TITLE Example 22.--Compare experimental CO2 solubilities at high CO2 pressure with
      Peng-Robinson calc`ns with fixed-volume gas_phase, 25, 50, 75, 100, 150 oC.
#PRINT; -reset false
SOLUTION 1
GAS_PHASE 1
    -fixed_volume
    CO2(g) 0
    H2O(g) 0
REACTION
 CO2 1; 0 31*1
INCREMENTAL_REACTIONS true

USER_GRAPH 1 Example 22,
  -chart_title "CO2 solubility at high pressures, 25 - 100 ºC"
 -plot_csv_file co2.dat
 -axis_titles "Pressure / atm" "CO2 / (mol/kgw)"
 -axis_scale x_axis 0 500
 -axis_scale y_axis 0 2
 -connect_simulations false
 10 graph_x PR_P("CO2(g)")
 20 graph_y TOT("C(4)")
 -end
USER_GRAPH 2 Example 22, P-Vm of CO2 gas, 25 - 150 ºC
 -chart_title "P-Vm of CO2 gas, data from Michels et al., 1935, 1937"
 -plot_csv_file co2_VP.dat
 -axis_titles "molar volume of CO2 gas / (L/mol)" "CO2 pressure / atm"
 -axis_scale x_axis 0 0.7
 -axis_scale y_axis 0 700
 -connect_simulations false
 10 graph_x gas_vm : graph_y gas_p
END

USE solution 1
USE gas_phase 1
REACTION
 CO2 1; 0 29*1
REACTION_TEMPERATURE 2
 50
END

USE solution 1
USE gas_phase 1
REACTION
 CO2 1; 0 27*1
REACTION_TEMPERATURE 3
 75
USER_GRAPH 2; -active false
END

USE solution 1
USE gas_phase 1
REACTION
 CO2 1; 0 28*1
REACTION_TEMPERATURE 4
 100
USER_GRAPH 2; -active true
END

USE solution 1
USE gas_phase 1
REACTION
 CO2 1; 0 26*1
REACTION_TEMPERATURE 5
 150
USER_GRAPH 1; -active false
USER_GRAPH 2; -active true
END