s_algorithm/text/tracer.tex

% $Header: $
% $Name: $

\section{Tracer equations}

The tracer equations are discretized consistantly with the continuity
equation to facilitate conservation properties analogous to the
continuum:
\begin{equation}
{\cal A}_c \Delta r_f h_c \partial_\theta
+ \delta_i U \overline{ \theta }^i
+ \delta_j V \overline{ \theta }^j
+ \delta_k W \overline{ \theta }^k
= {\cal A}_c \Delta r_f h_c {\cal S}_\theta + \theta {\cal A}_c \delta_k (P-E)_{r=0}
\end{equation}
The quantities $U$, $V$ and $W$ are volume fluxes defined:
\marginpar{$U$: {\bf uTrans} }
\marginpar{$V$: {\bf vTrans} }
\marginpar{$W$: {\bf rTrans} }
\begin{eqnarray}
U & = & \Delta y_g \Delta r_f h_w u \\
V & = & \Delta x_g \Delta r_f h_s v \\
W & = & {\cal A}_c w
\end{eqnarray}
${\cal S}$ represents the ``parameterized'' SGS processes and
physics associated with the tracer. For instance, potential
temperature equation in the ocean has is forced by surface and
partially penetrating heat fluxes:
\begin{equation}
{\cal A}_c \Delta r_f h_c {\cal S}_\theta = \frac{1}{c_p \rho_o} \delta_k {\cal A}_c {\cal Q}
\end{equation}
while the salt equation has no real sources, ${\cal S}=0$, which
leaves just the $P-E$ term.

The continuity equation can be recovered by setting ${\cal Q}=0$ and
$\theta=1$. The term $\theta (P-E)_{r=0}$ is required to retain local
conservation of $\theta$. Global conservation is not possible using
the flux-form (as here) and a linearized free-surface
(\cite{Griffies00,Campin02}).


1	% $Header: $
2	% $Name: $
3
4	\section{Tracer equations}
5
6	The tracer equations are discretized consistantly with the continuity
7	equation to facilitate conservation properties analogous to the
8	continuum:
9	\begin{equation}
10	{\cal A}_c \Delta r_f h_c \partial_\theta
11	+ \delta_i U \overline{ \theta }^i
12	+ \delta_j V \overline{ \theta }^j
13	+ \delta_k W \overline{ \theta }^k
14	= {\cal A}_c \Delta r_f h_c {\cal S}_\theta + \theta {\cal A}_c \delta_k (P-E)_{r=0}
15	\end{equation}
16	The quantities $U$, $V$ and $W$ are volume fluxes defined:
17	\marginpar{$U$: {\bf uTrans} }
18	\marginpar{$V$: {\bf vTrans} }
19	\marginpar{$W$: {\bf rTrans} }
20	\begin{eqnarray}
21	U & = & \Delta y_g \Delta r_f h_w u \\
22	V & = & \Delta x_g \Delta r_f h_s v \\
23	W & = & {\cal A}_c w
24	\end{eqnarray}
25	${\cal S}$ represents the ``parameterized'' SGS processes and
26	physics associated with the tracer. For instance, potential
27	temperature equation in the ocean has is forced by surface and
28	partially penetrating heat fluxes:
29	\begin{equation}
30	{\cal A}_c \Delta r_f h_c {\cal S}_\theta = \frac{1}{c_p \rho_o} \delta_k {\cal A}_c {\cal Q}
31	\end{equation}
32	while the salt equation has no real sources, ${\cal S}=0$, which
33	leaves just the $P-E$ term.
34
35	The continuity equation can be recovered by setting ${\cal Q}=0$ and
36	$\theta=1$. The term $\theta (P-E)_{r=0}$ is required to retain local
37	conservation of $\theta$. Global conservation is not possible using
38	the flux-form (as here) and a linearized free-surface
39	(\cite{Griffies00,Campin02}).
40
41
42
43