anchor.diff:  
         real(kind=fPrec) emr_cry(4)
   
  -!      common /Par_Cry1/ Cry_length, Rcurv,C_xmax,C_ymax,Alayer,C_orient
  -      common /Par_Cry1/ Rcurv,C_xmax,C_ymax,Alayer,C_orient
  -      common /miscut/ miscut
  -      common /Proc2/PROC
         COMMON/NPC/     NAM,ZN
         COMMON/CRYS/    DLRI,DLYI,AI,DES
         COMMON/eUc/     eUm  !
  
positive.code: 
        function ruth_cry(t_cry)
          use floatPrecision
        implicit none
        integer nmat_cry,mcurr_cry
        parameter(nmat_cry=4)
  !      real(kind=fPrec) z(4),emr_cry(4),tlcut_cry,hcut_cry(4)
        real(kind=fPrec) tlcut_cry,hcut_cry(4)
  
        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)
  
        real(kind=fPrec) cgen_cry(200,4)
  
        common/ion/rho,z,Ime,k,re,me,mp,anuc_cry2,emr_cry
  
  !      common/ion/z,emr_cry
        common/ruth_scat_cry/tlcut_cry,hcut_cry
        common/ruth_scat_cry/cgen_cry,mcurr_cry
  
        real(kind=fPrec) ruth_cry,t_cry
        real(kind=fPrec) cnorm,cnform
  
negative.code: !        PARAMETER         (MAX_NPART=1500000)
  !
           real(kind=fPrec)            p0
           real(kind=fPrec)            zlm
           real(kind=fPrec)            x,xp
           real(kind=fPrec)            shift
  