anchor.diff:    character(len=:), allocatable :: lnSplit(:)
     character(len=mInputLn) inLine
     integer i, j, nSplit, lineNo, fUnit
  -  logical spErr, cErr, isFirst
  +  logical spErr, cErr
   
     write(lout,"(a)") "DIST> Reading particles from '"//trim(dist_distFile)//"'"
   
  -  dist_readCols = dist_nColumns
  -
  -  isFirst = .true.
  -  lineNo  = 0
  -  cErr    = .false.
  -  j       = dist_numPart ! PARTICLE definitions in the block are saved first
  +  lineNo = 0
  +  cErr   = .false.
  +  j      = dist_nParticle ! PARTICLE definitions in the block are saved first
   
     if(dist_hasFormat .eqv. .false.) then
       call dist_setColumnFormat("OLD_DIST",cErr)
  
positive.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
      call dist_appendFormat(dist_fmtCHARGE, one)     ! Particle charge
      call dist_appendFormat(dist_fmtPDGID, one)      ! Particle PDG ID
  
    case default
      write(lerr,"(a)") "DIST> ERROR Unknown column format '"//trim(fmtName)//"'"
  
negative.code: 
    integer i, j
  
    do i=1,dist_nColumns
      select case(dist_colFormat(i))
      case(dist_fmtPX)
        yv1(1:napx) = (yv1(1:napx)/ejfv(1:napx))*c1e3
  