INTEGER, PARAMETER :: DP = selected_real_kind(14,200)
INTEGER, PARAMETER :: SP = selected_real_kind(5,20)