open e:/Merra2/Europe/phi_av_CLimate.ctl
open e:/Merra2/Europe/uvt_av_CLimate.ctl
reset
set t 1
set font 1
set string 1 c 6 0
set strsiz 0.16
*draw string  4.4 8.1 Wave activity flux during SSW
*draw string  4.4 8.1 Wave activity flux 15d before SSW
draw string  4.4 7.7 Dev. h(m) and wind(m/s), 4 km, Jan, Clim
*draw string  4.4 8.1 Fz, 20km, 10 Feb (after SSW)
*draw string  4.4 8.1 Fz, 20km, during SSW
*draw string  4.4 8.1 Fz, 20km, after SSW

set lon 00 60
set lat 0 90
set lev 4
set clopts 1 6 0.2
set xlopts 1 4 0.14
set ylopts 1 4 0.14
*set arrscl 0.5
*set arrscl 0.5
set arrowhead 0.07
set vpage 0.5 8. 1.5 8
set grid on
set grads off
set clab forced
set ylint 15
set xlint 30
set cterp on
set xlab on
set cthick 10
*set mproj nps
set gxout shaded
d ave(h-ave(h,lon=0,lon=355),t=489,t=489+248)
*d t
set strsiz 0.5
set arrscl 1 20
set ccolor 0
set lev 4
*d skip(u,1,1);v
d skip(ave(u.2-ave(u.2,lon=0,lon=355),t=489,t=489+248),1,1);ave(v.2-ave(v.2,lon=0,lon=355),t=489,t=489+248)
cbarn 1 1 7.8 3.7
set font 1
draw xlab Longitude (deg)
draw ylab Latitude (deg)
enable print xy.prn
print 
disable print
close 2
close 1