ind2 = np.flatnonzero((data1[iy, ind] > data1[iy - 1, ind]) &
                                      (data1[iy, ind] > data1[iy + 1, ind]))