diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-02-11 16:30:02 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-02-11 16:30:02 +0000 |
| commit | b770b656fbd9926096a99da9b72dc39c44a4fcf9 (patch) | |
| tree | 9b98b4a36ad9f8b6b8f28a7d362a089be74c4118 | |
| parent | 426c5ef849c66987a4df7143e097ebc18efe4120 (diff) | |
GRHydro: disable in-AH polytrope call by default for Avec
From: Roland Haas <rhaas@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@472 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
| -rw-r--r-- | src/GRHydro_Con2PrimAM.F90 | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/src/GRHydro_Con2PrimAM.F90 b/src/GRHydro_Con2PrimAM.F90 index 788e3b8..5bdf02a 100644 --- a/src/GRHydro_Con2PrimAM.F90 +++ b/src/GRHydro_Con2PrimAM.F90 @@ -373,7 +373,7 @@ subroutine Conservative2PrimitiveAM(CCTK_ARGUMENTS) local_pgam=log(xpress(1)/local_K)/log(xrho(1)) sc = local_K*dens(i,j,k) - if(sdet.ge.sqrtdet_thr) then + if(sqrtdet_thr.gt.0d0 .and. sdet.ge.sqrtdet_thr) then GRHydro_C2P_failed(i,j,k) = 0 rho_tmp = rho(i,j,k) |