User: rhaas Date: 2012/06/20 10:51 PM
Modified: /trunk/ param.ccl /trunk/src/ GRHydro_Con2Prim.F90
Log: * introduce parameters GRHydro_Y_e_min and GRHydro_Y_e_max to fix issues with EOS tables that don't extend sufficiently far. Enforcement happens in Con2Prim.
File Changes:
Directory: /trunk/src/ ======================
File [modified]: GRHydro_Con2Prim.F90 Delta lines: +2 -1 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2012-06-21 03:51:16 UTC (rev 362) +++ trunk/src/GRHydro_Con2Prim.F90 2012-06-21 03:51:18 UTC (rev 363) @@ -150,7 +150,8 @@ endif
if(evolve_Y_e.ne.0) then - Y_e(i,j,k) = Y_e_con(i,j,k) / dens(i,j,k) + Y_e(i,j,k) = max(min(Y_e_con(i,j,k) / dens(i,j,k),GRHydro_Y_e_max),& + GRHydro_Y_e_min) endif
if ( dens(i,j,k) .le. sqrt(det)*GRHydro_rho_min*(1.d0+GRHydro_atmo_tolerance) ) then
Directory: /trunk/ ==================
File [modified]: param.ccl Delta lines: +13 -2 =================================================================== --- trunk/param.ccl 2012-06-21 03:51:16 UTC (rev 362) +++ trunk/param.ccl 2012-06-21 03:51:18 UTC (rev 363) @@ -288,14 +288,25 @@
REAL GRHydro_hot_atmo_temp "Temperature of the hot atmosphere in MeV" STEERABLE=ALWAYS { - 0:* :: "Larger than 0 MeV" + 0.0:* :: "Larger than 0 MeV" } 0.5e0
REAL GRHydro_hot_atmo_Y_e "Y_e of the hot atmosphere" STEERABLE=ALWAYS { - 0:* :: "Larger than 0" + 0.0:* :: "Larger than 0" } 0.5e0
+REAL GRHydro_Y_e_min "minimum allowed Y_e" STEERABLE=ALWAYS +{ + 0.0:* :: "Larger than or equal to zero" +} 0.0 + +REAL GRHydro_Y_e_max "maximum allowed Y_e" STEERABLE=ALWAYS +{ + 1.0:* :: "Larger than or equal to zero; 1 is default" +} 1.0 + + ####################### Other Parameters ##############################
#Parameters controlling conservative <-> primitive change.