character*40 tline
          character*14 inom(6)
	    integer nlin(6)
	    real xval(12)
	    real xsoil(6,8355,12)  !day taken, layer
	    integer zcnt(6,60)
	    integer dayrec(6,270)
	    integer lomo(23,20),lody(23,20),lopnt(23)
	    character*1 xmark(23,240)
	    
	    do i=1,23
	     do j=1,240
	      xmark(i,j)=' '
	     enddo
	    enddo
	    
	    data inom/'Iowa.W2-1.dat','Iowa.W2-3.dat',
     &  'Iowa.W2-4.dat','Iowa.W3-11.dat','Iowa.W3-14.dat',
     &  'Iowa.W3-16.dat'/
          data nlin/270,270,270,269,269,269/

15        format(a40)
20	    format(i4,x,i2,x,i2,x,i3,x,i5,6x,i2,
     &        12(x,f5.1))	         
     
          undef=-99.99
	    do i=1,6
           do j=1,8355
	      do k=1,12
	       xsoil(i,j,k)=undef
	      enddo
	     enddo
          enddo
     
          do i=1,6
	     jdmax=0
	     izpnt=0
	     open(10,file=inom(i),status='old')
	     do j=1,9
  	       read(10,15) tline
	     enddo
	     
	     do j=1,nlin(i)
	       read(10,20) iyr,imn,idy,idoy,idtt,isit,
     &            (xval(kk),kk=1,12)
              if(i.eq.1) then
		    ipt=iyr-1971
		    jpt=idoy-90
		    xmark(ipt,jpt)='X'
		    lopnt(ipt)=lopnt(ipt)+1
		    lomo(ipt,lopnt(ipt))=imn
		    lody(ipt,lopnt(ipt))=idy
		  endif
              idif=idtt-izpnt
		  if(idif.lt.61) zcnt(i,idif)=zcnt(i,idif)+1
		  izpnt=idtt
		  dayrec(i,j)=idtt
             do k=1,12
 		   xsoil(i,idtt,k)=xval(k)
		 enddo
		 if(idtt.gt.jdmax) jdmax=idtt
	     enddo  !the j loop
	     print*,'maximu day is',jdmax
	    enddo  !the i loop		 
	
	
45        format(i3,x,6i6)	
	    do i=1,60
c	      write(*,45) i,(zcnt(jj,i),jj=1,6)
	    enddo
	    write(*,*)
	    do j=1,270
c	     write(*,45) j,(dayrec(ii,j),ii=1,6)
	    enddo
c	    open(30,file='iowa.bdat',status='unknown',
c     &      access='direct',recl=1)
          ire=0
	    
c 	    do j=1,8355
 	     do i=1,6
	      do k=1,12
	         ire=ire+1
c 		   write(30,rec=ire) xsoil(i,j,k)
 		enddo
 	     enddo
c	    enddo
65        format(i4,x,14(i3,'/',i2,2x))
          do i=1,23
	     iyr=i+1971
c	     write(*,*) 'o pnt is',lopnt(i)
 	     write(*,65) iyr,(lomo(i,jj),lody(i,jj),jj=1,lopnt(i))
	    enddo
75	    format(240a1)
	    do i=1,23
c	     write(*,75) (xmark(i,jj),jj=1,240)
	    enddo
	    stop
	    end