anchor.diff:  
       enddo ! iij2 N
   
  -    do iij2=1,N/2
  -      nn (iij1,N-iij2,1)=  nn(iij1,iij2,1)
  -      nn (iij1,N-iij2,2)= -nn(iij1,iij2,2)
  -    enddo ! iij2 N/2
  +    nn(iij1,N-1:N/2,1) = nn(iij1,N-1:N/2,1)
   
       nn (iij1,N/2,1)=0
       nn (iij1,N/2,2)=0
  
positive.code:   end do ! ii
  
    write(argv,'(A,I0.3,A,I0.3)') odir//"/output1_n",beami,"_v", beamj
    open(newunit=u,file=trim(argv),status='unknown',action='write')
  
    write(u,*) output1
  
    close(u)
    print *, "Wrote to", trim(argv)
  
  
  do realization=1,QW
  
negative.code:   call random_seed(put=seed)
  
    write(argv,'(A,I0.3,A,I0.3,A,I0.3)') odir//"/EE",realization,"_n",beami,"_v",beamj
    open(newunit=uEE,file=trim(argv), status='new',action='write',access='stream')
    print *,'writing to ',trim(argv)
  
    write(argv,'(A,I0.3,A,I0.3,A,I0.3)') odir//"/nn",realization, "_n" ,beami, "_v",beamj
    open(newunit=uNN,file=trim(argv), status='new',action='write',access='stream')
    print *,'writing to ',trim(argv)
  
  !   main loops
  
  