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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/02/11 10:30 AM

Modified:
 /trunk/src/
  GRHydro_Con2PrimAM.F90

Log:
 GRHydro: disable in-AH polytrope call by default for Avec
 
 From: Roland Haas <rhaas at tapir.caltech.edu>

File Changes:

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

File [modified]: GRHydro_Con2PrimAM.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/GRHydro_Con2PrimAM.F90	2013-01-20 18:20:48 UTC (rev 471)
+++ trunk/src/GRHydro_Con2PrimAM.F90	2013-02-11 16:30:02 UTC (rev 472)
@@ -373,7 +373,7 @@
                 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)