Neale and Hoskins (2000) による水惑星実験 (軸対称計算)

Neale and Hoskins (2000) によって提案された, 惑星表面が全て水であるような条件 (水惑星) で, 軸対称実験を行う方法を説明します.

この計算には以下の物理過程を用いています.

概要

本実験では, 軸対称モデルを用いて計算を行います. 軸対称版力学過程を使うためには, FFLAGS に -DAXISYMMETRY または -DAXISYMMETRY_SJPACK を 与えてビルドしなおす必要があることに注意してください.

本実験では, 海表面温度のデータが必要です. したがって, 実験の実行は以下の 4 つのステップで行います.

実験用ディレクトリ準備

dcpam5 の数値実験はソースツリー内部では行わず ソースツリーとは別の外部ディレクトリにて行うことを推奨いたします.

まず dcpam5 ソースのトップディレクトリ(以下の例では dcpam5-YYYYMMDD とする)に移動してください. 以下では dcpam5 ソースディレクトリの隣に ../dcpma5-exp/ape-nh00-axisym-exp ディレクトリを作成し, そこで実験を行うことにします. 次のように ../dcpam5-exp/ape-nh00-axisymexp ディレクトリを作成し, 移動してください.

$ mkdir -p ../dcpam5-exp/ape-nh00-axisym-exp
$ cd ../dcpam5-exp/ape-nh00-axisym-exp

次に, このディレクトリに実行ファイルと設定ファイル置き場を作成します.

$ mkdir bin
$ mkdir conf

最後に作成したディレクトリに実行ファイルと設定ファイルをコピーします.

$ cp ../../dcpam5-YYYYMMDD/src/main/dcpam_main bin
$ cp ../../dcpam5-YYYYMMDD/src/main/dcpam_init_data  bin
$ cp ../../dcpam5-YYYYMMDD/src/main/dcpam_init_data_surface   bin
$ cp ../../dcpam5-YYYYMMDD/exp_setup_files/*.conf conf

なお, 実行ファイルと設定ファイル (NAMELIST ファイル) があれば, どのディレクトリにおいても計算を行うことができます.

初期値データファイルの作成

dcpam_init_data と init_data_APE_NH00_Axisym_T42L22.conf を用いて初期値ファイル init_T42L22.nc を作成します.

$ ./bin/dcpam_init_data -N=./conf/init_data_APE_NH00_Axisym_T42L22.conf

 *** MESSAGE [init_data] ***  Run: Initial data generation
 *** MESSAGE [init_data] ***  -- version = $Name: dcpam5-20140228 $$Id: exp-APE-NH00-Axisym.rd,v 1.4 2013-10-06 13:46:11 yot Exp $
 *** MESSAGE [namelist_util] ***  ----- Initialization Messages -----
 *** MESSAGE [namelist_util] ***    MaxNmlArySize = 256
 *** MESSAGE [namelist_util] ***  -- version = $Name: dcpam5-20140228 $$Id: exp-APE-NH00-Axisym.rd,v 1.4 2013-10-06 13:46:11 yot Exp $
 *** MESSAGE [timeset] ***  NAMELIST group "timeset_nml" is loaded from "./conf/init_data_APE_NH00_Axisym_T42L22.conf".
                       : 
 *** MESSAGE [initial_data] ***    PsAvr   = 101080.
 *** MESSAGE [initial_data] ***    QVapAvr = 0.
 *** MESSAGE [initial_data] ***    Ueq     = 0.
 *** MESSAGE [initial_data] ***
 *** MESSAGE [initial_data] ***  -- version = $Name: dcpam5-20140228 $$Id: exp-APE-NH00-Axisym.rd,v 1.4 2013-10-06 13:46:11 yot Exp $
 *** MESSAGE [HistoryClose] ***  "init_T42L22.nc" is closed

 ############## CPU TIME SUMMARY ################
 others                 0.400020E-01
 ------------------------------------------------
        TOTAL TIME =    0.400020E-01

実験用データの準備

実験用データとして, オゾン分布のデータを 以下に用意してあります.

上記のページから "O3_NH00_lon1_T042.nc" をダウンロードして, 実験を行うディレクトリに置いてください.

SST データファイルの作成

次に, dcpam_init_data_surface と sst_data_APE_NH00_Axisym_T42.conf を用いて, 海水面データ sst_T42.nc を作成します.

$ ./bin/dcpam_init_data_surface -N=./conf/sst_data_APE_NH00_Axisym_T42.conf

 *** MESSAGE [sst_data] ***  Run: SST data generation
 *** MESSAGE [sst_data] ***  -- version = $Name: dcpam5-20140228 $$Id: exp-APE-NH00-Axisym.rd,v 1.4 2013-10-06 13:46:11 yot Exp $
 *** MESSAGE [namelist_util] ***  ----- Initialization Messages -----
 *** MESSAGE [namelist_util] ***    MaxNmlArySize = 256
 *** MESSAGE [namelist_util] ***  -- version = $Name: dcpam5-20140228 $$Id: exp-APE-NH00-Axisym.rd,v 1.4 2013-10-06 13:46:11 yot Exp $
 *** MESSAGE [timeset] ***  NAMELIST group "timeset_nml" is loaded from "./conf/sst_data_APE_NH00_Axisym_T42.conf".
                       :
 *** MESSAGE [surface_data] ***    SoilHeatCap      = 2100000.
 *** MESSAGE [surface_data] ***    SoilHeatDiffCoef = 1.2
 *** MESSAGE [surface_data] ***  -- version = $Name: dcpam5-20140228 $$Id: exp-APE-NH00-Axisym.rd,v 1.4 2013-10-06 13:46:11 yot Exp $
 *** MESSAGE [HistoryClose] ***  "sst_T42.nc" is closed

 ############## CPU TIME SUMMARY ################
 others                 0.400000E-02
 ------------------------------------------------
        TOTAL TIME =    0.400000E-02

実験の実行

最後に dcpam_main, dcpam_APE_NH00_Axisym_T42L22.conf を用いて実験を実施します. 解像度 T42L21, 時間ステップ 12 分で 2 日分計算を行います. リスタートデータ, 地表面リスタートデータといくつかのヒストリデータ ファイルが出力されます.

$ ./bin/dcpam_main -N=./conf/dcpam_APE_NH00_Axisym_T42L22.conf | tee ape-nh00.log

初期値データ, 地表面データ, 実験設定を変更する場合には, 上記の実行に際して conf ディレクトリ内の init_data_APE_NH00_Axisym_T42L22.conf, sst_data_APE_NH00_Axisym_T42.conf, dcpam_APE_NH00_Axisym_T42L22.conf を変更してください.

計算結果の解析・可視化の方法についてはこちらを参照してください.

最後に

実験ディレクトリに作成された計算結果を保存しておきたい場合には データ用ディレクトリ(例えば data01)を作成し, ファイルを移動させて整理しておきましょう.

$ mkdir data01
$ mv *.nc *.log data01/

なお, 計算結果を計算実行時にサブディレクトリへ出力することも可能です. 実験設定の変更の出力設定の変更を参照し, NAMELIST "&gtool_historyauto_nml" の変数 "FilePrefix" による設定を試してみて下さい.

参考文献


$Id: exp-APE-NH00-Axisym.rd,v 1.4 2013-10-06 13:46:11 yot Exp $