A method to perform an aqua planet experiment (planetary surface is covered by watar) proposed by Neale and Hoskins (2001) in an axisymmetric condition is described.
Following physical processes are used in this experiment.
This experiment uses an axisymmetric version of dynamical process. The axisymmetric version can be built by compiling the code with setting FFLAGS to -DAXISYMMETRY or -DAXISYMMETRY_SJPACK.
This experiment requires data of distribution of sea surface temperature. So, the experiment is performed with the following 4 steps:
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-nh01-exp-axisym" in the following example). Input 'Enter' to the question after that.
$ make ****** Setup a directory for a experiment ****** Enter directory name [testXX]: ape-nh01-exp-axisym Directory in which executable files are prepared [../src/main]: :
Move the directory
$ cd ape-nh01-exp-axisym
Note that you can perform an experiment in any directory by using executable files and configuration (namelist) files.
Create initial data file "init_T21L22.nc" using "init_data" and "init_data_ape_nh01_axisym_T21L22.nml"
$ ./init_data -N=init_data_ape_nh01_axisym_T21L22.nml *** MESSAGE [init_data] *** Run: Initial data generation *** MESSAGE [init_data] *** -- version = $Name: dcpam5-20110615 $$Id: exp-ape-nh01-axisym.rd,v 1.3 2011-02-25 06:53:37 yot Exp $ *** MESSAGE [namelist_util] *** ----- Initialization Messages ----- *** MESSAGE [namelist_util] *** MaxNmlArySize = 256 *** MESSAGE [namelist_util] *** -- version = $Name: dcpam5-20110615 $$Id: exp-ape-nh01-axisym.rd,v 1.3 2011-02-25 06:53:37 yot Exp $ *** MESSAGE [timeset] *** NAMELIST group "timeset_nml" is loaded from "init_data_ape_nh01_T21L22.nml". : *** MESSAGE [initial_data] *** PsAvr = 101080. *** MESSAGE [initial_data] *** QVapAvr = 0. *** MESSAGE [initial_data] *** Ueq = 0. *** MESSAGE [initial_data] *** *** MESSAGE [initial_data] *** -- version = $Name: dcpam5-20110615 $$Id: exp-ape-nh01-axisym.rd,v 1.3 2011-02-25 06:53:37 yot Exp $ *** MESSAGE [HistoryClose] *** "init_T21L22.nc" is closed ############## CPU TIME SUMMARY ################ others 0.400020E-01 ------------------------------------------------ TOTAL TIME = 0.400020E-01
This experiment requires additional data of ozone distribution. Those data are prepared at following site.
Please download "O3_NH01_lon1_T021.nc" from above site, and place those files at experimental directory.
Next, create SST data "sst_T21.nc" using "sst_data" and "sst_data_ape_nh01_axisym_T21.nml".
$ ./sst_data -N=sst_data_ape_nh01_axisym_T21.nml *** MESSAGE [sst_data] *** Run: SST data generation *** MESSAGE [sst_data] *** -- version = $Name: dcpam5-20110615 $$Id: exp-ape-nh01-axisym.rd,v 1.3 2011-02-25 06:53:37 yot Exp $ *** MESSAGE [namelist_util] *** ----- Initialization Messages ----- *** MESSAGE [namelist_util] *** MaxNmlArySize = 256 *** MESSAGE [namelist_util] *** -- version = $Name: dcpam5-20110615 $$Id: exp-ape-nh01-axisym.rd,v 1.3 2011-02-25 06:53:37 yot Exp $ *** MESSAGE [timeset] *** NAMELIST group "timeset_nml" is loaded from "sst_data_ape_nh01_axisym_T21.nml". : *** MESSAGE [surface_data] *** SoilHeatCap = 2100000. *** MESSAGE [surface_data] *** SoilHeatDiffCoef = 1.2 *** MESSAGE [surface_data] *** -- version = $Name: dcpam5-20110615 $$Id: exp-ape-nh01-axisym.rd,v 1.3 2011-02-25 06:53:37 yot Exp $ *** MESSAGE [HistoryClose] *** "sst_T21.nc" is closed ############## CPU TIME SUMMARY ################ others 0.400000E-02 ------------------------------------------------ TOTAL TIME = 0.400000E-02
At the last, run a experiment using "dcpam_main" and "dcpam_ape_nh01_axisym_T21L22.nml". Resolution is T21L22, 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_nh01_axisym_T21L22.nml | tee ape-nh01-axisym.log
If initial data, surface data, settings of experiments want to be changed, edit "init_data_ape_nh01_axisym_T21L22.nml", "sst_data_ape_nh01_axisym_T21.nml", "dcpam_ape_nh01_axisym_T21L22.nml".