In Coles book 
describes how to estimate the standard errors of the return periods based on the delta method. He does not have the fill formula though. Are there any packages that do that automatically? Although I would prefer to do it myself. Has anyone implemented that? Preferably in R:)Thank you!