anchor.diff:    sRot   = sin_mb(c_rotation)
     cRRot  = cos_mb(-c_rotation)
     sRRot  = sin_mb(-c_rotation)
  +
  +  !Set energy and nucleon change variables as with the coupling
  +  nnuc0 = 0
  +  ien0  = zero
  +  nnuc1 = 0
  +  ien1  = zero
   
     do j=1,napx
   
  
positive.code:       xp_in(j) = xp*cRRot + zp*sRRot
        yp_in(j) = zp*cRRot - xp*sRRot
  
        if(icoll == ipencil .and. iturn == 1 .and. pencil_distr /= 3) then
          x00      = mirror * x00
          x_in(j)  = x00*cRRot + z00*sRRot
  
negative.code:   x_in, xp_in, y_in, yp_in, p_in, s_in, enom, lhit_pos, lhit_turn, part_abs_pos_local,             &
    part_abs_turn_local, impact, indiv, lint, onesided, nhit_stage, j_slices, nabs_type, linside)
  
    use, intrinsic :: iso_fortran_env, only : int16
    use parpro
    use crcoall
    use coll_db
    use coll_common
    use coll_crystal, only : cry_doCrystal
    use coll_materials
    use mod_common, only : iexact, napx, unit208
    use mod_common_main, only : partID, naa
    use mathlib_bouncer
    use mod_ranlux
  #ifdef HDF5
    use hdf5_output
  #endif
  #ifdef ROOT
    use root_output
  #endif
  
    integer,          intent(in)    :: icoll        ! Collimator ID
  