!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 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 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 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


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)"


go isquare3x3_ext.jnl
go prepareFigure nfct_nppexpo2_overlay_diff
go transparency
set win 2
set vie v1
sha/nolab/set/pal=blue_orange/lev=($iscalennpdiff)/nokey nofe_nppdiff[($oppa)]*($s2yC)
ppl ylab "diff NPP 2100-1800"
ppl axlabp,1,-1
ppl xaxis, 245,300,10
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($iscalennp)/nokey nofe_npp_1800[($oppa)]*($s2yC)
go land
label 235,22,-1,0.2 (a)
set vie v2
sha/nolab/set/pal=blue_orange/lev=($iscalennpdiff)/nokey femask_nppdiff[($oppa)]*($s2yC)
ppl axlabp,1,0
ppl xaxis, 245,300,10
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($iscalennp)/nokey femask_npp_1800[($oppa)]*($s2yC)
go land
label 235,22,-1,0.2 (b)
set vie v3
sha/nolab/set/pal=blue_orange/lev=($iscalennpdiff) fedyn1_nppdiff[($oppa)]*($s2yC)
ppl axlabp,1,0
ppl xaxis, 245,300,10 
ppl shade
contour/overlay/size=0.1/pen=1/nolab/lev=($iscalennp) fedyn1_npp_1800[($oppa)]*($s2yC)
go land
label 235,22,-1,0.2 (c)
set vie v4
sha/nolab/set/pal=blue_orange/lev=($iscaleexpdiff)/nokey nofe_expdiff[($oppa)]*($s2yC)
ppl ylab "diff EXP 2100-1800"
ppl axlabp,0,-1
ppl xaxis, 245,300,10 
ppl shade
contour/nolab/overlay/size=0.1/pen=1/lev=($iscaleexp)/nokey nofe_exp_1800[($oppa)]*($s2yC)
go land
label 235,22,-1,0.2 (d)
set vie v5
sha/nolab/set/pal=blue_orange/lev=($iscaleexpdiff)/nokey femask_expdiff[($oppa)]*($s2yC)
ppl axlabp,0,0
ppl xaxis, 245,300,10 
ppl shade
contour/nolab/overlay/size=0.1/pen=1/lev=($iscaleexp)/nokey femask_exp_1800[($oppa)]*($s2yC)
go land
label 235,22,-1,0.2 (e)

set vie v6
sha/nolab/set/pal=blue_orange/lev=($iscaleexpdiff) fedyn1_expdiff[($oppa)]*($s2yC)
ppl axlabp,0,0
ppl xaxis, 245,300,10  
ppl shakey , , , ,1
ppl shade
contour/nolab/overlay/size=0.1/pen=1/lev=($iscaleexp) fedyn1_exp_1800[($oppa)]*($s2yC)
go land
label 235,22,-1,0.2 (f)
set vie v7
sha/nolab/set/pal=blue_orange/lev=($iscaleo2diff)/nokey nofe_o2diff[($oppa)]
ppl ylab "diff O_2 2100-1800"
ppl xlab "NoFe"
ppl axlabp,-1,-1
ppl xaxis, 245,300,10 
ppl shade
contour/overlay/nolab/size=0.1/pen=1/lev=($iscaleo2)/nokey nofe_o2_1800[($oppa)]
go land
label 235,22,-1,0.2 (g)
set vie v8
sha/nolab/set/pal=blue_orange/lev=($iscaleo2diff)/nokey femask_o2diff[($oppa)]
ppl xlab "Femask"
ppl axlabp,-1,0
ppl xaxis, 245,300,10 
ppl shade
contour/overlay/nolab/size=0.1/pen=1/lev=($iscaleo2)/nokey femask_o2_1800[($oppa)]
go land
label 235,22,-1,0.2 (h)

set vie v9
sha/nolab/set/pal=blue_orange/lev=($iscaleo2diff) fedyn1_o2diff[($oppa)]
ppl xlab "FeDyn"
ppl axlabp,-1,0
ppl xaxis, 245,300,10 
ppl shade
contour/overlay/nolab/size=0.1/pen=1/lev=($iscaleo2) fedyn1_o2_1800[($oppa)]
go land
label 235,22,-1,0.2 (i)
go finalizeFigure_psthicken


let/bad=-999 dfe_pa= if g_mskhr eq 2 then o_dfe[d=3] else (-999)
let/bad=-999 dfe_at= if g_mskhr eq 1 then o_dfe[d=3] else (-999)
let/bad=-999 dfe_in= if g_mskhr eq 3 then o_dfe[d=3] else (-999)
set win 5
sha dfe_pa[k=1:2@ave,($oppa)]
list dfe_pa[k=1:2@ave,x=115w:60w@ave,y=8s:20n@ave,l=300]/dfe_pa[k=1:2@ave,x=115w:60w@ave,y=8s:20n@ave,l=1]
list dfe_at[k=1:2@ave,x=0:20E@ave,y=20s:0@ave,l=300]/dfe_pa[k=1:2@ave,x=0:20E@ave,y=20s:0@ave,l=1]

list dfe_in[k=1:2@ave,x=75e:110E@ave,y=0:20n@ave,l=300]/dfe_in[k=1:2@ave,x=75e:110E@ave,y=0:20n@ave,l=1]