|
|
(93 intermediate revisions by 7 users not shown) |
Line 1: |
Line 1: |
− | == Nested Grid Frequently Asked Questions ==
| + | This content has been moved to our [https://geos-chem.readthedocs.io/en/latest/supplemental-guides/nested-grid-guide.html '''Run nested grid simulations'''] Supplemental Guide at [https://geos-chem.readthedocs.io <tt>geos-chem.readthedocs.io</tt>]. |
− | | + | |
− | 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:
| + | |
− | | + | |
− | * [http://www.as.harvard.edu:16080/ctm/geos/doc/man_test/appendix_2.html#A2.5.1 GEOS-3 1x1, China nested grid]
| + | |
− | * [http://www.as.harvard.edu:16080/ctm/geos/doc/man_test/appendix_2.html#A2.5.2 GEOS-3 1x1, North America nested grid]
| + | |
− | | + | |
− | | + | |
− | == 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 <tt>define.h</tt> 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 <tt>input.geos</tt> 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 ====
| + | |
− | | + | |
− | # 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.
| + | |
− | # There needs to be two run directories: one for the 4x5 global run and one for the 1x1 nested run
| + | |
− | # 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
| + | |
− | # 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 ====
| + | |
− | # Change grid size in define.h
| + | |
− | # 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
| + | |
− | # 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
| + | |
− | # 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
| + | |
− | # 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)
| + | |