Coupling MOZART-4 with Meso-NH for initial and boundary conditions
The following procedure is a beta version of the coupling of Meso-NH with the MOZART-4 CTM. The fortran routines are compatible with the MASDEV_4_10 and MASTER4_9 versions.
The MOZART-4 CTM model can be used to prescribe chemical initial and boundary conditions for Meso-NH. A full description of the MOZART-4 model is given in Emmons et al. (2010) (http://www.geosci-model-dev.net/3/43/2010/gmd-3-43-2010.html)
STEP 1: download the MOZART-4 simulations from http://www.acom.ucar.edu/wrf-chem/mozart.shtml
DOWNLOAD ONLY ONE DAY PER FILE!
Each daily file contains 4 time slots at 0000, 0600, 1200 and 1800 UTC
MOZART-4 simulations use the meteorological model GEOS-5 and are available from Jan 1, 2007 to Dec 31, 2016 (1.9x2.5 deg and 56 levels)
Change the name of the MOZART file (the default name is too long for Meso-NH)
> cp mozart4-20120206041851846590.nc mozart4-20090101.nc
STEP 2: create an user library (see the A-INSTALL file) and compile the three routines:
read_chem_data_netcdf_case.f90
STEP 3: Run PREP_REAL_CASE
Specify the type of the chemical file: HCHEMFILETYPE='NETCDF'
Specify the hour (UTC) you need from the MOZART-4 daily file in CDUMMY1 (CDUMMY1="00" or "06" or "12" or "18").
- Be aware that there is no checking of the coherency between the MOZART date and hour and those from the meteorological model used to initialized Meso-NH (ARPEGE, ECMWF, ...).
Specify the type of MOZART-4 simulation (CDUMMY2="NCEP" or CDUMMY2="GEOS5")
Specify the correspondence between the Meso-NH chemical species (including aerosol if necessary) and the MOZART chemical species in an ascii file MOZ1.nam.
Example of PRE_REAL1.nam:
&NAM_FILE_NAMES HATMFILE ='ecmwf.O3RIVER.20040810.18', HATMFILETYPE='GRIBEX',
- HCHEMFILE='mozart4-20040810.nc', HCHEMFILETYPE='NETCDF', HPGDFILE ='ICARTT1008_PGD_60km' , CINIFILE='CPLCH20040810.60' /
&NAM_VER_GRID NKMAX=50, YZGRID_TYPE='FUNCTN',
- ZDZGRD=60., ZDZTOP=650., ZZMAX_STRGRD=2000., ZSTRGRD=4., ZSTRTOP=6. /
&NAM_BLANK CDUMMY1="18", CDUMMY2="NCEP" /
The file MOZ1.nam must be pasted at the end of the PRE_REAL1.nam. The namelist NAM_AERO_CONF must be added for activating the ORILAM module.
Example of MOZ1.nam:
MOZ1.nam for a simulation with ReLACS3 and ORILAM-SOA module activated for aerosols: MOZ1.nam
News of the current version:
- The coupling code is compatible whatever the chemical scheme used with the model.
- Aerosols are accounted for in the new version as MOZART-4 provides 3D fields for sulfate, ammonium, nitrate, SOA, black carbon (hydrophobic, hydrophylic), organic carbon (hydrophobic, hydrophylic), sea salt (4 bins) and dust (5 bins).
Illustrations
Source: S. Strada