[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 527)

[rhaas at tapir.caltech.edu]

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 at 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