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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/01/14 08:23 AM

Modified:
 /trunk/src/
  GRHydro_Con2PrimM.F90

Log:
 GRHydro: correct tau when resetting to atmosphere
 
 before a term rho*b2 was present when the correct one is only b2
 
 From: Roland Haas <roland.haas at physics.gatech.edu>

File Changes:

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

File [modified]: GRHydro_Con2PrimM.F90
Delta lines: +1 -5
===================================================================
--- trunk/src/GRHydro_Con2PrimM.F90	2013-01-14 14:23:45 UTC (rev 461)
+++ trunk/src/GRHydro_Con2PrimM.F90	2013-01-14 14:23:47 UTC (rev 462)
@@ -656,12 +656,8 @@
 
               !!$ tau does need to take into account the existing B-field
               !!$ with w_lorentz=1, we find tau = sqrtdet*(rho (1+eps+b^2/2)) - dens  [Press drops out]
-              tau(i,j,k)  = sdet * (rho(i,j,k)*(1.0+eps(i,j,k)+b2/2.0)) - dens(i,j,k)
+              tau(i,j,k)  = sdet * (rho(i,j,k)*eps(i,j,k)+b2/2.0)
 
-              if(tau(i,j,k).le.sdet*b2*0.5d0)then
-                tau(i,j,k) = GRHydro_tau_min + sdet*b2*0.5d0
-              endif
-
            end if
 
 !           ! Again, reset 3-metric only for con2prim inversion. Restoring