program readaari

      parameter (width = 721, height = 721, scale = 0.01)

      real*4 lat(width,height), lon(width,height)
      integer*2 datum, column, row

      open(unit=1, file='aarilats.bin', form='unformatted',
     1     access='direct', recl=2)

      open(unit=2, file='aarilons.bin', form='unformatted',
     1     access='direct', recl=2)

C     note switch to column,row here

      do 200 row = 1, height
        do 100 column = 1, width
          read(1) datum
          lat(column,row) = datum*scale
          read(2) datum
          lon(column,row) = datum*scale
 100    continue
 200  continue

      close(1)
      close(2)


      write(*,*) 'column, row, lat, lon'
      write(*,*) 721, 721, lat(721,721), lon(721,721)
      write(*,*) 1, 1, lat(1,1), lon(1,1)
      write(*,*) 1, 721, lat(1,721), lon(1,721)
      write(*,*) 721, 1, lat(721,1), lon(721,1)
      write(*,*) 361, 361, lat(361,361), lon(361,361)

      end