A multiple threshold method for fitting the generalized Pareto distribution to rainfall time series