Fixed an array indexing error in tpx::Water::ldens

This commit is contained in:
Ray Speth 2012-02-27 18:15:27 +00:00
parent 236756fbf0
commit 1e5d4d50e1

View file

@ -191,8 +191,8 @@ double water::ldens()
if ((T < Tmn) || (T >= Tc)) {
set_Err(TempError); // Error("water::ldens",TempError,T);
}
for (i=0; i<=8; i++) {
sum+=D[i-1]*pow(1.0 - T/Tc, double(i)/3.0);
for (i=0; i<8; i++) {
sum+=D[i]*pow(1.0 - T/Tc, double(i+1)/3.0);
}
double density = Roc*(1+sum);
return density;