      Module numerics
! ---------------------------------------------------------------------
      INTEGER,  PARAMETER  ::  sp = kind(1.0  )     
      INTEGER,  PARAMETER  ::  dp = kind(1.0D0)
!
      INTEGER,  PARAMETER  ::  rp = dp
!
      INTEGER   ::  ierr = 0
!
      REAL(rp), PARAMETER  ::  un   = 1.0_rp
      REAL(rp), PARAMETER  ::  zero = 0.0_rp
! ---------------------------------------------------------------------
      End Module numerics
