uv_depths.F90 3.66 KB
 kbk committed Mar 01, 2006 1 !$Id: uv_depths.F90,v 1.12 2006-03-01 15:54:07 kbk Exp$  gotm committed May 02, 2002 2 3 4 5 #include "cppdefs.h" !----------------------------------------------------------------------- !BOP !  hb committed Feb 04, 2006 6 ! !IROUTINE: uv_depths - calculate depths in u and v points.\label{sec-uv-depth}  gotm committed May 02, 2002 7 8 ! ! !INTERFACE:  kbk committed Jan 05, 2004 9  subroutine uv_depths(vel_depth_method)  gotm committed May 02, 2002 10 11 12 ! ! !DESCRIPTION: !  hb committed Feb 04, 2006 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 ! 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 {\tt vel\_depth\_method}: ! ! ! H^u_{i,j} = \left\{ ! \begin{array}{ll} ! \displaystyle ! \frac12 \left(H_{i,j}+H_{i+1,j}\right), & ! \displaystyle ! \mbox{ for {\tt vel\_depth\_method}} =0, \\ \\ ! \displaystyle ! \min\left\{H_{i,j}+H_{i+1,j}\right\}, & ! \displaystyle ! \mbox{ for {\tt vel\_depth\_method}} =1, \\ \\ ! \displaystyle ! \min\left\{H_{i,j}+H_{i+1,j}\right\}, & ! \displaystyle ! \mbox{ for {\tt vel\_depth\_method}} =2 \mbox{ and } \min\{H_{i,j}i,H_{i+1,j}\}