2 |
% $Name$ |
% $Name$ |
3 |
|
|
4 |
\section{Spatial discretization of the dynamical equations} |
\section{Spatial discretization of the dynamical equations} |
5 |
|
\begin{rawhtml} |
6 |
|
<!-- CMIREDIR:spatial_discretization_of_dyn_eq: --> |
7 |
|
\end{rawhtml} |
8 |
|
|
9 |
Spatial discretization is carried out using the finite volume |
Spatial discretization is carried out using the finite volume |
10 |
method. This amounts to a grid-point method (namely second-order |
method. This amounts to a grid-point method (namely second-order |
150 |
The model domain is decomposed into tiles and within each tile a |
The model domain is decomposed into tiles and within each tile a |
151 |
quasi-regular grid is used. A tile is the basic unit of domain |
quasi-regular grid is used. A tile is the basic unit of domain |
152 |
decomposition for parallelization but may be used whether parallelized |
decomposition for parallelization but may be used whether parallelized |
153 |
or not; see section \ref{sect:tiles} for more details. Although the |
or not; see section \ref{sect:domain_decomposition} for more details. |
154 |
tiles may be patched together in an unstructured manner |
Although the tiles may be patched together in an unstructured manner |
155 |
(i.e. irregular or non-tessilating pattern), the interior of tiles is |
(i.e. irregular or non-tessilating pattern), the interior of tiles is |
156 |
a structured grid of quadrilateral cells. The horizontal coordinate |
a structured grid of quadrilateral cells. The horizontal coordinate |
157 |
system is orthogonal curvilinear meaning we can not necessarily treat |
system is orthogonal curvilinear meaning we can not necessarily treat |
370 |
The above grid (Fig.~\ref{fig:vgrid}a) is known as the cell centered |
The above grid (Fig.~\ref{fig:vgrid}a) is known as the cell centered |
371 |
approach because the tracer points are at cell centers; the cell |
approach because the tracer points are at cell centers; the cell |
372 |
centers are mid-way between the cell interfaces. |
centers are mid-way between the cell interfaces. |
373 |
This discretisation is selected when the thickness of the |
This discretization is selected when the thickness of the |
374 |
levels are provided ({\bf delR}, parameter file {\em data}, |
levels are provided ({\bf delR}, parameter file {\em data}, |
375 |
namelist {\em PARM04}) |
namelist {\em PARM04}) |
376 |
An alternative, the vertex or interface centered approach, is shown in |
An alternative, the vertex or interface centered approach, is shown in |
476 |
|
|
477 |
|
|
478 |
\section{Continuity and horizontal pressure gradient terms} |
\section{Continuity and horizontal pressure gradient terms} |
479 |
|
\begin{rawhtml} |
480 |
|
<!-- CMIREDIR:continuity_and_horizontal_pressure: --> |
481 |
|
\end{rawhtml} |
482 |
|
|
483 |
|
|
484 |
The core algorithm is based on the ``C grid'' discretization of the |
The core algorithm is based on the ``C grid'' discretization of the |
485 |
continuity equation which can be summarized as: |
continuity equation which can be summarized as: |
518 |
evaporation and only enters the top-level of the {\em ocean} model. |
evaporation and only enters the top-level of the {\em ocean} model. |
519 |
|
|
520 |
\section{Hydrostatic balance} |
\section{Hydrostatic balance} |
521 |
|
\begin{rawhtml} |
522 |
|
<!-- CMIREDIR:hydrostatic_balance: --> |
523 |
|
\end{rawhtml} |
524 |
|
|
525 |
The vertical momentum equation has the hydrostatic or |
The vertical momentum equation has the hydrostatic or |
526 |
quasi-hydrostatic balance on the right hand side. This discretization |
quasi-hydrostatic balance on the right hand side. This discretization |