anchor.diff:    ! Ion Columns
     integer, parameter :: dist_fmtMASS        = 51 ! Particle mass
     integer, parameter :: dist_fmtCHARGE      = 52 ! Particle Charge
  -  integer, parameter :: dist_fmtIonA        = 53 ! Ion atomic number
  -  integer, parameter :: dist_fmtIonZ        = 54 ! Ion atomic charge
  +  integer, parameter :: dist_fmtIonA        = 53 ! Ion atomic mass
  +  integer, parameter :: dist_fmtIonZ        = 54 ! Ion atomic number
     integer, parameter :: dist_fmtPDGID       = 55 ! Particle PDG ID
   
   contains
  
positive.code:     call dist_appendFormat(dist_fmtE, one)           ! Particle energy
    case("P")
      call dist_appendFormat(dist_fmtP, one)           ! Particle momentum
    case("DE/E0")
      call dist_appendFormat(dist_fmtDEE0, one)        ! Relative particle energy (to reference particle)
    case("DP/P0")
      call dist_appendFormat(dist_fmtDPP0, one)        ! Relative particle momentum (to reference particle)
  
    case("X_NORM")
  
negative.code:     call dist_appendFormat(dist_fmtE_NORM, one)      ! Normalised particle energy
    case("P_NORM")
      call dist_appendFormat(dist_fmtP_NORM, one)      ! Normalised particle momentum
    case("DE/E0_NORM","DEE0_NORM")
      call dist_appendFormat(dist_fmtDEE0_NORM, one)   ! Normalised relative particle energy (to reference particle)
    case("DP/P0_NORM","DPP0_NORM")
      call dist_appendFormat(dist_fmtDPP0_NORM, one)   ! Normalised relative particle momentum (to reference particle)
  
    case("MASS","M")
  