--- manual/s_algorithm/text/tracer.tex	2004/10/14 19:53:04	1.17
+++ manual/s_algorithm/text/tracer.tex	2006/06/28 16:55:52	1.21
@@ -1,8 +1,11 @@
-% $Header: /home/ubuntu/mnt/e9_copy/manual/s_algorithm/text/tracer.tex,v 1.17 2004/10/14 19:53:04 jmc Exp $
+% $Header: /home/ubuntu/mnt/e9_copy/manual/s_algorithm/text/tracer.tex,v 1.21 2006/06/28 16:55:52 jmc Exp $
 % $Name:  $
 
 \section{Tracer equations}
 \label{sect:tracer_equations}
+\begin{rawhtml}
+<!-- CMIREDIR:tracer_equations: -->
+\end{rawhtml}
 
 The basic discretization used for the tracer equations is the second
 order piece-wise constant finite volume form of the forced
@@ -16,6 +19,9 @@
 
 \subsection{Time-stepping of tracers: ABII}
 \label{sect:tracer_equations_abII}
+\begin{rawhtml}
+<!-- CMIREDIR:tracer_equations_abII: -->
+\end{rawhtml}
 
 The default advection scheme is the centered second order method which
 requires a second order or quasi-second order time-stepping scheme to
@@ -621,10 +627,10 @@
 \tau^{n+1/3} & = & \tau^{n}
 - \Delta t \left( \frac{1}{\Delta x} \delta_i F^x(\tau^{n})
            + \tau^{n} \frac{1}{\Delta x} \delta_i u \right) \\
-\tau^{n+2/3} & = & \tau^{n}
+\tau^{n+2/3} & = & \tau^{n+1/3}
 - \Delta t \left( \frac{1}{\Delta y} \delta_j F^y(\tau^{n+1/3})
            + \tau^{n} \frac{1}{\Delta y} \delta_i v \right) \\
-\tau^{n+3/3} & = & \tau^{n}
+\tau^{n+3/3} & = & \tau^{n+2/3}
 - \Delta t \left( \frac{1}{\Delta r} \delta_k F^x(\tau^{n+2/3})
            + \tau^{n} \frac{1}{\Delta r} \delta_i w \right)
 \end{eqnarray}
@@ -660,6 +666,10 @@
 
 
 \section{Comparison of advection schemes}
+\label{sect:tracer_advection_schemes}
+\begin{rawhtml}
+<!-- CMIREDIR:comparison_of_advection_schemes: -->
+\end{rawhtml}
 
 \begin{table}[htb]
 \centering
@@ -668,15 +678,18 @@
    Advection Scheme & code & use  & use Multi- & Stencil & comments \\
                     &      & A.B. & dimension & (1 dim) & \\
    \hline \hline
+   $1^{rst}$order upwind  & 1 &  No & Yes & 3 pts & linear/$\tau$, non-linear/v\\
+   \hline
    centered $2^{nd}$order & 2 &  Yes & No & 3 pts & linear \\
    \hline
-   $3^{rd}$order upwind   & 3 &  Yes & No & 5 pts & linear/tracer\\
+   $3^{rd}$order upwind   & 3 &  Yes & No & 5 pts & linear/$\tau$\\
    \hline
    centered $4^{th}$order & 4 &  Yes & No & 5 pts & linear \\
    \hline \hline
-%  Lax-Wendroff       & 10 &  No & Yes & 3 pts & linear/tracer, non-linear/flow\\
-%  \hline
-   $3^{rd}$order DST & 30 &  No & Yes & 5 pts & linear/tracer, non-linear/flow\\
+   $2^{nd}$order DST (Lax-Wendroff)  & 20 &
+                         No & Yes & 3 pts & linear/tracer, non-linear/flow\\
+   \hline
+   $3^{rd}$order DST & 30 &  No & Yes & 5 pts & linear/$\tau$, non-linear/v\\
    \hline \hline
    $2^{nd}$order Flux Limiters & 77 &  No & Yes & 5 pts & non-linear \\
    \hline
@@ -686,7 +699,7 @@
  \caption{Summary of the different advection schemes available in MITgcm.
           ``A.B.'' stands for Adams-Bashforth and ``DST'' for direct space time.
           The code corresponds to the number used to select the corresponding
-          advection scheme in the parameter file (e.g., {\em tempAdvScheme=3} in
+          advection scheme in the parameter file (e.g., {\bf tempAdvScheme}=3 in
           file {\em data} selects the $3^{rd}$ order upwind advection scheme 
           for temperature).
    }
