anchor.diff:              x = x + (half*Rlength)*XP
               y = y + (half*Rlength)*YP
   
  -            call calc_ion_loss_cry(is,pc,rlength,dest)
  +            call cry_calcIonLoss(is,pc,rlength,dest)
               pc = pc - dest*Rlength
             else
               Dxp = (Length-Lrefl)/c_rcurv
  
positive.code: !  Subroutine for the movements of the particles in the crystal
  !  Simple tranport protons in crystal 2
  ! ================================================================================================ !
  subroutine cryst(is,x,xp,y,yp,pc,length)
  
    use mod_ranlux
    use mod_funlux
  
negative.code:           x = x + (half*L_chan)*XP
            y = y + (half*L_chan)*YP
  
            call cry_calcIonLoss(is,pc,length,dest)
            pc = pc - dest*length ! energy loss to ionization [GeV]
          else
            Dxp = L_chan/c_rcurv + (half*ran_gauss(one))*xpcrit ! Change angle[rad]
  