anchor.diff:    c_orient = cdb_cryOrient(icoll)
     c_miscut = cdb_cryMiscut(icoll)
     cry_bend = cry_length/c_rcurv
  +  c_cBend  = cos_mb(cry_bend)
  +  c_sBend  = sin_mb(cry_bend)
     c_cpTilt = cos_mb(cry_tilt)
     c_spTilt = sin_mb(cry_tilt)
     c_cnTilt = cos_mb(-cry_tilt)
  
positive.code:   real(kind=fPrec), private, save :: c_spTilt = zero ! Sine of positive crystal tilt
    real(kind=fPrec), private, save :: c_cnTilt = zero ! Cosine of negative crystal tilt
    real(kind=fPrec), private, save :: c_snTilt = zero ! Sine of negative crystal tilt
  
    ! Rutherford Scatter
    real(kind=fPrec), parameter     :: tlcut_cry = 0.0009982_fPrec
  
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
  