#! /bin/sh
#
#    Last change: Norel Rimbu 2004
#
if [ $# -lt 4 ]
then
  echo " "
  echo "           *****     CORWINDOW     *****      "
  echo " "
  echo "      The program performs window correlations"
  echo " "
  echo " Syntax: corwindow <file1> <file2> <file3> <win_length>"
  echo " "
  echo "    <file1>      - input series 1 (one record)"
  echo "    <file2>      - input series 2 (one record)"
  echo "    <file3>      - input field"
  echo "    <win_length> - window length"
  echo " "
  exit 0
else
 echo " "
 echo "     *****     CORWINDOW     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 echo "file3="$3
 echo "Window length="$4
 cp $1 fort.1
 cp $2 fort.2
fi
cat > corwindow.f <<EOF
C
C       Correlates two series in one GIVEN WINDOW file 
C
      program corwindow
      dimension x(20000),y(20000),r(20000)
      integer ia(4),ib(4)

      read(*,'(i10)')nw

C  I read files of an input

      read(1)(ia(l),l=1,4)
      n=ia(4)
      read(1)(x(l),l=1,n)
      read(2)(ib(l),l=1,4)
      n1=ib(4)
      read(2)(y(l),l=1,n1)

C  Has stopped to read files of an input 

C  The algorithm follows

             
C I calculate average 


      nf=n-nw
C      write(*,*)nf
C      read(*,*)

      do 99 k=1,nf
        s1=0.0
        s2=0.0
        nnw=k+nw-1
        do 1 l=k,nnw
          s1=s1+x(l)
          s2=s2+y(l)
    1   continue
        a=nw
        xm=s1/a
        ym=s2/a
             
C  I calculate STANDARD DEVIATION
        s1=0.0
        s2=0.0
        do 2 l=k,nnw
          s1=s1+(x(l)-xm)*(x(l)-xm)
          s2=s2+(y(l)-ym)*(y(l)-ym)
    2   continue
        xd=sqrt(s1/a)
        yd=sqrt(s2/a)

c   I calculate the cross sum 

        s1=0.0
        do 3 l=k,nnw
          s1=s1+(x(l)-xm)*(y(l)-ym)
    3   continue
        sxy=s1/a

c I calculate factor of correlation 

        r(k)=sxy/(xd*yd)
C        write(*,*)r(k)
   99 continue

C I write results 

      open(3,form='unformatted')

      ia(4)=nf
      write(3)(ia(l),l=1,4)
      write(3)(r(l),l=1,nf)
                         
      close(1)
      close(2)
      close(3)
      print *,'REGULAR END!'
      end
EOF
f77 corwindow.f -o corwindow.x
corwindow.x $1 $2 $3 << M
$4
M
cp fort.3 $3
rm fort*
rm corwindow.x corwindow.f
exit