/[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.4 by cnh,
Thu Oct 25 00:55:28 2001 UTC
+revision 1.12 by jmc,
Fri Jan 18 21:26:45 2008 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 145 
 useEnergyConservingCoriolis} to {\em tru
+Line 149 
 useEnergyConservingCoriolis} to {\em tru
  {\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} }
-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 257 
 represent the an-isotropic cosine scalin
+Line 261 
 represent the an-isotropic cosine scalin
  ``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 276 
 viscA4}), has units of $m^4 s^{-1}$.
+Line 280 
 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
- 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 307 
 G_v^{side-drag} & = &
+Line 311 
 G_v^{side-drag} & = &
  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.4
 


changed lines


 
Added in v.1.12
 Legend:



Removed from v.1.4
 


changed lines


 
Added in v.1.12
-Removed from v.1.4
+Added in v.1.12

	ViewVC Help
Powered by ViewVC 1.1.22