MODULE blas_kinds
      PUBLIC
      INTEGER, PARAMETER :: dp = KIND(1.0D0)
      INTEGER, PARAMETER :: sp = KIND(1.0)
    END MODULE blas_kinds