For computational efficiency, the mode splitting technique (Blumberg and Mellor, 1987) is applied with a barotropic time step of 25 seconds, based on the Courant-Friederichs-Levy (1928) computational stability (CFL) condition and the external wave speed; and a baroclinic time step of 900 seconds, based on the CFL condition and the internal wave speed.