From b770b656fbd9926096a99da9b72dc39c44a4fcf9 Mon Sep 17 00:00:00 2001 From: rhaas Date: Mon, 11 Feb 2013 16:30:02 +0000 Subject: GRHydro: disable in-AH polytrope call by default for Avec From: Roland Haas git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@472 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45 --- src/GRHydro_Con2PrimAM.F90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) 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) -- cgit v1.2.3