anchor.diff:    nn(:,N-N/2+1:N,2) = -nn(:,:N/2,2) ! yes minus
     nn(:,N/2,:)=0.0_wp
   
  -  counter1=0
     do tt1=1,TT
   
   !		int c0=(tt1-1) % 3;
  
positive.code:       cte1=1.5*omegae*(lambdaD*k(pp))*(lambdaD*k(pp))
  			!cte1=1.5*Kb*Te/me/omega_off*k(pp)*k(pp)-(pow(omega_off,2)-pow(omegae,2))/2.0/omega_off;
        k4(pp,1)=Tstep*(cte1*(EE(c1,pp,2)+k3(pp,2)-SSE(pp,1)*Tstep)-nuE(pp)*(EE(c1,pp,1)+k3(pp,1)+SSE(pp,2)*Tstep)+cte2*CC(2));
        k4(pp,2)=Tstep*((-1)*cte1*(EE(c1,pp,1)+k3(pp,1)+SSE(pp,2)*Tstep)-nuE(pp)*(EE(c1,pp,2)+k3(pp,2)-SSE(pp,1)*Tstep)-cte2*CC(1));
  
        EE(c2,pp,1)=EE(c1,pp,1)+(k1(pp,1)+2.0*k2(pp,1)+2.0*k3(pp,1)+k4(pp,1))/6.0+SSE(pp,2)*Tstep
        EE(c2,pp,2)=EE(c1,pp,2)+(k1(pp,2)+2.0*k2(pp,2)+2.0*k3(pp,2)+k4(pp,2))/6.0-SSE(pp,1)*Tstep
        EE(c2,N/2,:)=0.0_wp
  
  
  
negative.code: 
  character(:), allocatable :: odir, ofn
  character(256) :: argv
  integer :: argc,i, ii,iij1, q,tt1,c1,c2,beami, beamj, realization,u,uEE,uNN
  real(wp) :: tic,toc
  
  !---- main loop variables
  