!NPZD inventory for 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 basin mask file.

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





!!!!! 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/PHYTmm.nc" !2
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/DIAZmm.nc"		!3
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/ZOOPmm.nc"	!4
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/DETRmm.nc"	!5
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/NO3mm.nc"	!6



! FeMAsk

use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/PHYTmm.nc" !7
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/DIAZmm.nc"		!8
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/ZOOPmm.nc"	!9
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/DETRmm.nc"	!10
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/NO3mm.nc"	!11



! FeDyn1

use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/PHYTmm.nc" !12
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/DIAZmm.nc"		!13
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/ZOOPmm.nc"		!14
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/DETRmm.nc"	!15
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/NO3mm.nc"	!16







let dummy = compressl(O_no3[d=1,l=1,k=1]*0.)
def sym oper = gz=zax,x=@din,y=@din,l=@ave



!list data out

list/nohead/file=4surfaceNPZD.txt phyt[d=2,($oper),k=1:2@din]+diaz[d=3,($oper),k=1:2@din],zoop[d=4,($oper),k=1:2@din],detr[d=5,($oper),k=1:2@din],no3[d=6,($oper),k=1:2@din]
list/nohead/file=4surfaceNPZD.txt/append phyt[d=7,($oper),k=1:2@din]+diaz[d=8,($oper),k=1:2@din],zoop[d=9,($oper),k=1:2@din],detr[d=10,($oper),k=1:2@din],no3[d=11,($oper),k=1:2@din]
list/nohead/file=4surfaceNPZD.txt/append phyt[d=12,($oper),k=1:2@din]+diaz[d=13,($oper),k=1:2@din],zoop[d=14,($oper),k=1:2@din],detr[d=15,($oper),k=1:2@din],no3[d=16,($oper),k=1:2@din]