!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=k=1:2@din
def sym op5=k=2@din
def sym op6=Z=350
def sym op7=k=3

def sym s2y=3600*24*365 !*6.625*12
def sym s2yN=3600*24*365*14.0067
def sym s2yC=3600*24*365*6.625*12

def sym s2y=3600*24*365
def sym molN2PgN=14.0067/1e15
def sym molN2TgN=14.0067/1e12
def sym opind=x=75e:110E,y=0:20n
def sym oppa=x=115w:60w,y=8s:20n
def sym opatl=x=0:20E,y=20s:0

let nofe_exp_1800=o_DETRexp[d=1,l=1,($op5)]
let nofe_exp_2100=o_DETRexp[d=1,l=300,($op5)]
let nofe_expdiff=nofe_exp_2100-nofe_exp_1800
let nofe_npp_1800=o_diaznpp[d=1,l=1,($op4)]+o_phytnpp[d=1,l=1,($op4)]
let nofe_npp_2100=o_diaznpp[d=1,l=300,($op4)]+o_phytnpp[d=1,l=300,($op4)]
let nofe_nppdiff=nofe_npp_2100-nofe_npp_1800
let nofe_o2_1800=o_o2[d=1,l=1,($op6)]*1000
let nofe_o2_2100=o_o2[d=1,l=300,($op6)]*1000
let nofe_o2diff=nofe_o2_2100-nofe_o2_1800
let nofe_no3_1800=o_no3[d=1,l=1,($op7)]
let nofe_no3_2100=o_no3[d=1,l=300,($op7)]
let nofe_no3diff=nofe_no3_2100-nofe_no3_1800
let nofe_no3fix_1800=o_no3fix[d=1,l=1,($op4)]
let nofe_no3fix_2100=o_no3fix[d=1,l=300,($op4)]
let nofe_no3fixdiff=nofe_no3fix_2100-nofe_no3fix_1800
let/bad=-999 nofe_no3fix_1800_no0= if (nofe_no3fix_1800 lt 0 or nofe_no3fix_1800 gt 0) then nofe_no3fix_1800 else (-999)
let/bad=-999 nofe_no3fixdiff_no0= if (nofe_no3fixdiff lt 0 or nofe_no3fixdiff gt 0) then nofe_no3fixdiff else (-999)
let/bad=-999 nofe_no3fix_2100_no0= if (nofe_no3fix_2100 lt 0 or nofe_no3fix_2100 gt 0) then nofe_no3fix_2100 else (-999)

let femask_exp_1800=o_DETRexp[d=2,l=1,($op5)]
let femask_exp_2100=o_DETRexp[d=2,l=300,($op5)]
let femask_expdiff=femask_exp_2100-femask_exp_1800
let femask_npp_1800=o_diaznpp[d=2,l=1,($op4)]+o_phytnpp[d=2,l=1,($op4)]
let femask_npp_2100=o_diaznpp[d=2,l=300,($op4)]+o_phytnpp[d=2,l=300,($op4)]
let femask_nppdiff=femask_npp_2100-femask_npp_1800
let femask_o2_1800=o_o2[d=2,l=1,($op6)]*1000
let femask_o2_2100=o_o2[d=2,l=300,($op6)]*1000
let femask_o2diff=femask_o2_2100-femask_o2_1800
let femask_no3_1800=o_no3[d=2,l=1,($op7)]
let femask_no3_2100=o_no3[d=2,l=300,($op7)]
let femask_no3diff=femask_no3_2100-femask_no3_1800
let femask_no3fix_1800=o_no3fix[d=2,l=1,($op4)]
let femask_no3fix_2100=o_no3fix[d=2,l=300,($op4)]
let femask_no3fixdiff=femask_no3fix_2100-femask_no3fix_1800
let/bad=-999 femask_no3fix_1800_no0= if (femask_no3fix_1800 lt 0 or femask_no3fix_1800 gt 0) then femask_no3fix_1800 else (-999)
let/bad=-999 femask_no3fixdiff_no0= if (femask_no3fixdiff lt 0 or femask_no3fixdiff gt 0) then femask_no3fixdiff else (-999)
let/bad=-999 femask_no3fix_2100_no0= if (femask_no3fix_2100 lt 0 or femask_no3fix_2100 gt 0) then femask_no3fix_2100 else (-999)

let fedyn1_exp_1800=o_DETRexp[d=3,l=1,($op5)]
let fedyn1_exp_2100=o_DETRexp[d=3,l=300,($op5)]
let fedyn1_expdiff=fedyn1_exp_2100-fedyn1_exp_1800
let fedyn1_npp_1800=o_diaznpp[d=3,l=1,($op4)]+o_phytnpp[d=3,l=1,($op4)]
let fedyn1_npp_2100=o_diaznpp[d=3,l=300,($op4)]+o_phytnpp[d=3,l=300,($op4)]
let fedyn1_nppdiff=fedyn1_npp_2100-fedyn1_npp_1800
let fedyn1_o2_1800=o_o2[d=3,l=1,($op6)]*1000
let fedyn1_o2_2100=o_o2[d=3,l=300,($op6)]*1000
let fedyn1_o2diff=fedyn1_o2_2100-fedyn1_o2_1800
let fedyn1_no3_1800=o_no3[d=3,l=1,($op7)]
let fedyn1_no3_2100=o_no3[d=3,l=300,($op7)]
let fedyn1_no3diff=fedyn1_no3_2100-fedyn1_no3_1800
let fedyn1_no3fix_1800=o_no3fix[d=3,l=1,($op4)]
let fedyn1_no3fix_2100=o_no3fix[d=3,l=300,($op4)]
let fedyn1_no3fixdiff=fedyn1_no3fix_2100-fedyn1_no3fix_1800
let/bad=-999 fedyn1_no3fix_1800_no0= if (fedyn1_no3fix_1800 lt 0 or fedyn1_no3fix_1800 gt 0) then fedyn1_no3fix_1800 else (-999)
let/bad=-999 fedyn1_no3fixdiff_no0= if (fedyn1_no3fixdiff lt 0 or fedyn1_no3fixdiff gt 0) then fedyn1_no3fixdiff else (-999)
let/bad=-999 fedyn1_no3fix_2100_no0= if (fedyn1_no3fix_2100 lt 0 or fedyn1_no3fix_2100 gt 0) then fedyn1_no3fix_2100 else (-999)

def sym iscalennp="(-inf)(0,1000,50)(inf)"
def sym iscalennpdiff="(-inf)(-250,250,25)(inf)"
def sym iscaleexp="(-inf)(0,50,2)(inf)"/sigdig=2
def sym iscaleexpdiff="(-inf)(-12,12,2)(inf)"
def sym iscaleo2="(-inf)(0,100,5)(inf)"/sigdig=2
def sym iscaleo2diff="(-inf)(-15,15,2.5)(inf)"
def sym iscaleno3="(-inf)(0,16,2)(inf)"
def sym iscaleno3diff="(-inf)(-5,5,1)(inf)"
def sym iscalenfix="(-inf)(0,5,0.25)(inf)"
def sym iscalenfixdiff="(-inf)(-3,3,0.3)(inf)"

go isquare3x3_ext.jnl
go prepareFigure nfct_nfixmap
go transparency
set win 2
set vie v1
sha/nolab/set/pal=rainbow/lev=($iscalenfix)/nokey nofe_no3fix_1800_no0[X=30E:30E]*($s2yN)
ppl ylab "N-fiaxation 1800"
ppl axlabp,1,-1
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalenfix)/nokey nofe_no3_1800[($oppa)]*1000
go land
label -30,110,-1,0.2  (a)
set vie v2
sha/nolab/set/pal=rainbow/lev=($iscalenfix)/nokey femask_no3fix_1800_no0[X=30E:30E]*($s2yN)
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($iscalenfix)/nokey femask_npp_1800[($oppa)]*1000
go land
label -30,110,-1,0.2   (b)
set vie v3
sha/nolab/set/pal=rainbow/lev=($iscalenfix) fedyn1_no3fix_1800_no0[X=30E:30E]*($s2yN)
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalenfix) fedyn1_no3_1800[($oppa)]*1000
go land
label -30,110,-1,0.2   (c)

set vie v4
sha/nolab/set/pal=rainbow/lev=($iscalenfix)/nokey nofe_no3fix_2100_no0[X=30E:30E]*($s2yN)
ppl ylab "N-fiaxation 2100"
ppl axlabp,1,-1
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalenfix)/nokey nofe_no3_2100[($oppa)]*1000
go land
label -30,110,-1,0.2  (d)
set vie v5
sha/nolab/set/pal=rainbow/lev=($iscalenfix)/nokey femask_no3fix_2100_no0[X=30E:30E]*($s2yN)
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalenfix)/nokey femask_npp_2100[($oppa)]*1000
go land
label -30,110,-1,0.2   (e)
set vie v6
sha/nolab/set/pal=rainbow/lev=($iscalenfix) fedyn1_no3fix_2100_no0[X=30E:30E]*($s2yN)
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalenfix) fedyn1_no3_2100[($oppa)]*1000
go land
label -30,110,-1,0.2   (f)


set vie v7
sha/nolab/set/pal=blue_orange/lev=($iscalenfixdiff)/nokey nofe_no3fixdiff_no0[X=30E:30E]*($s2yN)
ppl ylab "@SGD@SR N-fixation 2100-1800"
ppl xlab "NoFe"
ppl axlabp,-1,-1
ppl xaxis, 30,390,60 
ppl shade
!contour/nolab/overlay/size=0.1/pen=1/lev=($iscalenfixdiff)/nokey nofe_exp_1800[($oppa)]*($s2yC)
go land
label -30,110,-1,0.2   (g)
set vie v8
sha/nolab/set/pal=blue_orange/lev=($iscalenfixdiff)/nokey femask_no3fixdiff_no0[X=30E:30E]*($s2yN)
ppl axlabp,-1,0
ppl xlab "Femask"
ppl xaxis, 30,390,60 
ppl shade
!contour/nolab/overlay/size=0.1/pen=1/lev=($iscalenfixdiff)/nokey femask_exp_1800[($oppa)]*($s2yC)
go land
label -30,110,-1,0.2  (h)

set vie v9
sha/nolab/set/pal=blue_orange/lev=($iscalenfixdiff) fedyn1_no3fixdiff_no0[X=30E:30E]*($s2yN)
ppl axlabp,-1,0
ppl xlab "FeDyn"
ppl xaxis, 30,390,60  
ppl shakey , , , ,1
ppl shade
!contour/nolab/overlay/size=0.1/pen=1/lev=($iscalenfixdiff) fedyn1_exp_1800[($oppa)]*($s2yC)
go land
label -30,110,-1,0.2   (i)

go finalizeFigure_psthicken