open e:/Merra2/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. of T(K), 1.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 60 180
set lat 0 90
set lev 1.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(t-ave(t,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-ave(u,lon=0,lon=355),t=489,t=489+248),1,1);ave(v-ave(v,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 e:/Merra2/xy.prn
print
disable print
close 1