=begin JA = 火星サーマル上昇実験 # * 山下 達也 (yamasita), 杉山 耕一朗 (sugiyama), 小高 正嗣 (odakker) # * $Id: mars-plumetest.rd,v 1.1 2011-02-24 09:25:48 yamasita Exp $ =end JA =begin EN = A Numerical experiment of ascending hot plume # * Tatsuya YAMASHITA (yamasita), Ko-ichiro SUGIYAMA (sugiyama), Masatsugu ODAKA (odakker) # * $Id: mars-plumetest.rd,v 1.1 2011-02-24 09:25:48 yamasita Exp $ =end EN =begin JA 単独のホットプリュームを上昇させる実験を行なう方法を説明します. #この計算には以下の物理過程を用いています. # # * 放射 (AGCM5 デフォルト) # * 長波放射 (水蒸気 3 領域, 乾燥大気 1 領域) # * 短波放射 (水蒸気 1 領域) # * 鉛直拡散 (Mellor and Yamada, 1974, レベル 2) # * 積雲パラメタリゼーション: 湿潤対流調節 (Manabe et al., 1965) # * 診断型大規模凝結 # * 乾燥対流調節 =end JA =begin EN A method to perform an aqua planet experiment (planetary surface is covered by watar) is described. #Following physical processes are used in this experiment. # # * Radiation (Band model) # * Longwave radiation (moisture 3 band, dry air 1 band) # * Shortwave radiation (moisture 1 band) # * Vertical diffusion (Mellor and Yamada, 1974, level 2) # * Cumulus parameterization: Convective adjustment scheme (Manabe et al., 1965) # * Large scale condensation # * Dry convective adjustment =end EN =begin JA == 概要 本実験では, 海表面温度のデータが必要です. したがって, 実験の実行は以下の 4 つのステップで行います. * 実験ディレクトリの準備 * 初期値の準備 * 実験用データ (海表面温度分布) の準備 * 実験の実行 =end JA =begin EN == Overview This experiment requires data of distribution of sea surface temperature. So, the experiment is performed with the following 4 steps: * Preparation of directory for experiments * Preparation of initial condition * Preparation of additional data (sea surface temperature) for experiments * Execution of experiments =end EN =begin JA == 実験用ディレクトリ準備 practice ディレクトリに移動してください. $ cd practice そのディレクトリ内で, make コマンドを実行してください. 自動的に実行ファイルや NAMELIST ファイルがコピーされます. make コマンド入力後, まずディレクトリ名を適宜入力してください. (下記の例では ape-exp としています). それ以降の質問はとりあえずそのまま Enter キーを入力してください. $ make ****** Setup a directory for a experiment ****** Enter directory name [testXX]: ape-exp Directory in which executable files are prepared [../src/main]: : そしてそのディレクトリに移動します. $ cd ape-exp なお, このディレクトリ内以外でも, 実行ファイル, 設定ファイル (namelist ファイル) があれば計算を行うことはできます. == 初期値データファイルの作成 init_data と init_data_T21L16.nml を用いて初期値ファイル init_T21L16.nc を作成します. $ ./init_data -N=init_data_T21L16.nml : *** MESSAGE [restart_file_io] *** ----- Initialization Messages ----- *** MESSAGE [restart_file_io] *** Input:: *** MESSAGE [restart_file_io] *** InputFile = *** MESSAGE [restart_file_io] *** Output:: *** MESSAGE [restart_file_io] *** OutputFile = init_T21L16.nc *** MESSAGE [restart_file_io] *** IntTime = 1. [day] *** MESSAGE [HistoryCreate3] *** "init_T21L16.nc" is created (origin=0. [day]) *** MESSAGE [initial_data] *** NAMELIST group "initial_data_nml" is loaded from "init_data_T21L16.nml". *** MESSAGE [initial_data] *** NAMELIST group "initial_data_nml" is loaded from "init_data_T21L16.nml". *** MESSAGE [initial_data] *** ----- Initialization Messages ----- *** MESSAGE [initial_data] *** Pattern = Small Disturbance of Temperature *** MESSAGE [initial_data] *** TempAvr = 280. *** MESSAGE [initial_data] *** PsAvr = 100000. *** MESSAGE [initial_data] *** QVapAvr = 0.1E-09 *** MESSAGE [HistoryClose] *** "init_T21L16.nc" is closed == SST データファイルの作成 次に, sst_data と sst_data_T21.nml を用いて, 海水面データ sst_T21.nc を作成します. $ ./sst_data -N=sst_data_T21.nml : *** MESSAGE [restart_surftemp_io] *** ----- Initialization Messages ----- *** MESSAGE [restart_surftemp_io] *** Input:: *** MESSAGE [restart_surftemp_io] *** InputFile = *** MESSAGE [restart_surftemp_io] *** InputName = SurfTemp *** MESSAGE [restart_surftemp_io] *** Output:: *** MESSAGE [restart_surftemp_io] *** OutputFile = sst_T21.nc *** MESSAGE [restart_surftemp_io] *** IntTime = 1. [day] *** MESSAGE [HistoryCreate3] *** "sst_T21.nc" is created (origin=0. [day]) *** MESSAGE [surface_data] *** NAMELIST group "surface_data_nml" is loaded from "sst_data_T21.nml". *** MESSAGE [surface_data] *** NAMELIST group "surface_data_nml" is loaded from "sst_data_T21.nml". *** MESSAGE [surface_data] *** ----- Initialization Messages ----- *** MESSAGE [surface_data] *** Pattern = Hosaka et al. (1998) *** MESSAGE [surface_data] *** SurfTemp = 302. *** MESSAGE [surface_data] *** Albedo = 0.15 *** MESSAGE [surface_data] *** HumidCoeff = 1. *** MESSAGE [surface_data] *** RoughLength = 0.1E-03 *** MESSAGE [surface_data] *** HeatCapacity = 0. *** MESSAGE [surface_data] *** TempFlux = 0. *** MESSAGE [surface_data] *** SurfCond = 0 *** MESSAGE [HistoryClose] *** "sst_T21.nc" is closed == 実験の実行 最後に dcpam_main, dcpam_ape_T21L16.nml を用いて実験を実施します. 解像度 T21L16, 時間ステップ 20 分で 10 日分計算を行います. リスタートデータ, 地表面リスタートデータといくつかのヒストリデータ ファイルが出力されます. $ ./dcpam_main -N=dcpam_ape_T21L16.nml | tee ape.log 初期値データ, 地表面データ, 実験設定を変更する場合には, 上記の実行に際して, init_data_T21L16.nml, sst_data_T21.nml, dcpam_ape_T21L16.nml を変更してください. == リスタートデータからの実行 上記の初期値データを元に, 1000 日分の計算を行った際のリスタートデータを 以下に用意してあります. このデータを用いてリスタートを行ってみましょう. * ((<サンプルリスタートデータ|URL:http://www.gfd-dennou.org/library/dcpam/sample/2009-03-19_tutorial/ape-T21L16>)) 上記のページから "ape-T21L16-rst.nc" と "ape-T21L16-rst-rad.nc" をダウンロードしてください. そして, "dcpam_ape_T21L16.nml" の ×et_nml, &restart_file_io, &radiation_band_nml を以下のように変更してください. ×et_nml RestartTimeValue = 1000.0, ! 計算開始時刻. ! Start time of calculation RestartTimeUnit = 'day', ! 計算開始時刻の単位. ! Unit of start time of calculation IntegPeriodValue = 10.0, ! 計算終了時刻. ! End time of calculation IntegPeriodUnit = 'day', ! 計算開始時刻の単位. ! Unit of end time of calculation : &restart_file_io_nml InputFile = 'ape-T21L16-rst.nc', ! 入力するリスタートデータのファイル名 ! filename of input restart data : &radiation_band_nml RstInputFile = 'ape-T21L16-rst-rad.nc', ! 入力するリスタートデータのファイル名 ! Filename of input restart data 以下のコマンドで実行を行ってください. $ ./dcpam_main -N=dcpam_ape_T21L16.nml | tee ape.log この結果得られる Vor.nc の変数 Vor を GAVE で開くと 以下のような渦度の分布が表示されます. ((<"IMG:images/gave_ape_vor.png">)) =end JA =begin EN == Preparation of a directory for an experiment Move "practice" directory: $ cd practice In the directory, execute 'make' command. The execution files and the NAMELIST files will be automatically copied. Input directory name first of all after inputting 'make' command. (It is assumed "ape-exp" in the following example). Input 'Enter' to the question after that. $ make ****** Setup a directory for a experiment ****** Enter directory name [testXX]: ape-exp Directory in which executable files are prepared [../src/main]: : Move the directory $ cd ape-exp Note that you can perform an experiment in any directory by using executable files and configuration (namelist) files. == Create initial data file Create initial data file "init_T21L16.nc" using "init_data" and "init_data_T21L16.nml" $ ./init_data -N=init_data_T21L16.nml : *** MESSAGE [restart_file_io] *** ----- Initialization Messages ----- *** MESSAGE [restart_file_io] *** Input:: *** MESSAGE [restart_file_io] *** InputFile = *** MESSAGE [restart_file_io] *** Output:: *** MESSAGE [restart_file_io] *** OutputFile = init_T21L16.nc *** MESSAGE [restart_file_io] *** IntTime = 1. [day] *** MESSAGE [HistoryCreate3] *** "init_T21L16.nc" is created (origin=0. [day]) *** MESSAGE [initial_data] *** NAMELIST group "initial_data_nml" is loaded from "init_data_T21L16.nml". *** MESSAGE [initial_data] *** NAMELIST group "initial_data_nml" is loaded from "init_data_T21L16.nml". *** MESSAGE [initial_data] *** ----- Initialization Messages ----- *** MESSAGE [initial_data] *** Pattern = Small Disturbance of Temperature *** MESSAGE [initial_data] *** TempAvr = 280. *** MESSAGE [initial_data] *** PsAvr = 100000. *** MESSAGE [initial_data] *** QVapAvr = 0.1E-09 *** MESSAGE [HistoryClose] *** "init_T21L16.nc" is closed == Create SST data file Next, create SST data "sst_T21.nc" using "sst_data" and "sst_data_T21.nml". $ ./sst_data -N=sst_data_T21.nml : *** MESSAGE [restart_surftemp_io] *** ----- Initialization Messages ----- *** MESSAGE [restart_surftemp_io] *** Input:: *** MESSAGE [restart_surftemp_io] *** InputFile = *** MESSAGE [restart_surftemp_io] *** InputName = SurfTemp *** MESSAGE [restart_surftemp_io] *** Output:: *** MESSAGE [restart_surftemp_io] *** OutputFile = sst_T21.nc *** MESSAGE [restart_surftemp_io] *** IntTime = 1. [day] *** MESSAGE [HistoryCreate3] *** "sst_T21.nc" is created (origin=0. [day]) *** MESSAGE [surface_data] *** NAMELIST group "surface_data_nml" is loaded from "sst_data_T21.nml". *** MESSAGE [surface_data] *** NAMELIST group "surface_data_nml" is loaded from "sst_data_T21.nml". *** MESSAGE [surface_data] *** ----- Initialization Messages ----- *** MESSAGE [surface_data] *** Pattern = Hosaka et al. (1998) *** MESSAGE [surface_data] *** SurfTemp = 302. *** MESSAGE [surface_data] *** Albedo = 0.15 *** MESSAGE [surface_data] *** HumidCoeff = 1. *** MESSAGE [surface_data] *** RoughLength = 0.1E-03 *** MESSAGE [surface_data] *** HeatCapacity = 0. *** MESSAGE [surface_data] *** TempFlux = 0. *** MESSAGE [surface_data] *** SurfCond = 0 *** MESSAGE [HistoryClose] *** "sst_T21.nc" is closed == Run the experiment At the last, run a experiment using "dcpam_main" and "dcpam_ape_T21L16.nml". Resolution is T21L16, time step is 20 minutes. Integration time is 7 days. A restart data file, a restart data file of surface data, and some history data files are output. $ ./dcpam_main -N=dcpam_ape_T21L16.nml | tee ape.log If initial data, surface data, settings of experiments want to be changed, edit "init_data_T21L16.nml", "sst_data_T21.nml", "dcpam_ape_T21L16.nml". =end EN =begin JA == 参考文献 =end JA =begin EN == References =end EN =begin * Manabe, S., J. Smagorinsky, and R. F. Strickler, 1965: Simulated climatology of a general circulation model with a hydrologic cycle, ((|Mon. Wea. Rev.|)), ((*93*)), 769--798. * Mellor, G. L., and T. Yamada, 1974: A hierarchy of turbulence closure models for planetary boundary layers, ((|J. Atmos. Sci.|)), ((*31*)), 1791--1806. =end =begin HTML