!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 s2D=3600*24
def sym s2DN=3600*24* 14.0067
def sym s2DC=3600*24*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 dummy=o_phyt[d=1,l=1]*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_phytnpp_1800=o_phytnpp[d=1,l=1,($op4)]
let nofe_phytnpp_2100=o_phytnpp[d=1,l=300,($op4)]
let nofe_phytnppdiff=nofe_phytnpp_2100-nofe_phytnpp_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 nofe_diazgraz_1800=if o_diaznpp[d=1,l=1,k=1:2@ave] gt 0 then o_diazgraz[d=1,l=1,k=1:2@ave]/o_diaz[d=1,l=1,k=1:2@ave]*($s2D) else dummy[k=1:2@ave,l=1]
let nofe_diazgraz_2100=if o_diaznpp[d=1,l=300,k=1:2@ave] gt 0 then o_diazgraz[d=1,l=300,k=1:2@ave]/o_diaz[d=1,l=300,k=1:2@ave]*($s2D) else dummy[k=1:2@ave,l=300]
let nofe_phytgraz_1800=if o_phytnpp[d=1,l=1,k=1:2@ave] gt 0 then o_phytgraz[d=1,l=1,k=1:2@ave]/o_phyt[d=1,l=1,k=1:2@ave]*($s2D) else dummy[k=1:2@ave,l=1]
let nofe_phytgraz_2100=if o_phytnpp[d=1,l=300,k=1:2@ave] gt 0 then o_phytgraz[d=1,l=300,k=1:2@ave]/o_phyt[d=1,l=300,k=1:2@ave]*($s2D) else dummy[k=1:2@ave,l=300]
let nofe_diazh=if nofe_diazgraz_1800+nofe_diazgraz_2100 gt 0 then 1 else dummy[k=1:2@ave]
let nofe_prefgraz_1800= if nofe_diazh gt 0 then nofe_phytgraz_1800-nofe_diazgraz_1800 else 0
let nofe_prefgraz_2100= if nofe_diazh gt 0 then nofe_phytgraz_2100-nofe_diazgraz_2100 else 0
let/bad=-999 nofe_prefgrazdiff= if (nofe_prefgraz_2100-nofe_prefgraz_1800) gt 0 or (nofe_prefgraz_2100-nofe_prefgraz_1800) lt 0  then (((nofe_prefgraz_2100)^2)^0.5-((nofe_prefgraz_1800)^2)^0.5) else (-999)
let nofe_diaznpp_1800=o_diaznpp[d=1,l=1,($op4)]
let nofe_diaznpp_2100=o_diaznpp[d=1,l=300,($op4)]
let/bad=-999 nofe_diaznppdiff= if nofe_diaznpp_1800 gt 0 or nofe_diaznpp_2100 gt 0 then nofe_diaznpp_2100-nofe_diaznpp_1800 else (-999)
let nofe_diaz_1800=o_diaz[d=1,l=1,($op4)]
let nofe_diaz_2100=o_diaz[d=1,l=300,($op4)]
let/bad=-999 nofe_diazdiff= if nofe_diaz_1800 gt 0 or nofe_diaz_2100 gt 0 then nofe_diaz_2100-nofe_diaz_1800 else (-999)
let nofe_nsource_1800=nofe_no3fix_1800/nofe_diaznpp_1800
let nofe_nsource_2100=nofe_no3fix_2100/nofe_diaznpp_2100
let nofe_nsourcediff=nofe_nsource_2100-nofe_nsource_1800
let nofe_prefgrazdiff_w0=(nofe_prefgraz_2100-nofe_prefgraz_1800)
let nofe_diaznppdiff_w0=nofe_diaznpp_2100-nofe_diaznpp_1800
let nofe_diazdiff_w0=nofe_diaz_2100-nofe_diaz_1800 


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_phytnpp_1800=o_phytnpp[d=2,l=1,($op4)]
let femask_phytnpp_2100=o_phytnpp[d=2,l=300,($op4)]
let femask_phytnppdiff=femask_phytnpp_2100-femask_phytnpp_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 femask_diazgraz_1800=if o_diaznpp[d=2,l=1,k=1:2@ave] gt 0 then o_diazgraz[d=2,l=1,k=1:2@ave]/o_diaz[d=2,l=1,k=1:2@ave]*($s2D) else dummy[k=1:2@ave,l=1]
let femask_diazgraz_2100=if o_diaznpp[d=2,l=300,k=1:2@ave] gt 0 then o_diazgraz[d=2,l=300,k=1:2@ave]/o_diaz[d=2,l=300,k=1:2@ave]*($s2D) else dummy[k=1:2@ave,l=300]
let femask_phytgraz_1800=if o_phytnpp[d=2,l=1,k=1:2@ave] gt 0 then o_phytgraz[d=2,l=1,k=1:2@ave]/o_phyt[d=2,l=1,k=1:2@ave]*($s2D) else dummy[k=1:2@ave,l=1]
let femask_phytgraz_2100=if o_phytnpp[d=2,l=300,k=1:2@ave] gt 0 then o_phytgraz[d=2,l=300,k=1:2@ave]/o_phyt[d=2,l=300,k=1:2@ave]*($s2D) else dummy[k=1:2@ave,l=300]
let femask_diazh=if femask_diazgraz_1800+femask_diazgraz_2100 gt 0 then 1 else dummy[k=1:2@ave]
let femask_prefgraz_1800= if femask_diazh gt 0 then femask_phytgraz_1800-femask_diazgraz_1800 else 0
let femask_prefgraz_2100= if femask_diazh gt 0 then femask_phytgraz_2100-femask_diazgraz_2100 else 0
let/bad=-999 femask_prefgrazdiff= if (femask_prefgraz_2100-femask_prefgraz_1800) gt 0 or (femask_prefgraz_2100-femask_prefgraz_1800) lt 0  then (((femask_prefgraz_2100)^2)^0.5-((femask_prefgraz_1800)^2)^0.5) else (-999)
let femask_diaznpp_1800=o_diaznpp[d=2,l=1,($op4)]
let femask_diaznpp_2100=o_diaznpp[d=2,l=300,($op4)]
let/bad=-999 femask_diaznppdiff= if femask_diaznpp_1800 gt 0 or femask_diaznpp_2100 gt 0 then femask_diaznpp_2100-femask_diaznpp_1800 else (-999)
let femask_diaz_1800=o_diaz[d=2,l=1,($op4)]
let femask_diaz_2100=o_diaz[d=2,l=300,($op4)]
let/bad=-999 femask_diazdiff= if femask_diaz_1800 gt 0 or femask_diaz_2100 gt 0 then femask_diaz_2100-femask_diaz_1800 else (-999)
let femask_nsource_1800=femask_no3fix_1800/femask_diaznpp_1800
let femask_nsource_2100=femask_no3fix_2100/femask_diaznpp_2100
let femask_nsourcediff=femask_nsource_2100-femask_nsource_1800
let femask_prefgrazdiff_w0=(femask_prefgraz_2100-femask_prefgraz_1800)
let femask_diaznppdiff_w0=femask_diaznpp_2100-femask_diaznpp_1800
let femask_diazdiff_w0=femask_diaz_2100-femask_diaz_1800 


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_phytnpp_1800=o_phytnpp[d=3,l=1,($op4)]
let fedyn1_phytnpp_2100=o_phytnpp[d=3,l=300,($op4)]
let fedyn1_phytnppdiff=fedyn1_phytnpp_2100-fedyn1_phytnpp_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 fedyn1_diazgraz_1800=if o_diaznpp[d=3,l=1,k=1:2@ave] gt 0 then o_diazgraz[d=3,l=1,k=1:2@ave]/o_diaz[d=3,l=1,k=1:2@ave]*($s2D) else dummy[k=1:2@ave,l=1]
let fedyn1_diazgraz_2100=if o_diaznpp[d=3,l=300,k=1:2@ave] gt 0 then o_diazgraz[d=3,l=300,k=1:2@ave]/o_diaz[d=3,l=300,k=1:2@ave]*($s2D) else dummy[k=1:2@ave,l=300]
let fedyn1_phytgraz_1800=if o_phytnpp[d=3,l=1,k=1:2@ave] gt 0 then o_phytgraz[d=3,l=1,k=1:2@ave]/o_phyt[d=3,l=1,k=1:2@ave]*($s2D) else dummy[k=1:2@ave,l=1]
let fedyn1_phytgraz_2100=if o_phytnpp[d=3,l=300,k=1:2@ave] gt 0 then o_phytgraz[d=3,l=300,k=1:2@ave]/o_phyt[d=3,l=300,k=1:2@ave]*($s2D) else dummy[k=1:2@ave,l=300]
let fedyn1_diazh=if fedyn1_diazgraz_1800+fedyn1_diazgraz_2100 gt 0 then 1 else dummy[k=1:2@ave]
let fedyn1_prefgraz_1800= if fedyn1_diazh gt 0 then fedyn1_phytgraz_1800-fedyn1_diazgraz_1800 else 0
let fedyn1_prefgraz_2100= if fedyn1_diazh gt 0 then fedyn1_phytgraz_2100-fedyn1_diazgraz_2100 else 0
let/bad=-999 fedyn1_prefgrazdiff= if (fedyn1_prefgraz_2100-fedyn1_prefgraz_1800) gt 0 or (fedyn1_prefgraz_2100-fedyn1_prefgraz_1800) lt 0  then (((fedyn1_prefgraz_2100)^2)^0.5-((fedyn1_prefgraz_1800)^2)^0.5) else (-999)
let fedyn1_diaznpp_1800=o_diaznpp[d=3,l=1,($op4)]
let fedyn1_diaznpp_2100=o_diaznpp[d=3,l=300,($op4)]
let/bad=-999 fedyn1_diaznppdiff= if fedyn1_diaznpp_1800 gt 0 or fedyn1_diaznpp_2100 gt 0 then fedyn1_diaznpp_2100-fedyn1_diaznpp_1800 else (-999)
let fedyn1_diaz_1800=o_diaz[d=3,l=1,($op4)]
let fedyn1_diaz_2100=o_diaz[d=3,l=300,($op4)]
let/bad=-999 fedyn1_diazdiff= if fedyn1_diaz_1800 gt 0 or fedyn1_diaz_2100 gt 0 then fedyn1_diaz_2100-fedyn1_diaz_1800 else (-999)
let fedyn1_nsource_1800=fedyn1_no3fix_1800/fedyn1_diaznpp_1800
let fedyn1_nsource_2100=fedyn1_no3fix_2100/fedyn1_diaznpp_2100
let fedyn1_nsourcediff=fedyn1_nsource_2100-fedyn1_nsource_1800
let fedyn1_prefgrazdiff_w0=(fedyn1_prefgraz_2100-fedyn1_prefgraz_1800)
let fedyn1_diaznppdiff_w0=fedyn1_diaznpp_2100-fedyn1_diaznpp_1800
let fedyn1_diazdiff_w0=fedyn1_diaz_2100-fedyn1_diaz_1800 


def sym iscalennp="(-inf)(0,1000,50)(inf)"
def sym iscalennpdiff="(-inf)(-250,250,25)(inf)"
def sym iscalediaznnp="(-inf)(0,100,5)(inf)"
!def sym iscalediaznppdiff="(-inf)(-25,25,2.5)(inf)"
def sym iscalediaznppdiff="(-inf)(-40,40,8)(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)"
!def sym iscaleperf="(-inf)(-0.2,0.2,0.02)(inf)"
def sym iscaleperf="(-inf)(-0.3,0.3,0.025)(inf)"

!def sym iscalediaznppdiffc="(-inf)(-25,25,5)(inf)"
def sym iscalediaznppdiffc="(-inf)(-40,40,8)(inf)"
def sym iscalensource="(-inf)(0,100,5)(inf)"
go isquare3x3_ext.jnl


go prepareFigure nfct_grazingcontrol_nsource
go transparency
set win 2
set vie v1
sha/nolab/set/pal=blue_orange/lev=($iscalediaznppdiff)/nokey nofe_diaznppdiff[X=30E:30E]*($s2yC)
ppl ylab "@SGD@SR diaz. NPP 2100-1800"
ppl axlabp,1,-1
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey nofe_diaznppdiff_w0*($s2yC)
go land
label -30,110,-1,0.2 (a)
set vie v2
sha/nolab/set/pal=blue_orange/lev=($iscalediaznppdiff)/nokey femask_diaznppdiff[X=30E:30E]*($s2yC)
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey femask_diaznppdiff_w0*($s2yC)
go land
label -30,110,-1,0.2  (b)
set vie v3
sha/nolab/set/pal=blue_orange/lev=($iscalediaznppdiff) fedyn1_diaznppdiff[X=30E:30E]*($s2yC)
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc) fedyn1_diaznppdiff_w0*($s2yC)
go land
label -30,110,-1,0.2  (c)
let/bad=-999 nofe_prefgraz_2100n0= if nofe_prefgraz_2100 gt 0 or nofe_prefgraz_2100 lt 0 then nofe_prefgraz_2100 else (-999)
let/bad=-999 femask_prefgraz_2100n0=if femask_prefgraz_2100 gt 0 or femask_prefgraz_2100 lt 0 then femask_prefgraz_2100 else (-999)
let/bad=-999 fedyn1_prefgraz_2100n0=if fedyn1_prefgraz_2100 gt 0 or fedyn1_prefgraz_2100 lt 0 then fedyn1_prefgraz_2100 else (-999)
set vie v4
sha/nolab/set/pal=blue_orange/lev=($iscaleperf)/nokey nofe_prefgraz_2100n0[X=30E:30E]
ppl ylab "graz. pref. on phyt. 2100"
ppl xlab "NoFe"
ppl axlabp,0,-1
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey nofe_diaznppdiff_w0*($s2yC)
go land
label -30,110,-1,0.2 (d)
set vie v5
sha/nolab/set/pal=blue_orange/lev=($iscaleperf)/nokey femask_prefgraz_2100n0[X=30E:30E]
ppl axlabp,0,0
ppl xlab "Femask"
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey femask_diaznppdiff_w0*($s2yC)
go land
label -30,110,-1,0.2  (e)

set vie v6
sha/nolab/set/pal=blue_orange/lev=($iscaleperf) fedyn1_prefgraz_2100n0[X=30E:30E]
ppl axlabp,0,0
ppl xlab "FeDyn"
ppl xaxis, 30,390,60  
ppl shakey , , , ,1
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc) fedyn1_diaznppdiff_w0*($s2yC)
go land
label -30,110,-1,0.2  (f)


set vie v7
sha/nolab/set/pal=blue_orange/lev=($iscalensource)/nokey nofe_nsource_2100[X=30E:30E]*100
ppl ylab "diaz dependence on N fixation"
ppl xlab "NoFe"
ppl axlabp,-1,-1
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey nofe_diaznppdiff_w0*($s2yC)
go land
label -30,110,-1,0.2  (g)
set vie v8
sha/nolab/set/pal=blue_orange/lev=($iscalensource)/nokey femask_nsource_2100[X=30E:30E]*100
ppl xlab "Femask"
ppl axlabp,-1,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey femask_diaznppdiff_w0*($s2yC)
go land
label -30,110,-1,0.2  (h)
set vie v9
sha/nolab/set/pal=blue_orange/lev=($iscalensource) fedyn1_nsource_2100[X=30E:30E]*100
ppl xlab "FeDyn"
ppl axlabp,-1,0
ppl xaxis, 30,390,60 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc) fedyn1_diaznppdiff_w0*($s2yC)
go land
label -30,110,-1,0.2  (i)


go finalizeFigure_psthicken

go prepareFigure nfct_nppmap
go transparency
set win 3

set vie v1
sha/nolab/set/pal=blue_orange/lev=($iscalennp)/nokey nofe_npp_1800[X=30E:30E]*($s2yC)
ppl ylab "NPP 1800"
ppl axlabp,1,-1
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey nofe_diaznppdiff_w0
go land
label -30,110,-1,0.2 (a)
set vie v2
sha/nolab/set/pal=blue_orange/lev=($iscalennp)/nokey femask_npp_1800[X=30E:30E]*($s2yC)
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey femask_diaznppdiff_w0
go land
label -30,110,-1,0.2  (b)
set vie v3
sha/nolab/set/pal=blue_orange/lev=($iscalennp) fedyn1_npp_1800[X=30E:30E]*($s2yC)
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc) fedyn1_diaznppdiff_w0
go land
label -30,110,-1,0.2  (c)


set vie v4
sha/nolab/set/pal=blue_orange/lev=($iscalennp)/nokey nofe_npp_2100[X=30E:30E]*($s2yC)
ppl ylab "NPP 2100"
ppl xlab "NoFe"
ppl axlabp,0,-1
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey nofe_diaznppdiff_w0
go land
label -30,110,-1,0.2 (d)
set vie v5
sha/nolab/set/pal=blue_orange/lev=($iscalennp)/nokey femask_npp_2100[X=30E:30E]*($s2yC)
ppl axlabp,0,0
ppl xlab "Femask"
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey femask_diaznppdiff_w0
go land
label -30,110,-1,0.2  (e)

set vie v6
sha/nolab/set/pal=blue_orange/lev=($iscalennp) fedyn1_npp_2100[X=30E:30E]*($s2yC)
ppl axlabp,0,0
ppl xlab "FeDyn"
ppl xaxis, 30,390,60  
ppl shakey , , , ,1
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc) fedyn1_diaznppdiff_w0
go land
label -30,110,-1,0.2  (f)

set vie v7
sha/nolab/set/pal=blue_orange/lev=($iscalennpdiff)/nokey nofe_nppdiff[X=30E:30E]*($s2yC)
ppl ylab "@SGD@SR NPP 2100-1800"
ppl xlab "NoFe"
ppl axlabp,-1,-1
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey nofe_diaznppdiff_w0
go land
label -30,110,-1,0.2  (g)
set vie v8
sha/nolab/set/pal=blue_orange/lev=($iscalennpdiff)/nokey femask_nppdiff[X=30E:30E]*($s2yC)
ppl xlab "Femask"
ppl axlabp,-1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey femask_diaznppdiff_w0
go land
label -30,110,-1,0.2  (h)
set vie v9
sha/nolab/set/pal=blue_orange/lev=($iscalennpdiff) fedyn1_nppdiff[X=30E:30E]*($s2yC)
ppl xlab "FeDyn"
ppl axlabp,-1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc) fedyn1_diaznppdiff_w0
go land
label -30,110,-1,0.2  (i)

go finalizeFigure_psthicken

go prepareFigure nfct_diaznppmap
go transparency
set win 4

set vie v1
sha/nolab/set/pal=blue_orange/lev=($iscalediaznnp)/nokey nofe_diaznpp_1800[X=30E:30E]*($s2yC)
ppl ylab "Diaz PP 1800"
ppl axlabp,1,-1
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey nofe_diaznppdiff_w0
go land
label -30,110,-1,0.2 (a)
set vie v2
sha/nolab/set/pal=blue_orange/lev=($iscalediaznnp)/nokey femask_diaznpp_1800[X=30E:30E]*($s2yC)
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey femask_diaznppdiff_w0
go land
label -30,110,-1,0.2  (b)
set vie v3
sha/nolab/set/pal=blue_orange/lev=($iscalediaznnp) fedyn1_diaznpp_1800[X=30E:30E]*($s2yC)
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc) fedyn1_diaznppdiff_w0
go land
label -30,110,-1,0.2  (c)

set vie v4
sha/nolab/set/pal=blue_orange/lev=($iscalediaznnp)/nokey nofe_diaznpp_2100[X=30E:30E]*($s2yC)
ppl ylab "Diaz PP 2100"
ppl xlab "NoFe"
ppl axlabp,0,-1
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey nofe_diaznppdiff_w0
go land
label -30,110,-1,0.2 (d)
set vie v5
sha/nolab/set/pal=blue_orange/lev=($iscalediaznnp)/nokey femask_diaznpp_2100[X=30E:30E]*($s2yC)
ppl axlabp,0,0
ppl xlab "Femask"
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey femask_diaznppdiff_w0
go land
label -30,110,-1,0.2  (e)

set vie v6
sha/nolab/set/pal=blue_orange/lev=($iscalediaznnp) fedyn1_diaznpp_2100[X=30E:30E]*($s2yC)
ppl axlabp,0,0
ppl xlab "FeDyn"
ppl xaxis, 30,390,60  
ppl shakey , , , ,1
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc) fedyn1_diaznppdiff_w0
go land
label -30,110,-1,0.2  (f)

set vie v7
sha/nolab/set/pal=blue_orange/lev=($iscalediaznppdiff)/nokey nofe_diaznppdiff[X=30E:30E]*($s2yC)
ppl ylab "@SGD@SR Diaz PP 2100-1800"
ppl xlab "NoFe"
ppl axlabp,-1,-1
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey nofe_diaznppdiff_w0
go land
label -30,110,-1,0.2  (g)
set vie v8
sha/nolab/set/pal=blue_orange/lev=($iscalediaznppdiff)/nokey femask_diaznppdiff[X=30E:30E]*($s2yC)
ppl xlab "Femask"
ppl axlabp,-1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey femask_diaznppdiff_w0
go land
label -30,110,-1,0.2  (h)
set vie v9
sha/nolab/set/pal=blue_orange/lev=($iscalediaznppdiff) fedyn1_diaznppdiff[X=30E:30E]*($s2yC)
ppl xlab "FeDyn"
ppl axlabp,-1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc) fedyn1_diaznppdiff_w0
go land
label -30,110,-1,0.2  (i)

go finalizeFigure_psthicken

go prepareFigure nfct_phytnppmap
go transparency
set win 5

set vie v1
sha/nolab/set/pal=blue_orange/lev=($iscalennp)/nokey nofe_phytnpp_1800[X=30E:30E]*($s2yC)
ppl ylab "Phyt NPP 1800"
ppl axlabp,1,-1
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey nofe_diaznppdiff_w0
go land
label -30,110,-1,0.2 (a)
set vie v2
sha/nolab/set/pal=blue_orange/lev=($iscalennp)/nokey femask_phytnpp_1800[X=30E:30E]*($s2yC)
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey femask_diaznppdiff_w0
go land
label -30,110,-1,0.2  (b)
set vie v3
sha/nolab/set/pal=blue_orange/lev=($iscalennp) fedyn1_phytnpp_1800[X=30E:30E]*($s2yC)
ppl axlabp,1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc) fedyn1_diaznppdiff_w0
go land
label -30,110,-1,0.2  (c)

!contour/overlay/size=0.1/pen=1/nolab/lev=($iscaleperf)/nokey nofe_prefgrazdiff[X=30E:30E]
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscaleperf)/nokey femask_prefgrazdiff[X=30E:30E]
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscaleperf) fedyn1_prefgrazdiff[X=30E:30E]

set vie v4
sha/nolab/set/pal=blue_orange/lev=($iscalennp)/nokey nofe_phytnpp_2100[X=30E:30E]*($s2yC)
ppl ylab "Phyt NPP 2100"
ppl xlab "NoFe"
ppl axlabp,0,-1
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey nofe_diaznppdiff_w0
go land
label -30,110,-1,0.2 (d)
set vie v5
sha/nolab/set/pal=blue_orange/lev=($iscalennp)/nokey femask_phytnpp_2100[X=30E:30E]*($s2yC)
ppl axlabp,0,0
ppl xlab "Femask"
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey femask_diaznppdiff_w0
go land
label -30,110,-1,0.2  (e)

set vie v6
sha/nolab/set/pal=blue_orange/lev=($iscalennp) fedyn1_phytnpp_2100[X=30E:30E]*($s2yC)
ppl axlabp,0,0
ppl xlab "FeDyn"
ppl xaxis, 30,390,60  
ppl shakey , , , ,1
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc) fedyn1_diaznppdiff_w0
go land
label -30,110,-1,0.2  (f)

set vie v7
sha/nolab/set/pal=blue_orange/lev=($iscalennpdiff)/nokey nofe_phytnppdiff[X=30E:30E]*($s2yC)
ppl ylab "@SGD@SR Phyt NPP 2100-1800"
ppl xlab "NoFe"
ppl axlabp,-1,-1
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey nofe_diaznppdiff_w0
go land
label -30,110,-1,0.2  (g)
set vie v8
sha/nolab/set/pal=blue_orange/lev=($iscalennpdiff)/nokey femask_phytnppdiff[X=30E:30E]*($s2yC)
ppl xlab "Femask"
ppl axlabp,-1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc)/nokey femask_diaznppdiff_w0
go land
label -30,110,-1,0.2  (h)
set vie v9
sha/nolab/set/pal=blue_orange/lev=($iscalennpdiff) fedyn1_phytnppdiff[X=30E:30E]*($s2yC)
ppl xlab "FeDyn"
ppl axlabp,-1,0
ppl xaxis, 30,390,60 
ppl shade
!contour/overlay/size=0.1/pen=1/nolab/lev=($iscalediaznppdiffc) fedyn1_diaznppdiff_w0
go land
label -30,110,-1,0.2  (i)

go finalizeFigure_psthicken