User: cott Date: 2012/12/29 03:37 PM
Modified: /trunk/src/ EOS_Omni_Module.F90
Log: * update unit conversions based on new conversion factors
File Changes:
Directory: /trunk/src/ ======================
File [modified]: EOS_Omni_Module.F90 Delta lines: +14 -15 =================================================================== --- trunk/src/EOS_Omni_Module.F90 2012-12-29 12:57:38 UTC (rev 71) +++ trunk/src/EOS_Omni_Module.F90 2012-12-29 21:37:30 UTC (rev 72) @@ -2,25 +2,24 @@
implicit none
- ! Magic numbers in convenient units, calculated by someone whom we - ! trust to have worked very carefully - - real*8,parameter :: rho_gf = 1.61620075314614d-18 - real*8,parameter :: press_gf = 1.7982953469278d-39 - real*8,parameter :: eps_gf = 1.11265006d-21 - real*8,parameter :: time_gf = 2.03001708d+05 - real*8,parameter :: mass_gf = 5.02765209d-34 - real*8,parameter :: length_gf = 6.77140812d-06 + ! 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.18735016707159d17 - real*8,parameter :: inv_press_gf = 5.56082181777535d38 - real*8,parameter :: inv_eps_gf = 8.98755175549085d20 - real*8,parameter :: inv_time_gf = 4.92606692747629d-6 - real*8,parameter :: inv_mass_gf = 1.98899999860571d33 - real*8,parameter :: inv_length_gf = 147679.77092481d0 + 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