c  written by Vern Kousky 8/18/09
c  Program to compute anomaly correlations    
c23456789012345678901234567890123456789012345678901234567890123456789012
      parameter (num=930)
      dimension rcor(num),acor(num,num),idate(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-correlation-counts-jan1980-2009'
     *,form='formatted')
ccc   compute idate
      n=0
      do 99 iyr=1,30
      iyear=iyr+1979
      imon=2
      do 98 iday=1,28
      n=n+1
      idate(n)=iyear*10000+imon*100+iday
 98   continue
 99   continue
      do 1 n=1,num
      read(11,rec=n)rcor
      do 2 k=1,num
      acor(n,k)=rcor(k)
 2    continue
 1    continue
ccc   count the number of correlations greater or less than selected thresholds
      do 3 n=1,num
      fcnt1=0.
      fcnt2=0.
      fcnt3=0.
      fcnt4=0.
      fcnt5=0.
      fcnt6=0.
      fcnt7=0.
      fcnt8=0.
      do 4 k=1,num
      if(acor(n,k).le.-70.)fcnt1=fcnt1+1.
      if(acor(n,k).le.-60.)fcnt2=fcnt2+1.
      if(acor(n,k).le.-50.)fcnt3=fcnt3+1.
      if(acor(n,k).le.-40.)fcnt4=fcnt4+1.
      if(acor(n,k).ge.40.)fcnt5=fcnt5+1.
      if(acor(n,k).ge.50.)fcnt6=fcnt6+1.
      if(acor(n,k).ge.60.)fcnt7=fcnt7+1.
      if(acor(n,k).ge.70.)fcnt8=fcnt8+1.
 4    continue
      write(51,200)idate(n),fcnt1,fcnt2,fcnt3,fcnt4,fcnt5,fcnt6,fcnt7,fc
     *nt8
 200  format(i10,8f6.0)
 3    continue
      stop
      end