--- manual/s_algorithm/text/tracer.tex 2008/01/15 23:52:12 1.23 +++ manual/s_algorithm/text/tracer.tex 2016/10/26 17:30:23 1.27 @@ -1,8 +1,8 @@ -% $Header: /home/ubuntu/mnt/e9_copy/manual/s_algorithm/text/tracer.tex,v 1.23 2008/01/15 23:52:12 jmc Exp $ +% $Header: /home/ubuntu/mnt/e9_copy/manual/s_algorithm/text/tracer.tex,v 1.27 2016/10/26 17:30:23 heimbach Exp $ % $Name: $ \section{Tracer equations} -\label{sect:tracer_equations} +\label{sec:tracer_equations} \begin{rawhtml} \end{rawhtml} @@ -18,7 +18,7 @@ described here. \subsection{Time-stepping of tracers: ABII} -\label{sect:tracer_equations_abII} +\label{sec:tracer_equations_abII} \begin{rawhtml} \end{rawhtml} @@ -129,13 +129,13 @@ \section{Linear advection schemes} -\label{sect:tracer-advection} +\label{sec:tracer-advection} \begin{rawhtml} \end{rawhtml} \begin{figure} -\resizebox{5.5in}{!}{\includegraphics{part2/advect-1d-lo.eps}} +\resizebox{5.5in}{!}{\includegraphics{s_algorithm/figs/advect-1d-lo.eps}} \caption{ Comparison of 1-D advection schemes. Courant number is 0.05 with 60 points and solutions are shown for T=1 (one complete period). @@ -153,7 +153,7 @@ \end{figure} \begin{figure} -\resizebox{5.5in}{!}{\includegraphics{part2/advect-1d-hi.eps}} +\resizebox{5.5in}{!}{\includegraphics{s_algorithm/figs/advect-1d-hi.eps}} \caption{ Comparison of 1-D advection schemes. Courant number is 0.89 with 60 points and solutions are shown for T=1 (one complete period). @@ -359,6 +359,7 @@ \section{Non-linear advection schemes} +\label{sec:non-linear_advection_schemes} \begin{rawhtml} \end{rawhtml} @@ -549,7 +550,7 @@ \subsection{Multi-dimensional advection} \begin{figure} -\resizebox{5.5in}{!}{\includegraphics{part2/advect-2d-lo-diag.eps}} +\resizebox{5.5in}{!}{\includegraphics{s_algorithm/figs/advect-2d-lo-diag.eps}} \caption{ Comparison of advection schemes in two dimensions; diagonal advection of a resolved Gaussian feature. Courant number is 0.01 with @@ -570,7 +571,7 @@ \end{figure} \begin{figure} -\resizebox{5.5in}{!}{\includegraphics{part2/advect-2d-mid-diag.eps}} +\resizebox{5.5in}{!}{\includegraphics{s_algorithm/figs/advect-2d-mid-diag.eps}} \caption{ Comparison of advection schemes in two dimensions; diagonal advection of a resolved Gaussian feature. Courant number is 0.27 with @@ -591,7 +592,7 @@ \end{figure} \begin{figure} -\resizebox{5.5in}{!}{\includegraphics{part2/advect-2d-hi-diag.eps}} +\resizebox{5.5in}{!}{\includegraphics{s_algorithm/figs/advect-2d-hi-diag.eps}} \caption{ Comparison of advection schemes in two dimensions; diagonal advection of a resolved Gaussian feature. Courant number is 0.47 with @@ -665,20 +666,21 @@ \end{minipage} } \begin{figure} -\resizebox{3.5in}{!}{\includegraphics{part2/multiDim_CS.eps}} +\resizebox{3.5in}{!}{\includegraphics{s_algorithm/figs/multiDim_CS.eps}} \caption{Muti-dimensional advection time-stepping with Cubed-Sphere topology \label{fig:advect-multidim_cs} } \end{figure} \section{Comparison of advection schemes} -\label{sect:tracer_advection_schemes} +\label{sec:tracer_advection_schemes} \begin{rawhtml} \end{rawhtml} \begin{table}[htb] \centering +{\small \begin{tabular}[htb]{|l|c|c|c|c|l|} \hline Advection Scheme & code & use & use Multi- & Stencil & comments \\ @@ -696,12 +698,31 @@ No & Yes & 3 pts & linear/$\tau$, non-linear/v\\ \hline $3^{rd}$order DST & 30 & No & Yes & 5 pts & linear/$\tau$, non-linear/v\\ + \hline + $2^{nd}$order-moment Prather & 80 & No & Yes & ~ & ~ \\ \hline \hline $2^{nd}$order Flux Limiters & 77 & No & Yes & 5 pts & non-linear \\ \hline $3^{nd}$order DST Flux limiter & 33 & No & Yes & 5 pts & non-linear \\ \hline + $2^{nd}$order-moment Prather w. limiter & 81 & No & Yes & ~ & ~ \\ + \hline + piecewise parabolic w. ``null'' limiter & 40 & No & Yes & ~ & ~ \\ + \hline + piecewise parabolic w. ``mono'' limiter & 41 & No & Yes & ~ & ~ \\ + \hline + piecewise quartic w. ``null'' limiter & 50 & No & Yes & ~ & ~ \\ + \hline + piecewise quartic w. ``mono'' limiter & 51 & No & Yes & ~ & ~ \\ + \hline + piecewise quartic w. ``weno'' limiter & 52 & No & Yes & ~ & ~ \\ + \hline + $7^{nd}$order one-step method & 7 & No & Yes & ~ & ~ \\ + with Monotonicity Preserving Limiter & ~ & ~ & ~ & ~ & ~ \\ + \hline + \end{tabular} + } \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