anchor.diff:  
         ! Rotate beam 2 in x,y plane so the crossing angle is only in x
         ! elRot = -atan2_mb(yFac-one,xFac-one)
  -      elRot = -atan2_mb(sigmaYeff-sigmaY,sigmaXeff-sigmaX)
  +      elRot = -atan2_mb(sigmaYeff - sigmaY, sigmaXeff - sigmaX)
         sRot  = sin_mb(elRot)
         cRot  = cos_mb(elRot)
   
  
positive.code: ! =================================================================================================
  !  K. Sjobak, V.K. Berglyd Olsen, BE-ABP-HSS
  !  Created: 2017-08
  !  Updated: 2018-11-12
  ! =================================================================================================
  subroutine scatter_init
  
  
negative.code: ! =================================================================================================
  !  V.K. Berglyd Olsen, K. Sjobak, BE-ABP-HSS
  !  Created: 2017-08
  !  Updated: 2018-11-20
  ! =================================================================================================
  subroutine scatter_getDensity(idPro, xPos, yPos, sPos, xProj, yProj, rhoOut)
  
  