!script by bencat.xuan@gmail.com



use ../model_and_output_data/noIron_spinup_fct_drift/1800_2100/tavg.01801.01.01.nc
use ../model_and_output_data/ironMask_spinup_fct_drift/1800_2100/tavg.01801.01.01.nc
use ../model_and_output_data/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 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,($op4)]*1000
let nofe_no3_2100=o_no3[d=1,l=300,($op4)]*1000
let nofe_po4_1800=o_po4[d=1,l=1,($op4)]*1000
let nofe_po4_2100=o_po4[d=1,l=300,($op4)]*1000
let nofe_nlim_1800=o_no3[d=1,l=1,($op4)]*1000/(k1n+o_no3[d=1,l=1,($op4)]*1000)
let nofe_nlim_2100=o_no3[d=1,l=300,($op4)]*1000/(k1n+o_no3[d=1,l=300,($op4)]*1000)
let nofe_plim_1800=o_po4[d=1,l=1,($op4)]*1000/(k1p+o_po4[d=1,l=1,($op4)]*1000)
let nofe_plim_2100=o_po4[d=1,l=300,($op4)]*1000/(k1p+o_po4[d=1,l=300,($op4)]*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_pstardiff=nofe_pstar_2100-nofe_pstar_1800

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,($op4)]*1000
let femask_no3_2100=o_no3[d=2,l=300,($op4)]*1000
let femask_po4_1800=o_po4[d=2,l=1,($op4)]*1000
let femask_po4_2100=o_po4[d=2,l=300,($op4)]*1000
let femask_nlim_1800=o_no3[d=2,l=1,($op4)]*1000/(k1n+o_no3[d=2,l=1,($op4)]*1000)
let femask_nlim_2100=o_no3[d=2,l=300,($op4)]*1000/(k1n+o_no3[d=2,l=300,($op4)]*1000)
let femask_plim_1800=o_po4[d=2,l=1,($op4)]*1000/(k1p+o_po4[d=2,l=1,($op4)]*1000)
let femask_plim_2100=o_po4[d=2,l=300,($op4)]*1000/(k1p+o_po4[d=2,l=300,($op4)]*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_pstardiff=femask_pstar_2100-femask_pstar_1800

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,($op4)]*1000
let fedyn1_no3_2100=o_no3[d=3,l=300,($op4)]*1000
let fedyn1_po4_1800=o_po4[d=3,l=1,($op4)]*1000
let fedyn1_po4_2100=o_po4[d=3,l=300,($op4)]*1000
let fedyn1_nlim_1800=o_no3[d=3,l=1,($op4)]*1000/(k1n+o_no3[d=3,l=1,($op4)]*1000)
let fedyn1_nlim_2100=o_no3[d=3,l=300,($op4)]*1000/(k1n+o_no3[d=3,l=300,($op4)]*1000)
let fedyn1_plim_1800=o_po4[d=3,l=1,($op4)]*1000/(k1p+o_po4[d=3,l=1,($op4)]*1000)
let fedyn1_plim_2100=o_po4[d=3,l=300,($op4)]*1000/(k1p+o_po4[d=3,l=300,($op4)]*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_pstardiff=fedyn1_pstar_2100-fedyn1_pstar_1800

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)

def sym iscale="(-inf)(-0.5,0.5,0.1)(inf)"
def sym cscale="(-inf)(-0.6,0.6,0.1)(inf)"




go ilandscape3x1_ext.jnl




let nofe_pstar_diff=nofe_pstar_2100-nofe_pstar_1800
let femask_pstar_diff=femask_pstar_2100-femask_pstar_1800
let fedyn1_pstar_diff=fedyn1_pstar_2100-fedyn1_pstar_1800

def sym iscale="(-inf)(-0.5,0.5,0.1)(inf)"
def sym cscale="(-inf)(-0.6,0.6,0.1)(inf)"

go prepareFigure Fig5
go transparency
set win 2
set vie v1
sha/nolab/set/pal=blue_orange/lev=($iscale)/nokey nofe_pstar_diff[X=30E:30E]
ppl ylab "@SGD@SR P^@PM22@SR 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=($cscale) nofe_pstar_diff[X=30E:30E]
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]
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]
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]
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]
go land
label -30,110,-1,0.2 (c)
go finalizeFigure_psthicken