OrthoHydrogen#
References#
Equation of State#
J.W. Leachman, R.T. Jacobsen, S.G. Penoncello, and E.W. Lemmon. Fundamental Equations of State for Parahydrogen, Normal Hydrogen, and Orthohydrogen. J. Phys. Chem. Ref. Data, 38(3):721–748, 2009. doi:10.1063/1.3160306.Aliases#
Orthohydrogen, orthohydrogen, ORTHOHYDROGEN, ORTHOHYD
Fluid Information#
Parameter, Value |
|
|---|---|
General |
|
Molar mass [kg/mol] |
0.00201594 |
CAS number |
1333-74-0o |
ASHRAE class |
UNKNOWN |
Formula |
N/A |
Acentric factor |
-0.21859198995573192 |
InChI |
N/A |
InChIKey |
N/A |
SMILES |
N/A |
ChemSpider ID |
-1 |
Limits |
|
Maximum temperature [K] |
1000.0 |
Maximum pressure [Pa] |
2000000000.0 |
Triple point |
|
Triple point temperature [K] |
14.008 |
Triple point pressure [Pa] |
7559.882341057205 |
Critical point |
|
Critical point temperature [K] |
33.21981462174088 |
Critical point density [kg/m3] |
31.132951755877347 |
Critical point density [mol/m3] |
15443.39204335315 |
Critical point pressure [Pa] |
1309791.9204848697 |
REFPROP Validation Data#
Note
This figure compares the results generated from CoolProp and those generated from REFPROP. They are all results obtained in the form \(Y(T,\rho)\), where \(Y\) is the parameter of interest and which for all EOS is a direct evaluation of the EOS
You can download the script that generated the following figure here: (link to script), right-click the link and then save as… or the equivalent in your browser. You can also download this figure as a PDF.
Consistency Plots#
The following figure shows all the flash routines that are available for this fluid. A red + is a failure of the flash routine, a black dot is a success. Hopefully you will only see black dots. The red curve is the maximum temperature curve, and the blue curve is the melting line if one is available for the fluid.
In this figure, we start off with a state point given by T,P and then we calculate each of the other possible output pairs in turn, and then try to re-calculate T,P from the new input pair. If we don’t arrive back at the original T,P values, there is a problem in the flash routine in CoolProp. For more information on how these figures were generated, see CoolProp.Plots.ConsistencyPlots
Note
You can download the script that generated the following figure here: (link to script), right-click the link and then save as… or the equivalent in your browser. You can also download this figure as a PDF.
Flash consistency (HEOS): 173 inconsistent, 812 exceptions, 0 bad-phase across 11 input pair(s).
Download full failure list (CSV)
Failing state points (sample, up to 20 per pair/class)
Pair |
Class |
Region |
P [Pa] |
T [K] |
In1 |
Val1 |
In2 |
Val2 |
Error |
|---|---|---|---|---|---|---|---|---|---|
DmolarHmolar |
EXCEPTION |
1phase |
1.54741e+08 |
15.108 |
Dmolar |
55111.6 |
Hmolar |
3890.64 |
D < DLtriple 3890.64 3916.69 |
DmolarHmolar |
EXCEPTION |
1phase |
4.04008e+08 |
15.108 |
P |
4.04008e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001705,-0.967190] |
DmolarHmolar |
EXCEPTION |
1phase |
5.56311e+08 |
15.108 |
P |
5.56311e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001238,-0.976173] |
DmolarHmolar |
EXCEPTION |
1phase |
5.56311e+08 |
40.3616 |
Dmolar |
70056.3 |
Hmolar |
10295.5 |
D < DLtriple 10295.5 13121.2 |
DmolarHmolar |
EXCEPTION |
1phase |
7.6603e+08 |
15.108 |
P |
7.6603e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000899,-0.982696] |
DmolarHmolar |
EXCEPTION |
1phase |
7.6603e+08 |
40.3616 |
Dmolar |
74831.9 |
Hmolar |
13201.4 |
D < DLtriple 13201.4 17905.6 |
DmolarHmolar |
EXCEPTION |
1phase |
1.05481e+09 |
15.108 |
P |
1.05481e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000653,-0.987433] |
DmolarHmolar |
EXCEPTION |
1phase |
1.05481e+09 |
40.3616 |
Dmolar |
79965.8 |
Hmolar |
16988.4 |
D < DLtriple 16988.4 24254.6 |
DmolarHmolar |
EXCEPTION |
1phase |
1.05481e+09 |
65.6153 |
Dmolar |
79945.5 |
Hmolar |
16942.3 |
D < DLtriple 16942.3 24227 |
DmolarHmolar |
EXCEPTION |
1phase |
1.45245e+09 |
15.108 |
P |
1.45245e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000474,-0.990874] |
DmolarHmolar |
EXCEPTION |
1phase |
1.45245e+09 |
40.3616 |
Dmolar |
85468.9 |
Hmolar |
21912.8 |
D < DLtriple 21912.8 32607.4 |
DmolarHmolar |
EXCEPTION |
1phase |
1.45245e+09 |
65.6153 |
Dmolar |
85861.6 |
Hmolar |
21655.7 |
D < DLtriple 21655.7 33269.5 |
DmolarHmolar |
EXCEPTION |
1phase |
1.45245e+09 |
90.8689 |
Dmolar |
85176.2 |
Hmolar |
21893.3 |
D < DLtriple 21893.3 32120 |
DmolarHmolar |
EXCEPTION |
1phase |
2e+09 |
15.108 |
P |
2e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000344,-0.993372] |
DmolarHmolar |
EXCEPTION |
1phase |
2e+09 |
40.3616 |
Dmolar |
91348.6 |
Hmolar |
28304.2 |
D < DLtriple 28304.2 43509.9 |
DmolarHmolar |
EXCEPTION |
1phase |
2e+09 |
65.6153 |
Dmolar |
92130.7 |
Hmolar |
27766.5 |
D < DLtriple 27766.5 45126.6 |
DmolarHmolar |
EXCEPTION |
1phase |
2e+09 |
90.8689 |
Dmolar |
91717.3 |
Hmolar |
27878.5 |
D < DLtriple 27878.5 44267 |
DmolarHmolar |
EXCEPTION |
1phase |
2e+09 |
116.123 |
Dmolar |
90912.1 |
Hmolar |
28252 |
D < DLtriple 28252 42625 |
DmolarHmolar |
INCONSISTENT |
1phase |
5.56311e+08 |
90.8689 |
Dmolar |
67559.3 |
Hmolar |
10970.4 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
7.6603e+08 |
116.123 |
Dmolar |
71620.4 |
Hmolar |
14227.8 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
2e+09 |
242.391 |
Dmolar |
85995.2 |
Hmolar |
31505.8 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
1.45245e+09 |
191.883 |
Dmolar |
80731.6 |
Hmolar |
24171 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
4.04008e+08 |
65.6153 |
Dmolar |
64251.2 |
Hmolar |
8382.06 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
1.05481e+09 |
141.376 |
Dmolar |
76409.5 |
Hmolar |
18325.4 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
2e+09 |
217.137 |
Dmolar |
86967 |
Hmolar |
30751.3 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
1.45245e+09 |
166.63 |
Dmolar |
81873.1 |
Hmolar |
23495.2 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
2.93401e+08 |
40.3616 |
Dmolar |
61515.2 |
Hmolar |
6374.71 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
2e+09 |
141.376 |
Dmolar |
89962.8 |
Hmolar |
28763.2 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
7.6603e+08 |
90.8689 |
Dmolar |
73102.4 |
Hmolar |
13711.3 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
1.45245e+09 |
116.123 |
Dmolar |
84155.6 |
Hmolar |
22330.8 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
5.56311e+08 |
65.6153 |
Dmolar |
69151.4 |
Hmolar |
10521.5 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
1.05481e+09 |
116.123 |
Dmolar |
77725.8 |
Hmolar |
17756.9 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
2e+09 |
166.63 |
Dmolar |
88964.9 |
Hmolar |
29364.9 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
2e+09 |
191.883 |
Dmolar |
87959.9 |
Hmolar |
30033.2 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
1.05481e+09 |
90.8689 |
Dmolar |
78972.2 |
Hmolar |
17273.2 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
7.6603e+08 |
65.6153 |
Dmolar |
74378.1 |
Hmolar |
13312 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
4.04008e+08 |
40.3616 |
Dmolar |
65623.9 |
Hmolar |
8071.58 |
|
DmolarHmolar |
INCONSISTENT |
1phase |
1.45245e+09 |
141.376 |
Dmolar |
83026.3 |
Hmolar |
22876 |
|
DmolarP |
EXCEPTION |
1phase |
4.04008e+08 |
15.108 |
P |
4.04008e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001705,-0.967190] |
DmolarP |
EXCEPTION |
1phase |
5.56311e+08 |
15.108 |
P |
5.56311e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001238,-0.976173] |
DmolarP |
EXCEPTION |
1phase |
7.6603e+08 |
15.108 |
P |
7.6603e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000899,-0.982696] |
DmolarP |
EXCEPTION |
1phase |
1.05481e+09 |
15.108 |
P |
1.05481e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000653,-0.987433] |
DmolarP |
EXCEPTION |
1phase |
1.45245e+09 |
15.108 |
P |
1.45245e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000474,-0.990874] |
DmolarP |
EXCEPTION |
1phase |
2e+09 |
15.108 |
P |
2e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000344,-0.993372] |
DmolarP |
INCONSISTENT |
1phase |
2e+09 |
520.181 |
Dmolar |
76962.4 |
P |
2e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
2.13076e+08 |
15.108 |
Dmolar |
57791.7 |
P |
2.13076e+08 |
|
DmolarP |
INCONSISTENT |
1phase |
1.45245e+09 |
318.152 |
Dmolar |
75506 |
P |
1.45245e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
2e+09 |
217.137 |
Dmolar |
86967 |
P |
2e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
2.93401e+08 |
15.108 |
Dmolar |
60739.3 |
P |
2.93401e+08 |
|
DmolarP |
INCONSISTENT |
1phase |
2e+09 |
494.927 |
Dmolar |
77672.3 |
P |
2e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
1.45245e+09 |
141.376 |
Dmolar |
83026.3 |
P |
1.45245e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
1.05481e+09 |
191.883 |
Dmolar |
73817 |
P |
1.05481e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
2e+09 |
191.883 |
Dmolar |
87959.9 |
P |
2e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
1.05481e+09 |
267.644 |
Dmolar |
70249.6 |
P |
1.05481e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
1.12377e+08 |
15.108 |
Dmolar |
52680.8 |
P |
1.12377e+08 |
|
DmolarP |
INCONSISTENT |
1phase |
2e+09 |
469.674 |
Dmolar |
78401.4 |
P |
2e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
1.45245e+09 |
292.898 |
Dmolar |
76481 |
P |
1.45245e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
2e+09 |
444.42 |
Dmolar |
79150.6 |
P |
2e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
1.45245e+09 |
444.42 |
Dmolar |
71096.9 |
P |
1.45245e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
2e+09 |
898.985 |
Dmolar |
68092.8 |
P |
2e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
1.05481e+09 |
166.63 |
Dmolar |
75097 |
P |
1.05481e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
2e+09 |
570.688 |
Dmolar |
75596.2 |
P |
2e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
2e+09 |
772.717 |
Dmolar |
70735.5 |
P |
2e+09 |
|
DmolarP |
INCONSISTENT |
1phase |
2e+09 |
848.478 |
Dmolar |
69118.6 |
P |
2e+09 |
|
DmolarSmolar |
EXCEPTION |
1phase |
2.13076e+08 |
15.108 |
Dmolar |
57791.7 |
Smolar |
25.3638 |
D < DLtriple 25.3638 25.38 |
DmolarSmolar |
EXCEPTION |
1phase |
4.04008e+08 |
15.108 |
P |
4.04008e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001705,-0.967190] |
DmolarSmolar |
EXCEPTION |
1phase |
5.56311e+08 |
15.108 |
P |
5.56311e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001238,-0.976173] |
DmolarSmolar |
EXCEPTION |
1phase |
5.56311e+08 |
40.3616 |
Dmolar |
70056.3 |
Smolar |
26.7907 |
D < DLtriple 26.7907 58.1133 |
DmolarSmolar |
EXCEPTION |
1phase |
7.6603e+08 |
15.108 |
P |
7.6603e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000899,-0.982696] |
DmolarSmolar |
EXCEPTION |
1phase |
7.6603e+08 |
40.3616 |
Dmolar |
74831.9 |
Smolar |
27.1333 |
D < DLtriple 27.1333 83.931 |
DmolarSmolar |
EXCEPTION |
1phase |
7.6603e+08 |
65.6153 |
Dmolar |
74378.1 |
Smolar |
28.8987 |
D < DLtriple 28.8987 81.1012 |
DmolarSmolar |
EXCEPTION |
1phase |
1.05481e+09 |
15.108 |
P |
1.05481e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000653,-0.987433] |
DmolarSmolar |
EXCEPTION |
1phase |
1.05481e+09 |
40.3616 |
Dmolar |
79965.8 |
Smolar |
28.6133 |
D < DLtriple 28.6133 121.969 |
DmolarSmolar |
EXCEPTION |
1phase |
1.05481e+09 |
65.6153 |
Dmolar |
79945.5 |
Smolar |
27.2519 |
D < DLtriple 27.2519 121.796 |
DmolarSmolar |
EXCEPTION |
1phase |
1.05481e+09 |
90.8689 |
Dmolar |
78972.2 |
Smolar |
31.442 |
D < DLtriple 31.442 113.712 |
DmolarSmolar |
EXCEPTION |
1phase |
1.45245e+09 |
15.108 |
P |
1.45245e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000474,-0.990874] |
DmolarSmolar |
EXCEPTION |
1phase |
1.45245e+09 |
40.3616 |
Dmolar |
85468.9 |
Smolar |
31.6374 |
D < DLtriple 31.6374 176.113 |
DmolarSmolar |
EXCEPTION |
1phase |
1.45245e+09 |
65.6153 |
Dmolar |
85861.6 |
Smolar |
26.0692 |
D < DLtriple 26.0692 180.547 |
DmolarSmolar |
EXCEPTION |
1phase |
1.45245e+09 |
90.8689 |
Dmolar |
85176.2 |
Smolar |
29.032 |
D < DLtriple 29.032 172.859 |
DmolarSmolar |
EXCEPTION |
1phase |
1.45245e+09 |
116.123 |
Dmolar |
84155.6 |
Smolar |
33.2526 |
D < DLtriple 33.2526 161.852 |
DmolarSmolar |
EXCEPTION |
1phase |
2e+09 |
15.108 |
P |
2e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000344,-0.993372] |
DmolarSmolar |
EXCEPTION |
1phase |
2e+09 |
40.3616 |
Dmolar |
91348.6 |
Smolar |
36.7056 |
D < DLtriple 36.7056 250.999 |
DmolarSmolar |
EXCEPTION |
1phase |
2e+09 |
65.6153 |
Dmolar |
92130.7 |
Smolar |
25.5398 |
D < DLtriple 25.5398 262.38 |
DmolarSmolar |
EXCEPTION |
1phase |
2e+09 |
90.8689 |
Dmolar |
91717.3 |
Smolar |
26.8521 |
D < DLtriple 26.8521 256.322 |
DmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
621.195 |
Dmolar |
74296 |
Smolar |
78.9045 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
595.942 |
Dmolar |
74938.2 |
Smolar |
77.5869 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
570.688 |
Dmolar |
75596.2 |
Smolar |
76.2094 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
545.434 |
Dmolar |
76270.6 |
Smolar |
74.7664 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
520.181 |
Dmolar |
76962.4 |
Smolar |
73.2514 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
494.927 |
Dmolar |
77672.3 |
Smolar |
71.6569 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
469.674 |
Dmolar |
78401.4 |
Smolar |
69.9747 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
1.45245e+09 |
393.913 |
Dmolar |
72773.9 |
Smolar |
67.365 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
444.42 |
Dmolar |
79150.6 |
Smolar |
68.195 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
1.45245e+09 |
368.659 |
Dmolar |
73654.1 |
Smolar |
65.2498 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
419.166 |
Dmolar |
79921 |
Smolar |
66.307 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
1.45245e+09 |
343.405 |
Dmolar |
74564.2 |
Smolar |
62.9853 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
393.913 |
Dmolar |
80713.6 |
Smolar |
64.2981 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
1.45245e+09 |
318.152 |
Dmolar |
75506 |
Smolar |
60.5544 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
368.659 |
Dmolar |
81529.7 |
Smolar |
62.1547 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
1.05481e+09 |
267.644 |
Dmolar |
70249.6 |
Smolar |
58.2091 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
343.405 |
Dmolar |
82370.3 |
Smolar |
59.8612 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
1.45245e+09 |
292.898 |
Dmolar |
76481 |
Smolar |
57.939 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
1.05481e+09 |
242.391 |
Dmolar |
71391.6 |
Smolar |
55.1953 |
|
DmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
318.152 |
Dmolar |
83236.5 |
Smolar |
57.401 |
|
DmolarT |
EXCEPTION |
1phase |
4.04008e+08 |
15.108 |
P |
4.04008e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001705,-0.967190] |
DmolarT |
EXCEPTION |
1phase |
5.56311e+08 |
15.108 |
P |
5.56311e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001238,-0.976173] |
DmolarT |
EXCEPTION |
1phase |
7.6603e+08 |
15.108 |
P |
7.6603e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000899,-0.982696] |
DmolarT |
EXCEPTION |
1phase |
1.05481e+09 |
15.108 |
P |
1.05481e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000653,-0.987433] |
DmolarT |
EXCEPTION |
1phase |
1.45245e+09 |
15.108 |
P |
1.45245e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000474,-0.990874] |
DmolarT |
EXCEPTION |
1phase |
2e+09 |
15.108 |
P |
2e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000344,-0.993372] |
DmolarUmolar |
EXCEPTION |
1phase |
2.13076e+08 |
15.108 |
Dmolar |
57791.7 |
Umolar |
1263.93 |
D < DLtriple 1263.93 1264.15 |
DmolarUmolar |
EXCEPTION |
1phase |
4.04008e+08 |
15.108 |
P |
4.04008e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001705,-0.967190] |
DmolarUmolar |
EXCEPTION |
1phase |
5.56311e+08 |
15.108 |
P |
5.56311e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001238,-0.976173] |
DmolarUmolar |
EXCEPTION |
1phase |
5.56311e+08 |
40.3616 |
Dmolar |
70056.3 |
Umolar |
2354.58 |
D < DLtriple 2354.58 2977.31 |
DmolarUmolar |
EXCEPTION |
1phase |
7.6603e+08 |
15.108 |
P |
7.6603e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000899,-0.982696] |
DmolarUmolar |
EXCEPTION |
1phase |
7.6603e+08 |
40.3616 |
Dmolar |
74831.9 |
Umolar |
2964.68 |
D < DLtriple 2964.68 4122.65 |
DmolarUmolar |
EXCEPTION |
1phase |
1.05481e+09 |
15.108 |
P |
1.05481e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000653,-0.987433] |
DmolarUmolar |
EXCEPTION |
1phase |
1.05481e+09 |
40.3616 |
Dmolar |
79965.8 |
Umolar |
3797.64 |
D < DLtriple 3797.64 5720.87 |
DmolarUmolar |
EXCEPTION |
1phase |
1.05481e+09 |
65.6153 |
Dmolar |
79945.5 |
Umolar |
3748.17 |
D < DLtriple 3748.17 5713.77 |
DmolarUmolar |
EXCEPTION |
1phase |
1.45245e+09 |
15.108 |
P |
1.45245e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000474,-0.990874] |
DmolarUmolar |
EXCEPTION |
1phase |
1.45245e+09 |
40.3616 |
Dmolar |
85468.9 |
Umolar |
4918.84 |
D < DLtriple 4918.84 7909.19 |
DmolarUmolar |
EXCEPTION |
1phase |
1.45245e+09 |
65.6153 |
Dmolar |
85861.6 |
Umolar |
4739.53 |
D < DLtriple 4739.53 8085.64 |
DmolarUmolar |
EXCEPTION |
1phase |
1.45245e+09 |
90.8689 |
Dmolar |
85176.2 |
Umolar |
4841 |
D < DLtriple 4841 7779.51 |
DmolarUmolar |
EXCEPTION |
1phase |
2e+09 |
15.108 |
P |
2e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000344,-0.993372] |
DmolarUmolar |
EXCEPTION |
1phase |
2e+09 |
40.3616 |
Dmolar |
91348.6 |
Umolar |
6410.1 |
D < DLtriple 6410.1 10854.8 |
DmolarUmolar |
EXCEPTION |
1phase |
2e+09 |
65.6153 |
Dmolar |
92130.7 |
Umolar |
6058.24 |
D < DLtriple 6058.24 11297.6 |
DmolarUmolar |
EXCEPTION |
1phase |
2e+09 |
90.8689 |
Dmolar |
91717.3 |
Umolar |
6072.42 |
D < DLtriple 6072.42 11062 |
DmolarUmolar |
EXCEPTION |
1phase |
2e+09 |
116.123 |
Dmolar |
90912.1 |
Umolar |
6252.73 |
D < DLtriple 6252.73 10613.1 |
DmolarUmolar |
INCONSISTENT |
1phase |
1.05481e+09 |
141.376 |
Dmolar |
76409.5 |
Umolar |
4520.68 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
2.93401e+08 |
40.3616 |
Dmolar |
61515.2 |
Umolar |
1605.14 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
2e+09 |
217.137 |
Dmolar |
86967 |
Umolar |
7754.09 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
1.45245e+09 |
166.63 |
Dmolar |
81873.1 |
Umolar |
5754.87 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
7.6603e+08 |
90.8689 |
Dmolar |
73102.4 |
Umolar |
3232.46 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
2e+09 |
141.376 |
Dmolar |
89962.8 |
Umolar |
6531.79 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
5.56311e+08 |
65.6153 |
Dmolar |
69151.4 |
Umolar |
2476.67 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
1.45245e+09 |
116.123 |
Dmolar |
84155.6 |
Umolar |
5071.64 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
1.05481e+09 |
116.123 |
Dmolar |
77725.8 |
Umolar |
4185.96 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
2e+09 |
166.63 |
Dmolar |
88964.9 |
Umolar |
6884.17 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
2e+09 |
191.883 |
Dmolar |
87959.9 |
Umolar |
7295.58 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
1.05481e+09 |
90.8689 |
Dmolar |
78972.2 |
Umolar |
3916.48 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
4.04008e+08 |
40.3616 |
Dmolar |
65623.9 |
Umolar |
1915.17 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
7.6603e+08 |
65.6153 |
Dmolar |
74378.1 |
Umolar |
3012.85 |
|
DmolarUmolar |
INCONSISTENT |
1phase |
1.45245e+09 |
141.376 |
Dmolar |
83026.3 |
Umolar |
5382.07 |
|
HmolarP |
EXCEPTION |
1phase |
2.93401e+08 |
65.6153 |
Hmolar |
6747.09 |
P |
2.93401e+08 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: The temperature of -50.475369 K is below the minimum of 0.000000 K |
HmolarP |
EXCEPTION |
1phase |
4.04008e+08 |
15.108 |
P |
4.04008e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001705,-0.967190] |
HmolarP |
EXCEPTION |
1phase |
4.04008e+08 |
40.3616 |
Hmolar |
8071.58 |
P |
4.04008e+08 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.002698,-0.970283] |
HmolarP |
EXCEPTION |
1phase |
5.56311e+08 |
15.108 |
P |
5.56311e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001238,-0.976173] |
HmolarP |
EXCEPTION |
1phase |
5.56311e+08 |
40.3616 |
Hmolar |
10295.5 |
P |
5.56311e+08 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.001959,-0.978419] |
HmolarP |
EXCEPTION |
1phase |
7.6603e+08 |
15.108 |
P |
7.6603e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000899,-0.982696] |
HmolarP |
EXCEPTION |
1phase |
7.6603e+08 |
40.3616 |
Hmolar |
13201.4 |
P |
7.6603e+08 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.001423,-0.984327] |
HmolarP |
EXCEPTION |
1phase |
1.05481e+09 |
15.108 |
P |
1.05481e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000653,-0.987433] |
HmolarP |
EXCEPTION |
1phase |
1.05481e+09 |
40.3616 |
Hmolar |
16988.4 |
P |
1.05481e+09 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.001033,-0.988618] |
HmolarP |
EXCEPTION |
1phase |
1.45245e+09 |
15.108 |
P |
1.45245e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000474,-0.990874] |
HmolarP |
EXCEPTION |
1phase |
1.45245e+09 |
40.3616 |
Hmolar |
21912.8 |
P |
1.45245e+09 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.000750,-0.991734] |
HmolarP |
EXCEPTION |
1phase |
2e+09 |
15.108 |
P |
2e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000344,-0.993372] |
HmolarP |
EXCEPTION |
1phase |
2e+09 |
40.3616 |
Hmolar |
28304.2 |
P |
2e+09 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.000545,-0.993997] |
HmolarSmolar |
EXCEPTION |
1phase |
4.04008e+08 |
15.108 |
P |
4.04008e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001705,-0.967190] |
HmolarSmolar |
EXCEPTION |
1phase |
5.56311e+08 |
15.108 |
P |
5.56311e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001238,-0.976173] |
HmolarSmolar |
EXCEPTION |
1phase |
7.6603e+08 |
15.108 |
P |
7.6603e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000899,-0.982696] |
HmolarSmolar |
EXCEPTION |
1phase |
1.05481e+09 |
15.108 |
P |
1.05481e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000653,-0.987433] |
HmolarSmolar |
EXCEPTION |
1phase |
1.45245e+09 |
15.108 |
P |
1.45245e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000474,-0.990874] |
HmolarSmolar |
EXCEPTION |
1phase |
2e+09 |
15.108 |
P |
2e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000344,-0.993372] |
HmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
40.3616 |
Hmolar |
28304.2 |
Smolar |
36.7056 |
|
HmolarSmolar |
INCONSISTENT |
1phase |
1.45245e+09 |
40.3616 |
Hmolar |
21912.8 |
Smolar |
31.6374 |
|
HmolarSmolar |
INCONSISTENT |
1phase |
2.93401e+08 |
15.108 |
Hmolar |
6357.87 |
Smolar |
28.8472 |
|
HmolarSmolar |
INCONSISTENT |
1phase |
1.05481e+09 |
40.3616 |
Hmolar |
16988.4 |
Smolar |
28.6133 |
|
HmolarSmolar |
INCONSISTENT |
1phase |
2e+09 |
65.6153 |
Hmolar |
27766.5 |
Smolar |
25.5398 |
|
HmolarSmolar |
INCONSISTENT |
1phase |
7.6603e+08 |
40.3616 |
Hmolar |
13201.4 |
Smolar |
27.1333 |
|
PSmolar |
EXCEPTION |
1phase |
4.04008e+08 |
15.108 |
P |
4.04008e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001705,-0.967190] |
PSmolar |
EXCEPTION |
1phase |
4.04008e+08 |
40.3616 |
P |
4.04008e+08 |
Smolar |
27.2599 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.002698,-0.970283] |
PSmolar |
EXCEPTION |
1phase |
5.56311e+08 |
15.108 |
P |
5.56311e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001238,-0.976173] |
PSmolar |
EXCEPTION |
1phase |
5.56311e+08 |
40.3616 |
P |
5.56311e+08 |
Smolar |
26.7907 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.001959,-0.978419] |
PSmolar |
EXCEPTION |
1phase |
7.6603e+08 |
15.108 |
P |
7.6603e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000899,-0.982696] |
PSmolar |
EXCEPTION |
1phase |
7.6603e+08 |
40.3616 |
P |
7.6603e+08 |
Smolar |
27.1333 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.001423,-0.984327] |
PSmolar |
EXCEPTION |
1phase |
1.05481e+09 |
15.108 |
P |
1.05481e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000653,-0.987433] |
PSmolar |
EXCEPTION |
1phase |
1.05481e+09 |
40.3616 |
P |
1.05481e+09 |
Smolar |
28.6133 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.001033,-0.988618] |
PSmolar |
EXCEPTION |
1phase |
1.45245e+09 |
15.108 |
P |
1.45245e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000474,-0.990874] |
PSmolar |
EXCEPTION |
1phase |
1.45245e+09 |
40.3616 |
P |
1.45245e+09 |
Smolar |
31.6374 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.000750,-0.991734] |
PSmolar |
EXCEPTION |
1phase |
2e+09 |
15.108 |
P |
2e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000344,-0.993372] |
PSmolar |
EXCEPTION |
1phase |
2e+09 |
40.3616 |
P |
2e+09 |
Smolar |
36.7056 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.000545,-0.993997] |
PSmolar |
INCONSISTENT |
1phase |
2.93401e+08 |
40.3616 |
P |
2.93401e+08 |
Smolar |
28.2853 |
|
PT |
EXCEPTION |
1phase |
4.04008e+08 |
15.108 |
P |
4.04008e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001705,-0.967190] |
PT |
EXCEPTION |
1phase |
5.56311e+08 |
15.108 |
P |
5.56311e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001238,-0.976173] |
PT |
EXCEPTION |
1phase |
7.6603e+08 |
15.108 |
P |
7.6603e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000899,-0.982696] |
PT |
EXCEPTION |
1phase |
1.05481e+09 |
15.108 |
P |
1.05481e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000653,-0.987433] |
PT |
EXCEPTION |
1phase |
1.45245e+09 |
15.108 |
P |
1.45245e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000474,-0.990874] |
PT |
EXCEPTION |
1phase |
2e+09 |
15.108 |
P |
2e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000344,-0.993372] |
PUmolar |
EXCEPTION |
1phase |
4.04008e+08 |
15.108 |
P |
4.04008e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001705,-0.967190] |
PUmolar |
EXCEPTION |
1phase |
4.04008e+08 |
40.3616 |
P |
4.04008e+08 |
Umolar |
1915.17 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.002698,-0.970283] |
PUmolar |
EXCEPTION |
1phase |
5.56311e+08 |
15.108 |
P |
5.56311e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001238,-0.976173] |
PUmolar |
EXCEPTION |
1phase |
5.56311e+08 |
40.3616 |
P |
5.56311e+08 |
Umolar |
2354.58 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.001959,-0.978419] |
PUmolar |
EXCEPTION |
1phase |
7.6603e+08 |
15.108 |
P |
7.6603e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000899,-0.982696] |
PUmolar |
EXCEPTION |
1phase |
7.6603e+08 |
40.3616 |
P |
7.6603e+08 |
Umolar |
2964.68 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.001423,-0.984327] |
PUmolar |
EXCEPTION |
1phase |
1.05481e+09 |
15.108 |
P |
1.05481e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000653,-0.987433] |
PUmolar |
EXCEPTION |
1phase |
1.05481e+09 |
40.3616 |
P |
1.05481e+09 |
Umolar |
3797.64 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.001033,-0.988618] |
PUmolar |
EXCEPTION |
1phase |
1.45245e+09 |
15.108 |
P |
1.45245e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000474,-0.990874] |
PUmolar |
EXCEPTION |
1phase |
1.45245e+09 |
40.3616 |
P |
1.45245e+09 |
Umolar |
4918.84 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.000750,-0.991734] |
PUmolar |
EXCEPTION |
1phase |
2e+09 |
15.108 |
P |
2e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000344,-0.993372] |
PUmolar |
EXCEPTION |
1phase |
2e+09 |
40.3616 |
P |
2e+09 |
Umolar |
6410.1 |
unable to solve 1phase PY flash with Tmin=14.007, Tmax=1500 due to error: Inputs in Brent [37708.855316,41898.876838] do not bracket the root. Function values are [-1.000545,-0.993997] |
SmolarT |
EXCEPTION |
1phase |
4.04008e+08 |
15.108 |
P |
4.04008e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001705,-0.967190] |
SmolarT |
EXCEPTION |
1phase |
5.56311e+08 |
15.108 |
P |
5.56311e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.001238,-0.976173] |
SmolarT |
EXCEPTION |
1phase |
7.6603e+08 |
15.108 |
P |
7.6603e+08 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000899,-0.982696] |
SmolarT |
EXCEPTION |
1phase |
1.05481e+09 |
15.108 |
P |
1.05481e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000653,-0.987433] |
SmolarT |
EXCEPTION |
1phase |
1.45245e+09 |
15.108 |
P |
1.45245e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000474,-0.990874] |
SmolarT |
EXCEPTION |
1phase |
2e+09 |
15.108 |
P |
2e+09 |
T |
15.108 |
Inputs in Brent [37403.641030,41898.876838] do not bracket the root. Function values are [-1.000344,-0.993372] |
SmolarT |
EXCEPTION |
1phase |
2e+09 |
65.6153 |
Smolar |
25.5398 |
T |
65.6153 |
Even by increasing rhoc, not able to bound input; input 25.5398 is not in range 1643.73,339.177 |
SmolarT |
INCONSISTENT |
1phase |
2.93401e+08 |
15.108 |
Smolar |
28.8472 |
T |
15.108 |
|
SmolarT |
INCONSISTENT |
1phase |
2e+09 |
40.3616 |
Smolar |
36.7056 |
T |
40.3616 |
|
SmolarT |
INCONSISTENT |
1phase |
1.45245e+09 |
40.3616 |
Smolar |
31.6374 |
T |
40.3616 |
|
SmolarT |
INCONSISTENT |
1phase |
2.13076e+08 |
15.108 |
Smolar |
25.3638 |
T |
15.108 |
|
SmolarT |
INCONSISTENT |
1phase |
1.05481e+09 |
40.3616 |
Smolar |
28.6133 |
T |
40.3616 |
|
SmolarT |
INCONSISTENT |
1phase |
1.54741e+08 |
15.108 |
Smolar |
23.5365 |
T |
15.108 |
|
SmolarT |
INCONSISTENT |
1phase |
7.6603e+08 |
40.3616 |
Smolar |
27.1333 |
T |
40.3616 |
|
SmolarT |
INCONSISTENT |
1phase |
1.12377e+08 |
15.108 |
Smolar |
22.8137 |
T |
15.108 |
Superancillary Plots#
The following figure shows the accuracy of the superancillary functions relative to extended precision calculations carried out in C++ with the teqp library. The results of the iterative calculations with REFPROP and CoolProp are also shown.
Note
You can download the script that generated the following figure here: (link to script), right-click the link and then save as… or the equivalent in your browser. You can also download this figure as a PDF.
