open d:/Muam_NF_wQBO/1_3/m1_gh.ctl
open d:/Muam_NF_eQBO/1_3/m1_gh.ctl
reset
set t 1
set font 1
set string 1 c 6 0
set strsiz 0.22
draw string  4.7 8. Zonal wind, m/s, January
set lon 0
set lat -90 90
set lev 0 100
set clopts 1 5 0.16
set xlopts 1 5 0.16
set ylopts 1 5 0.16
set vpage 0.5 8. 1. 8.
set vrange 0 100
set yaxis  0 100 10
set grid off
set grads off
set cthick 10
set cint 50
*d ave(u,t=1,t=1201)-ave(u.2,t=1,t=1201)
*d ave(a.2,t=361,t=1069)
d ave(a,t=361,t=1069)-ave(a.2,t=361,t=1069)
set font 1
draw xlab Latitude (deg)
draw ylab Geopotential Height (km)
enable print uv.prn
print
disable print
close 2
close 1