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