my subplot

*
 itot=3
 jtot=3
 dx=11.0/itot
 dy=8.1/jtot
 dxovl=0.15
 dyovl=0.15
 dxd=dx-dxovl
 dyd=dy-dyovl
 xsi=0.2
 ysi=0.2

i=0
  jx=jtot
**********
 while(jx>=1)
  ix=1
**********
 while(ix<=itot)
 i=i+1
*
 xs=xsi+(ix-1)*dxd
 xe=xs+dx
 ys=ysi+(jx-1)*dyd
 ye=ys+dy
 'set vpage 'xs' 'xe' 'ys' 'ye
*
 say ix ' 'jx' 'i
 say 'set vpage 'xs' 'xe' 'ys' 'ye
*

'set t 'i
pull dummy
*
'set grads off'
'set gxout grfill'
'set clevs 0 0.05 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 '
'set ccols 0 41 42 44 46 33 35 37 23 25 27 29'
'define a'i'=ave(t(z='i'),time=jan1860,time=dec1869)'
*'define b'i'=ave(t(z='i')',time=jan1990,time=dec1999)'
'd a'i
'cbarn'
ix=ix+1
endwhile

jx=jx-1
endwhile