changed dispersion implementation to that using NonlinearSolve, which greatly enhanced accuracy for super shallow water waves