c                                                                      
c     this program reads the daily gridded data for real-time OMI
c     data for January 31, 2007.  variables lat and lon contain the 
c     latitudes and longitudes of the center of each of the grid    
c     cells in the array ozone.                                    
c                                                                 
      character*80   header                                      
      integer*4      ozone                                      
      real*4         lat(180), lon(360)                        
                                                              
      dimension      ozone(360,180)                          
c                                                           
c     calculate latitudes and longitudes                   
c                                                         
      dlat = 1.0                                         
      do 10 i=1,180                                     
      lat(i) = -89.5 + (i-1)*dlat                      
10    continue                                        
      dlon = 1.0                                    
      do 20 i=1,360                                 
      lon(i) = -179.5 + (i-1)*dlon               
20    continue                                    
c                                                
c     open the input file                       
c                                              
      open(1,file='L3_ozone_omi_20070131.txt',status='old')               
c                                                           
c     read in the header lines                             
c                                                         
      read(1,'(a80)') header                             
      read(1,'(a80)') header                            
      read(1,'(a80)') header                           
c                                                     
c     read in the data into the array ozone          
c                                                   
      do 30 i=1,180                                
      read(1,'(1x,25i3)') (ozone(j,i),j=1,360)    
30    continue                                                                
c                                                             
c     close the input file                                   
c                                                           
      close(1)                                             
c                                                         
c     process/print the ozone data                       
c                                                       
      stop                                             
      end