anchor.diff:    logical, private, save :: dist_norm4D     = .false.
     logical, private, save :: dist_norm6D     = .false.
   
  +  ! Arrays to send to DISTlib
  +  real(kind=fPrec), allocatable, private, save :: dist_partCol1(:)    ! Ends up in array xv1
  +  real(kind=fPrec), allocatable, private, save :: dist_partCol2(:)    ! Ends up in array yv1
  +  real(kind=fPrec), allocatable, private, save :: dist_partCol3(:)    ! Ends up in array xv2
  +  real(kind=fPrec), allocatable, private, save :: dist_partCol4(:)    ! Ends up in array yv2
  +  real(kind=fPrec), allocatable, private, save :: dist_partCol5(:)    ! Ends up in array sigmv
  +  real(kind=fPrec), allocatable, private, save :: dist_partCol6(:)    ! Ends up in array dpsv, evj, and ejfv
  +  integer,                       private, save :: dist_partFmt(6) = 0 ! The format used for each column
  +
   contains
   
   ! ================================================================================================ !
  
positive.code:       unitScale = c1e9
      case("[ps]")
        unitScale = c1e12
      case("[none]")
        unitScale = one
      case default
        goto 100
  
negative.code: subroutine dist_generateDist
  
    use crcoall
    use mod_alloc
    use mod_common
    use mod_particles
    use numerical_constants
  
    call alloc(dist_partCol1, napx, zero, "dist_partCol1")
    call alloc(dist_partCol2, napx, zero, "dist_partCol2")
    call alloc(dist_partCol3, napx, zero, "dist_partCol3")
    call alloc(dist_partCol4, napx, zero, "dist_partCol4")
    call alloc(dist_partCol5, napx, zero, "dist_partCol5")
    call alloc(dist_partCol6, napx, zero, "dist_partCol6")
  
    if(dist_nParticle > 0) then
      call dist_parseParticles
  