放射MTGメモ(2012/05/07)

参加者

  • 倉本圭、石渡正樹、はしもとじょーじ、 高橋康人、三上峻、大西将徳

研究の方針・進捗

  • 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/)に置く

次回の日程

  • 5/14(月) 9:00-