For the calculation of steep water waves, a Lagrangian method is presented. This makes it easier to derive and understand the complicated non-linear structure of the equations of motion ( EOM ), for the behaviour of water surfaces. In addition, this formulation offers the possibility to model an EOM which describes water waves with deference to variation in depth.