{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## using obspy + yt to extract properties along a seismic ray path \n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import yt\n",
"import xarray as xr\n",
"import numpy as np\n",
"import pandas as pd \n",
"import matplotlib.pyplot as plt\n",
"from obspy.taup import TauPyModel\n",
"from obspy import UTCDateTime\n",
"from obspy.clients.fdsn import Client\n",
"from pygeodesy.sphericalNvector import LatLon\n",
"import cartopy.crs as ccrs\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## load the data into yt"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"yt : [INFO ] 2021-01-15 15:03:21,902 Parameters: current_time = 0.0\n",
"yt : [INFO ] 2021-01-15 15:03:21,902 Parameters: domain_dimensions = [112 181 361]\n",
"yt : [INFO ] 2021-01-15 15:03:21,903 Parameters: domain_left_edge = [ 0.00000000e+00 -8.99899979e+01 1.00000005e-03]\n",
"yt : [INFO ] 2021-01-15 15:03:21,903 Parameters: domain_right_edge = [ 2900. 89.98999786 359.99899292]\n",
"yt : [INFO ] 2021-01-15 15:03:21,903 Parameters: cosmological_simulation = 0.0\n"
]
}
],
"source": [
"\n",
"ddir = \"/home/chavlin/hdd/data/yt_data/IRIS_models/\"\n",
"modelfile='GYPSUMS_kmps.nc'\n",
"xr_ds = xr.open_dataset(ddir+modelfile)\n",
"bbox = np.array([\n",
" (xr_ds[dim].values[:].min(),xr_ds[dim].values[:].max()) for dim in xr_ds.vs.dims\n",
"])\n",
"bbox[1][0]=-89.99\n",
"bbox[1][1]=89.99\n",
"bbox[2][0]=0.001\n",
"bbox[2][1]=360 - 0.001\n",
"\n",
"data = {'vs':xr_ds.vs.values}\n",
"ds = yt.load_uniform_grid(data,data['vs'].shape,1.0,bbox = bbox, \n",
" geometry=(\"internal_geographic\",\n",
" ['depth','latitude','longitude']\n",
" )\n",
" )"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"yt : [INFO ] 2021-01-15 15:03:21,909 Setting origin='native' for internal_geographic geometry.\n",
"yt : [INFO ] 2021-01-15 15:03:22,040 xlim = 0.001000 359.998993\n",
"yt : [INFO ] 2021-01-15 15:03:22,040 ylim = -89.989998 89.989998\n",
"yt : [INFO ] 2021-01-15 15:03:22,041 xlim = 0.001000 359.998993\n",
"yt : [INFO ] 2021-01-15 15:03:22,041 ylim = -89.989998 89.989998\n",
"yt : [INFO ] 2021-01-15 15:03:22,043 Making a fixed resolution buffer of (('stream', 'vs')) 800 by 800\n",
"/home/chavlin/miniconda3/envs/yt_dev/lib/python3.7/site-packages/cartopy/mpl/geoaxes.py:388: MatplotlibDeprecationWarning: \n",
"The 'inframe' parameter of draw() was deprecated in Matplotlib 3.3 and will be removed two minor releases later. Use Axes.redraw_in_frame() instead. If any parameter follows 'inframe', they should be passed as keyword, not positionally.\n",
" inframe=inframe)\n"
]
},
{
"data": {
"text/html": [
"
"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"\n",
"target_depth = 150\n",
"center = ds.domain_center\n",
"p = yt.SlicePlot(ds,\"depth\",'vs',center=[target_depth,center[1],center[2]]) \n",
"p.set_cmap('vs', 'dusk_r')\n",
"p.set_log('vs', False)\n",
" \n",
"p.set_mpl_projection(('Mollweide', (), {'central_longitude':-110}))#{'central_longitude':-100, 'central_latitude':40}))\n",
"#p.set_mpl_projection(('Mollweide', (), {'central_longitude':0}))#{'central_longitude':-100, 'central_latitude':40}))\n",
"\n",
"p._setup_plots()\n",
"\n",
"p.plots['vs'].axes.coastlines()\n",
"p.show()\n",
"\n",
"\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## use obspy + pygeodesy to get event info and generate a ray path"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[Event:\t2011-03-11T05:46:23.200000Z | +38.296, +142.498 | 9.1 MW\n",
"\n",
"\t resource_id: ResourceIdentifier(id=\"smi:service.iris.edu/fdsnws/event/1/query?eventid=3279407\")\n",
"\t event_type: 'earthquake'\n",
"\t preferred_origin_id: ResourceIdentifier(id=\"smi:service.iris.edu/fdsnws/event/1/query?originid=18196524\")\n",
"\t preferred_magnitude_id: ResourceIdentifier(id=\"smi:service.iris.edu/fdsnws/event/1/query?magnitudeid=176443936\")\n",
"\t ---------\n",
"\t event_descriptions: 1 Elements\n",
"\t origins: 1 Elements\n",
"\t magnitudes: 1 Elements]"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"client = Client()\n",
"cat = client.get_events(eventid=3279407) # Tohoku\n",
"cat.events"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[197.0, 38.2963, 142.498]"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"event = cat.events[0]\n",
"depth_km = event.origins[0].depth / 100\n",
"lat = event.origins[0].latitude\n",
"lon = event.origins[0].longitude\n",
"[depth_km,lat,lon]"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Station AIN (Ainapo Ranch, Hawaii)\n",
"\tStation Code: AIN\n",
"\tChannel Count: 0/7 (Selected/Total)\n",
"\t2007-04-01T12:00:00.000000Z - 2011-08-12T00:00:00.000000Z\n",
"\tAccess: open \n",
"\tLatitude: 19.37, Longitude: -155.46, Elevation: 1524.0 m\n",
"\tAvailable Channels:\n"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"stat = client.get_stations(network='HV',station='AIN')[0][0]\n",
"stat"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"56.55070260333306"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"eq_start = LatLon(lat,lon)\n",
"station_loc = LatLon(stat.latitude,stat.longitude)\n",
"gc_dist_km = eq_start.distanceTo(station_loc) / 1000\n",
"gc_dist_deg = gc_dist_km /111. # approx 111 km /deg\n",
"gc_dist_deg"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [],
"source": [
"model = TauPyModel('iasp91')\n",
"arrivals = model.get_ray_paths(source_depth_in_km=depth_km,\n",
" distance_in_degree=gc_dist_deg,\n",
" phase_list=[\"P\"])"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAOcAAADuCAYAAADC+xDoAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAABidElEQVR4nO39eXxbZ5n3j39ubZa8Sd4keZH3NU5iJ3H2pEBL92laljJAS1sedgoUmIGBMlC2mTIPPFDgNy10gGlLCwMDv7ZpS2kBZ3dsx4nt2Il32dZuybIlWbt0zv39wz0HyZZjJ7HjI0fv10uvRLJ9zi3pfM513dd13ddNKKVIkSKF8BCt9wBSpEiRmJQ4U6QQKClxpkghUFLiTJFCoKTEmSKFQEmJcwNx6NCh+kOHDtWv9zhSrA4klUrZGBw6dEgsItQBAsqypODw4cPseo8pxdWRspwbhxaWkhyWJbkAWtZ7MCmunpQ4Nw63AxTzD9y+zmNJsQqkxLlBIITeUVPlo9WVPpYQesd6jyfF1ZMS5wbg0KFDBZSipaVplrQ0u0SUYuehQ4fy13tcKa6OlDg3AJTSWwBCtje7sL3JBYCQ+ddSJDOS9R5AingIIQSAEkDhW4+i4uLirZs2bXqIUqpgWVZGKRUTQvgbq0wmIyoVUF3hAwWQkR6By40Xbrrppuff+hVKKaUikShKCAkRQgK9vb1P2+32iwAsAKwArJRS77V9tykuRSqVso4QQjIANGdmZu4tLCy8OxQKbRGJRIqsrCxRbm4uycnJEeXk5BC5XA6TyQQAyM6KoLQksOhYN79jCm/f7wQAHD2Vh78c0Sz6HYNJAc+cFABQVFSEYDCI2dlZOjMzw87MzFCfz8cA8Eml0jMGg+HlcDjcCaCfUhpaq88gxdKkxHmNIIRkAtj2lhDvCQaDjenp6YrKykpxVVWVqK6uDps3b0Z+fj4yMzORmZkJuVwOuVwOiUSC/v5+/OTHP4JtyoE7b7HigX80Qi5fWSozGBThud+V4rU3tdBqCvC5R76AxsZGRKNRBINBBAIBzM3NwefzwWq1or+/H8PDwxgbG2MnJiaYSCQyJ5PJzppMppeDwWAHgL6UYNeelDjXCEKIBMC+srKyj1JKDykUivSqqipJVVUV4YRYWFiIgoICKJVKiESXnv5TSuH1evHss8/izTffhEYdxuc+PoLNDXOX/Lv+gSz85OkaTNlluPXWW/Hggw8iMzNz2fFHo1G4XC7Y7XaYzWb09fVhZGQEo6Oj7MTEBBONRt2hUOh5q9X6HIAemrqQVp2UOFcRQkg2IeTWqqqqR0KhUEtjY6Ns165d5ODBgygtLYVarUZ2dnZCIbIsC6/XC7fbjWAwGPdgGAYAIJFIIJfLYbFY8OKLf4DT6cLDHx3DLe9wJBzPm61q/OcvK5GXp8I997yHd2W540mlUsjlcqSlpUEul0OhUECpVCIjIwPzU994otEo3G43pqamMDY2huPHj+PMmTPs2NhYUCaT/U2v1z8J4EjKqq4OKXFeJYSQMqVSeW9eXt6nxWJxyY4dO6R79uzB3r17UVVVhdzc3EUXOsuymJubg8vlgtvthtvtBsMwyMrKQnZ2NhQKBe/Scm7tQvx+Pz70oftxcI8Vj3xyLOHYnvhZFU62F+LXv34e6enpcT+jlPJuLffw+/1wu93w+XyQyWRQKpVQqVRQqVQJBcuyLKampjA6OoqTJ0+io6OD9vT0RKRS6aDBYPhJOBw+TClNfOdIsSypaO0VQAgpyM/Pfzg9Pf1z27Zty969e7d4//792LZtG8rKyha5jYFAAHa7HS6XCx6PhxeiSqVCUVERGhoaIJVKL2sMJpMJkUgUzVtcS/5O8xYXjpwogNlsRk1NzcL3AKlUCqlUiqysrEV/GwqF4Ha74XK5YLFYFglWrVZDJpOhsLAQhYWFOHDgANxuNxkfH5edOXNma1tb2y/OnDkTraystNvt9u/4fL7nU9HgyyNlOVfIW6mLG6uqqh4XiUTNN910k+SOO+7A5s2bUVJSEicuSincbjdsNhvsdjskEgnUajVycnKgVCoTWsJYotEoQqHQIveWey0cDqOtrQ0dHafx3FNnkZ0VRTAowh8OFwEA3nvIArmchWdOggc+tQN79uzD3r17efc11pWNfS4Wiy85rnA4DJfLhdnZWdjtdhBCoNVqodVqF92QgsEgDAYDzp07h9dff50eO3YsIhKJ/jo+Pv4NSunZK/sWri9S4lwGQoiqoKDgEblc/s+NjY0Zt912G7nppptQU1ODtLQ0/vcYhsH09DRsNhtmZmaQnZ0NrVYLtVp9SasYCAR4C+VyuRAIBCCRSBaJJ1ZQMpkM//zP/wTC9uL73+5D/0AWfvp0LWz2+fPEBov+6etbIJI04fvf/39xgk8kfpZlkZmZCZVKxVtImUy25NiDwSCmpqZgs9kQCARQUFAArVa7yJX3er0YGBjAG2+8gTfeeIM1Go0zTqfza16v9zlKaXAVvqYNSUqcS0AIqSsvL/8eIeTOm2++WXrPPfdgx44dKCgo4C+8aDQKi8UCm80Gv9+P/Px8/uJMFPQJBAL8PJMTolwuj5vbKRSKhMGYWNxuNx544EO45w4LwhFRXIoEQFzKRSalePG1Ijz//PPIzs6+5HG5iDA3PrfbjXA4jIyMDH58SwmWYRg4HA7YbDbMzs5CqVSisLAQGo2G/yxYloXZbEZ7ezteeukleuLEiRDLsr8ym83/Rim1rOiLuY5IiXMBhJDdFRUVz6hUqpq77rpL/A//8A/YvHkzFAoFgPkL2OVyYXJyEi6XC4WFhSgqKkJmZmbCgMns7CxsNhscDgfS0tLirNJKhJiIo0eP4oc//CFkUhbhiAh33XUXPvShD0EulwOYt2jPPfccXn31VUilLCIREfbe8yC++n/ec9nnopTC5/Pxlt3lcoFlWWg0Gmi1WmRnZy96D9xnZLFYYLfboVarF83FPR4Puru78fLLL+NPf/pTNBwOt42Pj3+UUjpy2YPcoKTE+RaEkPry8vLnCwoKtj3wwAOi2267DRUVFfw8LBqNwmg0wmg0Ij09HWVlZcjPz190YUYiETgcDlitVng8HuTk5ECr1aKgoGDZOV0iIpHIIhf0t7/9Lc6cOYOcnGzcfvs/QKfTgVIKp3O+QigvLw+EEBiNRvz/X3kFwbk5zOTU48P3fwC3bVLHucpXMqZwOAy73Q6bzYa5uTnk5uaisLAQeXl5i47HsixsNhsmJyfBsizKy8tRWFjIW9NIJILBwUEcPnwYL7zwAuPz+V41GAyfopRaL3tgG4zrXpyEkJKysrKfZ2Zm3vqhD31I/O53vxvV1dW86Px+P/R6PRwOB0pKSlBaWho31wTm3VWr1YqpqSmEw2Go1WpotVqoVKoVWUbOnYxNrUQiEQB/z23GzkGj0Sg6Ojpw2223QaFQQCQSJbReLMvC7/fjG0/9Fn+a1YCVyPGZvQW4uVweN88EgLS0tMt2r4F58TmdTthsNkxPTyMjI4MPEi10f30+HyYmJmC321FcXIyysjL+s4xGo7hw4QJ+97vf4Xe/+10kFAr90mw2f4VS6l52EBuU61achJBcnU73/2Qy2f0f+MAHJO973/uwadMm/s4/MzOD0dFRhMNhVFRUxN3tgcUWobCwEFqtdlE+cSGxQuTEGI1G+XndSgIxlwulFE8eHcP33xwCADxyUw0+f1MNLz5KKYLBYNxc0+/38274SgVLKcXc3BxsNhusVisyMjISehjRaBQmkwkTExNQKpWoqanhXd5wOIyzZ8/iN7/5DQ4fPhwKBoPfs9vt37seA0fXnTgJIelarfZf09LS/umee+6RffCDH8S2bdv4iOrc3BwuXrwIAKitrUVOTk7c38fe/RPNpRIRG9X0+/18RJQT42oK8VL8zxkDHn2xDywFPrBTh2/fvRlS8dJlg8FgMC6A5ff7VxyFXjg3LyoqirOU3O84HA4MDQ0hKysLdXV1/Nze7/ejvb0dzz33HFpbW/1ut/uLHo/nF5RSZvU+EWFz3YiTEEKUSuVD2dnZ/3nLLbco7rvvPuzZs4e/GAKBAAYHB+Hz+dDQ0IC8vDz+bymlmJ6ehl6vRzQaRUVFBbRa7ZL1sJwFsVqtsNvtEIlEfABlJXWta8kbF2z47P90Ixxlsb8qD0/etwNKxcoKIBLlb7n3lZGRseTfRSIRmM1mTE5OIisrC1VVVVAqlXHHtVqtGB4ehlqtRk1NDS98t9uNY8eO4bnnnkNnZ+es1Wq9NxKJ/O3qPoXk4LoQJyGksLS09LX6+vrmz3zmM+Ttb387XxUTDocxPDwMp9OJuro6aDQa3gVjWRZGoxETExPIzs5GZWVl3EUVC8uyfJ7T6XQiMzMTWq0WGo3mmlnGlXLOMIuP//ospr0hVOZn4JcP7kRF/tLiWopAIACbzYapqSmEQiE+z5mTk5PQ/V14k6usrIRWq130eev1epSUlKCyspKfZkxPT+PPf/4zfvzjH7N2u/0lg8Hw4EavONrQ4uSspUqleur//J//k/bAAw+gvLwcwPy8R6/Xw2w2o6qqCjqdLm4OZjabMTo6Cq1Wi4qKikVBIA6fz4fJyUlMTU0hLy8PWq0W+fn5y64yWW/MrgA++lwXBqweKBVSPHXfduyruvLOJtFolI/gejweFBcXQ6fT8emdhfh8PoyOjsLj8aC+vh4FBQX8zxiGwfj4OIxGIyorK6HT6SASiUApxYULF/D000/jpZdemjObzYcYhjl6xYMWOBtWnIQQbWlp6Wt1dXXbPve5z5F3vvOdkMvloJTCYDBAr9ejtLQU5eXl/N2ZUgq73Y6hoSHk5OSgtrY2oShZloXVasXk5CQAoKysDFqt9orSEuuJLxTFI7/rwV8HpiAREXztjgY8tK/8inKvsXBurNFoRFpaGsrLy+OKN+LG4PNhYGAA4XAYmzZtgkqlijvOyMgI7HY7GhoaoNHMLyD3eDw4fPgwnnjiCdbhcLxoMBge2ohWdMOJkxBCsrKyHsjNzf35hz/84bQHH3yQt5aBQAA9PT3IyMhYVGw+MzODgYEBKBQK1NfXJ4y6hkIhTE5Owmw2Q61Wo7y8/JJzrWSAYSn+7xuD+PlxPQDgzi2F+I/3bEVm2uqsiXC73RgfH4fL5UJZWRl0Ol3C2mK3242LFy9CIpGgoaEhbm4eCARw/vx5yGQybN68GVKplLeiP//5z/HSSy/NWSyWuxiGObYqgxYIG0qchBBNWVnZazU1NdsfeeSRhNZy8+bNcS7U3NwcBgYGQClFQ0NDwhI3r9eLkZEReDwelJWVoaSkZNni9WTjtT4rvvyHXvjCDKoKMvCz+3agRrN4tcqVEg6HMTExwd/YqqurE3ol09PTGBgYQHZ2Nurq6ni3mFIKk8mE0dFRNDY2Qq1WA/i7Ff3Rj37ETk9P/9FgMHyYUupbtYGvIxtGnOnp6YfUavXvH3roobSHHnoozlr29vZCoVCgsbGRF1UkEsGFCxfg9XoXRWc5gsEghoaG4PF4UFdXt6RrtlEYc3jxqRfOYnjKi3SZGI+/awvubi5e1XOwLAuLxcLP56uqqhalZGKjt1qtFjU1NfyUIRgMore3F2lpaWhsbOStaH9/PzcX9ZhMprdTSrtXdeDrQNKLkxBCCgsL/02j0fzLN77xDdHtt98eZy3Hx8fR2NgYZy2npqZw8eJF1NTUoLi4OGEJ3ujoKKamplBbW4vCwsINLcpY/OEoHn2xDy/1zNehv2d7Cb51qHHV3FwOlmUxOTmJiYmJRXP/2N/hgnZNTU38fPRSVvQPf/gDHn/88ajdbn/I7Xa/sKqDvsYktTgJIYrS0tLXt27desO//uu/kl27doEQEmctN23axN+ZI5EI+vr6EI1GsXXr1kWRxKWihNcblFL8ptOA77x2EcEIi7LcdDzxj83YVpqz/B9fJtFoFGNjY7BaraiqqkJJScmiG+Hc3Bx6e3uRl5eHuro6/jtJZEWj0SjefPNNfOtb32LNZvNTZrP5c5TSpNzUKWnFSQgpLikp6bznnnuKvvCFL6CyshIAYLPZMDAwsGhueSlrGTsnXZhfu54Ztc/hs//TgwGrB2IRwedurMGn3151yaqiK+VS+WZg/jsaGxtb0oqOjY2hubkZKpUKlFJ0d3fj8ccfR2dnZ5fBYHhHMkZzk1KchJDdOp3uyMMPP6z4yEc+gvz8fFBKMTIygunpabS0tPCJ/0gkgv7+fkQikYTW0u1283flmpoawRUMrDehKIPvvzGEX5wcBwA0FmXj++9twqbCS68NvVL8fj8GBwcRCoXQ1NS0KGrOWdH8/HzU1tbyVtTn86Grq4u3vgBgMBjwk5/8BL/73e8cJpNpN6V0fE0GvUYknThzcnI+plarn/ra174mvvfee6FQKBCNRtHT08O7N9wXdilrybIshoeH4XA40NTUtOxC5OudtrFpfPmP52GaDUAiInj4HdV4+O3VkEnWxu2fnp5Gf38/KioqUFpausiKjo6OwmKxxFnRSCSCc+fOISsrCw0NDSCEYHZ2Fs888wyeeOKJkMViuT0SiRxZkwGvAUkjTkKIuKSk5JelpaUPPPbYY+Smm26CWCxGIBDAmTNnUFZWhrKyMgDzX97AwAA8Hg+am5uXtJZarRbV1dXX5bzySvCFovi/bwzi2dPzxRf12iz84L1N2FycuKTxaolGo7h48SJ8Ph+am5v5OmiOubk5dHd3J/zu5+bmsH37dkilUoRCIbz44ov4zne+wzgcjq/Y7fYfrMmAV5mkECchJE2n0504cODAzq985SvYsmULCCGYmZlBb28vtm7dyqdCuLtndnY26uvrF1lLruIkZS2vnI5xJ778h/OYnPFDRIAH9pTji7fUIlt+eR0EV4rD4cCFCxcSWlGGYdDd3c0XKHA3Wm4e2tLSgoyMDLAsixMnTuCb3/wmHRsbe85oNH5Y6I2wBS/OtyKyHbfffvuWr33ta9DpdADm5xMTExNoaWnh5yXcvKO6uhrFxfH5uZS1XF0CYQb/7y9D+O+2CTAsRX5mGh69vR7v2rY4NbUacIuxA4EAmpqa4qzoUvEGl8uF7u7uuOBgX18fvvrVr6K3t/f3JpPpA0KO5ApanISQDJ1Od/buu++u+8pXvoLi4mK+bCsQCGDbtm18UYHD4UB/fz+2bdsWV5/JzU9sNlvKWq4BA1YPvnG4H2cmZgEAu8pz8e27G1GvXZvPmbOiXBwhFpvNhsHBQWzfvp3/noPBILq6ulBUVMRH9AcGBvCv//qvOHPmzKtGo/Eeoa4RFaw4CSFZOp2u973vfW/Fl7/8ZWi1WrAsi+7ubqSnp/MuK6UU4+PjsFgsaGlpiZtfci7PwkBRitWFUoo/njPj8dcH4PSFISLA+3eW4gvvrEVBVuLVPFfDpaYuHo8H586dQ319PbRaLYC/XweZmZmoq6sDIQRDQ0P45je/iVOnTrUajcZbKaXRVR/oVSJIcb5lMfs+8IEPVHzxi1+ERqMBy7Lo6upCTk4O372cZVmcP38eLMuiqakpLjcZCATQ1dWF0tJSPliQYm1xByL44V+G8HyHAQxLkSET45Nvq8JHD1RCIVvdvHGiwA9HOBxGV1cXCgoK+H5QlFI+os9FcvV6PR577DF67NixI0aj8RahWVDBiZMQotDpdD333ntv7b/8y79ArVaDYRh0dXUhPz8fVVVVAObvnp2dndBqtaisrIy7eyYKFKW4dozavfjenwfx14EpAIAmOw3/fEsd3r2tBGLR6s5HFwZ+OFiWRX9/P6LRKJqbm/n1oOfPn4dIJMLmzZt5gT766KNoa2t71Wg03i2kOaigxPlWVPbsPffc0/joo49Cq9WCYRh0dnaisLCQL2YPh8Po6OhAVVUVioqK4o6RKFCUYn04rXfi3/80gD7zfAO9em0WvvDOWtyySbOqQaPZ2Vn09PQsqgoDgOHhYbjdbuzYsYMXaH9/P1iWxdatW3kX99FHH0VnZ+fvTSbT+4USxRWMOAkh0tLS0vY77rhj+9e//nUUFRWBZVneOnLCDIVC6OjoQG1tLT+nALBkoCjF+sKyFId7Lfj+m0Mwu+Z35N5clI0vvLMWN9arV02kXOCnuLgYFRUVcT8bGxvjI7lisZgXKADegg4MDODLX/4yent7nzUYDA+tyqCuEsGIs7S09H/e9ra3ve/f//3fiU6n4+eYeXl5vCsbDAbR0dGBhoYGfiUCMB9m5+ajtbW1180KkmQiGGHwP2cM+M+jY3DMzW/fubVYiU+/oxq3NGggWgV3l2EY9Pb2QiwW81aRY3x8HDabDbt27eIFev78eUilUn4O2tfXh8997nP0woULXxZCoYIgxJmfn/+5xsbGJ5566imyadMmUEpx9uxZvqcpMG8x29vbsWnTpjjXJRKJoKOjg19ln0LYBCMMnu+YxM+OjWHaGwYAVBVk4GMHK3FPczHk0qsLHFFKMTQ0BJ/Ph23btsVF6CcnJ2GxWOIE2tPTg/T0dNTV1QGY3+ri05/+NDMyMnLz5Zb6HTp0iAB4HsCFw4cP//tVvREIQJxisfhgbW3tkSeffFL89re/HYQQviVFfX09gPk5Znt7O+rr6+Ms5qXmnimETSDM4PddRjx9Qs+7u/mZaXhoXxnu310GVfrVLUBYONfkGB8fh91ux86dO/k56NmzZ5GXl4eKigqwLIvf//73+Jd/+ZeQwWBouJxi+UOHDr0TwF8AUAA3HD58+OTVvId1FSchpEyn0w0+/vjj8g984AMQiUSYmJiA0+nE9u3bQQhBJBJBe3s7ampq4uaYS809UyQXEYbFn/qs+PlxPS5aPQAAhVSMd28vxn27y65q9cvCuWbs6zMzM7xwGYbB6dOn+W4XwWAQTz31FH74wx86TCZTFaV0brlzdb7nxvRfRhXWaSpiIyCzAKIAmg8fPuy/0vGvW1aeEJJRUlLS+cgjj8jf8573QCQSYXp6GkajEc3NzSCE8JHaqqqqOAFyFjM20ZwiOZGKRbi7uRivffYAnv/IbhysyUcgwuCFDgPu+MkJ3PPkKfzvWSMC4ctPQVZVVUGtVqOrq4vfE4Z7XalUoqenB5RSiMVitLS0oL+/Hz6fD3K5HA899BD+8R//sUCn07W+tXHyJeliJN+2UnH220ThzCoS/QaAGgDfvuxBx7AulpMQIiotLT317ne/e883vvEN5OTkwO/3o6Ojg+/CTinFuXPnkJubGxd9W8qSptg4DE/N4TcdBvyx24S54HzhTrZcgndvL8F9u0ovu/HYQkvJcf78eSgUCj6u4XK50Nvbi3379kEqlcJgMOCf/umfcPr06f80mUyfWer4ne+5UfxCVD5no2LFZyU+VkZw6ruRjIsA+RiA/YcPH26/go9hfcRZXFz8o127dn3+Rz/6EcrLyxGNRnHq1Cls3bqV35tkdHQUPp8PTU1N/N9Fo1G0t7ejsrIyNce8DgiEGbxy3oIXOgzoNbn417cWK3GouQh3bS2CJjtx0+qFDA8P89VEsR3mueuJu9GbTCaYzWZwLW/OnTuHT33qU+zg4OADS/UkeuqeW770Oiv/v28ThXBQPL87nJuSL/00mvE5AD4A2w4fPnzZGzFdc3FmZWW9r66u7rc/+9nPRC0tLaCU4syZMygsLOSjrVNTUxgdHcXevXvjdkXu6OiATqfjV7qnuH7ot7jxmw4DDvda4A3NW1NCgL2Vebi7qQi3bS5cds+X2A4LHOFwGG1tbXHF8hcvXgQhBA0NDaCU4vXXX8cjjzwSHR0d3bWwq1/ne24kv43KnSYqzvmsxAf5W5sGAAj+Oir/5CSVPAvge4cPH/7q5b7naypOQkhZeXn5yE9/+lPpnXfeCUIIBgcHwTAMGhsbAcwvoD179iz27t0b19d0oQuS4vokGGFwZMiOl3ssaB20I8zMzyVlYhEOVOfj5k0a3FSvhjqBReV6C+Xk5MRNlbhi+X379kEmk4FSis7OTpSUlKC4uBjRaBTPPfccHnvsMY/JZCqK7Yv7n/fcev8brOzX+0QR3CgOx56OBXDsu5EMPUA+DGDP4cOHz1zOe71m4iSEkNLS0sFPf/rTtV/84hchlUphsVhgMBiwe/duEEIS3sUALIrgpkgBzBfav3HBhsO9FrSNTYONuZSbSlS4ZZMGb68rQIM2my9yWBiZ5bDZbNDr9dizZw9EIhEikQja2tr4NigejwePPvooXnnllf+dnJx8H/d3/3T3HcZxKi75rMSPDLJYS7OUfP4/oxlfAuACsOPw4cOhlb6/aybOvLy8z+/fv/+HTz31FCkuLobP58OZM2ewf/9+SKVS3m3ltiXn4DqA79u3L9URL8WSOOZC+NvgFP46MIUTI9MIRf8enc3NkGFvZR4OVOfjQHU+CtJFOH36NHbt2hVXLD88PIxQKIQtW7YAmF+839nZiQMHDkAqlaKvrw8PPfQQ29PTcyPDMMd+fM+tNx1hZX/dIYrg1niryUEB+P87qnjYTMXPAPju4cOHv77S93RNxEkIKa+oqBh++umnpe985ztBKUVbWxsaGhqQm5sLAOjv74dUKuUrNQAsiuCmSLESAmEGJ0Yd+OvAFE6OTMPijo/FFCnlaNRmIB9zuPtAEzaV5CDrrabZZ8+eRUFBAb/M0GAwYGZmBs3NzaCU4vnnn8fXvvY1j9FoLPrkoTvGzFSs2SEKQw4gAiBMCSIAIiAIv/Wvj5I5F0QZmBdryeHDh20reR9rLk7OnX344Ydrv/CFL0AqlWJsbAzBYJCfZ1osFpjNZrS0tPBuazQaRVtbG7Zs2bJod+kUKVYKpRQTTj9Ojk7j1Og02sam4QkuXledIRNDo5RDk5UGxu+BNj8XqkwFRITAbLVBkZkJiTQNHn8Q7a+8gJu9/c6xtCx+PSIBhQyAFBTS+X+plICVgLIUxGGg4lYAZgDfXmlhwpqLMz8//4v79+//wZNPPkmKi4vh9Xpx9uxZHDhwAGKxGKFQCG1tbdi/fz/f+yVRBPd6h2VZ+Hw+BINBhEIhBINB/hGNRvntC0UiEfLz80EIgVQqhVwuX/RIT0+/bufuDEsx5vCix+hCt2EW7aNTsMxFEIquXAeUUnyw/XHMpKcPNMrI/VpCpzSEmVPMG8vwrj+2rsqa0DUVJyGkoqKiYui//uu/pDfddBPvzm7atIm3hmfOnIFOp4srKBgaGkI0GuUt6/UGwzCYm5uDy+WCy+WCx+MBy7LIzMzkBZaWlsb/XyKR8HWiAEAIAcuyiEQicSLmHn6/HxKJBEqlEkqlEiqVCllZWdelYCml6OjogEpdCFFGLqzuINyBCCZMFvjDURQUaCAVE8x5XAj6faivrkRuuhRtrX/GE99+1GsyGQvXqpv8mi16fKsK6M1PfepT0htuuAEAoNfroVKpeGGazWaIxeI4Yc7MzGB6ehr79u1bq6EJkkAggKmpKdhsNgSDQWRnZ0OlUqG0tBTZ2dlXtD5VoVAs2dAsEonA5XLB7XbzCfqsrCxotVqo1errpvM9IQQ7duzAqVOnsHu3GrWa+Qgu3VqIkydPYsvmHH6Lh87OTlRoxFCr89FYcAgT5zsyX3nllecAvHtNBkcpXZNHfn7+l+6++25qMpkopZTOzc3Ro0eP0mg0SimlNBgM0tbWVhoKhShHNBqlR44coV6vl14PuN1uOjg4SI8dO0ZPnjxJR0dH1+29syxLXS4XHRwcpMePH6cnT56kIyMj1O/3r8t4rjV2u52ePn2asizLv+bxeOjRo0cpwzCUUkoDgQBtbW2l4XCYUkppX18f3bFjByMSiW6ka6ChNXFrCSG5ZWVl1l/84hcyLjp76tQpNDY2Iicnh59TlpaWxlnN/v5+pKen8y0MNyLRaJTfkl0ikaCkpAQajWbRHpXrTTAYhM1mg8lkgkQiQVlZGTQazYbuYNjb2wuVShXXEG5kZATRaBQNDQ0AAKPRiOnpaWzbtg2UUrzwwgt49NFHZ41GYwFd5QZha/JJ63S6H77rXe+S7d+/H8C8O5ubm8u7sxaLBRKJZJE763a7F7WY2CgEAgFcuHABJ06cQDAYxI4dO7Bnzx6UlJQITpgAIJfLUV5ejgMHDqCxsRHT09M4duwYRkZGEA4nzOklPY2NjRgfH0cgEOBfq66uhtPphMvlAgCUlJQgHA5jamoKhBDcdttt2LVrV052dvZHV31Aq22KAZRUVVVFOjs7KaWUhkIh2trayruzC10DSje2OxsIBGhvby89duwYNRqNvIuUjEQiETo2NkZbW1vp4OAgjUQi6z2kVcfhcNC2trZLurd+v58eOXKEf/7yyy/TkpISHwA5XUUtrbrlLC8v/68PfvCDkubmZgDzVRfV1dV8dU9fX1/chrbAfAfusrKyuGqNZCcSiWBgYADt7e3Iy8vDwYMHUVJSktRuoUQiQWVlJd72trdBKpXixIkT0Ov1cWslk538/HxkZGTAYDDwr2VlZaG4uBjDw8MA5gNtarUaRqMRAPC2t70NN910U3pBQcGXV3Msq3qlEEIaMjMzb7733nshlUrh9/vhdDr5VSQzMzNgWRYajYb/G6fTCY/Hw3fXS3ZYlsXY2BhOnjwJhUKBG264IeHW9smMSCRCZWUlDh48iEgkguPHj8NkMvGpnGRn06ZNi9zbyspKPpIOADU1NdDr9WAYBkqlEvfffz8UCsWjhJBV23JtVcVZXl7+wgMPPCDetGkTgPl8Jdf+ntL5Dt3cxBqYz+f19fWhqalpQ1y8Ho8HJ0+eRCQSwcGDB1FeXp7UlnI5JBIJ6urqsHfvXkxPT6OjoyPugk5WJBIJNm/ezHdKAOZvSLW1tRgaGgIASKVSlJSUYHx8vsXQvn37cNddd6UVFxc/vlrjWLUrhxCyR61WN73rXe+CWCyGx+OB3+/nraTdbl+UdxsfH0dxcXHSu7Msy2JoaAg9PT1obm5GfX39ddU3Ny0tDc3NzaiqqkJ7ezsMBkPSW9H8/HykpaXBbrfzrxUWFsLj8cDrna85qKyshNFoRCQSQXp6Ou677z6kpaV9lBBSuNRxL4dVESchhJSVlf32wQcfFHE9ZjkryVnNwcFBvpseMD8nMxqNSZ828Xg8OHXqFADgwIED1/UuZgUFBThw4ABmZ2c3hBWtr6/H4OBgXOUV9xoAiMViVFZWYmRkBACwY8cOvO9975OWlZU9uRrnXxVxSiSS26uqqsruuusuflNbAPyKE5PJhLy8vLjtEUZGRlBZWZnUy8CMRiO6u7uxdetW1NXVbWgXdqVIpVI0NTWhsrIS7e3tmJ6eXu8hXTHp6enIy8uDyWTiXysoKOCrqwBAp9PBbrcjEAhAJpPh/e9/P9LT0+8ihFRf7fmv+moihJCioqL/fvDBB0lJScmiuSUXIIntYBAMBmG325O2qJ3S+a0frFYr9u/fD6VybbZdT2bUajX27t2LgYEBTExMrPdwrpja2lqMjY3FRaQbGhowMDAA4O9zUS6S29jYiPvuu09cUVHxy6s992rc6g9UVVXl33TTTSCEYGpqKm5uOTExgcLCwriWI0NDQ6itrU1KS8N1mBeLxdi5c+d1Nbe8XORyOfbt2wen04ne3t6kTLnIZDIUFhbG3WBUKhWkUikcDgeA+LmoRCLBoUOHIJPJ9l3t3POq1VFZWfnju+++W8R1wxsdHeUXTEciEUxOTvJ7nQCA1+uFx+OJ63aQLHi9Xpw6dQqlpaWLNm1NkRixWIzt27cjIyMDp0+fTsrqoqqqKkxOTiISifCv1dfXY2BgYL4GlhDU1dVhdHSU/9kdd9whKS4u/trVnPeqxEkIKUtPT996yy23gBCC2dlZyGQyPvo6Pj6OsrKyOOsSGyhKJjweD86cOYPm5uZUW87LhBCC6upqPpobCq24jY4gkEgkKC8vx9jYGP9aZmYmlEolbLb5pgYFBQVwuVwIhUKQSqU4dOgQRCLRRwghK+vdmYCrEmdpaem37rjjDjE3n9Tr9Xz0lWVZmM1mlJaW8r/vcrkQjUaRn59/Nae95rjdbpw9exY7duyASqVa7+EkLVqtFvX19Whvb+eT+clCWVkZbDZb3I2luroaer0ewPwNqKysDJOTkwDmI7cHDhxIy8zM/NCVnvOKxUkIUQB4/6FDhyCRSBAMBuHz+fidpC0WCzQaTZzVHBoaiitCSAa4tok7d+68rtMkq4VarUZjYyM6OjqSyoKKRCJUV1fzrisAZGRkQCKRwO2e3xxYp9PBbDaDZVlkZWXhXe96F8nLy7vi3cauWJxyufwf9+/fn7Z161YA84Gf8vJy3l0dHx+PW2Hi9XoRjUaTyvJwLVVaWlqQmZm53sPZMOTn56OhoQHt7e1JNQctKiqC3W5HNPr3HkSVlZW8uyuRSKBWq2G1WgEAe/bsgU6nyyWE7LiS812xOLVa7XfuuOMOZGVlgWVZWK1WFBcXA5ivoVUoFHEd88bHx5Oq4CAcDqOrqwvbt29HVtbl7c2RYnnUajXq6upw5syZpIniikQiFBcXx+U98/PzMTc3x3sB5eXlfGS3uLgYt9xyi6iiouJbV3S+K/kjQkh9Xl5eIddKxG63Iz8/ny8omJiYWLT5kNPpTJqNh7hdtevr61M5zDVEq9VCo9Ggr69vvYeyYjjxxVYNxc41MzIyQAiB1+uFSCTCbbfdBpZlbyaEXLbrdUXi1Ol0j956661ibsX45OQkv3o8HA7D4/Hw1UHAfO9PnU6XNBHaCxcuID8/P2luJslMVVUVotEoX0AudGQyGVQqFZ/jBOYtpMVi4QUbK9aGhgbccMMN0szMzPsu91yXLU5CiBjAvbfffjvEYjECgQAikQgfLDGZTHFCpJTCaDQmTTXQ5OQkQqFQak+WawQhBM3NzTCZTElT6hfrugLzJYsqlYoff2FhIex2O98x8Y477iAFBQXfuNzzXInl3N3Q0JDGFbHHCo8TYuwuYDMzM8jKykqKbm6zs7MwGAzYtm1b0lj5jQBXbdXX15cUxfIqlYpvM8oRay1FIhHUajWmpqYAANu3b4dCoVATQgoSHnAJLlucpaWlH921axfhUiY2m42v9nG5XMjIyIgr1ZucnEyKhdQMw6C3txfbtm1L6mL8ZEUul6OxsRG9vb1JsdystLSU74QAzAvW5/Px0eeioiI+altaWoqWlhZJWlraXZdzjiuxnO8+ePAgCCEIBAIQi8W8VbRYLHzEFpgPBC2cfwqVwcFB6HS6VMpkHVGr1VAoFHEXvVApLi6G2WyOCwwVFhbyFUMqlQputxuUUr7GWKfTfe5yznFZ4iSEVGs0mkzOpZ2amooLmjgcjrht1bj0itBdxJmZGbhcrqRK9WxUNm3ahLGxMcG7t1KpFNnZ2ZidneVf02q1vCtLCEFOTg6/fHLnzp2IRCINhJC0hAdMwGWJMy8v7x937dol5mpLbTYbL06v18tvDcAR+3OhwjAMzp8/v2FapSQ7Uqk0adzbWDEC843AvF4vGIbhf85Z0srKSjQ1NUkBvH2lx78scWZnZ3/i4MGDkEgkiEajCAaDfJH7QivKMAz8fr/g3cTh4eGUOysw1Go15HI5zGbzeg/lkqjV6rg2JoQQ5Ofnw+l0ApgvUOBSLkqlEnv27CGVlZWfXunxVyxOQkhOWlpaYVNTE4DFLqzNZovrqjc9Pc3vdiVUgsEgpqamNmwj62SmoaEBo6Ojgq4ekkgkkMlk8Pv/vqNfrLWUSCRQKBTwer0ghODAgQMIh8PvJCsUxYrFKZVK79yxYwdfeBAbpQ2Hw2AYJq5cLxlc2mRe9L3RSUtLW7TIWYjEihEA8vLy4HQ6eZc89ufV1dWorKyUA2haybFXfFWWlZU9sm/fPqJQKEApxezsLL+9gt1uj7OalFLMzMwIOkqbzIu+rxe4Rc6xheZCY6E4RSIRsrKy4PF4AAAajYafl6rVauzatUtUWFj4wEqOvSJxEkLE4XB4y65duwDMJ+tVKhXvsi60km63G9nZ2YK2SFw3QCG73dc73AZKsYuchYZCoUA0Go3rkhArWLlcDpZlEQ6HIRaLccMNNyAtLe3+lRx7peqpKysrk3KVP3a7nRcjy7KYm5uLW+sodJfW7XYjEonEzZlTCJPy8nJYrVZBLy1bGBiKrQ7innOBoaqqKojFYtVKUiortZw7qqurRZybGuvScv+PtUB2ux1qtXolh14XxsbGUF191Z0LU1wDRCIRSktL+dI4IbLQtZXJZBCJRPwNJScnh8+HqtVqVFRUiAFsWe64KxJnRUXFoZqaGshkMlBKEQwGIZfPt0ZxuVy8UAEgFApBLBYLcls7YH58c3NzSdcq5XpGp9PBZDIJNnKrVCrh8Xji8rIqlYrvbctVC3H/r6qqEslksp3LHXdF4mRZ9sDmzZsBgBcmZyldLldcd4OFz4UGt7wtNddMHqRSKfLz8+NcRSFBCEF6enpcVVOsOGUyGcLhMCilkEgkqK2tRWlp6d3LHXdZcb61RCyXa9S1UHxzc3NxnQKELE5KKcxmc9yqmY0ApRTRaBR+vx+BQADRaFTw1TWXS3l5uaBdW6VSyYuRe85ZS2B+EbbP5wMAbN68GZFIZFnLuZKOyLUlJSViLlUSK75oNAqRSBQXlXW5XIK9+Kenp5GTk5P0jaC5rRVdLhcf3JJKpXA6nRCLxVAqlYhGo0hLS4NKpYJKpUJeXl7caqFkIysrC9FoFIFAIC6fLhRUKhVmZmb4tqmZmZn8hkfcz10uFzIzM1FcXAyxWJxFCJFRSpeMdC17lRJCdlRVVYm5YJDb7eaXgLnd7kVtPPx+f9yeKELCbDYnzaLvhXB9miYnJ0EphVqthlarRV1d3ZJrZYPBIFwuF1wuF0ZHRyGXy1FRUSH4yq2lKCkpgcViiWtSLhRUKlVcNwdCCKRSKcLhcFz3hJKSEj4opNfrNwM4t9Qxl3Vry8vL74oNBgUCgbhgUKwLGwqFIJPJBPnFJ0NhRCI4V/zYsWNwuVxoamrC/v37UVNTg4KCgksuYpfL5Xyv2BtuuAF1dXUwmUw4ceJE0nQdiEWj0cRFRYWEXC5HMBhcMigU6/bm5ORwQaFdlzrmsuKklB5cKhi00HIKeb7JjU2IN46lCIVCOHPmDKamprB//340NjZe1V6mSqUS27ZtQ0tLC4aHh9HX1yfo6puFKBQKMAwj2JynQqGICwrFzjvT0tLigkI1NTXQ6XSHLnW8S4qTECKilOZxbuxC8Xk8nrjiAyGLU+iFEQuxWCxoa2tDaWkptm/fvqptXtLT07F3715kZWXh5MmT/JrDZECj0cQl/IVErKVM9Dw2KLRlyxZEo9GrspxFarVazOUEPR4PbykTBYPcbrdgxWm325OmImh8fBwTExPYv3//mt1QCCEoLy/Hrl270NfXJ9g0xUIWJvyFRGw+E5gPCnFiBOIt6Vs13Zdcp7icOAtzc3MJ50rFBnt8Pt+iNZBCDQZxG5sKtTAiFr1eD7vdjt27d1+TpmicFR0aGkoKgWZnZ2Nubk6QqaKF6RRCCAghfPFEbC40IyMD6enpYkLIkvOUlYhTxAWAQqEQHwyKrRICwG+FJsQ53ezsbFIEgsxmM6amptDS0nJNm4zJZDLs3r0bg4ODgndxCSHIysrC3Nzceg9lEXK5fNH+L1ygaOHPMzMzkZubSwAsuSzqkuKUSCTFOTk5fH4sFArx/w8Gg3F5My5SK0RcLpfgO7f7/X6MjIxg586d69L9Ly0tDTt37sT58+cFHyRamOAXCokMU6w409LS+P8rFArk5uaKcKXiLCwsbMzNzY1rEM39f6HlXPhcSAg5UAXMf669vb3YvHnzuhZIpKeno7y8nN9SXagsDLQICZFIxPcQAuIFudCK5uTkEABLbvZ6SXHK5fI6LhjEsmzcnSHWxQWEK05K6aKxCo3JyUlkZGQIohi/rKwMc3NzfB8cIbIw8CIkYgXIPedcWalUyq/7lEqlyM3NRUZGRtlSx1puzlnGle3FurRA8ljOQCAgyCAVB7dPyKZNm9Z7KAD+vj3ChQsXBBl0AeYvbKHWDycSJ/c81rgRQpCXlwe1Wt241LEuKc5QKFTAhfKXs5QLxSsUFuZihYbBYEBJSYmg6n3T09ORmZkp6Cqi9PT0uMZaQiEtLS0uKLRQrGKxmJ/T5+fnQyqVLrkpz3JFCGncipOFASCWZeMCF0K1nEIdFzDvchsMBkFuV1FVVSXonb8WXvRC4VKWE4gXr0ajQTQaXXKVyJLiJISIxWKxhHMJYy/yRO6EUEUg1HEB8ymerKwsQeZflUolAoGAYLeGTxZxXsqSajQaRCIR1VLHupTlzFepVCKu0CD2IueWKMUiVBEIORjkcDgEXVKo0WgEGxhKFnEuTK/E/jw9PR1isXjJO/OlxJkul8t5ETIMw8+LEomTUirIbntCvWkAwk/xCDllIVRxcl0PloLbLYH7PyFkSdFcSk0SsVjMzytjUylCFWIiFs6VhYRQyx05hCzOhe6iUBCJRJeMIhNC+J+/pa0lS+qWEyfhRBgryIU5TyEj1BuJkNe+cnA5OiGmLJazUOtFbC1tLNxnGCtertx1qe0ZlrWciaqDYv+f4srw+/1XtTbzWiGRSOIqXoRCrAUSEonGFftarHhjVnUl1OGKxZmsbq1QYVk2KT5DkUgkyJaUy7mP60WicS20lrFurVgsJliiXdClro6UaVxjksH7EKqFAhKn9IRAonHFijOWt55ftuWMMgyTUPFL+dUpVs7CAmmhwjDMuqySWQ6hem+JPKLY12I9UJZlwTAMBRBZeBzgMsS5lGlOcWXErlYQMtFoVJDiFGpQMlE8ZmG8Jlaob92gE96llxMn5SxkrCCF6u8nQqg3koXNoIRIOByGVCoVpAgS5dqFwFIWPVEwlWVZUEpBl7hALyXOCMMwiBVnIqEKHZlMJsh8GCFE8NYzUV9ioSDUhRbLWfSFbi2AJYV0KXGGwuEwX80gFov5OZJEIonbjxAQ7jxUqJUkgLCT/ICwK5iEWvm1nEWPncO/NW28InE63G43y3UPi73IEyWAhVqxIWRxCnlzHmC+Y2FeXt56DyMhQhVnot5aS/3c5/OBZdmEwSDgEuKklEbC4TDDiTNWfMv1ShESQh0XABQUFGBmZkaQUVufz8fvniVEkkWc3HYMHLHu+Ft72yzZqeySsWiWZUNcl7NE1faxbmzKcl4+hBAUFhbCZDKt91AWMT4+jrKyJTtorDtCnXNeTscQq9UKqVRqXupYlxRnWlraDNdde7lFpEIVQWZmpiDbKHJUVFRgfHxcUPP1cDgMh8PBNT4WJF6vV5Dlj8u174lGo/ycdHp6GtFodHSpY11SnCKRyMB1105mccZuxSY00tLSoNVqMTY2tt5D4env70dtba0gk/xA4t0GhEIicS5l4Z1OJ2ZnZy8udazlegiNcH1kFuY2E63wFqJbG7sVm1Cpra2FxWKBx+NZ76HAZrMhGo3y+0wKkdhtQYTGwsX9sc8jkUjcmuiZmRnMzs7qlzrWJcVps9n6nU5nnCi5/yeL5QSE24SYQyQSobm5Gb29vevq3obDYQwMDKCpqUmQhQccQm4SHo1G45q1xVrSWKGGQiHMzs5SAJaljrWc5TTOzs5SzurEplAS9UoRqjiFnk8E5m8gWq0W58+fX5cCD4Zh0NXVhfr6ekEGWmIRav51ud5asS5uIBDAzMwMC8C61PGWc9qts7OzdKmO1bFi5JYWCbFySKVSYXZ2dr2HsSzV1dUQi8XXvGcsJ8zCwkJBB4E4hNrudLneWrH/93q9mJmZobgacc7MzFCuP2hsPWjsXoMcQnVtua3YhJhPjIUQgs2bN/PbM1wLFzcajaKzsxMFBQWoqKhY8/NdLT6fD3K5XJDBoIXuNqU0riIodld4r9cLr9fLAlgy0LDcOzRPTU2xXAe27Oxsfu6WqOu2UN1HQgjy8/MF20kuFk6gGRkZOHXq1JoGiaanp3Hy5EmUlpaisrJyzc6zmgh5E+SF7vbC3QZia5Wnp6fBMEzgSsv3QCmNsCzr5pLkiXbqjU1TCDnwIuRNVxdCCEFNTQ22bt2Knp4ejIyMrKoVjUajOH/+PEZGRrB7924UFxev2rHXGpvNBm6LEKGxcKHAQkvq9Xr5PW37+/shl8vPX+p4y/oGUqm0s6+vD8DiFvgLxZqTkyNIywkAeXl5WBh5FjpKpRIHDhwAwzA4fvw4Jicnr8o1j0QiGBkZwYkTJ6BUKrFnzx4oFIpVHPHaEg6HwTCMYMe8cEPpWEsaiUQgFov5RfZDQ0OwWq0vX+p4y4pzcnLy5ZGREUSjURBC4pZgLRSnXC4X7BpFkUiErKwsQeQSLweRSIT6+nrs3bsXwWAQx48fR29vL2w227KfNaUUPp8PZrMZXV1dOHXqFMRiMQ4ePIiysjJBp0sSYbfbBWs1E619jRVnrFX1eDwYGxtjvV7v6Usdc9ndcyKRyJmxsTF2dnZWVFBQwAtSo9FAqVRiaGiI/11CCC9QId7dioqKYDKZBJsjuxRpaWmoq6tDbW0tHA4HnE4nxsfHEQwGoVAoIJfLIRaLQSkFy7Lw+XwIh8NIT0+HSqVCVVUVVCpV0gkyFpPJJJjd2BaSKBgUm9eMFardboder2cA9FzqmCvZ2urCxMQEY7fbF4lTJpMhEonEre7m9k4Uoji1Wi0GBwfR0NAgyGjfSiCEQK1WQ61WA5i/CILBIILBID8vFYvFUCgUgs9XXg6BQADRaFSQKRRg3jIuDAbFasDtdqOmZn5DMZvNhmAwGKCUXnKbtGWvUEppmGEYDxcUWhj0WVi7qlQqBTvvFIlE0Gg0sFqXTC0lHYQQKBQK5OTkIC8vD3l5eVCpVBtKmMD8BsOlpaXrPYwlWWg5F0Zu5+bm+PloX18f5HJ533LHXJH5kEqlZ/r7+wEszm8unHcKNZ3CUV5ejomJifUeRorLgGVZWK1WQUeVYyOxwNLBIJZlMTQ0BIvFcslgELBCcRoMhpeHh4fjgkJcGd9CS6lQKBAMBgWb8M/IyIBEIhFsyifFYiwWCzQajSC7AALzwkxPT18UDIoNAHH/d7vd0Ov1dLlgELBCcYbD4c6xsTGWE2FsOVxubm5ccp9L+At5V+SqqiqMji65jC6FgKCUQq/XC7p6aWFhRDQaRTgc5uecs7OzvDgdDgfGxsaiWCYYBKxQnAD69Xo9Y7HMF9AXFBTwvW8kEgnkcnncvFPoCf+8vDyEQqGU9UwCLBYLVCqVIAOMHFNTU3EpnunpaRQUFMT9nAvgGQwGBAIBP6V02UXGKxInpTQsEonGz549CwB8KRyX0NdqtXGNqnJzczEzMyPYhD8hBA0NDRgYGFjvoaS4BCzLYmRkBLW1tes9lCUJh8NgWTZuDWesJY1EInzhBKUUx48fB4CXVnLsFecTbDbbT06dOoVwOAyRSITMzEw+oa/RaOIspUgkiqvDFSI5OTkQi8WCdr+vdyYnJ6HVagXZyItjodWklGJ2dha5ubn8zzmr6XQ60dnZSQ0Gwy9XcuwVi9Pn87145swZxmg0Aoh3XRUKBRiGies2IHTXFgBvPYVq4a9notEoJiYmUFVVtd5DuSQL55tcIIgLDk1NTfE/Hxsbw8WLF8MAlg0GAZchTkqpxe12z168ON/yRKPRgGv+lei5Wq2Oey5EMjMzkZWVBW4unUI4jI6OorS0VJBbLnCwLAuv14usrCz+tVixsiwbF6k9ffo00tLSOiml0ZUc/7LKZEKh0LPHjx8HwzCQyWQQiUT8+s2FllIqlUIikQi21pZj06ZNGB4eFmT/o+sVt9sNu90u6AgtMB/4yc/Pj0uh2O12PhjkdDqRm5sLQgi8Xi/a29sxNjb2k5Ue/7LEabPZnu/s7GQ5i6jRaPhAUHZ2Nubm5uKWNiWDayuTyVBfXw9u5U2K9YVlWfT29qK5uVnwJZYLXVq/3w+ZTMZb+1iXdnJyEufOnYtQSt9Y6fEv99336vX64MjICIB48RFC+GVZHFyhudApLCwEISTl3gqAkZERaLVawdbQcjAMA6fTGbddRaxYKaW8ZQWAc+fOgWEYI6V0xU2UL0uclFIqk8n+curUKVBKkZmZyRckA/MXeewFLpfLIZPJkmKZ1pYtWzA0NJRyb9cRzp2trq5e76EsC7foO9a6W61WXpxc1ZBYLEY4HMapU6fgcDj+f5dzjsv2G/R6/VPt7e2US5MUFBTwgZ/8/HzMzMzwYgWAsrKypKhllclkaGhoQE9PTyp6uw5Eo1H09PSgqalJ8O4sMO+mxm5XMTc3B5FIxBdLWCwWXqhmsxmdnZ3M3NzcHy/nHFfyKRzt7e2NcOVvpaWlMBgMAOZd26KiIpjNf9/+QaPRwOl0CrbWNhbOnRocHFzvoVxXUEpx7tw5VFZWCt6dBcAv/IjdDsJgMPBipZTCarXyjbn7+/sxMzPjoZQaLuc8ly1OSmlIJBId+8tf/gJKKbKyssAwDN++pLS0FFwuFPj7Zj2xghUy9fX18Hg8STFX3igMDQ0hIyMDOp1uvYeyIiYmJuKsJsuysNvtvKV0OBzIycmBRCJBKBTCn/70J/h8vicu9zxX5D+Mj49//Y033mC5+WWs9VQoFJBKpXHzTG6ZVjK4i4QQbN++HWNjY4Je+rZRMJvNcLlcgu1wsJBoNAq73R7X39dqtUKtVvPueKzLOzQ0hNbW1uj09PSTl3uuK3XuO8fHx92dnZ0A5qOyVquVF195eTnGx8f5X5bL5cjIyMDMzMwVnu7aIpVK0dLSgu7ubsHnaZMZl8uF0dFR7NixI2nap5hMJhQXF8fNiycmJlBeXg5gfpuFQCAAlUoFSimOHDkChmG6KaWXXSd6ReKklFKPx/P4a6+9hmAwCLFYjPz8fD6tolarMTs7G1fOV1lZKaidtJYjIyMDW7duRWdnZyqCuwZ4PB50d3ejpaVF0FVAsVBK44QIzEeYJRIJP/+cnJzk3XOHw4HXX3+djo2NffVKznfFYTGXy/WLI0eORLiVHdw+k8C8a1haWorJyUn+93NychAOhxd1iRcyeXl5aGhoQEdHh6B3KUs2vF4vzp49i5aWFkHusbkUDocDKpUqbqfqsbExviE3y7Iwm828OM+ePYsLFy54ARy5kvNdsTgppbOU0iOvvfYaWJZFZmYmRCIRP9csLS2FyWSKqxiqqamJ69aXDKjVatTV1aG9vT1lQVcBj8eDM2fOYMeOHXE1qUKHUoqhoaG4HGwoFMLc3BxfaMDNPSUSCfx+P1588UUEg8HvU0qvqCP4VSWUxsfH//mVV15hOAtZWVkJvX5+u0GJRIKCgoK4ogS1Wo1AIJAURQmxaDQa1NfXo729PTUHvQrcbjdvMZMhZRKL1WpFdnZ2XJ+g8fFxlJeX8/Pl8fFxvh74woUL+Mtf/hKZnp7+6ZWe86rESSnts9vtxiNH5q12QUEBXC4Xb2Gqq6sxOjrKW89kXuSsVquxefNmtLe3J8WOZULDarWiu7sbO3fuTCqLCcy7q8PDw6irq+NfC4VCsNlsvAs7OzsLmUyG9PR0MAyDV155BQBepZS6rvS8V12KYTQav/DSSy/R6elpEEJQWVnJ9+eRy+XQaDR8mgUAX6WfDJsKLSQvLw+7du3C+fPnU3nQFUIpxfDwMMbHx7Fv3744y5MsGAwGaDSauEXfw8PDqK6u5qO23HMA0Ov1ePXVV5mJiYkrCgRxXLU4GYZ5pbu723vixAkAgE6ng8Ph4N2/6upqjI+Px5X0JfMi54yMDOzbtw9msxkXL15MyvdwrYhGozh79ixCoRD27NkTF0hJFqLRKMbHx+Pmmn6/H7Ozs3yrTqfTCUIIcnNzwbIs3njjDbhcrhFK6VUFWK5anJRSxuVyffHXv/41HA4HCCGoq6vjAz9SqRSlpaX8XBQAsrKykJmZGdd3KJmQSqXYtWsXCCGpVMsSeL1enD59GgUFBdiyZUtS1MsmYnx8fNGi74GBAdTX14MQAkopBgYG0NDQAGA+evvss8+y4+PjD17tuVflE/N6vf/d1dU188Ybb4BSCq1WC6/Xi7m5+dUx5eXlMJvNcekITsDJanm4+XNZWRna2tpSy83eglKKsbExnD17Fps3b44rc0s2wuEwTCbTorxmKBTiF1TbbDZkZGTwZax//OMfMT09fY5S2nm1518VcVJKGYvF8uFnn32Wms3mRYEfsViMqqoqcOtAgfkyv/z8/Lg63GREq9Vi//79sFqtOHPmzHVtRb1eL9ra2hAMBnHgwAHk5OSs95CuipGREVRVVcU1s+asJGc1h4eHUV9fD2A+QvvrX/+amZiYuH81zr9qvgbDMK+MjIyYXnvtNVBKkZeXB5Zl+cimTqfD9PR0XCqipqYGY2NjSZ/gl8lk2LFjB0pKStDW1gaz2Zy0HsGVwLIsxsbG0NXVhYaGBjQ2Ngq2O/tK8Xg8cDqdccX4TqcTIpGIv+kYjUbk5+dDoVAgEong97//PXw+3xtXO9fkoZSu2gPA/r179zKjo6OUUkpdLhc9efIkZVmWUkqpxWKh3d3dNBaz2Uy7urroRiEUCtHu7m56/Phx6nA41ns4awrLstRkMtEjR47QgYEBGo1G13tIqwLDMPT48ePU5XLxr7EsS0+cOEE9Hg+llNJoNEpbW1tpKBSilFLa0dFBq6qqIgCK6SrpaVVn6ZTSU1ar9cJLL70EhmGgVCohl8vhcDgAzLuAPp8vLk9YVFTEr3/bCMhkMjQ3N6O5uRnj4+Nob28XdP/eK8Vut+PEiRNwOp3Ys2cP6uvrk95acoyOjkKtVsftGmY2m5GVlcXnaCcmJlBcXAyZTIZAIIAXXngB4XD4OUrp6q2NXC2Vcw8Am7du3Rrt7++nlFLq9XrpsWPHeOs5NzdHjx49GneXDQaDcXehjcTMzAw9deoU7erqorOzs+s9nKuCZVlqt9vpqVOn6JkzZ6jX613vIa06brebHjt2jDIMw78WCARoa2srDYfDlFJKw+EwbW1tpZFIhFJKaWtrKy0rKwsByKGrqKVVj29TSvs9Hs/RP/zhD4hEIsjIyEBeXh5fFJ+ZmYni4uK4Gltu1+aN2AEvJycHe/fuRVlZGYaGhnDy5EmYzea4mmOhE4lEMDExgePHj8NkMqGxsTHpitZXAtf5L7ZVCqUU58+fx6ZNm/h0ysWLF1FZWQmJRIK5uTk8//zzCIVC36eUrmrp2JoknyYmJj7ym9/8JtrT0wNgvruAwWDgV6RUVVVhZmZmQ7u3sXA7r+3evRvbtm2Dy+XC0aNHceHCBbjdbkEGj1iWhdPpRE9PD7htOPbs2YNt27bFuXsbiaXcWalUym+54HA44Pf7+Y18jx49ijfffNNvs9n+fdUHtJpmOPZRXFz85Ec+8hHqdDoppZQ6nc644ND15t4uhGEYajab6ZkzZ2hrayvt6emhU1NT6xpUCYfD1Gw207Nnz9LW1lba3d1N7XY7/51tZC7HnfX7/ZRSSsfHx+nNN9/MKpXKD9M10BCha3TXJoSk6XQ689e//vW8j3zkIxCJRLhw4QIUCgW//m1sbAyhUCiuRYXFYoHFYkFLS8uajEuIcFbKZrNhenoaGRkZyMnJgVKpXLR+cDUJBoNwuVxwuVyYmZlBJBKBWq2GVquFSqVKmu4EVwvLsjh16hSampr41TKUUnR2dqKiooLfiKi3txcqlQplZWUIBoP4j//4D/zqV786ZzAYWugaCGnNxAkAhJAdmzZt6nj22WfFLS0tYBgGJ06cwM6dO5GRkQFKKU6dOoXGxsa4hHV3dzdUKpXg2/GvBZRSPqLtdrvhcrkQiUT4fV3kcjn/SEtLQ1pa2pKlcSzLIhgMIhgMIhQKIRgM8kv2AoEA0tLSoFKpoFKpkJOTI+g9MNeS8+fPQ6FQoKamhn/NaDRienoa27ZtA8Bveovdu3eDEIJXX30Vn/nMZ0KTk5PllNI12dZgTcUJACUlJU/ecsstn/rBD37A79s5MDCAffv28XtInD17FgcOHOBD8QzDoK2tDQ0NDfxC1usZSim8Xi+8Xu8isYVCoSXnrISQRWKWy+XIzs6GQqG4bizjpZiYmIDT6cT27dv5zyMYDOL06dM4cOAApFIpIpEITp48iT179kChUGBiYgIf+9jHaEdHx4c9Hs+zaza4tfCVYx8A0kpKShw///nPeX++v7+fcoUKlFI6OjpK+/r64uYAnL+/EcP1KYTB9PQ0PX78eNw8n2VZevr0aTo1NcW/1tPTQycnJyml89flY489RsvKyrrwlnFbq8eaLxWglIZMJtNtP/7xjxluZ+z6+noYjUZ+q/rKykr4/f64Olu5XI7m5mZ0dXUhEoms9TBTXGf4/X6cP38eLS0ti2pns7Oz+Xmm3W5HIBCATqcDpRR//etf8cwzz4QmJyf/gdK1dTuvyToeSulZt9v9syeffBIzMzMQi8XYunUrenp6wLIs3ytWr9fHpVdycnJQVVWFc+fOCTLdkCI5iUaj6OrqQlNTU9w822QyYW5ujl/+FQ6HceHCBTQ1NYEQgomJCTzxxBN0dnb2Y3SN5pmxXLNFdmaz+Yt/+9vfpv/3f/8XLMsiNzcXWq2WLzyQSCTYuXMnenp64orjS0pKkJWVldoiIcWqQClFd3c3Kioq+K3hgfkeumNjY/zck2VZvpBfoVAgGAzimWeewejoaJfH43n+Woz1momTUho2Go23/PjHP2ZOn57fdbuqqgoMw/DVQ+np6diyZQu6urri9lZpaGiAx+NJmi0dUgiX4eFhpKenx602CQaDi3ro9vf3o6CgAFqtFpRSvP7665w7e9dau7Mc13R5OqW0e3Z29tvf+c536PDwMAghaGpqgslkwvT0fEPs/Px86HS6uN2+OLd3dHQ0KXsPpRAGRqMRs7OzcXl1hmHQ1dWFzZs38+WIExMTiEQifGuSjo4OfOtb32Knp6fvpZRes/Yd17x3hM1m+87AwMBffvSjH8HhcEAsFqOlpQV9fX18eV95eTmkUinfKAz4e2uQvr6+VPe7FJeN2WyGwWBAS0sLnzKhlKK3txfFxcVxW8UbjUY0NzeDEIKxsTF897vfpQ6H43Gfz/fKtRzzNRcnpZQaDIZDr732mvHpp5+G3++HQqFYFJndvHkzHA5HXJ8hhUKBXbt2obe3d0Muw0qxNlitVuj1euzatQsSiYR/Xa/XQywW821IFkZwp6en8cMf/hB9fX1HLRbL16/5wNcyT3OpBwCtTqfz/epXv+LzTAaDgba3t/O1nKFQiB45coTOzMzE5ae8Xi9tbW2lbrd7yRxWihSUUmqz2eixY8f4+lgOk8lET548yefeI5EIPXbsGF8L7vf76eOPP051Op0JgJyuh0bW46T8yYEd1dXVkddff50XZH9/P7148SL/Ifp8PnrkyJG4VemUUurxeGhra+ui11Ok4LBarfTYsWM0GAzGvW6xWOjx48d5wbIsSzs7O/lCg2g0Sp955hmq0+l8AIroeuljvU7MPZRK5QO7d+9mufYlXIWGwWDgP8ylLKXX601oWVOkMJvN9Pjx44tWOCWypBcvXuQr1FiWpX/+859pTU1NFMAuup7Gaz1Pzj2Ki4t//r73vY8XZCQSoSdOnKBms5n/ADlLyfVw4eAs6/T09FLfU4rrDIPBQE+ePLnIlbXb7fTo0aNxlnRoaIh2dXXxnltvby/du3cvm5OT8xG6zrpYd2FSSgFApNPpzn3pS1/iff5wOEyPHz9OrVYr/0G63e6EFtTv99OjR49Sk8m0xNeV4nqAZVk6PDxM29ra+BYiHFNTU/To0aM0EAjwr42OjtLOzk5+3mkwGOj73/9+WlJS8isqBF2s9wD4gQDZJSUllscee4yfR4ZCIXrs2LG4ImTOgi7sxxMOh2l7ezu9ePHidbE4OEU80WiUdnV10d7e3rgF05Qmnnvq9Xra3t7O/67JZKIf//jHqU6n6wQgpkLQxHoPIG4wQF5JSYn9u9/9Lm8dg8EgPXr0aJxAl5prsixLL1y4QNvb2xe5NCk2Ln6/nx4/fpyOj48v+lmiuef4+Dhta2vjswRWq5U+/PDDVKfT9QJIowLQAhWaOOm8QNUlJSXO733ve3Rubo5SOi/QY8eOxbm43Fwz9jUOo9FIjx49mlpudh3gdDppa2vropgDy7JUr9cvmnuOjo7S9vZ2Xpg2m41+4QtfoDqdbmC9UiZLPdZ9AAkHBRSVlJS4fvCDH/ACDYVC9Pjx43FBolAoRE+ePEmHh4cXubIzMzO0tbV1wzd2vp4xGAz02LFj1Ofzxb3OMAzt6emhZ8+ejVurOTw8TDs6OnhX1maz0S9/+ctUp9ONAcigArj2Yx/rPoAlBwbodDqd63vf+x7v4obDYXrixIm4NAvDMLS3t5d2dXUtao4VCAToiRMnqF6vpyk2DizL0v7+ftrZ2bko8BMMBunJkyfpyMgIf8NmWZYODAzQrq4uXphWq5WzmGMAsqgArvmFjzVvU3I1EEKKSkpK+j/5yU/mPPzww1CpVPxaPKVSyW/DBswXKxuNRrS0tMSt0WMYBufPnwfDMNi6dWtS7hGZ4u/4fD709PQgPz8ftbW1ca1W3G43uru70dDQwLeyZBgGPT09kEql2LJlCwghMJvN+N73voeXX355yGg07qCU+tbr/VyS9b47LPfA/BzU/thjj/FpFi7w09HREXfndDgctLW1NWFRgsVioa2trXFucYrkgWVZOjY2Ro8cOcJfB7GYzWZ65MiRuDx4okCRwWCgn/jEJ6hOpzsPQEEFcI0v9Vj3AaxokPNRXMuXvvSlOJc2UeDH5/PRo0ePxv0eRygUomfOnKFdXV3XRW/cjYLX66UnT56k/f39i6YuLMvSwcFB2tbWFhf4WRgo4lzh+++/n9XpdGcgoKjsUo91H8CKBwpk63S6c/feey/t7u7m5xOJAj+RSIR2dHTQvr6+hE2azWYzbW1tpRaLZdHPUgiH5axlKBTiv+fYgODk5GRcoIhhGPq3v/2NHjx4kC0pKfktAAkVwDW93GPdB3BZgwVEJSUlP9u9ezf7+uuv88ILBAL0+PHjcYEflmXp6OgoPXr0aEI3NxgMpqyogOGs5VI32ETTFM46xk53AoEAff7552ldXV00Nzf3s1QA1/FKH4IOCC2FSqV6sKCg4Bdf/epXJe9///uRnp7OT/wlEgm2bNnCN1r2er3o6elBbm4u6urqFm1TZ7FYMDQ0hKqqKuh0ulQv13WGYRjo9XqYzWZs3bo1rs8PMN90q6+vDyzLYuvWrUhLSwMwv9nS2bNnoVKpUFdXB0IIZmZm8Ktf/Qo//elPgyaT6WaGYU6ux3u6UpJSnMB8N3mdTnfsE5/4RMbHP/5xFBQUgFKKsbEx2Gw2NDc3IzMzE8C8d6DX62EymbB169ZF26GHw2GMjIxgenoatbW10Gq1KZFeY1iWhdFohF6vh06nQ0VFxaIbqc1mw8DAAGpra1FcXMy/PjMzg/Pnz6O2thZFRUUAwHXKwx//+Mcpk8m0i1JquKZvaBVIWnECACFEW1JS0nHnnXeWfvGLX0RtbS2Av39ZOp0OlZWVvNCWs6KBQABDQ0Pwer1oaGhAXl7eNX9P1xuUzu8sNzw8DI1Gg+rqar7JFkc4HEZ/fz+fDuOsJcMwGBwchMvlQnNzM7/FR1dXF/7t3/4N3d3dpw0Gw81UqKmSZUhqcQL8hkkvb9q06Zavf/3rZO/evRCJRAm/OAC8dTWbzWhqaoJKpVp0zLm5OQwMDIBlWWzatInf3CbF6uJwODA4OIjs7GzU1dVBLpcv+p3lrGVpaSkqKipACEEkEsGf//xnfPOb32SnpqaeMJvN/0yT+AJPenECACGEaLXar+fk5HzjkUceEd977738XCXRlwjMC/D8+fNIT09HfX19wk18ZmdncfHiRcjlctTX12+4zWLXC5fLhYGBAUilUtTX1/PTj1g8Hg8GBgYgEomWtZbAfAOv5557Dk8//XRkenr6vrm5uf+9pm9qDdgQ4uQghGzT6XRv3njjjfkPP/wwduzYsawVnZqawtDQEPLz81FTU7OogohSCofDgeHhYchkMlRWViIvLy81J71MKKWw2WzQ6/WQSCSoq6tL6LX4/X4MDg4iEAigoaEhLiA0OzuL3t7euOlKJBLB0aNH8ZOf/ISeP39+yGAw3JqM88tEbChxAgAhRFZcXPxDpVL5yUceeUT83ve+d1krSimF0WjE2NgYiouL+S3FF8J1Bff5fCgvL0dxcfGieWuKeMLhMAwGA4xGI/Lz81FRUZHQUoZCIQwPD2N2dhZ1dXVQq9X89xN7c21qauL/3mw249e//jWefvrpiMfj+ZLT6fxJMruxC9lw4uQghGwvKSl5453vfGf+ww8/jO3bt8dZUafTmfAimJiYgMFgQHl5OcrKyhLufRkMBjExMQGr1Yq8vDyUlZVt2K3YrwRKKZxOJyYnJzE3NwedTofS0tJFgR5gft+SsbExWK1WVFdXo7i4OOFNs7S0dJG1/OlPf0p7e3tHDAbDLZTSyWv9PteaDStOgLeiP1IqlZ/4/Oc/L37Pe97DW1G/34+BgQEEg8FF7lMkEsHo6CimpqZQXV2NoqKihCKllMJut2NiYgLhcBjFxcUoLCy8LjehpXR+D1Gr1QqLxYLs7GyUl5cjJycn4RQgGo3CYDBgcnISZWVlKC8v5z/jS003zGYznn/+efz85z+PuN3uL8/MzPx4I1nLWDa0ODneyom+cfPNN+d94hOfwLZt2/i7uNvt5gMPDQ0NyMrK4v8uGAxidHQUDocDJSUlKCsrW3JVSzAYhMVigc1mQzQa5bdvVyqVG3Z+yrIsZmZmYLPZMD09jfT0dGi1WhQWFia0ksB8ukqv18NutyecQjidTgwMDCAjIyMuUBcIBHDq1Cn89Kc/pT09PSNvzS0nrsX7XC+uC3EC8ymX4uLiH6Wnp3/sgx/8oOTee+9FfX09P2fkLorMzEzU1dXFWT/uLm8wGJCTk4OKiopLplcikQjsdjtsNhs8Hg+/o1p+fn7Sz1Gv5L1RSjEzMwO9Xo9gMIiKiopF3ggXnQXmN67iPt9IJILu7m785je/wUsvvRT2er2POp3OH25UaxnLdSNODkJIWXl5+a+ys7Pf9sADD4jvuecefi7DRRSHh4eRn5+P6upqPoQPxLuxkUgEOp0OxcXFCYNHHAuti1wuh0qlglKphEqlglwuF6xlpZTC7/fD5XLxD4ZhVuwVhEIhGI1GmEwm3s1dWI7n8/kwNDS0KDrLMAwuXryI3//+9/jtb38bDYfDvzYajf9MKZ1Z0zctIK47cXIQQjaXl5e/oNVqNz/00EOiO++8kw9GcIEIvV4PpVKJqqqqRZYyEAjAYDDAYrEgNzcXRUVFyMvLSzg3Xfh3sRd7MBiEQqHgxbpegqWUwufzwe1282MLh8PIyMiIG1vszSoR0WgUDocDZrMZPp8POp0OOp0uzs3lAkZ6vR7hcBg1NTV8YI5SitHRUbz44ot47rnnGK/X++bk5OTHKaWmtf4MhMZ1K04OQsj+srKy39bU1JQ8+OCD5NZbb+V3nOJynHq9HgzDoLKyclHdLaUU09PTsFqtcDqdyMrKglarhVqtXlHXBUopgsHgIsGKRCKkpaVBLpfzj9jnaWlpy94IOBiGQSgUQjAY5B+xz0OhEID5/VE5ESqVymWFyBEIBDA1NQWbzYZgMAi1Wo3CwkKoVKq4z4phGJjNZoyPjyMzMxNVVVV8rpO7Ib766qt45plnWIfD0T0xMXE/pfS63TX5uhcnMF9hJBaL/6GoqOiZHTt25D7wwAM4ePAg8vPz+d/xer0YHx/H9PQ0CgsLUVpaivT09LjjUErh8Xhgs9lgt9shEomg1Wqh0WgS5vYuBcuyywpqpd9dIqHHij0tLe2yLDWlFG63G1NTU5iamoJYLIZWq4VWq01YReXxeDA5Ocl/duXl5XypHqUUFosFf/3rX/Hss89SvV4/OTk5+X5KaceKB7RBSYkzBkKIODMz8yGVSvWj5ubmzHvuuYe84x3viAvzMwwDi8UCg8EAsViM0tJSqNXqhPPOYDAIm83GW5S8vDzeMmVmZgp2rrkQlmXh8XjgdrsxOzuL2dlZ3kPQaDQJI7PhcBhWqxUGgwEymQxlZWVQq9X858ilq9544w28/PLL7Pj4+LTBYHiQUvrG9RDsWQkpcSZg3pCK/6G0tPSJ/Px83aFDh8R33nnnohrcubk5GI1G2O12yOVy3nokKuCORqOYnZ3lXVefzweJRMLP55RKJbKystZdsJwQXS4X3G433G43WJZFdnY2lEolcnJyoFKpErrUXq8XNpsNU1NTYFkWGo0GOp1u0Wd2/vx5HD58GK+++ioTDAZ79Hr95wGcSokynpQ4l4EQsqWiouL7hJAbb7zxRukdd9yBXbt2oaioKE5IXq+Xn3exLMtHNLOzs5cUXCQS4cXqdrvh9XohFouRmZm5aI7JPb/aVEw0Gl3kLnMPr9cLAMjKyuItfHZ29pLRaC5FYrPZ4HA4oFAoeGsae4NiGAaTk5Noa2vDa6+9hra2thCl9DdGo/GbG6UOdi1IiXOFEEJyVCrVR7Kzs79aXl6uuvXWW0W33HILqqurF6UUwuEwnwucm5tDXl4eNBrNiqKdkUgEPp9v0RyTe86yLABAJpNBIpFAJBKBEAJCCEQiUVybC5ZlQSlFJBLhdwwXi8VLzj8zMjIumRbiglczMzOYmpqCy+WKy3PG/i0XKBsaGsJrr72Gv/zlL8zMzIzVZrN9LRAI/C+lNHA138f1QEqclwmZV+HOysrK70aj0bc1NzdLd+/eTQ4cOMCnBGKtG8uycDqdsNvtV5yeWAilFOFwGNFodJEQAcQJlhACqVQKqVR62UEfLu3DpVe4tI9KpYJarV5UmheJRGCxWDA4OIgTJ06gs7OTDg4OhgghvzcYDP9OKR26rDd6nZMS51VACEkD8LbKyspPhcPhWyorK+U7d+4U3XDDDWhsbERJScki4SVK7IfDYT6NoVQqoVAoIJfLIZPJ1nwOGhsVjhVjrBAvVTDh9/sxOTmJ3t5eTpCM3W73Ukr/YDQafwmgk1LKrOmb2KCkxLlKvGVRtxYWFt4vk8keUiqVOTt37hTv378fO3bsQGFhIXJychK6jbGC9Xg8CAQCCIVCfP6REJIw1ykWi3l3lrOS3PE4S8q5tYlcZO7Y3DG5YgjuBpGIcDiMmZkZmEwmdHR04PTp07Srq4uJRCIWl8v185mZmd9RSsfW5lO+vkiJc40ghBRlZGS8S6vVfpZl2YrKykpJVVWVqKamBlu2bEFxcTE0Gs2Sgo0lUc4zGAyCZdk4EbIsG+fOcqKVSCQJg0vLWWVOiFNTU5iYmMCFCxcwMjKCsbExZnJyMiqTyc6Pj48/wTDM65TS2dX8/FKkxHlNIITIADTKZLJdpaWld0ej0V0ikSirvLxcUllZyQtWp9MhLy8P6enp6+LW+v1+OBwOTExMoL+/H6OjoxgdHWVMJhNDCHESQk5MTEy8Qik9C2A45a6uLSlxrhOEECmARqlUuvMtwe4WiUTZarVanJeXR3JyckhOTg7Jzc1FQUEBn5rJyspCeno6MjMzIRaLIRKJ+MdCt5Z7RCIReL1e+P1+eDweWK1WTE1Nwel0YmZmhissYJ1OJ3U4HIxIJHICOL5AiOx6fl7XIylxCoi3BFsEoBBAoVwu12k0msa0tLRalmVLw+FwHiFEJpVKJUqlkshkMiIWiyEWi4lEIoFYLAalFAzDgGEY+ta/CAaDcLvdLMMwUQBBmUzmADDh9/uHrFbrBZZlLQCs3CMlRGGQEmcS8paI8wGkAZACkLz1kAJgAUQXPAIAplNuaHKREmeKFAJlZWuOUqRIcc1JiTNFCoGSEmeKFAIlJc4UKQRKSpwpUgiU/w+Xs8ucDGSg7QAAAABJRU5ErkJggg==\n",
"text/plain": [
"