| 1 | real*8 function psvco2(T) |
|---|
| 2 | ! solid-vapor transition for CO2 |
|---|
| 3 | ! returns saturation pressure in Pascal, input is temperature in Kelvin |
|---|
| 4 | implicit none |
|---|
| 5 | real*8 T |
|---|
| 6 | |
|---|
| 7 | psvco2 = exp(23.3494 - 3182.48/T)*100. |
|---|
| 8 | end |
|---|
| 9 | |
|---|
| 10 | |
|---|
| 11 | real*8 function tfrostco2(p) |
|---|
| 12 | ! the inverse of function psvco2 |
|---|
| 13 | ! input is pressure in Pascal, output is temperature in Kelvin |
|---|
| 14 | implicit none |
|---|
| 15 | real*8 p |
|---|
| 16 | tfrostco2 = 3182.48/(23.3494+log(100./p)) |
|---|
| 17 | end |
|---|
| 18 | |
|---|
| 19 | |
|---|
| 20 | |
|---|
| 21 | |
|---|
| 22 | !------------------------------------------------------------------------------ |
|---|
| 23 | ! Antoine equation parameters from NIST, 154K-196K |
|---|
| 24 | ! based on Giauque and Egan (1937) |
|---|
| 25 | ! A=6.81228, B=1301.679, C=-3.494 |
|---|
| 26 | ! p = 10**(A-(B/(T+C)))*1.e5 |
|---|
| 27 | |
|---|
| 28 | ! Expressions from Int. Crit. Tabl. 3.207, based on many references |
|---|
| 29 | ! mm2Pa = 133.32 |
|---|
| 30 | ! -135C to -56.7C (138K-216K) |
|---|
| 31 | ! p = 10**(-0.05223*26179.3/T + 9.9082)*mm2Pa |
|---|
| 32 | ! -183C to -135C (90K-138K) |
|---|
| 33 | ! p = 10**(-1275.62/T + 0.006833*T + 8.3071)*mm2Pa |
|---|
| 34 | ! Expressions from Int. Crit. Tabl. 3.208, based on Henning |
|---|
| 35 | ! -110 to -80C (163K-193K) |
|---|
| 36 | ! p = 10**(- 1279.11/T + 1.75*log10(T) - 0.0020757*T + 5.85242)*mm2Pa |
|---|
| 37 | ! p = 10**(- 1352.6/T + 9.8318)*mm2Pa |
|---|
| 38 | |
|---|
| 39 | ! Mars book (1992), p959, fit by chapter authors |
|---|
| 40 | ! p = exp(23.3494 - 3182.48/T)*100. ! 120-160K |
|---|
| 41 | ! p = exp(25.2194 - 3311.57/T - 6.71e-3*T)*100 ! 100-170K |
|---|
| 42 | ! Mars book (1992), p960, based on Miller & Smythe (1970) |
|---|
| 43 | ! p = exp(26.1228 - 3385.26/T - 9.4461e-3*T)*100 ! 123-173K |
|---|
| 44 | |
|---|
| 45 | ! Fray & Schmitt, PSS 57, 2053 (2009) |
|---|
| 46 | ! A0=1.476e1, A1=-2.571e3, A2=-7.781e4, A3=4.325e6, A4=-1.207e8, A5=1.350e9 |
|---|
| 47 | ! p = exp(A0+A1/T+A2/T**2+A3/T**3+A4/T**4+A5/T**5)*1e5 ! 40K-194.7K |
|---|
| 48 | ! A0=1.861e1, A1=-4.154e3, A2=1.041e5 |
|---|
| 49 | ! p = exp(A0 + A1/T + A2/T**2)*1e5 ! 194.7K-216.58K |
|---|
| 50 | !------------------------------------------------------------------------------ |
|---|
