anchor.diff:      if(iseed2 < 0) iseed2 = iseed2+2147483399
       iz = iseed1-iseed2
       if(iz < 1) iz = iz+2147483562
  -    r(j) = real(iz,fPrec)*4.656613e-10_fPrec
  +    r(j) = real(iz,fPrec)*4.656613e-10_fPrec ! Note: this is still single precision
     end do
   
     if(mode == 1) then
  
positive.code: ! ================================================================================================ !
  module mod_ranecu
  
    implicit none
  
    integer, private, save :: iseed1 = 12345
    integer, private, save :: iseed2 = 67890
  
  contains
  
  
negative.code:     if(iseed2 < 0) iseed2 = iseed2+2147483399
      iz = iseed1-iseed2
      if(iz < 1) iz = iz+2147483562
      rvec(j) = real(iz,fPrec)*rScale
    end do
  
  end subroutine ranecuu
  