Density Current Experiment by Straka et al. (1993)

A method to perform a density current experiment by Straka et al. (1993) is described.

Overview

This 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 deepconv/arare5 src tree (assuming arare5-YYYYMMDD in the following example). Here, we perform an experiment in ../arare5-exp/s93-exp directory. Please create the directory and enter there as follows:

$ mkdir -p ../deepconv-exp/s93-exp
$ cd ../deepconv-exp/s93-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 ../../arare5-YYYYMMDD/src/main/arare bin
$ cp ../../arare5-YYYYMMDD/src/main/arare_init-data bin/
$ cp ../../arare5-YYYYMMDD/exp_setup_files/04_density-current/density-current_512x1x64_Center4_Center2_PTemp.nc 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 "density-current_512x1x64_Center4_Center2_ini.nc" using "bin/arare_init-data" and density-current_512x1x64_Center4_Center2.conf.

$ ./bin/arare_init-data -N=conf/density-current_512x1x64_Center4_Center2.conf

 *** MESSAGE [main] ***  Namelist file is 'density-current_512x1x64_Center4_Center2.conf'
 *** MESSAGE [argset_init] ***  NAMELIST FILE = density-current_512x1x64_Center4_Center2.conf
 *** MESSAGE [namelist_util] ***  ----- Initialization Messages -----
 *** MESSAGE [namelist_util] ***    MaxNmlArySize = 256
 *** MESSAGE [namelist_util] ***  -- version = $Name:  $$Id: exp-s93.rd,v 1.3 2016/02/13 04:42:06 odakker Exp $
 *** MESSAGE [gridset_init] ***  xsub = 1
                       : 
 *** MESSAGE [main] ***  Making Initial data....
 *** MESSAGE [main] ***  Making Initial data (basic) named DRY...
 *** MESSAGE [main] ***  Making Initial data (disturb) named CosXZ...
 *** MESSAGE [main] ***  Output variables into netCDF file...
 *** MESSAGE [restartfileioIO_init] ***  InitialFile  = density-current_512x1x64_Center4_Center2_ini.nc
 *** MESSAGE [restartfileioIO_init] ***  InputFile  =
 *** MESSAGE [restartfileioIO_init] ***  OutputFile = density-current_512x1x64_Center4_Center2_res.nc
 *** MESSAGE [restartfileioIO_init] ***  INPUT  =
 *** MESSAGE [restartfileioIO_init] ***  OUTPUT = density-current_512x1x64_Center4_Center2_ini.nc
 *** MESSAGE [GDNcFileOpen] ***  "density-current_512x1x64_Center4_Center2_ini.nc" is overwritten.

Run the experiment

Using an executable files 'arare' and a NAMELIST file density-current_512x1x64_Center4_Center2.conf, execute 'arare' as follows. This program will be finished in few minutes - tens of minutes.

$ ./bin/arare -N=conf/density-current_512x1x64_Center4_Center2.conf | tee s93.log

 *** MESSAGE [main] ***  Namelist file is 'conf/density-current_512x1x64_Center4_Center2.conf'
 *** MESSAGE [argset_init] ***  NAMELIST FILE = conf/density-current_512x1x64_Center4_Center2.conf
 *** MESSAGE [namelist_util] ***  ----- Initialization Messages -----
 *** MESSAGE [namelist_util] ***    MaxNmlArySize = 256
 *** MESSAGE [namelist_util] ***  -- version = $Name:  $$Id: exp-s93.rd,v 1.3 2016/02/13 04:42:06 odakker Exp $
 *** MESSAGE [timeset_init] ***  DelTimeLong  = 1.
                   :
  *** MESSAGE [HistoryClose] ***  "density-current_512x1x64_Center4_Center2_PTemp.nc" is closed

  ############## CPU TIME SUMMARY ################
  initialization         0.176000E+00
  time-integration       0.672216E+03  (11.20 minutes)
  ------------------------------------------------
         TOTAL TIME =    0.672392E+03  (11.21 minutes)

In this case, about 10 minites integration is performed. History data are output to 'VelX.nc' and 'PTemp.nc' etc., and a restart data is output to 'density-current_512x1x64_Center4_Center2_res.nc'.

Visualization

Please see First step analysis and visualization.

References

• Straka, J. M., Wilhelmson, R. B., Wicker, A. L., Anderson, J. R., Droegemeier, K. K., 1993 : Numerical solutions of a non-linear density current : A benchmark solution and comparisons. International Journal for Numerical Methods in Fluids , 17, 1--22

---

$Id: exp-s93.rd,v 1.3 2016/02/13 04:42:06 odakker Exp $