! plots of mixed layer height and net primary production
! comparing 5 scenarios (CONTROL, 2xCO2, 2xCO2hosing, 4xCO2, 4xCO2hosing) and 3 regions (East equatorial Pacific, Southern Ocean and Indian Ocean)

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

use "/p/projects/climber3/hofmann/RL_NEW/CO2X4/snap_hmxl_1000.nc" 11
use "/p/projects/climber3/hofmann/RL_NEW/CO2X2/snap_hmxl.nc" 12
use "/p/projects/climber3/hofmann/RL_NEW/CO2X4_HOSING/snap_hmxl_hosing_1000.nc" 13
use "/p/projects/climber3/hofmann/RL_NEW/CONTROL/snap_hmxl_1000.nc" 14
use "/p/projects/climber3/hofmann/RL_NEW/CO2X2_HOSING/snap_hmxl_hosing.nc" 15

!----------------GRAFIKEINSTELLUNGEN-------------------------------

! for 6 plots on a page (3 columns)
define viewport/xlimits=0,0.39/ylim=0.08,0.53 leftlow3
define viewport/xlimits=0,0.39/ylim=0.54,0.99 leftup3
define viewport/xlimits=0.285,0.675/ylim=0.08,0.53 midlow3
define viewport/xlimits=0.285,0.675/ylim=0.54,0.99 midup3
define viewport/xlimits=0.57,0.96/ylim=0.08,0.53 rightlow3
define viewport/xlimits=0.57,0.96/ylim=0.54,0.99 rightup3

!----------------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=1920/ypix=1080 1
!
set view leftup3
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/nolabels/thickness=2/color=red/hlimits=-50:1000:100/vlimits=7:11.5 glb_xpp[gt=ntime1@asn,x=180w:80w@din,y=20s:20n@din,d=2]*1.e-15
plot/over/nolabels/thickness=2/color=(100,55,0) glb_xpp[gt=ntime1@asn,x=180w:80w@din,y=20s:20n@din,d=4]*1.e-15
plot/over/nolabels/thickness=2/color=blue glb_xpp[gt=ntime1@asn,x=180w:80w@din,y=20s:20n@din,d=6]*1.e-15
plot/over/nolabels/thickness=2/color=black glb_xpp[gt=ntime1@asn,x=180w:80w@din,y=20s:20n@din,d=8]*1.e-15
plot/over/nolabels/thickness=2/color=lightblue glb_xpp[gt=ntime1@asn,x=180w:80w@din,y=20s:20n@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 4.2,3.1,0,0,.2 "@CRa)"
label/nouser 2.2,-0.8,0,0,.2 "@CR EEP"
label/nouser -1,1.6,0,90,.2 "@CRNPP (GtC/a)"
!
set view midup3
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/nolabels/thickness=2/color=red/hlimits=-50:1000:100/vlimits=2.4:5.0 glb_xpp[gt=ntime1@asn,x=20e:100e@din,y=30n:40s@din,d=2]*1.e-15
plot/over/nolabels/thickness=2/color=(100,55,0) glb_xpp[gt=ntime1@asn,x=20e:100e@din,y=30n:40s@din,d=4]*1.e-15
plot/over/nolabels/thickness=2/color=blue glb_xpp[gt=ntime1@asn,x=20e:100e@din,y=30n:40s@din,d=6]*1.e-15
plot/over/nolabels/thickness=2/color=black glb_xpp[gt=ntime1@asn,x=20e:100e@din,y=30n:40s@din,d=8]*1.e-15
plot/over/nolabels/thickness=2/color=lightblue glb_xpp[gt=ntime1@asn,x=20e:100e@din,y=30n:40s@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 4.2,3.1,0,0,.2 "@CRb)"
label/nouser 2.2,-0.8,0,0,.2 "@CRIndian Ocean"
!
set view rightup3
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/nolabels/thickness=2/color=red/hlimits=-50:1000:100/y=40s:80s/vlimits=13.7:18 glb_xpp[gt=ntime1@asn,i=@din,j=@din,d=2]*1.e-15
plot/over/nolabels/thickness=2/color=(100,55,0)/y=40s:80s glb_xpp[gt=ntime1@asn,i=@din,j=@din,d=4]*1.e-15
plot/over/nolabels/thickness=2/color=blue/y=40s:80s glb_xpp[gt=ntime1@asn,i=@din,j=@din,d=6]*1.e-15
plot/over/nolabels/thickness=2/color=black/y=40s:80s glb_xpp[gt=ntime1@asn,i=@din,j=@din,d=8]*1.e-15
plot/over/nolabels/thickness=2/color=lightblue/y=40s:80s 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 2.2,-0.8,0,0,.2 "@CRSouthern Ocean"
label/nouser 4.2,3.1,0,0,.2 "@CRc)"
!
set view leftlow3
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/nolabels/thickness=2/color=red/hlimits=-50:1000:100/vlimits=52.5:57.5 hmxl[gt=ntime1@asn,x=180w:80w@ave,y=20s:20n@ave,d=11]/100 !/vlimits=57:64
plot/over/nolabels/thickness=2/color=(100,55,0) hmxl[gt=ntime1@asn,x=180w:80w@ave,y=20s:20n@ave,d=12]/100
plot/over/nolabels/thickness=2/color=blue hmxl[gt=ntime2@asn,x=180w:80w@ave,y=20s:20n@ave,d=13]/100
plot/over/nolabels/thickness=2/color=black hmxl[gt=ntime1@asn,x=180w:80w@ave,y=20s:20n@ave,d=14]/100
plot/over/nolabels/thickness=2/color=lightblue hmxl[gt=ntime2@asn,x=180w:80w@ave,y=20s:20n@ave,d=15]/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.2,3.1,0,0,.2 "@CRd)"
!label/nouser 2.2,-0.8,0,0,.2 "@CR HMXL EEP"
label/nouser -1,1.6,0,90,.2 "@CRmixed layer depth (m)"
label/nouser 2.2,-0.8,0,0,.2 "@CR Year of simulation"
!
set view midlow3
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/nolabels/thickness=2/color=red/hlimits=-50:1000:100/vlimits=76:82 hmxl[gt=ntime1@asn,x=20e:100e@ave,y=30n:40s@ave,d=11]/100
plot/over/nolabels/thickness=2/color=(100,55,0) hmxl[gt=ntime1@asn,x=20e:100e@ave,y=30n:40s@ave,d=12]/100
plot/over/nolabels/thickness=2/color=blue hmxl[gt=ntime2@asn,x=20e:100e@ave,y=30n:40s@ave,d=13]/100
plot/over/nolabels/thickness=2/color=black hmxl[gt=ntime1@asn,x=20e:100e@ave,y=30n:40s@ave,d=14]/100
plot/over/nolabels/thickness=2/color=lightblue hmxl[gt=ntime2@asn,x=20e:100e@ave,y=30n:40s@ave,d=15]/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.2,3.1,0,0,.2 "@CRe)"
label/nouser 2.2,-0.8,0,0,.2 "@CR Year of simulation"
!
set view rightlow3
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/nolabels/thickness=2/color=red/hlimits=-50:1000:100/vlimits=90:330/y=40s:80s hmxl[gt=ntime1@asn,i=@ave,j=@ave,d=11]/100
plot/over/nolabels/thickness=2/color=(100,55,0)/y=40s:80s hmxl[gt=ntime1@asn,i=@ave,j=@ave,d=12]/100
plot/over/nolabels/thickness=2/color=blue/y=40s:80s hmxl[gt=ntime2@asn,i=@ave,j=@ave,d=13]/100
plot/over/nolabels/thickness=2/color=black/y=40s:80s hmxl[gt=ntime1@asn,i=@ave,j=@ave,d=14]/100
plot/over/nolabels/thickness=2/color=lightblue/y=40s:80s hmxl[gt=ntime2@asn,i=@ave,j=@ave,d=15]/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 2.3,-0.8,0,0,.2 "@CR HMXL SO"
label/nouser 4.2,3.1,0,0,.2 "@CRf)"
label/nouser 2.2,-0.8,0,0,.2 "@CR Year of simulation"
frame/file="figures/COMFORT/paper/COMFORT_NPP-HMXL_new.pdf"