User: cott
Date: 2013/03/23 03:27 PM
Modified:
/trunk/src/
EOS_Omni_MultiVarCalls.F90, EOS_Omni_SingleVarCalls.F90
Log:
* fix an issue with coldEOS+gamma-law (thermal pressure was not added)
* improve a multi-var EOS call for dpdrho and dpdeps
File Changes:
Directory: /trunk/src/
======================
File [modified]: EOS_Omni_MultiVarCalls.F90
Delta lines: +2 -2
===================================================================
--- trunk/src/EOS_Omni_MultiVarCalls.F90 2013-03-23 20:10:20 UTC (rev 78)
+++ trunk/src/EOS_Omni_MultiVarCalls.F90 2013-03-23 20:27:23 UTC (rev 79)
@@ -179,8 +179,8 @@
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
File [modified]: EOS_Omni_SingleVarCalls.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/EOS_Omni_SingleVarCalls.F90 2013-03-23 20:10:20 UTC (rev 78)
+++ trunk/src/EOS_Omni_SingleVarCalls.F90 2013-03-23 20:27:23 UTC (rev 79)
@@ -151,7 +151,7 @@
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