open d:/LH_ENSO/El/s10.ctl
set lat -40 40
set lev 0 15
set t 1
set vpage 0.5 4.5 7.5 10.5
*Set clopts 1 4 0.16
*Set xlopts 1 4 0.20
*Set ylopts 1 4 0.20
set cint 0.2
set cthick 10
set grads off
set vrange 0 15
set yaxis  0 15 5

set ylab on
d ave(lh, t=1, t=248)
set annot 1 5
set font 1
set string 1 c 6 0
set strsiz 0.2
draw string  4.7 6.2 a) Zonally averaged heating, K/day, EN
set strsiz 0.22
set font 1
set string 1 c 5 0
*draw string 4.8 0.3 Latitude (deg)
set string 1 c 5 90
draw string 0.7 3.5 log-pressure height (km)
set string 1 c 6 0
set strsiz 0.3
draw string  2.0 5.5 
close 1

open d:/LH_ENSO/El/wave1.ctl

set vpage 0.5 4.5 4.5 7.5
set lon 2
set lat -40 40
set lev 1 15
set t 1

set cint 0.2
set cthick 10
set grads off
set yaxis 0 15 5
set ylab on
d Aw
set annot 1 5
set font 1
set string 1 c 6 0
set strsiz 0.2
draw string  4.7 6.2 c) SPW1 forcing, K/day, El-Nino
set strsiz 0.22
set font 1
set string 1 c 5 0
*draw string 4.8 0.3 Latitude (deg)
set string 1 c 5 90
draw string 0.7 3.5 log-pressure height (km)
set string 1 c 6 0
set strsiz 0.3
draw string  2.0 5.5 

set font 1
*draw xlab Latitude (deg)
*draw ylab log-pressure height (km)
close 1

open  d:/LH_ENSO/El/wave2.ctl
set lon 2
set lat -40 40
set lev 1 15
set t 1
set vpage 0.5 4.5 1.5 4.5

*set font 1
*draw xlab Latitude (deg)
*draw ylab log-pressure height (km)
set font 1
set strsiz 0.22 
set string 1 c 5 90

draw string 0.7 3.5 log-pressure height (km)

set cint 0.2
set cthick 10
set grads off
set yaxis 0 15 5
set ylab on
d Aw
set annot 1 5
set font 1
set string 1 c 6 0
set strsiz 0.2
draw string  4.7 6.2 e) SPW2 forcing, K/day, El-Nino
set strsiz 0.22
set font 1
set string 1 c 5 0
draw string 4.8 0.3 Latitude (deg)
set font 1
set strsiz 0.22 
set string 1 c 5 90

draw string 0.7 3.5 log-pressure height (km)
set string 1 c 6 0
set strsiz 0.3
draw string  2.0 5.5 

*set font 1
*draw xlab Latitude (deg)
*draw ylab log-pressure height (km)
close 1

open d:/LH_ENSO/La/s10.ctl

set lat -40 40
set lev 0 15
set t 1
set vpage 4.0 8.0 7.5 10.5

Set clopts 1 4 0.16
*Set xlopts 1 4 0.20
Set ylopts 1 4 0.20
set cint 0.2
set cthick 10
set grads off
set vrange 0 15
set yaxis  0 15 5
*set xaxis  30N 90N 10N
*set xlab on
set ylab on
set clab forced
d ave(lh, t=1, t=248)
set annot 1 5
set font 1
set string 1 c 6 0
set strsiz 0.2
draw string  4.7 6.2 b) Zonally averaged heating, K/day, LN

set strsiz 0.22
set font 1
set string 1 c 5 0
*draw string 4.8 0.3 Latitude (deg)
set string 1 c 5 90
*draw string 0.7 4.5 log-pressure height (km)
set string 1 c 6 0
set strsiz 0.3
draw string  2.0 5.5 
close 1

open d:/LH_ENSO/La/wave1.ctl
set lon 2
set lat -40 40
set lev 1 15
set t 1

set vpage 4.0 8.0 4.5 7.5
set font 1
*draw xlab Latitude (deg)
*draw ylab log-pressure height (km)
*set font 1
*set strsiz 0.22 
*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.2
set cthick 10
set grads off
set yaxis 0 15 5
set ylab on
set clab forced
d Aw
set annot 1 5
set font 1
set string 1 c 6 0
set strsiz 0.2
draw string  4.7 6.2 d) SPW1 forcing, K/day, La-Nina
set strsiz 0.22
set font 1
set string 1 c 5 0
*draw string 4.8 0.3 Latitude (deg)
set font 1
*draw xlab Latitude (deg)
*draw ylab log-pressure height (km)
set string 1 c 6 0
set strsiz 0.3
draw string  2.0 5.5 
close 1

open  d:/LH_ENSO/La/wave2.ctl

set lon 2
set lat -40 40
set lev 1 15
set t 1

set vpage 4.0 8.0 1.5 4.5
Set clopts 1 4 0.16
Set xlopts 1 4 0.20
Set ylopts 1 4 0.20
set cint 0.2
set cthick 10
set grads off
set vrange 0 15
set yaxis  0 15 5
*set xaxis  30N 90N 10N
set xlab on
set ylab on
d Aw
set annot 1 5
set font 1
set string 1 c 6 0
set strsiz 0.2
draw string  4.7 6.2 f) SPW2 forcing, K/day, La-Nina
set strsiz 0.22
set font 1
set string 1 c 5 0
draw string 4.8 0.3 Latitude (deg)
set string 1 c 5 90
*draw string 0.7 4.5 log-pressure height (km)
set string 1 c 6 0
set strsiz 0.3
draw string  2.0 5.5 

enable print Fig2.prn 
print
*print d:/LH_ENSO/La/Fig2.ps
disable print

close 1