anchor.diff:  !		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+1,tt1*100.0/TT,"% complete.  n",beami,"v",beamj
  +    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
   
  
positive.code: 
  type params
  
  
      real(wp) :: pi, me,electroncharge,mi, Kb,eV,epsilon0, Te, Ti, nuic, nuec, n0, nbeam, &
              vbeam_ev, vbeam, tetabeam, endTime, Tstep, TT, res, TT_res, L, N, Xstep, QW, &
              SEED, eta, ve, Cs, omegae, lambdaD
      integer :: Z
  
    	! nbeam.at(beami); vbeam_ev.at(beamj); vbeam.at(beamj); tetabeam.at(beamj);
  !  	//have to include other parameters regarding the Kappa distribution
  
  end type params
  
negative.code: do realization=1,QW
    cte2=omegae/2.0_wp/n0/1
  
    call system_clock(clock)
    seed = clock + 37 * [ (i - 1, i = 1, nseed) ]
    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='unknown',action='write')
    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='unknown',action='write')
    print *,'writing to',trim(argv)
  
  !   main loops
  
  