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