Mars experiment (Vertical 1-D experiment)

A method to perform a vertical 1-dimensional Mars experiment is described.

Following physical processes are used in this experiment.

  • Radiation (improved and/or simplified version of Takahashi et al., 2003; 2006)
  • Vertical diffusion (Mellor and Yamada, 1974, level 2)
  • Dry convective adjustment
  • Bucket model

Overview

This experiment requires data of distribution of surface properties and radiative property. So, the experiment is performed with the following 4 steps:

  • Preparation of directory for experiments
  • Preparation of initial condition
  • Preparation of additional data for experiments
  • Execution of experiments

Preparation of a directory for an experiment

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/Mars-1D-exp directory. Please create the directory and enter there as follows: 

$ mkdir -p ../dcpam5-exp/Mars-1D-exp
$ cd ../dcpam5-exp/Mars-1D-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/src/main/dcpam_init_data_surface  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.

Preparation of initial data

Create initial data file "init_1D_L36.nc" using "bin/dcpam_init_data", "conf/init_data_M_1D_L36.conf", and "bin/dcpam_init_data_surface", "conf/surface_data_M_1D.conf"

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

 *** MESSAGE [dcpam_init_data] ***  Run: Initial data generation
 *** MESSAGE [dcpam_init_data] ***  -- version = $Name:  $$Id: exp-Mars-1D.rd,v 1.3 2015/02/14 07:21:31 yot Exp $
 *** MESSAGE [namelist_util] ***  ----- Initialization Messages -----
 *** MESSAGE [namelist_util] ***    MaxNmlArySize = 256
 *** MESSAGE [namelist_util] ***  -- version = $Name:  $$Id: exp-Mars-1D.rd,v 1.3 2015/02/14 07:21:31 yot Exp $
 *** MESSAGE [timeset] ***  NAMELIST group "timeset_nml" is loaded from "init_data_M_1D_L36.conf".
&TIMESET_NML
 CAL_TYPE="noleap                          ",
                           :
 *** MESSAGE [restart_file_io] ***  Output::
 *** MESSAGE [restart_file_io] ***    OutputFile = init_1D_L36.nc
 *** MESSAGE [restart_file_io] ***    IntTime    = 365. [day]
 *** MESSAGE [restart_file_io] ***  -- version = $Name:  $$Id: exp-Mars-1D.rd,v 1.3 2015/02/14 07:21:31 yot Exp $
 *** MESSAGE [GDNcFileOpen] ***  "init_1D_L36.nc" is overwritten.
 *** MESSAGE [HistoryCreate1] ***  "init_1D_L36.nc" is created (origin=0.E+00 [sec])
 *** MESSAGE [HistoryClose] ***  "init_1D_L36.nc" is closed

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

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

 *** MESSAGE [dcpam_init_data_surface] ***  Run: SST data generation
 *** MESSAGE [dcpam_init_data_surface] ***  -- version = $Name:  $$Id: exp-Mars-1D.rd,v 1.3 2015/02/14 07:21:31 yot Exp $
 *** MESSAGE [namelist_util] ***  ----- Initialization Messages -----
 *** MESSAGE [namelist_util] ***    MaxNmlArySize = 256
 *** MESSAGE [namelist_util] ***  -- version = $Name:  $$Id: exp-Mars-1D.rd,v 1.3 2015/02/14 07:21:31 yot Exp $
 *** MESSAGE [timeset] ***  NAMELIST group "timeset_nml" is loaded from "surface_data_M_1D.conf".
&TIMESET_NML
 CAL_TYPE="noleap                          ",
                 :
 *** MESSAGE [GDNcFileOpen] ***  "surface_1D.nc" is overwritten.
 *** MESSAGE [HistoryCreate1] ***  "surface_1D.nc" is created (origin=0.E+00 [day])
 *** MESSAGE [HistoryClose] ***  "surface_1D.nc" is closed

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

Preparation of additional data

This experiment requires additional data, radiation parameter. Those data are prepared at following site.

Please download "lnactbl15m_2006-07-13.nc" from above site, and place those files at experimental directory.

Run the experiment

At the last, run a experiment using "dcpam_main" and "dcpam_M_1D_L36.conf". Resolution is L36. Integration time is 67 days. A restart data file, a restart data file of surface data, and some history data files are output. 

$ ./dcpam_main -N=dcpam_M_1D_L36.conf | tee Mars.log

If you want to change experimental settings, edit "dcpam_M_1D_L36.conf" in "conf" directory.

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 "&gtool_historyauto_nml".

$Id: exp-Mars-1D.rd,v 1.3 2015/02/14 07:21:31 yot Exp $