'open /cpc/home/wd52vk/samer/satellite/gridsat.ctl'
'set vpage 0.0 8.5 0.0 11.0'
'set parea 0.4 7.3 3.5 10.4'
'set annot 1 5'
'set cterp off'
'set csmooth on'
'set lat -50 10'
'set lon -100 -20'
'set time 00Z12jan2005'
'q dims'
temp=sublin(result,5)
beg=subwrd(temp,9)
end = beg+63
it=beg
n=1
while (it <= end)
'set t 'it
cdate = mydate()
chr = subwrd(cdate,1)
cdy = subwrd(cdate,2)
cmo = subwrd(cdate,3)
cyr = subwrd(cdate,4)
'set grads off'
'set mpdset hires'
'set map 3 1 8'
'set string 1 tc 7'
'set strsiz 0.15 0.15'
'draw string 4.25 10.9 IR' 
'draw string 4.25 10.6 'chr''cdy''cmo''cyr
'run /cpc/home/wd52vk/grads/util/bw2.gs'
'set clevs 180 190 200 210 220 230 240 250 260 270 280 290 300'
'set ccols 55 53 51 49 47 45 43 41 39 37 35 33 31 30'
'set gxout shaded'
'd irwin'
'run /cpc/home/wd52vk/grads/util/cbarhovn.gs'
'printim IR_'n'.gif white'
'disable print'
'clear'
n=n+1
it=it+1
endwhile
function mydate
'query time'
sres = subwrd(result,3)
i=1
while (substr(sres,i,1)!='Z')
i=i+1
endwhile
hour = substr(sres,1,i)
isav = i
i = i + 1
while (substr(sres,i,1)>='0' & substr(sres,i,1)<='9')
i = i + 1
endwhile
day = substr(sres,isav+1,i-isav-1)
month = substr(sres,i,3)
year = substr(sres,i+3,4)
return (hour' 'day' 'month' 'year)