anchor.diff:  
     character(len=:), allocatable :: lnSplit(:)
     character(len=mInputLn) inLine
  -  integer i, nSplit, lineNo, fUnit, nRead
  -  logical spErr, cErr
  +  real(kind=fPrec) fmtFac
  +  integer i, nSplit, lineNo, fUnit, nRead, fmtID, fmtCol
  +  logical spErr, cErr, isValid
   
     write(lout,"(a)") "DIST> Reading particles from '"//trim(dist_distFile)//"'"
   
  
positive.code:   logical,          intent(inout) :: iErr
  
    character(len=:), allocatable   :: lnSplit(:)
    integer nSplit, i
    logical spErr, cErr
  
    call chr_split(inLine, lnSplit, nSplit, spErr)
    if(spErr) then
  
negative.code:     call dist_appendFormat(dist_fmtPY,       one,  4)
      call dist_appendFormat(dist_fmtZETA,     one,  5)
      call dist_appendFormat(dist_fmtDPP0,     one,  6)
      isValid = .true.
  
    case("NORM") ! 6D normalised coordinates
      call dist_appendFormat(dist_fmtXN,       one,  1)
  