anchor.diff:      end if
       ichoix = i
   
  -    ! Do the interaction 
  +    ! Do the interaction
       select case(ichoix)
       case(1) ! deep inelastic, impinging p disappeared
         iProc = proc_ch_absorbed
  -    
  +
       case(2) ! p-n elastic
         iProc  = proc_ch_pne
         bn(IS) = bnref(is2) * cs(0,IS) / csref_tot_rsc
  
positive.code:       length_cry = 1
        call funlux(cgen_cry(1,is),xran_cry,length_cry)
        t = xran_cry(1)
      
      end select
  
      ! Calculate the related kick -----------
  
negative.code:   ! Can nuclear interaction happen?
    ! Rescaled nuclear collision length
    if(avrrho == zero) then
      nuc_cl_l = c1e6
    else
      nuc_cl_l = collnt_cry(is)/avrrho
    end if
  