In this thesis I describe a procedure for forecasting the propagation of ocean surface waves. I start by developing the third order Schrödinger equations in chapter 1, show how one can extract the initial state from surface data using the Fourier transform in chapter 2, and then implement and test the procedure in chapters 3-5.