From da638564663c0b9d6217fd90eb7fcb0ac002acf8 Mon Sep 17 00:00:00 2001 From: L-Nafaryus Date: Fri, 4 Mar 2022 17:19:56 +0500 Subject: [PATCH] Mod: data files --- data/analyze.ipynb | 341 +++++++++++++++++++++++++++++++++++++++++++++ data/anisotropy.db | Bin 0 -> 188416 bytes 2 files changed, 341 insertions(+) create mode 100644 data/analyze.ipynb create mode 100644 data/anisotropy.db diff --git a/data/analyze.ipynb b/data/analyze.ipynb new file mode 100644 index 0000000..433192e --- /dev/null +++ b/data/analyze.ipynb @@ -0,0 +1,341 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [], + "source": [ + "from anisotropy.database import Database, tables\n", + "import pathlib\n", + "import peewee as pw" + ] + }, + { + "cell_type": "code", + "execution_count": 14, + "metadata": {}, + "outputs": [], + "source": [ + "db = Database(pathlib.Path(\"anisotropy.db\").resolve())\n", + "savefig = False" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [], + "source": [ + "execution = 5\n", + "data = \"flowRate\"\n", + "\n", + "if not db.getExecution(execution):\n", + " print(\"Execution not found\")\n", + "\n", + "for model in db.tables:\n", + " try:\n", + " column = getattr(model, data)\n", + " \n", + " except AttributeError:\n", + " pass\n", + "\n", + " else:\n", + " break" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": {}, + "outputs": [], + "source": [ + "query = model.select(tables.Shape.alpha, column, tables.Shape.direction, tables.Shape.label)\n", + "idn = db.tables.index(model)\n", + "\n", + "for table in reversed(db.tables[ :idn]):\n", + " query = query.join(table, pw.JOIN.LEFT_OUTER)\n", + " \n", + "query = query.switch(tables.Shape)\n", + "query = query.where(\n", + " tables.Shape.exec_id == execution,\n", + "# tables.Shape.label == structure,\n", + ")\n", + "query = query.order_by(tables.Shape.label, tables.Shape.direction, tables.Shape.alpha)" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": {}, + "outputs": [], + "source": [ + "with db:\n", + " if query.exists():\n", + " table = []\n", + " for row in query.dicts():\n", + " for k in row.keys():\n", + " if type(row[k]) == list:\n", + " row[k] = str(row[k])\n", + "\n", + " table.append(row)\n", + " \n", + " else:\n", + " table = None\n", + "\n", + "if not table:\n", + " print(\"Results not found\")" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": {}, + "outputs": [], + "source": [ + "from pandas import DataFrame\n", + "\n", + "df = DataFrame(table)" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/tmp/ipykernel_33356/2997525597.py:2: UserWarning: Boolean Series key will be reindexed to match DataFrame index.\n", + " df[df.label == \"simple\"][df.direction == direction][\"alpha\"].reset_index(drop = True)\n", + "/tmp/ipykernel_33356/2997525597.py:6: UserWarning: Boolean Series key will be reindexed to match DataFrame index.\n", + " df[df.label == \"simple\"][df.direction == direction][\"flowRate\"].reset_index(drop = True)\n" + ] + } + ], + "source": [ + "a_s = [\n", + " df[df.label == \"simple\"][df.direction == direction][\"alpha\"].reset_index(drop = True)\n", + " for direction in ['[1.0, 0.0, 0.0]', '[0.0, 0.0, 1.0]', '[1.0, 1.0, 1.0]']\n", + "]\n", + "k_s = [ \n", + " df[df.label == \"simple\"][df.direction == direction][\"flowRate\"].reset_index(drop = True)\n", + " for direction in ['[1.0, 0.0, 0.0]', '[0.0, 0.0, 1.0]', '[1.0, 1.0, 1.0]']\n", + "]\n", + "k_s2 = 2 * k_s[1] / k_s[0]\n", + "k_s3 = 2 * k_s[2] / k_s[0]" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/tmp/ipykernel_33356/723496196.py:2: UserWarning: Boolean Series key will be reindexed to match DataFrame index.\n", + " df[df.label == \"bodyCentered\"][df.direction == direction][\"alpha\"].reset_index(drop = True)\n", + "/tmp/ipykernel_33356/723496196.py:6: UserWarning: Boolean Series key will be reindexed to match DataFrame index.\n", + " df[df.label == \"bodyCentered\"][df.direction == direction][\"flowRate\"].reset_index(drop = True)\n" + ] + } + ], + "source": [ + "a_b = [\n", + " df[df.label == \"bodyCentered\"][df.direction == direction][\"alpha\"].reset_index(drop = True)\n", + " for direction in ['[1.0, 0.0, 0.0]', '[0.0, 0.0, 1.0]', '[1.0, 1.0, 1.0]']\n", + "]\n", + "k_b = [ \n", + " df[df.label == \"bodyCentered\"][df.direction == direction][\"flowRate\"].reset_index(drop = True)\n", + " for direction in ['[1.0, 0.0, 0.0]', '[0.0, 0.0, 1.0]', '[1.0, 1.0, 1.0]']\n", + "]\n", + "k_b2 = 2 * k_b[1] / k_b[0]\n", + "k_b3 = 2 * k_b[2] / k_b[0]" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/tmp/ipykernel_33356/1402688131.py:2: UserWarning: Boolean Series key will be reindexed to match DataFrame index.\n", + " df[df.label == \"faceCentered\"][df.direction == direction][\"alpha\"].reset_index(drop = True)\n", + "/tmp/ipykernel_33356/1402688131.py:6: UserWarning: Boolean Series key will be reindexed to match DataFrame index.\n", + " df[df.label == \"faceCentered\"][df.direction == direction][\"flowRate\"].reset_index(drop = True)\n" + ] + } + ], + "source": [ + "a_f = [\n", + " df[df.label == \"faceCentered\"][df.direction == direction][\"alpha\"].reset_index(drop = True)\n", + " for direction in ['[1.0, 0.0, 0.0]', '[0.0, 0.0, 1.0]', '[1.0, 1.0, 1.0]']\n", + "]\n", + "k_f = [ \n", + " df[df.label == \"faceCentered\"][df.direction == direction][\"flowRate\"].reset_index(drop = True)\n", + " for direction in ['[1.0, 0.0, 0.0]', '[0.0, 0.0, 1.0]', '[1.0, 1.0, 1.0]']\n", + "]\n", + "k_f2 = 2 * k_f[1] / k_f[0]\n", + "k_f3 = 2 * k_f[2] / k_f[0]" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": {}, + "outputs": [], + "source": [ + "k_f3[20] = (k_f3[18] + k_f3[21]) * 0.5\n", + "k_f3[19] = (k_f3[20] + k_f3[18]) * 0.5\n", + "k_f2[19] = (k_f2[20] + k_f2[18]) * 0.5\n", + "\n", + "k_b2[24] = (k_b2[25] + k_b2[23]) * 0.5\n", + "k_b3[24] = (k_b3[25] + k_b3[23]) * 0.5" + ] + }, + { + "cell_type": "code", + "execution_count": 12, + "metadata": {}, + "outputs": [], + "source": [ + "import matplotlib.pyplot as plt" + ] + }, + { + "cell_type": "code", + "execution_count": 17, + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAqXUlEQVR4nO3de7wcdX3/8dcbCARJCJBAKgQIQUAwBSFBQWlJjOUetGopWCkgGkEe8rO2PqhAQRFv9UJBxIgKEUUOUZACCgiYgIrUEiABjZQYQILKJUjgBAghfH5/zJywOdnZM3POzl7OvJ+Pxz7O7tz2893dM5/5XmZGEYGZmVXXBu0OwMzM2suJwMys4pwIzMwqzonAzKzinAjMzCrOicCsTSSNaHcMZuBEYNYyksZKukjSg5KWA9e0OyYzgI3aHYANL5IeBsYDa2ombwTcFREHtCWoDpAe/c8DbgH2jYhn2huR2atcI7AyzIyIUX0P4KR2B9QB/gl4LCI+5iRgncaJwFpO0sOSPiHpt5L+IulSSSNr5n9Q0hJJT0u6VtK2NfPeIOnmdN7jkk6XtL+k3vSxWtJLNa93kLSzpJ9JWi7pKUmXS9qiZpsHSLpX0nPpOq9ImpYR+3xJH0if7yzpUUlHSto4jemva5bdRtLzkrYG3gSskfQbSSvS7eze7zN5e/p8VFq2X6Svp0laVrNs/9fbSrpK0pOSHpJ0as28DdPP6Pdp+RZI2j6dF5Jelz4/UtIfJE1KXx8u6R5Jz6Zl/GS/z+FDkh5JP6+VknyJgi7mRGDt8k/AwcDOwK7AmQCS3gZ8DjgKeC3wCNCTzhtN0rRyI7At8Drg1oj4VU3t43LgP2tqJH8AlG5zW2B3YHvgkzWxfAn4EbB5uo0/DhS8pNcCNwGnR8S1EfFSGuf7ahY7Jo3vSeA1wN8BHwXGAT8BrpO0cZ3NfxxYXfP6FTL+VyVtAFwHLAS2A2YAH5V0cLrIx9I4DgM2B94PPN9vGwcCs4HDI2JpOnkl8M/AFsDhwMmS3pkuvxlwEXBc+nntVS826x5OBNYuF0bEoxHxNPAZkp0VJAnikoi4OyJWAZ8A9pc0ETgC+HNEfDkiXoyI5yLifwZ6o4hYEhE3R8SqdKf8FeDAfottSJIw8tiSJAlcHhHfrZn+HeAYSX3bORaonf/faRyrSZLPpsBbajcs6a+AE9MY+zwKbCOp3g53X2DriDgnIl5Kd+TfBI5O538AODMiHojEwohYXrP+3sC1wD9FxH19EyNifkTcFxGvRMQi4Ape/cw2IElO7mMcJpwIrF0erXn+CMnROunfR/pmREQvsJzkaHd74PdF30jSeEk9kh6T9CzwPZKj8j4fAY4EXpT0TE0sWT4F9AJvq9npkyal54Fpkl5PUmO5Np29ql+5XiH5DLbrt+2zga8CT9cs+xBwDnBzGt/1NcvvCGwr6Zm+B3A6SYc9DPyZfQt4kKS2spakN0ualzY3rSDp5xmXxvMcSbK6TNLzwN0Ntm9dwInA2mX7muc78GpzzB9Jdm7A2maIscBjJDvOSYN4r88CAfx1RGxO0nxTuwP/X5Kd9BkRsQUDNw3NBQ5It/GRfvO+k27/WOCHEfFiOv0P/colks/gsZp1dyVpLju//xumR/zbpPEdUTPrUeChiNii5jE6Ig6rmb9zg7J8NN3eiZL2qZn+fZIktn1EjCFpOqqtMV1D0nw1Hahdz7qQE4G1yymSJkjaCjgDuDKdfgVwgqQ3StqEZCf+PxHxMMmR8GslfVTSJpJGS3pzjvcaTXIEv0LSdiRt8GtJOookGZ2XM/ZfpEf07wfO6utgTX0P+HuSZHBZzfS5wOGSZigZSvqvJLWEO2qWORM4pyZ55PFr4DlJp0naNO0cnixp33T+t4BPS9pFiT0lja1Z/+cR8Wfg34BL9epJbqOBpyPiRUlvAt7b730/D1ybp2nOOp8TgbXL94GfAktJmi7OBYiIW4D/AK4C/kRyNHt0Ou85kiaMmcCfSZo0pud4r0+RHLWuAH4MXN03Q9KWJAnggxHxcpECRMT/kewQv9XXRBQRj5I0lQTw85plHySpJXwVeCotw8y0k7nPU6ybPPLEsIbkiP6NwEPpNr4FjEkX+QpJEvop8CzwbZK+if7b+S5J7eH0dNKHgXMkPQeclW4DAElvJelAPr3/dqw7yTemsVZTctLZB9Kd/rAj6RLgjxFxZrtjMcvDvf5mTZSObnoXyWgcs67gpiGzJpH0aeB+4IvpSB+zruCmITOzinONwMys4rquj2DcuHExceLEta9XrlzJZptt1r6AWqQq5YTqlLUq5YTqlLWTy7lgwYKnImLrevO6LhFMnDiRu+66a+3r+fPnM23atPYF1CJVKSdUp6xVKSdUp6ydXE5Jj2TNc9OQmVnFORGYmVWcE4GZWcV1XR9BPatXr2bZsmW8+GKRS7R0lzFjxrB48eIBlxs5ciQTJkxgxAjfF93M8hkWiWDZsmWMHj2aiRMnUnNV4GHlueeeY/To0Q2XiQiWL1/OsmXL2GmnnVoUmZl1u2GRCF588cVhnQTyksTYsWN58skn2x2KmTXJ1HNv5qnel9abPm7Uxtx15t/VWaO4YdNHUPUk0Mefg9nwUi8JNJo+GMMmEZiZ2eA4EZiZVZwTgZlZxTkRNNGtt97KscceW3i9k046iV/+8pdN256ZWRHDYtRQEWX2wC9cuJC99y5+P5I777yTr33ta03bnpkNH+NGbZy5z2qWyiWCMnvgFy5cyPHHH8+qVav40Ic+xLbbbstnPvOZhiN5Fi9ezK677sqGG27YcHsnnXQSO+6444DbM7PhpVlDRBtx01ATLVq0iG222YaDDz6Yt7/97Xz2s58dcKd9ww03cMghhwy4vWnTpuXanplZUZWrEZRl9erVLF26lGOOOYZvfOMb7L///lxzzTX8+Mc/5tlnn+XEE0/koIMOWm+9m266iUsvvXTA7U2ePBmApUuX8pnPfIYVK1bwwx/+sPRymdnw5xpBkyxevJh9992XjTbaaG0zzzvf+U6++c1vMnv2bK688sr11nn++ed55pln2HbbbXNtD2DSpEl8+9vfLq8gZlY5pSUCSdtLmifpt5J+I+n/1VlGki6QtETSIkn7lBVP2RYuXMhb3vIWenp6OOGEE3j88cfXzjv33HM55ZRT1ltn3rx5TJ8+Pdf2nnjiidJiN7NqK7NG8DLwrxGxB7AfcIqkPfotcyiwS/qYBXy9xHiA7J72ofbAL1y4kMmTJ7PrrrvyhS98gaOOOorVq1dz2mmnceihh7LPPuvnuEb9A/23d9xxx7F69eohxWhmVk9pfQQR8SfgT+nz5yQtBrYDfluz2DuAyyIigDslbSHptem6pSirB/5LX/rS2udHHHEERxxxBBdccAG33HILK1asYMmSJZx00knrrHPHHXdw3nnn5dregQceyIgRI1i+fDlnnHEG99xzD5/73Of4xCc+UUp5zKw6lOyDS34TaSJwOzA5Ip6tmX498PmI+EX6+lbgtIi4q9/6s0hqDIwfP35KT0/P2nm9vb1st912vO51ryu9HO20Zs2aukNM61myZAkrVqwoOaLy9Pb2MmrUqHaHUbqqlBOqU9ZOLuf06dMXRMTUevNKHzUkaRRwFfDR2iRQRERcDFwMMHXq1Ki9OfT8+fMZOXLkgNfq73Z57kfQZ+TIkV19Ilon3wC8mapSTqhOWbu1nKWOGpI0giQJXB4RV9dZ5DFg+5rXE9JpZmbWImWOGhLwbWBxRHwlY7FrgX9ORw/tB6wos3/AzMzWV2bT0FuBY4H7JN2bTjsd2AEgImYDPwEOA5YAzwMnlBiPmZnVUeaooV8ADa+HkI4WWn+AvZmZtYzPLDYzqzgnAjOzinMiMDOrOCcCM7OKcyJoIt+q0sy6UTUTwaK5cN5k+OQWyd9Fc5uy2aHcqnK//fZr2vbMzIqoXiJYNBeuOxVWPApE8ve6U5uSDPp23KtWreL444/n9NNPZ6BrOQ10q8q+7Z100km5tmdmVlT1EsGt58DqF9adtvqFZPoQ+VaVZtaNBjyhTNJ9QO1hqEjOBduztKjKtGJZsek51btV5eLFizn//PN56qmnmDFjBieffPJ66/lWlWbWbnlqBAuBNcAZwEzgiPRvdxozodj0nOrdWnL33Xdn9uzZzJ07t25nsG9VaWadYMBEEBHvI7lm0InAucAGEfFI2YGVZsZZMGLTdaeN2DSZPgRZt6q89tprOfzwwznssMPWW8e3qjSzTjBgIpC0FfBH4P3AXOAHki4sO7DS7HkUzLwAxmwPKPk784Jk+hBk3aryyCOP5IYbbuDyyy9fbx3fqtLMOkGei84t4NU+gr6eyvUPb7vJnkcNecffX71bVc6fP5+rr76aVatW1a0R+FaVZtYJBkwEEbFTKwIZjqZNm9bwbkV333134W2OHTuW2bNnDyEqM7N15WkaKr63MjOzrpFn1JAHrpuZDWN5+gh2k7So5nV3n0dgZmbryJMIHqILzhuICJ91C74EhZkVlicRvNTp5w2MHDmS5cuXM3bs2Eong4hg+fLljBw5st2hmFkXyZMIPlJ6FEM0YcIEli1bxpNPPtnuUErz4osv5trBjxw5kgkThnaWtJm13tRzb+ap3pfWmz5u1MbcdebflfreeRLBYZLuj4hnACRtCfxrRJxZamQFjBgxgp12Gt6jXOfPn+9LUpsNY/WSQKPpzZRn1NChfUkAICL+Qo4TyiRdIukJSfdnzN9S0o8kLZL0a0mTc0dtZmZNkycRbChpk74XkjYFNmmwfJ85QP3rJyROB+5NRx/9M3B+jm2amVmT5UkElwO3SjpR0onAzcB3BlopIm4Hnm6wyB7Az9JlfwdMlDQ+RzxmZtZEyjPcUNIhwNvTlzdHxE25Ni5NBK6PiPWafSR9Ftg0Iv5F0puAO4A3R8SCOsvOAmYBjB8/fkpPT8/aeb29vYwaNSpPOF2tKuWE6pS1KuWE6pR1KOU8/saVmfPmHLLZYENaa/r06QsiYmq9eXk6iwHuAUaQXHzuniFHlPg8cL6ke4H70u2uqbdgRFwMXAwwderUqL1+z/z58xtez2e4qEo5oTplrUo5oTplHUo5x/0ie9RQ2Z9dnjuUHQV8EZhPclbxVyV9PCKGdHusiHgWOCF9D5GcuLZ0KNs0M+tWZQ8RbSRPjeAMYN+IeAJA0tbALcCQEoGkLYDnI+Il4APA7WlyMDOzFsqTCDboSwKp5eS7aukVwDRgnKRlwNkkzUtExGxgd+A7kgL4Dckd0MzMrMXyJIIbJd0EXJG+/kfgJwOtFBHHDDD/V8CuOd7fzMxKlOfGNB+X9G7gremkiyPiR+WGZWZmrZJr1FBEXAVcVXIsZmbWBnlGDb0C9NZOIrkfwealRWVmNky18+JyWfLUCL4I7A3cCHyvX8exmZkV0M6Ly2XJ00dwmqQNgEOBr0n6S0TMKj80MzNrhTzXGoJkdM+BwBjg1+WFY2ZmrZanj+AOkn6BS4GedNo+EXF3ybGZmVkL5LpVJck1ho4BjibtLAbeVmJcZmbWInn6CKa1IA4zs0oYN2rjzFFD7ZKnaWg88Flg24g4VNIewP4R8e3SozMzG2baeXG5LHk6i+cANwHbpq//D/hoSfGYmVmL5UkE4yJiLvAKQES8TMZ9A8zMrPvkSQQrJY0l6SBG0n7AilKjMjOzlskzauhjwLXAzpJ+CWwNvKfUqMzMrGXyjBq6W9KBwG4kQ0cfiIjVpUdmZmYtkfeexa8H9gBGAvtIIiIuKy8sMzNrlTzDR88mudPYHiQ3pDkU+AXgRGBmNgzk6Sx+DzAD+HNEnADsRXLNITMzGwbyJIIXIuIV4GVJmwNPANuXG5aZmbVKnj6CuyRtAXwTWEByk5pflRmUmZm1Tp5RQx9On86WdCOweUQsGmg9SZcARwBPRMTkOvPHAN8Ddkjj+FJEXFok+KHoxLsEmZm1Q95RQwBExMMAkv6d5B4FV0XEjzMWnwNcSHan8inAbyNipqStgQckXR4RLblNT6O7BE389/WL1HfJ1f6cOMys2+UZNfQQ6+4DBYwnGVL6TNZ6EXG7pIkNNh3AaEkCRgFPAy8PHHIxWUf+RdVLAtDe28uZmTWDIrJ2cekCyeUl1pkEXBcR+w+48SQRXJ/RNDSa5Izl1wOjgX/Mql1ImgXMAhg/fvyUnp6etfN6e3sZNWpUZgzH37hyoDBLkVWD2HxjuOBtmxXe3kDlHE6qUtaqlBOqU9ZOLuf06dMXRMTUevPy9BEs7z9NUjPOLD4YuJfkBjc7AzdL+nlEPFsnhouBiwGmTp0a06ZNWztv/vz51L5ez41ZLVflykqvz75E43gzDFjOYaQqZa1KOaE6Ze3Wcg44fFTSVv0fFOxbyHACcHUklgAPkdQOzMyshfLs0BeQHOCqZlrj9qR8/kByotrP05vf7AYsbcJ2O169zmh3OpsNL900MjFP09BOg9mwpCtILk0xTtIy4GxgRLrN2cCngTmS7iNJMqdFxFODea/Bevjzh683LevLy2rzb5Z677leLGkzV7NGMHXTD9Ws2zQamdhp8owaele96RFxdaP1IuKYAeb/EThooPcfqqL3By26A6x3dN8sWT+YRiOYitQ2uumHWoZuToTdHLt1njxNQ1cCvyVpIuprHgqgYSLoFGX/U2QlmsHUIMpKKlXZsRfVzYmwm2O3zpMnEUwmacYZBfxHRDxQbkjdpZNqEJ34vp2gWeeSmA1XefoIHgCOkjQF+IqkPwKfjIjHSo/Ohp12NGk0KwkUid1NN+3hz31w8vQRfJVXWzmWAgcCDwKvKTGuYSurKanT5K1BFP0Ha0aTRrv+2YvEnvsSJk0eAADdOyqtGYk2Szv+54r2T7ZTrquPDvDaCqj3z1i02absEUxFZP2DDaY5psgOMiuWZjWBtboprexLmHTDwUczEm0n6fTEWytP09B3WhGI1dc3xDXPGYuD2fnWG0JbdCdY5k6zUxJeOxW5CGLR7XRDTQGq3cfVCoO96FxExKTSoqqYZlUhs/6h/U+UrRmJsB2C5sTeDUfWnS7rfJ9uSbKQr2loKsnO/2fA9HLDqaZ2DXHtxLbKVqp6+S1bM04s7aYkm/uic5JerncBOut83XJUUqZ6R8/NIrqjFlGW3DvNNh4pZyX9os2pw7WpMk/T0Fbp0w0lbUl6UllEPF1mYNY+ZY9s6tbmGCgee//ly76ESTtGpQ3mDPh6mhF7VpKZeu7NXfMba4eiF527O50WgPsIhqkiR2tF/7myjsyKnqFddOx+0RjLakqrjTnPAICin+9g+on6z2u0M21WksnquO6fOIsk2Ua6qZmmHUq76JxVQ6OdZpGEUnQHOdA2hqKTmtKalZSKHG236/IVWeV0/1b58jQNjQBOBv42nTQf+EZENOPmNNblOmmnORyVmdy6oamkm39f3ZSs8jQNfZ3k8tEXpa+PTad9oKygzKy9uiFJtMNATZXdeoeyPIlg34jYq+b1zyQtLCsgs27j5ov6OukM+KK66RyAZsiTCNZI2jkifg8gaRKwptywzLpHlXYYWQbquK09Um5HbaNZfVnDVZ5E8HFgnqSlJEl+R5L7DZtZF2vWUNNOGpWVxTv7xvKMGrpV0i4k9xQGeCAiVpUblpmVbTBDTZtxYp53yp0nz6ihkcCHgQNImvx+Lml2RLxYdnBmZla+DXIscxnwBuCrwIXp8++WGZSZtU+jk/5seMp1q8qI2KPm9TxJvy0rIDNrLzfdVE+eGsHdkvbreyHpzeS4OY2kSyQ9Ien+jPkfl3Rv+rhf0pqa6xqZmVmL5EkEU4A7JD0s6WHgV8C+ku6TtKjBenOAQ7JmRsQXI+KNEfFG4BPAbb6QnZlZ6+VpGsrcmTcSEbdLmphz8WOAKwbzPmZmNjSKaHzun6Qd6k2PiD8MuPEkEVwfEZMbLPMaYBnwuqwagaRZwCyA8ePHT+np6Vk7r7e3l1GjRg0USterSjmhOmWtSjmhOmXt5HJOnz59QURMrTcvT43gd8AS0vsQpALYswmxAcwEftmoWSgiLgYuBpg6dWrUXsujW6/tUVRVygnVKWtVygnVKWu3ljNPInggIvYuMYajcbOQmVnb5OksLu26UZLGAAcC/13We5iZWWN5agRbSHpX/4kRcXWjlSRdAUwDxklaBpxNcjlrImJ2utjfAz+NiJVFgjYzs+bJkwhuI2nHrxVAw0QQEccMtOGImEMyzNTMzNokz0XnfKVRM7NhLE8fgZmZDWNOBGZmFedEYGZWcXnuR3BWvekRcU7zwzEzs1bLUyOYBfQCK4EPpn893NPMbJjIM3z0yYj4CoCkk4DLI+LP5YZlZmatkqdGIEl/L+k4YBXwI0nTSo3KzMxaJk+N4IPAacAa4B+Ax4FLgfnlhWVmZq2S54Sy/wXe02/yO8oJx8zMWi3PqKGP1Zve129gZmbdLU/T0H8AjwA/KjkWMzNrgzyJYGeSewrPAM6JiFvKDcnMzFppwFFDEfF0RHyc5AYy/yDpRkn7lh+amZm1Qp4+gut49eY0AnYA7gQ2LDEuMzNrkTxNQ18qPQozM2ubPMNHb2tFIGZm1h55moZeBp6vnQRERGxeWlRmZtYyeZqG7ouIvUuPxMzM2iLPtYZi4EXMzKxb5akRvEbS3iRNQmtFxN2NVpJ0CXAE8ERETM5YZhrwX8AI4KmIODBHPGZm1kR5EsGfgP6XkwjgbQOsNwe4ELis3kxJWwAXAYdExB8kbZMjFjMza7I8o4amD2bDEXG7pIkNFnkvcHVE/CFd/onBvI+ZmQ2NIhp3AUgaA5wN/G066TaSS02sGHDjSSK4vl7TkKT/ImkSegMwGjg/IrJqD7NI7pTG+PHjp/T09Kyd19vby6hRowYKpetVpZxQnbJWpZxQnbJ2cjmnT5++ICKm1p0ZEQ0fwFXAp4BJ6eNskiP5POtOBO7PmHchyRnKmwHjgAeBXQfa5pQpU6LWvHnzogqqUs6I6pS1KuWMqE5ZO7mcwF2RsV/NddG5iHh3zetPSbq3cDpa3zJgeUSsBFZKuh3YC/i/JmzbzMxyyjN89AVJB/S9kPRW4IUmvPd/AwdI2kjSa4A3A4ubsF0zMysgT43gZOA7aV+BgKeB4wZaSdIVwDRgnKRlJE1KIwAiYnZELJZ0I7AIeAX4VkTcP6hSmJnZoOUZNXQvsJekzdPXz+bZcEQck2OZLwJfzLM9MzMrx4BNQ5LGSrqA5Gb18ySdL2ls6ZGZmVlL5Okj6AGeBN5NchP7J4ErywzKzMxaJ08fwWsj4tM1r8+V9I9lBWRmZq2Vp0bwU0lHS9ogfRwF3FR2YGZm1hp5EsEHge8Dq9JHD/AhSc9JytVxbGZmnSvPqKHRrQjEzMzaI8+ooetbEYiZmbVHnqahbUuPwszM2ibPqKFJkq7tPzEijiwhHjMza7E8ieBJ4MtlB2JmZhkWzYVbz4EVy2DMBJhxFux5VNM2nycR9EbEbU17RzMzy2/RXLjuVFidXutzxaPJa2haMsjTR/C5pryTmZkVd+s5ryaBPqtfSKY3SZ4awTWSTmXdO5TNjojVTYvCzMzqW7Gs2PRByFMjuAiYkv69CNgH+HrTIjAzs2xjJhSbPgh5agT7RsReNa9/Jmlh0yIwM7NsM85at48AYMSmyfQmyVMjWCNp574XkiYBa5oWgZmZZdvzKJh5AYzZHlDyd+YFLR819HGS+xAsTaJgR+CEpkVgZmaN7XlUU3f8/eW51tCtknYBdksnPRARq0qLyMzMWirPtYZ2Bc4iuWH9QcDsdJqZmQ0DefoIvpv+vRnYEFgIXFpaRGZm1lJ5EsFGEXEG8EpEfCEi/gsYVW5YZmbWKnkSwSvp31NqpmmglSRdIukJSfdnzJ8maYWke9NH88ZCmZlZbnlGDb0fICJuAJC0OUmfwUDmABcClzVY5ucRcUSObZmZWUkGrBFExH39Xj8bEdfkWO924OnBh2ZmZq2giGi8gPQc0LdQX5NQRMTmA25cmghcHxGT68ybBlwFLAP+CPxbRPwmYzuzgFkA48ePn9LT07N2Xm9vL6NGDf8ui6qUE6pT1qqUE6pT1k4u5/Tp0xdExNS6MyOi4QM4F7gDOHygZeusOxG4P2Pe5sCo9PlhwIN5tjllypSoNW/evKiCqpQzojplrUo5I6pT1k4uJ3BXZOxX8zQNnQm8AzhY0q2S3jrUzJRu99mI6E2f/wQYIWlcM7ZtZmb5DdhZLGmf9OkcYCfgIkmPxhA7eSX9FfB4RISkN5H0VywfyjbNzKy4PKOG+t+m8mlgs4FWknQFMA0YJ2kZcDYwAiAiZgPvAU6W9DLJWctHp9UXMzNroTzXGpo+mA1HxDEDzL+QZHipmZm1UZ6mobrnDERE8+6TZmZmbZPnzOKV6eODNc9XlhmUmZm1Tp6moS8DSHpf33MzMxs+8tQI+rgj18xsGMrTR3AdSRKYJOnavukRcWSZgZmZWWvkGT76pfSvm4XMzIahPH0Et9W+lnQAcAxwW/01zMysm+SpESBpb+C9wD8AD5FcLM7MzIaBzESQ3pf4mPTxFHAlydVKB3WCmZmZdaZGNYLfAT8HjoiIJQCS/qUlUZmZWcs0Gj76LuBPwDxJ35Q0gxy3qDQzs0FaNBfOmwyf3CL5u2huS942MxFExDURcTTwemAe8FFgG0lfl3RQS6IzM6uKRXPhulNhxaNAJH+vO7UlySDP/QhWRsT3I2ImMAG4Bzit9MjMzKrk1nNg9QvrTlv9QjK9ZEXOLCYi/hIRF0fEjLICMjOrpBXLik1vokKJwMzMSjJmQrHpTeREYGbWCWacBSM2XXfaiE2T6SXLdUKZmZk10aK5Sdv/imXJEf+Ms2DPo5J5WdNL5ERgZtZKfaOD+jqG+0YHQbLTb8GOvz83DZmZtVIbRwdlcSIwMytiqCd9tXF0UJbSEoGkSyQ9Ien+AZbbV9LLkt5TVixmRrlnrbbpjNiWx9KMk77aODooS5k1gjnAIY0WkLQh8AXgpyXGYWZlnrXaxjNiWx5L0WadekmpjaODspSWCCLiduDpARb7CMklrZ8oKw6zyqm38ymzXbqT2rwbxdKMmkKRZp2spAQw8wIYsz2g5O/MC9rSSdxHEeXdiljSROD6iJhcZ952wPeB6cAl6XI/zNjOLGAWwPjx46f09PSsndfb28uoUaOaH3yHqUo5YfiXdZvHb2PS0u+yyaonWbXJ1iyddCxPjD9wCNt5ilWbjGPppGMB2O2Br7HhK6vWLrdmg03Y4JVVda8YGcCqTbZeZxuNYqn3nrsvPg/VuaV5IG6bdg3Quu/0wPnvzIgFXtlgk/U+lwd2OwVgvTJlfQb7/eoDjFz15HrTX9xka+7c/1vrlHOgZVtt+vTpCyJiar157UwEPwC+HBF3SppDg0RQa+rUqXHXXXetfT1//nymTZvWtJg7VVXKCcO8rP2HDkLSLDDzguR53jHkWdvZaFN4oU5FXBtCrKmzIUHtjrNRLFDsPcdsD/+SdBHm/k4bja/Ps+yt56RH4P2LmVH+TbeCl1+o/33Ue9+sz32v98KDPyVWLEN9sVw9C+okJRB88pkGH0I5JGUmgnaeRzAV6JEEMA44TNLLEXFNyyIo8qMza4asposbTlt3h1TbjNC3Xv8dXr3t9J/WJ9YkO6x15vdLAgPFslH/9dPlN9p0/W33tXmn/2MHrlgG99T8j9X734PG4+trZY3F3+u9sPD768eS9bnUS2C1zVp5Tvra5aC176naWDbdMiNBtq9TOEvbagT9lptDmTWCPD86GPjIrGjiaGKiGdZHyX3Sz2udo6pmJeas76Lo9KLb7++TW1D/KDFD1hFr1o4ty5jta46Y0xjrHTkPiuBdF+f/H8vaWeeoWax13uT68dcrZ6OaQiP1klu9mkJWLEVrGyVrVCMoLRFIugKYRnK0/zhwNjACICJm91t2DmUlgqJV6KwvL+vHW7Qa3Wj5BjuSUqrWg1m+iCI7Wcj+vIYaT6PqfL3vtFnfddpcMOQdUj3NaOrI2oEVVW9H3Wj7mc1UA7xH7edYtNmlWU1p9cqamdwzEmSbWh3akgjKUjgRNOvHXvQfr0ii2WAESLDmpXW3UbPjWe8ouRk700bt1a3e+RY5Guzbft4kU7TduFnfdb3296LlbyTriBXy7XyadZCU9XspWvvJVOdzLPp7gWL/M5k1rjqJplHtJCuWNqh2Imjaj7ENitZOGv1zFKkuZy0P5e18M9X5xyuaZIo2o5Sp5rNdJ7lD8Y7Yet9R0QTerGbTeorWCOr93uv1Y2QtO9iDmCK/33o79zIPqJqo2omgaPtd4VEXbTCYWOodPTbaQfafn1VrKXvnW2+H99LKYt9Rs6Y3xauJbb3fbqubzBppRpNh0YRdL9Fk1uZLbnYpunMvs3+rSTp11FBrzDir/hd66BeS52V1cBVuRiig6E5KG9Yf7dFoR9h/+VdWr7/c6hdgwZz1tzHQtos0u+xy0PqjQ7JkfS71RswMpo+gyHeddSTbaMRIoytPtrqduRlXwawZYbPeDnKH/bLLVPu+mc0uE8q9UmfRS0KnsdzWpYM6hn8iGOgLLfKPV+/HC0NPNFlH20VrJ0VHmGTtIIsczTdr51vzea2z06g3TDJLow6+rGaUrB3SUL/rmiGF6yw7mMsItOnSxE2RtYPMW6asA7lWXI6hmz/3goZ/IoDiX2jW8oM5Ysu7fNa0IrWTrORTtC+gyMiWZu58089rnZ3G1bPyxdHoc+nbfpHvtBnfdaNyWj5tvFlLlVQjEZStmYmmnqFUrftkHVVlxZK31tLMnW89We3Em24FG29WLMk0Q5HYK3REWSp/jqVzIuh0Q61a9y0Lxdo76y2ftY0yd76N+niGmmTMDHAiqI4yay2d1GlnZoU5EVjn81G+Wak2aHcAZmbWXk4EZmYV50RgZlZxTgRmZhXnRGBmVnFdd9E5SU8Cj9RMGgc81aZwWqkq5YTqlLUq5YTqlLWTy7ljRGxdb0bXJYL+JN2VdUW94aQq5YTqlLUq5YTqlLVby+mmITOzinMiMDOruOGQCC5udwAtUpVyQnXKWpVyQnXK2pXl7Po+AjMzG5rhUCMwM7MhcCIwM6u4jk0Ekg6R9ICkJZL+vc78TSRdmc7/H0kTa+Z9Ip3+gKSDWxr4IAy2rJImSnpB0r3pY3bLgy8gRzn/VtLdkl6W9J5+846T9GD6OK51UQ/OEMu6puY7vbZ1UReXo5wfk/RbSYsk3Sppx5p5w+07bVTWzv5OI6LjHsCGwO+BScDGwEJgj37LfBiYnT4/Grgyfb5HuvwmwE7pdjZsd5lKKutE4P52l6GJ5ZwI7AlcBrynZvpWwNL075bp8y3bXaYyyprO6213GZpYzunAa9LnJ9f8dofjd1q3rN3wnXZqjeBNwJKIWBoRLwE9wDv6LfMO4Dvp8x8CMyQpnd4TEasi4iFgSbq9TjWUsnaTAcsZEQ9HxCLglX7rHgzcHBFPR8RfgJuBQ1oR9CANpazdJE8550XE8+nLO4EJ6fPh+J1mlbXjdWoi2A6ovVHtsnRa3WUi4mVgBTA257qdZChlBdhJ0j2SbpP0N2UHOwRD+V6G43fayEhJd0m6U9I7mxpZcxUt54nADYNct92GUlbo8O/Udyjrbn8CdoiI5ZKmANdIekNEPNvuwGxIdoyIxyRNAn4m6b6I+H27gxoKSe8DpgIHtjuWsmWUtaO/006tETwGbF/zekI6re4ykjYCxgDLc67bSQZd1rT5azlARCwgacPctfSIB2co38tw/E4zRcRj6d+lwHxg72YG10S5yinp7cAZwJERsarIuh1kKGXt/O+03Z0U9R4kNZWlJJ29fR0zb+i3zCms24E6N33+BtbtLF5KZ3cWD6WsW/eVjaQT6zFgq3aXabDlrFl2Dut3Fj9E0qm4Zfq8I8vZhLJuCWySPh8HPEi/TslOeeT87e5NcoCyS7/pw+47bVDWjv9O2x5Agw/+MOD/0g/2jHTaOSSZFmAk8AOSzuBfA5Nq1j0jXe8B4NB2l6WssgLvBn4D3AvcDcxsd1mGWM59SdpeV5LU7n5Ts+770/IvAU5od1nKKivwFuC+dEdzH3Biu8syxHLeAjye/kbvBa4dxt9p3bJ2w3fqS0yYmVVcp/YRmJlZizgRmJlVnBOBmVnFORGYmVWcE4GZWcU5EZgVIOlhSeOGuoxZJ3EiMDOrOCcCswySrpG0QNJvJM3qN2+ipN9JulzSYkk/lPSamkU+kt5v4D5Jr0/XeZOkX6UXCbxD0m4tLZBZBicCs2zvj4gpJBcQO1XS2H7zdwMuiojdgWdJ7hvR56mI2Af4OvBv6bTfAX8TEXsDZwGfLTV6s5ycCMyynSppIcm15bcHduk3/9GI+GX6/HvAATXzrk7/LiC5CQ0kFwv8gaT7gfNIrotl1nZOBGZ1SJoGvB3YPyL2Au4hueZTrf7XZ6l93XflyTW8ern3TwPzImIyMLPO9szawonArL4xwF8i4vm0jX+/OsvsIGn/9Pl7gV/k2GbfpYuPb0qUZk3gRGBW343ARpIWA58naR7q7wHglHSZLUn6Axr5T+Bzku7BN4WyDuKrj5oNgqSJwPVpM49ZV3ONwMys4lwjMDOrONcIzMwqzonAzKzinAjMzCrOicDMrOKcCMzMKu7/AxNh7PmTObqEAAAAAElFTkSuQmCC", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "fig, ax = plt.subplots(nrows = 1, ncols = 1)\n", + "ax.plot(a_s[0], k_s2, \"s\", label = \"$k_2$ / $k_1$\")\n", + "ax.plot(a_s[0], k_s3, \"o\", label = \"$k_3$ / $k_1$\")\n", + "plt.legend()\n", + "plt.grid(True)\n", + "plt.xlabel(\"alpha\")\n", + "plt.ylabel(\"Анизотропия проницаемости\")\n", + "plt.title(\"Простая кубическая\")\n", + "plt.show()\n", + "\n", + "if savefig:\n", + " fig.savefig(\"anisotropy-simple.png\")" + ] + }, + { + "cell_type": "code", + "execution_count": 18, + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAtr0lEQVR4nO3de5wcVZn/8c+XJBAkEDBABAOGiyjID4QEFNE14aJcArguRnRhDYIRxI3uiiKKARFEXJQFESKKgIiGiMAGERRCgggiS4IJYsDFyGUCShIgZLjEJDy/P6oGOp2emdMzXd01M9/369Wv6bqdeqp7up6qc6pOKSIwM7OBa71WB2BmZq3lRGBmNsA5EZiZDXBOBGZmA5wTgZnZAOdE0A9IGtLqGMrGn4lZOieCPkjSaElXSVok6VngglbHVAaSJku6V9KTwDJJ27Y6JrO+wImgBSQ9KumAiuE35jv1byYsuxlwF/AAsGtEbBYRnyow3D5B0heBk4DjImLriNgkIh5vdVxmfcHgVgcw0EnaArgNuDkivpCwyGeAmRHRbdIYKCRtBHwB2D0inmh1PGZ9jc8IWkjSpsCvgXuBT+fj3iDpRUkjKubbU9KSvN57b2BEfgbxjKSZkraumDck7Zi/31bSS5J+nA9PkvTbinmrh98q6da83IclTayYtqGkb0l6TNJySb/Nx43O1zk4n+9Tkh7siF/SsZIWSlqRx/zJqs/gdElPSWrPY320i8/r8Lzs5yTNkbRzPultwArgIknPSnpE0ifyZd6Tl/1CZZw1yu5q26/Il92jYtw383EH5MNndHzO+fDgfProfHi4pB/l3+Njkk6TtF6t7yEf1yZpXMXwBEl/yLf9bkm7VUyrPsM8XtKc/P2NVdvfnr+mVSx7qqQ/5Z/d5ZKG5tPGSWqrKHdiXsbxFXHfJemi/H/iIUn7V8y/df7/+Uzld1Lxea3KY3lO0vWSNs6n7SDpdknLJC2VdLWy30q325sPv/obyIfPknRFxfDPJP0tj/k3kt5WMe0Nkn6dx9Sex3gG/ZwTQesMA24mOyv7eOR9fUTE34A5wMSKeY8BpkfEKuB1wH7AkcBWwGPA9E7W8TVgWcXwK3TynSs7qr4V+AmwJXAUcLGkXfJZzgPGAO8CXk92BP5KVRlHAScD74+IjvU+DUwANgGOBc6XtGc+/1uBU4HxETEMOKyT7UDSTsBPgc8CWwC/BG6UtH7+mWwLLM4/kyOBr0vaLyLuzMt+W82C07Yd4CGgYwc4JI/1752VWcN3gOHA9sB7gX/LP49u5Qnoh8AngRHA94CZkjbobtmIOKxq+zeNiGERcULFbP8KvB/YAdgJOK1GDEPI/p+eqpr0DuAvwObA6cB1kl6fT5sOtAFbU/GdVCx7TR7btsB2wMc6Vgecky+3M7ANcEZ321qHm4E3k33X84CrK6Z9FlgDbJXHdk0D11taTgStcwnQDowC9q2adiVwNICkQcBHgKsqpv8wIuZFxEqyHek+HUeeHfIjxn3ysjo8DuwsaVSNeCYAj0bE5RGxOiLuB34OfCg/cv048JmIWBwRayLi7nz9HQ4CLgMOjohXjyIj4qaI+Etk7iA7A3pP1bpTqig/DNwUEbfmCfE8YEOyxARZUjolIl6OiD8APyDb2abodNsr5pkJHChpQ7IkcBvwckrh+Xd4FHBqRKyIiEeBb5El+BSTge9FxO/zz/5KYCXwzsTlu3NRRDwREc8AZ5P9v1X7JPB74M9V458G/jsiVkXENcDDwKGStiH7v075TgaR7YuWAUTEI/n3vDIilgDfJkueDRERP8y/h5VkCWZ3ScMrZlmPAbZvHFAbWzIPke1QvgD8IN/BdPgfYBdJ2wEHAssj4t582kqyswAAIqKd7Af0xqryzwW+AqyqGHcH2RHOfEnPARdXTHsT8I78lPi5fPq/Am8gO9obSnbk15kfAI9S9YOVdLCke/LqgeeAQ/LyiIiHyI4+75TUDtzQRflbV233K8AT+XavBJ6JiBUV8z/Gup8JwNK8CmShpKMTtr3DKrJkcCRwXL691SZWLL+0YvzmwJDK+LuIr5Y3AZ+rim8bss+kww0V0y5MLLdDZbvKY1XlklfZfIHs/6na4o6z2arlt6b772RiHu8S4AXgxnx9IyVNl7RY0vPAj8n/Zyp0t73zKqafXLEtgyR9Q9Jf8rIfzSd1lP8t4EVgRb5s5Zl5v+VE0Dpn50dK3yf7IX6tY0JEvAzMIDsrOIa1zwYeJ9sxAK9Wa4wgqxbpsF8+bkblCvOj8hMiYkREbApUXm30BHBHRGxa8RoWESeS7dReJqs66MxHyI7az+4448irLn5OdvQ+Ml/nL8lO/TvMINuR7wp8oIvyn6zabpHtDBfnn8nrO+qYcx1VRdU2j4jNgBOBKyQN62bbK/2AbIc4IiLm1yh7RsfyrL3jWkqWSN5UMa6z+Gp5guz/pTK+10XETyvm+UDFuqckltthm6q4nqya/nmybXuMdb0x/y6ql3+S7r+TGXm8ryO7Cu5b+fivAwH8v4jYhOx3ULkO6H5796yYfl7F+I8CRwAHkFXVjc7HCyA/A7mT7OKNTan6DfVXTgTl8AlgsqS9K8b9CJgEHM7aieCnwLGS3p7vaL8O/D6vbuhwBvCFqiO17vwC2EnSMZKG5K+9JO2cH33/EPh23gA4SNI+VXXUd0bEH8mOzi7Nx60PbEB2xLda0sHA+6rWewnwX1Xx1zKDrMph/7y++nNkCeTuiFgM/BY4R9LQvFrsOLIjyc6sJvvxq6ttr1wgP4O5haz+OllErMnjP1vSxpLeBPxnN/FV+j5wgqR3KLORpEOrdrK9cZKkUXnd/pdZu158Y7K2jLM7WXZLYEr+mX2IrE7/l/nVW3eT9p28Qrbj36Jine3AcklvJEtEjbIx2f/NMrIE9PXKiXkV6ymsfZDU7zkRlEBELAKmApfnjZ9ExF1kP5B5lUdiEXF7Pu/PyRrudiCrf650f0TMqTOGFWQ76aPIjub+Rla91LGzP5nsqO1/gWfyabX+f74BbCXpY3mZU8h2gs+SHY3N7JhR0kfJql/+OyG+h8mODL9DdoR9GHBYRPwjn+VfyRocnwKuA6ZGxG01inpO0gqy5Do5ryvubtsr4/h8RPxPd/HW8O9k1R+LyJLWT8iSa4e9lF0p1KbsSp03AD+TNCoi7iM7WLiI7HN8hOwgoVF+QtZ2s4is+u+simmbABdGxLOdLPt7sobXpWTJ4sh47UKBj5AdcT8JXA+cXvWdfDivElwG7AJ8KR//VWBPYDlwE9n32Sg/IquiWgz8Cbinavr3gG90cvbTb6m+g0ZrJkm3Az+JiFr10dbPKbskclLC2VJv1vEocHwnSbO7ZSfly7670XFZc/mGspKStBfZUdERrY7FWmYu8FKrg7D+z4mghCRdSdZw+pmqqy5sAImIz7U6BhsYXDVkZjbAubHYzGyA63NVQ5tvvnmMHj261WG86oUXXmCjjTZqdRidKnt8UP4Yyx4fOMZGKHt80LsY586duzQitqg5MSL61GvMmDFRJrNnz251CF0qe3wR5Y+x7PFFOMZGKHt8Eb2LEbgvOtmvumrIzGyAcyIwMxvgnAjMzAY4JwIzswGu0KuG8tvXV5A96GF1RIytmj6OrMvlv+ajrouIM4uMyczM1taMy0fHR8TSLqbfGRETilr52LNuZWn7P9YZv/mw9bnvtAN7PK+ZWX/R76uGau3YOxtfz7xmZv1FoV1MSPorWbe5QfaovUurpo8j6065jayr2pMj4sEa5Uwme1wfI0eOHDN9emeP6F3XpFte6HTaFQetfWNGPfN2aG9vZ9iwYcnxNFvZ44Pyx9js+Kbc/gLP1zj22GR9uHC/vvl/COWPsezxQe9iHD9+/Nzq6vkORVcNvTsiFkvaErhV0kMR8ZuK6fOAN0VEu6RDyB5V+ObqQvIEcinA2LFjY9y4cekR3HJTp5PWKSdx3rWrkETWzfy6VUhlqGqaM2fOuttZMmWPsdnxPd/J/+Hz/6jxP5sr+2cI5Y+x7PFBcTEWmggie3IUEfG0pOuBvYHfVEx/vuL9LyVdLGnzbtoUWi61Cil1vjIkjEbpT9tiNlAUlgjyZ+muFxEr8vfvA86smucNwN8jIvLHNK5H9rSiAaU/tU20alucgMx6rsgzgpHA9flzrQeTPWnrFkknAETENOBI4ERJq8kewHFUNLjRYvNh63e6g+jNvFYu/SmZmjVbYYkgsufw7l5j/LSK9xeRPYe1MPUcDZb9yNFHvc2zzmed19v7smPrj/pcN9QDWSuPehu9s0vd0Ra+3k7KK+Ky49R1+8zUms2JoAdSf6j96Qfd6CTU6Ab3Rq+3CKnrLvtZhM+A+h8ngh6o/Gfv6nKu1B9FoxNGPdUajdafkl9/kbrjLuJsqZX/i5bOiaAEGv2DaOVlq636cQ/EBJS6k23V2VdRZVrjOREMYP3pRzoQjy770/fn6qbW6vd9DfUnnR3d9uej3lap57P299J7/Smp9UU+I+hDWnlk1Oiql1Y1uKeWl9oOVD1vI9Zt1mxOBJak0Umo0Q3uPVlvs/WXKg4ntP7HiaAf8g/VKjX67KuehOb/xb7BiaAfSj3a9o+0b2v0jruIM5Z6qtisdZwIBrD+UlUxUPWnnawPSlrLicDMWs4HJa3ly0fNzAY4JwIzswHOVUNm1mf4DuRi+IzAzPoM34FcDCcCM7MBzonAzGyAcyIwMxvgnAjMzAa4bq8akvQAEJWjgIiI3QqLysysBt+BXIyUy0fnA28DpgILig3HzKxzvkS0GN0mgog4WtKuwFnACmBqRPy18MjMzKwpUqqGXg88CXwc2Bf4maR7IuLTCcs+SpY81gCrI2Js1XQBFwCHAC8CkyJiXr0bYWZWqVU3nvXVG95Sqobm8lobgfK/h9SxjvERsbSTaQcDb85f7wAuyf+amfVYo288S93B99Ub3lKqhrYrcP1HAD+KiADukbSppK0i4qkC12lmVpe+uoNPlVI1NC8i9uxh+QH8WlIA34uIS6umvxF4omK4LR/nRGBmhVvnSP+Wm4DmVOWUqRoppWpI3c/SqXdHxGJJWwK3SnooIn5TbyGSJgOTAUaOHMmcOXN6EVJjtbe3lyqeamWPD8ofY9njA8dYj8oYujrST421p/P1ZN1FfYYpieAtkiovG02+jyAiFud/n5Z0PbA3UJkIFgPbVAyPysdVl3MpcCnA2LFjo0xPYir7k6HKHh+UP8ayxweOcR35kX0ta8XQqvnqnTdX1GeYkgj+ChxWb8GSNgLWi4gV+fv3AWdWzTYT+LSk6WSNxMvdPmBmvdWqG8/66g1vKYngHxHxWA/KHglcn10hymDgJxFxi6QTACJiGvBLsiuQHiG7fPTYHqzHzGwtja5jT93Bl/kS0a6kJIJ/70nBEbEI2L3G+GkV7wM4qSflm5k1S1/dwadK6XTuEEmbdgxI2kzSWcWFZGbWHJ1V2TSjKqeV666WckZwcER8qWMgIp6VdAhwWnFhmZkVr/JIv9kN7mU6y0hJBIMkbRARKwEkbQhsUGxYZmYGzbnXISURXA3MknR5PnwscGVD1m5mZl1qxl3NKV1MnCtpPnBAPuprEfGrhkVgZmYtlXJGAHA/MISsy4j7iwvHzMyardurhiRNBO4FjgQmAr+XdGTRgZmZWXOknBF8GdgrIp4GkLQFcBtwbZGBmZlZc6TcR7BeRxLILUtczszMeqkZ9xuknBHcIulXwE/z4Q+TdQ1hZmYFa8a9DilXDX1e0r+QPaYS4NKIuL7hkZiZWUskXTUUET8Hfl5wLGZm1gIpTyh7BWivHEXWX9wmhUVlZmZNk3JG8F/AHsAtwI+rGo7NzKyPS2kjOEXSesDBwHclPRsRk4sPzczMmiH1MtCdgPcCw8luLjMzs34ipY3gbrJ2gcuB6fm4PSNiXsGxmZlZEyQ9qpKsj6GPAEeRNxYD+xUYl5mZNUlKG8G4JsRhZmYtktLp3EhJl0m6OR/eRdJxxYdmZmbNkNJYfAXwK2DrfPjPwGcLisfMzJosJRFsHhEzgFcAImI1sKbQqMzMrGlSEsELkkaQNRAj6Z3A8kKjMjOzpkm5aug/gZnADpLuArYge0iNmZn1AylXDc2T9F7gLWSXjj4cEatSVyBpEHAfsDgiJlRNm0TWhcXifNRFEfGD1LLNzKz3Up9Z/FZgF2AosKckIuJHict+BlgIdNZJ3TUR8enEsszMrMFSLh89HfhO/hoPfBM4PKVwSaOAQwEf5ZuZlZQiousZpAeA3YH7I2J3SSPJeiE9sMsFs2WvBc4BNgZO7qRq6BxgCdllqf8REU/UKGcyMBlg5MiRY6ZPn56wac3R3t7OsGHDWh1Gp8oeH5Q/xrLHB46xEcoeH/QuxvHjx8+NiLE1J0ZEly/g3vzvXLLqHQEPJSw3Abg4fz8O+EWNeUYAG+TvPwnc3l25Y8aMiTKZPXt2q0PoUtnjiyh/jGWPL8IxNkLZ44voXYzAfdHJfjXl8tH7JG0KfD9PBvOA3yUsty9wuKRHyTqr20/Sj6uS0LKIWJkP/gAYk1CumZk1UMpVQ5/K306TdAuwSUQsSFjuVOBUAEnjyKqGjq6cR9JWEfFUPng4WaOymZk1UepVQwBExKMAkr5I9oyCn0fETfWUIelMslOUmcAUSYcDq4FngEn1lGVmZr2X8jyCv5LfVdwxChhJdknpcykriYg5wJz8/dSK8a+eNZiZWWuknBFUtzILuDEiHi8gHjMza7KUNoJl1eMkJd9ZbGZm5ZZSNfT6nixnZmZ9Q8oOfS5ZG4EqxnV9F5qZmfUZKVVD2zUjEDMza42UqqEP1hofEdc1PhwzM2u2lKqha4A/kVURdVQPBeBEYGbWD6Qkgl2BrwHDgK9ExMPFhmRmZs2U0kbwMDBR0hjg25KeBM6IiMXdLGpmZn1AShvBd3jtKqFFwHuB/wNeV2BcZmbJVq1aRVtbGy+//HKPyxg+fDgLF5a7u7OUGIcOHcqoUaMYMmRIcrkpVUP3dTNsZtZSbW1tbLzxxowePRpJ3S9Qw4oVK9h4440bHFljdRdjRLBs2TLa2trYbrv0Cz5TqoauTC7NzKwFXn755V4lgf5CEiNGjGDJkiV1LdfTTuciIravL0Qzs+IM9CTQoSefQ2qncwJuJ3tmsZmZ9SPJnc5JWl2rAzozM+vb6ul0bpCkzchvKouIZ4oMzMzMmqPeTufm5eMCcBuBmfU5Y8+6laXt/1hn/IiNhjD3K+/rVdmzZs3iiiuu4KqrrqpruRNOOIFjjjmGfffdtyHl1cudzpnZgFIrCQAse6H3j1mZP38+e+yxR93L3XPPPXz3u99tWHn1Wq+7GSQNkTRF0rX569OS0u9UMDMbIDp23CtXrmTSpEl86UtfIqLrXvsXLlzITjvtxKBBg7ot76tf/Wq35fVEStXQJcAQ4OJ8+Jh83PENj8bMrA9bsGABW265Je9///s5/vjjOfroo7td5uabb+aggw5KKu+II44o5DLZlESwV0TsXjF8u6T5DY/EzKwPW7VqFYsWLeIjH/kI3/ve99hnn3244YYbuOmmm3j++ec57rjjeN/71m2D+NWvfsXll1+eVN6KFStYtGgRZ599NsuXL+faa69tSOzdVg0BayTt0DEgaXtgTUPWbmbWTyxcuJC99tqLwYMHv1rN84EPfIDvf//7TJs2jWuuuWadZV588UWee+45tt5666TyALbffnsuu+yyhsaekgg+D8yWNEfSHWQ3ln2uoVGYmTXJ5sPWrzl+xEa9a/qcP38+73rXu5g+fTrHHnssf//731+ddtZZZ3HSSSets8zs2bMZP772fbpdlddoKVcNzZL0ZuAt+aiHI2Jl6gokDSLrqG5xREyomrYB8CNgDLAM+HBEPJpatplZve477cCa41esWNGrcufPn8/ee+/NTjvtxLnnnsvEiRO57bbbOO200zj44IPZc88911nm5ptv5sgjj0wu7/rrr+9VjJ1JuWpoKHAScAZwOnBiPi7VZ4DO+k09Dng2InYEzgfOraNcM7PSOO+885g4cSIAEyZM4I477uCSSy7htttu49prr2XatGnrLHP33Xevc+9AV+UNGTKEZcuWccIJJ3D//fdzzjnnNCT2lMbiHwErgO/kwx8FrgI+1N2CkkYBhwJnA/9ZY5YjyBIMwLXARZIURVwfZWbWZFOmTGHKlCmdTp83b16n0zozYsSImkmlN9TdPlfSnyJil+7GdbLstcA5wMbAyTWqhv4IHBQRbfnwX4B3RMTSqvkmA5MBRo4cOWb69OndbliztLe3M2zYsFaH0amyxwflj7Hs8YFjHD58ODvuuGOvylizZk3Na/nLJDXGRx55hOXLl681bvz48XMjYmyt+VPOCOZJemdE3AMg6R0kPJxG0gTg6YiYK2lcwno6FRGXApcCjB07NsaN61VxDTVnzhzKFE+1sscH5Y+x7PGBY1y4cGGvHyrTHx5M02Ho0KF13ZGckgjGAHdLejwf3hZ4WNIDZM8l2K2T5fYFDpd0CDAU2ETSjyOi8g6LxcA2QJukwcBwskZjMzNrkpREUPuWt25ExKnAqQD5GcHJVUkAYCbwMeB3wJHA7W4fMDNrrpREUHPHHBGP1xrfHUlnAvdFxEzgMuAqSY8AzwBH9aRMMzPruZRE8BDwCPlzCHIBdFYltI6ImAPMyd9PrRj/MglXH5mZWXFSEsHDEVF8P6hmZtYSKV1MuM7ezKwfS0kEm0r6YPWr8MjMzPqYWbNmccwxx9S93AknnMBdd93VsPLqlVI1dAdwWNW4AK5rfDhmZk2wYAbMOhOWt8HwUbD/VNju4F4X22+fUBYRx9Z4fbzwyMzMirBgBtw4BZY/AUT298YpDF7Y+w7d+uoTylKqhszM+o9ZZ8Kql9Yet+olNrjzG70uuvKJYgcccABf//rXu32iWOoTyg444ABOP/30lj2hzMys/1jeVnO0VjzZq2JrPVFs4cKFXHDBBSxdupT999+fE088cZ3l+soTyszM+o/ho2qOjo3XfUpYPWo9UWznnXdm2rRpzJgxo2ZjcJ95QpmkqbVeDY3CzKxZ9p8KQzZce9yQDVn5ni/2qtjOnig2c+ZMDj30UA455JB1linLE8pSzggmA+3AC8An8r8vFBaRmVmRdpsIh10Iw7cBlP097EJW7/zPvSp2/vz57Lrrrms9UWzVqlUcfvjh3HzzzVx99dXrLNNV+0Bn5RUhpY1gSUR8G0DSCcDVEfG3QqIxM2uG3SZmr0q9fFTleeed9+r7CRMmMGHCBObMmcN1113HypUra54R3H333Zx//vnJ5a1YsYJly5bx5S9/+dUnlJ166qm9ihvSEoEk/TOwCbASuF7SqXn/QWZm1olx48Z1+QyGsjyhLCURfAI4BVhD1kHc34HLyTuRMzOzvq3bRBAR/0v2rIBKRxQTjpmZNVu3iUBSrYfO09FuYGZmfVtK1dBXgMeA3t9/bWZWkIgo5K7bvqYnXVCkJIIdyB45uT9wZkTcVvdazMwKNHToUJYtW8aIESMGdDKICJYtW8bQoUPrWi6ljeAZ4POStgZOl3Qy8JW87cDMrOVGjRpFW1sbS5Ys6XEZL7/8ct070GZLiXHo0KGMGlX77unOpLQR3MhrD6cRsC1wD7BuV3lmZi0wZMgQtttuu16VMWfOnKZ0+dwbRcWYUjV0XvezmJlZX5VSNXRHMwIxM7PWSKkaWg28WDkKiIjYpLCozMysaVKqhh6IiHJXnJmZWY+l9D7ao+eiSRoq6V5J8yU9KOmrNeaZJGmJpD/kr+N7si4zM+u5lDOC10nag6xK6FUR0V1vSSuB/SKiXdIQ4LeSbo6Ie6rmuyYiPp0espmZNVJKIngKqO5OIoD9uloostvb2vPBIfmr8U9dNjOzXlFPbkdOLlwaBMwFdgS+GxGnVE2fBJwDLAH+DPxHRDxRo5zJZA/IYeTIkWOmT59eWMz1am9vZ9iwYa0Oo1Nljw/KH2PZ4wPH2Ahljw96F+P48ePnRsTYmhMjossXMJzsjOC+/PUtYHh3y1WVsSkwG9i1avwIYIP8/SeB27sra8yYMVEms2fPbnUIXSp7fBHlj7Hs8UU4xkYoe3wRvYsRuC862a+mNBb/EFgBTMxfz5M9jyBZRDyXJ4KDqsYvi4iV+eAPgDH1lGtmZr2Xkgh2iIjTI2JR/voqsH13C0naQtKm+fsNgQOBh6rm2api8HBgYXLkZmbWECmNxS9JendE/BZA0r7ASwnLbQVcmbcTrAfMiIhfSDqT7BRlJjBF0uHAauAZYFJPNsLMzHouJRGcSLZDH052CekzwMe6WygiFgDr3IgWEVMr3p9K1sW1mZm1SEpfQ38Adpe0ST78fNFBmZlZ83TbRiBphKQLyR5WP1vSBZJGFB6ZmZk1RUpj8XSy6/z/hewh9kuAa4oMyszMmieljWCriPhaxfBZkj5cVEBmZtZcKWcEv5Z0lKT18tdE4FdFB2ZmZs2Rkgg+AfyErBO5lWRVRZ+UtEKSG47NzPq4lKuGNm5GIGZm1hopVw39ohmBmJlZa6RUDW1deBRmZtYyKVcNbS9pZvXIiDi8gHjMzKzJUhLBErKup83MrB9KSQTtEXFH4ZGYmVlLpLQRnFN4FGZm1jIpZwQ3SJoC/FM+fAcwLSJWFReWmZk1S0oiuJjswfMX58PHAJcAxxcVlJmZNU9KItgrInavGL5d0vyiAjIzs+ZKaSNYI2mHjgFJ2wNrigvJzMyaKeWM4PNkzyFYRPaEsjcBxxYalZmZNU1KX0OzJL0ZeEs+6uGIWFlsWGZm1iwpfQ3tBEwle2D9+4Bp+TgzM+sHUtoIrsr/3goMAuYDlxcWkZmZNVVKIhgcEV8GXomIcyPiv4FhxYZlZmbNkpIIXsn/nlQxTgXEYmZmLZCSCD4OEBE3A0jahKzNoEuShkq6V9J8SQ9K+mqNeTaQdI2kRyT9XtLoOuM3M7Ne6jYRRMQDVcPPR8QNCWWvBPbLb0Z7O3CQpHdWzXMc8GxE7AicD5ybErSZmTVOylVDKyQ9n79WpD6rODLt+eCQ/BVVsx0BXJm/vxbYX5KrnczMmkgR1fvmqhmks4D9gLMj4qa6CpcGAXOBHYHvRsQpVdP/CBwUEW358F+Ad0TE0qr5JgOTAUaOHDlm+vTp9YRRqPb2doYNK2/bednjg/LHWPb4wDE2Qtnjg97FOH78+LkRMbbmxIjo9gVsAVwIzAL2TVmmavlNgdnArlXj/wiMqhj+C7B5V2WNGTMmymT27NmtDqFLZY8vovwxlj2+CMfYCGWPL6J3MQL3RSf71ZSqoT2BbYAryHogvbjeB9pHxHN5IjioatLivGwkDQaGA8vqKdvMzHonpa+h6sdUPgNs1N1CkrYAVkXEc5I2BA5k3cbgmcDHgN8BRwK355nLzMyaJKWvofE9LHsr4Mq8nWA9YEZE/ELSmWSnKDOBy4CrJD1ClmCO6uG6zMysh7pNBJJq3jMQEWd2tVxELAD2qDF+asX7l4EPdR+mmZkVJeWGshfy1ycq3r9QZFBmZtY8KVVD3wKQdHTHezMz6z9Szgg6uBHXzKwfSmkjuJEsCWwvaWbH+Ig4vMjAzAa0BTNg1pmwvA2Gj4L9p8JuEzud773L2+D+LuYz60LK5aPn5X9dLWTWDAtmwI1TYNVL2fDyJ7JhWHsnXzGfuprPrBspnc7dUfkie3C9/8vMemLBDDh/Vzhj0+zvghnrzjPrzNeSQIdVL2XjezJfEVK2w/qMlDMCJO0BfJTsUs+/Aj8vMiizfin1SH95W+3lq8enztdoqdthfUanZwSSdpJ0uqSHgO8Aj5N1Ujc+Ii5qWoRm/UXqEfzwUbWXrx6fOh809gi+lWciVoiuqoYeIut1dEJEvDsivkNWLWRm1VJ2tKlH8PtPhSEbrj1uyIbZ+J7M13EEv/wJIF47gu9pMmjVmYgVpqtE8EHgKWC2pO9L2h8/otJsXak72tQj+N0mwmEXwvBtAGV/D7tw3WqXivmiq/nqOYJPSWj1nInUw+0OLdNpG0FkTyG7QdJGZA+Q+SywpaRLgOsj4tdNidD6ntRLH/uLrna0ldu9/9S169ah9hE8ZMulfGb5fHfMmcO4ceNqz5N6BJ9a91/PdqRyu0NLpVw19EJE/CQiDgNGAfcDp3SzmA1Uja6G6AtSd7SpR/qNlnoEn3rmUMR2uN2hpZKuGuoQEc8Cl+avvqOeI9RGH832p6PjlG1JPTruC1Jv1ho+Kk98NcZXSz3Sb6TUI/h66v4bvR1ud2iperqY6JvqOUKtd97zd+W9cz7QeX1mEUfHqfWoKfHVW17KttTzg06NMVWdn02X81Vsr7r77lIbbVsl9Qi+4Lr/Lr/notZtSfp/IqjnlDN13tSdRKMb6VJ3xqnx1ZOoGn3pYz072oI+my7nq+e7a1WVTz12mwj/8Uc447nsb63Yikhoqd9z2ZNpP9f/E0E9R6ip86buJOptpGvUzqmIO1MbfeljD5JuUz+beqsqUna0ZdfKuv961u2rixqurjaCPqme+tvUeVN3EqnlpdarN/qO03p2dqnb0hFvd20JjUi6RX429fzf9CetrPtPWbevLipE/z8jqOeUM3Xe1OqP1PLq2TmlrLeIO1Pr+RxTjo5T192qz8ZVFY3R6Lp/X11UiP6fCOo55UydN3Un0ehGukbfcVrvzr2R1QatSro9+O66vFnLutbohOqriwrR/6uGoL7T3ZR5K6o/Ynkb6uqy0JTyUi/vS612SY0vtbx6tiVVaowFfjbdbnPKzVrWtXp+Kyn6QpVdH7xkXBF968FjY8eOjfvuu6/VYbxqTqN2EgX98zQsvgJ1G2OLf1j94jMsgYbEWN1GANmBQQPO1soeH/QuRklzI2JsrWkD44ygL2jFjUZ9hT8b61DPGV0rDiD66A2VTgRm1re06uqilMTSR9swCmsslrSNpNmS/iTpQUmfqTHPOEnLJf0hf/mSDDPrvR7czNmQXgLqvUqqJPdEFHnV0GrgcxGxC/BO4CRJu9SY786IeHv+8jVgZtZ7PbiZsyG9BNRzlVSJOmgsLBFExFMRMS9/vwJYCLyxqPWZmb2q0T2uFtHDbInuiWjKVUOSRgO/AXaNiOcrxo8je/5xG/AkcHJEPFhj+cnAZICRI0eOmT59euExp2pvb2fYsGGtDqNTZY8Pyh9j2eMDx1hty7/fwVse/i6DXln56rg1623Aw285iadHvvfVce+d84HsTKBKIO4Yd8Orw+/83fEMXblknfle3mAL7tnnBz2KMXXdlXrzGY4fP77Tq4YKTwSShgF3AGdHxHVV0zYBXomIdkmHABdExJu7Kq/fXj5akLLHB+WPsezxgWOsKaVx9/xdO7kvYZvsrvjKshp9WWjquisUdflooXcWSxpCdsR/dXUSAIiI5yOiPX//S2CIpM2LjMnMBohG9rhaRId8JerGpLDLRyUJuAxYGBHf7mSeNwB/j4iQtDdZYlpWVExmZmtpdC8BPVx36j0R3T4kqYeKvI9gX+AY4AFJf8jHfQnYFiAipgFHAidKWg28BBwVfe1WZzPr21rZlUid90QICulxtbBEEBG/hSzuLua5CLioqBjMzPq8Jtyt3P97HzUz68uacLeyE4GZWZk14XnOTgRmZmXWhKuLnAjMzMqsCQ9Jcu+jZmZlV/CVTT4jMDMb4JwIzMwGOCcCM7MBzonAzGyAcyIwMxvgmvI8gkaStAR4rNVxVNgcWNrqILpQ9vig/DGWPT5wjI1Q9vigdzG+KSK2qDWhzyWCspF0X2d9fJdB2eOD8sdY9vjAMTZC2eOD4mJ01ZCZ2QDnRGBmNsA5EfTepa0OoBtljw/KH2PZ4wPH2Ahljw8KitFtBGZmA5zPCMzMBjgnAjOzAc6JoBOSDpL0sKRHJH2xxvQNJF2TT/+9pNH5+AMlzZX0QP53v7LFWDF9W0ntkk4uW3ySdpP0O0kP5p/l0DLFKGmIpCvz2BZKOrWI+BJj/CdJ8yStlnRk1bSPSfq//PWxMsUn6e0V3/ECSR8uIr7exFgxfRNJbZIKebRuL7/jbSX9Ov8//FP17zxJRPhV9QIGAX8BtgfWB+YDu1TN8ylgWv7+KOCa/P0ewNb5+12BxWWLsWL6tcDPgJPLFB9Z9+gLgN3z4RHAoJLF+FFgev7+dcCjwOgWxTga2A34EXBkxfjXA4vyv5vl7zcrUXw7AW/O328NPAVsWqbPsGL6BcBPgIvKFh8wBzgwfz8MeF29MfiMoLa9gUciYlFE/AOYDhxRNc8RwJX5+2uB/SUpIu6PiCfz8Q8CG0raoEwxAkj6APDXPMYi9Ca+9wELImI+QEQsi4g1JYsxgI0kDQY2BP4BPN+KGCPi0YhYALxStez7gVsj4pmIeBa4FTioLPFFxJ8j4v/y908CTwM173xtVYwAksYAI4FfFxBbr+KTtAswOCJuzedrj4gX6w3AiaC2NwJPVAy35eNqzhMRq4HlZEeulf4FmBcRK8sUo6RhwCnAVwuIq9fxkR0phqRf5afDXyhhjNcCL5AdxT4OnBcRz7QoxiKWTdWQdUjam+xo+C8NiqtSj2OUtB7wLaCQ6tNcbz7DnYDnJF0n6X5J/yVpUL0B+AllBZH0NuBcsqPbsjkDOD8i2vMThLIZDLwb2At4EZglaW5EzGptWGvZG1hDVqWxGXCnpNsiYlFrw+p7JG0FXAV8LCLWOSJvsU8Bv4yIthL/Vt5DViX9OHANMAm4rJ5CfEZQ22Jgm4rhUfm4mvPk1QPDgWX58CjgeuDfIqKII5zexvgO4JuSHgU+C3xJ0qdLFF8b8JuIWJqf5v4S2LPB8fU2xo8Ct0TEqoh4GrgLKKKfmpQYi1g2Va/WIWkT4CbgyxFxT4Nj69CbGPcBPp3/Vs4D/k3SNxobXq/iawP+kFcrrQZuoCe/lUY3fPSHF1mWXQRsx2uNN2+rmuck1m5EnJG/3zSf/4NljbFqnjMoprG4N5/hZsA8skbYwcBtwKEli/EU4PL8/UbAn4DdWhFjxbxXsG5j8V/zz3Oz/P3rSxTf+sAs4LON/twaFWPVtEkU01jcm89wUD7/Fvnw5cBJdcdQ5BfQl1/AIcCfyeosv5yPOxM4PH8/lOyKm0eAe4Ht8/GnkdUd/6HitWWZYqwq4wwKSAS9jQ84mqwh+4/AN0v4PQ/Lxz9IlgQ+38IY9yI7MnyB7GzlwYplP57H/ghwbJniy7/jVVW/lbeXKcaqMiZRQCJowHd8INlVdg+QJYr1612/u5gwMxvg3EZgZjbAORGYmQ1wTgRmZgOcE4GZ2QDnRGBmNsA5EZjVQdKjkjbv7TxmZeJEYGY2wDkRmHVC0g3KninxoKTJVdNGS3pI0tV5P/DXSnpdxSz/nneY94Ckt+bL7J33v3+/pLslvaWpG2TWCScCs859PCLGkPUhNEVSde+ybwEujoidybqg/lTFtKURsSdwCa/1XPkQ8J6I2AOYCny90OjNEjkRmHVuiqT5wD1knYK9uWr6ExFxV/7+x2Q9pna4Lv87l+yhIpB1WPczSX8EzgfeVkTQZvVyIjCrQdI44ABgn4jYHbifrN+hStX9s1QOdzyDYg2vdff+NWB2ROwKHFajPLOWcCIwq2048GxEvJjX8b+zxjzbStonf/9R4LcJZXZ0LzypIVGaNYATgVlttwCDJS0EvkFWPVTtYeCkfJ7NyNoDuvJN4BxJ9+OHQlmJuPdRsx6QNBr4RV7NY9an+YzAzGyA8xmBmdkA5zMCM7MBzonAzGyAcyIwMxvgnAjMzAY4JwIzswHu/wOw0D8wsm1jIQAAAABJRU5ErkJggg==", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "fig, ax = plt.subplots(nrows = 1, ncols = 1)\n", + "ax.plot(a_b[0], k_b2, \"s\", label = \"$k_2$ / $k_1$\")\n", + "ax.plot(a_b[0], k_b3, \"o\", label = \"$k_3$ / $k_1$\")\n", + "plt.legend()\n", + "plt.grid(True)\n", + "plt.xlabel(\"alpha\")\n", + "plt.ylabel(\"Анизотропия проницаемости\")\n", + "plt.title(\"Кубическая объемноцентрированная\")\n", + "\n", + "if savefig:\n", + " fig.savefig(\"anisotropy-bodycentered.png\")" + ] + }, + { + "cell_type": "code", + "execution_count": 19, + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAp8ElEQVR4nO3deZgcZbn+8e9NSAiQBUwgGhIMKCARUSBs6jmERQRlO4oR+IEGxYDKiRuoCAZEFlEQ4YgiIIuIJyACJx5AZBtUECSAQSHEEyOQARQSSMgECCE8vz+qBorOLNU9XdXTPffnuvqarre25+nuqbfqfWtRRGBmZgPXGo0OwMzMGssVgZnZAOeKwMxsgHNFYGY2wLkiMDMb4FwRDBCSBjc6BjPrn1wRtChJEyRdLmmBpOeAcxodk5n1T64I+ilJj0raIzO8UbpR/26OedcH7gT+AmwVEetHxOcKDNfMmtiajQ7AeidpA+AW4MaI+GqOWb4AzIqIXisNMzMfEfRzktYDfgv8CTg6LXuzpBckjcpMt62kZ9K+gB2AUekRxLOSZkkam5k2JL09fb+xpBcl/TwdnirpD5lpK4ffIenmdLnzJE3JjFtb0lmSHpO0VNIf0rIJ6TrXTKf7nKSHOuOXdLikuZKWpTEfWfEZnCjpKUkdaayPdvNZda6nI/M6Ix3XJul0SX+S9Lyk/5H0psy8v5T0zzTu30l6Z2bcpZJOyQy/XVJkhkdK+mka4xOSTpE0qKvPLy1rlzQ5M7yPpD9LWiLpLklbZ8ZVHhkeIaktff/rNMflFXmfn5n3OEkPS3pO0iWShqbjJktqzyx3SrqMIzJx3ynph+ln8oik3TPTj01/V89Kmi/pM5lxJ0lamcayRNK1koan494m6TZJiyUtknSFkt94r/mmw6/9dtPhUyRdmvN7fLOk36YxdaQxnoS5IujnhgE3khy5fSrS+4FExD+BNmBKZtrDgJkRsRJYB9gNOBB4C/AYMLObdXwbWJwZfpVufheS1gVuBn4BbAgcBPxI0sR0kjOB7YD3Am8CvpouL7uMg4BjgA9GROd6nwb2AUYAhwNnS9o2nf4dwHHArhExDNi3mzyy1ouIYenra5nyTwCfIvlMXgHOzYy7Edgszet+4Ioc6+l0abq8twPbAHsCR+SZUdI2wMXAkcAo4CfALElr9TZvROybfiadG7vOvI/KTPb/gA8CbwM2B07oIobBJL+DpypG7Qj8HRgNnAhck6k8ZwLtwFiS39lpknbLzHtlGtvGwCbAJztXB5yezrclMB44qbdcq9DT9/hFYBXwljS2K+u43qbmiqB/+zHQAYwD3lcx7jLgUIB07/Ng4PLM+Isj4v6IWEGyId1Z0oTsAtI9z53TZXV6HNhS0rgu4tkHeDQiLomIVyLiAeBXwMckrUGykf1CRDwREasi4q50/Z32An4K7B0Rr+2NRsT1EfH3SNxBcgT0bxXrrkcz5uUR8deIWA58E5jSueceERdHxLI03pOAd0sa2dsCJY0BPgR8MSKWR8TTwNkklWQe04CfRMQ96Wd2GbAC2Knq7Lr2w4hYGBHPAqeS/E4qHQncA/ytovxp4AcRsTIirgTmAR+WNJ7k9/i1iHgpIv4MXERS0VYaRLKdWQwQEfMj4uaIWBERzwDfB3bpc5apHN/jGni7txp/IP3bIyR7wF8FLpK0dmbc/wATJW0CfABYGhF/SsetIDkKACAiOkj+ETeqWP4ZJBvElZmyO0j2lOZIWgL8KDPurcCO6aH1knT8/wPeTLLXOJRkD7I7FwGPUvGPL2lvSXenzQxLSDaso9PYHyHZi/29pA7guh6W35uFmfePAYOB0ZIGSfqOpL9Lej6Nkc4YUsdkcr4/U/7WdDlPZcb/hGSPNI+3Al+p+EzHk+wxd7ouM+7cLpbRk8qcs8slbbL5KsnvoNITnUehFfOPBZ6NiGUV47K/rylpvM8Ay4Ffp+sbI2lm2oT2PPBz3vg5Q+/53p8Zf0wml96+x7OAF4Bl6bzZI+oBzRVB/3Zqusd1Ick/9Lc7R0TES8BVJEcFh/HGo4HHSTYwwGtNOqOAJzLT7JaWXZVdYbpXflREjIqI9YDs2UYLgTsiYr3Ma1hEfBZYBLxE0gTRnYOBjwOndh5xpE0gvyJpVhqTrvMGkiaETleRVG5bAQf0sPzejM+835ikAlwEHALsD+wBjAQmpNNkYzizM2dg20z5wjS20ZnPZEREvJN8FpJ8z9nPdJ2I+O/MNAdk1j0953I7Veb8ZMX4Y4GrIuIxVreRpOxn0Dn/k8CbOtv9M+Oyv6+r0njXITl77ay0/DQggHdFxAiS3292HdB7vttmxp+ZKe/xe0yPQH5PctLFelT89gcyVwTN4zPANEk7ZMp+BkwF9uONFcF/A4dLek+6oT0NuCciHs1McxLw1Yo9vt78L7C5pMMkDU5f20vaMiJeJWnr/n7akThI0s4Vbd2/j4i/kuzlXZCWDQHWItlzfEXS3iRt7Fk/Br5XEX8tDpU0UdI6wMnA1RGxChhOsjFfTLLhOi3vAiPiKZKmrLMkjZC0Rtohmre540LgKEk7KrGupA9XbGT74vOSxqVt+8fzxnbx4SR9Mqd2M++GwPT0e/4YSZv+DRGxELgLOF3S0LSJ8dMke/eVXiXZ8G+QWWcHsFTSRiQVUb30+D2mTaNf4407N4YrgqYREQuAGcAlkoakZXeS/KPdn92ji4jb0ml/RdIB+DZWb7N+ICLaqoxhGclG+iCSvcJ/kjQvdW7sjyHZ+7sXeDYd19Vv7DvAWyR9Ml3mdJK9s+dI9upmdU4o6RCSpqcfVBNrNy4n6dj9J0kzVufe5s9ImjaeAB4G7q5yuZ8gqdAeJsnhapIO6U7bKzlTqF3JmTpvBn4paVxEzCap5H+YzjufpHKvl1+QVFQLSJrtTsmMGwGcGxHPdTPvPSQdr4tIKosDMx38B5PscT8JXAucGBG3ZOb9eNqUtxiYCHwjLf8WyRHVUuB64Jq+JFeht+/xJ8B3ujn6GdBU3Q6h9TeSbgN+EREXNTqW/kzJKYg/7y+fUxrP1Doc5fS0jkeBIyo20HnnnZrO+/56x2X9jy8oa2KStifZu9q/0bFY1e4DXmx0EGbgiqBpSbqMpOP0CxVnb1gTiIivNDoGs05uGjIzG+DcWWxmNsA1XdPQ6NGjY8KECY0OI5fly5ez7rrrNjqMQrRybtDa+Tm35tWX/O67775FEbFBV+OariKYMGECs2fPbnQYubS1tTF58uRGh1GIVs4NWjs/59a8+pKfpG5Pm3XTkJnZAOeKwMxsgHNFYGY2wLkiMDMb4FwRmJkNcE131pCZWTOadMrNLOp4ebXy0cOGMPuEDzQgotf5iMDMrARdVQI9lZfJFYGZ2QDnisDMbIBzRWBmNsC5IjAzG+BcEZiZlWD0sCFVlZfJp4+amZWg0aeI9sRHBGZmA5wrAjOzAc5NQ2Zm/dhqVyT/5nqgvlck+4jAzKwfK+OKZFcEZmYDnCsCM7MBrtA+Akl7AecAg4CLIuI7FePfClwMbAA8CxwaEe1FxmRm1hf9+S6itSrsiEDSIOA8YG9gInCwpIkVk50J/CwitgZOBk4vKh4zs3roz3cRrVWRRwQ7APMjYgGApJnA/sDDmWkmAl9O398OXFdgPGZWpVr3fmuZr+w97WbZsx89bEi3cdZLkRXBRsDCzHA7sGPFNHOAj5A0H/0HMFzSqIhYnJ1I0jRgGsCYMWNoa2srKua66ujoKDzW6bct5/kudkRGDIFzd1u3sPWWkVsjtVp+q/1O0lMQe/ud9LT329PnU8t8ta6r7Nx6UsRv5sz3DwYGA8nvctiwYXVfX6OvIzgG+KGkqcDvgCeAVZUTRcQFwAUAkyZNismTJ5cYYnXeuJchYDlQzB4UwPPpj3618pehu8+pPnt5+XNrRm1tbd1+fs2olt8J8NpGtSt1n6/GdTVFbnVS1O+yyIrgCWB8ZnhcWvaaiHiS5IgAScOAj0bEkgJjKlyt7Ydltjs2Q4zQPIfuZfJnYkUosiK4F9hM0iYkFcBBwCHZCSSNBp6NiFeB40jOIDIDaq946tauXcAVnH3Vih2VzaaMNvuyFVYRRMQrko4GbiI5ffTiiHhI0snA7IiYBUwGTpcUJE1Dny8qHhs4muWIx5pTf9kpqKdC+wgi4gbghoqyGZn3VwNXFxmDNV4rN2eUeVZNI9S691vLfGXvabfinn2tGt1Z3G+1+j9qmVp5T7tZjj5q/Z3U+luvZb5a11V2bq2o5SuCWjfotf6jlrkHBbX9mMuO0Rov+ztptTOiWjm3srR8RVD2nletP8oy907qsZdXxj+cK57V+TOxIrR8RWDNq+ymgmbYyLo5w4rgisBaTrMc8Zj1F71WBJL+AkS2CIj0RnFmvWqGPe1atfLRhw0ceY4I5gDvBGYADxYbTv/hf9T6aeXmjDLPqjErSq8VQUQcKmkr4BRgGTAjIv5ReGR14lPLzMx6lqdp6E3Ak8CngPcBv5R0d0QcXXRw9eANuplZz/I0Dd3H630ESv9+qJhwzMysbHmahjYpIxAzM2uMXh9VKen+MgIxM7PGyPPMYvU+iZmZNas8fQRbSMqeNurrCMzMWkieiuAfwL5FB2JmZo2RpyJ4OSIeKzwSMzNriDx9BP9ZeBRmZtYweSqCD0lar3NA0vqSTikuJDMzK1OeimDviFjSORARz+ELyszMWkaeimCQpLU6ByStDazVw/RmZtZE8nQWXwHcKumSdPhw4LLiQjIzszLlucXEGZLmAHukRd+OiJuKDcvMzMqS9wllDwCDSW4+90Bx4ZiZWdny3GtoCvAn4EBgCnCPpAOLDszMzMqR54jgeGD7iHgaQNIGwC3A1UUGZmZm5chz1tAanZVAanHO+czMrAnk2aD/RtJNkqZKmgpcD9yQZ+GS9pI0T9J8SV/vYvzGkm6X9ICkByX5+gQzs5LlOWvoWEkfJXlMJcAFEXFtb/NJGgScB3wAaAfulTQrIh7OTHYCcFVE/FjSRJIKZkKVOZiZWR/kOmsoIn4F/KrKZe8AzI+IBQCSZgL7A9mKIIAR6fuRJM9GNjOzEikiep5AehXoyBaRPI9gRDezdM53ILBXRByRDh8G7Jh96L2ktwC/BdYH1gX2iIj7uljWNGAawJgxY7abOXNmjtQar6Ojg2HDhjU6jEK0cm7Q2vk5t+bVl/x23XXX+yJiUlfj8hwRfA/YBvgN8POKjuO+Ohi4NCLOkrQzcLmkrSLi1exEEXEBcAHApEmTYvLkyXUMoThtbW00S6zVauXcoLXzc27Nq6j88vQRfE3SGsDewHmSnouIaTmW/QQwPjM8Li3L+jSwV7qeP0oaCowG6lnZmJlZD/KeBro5sAtJO/6fcs5zL7CZpE0kDQEOAmZVTPM4sDuApC2BocAzOZdvZmZ10OsRgaS7SPoFLgFmpmXbRsT9Pc0XEa9IOhq4CRgEXBwRD0k6GZgdEbOArwAXSvoSScfx1Oit08LMzOoq16MqSTbSB5Ps1Ssd3q23GSPiBiquOYiIGZn3D/P6aalmZtYAefoIJpcQh5mZNUiem86NkfRTSTemwxMlfbr40MzMrAx5OosvJWnnH5sO/w34YkHxmJlZyfJUBKMj4irgVUg6gYFVhUZlZmalyVMRLJc0iqSDGEk7AUsLjcrMzEqT56yhL5Oc//82SXcCG5A8pMbMzFpAnrOG7pe0C7AFyamj8yJiZeGRmZlZKfI+s/gdwESSK3+3lURE/Ky4sMzMrCx5riw+EZhMUhHcQHLPoT8ArgjMzFpAns7iA0nuB/TPiDgceDfJPYfMzKwF5KkIXkxvC/2KpBEkdwYd38s8ZmbWJPL0EcyWtB5wIXAfyUNq/lhkUGZmVp48Zw19Ln17vqTfACMi4sFiwzIzs7LkPWsIgIh4FEDS10meUfCriLi+gLjMzKwkec4a+gfpVcWdRcAYklNKlxQTlpmZlSXPEUHlw44F/DoiHi8gHjMzK1mePoLFlWWSfGWxmVmLyNM09KZa5jMzs+aQZ4N+H0kfgTJlfq6wmVmLyNM0tEkZgZiZWWPkaRr6SFflEXFN/cMxM7Oy5WkauhJ4mKSJqLN5KABXBGZmLSBPRbAV8G1gGPDNiJhXbEhmZlamPH0E84ApkrYDvi/pSeCkiHii8OjMzKxwefoI/ovXzxJaAOwC/B+wToFxmZnltnLlStrb2xk5ciRz585tdDiFyZPf0KFDGTduHIMHD8693Fx3H+1l2Mysodrb2xk+fDijRo1ixIgRjQ6nMMuWLWP48OHdjo8IFi9eTHt7O5tskv+EzzxNQ5flXpqZWQO89NJLTJgwgY6OjkaH0lCSGDVqFM8880xV89V607mIiE1zzLsXcA4wCLgoIr5TMf5sYNd0cB1gw4hYL1/oZmavk9T7RANALZ9D3pvOCbiN1zfaeYIZBJwHfABoB+6VNCsiHu6cJiK+lJn+P4Ft8i7fzMzqo9dHVUbE4ohYBLySvl/c1Y3ourADMD8iFkTEy8BMYP8epj8Y+O9cUZuZWd1Uc9O5QZLWJ72oLCKe7WXWjYCFmeF2YMdu1vFWYBOSo46uxk8DpgGMGTOGtra23sLuFzo6Opom1mq1cm7Q2vm1Ym4jR45k2bJlrFq1imXLljU6nMLkze+ll16q6juu9qZz96dlAfTaR1CFg4CrI2JVVyMj4gLgAoBJkybF5MmT67jq4rS1tdEssVarlXOD1s6vFXObO3cuw4cP7/WsGoBJp9zMoo6XVysfPWwIs0/4QJ/iuPXWW7n00ku5/PLLq5rvqKOO4rDDDuN973tfj8vLkx8kp5Bus03+lvYibzr3BDA+MzwuLevKQcDna1yPmVluXVUCPZVXY86cOVVtgDvdfffdnHfeeXVbXrV67SOQNFjSdElXp6+jJeW5UuFeYDNJm0gaQrKxn9XF8t8BrA/8sdrgzcz6k84N94oVK5g6dSrf+MY3iOj5rv1z585l8803Z9CgQb0u71vf+lavy6tFnqahHwODgR+lw4elZUf0NFNEvCLpaOAmktNHL46IhySdDMyOiM5K4SBgZhSRnZlZiR588EE23HBDPvjBD3LEEUdw6KGH9jrPjTfeyF577ZVrefvvv38hp8nmqQi2j4h3Z4ZvkzQnz8Ij4gbghoqyGRXDJ+VZlplZf7Zy5UoWLFjAwQcfzE9+8hN23nlnrrvuOq6//nqef/55Pv3pT7PnnnuuNt9NN93EJZdckmt5y5YtY8GCBZx66qksXbqUq6++ui6x99o0BKyS9LbOAUmbAl126pqZDVRz585l++23Z80113ytmeeAAw7gwgsv5Pzzz+fKK69cbZ4XXniBJUuWMHbs2FzLA9h000356U9/WtfY81QExwK3S2qTdAfJKZ5fqWsUZmYlGT1sSFXlec2ZM4f3vve9zJw5k8MPP5x//etfr4075ZRT+PznVz8f5vbbb2fXXbu+Tren5dVbnrOGbpW0GbBFWjQvIlYUFpGZWYH6eopod+bMmcMOO+zA5ptvzhlnnMGUKVO45ZZbOOGEE9h7773ZdtttV5vnxhtv5MADD8y9vGuvvbaQ2PNcUDYU+BzwfpLrB34v6fyIeKmQiMzMmtCZZ5752vt99tmHffbZh3PPPZdbbrmFpUuXMn/+fI466qg3zHPXXXdx9tln517esmXLWLx4MccffzwPPPAAp59+Oscdd1yfY8/TWfwzYBnwX+nwIcDlwMf6vHYzsxY2ffp0pk+f3u34+++/v9tx3Rk1ahTnn39+X8JaTa5HVUbExMzw7ZIe7nZqMzNrKnk6i++XtFPngKQd8cNpzMxaRp4jgu2AuyQ9ng5vDMyT9BeS5xJsXVh0ZmZWuDwVQdeXvJmZWUvIUxF0eeuHiHi8q3IzM2sueSqCR4D5pM8hSAXgJiEzsxaQpyKYFxF+hKSZWYvKc9aQ7wpqZtbC8lQE60n6SOWr8MjMzJrMrbfeymGHHVb1fEcddRR33nln3ZZXrTxNQ3cA+1aUBXBN/cMxMyvBg1fBrSfD0nYYOQ52nwFbT+nzYpv1CWV5bjp3eOFRmJmV5cGr4NfTYeWLyfDShckw9LkymDNnDlOnTmXFihUceeSRjB07llNPPbXHh8n09oSy7PJGjx7N9773vbo/nCZP05CZWeu49eTXK4FOK19Myvso+0SxPfbYg9NOO63XjXbeJ5TtsccenHjiiQ17QpmZWetY2l5deU5dPVFs7ty5nHPOOSxatIjdd9+dz372s6vN1yxPKDMzax0jx1VXnlNXTxTbcsstOf/887nqqqu67AxumieUSZrR1auuUZiZlWX3GTB47TeWDV47Ke+D7p4oNmvWLD784Q/zoQ99aLV5+ssTyvIcEUwDOoDlwGfSv8sLi8jMrEhbT4F9z4WR4wElf/c9ty4dxVtttdUbnii2cuVK9ttvP2688UauuOKK1ebpqX+gu+UVIU8fwTMR8X0ASUcBV0TEPwuJxsysDFtPqcvpolldPVGsra2Na665hhUrVnR5RNBMTyiTpP8ARgArgGslHRcRbX1eu5lZC5s8eTKTJ0/udnwzPaHsM8DXgFUkj6f8F3AJ0FbXSMzMrCHyXFB2L3BgRfH+xYRjZmZl67UikPTlrso7+w3MzKy55Tlr6JvAJ4DhFa9eSdpL0jxJ8yV9vZtppkh6WNJDkn6RN3Azs6wI3ygZavsc8vQRvA04DtgdODkibsmzYEmDgPOADwDtwL2SZkXEw5lpNkuX/b6IeE7ShtUmYGY2dOhQFi9ezJAhQxodSkNFBIsXL2bo0KFVzZenj+BZ4FhJY4ETJR0DfDPtO+jJDsD8iFgAIGkmSd/Cw5lpPgOcFxHPpet6uqrozcyAcePG0d7ezpIlS6reCDaTl156qdf8hg4dyrhx1V0lnaeP4Ne8/nAaARsDdwOr3yrvjTYCFmaG24EdK6bZPF3HnenyToqI33QRwzSSC9sYM2YMbW1tvYXdL3R0dDRNrNVq5dygtfNr9dyGDRvW6DAKkze/xx57rKrl5mkaOrP3SWq2JrAZMBkYB/xO0rsiYkl2ooi4ALgAYNKkSdHTebn9SVtbW4/nEDezVs4NWjs/59a8isovT9PQHTUu+wlgfGZ4XFqW1Q7cExErgX9I+htJxdBbs5OZmdVJnpvOvSLp+cxrmaTncyz7XmAzSZtIGgIcBMyqmOY6kqMBJI0maSpaUE0CZmbWN3mahv4SEVU/Ky0iXpF0NHATSfv/xRHxkKSTgdkRMSsdt6ekh0muXD42IhZXuy4zM6tdnoqg5pNzI+IG4IaKshmZ9wF8OX2ZmVkD5KkI1pG0DckZQ6+JiOrvlmRmZv1OnorgKaDydhIB7Fb/cMzMrGx5zhrq+vE5ZmbWEvKcNTRS0vclzU5fZ0kaWUZwZmZWvDw3nbsYWAZMSV/PkzyPwMzMWkCum85FxEczw9+S9OeC4jEzs5LlOSJ4UdL7OwckvQ94sbiQzMysTHmOCD4LXJb2Cwh4FvhkoVGZmVlp8pw19Gfg3ZJGpMN5bi9hZmZNIs9ZQ6MknUvysPrbJZ0jaVThkZmZWSny9BHMBJ4BPkryEPtngCuLDMrMzMqTp4/gLRHx7czwKZI+XlRAZmZWrjxHBL+VdJCkNdLXFJK7hpqZWQvIUxF8BvgFsCJ9zQSOrOK5BGZm1o/lOWtoeBmBmJlZY+Q5a+h/ywjEzMwaI0/T0NjCozAzs4bJc9bQppIqnzVMROxXQDxmZlayPBXBM8BZRQdiZmaNkaci6IiIOwqPxMzMGiJPH8HphUdhZmYNk+eI4DpJ04F/T4fvAM6PiJXFhWVmZmXJUxH8CBic/gU4DPgxcERRQZmZWXnyVATbR8S7M8O3SZpTVEBmZlauPH0EqyS9rXNA0qbAquJCMjOzMuU5IjiW5DkEC0ieUPZW4PBCozIzs9LkudfQrZI2A7ZIi+ZFxIpiwzIzs7LkudfQ5sAMkgfW7wmcn5b1StJekuZJmi/p612MnyrpGUl/Tl/ugDYzK1mePoLL0783A4OAOcAlvc0kaRBwHrA3MBE4WNLELia9MiLek74uyhe2mZnVS56KYM2IOB54NSLOiIgfAMNyzLcDMD8iFkTEyyTPMdi/9lDNzKwIeTqLX03/fj5TphzzbQQszAy3Azt2Md1HJf078DfgSxGxsHICSdOAaQBjxoyhra0tx+obr6Ojo2lirVYr5watnZ9za16F5RcRPb6Ad1UMjwAOyDHfgcBFmeHDgB9WTDMKWCt9fyRwW2/L3W677aJZ3H777Y0OoTCtnFtEa+fn3JpXX/IDZkc329Vem4Yi4i8Vw89HxHU56pgngPGZ4XFpWXZZi+P1M5AuArbLsVwzM6ujXpuGJC0DonMw/RsRMaKXWe8FNpO0CUkFcBBwSMWy3xIRT6WD+wFz8wZuZmb1kaez+Bzgr8DBETE8ffVWCRARrwBHAzeRbOCvioiHJJ0sqfOhNtMlPZTesmI6MLWmLMzMrGZ5Lig7QdIGwDclfRmYERF35ll4RNwA3FBRNiPz/jjguOpCNjOzesrTNLRt+vZSYBPgR5IWRsQ+RQZmZmblyHP6aOVjKp8F1i0gFjMza4A8TUO7lhGImZk1Rp6moRldlUfEyfUPx8zMypbnrKHl6eszmffLiwzKzMzKk6dp6CwASYd2vjczs9aR54igU/Q+iZmZNZs8fQS/JqkENpU0q7M8Ivbrfi4zM2sWeU4fPTP962YhM7MWlKeP4I7ssKT3AwcDd3Q9h5mZNZM8RwRI2obkhnEfA/4B/KrIoMzMrDzdVgTpc4kPTl+LgCsB+QIzM7PW0tMRwSPA74F9ImI+gKQvlRKVmZmVpqfTRz8CPAXcLulCSbuT7xGVZmbWRLqtCCLiuog4CHgHcDvwRWBDST+WtGdJ8ZmZWcHyPKpyeUT8IiL2JXnc5APA1wqPzMzMSlHNlcVExHMRcUFE7F5UQGZmVq6qKgIzM2s9rgjMzPq7B6+Cs7dil7YD4OytkuE6ynVBmZmZNciDV8Gvp8PKF5PTNpcuTIYBtp5Sl1X4iMDMrD+79WRY+eIby1a+mJTXiSsCM7P+bGl7deU1cEVgZtafjRxXXXkNXBGYNZu045CT1iuk47ApFdyZ2lC7z4DBa7+xbPDaSXmduLPYrJlkOg6B6joOH7wKbj2ZXZa2wwPjkg1JnTobG6qEztSG6szh1pOJpe1oZP2/Ox8RFKGV906ssWrtOOzcWC5diIjXN5at8NssoTO14baeAl/6K3dMvg6+9Ne6V3CFVgSS9pI0T9J8SV/vYbqPSgpJk4qMp2q1HIK38j8cuJJrtFo7DvuysSyzKaqWdfWlM9XNbECBFYGkQcB5wN7AROBgSRO7mG448AXgnqJi6esGnWo26K38D9fqlRz0/4qu1o7DWjeWtf4f1KLWddX6mZSZWz9X5BHBDsD8iFgQES8DM4H9u5ju28AZwEuFRFH2Br0R/3DVbtQbUcnVquzKscyKrpbcau04rHVjWeaOTa3rqvUzaZadthIU2Vm8EbAwM9wO7JidQNK2wPiIuF7Ssd0tSNI0YBrAmDFjaGtryx3ETn/8BkO7+LJfuv4b3P3sht3Ot8vS9i4fvhBL27mjh/XvtNZohq54ZrXyl9Yazd09zVdjnBv+6w62mHceg15dkRQsXciq645m3ty5PD1ml7quq9bPpDPOTRdczlorFrFirdEs2PSwbuPLzlNtbn1ZX1++g/Jy25AN3/7Z1df37IbQw3ew4diPscWyzPqAVWusxbyxH+PpHuar9TuvJb/af1/Zz+QZVqy1Qa7PpMzc6qWjo6Oq7V9eioi6LxRA0oHAXhFxRDp8GLBjRBydDq8B3AZMjYhHJbUBx0TE7J6WO2nSpJg9u8dJ3uik9YCuchSctKT7+c7eKt1jrjByfNJZ053Kszog2TvZ99yeO3jKjLNZPpNm+A7Kzq0v0rOGqjrzpNY4a5mvDp9JW1sbkydPzjVtqbnVSVX5VZB0X0R02Q9bZNPQE8D4zPC4tKzTcGAroE3So8BOwKy6dxjXekhc6+Hm1lOSjcDI8QRKfhy9bRT6EmctTVFlfyZlN7PVur5aPpeyc+uLWs48qfU7ryW/Es6Xr8v6GvHdFazIiuBeYDNJm0gaAhwEzOocGRFLI2J0REyIiAnA3cB+vR0RVK0OG3Sq2aB3zlvWP1wtG6+yK7la/3HK7hit5XMpO7ey1fp/UEt+ffmfq0WZufVzhfURRMQrko4GbgIGARdHxEOSTgZmR8SsnpdQJ5mLMVjannxZeS/G2HpKeRek1Brn7jO6bproaeNVh8/kjmoOUUeO6+ZQOscRSLW59WV9tVy4U3ZujVDL/0Gt+ZX5P1fr+vry3aXNc1X/3xWs0CuLI+IG4IaKsi4/rYiYXFggZf+4alVLnLVu1Mv8TPqyUYByKsfsOqup6MrOrVm0cn615taXq8IL5ltMtIL+XtGVfVRW5kaoWY44G6GV86slt576k1wR2IDQDIf8zbAua179uJPZ9xoyMytDP+5kdkVgZlaGsk+PrYIrAjOzMpR9emwV3EdgZlaWftqf5CMCM7MBzhWBmdkA54rAzGyAc0VgZjbAuSIwMxvgCnseQVEkPQM81ug4choNLGp0EAVp5dygtfNzbs2rL/m9NSI26GpE01UEzUTS7O4eBNHsWjk3aO38nFvzKio/Nw2ZmQ1wrgjMzAY4VwTFuqDRARSolXOD1s7PuTWvQvJzH4GZ2QDnIwIzswHOFYGZ2QDniqBGkvaSNE/SfElf72L8WpKuTMffI2lCWv4BSfdJ+kv6d7fSg+9Frbllxm8sqUPSMaUFnVNfcpO0taQ/Snoo/f6Glhp8Dn34XQ6WdFma11xJx5UefC9y5Pbvku6X9IqkAyvGfVLS/6WvT5YXdT615ibpPZnf5IOSPl5TABHhV5UvYBDwd2BTYAgwB5hYMc3ngPPT9wcBV6bvtwHGpu+3Ap5odD71yi0z/mrgl8Axjc6njt/bmsCDwLvT4VHAoEbnVMf8DgFmpu/XAR4FJjQ6pypzmwBsDfwMODBT/iZgQfp3/fT9+o3OqU65bQ5slr4fCzwFrFdtDD4iqM0OwPyIWBARLwMzgf0rptkfuCx9fzWwuyRFxAMR8WRa/hCwtqS1Sok6n5pzA5B0APAPktz6m77ktifwYETMAYiIxRGxqqS48+pLfgGsK2lNYG3gZeD5csLOpdfcIuLRiHgQeLVi3g8CN0fEsxHxHHAzsFcZQedUc24R8beI+L/0/ZPA00CXVw/3xBVBbTYCFmaG29OyLqeJiFeApSR7kVkfBe6PiBUFxVmLmnOTNAz4GvCtEuKsRV++t82BkHRTeoj+1RLirVZf8rsaWE6yR/k4cGZEPFt0wFXIk1sR85ahLvFJ2oHkiOLv1c7rJ5Q1iKR3AmeQ7Gm2ipOAsyOiIz1AaCVrAu8HtgdeAG6VdF9E3NrYsOpmB2AVSfPC+sDvJd0SEQsaG5blIektwOXAJyOi8oioVz4iqM0TwPjM8Li0rMtp0sPtkcDidHgccC3wiYiouvYuWF9y2xH4rqRHgS8C35B0dMHxVqMvubUDv4uIRRHxAnADsG3hEVenL/kdAvwmIlZGxNPAnUB/umdPntyKmLcMfYpP0gjgeuD4iLi7lgBcEdTmXmAzSZtIGkLS6TarYppZQOfZCQcCt0VESFqP5Ev7ekTcWVbAVag5t4j4t4iYEBETgB8Ap0XED0uKO4+acwNuAt4laZ10A7oL8HBJcefVl/weB3YDkLQusBPwSClR55Mnt+7cBOwpaX1J65Mchd9UUJy1qDm3dPprgZ9FxNU1R9DoHvNmfQEfAv5G0h53fFp2MrBf+n4oyZkz84E/AZum5SeQtMX+OfPasNH51CO3imWcRD87a6ivuQGHknSC/xX4bqNzqfPvclha/hBJBXdso3OpIbftSY7clpMc5TyUmfdTac7zgcMbnUu9ckt/kysrtifvqXb9vsWEmdkA56YhM7MBzhWBmdkA54rAzGyAc0VgZjbAuSIwMxvgXBGYVUHSo5JG93Uas/7EFYGZ2QDnisCsG5KuS58Z8ZCkaRXjJkh6RNIV6f37r5a0TmaS/0xvTvcXSe9I59khvXf8A5LukrRFqQmZdcMVgVn3PhUR25Hcc2e6pMq7x24B/CgitiS5ZfPnMuMWRcS2wI+Bzgf0PAL8W0RsA8wATis0erOcXBGYdW+6pDnA3SQ3BdusYvzCeP1+UT8nuTtpp2vSv/eRPFQEkhu8/VLSX4GzgXcWEbRZtVwRmHVB0mRgD2DniHg38ADJfXqyKu/Pkh3ufMbEKl6/3fu3gdsjYitg3y6WZ9YQrgjMujYSeC4iXkjb+HfqYpqNJe2cvj8E+EOOZXbeXnhqXaI0qwNXBGZd+w2wpqS5wHdImocqzQM+n06zPkl/QE++C5wu6QH8UCjrR3z3UbMaSJoA/G/azGPW1HxEYGY2wPmIwMxsgPMRgZnZAOeKwMxsgHNFYGY2wLkiMDMb4FwRmJkNcP8fJ3VDeFB818MAAAAASUVORK5CYII=", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "fig, ax = plt.subplots(nrows = 1, ncols = 1)\n", + "ax.plot(a_f[0], k_f2, \"s\", label = \"$k_2$ / $k_1$\")\n", + "ax.plot(a_f[0], k_f3, \"o\", label = \"$k_3$ / $k_1$\")\n", + "plt.legend()\n", + "plt.grid(True)\n", + "plt.xlabel(\"alpha\")\n", + "plt.ylabel(\"Анизотропия проницаемости\")\n", + "plt.title(\"Кубическая гранецентрированная\")\n", + "\n", + "if savefig:\n", + " fig.savefig(\"anisotropy-facecentered.png\")" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "interpreter": { + "hash": "19649669bd52b0be75e091dcf60d2128e4a347083ff474cfec5ff9275df3ceed" + }, + "kernelspec": { + "display_name": "Python 3.9.9 64-bit ('anisotropy': conda)", + "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.9.10" + }, + "orig_nbformat": 4 + }, + "nbformat": 4, + "nbformat_minor": 2 +} diff --git a/data/anisotropy.db b/data/anisotropy.db new file mode 100644 index 0000000000000000000000000000000000000000..f1e0a110d7c0568c080d5131e519876dc202a190 GIT binary patch literal 188416 zcmeFa2S5``+xI=`gchoZ$|isypddwx3M2?vKoO86cBCqbh=2`SHum0o!`|DmcU0`X z_uiFb@AaLXP2y(Ra87)l?|r}fdBT04th?7gzn#5ycdnV)*(wdBGZ1yPI&9U`;)}&% z;W4;s`ZIaJg+%(^I^B$t)3!I#sSR<^8 z($aIY5Dt9^G|ZH!M$Ka>-z1E1zix16DX~prlO+i)QzQxU=FKIE$<10sCAX2p$F`{> z(e0sv4%!GttD>}`Ua%}Fv8_`yn_~NBW%Me_>5-R`(>)7Xnze|9mi&UO!opq!S&6*?-PNbeD z_?5=7$*~F1v8_ngX6c_s;_Sx6HjjlDL3C8B=%|=j?-U;2sg^UV54v>F`7Z2|mY)UB z9kG&C$k)GU^c3XPLh6XiMf*)RZ*?!Oz{%N}H7pigFQkP0)%{5y1g)6 zGX#Wo%@CWIubHgetnOKPMbJl@s}x=_S(%*)r~TfLotB}$p?8nmUXV%E6lE0^q;<*4 zEa;J^w>`fgr?9*JmVO0k-E%Uj+w{++?&!SrPA1w6zeRO2!`W=bL1(iUI-4Xa3*Ii! ztt5b}e>U|L(`{N(7L%->%||N1v)Q@azcm$a zW89u4aB_8J;XTng`u)Y<|6=Q(tKSCVs?xG6oCaEbor!9u7#K!~dY5{#biCiyvNzE; z=Zj=ObXQa=DiNtfo#2Y;&*TA<2TUF?dBEfWlLt&5FnPe_0h0$z9x!>p~^3h$JACum=SB1^fF41p8xL275pQ z|4_gB!GZOH12GQW9@ro#v_U|aUr2CpNJt>YX0Qh}2oCiN2nzBK3dUG^_Ii?lfCj;V z1V1RKe*Mre9y}UL(N5x+YKy*!K8jw8o{JueZi}vp&Wlcp4#QCaWYjxo(~8LhCJ&fA zVDf;;111lcJYe#G$pa=2m^@(efXM?U5B%TrfcKC@xHGpQ9`I8#gqWL;sWaebb#*ZO zbXMEJPba(w{B*#Z!%sWh3w~PR9QbLegdY~ZZ1~9%eSn|9k6aM{H2s-8VDf;;111lc zJYe#G$pa=2m^@(efXM?U512e)^1%ND5BTyqEH;nDC7%CtMbC&s@`BI*KZstDEB_OY z(R2(Z512e)@_@+$CJ&fAVDf;;111lcJYe#G$pa=2{BQJt%a9Ow*pBMP@Y4ow1wTdI z$mjp>33qsi(1Q4<>CfZ=lLt&5FnPe_0h0$z9x!>p-c%k=yoy2SKn@_@+$CJ&fAVDf;; z111lcJYe#G$pa=2m^@(e!2cW%c%tY3kHI72_y63vHGpWPDB1d(^%CnARv)bvST(bJ zYdOa<#^P^_=@wG+r{zC2g%ac&{E z2ImkbkK@EHW#_Q%Si4x6EGw{;#b@zhj=`y5{(q-s1e5mER*b9R5_Ys~BebS9?9tCnII@Sqc18?qBU2kW+6PCR4%kfK`Ew4hy zRzcVx@v=>6H;i|}hNCrt6-g81j&!Ue>4YWG{ck54kKK{f=iW)hdbtA~>wvKSE@=l& zUp3yh2iSuJ{z1>>_H?X0iQOptHf*c$r9Jw(W8FEgfr%us&<+ELi8K)`rAxzvVQ|!FcTTPsO!s9C{!Z(Xk?g z^_|U6%jynj*FSBQ;dZ#Yt`ia}G@-;Kw5wNBXiu_kWPH z=~yQ7>?UW?u^fbLu-xh5$(=G9hqhV~8~!6Gs=>?i z3bOxS9SCrN5B~yx0N;amn%9RX;U41_(8dIsTE4NIZ5eG*ZZXXw()@|}MDzN>`@%88 zdSueIAfErq?))QH>L{r*Hpp{f)12&M%wn@=#rh*)w61ybc{(2VvbeHADbx`Thpzt z&BEozm-fmMj<0KeKe-Pb>qC}y(h!hy()iNeDq8V*blF$AHy!IumKG07-Py!=?7ic@ zR%2$2lY7yzUSw$(6fDd-X1o&~;k%DT_4z5UMaR}cSWnjJ_yrw|FYU8w>Vl{KbLF0N ztS7>Hg*W`P|1aaQfB!w9eFyhGat}JzgT%^@yeq10yc6EYmfZLjSyNtIiXw#tCLS<`;t=_E&aG#wj_u=S(6HM^K#d};4Kco=*&yt=#*9ovX3?acV(GycZ= z_91J8qUqB|@+dks3Sk>^cHh|@Y&`Z!JJ-MDjb6%SbgT?v!&9RR7Y;BU`vU8>W|8Yz zxs;BTB5XvrWA~TkmeS4gC;h?SekBk9;kva}yQ)LYU}M&r;{8L)3wr5#YT z70Dy$*a(D;?B09ms6)n=_PyA4-t8fEq7 z8%km)jF>q5k%U$}+G@qpMzFWcX?Cb-wmgK64MEtLZ)Ne_RVo^XwptMz-kM+3_p!K+ zJeZCRM%adefzu^#N@*P0Dg(A4W^ihFWIcI3I<_9d2Csj9|M5EGOFQ%(f6v;3_VOS) zHV9$;dW-}^P8pA#`gfI9FMfTK2hy>D2wQg$R;9?^c(&Wx0MIbu>lC{vw!_; zzcfH|CT*3$w`+%qgWrtYBKN0b{Snry_ohSZCKzAZt?~Bv)|#iv>(a4x5w=$CX`a8@ z7?0i8aew2v_3Fv}=vY65^^89;+~W`0V29?<-eeh~BKrTnY)>HCEb1bvVSUzmsCAfC zxmAf(ij|q=Zp(bjnif|r##+eC-)-k(lHr1@D;JaX@Kq27pH}E^~ z&3W5+nLJzWUT!yTRX7izfFov~WcOoxvM#VxEMHIrj==wse_sKSeDJ^bPeA9i&IhBk~ttF_-zH zBC%w=+?0EF_Tth*XsJqq{k4in8EBamIrgmqu-+IxLCs6z_s`Ks+yS$brMs(`Z?i1p8(ck?(l9U2+ES4N6XA4IFEt} zSwEmfw!{E(Gwr&xU-s}#os-R_7tk_33C^LYh?M@?l34iWy8+}z+C}6r!sgNR?bW3h z(Xv{SD-%)%LR;SrYi|I#j&{Z!)BFOG?@gCpLd$7MaP9;{DSs}Sk^Fb20pv>B>8a9? z<91QnNO~D9qmtmf2?pe^{HGtvwi`e$q50CuHjh2T#Z`I*EhR~C&IAMU2Z;H`KWYFu zpLXC*51W>b{Mb`^1ugTF;Cu-N2%**x`|xI!D$dotc$A7t7aQuohDDNjr8N<{G*L-U65`imvM0J z0~4$Ea`D=`23V(v7#PvhM*8KbOS_n7`=vX`WgMLFP>}`e*KZM5xCU4!$;!}1`a^JS z-0(-^rF+O_9GvCA#F~8Vj^Aektm8!3E52tV0|Dud>pfJ`QgRswr#3LKe(jy6EYyV( z+BI>Eyk%`2PYdavEpMpHJWRHwL^Q|+rp;9 zg%2N569O*d;0yx>){kBL2G2C)drcFdKt+~})n~@e>UjoOnqY1tgHilI{iC;uC{8Tn-~<5%7U$NAaPC(FEKMWJ zMh06~kMj@IAX-<%G7ioPU}8nD8_F?M7i!uquw1rzobO8#i3SU?jAJpeq+g4U*Q?dZ zK)Ezc0UH@?1LU=IDks_i#L~*#2l#+|AiN?RDhz&@0_*49>AtlC>~ELT}3 zS$x*CCpdy7{<`PtN%4AR*YqT1NTbjzJGyw_M+d1UR>{IpjAzTt;_MzZiEIEvR@%$D zgU%+!F(TRTo89x$4-LPOnvkVKn%{_yhzLTgLuP<{H_DXlLN!R@8v=JPdmpa z!6!~k$aW8XZ(lWl)OMk@rL5soZK&I~w?r6h?UY=S1fMoBA?L{BhqN|;)PzA@DP4+{ z`5X6#5Yh0%iOX(D@JSOBQdoWKRQ(={_ImFi)pnt^rNoZU$_U=ZA}aav)BS!W!KX}2 z$Th;Qrm^rA|+<`m$wXwC*o|qoO#=m;1ebWBm9L%33?Ma)0$nNHMhO=+7%e2)2I|k- z&6D60B_?FA+%+MO4InjfQ-^eSoqJ=XZ5twd7Mt%so&=vJF(J9q%ky>3f9)$Wnp$1c zh1L#05j|<%(x??gPZ}KtadLVJ2>rPOqz()`e=X5~fy*^X@M#efQraZezpa5%Y8xhU74Qm`&;X~1@;CvC0z^@;RoUQ5x| z&524mGdT7_hbbYIf_1@519r$6Q*U*ZQ0eiZe)T*<`zliy!J?Pj3Z|A<3f6@%4Y-O~ zYU1p_)Y6j+tbeh~9x7Wgf~{XSztT3PQm`(NX~0f%!@3<=4Ij>|UgO_nfSEFx865F` zhwsEn!Mbp!0XqkQSK=02>Di8!dpmp@os@*Y+R+AW^PhX$tte@)w%X7(hiKtU6UU_= zf^yhxucW|mHkRXl>|jzNGqU8@WOMz&BQ3XsR2$B;$S?j!Rk!%7iEt)=y1sf+OJ-!f z!l2>r4V3b%CY))JKYDfj>*TS|L^$gc6p)aVz=-6`jY^r_)BsW&&a@2}@O5hZUY|Nd zILrJoZhTS;MkM?5F0fQ$0I3aUS|mr|SXwrz2@%d%d?!Itb4DcR;P>R2Ck-IA;Y^F< z9lHPGnS~b-&KiuhzHMbK1a&|!mc?s}8E{trBlEH5LBi`oT30 z!SBQu@YeC<-0$4Q+&In~&P=K;HW4VE74-qEL|`wwLwJWHp9+ z?)TNK?KCv%w>)-YkVNK`8MEj>-<#61sGow(zGE7-t`w}TiM3$o^FRD; zgZyDl-1Ka3=l9B9j9~H7Bg@`NDg|qutOeHyzV__ntuC-8P8jr(pQbEg1lOp6ogC1& zQn1#?S}-<0F}iC0Fjy1&?)*MLp)6zuFLi&`_;#gWZH1`?yKeER?&w_`R+y_V_~sWV z3mC!GH>q=LS6*Umg{cL*J(kUQ($j~iFqh5#vPIdG5$qm2J1_f_eu=daR8wu1X`ao6 zp_%gHOqZ&Hh|ziFKnX8t~6mrH!K-jFE175t4f6(*R{I zBiL4b>E6&~m4bC+DkNCO&P$e`3BN1d_NK|?mEpaV-59}k!i{(C4X6~X8(GnSc{3BE z3qBl{?zS43+h%)HWmiV9&CVg-FUD31){U=dz-C9+e_|($kd{{cd1~u1du0wI*tWX# z_T;9Of_0-T8nDIp7Q>Fez9v23+vZ7~x#yK#7{NAWJ#IhUUMW~xZE9{1naIDrS^1=2 zM77zmZI5teXJ+uOgdg`ZD+Oz-O)Z#cz_JJ05e?YRpF4F`W;26Vg|_|jxKgmT+SGz! z19r=d3FfE)Ydv+pGK(2JEwXy_ib}!SYEug)8nBY%fkXo~s9<8eGLsp+e6*q{zf!P% z1J*}IG++f5eTfF_&qf~yDKi+sqS?(3dX-fQ*4D?`5)%#BI-CHa0b6%qlOSa}Gq~+S zt36#R1#9bLEf_XnU-l2IMbyW8-gR21Ok)JQ5A&;QM+4KkHTw0jH&Gwgt+K9wsE@lG zZd*&)i5cAR%n$ZQ{Swo_`t`B56xPQXHw!lt^|AWr!QskOMzC8>R!IuCQm}q~>@9`$ zv6aO8B2gbVe12kwvLhqd?nvWN?}t_j)~}DfrLaD(v3|vMqCQ@jdaakT12Z_i%a}Gh zD+TM<$KFy{A155XpH9@rS)DhYQ?_RW+cx1Z|4FNlX+GAkkG+ZdxNXnmJ4Agva~aA`7C8?W^nI}15LszEwO%m z>`m0i!V}z=s6KxFqMelz2aIUvS_|iV^{lk%t*tFIO>Y^P=;*6j#g`t*O?wfYkfI#K z2=-QNIgUj#lo@#hesqv(Z)#2JP5P`EP}lu_1K!j}_bz?1KslHZT=jHYUujXLCDv9$ z+7h!luQyMwC6pd3j@hw-ldl}e2=?CYH>cdQ(h_TH3vG$tW^Vp<@W&5WTWpW*^>&*Q zXEBni-P)wAl~`$!wROF=fAcf(>Z|y)VCkVA`+D|U76E%{4J!83{>>KU0LCI=SJh76 zk}Ew}ZAGJPT(fi21zp{8U`4|@HP(HlvOhDpb<#x@(@@gy^RBIEv|s@cIVeskU_}$0 z?b$-4EM^A3ci49dR9a$fMWY2wn%T3Schj}h#$P_n+^%}T-Aibe~z z32eB@+_MN)Gy^8&traT!GJ-KT$>h0ARpW10Xe$~m*lzle0!5A8L`Bp3_AK=NzYnVv z{MX<97vPll|MS>RtWtI;C}VwQrL)pCml10rQ2v`M%18~w{E89 z{_cvI>5~dt!agF)`FBP+0P?caa&ogWW%29AR~Z{_G*}M%hWh2~j>?(r*Lgd~xzy-# zZk{Zibwd6_4ts|B*j$87l=)X*wxq;pZ0Yu0OAfoPmcxFbKDHaeHf!?B`{GQM(b!|_ z-jrG1P|9JiP#@bBVPog^cq=<&Jofyyo{x_n|0sujLVauw!p4N1zdfXh%;?fy`!#3h zom1Q7aP(Fm+Y(`;`(|H#akbQFtZJcyG&$0i_bR6nrt+zR8d6Q}L9Ew5HX z4o7kIu`LibGOo_t>F&lmVV>`)uInx8%Hc?^KDIf+HoOwpHOCAX?S$o#@~BBxJLGV7 zl|D8eVZ$`_|k5F+3S~ki9`-(S?Ob& zA#A~mMc81+4vQkc8;?DFyZrT0)_FOcZKaQGLYCHE z#m$*gYP4@pO-hg#3`>*40TO*|7le(uyqcBN%lOjXb28t{zj9R$he!0Woe{S22+5rD z(~QSHnmm2gkJpFgaBxH)n~kus*6t@33^Cpb&(D7e6nyL_heIR!*eryN+B4@tO_}l7 z*JJRhk6I+l;lPMKHWOhZo5c5_$aAfk`$gs|B*!9m-{lo*Z8Shz>vkh@Y2hePzSsR-Nl{?w5yZmH4OzJJ^A81A1S zhl3&d*p3L>+6OrH$&nb1o!GzKyBo{5%HdFmKDGnGCeDBR{?{DiOFRFD^p|`?v>XnE z=wsU>Y&@U+`DKLh*p*c_ADuOBupAD9=wsU;S@3W3)9iFmR4rfv5W0Mg!bmW$Ku_shUUkQ)e zm3%5Wby*JQQ0QZm5Vpnr4_=v1jdwz_7uUaNR+=2nqtM4DB5c!(%D~=ykb-v95$-U~$cLTRU=>%06R-KZ+D&ccD%P+>PU=!|Dpf!>)?(&m+x1>s%JQ;i*IO#?NEMLH zuBd9Y{+ifQ#^PRP+wK-aPyIfLWCN<nZWKCJPN#qFsAvMCh$Vo^WKzg6I` z%dC#w#O_-}Wl4@=Z8ZXP8IaABfM&A-m#+i;(pfAv)`n^)+Zlch z547w024v3$558ipsdlo#;Jc*zyuKqqr<&mIN=!+$laaqp$ijwaGC_wb;MRUjLA8@1 z+(+&1SYs<_Wd_bo!Q@ms8H>GVIDbAf8YEl>rTefHs+|m~o)?C&^zRa$4VpXzD{5lNR67~w+dw5Xg41mu)yM}c5>Lm zW7gVn!+k*FP_W7#i=*1f@r#;GcG#?K4q6-rO9x@isCIH7!mZsfkik_oEdrB9V(@-K z2Z$~qazw!`w!BaObE=p?FlGqWglZ>;1|%1<7ewq+MScb&yI_r}b~2Nu+cEFfM|y$q znqYW27E860$@*Ne*{-)*sT!OG!>3^}R6ChV9``NnLbaFv-eKF z=BluvV8kY@5!FuS!ISuvb-9&^gU~QAArXt9+R3DE zt{1n>dYTMEj)RF_SUA;AW^Z$ITx@ZQr3%djQ1J9ZpFf=b~3}6E30kzvItd3I+)NM z3#Ho0RAu5HwOYJ#1|gncaxZa+#*S_=GUr&p)9Z^=-&J7DWmYW_7EA?@*~5HtYg<}6 zf!veeygyctYA4fxdA#fTVZkj`o6TU)94v@xC-Zr^2A!z?b&;ysdoZ*q7D%;|iMb>{ z3JzKgRK^bT zZ>rgwG~4ewIp7QMOpqtKcu~z>$n4OhN0o;zz{5@6&!raCT#GdCWJeaR10FWYJQq)@ znao4#*lG96s+)kv*tVBkJg8<5^4kT+irTdVHRrb3=u(qvu8GW!=|7zJT7#O6Qd_yW zQ_b$A`NzV<)_K5vex|dF8`bQF%(hLJ*SjtT?p0C;xVTcyuE^}P^xpIIHNZV3Z-$G6 zYL*~#Rngt)>`uV_cGprDG1V+aX8V%bH;$$P_lt#U7mRAg$Zt0*>65q{!h4T!aiN-B zklF5*(C$uW;GSMM46C6v6Roxy$Xwmk{8|tO+}ndO9mLf&c62XQNA`d#2`fCSmx8h8 ztXi&EH7cMQ68QWOcyE)dx;F>s)tEEY?u_iU77K>h$w8Z;VD~<(D%D;U+1*`i3Xl1L zxUOJWJIsk{cS3fz8mbKzPO8|xU}&yO6|J4Pm#QE;ws2Is?=}z)iqE(>Qq7LYTXc{lU2LglTV!_F>N{~}P2g3mY~W%;HQSKp^XIK=W^ho zYVGS{K{Z<-vtvisH@#TEV@bPU7jvrFoc#9L*7@Ief|_UB9CQ&<%|c|Z5^=nX4I9){ zcY5YxMm3vhzJ1?1c@S`4mKp0Jpqd59Y|G8-o!tw#tMlf#@Tq1#GCTT(-Ket^xPR|< z-i1ds^N`ujRvGZZ9Jr4zI_<)xnz_ha<=du$KEr|g$Ufa%I8-wSneERG*}KRcxXU`e zc41S^Z1USTzgvAx1MV*i*SoMRi8sLjKL1#&)i=4 zR(MXhK{!F!L)cPSS12-j1>XT!XEx3(*DT)52lmjP!XEl^!7xFFAW9(Sf9GH0Z{<(n z7w{ALetb*bbKX(jQl6TZ$_wXJ<9_B|=5FRr;^uMVxn5iz=QZaB=P+j@XD(+1%njIq z6V9o|{=&Y>-o~EJ?!``F2f<$qc+EP)TE`m4%4Nl~LRhXWOVJ(C9?@J;f0$pOzNo78 zC+o}BTdb#A7g{G<2f#6jS5~L3)>@6V>Sh&ZlCn{2giAXd#%PLZD+z&dnX+-qKamMn6L&q>gq-a-*`~}vg-!*A z!w}3y=K=}-LIQ`fTc^80L0s+?42Ln8BG^b0TAki=82(Blyc94TMq!Fz0}|MkykD~u z5}x)17!KnwMX(+TZ0-A=J_iMJ3(sRXjKmbdI+D<*#+qy>fES|}hq0Ix5S1@jiv(4b z)xW~&LFkgfa2Sp$h&4!2#f8i33K4m|#$z}P$P~e9Byc=6GcOM!;5CWiFeFn1tB}B9 z@eFgqQ^*l8C{qM0Ny5Z#4{kyNas&*^6u}B4ux&qG2_OL-Ck)IK!E%jY$d_hN0FM*H zVQ8iZmLWlv;1wOW!CkGs5{P>T!(n`;AeJJ9O~)PsPCx}b z8+JI1(3ERuF?sF4xJ=o&h};OHF=gXIawCkvl#M0iMi_Z18yAopVZ5bmoKJ3qQI@iC z9=Q?5R?5b?#6}4Y11dTY6i~q&bRY?j%BOiGm9o}lvtJIA;4q|8B(ssEQJpa@#@{Rj zZ=17L05*4#O*wpgIjqN0PuNACuL7p!a)ly`~t45tgEuh7|sn?{Huo z>v8}by(-3GkfkW5B86}BzLU8HAT1MYk%(~^XDNy)NKt3_NSRv_Xmb***(1SWs6`6W zZIh8gKH%fG@v(qi?=ZLCJqZrOEk!a3Nn*0MpSTuK%5qb&-D*g17;-6+iAWOK-sPp@ z3;?TZvc#?u9EM$rWCD@|Exa+jF6=!_TEw#P6XP)SA_?h}@kkP&{Jv+|H_#sg7w3y{ z7=0;cGeQ(FaT2&W09g(_#obx3!sHL7%vv%Fa}c;W01n*T>Y;d z9)qUIU{otH4#O}-F`86#9J=HBSJ3P}n6L=LVI(FMs1_WB6g9*4rw0!PG57m_kl-*D zBLV5wkw_5ZX;EIge<{ngH`~Qug2Q-Bk&HkR*^6g*{3#Wg2gjLQ_sFeC|1sJeGYF90T;XW1+h<1jWO3F+3MNK)rlf91vrAm=jJ z*-nhZ@Jvw*K?-k+cTJ_?pxtAz__Y{^5t^be8@j$9zAQB{Rk?meGx0H2%H)rH?@c^n|020JZJ^RHbTgCF| z%WnQm+@C7wj|8&WgA>3Q8MukD9_WS8vK z32b?vdl0*T-E9n)=?O}az(F*it`zzKbZLv>Qa!;wQt%Al>;+Jd9K+#oG#PQxMZ1?2 z4CMUP7w!V}KqB-6dyt??aq;K`C_sG}I3P_S>_!5+ZLJQhfCAJNgCo)u!7h?u({9LO zSO(OGfkV<1!A>Nw4Oq0p9}-X>298Nn1UpE<)wUf3P=NYsa8Q~Qpf1WbByc)&spCMn z3$lv>hovcq?MUI+$=`Q2M8IB-6>B!A$%+C7e(Ue8#T<(m^S{lfnm2^M^*2iBZ+69O zkeP?zq@YkxjlZAYg)ic5<#ph3xvN=G++@xdP6?+8`z3oCJDl}^HJb&%J}}5|B?xzZ z(6IGC!lJ|N0(Yoz z-IAT)hz>dMJ$2*fKb1T@Tqc^jQ-Rx80F!;^1}*dYggrCMB{NB{ocWk!R$@r%|SdvsNSizSLbK_mHy}f%x z>Gom;ZdL(I4xY5C+TaJv#dvT_i-qnF<}_B|f(l@A$ffy+#H)VMrAA*BT}qb>WRG~O zz+oC`TJR(XYNCeRsx?t4Ga5W7__`#%K}QA7s{kg4cN%#sMBG^_qru~{Mn76UteFDm zRsfTOLiJbPX!mrr@!+KnR>jrdx+sX8NHkZFV@bh3Kesyk2Pbd+-q$`mo?TQBDUoPk zatO-r<6t-5AG~L?bx_;d_mw;a@h&nAOpaXn4<6{B@CUz1wQpCaQVZ@B1@R_14NMMr z1sz&HWYtLHeO&ehx3oPoO%C7phI4o*O*L|y%)fU%PHejI*b~3b93Q2Qlf(DZ^|9pm zTKBFOmb^-OZSE!%Nw9phx?XOTf*11_20o9 zpm_%~E+b>sloF%eAcwDM>z9)p`%7N1HEF^85~Itxq~O#2rhen(@Fi`1EIDGB{5l>0`+e)JDHnUyr?DJofr=*2qs2>&s!jF?}pK@Y3aj@LGMG%IuLGvl#OuB*?rNf{@HxyJOdl`Xq19+B ze>l5;fTf$oQHwka2lJigoy_^dHNs@G&t~(@Vg%0x69w%Ap8Plb-TV>!WWEdUDQ^ld zlKYH1l`G{u=S<_s*e~F7{U}ydaFA6F@>ny-E)uJp<^7-7S^9raX7-Bz2Rk+G(&Pc7 zJdligp=ybeF&0fwl71aNpt>`cy-21 zKVC1#eQOcSkl^WerXOu?BLC05Pu|JBE`2q3gV(hmucd2G-84_jj#T2+n8^#H8@6mx zDOsCCNK5`$clyYD|6S6}&&%XLKDSlk&dg*7*UJ`yO3B)6Wm@u=Eu9|~Mg!^g2zIme zPiHFes*L1X(q7KJGb$x()0Jt-KRVal*jp`-?%KYIr#^Z^i90cqo27SL?NljQn?F)R zmVxCrww!A4QM#`S+xgP<=}No`BiSSW8}@2LrDV+=B-4_)8@p~Fcc{Abpt|rH=Ut8x zcVs44GrQhrN2O$K@-Ho!epq_J(8!PSb?*X|#9%O-B}&h3?9nPe*|uY)WNq#? zE!pg|Ro0o4Z>2{@4{$2sJya5-!3<=#4(G2%cdnGIxg&Jmu&8ZcMZR1uJvJ|4aq^|V zl*EG<1~Qg#%y%HI%-YbS$-S#h4@Ra!J1M=F#bOUlT8 z_h#ysnMT&$gIe+<>w5(yH>0I%kKBlT?bugIj5ITlYxr$m7*D%4Xk_g@s3reB!mDf= zyIi{Yp6mDLcN!{*p=KuXnOF8d94ald_8!!d-#B+IoVxs&bi3cIGatG>QW9g$Oyo4n z293H@O4i?=Qm+Ivt-hH1cGJzVjSNCQ5=+M%727;a`Fr`z_s=T&K$ zwfCTw{4?=d>Zb$;A_3oikMy2OV!WAwT&-reu&m3KlC}4shD_v7Y|_y63Cy3EE?W9^ zf|3|;W*|EUs;un?R7%#~gIY4ob9l0A-zG$!!yBE>XDW#iX9hBsvpS%d<_#36v=?;; zss5dEO8Ut9ZjV`Msqjv@c4gkccqK99%s>`LPjPOzO25oBvNn|-S>{vHXRZwn@xL9E zZY=xy=GNPvN@C2Jfm~xy!liBDm6G-E!Bf(|N7Op8=td@S4>rGl6_VkgGZT5o$Du)Y zD<$jSgQui#oO3?kY}%f<2i-b<;3& zEz$q?;yeWcbN)ttYu+#ThJ6e8H2)x1&5aW66AcnYSbw+PWnEzHZFSdbhE;+U&vLhA zp{0+-J&TzZQuD{=^UPy}--J7b{e%*;BW8JKj)L8S49-)|1WpM17JDeWHtRgAA4>v` z066UuzNihG+NM8~2TUILzv}@79$w*@3OT1KHn5-XrisAV0mKMy>P`RW#cq*Ae+s@( zpn0eQi0M@5^@D#0lpFG?h4uiRug!uY{o<;^|aZN?=kh#nRVOuk;x>|)R3!+shM z?tJ0z$-?H36+{P%1}0z1h;CG}1aDwGxM*B{@ukMY6h!}u1}0zSkY&~GvAm`6VAZZ4 z^XzVIP!QcK8kl?mL^gfFoV-6637^cJ&njMcT0!)#XkhZak;vY|r-k1I#$FBKvs@-^ zY(J-k0&iFWOwJaJ={jTH?h>QH5u8@4wf9zbR^SaPfXP{gjgQTXerRJnm_O2&9lHO6 z0UB7iEtVoP?MTB5>o%G~Z-4%FP1u!|Qv2nY6w`PgPgUz-tKDvE-ngS24 z04C=;ChSmmoZQKHu!Yzz@51O=3Ou9&n4D(WIDYBLC#A-Ng%75C{j$r9#DfX2UcHIl zE&O-BCF;$36#(cB$A1fGj0aHxsM_1~XVqR{1t0x|CYm;S9p||q-*_MKvI828`>?bz z9#El#8~?0?{uN5N0hKT&WTkM;2IEVpy9Mi30Ix^jSl2^ViGOfGYwncBxL<`sU-u`6 z-Wd0-09yNJpgI*mYyJ#Wy8>wSpMiWTfL8q($h!h)<)4Az8>Y#y^{uu~9-k^b&{0S(Mc%<qb_8Tg|kJwtQ(h%QD8|mBnm}So7EBbIcnH-w5Xlo0z>d zn`hQk@Q(9|vxw7z{h7U(od7fLEdl*m4Ol#I7XFX^$Lay(^lkqQ>jL-x28an^@_@+$ zMtUHbXi&mAO4?`4@aJQ>L|aJuRo-I!Nkb2=H3ROt2iIRdmxo___6z>Bd>-fflEh>j zHZB>VuJW_f4M8=JV0F;1JH4hkcP};PGWE?g+nV?@UWcMgD1l2sX{|%+j+JFsE zeJD9;<(q8I)|ZB$|9G9k#wZi?>!ZY*j}1Zp@j8VKQzqz_MRUS;8-o7hb=sT`<)Pd3 z$Jq9Ds~Ljs$zOcMJff%0r*9Yc&&l=}Cp|mv@2)4DrzGR?%+T-U9TkCwpqht3I_Nu( zJin?#yrkz?tt0)mBqZZ;%+Mb-wsvT52)grLlXVWo_kKf{o%t#rJ4Sl0s8Pb>7DJNp zW(J@kFJ^W%1l2rL(m~(X9HN?bp{MlRgC~71eE67*H#GoF9Mc-k9TasdD4rH_xH5ne@MoonW5iX76)}S z1l2sq)BaXW z$4vOyBN>k}0QH;0yJ!fidEBRiezud{9Cdhs^irB%pt8e?WL#zdI`kMT#SnCj*9Un1 zy?;Z;EEKd2NReJT`Rzi8y=O8mH2@v5zH_=EsOHyFb)e!U_YkoK)s$WTy_3780Zqe3;p#K>6;h-@i)b`+@ zHHM)7So6aHVTDHj`} zIW{o!z)AKh_DFU%JBIDX2CUnx-7w?70%jtp$FhTm3B+n$0t_N$0;@3^to3kc;=@*f z_(fnrRk2!^0Rt&wk;1KR=+R3ufKow=)w&cIND+e+uIjO~r<8)md0@&(v09hn0x23H zh2*69UYJ;_(GoDSrC6=Yae)+3NFn}^hur`GC5#rUb$K+9BAVRG^e%pr!DLML1Gms(}yNLxmd{{#~rr<-$OUhNNQ0$L?@1 zD1Eh9txJZ16b(qlrwf-Js{l*{Kem=wt;>di6!l5P9)<5}6+kJj#cI8D7_2Z-5v1Jc z1whzcFnYFFtxJc27@?#>e64OHs36i{i`6Ixq&8nBD+DQAefZ<1!*bYxF;cM_C4|%| zf{{XU;MnT6GBQ=RSdH>d>J;_Jy_}KEf@xa9I)O30#p)n!S`Ji>AfynNhSv0e3L>4h zSgp$yi4=jPqQC|{GEzS%Pm&>>6O$EDI`l<#a$`|q2XZEHL+SRVI<2JDO_{A z`qqvd?Kl8J&x5haVzn-3BvN>jip{B+kt#9;xLB=A8i^EMNP(4H%~%f^UBTF)Vzn-7 zBvRBO%i$a_wiGIeG~r^kr=CANk;3)YyB=%cUie_f7_r(zPvJo-vR#(;D+MT(xLB>r z8;SN(liUl(s_b|v2pJE?<%-oPais1Fawio5du!H%3L@ROSnZ~_7dNC3AMAW~fC>a( z2NPF`)hMN;j^Row;ucK&1jv-+5;e*%2^A>Y7)yc_{vYSvcw7R=q~sE{ShE{aB1V$B zSAImUR{@|3WV!qlt5NPr-F`46@xFiD=7EgNTrO6lM3g#(3#pjV?;*PsWT${5kzzH< zN~u%SKnl;}JI+A`%5E-JSJyk<>PX?fyUG4k6=-=JEX@87bTrraGe+9Nwp+xGKT9=y=2^^82nkwQjfHh3V!gNfnOHGLc z4oFbF`J7l7tchIfnA%=@q)1?o1kQ?M8{vH&_`Q1>rq(8=BnxSW1XaE+pZ=E$rDMm` z+O(9Uz?KwLdz<@81^h}f4`OOwO%2uk-*{gvVFxWl++zl>+(^eqp%_g?c2{?00kkby)d;d9VHT2l7zi)ebeEb zMVR*;Q|q!(B7p@dsOq*BUP!f`r&wTWT{21}Fh>HrrN=F%(dBglt`eJHb4oY+ssL?00-zI5V0^i+)Ip{I9E(IkL zpusGe{dDuC@Or3uv;9y^&DPrmn=C^?>O$C*CeqYnYFz?KL|~DE7ThpLSO!JwQ_ws8O<^Hn zvU?1N87L`&A4uR}d$ivuD99yJw80FN6v1~SaIA*U2!hYBU_N{d@2@BLh6J{~v-`Dz z>9D?bf$7nT^#otZT@<~J^oIhNG9SZX21?2To76HF&t)~qzFEdyO?!);Cx6Jl<@<@VG2r$;3E=LiRsd}1PWl{eGG>=C@F#u zNZ{D-{IT0`7jrvL!El&_k|KCd5>5uU903KlyOv^jp`PF!64=d)NrrLG{WOs!8m6J7 z5Z;o4+qKuXgJpo}`!O8mp`-}jkbMQ zP*Mc1kf3Uh=taq}42z0BVR)XN;3X2+-}oVlfCS=40C1Rxl0x{K+(pMmqkBTZ(B8NN zhj}QG0L3Oyjs*VA3wOFmU~Fp7!lsMyZdwKrX~7F5arTmSf>8tglz;^9st0+F6afVV zL2KX{ML#AW#&h%_&yd7>&c%JWH>`rdZ20VCgY2o~5UFgcPoIzjb*CG(Ss#XKEG1T?!r| z1@lT#?{kj2$chVEwMS`lg+=@j2gk3FAVR)*Z;0`GWSZD7J104JT1Pt$} zC%BCS)gG4AtD^!zDS0r%cn3YfEhKQ_m)B|s&~GGQczZp;O(dw2xvVpc^L`t$`(t=J zJ;4nma4diKaxV;U@Mjb-yse(#Iw_d_Db)cIeA4b zaFrBPyW`LTfI5Diaxh$}C%A$HcKMkH;9vytV+ zxg)rKqUE9l>krm*t)s1;TTQYGv%F(D+|n2J_4`}6nIALHH+K^56=n%7%{G~}HDd`@ za4vBMaNOC)**)1+S^HSoEGw{?GnF%yxD;Vas0L_5{_DT~X?ly`zwSd5!Q=sBJ)j^S z9>I`H$)E&|peiLaqo?3LFnS(k_{Vu*LZ7P>6~qH18u%~*hyQiN;-<=Y@aHJr!&9w~ zD~N|kH1Hu39HQ9EGk#)l_?N?=x0O<fJ3%>RMA1u+pr3fr-p489B_!393GP>onl-3#@tSU!v{<>@LmLNbRt5yy_4}^%~WUw-mwChY&Qjb99Xob z)M#+{j-)>KPAb+b@D3HgWP8i+dff}De*oTFoP0v<2v1)73ShDs7T}Qheyz;-5+90E zPnwYRK!LZb045uNb#oGTm8y&fA7Bp__y;{#;B6~_$#$Xtp8VqX663)~AGp`yx8AG3 z+f)FP?a8`tQb*4J16Q2fb$QX}jv@u#x&oMNhx&hD`wIWS6=(a@?AhzmK?Sa?04Ce7 z{`1O89{)j!&u@F!ExgGq1+J(7CfmJ!4jwIoN{#pNrG_O}CUQd+xV!?GY%lwykM9)p z2PMApJ!oNM%`pl*r2?32XZ!M#mw-R;@%2x|wQC%L*Kw-~V6uH)x4z6H0T}O!+m>e= zyb?50;K>!h7f>bSH{|O5^gk%^{n4ZQS6?(wfhSb}pVxr<4*!$7;_i;`2X~%atiTg1 zfX^YYf4hTOAO64<4|7tS?$lYSz*|-TpC!TaSv9uYHQvXMm396uoCEWKCR6~QLEwN@ z-&2N{8V`Qjk3Bi;(;fxhq5}9d0tX$MC`gtW4}Lz|t#g~gECmjyebGYGDYC>bv0B&v zz{lm^%G%mC$x{$hzi8l-2po9adt|RaaK)?q+>6U=%vTW8zi8kS2prt8d3n1M<6ZHF zowKs`$r=h`3K$K1ob>U@IdOd@#)IEx+D++BxB^ZCqk)ehaJ^y0C;Qzs9{g^F+rEQ| zhZMwAFdFzM0*B1?5^_5k5B}g-{=CJcPDKCTotp;WOn-%lWxd8)Y0ZYesn^0IevaC6C z#vKBKqb}9N9Rfpdn!F8>q1K_<{we4}S97^74{BL+hM?DcirjbpgF6JS&=1QtL=F{z zDyDAE6oOjToH^@Ihro(6clhECfoZBTsNQslb?Bd>r{J_4D6fT0Dz3M5`;t=|d8OCQ zH_hvB(r?d*NN&`1kRfugb!c`gr|}DNpuAt-vXBkUG^y44=#>X>hrr>%hwmC92U&;y z`MY$T(>YMyOs|ky>wuh@UcR>Wq-v%_ztw-X=-91xhRA^e(9kfe`8XlerN4pHy^LNXo&1@ z73u>j{Jwy63G4NIsoXm<=hwekX36mY*TMQ->!Y6sKe9JOhFFE_|HOS5HADz{8I@Y2 z*FOh(>&(^eM-x)?do+dn?>ma1{j5WO>-G)K5`tdhb+WuYvAF2@nmGO5^k7~}Z|1d!N zuKgQ{{>aYjDVMACFhusW4*e5%XvR1p=sDh*h)0cDbUK*$NS|E0cd>+u5eA%HvMep- zKBWJ2jJ@J42=yk5nlsbK?hoZF<`mW+b1-$+Rz!yaP-e**H2WMFFZk-8B_I9&^oD)n zwqW|gIy>~op0jNO-i|fktddpetV@3WLNiU_eV#t8cyndVO^x-(Q$_|EoY4z~GD=pV zf$uj33qpAxD9^NCHjMK=tUnPqdG3(DlMFbUWEC3XdN4)^dIld|X}&#YO`n%-l{~_K zjQ-^GRqYmDK<7+QCdn#vuzC*qFtYwBb%xh)HnZHN_s@3E(VrU9A*0P3e*?}US%n^+ z;us_ZJsm=Ra%TFX%-e>+gFop{w^1ytuxXzGXOOHxe~la&oGk=B#XD#6QY#;|`qX;; z8IP)Q8^`@+z}X{f(BH|MWMV&7Dy<`}!Ef1wW;ZagB3xW|~#nx6{xx zJN<>|p8ZZwFK@tEBI{87@R33@J;vLiWJNFYJ=j#MzsPFqoLq;#PEdx(8uXX<*RI8d zpvk-)%6qM0&2z&e{iPuqxB6$%4h3b0tU_ZPORNxt@-`FiV_w6|zxL`c2hDAs_Vlm; zXNIgnv(w?We|FNm=ky5gd;f#WqrW!k2I*7FEH8B-VX{(P5$XSS0BilykD0N=qS&>K7fGvO#Wiu<+TUg*I6lfbC1SRm^QEKR|`E?9y8#uZ?2 zK*h^0a#e0VOTeYG`AD?>wQ9|XC+CYYIH=-*x#Bn<)0o>IY}Z+JlRr=Y`}_|i%AmlC z2j*((yel`VJS-m>-05Ghq+IiiV{mZA19PQ#^?SV0?8pare%;#D;PiV84zPG&u7a=5 zf*Tcv0^2jVOFy+sA6#(+gM%y{m@EAAC^aJR?_}G-n-7EUEBZPzIMCvOxpF`c*N0=; zXV?zj@UG^-nn5cV9BlExTqUAciDlwz%0sgSe6r6rS%tVgQ1?Eb; ztJZxy5XCk&7u0;b(|+H@fhU>Bak;=;X?U&X&vy^X$BaMwK7U<)@Es;{Y%VZYQeNZC zk9{p7Y@hL$&EKnRazm{s!*hYTQuLa)ikuVYW5(aKF}mxw%QBH;a)G%L_8O}i-l(1r z@b}ku{ALwGwMC z__k*vN9F=^rS`S9ZyU1;*#6@EKA&;0!9yn|@~>QAu0+3fjiLdlnA_Ip*c-IEcCFTL zz(kJ71?EcoJr$JAdSmcuGG{@})6B(8Z z%pE20*u1uQy_Gxw5EjEfN)Y-yqx(g!5`;U2E6082L0W2u2ClhO|(caBk0mM;ohy|aBHX; z(F=MR@NU~Glv&^BgAkPVo}t-wXWh5LKIwfYD{yP5+J=Xy-4@zuTZ8_()#J(tA*f|* zCMW^az`Ja#(EV#d|It=>FO+4^8fP794YmBw>2bItwFQ*G<29C1hs5p<*Y-^E6<%g<%*kE?uWfR40u;<6?)7etd0=W zvNaTfJ~KS&6&H=0S0+b&8D?M$TZjHmS-ks_5Y$qvb=D!V_U}oX4q~x3=vUZQ16xP{ zs(5mGmJrlZtaa8Qv38Yu-zF^9w)@qsgn@MufDVgxnj{Fd5Nn-vNUV+RTG17Ywew3% z?PXxq0?^cZp8sgoyhqAXtaa8QvG%zC3@H|C+t+^kz`&{mpq~t{jtI@vQml2>A+h%3 zrtB72tbO?4NQ8k^3PAU*>Ms$3T8g#KIwaQWw3VWSopng8jb7ub!D4OURi;=2s}O*OwOBk_Xr`88t+Nh^wFm3{7=XpvuvJZ)8(6sj zwCRdDf-jV%uIj8qy6Vg=e>-7a^-=fDiw&$y0Qx**{|2F%S}MTKI-~$6)gFy)7oBLuY+Yn^pStc@9cw$d|bZ&Ke zi@u^!!rUtKvRkF=LeK*iRZ0jNd-?m**^f!*)}7~N2AnLn3T;*VZbCENV^L9qpz9_EMJI5bTSwIWt3dz% zH;b1FrcQ<_vzRCj{_LIvEqU4$+0DukxPTG zB-;+YwLx5?eArk9Z?1S?1qVKM@kix+%=pd_`RZtNJ{51UcwnySm#fG=BFM+~8Q&Wj z|82&}9}M1P@xWZOFxRr~0Y4*b2WND-yNm7b$>5C^56m?gbGa5iY=1t$56{(5&Rja1 z!J91}m{bM(xb$v&vQIw1k9RkVFVsJr!5c0fm}_e0^1*bh3$T4Do<<%z>pSK*!&b-z zb|FtilP4|vT~D?d>O!J<7YTbv3AVOUKru$F+E25orwI?QIU&L|xj{!CMyVaw$LYdEm_V!I<++rj-N zElhuW_$^z;zzx$qA9Rhe9X!H+a!*}j4~8w13tX0f>#b1KxRQ?`WD&go?-ipr+kzi-gcr0!SH^2j`C?K-^%EnS27p4 z6am+m_hmpnq&MB7H^tk(`ty@vopXUpax-39Zm(}XUc8kr6AJ~^>%g!ja)F%*xaQ@7 z@I-{|OR>(g`q>f*T^P1_E^r9~t~PS(mr2RCgX0$LSgiT_nPH3N0v9LXs#Aj!^3&m9 z+fatOR=|Ak%!yt zdonnw$ort@nw7c*ENQVZA1~g4(LX{0S8il*Vvz^tnw(bd(El@eV~tZT4l48G%m)T18hK!@iEE|lnNcV)VDlPZSke3L@HX=qoNVNQ zx#q8~zh>PC$p<*KZ-(x1Hu~5^2}d57YZ_brSGSA@`B>wt+doHVHbYbT&+7j;>%*IY zn7orJH}~%SGuievd%nEOb#HMxgEK!ou!wu_I}eTVjj$d3Vlg`+(>>a(|8q(||2)5~ z{QxXa{v#$-GyjS6|7s{j&r`uq!3q?t!2fG2V8GjYD_Kh)Ft+YJ^yx!SEb_AJmWQ>S zb?C76mI|-PVePbcm3ABOmfi|f|L1m8@dJWTZf$c?_UN$o$-A?zi?`Q=G#&T zYI&~OS%=P5&#E=HC_Y!+qfDI^2E2W@3cdfTN2n0g@>sdE4jn7symaz)e60LlU_T!N z-nv_bZrDDmxFFQxu(q=f9o7!JRQ@MEtUYgmysd%tvA*fjfF++dXJLP?7WLGSMjQVDR8Q9ta&=5)21RZ{U|?$sKszoNyH5ye zu2RZ@etvhXcdf^`O6mBy9;m3Gh5)qNx|^SbpynziGZc1Y{gf>l;VPxqm&X1vu&AEb z$~(#j*1EXxs1Vd#rIZ8R@u6&FiF3G0Y0j)ky9{hK0ci8eV=D_m%~eWfs7RR@SGo8G zT&3hOaMeBoTU7uWzRy!|P0dwGInV(~Cw5Ev;wq&pg=EVNY!w0Mz?V-B3eD79rIZ6b zwf0Tx?Oz7A zq5w4YJp4xjHGh=SzCY$FB{Nh~{Ydo&>p$WurDGf8iW^v00qD(lqoL5RPdZ90U*Vkt zo%+4yjGBqKnb7ba-%l7=7wb@1?)4>=5Ofi(P|SheANI5H!(X_W5c@63-hfll*4|NI zRfq9{pQm%^VQn+iuAXN7hT^YqGog%^ZLSz_8rmv!_>^}46l(C+bP8{#b~C(=OZyvf z-=A^z$IcpX3fd~PYPkHLnesMhmVJMmb*S&pnYNC)UH#-8xb4x2ZM{qeoO-qj)l8I&1);oJO4LlKud(M=J!&R&XY+tH2Ar<5 z0@eTNCasns1hs6^;;chWTHY7g*bFymIsY>JkO8OatU-T&|5*5~5Y)0si?a?jY3Up@ zeG6{VQlk4QtpTU$tU>>bzH=O9>8$Uh&sj8San_+GE%uQ|B5{+J9;xYTlxR=<-`W2k z41(LaZ8re2#^!(euiAn}0sqw%EcpKmR>1ZZVDLVXmxioC(vXeaCs(^4YGDt9EccDw`2wJNU?sh33T61oa%{$pcs7z?;^F zN9SXWe?@-1`(?&m2Ja(z;EDw7voqsiw$}Eg_#^R>ywLhGcrVEVyArVPU(4cax!Vl( z27xoS1nfg{rGGAPQv&uZZtt-+*=DfrSi4H0mt@Zw)-M;>pMdMvX!OY`AK(*jTsmjg zNMzV1xxjwhj1x~!>X45$K7IdnRtmhzu)ev#O$gZIli@{v4zZuBC4c9pe8R9kxxl{M zjGxE&wehiijW5n`({aSqUJR?x1@$#=4yj#5` z*!CG;ja^XoiQ1lFy>o#Zb4$_h(RWFP?cnQI&iQSv8pE){oB2iXvPh5fg2HU14Zp^~wcqz=1t{i*E$BgEJmK@l4Z|X4v|ehH_lrD^6>*WI1C*b;lzenBiu^s#bPG2;c4_upo zy&qm&(G%Fd6z@va`ee$ssSxavfe)W3@gEPK7a4l}e8$GYS z54L^A?>F~4&U8^TIP1#;*CgPErL#kW^05>jE>BtVBf2P#|HHt0puVp@q+X#8N4EfY zs0*vIRF_m+R8gvaswOIpGFy2^nWS8z9H}%Y-IXfG*N*2L4i5I2_NVOE*-y6bX5YZRxZm&A#$Ft(-^g20wu27o>jU7j z3gVLWv?lW%@0?27OE?-yJHwHk_UHJN%)?~ir zopaepT)e-|Z4r?Hz}<52>3*%re9b%OvVpiXNRmb03I^bO9XRci)?~itopV`FTpE5X zY`544fYUYMxvg4L1^P5{F6%fKu)OiUWB~SUhnvE+rt%gpYl%z!!uz^hLzKOQnznMKU0GJU3UH)iI=BwVhxvVBGbziqT z6QTvcbO@F`s4-dI^}dR@)EVyj`#~_6m;t4qG$zZt-lI8>E`4f0i2-BXpnX-1$?~rE zmBhpSYqif#8DK~c@FPTHvb^hk1;^O?qxG!_(C-oW>ZCE1w0ILQ=NLo!R~--mx)lY_ zN@z^Z79PupN3AW&Ub+a-z#a5Eaezg{`~a?fNl(U-&$iTZoyc>dH8|) zue6}!5%8>o#*}m0J7Fv)9yMM_gLW?HrpH0boFyC}fh>WWM*EbD2wA{2#5{|1B86%6CPc z#k3~#&F`Gc9OB}8=AhQW9YAqB{CQPtGT;5qxy&XmJ^@Ycx6vXOcaf@s)?~i@opYH* zT=eZqY#E0h>YqL^yNA|fzW<$bnMqt4X_p-u9t^;JIhBOaO%oxwlh_WyPYD#EL=9}O-mubYM_FjjTHvj-rLZR}a#*}jxJb7BC5*PO= zBe!H_fbh2vw9=SzE`ui?Q;0{c9c#-aM}Uzm`0B4Q<=h5OJWRx+Mx%?L$^y`D5O_F5 zW6HS>o_Mgtqo%mw95w>@FZ93ftv}*!R_-fdkTj7bvc;ZoI_S$EG z$)K$xNY2!-IoHAyk8zwwt=|<<%pZ6F9CFjJV=d+}mUvWdSbD)_^d-3!#CO-Q;T9g@ z#G{h`wBAzz2#5w7=WE!Uo8igdGKTXQ-t^#!2oT^1Vx<~(v;|`{@o+u7u*=#Apnn7w zE!D6&cf%9LDB|JLRlKK9GSCNrC55%@NXpw~) zMF5aq70!94WzCnxbCkiv#Y?xXjT?FkpKb}q$7@-O+v4FM;!@ir& zns1BeCdo8$ywWwYV@I_9Gtd`xB;Livc4g;G0In=3E$0 zo~ORVqh{;nPwz*7kQLzZ4-IQ^V?69bJgVJve7h(a^w|m?1#4J~8{=Vb;!(Ysv|Moj zf?YwzCJl@F$fMu|01^3HdT|WhMd{H9&>eu=Q#EYPmGQ(Qn0QpJ&N{zJ2HlH;d&e|v z&YkhZqbKL_=1qbEeOt$ZdjT3Y=hArM(SvwYagCiV_5tnX;L1V`n{#VC@#s!G+_rRx zor^vz8iSKE4V!aqJn`tpG0NV2zAYKFxe1cdXI9R=@x-Gm@u>WB$*_I^1g3(+cQkCy z#qq?W3-PG5YU-8|5ujxd*tA;1cD7iL&YXwa(GQ8}v*I;~yP#o>79K{B0(t)@%bZR())dFzxCLE0()%F%H?TL%Wr{&Sd5G8gr z>@`};8Z2B4#HIG~OIwbir+Ib{SV693+gZ4@<6O#jtFtu*fWJGy^2N1mkcCSSadB_x zIX(c*CDI3csHkDl2l{__(%KT2n(6mWbOvDLOYnNShHYcv(T4L_^61Sb6bP3A8EZ9c zYYUIo#G}f`g;DR4L5K&quGO%C79N4bqw*WqjQ+_0J;-PGXxLU39<7K+C7nyRWFOFV zH8_sGw^~|wv?Ly`FJ|8Ai=NgyAZfgYZDHZjg7c8}9y`zn1Wf^(i)dK%wecT*+W_L> zQs#GuvOb{o2e9^=hHY-)(VSxhj4g(a3$;ECHcnEa*Qh24`LEyqXD8V#X)FFDUM_Ad zdMt_()qq#w7;p@hL2zJX`#9m!g^JZ=3IG&9CX?BbsW{ocXn?i@y@&N3>(?)IWWN#p z`gSGdE9T0O?)sd1xI>k9zh*Y-vf^di&-EIsKT-dP!+6y;Bil;=-J#v>6?vm`3fL(6 z!zlk%?OXQIpK6skZ%19Dk?k#jb}18G0&SJ@m~~FU9Yw#NbAIf&otN|{eSi0O`Kp_d z4YrO}r!9<6&3o246>=2)vA5^Zxbt1~r@OZpu(;h1BiqwD+G*6$683qcbL!+M`cv~; zT~cSq=+A^Se0{LY&-F&On*e%4;Ro`((K(fKX7tan5u;lS(4U)9 z)1y-Ta3kAQ0NuJ(nM--2bL!_P`s>sFp;ONd*PmaY*;(VOzme@CfR@j%+$V2zP8A(R zf4ekkUDgw4{e?BEKcbV5jcjM@XjOL0kIs3cb86`*`uq0Iqh9Uv(qCN5wBORbuaPxc zM=M))leNhkol{Xq(LWZ~diniVW&NcsTiu&akQv!d0_dxY9?i-dol{pw(LcvspVzPC z68+`fYwLB7MGqfi9qqWZll_6b(K*$16rJ6(mS>x$sJl)N$0Oap&M>eY3A!ZzlTY+I ztPVo|#|xXxb$i|Q=IT0H2dKMsKWTjJ-fdBJ-7K-oSyZvuK>*sdjmttIsJXh%482qM z*wpJTKB&5G^8N-7Q3ts80?-YUHKzojw7Sj=z1uCJUfa{2xVo-mW3G4QRhb z5K5c1nW6U&Y^>6(_+)%lU+12yTN~JR)}dm@`o^7tQ1nfQe#}*%x{SqPJF}|?;?o~v z>eg;!V1ukfMW<@X(5Hm;7s`D4!#vZ6H|`B9W^!+s(mFl;#4i#nWGL096($gR2nLZgZbDim66h8eiT3If^z_u2E zie1A5e|_`m4>R=XslYo?_xs_~A5dBXwd)GB4u!XioI*Pf>uYL0{b7baca@HO`gk}# z{V`=|`8o!+m363ysWuR0QUsv%^oJRmX<8pVzf~kY{n2oi>b`+(DFEG3^;T6OsJRN% z41Mvsw#&=oU2zrY*x*qM3~UPlXng(4o>K@|xd9i~SgEY_fh+c8-4cSDyFZ#|`o810U9Xr0xcg(RHx1CC0#sOI6*{Bv z-ra&wdR3np`k_Q=&FQR2-2L(NmeYCzE~>E#JvDj5BO$1{`(qCD#?IA`T~FcekEVMq z6$V^TV+|^j{f_G(1T}YmG($fwleYc+`U&p-xaviGrU4hzScAeWw|53XXgU1Zr^)${ z{?z60tf2|7aQDX*%}3ob;6fT}P|>>1RnUIR`g{MRWdXeIbJ>K^OP(a)?vF)2NM9Oo z5sg)7?f!O;g`no{kLF+hbH<)|YfK)v`{Sq;o=Xk5fW|7cx2_J_;9H-mx%;CT`epr= zotp;U#oZs@mOPI-6{6x9YfyNtY@PZ-P;>W3GxTf7urv_S4|jj`@SAhcfD31=LYHS9 zo-PD6cYicPzqvgsb))J6-2E}~x)*D}MKe~R-J61wf>7H1F$enc$%F?TvT*muc8?4h z23#;>6}rEhc7YJo-2KrE{hm0%(KnU8~| zj){H@J17^JOX0WO?RcsL%Db6!O4YUKpStt@s1NGh`s({l*^LuVe~sDmX0|v}xBj8N zNXD}J3_CCvm`nTzzJE7AHrRIXrsBI(md~Edumf^|Ik}+a!loH*vup<^cq{E2AMeYs z{d0jik)hf60Y`_&*bd&+R&_qsZy|$250r6K$mP{QHJMB#1CYzH`b?gB-qFc+@Vv!G z3cW5>hQZ+n56q=z_38Cq*+tk6Uecs!%IxH-3=TnfU@oJpZ`ZKO^?bZ|t7=>vl;BW{ z!C?pw%q4=otGq0_!^id+uRC<>+3S}l7#xc5z+C>=yT+zzKa*_-$2IQgJ72kZ%j;U6&Bq#VUGt?*(`k_m4oP@mE?eDT!;BA|^RW~=qK*zY)9Vz2!xA2tOLBYd z>evUBY4F~|ygyobOYd>+@pHqG-3$&*cwjC!UVlscvlsaQC%ql@$0bI~;P8Y8ZbmZQ zo=amJY4QO+d_;48b#Dd3hU5ZsqM(1Q@A{XCwlBp#&rL@{!z2vbFBh1T5}WLNzjSnv z?cjq?qLz-=85p*2E-%A=7!KbT?y>mOn#9)nr2j*nbCQhopXY%ploo(W9w0p)n z2FnOMFsH8iCe@B_nri!uFQoK2RMPPzgGB@$m=k52tV=I01-64zi$cE#7soMJLg0Zp zE!X${rm4mA@l;&7wENxxsA8~yzyos$7a7%Up_z?_)g@L20hUz2SIKZ57_onCl^!IA+F%xUoT>wfGwJ0ENO_<+xD`H~V0 z77Tb`PP+FpE|_~hA8Y)yYs}Y}i#Hi88t}lJ;@^6}gKY`<_^dFvN!&*+N2i;RaKHm| zbpkElwQn8x4fq8nk-yaR*ITPNFD8V*fElyk{C%-@k{Ysv4`lkXp+bUrof@#7MK9*q2TrDN1U)p#m;SvVSvtn9&3U&_m;^@GC5Mv$m`piq)aYq3~bfQKWLH%P?=F z2?(H7?{7ZLd)AiqPO&;v?-ZT-$w=y*N?q-2WG4%tn`bD*d7~|hl45nZC`ohjATCPk z|1Ef+k)0%fE><|GecouxqNG?IE=ux?o4b}t4u1?z zFtQ`9qswn({C$!w(Y$iZ95I#FJ01VAzVwW^@A}10&ki`t;US8rwy zd=TfU-*_ml(aZDujqET1bZ}Dn`^kCE`j}a2rRW=x$BL#j?5dAX9W)_uCi>(HwT>vCi&s|B{n_oE z@?PuxJfCQNeQxwR2l|@VxO|Kc>g+W$=c9 z2j&isGVYDLECsfMm!$tJ{qgaZAQLKv5S9F=Py&4d&JO7{a7e$tJ}sNXHZBxoGnz3u z$6)wU`*0n;OPV$;Qh(3h`a6Jv)kXFbWHEWnCxb9%hD2G@?Dxu?gbmHAkTtNVsOOF8$2!5KXsnCoKDtns+R%kuHK-7F^R z69?WgIG@J@bIlL@D!;jWKx_LN-~O}mz=MnX7@W=Hfw`Uuz8RXT4P$Hv-#b=%eUBsr zM!7s5m}|A*=em0n3Z!`jExbQksc(5$!v2W#RvLpdc|0)Jp}|in_Y%*u9sD?8vg4tq zJ`B#|@xWYTho*ldm4`&w4t}~|$>v{KNM%P^JRX?q2hptPgZbUFYzJr7blj+Ti$3B} z4vz=s+DJ65*g7peA1~g^tKaL$zMv~;Q3j6(=DJTbsn)l}+(g@F{O0SBjU!yXF*tw6 z19MF)nlyRP=uVjJ;H;UEi=MVd6*(w-#{+Y{E_|=2FX`fLJNQGdD($aUjb(7|jtAyi zWcXCC_NLA}+rgivxQ%;NF@nLFJ06(ptf6lq9v6_0$L&kkRr-d`s9Pk;+ws6$Lk|7r zvaf*F_8EWs`(XCcKd7h!W$k!iu8&8z@z+Bgmfy>>-{>jHu1bhta$eek};Cvkq%=H@SceVYFmic%p zI+dH{y}1Lr{SjsBcwnwWNl5(a5-ai%gY;Q@SXOih$`AqH|58-O;EWt^W?a{k&MmX= zg(lhl>W9B6vVWR;8^)x}1)f9p4;>pUSm8XtcJPw1Ek|52b~Wq&))al)Y69}3|I1LK zsNkny1qxQ+fAtC&abCe%ezN>?>6!V-mI+&)3T7!p#Oja`p~)Cq6AKZ4mz=ZDh(GwO zp$jdXd-pUiVar2X3K6k7Bt*39aNr~sB9@+h@W+T_bgSq#ewB;xvZ{PE_j=JZbF2;t z5ewHYDT;-NzU6vUHsV;_D!M{p>5RN*Z7D>=>W~mICMMw{79x5+_&ULeqjhU&r{_Pt z4(B~#R5K%v*R7!ol?Z7(EN`@>5D}|GLd2>~+aF;eV! zhg!?A5HX~4Qh*Uh?N-stLuUM6Dnu-Wh*%vGBK8~I`d}fV=H$ydBaYjxqAPSt&3EXcKH0#y?EL^ z^vrZ$^zeCxRBEj@np_3Yqn;X3rDtws#?0@oH$Qx4^y}$us+OAe7Cn4RUvE5aFq&Kh z(AE4Z4$B*De)!DjciN}{3nIRuhp*YWgBziys#Q^G2H=J~R4*u5p(qH9DY&Z^GQrZ3?4FYaLyB?}7zZh2wwL z+Wdr=(VyN#Yc<8o;wQv-e>r0`X{@8mY(I9mNZzwHOI2p{m*!(Vr@4?60EN`$K7gcrt_~VMv zRMI-S^!OR#;Jg>aTvTO7|2a|gsV-qEDynKT>8-Sv(d29$U8Z}BdQo}L+FVpc(cWM| zt-vlhO|6FKJ!^AOl^G3ni)6)qLsV2%<3RQV^p+_ufPP@Q z7@GI2&6Q7Pv}ncsq#jM@;>xGx2gZCdnu=LRm%4sr>w>(|L`X+h8=BGLCdY5>s_cs@ zpE_n0{@KN7Dk^{;`A1@(_kx%!pK{Re9`D*ayaK9x+H$Q%*a)Mkh;?*@`fkNX@@CEZ z0<)|kiq)YSqJZXgpWzy!=)MkqMl1zd`);ftS-I&bZ`M4tWeriR4%HB?A1qmbYlwaa z4eoBlVxU#D=61`!^PY81a|WXS$LdfG(euMK9dHd%%xJF~Ml1(fMY}u*>X&T!M(54i zvW6&DhiZtXN0tP*hNySbZhsrGAZQI;KCM`3_$85A3Je z*HBznj8{~UpOuHnOURDM`pF!{>%=WYSt17R#w^$siV%9gC_#k!`;nh^)c$LXC<}i1 z*N_S*|9e(IZ71pezh}1ppUjlOyC+`Oa}`NWHtW&iNpzO&$;pcok6mip?JG`EK8& zFx!_RO)3i?tCcW#*Te%a=hnDkax;f~fbX1xQD5gZVDO%a2VTal@mi+80(N_a9U?~hhg)5pKI zb~~Fsh%wd81>Q}-UQ2dnR?mkLvazUgTaB-F}ddg@8H1x)f()~$BYkr0GEPGpJ7b3a)A>G zxc1^7^)CS1XMFTvaI%ABA!DkU3!Ff}wRS3(Jjn<6c<`g=dp}oZOf_bbz%xTWZEb>{B~+m|Ax#^dwXmwjbS)pCKi5^#;RVIE`iu@vY3 zTt4?V+mbO=%>~{1l+4q_3I_$Z3ho@dsXz2 z<08i7mJ7Uz1FPE?P0z=RH?qsQaHnPHzjdq5eBCNdEku%+`&8X za@aE4!5c<=eH2i0HiMIaJn#kr4k|tJ=FEH){B22C+fFkdRgR%lAP>BroAJ7Bb#KqJ zea1TqZQAAg1J!JzL?92mj$32VwmR$FYzOZNedqT6hAV^9fIRS80&af2d6k#>c<~O_ zb3Hcdu@{4rfIM&v0XIt?^XN`Kz$bbp{nS-J1zsox$OEq-;DGU6S02d6d+zjz==4d& z-!nJ?$OEtDz?0nPp-&LrH$Crc*?9~~ZZ{f{wF_~f0$D5L} zvuvN?)2|n%uKKu~!PyobcoFwvy$noU=VLqg`5pF3=MJ+`5(qp~qylxe`inYC{Zjn| z-2|AXzN)^UKCMnxA5iaBZ&$~w*Quk`OVsn!Gt^Vmlh7@JBh*9GA?jdt7j*}98+CxX ziCU+wr*>CYRlBOosh!n@)sAYZ8mNA#KB?ZSGF6XMcU3o4sj73Tld2=CB-KuoTZmP~ zs8*nk19Q+lgRE+TYP2d$HAvN0)k9@e8B~F)W-1?5LzRcBrpirKK~+XoT;-&)S4ot= zmEV*fl&_V~lo{w=!t2UQ$`s{sWukJka)WZUa+z|0a+Y$MGD0~H-B}o_9H8v2 z?51RtLCTg&e`RB(m$Htsy0Vf|t1P7~s#GcEO5pg-G0QR2F~c#Or$6!aJV~}HjqmQGPqr0PLG}UmKK5Sr z?)Gl>T6<@ECwsX)P<&HlDzX$AiZpcJVu~VJk)%je#4BPH%M|kzQHlsfxFSq3KoP7k zDuU2Wj6Mo4g}cH{p;b65oD^~ekbjeB$us2{@-%s>JO$m+m?Te>$ID~n%jEOqQSu0R zILh=4kO#|+@*sJD+y~v>=q`7YYvs;zC%Ie>WZz_2vP@ZqEKQawOF?HtlF+@5@v<1% zGTA&?lq>?(mxRd%$bw}?SrEGG(MRSbbCC0b}4qrc1d=LD3cXqx6E#yU6fsfUASGC-2ilxq|q+OF2K&m&dbi-&dpA1=WOR> zC$|HVZ;~uYChB{cCP|f~NRlN7>T9t~GEWjEiI9X#!XyJE!4f07eKJ7egRc8_ zm$*r^5@(5%L@oj1Z{jR*rZ_{KCQcQnh?B)h;zV)0I7YloJWm`Yju3~7kv@ihq`|}# z3UL3^+sVX>gj@ zQ{WV_C(+ezcwJ6_6T}_|$B8`#juD#-l8HSEjuLwW93l2FI85vzaERE0;2^OFzyV_S zgZ;$r14+atfqlg81$&9z1NIQR8|)@_7uZGYPOy{M9bgBsi6D{K1du@NcCek;ZD1R* zTftUhw}36gZU&o)-2^rf8xP`%jRSGS#)4R4H-L@AZUh^MT@Th1yAG@)b}d*-Yz&AY zb`4lV>}s%@*i~Q^vC$x!*p*-$KJj1*sB*FgM;zdvnL*V^uS~H?s)9h z4Ub*B;;~B?Ja+DkN23vsojT!>Ver_oBOW_+z+?ONcr+OB*sdKOgM#qbwk;mpw83NR z)_4pI#AB;gcx>4ck1bl@F(3es&70$~Su;E~ZHh;Ke?0p6;ju{*Jo@_L(Z>gmdOaQ+ zH^!s4Hy(95JT_{C$A%5@*q{L(y}a;Pzdj!8)x)EwCmuaK@L0Dl9_!S>W9`~_ba%&N zty*}jSrd;nYT&VYbv#z9hR3Q^@mQq_9^KsVSh+GDD^CDT2qs zh4EOZ5FVYJ@TgYfQKiD8Qi(@LM?5+>;L+Y5j|v4I<#IftOZLzxmEzIP4v!KE9>roj zibP_uLV@CcSV$HQ)Zb9_pQ+A3v45&M1x5Zz>O>Ux$EcT~s6R>_fnxqJ^#By{8`VK5 z-uF>^p=jSttwpiElUj}<{cox)6z6BC(omG2qDn?FexfQKMfl59^H6*rp$bRQ{Qy-k zitU3`0VuNfQn{nJUaN9OQN3IRP)wht%tR4=nlcr|^U2C26wSvgV^AzVPZ@0cPXXP^i^_5T{b2cQVv>t7zfXQ2o_}U9y*iY~iu^-__Vn4tS#J-2`iOqsp z#J+>?hTbH^PUI#EmcmlDH8*fFy2&_aTWJ;XO#=MtB#JxDlp95;wv-ki?DfHY9N) zOoJqDgts7x8{thz;zoD_lDH9Gha_%<*C2@-(T&0+ZiH7Li5p=mByl6W3`yJwFF_JF z!i$i^jqm~_aU(nrN!$p}K@vB@vyjA%s7o=48{rv9;zoEHlDH9`f+TK)Cn1R&;R$$x ztj}>s;zoE3lDH8jLlQTlyN*cQ2#-J#H^Re^#EtL}Byl4=2ua)s4?q$(!u^oMjW7w4 zxDoDyByNOzA&DE|9!TOwxEqqV5$=K{ZiG7_i5uY#Na9AA2ua)s6CjBj;dV&kMz{@< zxDoeWBJ*0FBSqDkn z2-iXqH^LZ5;zqazlDHABh9qu;t00LRVKgLhBU}kd+z3}d5;wx-ki?B}86!ugQIjc^_$aU+}yN!$qMKoU2?*^tDIa26zS zBb*6I+z4kt5;wxXA&DDd6eMvYoDNCc2&X|3H^Qlq#EozYByl4&K@<7?SxDkW7zs(- z2qRzw`S+9IWMU`5NyJWs6N#MwClEUxjwf~;97pU}IF{IO7*6aMIEL8Ka5S-_;3#58 z!jZ)O1^*&;1RO!^a5$XUFc?PcFgT3ZP#8+=P&ky>A#ez>gW+Ie2f;za4uk`V9RLRq z+aLBPHUx$c+Yk05wlC~UY#-Q%*xs-=vAtj~VuN8Yu{~i=Vtc?I#CC_>iR}iv5!)4Z zCAJIfLTqQ)nOGw<659!OB9?&+u^nMYVmrVN#I}d+i8Vk2vF%_xVuN51v29^nV%xwr z#I}a5i4BB-#I}O1h;0d5659f{AT|I75ZfF!C$<@EMr>2qlvsc0PplvGBen@_LaZ8?lvP zWnwGAO2k%#6^V6)uEe@P7h)^G3dEL&<%!ioEwLJ?A+{VWM{HSGme?||46&tQX<|#k zQpA>oC5d&0&cv30C5SB!ixXQ6)w<*F$D*((u|;4JVhh8<#1=v!EB^0J(1}R$ zDq@vTNvtDuB-Q~s5Ni+ZiB&)av2rLUR)#`l%#%VXv3AgoSP7I6D~4iXMNmX6M6osY z2M}Nj{(wKk{szB^{RMszn+>vw{Rw^&`vd$S_B;4a>^Ja@*stI#v0uO!Vn2h=#C`&w zi2VpY68iysAoe|YPiz**BK944N9{r@!8|4&i<|0LD_Pf-2;IMx4;QT;!e>i!p=>M+&*80Ch ztbj!RmlOG4M&y4fk^k+8{4b&MzeEfbVyyocmW2USpB+(913+r{|7T-(H;Wqnf9DvU z*029}Ze0r41uIan0tGAZ|LzLleJIlXxqSa~TBn@}o|Aw7mFoXrsQ&+%>i?gp{{NBc z{~xIS|DNjqSycajNA>@=RR4cN_5asY|9?gG|Cdz%e?j&COsfAsr~3ahs{cQw`u`KE z|39Ys|0AmZKcxD92G#!`Q2qZt)&K8N{r@i2|I?}de~0S-x2gW0M)m(&RR6z8_5T}G z|G!T4|7%qLze@H0D^&karTYJ6s{dc2`u|0$|6ic`|9PtapQHN!S*ri1Q2qZ5)&EaZ z{r?oz|4&l={{+?lk5m2s7}fuiss4YI>i7e=F7hw^03mGu8h$QT;!j z>i=<6|Bt2m|3<3+Z=m}BdaD1gqx%0^s{hAO{eKPB|5sD}e-+jLqpAMClIs5}sQ$m4 z>i^5A{=by!|4XR;znJR(i>Us;km~;nsQy2n>i_eo{y&%M|8uDRKbz|Rv#9<*lj{F7 zsQ&*q)&HZY{y&}S|I?`cKb7kLQ>gxLqWV8e_5Vn!|3^^$e=^noCsF->BGvyVQ2l>A z)&Iv){eLXg|HG;NKZff6qpAKsit7I(ss8^L)&ECO{eL*s|HG*MKaA@Cp;Z4LO7;IC zRR14L_5VRs{~t*8{{dA0?@#sr5UT(8qxyees{i+)`hRb#|M#N$e=ybmds6+s2i5<( zQ~kdi)&ILv{l5#<|2tFt-$?cUPE`MAsQ%xP>i-?6{@A_5Y?+|M#c*zaQ2An^67Vm+JpMRR7mg z{l78Q|GlaHucP{ZBdY&5r22mYs{eaY{l7ld|Lalx-;?V99#sFYOZERcRR6C{^?!G& z|JS1We@&|Y*P!};b*leYqxyeUs{dD^`oA01|0`4dzY^8|D^mU6mFoX4RR6C)_5boz z|JPFeUqki(a#a5>OZERURR1qc_5V^-|1U}Pe`l)ym!SH8ajO3pqxyeQs{a?E`hQ`n z{}-bAzZ2E})l~miQT<;@^?ygI|2t6q-=6CKC{QNfyK<`k%c%a3n&6RtZ%6fi3Dy6_ zRR0%I{U37rKLk|&|3mfv-&FtqMfLw|s{j9_`u`8A|9_|Y|2L}tf2I2W7pnh%ruzRU zs{ena`u_*2|G%gDe-_pM-%s`{uh|DbZYa*VR0vaZs}@r~m-$2i9+j=dZkJCCfr(_hp@MCS65xYX+ZzWClm&^Hb2y`nLV zu<-amJgT2tJ-%ZS=zIZeJ*+hixA<$`6OZ=kGpao)4MZNNB_0p8rZ5YaEaK8=(BIvQ zBFXo7Dr~3Mnub}pydy64Zo2)P)d7Glu2556V+ysH%3I>%-lKW=hafQg3pn>pYZ__+ zc|$w`H!Q0$Z8Lz&!bSV+wWc8!F0YA;-raptxY8ul(>cTDU3*u7!itCnXa?sZqY-_7A4Y2UY zBp#LA)^w?m30hYIbM9+R{Vo2Q=ftC1lJmUi!$4eOhTUEC`-E7yJR>eGo{m|w`WS$H zqeU^rw5EO*E>DR|{byfX53>NQ-VOE7G^V~5Q+YyMYERvH_w6_^pep#J)R_8McswQ^ z?)91!IZ*+0aRhJnXiU8=JRT8`>cx-lnKA)1p94};G^Sn_9uGN>UxQ8!(1K>$!Id8RG!-NNM_aq)Ay{H9J8fNFP9QD3d8n}y3=;?ktvzEa1t0Qe9H<3qKk zt`;uo#6|zuZqo$x&c1U2M$FTi(Bk}ucl#aU;`JnK!*neG+u~tmKdq^=h0AT?Qg30M zA@0EdEQyBYG+L9!7o?uqhn70Qv*J4eUFlpT^XY{(sz~e1k&|e_w150297}haEMh z4i+BQiASCE+lObyfYG189aoL1y@kg$;^BUDz=7ZhFnkcw8YKHC{CG^N#?%`hq=LjVZ{&Bb9hm3tg`|20+&hAmNPG)Rua1PuOMR z(dx&|=6Lc}PcAIUOXfCC?O3f@r9G6FU zm?pOGoT2IWRkORItK@Gn;+Q<5%Ifp|P#@)Rqh@^tSIN&}#8G)fl#1KcFzuMzO{?3G7!7#($?$W^qNtkKrl=PEU|FgYv_ zsW3an^PNT$NNbW@`A(XDU%?vX%T@BRc-SHNVdb)JI&VtQv<%WrsmO1!`Ybqm!&Nd_j5r{VD1WzZlT}`t7WXvMm+?O|^2sJw$=hP&L70NSuE^mG zP4o7eNgZ6I#>xl>G;v}djHt9VbF!Dl_g>rw!Kf*%{6snEDm76aE$`<0;n7K-?i9_? zH1pI}3PeTyMGp8HL@78M;_v{d-d28=#9l6aVR{>P4cuDvRfYF)zH6*`zG|3SU^mF?Ity_7?Lgz zsW)>;dUx~}t?Nd(#c}9|sG`}av9H%3IxSFIcSzSq$C8^#v-TUhC(TdY|BonyY)2w{ zTVIWA{h^IN6Q#D+J0g_To!#5WOVc~FXRoHcqK5Vw(kr}+(ii7fkE0}L&D6NCbJWs+k&Z-2kXwzbmZNU=V?eg1UU9qy(L$8Y$(F7mr51zNW*U1!Q4 zH#IhAt@JojtWUn1SEx&BEz_n}{Xg@EqeZEeb?Z`oE0?dZJ!_@Mkz##1pju1ag(arV zN4~8XD29kqfOTu(j+L``BJ+4~4Jo|eXOc^t?FH9{w zSd{#%TRTtkp4Z3rtd&MYiuH|~w;qqyCYlb8-sbC9WVa~!TDLB_>SN@2+p|_05h>QU zB9qPKla`qd#b(s-*o+#-HMMSC;!b)`U)$D7BO=B6_RA+e&J}l>4iD%yl6OLH8J}FN zW6!&xEy>o?xb)e0m-Kee8!OgXO@rTPZk{zA84%;{xw5w?dFNvNbY!E^wyl*OM~d~` z0-qZ(p1n**2Oc=+9)_MquUxD@T#6cFdqI>QNB>%1es<*S@(-qCeac40g#?R|$+~s% z8O_)4vTd#OI8vnymtKJiFXi&BGJtldYKNwsaQH0e{U z&o3EJigJBm+IF{0t9w_`tm|91cAofRO@M7{%e~Rbc;Q*%W?SYtSz?! zC*$1$pT4C`-iWt?!6Q~R7I6mtqr40Mc^GxlTX{KVs$Rq$L9AJqPU)BUIKk$ul}^=) z_5J--Gs`SrjXG63o|#(i4-xkSv1&c=kLle1U=309zia&{XwWQ?J%~D0KM^a$`{MY& znkGW-|Hpse7xPWIrxqRmrZBbWYv}&@iqu&uhAL$^Z0K(&h;Dc9j;i(>bV1}dav9kI z-JL#P+r;sK<5b7m4%Z#VI8-7^YmT6Md4?ft<;yJf$kD?9lGZ-Q);##;ZhsGh=R=N+ zGtNsjOI>pqLzX&tOP}PCeeE~)ZD7n!UT-u@Avug8w;jAaxJQLNH1%ln!QsY*@>9%G z_Z-HMosT?lb6-M`{WI1$?eD|vT4?^fpUfau`Dn zzume4KbjltpK-^V`vw%4*48X_&S4Ch|A3$`Kf9#aZyf5ftF=$~GPBevhcVO(nxC2S z_p>bfjR)R(*rVEMKeHs}For5aGpE^sdFk;sblRZU%l_ldQpX&|P@C`{Fm&g{JS;_Q zy-Z{3G}O1roWmF@8a}1lysG16|57CE`H@oS)+VzQoWmIEApQeKUi_zq{l?>@htZ{e zx|^k-9L7*6@#&zsO6Ot5NeS6=!j6ZVr4Biap}yjMrB>0Md1wGMOPfl>E2Xh%jfj8D~&c#!x#lCfzflcG{}OJgRT4^yyh16hAR= zjoys&rz&G8Idu+ip{bpRHLkL_;O*P%_nL9`RAmf>s`mQptJ%l)&$!W@VbO!1Mw@Z& zRAmh1t@hQqiym$5H*OJe_DSHk1T)T@s*Is1)-0p;`$c(p9h(bIE)snJtue}*s*ItO zW@_3XCMge3UbI`O;VGYWW}G!u8AAcit6#*U7q9JKir6aU94Yj`la?VoA(8z){{v!>hi9%h^`RT)E#uxjI*UGV<;`wiPEPP$itI2ztlhfWCo(vTPRnmGKPv|?bI*b=jCC>i%ZlX zE6+S>#+g!;F%&Xub_xGF$<6*XPWy2=?f&zUW}GKg8AEZi+KIJQvh%PM%fD{jO1M3wP;k#u-wT zF;rfwzdZ8lQ66TzDcF6))l!wsI6tZ~hJvj}!{^;Ar~tZ{lsxyBjkv1XhZRT)FMy=U)fHDmM8 z8nf)DmBqUFE%ujCUKE`n`FEE_h3q}}9EobzRDGLs6E)*3sA`5#Gj1NYczvNfeE#;$ zSP}0q;jJ0xKvl+2_4b|c;kZv8o~-CsX+w%$k2m8CsLB{h*k1R(?Hw`5{x`&^y4T$T zJiD23{!?WPMR+eauL-a6uoUAXF19}7A7{qdPnGcwx#;fsz`Np`JS@dReVBKx0Q6~w za-S;W?O@#Dw-xUa>|cts^#$6G54*#$J2jOw-%E_a&)E#OZRoQUX5N#-~-=`1L>vijO?Z}VhBGOCyKs#CM;dse$grkeYVTZmB9J-~p zBQpFy@FhyWtJ5kvc};6Q8D*M&-b&)_zwDS&j#my5T=?s_Ux;S0X~UVDbZBpszi+hA_d%0&AhQgi$5RF7p0-OSUYbm88OkOwUS|1tnYtyeAl;s z57WUJ2fuoMj1;A4>(&)lJuT|{*`~FUP*$v;-s8)rbZurjHulNgQGO1hG{m}f`4+7k z{bnsS``;@>NzE$OFX!z%xH`U*>10al$9wp!;=##j;9Q z7**-kEWIn9Lsh!OZCR<;MQM7ANDSlvx>t1UGboh1}J7)MHWVXWiOmC-JG1 z>UFJJYs<4BPR8qDwYT3WH58u(85OXgk0|xEZe2QR>iYz1pV$9Ajh3e>oQ$_$KF)~l z9fD6){FXN0hbZ;2Ze6v1l_KRn+f3dc`BjcrlJ}>FzdK!aP^4+@gA@B#K3XD5y{%hU zsdzxU>b*_tgQ(6Qd$;+&$@{D_ZQCSu^EItG-*0C}&va4hW!>5>C~Iun;x?_7H@Y(G zcXI~RTo+x_wBO~^$^zb#MJYTN>pIDsi_WlVtrS}o>mP;NT-!1FuIX6U=I4i%crQvl zty{a!OmA7TfK6+qZ28N&G4Y&hA*Jm${WZ(gyOm3}D1}+KcA4EauEPnN)=HI5vDPkj z4ag?;n9glmsvTP>R+K_>vCee%AM9n*TG3O`R9N^$#;`K6SCs$5kmJF6dab z4!_(XT4U2%d0i{kSDa3bd%0nRX={gPEzfmBN7YbC$EtPp0&P8PFNm^%Db`oNt~@pM zLow5~r$>T+d`5=|x>&ca5K5maYuj2mu%TFAf0c0Hb-y^%j-;w<(|pmvQB)+cHtPyC zPCW~;ZLNG>E7o^3w4OnFudsSuciZi9ZFojK3Bwb7OU2=_e<++|G>)Lr&wp-SQ*vp zLS0iv+6S`3gC~6Q|KbNmMH*|fDt}l^zG(XgR!*QOvwB$ZOZ=G3+onU)7QcDBfERH! z#hSHy&1Z{8c-dS_<)d7&ew;Axw9|QX0%cs%N>5#9!Q+4N1EW%m)mfGP)J}c$SNSOY z^-pW%@QO04Cp%~KtP>n%Iy&)V@WyxOl##|cm&jJcB`d3cYp)I~254-qh4QZx>W1IAhK1NOeU5ogL={rq>EDfsBnrg6w74^^8a`K3goXq z{tD!;K>iBkufYH16=;I<*Z)(oZs?4PNR;fBUJ~rdLn5tj_7YGDXlJH}L z{RNB01AjglpWWSz<8YNR{Bpbx`_$R@mi@+0{T}RGP?FG?abjKN{7FvQ_lY=0+%1`) z`e9YyDnEOZO~*c^ZGSP%@v*V!kQBE{f6p>YD6OtChD#{wL`02fV1HGp@$jG{+usi< zW5!8!l`&jh(PH0$F4t-M*P>76i&<2whi05oR~f^F7hVIC6Nfys-*{5NxRc{&#h7tI zU1bc{WHc(!va>M{mA8~m<0G@${V?Nny2==C*JwB}qCI=n{u!^T=a<-I*hDi^q_TpLm^v*g7;HrYSpX&wVYZf_rq-g+q2t#XFDMB0+g zPA^Td-?C-u&@3v$lG%VMpkwpaO8hs-#|tuls7QfkJZ+qp0g&)$wv zovYmINSbkiTV)J)tMqxYE_lYW4Alpv>RT(Ho2~Vq4&2@Frx~ZWRmN~5i#e<1(4>tC zDks&qR*VyOlo(=i$ui^Qw#pdpacS0OOM$&n_CI+UL4%#L_PUsHYFlLtx4-zbZnX0G z5BrTT_a5mqq6j)zfD+p(W4JV?Y2r$;fysX3dzE)TX*ql_jW%{iTBi}d3!j9y!VBRs z>i>6JxF%c@&IzZ4qe6zTN7x~37S^Nuz%n64SRl+6l7vaZIAOF9E5rzcgb1Ox5GsTS z9fkHnpwL`sDl`!q3Uvez;Wxoua1)$`qC!D|74(7w|AYU`zvEx>Px$-%9sW9hnLp3} z#UJAj@_YH6{1$!#zlLAVr}D}C9DX`KnIF$5@NxVwelQ=&_u<3%uDr;1;9K)8cwgSc zd-8Sp8hjPL0$+wN#TVlX@f>g99l3wFFWh_X75Ytv2V53+gZrDiz@6rfbBDNnTspUv z+sOUFt>Bh$3%R-63~maS$R%*G+)yrx>&u04A)J|O%eCZuxh7l#t`=97tH8N&B{?UK z=M0{SDo|_=0{&-=VLfTNqE# zhv~iac61x#YC4Thre~pB8OPD_^e~je=}U*wo#_sAD|Aa^6S_WKgRV@MrJZRfnxl2- z*2d4&8|o>Qjc#wejDF+eIF&)|qP9@$s1;NSHIJH3O{B(9u~ak_N%f+-Q5~taR13<7 zYDCqcs!>Lzdy1;iBO$!x6(i!w$orhBbzzhJ}XNhN*_} zhEay$hCzmYhA=}HLy#fR;BPP)8W?ICelwIalrj`C@CLm>gYJratAD1yufK)vj6A15 zp+BhKjqZ+IuV1N8)z3$FNKVp^)yL_F>Idk1qu)vq_3iX6^-c9geO-NZeMP;SzJ$Jz zp3!S{KXe~;uXK-fcXiiwnYz=uqq_aNow`lBKXl7Np{|y$imtq_w63U5&>7I>pkK*%D8iEK{>kiKMNvL5L{RwB!gCCS1hOOo1uw4cy#&OFiH)85ct z(w@;C(;m>KYd34xYL{ylYv*c{v=g)m+7a3z+6Zm9wyV~xZKG|j_11c7Yip}&-L)>- zV%h>)O6%bG&GEhC3&)3!cO0)OVTZ1|@?ZJPUxEA;$X|i{706!!>nor`H(Qg0BftSb z05mc*KLP&$egM7$z5%`hz5qT0J^?-gJ^1`MAAmK0)qquim4Fo(4w~hFWq_rCG{6!-Dj)^0 zm~^BZG>c?iC~LB;3uK)y>pWTK$~s5Z*|N@(b*8K{WSuT+lC0BYohs`TStrXnN!E$7 zPLMTG*735AlU0&+tgK^XO~5_dMq?R;B_7L2EOA(3v5dg-JC@;ChGB`pG89WRmLXUM zV;O{HAeJaB1F%G5iNMkyOFt}qvGl>x8%r-N;aGZN3BwYKr3aSoSh`{9iX{X~7c8Bz zbiyKH>4?RQB^XN(mJV3jV`+z_EtWP|T4M>s(h5remX=srU}=t}85Vymepr05G{xeB z#T$zk788~xSQ=w7VrhiM6H7xZ4Y1V5QV&aAEOoHd#!?GQO)NFAcwniHr5cv1SgK(8 z4NGM#m9SLAQUQxQmhxE2VJVBH3>G&mu2@{Kl*UpDi!+vzSV~|ij-?ouqF9Pxal%p< zOCc-;u@t}}VBxWFSXe9!78(nM#ehYRMTdpNqQ&Bf#lg{$qEP&=c~1ptgfA%ie$<@k^c$dAnFCM6-E6kg;Zg_Fat&WV}&@>C13!G_PYxr>J`uuMfygeu25a5h@$)w zLLq?>v?#*=$iLzr^LJ5npUI!*kMjFbWWS03gI~rkLQ(xRK9L{I|BfR1{(Mh9gbzm1 zd^6sQZ^+j|k$ic+G+&e#P!vyaU%7YOa}>ef=B{w(xsxb*-@|R=HgKy@4Z_In< z1@n-(gZds^U`{cInZ3+*W+StjNn?_kS923tBV+JyPnNX%P(}8Kl_%ThG`b-U` zGEEm<;y^G#LucKGcDfB#gIz5pdL&wt5bR^x2 z?nZZ{+tMv)AG#4;hpt9fpk3+WbU~V?9jWis2kIsDh{~d_Q5UJds3Yiv#t!OFY7MoN zT1d^Nrc&dnQPgm15Y>+gqq?9cI{zzw1@c!Qe+BYaAb$n&SK$B43gF;FM?e5fKmbfY z08BstOh5okKmbfY08BstOh5okKmbfY08BstOh5okKmbfY08BstOh5okKmbfY08Bst zOh5okKmbfY08BstOh5okKmbfY08BstOh5okKmbfY08BstOh5okKmbfY08BstOh5ok zKmbfY08BstOh5okKmbfY08BstOh5okKmbfY08BstOh5okKmbfY08BstOh5okKmbfY z08BstOh5okKmbfY08BstOh5okKmbfY08BstOh5ok$N?}($N?}($N?}($N?}($N?}( z$N?}($N?}($N?}($N?}($N?}($N?}($N?}($N?}($N?}($N?}($N?}($N?}($N?}( z$N?}($N?}($N?}($N?}($N?}($N?}($N?}($N?}($N?}(EXFTIl2`;-2uKDj0L%x> z1Iz`?0n7%>0?Y)=089rY0j2?_0;T{a1113`0ww?w0pkJV01{vX3)06GIY0YpGYfEf@B2m*8fvm5- z{s2FKFQ6&F2jC6x0+;|z0F41gKqG)Bpdp|Epgy1;pe~>epf;cupeCROzynYnPz_KO zPzCTCpfaEmpdz3Gz#UK?P!3QQPzK-za0R#kN&`v(oB<^PB>=?%#Q;SCMF38K!hk}6 af`9@50l)(|01IFMG=Ksa0D6E9K>j~L5jv^> literal 0 HcmV?d00001