parameter(it=162,jt=82,kt=130)
      dimension xc(it,jt,kt),yc(it,jt,kt),zc(it,jt,kt)

      open(unit=18,file='mesh_3dhill.dat',form='formatted')
      read(18,*) ni,nj,nk
      write(6,*)'ni,nj,nk',ni,nj,nk
      nim1=ni-1
      njm1=nj-1
      nkm1=nk-1
      do k=1,nkm1
         do j=1,njm1
            do i=1,nim1
               read(18,*)xc(i,j,k),yc(i,j,k),zc(i,j,k)
           end do
        end do
      end do
      close(18)


      stop
      end