program rdtest

      real :: amp(72,59,55)
      real :: glat(59),ht(55),glon(72),temp(6)
      integer :: i,j,k,m
C  #23-30=90.57 to 120.019km so save these 8 hts
      data ht/0.,4.161,8.188,12.201,16.218,20.239,24.264,28.295,32.332,
     |  36.378,40.434,44.503,48.588,52.691,56.819,60.975,65.167,
     |  69.403,73.695,78.056,82.129,86.293,90.57,94.58,98.736,
     |  103.075,107.177,111.514,115.619,120.019,124.175,128.688,
     |  132.911,137.539,141.772,146.414,150.470,154.871,159.653,
     |  164.847,169.021,173.437,178.086,182.957,188.030,193.283,
     |  198.689,204.224,209.865,215.589,221.379,227.221,233.102,
     |  239.013,244.946/
      

      do i=1,72
        glon(i) = (i-1)*5
      enddo
      do j=1,59
        glat(j) = -90 + j*3
      enddo

      open (10,file='temper_sep_all_Dam.dat',status='old')

      do k=1,55
      do j=1,59
       do i=1,72,6
         read (10,*) temp
          do m=1,6
           amp(i+m-1,j,k) = temp(m)
          enddo
       enddo
      enddo
      enddo
    
      do j=1,59
        write (6,"(1x,'j glat amp(3,j,3,31,55)=',i3,f6.1,3f11.6)")
     |   j,glat(j),amp(3,j,3),amp(3,j,31),amp(3,j,55)
      enddo

      stop
      end