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

[cott at tapir.caltech.edu]

User: cott
Date: 2011/12/07 04:33 PM

Modified:
 /trunk/src/
  GRHydro_Eigenproblem.F90

Log:
 * introduce check that cs2 > 0.

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_Eigenproblem.F90
Delta lines: +5 -0
===================================================================
--- trunk/src/GRHydro_Eigenproblem.F90	2011-12-06 05:51:39 UTC (rev 308)
+++ trunk/src/GRHydro_Eigenproblem.F90	2011-12-07 22:33:23 UTC (rev 309)
@@ -84,6 +84,11 @@
   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