'open atoll_rain_mon_v0301.lnx.ctl' 
'set gxout shaded' 
'set grads off'  
'set grid off'  
'set x 1 48' 
'set y 1 32' 
'set z 1' 

kyr = 1991  
kmn = 1     
kchk=1 

while (kchk=1)  

  nn  = (kyr-1971)*12+kmn    

  'set vpage 0.0 8.5 0.0 11.0'   
  'set gxout grfill'  
  'set parea 1.0 7.5 6.0 10.0'   
  'set grads off'   
  'set grid off' 
  'set rgb 30 150 150 150' 
  'set rgb 31 250 250 250' 
  'set rgb 32 0 225 255' 
  'set rgb 33 0 165 185' 
  'set rgb 34 0 185 0' 
  'set rgb 35 0 255 45' 
  'set rgb 36 255 255 0' 
  'set rgb 37 225 205 0' 
  'set rgb 38 160 120 0' 
  'set rgb 39 160 40 0' 
  'set rgb 40 185 0 0' 
  'set rgb 41 255 0 155' 
  'set rgb 42 185 0 255' 
  'set rgb 43 125 0 205' 
  'set rgb 44 75 0 160' 
  'set ccols  31 32 33  34  35  36  37   38  39  40   41  42  43 44' 
*  'set ccols   32   33   34   35  36  37  38  39  40  41  42  43  44  ' 
  'set clevs  0.0  0.2  0.5   1   2   4    6   8  10  12  14 16'  

  'set t 'nn 
  'd rain'   
  'draw title atoll rainfall'   
  'run /export/lnx193/xping/GrADS/utl/cbarhov1.gs' 


  'set parea 1.0 7.5 1.0 5.0' 
  'set ccols  31 32 33 34 35 36 37 38 39 40  41  42  43 44'
  'set clevs   0  1  2  3  4  5  6  7  8  9  10 '  
  'd numb' 
  'draw title number of gauges'   
  'run /export/lnx193/xping/GrADS/utl/cbarhov1.gs' 

  'set string 1 bl 6'
  'set strsiz 0.14'
  'draw string 1.0 0.1 'kyr  
  'draw string 2.0 0.1 'kmn 

  say 'Do you want to continue ? (y/n)' 
  pull  con 
  if  (con != n)   
    kmn=kmn+1 
    if  (kmn=13)   
      kyr=kyr+1 
      kmn=1   
    endif 
    'clear' 
  else 
    kchk=0 
  endif 
endwhile