[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 282)

[roland.haas at physics.gatech.edu]

User: rhaas
Date: 2011/10/14 03:35 AM

Modified:
 /trunk/src/
  GRHydro_EoSChangeGamma.F90, GRHydro_Prim2Con.F90

Log:
 correct formula used to compute eps in EoSChangeK
 
 this only shows up if we use a general eos (not a polytype) during evolution

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_EoSChangeGamma.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/GRHydro_EoSChangeGamma.F90	2011-10-03 16:05:30 UTC (rev 281)
+++ trunk/src/GRHydro_EoSChangeGamma.F90	2011-10-14 08:35:29 UTC (rev 282)
@@ -195,7 +195,7 @@
   end if
   
   press = local_k * rho**local_gamma
-  eps = (local_gamma - 1.d0) * local_k * rho**local_gamma
+  eps = local_k / (local_gamma - 1.d0) * rho**(local_gamma-1)
 
   call Primitive2ConservativePolyCells(CCTK_ARGUMENTS)
 

File [modified]: GRHydro_Prim2Con.F90
Delta lines: +3 -3
===================================================================
--- trunk/src/GRHydro_Prim2Con.F90	2011-10-03 16:05:30 UTC (rev 281)
+++ trunk/src/GRHydro_Prim2Con.F90	2011-10-14 08:35:29 UTC (rev 282)
@@ -541,14 +541,14 @@
         avg_detl = SPATIAL_DETERMINANT(g11l,g12l,g13l,g22l, g23l,g33l)
         avg_detr = SPATIAL_DETERMINANT(g11r,g12r,g13r,g22r, g23r,g33r)
 
-        call prim2conpolytype(GRHydro_eos_handle, g11l,g12l,g13l,&
+        call prim2conpolytype(GRHydro_polytrope_handle, g11l,g12l,g13l,&
              g22l,g23l,g33l, &
              avg_detl,densminus(i,j,k),sxminus(i,j,k),&
              syminus(i,j,k),szminus(i,j,k),tauminus(i,j,k),rhominus(i,j,k), &
              velxminus(i,j,k),velyminus(i,j,k),velzminus(i,j,k),&
              epsminus(i,j,k),pressminus(i,j,k),w_lorentzminus(i, j, k))
  
-        call prim2conpolytype(GRHydro_eos_handle, g11r,g12r,g13r,&
+        call prim2conpolytype(GRHydro_polytrope_handle, g11r,g12r,g13r,&
              g22r,g23r,g33r, &
              avg_detr, densplus(i,j,k),sxplus(i,j,k),&
              syplus(i,j,k),szplus(i,j ,k),tauplus(i,j,k),&
@@ -705,7 +705,7 @@
         det = SPATIAL_DETERMINANT(g11(i,j,k),g12(i,j,k),g13(i,j,k),\
              g22(i,j,k),g23(i,j,k),g33(i,j,k))
 
-        call prim2conpolytype(GRHydro_eos_handle,g11(i,j,k),g12(i,j,k),&
+        call prim2conpolytype(GRHydro_polytrope_handle,g11(i,j,k),g12(i,j,k),&
              g13(i,j,k),g22(i,j,k),g23(i,j,k),g33(i,j,k),&
              det, dens(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k),&
              tau(i,j,k),rho(i,j,k),velx(i,j,k),vely(i,j,k),velz(i,j,k),&