anchor.diff:           real(kind=fPrec)      C_TILT(2)      !Tilt in radian
            real(kind=fPrec)      C_TILT0(2)      !Tilt in radian
            real(kind=fPrec)      cry_bend
  +         real(kind=fPrec)      cRot, sRot, cMRot, sMRot
   
   !
   !
  
positive.code: 
             IF (.NOT. CHANGED_TILT1(ICOLL) .AND. MIRROR.GT.0.) THEN
                     WRITE (*,*) 'NEVER!!!'
                     C_TILT(1) = XP_PENCIL0*COS(C_ROTATION)+SIN(C_ROTATION)*YP_PENCIL0
                     WRITE(*,*) 'INFO> Changed tilt1  ICOLL  to  ANGLE  ',ICOLL, C_TILT(1)
  !
                     CHANGED_TILT1(ICOLL) = .true.
             ELSEIF (.NOT. CHANGED_TILT2(ICOLL) .AND. MIRROR.LT.0.) THEN
                     C_TILT(2) = -1.*(XP_PENCIL0*COS(C_ROTATION)+SIN(C_ROTATION)*YP_PENCIL0)
                     WRITE(*,*) 'INFO> Changed tilt2  ICOLL  to  ANGLE  ',ICOLL, C_TILT(2)
  !
                     CHANGED_TILT2(ICOLL) = .true.
  
negative.code:       n_chan  = 0          !valentina :initialize to zero the counters for crystal effects
        n_VR    = 0          !
        n_amorphous = 0      !
  
        cRot  = cos(c_rotation)
        sRot  = sin(c_rotation)
        cMRot = cos(-c_rotation)
        sMRot = sin(-c_rotation)
  
        do j = 1, nev
  !
  