anchor.diff:    ! Check that particle hit the crystal
     if(x >= zero .and. x < c_xmax) then
   
  -    call cry_interact(mat,x,xp,z,zp,p,cry_length)
  +    ! MISCUT first step: P coordinates
  +    s_P = (c_rcurv-c_xmax)*sin_mb(-c_miscut)
  +    x_P = c_xmax + (c_rcurv-c_xmax)*cos_mb(-c_miscut)
  +    write(*,*) "P front", s_P, x_P
  +
  +    call cry_interact(mat,x,xp,z,zp,p,cry_length,s_P,x_P)
       s   = c_rcurv*c_sBend
       zlm = c_rcurv*c_sBend
       if(iProc /= proc_out) then
  
positive.code:   real(kind=fPrec) x_temp,x_shift,x_rot,x_int
    real(kind=fPrec) xp_temp,xp_shift,xp_rot,xp_int,xp_tangent
    real(kind=fPrec) tilt_int,shift,delta,a_eq,b_eq,c_eq
  
    s_temp     = zero
    s_shift    = zero
  
negative.code:   real(kind=fPrec), intent(inout) :: yp
    real(kind=fPrec), intent(inout) :: pc
    real(kind=fPrec), intent(in)    :: length
    real(kind=fPrec), intent(in)    :: s_P
    real(kind=fPrec), intent(in)    :: x_P
  
    integer nam,zn                        ! Switch on/off the nuclear interaction (NAM) and the MCS (ZN)
    real(kind=fPrec) ymax,ymin            ! Crystal geometrical parameters
  