Praj-Paramit

'define MOclm=ave('varm'*'unitm',t+'tt1',t='tt2',12)'
'modify MOclm seasonal'
*
'define MOpow=ave(pow('varm'*'unitm'-MOclm,2),t+'tt1',t='tt2',12)'
'define MOstd=sqrt(MOpow)'
'modify MOpow seasonal'
'modify MOstd seasonal'

while  ( i <= 20 )

years=1900+(i-1)*5
yeare=1900+i*5-1

is=(i-1)*12*5+1
ie=i*12*5
say is ' 'ie
*---------------------------
  if(i<=25) ; ii=5 ; jj=i-5*4 ; endif
  if(i<=20) ; ii=4 ; jj=i-5*3; endif
  if(i<=15) ; ii=3 ; jj=i-5*2 ; endif
  if(i<=10) ; ii=2 ; jj=i-5*1 ; endif
  if(i<=5)  ; ii=1 ; jj=i ; endif
  if(jj=4) ; jj=6 ; endif
  if(jj=5) ; jj=7 ; endif

*
  x0=8.6
  y0=5.0
  x1 = x0
  x2 = x0+0.5
  x3 = x0+0.8
  y1=y0-i*0.2+0.05
  y2=y0-i*0.2
*
'set ccolor 'ii
'set cstyle 'jj
'set cmark 0'

'set line ' ii ' 'jj
'draw line 'x1' ' y1' ' x2' ' y1
'draw string 'x3 ' 'y2' ' years'-'yeare
*
i=i+1
endwhile

*
 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