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: 
    !  Column formats
    ! ================
    integer, parameter :: dist_fmtNONE        = 0
    integer, parameter :: dist_fmtPartID      = 1
    integer, parameter :: dist_fmtParentID    = 2
  
    ! Physical Coordinates
    integer, parameter :: dist_fmtX           = 11 ! Horizontal position
  
negative.code:   cErr = .false.
  
    if(dist_hasFormat .eqv. .false.) then
      call dist_setColumnFormat("OLD_DIST",cErr)
    end if
  
    call f_open(unit=dist_readUnit,file=dist_readFile,mode='r',err=cErr,formatted=.true.,status="old")
  