[Commits] [svn:einsteintoolkit] incoming/EOS_Omni/src/ (Rev. 20)

[cott at tapir.caltech.edu]

User: cott
Date: 2010/11/16 09:19 AM

Modified:
 /EOS_Omni/src/
  EOS_Omni_SingleVarCalls.F90
 /EOS_Omni/src/nuc_eos/
  nuc_eos.F90, readtable.F90

Log:
 * minor mods related to hot EOS

File Changes:

Directory: /EOS_Omni/src/nuc_eos/
=================================

File [modified]: nuc_eos.F90
Delta lines: +6 -2
===================================================================
--- EOS_Omni/src/nuc_eos/nuc_eos.F90	2010-11-14 06:17:54 UTC (rev 19)
+++ EOS_Omni/src/nuc_eos/nuc_eos.F90	2010-11-16 15:19:04 UTC (rev 20)
@@ -195,7 +195,9 @@
   endif
 
   if(xye.lt.eos_yemin) then
-     stop "nuc_eos: ye < yemin"
+     keyerr = 102
+     write(6,*) "ye: ", xye
+     stop "nuc_eos ye < yemin"
   endif
 
   if(keytemp.eq.1) then
@@ -298,7 +300,9 @@
   endif
 
   if(xye.lt.eos_yemin) then
-     stop "nuc_eos: ye < yemin"
+     keyerr = 102
+     write(6,*) "ye: ", xye
+     stop "nuc_eos ye < yemin"
   endif
 
   if(keytemp.eq.1) then

File [modified]: readtable.F90
Delta lines: +0 -1
===================================================================
--- EOS_Omni/src/nuc_eos/readtable.F90	2010-11-14 06:17:54 UTC (rev 19)
+++ EOS_Omni/src/nuc_eos/readtable.F90	2010-11-16 15:19:04 UTC (rev 20)
@@ -210,7 +210,6 @@
     stop "Problem reading EOS table file"
   endif
 
-
   call h5fclose_f (file_id,error)
 
 ! not a good idea to call this function as it will screw recovery

Directory: /EOS_Omni/src/
=========================

File [modified]: EOS_Omni_SingleVarCalls.F90
Delta lines: +7 -0
===================================================================
--- EOS_Omni/src/EOS_Omni_SingleVarCalls.F90	2010-11-14 06:17:54 UTC (rev 19)
+++ EOS_Omni/src/EOS_Omni_SingleVarCalls.F90	2010-11-16 15:19:04 UTC (rev 20)
@@ -96,6 +96,13 @@
            xtemp = temp(i)
            xye = ye(i)
            xenr = eps(i) * inv_eps_gf
+           if(xye.lt.0.1) then
+              write (warnstring,"(A6,1P10E16.7)") "blah: ",xye,xrho
+              call CCTK_WARN(1,warnstring)
+              keyerr(i) = 102
+              anyerr = 1
+              return
+           endif
            call nuc_eos_press_eps(xrho,xtemp,xye,xenr,xprs,& 
                 keytemp,keyerr(i),rf_precision)