SUBROUTINE SUPERBEEUV

The subroutine SUPERBEEUV is called from the main program after the call to INTERNAL. SUPERBEEUV estimates contributions to $\Delta U$ and $\Delta V$ due to the terms in equations (56) and (57) describing advection of momentum and add these to the contributions estimated by INTERNAL. SUPERBEEUV applies an combination of a 1th orher upwind method and a 2nd order Lax-Wendroff method to compute the advective fluxes. The method is implemented following the description in Yang and Przekwas [46].

If RIVLOG is .TRUE. , fluxes of the momentum from rivers are also included.