! global maps of net primary production (NPP) and limitation factors of PO4 and Fe
! NPP maps showing difference between all forcing scenarios (2xCO2, 4xCO2, 2xCO2hosing, 4xCO2hosing) and CONTROL scenario after 500 years
! limitation factor maps (L_PO4, L_Fe) showing difference between strong forcing scenarios (4xCO2, 4xCO2hosing) and CONTROL scenario after 500 years

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

!----------------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 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])
!
!----------------REGION----------------------------------------
!
SET REGION/X=22E:21E
!
!----------------GRAFISCHE AUSGABE---------------------------------

SET TEXT/FONT=symbol
set window/xpix=1920/ypix=1080 1
!
PALETTE mypal_c.spk
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
go mp_hammer
go mp_grid glb_xpp[d=2]
go mp_aspect
shade/noaxis/nolabel/levels=(-inf)(-200,200,10)(inf) (glb_xpp[l=546:555@ave,d=2]-glb_xpp[l=546:555@ave,d=8])*mp_mask, x_page, y_page
go mp_fland
go mp_graticule
label/nouser 9.9,-0.25,0,0,.15 "<&Delta>NPP (g/m<sup>2</sup>)"
label/nouser 4.8,-0.8,0,0,.3 "4xCO<sub>2</sub>"
label/nouser 0.25,4.4,0,0,.3 "a)"
frame/file="figures/COMFORT/paper/COMFORT_primprod_diff500_curve_CO2X4.pdf"

set window/xpix=1920/ypix=1080 2
!
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
go mp_hammer
go mp_grid glb_xpp[d=2]
go mp_aspect
shade/noaxis/nolabel/levels=(-inf)(-200,200,10)(inf) (glb_xpp[l=546:555@ave,d=4]-glb_xpp[l=546:555@ave,d=8])*mp_mask, x_page, y_page
go mp_fland
go mp_graticule
label/nouser 9.9,-0.25,0,0,.15 "<&Delta>NPP (g/m<sup>2</sup>)"
label/nouser 4.8,-0.8,0,0,.3 "2xCO<sub>2</sub>"
label/nouser 0.25,4.4,0,0,.3 "b)"
frame/file="figures/COMFORT/paper/COMFORT_primprod_diff500_curve_CO2X2.pdf"

set window/xpix=1920/ypix=1080 3
!
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
go mp_hammer
go mp_grid glb_xpp[d=2]
go mp_aspect
shade/noaxis/nolabel/levels=(-inf)(-200,200,10)(inf) (glb_xpp[l=546:555@ave,d=6]-glb_xpp[l=546:555@ave,d=8])*mp_mask, x_page, y_page
go mp_fland
go mp_graticule
label/nouser 9.9,-0.25,0,0,.15 "<&Delta>NPP (g/m<sup>2</sup>)"
label/nouser 4.8,-0.8,0,0,.3 "4xCO<sub>2</sub>HOSING"
label/nouser 0.25,4.4,0,0,.3 "c)"
frame/file="figures/COMFORT/paper/COMFORT_primprod_diff500_curve_CO2X4HOSING.pdf"

set window/xpix=1920/ypix=1080 4
!
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
go mp_hammer
go mp_grid glb_xpp[d=2]
go mp_aspect
shade/noaxis/nolabel/levels=(-inf)(-200,200,10)(inf) (glb_xpp[l=546:555@ave,d=10]-glb_xpp[l=546:555@ave,d=8])*mp_mask, x_page, y_page
go mp_fland
go mp_graticule
label/nouser 9.9,-0.25,0,0,.15 "<&Delta>NPP (g/m<sup>2</sup>)"
label/nouser 4.8,-0.8,0,0,.3 "2xCO<sub>2</sub>HOSING"
label/nouser 0.25,4.4,0,0,.3 "d)"
!label/nouser 5,-4.8,0,270,.2 "@CR Fig. 1: longterm NPP change (NPP(t=500)-NPP(t=0)) in g/m^2"
frame/file="figures/COMFORT/paper/COMFORT_primprod_diff500_curve_CO2X2HOSING.pdf"

set window/xpix=1920/ypix=1080 5
!
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
go mp_hammer
go mp_grid plimit1
go mp_aspect
shade/noaxis/nolabel/levels=(-inf)(-8,8,0.4)(inf) (plimit1[l=546:555@ave,z=100]-plimit7[l=546:555@ave,z=100])*100*mp_mask, x_page, y_page
go mp_fland
go mp_graticule
label/nouser 9.8,-0.25,0,0,.15 "L<sub>PO4</sub>(x100)"
label/nouser 4.8,-0.8,0,0,.3 "4xCO<sub>2</sub>"
label/nouser 0.25,4.4,0,0,.3 "b)"
frame/file="figures/COMFORT/paper/COMFORT_LPO4_diff500_curve_CO2X4.pdf"

set window/xpix=1920/ypix=1080 6
!
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
go mp_hammer
go mp_grid plimit5
go mp_aspect
shade/noaxis/nolabel/levels=(-inf)(-8,8,0.4)(inf) (plimit5[l=546:555@ave,z=100]-plimit7[l=546:555@ave,z=100])*100*mp_mask, x_page, y_page
go mp_fland
go mp_graticule
label/nouser 9.8,-0.25,0,0,.15 "L<sub>PO4</sub>(x100)"
label/nouser 4.8,-0.8,0,0,.3 "4xCO<sub>2</sub>HOSING"
label/nouser 0.25,4.4,0,0,.3 "d)"
frame/file="figures/COMFORT/paper/COMFORT_LPO4_diff500_curve_CO2X4Hosing.pdf"

set window/xpix=1920/ypix=1080 7
!
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
go mp_hammer
go mp_grid fvalue1
go mp_aspect
shade/noaxis/nolabel/levels=(-inf)(-8,8,0.4)(inf) (fvalue1[l=546:555@ave,z=100]-fvalue7[l=546:555@ave,z=100])*100*mp_mask, x_page, y_page
go mp_fland
go mp_graticule
label/nouser 9.8,-0.25,0,0,.15 "L<sub>Fe</sub>(x100)"
label/nouser 4.8,-0.8,0,0,.3 "4xCO<sub>2</sub>"
label/nouser 0.25,4.4,0,0,.3 "a)"
frame/file="figures/COMFORT/paper/COMFORT_LFe_diff500_curve_CO2X4.pdf"

set window/xpix=1920/ypix=1080 8
!
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
go mp_hammer
go mp_grid fvalue5
go mp_aspect
shade/noaxis/nolabel/levels=(-inf)(-8,8,0.4)(inf) (fvalue5[l=546:555@ave,z=100]-fvalue7[l=546:555@ave,z=100])*100*mp_mask, x_page, y_page
go mp_fland
go mp_graticule
label/nouser 9.8,-0.25,0,0,.15 "L<sub>Fe</sub>(x100)"
label/nouser 4.8,-0.8,0,0,.3 "4xCO<sub>2</sub>HOSING"
label/nouser 0.25,4.4,0,0,.3 "c)"
frame/file="figures/COMFORT/paper/COMFORT_LFe_diff500_curve_CO2X4Hosing.pdf"