************************************************************************
***** 標高データを配列に代入するプログラム(import) *********************
c      subroutine(a,e,n)
      real    a(4800,4800),k(1600,1600)
      integer i, j, n,x,y
      character mn*4
************************************************************************
      do n=1,9
      write(*,*) '地域番号を入れてください。'
      read(*,*)mn
      open(unit=77,file=mn//'.bin',access='direct',
     +                 form='unformatted',  recl=6400)
       do i=1,1600
       read(77,rec=i)(k(i,j),j=1,1600)   
       end do
      close(77)
************************************************************************


       if(n.ge.1.and.n.le.3)y=0
       if(n.ge.4.and.n.le.6)y=1
       if(n.ge.7.and.n.le.9)y=2
       if(n.eq.1 .or. n.eq.4 .or. n.eq.7)x=0
       if(n.eq.2 .or. n.eq.5 .or. n.eq.8)x=1
       if(n.eq.3 .or. n.eq.6 .or. n.eq.9)x=2
       do i=1,1600
        do j=1,1600
        a(i+1600*y,j+1600*x)=k(i,j)
        end do
       end do
      end do
      
      open(unit=7,file='data.terrain.bin',access='direct'
     +                  ,form='unformatted',recl=19200)
       do n=1,4800
         write(7,rec=n)(a(n,m),m=1,4800)
       end do
      close(unit=7)
      
      
c      write(*,*)a(3000,4000)
************************************************************************
c      close(78)
************************************************************************
      end
