OpenSWMM Engine  6.0.0-alpha.4
Data-oriented, plugin-extensible SWMM Engine (6.0.0-alpha.4)
Loading...
Searching...
No Matches
xsect_tables.hpp
Go to the documentation of this file.
1// SPDX-License-Identifier: Apache-2.0
2//
3// Copyright 2026 Caleb Buahin
4//
5// Licensed under the Apache License, Version 2.0 (the "License");
6// you may not use this file except in compliance with the License.
7// You may obtain a copy of the License at
8//
9// http://www.apache.org/licenses/LICENSE-2.0
10//
11// Unless required by applicable law or agreed to in writing, software
12// distributed under the License is distributed on an "AS IS" BASIS,
13// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14// See the License for the specific language governing permissions and
15// limitations under the License.
16
39
40#ifndef OPENSWMM_XSECT_TABLES_HPP
41#define OPENSWMM_XSECT_TABLES_HPP
42
43namespace openswmm::xsect_tables {
44
45// ============================================================================
46// Amax — ratio of area at max flow to full area (from xsect.c)
47// ============================================================================
48
49static constexpr double Amax[26] = {
50 1.0, // 0 - DUMMY
51 0.9756, // 1 - CIRCULAR
52 0.9756, // 2 - FILLED_CIRCULAR
53 0.97, // 3 - RECT_CLOSED
54 1.0, // 4 - RECT_OPEN
55 1.0, // 5 - TRAPEZOIDAL
56 1.0, // 6 - TRIANGULAR
57 1.0, // 7 - PARABOLIC
58 1.0, // 8 - POWERFUNC
59 0.98, // 9 - RECT_TRIANG
60 0.98, // 10 - RECT_ROUND
61 0.96, // 11 - MOD_BASKET
62 0.96, // 12 - HORIZ_ELLIPSE
63 0.96, // 13 - VERT_ELLIPSE
64 0.92, // 14 - ARCH
65 0.96, // 15 - EGGSHAPED
66 0.96, // 16 - HORSESHOE
67 0.96, // 17 - GOTHIC
68 0.98, // 18 - CATENARY
69 0.98, // 19 - SEMIELLIPTICAL
70 0.96, // 20 - BASKETHANDLE
71 0.96, // 21 - SEMICIRCULAR
72 1.0, // 22 - IRREGULAR
73 0.96, // 23 - CUSTOM
74 0.9756, // 24 - FORCE_MAIN
75 1.0 // 25 - STREET_XSECT
76};
77
78// ============================================================================
79// CIRCULAR
80// ============================================================================
81
82static constexpr int N_A_Circ = 51;
83static constexpr double A_Circ[51] = {
84 0.0,.00471,.0134,.024446,.0374,.05208,.0680,.08505,.1033,.12236,
85 .1423,.16310,.1845,.20665,.2292,.25236,.2759,.29985,.3242,.34874,
86 .3736,.39878,.4237,.44907,.4745,.500,.5255,.55093,.5763,.60135,
87 .6264,.65126,.6758,.70015,.7241,.74764,.7708,.79335,.8154,.83690,
88 .8576,.87764,.8967,.91495,.9320,.94792,.9626,.97555,.9866,.99516,1.000
89};
90
91static constexpr int N_R_Circ = 51;
92static constexpr double R_Circ[51] = {
93 .0100,.0528,.1048,.1556,.2052,.254,.3016,.3484,.3944,.4388,
94 .4824,.5248,.5664,.6064,.6456,.6836,.7204,.7564,.7912,.8244,
95 .8568,.888,.9176,.9464,.9736,1.00,1.024,1.048,1.070,1.0912,1.110,1.1272,
96 1.144,1.1596,1.174,1.1848,1.194,1.2024,1.210,1.2148,
97 1.217,1.2172,1.215,1.2104,1.203,1.192,1.178,1.1584,1.132,1.094,1.000
98};
99
100static constexpr int N_Y_Circ = 51;
101static constexpr double Y_Circ[51] = {
102 0.0, 0.05236, 0.08369, 0.11025, 0.13423, 0.15643, 0.17755, 0.19772, 0.21704,
103 0.23581, 0.25412, 0.27194, 0.28948, 0.30653, 0.32349, 0.34017, 0.35666,
104 0.37298, 0.38915, 0.40521, 0.42117, 0.43704, 0.45284, 0.46858, 0.4843,
105 0.50000, 0.51572, 0.53146, 0.54723, 0.56305, 0.57892, 0.59487, 0.61093,
106 0.62710, 0.64342, 0.65991, 0.67659, 0.69350, 0.71068, 0.72816, 0.74602,
107 0.76424, 0.78297, 0.80235, 0.82240, 0.84353, 0.86563, 0.88970, 0.91444,
108 0.94749, 1.0
109};
110
111static constexpr int N_S_Circ = 51;
112static constexpr double S_Circ[51] = {
113 0.0, 0.00529, 0.01432, 0.02559, 0.03859, 0.05304, 0.06877, 0.08551, 0.10326,
114 0.12195, 0.14144, 0.16162, 0.18251, 0.2041, 0.22636, 0.24918, 0.27246,
115 0.29614, 0.32027, 0.34485, 0.36989, 0.39531, 0.42105, 0.44704, 0.47329,
116 0.4998, 0.52658, 0.55354, 0.58064, 0.60777, 0.63499, 0.66232, 0.68995,
117 0.7177, 0.74538, 0.77275, 0.79979, 0.82658, 0.8532, 0.87954, 0.90546,
118 0.93095, 0.95577, 0.97976, 1.00291, 1.02443, 1.04465, 1.06135, 1.08208,
119 1.07662, 1.0
120};
121
122static constexpr int N_W_Circ = 51;
123static constexpr double W_Circ[51] = {
124 0.0, .2800, .3919, .4750, .5426, .6000, .6499, .6940, .7332,
125 .7684, .8000, .8285, .8542, .8773, .8980, .9165, .9330,
126 .9474, .9600, .9708, .9798, .9871, .9928, .9968, .9992,
127 1.000, .9992, .9968, .9928, .9871, .9798, .9708, .9600,
128 .9474, .9330, .9165, .8980, .8773, .8542, .8285, .8000,
129 .7684, .7332, .6940, .6499, .6000, .5426, .4750, .3919,
130 .2800, .0
131};
132
133// ============================================================================
134// EGG SHAPE
135// ============================================================================
136
137static constexpr int N_A_Egg = 26;
138static constexpr double A_Egg[26] = {
139 .0000,.0116,.0316,.0550,.0850,.1200,.1555,.1900,.2250,.2750,
140 .3200,.3700,.4200,.4700,.5150,.5700,.6200,.6800,.7300,.7800,
141 .8350,.8850,.9250,.9550,.9800,1.000
142};
143
144static constexpr int N_R_Egg = 26;
145static constexpr double R_Egg[26] = {
146 .0100,.1300,.2440,.3020,.3860,.4650,.5360,.6110,.6760,.7350,
147 .7910,.8540,.9040,.9410,1.008,1.045,1.076,1.115,1.146,1.162,
148 1.186,1.193,1.186,1.162,1.107,1.000
149};
150
151static constexpr int N_Y_Egg = 51;
152static constexpr double Y_Egg[51] = {
153 0.0, 0.04912, 0.08101, 0.11128, 0.14161, 0.16622, 0.18811, 0.21356, 0.23742,
154 0.25742, 0.27742, 0.29741, 0.31742, 0.33742, 0.35747, 0.37364, 0.4,
155 0.41697, 0.43372, 0.45, 0.46374, 0.47747, 0.49209, 0.50989, 0.53015,
156 0.55, 0.56429, 0.57675, 0.58834, 0.6, 0.61441, 0.62967, 0.64582,
157 0.66368, 0.68209, 0.7, 0.71463, 0.72807, 0.74074, 0.75296, 0.765,
158 0.77784, 0.79212, 0.80945, 0.82936, 0.85, 0.86731, 0.88769, 0.914,
159 0.95, 1.0
160};
161
162static constexpr int N_S_Egg = 51;
163static constexpr double S_Egg[51] = {
164 0.0, 0.00295, 0.01331, 0.02629, 0.04, 0.05657, 0.075, 0.09432, 0.11473,
165 0.13657, 0.15894, 0.1803, 0.20036, 0.22, 0.23919, 0.25896, 0.28,
166 0.30504, 0.33082, 0.35551, 0.37692, 0.39809, 0.42, 0.44625, 0.47321,
167 0.5, 0.52255, 0.54481, 0.56785, 0.59466, 0.62485, 0.65518, 0.68181,
168 0.70415, 0.72585, 0.74819, 0.77482, 0.80515, 0.83534, 0.86193, 0.88465,
169 0.9069, 0.93, 0.95866, 0.98673, 1.01238, 1.03396, 1.05, 1.06517,
170 1.0538, 1.0
171};
172
173static constexpr int N_W_Egg = 26;
174static constexpr double W_Egg[26] = {
175 .0, .3250, .4270, .5080, .5820, .6420, .6960, .7460, .7910,
176 .8360, .8660, .8960, .9260, .9560, .9700, .9850, 1.000,
177 .9850, .9700, .9400, .8960, .8360, .7640, .6420, .3100,
178 .0
179};
180
181// ============================================================================
182// HORSESHOE
183// ============================================================================
184
185static constexpr int N_A_Horseshoe = 26;
186static constexpr double A_Horseshoe[26] = {
187 .0000,.0181,.0508,.0908,.1326,.1757,.2201,.2655,.3118,.3587,
188 .4064,.4542,.5023,.5506,.5987,.6462,.6931,.7387,.7829,.8253,
189 .8652,.9022,.9356,.9645,.9873,1.000
190};
191
192static constexpr int N_R_Horseshoe = 26;
193static constexpr double R_Horseshoe[26] = {
194 .0100,.1040,.2065,.3243,.4322,.5284,.6147,.6927,.7636,.8268,
195 .8873,.9417,.9905,1.036,1.077,1.113,1.143,1.169,1.189,1.202,
196 1.208,1.206,1.195,1.170,1.126,1.000
197};
198
199static constexpr int N_Y_Horseshoe = 51;
200static constexpr double Y_Horseshoe[51] = {
201 0.0, 0.04146, 0.07033, 0.09098, 0.10962, 0.12921, 0.14813, 0.16701, 0.18565,
202 0.20401, 0.22211, 0.23998, 0.25769, 0.27524, 0.29265, 0.3099, 0.32704,
203 0.34406, 0.36101, 0.3779, 0.39471, 0.41147, 0.42818, 0.44484, 0.46147,
204 0.47807, 0.49468, 0.51134, 0.52803, 0.54474, 0.56138, 0.57804, 0.59478,
205 0.61171, 0.62881, 0.64609, 0.6635, 0.68111, 0.69901, 0.71722, 0.73583,
206 0.7549, 0.77447, 0.79471, 0.81564, 0.83759, 0.86067, 0.88557, 0.91159,
207 0.9452, 1.0
208};
209
210static constexpr int N_S_Horseshoe = 51;
211static constexpr double S_Horseshoe[51] = {
212 0.0, 0.00467, 0.01237, 0.02268, 0.03515, 0.04943, 0.06525, 0.08212, 0.10005,
213 0.11891, 0.13856, 0.15896, 0.18004, 0.20172, 0.22397, 0.24677, 0.27006,
214 0.2938, 0.3179, 0.34237, 0.3672, 0.39239, 0.41792, 0.44374, 0.46984,
215 0.49619, 0.52276, 0.5495, 0.5764, 0.60345, 0.63065, 0.65795, 0.68531,
216 0.71271, 0.74009, 0.76738, 0.79451, 0.82144, 0.84814, 0.8745, 0.90057,
217 0.92652, 0.95244, 0.97724, 0.99988, 1.02048, 1.03989, 1.05698, 1.07694,
218 1.07562, 1.0
219};
220
221static constexpr int N_W_Horseshoe = 26;
222static constexpr double W_Horseshoe[26] = {
223 .0, .5878, .8772, .8900, .9028, .9156, .9284, .9412,
224 .9540, .9668, .9798, .9928, .9992, .9992, .9928, .9798,
225 .9600, .9330, .8980, .8542, .8000, .7332, .6499, .5426,
226 .3919, .0
227};
228
229// ============================================================================
230// GOTHIC
231// ============================================================================
232
233static constexpr int N_Y_Gothic = 51;
234static constexpr double Y_Gothic[51] = {
235 0.0, 0.04522, 0.07825, 0.10646, 0.12645, 0.14645, 0.16787, 0.18641, 0.20129,
236 0.22425, 0.24129, 0.25624, 0.27344, 0.29097, 0.30529, 0.32607, 0.33755,
237 0.35073, 0.36447, 0.37558, 0.4, 0.4181, 0.43648, 0.45374, 0.46805,
238 0.48195, 0.49626, 0.51352, 0.5319, 0.55, 0.56416, 0.57787, 0.59224,
239 0.6095, 0.62941, 0.65, 0.67064, 0.69055, 0.70721, 0.72031, 0.73286,
240 0.74632, 0.76432, 0.78448, 0.80421, 0.82199, 0.84363, 0.87423, 0.90617,
241 0.93827, 1.0
242};
243
244static constexpr int N_S_Gothic = 51;
245static constexpr double S_Gothic[51] = {
246 0.0, 0.005, 0.0174, 0.03098, 0.04272, 0.055, 0.0698, 0.0862, 0.10461,
247 0.12463, 0.145, 0.16309, 0.18118, 0.2, 0.22181, 0.24487, 0.26888,
248 0.2938, 0.31901, 0.34389, 0.36564, 0.38612, 0.4072, 0.43, 0.45868,
249 0.48895, 0.52, 0.55032, 0.5804, 0.61, 0.63762, 0.66505, 0.6929,
250 0.72342, 0.75467, 0.785, 0.81165, 0.83654, 0.86, 0.88253, 0.90414,
251 0.925, 0.94486, 0.96475, 0.98567, 1.00833, 1.03, 1.0536, 1.065,
252 1.055, 1.0
253};
254
255static constexpr int N_W_Gothic = 21;
256static constexpr double W_Gothic[21] = {
257 0.0, 0.286, 0.643, 0.762, 0.833, 0.905, 0.952, 0.976, 0.976,
258 1.0, 1.0, 0.976, 0.976, 0.952, 0.905, 0.833, 0.762,
259 0.667, 0.524, 0.357, 0.0
260};
261
262// ============================================================================
263// CATENARY
264// ============================================================================
265
266static constexpr int N_Y_Catenary = 51;
267static constexpr double Y_Catenary[51] = {
268 0.0, 0.02974, 0.06439, 0.08433, 0.10549, 0.12064, 0.13952, 0.1556, 0.17032,
269 0.18512, 0.20057, 0.21995, 0.24011, 0.25892, 0.27595, 0.29214, 0.30802,
270 0.32372, 0.33894, 0.35315, 0.36557, 0.37833, 0.3923, 0.4097, 0.42982,
271 0.45, 0.46769, 0.48431, 0.5, 0.51466, 0.52886, 0.54292, 0.55729,
272 0.57223, 0.5878, 0.60428, 0.62197, 0.64047, 0.6598, 0.67976, 0.7,
273 0.71731, 0.73769, 0.76651, 0.8, 0.8209, 0.84311, 0.87978, 0.91576,
274 0.95, 1.0
275};
276
277static constexpr int N_S_Catenary = 51;
278static constexpr double S_Catenary[51] = {
279 0.0, 0.00605, 0.01455, 0.0254, 0.03863, 0.0543, 0.07127, 0.08778, 0.10372,
280 0.12081, 0.14082, 0.16375, 0.18779, 0.21157, 0.23478, 0.25818, 0.28244,
281 0.30741, 0.33204, 0.35505, 0.37465, 0.39404, 0.41426, 0.43804, 0.46531,
282 0.49357, 0.52187, 0.54925, 0.57647, 0.60321, 0.62964, 0.65639, 0.68472,
283 0.71425, 0.74303, 0.76827, 0.79168, 0.815, 0.84094, 0.86707, 0.89213,
284 0.91607, 0.94, 0.96604, 0.99, 1.00714, 1.02158, 1.03814, 1.05,
285 1.05, 1.0
286};
287
288static constexpr int N_W_Catenary = 21;
289static constexpr double W_Catenary[21] = {
290 0.0, 0.6667, 0.8222, 0.9111, 0.9778, 1.0000, 1.0000, 0.9889, 0.9778,
291 0.9556, 0.9333, 0.8889, 0.8444, 0.8000, 0.7556, 0.7000, 0.6333,
292 0.5556, 0.4444, 0.3333, 0.0
293};
294
295// ============================================================================
296// SEMI-ELLIPTICAL
297// ============================================================================
298
299static constexpr int N_Y_SemiEllip = 51;
300static constexpr double Y_SemiEllip[51] = {
301 0.0, 0.03075, 0.05137, 0.07032, 0.09, 0.11323, 0.13037, 0.14519, 0.15968,
302 0.18459, 0.19531, 0.21354, 0.22694, 0.23947, 0.25296, 0.265, 0.27784,
303 0.29212, 0.3097, 0.32982, 0.35, 0.36738, 0.3839, 0.4, 0.41667,
304 0.43333, 0.45, 0.46697, 0.48372, 0.5, 0.51374, 0.52747, 0.54209,
305 0.5595, 0.57941, 0.6, 0.62, 0.64, 0.66, 0.68, 0.7,
306 0.71843, 0.73865, 0.76365, 0.7926, 0.82088, 0.85, 0.88341, 0.90998,
307 0.93871, 1.0
308};
309
310static constexpr int N_S_SemiEllip = 51;
311static constexpr double S_SemiEllip[51] = {
312 0.0, 0.00438, 0.01227, 0.02312, 0.03638, 0.05145, 0.06783, 0.085, 0.10093,
313 0.11752, 0.1353, 0.15626, 0.17917, 0.20296, 0.22654, 0.24962, 0.27269,
314 0.29568, 0.31848, 0.34152, 0.365, 0.38941, 0.41442, 0.44, 0.46636,
315 0.49309, 0.52, 0.54628, 0.57285, 0.6, 0.62949, 0.65877, 0.68624,
316 0.71017, 0.73304, 0.75578, 0.77925, 0.80368, 0.83114, 0.8595, 0.88592,
317 0.90848, 0.93, 0.95292, 0.97481, 0.99374, 1.01084, 1.02858, 1.04543,
318 1.05, 1.0
319};
320
321static constexpr int N_W_SemiEllip = 21;
322static constexpr double W_SemiEllip[21] = {
323 0.0, 0.7000, 0.9800, 1.0000, 1.0000, 1.0000, 0.9900, 0.9800, 0.9600,
324 0.9400, 0.9100, 0.8800, 0.8400, 0.8000, 0.7500, 0.7000, 0.6400,
325 0.5600, 0.4600, 0.3400, 0.0
326};
327
328// ============================================================================
329// BASKETHANDLE
330// ============================================================================
331
332static constexpr int N_A_Baskethandle = 26;
333static constexpr double A_Baskethandle[26] = {
334 .0000,.0173,.0457,.0828,.1271,.1765,.2270,.2775,.3280,.3780,
335 .4270,.4765,.5260,.5740,.6220,.6690,.7160,.7610,.8030,.8390,
336 .8770,.9110,.9410,.9680,.9880,1.000
337};
338
339static constexpr int N_R_Baskethandle = 26;
340static constexpr double R_Baskethandle[26] = {
341 .0100,.0952,.1890,.2730,.3690,.4630,.5600,.6530,.7430,.8220,
342 .8830,.9490,.9990,1.055,1.095,1.141,1.161,1.188,1.206,1.206,
343 1.206,1.205,1.196,1.168,1.127,1.000
344};
345
346static constexpr int N_Y_BasketHandle = 51;
347static constexpr double Y_BasketHandle[51] = {
348 0.0, 0.04112, 0.0738, 0.1, 0.12236, 0.14141, 0.15857, 0.17462, 0.18946,
349 0.20315, 0.21557, 0.22833, 0.2423, 0.25945, 0.27936, 0.3, 0.3204,
350 0.34034, 0.35892, 0.37595, 0.39214, 0.40802, 0.42372, 0.43894, 0.45315,
351 0.46557, 0.47833, 0.4923, 0.50945, 0.52936, 0.55, 0.57, 0.59,
352 0.61023, 0.63045, 0.65, 0.66756, 0.68413, 0.7, 0.71481, 0.72984,
353 0.74579, 0.76417, 0.78422, 0.80477, 0.82532, 0.85, 0.88277, 0.915,
354 0.95, 1.0
355};
356
357static constexpr int N_S_BasketHandle = 51;
358static constexpr double S_BasketHandle[51] = {
359 0.0, 0.00758, 0.01812, 0.03, 0.03966, 0.04957, 0.0623, 0.07849, 0.09618,
360 0.11416, 0.13094, 0.14808, 0.16583, 0.18381, 0.20294, 0.225, 0.2547,
361 0.28532, 0.31006, 0.32804, 0.34555, 0.36944, 0.40032, 0.43203, 0.46004,
362 0.47849, 0.49591, 0.51454, 0.5381, 0.56711, 0.6, 0.64092, 0.68136,
363 0.71259, 0.73438, 0.755, 0.78625, 0.8188, 0.85, 0.8679, 0.88483,
364 0.90431, 0.9369, 0.97388, 1.00747, 1.033, 1.05, 1.05464, 1.06078,
365 1.055, 1.0
366};
367
368static constexpr int N_W_BasketHandle = 26;
369static constexpr double W_BasketHandle[26] = {
370 0.0, 0.49, 0.667, 0.82, 0.93, 1.00, 1.00, 1.00, 0.997,
371 0.994, 0.988, 0.982, 0.967, 0.948, 0.928, 0.904, 0.874,
372 0.842, 0.798, 0.75, 0.697, 0.637, 0.567, 0.467, 0.342,
373 0.0
374};
375
376// ============================================================================
377// SEMI-CIRCULAR
378// ============================================================================
379
380static constexpr int N_Y_SemiCirc = 51;
381static constexpr double Y_SemiCirc[51] = {
382 0.0, 0.04102, 0.07407, 0.1, 0.11769, 0.13037, 0.14036, 0.15, 0.16546,
383 0.18213, 0.2, 0.22018, 0.2403, 0.25788, 0.27216, 0.285, 0.29704,
384 0.30892, 0.32128, 0.33476, 0.35, 0.36927, 0.38963, 0.41023, 0.43045,
385 0.45, 0.46769, 0.48431, 0.5, 0.51443, 0.52851, 0.54271, 0.55774,
386 0.57388, 0.59101, 0.60989, 0.63005, 0.65, 0.66682, 0.68318, 0.7,
387 0.71675, 0.73744, 0.76651, 0.8, 0.8209, 0.84311, 0.87978, 0.91576,
388 0.95, 1.0
389};
390
391static constexpr int N_S_SemiCirc = 51;
392static constexpr double S_SemiCirc[51] = {
393 0.0, 0.00757, 0.01815, 0.03, 0.0358, 0.04037, 0.04601, 0.055, 0.07475,
394 0.09834, 0.125, 0.1557, 0.18588, 0.20883, 0.223, 0.23472, 0.24667,
395 0.26758, 0.29346, 0.32124, 0.35, 0.3772, 0.4054, 0.43541, 0.46722,
396 0.5, 0.53532, 0.56935, 0.6, 0.61544, 0.62811, 0.6417, 0.66598,
397 0.7001, 0.73413, 0.76068, 0.78027, 0.8, 0.82891, 0.85964, 0.89,
398 0.9127, 0.93664, 0.96677, 1.0, 1.02661, 1.04631, 1.05726, 1.06637,
399 1.06, 1.0
400};
401
402static constexpr int N_W_SemiCirc = 21;
403static constexpr double W_SemiCirc[21] = {
404 0.0, 0.5488, 0.8537, 1.0000, 1.0000, 0.9939, 0.9878, 0.9756, 0.9634,
405 0.9451, 0.9207, 0.8902, 0.8537, 0.8171, 0.7683, 0.7073, 0.6463,
406 0.5732, 0.4756, 0.3354, 0.0
407};
408
409// ============================================================================
410// HORIZONTAL ELLIPSE
411// ============================================================================
412
413static constexpr int N_A_HorizEllipse = 26;
414static constexpr double A_HorizEllipse[26] = {
415 .0000,.0150,.0400,.0650,.0950,.1300,.1650,.2050,.2500,.3000,
416 .3550,.4150,.4800,.5200,.5850,.6450,.7000,.7500,.7950,.8350,
417 .8700,.9050,.9350,.9600,.9850,1.000
418};
419
420static constexpr int N_R_HorizEllipse = 26;
421static constexpr double R_HorizEllipse[26] = {
422 .0100,.0764,.1726,.2389,.3274,.4191,.5120,.5983,.6757,.7630,
423 .8326,.9114,.9702,1.030,1.091,1.146,1.185,1.225,1.257,1.274,
424 1.290,1.282,1.274,1.257,1.185,1.000
425};
426
427static constexpr int N_W_HorizEllipse = 26;
428static constexpr double W_HorizEllipse[26] = {
429 .0,.3919,.5426,.6499,.7332,.8000,.8542,.8980,.9330,.9600,
430 .9798,.9928,.9992,.9992,.9928,.9798,.9600,.9330,.8980,.8542,
431 .8000,.7332,.6499,.5426,.3919,.0
432};
433
434// ============================================================================
435// VERTICAL ELLIPSE
436// ============================================================================
437
438static constexpr int N_A_VertEllipse = 26;
439static constexpr double A_VertEllipse[26] = {
440 .0000,.0100,.0400,.0700,.1000,.1400,.1850,.2300,.2800,.3300,
441 .3800,.4300,.4800,.5200,.5700,.6200,.6700,.7200,.7700,.8150,
442 .8600,.9000,.9300,.9600,.9900,1.000
443};
444
445static constexpr int N_R_VertEllipse = 26;
446static constexpr double R_VertEllipse[26] = {
447 .0100,.1250,.2436,.3536,.4474,.5484,.6366,.7155,.7768,.8396,
448 .8969,.9480,.9925,1.023,1.053,1.084,1.107,1.130,1.154,1.170,
449 1.177,1.177,1.170,1.162,1.122,1.000
450};
451
452static constexpr int N_W_VertEllipse = 26;
453static constexpr double W_VertEllipse[26] = {
454 .0,.3919,.5426,.6499,.7332,.8000,.8542,.8980,.9330,.9600,
455 .9798,.9928,.9992,.9992,.9928,.9798,.9600,.9330,.8980,.8542,
456 .8000,.7332,.6499,.5426,.3919,.0
457};
458
459// ============================================================================
460// ARCH
461// ============================================================================
462
463static constexpr int N_A_Arch = 26;
464static constexpr double A_Arch[26] = {
465 .0000,.0200,.0600,.1000,.1400,.1900,.2400,.2900,.3400,.3900,
466 .4400,.4900,.5400,.5900,.6400,.6900,.7350,.7800,.8200,.8600,
467 .8950,.9300,.9600,.9850,.9950,1.000
468};
469
470static constexpr int N_R_Arch = 26;
471static constexpr double R_Arch[26] = {
472 .0100,.0983,.1965,.2948,.3940,.4962,.5911,.6796,.7615,.8364,
473 .9044,.9640,1.018,1.065,1.106,1.142,1.170,1.192,1.208,1.217,
474 1.220,1.213,1.196,1.168,1.112,1.000
475};
476
477static constexpr int N_W_Arch = 26;
478static constexpr double W_Arch[26] = {
479 .0,.6272,.8521,.9243,.9645,.9846,.9964,.9988,.9917,.9811,
480 .9680,.9515,.9314,.9101,.8864,.8592,.8284,.7917,.7527,.7065,
481 .6544,.5953,.5231,.4355,.3195,.0
482};
483
484// ============================================================================
485// Standard ellipse and arch size tables
486// ============================================================================
487
488static constexpr int NumCodesEllipse = 23;
489
490static constexpr double MinorAxis_Ellipse[23] = {
491 14,19,22,24,27,29,32,34,38,43,48,53,58,63,68,72,77,82,87,92,97,106,116
492};
493
494static constexpr double MajorAxis_Ellipse[23] = {
495 23,30,34,38,42,45,49,53,60,68,76,83,91,98,106,113,121,128,136,143,151,166,180
496};
497
498static constexpr double Afull_Ellipse[23] = {
499 1.80,3.30,4.10,5.10,6.30,7.40,8.80,10.20,12.90,16.60,20.50,24.80,29.50,34.60,
500 40.10,46.10,52.40,59.20,66.40,74.00,82.00,99.20,118.60
501};
502
503static constexpr double Rfull_Ellipse[23] = {
504 0.367,0.490,0.546,0.613,0.686,0.736,0.812,0.875,0.969,1.106,1.229,1.352,
505 1.475,1.598,1.721,1.845,1.967,2.091,2.215,2.340,2.461,2.707,2.968
506};
507
508static constexpr int NumCodesArch = 102;
509
510static constexpr double Yfull_Arch[102] = {
511 11,13.5,15.5,18,22.5,26.625,31.3125,36,40,45,54,62,
512 72,77.5,87.125,96.875,106.5,
513 13,15,18,20,24,29,33,38,43,47,52,57,
514 31,36,41,46,51,55,59,63,67,71,75,79,83,87,91,
515 55,57,59,61,63,65,67,69,71,73,75,77,79,81,83,85,87,89,91,
516 93,95,97,100,101,103,105,107,109,111,113,115,118,119,121,
517 112,114,116,118,120,122,124,126,128,130,132,134,
518 136,138,140,142,144,146,148,150,152,154,156,158
519};
520
521static constexpr double Wmax_Arch[102] = {
522 18,22,26,28.5,36.25,43.75,51.125,58.5,65,73,88,102,
523 115,122,138,154,168.75,
524 17,21,24,28,35,42,49,57,64,71,77,83,
525 40,46,53,60,66,73,81,87,95,103,112,117,128,137,142,
526 73,76,81,84,87,92,95,98,103,106,112,114,117,123,128,131,137,139,
527 142,148,150,152,154,161,167,169,171,178,184,186,188,190,197,199,
528 159,162,168,170,173,179,184,187,190,195,198,204,
529 206,209,215,217,223,225,231,234,236,239,245,247
530};
531
532static constexpr double Afull_Arch[102] = {
533 1.1,1.65,2.2,2.8,4.4,6.4,8.8,11.4,14.3,17.7,25.6,34.6,
534 44.5,51.7,66,81.8,99.1,1.1,1.6,2.2,2.9,4.5,6.5,8.9,11.6,14.7,18.1,
535 21.9,26,7,9.4,12.3,15.6,19.3,23.2,27.4,32.1,37,42.4,48,54.2,60.5,
536 67.4,74.5,22,24,26,28,31,33,35,38,40,43,46,49,52,55,58,61,64,67,
537 71,74,78,81,85,89,93,97,101,105,109,113,118,122,126,131,97,102,
538 105,109,114,118,123,127,132,137,142,146,151,157,161,167,172,177,
539 182,188,194,200,205,211
540};
541
542static constexpr double Rfull_Arch[102] = {
543 0.25,0.3,0.36,0.45,0.56,0.68,0.8,0.9,1.01,1.13,1.35,
544 1.57,1.77,1.92,2.17,2.42,2.65,0.324,0.374,0.449,0.499,0.598,0.723,
545 0.823,0.947,1.072,1.171,1.296,1.421,0.773,0.773,1.022,1.147,1.271,
546 1.371,1.471,1.570,1.670,1.770,1.869,1.969,2.069,2.168,2.268,1.371,
547 1.421,1.471,1.520,1.570,1.620,1.670,1.720,1.770,1.820,1.869,1.919,
548 1.969,2.019,2.069,2.119,2.168,2.218,2.268,2.318,2.368,2.418,2.493,
549 2.517,2.567,2.617,2.667,2.717,2.767,2.817,2.866,2.941,2.966,3.016,
550 2.792,2.841,2.891,2.941,2.991,3.041,3.091,3.141,3.190,3.240,3.290,
551 3.340,3.390,3.440,3.490,3.539,3.589,3.639,3.689,3.739,3.789,3.838,
552 3.888,3.938
553};
554
555} // namespace openswmm::xsect_tables
556
557#endif // OPENSWMM_XSECT_TABLES_HPP
openswmm::xsect_tables
Definition xsect_tables.hpp:43