地球的惑星の実験を行う方法を説明します.
この計算には以下の物理過程を用いています.
本実験では, 海表面温度, 地面特性, 海表面密度, オゾン分布のデータが必要です. したがって, 実験の実行は以下の 4 つのステップで行います.
practice ディレクトリに移動してください.
$ cd practice
そのディレクトリ内で, make コマンドを実行してください. 自動的に実行ファイルや NAMELIST ファイルがコピーされます. make コマンド入力後, まずディレクトリ名を適宜入力してください. (下記の例では Earth-exp としています). それ以降の質問はとりあえずそのまま Enter キーを入力してください.
$ make ****** Setup a directory for a experiment ****** Enter directory name [testXX]: Earth-exp Directory in which executable files are prepared [../src/main]: :
そしてそのディレクトリに移動します.
$ cd Earth-exp
なお, このディレクトリ内以外でも, 実行ファイル, 設定ファイル (namelist ファイル) があれば計算を行うことはできます.
init_data と init_data_E_T21L22.nml を用いて初期値ファイル init_T21L22.nc を作成します.
$ ./init_data -N=init_data_E_T21L22.nml *** MESSAGE [init_data] *** Run: Initial data generation *** MESSAGE [init_data] *** -- version = $Name: dcpam5-20110327 $$Id: exp-Earth.rd,v 1.4 2011-02-21 07:19:54 yot Exp $ *** MESSAGE [namelist_util] *** ----- Initialization Messages ----- *** MESSAGE [namelist_util] *** MaxNmlArySize = 256 *** MESSAGE [namelist_util] *** -- version = $Name: dcpam5-20110327 $$Id: exp-Earth.rd,v 1.4 2011-02-21 07:19:54 yot Exp $ *** MESSAGE [timeset] *** NAMELIST group "timeset_nml" is loaded from "init_data_E_T21L22.nml". : *** MESSAGE [initial_data] *** Ueq = 0. *** MESSAGE [initial_data] *** *** MESSAGE [initial_data] *** -- version = $Name: dcpam5-20110327 $$Id: exp-Earth.rd,v 1.4 2011-02-21 07:19:54 yot Exp $ *** MESSAGE [HistoryClose] *** "init_T21L22.nc" is closed ############## CPU TIME SUMMARY ################ others 0.120000E-01 ------------------------------------------------ TOTAL TIME = 0.120000E-01
実験用データとして, 海表面温度, 地面特性, 海表面密度, オゾン分布のデータを 以下に用意してあります.
上記のページから "sst_amipII_bc_clim_T021.nc", "sp_for_Earth_T021.nc", "sic_amipII_bc_clim_T021.nc", "O3_CMIP5_climatology_zonalmean_T021.nc" をダウンロードして, 実験を行うディレクトリに置いてください.
最後に dcpam_main, dcpam_E_T21L22.nml を用いて実験を実施します. 解像度 T21L22, 時間ステップ 30 分で 2 日分計算を行います. リスタートデータ, 地表面リスタートデータといくつかのヒストリデータ ファイルが出力されます.
$ ./dcpam_main -N=dcpam_E_T21L22.nml | tee Earth.log
初期値データ, 地表面データ, 実験設定を変更する場合には, 上記の実行に際して, init_data_E_T21L22.nml, dcpam_E_T21L22.nml を変更してください.