[Commits] [svn:einsteintoolkit] GRHydro/trunk/ (Rev. 363)
rhaas at tapir.caltech.edu
rhaas at tapir.caltech.edu
Wed Jun 20 22:51:18 CDT 2012
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.
More information about the Commits
mailing list