ind = np.flatnonzero(1 - np.isfinite(Hw))