AI concentration and wind field slurper and plotter

ai
(ver. 5/9/06 )

Reads CAMx concentration and wind-field files for both
coarse and fine grids, and fiddles around with them

5/ 9/06 14:55:30

Hours are limited to 1968; model grid to 97 columns by 90 rows; layers to 1; multinested cells to 3300000.

Cell of highest emission for each species and each category
Reasonableness-checking of the chemical emission rates
National monitoring data
Ohio monitors and design values
Attainment (or not) Demonstration
Monitor-to-Model Relative Reduction Factors
Monitor location needs analysis
Back-trajectory analysis
Human exposure modeling
List of data sources

Summary of concentration data here

Day-by-day fine grid specifications here. Day-by-day fine grid specifications here. Reading of a second set of concentration files completed, from 2009_23_summer12.conf (but the values listed below are from the first set, unless otherwise indicated)

572 ozone monitors exist nationwide.

Header lines (3) from O3_design_values_5yr :
three header lines expected
don't know
just why

Results of the recalculation of design values here.

Here are the Ohio monitors present in the long list but absent from the short one:

0 found.

Attainment (or not) Demonstration

demonstration interval: 6/12/02 hr 1 Wed to 8/31/02 hr 24 Sat

Design values, relative reduction factors, and prospective design values for a bunch of monitors.

Old case: K23_summer12_w.conf
New case: 2009_23_summer12.conf

Only days with baseline modeled 8-hr concentration greater than 0.0700
participate in the calculation of rrf, unless otherwise noted.

					old	new	cell	lowest
station	state	county	rrf	days	design val.	design val.	size	accepted
390030002	Ohio	Allen	0.89007	10	88.0	78.33	12.0	0.0712
390071001	Ohio	Ashtabula	0.86741	23	95.6	82.92	12.0
390170004	Ohio	Butler	0.91512	15	89.6	81.99	12.0
390171004	Ohio	Butler	0.91604	14	87.6	80.25	12.0
390250022	Ohio	Clermont	0.91762	12	89.3	81.94	12.0
390610006	Ohio	Hamilton	0.94835	10	90.3	85.64	12.0
390610010	Ohio	Hamilton	0.90944	16	85.6	77.85	12.0
390610040	Ohio	Hamilton	0.96205	10	86.0	82.74	12.0	0.0813
391650006	Ohio	Warren	0.90794	13	90.6	82.26	12.0
390230001	Ohio	Clark	0.89195	10	88.6	79.03	12.0	0.0844
390230003	Ohio	Clark	0.88843	10	85.3	75.78	12.0	0.0821
390570006	Ohio	Greene	0.89451	10	87.6	78.36	12.0	0.0837
391090005	Ohio	Miami	0.88545	10	87.0	77.03	12.0	0.0791
391130019	Ohio	Montgomery	0.91479	10	86.5#	79.13	12.0	0.0804
391351001	Ohio	Preble	0.88344	10	80.0	70.67	12.0	0.0771
390271002	Ohio	Clinton	0.88080	10	94.3	83.06	12.0	0.0838
390350034	Ohio	Cuyahoga	0.89271	15	78.3	69.90	12.0
390350064	Ohio	Cuyahoga	0.87660	10	81.6	71.53	12.0	0.0839
390355002	Ohio	Cuyahoga	0.92574	10	88.0	81.47	12.0	0.0841
390550004	Ohio	Geauga	0.88391	12	99.0	87.51	12.0
390850003	Ohio	Lake	0.91543	10	92.6	84.77	12.0	0.0849
390853002	Ohio	Lake	0.91067	17	83.0	75.59	12.0
391331001	Ohio	Portage	0.88801	10	91.0	80.81	12.0	0.0809
391030003	Ohio	Medina	0.85803	10	87.0	74.65	12.0	0.0768
391530020	Ohio	Summit	0.90250	10	93.3	84.20	12.0	0.0762
390410002	Ohio	Delaware	0.86636	10	89.0	77.11	12.0	0.0759
390490028	Ohio	Franklin	0.88457	10	86.0#	76.07	12.0	0.0821
390490029	Ohio	Franklin	0.90067	10	93.0#	83.76	12.0	0.0792
390490037	Ohio	Franklin	0.90576	10	86.5#	78.35	12.0	0.0768
390490081	Ohio	Franklin	0.91969	10	84.0	77.25	12.0	0.0750
390830002	Ohio	Knox	0.85922	10	87.0	74.75	12.0	0.0807
390890005	Ohio	Licking	0.86323	10	88.0	75.96	12.0	0.0798
390970007	Ohio	Madison	0.88291	10	88.6	78.23	12.0	0.0797
390810016	Ohio	Jefferson	0.87535	10	86.0#	75.28	12.0	0.0790
390870006	Ohio	Lawrence	0.86608	10	82.0	71.02	12.0
390870011	Ohio	Lawrence	0.81399	10	76.6	62.35	12.0	0.0845
390930017	Ohio	Lorain	0.88094	10	87.0	76.64	12.0
390950024	Ohio	Lucas	0.94150	10	83.3	78.43	12.0	0.0704
390950027	Ohio	Lucas	0.88482	10	82.0	72.56	12.0	0.0737
390950034	Ohio	Lucas	0.93780	8	90.0	84.40	12.0	0.0728
390950081	Ohio	Lucas	0.94150	10	88.6	83.42	12.0	0.0704
391730003	Ohio	Wood	0.88135	10	87.6	77.21	12.0	0.0739
390990013	Ohio	Mahoning	0.86961	10	87.0	75.66	12.0	0.0845
391550011	Ohio	Trumbull	0.88099	10	92.5#	81.49	12.0	0.0808
391550009	Ohio	Trumbull	0.84272	10	88.0	74.16	12.0	0.0838
391510016	Ohio	Stark	0.85214	10	88.0	74.99	12.0	0.0774
391510021	Ohio	Stark	0.83986	10	86.3	72.48	12.0	0.0758
391511009	Ohio	Stark	0.89849	10	88.3	79.34	12.0	0.0767
391514005	Ohio	Stark	0.86482	10	88.3	76.36	12.0	0.0804
391670004	Ohio	Washington	0.79968	12	85.6	68.45	12.0
170310042	Illinois	Cook	1.01476	6	88*	89.30	12.0	0.0702
170317002	Illinois	Cook	1.02155	10	86*	87.85	12.0	0.0717
170971007	Illinois	Lake	0.95768	10	87*	83.32	12.0	0.0795
260990009	Michigan	Macomb	0.92552	10	97*	89.78	12.0	0.0755
260991003	Michigan	Macomb	0.95354	9	95*	90.59	12.0	0.0706
261630019	Michigan	Wayne	0.95987	3	91*	87.35	12.0	too few
420030010	Pennsylvania	Allegheny	0.88828	10	93*	82.61	12.0	0.0830
420070002	Pennsylvania	Beaver	0.87658	10	94*	82.40	12.0	0.0798
420490003	Pennsylvania	Erie	0.90192	12	92*	82.98	12.0

All design values are based on five years of data,
except those marked with "*" which are based on three, or"#" which are based on four.

The threshold is 85.0 ppb unless the number of qualifying days is less than 10,
in which case the threshold is adjusted downward as needed to qualify 10 days,
to a lower limit of 70.0 ppb, after which the acceptable number of days is reduced to as low as 5.

4516 coarse grid rrf's calculated; 4214 cells are uncalculable.
The highest rrf is 1.0446 at ( 71, 83); the lowest is 0.7810 at ( 57,116)

rrf calculations attempted for 1 nests totalling 16641 cells, over 81 days; 13382 cells calculated, 3259 uncalculable.
minimum fine-grid rrf = 0.713 in nest 1 at ( 3, 2) in coarse cell ( 40, 36); maximum fine-grid rrf = 1.071 in nest 1 at ( 2, 2) in coarse cell ( 61, 35)

Maximum red = 1.071
Deepest blue = 0.713

Attainment (or not) Demonstration

demonstration interval: 6/12/02 hr 1 Wed to 8/31/02 hr 24 Sat

Design values, relative reduction factors, and prospective design values for a bunch of monitors.

Old case: K23_summer12_w.conf
New case: 2009_23_summer12.conf

Only days with baseline modeled 8-hr concentration greater than 0.0700
participate in the calculation of rrf, unless otherwise noted.

A Timin collar will be employed. < ">390030002 ">390071001 ">390170004 ">390171004 ">390250022 ">390610006 ">390610010 ">390610040 ">391650006 ">390230001 ">390230003 ">390570006 ">391090005 ">391130019 ">391351001 ">390271002 ">390350034 ">390350064 ">390355002 ">390550004 ">390850003 ">390853002 ">391331001 ">391030003 ">391530020 ">390410002 ">390490028 ">390490029 ">390490037 ">390490081 ">390830002 ">390890005 ">390970007 ">390810016 ">390870006 ">390870011 ">390930017 ">390950024 ">390950027 ">390950034 ">390950081 ">391730003 ">390990013 ">391550011 ">391550009 ">391510016 ">391510021 ">391511009 ">391514005 ">391670004 ">170310042 ">170317002 ">170971007 ">260990009 ">260991003 ">261630019 ">420030010 ">420070002 ">420490003 those marked with "*" which are based on three, or"#" which are based on four. days is less than 10,
in which case the threshold is adjusted downward as needed to qualify 10 days,
to a lower limit of 70.0 ppb, after which the acceptable number of days is reduced to as low as 5. .
">( 71, 83); the lowest is 0.7810 at ( 57,116)
rrf calculations attempted for 1 nests totalling 16641 cells, over 81 days; 14337 cells calculated, 2304 uncalculable.
minimum fine-grid rrf = 0.680 in nest 1 at ( 2, 1) in coarse cell ( 40, 36); maximum fine-grid rrf = 1.018 in nest 1 at ( 3, 3) in coarse cell ( 61, 32)

Maximum red = 1.045
Deepest blue = 0.680
4516 coarse grid rrf to rrf ratios's calculated; 4214 cells are uncalculable.
The highest ratio is 1.0000 at ( 12, 68); the lowest is 1.0000 at ( 12, 68)
ratio of rrf's computed with ttt.rrf in the denominator

Maximum red = 1.105
Deepest blue = 0.800
Analysis performed 5/ 9/06 14:55:30
Culpability-related analysis of monitors of particular interest not performed because H2O2 and HNO3 concentrations not found in concentration files (needed to identify VOC and NOx-limited regimes) (ih2o2= -1; ihno3= -1
Monitor vs. model comparisons at 572 locations

old new cell lowest
station state county rrf days design val. design val. size accepted
Ohio Allen 0.88568 10 88.0 77.94 12.0 0.0764
Ohio Ashtabula 0.86973 29 95.6 83.15 12.0
Ohio Butler 0.92226 21 89.6 82.63 12.0
Ohio Butler 0.91318 21 87.6 79.99 12.0
Ohio Clermont 0.92105 25 89.3 82.25 12.0
Ohio Hamilton 0.94028 25 90.3 84.91 12.0
Ohio Hamilton 0.91929 23 85.6 78.69 12.0
Ohio Hamilton 0.93820 26 86.0 80.69 12.0
Ohio Warren 0.91077 22 90.6 82.52 12.0
Ohio Clark 0.88862 13 88.6 78.73 12.0
Ohio Clark 0.90088 15 85.3 76.84 12.0
Ohio Greene 0.89618 15 87.6 78.51 12.0
Ohio Miami 0.88311 10 87.0 76.83 12.0
Ohio Montgomery 0.89904 16 86.5# 77.77 12.0
Ohio Preble 0.88043 10 80.0 70.43 12.0 0.0814
Ohio Clinton 0.88831 16 94.3 83.77 12.0
Ohio Cuyahoga 0.89117 27 78.3 69.78 12.0
Ohio Cuyahoga 0.88365 18 81.6 72.11 12.0
Ohio Cuyahoga 0.90195 26 88.0 79.37 12.0
Ohio Geauga 0.89994 16 99.0 89.09 12.0
Ohio Lake 0.89641 29 92.6 83.01 12.0
Ohio Lake 0.89464 29 83.0 74.25 12.0
Ohio Portage 0.87578 12 91.0 79.70 12.0
Ohio Medina 0.85779 10 87.0 74.63 12.0 0.0814
Ohio Summit 0.87340 10 93.3 81.49 12.0 0.0834
Ohio Delaware 0.87527 10 89.0 77.90 12.0 0.0794
Ohio Franklin 0.89211 10 86.0# 76.72 12.0 0.0850
Ohio Franklin 0.89105 10 93.0# 82.87 12.0
Ohio Franklin 0.89182 11 86.5# 77.14 12.0
Ohio Franklin 0.87950 10 84.0 73.88 12.0 0.0842
Ohio Knox 0.86581 11 87.0 75.33 12.0
Ohio Licking 0.87344 10 88.0 76.86 12.0 0.0833
Ohio Madison 0.87986 10 88.6 77.96 12.0 0.0837
Ohio Jefferson 0.83677 12 86.0# 71.96 12.0
Ohio Lawrence 0.87391 17 82.0 71.66 12.0
Ohio Lawrence 0.84600 20 76.6 64.80 12.0
Ohio Lorain 0.88145 23 87.0 76.69 12.0
Ohio Lucas 0.91372 10 83.3 76.11 12.0 0.0792
Ohio Lucas 0.90294 10 82.0 74.04 12.0 0.0763
Ohio Lucas 0.89582 15 90.0 80.62 12.0
Ohio Lucas 0.91372 10 88.6 80.96 12.0 0.0792
Ohio Wood 0.88322 10 87.6 77.37 12.0 0.0780
Ohio Mahoning 0.86882 12 87.0 75.59 12.0
Ohio Trumbull 0.86352 13 92.5# 79.88 12.0
Ohio Trumbull 0.84835 10 88.0 74.66 12.0 0.0845
Ohio Stark 0.86151 10 88.0 75.81 12.0 0.0822
Ohio Stark 0.83985 10 86.3 72.48 12.0 0.0786
Ohio Stark 0.86550 10 88.3 76.42 12.0 0.0827
Ohio Stark 0.86639 11 88.3 76.50 12.0
Ohio Washington 0.78983 26 85.6 67.61 12.0
Illinois Cook 0.98434 16 88* 86.62 12.0
Illinois Cook 0.98255 18 86* 84.50 12.0
Illinois Lake 0.94853 17 87* 82.52 12.0
Michigan Macomb 0.92494 10 97* 89.72 12.0 0.0828
Michigan Macomb 0.93331 10 95* 88.66 12.0 0.0752
Michigan Wayne 0.92423 10 91* 84.11 12.0 0.0782
Pennsylvania Allegheny 0.88464 18 93* 82.27 12.0
Pennsylvania Beaver 0.88022 10 94* 82.74 12.0
Pennsylvania Erie 0.88189 26 92* 81.13 12.0

>>>>

monitor model
matched mod.-mon.

station pairs mean err. RMS err. mon. mod.

Here are all the coarse-grid cells with at least one predicted 8-hr avg. greater than 0.1650

cell	highest 8-hr	begins	2nd highest 8-hr	begins
( 42, 76)	0.1753	7/ 4/02 hr 11 Thu	0.1749	7/ 4/02 hr 10 Thu
( 52, 79)	0.1660	6/24/02 hr 12 Mon	0.1638	6/24/02 hr 11 Mon
( 53, 80)	0.1851	7/29/02 hr 12 Mon	0.1825	7/ 3/02 hr 12 Wed
( 53, 81)	0.1852	7/29/02 hr 13 Mon	0.1834	7/29/02 hr 12 Mon
( 53, 82)	0.1669	7/30/02 hr 12 Tue	0.1665	7/29/02 hr 13 Mon
( 54, 82)	0.1718	7/29/02 hr 14 Mon	0.1702	7/29/02 hr 15 Mon
( 54, 83)	0.1656	7/29/02 hr 14 Mon	0.1642	7/29/02 hr 15 Mon
( 55, 83)	0.1721	7/29/02 hr 14 Mon	0.1713	7/29/02 hr 13 Mon
( 55, 84)	0.1685	7/29/02 hr 14 Mon	0.1666	7/29/02 hr 15 Mon
( 56, 84)	0.1668	7/29/02 hr 15 Mon	0.1662	7/29/02 hr 14 Mon

10 cells are affected.

Here are all the fine-grid cells with at least one predicted 1-hr avg. greater than 0.14500 (bypass this logic for now)

Here are all the cells from the coarse grid with at least one predicted 8-hr avg. conc. of O3 greater than 0.16000
that are a local peak, or that are within 420.km of central Ohio and have an 8-hr avg. greater than 0.11500
and are a local peak.

cell	highest 8-hr	begins	2nd highest 8-hr	begins
( 42, 76)	0.1753	7/ 4/02 hr 11 Thu	0.1749	7/ 4/02 hr 10 Thu
( 46, 75)	0.1633	7/ 3/02 hr 11 Wed	0.1562	7/ 3/02 hr 12 Wed
( 52, 79)	0.1660	6/24/02 hr 12 Mon	0.1638	6/24/02 hr 11 Mon
( 53, 81)	0.1852	7/29/02 hr 13 Mon	0.1834	7/29/02 hr 12 Mon
( 54, 82)	0.1718	7/29/02 hr 14 Mon	0.1702	7/29/02 hr 15 Mon
( 55, 83)	0.1721	7/29/02 hr 14 Mon	0.1713	7/29/02 hr 13 Mon
( 60, 84)	0.1612	7/ 3/02 hr 13 Wed	0.1610	7/ 3/02 hr 12 Wed

7 cells are affected.

The highest coarse gridwide average concentration is 0.04674 for hour 232( 6/21/02 hr 16 Fri)

The highest cell concentration is 0.21154 at ( 53, 80) for hour 1143( 7/29/02 hr 15 Mon)

The highest cell concentration within a 420.0 km radius of ( 1188.8, 83.1) (which corresponds to Columbus, Ohio) is 0.1261 at ( 53, 49) for hour 975( 7/22/02 hr 15 Mon)

back-trajectory analysis declined.

"Magno conatu nugas"

Nothing convinces persons of a weak understanding so effectually as what they do not comprehend.

"Political Maxims from Cardinal de Retz"

					monitor	model
	matched	mod.-mon.
station	pairs	mean err.	RMS err.	mon.	mod.

ai(ver. 5/9/06 ) Reads CAMx concentration and wind-field files for both coarse and fine grids, and fiddles around with them

ai

Reads CAMx concentration and wind-field files for both coarse and fine grids, and fiddles around with them

5/ 9/06 14:55:30

Contents

Attainment (or not) Demonstration

Attainment (or not) Demonstration

Monitor vs. model comparisons at 572 locations

Here are all the coarse-grid cells with at least one predicted 8-hr avg. greater than 0.1650

Here are all the fine-grid cells with at least one predicted 1-hr avg. greater than 0.14500 (bypass this logic for now)

Here are all the cells from the coarse grid with at least one predicted 8-hr avg. conc. of O3 greater than 0.16000 that are a local peak, or that are within 420.km of central Ohio and have an 8-hr avg. greater than 0.11500and are a local peak.

ai
(ver. 5/9/06 )

Reads CAMx concentration and wind-field files for both
coarse and fine grids, and fiddles around with them

Reads CAMx concentration and wind-field files for both
coarse and fine grids, and fiddles around with them

Here are all the cells from the coarse grid with at least one predicted 8-hr avg. conc. of O3 greater than 0.16000
that are a local peak, or that are within 420.km of central Ohio and have an 8-hr avg. greater than 0.11500
and are a local peak.