01 の GCM セミナーのメモ書き

ToDo リスト

一番上にはギャラリーを持ってくる. ソースコードはその後ろ.

rdoc-f95 の書式に合わせてソースコードのコメントを記述する.

ソースコードに物理過程が加えられたのに伴い, ドキュメントにも物理過程部分を追加する.

例: rikigaku -> dynamics

agcm5 の第3章「コード解説」に相当する部分の挿入. 変数に関しては dcpam に置き換える.

モジュール名の末尾の "_mod" はやめる
- 文字数が不足しつつあるため
  - "_mod" に 4 文字使うよりも, よりモジュールの内容を示すことにその文字数を使うほうが良い. (F90/95 規格では言語要素の文字数は 31 文字まで)
- モジュール名が登場するのは定義の部分と参照する部分であり, "module" や "use" という表記からそれがモジュールであることは自明であるため, あえて名前そのものにモジュールであることを示すサフィックスをつける必要は無い.
力学過程, 物理過程であることを示すプレフィックスは短く
- 力学過程: dyn_
- 物理過程: phys_
- 理由
  - よく知られた「力学過程」, 「物理過程」を示す文字列を長くするより, より詳細な部分を示すことに文字列を使いたい.

現在の explicit を semi-implicit に改造する.

以下の作業はまだこれから.

init の中身を dcpam.f90 内に入れてしまう. ただし, 他の初期値を利用したい場合にはその初期値を入力できるよう NAMELIST を読み込むようにしておく.
データ I/O の厄介な部分は gt4 で吸収するようにする.

リスタートファイルを組み込む機構を取り入れる際には, spmodel, arare などと合わせる.

T63L20 で初期値 300K, 250K の 2 つでの計算は終了. 「計算結果」置き場にすべて置いた.

時間があれば初期値をより変えたものを走らせておいてもいいかもしれないが, まずはリスタートファイルを入出力できるようにするのが急務

dcpam に対して

'-w -sx6 -Wf"-pvctl nodivloop noloopchg -O"'

というオプションをつけているが, 単に計算速度を低下させているだけかもしれない. (少なくとも AGCM5 にはこれを指定せずともうまく動いている. ISPACK には必要なオプションだが)

DCPAM のコンパイルのみこれをはずしてみて, 計算速度を計測してみたが, 数％ぐらいしか違わなかった.

今度は, ispack, gt4f90io, spmodel それぞれでも上記のオプションをはずして計算してみて速度を計測する.

なお, ispack はループ入れ替えなどの最適化をすると実行速度が極端に落ちるらしい. (ispack の README 参照).

ISPACK と spmodel の OpenMP による並列化を利用し, DCPAM の並列計算を行なった. 計算効率は以下の通り.

CPU 数     CPU 時間

1          1.0

2          1.4

4          1.5

残念ながら (ISPACK や spmodel を利用していない部分は全然チューニングしてないので当然かもしれないが), 今のところ並列化効率は全然良くない. 全体の並列化効率を向上させるには, ISPACK 以外の部分の並列化も考慮すべきかも.