anchor.diff:        double precision s_length             !element length along s
         double precision Alayer               !amorphous layer [m]
         integer C_orient                      !crystal orientation
  -      integer IS                            !index of the material
  +      integer IS,j                            !index of the material
   !      integer counter
         double precision  DLRI(4),DLYI(4),AI(4),DES(4)!cry parameters:see line~270
         double precision DESt                  ! Daniele: changed energy loss by ionization now calculated and not tabulated
  
positive.code:         lint(j)   = -1.
          indiv(j)  = -1.
  
          idx_proc = name(j)-100*samplenumber
  
          if (ITURN .eq. 1) then                                                          !daniele
                  bool_proc_old(idx_proc)=-1                                                     !daniele
  
negative.code: !.**************************************************************************
  !     SUBROUTINE FOR THE MOVEMENTS OF THE PARTICLES IN THE CRYSTAL
  !.**************************************************************************
        SUBROUTINE CRYST(IS,x,xp,y,yp,PC,Length,j)
  
        use mod_ranlux
        use mod_funlux
        use mod_common_main
  !
  !     Simple tranport protons in crystal 2
  !-----------------------------------------------------------C
  