dcpam (火星モデル) 作業ミーティング記録 (2014/08/19)

参加者

20130302 版と 20140630-2 版のソースで同期回転惑星計算の結果がどの程度違うかを見てみた.
- 図は以下の通り.
  - 20140630-2 版を用いた結果
  - 20130302 版を用いた結果
- 結果
- OLR, PRCP などの値は 10 W/m**2 程度変わっている.
- 循環の構造は騒ぐほど変わってない.
- 高橋コメント
- 雲の寿命は, 放射収支にあう雲の寿命は設定によって異なる. 全球フラックスが違っているのはこのための可能性がある.
- 鉛直総数も違う. 地球計算で放射収支があうように雲の寿命を決めたときの鉛直層数は L26
次は, 雲の寿命と鉛直総数を変えた計算をやってみる.
現在の dcpam5 における雲の問題
- 地球条件で西風が早すぎる問題に関連して放射加熱率の分布を見てみた.
- APE atlas の長波放射加熱率 (緯度-圧力分布) と比較した.
  - APE のモデル達では, 赤道上の 100 hPa で加熱, それ以外の領域では冷却.
  - dcpam5 では, 熱帯の圏界面の下も加熱になってしまっている. 雲が厚すぎるのか, 放射の計算がおかしいのかどちらか. モデルの雲量は 1 になっているのでそれが悪いのかもしれない.

地表面特性を水平一様にした場合の計算をおこなった.
結果の図を並べたページは
- 火星ダストと巻き上げ実験 DustScenario Const 0.2 DD巻き上げあり
熱慣性, アルベド, 地形を全て一定. 全て全球平均値を与える. Ls 30 ごとの平均 (季節平均) 季節平均場では経度方向の変化がほとんど無くなる.
スナップショットと動画を見てみた.

やはりダストデビルによる巻き上げフラックスは結構複雑に決まっているようだ. 地表面顕熱フラックスも熱効率も複雑な分布をしている. 単純に巻き上げフラックスのピークが地表面顕熱フラックスのピークもしくは熱効率のピークに対応しているわけではなさそう. 線グラフを作って比較してみる必要がある (時刻と緯度を決めて, DDAF, 顕熱, 熱効率の 3本のグラフを書く).
地表面特性を一定にしない場合に比べて巻き上げ量がどの程度変わるかは, これからちゃんと比較できる図を作る予定.