Skip to content

Instantly share code, notes, and snippets.

Created March 28, 2017 13:38
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save anonymous/3a2afd006be4d1918f28d98ef2a1e077 to your computer and use it in GitHub Desktop.
Save anonymous/3a2afd006be4d1918f28d98ef2a1e077 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"import yt\n",
"\n",
"ds = yt.load('IsolatedGalaxy/galaxy0030/galaxy0030')\n",
"\n",
"profile = yt.create_profile(ds.all_data(), ['density', 'temperature'], 'cell_mass', n_bins=128, weight_field=None)"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[[ 0. 0. 0. ..., 0. 0. 0.]\n",
" [ 0. 0. 0. ..., 0. 0. 0.]\n",
" [ 0. 0. 0. ..., 0. 0. 0.]\n",
" ..., \n",
" [ 0. 0. 0. ..., 0. 0. 0.]\n",
" [ 0. 0. 0. ..., 0. 0. 0.]\n",
" [ 0. 0. 0. ..., 0. 0. 0.]] g\n",
"(128, 128)\n",
"1.3172140179324272e+34 g\n",
"1.388760359778708e+43 g\n"
]
}
],
"source": [
"# the 2D histogram in density-temperature space\n",
"histogram = profile.field_data['gas', 'cell_mass']\n",
"\n",
"print(histogram)\n",
"\n",
"print(histogram.shape)\n",
"\n",
"print(histogram[histogram > 0].min())\n",
"\n",
"print(histogram[histogram > 0].max())"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"1.001749126507666e+39 g\n",
"1.388760359778708e+43 g\n"
]
}
],
"source": [
"histogram[histogram < 1e39] = 0\n",
"\n",
"print(histogram[histogram > 0].min())\n",
"\n",
"print(histogram[histogram > 0].max())"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"1.001749126507666e+39 g\n",
"1.388760359778708e+43 g\n"
]
}
],
"source": [
"print(profile['cell_mass'][profile['cell_mass'] > 0].min())\n",
"print(profile['cell_mass'][profile['cell_mass'] > 0].max())"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"phase = yt.PhasePlot.from_profile(profile)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"<img style=\"max-width:100%%;max-height:100%%;\" src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAABDMAAAOYCAYAAADIS7gYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAAEt0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMC4wLnBvc3QzNjg3K2dmMjg3YWM2LCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvyjh8cAAAIABJREFUeJzs3XuUU/W9//9XJjOTASSpc0AUGy62RyQoIiWKXDTgBamWw1REbG3FCyAdz4+DipRVvF/SWkWrDUcRy1DPKeUgoFSwKBRdIkWGL0wBsQqCmGFGLooJzHCbZP/+wEkNk5nMQGZn7+H5WGuvlez92fvzDjvDWvudz+fzdhiGYQgAAAAAAMAmcrIdAAAAAAAAQFOQzAAAAAAAALZCMgMAAAAAANgKyQwAAAAAAGArJDMAAAAAAICtkMwAAAAAAAC2QjIDAAAAAADYCskMAAAAAABgKyQzAAAAAACArZDMAAAAAAAAtkIyAwAAAAAA2ArJDAAAAAAAYCskMwAAAAAAgK2QzAAAAAAAALZCMgMAAAAAANgKyQwAAAAAAGArJDMAAAAAAICtkMwAAAAAAAC2QjIDAAAAAADYCskMAAAAAABgKyQzAAAAAACArZDMAAAAAAAAtkIyAwAAAAAA2ArJDAAAAAAAYCskMwAAAAAAgK2QzAAAAAAAALZCMgMAAAAAANgKyQwAAAAAAGArudkOAE0Xj8dVUVGhtm3byuFwZDscAAAAAC2YYRjas2eP2rVrp5wcc38Pz8/PV0FBgal9wh5IZthQRUWFvF5vtsMAAAAAcIo4s0O+vth1xPx+zzxT27dvJ6GBOkhm2FDbtm0lSeFwWG63W9FoVF6vN/HebH6/X6Wlpbbt90Suk+4c7on512nMOdm8L9m6J5nqm3tizb6beh2r3xOJ/79S4Z6Yfx2r/63w/1dq/K0033V27twpn8+nzz8aKHdb8x4ho/tr1Kn7ezpy5AjJDNRBMsOGaqeWuN3upP+oj39vFqfTaet+T+Q6jT2He2LedZpyTjbuS7buSab65p5Ys++mXsfq90Ti/6+GcE/Mu47V/1b4/6th/K1k/jrRaFSS5G6bK7ebR0hYAwuA4qQVFxfbut8TuU62PnNjcU+sJ5vxZaJv7ok1+27qdax+TyT+/7Ii7on18P+XNZ2KfytANjkMwzCyHQSaJhqNyuPxKBKJJKaZfPs9so97Yk3cF+vhnlgP98R6uCfWxH2xHu5J8ykvL5fX69VX4QGmjsyIRmtU6F3JPUVKjMxoAVwulx588EG5XK5sh4JvcE+siftiPdwT6+GeWA/3xJq4L9bDPQFOLYzMsCGyzgAAAADMwsgMWBGrtwAAAAAA0oqpRjGT+wPqwzQTAAAAAABgKyQzbMzv98vn8ykUCmU7FAAAAAAtTCgUks/nUyAQyHYoQB2smWFDrJkBAAAAwCy1a2bsCvcxfc2MDt61PPcgJUZmAAAAAAAAWyGZAQAAAAAAbIVqJgAAAACAtGKKmVzNxMzeYDeMzAAAAAAAALbCyAwAAAAAAE7CwoUL9de//lUvvviiJOnw4cMqKSnRl19+qQMHDmjTpk36xS9+oWuuuSbLkbYcJDMAAAAAAGkxzSS1vXv36s4779TQoUMT+yZNmqR169bpb3/7m/Lz87Vw4UJdd911WrFihQYOHJjFaFsOppkAAAAAAHCCnn32WXXp0iVpX35+vnbu3KmamhpJ0nnnnadYLKZVq1ZlIcKWiWQGAAAAAAAnYN68eRo6dKhatWqVtP+pp57S9u3b1bp1a0nSli1bJEl+v9/0GFsqppkAAAAAANKKKa6YHKb2Z2W7du3S1q1bdcMNNzTYLhaL6emnn9bdd9+twYMHmxRdy8fIDAAAAABAi1VcXKzHH3+8wTYfffSRRowYoQsvvFCDBg2S3+/Xiy++KMMw6j3nmWee0cSJE+s9Ho/H9dhjj2nEiBHq0qVL2hjQNCQzAAAAAAAtSk1NjZYtW6ZAIKDp06fr6NGj9bb9+OOPdemll6qwsFDr16/XihUrNHv2bP3yl7/UuHHjUp4zZ84cDR8+XAUFBfVeNycnR1OnTtXChQvVq1cvXXDBBdq+fftJfzYcQzIDAAAAAJBWXPFvppqYs8VPcJrJXXfdpUAgoDlz5ui0005rsK1hGLrjjjsUi8U0bdo05eQce0T2+XyaOHGiXnrpJS1evDjpnIqKCu3YsUN9+/ZtdExjx47Vjh07NHbs2KZ/IKREMgMAAAAA0GL8/ve/18qVK/Xyyy+rT58+DbZdvXq1Vq5cqaKiojqJj5tvvlmSNG3atKT9ixcv1oYNGzR69OjEtmHDBq1cuVKjR4/Wu+++q44dO+qRRx5JnNOmTRu1a9dO77//foY+JVgAFAAAAABwSlqwYIGk1FVGzjnnHBUWFuqdd95RNBqV2+2WJI0ZM0ZjxoxJatulSxcNGDBAJSUl+sc//qHKykp99dVXieOxWEz79u3T9773vWb8NKcWRmYAAAAAANKqkWH61tzKysokSZ07d055vFOnTorH41q/fn2D14nFYorFYpKkCy64QFdffbWKi4sTxxctWqTDhw/rsccey1DkYGQGAAAAAOCUtHXrVklSYWFhyuO1+ysrK1MeLykpUUlJicrLyzV//nzt27dPDz74oP785z/riSee0KFDh3T06FF99tlnevPNNzVkyJDm+SCnIJIZAAAAAADLikajSe9dLpdcLldGrl07FaS+qiS1+/ft25fyeO2aGan89re/PfkAUS+mmQAAAAAALMvr9crj8SS2YDCYsWtXVVVJknJzU//O73Q6JUlff/11xvpEZjAyAwAAAACQVkyGYiasY/Ht/iQpHA4nFt+UlLFRGZKUk5OTWOsilerq6kQ7WAvJDAAAAACAZbnd7qRkRia1bds2qerI8WqTGceXbUX2kV4CAAAAAJySOnXqJEk6cuRIyuO101C++93vmhYTGodkBgAAAAAgrVgWtubWo0cPSdLu3btTHt+1a5ckqVu3biZEg6YgmQEAAAAAOCVddtllkqRNmzbVObZ7927t2rVLHTt2JJlhQSQzAAAAAACnpKKiIuXm5uqNN96oc2zhwoWSpBtvvFEOh8Ps0JAGyQwAAAAAQFo1cpi+SVI0Gk3aDh8+3Kh4DcNQeXm5JOnzzz+XYdStxNK+fXuNHz9e77//vt55553E/lgspunTp6tDhw6aOnXqyf/jIeNIZgAAAAAALMvr9crj8SS2YDDYYPtgMKiLL75YZ5xxhl5++WVJ0qxZs9S+fXv5/f465//mN7/R9ddfr6KiIv3+97/XX/7yFw0bNkyVlZWaN2+eCgsLm+2z4cQ5jFTpKZiqb9+++uCDD+rs9/v9WrNmTZ390WhUHo9HkUik2UoUAQAAAIAklZeXy+v1qizcXW3dTtP63R+NqZf3I4XD4aTnHpfLJZfLlfH+ysrKtGrVKkUiEXXt2lXXXnut2rZtm/F+kBm52Q7gVLdhwwa1a9dOr776atIf5HPPPZdYjAYAAAAAsi0ucyqMfLs/SXK73ab8iNurVy/16tWr2ftBZpDMyLLPPvtM8+fPT0pkGIah8ePHa9q0aVmMDAAAAAAAa2LNjBOwYsUKLVq0KCPXGjZsWJ0hUmvWrFFBQYHOP//8jPQBAAAAAEBLQjKjCdatW6chQ4Zo8ODBWrduXYNtP/30U/385z/XRRddpH79+ql3796aMWNGo/qZP3++ioqKMhEyAAAAAGREtqqZAKkwzaQRvv76a82YMUO7d+/W6tWr07YvKyvT5ZdfrmHDhmnNmjXKy8vT+++/ryFDhqi0tFQvvfRSg+cvWLBAr7zySqbCBwAAAACgRWFkRiN4PB7dd999euqppzR06NAG20YiEQ0bNky5ubmaPn268vLyJEn9+/fXpEmTNHPmTJWUlNR7fllZmaqrq9W3b99MfgQAAAAAAFoMkhmN4HD8a3hTQUFBg21feOEFhcNhjRgxok4Zn9tuu02SNHXqVB09ejTl+fPnz9d//Md/JPUJAAAAANlWoxzTN6A+fDsybObMmZKkQCBQ55jX69U555yjnTt3atmyZSnPZ70MAAAAAPiXaDSatB0+fDjbIcECSGZk0J49e7R161ZJks/nS9mmtkLJW2+9VefYRx99pIqKCg0aNKj5ggQAAAAAG/F6vfJ4PIktGAxmOyRYAAuAZtCHH36YeH322WenbNOxY0dJ0qZNm+ocmz9/vq677rrEOhsAAAAAYBUx5Shm4u/hMRmSpHA4LLfbndjvcrlMiwHWRTIjg7766qvE6+PXy6hV+0e4e/fuOsfmz5+vX/3qV43uLxqNJr13uVz8YQMAAAA4KYcPH06ayrF///4sRnPsGerbyQxAYppJRlVXVyde5+fnp2xTm2yoqqqqc2z9+vUaMWJEo/tjuBUAAACATAsGg0nPGfVNoQeyiZEZGdSqVavE66NHj6ZMaNRWMWnduvVJ98dwKwAAAACZNmXKFN19992J9zt37pTP51ONHKqReVUXzewL9kMyI4POPPPMxOsDBw6osLCwTpsDBw5Ikjp06HDS/THcCgAAAECmHT99/fjp7YAVMM0kg7p37y6H41j2sLKyMmWb2v09evQwLS4AAAAAAFoSkhkZVFhYqJ49e0qSNm/enLJNbcUTyq8CAAAAAHBiSGZk2KhRoyRJ7733Xp1ju3bt0ieffKLCwkJdffXVZocGAAAAACcsphzVmLiZWQYW9sO3o4kOHjwoSYrH4ymPjx07Vu3bt9fcuXOTqptI0qxZsxSPx3XPPfckLRZ6ovx+v3w+n0Kh0ElfCwAAAAC+LRQKyefzKRAIZDsUoA6HYRhGtoOwiy+//FI9e/ZURUWFBg8erKVLlyo3t+4aqkuXLlVRUZFGjRqlF198UXl5eSotLdWVV16pAQMGaNGiRXI6nSccRzQalcfjUSQSYQFQAAAAAM2qvLxcXq9Xy8IXqY37xJ9jmqoqGtOV3vUpqzhSyRGMzGiEWCymPn36qHPnzqqoqJAk/e1vf1PHjh1VVFRUp/2QIUO0evVqVVVVqW/fvho4cKDGjRunRx555KQTGQAAAACQDcemfzhN3I49rnq9Xnk8nsQWDAaz/C8BK6A0ayM4nU6tXbu2Sef07NlTc+fObaaIAAAAAODUkGpkBkAyAwAAAABgWW63m+n1qINkBgAAAAAgrZjJFUaoZoKG8O0AAAAAAAC2QjLDxijNCgAAAKC5UJoVVkZpVhuiNCsAAAAAs9SWZn09fInauM1bqaAqWqP/8H7Acw9SYmQGAAAAAACwFZIZAAAAAADAVqhmAgAAAABIq0ZO1chpYn+siID6MTIDAAAAAADYCskMAAAAAABgK0wzAQAAAACkFVOOakz8PTz2TV/RaDRpv8vlksvlMi0OWBMjM2zM7/fL5/MpFAplOxQAAAAALUwoFJLP51MgEMhqHF6vVx6PJ7EFg8GsxgNrcBiGwaoqNhONRuXxeKi3DAAAAKDZlZeXy+v16v/CA9Tabd7g/upojUZ6VyocDic99zAyAxLTTAAAAAAAjRAznIoZ5lUziX3zu7vb7eZHXNTBNBMAAAAAAGArJDMAAAAAAICtMM0EAAAAAJBWjcnVTMzsC/bDtwMAAAAAANgKyQwAAAAAAGArJDMAAAAAAICtkMywMb/fL5/Pp1AolO1QAAAAALQwoVBIPp9PgUBAklQjp+kbUB+HYXxTvBe2EY1G5fF4FIlEqLcMAAAAoFmVl5fL6/Vq1udXqrXbvBoS1dEa3dppGc89SImRGQAAAAAAwFYozQoAAAAASOtYaVbzpn7UKG5aX7AfRmYAAAAAAABbIZkBAAAAAABshWkmAAAAAIC0YspRjYm/h8e+6SsajSbtd7lccrlcpsUBa2JkBgAAAADAsrxerzweT2ILBoPZDgkWwMgMAAAAAIBlhcPhpNKsjMqARDIDAAAAANAIMcOpmGFeNZOYcayaidvtTkpmABLTTGzN7/fL5/MpFAplOxQAAAAALUwoFJLP51MgEMh2KEAdDsMwjGwHgaaJRqPyeDyKRCJkKAEAAAA0q/Lycnm9XoV2XKtW7jzT+j0YParizot57kFKTDMBAAAAAKRVY3I1EzP7gv3w7QAAAAAAALZCMgMAAAAAANgK00wAAAAAAGnVyKkamVfNpEZx0/qC/TAyAwAAAAAA2ArJDAAAAAAAYCtMMwEAAAAApFVjOFVjmDjNxGCaCerHyAwAAAAAAGArJDMAAAAAAICtMM0EAAAAAJBWXDmKmfh7eJzf3tEAvh0AAAAAAMBWSGbYmN/vl8/nUygUynYoAAAAAFqYUCgkn8+nQCCQ7VCAOhyGYRjZDgJNE41G5fF4FIlE5Ha7sx0OAAAAgBasvLxcXq9Xwc9uUIE7z7R+D0WPakqXeQqHw0nPPS6XSy6Xy7Q4YE2MzAAAAAAApFWjHNXIaeJ27HHV6/XK4/EktmAwmOV/CVgBC4ACAAAAACwr1cgMgGQGAAAAAMCy3G430+tRB8kMAAAAAEBaNUaOagzzViowsy/YD98OAAAAAABgKyQzAAAAAACArTDNBAAAAACQVm2VETP7A+rDyAwAAAAAAGArJDMAAAAAAICtMM0EAAAAAJBWTDmKmVhhJMZv72gA3w4AAAAAAGArJDMAAAAAAICtMM0EAAAAAJBWzORqJjGqmaABjMwAAAAAAAC2QjIDAAAAAADYCskMG/P7/fL5fAqFQtkOBQAAAEALEwqF5PP5FAgEJEkxI8f0DaiPwzAMI9tBoGmi0ag8Ho8ikYjcbne2wwEAAADQgpWXl8vr9erubbfL1TbftH4P7z+iaee8zHMPUiLVBQAAAAAAbIVqJgAAAACAtGrklNPECiNmVk6B/TAyAwAAAAAA2ArJDAAAAAAAYCtMMwEAAAAApBUzHKZWGIkZDtP6gv0wMgMAAAAAANgKyQwAAAAAAGArTDMBAAAAAFhWNBpNeu9yueRyubIUDayCkRkAAAAAgLRq5DR9kySv1yuPx5PYgsFglv8lYAWMzAAAAAAAWFY4HJbb7U68Z1QGJJIZAAAAAAALc7vdSckMQCKZAQAAAABohLiRY2pp1riJfcF++HYAAAAAAABbIZkBAAAAAABshWkmAAAAAIC0aowc5Zg49aOGaSZoAN8OAAAAAABgKyQzAAAAAACArTDNBAAAAACQVkxOxeQ0tT+gPozMAAAAAAAAtkIyAwAAAAAA2ArTTAAAAAAAacWMHMVMrDBiZl+wH74dNub3++Xz+RQKhbIdCgAAAIAWJhQKyefzKRAIZDsUoA6HYRhGtoNA00SjUXk8HkUiEbnd7myHAwAAAKAFKy8vl9fr1U8/uUf5bV2m9Xtk/2H977lP89yDlJhmAgAAAABIq8bIUY6JUz9qmGaCBvDtAAAAAAAAtkIyAwAAAAAA2ArTTAAAAAAAacWVo5iJv4fH+e0dDeDbAQAAAAAAbIVkBgAAAAAAsBWmmQAAAAAA0ooZOcoxnKb2B9SHbwcAAAAAALAVkhkAAAAAAMBWSGYAAAAAAABbYc0MAAAAAEBax9bMMO/3cNbMQEP4dgAAAAAAAFshmQEAAAAAAGyFaSYAAAAAgLRiylGOib+Hx/jtHQ3g2wEAAAAAAGyFZAYAAAAAALAVppkAAAAAANKimgmshG8HAAAAAACwFUZmAAAAAAAsKxqNJr13uVxyuVxZigZWwcgMAAAAAEBaMSPH9E2SvF6vPB5PYgsGg1n+l4AVMDIDAAAAAGBZ4XBYbrc78Z5RGZBIZgAAAAAALMztdiclMwCJZAYAAAAAoBHi35r6YVZ/QH34dgAAAAAAAFthZIYFHThwQGvXrtXRo0d11VVXZTscAAAAAAAshWSGhWzfvl1TpkxRPB7XhAkTdNlll2U7JAAAAACQdKyaicPEqR9mTmmB/fDtsIiVK1eqT58+uuyyy/R///d/6t+/v3JyuD0AAAAAAByPp2UL2LZtm370ox/ptttu0y9+8YtshwMAAAAAgKUxzcQCJk+erCNHjuhXv/pVtkMBAAAAgJTiylHcxN/DzewL9sO34wSsWLFCixYtysi1vv76ay1cuFBdu3bV5MmT1adPH7Vr10533nmnDh48mJE+AAAAAABoSUhmNMG6des0ZMgQDR48WOvWrWuw7aeffqqf//znuuiii9SvXz/17t1bM2bMqNNuzZo1isViOv/88/X8889r7dq1WrhwoUpKSjRp0qTm+igAAAAAANgW00wa4euvv9aMGTO0e/durV69Om37srIyXX755Ro2bJjWrFmjvLw8vf/++xoyZIhKS0v10ksvJdru3btXkjRmzBjl5+dLkgYOHKjrr79eL774op5++mm5XK7m+WAAAAAA0EhUM4GV8O1oBI/Ho/vuu09PPfWUhg4d2mDbSCSiYcOGKTc3V9OnT1deXp4kqX///po0aZJmzpypkpKSRPvCwkJJktPpTLpOr169VFNTo88++yyjnwUAAAAAALsjmdEIDocj8bqgoKDBti+88ILC4bBGjBihtm3bJh277bbbJElTp07V0aNHJUn//u//LknavXt3Ulu32y1Jda4BAAAAAMCpjmRGhs2cOVOSFAgE6hzzer0655xztHPnTi1btkyS9L3vfU8/+MEP9O677ya1raio0He/+1117Nix2WMGAAAAAMBOSGZk0J49e7R161ZJks/nS9nm/PPPlyS99dZbiX1Tp07V3LlztWfPHknSkSNH9Oqrr+rRRx9t5ogBAAAAoHHihkNxI8fEzZE+KJyyWAA0gz788MPE67PPPjtlm9qRFps2bUrsGz58uKqqqnTLLbfooosuUnl5uSZMmKDRo0c3a7wAAAAAANgRyYwM+uqrrxKv61vronYtjOPXyPjpT3+qn/70p03qLxqNJr13uVxUPgEAAABwUg4fPqzDhw8n3u/fvz+L0QCpMc0kg6qrqxOva8usHq822VBVVXXS/Xm9Xnk8nsQWDAZP+poAAAAATm3BYDDpOaN2Cn3McJi+AfVhZEYGtWrVKvH66NGjKRMatVVMWrdufdL9hcPhxEgPSYzKAAAAAHDSpkyZorvvvjvxfufOnfWuCQhkC8mMDDrzzDMTrw8cOKDCwsI6bQ4cOCBJ6tChw0n353a7k5IZAAAAAHCyjp++fvz0dsAKSGZkUPfu3eVwOGQYhiorK1MmMyorKyVJPXr0MDs8AAAAADhhtVVGzOwPqA/fjgwqLCxUz549JUmbN29O2aa24smgQYNMiwsAAAAAgJaEZEaGjRo1SpL03nvv1Tm2a9cuffLJJyosLNTVV19tdmgAAAAAALQIJDOa6ODBg5KkeDye8vjYsWPVvn17zZ07N6m6iSTNmjVL8Xhc99xzT9JioSfK7/fL5/MpFAqd9LUAAAAA4NtCoZB8Pp8CgYCkf00zMXMD6uMwDMPIdhB28eWXX6pnz56qqKjQ4MGDtXTpUuXm1l12ZOnSpSoqKtKoUaP04osvKi8vT6Wlpbryyis1YMAALVq0SE6n84TjiEaj8ng8ikQiLAAKAAAAoFmVl5fL6/XqB+ueUW7bk/9RtrFq9h/U/+s9kecepESqqxFisZj69Omjzp07q6KiQpL0t7/9TR07dlRRUVGd9kOGDNHq1atVVVWlvn37auDAgRo3bpweeeSRk05kAAAAAABwqqOaSSM4nU6tXbu2Sef07NlTc+fObaaIAAAAAMBchsyd+mHw2zsawLcDAAAAAADYCskMAAAAAABgK0wzAQAAAACkFTccchgOU/sD6sPIDBujNCsAAACA5nJ8aVbASijNakOUZgUAAABgltrSrL3+33NynmZeadbYgYMq+8H/x3MPUmKaCQAAAAAgrbiRI4eJ1UzMrJwC++HbAQAAAAAAbIVkBgAAAAAAsBWmmQAAAAAA0jKMHBkmTv0wsy/YD98OAAAAAABgKyQzAAAAAACArZDMsDG/3y+fz6dQKJTtUAAAAAC0MKFQSD6fT4FAINuhAHU4DMMwsh0EmiYajcrj8VBvGQAAAECzKy8vl9frVY/S/5bztFam9Rs7cFAf+sfz3IOUGJkBAAAAAABshWQGAAAAAACwFUqzAgAAAADSojQrrIRvBwAAAAAAJ2HhwoUaN25c4v2hQ4f0xBNP6L/+67903XXX6YorrtCqVauyGGHLQzIDAAAAAIATtHfvXt155506fPhwYt/jjz+um266Sc8++6zeeOMNXX/99Ro0aJBWr16dxUhbFpIZAAAAAIC0DMNh+mYHzz77rLp06ZJ4f/jwYT377LOaPXt2Yt/48eNVWFioJ554IgsRtkwkM2zM7/fL5/MpFAplOxQAAAAALUwoFJLP51MgEMh2KJY1b948DR06VK1a/atkbU1NjQoLC1VdXZ3Y53A41LlzZ23ZsiUbYbZILABqY6WlpdRbBgAAANAsiouLVVxcrPLycnm93myHYzm7du3S1q1bdcMNNyTtb9OmjXbs2FGn/WeffaY+ffqYFV6Lx8gMAAAAAEBax6Z+5Ji4ZWaaSXFxsR5//PEG23z00UcaMWKELrzwQg0aNEh+v18vvviiDMOo95xnnnlGEydObFQMS5cu1d69ezV16tQmxY76kcwAAAAAALQoNTU1WrZsmQKBgKZPn66jR4/W2/bjjz/WpZdeqsLCQq1fv14rVqzQ7Nmz9ctf/jKpQsm3zZkzR8OHD1dBQUHaWA4dOqR7771XTzzxhPr27XvCnwnJSGYAAAAAAFqMu+66S4FAQHPmzNFpp53WYFvDMHTHHXcoFotp2rRpysk59ojs8/k0ceJEvfTSS1q8eHHSORUVFdqxY0ejExPjx4/X6NGjdd99953YB0JKJDMAAAAAAGnZpZrJ73//e61cuVIvv/xy2jUqVq9erZUrV6qoqKhO4uPmm2+WJE2bNi1p/+LFi7VhwwaNHj06sW3YsEErV67U6NGj9cEHHyTaPvbYYxo0aJDuueceSdK2bdtO6DOhLhYABQAAAACckhYsWCDpWKXI451zzjkqLCzUO++8o2g0mii+MGbMGI0ZMyapbZcuXTRgwACVlJQk9v3xj3/Uueeeq5EjRyb2/eEPf9Bjjz3WDJ/k1EMyAwAAAABwSiorK5Mkde7cOeXxTp066auvvtJ6b10DAAAgAElEQVT69et1+eWX13udWCymWCyWeL9kyRL98Y9/1JVXXqlf//rXiTZff/11BqM/tZHMAAAAAACkdTJTP060v+a2detWSVJhYWHK47X7KysrUx4vKSlRSUmJysvLNX/+fO3bt0/333+/Ro4cqaqqKi1fvjyp/W9+85sMRn9qI5kBAAAAALCsaDSa9N7lcsnlcmXk2l999ZUk1VuVpHb/vn37Uh6vXTPjeAcOHMhIfKgfC4ACAAAAACzL6/XK4/EktmAwmLFrV1VVSZJyc1P/zu90OiWJ6SEWxMgMG/P7/XI6nSouLlZxcXG2wwEAAADQgoRCIYVCIR05cuTYDiPn2GaWb/oKh8OJxTclZWxUhiTl5OQkrXVxvOrq6kQ7WAvJDBsrLS1N+qMGAAAAgEyp/dG0vLxcXq83a3G43e5me+5p27ZtYqpJKrXJjOPLtiL7SC8BAAAAAE5JnTp1kqR/jT45Tu00lO9+97umxYTGIZkBAAAAAEjPcJi/NbMePXpIknbv3p3y+K5duyRJ3bp1a/ZY0DQkMwAAAAAAp6TLLrtMkrRp06Y6x3bv3q1du3apY8eOJDMsiGQGAAAAAOCUVFRUpNzcXL3xxht1ji1cuFCSdOONN8rhaP5RImgakhkAAAAAgBbHMAyVl5dLkj7//HMZhlGnTfv27TV+/Hi9//77eueddxL7Y7GYpk+frg4dOmjq1KlmhYwmoJoJAAAAACAtw3DIMGEdi2/3J0nRaDRpv8vlarA8azAY1MKFC7V9+3bt3btXkjRr1iwtWrRIXbt21Y9//GNNmTIl0f43v/mNKioqVFRUpEcffVSdO3fWCy+8oMrKSs2fP1+FhYXN8OlwshxGqvQULC0ajcrj8SgSiVCaFQAAAECzqi3N+u/vzZbztNam9Rs7UK0tA2+ps//BBx/UQw89lPH+ysrKtGrVKkUiEXXt2lXXXnut2rZtm/F+kBmMzAAAAAAAWFY4HE76EbehURkno1evXurVq1ezXBuZRzIDAAAAAJCekXNsM7M/SW63mxHpqIMFQAEAAAAAgK2QzAAAAAAAALZCMsPG/H6/fD6fQqFQtkMBAAAA0MKEQiH5fD4FAoFjOwyH+RtQD6qZ2BDVTAAAAACYJVHN5N3/Mb+ayeU389yDlBiZAQAAAAAAbIVqJgAAAACA9Mye+vFNX9FoNGm3y+VqtvKssA9GZgAAAAAALMvr9crj8SS2YDCY7ZBgAYzMAAAAAABYVjgcTlozg1EZkEhmAAAAAAAaI0vTTNxuNwuAog6mmQAAAAAAAFshmQEAAAAAAGyFaSYAAAAAgLQccYcccfOmmZjZF+yHkRkAAAAAAMBWSGYAAAAAAABbYZoJAAAAAKARTK5mIqaZoH6MzAAAAAAAALZCMgMAAAAAANgK00wAAAAAAGk5DIccJk4zqe0rGo0m7Xe5XHK5XKbFAWtiZIaN+f1++Xw+hUKhbIcCAAAAoIUJhULy+XwKBAJZjcPr9crj8SS2YDCY1XhgDQ7DMIxsB4GmiUaj8ng8ikQicrvd2Q4HAAAAQAtWXl4ur9erbsvmytmmtWn9xqqq9fGVNyocDic99zAyAxLTTAAAAAAAFuZ2u/kRF3WQzAAAAAAApGeYXJrV1DKwsBvWzAAAAAAAALZCMgMAAAAAANgK00wAAAAAAGk5jBw5DPN+DzezL9gP3w4AAAAAAGArJDMAAAAAAICtMM0EAAAAAJCeIZOrmZjXFeyHkRkAAAAAAMBWSGYAAAAAAABbYZoJAAAAACAth+GQw8RpJrV9RaPRpP0ul0sul8u0OGBNjMwAAAAAAFiW1+uVx+NJbMFgMNshwQIYmQEAAAAAsKxwOCy32514z6gMSCQzAAAAAACN4DCObWb2J0lutzspmQFITDMBAAAAAAA2QzIDAAAAAADYCtNMAAAAAADpGY5jm5n9AfVgZAYAAAAAALAVkhkAAAAAAMBWmGYCAAAAAEjLYTjkMHHqh5l9wX4YmQEAAAAAAGyFZAYAAAAAALAVkhk25vf75fP5FAqFsh0KAAAAgBYmFArJ5/MpEAhI+tc0EzM3oD4OwzCMbAeBpolGo/J4PIpEInK73dkOBwAAAEALVl5eLq/Xq/PfeE3ONm1M6zdWVaVN1w3nuQcpMTIDAAAAAADYCtVMAAAAAACWFY1Gk967XC65XK4sRQOrYGQGAAAAACA9w2H+Jsnr9crj8SS2YDCY5X8IWAEjMwAAAAAAlhUOh5PWzGBUBiSSGQAAAAAAC3O73SwAijpIZgAAAAAA0jK7XCqlWdEQ1swAAAAAAAC2QjIDAAAAAADYCtNMAAAAAABpMc0EVsLIDAAAAAAAYCskMwAAAAAAgK0wzQQAAAAAkBbTTGAljMwAAAAAAAC2QjIDAAAAAADYCtNMAAAAAADpGY5jm5n9AfVgZAYAAAAAALAVkhkAAAAAAMBWmGYCAAAAAEjLYRzbzOwPqA8jMwAAAAAAgK2QzAAAAAAAALbCNBMAAAAAQFoOwyGHiRVGavuKRqNJ+10ul1wul2lxwJoYmQEAAAAAsCyv1yuPx5PYgsFgtkOCBTAyAwAAAABgWeFwWG63O/GeURmQSGYAAAAAABrDcBzbzOxPktvtTkpmABLTTCzLMAxt2bIl22EAAAAAAGA5JDMsYvPmzcrJyZHD4ZDD4VBOTo5KS0uzHRYAAAAAAJbDNBOLeOaZZxQKhdSqVStJUqtWrXT99ddnOSoAAAAAAKzHlsmMQ4cOqaKiQpFIRAcPHlSbNm10xhln6Kyzzsp2aCfkiy++UCQS0fjx47MdCgAAAACk5JDJpVll4vocsB1bJDMOHTqkN998U6+99po++OADffrpp4rH43XaFRQUqE+fPho0aJBGjhwpn8+XhWib7vnnn9e8efPUpUsX/fCHP9SECRPUrVu3bIcFAAAAAIAlOQzDMLIdRH2i0aimTZum3/3ud4rFYrr44ovVu3dvde3aVWeddZZat26tvLw8VVdXKxKJKBwOa/PmzSotLdU///lP9evXT/fff7+GDBmS0bhWrFih/fv3a9iwYRm53iuvvKKysjItX75c//jHP1RQUKC5c+fWe/1oNCqPx6NIJMKqvgAAAACaVXl5ubxer3q/+lc5W7cxrd9YdZXWjbiG5x6kZNmRGYsXL9btt9+ufv36ac6cORo8eLDy8/Mbff7nn3+uOXPm6M4771SvXr00Y8YMtW/f/qRiWrdunaZMmaK33npLDz74YIPJjE8//VQPP/ywNm7cqFatWunQoUO68847NXbs2Dptf/azn+lnP/uZJGnRokW69dZb9ZOf/ESfffaZ2rVrd1IxAwAAAEAmOOIOOeImTjMxsS/YjyWrmTzwwAN6+umn9fbbb2vBggW65pprmpTIkKROnTpp8uTJ2rJliwKBgK644gpt3rz5hOL5+uuv9eSTT+pPf/qTVq9enbZ9WVmZevfuLcMwtGbNGq1atUrPP/+87r77bo0ZM6bBc4cNG6b58+erqqpKixYtOqF4AQAAAABoySyXzAgGg3I6nVq+fLkuuOCCk75ebm6uJkyYoAULFmjChAnavn17k6/h8Xh033336amnntLQoUMbbBuJRDRs2DDl5uZq+vTpysvLkyT1799fkyZN0syZM1VSUtLgNQKBgL73ve9p165dTY4VAAAAAICWznLJjHPPPVcPPvigHI7MDin6/ve/r3nz5mnHjh1NPvfbsRQUFDTY9oUXXlA4HNaIESPUtm3bpGO33XabJGnq1Kk6evRog9cpLCxUhw4dmhwrAAAAADQLIwsbUA/LJTOuv/76Zrv2d77zHQUCgWa7viTNnDlTklL24/V6dc4552jnzp1atmxZvdfYv3+/tm/frmuuuaa5wgQAAAAAwLYsl8yQjj3Mn4zXXnstQ5E0zZ49e7R161ZJqrcs7Pnnny9JeuuttyRJa9asUVFRkZYsWSJJMgxD9957rx5++GF17NjRhKgBAAAAALAXS1YzGTRokNauXXvC5z/00EMaPnx4BiNqnA8//DDx+uyzz07ZpjZBsWnTJkmS2+1WOBzWj3/8Y1177bXq2rWrRo0apUGDBjV/wAAAAADQWIZDDsPECiNm9gXbsWQyY9OmTaqqqlKbNk2vYbxmzZqkpIKZvvrqq8Tr49fLqFVbH3n37t2SpPPOO++EEzfRaDTpvcvlksvlOqFrAQAAAIAkHT58WIcPH068P9mR80BzsOQ0kyNHjmjatGlNPm/27Nm64oorFI/HmyGq9KqrqxOv6yslW5tsqKqqOun+vF6vPB5PYgsGgyd9TQAAAACntmAwmPScUd8UeiCbLDkyQ5IeffRRXXLJJbr66qvTtq2srNTYsWO1ZMkSGYaR8UoojdWqVavE66NHj6ZMaNRWMWnduvVJ9xcOhxMjPSQxKgMAAADASZsyZYruvvvuxPudO3fK5/PJYfI0E1OntMB2LDkyo02bNvrzn/+sJ598Uu+++26Dbf/4xz/q/PPP15IlS+R0OjVw4ECToqzrzDPPTLw+cOBAyja1+zNRdtXtdidtJDMAAAAAnCyXy5X0nFHfFHogmyw5MmPRokUaNGiQrrrqKv34xz9Wfn6+Lr300qQ2x4/GuPDCCzV79mz17NlTPXv2zErc3bt3l8PhkGEYqqysVGFhYZ02lZWVkqQePXqYHR4AAAAAAC2CJUdm1FbyaNu2rV599VVNnTo1aZHMb4/GyM3N1UMPPaTS0tJEEiPdaI7mUlhYmIhh8+bNKdvULk5KtRIAAAAAdlI7zcTMDaiPJZMZ3+bxeDRv3jxNmjRJS5Ys0Y9+9CPdeuut2rdvny688EKVlpbqgQceUG7uvwaZnH766VmLd9SoUZKk9957r86xXbt26ZNPPlFhYWGj1gIBAAAAAAB1WT6ZIR0b8TBv3jxNnjxZixcvVm5urh5++OGk0RhmOXjwoCTVWzFl7Nixat++vebOnZtU3USSZs2apXg8rnvuuSdpsdAT5ff75fP5FAqFTvpaAAAAAPBtoVBIPp9PgUAgq3FEo9Gk7dtlY3HqsmQyY9WqVXX2tWvXTsuWLdN5552n//3f/9X9998vp9PZ6PMz4csvv9TKlSslSe+//75qamrqtCksLNQrr7yi/fv366677kpULyktLVUwGNQPf/hDTZ48OSPxlJaWavPmzSouLs7I9QAAAACgVnFxsTZv3qx33nnn2A7DYf4myev1JpWKDQaD2ftHgWU4DMMwsh3E8fr27avVq1enPFZZWambbrpJf/jDH3TOOeekbHPppZfq73//e8biicViuuSSS/TPf/5TVVVVif3t27dX//79tXDhwjrnbNiwQY8//ri2bt2q1q1bq6qqSrfccovuuuuuepMwjRWNRuXxeBSJRJJKswIAAABAppWXl8vr9cr/v8uV27qNaf3WVFep9KdXKBwOJz33uFwuKjnCmtVMtmzZolgslvKh/6yzztL//M//6Oabb1ZJSYm6dOmSdLympkZbtmzJaDxOpzNpAdLG6Nmzp+bOnZvROAAAAADgVFNbIhb4NksmM/bt26fevXvr3/7t3+pts3v3bl166aXq3r170v4vv/xS+/bta+4QAQAAAOCU4jCObWb2B9THkskMSdq0aZMaMwNm165ddfY5HJTwAQAAAACgpbJsMmPIkCHq0KFDk8+rrKzU22+/3QwRAQAAAACQfYcOHVJFRYUikYgOHjyoNm3a6IwzztBZZ52V7dBMY8lkxhlnnKElS5ac0LnxeFwdO3bMcETW5Pf75XQ6VVxcTEUTAAAAABkVCoUUCoV05MiRbIdyyjt06JDefPNNvfbaa/rggw/06aefKh6P12lXUFCgPn36aNCgQRo5cqR8Pl8WojWHJauZnGw1kkxXM7EaqpkAAAAAMEttNZOLX1mu3NanmdZvTfUBrfnZFaf0c080GtW0adP0u9/9TrFYTBdffLF69+6trl276qyzzlLr1q2Vl5en6upqRSIRhcNhbd68WaWlpfrnP/+pfv366f7779eQIUOy/VEyzpIjM2bPnp3V8wEAAAAAyKbFixfr9ttvV79+/TRnzhwNHjxY+fn5jT7/888/15w5c3TnnXeqV69emjFjhtq3b9+MEZsrJ9sBpHLuuedm9XwAAAAAALLlgQce0NNPP623335bCxYs0DXXXNOkRIYkderUSZMnT9aWLVsUCAR0xRVXaPPmzc0UsfkslcyIRCK2vj4AAAAAtFQOw2H6dioKBoNyOp1avny5LrjggpO+Xm5uriZMmKAFCxZowoQJ2r59ewaizD5LTTN5++23tWHDBj3yyCMZv/bGjRv18MMP69VXX834tQEAAADA7u72pF6XYn+KhSbRfM4991xdf/31Gb/u97//fc2bN09lZWXq2rVrxq9vNkuNzBgxYoTcbreGDx+uvXv3Zuy6M2bM0Lhx4zRjxoyMXRMAAAAAgExrKJFRU1OjmpqaE772d77zHQUCgRM+30oslcyQpHvvvVcjR47UhRdeqEcffVR79uw5oesYhqHXX39dfr9fS5Ys0dKlS1VYWJjhaLPL7/fL5/MpFAplOxQAAAAALUzZkSOafaBKrx6oksQ0Eyvo16+fbrvttrTtvvjiC91xxx266KKLNHz4cO3YscOE6MxlqWkmtX7yk5+of//+mjhxojp16qTBgwfrqquuUp8+fXTeeefp9NNPl9PpTDrn0KFD2r59u9avX693331Xb7zxhgoKCvTEE0/oxhtvzNInaV6lpaWnbIkiAAAAAM2rV36+euXna388rpnfJDSQXbm5ufrlL3/ZYJvq6moFAgGFw2GNHj1aX3zxha6++mpt3LixyYuIWpklkxmS1LlzZy1YsEAbN27USy+9pCeffFJffPGFHI5j2bnTTjtNHo9H8XhcVVVVikajkqS8vDxdeeWVmjZtmkaMGFEn6QEAAAAAgB117dpVHTp00P3336/169frBz/4ge65556kH7n/9Kc/6ZNPPlEwGNTkyZMlSePGjdMrr7yi22+/PVuhZ5xlkxm1LrjgAj333HN67rnn9PHHH2vt2rXatm2b9uzZo+rqauXn58vj8ahLly7q0aOH/H6/XC5XtsMGAAAAgJbF+GYzsz8k+cUvfqFzzz1X+/btkyQtWbJEixYt0t///ncVFBRIklasWCGHw6Grrroqcd6UKVP0n//5nyQzsqVbt27q1q1btsMAAAAAAMB0f/3rX9W9e3eNGzdOHTp00I4dO/S73/1Ozz//vCZNmiRJ2rZtm6Rj1UtqdenSRRUVFVmJubnYKpkBAAAAAMCpasWKFVq+fHnSbITrrrtON910UyKZsX//fklS27Ztk85t06aNeYGagGQGAAAAAEDTItGU+8vLyzXT6zW9wgjVTOrKy8urs6zCWWedpZycfxUqjcVikpRYb7JWdXV18wdoIsuVZgUAAAAAAHUdPXpU4XA4ad/GjRtVWFiYeB+Px+ucd/DgQeXl5TV7fGZiZAYAAAAAADZwww03qHfv3ho+fLjOOOMMff755/rLX/6iV199VZJUVlam7du3S5I++ugjde/eXdKxCicDBgzIWtzNwWEYBmvE2kw0GpXH49G5554rp9Op4uJiFRcXZzssAAAAAC1IKBRSKBTSkSNH9Omnn6rfH95VbuvTTOu/pvqAVt12uSKRSFLp0VPZkSNHVFRUpDfffDOx74477lAkEtFnn32m9evX65JLLtGIESP08ssv6+WXX9auXbs0ZswYrVixQuedd14Wo88sRmbYWGlpKX/UAAAAAJpF7Y+m5eXl8nq92Q4HkvLz8xMjMbZv366BAweqX79+ko6tiXHo0KHElJOysjJdcsklysnJ0ZNPPtmiEhkSyQwAAAAAAGwjJydHI0eOrLO/devWat26deL9rFmzdN9998ntduvss882M0RTkMwAAAAAAKRnOI5tZvaHk1K7ZkZLRDUTAAAAAAAsIBKJ2Pr6ZiKZAQAAAACABbz99tt64IEHmuXaGzdu1O23394s184GkhkAAAAAgLQccYfp26lmxIgRcrvdGj58uPbu3Zux686YMUPjxo3TjBkzMnbNbCOZAQAAAACARdx7770aOXKkLrzwQj366KPas2fPCV3HMAy9/vrr8vv9WrJkiZYuXZqodNISsAAoAAAAAAAW8pOf/ET9+/fXxIkT1alTJw0ePFhXXXWV+vTpo/POO0+nn366nE5n0jmHDh3S9u3btX79er377rt64403VFBQoCeeeEI33nhjlj5J8yGZAQAAAACAxXTu3FkLFizQxo0b9dJLL+nJJ5/UF198IYfj2PSb0047TR6PR/F4XFVVVYpGo5KkvLw8XXnllZo2bZpGjBhRJ+nRUpDMAAAAAACk5TAccphYLtXMvqzsggsu0HPPPafnnntOH3/8sdauXatt27Zpz549qq6uVn5+vjwej7p06aIePXrI7/fL5XJlO+xmRzLDxvx+v5xOp4qLi1VcXJztcAAAAAC0IKFQSKFQSEeOHMlqHLUjDmq5XK5T4mE9lW7duqlbt27ZDsMSHIZhGNkOoqkOHz6sv//97zp48KCGDh0qSdqwYYO8Xq9OP/30LEfX/KLRqDwejyKRiNxud7bDAQAAANCClZeXy+v1asCMlcptfZpp/dZUH9DKsQPq7H/wwQf10EMPmRYHrMl2IzNKSkp07733at++fercubO2bdsm6dhiJ7feequKiop0yy23ZDlKAAAAAGhZHJIcJv4UXjvJJBwOJ/2Ie6qOykAyW5VmXbx4se644w71799f06ZNU7t27RLHLr74Yr322mv6xz/+oeXLl2cxSgAAAABAprjd7qSNZAYkm43M+O1vf6tXXnlFN910kyTp9ddfr9Pm17/+tUaNGqUrrrjC7PAAAAAAwDLu9qSeku7ILUi535nbOuX+aE1NxmICMsVWyYyqqqpEIqM++fn5OnjwoEkRAQAAAMApwnAc28zsD6iHraaZfOc732lUuy+//LKZIwEAAAAAANliq2RGPB7X1q1bE+9TFWL57//+b3Xu3NnMsP5/9u48Pqrq/v/4e2aSTBJIAgECBAJhVUPFjaA2VVFUlloQxCqgiKgojXVFkVabqFQWAeWrYF0q1Fgt7iLUpVjxJ1YkFhARERDRIBC2kITss/z+oERjzg0Ekntnwuv5eMzD5HPPveczBMPMZ849HwAAAAAAYKOwus0kMzNTF1xwgWbMmKGhQ4fK5Tq47Gjfvn1asWKFnn32WS1atEjLly93OFMAAAAAaFpcQZdcNt76YedcCD9hVcwYPny4PvzwQ40aNUput1tut1sxMTGqrKysHjNr1iz17dvXwSwBAAAAAEBjCqvbTCRpzpw5WrBggbp27Sqfz6eKigoFg0H17t1bb731lm677TanUwQAAAAAAI3IFTRtPBEmtm7dqvz8fHXo0EEdO3Z0Oh3bFBUVKSEhQYWFhYqPN7dbAgAAAICGsG3bNqWkpOi8Jz5RRExz2+b1lR3QhxPO5n0PjMLqNhNJ8vv98ng8kqTU1FSlpqY6mxAAAAAAALBVWN1mMnv2bMXExOh3v/ud06mEhPT0dKWlpWnu3LlOpwIAAACgiZk7d67S0tLUr18/p1MBagmrlRlz5syR3+9XcXGx06mEhNzcXJZbAQAAAGgUmZmZyszMrL7NhG4mCCVhtTIjNTVVu3fv1nPPPVfnuPHjx9uUEQAAAAAAsFtYrcy49dZblZWVpdmzZysyMtJy3EcffWRjVgAAAAAQeu5JTjPGXZ5oYzy62QnGeGFVaYPlBDSUsCpmtGnTRh07dtRpp52mjIwM9erVSwkJCXK5flx+lJ+fr40bNzqYJQAAAAA0QQHXwYed86GGF154QT6fT5J04oknqm/fvsrPz9f48eP1/vvvq3Xr1srKytK1117rcKaNL6yKGUOHDlVhYaGCwaDWr19fHf9pMSMYDNb4HgAAAACApmD16tV64oknNGHCBHXr1k2SNHjwYK1fv14PPPCAunTpolmzZqlly5a69NJLHc62cYVVMaNt27bq0aOHxo8fr4gIc+qFhYWaNGmSzZkBAAAAANC4IiIi9O677yojI0OStHTpUq1evVpTpkzRXXfdJUk655xzNHr0aIoZoSQpKUmTJ0/WwIED6xyXk5NjU0YAAAAAcHygm4nzVq5cqalTp1Z//+qrr8rlcmnMmDHVsbZt2yoYDDqRnq3Cqpgxbdo09ezZ87DjZs2aZUM2AAAAAADYJxAI1Ph+9erVSkpKUseOHWvES0ub/qatYVXMOPvssw875rvvvtP27dttyAYAAAAAnHV321TLY25PrH2JwBZVVVXy+/3yeDzaunWrPvvsM11xxRU1xmzatEnt2rVzKEP7hFUx40isXr1ajz/+uK688kqnUwEAAAAAxwT8pZqRv7VWfFL7E43jg0G/+UJWcdhuwIABmjhxosaNG6f77rtPwWBQ48aNqz6+bt06XXXVVXrhhRcczNIeYVXMeOCBB+o87vP59NJLLx0X9wcBAAAAgJ3YM8N5kyZN0qhRo3TKKacoMjJS06dPV//+/bVx40aNHTtWn376qYLBoO6++24tXrzY6XQbVVgVM7Kzs+VyuSyLFYeOtWjRwubMAAAAAABoXFFRUXrllVdUXFwst9utZs2aSZK6du2qF198sca4pi6sihkul0t/+tOflJqaWiPu9/u1detWvfnmm7r44ov1m9/8xpkEAQAAAABoZHFxcTW+j4iIUOfOnR3KxhlhVcxIT09XVlaW5fH7779fU6ZMUUFBgY1ZAQAAAEDTx20mznvhhRfk8/kkSSeeeKL69u2r/Px8jR8/Xu+//75at26trKwsXXvttQ5n2vjcTidQH4899lidx10ul+677z49/fTTNmUEAAAAAIA9Vq9erd/97nf64osvVFVVJUkaPHiw3nvvPWVlZWnmzJl66qmn9MYbbzicaeMLu5UZR6KioqKRMwkN6WQqs9IAACAASURBVOnp8ng8yszMVGZmptPpAAAAAGhCVu7frNz9m+UPBJxOBf8TERGhd999VxkZGZKkpUuXavXq1ZoyZYruuusuSdI555yj0aNH69JLL3Uy1UYXVsWMI3HgwAEVFRU5nYYtcnNzFR8f73QaAAAAABziryqu9znuiNgjGte3RXf1bdFdRVWleuTbJdxmEgJWrlypqVOnVn//6quvyuVyacyYMdWxtm3bHhcdPsOqmPHT/rkmpaWl+uSTT3TZZZfZlBEAAAAAAPYI/GyVzOrVq5WUlKSOHTvWiJeWltqZliPCqpixYMGCw4459dRT9Yc//KHxkwEAAAAAwEZVVVXy+/3yeDzaunWrPvvsM11xxRU1xmzatEnt2rVzKEP7hFUxIzIyUk8++aS6dOlS65jL5VLr1q2VlpbmQGYAAAAA0MQFJJed22ewVUctAwYM0MSJEzVu3Djdd999CgaDNe5gWLduna666iq98MILDmZpj7AqZvTp00djx451Og0AAAAAAGw3adIkjRo1SqeccooiIyM1ffp09e/fXxs3btTYsWP16aefKhgM6u6779bixYudTrdRhVUxY/78+U6nAAAAAKAJm9i6fZ3HZ+7ZYVMmDWNS+xNrxY50A1CEnqioKL3yyisqLi6W2+1Ws2bNJEldu3bViy++WGNcUxdWxYyePXsa47m5uaqoqFBycrK6du1qc1YAAAAA4AxPZJzlMZc72hivLN1mjEdEtTLG/f7yg9ejm0nIiIur+XOPiIjQpk2btGPHDnXv3l1nn322Q5nZx+10AvVx3XXXGeMrVqzQ0qVLNXXqVA0ePFgbNmywOTMAAAAAAJwTDAYVCAT09ttv67HHHnM6nUYXVisztmzZYoz//ve/r/56//79GjlypN5++2270gIAAAAAwFEXXXSRJGnr1q0aMmRIjffJTVFYrcxwuQ6/zCghIUE7d+60IRsAAAAAOH64gj/eamLPw+lnHJpyc3N19tlny+v1yuPx1Hp069ZNMTExTqfZ6EJ2Zcbnn3+uzz//vEZs586dysnJUTBY+291MBjU7t279cYbbygxMdGuNAEAAAAAsMW3336rCy+8UBUVFWrfvr3y8/PVrl07SVJlZaV27typSy65RE888YTDmTa+kC1m+Hw+bd26VcuXL9f7779fHa+rNWswGFTr1q21dOlSGzIEAAAAAOnOVknGuLdZF8tzHvr+03rN8YdOZ9ZrPJqmhx56SBdffLGefvpptWjRQuedd54+/PDD6uObNm3STTfdpP379ys5OdnBTBtfyBYzzjjjDJ1xxhmSpK+++kqXXXaZysvLdc011xjHu91ude7cWZdccgkrMwAAAAAclXBrvRqTcIrlMV/lnvpdq/npxnhlZZGk5XQzCQG5ublasWKFoqMPdqqJiIhQQUGBWrZsKUnq0aOHcnJydOutt+rll192MtVGF7LFjJ866aST9Oyzz2ry5MnKyspyOh0AAAAAAGzXqlWr6kKGJPXq1UvvvPOORo4cWR1LTk7W/v37nUjPVmGzAWh6erqGDRvmdBoAAAAAADiiRYsW8vl8KisrkyT1799fM2bMkM/nqx4TDAa1bds2p1K0TdgUMzwej2655RbL4//617+Uk5OjTz75xMasAAAAAOD44Aq4bH+gpi5duqhv377q3r27SktLNXDgQG3fvl0jRozQ5s2bVVRUpIkTJyo+Pt7pVBtd2BQzjkQgENDbb7+txx57zOlUAAAAAADHiddff1033nhjrfh///tf9e3bV3l5eQ0yz8iRI/XVV19p//79qqyslNfr1cyZM7Vo0SKdcMIJatmypR599NE6G2c0FWGxZ8aRuOiiiyRJW7du1ZAhQ/T73//e4YwAAAAAAE3dnj17dNNNN2nQoEHVsVdeeUWLFy+Wy+VSbm6u/H5/g8x1xhlnaMuWLSovL1eLFi0kSVdffbUKCgo0Y8YMFRcXa/To0ZowYUKDzBfKwq6YkZubq1tuuUWrVq2qcV/QT/Xp08fmrAAAAACEojsSrJfbzy4sqhWb2Lp9nddrqG4nf0zNMMaD/grzCS5Pg8yLhvfoo48qNTW1RmzEiBEaMWKEli1bpgULFjTofO3b1/47esstt9S5LUNTFFbFjG+//VYXXnihKioq1L59e+Xn56tdu3aSpMrKSu3cuVOXXHKJnnjiCYczBQAAAHC8mLV3lzF+Z2IrYzw6/sR6zxERZb6W29Pc8pxIr9cYd1kURqKiepmvo30Hz5PNrVkV+ntmvPzyyxo0aJCWL1/eoNctLi5WXFycY+eHg7DaM+Ohhx7SxRdfrJ07d2rr1q3q27evvv32W3377bf64YcftGHDBh04cOC4aEMDAAAAAHBOfn6+Nm/erIwM8yqbYzFw4MBjOv+nt7w0VWFVzMjNzVVOTk71vUEREREqKCioPt6jRw/l5OQoKyvLqRQBAAAAACEkMzNTf/7zn+sc89VXX2nEiBE65ZRTdP755ys9PV1PPvmkgsGg5TmPPPKIbr/99oZOV5KOubVqQ204GsrC6jaTVq1aKTo6uvr7Xr166Z133tHIkSOrY8nJyazMAAAAAIAG5grafJvJMczl8/m0bNkyTZkyRR9++GGdH3h//fXXOvvss/Xb3/5Wq1evltvt1vr165WRkaH//ve/euqpp2qd8+KLL+rSSy+t8f60IeXl5el3v/td9bYK9bFjx45jLoaEg7AqZrRo0UI+n09VVVWKiYlR//79lZ2drcsvv1wREQefSjAYPC5+cJKUnp4uj8ejzMxMZWZmOp0OAAAAgCbk411L9Z9d78sXMDdeCFU333yz1qxZoxNOOEHNm1vvKSIdfP94/fXXy+/3a/bs2XK7D968kJaWpttvv11ZWVkaOnSofv3rX1efs337dn333Xc1PlRvDE8++eRRnRcMBuVyhf5+I8cqrIoZXbp0Ud++fZWfn69NmzZp4MCBGj9+vEaMGKGZM2cqKSlJ999/v+LjrXcsbkpyc3OPm+cKAAAAhKq723ZzOoVGkZF0oTKSLtT+yn3689rGuZ2iMTz++OPVX2dnZ2vJkiWWY1esWKHly5fr6quvrlX4uOqqq5SVlaXZs2fXKGYsWbJEa9eu1dixY6tja9eu1bZt2zR27FhNmDBBZ5555jE9B6/Xq4yMjKPaxDM/P1+ffvrpMc0fDsKqmDFy5EjNnTtXbrdblZWVatGihWbOnKlrrrlGb731VvW4n/7lBQAAAHD8MrVfrYvbE9tgc3siW2hG/je14vd26Wd5TlSzTsZ4wH/AGA/6Sy2vFentaJ7Da55Dwao6466AS66AjbeZ2DDXa6+9Jungqvef69q1qxITE7Vs2TIVFRVVf5B8ww036IYbbqgxNjU1Vb/61a8arA3rqlWr9MILL6i8vFxDhgzROeecc8Tn+v3+o7o9JdyEVTHjjDPO0JYtW1ReXl69CejVV1+tgoICzZgxQ8XFxRo9erQmTJjgcKYAAAAAgFC3Zs0aSVLnzp2Nxzt16qR9+/Zp9erVOu+88yyv4/f75ff7a8UDgUCN/x6pk046SQ8++KB8Pp8WLVqk2267TR06dNDo0aOVnJxc57kej0fdujXN1UI/FVbFDElq3759rdgtt9yiW265xYFsAAAAAADhavPmzZKkxMRE4/FD8R07dhiPL1iwQAsWLNC2bdv06quvqqCgQFlZWSooKNBTTz1VfbvHwIEDddppp2nSpEk6/fTTjzi/iIgIDR8+XMOHD9f27duVk5OjHTt26Nxzz9WQIUOq9478uXnz5h3xHOEq7IoZfr9fHo/H6TQAAAAA4LjiDrrktrGbyaG5iopq3irk9Xrl9XobZI59+/ZJkmVXkkPxgoIC4/GxY8fW2Dvjpy6++OJjT/AnkpOTNWnSJEnSsmXLdPfdd6tZs2a68sor1atXrxpj61MwCVdupxOoj9mzZysmJka/+93vnE4FAAAAAGCDlJQUJSQkVD+mTp3aYNcuKSmRJMsVDoc+SN+/f3+DzdkQ+vXrp9mzZ+uuu+7S8uXLlZmZqaefflrFxcVOp2absFqZMWfOHPn9/uPqBwQAAHBnqyRjPCLKvCza5arj8yqXxQrXYO17vSXJV1VYZ24zd/9Q53HgeJfdc5jTKYS9vLy8Gl0cG2pVhiS53W7jXheHlJaWVo8LRfHx8brxxhslSe+++65OPvlk9evXTzfccIMyMjIczq5xhVUxIzU1VatXr1bLli3rHDd+/Hg99dRTNmUFAAAAoKmIiKr7vYaJt3lXY9zjaW6MR0RazxERVXuPQEmKiGhrjFdVfmt5rSjvyca4yxVpjHu8vzDHlX/wvKBLLhtvMzk0V3x8fI1iRkOKi4urvtXE5FAx4+dtW0OF3+/XokWLNH/+fL3zzjvy+Xx67rnnlJ+fr7ffftvp9BpVaJaXLNx6663KysqSz+erc9xHH31kU0YAAAAAgHDVqdPBNrWVlZXG44duQ+nY0dzm1ilffvml7rzzTiUnJ2vEiBFavHixkpKSdM8992jjxo1NvpAhhdnKjDZt2qhjx4467bTTlJGRoV69eikhIUEu14/Vwfz8fG3cuNHBLAEAAAAA4aBXr15as2aNdu3aZTyen39wVcoJJ5xgZ1pGRUVFevHFF/Xss8/qs88+UzAYVGRkpC699FKNGzdOgwYNCtnbYRpDWBUzhg4dqsLCQgWDQa1fv746/tNiRjAYrPE9AAAAAODYuQJuuQL2vVm2Y65zzz1Xf//737Vu3TpdcsklNY7t2rVL+fn5Sk5OdrSY8e9//1vPPvusXn/9dZWXlysYDCotLU3jxo3T1VdfrTZt2jiWm5PCqpjRtm1b9ejRQ+PHj7fcbbawsLC6XQ0AAAAAAFaGDRumzMxMLV68WPfcc0+NY6+//rok6YorrrD9A/Pvv/9eCxYs0IIFC/Tdd98pGAwqLi5O1113na677jqdeeaZdZ5/1VVX6fnnn7cpW2eEVTEjKSlJkydP1sCBA+scl5OTY1NGAAAAAIBQFAwGtW3bNkkHiwOmVfxt2rTRhAkT9Nhjj2nZsmXq16+fpIMba86bN09t27bVvffea3fq6tmzp6qqqhQMBnXOOedo3LhxuvzyyxUbG3vYc4PBoN5//30bsnRWWBUzpk2bpp49ex523KxZs2zIBgAAHHJvl37GeLOEXxnjXm+6MR6IP81yjqKkVsb43nYBY/xAfIXltcpizBu9BTzW7fmsuP3mVqcev3l5tFU8osqiZaok8zMPDbP6fV07WLrWOLa8/BNz/MDnltevLDO3fvXX0TJ25p4dlscQ3u5MNP/fYNWmWJICAevfBQ/v+r7eOWT1HFrvc5oKp7qZFBUV1Yh7vd4627NOnTpVr7/+ur799lvt2bNHkjR//nwtWrRIXbp00fDhwzV58uTq8dOnT9f27ds1bNgwPfjgg+rcubP+8pe/aMeOHXr11VeVmGj996uxVFZW6uyzz9a4cePUvXt3uVwuffbZZ4c9LxgM6tVXX7XcA6QpCatixtlnny1Jqqio0CeffKKysjINGjRIkrR27VqlpKSoZcuWOvfcc51MEwAAoEEV3WgoGEiKtCiAeHzWbzaiys3FlOhS8znx+63fCHoLm/6LZRx/3BHWn3x7LFqqBgNVxrgrMtoYj/R2sp7fZT7H7ze3D42ONheH6+L2tDPGg27z3FZxu6SkpNT4PisrS9nZ2ZbjJ0+eXKNYcTgxMTF65ZVXtGbNGv3nP//RunXrdPXVV+sf//iH4uLijjbtY9a7d29t27atenXJkQgEAlq5cmUjZhU6wqqYIUkLFizQxIkTVVBQoM6dO2vLli2SpPLycl177bUaNmyYrrnmGoezBAAAAAA0hLy8PMXHx1d/X9eqjGNx6qmn6tRTT22Ua9dXx44d9cQTTxzVuRMnTlRycnIDZxR6wqpvy5IlS3T99dcrIyNDs2fPVuvWrauP9e3bV2+88YY+//zz4+L+IAAAAACwkzvgsv0hSfHx8TUejVXMCCUdO3Y86nPj4uJ00kknNWA2oSmsihkPP/ywcnJy9Oabb+rWW29V8+bNa42ZNm2a5s6d60B2AAAAAAAcu3feeeeYzl+6dGkDZRK6wqqYUVJSopEjR9Y5JioqSmVlZTZlBAAAAABAw/rpbTVOnB8OwmrPjBYtWhzRuL179zZyJgCApmJi6/bGeFRMB2O8eWJ/y2tFx/Qzxks6ZBjj36b5jPH8ZOt/x4pbFRnjVdHm7hz+yPp353D5zRtBenzWLxsuqPcsQOO5r9tFxrjlJo0u80aqLlek5RyBQLnVEWPUqvuKr7LAcg7p6LpthJo7WyUZ495Y600wJemhvMN3bjhWf0yt/fvZ6u8DgNASVsWMQCCgzZs3q3v37pIOtp35uSeeeEKdO3e2OzUAAI5rHz480xh3W7QhdfvMbxYiqoot54gsN795jC2JMcYTd1vvQN9tnfl+6+bb/mOMlx1YbHmtwl3/tDw2I/8by2NOmZxyujHePHGAMe6NPsv6YlGpxnAg4UzztSzisXHbLafwVaw2xstLzW90KyvC/80/rMUnXWyMlxaa/55IUsBv/t0RHdfDGI+ItG6G7PeZC8qR0eYCeESEuaWnx2Pd6tMl8+/HyIie5vGu2rfeH04w9gRjvKRVG3O89GBhz6nWrIBJWN1mkpmZqQsuuED/+Mc/VFZWJpfr4F/uffv26Z///KdGjBihW2+9VXfddZfDmQIAAAAAgMYSViszhg8frg8//FCjRo2S2+2W2+1WTEyMKit/XFo7a9Ys9e3b18EsAQAAAABAYwqrlRmSNGfOHC1YsEBdu3aVz+dTRUWFgsGgevfurbfeeku33Xab0ykCAAAAQJPjCrhsf0hSUVFRjUdFRYXDfxIIBWFXzJCkMWPGaOPGjdqyZYs++eQTff/991q9erUGDx7sdGoAAAAAgAaUkpKihISE6sfUqVOdTinkvfbaa06n0OjC6jaTn0tNTVWHDh0UGWm90zRwNB5Mu84Y90S0NMaDQfPu6AeZdzUPBs0dBuq+lpnVrttWc9S1O7vbHWs+xx1tcS1zXJI8bnMHIrentcW1zJv4uTzmP3dJCkS1Ncb9XnM7qqpY8/OTpPLm5l+JJRadrcpjrbtEFMeZPzGoiDb/fMtjzJuTVVqMP3jMPIfPaz7HH2nunCFJ/gjzc7HqhBGwGC9Jfo/FORZxuWpv5ixJR7Pnl8Wl6mTeSg4AQsOMM+YYolavLcxdjSTrTSUDgRKLuPVmvH6/uTNLwH/A8pyj8aduA2vFPJHWrwnQdOXl5dVoNer1mjdxxo+mTJmi4cOHO51GowrLYsbatWs1ffp0vf322yosLFSrVq00dOhQTZ48WV27dnU6vWNy7bXXqnPnzsrOznY6FQA4Lix9Jcd8wKKa4g5YL2q02nXdFTDvsu+pMv8zHFFhXXCMLWpmjHfeYG5xePLH1kWs0q23Wh7L/voly2NNm7mjh3Vckh5rjEQazdS8VY7N/YdO5m4mLZPGWp7jiUgxxpvFJRvjMbF7LK/l9+8yxoNBc0E5KOvCrWVxwOJabre5u47VdSTJZfEBQ7iJa93PGK+z+BE0F2ysPnSJ8Fp3IIlubu5a4veZCzbBoPXvTW+s+VpuT4IxbtW1pK4Pr6w/QDJ3LXFFmv9fkKSgReHJSlSZ+YOSqPKDhTJ30CW3jR1GDs0VHx9fo5jR1B04cECzZ88+6vOLior0xRdfNGBGoSnsihnz58/XhAkTamz6uWfPHv31r3/Vyy+/rFdeeUUXXnihgxkevbfeekvPPfec7rvvPqdTAQAAAAA4JDs7Wy6XS8HgUSw7lao7fzZlYVXMWLlypW688UZ17NhRN954o84991y1atVKeXl5+vjjjzVv3jz99re/1Zo1a9Spk/kTqlC1d+9evf3220pJMX/6AAAAAABo+po3b67IyEjddtttGjRoUL0LE4WFhRo9enQjZRc6wqqY8dBDD+myyy7TggULatwn1bNnT/Xv31+33nqrRowYoZkzZ+r//u//HMy0/rKyspSdna1//vOfTqcCAAAAALX9pMOIXfMdrxITE/Xggw8qKirqqM7v1atXA2cUesKqmLFx40atXr3acsOXhIQELVy4UOeff36j5fDBBx+ouLhYQ4YMabBrvvTSSzr33HOVlJTUYNcEAAAAAISnpUuXHnUhQ5Kee+65BswmNIVVMaNdu3aH3bk2MTFRLVqYOygci1WrVmny5Ml67733lJWVVWcx45tvvtH999+vL774QjExMSovL9dNN92k8ePH1xq7c+dOLVu2TPPmzWvwnHHQL9/41Bgf8kw3mzMBAAAAgMM71pUV69atU8+ePRsom9AUVsWM+Ph4+f1+eTzWuz4Hg0G53bV3ml+2bJn69etX7zn379+vp556Srt27dKKFSsOO37NmjU677zzNGTIEK1cuVKRkZH6+OOPNWDAAOXm5urpp5+uMf7ee+/VtGnT6p0XGlezhJHGeMC33RwPFFhey2qHdOsTzLt318Vqt3OrHdJdbnM3BKmOnbI95h2kXRa7dEuSP8bcNrU03lxwLGlh/pW0v5X1jt9FCeY/39LmZcZ4eWyR5bWsWp1atSetq9VpwGPRNs9d/5+vFZdFVw233xy36pwhSd4D5ra40SXmHdWbFVv/HUrZav75Jm/ca4z7975ojO/Oq/t2wZm7f6jzOI7E8dqxBE556HvzBwx1ubttqjEe16q/MR4Z1bGOq5l/P1p1IAkGrP8NDwbN/2ZYtUWv8n9vzqiOjiXuCHOHDFl0wnC56vjQz6KNu1W79ohI633cvNGnWVzL/G+JLF53WI5XXW1mzX++zROtX48ELbqDBCPMefmizZ1nJMkfaf63tDLW/OdYGW3+O+dquJcD8vjqv0FkTLH572+V1+J5BA+uEHCqmwnqj9asIWbMmDGaO3eubrnlFssxTz75pEaMGFErfscdd2jVqvq3IktISNDdd98tSdq2bZsWLlxoObawsFBDhgxRRESE5s2bp8jIg78MMjIydNdddyk7O1sZGRkaO3asJGnBggUaPHiwWrduXe+8AAAAAABND61Zj0xYFTNKSkq0YMECFRQUqEuXLrWOr1+/XosWLdKkSZNq3CO0Z88erVu37qjm/OnOsdHR5qruIX/5y1+Ul5en8ePHKy6uZkV33Lhxys7O1r333qvRo0crMjJSf/vb35Sbm6sxY8ZUjystLdVDDz2kmTNn6sknnzwudqEFAAAAACtFRTVX1nq93sNuPxDuaM16eGFVzLjtttu0f/9+rVmzxvjDOfSDHjduXK24HT/MZ555RpKMt7OkpKSoa9eu2rJli5YuXapBgwbp+eefV1lZzaXw/fr10/Dhw3XLLbeobVvzEn0AAAAAsJsraG83E9f/bjNJSal529OhTpBNFa1Zj0xYFTNatWqlnj17asCAAfX6gRYUFOjxxx9vxMyk3bt3a/PmzZKktLQ045hf/OIX2rJli9577z0NGjRIHTp0qDUmIiJCiYmJ6t69e6PmCwAAAADhIC8vT/HxP+7f1tRXZUi0Zj0SYVfMeOSRR3TWWWfV+9xly5Y1fEI/8eWXX1Z/bSpSSFJycrIkHfUtLz93PC63AgAAANC4fP4K+QI/bhJaUXXAwWwONoL4aTHjeEBr1sMLq2LGAw88YLnq4adefvllXX755TVi06dPb6y0JEn79u2r/vrn+2Uccuh/wF27dlleZ+vWrUc85/G23AoAAABA4/towywt+6p2x0V3QHLbeJtJAzaBCzuHW1lRUVGhTz75RGVlZRo0aJAkae3atUpJSVHLli2bfFtWKcyKGRdddNERjZs6dWqtYsaAAQMaI6VqpaWl1V9bVdAOrZooKSlpkDmPx+VWR+Pj/veZD/SXXHHvGQ/d0eJyYzwy2tx5xuWx3hzWE2nR0tSiFZonwlwMc7mt25d5LFq3uT3ma9XVBs7jaWWOq705L5d1xdhVZX7uUWXm51IVbW4bF1tS168q8599pM/iWhYtSCXJZ9GCNegyb7zksWiBKkmRleacI6vMeUWXm8dHldc1h/nFRIzFhydxu/dbXsu9/z/G+IGi1yzPydpgbqlqzaqd66R6xgEcT2bkb3U6hSZhYhvzymGp/q2u70iw/oR8dmHtFuh/TM0wjnV7rF+PWL1O8sZ0M8b9/mLLa0VEWLx+s2jZ6nZZtx/3RnYyxqM9Ft0JI9oZw8E6Xjv6YswtzstamPPyR1gXF0qbm+MHWpjfN3j+13X+1F/epZPP/LGD5IGSHZr/8hmW88BeCxYs0MSJE1VQUKDOnTtry5YtkqTy8nJde+21GjZsmK655hqHs2x8YVXMkKQdO3ZoyZIlys/Pl99f+42HU21oYmJ+fINUVVVlLGhUVR3sCR4ba/2Luz6Ox+VWAAAAABpXhMerCM+PBY/KSutiEey1ZMkSXX/99fr1r3+tCy64QH//+9+rj/Xt21dvvPGG7rjjDr3//vvq37+/g5k2vrAqZixfvlyDBw9WSUlJnS1qnGhD067dj1XXAwcOKDExsdaYAwcOflxKlxIAAAAA4cYVPPiwcz7U9PDDDysnJ0cjR46UJL355pu1xkybNk1XXnklxYxQ8oc//EEnnHCChg8frrZt28rjqb1cu7CwUJMm2b80+aSTTqruA7xjxw5jMWPHjh2Sjo+dZQEAAAAADaukpKS6kGElKipKZWVlNmXknLAqZuzcuVPr169XRETdaefk5NiU0Y8SExPVu3dvff7551q/fr2xYHGo48n5559vd3oAAAAAgDDXooV5T5Wf27t3byNn4jzrneVCUI8ePQ5byJCkWbNm2ZBNbVdeeaUk6aOPPqp1LD8/Xxs3blRiYqIuvvhiu1MDAAAAAIS5QCCgzZs3V39v2n7hiSeeUOfOne1MyxFhtTIjOTlZO3furLE/hUl+fn6jzH9oqU4gYO4RNH78eM2ePVsLFy7UtGnTamz0C5wHcgAAIABJREFUOX/+fAUCAd155501Ngs9Funp6fJ4PMrMzFRmZmaDXLMpsupYUpfZ+wuM8YmtzR09zP0pDgq4zbtxu61OcplrjG6LHbcl6x2/rbqWeDxtLK/l8Zj3dHFHJBvjAa/5z0Sy3o27PM78/0BltHm/G1fA+obJ6DLzH6S3whz3+Kz31Dm0g/fPWXUNiT5g7n4iSVFlFcZ4ZKm5o4iryryxlstfaDlHwJdnjJeXrTTG9+z+p+W1JGn6zk2GqLmzDwAgvNS3Y0ldTB1L6vLnrR/Xe457OpxsjJfuX22Me5t3tbxWpWuLMW7VDc7lsn6L5InYbIy7Pea2IVav0SIjUy3ncBVXGePNd1t0yHNbNwRItOjkIpfV+5Gar22W73hVH+98Xb7AwRdJ7oDL5tas9u+FGOoyMzN1wQUXaMaMGRo6dGj1fpH79u3TihUr9Oyzz2rRokVavny5w5k2vrAqZjzwwAO69dZb9fDDD6tTJ3NbJMncmvVY7d27t/ovxMcffyyfz1drlUhiYqJycnI0bNgw3XzzzXryyScVGRmp3NxcTZ06VYMHD27Q/Txyc3PpZgIAAACgUfyq/WX6VfvLtL9ilx7473Cn04Gk4cOH68MPP9SoUaPkdrvldrsVExOjysrK6jGzZs1S3759HczSHmFVzGjfvr3uv/9+9evXT3369FHPnj1rtUBt6Nasfr9fZ555pjZs2KCSkhJJ0r///W8lJycrIyNDr7/+eo3xAwYM0IoVK/TnP/9ZZ511lmJjY1VSUqIHHnhAN998s3HTUgAAAACAWVFRzdVAXq9XXq/XYnTTN2fOHJ1xxhmaMmWKNm/eLJ/v4MqZU045RX/+8581ePBghzO0R1gVM1auXKmBAwdq//792rp1q+W4hmzN6vF49Nlnn9XrnN69e2vhwoUNlgMAAAAAOM0VcMll460fh+ZKSUmpEc/KylJ2drZteYSiMWPGaMyYMdq6davy8/PVoUMHdezY0em0bBVWxYx77rlH3bt31/Dhw5WUlBRSrVkBAAAAAA0vLy+vxu31x/OqjJ9LTU1Vamqq02k4IqyKGfn5+Vq7du1hb9V4/vnnbcoIAAAAANCY4uPjj6u9At955x299NJLxmO//e1vNXDgwOrvN27cqHnz5iktLU1Dhgw5bLOMpiSsihldunQ5oj0n5s6da0M2zqObif1m7tlhjN/ZKsnynEiLbiZBi92tFTR3y7HqciLV0c3E1cwcr2PXa5fbvBu33Ba7XlvlK8ntqzTGY/ab24bEmBt91MlTecAYd/lKzfE6uoME/eYuNn7/dou4df/ugN/8ZMoC5nx9VXuM8crSrZZzSNJDeabb4EZbjJ5T57UAAAgV035ouD3wJrXrYYz7POZ/9yO91q/rgjK/7glWml8jeiJbGuNuj7njmyT5/fssJjd3UYuMTDHGJam8/BNjPMKis53rZ53w/rP7I32yezndTGx27rnn6oUXXtDzzz8vl8ulqKgoXXHFFRo4cKB69+5dY2y3bt103XXXad26dZowYYLS0tL04IMPyu22fu/QVIRVMaN79+767rvvDtsz94svvlB6erpNWTmHbiYAAAAAGssv25yjX7Y5R/sr92vqumyn0zluxMbG6rnnntPmzZvVpUsXPfLII0pKMhfZPB6PTj75ZJ188skaOXKkFi9erDFjxignJ6dB95IMRWFVrrn//vt17733atWqVXWOe/TRR23KCAAAAACAhvXYY49pyJAh+vvf/25ZyDC55JJLdPXVV2v69OmNmF1oCKuVGc8884y6deumSy65RB07djS2Zi0uLtaXX37pUIYAAAAA0DS5gy65gzbeZmLjXKGkrKxM7733nt56662jOn/AgAF6/vnnVVpaqtjY2MOfEKbCqpgxY8YM7dmzR8FgUDt37rRsmdrUl9MAAAAAAJqmf/3rXxoyZMgxXePCCy/Uv/71Lw0dOrSBsgo9YVXMSEpKUp8+fXTHHXcoIsKcemFhoUaPttr8DgAAAACA0LVx40adeuqpx3SN5ORkrV27toEyCk1hVcxo06aNbr/9dvXv37/Ocb169bIpIwAAAAA4PrgCLrls7DBi51yhJCoqSrt37z6ma+zevVuRkRbdE5uIsCpmPP3002rbtu1hxz333HM2ZAP8KOgrtzzmc5tbgUZFxJmvFagyx4PmeF2sWrbW1U5VQXM71aDP3DrUZdU+TJK7tMh8rWCZRdz85+j3/WA5R5VFC9QqixZpvsqd1tcq32U+YNEKLRCosLyWJD286/s6jwMAgMYXtHjdY/U22VdlbtkqSdHedsZ4VVWeMR4R2doYLyvJtZzDG3OiMR60eD1SeuBjy2tFRnexOOc/FnOfYIz7q4ot50DD69SpkxYvXqyRI0ce9TWWLVumwYMHN2BWoSesupl069ZNzZs3V0VFhZYtW6a33367+tjatWtVUHDwF0/Pnj2dStFW6enpSktL09y5c51OBQAAAEATs2Lvaj26cb6e3vKS06kcV84//3y9+eab2rnT+kO4uuzYsUNvvPGGLrjgggbOLLSE1coMSVqwYIEmTpyogoICde7cWVu2bJEklZeX69prr9WwYcN0zTXXOJylPXJzcxUfH+90GgAAAACaoLNanaazWp2mwspizfj6KbkDkruOBb4Nzc65QklCQoL69++va665RosXL67X7SJVVVUaN26cLrrooib/XjGsVmYsWbJE119/vTIyMjR79my1bv3jsq2+ffvqjTfe0Oeff67333/fwSwBAAAAADh6U6dO1UcffaRBgwZp27ZtR3TODz/8oEsuuUQfffSRpk6d2sgZOi+sVmY8/PDDysnJqb536M0336w1Ztq0abryyisPu0koAAAAACD0FRXV3IfN6/XK6/U6lI09unTpovnz52vUqFE66aSTNGLECP3mN79R79691b59e8XGxqqsrEw7d+7UunXrtGjRIr300ksqLS3VwoUL1alTJ6efQqMLq2JGSUnJYTdBiYqKUlmZeXNBAAAAAMDRcQVccjvQzSQlJaVGPCsrS9nZ2bbl4ZQrrrhCfr9f48eP13PPPVdno4tgMKhmzZrphRde0GWXXWZjls4Jq2JGixYtjmjc3r17GzkToKbZheauHZJ0Z2IrY9yqO4lV3LIzSR3n+APmTh91dUbx+czL2AIB8y7Wfr+5W4tk3TmkqtzcaaSyzLprycw95nMAAADq4veZX8MEAtbd6GbuNr8muSvJ/Gl3RFRLY9znMXeDi4hKspy7svwbiyPmHQIivfX/BD4quqsxXt/Xp3bJy8ursf9DU1+V8VOjRo1Senq6/vCHP+i1115TMBisNcblcum3v/2tpkyZou7duzuQpTPCqpgRCAS0efPm6h+Q6Qf5xBNPqHPnznanBgAAAABoBPHx8U1+M8u69OjRQy+//LJ27typd955Rxs2bFBBQYESExN14oknauDAgWrbtq3TadourIoZmZmZuuCCCzRjxgwNHTpULtfBZUf79u3TihUr9Oyzz2rRokVavny5w5kCAAAAQNPiCkouGzuMuGp/dn1ca9euncaOHet0GiEjrIoZw4cP14cffqhRo0bJ7XbL7XYrJiZGlZWV1WNmzZqlvn37OpglAAAAAABoTGHVmlWS5syZowULFqhr167y+XyqqKhQMBhU79699dZbb+m2225zOkXbpKenKy0tTXPnznU6FQAAAABNzIq9n2vOphw98+1rTqcC1BJWKzMOGTNmjMaMGaOtW7cqPz9fHTp0UMeOHZ1Oy3a5ubnH9b1jAAAAQGOYnHK60ymEhLNanaKzWp2iwspiPbxxvtPpADWEfDHD7/dr9+7d8vv96tChQ41jqampSk1NdSYx4BgFfKXmuCfGGPf7Cqyv5TpgjFdVmnfi9ldZX8tXaT5WeeA7y3Pq6uYCAADgpKDf3IkjNuFkY9xXx+skb7NUY9zqdV158XpjPDLGuhtc0F9hjsu8WYU3upvltay6zll1ybN6ven3lUiS3Da3ZrVzLoSfkCxm9OvXT7t27VJhYaFSUlLUtWtXDRo0SFdffbXTqQEAAAAAAIeFZDHj//2//6ff/OY3ysnJ4TYKAAAAAABQQ0gWMyIiIvTXv/6VQgYAAAAAhAhXwObWrDbOhfATkt1MTjzxRLVu3fqoz7/yyisbMBsAAAAAABBKQrKYkZiYeEznf/XVVw2UCQAAAAAACDUheZvJsdq7d6/TKQAAAABAk+IOHHzYOR9gJSSLGV9++aXeeecdxcbG1vvcbdu2aceOHY2QFXB0Zu0zF9fubJVkjAcD5vZhVWXmNquS5K8qrvf8AAAAx4PZ+83tRu9K6mR5zsO7vjfG/9DpTGPc5Ta3OnV7zO9nggGf5dxWvLHdjfGK8m+U/fVLxmN/7n2HMR7wm187WrVyDQTMrWcBJ4VkMWPv3r369a9/7XQaIS89PV0ej0eZmZnKzMx0Oh0AAAAATUju/m+VW/it/AGWSCD0hGQxQ5KCweBRn+tyuRowk9CVm5tLxxcAAAAAjSK9RRelt+iioqpSPbp1qdwBl9wB+95r2TkXwk9IFjNSUlL0r3/9S16vt97nbtu2Teedd14jZAUAAAAAAEJBSBYzunTpop49ex7VuZ07d1aHDh0aOCMAAAAAgBOKiopqfO/1eo/qg280LSFZzDhWLVq0cDoFAAAAALDFHzufbYzHJpjjR8sVOPiwy6G5UlJSasSzsrKUnZ1tXyIISSFZzNi9e/cxnX/uuec2UCZA45m1d5fTKQAAABy3XC5zB5K6PPT9p/Uaf28X8/uSqjLr7osR3tbGeEXJ18a4VccUSSrY/qIx3qxlH2M80ptijEfogOUcdsjLy6uxVyCrMiCFaDFj48aNysvLq1WBO1KPPfZYA2cEAAAAAHBCfHw8jQ9Qi9vpBEz8fr+uuOIKbdiwwelUAAAAAACS3AH7H4CVkFyZMXr0aO3Zs0dXXXWV3G63EhISNHjwYN1+++1OpwYAAAAAABwWksWMnJycWrGqqioHMgEAAAAAAKEmJIsZJpGR9d+gBwAAAADqa2Lr9sa4y21+TxIV06Fe48OVU91MAJOwKWYAAAAAQEPxVewzxt0R0Q02R1W5Rfc6l0fTd5j3B/xjaoYxXmnRASW6eTfL+WMSfmGMuz3NjfGqyh+McV9VqeUcgFNCcgNQAAAAAAAAK6zMCGPp6enyeDzKzMxUZmam0+kAAAAAaEJy93+r3MJv5Q8cvN/DHQjKHQjaNr+dcyH8UMwIY7m5ufRbBgAAANAo0lt0UXqLLiqqKtWjW5c6nQ5QA7eZAAAAAACAsMLKDAAAAACw2Z2tkozx6LgeNmcChCeKGQAAAACOO7P3F9T7nDsTWxnj/qpCc9xXbIwHA1WWc1SVmTuKREQl1mvuuuYJ+M3dScoKvjDGD/xvzwxX0ObWrGyZgTpwmwkAAAAAAAgrFDMAAAAAAEBY4TYTAAAAAMBh0ZoVoYSVGQAAAAAAIKxQzAAAAAAAAGGF20wAAAAAAIcXsLebieycC2GHYgYAAAAAHIFZ+/bWa/zklNONcbcnxvIcX4X1HNN3bKgV+1O3gZbjy4rXG+NRsSnGeER0G2Pc4/NJKrGcB3ACt5kAAAAAAICwwsoMAAAAAMBhuQJBuWzsMGLnXAg/FDMAAAAAACGrqKioxvder1der9ehbBAquM0EAAAAABCyUlJSlJCQUP2YOnWq0ykhBLAyAwAAAECTdU+Hky2PTfvhCxszCX9O3WaSl5en+Pj46jirMiBRzAhr6enp8ng8yszMVGZmptPpAAAAAPiJqXmrGn0Ov6/Q8pjbE2uMR0QmHtG1Pyveof8W75Q/4GyP1Pj4+BrFDECimBHWcnNz+Z8aAAAAQKPoE9defeLaq8hXocd/+K/T6QA1UMwAAAAAAByWOxCU28bbTOycC+GHDUABAAAAAEBYoZgBAAAAAADCCreZAAAAAECImNSuh+Wx6Ts32ZhJba5AQC4bNwO1cy6EH4oZAAAAAJqspt5+taxwQ53HZ+8vqBWbnHK6cazbE9MgOQF24DYTAAAAAAAQVliZAQAAAAA4LFcwKFfQvg4jds6F8MPKDAAAAAAAEFYoZgAAAAAAgLBCMQMAAAAAAIQV9swAAAAAABvdldTJ8pjb7W3QuSa26VArFhnd9qiu5QoEbW7Nyp4ZsEYxAwAAAABChDsitl7jo2LbWx7zVRUa497Y7sa4y+UxxiOrSiX9t155AY2N20wAAAAAAEBYYWUGAAAAAOCwDt5mYmNrVm4zQR1YmQEAAAAAAMIKxQwAAAAAABBWuM0EAAAAAHAEAlLQvm4mkp1zIdxQzAAAAAAAGz2863vLY3e3Ta3XtaZtX295zKoFbFXFD+Z4mTle6vPVKyfADtxmAgAAAAAAwgorMwAAAAAAh+UKBOQK2Hfrh51zIfywMgMAAAAAAIQVihlhLD09XWlpaZo7d67TqQAAAABoYlaVFuvZvdu1sGCn06kAtXCbSRjLzc1VfHy802kAAAAAMLgjwfxaPTKm7VFdz2pDz7o2FD0Wp8fG6fTYOBX5fHpy3w9yBQNy2djNxM65EH4oZgAAAABAiPD7So1xtye63teyKnJMbN3e8pyZe3bUim3btk1PpqTUe36gMXGbCQAAAAAACCuszAAAAAAAHJYr6Jcr6Ld1PsAKKzMAAAAAAEBYoZgBAAAAAADCCreZAAAAAAAOLxA4+LBzPsACxQwAAAAAaILuSU6r9zl/6jawVqzIV94Q6Ry1oqKiGt97vV55vV6HskGooJgBAAAAAI1gdmGRMX5nq6Q6z5u1d1djpFPNE9HcGPc269ao8x6tlJ+1hc3KylJ2drYzySBkUMwAAAAAAByWKxiQK2jfrR+H5srLy1N8fHx1nFUZkChmAAAAAABCWHx8fI1iBiDRzQQAAAAAAIQZihkAAAAAACCscJsJAAAAAIQQqw1CG3tj0MMKBqSg3975AAsUMwAAAADARnUVJQ7X6aQ+pm1fb4xPatfDGK8o+UbTd26qFd+2bZvm/KyjCOA0bjMBAAAAAABhhZUZAAAAAIDDCwbsvfWD20xQB1ZmAAAAAACAsEIxAwAAAAAAhBVuMwEAAAAAG92REG95zBURbWMmtU3ueEqtWJGvUpLkCvrlsrGbiZ1zIfxQzAAAAACAEBH0lWt2YVGjzuGOjLOYnOIBwge3mQAAAAAAgLDCygwAAAAAwOEF/fau3mClCOrAygwAAAAAABBWKGYAAAAAAICwwm0mAAAAABBCrLqdNPbGoIcVDBx82DkfYIFiBgAAAADYqK6iRF1tWxtKZek2y2Oz9u6qFdu2bZvmpaQ0ZkpAvXGbSQgIBAJ66KGHlJqaqri4OF100UX65ptvnE4LAAAAAICQRDEjBEydOlUtW7bUp59+qg8++EDff/+9rr76aqfTAgAAAICf8DvwAMwoZjisrKxMJ554oiZMmKC2bduqT58+yszM1Nq1a51ODQAAAACAkEQx4yh88MEHWrRoUYNcKyYmRpdddlmNWEFBgX796183yPUBAAAAAGhqKGbUw6pVqzRgwABdcMEFWrVqVZ1jv/nmG40ZM0annXaafvnLX+r000/XU089ddg5Nm3apLVr1+qZZ55pqLQBAAAA4Ji5gn7bH4AVupkcgf379+upp57Srl27tGLFisOOX7Nmjc477zwNGTJEK1euVGRkpD7++GMNGDBAubm5evrpp2uds3v3bk2fPl1/+ctfVFpaqpSUFD366KON8XQAAAAANBGTU06v8/jUvLo/hAXCFcWMI5CQkKC7775b0sG2RAsXLrQcW1hYqCFDhigiIkLz5s1TZGSkJCkjI0N33XWXsrOzlZGRobFjx9Y4r02bNpoyZYpGjRqlrKwszZkzRxdeeKEuueSSRnteAAAAAEJLXW1bG0pEVGKjzwE0Nm4zOQIul6v66+jo6DrH/uUvf1FeXp5GjBihuLi4GsfGjRsnSbr33ntVVVVV69zo6Gidfvrpev3119W5c2f9+9//boDsAQAAAKABBANS0G/jI+D0M0YIo5jRwA7tddGvX79ax1JSUtS1a1f98MMPWrp0qeU1IiIi1Ldv3+pVHQAAAAAA4EcUMxrQ7t27tXnzZklSWlqaccwvfvELSdJ7771X57V++OEHXXjhhQ2bIAAAAAAATQDFjAb05ZdfVn/doUMH45jk5GRJ0rp16yQdbMN6xx136NNPP60e89Zbbyk1NVUXXXRRI2YLAAAAAEB4YgPQBrRv377qr3++X8Yh8fHxkqRdu3ZJkiorK/XBBx9o7ty5GjBggPr27atOnTrpueeeO+x8RUU1Nwfyer3yer1Hmz4AAACAJuaupE61Ym5PbJ3nVFRUqKKiovr74uLi/33l/9/DLrRmhTWKGQ2otLS0+uuoqCjjmEPFhpKSEklS27ZttXr16qOaLyUlpcb3WVlZys7OPqprAQAAAAg//sp9lseCFhto+qsK67zm1KlTdf/99x9TXkBjo5jRgGJiYqq/rqqqMhY0DnUxiY2tuxp6JPLy8qpXekhiVQYAAACAYzZ58mTdcccd1d//8MMPlnsCAk6hmNGA2rVrV/31gQMHlJhYu3/zgQMHJB1ckXGs4uPjaxQzAAAAAOBY/fz29erb2w+1TLWLnXMh7LABaAM66aST5HK5JEk7duwwjjkU79Wrl215AQAAAADQlFDMaECJiYnq3bu3JGn9+vXGMYc6npx//vm25QUAAAAAQFPCbSYN7Morr9Tnn3+ujz76SJdffnmNY/n5+dq4caMSExN18cUXO5QhAAAAgFA1qf2Jlsem79jQ6PPf3Ta1VqzI5/vfV3QzQeigmFFPZWVlkqRAwLwz8Pjx4zV79mwtXLhQ06ZNq7HR5/z58xUIBHTnnXfW2Cz0aKWnp8vj8SgzM1OZmZnHfD0AAAAA4WVG/tZ6n3NnYitj3BMZV+P7VaXFWlNWLL/Fex/ASdxmUg979+7V8uXLJUkff/yxfNUVyh8lJiYqJydHxcXFuvnmm6u7l+Tm5mrq1KkaPHiwJk2a1CD55Obmav369RQyAAAAADS402PjNK5Vsq5o2e7wgwGbUcw4An6/X3369FHnzp21fft2SdK///1vJScna9iwYbXGDxgwQCtWrFBJSYnOOussnXPOObrxxhv1wAMPaNGiRfJ4PHY/BQAAAAA4JsGg3/YHYIXbTI6Ax+PRZ599Vq9zevfurYULFzZSRgAAAAAAHL9YmQEAAAAAAMIKKzMAAAAAAEeAbiYIHRQzAAAAACBM3dGiZZ3HZ+8vsCkTwF4UM8IYrVkBAACApmX6jg2NPsesfXuN8ckdT6nx/X8P7NGqkr20ZkVIopgRxnJzcxUfH+90GgAAAACaoDOat9YZzVuryFepefkbpGBQCtpY2AgG7ZvrGL3++ut655139OSTT1bH5s2bp7Vr16pjx476/PPPNWnSJPXp08fBLJsWihkAAAAAABylPXv26KabbtKgQYOqY3/729/0zDPPaOXKlYqIiFBubq4uuugiff3110pKSnIw26aDbiYAAAAAABylRx99VKmpqTViU6ZM0ZgxYxQRcXD9QHp6ujp16qQFCxbYn2ATRTEDAAAAAHAE/A48QtvLL7+sQYMGKSYmpjr2zTffaPPmzfrFL35RY+zJJ5+sd9991+4UmyyKGQAAAAAA1FN+fr42b96sjIyMGvFNmzZJUq39DePi4rRhQ+Nv8Hq8YM8MAAAAAAhT3rguDXatqds+N8a3bdumeSkpDTaP3TIzM5WcnKw//vGPlmO++uor3Xfffdq0aZMSExN14MABXX/99Ro/frxcLpfxnEceeUTZ2dm14vv27ZMkNWvWrEa8WbNmKiigVW5DYWUGAAAAAOAI+KWgjY9juM3E5/Np6dKl6tevn+bNm6eqqirLsV9//bXOPvtsJSYmavXq1frggw/0t7/9Tffcc49uvPFG4zkvvviiLr30UkVHR9c6dmifDI/HUysnn8931M8JNVHMCGPp6elKS0vT3LlznU4FAAAAQBMzd+5cpaWlqV+/fk6nUi8333yz+vXrpxdffFHNmzevc2wwGNT1118vv9+v2bNny+0++BY5LS1Nt99+u55++mktWbKkxjnbt2/Xd999p7POOst4zTZt2kiSAoGabWxLSkrUokWLo31a+BmKGWEsNzdX69evV2ZmptOpAAAAAGhiMjMztX79ei1btszpVOrl8ccf1/Lly/XXv/5Vffr0qXPsihUrtHz5cg0bNqxW4eOqq66SJM2ePbtGfMmSJVq7dq3Gjh1b/Vi7dq2WL1+usWPHVhdE8vPza5y3b98+9ejR41ifHv6HPTMAAAAAAMel1157TdLBVe8/17VrVyUmJmrZsmUqKiqq3tDzhhtu0A033FBjbGpqqn71q19Vt1496aSTtGHDBp1//vnVY1atWqWxY8c2zhM5DlHMAAAAAIAQd2erJGM8KrZjnefdk5xmjE/bvr5WbHLK6caxRb7K/31ld7vUxp9rzZo1kqTOnTsbj3fq1En79u3T6tWrdd5551lex+/3y+//Md8bbrhB8+fP10033SSXy6V3331XRUVFuummmxr2CRzHKGYAAAAAQJiqLN2mWXt3GY9ZFTLwo82bN0uSEv9/e3cfJ3O9/3/8ubsYZGe1aKPGIl0NRzhWhGyiTRxFDuJISrJt36SoOIoS24VycRq5qKTS1bFtKTl1qtWPOq5CyFFHrnaRazvWxWLn8/vDmTnG7tplZz/z+ew+7rfb52bmM595v18fr9tezGvfF7GxBb7uP79r164CX3/rrbf01ltvKSsrS2lpaTp48KDGjBmjoUOHyuv16r777lP9+vW1atUq/fOf/1RcXFzp3Eg5RDEDAAAAAGBZXq836LnD4ZDD4QhJ2/5tVAvaleTM84VtqepfM6MgY8aMKXmAKBQLgAIAAAAAimQYeaYfkuRyuRQTExM4UlNTQ3ZpfjFNAAAgAElEQVRPR44ckfS/7VTP5t9e9dChQyHrE6HByAwAAAAAgGVlZmYGFt+UFLJRGZIUGRkZtNbF2Y4ePRq4DtZCMQMAAAAAYFlOpzOomBFK0dHRgakmBfEXM87ethXhRzHDxhISEhQVFaWUlBSlpKSEOxwAAAAAZciPOXu1Kmef8ny+/54pe7uZ+HcrOXHiRIGv+6ehXH75uXeNgfkoZtjYihUrSq1CCQAAAMB8w2vWPu/3PBZbo8DzFSuXbOeMP1arpT9WqyXvqROa9vvPJWrLqho1aqQ1a9Zoz56Cd4TZvXu3JOnqq682MywUA8UMAAAAALC4iIgoTdyXf3vQwgoZknTy+G69fGB/sftIzVxV4PmsrCxNc7mK3Y6d3HjjjZo7d67Wr1+vrl27Br22Z88e7d69W3Xq1KGYYUGsYgIAAAAAKIa8MBylq3v37qpQoYI+//zzfK+lp6dLknr37q2IiIhSjwXnh2IGAAAAAKDMMQxDWVlZkqTt27fLMIx819SqVUvJycn6/vvvtWjRosD5vLw8TZs2TXFxcRo9erRZIeM8MM0EAAAAAGBZXq836LnD4Tjn9qypqalKT0/Xli1btG/fPknS7NmzNX/+fNWvX189evTQyJEjA9e/8MIL2rlzp7p3765x48YpPj5e06dP165du5SWlqbY2NjSuTGUCMUMAAAAAECRDMMnwzBvNxPDOL2Liuus9TrGjBmjsWPHFvq+kSNHBhUrilKlShXNmzdPa9as0Q8//KD169erf//++uCDDxQdHX1BsaP0UcwAAAAAgHJkpKt5oa8VtghoOGVmZgbt4niuURkl0bRpUzVt2rRU2kboUcwAAAAAAIsoaMeSczmf3Ursyul0BhUzAIliBgAAAACgWMzZYSS4P6Bg7GZiYwkJCXK73fJ4POEOBQAAAEAZ4/F45Ha7lZiYGO5QgHwYmWFjK1asYLgVAAAAgFKRkpKilJQUZWVl5VuEEwg3ihkAAAAAgKIZvtOHmf0BhWCaCQAAAAAAsBVGZgAAAABAOWLF7VeB80UxAwAAAABQJEN5MkzcYcTfl9frDTrvcDjkcDhMiwPWxDQTAAAAAIBluVwuxcTEBI7U1NRwhwQLYGQGAAAAAMCyMjMzg3ZxZFQGJIoZAAAAAAALczqdQcUMQKKYAQAAAAAolrz/Hmb2BxSMNTMAAAAAAICtUMwAAAAAAAC2wjQTAAAAAEDRjLzTh5n9AYVgZIaNJSQkyO12y+PxhDsUAAAAAGWMx+OR2+1WYmJiuEMB8mFkho2tWLGCVX0BAAAAlIqUlBSlpKQoKytLLpcr3OEAQShmAAAAAACKZMgnw8QdRgz5TOsL9kMxAwAAAABgWV6vN+i5w+GQw+EIUzSwCtbMAAAAAABYlsvlUkxMTOBITU0Nd0iwAEZmAAAAAACKwfffw8z+pMzMzKC1AhmVAYliBgAAAADAwpxOJxsfIB+mmQAAAAAAAFthZAYAAAAAoGhG3unDzP6AQjAyAwAAAAAA2ArFDAAAAAAAYCtMMwEAAAAAFMlQngyZN/XDzL5gP4zMAAAAAAAAtkIxAwAAAAAA2ArTTAAAAAAARTN8pw8z+wMKwcgMAAAAAABgK4zMAAAAAABYltfrDXrucDjkcDjCFA2sgpEZNpaQkCC32y2PxxPuUAAAAACUMR6PR263W4mJiZL+t5uJmYckuVwuxcTEBI7U1NQw/q/AKiIMwzDCHQTOj9frVUxMjLKzs+V0OsMdDgAAAIAyLCsrSy6XS2P+8FdVjqpsWr/H847rmXXjlZmZGfS5h5EZkJhmAgAAAACwMKfTyR9xkQ/TTAAAAAAAgK0wMgMAAAAAUDQj7/RhZn9AIRiZAQAAAAAAbIViBgAAAAAAsBWmmQAAAAAAisGQ5DO5P6BgjMwAAAAAAAC2QjEDAAAAAADYCtNMAAAAAABFMpQnQ+btMGJmX7AfihkAAAAAYHEjLqlb6Gsv7dluYiSANTDNBAAAAAAA2AojMwAAAAAARTPyTh9m9gcUgmIGAAAAAMCyvF5v0HOHwyGHwxGmaGAVTDMBAAAAAFiWy+VSTExM4EhNTQ13SLAARmYAAAAAAIpkyCdDPlP7k6TMzEw5nc7AeUZlQKKYAQAAAAC29nhcvUJfe3H3VtPiKC1OpzOomAFIFDMAAAAAwPIK2371XIUMoCyjmAEAAAAAKBq7mcBCWAAUAAAAAADYCsUMAAAAAABgK0wzAQAAAAAUKVy7mQAFYWSGjSUkJMjtdsvj8YQ7FAAAAABljMfjkdvtVmJiYrhDAfJhZIaNrVixgi2KAAAAAJSKlJQUpaSkKCsrSy6XK9zhAEEoZgAAAACATeWdPGxeZ+xmAgthmgkAAAAAALAVihkAAAAAAMBWKGYAAAAAAABbYc0MAAAAAEAxGJKp26UaJvYFu2FkBgAAAAAAsBVGZgAAAACATb18YH+4QwDCgmIGAAAAAKBIhpEnw8TtUs3sC/ZDMQMAAAAAYFlerzfoucPhkMPhCFM0sArWzAAAAAAAWJbL5VJMTEzgSE1NDXdIsABGZgAAAAAAisEnc3czOd1XZmamnE5n4CyjMiBRzAAAAAAAWJjT6QwqZgAS00wAAAAAAIDNMDIDAAAAAFAkdjOBlTAyAwAAAAAA2ArFDAAAAAAAYCtMMwEAAAAAFM3wnT7M7A8oBCMzAAAAAACArVDMAAAAAAAAtsI0EwAAAABAkQzlyZCJu5mY2Bfsh5EZAAAAAADAVihmAAAAAAAAW2GaCQAAAACgaOxmAgthZAYAAAAAALAVihkAAAAAAMBWmGYCAAAAACiSIZ/Ju5kwzQSFY2RGmB05ckTDhg1TfHy8qlevrltvvVUbN24Md1gAAAAAAFgWxYwwS05OltPp1KRJk5ScnKyMjAx17NhRhw4dCndoAAAAAABYEtNMwmj9+vW67rrr9Nhjj0mSevTooRo1amjEiBH66quv1KtXrzBHCAAAAACA9TAy4zxlZGRo/vz5IWnr8OHDSklJCTqXlJQkScrOzg5JHwAAAAAQEv6tWc08gEJQzCimVatWKSkpSR06dNCqVavOee1vv/2mu+++W82aNdMNN9yg5s2ba+bMmfmua926tSpXrhx07vjx44HXAAAAAABAfkwzKcKhQ4c0c+ZM7dmzR0uXLi3y+jVr1qh9+/bq1q2bli9frooVK+r7779XUlKSVqxYoVmzZp3z/RkZGercubMaN24cqlsAAAAAAKBMoZhRhJiYGD3++OOSpKysLH344YeFXpudna1u3bqpQoUKmjZtmipWrChJatOmjUaMGKGxY8eqTZs2uueeewp8/5EjR/T+++8rPT095PcBAAAAACVi5MkwTBzcb5zeBtbr9Qaddjgccjgc5sUBS2KaSREiIiICj8+eEnK26dOnKzMzUz179lR0dHTQa/fee68kafTo0Tp58mSB7x8xYoQmT56sevXqlSxoAAAAACgjXC6XYmJiAkdqamq4Q4IFMDIjhF5//XVJUmJiYr7XXC6XGjRooM2bN+vrr79W586dg16fOHGikpKS1L59ezNCBQAAAABbyMzMlNPpDDxnVAYkRmaEzN69e7Vp0yZJktvtLvAa/zoYX331VdD52bNnKy4uTrfffnvgXE5OjrZu3Vo6wQIAAADAefOF4ZCcTmfQQTEDEiMzQubnn38OPL7ssssKvKZOnTqSpPXr1wfOpaenKy0tTffcc4/mzZsn6fTaGZ988onee++9UowYAAAAAAB7opgRIgcOHAg8Pnu9DD//0Kg9e/ZIkpYsWaK77rpLubm5WrBgQdC1ycnJqlKlyjn7ZCEcAAAAAKGWm5ur3NzcwPPDhw+HMRqgYBQzQuTo0aOBx5UqVSrwGn+h4ciRI5Kktm3b6vjx4xfcp8vlCno+ZswYjR079oLbAwAAAIDU1FQ988wz+V8w8qQw7GYCFIRiRoicOYri5MmTBRY0/LuYVK1aNSR9shAOAAAAgFAbOXKkHn300cDzHTt2FLouIBAuFDNC5NJLLw08zsnJUWxsbL5rcnJyJElxcXEh6dO/AA4AAAAAhMrZ09fPnt4OWAHFjBC59tprFRERIcMwtGvXrgKLGbt27ZIkNWrUyOzwAAAAAKBEDPlk/HeHEbP6AwrD1qwhEhsbqyZNmkiSNmzYUOA1/h1PbrrpJtPiAgAAAACgrKGYEUJ9+vSRJC1evDjfa7t379avv/6q2NhY3XLLLWaHBgAAAABAmUEx4zwcO3ZMkuTzFTzcafDgwapVq5Y+/PDDoN1NJGn27Nny+Xx67LHHitxytbgSEhLkdrvl8XhC0h4AAAAA+Hk8HrndbiUmJv73jE8yTDyYZoJziDAMwwh3EHawf/9+NWnSRDt37lSHDh305ZdfqkKF/EuOfPnll+revbv69OmjGTNmqGLFilqxYoU6duyotm3bav78+YqKiipRLF6vVzExMcrOzmYBUAAAAAClKisrSy6XS080SJIjsqJp/eb6TuqFzV/yuQcFYmRGEfLy8tSiRQvFx8dr586dkqRvv/1WderUUffu3fNdn5SUpKVLl+rIkSNq1aqV2rVrpwceeEDPPvtsSAoZAAAAAACUd+xmUoSoqCitXLnyvN7TpEkTffjhh6UUEQAAAACYzzDyZBjm/T3cMPJM6wv2w8gMAAAAAABgKxQzAAAAAACArTDNBAAAAABQNMP47y4jJvYHFIKRGTbG1qwAAAAASkv+rVkB62BrVhtia1YAAAAAZvFvzfp4vZvliDRvcH+u75Re3PoNn3tQIEZmAAAAAAAAW2HNDAAAAABAkQzlyVCEqf0BhWFkBgAAAAAAsBWKGQAAAAAAwFaYZgIAAAAAKJrhM3lrVhP7gu0wMsPG2JoVAAAAQGlha1ZYGVuz2hBbswIAAAAwi39r1hHxN5q+NetL2/4fn3tQIKaZAAAAAACKZBh5MgwTdzMx2M0EhWOaCQAAAAAAsBWKGQAAAAAAwFaYZgIAAAAAKAbffw8z+wMKxsgMAAAAAABgKxQzAAAAAACArTDNBAAAAABQpHDtZuL1eoPOOxwOORwO0+KANTEyw8YSEhLkdrvl8XjCHQoAAACAMsbj8cjtdisxMTGscbhcLsXExASO1NTUsMYDa4gwDMMIdxA4P16vVzExMcrOzpbT6Qx3OAAAAADKsKysLLlcLj3mailHpHmD+3N9p/Ry5nJlZmYGfe5hZAYkppkAAAAAAIrD8J0+zOxPktPp5I+4yIdpJgAAAAAAwFYoZgAAAAAAAFthmgkAAAAAoEiGYcgwcZoJyzviXBiZAQAAAAAAbIViBgAAAAAAsBWmmQAAAAAAimbkSWbO/DDyTOwMdsPIDAAAAAAAYCsUMwAAAAAAgK1QzLCxhIQEud1ueTyecIcCAAAAoIzxeDxyu91KTEwMdyhAPhEG+93YjtfrVUxMjLKzs+V0OsMdDgAAAIAyLCsrSy6XS8PqNJEjMsq0fnN9eZq0cy2fe1AgRmYAAAAAAABboZgBAAAAAABsha1ZAQAAAABFY2tWWAgjMwAAAAAAgK1QzAAAAAAAALbCNBMAAAAAQJEMwyfDiDC1P6AwjMwAAAAAAAC2QjEDAAAAAADYCtNMAAAAAABFYzcTWAgjMwAAAAAAgK1QzAAAAAAAALZCMcPGEhIS5Ha75fF4wh0KAAAAgDLG4/HI7XYrMTFRkn83E3MPoDARhmGYOesJIeD1ehUTE6Ps7Gw5nc5whwMAAACgDMvKypLL5dLQSxrKERllWr+5vjxN2bOJzz0oECMzAAAAAACArbCbCQAAAACgaIbP5N1MmGaCwjEyAwAAAAAA2ArFDAAAAAAAYCtMMwEAAAAAFOn0DiPm9gcUhpEZAAAAAADAVihmAAAAAAAAW2GaCQAAAACgGPJk7nYmTDNB4RiZAQAAAAAAbIViBgAAAAAAsBWKGQAAAAAAwFZYMwMAAAAAUCS2ZoWVMDIDAAAAAADYCsUMAAAAAABgKxQzbCwhIUFut1sejyfcoQAAAAAoYzwej9xutxITEyVJhpFn+gEUJsIwzJz1hFDwer2KiYlRdna2nE5nuMMBAAAAUIZlZWXJ5XLpodhL5Yg07+/huT6fXj3wO597UCBGZgAAAAAAAFthNxMAAAAAQNEMn2TmuH52M8E5MDIDAAAAAADYCsUMAAAAAABgK0wzAQAAAAAUyVCezNw/whDTTFA4RmYAAAAAAABboZgBAAAAAABshWkmAAAAAICiGT6Zup2JiVNaYD+MzAAAAAAAALZCMQMAAAAAANgK00wAAAAAAEUyfD4ZERHm9cc0E5wDIzPKgNzcXI0dO1a5ubnhDgX/RU6sibxYDzmxHnJiPeTEmsiL9ZAToHyJMCh32Y7X61VMTIyys7PldDrzPUf4kRNrIi/WQ06sh5xYDzmxJvJiPeSk9GRlZcnlculBZ4wcJo7MyDUMTfNmk1MUiJEZKDGPx2Prfi+knXDdc3GRE+sJZ3yh6JucWLPv823H6jmR+P5lReTEevj+ZU3l4mvFyDP/AArByAwbstrIDLfbrQ0bNti23wtpp6j3kBPz2ynOe8KZl3DlJFR9kxNr9n2+7Vg9JxLfvwpCTsxvx+pfK3z/KhhfK6XXTmBkRnQ180dmHM5hZAYKxAKgNuSvP3m93gL/NVteXl5Y+g5VvxfSTlHvISfmt1Oc94QzL+HKSaj6JifW7Pt827F6TiS+fxWEnJjfjtW/Vvj+VTC+VkqvncOHD0uSTpj8d3Cz+4O9MDLDhvyVUQAAAAAoyy699FJt2bJFlStXDncosBiKGTbk8/m0c+dORUdHK8LEYV4AAKBsyMnJ0erVq3Xy5El16NAh3OHgDIZhaPPmzbriiivCHQoQYBiG9u7dq5o1ayoy0txlFytVqkQhAwWimFGO+Hw+Pf/885o5c6b279+vVq1aafr06fl+WG7evFkTJkzQHXfcoa5du4YpWgAAUFxHjhzR6NGj9fHHHys7O1utWrXS5MmTdc011wRdt2XLFo0cOVI+n09Dhw5V69atTf9gUp4UJy8bNmxQ48aNdeav5HPnzlXfvn3DEXKZV9yvFb+BAwcqPj5eY8eONTdQAEXip1c5kpqaqosvvljLli1TRkaGtm/frv79+wde9/l8Sk9P10MPPaQ33nhDx48fD2O0ZZ/P59OECRNUr149RUdHq1OnTvrtt9+Crnnrrbf0yCOP6Pnnn1evXr20bt26MEVbfhw5ckTDhg1TfHy8qlevrltvvVUbN24MvJ6bm6snnnhCTz31lIYNG6aePXtq+/btYYy47CsqJz6fTy+//LKefvppPfvss+rRo4c++uijMEZc9hWVk1atWikiIiLf0bJlyzBGXbYlJyfL6XRq0qRJSk5OVkZGhjp27KhDhw4FrlmyZIlatGihG2+8UR999JHatGlDIaOUFScvkyZNksfj0ezZszV79mx98MEH6tWrVxijLtuKkxO/zz77TG+//XYYogRQHCwAWk4cO3ZM11xzje68805JUlxcnFJSUjRq1KjANZGRkerevbucTqcWLlwYrlDLjdTUVMXGxmrZsmXKzMxUv3791L9/f/3www+SpFdffVVz5szR0qVLFRUVpR9//FEdO3bU6tWrVadOnTBHX3YlJyerfv36mjRpklasWKFXXnlFHTt21Pr161W9enUNGDBAbrdbTz/9tCRp6tSpat++vX766SdW2S4lReVk3LhxOnDggKZMmSJJOnjwoBo2bKj4+Hhdf/31YY6+bDpXTrZv366aNWtq3rx5cjgcgfdMnTpVN954YxijLrvWr1+v6667To899pgkqUePHqpRo4ZGjBihr776Sr169dLmzZv1pz/9SYMGDdKDDz4Y5ojLh+Lk5ffff1d2draSk5PDHG35UJyc+O3fv18LFy5knTrAygyUW2PHjjV69eqV73xGRoYhyfj73/8ehqjKh6NHjxrz5s0LOjdlyhTjoosuMgzDMA4dOmRER0cbM2bMCLqmZcuWxoMPPmhanOXNunXrjIkTJwade+mllwxJxocffmisXr3akGQsW7Ys8Prhw4eNChUqGE8//bTZ4ZYLReXEMAzjuuuuM6ZMmRJ0TcuWLfO9D6FRVE4+/fRT4/jx40Gv+3w+4/LLLzfWrVtnZqjlxg8//GAcO3Ys6NzatWsNScbMmTMNwzCMnj17GlWrVjUOHjwYjhDLpeLkZdSoUYYkIz4+3khOTjY2btwYjlDLjeLkxC8lJcXYvXu3ER8fb4wZM8bEKAEUF2MLLSAjI0Pz5883tc///Oc/Wrt2rV5//XVT+7UDM/JRpUqVwCgZv4MHD6pLly6SpIULF+rw4cNq1apV0DWtW7fWBx98EDSvtrwwIy+HDx9WSkpK0LmkpCRJUnZ2tr777jtJ0mWXXRZ4vVq1arriiiuUnp5eqrFZkRVyIkk1a9bUnDlzlJubK+n0FIhff/1Vf/zjH0s1NiuyQk66desWNCJDkpYvX67KlSurcePGpRqbFZmRk9atW+dbHM8/VbR169Y6dOiQ0tPTVb9+fT3xxBNq0aKFatasqSFDhujYsWOlGptVWSEvknTNNdfo0UcfVfXq1fXaa6+padOmpv9OaBVWyYkkffTRR7rxxht1ySWXlGo8AEqGYkYYrVq1SklJSerQoYNWrVp1zmt/++033X333WrWrJluuOEGNW/eXDNnzjzvPvfu3avhw4erWbNmSk9P11NPPXWh4Zc54ciH39nFpbVr10qS6tatG3Rd3bp1deDAAW3ZsuWC+7IbM/NS1C85J0+elPS/D9F+NWrUyLfeSVlmpZxI0pgxY7Ru3Tp16tRJq1ev1sCBA/XCCy8oMTHx/G7MxqyWk7OlpaWpe/fuxe6jLAjnzxTp9AfDzp07q3Hjxlq+fLny8vLUuHFj/e1vf9PKlSuVnp6ut956SyNGjChRP3ZjpbxIUv/+/fXyyy9rzZo1+vTTT1W1alX17dtX+/btK1E/dmK1nPz+++9atGgR65YANsCaGWFw6NAhzZw5U3v27NHSpUuLvH7NmjVq3769unXrpuXLl6tixYr6/vvvlZSUpBUrVmjWrFnF7rtWrVp67rnn1LdvX40ZM0ZTpkxRx44dy/WuJeHMx969e/XCCy9o+vTpOnr0qFwulyZPnqz9+/dLkqKjo4Ou9z/fs2ePGjRocB53aT/hzMuZzvwlZ+vWrZKkf/3rX3K73YFroqKiLqhtu7FiTiSpXbt2Sk9PV69evdS8eXONGjVKgwcPvqC27caqOTnbxx9/rHfeeeeC2rYbK+TkyJEjev/99wMjxvwfjO+//35VqlRJ0umvmzvvvFMzZszQyy+/nG80TVljxbycrVu3bkpLS9NNN92k+fPn69577z3vPuzEqjkZPXq0nn/++fNuC0AYhHueS3nk8/kCj3v37m1IKnQu3qFDhwyXy2XExsYaXq836LWxY8cakozZs2cHzr377rvGRRddFDjcbnehcZw8edKIj483hg0bFnS+vK2ZUZr5KI5jx44ZP/74o9G1a1dDkvHZZ58ZDz30kCHJyMvLC7r29ddfNyQZq1atOq8+7CjceTEMw8jJyTGaNm1qbNmyxTAMwzh16pTRqFEjo0GDBsb27dsNwzCMJUuWGE6n02jYsOF5t283VsyJ3+TJk40hQ4YYjRs3NiQZffr0MU6dOnXe7duNlXPit3r1aqN27dpBsZZlVvgZn5ycbCxatCjwfOHChYYkIyMjI+i6F1980ZBULtZpsGJeCnPFFVcYEyZMKPqmbM6KOZk9e7aRlpYWdA1rZgDWxTSTMIiIiAg8Pnuo7tmmT5+uzMxM9ezZM99f6f0V+9GjRweGv3fr1k1r1qwJHF988UWhbVeoUEEtW7ZUxYoVL/RWyoTSzMfcuXNVrVq1wNGoUaN8bVauXFnNmzdXenq64uPj9e2336p27dqSTs9NP1NOTo6k4DUbyqpw50WSRowYocmTJ6tevXqSTo/AyMjIUPv27dW1a1f169dP69ev18mTJ8vFlAYr5kSSpkyZok8//VSvvfaaVq5cqYEDB+qDDz7QK6+8ciG3aStWzcmZ0tLSdPvttwfFWpaF+2f8xIkTlZSUpPbt2wfOXXnllZJOj+o7k38HprP7LousmJfCxMbGKi4ursjr7M6KOZkzZ47uvvvuoO9927dv14QJE1StWjXNnTv3gu8XQOgxzcTi/GsoFPRByeVyqUGDBtq8ebO+/vprde7cWdHR0ef1S8mOHTt0//33hyrcMu9889GtW7egrSHPVTg6s7jUrFkzSVJWVpZiYmIC12zfvl2XXnqpatWqFaI7KhtKIy+F/eJZq1Ytvfnmm4HnHo9Hx48f1wMPPBCiuykbzMzJiy++qCeffFKS5HA49MYbb+jnn3/WRx99VO7WAzgXM3NyprS0NE2ePLnkN1AGhfpn/OzZsxUXF6fbb789cC4nJ0dRUVH64x//qO+++y5oHYCdO3fq8ssvZ7vvs5iVl3379uUrAh4+fFhbtmzRrbfeGpJ7KSvMyslzzz2Xr5CUmJioHj166OGHHy4XRSbAThiZYWF79+7Vpk2bJClofv6Z/POTv/rqq3O2dfDgQT366KNatmxZ4Nxnn32mevXqqVOnTkHX+ncE8Pl8Fxx7WXQh+YiOjlbDhg0DR3x8/Dn72LFjhzp27Kibb75ZMTExQfmSpKVLl6p3797l5i+cxVEaeSnslxz/mhl+W7du1VNPPaWhQ4eqRYsWobol2zM7J7m5uTp16lTgfEREhNq3b19u1jIpjnB9nfz73//Wzp07ddNNN4XydsqEUP6Ml6T09HSlpaWpSpUqmjdvnubNm6c5c+aof//+iouL0+jRo/Xhhx9q7969kqQTJ05o3rx5GjduXIjuqGwwMy/btm1T9+7dA6MGDMPQ8OHD9cwzz1BgOoOZOWnevHnQ972GDRuqQoUKio2NVcOGDcvFKCbAThiZYWE///xz4IgTR/oAABStSURBVHFh0wr8P+zWr19/zrZOnDihjIwMeTweJSUlqWXLlqpbt67efvvtoOs+//xzTZw4UZI0c+ZMVapUSXfccUdJbqPMCGU+Dh48qHHjxql3796Bv3yeXVx6/PHHNWvWLN1zzz2KjIzUmjVrtH79+nKziF5xhTIv0v9+ybnnnns0b948SacXCPvkk0/03nvvBa7btGmTevToob59++rll18uyS2UOWbnpG/fvpo3b54efvjhQAFj2bJluuuuu0p6K2VGuL5O0tLS1LVr13I/nbEgoczJkiVLdNdddyk3N1cLFiwIei05OVlVqlTRHXfcoSNHjmjAgAFq1qyZsrKyNHToUN1zzz0lu5Eyxsy8xMXFKTMzUz169FCXLl1Uv3599enTh+LfWcz+WgFgHxQzLOzAgQOBx4VVgv3zXc+eB3u2uLg4rV69usg+u3btWq53NjmXUOajOMWlUaNGqVKlSkpOTlaDBg30448/6uuvvy7zu5icr1DmpTi/5Pz73//WvHnz9N133yk1NVVdunQp4R2UPWbnZOLEiRo/frz69Omja6+9Vnv27NGdd96phx9+uIR3UnaYnRO/tLQ0/fWvf73QsMu0UOakbdu2ge1xz6Vfv37q16/feURZ/pidl5UrV55nhOVPOL5WznT2qEwA1kExw8KOHj0aeOzfSu1s/q3Ujhw5YkpM5Vko81Hc4tLw4cPPI8LyKZR5KeqXnNzcXG3ZskWDBw/WU089dQHRlg9m5sTfxzPPPHOeUZYvZufErzjf58orfsZbE3mxHnICoDCsmWFhZ/51y78689n856tWrWpKTOUZ+bAmM/PicDh02223sQBYEfhasR5yYj3kxJrIi/WQEwCFoZhhYZdeemngsX9LzrP5z/PhqvSRD2siL9ZDTqyHnFgPObEm8mI95ARAYShmWNi1114b2LVi165dBV7jP9+oUSPT4iqvyIc1kRfrISfWQ06sh5xYE3mxHnICoDAUMywsNjZWTZo0kSRt2LChwGv8Kzyz8nXpIx/WRF6sh5xYDzmxHnJiTeTFesgJgMJQzLC4Pn36SJIWL16c77Xdu3fr119/VWxsrG655RazQyuXyIc1kRfrISfWQ06sh5xYE3mxHnICoCAUM8Ls2LFjkiSfz1fg64MHD1atWrX04YcfBq3mLEmzZ8+Wz+fTY489xr7YIUI+rIm8WA85sR5yYj3kxJrIi/WQEwAXxEDY7Nu3z6hTp44hyejQoYNx8uTJAq/7xz/+YVSpUsUYOHCgceLECcMwDGP58uWG0+k0brvtNuPUqVNmhl1mkQ9rIi/WQ06sh5xYDzmxJvJiPeQEwIWKMAzDCF8ppXzKy8vT9ddfr40bNwbth12rVi21adNG6enp+d6zdu1ajR8/Xps2bVLVqlV15MgRDRgwQA899JCioqLMDL/MIR/WRF6sh5xYDzmxHnJiTeTFesgJgJKimAEAAAAAAGyFNTMAAAAAAICtUMwAAAAAAAC2QjEDAAAAAADYCsUMAAAAAABgKxQzAAAAAACArVDMAAAAAAAAtkIxAwAAAAAA2ArFDAAAAAAAYCsUMwAAAAAAgK1QzAAAAAAAALZCMQMAAAAAANgKxQwAAAAAAGArFDMAAAAAAICtVAh3AAAAAEXZvXu3Ro0apfr168vr9SoqKkqVKlXS0qVL9dBDD+lPf/pTuEMEAAAmijAMwwh3EAAAoOz7+9//ro4dO+riiy8+7/d26tRJ119/vZ577jlJUuvWrTVw4EA1b95cTqdTV111VajDBQAAFsY0EwAAUOpeeukl/fTTTxdUyJCkFStWyOVyBZ7Hx8fr888/V4sWLUpUyJg3b54mTJhwwe8HAADhQTEDAFAm/PLLL2rVqpXi4uIUERGhiIgIVaxYUTfccIMSExOVmJio1q1b6+qrr1bTpk317LPPavfu3eEOu1j+/Oc/Ky4uTt9//324Q7kg7733nhYvXhwYVXEh/vKXvwTuPzc3Vz/++KN69OhR4th69uypEydO6KmnnipxWwAAwDxMMwEAlCk+n09dunTRP/7xDw0aNEizZs3Kd83ixYv16KOPavPmzZo1a1ZIPhSXlv3796tmzZqSpIcfflhTpkzJd83KlStVoUIFNW3a1OzwirRt2za1bt1aq1evVlxc3AW34/F45PP55PV6dfjwYbVr105dunQJSYw+n08dOnTQ8OHD1bVr15C0CQAAShcLgAIAypTIyMjAh+aLLrqowGvatWunH374QV26dFGvXr00d+5c9e7d28wwi61GjRoaOXKkFi9erCFDhhR4zYIFCxQfH2/JYsYjjzyiwYMHl6iQIUmrVq3SG2+8EaKogkVGRmrq1Knq1KmTfvnlF1WvXr1U+gEAAKHDNBMAQLlUsWJFzZ07VxdddJEGDRqkLVu2hDukQk2YMEGLFy/WtddeW+DrS5YsMTmi4lm3bp2++OILJScnl7itvLw8ffTRR/r999+1f/9++Xy+EET4P02aNFHr1q01adKkkLYLAABKB8UMAEC5VatWLd13333KycnR6NGjwx3OBfn3v/+tb775JtxhFGjmzJnq0KFDiUdlSNKAAQPUu3dv1a5dWzVr1lTFihV15ZVX6p133glBpKc98MADmjJlio4ePRqyNgEAQOmgmAEAKNduvfVWSae3Dd23b1+YoylYXl6e9u3bp4MHDwbOnTp1SkuXLtWdd94pqy5/NW/ePHXq1KnE7SxevFhPPfWUtm3bJsMwlJeXp71792ry5Ml6+OGHtXXr1pIHK6lDhw7Ky8vTwoULQ9IeAAAoPRQzAADl2tVXXy1JOnnyZIEjHD7++GN16NBBN998s1q0aKEbb7xRn3/+eeD1TZs26aabblLjxo11ySWX6KWXXtKhQ4f0yCOPqGvXrmrUqJESExO1dOnSfG3n5ubq2WefVZs2bdS+fXvddNNNevLJJzV8+HC9++67OnjwoGJjY1WxYkXVqlUraPFPj8ejJ598MjCK4Pnnnw/s2jJw4MDAdbNnz9a1114b2OElOjpaQ4cODbz+/vvvq3LlyoqIiFCVKlU0bNiwkv+nStq4caN+//13tWjRosRtbd68WbVr1w5szRoZGanY2Fi1b99el1xyiSIiIkrchyQ5HA61bNlSn3zySUjaAwAApYcFQAEA5VqNGjUCj3/66aeghUCHDh2qtLQ0ffPNN4Gix3vvvafbb79ds2bN0r333qu4uDiNHDlSkydP1s8//6ycnBylpKToueeeU/369XXy5EndcMMN6ty5szZu3Bg05aJfv3667LLLtGjRIlWsWFHS6fUvkpKS5PF4dPHFF2vr1q0aP368XnzxxaC4hw4dqqFDh2rs2LF65pln9OSTT+qee+7Jd38DBw7UwIEDNWzYME2ePFljx47VY489Fnj9rrvu0rFjx5SamqpFixbpsssuC8n/68qVKyVJ11xzTYnbGjBggKpUqaKHHnpIderUUVRUlI4fP64dO3Zozpw5io+PL3Effs2bN9enn34asvYAAEDpoJgBACjXqlSpEni8f//+wOO3335bU6dO1bvvvhsoZEhS37599cEHH2j48OG66667FB0drVtuuUWRkZFauHChZsyYoc8++0z169eXdHqh0R49emjUqFH66quv1L9/f0mS1+tVWlqa1qxZEyhkSFLbtm2Ddi1xOp265ZZb8hUzztcjjzyiqVOn6qOPPgoqZvjve/z48SErZEinR1NI0sUXXxyS9nr16qVevXqFpK1zqV+/vn777TcdPXpUVatWLfX+AADAhWGaCQCgXDtx4kTgcXR0dODxc889p4iICN1xxx353tOhQwcdPHhQ//rXvwLnKlQ4/feByy+/XAkJCUHXX3LJJZKkbdu2Bc5FRkaqYsWKevrpp/PtpHLHHXcEjWiIioq6kFsLEh8fr6SkJC1fvlyrV68OnDcMQ/Pnz1f37t1L3MeZ9u7dq2rVqgUVauygdu3a8vl82rlzZ7hDAQAA58DIDABAuXbo0KHA4yuvvFKStGXLFv3nP/9RVFSUunTpku89OTk5io+P17Fjx/K9dtVVV+U75y9GnDp1KnCuWrVqGj9+vB5//HF9/vnnatGihZo1a6Y777xTHTt2DNk6EGcaMmRIYPTI9OnTJUnffvut2rZtG/Kiw9GjRxUTExPSNs3gH42RnZ0d5kgAAMC5UMwAAJRrv/32W+Cxf+eN3bt3SzpdhFi0aNF5tVepUqViXztixAg1bdpUb775pr799lstX75cM2bMULNmzbRgwQLVrl37vPouSpcuXXT55Zdr7ty5eumllxQdHa1Zs2ZpwoQJIe1HOv3/4PP5Qt5uafMXdQoqVAEAAOtgmgkAoFz75z//KUnq3LmzGjRoIEmBtSNOnDhRatu1njx5Ur/88os6deqk999/X7t379a2bds0btw4/fzzzxo8ePAFtz116lQtW7Ys3/moqCgNGjRIOTk5eu+997Rnzx7l5OQE7juUYmJidOTIkZC3W9r8u8NcdNFFYY4EAACcC8UMAEC55fV6NXPmTDkcDj3//POB8y6XS9ddd50kFbilqiQtW7ZMa9asueC+d+zYkW/3kbp162r06NEaN25cgdvEFiQy8vSPcsMwAuf279+v3NzcAq8fNGiQoqKiNH36dL355ptB27iGUr169eT1egPFgbP5t4oN11EY//QS/zonAADAmihmAADKpVOnTmngwIHavXu3Xn31VTVp0iTo9fHjxysiIkJvvvlmge+fMGGCXC5XiWJYtmyZ1q5dm+98gwYNVLNmzWK14d+WdO/evYFzO3bsUJ06dQq8/rLLLlPXrl21Zs0azZkzR926dbuAyIvWqFEjScq3uKmfYRhhPQqTmZmpypUrF/r/BwAArIFiBgCgTPH5fPr9998lqdDRCd9//73atWunb775Ru+8844GDRqU75ouXbro5Zdf1vz58+XxeAIfgHNzczVs2DDdcccdqlGjRuB6r9crSTp8+HC+tvyv+f/1MwxDf/nLX/Trr78Gzu3evVsTJ07U008/HTjnHy1QUNtdu3ZVzZo19cknn8jn8+nQoUPKyspSw4YNC7x3SYGtX3v06FFqu40kJCSoUqVKQTunmGnBggV65513NGbMmMBip8WxYcMGtWjRolQWYAUAAKETYZzrzxMAANjEL7/8ogEDBmjLli3as2ePJMnhcKh58+aBRTlPnDih/fv3q3LlyurWrZsefPDBIhfZXLp0qSZOnKitW7eqRo0aqlKlilJSUgKLhf7yyy+6++67tX79+sCUikaNGumRRx5RkyZNlJKSorVr1+rEiROKjIzUH/7wB40cOVLt2rXToEGD9OKLL2r06NHatWtXYJ2GlJQU3Xnnndq0aZP69u2r9evX69ixY4qMjNR1112n6dOnq2XLloEY16xZo2HDhunAgQO69NJLNWnSJLnd7kLvyefzKT4+XosXL1a9evUu+P+8KF26dFGdOnU0a9asUuujINnZ2YqLi9OBAwdUuXJl1axZU6tWrSrWvV577bX685//rGeffbb0AwUAABeMYgYAAOXMjh07NGTIEH322Wel2s8HH3ygYcOGKSsrK7A9rVnWrVunP/zhD5Ikp9Op1atX64orrjjne7Zt26Z69epp7dq1gfcCAABrYpoJAABlWF5enr755pugnUVmzZql++67r9T77tmzpypXrqz58+eXel9n8xcjFi9erLZt2xZZyJBOF19at25NIQMAABugmAEAQBk2a9YsdezYUSNGjJB0et2NjIyMUlv480wVKlRQamqqxo8ff85FN0vLe++9p2nTpumll14q8tqTJ0/qtdde0+jRo02IDAAAlBTFDAAAyrAaNWrI4XCoTZs2On78uAYPHqwxY8YEtnQtbX369FGdOnU0c+ZMU/o7U9++fTVnzhz1799fP/744zmvnTZtmq677jrddtttJkUHAABKgmIGAABlWM+ePfXggw9qxowZ6tChg/r166cOHTqYGsO7776r2bNn66effjKlvy+++ELXX3+9JKlSpUq6+uqrC91iV5I2bdqkWbNmhaXgAgAALkyFcAcAAABKT0REhF555ZWwxuB0OrVgwQLdf//9euONN3TxxReXan8RERFq06ZN4Pm2bdsKLeAcPHhQQ4YM0fvvv6+4uLhSjQsAAIQOu5kAAICwef311/Xtt9+qbt26ioqK0pVXXqnq1atrzJgx2rRpkzZs2KD3339fhmHou+++0xNPPKEaNWpowYIF8vl8+uabbzR58mQ1adIkqN0333xThw8f1r59++RwOApdCyMjI0N169Yt1gKhAADAOhiZAQAAwiI1NVXz58/XkiVLFBUVpbffflsDBw5Udna26tevr6ZNm+qNN97Q2LFjFRkZqdq1a+u+++7TkCFDNHLkSEnSqVOnNGLECH355ZdBbd97773FiuGmm24K+X0BAIDSRzEDAACYbv/+/Ro7dqxmzpypqKgoSVJiYqImT56satWqKSYmRpJ08803BxYrrVu3rrZs2aLevXsH2nG5XNq0aZP5NwAAAMKKBUABAIDplixZohMnTuiqq64KnKtbt67+7//+L+i6+Pj4wGN/UePsc3l5eaUcLQAAsBqKGQAAwHQ+ny/cIQAAABujmAEAAEx3ww03qEKFClq/fn3Q+b1792rbtm1higoAANgFxQwAAGC6uLg4jRo1Sq+++qqOHTsWOD9lyhRVr15d/s3Wztx0rbBzbMwGAED5w9asAAAgbKZNm6aFCxfK7XarQoUK6tGjhzIzMzVhwgStWLFCCQkJeuKJJ7R371797W9/04YNG9S+fXuNGzdOGRkZmj17trZt26Zbb71Vr7zyiq655ppw3xIAADABxQwAAAAAAGArTDMBAAAAAAC2QjEDAAAAAADYCsUMAAAAAABgKxQzAAAAAACArVDMAAAAAAAAtkIxAwAAAAAA2ArFDAAAAAAAYCsUMwAAAAAAgK1QzAAAAAAAALZCMQMAAAAAANgKxQwAAAAAAGArFDMAAAAAAICtUMwAAAAAAAC2QjEDAAAAAADYCsUMAAAAAABgKxQzAAAAAACArVDMAAAAAAAAtvL/AaUe9pk/l2IWAAAAAElFTkSuQmCC\"><br>"
],
"text/plain": [
"<yt.visualization.profile_plotter.PhasePlot at 0x10bb6b080>"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"phase"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"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.5.2"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment