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