anchor.diff:          lint(j)   = -1.
           indiv(j)  = -1.
   
  -        idx_proc = name(j)-100*samplenumber
  +        idx_proc = j
   
           if (ITURN .eq. 1) then                                                          !daniele
                   bool_proc_old(idx_proc)=-1                                                     !daniele
  
positive.code: !           write(*,*)'impact at s,x = ', s_impact,x_in0(j)
  !           write(*,*)'with angle xp = ',xp
  !           write(*,*)'s before', s
             CALL CRYST(mat-7,X,XP,Z,ZP,p,cry_length)
  !           write(*,*) "xp after crystal coll routine exit", XP
             s=Rcurv*sin(cry_bend)
             zlm=Rcurv*sin(cry_bend)
  
negative.code: !      write(*,*) 'Length [m]:', Length
  !      write(*,*) 'Random:', rndm4()
  !      write(*,*)'xp',xp,'x',x , 's', s
  !      DESt = 0
        s=0
  !      write(*,*) "s_length", Rcurv, length
        s_length=Rcurv*(sin(length/Rcurv)) !
  