anchor.diff:          x = x + xp*Srefl
           y = y + yp*Srefl
   
  -        TLdech2 = (const_dech/c1e1)*pc*(one-one/ratio)**2       ! Updated typical dechanneling length(m)
  -        Ldech   = TLdech2*(sqrt(c1m2 - log(rndm4())) - c1m1)**2 ! Updated DC length
  +        TLdech2 = (const_dech/c1e1)*pc*(one-one/ratio)**2          ! Updated typical dechanneling length(m)
  +        Ldech   = TLdech2*(sqrt(c1m2 - log_mb(rndm4())) - c1m1)**2 ! Updated DC length
           tdech   = Ldech/Rcurv
  -        Sdech   = Ldech*cos(xp + half*tdech)
  +        Sdech   = Ldech*cos_mb(xp + half*tdech)
   
           if(Ldech < Length-Lrefl) then
   
             iProc = proc_DC
             Dxp   = Ldech/Rcurv + half*ran_gauss(one)*xpcrit
  -          x     = x + Ldech*(sin(half*Dxp+xp)) ! Trajectory at channeling exit
  +          x     = x + Ldech*(sin_mb(half*Dxp+xp)) ! Trajectory at channeling exit
             y     = y + Sdech*yp
             xp    =  Dxp
             Red_S = s_length - Srefl - Sdech
  
positive.code:            X  = X + tiltangle*C_LENGTH
             XP = XP + tiltangle
  !
             X  = X - SIN(tiltangle) * C_LENGTH
           ENDIF
  !
  !++  Transform back to particle coordinates with opening and offset
  
negative.code: !               write(*,*) "debug - S cryRF 2" ,  S_rot
  !               write(*,*) "debug - X cryRF 2" ,  X_rot
  !               write(*,*) "debug - XP cryRF 2" ,  XP_rot
                 S_shift=S_rot *cos_mb(-Cry_tilt)+X_rot*sin_mb(-Cry_tilt)
                 X_shift=-S_rot*sin_mb(-Cry_tilt)+X_rot*cos_mb(-Cry_tilt)
                 XP_shift=XP_rot + Cry_tilt
  !     2nd: shift back the reference frame
                 if (Cry_tilt .lt. zero) then
  