GEOS-Chem nested grid simulations

From Geos-chem
Revision as of 13:12, 13 May 2008 by Bmy (Talk | contribs) (Nested Grid Frequently Asked Questions)

Jump to: navigation, search

Nested Grid Frequently Asked Questions

Q. Can I run a nested grid for GEOS-5 in v8-01-01?

A. Not yet. However, Yuxuan Wang and Dan Chen have been working on this capability and we anticipate that it will be ready in a near future release of GEOS-Chem.


Q. Which nested grids are available for me to use?

A. At present, you can use the following nested grids:


Dylan Millet (dbm@umn.edu) wrote:

Do you know if the 0.5x0.667 nested simulation that's being developed uses a 4x5 or a 2x2.5 grid outside of the nested domain?

Yuxuan Wang (wang3@fas.harvard.edu) replied

The nested simulation uses 4x5 simulation results only at the lateral boundaries of the nested domain. For example, the East Asia nested domain covers 10S to 55N and 70E to 150E, one needs to save 4x5 simulation results (tracer concentrations only) for each vertical layer, at a frequency of every 3 hours, at the following four boundaries:
lat = 10S, lon = 70E
lat = 10S, lon = 150E
lat = 55N, lon = 70E
lat = 55N, lon = 150E
One can easily change the code to have the nested-grid simulation use 2x2.5 simulation as boundary conditions, but the default has been to take 4x5 results since it's much faster to run a 4x5 simulation.

--Bob 09:12, 13 May 2008 (EDT)

Setting up a China nested-grid simulation

Here is a check list and input.geos file that you need to run the nested-grid China simulation (supplied by Yuxuan Wang):

define.h

These are the settings you need to make in the GEOS-Chem define.h file:

#define GEOS_3      'GEOS_3'
#define GRID1x1     'GRID1x1'
#define NESTED_CH   'NESTED_CH'
#define GRID30LEV   'GRID30LEV'

input.geos

These are the settings you need to make in the GEOS-Chem input.geos file:

Saving out 4 x 5 boundary conditions

You will have to run the global model at 4 x 5 resolution to save out the transport boundary conditions. Here are the relevant settings in input.geos:

%%% SIMULATION MENU %%% :
Start YYYYMMDD, HHMMSS  : 20010101 000000
End   YYYYMMDD, HHMMSS  : 20020101 000000
Run directory           : ./
Input restart file      : restart.YYYYMMDDhh
Make new restart file?  : F
Output restart file(s)  : restart.YYYYMMDDhh
Root data directory     : /data/ctm/GEOS_4x5/
 => GCAP       subdir   : AGRID/YYYY/MM/
 => GEOS-3     subdir   : GEOS_3/YYYY/MM/
 => GEOS-4     subdir   : GEOS_4_v4/YYYY/MM/
 => GEOS-5     subdir   : GEOS_5/YYYY/MM/
Dir w/ 1x1 emissions etc: /data/ctm/GEOS_1x1/
Temporary directory     : TEMP/
Unzip met fields?       : T
Wait for met fields?    : T
Use variable tropopause : N
Global offsets I0, J0   : 0 0

 ... 

%%% TRANSPORT MENU %%%  :
Turn on Transport       : T
 => Use Flux Correction?: F
 => Fill Negative Values: T 
 => IORD, JORD, KORD    : 3  3  7
Transport Timestep [min]: 30  
Use strat O3/NOy BC's   : T

...

%%% CONVECTION MENU %%% :
Turn on Cloud Conv?     : T
Turn on PBL Mixing?     : T
Convect Timestep (min)  : 30

...

%%% NESTED GRID MENU %%%:  
Save TPCORE 4x5 BC's    : T
TPCORE 4x5 BC directory : /as/home/ctm/yxw/bc/bc_ozone_std/
4x5 BC timestep [min]   : 180
LL box of 4x5 BC region :  51  21
UR box of 4x5 BC region :  67  37
1x1 offsets I0_W, J0_W  :   3   3

1x1 nested grid simulation

Then once you have saved out the 4 x 5 boundary conditions you can do the 1 x 1 nested grid simulation. Here are the appropriate settings in input.geos:

%%% SIMULATION MENU %%% :
Start YYYYMMDD, HHMMSS  : 20010301 000000
End   YYYYMMDD, HHMMSS  : 20010302 000000
Run directory           : ./
Input restart file      : restart.YYYYMMDDhh
Make new restart file?  : F
Output restart file(s)  : restart.YYYYMMDDhh
Root data directory     : /data/geos/GEOS_1x1_CH/
 => GCAP       subdir   : AGRID/YYYY/MM/
 => GEOS-3     subdir   : GEOS_3/YYYY/MM/
 => GEOS-4     subdir   : GEOS_4_v4/YYYY/MM/
 => GEOS-5     subdir   : GEOS_5/YYYY/MM/
Dir w/ 1x1 emissions etc: /data/ctm/GEOS_1x1/
Temporary directory     : TEMP/
Unzip met fields?       : T
Wait for met fields?    : T
Use variable tropopause : F
Global offsets I0, J0   : 250 79       # need to be changed for NA; index refer to 1x1 resolution

 ...

%%% TRANSPORT MENU %%%  :
Turn on Transport       : T
 => Use Flux Correction?: F
 => Fill Negative Values: T 
 => IORD, JORD, KORD    : 3  3  7
Transport Timestep [min]: 10  
Use strat O3/NOy BC's   : T

...

%%% CONVECTION MENU %%% :
Turn on Cloud Conv?     : T
Turn on PBL Mixing?     : T
Convect Timestep (min)  : 10

... 

%%% NESTED GRID MENU %%% 
Save TPCORE 4x5 BC's    : T
TPCORE 4x5 BC directory : /as/home/strat/yxw/bc/  # change to be your own directory
4x5 BC timestep [min]   : 180
LL box of 4x5 BC region :  51  21                 # need to be changed for NA; index refer to 4x5 resolution
UR box of 4x5 BC region :  67  37                 # need to be changed for NA; index refer to 4x5 resolution
1x1 offsets I0_W, J0_W  :   3   3                 # do not change for NA

NOTES

General

  1. If not NESTED_CH (East Asia window) or NESTED_NA (north America window), need to first obtain
    • The 1x1 meteorological field for the nested domain;
    • Tll the input files (mostly emissions) at 1x1 resolution for the nested domain.
  2. There needs to be two run directories: one for the 4x5 global run and one for the 1x1 nested run
  3. Only one code directory is necessary. Code version v7-01-02 and higher can be compiled for both the global run and the nested run
  4. Run the 4x5 global simulation first. With correct switches in the 4x5 input.geos file, boundary conditions will
    • be saved out automatically into the designated directory and will be read later by the nested run.
    • The code automatically takes care of formatting.

Code

  1. Change grid size in define.h
  2. Change scaling factor in lightning_nox_mod.f (if not Chinese window)
    • search for 1x1 in the file to locate the lines of code where changes should be made
  3. Change IGLOB, JGLOB in CMN_SIZE (if not Chinese window)
    • IGLOB and JGLOB refers to the size of the nested window measured at 1x1 resolution
  4. Change PARAMETER( NLAND=2681, NPULSE=3 ) in commsoil.h (if not Chinese window)
    • changes are already made in the code; just uncomment the NA part
  5. Change datadir for O3 column in toms_mod.f --- if not CHINA or NA windows
    • need to create TOMS_O3col_2001.geos file at 1x1 resolution for the nested region
    • use the IDL code: regrid_4to1_O3_column.pro (note that line 33 and 34 needs to be modified for NA)