計算例: MPDATA

MPDATA ((22)$\sim$ (25)式)を用いて同様の 計算を行なった. なお $\Delta t=0.1$ に戻してある. 途中必要となる 平均値 計算は(26)$\sim$ (28)に従い求めた. これらの平均操作については Smolarkiewicz(1982) の記述を参考にしてある.

1 回転後(628ステップ)と 3 回転後(1884ステップ)後の結果 をFig.6とFig.7にそれぞれ示す. 初期分布は比較的 よく保たれている. 2 次中心差分で発生したような波も生じていない.

図 6: MPDATA による計算. 1 回転(628ステップ)後の結果.
図 7: MPDATA による計算. 3 回転(1884ステップ)後の結果.