f21[mask] = Z[Ly0[mask].astype("int"), Lx1[mask].astype("int")]