15 |
|
|
16 |
\subsection{The finite volume method: finite volumes versus finite difference} |
\subsection{The finite volume method: finite volumes versus finite difference} |
17 |
\begin{rawhtml} |
\begin{rawhtml} |
18 |
<!-- CMIREDIR:finite_volume --> |
<!-- CMIREDIR:finite_volume: --> |
19 |
\end{rawhtml} |
\end{rawhtml} |
20 |
|
|
21 |
|
|
366 |
|
|
367 |
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 |
368 |
approach because the tracer points are at cell centers; the cell |
approach because the tracer points are at cell centers; the cell |
369 |
centers are mid-way between the cell interfaces. An alternative, the |
centers are mid-way between the cell interfaces. |
370 |
vertex or interface centered approach, is shown in |
This discretisation is selected when the thickness of the |
371 |
|
levels are provided ({\bf delR}, parameter file {\em data}, |
372 |
|
namelist {\em PARM04}) |
373 |
|
An alternative, the vertex or interface centered approach, is shown in |
374 |
Fig.~\ref{fig:vgrid}b. Here, the interior interfaces are positioned |
Fig.~\ref{fig:vgrid}b. Here, the interior interfaces are positioned |
375 |
mid-way between the tracer nodes (no longer cell centers). This |
mid-way between the tracer nodes (no longer cell centers). This |
376 |
approach is formally more accurate for evaluation of hydrostatic |
approach is formally more accurate for evaluation of hydrostatic |
377 |
pressure and vertical advection but historically the cell centered |
pressure and vertical advection but historically the cell centered |
378 |
approach has been used. An alternative form of subroutine {\em |
approach has been used. An alternative form of subroutine {\em |
379 |
INI\_VERTICAL\_GRID} is used to select the interface centered approach |
INI\_VERTICAL\_GRID} is used to select the interface centered approach |
380 |
but no run time option is currently available. |
This form requires to specify $Nr+1$ vertical distances {\bf delRc} |
381 |
|
(parameter file {\em data}, namelist {\em PARM04}, e.g. |
382 |
|
{\em verification/ideal\_2D\_oce/input/data}) |
383 |
|
corresponding to surface to center, $Nr-1$ center to center, and center to |
384 |
|
bottom distances. |
385 |
|
%but no run time option is currently available. |
386 |
|
|
387 |
\fbox{ \begin{minipage}{4.75in} |
\fbox{ \begin{minipage}{4.75in} |
388 |
{\em S/R INI\_VERTICAL\_GRID} ({\em |
{\em S/R INI\_VERTICAL\_GRID} ({\em |
400 |
|
|
401 |
|
|
402 |
\subsection{Topography: partially filled cells} |
\subsection{Topography: partially filled cells} |
403 |
|
\begin{rawhtml} |
404 |
|
<!-- CMIREDIR:topo_partial_cells: --> |
405 |
|
\end{rawhtml} |
406 |
|
|
407 |
\begin{figure} |
\begin{figure} |
408 |
\begin{center} |
\begin{center} |