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: !		long double omega_off=omegae+2*pi*300000;
  
  		! update display every 50th iteration
      if (mod(tt1,50) == 0) print '(A,I0.3,F5.2,A,I0.3,A,I0.3)',"Realization: ",&
          realization,tt1*100.0/TT,"% complete.  n",beami," v",beamj
  
      do pp=1,N
  
  
negative.code: 
      enddo ! iij2 N
  
      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
  