stresses_3d - bottom and surface stresses (Source File: stresses_3d.F90)

INTERFACE:

    subroutine stresses_3d
DESCRIPTION:

As preparation of the call to do_turbulence in the routine gotm, see section 8.13.15, the normalised surface and bottom stresses, $\tau_s/\rho_0$ (variable taus) and $\tau_b/\rho_0$ (variable taub), respectively, are calculated and interpolated to the T-points. Input parameters to this routine are rru and tt rrv, which contain $r\sqrt{u^2+v^2}$ for the U- and V-points, respectively. The modules of the surface and bottom stress vectors are calculated then by means of taking the square root of the sum of the squares of the stess components. In a similar way also the $x$- and $y$-components of the bottom stress are computed for output. USES:

    use parameters, only: rho_0
    use domain, only: az,au,av,imin,imax,jmin,jmax
    use variables_3d, only: kumin,kvmin,uu,vv,hun,hvn,rru,rrv
    use variables_3d, only: taus,taubx,tauby,taub
    use meteo, only: tausx,tausy
    use halo_zones, only : update_2d_halo,wait_halo,z_TAG
    use getm_timers, only: tic, toc, TIM_STRESSES3D, TIM_STRESSES3DH
  $ use omp_lib
    IMPLICIT NONE
REVISION HISTORY:
    Original author(s): Hans Burchard & Karsten Bolding
LOCAL VARIABLES:
    integer                   :: i,j,k,ku1,ku2,kv1,kv2
    REALTYPE                  :: rho_0i