anchor.diff:  
   ! For calculation of dE/dX due to ionization
   
  -      double precision rho(4),z(4),Ime(4),anuc_cry2(4) !Daniele: material parameters for dE/dX calculation
  -      double precision k,re,me,mp !Daniele: parameters for dE/dX calculation (const,electron radius,el. mass, prot.mass)
  -
  -      double precision emr_cry(4) !nuclear radius for Ruth scatt.
  +      real(kind=fPrec) rho(4),z(4),Ime(4),anuc_cry2(4) !Daniele: material parameters for dE/dX calculation
  +      real(kind=fPrec) k,re,me,mp !Daniele: parameters for dE/dX calculation (const,electron radius,el. mass, prot.mass)
  +
  +      real(kind=fPrec) emr_cry(4) !nuclear radius for Ruth scatt.
   
         common/ion/rho,z,Ime,k,re,me,mp,anuc_cry2,emr_cry
   
  
positive.code:       data dP/1.92d-10/   !distance between planes (110) [m]
        data u1/0.075d-10/  !thermal vibrations amplitude
  
        double precision tlcut_cry,hcut_cry(4)   !param. for Ruth scatt.
        double precision cgen_cry(200,4)
        integer mcurr_cry
  
        common/ruth_scat_cry/tlcut_cry,hcut_cry
  
negative.code: 
        use mod_ranlux
        use mod_funlux
        use floatPrecision
  
  !. Moving in amorphous substance...........................
        IMPLICIT none
        integer IS,NAM
        real(kind=fPrec) DZ,DEI,DLY,DLr, XP,YP,PC
        real(kind=fPrec) DLAI(4),SAI(4),DES(4)
        real(kind=fPrec) DLRI(4),DLYI(4),AI(4)
        real(kind=fPrec) AM(30),QP(30),NPAI
        real(kind=fPrec) Dc(4),eUm(4)
        real(kind=fPrec) DYA,W_p
  !      REAL RNDM4
  !         REAL      RAN_GAUSS
        CHARACTER(LEN=50) PROC              !string that contains the physical process
  