diff options
| author | cott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9> | 2010-11-16 15:19:04 +0000 |
|---|---|---|
| committer | cott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9> | 2010-11-16 15:19:04 +0000 |
| commit | 44ff1e4799a014f4bac2b2fcb2cc34d75a2c5edc (patch) | |
| tree | 0f66d37119395fe08bb81eb79e8a9dde9219a1ee /src | |
| parent | 68fef2fe3544197073744ccb7e69be1a34763225 (diff) | |
* minor mods related to hot EOS
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOS_Omni/EOS_Omni@20 8e189c6b-2ab8-4400-aa02-70a9cfce18b9
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 |