48 |
C The algorithm is as follows: |
C The algorithm is as follows: |
49 |
C \begin{itemize} |
C \begin{itemize} |
50 |
C \item{$\theta^{(n+1/3)} = \theta^{(n)} |
C \item{$\theta^{(n+1/3)} = \theta^{(n)} |
51 |
C - \Delta t \partial_x (u\theta) + \theta \partial_x u$} |
C - \Delta t \partial_x (u\theta^{(n)}) + \theta^{(n)} \partial_x u$} |
52 |
C \item{$\theta^{(n+2/3)} = \theta^{(n+1/3)} |
C \item{$\theta^{(n+2/3)} = \theta^{(n+1/3)} |
53 |
C - \Delta t \partial_y (v\theta) + \theta \partial_y v$} |
C - \Delta t \partial_y (v\theta^{(n+1/3)}) + \theta^{(n)} \partial_y v$} |
54 |
C \item{$\theta^{(n+3/3)} = \theta^{(n+2/3)} |
C \item{$\theta^{(n+3/3)} = \theta^{(n+2/3)} |
55 |
C - \Delta t \partial_r (w\theta) + \theta \partial_r w$} |
C - \Delta t \partial_r (w\theta^{(n+2/3)}) + \theta^{(n)} \partial_r w$} |
56 |
C \item{$G_\theta = ( \theta^{(n+3/3)} - \theta^{(n)} )/\Delta t$} |
C \item{$G_\theta = ( \theta^{(n+3/3)} - \theta^{(n)} )/\Delta t$} |
57 |
C \end{itemize} |
C \end{itemize} |
58 |
C |
C |