anchor.diff:        teta = sqrt(t)/pc
       end if
   
  -    tx = (teta*ran_gauss(one))*c1e3
  -    tz = (teta*ran_gauss(one))*c1e3
  +    tx = (teta*ran_gauss(zero))*c1e3
  +    tz = (teta*ran_gauss(zero))*c1e3
   
       ! Change p angle
       xp = xp + tx
  
positive.code:         if(Ldech < Length-Lrefl) then
  
            iProc = proc_DC
            Dxp   = Ldech/c_rcurv + (half*ran_gauss(one))*xpcrit
            x     = x + Ldech*(sin_mb(half*Dxp+xp)) ! Trajectory at channeling exit
            y     = y + Sdech*yp
            xp    =  Dxp
  
negative.code:         x     = x + xp*Srefl
          y     = y + yp*Srefl
          Dxp   = Ang_avr
          xp    = xp + Dxp + Ang_rms*ran_gauss(zero)
          x     = x  + (half*xp)*(s_length - Srefl)
          y     = y  + (half*yp)*(s_length - Srefl)
  
  