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@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