Aqua planet experiment

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

1. Overview
2. Preparation of a directory for an experiment
3. Create initial data file
4. Create SST data file
5. Run the experiment
6. References

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

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

References

• 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.

DCPAM Development Group / GFD Dennou Staff