Experiment for Venus-like planet

A method to perform an experiment for Venus-like planet by using forcing by Yamamoto and Takahashi (2003) is described.

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

Overview

This experiment does not require additional 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

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

$ cd practice
$ mkdir Venuslike-exp
$ cd Venuslike-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_T10L50.nc" using "init_data" and "init_data_Vlp_T10L50.nml"

$ ./init_data -N=init_data_Vlp_T10L50.nml

 *** MESSAGE [init_data] ***  Run: Initial data generation
 *** MESSAGE [init_data] ***  -- version = $Name: dcpam5-20120220 $$Id: exp-Venuslike.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-Venuslike.rd,v 1.3 2012-02-01 06:27:27 yot Exp $
 *** WARNING [timeset] ***  NAMELIST group "timeset_nml" is not found in "init_data_Vlp_T10L50.nml" (iostat=186).
 *** MESSAGE [timeset] ***  ----- Initialization Messages -----
                       : 
 *** MESSAGE [initial_data] ***  ----- Initialization Messages -----
 *** MESSAGE [initial_data] ***    Pattern = SmallVapAvr = 0.
 *** MESSAGE [initial_data] ***    Ueq     = 0.
 *** MESSAGE [initial_data] ***
 *** MESSAGE [initial_data] ***  -- version = $Name: dcpam5-20120220 $$Id: exp-Venuslike.rd,v 1.3 2012-02-01 06:27:27 yot Exp $

Run the experiment

Using an executable files 'dcpam_main' and a NAMELIST file dcpam_Vlp_T10L50.nml, execute 'dcpam_main' as follows. This program will be finished in few minutes - tens of minutes.

$ ./dcpam_main -N=dcpam_Vlp_T10L50.nml | tee Vlp.log

 *** MESSAGE [namelist_util] ***  ----- Initialization Messages -----
 *** MESSAGE [namelist_util] ***    MaxNmlArySize = 256
 *** MESSAGE [namelist_util] ***  -- version = $Name: dcpam5-20120220 $$Id: exp-Venuslike.rd,v 1.3 2012-02-01 06:27:27 yot Exp $
 *** MESSAGE [dcpam_main] ***  NAMELIST group "dcpam_main_nml" is loaded from "dcpam_Vlp_T10L50.nml".
 &DCPAM_MAIN_NML FLAGDYNAMICS=T,FLAGFULLPHYSICS=F,FLAGHS94=F,FLAGVENUSSIMPLE=T,FLAGJUPITERSIMPLE=F,FLAGRADIATIONDENNOUAGCM=T,FLAGRADIATIONDCPAMEV1=F,FLAGRADIATIONDCPAMMV1=F,FLAGRADIATIONSL09=F,FLAGPHYIMP1LAYMODEL=T,FLAGPHYIMPSOILMODEL=F,FLAGPHYIMPSOILMODELSO=F,FLAGPHYIMPATMONLY=F,FLAGVERTICALFILTER=F/
 *** MESSAGE [dcpam_main] ***  FlagDynamics=<T>.
 *** MESSAGE [dcpam_main] ***  FlagRadiationDennouAGCM=<T>.
 *** MESSAGE [dcpam_main] ***  FlagRadiationDcpamEV1  =<F>.
 *** MESSAGE [dcpam_main] ***  FlagRadiationDcpamMV1  =<F>.
                   :
 *** MESSAGE [HistoryClose] ***  "VUBalance.nc" is closed

 ############## CPU TIME SUMMARY ################
 surface_properties     0.000000E+00
 auxiliary              0.108007E+00
 vdiffusion_my1974      0.520025E-01
 dynamics_hspl_vas83    0.187212E+01
 mass_fixer             0.280015E-01
 check_prog_vars        0.800041E-02
 timefilter_asselin19   0.280015E-01
 others                 0.552037E+00
 ------------------------------------------------
        TOTAL TIME =    0.264816E+01

In this case, 5 days integration is performed. Resolution is T10L50 (Number of grid points of longitude, latitude is 32, 16, and number of vertical layer is 50). 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 distribution after integration for 5 days is shown below.

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.set_fig( 'itr'=> 2 )
7: GGraph.tone( gphys.mean('lon').cut('time'=>5) )
8: GGraph.color_bar
9: DCL.grcls

Type quit to exit irb. 

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

References

• Yamamoto, M, and M. Takahashi, 2003: The fully developed superrotation simulated by a general circulation model of a Venus-like atmosphere, J. Atmos. Sci., 60, 561--574.

DCPAM Development Group / GFD Dennou Staff