/[ascend]/trunk/models/johnpye/fprops/helmholtz.c
ViewVC logotype

Diff of /trunk/models/johnpye/fprops/helmholtz.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1993 by jpye, Tue Feb 3 23:52:53 2009 UTC revision 1994 by jpye, Wed Feb 4 03:52:36 2009 UTC
# Line 246  double helmholtz_w(double T, double rho, Line 246  double helmholtz_w(double T, double rho,
246      double phir_dt = helm_resid_deltau(tau,delta,data);      double phir_dt = helm_resid_deltau(tau,delta,data);
247            
248      /* note similarities with helmholtz_cp */      /* note similarities with helmholtz_cp */
249      double temp1 = 1 + 2*delta*phir_d + SQ(delta)*phir_dd;      double temp1 = 1. + 2.*delta*phir_d + SQ(delta)*phir_dd;
250      double temp2 = 1 + delta*phir_d - delta*tau*phir_dt;      double temp2 = 1. + delta*phir_d - delta*tau*phir_dt;
251      double temp3 = SQ(tau)*(helm_ideal_tautau(tau,data->ideal) + helm_resid_tautau(tau,delta,data));      double temp3 = SQ(tau)*(helm_ideal_tautau(tau,data->ideal) + helm_resid_tautau(tau,delta,data));
252    
253      return sqrt(data->R * T *(temp1 - SQ(temp2)/temp3));      return sqrt(data->R * T * (temp1 - SQ(temp2)/temp3));
254    
255  }  }
256    

Legend:
Removed from v.1993  
changed lines
  Added in v.1994

john.pye@anu.edu.au
ViewVC Help
Powered by ViewVC 1.1.22