! global maps of PO4 and 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_HOSING/snap_ocean_hosing_1000.nc" 2
use "/p/projects/climber3/hofmann/RL_NEW/CONTROL/snap_ocean_1000.nc" 3

!----------------REGION----------------------------------------
!
SET REGION/X=22E:21E
!
!----------------GRAFISCHE AUSGABE---------------------------------


SET TEXT/FONT=symbol
set window/xpix=1920/ypix=1280 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 ! eckert_iii mp_wagner_vii
go mp_grid tracer_13[d=1]
go mp_aspect
shade/noaxis/nolabel/levels=(-inf)(-0.8,0.8,0.04)(inf) (tracer_13[l=546:555@ave,z=100,d=1]-tracer_13[l=546:555@ave,z=100,d=3])*1000*mp_mask, x_page, y_page ! 
go mp_fland
go mp_graticule
label/nouser 9.9,-0.25,0,0,.15 "Fe (nmol/l)"
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_Fe_diff500_curve_CO2X4.pdf"

set window/xpix=1920/ypix=1200 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 tracer_13[d=1]
go mp_aspect
shade/noaxis/nolabel/levels=(-inf)(-0.8,0.8,0.04)(inf) (tracer_13[l=546:555@ave,z=100,d=2]-tracer_13[l=546:555@ave,z=100,d=3])*1000*mp_mask, x_page, y_page
go mp_fland
go mp_graticule
label/nouser 9.9,-0.25,0,0,.15 "Fe (nmol/l)"
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_Fe_diff500_curve_CO2X4Hosing.pdf"

set window/xpix=1920/ypix=1440 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 tracer_13[d=1]
go mp_aspect
shade/noaxis/nolabel/levels=(-inf)(-0.8,0.8,0.04)(inf) (tracer_06[l=546:555@ave,z=100,d=1]-tracer_06[l=546:555@ave,z=100,d=3])*mp_mask, x_page, y_page
go mp_fland
go mp_graticule
label/nouser 9.5,-0.25,0,0,.15 "PO<sub>4</sub> (mmol/l)"
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_PO4_diff500_curve_CO2X4.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 tracer_13[d=1]
go mp_aspect
shade/noaxis/nolabel/levels=(-inf)(-0.8,0.8,0.04)(inf) (tracer_06[l=546:555@ave,z=100,d=2]-tracer_06[l=546:555@ave,z=100,d=3])*mp_mask, x_page, y_page
go mp_fland
go mp_graticule
label/nouser 9.5,-0.25,0,0,.15 "PO<sub>4</sub> (mmol/l)"
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_PO4_diff500_curve_CO2X4Hosing.pdf"