CHIMERE Training Course - LMD - Ecole Polytechnique

CHIMERE Training Course General CHIMERE Structure Dmitry Khvorostyanov, Laurent Menut, Myrto Valari Laboratoire de Météorologie Dynamique IPSL Ecole Polytechnique, 91128 Palaiseau Cedex Contact: [email protected]

CHIMERE Training Course, D. Khvorostyanov, L. Menut, M. Valari

ChemistryTransport Model (CTM)

Mass outflow Mass inflow


CHIMERE CORE


CHIMERE Core Routines


CHIMERE Preprocessing •

scripts/makechemistry.sh – the chemical preprocessor chemprep/

•

chimeredomain.sh – to create domainspecific data in domains/ :

–

Domain coordinates: 

domains/HCOORD/ COORD_CONT5



domains/VCOORD/VCOORD_8_995_500

–

Landuse data : domains/CONT5/LANDUSE_GLCF_CONT5

–

Parameters for MEGAN biogenic emissions (LAI, emission factors for each species): EFMAP_LAI_CONT5.nc

–

Dust emissions input parameters (soil properties, z0, erodibility)

→ makechemistry results in chemprep/inputdata.xxx/ directory → chimeredomains results in domains/CONT5/ directory

Specified in chimere.par: # CHIMERE Simulation Domain [$nested] Nested run? (yes/no) : no, yes, no [$dom] CHIMERE domain : CONT5, FRAGN15 [$nzdoms] Number of parallel zonal subdomains : 4 [$nmdoms] Number of parallel merid subdomains : 4 [$nlevels] Number of vertical layers : 8, 15 [$top_chimere_pressure] Top layer pressure : 500 [$first_layer_pressure] First layer pressure : 997


Chemical Preprocessor: chemprep/

•

Creates the directory chemprep/inputdata.211001021.10/ where '211001021.10' correspond to parameters specified in the parameter file chimere.par

•

Calculates aerosol size distribution (src/chemprepdistrib.f)

•

Defines families (src/chemprepfamilies.f)

•

The resulting inputdata.211001021.10/ contains data files for running CHIMERE: parameters and constants for gas & aerosol chemistry, depositions, etc.

Parameters specified in chimere.par: # Chemistry Options [$mecachim] Chemistry mechanism (0..2) :2 [$aero] Chemically-active aerosols (1/0) :1 [$nbins] Number of aerosol size sections :9 [$seasalt] Include sea salts? (1/0) :0 [$dust] Include dust? (1/0) :1 [$carb] Allow carbon speciation? :1 [$soatyp] Secondary organic aerosol scheme (0..2) : 2 [$gtrc] Include gasesous tracers? :0 [$ptrc] Include particulate tracers? :0 [$nequil] ISORROPIA coupling (1) / tabulation (0) : 0 [$npeq] - equilibrium calc freq (phys steps, < 5) :1


data/ File Name

Parameters

Usage

AEROSOL

Aerosol size distribution, type, molar mass, origin (bio/anthropogenic)

Aerosol emissions and transport

Names of emitted anthropogenic species

Anthropogenic emissions & transport

Names of emitted biogenic species

Biogenic emissions & transport

DEPO_PARS

Dry deposition parameters depending on vegetation type

Dry deposition

DEPO_SPEC

Speciesdependent parameters to calculate resistances for dry deposition

Dry deposition initio/inidepo.f90

WETD_SPEC

Speciesdependent wet deposition parameters (wetdk)

Wet deposition (inidepo.f90 → depwet.f90)

SEA

Coefficients for exchange between air & sea water (dKh0, dAh, dTh, ksea)

Exchanges between air and sea water (inidepo.f90 → locvalues → depvel)

SOIL

Coefficients for exchange between air & soil (koc, kdoc, diffair, diffwat, rdiffwat, ksoil)

Exchanges between air & soil (inidepo.f90 → depvel.f90)

VEGETATION

Coefficients for exchange between air & vegetation (dKoa0, dAoa, dToa, dmmm, kveget)

Exchanges between air & vegetation (inidepo.f90 → locvalues → air2veget.f90

FAMILIES

Families and their species

Output (outprint.f90)

LANDPAR

Seasonal cycle of momentum roughness height, depending on landuse type

Surface layer parameters diagnostics: mean z0 , SHF, urban T & U correction (diagmet_science.f90)

MAKEPRI

Used to generate PRIMARY file

See below (inputdata.xxx/ dir)

Description of gas chemical reactions used in CHIMERE

Model gasphase chemistry

Reaction rate constants

Reaction rate calculation for gases (rates.f90)

SEMIVOL

Parameters for semivolatile species

Absorption of semivolatile species

TRACER

Tracer sources characteristics

Tracer transport

Default hourly tracer emission profile (one tracer per line, one hour per column)

Tracer transport

ANTHROPIC BIOGENIC

REACTIONS REACTION_RATES

TRACER_PROFILE


AEROSOL # Number of sections D(2.5) D(10.0) 8 6 8 0.3906250E07 0.7812500E07 0.1562500E06 0.3125000E06 0.6250000E06 0.1250000E05 0.2500000E05 0.5000000E05 0.1000000E04 # Name of aerosol components, properties (no blank line a the end!) # NAME TYPE Molar mass Origin [0=bio,1=anth] BCAR 1 100. 1 DUST 1 100. 0 OCAR 1 100. 1 PPM 1 100. 1 SALT 1 58.5 0 AnA1D 2 168. 1 AnBmP 2 152. 1 BiA1D 2 170. 0 BiBmP 2 236. 0 H2SO4 3 96. 1 HNO3 3 62. 1 NH3 3 18. 1 WATER 4 18. 1

After the chemprep run: chemprep/inputdata.xxxx/ File Name

Parameters

Usage

List of all active species with corresponding horizontal transport scheme choices

Gas and aerosol transport

PRIMARY

Coefficients of decomposition of emitted species mass into CHIMERE size distribution

Aerosol emissions transport

AEROMIN.bin

Lookup table to compute aerosol solid/liquid Aerosol solid/liquid phase partitioning phase partitioning versus T, hum (based on ISORROPIA thermodynamic equilibrium model)

ACTIVE_SPECIES

AEROORG.bin

Aerosol parameters for the SOA scheme

SOA scheme

GOCART_builder

Boundary conditions: aggregation parameters GOCART ↔ CHIMERE species

Boundary conditions from global models

LMDZAERO_builder

Boundary conditions: aggregation parameters LMDZ/INCA ↔ CHIMERE species

Boundary conditions from global models


Domains: domains/

domains/domainlist.nml domain NX NY DX DY XMIN YMIN CONT1 65 33 0.5 0.5 10.0 40.5 CONT3 67 46 0.5 0.5 10.5 35.0 CONT4 70 44 0.5 0.5 10.5 36.0 CONT5 79 47 0.5 0.5 14.0 35.0 FRA10 101 111 0.15 0.1 5.0 41.0

makeCOORDdomains.sh

domains/HCOORD/ COORD_CONT5 domains/VCOORD/ VCOORD_8_995_500

in chimere.par nlayer_chimere top_chimere_pressure first_layer_pressure

define_geom.sh geom.awk


Meteorological preprocessor ➢ Can be used with any meteorological model. WRF shown below as an example ➢ Interfaces provided with CHIMERE: WRF, MM5, ECMWF (the latter one not supported)

# Meteorology Driver [$meteo] Meteo driver (WRF,MM5,ecm) : WRF, ecm [$metdom] Meteo driver domain : d01, IDF15, IDF5 # CHIMERE Subgrid Processes [$ideepconv] Deep convection activation (1/0) : 1 [$urbancorr] Urban correction [1] or not [0] : 0 [$resusp] Resuspension [1] or not [0] : 0 [$ilidar] Lidar profiles [1] or not [0] : 0 .......... [$meteo_DIR] Meteo driver output dir : $bigfiles/METEO [$meteo_file] Meteo file : # ${meteo_DIR}/wrfout_d01_$di
