GMT script to generate Figure 4 of Park and Royer (2011)

#!/bin/csh

gmtset ANNOT_FONT_SIZE 12p LABEL_FONT_SIZE 15p HEADER_FONT_SIZE 15p

gmtset ANNOT_FONT Times-Roman LABEL_FONT Times-Roman HEADER_FONT Times-Roman

awk '{print $1, $2, $3/1000}' out_ppdf_control.dat >! blob.dat

surface blob.dat -R0.2/10.6/0.025/0.135 -I0.4/0.01 -Gblobg -T0.25

grdimage blobg -JX3/2 -Cppdf.cpt -R -Ba2f1:'CO@-2@- Sensitivity (degC)':/0.02:'ACT'::'.Fraction of ModelRuns within 3-@~s@~':/SWen -X1.0 -Y1.0 -K >! aaplot

awk '{print $1, $5, $6/1000}' out_ppdf_control.dat >! blob.dat

surface blob.dat -R0.2/10.6/0.15/0.85 -I0.4/0.05 -Gblobg -T0.25

grdimage blobg -JX3/2 -Cppdf.cpt -R -Ba2f1:'CO@-2@- Sensitivity (degC)':/0.1:'FERT'::'.Fraction of ModelRuns within 3-@~s@~':/SWen -X0.0 -Y3.5 -K -O >> aaplot

awk '{print $1, $8, $9/1000}' out_ppdf_control.dat >! blob.dat

surface blob.dat -R0.2/10.6/0.05/0.55 -I0.4/0.025 -Gblobg -T0.25

grdimage blobg -JX3/2 -Cppdf.cpt -R -Ba2f1:'CO@-2@- Sensitivity (degC)':/0.1:'LIFE'::'.Fraction of ModelRuns within 3-@~s@~':/SWen -X4.0 -Y0.0 -K -O >> aaplot

awk '{print $1, $11, $12/1000}' out_ppdf_control.dat >! blob.dat

surface blob.dat -R0.2/10.6/0.45/1.25 -I0.4/0.05 -Gblobg -T0.25

grdimage blobg -JX3/2 -Cppdf.cpt -R -Ba2f1:'CO@-2@- Sensitivity (degC)':/0.2:'GYM'::'.Fraction of ModelRuns within 3-@~s@~':/SWen -X0.0 -Y-3.5 -K -O >> aaplot

psscale -D3.5/1/2/0.25c -Cppdf.cpt -Ba0.2f0.1:'.Fraction <3-@~s@~': -O >> aaplot

gv aaplot &