%%(8) sensitivity tests
%% need to copy no.nc, nosh.nc and nosst.nc into this folder

%neither
ferret -nojnl
cancel data/all
set mode metafile
sp rm metafile.plt*

go margins 1.5 1.5 1.5 2.0

use inf_omega.nc
use 10_omega.nc

sha/set/nolab/palette=blue_orange/lev=(-inf)(1.2,9.6,0.4)(inf) omega[x=@ave,d=1]

ppl txnmtc -1
ppl labset,0.16,0.2,0.2,0.2
ppl shakey 1 1 .2
ppl axlsze 0.18,0.18
ppl txlsze 0.14,0
ppl pen 0 1 
ppl shade

contour/nolab/ov/lev=(9.0)/color=black/size=0.2/spacing=2 omega[x=@ave,d=1]


cancel mode metafile
sp Fprint -l cps -o inf_omega.ps -R metafile.plt

exit

ps_thicken inf_omega.ps 5 > inf_omega_1.ps
rm inf_omega.ps






ferret -nojnl
cancel data/all
set mode metafile
sp rm metafile.plt*

go margins 1.5 1.5 1.5 2.0

use inf_omega.nc
use 10_omega.nc

sha/set/nolab/palette=blue_orange/lev=(-inf)(1.2,9.6,0.4)(inf) omega[x=@ave,d=2]

ppl txnmtc -1
ppl labset,0.16,0.2,0.2,0.2
ppl shakey 1 1 .2
ppl axlsze 0.18,0.18
ppl txlsze 0.14,0
ppl pen 0 1 
ppl shade

contour/nolab/ov/lev=(9.0)/color=black/size=0.2/spacing=2 omega[x=@ave,d=2]


cancel mode metafile
sp Fprint -l cps -o 10_omega.ps -R metafile.plt

exit

ps_thicken 10_omega.ps 5 > 10_omega_1.ps
rm 10_omega.ps




ferret -nojnl
cancel data/all
set mode metafile
sp rm metafile.plt*

go margins 1.5 1.5 1.5 2.0

use oliv10_9_tavg.nc
use olivinf_9_tavg.nc

let omega=o_oaragsur[d=1]-o_oaragsur[d=2]

sha/set/nolab/palette=dynamic_cmyk/lev=(-inf)(0,3,0.2)(inf) omega[l=80]

ppl txnmtc -1
ppl labset,0.16,0.2,0.2,0.2
ppl shakey 1 1 .2
ppl axlsze 0.18,0.18
ppl txlsze 0.14,0
ppl pen 0 1 
ppl shade

let/bad=-999 msk = if g_mskt[d=1] eq 0 then 0 else (-999)
shade/ov/nokey/palette=grey_dark/nolab/levels="(-10,1,0.1)" msk


cancel mode metafile
sp Fprint -l cps -o omega_dif_shade.ps -R metafile.plt

exit

ps_thicken omega_dif_shade.ps 5 > omega_dif.ps
rm omega_dif_shade.ps