[ 地球流体電脳倶楽部 / dcrtm / SIGEN ]
- 倉本圭、石渡正樹、はしもとじょーじ、 高橋康人、三上峻、大西将徳
- Nakajima et al. 1992の再現
- 進捗
- Nakajima et al. 1992のFig.3を再現
- 非圧縮性気体の圧力(Pn)が負になる問題について、グリッドを細かくすることで解消することを確認
- net radiation flux concergenceが対流圏下層で正になることについては、現在のプログラムではconvergenceの正負が反転する高度が数点存在することを確認
- Fortran90での実装
- Tsを与えて、Ttp、τtpを計算する部分まで実装
- dennouの作法にはよっていない(fortranでプログラムを書いたのち、dennouスタイルで書き直す)
- 議論、再考などが必要な点
- 上空の非圧縮性気体の圧力(Pn)が負になってしまうことについて
- 精度を考慮した積分の仕方を考えるべき
- Nakajima et al. 1992の式(4) (断熱減率の式)の値は、grid上で与えられるので、層の中は一様ではない
- 適当な関数形を与える、ルンゲクッタ法を使うなどして積分を行う
- net radiation flux convergenceが対流圏下層で正になることについて
- 問題のありかを明らかにするために、水蒸気のない大気など、解のわかっているもので検討しなおす
- ToDo
- 精度を考えた積分について
- 断熱減率の積分をルンゲクッタ法で試してみる
- 積分法の教科書の例:数値計算の基礎, 川上一郎 http://www7.ocn.ne.jp/~kawa1/
- 放射の計算がうまくいっているかの確認について
- 水蒸気のない大気(放射平衡など)について手計算とプログラムの結果を比較してみる
- 現在のプログラムを公開場(/GFD_Dennou_Club/dc-arch/dcrtm/model/onishi/)に置く
dcmodel Development Group / GFD Dennou Staff
Last Updated: 2012/05/07, Since: 2012/05/07