125 |
\section{Linear advection schemes} |
\section{Linear advection schemes} |
126 |
\label{sect:tracer-advection} |
\label{sect:tracer-advection} |
127 |
\begin{rawhtml} |
\begin{rawhtml} |
128 |
<!-- CMIREDIR:linear_advection_schemes --> |
<!-- CMIREDIR:linear_advection_schemes: --> |
129 |
\end{rawhtml} |
\end{rawhtml} |
130 |
|
|
131 |
\begin{figure} |
\begin{figure} |
353 |
|
|
354 |
|
|
355 |
\section{Non-linear advection schemes} |
\section{Non-linear advection schemes} |
356 |
|
\begin{rawhtml} |
357 |
|
<!-- CMIREDIR:non-linear_advection_schemes: --> |
358 |
|
\end{rawhtml} |
359 |
|
|
360 |
Non-linear advection schemes invoke non-linear interpolation and are |
Non-linear advection schemes invoke non-linear interpolation and are |
361 |
widely used in computational fluid dynamics (non-linear does not refer |
widely used in computational fluid dynamics (non-linear does not refer |
661 |
|
|
662 |
\section{Comparison of advection schemes} |
\section{Comparison of advection schemes} |
663 |
|
|
664 |
|
\begin{table}[htb] |
665 |
|
\centering |
666 |
|
\begin{tabular}[htb]{|l|c|c|c|c|l|} |
667 |
|
\hline |
668 |
|
Advection Scheme & code & use & use Multi- & Stencil & comments \\ |
669 |
|
& & A.B. & dimension & (1 dim) & \\ |
670 |
|
\hline \hline |
671 |
|
centered $2^{nd}$order & 2 & Yes & No & 3 pts & linear \\ |
672 |
|
\hline |
673 |
|
$3^{rd}$order upwind & 3 & Yes & No & 5 pts & linear/tracer\\ |
674 |
|
\hline |
675 |
|
centered $4^{th}$order & 4 & Yes & No & 5 pts & linear \\ |
676 |
|
\hline \hline |
677 |
|
% Lax-Wendroff & 10 & No & Yes & 3 pts & linear/tracer, non-linear/flow\\ |
678 |
|
% \hline |
679 |
|
$3^{rd}$order DST & 30 & No & Yes & 5 pts & linear/tracer, non-linear/flow\\ |
680 |
|
\hline \hline |
681 |
|
$2^{nd}$order Flux Limiters & 77 & No & Yes & 5 pts & non-linear \\ |
682 |
|
\hline |
683 |
|
$3^{nd}$order DST Flux limiter & 33 & No & Yes & 5 pts & non-linear \\ |
684 |
|
\hline |
685 |
|
\end{tabular} |
686 |
|
\caption{Summary of the different advection schemes available in MITgcm. |
687 |
|
``A.B.'' stands for Adams-Bashforth and ``DST'' for direct space time. |
688 |
|
The code corresponds to the number used to select the corresponding |
689 |
|
advection scheme in the parameter file (e.g., {\em tempAdvScheme=3} in |
690 |
|
file {\em data} selects the $3^{rd}$ order upwind advection scheme |
691 |
|
for temperature). |
692 |
|
} |
693 |
|
\label{tab:advectionShemes_summary} |
694 |
|
\end{table} |
695 |
|
|
696 |
|
|
697 |
Figs.~\ref{fig:advect-2d-lo-diag}, \ref{fig:advect-2d-mid-diag} and |
Figs.~\ref{fig:advect-2d-lo-diag}, \ref{fig:advect-2d-mid-diag} and |
698 |
\ref{fig:advect-2d-hi-diag} show solutions to a simple diagonal |
\ref{fig:advect-2d-hi-diag} show solutions to a simple diagonal |
699 |
advection problem using a selection of schemes for low, moderate and |
advection problem using a selection of schemes for low, moderate and |