ind = np.where(((y > Ly/4-0.5*dy) & (y < Ly/4+0.5*dy)))