From ec464f0fb0bb2d9faed7b9de03fb6e6a48a38afe Mon Sep 17 00:00:00 2001 From: cott Date: Sat, 29 Dec 2012 21:37:30 +0000 Subject: * update unit conversions based on new conversion factors git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOS_Omni/trunk@72 8e189c6b-2ab8-4400-aa02-70a9cfce18b9 --- src/EOS_Omni_Module.F90 | 29 ++++++++++++++--------------- 1 file changed, 14 insertions(+), 15 deletions(-) diff --git a/src/EOS_Omni_Module.F90 b/src/EOS_Omni_Module.F90 index 1cc6719..2fbef68 100644 --- a/src/EOS_Omni_Module.F90 +++ b/src/EOS_Omni_Module.F90 @@ -2,25 +2,24 @@ module EOS_Omni_Module 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 -- cgit v1.2.3