structure_friction_3d - (Source File: structure_friction_3d.F90)

INTERFACE:

    subroutine structure_friction_3d()

DESCRIPTION:

Here, the quadratic friction term resulting from a structure in the water column is calculated. This term will be added as additional forcing to the three-dimensional momentum equations, where it is treated numerically implicitly. Therefore here, only the following terms is calculated:

$$= C(z) \sqrt{u(z)^2+v(z)^2},$$ (134)

with the friction coefficient $C$ bearing the physical unit [1/m]. USES:

    use domain, only: imin,imax,jmin,jmax,kmax
    use variables_3d, only: uu,vv,sf,huo,hvo
    #define CALC_HALO_WW
 #ifndef CALC_HALO_WW
    use domain, only: az
    use halo_zones, only: update_3d_halo,wait_halo,z_TAG
 #endif
    use getm_timers, only: tic, toc, TIM_STRCTFRICT
    IMPLICIT NONE

REVISION HISTORY:

    Original author(s): Hans Burchard & Karsten Bolding

LOCAL VARIABLES:

    REALTYPE                  :: dtm1
    integer                   :: i,j,k
 #ifdef STRUCTURE_FRICTION
    REALTYPE                  :: cds(I2DFIELD)
 #endif