Home Contact Us Site Map  
 
       
    next up previous contents
Next: 2.16.4 First order upwind Up: 2.16 Linear advection schemes Previous: 2.16.2 Third order upwind   Contents

2.16.3 Centered fourth order advection

Centered fourth order advection is formally the most accurate scheme we have implemented and can be used to great effect in high resolution simulation where dynamical scales are well resolved. However, the scheme is noisy like the centered second order method and so must be used with some finite amount of diffusion. Bi-harmonic is recommended since it is more scale selective and less likely to diffuse away the well resolved gradient the fourth order scheme worked so hard to create.

The centered fourth order fluxes are discretized:

$\displaystyle F_x$ $\displaystyle =$ $\displaystyle U \overline{\tau - \frac{1}{6} \delta_{ii} \tau}^i$ (2.181)
$\displaystyle F_y$ $\displaystyle =$ $\displaystyle V \overline{\tau - \frac{1}{6} \delta_{ii} \tau}^j$ (2.182)
$\displaystyle F_r$ $\displaystyle =$ $\displaystyle W \overline{\tau - \frac{1}{6} \delta_{ii} \tau}^k$ (2.183)

As for the third order scheme, the best discretization near boundaries is under investigation but currently $ \delta_i \tau=0$ on a boundary.

\fbox{ \begin{minipage}{4.75in} {\em S/R GAD\_C4\_ADV\_X} ({\em gad\_c4\_adv\_x.... ...bf rTrans} (argument) \par $\tau$: {\bf tracer} (argument) \par \end{minipage} }

next up previous contents
Next: 2.16.4 First order upwind Up: 2.16 Linear advection schemes Previous: 2.16.2 Third order upwind   Contents
mitgcm-support@dev.mitgcm.org
Copyright © 2002 Massachusetts Institute of Technology