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:

$$H^u_{i,j} = \left\{ \begin{array}{ll} \displaystyle \frac12 \... ...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