! plots atmospheric CO2, freshwater forcing, net primary production, Atlantic meridional overturning circulation and limitation factors for PO4 and Fe


cancel data/all
cancel window/all
set memory/size=200

set window/quality=high/thicken=1.2/aspect=1.25:axis 1


!----------------SIMULATIONSDATEN LADEN----------------------------
use "/p/projects/climber3/hofmann/RL_NEW/CO2X4/snap_ocean_1000.nc" 1 
use "/p/projects/climber3/hofmann/RL_NEW/CO2X4/history.nc" 2
use "/p/projects/climber3/hofmann/RL_NEW/CO2X2/snap_ocean.nc" 3
use "/p/projects/climber3/hofmann/RL_NEW/CO2X2/history.nc" 4
use "/p/projects/climber3/hofmann/RL_NEW/CO2X4_HOSING/snap_ocean_hosing_1000.nc" 5 
use "/p/projects/climber3/hofmann/RL_NEW/CO2X4_HOSING/history.nc" 6
use "/p/projects/climber3/hofmann/RL_NEW/CONTROL/snap_ocean_1000.nc" 7
use "/p/projects/climber3/hofmann/RL_NEW/CONTROL/history.nc" 8
use "/p/projects/climber3/hofmann/RL_NEW/CO2X2_HOSING/snap_ocean_hosing.nc" 9 
use "/p/projects/climber3/hofmann/RL_NEW/CO2X2_HOSING/history.nc" 10

!----------------GRAFIKEINSTELLUNGEN-------------------------------
! for 6 plots on a page (2 columns)
define viewport/xlimits=0,0.55/ylim=0.05,0.35 leftlow
define viewport/xlimits=0.45,0.99/ylim=0.05,0.35 rightlow
define viewport/xlimits=0,0.55/ylim=0.36,0.66 leftmid
define viewport/xlimits=0.45,0.99/ylim=0.36,0.66 rightmid
define viewport/xlimits=0,0.55/ylim=0.67,0.97 leftup
define viewport/xlimits=0.45,0.99/ylim=0.67,0.97 rightup

!----------------FORMELN---------------------------------
!
let argT1 =  32.0 - temp[d=1]
let flimit1 = 1.e-3 * (max(0.08,argT1/(argT1 + 20.0)))
let fvalue1 = tracer_13[d=1]/(tracer_13[d=1] + flimit1)
!
let plimit1 = tracer_06[d=1]/((1./122.)+tracer_06[d=1])
!
let argT3 =  32.0 - temp[d=3]
let flimit3 = 1.e-3 * (max(0.08,argT3/(argT3 + 20.0)))
let fvalue3 = tracer_13[d=3]/(tracer_13[d=3] + flimit3)
!
let plimit3 = tracer_06[d=3]/((1./122.)+tracer_06[d=3])
!
let argT5 =  32.0 - temp[d=5]
let flimit5 = 1.e-3 * (max(0.08,argT5/(argT5 + 20.0)))
let fvalue5 = tracer_13[d=5]/(tracer_13[d=5] + flimit5)
!
let plimit5 = tracer_06[d=5]/((1./122.)+tracer_06[d=5])
!
let argT7 =  32.0 - temp[d=7]
let flimit7 = 1.e-3 * (max(0.08,argT7/(argT7 + 20.0)))
let fvalue7 = tracer_13[d=7]/(tracer_13[d=7] + flimit7)
!
let plimit7 = tracer_06[d=7]/((1./122.)+tracer_06[d=7])
!
let argT9 =  32.0 - temp[d=9]
let flimit9 = 1.e-3 * (max(0.08,argT9/(argT9 + 20.0)))
let fvalue9 = tracer_13[d=9]/(tracer_13[d=9] + flimit9)
!
let plimit9 = tracer_06[d=9]/((1./122.)+tracer_06[d=9])
!
!----------------ZEITACHSE----------------------------------------
!
define axis/units=years/t=-50:1000:1 nt1
define grid/T=nt1 ntime1
!
let tt=T[gt=ntime1@asn]
let t0=tt[T="01-jan-1901"]
let t1=tt[T="01-jan-2040"]
let t2=tt[T="01-jan-2180"]
let t3=tt[T="01-jan-1970"]
!
!----------------GRAFISCHE AUSGABE---------------------------------

SET TEXT/FONT=symbol
set window/xpix=1680/ypix=1920 1
!
set view leftup
GO margins  0.5 1.0 2 1.0
ppl labset, 0.12,0.12,0.12,0.4
ppl axlsze 0.14,0.14
ppl axset 0,1,1,0
plot/nolabels/thickness=2/color=red/vlimits=200:1250/hlimits=-50:1000:100 pco2diag[gt=ntime1@asn,d=2]
plot/over/nolabels/thickness=2/color=(100,55,0) pco2diag[gt=ntime1@asn,d=4]
plot/over/nolabels/thickness=2/color=black pco2diag[gt=ntime1@asn,d=8]
let tt=T[gt=ntime1@asn]
let t0=tt[T="01-jan-1901"]
let t1=tt[T="01-jan-2040"]
let t2=tt[T="01-jan-2180"]
let t3=tt[T="01-jan-1970"]
plot/over/vs/line/color=black/nolabel {`T0`,`T0`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T1`,`T1`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T2`,`T2`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T3`,`T3`},{($YAXIS_MIN),($YAXIS_MAX)}
label/nouser 4.3,2.9,0,0,.2 "a)"
label/nouser -0.9,1.5,0,90,.2 "pCO<sub>2</sub> (ppmv)"
!
set view rightup
GO margins  0.5 1.0 2 1.0
ppl labset, 0.12,0.12,0.12,0.4
ppl axlsze 0.14,0.14
ppl axset 0,1,1,0
let FWT4 = (0.05/4.) * (pco2diag[gt=ntime1@asn,d=6]/280.)
let val4 = if (ll GT 50 AND ll LT 330) then FWT4 else 0
let FWT2 = (0.05/4.) * (pco2diag[gt=ntime1@asn,d=10]/280.)
let val2 = if (ll GT 50 AND ll LT 190) then FWT2 else 0
plot/nolabels/thickness=2/color=blue/hlimits=-50:1000:100 val4 ! [l=1:500]
plot/over/nolabels/thickness=2/color=LightBlue val2 ! [l=1:500]
plot/over/vs/line/color=black/nolabel {`T0`,`T0`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T1`,`T1`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T2`,`T2`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T3`,`T3`},{($YAXIS_MIN),($YAXIS_MAX)}
label/nouser 4.3,2.9,0,0,.2 "@b)"
label/nouser -1.0,1.5,0,90,.15 "NA freshw. add. F<sub>FWT</sub> (Sv)"
!
set view leftmid
GO margins  0.5 1.0 2 1.0
ppl labset, 0.12,0.12,0.12,0.4
ppl axlsze 0.14,0.14
ppl axset 0,1,1,0
plot/nolabels/thickness=2/color=red/hlimits=-50:1000:100/vlimits=0:24 atl_over[gt=ntime1@asn,y=30n:60n@max,z=500:1500@max,d=2]
plot/over/nolabels/thickness=2/color=(100,55,0) atl_over[gt=ntime1@asn,y=30n:60n@max,z=500:1500@max,d=4]
plot/over/nolabels/thickness=2/color=blue atl_over[gt=ntime1@asn,y=30n:60n@max,z=500:1500@max,d=6]
plot/over/nolabels/thickness=2/color=black atl_over[gt=ntime1@asn,y=30n:60n@max,z=500:1500@max,d=8]
plot/over/nolabels/thickness=2/color=lightblue atl_over[gt=ntime1@asn,y=30n:60n@max,z=500:1500@max,d=10]
plot/over/vs/line/color=black/nolabel {`T0`,`T0`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T1`,`T1`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T2`,`T2`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T3`,`T3`},{($YAXIS_MIN),($YAXIS_MAX)}
label/nouser -0.9,1.7,0,90,.175 "Atlantic Overturning (Sv)"
label/nouser 4.3,2.9,0,0,.2 "c)"
!
set view rightmid
GO margins  0.5 1.0 2 1.0
ppl labset, 0.12,0.12,0.12,0.4
ppl axlsze 0.14,0.14
ppl axset 0,1,1,0
plot/nolabels/thickness=2/color=red/hlimits=-50:1000:100 glb_xpp[gt=ntime1@asn,i=@din,j=@din,d=2]*1.e-15
plot/over/nolabels/thickness=2/color=(100,55,0) glb_xpp[gt=ntime1@asn,i=@din,j=@din,d=4]*1.e-15
plot/over/nolabels/thickness=2/color=blue glb_xpp[gt=ntime1@asn,i=@din,j=@din,d=6]*1.e-15
plot/over/nolabels/thickness=2/color=black glb_xpp[gt=ntime1@asn,i=@din,j=@din,d=8]*1.e-15
plot/over/nolabels/thickness=2/color=lightblue glb_xpp[gt=ntime1@asn,i=@din,j=@din,d=10]*1.e-15
plot/over/vs/line/color=black/nolabel {`T0`,`T0`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T1`,`T1`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T2`,`T2`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T3`,`T3`},{($YAXIS_MIN),($YAXIS_MAX)}
label/nouser -0.9,1.5,0,90,.2 "NPP (GtC/a)"
label/nouser 4.3,2.9,0,0,.2 "d)"
!
set view leftlow
GO margins  0.5 1.0 2 1.0
ppl labset, 0.12,0.12,0.12,0.4
ppl axlsze 0.14,0.14
!plot plimit1[x=@ave,y=@ave,d=1,z=100]
plot/nolabels/thickness=2/color=red/vlimits=0.88:0.98:0.01/hlimits=-50:1000:100 plimit1[gt=ntime1@asn,i=@ave,j=@ave,z=100]
plot/over/nolabels/thickness=2/color=(100,55,0) plimit3[gt=ntime1@asn,i=@ave,j=@ave,z=100]
plot/over/nolabels/thickness=2/color=blue plimit5[gt=ntime1@asn,i=@ave,j=@ave,z=100]
plot/over/nolabels/thickness=2/color=black plimit7[gt=ntime1@asn,i=@ave,j=@ave,z=100]
plot/over/nolabels/thickness=2/color=lightblue plimit9[gt=ntime2@asn,i=@ave,j=@ave,z=100]
plot/over/vs/line/color=black/nolabel {`T0`,`T0`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T1`,`T1`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T2`,`T2`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T3`,`T3`},{($YAXIS_MIN),($YAXIS_MAX)}
label/nouser 4.3,2.9,0,0,.2 "e)"
label/nouser -0.9,1.5,0,90,.2 "P availability L<sub>PO4</sub>"
label/nouser 2.8,-0.8,0,0,.2 "Year of simulation"
!
set view rightlow
GO margins  0.5 1.0 2 1.0
ppl labset, 0.12,0.12,0.12,0.4
ppl axlsze 0.14,0.14
!plot fvalue1[x=@ave,y=@ave,d=1,z=100]
plot/nolabels/thickness=2/color=red/vlimits=0.47:0.62:0.01/hlimits=-50:1000:100 fvalue1[gt=ntime1@asn,i=@ave,j=@ave,z=100]
plot/over/nolabels/thickness=2/color=(100,55,0) fvalue3[gt=ntime1@asn,i=@ave,j=@ave,z=100]
plot/over/nolabels/thickness=2/color=blue fvalue5[gt=ntime1@asn,i=@ave,j=@ave,z=100]
plot/over/nolabels/thickness=2/color=black fvalue7[gt=ntime1@asn,i=@ave,j=@ave,z=100]
plot/over/nolabels/thickness=2/color=lightblue fvalue9[gt=ntime1@asn,i=@ave,j=@ave,z=100]
plot/over/vs/line/color=black/nolabel {`T0`,`T0`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T1`,`T1`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T2`,`T2`},{($YAXIS_MIN),($YAXIS_MAX)}
plot/over/vs/line/color=black/nolabel {`T3`,`T3`},{($YAXIS_MIN),($YAXIS_MAX)}
label/nouser 4.3,2.9,0,0,.2 "@CRf)"
label/nouser -0.9,1.5,0,90,.2 "@CR Fe availability L<sub>Fe</sub>"
label/nouser 2.8,-0.8,0,0,.2 "@CR Year of simulation"
!
define view/x=0.35:0.45/y=0.77:0.87/axes legend
set view legend
plot/vs/sym=dot/nolab/noax/vlim=0:1/hlim=0:1/color=white {0,1}, {0,1}
plot/over/vs/nolab/line=1/thickness=3 {0.1,0.8}, {0.95,0.95}
annotate/norm/xpos=.1/ypos=0.85/halign=-1/valign=0/size=0.2 "CONTROL"
plot/over/vs/nolab/line/color=(100,55,0) {0.1,0.8}, {0.75,0.75}
annotate/norm/xpos=.1/ypos=0.65/halign=-1/valign=0/size=0.2 "2xCO2"
plot/over/vs/nolab/line=2/thickness=3 {0.1,0.8}, {0.55,0.55}
annotate/norm/xpos=.1/ypos=0.45/halign=-1/valign=0/size=0.2 "4xCO2"
plot/over/vs/nolab/line=5/thickness=3 {0.1,0.8}, {0.35,0.35}
annotate/norm/xpos=.1/ypos=0.25/halign=-1/valign=0/size=0.2 "2xCO2HOSING"
plot/over/vs/nolab/line=4/thickness=3 {0.1,0.8}, {0.15,0.15}
annotate/norm/xpos=.1/ypos=0.05/halign=-1/valign=0/size=0.2 "4xCO2HOSING"
!frame/file="figures/COMFORT/paper/COMFORT_means_forcings-limits-NPP.pdf"