!surface fluxes for 3 calibrations.
! fluxes include : NPP, Expo, remain, scret, recy.

!!! 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/O_phytnpp.nc" !2
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/O_diaznpp.nc"		!3
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/O_detrexp.nc"	!4
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/O_detrremi.nc"	!5
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/O_phytrecy.nc"	!6
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/O_excret.nc"	!7
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20181006_i90_2_norion_046639_extradiag/NO3mm.nc"	!8


! FeMAsk
use "/Volumes/BM_2017_36/UVicWorking/2.9_ln_iron/data.100.100.19/O_fe_dissolved.nc" !9
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/NO3mm.nc" !10
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/O_phytnpp.nc"		!11
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/O_diaznpp.nc"	!12
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/O_detrexp.nc"	!13
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/O_detrremi.nc"	!14
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/O_phytrecy.nc"	!15
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180828_i116_5_0462_FeMask_extradiag/O_excret.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/O_phytnpp.nc"		!18
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/O_diaznpp.nc"		!19
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/O_detrexp.nc"	!20
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/O_detrremi.nc"	!21
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/O_phytrecy.nc"	!22
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/O_excret.nc"	!23
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/O_avej_D.nc"	!24
use "/Volumes/BM_2017_36/HLRN_UVOKWorking/UVOKWorking/best_pops/pop_hlrn_20180829_i112_5_04427_FeDyn1_extradiag/O_deffe.nc" !25






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=4surfaceflux.txt o_phytnpp[d=2,($oper),k=1:2@din]+o_diaznpp[d=3,($oper),k=1:2@din],o_detrexp[d=4,($oper),k=2]*80,o_detrremi[d=5,($oper),k=1:2@din],o_phytrecy[d=6,($oper),k=1:2@din],o_excret[d=7,($oper),k=1:2@din]
list/nohead/file=4surfaceflux.txt/append o_phytnpp[d=11,($oper),k=1:2@din]+o_diaznpp[d=12,($oper),k=1:2@din],o_detrexp[d=13,($oper),k=2]*80,o_detrremi[d=14,($oper),k=1:2@din],o_phytrecy[d=15,($oper),k=1:2@din],o_excret[d=16,($oper),k=1:2@din]
list/nohead/file=4surfaceflux.txt/append o_phytnpp[d=18,($oper),k=1:2@din]+o_diaznpp[d=19,($oper),k=1:2@din],o_detrexp[d=20,($oper),k=2]*80,o_detrremi[d=21,($oper),k=1:2@din],o_phytrecy[d=22,($oper),k=1:2@din],o_excret[d=23,($oper),k=1:2@din]