! shows marine tracer dispersion


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/INK/CO2X4_HOSING/snap_ocean_ink.nc" 1
use "/p/projects/climber3/hofmann/RL_NEW/INK/CONTROL/snap_ocean_ink.nc" 2
use "/p/projects/climber3/hofmann/RL_NEW/CO2X4_HOSING/history.nc" 3
use "/p/projects/climber3/hofmann/RL_NEW/CONTROL/history.nc" 4
use "/p/projects/climber3/hofmann/RL_NEW/INK/CONTROL/snap_ink_indian.nc" 5
use "/p/projects/climber3/hofmann/RL_NEW/INK/CO2X4/snap_ocean_ink.nc" 6
!
!----------------GRAFIKEINSTELLUNGEN-------------------------------
! for 12 plots on a page (4 columns)
define viewport/xlimits=0,0.35/ylim=0.67,0.99 pos11
define viewport/xlimits=0.21,0.56/ylim=0.67,0.99 pos12
define viewport/xlimits=0.42,0.77/ylim=0.67,0.99 pos13
define viewport/xlimits=0.63,0.98/ylim=0.67,0.99 pos14
define viewport/xlimits=0,0.35/ylim=0.36,0.68 pos21
define viewport/xlimits=0.21,0.56/ylim=0.36,0.68 pos22
define viewport/xlimits=0.42,0.77/ylim=0.36,0.68 pos23
define viewport/xlimits=0.63,0.98/ylim=0.36,0.68 pos24
define viewport/xlimits=0,0.35/ylim=0.05,0.37 pos31
define viewport/xlimits=0.21,0.56/ylim=0.05,0.37 pos32
define viewport/xlimits=0.42,0.77/ylim=0.05,0.37 pos33
define viewport/xlimits=0.63,0.98/ylim=0.05,0.37 pos34

! for 8 plots on a page (4 columns)
define viewport/xlimits=0,0.35/ylim=0.54,0.99 up1
define viewport/xlimits=0.21,0.56/ylim=0.54,0.99 up2
define viewport/xlimits=0.42,0.77/ylim=0.54,0.99 up3
define viewport/xlimits=0.63,0.98/ylim=0.54,0.99 up4
define viewport/xlimits=0,0.35/ylim=0.08,0.53 low1
define viewport/xlimits=0.21,0.56/ylim=0.08,0.53 low2
define viewport/xlimits=0.42,0.77/ylim=0.08,0.53 low3
define viewport/xlimits=0.63,0.98/ylim=0.08,0.53 low4
!
!----------------GRAFISCHE AUSGABE---------------------------------

SET TEXT/FONT=symbol

! vertical cut at 30W through Atlantic ocean, comparing Atlantic overturning (streamfunction) with Tracer concentrations
set window/xpix=1920/ypix=1280 1
!
set view pos11
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
PALETTE blue_darkorange
shade/nolabel/levels=(-inf)(0.0,400.0,10)(inf)/nokey tracer_16[d=2,l=1:10@ave,x=30W,y=30S:80N]
contour/overlay/nolabel/noaxis/levels=(-inf)(-10,20.0,2)(inf) atl_over[d=4,l=1:10@ave,x=30W,y=30S:80N]
label/nouser 1.7,-0.8,0,0,.15 "t=5"
label/nouser -0.9,1.5,0,90,.2 "@CR CONTROL"
!
set view pos12
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,0,0
shade/nolabel/levels=(-inf)(0.0,400.0,10)(inf)/nokey tracer_16[d=2,l=136:145@ave,x=30W,y=30S:80N]*1000
contour/overlay/nolabel/noaxis/levels=(-inf)(-10,20.0,2)(inf) atl_over[d=4,l=136:145@ave,x=30W,y=30S:80N]
label/nouser 1.7,-0.8,0,0,.15 "t=140"
!
set view pos13
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,0,0
shade/nolabel/levels=(-inf)(0.0,400.0,10)(inf)/nokey tracer_16[d=2,l=276:285@ave,x=30W,y=30S:80N]*1000
contour/overlay/nolabel/noaxis/levels=(-inf)(-10,20.0,2)(inf) atl_over[d=4,l=276:285@ave,x=30W,y=30S:80N]
label/nouser 1.7,-0.8,0,0,.15 "t=280"
!
set view pos14
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,0,0
shade/nolabel/levels=(-inf)(0.0,400.0,10)(inf) tracer_16[d=2,l=416:425@ave,x=30W,y=30S:80N]*1000
contour/overlay/nolabel/noaxis/levels=(-inf)(-10,20.0,2)(inf) atl_over[d=4,l=416:425@ave,x=30W,y=30S:80N]
label/nouser 1.7,-0.8,0,0,.15 "t=420"
label/nouser 4.6,-0.18,0,0,.12 "@CR Fe (nmol/l)"
!
set view pos21
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
shade/nolabel/levels=(-inf)(0.0,400.0,10)(inf)/nokey tracer_16[d=6,l=1:10@ave,x=30W,y=30S:80N]*1000
contour/overlay/nolabel/noaxis/levels=(-inf)(-10,20.0,2)(inf) atl_over[d=3,l=1:10@ave,x=30W,y=30S:80N]
label/nouser 1.7,-0.8,0,0,.15 "t=5"
label/nouser -0.9,1.5,0,90,.2 "@CR 4xCO_2"
!
set view pos22
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,0,0
shade/nolabel/levels=(-inf)(0.0,400.0,10)(inf)/nokey tracer_16[d=6,l=136:145@ave,x=30W,y=30S:80N]*1000
contour/overlay/nolabel/noaxis/levels=(-inf)(-10,20.0,2)(inf) atl_over[d=3,l=136:145@ave,x=30W,y=30S:80N]
label/nouser 1.7,-0.8,0,0,.15 "t=140"
!
set view pos23
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,0,0
shade/nolabel/levels=(-inf)(0.0,400.0,10)(inf)/nokey tracer_16[d=6,l=276:285@ave,x=30W,y=30S:80N]*1000
contour/overlay/nolabel/noaxis/levels=(-inf)(-10,20.0,2)(inf) atl_over[d=3,l=276:285@ave,x=30W,y=30S:80N]
label/nouser 1.7,-0.8,0,0,.15 "t=280"
!
set view pos24
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,0,0
shade/nolabel/levels=(-inf)(0.0,400.0,10)(inf) tracer_16[d=6,l=416:425@ave,x=30W,y=30S:80N]*1000
contour/overlay/nolabel/noaxis/levels=(-inf)(-10,20.0,2)(inf) atl_over[d=3,l=416:425@ave,x=30W,y=30S:80N]
label/nouser 1.7,-0.8,0,0,.15 "t=420"
label/nouser 4.6,-0.18,0,0,.12 "@CR Fe (nmol/l)"
!
set view pos31
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
shade/nolabel/levels=(-inf)(0.0,400.0,10)(inf)/nokey tracer_16[d=1,l=1:10@ave,x=30W,y=30S:80N]*1000
contour/overlay/nolabel/noaxis/levels=(-inf)(-10,20.0,2)(inf) atl_over[d=3,l=1:10@ave,x=30W,y=30S:80N]
label/nouser 1.7,-0.8,0,0,.15 "t=5"
label/nouser -0.9,1.5,0,90,.2 "@CR 4xCO_2HOSING"
!
set view pos32
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,0,0
shade/nolabel/levels=(-inf)(0.0,400.0,10)(inf)/nokey tracer_16[d=1,l=136:145@ave,x=30W,y=30S:80N]*1000
contour/overlay/nolabel/noaxis/levels=(-inf)(-10,20.0,2)(inf) atl_over[d=3,l=136:145@ave,x=30W,y=30S:80N]
label/nouser 1.7,-0.8,0,0,.15 "t=140"
!
set view pos33
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,0,0
shade/nolabel/levels=(-inf)(0.0,400.0,10)(inf)/nokey tracer_16[d=1,l=276:285@ave,x=30W,y=30S:80N]*1000
contour/overlay/nolabel/noaxis/levels=(-inf)(-10,20.0,2)(inf) atl_over[d=3,l=276:285@ave,x=30W,y=30S:80N]
label/nouser 1.7,-0.8,0,0,.15 "t=280"
!
set view pos34
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,0,0
shade/nolabel/levels=(-inf)(0.0,400.0,10)(inf) tracer_16[d=1,l=416:425@ave,x=30W,y=30S:80N]*1000
contour/overlay/nolabel/noaxis/levels=(-inf)(-10,20.0,2)(inf) atl_over[d=3,l=416:425@ave,x=30W,y=30S:80N]
label/nouser 1.7,-0.8,0,0,.15 "t=420"
label/nouser 4.6,-0.18,0,0,.12 "@CR Fe (nmol/l)"
!frame/file="figures/COMFORT/paper/COMFORT_ink_30W.pdf"


! global maps of tracer dispersion injected in Arabian Sea (Indian Ocean)
set window/xpix=1920/ypix=1080 2
!
set view up1
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
PALETTE cmocean_dense
shade/nolabel/levels=(-inf)(0.0,100.0,5.0)(inf)/nokey tracer_16[d=5,l=1:10@ave,z=100]*1000 ! /levels=(-inf)(-2.0,2.0,0.1)(inf)
go fland
label/nouser 1.7,-0.8,0,0,.15 "Ink (t=5)"
label/nouser -0.9,1.5,0,90,.2 "@CR CONTROL"
!
set view up2
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
!PALETTE no_green_centered
shade/nolabel/levels=(-inf)(0.0,100.0,5.0)(inf)/nokey tracer_16[d=5,l=36:45@ave,z=100]*1000 ! /levels=(-inf)(-2.0,2.0,0.1)(inf)
go fland
label/nouser 1.7,-0.8,0,0,.15 "Ink (t=40)"
!
set view up3
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
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,100.0,5.0)(inf)/nokey tracer_16[d=5,l=76:85@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "Ink (t=80)"
!
set view up4
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
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,100.0,5.0)(inf) tracer_16[d=5,l=116:125@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "Ink (t=120)"
!
set view low1
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
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,100.0,5.0)(inf)/nokey tracer_16[d=5,l=156:165@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "Ink (t=160)"
!
set view low2
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
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,100.0,5.0)(inf)/nokey tracer_16[d=5,l=196:205@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "Ink (t=200)"
!
set view low3
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
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,100.0,5.0)(inf)/nokey tracer_16[d=5,l=236:245@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "Ink (t=240)"
!
set view low4
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
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,100.0,5.0)(inf) tracer_16[d=5,l=271:280@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "Ink (t=275)"
!frame/file="figures/COMFORT/paper/COMFORT_ink_indian.pdf"


! global maps of tracer dispersion injected near Westafrican coast (North Atlantic), comparing three scenarios: CONTROL, 4xCO2 and 4xCO2 with NA freshwater forcing
set window/xpix=1920/ypix=1280 3
!
set view pos11
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
PALETTE cmocean_dense ! blue_darkorange cmocean_solar no_green_centered orange_purple
shade/nolabel/levels=(-inf)(0.0,20.0,1)(inf)/nokey tracer_16[d=2,l=1:10@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "t=5"
label/nouser -0.9,1.5,0,90,.2 "@CR CONTROL"
!
set view pos12
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,0,0
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,20.0,1)(inf)/nokey tracer_16[d=2,l=136:145@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "t=140"
!
set view pos13
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,0,0
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,20.0,1)(inf)/nokey tracer_16[d=2,l=276:285@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "t=280"
!
set view pos14
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,0,0
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,20.0,1)(inf) tracer_16[d=2,l=416:425@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "t=420"
!label/nouser 4.6,-0.18,0,0,.12 "@CR Fe (nmol/l)"
!
set view pos21
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
!PALETTE no_green_centered
shade/nolabel/levels=(-inf)(0.0,20.0,1)(inf)/nokey tracer_16[d=6,l=1:10@ave,z=100]*1000 
go fland
label/nouser 1.7,-0.8,0,0,.15 "t=5"
label/nouser -0.9,1.5,0,90,.2 "@CR 4xCO_2"
!
set view pos22
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,0,0
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,20.0,1)(inf)/nokey tracer_16[d=6,l=136:145@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "t=140"
!
set view pos23
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,0,0
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,20.0,1)(inf)/nokey tracer_16[d=6,l=276:285@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "t=280"
!
set view pos24
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,0,0
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,20.0,1)(inf) tracer_16[d=6,l=416:425@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "t=420"
!
set view pos31
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
!PALETTE no_green_centered
shade/nolabel/levels=(-inf)(0.0,20.0,1)(inf)/nokey tracer_16[d=1,l=1:10@ave,z=100]*1000 
go fland
label/nouser 1.7,-0.8,0,0,.15 "t=5"
label/nouser -0.9,1.5,0,90,.2 "@CR 4xCO_2HOSING"
!
set view pos32
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,0,0
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,20.0,1)(inf)/nokey tracer_16[d=1,l=136:145@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "t=140"
!
set view pos33
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,0,0
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,20.0,1)(inf)/nokey tracer_16[d=1,l=276:285@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "t=280"
!
set view pos34
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,0,0
!PALETTE centered
shade/nolabel/levels=(-inf)(0.0,20.0,1)(inf) tracer_16[d=1,l=416:425@ave,z=100]*1000
go fland
label/nouser 1.7,-0.8,0,0,.15 "t=420"
!frame/file="figures/COMFORT/paper/COMFORT_ink_100m.pdf"