!surface difference in NPP, Export and N concentration between 3 calibrations.

!!! the output files form TMM may have no depth_edge, which leads to errors in @DIN and @AVE.so here we define the right boundary for the depth boxes.

let coordin={17.5, 82.5, 177.5, 302.5, 457.5, 642.5, 857.5, 1102.5, 1377.5,1682.5, 2017.5, 2382.5, 2777.5, 3202.5, 3657.5, 4142.5, 4657.5, 5202.5,5777.5}
let bond={ 0, 50, 130, 240, 380, 550, 750, 980, 1240, 1530, 1850, 2200,2580, 2990, 3430, 3900, 4400, 4930, 5490, 6080 }
define axis/z/depth/bounds zax=coordin,bond
show axis/all

!!! load in the grid mask file.

use "/Volumes/BM_2017_36/UVicWorking/2.9_ln_iron/quickerTMMphysicOnline/tavg.08100.01.01.nc" !1

!!! load in observations
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/update/Forcing_newiniforopt_20171219/obs_for_calibration/O_no3.nc"	!2
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/update/Forcing_newiniforopt_20171219/obs_for_calibration/O_po4.nc"	!3
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/update/Forcing_newiniforopt_20171219/obs_for_calibration/O_o2.nc"	!4



!!!!! load in data from 4 calibrations

! noFe

use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/NO3mm.nc"	!5
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/PO4mm.nc"	!6
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/O2mm.nc"	!7
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/O_phytnpp.nc"	!8
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/O_diaznpp.nc"	!9
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/O_detrexp.nc"	!10

! FeMAsk
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/NO3mm.nc"	!11
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/PO4mm.nc"	!12
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/O2mm.nc"	!13
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/O_phytnpp.nc"	!14
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/O_diaznpp.nc"	!15
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/O_detrexp.nc"	!16


! FeDyn1
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/NO3mm.nc"	!17
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/PO4mm.nc"	!18
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/O2mm.nc"	!19
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/O_phytnpp.nc"	!20
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/O_diaznpp.nc"	!21
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/O_detrexp.nc"	!22
let possibleTragname={"NO3","PO4","O2","DFE"}
let possibleObsVariname={"O_NO3","O_PO4","O_O2","dfe_obs_lt10"}
let possibleXlims={"-8:8:2","-0.55:0.55:0.15","-0.055:0.055:0.02","0.:0.0011:0.00025"}
let possibleScale={"*1000","*1000","*1","*1000"}
let possibleUnite={"mMol/M^3","mMol/M^3","Mol/M^3","mMol/M^3"}


let griddata=o_dfe[d=1]
def sym opt1=k=1:2@din,l=@ave,gxyz=griddata
def sym opt2=k=2@din,l=@ave,gxyz=griddata
def sym opt3=k=1:2@ave,l=@ave,gxyz=griddata
def sym s2y=3600*24*365*6.625*12


go ilandscape3x2.jnl
go prepareFigure surfaceNPPEXPONdiff
set win 1
go transparency

set vie v1
sha/nolab/pal=blue_orange/lev="(-inf)(-400,400,50)(inf)"/set ((o_phytnpp[d=20,($opt1)]+o_diaznpp[d=21,($opt1)])-(o_phytnpp[d=8,($opt1)]+o_diaznpp[d=9,($opt1)]))*($s2y)
ppl ylab "Latitude"
ppl shade
go land
label `360/2`, `-125`,0,0,.15 @SGd@SRNPP FeDyn-NoFe
label `360/2`, `-135`,0,0,.15 (g C m^-^2yr^-^1)
label `360/2`, `100`,0,0,.15 Longitude
set vie v2
sha/nolab/pal=blue_orange/lev="(-inf)(-60,60,10)(inf)" (o_detrexp[d=22,($opt2)]-o_detrexp[d=10,($opt2)])*($s2y)
ppl shade
go land
label `360/2`, `-125`,0,0,.15 @SGd@SRExport FeDyn-NoFe
label `360/2`, `-135`,0,0,.15 (g C m^-^2yr^-^1)
label `360/2`, `100`,0,0,.15 Longitude
set vie v3
sha/nolab/pal=blue_orange/lev="(-inf)(-9,9,1)(inf)" (no3[d=17,($opt3)]-no3[d=5,($opt3)])
ppl shade
go land
label `360/2`, `-125`,0,0,.15 @SGd@SRNO_3 FeDyn-NoFe
label `360/2`, `-135`,0,0,.15 (mmol m^-^3)
label `360/2`, `100`,0,0,.15 Longitude

!go transparency
set vie v4
sha/nolab/pal=blue_orange/lev="(-inf)(-400,400,50)(inf)"/set ((o_phytnpp[d=14,($opt1)]+o_diaznpp[d=15,($opt1)])-(o_phytnpp[d=8,($opt1)]+o_diaznpp[d=9,($opt1)]))*($s2y)
ppl ylab "Latitude"
ppl shade
go land
label `360/2`, `-125`,0,0,.15 @SGd@SRNPP FeMask-NoFe
label `360/2`, `-135`,0,0,.15 (g C m^-^2yr^-^1)
set vie v5
sha/nolab/pal=blue_orange/lev="(-inf)(-60,60,10)(inf)" (o_detrexp[d=16,($opt2)]-o_detrexp[d=10,($opt2)])*($s2y)
ppl shade
go land
label `360/2`, `-125`,0,0,.15 @SGd@SRExport FeMask-NoFe 
label `360/2`, `-135`,0,0,.15 (g C m^-^2yr^-^1)

set vie v6
sha/nolab/pal=blue_orange/lev="(-inf)(-9,9,1)(inf)" (no3[d=11,($opt3)]-no3[d=5,($opt3)])
ppl shade
go land
label `360/2`, `-125`,0,0,.15 @SGd@SRNO_3 FeMask-NoFe 
label `360/2`, `-135`,0,0,.15 (mmol m^-^3)
go finalizefigure_psthicken