anchor.diff:    logical,                  private, save :: dist_echo      = .false. ! Echo the read distribution?
     logical,                  private, save :: dist_hasFormat = .false. ! Whether the FORMAT keyword exists in block or not
     logical,                  private, save :: dist_libRead   = .false. ! Read file with dist library instead of internal reader
  -  integer,                  private, save :: dist_updtEFrom = 3       ! The parameter sent to part_updatePartEnergy after DIST
     character(len=mFileName), private, save :: dist_distFile  = " "     ! File name for reading the distribution
     character(len=mFileName), private, save :: dist_echoFile  = " "     ! File name for echoing the distribution
   
  
positive.code:   end if
  
    pstop(1:napx) = .false.
    ejf0v(1:napx) = ejfv(1:napx)
    mtc(1:napx)   = (nqq(1:napx)*nucm0)/(qq0*nucm(1:napx))
  
    ! If we have no energy arrays, we set all energies from deltaP = 0, that is reference momentum/energy
    call part_updatePartEnergy(dist_updtEFrom,.false.)
  
    if(dist_readIonA .and. dist_readIonZ .and. .not. dist_readPDGID) then
      do j=1,napx
  
negative.code: !  V.K. Berglyd Olsen, BE-ABP-HSS
  !  Created: 2019-07-05
  !  Updated: 2019-07-10
  !
  !    This routine splits the unit from the format description and adds the format and the scaling
  !  factor to the list of format columns for later file parsing.
  ! ================================================================================================ !
  subroutine dist_setColumnFormat(fmtName, fErr)
  
  