EBAS example¶
import pyaerocom as pya
# if you are tired of warnings and output:
pya.change_verbosity('critical', log=pya.print_log)
from warnings import filterwarnings
filterwarnings('ignore')
pya.__version__
'0.10.0rc1'
Make sure you have version 0.10.0rc1 or later…
Register EBAS location (and reading interface)¶
This is needed in order to be able to read EBAS data.
pya.const.add_ungridded_obs(obs_id='EBAS',
data_dir=f'{pya.const.HOMEDIR}/shared-ns1000k/obsdata/EBAS',
reader=pya.io.ReadEbas,
check_read=True)
Read all available pm10 and pm2.5 mass concentration measurements from EBAS¶
reader = pya.io.ReadUngridded()
The following command will take a while! Have a coffee :)
data = reader.read('EBAS', ['concpm10', 'concpm25'],
start_date='2010-01-01', stop_date='2012-12-31')
Reading EBAS data
100%|██████████| 881/881 [00:22<00:00, 39.58it/s]
Set invalid data points to NaN¶
In EBAS, each measurement datapoint contains up to three flag values, some of which correspond to invalid measurements. By default, these values are not invalidated during reading of the data. Instead, the flags are evaluated and can be applied to the data object via:
data.set_flags_nan(inplace=True);
Compute coarse PM - colocating different ungridded measurements to infer new variables¶
Basically, we want to subtract pm2.5 from pm10 measurements for the EBAS data. However, you have to keep in mind that observations are sparse and you first have to identify common locations and times / sampling frequencies. Luckily, this can be done on the UngriddedData object:
coarse_pm = data.colocate_vardata(var1='concpm10',
data_id1 = 'EBAS',
var2='concpm25',
data_id2='EBAS',
match_stats_how='closest',
match_stats_tol_km=1,
merge_how='eval',
merge_eval_fun='EBAS;concpm10-EBAS;concpm25',
var_name_out='concpmgt25',
data_id_out='EBAS',
var_unit_out='ug cm-3',
resample_how='mean',
apply_time_resampling_constraints=True)
See what is in there¶
print(coarse_pm)
Pyaerocom UngriddedData
-----------------------
Contains networks: ['EBAS']
Contains variables: ['concpm10', 'concpm25', 'concpmgt25']
Contains instruments: ['filterpack_US2098', 'filterpack_US0065', 'lvs_02;teom_02;hvs_pm10_02', 'filterpack_US0025', ';b-attenuation_01;Eberline_FH_62_IR;b_att_01', 'filterpack_US0009', 'filterpack_US0076', 'filterpack_US3900', 'TEOM_Asp_12;sm200_Asp;opsis_sm200_12;leckel_12;lvs_se12', 'filterpack_US0029', 'tf_48;PM2501;lvs_48;PM25G01;PM1001;PM10G01;thermo_48', 'filterpack_US0072', 'filterpack_KR0001_BYIS1;filterpack_KR0001_BYIS1_2;filterpack_KR0001_BYISX_2', 'filterpack_US0043', 'filterpack_US0062', 'hvs_11;ES0011R_D-496;abs_11;ES0011R_D-492', 'filterpack_CA0107', 'filterpack_US4504', 'Grimm_190_BIR_dry;kfg_no01;dicho_d_0001;kfg_no02_pm25;kfg_no02_pm10;RP_TEOM;hvs_no01', 'filterpack_US0010', 'filterpack_US0019', 'filterpack_US0050', 'filterpack_US0035', 'filterpack_US0063', 'filterpack_US0111', 'filterpack_US0048', 'PM25_Bred;PM25_Bred_05;Grimm_Bre;PM_Bred', 'filterpack_US0041', 'filterpack_US0112;filterpack_US0112_2', 'filterpack_US0089', 'filterpack_US2002', 'filterpack_US0061', 'filterpack_US0087', 'filterpack_US4699', 'BAM1020_0644;Derenda-PNS16T;lvs_644;FAG_IR_EMS_nl0644', 'filterpack_US0055', 'hvs_10;ES0010R_A-120;abs_10;ES0010R_D-503', 'filterpack_US0040', 'ES0009R_1042;hvs_9;hvs_09;abs_9;P5150;ES0009R_P5150;ES0009R_140AB264500612;ES0009R_D-501', 'filterpack_US0011', 'filterpack_US0077', 'filterpack_US0042', 'filterpack_US4205', 'filterpack_US0079', 'filterpack_US0113', 'filterpack_US2096', 'filterpack_US0038', 'filterpack_US1005', 'digitel_4', 'filterpack_US0068', 'filterpack_US0036', 'tf_36;lvs_36', 'filterpack_US0398', 'filterpack_US0088', 'filterpack_US0028', 'filterpack_US0084', 'filterpack_US0083', 'filterpack_US0066', 'filterpack_US0070', 'filterpack_US1303', 'filterpack_US0059', 'BAM1020_0091;lvs_91_1;Derenda-PNS16T;FAG_IR_EMS_nl0091;FAG_IR_EMS_91', 'filterpack_US0020', 'filterpack_US0203', 'filterpack_US0047', 'Digitel_04;Digitel_4', 'HVPM25;HVPM10', 'filterpack_US0086', 'filterpack_US0022', 'filterpack_US0017', 'filterpack_CA0011', 'filterpack_US0051', 'filterpack_US0067', 'ES0014R_149;hvs_14;abs_14;ES0014R_D-513', 'filterpack_US1011', 'PM2', 'filterpack_US0115', 'filterpack_US0053', 'filterpack_US0116', 'filterpack_US0056', 'filterpack_US0018', 'filterpack_US0044', 'filterpack_US0117', 'filterpack_US0054', 'filterpack_US3201', 'filterpack_US2705', 'filterpack_US2905', 'filterpack_US0060', 'filterpack_US0049', 'filterpack_US0027', 'filterpack_US0073_GRSM1', 'filterpack_US0074', 'Scientific_SHARP_Monitor_Model_503', 'TEOM_FDMS_PM25_FR24;TEOM_FDMS_PM10_FR24', 'filterpack_US0030', 'filterpack_US0032', 'filterpack_US0013', 'TEOM_36;tf_36;PM2502;tf_36_dup;lvs_36;PM1002;thermo_36', 'filterpack_US0031', 'filterpack_US0071', 'filterpack_US0052', 'filterpack_US0016', 'filterpack_US0023', 'filterpack_US0012', 'filterpack_US0075;filterpack_US0075_2', 'kfg_no56_pm25;kfg_no56_pm10', 'DEKATI_PM10;DekatiPM10impactor', 'filterpack_US0015', 'hvs_02;DA80H_448;lvs_at02', 'filterpack_US0014', 'Leckel_pm25;lvs_08;LVS3D_D_01', 'filterpack_US0045', 'filterpack_US0046', 'f1p_04;JRC_PM25_01;f1p_4;LVS_RPP-2;LVS_KFG;leckel_31;lvs_it04', 'filterpack_US0080', 'filterpack_US0064;filterpack_US0064_2', 'filterpack_US0081', 'filterpack_US0567;filterpack_US0567_2', 'HVS;PM4;hvs_2;f1p_2;beta_02;PM1;AM;beta_ray_abs_2;PM3', 'filterpack_US0024', 'filterpack_US0069', 'filterpack_US0118', 'BAM1020_0009;lvs_09_1;Derenda-PNS16T;FAG_IR_EMS_nl0009;FH62_I_R;FAG_IR_EMS_9;lvs_nl09', 'Environnement_MP101_pm25;pm_FH_62_I-R;air_sampler_03;f1p_3_1;fp1_3;rad_cz03;Environnement_MP101_pm10;pm_FH_62_I-R_dup;lvs_cz03', 'kfg_no39_pm25;kfg_no39_pm10', 'TEOM_FDMS_PM25_FR18;MP101M_PM10_FR18;TEOM_FDMS_PM10_FR18', 'TEOM_FDMS_PM25_FR15;TEOM_FDMS_PM10_FR15;TEOM_15', 'EE01_ra2;hs1_09;f1p_09;hvs2', 'filterpack_US0033', 'filterpack_US0034', 'filterpack_US0119', 'filterpack_US0596', 'filterpack_US0026', 'Fai_DUAL.Channel_ECO_pm25;Fai_DUAL.Channel_ECO_pm10', 'teom_41;SM200;lvs_41', 'filterpack_US0057', 'filterpack_US0037', 'filterpack_US0120', 'filterpack_US0078', 'filterpack_US0058', 'filterpack_US0082', 'filterpack_US0121', 'filterpack_US9021', 'PM2.5;thermo_31;lvs_ie31', 'ES0006R_Digitel-151;DH-77-151;hvs_06;hvs_6;P7200TSPBLX;ES0006R_P7200TSPBLX;hvs_06_dup', 'filterpack_US9028', 'filterpack_US0122', 'filterpack_US1810', 'filterpack_US9006', 'filterpack_US0094', 'filterpack_US9014', 'digitel_44;Metrohm;hvs_44;DIONEX_3000;Andersen_PM10_Melpitz;thermo_44', 'filterpack_US0699_2', 'filterpack_US9011', 'filterpack_US9003', 'PARAMCD39_IDSTAZ63;PARAMCD25_IDSTAZ63', 'lvs_QBQ_1;lvs_t_01;LVS_KFG;abst_1_1;bgps_01;lvs_d_0001', 'GRIMM_1107;hvs_pm25_1178;hvs_17;hvs_17_dup;ES1778R_HVS_DIG_PM25_02;ES1778R_HVS_MCV_PM25_01;hvs_pm10_1178;thermo_1778', 'filterpack_US9012', 'filterpack_US9010', 'MP101M_PM25_FR10;TEOM_FDMS_PM25_FR10;MP101M_PM10_FR10;TEOM_FDMS_PM10_FR10', 'filterpack_US0091', 'filterpack_US9020', 'filterpack_US4999', 'filterpack_US9002', 'filterpack_US9015', 'UBA_Ng_HVS_0002;Digitel_07;UBA_Ng_HVS_0001;Digitel_7', 'filterpack_US9018', 'ES0008R_148;hvs_08;hvs_8;ES0008R_P7721X;abs_8', 'filterpack_US0123', 'filterpack_US9029', 'filterpack_US9027', 'filterpack_US9013', 'sm200_nor;TEOM_nor', 'abs_16;hvs_16;ES0016R_D-461;ES0016R_D-515;ES0016R_140AB264650612', 'filterpack_US1109', 'filterpack_US9009', 'filterpack_US9024;filterpack_US9024_2', 'filterpack_US9016', 'filterpack_US0306', 'filterpack_US0125', 'filterpack_US9017', 'filterpack_US9026', 'digitel_22;leckel_2', 'hvs_13;D-518;ES0013R_D-518;abs_13;D-516;ES0013R_D-516;ES0013R_140AB264340611', 'filterpack_US0126', 'filterpack_US0090', 'filterpack_US0095', 'TEOM_FDMS_PM25_FR13;TEOM_FDMS_PM10_FR13;TEOM_13', 'filterpack_US0092', 'filterpack_US0093;filterpack_US0093_2', 'filterpack_US0566', 'filterpack_US9022', 'filterpack_US0098', 'filterpack_US9008', 'filterpack_US9023_2', 'filterpack_US9025', 'PTLWS_pm25', 'filterpack_US9007;filterpack_US9007_2', 'filterpack_US9019', 'filterpack_US0096', 'filterpack_US0097', 'BAM_PM25_FR09;TEOM_FDMS_PM25_FR09;TEOM_FDMS_PM10_FR09;BAM_PM10_FR09;TEOM_09;thermo_09', 'digitel_52', 'filterpack_US0127', 'hvs_15;abs_15', 'filterpack_US9000', 'filterpack_US9001', 'filterpack_US0099', 'opsis_sm200_10_dup', 'PM25_Rao;PM_R_205_14;PM25_R_205_14;Grimm_Rao;PM_Rao', 'filterpack_US9057_2', 'filterpack_US9038', 'filterpack_US9039', 'MP101M_PM25_FR23;MP101M_PM10_FR23', 'filterpack_US9056', 'filterpack_US9068', 'filterpack_US9041', 'filterpack_US9043', 'ES0001R_Digitel-150;DH77-150;hvs_01;ES0001R_D-207;ES0001R_DH77-150;ES0001R_TISCH_595;ES0001R_595', 'filterpack_US9070', 'filterpack_US9046', 'filterpack_US9055', 'Digitel_03;UBA_Sc_HVS_0002;UBA_Sc_HVS_0001', 'UBA_Sm_HVS_0002;Digitel_08;UBA_Sm_HVS_0001;f1p_08', 'filterpack_US9060', 'filterpack_US9045', 'filterpack_US9047', 'filterpack_US4828', 'filterpack_US9065', 'filterpack_US9037', 'filterpack_US9059', 'filterpack_US9030', 'filterpack_US9033', 'filterpack_US9077', 'filterpack_US9040', 'filterpack_US9079', 'filterpack_US9052_2', 'filterpack_US3818', 'filterpack_US0128', 'filterpack_US9064', 'filterpack_US9066;filterpack_US9066_2', 'filterpack_US9035', 'filterpack_US0129', 'filterpack_US9058', 'filterpack_US3500', 'filterpack_US9074', 'filterpack_US9080', 'filterpack_US9036', 'filterpack_US9032;filterpack_US9032_2', 'filterpack_US9044', 'filterpack_US9031', 'filterpack_US9063', 'filterpack_US9034', 'TEOM_Vav;PM25_VAV_D_11;TEOM_Vav_11', 'filterpack_US9053', 'TEOM_Vin', 'filterpack_US5501', 'FMI_FH62IR_0846;lvs_17;FMI_FH62IR_0632;lvs_fi17', 'FMI_FH62IR_0846;FMI_FH62IR_0632', 'filterpack_US9061', 'filterpack_US9062', 'BAM1020_0010;Derenda-PNS16T;lvs_10;FAG_IR_EMS_nl0010;FAG_IR_EMS_10;FH62_I_R', 'abs_7;hvs_07;ES0007R_MCV_D-497;ES0007R_D-497;ES0007R_Tisch-0393;ES0007R_P7774X;hvs_7', 'Digitel_02;UBA_Wa_HVS_0002;UBA_Wa_HVS_0001;Digitel_2', 'filterpack_US9072', 'filterpack_US9051', 'filterpack_US9049', 'filterpack_US9071', 'filterpack_US9069', 'filterpack_US9078', 'filterpack_US9048', 'filterpack_US9073', 'filterpack_US4304;filterpack_US4304_2', 'filterpack_US0131', 'filterpack_US9081_YELL2;filterpack_US9081_YELL1', 'filterpack_US9042_YOSE1;filterpack_US9042_YOSEX_2;filterpack_US9042_YOSE1_2', 'ES0012R_C-280;hvs_12;abs_12;C-279;ES0012R_C-331;ES0012R_C-279;ES0012R_140AB264470611;ES0012R_D-492', 'HVPM25-Z;HVPM10-Z', 'filterpack_US9075', 'filterpack_US9076', 'hvs_48', 'opsis_sm200_16']
Total no. of meta-blocks: 283
coarse_pm.unique_station_names
['Acadia National Park-McFarland Hill (ME98)',
'Addison Pinnacle',
'Agia Marina Xyliatou ; Cyprus Atmospheric Observatory',
'Agua Tibia',
'Aliartos',
'Ambler',
'Arches NP',
'Arendtsville (PA00)',
'Aspvreten',
'Atlanta',
'Auchencorth Moss',
'Badlands NP',
'Baengnyeong Island',
'Baltimore',
'Bandelier NM',
'Barcarrota',
'Barrier Lake',
'Big Bend National Park-K-Bar (TX04)',
'Birkenes II',
'Birmingham',
'Bliss SP (TRPA)',
'Blue Mounds',
'Bondville',
'Bosque del Apache',
'Boulder Lake',
'Boundary Waters Canoe Area',
'Bredkälen',
'Breton',
'Breton Island',
'Bridger Wilderness',
'Bridgton (ME02)',
'Brigantine NWR',
'Brooklyn Lake',
'Bryce Canyon National Park-Repeater Hill (UT99)',
'Cabauw Wielsekade',
'Cabinet Mountains',
'Cabo de Creus',
'Cadiz',
'Campisabalos',
'Caney Creek',
'Canyonlands NP',
'Cape Cod',
'Cape Romain National Wildlife Refuge',
'Capitol Reef NP',
'Carlsbad Caverns',
"Casco Bay-Wolfe's Neck Farm (ME96)",
'Cedar Bluff',
'Chassahowitzka National Wildlife Refuge (FL05)',
'Chaumont',
'Cherokee Nation',
'Chicago',
'Chilbolton Observatory',
'Chiricahua (AZ98)',
'Cloud Peak',
'Cohutta',
'Columbia Gorge #1',
'Columbia River Gorge',
'Connecticut Hill',
'Crater Lake NP',
'Craters of the Moon National Monument (ID03)',
'Crescent Lake',
'De Zilk',
'Death Valley NP',
'Denali National Park-Mt. McKinley (AK03)',
'Detroit',
'Deuselbach',
'Diabla Gora',
'Dolly Sods Wilderness',
'Dome Lands Wilderness',
'Douglas',
'Egbert',
'El Dorado Springs',
'Ellis',
'Els Torms',
'Everglades National Park-Research Center (FL11)',
'Farkasfa',
'Flat Tops',
'Flathead',
'Forest County Potawatomi Community',
'Fort Peck (IMPROVE)',
'Fresno',
'Frostberg Reservoir (Big Piney Run)',
'Gates of the Arctic NP',
'Gates of the Mountains',
'Gila Cliff Dwellings National Monument (NM01)',
'Glacier National Park-Fire Weather Station (MT05)',
'Great Basin National Park-Lehman Caves (NV05)',
'Great Gulf Wilderness',
'Great River Bluffs',
'Great Sand Dunes NM',
'Great Smoky Mountains NP',
'Guadalupe Mountains NP',
'Gual Pahari',
'Guipry',
'Haleakala Crater',
'Haleakala NP',
'Hance Camp at Grand Canyon NP',
'Harwell',
'Hawaii Volcanoes NP',
'Hells Canyon',
'Hercules-Glades',
'Hillside',
'Hoover',
'Hopi Point #1',
'Houston',
'Hurdal',
'Hyytiälä',
"Ike's Backbone",
'Illmitz',
'Indian Gardens',
'Iskrba',
'Isle Royale NP',
'Isle Royale NP - code reserved pre 2014',
'Ispra',
'James River Face Wilderness',
'Jarbidge Wilderness',
'Jefferson NF',
'Joshua Tree National Park-Black Rock',
'K-puszta',
'Kaiser',
'Kalmiopsis',
'Kenai Peninsula Borough',
'Kollumerwaard',
'Kosetice (NOAK)',
'Kårvatn',
'La Coulonche',
'La Tardière',
'Lahemaa',
'Lake Sugema 1',
'Lake Sugema 2',
'Lake Tahoe Community College',
'Lassen Volcanic National Park-Manzanita Lake',
'Lava Beds NM',
'Lecce (ECO)',
'Lille Valby',
'Linville Gorge',
'Livonia',
'Londonderry',
'Lone Peak Wilderness',
'Lostwood',
'Lye Brook Wilderness',
'Lye Brook Wilderness (IMPROVE)',
'M.K. Goddard',
'Mace Head',
'Mahón',
'Makah Tribe',
'Makah Tribe Site #2',
'Mammoth Cave National Park-Houchin Meadow',
"Martha's Vineyard",
'Meadview',
'Medicine Lake',
'Melpitz',
'Mesa Verde National Park-Chapin Mesa (CO99)',
'Mingo',
'Mohawk Mt.',
'Monte Martano',
'Montelibretti',
'Montseny',
'Monture',
'Moosehorn NWR',
'Morvan',
'Mount Baldy',
'Mount Hood',
'Mount Rainier National Park-Tahoma Woods (WA99)',
'Mount Zirkel Wilderness',
'Nebraska NF',
'Neuglobsow',
'New York City',
'Niembro',
'Nogales',
'North Absaroka',
'North Cascades',
'Northern Cheyenne',
'Norunda Stenen',
'O Saviñao',
'Okefenokee National Wildlife Refuge (GA09)',
'Old Town',
'Olympic',
'Omaha',
'Organ Pipe Cactus National Monument',
'Owens Valley',
'Pack Monadnock Summit',
'Pasayten',
'Payerne',
'Penausende',
'Penobscot',
'Petersburg',
'Petrified Forest NP',
'Peyrusse Vieille',
'Phoenix',
'Phoenix Colocated Sampler',
'Pinnacles National Monument-Bear Valley (CA66)',
'Pittsburgh',
'Point Reyes National Seashore',
'Presque Isle',
'Proctor Maple R. F.',
'Puget Sound',
'Puntijarka',
'Quabbin Summit',
'Quaker City',
'Queen Valley',
'Redwood NP',
'Revin',
'Rigi',
'Ripple Creek',
'Risco Llamo',
'Rocky Mountain NP',
'Rocky Mountain NP HQ',
'Rubidoux',
'Rucava',
'Råö',
'Sac and Fox',
'Saguaro NM',
'Saguaro West',
'Saint-Nazaire-le-Désert',
'Salmon NF',
'Salt Creek',
'San Gabriel',
'San Gorgonio Wilderness',
'San Pablo de los Montes',
'San Pedro Parks',
'San Rafael',
'Sawtooth NF',
'Schauinsland',
'Schmücke',
'Seney',
'Sequoia NP',
'Shamrock Mine',
'Shenandoah National Park-Big Meadows (VA28)',
'Shining Rock Wilderness',
'Sierra Ancha',
'Sikes',
'Simeonof',
'Sipsy Wilderness',
'Snoqualmie Pass',
'South Lake Tahoe',
'Spokane Res.',
'St. Marks',
'Starkey Experimental Forest (OR18)',
'Stilwell (IMPROVE)',
'Sula Peak',
'Swanquarter',
'Sycamore Canyon',
'Sycamore Canyon 2',
'Tallgrass',
'Theodore Roosevelt National Park-Painted Canyon',
'Three Sisters Wilderness',
'Thunder Basin',
'Tonto NM',
'Trapper Creek',
'Trinity',
'Tuxedni',
'UL Bend',
'Upper Buffalo Wilderness',
'Vavihill',
'Viking Lake',
'Vindeln',
'Virgin Islands National Park-Lind Point (VI01)',
'Virolahti II',
'Virolahti III',
'Voyageurs NP #1',
'Voyageurs NP #2',
'Vredepeel',
'Víznar',
'Waldhof',
'Walker River Paiute Tribe',
'Washington D.C.',
'Weminuche Wilderness',
'Wheeler Peak',
'White Mountain',
'White Pass',
'White River NF',
'Wichita Mountains',
'Wind Cave National Park-Elk Mountain',
'Wrightwood',
'Yellowstone NP 1',
'Yosemite NP',
'Zarra',
'Zielonka',
'Zion',
'Zion Canyon',
'Zoebelboden',
'Zoseni']
sitedata = coarse_pm.to_station_data('Petersburg', vars_to_convert='concpmgt25')
sitedata
StationData([('dtime',
array(['2004-07-05T00:00:00.000000000', '2004-07-12T00:00:00.000000000',
'2004-07-19T00:00:00.000000000', '2004-07-26T00:00:00.000000000',
'2004-08-02T00:00:00.000000000', '2004-08-09T00:00:00.000000000',
'2004-08-16T00:00:00.000000000', '2004-08-23T00:00:00.000000000',
'2004-08-30T00:00:00.000000000', '2004-09-06T00:00:00.000000000',
'2004-09-13T00:00:00.000000000', '2004-09-20T00:00:00.000000000',
'2004-09-27T00:00:00.000000000', '2004-10-04T00:00:00.000000000',
'2004-10-11T00:00:00.000000000', '2004-10-18T00:00:00.000000000',
'2004-10-25T00:00:00.000000000', '2004-11-01T00:00:00.000000000',
'2004-11-08T00:00:00.000000000', '2004-11-15T00:00:00.000000000',
'2004-11-22T00:00:00.000000000', '2004-11-29T00:00:00.000000000',
'2004-12-06T00:00:00.000000000', '2004-12-13T00:00:00.000000000',
'2004-12-20T00:00:00.000000000', '2005-01-10T00:00:00.000000000',
'2005-01-17T00:00:00.000000000', '2005-01-24T00:00:00.000000000',
'2005-01-31T00:00:00.000000000', '2005-02-14T00:00:00.000000000',
'2005-02-21T00:00:00.000000000', '2005-02-28T00:00:00.000000000',
'2005-03-07T00:00:00.000000000', '2005-03-14T00:00:00.000000000',
'2005-03-21T00:00:00.000000000', '2005-03-28T00:00:00.000000000',
'2005-04-04T00:00:00.000000000', '2005-04-11T00:00:00.000000000',
'2005-04-18T00:00:00.000000000', '2005-04-25T00:00:00.000000000',
'2005-05-02T00:00:00.000000000', '2005-05-09T00:00:00.000000000',
'2005-05-16T00:00:00.000000000', '2005-05-23T00:00:00.000000000',
'2005-05-30T00:00:00.000000000', '2005-06-06T00:00:00.000000000',
'2005-06-13T00:00:00.000000000', '2005-06-20T00:00:00.000000000',
'2005-06-27T00:00:00.000000000', '2005-07-04T00:00:00.000000000',
'2005-07-11T00:00:00.000000000', '2005-07-18T00:00:00.000000000',
'2005-07-25T00:00:00.000000000', '2005-08-01T00:00:00.000000000',
'2005-08-08T00:00:00.000000000', '2005-08-15T00:00:00.000000000',
'2005-08-22T00:00:00.000000000', '2005-08-29T00:00:00.000000000',
'2005-09-05T00:00:00.000000000', '2005-09-12T00:00:00.000000000',
'2005-09-19T00:00:00.000000000', '2005-09-26T00:00:00.000000000',
'2005-10-03T00:00:00.000000000', '2005-10-10T00:00:00.000000000',
'2005-10-17T00:00:00.000000000', '2005-10-24T00:00:00.000000000',
'2005-10-31T00:00:00.000000000', '2005-11-07T00:00:00.000000000',
'2005-11-14T00:00:00.000000000', '2005-11-21T00:00:00.000000000',
'2005-11-28T00:00:00.000000000', '2005-12-05T00:00:00.000000000',
'2005-12-12T00:00:00.000000000', '2005-12-19T00:00:00.000000000',
'2005-12-26T00:00:00.000000000', '2006-01-02T00:00:00.000000000',
'2006-01-09T00:00:00.000000000', '2006-01-16T00:00:00.000000000',
'2006-01-23T00:00:00.000000000', '2006-01-30T00:00:00.000000000',
'2006-02-06T00:00:00.000000000', '2006-02-13T00:00:00.000000000',
'2006-02-20T00:00:00.000000000', '2006-02-27T00:00:00.000000000',
'2006-03-06T00:00:00.000000000', '2006-03-13T00:00:00.000000000',
'2006-03-20T00:00:00.000000000', '2006-03-27T00:00:00.000000000',
'2006-04-03T00:00:00.000000000', '2006-04-10T00:00:00.000000000',
'2006-04-17T00:00:00.000000000', '2006-04-24T00:00:00.000000000',
'2006-05-01T00:00:00.000000000', '2006-05-08T00:00:00.000000000',
'2006-05-15T00:00:00.000000000', '2006-05-22T00:00:00.000000000',
'2006-05-29T00:00:00.000000000', '2006-06-05T00:00:00.000000000',
'2006-06-12T00:00:00.000000000', '2006-06-19T00:00:00.000000000',
'2006-06-26T00:00:00.000000000', '2006-07-03T00:00:00.000000000',
'2006-07-10T00:00:00.000000000', '2006-07-17T00:00:00.000000000',
'2006-07-24T00:00:00.000000000', '2006-07-31T00:00:00.000000000',
'2006-08-07T00:00:00.000000000', '2006-08-14T00:00:00.000000000',
'2006-08-21T00:00:00.000000000', '2006-08-28T00:00:00.000000000',
'2006-09-04T00:00:00.000000000', '2006-09-11T00:00:00.000000000',
'2006-09-18T00:00:00.000000000', '2006-09-25T00:00:00.000000000',
'2006-10-02T00:00:00.000000000', '2006-10-09T00:00:00.000000000',
'2006-10-16T00:00:00.000000000', '2006-10-23T00:00:00.000000000',
'2006-10-30T00:00:00.000000000', '2006-11-06T00:00:00.000000000',
'2006-11-13T00:00:00.000000000', '2006-11-20T00:00:00.000000000',
'2006-11-27T00:00:00.000000000', '2006-12-04T00:00:00.000000000',
'2006-12-11T00:00:00.000000000', '2006-12-18T00:00:00.000000000',
'2006-12-25T00:00:00.000000000', '2007-01-08T00:00:00.000000000',
'2007-01-15T00:00:00.000000000', '2007-01-22T00:00:00.000000000',
'2007-01-29T00:00:00.000000000', '2007-02-05T00:00:00.000000000',
'2007-02-12T00:00:00.000000000', '2007-02-19T00:00:00.000000000',
'2007-02-26T00:00:00.000000000', '2007-03-05T00:00:00.000000000',
'2007-03-12T00:00:00.000000000', '2007-03-19T00:00:00.000000000',
'2007-03-26T00:00:00.000000000', '2007-04-02T00:00:00.000000000',
'2007-04-09T00:00:00.000000000', '2007-04-16T00:00:00.000000000',
'2007-04-23T00:00:00.000000000', '2007-04-30T00:00:00.000000000',
'2007-05-07T00:00:00.000000000', '2007-05-14T00:00:00.000000000',
'2007-05-21T00:00:00.000000000', '2007-05-28T00:00:00.000000000',
'2007-06-04T00:00:00.000000000', '2007-06-11T00:00:00.000000000',
'2007-06-18T00:00:00.000000000', '2007-06-25T00:00:00.000000000',
'2007-07-02T00:00:00.000000000', '2007-07-09T00:00:00.000000000',
'2007-07-16T00:00:00.000000000', '2007-07-23T00:00:00.000000000',
'2007-07-30T00:00:00.000000000', '2007-08-06T00:00:00.000000000',
'2007-08-13T00:00:00.000000000', '2007-08-20T00:00:00.000000000',
'2007-08-27T00:00:00.000000000', '2007-09-03T00:00:00.000000000',
'2007-09-10T00:00:00.000000000', '2007-09-17T00:00:00.000000000',
'2007-09-24T00:00:00.000000000', '2007-10-01T00:00:00.000000000',
'2007-10-08T00:00:00.000000000', '2007-10-15T00:00:00.000000000',
'2007-10-22T00:00:00.000000000', '2007-10-29T00:00:00.000000000',
'2007-11-05T00:00:00.000000000', '2007-11-12T00:00:00.000000000',
'2007-11-19T00:00:00.000000000', '2007-11-26T00:00:00.000000000',
'2007-12-03T00:00:00.000000000', '2007-12-10T00:00:00.000000000',
'2007-12-17T00:00:00.000000000', '2007-12-24T00:00:00.000000000',
'2007-12-31T00:00:00.000000000', '2008-01-07T00:00:00.000000000',
'2008-01-14T00:00:00.000000000', '2008-01-21T00:00:00.000000000',
'2008-01-28T00:00:00.000000000', '2008-02-11T00:00:00.000000000',
'2008-02-18T00:00:00.000000000', '2008-03-03T00:00:00.000000000',
'2008-03-10T00:00:00.000000000', '2008-03-17T00:00:00.000000000',
'2008-03-24T00:00:00.000000000', '2008-03-31T00:00:00.000000000',
'2008-04-07T00:00:00.000000000', '2008-04-14T00:00:00.000000000',
'2008-04-21T00:00:00.000000000', '2008-04-28T00:00:00.000000000',
'2008-05-05T00:00:00.000000000', '2008-05-12T00:00:00.000000000',
'2008-05-19T00:00:00.000000000', '2008-05-26T00:00:00.000000000',
'2008-06-09T00:00:00.000000000', '2008-06-16T00:00:00.000000000',
'2008-06-23T00:00:00.000000000', '2008-06-30T00:00:00.000000000',
'2008-07-07T00:00:00.000000000', '2008-07-14T00:00:00.000000000',
'2008-07-21T00:00:00.000000000', '2008-07-28T00:00:00.000000000',
'2008-08-04T00:00:00.000000000', '2008-08-11T00:00:00.000000000',
'2008-08-18T00:00:00.000000000', '2008-08-25T00:00:00.000000000',
'2008-09-01T00:00:00.000000000', '2008-09-08T00:00:00.000000000',
'2008-09-15T00:00:00.000000000', '2008-09-22T00:00:00.000000000',
'2008-09-29T00:00:00.000000000', '2008-10-06T00:00:00.000000000',
'2008-10-13T00:00:00.000000000', '2008-10-20T00:00:00.000000000',
'2008-10-27T00:00:00.000000000', '2008-11-03T00:00:00.000000000',
'2008-11-10T00:00:00.000000000', '2008-11-17T00:00:00.000000000',
'2008-11-24T00:00:00.000000000', '2008-12-01T00:00:00.000000000',
'2008-12-08T00:00:00.000000000', '2008-12-22T00:00:00.000000000',
'2008-12-29T00:00:00.000000000', '2009-01-12T00:00:00.000000000',
'2009-01-19T00:00:00.000000000', '2009-01-26T00:00:00.000000000',
'2009-02-02T00:00:00.000000000', '2009-02-09T00:00:00.000000000',
'2009-02-16T00:00:00.000000000', '2009-02-23T00:00:00.000000000',
'2009-03-02T00:00:00.000000000', '2009-03-09T00:00:00.000000000',
'2009-03-16T00:00:00.000000000', '2009-03-23T00:00:00.000000000',
'2009-03-30T00:00:00.000000000', '2009-04-06T00:00:00.000000000',
'2009-04-13T00:00:00.000000000', '2009-04-20T00:00:00.000000000',
'2009-05-11T00:00:00.000000000', '2009-05-18T00:00:00.000000000',
'2009-05-25T00:00:00.000000000', '2009-06-01T00:00:00.000000000',
'2009-06-08T00:00:00.000000000', '2009-06-15T00:00:00.000000000',
'2009-06-22T00:00:00.000000000', '2009-06-29T00:00:00.000000000',
'2009-07-06T00:00:00.000000000', '2009-07-13T00:00:00.000000000',
'2009-07-20T00:00:00.000000000', '2009-07-27T00:00:00.000000000',
'2009-08-17T00:00:00.000000000', '2009-08-31T00:00:00.000000000',
'2009-09-07T00:00:00.000000000', '2009-09-14T00:00:00.000000000',
'2009-09-21T00:00:00.000000000', '2009-09-28T00:00:00.000000000'],
dtype='datetime64[ns]')),
('var_info',
BrowseDict([('concpmgt25',
OrderedDict([('ts_type', 'weekly'),
('units', 'ug cm-3'),
('merge_how', 'eval'),
('merge_eval_fun',
'EBAS;concpm10-EBAS;concpm25'),
('overlap', False)]))])),
('station_coords',
{'latitude': 56.61076,
'longitude': -132.81234,
'altitude': 0.0}),
('data_err', BrowseDict()),
('overlap', BrowseDict()),
('numobs', BrowseDict()),
('data_flagged', BrowseDict()),
('filename',
'US0090R.20060102090000.20191216000000.filterpack..pm25.52w.1w.US03L_filterpack_US0090.US03L_IMPROVE_Davis.lev2.nas;US0090R.20050101090000.20191216000000.filterpack..pm25.52w.1w.US03L_filterpack_US0090.US03L_IMPROVE_Davis.lev2.nas;US0090R.20070103090000.20191216000000.filterpack...1y.3d.US03L_filterpack_US0090.US03L_IMPROVE_Davis.lev2.nas;US0090R.20080101090000.20191216000000.filterpack...52w.3d.US03L_filterpack_US0090.US03L_IMPROVE_Davis.lev2.nas;US0090R.20090101090000.20191216000000.filterpack...9mo.3d.US03L_filterpack_US0090.US03L_IMPROVE_Davis.lev2.nas;US0090R.20040702090000.20191216000000.filterpack..pm25.6mo.1w.US03L_filterpack_US0090.US03L_IMPROVE_Davis.lev2.nas'),
('station_id', 'US0090R'),
('station_name', 'Petersburg'),
('instrument_name', 'filterpack_US0090'),
('PI', 'MISSING INFORMATION'),
('country', None),
('country_code', None),
('ts_type', 'weekly'),
('latitude', 56.61076),
('longitude', -132.81234),
('altitude', 0.0),
('data_id', 'EBAS'),
('dataset_name', None),
('data_product', None),
('data_version', None),
('data_level', 2),
('revision_date', numpy.datetime64('2019-12-16T00:00:00')),
('website', None),
('ts_type_src', 'weekly'),
('stat_merge_pref_attr', None),
('data_revision', None),
('concpmgt25', 2004-07-05 1.249500
2004-07-12 1.701500
2004-07-19 2.634500
2004-07-26 1.850000
2004-08-02 1.723000
...
2009-08-31 1.526500
2009-09-07 1.987333
2009-09-14 1.102500
2009-09-21 1.452000
2009-09-28 0.490000
Length: 260, dtype: float64)])
sitedata['concpmgt25'].plot()
<matplotlib.axes._subplots.AxesSubplot at 0x7fe4d907c630>
Plot all coordinates¶
ax = coarse_pm.plot_station_coordinates(markersize=10)
coarse_pm.plot_station_coordinates(markersize=5, var_name='concpmgt25', start=2010, stop=2012, color='g', ax=ax);
f'No. of sites in data object:{len(coarse_pm.unique_station_names)}'
'No. of sites in data object:283'
model_id = 'ECMWF-IFS-CY45R1-CAMS-CTRL-met2010_AP3-CTRL'
home = pya.const.HOMEDIR
home
'/home/065403d4-2d5a02-2d46ec-2da0dd-2db878e4b7c040/'
pya.const.add_data_search_dir(f'{home}shared-ns1000k/aerocom_modeldata')b
pya.browse_database(model_id)
Pyaerocom ReadGridded
---------------------
Data ID: ECMWF-IFS-CY45R1-CAMS-CTRL-met2010_AP3-CTRL
Data directory: /home/065403d4-2d5a02-2d46ec-2da0dd-2db878e4b7c040/shared-ns1000k/aerocom_modeldata/ECMWF-IFS-CY45R1-CAMS-CTRL-met2010_AP3-CTRL/renamed
Available experiments: ['AP3']
Available years: [1850, 2010]
Available frequencies ['3hourly' 'daily' 'monthly']
Available variables: ['abs550aer', 'clh', 'cll', 'clm', 'clt', 'depbc', 'depdu', 'depoa', 'depso2', 'depso4', 'depss', 'drybc', 'drydu', 'dryoa', 'dryso2', 'dryso4', 'dryss', 'emibc', 'emidu', 'emioa', 'emiso2', 'emiss', 'iwp', 'loadbc', 'loaddu', 'loadlt10du', 'loadlt10ss', 'loadlt1du', 'loadlt1ss', 'loadlt25du', 'loadlt25ss', 'loadoa', 'loadso2', 'loadso4', 'loadss', 'lwp', 'od440aer', 'od550aer', 'od550bc', 'od550du', 'od550lt10du', 'od550lt10ss', 'od550lt1du', 'od550lt1ss', 'od550lt25du', 'od550lt25ss', 'od550oa', 'od550so4', 'od550ss', 'od865aer', 'prodso4', 'ps', 'sconcbc', 'sconcdu', 'sconcoa', 'sconcpm1', 'sconcpm10', 'sconcpm25', 'sconcso2', 'sconcso4', 'sconcss', 'wetbc', 'wetdu', 'wetoa', 'wetso2', 'wetso4', 'wetss', 'ang4487aer']
['ECMWF-IFS-CY45R1-CAMS-CTRL-met2010_AP3-CTRL']
reader = pya.io.ReadGridded(model_id)
model_pm10 = reader.read_var('sconcpm10', ts_type='daily', start=2010)
#lowres = model_pm10.regrid(lat_res_deg=5, lon_res_deg=10)