{ "cells": [ { "cell_type": "markdown", "id": "6eb96853", "metadata": {}, "source": [ "# Ground Models" ] }, { "cell_type": "markdown", "id": "6751ef47", "metadata": {}, "source": [ "The ground models in PyWake are used to model the effects that the ground has on the inflow and wake." ] }, { "cell_type": "markdown", "id": "25037577", "metadata": {}, "source": [ "**Install PyWake if needed**" ] }, { "cell_type": "code", "execution_count": 1, "id": "a8c3019b", "metadata": {}, "outputs": [], "source": [ "# Install PyWake if needed\n", "try:\n", " import py_wake\n", "except ModuleNotFoundError:\n", " !pip install git+https://gitlab.windenergy.dtu.dk/TOPFARM/PyWake.git" ] }, { "cell_type": "markdown", "id": "8a11a3c3", "metadata": {}, "source": [ "### Mirror\n", "The Mirror ground model lets the ground mirror the wake deficit. It is implemented by adding wakes from underground mirrored wind turbines" ] }, { "cell_type": "code", "execution_count": 2, "id": "78e44d5d", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 6, "id": "58d63f34", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'y [m]')" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAksAAAGwCAYAAAC5ACFFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy89olMNAAAACXBIWXMAAA9hAAAPYQGoP6dpAABiPElEQVR4nO3deVxUVf8H8M/MwAz7IMgiiIqauK8pobn0ZC5ZPZaPlVlqLpmBZrQolWsZ7lpWWj0KlfmzzaXsycIFjdxyK1FzV5DFHUZAZrn3/v5AJkeGgVlgYPi8X6/7aubec8/9ngnhO+eee45MkiQJRERERGSW3NkBEBEREdVkTJaIiIiILGCyRERERGQBkyUiIiIiC5gsEREREVnAZImIiIjIAiZLRERERBa4OTuA2kYURWRnZ8PX1xcymczZ4RARUQ0mSRJu3ryJsLAwyOVV1z9RXFwMnU5ndz1KpRIeHh4OiMi1MFmyUnZ2NiIiIpwdBhER1SKZmZlo2LBhldRdXFyMJpGRuJSba3ddoaGhOHfuHBOmuzBZspKvry8A4MjJ8/D19XNo3aLo/MnUq2M+d8HBF3HE5ybeFZPKTQF3NzmUbrxTTUS2u3lTg3siI4x/O6qCTqfDpdxcnDybCT8/2/8uaTQatGgaAZ1Ox2TpLkyWrFR6683X18+uH0pzBCZLNqmKZKmUwl0BlbscblXYfU5Erq86hm34+Tn+7xKV4F8AIguK9QLyi/TQGURnh0JERE7CniWiSrhZrIebQg5fDzfIObCfiGogUZLK7SWv7PlkHnuWiCrJIIi4UahDodbg7FCIiKgasWeJyErFegEGUYLKTQ4Pd4WzwyEioirGZInIBgZBhM4gQhAlDgAnInJx/A1PZIcinYAbhXoU6wVnh0JEdVzpmCV7NmsIgoBp06YhMjISnp6eaNasGd555x1IFupJS0tDjx49EBgYCE9PT7Rs2RJLliwxKZOYmIiuXbvC19cXwcHBGDx4ME6cOGHTZ+Io7FkicoCbxQbc0gnw9XRjLxMR1Qnz5s3D8uXL8fnnn6NNmzbYv38/nn/+eajVakyaNMnsOd7e3oiLi0P79u3h7e2NtLQ0jB8/Ht7e3njhhRcAADt27EBsbCy6du0Kg8GAN998E/369cOxY8fg7e1dnU00kkmWUkAqQ6PRQK1W43zOdc6zZKPaNM+S5XPM7/d0V8DHg99DiKjkb0ZofTXy8/OrbA6k0r9LFy/fsHtSyobB9Sod6yOPPIKQkBCsXLnSuG/IkCHw9PTE6tWrK33dJ554At7e3vjyyy/NHr9y5QqCg4OxY8cO9OrVq9L1OhK/AhM52C29gBuFOt6aI6JaSaPRmGxardZsue7du2Pr1q04efIkAODPP/9EWloaBg4cWOlrHTp0CLt27ULv3r3LLZOfnw8ACAgIsKIVjsWvv0RVwCBKuFlsgEGQ4KHkAHAiqnqiWLLZcz6AMuufzpgxAzNnzixTfurUqdBoNGjZsiUUCgUEQcCcOXMwfPjwCq/VsGFDXLlyBQaDATNnzsTYsWPLiUnE5MmT0aNHD7Rt29bqNjkKkyWiKnRLL+CWXoCvhxunGSCiWiEz03SNOZVKZbbcN998g6+++gpr1qxBmzZtcPjwYUyePBlhYWEYOXKkxWv89ttvKCgowJ49ezB16lQ0b94cw4YNK1MuNjYW6enpSEtLs69RdmKyRFQNOACciGqLyq4x9/rrr2Pq1Kl4+umnAQDt2rXDhQsXkJiYWGGyFBkZaTzn0qVLmDlzZplkKS4uDps2bcLOnTvRsGFDG1vjGPytTVRNdAYR127qUFDMGcCJqPYrKiqC/K4vfwqFAqKV9wJFUTQZFyVJEuLi4rB+/Xps27bNmFg5E3uWiKpZodYArUGAt4q35ojIcap7bbhHH30Uc+bMQaNGjdCmTRscOnQIixcvxujRo41lEhISkJWVhS+++AIA8NFHH6FRo0Zo2bIlAGDnzp1YuHChyVQDsbGxWLNmDTZu3AhfX1/k5uYCANRqNTw9PW1unz2YLBE5gUGQkF+kh1Ypwlul4K05Iqp1li1bhmnTpuGll17C5cuXERYWhvHjx2P69OnGMjk5OcjIyDC+F0URCQkJOHfuHNzc3NCsWTPMmzcP48ePN5ZZvnw5AKBPnz4m10tKSsKoUaOqtE3l4TxLVuI8S/Zz9XmWyi1v4QQfDzd4q/jdhcjVVOc8S/b+XdJoNGjSIKBKY62t+HWWqAYoKDbgWoEWOoMdz/0SEVGV4FdZohrCIEi4UaiDh1IBXw83yGUyZ4dERLWIKEp23aFwRC+9q2LPElENU6wTcK2AM4ATEdUUTJaIaiC9IOJ6oQ43CnUw2DMlLxER2Y3JElENdksv4LJGy7mZiIiciGOWiGoBTbGey6YQkUWSZN8TzXw2vnzsWSKqJUpvzeXf0ts18RwREVmHyRJRLVOoNSA3vxiFWt6aIyKqDkyWiGqpG4V6XLmp5QBwIqIqxjFLRLWYVi8iN08LX083qD3dnR0OETmRIEl2rZDg6NUVXAmTJSIXcPOWATqDCC+lgsumEBE5GH+rErkIrV6EVi9CZxDh6+nGxXmJiByEv02JXEyhVkBunpYDwImIHIQ9S0Qu6nqBDgXFBgT6KtnLRFQHiKJk1/puXBuufPwNSuTCdAYROTeKcaNQ5+xQiIhqLfYsEdUBBcUlA8C9VW7w8eA/eyIia/C3JlEdoTOI0Bl00BoE+Hq4Q+nGjmUiospgskRUxxRpBdwsNqC+j4q9TEQuRJQku5ZC4jJK5eNXS6I66mqBFhev34LOwBnAiYgsYbJEVIcZRBHZebdw5abW2aEQEdVY7IMnIhRqDdALJQPAuWwKEZEpJktEBOCfAeCG20mTh7vC2SERkRVEqWSz53wyj7fhiMhEfpEB2TeKUVBs4IBPIiIwWSKiclzWaJHFAeBEREyWiKh8ekHCxeu3cCm/2NmhEBE5DccsEVGFbhYboBduwVulQD1vpbPDISIzRMnOteF4271cTJaIqFKK9QKK9QIMogQvpQLeKv76IKK6gbfhiMgq+UV65OQVI/+WHgaR45mIyPUxWSIim1zRaJFxlQPAicj1MVkiIpsZRBHnrxbi4vVbzg6FqM4TJMnuzRpNmjSBTCYrs8XGxpotr9frMXv2bDRr1gweHh7o0KEDNm/ebFImMTERXbt2ha+vL4KDgzF48GCcOHHC5s/EUZgsEZHdinQGnL9ayGVTiOqQP/74Azk5OcYtJSUFADB06FCz5d9++2188sknWLZsGY4dO4YXX3wRjz/+OA4dOmQss2PHDsTGxmLPnj1ISUmBXq9Hv379UFhYWC1tKo9Mkjj83RoajQZqtRrnc67Dz8/PoXULNWD61Or4abD220tF7Hn6w1iHDTFZe1lr4qzMZ1TRcCFL/7Qt1W+u3vI+H3P71Z7u8FAquGwKEUr+ZoTWVyM/P9/hfzPuvIZarcahM7nw9bX9GjdvatCpWajNsU6ePBmbNm3CqVOnIJPJyhwPCwvDW2+9ZdLzNGTIEHh6emL16tVm67xy5QqCg4OxY8cO9OrVy+qYHIWPsxCRQ90o1AOFegCAp7sCSjd2YBPVJhqNxuS9SqWCSqWyeI5Op8Pq1asRHx9vNlECAK1WCw8PD5N9np6eSEtLK7fe/Px8AEBAQEBlQq8y/C1GRFUi50Yxzl4uRLFecHYoRHWCJJb0DNu6Sbd7lSMiIqBWq41bYmJihdfesGED8vLyMGrUqHLL9O/fH4sXL8apU6cgiiJSUlKwbt065OTkmC0viiImT56MHj16oG3btrZ8JA7DniUiqlJnLxfCS6VAk/rezg6FiCohMzPT5DZcRb1KALBy5UoMHDgQYWFh5ZZ5//33MW7cOLRs2RIymQzNmjXD888/j1WrVpktHxsbi/T0dIs9T9WFPUtEVOWKtALOXi5Edh6fmiOq6fz8/Ey2ipKlCxcuYMuWLRg7dqzFckFBQdiwYQMKCwtx4cIF/P333/Dx8UHTpk3LlI2Li8OmTZuwfft2NGzY0K72OAJ7loioWpTOAA6UjGXisilEriEpKQnBwcEYNGhQpcp7eHggPDwcer0e33//PZ588knjMUmSMHHiRKxfvx6pqamIjIysqrCtwmSJiKpVXqEe10UdAMBTqYCHu8LJERG5BkmSLD4FW5nzrSWKIpKSkjBy5Ei4uZmmFCNGjEB4eLhxzNPevXuRlZWFjh07IisrCzNnzoQoinjjjTeM58TGxmLNmjXYuHEjfH19kZubCwBQq9Xw9PS0uW32qlW34Xbu3IlHH30UYWFhkMlk2LBhg8nxUaNGlZkca8CAASZlrl+/juHDh8PPzw/+/v4YM2YMCgoKqrEVRAQAWddv4XRuAQeAE9ViW7ZsQUZGBkaPHl3mWEZGhsng7eLiYrz99tto3bo1Hn/8cYSHhyMtLQ3+/v7GMsuXL0d+fj769OmDBg0aGLevv/66OppTrlrVs1RYWIgOHTpg9OjReOKJJ8yWGTBgAJKSkozv777XOnz4cOPkWXq9Hs8//zxeeOEFrFmzpkpjJyLzTucWQOkuR4tQX2eHQkRW6tevX7k9UqmpqSbve/fujWPHjlmsr6ZO/VirkqWBAwdi4MCBFsuoVCqEhoaaPXb8+HFs3rwZf/zxB+69914AwLJly/Dwww9j4cKFFkfxE1HVKdaJOJl7E57uCkQEejk7HCIiE7XqNlxlpKamIjg4GFFRUZgwYQKuXbtmPLZ79274+/sbEyUA6Nu3L+RyOfbu3Wu2Pq1WC41GY7IRkeMV60TcKNQj81oRl00hskF1rw1Xl7hUsjRgwAB88cUX2Lp1K+bNm4cdO3Zg4MCBEISSMRG5ubkIDg42OcfNzQ0BAQHGQWR3S0xMNJmcKyIiosrbQVSXXSvQIev6LVwr0KKg2ODscIiIatdtuIo8/fTTxtft2rVD+/bt0axZM6SmpuLBBx+0qc6EhATEx8cb32s0GiZMRNUg81rJnEz3hPrAU6mAvJwlFIiIqppL9SzdrWnTpqhfvz5Onz4NAAgNDcXly5dNyhgMBly/fr3ccU4qlarMBF1EVH2OZ9/EsYu8/U1EzuPSydLFixdx7do1NGjQAAAQExODvLw8HDhwwFhm27ZtEEUR0dHRzgqTiCpQbBBxLEuD05c4zQdReexZF650I/Nq1W24goICYy8RAJw7dw6HDx9GQEAAAgICMGvWLAwZMgShoaE4c+YM3njjDTRv3hz9+/cHALRq1QoDBgzAuHHjsGLFCuj1esTFxeHpp5/mk3BENVyRToBGq4f8sgweSjnC/J03QR0R1S21Klnav38/HnjgAeP70rFEI0eOxPLly/HXX3/h888/R15eHsLCwtCvXz+88847JnMtffXVV4iLi8ODDz4IuVyOIUOG4IMPPqj2thCRbS4XFAMAFDIZlG5yLptCRFWuViVLffr0sThh1S+//FJhHQEBAZyAksgFnLlackuufbg/lG5yKN1celQBETlRrUqWiIju9ldWHuQyGe5rGujsUIicSpQkiHbMlWTPua6OX8WIqNbTCSL+ysjHsSw+NUdEjseeJSJyCXlaHQDg9KUCuClkaFLf28kREZGrYLJERC4lI68IAKB0k8NdIUeQr6qCM4iILGOyREQuKT03HwDQvUl9uClk8HBXODkioqrFMUtVh2OWiMil7Tp/FTvPXHF2GERUizFZIiKXJ0rAnrPXsP/cDWeHQkS1EG/DEVGdkKfVwyCKOJalgJtChhahvs4OiYhqCSZLRFSnnM0rgEGS4OmugEwmQ8MALptCrsHe9d24Nlz5mCwRUZ10IOcGBFGCtyoY7go5fDz465CIzOOYJSKq07aduYz/ncjhk0BEVC4mS0REALacuIRtJy47OwwiqoHY70xEhNIB4BIOXyhZa659I7WzQyKyCudZqjpMloiI7nDsugZ6QYSXSgG5TIamwVw2haiuY7JERGTG75lXoRcl+HqGQS4DAn24bApRXcUxS0REFmw8no3v0rNh4HPVRHUWe5aIiCph3ZEsAMCTHSKcHAmReaIoQRTtGLNkx7mujskSEVElXC7QAwC2nbgMN7kMve4JcnJERFRdmCwREVnhz0s3YRAlBHgoIZfL0Drcz9khEVEVY7JERGSDn89cgUGU4OfpDneFDCFqD2eHRERVhAO8iYjskHQgE5/uy0Ch1uDsUKiOkyRAtGOzdpqlJk2aQCaTldliY2PNlk9OTi5T1sPjny8Zer0eU6ZMQbt27eDt7Y2wsDCMGDEC2dnZ9nwsDsGeJSIiB/h493m4yWV4pVczZ4dCVC3++OMPCIJgfJ+eno6HHnoIQ4cOLfccPz8/nDhxwvheJpMZXxcVFeHgwYOYNm0aOnTogBs3buDll1/GY489hv3791dNIyqJyRIRkQPc1ArQCxK++TMTXm4KPNImzNkhEVWpoCDThxzmzp2LZs2aoXfv3uWeI5PJEBoaavaYWq1GSkqKyb4PP/wQ3bp1Q0ZGBho1amR/0DZiskRE5ED7MwuMA8C93NzQsbG/s0MisopGozF5r1KpoFJZnpRVp9Nh9erViI+PN+ktultBQQEaN24MURTRuXNnvPfee2jTpk255fPz8yGTyeDv729VGxyNY5aIiKrAd0cvY9XBizh9qQCZ14qcHQ7VAaVrw9mzAUBERATUarVxS0xMrPDaGzZsQF5eHkaNGlVumaioKKxatQobN27E6tWrIYoiunfvjosXL5otX1xcjClTpmDYsGHw83PuU6fsWSIiqkKL087BIIiYN6gV6nkrnR0OUYUyMzNNkpOKepUAYOXKlRg4cCDCwsq//RwTE4OYmBjj++7du6NVq1b45JNP8M4775iU1ev1ePLJJyFJEpYvX25DKxyLyRIRUTWY8tNx+Hi4Y/FjrZ0dCpFFfn5+VvXkXLhwAVu2bMG6deusuo67uzs6deqE06dPm+wvTZQuXLiAbdu2Ob1XCeBtOCKiaiGKEq4XaLEg9TQ+2XPO2eEQOUxSUhKCg4MxaNAgq84TBAFHjhxBgwYNjPtKE6VTp05hy5YtCAwMdHS4NmHPEhFRNfozMx8GQUKIjxJ+7u74V1Sws0MiF2GQJBjsWN/NYO1ESwBEUURSUhJGjhwJNzfTlGLEiBEIDw83jnmaPXs27rvvPjRv3hx5eXlYsGABLly4gLFjxwIoSZT+85//4ODBg9i0aRMEQUBubi4AICAgAEql825jM1kiInKCr/ZlQxBF1PdUwUulQPMQH2eHRGS1LVu2ICMjA6NHjy5zLCMjA3L5Pzewbty4gXHjxiE3Nxf16tVDly5dsGvXLrRuXXJrOisrCz/88AMAoGPHjiZ1bd++HX369KmydlREJkk2pJJ1mEajgVqtxvmc6w6/jyrUgBWfq+OnQXDwRRyxUrZoyzcqK0+xJs7KfEaiaPm4pX/aluo3V295n4+5/Xefb76MdFeZ8svf+U3ZIIlmy5mUuf3atGzpMfH2sZIdpf/mSj+P0nP1Qkk5felxEdALpcf++fZuECUYSveLEvSCBEG847gg3i4jQry93yCIMNwuJ4glrw2CiP8+0wlh/p5lPiuq3TQaDULrq5Gfn19lY29K/y79fOA8vH1sv0ZhgQYDuzSp0lhrK/Ys2cgRf6CJiEqNXn0QXh5uWDc22tmhUC1lkCSTLwm2nE/mcYC3jSSpZvQEEZHryNNoMXbtYcT/cMzZoRDRHdizZAdJAgyCBJkMUMjLn7GUiKiyTmXmQRQlLPBTIlytxDOdnLfEAxGVYLLkAHf2MjFpIiJH+O638xAEEfU9lQj18kT7Rmpnh0RUZ/E2nINI0j9JE2/PEZGjTPnqTwz76HecvlTg7FCqxP9+2oQObaLQrvU9SFr1X2eHU6s5arkTKos9Sw5W+rMmiBJ7mYjIYR5dkAq5Qo5DcwZA6eYa33MNBgOmvh6Pn1O2Q61Wo0d0Fzz278drzESERKVc419cDVQ6nomIyFF0xTr0mZ+KJ/6719mhOMT+P/ahVes2CA8Ph4+PD/oNGIitKb86OyyiMtizVMU4AJyIHOlSdh6yLlzFKC8lGvh7IvHhls4OyWY52dkICw83vg8LC0d2dpYTIyIyj8lSNeBTc0TkaL+mHINBb0C4Wokm/h54pE35q71T3WAQ7VzuhONty8XbcNWIA8CJyNFmfrgDI6ZvxP5zN5wditUahIUhO+ufnqTs7Cw0aMCkj2oeJkvVjE/NEVFV6PdyMgIGvIcbhTpnh1Jp93bthmPH0pGVlYWCggL8uvln9O3X39lhEZXBZMlJ7kyaiIgcpdnzX6Bb7AJ4KWX44vMkZ4djkZubGxLnLcLAhx7AfV07YtIrr/JJOKqROGbJye4cz0REZLdLZyH4eQAANPn5Tg6mYo88+hgeefQxZ4fhEgQ7xyzxy3v52LNUQ5T2NBEREVHNwmSphmHCREREVLMwWaqB+PgmERFRzcExSzVUacLkxnmZiIioEgySCIMk2nU+mceepRrOIEq8NUdERORETJZqAUGSOACciIjISXgbrpYQbmdKCpTcluNUA0RERNWjVvUs7dy5E48++ijCwsIgk8mwYcMGk+OSJGH69Olo0KABPD090bdvX5w6dcqkzPXr1zF8+HD4+fnB398fY8aMQUFBQTW2wj6CJBl7moiIiEqJkv0bmVerkqXCwkJ06NABH330kdnj8+fPxwcffIAVK1Zg79698Pb2Rv/+/VFcXGwsM3z4cBw9ehQpKSnYtGkTdu7ciRdeeKG6muAwgmTf5GNERERUObXqNtzAgQMxcOBAs8ckScLSpUvx9ttv49///jcA4IsvvkBISAg2bNiAp59+GsePH8fmzZvxxx9/4N577wUALFu2DA8//DAWLlyIsLCyCzhqtVpotVrje41GUwUtsx2fmiMiIqpatapnyZJz584hNzcXffv2Ne5Tq9WIjo7G7t27AQC7d++Gv7+/MVECgL59+0Iul2Pv3r1m601MTIRarTZuERERVdsQG/GpOSIioqrhMslSbm4uACAkJMRkf0hIiPFYbm4ugoODTY67ubkhICDAWOZuCQkJyM/PN26ZmZlVEL1j8Kk5IqK6yyCKdm9kXq26DecMKpUKKpXK2WFUGp+aIyIiciyX6VkKDQ0FAFy6dMlk/6VLl4zHQkNDcfnyZZPjBoMB169fN5ZxFaVPzREREZF9XCZZioyMRGhoKLZu3Wrcp9FosHfvXsTExAAAYmJikJeXhwMHDhjLbNu2DaIoIjo6utpjrg4GkU/NERER2aNW3YYrKCjA6dOnje/PnTuHw4cPIyAgAI0aNcLkyZPx7rvv4p577kFkZCSmTZuGsLAwDB48GADQqlUrDBgwAOPGjcOKFSug1+sRFxeHp59+2uyTcK6ET80REbk2gyTBYMcdBXvOdXW1Klnav38/HnjgAeP7+Ph4AMDIkSORnJyMN954A4WFhXjhhReQl5eH+++/H5s3b4aHh4fxnK+++gpxcXF48MEHIZfLMWTIEHzwwQfV3hZnMYgSFDIZxzIREZFdmjRpggsXLpTZ/9JLL5mdD/Ho0aOYPn06Dhw4gAsXLmDJkiWYPHmySZnExESsW7cOf//9Nzw9PdG9e3fMmzcPUVFR5cZRmgtY4+2330ZAQECly9eqZKlPnz6QLGS+MpkMs2fPxuzZs8stExAQgDVr1lRFeLWGXhDhrii5A8ukiYiIbPHHH39AEATj+/T0dDz00EMYOnSo2fJFRUVo2rQphg4dildeecVsmR07diA2NhZdu3aFwWDAm2++iX79+uHYsWPw9vY2e87SpUsRExMDpVJZqbjT0tIQFxfnuskSOY5eKHlElEkTERHd6e7Jl8t7KjwoKMjk/dy5c9GsWTP07t3bbL1du3ZF165dAQBTp041W2bz5s0m75OTkxEcHIwDBw6gV69e5ca8fv36MlMDlcfX17dS5e7kMgO8yTZ6QeRTc0RELkAUJQh2bOLtsa0REREmkzEnJiZWeG2dTofVq1dj9OjRkDnw23d+fj4AWOwFSkpKglqtrnSdn3zySZk5GSvCniWCKEow3H6tYBcTEVGdlpmZCT8/P+P7ysw1uGHDBuTl5WHUqFEOi0MURUyePBk9evRA27Ztyy03cuRIq+p95plnrI6FyRIBgPEbRWlfI5MmIqK6yc/PzyRZqoyVK1di4MCBDn2yPDY2Funp6UhLS6v0OZmZmZDJZGjYsCEAYN++fVizZg1at26NF154weZYeBuOTIh3dMUSERFV5MKFC9iyZQvGjh3rsDrj4uKwadMmbN++3Zj4VMYzzzyD7du3AyhZ4uyhhx7Cvn378NZbb1l8+KsiTJbILL0gGgeBExFRzVe6coM9my2SkpIQHByMQYMG2d0GSZIQFxeH9evXY9u2bYiMjLTq/PT0dHTr1g0A8M0336Bt27bYtWsXvvrqKyQnJ9scF2/DkUV3PzVHRERUShRFJCUlYeTIkXBzM00pRowYgfDwcOMAcZ1Oh2PHjhlfZ2Vl4fDhw/Dx8UHz5s0BlNx6W7NmDTZu3AhfX1/jIvdqtRqenp4VxqPX641jrLZs2YLHHnsMANCyZUvk5OTY3E7+BaRKYU8TERHdbcuWLcjIyMDo0aPLHMvIyDBJULKzs9GpUyd06tQJOTk5WLhwITp16mRy+2758uXIz89Hnz590KBBA+P29ddfVyqeNm3aYMWKFfjtt9+QkpKCAQMGGK8dGBhoczvZs0RWKe2m5QBwIiLq169fuZNFp6ammrxv0qSJxYmlAVR4vCLz5s3D448/jgULFmDkyJHo0KEDAOCHH34w3p6zBZMlsgqfmiMiqpnsXTi9Ni+6XlRUBC8vL/Tp0wdXr16FRqNBvXr1jMdfeOEFeHl52Vw/b8ORTfjUHBER1RT169fHI488gk8//RRXr141SZSAkl6tys7wbQ6TJbKL1iBwLBMRETnV33//jf79++Obb75B48aNER0djTlz5uDIkSMOqZ/JEtlNlCQOACciIqdp1KgRJk6ciC1btuDSpUuYPHkyjhw5gp49e6Jp06aYPHkytm3bZrLwrzWYLJFDiJLEpImIyIn0omj8HWzTJrrG7261Wo1hw4Zh7dq1uHLlClasWAFBEPD8888jKCgIX331ldV1coA3OZR4+0kGPjVHRETO5u7ujn79+qFfv35YtmwZDh06BIPBUPGJd2GyRFXCUNq7dHsySyZNRERUHYqLi/HXX3/h8uXLEO/oLZPJZHj00UdtqpPJElWpO5MmhUxm83T6REREFdm8eTNGjBiBq1evljkmk8k4ZolqNoMgQmuw7YeUiIgqphclu7fabuLEiRg6dChycnIgiqLJZmuiBDBZomrG+ZmIiKiqXLp0CfHx8QgJCXFovUyWqFqVPjXHpImIiBztP//5T5llVhyBY5bIKUqfmpODA7+JiMgxPvzwQwwdOhS//fYb2rVrB3d3d5PjkyZNsqleJkvkVKUDwN0U7OQkIrKHKAL2THPnCtMs/d///R9+/fVXeHh4IDU1FbI7nsSWyWRMlqj2EiUmTUREZL+33noLs2bNwtSpUyGXO+7vCZMlqhFKhy8xaSIiIlvpdDo89dRTDk2UAA7wphpGlG5vHABORERWGjlyJL7++muH18ueJaqRDLcTpdIfULmcA8GJiCzRCxIUgu1fMvV2nFtTCIKA+fPn45dffkH79u3LDPBevHixTfUyWaIa7c6kiQkTERFZcuTIEXTq1AkAkJ6ebnJMZseyW0yWqFYw3L4/p3TjnWMiIjJv+/btVVIvkyWqVXQGEXIZB4AT3UkURUB7E9AXAwYtRH3Jvw+tTgdJkuz6Rk1EHOBNtVDpVAMGQeQgcCIA0BVAt30GdGmJ0P2xGCf+bw4AYNqbU1BUVOTk4Ki66AXJ7q02euKJJ6DRaCpdfvjw4bh8+bJV12DPEtVKpU/MuSkAOb81ExHVWRs3bsSVK1cqVVaSJPz444945513EBwcXOlrMFmiWs0gSAA4lonqOIXS+FLZZRIiW/jixP/NwTvvzYOXl5cTAyOqepIkoUWLFlV6DSZL5BJ0BhFyuQxufGKO6iCTJR0U7pC7lyRPKqWS45XI5dkyqDs8PNyq8kyWyGWIogTD7ddMmoiorjGIknG6FVvPr4169+5d5ddgskQupXTAtwGAXMbxTEREZD8mS+SSRFGCCMBNUfKeSRMREdmKyRK5NMPtR2H51BwREdmKjxBRnWAQJOgMorPDICKqMoL4z7glWzbBhl+RWVlZePbZZxEYGAhPT0+0a9cO+/fvL7f8qFGjIJPJymxt2rQxlklMTETXrl3h6+uL4OBgDB48GCdOnLDlI3EYJktUZwhSScJUWwcxEhHVJDdu3ECPHj3g7u6On3/+GceOHcOiRYtQr169cs95//33kZOTY9wyMzMREBCAoUOHGsvs2LEDsbGx2LNnD1JSUqDX69GvXz8UFhZWR7PMqtRtuICAAKsqlclkOHjwIBo3bmxTUERVRZAkCIIEUZRxqgEiIjvMmzcPERERSEpKMu6LjIy0eI5arYZarTa+37BhA27cuIHnn3/euG/z5s0m5yQnJyM4OBgHDhxAr169ytTZqVOnSk+RcfDgwUqVu1ulkqW8vDwsXbrUpIHlkSQJL730EgRBsCkgoupQmjQBcjBfIiL6x91Lh6hUKqhUqjLlfvjhB/Tv3x9Dhw7Fjh07EB4ejpdeegnjxo2r9LVWrlyJvn37Wuxcyc/PB1B+x83gwYONr4uLi/Hxxx+jdevWiImJAQDs2bMHR48exUsvvVTpuO5W6QHeTz/9dKWnBp84caLNARFVJ/3tm/TuCiZNRFS7GUQJCjvWdysdohAREWGyf8aMGZg5c2aZ8mfPnsXy5csRHx+PN998E3/88QcmTZoEpVKJkSNHVni97Oxs/Pzzz1izZk25ZURRxOTJk9GjRw+0bdvWbJkZM2YYX48dOxaTJk3CO++8U6ZMZmZmhTGVp1LJkihaN+rr5s2bNgVD5Cx6QYRCJoObghkTEdVtmZmZ8PPzM74316sElOQG9957L9577z0AJbfD0tPTsWLFikolS59//jn8/f1NeobuFhsbi/T0dKSlpVUq9m+//dbsAPNnn30W9957L1atWlWpeu7GAd5EtwmSBK1BNE43QERUF/n5+Zls5SVLDRo0QOvWrU32tWrVChkZGRVeQ5IkrFq1Cs899xyUSqXZMnFxcdi0aRO2b9+Ohg0bVip2T09P/P7772X2//777/Dw8KhUHebYNM9SdnY20tLScPny5TK9TpMmTbI5GKKaQJAkQADkchlvzRERlaNHjx5lHuk/efJkpR7u2rFjB06fPo0xY8aUOSZJEiZOnIj169cjNTW1wkHjd5o8eTImTJiAgwcPolu3bgCAvXv3YtWqVZg2bVql67mb1clScnIyxo8fD6VSicDAQNMFHGUyJkvkEvSCBAgS3BUyJk1EVCvoRQlyO6ZG0Vt57iuvvILu3bvjvffew5NPPol9+/bh008/xaeffmosk5CQgKysLHzxxRcm565cuRLR0dFmxyHFxsZizZo12LhxI3x9fZGbmwug5Ek6T09PizFNnToVTZs2xfvvv4/Vq1cDKOntSkpKwpNPPmlV++5kdbI0bdo0TJ8+HQkJCZDLeRePXFtp0qRy5wBwIqI7de3aFevXr0dCQgJmz56NyMhILF26FMOHDzeWycnJKXNbLj8/H99//z3ef/99s/UuX74cANCnTx+T/UlJSRg1alSFcT355JN2JUbmWJ0sFRUV4emnn2aiRHWKVi9CLudTc0REd3rkkUfwyCOPlHs8OTm5zD61Wo2ioqJyz5Ek+8aN5uXl4bvvvsPZs2fx2muvISAgAAcPHkRISAjCw8NtqtPqjGfMmDH49ttvbboYUW0miiVJk5bLphAR1Uh//fUXWrRogXnz5mHBggXIy8sDAKxbtw4JCQk212t1z1JiYiIeeeQRbN68Ge3atYO7u7vJ8cWLF9scDFFtIIqA1lAy1YCc3UxEVEPoBQlyO57m1bvAk8Dx8fEYNWoU5s+fD19fX+P+hx9+GM8884zN9dqULP3yyy+IiooCgDIDvInqAlEEBEmEQiGDu4K3pImIaoI//vgDn3zySZn94eHhxoHitrD6t/yiRYuwatUqHD9+HKmpqdi+fbtx27Ztm82BOMLMmTPLrGTcsmVL4/Hi4mLExsYiMDAQPj4+GDJkCC5duuTEiKm2EwQJxToBIhfnJSJyOpVKVWa5FqBkSoOgoCCb67U6WVKpVOjRo4fNF6xqbdq0MVnR+M5ZP1955RX8+OOP+Pbbb7Fjxw5kZ2fjiSeecGK05Cp0BhFavcikiYjIiR577DHMnj0ber0eQMkdr4yMDEyZMgVDhgyxuV6rk6WXX34Zy5Yts/mCVc3NzQ2hoaHGrX79+gBKHlVcuXIlFi9ejH/961/o0qULkpKSsGvXLuzZs6fc+rRaLTQajclGZI4kSdAZRON6c0RE1UkUJQh2bK7wZW/RokUoKChAcHAwbt26hd69e6N58+bw9fXFnDlzbK7X6jFL+/btw7Zt27Bp0ya0adOmzADvdevW2RyMI5w6dQphYWHw8PBATEwMEhMT0ahRIxw4cAB6vR59+/Y1lm3ZsiUaNWqE3bt347777jNbX2JiImbNmlVd4ZMLEAQJovjPVANERFQ91Go1UlJSkJaWhr/++gsFBQXo3Lmzyd9+W1idLPn7+9fYW1fR0dFITk5GVFQUcnJyMGvWLPTs2RPp6enIzc2FUqmEv7+/yTkhISEWB30lJCQgPj7e+F6j0ZRZkZnobpIkQRBKFposnQWciIiqx/333497770XKpXKIQ+fWZ0sJSUl2X3RqjJw4EDj6/bt2yM6OhqNGzfGN998U+EU6eVRqVTlLiJIVBGDKMIgAko3ORR8WpSIqEqJoog5c+ZgxYoVuHTpEk6ePImmTZti2rRpaNKkidm16CrDpe8R+Pv7o0WLFjh9+jRCQ0Oh0+mME1SVunTpEkJDQ50TINUZOoMInUEsWaSXiKgKGETJ7q22e/fdd5GcnIz58+dDqVQa97dt2xb//e9/ba63UslS586dcePGjUpXev/99yMrK8vmoByloKAAZ86cQYMGDdClSxe4u7tj69atxuMnTpxARkYGYmJinBgl1RWCcQB47f+FRERUE33xxRf49NNPMXz4cCgUCuP+Dh064O+//7a53krdhjt8+DD+/PNPBAQEVKrSw4cPQ6vV2hyUrV577TU8+uijaNy4MbKzszFjxgwoFAoMGzYMarUaY8aMQXx8PAICAuDn54eJEyciJiam3MHdRFWh9Nacm1wOdwVvzREROUpWVhaaN29eZr8oisbpBGxR6TFLDz74YKUXt3PWTN4XL17EsGHDcO3aNQQFBeH+++/Hnj17jBNRLVmyBHK5HEOGDIFWq0X//v3x8ccfOyVWIoMoQpJkUChkHM9EROQArVu3xm+//YbGjRub7P/uu+/QqVMnm+utVLJ07tw5qytu2LCh1efYa+3atRaPe3h44KOPPsJHH31UTRERWaYXJOgFCUo39jIRkX0Mggi5HfO8GVxgjrjp06dj5MiRyMrKgiiKWLduHU6cOIEvvvgCmzZtsrneSiVLd2doRORYOsM/t+bkLv3YBRFR1fn3v/+NH3/8EbNnz4a3tzemT5+Ozp0748cff8RDDz1kc71WTx1ARFVDFAHd7XmZFOxlIiKySc+ePZGSkuLQOpksEdUwd96aYy8TEZH19u/fj+PHjwMoGcfUpUsXu+pjskRUQxXrBbgpZLw1R0SVYhAlyO2YK8kV5lkqfdDr999/N67YkZeXh+7du2Pt2rU2j6fmr2CiGswgSCjWCxA4NxMRUYXGjh0LvV6P48eP4/r167h+/TqOHz8OURQxduxYm+u1OlkaOXIkdu7cafMFich6WoOIYp14e5HekvFNRERkaseOHVi+fDmioqKM+6KiorBs2TK7cherk6X8/Hz07dsX99xzD957770aMVM3UV0gShK0BhEGZkpERGZFRESYnXxSEASEhYXZXK/VydKGDRuQlZWFCRMm4Ouvv0aTJk0wcOBAfPfdd3bNjklElaMziCjWC0yaiMiEQRDt3mq7BQsWYOLEidi/f79x3/79+/Hyyy9j4cKFNtdr05iloKAgxMfH488//8TevXvRvHlzPPfccwgLC8Mrr7yCU6dO2RwQEVVMvL3OXLFeqPTM+kRErm7UqFE4fPgwoqOjoVKpoFKpEB0djYMHD2L06NEICAgwbtaw62m4nJwcpKSkICUlBQqFAg8//DCOHDmC1q1bY/78+XjllVfsqZ6IKiBKEm7pBbjJ5VC68XkNIqrbli5dWiX1Wp0s6fV6/PDDD0hKSsKvv/6K9u3bY/LkyXjmmWfg5+cHAFi/fj1Gjx7NZImomhhEEaJeghvXmSOiOmzkyJFVUq/VyVKDBg0giiKGDRuGffv2oWPHjmXKPPDAA8b5DYioehgECQahJGFSsZeJqM6RJAmiHXMlucIt/YMHD8Ld3R3t2rUDAGzcuBFJSUlo3bo1Zs6cCaVSaVO9Vv9GXbJkCbKzs/HRRx+ZTZQAwN/f36bFd4nIfgZBQqFWgJ5zMxFRHTN+/HicPHkSAHD27Fk89dRT8PLywrfffos33njD5nqtTpaee+45eHh42HxBIqoeOoOIW3oBggt8WyQiqoyTJ08aO3K+/fZb9O7dG2vWrEFycjK+//57m+tlXz2RCzMIEoq0ArSG2v9IMBFRRUpuRZb8vtuyZQsefvhhACXzL129etXmepksEdUBOoOIwttJE3uaiFyTQZTs3qyVlZWFZ599FoGBgfD09ES7du1M5ji6W1paGnr06GEs37JlSyxZssSkTGJiIrp27QpfX18EBwdj8ODBOHHiRKXiuffee/Huu+/iyy+/xI4dOzBo0CAAwLlz5xASEmJ1+0pxIV2iOqJkbiYJosTFeYnIfjdu3ECPHj3wwAMP4Oeff0ZQUBBOnTqFevXqlXuOt7c34uLi0L59e3h7eyMtLQ3jx4+Ht7c3XnjhBQAlS5bExsaia9euMBgMePPNN9GvXz8cO3YM3t7eFmNaunQphg8fjg0bNuCtt95C8+bNAQDfffcdunfvbnNbZZIrDH+vRhqNBmq1GmcuXoPv7akSHK2qv/nb87REmbocEKut4VjbDms+18pMjl2ZfzoVXdPSdcr7bMvdb6au8stK8FAq4KYoO82AVfWY2X/3+ebLSHeVKb/8nd92DZJotpxJmduvTcuWHhNvHyvZIdw+UPr/qfRc/e2ZjPWlx0UYB8zrhX++gRvEkicQS8vqBQnCHd/QDYJ4u4wIsfTbuyDCcLucIIq3n2IsOa7TCRDuei0IJa9FQYRBbzD+V6/TG98LRfnQ/lwyVctvR7PRsWkwCgsL4eXlBTc3fid2Jo1Gg9D6auTn5xun16mKa6jVajy+fCfcPX1srkd/qwDrJ/SqdKxTp07F77//jt9++83mawLAE088AW9vb3z55Zdmj1+5cgXBwcHYsWMHevXqZdM1iouLoVAo4O7ubtP5/G5JVEcV6wQUaQUuyutiWob5QaFQwM/Pj4kS2USj0ZhsWq3WbLkffvgB9957L4YOHYrg4GB06tQJn332mVXXOnToEHbt2oXevXuXWyY/Px8ArJ51+04eHh42J0oAkyWiOk0UJRRpDSjWiRDFyvWoUc3y5H/uxWsv/svZYVAN4Ki14SIiIqBWq41bYmKi2eudPXsWy5cvxz333INffvkFEyZMwKRJk/D5559XGGvDhg2hUqlw7733IjY2FmPHjjVbThRFTJ48GT169EDbtm1t/3DsxK8dRASdQYQolUxo6cbBTLXCPa3DofbzwMdD2qGwsBBznB0QuYzMzEyT23AqlcpsOVEUce+99+K9994DAHTq1Anp6elYsWJFhTNp//bbbygoKMCePXswdepUNG/eHMOGDStTLjY2Funp6UhLS7OjRfZjskREAP6ZAVzpBrgpZJBz2ZQaK7RhAHa83sfZYZCL8vPzq9SYpQYNGqB169Ym+1q1alWp+YwiIyMBAO3atcOlS5cwc+bMMslSXFwcNm3ahJ07d6Jhw4ZWtMDxmCwRkQmdQYTOAHgo5XCTyx0yiJ8cZ8fsgWgY4OnsMIjQo0ePMo/0nzx5Eo0bN7aqHlEUTcZFSZKEiRMnYv369UhNTTUmVs7EZImIzCrSCpDLRHipFM4OhQCsnHAfAjxVTJSoXAZBgsyOZY4MVp77yiuvoHv37njvvffw5JNPYt++ffj000/x6aefGsskJCQgKysLX3zxBQDgo48+QqNGjdCyZUsAwM6dO7Fw4UJMmjTJeE5sbCzWrFmDjRs3wtfXF7m5uQAAtVoNT8+yP//x8fGVjnnx4sVWtbEUkyUiKpcoSSgoNsBNLoOHkkmTM7w4qAXCfVXo3qy+s0MhMtG1a1esX78eCQkJmD17NiIjI43zHJXKyclBRkaG8b0oikhISMC5c+fg5uaGZs2aYd68eRg/fryxzPLlywEAffr0MbleUlISRo0aVSaOQ4cOmbw/ePAgDAYDoqKiAJT0dikUCnTp0sXmtjJZIqIKGcSSZVPcFDIo3TgAvDp0aR0Ctac7xt/n/FsQROV55JFH8Mgjj5R7PDk52eT9xIkTMXHiRIt1Wjv94/bt242vFy9eDF9fX3z++efGyTFv3LiB559/Hj179rSq3jsxWSKiStEJInRCyXQDbgrOAF5VPD3d4eulxMdD2jk7FKJaZ9GiRfj1119NZhGvV68e3n33XfTr1w+vvvqqTfUyWSIiqxQbRMAgwqucGcDJPv/3fFfU81Y6OwyqhQRRgtyOFRoEB67u4CwajQZXrlwps//KlSu4efOmzfUyWSIim2i0erjJ5PBRubGXyQGWPtEe3ioFEyUiOzz++ON4/vnnsWjRInTr1g0AsHfvXrz++ut44oknbK6XyRIR2cwgicgr1sFDoeBTczZ6qVcT+Cnd0DLM19mhENV6K1aswGuvvYZnnnkGer0eAODm5oYxY8ZgwYIFNtfLZImI7FZkMMAgSlAq5BwAXkm9owIR4OWGAa1CnR0Kkcvw8vLCxx9/jAULFuDMmTMAgGbNmsHb29uuepksEZFD6EQBOlGAl+QGN7kMbgomTeYE+Kjg7+WOuB5NnR0KuRhRlCDYscCj6AJjlkp5e3ujffv2DquPyRIROVRpL5O/SskB4GYkPtwS3ir+6iWqCoWFhZg7dy62bt2Ky5cvQ7wreTx79qxN9fJfLBFViTytDgAQ4KGEXM6kacZDLaBQyJgoEVWhsWPHYseOHXjuuefQoEEDyBy0xiX/1RJRlbpcpIWHmxz+HnXzKa/RnRvCS6lARKCXs0Mhcnk///wzfvrpJ/To0cOh9TJZIqIqV2wQcbVIC6VCDh/3uvFrp1dTNTwUCtwbWa/iwkQOYBAkoBrXhquJ6tWrh4CAAIfXyxGYRFQtDKIEjU6PAr2hZGJLF9XAT4moYE8MbhvOJ92Iqtk777yD6dOno6ioyKH11o2veERUY2h0JXOfBEAFN7kMcgeNKagpJt4fCTfO0knkFIsWLcKZM2cQEhKCJk2awN3d3eT4wYMHbaqXyRIROcXVYi0AINTLwyUSptj7GkPpJmeiROREgwcPrpJ6mSwRkVPlFhVDLgNCvTydHYpNhrZpALkcCFF7ODsUquMMoggItt/iNtgxR1NNMWPGjCqpl8kSETmdThCRW3QLbjI5AmrJU3M9Ivwhl8m4TAlRHcBkiYhqBJ0gokgSAABuchn8lO4VnOEcjfxVUMhl6N6svrNDISIAAQEBOHnyJOrXr4969epZnFvp+vXrNl2DyRIR1Sia24tfAiVJk1JesxboHdK+obNDIKI7LFmyBL6+vsbXjpqI8k5MloioRrpyq2QAeLiPJ9xkchgk546nGN6xIWcipxpNFCW71nerrWvDjRw50vh61KhRVXINJktEVKNl3CyCUiFHqJdzBlD3bRoMhUyGet61YywVUV02YsQIPPDAA+jVqxeaNWvmsHr5jCsR1XglA8CLcfV2b1N16BJSD/c1CECT+t5cqoSollAqlUhMTMQ999yDiIgIPPvss/jvf/+LU6dO2VUve5aIqFbQCSKKbt8mcJPL4aesml9foV4ecJPL0DZCXSX1E1HV+e9//wsAyMrKws6dO7Fjxw4sWrQI48ePR4MGDXDx4kWb6mWyRES1ikZngF4QAXhALgM8FI4bAK5SyPGvqGCH1UdUnQRBhMyOeZYEO86taerVq4fAwEDUq1cP/v7+cHNzQ1BQkM31MVkiolrpUmEx9KKEpmpvuDlg4PXDUQ3gpuAAbqLa7M0330RqaioOHTqEVq1aoXfv3pg6dSp69eqFevVsX9SayRIR1Wpn8wsBAC3q+dh0fkyjQChkMvh48NchUW03d+5cBAUFYcaMGXjiiSfQokULh9RbZwd4f/TRR2jSpAk8PDwQHR2Nffv2OTskIrKRIALn8gtxseBWpc9pXd8PHUL8EebvyaVKiFzEoUOH8NZbb2Hfvn3o0aMHwsPD8cwzz+DTTz/FyZMnba63Tn6V+vrrrxEfH48VK1YgOjoaS5cuRf/+/XHixAkEB3O8AlFtVGyQcFNrAFAMN7kM/irzM4CHenlCqZCjRSiXKSHXIggSZILtcyUJdpxbU3To0AEdOnTApEmTAAB//vknlixZgtjYWIiiCEEQbKq3TiZLixcvxrhx4/D8888DAFasWIGffvoJq1atwtSpUytVR2FRIeQOHFhaEVGquh9iR05E5og4bQnH2jYIVsRZ2bUlJVius6JrWrpOeZ+rVM5+c3WZq8Pc53b3rrvPM9xV4O7JIu8sL9xR9s7zDHeUubO60roMZuoorffOesTbA1JL94kicFUHGAQJet+SeZG83f75NaeUy9CqSclTboWFhXAVd7bFldrlCor4/6NaSZKEQ4cOITU1FampqUhLS4NGo0H79u3Ru3dvm+utc8mSTqfDgQMHkJCQYNwnl8vRt29f7N69u0x5rVYLrfafuV00Gg0AoH2LRlUfLBGRlZo0DHF2CEROExAQgIKCAnTo0AG9e/fGuHHj0LNnT/j7+9tVb51Llq5evQpBEBASYvoLJSQkBH///XeZ8omJiZg1a1Z1hUdEREQ2Wr16NXr27Ak/Pz+H1lvnkiVrJSQkID4+3vheo9EgIiICx85kwtfXsf8zhBqwLk8V3u0zcvQtRUfcRrQlJmsva02clblNWNHtQUu3BS3VX9lbeID524B3n1+Z2393vrV02+/OW34V3e67+1ZfoIcSSoXcJSeblCQJRUVFAAAvL68qWUiUbHNTo0HTxmHVci1REO2aK0m08tyZM2eW6UyIiooy2/EAAHq9HomJifj888+RlZWFqKgozJs3DwMGDDCWSUxMxLp16/D333/D09MT3bt3x7x58xAVFVWpmAYNGmRVGyqrziVL9evXh0KhwKVLl0z2X7p0CaGhoWXKq1QqqFSqMvu9vLzh7e3t0NjqSrJkzXihymCyZF5545kqqt+aZMl8IlSZMo5Nlgxmk6U7ywKFAPIFEVk3RchlMjQPsW2qgZrKx8e12uMqbB1QXFu0adMGW7ZsMb53cys/rXj77bexevVqfPbZZ2jZsiV++eUXPP7449i1axc6deoEANixYwdiY2PRtWtXGAwGvPnmm+jXrx+OHTvm8L+51qhzyZJSqUSXLl2wdetWDB48GAAgiiK2bt2KuLg45wZHRFXu5LWbMEgSvFQKKGQyThtAdJfSsbmlyus0AEqSI3MdDeZ8+eWXeOutt/Dwww8DACZMmIAtW7Zg0aJFWL16NQBg8+bNJuckJycjODgYBw4cQK9evaxtisPUyXmW4uPj8dlnn+Hzzz/H8ePHMWHCBBQWFhqfjiMi17fv4nWkXbiKgmIDCrUGZ4dDZDdRlOzeACAiIgJqtdq4JSYmlnvNU6dOISwsDE2bNsXw4cORkZFRblmtVgsPD9MvJ56enkhLSyv3nPz8fAAlA7edqc71LAHAU089hStXrmD69OnIzc1Fx44dsXnz5jKDvonI9f1yMhcKuQyD24Y7OxSiGiEzM9NkgHR5vUrR0dFITk5GVFQUcnJyMGvWLPTs2RPp6enw9S07j1n//v2xePFi9OrVC82aNcPWrVuxbt26cm9ViqKIyZMno0ePHmjbtq1jGmcjmWRpYAOVodFooFarcT7nusNH23PMkm04Zsk8jlkyHbNUcsx0HqfSf3OCJKG+hwpuchl63WP7YptEd9NoNAitr0Z+fr7D/2bceQ21Wo2Ob/0IhYft43qE4kIcnvOozbHm5eWhcePGWLx4McaMGVPm+JUrVzBu3Dj8+OOPkMlkaNasGfr27YtVq1bh1q2ys+9PmDABP//8M9LS0tCwYUOb2uQodfI2HBHR3a4Wa3Gx4BbSM/ORnpnv7HCIah1/f3+0aNECp0+fNns8KCgIGzZsQGFhIS5cuIC///4bPj4+aNq0aZmycXFx2LRpE7Zv3+70RAlgskREZOKvq/k4cOkGzl8tROa1ImeHQ1RpgiDZvdmjoKAAZ86cQYMGDSyW8/DwQHh4OAwGA77//nv8+9//Nh6TJAlxcXFYv349tm3bhsjISLtichQmS0REZuw4fwW/nLmEG4U65N/SOzscohrntddew44dO3D+/Hns2rULjz/+OBQKBYYNGwYAGDFihMlqGXv37sW6detw9uxZ/PbbbxgwYABEUcQbb7xhLBMbG4vVq1djzZo18PX1RW5uLnJzc83epqtOdXKANxFRZX1/NAuCCIyLbgI5J3skMrp48SKGDRuGa9euISgoCPfffz/27NmDoKCScX8ZGRmQy//pkykuLsbbb7+Ns2fPwsfHBw8//DC+/PJLk6VIli9fDgDo06ePybWSkpIwatSoqm5SuZgsERFVwvojWXxqjugOa9eutXg8NTXV5H3v3r1x7Ngxi+fU1GfOmCwREVVCzk0dAGDXmauQy2S4r2mgkyMiMiWIImDHcidCRY/Y1mFMloiIrLA3Kx8GUYK/hxJuCtdbNoWIymKyRERkgx9PXoJBlDDaIwJuChkCfcxP3EdEtR+fhiMissOn+zLw8e4L0Bl4C4PIVbFniYjIAT7dex7uChnG31cz5oWhukcQ7ByzZMe5ro7JEhGRA1wu0EMvSNiQngUPhQIDWlVuJXYiqvmYLBEROdCu8xoYRAmhXh5QusnROrxq1gMjourDZImIqAp88Wc2DIKI+J6R8HRXIETt4eyQiMhGTJaIiKrQ/B1nYRBELPl3G3ir+CuXqo4gSIAd67vZuzacK+PTcERE1WDaLycx/ZcTzg6DiGzAZImIqBpoinTIuFaED38/izWHMpwdDhFZgX3CRETVaM/ZGzAIEgI8lPBXuXPZFKJagMkSEZETrNydCUEU8a7KHR7uCjQN9nZ2SFTLiXaOWRI5ZqlcvA1HROREUzcdQ9x3f+JagdbZoRBROZgsERHVAM99cQBPJ+93dhhEZAaTJSKiGkAUJVy+XoRJ69P51BxRDcMxS0RENcifp69CFCV8olaivpc7hrRv6OyQqJYQBAmSHeu7ccxS+ZgsERHVQKtSzkAQRDTw9kCApwotw3ydHRJRncXbcERENVhs8gEMef83XLx+y9mhENVZTJaIiGqBh97bik7TfnF2GER1Em/DERHVEgWaW/jX4p0IC/LG6ue6ODscqmFEQQTsGrNk+7mujskSEVEtcuHsVZw5kYtJPipE1PPA632aOzskIpfHZImIqBb67oc/YdAbEK5WoomfF7o3q+/skIhcFpMlIqJabNL8rRB1xUj7YDhah/s5Oxwil8RkiYjIBdw/YSWgcEf2NxPg4a5wdjjkBKIkAqIdY5YkjlkqD5+GIyJyFbc0aPrCWjz88S5nR0LkUpgsERG5kOKzR50dApHLYbJEREREZAHHLBEREbkAzrNUddizRERERFabOXMmZDKZydayZctyyycnJ5cp7+HhYTyu1+sxZcoUtGvXDt7e3ggLC8OIESOQnZ1dHc2xiD1LREREZJM2bdpgy5YtxvdubpbTCj8/P5w4ccL4XiaTGV8XFRXh4MGDmDZtGjp06IAbN27g5ZdfxmOPPYb9+/c7PngrMFkiIiIiI41GY/JepVJBpVKZLevm5obQ0NBK1y2Tycotr1arkZKSYrLvww8/RLdu3ZCRkYFGjRpV+jqOxttwRERELkAURLs3AIiIiIBarTZuiYmJ5V7z1KlTCAsLQ9OmTTF8+HBkZGRYjLGgoACNGzdGREQE/v3vf+PoUctPb+bn50Mmk8Hf39/qz8OR2LNERERERpmZmfDz+2c2+PJ6laKjo5GcnIyoqCjk5ORg1qxZ6NmzJ9LT0+Hr61umfFRUFFatWoX27dsjPz8fCxcuRPfu3XH06FE0bNiwTPni4mJMmTIFw4YNM4nHGZgsERERkZGfn1+lkpOBAwcaX7dv3x7R0dFo3LgxvvnmG4wZM6ZM+ZiYGMTExBjfd+/eHa1atcInn3yCd955x6SsXq/Hk08+CUmSsHz5cjta4xhMloiIiMhu/v7+aNGiBU6fPl2p8u7u7ujUqVOZ8qWJ0oULF7Bt2zan9yoBHLNERETkEgRBsHuzR0FBAc6cOYMGDRpUOt4jR46YlC9NlE6dOoUtW7YgMDDQrpgchckSERERWe21117Djh07cP78eezatQuPP/44FAoFhg0bBgAYMWIEEhISjOVnz56NX3/9FWfPnsXBgwfx7LPP4sKFCxg7diyAkkTpP//5D/bv34+vvvoKgiAgNzcXubm50Ol0TmljKd6GIyIiIqtdvHgRw4YNw7Vr1xAUFIT7778fe/bsQVBQEAAgIyMDcvk/fTI3btzAuHHjkJubi3r16qFLly7YtWsXWrduDQDIysrCDz/8AADo2LGjybW2b9+OPn36VEu7zGGyRERERFZbu3atxeOpqakm75csWYIlS5aUW75JkyaQJMkRoTkckyUiIiIXIIkSJNH2ZMOec10dxywRERERWcBkiYiIiMgCJktEREREFnDMEhERkQsQBQGSHXMl2XOuq2PPEhEREZEFLpUsNWnSBDKZzGSbO3euSZm//voLPXv2hIeHByIiIjB//nwnRUtERES1gcvdhps9ezbGjRtnfH/nyscajQb9+vVD3759sWLFChw5cgSjR4+Gv78/XnjhBWeES0RERDWcyyVLvr6+CA0NNXvsq6++gk6nw6pVq6BUKtGmTRscPnwYixcvZrJERES1miiKkAmizedLou3nujqXug0HAHPnzkVgYCA6deqEBQsWwGAwGI/t3r0bvXr1glKpNO7r378/Tpw4gRs3bpitT6vVQqPRmGxERERUd7hUz9KkSZPQuXNnBAQEYNeuXUhISEBOTg4WL14MAMjNzUVkZKTJOSEhIcZj9erVK1NnYmIiZs2aVfXBExERUY1U43uWpk6dWmbQ9t3b33//DQCIj49Hnz590L59e7z44otYtGgRli1bBq1Wa/P1ExISkJ+fb9wyMzMd1TQiIiKqBWp8z9Krr76KUaNGWSzTtGlTs/ujo6NhMBhw/vx5REVFITQ0FJcuXTIpU/q+vHFOKpUKKpXK+sCJiIiqkShIkMntGLMkcG248tT4ZCkoKAhBQUE2nXv48GHI5XIEBwcDAGJiYvDWW29Br9fD3d0dAJCSkoKoqCizt+CIiIiIavxtuMravXs3li5dij///BNnz57FV199hVdeeQXPPvusMRF65plnoFQqMWbMGBw9ehRff/013n//fcTHxzs5eiIiIqqpanzPUmWpVCqsXbsWM2fOhFarRWRkJF555RWTREitVuPXX39FbGwsunTpgvr162P69OmcNoCIiIjK5TLJUufOnbFnz54Ky7Vv3x6//fZbNURERERUfURRhMyOuZI4z1L5XOY2HBEREVFVYLJEREREZAGTJSIiIiILXGbMEhERUV0mCqKd8yxxzFJ52LNEREREZAGTJSIiIiILmCwRERERWcAxS0RERC5AFATI5ILN50uC7ee6OvYsEREREVnAZImIiIjIAiZLRERERBYwWSIiInIFgt7+zQozZ86ETCYz2Vq2bFlu+aNHj2LIkCFo0qQJZDIZli5dWqZMYmIiunbtCl9fXwQHB2Pw4ME4ceKEtZ+EwzFZIiIiIpu0adMGOTk5xi0tLa3cskVFRWjatCnmzp2L0NBQs2V27NiB2NhY7NmzBykpKdDr9ejXrx8KCwurqgmVwqfhiIiIyCZubm7lJj5369q1K7p27QoAmDp1qtkymzdvNnmfnJyM4OBgHDhwAL169bIvWDuwZ4mIiIiMNBqNyabVasste+rUKYSFhaFp06YYPnw4MjIyHBpLfn4+ACAgIMCh9VqLyRIREZErEAX7NwARERFQq9XGLTEx0ezloqOjkZycjM2bN2P58uU4d+4cevbsiZs3bzqmOaKIyZMno0ePHmjbtq1D6rQVb8MRERGRUWZmJvz8/IzvVSqV2XIDBw40vm7fvj2io6PRuHFjfPPNNxgzZozdccTGxiI9Pd3iOKjqwmSJiIiIjPz8/EySpcry9/dHixYtcPr0abtjiIuLw6ZNm7Bz5040bNjQ7vrsxdtwREREZLeCggKcOXMGDRo0sLkOSZIQFxeH9evXY9u2bYiMjHRghLZjzxIREZErEA1Wz5VU5nwrvPbaa3j00UfRuHFjZGdnY8aMGVAoFBg2bBgAYMSIEQgPDzeOedLpdDh27JjxdVZWFg4fPgwfHx80b94cQMmttzVr1mDjxo3w9fVFbm4uAECtVsPT09P2ttmJyRIRERFZ7eLFixg2bBiuXbuGoKAg3H///dizZw+CgoIAABkZGZDL/7mBlZ2djU6dOhnfL1y4EAsXLkTv3r2RmpoKAFi+fDkAoE+fPibXSkpKwqhRo6q0PZYwWSIiIiKrrV271uLx0gSoVJMmTSBJksVzKjruLByzRERERGQBe5aIiIhcgWAAZNaNOypzPpnFniUiIiIiC5gsEREREVnAZImIiIjIAo5ZIiIicgWCHpAp7DufzGLPEhEREZEFTJaIiIiILGCyRERERGQBxywRERG5AlEo2ew5n8xizxIRERGRBUyWiIiIiCxgskRERERkAccsERERuQKDAYAdcyUZuDZcedizRERERGQBkyUiIiIiC5gsEREREVnAMUtERESuQDDYuTYcxyyVhz1LRERERBYwWSIiIiKygMkSERERkQUcs0REROQKRD0g2NEHItoxR5OLY88SERERkQVMloiIiIgsYLJEREREZAHHLBEREbkCUQBkgn3nk1nsWSIiIiKygMkSERERkQVMloiIiIgsqDXJ0pw5c9C9e3d4eXnB39/fbJmMjAwMGjQIXl5eCA4Oxuuvvw6DwXStm9TUVHTu3BkqlQrNmzdHcnJy1QdPRERU1QShZH03mzf7xizNnTsXMpkMkydPLrdMnz59IJPJymyDBg0CAOj1ekyZMgXt2rWDt7c3wsLCMGLECGRnZ9sVm71qTbKk0+kwdOhQTJgwwexxQRAwaNAg6HQ67Nq1C59//jmSk5Mxffp0Y5lz585h0KBBeOCBB3D48GFMnjwZY8eOxS+//FJdzSAiInI5f/zxBz755BO0b9/eYrl169YhJyfHuKWnp0OhUGDo0KEAgKKiIhw8eBDTpk3DwYMHsW7dOpw4cQKPPfZYdTSjXLXmabhZs2YBQLk9Qb/++iuOHTuGLVu2ICQkBB07dsQ777yDKVOmYObMmVAqlVixYgUiIyOxaNEiAECrVq2QlpaGJUuWoH///tXVFCIiIpdRUFCA4cOH47PPPsO7775rsWxAQIDJ+7Vr18LLy8uYLKnVaqSkpJiU+fDDD9GtWzdkZGSgUaNGjg2+kmpNslSR3bt3o127dggJCTHu69+/PyZMmICjR4+iU6dO2L17N/r27WtyXv/+/S12GWq1Wmi1WuP7/Px8AMDNmxrHNgCAKEoOr9NaUjWEIDj4Io743EQbYrL2stZcozKfkSRWcNxCHZbqF83UW17s5vbffb75Mqb77ixyd3nDHe/vfH1nOcMd9Qm3Xxvu+IBKDxtuB1daT2kcpZ9HaT362+X0Ymk5QC/cPiZIxnKCeOc5EvSCBFH857hBEGEQJRgEEZIkGV8bBAmCeLusIMEgiiWvdQJEoeS1oBcgCCIEQYKoE0v2Gwwl/9UbIOn1kPQCJEGEpC8GDFrAUPJfQ3EhNBrH/44i65X+rbD079FhBB3suoqgA4AyPzsqlQoqlarc02JjYzFo0CD07du3wmTpbitXrsTTTz8Nb2/vcsvk5+dDJpOVOwSnOrhMspSbm2uSKAEwvs/NzbVYRqPR4NatW/D09CxTb2JiorFX607tWjRxUORERI6V9sdihE51dhR0p5s3b0KtVldJ3UqlEqGhocg9mmR3XT4+PoiIiDDZN2PGDMycOdNs+bVr1+LgwYP4448/rL7Wvn37kJ6ejpUrV5Zbpri4GFOmTMGwYcPg5+dn9TUcxanJ0tSpUzFv3jyLZY4fP46WLVtWU0RlJSQkID4+3vheFEVcv34dgYGBkMlkTovLXhqNBhEREcjMzHTqD6C9XKEdrtAGgO2oaVyhHa7QBkmScPPmTYSFhVXZNTw8PHDu3DnodDq765IkqczftvJ6lTIzM/Hyyy8jJSUFHh4eVl9r5cqVaNeuHbp162b2uF6vx5NPPglJkrB8+XKr63ckpyZLr776KkaNGmWxTNOmTStVV2hoKPbt22ey79KlS8Zjpf8t3XdnGT8/P7O9SoD57kdndgU6mp+fX639JXQnV2iHK7QBYDtqGldoR21vQ1X1KN3Jw8PDpoTFHgcOHMDly5fRuXNn4z5BELBz5058+OGH0Gq1UCgUZs8tLCzE2rVrMXv2bLPHSxOlCxcuYNu2bU7//+/UZCkoKAhBQUEOqSsmJgZz5szB5cuXERwcDABISUmBn58fWrdubSzzv//9z+S8lJQUxMTEOCQGIiKiuuLBBx/EkSNHTPY9//zzaNmyJaZMmVJuogQA3377LbRaLZ599tkyx0oTpVOnTmH79u0IDAx0eOzWqjVjljIyMnD9+nVkZGRAEAQcPnwYANC8eXP4+PigX79+aN26NZ577jnMnz8fubm5ePvttxEbG2vsGXrxxRfx4Ycf4o033sDo0aOxbds2fPPNN/jpp5+c2DIiIqLax9fXF23btjXZ5+3tjcDAQOP+ESNGIDw8HImJiSblVq5cicGDB5dJhPR6Pf7zn//g4MGD2LRpEwRBMI47DggIgFKprMIWla/WJEvTp0/H559/bnzfqVMnAMD27dvRp08fKBQKbNq0CRMmTEBMTAy8vb0xcuRIky6+yMhI/PTTT3jllVfw/vvvo2HDhvjvf/9bJ6cNUKlUmDFjhsUnHGoDV2iHK7QBYDtqGldohyu0oa7LyMiAXG46peOJEyeQlpaGX3/9tUz5rKws/PDDDwCAjh07mhwr/XvvDDKpWp5nJCIiIqqdas0M3kRERETOwGSJiIiIyAImS0REREQWMFkiIiIisoDJkgs7f/48xowZg8jISHh6eqJZs2aYMWOGySyv58+fh0wmK7Pt2bPHpK5vv/0WLVu2hIeHB9q1a1dmvipntwMA/vrrL/Ts2RMeHh6IiIjA/Pnzy9TlzHYAwJw5c9C9e3d4eXmVO7mpuf8fa9euNSmTmpqKzp07Q6VSoXnz5uUuMF1VKtOOjIwMDBo0CF5eXggODsbrr78Og8FgUsbZ7bhbkyZNynz2c+fONSlTmZ8zZ/voo4/QpEkTeHh4IDo6usyEvTXNzJkzy3zud67cUFxcjNjYWAQGBsLHxwdDhgwpM8EwUZWSyGX9/PPP0qhRo6RffvlFOnPmjLRx40YpODhYevXVV41lzp07JwGQtmzZIuXk5Bg3nU5nLPP7779LCoVCmj9/vnTs2DHp7bffltzd3aUjR47UmHbk5+dLISEh0vDhw6X09HTp//7v/yRPT0/pk08+qTHtkCRJmj59urR48WIpPj5eUqvVZssAkJKSkkz+f9y6dct4/OzZs5KXl5cUHx8vHTt2TFq2bJmkUCikzZs3V1MrKm6HwWCQ2rZtK/Xt21c6dOiQ9L///U+qX7++lJCQUKPacbfGjRtLs2fPNvnsCwoKjMcr83PmbGvXrpWUSqW0atUq6ejRo9K4ceMkf39/6dKlS84OrVwzZsyQ2rRpY/K5X7lyxXj8xRdflCIiIqStW7dK+/fvl+677z6pe/fuToyY6homS3XM/PnzpcjISOP70mTp0KFD5Z7z5JNPSoMGDTLZFx0dLY0fP76qwqzQ3e34+OOPpXr16klarda4b8qUKVJUVJTxfU1qR1JSksVkaf369eWe+8Ybb0ht2rQx2ffUU09J/fv3d2CElVNeO/73v/9Jcrlcys3NNe5bvny55OfnZ/x/VJPaUapx48bSkiVLyj1emZ8zZ+vWrZsUGxtrfC8IghQWFiYlJiY6MSrLZsyYIXXo0MHssby8PMnd3V369ttvjfuOHz8uAZB2795dTRFSXcfbcHVMfn4+AgICyux/7LHHEBwcjPvvv984IVip3bt3o2/fvib7+vfvj927d1dprJbc3Y7du3ejV69eJrO79u/fHydOnMCNGzeMZWpaO8oTGxuL+vXro1u3bli1ahWkO6ZDqw3t2L17N9q1a4eQkBDjvv79+0Oj0eDo0aPGMjWxHXPnzkVgYCA6deqEBQsWmNw6rMzPmTPpdDocOHDA5HOVy+Xo27ev0z/Xipw6dQphYWFo2rQphg8fjoyMDAAl64/p9XqTNrVs2RKNGjWq8W0i11FrZvAm+50+fRrLli3DwoULjft8fHywaNEi9OjRA3K5HN9//z0GDx6MDRs24LHHHgMA5ObmmvzRA4CQkBDjFPTVzVw7cnNzERkZaVKuNObc3FzUq1evxrWjPLNnz8a//vUveHl54ddff8VLL72EgoICTJo0CUD5/z80Gg1u3bpV7qLQ1am8GEuPWSrjzHZMmjQJnTt3RkBAAHbt2oWEhATk5ORg8eLFxpgr+jlzpqtXr0IQBLOf699//+2kqCoWHR2N5ORkREVFIScnB7NmzULPnj2Rnp6O3NxcKJXKMmPjauK/XXJd7FmqhaZOnWp2EPCd292/GLOysjBgwAAMHToU48aNM+6vX78+4uPjER0dja5du2Lu3Ll49tlnsWDBglrVDmeypR2WTJs2DT169ECnTp0wZcoUvPHGGzX2/0dtYE274uPj0adPH7Rv3x4vvvgiFi1ahGXLlkGr1Tq5Fa5t4MCBGDp0KNq3b4/+/fvjf//7H/Ly8vDNN984OzQiAOxZqpVeffVVjBo1ymKZpk2bGl9nZ2fjgQceQPfu3fHpp59WWH90dDRSUlKM70NDQ8s8eXLp0iWEhoZaF/hdHNmO8mIsPWapTHW3w1rR0dF45513oNVqoVKpym2Hn5+fXb0xjmxHaGhomSewKvv/w9523M2edkVHR8NgMOD8+fOIioqq1M+ZM9WvXx8KhaJKfs6rk7+/P1q0aIHTp0/joYcegk6nQ15enknvUm1rE9VuTJZqoaCgIAQFBVWqbFZWFh544AF06dIFSUlJZRY0NOfw4cNo0KCB8X1MTAy2bt2KyZMnG/elpKQgJibG6tjv5Mh2xMTE4K233oJer4e7u7sxxqioKOOtkZrQDlscPnwY9erVMy4mGhMTU2bKg5rWjpiYGMyZMweXL19GcHAwgJIY/fz80Lp1a2OZqmjH3exp1+HDhyGXy41tqMzPmTMplUp06dIFW7duxeDBgwEAoihi69atiIuLc25wVigoKMCZM2fw3HPPoUuXLnB3d8fWrVsxZMgQACULsWZkZDj8Z4WoXM4eYU5V5+LFi1Lz5s2lBx98ULp48aLJY7mlkpOTpTVr1kjHjx+Xjh8/Ls2ZM0eSy+XSqlWrjGV+//13yc3NTVq4cKF0/PhxacaMGdX6yH1l2pGXlyeFhIRIzz33nJSeni6tXbtW8vLyKjN1gDPbIUmSdOHCBenQoUPSrFmzJB8fH+nQoUPSoUOHpJs3b0qSJEk//PCD9Nlnn0lHjhyRTp06JX388ceSl5eXNH36dGMdpY/cv/7669Lx48eljz76qNofua+oHaVTB/Tr1086fPiwtHnzZikoKMjs1AHObMeddu3aJS1ZskQ6fPiwdObMGWn16tVSUFCQNGLECGOZyvycOdvatWsllUolJScnS8eOHZNeeOEFyd/f3+TJxJrm1VdflVJTU6Vz585Jv//+u9S3b1+pfv360uXLlyVJKpk6oFGjRtK2bduk/fv3SzExMVJMTIyTo6a6hMmSC0tKSpIAmN1KJScnS61atZK8vLwkPz8/qVu3biaP6Jb65ptvpBYtWkhKpVJq06aN9NNPP9WodkiSJP3555/S/fffL6lUKik8PFyaO3dumbqc2Q5JkqSRI0eabcf27dslSSqZU6pjx46Sj4+P5O3tLXXo0EFasWKFJAiCST3bt2+XOnbsKCmVSqlp06ZSUlJSjWqHJEnS+fPnpYEDB0qenp5S/fr1pVdffVXS6/U1qh13OnDggBQdHS2p1WrJw8NDatWqlfTee+9JxcXFJuUq83PmbMuWLZMaNWokKZVKqVu3btKePXucHZJFTz31lNSgQQNJqVRK4eHh0lNPPSWdPn3aePzWrVvSSy+9JNWrV0/y8vKSHn/8cZMvS0RVTSZJdzyTTEREREQm+DQcERERkQVMloiIiIgsYLJEREREZAGTJSIiIiILmCwRERERWcBkiYiIiMgCJktEREREFjBZIiIiIrKAyRIRWe38+fOQyWSQyWTo2LGj3fWV1nXnQqlERDUFkyUistmWLVuwdetWu+vJycnB0qVL7Q+IiKgKMFkiIpsFBgYiMDDQ7npCQ0OhVqsdEBERkeMxWSKq465cuYLQ0FC89957xn27du2CUqm0utdo1KhRGDx4MN577z2EhITA398fs2fPhsFgwOuvv46AgAA0bNgQSUlJjm4GEVGVcXN2AETkXEFBQVi1ahUGDx6Mfv36ISoqCs899xzi4uLw4IMPWl3ftm3b0LBhQ+zcuRO///47xowZg127dqFXr17Yu3cvvv76a4wfPx4PPfQQGjZsWAUtIiJyLPYsEREefvhhjBs3DsOHD8eLL74Ib29vJCYm2lRXQEAAPvjgA0RFRWH06NGIiopCUVER3nzzTdxzzz1ISEiAUqlEWlqag1tBRFQ12LNERACAhQsXom3btvj2229x4MABqFQqm+pp06YN5PJ/voeFhISgbdu2xvcKhQKBgYG4fPmy3TETEVUH9iwREQDgzJkzyM7OhiiKOH/+vM31uLu7m7yXyWRm94miaPM1iIiqE3uWiAg6nQ7PPvssnnrqKURFRWHs2LE4cuQIgoODnR0aEZHTsWeJiPDWW28hPz8fH3zwAaZMmYIWLVpg9OjRzg6LiKhGYLJEVMelpqZi6dKl+PLLL+Hn5we5XI4vv/wSv/32G5YvX+7s8IiInI634YjquD59+kCv15vsa9KkCfLz862uKzk5ucy+1NTUMvvsGRNFRFTdmCwRkc26d++Ojh07YteuXXbV4+PjA4PBAA8PDwdFRkTkOEyWiMhqDRs2xKlTpwDA5ikG7nT48GEAJdMKEBHVNDJJkiRnB0FERERUU3GANxEREZEFTJaIiIiILGCyRERERGQBkyUiIiIiC5gsEREREVnAZImIiIjIAiZLRERERBYwWSIiIiKy4P8BdzRxPDs2qcUAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from py_wake.ground_models import Mirror\n", "from py_wake import NOJ\n", "from py_wake.flow_map import YZGrid\n", "\n", "from py_wake.examples.data.hornsrev1 import V80, Hornsrev1Site\n", "\n", "site = Hornsrev1Site()\n", "windTurbines = V80()\n", "wt_x, wt_y = site.initial_position.T\n", "\n", "wfm = NOJ(site, windTurbines, k=.5, groundModel=Mirror())\n", "\n", "plt.figure()\n", "wfm([0], [0], wd=0).flow_map(YZGrid(x=0, y=np.arange(-300, 100, 1) + .1, z=np.arange(-100, 200))).plot_wake_map()\n", "plt.xlabel('x [m]')\n", "plt.ylabel('y [m]')" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.13" } }, "nbformat": 4, "nbformat_minor": 5 }