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

**INTERFACE:**

subroutine stresses_3d

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

Original author(s): Hans Burchard & Karsten Bolding

integer :: i,j,k,ku1,ku2,kv1,kv2 REALTYPE :: rho_0i

kklingbe 2017-10-02