[Commits] [svn:einsteintoolkit] EOS_Omni/trunk/src/ (Rev. 72)

[cott at tapir.caltech.edu]

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