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

Contents of /trunk/models/johnpye/fprops/python/spinodal.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 2118 - (show annotations) (download) (as text)
Wed Dec 9 05:00:17 2009 UTC (10 years ago) by jpye
File MIME type: text/x-python
File size: 423 byte(s)
Fix build for Linux.
1 import fprops
2
3 D = fprops.helmholtz_data_carbondioxide;
4
5 from pylab import *
6 hold(1)
7
8 TT = array([217, 230,240,250,260,280,300],'float')
9 L = []
10
11 for T in TT:
12 rr = logspace(log10(10),log10(1200), 200)
13 print rr
14 pp = [max(fprops.helmholtz_p(T,r,D),1.) for r in rr]
15
16 semilogy(rr,pp)
17 L.append("T = %f K" % T)
18
19 xlabel('Density')
20 ylabel('Pressure')
21 #legend(L)
22
23 axis([10,1200,1.,1e6 * D.p_c])
24 #axis([0,1000,0,100e6])
25 show()
26

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