open e:/Merra2/uvt_av_wQBO.ctl
open e:/Merra2/uvt_av_eQBO.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 wQBO-eQBO T(K) and wind(m/s), 4 km, Jan
*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 30 150
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 xaxis 30 150 30
set clab forced
set ylint 15
set xlint 30
set cterp on
set xlab on
set cthick 10
*set mproj nps
set gxout shaded
define tw=ave(t-ave(t,lon=0,lon=355),t=489,t=489+248)
define te=ave(t.2-ave(t.2,lon=0,lon=355),t=489,t=489+248)
*d ave(t-ave(t,lon=0,lon=355),t=489,t=489+248)
d tw-te
set strsiz 0.5
set arrscl 1 20
set ccolor 0
set lev 4
*d skip(u,1,1);v
define uw=ave(u-ave(u,lon=0,lon=355),t=489,t=489+248)
define vw=ave(v-ave(v,lon=0,lon=355),t=489,t=489+248)
define ue=ave(u.2-ave(u.2,lon=0,lon=355),t=489,t=489+248)
define ve=ave(v.2-ave(v.2,lon=0,lon=355),t=489,t=489+248)
define du=uw-ue
define dv=vw-ve
*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)
d skip(du,1,1);dv
cbarn 1 1 7.9 3.7
set font 1
draw xlab Longitude (deg)
draw ylab Latitude (deg)
enable print g:/Merra2/uv.prn
print
disable print
close 2
close 1