/[ascend]/trunk/models/johnpye/fprops/python/solve_ph_array.py
ViewVC logotype

Diff of /trunk/models/johnpye/fprops/python/solve_ph_array.py

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

revision 2265 by jpye, Thu Aug 5 06:28:01 2010 UTC revision 2266 by jpye, Thu Aug 5 13:53:01 2010 UTC
# Line 3  from pylab import * Line 3  from pylab import *
3    
4  D = helmholtz_data_water;  D = helmholtz_data_water;
5    
6  TT = linspace(273.16, 1000, 40);  TT = linspace(273.16, 1000, 100);
7  rr = logspace(log10(0.01), log10(900), 40);  rr = logspace(log10(0.01), log10(900), 100);
8    
9  goodT = []  goodT = []
10  goodv = []  goodv = []
# Line 22  for T in TT: Line 22  for T in TT:
22    
23          res, T1, rho1 = fprops_solve_ph(p,h,0,D);          res, T1, rho1 = fprops_solve_ph(p,h,0,D);
24          if res:          if res:
25              print "   +++ BAD RESULT"              print "   +++ BAD RESULT T1 = %f, rho1 = %f" % (T1, rho1)
26              if not isnan(T) and not isnan(rho):              if not isnan(T) and not isnan(rho):
27                  badT.append(T); badv.append(1./rho)                  badT.append(T); badv.append(1./rho)
28          else:          else:
# Line 50  for T in TT: Line 50  for T in TT:
50    
51  semilogx(vf1,TT1,"b-")  semilogx(vf1,TT1,"b-")
52  semilogx(vg1,TT1,"b-")  semilogx(vg1,TT1,"b-")
53    xlabel("specific volume")
54    ylabel("temperature")
55    
56  show()  show()
57    

Legend:
Removed from v.2265  
changed lines
  Added in v.2266

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