aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorcott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9>2013-03-23 20:27:23 +0000
committercott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9>2013-03-23 20:27:23 +0000
commitc16b69c8789baabd6f2fc3a9263c5e009b99a908 (patch)
treec399ca8dd2418d0f128dbcfddebf40e4397654f0
parent4530e894a8002d16f7eb7c822a79c4038c6b7682 (diff)
* fix an issue with coldEOS+gamma-law (thermal pressure was not added)
* improve a multi-var EOS call for dpdrho and dpdeps git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOS_Omni/trunk@79 8e189c6b-2ab8-4400-aa02-70a9cfce18b9
-rw-r--r--src/EOS_Omni_MultiVarCalls.F904
-rw-r--r--src/EOS_Omni_SingleVarCalls.F902
2 files changed, 3 insertions, 3 deletions
diff --git a/src/EOS_Omni_MultiVarCalls.F90 b/src/EOS_Omni_MultiVarCalls.F90
index 417ae99..e6a37b3 100644
--- a/src/EOS_Omni_MultiVarCalls.F90
+++ b/src/EOS_Omni_MultiVarCalls.F90
@@ -179,8 +179,8 @@ subroutine EOS_Omni_EOS_dpderho_dpdrhoe(eoskey,keytemp,rf_precision,npoints,&
xtemp = temp(i)
xye = ye(i)
xenr = eps(i) * inv_eps_gf
- call nuc_eos_dpdr_dpde(xrho,xtemp,xye,xenr,xdpdrhoe,&
- xdpderho,&
+ call nuc_eos_dpdr_dpde(xrho,xtemp,xye,xenr, &
+ xdpderho, xdpdrhoe,&
keytemp,keyerr(i),rf_precision)
if(keyerr(i).ne.0) then
diff --git a/src/EOS_Omni_SingleVarCalls.F90 b/src/EOS_Omni_SingleVarCalls.F90
index f789806..fc7a90b 100644
--- a/src/EOS_Omni_SingleVarCalls.F90
+++ b/src/EOS_Omni_SingleVarCalls.F90
@@ -151,7 +151,7 @@ subroutine EOS_Omni_EOS_Press(eoskey,keytemp,rf_precision,npoints,&
press_cold = coldeos_kappa * rho(i)**gamma
press_th = coldeos_thfac*(coldeos_gammath - 1.0d0)*rho(i)*eps_th
- press(i) = press_cold !+ press_th
+ press(i) = press_cold + press_th
enddo
case DEFAULT