diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-05-29 16:34:58 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-05-29 16:34:58 +0000 |
| commit | dabcbc8fd072964f70c41d2039f9e010457f8edc (patch) | |
| tree | c368a7ab0b85aa90cea03337a83b5c8775fdaf3e | |
| parent | ad183a5ea2629e9a40b8d82f57b18531889a88a6 (diff) | |
GRHydro: add error checking after call for cs2 in Eigenproblem
From: Christian David Ott <cott@zwicky-b.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@527 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
| -rw-r--r-- | src/GRHydro_Eigenproblem.F90 | 11 |
1 files changed, 11 insertions, 0 deletions
diff --git a/src/GRHydro_Eigenproblem.F90 b/src/GRHydro_Eigenproblem.F90 index 90c46f8..0b7b3c5 100644 --- a/src/GRHydro_Eigenproblem.F90 +++ b/src/GRHydro_Eigenproblem.F90 @@ -146,6 +146,8 @@ subroutine eigenvalues_hot(handle,keytemp,ii,jj,kk,rho,velx,vely,velz,eps, & CCTK_REAL lam1,lam2,lam3,lamm,lamp,lamm_nobeta,lamp_nobeta CCTK_INT handle,ii,jj,kk + character(len=512) :: warnline + ! begin EOS Omni vars integer :: n,keytemp,anyerr,keyerr(1) real*8 :: xpress,xeps @@ -159,6 +161,15 @@ subroutine eigenvalues_hot(handle,keytemp,ii,jj,kk,rho,velx,vely,velz,eps, & !!$ Set required fluid quantities call EOS_Omni_cs2(handle,keytemp,GRHydro_eos_rf_prec,n,& rho,eps,temp,ye,cs2,keyerr,anyerr) + if(anyerr.ne.0) then + !$OMP CRITICAL + call CCTK_WARN(1,"EOS ERROR in eigenvalues_hot") + write(warnline,"(A10,i5,A10,i5)") "keyerr: ", keyerr, "keytemp: ",keytemp + call CCTK_WARN(1,warnline) + write(warnline,"(1P10E15.6)") rho,eps,temp,ye,cs2 + call CCTK_WARN(0,warnline) + !$OMP END CRITICAL + endif vlowx = gxx*velx + gxy*vely + gxz*velz vlowy = gxy*velx + gyy*vely + gyz*velz |