diff options
Diffstat (limited to 'src/nuc_eos/nuc_eos.F90')
| -rw-r--r-- | src/nuc_eos/nuc_eos.F90 | 29 |
1 files changed, 16 insertions, 13 deletions
diff --git a/src/nuc_eos/nuc_eos.F90 b/src/nuc_eos/nuc_eos.F90 index 18928cf..b576ce0 100644 --- a/src/nuc_eos/nuc_eos.F90 +++ b/src/nuc_eos/nuc_eos.F90 @@ -17,6 +17,9 @@ ! 1 -> coming in with temperature ! 2 -> coming in with entropy ! + +#include "cctk.h" + subroutine nuc_eos_full(xrho,xtemp,xye,xenr,xprs,xent,xcs2,xdedt,& xdpderho,xdpdrhoe,xxa,xxh,xxn,xxp,xabar,xzbar,xmu_e,xmu_n,xmu_p, & xmuhat,keytemp,keyerr,rfeps) @@ -45,11 +48,11 @@ subroutine nuc_eos_full(xrho,xtemp,xye,xenr,xprs,xent,xcs2,xdedt,& keyerrr = 0 if(xrho.gt.eos_rhomax) then - stop "nuc_eos: rho > rhomax" + call CCTK_WARN(0, "nuc_eos: rho > rhomax") endif if(xrho.lt.eos_rhomin) then - stop "nuc_eos: rho < rhomin" + call CCTK_WARN (0, "nuc_eos: rho < rhomin") endif if(xye.gt.eos_yemax) then @@ -64,11 +67,11 @@ subroutine nuc_eos_full(xrho,xtemp,xye,xenr,xprs,xent,xcs2,xdedt,& if(keytemp.eq.1) then if(xtemp.gt.eos_tempmax) then - stop "nuc_eos: temp > tempmax" + call CCTK_WARN (0, "nuc_eos: temp > tempmax") endif if(xtemp.lt.eos_tempmin) then - stop "nuc_eos: temp < tempmin" + call CCTK_WARN (0, "nuc_eos: temp < tempmin") endif endif @@ -84,7 +87,7 @@ subroutine nuc_eos_full(xrho,xtemp,xye,xenr,xprs,xent,xcs2,xdedt,& !need to find temperature based on xeps call findtemp(lr,lt,y,leps,keyerrt,rfeps) if(keyerrt.ne.0) then - stop "Did not find temperature" + call CCTK_WARN (0, "Did not find temperature") endif xtemp = 10.0d0**lt @@ -206,7 +209,7 @@ subroutine nuc_eos_short(xrho,xtemp,xye,xenr,xprs,xent,xcs2,xdedt,& keyerrr = 0 if(xrho.gt.eos_rhomax) then - stop "nuc_eos: rho > rhomax" + call CCTK_WARN (0, "nuc_eos: rho > rhomax") endif if(xrho.lt.eos_rhomin*1.2d0) then @@ -227,7 +230,7 @@ subroutine nuc_eos_short(xrho,xtemp,xye,xenr,xprs,xent,xcs2,xdedt,& if(keytemp.eq.1) then if(xtemp.gt.eos_tempmax) then - stop "nuc_eos: temp > tempmax" + call CCTK_WARN (0, "nuc_eos: temp > tempmax") endif if(xtemp.lt.eos_tempmin) then @@ -329,7 +332,7 @@ subroutine nuc_eos_press_eps(xrho,xtemp,xye,xenr,xprs,& keyerrt = 0 if(xrho.gt.eos_rhomax) then - stop "nuc_eos: rho > rhomax" + call CCTK_WARN (0, "nuc_eos: rho > rhomax") endif if(xrho.lt.eos_rhomin*1.2d0) then @@ -349,7 +352,7 @@ subroutine nuc_eos_press_eps(xrho,xtemp,xye,xenr,xprs,& if(keytemp.eq.1) then if(xtemp.gt.eos_tempmax) then - stop "nuc_eos: temp > tempmax" + call CCTK_WARN (0, "nuc_eos: temp > tempmax") endif if(xtemp.lt.eos_tempmin) then @@ -361,7 +364,7 @@ subroutine nuc_eos_press_eps(xrho,xtemp,xye,xenr,xprs,& keyerr = 0 if(keytemp.gt.1) then - stop "eos_nuc_press does not support keytemp other than 0 and 1" + call CCTK_WARN (0, "eos_nuc_press does not support keytemp other than 0 and 1") endif lr = log10(xrho) @@ -408,7 +411,7 @@ subroutine nuc_eos_dpdr_dpde(xrho,xtemp,xye,xenr,xdpdrhoe,& keyerrt = 0 if(xrho.gt.eos_rhomax) then - stop "nuc_eos: rho > rhomax" + call CCTK_WARN (0, "nuc_eos: rho > rhomax") endif if(xrho.lt.eos_rhomin*1.2d0) then @@ -428,7 +431,7 @@ subroutine nuc_eos_dpdr_dpde(xrho,xtemp,xye,xenr,xdpdrhoe,& if(keytemp.eq.1) then if(xtemp.gt.eos_tempmax) then - stop "nuc_eos: temp > tempmax" + call CCTK_WARN (0, "nuc_eos: temp > tempmax") endif if(xtemp.lt.eos_tempmin) then @@ -455,7 +458,7 @@ subroutine nuc_eos_dpdr_dpde(xrho,xtemp,xye,xenr,xdpdrhoe,& xtemp = 10.0d0**lt elseif(keytemp.gt.1) then - stop "eos_nuc_press does not support keytemp > 1" + call CCTK_WARN (0, "eos_nuc_press does not support keytemp > 1") endif ! have temperature, proceed: |