diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:47 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:47 +0000 |
commit | cb186e758987c714b7d77709b45d6b90b6e4a280 (patch) | |
tree | 50b77aaff1f8e70979a58ec1ded18d739e2e34fb /src | |
parent | 5e3fd99e900f957b320e9dd428fa062fc1713980 (diff) |
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@physics.gatech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@462 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src')
-rw-r--r-- | src/GRHydro_Con2PrimM.F90 | 6 |
1 files changed, 1 insertions, 5 deletions
diff --git a/src/GRHydro_Con2PrimM.F90 b/src/GRHydro_Con2PrimM.F90 index c332aff..99260c0 100644 --- a/src/GRHydro_Con2PrimM.F90 +++ b/src/GRHydro_Con2PrimM.F90 @@ -656,11 +656,7 @@ subroutine Conservative2PrimitiveM(CCTK_ARGUMENTS) !!$ 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) - - if(tau(i,j,k).le.sdet*b2*0.5d0)then - tau(i,j,k) = GRHydro_tau_min + sdet*b2*0.5d0 - endif + tau(i,j,k) = sdet * (rho(i,j,k)*eps(i,j,k)+b2/2.0) end if |