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. The vertical averaged currents (barotropic) are defined by