c RDATA.F 
c 
c PURPOSE:  TO CALCULATE ARCTIC MODEL GRID COORDINATES, AND READ 
c           IN LAND MASK 

        subroutine rdata(xt,yt,zt,kmt,imt,jmt,km) 
c 
        real xt(imt),yt(jmt),zt(km) 
        integer kmt(imt,jmt) 
c 
c assign ARCTIC MODEL data 
c xt,yt,zt 
        xt(1) = -21. 
        yt(1) = -14. 
        do i=2,imt 
                xt(i) = xt(i-1) + 0.5 
        enddo 
        do j=2,jmt 
                yt(j) = yt(j-1) + 0.5 
        enddo 
        do k=1,km 
                zt(k) = float(k)*10.   ! dummy 
        enddo 
c 
c kmt, setting max depths to 12 
        open(31,file='kmt.arc',status='old',iostat=ierr) 
          do j=jmt,1,-1 
            do i=1,imt 
              read(31,'(1x,i3)') kmt(i,j) 
              if(kmt(i,j).gt.0) kmt(i,j)=km 
            enddo 
          enddo 
        close(31) 
c 
        return 
        end