c  written by Vern Kousky 8/18/09
c23456789012345678901234567890123456789012345678901234567890123456789012
      parameter (num=930)
      dimension rcor(num)
      character path*33
      data path/'/cpc/home/wd52vk/samer/satellite/'/
      open (11,file='pattern-correlations-comp-sebr-wet-new-jan1980-2009
     *',access='direct',form='unformatted',recl=num*4)
      open (51,file='pattern-correlations-comp-rcor60-6jan1986'
     *,form='formatted')
      read(11,rec=1)rcor
      do 1 k=1,num
      if(rcor(k).lt.60.)go to 1
ccc   compute date
      krec=0
      do 13 iyr=1,30
      nyear=iyr+1979
      imon=1
      do 14 iday=1,31
      krec=krec+1
      kdate=nyear*10000+imon*100+iday
      if(krec.eq.k)go to 15
 14   continue
 13   continue
 15   write(6,206)kdate,krec,rcor(k)
 206  format(2i10,f8.1)
      write(51,206)kdate,krec,rcor(k)
 1    continue
      stop
      end