open  h:/lh_new/heating/jul/wave2.ctl
set lon 0
set lat -40 40
set lev 1 15
set t 1
set vpage 0.5 4.5 5.5 9.5
Set clopts 1 4 0.16
Set xlopts 1 4 0.20
Set ylopts 1 4 0.20
set cint 0.1
set cthick 10

set grads off
set vrange 1 15
set yaxis  1 15 5
set ylab on
d Aw
set annot 1 5
*set font 6
* draw title Phase (deg)
draw title LH m=2 24h westward (K/day), Jul
* draw title SPW, geopotential height, July
*draw title 2d (54h) geopotential height (m), July
draw ylab Altitude (km) 
*set font 6
set strsiz 0.22 
set string 1 c 5 90
set vpage 4. 8. 5.5 9.5
set cint 0.1
set grads off
*set yaxis 0 70 10
set ylab off
set cthick 10
set lon 1
d Aw
set annot 1 5
set font 1
draw title LH m=2 12h westward (K/day), Jul
*draw title 12-hour, m=1 nonmigrating tide
*draw title Phase (deg)
*close 1
*open  c:/s2/wave2.ctl
set lon 0
set lat -40 40
set lev 1 15
set t 1
set vpage 0.5 4.5 1.5 5.5
Set clopts 1 4 0.16
Set xlopts 1 4 0.20
Set ylopts 1 4 0.20
set cint 0.1
set grads off
set vrange 1 15
set yaxis 1 15 5
set ylab on
set cthick 10

d Ae
set annot 1 5
set font 1
draw title LH m=2 24h eastward (K/day), Jul
*draw title GH amplitude (m), SPW2, Jul2009

*draw title SDT meridional wind (m/s), January
set strsiz 0.22
set font 1
set string 1 c 5 0
draw string 4.8 0.2 Latitude (deg)
set string 1 c 5 90
draw string 0.5 4.5 Altitude (km)
set vpage 4. 8. 1.5 5.5
set cint 0.1
set grads off
*set yaxis 0 70 10
set ylab off
set cthick 10
set lon 1
d Ae
set font 1
set annot 1 5
*draw title m=1 12-hour, m=1 nonmigrating tide
draw title LH m=2 12h eastward (K/day), Jul
*draw title Phase (deg)
set font 1
set strsiz 0.22
set string 1 c 5 0
draw string 4.8 0.2 Latitude (deg)
enable print xy.prn
print
disable print
close 1