放射MTGメモ(2012/10/16)

参加者

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

Nakajima et al. 1992の再現(大西)

• dennouスタイルでのプログラムのリライト
  • 以下の3点をdennouスタイルに合わせる作業中
    • 1. 変数の書き換え
    • 2. コメントの書き方をそろえる(RDOC処理ができる形式)
    • 3. データの入出力(gtoolを使ったnetCDF出力)
      • gtool5のインストールは完了

同期回転惑星の放射計算

• 納多さんの同期回転惑星のGCMの計算結果を議論する材料として、同期回転惑星の条件で射出限界の計算を行う
  • 計算条件(石渡さんより)
    • 基本的に Nakajima et al. が用いた標準的なパラメータと同じ値
      • g=9.8, kappa_v=0.01 m^2 Kg^{-1}, kappa_n=0.0, p_{n0}=10^5 Pa,
    • ただし, 分子量, 比熱を変える.
      • m_n = 28.964 \times 10^{-3} Kg mol^{-1}
      • m_v = 18.01528 \times 10^{-3} Kg mol^{-1}
      • c_{pn} = 1004.6 J Kg^{-1} K^{-1}
      • c_{pv} = 1810.0 J Kg^{-1} K^{-1}
  • 計算結果
    • 地表面温度を250, 290, 300, 310, 330, 350, 400, 450, 500[K]で計算(450, 500[K]は計算中)
    • 上向きフラックスの最大値は381[W/m2](地表面温度が300[K]で最大)
    • 結果のグラフ(地表面温度に対する上向きフラックス)は /GFD_Dennou_Club/dc-arch/dcrtm/model/onishi/syncrot/ に格納(計算プログラムも)
  • 追加計算など
    • 今回の結果を受け、どのような計算が追加で必要かを石渡さんに確認
      • 地表面温度300[K]付近の温度を細かく計算する、とか、湿度を80[%]で計算する、など

木星のガスの吸収線について(高橋)

• 線吸収の計算
  • 以前の計算(バンドの吸収を与えた計算)に比べ光学的厚さが大きく計算される問題があったが、定式化に誤りがあったことが判明。修正を行う。
• Fluxの計算
  • 計算時間を短縮するために計算の仕方を変更したが(全層を逐一計算するコードから、ひとつ前の層のフラックスから計算するプログラムへ変更)、散乱の角度分布を計算する現在のアルゴリズム(AH)ではうまくいかない。そのため、係数(5/3)で近似したアルゴリズムで計算を試みる
• 惑星科学会の発表スライド
  • 結果を説明するようにイントロをつくる

次回の日程

• 10/23(火) 9:00-