/[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 1869 by jpye, Wed Sep 17 06:37:45 2008 UTC revision 1870 by jpye, Wed Sep 17 06:45:10 2008 UTC
# Line 143  double helmholtz_s(double T, double rho, Line 143  double helmholtz_s(double T, double rho,
143      double tau = data->T_star / T;      double tau = data->T_star / T;
144      double delta = rho / data->rho_star;      double delta = rho / data->rho_star;
145    
146  #ifdef TEST  #ifdef ENTROPY_DEBUG
147      assert(data->rho_star!=0);      assert(data->rho_star!=0);
148      assert(T!=0);      assert(T!=0);
149      assert(!isnan(tau));      assert(!isnan(tau));
# Line 187  double helmholtz_a(double T, double rho, Line 187  double helmholtz_a(double T, double rho,
187      fprintf(stderr,"multiplying by RT = %f\n",data->R*T);      fprintf(stderr,"multiplying by RT = %f\n",data->R*T);
188  #endif  #endif
189    
190      return data->R * T * (helm_ideal(tau,delta,data->ideal) - helm_resid(tau,delta,data));      return data->R * T * (helm_ideal(tau,delta,data->ideal) + helm_resid(tau,delta,data));
191  }  }
192    
193    

Legend:
Removed from v.1869  
changed lines
  Added in v.1870

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