In order to save computational time, the equations can be solved by the mode splitting method [13,17]. This method involves
separating out the external (vertically averaged) and the internal (three dimension) modes and solving each of them
separately at appropriate time steps.