Held and Suarez (1994) dynamcal core experiment
A method to perform dynamical core experiment proposed by Held and Suarez (1994) is described.
Following physical processes are used in this experiment.
- Newtonian relaxation of the temperature field to a zonally symmetric state
- Rayleigh damping of low-level winds
- processes related to moisture are excluded.
- Dry convective adjustment
Overview
This experiment does not require additional data, such as SST data. So, the experiment is performed with the following 3 steps:
- Preparation of directory for experiments
- Preparation of initial condition
- Execution of experiments
Preparation of a directory for experiments
Let us move to the top directory of dcpam5 src tree (assuming dcpam5-YYYYMMDD in the following example). Here, we perform an experiment in ../dcpam5-exp/hs94-exp directory. Please create the directory and enter there as follows:
$ mkdir -p ../dcpam5-exp/hs94-exp $ cd ../dcpam5-exp/hs94-exp
Then, pleaase create the directories for executable files and configuration files as follows:
$ mkdir bin $ mkdir conf
Finally, executable files and configuration files are copied as follows:
$ cp ../../dcpam5-YYYYMMDD/src/main/dcpam_main bin $ cp ../../dcpam5-YYYYMMDD/src/main/dcpam_init_data bin $ cp ../../dcpam5-YYYYMMDD/exp_setup_files/*.conf conf
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_T21L20.nc" using "dcpam_init_data" and "init_data_hs94_T21L20.conf"
$ ./bin/dcpam_init_data -N=./conf/init_data_hs94_T21L20.conf
:
*** 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_T21L20.nc
*** MESSAGE [restart_file_io] *** IntTime = 1. [day]
*** MESSAGE [HistoryCreate3] *** "init_T21L20.nc" is created (origin=0. [day])
*** MESSAGE [initial_data] *** NAMELIST group "initial_data_nml" is loaded from "./conf/init_data_T21L16.conf".
*** 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_T21L20.nc" is closed
Run the experiment
Using an executable files 'dcpam_main' and a NAMELIST file dcpam_hs94_T21L20.conf, execute 'dcpam_main' as follows. This program will be finished in few minutes - tens of minutes.
$ ./bin/dcpam_main -N=./conf/dcpam_hs94_T21L20.conf | tee hs94.log
*** MESSAGE [dcpam_main] *** Run: Held and Suarez (1994) benchmark test
:
*** MESSAGE [timeset] *** ----- Initialization Messages -----
*** MESSAGE [timeset] *** StartTime = 0. [day]
*** MESSAGE [timeset] *** EndTime = 10. [day]
*** MESSAGE [timeset] *** DelTime = 30. [min]
*** MESSAGE [timeset] *** = 1800. [sec]
:
:
*** MESSAGE [fileset] *** ----- Initialization Messages -----
*** MESSAGE [fileset] *** FileTitle = Held and Suarez (1994) benchmark test
*** MESSAGE [fileset] *** FileSource = dcpam5 (http://www.gfd-dennou.org/library/dcpam)
*** MESSAGE [fileset] *** FileInstitution = GFD Dennou Club (http://www.gfd-dennou.org)
*** MESSAGE [fileset] *** -- version = ....
*** MESSAGE [gridset] *** ----- Initialization Messages -----
*** MESSAGE [gridset] *** nmax = 21
*** MESSAGE [gridset] *** imax = 64
*** MESSAGE [gridset] *** jmax = 32
*** MESSAGE [gridset] *** kmax = 20
*** MESSAGE [gridset] *** -- version = ....
:
########## PREDICTION OF CALCULATION ###########
Start Date 2009-03-06T08:26:41+09:00
Current Date 2009-03-06T08:26:45+09:00
Progress 20.00% [***** ]
Remaining CPU TIME 0.160000E+02
Completion Date 2009-03-06T08:27:01+09:00
:
############## CPU TIME SUMMARY ################
held_suarez_1994 0.165611E+01
dynamics_hspl_vas83 0.191452E+02
timefilter_asselin19 0.264027E+00
others 0.896038E+00
------------------------------------------------
TOTAL TIME = 0.219614E+02
In this case, 10 days integration is performed. Resolution is T21L20 (Number of grid points of longitude, latitude is 64, 32, and number of vertical layer is 20). History data are output to 'U.nc' and 'Temp.nc' etc., and a restart data is output to 'restart.nc'.
Visualization of result
Zonal mean temperature at 10 day is shown in following figures.
This figure is drawn by using Gphys, as follows:
Please type irb.
$ irb
Following prompt of irb is shown.
irb(main):001:0>
Type following lines. Please note that the left-most number is line number, and it has not been typed.
1: require "numru/ggraph"
2: include NumRu
3: gphys = GPhys::IO.open('Temp.nc', 'Temp')
4: DCL.gropn(1)
5: DCL.sgpset('lcntl', false) ; DCL.uzfact(0.7)
6: GGraph.tone( gphys.mean('lon').cut('time'=>11) )
7: GGraph.color_bar
8: DCL.grcls
Type quit to exit irb.
In Held and Suarez (1994), integration period is 1200 days and resolution is T63L20 (Number of grid points of longitude, latitude is 192, 96, and number of vertical layer is 20). The result of such integration by using the dcpam is shown here (Japanese page only).
Final remark
If you would like to keep the results of the numerical experiment, please create a sub directory ("data01" for example), move the files, and keep the current directory neatly.
$ mkdir data01 $ mv *.nc *.log data01/
You can output the results of the numerical experimet concurrently with the execution. Please refer to "Change output settings" in Change experimental setup and try to configure "FilePrefix" variable in NAMELIST ">ool_historyauto_nml".
References
$Id: exp-hs94.rd,v 1.14 2013-10-06 13:46:12 yot Exp $