uv_depths - calculate depths in u and v points.

INTERFACE:

    subroutine uv_depths(vel_depth_method)
DESCRIPTION:

In this routine which is called once during the model initialisation, the bathymetry value in the U- and the V-points are calculated from the bathymetry values in the T-points. The interpolation depends on the value which is given to vel_depth_method:

$\displaystyle H^u_{i,j} = \left\{ \begin{array}{ll} \displaystyle \frac12 \left...
...2 \mbox{ and } \min\{H_{i,j},H_{i+1,j}\}\geq D_{crit} \ \ \end{array} \right.$ (56)

The calculation of $ H^v_{i,j}$ is done accordingly.

The options 1 and 2 for vel_depth_method may help to stabilise calculations when drying and flooding is involved. USES:

    use exceptions
    use domain, only: imin,imax,jmin,jmax,az,au,av,H,HU,HV
    use getm_timers, only: tic,toc,TIM_UVDEPTHS
    IMPLICIT NONE
INPUT PARAMETERS:
    integer, intent(in)                 :: vel_depth_method
REVISION HISTORY:
    Original author(s): Hans Burchard & Karsten Bolding
LOCAL VARIABLES:
    integer                   :: i,j
    REALTYPE                  :: d_crit=2.0



kklingbe 2017-10-02