After INTERNAL and SUPERBEEUV have estimated and , MODESPLIT is called to update the U, V , W and ETA fields including the effects of all terms in the momentum equations and the equation of continuity, except the atmospheric pressure term.

In MODESPLIT the following steps are taken:

- The contributions in and are vertically integrated to produce terms and according to equation (55).
- The and fields are split into 2-D external or barotropic and 3-D internal or baroclinic parts according to equation (51).
- Then for each call to MODESPLIT N2D steps are used to propagate the depth
integrated velocities and and the water level with a
-method using DTE = DT/N2D as time step.
In each 2D time step the following actions are taken:
Estimate the field through a call to SMAGOR2D, if SMAG is true.
Include the effects of horizontal viscosity on the and
fields. (Calls to HORVISCUV2DPOM or HORVISCUV2D)
In the first substep, a forward method is used as a predictor method.
That is: = + DTE
,
where contain the spatially discretized terms of the -equation.
In the first substep, and are computed correspondingly with the forward method.
In each of the following 2D sub-steps, the leapfrog method is used as predictor for ,
, and .
That is: = +
.
The final values at time step are computed with a -method according to:
= +
must be in the range
to ensure stability.
- After the 2D steps, the velocity fields, that have been used to change the water level from ETAP to ETA over the N2D 2D time steps, are saved in the UADV and VADV arrays. It is essential that these non-divergent fields are used in the advection of scalar fields and momentum to maintain the mass balance.
- Estimate the and fields through a call to SMAGOR, if SMAG is true.
- Estimate the horizontal viscosity fluxes in equations (58) and (59) and store these through a call to HORVISCUVPOM if VISCPOM is true, and a call to HORVISCUV otherwise.
- Update the , , and fields taking into account the effects of the Coriolis terms, the horizontal viscosity terms, the , , and terms of equations (58) and (59).
- Call VERTVISCUVB to include the effects of vertical viscosity on the momentum.
- Include the effects of surface and bottom stress on the , , and fields.
- Compute the sigma coordinate vertical velocities from equation (36) using fluxes produced by UADV and VADV and the water levels ETA and ETAP.

If we are running experiments with river runoff, RIVLOG = .TRUE. , the effects of the runoff on the U, V, W and ETA fields are also included in modesplit.f90.

- SUBROUTINE SMAGOR2D
- SUBROUTINES HORVISCUV2D and HORVISCUV2DPOM
- SUBROUTINE SMAGOR
- SUBROUTINES HORVISCUV and HORVISCUVPOM
- SUBROUTINE STEPU3D and SUBROUTINE STEPV3D
- SUBROUTINE VERTVISCUVB