!script by bencat.xuan@gmail.com



use /Volumes/BM_2017_36/gauss/noIron_spinup_fct_drift/1800_2100/tavg.01801.01.01.nc
use /Volumes/BM_2017_36/gauss/ironMask_spinup_fct_drift/1800_2100/tavg.01801.01.01.nc
use /Volumes/BM_2017_36/gauss/ironDyn1_spinup_fct_drift/1800_2100/tavg.01801.01.01.nc




set mem/size=2048

def sym op1=x=@din,y=@din,k=@din
def sym op2=x=@din,y=@din,k=1:2@din
def sym op3=x=@din,y=@din,k=2@din
def sym op4=z=350
def sym op6=k=3
def sym op7=k=1


def sym op5=k=@din

let k1n=0.7
let redptn=1/16
let k1p=k1n*redptn



let nofe_no3_1800=o_no3[d=1,l=1]*1000
let nofe_no3_2100=o_no3[d=1,l=300]*1000
let nofe_po4_1800=o_po4[d=1,l=1]*1000
let nofe_po4_2100=o_po4[d=1,l=300]*1000
let nofe_nlim_1800=o_no3[d=1,l=1]*1000/(k1n+o_no3[d=1,l=1]*1000)
let nofe_nlim_2100=o_no3[d=1,l=300]*1000/(k1n+o_no3[d=1,l=300]*1000)
let nofe_plim_1800=o_po4[d=1,l=1]*1000/(k1p+o_po4[d=1,l=1]*1000)
let nofe_plim_2100=o_po4[d=1,l=300]*1000/(k1p+o_po4[d=1,l=300]*1000)
let nofe_no3diff=nofe_no3_2100-nofe_no3_1800
let nofe_po4diff=nofe_po4_2100-nofe_po4_1800
let nofe_nlimdiff=nofe_nlim_2100-nofe_nlim_1800
let nofe_plimdiff=nofe_plim_2100-nofe_plim_1800
let nofe_pstar_1800=nofe_po4_1800-nofe_no3_1800/16
let nofe_pstar_2100=nofe_po4_2100-nofe_no3_2100/16
let nofe_pstar_diff=nofe_pstar_2100-nofe_pstar_1800
let nofe_deni_1800=o_denitr[d=1,l=1,k=@din]
let nofe_deni_2100=o_denitr[d=1,l=300,k=@din]
let nofe_denidiff=o_denitr[d=1,l=300,k=@din]-o_denitr[d=1,l=1,k=@din]
let/bad=-999 nofe_denidiff_n0= if nofe_denidiff gt (1e-15) or nofe_denidiff lt (-1e-15) then nofe_denidiff else (-999)


let nofe_pstar=o_po4[d=1,($op5)]*1000-o_no3[d=1,($op5)]*1000/16
let nofe_pstar_l=nofe_pstar-(o_no3[d=1,($op5),l=1]*1000-o_po4[d=1,($op5),l=1]*1000*16)

let femask_no3_1800=o_no3[d=2,l=1]*1000
let femask_no3_2100=o_no3[d=2,l=300]*1000
let femask_po4_1800=o_po4[d=2,l=1]*1000
let femask_po4_2100=o_po4[d=2,l=300]*1000
let femask_nlim_1800=o_no3[d=2,l=1]*1000/(k1n+o_no3[d=2,l=1]*1000)
let femask_nlim_2100=o_no3[d=2,l=300]*1000/(k1n+o_no3[d=2,l=300]*1000)
let femask_plim_1800=o_po4[d=2,l=1]*1000/(k1p+o_po4[d=2,l=1]*1000)
let femask_plim_2100=o_po4[d=2,l=300]*1000/(k1p+o_po4[d=2,l=300]*1000)
let femask_no3diff=femask_no3_2100-femask_no3_1800
let femask_po4diff=femask_po4_2100-femask_po4_1800
let femask_nlimdiff=femask_nlim_2100-femask_nlim_1800
let femask_plimdiff=femask_plim_2100-femask_plim_1800
let femask_pstar_1800=femask_po4_1800-femask_no3_1800/16
let femask_pstar_2100=femask_po4_2100-femask_no3_2100/16
let femask_pstar_diff=femask_pstar_2100-femask_pstar_1800
let femask_deni_1800=o_denitr[d=2,l=1,k=@din]
let femask_deni_2100=o_denitr[d=2,l=300,k=@din]
let femask_denidiff=o_denitr[d=2,l=300,k=@din]-o_denitr[d=2,l=1,k=@din]
let/bad=-999 femask_denidiff_n0= if femask_denidiff gt (1e-15) or femask_denidiff lt (-1e-15) then femask_denidiff else (-999)

let femask_pstar=o_po4[d=2,($op5)]*1000-o_no3[d=2,($op5)]*1000/16
let femask_pstar_l=femask_pstar-(o_no3[d=2,($op5),l=1]*1000-o_po4[d=2,($op5),l=1]*1000*16)


let fedyn1_no3_1800=o_no3[d=3,l=1]*1000
let fedyn1_no3_2100=o_no3[d=3,l=300]*1000
let fedyn1_po4_1800=o_po4[d=3,l=1]*1000
let fedyn1_po4_2100=o_po4[d=3,l=300]*1000
let fedyn1_nlim_1800=o_no3[d=3,l=1]*1000/(k1n+o_no3[d=3,l=1]*1000)
let fedyn1_nlim_2100=o_no3[d=3,l=300]*1000/(k1n+o_no3[d=3,l=300]*1000)
let fedyn1_plim_1800=o_po4[d=3,l=1]*1000/(k1p+o_po4[d=3,l=1]*1000)
let fedyn1_plim_2100=o_po4[d=3,l=300]*1000/(k1p+o_po4[d=3,l=300]*1000)
let fedyn1_no3diff=fedyn1_no3_2100-fedyn1_no3_1800
let fedyn1_po4diff=fedyn1_po4_2100-fedyn1_po4_1800
let fedyn1_nlimdiff=fedyn1_nlim_2100-fedyn1_nlim_1800
let fedyn1_plimdiff=fedyn1_plim_2100-fedyn1_plim_1800
let fedyn1_pstar_1800=fedyn1_po4_1800-fedyn1_no3_1800/16
let fedyn1_pstar_2100=fedyn1_po4_2100-fedyn1_no3_2100/16
let fedyn1_pstar_diff=fedyn1_pstar_2100-fedyn1_pstar_1800
let fedyn1_deni_1800=o_denitr[d=3,l=1,k=@din]
let fedyn1_deni_2100=o_denitr[d=3,l=300,k=@din]
let fedyn1_denidiff=o_denitr[d=3,l=300,k=@din]-o_denitr[d=3,l=1,k=@din]
let/bad=-999 fedyn1_denidiff_n0= if fedyn1_denidiff gt (1e-15) or fedyn1_denidiff lt (-1e-15) then fedyn1_denidiff else (-999)

let fedyn1_pstar=o_po4[d=3,($op5)]*1000-o_no3[d=3,($op5)]*1000/16
let fedyn1_pstar_l=fedyn1_pstar-(o_no3[d=3,($op5),l=1]*1000-o_po4[d=3,($op5),l=1]*1000*16)




go isquare3x3_ext.jnl

use "/Users/wyao/Desktop/UVicWorking/2.9_ln_iron/dataprocessing/WOA2013/O_no3.nc"
use "/Users/wyao/Desktop/UVicWorking/2.9_ln_iron/dataprocessing/WOA2013/O_o2.nc"
use "/Users/wyao/Desktop/UVicWorking/2.9_ln_iron/dataprocessing/WOA2013/O_po4.nc"


def sym iscale="(-inf)(-0.5,0.5,0.1)(inf)"
def sym cscale="(-inf)(-0.6,0.6,0.1)(inf)"
def sym iscaledeni="(-inf)(-16,16,2)(inf)"
def sym cscaledeni="(-inf)(-16,16,4)(inf)"

go prepareFigure nfct_pstar_denitr
go transparency
set win 1
set vie v1
sha/nolab/set/pal=blue_orange/lev=($iscaledeni)/nokey nofe_denidiff_n0[X=30E:30E]*3600*24*365*14
ppl ylab "@SGD@SR denitrification 2100-1800"
ppl axlabp,1,-1
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscaledeni) nofe_denidiff_n0[X=30E:30E]*3600*24*365*14
go land
label -30,110,-1,0.2  (a)
set vie v2
sha/nolab/set/pal=blue_orange/lev=($iscaledeni)/nokey femask_denidiff_n0[X=30E:30E]*3600*24*365*14
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscaledeni) femask_denidiff_n0[X=30E:30E]*3600*24*365*14
go land
label -30,110,-1,0.2  (b)
set vie v3
sha/nolab/set/pal=blue_orange/lev=($iscaledeni) fedyn1_denidiff_n0[X=30E:30E]*3600*24*365*14
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscaledeni) fedyn1_denidiff_n0[X=30E:30E]*3600*24*365*14
go land
label -30,110,-1,0.2  (c)


set vie v4
sha/nolab/set/pal=blue_orange/lev=($iscale)/nokey nofe_pstar_diff[X=30E:30E,z=150]
ppl ylab "@SGD@SR P* 2100-1800 150m"
ppl axlabp,0,-1
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) nofe_pstar_diff[X=30E:30E,z=150]
go land
label -30,110,-1,0.2  (d)
set vie v5
sha/nolab/set/pal=blue_orange/lev=($iscale)/nokey femask_pstar_diff[X=30E:30E,z=150]
ppl axlabp,0,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) femask_pstar_diff[X=30E:30E,z=150]
go land
label -30,110,-1,0.2  (e)
set vie v6
sha/nolab/set/pal=blue_orange/lev=($iscale) fedyn1_pstar_diff[X=30E:30E,z=150]
ppl axlabp,0,0
ppl shakey , , , ,1
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) fedyn1_pstar_diff[X=30E:30E,z=150]
go land
label -30,110,-1,0.2  (f)


set vie v7
sha/nolab/set/pal=blue_orange/lev=($iscale)/nokey nofe_pstar_diff[X=30E:30E,z=350]
ppl ylab "@SGD@SR P* 2100-1800 euphotic 350m"
ppl xlab "NoFe"
ppl axlabp,-1,-1
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) nofe_pstar_diff[X=30E:30E,z=350]
go land
label -30,110,-1,0.2  (g)
set vie v8
sha/nolab/set/pal=blue_orange/lev=($iscale)/nokey femask_pstar_diff[X=30E:30E,z=350]
ppl xlab "Femask"
ppl axlabp,-1,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) femask_pstar_diff[X=30E:30E,z=350]
go land
label -30,110,-1,0.2  (h)
set vie v9
sha/nolab/set/pal=blue_orange/lev=($iscale) fedyn1_pstar_diff[X=30E:30E,z=350]
ppl xlab "FeDyn"
ppl axlabp,-1,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) fedyn1_pstar_diff[X=30E:30E,z=350]
go land
label -30,110,-1,0.2  (i)
go finalizeFigure_psthicken



go prepareFigure nfct_pstar_layers
go transparency
set win 2
set vie v1
sha/nolab/set/pal=blue_orange/lev=($iscale)/nokey nofe_pstar_diff[X=30E:30E,k=1:2@ave]
ppl ylab "@SGD@SR P* 2100-1800 euphotic"
ppl axlabp,1,-1
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) nofe_pstar_diff[X=30E:30E,k=1:2@ave]
go land
label -30,110,-1,0.2  (a)
set vie v2
sha/nolab/set/pal=blue_orange/lev=($iscale)/nokey femask_pstar_diff[X=30E:30E,k=1:2@ave]
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) femask_pstar_diff[X=30E:30E,k=1:2@ave]
go land
label -30,110,-1,0.2  (b)
set vie v3
sha/nolab/set/pal=blue_orange/lev=($iscale) fedyn1_pstar_diff[X=30E:30E,k=1:2@ave]
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) fedyn1_pstar_diff[X=30E:30E,k=1:2@ave]
go land
label -30,110,-1,0.2  (c)


set vie v4
sha/nolab/set/pal=blue_orange/lev=($iscale)/nokey nofe_pstar_diff[X=30E:30E,z=150]
ppl ylab "@SGD@SR P* 2100-1800 150m"
ppl axlabp,0,-1
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) nofe_pstar_diff[X=30E:30E,z=150]
go land
label -30,110,-1,0.2  (d)
set vie v5
sha/nolab/set/pal=blue_orange/lev=($iscale)/nokey femask_pstar_diff[X=30E:30E,z=150]
ppl axlabp,0,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) femask_pstar_diff[X=30E:30E,z=150]
go land
label -30,110,-1,0.2  (e)
set vie v6
sha/nolab/set/pal=blue_orange/lev=($iscale) fedyn1_pstar_diff[X=30E:30E,z=150]
ppl axlabp,0,0
ppl shakey , , , ,1
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) fedyn1_pstar_diff[X=30E:30E,z=150]
go land
label -30,110,-1,0.2  (f)


set vie v7
sha/nolab/set/pal=blue_orange/lev=($iscale)/nokey nofe_pstar_diff[X=30E:30E,z=350]
ppl ylab "@SGD@SR P* 2100-1800 euphotic 350m"
ppl xlab "NoFe"
ppl axlabp,-1,-1
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) nofe_pstar_diff[X=30E:30E,z=350]
go land
label -30,110,-1,0.2  (g)
set vie v8
sha/nolab/set/pal=blue_orange/lev=($iscale)/nokey femask_pstar_diff[X=30E:30E,z=350]
ppl xlab "Femask"
ppl axlabp,-1,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) femask_pstar_diff[X=30E:30E,z=350]
go land
label -30,110,-1,0.2  (h)
set vie v9
sha/nolab/set/pal=blue_orange/lev=($iscale) fedyn1_pstar_diff[X=30E:30E,z=350]
ppl xlab "FeDyn"
ppl axlabp,-1,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($cscale) fedyn1_pstar_diff[X=30E:30E,z=350]
go land
label -30,110,-1,0.2  (i)
go finalizeFigure_psthicken