c*********************************************************                      
c* fortran program to read and write wcgrid.asc (file 5) *                      
c*********************************************************                      
      integer nlin                                                              
      integer id, elnum, gl, gm, slyr, ernum, trnum                             
      real    elavg, elmax, elmin, slr, slg, slc, sls                           
      real    eravg, ermax, ermin, travg, trmax, trlvl                          
      real    whavg, whmax, whsd                                                
c*********************************************************                      
c* initialize a counter and open files for input/output  *                      
c*********************************************************                      
      nlin=0                                                                    
      open(unit=5,file='wcgrid.asc',readonly,status='old')                      
      open(unit=6,file='wcgrid.out',status='new')                               
c*********************************************************                      
c* read/write the grid cell id and the 22 data variables *                      
c*********************************************************                      
   10 read(5,100,end=999) id, elavg, elmax, elmin, elnum,                       
     1    gl, gm, slr, slg, slc, sls, slyr                                      
      read(5,110) eravg, ermax, ermin, ernum, travg,                            
     1    trmax, trlvl, trnum, whavg, whmax, whsd                               
                                                                                
      if (nlin.gt.32) nlin=0                                                    
      if (nlin.eq.0) write (6,120)                                              
      if (nlin.eq.0) write (6,130)                                              
      nlin=nlin+1                                                               
                                                                                
      write(6,105) id, elavg, elmax, elmin, elnum,                              
     1     gl, gm, slr, slg, slc, sls, slyr                                     
      write(6,115) eravg, ermax, ermin, ernum, travg,                           
     1     trmax, trlvl, trnum, whavg, whmax, whsd                              
   20 continue                                                                  
      go to 10                                                                  
  100 format (i5,3f8.2,3i4,4f8.2,i4)                                            
  105 format (1x,i5,3f8.2,3i4,4f8.2,i4)                                         
  110 format (3f8.2,i4,3f8.2,i4,3f8.2)                                          
  115 format (1x,3f8.2,i4,3f8.2,i4,3f8.2)                                       
  120 format (4x,'id',1x,'elavg',3x,'elmax',3x,'elmin',1x,                      
     1    'elnum',2x,'gl',2x,'gm',1x,'slr',5x,'slg',5x,'slc',                   
     1    5x,'sls',4x,'slyr')                                                   
  130 format (2x,'eravg',3x,'ermax',3x,'ermin',1x,'ernum',                      
     1    1x,'travg',3x,'trmax',3x,'trlvl',1x,'trnum',1x,                       
     1    'whavg',3x,'whmax',3x,'whsd')                                         
c***********************************************************                    
c*****           close files and exit                 ******                    
c***********************************************************                    
  999 close(unit=5)                                                             
      close(unit=6)                                                             
      stop                                                                      
      end