def num_above(values, level):
    return np.sum(((values >= level) & np.isfinite(values)), dtype=np.float)
def num_above(values, level):
    return np.sum(((values >= level) & np.isfinite(values)), dtype=np.float)