! Script to plot fluxes and overlay trap data on it

set mem/size=100

! ! ! ! !  CONVERSION FACTORS  ! ! ! ! ! ! ! 
! convert from per second to per day
let mult = 60*60*24

! convert from per day to per year
let mult1 = mult*365

! convert nitrogen to carbon
let ntc = 6.625

! convert mol C to grams C
let mtg = 12.011
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! model data
use tavg.01775.01.01.nc



go prepareFigure Fig22
go "$HOME/opt_postprocessing_scripts/go/ilandscape3x2.jnl"
ppl axlsze,0.15,0.15

set view v1
shade/set/pal=blue_darkred/levels=(-inf)(-0.1,0.1,0.02)/nolab (o_po4[k=1,l=25]-o_po4[k=1,l=20])*1000
ppl shakey 1,0,0.15, , , , , ,`($ppl$ylen)*(0.15)`,`($ppl$ylen)*(0.2)`
ppl shade
go fland;go land black
label/nouser `($ppl$xlen)*.5` `($ppl$ylen)*(1.1)`  0,0,0.2 "PO_4^3^- (mmol m^-^3)"

set view v2
shade/set/pal=blue_darkred/levels=(-14,14,2)/nolab (o_sil[k=1,l=25]-o_sil[k=1,l=20])*1000
ppl shakey 1,0,0.15, , , , , ,`($ppl$ylen)*(0.15)`,`($ppl$ylen)*(0.2)`
ppl shade
go fland;go land black
label/nouser `($ppl$xlen)*.5` `($ppl$ylen)*(1.1)`  0,0,0.2 "Si (mmol m^-^3)"

set view v3
shade/set/pal=blue_darkred/levels=(-inf)(-100,100,20)(inf)/nolab (o_dfe[k=1,l=25]-o_dfe[k=1,l=20])*1e9
ppl shakey 1,0,0.15, , , , , ,`($ppl$ylen)*(0.15)`,`($ppl$ylen)*(0.2)`
ppl shade
go fland;go land black
label/nouser `($ppl$xlen)*.5` `($ppl$ylen)*(1.1)`  0,0,0.2 "Fe (nmol m^-^3)"

set view v4
shade/nokey/pal=blue_darkred/levels=(-inf)(-0.1,0.1,0.02)/nolab (o_po4[z=300,l=25]-o_po4[z=300,l=20])*1000
go fland;go land black
label/nouser `($ppl$xlen)*.5` `($ppl$ylen)*(1.1)`  0,0,0.2 "PO_4^3^- (mmol m^-^3)"

set view v5
shade/nokey/pal=blue_darkred/levels=(-14,14,2)/nolab (o_sil[z=300,l=25]-o_sil[z=300,l=20])*1000
go fland;go land black
label/nouser `($ppl$xlen)*.5` `($ppl$ylen)*(1.1)`  0,0,0.2 "Si (mmol m^-^3)"

set view v6
shade/nokey/pal=blue_darkred/levels=(-inf)(-100,100,20)(inf)/nolab (o_dfe[z=300,l=25]-o_dfe[z=300,l=20])*1e9
go fland;go land black
label/nouser `($ppl$xlen)*.5` `($ppl$ylen)*(1.1)`  0,0,0.2 "Fe (nmol m^-^3)"

go finalizeFigure_psthicken