anchor.diff:        CC(:)=0.0_wp
   
         do q=LL,UU
  +        ! no vect
           CC(1)=CC(1)+(EE(c1,q+N/2,1)+k3(q+N/2,1))*nn(c1,p(pp)-q+N/2,1)-(EE(c1,q+N/2,2)+k3(q+N/2,2))*nn(c1,p(pp)-q+N/2,2);
           CC(2)=CC(2)+(EE(c1,q+N/2,1)+k3(q+N/2,1))*nn(c1,p(pp)-q+N/2,2)+(EE(c1,q+N/2,2)+k3(q+N/2,2))*nn(c1,p(pp)-q+N/2,1);
         end do
  
positive.code: 
  character(:), allocatable :: odir, ofn
  character(256) :: argv
  integer :: argc,i, ii,iij1,counter1, q,tt1,c1,c2,beami, beamj, realization,u,uEE,uNN
  real(wp) :: tic,toc
  
  !---- main loop variables
  
negative.code:   nn(:,N-N/2+1:N,2) = -nn(:,:N/2,2) ! yes minus
    nn(:,N/2,:)=0.0_wp
  
    do tt1=1,TT
  
  !		int c0=(tt1-1) % 3;
  