使っていて、1ヶ月以上にわたり悩まされた問題について、メモ代わりに残しておきたい。書いておくと、自分は初期値として水平に密度勾配があるという条件について数値実験していた。このことが問題となっていた。この海洋モデルでは座標系のうち鉛直方向を、一般的な直交座標系ではなくシグマ座標系というのを採用している。シグマ座標についてはほかにいくらでも書かれているので、ここでは詳細に記載しない。メリットもあり、デメリットもある。海洋において圧力とは、主に水の密度である。この密度が水平に異なるとき、高密度の方から低密度の方へ流速が起きるような力が生まれる。これを圧力傾度力と呼ぶ。シグマ座標系において圧力傾度力の計算で起きる密度の打ち切り誤差は、圧力傾度力に対して過大に影響する。シグマ座標系の問題として長く議論されている。することで、密度の絶対値ではなく相対値によって圧力傾度力が計算され、打ち切り誤差が軽減される。