open e:/Merra2/uvt_av_wQBO.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 T(K), wind(m/s), 1.4 km, January, Climate
*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
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 10
*set xlint 45
set xlab on
set cthick 10
*set mproj nps
set gxout shaded
define tav=ave(t-ave(t,lon=0,lon=355),t=489,t=489+248)
*d ave(t-ave(t,lon=0,lon=355),t=489,t=489+248)
d tav
set strsiz 0.5
set arrscl 1 10
set ccolor 0
set lev 1
*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.9 3.7
set font 1
draw xlab Longitude (deg)
draw ylab Latitude (deg)
enable print uv.prn
print
disable print
close 1