c  written by Vern Kousky 8/18/09
c  Program to compute anomaly correlations    
c23456789012345678901234567890123456789012345678901234567890123456789012
      parameter (num=930)
      dimension rcor(num),idate(num)
      character path*33
      data path/'/cpc/home/wd52vk/samer/satellite/'/
      open (11,file='pattern-correlations-comp-sebr-dry-jan1980-2009'
     *,access='direct',form='unformatted',recl=num*4)
      open (51,file='pattern-correlation-comp-dates-sebr-dry'
     *,form='formatted')
ccc   compute idate
      n=0
      do 99 iyr=1,30
      iyear=iyr+1979
      imon=1
      do 98 iday=1,31
      n=n+1
      idate(n)=iyear*10000+imon*100+iday
 98   continue
 99   continue
      read(11,rec=1)rcor
ccc   count the number of correlations greater or less than selected thresholds
      fcnt1=0.
      fcnt2=0.
      fcnt3=0.
      fcnt4=0.
      fcnt5=0.
      fcnt6=0.
      fcnt7=0.
      fcnt8=0.
      do 4 k=1,num
      if(rcor(k).le.-70.)fcnt1=fcnt1+1.
      if(rcor(k).le.-60.)fcnt2=fcnt2+1.
      if(rcor(k).le.-50.)fcnt3=fcnt3+1.
      if(rcor(k).le.-40.)fcnt4=fcnt4+1.
      if(rcor(k).ge.40.)fcnt5=fcnt5+1.
      if(rcor(k).ge.50.)fcnt6=fcnt6+1.
      if(rcor(k).ge.60.)fcnt7=fcnt7+1.
      if(rcor(k).ge.70.)fcnt8=fcnt8+1.
      if(rcor(k).ge.70.)write(51,201)idate(k),rcor(k)
 201  format(i10,2x,f5.1)
 4    continue
      write(6,200)fcnt1,fcnt2,fcnt3,fcnt4,fcnt5,fcnt6,fcnt7,fcnt8
 200  format(8f6.0)
 3    continue
      stop
      end