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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/01/14 08:23 AM

Modified:
 /trunk/
  param.ccl
 /trunk/src/
  GRHydro_Con2PrimM.F90

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

File Changes:

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

File [modified]: GRHydro_Con2PrimM.F90
Delta lines: +2 -2
===================================================================
--- trunk/src/GRHydro_Con2PrimM.F90	2013-01-14 14:23:23 UTC (rev 451)
+++ trunk/src/GRHydro_Con2PrimM.F90	2013-01-14 14:23:24 UTC (rev 452)
@@ -270,7 +270,7 @@
            if(evolve_temper.eq.0) then
             
 
-              if(sdet.ge.sqrtdet_thr) then
+              if(sqrtdet_thr.gt.0d0 .and. sdet.ge.sqrtdet_thr) then
                     d2 = dens(i,j,k)**2
                     s2 = uxx*scon(i,j,k,1)**2 + uyy*scon(i,j,k,2)**2 &
                                               + uzz*scon(i,j,k,3)**2 &
@@ -359,7 +359,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)

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

File [modified]: param.ccl
Delta lines: +2 -1
===================================================================
--- trunk/param.ccl	2013-01-14 14:23:23 UTC (rev 451)
+++ trunk/param.ccl	2013-01-14 14:23:24 UTC (rev 452)
@@ -354,7 +354,8 @@
 REAL sqrtdet_thr "Threshold to apply cons rescalings deep inside the horizon" STEERABLE=ALWAYS
 {
  1.0:    :: "Larger values guarantees this sort of rescaling only deep inside the horizon"
-} 30.0
+ -1.0    :: "Do not apply limit"
+} -1.0
 
 REAL max_magnetic_to_gas_pressure_ratio "consider pressure to be magnetically dominated if magnetic pressure to gas pressure ratio is higher than this" STEERABLE=ALWAYS
 {