MODULE blas_precision
      INTEGER, PARAMETER :: sp = KIND(0.0E0)
      INTEGER, PARAMETER :: dp = KIND(0.0D0)
    END MODULE blas_precision