CHIMERE Training Course

Transcription

1 CHIMERE Training Course General CHIMERE Structure Dmitry Khvorostyanov, Laurent Menut, Myrto Valari Laboratoire de Météorologie Dynamique IPSL Ecole Polytechnique, 928 Palaiseau Cedex Contact:

2 Chemistry Transport Model (CTM) Mass outflow Mass inflow

3 CHIMERE CORE

4 CHIMERE Core Routines

5 CHIMERE Pre processing scripts/make chemistry.sh the chemical pre processor chemprep/ chimere domain.sh to create domain specific data in domains/ : Domain coordinates: domains/hcoord/ COORD_CONT5 domains/vcoord/vcoord_8_995_500 Land use 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) make chemistry results in chemprep/inputdata.xxx/ directory chimere domains results in domains/cont5/ directory Specified in chimere.par: # CHIMERE Simulation Domain [$nested] Nested run? (yes/no) : no, yes, no [$dom] CHIMERE domain : CONT5, FRAGN5 [$nzdoms] Number of parallel zonal subdomains : 4 [$nmdoms] Number of parallel merid subdomains : 4 [$nlevels] Number of vertical layers : 8, 5 [$top_chimere_pressure] Top layer pressure : 500 [$first_layer_pressure] First layer pressure : 997

6 Chemical Preprocessor: chemprep/ Creates the directory chemprep/inputdata / where ' ' correspond to parameters specified in the parameter file chimere.par Calculates aerosol size distribution (src/chemprep distrib.f) Defines families (src/chemprep families.f) The resulting inputdata / 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 (/0) : [$nbins] Number of aerosol size sections :9 [$seasalt] Include sea salts? (/0) [$dust] Include dust? (/0) : [$carb] Allow carbon speciation? : [$soatyp] Secondary organic aerosol scheme (0..2) : 2 [$gtrc] Include gasesous tracers? [$ptrc] Include particulate tracers? [$nequil] ISORROPIA coupling () / tabulation (0) : 0 [$npeq] - equilibrium calc freq (phys steps, < 5) :

7 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 Species dependent parameters to calculate resistances for dry deposition Dry deposition initio/inidepo.f90 WETD_SPEC Species dependent 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 land use 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 gas phase chemistry Reaction rate constants Reaction rate calculation for gases (rates.f90) SEMIVOL Parameters for semi volatile species Absorption of semi volatile 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

8 AEROSOL # Number of sections D(2.5) D(0.0) E E E E E E E E E 04 # Name of aerosol components, properties (no blank line a the end!) # NAME TYPE Molar mass Origin [0=bio,=anth] BCAR 00. DUST OCAR 00. PPM 00. SALT AnAD AnBmP BiAD BiBmP H2SO HNO NH WATER 4 8. 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 Look up 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

9 Domains: domains/ domains/domainlist.nml domain NX NY DX DY XMIN CONT CONT CONT CONT FRA YMIN makecoorddomains.sh domains/vcoord/ VCOORD_8_995_500 in chimere.par nlayer_chimere top_chimere_pressure first_layer_pressure domains/hcoord/ COORD_CONT5 define_geom.sh geom.awk

10 Meteorological pre processor 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) [$metdom] Meteo driver domain # CHIMERE Subgrid Processes [$ideepconv] Deep convection activation (/0) [$urbancorr] Urban correction [] or not [0] [$resusp] Resuspension [] or not [0] [$ilidar] Lidar profiles [] or not [0]... [$meteo_dir] Meteo driver output dir [$meteo_file] Meteo file : WRF, ecm : d0, IDF5, IDF5 : : $bigfiles/meteo : # ${meteo_dir}/wrfout_d0_$di