k_lorene2cu_pwp Function

public pure function k_lorene2cu_pwp(gamma0)

Uses

    • constants
  • proc~~k_lorene2cu_pwp~~UsesGraph proc~k_lorene2cu_pwp k_lorene2cu_pwp constants constants proc~k_lorene2cu_pwp->constants

Compute the constant to convert the polytropic constant from the units used in \(\texttt{LORENE}\) for the piecewise polytropic \(\mathrm{EOS}\), to the units used in \(\texttt{SPHINCS}\)

FT xx.xx.2020


Arguments

Type IntentOptional Attributes Name
double precision, intent(in) :: gamma0

Polytropic exponent

Return Value double precision


Called by

proc~~k_lorene2cu_pwp~~CalledByGraph proc~k_lorene2cu_pwp k_lorene2cu_pwp proc~print_bns_properties print_bns_properties proc~print_bns_properties->proc~k_lorene2cu_pwp proc~print_diffstar_properties print_diffstar_properties proc~print_diffstar_properties->proc~k_lorene2cu_pwp proc~read_bns_properties read_bns_properties proc~read_bns_properties->proc~k_lorene2cu_pwp interface~print_bns_properties print_bns_properties proc~read_bns_properties->interface~print_bns_properties proc~read_diffstar_properties read_diffstar_properties proc~read_diffstar_properties->proc~k_lorene2cu_pwp interface~print_diffstar_properties print_diffstar_properties proc~read_diffstar_properties->interface~print_diffstar_properties program~write_par_eos write_par_eos program~write_par_eos->proc~k_lorene2cu_pwp interface~print_bns_properties->proc~print_bns_properties interface~print_diffstar_properties->proc~print_diffstar_properties interface~read_bns_properties read_bns_properties interface~read_bns_properties->proc~read_bns_properties interface~read_diffstar_properties read_diffstar_properties interface~read_diffstar_properties->proc~read_diffstar_properties proc~construct_bnslorene construct_bnslorene proc~construct_bnslorene->interface~read_bns_properties proc~construct_diffstarlorene construct_diffstarlorene proc~construct_diffstarlorene->interface~read_diffstar_properties interface~construct_bnslorene construct_bnslorene interface~construct_bnslorene->proc~construct_bnslorene interface~construct_diffstarlorene construct_diffstarlorene interface~construct_diffstarlorene->proc~construct_diffstarlorene

Contents

Source Code


Source Code

  PURE FUNCTION k_lorene2cu_pwp( gamma0 )

    !****************************************************************
    !
    !# Compute the constant to convert the polytropic constant \(K\)
    !  from the units used in |lorene| for the piecewise polytropic
    !  |eos|, to the units used in |sphincs|
    !
    !  FT xx.xx.2020
    !
    !****************************************************************

    USE constants,  ONLY: m2cm

    IMPLICIT NONE

    DOUBLE PRECISION, INTENT(IN) :: gamma0
    !! Polytropic exponent \(\gamma_0\)
    DOUBLE PRECISION :: k_lorene2cu_pwp

    ! LORENE has K0 in units of (g cm^{-3})^{1-gamma0} for the piecewise
    ! polytropes. This factor writes it in SPHINCS units

    k_lorene2cu_pwp= ( MSun/((MSun_geo*km2m*m2cm)**3) )**(gamma0 - one)


  END FUNCTION k_lorene2cu_pwp