anchor.diff:  
         use mod_ranlux
         use mod_funlux
  +      use floatPrecision
   
         IMPLICIT none
         integer IS
  -      double precision PC,DZ,EnLo
  -      double precision rho(4),z(4),Ime(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 enr,mom,betar,gammar,bgr !Daniele: energy,momentum,beta relativistic, gamma relativistic
  -      double precision Tmax,plen !Daniele: maximum energy tranfer in single collision, plasma energy (see pdg)
  -      double precision anuc_cry2(4)
  -      double precision emr_cry(4)
  -      double precision thl,Tt,cs_tail,prob_tail
  -      double precision ranc
  +      real(kind=fPrec) PC,DZ,EnLo
  +      real(kind=fPrec) rho(4),z(4),Ime(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) enr,mom,betar,gammar,bgr !Daniele: energy,momentum,beta relativistic, gamma relativistic
  +      real(kind=fPrec) Tmax,plen !Daniele: maximum energy tranfer in single collision, plasma energy (see pdg)
  +      real(kind=fPrec) anuc_cry2(4)
  +      real(kind=fPrec) emr_cry(4)
  +      real(kind=fPrec) thl,Tt,cs_tail,prob_tail
  +      real(kind=fPrec) ranc
   !      REAL RNDM4
   
         common/ion/rho,z,Ime,k,re,me,mp,anuc_cry2,emr_cry
  
positive.code:             !           write(*,*)'s tan',Rcurv*sin(XP_tangent)
              !           write(*,*) 's tot', c_length,Rcurv*sin(cry_bend)
                         if ( XP .ge. XP_tangent  ) then
              
              ! if it hits the crystal, calculate in which point and apply the
              ! transformation and drift to that point
                           a_eq=(1.+xp**2)
  
negative.code:             !           write(*,*)'s tan',Rcurv*sin(XP_tangent)
              !           write(*,*) 's tot', c_length,Rcurv*sin(cry_bend)
                         if ( XP .le. XP_tangent  ) then
  
              ! if it hits the crystal, calculate in which point and apply the
              ! transformation and drift to that point
                           a_eq=(1.+xp**2)
  