anchor.diff:    end if
   
     if(hasParticles) then
  +    call dist_postprParticles
       call dist_finaliseDist
       call part_applyClosedOrbit
       if(dist_echo) then
  
positive.code:     nqq(1:napx) = nzz(1:napx)
    end if
  
    mtc(1:napx)   = (nqq(1:napx)*nucm0)/(qq0*nucm(1:napx))
    pstop(1:napx) = .false.
    ejf0v(1:napx) = ejfv(1:napx)
  
    if(dist_readIonA .and. dist_readIonZ .and. .not. dist_readPDGID) then
      do j=1,napx
        call CalculatePDGid(pdgid(j), naa(j), nzz(j))
      end do
    end if
  
  end subroutine dist_postprParticles
  
  ! ================================================================================================ !
  !  A. Mereghetti and D. Sinuela Pastor, for the FLUKA Team
  !  V.K. Berglyd Olsen, BE-ABP-HSS
  !  Last modified: 2018-10-31
  ! ================================================================================================ !
  subroutine dist_finaliseDist
  
    use parpro
    use mod_common
    use mod_common_track
    use mod_common_main
    use numerical_constants
  
    implicit none
  
    integer          :: j
    real(kind=fPrec) :: chkP, chkE
  
    ! Check existence of on-momentum particles in the distribution
    do j=1, napx
  
negative.code: 
  end subroutine dist_appendFormat
  
  ! ================================================================================================ !
  !  Parse PARTICLE Lines from DIST Block
  !  V.K. Berglyd Olsen, BE-ABP-HSS
  !  Created: 2019-07-08
  !  Updated: 2019-07-08
  ! ================================================================================================ !
  subroutine dist_parseParticles
  
    use mod_common
  