ihpb(i)=0
jhpb(i)=0
enddo
+C Initialize correlation arrays
+ do i=-maxtor,maxtor
+ 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
+ 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
C
C Initialize timing.
C