c
c  This program is to read Mongolian soil moisture data.
c       Created by Mingquan, Mu in April 19, 2001. 
c
         parameter (num=42)
         character*12 name(num),agri(num)
         character*25 info(num)
         integer id(num),ibeg(num),iend(num)
         real lat(num),lon(num),alt(num)
	 real sst(num,11,21,1964:1993) !station,layer,obs,year

         undef=-9.9
         
         do iy=1964,1993
          do k=1,21
           do li=1,11
            do i=1,num
		 sst(i,li,k,iy)=undef
            enddo
           enddo
          enddo
         enddo
 
C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
         open(11,file='soilm_mong.dat',status='old')
        
         do i=1,9
           read(11,*)
         enddo
         do 101 i=1,num
           read(11,21) name(i),id(i),ibeg(i),iend(i),
     &             lat(i),lon(i),alt(i),agri(i),info(i)
           write(*,21) name(i),id(i),ibeg(i),iend(i),
     &             lat(i),lon(i),alt(i),agri(i),info(i)
101      continue

         do 201 i=1,num
          print*,'i is',i

           do ii=1,4
              read(11,*)
           enddo

           do 201 iy=ibeg(i),iend(i)
             read(11,*) 
             read(11,*) 
             do 201 li=1,11
               read(11,31) (sst(i,li,k,iy),k=1,21) 
               write(*,31) (sst(i,li,k,iy),k=1,21) 
 201     continue
         close(11)

21       format(4X,A12,4X,I5,3X,I4,1X,I4,2X,F5.2,3X,
     +        F6.2,3X,F6.1,2X,A7,3X,A25)
31       format(5X,21(x,f5.1))

         print*,'done '
         stop
         end