anchor.diff:    use mathlib_bouncer
     use numerical_constants
   #ifdef HDF5
  +  use mod_alloc
     use hdf5_output
   #endif
   
  
positive.code:   iLost    = 0
    isDiff   = .false.
    updateE  = .false.
  
    ! Check if we're dumping the density diagnostics data on this element/turn
    densDump = scatter_dumpDensity(1) == iElem .and. scatter_dumpDensity(2) == nTurn
  
negative.code: ! =================================================================================================
  !  K. Sjobak, V.K. Berglyd Olsen, BE-ABP-HSS
  !  Created: 2017-08
  !  Updated: 2019-11-20
  ! =================================================================================================
  subroutine scatter_thin(iStru, iElem, nTurn)
  
    use crcoall
    use mod_time
    use mod_units
    use mod_ranecu
    use mod_common
  