def dampHighCoverage(x):
cutoff = min(500, utl.ceilto(x[x.name].D.quantile(0.999), 10))
x.loc[(x[x.name].D > cutoff).values, (x.name[0], x.name[1], x.name[2]) + ("C",)] = 0
x.loc[(x[x.name].D > cutoff).values, (x.name[0], x.name[1], x.name[2]) + ("D",)] = 0
return x
