program datacut
c
      integer i,j,im,jm,irec,irec2
      parameter (im=176, jm=137)
      real var(im,jm)
c----------------- atom -----------------------------------------
c  uu:1-23 vv:24-46 tt:47:69 qv:70-92 pt:93-115 pp:116-138 pw:139
c----------------- surf -----------------------------------------
c 1.www 2.rain 3.evap 4.runoff 5.SWd 6.SWu 7.LWd 8.LWu 
c 9.Rnet 10.Qle 11.Qh 12.tau 13.Ts 14.Td 15.W1 16.W2
c 17.W3 18.Tm 19.em 20.NDVI 
c-----------------------------------------------------
c
c  0621 00UTC -  0721 00UTC  (hourly)
c
c----------------------------- 
      open(11,file='HUBEXatom0621-0721v15.dat'
     * ,form='unformatted',access='direct',recl=im*jm)
      open(12,file='HUBEXsurf0621-0721v15.dat'
     * ,form='unformatted',access='direct',recl=im*jm)
      open(21,file='HUBEXatom070100-070206.dat'
     * ,form='unformatted',access='direct',recl=im*jm)
      open(22,file='HUBEXsurf070100-070206.dat'
     * ,form='unformatted',access='direct',recl=im*jm)
c
c====== atmosphere ======
      irec=139*24*10
      irec2=0
c
      do ihour=1,24  ! 0701 00UTC - 0701 23UTC
      do ivar=1,139
      irec=irec+1
      irec2=irec2+1
      read(11,rec=irec) ((var(i,j),i=1,im),j=1,jm) 
      write(21,rec=irec2) ((var(i,j),i=1,im),j=1,jm) 
      enddo
      enddo
c
      do ihour=1,7   ! 0702 00UTC - 0702 06UTC
      do ivar=1,139
      irec=irec+1
      irec2=irec2+1
      read(11,rec=irec) ((var(i,j),i=1,im),j=1,jm) 
      write(21,rec=irec2) ((var(i,j),i=1,im),j=1,jm) 
      enddo
      enddo

c====== surface ======
      irec=20*24*10
      irec2=0

      do ihour=1,24  ! 0701 00UTC - 0701 23UTC
      do ivar=1,20
      irec=irec+1
      irec2=irec2+1
      read(12,rec=irec) ((var(i,j),i=1,im),j=1,jm) 
      write(22,rec=irec2) ((var(i,j),i=1,im),j=1,jm) 
      enddo
      enddo
c
      do ihour=1,7   ! 0702 00UTC - 0702 06UTC
      do ivar=1,20
      irec=irec+1
      irec2=irec2+1
      read(12,rec=irec) ((var(i,j),i=1,im),j=1,jm) 
      write(22,rec=irec2) ((var(i,j),i=1,im),j=1,jm) 
      enddo
      enddo
c
      close(11)
      close(21)
      close(12)
      close(22)
c
      stop
      end
c