reinit
sdfopen /r1a/ncep/surface/monthly-climate/nc/air.mon.ltm.nc
sdfopen  /r1a/ncep/surface/daily/nc/air.sig995.2003.nc 
set lon -20 160
set lat 25 90
set grads off
set mproj nps

set gxout shaded
set clevs -9 -3 3 9
set ccols 9 4 0 8 2
d air.2(time=9jul2003)-air(t=7)-273

close 2
close 1
sdfopen /r1a/ncep/plev/monthly-climate/nc/hgt.mon.ltm.nc
sdfopen  /r1a/ncep/plev/daily/nc/hgt.2003.nc 

set grads off
set lon -20 160
set lat 25 90
set mproj nps
set lev 300
set gxout contour
*set cint 100
set clevs -450 -350 -250 -150 -50 50 150 250 350 450
set cthick 5
set clopts 1 1 0.1
d hgt.2(time=9jul2003)-hgt.1(t=7)

draw title za300(cntr) & tsa(shd) (2003-climate) 9jul2003

enable print za300-tsa.9jul.ncep.gm
print
disable print


!gxeps -2 -c -i za300-tsa.9jul.ncep.gm -o za300-tsa.9jul.ncep.eps
!convert -rotate 90 -geometry 600x400  za300-tsa.9jul.ncep.eps za300-tsa.9jul.ncep.gif