ind = np.where(((x > Lx/4-0.5*dx) & (x < Lx/4+0.5*dx)))