###

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,
(variable taus) and
(variable taub), respectively, are calculated and interpolated to the T-points.
Input parameters to this routine are rru and tt rrv, which
contain
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 - and -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