anchor.diff:      end if
   
       if (counter1==20000) then
  -      write(uEE,*) total_EE
  -      write(unn,*) total_nn
  +      write(uEE) total_EE
  +      write(unn) total_nn
         counter1=0
       end if
   
     end do ! tt1
   
     if (counter1>0) then
  -    write(uEE,*) total_EE
  -    write(uNN,*) total_nn
  +    write(uEE) total_EE
  +    write(uNN) total_nn
     end if
   
   
  
positive.code:       vv(c2,pp,2)=vv(c1,pp,2)+(kv1(2)+2*kv2(2)+2*kv3(2)+kv4(2))/6+SSn(2)*Tstep;
        nn(c2,pp,1)=nn(c1,pp,1)+(kn1(1)+2*kn2(1)+2*kn3(1)+kn4(1))/6;
        nn(c2,pp,2)=nn(c1,pp,2)+(kn1(2)+2*kn2(2)+2*kn3(2)+kn4(2))/6;
        nn(c2,N/2,1)=0.0;
        nn(c2,N/2,2)=0.0;
  
        if (pp>=1) then
          nn(c2,N-pp,1)=nn(c2,pp,1);
  
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
  
  