dcmodel ミーティング記録(2012/12/19)

参加者

経過報告
- 水平/鉛直セミラグランジアン計算
- 上流点探索
- 補間 : エルミート 3 次(by Enomoto 2008)
  - 水平スペクトル微分値
  - 鉛直 4 次中央差分
- 非負値保証のフィルター
- dcpam に仮導入
短波問題
- 非負値はフィルターで保証されている
- ギブス現象により短波が発生
- Sun et al. (1996) フィルタ
  - 4 点エルミート補間後に極値を持つ場合は Global な最大・最小で押さえる
  - 単調な場合は local な最大・最小で押さえる
  - 長時間計算で振幅が減少
  - Global 最大・最小で押さえるのをやめると振幅減少が抑えられる
さらなる高精度化
- フィルターを入れることで補間に要する点が増える
- そこで初めから 4 点の値と微分値を用いた高精度補間とフィルターの組み合わせを試した
- 5 次から 7 次はあまり変わらない, 7 次は計算コストに見合わない
- 球面での定式化
鉛直移流
- 微分値が中央差分で補間精度が上がらない
- 境界の扱い, 境界直上・直下の鉛直流が悪さをしている?
- 境界でのセル内の分布をどう仮定するか?
低解像度水惑星実験(T21L16)
- セミラグのほうがシャープな水蒸気分布
- 単にセミラグの分布をなめしたのがオイラーの分布になっているわけではない

dcpam5 へのダスト巻き上げスキーム(Kahre et al. 2006)の組込み
トレーサ移流実験
- T21L32, 3 火星年
- 放射計算するためのダスト分布は外から与える. 水平一様の分布
- 結果は dennou dcpam のページに載せてある