subroutine structure_friction_3d()

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:

sf | (134) |

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

Original author(s): Hans Burchard & Karsten Bolding

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