[dr_zoidberg]

And for this particular use case, numpy has .clip(), that takes:

    np.clip(ndarray, lower_bound, upper_bound)

And handles NaNs correctly.