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-irwin-SA-jan1980-2009' *,access='direct',form='unformatted',recl=num*4) open (51,file='pattern-correlations-rcor60-17jan2006' *,form='formatted') nrec=0 idate=20060117 ccc compute record number - n do 10 iyr=1,30 iyear=iyr+1979 imon=1 do 11 iday=1,31 nrec=nrec+1 ndate=iyear*10000+imon*100+iday if(ndate.eq.idate)go to 12 11 continue 10 continue 12 n=nrec write(51,210)idate,nrec write(6,210)idate,nrec 210 format(2i10) read(11,rec=n)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