[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