[MITgcm-devel] multidimadvection and manual
Jean-Michel Campin
jmc at ocean.mit.edu
Thu Apr 20 09:11:37 EDT 2006
Hi Martin,
Thanks for the report (I read this part.2 last time
and did not see it)
I added this to the "manual_todo.txt" file
Cheers,
Jean-Michel
On Thu, Apr 20, 2006 at 11:36:50AM +0200, Martin Losch wrote:
> 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
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel
More information about the MITgcm-devel
mailing list