diff options
Diffstat (limited to 'src')
| -rw-r--r-- | src/EOS_Omni_SingleVarCalls.F90 | 7 | ||||
| -rw-r--r-- | src/nuc_eos/nuc_eos.F90 | 8 | ||||
| -rw-r--r-- | src/nuc_eos/readtable.F90 | 1 |
3 files changed, 13 insertions, 3 deletions
diff --git a/src/EOS_Omni_SingleVarCalls.F90 b/src/EOS_Omni_SingleVarCalls.F90 index 72abe97..3cd6c38 100644 --- a/src/EOS_Omni_SingleVarCalls.F90 +++ b/src/EOS_Omni_SingleVarCalls.F90 @@ -96,6 +96,13 @@ subroutine EOS_Omni_EOS_Press(eoskey,keytemp,rf_precision,npoints,& 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) diff --git a/src/nuc_eos/nuc_eos.F90 b/src/nuc_eos/nuc_eos.F90 index bc0f912..e84ac82 100644 --- a/src/nuc_eos/nuc_eos.F90 +++ b/src/nuc_eos/nuc_eos.F90 @@ -195,7 +195,9 @@ subroutine nuc_eos_short(xrho,xtemp,xye,xenr,xprs,xent,xcs2,xdedt,& 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 @@ subroutine nuc_eos_press_eps(xrho,xtemp,xye,xenr,xprs,& 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 diff --git a/src/nuc_eos/readtable.F90 b/src/nuc_eos/readtable.F90 index f332fab..1eee600 100644 --- a/src/nuc_eos/readtable.F90 +++ b/src/nuc_eos/readtable.F90 @@ -210,7 +210,6 @@ subroutine nuc_eos_readtable(eos_filename) 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 |