anchor.diff:  
   !      ichoix=3
   
  -      if (ichoix.eq.1) then
  -
  -        PROC = 'ch_absorbed'            !deep inelastic, impinging p disappeared
  +      if (ichoix.eq.1) then  !deep inelastic, impinging p disappeared
  +
           iProc = proc_ch_absorbed
   
         endif
   
         if (ichoix.eq.2) then          ! p-n elastic
  -        PROC = 'ch_pne'
           iProc = proc_ch_pne
   !        write(*,*) "bn(IS)"
           bn(IS) = bnref(is2) * cs(0,IS) / csref_tot_rsc
  
positive.code:             !              write(*,*) "p after exit", p
                            s=Rcurv*sin(cry_bend-tilt_int)
                            zlm=Rcurv*sin(cry_bend-tilt_int)
                            if (PROC(1:3).ne.'out')then
                               X_rot=X_int
                               S_rot=S_int
                               XP_rot=XP_int
                               S_shift=S_rot*cos(-Cry_tilt)+X_rot*sin(-Cry_tilt)
  
negative.code:       endif
  
        if ( ichoix .eq. 5 ) then
          iProc = proc_ruth
             length_cry=1
             call funlux( cgen_cry(1,IS) ,xran_cry,length_cry)
  