anchor.diff:  
     end if
   
  -  if(C_orient == 2) then
  +  if(c_orient == 2) then
       Ang_avr = Ang_avr*0.93_fPrec
       Ang_rms = Ang_rms*1.05_fPrec
       xpcrit  = xpcrit*0.98_fPrec
  
positive.code:     indiv(j)  = -one
  
      ! CRY ---------------------
      iProc       = proc_out
      cry_proc(j) = proc_out
      ! CRY ---------------------
  
      ! CRY ---------------------
      nabs = 0
      s    = 0
      ! CRY ---------------------
  
      x   = x_in(j)
  
negative.code:   real(kind=fPrec) xpcrit               ! Critical angle for curved crystal[rad]
    real(kind=fPrec) xpcrit0              ! Critical angle for str. crystal [rad]
    real(kind=fPrec) Rcrit                ! Critical curvature radius [m]
    real(kind=fPrec) ratio                ! X=c_rcurv/Rcrit
    real(kind=fPrec) TLdech2              ! Tipical dechanneling length(1) [m]
    real(kind=fPrec) TLdech1              ! Tipical dechanneling length(2) [m]
    real(kind=fPrec) tdech, Ldech,Sdech   ! Angle, lenght, and S coordinate of dechanneling point
    real(kind=fPrec) Rlength, Red_S       ! Reduced length/s coordinate (in case of dechanneling)
    real(kind=fPrec) Am_length            ! Amorphous length
    real(kind=fPrec) Length_xs, Length_ys ! Amorphous length
    real(kind=fPrec) xp_rel               ! Xp-c_miscut angle in mrad
    real(kind=fPrec) alpha                ! Par for new chann prob
    real(kind=fPrec) Pvr                  ! Prob for VR->AM transition
  
  