anchor.diff:      output1(ii,12) = Source_factor_n(ii)
     end do ! ii
   
  -  write(argv,'(A,I0.3,A,I0.3)') odir//"output1_n",beami,"_v", beamj
  +  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
  
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+1,tt1*100.0/TT,"% complete.  n",beami,"v",beamj
  
      do pp=1,N
  
  
negative.code: real(wp), parameter :: Tstep=0.5e-7_wp ! simulation time steps
  integer, parameter :: TT= floor(endTime / Tstep)              ! floor(endTime/Tstep)+2;
  integer, parameter :: res=20
  integer, parameter :: TT_res = floor(endTime/Tstep/res)  !floor(endTime/Tstep/res);
  real(wp), parameter :: L=70.0_wp           ! simulation box length (meter)j
  integer, parameter :: N=2046          ! number of samples in L; should be devidable by 6
  real(wp), parameter :: Xstep= L / N
  