anchor.diff:      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)
  
positive.code: 
    integer,                       private, save :: dist_numPart   = 0  ! Number of actual particles generated or read
  
    !  Column formats
    ! ================
    integer, parameter :: dist_fmtNONE        = 0  ! Column ignored
    integer, parameter :: dist_fmtPartID      = 1  ! Paricle ID
  
negative.code:     call dist_appendFormat(dist_fmtJZ,       one,  5)
      call dist_appendFormat(dist_fmtPhiZ,     one,  6)
      dist_distLib = .true.
      isValid = .true.
  
    case("IONS") ! The ion columns
      call dist_appendFormat(dist_fmtMASS,     c1e3, 0)
  