contour_hfoucolor.ppi

# produces a contour plot for the Fe-H2O system with FeT = 0.01 mol/kgw and Fe(OH)3(a) as a possible mineral phase
# select contours from user-defined set of values

SPECIATION
  calculationType                      "contour"
  calculationMethod                    1
  contourZvariable                     FeT                                     # plot is for FeT defined in USER_PUNCH
  xmin                                 2.0                                     # pH range
  xmax                                 12.0
  ymin                                 -90.0                                   # fO2(g) range (controls the redox)
  ymax                                 0.0
  resolution                           50                                      # calculate a 50 x 50 grid

PLOT
  plotTitle                            "Fe solubility<br> user contours and colours"
  xtitle                               pH
  ytitle                               "log <i>f</i> O<sub>2</sub>(g) (atm)"

  contours                             -7  -6  -5  -4  -3 -2.5 -2.1 -2.0                            # n user-defined levels
  contourFillcolor                     blue3 sky3 cyan3 spring3  green3 lawn3 yellow3 orange3 red3  # n+1 user-defined colours
  contourLineWidth                     auto
  contourLineColor                     auto
  contourShiftLabel                    c
  contourLabelSize                     2
  contourLabelFont                     "Helvetica"
  contourLabelColor                    auto
  contourlabelFigs                     auto

  extraText                            extratextcontour_hfo.dat

CHEMISTRY

# first simulation
PHASES
Fix_H+; H+ = H+; log_k 0.

SELECTED_OUTPUT
-reset FALSE
-high_precision TRUE

USER_PUNCH
-headings  FeT
10 PUNCH log10(TOT("Fe"))

SOLUTION 1
   pH        1.8
   units     mol/kgw
   Fe(3)     1e-2
   Na        1e-1
   Cl        1e-1
END

# second simulation - loop on this
USE solution 1
EQUILIBRIUM_PHASES 1
   Fix_H+ -<x_axis> NaOH 10          # N.B. this works because <x_axis> is substituted without any leading spaces
     -force_equality true            #   for negative values of <x_axis> would have to form a new tag to avoid --value
   O2(g)  <y_axis> 0.1
   Fe(OH)3(a) 0 0
END