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
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
double precision, | intent(in) | :: | gamma0 |
Polytropic exponent |
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