f22[mask] = Z[Ly1[mask].astype("int"), Lx1[mask].astype("int")]