module precision_module
    ! define SP as single precision kind
  integer, parameter, public :: SP=selected_real_kind(4)
    ! define DP as double precision kind
  integer, parameter, public :: DP = selected_real_kind(15)
end module precision_module