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:

\begin{displaymath}H^u_{i,j} = \left\{
\begin{array}{ll}
\displaystyle
\frac12 \...
...min\{H_{i,j},H_{i+1,j}\}\geq D_{crit} \\ \\
\end{array}\right.\end{displaymath} (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