2010/03/31 の dcmodel ネットミーティングのメモ書き

参加者

  • 北大
    • 石渡 正樹, 杉山 耕一朗, 山下 達也
  • 神戸大
    • 高橋 芳幸, 西澤 誠也, 納多 哲史, 今関 翔, 北野 太朗, 島津 通, 井谷 優花, 黒田 美紀, 河合 佑太
  • 九州大
    • 中島 健介

次回日程

  • 日時
    • 04 月 07 日(水) 14:00 - 18:00
  • 場所
    • 神戸大: 自然科学 3 号館 508
    • 北大 : 理学部 8 号館コスモスタジオ (?)
    • 九大 : 理学部 3 号館 3605

deepconv, dcpam に関して新たに行うゼミ

  • deepconv について
    • コード読み
  • dcpam について
    • コード読み
    • 論文読み
    • GCM について
  • 時間の使い方
    • dcpam は別の曜日に行う
    • deepconv は dcmodel ミーティング後に行う

徳永さんの修論

  • 石渡さんに話してもらう
    • 04/14 (水) 16:00 -

スパコン

  • 環境研
    • 林, 中島, 高橋, 納多は山下に誓約書を送る
    • 説明会には山下が出席
    • ジョブは一つしか入らない
  • スパコンの使用者の確認
    • ES: 高橋
      • 結果がまだ来ない
    • NIES: 石渡, 高橋, 納多, 山下, 島津(?)
      • deepconv だけでなく dcpam も走らせなくてはならない
    • 宇宙研: 小高, 杉山, 山下, はしもと(? DCPAM走らせてもらう), 北野(?), 今関(?)
    • 天文台: 杉山, 佐々木
  • 2010 年度に環境研と宇宙研の CPU を主に使う人
    • 納多, 山下
    • もしかしたら今関 (宇宙研), 北野 (環境研 or 宇宙研), 島津 (宇宙研), 馬場 (環境研 or 宇宙研)
    • 山下は夏までは環境研を主に使用
      • 宇宙研に替えるときに馬場・北野のどちらを環境研と宇宙研とに振り分けるかはまた相談する
    • 納多はメインは宇宙研, 足りなくなったら環境研

定式化チェックの作戦

  • 必須メンバー
    • 中島, 石渡, 杉山, 山下
  • 4/7 のミーティング時にチェックをする

進捗状況の報告

納多

  • 熱収支の計算
    • 放射について変更した
      • 補正前の値を使うようにした
        • 3 桁くらいは合っていた
  • 熱輸送を解析するプログラムの完成
  • * 鉛直層数を 32 層にした計算 (地球の年平均日平均日射分布)
    • 太陽定数は地球と同じ
    • フィルタの入れ方がよくないせいか, 平衡状態は 得られていない

杉山

  • deepconv
    • 微分演算の書き方
      • 部分配列の引き算をしている部分を, 添え字でループを回すように変更した
        • 3 % くらい速くなった

高橋

  • mpi の並列化をやり直してみた
    • 3 並列以上でもエラーにならない
    • 入力ファイルは各プロセス用に分割しなければならない
    • 出力ファイルは各プロセスごとに分割される
    • 今後もテストが必要

山下

  • 定式化ドキュメント
    • 付録の辺りを電子化
      • 混合比で書かれていた部分を比湿に書き換えた
      • 乱流エネルギーの式を直した
  • 離散化ドキュメント
    • 修正に若干時間が必要

今関

  • 地表面の熱フラックスについて書き換え中
    • 結果を Web に置く
  • 水蒸気の混合比の初期場の値を検算したところ, 値が一桁違った
    • eccm.f90 での初期場の作成がうまくいってないかもしれない

deepconv のコードについて

  • 乱流拡散項
    • x 方向の微分については \rho(z) はいらないが z 方向の微分では必要かもしれない
    • CReSS ではどうなっているか確認が必要
      • ざっと見た感じでは\rho(z) は残っていそう
    • Klemp and Wilhelmson
      • \rho は入ってない
      • deepconv でも同じ式を使用
    • 書き方
      • 地面におけるフラックスの寄与と, それよりも 上の層のフラックスの項・収束・発散が別れて いる
      • deepconv の書き換え時にまた考える
  • HeatFlux モジュールの xza_MixRtFluxDiff の計算式
    • / ( ( DelZ * 5.0d-1 ) ** 2.0d0 ) について
      • 下部境界の直上の点 (1/2 * DelZ) で計算して いる
      • 1 / DZ ( F - kappa * (MixRt - 0 ) / (DZ/2) ) と考えると, 1/2 を一つ多くかけている?
  • HeatFlux_N1994 モジュールの xz_MixRtFluxBulk について
    • バルク法について勉強して書き換える (今関)
      • 大気科学講座
      • 気象力学通論など