anchor.diff:        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
   !
  
positive.code: !
  
  !          write(*,*) "debug - rotazione" ,c_rotation
            x  = x_in(j)*cos(c_rotation) +sin(c_rotation)*y_in(j)
            z  = y_in(j)*cos(c_rotation) -sin(c_rotation)*x_in(j)
            xp = xp_in(j)*cos(c_rotation)+sin(c_rotation)*yp_in(j)
            zp = yp_in(j)*cos(c_rotation)-sin(c_rotation)*xp_in(j)
  
  !          write(*,*) x, xp, z, zp
  
  
negative.code:          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
  
  !
  !
  