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

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

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

revision 2300 by jpye, Mon Aug 16 08:46:28 2010 UTC revision 2301 by jpye, Sat Aug 21 13:28:35 2010 UTC
# Line 41  Line 41 
41    
42  #define SQ(X) ((X)*(X))  #define SQ(X) ((X)*(X))
43    
44  #define THROW_FPE  //#define THROW_FPE
45    
46  #ifdef THROW_FPE  #ifdef THROW_FPE
47  #define _GNU_SOURCE  #define _GNU_SOURCE
# Line 403  int fprops_sat_hf(double hf, double *Tsa Line 403  int fprops_sat_hf(double hf, double *Tsa
403      while(i++ < 60){      while(i++ < 60){
404          assert(T1 >= d->T_t - 1e-4);          assert(T1 >= d->T_t - 1e-4);
405          assert(T1 <= d->T_c);          assert(T1 <= d->T_c);
406            MSG("T1 = %f\n",T1);
407          res = fprops_sat_T(T1, &p, &rhof, &rhog, d);          res = fprops_sat_T(T1, &p, &rhof, &rhog, d);
408          if(res){          if(res){
409              ERRMSG("Failed to solve psat(T = %.12e) for %s",T1,d->name);              ERRMSG("Failed to solve psat(T = %.12e) for %s",T1,d->name);

Legend:
Removed from v.2300  
changed lines
  Added in v.2301

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