Skip to content
Snippets Groups Projects
postprocess_cascade_case.ipynb 34.7 KiB
Newer Older
many's avatar
many committed
{
 "cells": [
  {
   "cell_type": "markdown",
   "source": [
    "Well use a GenericCascadeCase and we will load the example dataset, generated with OpenFOAM.\n",
    "We don't necessarily need a complete dataset for postprocessing procedures."
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%% md\n"
    }
   }
  },
  {
   "cell_type": "code",
Malte Nyhuis's avatar
Malte Nyhuis committed
   "execution_count": 5,
Malte Nyhuis's avatar
Malte Nyhuis committed
   "outputs": [],
many's avatar
many committed
   "source": [
Malte Nyhuis's avatar
Malte Nyhuis committed
    "from ntrfc.cascade_case.solution import GenericCascadeCase\n",
many's avatar
many committed
    "import pyvista as pv\n",
Malte Nyhuis's avatar
Malte Nyhuis committed
    "import importlib.resources\n",
    "\n",
many's avatar
many committed
    "\n",
    "pv.set_jupyter_backend('trame')\n",
many's avatar
many committed
    "\n",
Malte Nyhuis's avatar
Malte Nyhuis committed
    "inlet_path =  importlib.resources.files(\"ntrfc\") / \"data/openfoam_cascade_case/boundary/inlet.vtp\"\n",
    "outlet_path = importlib.resources.files(\"ntrfc\") / \"data/openfoam_cascade_case/boundary/outlet.vtp\"\n",
    "blade_path = importlib.resources.files(\"ntrfc\") / \"data/openfoam_cascade_case/boundary/blade_wall.vtp\"\n",
    "fluid_path = importlib.resources.files(\"ntrfc\") / \"data/openfoam_cascade_case/internal.vtu\"\n",
many's avatar
many committed
    "\n",
    "case = GenericCascadeCase()\n",
    "\n",
Malte Nyhuis's avatar
Malte Nyhuis committed
    "case.read_meshes(inlet_path,\"inlet\")\n",
    "case.read_meshes(outlet_path,\"outlet\")\n",
    "case.read_meshes(blade_path,\"blade\")\n",
    "case.read_meshes(fluid_path,\"fluid\")\n"
many's avatar
many committed
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
Malte Nyhuis's avatar
Malte Nyhuis committed
   "execution_count": 6,
Malte Nyhuis's avatar
Malte Nyhuis committed
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/many/miniconda3/envs/NTRfC/lib/python3.10/site-packages/pyvista/core/filters/poly_data.py:2848: PyVistaFutureWarning: The default value of the ``capping`` keyword argument will change in a future version to ``True`` to match the behavior of VTK. We recommend passing the keyword explicitly to prevent future surprises.\n",
      "  warnings.warn(\n"
     ]
    }
   ],
many's avatar
many committed
   "source": [
    "\n",
    "case.set_active_blade_slice()\n",
    "case.compute_domainparams_from(0.01)"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
Malte Nyhuis's avatar
Malte Nyhuis committed
   "execution_count": 8,
many's avatar
many committed
   "outputs": [
Malte Nyhuis's avatar
Malte Nyhuis committed
    {
     "data": {
      "text/plain": "0"
     },
Malte Nyhuis's avatar
Malte Nyhuis committed
     "execution_count": 8,
Malte Nyhuis's avatar
Malte Nyhuis committed
     "metadata": {},
     "output_type": "execute_result"
    },
many's avatar
many committed
    {
     "data": {
      "text/plain": "<Figure size 640x480 with 1 Axes>",
Malte Nyhuis's avatar
Malte Nyhuis committed
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAHJCAYAAACmFmJFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABdKElEQVR4nO3deVxUVf8H8M/MyC4DAiGoiIgrivqoYWrmnrj7uJbm1vqo5PZUWv0MzTI1NcvKdq20zKx8XJDSTM0Vd0XUXBBNQcMFVBaHmfv7Y5yRgQHmDnNn7sx83q+Xr5rrmTtnjgN8Oed7vkchCIIAIiIiIhlTOroDRERERBVhwEJERESyx4CFiIiIZI8BCxEREckeAxYiIiKSPQYsREREJHsMWIiIiEj2GLAQERGR7DFgISIiItljwELkImbOnAmFQoHs7OwK29apUwdjxoyx2Wtv27YNCoUC27Zts8n9Lly4AIVCgeXLl9vkfmIsX74cCoUCFy5cMF7r1KkTOnXqZPe+ENEDDFiIiIhI9qo4ugNERHL322+/OboLRG6PAQsRUQU8PT0d3QUit8clISIXk52djaFDh0KtViM4OBiTJk1CQUFBuc+5ceMGXnrpJcTGxqJq1apQq9Xo2bMnjh49Wqrt33//jQEDBsDPzw+hoaGYMmUKCgsLzd533759iI+PR0BAAHx9fdGxY0fs2rXL6ve2detWdOjQAX5+fggMDET//v1x8uRJkzYZGRkYP348GjZsCB8fHwQHB2PIkCEmOSkGJ06cQJcuXeDj44NatWrhrbfegk6nK9WuZA6LIWdn9erVePvtt1GrVi14e3uja9euOHv2bKnnf/TRR6hbty58fHwQFxeHP//8k3kxRCJxhoXIxQwdOhR16tTBO++8g7179+KDDz7AzZs38c0335T5nPPnz2Pt2rUYMmQIoqKicPXqVXz66afo2LEj0tLSUKNGDQBAfn4+unbtiosXL2LixImoUaMGvv32W2zdurXUPbdu3YqePXuiVatWSExMhFKpxLJly9ClSxf8+eefiIuLE/W+tmzZgp49e6Ju3bqYOXMm8vPzsWTJErRv3x6HDh1CnTp1AAD79+/H7t278cQTT6BWrVq4cOECli5dik6dOiEtLQ2+vr4AgKysLHTu3BlFRUWYPn06/Pz88Nlnn8HHx8fiPs2dOxdKpRIvvfQScnJyMH/+fIwYMQL79u0ztlm6dCkSEhLQoUMHTJkyBRcuXMCAAQNQrVo11KpVS9QYELk1gYhcQmJiogBA6Nevn8n18ePHCwCEo0ePGq9FRkYKo0ePNj4uKCgQtFqtyfPS09MFLy8v4c033zReW7x4sQBAWL16tfHa3bt3hXr16gkAhD/++EMQBEHQ6XRC/fr1hR49egg6nc7YNi8vT4iKihK6d+9e7ntJT08XAAjLli0zXmvRooUQGhoqXL9+3Xjt6NGjglKpFEaNGmXyGiXt2bNHACB88803xmuTJ08WAAj79u0zXrt27ZoQEBAgABDS09ON1zt27Ch07NjR+PiPP/4QAAiNGzcWCgsLjdfff/99AYBw/PhxQRAEobCwUAgODhYefvhhQaPRGNstX75cAGByTyIqH5eEiFzMhAkTTB6/+OKLAICkpKQyn+Pl5QWlUv/tQKvV4vr166hatSoaNmyIQ4cOGdslJSUhPDwcgwcPNl7z9fXF888/b3K/I0eO4MyZMxg+fDiuX7+O7OxsZGdn4+7du+jatSt27NhhdumlLJmZmThy5AjGjBmDoKAg4/VmzZqhe/fuJu+t+AyJRqPB9evXUa9ePQQGBpZ6L4888ojJTM9DDz2EESNGWNyvsWPHmuS3dOjQAYB+xgoADhw4gOvXr+O5555DlSoPJrRHjBiBatWqWfw6RMQlISKXU79+fZPH0dHRUCqVZnM4DHQ6Hd5//318/PHHSE9Ph1arNf5dcHCw8f8zMjJQr149KBQKk+c3bNjQ5PGZM2cAAKNHjy7zNXNyciz+oZ2RkWH2dQCgcePG+PXXX3H37l34+fkhPz8f77zzDpYtW4bLly9DEAST1yx+zzZt2pS6n7nXKEvt2rVNHhvez82bN036Xa9ePZN2VapUMS5hEZFlGLAQubiSwYU5c+bMwYwZM/D0009j9uzZCAoKglKpxOTJk0XNhBgYnvPuu++iRYsWZttUrVpV9H0t8eKLL2LZsmWYPHky2rZti4CAACgUCjzxxBNWvZfyqFQqs9eLB0lEZBsMWIhczJkzZxAVFWV8fPbsWeh0unJ/o1+zZg06d+6ML7/80uT6rVu3EBISYnwcGRmJ1NRUCIJgEgidPn3a5HnR0dEAALVajW7dulXm7Rhf19zrAMCpU6cQEhICPz8/43sZPXo0Fi5caGxTUFCAW7dulbqnYSaoOHOvUdl+nz17Fp07dzZeLyoqwoULF9CsWTObvRaRq2MOC5GL+eijj0weL1myBADQs2fPMp+jUqlKzQr8+OOPuHz5ssm1Xr164cqVK1izZo3xWl5eHj777DOTdq1atUJ0dDQWLFiAO3fulHq9f/75x7I3c194eDhatGiBr7/+2iTwSE1NxW+//YZevXqV+16WLFlissxleC979+5FSkqKSb9Wrlwpqm/lad26NYKDg/H555+jqKjIeH3lypXGZSMisgxnWIhcTHp6Ovr164f4+Hjs2bMHK1aswPDhw9G8efMyn9OnTx+8+eabGDt2LNq1a4fjx49j5cqVqFu3rkm75557Dh9++CFGjRqFgwcPIjw8HN9++61xq7CBUqnEF198gZ49e6JJkyYYO3YsatasicuXL+OPP/6AWq3G+vXrRb2vd999Fz179kTbtm3xzDPPGLc1BwQEYObMmSbv5dtvv0VAQABiYmKwZ88ebNmyxSQXBwBeeeUVfPvtt4iPj8ekSZOM25ojIyNx7NgxUX0ri6enJ2bOnIkXX3wRXbp0wdChQ3HhwgUsX74c0dHRFi3XEdF9Dt2jREQ2Y9jWnJaWJgwePFjw9/cXqlWrJiQkJAj5+fkmbc1ta/7vf/8rhIeHCz4+PkL79u2FPXv2lNrOKwiCkJGRIfTr10/w9fUVQkJChEmTJgnJyckm25oNDh8+LAwcOFAIDg4WvLy8hMjISGHo0KHC77//Xu57MbetWRAEYcuWLUL79u0FHx8fQa1WC3379hXS0tJM2ty8eVMYO3asEBISIlStWlXo0aOHcOrUqVLvWRAE4dixY0LHjh0Fb29voWbNmsLs2bOFL7/80uJtzT/++KNF/f7ggw+EyMhIwcvLS4iLixN27doltGrVSoiPjy93HIjoAYUgMDuMiMiedDodHnroIQwcOBCff/65o7tD5BSYw0JEJKGCgoJSOTXffPMNbty4wdL8RCJwhoWISELbtm3DlClTMGTIEAQHB+PQoUP48ssv0bhxYxw8eJAHKxJZiEm3REQSqlOnDiIiIvDBBx/gxo0bCAoKwqhRozB37lwGK0QicIaFiIiIZI85LERERCR7DFiIiIhI9lwih0Wn0+HKlSvw9/dnISYiIiInIQgCbt++jRo1ahhPjC+LSwQsV65cQUREhKO7QURERFa4dOkSatWqVW4blwhY/P39AejfsFqtlvS1NBoNfvvtNzz++OPw8PCQ9LXcFcdYehxj++A4S49jLD0pxzg3NxcRERHGn+PlcYmAxbAMpFar7RKw+Pr6Qq1W84tDIhxj6XGM7YPjLD2OsfTsMcaWpHMw6ZaIiIhkjwELERERyZ4sAhatVosZM2YgKioKPj4+iI6OxuzZs0udv0FERETuSRY5LPPmzcPSpUvx9ddfo0mTJjhw4ADGjh2LgIAATJw40SavIQgCioqKoNVqK3UfjUaDKlWqoKCgoNL3IvPkMMYqlQpVqlThNnkiIpmQRcCye/du9O/fH7179wagP3vj+++/R0pKik3uf+/ePWRmZiIvL6/S9xIEAWFhYbh06RJ/mElELmPs6+uL8PBwnvdCRCQDsghY2rVrh88++wx//fUXGjRogKNHj2Lnzp1YtGiR2faFhYUoLCw0Ps7NzQWg/81co9GYtNXpdEhPT4dKpUJ4eDg8PDwq9UNQEATcvXsXfn5+DFgk4ugxFgQBGo0G//zzD86fP4+oqKgKCxo5G8PXScmvF7ItjrP0OMbSk3KMxdxTFocf6nQ6vPbaa5g/fz5UKhW0Wi3efvttvPrqq2bbz5w5E7NmzSp1/bvvvoOvr6/JtSpVqiAsLAy1atWCl5eXJP0n11RYWIi///4bmZmZXP4jIpJAXl4ehg8fjpycnArLksgiYFm1ahVefvllvPvuu2jSpAmOHDmCyZMnY9GiRRg9enSp9uZmWCIiIpCdnV3qDRcUFODSpUuoU6cOvL29K91XQxlhHgMgHbmMcUFBAS5cuICIiAibfHbkRKPRYPPmzejevTtrV0iI4yw9jrH0pBzj3NxchISEWBSwyGJJ6OWXX8b06dPxxBNPAABiY2ORkZGBd955x2zA4uXlZXa2xMPDo9RgarVaKBQKKJVKm0zr63Q6ADDek2xPLmOsVCqhUCjMfq5chSu/NznhOEuPYyw9KcZYzP1kEbDk5eWV+sGkUqmMP7iIiGxKpwUydgN3rgJ+DwGCANy9Btz9R/+4anVx1/KyH9wnL1t/PbKdo98lkUuRRcDSt29fvP3226hduzaaNGmCw4cPY9GiRXj66acd3TW306lTJ7Ro0QKLFy92yOtfuHABUVFR2LFjB9q3b2+2zbZt29C5c2fcvHkTgYGB9u0gOZ+SwUnGbiDlEyD/lrSv6xsCxD4J4GF9H8Df/okqQxYBy5IlSzBjxgyMHz8e165dQ40aNfDCCy/gjTfecHTXXFZZP/R//vlnh06rRkRE4PLly9xKTNYzBCi3M4Fz24DTG4CCHPv3Iy8b2P850Pxh4IPmQOsxQO22pjMwSpX9+0XkpGQRsPj7+2Px4sUO+63eUlqdgH3nryPj2i1EhmrQpm4IVErXSrwNCgpy6OurVCqEhYUZt6oTWcQQpJxOAo79AORdd3SPTOXfArbPNb3mGwI0Gwo07MXghcgCzBq1UHJqJh6dtxXDv0jBq+v+wvAvUvDovK1ITs2U7DXXrFmD2NhY+Pj4IDg4GN26dcPdu3cB6JduJk+ebNJ+wIABGDNmjPFxYWEhpk2bhoiICHh5eaFevXr48ssvceHCBXTu3BkAUK1aNSgUCuPzSt735s2bGDVqFKpVqwZfX1/07NkTZ86cMf798uXLERgYiF9//RWNGzdG1apVER8fj8zMssfl5s2bGDFiBB566CH4+Pigfv36WLZsGQD9kpBKpcLx48eN7ZOSktCgQQP4+Pigc+fOuHDhQql77ty5Ex06dICPjw8iIiIwceJE41iRC9NpgW3zgHfrAl/3AfZ+LL9gpSx52fr+ft0HWNAASH4VOL9dPyt0fA2Q/uf9pSQiAmQywyJ3yamZGLfiEEru/87KKcC4FYew9KmWiG8abtPXzMzMxJNPPon58+fj3//+N27fvo0///xT1PlKo0aNwp49e/DBBx+gefPmSE9PR3Z2NiIiIvDTTz9h0KBBOH36NNRqNXx8fMzeY8yYMThz5gzWrVsHtVqNadOmoVevXkhLSzMuHeXl5WHBggX49ttvoVQq8dRTT+Gll17CypUrzd5zxowZSEtLw6ZNmxASEoKzZ88iPz/fbNtLly5h4MCBmDBhAp5//nkcOHAA//3vf03anDt3DvHx8Xjrrbfw1Vdf4Z9//kFCQgISEhKMgRC5iJLLPSd+AYoqX8Ha4QzBy96PTa+rawDx84CYfo7pF5GMMGCpgFYnYNb6tFLBCgAIABQAZq1PQ/eYMJsuD2VmZqKoqAgDBw5EZGQkAP12b0v99ddfWL16NTZv3oxu3boBAOrWrWv8e8PST2hoaJmJq4ZAZdeuXWjXTr/jYeXKlYiIiMDatWsxZMgQAPo9+p988gmio6MBAAkJCXjzzTfL7NvFixfxr3/9C61btwagP4qhLEuXLkV0dDQWLlwIAGjYsCGOHz+OefPmGdu88847GDFihHFmqH79+vjggw/QsWNHLF261OVqqLittHVA8jQg94qje2I/uVeA1SOBTq8Bj73EZSNya1wSqkBK+g1k5hSU+fcCgMycAqSk37Dp6zZv3hxdu3ZFbGwshgwZgs8//xw3b960+PlHjhyBSqVCx44dre7DyZMnUaVKFbRp08Z4LTg4GA0bNsTJkyeN13x9fY3BCgCEh4fj2rVrZd533LhxWLVqFVq0aIFXXnkFu3fvLrcPxV8fANq2bWvy+OjRo1i+fDmqVq1q/NOjRw/jsQzk5AzLPqtHulewUty2OcC79fXLRlwqIjfFgKUC126XHaxY085SKpUKmzdvxqZNmxATE4MlS5agYcOGxh/ASqWy1PJQ8TMZylrikULJXUUKhaLcpauePXsiIyMDU6ZMwZUrV9C1a1e89NJLVr/+nTt38MILL+DIkSPGP0ePHsWZM2dMAilyMsXzU7bNcXRvHC//eumcFwYv5EYYsFQg1N+y5QRL24mhUCjQvn17zJo1C4cPH4anpyd++eUXAMBDDz1kktiq1WqRmppqfBwbGwudToft27ebvbdh23B5Z+Q0btwYRUVF2Ldvn/Ha9evXcfr0acTExFTqvT300EMYPXo0VqxYgcWLF+Ozzz4rsw8lT+3eu3evyeOWLVsiLS0N9erVK/WH26OdkCFQeSdCH6hIXS/FGZlL2GXwQi6OOSwViIsKQniAN7JyCszmsSgAhAV4Iy7KttuB9+3bh99//x2PP/44QkNDsW/fPvzzzz9o3LgxAKBLly6YOnUqNm7ciOjoaCxatAi3bt0yPr9OnToYPXo0nn76aWPSbUZGBq5du4ahQ4ciMjISCoUCGzZsQK9eveDj44OqVaua9KF+/fro378/nnvuOXz66afw9/fH9OnTUbNmTfTv39/q9/bGG2+gVatWaNKkCQoLC7Fhwwbj+yrpP//5DxYuXIiXX34Zzz77LA4ePIjly5ebtJk2bRoeeeQRJCQk4Nlnn4Wfnx/S0tKwefNmfPjhh1b3k+xMpwV2LAB2vQ9o7LTDyycIiHteXx/FVpVuz/wKHPoWuHfbPu+heMKuugbQcgwQHM1aL+RyGLBUQKVUILFvDMatOAQFYBK0GFJsE/vG2Lwei1qtxo4dO7B48WLk5uYiMjISCxcuRM+ePQEATz/9NI4ePYpRo0ahSpUqmDJlinGrssHSpUvx2muvYfz48bh+/Tpq166N1157DQBQs2ZNzJo1C9OnT8fYsWMxatSoUoEAACxbtgyTJk1Cnz59cO/ePTz22GNISkqqVHE5T09PvPrqq7hw4QJ8fHzQoUMHrFq1ymzb2rVr46effsKUKVOwZMkSxMXFYc6cOSZVkJs1a4bt27fj9ddfR4cOHSAIAqKjozFs2DCr+0h2ZK9ApXhwImXxtuhOwONv6d/T/q9se++K5F4xXT7jLiNyIbI4rbmycnNzERAQYPa0x4KCAqSnpyMqKqpSu0WSUzMxa32aSQJueIA3EvvG2HxLs7vT6XTIzc2FWq126OGHtvrsyJFGo0FSUhJ69erluMrGUgYqnv7Av54CAiP0Mx/+4Q6ZbdAUFiAp+Vf0qqOBh3+xGRh7F7dz4V1Gsvgsuzgpx7i8n98lcYbFQvFNw9E9Jgz7zmffr3Qb6JKVbonsIm0dsH4SkG/b3XXwrAq0myifH86GPsT0Bwzf6A0zMPaszLttDrDvU6D5MFbWJafFgEUElVKBR+oGIybE4/5v/wxWiERLXQusGW3be/oEAW3+I59ApSJKFRDVQf/HELwUP5yxZBl/WzDsMtr7MY8FIKfEgIWI7EOnBbbPB7bPq7itJVoMB+p2dthyj80YgheDuh2B6k2kLZJXMlGXeS7kBBiwEJG0bJ2r4hME9H3ftX/AxvQDGvW2z7KRoZruI+M540KyxoCFiKRjy1wVueWnSM3cspGUwQtnXEjmGLAQkTRslaviboGKOfYMXnh+EckUAxYisi1b5aowUDGvrITd6+eAg8v0J1nbwrY5+vu1GstCdCQLDFiIyHZssQTEQMVyJRN2H3tJny9kq7OXbmeyEB3JBgMWIrKNtHX6pQSrKYCO04COrzBQsZZSBXSaBoQ2lmaXEZeLXJ5WJyAl/Qau3S5AqL/tj52pDAYs5JSWL1+OyZMnm5yfVNLMmTOxdu1aHDlyxG79cls6LbB+YuXuMWQ50GSALXpDUu8yMiwX9ZzP2RYXUlZF9zd6N3Rgrx7gac3klIYNG4a//vrL0d0gg+3zgfyb1j3XJwgY+i2DFVszLBfFvwO8dAYYvUG/ddk32Db3v52pn23hSdEuITk1E+NWHDIJVgAgK6cAU3444phOlcAZFjF0WuDCLnhcSwdCo4A67WU3JXrv3j14eno6uhsmpOiTj48PfHx8bHpPsoIxwdaKyqzMVbGfsnYZ7f248vdm9Vynp9UJmLU+DeYOFhRKtHPkaU2cYbFU2jpgcVMov+kLv+SJUH7TF1jcVH9dIp06dUJCQgISEhIQEBCAkJAQzJgxA8XPq6xTpw5mz56NUaNGQa1W4/nnnwcA7Ny5Ex06dICPjw8iIiIwceJE3L37oGjXxx9/jPr168Pb2xvVq1fH4MGDjX+3Zs0axMbGwsfHB8HBwejWrZvxuZ06dcLkyZNN+jlgwACMGTOm0n0q6ejRo+jcuTP8/f2hVqvRqlUrHDhwAIB+SSgwMNCk/dy5c1G9enX4+/vjmWeeQUFBQal7fvHFF2jcuDG8vb3RqFEjfPyxDb5hu6u0dcC79awIVhRAx+nA9Iv6fAv+YLOv4jMvQ7/VJ9LagqF67td9JP/eSLaVkn6j1MxKcYafOAczrJxFtREGLJZIWwesHlU6gS03U39dwi/Mr7/+GlWqVEFKSgref/99LFq0CF988YVJmwULFqB58+Y4fPgwZsyYgXPnziE+Ph6DBg3CsWPH8MMPP2Dnzp1ISEgAABw4cAATJ07Em2++idOnTyM5ORmPPfYYACAzMxNPPvkknn76aZw8eRLbtm3DwIEDIfZQb7F9MmfkyJGoVasW9u/fj4MHD2L69OllnhS6evVqzJw5E3PmzMGBAwcQHh5eKhhZuXIl3njjDbz99ts4efIk5syZgxkzZuDrr78W9d4I+horq0datxtoyHKg86sMVOQgph8wOVWfRGtLhuRcLhc5hWu3yw5Wisu+UyhxT8rHJaGK6LT6bPsyJ8sUQPJ0fYKbBN+AIyIi8N5770GhUKBhw4Y4fvw43nvvPTz33HPGNl26dMF///tf4+Nnn30WI0aMMM6E1K9fHx988AE6duyIpUuX4uLFi/Dz80OfPn3g7++PyMhI/Otf/wKgD1iKioowcOBAREZGAgBiY2NF91tsn7y9vUvd4+LFi3j55ZfRqFEj43PKsnjxYjzzzDN45plnAABvvfUWtmzZYjLLkpiYiIULF2LgwIEAgKioKKSlpeHTTz/F6NE2PozPlaX+DKx5WvzzPP2BAR8zSVNupNxZxOq5TiHUv/T3X3MyrudJ3JPycYalIhm7K/gCFoDcy/p2EnjkkUegUDw4Fbpt27Y4c+YMtNoHv7G0bt3a5DlHjx7F8uXLUbVqVeOfHj16QKfTIT09Hd27d0dkZCTq1q2LkSNHYuXKlcjL038Qmzdvjq5duyI2NhZDhgzB559/jps3xU8Diu2TOVOmTMGzzz6Lbt26Ye7cuTh37lyZr3fy5Em0adPG5Frbtm2N/3/37l2cO3cOzzzzjEkf3nrrrXLvSyX8NgNYMxbmA/hyeKmBV87zB5acGWZbRm8ABn2pn3XxD7fNvXOZoCtncVFBCA/whqKCdh9vO4vkVBsVJrQCA5aK3Llq23YS8PPzM3l8584dvPDCCzhy5Ijxz9GjR3HmzBlER0fD398fhw4dwvfff4/w8HC88cYbaN68OW7dugWVSoXNmzdj06ZNiImJwZIlS9CwYUNjUKFUKkstD2k0mkr3yZzExEScOHECvXv3xtatWxETE4NffvnFqjG6c+cOAODzzz836UNqair27t1r1T3dTurPwO4PrHtu/4+AKvJKBiczDPktsYP1sy5TTthouej+9wzmuMiSSqlAYt8Yi34NmbU+DVqdyF9YbIQBS0WqVrdtO5H27dtn8njv3r2oX78+VKqyl59atmyJtLQ01KtXr9Qfw26dKlWqoFu3bpg/fz6OHTuGCxcuYOvWrQAAhUKB9u3bY9asWTh8+DA8PT2NgcJDDz2EzMwHEbZWq0VqamqF78OSPpnToEEDTJkyBb/99hsGDhyIZcuWmW3XuHFjs2NlUL16ddSoUQPnz58v9fpRUVEV9t/tnVgL/PSM+OcpVMDgrzmz4qwMy0W2TM4FOOMiQ/FNwzGlW9nL7oA+7MzMKUBKug0OM7UCc1gqEtlO/4Wamwnz0+AK/d9HtpPk5S9evIipU6fihRdewKFDh7BkyRIsXLiw3OdMmzYNjzzyCBISEvDss8/Cz88PaWlp2Lx5Mz788ENs2LAB58+fx2OPPYZq1aohKSkJOp0ODRs2xL59+/D777/j8ccfR2hoKPbt24d//vkHjRs3BqDPTZk6dSo2btyI6OhoLFq0qNzibZb2qaT8/Hy8/vrrGDJkCKKiovD3339j//79GDRokNn7T5o0CWPGjEHr1q3Rvn17rFy5EidOnEDdunWNbWbNmoWJEyciICAA8fHxKCwsxIEDB3Dz5k1MnTq1wvfgtipziOHgr1hfxRWULERX6e3QxWZcmOMiG3VC/CpuBMuTdG2NAUtFlCr9F9LqUQAUMA1a7q/4xc+VbMfDqFGjkJ+fj7i4OKhUKkyaNMm4TbgszZo1w/bt2/H666+jQ4cOEAQB0dHRGDZsGAAgMDAQP//8M2bOnImCggLUr18f33//PZo0aYKTJ09ix44dWLx4MXJzcxEZGYmFCxeiZ8+eAICnn34aR48exahRo1ClShVMmTIFnTt3rvB9VNSnklQqFa5fv45Ro0bh6tWrCAkJwcCBAzFr1iyz7YcNG4Zz587hlVdeQUFBAQYNGoRx48bh119/NbZ59tln4evri3fffRcvv/wy/Pz8EBsbW2qbNt1XmUMMfYKAvu/zB5ArKV7LpXZb2ybosuS/LFiafGtpO1tTCGL3q8pQbm4uAgICkJOTA7VabfJ3BQUFSE9PR1RUlNmdKBZLW1f6C1RdUx+sSPRNuVOnTmjRogUWL14syf3lSqfTITc3F2q1Gkql41YtbfbZkSGNRoOkpCT06tXL/Fbxyhxi2HE6zwO6r8JxdmY6rTSl//3DRZX8d+kxtjOtTsCj87YiK6fA5FdzL5WA+XFaTEtRoVpVH+yc1gUqZUUpupYp7+d3SZxhsdT9KVHdhV3Iv5YOn9AoKGVY6Zao0ipziOHgZUDTgbbtD8mTVNVzDSX/HxnPqrl2Zki+HbfiUKn1BIPEvjE2C1bEYtKtGEoVUOdRaBr1B+o8yi8icj2VOcSw3UQGK+5Kiuq53FHkEPFNw7H0qZYICyg9q/zesBaIb2qjre5W4AyLjG3bts3RXSB3Y80hhgqlvm4HgxUCbJ+gm8sZF3uLbxqO7jFhSEm/gWu3CxDiWwXZJ/eiW2NpdsNaigELEemdWGtdgu3gZdwJRKZsmqDLHUWOoFIq0DZaf7K3RqNB0kkHdwhutCTkArnFZGdu9ZlJXQv8OBqiKtgaaqwwWKHyFK+g+8j4yt+PNVzclssHLIascUPpeSJLGT4zLr/z4MRa4Kex4p83+Cug6QBb94ZckU1zXEpUzf34EZt0keTP5ZeEVCoVAgMDce3aNQCAr6+vydk8Yul0Oty7dw8FBQUO3XLryhw9xoIgIC8vD9euXUNgYGC5VYWd3smNwM8ii8LxEEOqDEOOy44FwLY5lb/f7Sz9f7fMAhr1YI6LC3P5gAUAwsLCAMAYtFSGIAjIz8+Hj49PpQIfKptcxjgwMND42XFZ/xsnrr2XGnj5HM8Fosqx6QnR92dc9n8O7FvCHBcX5hYBi0KhQHh4OEJDQ80e1CeGRqPBjh078Nhjj7n+UoGDyGGMPTw8XHtm5VSS/r+CTtzzeIgh2ZLNS/7jfo7LKKDTq0BwtP6cN866uAS3CFgMVCpVpX8IqVQqFBUVwdvbmwGLRDjGEtNpgeRXgAbzLX+OQgUM+oq/tZLt2bzk//0Zl+LLTZx1cQlMwiByN9vnA/m3xD2HCbZkD7beUWTAnUUugQELkTsRW2uFW5fJ3qSomltyZxGr5zolBixE7sKaWiucWSFH4owLFcOAhcgdiK61ouDMCskDZ1zoPgYsRK4ubZ1+ZkXMjqCO0zizQvIj6YzLKAYtMudWu4SI3I5OC2x6RdxzfIKAjiKfQ2QvJXcV/fqGDW56f8Zl/SR9raGoDtwGLUOcYSFyZTsWALczxT2n7/v8Zk3OIaYfMH6v/v/7fwR0eu3+X1hZcDL/BvBtfy4RyRRnWIhcVdo6caXPWWuFnJEhuI7pD3h42KZ6LovPyRIDFiJXpNMC6yeKe87gr5hkS87PJtVzWXxOjrgkROSKts8H8m9a1pa1VsjVSLGziFuhHY4BC5GrEVscjrVWyJXZbGcRt0I7GpeEiFxJ6lpgzWjL23d4iTMr5Ppsfl4RmOfiAAxYiFyF6OJwANq9KElXiGTLkONyYSfw4yjx52oZMc/F3rgkROQKrCkOB/C3QXJPShVQtyPQdwn0W6Ct3AZdEgvQSYoBC5Gzs6o4XDVp+kLkTGL6AUO/AdThNrqhoP+zfhJwbhsTc22MS0JEzs6a4nDx84Dz0nSHyKkU3wZ95ypw/dz9ZR4FRB0UWpyhAB2XiGyKAQuRM7O2OFzDXsD5JOn6ReRMDEm5BrYoPgcwMdfGGLAQOavKFIfTaCTpEpFLsEnxOYCJubbFgIXIWYktDjeIlWyJLCbFVmiAsy6VwICFyBlZUxyOwQqRdWya5+I8sy5anYCU9Bu4lnPX+NjDgf1hwELkbAxbmC3V6TUGK0SVJVWeC/Bg1mXoN7IJWpJTMzFrfRoycwrgo9Jibhwwd+Fc9HqkGeI69XXIjBC3NRM5E7FbmH2qAY+9JF1/iNyVoeT/qHWAT2Alb3Z/1mXDFODYaoefVZScmolxKw4hM6cAPZQp+M1T/z3n9aKPELdjNPLfjXFIrRkGLETOROwW5jbjuDZOJBWbFqATgLxs4OfnHHpWkVYnYNb6NAgAeihTsNRjMUJhmivnlZ8FwQEF8hiwEDmL1LXitjD7BHF2hcgebF6ADg+WibbNA46vsdusS0r6DWTmFEAJHRI9vgEAKEvEYcbAIXm6XWeCmMNC5AysOSeo7/ucXSGyF5sXoHNMcu612wUAgDjlKdRQ3AAAmAtJFBCA3Mv691s8t0dCsplhuXz5Mp566ikEBwfDx8cHsbGxOHDggKO7ReR4qWvFnROkUAGDv5ZN8h6R2zAk5sYOBjpNA4Z+K82si4RLMaH+3vr/4pZlT7hzVbK+lCSLGZabN2+iffv26Ny5MzZt2oSHHnoIZ86cQbVqPO+E3Jw1MyvcwkwkD5LMuij0yblFBYB/uM1ruMRFBSE8wBvXbgda9oSq1W322hWRRcAyb948REREYNmyZcZrUVFRDuwRkQyI3b4McAszkdzYfDt0seRcwObLRCqlAol9YzBhRR5uCFVRDXfKbqyuqQ+Y7EQWAcu6devQo0cPDBkyBNu3b0fNmjUxfvx4PPfcc2bbFxYWorCw0Pg4NzcXAKDRaKCRuOS44f5Sv4474xhDn8j26xuA0tvy5/jXANpOsqjsPsfYPjjO0nO6Ma7fE4h+HLiUAtzJArYkAnk3rL/fnZvAmueBDqeBoCjALxSIiKvUrEvXhiFY3ekGqqYUQSN4o+j+9yFNye9HTYYAWp3+j5XE/LspBEGw8jhK2/H21g/C1KlTMWTIEOzfvx+TJk3CJ598gtGjS/+GOXPmTMyaNavU9e+++w6+vr6S95eIiIgqLy8vD8OHD0dOTg7UanW5bWURsHh6eqJ169bYvXu38drEiROxf/9+7Nmzp1R7czMsERERyM7OrvANV5ZGo8HmzZvRvXt3eHg4skix6+IYA9icCBz40rK2ChXQ/2OgcW+Lb88xtg+Os/RcYoxPJelnWsTUWCrX/X3IHf5r3axLxh7guyHGhxqlNzbHfoDuxyfCQ1dg2nb4j0BkW6t7mpubi5CQEIsCFlksCYWHhyMmJsbkWuPGjfHTTz+Zbe/l5QUvL69S1z08POz2gbXna7krtx3jE2uB/Ust3xU05Gur81bcdoztjOMsPace49j+QJM+NkzOvW/7Ww/+X0Suiy7vKpQlAxMAHrqC0gFL/jWgEuMu5t9MFgFL+/btcfr0aZNrf/31FyIjIx3UIyIHSV0LrLE00VYJDF7GJFsiVyDlWUWAqPOKTt72RRNL7+tuu4SmTJmCdu3aYc6cORg6dChSUlLw2Wef4bPPPnN014jsR+wW5jbPA00HSNUbInKk4luib2fqq8rmXa/EDS3fEn3WNxZBQhDCcAOKMk4bEARAUCigjGhTiT6JI4vCcQ8//DB++eUXfP/992jatClmz56NxYsXY8SIEY7uGpF9GLYwW7oMBACN+kjXHyJyPMOsS7OhQJ/FsNd5RaFqP6RoG5QZrACAQgEoIQCX9lWiP+LIYoYFAPr06YM+ffgNmNyQ2BOYAbvXPyAiBzOcVyTFMlGnV4HgaP3yTmQ7tKrpB63qoGX3sFmicMVkE7AQuS2xJzADQPxcnhNE5G7scV6RTxBU9/LgqbCwPsqda1a8rnUYsBA5ktgTmBUqYNBXPCeIyF1JnZybf0Ncrkhlit6JxICFyFF4ThARVZbNZ11EZsko7ZcKy4CFyBHEnhNkmFlhsEJEJUk961KeyEelf437GLAQ2ZtOC6yfKO45nFkhIkvZfEt0GTz9TQMliTFgIbK3HQuA/JuWt+cJzEQkVvFZlyre+t1AACpVObfYHRQA0P9Duyb/y6IOC5Hb0GmBvR9b3l5dE3jsJen6Q0Suz7AlWh1uk9spAKDdRLv/IsUZFiJ72rEAKLhleXtuXyYiW7BVcq6nP9DvQ4dU2WbAQmQvorYwK4DBy7l9mYhsx5LkXJ8gAILpsrWXv/6/T6wC6j3msF+iGLAQ2YPYLcwdp/GcICKSVolZF61fKFK0jXDtTgHq5R1HY/88KP3DgBoPA8m/AlGPOnTGlwELkdREncAM/W84HUWW6icissb9WZfk1EzM+iENmTn7jX8VHhCExL710FUmy9JMuiWSkjXF4fq+z7wVIrKb5NRMjFtxCJk5BSbXs3IKMG7FIWw5edVBPTPFgIVIKqlrxZ/A3Ok15q0Qkd1odQJmrU8zm3ZruDZ30yl7dqlMDFiIpGDNzAq3MBORnaWk3yg1s1KcACArt+y/tyfmsBDZmtiy+wbcwkxEdnbttjyCEUtwhoXIlqwpu69QAYO/5lIQEdldqL+3o7tgMQYsRLa0fb64svuA/pwgbmEmIgeIiwpCeIB3mSc0KwCEqeUR1DBgIbKVE2uB7fMsb2+YWeE5QUTkICqlAol9YwCgVNBieDy9ZyO79qksDFiIbMGwI0hMiWvOrBCRDMQ3DcfSp1oiLMB0JiUswBtLn2qJbo2rO6hnpph0S1RZoncE3S+7z5kVIpKJ+Kbh6B4ThpT0G7h2uwCh/t6IiwqCSqmARqNxdPcAMGAhqhyxVWwBlt0nIllSKRVoGx3s6G6UiUtCRNayptYKy+4TEVmFMyxE1rBmZgVg2X0iIitxhoVILGtmVlhrhYioUjjDQiSGtVVsB3/FJFsiokrgDAuRpSpTxZbBChFRpXCGhchS1laxZbBCRFRpnGEhsgSr2BIRORRnWIgqYs2OIM6sEBHZFGdYiMpjVRVbzqwQEdkaZ1iIysIqtkREssEZFiJzWMWWiEhWOMNCVJK1tVZYxZaISDKcYSEqrjK1VljFlohIMpxhISqOtVaIiGSJMyxEBqy1QkQkW5xhIQJYa4WISOY4w0LEWitERLLHGRZyb6y1QkTkFDjDQu6LtVaIiJwGZ1jIPZ1Yy1orREROhDMs5H6O/wysETmzwlorREQOxRkWci+/zQB2fyD+edwRRETkUJxhIfeR+rP4YIW1VoiIZIEzLOQeTqwFfnpG/PM4s0JEJAsMWMj1WbN1WaECBjFYISKSCwYs5Lp0Wv3ZQGLK7RtwZoWISFYYsJBrSlsHrJ8E5N8Q9zyFEhi8jMEKEZHMMGAh15O2Dlg9CoAg/rkDv2SwQkQkQ9wlRK6l6B6wYTKsClbaTQRiB9q6R0REZAOcYSHXkbYO2DAFyLsu7nkKJTDoS6ApgxUiIrliwEKuwZqdQAbMWSEikj0GLOTcKrMTiFuXiYicBgMWcl6pa4F1CcC929Y9n1uXiYicBgMWck7WngkEAD5B+lOXeZAhEZHTYMBCzseaM4EAwNMfaPci8NhLgFJl+34REZFkGLCQ86hMvopvCDD1JFDF0/b9IiIiyTFgIedQ2XyVPu8xWCEiMkOrE5CSfgPXbhcg1N8bcVFBUCkVju5WKQxYSP4qk69i2AnEfBUiolKSUzMxa30aMnMKjNfCA7yR2DcG8U3DHdiz0ljpluTN2nwVg8FfAU0H2Kw7RESuIjk1E+NWHDIJVgAgK6cA41YcQnJqpoN6Zh4DFpKvPxcBa5627rne1YCh33LbMhGRGVqdgFnr08weYmK4Nmt9GrQ6K445kQgDFpKfkxv1/925CFadCdRkIPDKOS4DERGVISX9RqmZleIEAJk5BUhJF3nivYQYsJB86LTAj2OBtS9Y93yFUl9mf8gyblsmIirHtdtlByvWtLMHJt2SPKStA9ZPBPJvAkpv6+7BM4GIiCwS6m/Z91lL29kDZ1jI8VLXAqtH6oMVazBfhYhIlLioIIQHeKOszcsK6HcLxUUF2bNb5ZJlwDJ37lwoFApMnjzZ0V0hqaX+DKwZY/3zma9CRCSaSqlAYt8YACgVtBgeJ/aNkVU9FtkFLPv378enn36KZs2aOborJCVDvsqasbAqsZb5KkRElRLfNBxLn2qJsADTZZ+wAG8sfaql7OqwyCqH5c6dOxgxYgQ+//xzvPXWW47uDklBpwV2LAB2vQ9o7lp/H+arEBFVWnzTcHSPCWOlW7EmTJiA3r17o1u3buUGLIWFhSgsLDQ+zs3NBQBoNBpoNBpJ+2i4v9Sv45JObgSSXnpQXr+M5FrN/esac3/vHQj0nA806AXw38Bq/BzbB8dZehxj22hdWw1ADQDQaYug0z74OynHWMw9FYIgyKIqzKpVq/D2229j//798Pb2RqdOndCiRQssXry4VNuZM2di1qxZpa5/99138PX1tUNviYiIqLLy8vIwfPhw5OTkQK1Wl9tWFgHLpUuX0Lp1a2zevNmYu1JewGJuhiUiIgLZ2dkVvuHK0mg02Lx5M7p37w4PDw9JX8slnEoCkl8B8m9Z/BSN0hubYz9A9+MT4aErADz9gd4LgEa9peunm+Hn2D44ztLjGEtPyjHOzc1FSEiIRQGLLJaEDh48iGvXrqFly5bGa1qtFjt27MCHH36IwsJCqFQPEiu9vLzg5eVV6j4eHh52+8Da87WcVupa4KfRVj/dQ1cAjw6TgY6vMLFWIvwc2wfHWXocY+lJMcZi7ieLgKVr1644fvy4ybWxY8eiUaNGmDZtmkmwQk5ApwW2zwe2z6vcfQYsBZoNtE2fiIjIqckiYPH390fTpk1Nrvn5+SE4OLjUdZK51LXAuoQHibWV0bhv5e9BREQuQRYBC7kAnRb46VngxM+Vu4+nP9D3Q+C8bbpFRESuQbYBy7Zt2xzdBbJU8XOArOXhB7SfBDz2EqDVAeeTbNc/IiJyerINWMhJpK4F1lifWAtAX15/0BcPEmu1ukp3i4iIyqfVCU5RMM6AAQtZxxaJtZ7+QP8PWbGWiMjOklMzMWt9GjJzCozXwgO8kdg3RnYl+Q0YsJA4NimtrwA6TuN2ZSIiB0hOzcS4FYdKneKWlVOAcSsOyfIcIYABC1nKVmcAAcCQ5ZxVISJyAK1OwKz1aWaPnBWgP6l51vo0dI8Jk93yEAMWqljaOmD9JCD/RuXu410N6PcBENPPNv0iIiJRUtJvmCwDlSQAyMwpQEr6DbSNDrZfxyzAgIXKZ4ukWqB0Yi0REdndtdtlByvWtLMnBixknq2q1TKxlohINkL9vW3azp4YsJApm+WqMLGWiEhu4qKCEB7gjaycArN5LAoAYQH6Lc5yo3R0B0hG0tYB79YDts2xTWJt51cZrBARyYhKqUBi3xgA+uCkOMPjxL4xsku4BRiwkEHqWmD1yMon1voEAUO/5RIQEZFMxTcNx9KnWiIswHTZJyzAW7ZbmgEuCbk3nRbI2A2c2gDs+7Ry9/KsCrSbqC+tz1kVIiJZi28aju4xYax0S04gbR2QPA3IvVK5+zBQISJySiqlQnZbl8tjdcBy6NAhLFiwALdu3UJsbCz+85//ICoqypZ9I6nYZKsyk2qJiMh+rM5hGTZsGPr06YO3334bDRo0wMCBA/Hbb7/Zsm9kazot8Mc7wJoxlb8Xk2qJiMiOrJ5hqVatGoYPHw4A+Ne//oUBAwagS5cuePzxx23WObIhW1Wr9QkC+r7ParVERGRXVs+wREdHY+HChRAE/U7uwMBAW/WJbM0WO4A8qwKdXgNePstghYiI7M7qGZbCwkIsXboUixcvRpMmTZCRkYFevXrh8uXLqFmzpi37SNaw2Q4g5qoQEZHjWRyw7N69G2q1Gk2bNgUA/PzzzwCAu3fv4vjx4zh27BiOHTuGJ598EleuXMHZs2el6TFVzFY7gACerExERLJgccAyYcIEJCQkGAMWg6ysLDRp0gSPPPKIzTtHVrDVYYXqmkD8XC7/EBGRLFgcsJw+fRqdOnUqdX3Lli1Yv349NmzYYMt+kTVOrAV+Gmv98xVKIO55oFEfILIdl4CIiEg2LE66VavVuHnzZqnrHTp0wN69e23aKbJC6lrgx9GAoLP+HoOXAT3nAVEdGKwQEZGsWBywxMfHY8GCBaVvoFTi3r17Nu0UiVTZmRWe/0NERDJnccAye/ZsbN++HYMGDcLx48cBAAUFBZg3bx6aNWsmWQepHDotsG1eJWZWFEDH6dyqTEREsmdxDktERAT27t2LcePGoXnz5vDy8kJRURECAgKwfv16KftI5thiJxB3ABERkZMQVYclMjISSUlJuHjxIo4cOQIPDw+0adMGQUFBUvWPzKnsTiDuACIiIidjVeG42rVro3bt2rbuC1nC2nwV7gAiIiInZnWlW3KAysysDF7G5R8iInJaVp8lRHZm9cyKChj8NYMVIiJyapxhkTudFtixANg2x7rnD/6KwQoRETk9BixyVpmdQP419EXgmFhLROT2tDoBKek3cO12AUL9vREXFQSVUuHobonCgEWuKpOv0uk14LGXmFhLRERITs3ErPVpyMwpMF4LD/BGYt8YxDcNd2DPxGEOixxVNl+l0zQGK0REhOTUTIxbccgkWAGArJwCjFtxCMmpmQ7qmXgMWOQmbZ31lWsHfwU0HWDzLhERkfPR6gTMWp8GwczfGa7NWp8Grc5cC/lhwCInOi2w6RXxz+NOICIiKiEl/UapmZXiBACZOQVISb9hv05VAnNY5GTHAuC2FdNz3AlEREQlXLtddrBiTTtH4wyLXKSuFb912b8GT1kmIiKzQv29bdrO0TjDIgfWJNlyJxAREZUjLioI4QHeyMopMJvHogAQFqDf4uwMOMPiaGKTbLkTiIiILKBSKpDYNwaAPjgpzvA4sW+M09RjYcDiSEX3gA2TxT2HO4GIiMhC8U3DsfSplggLMF32CQvwxtKnWjpVHRYuCTlK2jpgwxQg77rlz+n0GvNViIhIlPim4egeE8ZKt2SFtHXA6lGA2VXFMqhr6nNWiIiIRFIpFWgbHezoblQKl4TszVhrRWShnvi5zFkhIiK3xRkWexNba0WhAgZ9xUMMiYjIIq5w0KE5DFjsKW2d+ForLApHREQWcpWDDs3hkpC9iC277xvConBERGQxVzro0BwGLPYiZinINwSYepLLQEREZBFXO+jQHAYs9iB2KajPe0AVT+n6Q0RELsXVDjo0hwGL1MQuBXV6jTMrREQkiqsddGgOAxapiVkKYq0VIiKygqsddGgOAxYpiV0KYq0VIiKyguGgw7I2Lyug3y3kLAcdmsOARSo6LZA8zfL2XAoiIiIrudpBh+YwYJFKxm4g94plbbkUREREleRKBx2aw8JxUjmdZHlbLgUREZGVSla23f5yZxzMuMlKt2SBtHXA3o8ta8ulICIislJ5lW37t6jpwJ7ZHpeEbE3MNmYuBRERkZVcvbJtSQxYbE3MNmYuBRERkRXcobJtSQxYbEnMNuZHxnMpiIiIrOIOlW1LYsBiK2K3MTfsJV1fiIjIpblDZduSGLDYithtzJHtpO0PERG5LHeobFsSAxZb4TZmIiKyE3eobFsSAxZb4DZmIiKyI3eobFsSA5bKEpO7wm3MRERkI65e2bYkFo6rLDG5K1wKIiKiSihZ1bZ7TBi6x4SZXHOVyrYlMWCpLEtzV7iNmYiIKqG8qrauNptiDpeEKkNM7gq3MRMRkZXcraqtObIIWN555x08/PDD8Pf3R2hoKAYMGIDTp087ulvlE5u7wm3MRERkBXesamuOLAKW7du3Y8KECdi7dy82b94MjUaDxx9/HHfv3nV018p2KYW5K0REJDl3rGprjixyWJKTk00eL1++HKGhoTh48CAee+wxB/WqAmd+s6wdc1eIiKgS3LGqrTmyCFhKysnJAQAEBZkveFNYWIjCwkLj49zcXACARqOBRqORtG+G+2sOfgsoLaggWK8nIHGfXI1xjDlukuEY2wfHWXruMMYhvlXgpap4uSfEt4ok4yDlGIu5p0IQBFkteul0OvTr1w+3bt3Czp07zbaZOXMmZs2aVer6d999B19fX6m7SERERDaQl5eH4cOHIycnB2q1uty2sgtYxo0bh02bNmHnzp2oVauW2TbmZlgiIiKQnZ1d4RuuLM353dh86ha6H58ID10F02///gxoxN1BYmk0GmzevBndu3eHh4eHo7vjkjjG9sFxlp6rjvGWk1cxd9MpZOWW/3PGUG3lvWEt0K1xdUn6IuUY5+bmIiQkxKKARVZLQgkJCdiwYQN27NhRZrACAF5eXvDy8ip13cPDQ/oPbME/ADzgoSsoP2B5ZDwQ21/avrg4u/x7ujmOsX1wnKXnSmOcnJqJ8d8dvb8DqPwCcPaswyLFGIu5nywCFkEQ8OKLL+KXX37Btm3bEBUV5egule1GOoAGFbdj3RUiIhKpvC3MBkF+HpjRpwnC1K5b1dYcWWxrnjBhAlasWIHvvvsO/v7+yMrKQlZWFvLz8x3dNVM6LXBkZcXt/Guw7goREYlW0RZmALhxV4MwtTfaRge7TbACyCRgWbp0KXJyctCpUyeEh4cb//zwww+O7pqpjN3AbQuqCbYaw7orREQkGrcwl002S0JOwdJzg4Kjpe0HERG5pFB/C8pliGjnSmQxw+IUxJwbVFWaTG0iInJtcVFBCA/wLjPVVgF9om1clPk6Za6MAYsleG4QERFJSKsTsOfcdWw4dgVPPFwbQOn9QYbHiX1j3Cp3xUAWS0Kyl7Gb5wYREZEkklMzMWt9mkmybaCvfrvvrbwHlWDD7LiFWY4YsFjizlXL2vHcICIiEiE5NRPjVhwqtY05J08DAcCUbvVRJ8QPof7utYXZHAYslrA0J4W1V4iIyELl1VwRoF8CWrX/EnZO6+LWgYoBc1gskXe94jbMXSEiIhEqqrkiAMjMKUBK+g37dUrGGLBURKcFfn214naPz2HuChERWYw1V8RhwFIRSxNu/YKl7wsREbkErU5A9u3CihvCPWuumMMclopYmnBraTsiInJr5nYFmaOAfmeQO9ZcMYcBS0UsTbhlsTgiIqpAWbuCSnL3mivmMGCpCBNuiYjIBiw5idnA3WuumMOApTxMuCUiIhux5CRmAJjRuzHGtI/izEoJDFjKw4RbIiKqJK1OQEr6DWxKzbSofYi/F4MVMxiwlIcJt0REVAmWJtgWx11B5jFgKc/1c5a1Y8ItERGVYGmCrQF3BZWPdVjKotMC+5ZW3M6/BhNuiYjIhJgEW4C7gizBGZay7FgA5N+suF2rMUy4JSIiE5Ym2BpwV1DFGLCYY+nsCgAER0vbFyIichpiE2xHtY1Ez6bhbn8SsyUYsJiTsduy2RWA+StERATAugTbnk3D0TaaO00twYDFHEt3/fhUY/4KERExwdYOmHRrjqWzJm3GMX+FiMjNMcHWPhiwmBPRBlBUNDQK4NEpdukOERHJlzUJtkufaskEW5G4JGTOpX2AoKugkaBvF9XBLl0iIiJ5YYKtfTFgMYcVbomIqBxMsLU/BizmsMItERGVgQm2jsEclpJ0WuDQ8orbVQ3nDiEiIjfDBFvH4QxLSZae0PyvEdwhRETkRrQ6Act3pbOCrYMwYCnJ0ryUoChp+0FERLIhNmeFCba2x4ClJEvzUvxCpe0HERHJgticFYAJtlJgDktJFtVgAVCzlfR9ISIih7ImZyWcCbaS4AxLSRbVYAFw+SBQ7zHp+0NERHZnqLGy6+w/Fi8DMcFWWgxYSrI0h+XuNWn7QUREDmFNjRWACbZSY8BSEnNYiIjcljX5KgAwo3djjGkfxZkVCTFgKSmyHaCuAeRmAmY/svc/jBFx9uwVERFJxLD8k5WTj9kbT4oKVgxF4RisSI8BS0lKFRA/D1g9sowGwoN2RETk1Kxd/gGYs2JvDFiIiMgtWbv8Y8CcFftiwFKSTgskTyungeJBO3jYo0dERGRjYrcrF5fQuR7a1wthUTg7Y8BSUoWl+e9/vC+lcFszEZGTsWa7soEhX2VK9wYMVByAAUtJ3NZMROSSmK/i3BiwlMRtzURELof5Ks6PAUtJ3NZMROQytDoBe89dx/SfjosOVoL8PDCjTxOEqb2ZryIDDFhKsnJbs2Fd9NrtAoT688NNRORo1i4BGb5zz/l3LGdUZIQBiw2Y+6II5/QhEZFdFf/F8UJ2HhZv+cuqJSAu/8gTA5aSRG5rTjqWifHfHSrVKjOnAP9ZcQifPNWSH3oiIoltOXkVb248bVVCrQG3K8sbA5aSRGxrTsqrh4TvD5d7u+k/H0f3mDB++ImIJDTlhyMo0Fr3fZbblZ2D0tEdkB0LtzWnnv4L4787DF0F84238jT4cOtZG3SMiIgMtDoBe85dx4aj+l8wrd39w+3KzoMzLCVZuK35q2N5Ft/yq13nkdClHr8YiIhsoHjeoJdKwPxKbNpkvorzYMBSkoXbmrfcqWPxLXPyi/Dh1rOY1K2+TbpIRORObJVMW1ygjwc+GtESj9QN5i+TToIBS0nGbc2joA9Oin9ZPPhQ60Supr235S80DKvKKJ6ISITKVKc1x/BdfO6gWLSvF2KTe5J9MIfFnJh+wNBvAHWJ4EJdA/j3p1bfdua6E9BWlPRCROTGDLkp/ztyGe9vOYNxKw7ZLFgB9EtAS7l70ylxhqUsMf2ARr31u4buXNXntkS2A7Q64HwSwtTeuHizUNS0ZFZuIZeGiIjKIMVsigBgSrf6qBPix6KeTo4BS3mUKiCqg+k1rQ4AML1nI4z/7mipRaOKcGmIiEhPityU4phQ61oYsFipW+Pq+Gh4S/zf/1Jx4+49Uc+due4Ea7MQkVuz9WyKAc//cV0MWKy05eRVzN542iRYUSgAwYJfD7g0RETuRurZFMX9Pzz/x3UxYLGSuaqKlgQrBu9t+QuAgIQu9fkbABG5lJKHwd68ew+zN9p+NqW46mpvvNq7CYMVF8aARSTDLh9b/Gbw3pYz+D7lEmb24xorEbkGqZZ6ijNJpg3yBi4dxq+TH4O3l6dkr0mOx4BFpIMZN216v6xc/SGJz7Svg24xYVxzJSKnIvVSjznFk2k1Gg2SLh3m9003wIBFpOw7hZLc98tdF/DlrgsIv/+F2D0mzGRKlYEMETmaI5Z6uDWZDBiwiBRS1QvZFrQb3LIW1hz6W/T9s3L0My6Bvh64lacxXg/y88C/W9TkLAwR2YUjghNzuDWZDBiwiNQqshp+PVm8SL8pwzHlcwbGYufZbGTlivviNkylFg9WAODGXY1xFiZM7YUn42rztw0isgm5BCecTaHyMGARqfgXTlknDSX2jYFnFSVm9ovBuBWHbL6em5VbiPe2nDE+Lj770iqyGg5m3ORSEhGZJZfgxBzOplB5GLBY6b1hLfDmxtMmX+Qlv9jim4Zj6VMtMXPdCWTlSpP7ApjOvigVQPHjikrOxjCgIXIPJQOTuKggbE7LknwHj6U4m0JiMWCxUrfG1fF405oVJsbGNw1H95gwfLj17P3aK9IqebZiydmYigKauKggAGDCL5ETsWTWpGRenKNxNoXEYsBSCSqlAm2jgy1qN6lbfTQMq+rw324qCmgCfT0AmObQWBrUlLzG2Rwi27N2SceRwUp4gDdm9G6Man5e/H5AVmPAYkeG2ZaU9BvYnJaFr3ZdEH14otTMfVOzJKgxd63kbE5Z37QA00DnX7X8jdey84r4zY1cnrnlG5VSIet8k7JwqYekwoDFzgyzMm2jgxEXFVRqxkVu07bmmOufuWslZ3Mycwow/rvDJtfMBTrVq1bBa7HA01/vR+H94w/K+w2trG/2RHJkSRASHuCNfs3Dse5opqyDE3O41ENSkVXA8tFHH+Hdd99FVlYWmjdvjiVLliAuLs7R3TJLqxNw4Nz1Sv2QLD7jItfEOKmZDX7yS18zF+yU9U3d0pkcBjxkDXOfF6DszxYAJB3PRGiAn8UzJJk5Bfh0R7rk76WyuNRD9iSbgOWHH37A1KlT8cknn6BNmzZYvHgxevTogdOnTyM0NNTR3Sulx+IdyLj5YOdPuJW/VZjLgym5dLT2yBWTU6FJr6xv6pbO5IgNeMoKbgDLAyExbflN334s/bcy98tEeZ+t5ONXMLUR8MpPx4yzhc6KwQk5mmwClkWLFuG5557D2LFjAQCffPIJNm7ciK+++grTp083aVtYWIjCwgfBQm5uLgBAo9FAo5F2OWXziSsAgJt38uGlevCFevNOPiZ/fxDvDWuBbo2r2+S1WtdWo3VtNV55vD4OZtxE9p1CZFzPw5qDf+PqbdedffFSCib/tYX8Qn3A56V6cO3GnXws33Xe7PUpPxwyeX6Y2hu9mlZHUupVk2KAgT73f1jla2zaNkztjek9G6Fzw1Djv31IVS+0iqwGAKWuGYIeS6/rtEUAIPnXi62Zey+A5eNhru0fp69h7qZTFf5bBfp4GB8X/7yU99mS4rNsD2Fqb7zSoyECfT1LjaueGgCg0xZBp3VcP4EHn2Fn+yw7EynHWMw9FYIgOPwr6d69e/D19cWaNWswYMAA4/XRo0fj1q1b+N///mfSfubMmZg1a1ap+3z33Xfw9fWVurtERERkA3l5eRg+fDhycnKgVqvLbSuLGZbs7GxotVpUr246M1G9enWcOnWqVPtXX30VU6dONT7Ozc1FREQEHn/88QrfsCinkoAticDtTOOly4owHG32JmYcUKJQZ34q9KvRDxunlKVW8rfIFhGBOHLpVpmzMSV37siRl1LA7Na6cseYKqfkGJc302OL2SJbzE6Zy22SOzl9lsubtXPmZR2NRoPNmzeje/fu8PDwcHR3XJKUY2xYIbGELAIWsby8vODl5VXquoeHh+0GM20d8NNolNx0HK7IwlEAHXAIG7WtzT41O6/Ibl84HgDaNzAN9Io/ntC1YZm1US5k5+H7lIumPxjMrMc7SqFO4fTr/nJnGOOMm4VY+ufF+1cfjHlZ16/eKbK4rS3uoW/rvJ8Fe3+WK8o3Kfk9wxXY9Ps/mSXFGIu5nywClpCQEKhUKly9etXk+tWrVxEWFmb/Dum0QPI0mKuQYvhFZHqVVdikaQkdlKXahPp7S9xBy5lL6i3+OKFLvXKTDS0Naiypw0JEtsdkWHIXsghYPD090apVK/z+++/GHBadTofff/8dCQkJ9u9Qxm4g90q5TcIUNxCnPIW9uhjjNcNJzfZaDrKFsqr1iglqLK10W1a5cKBE8OPjAcDBmXxEDmTNDjYiVyeLgAUApk6ditGjR6N169aIi4vD4sWLcffuXeOuIbu6c7XiNgBCccv4/8VPana1bx6WBDVlXSv5uEfT0nVngNKVbn9N3oSvRj9srHTrasW1iAzKC0JeiW/M7e5E98kmYBk2bBj++ecfvPHGG8jKykKLFi2QnJxcKhHXLqpa9prXEGj8f1Z3tIwlwY9hm1tcVJDJ+qa5YMfcN3VLZ3IY8JAYhpLzJatRV1SHBbhrcl3MDIml55URuQPZBCwAkJCQ4JgloJIi2wHqGkBuJso66Ufwr4FJT47Gk3c1/M3HTsr65m3uuiUzOWICnrKCGzGBkNi2t/I0sjtrytWICUIMv5SYq04NmP9sTelaD78mb8L8Qc0QGuDH7xNElSCrgEU2lCogfh6wehRQ6keG/puNottMtK0vvwq8pCdmGcvSgKesKXrAskBIbFuxVVXN5Tu4ygySNbMblgSM1gQhgGXLoQCM7XvFhnMHC1ElMWApS0w/YOg3+t1CxRNw/e8v+TTq5Zh+kV2Imc0BLA+ExLQt66wpwPKgp6wZpHmbTqDkUoWYoMcWgZOYe1gTWIgJGMv6d+FyDJF8MGApT0w/oFFv/a6hO1f1uS01HgaSf3V0z8hN2CJAMne9S8Ngs0sVlgY9tgicrLlHWe+9sgEjEckfA5aKKFVAVIcHj3leBbmAspYqpJxZssU9iMh9MWCxUkr6DeOWWybSERERSYsBi0hbTuprtDz99X5jqe1wbmkmIiKSVOm68lSm5NRMTPnhSKnrWTkFGLfiEJJTM0s/iYiIiCqNAYuFtDoBs9anma2JYbg2a30atDw8h4iIyOYYsFgoJf1GubUsBACZOQVISb9hv04RERG5CQYsFrp227LCW5a2IyIiIssxYLFQqL+3TdsRERGR5RiwWCguKgjhAd4oa/OyAvrdQoZiWERERGQ7DFgspFIqkNg3xuzfGYKYxL4xrMdCREQkAdZhsYROC2TsRrxwFcu7VMHVO6ZBSRjrsBAREUmKAUtF0taZHIDYSumNpOaf4ZeuOTgT1ImVbomIiOyAS0LlSVsHrB5lelrzffV3/hf9PQ+ibXQwgxUiIiKJMWApi06rn1kxWyruvuTp+nZEREQkKQYsZcnYbXZm5QEByL2sb0dERESSYsBSljtXbduOiIiIrMaApSxVq9u2HREREVmNAUtZItsB6hpAeaXi1DX17YiIiEhSDFjKolQB8fPuPygdtAgQcDFqCE9nJiIisgMGLOWJ6QcM/QZQly4IpwBQ++hiZL/VAId//dr+fSMiInIjDFgqEtMPmJyKMzEvQhAAocSEykPCdTTfPZFBCxERkYQYsFhAqxOgTvsOAgBFidUhQ8248D2zoC0qsnvfiIiI3AEDFguc2vcrquM6yipoq1QAYbiOU/t+tW/HiIiI3AQDFgvk37xs03ZEREQkDgMWC/hUq2nTdkRERCQOAxYLNGrTA7fgXyrh1kAnAFkIRqM2PezbMSIiIjfBgMUCqr+SEIDbZv9OEPRbnDPbJkJVpYp9O0ZEROQm+BO2IvdPbVYA5oveKoB7noH4V/en7NwxIiIi98EZlopUcGqzAoDXvVs8tZmIiEhCDFgqwlObiYiIHI4BS0V4ajMREZHDMWCpCE9tJiIicjgGLBUp99Tm+4/j5+rbERERkSQYsFiirFOb1TX012P6OaZfREREboLbmi0V0w9o1Bs4vws4cRMY/iNQtz1nVoiIiOyAMyxiKFVAZFv9/0e2ZbBCRERkJwxYiIiISPYYsBAREZHsMWAhIiIi2WPAQkRERLLHgIWIiIhkjwELERERyR4DFiIiIpI9BixEREQkey5R6VYQBABAbm6u5K+l0WiQl5eH3NxceHh4SP567ohjLD2OsX1wnKXHMZaelGNs+Llt+DleHpcIWG7fvg0AiIiIcHBPiIiISKzbt28jICCg3DYKwZKwRuZ0Oh2uXLkCf39/KBQlT1S2rdzcXERERODSpUtQq9WSvpa74hhLj2NsHxxn6XGMpSflGAuCgNu3b6NGjRpQKsvPUnGJGRalUolatWrZ9TXVajW/OCTGMZYex9g+OM7S4xhLT6oxrmhmxYBJt0RERCR7DFiIiIhI9hiwiOTl5YXExER4eXk5uisui2MsPY6xfXCcpccxlp5cxtglkm6JiIjItXGGhYiIiGSPAQsRERHJHgMWIiIikj0GLERERCR7DFiIiIhI9hiwmPHRRx+hTp068Pb2Rps2bZCSklJu+x9//BGNGjWCt7c3YmNjkZSUZKeeOi8xY/z555+jQ4cOqFatGqpVq4Zu3bpV+G9C4j/HBqtWrYJCocCAAQOk7aALEDvGt27dwoQJExAeHg4vLy80aNCA3y8sIHacFy9ejIYNG8LHxwcRERGYMmUKCgoK7NRb57Njxw707dsXNWrUgEKhwNq1ayt8zrZt29CyZUt4eXmhXr16WL58ueT9hEAmVq1aJXh6egpfffWVcOLECeG5554TAgMDhatXr5ptv2vXLkGlUgnz588X0tLShP/7v/8TPDw8hOPHj9u5585D7BgPHz5c+Oijj4TDhw8LJ0+eFMaMGSMEBAQIf//9t5177jzEjrFBenq6ULNmTaFDhw5C//797dNZJyV2jAsLC4XWrVsLvXr1Enbu3Cmkp6cL27ZtE44cOWLnnjsXseO8cuVKwcvLS1i5cqWQnp4u/Prrr0J4eLgwZcoUO/fceSQlJQmvv/668PPPPwsAhF9++aXc9ufPnxd8fX2FqVOnCmlpacKSJUsElUolJCcnS9pPBiwlxMXFCRMmTDA+1mq1Qo0aNYR33nnHbPuhQ4cKvXv3NrnWpk0b4YUXXpC0n85M7BiXVFRUJPj7+wtff/21VF10etaMcVFRkdCuXTvhiy++EEaPHs2ApQJix3jp0qVC3bp1hXv37tmriy5B7DhPmDBB6NKli8m1qVOnCu3bt5e0n67CkoDllVdeEZo0aWJybdiwYUKPHj0k7JkgcEmomHv37uHgwYPo1q2b8ZpSqUS3bt2wZ88es8/Zs2ePSXsA6NGjR5nt3Z01Y1xSXl4eNBoNgoKCpOqmU7N2jN98802EhobimWeesUc3nZo1Y7xu3Tq0bdsWEyZMQPXq1dG0aVPMmTMHWq3WXt12OtaMc7t27XDw4EHjstH58+eRlJSEXr162aXP7sBRP/dc4rRmW8nOzoZWq0X16tVNrlevXh2nTp0y+5ysrCyz7bOysiTrpzOzZoxLmjZtGmrUqFHqC4b0rBnjnTt34ssvv8SRI0fs0EPnZ80Ynz9/Hlu3bsWIESOQlJSEs2fPYvz48dBoNEhMTLRHt52ONeM8fPhwZGdn49FHH4UgCCgqKsJ//vMfvPbaa/boslso6+debm4u8vPz4ePjI8nrcoaFnMrcuXOxatUq/PLLL/D29nZ0d1zC7du3MXLkSHz++ecICQlxdHdclk6nQ2hoKD777DO0atUKw4YNw+uvv45PPvnE0V1zKdu2bcOcOXPw8ccf49ChQ/j555+xceNGzJ4929Fdo0riDEsxISEhUKlUuHr1qsn1q1evIiwszOxzwsLCRLV3d9aMscGCBQswd+5cbNmyBc2aNZOym05N7BifO3cOFy5cQN++fY3XdDodAKBKlSo4ffo0oqOjpe20k7HmcxweHg4PDw+oVCrjtcaNGyMrKwv37t2Dp6enpH12RtaM84wZMzBy5Eg8++yzAIDY2FjcvXsXzz//PF5//XUolfw9vbLK+rmnVqslm10BOMNiwtPTE61atcLvv/9uvKbT6fD777+jbdu2Zp/Ttm1bk/YAsHnz5jLbuztrxhgA5s+fj9mzZyM5ORmtW7e2R1edltgxbtSoEY4fP44jR44Y//Tr1w+dO3fGkSNHEBERYc/uOwVrPsft27fH2bNnjcEgAPz1118IDw9nsFIGa8Y5Ly+vVFBiCBIFnvVrEw77uSdpSq8TWrVqleDl5SUsX75cSEtLE55//nkhMDBQyMrKEgRBEEaOHClMnz7d2H7Xrl1ClSpVhAULFggnT54UEhMTua25AmLHeO7cuYKnp6ewZs0aITMz0/jn9u3bjnoLsid2jEviLqGKiR3jixcvCv7+/kJCQoJw+vRpYcOGDUJoaKjw1ltvOeotOAWx45yYmCj4+/sL33//vXD+/Hnht99+E6Kjo4WhQ4c66i3I3u3bt4XDhw8Lhw8fFgAIixYtEg4fPixkZGQIgiAI06dPF0aOHGlsb9jW/PLLLwsnT54UPvroI25rdpQlS5YItWvXFjw9PYW4uDhh7969xr/r2LGjMHr0aJP2q1evFho0aCB4enoKTZo0ETZu3GjnHjsfMWMcGRkpACj1JzEx0f4ddyJiP8fFMWCxjNgx3r17t9CmTRvBy8tLqFu3rvD2228LRUVFdu618xEzzhqNRpg5c6YQHR0teHt7CxEREcL48eOFmzdv2r/jTuKPP/4w+z3WMK6jR48WOnbsWOo5LVq0EDw9PYW6desKy5Ytk7yfCkHgHBkRERHJG3NYiIiISPYYsBAREZHsMWAhIiIi2WPAQkRERLLHgIWIiIhkjwELERERyR4DFiIiIpI9BixEREQkewxYiIiISPYYsBCRrIwdOxb/93//5+huEJHMVHF0B4iIDLRaLTZs2ICNGzc6uitEJDOcYSEiu/j+++/h4+ODzMxM47WxY8eiWbNmyMnJAQDs3r0bHh4eePjhh41tLl68iNGjR6N69erw8fFB8+bNsXPnTrv3n4gciwELEdnFE088gQYNGmDOnDkAgMTERGzZsgWbNm1CQEAAAGDdunXo27cvFAoFACAjIwNxcXHIz8/HunXrcOzYMSQkJECtVjvsfRCRY/C0ZiKymw0bNmDw4MGYMWMGFi5ciD///BNNmjQx/n2DBg3w3nvvoXfv3gCAXr16wdPTE2vXrnVQj4lILhiwEJFdtWzZEidOnMBvv/2Gjh07Gq+fPHkSDz/8MLKzs+Ht7Y2MjAzUqVMHhw8fRosWLRzXYSKSBS4JEZHdJCcn49SpU9BqtahevbrJ361btw7du3eHt7c3AODIkSPw9PQsM1hZsWIF4uLiEBsbi969e6OwsBAA0L59e+zbtw8A8Mwzz+C9996T7g0Rkd0wYCEiuzh06BCGDh2KL7/8El27dsWMGTNM/v5///sf+vfvb3zs4eGBoqIi5OXlmb1fz549kZKSguPHj6NGjRrYtm0bAGDGjBmYO3cuFi1aBKVSiSlTpkj2nojIfrgkRESSu3DhAtq2bYtJkyZh+vTp2LdvH9q2bYsDBw6gZcuWuHbtGmrVqoUrV64gJCQEAHD9+nXUr18fffv2xfTp0yEIAnbs2IGuXbuiXr16mDdvHn766Sfcu3cPly5dwi+//GJcYmrVqhUCAwORnJwMDw8PR751IrIRzrAQkaRu3LiB+Ph49O/fH9OnTwcAtGnTBj179sRrr70GAFi/fj3i4uKMwQoABAcHY/369Thz5gwefvhhPProo1i3bh1CQ0OxfPlynDp1Cjt27MDRo0dRrVo1xMTEAAD279+PGzduICAggMEKkQth4TgiklRQUBBOnTpV6nrx4nD/+9//0K9fv1Jt2rdvj927d5e6fuLECbRv3x4+Pj746KOPkJeXh4ceegiXL1/Gs88+i61bt2LQoEFITU1F06ZNbfuGiMghOMNCRA736KOP4sknn7S4/ciRIzF//nw88sgjSE9PR2xsLPLz8zFkyBAsWbIEUVFRePXVVzF79mwJe01E9sQcFiIiIpI9zrAQERGR7DFgISIiItljwEJERESyx4CFiIiIZI8BCxEREckeAxYiIiKSPQYsREREJHsMWIiIiEj2GLAQERGR7DFgISIiItljwEJERESy9/92hhHJQyIktgAAAABJRU5ErkJggg=="
many's avatar
many committed
     },
     "metadata": {},
     "output_type": "display_data"
many's avatar
many committed
    }
   ],
many's avatar
many committed
   "source": [
    "case.blade_loading()"
   ],
many's avatar
many committed
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
Malte Nyhuis's avatar
Malte Nyhuis committed
   "execution_count": 4,
   "outputs": [
    {
Malte Nyhuis's avatar
Malte Nyhuis committed
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2.6096421199420345\n"
     ]
    }
   ],
   "source": [
Malte Nyhuis's avatar
Malte Nyhuis committed
    "avdr =case.compute_avdr_inout_massave()\n",
Malte Nyhuis's avatar
Malte Nyhuis committed
    "print(avdr)\n"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
many's avatar
many committed
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}