dcpam (火星モデル) 作業ミーティング記録 (2013/10/29)

参加者

北大
- 石渡, 荻原
神戸大
- 高橋(芳), 関

荻原

20131008 版で火星計算を試そうとしている.

関

現状報告
- 水過程を入れた計算で図を書いた.
  
  <URL:http://dennou-k.gfd-dennou.org/library/dcpam/sample/2013-10-24_tbseki/> 以下に図を置きつつある.
  - セミラグで単調性を保証するオプションをつけない場合の結果
  - セミラグで単調性を保証するオプションをつけた場合の結果
  両者の結果を比較すると, 今後はオプションをつけた場合でやっていくのが良さそう.
- 上記の結果と水蒸気分布について過去の計算と比較
```
Montmessin et al (2004, JGR 109, E10004) の Ls=90 の図と
計算で得た Ls=90-120 平均の図を比較した.

パターンは大雑把にはあっているみたいに見えるけど, 
水蒸気量は 2 桁くらい小さい.
蒸発量が足らないことが原因の 1 つであろう.
蒸発量が足りない理由の 1 つは, 地表面アルベドの取り扱いが違うこと.
dcpam5 では地球の氷のアルベドの値が使われている.
```
水循環過程をまとめてくる宿題の進捗状況

いくつか論文をみて粒子半径の典型的な数値をチェック中. particle size と書かれている場合に半径か直径かはっきりしないところもあるが, 粒子半径の典型的な値は 2μm ということで良さそう.

重力沈降の式について勉強中.
次の予定
- 作図に関する宿題
  - sufrace snow amount の図は全ての月でも同じに見えるのでチェックする.
  - 今地球の海岸線が図に引かれてしまっているが, これを火星の地形にする.
  - 水氷量の鉛直積分をプロットする.
- 飽和蒸気圧曲線のチェック
  
  今の火星計算で使われている式 (AGCM5 で使っていた式) で出てくる値と Tetens の式で出てくる値とを比較してみる.
- 次の計算
  
  地面の積雪のアルベドを変えた計算をやってみる. constants_snowsea.f90 で定義されている constants_snowseaice_nml の SnowAlbedo の値を変えた計算をやってみる. Montmessin et al (2004) では 0.4 が使われていた.

高橋

dcpam5-0131008 版のバグ報告
- 重力沈降しているルーチンで 0 で割るバグがあった. 計算には使っていないので走るのなら, 結果に影響はしない.
  
  気圧で割る計算のところで最上層の気圧の値(ゼロ)でも割ってしまっている.

次回日程

2013 年 11 月 06 日(水) 14:00 より. TV 会議 (dcmodel meeting2 の部屋) にて