integer, parameter :: dp = kind(0.d0)