/[MITgcm]/manual/s_algorithm/text/mom_fluxform.tex

Diff of /manual/s_algorithm/text/mom_fluxform.tex

Parent Directory | Revision Log | View Revision Graph Revision Graph | View Patch Patch

-revision 1.2 by adcroft,
Tue Sep 25 20:13:42 2001 UTC
+revision 1.13 by jmc,
Thu Feb 18 21:26:26 2010 UTC
 Line 2
  % $Name$
  \section{Flux-form momentum equations}
+ \label{sect:flux-form_momentum_equations}
+ \begin{rawhtml}
+ <!-- CMIREDIR:flux-form_momentum_eqautions: -->
+ \end{rawhtml}
  The original finite volume model was based on the Eulerian flux form
  momentum equations. This is the default though the vector invariant
-Line 25 
 In the hydrostatic limit, $G_w=0$ and $\
+Line 29 
 In the hydrostatic limit, $G_w=0$ and $\
  vertical momentum to hydrostatic balance.
  These terms are calculated in routines called from subroutine {\em
- CALC\_MOM\_RHS} a collected into the global arrays {\bf Gu}, {\bf Gv},
+ MOM\_FLUXFORM} a collected into the global arrays {\bf Gu}, {\bf Gv},
  and {\bf Gw}.
  \fbox{ \begin{minipage}{4.75in}
- {\em S/R CALC\_MOM\_RHS} ({\em pkg/mom\_fluxform/calc\_mom\_rhs.F})
+ {\em S/R MOM\_FLUXFORM} ({\em pkg/mom\_fluxform/mom\_fluxform.F})
  $G_u$: {\bf Gu} ({\em DYNVARS.h})
-Line 86 
 conserves kinetic energy.
+Line 90 
 conserves kinetic energy.
  {\em S/R MOM\_U\_ADV\_WV} ({\em mom\_u\_adv\_wv.F})
- $uu$, $uv$, $vu$, $vv$: {\bf aF} (local to {\em calc\_mom\_rhs.F})
+ $uu$, $uv$, $vu$, $vv$: {\bf aF} (local to {\em mom\_fluxform.F})
  \end{minipage} }
-Line 106 
 discretized:
+Line 110 
 discretized:
  \end{eqnarray}
  where the Coriolis parameters $f$ and $f'$ are defined:
  \begin{eqnarray}
- f & = & 2 \Omega \sin{\phi} \\
+ f & = & 2 \Omega \sin{\varphi} \\
- f' & = & 2 \Omega \cos{\phi}
+ f' & = & 2 \Omega \cos{\varphi}
  \end{eqnarray}
- where $\phi$ is geographic latitude when using spherical geometry,
+ where $\varphi$ is geographic latitude when using spherical geometry,
  otherwise the $\beta$-plane definition is used:
  \begin{eqnarray}
  f & = & f_o + \beta y \\
-Line 132 
 G_w^{Cor} & = &
+Line 136 
 G_w^{Cor} & = &
  \marginpar{Need to change the default in code to match this}
  where the subscripts on $f$ and $f'$ indicate evaluation of the
  Coriolis parameters at the appropriate points in space. The above
- discretization does {\em not} conserve anything, especially energy. An
+ discretization does {\em not} conserve anything, especially energy and
- option to recover this discretization has been retained for backward
+ for historical reasons is the default for the code. A
- compatibility testing (set run-time logical {\bf
+ flag controls this discretization: set run-time logical {\bf
- useNonconservingCoriolis} to {\em true} which otherwise defaults to
+ useEnergyConservingCoriolis} to {\em true} which otherwise defaults to
- {\em false}).
+ {\em false}.
  \fbox{ \begin{minipage}{4.75in}
  {\em S/R MOM\_CDSCHEME} ({\em mom\_cdscheme.F})
-Line 145 
 useNonconservingCoriolis} to {\em true}
+Line 149 
 useNonconservingCoriolis} to {\em true}
  {\em S/R MOM\_V\_CORIOLIS} ({\em mom\_v\_coriolis.F})
- $G_u^{Cor}$, $G_v^{Cor}$: {\bf cF} (local to {\em calc\_mom\_rhs.F})
+ $G_u^{Cor}$, $G_v^{Cor}$: {\bf cF} (local to {\em mom\_fluxform.F})
  \end{minipage} }
  \subsection{Curvature metric terms}
  The most commonly used coordinate system on the sphere is the
- geographic system $(\lambda,\phi)$. The curvilinear nature of these
+ geographic system $(\lambda,\varphi)$. The curvilinear nature of these
  coordinates on the sphere lead to some ``metric'' terms in the
  component momentum equations. Under the thin-atmosphere and
  hydrostatic approximations these terms are discretized:
  \begin{eqnarray}
  {\cal A}_w \Delta r_f h_w G_u^{metric} & = &
-   \overline{ \frac{ \overline{u}^i }{a} \tan{\phi} {\cal A}_c \Delta r_f h_c \overline{ v }^j }^i \\
+   \overline{ \frac{ \overline{u}^i }{a} \tan{\varphi} {\cal A}_c \Delta r_f h_c \overline{ v }^j }^i \\
  {\cal A}_s \Delta r_f h_s G_v^{metric} & = &
- - \overline{ \frac{ \overline{u}^i }{a} \tan{\phi} {\cal A}_c \Delta r_f h_c \overline{ u }^i }^j \\
+ - \overline{ \frac{ \overline{u}^i }{a} \tan{\varphi} {\cal A}_c \Delta r_f h_c \overline{ u }^i }^j \\
  G_w^{metric} & = & 0
  \end{eqnarray}
  where $a$ is the radius of the planet (sphericity is assumed) or the
  radial distance of the particle (i.e. a function of height).  It is
  easy to see that this discretization satisfies all the properties of
  the discrete Coriolis terms since the metric factor $\frac{u}{a}
- \tan{\phi}$ can be viewed as a modification of the vertical Coriolis
+ \tan{\varphi}$ can be viewed as a modification of the vertical Coriolis
- parameter: $f \rightarrow f+\frac{u}{a} \tan{\phi}$.
+ parameter: $f \rightarrow f+\frac{u}{a} \tan{\varphi}$.
  However, as for the Coriolis terms, a non-energy conserving form has
  exclusively been used to date:
  \begin{eqnarray}
- G_u^{metric} & = & \frac{u \overline{v}^{ij} }{a} \tan{\phi} \\
+ G_u^{metric} & = & \frac{u \overline{v}^{ij} }{a} \tan{\varphi} \\
- G_v^{metric} & = & \frac{ \overline{u}^{ij} \overline{u}^{ij}}{a} \tan{\phi}
+ G_v^{metric} & = & \frac{ \overline{u}^{ij} \overline{u}^{ij}}{a} \tan{\varphi}
  \end{eqnarray}
- where $\tan{\phi}$ is evaluated at the $u$ and $v$ points
+ where $\tan{\varphi}$ is evaluated at the $u$ and $v$ points
  respectively.
  \fbox{ \begin{minipage}{4.75in}
-Line 184 
 respectively.
+Line 188 
 respectively.
  {\em S/R MOM\_V\_METRIC\_SPHERE} ({\em mom\_v\_metric\_sphere.F})
- $G_u^{metric}$, $G_v^{metric}$: {\bf mT} (local to {\em calc\_mom\_rhs.F})
+ $G_u^{metric}$, $G_v^{metric}$: {\bf mT} (local to {\em mom\_fluxform.F})
  \end{minipage} }
-Line 193 
 $G_u^{metric}$, $G_v^{metric}$: {\bf mT}
+Line 197 
 $G_u^{metric}$, $G_v^{metric}$: {\bf mT}
  For the non-hydrostatic equations, dropping the thin-atmosphere
  approximation re-introduces metric terms involving $w$ and are
- required to conserve anglular momentum:
+ required to conserve angular momentum:
  \begin{eqnarray}
  {\cal A}_w \Delta r_f h_w G_u^{metric} & = &
  - \overline{ \frac{ \overline{u}^i \overline{w}^k }{a} {\cal A}_c \Delta r_f h_c }^i \\
-Line 219 
 G_w^{metric} & = &
+Line 223 
 G_w^{metric} & = &
  {\em S/R MOM\_V\_METRIC\_NH} ({\em mom\_v\_metric\_nh.F})
- $G_u^{metric}$, $G_v^{metric}$: {\bf mT} (local to {\em calc\_mom\_rhs.F})
+ $G_u^{metric}$, $G_v^{metric}$: {\bf mT} (local to {\em mom\_fluxform.F})
  \end{minipage} }
-Line 240 
 tensor such as symmetry.
+Line 244 
 tensor such as symmetry.
  The lateral viscous stresses are discretized:
  \begin{eqnarray}
- \tau_{11} & = & A_h c_{11\Delta}(\phi) \frac{1}{\Delta x_f} \delta_i u
+ \tau_{11} & = & A_h c_{11\Delta}(\varphi) \frac{1}{\Delta x_f} \delta_i u
-                -A_4 c_{11\Delta^2}(\phi) \frac{1}{\Delta x_f} \delta_i \nabla^2 u \\
+                -A_4 c_{11\Delta^2}(\varphi) \frac{1}{\Delta x_f} \delta_i \nabla^2 u \\
- \tau_{12} & = & A_h c_{12\Delta}(\phi) \frac{1}{\Delta y_u} \delta_j u
+ \tau_{12} & = & A_h c_{12\Delta}(\varphi) \frac{1}{\Delta y_u} \delta_j u
-                -A_4 c_{12\Delta^2}(\phi)\frac{1}{\Delta y_u} \delta_j \nabla^2 u \\
+                -A_4 c_{12\Delta^2}(\varphi)\frac{1}{\Delta y_u} \delta_j \nabla^2 u \\
- \tau_{21} & = & A_h c_{21\Delta}(\phi) \frac{1}{\Delta x_v} \delta_i v
+ \tau_{21} & = & A_h c_{21\Delta}(\varphi) \frac{1}{\Delta x_v} \delta_i v
-                -A_4 c_{21\Delta^2}(\phi) \frac{1}{\Delta x_v} \delta_i \nabla^2 v \\
+                -A_4 c_{21\Delta^2}(\varphi) \frac{1}{\Delta x_v} \delta_i \nabla^2 v \\
- \tau_{22} & = & A_h c_{22\Delta}(\phi) \frac{1}{\Delta y_f} \delta_j v
+ \tau_{22} & = & A_h c_{22\Delta}(\varphi) \frac{1}{\Delta y_f} \delta_j v
-                -A_4 c_{22\Delta^2}(\phi) \frac{1}{\Delta y_f} \delta_j \nabla^2 v
+                -A_4 c_{22\Delta^2}(\varphi) \frac{1}{\Delta y_f} \delta_j \nabla^2 v
  \end{eqnarray}
- where the non-dimensional factors $c_{lm\Delta^n}(\phi), \{l,m,n\} \in
+ where the non-dimensional factors $c_{lm\Delta^n}(\varphi), \{l,m,n\} \in
  \{1,2\}$ define the ``cosine'' scaling with latitude which can be
  applied in various ad-hoc ways. For instance, $c_{11\Delta} =
- c_{21\Delta} = (\cos{\phi})^{3/2}$, $c_{12\Delta}=c_{22\Delta}=0$ would
+ c_{21\Delta} = (\cos{\varphi})^{3/2}$, $c_{12\Delta}=c_{22\Delta}=1$ would
  represent the an-isotropic cosine scaling typically used on the
  ``lat-lon'' grid for Laplacian viscosity.
  \marginpar{Need to tidy up method for controlling this in code}
- It should be noted that dispite the ad-hoc nature of the scaling, some
+ It should be noted that despite the ad-hoc nature of the scaling, some
  scaling must be done since on a lat-lon grid the converging meridians
  make it very unlikely that a stable viscosity parameter exists across
  the entire model domain.
-Line 275 
 viscA4}), has units of $m^4 s^{-1}$.
+Line 279 
 viscA4}), has units of $m^4 s^{-1}$.
  {\em S/R MOM\_V\_YVISCFLUX} ({\em mom\_v\_yviscflux.F})
- $\tau_{11}$, $\tau_{12}$, $\tau_{22}$, $\tau_{22}$: {\bf vF}, {\bf
+ $\tau_{11}$, $\tau_{12}$, $\tau_{21}$, $\tau_{22}$: {\bf vF}, {\bf
- v4F} (local to {\em calc\_mom\_rhs.F})
+ v4F} (local to {\em mom\_fluxform.F})
  \end{minipage} }
  Two types of lateral boundary condition exist for the lateral viscous
-Line 290 
 handled using the lopped cells.
+Line 294 
 handled using the lopped cells.
  The no-slip condition defines the normal gradient of a tangential flow
  such that the flow is zero on the boundary. Rather than modify the
  stresses by using complicated functions of the masks and ``ghost''
- points (see \cite{Adcroft+Marshall98}) we add the boundary stresses as
+ points (see \cite{adcroft:98}) we add the boundary stresses as
  an additional source term in cells next to solid boundaries. This has
  the advantage of being able to cope with ``thin walls'' and also makes
  the interior stress calculation (code) independent of the boundary
-Line 298 
 conditions. The ``body'' force takes the
+Line 302 
 conditions. The ``body'' force takes the
  \begin{eqnarray}
  G_u^{side-drag} & = &
  \frac{4}{\Delta z_f} \overline{ (1-h_\zeta) \frac{\Delta x_v}{\Delta y_u} }^j
- \left( A_h c_{12\Delta}(\phi) u - A_4 c_{12\Delta^2}(\phi) \nabla^2 u \right)
+ \left( A_h c_{12\Delta}(\varphi) u - A_4 c_{12\Delta^2}(\varphi) \nabla^2 u \right)
  \\
  G_v^{side-drag} & = &
  \frac{4}{\Delta z_f} \overline{ (1-h_\zeta) \frac{\Delta y_u}{\Delta x_v} }^i
- \left( A_h c_{21\Delta}(\phi) v - A_4 c_{21\Delta^2}(\phi) \nabla^2 v \right)
+ \left( A_h c_{21\Delta}(\varphi) v - A_4 c_{21\Delta^2}(\varphi) \nabla^2 v \right)
  \end{eqnarray}
  In fact, the above discretization is not quite complete because it
  assumes that the bathymetry at velocity points is deeper than at
- neighbouring vorticity points, e.g. $1-h_w < 1-h_\zeta$
+ neighboring vorticity points, e.g. $1-h_w < 1-h_\zeta$
  \fbox{ \begin{minipage}{4.75in}
  {\em S/R MOM\_U\_SIDEDRAG} ({\em mom\_u\_sidedrag.F})
  {\em S/R MOM\_V\_SIDEDRAG} ({\em mom\_v\_sidedrag.F})
- $G_u^{side-drag}$, $G_v^{side-drag}$: {\bf vF} (local to {\em calc\_mom\_rhs.F})
+ $G_u^{side-drag}$, $G_v^{side-drag}$: {\bf vF} (local to {\em mom\_fluxform.F})
  \end{minipage} }
-Line 324 
 Vertical viscosity terms are discretized
+Line 328 
 Vertical viscosity terms are discretized
  to the variable grid lengths introduced by the finite volume
  formulation. This reduces the formal accuracy of these terms to just
  first order but only next to boundaries; exactly where other terms
- appear such as linar and quadratic bottom drag.
+ appear such as linear and quadratic bottom drag.
  \begin{eqnarray}
  G_u^{v-diss} & = &
  \frac{1}{\Delta r_f h_w} \delta_k \tau_{13} \\
-Line 342 
 In the interior the vertical stresses ar
+Line 346 
 In the interior the vertical stresses ar
  \tau_{33} & = & A_v \frac{1}{\Delta r_f} \delta_k w
  \end{eqnarray}
  It should be noted that in the non-hydrostatic form, the stress tensor
- is even less consistent than for the hydrostatic (see Wazjowicz
+ is even less consistent than for the hydrostatic (see
- \cite{Waojz}). It is well known how to do this properly (see Griffies
+ \cite{wajsowicz:93}). It is well known how to do this properly (see
- \cite{Griffies}) and is on the list of to-do's.
+ \cite{griffies:00}) and is on the list of to-do's.
  \fbox{ \begin{minipage}{4.75in}
  {\em S/R MOM\_U\_RVISCLFUX} ({\em mom\_u\_rviscflux.F})
  {\em S/R MOM\_V\_RVISCLFUX} ({\em mom\_v\_rviscflux.F})
- $\tau_{13}$: {\bf urf} (local to {\em calc\_mom\_rhs.F})
+ $\tau_{13}$: {\bf urf} (local to {\em mom\_fluxform.F})
- $\tau_{23}$: {\bf vrf} (local to {\em calc\_mom\_rhs.F})
+ $\tau_{23}$: {\bf vrf} (local to {\em mom\_fluxform.F})
  \end{minipage} }
-Line 389 
 dimensionless with typical values in the
+Line 393 
 dimensionless with typical values in the
  {\em S/R MOM\_V\_BOTTOMDRAG} ({\em mom\_v\_bottomdrag.F})
- $\tau_{13}^{bottom-drag}$, $\tau_{23}^{bottom-drag}$: {\bf vf} (local to {\em calc\_mom\_rhs.F})
+ $\tau_{13}^{bottom-drag}/\Delta r_f$, $\tau_{23}^{bottom-drag}/\Delta r_f$:
+ {\bf vf} (local to {\em mom\_fluxform.F})
  \end{minipage} }
  \subsection{Derivation of discrete energy conservation}
-Line 401 
 KE = \frac{1}{2} \left( \overline{ u^2 }
+Line 406 
 KE = \frac{1}{2} \left( \overline{ u^2 }
  \epsilon_{nh} \overline{ w^2 }^k \right)
  \end{equation}
+ \subsection{Mom Diagnostics}
+ \label{sec:pkg:mom_common:diagnostics}
+ \begin{verbatim}
+ ------------------------------------------------------------------------
+ <-Name->|Levs|<-parsing code->|<--  Units   -->|<- Tile (max=80c)
+ ------------------------------------------------------------------------
+ VISCAHZ | 15 |SZ      MR      |m^2/s           |Harmonic Visc Coefficient (m2/s) (Zeta Pt)
+ VISCA4Z | 15 |SZ      MR      |m^4/s           |Biharmonic Visc Coefficient (m4/s) (Zeta Pt)
+ VISCAHD | 15 |SM      MR      |m^2/s           |Harmonic Viscosity Coefficient (m2/s) (Div Pt)
+ VISCA4D | 15 |SM      MR      |m^4/s           |Biharmonic Viscosity Coefficient (m4/s) (Div Pt)
+ VAHZMAX | 15 |SZ      MR      |m^2/s           |CFL-MAX Harm Visc Coefficient (m2/s) (Zeta Pt)
+ VA4ZMAX | 15 |SZ      MR      |m^4/s           |CFL-MAX Biharm Visc Coefficient (m4/s) (Zeta Pt)
+ VAHDMAX | 15 |SM      MR      |m^2/s           |CFL-MAX Harm Visc Coefficient (m2/s) (Div Pt)
+ VA4DMAX | 15 |SM      MR      |m^4/s           |CFL-MAX Biharm Visc Coefficient (m4/s) (Div Pt)
+ VAHZMIN | 15 |SZ      MR      |m^2/s           |RE-MIN Harm Visc Coefficient (m2/s) (Zeta Pt)
+ VA4ZMIN | 15 |SZ      MR      |m^4/s           |RE-MIN Biharm Visc Coefficient (m4/s) (Zeta Pt)
+ VAHDMIN | 15 |SM      MR      |m^2/s           |RE-MIN Harm Visc Coefficient (m2/s) (Div Pt)
+ VA4DMIN | 15 |SM      MR      |m^4/s           |RE-MIN Biharm Visc Coefficient (m4/s) (Div Pt)
+ VAHZLTH | 15 |SZ      MR      |m^2/s           |Leith Harm Visc Coefficient (m2/s) (Zeta Pt)
+ VA4ZLTH | 15 |SZ      MR      |m^4/s           |Leith Biharm Visc Coefficient (m4/s) (Zeta Pt)
+ VAHDLTH | 15 |SM      MR      |m^2/s           |Leith Harm Visc Coefficient (m2/s) (Div Pt)
+ VA4DLTH | 15 |SM      MR      |m^4/s           |Leith Biharm Visc Coefficient (m4/s) (Div Pt)
+ VAHZLTHD| 15 |SZ      MR      |m^2/s           |LeithD Harm Visc Coefficient (m2/s) (Zeta Pt)
+ VA4ZLTHD| 15 |SZ      MR      |m^4/s           |LeithD Biharm Visc Coefficient (m4/s) (Zeta Pt)
+ VAHDLTHD| 15 |SM      MR      |m^2/s           |LeithD Harm Visc Coefficient (m2/s) (Div Pt)
+ VA4DLTHD| 15 |SM      MR      |m^4/s           |LeithD Biharm Visc Coefficient (m4/s) (Div Pt)
+ VAHZSMAG| 15 |SZ      MR      |m^2/s           |Smagorinsky Harm Visc Coefficient (m2/s) (Zeta Pt)
+ VA4ZSMAG| 15 |SZ      MR      |m^4/s           |Smagorinsky Biharm Visc Coeff. (m4/s) (Zeta Pt)
+ VAHDSMAG| 15 |SM      MR      |m^2/s           |Smagorinsky Harm Visc Coefficient (m2/s) (Div Pt)
+ VA4DSMAG| 15 |SM      MR      |m^4/s           |Smagorinsky Biharm Visc Coeff. (m4/s) (Div Pt)
+ momKE   | 15 |SM      MR      |m^2/s^2         |Kinetic Energy (in momentum Eq.)
+ momHDiv | 15 |SM      MR      |s^-1            |Horizontal Divergence (in momentum Eq.)
+ momVort3| 15 |SZ      MR      |s^-1            |3rd component (vertical) of Vorticity
+ Strain  | 15 |SZ      MR      |s^-1            |Horizontal Strain of Horizontal Velocities
+ Tension | 15 |SM      MR      |s^-1            |Horizontal Tension of Horizontal Velocities
+ UBotDrag| 15 |UU   129MR      |m/s^2           |U momentum tendency from Bottom Drag
+ VBotDrag| 15 |VV   128MR      |m/s^2           |V momentum tendency from Bottom Drag
+ USidDrag| 15 |UU   131MR      |m/s^2           |U momentum tendency from Side Drag
+ VSidDrag| 15 |VV   130MR      |m/s^2           |V momentum tendency from Side Drag
+ Um_Diss | 15 |UU   133MR      |m/s^2           |U momentum tendency from Dissipation
+ Vm_Diss | 15 |VV   132MR      |m/s^2           |V momentum tendency from Dissipation
+ Um_Advec| 15 |UU   135MR      |m/s^2           |U momentum tendency from Advection terms
+ Vm_Advec| 15 |VV   134MR      |m/s^2           |V momentum tendency from Advection terms
+ Um_Cori | 15 |UU   137MR      |m/s^2           |U momentum tendency from Coriolis term
+ Vm_Cori | 15 |VV   136MR      |m/s^2           |V momentum tendency from Coriolis term
+ Um_Ext  | 15 |UU   137MR      |m/s^2           |U momentum tendency from external forcing
+ Vm_Ext  | 15 |VV   138MR      |m/s^2           |V momentum tendency from external forcing
+ Um_AdvZ3| 15 |UU   141MR      |m/s^2           |U momentum tendency from Vorticity Advection
+ Vm_AdvZ3| 15 |VV   140MR      |m/s^2           |V momentum tendency from Vorticity Advection
+ Um_AdvRe| 15 |UU   143MR      |m/s^2           |U momentum tendency from vertical Advection (Explicit part)
+ Vm_AdvRe| 15 |VV   142MR      |m/s^2           |V momentum tendency from vertical Advection (Explicit part)
+ ADVx_Um | 15 |UM   145MR      |m^4/s^2         |Zonal      Advective Flux of U momentum
+ ADVy_Um | 15 |VZ   144MR      |m^4/s^2         |Meridional Advective Flux of U momentum
+ ADVrE_Um| 15 |WU      LR      |m^4/s^2         |Vertical   Advective Flux of U momentum (Explicit part)
+ ADVx_Vm | 15 |UZ   148MR      |m^4/s^2         |Zonal      Advective Flux of V momentum
+ ADVy_Vm | 15 |VM   147MR      |m^4/s^2         |Meridional Advective Flux of V momentum
+ ADVrE_Vm| 15 |WV      LR      |m^4/s^2         |Vertical   Advective Flux of V momentum (Explicit part)
+ VISCx_Um| 15 |UM   151MR      |m^4/s^2         |Zonal      Viscous Flux of U momentum
+ VISCy_Um| 15 |VZ   150MR      |m^4/s^2         |Meridional Viscous Flux of U momentum
+ VISrE_Um| 15 |WU      LR      |m^4/s^2         |Vertical   Viscous Flux of U momentum (Explicit part)
+ VISrI_Um| 15 |WU      LR      |m^4/s^2         |Vertical   Viscous Flux of U momentum (Implicit part)
+ VISCx_Vm| 15 |UZ   155MR      |m^4/s^2         |Zonal      Viscous Flux of V momentum
+ VISCy_Vm| 15 |VM   154MR      |m^4/s^2         |Meridional Viscous Flux of V momentum
+ VISrE_Vm| 15 |WV      LR      |m^4/s^2         |Vertical   Viscous Flux of V momentum (Explicit part)
+ VISrI_Vm| 15 |WV      LR      |m^4/s^2         |Vertical   Viscous Flux of V momentum (Implicit part)
+ \end{verbatim}

 Legend:



Removed from v.1.2
 


changed lines


 
Added in v.1.13
 Legend:



Removed from v.1.2
 


changed lines


 
Added in v.1.13
-Removed from v.1.2
+Added in v.1.13

	ViewVC Help
Powered by ViewVC 1.1.22