[地球流体電脳倶楽部]
[dcmodel |
dcmodel-tools]
[numexp]
[spmodel]
SPMODEL 球面浅水モデルを用いた
Showman and Polvani (2010) の計算 (並列計算版)
このページには SPMODEL 球面浅水モデル (shallow-zd) を用いた
Showman and Polvani (2010) の結果を示している.
藤林さんのプログラムを参考に並列計算用にプログラムを書きかえ,強制項が加わっても並列計算ができることを確かめる.
- 参考文献
- 計算に使用した計算機
- nozomi
- CPU : Intel(R) Core(TM) i7-9700K CPU @ 3.60GHz
方程式
\begin{align}
\frac{D \boldsymbol{u}}{D t} + g' \nabla_z h + f \boldsymbol{k} \times \boldsymbol{u} &= \boldsymbol{R} - \frac{\boldsymbol{u}}{\tau_{drag}} - \boldsymbol{D}_u, \\
\frac{D h}{D t} + h (\nabla_z \cdot \boldsymbol{u}) &= S - \frac{h-H}{\tau_{rad}} - D_h . \\
\end{align}
もしくは
\begin{align}
\frac{\partial \boldsymbol{u}}{\partial t} + (\boldsymbol{u} \cdot \nabla_z)\boldsymbol{u} &= - f \boldsymbol{k} \times \boldsymbol{u} - g' \nabla_z h + \boldsymbol{R} - \frac{\boldsymbol{u}}{\tau_{drag}} - \boldsymbol{D}_u, \\
\frac{\partial h}{\partial t} + \nabla_z \cdot (\boldsymbol{u} h) &= S - \frac{h - H }{\tau_{rad}} - D_h = Q - D_h .
\end{align}
運動量の式には下層から上層への運動量の移流を表す\(R\), レイリー摩擦,
連続の式には質量のソースとシンクを表す\(S\), ニュートン冷却の項が含まれる.
標準実験
ジオポテンシャル
spml2 T170 8並列
流線関数
spml2 T170 8並列