jhpb(i)=0
enddo
C Initialize correlation arrays
- do i=-maxtor,maxtor
+ do i=1,maxres
do k=1,2
b1(k,i)=0.0
b2(k,i)=0.0
b1tilde(k,i)=0.0
c b2tilde(k,i)=0.0
do j=1,2
+C CC(j,k,i)=0.0
+C Ctilde(j,k,i)=0.0
+C DD(j,k,i)=0.0
+C Dtilde(j,k,i)=0.0
+ EE(j,k,i)=0.0
+ enddo
+ enddo
+ enddo
+ do i=-maxtor,maxtor
+ do k=1,2
+ do j=1,2
CC(j,k,i)=0.0
Ctilde(j,k,i)=0.0
DD(j,k,i)=0.0
Dtilde(j,k,i)=0.0
- EE(j,k,i)=0.0
enddo
- enddo
+ enddo
enddo
C
C Initialize timing.