gum = ws.gumbel_r.fit(Hs)