slow_bottom_friction - slow bed friction (Source File: slow_bottom_friction.F90)

INTERFACE:

    subroutine slow_bottom_friction
DESCRIPTION:

This routine basically calculates the bed friction, as it would come out if the vertically and macro timestep averaged velocity would be used. The output of this subroutine is thus $R\sqrt{u^2+v^2}$ on the U-points (see variable ruu) and on the V-points (see rvv) with the vertically and macro timestep averaged velocity components on the old time step, $u$ and $v$, which are in the code denoted by Ui and Vi, respectively. The drag coefficient $R$ is given by eq. (71) on page [*]. The results for the variables ruu and rvv will then be used in the routine slow_terms described on page [*] for the calculation of the slow terms $S^x_F$ and $S^y_F$, see section 7.1.

USES:

    use parameters, only: kappa
    use domain, only: imin,imax,jmin,jmax,HU,HV,min_depth,au,av
    use variables_2d, only: zub,zvb,ru,rv,Uinto,Vinto
    use variables_3d, only: ssuo,ssun,ssvo,ssvn
    use getm_timers, only: tic, toc, TIM_SLOWBFRICT
    use exceptions, only: getm_error
  $ use omp_lib
    IMPLICIT NONE
REVISION HISTORY:
    Original author(s): Hans Burchard & Karsten Bolding
LOCAL VARIABLES:
    integer                   :: i,j
    REALTYPE                  :: uloc,vloc,HH
    logical,save              :: first=.true.
    REALTYPE                  :: Ui(I2DFIELD)
    REALTYPE                  :: Vi(I2DFIELD)
    REALTYPE                  :: ruu(I2DFIELD)
    REALTYPE                  :: rvv(I2DFIELD)