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)