diff options
| author | reisswig <reisswig@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-02-03 23:09:46 +0000 |
|---|---|---|
| committer | reisswig <reisswig@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-02-03 23:09:46 +0000 |
| commit | f419a8fba9ea693f86fd59e05b4c58bdc3547814 (patch) | |
| tree | 68ba8d1299b33f72839ec1d327e6cf020a712486 /src/GRHydro_Eigenproblem.F90 | |
| parent | 8789e09dce73f0f74dc7fc578b0a8dd665079057 (diff) | |
Set cs2 to zero if cs2 < 0 and abs(cs2) < 1e-4
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@312 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Eigenproblem.F90')
| -rw-r--r-- | src/GRHydro_Eigenproblem.F90 | 9 |
1 files changed, 8 insertions, 1 deletions
diff --git a/src/GRHydro_Eigenproblem.F90 b/src/GRHydro_Eigenproblem.F90 index 35313fa..605f372 100644 --- a/src/GRHydro_Eigenproblem.F90 +++ b/src/GRHydro_Eigenproblem.F90 @@ -57,6 +57,7 @@ subroutine eigenvalues(handle,rho,velx,vely,velz,eps, & CCTK_REAL lam1,lam2,lam3,lamm,lamp,lamm_nobeta,lamp_nobeta CCTK_INT handle CCTK_REAL dpdrho,dpdeps,press + character*256 :: warnline ! begin EOS Omni vars integer :: n,keytemp,anyerr,keyerr(1) @@ -85,7 +86,13 @@ subroutine eigenvalues(handle,rho,velx,vely,velz,eps, & (1.0d0 + eps + press/rho) if(cs2.lt.0.0d0) then - call CCTK_WARN(0,"cs2 < 0! Check speed of sound calculation!") + if (abs(cs2) .gt. 1.0d-4) then + write(warnline,'(a50,6g16.7)') 'rho, dpdrho, press*dpdeps/rho**2, eps, press/rho: ', abs(cs2), rho, dpdrho, press * dpdeps / (rho**2), eps, press/rho + call CCTK_WARN(1,warnline) + call CCTK_WARN(0,"cs2 < 0! Check speed of sound calculation!") + else + cs2 = 0 + endif endif |