anchor.diff:  print *,"vbeam=",vbeam
   print *,"tetabeam=",tetabeam
   
  +call random_seed(size=nseed)
  +allocate(seed(nseed))
  +
   do beami=1,Nnbeam
   do beamj=1,Nvbeam
   
  -! TODO: may have blank spaces in filename
  +  parameters%pi = pi
  +  parameters%me = me
  +  parameters%electroncharge = electroncharge
  +  parameters%mi = mi
  +  parameters%Z = Z
  +  parameters%Kb = Kb
  +  parameters%eV = eV
  +  parameters%epsilon0= epsilon0
  +  parameters%Te = Te
  +  parameters%Ti = Ti
  +  parameters%nuic = nuic
  +  parameters%nuec= nuec
  +  parameters%n0 = n0
  +  parameters%nbeam = nbeam(beami)
  +  parameters%vbeam_ev = vbeam_ev(beamj)
  +  parameters%vbeam = vbeam(beamj)
  +  parameters%tetabeam = tetabeam(beamj)
  +  parameters%endTime = endTime
  +  parameters%Tstep = Tstep
  +  parameters%TT = TT
  +  parameters%res = res
  +  parameters%TT_res = TT_res
  +  parameters%L = L
  +  parameters%N = N
  +  parameters%Xstep = Xstep
  +  parameters%QW = QW
  +  parameters%eta = eta
  +  parameters%ve = ve
  +  parameters%Cs = Cs
  +  parameters%omegae = omegae
  +  parameters%lambdaD = lambdaD
  +
     write(argv,'(A,I0.3,A,I0.3)')  odir//"/parameters_n" , beami, "_v" , beamj
     ofn = trim(argv)
   
  -  open(newunit=u,file=ofn,status='unknown',action='write')
  +  open(newunit=u, file=ofn, form='formatted',status='unknown',action='write')
     write(u,*) parameters
  +  write(u,*) SEED
     close(u)
  -  print *, "Wrote",parameters,"to",ofn
  +  print *, "Wrote",parameters,SEED,"to",ofn
   
     do ii=1,N
       p(ii)=ii-N/2
  
positive.code: call get_command_argument(1,argv)
  odir = trim(argv)
  print *,'writing output to', odir
  
  allocate(beamev(argc-1))
  do i = 2,argc
  
negative.code:     output1(ii,12) = Source_factor_n(ii)
    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
  