--- manual/s_algorithm/text/spatial-discrete.tex 2004/10/16 03:40:12 1.17 +++ manual/s_algorithm/text/spatial-discrete.tex 2008/01/22 20:50:41 1.22 @@ -1,4 +1,4 @@ -% $Header: /home/ubuntu/mnt/e9_copy/manual/s_algorithm/text/spatial-discrete.tex,v 1.17 2004/10/16 03:40:12 edhill Exp $ +% $Header: /home/ubuntu/mnt/e9_copy/manual/s_algorithm/text/spatial-discrete.tex,v 1.22 2008/01/22 20:50:41 heimbach Exp $ % $Name: $ \section{Spatial discretization of the dynamical equations} @@ -13,8 +13,6 @@ representation of the position of the boundary. We treat the horizontal and vertical directions as separable and differently. -\input{part2/notation} - \subsection{The finite volume method: finite volumes versus finite difference} \begin{rawhtml} @@ -141,8 +139,8 @@ grid for all panels. a) The area of a tracer cell, $A_c$, is bordered by the lengths $\Delta x_g$ and $\Delta y_g$. b) The area of a vorticity cell, $A_\zeta$, is bordered by the lengths $\Delta x_c$ and -$\Delta y_c$. c) The area of a u cell, $A_c$, is bordered by the -lengths $\Delta x_v$ and $\Delta y_f$. d) The area of a v cell, $A_c$, +$\Delta y_c$. c) The area of a u cell, $A_w$, is bordered by the +lengths $\Delta x_v$ and $\Delta y_f$. d) The area of a v cell, $A_s$, is bordered by the lengths $\Delta x_f$ and $\Delta y_u$.} \label{fig:hgrid} \end{figure} @@ -150,8 +148,8 @@ The model domain is decomposed into tiles and within each tile a quasi-regular grid is used. A tile is the basic unit of domain decomposition for parallelization but may be used whether parallelized -or not; see section \ref{sect:tiles} for more details. Although the -tiles may be patched together in an unstructured manner +or not; see section \ref{sect:domain_decomposition} for more details. +Although the tiles may be patched together in an unstructured manner (i.e. irregular or non-tessilating pattern), the interior of tiles is a structured grid of quadrilateral cells. The horizontal coordinate system is orthogonal curvilinear meaning we can not necessarily treat @@ -188,7 +186,7 @@ Fig.~\ref{fig:hgrid}b shows the vorticity cell. The length of the southern edge, $\Delta x_c$, western edge, $\Delta y_c$ and surface area, $A_\zeta$, presented in the vertical are stored in arrays {\bf -DXg}, {\bf DYg} and {\bf rAz}. +DXc}, {\bf DYc} and {\bf rAz}. \marginpar{$A_\zeta$: {\bf rAz}} \marginpar{$\Delta x_c$: {\bf DXc}} \marginpar{$\Delta y_c$: {\bf DYc}} @@ -370,7 +368,7 @@ The above grid (Fig.~\ref{fig:vgrid}a) is known as the cell centered approach because the tracer points are at cell centers; the cell centers are mid-way between the cell interfaces. -This discretisation is selected when the thickness of the +This discretization is selected when the thickness of the levels are provided ({\bf delR}, parameter file {\em data}, namelist {\em PARM04}) An alternative, the vertex or interface centered approach, is shown in