1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
|
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Overview\n",
"\n",
"This a notebook that inspects the results of a WarpX simulation.\n",
"\n",
"# Instruction\n",
"\n",
"Enter the path of the data you wish to visualize below. Then execute the cells one by one, by selecting them with your mouse and typing `Shift + Enter`"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"# Import statements\n",
"import yt ; yt.funcs.mylog.setLevel(50)\n",
"import numpy as np\n",
"import scipy.constants as scc\n",
"import matplotlib.pyplot as plt\n",
"%matplotlib notebook"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Read data in the simulation frame"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/Users/mthevenet/anaconda3/lib/python3.7/site-packages/yt/units/yt_array.py:1394: RuntimeWarning: invalid value encountered in true_divide\n",
" out=out, **kwargs)\n"
]
},
{
"data": {
"text/html": [
"<img style=\"max-width:100%%;max-height:100%%;\" src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAw0AAAHGCAYAAAA2UWe+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzsvXu8XVV57/17xpxr7Z3snUASNCok8S0g1hCgRRQhYCUgpaD02OoHjsf3VYsYKhRrEUUEIpUoeKE9tVZPCyq1otijgKWAclO8BMgRBHpASy9CAGNzI9nJXpc5x/P+Ma5zrstee++1k315fnzCXmuu7xpjzDHmnGs8YzzPGMTMDJFIJBKJRCKRSCTqILWvCyASiUQikUgkEommt8RoEIlEIpFIJBKJRF0lRoNIJBKJRCKRSCTqKjEaRCKRSCQSiUQiUVeJ0SASiUQikUgkEom6SowGkUgkEolEIpFI1FViNIhEIpFIJBKJRKKuEqNBJBKJRCKRSCQSdZUYDSKRSCQSiUQikairxGgQiUQikUgkEolEXSVGg0gkEolEIpFIJOoqMRpEIpFIJBKJRCJRV4nRIBKJRCKRSCQSibpKjAaRSCQSiUQikUjUVWI0iEQikUgkEolEoq4So0EkEolEIpFIJBJ1lRgNIpFIJBKJRCKRqKvEaBCJRCKRSCQSiURdle7rAuxtXXHFFfjlL3+JL3/5y4XjN910E+666y4ccsgheOihh3DOOefg1FNPHTcjEolEIpFIJBLNNs0po+HBBx/EJz/5SZx99tmF49dddx3WrVuHJ598EkNDQ9i0aRNWrlyJ2267DatXr+6ZEYlEIpFIJBKJZqOImXlfF2JvaM+ePVi7di3+4z/+AwcffLCfadi1axeWL1+O8847D+vXr/f8WWedhZ///Od4+OGHe2LKIiKQSiZdbmYGEYF1DpACEdnjGmAu5DFexjQ9TzqdMiNlljJLmfdxmXXu/86YMs/EepYyS5lnYZknqze+YX9s2daYdDplHfDi43DHHXf0PV1R75ozMw2XX345LrvsMrznPe8pHL/jjjuwY8cOHHvssYXjxx57LL7xjW/g8ccfxxNPPDEmc/jhhxc+I5XgjL/4xaTL/bOvfwQHv/F83LvuBMw/4OVYc9n3AADfef8rAM59HiPbnhs387OvfwRPb/jHSadTZn7wqTPx23/0NzOqzDOxnqeyzD/41Jl44dknZlSZZ2I9T1WZ/+n9r8DvXH7fjCrzTKxnuQdnfj3L72Cxno88KwyMTlTPfekoPPT9Y8cGx6lj3rCl72mKxqc5EQh9++234+CDD8ahhx7a8tkjjzwCAFi+fHnh+IoVKwAADz30UE/MVGnpyjUAzGQQ51n4QJfJ8TNLV67pSzp7k5mqMs/Eep7KMk8VI/Us9+Bsqme5B2d+Pcs9GF6aeu6DmKG52fd/on2vWW80bN26FTfddBPOO++8tp9v2WIs1wULFhSODw8PAwA2b97cEzNVesmqNQDsdKGKpg2p7FU2fuYlq9b0JZ29yUxVmWdiPU9lmaeKkXqWe3A21bPcgzO/nuUeLNezSNRZs9496ZJLLsFVV13V8fOBgQEAaPHjU8rYU9VqtSemLGbGDz51pn+/4rizsOL4s1u43mRHBzSXD01LZsVxZ02r8ggzfmbFcWfh0W9cPm3KI8z4GFUZnFblEUbuwbnGzNzfwYnplz+6Eb/88dcBAAcvAXJkY3xDNBM1q42Gr3zlKzjppJPwspe9rCOzdOlSAMDIyEjhuHt/4IEHYnR0dEymLCLCiR+8ZeKFj2WnDwnx6ACKN/s0YlYcfzZGtjw7bcojzPiZFcefjUdvunzalEeY8TFJZXBalUeY8TNyD85sZsb+Dk5QK44/2w+MPnv9EdBiNMxKzWr3pC9/+cs455xzMDw87P/df//9+Id/+AcMDw9j/fr1OProowEAzzzzTOG7Tz/9NADgqKOO6omZUtlW4papRmGEmRpGR6wmNv8S9ow7JkzvTLt6bithhBFGmH3FiERdNKtnGr761a/6WQKnt7/97TjooINw9dVXY/HixVi4cCGWLFmCBx54AKeffrrnNmzYgGOOOQaHHXYYDj744DGZqZWK/m/Vcs8LI0z/GMXkGcVk/yaecceEGR8z3dt9XzLG0DL1RBQZWhw+N3+FmamM1zS+DoWZvBgs7kmzVGpsZObqwAMPxCGHHFL4N2/ePCxYsACHHHIIFi9ejDRNcdlll+GGG27wBsbzzz+PW2+9FZdfbvxJe2GmUjZ0ArEbYjy7KIwwwghTYAiAnfnQHM18WMa/n0ZMfF4Fn+sO5y7MzGP21vUzLe7BGcr0Q8Zo6P9/on2vWT3T0KsuvPBCDAwM4Nxzz8WqVauwceNG3HDDDTjjjDPGxUyVtFsSrYuJNxuYdu4cmdt7UBU/F0YYYbowDCC3vQFF0UyHPeTeCyPMLGKm+2/cTGBEom6ac0bDfffd1/b42rVrsXbt2q7f7YWZSrG/4wF02Mh7JjPGJYbD6+gvax25fAgjjDBdGQrGBLOOXEfsaC0V/wojzFxgYmNjuv4OThdmsspbN4nog8Ta2deac0bDzJS5+QpLvpICOBdmLzJmSjz3jP9Rskz44RJGmH3MMKBsZ4AQjcxaJhgfwgizbxlz/WrPeHV6Pk+Cmapn+JjlmVHM5MWYKqNBtK8lRsOMkLGuFUVWdsvowN5hzEOXPNN5xGf6MOaHavLnrpjgg11JdRkFFkaYfczQ9LoHhRGmGzPe5/M+v79KTC9lnjnM5MVgZFMRYS3a5xKjYYrEzPjZ1z+CpSvX9GGXRWOx63KwGe99xnSc2TNh5Ah2dLM4ujRdmOlUh8IIs1cZFkaY2cWE0X/DlEf/Zyoz5rlP8XMjGHGdFbt4lfWrx+7G5n+5G/uPmYpopkqMhikSEeHIs9YXjvVyQ7aTv7/je7WclJ4bjHO98K9RrJ8wkiWMMHOYYUDl4cYKo6XF98III8zYTPy746XbdJ7HYBRT17w6p2P/TvFvbjeDoBe9ZJUZJP3l9bcgl5mGWSkxGvaSJmowAABU+G7H0Ypk7jHt6qfjbIR/FpIwwsx+hqbPfSqMMHuDUUwIcbJR57fl92L8jLun+nmfTqY8U8r0QQz0dYHUu27fhbtu34UXXljSx1RFE5EYDXtJk7Hgne8ktUmn7F8515nOUpbp9nQURpjZyHSTMMLMDkZTtDkgijuyu8/N35nFmN+5qXluxPlNZmBzsjMUY+nk0xbg5NMW4C0nDU5pPqKxJUbDDFL3HTWFEUYYYdpBwggjzExljDEUZlY6uixOgCm7S03KI6JU/GwME1A0MyVGwxSqXzegC5TKdbfRivExXXIz/xdGGGGEEUYYYfYa03k23Xyfp5CJuX5I9m+enVJjI6J9Lhu8VFxbeXyMYvIBV72kA2GEEWZWMJN5bggjjDB7j/GzCpbxA4B787khEnWRzDRMofpltSvn88jUeQShB8aZiNStXMIII8ysYiCMMMLMWEYThyDwNi5ELYHiHZg4ULzTc2Oi3hHlPgeDxD1plkqMhhkhf6eHQy33tjDCCCNMK0PCCCPMnGfiQHHEcRIcPp+oyt9lAJk4ssxKidEwRWq3uduEb0o7gqALabBPUzFB2ftTF75YzE8YYYQRRhhhhJl7jGJCao0GjXh2gP3nk5Xb3G140imJpquIeQr2EBeBiFBd0J81hVlrNHdvB0ihOrwIANDYtRUAfB4tDAONEcsMW4YjZsimI4wwwsw6prl7B9L5C8f/3BBGGGFmPjOJZ0s/dPAhe/C177+yb+k5nXOSwsaNG/uerqh3idEwRVJJijP+4hfj+k6nmYiRbc/j+1ecgMEly7Fm3T1QTPjOhYcCrPGm//lvnrl33WrMW7IcJ19xLwAII4wwc5S589LX4HUX3Vx4bgDAbX/yCoA1Tv+rp3w6wggjzOxnFFNPz5Z+6KnrV+GGKTAa1p5EYjTsY4l70hSqX0uuumlH4uCSFPshxu5JiG3AUvbCCCPM3GNIGGGEEQa9PTf6IcZUxTRMQWFF45JEqswAuTWeucsyacIII4wwwggjjDCTYUSibpKZhhkhe6N3HR2wB4QRRhhhogO9PDeEEUaYucNM/UwDIUPS/4SRTUGaovFIjIYpVN/2abCbsimK9mAgABzl0WkTmfiBIIwwwswhhoQRRhhh7McMvzE1RXs5WKZf7tSwWebiyDIrJUbDPlSvN6nbtIXBYZOWsvz92SVNYYQRZg4xLIwwwgjTlum4CewE1E+DQzS9JUbDFKp/N5LrIVBINwqKNjd7kQFQnIJsk44wwggjjDDCCDP7GcUExYkluqUzfrUzOGRzt9kpMRqmSMyMR2+8FEsPPwlL7eZuE5fxRWCto5vTzCmG94EJapm/FEYYYYQRRhhh5hhjdoTOzWutowHHcjoTl9vcbQBTFdMg2tcSo2GKREQ46qz15s0k70fFbi21yHJvuc+N8cDCCCOMMBHTy3NDGGGEmXuMJvbMZDwj3ODlS1atwUtWrcGT19824bRE01tiNEyh+uWepO30YRzf1GKIaGGEEUaYmIEwwggjTEdGMXmmXwu3uGzFPWl2SoyGmSA30YD2Kx4oJkCRZ7zKowzCCCPMHGIgjDDCCGM+ZoKi8Nwor8TYTzEIOYt70myUGA1TqLLlXp55GOtzxzj3JEIUw9AyOkCe8WpJThhhhJlrDIQRRpg5z5iYhvCB729w+Lyfsw0y0zA7JUbDFCo2Akznv/sN2flzM6eoC6sclJ8YwggjjDCtDAsjjDDC9MBM1KW6n8aGaHpLjIa9pMndjG7eMZ52JBSXShNGGGGEacMoYYQRRpixmX51/qduR2jRvpYYDVOoXtyPxqdevi+MMMIIEyPCCCOMMONkxqFy36bfRsNPb9+Mn96xGS+8MNC3NEUTkxgNe1ETteKVG1SMv8/CCCOMMN0Y+5kwwggjTA/MdN3Z+bdPW4rfPm0pPrvmV/u6KHNeYjRMkZgZj3z9I+Pe3K2dYeH3aulicwgjjDDCFBnTAYiXbRdGGGGE6cbEGu9Ap9vcTUECoWerxGiYIhERjjj7qnF/r/0KSy7oses3hRFGGGFaGS2MMMIIMzYz2ZgGt7nbo9ffKUuuzlKJ0TCFmkhMQ9ub1t7ghU1byvOOwggjjDAFhoQRRhhhembGWha+k6arW5Oo/5r180fMjE9/+tM4+OCDMTQ0hNWrV+OnP/1pC3fTTTfh3HPPxTXXXIO3vvWtuPPOOyfEdJNbdrXbv/ZfNH+osFKKEkYYYYTpwrCDhBFGGGF6YyKZvR3G/tdOGVTf/4n2vWb9TMM111yDoaEhbNiwAVu3bsUf/dEf4c1vfjOeeeYZ/4N73XXXYd26dXjyyScxNDSETZs2YeXKlbjtttuwevXqnpmyxrK+e58KtDdLnByX5xeFEUYYYYQRRhhhJsN0Vy/9FllydfZqVptuo6OjOPTQQ3H++efjRS96EV75ylfiQx/6EJ599lls3boVALBr1y5cdNFFeMc73oGhoSEAwEEHHYTTTjsNF1xwQc/MRNSrxe58EYqflm9cYYQRRphWpkgJI4wwwnRjgnr2hhDNGRFznxfsneb6yle+gs9//vN44IEHAADf/OY38ba3vQ233HIL3vzmN3vuL/7iL/Cnf/qneOyxx/DEE0+MyRx++OGFfIgI1eElfSkzs0Zz93aAFKpDiwAAjRFj9Lg8hBFGGGEc09y9A+n8hYEZtswuyyyw6WgtjDDCCOOZfmjpwRofvOfUvqXn9JVT/hUbN27se7qi3jXr3ZNibd68Gd/61rfw7W9/2x975JFHAADLly8vsCtWrAAAPPTQQ3jqqafGZFqMBpXg1PUP9qXcI1uexb1XnojBxQfhlCvuBQB8508OAcA+D2GEEUYYx9x56Wvwug98C9//2OsxuPggrFl3DwDgtgsOBcA4Zf0DPh1hhBFGGMf0Qz+77jXIZ7cjy5zVnDAaduzYgWuvvRaf+9znsG3bNgwNDeGGG25AmqbYsmULAGDBggWF7wwPDwMwhkYvTDuNd0WBsQKhizOKhMKujsIII4wwwggjjDCTYUqaqEuSxDTMTs0Jo2G//fbDxRdfjLe85S349Kc/ja9+9as47rjjcP7552NgwGxLXliFBICyW65Wq9WemLKYGT+85vf9++XHn4Xlx5/VtZyauMMNSoU/NgNhhBFGGGGEEUaYPjJB4zEYfvmjG/HLH38dANBHTyfRNNOcMBqICENDQzjyyCPx93//9/j3f/933HPPPTj//POxdOlSAMDIyEjhO+79gQceiNHR0TGZdnme+MFbigc735tjyH6xsJMTlRIURhhhhBFGGGGEmQwTNB5viWWrz8Ky1WZg9JG/ey0yFvek2ag52arHH3+8nx04+uijAQDPPPNMgXn66acBAEcddVRPzNTKNBOpqLlaBgCEEUYYYdoxJIwwwgjTIzN5MYx7Ur//ifa95qTRsGnTJpxyyikAgJNPPhlLlizxqyk5bdiwAccccwwOO+ywnph26nVjlO7LrQJhecVu6ykLI4wwwrQyusMoojDCCCPMWOplU1pZlnXuaFYbDTt27MAHPvAB/PjHP/bHHnzwQWzfvh3vfOc7AQBpmuKyyy7DDTfc4N2Qnn/+edx66624/PLLe2amVvYmjPsHnXwVhRFGGGEihoQRRhhhemYmLwYhZ9X3f6J9r1kd09BoNHDffffhc5/7HE455RSccMIJOOCAA3DLLbcgScJU14UXXoiBgQGce+65WLVqFTZu3IgbbrgBZ5xxxriYsvptcRenFBXAuTDCCDMOhpg8QxwdA0BK+dcxw9Es4HQ9r9nAMLGva1II9W4Z9z5m9nWZhRFmtjJO410F0knciWanZrXR8OIXvxg//elPe2LXrl2LtWvXTpqZShUcETqMDggjzFxgiAnut4xBobNvmWAQdGaihKP/9yeviZ7XbGGY2DPlzr6OXpeZkE4vdUgd0xkrL4oGdKZrHQojzL5kRKJ2mtVGw75W3/ZpsCsdUPFOF0aYWckQk/+MNKIOXrGT7r9SWAmklNc+YtqVmYmnVT33i2FiMJsPiWMDoEM6RetsnzBjldkbFtOonoURZq8wJU3EY8K5J4lmn8RomAlSpb8AOm7aIoww05ghprDfiYoMAMuE9z3kNcOYdufOxHuhLbhvTHAZMox/79wdqFs6M4thYrDzsOhw7sQ0JfUsjDD7jOmDGOKeNFslRsMUiZnx2I2XYunKNXjJqjWTTE1F//cZCCPMtGKMG09gOrnxuN5197zmBkNcNCpMx7S/7YWurj7dGRNnQF2Z4mdziynXT2dXqMm3hTDC7B2mqPF4TGx+7G5sfvweQGYaZq3EaJgiERFWnX0VAPS8vFnnaUAzl5h3vbmFEWbvMVS6VsvuOLprOv1nSJfKY9+7aXitMaUMK+7LeRnDq5d0em8vHgcTxxaE15bpYSnHTgwrBqvAsOKWz83fqWHKbdav8+qFMYZFHspDrfdPIZ0+t6kwwoyPKWo87kkvPfxkvPTwk/Hgda/r+TuimSUxBWeETDN51wYAxaVDhBFmahnTmTXTzWovl4e0+ecY9560LQ/GTideLGQqmHJ52pW5fF5j5TV9mKQrYwwCl457z54JHfvAjF2e/jLl8rQrc/t0up/7RBnjlqbC60KZp0u7CzM3mcmLIZu7zVbJTMMMUtHe7zBSKYwwfWDMSChHrzul1GnEqjcmdKZtXn5kv41bSMeR3JnBeEPCvBszHcQrNZUDqtuk0jmd9ky8xKl73z4dLv6JZlXKI/3gHs5rmjLtzsufHyOaseheP+3z6s6U22LMMgsjzF5lxicJhJ69EqNhH6r3aT83pRgfK39XGGEmxsQuMWC0xB6U0ykmND7GdJy5KzORvMplDsug8l5hTIdvnGX2bk0TP/de2h09Mb2VORgJHRi0T6edi07MhKVSp4YxbTa+Mrdj4vYKRsTk74sy0+68zDmMp02FEWYiTOkbJWNWdnye2xKjYR+q95uxzZRihwAnYYQZiykH3xrZv13T6c4Yg0B5xo+qu7zKo+yd0mECIUqn074IEdOxzNTDefWB6aXMpiNYTIc0FeqnYETs1WtDtWVYsffJ6sgQh34HqY57JxC6pzOVTLt9I8Yqs2nT0rVWPvdS/YSZiWBkmOu+ez0bjdEWUdD1dH22CDMbmNI3JrLkKgOZzDTMSonRMCNkRgc0a2GEmRTTPgATcENP3dMJTKeA4XGn0ymg2gUD6/bpxHs5aN05MHtvMXGntGuZGZ4puxqZfKLz2pvXhpu5iNyNvIuOay/olvL6kX03Et8hLyaGJu2ZTunsLYaYxlVmd+7t7p+4ftqmEwWBa+gxA7M7pVNmxgqono7PH2FmAtMPEXKJQZiVElNwJsje38UAVGGEGZsxnWsaOx0eIx1NoHz86RBTGDhnKpYHvaXjvhenU2Si73VMJ5wHotfmvMz7hMKxmOklr3I9j13m9vXTjhkrr47nXmLiGQHfkXYxJB3awmyCZmcKOuRVZNyIeLu84Jm2ZY4MLyI27wuzGNzKdDz37nm1K3P3dNrXYS/10zYdFaWjemuLdky3Nh3PtSGMMCJRL5KZhpkga9p1XfZPGGGAwshjeaO0rumUmHgk1L92TIfp7Nido+j73z4voMvI/hi/ZKSp0MnuuFRqFybA5Irct7zaphPXTzdF9dxxedeovdrNWACA28CpU5AtKwanEVMOavbpmD9+NqJd8LRnXOc77tDb1wkHplyelnS6KGY6pTNWXtGmbK7M7c6LfTpdRvaj+uk0qzHWvVxui45S3RmzIZ0dTZ5EOsLMYaYPYkACoWepxGiYIjEzHr3xUiw9/CQsnfTmbk6xr6swwhSZsTqt4PbpFFb1YdUyut6pPCY2AuG1edG1zKbjjO6MDm48BGopTzBiwgivKo/S2/cuB8LUMtqNgtvytCuz6exzSMm7NbWvH+ql3Tu0aa9MGMSPRq4tE2IRwnkG3/84Efaj/goUzjFi4nNXRIVVifznNq2pYEzsyDjKbA0Lf+5oXz9xexV8v+Nzj6+NuAxd2qLdeRUM+R7anVgFo2WKrh9hZjEzAf3qsbux+V/uBoMkpmGWSoyGKRIR4Qi7uVsfUov+71QeixNmrjHxDsw9pUMTZ+JAX0I8o1FORxW+EzPhfdQF1yXGGwQ2gNWWS0Urx4SOvPGZ5cLob1RmHfaW6JbXRBjTqQt5KaZCebwhMVZ5VLF+4h2hi5p4m6Ir406gPVMIclaIer3O38YZFiGvMGNRjJOIzcWWjePauAdNlGk3oxO3VzAs0JJeoX6IQ/0SwgxFl/ppWeUoJFZg2pXZq0NbFNIpWB/dmXhlr5aVpCZ9/Qgz+5iJ6yWr1uAlq9bgR3+3Grl4v89KEY9na0BRzyIiVBcs6Q71WPPMGs3d2wFSqA4tAgA0RrYCAKrDS4QRZgoYQmNki2UOsEweMYtBAOqWGbCM1jmaewwzMLQYmCDDOkfDMfMts9syQ6E8M4qh0nn1UD9VyxTbgsds0+buHUjnLwzMsGV2WcY+m1jrHpgczd07PMMAmpapLFgCgAHWaI4YpjK8PwCguWubZRbbc9fIHLPAMjsts3CxL8+UMFwqj4FKZabSeQGIzr0yvAjUQ/30UocTbwthhJlaZkz10G8Z/o0qfveO9/SW3jj08zfdho0bN/Y9XVHvEqNhiqSSFL/3lz/vzvS4lNnIlmdx75UnYnDxMpyy7j4AwHf+5BAAjDf9z38TZhYyxISRrc/initPwLzFy3DKFd8HANx64cEAGG/+y3836YzBkDbp3P3x1Zi3aDneeLlhbvnTgwFonHntf/i87rrqeMxbtAKnXmaYmz/wGwA0fv+z/+mZ760/DvMWrcDvfvQHAIBv/9n/A0Djv33mlzavTfjuJw1z2kfuBwB864MvB6Dxlk89DcWGuePq12H+ohX4vUt+CNKEb354OQCNt35ykzmvbZtw+zXHYv6iFTj9Qz8CMeGmS5YB0HjbJ541zPZN+OdrXov5i16ON1/0YwDA1y89CIDGWVc955l/+vRrppzRigvlOePiH4GJC+fFilvOXRMX6scwoZ7btQUTF9r0lHX3gYkL18+dl74Gr/vAt/H9j52IgcXLsObKe8CKcfsfHwaAcdrnf2HS2fIcfnjZGzCwZBne8PG7wMS4c+1vAmCc+oUnAcv86NI1GDhgGV7/iTsBAN99z0oAjDf+7f816fzXc/jJJaegesAynHj17QCAu/5oFQDGydc9bpnn8cCH3ojqActwwjX/DAC4+91HAGCsuf4xU89TxBBToTyGeQ4bPmyY13/yDgDA9845HADjlL/7F8+48/qdNudOmgr18ztXfRcAcOd7X2Xq8ItPAGyZj55kmI9/zzBRPZNl4rYAgDvOeyUAxu/+zc89c//lv4OBJctw0p/fDQCFNnXl+cG61/t2B4B/Pv8VABi/97l/Nec1BuPyctfPyR+7FwBw2wWHAmCc/ldP+XSEmV3MWOql3/LDvzsBb7z9vT2lNx499eZbp8xo+NKXvoQkSfDUU0/h5S9/Od797ndPST4zXeKetA/V8z4Nfh3w6HNSAOfCzGbGu0b0zpiVfwxTDmAuTEkTFUeMPMNjMkSI0jZM4b1F/fVsGWVX/SkEa5fSMUUgxDv5EhOUppCOcxnKQ15Uysu/j1ab8asOuXRy626VKc+4YzHDilvSKeeldLE8ShO0KpVHo+XcVcQoJuiIcedlXGA6tVd8rHyN2T9dr0PHIPjeKwXoHLDB0S5AF0AU+GwYv4JQ4orDHRiUGJQYmPxSFxcSBWY7xgVxxky7dNAuLy6Up5VpLbMJYA5t0Y6BQrF+XOyBq0OY83KB2eTec4kptUU41olp36ZxeynVQ7v3wvR0/Qgza5iSJrqZ20wKhN6wYQNuvfVWfPvb3wYAHH744TjyyCNx9NFH7+OSTT/NnFYVodBD6DhBJMxMZLpuuNY1nd4ZipcCLedFgWlXHtLx5mVR0lFeignKbUzl0NLmZcQEosivv015lCYkEaN0GyZXSMiMebBLu5RXkiuk5MZFGEmuzA+gZZSNPUiSxDPuWMyU02mXFzFBqSSUxxk6Lp1O515Kp20ddrk2ikHvopUcAAAgAElEQVTXnSJi0bLpGGA6wa7v4GJHjIHE/jUTR1cYhZWYHGPf+yVtQfA2n2cAECOJbVi3yo9jErN/QeLP2eZfLk+ZSXRLOm3zajmv0HHqdl4UpdORKbRIax0CIZDaU1E6KOXF8Q2mi+0Fz0RGTMvGdoEJX+7BoaBlI7lS7Ee3dISZfcwc0s0334yDDz7Yv1+9ejWuv/76fVii6SuZadiH6tWCd6NBnTYhEma2M+0f6qQJie1ZtlueEzCdzMR2QnWXpVJdOrrLUqluhEEzt1y7bhYhsWVgza2zCs4gcAzYzxr4dFwHPGLifQnceZvOvDNwOpdHZeQK3XGp1PEySR7t/GvzQk5I3YyFbYvCOlG5AqBbzl1HwzZuwzF3yLVpu6VQ3bXRqU1jhnWXToEK59WyNCls59H1wBFG6V2X1C3Vqv0+Tjp00h1jlwB16TAYOtEt6QCAducVj/6Pl1HFvEL9hY486ejcoaOJrcjQQZhBQWE51WJ5Qv10YBRCOtxlc7eoLTot3Upxe0Vl6MSUrx2Xdy/XhjBzkymrZ4+ISAzMqEDo0dFR7Nmzx79fsmQJHn744X1YoukrMRpmgNwGs6TiHgYVhpKEmTkMMYHdRmkqWiq1QzpFF4IODPXOmA3FDOPy1vZHxDORG4/rpMOVGe0ZYgI7Nx6KOv+RW5Fx03HphFmM2B1IMUHZDriCGemPz8sZBPAdL4UkKzLufbwaUzkd9967KfXIKJd3lJeOXJgA1dY1i9ucu9Io1A9TuDaoUz1HTLc2zd31Zq/D2M3JuSJ5o8OedplhYuSWYTuL4OoTbNyBWDFy30kmP0Nh0mGfThN5YJwrlGPsjsmZX5UJ/tiYjD/39nlppaP6Yd+Rd+loCoaNZxINaPLnrgltzx0o1w+XGLZtEeoHgF2C1tazNWLKbRFfP65tcs90u9+pZ6b7cyww4bzCtUHWnW7sdISZSUw/xEx73T1p8+bN+OxnP4uXvvSleP/739/y+U033YS77roLhxxyCB566CGcc845OPXUUwEAJ598Mt773vdi27ZtWLx4MZ555hnMmzdvr5Z/pkiMhhmk4qh0+7tcmJnFFEeKOzNlF5RyvILboKq96xGHdLiYjkvLr9Nv8yq76JQ3OPO+/yW3It9Z1mjrVlTMy3TKYxedJDed8iRKP8lUIZ20mUBpQqUZ3IrK6aT2M/fXpROfl/t+JWLi1+4zTeyNEGZG2kxMpy8qMzGHMkMHI8Yx9n2ilT+c5Ap5oov1w8V69rMaUVvoNm1qNg0rXz+2o8qtI+AAQufepqcTx7kOro1FiN+XZwgSDa3MP3OUzet4BFzpgrtLcabB1pjdxyGkg1bGbVoWM6p43oEJebHShZrRSW46vEnw43Zpxu+JqFAGTrS970KdMahYP51mI/xfHYyY8mxE1BYot5df7tWXJlo+NarrNmUov/ZuXwAY8cxH5+unXToFUphZx8wk3X777fja176Gr371q7jiiitaPr/uuuuwbt06PPnkkxgaGsKmTZuwcuVK3HbbbVi9ejXe9KY34emnn8b73vc+rFy5Ek888QT+8A//cB+cyfTXzJk/mmFiZjx246X49aP3TDiQKMj+SBbu8/KNLsx0Zdr7pDsm+m65I+9cEAr9mYjhwBSDbMtMOBw2cgvGh3Ij3gBgffrLTCEvGINARYw3EPzwdjSLEKWT5AqqGc4myVTbdDh3HRbyhkTMJLkCuXQ0QryCY5iQ5gqUBSbNFZIsMElG5n0z1KE/FjFt08lVYHIz00GNUIdJB0Z5JtRHzChNUHnxemrfFiEdrzGuQ7Mxmfavzb8xmESDXedam2Naac9opa1RkNnimE66VhwxtpOqXCfdduSjdNgaHqxCOt4YifIqM+3SATGI4rxs/p6x5+5mI2w6iM4dvn7yYnmSqDw2foLczWMnPcxl6OoQhXp2ddi9LahtmY0xkYembpNOod05zovbM67dI8O0ZZ+NDs+6Qv1M02evMBNhxi9NDE2M5x+/C498/SMATCB0v/910mmnnYZ169a1/WzXrl246KKL8I53vANDQ0MAgIMOOginnXYaLrjgAs+9733vw4033oizzz4b//Zv/4Z3vetdk66X2SgxGqZIRIRVZ1+FFx9xUotP4PhlXSTK083CTHumHGg8sXQwJuP8mL0bhnvdhmmbjg+EjmcxAqN02LysEFBNzo/fxhDYQF9FUcfZMn4WwQYMKxsMXGaSTKFiJ0GVDfRtl07q0lEUdfYtY99X4cqj7DEVMeW8lD9WZNqlU8wrzRWqFMrsZz5c/bhYjSQZk3HB0q6e27ZX1Ka9XGMFt5mYoPDCzSw4xnV43fVHikLn2h8030kSl5eyHV7tGWcQuBhwZdMpM5wwVGoKlCQEThichLzc+wLTJp1yXtrOdCBy02LFSGw6Sikw2dkJy2h7Di5GXpEKnfQoHZ3osDKSIh+I7evZ1o+/P21eHDGubVR0u3ZqL+eyVlhhKWLiYOkiU742XDr9e9YVXJgmkY4w+5oZv9wCEi89/GQcddZ6MGivGg0AogUtirrjjjuwY8cOHHvssYXjxx57LB555BE8/vjj/tiOHTvwzne+EzfccANe/OIXT7peZqPEPWlGyDyMtW45JMw0Y0xnus0MQUs6pdkITZ4p70pcnLGIXheYeMYiKk88Y5FHAcOW8a5HuTuvaEYkyiueRdC6vbtUkofRdp3bWYzSCkJJrkB1yzBa3HhImw64qsNDaa7MyLp39TGdf6pZJreuRzqcWKVhO2R1exKsUWm6Hx1zLHWxAw03q6H9sQKTU8iLtckfYfS/0lDIUwZsOqY8NtjUuUu5mAibjtYhhiGuH9IEZKGefbB03BalNvXqdK264/HghYtXcF9ywbdtmMwnzMGwcF+znevMjdprHdyTXNUrDZ0wMjsboZmhEy4wOtW26gLjYxqidMqMTt3wvk0nYQDFvDhh40IU5UU5+bgHzbp9Xlqh6e5lZmjF5t4uMS31k4Q61ImGyhW0H9LX3q3Jt5I1IuK2MMeLo/pQwY3K1A+HGSfbXuDgLtUp1jWe1Si4R8YvVaj/Tkx8wOXFcbuX3KcK5RnrWhVmHzLTS898bQP+/Qv3IdtZw/7p/HF//5FHHgEALF++vHB8xYoVAICHHnoIy5Ytw6233ooHHngAf/VXf4Wjjjpq8gWfpRKjYQo1ebckK9dpbLOKijDTkXFDfZNl7F8VGwlBxOT3Bii6A9lDumgg+DXnS0zsR6/iAGb3NduZ9fsjELVdrUhpQuqCd4FibAJC/ELiRjkROtNOzoUo5ZBOGNk3qjSMK1LFM4RKs1ieJDflqXKYCUkyKnSy0ob5TsWdF5M/FjOsEM6LCZWm8j7zPi8mVJ1hAGPUxCvXuJkJVx4FFOrG1Y9WjNQx7NoiqucObdpW7lqNDYH4JXFYmUuZIOfYh18n1vXIBQwDxTgIIIz0O+OV7IyA0j4znepCR5WBKJ7AMkojDn5moGWFJfe+wKgyY+IVynkVYgiUtsHY4Vy5lE4IwnbpsDcSQv0BKNdPm1WPWLHvyDOptp1948EYt4XreYd0AARPkrjd4ywjVyQQ2qbj9gyJmbYrpvnrp/WjVqaHZ6YwM4MZQ714TjD6t0/Dy84+Di87+zgAwLY/+Nq4v79lyxYAwIIFCwrHh4eHAZgA6mq1irPPPhvveMc7Jlna2S8xGmaCVMsLRAOdwsxwphDU7H7vmYp7JoyRTry8Z8woJr/akQKKrkeO0VQ47uMToo6yigwLQjSLYBkXrOxX/ondb6JVgNJc+ZWRwLEx4NJxcQYuL4VK0840OMYaBOS+a9NREVNpmPeqFuqn0kAhnYqdFSkzcXkqDYBVK2OW/Qx55SkK5XEzFJ5pKgAayMJ5BRcnd142MDoLbWE2CgvlcW1RblPj5gSAQxv7TccU4JdBtYzbq4DjGX0fROs6uPY7zv2G4iDn0LnWqfYzBWQ7/nFgsYlDYLBjFNoHORODU7chGgqByoaxsxkFpjUdYirkxSoyjmx5iMkHdRMhzFj4dBgKGkiK52UUGRZaAW6TODKGRSE42daZdm1B2s8ahHTYzBBEbVHYI4Gjjr0qMQqFNg3XS6mz3+m54W9F7pyXitLu9BwruDl1YKbhs1eYNkyfpKeJ9/vAwAAAgErWr7JuhtVqVVZKGofEaJhCxRb55GYd3HfjYcLyXS7MvmC6b8rWPh3TyVNdGSPrYoPWVY/CSkQhnY6rHvn1NDnMIhRW/gmxCKTZL3vqyubiDFL3uLB7FcTByS5+YYArnqk0E9tJN0ylaWMIcpsXODIIDFOtm870QB5WMUobhDQLzECNkGTAvJrrSDMG/RLbNp0akGSM+d6tiDGwh0sMgzRjfs0Nw2pUa60MgEI61RpboyHklafAQCOUp1qjUjrGQKjaDQ2I2cyY6NABrzQU8kSh4mJHWJvVkzgwSa6QFVayMi5nuX3dVhphhiAKiOVEQ6WulOxnFvyotNJmlL7i0iXTuS65FeVpDgy6dIE8zQvp6DQ3MzODIa88zRG7QuUVG9w84OrNpBuP3GeVzHT2Y6aSIVZeycy9FeXlvufKk1fMbulkGR0bH1GZWREw4M7LnDtpFM6dcgalUf1Yo8UzdsWlJGJYaXDMKNO9itvCGRKhveysSrRJX3lzN29E+NudQx372Qc75+EuHx2twtTpWeeqphQsHQ8IFL7bEizt3um2TCGvafqcn1vM+FXu37iYhumgpUuXAgBGRkYKx937Aw88cK+XaSZLjIa9pHZTeuM1JHqhhdl7TKdN0LxB0TWFIlPeBM2+ANDdA8V/hxCtemQ6tGXXIzeL0PLdeBlUildPKn4e0on2RYiY2I3I7FSMApdkCpVmvAlaqztQkpmR/KrrgGvzPon6hmnDGARuVoAYSBvadsCNqrUcpBlpw3RUFANpIy+sAFWpZTY9x5A/FjOsqJBOpWY7lFFeWVWh6soD2NmI+LyMq1PaDC5MfvYhOvckV0icYQFjmBU8WVxbtB1OLsrv8+D2Bij5yJsVhMJ16GcWHKNM7EHsiqWT3F78ADhyYaKQjk60dzcCA3lijI/czhSYkf3yLIJJI/ezCBReRwyAAtN2xkLpYl7OPcmXOQeR8rMhiqJZjYhRSHyH150XRSPynlGt9erLozQUJ4V+mVZccDlDu7boEK8Qkia0c4XyLEz7d9xIzqVC7Tt2hfiWHgJly6O4nZiOrlCOGTMVYfYW00l9c7veSzr66KMBAM888wyOOOIIf/zpp58GAIlfGKfEaJhC9e/m4uj/TuU5RWGmmikEOXN7ppAOI/rBDh35MtMxL9tpcP3cskFg4ghsOsW+k5dx4wlM+Qfbux5lYUQyjl9wacSBx8TcIRZBAbbjzNoEJ6eRi06aF5c4JWakGVpchioNBtXdMp8alQZDZRzSaWgkmQYaTZsOWjr7SSOHynNQ3RY6Z6T1JiiKYE/rprCB0f5YzLBShXSSRgaOlrRKrDFCNVseDaQNBkduRZWGKT/VXKNqJFnRzSmx9UN1dx26vR9Cmd1St6FN45237bVh21j7YNeoM1vaBM1tcMZAcCuyjFtO1QUee//8mEly5IlGxubcNUxn2rgbWSY17zN7cWiKZhosk1cyMDGa9iJj4jCLEDEAxmRIq2JeldzsK+LLbDr/GTd8/eSlzd3ce5eOW3JVIaRj3J6i+kEUBO03iTNM5paAtUHXOkpHK4biUlt416OQF5iQRYaBaeZimxKTD1z2cRSKC+mQRoj5iGcpSulwu2dUFLhfDJZG5MJEpS+xZzrl5RnEav9cFWaqme7qaTVIBvQ0mWk4+eSTsWTJEjzwwAM4/fTT/fENGzbgmGOOwWGHHbYPSzfzJEbDTFC8oo0w+5gxfwoxBONk4oBlAnXOz/UPQG1HC8OeC4Epz1hQ5P9O1kc+Nma965EOo//tApgVE1I/mlgMTAbcsqNxALMJTo73PEgbdhbB9e+ZkDZQctExswPVpvMlV6jUcqgs+IVUahlUnocOPgNp3S23ZJTW66A8Q2I7+wRGUq+BIusqaZjPVMMx4VjMMKlCOmmjDo5GaNN6HbmuIGmEDSjcDEV4nwNIMOCNM2VnSiKmQWhWorbQqm2wdKFNoTpfP3E/uOAmAuP5ojSY3Wi+3RsgCW3BSeSyZL/jZxFspzBPtDEK/Oi/MQh0tGmdeZ8jT4Nhp9Os4C6VpxlADF1iTMZceB8zeYnJU2M05CWGlQqMNSzyio2RUMawMb11m5edefAbybk4jKjMPsbCTRuoDkHXCuCoMULcg6t6bm2LdjMWiNIBodMqTL6/Tu1nIwhUCIR2sxgtKl0/JoEonSigujDRUO5TuucqemCmzXN+DjN9EAPI+5jHtrsfw7a7H8OCF17oyNTtb4DWRZMoTVNcdtlluPbaa3HJJZdg3rx5eP7553Hrrbfixhtv7FsZ54rEaJgiMTN+9vWPYOnKNXjJqjWTSywOVnMqb/0uzF5l0IExMQTF2YS26ZR+hP3+ChyMCp8FR7NWlmmZxeLQ4feMMyrcCkscuTC5vPymbOSTbxfAnGQEOJchN9MRBydb1ySyAcwKwdWoPIuQ1N0PV3ArIhAYdhahkYMaoTOYNIyR4Jik2YRqNqBsR56hkTRdp94xNZBuQmV1yzCSbBTgkI7K9gCsoZq7fToq21NIR2V7AEqK6TR2g1WlmBdrJA3LsEbSzMBEfgzPzE4QlHVzImakDRcbAVuHZjaGvIsX2xmfYrB0oU0RrajU6VplN7rLPh0XGM3uF8AG6JoOsCm1WxnJBwzDziKUZxrSHHkaVhnS1kgIMw2Z+V7kbpRXmqGMuTUEiH1ANWANCNKBsVHpLUwpHdKqkJdOM7AOI/t5mkHlSWF1pjzNo3swlMP9ZXuuTFE6bvfpKK8QdG3qUCsNJOSD0k0dwy7BatrCuzD5QGgdGRZRe1G8nKtGcGGKmWDEEEed+9i4UCEdsu9b0gFCXtzlOVYKzQqzGkA5MLugcT5XhdmLzCS0+bG7sfnxe/qTWKTFa1Zh8ZpVqL3tS20/v/vuu/H5z38eAPDNb34Tv/mbv4nTTz8dCxcuBABceOGFGBgYwLnnnotVq1Zh48aNuOGGG3DGGWf0vayzXWI0TJGICEeetb5PqTk/6G53tjD9ZNq58dhX9m+3dKjwpz2ifJYdR308Ezrv/qPSRmCFVY8iQ0Rx2EyNKKxy5JiWIGdYv/p4idOmQqWhMOB6mJowUHcBzEbVuhk5n+d89hkYqBk3Hc/UGJVahnmjoVNT3dMAMfvarIw2kNRrqNbsLAAzKqO7QTqM4qe1XaC8jlSbDn7CGqqxA+Ac7sRUfSvATSi906ajQY2tBYYaWwHkUPyCZXJQ49c2F8f8GqBqSzpElUJelM1Dmo348lRGRwqzEZXR3SCej4FR68IERnVPE6TDI3hgj0aeKgxmdgM4DaQN5TvkgDHEym1aiDdBdD35S9Z1JtmXmZV1t/EdVQpBzVZ5miOrZGAfVAxk1aywnGpezZBXmuB5NoZAAVm1AU6bgRmoAUkOGrTnrjQwMBpXMzAwCqVypIM1z6jB3QVGDYyCWSEZqHdNh3VSyIurdWM0ONeMah2cp+DBUOa82kC8AUVWbSBRKdgGeJNiZNXMuDm586pkACWADcxmW2dxf0ynJgaEqraWiYzRQkkoj21f5W9BO+vD4cS0Mqs5+f2riIK7VNSmYPKrwgDxbEUwWtw5t6r4TKPoOVZ2I2pl4uPlZ2O7Z294bWZVwoBE57iH6ft7MbuYiWvpqjVYumoN7vviSXvVPWnNmjVYs6b74OzatWuxdu3avVSi2SsxGmaA3G9AYdKt5MQuzN5iXAch6shbxnXqPcNRB84xtoOXuJWR3HRxKR3FYc8DaN2ymRpgYw3c8pp2ZSSKVjRy7kvKjkprHa2e5Dq8dhnUxG18ZjcmMwHMhkkyhSQPe0IwM5IMhc3U0szOIth4Bc7NSHraCJ30tJEjaWZQmfUT1xppow6zU57Nq15Dko2Cmnssk0Nle0B5zbtmqOYuQNeAxi57XhrIdwKcw5sf+U6AG+DmdsNwE8i3gSOG823ms2yHZXJw/oJtX8e8AKIE3NwapfMCQNVCXsQ50Njty6Oau8CURGUeQUIK1BgN9VOvWxemqH4yBbsvGXLNNiak2BYtbaqL15iT24VY21V0WpZTTRjOZmHrax8bBGb1JA2tnCGhratPWBkpTzPoSgPazhxo1mAXue4u1bQJlTahbOS61hqVSt0EFVumWq2BSINtOsyGsY0BAOG9jThn1qhWa4WViCqVOnKdFPJSlTq0TkJ5Kg1zrq7MOjczFvGKRnYGw7k5aXYzH0nEmNmJuH50kheMjzzRAFPYlE1rHwsRYgM0OKFCW7hlYQsMK7+9i9Y6tGkUl0DkYmlMmX3cQ8y4GRXXXi5+orAyEvyu2jr6bmGlJsAvN6vj+IR49SQVVpcqPleLndbp95yfu0y/NF1iGkT91axv1WaziY9//OM49NBDMTw8jFe/+tX4p3/6pxbupptuwrnnnotrrrkGb33rW3HnnXdOiJkKORc9FQ9dl1arEGbvMoWjpWE7zxTWLp8Aw4Hxo8nR3g3Gp9gaKva4ilwRFJtZA2X94hXIdzoRGT2ko4FCZTdpKwUwVxqE1LoeJTBuR2nkepQ2jGtRWte+PGkjR9LI/Ghk0rAxBrW6LQ+gmqPGSLDZG4NhJ8i6DBE0qLnTdMwdle80Hfd8p60f09ln/YJnWL8A1iNg7dyNANYjAIe8wKMAjwCeYXusyLAeBetaIZ1yXtA7gXyXLY8G8t2gbMQzlI0gyUahmvbciZE0a0iaUf00M1QabIOmTR26WJG4LeJ9NRSi1a7IxUJExgIAshuKceSe5FyRtPOjV8ZNJncBzAjxCpldukqT6UybzrPNI82AtAnttvBWANImkkrD3yRJpYGBSg2V1M4eJUA1baCSND1TSZqoJk0kynT2lSJUkyYG07pnBtN6CzOQ1lFNQ17VtDWvSqVuDA6XV6UOVakD6Wgoc6Vu/Odcw1ca0GmGLAlLYuk0s7EXtn7SDHmaoxnXj10pytePDQrPKV7NSdtVnhyjW9rCGwTRyldMQO7aVMHvxu0Zux+G21naGY3mkRBeM7Hb3N1cGx1WSQrPqOi5NCmmh+eqMPucEYm6adbPNFx88cVQSuHqq6/Gf/7nf+Kaa67BmWeeiXvvvRcnnngiAOC6667DunXr8OSTT2JoaAibNm3CypUrcdttt2H16tU9M2X1b58GozjordRtFWaSTNHFY+x0Oh2P0+llkQkXwBxfH35/hSjIGUDLikaFwGcOMwvOkd6tghR83uF3bnYLnLhVj9zSqC4WIckCY96HZU9NLILZ48Bm5WMRVMOtEmNWK0qaYTWitFFHko2ikoVgY9XcBXBwifEGgu2AgxnQO8F61DOcvwBwA6Qj9xUursENHgVzE8RuAwcGc83mZUrN9vvE9cC4NGOGklI6u0GoFPPKc5B1YQIzkG+D90kBAL0TlFWhMnceGpTtQVIPj+CkXkPaqPrlXYnhY0IKbZFHAdKEgttaO5FbLSlaBtWtlBQ2LIs6vP59bjvGpuGJ2LgdxeuFWgMhsaP2ihjVas0aBCazwUoNqcqQK5eOxmBSD9UMYDCpQ5FGnjijqpWpKpNHzLhjcToZJ8W80joynfryVNMGMpUjT0OZk0rDzFhYRlXqZjTWMkQAVxpgCu2uU7P/g5uVUBRWbwqMBnEegqaJbLB0qELTFlRqCw4hKAy/CpNf3tX4pSFaNCssoxu7JMUGAcN/J14i15SLPeNmHgJj0ynFRsTyIVolw7TgZlT4DexwzXL4bj/SEWYSzDjUbp8GmWmYnZrVRsNzzz2HefPmYf36EFtwwgkn4DWveQ2uv/56nHjiidi1axcuuuginHfeeRgaGgIAHHTQQTjttNNwwQUX4OGHH+6JGUuT2afBTSnGK790mnYUZvIMemQ6uYUo+8Pu112PpuydC5FzTQ7e/ChM2RsXJvujyRzFK4R0FBMS9yML9gaB39xNl4wPbX3iNfmlFt1ItjMsWJulPysNeMYFMFec6xGHPQ9c+dNGjrReR2rdbwiMpNmAao5GLjqjoGw3kNsOOGdAvhvgRsHVh/UImJ1BoO0MQcTwKJhrYLh0crCuwYyxOjeQ3QByMNUDwyYQ2jHMpmw6Toet+1DEAElLOsyNQl5EGoxaVOZRAKMlN6eqmZEwXwJlO5Goik8nyUaRNIagmm4pWfaGWtwWJhDaXieag1FYcvmgqPPoV0+KNwtLtN85mQG7ilBwmzE7PWeBIVMApfJwjVXqGKzUoGznGqz9rIEb1qwmTQyoBpCE8xpwI/h2Om0gaSAhDR1OGIP2mMtrKK0hZ1Vg2qVTYVXIq6qaSCn35RlM6mhQxQd4gxnVtG46OZappA00AaRpWE6V0qZpTVfPaROaCbCb35mZhtwYVYXVnCj80hLblaXCfepcmFw6pi2KjFuFybn6MIcA6rhNmRgqid77zeBsG9q4B/LPqLDHhjt310kPqwo7w4PDs84GS5MPqO7l2ZuH4Ohp8JwXpj0zlnrqtzDMvSGadZrVRsPmzZtx8cUXF44dc8wxWLx4MbZvN77Od9xxB3bs2IFjjz22wB177LH4xje+gccffxxPPPHEmMzhhx/ekn+/9mnwK4iVp4CjToIwU8+QJrALEO6aDvnjhaDUqE/nXY8QZhZ8cs69xB9A28DA2M3IuLKolqVS01whtSM+RIROm7Kl0axGpUG2s2qOmY3JNCrO9Qhk3I6ic3axCiq3HTptjASzGpEdvc32WJci53oEQL8QLTUJayCMAMq6DCltjARdi5jdYDRCR57YGwCesbMKmt3mEs6QCDKf5aGzT4Dm3fCO2o5hQLuZBoIpCwWG2e0bEHab41I6zDUgfwHszp00kL8AygZ9/VC2G2l9ASr1sP9E2nAjrlFb5NGu0QxK+XAAACAASURBVEr5JXGdFJswR98ZtKPYrSsjsT8NUmFmwV2ruXU9oiSMtlcqdWM0WKZSMS5DHLkMDSZ1VFTm74sB1cBg0kDmZggUGYMAYeZjKDF1mFFIxx1zGlQN5KxamBzKpzOYNKBZteSVswrlSRpQpAFlDMREGcPGnySzcYsCwNYoSZSJscgolFlV6tBM/tyJyK/eFGYIcug0Ws2J7OZ08bMgyaG0gtvt2bSFju5vE7ROSvkhfaXCHhCFNmXArVBANpYi3sPDLLkaiQhMNtnoOWY2m3NGRLRJXOlZ52Ky2Lo5xStQ+fNztpmNum/dN4IiZuznszBTy4ylXvZpYPQ3puGFex7BC/f+DIu6LLkq2jsi5nFcLbNEQ0NDuOKKK3DxxRfj0ksvxfr16/Hwww8Xdgb89re/jbe85S24/vrr8dRTT43JvOtd7yrkQUSoDi/pS3mZNZq7twOkUB1aBABojJjATJeHMPuKITRGtljmABAArXM09xhmYGgxAKBumYHhA0w6OkcjZhio77bMkGU4YuYvBgGoWWbQMVqjProNIIXB+YsBJtT2/Jdh5r/IDA5yjlrNMPMGzXmMjhpm3jzDaNao1beCYBhiYHfdMEMDLzI/8Kyxp7kNBIWhijmvkaYpz3DlAAAamjX2ZDsssz/AGiM20Hg43R8AGybfBQJhKF1g0sl2WmZhxOy2zLBldllmARi2PPkeEAjz0/kAgN2Z6ZQPpUO+vfbko3uBoUJ52pUZIHteI9G5q2L9UGLrcDsICvOrSwAqtoVWUZtCYXD+EjBx4fqpj25DZf5CNEd2AKRQWbA/QIzmC2awpLLfItMx1Br5zp2AUkgXLjT+7tvND3OyaD9jAGkNvWMEUIR00TAIQHOrOa/KkgWmx5lrNLfvARShsngIBEZji6mf6gFDpsuaazS2jwKKMLDY1GHdMgOW0REzuHg+AEZtizHYBg+Yb+pwjHTQJi8GCuUBCByXeZH5XnPriD2vYcNojWz7bn/uAJDZc0+XLAAzFeon2d8s8Zhvs3W4eD9jPGiGfmGnYfbbzzDbd9h63t8O4sdtYZhsh2mvdP9FNh2NfNcLpk0tE7cpmMCcI3PMgv0Ns9ME/VcWmvuWtUYWXxtjMcPmedjcZZ6HlQVLTJvqHM3dhqlapmGZ6oIlPh3/zBRm2jP9ULJ8GK/835f3LT0n9fbPY+PGjX1PV9S7ZvVMQztt2LABQ0NDOOeccwAAW7aYH9oFCxYUuOFh8wOxefPmnpiySCU45RMPdC1LrzMRI1uexb1XnojBxQfhlCvuBQB8508OBaBx6voHhZkAQ0y49cJDAGj87lUPGWbrs7jnyhMwuHgZ3nj5vSBNuOVPDwagcdqfb/TM3R9fjXmLl+GNH70PxISbP/AbADR+78qNICaMbH0W31t/HOYtWobTLrkfAPCtD74cgMbpV/wUigkjWzfhjqtfh/mLluP0i38IYsJNlywDoHHmR3+GJFcY2b4J3/nMMRjafzn+2/sfAGnC33/sJQA03nrR/0WlmWBkxybc9PkjMbRwOc4676dIc4UvfmoxAI3/733/impNYefOTfjKV16FBcMr8J63PwbSwGf+10IAGn/8jv/A4B5gx8gm/PXNr8TCoRW48E2PolLLcNnNSwBofPi0X6C6Zw+2jTyDq398LBYNLMclx90F0k1c9P2VADTWHXc/qLEVW0aewyceexsWVV+CS4++EZy/gIsePBMAY91v3wDWO7F1z6/wyf97ERZVXoxLjvgEGDku/j/nAmBcfuRnwHo3to3+F65+8iosShfjQ6suA3MNH37kwwCAjx5+KZgb2F7bjmt+fi32T/bHB1deCCDHR3725wCAS1d+AGCN7bUd+NQv/hr7J/vholeZ5fYuffSTAICPvOoCgHNsr+3Ap//1b7F/shAX/ea5ACVFBjDML76A/ZOF+OCrzgdIFfIiGsT20e245uefxv7p/vjQqksAJPjww39myrxqHZRaiK17fo2rn7gUiypLcMmRnwGphbjowbfY+vkSaOBgbBl5Fp94+DQsqh6ID67ZiGygWmiLXy9LQpsuWIH/vvZhbFsygm9+eDkAjTdf9ghuvuoIvPoj38BPLjkFA4sPxDGf/xbySo6f/MFxAGu8+ku3Y89+I6j96lf4+dv/B5KlS/HKb3wJ2fzdePI1fwhoxqHfux7zhrej9vwWPPGGDyM9cDFee+/HUFEZ7jvk/YAGjn/oz7Eg3YORZ7fj3td9CoPL9sPv/vAiVFWGm5atAzTw+z/7IOarOkae3YF/fPXfYGjZQpy1YS1S0vjiyz4FaOCdj/+xub+e34l/+K0vYnjZAvy/D74HAPD5l34G0MC7/+U85KwKzNsfeC8yVrjuwKsBDfyPx85HDlXI6y0/uQA5K3xt2ZWABv7wZ3+GOqcYeXY77nzttZi3bD+84X5zbf3zyz8MaOCUhz+KPdkgRp7bio2r/xzVgxbhtfeuQ8YpfvKK9wGaccwDV2NPcxC157fg5yddjOTAJTjstr8Es8ITv/XfAc14xb1fANfnofarX+E/37wWyUuX4hX/+39BZSn+5cTTAc145Xe+hrReLbTFEV/+Gkgr/J/fewPAGkfeeDPSZoo9v/oVHjv3D1Bd+lL81hdvgsqTQpsmzQS7N/8KD//Jm1B90Uvx2r+8FcTAD95+DMAar/vC96DyBCO/fh4b/+w0VF/0Mhz/6dtAucI97z4SYI0T/vL7UJow8l/PYcOH34iBxQfixKvvADHhu+ceDrDGGz57P0gTRrY8hx9dugaDSw7C66+6E2DCnee9yjCf+pF5Hm55Dj+87A0YXHIQ3nDlXQCA29/3SoA11lz9E/t8fg73X/47GFx8EE762N2mLc4/DIDGyZ/Y4JkfrHs9BhcfhDXr7gEA3HbBKwBonLL+Ac98/2MnCjMJZiz10m+564trwBLTMCs154yG9evX42//9m+xeLEZTRkYGACA4moCgF/rulqt9sSUxcz44TW/798vP/4sLD/+rIkV2q3TH08sl9xdhOmd6eQy5FcmAoK7UAsTvSxtuNaOKaTtXJIKAczRMqiWad2UzW2m1so4dynn5mQCmBWYtQ2ajTYLY+cTDxAUGNqugqSjTcfC0qhh8zLrepQFX3KV7QF0w1cPZbsBPQLy8QFsA5ZHQhXyiI0RcOnk1h0oVDNzzcYruH0ZGJp3m6BmEBhsXYaa3vWIKApqtgzrGpibyPWuiKm1MtDIbRyDcU+qAxwYrfeAoKDdCktkXJaIVSEdKEAjxHM4FybPcA3MVQBuF2s2rlg2coUBsN4NynaCchfUnUM1G0iUKrRFkqeFNk2yaGUtttclwW8WRgrBPcm5zajiTs5mVZ8cVAmB61RpYDBthEBfBcxPaybOwAU5Jw3MSxohqBjAfNUwAcqWma/qmK8ayFxQMYAhVUdCGkSmfuarOlLkYOuWloAxRK59DTNENWSUoKlqUV71ArMwGUWdK4W8FiSjqEeB0PNUHQlr5NYlKQFjOKkFFyYw5qkGdKKQJyHIeSitoZaHZ/2AakCnZPekMEzFLe9qlaYNNPMUqLhYAA12sRDO9chufFdsC21iUmx76SQ3bo+WYZjYCNfWrk1ZEeD39TB+R213jQ5F9Eu1hoa3wdKO8Stvlb4DQlh/sYvDggpIxz0Yov5lcIUqJRsz3hWqE9Ptt0CYrswE9csf3Yhf/vjr5s1CQLf9IZyckrER0RRrThkNf/3Xf401a9bgzDPP9MeWLl0KABgZKa6+4t4feOCBGB0dHZMpi4hw4gdvKR6c9I0ZJdDRq0yYqWKICdRuc5wokMzsCB0YPypjGcWlTdmitOO9HBQTkiQwIRA6MGmuUKHU5mU7j5qsnzN8vELF9x6NTzxpIATfGiNhwG/CxlCZRtIMgb5JswHK9oB8J10Dmd0HwdWDDWCG6zhTiEXwwcDaGgSuMwi2wck6YnZDR0YDI7cd7mZg7HuQY9h09u1rwHX+myC71KWGBnMd8T4N3hBgG+zK2hsWcZlBSTEdXbN7MIS8SGsAoQzMdQCqVOYK4OIwYIwPxdWIGTWGF0XB5FkxMNvsYzGv0KbeoHTXmF+G17WoWbufFftgGlZmZ2R2G5OR2Yegkjb8NZ+mDQyoBliFoOJBZYOTbTqDqoH5qg7tA4815qk6hlTdMwvVKAaoCe3iHlgHI8EyQ2SMEdeRZ8uYNjCMNxCSEJTuDIs4nRR5Ia8hqiFVFV+eIVVHytqvCgU2RkvOiWfm2++71ZOUPXfHA/CxEtqu1KRgYiN0FCheTRvIdYI0DUYwuZWaXHBy2oTOE6DqBiGsQRD5mutEQ+cMVPzubmYZ3VKbsgIoUSGdNvs9QBPSeEyEbDrRPg1g8ss4O5X3aWDFfrVOhuvsh/pxKy75ZYfBbdMxz1UUmLieyyqUqiPTSzrCdGM6qVNMw7LVZ2HZajMwes8XTh5XmqKZozljNNx8883YuXMnLrnkksLxo48+GgDwzDPP4IgjjvDHn376aQDAUUcdhf2sD2o3ZmplbtLudrswnRizDCr6lFdxpYm2O/H2shqFX4WkuJSqT0sTKHOFVmE2wkq5/RWyMKoY9lewnb6ckGaAypwRASQZQ2UcmIY2swpu/wAmJM2m34ANAFTWMHsONN2KRgzoETA3fDpmxaMawKHTbWYMQjqMhjUSXOcapvMfB0JzwwYWu46zYYqB0LXSMS587hhAg61BQGDoMsNNY7CwC2TVYF1HHOTO3ARzE1rbDqs1PhBNvbuZBuZ4dsTMRhTKjApYu/rQADcBDvUDbtg9IEY8Q3kNKiqPyutIM9OOgGnTNO9wxboAfii/xKqTVhqc5GAVOs5IMrPngVU1baCiMjQpBLdXVYYB9x6mcz1fNZD5YGDGkKpjgMIO3gPUxHxVDwHMZDvy0ZbiQ9YgqCMrMLGGqIYcSZFRdWRRPQ9QEwnpQl7z7cxCKLN5n1EwggcoAwrnVTdB1854ZUZVhc8BYzA1VIrMfo+1WalJR0Pi1aSJLE18QDXAZvYhWpYVaRM6q4RlWdnMPrh9FgCYmYg8rBylYTeAixmloVRYlpXdpn5Rmd2SrO5a8EzUEXTB0i7ckbvMEPgV4jo8NU1+IaB6THVlig/xljJHm+p1lzB7Q+KeNDs1J1r1u9/9Lh599NEWg+HRRx/FySefjCVLluCBB4r+fBs2bMAxxxyDww47rCemnTRx13/jlY6/Qu2bTpgxGN2eiZcm7cQA8M/cXhiO91uINt5S0Y6+rOGXYS0zzgUFWoeVkSyT5MosnepdmML+Cu62dvsrRP03YyRk2s+GJHZWwS3zCa2RNGpmZsHNmGR7AL0bnO+2WbllUEdDZ4FHwTwaOubarmika55hXQPrGrTfq8DNItQDw3WwriHP7I7QcK5Gdd8tMTMGdeS5W+KUwboO1hGj69C6hszuCcGdGK4hy+ue0VyzLkmuOvaAdR155hhA61pLOqbMo1F5amYGwqVjy6x1mDnQXDPGhaW0c82y5dGagXy32RHbtoXK9piVrZrhYg/uSba9bOcus88YhtuTgf30Ayuz50KehBWokNqlUi1TTZoYdK5GAEhpzFcNzFN1z8xTdQxRDfPt7IhiWHekur8t5ltmXjSDMmSPOWaIatiP9mCYzDWWABimOoYppDNM5juDLi/7vaEor6F2eVGtJS9T5rDJoCszojLPU3UM2PpJCN7Nyp37fLsqVMXOjihlVmYaiJgB1TAb0tnOPiljjKk01LNKm6Yt/CyUmfUpbxLHipFbA5usW1O5TbXS3tABAToxRoE/eTtDkLnOPsFsEhcx5jVDu1WYCGFvh4jh6HeMERkepechl/52ZdyzU0V5KS6y3Z69pXSEmTjTScrOhHf7xyCzSECf/4n2vWb9TMP999+PK6+8Eu9///vxj//4jwBM5+DHP/4xTj31VBxxxBG47LLLcO211+KSSy7BvHnz8Pzzz+PWW2/FjTfeCABI03RMZkrlBqXje6Y8pShMd6bd71knxnsJhGn9MItAnvGd/WhJVr9XAuD9hFv2aYim7JV1L0ly5Rm3eZdbKhVkXJHiqf8kI1SaCgPadRTD/gruRCoNoFLTGKhH8Qq1JlQeXHQqo3Ukjd2ghluO9P9n793DLb3qOs/Pb613X8696lRVKiE3QgIRQkEiidyJiBCufeGZngd5RmecdmiwpSO2OkMjGgcvbTvCaCM6OrRCO9jo0z0OtI22xAsKhE4UEMgTMEhMSCWVqjpVp85t7/2+a635Y13e9e69z8lJcsokVe+3nvPU3u/+7t9a72W/7++3fjeLKteR8lS9iliewpkVrFkJx8Dg3Ib/Px46uxH+6lVyG7elcKAzWLdJKnEKWLPmjZDAMXYtGBcxFMRi7TrOlmksazb96r8ZZWNFIyNssZ4TDR1xgZPN2cYGciaGKTmcacpxZtOHIgVlX5ytt2VjiSuxJnpZ/FiScZzdxKKSHJzDhnyL+visodQc1oU8DByuOu7r20dJ1ZmJc5p6N2TXYb4TzvnwlhSOQ+gf0Bkh0bMg0OttMqMHSc6MHrCv2KTQIc9AYL9eb/ROWNRbLKpNnIrKtWVJbTIng/TzOqDW6ElJFUuc4liSDQpsmuOSbFJgGCWjwXPy/YjvrfhrrMCxJJsNzgFZYyCd5lhqk6HrpPksqU0GroNV0fhwLKs1KnSSs6j8fArtx1Ri2a832JRe2vcF7b9f5HkYsQRs4MyG3hJVzMNwzjet69ZznulusWk10q27atvOCKxOJNMpEauoQs6JE8F0TGP1P/aEiCv7xK7fUq/A+z4cksLdrFCXao1eQ2WxWmqDQIKsLPg95UaM3UO9gu/S6/yeJcHQGPcINOQ0ut6N3+flEXBoOY+W06LFDjinjYYvfOELvO51r2NtbY1Pf/rTjc8uv/xyfuEXfgGAm2++mV6vx1ve8haOHDnCHXfcwYc//GFe//rXJ/5uOGcNQXds3KDHPbEtJ3HyBLlG0jPjKyxjcpKuNT1ZOo+9bVQ730YO28mxgkRPQ9yFvOfCeJKzA7GTHF0JlHGFvvYsRI7v4mxTCJNYhy5LxNk0HWUqHwITVrexDqnWQ1JzUFVj6FFcwXSW2CgtynGx4VoWohMTmhMHH1Jkk+JsfKhPfnhcGbwEWQiTLcOKvH9v3QBsmRRwF5q75fDeiSrJsXhPQ4PjKpwzGJt5GsY41g4R0Q05ftu08KQ8yXmIE1XP2Q4Q6aTcC5/PUYIrx/Y9P4Y+LwS7mp2L1clzWoWKJkFQzKNxsREYro5/z+LWO51hyiEA4z0LekSe5NyVkiJ6GrD0pPShR+Gyjiv9w2yVfE4G9KRMuf89KZlnkHFCeFLs00CdizBQZYMTBk6cCs1MDK3BTeUUzqSxBMc8AwoxaT5x7DKEMAnO7xdlg2NE0ZU6fGtO/HmOY82qIUNX1D0esD4pfCyhutQFJh5D8d6I8RCmUTGiSOFP1udSmOymoiuc6mQdoW3wNOjGOXWhf0OYsn/vakXeKuvD3aIcmzX+i2cwGgRJTm2EJMSOzSnJ2W0bwpTuz7FZ/LSwosix2baJ+3y0grNtGcfPOTYi2Ol54VrOTpy9gIM2POncxDltNFx77bWcOXNmV9y3vvWtvPWtb33MnLODuHItE9taTpNTP7S240xum+RMtQDCW7WjHJ8I3UxgzjkqeiNCkrMgKTk6N3SUE4qYCC1SV1TKxtEVdMONWUGKdY/PAF059KiiO6obS6lqhNiynk+5BWYDiQqUAG4QPAZeUuzIrJLig1/Fz+LxY6iNkzGjwQ4yOV4hlrQSGuXUOQ3WejlkCl00EmpSzDGIq+SSDIBaSR+Cswixsg0TOQ3WDhqcOJ8cMe8hlxMNglrOMOiSsZJPNAgyOXaIlU6S40I1JyuddHysG6AYIUlJD+PbrewYbk2cU125RqO/eJ3oELaixTf5MoWJJw9TeCPBxcpIIhTK1Eoy0JWKOTVM1ZMEnzDcl9rQieFAZfBYqKCAz2fepHl8ONAo7JcWmKdpLPYpKTDMSKwgJPQp0dQG7hxDDIpRylfwsqusAlUf3/05jqXEz0dTXxvR0OjFsCIc89I0FudliFGaKlZYCnLysXyuxsgf1zCflFAdMBvzHmJ4kvjciNJW6bZRSMgl6dTXmNKVX+mPXoKixBYd6IaFESU4bSbOqao00qk9lCkROmum5hSoItzHYsO4rFJTfF+Em5xGZSFCwSBVDmddSpYWJSlkKXFChSUVr8uGDtm8r6ZkaSXbch7ps6A2Yp4Yz6YnD2dv0IYTnZs4p42Gcw2N1fZtXIotZzuOa3CmhnNEjt1BTlzVs6RwpLgtVUqytvHVaWFOrqo5dfUbz1fx/aiej66kUR+7MIKyYGMitKs7N8cZ+yTnEsqgXFufSIutKxGJGXpPQumNa+sM2HUgq4wUEpiN2QpiDONVjwir5qn6kLM+hMgN6/CbsVwE62xS3ONIzvr52ZCI7SsaDZvhSXaIcyXGRE6d5JzkuBJny5TTYI3BWW8ApFMZvCspX8GatC3nOFFNOa4kTwfzidLUHFeXXG3Mx5VUZmy/pJOOT8r5iMcH79ERUY1zocsSRtEzBKpyE9cPgAnGryWEqei6qo/VhkJVlCk52tILSc0ptCZURiqzalcxgTnOpiclfUo6tjbyopEQMScDr8zb2hPTC0ZClDOHl7tFLKdqmJPa2wTQDwZpzunhPQS5nCGdxljzDBjQqeXgPSWDWA0seDWq7Jz2ZMQcilFKyLehIlQ9n+iNqKhD13pSUVGHb3Wl9MnjWcJ5X40wug5h6usRI9tpnIuiKDG2lqOLEqOrlPdgrT+fMnZOrTZUrvbo+OpK1PetELJkYsUwLFY7xLomRxyl1B6d6WFFWd5MnmiQ3TOdcpjMMxQrKo3fV+sKYtPvvclT5vdqKieH3cXzouXszHkiYe1P/5K1P/0rDrYdoR93tEbDkwJB0Wz8uMd/6C1nZ47aPaexrX7n+ytEjp3OcbWnwZG765tyVOQ4l3kRmqVSi7BiKbbuwRDlxFCkjk2xWaEykk1yVGXH+iv4ajzYuuqRmK3gSahj/53dCsp/rczakLAMXnzscRARPQ116VOblOB0lOwgJBFHI2ZayJD3IuQ5DdFoqOUMwwr8qJ7zNGU/4zA1PKkMyntW9WhqeJJqyLF20PAmOTP0mbuxUpOzKaypMWfp1B4T65Ih0ZzPsBnm5Aa4LLbd2YH3Fpm610UxirHtdTleP0QyWbDaV0tKcrShq0rKrIxoX4ck5/C1mVAFqa/q5Nt+DE8KnJhUvBXLoQYFvE/N6VMyJ8PEURjmZYjONPA5GVFg2QjeI7GOudTTInAYUaHohepJYh1zMsYJctJYzhsf2mVjMaTAMpDmnCtU0p9iKFR3PFxKZWOpIRW1N0LEpapLOWfkOqkErDjnvRMZp69GDFS3cS46ukSrqjY+iiFbug8p5AwfsqSzayMkRrtQT9USwpFyAzckS0dniHMQcx+S/a8sztZxWJYsyTkzLIKvMszH1cnSZBwnpEpjqdBbzUlGRB1jSo3p9/ldcWg5j56zN9hLT8P8jdczf+P1zH7v/7FnMls8OrRGw1mCc46//u13cfjZ38HhI6/YE5miVP6GvFRly9kdRyk9lRMTlAF0XoUmcMZjdbVEhT7jxMpIsU5+rD60EyeEVCkXm7IZ/zprACcq9mAIbndHKp+aSm8iqMqFpmwKh2/Qpsww5SuIAswGuBGx6RjWvycqryLEUqmxAZwjfB6VYgFHDNsJcqIC7KICJUEBrpuyRSU5lUFVQiy5Gs2hVE41n48tg5zAcaUPY4oKvkTZdSiLD08yWFs3zPJyMk6qUlQ2OP54xvmUIKopxwxBdRrzwVAbLkIW1lTLcareL1ES+j106uNjvcGQjAbBv5dO41woM6QwtRxVGW94xuvHxvCQKEb5Ov/ahOvQv+6oKsXRKy10xK+IixKccfRU5UOPgtdAU+cDKPE6YJ8y/cXrORoJdX7AkB4VPerQo2gQKAED9KnQODqBoxSJHzk9DD0Ma3Es5b83LscbH2Ump6IQm+YzL94w6o3NuUI19sswYBC8GlriNp3k9CkZiv/Lj8+QTn0MpfTehlBOVSvoqpKuK8KkfbjS+Lno6pLKar8DxlIog4pVlwCl8SFLpmic00aTOCXew5QuXG8waNFJd/TVk6y/dyrxxqzyv6nIUUpIpVsjJ4U0hftPuB86cX7O1iSPRTQ4RGf30MCJ30kckaly8jCjPJ9o22dBy3nMnEeDB790K8e+civ+5tXmNJyLaI2GswQR4Tnf9dN7Iivpwo3Fbdty8Iq9jvGZu5DjrKvDkbbhWJd5CMZct6lhlssbtzU5KVzUMTWEKTdQoE5qjl4QXfl+C7EsqzhJnZwjx1dGMnRibLuDzqBCD+tmano4QKoNJFYEMg7MqldE42qgWcV3U45JxRXOrfvKR7HCit3wXZklVhbKKxG5xHF2M3kjXAg9snYz42z6ykOx0pGzWLNJHp5kzVb4Xp0sHasl1ZxNr4Tb2JTNYavxqkexl0OW5Gw2vWchkwNkHFtvyzgiqinHDpE8XKraRFQHa0cNOYJujiWqLrlqrX+dN4mzm6HqUywla3F2w0eXZOdCVxve+AvntBgN0Wa2vn5MyHPp1Cu+tqh8cm285ouSOT1AsqZsC3rgKyEFOYtqkyW1iVN1P4NY5She8kuywZJsUmYr+0uyxRyjmsOAHpZhVIqtYw7jnTOBs4ihAEbBplfWsUQzzC++H6qas4hpcBYxDLHMxnKhFpYYMcy9CPgGdWXMe7B+vyqnsrE2KbJ+D+JixafmvvvjO0jHcE4N6bmycQwrNOg6hGlR1RW0wFdhGtmiPhdYZvSAUop0vvp6yKjTSTkoYOl0hr7bdHZOTVHgYvdpLKZjUMalsWK/h9ih2gUvlMpCoXxPD0leDeuoPQsxdDOUblVZ0nXyGtht7r2Outvz7aSrhgAAIABJREFUNhzyXIQJTn2fTXK2e15Ex4i4xIlyxzn+9c5yzifOY8GFR17BhUdewR/96ivrVYsW5xRao+FJgHjv3KmccsthLImuRqMS0Tac7eXUd1flhFiyrlGWNYb6RGMglWV1yWuQz0XZLDlQVOq5kMaJ4UmpPJ4Ez0I9Kx2atqlRUI4QlDEoW6X5KFsF5TKUllSx6lHeUGwLxyglMNdJzllYUerIHFejoocgD60ZhiTmOqE65ifUBzh6GTKPhRvLM4iKvqtXpWN+QuKERGjX8GqUjbjq6GmovSMueRpyOf5FNp98vpEjqinH+I7Q+ViqIUe8QSPT9qtugha9Dw1OfpzFBY9O81yMn1NVjVB2LnH8tVSHgAg+3l2KmADvkKL0lX5Ufb66qpxoytan3qbE+ffUnDm8F6GbeRG8V6He9x6WORxz8bpUwuyYQtFzgkZSuJ0olZL8I+L7+L8oRW+KnEIc81KP1cNRUK+kzoXX3XgdBm9ELxM1JyMqp1MjOx29EdJMhO4R+jIEOXMyoMqujTk1ZOg62GB8aBW+J53sOFfM6BFWx7wQ6kZy4fbTUZWv0pQSs6HQJdXYOXWql0quCpmnIcBq6/syBAVcghdhvPGfEoXEMMfUebqWk7wIMSQudZWWjBN48TqMxsO48m8k5ZE1OhWPId2f41rNFCPiiftsevJwWrTYCa3R8KRAWM3Jb6gynkzWchrJeBOcsNpkH54TV9N8WJHn1N6JOFa2LZMjVqCqV99SknOMAI7lVGOSs7FZToMkGZInORsXyqe6xNGVQxkDpi6DqssRYgYp3EXMANwWhFV4aw240USoT55860KH5MlY+6wPQSpxaqdwQoiOM1Nj9vOGa9bayeRk6w0EExKhbcg5yHMNbMhxaCQwTzMInMmasrl6bvlYzlCV2b5vYzTkcvy2sTmLqhOzrZ00LMwQkU6SY52bagw1j2GowpRV7PFGVvOcihmE6cTrx389NgKzDpw2vsFYEKS07+w8pFZUe1L5ykgx/EZKejKiQ51f0qOc4MxRsu5i9STLHBU9XF1yFUfXKbqZZ64bckKiAthFU6BSSBLO0aNAh3AggFl6mBCiFDmzdKiwDTnG2cZYPScYyecDYJhJOSguhTnVoUcVw2AUeDmmNpaC6PmQdN2JhmAoSVs4m8nxRteA2gs1GwyJCG9EVAxd7dXoSIXSNt1b+nrElukzInafDiFLY+fUaIOR2MvBN4DL9UOnrE+EtnWna6tc8KAGQSGyxITfk3USwpHq+2E0CKrs/uzE+cpz8X4YKinF3JpmOdXs3pt1jW426NzuPr+bZ0HLefScPYCDtnrSuYnWaHgyIN1Jsx/h+G+85TwMJ1LqkJ/ISQ3XxsqjTpUTv7MDR6beK5sei1jOMFL9Qp9Lr5WlTrUVH4ohNls9tg5VVagqGg0OrA1hM5FTghvVq/8EpTOrWuNfG2KegZflOeSckHtQw5KXSnXOgLNZGXAJhkOVcarAiScjliYdNz4scdWeqKTnYzuDzzPJxh9T9n0+Q5nyFXxOhWnIqZOtQ4hLvi3jONFjcmyzSKELlZmyhGqHbZS196FchvpC3MHTkObjj49Ip3G+xs+pZOM25569UQat6n3XyqDF0InhQM7UPRiCmOhp6LpYmlToS0mBTTpGgaXA0Ykr1c43Xetlvy/vRVDB9PFzLlDoLJm8Rwclmk4ouaoQOtJt7I8SLyXn6PAvlzOkbIzVRTNyJs1HO0GLoxNX2/HekIp6zhoXPCa1MdGXisrpjGNDzkJ9jReYpHsDaDEUYtEpZt9zxj06Wgyd6B3B5zmMbJEEaSwKi4T5EErkjp9TowyozFM1xYuQ5xDEcKDxBObcQ5AwnsAM0++940jD7+L+HItaiEucONeo6LpdyNkN55HM57zi7AmC5dninENrNDwZkMJlsh969mBqOdmKf+LIGCf/8rTX2VthqgHhk5Obm8e9CJAv2oT+CnkSX+wYHfMVkGzRuh5LWWIUTxjbGw71WM43aYseDkCcDQpwOEDOBAW0VlRDxm4t2Bm/jVqOwzQ1ziDX5Yq8m86xSbZXnMcNC+dMWgH3g9mmhyCNZcfkNDnOWayztZixuNwoJ+6rCvs6MVb2v9pGjqTjFDlmKqdxaTjTbJjlS041VvimH2dbb3Nh3DHOxDlNyc/ZdSggUSkOCauSeUdEbFBEg4GiBC22Vnrxyq03CurjVOC/F39eOhgN6UHiXFLh859jgaJIFcO8sq8yT4wSjaZAB0nWSXodEd/nnFxGlFM40xirQGEyb4RGwpzjb9BR1KsBYb4u7Gs4hi7u6/jxMehUTMGlY5VzNAZN7dXwRoRNHB3PRSxNKmHbGEeJTTHpIt6rM35Ow0TjqWBaD4aG0eBs2hZRGxRpS/a51O+zsRrfo7ktNYmLJVyz+cReDun+HOVObe4W9zNfJd+Gw8NzGmg5LVrsCq3R8KRAUAZyHXfiR95yGvfDCXermuRMCFKZtG04atpdt+lF0Co2bttejtLx6ebS9yLH50FAJ5uPTPE0YGujwW+LJUbjKuAI3AhRuRchrnjHVTuvlDYexIlDzcGm6lTiO5xNcBwWiauaSqYqvLGiUNz1pNw3xq6VJYmVPaZ5GqIyG70IGVwwECTt5/SxwCtfkTNRRcR5r0wuJ25rjGXr6lLbeiOcjRKCYTHtGJqkPKK8sTF5LkwyCDzJBCM3XhtBGS4yDVMcSurxlVgKscQ+YMrWRkLTi2DoJkVMguLc5GhIK/uIUOBX9yOnG7wBhRRpOtFIqFftCxQddGzuhkLRiTsd9tS/jxzIDItMjsU0xvKj13PWKLRz6FQ+uS6STNovbwAU6TiH92JrjlgKZ+nHYyiScicax1AsRawcJKCDsVGPZZrnInzPZGNFIyL3jijxIUqRVCjDUJnGGoovpVovoESFXXXCcRNVGwAR4b0OpVsbjcDG7qvJYxoo413HxUrdAK7xo2mISZ7b5m8GxlkT2IbTnPN2nN3IOR85e4Q2POmcRGs0PCkQ4j3NTmXRWk4si7kjZ8dKEbvgpNXpHe60rlZmt6fYRM2ToJuzifvlgpehCclWpcdXyJuDubE5j4Ue0TQSPHlSuXbZCvi23gibeyMm5+QwmBiXvd0zJfNYpJCEh/E0bDfnxgymlBVseD52QC5nmkFA1jwsKkDTPCg2XRs7jBXPU9JXJw2mFL6Fnbpfja9bB2JRmdamxPnV6yzpejskwzT8X1B7VrxhUT9IZMp16sfT6bhga2/BOLtuLRG9P/VxVXSwlFkovKu37TCWn3Nt7EfPgx7j5D6Lpv+iaQTUcmzjM2utN2iyhf3kjYi/C+s9DePIz0WqgpXxoqdB4vWz471u7Pc+bhBkk3bOTvUQQObLsjZb/ZeG4ZAujfy8j3PG/t9pznZH1vh9bPtZW7vD/bDlPAxnj9AaDeck2qCzJxNyV/x2v8eW8zActQMnrsLuwInPYWmWTW1QsnCjaRyxmc4XlgVzPcKHIpF0+1hSNlfI0usYniQqKBS2VpjwQpyJlW1IxoGM7ZjLSyky+UBxuUEw5fOaF5ccM09G/CxLYPVDuWwFPn5k68/Ct7Y1iNIYbrri3JhzvR/jrxshQymkK+NMkTM5FZsUObeDIk/av2hITHpiXH5utzMIdlFWPSkGjbt887xHJc0hYTXbJIYW4xuwxTg5lSnN2/x2UAodQoLiamZU2KM3S9KKeietiis6wdNQr4BHz0LuZRznxG2SyRkfK4Ywjeca1RypjZ6Mo3Epb8Zz7ASnwCCp90Yue/zoxN+p/0CLzY5zUOTqm0tmMIyfr3ggJBkStTfWK/dxsUGQRu5Cmou4ZCjXfRFqTkxyTosnee8EqWUAmPgbDN6saagrI8mU0Kd4fMJusdP9OU5njzi7Ges85rRoMQ2tp+EsYU+bu9mJF5PLNuczJz2HdsNxO3DCw3EnzsN/kJSHpMiMlVOFbaKcxlCHzUy/q+eehnrJLwudSZrlmAKS9VdIpU3T6vj28xkPIZgIrcnyE6Qx7ricyVfjcA/DyD0EaT7T8hW2HWFyrF1znGFazcLpal4Tu3k+y9j/kzBAnvi9jVGVb07hSR7xtUrXxPaTTNfzjgfq4VeBa6+RJOV+kj25Lj3OqfXf6WNpClRqeRyU9Mx34F9X9VynrIpPjLWNBwUg9osUG4yqOhrIJ0I7G4yvJsYlqqlnvGb58KTa6JwGyXNXxhT8xuvdVMtJxuv0STupDXrZ6eJwEy+mjBWntcMKeHJ97HB/tlM+aDmPjPMYUDd3g9bTcG6iNRrOEvayuVudAJb9CKP/u+VkWtZj5UxRVcc54xOLMvN7cHq45fyxFcP4XN+Jkxkxkq/WMe07TVW8KXaKtwNpGEdpJTS9z7oZZ7kDAHaHsVLTuO04eTLwNg+r5gr89g+eeq7b42w9tsYNFNnlfKapp9sew6TPT56L6YXWx/Z2avJ/jQKT8gN2tmGn7NmESBlnjK3aFxSpN4GbygHfV8EzZFtOip/fkTO5ALA9x6PIOPWDURr/6SzfRo+biflto3lLQE/pqhUaOYfx6nOR5yn5CTWvlXzuU7EbZU3tYqVZPcznGWe7hQ3//V1opVnezMNxdrw/PxLOLp475yXnMSA1d/vAqxr9iR4r1v78v7H+57dzcHV1z2S2eHRow5P+nqCcTPztHpMP5MmVopazO47bnjNtSXWcM62l5tjqWJ2Amm9tclQRkzAbk268K2L4hShcUqbG5MQEy+yXPJFLEROYZSfOpOIzsRKqJsM8JjmyM0cUShcTcx6HTJnPBCcNsoNhsZNyFaAeBScPf4ivdzMfpca5k8dQxhQoQWWelywZnXFFuCkpKqpNSs2p0JgQ6pM/3yd+FrXlXm8cWxCO5z0/TnnYtKHySbljcrYLrW4M9Sg4Sk3OeZwz7XyNL3SrseptBpX2PXZAiMfHZRf0tnPOXo+PlebTIE3/nTZL77sxikxuH+fkx8dtxxnP7JjkpHvUDvdVLVPuh+NjTfvNjHGmGhTb3OdbzmPh7AwrbuJvEgJu7/4WXvJ8LnrnD7C0tPSI5tpi7yFu56yiFo8SIkJ34cDOpF0eeecs5cYpEEV3bj8Ao/WTAHTnD7ScR8URRusnAucggk+kLjc9pze3DMAwcHrzB8GBc4ZR5MwGzkbgzEU5ltHWCoiiHziDwOnPHUScYK1hOPCcmb6f49bWcQBmZg6lmPet4UkExWzvADjYGHnOXPcQ4hzWGTarUwiKuc5+wLJengJgvrOfmHy7adYRhLliHoD1as1zigV8SVPHptlEEGaLWXCODeMbiM3pWcBinWPLDhFgVvcBx0ZoRDane+B8rZ8t65vHzaouDsdm6CMwqzqAwzrHwPnY+Rnla9FvhnyA2aCoeI5fYZ8JSsdm0LJmw3vrHIOwxjkTFKXNcDubFe+Rcc4xwD8q+2HVdiv87mYy3arm+G1NjkzIkbGxQMbmM2XOoqbsuzSOjwTOlh2F4zwDSONciHTCOV3z57SzD1MUjetna3iSYm6Bam0VRKH3zyPiqE76814cWKAjBmctg5UBKJhdnkGLZe24P6cLh3poLNY6zpysEAX7DhQIjpXj/nwtH9IofLL+yklfEnT5gE9EPn7c50IcOhSSnq3jxEmDUnDwgPc6PHTc7/sFh/x7a+HEyfIxcRyuMdaBA378fD6eAycDZ/mAv+5OhP06eMj7EZyFU4Gz74DGIZwKcvYfKnBI4/gsHPC9Jc4c9300Fg91sfiOxxsnR75r9IEezgkbJ3xTwbmDfSyqcS76y75/9vDEBgC9g3Pe32Qs5alNUEJn/xwOaZxT5/xY9vSaz3tYWgTAnPIrtHr/EuIEZy3mzBlQimJxERCq0/6+Uezb7+911mLC9dNZXAIH5RnP6Sz6e6izlmr9tOcs7AOgPLMSOMvhvI9xHJRrgbOw7L2e1lBueE53PtzD18I9PDxHnc3u4S3nrHEeFrtxHF2ynyt+8727k/cIMPvPf5I77rhjz+W22D1ao+EsQemC1/7iV3fm7NLbsL7yAH9yy0voH7iMV/7EnwDw8X/xdMDyhl/6+nnNESesrzzAH//ki+kfuIxX/fifAvCxm68CLP/gF/8WsZ5z63teRH/5Mm56958hTvi9H3oaYPlH770nyDnKH/30C+kvP5XXvvNTAPynH3kqYHnjz9+LNor1U0f5xM99GzP7n8rr/9e/QBvFf3jXJYDlTT99lKLSrJ8+yu+97zrm9j2V/+7ttyNO+M2fugCw/E8/doLesGDt9FE+8uvPYn7xCr7nn36BTqn4xfcvAZabf2CNmXXh9No3+dX/+C0szT2NH3rN5xHr+PH/tAxY/vc3nqa/usrJ9fv5uc+9iOXuU3nniz+JGp7kX37mhYDlF170OVx5Hyc3jvGzd/5zljuHeedzfhbHiB/9y7cBjn/zvF/FmJOsbK3w81/9BfYX+/nRa96BtWu860s/D8BPH/kRrF1nZXCa9339t9mvF7j56f8QnOWWr34EgFuufjPWrHFysM4vf/NP2C8z/MDTvgNnh/zUNz4JwI9d8Z0Ys8bKYJ3/66GvsCRd3nbxEXCGf33/FwD43y6+Fms2OTHY4t+dvo8lhH924dNxtuTnH/oGAD9ywRVYM+TEYMCHNlZYAv7p8mFwhvee8sbZD+0/iHOWE4MB/36wyRLwP4dVqvcFF/c7lpZwznBiMOC3RtW2HEQ35exbRkQ1xlK6z4nBgN9cP8ES8JbDVyGi+DcPfg2AH73wGSg9y8nBJr9+8mssUfD9l16HKub56W/46/pdV7yconsBJwdneP/ffYL9apYfetb3oaTXOBe9/vMb5/RdN3yMv73hmfzyL+4HDP/85jP82m9ezrP+9Yf54tveQOeCp3DlH76Pudkz/OWz3gLG8byv/jrP33cXaw+e4d8/99dZuHyBt97+PVxRPMT3H/pdMPCBlX/C09SDrBzf5M1P/0sOX9bh4186Qp+K65b/GmPgr1efwyKGlZUh33rF3Vx6qeKzX34mM67gggNfxBg4efq5zEiflZUtLrnir7ns0oK/+fLzAZhd/jTGwHD1xSg6nDo14IKn3sZll3b5+pdfhKJDd/mPMAbM6iuxlKysbHH4is8ljqVsyDFUjbHu/NK1GAzLBz6PMbB6+jo2KVlZGXLN0/6GSy/T/PmXvgWAK5e/gjFwz+o1nEFz4viA77jqa1x8mebjXzrCkIIXL9+BNfDZ1es56eY5dXyT7376f+Pw5V1+8YsvxzjN/3Dwv+AM/D+nXstJt8Dq8U3effXHWb58lh/7/D9ixczzUxf+Ghj4seNvYcUuNM7FGz7zDioUv3PZLWDgv7//Fk6MFtl8cIXbXvhTdC9d5vl/egub1UzjnG5sLjI4dpx7Xv0DqIsPcfV//CBiFXe+5A1gHc/6zH+muznD4KGHuOvN34W+6CKOfPi36Iw63P76G8E5bvj9T1GMOmw+9BBf/Gf/kM7hp/Bt7/9/Eav4izc9H5zlJR+9HV1p1o8f4y9/8DV0D13MC9/3n1FW8Sf/43XgLC//8BdRlWL95DE+9yPfSffgpdz4r/8AgD/6X46As7zy//4KyijWVx7k0+98Ob0Dl/LtP/NfwQl/+LZngbPc9Kt3hfvzg/zFu26kf+AyXv4ef0/5xPd/C2B5zQe+lu7zn/rxl9E/cBnf8ZM+tv6//MDVgOW17/+b9Lz4s594Kf0Dl/GKW/4YgN9/+zMAy+v+7d3puXO+ch4Ou9Fb/ugDN/G033jfruQ9Esz8wC2t0fA4ow1PelIgZonVWybDLlrOrjjsxPEC8oTqCc6UxL9xjkuxCdtzUrlQAZuacG3DQbApHCULpRBVhy8rB6JBUiHJ9J5Ysz59osde+7/ESXKklgN1KVDwfFEZR4FoXCqNI3XYTrZPgia2/8pDWRp7LjrF6ORVPhpywsp9Y7/G9536DEh2PHKOxGOUy38YORL2dfw4uyRn+pxzjhIJx6+WU5+HOnREiGNl56IxZ3++nHLNWvmQYuOVAE6FHBlPsk4wTlGlhoRChaZCpVCRCuVDmFJ5VIVBUWX7F1+nNBWJYTz1b8hgsc5gbJ0kb5isRmUpsYwaB85SpuPjP/d/kWMpJ+XkY0mcZz2f+Hp8ucxMeV3nPdTHZHzfyc6gcZoq+335Y6oZmSgHKjd23cKUc6Ewrv59GaewTmFt8zc1fk5xCmwmP4Z4xG+MVTfzLRnitbXd/TD73oS+GDnZB9vKGc/VyDFZmWo7TiN8VG3DybEdJz//u3imnD+cPcIehifV13GLxxttIvSTACn+NI+93S5G9TzmyI4c2TUnvzltO59GCG/9xooDLeOUCTlpHq7+XuRYcT6HIcYvK3BK/F9qnCSgFKLr4+Oi4h6FSmi7FRVf57eJI+VGJOU/j9EPymudHO0V2WYpWhXkuOw74xwv22XfQRSi6vlE5T8p5aqD2DIrkxmNmKyCe3gvKhXKrPcxbhENCpS2wGbgqDGOH9fLGabxxzmIbsiJxkY+lohK14bDTchhGoedj08aa+JcaCQYgs4BqhteT79WHYD1SmdUFq1TVE5h03FzXnH1bdAAkhFhUiljR4VCZxHsVeguYDPFYwjo7Or3XUAMhM7AFq/Y20xNL90QJRobG7c5S8kWuu7VjKUMuRE1x1BlXci9nMZYDgxmwogx4qJujcMbCVWm1lYIAwpGsRyqBCPKZfvuVNj/+hjGY+bSMfTfQUXDS0J/aJ1x9Ni58EaCcSqd0/Q+GnDiz+H4OcUJSrLfhQtln5ORIODynIa4reZEBU12kdMQ82ZiVSaxUi+sBCNIC5PY9j6fbbfbcPLN23HYBaeReNVyWrTYDVqj4UmAVGmncfN1LWeMw46caStl0zmiSN1Mc45TDhNXucK913sJXHqtTGZziPhtth7LicMqhyvCamdYHByfj1XUT1vnwgJiPR+rBKsLXFAYRbzREN/7+Xa9QqvD6rYSRDrNvRaN0AUV69urpLzWHIVIJymnEt67fJU3rr5nq/9C+F46HB1EdVBBgRJFkGOaY6kOospwRCSMnVM64ExmoEjaVo/lP9NZ7w5R3UYPhqj8pw7M0csyfnxEN+WITkZBmo/UxpmISh6Cxr5LJyn7ItFYG5OLbvQV8HLqsZAOIl1E+3OhlOBET1yHOLLeH14JtFbXCqbVjFyH0gWDI3gVKqfSZVgFIyJ5ExD/uRRJzxhSMKBkK1zEFqFCMOISx4hj6CqGWVM/M9Y52YQOybU3S5JhETeVzht2keMQSjeakFNha8+CwJAKQzYfHEOimeh3d4DCINl+Nb0q4oQhBUM6GafDgA6DsF9+3/0xixzvqdGZHMfAdRi6TjrOQ9fBZMYHeC/QyNWP55ErsE7hMu9G5TQuMxpc9DKkG4rzyr+tOWJVMCQyj0UyLEj7CqCycyHTOE7qRf/4eTygNHlA8tpMGDGQNYdT2Q10/P4cB89v4ts9L3bDoeVM4ewJHOkaa3FuoT2rTwrEh2R+Nx4/dS1nZ06zbKhvVFSvkvm/piIzVc60cqFjHGdskzvGseKosvIqVoVVz5jgG579dZ9b542EQlLoiy0Ep1S9WurA6QJUpw6PUR2QGZzxiqF/5mqvvKaGW17hjU3ZHCB0ETrJGyIETgqfiKvrnTowIyjFzmUr/iGxNx2CEI5j48PESW2QZHIQjU3hUl4pb4b8eMMihkL5SJ/OhCIvqoPNV11FoTKOUp3gMVCNbROccTnT9kt1svAkmTofRGXr3S7IyYp7SuHlpJVr8QYCKjsXCtQM1oakYASne15PjNdzNGiDkuUsiCkwpkgcYwof8pJ5vIa2YOjqRmnDoNwOXe01GdKhcnXjtsophmiqOLbzCvewng5DvKJeZT1FhngjQiVOxZCSYTBErbMMKTGYxDGYBscF42BczgjDMHRPdM4xEstQXD0fcVRIHSDl8PuASnIqFANXMAy/C4tLRlR0SFZojPN/ad9dd+IYDlyHgY1VqlTyGkTF1jjFyHaoXPQeuWQkREHeM6QpbX2crVNUpkjV0SpTgNVJMXTBYBAndQW1sM2V9X0sGQSqVurFCbHptHOZ8j/mja2mrVJPhFlGo2b7+3O6L6fmh5LkxPmkFhO7ehbYXXCeKM+mJxpnL7AHoUhteNITEq2n4SzBOceXfvtdHL7G1y1+TIg1+vP783i9vvOZ43bDYdectHKnXOI45cDUK2TOOb/NZnwJ20JohCDeO6Fc4pjCoo2CIjMaxGEK2+CYQjJvhFB1QVd1bwVTCFWvgymLmtPrIabKODMUehGn18NIBaLmcW5UPzRkBiVzCIHjBFF9hJojqo+iRNRm2C+FqFlwtVkjqu8740rfvxdBSR+bN+nWsyhA6+ihEJSebZwHpfz3ta6bzinVq3MlAFXMYs0QoofABU6ei6L7iDOINo3v5f0VRPcRUuSIH033ySG6742CXI7uNQwdpWeDkRC1LBDdQ6h3XqkeSs+iVNx3hVJ9RNXjieoj0gNX56CI9BA1l52LOdCLoM+EsTS2MxO8C57jvQ7ir1cAEVRVYMpezSl7rJk+m6H6lUXYsj02XD+d0w3XZ50+68QQHcU6fQbU3qpVZimwbIU8CiuwQRcYJZt5A00lljVVK+CbMewncDbFejkSVu3F1a8DZ028VyFtx7E2xtkUH/S0FQ0nHBvB+1HPR7FBh41ojIiwQdeHHgXOquuzQZ/NZFMJq8w2js+qm2XVzbFKNLgV666XeQz8Mdy0PTatP87GwRk7yxk7k87pGTvDpu2yGQxBJ4qB6VLaIoWgb1R9RqbDyNQG5bDqYcte2nlb9lBVgTLBa4igqwJlak+DLguUVaiUYyEoo0JYUbjGrCBGUiNwJPMqZPdaZRWqisZs7lWI2n1Q9tNPJ/NG5Eq9y67V3BU8bpCklIacM/68yGQ+HCff3HJ2jeklVj2OfelWjn35jx+RvBZPLrRGw1mCiPDcN/11tREnAAAgAElEQVTM3giLse0NQ7sZOnNec2RnjlOOsQax28jZmWOzOOjJhMD4XZdiwl2oYT1ebcKqZi9lpxy5J9cpMNpRJYXFURXUYwNlFzqFokr9Hhy2UNhOHZ5kO11cOQfFTD1pNY+EkA4AUfPgthDV8+9FEJlBpA49iu/VmEHQWKGXHqgSrQZpLFG9pNd7Th+rLE6GiaNUr6nIqx5+0TyEnYggutfY+cTRSatJRkIaK8qNHXRRwbBpcvwLkzhpW84RPSkn9yLovs+LiPMRbxDkK3iie95wSPL98ZFsPP++nwwnz+mn4x6PoSvmoJhLHNPpYnQWsx9exzAnnEUZjak69eVcdRjYLkMb8yEkKK71fDaDEbEVjQbrleAiu3o3XJ8+JethedtZWHe9xu9jgy4VVZJjnLCRehxEjqLAsZl5LDaSdyVy/KsmZ1KOQdgkJu/DBkUjgXmDDhuuy3pYSjcWNlyPKlf26bNOj40w58pJMqLiHWOdPuuux1Y4HhbPMVnew0Y4hpvpNyds2N7Ecc7PhbXC0HSpnCYkDjGyHUZVj1HVTce5qjpQZaFrVQdlNKqKnjxbGwRxdBuNhOh58Mq/5L+vGMIUPQbGNhOmyV5H48wFYyO/ncbXKamhNkrSmoVtymHHZ0F8kckJnAk5DWwjZzdjnY+cx4DDR17B4SOv4Nb3v6Zp3LU4Z9AaDWcRO1nksPuSq4w9PFtODR9WxI6cHAqyfIXxoaTmjJ07J86HiESXOeINgoay73sQuKKekFMOk92QjbZ+7G7YpqTeFvag0hZdCPTitISq69BVzTEFlH2FDZ4GtKLsa8R2azm9Lno0i+vMJQ7FfrAj0oNC7we7haiFMJhGZC7sl+cotQjOoPSW38/kacjDb+ZQgNKx+o1CqdkUsgN+RR5AB44SjehZJPc06FmcqCTH5xzM4jIjxnM0RVCORWuUniGHKuZxrqQogoKrNaqYaXo1Cj+fonDp+MRtDY6ohhxdzDcMAqVnENWhiN4RpcO+5mFPs97TErwsKhwfpWYbHM/bSudCyZw37NJxnscUM9g4T62oeh2sqojnPVbjKuL1LApdFpiiIJ3TqmCz6jPMPA0bps8ZWx/HM3aGVTvLpo0eAs2q856F+GNbdXNoMaxJxmEW41TGmaFHyVoKgVOsMEOBTb/dM/ieEOsh+dYgrBIU5CgnvM85Z8KPJJczcAVrYZXAz7nvDYJsPhv0WEuJ8YrTzPowo8SZ9UZBCmFUrLo5bxAFrNpZzuTHxwln7GwK5wJ8KVUzk/p1WCes2RnWTT/paeumz1o5y7Aq07kYmB6V1el8DaoeZtTHmbkwZ40d9VFVETwB/pzqsqATwhEFjao0yqrE0VWBWEHHvCJAmeB5iHKMD2EqgnqgVBbmlJR0hRhFYWM1MJJnIXHC+/icE6S+74axgMCpjT/Jb+iNNZ16u9jpnHQCpZadPqmTzXh4nM+cnbEbvcWfx8c8VIsnIFqj4UmB6O7d6VfYcnZTBWK7tiROXApvmeiWHDnK4SRyrDcSspV0Kw4RR/QjOOumGh9OuTqEyFlMYcMDrQ5PKp1QFXEF01Hp2HnTc8qORXcVo2B8WCyjviC2qDn9AlXNUVb9dJRMdwmxZeLY7j6UK3ESE0oVopdDzcnw4JU5lDIo8cqsKEGpRVzmsVBqwR89iSFMGqUXcDbPD5gHFKKCp0EErRcayclKz+OkAypbtS/mwdZGg9bzWMnCgVDoYlyOD3uKSc7KhrHIw5W8Mh5Dj8RK2pZzBN2QE8OR0nyKhZCPkAp1ovSYYVHMo/VCtl8apedReiGb8wJKLaCCt8bnRiz5kKQoWS1S9WcoyxjKIpR9jdEj0vUTjCkXQ1CsQ1UaXXaow1Q6bFZ9RmVMLob1apbTepDGOl3Nsai22AhyrIWTdoEiS04+aUOzwBASUznhpFtoJAyfdAv0pGRduomz6oLRkzjzFBg2bMxpEFbcXIMT3+eck25+Qk6FYj0Yq5UTTrLgKxpFDvPeaxCU+wphxS42OW7Bc4x/RBornLT+exErdiEYDSpxTpk5b3zEkCozw+lqlq3SpDmfqWbZqPpEjXetnGVgeoyqYAA52BrNYE39eB4OZ5FRDxkFT6KzqFEPPerW53TURVcaXap03nWlG4nH3kDQiKnzMMT6bXHOymiUrRV877BrLuHnOQbx+E/zRnjvQ5Dj3Lgbl/EVH5cZErkcP4RLc87HHRM4yZmAG/u/5bRosVu0RsNZxO49CQ+HGB7Tch6OM9VDkCnkKeRV1bxxz4Nluhz/fzQaZKoc5zJDQYIXIbsOrHIYbTEx3EVohJYAVB2vYAx1rDbjGPUNKpvjqO+/s7kVOcJg1qV+DgBb8wL0GMYVW6UYzS0g1hIVltHcPrqicVVU5Atc/xJkNIdfKbdI51KojiGh3oyzCq0P4FyuyB9EpIuOc3YKrQ7gZJBxDvjQHl3H9Wu9nzwXQRfLODtABTniNEWx3Byrs4zYARKTOdHoYn/DINCd/eAs8WtO6bCt5hTF/vj1MB1db8s5oifk5DkNuvDv6wqVOmzTmZxllJpPngbnBF0so1RtNGi9jFYH0cmA06jiMKjF+lwUT2Gw2GMYQpacEgazMjU8SeroLYqyeavvDLpsbi0wHMbkf8WZwQKFVLUCPlqkqyq2gvJaoniw2o8WmzgPmGWGrsOG9fOpUBw1y8zJIHHut8v0pWQrGJkW4X4XugVHjlumwKQwHuNUzaHmANM5mRzjdGOso3bZexoC56g94MOGjDeCrVMcNfsbVY+OmmUfnhVW7Y1TPGCW2bTdJOfBah9rZoatMno+NA+VS3WiP3C8XGStnCX8TDFOcWq0wKDqpT07M5pnOJhjsBUMC6uwmwvNkKGtOYphj2oQ8hWsohj2KEb1ee0MOxSjDqNh8CI48ec9W/3XZYEuNTKKoWugSx3yHvymaFi4GB8pDjEq8QHEqNBZOo7ufJjTlApLLvcupO01R5zUN3jZpgqTlawa2HSMGzLbh8e4cfEt5xHg4SIo6qH2LjzpzKdv48xnbuNAaLrZ4vFDazQ8KRBdt9mPcMIt23J2xZHpHJ9XIJOccWSrx9MMC5/3UD9sc68CgA0KnYmJ0FOMhrRS3LVp6qNu1XgIDnsVUGBmXOIMewYrirhsN5ixWKXYDEqM08Lmvm7j4bu5r48tNEMJN2OtGSwepDOoV7er+UvQg3lsFaoDqT5SXIZg0kFUnaci1fG02q5UF11chLMbiaP1BYjbqMOKpIMuDjX2vSgOY+0mRZFzDjaNhuIgypXoTgi/UQW6c7App+M7RHcSp0MxzukexjmDDvq/koKie3iCA4oicLTqUHQuIC+nqjsHUdJryikONa4VXRxEqYUU5qRUh0Jf6EO/wvEp9EWozqVJtFId7MwV2E4dMlTNX8LmvLAVvQgKNudsuq6gvsZSooxTFKNOI3SkM+qgNxZQm7VxZteXWPHvAMfKljegBqVfAbdOcd/oAF2pEufoaD+bRZehjQq45pvlAXqq8sNbuKc6zJwMGBhvQFZO843qcD1FC9+oDlNgGBiv7JdS8I3qMFrqEKZ7zSGMUw3OfdVBr+xncjynHuu+6iBD10ly7ikvYOg6jFzkFHyzOtAIT7q/PMCa6TMovUekouDoaD+DFP7neHC0n42qz1YVvBFOcXywPxgNfkLHB/vZGs4yHMQQJs2ZzX24US/JGa7vQw/76K14LhTdzdlGzH53c4Zi1MGO/LWgrNDd6qGrOoSpM+x6j1JVJ0KrUjfuG7rU3htR1ffDOoQJcARPhKoNBVS9YBE5VnwidKx0EJR9sZI40ctQl950bBd6NB42My4nzmOC+UR4ppxznD3AmNfpsWLpRS9k6UUvpPMv37lnMls8OrRGw5MCtvFfy3lknLycKnYnb0QtJ3kQwg01Kv82q5BjYy5C5IjDFAaTJejGECMfwwtVYRBtKTtVFEPZrXw32MAZ9kq0UQy6YVVaYGt21JCzOTdk2C9Zi6ul2nJmaeC9W2E+pw5uUZSa1SIoR4Xl2FNGwQPmb+gPXVzSGSlOrYSyjV3hgSsKOuUi/Fc/5ANXLdIbLHHytIG/AtPrsHLltSjr4LN+sNNXfCudgWGw+jdwJ9hilq2LX0UxrPMnqgtehx4NQP0N8NOgFigW34BfrvxRAIqF1yHVOh3zt8CvIHqW3txrwW4BPwc4+rOvwtl1uuYbwL9HZI6ZmZeFY/5vARfeG3rmHuB3EGaZmX1p4PwKADOzL8U5S9/+HfCfEJmjP/PiUE615oCmb+4BfheYZWb2ZXj3RD2WqBm69hvAhxA1R2/uVSBd4KcAR2/utbjOPjR/A7wX1AJq3+up+guN47O1NMuZ1YPwZbCdGR58+oWYwjXOxenlQTqnpnCsLm8x6lbp2ogGZhEUQ2Wht9HHZhWg+msziBWKlaCoVorFhw5Qrtbz6d1/CadmDjEMq3vWdPjG0WfhlA2rxo67Hng2dAeUZ3z1raHp8tnjz6EoymQR/MXJI3R0SXXGh65t2S63nro2zMRzbj11LUos5Rl/PW+ZLreuRo4f67+e/lask8QZmC5/dPq6hpJ+66lrfSnSbKw/PnVtyA/wnD9beQ5l1aVa2/Bzdh1uO/bcsDLqObcfOwJVh+Gqr1JVVQV3PfDsOocAx733PwtVFY3jM7rvyoayZO57Gv1KI6dmw3HWzB59il/ZD/u1cOwQutLpXCijmV9ZbBj38ysLXpkP+45TzKzOeWU+/N77azPeQ7ARchoM9Db7fr8Cp7vZ9wZH8FCIFTrDTmPOxchXXOqGECZlBT3mqdJl4XMRgvEqRqGqpoGiwvdjZ3IfBtUMc4rvJVbJy8KicnjjI67wyIRCOlGVaUcluOX8fUQntYnQ5yZao+FJgfDwz1trTvzoW85j4XjDQhInGRau/tz/3+Q47RLHhlCktJim8lVflzjgqyMljrIh+9pvM4XBFIayGw0Ox6hXNuSMeiXihFE/hJcIbM0NGx1ZN2eHKCdsVdHggDP7NsMee87p/RuIE9b0ZhjbcfKC9UZI17GnrCNWWFvxCkvVcdx7ZZTjcc/Tt9BGsXLCwK1Q9uHrzzZ+1fMTnnP3czsUpsfJ44twB4xmNXdff5H/8FNeWbv7hovRleL4cYEvQzkzwz3Pu8YbOrcFZe36G9CVcPz4V+BOqPrzHH3eq72SdbsX98B1r0ZZWDnxFbgLzMwSx577D/yHt3s5x577D1AWTp34EtwF1ew+TjznH/lCTRkH4OSJLwfOEg9c97ogpx7LFI6HHjoAX4Fqdp57rr/Of/gZL+ee511D2bWcfOgMfB5GMx3uvv4glbaN41N1BqyeDCvpXXjgkvVgwNbnouxWbIZcBKf8eTaFaVw/ADrGxCvobnUboXS9zT561EGtejnaKeZWFoIcz1k4vg+rLcO1bpJ34O8u8r+FcFkfuOcp2MJSba7zdUCZgqW7r/SGRZAz+7Vn+LE3vZKuqi7qrmc35MhXr/H0wJGyS3XXkXCpBqPqriP+ugwcyi7uriPe/M3kFBlHVV2Ku57lH3RhPv2vXc2sVVSRU3bY/7dPDWEynrP8t09FjDBcW+UbgDIdDtzzlMDxgpbvO4xYlR0fzb4HfbhUHGvp2H7EKganY06MYuHEUiMcaOHEEjhJ50I5zczqbCOMZ+bMnL+3bPYCR9Hf6Dc4vc0+YoVqFJswajqDkAMR9isaCLqKVdak9kIFOcWoE45d9CIQKizVHB0+S4aBuNogiOstwdNQczJFMg9hguSN2IlTe4upsd19XrWcR89p0WJ7TC1Q1uKJhe6sj3u+8Lk3pW0LFz8TdLflPME4Thyd+ZoTG8dFTnyfc6z24SULF38L6G56XywGzrWvToZE5JjCUHUq1NJs4NxE1akoe6Ufq+hS9kqG/RFuv483P/ytNzGYHTKYHSbOYHbI1twAc8hzLrzuVWwsbLG+uMnCU74FKbqsL26yvrTJ6MIg53k3sbq8zuryOgtPuRopuqwur7Ny6AzDy3qB8ypWDp3h+EWnWbjoGUjR4/hFp3ngkhXWnuGV2Quuv4mjl61w9LKVxDl62QrffOpJzjzTcw7dcBP3Pu0k91x5InHuufIEX7/6OMev9crIwRe8iq9ffZy7n3mc+QuvRoo+dz/zOF+75jgPPs9zDr3wJr52jd8WOV+75jh3HTnOgzd4JemCF9zEXUci5xmJ87VrjnPseUXi3P3MMNZF9VjfePoJVp7rr4FDz/dzvvdpJ8Oc+9z7tJMcvWyFM88MYz3/Jh64ZMUfn6dcnY7PqYNrbF0ajuH1N7G+tOnPxcXPRIoe64ubDPsjWJ4J5+smyl6J1ZbFS54FRTcZpv2+Dy+66NrXUIw6dAZdFi95NhRdOoMuM+szLFgfHnbRc15Df32GudPzLF76bNAd5k7Ps3BykX0bFyTO4vF9LD20P3GWHtrPvgeX2b9yoeccebV/f/Qgi5ceAd1h/9GDLH/zEAceeornPPsm//6+C1i87DmgOxy47wIO3HcBB49dHDiv4uC9hzl47+HEOXjv4QlO/F4c68B9F0yOdf8hlu8/lDjL9x9i3wMH2H/ywiRn34PLLB3bz+JlYb+O7Wfpof0sn/FG7VOO3OT3+/i+JGfx+L7m8TlyE3OnFpg7tcDipdf4Y3hqgdnVORbL+jjPrs7RX59h4RLP6a/PhHNxMHBeTX9jht5m359T3aG32ae71WOOmtPd6oVz6jmdQZdi1GFGLwfOTRRlQVEWLF7yTNAdirJAV5p+Z1+Soyuf0xA5PqdBpevnwue82pdpdZI4MRehN7Pkr9Xnvjpt85xuUva7s4uBU98zIydid5yFCc529+fD17665TxKzt5AaJu7nZsQt105mRaPCaI0l73wn+xNc7cWLVq0eAT4w3d9G6/8mc893tNo0aJFhr0rjvLExINfupVjX7mV4yc+z9W/8it7Ll//6I9yxx137LncFrtHG550lrCnzd1atGjR4pHAnfsKSosWLZ5YuPCIXyT95Ptf83hPpcVZQms0tGjRokWLFi1atNgztInQ5ybOG6Ph2LFjvPe97+Wiiy7iB3/wByc+/53f+R0++clPctVVV3H77bfzfd/3fdx0002PmNOiRYsWLVq0aHG+YrzRX4tzB+eF0fCJT3yCj3zkI/zWb/0WP/ETPzHx+Qc/+EFuueUW7rrrLubm5vjmN7/JNddcw+///u/zkpe8ZNecFi1atGjRokWLFi3ORZwX1ZNe85rXcMstt0z9bG1tjR/+4R/mu7/7u5mb8w2tLrnkEl7zmtfw9re/fdecFi1atGjRokWLFrTVk85RnBdGA4DWeur2P/iDP+D06dO84AUvaGx/wQtewBe+8AW+/OUv74rTokWLFi1atGjRIoYo7e1fi8cf50V40k74whe+AMBll13W2H755ZcDcPvtt3P33Xc/LOfZz3722Z5qixYtWrRo0aLFExyTXbsfC07f9mlOf+4z7Avd11s8fjjvjYYTJ04AsLCw0Ng+Pz8P+ATq3XDG4azhD//Vt+35fFu0aNHi4VBunG7vPy1atHhcUFx0cE/l7XvBi9n3ghfj3vVDeyq3xSPHeW809Hq++6pI0ypWoc96t9vdFWcCopjZf1F6e/mL3sTlL/6uPZt3ixYtWmyHP/xX38ZNP/PfHu9ptGjR4jzB3336t/m7z/wHAErKNgfhHMV5bzQcPnwYgPX19cb2+P7iiy9ma2vrYTnjEBFe9iP/357Pt0WLFi1atGjR4omEy1/8XWlh9NZfeu1ZKbnapjU8/jhvEqG3w/Oe9zwA7rvvvsb2e++9F4Brr712V5wWLVq0aNGiRYsWLc5VnPdGw3d+53dy4MABPve5zzW233bbbdxwww1cffXVu+K0aNGiRYsWLVq0ALGy538tHn+cN0bDcDgEwFrb2F4UBe9+97v58Ic/nMKQHnjgAT72sY/x4z/+47vmtGjRokWLFi1atMDHEu31X4vHHedFTsOtt97KBz7wAQB+93d/l2c+85m87nWvY3FxEYCbb76ZXq/HW97yFo4cOcIdd9zBhz/8YV7/+tcnGbvhtGjRokWLFi1atGhxLkKcc639dhagdMHr/8+vPd7TaNGixXmItnpSixYtHi/88S+9lme/99/tudzqJ/8Fd9xxx57LfbxgjOGBBx5gfX2d4XBIr9djYWGBiy66KFXnfKLhvPA0tGjRokWLFi1atPh7gNvb5m7nCu666y4+/vGP85nPfIbPf/7z3H///VhrydfuRQSlFJdeeinXXXcdL33pS3nDG97AlVde+TjOvEZrNLRo0aJFixYtWrTYM7RGg4dzjg996EO8//3v58EHH+TlL3853/7t3873f//3c8EFF7Bv3z727dvH/Pw86+vrrK6usrq6yoMPPsidd97Jbbfdxs/93M/xtKc9jbe//e286U1velz3pzUazhKcc3zxP/wrDl/zCi488orHezotWrRo0aJFixZnDQ9+6VaOfeXWx3saTxh86lOf4h3veAfXX389v/RLv8SLXvSiHflLS0ssLS0BcOTIEV75yldy8803A/Bnf/ZnfPCDH+R973sf73//+7nhhhvO+vynoTUazhJEhOe+6Wce72m0aNGiRYsWLVqcdVx4xC+S/vEvvm5PO0Kv3P7nnLr9L1hcXd0zmWcb73nPe7j77rv52Mc+NrUB8CPFjTfeyI033sjdd9/ND//wD3PjjTfyjne8Yw9m+sjwxMy0aNGiRYsWLVq0aPGkhLi9+ztw/Uu56m3vTKvwT3S8853v5JprruFDH/rQnhgMOa666ip+7/d+j7m5Od7znvfsqezdoPU0tGjRokWLFi1atGjxGLG1tcUb3/jGsx4+9Ja3vIVPfepTZ3WMaWiNhhYtWrRo0aJFixZ7h/M0EXpmZuZRGQxaa4wxj+g7L3vZyx7xOI8VrdHQokWLFi1atGjRYm/g2upJjxRPlpZpbU5DixYtWrRo0aJFixZnGUePHp26XeTJYWS1noYWLVq0aNGiRYsWewKhbe4G3ntw5513cvLkybTtl3/5l/noRz/6OM7qsaE1Glq0aNGiRYsWLVrsHZ4c0TZnFW9729s4duwY+/btS9u++MUvPo4zeuxojYazhLa5W4sWLVq0aNHifEHb3K2Jl73sZbz5zW9ubPvYxz72sN/7jd/4DT73uc/x7ne/mze/+c3ceeed/Nqv/Rr/+B//47M11V2jzWk4S4jN3VqDoUWLFi1atGhxruPCI69ITW3FyZ7/PdmglGIwGDS2nT59+mG/d/LkSb73e7+X973vfXzkIx/hq1/9Kvfcc89ZmuUjQ+tpaNGiRYsWLVq0aLFneDIq+XuNxcVFnvGMZ3DFFVegtcY5x1133cX3fM/37Pi9yy+/nKNHj/K2t72Niy++mL/6q7/iqquu+nua9c5ojYYWLVq0aNGiRYsWewPHedunIccnP/lJPv7xj6ecBuccP/uzP/uw33vhC1/IZz/7Wa688koAPvrRj/LGN77xrM51txD3/7d359FR1ff/x1+TlZiEULYgQUIFBUkCKFuEABaClEWpFKpUqBUp4hGkfotQoGx+lV/dihsu3xapkQoirRhFQQ0ISE0gYoAgYFEMREKUJRAwLEnu7480IwMhmST3zp3l+ThnzmE+88nnvjPkzr3v+Wy+sjisj3E4HAqLamJ3GAAC0PnTRQqNbFRzRQAwWViz5rphwT9Ma+/oZxt0dNtGRR3Yqf/85z+mtWu1ZcuWafTo0S5lu3btUkJCwiV1a9rc7dNPP1VoaKi6deumtWvXat++fRozZoxiYmJMj7s6lvY0nD9/Xvv379fx48cVHh6u5s2bKzY2VsHBwVYe1is4goI1aMEWu8MAEIDWzuzB5w8AW6xfOEyOcvN6Gppef5OaXn+TTv1lvGltekJISIief/55JSQkOO9767rk6sqVKzV+/HitXLlSTz/9tF566SX93//9nx566CGzw66W6UnDoUOHtHjxYq1atUo7d+5UUFCQGjdurNLSUh0/flxhYWHq16+fhgwZorvuukvR0dFmhwAAAABb+ObEZbM99NBD6tOnjz777DNnWV2XXO3Tp4+uu+46TZ8+XQ899JASExO1Y8cOs0J1m2lJw7lz5zRr1iytX79et956q5566il17dr1kqTg2LFjyszM1IYNG9SvXz/dddddmjRpUkD0PgAAAMD/PfPMMxo+fLhL2bvvvlunto4fP66DBw8qKytLS5cu1ZkzZ3TkyBEzwqwVU5Zc/e677/SrX/1KvXr1UnZ2tubMmaObbrqpyl6Exo0ba8iQIXrssceUlZWliIgI3XHHHTp9+rQZoQAAAMAmDkkOw/yHr7k4YZCkYcOG1amt5ORkPfHEE3rzzTd14sQJzZo1S99//319Q6w1UyZCT506VTNmzFCTJnWb+PvVV1/ptdde07x58+obitdwBAWr9Y2j2NwNgMcxpwGAp1Vu7nb80E51e3iZ6e2ffGacsrOzTW/XG9Q0EdpbmNLT8OSTT7qVMHz66acqLy+/pLxt27Z+lTBIbO4GAAACx4Wbu+FH//znP3X27Fm7wzCFR3eELisr0/Tp022ZvAEAAADrsSN0hZKSEs2dO1d5eXl2h2IKy5KGZcuWqVmzZmrevLlGjRql9PR03XjjjXriiSe0bJn53VYAAACwH0lDhYiICOXm5uraa681pT27ey0s26dh5cqVeuWVV1RcXKz169dr3LhxCg0N1W233aZDhw5ZdVgAAADYyFdv8r1ZZa9FUlKSaUlIbVmWNHTr1k233HKLJOnXv/61XnjhBb399tvatGmT/vjHP1p1WAAAAMCvVPZa2MmypMEwDJWWliokpOIQoaGhGjlypEaOHGnVIQEAAGAnwyGZuCO0rzp9+rQiIyPtDsNUls1pmDhxoubOnauioiKrDgEAAAAvY+ZchiOff6zdr87TiRMn7P61auWee+5RaWmpS9nOnSYaHTEAACAASURBVDt14403ql+/ftq6dWu92jdhx4RasyxpWLVqlZ577jm1adNGw4cP19NPP13n7bMBAAAQeJp1+Zk6/ma+YmJi7A6lVm644QaNGDFCiYmJmjVrliTpscce07PPPqsNGzZo7dq19Wp/zpw5ZoRZK5YlDatXr9bSpUu1cOFCtWjRQi+88IKuv/56NWvWTA888IBVh/UahmFo+/KZOrwzw+5QAAAALHV4Z4a2L5/53x2hWT3p0KFDmjVrllatWqVevXrpr3/9q/bt26cbbrhBktSmTZs6tZuTk6Np06bp+eefNzFa91g2p+HGG2/UrbfeKkm6++67JUnffvutMjIy/Ga92upUbu4GAADg71okDVCLpAHa8NStPnmTb7Zrr71WPXv2lCRdffXV+tvf/qYffvhBwcHBkqTz58+73dbevXu1fPlyLV++XKdOndLgwYPVsWNHS+KujqUToc+ePavw8HBnWVxcnH7zm99YdUgAAADAdtHR0eratauuuOIKFRcXa8CAAWrUqJEWLVqkrl27qqCgoNqfP3DggJYvX65ly5apoKBAd955pwYOHKhnn31WkvTvf//bE7+GC8uShnvuuUczZ87U7Nmz1ahRI6sO41ErVqzQRx99pHbt2mnr1q0aP368Bg0aZHdYgM8od3h+4lZtBPHtGADUj+HZfRq89d5sxIgRGjJkiL755ht17NhRERER+vzzz/XFF1/o7bff1owZM6r9+WnTpmnDhg168skndfvttyskJESPP/648/VevXpZ/StcwrKkYe7cuVq0aJH+/ve/q3///urfv79+9rOfqUOHDlYd0lKLFy/WvHnztGfPHkVGRio/P18JCQlavXq1UlJS7A4PcIu337TbzW/eH4c5vwtJFIC6cJR75jjefG82adIkLVmyRE2aNJEkFRUVadKkSdq8ebNbP798+XKdPHlSH3zwgVauXKmrrrrKZcWk77//Xs2aNbMk9suxbCL0uXPntGPHDj3zzDOKiYnRE088oY4dOyouLk7Tpk2z6rCWKC4u1tSpUzV27FjnmrutWrXS4MGDNXnyZJujQyAodximPBC4ggxHrR8A4K28/d7sjTfe0MaNGyVV9Ib06NFDhYWFtWqjYcOGGjlypO644w61b99eV111lZYtW6b09HTdfvvtVoRdLct6GhITE7Vjxw4NGzZMY8aMkSTl5eUpIyND+/fvt+qwllizZo2KioqUnJzsUp6cnKw33nhDubm5SkxMtCk6eDNu1GG1Km/ujYpywzB07vRxhUc19nxgAAKUZ1Y78vZ7s9zcXG3btk2DBw9Wo0aNtGnTpnq117RpU/3617+WVDHf4fvvvzcjzFqxrKdh8uTJ6tq1q9LT051l8fHxGjdunP73f//XqsNaIicnR5LUunVrl/L4+HhJqvcGHfANfLMPb3Th31tp+Xl9lvagDKNMpedKlPXSPcqY18+WTYAABCZPLbnqbfdmZWVl+vLLL52P0tJSderUST179tTEiRNVVlamefPmmXKs1q1b65FHHjGlrdqwrKdBktq3b6/27ds7n+fm5sowDCUlJVl5WNMdOXJEUsVM+AtFRUVJUq27m2A/bujhL05/9402/nmYWvb4hQ5teVvl50skSRseu1XnTh1VxE9a6tR3+3X+hxNyOBz6SZsuNkcMwK95aCK0t92bnT59Wt27d3fuw3ChDRs2SJL27NmjF1980ZTjDR8+3JR2asOypOHrr7/Wtm3b1LdvXzVv3lxSxZKrGRkZWrFihYYOHXpJl5K3qlw21uFwPQmCgio6asLCwi75GaO8TO/+/tof64Y2UHBoAwujDBAMswZcGIah8vNnlL95uUv5D99/LUk6VXJCHz86UI7gUDmCghUSHinDKJdkyOEItiFiAP6m7PwZlZ8/I0lqGHedae3mfbJcX330fzpfUqxGV7jestbl3sxKDRs21MKFCzVu3LjL1nnttdeqLPeV3mBLk4Z7771XRUVFateunfr27as+ffqod+/eGjJkiB555BGfSRpiY2MlSadOnXIpr3weFxd3yc84goI17OkvrQ/OT/DNPy7H8JO/DSPI3N+jvKxUX615Sfvef8Gt+rGJA9Xzt4skSdtXztbZ4iPqcbc533hdzOzf1Q7+8ncHeNrGJ4bLUW7ON3xteo1Wm16jJUmHXv2Vy2t1uTez2rhx45SVleXc1O1iY8eOrbK8vLz2y01t3bpV3bt3r/XP1YdlcxpSU1N15MgRffbZZ5o0aZKKioo0ffp0XXPNNWrYsKESEhKsOrTpunbtKkk6ePCgS/mBAwckSV260N1fHcb+By7DYdT74dF4gwzLHm4dvxbvxw/fH9S+1c9K5aVutV248wPlZa6QJCXdNkfdxjxT45jhunKUO0x52MmqcdhAIPDEnAZvvTdbtWqVXn/9dUuPsXz5cr377ruWHqMqliUNUkWXUZcuXTR58mS9+eabKigoUHZ2tqZMmeKcAe4LUlNT1aRJE2VlZbmUZ2Zmqnv37i7zNgINCUFg8IUb/pqYdWNvakx1eL8q64ZEx9TuWOWlOvzFOjkMh4KDwhQcEn5JHW+7AfbmpMJdJBqANbz13uz//b//py+//FJjxowxfV7F4cOHNXbsWH399deaP3++qW27w9KkoSo33HCD/vjHP+qll17y9KHrLCQkRLNnz1ZaWppKSiomGRYUFCg9PV1z5syxOTprkRD4F3+4+a8rb7zprM/N4sm83Fof7/TRPO16989u1/e2vwO7kzxPIbmAbzO/l6Gqv3dvvjebN2+e7rrrLg0cOFD333+/tm3bVq/2duzYocmTJ2vw4MEaN26cZs6caVKktWPZnIZXX31V4eHh6tu3r1q2bOnyWrNmzVRUVGTVoS0xZcoUhYeHa8KECUpKSlJ2drbS0tI0bNgwu0OrE274/Y833NT5GrMSB7NuYGtzI9j02mSFRjfW+eJj7v1AaISu7X+fGrXuZMuwL5jPnb8XPhfgaQ7DvM/WmnjzvdnAgQO1bds2LV68WL/5zW905swZ9e/fXz179lTHjh115ZVXqlGjRmrYsKGCgoJkGIZOnjypEydO6PDhw/riiy+0ZcsWrVu3TmFhYZo8ebIWLlyokBBLFz6tlsOwaMr2yJEj9c4776i0tFRXX321+vbtq+TkZMXHx+uLL77Qd999pwULFlhxaK8QFBxi20RoEgLfx4UeNSk5dkjnS4q18/U/6UReTo31g0LD1efBtxQT18ED0cGf8HmE2tj0+C900/+8Y3q7B//xS2VnZ5verqds27ZN7777rj755BN9/vnnOnr0qHPlpwYNGujMmYrVpwzDUJMmTXT99derT58+GjZsmK6//no7Q3eyLF1ZuXKlSkpKtHnzZm3atEkbNmzQsmXLdObMGTkcDv385z/X0qVLNXDgQOcMeNSMhMD3cQH2fd4wNOTAxn/oxMFdimrRrsakITbpZvX4XcVKSYY8//fnDe8X6q6q/z8+x1CdKneqD3A33HCDyx4OJ0+e1Lfffqvi4mKdPXtW4eHhio6OVlxcnBo2bGhjpJdnaR9HRESEUlNTlZqaKkk6d+6csrKytHHjRm3cuFH33XeffvjhB/Xv318ffvihlaH4LJIE3+IvF1Ju8rxf+1um6vDOj7RrxVz1nJSmkuMF+mLVYyo9XTFcqWHrzorv9Ss5gsMU1+XntsbqL+eFp/jC+cfQKFTHF/6G7dawYUOvTQ4ux7LhSe4oKytTdna2vvrqK59aTckdjqBgtb5xlGITBqhF0oAq65AQ+DZvuyDyIY1T3+3X4e0faG/6Exr4SKYaNGxmd0iwmbd/Lnjb5yjq7vDODBXuylBx3m4NmPKe6e3vf2O4Tw9P8ge2Jg3+7OI5DSQIvsXTFzJvv7DDt3w4o4cGLdhqdxjwY1Z+ZpFI+LZPHruNpMFP1XvJ1QULFtRrHdqjR4/qwQcfrG8YXomlSX2HWUtLsgEUvAN/U7CWlUs18/no+zyx5Co8r95zGh544AFNmDBBP//5z3XnnXcqODjY7Z9duXKl/vGPf/jUng3wPe5cvPhAAgBzmdVjwPwJ3xNk4pKr3+7+QId2f6igEydMaxN1U++ehqioKL322msqKipSt27dNGfOHH3wwQc6UcV/7unTp7Vx40Y9/PDD6tGjh3bt2qXXX3+d1ZNgKb7BAADfFSgbUKJqcdfdrO4jnlBMTIzdoQQ8U1ZPCg4O1gMPPKC7775br776qhYuXKgNGzbIMAw1atRIDodDx48fV1lZmZKTkzVkyBClp6erRYsWZhweAAAAXoIv4/yTqUuuRkdHa9KkSZo0aZLKyspUWFiow4cPq6ysTM2aNVOLFi3UoEEDMw8JAAAAL+HJHaHhWZbt07B9+3bdcMMNatmypVWHAAAAAOAB9Z7TcDm7du3S008/7VI2d+5cvfee+ctwAQAAwBs4FGSY/4D9LEsaxo4dqy1btmjBggXOsvnz5+u5557T66+/btVhvYZhGNqxbJYKd2bYHQoAAIClDu/M0PblMyVVDE8y+wH7WZY0LFq0SG3atNGgQYP07LPPOsvHjBmjhx56yKrDeg2Hw6FOox9V7GV2gwYAAPAXLZIGqPMdC2quCJ9lWdIwa9Ysde/eXV27dtV1112npUuXSpKKi4vVtGlTqw4LAAAAmzgMMTypjnbu3Gl3CNWyLGlITU3V/v37JUkDBw5UcHCw3nnnHX355ZdatWqVVYcFAACAjRiedHmjRo267Gvz58/3YCS1Z1nSkJaWphMnTqioqEiSNHr0aB07dkzffPONfvrTn1p1WAAAANgoqNxh+sNffPTRRzp69Ogl5StXrtTq1attiMh9DsMwPLqV4ltvvaW9e/fqj3/8oycP63EOh0Nh0U3sDgNAADp/qkihkY3sDgNAAIqObqVb71lners73r9Z2dnZprfraWFhYerfv79efvllxcfH69y5c3rwwQf16quv6ty5czp37pzdIV6WZfs0XM5tt92mwsJCTx/W4xxBwRq4IMvuMAAEoA9n9NSgBVvsDgNAAPr3ghFymPh19IEv1+jAf9aq/NQJ8xq10YYNG3TDDTdo7ty5atOmjV588UVFRERo+/bt2rVrl93hVcvjSYMkxcbG2nFYAAAAWMghc4cTtWk3WG3aDdbnHww0rU07devWTaGhoWrXrp0mT56s6667Ths2bFB4eLjatm1rd3jVsiVpAAAAAAJNUlKSOnfurA0bNmjVqlVq3bq1HnzwQU2dOlXr1q3T+PHj7Q7xskgaLFK5uVtsYn/2agAAAH7t8M4MFe6q2NDWn1Y7MtuXX36puLg4ff7557ryyislSU899ZTGjh2rt956i6QhEFVu7gYAAODvWiQNUIukAfr00V/61WpHZktMTNSHH36ooKAfFzCNiIjQ0qVLvX54kmVLrgIAAAD40T333OOSMFRq0KCB7r77bhsich89DQAAADCFQwqYHZxr69SpU5oyZYrz+YsvvijDMNShQwd16NBBjzzyiI3R1YykAQAgSSo3c53EWuImA/AX/rWDc32NGDFChw4dUmJiovr166exY8c6X7vvvvu0f/9+vffeexozZox++OEH56bI3oikAYDP85cbTjtv2u1GwgL4CUMKKrc7CO+Rnp6uTz/9VN27d6/y9Z/+9Ke6//77lZqaqsTERA9HVzskDR5i1kUpkG8q4Hu4Gaudi98vznfPsPJ95hyAt+LzxTO6det22YThQu3bt1ePHj08EFHdMRHaQkGGw/kAfMWFf7f1fQCBrtxhmPIAfIVDUlC5w/SHr4qPj3f++8UXX1RqaqqCg4PVr18/zZo167J1vRE9DT4myHBwAQlQ3IQDgYveEJjJ0vsIQ3IwPMkpMjLS+e/77rtPEyZMUEJCgtavX3/JKkoNGjTwdHi1QtJgEcMwtH35TMUmVKxbbCZ3PuBJLLwXF2gA3qQu1ws+x7yXXdf/wp0ZKsxdZ3q7X+9/T/v3v68z506Y3rYdgoOD1a1btyqXXfV2JA0WcTgc6nzHAtuOz9hoe3AhBRAI3Lmm8HnoGd5yfY9NGqDYpAHa+vAoU4cTtYsfqnbxQ7Uxq59pbXrSsWPHLikLCwursu6RI0esDqdeSBqAKnCxA4D6qevNLJ+/vs0hhiddKD09XcHBwZeUv/rqqzZEUz8BkzQUFhbqL3/5i6688kr9/ve/v+T1FStW6KOPPlK7du20detWjR8/XoMGDap1HW9Fz4MrLkoA4J3oxaheoF+/fU2DBg00atSoGocjlZeX65///KeHoqqbgEga3n//fb3++utaunSp5s6de8nrixcv1rx587Rnzx5FRkYqPz9fCQkJWr16tVJSUtyu40vc/cD19g+nQL5wAECg8sd5GN5+vXUb+zS4SE1NdbtXwZs3dpMCZMnVwYMHa968eVW+VlxcrKlTp2rs2LHOGe6tWrXS4MGDNXnyZLfrVIUl9NzHkp4AACux3G39uXtdZsnVH/3ud79zu+748eMtjKT+AiJpkFTleDJJWrNmjYqKipScnOxSnpycrJycHOXm5rpVx19ZdSPPzT8AwBfUJYkI5OTE8d8lV81++KpbbrnFkrp2CJik4XJycnIkSa1bt3Ypr9xgY+vWrW7VqQt//cCghwDwPf70GQRYKZASAOBCATGnoTqVy1tFR0e7lEdFRUmqmEDtTh0AAAB/Utcv+ZjT4J98LmnIz89XaWlpjfXi4uIUGhpaY73w8HBJFfsqXKhylntYWJhbdS5mGIY+efwXzuete9+h1r3vqDEeAKjEDvAA7HTx58/lkoi8zcuU9+/lkqQIhfr0cCJcns8lDSkpKcrLy6ux3u7du9WhQ4ca68XGxkqSTp065VJe+TwuLk4lJSU11rmEUa6S4wXOp3vfXai97y6sMR4AcFGHL/rOny7S2pk9TG8XAGoSEda65krwST6XNOzdu1eGUfM3b5W9AzXp2rWrJOngwYPq1KmTs/zAgQOSpC5duigmJqbGOhdzBAVr0IItbsVwIXezegCBoS49DR/O6Fnj5w89GAAuZNZ+TtvmjlJQOZ8v/sjnkgZ3kwF3paamqkmTJsrKytLQoUOd5ZmZmerevbvat2+vtm3b1lgHAKzAxowAfI2Zw5O+zH9PX377noodJ8xrFHUSMKsnnT17VlLFjnsXCgkJ0ezZs5WWluYchlRQUKD09HTNmTPH7TpWYUUGAAAQqK5tNUTDej7vHPUB+/hcT0NdZGRk6IUXXpAkvfnmm7ruuus0dOhQNWzYUJI0ZcoUhYeHa8KECUpKSlJ2drbS0tI0bNgwZxvu1LECw5MAAIDPMMTwJD/lMNyZIIBaCwoO0bCnv7Q7DAB+xp0eR+Y0ALDLtjmjdO/Ajaa3+85XNyk7O9v0duG+gBmeBAAAAKBuAmJ4EgAAAKznkORgeJJfImmwiGEY2r58pmITBqhF0gC7wwEAALBM4c4MFeaukwyDpMFPkTRYxOFwKGn0o5KkchlMaAYAAH6n8v7mysRUXZmYqs/m/JKJ0H6KpMFDqpp0SCIBAAB8BfctgY2kAQAAAOYwJEe5ibu7wWuQNAAAAMAUDhlysJq/XyJpAAAfx5ABAIDVSBoAwIeQIADwdmYOT9pduFa7C9fqRPgJ09pE3ZA0AAAAwByGufs0dGx2szo2u1lvfjfYtDZRNyQNABAAqlrBDQAAd5E0WMQwDO1YNkuxif0VmzSAIQUAbEPCAMAM1S0ff3hnhgp3ZUgyJIPVk/wRSYNFHA6HOv13czcAAAB/VJlINO/UX8079VfOrBEsueqnSBosRO8CAAAIKIbkoKfBL5E0AIC/cTAkCQBgLpIGAAAAmKJic7cyu8OABUgaAAAAYB7mNPglkgaLGOVlWjuzh91hAAhA508X6cOZPe0OA0AAaq6WdocAi5A0WMQRFKxBC7bYHQYAL+GwcGEE46L5C2tn9dDABVmXfd1MVv5eAHzP9pm/YCK0nyJpAACTmXkj7Sivvi0jyHDreEZQ9YlDTce5bLtVJCQkEkAAMySZOKdh15F1+uLoOp2IPGFam6ibILsD8FeGYWj78pk6vDPD7lAA+Ch3buQvV8dwGD8+akgYpJqTCgCoTuHODO1YNsv0dhOa9teo9o8oJibG9LZRO/Q0WMThcKjzHQvsDgOADzOCDLd6Gqri8m1/uXU9DQAgSbFJAxSbNEDbdwxnR2g/RdIAACYzc8iOW70EbsxZsCopYCgSgIux5Kp/ImkAAA+wauy/WwmD4TBlMjQJAgAELpIGALCJlasaXYwbfgCeYZg6ERreg6QBAAAA5mFOg18iaQAAAIBJDEn0NPgjllwFAAAAUC16GgAAAGAOw2D1JD9F0mCRys3dYhMGqEXSALvDAQAAsEzhzgwV5q6reELS4JdIGizC5m4AACBQODd32z7M7lBgEZIGAAAAmMTcidC7jn2iXcc360T0CdPaRN34/UTo8+fP65FHHtE111yjqKgodevWTe++++4l9VasWKEJEybo8ccf16hRo7R27do61amvcofh8gAAALBakOFwedSLUWbaI+EnN+pXV09VTEyMOb8o6szvexqmTZumoKAgPfbYY/rmm2/0+OOPa/jw4Vq/fr369u0rSVq8eLHmzZunPXv2KDIyUvn5+UpISNDq1auVkpLidh0AAADAH/l1T8OhQ4cUERGhp556SiNGjND//M//6J133lF5ebleeeUVSVJxcbGmTp2qsWPHKjIyUpLUqlUrDR48WJMnT3a7jj+4uJeDng4AADzj4m/66/1tv20qhyeZ/YDd/DppKCws1LRp01zKunfvrsaNG+v48eOSpDVr1qioqEjJycku9ZKTk5WTk6Pc3Fy36vgidxIEkggAAMznToLgm0mEIcMoM/0B+/n18KTrr7++yvIzZ86od+/ekqScnBxJUuvWrV3qxMfHS5K2bt2qffv21VgnMTHRvMA9gAQAAADfEmQ4PHb9rt+xuMn3R36dNFQlMzNTkZGRGj9+vCTpyJEjkqTo6GiXelFRUZIqeivcqWMGKz8I6tq273yzAQCA77j4uuzu9fbienwJCE/xuaQhPz9fpaWlNdaLi4tTaGjoJeULFizQX//6VzVu3FiSFB4eLqliX4ULBQVVjNwKCwtzq87FDMPQxieGO5/H97pD8b1H1xg3AACAN3A3QTmwebkObF4uSYqWIRnllscGz/O5pCElJUV5eXk11tu9e7c6dOjgUrZo0SINGDBAw4f/eDMfGxsrSTp16pRL3crncXFxKikpqbHOxRwOh1KmrXIpK1f13wbwrT4AALBDVQnBxfcll7tPadNrtNr0qvhi9PMZg8TwJP/kc0nD3r17ZRg1d8VV9g5UWrVqlU6ePKkZM2a4lHft2lWSdPDgQXXq1MlZfuDAAUlSly5dnGsDV1fnYkZ5mT6c2dOdX+mCH6pd9WqZlX/Q6wn4nPOni2r/+QPAXmZdb838/rEOMTXXT0wMAN7E55KGi5MBd3zwwQfasWOH5syZ41K+Y8cOpaamqkmTJsrKytLQoUOdr2VmZqp79+5q37692rZtW2OdizmCgjVwQVat4jSrp8HM8Y30fgC+Z+2sHrX+/AFgLzOvt2bdB9Qlps//eHPFxmzwO3695Kokbdq0SQ8//LA6duyolStXauXKlXrzzTf14IMP6tChQwoJCdHs2bOVlpbmHIZUUFCg9PR0Z5LhTh07sSwqAACoZP9Srebtz7DreKbe3P+8Tpw44elfAhfxuZ6G2sjJydHQoUNVXFyszZs3u7wWHx+vp556SpI0ZcoUhYeHa8KECUpKSlJ2drbS0tI0bNgwZ3136tSWt3+L7+3xAQDgL8odhtddd+u6wpOZEn7SXQk/6a7XyxZ4/Nhw5ddJQ5cuXXTy5Em36k6cOFETJ06sdx1PsKo3wds+rAAACCRW3aR7dplWg83Y/JRfJw3+iOFHAACgPqzdJM4Qqyf5J5IGC5nxDQFJAgAAgcnK4UF16X0ItHuSJUuWKDg4WPv27VObNm00btw4u0OyFUmDRQzD0PblMxWbMEAtkga4/XOePCEZjgQAAKT6D2Eq3Jmhwtx1/33m+z0NmZmZSk9P11tvvSVJSkxMVOfOnZ1L9QcikgaLOBwOdb7DeyftkDAAAOBbvHGydKXYpAGKTRqg7Tmpps5p+KLoc31R9LlONPLs6kmrVq1S27Ztnc9TUlL0yiuvBHTS4PdLrgIAAMBTKuc0mPPo2KiTRra5y7nRrqeUlJS4LPPapEkT7d+/36MxeBuSBi9i5dAk+9dsBgAA9eWpvZl8/V6hsLBQ06dP19NPP13l6ytWrNCECRP0+OOPa9SoUVq7dq3L66mpqVq9erWOHTsmSTp48KAiIiIsj9ubMTzJRoE2oQgAAPiOOs9zMMotiMZ977//vl5//XUtXbpUc+fOveT1xYsXa968edqzZ48iIyOVn5+vhIQErV69WikpKZKkW265RQcOHND999+vhIQE7d69WyNHjvT0r+JV6Gmw0MXfBnhy52Z6FgAA8H923ltU9TBkyFCZ6Y/aGDx4sObNm1fla8XFxZo6darGjh2ryMhISVKrVq00ePBgTZ482aXu/fffr2XLlmn06NH66quvdPfdd9fpffMXJA0AAADwK8HBwVWWr1mzRkVFRUpOTnYpT05OVk5OjnJzc13Ki4qK9Nvf/lZpaWlq3ry5ZfH6AoYn+SF6FgAACEzescKS9y65mpOTI0lq3bq1S3l8fLwkaevWrUpMTNSJEyeUnp6urKwsPffcc+rSpYvHY/U2JA0+yP4PAwAA4K2s3BSuZoZk0pKrmUc26+PDH+lMWYkiG19hSptHjhyRJEVHR7uUR0VFSaqYQC1JYWFhGj16tMaOHWvKcf0BSYNFDMPQjmWzFJvYX7G12NytJiQMAACgNqqa62D2/cThnRkq3JVhapvJTXsruWlvSdIbjr+b0mZ4eLikiv20LhQUVDFiPywsTJICKQ7qVAAAFvdJREFUfqWkqjCnwSIOh0OdRj9a74SBCc0AAMBsZk+gbpE0wLmprScmQm/atElRUVEujwMHDtQYZ2xsrCTp1KlTLuWVz+Pi4ur7VvgtehoAAABgEkOS9UuuduvWzTk/oVLLli1r/LnKHZ0PHjyoTp06OcsrEw7mLlweSYOF6BkAAAC+wLTlWg2ZNqfBxUW3VBEREWrXrl2tm0lNTVWTJk2UlZWloUOHOsszMzPVvXt3tW/fvr6R+i2GJwEAAMCvnD17VpJUXu7a6xESEqLZs2crLS1NJSUlkqSCggKlp6drzpw5Ho/Tl9DTAAAAAJMYtd6MzWwZGRl64YUXJElvvvmmrrvuOg0dOlQNGzaUJE2ZMkXh4eGaMGGCkpKSlJ2drbS0NA0bNszOsL0eSQMAAADMY1g/p6E6AwYM0IAB1S9EM3HiRE2cONFDEfkHkgaLGOVlWjuzh91hAAhA508X6cOZPe0OA0AAaqYGdocAi5A0WMQRFKxBC7bYHQaAALR2Vg8NXJBldxgAAlDOtL62D0+CNUgaLGIYhrYvn6nYhAFqYeLmbgAAAN6mcGeGCnPXycwdoeFdSBos4nA4nJucAAAA+LPYpAGKTRqgnM/72B0KLELSAAAAABPZOxEa1iBpAAAAgEnsX3IV1iBpAAAAgHmY0+CX2BEaAAAAQLXoaQAAAIApDBkymNPgl0gaAAAAYB6GJ/klhicBAAAAqBY9DRZhczcAABAoLtzcjeFJ/omkwSJs7gYAAAJF5eZun2/rxfAkP8XwJAAAAADV8vukwTAMPfnkk2rbtq0iIyOVkpKibdu2XVJvxYoVmjBhgh5//HGNGjVKa9eurVMdAACAwGWoYkdosx+wm98PT3r88ccVGRmpzMxMHT16VPfcc49uvfVWHTx4UA6HQ5K0ePFizZs3T3v27FFkZKTy8/OVkJCg1atXKyUlxe06AAAAgc5geJJf8uuehpKSEl1zzTWaNGmSmjVrpg4dOmj69On69ttvdfToUUlScXGxpk6dqrFjxyoyMlKS1KpVKw0ePFiTJ092uw4AAADgr/w6aYiIiNCIESNcyo4fP64ePXqoadOmkqQ1a9aoqKhIycnJLvWSk5OVk5Oj3Nxct+oAAACA4Un+yu+HJ12osLBQ//rXv/TWW285y3JyciRJrVu3dqkbHx8vSdq6dav27dtXY53ExETL4gYAAPBGQYbjkjKGJ/mngEgaioqKtHDhQj3//PM6duyYIiMjlZaWppCQEB05ckSSFB0d7fIzUVFRkioSDXfqAAAA+LOqEoRLGZJBz4A/8rmkIT8/X6WlpTXWi4uLU2hoqCQpJiZG06ZN04gRI/Tkk09q6dKl6tWrlyZNmqTw8HBJck6KrhQUVDFyKywszK06FzMMQxufGO58Ht/rDsX3Hu3urwkAAOAx7iUEVcvbvEx5/14uSYo0KyB4HZ9LGlJSUpSXl1djvd27d6tDhw6SKm72IyMj1blzZ7322mv6+uuvtW7dOk2aNEmxsbGSpFOnTrn8fOXzuLg4lZSU1FjnYg6HQ30feruWvx0AAIC16pMgVCW+92jnF6OfTeshQwxP8kc+lzTs3btXhmHUWK+yd6AqvXv31oEDByRJXbt2lSQdPHhQnTp1ctapfL1Lly6KiYmpsQ4AAIDdzE4Iao/hSf7K55KG6pIBd+Xn52vgwIGSpNTUVDVp0kRZWVkaOnSos05mZqa6d++u9u3bq23btjXWuZhRXqa1M3vUO1YAqK3zp4v04cyedocBwGw1f2dqu6Z2BwDL+FzSUBtFRUV6+OGHNXLkSPXq1UuStGXLFh0/fly//e1vJUkhISGaPXu2Fi5cqBkzZigiIkIFBQVKT0/XsmXL3K5zMUdQsAYt2OKR3xMALrR2Vg8NXJBldxgAasH+HgJzZE/rzvAkP+XXScO5c+f08ccf6/nnn9fAgQPVp08fNW3aVG+//baCg4Od9aZMmaLw8HBNmDBBSUlJys7OVlpamoYNG1arOgAAAFXxl6SgRgbDk/yVw3BnggBqLSg4RMOe/tLuMAAEIHoaAM8JmGTATdkPddO9bQaZ3u47Uf9Rdna26e3CfX7d0wAAAFAfJAW1x+Zu/omkAQAA+D1u/j3FkMTwJH9E0mARwzC0fflMxSYMUIukAXaHAwCAT+Jm3zcc3pmhwl0ZFU/oafBLJA0WcTgc6nzHArvDABAAyh1MTYP34+bfv7VIqviSNPuztXaHAouQNACAj+NmLDDx/w7vZMgwcXjS3lOH9Z/ThTpphJrWJuqGpAEAgDrgph24DBOXXG0f2VztI5trdUyBaW2ibkgaAMDfGJ67obVyaBQ35QDgPUgaAAB1xo09AFcGS676KZIGAAAAmMMQO0L7qSC7AwAAAADg3ehpAAAAgEkMGWJ4kj8iabAIm7sBAIBA4bq5G8OT/BFJg0XY3A0AAAQK5+Zu2e/ZHQosQtIAAAAAUxisnuS3SBoAAABgIoYn+SOSBgAAAJiEngZ/xZKrAAAAAKpFTwMAAADMw+pJfomkAQAAAOYwDBkkDX6JpMEiRnmZ1s7sYXcYAALQ+dNFfP4AsEVjuwOAZUgarOIIUotOqWzuBsDj1s7soUELttgdBoAAUrm527mt70hMhPZLJA0WYXM3AAAQKCo3d9u6NV2GJUuusnaP3fgfAAAAAFAtehoAAABgHkuGJ/E9t91IGgAAAGASVk/yVyQNAAAAMA8Tof0SfT0AAAAAqkVPAwAAAExhsLmb3yJpAAAAgHlMHJ607+xp7Tt7Wj9ENTetTdQNw5MsYhiGti+fqcM7M+wOBQAAwFKHd2Zo+/KZprfbLjxSP2/YXDExMaa3jdqhp8EibO4GAAACReXmblu2vMXwJD9F0gAAAAATsXqSPwq4pGHu3LnKy8vT3//+d5fyFStW6KOPPlK7du20detWjR8/XoMGDap1HQAAAMDfBFTSsGXLFv35z3/W6NGjXcoXL16sefPmac+ePYqMjFR+fr4SEhK0evVqpaSkuF0HAAAgsLF6kr8KmInQP/zwg55//nn16NHDpby4uFhTp07V2LFjFRkZKUlq1aqVBg8erMmTJ7tdBwAAAJJhlJn+gP0CJmmYM2eOZs+ereDgYJfyNWvWqKioSMnJyS7lycnJysnJUW5urlt1AAAAAH8VEEnD+++/r7Zt2+qaa6655LWcnBxJUuvWrV3K4+PjJUlbt251qw4AAEDAMwzJKDf/Adv5/ZyGo0ePasWKFVqyZEmVrx85ckSSFB0d7VIeFRUlSSosLHSrDgAAAMRwIj/lc0lDfn6+SktLa6wXFxen0NBQzZgxQ48++uhl64WHh0uq2FfhQkFBFZ0wYWFhbtW5mGEY2vjEcOfz+F53KL736EvqAQAA+LK8zcuU9+/lkqQGMkzdERrew+eShpSUFOXl5dVYb/fu3crKylL//v3VsmXLy9aLjY2VJJ06dcqlvPJ5XFycSkpKaqxzMYfDob4PvV1jnP4ob/MyEiQfx/+hbys7f8buEFBPnIO+LdD+/+J7j3b+vlsevNbmaGAVn5vTsHfvXpWUlNT4aN++vf7+979r/PjxioqKcj42bdqkf/zjH4qKitKCBQvUtWtXSdLBgwddjnPgwAFJUpcuXdyqY5XDOzN8ru3KbxusYFXMvvg+WxmzVf+HvM+eabvcwqSB99n6diXOQU+0a2XbgXwdNCQZ5eWmP2A/n0sawsPD1aBBgxofDodDS5cuVU5OjsujW7duuvXWW5WTk6OJEycqNTVVTZo0UVZWlstxMjMz1b17d7Vv396tOlYp3GXdh6WVbVvFqph98X3m/88zbftizFbifba+XSvxPnuubat4/fts/Hd4ktkP2M5hGIZhdxCedNNNN6lNmzYuO0I/88wzWrhwoXbv3q2IiAgVFBTo2muv1bJlyzRs2DC361yoadOmatOmTb3jPXHihGJiYurdjifb/v7779WsWTPT25Wsi9kX32crY7bq/5D32TNtcw56pm3OQc+07YsxB/I5uGPHDnXq1MmEiFw1bdpUa9asMb1duM/n5jRYYcqUKQoPD9eECROUlJSk7OxspaWluSQD7tS5UOWKSwAAAICvC7ieBgAAAAC143NzGgAAAAB4FsOTUG+FhYX6y1/+oiuvvFK///3vL3l9xYoV+uijj9SuXTtt3bpV48eP16BBg2pdB/bZu3evpRP+YQ7OI9/HuebduN4hoBlAPbz33nvGmDFjDEnG3LlzL3n9b3/7m9GqVSvj1KlThmEYxsGDB42GDRsamzZtqlUdeM7evXuNoKAgQ/9dOU+SsWDBApc6b7zxhvG73/3OeOyxx4yRI0caa9assSlaVOI88j2ca76F6x0CHUkD6m3fvn1VfoiePHnSaNSokTFjxgyX8ttvv93o0qWL23XgWb/73e+M5557zliyZImxZMkS49VXXzVOnjzpfJ2LnvfhPPJNnGu+h+sdAhnDk1BvwcHBVZavWbNGRUVFSk5OdilPTk7WG2+8odzcXO3evbvGOomJiZbFDleHDh3S0aNHNWnSpCpfLy4u1tSpU3XfffcpMjJSktSqVSsNHjxYkydP1ueff+7JcPFf7pxrnEfehXPNN3G9QyBjIjQsk5OTI0lq3bq1S3l8fLwkaevWrW7Vgef85S9/0b/+9S9dddVVuueeey65ManuwpiTk6Pc3FxPhov/4jzyPZxr/oXrHQIBSQMsU7lXRXR0tEt5VFSUpIoJZe7Uged07txZ06ZNU8uWLfXKK6+oR48e+tvf/uZ8nYued+I88j2ca/6F6x0CAcOT4JSfn6/S0tIa68XFxSk0NLTGeuHh4ZIkh8PhUh4UVJGrhoWFuVUHnjN27Fjnvz/++GONHTtW9957r3r16qWOHTty0fNSnEe+h3PNv3C9QyAgaYBTSkqK8vLyaqy3e/dudejQocZ6sbGxkqRTp065lFc+j4uLU0lJSY11UHf1SQRvuukmffDBB+rcubPeeOMNzZ8/n4uel3LnXIP34lzzfVzvEAhIGuC0d+9eGW5sEF55MatJ165dJUkHDx5Up06dnOUHDhyQJHXp0kUxMTE11kHd1TcRvO6669SnTx8dPnxYEjen3sqdcw3ejXPNt3G9QyBgTgOcwsPD1aBBgxofF3/zdTmpqalq0qSJsrKyXMozMzPVvXt3tW/f3q06qLu9e/eqpKSkxkd173Pjxo2dNzAXXhgvxEXPXpxH/oFzzXdxvUMgIGlAvZ09e1aSVF5e7lIeEhKi2bNnKy0tzdktW1BQoPT0dM2ZM8ftOqi7+iaCpaWl2rZtm4YPHy6Jm1NvxXnk+zjXfAPXOwSy4Hnz5s2zOwj4royMDP35z3/W7t279d1336lJkyaKj493DmFKTk5WgwYN9MILL2jfvn165ZVXNGfOHOeF0d06sN7HH3+sBx98UNHR0br22mslSXPmzFGnTp105513SqoYTx0eHq6XX35Z9957r0JDQ1VQUKBJkybpueeec/4cPI/zyHdwrvkmrncIdA7DnUHsAPxeTk6O7r77bn3xxRcaMmSIrrnmGt10000aMmTIJXVfeuklbd68WUlJScrOztbo0aN122232RA14Hs41wD4IpIGAAAAANViTgMAAACAapE0AAAAAKgWSQMAAACAapE0AAAAAKgWSQMAAACAapE0AAAAAKgWSQMAAACAapE0AAAAAKgWSQMAAACAapE0AICPKC4u1qOPPmp3GACAAETSAAA+oLS0VBMmTNC4ceM8fuw//OEPKiws9PhxAQDeg6QBAHzAokWLNGjQIF155ZUeP/aMGTM0YcIEnTlzxuPHBgB4B4dhGIbdQQAALu/kyZPq06ePsrOzFRoaaksM//znP/Xpp5/qySeftOX4AAB70dMAAF7u1Vdf1S233GJbwiBJI0aM0Lp165Sfn29bDAAA+5A0AICXW716tVJSUmyNweFwaPjw4XrllVdsjQMAYA+GJwGATV5++WUdOXJEu3fv1l133aW8vDwdPnxYOTk5Wrhwoa666iqVlpYqOjpahw4d0k9+8pMq2/nlL3+pYcOG6e6773aWffjhh7r66qvVtm1bSdKzzz6r3NxchYWFaf78+Vq8eLEcDoc+/fRTDR06VL/4xS/00ksvyeFwKDs7WyNHjtSvf/1rl+OsX79ef/jDH7Rt2zbr3hQAgFcKsTsAAAhES5YsUadOnXTjjTfqk08+0dChQ7Vo0SL16NFDs2fP1j333KOrrrpKhYWFKi8vv2zC8MMPP+idd97RL3/5S2dZeXm5br/9dr300ktq27at9uzZo8aNG2vy5Mnq1KmTQkNDtWDBAkVEROjtt9/WmDFjtG/fPs2fP19hYWFKT0/XuHHjLkkaunbtqp07d+rs2bMKDw+39P0BAHgXhicBgA2OHj2qG2+8UZJUUFCgoKAgjR49WikpKfr44481ePBgSdLhw4fVqFGjy7bz73//W+fPn1ffvn2dZTt27NDx48fVp08fSVJmZqaGDBmi7du3q1GjRvrTn/6kiIgIZxwOh0NTp05VWFiYpIr9IKrSsGFDBQUFaf/+/fV/AwAAPoWeBgCwwdSpU53/3rhxo/r27avg4GBdccUV6tevn/O1srIyORyOy7azceNGXX311WrVqpVLWbt27ZzLs/72t7+VVDG8qH///mrSpImzbmZmpvr16+dStn79epcYLtS4cWMdP368dr8sAMDn0dMAADZbv369brrppipfa9asWbU36Rs2bHDpZZB+TEIutm7dOvXv3/+SY//sZz9zPj937pzeeustjRkzpsrjlZaWMjQJAAIQSQMAeFhZWZk+/PBDlZWV6bvvvtOuXbtcvtl/7LHHnP9u0aKFysrKqhwydPbsWW3ZskXJycku5Z988on69OmjY8eOaefOnZKkb775Rt98841LgpCfn699+/a5lKWnpyskJETDhg3Tpk2btHnzZudrhmHo+PHjat68ef3fBACATyFpAAAP++tf/6qbb75ZeXl5euONNxQREaHY2FhJ0qpVq9SxY0dn3YiICCUnJ+uzzz67pJ2srCydOXNG8fHxzrJt27apsLBQvXv3VkZGhsrLyyVV9DK0aNHCpe1169apcePG6ty5s7Ns48aNuu222xQaGqq3335bPXv2dL527NgxhYaGOmMFAAQOkgYA8LA+ffpo9OjRWrp0qa666iq9/PLL+tOf/qT58+fr5MmTuuWWW1zqDx06VB9//PEl7WzcuFGS9MUXX0iSvvvuOy1dulQOh0MtW7bU9u3bnQnB119/rTvvvNPl57/55huNHj1aQUE/XgpGjx6tQ4cOafr06RozZoxCQn6c+lY5/8HOTeYAAPZgnwYA8HIFBQW65ZZbtGXLFpcb/JtvvlnBwcHq3LmzGjZsqJiYGE2cOFGvv/66du3apVtvvVW9evUyLY5p06apVatWeuCBB0xrEwDgG0gaAMAHzJ07V4mJiRo1apSkignJjRo10qJFi3TXXXdZfvxz586pR48e2rBhg2JiYiw/HgDAuzA8CQB8wPTp07V06VJ9//33kqTPPvtMp0+f1s033+yR4y9ZskR33XUXCQMABCiSBgDwAVdccYVeeuklPfLII5Kkffv2qX///s69GKy0b98+ffDBB5o0aZLlxwIAeCeGJwEAqjVnzhw99NBDio6OtjsUAIBNSBoAAAAAVIvhSQAAAACqRdIAAAAAoFokDQAAAACqRdIAAAAAoFokDQAAAACqRdIAAAAAoFokDQAAAACqRdIAAAAAoFokDQAAAACqRdIAAAAAoFokDQAAAACqRdIAAAAAoFokDQAAAACqRdIAAAAAoFokDQAAAACqRdIAAAAAoFokDQAAAACqRdIAAAAAoFr/HzQukPNJAP3yAAAAAElFTkSuQmCC\"><br>"
],
"text/plain": [
"<yt.visualization.plot_window.AxisAlignedSlicePlot at 0x1221068d0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"ds = yt.load( './diags/plotfiles/plt00200' ) # Create a dataset object\n",
"sl = yt.SlicePlot(ds, 2, 'Ex', aspect=.2) # Create a sliceplot object\n",
"sl.set_xlabel(r'$x (\\mu m)$') # Set labels x\n",
"sl.set_ylabel(r'$z (\\mu m)$') # Set labels y\n",
"sl.annotate_particles(width=(10.e-6, 'm'), p_size=2, ptype='ions', col='black') # Plot particles species=ions\n",
"sl.annotate_particles(width=(10.e-6, 'm'), p_size=2, ptype='electrons', col='black')\n",
"sl.annotate_particles(width=(10.e-6, 'm'), p_size=2, ptype='beam', col='black')\n",
"sl.annotate_grids() # Show grids\n",
"sl.show() # Show the plot\n",
"\n",
"#############################\n",
"### OTHER USEFUL COMMANDS ###\n",
"#############################\n",
"# # List all fields in the datasert\n",
"# ds.field_list\n",
"# # Get All Data from the dataset\n",
"# # Then get some data. \".v\" converts arrays from units-aware yt arrays to numpy arrays.\n",
"# ad = ds.all_data()\n",
"# Bx = ad['boxlib', 'Bx'].v\n",
"# # Get All Data from the dataset, on a given level and given dimension.\n",
"# # Then get some data. \".v\" converts arrays from units-aware yt arrays to numpy arrays.\n",
"# # This is similar to the 2 lines above, except that F has the proper shape.\n",
"# all_data_level_0 = ds.covering_grid(level=0,left_edge=ds.domain_left_edge, dims=ds.domain_dimensions)\n",
"# Bx = all_data_level_0['boxlib', 'Bx'].v.squeeze()\n",
"# # particle \n",
"# # CAREFUL! For the moment, 2d WarpX simulations use (x, z) spatial coordinate \n",
"# # but they are stored as (particle_position_x, particle_position_y) in Yt\n",
"# x = ad['beam', 'particle_position_x'].v\n",
"# z = ad['beam', 'particle_position_y'].v\n",
"# # For 2d simulations, WarpX and Yt use (ux, uz)\n",
"# # ux/c should be the nirmalized momentum\n",
"# ux = ad['beam', 'particle_momentum_x'].v\n",
"# uy = ad['beam', 'particle_momentum_y'].v\n",
"# uz = ad['beam', 'particle_momentum_z']\n",
"# w = ad['beam', 'particle_weight'].v\n",
"# # Set figure size\n",
"# sl.figure_size = (9, 7)\n",
"# # Save image\n",
"# sl.save('./toto.pdf')\n",
"# # This returns the domain boundaries\n",
"# sl.bounds"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Read data back-transformed to the lab frame when the simulation runs in the boosted frame (example: 2D run)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [],
"source": [
"# read_raw_data.py is located in warpx/Tools.\n",
"import os, glob\n",
"import read_raw_data"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [],
"source": [
"# For the moment, the back-transformed diagnostics must be read with \n",
"# custom functions like this one.\n",
"# It should be OpenPMD-compliant hdf5 files soon, making this part outdated.\n",
"def get_particle_field(snapshot, species, field):\n",
" fn = snapshot + '/' + species\n",
" files = glob.glob(os.path.join(fn, field + '_*'))\n",
" files.sort()\n",
" all_data = np.array([])\n",
" for f in files:\n",
" data = np.fromfile(f)\n",
" all_data = np.concatenate((all_data, data))\n",
" return all_data"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Available info: t_snapshot field_names xmin ymin zmin xmax ymax zmax nx ny nz x y z\n",
"Available fields: ['Ex', 'Ey', 'Ez', 'By']\n"
]
},
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x1202290f0>]"
]
},
"execution_count": 25,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAADNCAYAAACmVDdwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzsvXucZMlV3/mNzJtZWV3VPd3qHqmHeSCBBvBIuEfSrCSvwcvLaIQXjY0QFlob4c/w0YIl42WtXUvGBiwQSDaGDyzPAWaR/LEZBsnAYAvJWIBhbYE0EmqkkSwYRq+Wppnunq6equ7KqryVsX/cOPeeezLuzaxHd1d2x68++al740bEjUfmOXEeccJ570lISEhISGhC52o3ICEhISFhfyMxioSEhISEViRGkZCQkJDQisQoEhISEhJakRhFQkJCQkIrEqNISEhISGhFYhQJCQkJu4Rz7n7n3BPOuY/OkPcLnXPvdc79qXPu951zt1yJNu4GiVEkJCQk7B6/DNw9Y94fBd7uvf+rwJuAH7lcjdorJEaRkJCQsEt47/8AeFKnOee+2Dn3bufcB51zf+ic+7Lw6A7gveH694B7rmBTd4TEKBISEhIuD+4D/pH3/gXA64GfCekngZeH678DHHTOHb0K7ZsZ2dVuQEJCQsK1BufcMvA/A7/mnJPkhfD/9cBPOee+HfgD4HNAfqXbuB0kRpGQkJCw9+gAK977O+0D7/3ngW+CkqG83Ht/4Qq3b1tIqqeEhISEPYb3/ingk865VwC4AifC9THnnNDeNwL3X6VmzozEKBISEhJ2CefcrwDvA77UOXfKOXcv8L8B9zrnTgKPUBmtvwr4hHPuz4BnAG++Ck3eFlwKM56QkJCQ0IYkUSQkJCQktCIxioSEhISEVlzXXk/Hjh3zz3zmMwHwvvgAiDebpOn0Tic8t5lRabqQzuNcPb/JarWA0WJN72goZPtgmy11l32yBXS9ofNN2aSeTkfVNR5Pti3UMx4Xj6WOTic8QpX1vqzUuw55XiVlWcgriVnGlu8wHkO3Cx3GVcUJCbvEBz/4wbPe+xt3U8eznfOXWp4/Du/x3s+6w/uK4bpmFM985jN5+P3vZ0xBgPLgyZyFUZE0nT4YQD8bT2aWArqgQPJkWXk9DsKczpobT2pbrMO4+R2RQuOsP9EH22ypO8ugk29OdlrXOxhAlrGZd6LZpJ7BQNU1HE62bTBgnPUZDovHo1GRvLgYxpdQdm2t+J9lsLzMZnaA8+dhY6NIOnYs5D17tsh37BhP5QdYX4eDB+EAl6r3KeG5ZCAJCduA63Y/vds61oHXtjz/53Bst++4HLiuGQVQMgkNTdtitHhMh06WTVLKNiKuuU+4tgQ89i6NvsxWjPprhpVNTqtkFaI8GkGvF8nexCgibV9fL5K3tooVfDS/5bQyboGJDYf1Vw0GVAxGHpZMqmASa2shn7xDGMrhw6yvF3kOHoyMS0LCVYZjPvX91/2vyK6I7cpbCJKmN8MhZFmHLOtXK2eNNuLU8iz2SK/4ow9biGGxcu7UsgpzaG2E1NnCubKsqGs0mmQSeQ79Ka+ZGYY5RZ9Lg9Q4DAbA8Lr/eifsQ8zjt3Ie23zZ0KhGEZVPULtAwSwKNVREspCCuuIIMZ9lsTuRxyY0cZfI4xjtrxVXqqFGBpTnZFm/zNbIT/QAyviY+uSxSHAlE9PcUV4SykrS0lJQAZp3ba3A8nJQfUkjExL2CZJEMafQdN4uTmu2CCiJZF3jE9RQkqAr0C+Z8v4p9B4IKi/GdRVOW8UUTK6QfqZrYkqVWky0UpA6bR/q/enQ0VygpV0ikTQyxcAs8mF1u7hInQEtL5dMvJAmhjXbRLJLJOwHOOaT6M5jmxMSEhLmFkmimGPoxWkpTUQ8djqMGQzqXj9Z1q9PfqhIDOWyMNar2g7jcrW7HXur9t4hqICmrZbleT8z5Yl5Q4UVuG2UydjJsppUIY+0eaMmbek6DMTWEYUx0oj0URq9VZ6hkjhYyctXNkpmCQlXGEmimFPEVE5ZFnmg/lt1TsksrMomEPIJIhyI1HaIlSXw5SsixL+sN6Y2UlSziWnU2hpTcRlm0fC4GC5rw6n9L8ZHM4ksA4aT7Ramu7VVmS0Y1rm7qKZmGdfELBKuBhKjmFNoGiYEqCQgsVV1IHyd8DzLOsqLM3hCBeO3MJQYz5motwktrq72ulEyidVvMjeZJPpZJFOtnn7UpCFpE+VnFZ+MiNfq8RQwIZUoN9z+df9NT9gvSKqnhISEhIRGzKtEMY/Mbc+gw0+Iyqm2+1nromIb0UqpobxlOKTwvmmQBPJcqXyaNrbZQpFb3ZwWB6VGtVGsXbG6JlRe1lah1DdNks5ExeqhqJOmtRFa7BhGUomp6RIS9gs6LZ/9inlkbnsGYRQ1lVOTLkfSi912E3sVJojsoFPWKfr+ugqqU30xZlANacSyT2h1mjhIzZe1vi2uVY1l6wsZbPgTIeZlecs9VMV5XuUXZhG1ryjVU7erqtAcOs/p9SqbkOxvyQbVBslkk0i42kgSxZxCYjfVdljbfRCSJohIAvaxRKDQnk0T1VrlfuOSfDJLTMCJolXcaOxOiQni2pK5dcU/Qztq+QVqwHo9s7M8soGvLK4cDywkLUkeCVcasuFu3iSKPWmbc+5u59wnnHOPOufeEHm+4Jz71fD8j51zz1TP3hjSP+Gce8m0Op1zrwtp3jl3TKU759xPhmd/6px7/vR2KymiiYDNRI0LaLq1vl4xCxDJomFTmWUQU94XdSCaNWODKgsqQj9RX4PaLVaPViOV3mMWkYGoSQq1CuJFhkOijKLXE4luUDZMhngnjGFM56p95P0J1xayls9+xa6/hc65LvDTwEuBO4Bvdc7dYbLdC5z33j8b+HHgraHsHcArgecAdwM/45zrTqnzvwFfB9hIji8Fbg+f1wA/u9u+JSQkJOwlditROOfud8494Zz7aMPzbS+YZ8FeLFdeCDzqvX/Me78JPEB1NqzgHuBt4fodwNc651xIf8B7v+G9/yTwaKivsU7v/Z947z8Vacc9wNt9gT8CDjvnbmpruHO0q0ViK/GoWFBBPxqN6iooWSuWMYp0oRl0SbGmzixNtDRU2woa1Ue23kibJozSsTLmvdG8Wk8X8sf8Cmr15Xkt9pR2GJjVI7dpVX+1oKWKJHFcGxAbxS4kil+mWFQ34bIsmPfiW3Yz8Fl1fyqkRfN473PgAnC0pewsde6kHXV4HydibQQ7amyYtPOKLn00MtWFmzJonX23rizAhkLftl2ioT8z2zhiTMz0P8ZgJpwDVFlN4LrdYrzKDXqRMS52xBdJMqaxjZC62Vqt19TXeSa4MUai0xP2J3bDKLz3fwA82ZJl2wvmWbAXajEXSbNHrzXlaUqPfdMjx7ltux04515DwWm57bbbJkvMaI+I6cfFHhETTAoG0qFjpZMmGA8hCyHMNmz4hNeQbUykTu15JLaCaod6cxMLQtxshogmmgfRMzFiGx3zvAgySIetrUik2ZAnU4c1iU+XDrcew7VEWJuYhYSMSZ5fVxczRI895px7WN3f572/bxuvaFowP76NOiawF4ziFHCrur8F+HxDnlPOuQy4gYIrtpWdVudO2kEY9PsA7rrrLh/dVtxEXA0Bsz9Ga5fWRLwkXlbNZFfPDe/WRbarHppAVsWhmqo2mlKffSTSQeyddhz1cMr1RARbo2/L835zc/IqBHph0J5kIOImq2NtXQ+IMZDENK48ZnCPPeu9v2uXr7CYtsieir34pXwAuN059yznXJ/COP2QyfMQ8Opw/c3A73rvfUh/ZfCKehaFXu39M9Zp8RDwbcGY82Lggvd+V1w0ISEhYa9xmd1jZ1owbxe7bluwObwOeA/wceBB7/0jzrk3OedeFrL9EnDUOfco8H8CbwhlHwEeBD4GvBt4rfd+q6lOAOfcdzvnTlEMwJ86534xvONdwGMUBvFfAP7hTB2wq/mYTr5FfdOk+46t+vMcxlm//j595Kdtj3l3mzTRukFuRtiT6qI2hgYbTa9XP1diQpUUkcpi5ggxRI/pVFKF8jPWzSnz1FRPVT2NUloCQM2uIfcJlxd7YMyehsuyYN6Ttnnv30VBqHXa96nrIfCKhrJvBt48S50h/SeBn4yke9rPLZ8NWiUk9wJF0bSBucnGK+dSW7rfj6lWzElu0QqnoHHfwjbqkvbO8toYYan1146jqtgefoQqIln7GVGGV6+yT4fKOKQDMRY6+SI9q7ZVpACBBnYeky3j8iKmG5q5rHO/AnwVhS3jFPD9QA/Ae/9zFDTzGygWzJeAf7Crxgakn4xGTIqwVF7p920RC62rr6+cs/oRobJqtsxC2Sc0xJ5gJYBaOy10nVm2E4Fjog4LzSQmDOERSamTZeUZGTHHqhqjUOXF4F7yoQaupm1C0wzaCQW0LSMZwPcejkDVdwjv/bdOeb43C2aDxCgSEhISrhDmNdbTPLb58iDm0xpbqeY59kCiWPGZ1OMxFZRILSZgXxu2oZ1qRPTc6pjYEXELtvdlWqxh1h6UZeEwpUKqWF9XkpgdEzq1IZMs5UjFXGRVZnFhTqqn2RCzXyTpYveYR7k2/WTaEFM/RR7Lf723wWhL4tD6KFE/qQL6JNYZPF6nY0omvadhuwRhaj81rIpvMAAK4i5qtdLeIvlC7CbLgGM2n5qLbM1eUqQndcr2ELNFpTHcGeZVophH5nbl0aLQF08k7Y1kj/aMeiTpDOXZnpPvEv4xiz1hph9ulrV6tzQS/IjtRAc6jHkw1crFPMtUx2Qc81zZOvRYZFm5O1tQjkmLl5h+r4yPZEmEbmdInlK7g2v57FfMI3O7PND6DJ0Wy2cQMy7rcxkmiuh3yUMxbkPN+Bpzie126xvkKu1MJ/6znSLa2N3RE4ReXjLxsNl7CVo2dwmTUPVZaSLLqM7OFu4wHAY1Vb9WVQ15PbaT9nwqMs+u0jPVXhHshRrxcqKNMSQpYzp2a8y+WtjnX8uEhISEawczhPDYl7i+GYUEBYztndBLuwkX2SqMROwMBy1ZNNoqYvsLjO4+y/qlAVY3qdttcY2VumPG+Mhy1WZtVT1pqA53zHNZdU4IISJJaF2aerkein42rhKURMFgUNog9Ps6SkrTeymsQVveI265FldKcmhC0/v3o6RhXWl1epIsmrEPp3Iq5rHNewdvQqBoymapXOSXGnOUEsSYRfTHoym1UfNIUrcbZxbRmEq63rZ7tskkrItXUxrUGKn872dUndCdMUxiMAh8QRuxhVvOAqViEpuH9ooSjyjNJHbCHPaCoWyH+Mfe17QeuBpItorZMK/G7Hls897BMgqNJvuESd/amvR+EtiggOUGMluXXY2rFXOWVcZdTStj4TLK1bN+aUO7Y91stadYxmAJtzxX1uYar9UGbGEUy8tVn6k20nXyzfrZ5FbSGkw2x0p90eCAql3TCOyVkCymvWMaA4h97/YD09BI0sUk9rPRugn77GuVkJCQcO0iGbPnEdpGAZPL1IidwobnFlgtDNQ9n6RsR6+QjatqkxduRCtVkzQ0ynfYCkx/2hy6yhVgkzSh7TqxRhLph0gS+iBxk2lx0aidpE4pY9qvV9TWl6k2XupGx3yywkoTrqTdosHTt/a8wQFtoszVVE3FJIkkXSRj9nzCnnAHzYr7QLQ0jdSGbPvjjdHQXFxYzS/X0k2pW6uuut2CiIrWRvLE1EVtzEIzphgBmfghNzEJu7kjywpVUjD26+RSlaTVTublYoroZ+PKLVZTckXZxXU4SgDLvNX4lS6ykWyx4jFcbmYR+95Mez5N5aS/T1eKaSRbRTuSjWIeMW4Jox251t48o1HdPiHp1i4teeo/1klD6nBYP2kuz2Fjo2qGPnlO6HTsx18xrMLWYQn/rld01s4gUPtAtNdRllEQfssodJvyTQaD0Na1tarTWprQDGM4JMsOTLFRFJelQdtIO5rwtkmHUscssIc/yX6XVg+1GTBhj1HpAisUW4ZwNe0ZySuqjnlkpYlRxCitEClF8cd0avRKslotTOwHrfeWxX70wng2NuIrRtsczVzkuf5fL9tyotk06mPzaiZhz9CQZyrURilNrK1NEnxTZydTA6l3qltVVfhkg3ZCp6PFlowir29kjGnQLFOInfw3q3TR9p2Aevv1e3bCWGKC3yxSxNVSTQmuN6aRJIqEhISEhFYkRjGPsKqniNpJq5u0bUBW/7FVqfbs1FLAwkLdAC1ltBqrKa5TltWPq5C8guBp2oiJVWOTLiImYcWkCStR6PLDIR1RRUlekSqsFLK2Fp+DLJsso9qj1VsTUImlenBQHwArOcgcRKqY6X47aLJj2bq1+nJWKcMOr6S1qSlj79wLpJhQcczjSCRGAZO/IKVusszAqss1wbaBAfVpcfqHKjYHIUyWDsd04hKRdmBULpqOWkZiulO9TP/HPGszXscYhe6YHgxdTqueYvXKdZbB4cNVx6y6SjoXmJDslcgyIGeCOuq9FGLgF4Il9h+Bvo8NUdN1bPjsdRMiX7nyWpfX3ZqlzjYmFvte7Ibp7RTXq1dUkijmEbLcjHkIhTShUevrk6t+LWFYwzYUxEeYgmUa8hqrso/ZKrKszlSags0Oh7C0VLiZtqKJWcSojN5NvbZWEf319eK5jiioOyf1r63BykpVDirLvE6DOpPQ0obtZMmk+jMxQTFoi71G7nWXYyHd7X/NwGO2i1mgJYNeb9K7TfdHhDLN/HX03DbEeLitX+fT97otszC8hO3BuZYtd22bgK8i0lcgISEh4UrBuXZuLwuwfYbrm1GMm8Vc6+W0tTUZqkhLEnrRG4OWAmRVqVeSti6rzterytGokhrk3q78olKFFV/aYFVOIk1YNVKWFcYXvSHOqp1WV+tikohIMhBbW1WQJy1J6AnQ9Yf3d2KHU9SW0FXMKeWQVRtn2Zti97JYV2U7ZLNo7mIr8ti1SJ4CGRY7/zEBUP+X53po9abPWaY8SQ+XGc7N5SDPX4v3EpoKGCO2VTlplbkQdauCgriaotcraF9M9ST/dT1WpQV1RjMYFMRlYaHqhiU2vV69TKn7nUXtpBuhmcT583DxYp2LLSxMUlCrrhJGIe/R+hfp4PJyoXrS6iZdjz4jVU/GDLqYGCO3phRZyGkmop83Dd2seyza1ExWY6fnWRAzZlvVVAxaRSYbNDVmtWfI/RzSuP2FTqf9O7u6euXasg1c39Nu3VzUr8BKEjE7rtAwa6Noki4so9ARYLW0YI3cekGdZRXNXF6urzb1e/W7ypDdTQ20ymxNpMXGsLpa/NeMwi6Zt7bqFFcYhHBb6bSGNuRoDyj9f329yCfMRk+C3nnYgKIKz0fOnOTO4ycYDl051r1e3EFBM2vpzjTJoon3NkkUQuiF4UPcK84ygtjXVuoaqBAl+iws24emcbLtj6Un7ALTVE/7FHvCKJxzdwM/AXSBX/Tev8U8XwDeDrwAOAf8Xe/9p8KzNwL3AlvAd3vv39NWp3PuWcADwNOADwF/33u/6Zz7duBfA58Lr/0p7/0vtjZc/+IUZRVvp9GoTo+mMQxNSDSBgUkjppUsrEGxiaZrzyZZUFvjtkgvwkgaqaC83DIJyyhWV+HcuYphaOlAq4wGg6qjGxuTaiorScg7brihaOzKyqR+TvLZ82WXlyclighFk6584sJJvvP3Xs79L3knt/buLJuljdn6Ndb9Wc9Fk1FbD6M0U6B3acdOFIxdLy9XasY2om1VVVKPbr9MkR3GaUygSYq4nNKF9n665jyh9kD1NAO9/Xa2SwunYNcuvc65LvDTwEuBO4Bvdc7dYbLdC5z33j8b+HHgraHsHcArgecAdwM/45zrTqnzrcCPe+9vB86HugW/6r2/M3x2NTAJ1xbueNoJfuFvvpPnHD1xtZuScD1DVE9NnymYkd7CHtPCvVgTvBB41Hv/GIBz7gHgHuBjKs89wA+E63cAP+UKH7F7gAe89xvAJ51zj4b6iNXpnPs48DXAq0Ket4V6f3ZHLbeHNIePLHQ3NibV7aIul2fWwN0mDQjM60rNS8y4beuR/CJFLC9XLrFa0lhdLZ4tL0M/JvLoayvOaJFJVE/6oxujG6L1Huvrk/aMmLFb0iXGk9xriULyS/mFhSK/7DKU90rnFaToxobj2Ut3srlZn0tZVes5hXaHBatuspKE7oY8s9fSbCuIxVRPYo/a2prcWKnVTVtbk+ouKyzaqdZ2i1i5GJIqahfYvUQxC73dc+wFo7gZ+Ky6PwW8qCmP9z53zl0Ajob0PzJlbw7XsTqPAive+zySH+Dlzrm/AfwZ8D3ee10HAM651wCvAbjtwIH6Q6V2ssZr6/zTxDDaCIvVVWdZfbe2/jQRJCkrC5D19eLdy8tFXQcPFnm2tqp2lgcmxTYL2HtRN0nhlZVC7XT2bPHRex+03mttrc4o9AA2vSPPK24mai2xwkOlr7EDKIOuAwjqgTWQ5mgbjy6uFwKxXfKWr8UWA7MyCr1oFNuEVi3J3ErdonpqW2zqNmoVl/DNpu+Uta1ZY3fsKxIzhl8OFdQ1p3IS7N5GMQu9hRlo4XawF1Mc2z1id4005WlKj6nE2vID/BbwK977Defcd1JIG18zkdn7+4D7AO664QYfkyaEkNiFtbbrWhW8pomacUBcIpBrTWvlXi+mpbytZzCoFvGiTx8MKiO35F1bg0PHB0XDY6t7/QLdeN3hlZU4o9ASRWwXYGwQJH1lpejoTTdVDW4SwYUaycCIO66mghFKbQm8VCtzJ4xC5vzixXqzZ3FcsPMTs1FYPibdOHiw3m09pNIuaa8wldhqXtso9DvFUK/D09u2CDOyQ93Ul91KE9d9KI/pEsUx59zD6v6+QLPKGiJlLL2diRZuB3vBKE4Bt6r7W4DPN+Q55ZzLgBuAJ6eUjaWfBQ4757IgVZT5vffnVP5fINhBEhISEvYNprnHwlnv/V0tz6fS28tBC/eCUXwAuD14I32Owjj9KpPnIeDVwPuAbwZ+13vvnXMPAf/eOfdjwBcAtwPvp+CaE3WGMr8X6ngg1PmbAM65m7z3j4f3vQz4+NSWa73AYFCqnUSzsb5e18CIA9CFC3XtjI2ibf/LqwTWLiHSgV6YN3lByepP9mbIYhyKVaeWTKTsZt4pw2zXVEExPYJ0RDqspYmzZ4u9FFJHr1fvhF522ih7WuwS6ebw4bqNQ0QkSbPL1+Xl4l0XLxbLcS1R6B+f6peW9GS69bYQqCSM1dW4O3RMDWVtSLr+mEQhKiE9TCLBDAZ1NZR0Rd4hwyL1atNazEYhz7TJRlRubcLa1lYVU8ziaqqhrinsXvU0ld7uiBZOwa6nNdgcXge8h8Jd637v/SPOuTcBD3vvHwJ+Cfi3wVj9JEXnCPkepDDE5MBrvfdbALE6wyv/KfCAc+6HgD8JdQN8t3PuZaGeJ4Fvn9p4+WUNBoyzfo34a4ZgbbqWWeg8mqjo+FAWVmuzHWah6zh8uK6OWlwsjNsSNknK9rU+I6YH0/oVHZ9JM4mzZ6sBEGhmAXW9mXVxXVlhPBpxKRQ9cOZMoYgYDgvCf+xY3Sgu7ZK2CiXSuhTRrdg+ZRmjUd3FWbq6ulp07cKFepdjqsSmBYDdX9FmoxA7lHV1FcP0wkK1k16rofSULC3VjdV6qDWjGgzq79GeyTINAst8hAFKW5sM4xaJQWwDuzRmz0hvt08Lp2BPptd7/y7gXSbt+9T1EHhFQ9k3A2+epc6Q/hiVZ5ROfyPwxm01XC3xhLAL0bh4saKXsgiOfVZXq7zCKLwfAZsUW0NG4X8do1GX0ahHcdpzl16vN8EstCeMXaALbZYftyzKV1aqLQagvpMxpqAzSEUXL1adE0YhUsWZM3DuHOPRCMv7Oup/RxqtGEU+GrEWRmVI8cUbUzCLbHUVbryxeP+RIxUFte0UiiQbDPQS3yros4ytrbqEIM25cKHo0vnzdbu4NuNo5i8MfzTaCq2WuSX8HwfGUZ/n0ah44XDYYTiU62K+AVZWCuaqeaN8B6T7UH0vtbFa8ub5ZMDIMOQ1RiG81fJgPcTC1GIh15sYRMI2MV31NBUz0Nvt08IpuL7XAeHXNM76rK3UjZzaliuMQtNMTUtXVqDghesUhGSdgkEIQZFfnZy6JhvPusAA6DMa9RmNBqyu9kumYdUImh5q1UKe18OQ62McxPulDMMtlVp3VOGUwijEcC2SRGASw9GIzdCzcfgIhElkQDYakY1GdMKrh8BKKJdTsEdJPzAcsvxZ5ZQhW881xRIqLxSvaSt86ONm3plwRhCcP185c2ltmDD9PK+YSDGvm2Y+NfPfMtca/ZAmzGErzHfRl9FoCdji/PlFoIdz3Zrnmkgtw2G1x1FLAVrykP+a4cQcLDTDjDEKGxVXSxYCe51UUNtAivWUkDAf8N5zpnsSzwniTiQJCZcJ13MIj3mF9x4Ggwkpwu4xGw7rKvqwuGZlBUajIXCRQoq4SLXiFKlCrzgF3fDpUKw2+4hkAQNGowGjkaQDdHCuO6Er1moK+e4tLRUL8vX1akGeZcDacNLooe0VUBRaWSl0M6dPFwNx5kzxWVlhGOwLssbOw0dLFZnqTRZ6mANPhY/k74d6BqEugOXTp6tO6ACAo1HRQVEVahfdBm+BmGszFFX9+epJ/uPSy3nu6J24s3fW5rzYS7EVWidSokgTIllsqVaPqUsUXXUvEkWPSooU6aIT6t8CloA+3ncYDvsMh0JIuqVzhdWuQeXIsLRUrfq12kqkCrnXBnlJ0wZu/Q65tjYYa7gXJDXUjEgSxfzBA5eGndK7qckO0WTT9V7I30WKH/0aFVERIiBKGg1RP/UpiIiQ1sVwr/8XRMX7btB5i9qiS54XxFQ8eaAwYkv8PCEwHcYVpZACULfAy73dM3HuHJw7x+ZwyBC4FD7CLLQKStieEP8+FTl8KoyOjMQwPBtQMI8OhbpqcOZMZVHtdis1k7RdmISOn22RZayvFqqk1dV6N/MceudP8Lwn38Hj57dYPfAAfOpbuHC+w3A4Uq1co2IW0lOZU70AiKkWrYoxpnrqhvpQoyWqyGJ0zp49QuHb0Z1QN0F9l7Zo6rSBvG1vD8QZgaTvR8J/TWzC2wMbxdXAdc8ohsPCTr0gAAAgAElEQVRiAS2eMJpBCL209PP8+S2qNbJlFpvqvyUo2jbRpWISwhSEWVwM/4cqn5QR8tvF+y3W1gblDm3t+mn1zxNLyiwrKIsuIBRVOrq6CmfOMF5dLZmEZhaXqCQLIY0D4EBIl5/DJSgN2WPqjGUYrkUSyYZDsosXK2u9uC2JhCGhyGMBHdX9xkYxr9bRYH0dnrrguLDieOzFL2W88CSsfxoef21oqUiFIiVqqUJsUMI4YFKi0BCJQpiDSA8yl1KHZRSLyOKgsGEcIs+L7454P9lgu9YmAfX79fVqiq1Pg2UQ8lyfBCi+DnNI4/YXkupp/uCpiMjFi5WRU3/kCIYzZzSTWGGSUaxREBetdtLGbL3KtIxCK2Lkv9RTGbwr4lKppbzvMhz22NgoXCL1rnL5kY/p0BEqIRtEZDmp90xo958zZ0rDtjAIzSTWqNbcwigy6iop+WiGQvgvJHJAXRLpA4dWVytVkz78o9uNx4+K7LuwXr6yR0L2T2x85gT9R9/N8Phvw10/A489H04/K8ylZhSxOZVFANQlCq16qqS/imUSysn3YFOV01KH5OmXda2uLgFw4UIVgffIkboaUguNMgbCGOz5KbJOsJq7mCbPnqNhn+u07WhVrstd2kn1lJCw/+HxXDjwYYaHofvInXBpCx77G3D6i6Hc4ZGQcJmQVE/zCW3IFP964xFaHscwqXI6T131JDYKUVk02Sdk9aglCr1Slf899UzsFrJ6lfKVul5WeXqTWbny026xVmkv6ibxD5UOB7uEVqbJtUgUl6hLFNq4vUndjkH4L3V1VLlBqEvUVtlwWKw3g1QxBjpaVNJ+wlqiyLLyGFuZV5EK8xw+s3GSh2//RkbP8tB5Ldz18/DgfcCqmkuRKETtZI3as9ooRCJYpNppMlJzp9WRfSZVVYvqWTHna2tL5XTqWFV6k6dAJAYRIK3twQbmnfW0vpir7Bwukq8e5nCw5q/FewjHpJpJ3wvjEIJTEA+rctoOoxCVg9geRMVkVVRyvRieD0K66Lc1sZlUC0BdrVCDtkVoPZVWNw2HjFdXy30OlkloNZSYfoXgW4zNtdXwE3mHGLc71qlfVFG6Y3qXYtCzCHHUCwAhpOOzJzh+6rc4Nfwf+K/9J/Cxvwmnb6WYP/nIgmCnjALqxmrN3CXPokqTdN3ffngnVDaqLYbDglF0u3WtnGwKFOLf61XPrMcSTH437DO7CW+7SHspGpAkivmDo+4Sqzcla8axuqrdYEV3vUaduDxF3egpHjKakGhJwh6C3I2kYdKkXF+ldWob8Hq96qNXkX2oJ0gcEh2XInRcG6i1rWHTfDQr1AZq2XWNuu9E6utE6s9VXUX3i/53RqPJmBVKohhnfTqBMon5RU5vlZhOW1uwMXSsn7oT//SPQv8CvOjtRU/e853wnHfDI3+VwnqljdkxuUrLTVBnBtrmYOdSz7Mw/b4pL5AlvrDXwvvJ+y3yvFtuwtPeS/Kx/FRLm1D/fogNQntAyXjtFolhGCQbRULCfCA/9mF46T+BzrhYLbzwQWALXvgbwHfBI7Hw/gkJe4Dk9TR/cEyqm6yEUfjhb1GtLkfUV5nr4aNXnNo9VkO7QcKk+kI093Jvy1vVRJ9er1tqXyRmkJYwahKFQHZx6Uh4wbBhpYTcXMdW/7Jwtaol0dJn5jnUDxwZm0+JbrfyfhJ3WX0EXJAmNukzXIMs6xduomv1Ll68WFW5vg4XLwFbGZz9UrjtJHz0K+HkV8OpW+CR51LNp/Zy0tKEDuGhlWhjk67dYLuRdD2f9rtgvakKe0aImVli2n4Hq0baDfZCwrjukVRP8wfnXOltKWp77UFaxdgZUScaYlPQtoXYDmwNvQELJlUQ9mPzyL3U0wMGNSaxtDR5EBIEAqGd4wU6yFGgApohTBDugJh5XmvnhUmIyqlj8mbmPvbJoOjIwYOVMl4YhzFiy9x1u5PbRYQXanWKf/xO+Njfghf/YpF4xx/CF30Q3n1v6Iln0mYEk/M7K+WMzXvHPJv2HYipJZshdiv9PditxqO7vSYkxJBUT3MIRWT0hiwdc65YwTUxAfvDFg+l2A/b7p3Q+yPs/17D/z7Vxrw+g0GXw4crJnHwIOX9RKhyfaMjrwpnDIrqLDji29W//kha3+TTtggh9lqisM8wefSn0+sVHVoqjLclo1haqs59DcxCDNexjWQyp70edLIxq7f9GuNPfTOcej6MO4X6KQMGl+DuXwL+HjzyAuoEXQf3EwlCJIImZlFtjKzkOT33Mq9246Wd704tX6/XLQ3Mml8KEbc0yDIJ2cMoaWLXAqoAklT3u8Uc0sTLiyRRzCF6vQnjnz2MpkCMIWjCrVUR8qMeUN+RPQuj0OE7JL1L4e20GP4PgCWcWyoZw7FjVXTuGKOoRY8V9Y2GcMZut+yVVjtZoi8ts86gHVWmT+HqKj0TWDWTuMYOVN4+VIdqCKMQSEflhJ/BgI21OKMQyJyu3/ZrnH7B36ObfxL+p38DW4FRyBRnW/D1vwbnbobTz1Av1SqkHnXPNMsorCRgPZz0nMqcL1GfZ8mn3WSL74SODquvFxeLPne7leAl3kp20WC1d/I/FhVlp4Q+MYgWzOHgzF+LExJ2iMOffwWbm/DEJ58NL/oRyPKKj+fAH74cLhwO7rKb7ZUlJOwESaKYQwSJQkPEbxHte71uWGlJqDtrlygrox4IUEsZVo3RV9d6HW0lChswcIlCmjjEsWNw9GghTYhUMRhUmhktVZTQm9PkXusegrrnwOoqHdXbsWqFxHKCaoOdQIzcEvNpmXoUWYFVaw2AQ6HeA1B06MYbi//6qD55Jh0M7ReJUM5O0N3NssoO3s86HD/7d7m0OubCZ++CL/ndSgDIgK98B+QLcO6L4PQXqPmwQR7FdgGT4TtEJWXjc6HqEylD5ldLFH11rdMLtZMMR5ZVpxvqg48OHizSBoNqWu1JigL5rsckjTbExngn6DC+/sJ4JBvFHEKpYOTY58XFwoB98GClshgOq1g7BTQBKFRB1aYs8VbSjKRJ9WSDAsbUUloNsUyv1+PYMWofYRgLC4VmxqqfSgjVXFioYlMLdxG9Q7dLZzDgwOoqg+Gw3HktP2e951ivuWVzXk7xpdJqp9LuoOqwjOIAgVn0enD8eMX9Dh+ud0L0a2K5D4zCqgw1YZRjQYVIHnjWSS4c/zg8dQyWz1Y6sE4e5uEQcJTmuF3aG80yCj3fbTYKPb91hlCpm4r0gwd7ZFnR9WPHqv4dPlxMm9bQ6dDjctKhZhB6jSB7J5rsGvJ8JwwkoQHJPXYO0e/XzqvW0KvRjQ0YDLqsrBxiNFqk2mwnzEECx7Udf6o32mmJwhoxuyatKDcYdEubhHw0szh8uOiHHDttbb5lhxYXi1//wYPVylyHIBWj8XBIZzhk+dw58hA9VtsnhClA5RIrkWCtnGQZhZRBjcIBYLnXg2c+E266qejQ8eOTZ2gvLBSdPHiw2GhnVqQ6sofwFLFfCOM80x/z+IFzkG1WjYDCXvHu7wWOAcepZCTN+O1u+yYbhQ4A2VdpTQuBelqv1yuFpqNHq4WMSBQy13leLQygfj6JDN3CQl2qgMlrgZUQkqfTHiOpnhIS5gObS38O3S0YZdBTyrMM+Os/DYMtePDHQqDAhIQ9RFI9zSGyrFyBaX3twkKxGltbqxbcEqr64sUea2tHWFs7wmhk4zpZlVNM9aRXlXVVlHO92srX6paXl+GGG5iQLOQjK05xldV66povpSxDjx6tluCSbo02Z8+SnT7N8pkzMBxGN9uNoaai0t5RNZdXVW0txEevV7Xn1lvhllsqPYuIRYLFxWoJHdROMn+6m1pVo7vY7cJz//IVnDnz/7L69P88+Z1YWoX//GNw+u7QI32infZ2EojKSaue9JzbA4xkQqpDqQaDLr1ezexSUx3ecEMlUeiDiSTMuMw91CUO+e7q6+XlKmSHtVvIOOmxFLuPVlnpsBw2XaftNeb+0CJIqqe5RK9XqsBlb1ds452cY2Cfra722NjoMRweKjd5yd4wqP+gxKCqCZkmXkII7PNut+YJWmqL9EfbJBYWKq2S9rUnp05pxBIqL9YHGWgcP15Q28cfh4sXyUajgrBD0dlwVmweyncIeyCkY9LwGtcyEA4n6iZhEIcPF51YXKzyinJeqCqVRk2utdrpllvsWHs+ffAdbCx8KN6Wg0/Cy17HLUtfjnv8+bXtJtJlvUNZ5tjGNBJiK/YR3TZZjGhiHfsuSLplIFKXVjcJ7dGqJ+H9ci0HD+n3WaIv9Ujbmlxpdb8TtoEkUcwhAqPQh76MRoVNQsfKkzS9P00/t8HYYqEOLKOAOFOwzETyWcKhddB69Sn5hL7WvpM60qo0amGhoKiaUWjF9MZGQW0l9jpUL5TDxIdDstXVOLeb2NBhoPPdcEPFMBYX4xZ5yav6oVfEerF2/Hh1PKi85tTWSd4//md85fh7eC/fW2X26n9/jS98/ie47cLza3MLk6G49VzrYdM7o21XYwsB3Q/NNKBg/Dpd6tRMQdL03gr7PdPPobZnseybnSLdJs0UbL9sWkIDrmcbhXPubuAnKGTpX/Tev8U8XwDeDrwAOAf8Xe/9p8KzNwL3Usju3+29f09bnc65ZwEPAE8DPgT8fe/9Zts7GtHrcfx4fMVod/daRmDToTmmTpux0DIEIQ6W0Fi6qyWNNrqcZUFkj+k0BoNKf6F2Z5cv1R2WYFjaeiqMQg5BsJwwxh0t9NJbuJtQ9ont5QGaUeR1oirjJtoqORRPqrnd/VXWF3+YR8/+WfEtEThpD0CHr37u7dzaa5/XtnmOpcUId9OcxxYTTfk13dGEXa/8bRnbHvveWH9iU3E5mINWMV0T6iaNPZAodkNvd4pdOzE757rATwMvBe4AvtU5d4fJdi9w3nv/bODHgbeGsncArwSeA9wN/IxzrjulzrcCP+69v53i5KB7296RkKDx6Nqf8kuf/md88OJ/ask15omtR69YmxKuI4iNoukztfjO6e1usBfrgRcCj3rvHwNwzj0A3AN8TOW5B/iBcP0O4Keccy6kP+C93wA+6Zx7NNRHrE7n3MeBrwFeFfK8LdT7s03v8N6LUmESWcbx49WtXiE2qRhmPfyl5ZWt901p0C5t2HJyXa7ItF5KrJoxUUhXYjsrYpdWcNvzq5s6Mm0V1bTclv/2HVlxkp0ISbFzDxYXC20WVGP2nO4J1gY/zL983/c0NuXp/Zt5/UtfUVMl7VQX39btWb4L1j01ZuaZZZhnacflPDtCNte1bbK75qSHGHavetoxvW2lhVOwF1+Lm4HPqvtTgA3oX+bx3ufOuQsUO5puBv7IlL05XMfqPAqseO/zSP6md5zVDXHOvQZ4DcCtt97GIZ6KU39HOyfYCZeAKuZS072kEUl3cj/tf3Wtf5Q5/dDMfhkhVhhiqWsP8YKsnWWYFx+AQQaZGHnN1uxak/X9NCIkKhJf+QXVemWJqkrqZ+OiR2YO+hkcWrT6QfjOO76OQ/5f8Pd/77V4Jn87T2x+jrWVd/G8o8+pN3AniM1vrRdt96g5b8kTmfdonkZOUKT1y+2SlwfCCK5p1dI0TFc9HXPOPazu7/Pe36fud0Nvz7JD7MW3wkXS7K+vKU9TemzJ0ZZ/1nYQBv0+gC//8rs8p041GxyaDBNN91IuBq23h7hCuk1k0Pr+mP9szFI6GJSnvm3mnZoXl42aK4KBbo4eDjlLHKoNfRsbRZqOAm69ZKyOvWlYmiTxJpuOmCn62bhyRdOQTgWvrNq8eM/7PvKfokwCwpfvj/8Y+o/XRUvLPWNzbTvaZOG2dhxd1hoVmqQtPXiCJiOFXOu8VnKzjGQ74m7CTPC++C224Kz3/q6W57uhtzvGXgRaOQXcqu5vAT7flMc5lwE3AE+2lG1KPwscDnXYdzW9IyGhhpNrj/Lvzv9+4/N7e18BwC4k9YSEKLxvdpSZUSmxG3q7Y+zF8uADwO3BG+lzFMbpV5k8DwGvBt4HfDPwu95775x7CPj3zrkfA74AuB14PwVHnKgzlPm9UMcDoc7fbHtHW8PznGJJnOfFUlv7wNo07TdrfWOtj2xZuUJseR2TApryaX9X2SQhu+pi/rHWiT40SU59EycmfXS2bqpAJI/HHy/yiXlDnJ5On66GRJocE36apAq9b0CcqY4dK9KPHJn0ju126+dulOMvopFkFmni9OnKWwvw3S7jS4/xqvFz+Rn+uKpYfVMuPPkZXn7pJ/gPn/sK7rx4sHlep7m3xdLaPlYakHsJ5CQDpQdCJsRKF03+sbrN+jsiz5r059pOZNOSlDEzxuPmLUszYsf0djcv3fUMBx3Y64D3ULhr3e+9f8Q59ybgYe/9Q8AvAf82GKufpOgcId+DFIaYHHitD2c9xuoMr/ynwAPOuR8C/iTUTdM72jAaUakm1tYqxrC+Tu3ouzxv2nFX31ixtTV9I4X9b3dZWSoLdeKvj7OTnXV6x53eVCeEgLomZm2t2hZx5gycO9f85RUmce5cUVZ3Z2ur8o4dDkVF0ymPZ9V8SzaYWQgNk+0TR48WoZ4WFoopWFqqlxsMiimSMv2MSUaRZVUnT50q93qQZZzMzvIy9ys8NR7WN1IrPLj0GX7kQ4c5cfKjMMorXZ18afT8ikux/NcDBJPEXeZeNiDqa81Vp+24E0YB9d3rkhZbLAhTkYGH6bY2KR9LtyqrxDRmwk4dI2B39HY32JNZ9d6/C3iXSfs+dT0EXtFQ9s3Am2epM6Q/RuUZpdMb39GEklEUsTnqZ0jrw7OFUUi+8D8fjcpwFvIf9b92QtxoRGc4nDgiNFOfDjQTkMGgoJr6hKLDhwu3nvX1apOaDctdrgT7tW7JYlsW3NJN6+114UJBa8+cKQi09zqiKlTByJ8KaV1GowGjUezEvjhl7oWNjxLBY2Oj6KL8typ42bid5xSWde2RJdAShdihsowTWZcf9V/Kvbc07MwGshF8xcMrXDy9MhGqxB4PO4Zy0Dpq8DqjUTHHwyGd1dWiXqofnISF7AAdmWfhisIEZOFw5EjdMCMDIbvrb7ih2rVXnN1bSZx5Xr+2BF17k9lrGWRN2XS+3VC8GSBeUtcS9kCi2BW93Smua/af5xREf2WlYhRCNYUxiMRx8SKsrJCPRgypojvpjyUqGppB6P/20xkOyYbDGvPICMQkFj5WJJmtrboKAWpSifzeRfgR4/TZs3D+fKV+EhoiwtXZs7C6OqKImCs9fyrc6/OkL6q0AdVpbZpZaEZRxcEajQacOXOIlZUjpY+AhNAWGikQLczhw0UbD4n6yRI5zShOny6fuSzjFac9v77R5cEviRO6jocLwOkwlzFmQWSOoX5UbNNZ4BIPqzw5METq7a+uVnOuGUdT+Fit+rThj0XCkAGVa6uO0szDMhDrNt2UZtVSSbJoxF4wiquBNJsJ1x1+7aYV3nF7w2p4E370t+CO082nYSck7AaXWRC7LLiuGcVoRLFsFguv1stoi+9wyDicyXCJal0dkyxiK0+7ooxJFDbSqs7Xp1h1DoZDspWVYoUZM6Bb10eRKPKcPO/X7L76c+FCIVVcvFhVdf48DIdbTJ678RTFhviVkFaOJoVEMaKQHg5RnXE3oIqmCvXouqKiKqLxfupTT5fjMIDKHgGFFmZtrVhcr60VC+RCPdEwuaur1TxKBcDtm2v0nwPD/mSxL38c/pdHih6KlKilCplfmFRB6XZk6r9ca0lC5lsObpL/5ZyPRgxWVorvyXBYqaPERmKlR/0dELVRr1dJGyIm2oBScrKR1GHrS9LBniJJFHOIzdoRbUovs75eqZ1WVhgHddMlKkahrzWziKmf2hiFZRB9c62PN9oEBqMRB06frgzngiavGfWttOonzTi07X51VRiEEH85mugp4AkKMno+5NGDKPaLbsh7hPrxngJ7+t8SYvMYjQb85V8eKhmE1XCIzUL8CPIc+tZDrHyNCvsLpQW+k29ycARDzbsC/scz4IPH4VmnmbA/tamfYoxCz7mkCW8SppBRHR87pJrzARVDWhaV1MWL1QvEhiWR/ISw6x33ki7BHzc2qu+MREtsYg4xv81pS+ErzFTm0YYh7rHzhuuaUYxGTH6xc+XlMhoxHo1KYiEfLU1YZiEk1dopYoZrzRhiZzhYoqEJ1IEzZ4oLHXRPez7JoRqma+vrkzRgNKqu19agIPLCKISoX1QfYRRPUZ0lLVRXTvkbhtaKhKGlCTlzWs6WFmZWSCIrK4eAemgObTs5fLjO2PqaSYrfrnTYWuhXVzlx+hLf9d/hh77G2Bl84SUr86qZf0yysIZtiNsoNNOQOd2kmuuc+jx31Htq35vRiL5IR4NBJSmJ23Q1gXVbTZYV32kxPNnnlkFMI/Yx5hH7HSVJZALeJ4kiIWHfw+N58MvG/NhXwtg6YbkiwkJsW2tCwl4gSRRziNGI+tF2DbvC7ArSrjCt51OTRGFVT/2QL/Y/C+WlNbG6BqurlQulDtGgbBN69acjj2gsLlZmmsKzbkj93G854W2o/msJQ/KJikkkBjnhbVHVY88U71F5ShX1en+RtbUlVlbqsai63WrRrPdEHtJ7DrSdRjqq9jecfNomr/+qnMyDy8HLAOdFc1/5B/CFp2GD+JxOkygETaonkRLG1NVVep4zkybS5ZjCVsVoVHzEA0psD3q+s6wSuWTgbJiZJoq1HXXTNOyBZDGPKqYmJBvFHCLPqR8fZ3c293o1FYJFEwPR95oYaAIRYwAWkl/erxlTH8jFuD0YVCqHlZXqvEulVtD7xDRPXFioul4Q5U31dqiIu/5IujAV/UyO+JS6htTVTiOqY2OhUkvV3QJGoz5rawWB17RPDNk6bhXLg/oZFlBXyamTeU482ee3Huzy88/Z4udfXI31l34OPnEb3HyhUD/ZeY39j9kqOurefm90mv7kqqz+rogaqtYO2bMhDEPUpDGbgkYT4Re9Y5OLbKyOhB0jqZ7mEN77gsAUW4urj4lMp/c1WBtD7IevoYm9hmUyTZu5OtSJkWVKmewQlw0GoqPXokOWlUxCx6iTc5e3tgovp8JJpktdMtDfan32t4YwD2npQOUXBtKlsksIs5ARsuWLOtfXK0YxGhU83R45OxwCx9VuY1lJa7uNFJKd0Z0+rzq5zgdugg/dVhR9/Cj87/8FvvqROEGXllrJbpZ1rp7/mBG87Xug7y0TaoVmDNMCViZcMSTV01xiXO3sssfWKcbRoRD5xfNIeyGJukhURvoHrFeDct+0ASvGeDRDkvKCUnoJO74RjxjZrdsQvElHBzl4sHIMknTnlvBeG6W1oboT/vcoJIcelRJtK1z3qDbZQcVUuqouC6m7U3su9lfN86xWJM+LaJz95eUyKmdfS4diEV9bgyzj5DPgG1+8waUOXFiAOz8Df3YjvPp98Mr/XpfY2lRMTfcy5zJ3UJ8/O8cQ/17EnkP9O1AMXUMcEoEE35LrhKuKeVU9TV2YJCRcK/Dh7zd/+xDf8/9BfwQfvg02e/DvXgSPHr/aLUy41rEH0WOvCq7zJYavVBM6MuxwWJc0ul36584xCK6yOXGXVa1X1i6OqOf6E9tk1zef2H4LjTHQ0cEI9TfOGHf1wlIW2qKlkegQWQarqwOGwx6FimhM5eEvITnk/zJ1F1rJs0yxNwLq+yVEyuipdN3bRSq1V9FY6VIRK7Jbth+q9q+twWDQKbtdussuLVUSFnCy8wTf/MWn+OEPHuHnT0BPfrQ9eOkH4ctPF72xUqG1Eel5FwN1ORemnPRQ0rQjg722btExSUTsZrVzLOyZFfJfSxuxfSZi9IlFmbV5Y/cJ20ayUcwlfH3Tkd1UIPrdsFFrcOZM1IAJdWIiP3jLKCSf9auPMYimjXjWcwqYjOQHjSoHMWBrjYwOISTBac+e7TIcdhmNtiiIvhialygYxyHqgQGFURwKeY5Q7ZGQHdzCAMbqXuJCLVHt4i5IZeykPU3LtJawtmF5EDJJbKzAKE4Mu7zjL16A5zz/6o/WeO1XrJF5yD28+FTVMzvOej71BybtCdomJUxEqyQ1E5Ad2XbvjN5oCfW5r0WitdGFxRAtH2EC8j/GOPR9bE9EYg57imSjmFcIIdGMwlKnwaDY0TsacWBlpdXtVbvLitEZmj1epjEKu/qcYBIwGb9buIFeLQboowuC2r7cl7e8XI94UUTqXgqrefFUGlF5KMmGuZG6XqZgFsIwZEe3lBEvJ2nzksp/CFjEuUFJ3yp0Swc1/Uycu6RP5SZK4X6HD6sDM7r8WQ/e+OxP8po/P0TPr0EO2RjuOAcHwthuqv8iS9k5nYVRWIkCKubQoTL5y3tlrvVGS2E0pVwoTgtyLdxRhxTXbtLi1maPH7Ru4bAzptDgUr6rOq9hzKuN4jqfRVf3irG7eCVNQiVQ/ICXQwwea3gWw3YsRpDAlompnOzKssngXaqihFDIKlPCVpsf6eJixQOtykmcv5aXq4XnwkLR/fPnJcSGMATxWOpSMAZhGhKOQ5iE7AxfpNjJvR5GQ1RXIlEcCZ+CWejjF2Q3tgRPFYlHFsDi+aTdfZExkXDsId7HSf84b1z+CN+4cgs/9+xP85N/dAT/1FO8/mu2OHUUnne6GnthDHIt82oZhV4cWunRShSyW8QyilqMJyYZxQDoyDkThw9XUq4+uEi+x1rKkMGw551A3X0YqgGNEf59qoKax70VSaJISNjnOLF1jB9ZezFvPPTfeOsnvohv+eRFOHOJzhje8LVbPOsc3H76arcy4VpGslHMJVz9ABgtUWj1k5zqE2ClCr0K1RKF3pyl0SRRaElCGzWbXCYzKFaVIknIiXf6aNQsYxz0/YNBsTrXi0h9ONrGRv0Uu6BxI8+7rK4uUd9YJ9LA0KQtUrc5SL4ORcRZ7SLboZA6DiNSxWDQq0kUWlq44YZ612DSdl9GlBVpYmWl9AN2ec635C/gS1cPc2I8wmnbdT8AACAASURBVB39FHQzvuWJdb70N87znNObbKg5tLGeYtGBZ1E9ZSpdz2mTRCHXB6Suw4erk+vkP9R1hnZn+mBQHWYkOjvZjCKwEQlixuztqJAuk6Qxj5JDE5LqaS7RYTw4QEfEdw29SUmHhVDHmB5YWSkCtVGRzDHxcA8Qt1O0MQqZnJLoaMW8qBLkSFTLJIQzZBkdxvR6nVo3RqM6wYU6UYaCplT8cxDCjmsmIRvntE1HRmOZyh6xRGXI1juy+1SqqqP0ekvceGP9+FOtUj9ypKD9tt16M7FEq+jr4IjqZDi3vs4JvpiTN36ME2vHcL0ejDbBQ/bUkN7FS4zDyYVtQQGbbBQaonrSNiWtOtSqJW2DKm0aN95YnJontpatrfoJhrLAkXnX32Oth1tYKAZGggeWDWywUewT1dK1iKR6mkM45xgO4YA+YlJ7O8n9ujp3QX5EgVD3V1fpr6wUIcBpZhR2051cW4+XDIrTzaw+WVyAhElovbM+9U4YhRCFoHseDPpllI88L4pqYqxdZeW1q6t1O//p00tB4BJmYHdUy4Y5ea4jyorbqxyjqr2dDuHcEY4fLw7tE0agebSWKLTEIc/0/XCoGMXiYsF5pMKFBU72n+Tl2X/lnZt/jTuzY5w88BQv/6IV3vmHN3Hn5w7RGQ4ZrK4yGA6j9iYtJbYxiphEoRcA2lhdMhORDKE4dwSqRQAUAySURh+LKoZtYQ6Li3WJQq5jgxbzcoo9V9+nVq+phEYk1VNCwpzgROcLeOel/5UTG1vABU5cOsg7P/RsTqzkFAHGExIuD5JEMYcQNcVm1qevRXe7XbJ0paFaVcmKNah++qur9C9eBLUpz6qeNCbsFFqVFJMorNuuliiOHauW4rLk1q5BeQ70a4KS1lBAtXrXrvaDQXXWjTiGnT27xGgk0kBtNMueOdetSSZbWwNGIwkWKCflLSJeToNBjxtvLLpw+HDl2avNRdqJSbRJuv3yXzbojbN+cc54lhWraVWhu3iROy9+CVw6DVkPB9x51sFNwyriYHCr7a+uFhFbt7YYj0Y7kihi9gpkzqGYc21b0B5OMlnyvVT7QsprUT/meVVWDiuC+oDZUxCtHk9/ITSaJIY219iECSQbxRyi262Mu/1BVvmg25O/NKOAOrOQzQfihrm1VQQRFOpqT6ITaGYgTEcbpq17q3XdlTq02kn7upqjMnX4J61RqIh5pVHQG35lC4I+72Ztrdo53aTmtkFIh8MeKytPDxv4ngodOMSNN3Y5erTaGyfd0N2WepeWiq4tLlbDpEwxtb0Uea52aAfXUb+wwIc3Pg0L69x5ww047U5adKz6Dugd+mFHX2c4pBPmt8bFgjpwPBpN7m8RY5C9hmpPhMy/XIuaSNRMTYxCq57sQscas5uYwiz7IJrUTJdZ3XQtGbE1kkQxZxC1f6Gb71RShf7BAaU7DUyGrxZGocN+CIOwBEW/2NoZdJhs++O0kfCkTXqZLRRWMwlBlrG+WlzaBaWGdqOXEzZNtPKSYMeGwnoi6R+EDE+edzl79kh5f8stBT2UZmsHHl1e3iGLb6F9ehOeHp48D0ZhRYRPDj/FN/7FD8LY8x+Pvo47Zb7seGpGoQ++0KFStIQne2zsRs2YjUmLbXIGtnTCtLdmrBYGptss14uLVWj5hYVKJLN1asZovwQxQ3aMQcSQJIiZcTklCufc04BfBZ4JfAr4Fu/9+Ui+LeAj4fYz3vuXTat7VzO8jYa9Gvjn4faHvPdvC+kvAH6ZQg/xLuAfe+99U73OOQf8BPANFMrkb/fefyjUtf3OZ9VCfTiEbLlTqCs0s5CMMKmXkR+fbHbSG/fsEZPy47axdaSeJopr/8u1pqrCLGKMIssYZ322tqrqdZdiYyIquV6v+K+PWpajl3XzmxiF5ZOSfuxY4bUqWrPjxwsGoI8G0dAMTOfRQydTMdEZ+Swvc+KGE/zW8k/ApUucuHgQOhfqUppUoucvNqcxxg2Vo4GeazvYehGgV/96YuRaSwySZndmy3vEmK2N1ppRtKmZYv+nLXuTymlHuMyqpzcA7/Xev8U594Zw/08j+da993dup+LdRo+Vht0OvDfc1xCI/vcDLwJeCHy/c+5IePyzwGuA28Pn7in1vlTlfU0oL1j33t8ZPlOZRML1B+cczzv8ZTzv0JdQrDkSEq4sLnP02HuAt4XrtwF/e9c1Bux2CXAP8FXh+m3A7zPJwV4C/I73/kkA59zvAHc7534fOOS9f19IfztFx367pd57gLd77z3wR865w865m7z3j++k8dYAWvrf29mzKyVR3ssK1EoTbd8Aq/vV0oT+2HL6Xksiei+FXOv6l5fJ83oMOanONsmOiVU9iTZDVkTyenHR12YRPSz2nXlemnMmTCxSl7RFlwndKYUm/WmcYCtGSQeXlurShDbINJ0cJ52xEkXTL1zr8uVaN9au+KWMliisxKElCm24lk0vUsb2WyauTaLQ3/UmA/dltktc67jM7rHPEFrovX/cOff0hnwD59zDFH4Zb/He/8a0infLKGZp2M3AZ9X9qZB2c7i26W31NtX1ODN23jn3GgpphOXl2yaMrlnWoSM/rCaCr4mK1WHHVBOW2ci1EAgbsE3KWb2NlNWbpzSj0PUIVR0MSiIviEVkjZ1tY49a1vRO7qUL4qwlzyX8twyP5bXr68VZS9JMvX9C03PLr4VmSgTchYUWjUeMAQvxlMBXYjvQcwqT8b+a5lb/j70/Nvd6/vO8sinY55ZRjEaTO7P1c82UYguRaaqntnbH8iVsGzO4xx4LdExwn/f+Prlxzv0XIHZyyvduoxm3ee8/75z7IuB3nXMf8d7/RVuBqTO/Bw2Lyfi+JX0ndcGMnQ+Dfh/ATTfd5aGiBfL7KKUKyyygoqryX//Y9SY9yW9XrVKHps6zMgmoR4VdWCgIniY0MGEVzogbiKHeBF1c6ISmj/LKjY36vT12HCpGsb5et2loSGQUkRLkPAydV/NgaZs4iMm7pV/C2GqwRFLPrTwfDqu5bGIUei6sXUKgGY8eXBlgsWHYTuqzJSwjkf9Slw7BYSUG2y9dV5ukYMcq9szm07+HhJkxg43irPf+rqaH3vuva3rmnPtL0bA4524Cnmio4/Ph/2NBs/M8oJVRTLVReO+/znv/3MjnN4G/DA2ipWGngFvV/S3A50P6LZF0WuptqqvWeQpV1fOm9S0hISHhSuMy2igeAl4drl8N/KbN4Jw74pxbCNfHgL8OfGxaxbtdDkjD3tLUMOA9wA8rA/bXA2/03j/pnFt1zr0Y+GPg24D/Z0q9DwGvc849QGEcvxC45xHgkvd+Q3X+X01rvD1uuFo5B/WT1u1aRXivV19eN0kUGjFJQqe3rWB1WStR2BWeVuTnOVnWLxfOMcScc2SBqkNgSRO0OkmrnnQoIVncau2cXaxKaHORSGJ7KKSctH15edJFdqra3HqWWZWMHntR8WhpEurusE0qp9g+F5Ei9CAK5DszTaKIqavscytR6LzTJAU9aU3PZk1PaMVl9np6C/Cgc+5e4DPAKwCcc3cB3+m9/w7grwA/75yT0HNv8d5fdkYxtWGBIfwg8IFQ5k1i2Aa+i8o99rfDp7FeChfabwAepXCP/QchfUedbz2XPqbX1s9EEW93Tc/CJCxltkxCMwjxjddU2BI+azTQO7OBTr7JODCLmOrZNqWNpmhtmzRFjkbQxu2trWqLiaiYLM3S21bE1HLkSNx+rLun9040aViiYy+FdPRDPbdNxmq7F0bP76zLwCZ7hQxqjJnECL1xfa71zw7EtEltUz01tT9hV7icjMJ7fw742kj6w8B3hOv/Dnz5duveFaOYpWHh/n7g/oZ8z91GvR54bSR9R53Xh3wJhEb0NSHRPxC7pJ5VdpxG1Ww9smVcKG9sd5v+8UqaDd8R3tMJsWunqZ4F4kRjmal0X+IkiieUEHoh4pYx6D0V9swcYRS9Xj02owQx1O+WaBR6CAS93qShvkZgY8xfu2rpMdP/YxJEbK5nkSh0m5qItn2u57qNUWi7WhOjsLATYZ/FRMGEXWEe+e11PfOdFgvNmLD5Tse90ExjViWjZRD2B2zLiq4mdm4A1KmklSSEUTT4i2p6MG2BaJsae6brlJW+jR4iecUDyjZX16/dY/Ww6HAdVpqQMA9ZVk1mydz0+GjmEGMYMYbQ9N9ebweWITRJGfo6xgignWnYNPtMvz/27oTLghTrKSEhISGhFSl67BxCVE8WMpH9jEl9sRi2rSI9tkRvkias2sjWE5MIrNrJrpZtrKiWb2NsIaufWWjhRo+brsNugGt7h7zHxnXKskqVleeFKkvMQKKukj0Ug0GQJhpekmVQni4UxmgcnPw6etya5kP/b7vWsKLaLIPQdB1T+cQkUl1G96NJqmhrg77XfUjYM6TzKOYQnU67QXuM8n7S25HtTGu9MNR1LlbNEdMJW7VVkworxgSEMivjdZlf/yd6G9U8xGiDJuA2v7b1aGah6abeliBldX4Lu5VFD0G5FaE2bv3JStSLxlm/TrOtKsqWic1n08BtF7Nw6JiaqK3MDHPfmDZrGxN2jSRRzCHERmFtjXViGFahmkFo6cJCqNs0JlFUXpWJSROWqgoBk49lEjHC1vCDb1rsNt1rVX6seVqSkNV+PwMGnXIXtd7UKMxFu9Xqruk6Y+HRS2kiArt/kSyr1VtrtOV8diXdxCT2gpBud5Xflm870susZRIuAzzeR44d2OdI35CEhISEKwZPcTTwfOG6ZhSdTuVhYzFpblAhyJsCGLWpnHR8C13OLsslzUoFesmt05oklXA/btl8v92FZJbVz3Cyi3GREsrVfp7THww4eLDDxkYx1vb4cVFp9bMxYzrlngstYGk1VmlXsmNB3Tu1TeIoy8wyAFbVp94XHYimOmbFdtRFbZJHUjHtY8zfgUzX9bdEVE9N9D7KLKRAE7G3FUaIuCbeHXlRzPXWvkMTR8t86g0t3xGzSU47OayNuWg1nWUUNRWe6kuW9Wunca6v1/P2s8JnsJNlDAb9Gh8WZiGMo2afaIBVO2nVUy3DrKqX2Jeh8YU7xCxt2M3z7eZLuEwYkySKOYSlAZpWx+hD3zKDNvehiK1iTGeSeDdJErpROrKprjPyXnnHVOebhpVyh/EEo6m1N3If27wohfuDjMXFThkcMIrACfoDyIPhWZiDdjgo35FT65S0c2IPRSiUh6ErmWRsHmOdlGdXisDuVvqwz5MNYh/C7grd/7iuv0GdziSd0IjR0TFKqmgzhEaIeYxJ1FbvWb+QMGwj7FmkLWoTzSRizSuL2Qy2o9SZpv4vdTSFMIoliFtrrJ4SgRFmg36tefpUT5iuVqpJHEq6mpA0WtrbyiGvFHZC6G2ZxCT2GTyQjNkJCQkJCY3wJIlizmBPw9SLSHvssX7et9KDPNSruIYVa8y0UV80dsr7flYZhScM5bYSlW6FhTatyjRDrZUmBE17AqMIdgrZuNcozATpqb+8DIMO589X9dYOVoqIS/mwdovebFd7V+tgzNIZWqWZbWGW1f5O8iQpYh8jeT3NHZyb1Bjled17RmK6RTUAMSNGC5MQ/bhWN1n7teYH8cSqvlijYgzCNGMygzXItNSpH+tsraonxSTlYDmxP4xGamOjvGg4JBsckMsyjtQs9Lk0dg+rzLFyhQpx3K5yasLVJMSJCcw5kupp7iAShfz2dIRTOSYA6gFAG5lFk064wViqF8VyLe8uecIwwiQELbprXbdu+4SNwv5vIUIxxy65bvSiMvXJrZZGhAn3JUNofCffpNfrc/FikU+H+7Bo89LS0kbJfNuIbZMx+0ohMYJrHEn1lJCQkJDQiqR6mmvYRaMcvhM7Ba8MFtim79ZL7og0oesTxyZ5b5ZF3DhtAfuuCOzxCBNlmgwQeQ5Zv3bbhKgHUmhP00o/y6qD5GqJSqIQu4Z+fxkIcDjZID3EtTYptVdpu9juqj2t8hP2DEmiuKagD8CZiB1k0cAQmtIsjdberxPqFUtAdwCt0io0L53ph6U31KGbpXdgTyDCFFXyxHXNUF8yinrXJ5hSyDDhyRpxjYUW5puQcMWQbBTXHOwJlROwDCFGBVukDrFLSJaFBZMtdqjOLjDVJV8ZtTtZRpZ1Zn/tlH7qR3aYait9JYlpRlEyKhuIMWYHUS+MNavVppGQcNmRJIq5hXVgago/PrEKhjgFbqKwSq2T55XkkmV1yWUzr44tjW7wm6hzl1OpPZ9a6pIIsLVyLeqwJjfdqDFcdFLhOsuqedjaCp5KDe7AoJi6emmNUVxJo3RCQhTJRpGQkJCQ0Ir5VD0lGVxBVvVQnapmn0NQXcT0J01ShbExiF1C2w60bWJikd5kLG/ph4W8q6w31t4pEOmnUROmG25W9baM2Dg6+WazqEF1RrbUA0zGvKI+lrWXBMTsE0kFlXDlIcbsps/+RPqlGGhmIfcxB6aSWcSYhKWKhkLqx93uFCbR1MgYggG4KctUzYsh9DuKMht5SYx/1AzhEQO13mzR7UZUgUFFJrGtRG1nPZ4m9o7EmEdCwhWDSBRNn53DOfcK59wjzrmxc+6ulnx3O+c+4Zx71Dn3hlnq3hWjcM49zTn3O865Pw//jzTke3XI8+fOuVer9Bc45z4SGvyTzhVb4Jrqdc59mXPufc65Defc6807tt15qNMpS/ebiG6NFjatzGPMI8vYzCsDca9XrZgto4iiTQqIFGo1xM8I+7qtpkWPbriRnqAu0UyUMZMwzvolAxDGXfZF2WpiUXJjHk/JPpGwf3BZJYqPAt8E/EFTBudcF/hp4KXAHcC3OufumFbxbiWKNwDv9d7fDrw33NuGPQ34fuBFwAuB71cM5WeB1wC3h8/dU+p9Evhu4EfNO3bUeXt+rWYMTYKCxrhQnjRniHAhfTszg9CNs9cGHcZTNUp5biQiW6dpjK1Pe2pFuajpyGhUMJipzIKK+FuX4ZrbsHpnm2tsjIEkdVPC1YUYs5s+u6jZ+4977z8xJdsLgUe994957zeBB4B7ptW921/NPcDbwvXbgL8dyfMS4He89096788DvwPc7Zy7CTjkvX+f994Db1flo/V675/w3n+ASRltR51PSEhIuPK4qjaKm4HPqvtTIa0Vu/V6eob3/nEA7/3jzrmnb6NhN4drmz5rvdPe8aKZe6Ewq25/wj/fLm3t/WDAmM7EIr5JktCSRt+upGdwh5W6e736Zrlo0Sa32DwnRGGKIs9VjKbIw8IGUK1FRG1Vmh8aOh9TBwJVED8jUYxGhQovy6h2bWfZZNW7dSFOSNg1xkyxRRxzzj2s7u/z3t8nN865/wIcj5T7Xu/9b87QABdJ89MKTf3ltDVshkZBc8N21OBtvmMyo3OvoVB3cdttt9VoutBKa7jWiKmqajcRnbuutDS40h7iQmMiuqplRJHy0hdhFo0ZZ6hLG/i151P0pVP6ImXLA6B0ZXkOWX+Cf2SZCblu1HiCmnE6K4IB6vFOjCLh6mPqPoqz3vtGQ7T3/ut22YBTwK3q/hbg89MKTf3ltDXMOfeXzrmbwqr/JuCJhoZ9lWnY74f0WxoaPEu99h0zdT5w5/sA7nrBC7xdpFtTgF5oTxhOp0EMEUGaKAmWWb3PXKdtbItkYZmFRaOkMgUT3kczMLBer2IwNWN4S4cn4lRpw4WaqDwv6iykier5hO3IvCt5PCVcPVxVN9gPALc7554FfA54JfCqaYV2a6N4CBAvplcDMdHnPcDXO+eOBCP21wPvCaqlVefci4O307ep8rPUq1F23jnXp+j8QzvtVEJCQsLlwWV1j/07zrlTwF8D/pNz7j0h/Qucc+8C8N7nwOso6PLHgQe9949Mq3u3svhbgAedc/cCnwFeERp2F/Cd3vvv8N4/6Zz7QQpiDvAm7/2T4fq7gF8GFoHfDp+2eo8DDwOHgLFz7v8A7vDeP+Wck853gftn6TzeBy+h+kFCMQejmbwrdSZx0wkK+bZznu2ifCqmSBOxbDGJqfHlZq9BllXhRBo3tsWkibyKAGuSJxFxqd3aMmozo4+aqEfqGAwmPZ6CVJeQcHVx+aLHeu9/Hfj1SPrngW9Q9+8C3rWdunfFKLz354CvjaQ/DHyHur8fuL8h33O3Ue9p6uoq/WzbnQdgOKSTZfSzLO6XH2B15hMEN2K8FhXIRJTVCIG3SY08IEb9WzBTtm1ww5mDJJpkG7a9NIQbCGMShlSqoGKuxiY+YLNrbHufEhKuHFKsp/mDHHEXzuXsDAboiKkxuhk7KzrKRdqYhGJKkjQT2piEqaTNQSrK6KY0pNXhSksUDe2ABvvGDHaNwp4SMRQRkXBCG6whOxpuPCHhiiOdRzF/cMpZKhAgHV475uFkNTU1lZKSImqIMAldp86iMZMUMAWWDs/MIEKa7UtEO9X8MpUsbroTZ3tEuFajI1aLh1aeQ189K3nWDFJXQsKVw3wGBUy/oISEhIQrhqR6mk9olYZSQRFRDen/E66uRp0Eyve/dRk+9VFze2cs3KKhqj9oaWur+qpJ9DKFZoo9FTb5iV2jTBoY8Svki/bB1hfVFyYkXC3Mn+ozMYqY7jt461jiKHS0ZncIBFEH+5ugkQ1qnF3pytULZt24N/FOS1RVw6XOmL1GsjV6cql67BYLZcKpjMyGwMfofbnpUM1XlvVLW0aeT3a4uE1f8YT9hDFJophnGINBBxgM+vHdv2bjlzAJK4EUgffiITBaCXasTag6FWYpWpN+mjADg5BsE8yyIaP2Imt1qGrxHhAmUHpJaWY+qJhKngPLA1hbg+GQLDtQ61eHcX3TY0LCVUMyZickJCQkNCIZs+cTMf166QGV0w8btWpeNCq/3Xthq5T7RjNF23LbFJy22tfdiKrJYoWMuqlNApCmlHXKSXMRvdSYTil45Xnl8RRtV4OfsEgTw2FhZujXxAfoLy8zDG1eW4OnHQ7lh0MGywfKPokUkaSJhKuP5B4737BGYkWQOgiRUvmDK2zTXoih3QwWyRMl4E0bH0yRNrpfK6ILTNFT2ayyR0G7tJbMUnOBlraur9fjO03stNYHT4T6cvrxfkrVNQZV5d3MO4V6ajikf3jMpWHBWPtmOJMKKuHqIXk9zSe05TQmXWiKn2W1IH+WSUR4TO15I2LMwrYNJuwdrXp/GuwgUxpmmQQ0hO3QmSP1TWtbIyxDniG7vG84DAx9OAzzdqCwbUTGPzGLhKuDxCjmEhMhvKFyp2xYtm93V3WrS+osBY0bVVuxmd81oxtsY9aa+1IkLQhG+hjTiXqE8cpghnGX5NGI0k02y9QDVUmWwfKyEsRmdIVNTCLh6sCT3GMTEhISElqQJIr5RmzZbNVSDSt7WZ3qKLS2mtb3autzU7uoAubZprUUqdc/JfPMWeViebkx0wxZKnHAPJTHR45YYSPklzIh7+HDqt7BcilVDNK3O2HfYT6N2a44rvr6hHNuFZh2GPn1gGPA2avdiKuMNAYF0jg0j8EXeu9v3E3Fzrl3h/qbcNZ7f/du3nE5cL0ziofbjh28XpDGIY2BII1DGoMY0kkuCQkJCQmtSIwiISEhIaEV1zujuO9qN2CfII1DGgNBGoc0BhO4rm0UCQkJCQnTcb1LFAkJCQkJUzCXjMI5d7dz7hPOuUedc2+IPF9wzv1qeP7HzrlnqmdvDOmfcM69RKXf75x7wjn30YZ3vt45551zx1TaVznnPuyce8Q591/3tpft2A9j4Jz7v0L/P+yc+6hzbss597S9720c+2QMbnDO/ZZz7mT4HvyDve9pO/bJOBxxzv26c+5PnXPvd849d+972owrOQbOuR9wzn1Offe/YVpdcw/v/Vx9gC7wF8AXURxRcBK4w+T5h8DPhetXAr8aru8I+ReAZ4V6uuHZ3wCeD3w08s5bgfcAnwaOhbTDwMeA28L906+3MTDPvxH43ettDIB/Brw1XN8IPAn0r8Nx+NfA94frLwPee62OAfADwOsj7Wisa94/8yhRvBB41Hv/mPd+E3gAuMfkuQd4W7h+B/C1zjkX0h/w3m947z8JPBrqw3v/BxQ/8hh+HPi/KbZVCl4F/Afv/WdC+Sd23bPZsV/GQONbgV/ZYX92gv0yBh44GOpdDmV3GhJxJ9gv43AH8N5Q9n8Az3TOPWO3nZsRV2MMYmisa94xj4ziZuCz6v5USIvm8d7nwAXg6Ixla3DOvQz4nPf+pHn0JcAR59zvO+c+6Jz7tu12ZBfYL2Mgzw8AdwPvnL0Lu8Z+GYOfAv4K8HngI8A/9t5fyahv+2UcTgLfFPK8EPhC4JbtdGQXuKJjEPC6oGa73zl3ZBvtmEvMI6NwkTS7ym3KM0vZqpKCAH4v8H2RxxnwAuBvAS8B/oVz7kua6tpj7JcxEHwj8N+899tZfe0W+2UMXgJ8GPgC4E7gp5xzh5rqugzYL+PwFoqF04eBfwT8CVdOsrpiYxDws8AXU8z348C/2UY75hLzyChOUehIBbdQrOaieZxzGXADhQg5S1mNL6bQNZ507v9v3/5ZGgbCOI5/Aw6ii07qJm5SQVdfgg7ODiIWcXLRwamL+CYc3HQQFBycxNHRwUE6+I+2i4PgJDiJPA53samEk3RIm+b3gcAlba7N0wtPck8atfz776IomvR9XZnZp5m9AzfAfJfHlFW/xCC2Sr7TTtA/MajipiDNzF6AJm6OPi99EQcz+zCzqpktAOu4ek2z24PKKM8YYGZvZvbt7xyPaE8vZe6rMHpdJMm64K7kG7gBGxeuKn/es01n4erMtyt0FpsaJIpNwDQpxbvE6y3axbtZ3JzsEDAC1IG5MsXAr8cn3GhJx8EhsO/bE8ArKcX+EsRhDF/EB7aA40GNATCVaO/i6hL/9lXkpedfoMuBsQw84Z4qqPltB8CKbw8D57hi0i0wk9i35vd7BJYS209xt5FfuCuDzZTP/T0x/Poe7smnOrBT0hhsxCdKGccBbsrpGlefmPJcQwAAAFhJREFUqANrJY3DIvAMPAAXwPigxgA48b/3PXBJZ+JI7avoi/6ZLSIiQUWsUYiISI6UKEREJEiJQkREgpQoREQkSIlCRESClChERCRIiUJERIKUKEREJOgHfjMksYZVbXAAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x216 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"species = 'beam'\n",
"iteration = 3\n",
"field = 'Ex'\n",
"\n",
"snapshot = './lab_frame_data/' + 'snapshot' + str(iteration).zfill(5)\n",
"header = './lab_frame_data/Header'\n",
"allrd, info = read_raw_data.read_lab_snapshot(snapshot, header) # Read field data\n",
"F = allrd[field]\n",
"print( \"Available info: \", *list(info.keys()) )\n",
"print(\"Available fields: \", info['field_names'])\n",
"nx = info['nx']\n",
"nz = info['nz']\n",
"x = info['x']\n",
"z = info['z']\n",
"xbo = get_particle_field(snapshot, species, 'x') # Read particle data\n",
"ybo = get_particle_field(snapshot, species, 'y')\n",
"zbo = get_particle_field(snapshot, species, 'z')\n",
"uzbo = get_particle_field(snapshot, species, 'uz')\n",
"\n",
"plt.figure(figsize=(6, 3))\n",
"extent = np.array([info['zmin'], info['zmax'], info['xmin'], info['xmax']])\n",
"plt.imshow(F, aspect='auto', extent=extent, cmap='seismic')\n",
"plt.colorbar()\n",
"plt.plot(zbo, xbo, 'g.', markersize=1.)\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Read back-transformed data with hdf5 format (example: 3D run)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"import h5py\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<Figure size 432x288 with 0 Axes>"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"text/plain": [
"<Figure size 432x288 with 0 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"f = h5py.File('HDF5_lab_frame_data/snapshot00003', 'r')\n",
"print( list(f.keys())\n",
"# plt.figure()\n",
"# plt.imshow(f['Ey'][:,,:])"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"anaconda-cloud": {},
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.1"
},
"widgets": {
"state": {
"11d243e9f5074fe1b115949d174d59de": {
"views": [
{
"cell_index": 6
}
]
}
},
"version": "1.2.0"
}
},
"nbformat": 4,
"nbformat_minor": 1
}
|