...
We tested with the CAMS reanalysis data for driving the WRF-Chem model. Please have a look at how to download the CAMS Reanalysis data via the ECMWF Web API so that you have a good understanding of the scripts below. Scripts below are using ECMWF's Web API.
Script for retrieving aerosols, parameter names are given in the table bellow:
Param | Short Name | Long Name |
---|---|---|
t | Temperature | |
1.210 | aermr01 | Sea Salt Aerosol (0.03 - 0.5 um) Mixing Ratio |
2.210 | aermr02 | Sea Salt Aerosol (0.5 - 5 um) Mixing Ratio |
3.210 | aermr03 | Sea Salt Aerosol (5 - 20 um) Mixing Ratio |
4.210 | aermr04 | Dust Aerosol (0.03 - 0.55 um) Mixing Ratio |
5.210 | aermr05 | Dust Aerosol (0.55 - 0.9 um) Mixing Ratio |
6.210 | aermr06 | Dust Aerosol (0.9 - 20 um) Mixing Ratio |
7.210 | aermr07 | Hydrophobic Organic Matter Aerosol Mixing Ratio |
8.210 | aermr08 | Hydrophilic Organic Matter Aerosol Mixing Ratio |
9.210 | aermr09 | Hydrophobic Black Carbon Aerosol Mixing Ratio |
10.210 | aermr10 | Hydrophilic Black Carbon Aerosol Mixing Ratio |
11.210 | aermr11 | Sulphate Aerosol Mixing Ratio |
Code Block | ||||||
---|---|---|---|---|---|---|
| ||||||
#!/usr/bin/env python
from ecmwfapi import ECMWFDataServer
server = ECMWFDataServer()
server.retrieve({
'dataset' : "cams_reanalysis | ||||||
Code Block | ||||||
| ||||||
#!/usr/bin/env python from ecmwfapi import ECMWFDataServer server = ECMWFDataServer() server.retrieve({ 'dataset' : "cams_reanalysis", 'class' : "mc", 'type' : "an", 'stream' : "oper", 'expver' : "eac4", 'levtype' : "ml", 'repres' : "gg", 'levellist' : "7/to/60", 'param' : "T/1.210/2.210/3.210/4.210/5.210/6.210/7.210/8.210/9.210/10.210/11.210", 'date' : "20100101/to/20100102", 'stepclass' : "0mc", 'time 'type' : "0/to/21/by/3an", 'grid' stream' : "0.75/0.75oper", 'target 'expver' : "AER_20100101_20100102_eac4.grib" }) | ||||||
Code Block | ||||||
| ||||||
#!/usr/bin/env python from ecmwfapi import ECMWFDataServer server = ECMWFDataServer() server.retrieve({ 'dataset' : "cams_reanalysis", 'class' : "mc", 'type' : "an", 'stream' : "oper", 'expver' : "eac4", 'levtype' : "ml, 'levtype' : "ml", 'repres' : "gg", 'levellist' : "7/to/60", 'param' : "T/1.210/2.210/3.210/4.210/5.210/6.210/7.210/8.210/9.210/10.210/11.210", 'date' : "20100101/to/20100102", 'repres' step' : "gg0", 'levellist' : "7/to/60", 'param 'time' : "203.210/123.210/27.217/121.210/13.217/6.217/124.210/122.2100/to/21/by/3", 'date' grid' : "20100101/to/201001020.75/0.75", 'steptarget' : "0", 'time' : "0/to/21/by/3", 'grid' : "0.75/0.75", 'target' : "GRG_20100101_20100102_eac4.grib" }) AER_20100101_20100102_eac4.grib" }) |
Script for retrieving global reactive gases, parameter names are given in the table bellow:
Param | Short Name | Long Name |
---|---|---|
203.210 | O3 | ozone |
123.210 | CO | carbonmonoxide |
27.217 | NO | nitrogen_monoxide |
121.210 | NO2 | nitrogen_dioxide |
13.217 | PAN | peroxyacetyl_nitrate |
6.217 | HNO3 | nitric_acid |
124.210 | CH2O | formaldehyde |
122.210 | SO2 | sulfur_dioxide |
Code Block | ||||||
---|---|---|---|---|---|---|
| ||||||
#!/usr/bin/env python
from ecmwfapi import ECMWFDataServer
server = ECMWFDataServer()
server.retrieve({
'dataset' : "cams_reanalysis",
'class' : "mc",
'type' : "an",
'stream' : "oper",
'expver' : "eac4",
'levtype' : "ml",
'repres' : "gg",
'levellist' : "7/to/60",
'param' : "203.210/123.210/27.217/121.210/13.217/6.217/124.210/122.210",
'date' : "20100101/to/20100102",
'step' : "0",
'time' : "0/to/21/by/3",
'grid' : "0.75/0.75",
'target' : "GRG_20100101_20100102_eac4.grib"
})
|
Script for retrieving volatile organic compounds, parameter names are given in the table bellow:
param | Short Name | Long Name |
---|---|---|
52.217 | CH3COCH3 | acetone |
45.217 | C2H6 | ethane |
9.217 | PAR | paraffins |
42.217 | CH3OH | methanol |
47.217 | C3H8 | propane |
46.217 | C2H5OH | ethanol |
10.217 | C2H4 | ethene |
12.217 | ALD2 | aldehydes |
11.217 | OLE | olefins |
16.217 | C5H8 | isoprene |
43.217 | HCOOH | formic_acid |
7.217 | CH3OOH | methylperoxide |
15.217 | ONIT | organic_nitrates |
Code Block | ||||||
---|---|---|---|---|---|---|
| ||||||
#!/usr/bin/env python from ecmwfapi import ECMWFDataServer server = ECMWFDataServer() server.retrieve({ 'dataset' : "cams_reanalysis", 'class' : "mc", 'type' : "an", 'stream' : "oper", 'expver' : "eac4", 'levtype' : "ml", 'repres' : "gg", 'levellist' : "7/to/60", 'param' : "52.217/45.217/9.217/42.217/47.217/46.217/10.217/12.217/11.217/16.217/43.217/7.217/15.217", 'date' : "20100101/to/20100102", 'step' : "0", 'time' : "0/to/21/by/3", 'grid' : "0.75/0.75", 'target' : "GRG_voc_20100101_20100102_eac4.grib" }) |
Script for retrieving h2o2:
Code Block | ||||||
---|---|---|---|---|---|---|
| ||||||
#!/usr/bin/env python from ecmwfapi import ECMWFDataServer server = ECMWFDataServer() server.retrieve({ 'dataset' : "cams_reanalysis", 'class' : "mc", 'type' : "an", 'stream' : "oper", 'expver' : "eac4", 'levtype' : "ml", 'repres' : "gg", 'levellist' : "7/to/60", 'param' : "3.217", 'date' : "20100101/to/20100102", 'step' : "0", 'time' : "0/to/21/by/3", 'grid' : "0.75/0.75", 'target' : "GRG_h2o2_20100101_20100102_eac4.grib" }) |
...