      Module numerics
! ---------------------------------------------------------------------
      Integer, Parameter  ::  sp = kind(1.0  )
      Integer, Parameter  ::  dp = kind(1.0D0)
!
      Integer, Parameter  ::  rp = dp
      Integer, Parameter  ::  ip = 8
!
      Integer             ::  ntx, nty, ntz
!
      Real(rp)            ::  hx,  hy,  hz
!
      Real(rp)            ::  rlambda
      Real(rp)            ::  dt
!
      Real(rp), Parameter ::  one  = 1.0_rp
      Real(rp), Parameter ::  zero = 0.0_rp
!
      Real(rp)            ::  machep
! ---------------------------------------------------------------------
      End Module numerics