%% This is going to generate shading figures in MS. There are quite a lot to generate


% (1)
%-------------------------alk distribution

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

go margins 1.5 1.5 1.5 2.0

use tavg.02021.01.01.nc
use mask_c.nc
use control.nc

let alk_add_add=f_alk[d=1]*mask[d=2]*365*24*3600
let alk_add=alk_add_add[l=@ave]*80
sha/set/nolab/palette=dynamic_cmyk/lev=(-inf)(0,1600,100)(inf) alk_add


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 1
contour/nokey/size=0.00000001/ov/lev=3V,1.0max/color=black msk


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

exit

ps_thicken f_alk_shade.ps 5 > cummlative_alk.ps
rm f_alk_shade.ps


%%% Can use "go ccbar" for changing color bar properties

%(2)
%-------------------------omega distribution

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

go margins 1.5 1.5 1.5 2.0

use tavg.02021.01.01.nc
use mask_c.nc
use control.nc

sha/set/nolab/palette=dynamic_cmyk/lev=(-inf)(0,4,0.4)(inf)  o_oaragsur[l=80,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


let/bad=-999 msk = if g_mskt[d=1] eq 0 then 0 else 1
contour/nokey/size=0.00000001/ov/lev=3V,1.0max/color=black msk


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

exit

ps_thicken omega_shade.ps 5 > omega_2099_oliv100.ps
rm omega_shade.ps


%(3)
%-------------------------omega_dif distribution
ferret -nojnl
cancel data/all
set mode metafile
sp rm metafile.plt*

go margins 1.5 1.5 1.5 2.0

use tavg.02021.01.01.nc
use mask_c.nc
use control.nc

let omega_dif=O_OARAGSUR[d=1]-O_OARAGSUR[d=3]
sha/set/nolab/palette=dynamic_cmyk/lev=(-inf)(0,3.8,0.4)(inf) omega_dif[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 1
contour/nokey/size=0.00000001/ov/lev=3V,1.0max/color=black 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_diff_oliv100.ps
rm omega_dif_shade.ps

%(4)
%-------------------------dic_dif distribution

ferret -nojnl

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

go margins 1.5 1.5 1.5 2.0

use tavg.02021.01.01.nc
use mask_c.nc
use control.nc
let dic_dif_dif=f_dic[d=1]*356*24*3600-f_dic[d=3]*356*24*3600
let dic_dif=dic_dif_dif[l=@ave]*80
sha/set/nolab/palette=dynamic_cmyk/lev=c/lev=(-inf)(-100,1700,100)(inf) dic_dif[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 1
contour/nokey/size=0.00000001/ov/lev=3V,1.0max/color=black msk

%contour/nolab/ov/lev=(0)/color=blue/size=0.13/spacing=1 dic_dif



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

exit

ps_thicken coa_induced_co2uptake_shade.ps 5 > induced_co2uptake.ps
rm coa_induced_co2uptake_shade.ps




shade/ov/nokey/palette=white/nolab/levels="(-10,1,0.1)" msk

contour/nolab/ov/lev=(3.4)/color=w/size=0.18/spacing=1 o_oaragsur[l=80,d=1]