anchor.diff:    character(len=mFileName), private, save :: dist_distFile  = " "     ! File name for reading the distribution
     character(len=mFileName), private, save :: dist_echoFile  = " "     ! File name for echoing the distribution
   
  -  integer,          allocatable, private, save :: dist_colFormat(:) ! The format of the file columns
  -  real(kind=fPrec), allocatable, private, save :: dist_colScale(:)  ! Scaling factor for the file columns
  -  integer,                       private, save :: dist_nColumns = 0 ! The number of columns in the file
  -  integer,                       private, save :: dist_readCols = 0 ! The number of columns read from file
  -
  -  character(len=:), allocatable, private, save :: dist_partLine(:)  ! PARTICLE definitions in the block
  -  integer,                       private, save :: dist_numPart = 0  ! Number of PARTICLE keywords in block
  +  integer,          allocatable, private, save :: dist_colFormat(:)   ! The format of the file columns
  +  real(kind=fPrec), allocatable, private, save :: dist_colScale(:)    ! Scaling factor for the file columns
  +  integer,                       private, save :: dist_nColumns = 0   ! The number of columns in the file
  +
  +  character(len=:), allocatable, private, save :: dist_partLine(:)    ! PARTICLE definitions in the block
  +  integer,                       private, save :: dist_nParticle = 0  ! Number of PARTICLE keywords in block
  +  
  +  integer,                       private, save :: dist_numPart   = 0  ! Number of actual particles generated or read
   
     !  Column formats
     ! ================
  
positive.code: 
    integer i, j
  
    do i=1,dist_readCols
      select case(dist_colFormat(i))
      case(dist_fmtPX)
        yv1(1:napx) = (yv1(1:napx)/ejfv(1:napx))*c1e3
  
negative.code:     call dist_appendFormat(dist_fmtMASS, c1e3)      ! Particle mass
      call dist_appendFormat(dist_fmtP, c1e3)         ! Particle momentum
      call dist_appendFormat(dist_fmtDT, one)         ! Time delay
  
    case default
      write(lerr,"(a)") "DIST> ERROR Unknown column format '"//trim(fmtName)//"'"
  