diff options
author | cott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-12-07 22:33:23 +0000 |
---|---|---|
committer | cott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-12-07 22:33:23 +0000 |
commit | 6c6dde52a11ba3b7841fdbbff4f7f3cfc120c077 (patch) | |
tree | 9de8a40f9e953c96d1d52d97cb5dda66f853c7ef /src/GRHydro_Eigenproblem.F90 | |
parent | a6ff91cfc7a50f74d5de38777f61a5bf7d037f3d (diff) |
* introduce check that cs2 > 0.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@309 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Eigenproblem.F90')
-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 |