anchor.diff:  ! ----- Maybe a warning would be nice that c_tilt is overwritten !!!!!
   ! changed xp_pencil0(icoll) to xp_pencil0 due to definition mismatch
   ! this has to be solved if necassary and understood
  -!                 c_tilt(1) = xp_pencil0(icoll)*cos_mb(c_rotation)+         &
  -!     &                       sin_mb(c_rotation)*yp_pencil0(icoll)
  -!                 c_tilt(1) = xp_pencil0*cos_mb(c_rotation)+                &
  -!     &                       sin_mb(c_rotation)*yp_pencil0
  +!                 c_tilt(1) = xp_pencil0(icoll)*cRot+         &
  +!     &                       sRot*yp_pencil0(icoll)
  +!                 c_tilt(1) = xp_pencil0*cRot+                &
  +!     &                       sRot*yp_pencil0
   !                 write(*,*) "INFO> Changed tilt1  ICOLL  to  ANGLE  ",  &
   !     &                   icoll, c_tilt(1), j
   !                 changed_tilt1(icoll) = .true.
  
positive.code:   nhit    = 0
    fracab  = zero
    mirror  = one
  
  !==> SLICE here
  
  !- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    do j = 1, nev
  ! SR-GRD (04-08-2005):
  !   Don't do scattering process for particles already absorbed
  
negative.code: !     &                                   .and. mirror.lt.0.) then
  ! changed xp_pencil0(icoll) to xp_pencil0 due to definition mismatch
  ! this has to be solved if necassary and understood
  !                 c_tilt(2) = -1.*(xp_pencil0(icoll)*cRot+    &
  !     &                       sRot*yp_pencil0(icoll))
  !                 c_tilt(2) = -1.*(xp_pencil0*cRot+           &
  !     &                       sRot*yp_pencil0)
  !                 write(*,*) "INFO> Changed tilt2  ICOLL  to  ANGLE  ",  &
  !     &                   icoll, c_tilt(2), j
  !                 changed_tilt2(icoll) = .true.
  