%-----Use control run to get a location map:
cancel data/all
set mode metafile
sp rm metafile.plt*

go margins 1.5 1.5 1.5 2.0

use daily_control.nc
set region/y=50N:62N
Set region/x=5W:10E
sha/set/nolab o_temp[k=1,l=@ave]

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

go land

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

ps_thicken location_shade.ps 5 > location.ps
rm location_shade.ps

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

go margins 1.5 1.5 1.5 2.0

use control.nc
yes? set region/y=58N:59N
yes? set region/x=1E:2E

plot/set/nolab o_temp[k=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 plot

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

exit

ps_thicken sst.ps 5 > sst_thick.ps
rm sst.ps




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

go margins 1.5 1.5 1.5 2.0

use control.nc
use inf.nc
use 10.nc

yes? set region/y=58N:59N
yes? set region/x=1E:2E

plot/set/nolab/color=1 o_oaragsur[d=2,x=@ave,y=@ave]

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 plot

plot/ov o_oaragsur[d=3,x=@ave,y=@ave]

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

exit

ps_thicken omegainf.ps 5 > omegainf_thick.ps
rm omegainf.ps