anchor.diff:      EnLo = EnLo/dz ! [GeV/m]
     end if
   
  -end subroutine calc_ion_loss_cry
  +end subroutine cry_calcIonLoss
   
   ! ================================================================================================ !
   !  Subroutine for the movement in the amorphous
   ! ================================================================================================ !
  -subroutine move_am(is,nam,dz,dei,dly,dlr,xp,yp,pc)
  +subroutine cry_moveAM(is,nam,dz,dei,dly,dlr,xp,yp,pc)
   
     use mod_ranlux
     use mod_funlux
  
positive.code:           x     = x + (half*xp)*Red_S
            y     = y + (half*yp)*Red_S
  
            call calc_ion_loss_cry(is,pc,srefl,dest)
            pc = pc - dest*Srefl ! "added" energy loss before capture
  
            call calc_ion_loss_cry(is,pc,sdech,dest)
            pc = pc - (half*dest)*Sdech ! "added" energy loss while captured
  
            call calc_ion_loss_cry(is,pc,red_s,dest)
            call move_am(is,nam,red_s,dest,dlyi(is),dlri(is),xp,yp,pc)
            x = x + (half*xp)*Red_S
            y = y + (half*yp)*Red_S
  
  
negative.code: !  Subroutine for the movements of the particles in the crystal
  !  Simple tranport protons in crystal 2
  ! ================================================================================================ !
  subroutine cry_interact(is,x,xp,y,yp,pc,length)
  
    use mod_ranlux
    use mod_funlux
  