diff options
| -rw-r--r-- | src/GRHydro_Eigenproblem.F90 | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/src/GRHydro_Eigenproblem.F90 b/src/GRHydro_Eigenproblem.F90 index e65a0e6..35313fa 100644 --- a/src/GRHydro_Eigenproblem.F90 +++ b/src/GRHydro_Eigenproblem.F90 @@ -84,6 +84,11 @@ subroutine eigenvalues(handle,rho,velx,vely,velz,eps, & cs2 = (dpdrho + press * dpdeps / (rho**2))/ & (1.0d0 + eps + press/rho) + if(cs2.lt.0.0d0) then + call CCTK_WARN(0,"cs2 < 0! Check speed of sound calculation!") + endif + + vlowx = gxx*velx + gxy*vely + gxz*velz vlowy = gxy*velx + gyy*vely + gyz*velz vlowz = gxz*velx + gyz*vely + gzz*velz |