1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
|
module EOS_Omni_Module
implicit none
! conversion factors between cgs and M_Sun = c = G = 1
! see EOS_Omni/doc/units.py
real*8,parameter :: rho_gf = 1.61887093132742d-18
real*8,parameter :: press_gf = 1.80123683248503d-39
real*8,parameter :: eps_gf = 1.11265005605362d-21
real*8,parameter :: time_gf = 2.03040204956746d05
real*8,parameter :: mass_gf = 5.02916918125126d-34
real*8,parameter :: length_gf = 6.77269222552442d-06
! Inverses of the numbers above, calculated manually instead of by
! the compiler
real*8,parameter :: inv_rho_gf = 6.17714470405638d17
real*8,parameter :: inv_press_gf = 5.55174079257738d38
real*8,parameter :: inv_eps_gf = 8.98755178736818d20
real*8,parameter :: inv_time_gf = 4.92513293223396d-6
real*8,parameter :: inv_mass_gf = 1.98840000000000d33
real*8,parameter :: inv_length_gf = 1.47651770773117d05
real*8,parameter :: clite = 2.99792458d10
real*8,parameter :: cliteinv2 = 1.11265005605362d-21
! These values are set by EOS_Omni_Startup
real*8 :: hybrid_k2 = 0.0d0
! stuff for the cold, tabulated EOS with a gamma law
! set by the reader routine
integer :: coldeos_nrho = 0
real*8 :: coldeos_gammath = 0.0d0
real*8 :: coldeos_rhomin = 0.0d0
real*8 :: coldeos_rhomax = 0.0d0
real*8 :: coldeos_kappa = 0.0d0
real*8 :: coldeos_thfac = 1.0d0
real*8 :: coldeos_dlrho = 1.0d0
real*8 :: coldeos_dlrhoi = 1.0d0
real*8, allocatable :: coldeos_logrho(:)
real*8, allocatable :: coldeos_eps(:)
real*8, allocatable :: coldeos_gamma(:)
real*8, allocatable :: coldeos_cs2(:)
end module EOS_Omni_Module
|