! plots scavenging rates under 4xCO2: compares normal 4xCO2-scavenging and scavenging without 4xCO2 impact (scavenging rate frozen at baseline)

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

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


!----------------SIMULATIONSDATEN LADEN----------------------------
use "/p/projects/climber3/hofmann/RL_NEW/CO2X4/snap_ocean.nc" 1 
use "/p/projects/climber3/hofmann/RL_NEW/CO2X4/history.nc" 2
use "/p/projects/climber3/hofmann/RL_NEW/CO2X4_FROZ_SCAV/snap_ocean_froz_scav.nc" 3 
use "/p/projects/climber3/hofmann/RL_NEW/CO2X4_FROZ_SCAV/history.nc" 4
!----------------GRAFIKEINSTELLUNGEN-------------------------------
!
! for 2 plots on a page
define viewport/xlimits=0,0.9/ylim=0.05,0.55 low
define viewport/xlimits=0,0.9/ylim=0.5,1.0 up
!
!----------------ZEITACHSE----------------------------------------
!
define axis/units=years/t=-50:500:1 nt1
define grid/T=nt1 ntime1
!
!----------------FORMELN---------------------------------
!
let arg1 = CACO312[d=1]+tracer_07[d=1]
let scav1 = 0.01*(arg1^0.58)
!
let arg3 = CACO312[d=3]+tracer_07[d=3]
let scav3 = 0.01*(arg3^0.58)
!
!----------------GRAFISCHE AUSGABE---------------------------------

SET TEXT/FONT=symbol

set window/xpix=1680/ypix=640 1
!
set view up
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:450:25 scav1[gt=ntime1@asn,k=@ave,x=@ave,y=@ave]
plot/over/nolabels/thickness=2/color=green/hlimits=-50:450:25 scav3[gt=ntime1@asn,k=@ave,x=@ave,y=@ave]
!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 -1,1,0,90,.2 "@CR Scavenging"
label/nouser 6.5,-0.8,0,0,.2 "year of simulation"
!
set view low
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:450:25 glb_xpp[gt=ntime1@asn,i=@din,j=@din,d=2]*1.e-15
plot/over/nolabels/thickness=2/color=green/hlimits=-50:450:25 glb_xpp[gt=ntime1@asn,i=@din,j=@din,d=4]*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 -1,1,0,90,.2 "@CR NPP"
label/nouser 6.5,-0.8,0,0,.2 "year of simulation"
!frame/file="figures/COMFORT/paper/COMFORT_plots_froz_scav.pdf"