anchor.diff:            IF (tiltangle.GT.0.) THEN
               XP = XP - tiltangle
             ELSEIF (tiltangle.LT.0.) THEN
  -            X  = X + SIN(tiltangle) * C_LENGTH
  +            X  = X + sin_mb(tiltangle) * C_LENGTH
               XP = XP - tiltangle
             ENDIF
   !
  
positive.code:            CALL CRYST(mat,X,XP,Z,ZP,p,cry_length,j)
  !           write(*,*) "p after exit", p
             write(*,*) "xp after crystal coll routine exit", XP
             s=Rcurv*sin(cry_bend)
             zlm=Rcurv*sin(cry_bend)
             if(iProc /= proc_out) then
               NHIT = NHIT + 1
               LHIT(j) = ie
  
negative.code:       use coll_common, only : cry_proc, xp_pencil0, yp_pencil0, x_pencil, y_pencil, pencil_dx, ipencil
        use floatPrecision
        use coll_crystal
        use mathlib_bouncer
  
        IMPLICIT NONE
  !
  