[MITgcm-devel] multidimadvection and manual

[Martin Losch]

Hallo,

I know I should not "complain" about the manual without contributing  
to it, but I think that the description of the multi dimension  
advection is slightly wrong:
http://mitgcm.org/r2_web_testing/latest/online_documents/node81.html
eq.2.201,2.202,2.203
should be replaced by what is written in the head of gad_advection.F  
(which appears to describe what really happens in gad_advection.F)
> C The algorithm is as follows:
> C \begin{itemize}
> C \item{$\theta^{(n+1/3)} = \theta^{(n)}
> C      - \Delta t \partial_x (u\theta^{(n)}) + \theta^{(n)}  
> \partial_x u$}
> C \item{$\theta^{(n+2/3)} = \theta^{(n+1/3)}
> C      - \Delta t \partial_y (v\theta^{(n+1/3)}) + \theta^{(n)}  
> \partial_y v$}
> C \item{$\theta^{(n+3/3)} = \theta^{(n+2/3)}
> C      - \Delta t \partial_r (w\theta^{(n+2/3)}) + \theta^{(n)}  
> \partial_r w$}
> C \item{$G_\theta = ( \theta^{(n+3/3)} - \theta^{(n)} )/\Delta t$}
> C \end{itemize}
(that is theta(n+2/3) = theta(n+1/3) - ... and theta(n+3/3)=theta(n 
+2/3) - ... etc. instead of theta(n+2/3) = theta(n) - ...)

Martin