anchor.diff:  
         ! If it hits the crystal, calculate in which point and apply the transformation and drift to that point
         a_eq  = one + xp**2
  -      b_eq  = (two*xp)*(x-c_rcurv)
  +      b_eq  = (two*xp)*(x - c_rcurv)
         c_eq  = -(two*x)*c_rcurv + x**2
         delta = b_eq**2 - four*(a_eq*c_eq)
         s_int = (-b_eq - sqrt(delta))/(two*a_eq)
   
  -      if(s_int < c_rcurv*sin_mb(cry_bend)) then
  +      if(s_int < c_rcurv*c_sBend) then
           ! Transform to a new reference system: shift and rotate
           x_int  = xp*s_int + x
           xp_int = xp
  
positive.code: ! ================================================================================================ !
  !  Subroutine for the movements of the particles in the crystal
  !  Simple tranport protons in crystal 2
  !   J         - number of element
  !   S         - longitudinal coordinate
  !   IS        - number of substance 1-4: Si,W,C,Ge(110)
  !   x,xp,y,yp - coordinates at input of crystal
  !   PC        - momentum of particle*c [GeV]
  !   W         - weigth of particle
  ! ================================================================================================ !
  subroutine cryst(is,x,xp,y,yp,pc,length,j)
  
    use mod_ranlux
    use mod_funlux
  
negative.code:         tilt_int = s_int/c_rcurv
          xp       = xp-tilt_int
  
          call cryst(mat,x,xp,z,zp,p,cry_length-(tilt_int*c_rcurv))
          s   = c_rcurv*sin_mb(cry_bend - tilt_int)
          zlm = c_rcurv*sin_mb(cry_bend - tilt_int)
          if(iProc /= proc_out) then
            x_rot    = x_int
            s_rot    = s_int
  