aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorcott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9>2011-06-11 00:32:35 +0000
committercott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9>2011-06-11 00:32:35 +0000
commit42da92056c93da30be8fd4b8286f71aa92f345be (patch)
treec9b62425c2031e9fcdf4385ae1e8deca14390375
parent16fb957603005e7dcc10a1864dcc2db7090ed8ba (diff)
* fix dpdrho issue
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOS_Omni/trunk@48 8e189c6b-2ab8-4400-aa02-70a9cfce18b9
-rw-r--r--src/EOS_Omni_SingleVarCalls.F904
1 files changed, 3 insertions, 1 deletions
diff --git a/src/EOS_Omni_SingleVarCalls.F90 b/src/EOS_Omni_SingleVarCalls.F90
index 7f1f825..783f2c1 100644
--- a/src/EOS_Omni_SingleVarCalls.F90
+++ b/src/EOS_Omni_SingleVarCalls.F90
@@ -203,7 +203,9 @@ subroutine EOS_Omni_EOS_DPressByDRho(eoskey,keytemp,rf_precision,npoints,&
press_gf * hybrid_k1_cgs * inv_rho_gf**hybrid_gamma1 * &
hybrid_rho_nuc**(hybrid_gamma1 - 1.d0)
- dpdrhoe(i) = hybrid_dp_poly + hybrid_dp_th1 + hybrid_dp_th2
+ dpdrhoe(i) = hybrid_dp_poly + max(0.0d0,hybrid_dp_th1 + hybrid_dp_th2)
+
+
enddo