[ English | Japanese ] [ GFD Dennou Club / DCPAM ] [ dcpam5 Documents / GOKURAKU dcpam5 ]

Aqua planet experiment by Neale and Hoskins (2001)

A method to perform an aqua planet experiment (planetary surface is covered by watar) proposed by Neale and Hoskins (2001) is described.

Following physical processes are used in this experiment.


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

Overview

This experiment requires data of distribution of sea surface temperature. So, the experiment is performed with the following 4 steps:

Preparation of a directory for an experiment

Here, we perform an experiment in practice/ape-nh01-exp directory. Please create the directory and enter there as follows:

$ cd practice
$ mkdir ape-nh01-exp
$ cd ape-nh01-exp

Then, executable files and configuration files are copied as follows:

$ cp ../../src/main/dcpam_main ../../src/main/init_data ../../src/main/sst_data .
$ cp ../../src/sample_nml/*.nml .

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_T21L22.nc" using "init_data" and "init_data_ape_nh01_T21L22.nml"

$ ./init_data -N=init_data_ape_nh01_T21L22.nml

 *** MESSAGE [init_data] ***  Run: Initial data generation
 *** MESSAGE [init_data] ***  -- version = $Name: dcpam5-20120220 $$Id: exp-ape-nh01.rd,v 1.3 2012-02-01 06:27:27 yot Exp $
 *** MESSAGE [namelist_util] ***  ----- Initialization Messages -----
 *** MESSAGE [namelist_util] ***    MaxNmlArySize = 256
 *** MESSAGE [namelist_util] ***  -- version = $Name: dcpam5-20120220 $$Id: exp-ape-nh01.rd,v 1.3 2012-02-01 06:27:27 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-20120220 $$Id: exp-ape-nh01.rd,v 1.3 2012-02-01 06:27:27 yot Exp $
 *** MESSAGE [HistoryClose] ***  "init_T21L22.nc" is closed

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

Preparation of additional data

This experiment requires additional data of ozone distribution. Those data are prepared at following site.

Please download "O3_NH01_T021.nc" from above site, and place those files at experimental directory.

Create SST data file

Next, create SST data "sst_T21.nc" using "sst_data" and "sst_data_ape_nh01_T21.nml".

$ ./sst_data -N=sst_data_ape_nh01_T21.nml

 *** MESSAGE [sst_data] ***  Run: SST data generation
 *** MESSAGE [sst_data] ***  -- version = $Name: dcpam5-20120220 $$Id: exp-ape-nh01.rd,v 1.3 2012-02-01 06:27:27 yot Exp $
 *** MESSAGE [namelist_util] ***  ----- Initialization Messages -----
 *** MESSAGE [namelist_util] ***    MaxNmlArySize = 256
 *** MESSAGE [namelist_util] ***  -- version = $Name: dcpam5-20120220 $$Id: exp-ape-nh01.rd,v 1.3 2012-02-01 06:27:27 yot Exp $
 *** MESSAGE [timeset] ***  NAMELIST group "timeset_nml" is loaded from "sst_data_ape_nh01_T21.nml".
                       :
 *** MESSAGE [surface_data] ***    SoilHeatCap      = 2100000.
 *** MESSAGE [surface_data] ***    SoilHeatDiffCoef = 1.2
 *** MESSAGE [surface_data] ***  -- version = $Name: dcpam5-20120220 $$Id: exp-ape-nh01.rd,v 1.3 2012-02-01 06:27:27 yot Exp $
 *** MESSAGE [HistoryClose] ***  "sst_T21.nc" is closed

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

Run the experiment

At the last, run a experiment using "dcpam_main" and "dcpam_ape_nh01_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_T21L22.nml | tee ape-nh01.log

If initial data, surface data, settings of experiments want to be changed, edit "init_data_ape_nh01_T21L22.nml", "sst_data_ape_nh01_T21.nml", "dcpam_ape_nh01_T21L22.nml".

Visualization of result

Zonal mean temperature distribution after integration for 2 days is shown below.

images/ape_nh01_irb_temp_meanlon_time02day.png

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'=>2) )
7: GGraph.color_bar
8: DCL.grcls

Type quit to exit irb. 

Above figure shows spin-up phase of the atmosphere, because the Earth's atmosphere does not reach quasi-equilibrium in 5 days.

References


$Id: exp-ape-nh01.rd,v 1.3 2012-02-01 06:27:27 yot Exp $
DCPAM Development Group / GFD Dennou Staff dcstaff@gfd-dennou.org