User: rhaas Date: 2013/05/29 11:34 AM
Modified: /trunk/src/ GRHydro_Eigenproblem.F90
Log: GRHydro: add error checking after call for cs2 in Eigenproblem
From: Christian David Ott <cott@zwicky-b.(none)>
File Changes:
Directory: /trunk/src/ ======================
File [modified]: GRHydro_Eigenproblem.F90 Delta lines: +11 -0 =================================================================== --- trunk/src/GRHydro_Eigenproblem.F90 2013-05-29 16:34:50 UTC (rev 526) +++ trunk/src/GRHydro_Eigenproblem.F90 2013-05-29 16:34:58 UTC (rev 527) @@ -146,6 +146,8 @@ 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 @@ !!$ 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