anchor.diff:           integer         Nev
            integer         j
            integer         nabs
  -         integer         NHIT
  +         integer         nhit
            integer         MAX_NPART
            integer         NP
   !        PARAMETER         (MAX_NPART=1500000)
  
positive.code: 
    end do
  
  end
  
  ! ================================================================================================ !
  !  Subroutine for the movements of the particles in the crystal
  
negative.code:          real(kind=fPrec)            p0
           real(kind=fPrec)            zlm
           real(kind=fPrec)            x,xp
           real(kind=fPrec)            x00
           real(kind=fPrec)            z
           real(kind=fPrec)            z00
  