From 14459883ddcfd434e9975d87b55074782459fff1 Mon Sep 17 00:00:00 2001 From: aidanjgriffiths Date: Wed, 10 Jan 2024 09:44:19 +1100 Subject: [PATCH] feat: new flip flop tutorial notebook --- tutorials/Flip_flop_Index.ipynb | 1813 +++++++++++++++++++++++++++++++ 1 file changed, 1813 insertions(+) create mode 100644 tutorials/Flip_flop_Index.ipynb diff --git a/tutorials/Flip_flop_Index.ipynb b/tutorials/Flip_flop_Index.ipynb new file mode 100644 index 00000000..5422e078 --- /dev/null +++ b/tutorials/Flip_flop_Index.ipynb @@ -0,0 +1,1813 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 1, + "id": "91abca47", + "metadata": { + "ExecuteTime": { + "end_time": "2023-11-09T04:58:29.280799Z", + "start_time": "2023-11-09T04:58:29.275796Z" + } + }, + "outputs": [], + "source": [ + "from scipy.stats import skewnorm\n", + "import numpy as np\n", + "import xarray as xr\n", + "import pandas as pd\n", + "\n", + "from scores.continuous import mse\n", + "from scores.continuous import flip_flop_index_proportion_exceeding" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "id": "ac3ef136", + "metadata": { + "ExecuteTime": { + "end_time": "2023-11-09T04:51:28.537658Z", + "start_time": "2023-11-09T04:51:28.521157Z" + } + }, + "outputs": [], + "source": [ + "# Create synthetic temperature observations\n", + "obs = 40 * np.random.random((100, 100))\n", + "obs = xr.DataArray(\n", + " data=obs, \n", + " dims=[\"time\", \"x\"],\n", + " coords={\"time\": pd.date_range(\"2023-01-01\", \"2023-04-10\"), \"x\": np.arange(0, 100)}\n", + ")\n", + "\n", + "# Create forecasts for 7 lead days\n", + "fcst1 = xr.DataArray(data=[1]*7, dims=\"lead_day\", coords={\"lead_day\": np.arange(1, 8)})\n", + "fcst1 = fcst1 * obs\n", + "fcst2 = fcst1.copy()\n", + "\n", + "# add some noise to the forecasts\n", + "# fcst 2 will be set up to flip-flop more than fcst 1\n", + "fcst1 += 2 * skewnorm.rvs(4, size=(7, 100, 100))" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "id": "593844ac", + "metadata": { + "ExecuteTime": { + "end_time": "2023-11-09T04:57:02.193648Z", + "start_time": "2023-11-09T04:57:02.179379Z" + } + }, + "outputs": [], + "source": [ + "flipflop_noise = []\n", + "for lead_day in np.arange(1, 8):\n", + " noise = 2 * skewnorm.rvs(4, size=(100, 100))\n", + " if lead_day % 2 == 0:\n", + " noise *= -1\n", + " flipflop_noise.append(noise)\n", + "flipflop_noise = np.stack(flipflop_noise) \n", + "fcst2 += flipflop_noise" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "id": "fadab7f3", + "metadata": { + "ExecuteTime": { + "end_time": "2023-11-09T04:57:44.034837Z", + "start_time": "2023-11-09T04:57:44.020020Z" + } + }, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "
<xarray.DataArray (lead_day: 7)>\n",
+       "array([3.93994022, 4.05165078, 3.93569259, 3.93866259, 3.95555321,\n",
+       "       3.8961911 , 4.05679397])\n",
+       "Coordinates:\n",
+       "  * lead_day  (lead_day) int64 1 2 3 4 5 6 7
" + ], + "text/plain": [ + "\n", + "array([3.93994022, 4.05165078, 3.93569259, 3.93866259, 3.95555321,\n", + " 3.8961911 , 4.05679397])\n", + "Coordinates:\n", + " * lead_day (lead_day) int64 1 2 3 4 5 6 7" + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "mse(fcst1, obs, preserve_dims=[\"lead_day\"])" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "id": "ef993c82", + "metadata": { + "ExecuteTime": { + "end_time": "2023-11-09T04:57:45.759437Z", + "start_time": "2023-11-09T04:57:45.739361Z" + } + }, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "
<xarray.DataArray (lead_day: 7)>\n",
+       "array([4.01171069, 4.01689514, 3.98698368, 3.95889378, 4.09018116,\n",
+       "       4.05946734, 3.94731563])\n",
+       "Coordinates:\n",
+       "  * lead_day  (lead_day) int64 1 2 3 4 5 6 7
" + ], + "text/plain": [ + "\n", + "array([4.01171069, 4.01689514, 3.98698368, 3.95889378, 4.09018116,\n", + " 4.05946734, 3.94731563])\n", + "Coordinates:\n", + " * lead_day (lead_day) int64 1 2 3 4 5 6 7" + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "mse(fcst2, obs, preserve_dims=[\"lead_day\"])" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "id": "30b958f2", + "metadata": { + "ExecuteTime": { + "end_time": "2023-11-09T04:57:59.012220Z", + "start_time": "2023-11-09T04:57:58.997135Z" + } + }, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "
<xarray.DataArray ()>\n",
+       "array(3.96778349)
" + ], + "text/plain": [ + "\n", + "array(3.96778349)" + ] + }, + "execution_count": 6, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "mse(fcst1, obs)" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "id": "e5208781", + "metadata": { + "ExecuteTime": { + "end_time": "2023-11-09T04:58:03.470020Z", + "start_time": "2023-11-09T04:58:03.459742Z" + } + }, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "
<xarray.DataArray ()>\n",
+       "array(4.01020677)
" + ], + "text/plain": [ + "\n", + "array(4.01020677)" + ] + }, + "execution_count": 7, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "mse(fcst2, obs)" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "id": "20ade6c3", + "metadata": { + "ExecuteTime": { + "end_time": "2023-11-09T05:07:17.160703Z", + "start_time": "2023-11-09T05:07:16.955346Z" + } + }, + "outputs": [ + { + "data": { + "text/plain": [ + "[]" + ] + }, + "execution_count": 8, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAGwCAYAAAB7MGXBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAABIDUlEQVR4nO3deXxU5aH/8c+ZSSZ7whJICATCvhNkbUTFhYobal2KS0Wxen/Xiy2a6hVqkVqrqFWvrXKl0rp0QagbLiDKRZGqKAKCIPuOQBLCkhWyzJzfHycTiARIQpJnlu/79ZrXnJzMyXwTfTlfz/Oc51i2bduIiIiIGOIyHUBERETCm8qIiIiIGKUyIiIiIkapjIiIiIhRKiMiIiJilMqIiIiIGKUyIiIiIkZFmA5QFz6fj71795KQkIBlWabjiIiISB3Ytk1RURFpaWm4XCc//xEUZWTv3r2kp6ebjiEiIiINsHv3bjp06HDS7wdFGUlISACcXyYxMdFwGhEREamLwsJC0tPTqz/HTyYoyoh/aCYxMVFlREREJMicboqFJrCKiIiIUSojIiIiYpTKiIiIiBilMiIiIiJGqYyIiIiIUSojIiIiYpTKiIiIiBilMiIiIiJGqYyIiIiIUSojIiIiYlS9y8iSJUsYM2YMaWlpWJbF3LlzT3vM4sWLGTRoEFFRUXTr1o1XXnmlAVFFREQkFNW7jJSUlJCZmcn06dPr9Prt27dz+eWXc8EFF7Bq1Sruuece7rjjDj788MN6hxUREZHQU+8b5V166aVceumldX79jBkz6Ny5M08//TQAvXv35rPPPuN//ud/GD16dH3fvnEV5YK37MT9rgiIiAa3p+o5KO4nKCIiEpSa/FN26dKljBo1qsa+0aNHc88995z0mLKyMsrKjpWEwsLCJslmz/kZ1vfLTv9Cyw2RseCJA08sRMZBVDxEt4CYFseek9KhdTdo3RViW8Np7lIoIiIizVBGcnJySElJqbEvJSWFwsJCjhw5QkxMzAnHTJs2jYcffripo7GnqJI2eIh0W7j8xcG2wVcJtvfYC20vlBc5j7qKToLW3SG1H6QOcB4pfZ0yIyIiItUCcvxh8uTJZGdnV39dWFhIenp6o75HaXklN5RP4fujR4j1uHnyugFcMSDt2At8Xqgsg8qjznNFKZSXHHsuK4Kjh+HIYee59CAc2gEHt0HBbjhaAHuWOw8/y+UUki4XQNcLoWMWREY36u8lIiISbJq8jKSmppKbm1tjX25uLomJibWeFQGIiooiKiqqSXPFeiJ49+5z+MVrK/l8ywHunvUNq3cf5oFLehHhdoHL7ZzFaMiZjIojcHA77F8POWucx75voSTv2Ndf/MmZj9LpbOh/PfS5yhkGEhERCTNNXkaysrKYP39+jX0LFy4kKyurqd/6tFrFeXh1/DCe+mgTMz7dysx/b2fNngKev2kQyfFnUIYiYyClj/Pod+2x/YX7YMdnsPVj2PYJFO1ztrd+DPP/G/pfC2eNg/aDNN9ERETChmXbtl2fA4qLi9myZQsAZ511Fs888wwXXHABrVq1omPHjkyePJk9e/bwt7/9DXAu7e3Xrx8TJkzg9ttv5+OPP+aXv/wl8+bNq/PVNIWFhSQlJVFQUEBiYmI9f8W6+WDNPu57fTUl5V56piQw75fnOGdImoptw/6NsOF9+OYfcGj7se+1y4TR0yBjRNO9v4iISBOr6+d3vcvI4sWLueCCC07Yf+utt/LKK69w2223sWPHDhYvXlzjmHvvvZd169bRoUMHpkyZwm233dbov8yZ2pJXxPUzlnKotIInrx3AT4c27jyVk/L5YOfn8M3fYd07zjwVgH7XwcWPQGLaqY8XEREJQE1WRkxorjICMHPJNh6dv560pGg+uf98oiLcTfp+Jyg5AJ/8Hpa/DNjOZcTn3QdZEyCiaefRiIiINKa6fn7r3jQ/cEtWJ1ITo9lbcJR/frmr+QPEtYYr/gf+YzF0GAYVJbDoYZh5IeRtaP48IiIiTUxl5AeiI91MHNUdgOmfbKG4rNJMkLSBcPuHcPUMZwG13LXw4kj4+i/OfBMREZEQoTJSi+sGdyCjdSwHSsp56bPtpz+gqbhcMPBGuGspdL3ImUsy71cw+yYoyTeXS0REpBGpjNQi0u0i++KegDOH5FBJudlACSlw8xvOFTZuD2ycDy+cDbnfmc0lIiLSCFRGTuKK/u3o0y6RorJKZny61XQc5yxJ1n/BnR9Dck8ozoW/XQ0HAiCbiIjIGVAZOQmXy+L+0c7ZkVe+2EFOwVHDiaqk9oeffwgp/Z0VXf92FRR8bzqViIhIg6mMnML5PdswNKMlZZU+Xv7c4NyRH4ppCbe85dwhuGC3U0iK80ynEhERaRCVkVOwLIufn9MZgHdX78XnC6CrWOLbwrh3ICkdDmyBv18DRw6ZTiUiIlJvKiOncX7PtiRERbCv4Chf7zhoOk5NSR2cQhKfArlrYNZYqAiQ4SQREZE6Uhk5jehIN5f0SwXgndV7DaepReuucMtciE6C3V/B+/doHRIREQkqKiN1cNXA9gDMX7OP8kqf4TS1SOkD178ClhtWvwZfPGc6kYiISJ2pjNRBVtfWtEmI4nBpBUs27Tcdp3ZdL4RLpjnbCx+CTR+ZzSMiIlJHKiN14HZZjBng3Dk3IIdq/Ib9Bwy6FbDhjdt1LxsREQkKKiN1dNVAp4wsXJdDian71ZyOZcFlT0GnEVBeBK/dAKUBNulWRETkB1RG6mhAhyQyWsdytMLHwnW5puOcXIQHfvo3aNERDm13JrSKiIgEMJWROrIsiyurJrK+s2qP4TSnEZcMP/27M6F13Tvw3dumE4mIiJyUykg9XJnpDNUs2ZzPgeIyw2lOI20gnJvtbM+7T3f5FRGRgKUyUg/d2sbTr30iXp/N/DX7TMc5vfPuhza9oTQfPvhv02lERERqpTJST1dl+odqAviqGr+IKLj6f53hmrVvwvr3TCcSERE5gcpIPY3JTMOyYPnOQ+QVBsHS6+0HwYhfOtvvZ+vqGhERCTgqI/WUmhRNr9REAFbsDJIb042cBMk9oSQPFkwynUZERKQGlZEGGNSxBQArdwVJGYmMhqumg+WCb+fA1k9MJxIREammMtIAgzu1BILozAhA+lAYeqezPf9+qCw3m0dERKSKykgDDOrolJG1ewopq/QaTlMPF/wa4trAgc3w5f+aTiMiIgKojDRIp9axtI7zUO71sXZPoek4dRfTAn78O2f70yehIMAXbxMRkbCgMtIAlmVxVtXZkZXBNFQDMOAGSB8OFSXw0YOm04iIiKiMNJR/3kjQTGL1c7mcm+lZLmeZ+G2LTScSEZEwpzLSQP4ysnznIWzbNpymntoNgCE/d7bn/7cms4qIiFEqIw00oEMSES6L/UVlfH/oiOk49XfhgxCbDPkb4asZptOIiEgYUxlpoOhIN33TnMXPgm6oBiCmJYya6mx/9j9QVmw2j4iIhC2VkTMQtJNY/TJvglZd4MhBWP5X02lERCRMqYycgWOTWA+bDdJQ7gg491fO9hfPQXmp2TwiIhKWVEbOgL+MrNtXSGl5peE0DTRgLLToCCX7YcUrptOIiEgYUhk5A2ktYkhNjMbrs/n2+wLTcRrGHQnnZDvbn/8RKoLgTsQiIhJSVEbOUFDep+aHBt4Eie2hOAe++bvpNCIiEmZURs7QWf47+AZzGYmIghH3ONufPat1R0REpFmpjJyh41diDbrFz4436BaIT4HC72H1a6bTiIhIGFEZOUN905LwRLg4VFrB9vwS03EaLjIGRkx0tj97BrxBOiFXRESCjsrIGfJEuBjQPgkI4kt8/QaPd1ZlPbQD1vzLdBoREQkTKiONICQmsQJ4YuHsu53tfz8NPq/ZPCIiEhZURhqBfxLrt98fNpqjUQy9w1kq/sAW566+IiIiTUxlpBH0TXOGaTblFlFWGeRnE6ISIGuCs/3pk+Dzmc0jIiIhT2WkEXRoGUOL2EgqvDabckLghnPD/gOik5w7+q5/x3QaEREJcSojjcCyLPpVnR1ZuzdIV2I9XnQSDL/L2f70Dzo7IiIiTUplpJH0bZ8IwNo9IVBGAH70n+BJgLzvYON802lERCSEqYw0kmNnRgoNJ2kkMS1h+P9ztj99AoJ5QTcREQloKiONpF/VWiPr9xVS4Q2RYY0f/RdExkHOt7DpQ9NpREQkRKmMNJJOrWJJiIqgvNLHlrwQmMQKENcaht3hbC95UmdHRESkSaiMNBKXy6JPWojNGwHI+gW4o2DPCti9zHQaEREJQSojjcg/VPNdqMwbAYhvAwOud7a/esFsFhERCUkqI42oX6hdUePnv8x33btQ8L3ZLCIiEnJURhpR/+POjHh9ITS/IrUfZJwLthe+/ovpNCIiEmJURhpR5+R4YiLdHKnwsj0/RCax+g3/T+d5xStQXmo0ioiIhBaVkUbkrjGJNYTmjQD0vBRadIIjh2DN66bTiIhICFEZaWT9QvGKGgCX27lnDcBXM3SZr4iINBqVkUbWt2reyJpQKyMAZ/3MWQQtbx1sX2I6jYiIhAiVkUbmn8S6bm8hvlCaxAoQ0wIG3uhsfzXDaBQREQkdDSoj06dPJyMjg+joaIYPH86yZadeDOvZZ5+lZ8+exMTEkJ6ezr333svRo0cbFDjQdWsbjyfCRVFZJbsOhuBET/9E1o0fwMFtZrOIiEhIqHcZmTNnDtnZ2UydOpWVK1eSmZnJ6NGjycvLq/X1s2bNYtKkSUydOpX169fz17/+lTlz5vDrX//6jMMHoki3i96pCQCs3RuCQzXJ3aHbKMCGr/9qOo2IiISAepeRZ555hjvvvJPx48fTp08fZsyYQWxsLC+99FKtr//iiy8YMWIEN910ExkZGVx88cXceOONpzybUlZWRmFhYY1HMAnpeSMAw6ru5vvNP6DiiNksIiIS9OpVRsrLy1mxYgWjRo069gNcLkaNGsXSpUtrPebss89mxYoV1eVj27ZtzJ8/n8suu+yk7zNt2jSSkpKqH+np6fWJaVy/tKrFz0Lt8l6/bhdBUkc4ehi+e9t0GhERCXL1KiP5+fl4vV5SUlJq7E9JSSEnJ6fWY2666SZ+97vfcc455xAZGUnXrl05//zzTzlMM3nyZAoKCqofu3fvrk9M4/yTWNfuLcAOxUtgXW4Ycpuzvbz2M2IiIiJ11eRX0yxevJjHHnuM//3f/2XlypW89dZbzJs3j0ceeeSkx0RFRZGYmFjjEUx6pMYT4bI4XFrBnsMhOoxx1i3gioTvv4Z935pOIyIiQaxeZSQ5ORm3201ubm6N/bm5uaSmptZ6zJQpU7jlllu444476N+/Pz/5yU947LHHmDZtGj6fr+HJA1hUhJseKc4k1pC6g+/x4ttC7zHOts6OiIjIGahXGfF4PAwePJhFixZV7/P5fCxatIisrKxajyktLcXlqvk2brcbIDSHMKr0SIkHYEteiN2j5nhDf+48f/svOBqipUtERJpcvYdpsrOzmTlzJq+++irr16/nrrvuoqSkhPHjxwMwbtw4Jk+eXP36MWPG8MILLzB79my2b9/OwoULmTJlCmPGjKkuJaGoW1unjGwN5TLSaQQk94SKEvh2juk0IiISpCLqe8DYsWPZv38/Dz30EDk5OQwcOJAFCxZUT2rdtWtXjTMhv/nNb7Asi9/85jfs2bOHNm3aMGbMGB599NHG+y0CkL+MbNkfwmXEsmDI7bDgAVj+Mgy9w9knIiJSD5YdBGMlhYWFJCUlUVBQEDSTWTfnFvHj/1lCnMfN2odHY4Xqh/SRw/BMb6gohds/hI4/Mp1IREQCRF0/v3VvmibSqXUcbpdFSbmXfQWhufQ94Nyvpt+1zrZWZBURkQZQGWkinggXnVrHAiE+iRWcoRqAdXOh9KDRKCIiEnxURppQtzZVk1hDed4IQPtBkDoAvOWw9k3TaUREJMiojDSh6kmsoX5mBCDzBudZV9WIiEg9qYw0obAqI/2uA8vlrMh6YKvpNCIiEkRURppQ13AZpgFISIGuFzrbOjsiIiL1oDLShLpWnRnJLy7ncGm54TTNYMBxQzWBf8W4iIgECJWRJhQfFUG7pGggTM6O9LocPPFwaAfs/sp0GhERCRIqI00srOaNeGKh95XO9urZZrOIiEjQUBlpYv55I2FRRgAyxzrP370NlWVms4iISFBQGWliXcPpzAhAxrmQkAZHD8OmD02nERGRIKAy0sT8C5+F9A3zjudyw4DrnW1dVSMiInWgMtLE/HNGvj90hKMVXsNpmon/qppNH2p5eBEROS2VkSaWHO8hKSYS24Zt+0tMx2keKX0gtT/4KuC7t0ynERGRAKcy0sQsyzp2RU24DNXAcWuO/MtsDhERCXgqI82ga5s4IIwmsQL0uwawnPVGDu82nUZERAKYykgz8J8ZCYuFz/wS06DT2c72d2+bzSIiIgFNZaQZVJeRcDozAlVnR9C8EREROSWVkWbQrU0CANvyS/D6wuieLb2vcu7ku/cb3clXREROSmWkGbRvGUNUhIvySh+7D5aajtN84ttA55HOtoZqRETkJFRGmoHbZdE5OQwnsQL0u9Z5XquhGhERqZ3KSDMJy0msAL2vAFck5H0HeRtMpxERkQCkMtJMwuruvceLaQndLnK2NZFVRERqoTLSTMJy4TO/6qGaN8EOowm8IiJSJyojzeT4MyN2uH0g97wUIqLhwBbIWWM6jYiIBBiVkWaS0ToOy4Kio5XkF5ebjtO8ohKg+8XO9to3zWYREZGAozLSTKIj3aS3jAXCcBIr1FwALdzODImIyCmpjDQj/z1qwrKMdB8NkXFweBd8v9x0GhERCSAqI82oaxv/svAlhpMY4ImFXpc726tfM5tFREQCispIM+oarmuN+A28yXle+wZUHDWbRUREAobKSDOqPjMSrmWk80hISoejBbDhfdNpREQkQKiMNCP/nJE9h49wpNxrOI0BLhdk3uhsr5plNouIiAQMlZFm1CrOQ4vYSGwbtueH4bwRgIFVZWTrx1Cwx2wWEREJCCojzciyLA3VtOoCnUYANnw723QaEREJACojzSysL+/1809k/eafWnNERERURprbsTMjYTpMA9DnamfNkYNbYfdXptOIiIhhKiPN7NhaI2F8ZiQqHvpc5Wyv+qfZLCIiYpzKSDPzrzWyLb8Yny+MhyjOutl5Xvs2lIfxWSIREVEZaW7pLWOIdFscrfCxt+CI6TjmdDwbWmZAeRGs15ojIiLhTGWkmUW4XWS09k9iDeMzAi4XZFZNZF31D7NZRETEKJURAzRvpErmDc7z9n9DUa7ZLCIiYozKiAFd2+ryXgBadoK0QYCt5eFFRMKYyogBYb/w2fH6XOk8r3/XbA4RETFGZcQArTVynN5VZWT7v6H0oNksIiJihMqIAV2qVmHdX1RGwZEKw2kMa90VUvqB7YUN80ynERERA1RGDEiIjiQlMQqAbRqqObYAmoZqRETCksqIIRqqOY5/qGbrJ3C0wGwWERFpdiojhmgS63Ha9oLkHuCrgE0fmk4jIiLNTGXEkOq794b7WiN+/rMj694xm0NERJqdyoghx+5Ro2Ea4NglvlsW6V41IiJhRmXEEP8wzc4DJVR4fYbTBIDUAdCiE1Qegc0LTacREZFmpDJiSGpiNLEeNxVem90HS03HMc+ytACaiEiYUhkxxOWyqtcb0RU1VXpXXeK76UOoOGo2i4iINBuVEYP8QzVbNInV0X4wJKRBeTFs+8R0GhERaSYqIwZ1Sa6axKrLex0uF/Qe42yv01CNiEi4UBkxSHfvrYW/jGz6ALxhvlS+iEiYaFAZmT59OhkZGURHRzN8+HCWLVt2ytcfPnyYCRMm0K5dO6KioujRowfz589vUOBQcvwqrLZtG04TIDpmQWxrOHIIdn5uOo2IiDSDepeROXPmkJ2dzdSpU1m5ciWZmZmMHj2avLy8Wl9fXl7Oj3/8Y3bs2MEbb7zBxo0bmTlzJu3btz/j8MGuc3IclgUFRyo4UFJuOk5gcEdAz8uc7fXvmc0iIiLNot5l5JlnnuHOO+9k/Pjx9OnThxkzZhAbG8tLL71U6+tfeuklDh48yNy5cxkxYgQZGRmMHDmSzMzMk75HWVkZhYWFNR6hKDrSTfsWMYBWYq3Bvxrrhnng0xosIiKhrl5lpLy8nBUrVjBq1KhjP8DlYtSoUSxdurTWY959912ysrKYMGECKSkp9OvXj8ceewyv13vS95k2bRpJSUnVj/T09PrEDCr+oRqtxHqcLiPBkwBF+2DPCtNpRESkidWrjOTn5+P1eklJSamxPyUlhZycnFqP2bZtG2+88QZer5f58+czZcoUnn76aX7/+9+f9H0mT55MQUFB9WP37t31iRlUqueN6MzIMRFR0ONiZ1sLoImIhLwmv5rG5/PRtm1bXnzxRQYPHszYsWN58MEHmTFjxkmPiYqKIjExscYjVOmKmpPwX1Wz4X3Q5F4RkZAWUZ8XJycn43a7yc3NrbE/NzeX1NTUWo9p164dkZGRuN3u6n29e/cmJyeH8vJyPB5PA2KHjuOvqJHjdPsxuKPg4DbIWwcpfU0nEhGRJlKvMyMej4fBgwezaNGi6n0+n49FixaRlZVV6zEjRoxgy5Yt+I6biLhp0ybatWsX9kUEqF4SfvehUo5WnHweTdiJioeuFzrbuqpGRCSk1XuYJjs7m5kzZ/Lqq6+yfv167rrrLkpKShg/fjwA48aNY/LkydWvv+uuuzh48CATJ05k06ZNzJs3j8cee4wJEyY03m8RxNrER5EQHYFtw84DumFeDf6hGpUREZGQVq9hGoCxY8eyf/9+HnroIXJychg4cCALFiyontS6a9cuXK5jHSc9PZ0PP/yQe++9lwEDBtC+fXsmTpzIAw880Hi/RRCzLIuubeJZtfswW/cX0zM1wXSkwNHzUrDckLvWGa5p1cV0IhERaQKWHQRLfxYWFpKUlERBQUFITmb91b9W8+bK7/nVj3vwi4u6m44TWF4dA9uXwI8fgRG/NJ1GRETqoa6f37o3TQDQFTWn4F8ATUM1IiIhS2UkAPjv3qsramrR63Ln+ftlULjPbBYREWkSKiMBoFvVmZFt+4t1w7wfSkyDDkOd7Q3vm80iIiJNQmUkAHRsFYfbZVFS7iW3sMx0nMDT5yrned07ZnOIiEiTUBkJAJ4IF51axQKaN1Ir/yW+Oz+H4v1ms4iISKNTGQkQXapXYlUZOUHLDGg3EGyfhmpEREKQykiA6Fq1EqtumHcSffxX1ejGeSIioUZlJED471GzLV9X1NSqd9W8ke1LoPSg2SwiItKoVEYCRPVaIzozUrvkbtC2L/gqYeMHptOIiEgjUhkJEP61RvYWHKWkrNJwmgDlv6pGQzUiIiFFZSRAtIzz0DrOuYvxdg3V1M4/b2Trx3C00GwWERFpNCojAaRLGy0Lf0ptekFyD/CWw6YPTacREZFGojISQLq20bLwp2RZx+5Vs26u0SgiItJ4VEYCSFetNXJ6/nkjW/4PyvR3EhEJBSojAURX1NRBan9nEbTKo7Bloek0IiLSCFRGAoj/iprt+SV4fbphXq0sS/eqEREJMSojAaRDyxg8bhdllT72HDpiOk7g8s8b2bwQKnVjQRGRYKcyEkAi3K7qK2o25RYZThPA0gZBQjsoL4Ztn5pOIyIiZ0hlJMD0SEkAYFOeyshJuVzQ63Jne8N7ZrOIiMgZUxkJMD1SnHkjm3M1ifWU/GVk4wfg85rNIiIiZ0RlJMB0958Z0TDNqWWcC9FJULIfdi8znUZERM6AykiA8Q/TbMkr1hU1p+KOhB6XONsb3jebRUREzojKSIDp2CqWqAjniprdB0tNxwlsva5wnte/B7aKm4hIsFIZCTBul0W3ts68kY0aqjm1bhdBRDQc3gm535lOIyIiDaQyEoD8QzWbVUZOzRMHXS90tjVUIyIStFRGAlD3qitqNumKmtOrHqpRGRERCVYqIwGoR1tdUVNnPS8Fyw25a+DQDtNpRESkAVRGAlDPVKeMbNtfQqXXZzhNgIttBZ3OdrY3zDObRUREGkRlJAC1bxFDTKSbcq+Pnbqi5vQ0VCMiEtRURgKQy2UdmzeSo6Ga0/KvxrprKRTvN5tFRETqTWUkQHWvnjeiSayn1SId2g0EbNiooRoRkWCjMhKg/Peo0Q3z6qjPlc7zunfM5hARkXpTGQlQPVK11ki99Lnaed72KZQeNBpFRETqR2UkQPkXPtueX0KFrqg5vdZdIaU/2F7YON90GhERqQeVkQCVlhRNfFQEFV6bHfklpuMEhz5XOc/fzTUaQ0RE6kdlJEBZ1rF71GgSax35y8i2xXDkkNEoIiJSdyojAcw/iVU3zKujNj2gbR/wVcDGBabTiIhIHamMBDDdMK8B/GdHdFWNiEjQUBkJYP4yonvU1IO/jGxdBEcLzGYREZE6URkJYP4ysuNAKWWVXsNpgkSbXpDcA7zlsOlD02lERKQOVEYCWEpiFAnREXh9Ntt1RU3dWNaxNUc0VCMiEhRURgKYZVnHDdXoipo68w/VbF4IZRriEhEJdCojAa6HbphXfyl9oVVX8JZpqEZEJAiojAQ4TWJtAMvSVTUiIkFEZSTA9VQZaZi+VzvPmxdCuebbiIgEMpWRAOe/Yd7Og6UcKdcVNXWWOgBadILKI04hERGRgKUyEuCS46NIjvdg27A5T2dH6syyoM+Vzvb6d81mERGRU1IZCQL+eSMbNYm1fvyX+G76ECqOGo0iIiInpzISBHqmqow0SNogSGwP5cWw9WPTaURE5CRURoKAfxKrbphXTy4X9B7jbGuoRkQkYKmMBIEeOjPScL2r5o1snA+V5WaziIhIrVRGgoB/zkheURmHSvSBWi8dfwRxbZ2b5u1YYjqNiIjUQmUkCMRHRdChZQygoZp6c7mh9xXOthZAExEJSCojQaJXqhY/azD/UM2GeeCtNJtFREROoDISJPxDNRs0b6T+Ms6BmJZQegB2fWE6jYiI/IDKSJDwX96rG+Y1gDsSel3ubGuoRkQk4DSojEyfPp2MjAyio6MZPnw4y5Ytq9Nxs2fPxrIsrr766oa8bVirXmsktwjbtg2nCUK9q26ct/598PnMZhERkRrqXUbmzJlDdnY2U6dOZeXKlWRmZjJ69Gjy8vJOedyOHTu47777OPfccxscNpx1SY4nwmVRdLSSfQVaTbTeuoyEqEQozoHv61aeRUSkedS7jDzzzDPceeedjB8/nj59+jBjxgxiY2N56aWXTnqM1+vl5ptv5uGHH6ZLly5nFDhceSJcdGkTB2i9kQaJiIIelzjbGqoREQko9Soj5eXlrFixglGjRh37AS4Xo0aNYunSpSc97ne/+x1t27bl5z//eZ3ep6ysjMLCwhoPgZ6piYAu722wPlVDNeve0VCNiEgAqVcZyc/Px+v1kpKSUmN/SkoKOTk5tR7z2Wef8de//pWZM2fW+X2mTZtGUlJS9SM9Pb0+MUNWz5R4QGdGGqzbReCJh8I9sGeF6TQiIlKlSa+mKSoq4pZbbmHmzJkkJyfX+bjJkydTUFBQ/di9e3cTpgwe1WdGVEYaJjLmuKGauUajiIjIMRH1eXFycjJut5vc3Nwa+3Nzc0lNTT3h9Vu3bmXHjh2MGTOmep+v6vR4REQEGzdupGvXriccFxUVRVRUVH2ihQX/DfO27C+m0usjwq0rs+utz1Ww9g1nqObi34NlmU4kIhL26vVp5vF4GDx4MIsWLare5/P5WLRoEVlZWSe8vlevXqxZs4ZVq1ZVP6688kouuOACVq1apeGXeurQMoZYj5vySh87DpSajhOcuv8YIuOgYDfsWWk6jYiIUM8zIwDZ2dnceuutDBkyhGHDhvHss89SUlLC+PHjARg3bhzt27dn2rRpREdH069fvxrHt2jRAuCE/XJ6LpdF95QEVu8+zKbcIrq1jTcdKfhExkCP0fDdW85QTYfBphOJiIS9ep/nHzt2LE899RQPPfQQAwcOZNWqVSxYsKB6UuuuXbvYt29fowcVh38Sq5aFPwPVV9XMBS0gJyJinGUHwXKehYWFJCUlUVBQQGJiouk4Rv31s+088v46Lumbyoxb9H/1DVJeCn/oChWlcOcn0H6Q6UQiIiGprp/fmgEZZPyTWLXWyBnwxEL3i51tLYAmImKcykiQ8d+jZseBEo5WeA2nCWIaqhERCRgqI0EmOd5DqzgPtg2bc4tNxwlePUZDRAwc2gH7VptOIyIS1lRGgoxlWfSqOjuybl+B4TRBzBPnXOYLGqoRETFMZSQI9W+fBMC336uMnBEN1YiIBASVkSDUv4NTRtbsURk5Iz0ugYhoOLgN9moBNBERU1RGglBmhxYArN9XSFmlJrE2WFQ89K66VcHKv5vNIiISxlRGglCHljG0iI2kwmuzKUeTWM/IWbc4z2vfdNYfERGRZqcyEoQsyzo2b2TPYbNhgl3GudCiE5QVaiKriIghKiNBakDVvJFvd2veyBlxuY6dHflGQzUiIiaojASp/u1bAPCtJrGeuYE3geWCnZ/Dga2m04iIhB2VkSCVme6cGdmUW6SVWM9UUnvoepGzrbMjIiLNTmUkSKUmRpMcH4XXZ7NuX6HpOMFvUNVQzarXwFtpNouISJhRGQlSlmVVzxtZo8XPzlyPSyE2GYpzYMtC02lERMKKykgQ819Rs/r7w2aDhIIID2Te4GxrzRERkWalMhLEdGakkfmvqtm0AIpyzWYREQkjKiNBzL8s/Jb9xZSUaZ7DGWvbCzoMBdsLq18znUZEJGyojASxtgnRtEuKxrbhu72axNoo/GdHVv1TN88TEWkmKiNB7tgdfA+bDRIq+v4EImIgfxPs0c3zRESag8pIkKteiVXzRhpHdCL0vsLZXj3LbBYRkTChMhLk+lfdwXeNVmJtPJk3Os9r34TKMrNZRETCgMpIkBtQNUyzPb+EgiMVhtOEiC7nQ0I7OHIINn1oOo2ISMhTGQlyLeM8pLeKAeA7nR1pHC43DPips62rakREmpzKSAgYoJvmNb7Mm5znzR9BSb7ZLCIiIU5lJAT01+Jnja9tL0g7C3yVsOYN02lEREKaykgI8M8bWbX7sNkgocZ/dkRX1YiINCmVkRDQv0MSLgv2HD5CXuFR03FCR79rwRUJ+1ZD7jrTaUREQpbKSAhIiI6kZ2oiACt2HjKcJoTEtYYeo51tTWQVEWkyKiMhYnCnFgAsVxlpXP41R779F3h1/x8RkaagMhIihnRqBaiMNLruF0NMKyjOgW2fmE4jIhKSVEZCxOBOLQFnrZGjFV7DaUJIhOfYmiMrXzWbRUQkRKmMhIgOLWNomxBFpc9mta6qaVyDbnWeN34AxXlms4iIhCCVkRBhWRZDMpyzIyt2aaimUaX0gQ5DnTVHVukyXxGRxqYyEkIGdawqIztURhrdoHHO88q/gW2bzSIiEmJURkLIkAxnEuuKXYfw+fSB2aj6XgOeeDi4FXZ8ZjqNiEhIURkJIX3aJRIV4eJwaQXb8ktMxwktUfHQ/zpne+XfzGYREQkxKiMhxBPhIjO9BQArdh40GyYU+SeyrnsHSvX3FRFpLCojIWZI1SW+yzVvpPGlnQUp/cFb5iyCJiIijUJlJMT41xvRFTVNwLJgcNXZkZWvaiKriEgjURkJMf4rarbtL+FgSbnhNCGo//UQEQN562DPCtNpRERCgspIiGkZ56FrmzgAVmpp+MYX0wL6Xu1sr3jFYBARkdChMhKCdJ+aJuZfc2Ttm3C0wGwWEZEQoDISgqrnjeiKmqbRMQva9IaKUlg923QaEZGgpzISggZXLQu/+vsCyit9htOEIMuCoT93tr/+iyayioicIZWRENQlOY6WsZGUV/pYu1fDCE0i8wZnRdb8TbB9iek0IiJBTWUkBFmWVT1Us3yHhmqaRFSCU0gAvp5pNouISJBTGQlRwzo7k1i/2HrAcJIQNvQO53nDfCjYYzaLiEgQUxkJUed2bwPAl9sOUFbpNZwmRLXtDZ3OAdury3xFRM6AykiI6pWaQHJ8FEcrfKzQ0vBNxz+RdeWrUKlF5kREGkJlJERZlsV53ZMB+PeWfMNpQljvMRCfAsW5sOE902lERIKSykgIO7dHVRnZvN9wkhDmjoTBtznby/5iNIqISLBSGQlhI7o5ZWTtnkIOFJcZThPCBt8Glht2fQG535lOIyISdFRGQljbhGh6t0sE4DMN1TSdxDTodbmzvXS62SwiIkFIZSTEVc8b2awy0qRGTHSeV8+GA1vNZhERCTIqIyHunO7H5o3YWra86XQYAt0vdi7zXfIH02lERIKKykiIG5rRiqgIF7mFZWzJKzYdJ7SdP8l5/nYO5G8xm0VEJIg0qIxMnz6djIwMoqOjGT58OMuWLTvpa2fOnMm5555Ly5YtadmyJaNGjTrl66VxRUe6q1djXaKhmqbVfjD0uARsH3z6hOk0IiJBo95lZM6cOWRnZzN16lRWrlxJZmYmo0ePJi8vr9bXL168mBtvvJFPPvmEpUuXkp6ezsUXX8yePVo+u7mcV7Uaqy7xbQbnT3ae174B+zeazSIiEiQsu54TCYYPH87QoUN5/vnnAfD5fKSnp/OLX/yCSZMmnfZ4r9dLy5Ytef755xk3blyd3rOwsJCkpCQKCgpITEysT1wBNuQUcsmz/yY60sXqqRcTFeE2HSm0vXYTbJwH/a6F614ynUZExJi6fn7X68xIeXk5K1asYNSoUcd+gMvFqFGjWLp0aZ1+RmlpKRUVFbRq1eqkrykrK6OwsLDGQxquZ0oCbRK0NHyz8c8dWfsW5K03m0VEJAjUq4zk5+fj9XpJSUmpsT8lJYWcnJw6/YwHHniAtLS0GoXmh6ZNm0ZSUlL1Iz09vT4x5Qcsy+JcLQ3ffNoNcJaJx9bcERGROmjWq2kef/xxZs+ezdtvv010dPRJXzd58mQKCgqqH7t3727GlKHp3O5aGr5Zjaw6O/Ld25C7zmwWEZEAV68ykpycjNvtJjc3t8b+3NxcUlNTT3nsU089xeOPP85HH33EgAEDTvnaqKgoEhMTazzkzBy/NHxe0VHDacJAaj/oc5WzveRJs1lERAJcvcqIx+Nh8ODBLFq0qHqfz+dj0aJFZGVlnfS4J598kkceeYQFCxYwZMiQhqeVBmubEE1megsAPlxbtyE1OUMjH3Cev5sLeRuMRhERCWT1HqbJzs5m5syZvPrqq6xfv5677rqLkpISxo8fD8C4ceOYPHly9eufeOIJpkyZwksvvURGRgY5OTnk5ORQXKwFuJrbmAHtAHhv9T7DScJESt9jc0d0dkRE5KTqXUbGjh3LU089xUMPPcTAgQNZtWoVCxYsqJ7UumvXLvbtO/Zh98ILL1BeXs51111Hu3btqh9PPfVU4/0WUieX9XfKyNc7D5JToKGaZuE/O7L2La07IiJyEvVeZ8QErTPSeK574QuW7zzElCv68PNzOpuOEx5m3wwb3of+18O1fzGdRkSk2TTJOiMS/K6oGqp5/9u9hpOEkZH/7TyvfRPyN5vNIiISgFRGwsxlA9rhsuCbXYfZfbDUdJzw0C4Tel7m3LNGd/QVETmBykiYaZsQzfDOrQGYt0YTWZuN/+zImtd1R18RkR9QGQlDV2RqqKbZpZ117I6+i6eZTiMiElBURsLQpf3a4XZZrN1TyI78EtNxwsf5kwHLuaPv9iWm04iIBAyVkTDUKs5TvSKrzo40o7SBMPTnzvb72VBZZjSOiEigUBkJU8euqtG8kWZ14RSIawsHNsPnfzKdRkQkIKiMhKnRfVKJdFtsyClic26R6TjhI6YFXFI1Z2TJH+DAVqNxREQCgcpImEqKjeS87m0AeE9nR5pXv2uhy/ngLYP590HgrzsoItKkVEbC2JUD0wB4c8X3eH36QGw2lgWXPQ1uD2z9GL5723QiERGjVEbC2Oi+qbSMjWTP4SMsWp9rOk54Se4G52Q72wsmw9ECs3lERAxSGQlj0ZFuxg7tCMDflu40nCYMnXMvtOoCxTnw0W9MpxERMUZlJMzdPLwjLgs+25LPlrxi03HCS2Q0XPmcs73yb7B5odk8IiKGqIyEufRWsVzUOwWAf3ypsyPNLuMcGH6Xs/3uL+DIIbN5REQMUBkRbs3KAOCNFd9TXFZpNkw4uughaNUVivY580dERMKMyogwoltrurSJo7iskrdXfm86TvjxxMLVL4DlgtWvwYb5phOJiDQrlRHBsqzqsyOvLt2JrXUvml/H4ZB1t7P93kQoPWg2j4hIM1IZEQCuGdSeOI+bLXnFLN16wHSc8HTBg9CmF5TkwbxfaTE0EQkbKiMCQEJ0JNcM6gDAq0t3mA0TriKjq4Zr3PDdW/DNP0wnEhFpFiojUm1cVicAFq7LZc/hI4bThKn2g+DCqjVH5t8PeRvM5hERaQYqI1Kte0oCZ3dtjc+GVz7fbjpO+BpxD3S9ECqPwOu3QXmp6UQiIk1KZURquPPcLgC8tmw3hUcrDKcJUy4X/OTPEJ8C+9fDgkmmE4mINCmVEanh/J5t6N42nuKySl77apfpOOErvi1c8yJgwcpXYc0bphOJiDQZlRGpwbIs7jzPOTvy8uc7KK/0GU4UxrqcD+fd52y/NxEObDUaR0SkqaiMyAmuGphG24QocgqP8t7qvabjhLeRk6Dj2VBeDK/dqLv7ikhIUhmRE0RFuLltRAYAM/+9TYugmeSOgOtegoQ0yN8Ir48Hr5bsF5HQojIitbp5WCdiPW425BSxZHO+6TjhLbEd3PgaRMbC1kXwoe5fIyKhRWVEapUUG8kNQzsCMHPJNsNphLSBzhU2AMtehGUzjcYREWlMKiNyUrefk4HbZfHZlnzW7tFcBeP6XOnc4Rfggwdg68dm84iINBKVETmpDi1jubx/O8CZOyIB4JxsyLwRbC/861bY+YXpRCIiZ0xlRE7pP6ou83139V5W7DxkOI1gWTDmj9DpHCgrhL9dDd/NNZ1KROSMqIzIKfVrn8R1gztg2zD5rW+17kggiIiCn70Bva4Ab5mzZPxXfzadSkSkwVRG5LQevKw3reM8bMotZsanWngrIETGwE//BkN+DtjwwX/DwqngU1kUkeCjMiKn1TLOw0Nj+gDw/Mdb2JJXbDiRAOByw+VPw4VTnK8/fxbe+S/w6p5CIhJcVEakTq7MTOP8nm0o9/r49Vtr8Pm0EFpAsCxnyfirXwDLDatfgzm3QMUR08lEROpMZUTqxLIsfn91P2I9bpbtOMjsr3ebjiTHG3gT3DALIqJh0wfwj2u1dLyIBA2VEamzDi1j+dXFPQGY9sF6cguPGk4kNfS8BH72FkQlws7P4ZUroHi/6VQiIqelMiL1ctvZGWR2SKLoaCW/mPUNZZVe05HkeBkj4Lb3ITYZcr6Fly+BvA2mU4mInJLKiNSL22Xx9E8zSYiKYNmOgzz49lrdSC/QtMuE2z+EpHQ4sAVeHOlc+qt/TiISoFRGpN66tU3guZvOwmXBGyu+58+6d03gSe4Gd/wfdL0IKo86l/7+41oo3Gc6mYjICVRGpEHO79mWh65wLvd9YsEGPvwux3AiOUFCKvzsTbj0D87E1q2L4IUsrdgqIgFHZUQa7NazM/jZjzpi23DP7FW6mV4gsiwY/h/wH59C6gA4cghevxVm3wyFe02nExEBVEbkDFiWxdQxfTmnWzJHKrzc8epytuQVmY4ltWnbC+5YBOfdD64I2PA+PD8Mls3Uqq0iYpzKiJyRSLeL6TcPolvbeHIKj3LtC0tZtv2g6VhSmwgPXPgb+H9LoP0QKC+C+ffBS6Nh37em04lIGFMZkTOWFBPJv/5fFmd1bEHBkQp+9tevmPetJkoGrJS+8POPnLkknnj4fhn8+Tx4+y4o2GM6nYiEIZURaRSt4jzMuuNHXNwnhfJKH3e/tpK/frbddCw5GZfbmUsy4Svody1gw+pZ8NxgWPQIlGm4TUSaj2UHwSIRhYWFJCUlUVBQQGJiouk4cgpen81v3/2Ov3+5E4CbhnfkoSv6EB3pNpxMTun75fDRb2DXUufr2NYw7P/BsDshtpXZbCIStOr6+a0yIo3Otm3+vGQbj3/grPzZvW08f7rxLHq30z+7gGbbsGEeLHwIDm519kXGwlm3QNYEaNnJbD4RCToqI2Lckk37+dXrq9lfVIbH7WLSpb0YPyIDy7JMR5NT8VbCurnw+R+dJeXBuSNwn6sg627oMNhoPBEJHiojEhAOFJfx3298y6INeQCM7NGGx67pT/sWMYaTyWnZNmz/1CklWz8+tj/9R86Zkl6XO3NPREROQmVEAoZt2/zjy538ft56yip9REe6uPuCbtx5XheiIvRhFhRy1sDS/4U1r4OvwtnXohNk3gD9roM2PczmE5GApDIiAWdzbhEPzl1bvQ5JRutYfntlX87v2dZwMqmzohxnobTlf3VWc/VLHQD9r3euzElqby6fiAQUlREJSLZt8+7qvfx+3nr2F5UBcG73ZG7NyuCCXm1xuzSfJCiUlzqTXde87tzzxlfp7Ldc0OUCOOtm6Hk5REabzSkiRqmMSEArOlrBH/9vMy9/sQOvz/lXMC0pmhuHdWTssHTaJuhDLGiUHHAmvK55A3Z9cWx/dAvofx10HgntB0Fie+deOSISNlRGJCjsPFDCP7/axb+W7+ZwqTMXIcJlMahTS0b2aMN53dvQNy0Rl86YBIeD22DVLOdR+IPVXOPaQNpZVY9BznNCipmcItIsVEYkqByt8DJ/zT7++dUuVuw8VON7reM8ZHVtzbDOrRjSqRU9UxM0nBPofF7YthjWvQN7V0LuOrC9J74uIc05a5LcA1pmOGuZtMyAxA7gjmjm0CLS2Jq0jEyfPp0//OEP5OTkkJmZyXPPPcewYcNO+vrXX3+dKVOmsGPHDrp3784TTzzBZZddVuf3UxkJL7sOlPLp5v0s2bSfL7bkU1Je80MsISqCQZ1a0rtdIl3axNG1TRxd28TTItZjKLGcVsURyFnrFJO9q5zn/RuBk/znx3JDYpoztJPUwXm0SK8qLJ0hKd258Z+IBLQmKyNz5sxh3LhxzJgxg+HDh/Pss8/y+uuvs3HjRtq2PfGqiC+++ILzzjuPadOmccUVVzBr1iyeeOIJVq5cSb9+/Rr1l5HQU+H1sXLnIb7afpCvdxxk5c5DJ5QTv1ZxHrq1iadr23i6tY2ne9t42reMoXWch8ToSA31BJqyYti3GvatcoZ3Du2EQzvg8C7wlp3mYMspKgkpzvBPXDLEtXWWsY9p4cxXiU5ytqMSwJMAUfEQEdXUv5WIHKfJysjw4cMZOnQozz//PAA+n4/09HR+8YtfMGnSpBNeP3bsWEpKSnj//fer9/3oRz9i4MCBzJgxo1F/GQl9lV4fG3KKWLnrEFvyitm6v5ht+0vYV3D0lMe5XRYtYyNpEeshLiqC2Eg3sR43sVERxHncxEVFEBcVQXyUs+1xu/BEuIhwuYh0W0S6XVUPiwi3C4/bhauOt5l0WVbVAyzLqvMcTv8xrqoDwmbup+3DVZJHRNEe3FWPiKLvcRd+T0ThTtwFu3BVHmnYj3Z78EXGYXvi8UXGY3visT1x2BHR2C4PtjsSXJHYrkhstwdcEVX7PNjuCGyXB9zO93FHYltuwHKuIrIsbMtVtV31oGqfKwJcbuf1rgjnmB9ms6wTjj2BZWHjvJf/dbb/tT/cd/x+kTpolZJOVHRso/7Mun5+12tQtry8nBUrVjB58uTqfS6Xi1GjRrF06dJaj1m6dCnZ2dk19o0ePZq5c+ee9H3KysooKzv2f0aFhYX1iSkhLMLtol/7JPq1T6qxv6Ssku35JWzJK2ZLXjGb84rYkldMXmEZRWWVeH02+cXl5BeXG0ouDZMA9Kp6+NkkU0gHaz/JVgHJVgGtKSTZKqCVVUQiJSRZJSRRQqJVQgJHiLacydGWtxy3txyOHqrtzUTC2oYr3qLXkIuMvHe9ykh+fj5er5eUlJoz4FNSUtiwYUOtx+Tk5NT6+pycnJO+z7Rp03j44YfrE03CXFxURK0lBaC80seh0nIOFJdzuLScknIvpeWVHCn3Ulq1XVzmpaSskuKqR4XX5zwqbSp8PsorfVR6nW3/fl8dTiraOKuq27bzel8dz0Pato1tU31MXd4rnBTRkvW0PPEb9g+eq7jxEscR4jh67Nk69nUUFURSSQReIqnEQyURViWRVFbv91BJBJVEVr0mkkqc8xQ2LnxY2Lirni3Asny48VV/PwIfbry4OXGY0YLqn2XhHHey11hVP8/t/Nt1XAbnWAtwVb2v/7UidWHVdjaumQTkdPXJkyfXOJtSWFhIenq6wUQSzDwRLlISo0lJ1NolIiIn09Pge9erjCQnJ+N2u8nNza2xPzc3l9TU1FqPSU1NrdfrAaKiooiK0kQzERGRcFDHKXgOj8fD4MGDWbRoUfU+n8/HokWLyMrKqvWYrKysGq8HWLhw4UlfLyIiIuGl3sM02dnZ3HrrrQwZMoRhw4bx7LPPUlJSwvjx4wEYN24c7du3Z9q0aQBMnDiRkSNH8vTTT3P55Zcze/Zsli9fzosvvti4v4mIiIgEpXqXkbFjx7J//34eeughcnJyGDhwIAsWLKiepLpr1y5cx13zePbZZzNr1ix+85vf8Otf/5ru3bszd+7cOq8xIiIiIqFNy8GLiIhIk6jr53e95oyIiIiINDaVERERETFKZURERESMUhkRERERo1RGRERExCiVERERETFKZURERESMUhkRERERo1RGRERExKh6Lwdvgn+R2MLCQsNJREREpK78n9unW+w9KMpIUVERAOnp6YaTiIiISH0VFRWRlJR00u8Hxb1pfD4fe/fuJSEhAcuyGvxzCgsLSU9PZ/fu3brHTTPQ37t56e/dvPT3bl76ezevxvp727ZNUVERaWlpNW6i+0NBcWbE5XLRoUOHRvt5iYmJ+pe5Genv3bz0925e+ns3L/29m1dj/L1PdUbETxNYRURExCiVERERETEqrMpIVFQUU6dOJSoqynSUsKC/d/PS37t56e/dvPT3bl7N/fcOigmsIiIiErrC6syIiIiIBB6VERERETFKZURERESMUhkRERERo8KmjEyfPp2MjAyio6MZPnw4y5YtMx0pJE2bNo2hQ4eSkJBA27Ztufrqq9m4caPpWGHj8ccfx7Is7rnnHtNRQtaePXv42c9+RuvWrYmJiaF///4sX77cdKyQ5PV6mTJlCp07dyYmJoauXbvyyCOPnPY+J1J3S5YsYcyYMaSlpWFZFnPnzq3xfdu2eeihh2jXrh0xMTGMGjWKzZs3N3qOsCgjc+bMITs7m6lTp7Jy5UoyMzMZPXo0eXl5pqOFnE8//ZQJEybw5ZdfsnDhQioqKrj44ospKSkxHS3kff311/z5z39mwIABpqOErEOHDjFixAgiIyP54IMPWLduHU8//TQtW7Y0HS0kPfHEE7zwwgs8//zzrF+/nieeeIInn3yS5557znS0kFFSUkJmZibTp0+v9ftPPvkkf/rTn5gxYwZfffUVcXFxjB49mqNHjzZuEDsMDBs2zJ4wYUL1116v105LS7OnTZtmMFV4yMvLswH7008/NR0lpBUVFdndu3e3Fy5caI8cOdKeOHGi6Ugh6YEHHrDPOecc0zHCxuWXX27ffvvtNfZdc8019s0332woUWgD7Lfffrv6a5/PZ6emptp/+MMfqvcdPnzYjoqKsl977bVGfe+QPzNSXl7OihUrGDVqVPU+l8vFqFGjWLp0qcFk4aGgoACAVq1aGU4S2iZMmMDll19e499zaXzvvvsuQ4YM4frrr6dt27acddZZzJw503SskHX22WezaNEiNm3aBMDq1av57LPPuPTSSw0nCw/bt28nJyenxn9XkpKSGD58eKN/fgbFjfLORH5+Pl6vl5SUlBr7U1JS2LBhg6FU4cHn83HPPfcwYsQI+vXrZzpOyJo9ezYrV67k66+/Nh0l5G3bto0XXniB7Oxsfv3rX/P111/zy1/+Eo/Hw6233mo6XsiZNGkShYWF9OrVC7fbjdfr5dFHH+Xmm282HS0s5OTkANT6+en/XmMJ+TIi5kyYMIG1a9fy2WefmY4Ssnbv3s3EiRNZuHAh0dHRpuOEPJ/Px5AhQ3jssccAOOuss1i7di0zZsxQGWkC//rXv/jnP//JrFmz6Nu3L6tWreKee+4hLS1Nf+8QE/LDNMnJybjdbnJzc2vsz83NJTU11VCq0Hf33Xfz/vvv88knn9ChQwfTcULWihUryMvLY9CgQURERBAREcGnn37Kn/70JyIiIvB6vaYjhpR27drRp0+fGvt69+7Nrl27DCUKbffffz+TJk3ihhtuoH///txyyy3ce++9TJs2zXS0sOD/jGyOz8+QLyMej4fBgwezaNGi6n0+n49FixaRlZVlMFlosm2bu+++m7fffpuPP/6Yzp07m44U0i666CLWrFnDqlWrqh9Dhgzh5ptvZtWqVbjdbtMRQ8qIESNOuFR906ZNdOrUyVCi0FZaWorLVfNjyu124/P5DCUKL507dyY1NbXG52dhYSFfffVVo39+hsUwTXZ2NrfeeitDhgxh2LBhPPvss5SUlDB+/HjT0ULOhAkTmDVrFu+88w4JCQnV44pJSUnExMQYThd6EhISTpiPExcXR+vWrTVPpwnce++9nH322Tz22GP89Kc/ZdmyZbz44ou8+OKLpqOFpDFjxvDoo4/SsWNH+vbtyzfffMMzzzzD7bffbjpayCguLmbLli3VX2/fvp1Vq1bRqlUrOnbsyD333MPvf/97unfvTufOnZkyZQppaWlcffXVjRukUa/NCWDPPfec3bFjR9vj8djDhg2zv/zyS9ORQhJQ6+Pll182HS1s6NLepvXee+/Z/fr1s6OiouxevXrZL774oulIIauwsNCeOHGi3bFjRzs6Otru0qWL/eCDD9plZWWmo4WMTz75pNb/Zt966622bTuX906ZMsVOSUmxo6Ki7IsuusjeuHFjo+ewbFtL2YmIiIg5IT9nRERERAKbyoiIiIgYpTIiIiIiRqmMiIiIiFEqIyIiImKUyoiIiIgYpTIiIiIiRqmMiIiIiFEqIyJygsWLF2NZFocPH27W933llVdo0aLFGf2MHTt2YFkWq1atOulrTP1+IlI7lRER4fzzz+eee+4xHUNEwpTKiIg0ivLyctMRRCRIqYyIhLnbbruNTz/9lD/+8Y9YloVlWezYsQOAFStWMGTIEGJjYzn77LPZuHFj9XG//e1vGThwIH/5y1/o3Lkz0dHRABw+fJg77riDNm3akJiYyIUXXsjq1aurj1u9ejUXXHABCQkJJCYmMnjwYJYvX14j04cffkjv3r2Jj4/nkksuYd++fdXf8/l8/O53v6NDhw5ERUUxcOBAFixYcMrfcf78+fTo0YOYmBguuOCC6t9PRAKDyohImPvjH/9IVlYWd955J/v27WPfvn2kp6cD8OCDD/L000+zfPlyIiIiTrh1+5YtW3jzzTd56623qudoXH/99eTl5fHBBx+wYsUKBg0axEUXXcTBgwcBuPnmm+nQoQNff/01K1asYNKkSURGRlb/zNLSUp566in+/ve/s2TJEnbt2sV9991XI+/TTz/NU089xbfffsvo0aO58sor2bx5c62/3+7du7nmmmsYM2YMq1at4o477mDSpEmN+ScUkTPV6PcBFpGgM3LkSHvixInVX/tvK/5///d/1fvmzZtnA/aRI0ds27btqVOn2pGRkXZeXl71a/7973/biYmJ9tGjR2v8/K5du9p//vOfbdu27YSEBPuVV16pNcfLL79sA/aWLVuq902fPt1OSUmp/jotLc1+9NFHaxw3dOhQ+7/+679s27bt7du324D9zTff2LZt25MnT7b79OlT4/UPPPCADdiHDh061Z9FRJqJzoyIyEkNGDCgertdu3YA5OXlVe/r1KkTbdq0qf569erVFBcX07p1a+Lj46sf27dvZ+vWrQBkZ2dzxx13MGrUKB5//PHq/X6xsbF07dq1xvv637OwsJC9e/cyYsSIGseMGDGC9evX1/o7rF+/nuHDh9fYl5WVVee/gYg0vQjTAUQkcB0/fGJZFuDM2fCLi4ur8fri4mLatWvH4sWLT/hZ/kt2f/vb33LTTTcxb948PvjgA6ZOncrs2bP5yU9+csJ7+t/Xtu3G+HVEJEDpzIiI4PF48Hq9Z/xzBg0aRE5ODhEREXTr1q3GIzk5ufp1PXr04N577+Wjjz7immuu4eWXX67Tz09MTCQtLY3PP/+8xv7PP/+cPn361HpM7969WbZsWY19X375ZT1/MxFpSiojIkJGRgZfffUVO3bsID8/v8bZj/oYNWoUWVlZXH311Xz00Ufs2LGDL774ggcffJDly5dz5MgR7r77bhYvXszOnTv5/PPP+frrr+ndu3ed3+P+++/niSeeYM6cOWzcuJFJkyaxatUqJk6cWOvr//M//5PNmzdz//33s3HjRmbNmsUrr7zSoN9PRJqGyoiIcN999+F2u+nTpw9t2rRh165dDfo5lmUxf/58zjvvPMaPH0+PHj244YYb2LlzJykpKbjdbg4cOMC4cePo0aMHP/3pT7n00kt5+OGH6/wev/zlL8nOzuZXv/oV/fv3Z8GCBbz77rt079691td37NiRN998k7lz55KZmcmMGTN47LHHGvT7iUjTsGwNxoqIiIhBOjMiIiIiRqmMiIiIiFEqIyIiImKUyoiIiIgYpTIiIiIiRqmMiIiIiFEqIyIiImKUyoiIiIgYpTIiIiIiRqmMiIiIiFEqIyIiImLU/wcMZbk8FJgViwAAAABJRU5ErkJggg==", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "flip_flop_index_proportion_exceeding(fcst1, \"lead_day\", np.arange(0.1, 10, 0.1)).plot()\n", + "flip_flop_index_proportion_exceeding(fcst2, \"lead_day\", np.arange(0.1, 10, 0.1)).plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "id": "c9f233c3", + "metadata": { + "ExecuteTime": { + "end_time": "2023-11-09T05:07:23.729926Z", + "start_time": "2023-11-09T05:07:23.484614Z" + } + }, + "outputs": [ + { + "data": { + "text/plain": [ + "[]" + ] + }, + "execution_count": 9, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAG0CAYAAADO5AZFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAABRbElEQVR4nO3dd3hUVeLG8e+dSQ9JCAQSAqEovTeJFMUSxYZtVVwLLK4dFcjiT1ERdVXEwuoqiuLae0VURDAKCtIhKNJ7KEkIJZW0mfn9cSEYDUjCJCeZeT/PM09ubu7MvJNV8+65955jeTweDyIiIiI+wmE6gIiIiIg3qdyIiIiIT1G5EREREZ+iciMiIiI+ReVGREREfIrKjYiIiPgUlRsRERHxKSo3IiIi4lNUbkRERMSnqNyIiIiITzFabn788UcGDx5MfHw8lmUxbdq0v3zOnDlz6NmzJ8HBwbRu3Zo33nij2nOKiIhI3RFg8s3z8/Pp1q0bN9xwA5dffvlfHr9lyxYuvPBCbr31Vt59911SUlK48cYbadKkCYMGDTqu93S73ezatYuIiAgsyzrRjyAiIiI1wOPxkJubS3x8PA7HscdmrNqycKZlWXz++edceumlRz3mnnvu4euvv2bVqlVl+66++moOHDjAzJkzj+t9duzYQUJCwonGFREREQPS0tJo1qzZMY8xOnJTWQsWLCApKancvkGDBjFq1KijPqeoqIiioqKy7w93ubS0NCIjI6slp4iIiHhXTk4OCQkJRERE/OWxdarcpKenExsbW25fbGwsOTk5HDx4kNDQ0D89Z8KECTz88MN/2h8ZGalyIyIiUscczyUlPn+31NixY8nOzi57pKWlmY4kIiIi1ahOjdzExcWRkZFRbl9GRgaRkZEVjtoABAcHExwcXBPxREREpBaoUyM3ffv2JSUlpdy+2bNn07dvX0OJREREpLYxWm7y8vJITU0lNTUVsG/1Tk1NZfv27YB9Smno0KFlx996661s3ryZ//u//2Pt2rW8+OKLfPTRR4wePdpEfBEREamFjJabpUuX0qNHD3r06AFAcnIyPXr04MEHHwRg9+7dZUUHoFWrVnz99dfMnj2bbt268cwzz/Dqq68e9xw3IiIi4vtqzTw3NSUnJ4eoqCiys7N1t5SIiEgdUZm/33XqmhsRERGRv6JyIyIiIj5F5UZERER8isqNiIiI+BSVGxEREfEpKjciIiLiU1RuvGnzXCjKM51CRETEr6nceMvOZfDulfC/c2DvJtNpRERE/JbKjbe43RBaHzJXw9QzYcN3phOJiIj4JZUbb0k4BW6eC81OgcJsePcK+OkZ8K8JoEVERIxTufGmyCbwj6+h5zDAAymPwMfDdB2OiIhIDVK58baAYLj4v3DRf8ARCKu/gNfOgwNpppOJiIj4BZWb6tL7BvjHVxDeCDJ+ta/DSVtsOpWIiIjPU7mpTs1PhZu+h9gukL8H3rgQVn5gOpWIiIhPU7mpbvWbww0zof1F4CqGz2+B7x7ShcYiIiLVROWmJgTXg6vehtPG2N/P+w98c48KjoiISDVQuakpDgecPQ4G/xewYPHLKjgiIiLVQOWmpvUaBhc/jwqOiIhI9QgwHcAv9bze/jr9TrvgAJw/ESzLXCYREREfoZEbU3peX34E57vxphOJiIj4BJUbk3peb0/4BzD/OUh932weERERH6ByY1rPoXD63fb2lyNhx1KzeUREROo4lZva4Iz7oN2F4CqCD66FnF2mE4mIiNRZKje1gcMBl78MjTpAXrpdcEoOmk4lIiJSJ6nc1BbBEfD39yE0GnYth+l36RZxERGRKlC5qU0atIKr3gLLCb9+BIunmk4kIiJS56jc1DatTodzH7W3Z90Pu1KNxhEREalrVG5qo1NvO3SBcTF8/A8ozDGdSEREpM5QuamNLAsueQGiEmD/FvhS19+IiIgcL5Wb2iqsAVzxOjgC4LfPYelrphOJiIjUCSo3tVnCKXD2oWUZZo6F3b+YzSMiIlIHqNzUdn3vgDaD7An+Ph4GBw+YTiQiIlKrqdzUdg4HXDbFvv5m32b47GZwu02nEhERqbVUbuqCsAYw5G1wBsOGb2HuE6YTiYiI1FoqN3VFfA8Y/Jy9PXcirJ1hNo+IiEgtpXJTl3T/O/S52d7+/BbI2mA2j4iISC2kclPXDHocmveDohz44BooyjWdSEREpFZRualrnIFw1ZsQEQ9Z62HmvaYTiYiI1CoqN3VRvcZwxf8AC1a8AxtTTCcSERGpNVRu6qoW/SDxFnv7y5E6PSUiInKIyk1ddvaDUL8FZKfB7AdNpxEREakVVG7qsqBwe4FNsNee2jzXbB4REZFaQOWmrmt1OvS+wd6eficU5ZnNIyIiYpjKjS9Iehgim8GBbZDyiOk0IiIiRqnc+IKQSLj40OzFi1+GrfPM5hERETFI5cZXtE6CnkPt7Wm36+4pERHxWyo3vuTcxyCquX16atY402lERESMULnxJSGRcOlke3vZ67DxO7N5REREDFC58TWtTofEW+3tL+6Eg/vN5hEREalhKje+6Ozx0OBkyN0F32jtKRER8S8qN74oKAwumwKWA375ANZ8aTqRiIhIjVG58VUJfaD/SHv7q9GQv9dsHhERkRqicuPLzhgLjTpA/h6Y8S/TaURERGqEyo0vCwiGy14Cywm/fQ6rPjOdSEREpNqp3Pi6+B5w2qFRm6//BXmZZvOIiIhUM5Ubf3D63RDbBQ7us6+/8XhMJxIREak2Kjf+ICDIPj3lCIC1X8GqT00nEhERqTYqN/4irgsMvMfe/vpfkJthNo+IiEg1UbnxJwNGQ5NuUHgAvr3PdBoREZFqYbzcTJ48mZYtWxISEkJiYiKLFy8+5vHPPvss7dq1IzQ0lISEBEaPHk1hYWENpa3jnIEw+Dl7cr9Vn8Cm700nEhER8Tqj5ebDDz8kOTmZ8ePHs3z5crp168agQYPIzKz4jp733nuPe++9l/Hjx7NmzRr+97//8eGHH3LffRqFOG7xPaDPzfb21/+CEhVDERHxLUbLzaRJk7jpppsYPnw4HTt2ZMqUKYSFhfHaa69VePzPP/9M//79ueaaa2jZsiXnnnsuf//73/9ytEf+4Mz7IaIJ7NsM8yaZTiMiIuJVxspNcXExy5YtIykp6UgYh4OkpCQWLFhQ4XP69evHsmXLysrM5s2bmTFjBhdccMFR36eoqIicnJxyD78XEgnnPWFvz/sPZG0wm0dERMSLjJWbrKwsXC4XsbGx5fbHxsaSnp5e4XOuueYaHnnkEQYMGEBgYCAnn3wyZ5xxxjFPS02YMIGoqKiyR0JCglc/R53V8RJofQ64iuHrZM19IyIiPsP4BcWVMWfOHB5//HFefPFFli9fzmeffcbXX3/Nv//976M+Z+zYsWRnZ5c90tLSajBxLWZZcMFTEBACW36EXz4ynUhERMQrAky9cUxMDE6nk4yM8vOtZGRkEBcXV+Fzxo0bx/XXX8+NN94IQJcuXcjPz+fmm2/m/vvvx+H4c1cLDg4mODjY+x/AFzRoBQP/D1IegVn3Q7vzICTKdCoREZETYmzkJigoiF69epGSklK2z+12k5KSQt++fSt8TkFBwZ8KjNPpBMCj0ypV0/dOaNjGXjl8zkTTaURERE6Y0dNSycnJTJ06lTfffJM1a9Zw2223kZ+fz/DhwwEYOnQoY8eOLTt+8ODBvPTSS3zwwQds2bKF2bNnM27cOAYPHlxWcqSSAoKOXFy8+GXIXGs2j4iIyAkydloKYMiQIezZs4cHH3yQ9PR0unfvzsyZM8suMt6+fXu5kZoHHngAy7J44IEH2LlzJ40aNWLw4ME89thjpj6Cb2iTBO0ugHUz4Jv/g6Ff2NfkiIiI1EGWx8/O5+Tk5BAVFUV2djaRkZGm49Qe+7bA5ERwFcFVb0PHi00nEhERKVOZv9916m4pqUYNWkH/u+ztb++HkoNm84iIiFSRyo0cMWA0RDaD7O0w/znTaURERKpE5UaOCAqHcw/NGTTvP7B/m9k8IiIiVaByI+V1ugxangalhfDdQ6bTiIiIVJrKjZRnWXDeBMCC3z6DtCWmE4mIiFSKyo38WVwX6H6tvT3rfq07JSIidYrKjVTsrAcgMAzSFsHqL0ynEREROW4qN1KxyCbQ79Ct4d+Nh9Iis3lERESOk8qNHF3/u6BeHOzfCounmk4jIiJyXFRu5OiCwuGs++3tH5+Egn1m84iIiBwHlRs5tu7XQmxnKMyGuU+aTiMiIvKXVG7k2BzOIxP7LZkK+zabzSMiIvIXVG7kr518FrROAncppPzbdBoREZFjUrmR45P0EGUT++1cZjqNiIjIUancyPGJ6wJdh9jbs8drYj8REam1VG7k+J11PziDYOtPsDHFdBoREZEKqdzI8avfHPrcbG9/Nx7cLrN5REREKqByI5Vz2r8gOAoyVsEvH5lOIyIi8icqN1I5YQ3gtNH29g+PQUmh2TwiIiJ/oHIjlZd4K0TEQ3aaPfeNiIhILaJyI5UXGApn3mdvz3sWivONxhEREfk9lRupmm5XQ3RLKMiCpa+bTiMiIlJG5UaqxhloX1wMMP85KDloNo+IiMghKjdSdV2vhqjmkJ8Jy94wnUZERARQuZETERAEpyXb2/Oe1Z1TIiJSK6jcyInpfg1ENoO8dFjxtuk0IiIiKjdyggKCYcAoe3vef6C0yGgcERERlRs5cT2H2vPe5OyEFe+YTiMiIn5O5UZO3J9Gb4qNxhEREf+mciPe0XMo1Iu1Zy3WnVMiImKQyo14R2AoDPw/e3vuE1CYbTaPiIj4LZUb8Z6ew6BhGyjYa98aLiIiYoDKjXiPMxDOedjeXvgiZO8wm0dERPySyo14V7sLoHk/KC2E7x8znUZERPyQyo14l2XBuY/a2yvfh/RfzeYRERG/o3Ij3tesF3S6HPDArHGm04iIiJ9RuZHqcfaD4AiEzT/Axu9MpxERET+iciPVo0Er6HOzvT17PLjdZvOIiIjfULmR6nP6GAiOhIxVsGa66TQiIuInVG6k+oQ1gMRb7e25EzV6IyIiNULlRqpX39shOAoyV8PqaabTiIiIH1C5keoVGm0XHDg0euMym0dERHyeyo1Uv1Nvg5Ao2LMWfvvcdBoREfFxKjdS/UKioO8d9rZGb0REpJqp3EjNSLwVQupD1npY9anpNCIi4sNUbqRmhERCv9+N3rhKzeYRERGfpXIjNafPLfYFxns3wqpPTKcREREfpXIjNSckEvrdaW//+LSuvRERkWqhciM165Sb7Gtv9m7QnVMiIlItVG6kZoVEQt8R9vaPT2nWYhER8TqVG6l5fW62Zy3es1ZrTomIiNep3EjNC60PibfY2xq9ERERL1O5ETNOvQ2C6tkrhq//xnQaERHxISo3YkZYA/v0FNjz3ng8ZvOIiIjPULkRc/reAYHhsHslbJhtOo2IiPgIlRsvSU07wD9eX8xd768wHaXuCG8Ip9xgb2v0RkREvETlxkscFsxZt4fv12ZS6tIFsset310QEAo7l8LGFNNpRETEB6jceEmn+CjqhwWSV1TKyh0HTMepO+o1hlP+aW//8JhGb0RE5ISp3HiJ02HR7+SGAMzbsNdwmjqm/ygIDINdy2H9t6bTiIhIHWe83EyePJmWLVsSEhJCYmIiixcvPubxBw4cYMSIETRp0oTg4GDatm3LjBkzaijtsfVvHQPA/I1ZhpPUMfUaQZ+b7O05j2v0RkRETojRcvPhhx+SnJzM+PHjWb58Od26dWPQoEFkZmZWeHxxcTHnnHMOW7du5ZNPPmHdunVMnTqVpk2b1nDyig04VG6Wb99PflGp4TR1TL+R9rw3u1fCutpRVkVEpG4yWm4mTZrETTfdxPDhw+nYsSNTpkwhLCyM1157rcLjX3vtNfbt28e0adPo378/LVu2ZODAgXTr1q2Gk1esRcNwEhqEUur2sGiLTk1VSnjDI/Pe/DBBsxaLiEiVGSs3xcXFLFu2jKSkpCNhHA6SkpJYsGBBhc+ZPn06ffv2ZcSIEcTGxtK5c2cef/xxXC7XUd+nqKiInJycco/qdHj0RtfdVEG/OyEoAjJ+hbVfmU4jIiJ1lLFyk5WVhcvlIjY2ttz+2NhY0tPTK3zO5s2b+eSTT3C5XMyYMYNx48bxzDPP8Oijjx71fSZMmEBUVFTZIyEhwauf44903c0JCGtgL8sAMEejNyIiUjXGLyiuDLfbTePGjXnllVfo1asXQ4YM4f7772fKlClHfc7YsWPJzs4ue6SlpVVrxn4nx2BZsC4jl8zcwmp9L5/U93Z7xfDM1bD6c9NpRESkDjJWbmJiYnA6nWRkZJTbn5GRQVxcXIXPadKkCW3btsXpdJbt69ChA+np6RQXF1f4nODgYCIjI8s9qlOD8CA6xdvvodGbKgiNtgsOwNwnwX30U44iIiIVMVZugoKC6NWrFykpR2aldbvdpKSk0Ldv3wqf079/fzZu3Ij7d6cr1q9fT5MmTQgKCqr2zMerv667OTGn3gYhUbBnLfym0RsREakco6elkpOTmTp1Km+++SZr1qzhtttuIz8/n+HDhwMwdOhQxo4dW3b8bbfdxr59+xg5ciTr16/n66+/5vHHH2fEiBGmPkKFBvzuuhuP5mypvJAoe1FN0OiNiIhUWoDJNx8yZAh79uzhwQcfJD09ne7duzNz5syyi4y3b9+Ow3GkfyUkJPDtt98yevRounbtStOmTRk5ciT33HOPqY9QoVNaNiAowEF6TiGb9uTTunE905HqnsRbYMELkLXOHr3pcoXpRCIiUkdYHj8bWsjJySEqKors7Oxqvf7m2lcXMn/jXh6+uBPD+rWstvfxaXOfgh8ehZi2cPtCcDj/+jkiIuKTKvP3u07dLVWXHL7u5qcNuqi4yhJvgZD6kLUeVn1mOo2IiNQRKjfV5PB1Nws376XUpflaqiQkEvodvvZmoq69ERGR46JyU006xUcRFRpIXlEpK3dkm45Td/W5xb49fO8G+PUT02lERKQOULmpJk6HRf/WDQGYu36P4TR1WEjkkTunfnwSXFqQVEREjk3lphqd1d6+62vWbxUvJyHHKfHw6M1GzXsjIiJ/SeWmGp3dvjFOh8Xa9Fy27y0wHafuCo6AUw/NWjz/OfCvG/xERKSSVG6qUXR4EH1aNgBg1mqN3pyQU26EwDB7xfBN35tOIyIitZjKTTU7t9PhU1MZf3GkHFNYA+g5zN6e/6zRKCIiUrup3FSzczra5Wbptn1k5RUZTlPH9R0BjgDY8iPsXG46jYiI1FJeKTcZGRls377dGy/lc5pFh9G5aSRuD6Ss0ejNCamfAJ0PLcMw/zmzWUREpNaqVLnJzc3luuuuo0WLFgwbNozi4mJGjBhBkyZNaNWqFQMHDiQnJ6e6stZZ53aMA3Rqyiv632V/XTMd9m4ym0VERGqlSpWb++67j2XLljFmzBi2b9/OVVddxY8//shPP/3EDz/8QFZWFhMnTqyurHXW4etuftqYRX6R5mk5IbGdoM254HHbC2uKiIj8QaXKzRdffMGLL77InXfeybvvvsv06dOZMGEC/fv35/TTT+fJJ5/k008/ra6sdVa72AhaNAyjuNStCf28of8o++uKdyEv02gUERGpfSpVbjIzM2ndujUA8fHxhIaG0rZt27Kfd+7cmbS0NO8m9AGWZXFuR03o5zUt+kGzU8BVBIteNp1GRERqmUqVm4YNG7Jnz5GRh0suuYT69euXfZ+Xl0dwcLDXwvmSczvZ192krM2kRAtpnhjLgv4j7e0lU6Eoz2weERGpVSpVbrp27cqSJUvKvn/vvfdo3Lhx2fdLliyhQ4cO3kvnQ3o2jyamXhC5haUs3LzXdJy6r92F0OAkKMyGle+bTiMiIrVIpcrNu+++y5AhQ47689jYWB577LETDuWLnA6LpA6a0M9rHI4jSzIsfBHcLrN5RESk1qhUuWnQoEG501B/dP7553PGGWecYCTfVTZb8ep03G6tj3TCul8DIfVh32ZYP9N0GhERqSWqNImf213xNSMej0eT+R1Dv5NjCA9ykpFTxC87s03HqfuCwqH3cHt7wYtms4iISK1RqXKTk5PDVVddRXh4OLGxsTz44IO4XEdOB2RmZtKqVSuvh/QVIYFOzmxvX6M0c5XumvKKPjfbSzJsmwe7VphOIyIitUClys24ceNYuXIlb7/9No899hhvvfUWl1xyCcXFxWXHeDw63XIs53W275qauWq3flfeEBkPnS63tzV6IyIiVLLcTJs2jZdffpkrrriCG2+8kaVLl7Jnzx4GDx5MUZG9KKRlWdUS1Fec0a4xQQEOtu4tYH2GbmH2ir6HLiz+7TPI3mk2i4iIGFepcrNnzx5atGhR9n1MTAzfffcdubm5XHDBBRQUFHg9oK+pFxzA6W1iAJ2a8pr4HtBiALhLYfErptOIiIhhlSo3zZs3Z82aNeX2RUREMGvWLA4ePMhll13m1XC+atChCf1marZi7+k7wv667HVN6ici4ucqVW7OPfdcXn/99T/tr1evHt9++y0hISFeC+bLkjrE4nRYrNmdw7a9+abj+Ia25x2Z1G/F26bTiIiIQZUqNw8//DAPPfRQhT+LiIhg9uzZfP/9997I5dOiw4NIbNUAgG81euMdDgf0u9Penv8clBSazSMiIsZUqtxER0fTqVOno/48IiKCgQMHnnAof3D4rqlvNVux93S/FiKbQu5uWP6W6TQiImJIlSbxO5qMjAweeeQRb76kzzq3o11ulm3bT2aORhm8IiAYTku2t+dN0uiNiIif8mq5SU9P5+GHH/bmS/qsuKgQejSvD8C3qzV64zU9rj8yeqNrb0RE/FJAZQ7+5ZdfjvnzdevWnVAYf3NepzhWbD/At6vSuf7UFn/9BPlrAcEwYDTMGAM/TbLLTqAudBcR8SeVKjfdu3fHsqwKZ9Y9vF+T+B2/QZ3imPDNWhZs3suBgmLqhwWZjuQbeg61i03uLnv0ps9NphOJiEgNqvSq4FOnTmXLli1/emzevJmvvvqqunL6pJYx4bSPi8Dl9vDdmkzTcXzH76+9+WkSlBaZzSMiIjWqUuWmV69e7Nq1ixYtWlT4aNq0qdZLqqTfrzUlXtRzKETE26M3unNKRMSvVKrc3HrrrbRs2fKoP2/evHmFk/zJ0Z3fuQkAP67PIrewxHAaH1Luzqn/aPRGRMSPVKrcXHbZZVx33XVH/Xl0dDTDhg074VD+pG1sPU5uFE6xy02KTk15V4/r7dGbnJ2Q+q7pNCIiUkO8ciu4y+UiNTWV/fv3e+Pl/IplWVzQxR69mfGrTk15VWAI9B9pb//0H3BpZExExB9UqdyMGjWK//3vf4BdbAYOHEjPnj1JSEhgzpw53sznFw6Xmznr95BXVGo4jY/pNQzCG0P2dlj5gek0IiJSA6pUbj755BO6desGwJdffsmWLVtYu3Yto0eP5v777/dqQH/QPi6CVjHhFJe6+X6tTk15VWAo9L/L3v7pGXCpPIqI+LoqlZusrCzi4uy7fGbMmMGVV15J27ZtueGGG/j111+9GtAf2KemDv0+f9GpKa/rfQOENYT9W2DVp6bTiIhINatSuYmNjWX16tW4XC5mzpzJOeecA0BBQQFOp9OrAf3F4VNTP6zLJF+nprwrKBz63mFv//Q0uF1m84iISLWqUrkZPnw4V111FZ07d8ayLJKSkgBYtGgR7du392pAf9GxSSQtGoZRVOrmh3U6NeV1fW6CkPqQtR5WTzOdRkREqlGVys1DDz3Eq6++ys0338z8+fMJDg4GwOl0cu+993o1oL/QXVPVLDgCTr3d3v7xaXC7zeYREZFqY3mqMKXw5s2bOemkk6ojT7XLyckhKiqK7OxsIiMjTccp59cd2Qx+YR6hgU6WjUsiLKhSS3/JXzl4AJ7tAkU5MOQd6DDYdCIRETlOlfn7XaWRm9atW3PmmWfyzjvvUFhYWKWQ8medm0aS0CCUgyUu5qzbYzqO7wmtD4m32NtzJ4KWChER8UlVKjfLly+na9euJCcnExcXxy233MLixYu9nc3vWJbFBZ11aqpanXo7BNWD9F9h/UzTaUREpBpUqdx0796d5557jl27dvHaa6+xe/duBgwYQOfOnZk0aRJ79mjUoaoOX3fz/dpMCkt0V4/XhTWwLy4Gjd6IiPioE1p+ISAggMsvv5yPP/6YiRMnsnHjRsaMGUNCQgJDhw5l926NPlRW12ZRNK0fSkGxix80oV/16HsHBIbBrhWwYbbpNCIi4mUnVG6WLl3K7bffTpMmTZg0aRJjxoxh06ZNzJ49m127dnHJJZd4K6ffsCyLi7rZozfTV+4ynMZHhcfAKf+0tzV6IyLic6pUbiZNmkSXLl3o168fu3bt4q233mLbtm08+uijtGrVitNOO4033niD5cuXezuvX7i4WzwAKWszyS3UYo/Vot9dEBAKO5fCpu9NpxERES+qUrl56aWXuOaaa9i2bRvTpk3joosuwuEo/1KNGzcuW1xTKqdjk0hObmSvNTXrtwzTcXxTvcb2sgyg0RsRER9TpXKzYcMGxo4dS5MmTY56TFBQEMOGDatyMH9mWRYXd2sK6NRUtep/FziDIW0RbJlrOo2IiHjJCV1zU1BQwNq1a/nll1/KPeTEXdzdPjU1b2MWe/OKDKfxURFx0Osf9vbcJ41GERER76lSudmzZw8XXnghERERdOrUiR49epR7yIlrFRNOl6ZRuNwezXlTnfqPBGcQbJsPW34ynUZERLygSuVm1KhRZGdns2jRIkJDQ5k5cyZvvvkmbdq0Yfr06d7O6LcOX1isU1PVKKop9Bxqb895wmwWERHxiiqVm++//55JkybRu3dvHA4HLVq04LrrruPJJ59kwoQJ3s7oty7q1gTLgiVb97PzwEHTcXzXgNGHRm/mafRGRMQHVKnc5Ofn07hxYwCio6PLZiTu0qWLbv/2oiZRoZzSsgEAX2n0pvpENdPojYiID6lSuWnXrh3r1q0DoFu3brz88svs3LmTKVOmHPMOKqk8nZqqIQOSNXojIuIjqlRuRo4cWba0wvjx4/nmm29ISEjgueee4/HHH/dqQH93QZcmBDgsftuVw6Y9eabj+K6optDz0NQFGr0REanTLI/nxGcvO3xLePPmzYmJifFGrmqTk5NDVFQU2dnZREZGmo5zXP7x+mLmrNvDyLPbMPqctqbj+K7snfDf7uAqhmFfQavTTCcSEZFDKvP3O+B4XzQ5Ofm4A0yaNOm4jwWYPHkyTz31FOnp6XTr1o3nn3+ePn36/OXzPvjgA/7+979zySWXMG3atEq9Z11ycbd45qzbw/SVuxiV1AbLskxH8k2HR2+WTLVHb1RuRETqpOMuNytWrCj3/fLlyyktLaVdu3YArF+/HqfTSa9evSoV4MMPPyQ5OZkpU6aQmJjIs88+y6BBg1i3bl3ZRcsV2bp1K2PGjOG003z/D9C5neIIDVzFlqx8lm8/QK8W0aYj+a4Bo2H5m0euvVHBERGpc477mpsffvih7DF48GAGDhzIjh07WL58OcuXLyctLY0zzzyTCy+8sFIBJk2axE033cTw4cPp2LEjU6ZMISwsjNdee+2oz3G5XFx77bU8/PDDnHTSSZV6v7qoXnAA53eJA+CTZTsMp/Fxv7/25ofHteaUiEgdVKULip955hkmTJhAdPSREYTo6GgeffRRnnnmmeN+neLiYpYtW0ZSUtKRQA4HSUlJLFiw4KjPe+SRR2jcuDH//Oc/qxK/TrqiVzPAviW8sMRlOI2PGzAaAkJg+8+w7hvTaUREpJKqVG5ycnLK5rb5vT179pCbm3vcr5OVlYXL5SI2Nrbc/tjYWNLT0yt8zrx58/jf//7H1KlTj+s9ioqKyMnJKfeoi05t1ZCm9UPJLSrl298q/t2Il0Q1hVNvs7e/Gw+uUrN5RESkUqpUbi677DKGDx/OZ599xo4dO9ixYweffvop//znP7n88su9nbFMbm4u119/PVOnTj3uu7ImTJhAVFRU2SMhIaHa8lUnh8Pib4dGb3RqqgYMGA1hDSFrvX0NjoiI1BlVKjdTpkzh/PPP55prrqFFixa0aNGCa665hvPOO48XX3zxuF8nJiYGp9NJRkZGuf0ZGRnExcX96fhNmzaxdetWBg8eTEBAAAEBAbz11ltMnz6dgIAANm3a9KfnjB07luzs7LJHWlpa5T9wLfG3nk0Be6Xw3dlajqFahUTBwHvs7TkToOj4RyRFRMSsKpWbsLAwXnzxRfbu3cuKFStYsWIF+/bt48UXXyQ8PPy4XycoKIhevXqRkpJSts/tdpOSkkLfvn3/dHz79u359ddfSU1NLXtcfPHFnHnmmaSmplY4KhMcHExkZGS5R13VomE4fVo1wOOBz5bvNB3H9/UaDg1Ogvw9MP+/ptOIiMhxqlK5OSw8PJyuXbvStWvXSpWa30tOTmbq1Km8+eabrFmzhttuu438/HyGDx8OwNChQxk7diwAISEhdO7cudyjfv36RERE0LlzZ4KCgk7k49QJhy8s/nTZDrww/6IcS0AQJD1kby94AXJ2G40jIiLH54TKjTcMGTKEp59+mgcffJDu3buTmprKzJkzyy4y3r59e9lSD2IvxxAa6GRzVj4r0g6YjuP7OlwMCYlQUgA/PGY6jYiIHAevLL9Ql9TF5Rf+KPmjVD5bvpNrEpvz+GVdTMfxfWmL4X/nABbcNh9iO5lOJCLidyrz99v4yI1U3uFTU19qzpuakdAHOl4CeOC7h0ynERGRv6ByUweVzXlTqDlvaszZ48ERABtmwZYfTacREZFjULmpgxwOq2z05v3F2w2n8RMNT7bvngKY/SC43WbziIjIUanc1FFX90nAYcHCzfvYmKk5WGrEwHsgqB7sWgGrPzedRkREjkLlpo5qEhXK2R3sO8reXaTRmxpRrxH0H2lvpzwCpcVm84iISIVUbuqw605tAdhz3hws1oXFNaLvCKgXC/u3wtKjr1wvIiLmqNzUYae1jiGhQSg5haV89csu03H8Q1A4nHGvvf3jk1BYNxdiFRHxZSo3dZjDYXFNH3v05h2dmqo5PYZCwzZQsBfmP2c6jYiI/IHKTR13Ze9mBDotVqYdYNXObNNx/IMz4HfLMkyG/duMxhERkfJUbuq4mHrBnNe5CaALi2tU+wuh5WlQehC+Tgb/muhbRKRWU7nxAdcmNgfgi9Sd5BaWGE7jJywLLnoWnMGw8Tv49RPTiURE5BCVGx+Q2KoBrRvXo6DYxbQVO03H8R8xrWHg3fb2zHsgf6/ZPCIiAqjc+ATLsspGb95ZuB0/WwvVrH4joXFH++LiWfebTiMiIqjc+IzLezYjJNDBuoxclm7bbzqO/wgIgoufByxY+T5s+t50IhERv6dy4yOiQgO5pFtTAN5aoLt3alSz3tDnZnv7q9FQXGA2j4iIn1O58SHX97XnvJm5ajeZuYWG0/iZs8dBZDN75uK5T5hOIyLi11RufEjnplH0bF6fEpeHDxanmY7jX4Ij4MJn7O0FkyHjN7N5RET8mMqNjxnatyUA7y3aTqnLbTaMv2l3HrS/CNyl8FUyuPX7FxExQeXGx5zfJY6G4UGk5xQye3WG6Tj+5/yJEBgOaQsh9R3TaURE/JLKjY8JDnBydZ8EQBcWGxHVDM4ca2/PflBz34iIGKBy44OuSWyBw4IFm/eyISPXdBz/k3grxHaGg/vtgiMiIjVK5cYHNa0fSlKHWADeXqjRmxrnDISL/mNvp74D2342m0dExM+o3PiowxcWf7Z8J3lFpWbD+KOEPtBzmL391WgoLTabR0TEj6jc+Kj+rRtyUqNw8opK+Xz5DtNx/FPSQxDWEPashUVTTKcREfEbKjc+yrIsrj/VntTv7YXbtN6UCWEN4JxH7O25EyFnt9k8IiJ+QuXGh/2tVzNCA52sz8hj8ZZ9puP4p27XQNPeUJyni4tFRGqIyo0PiwwJ5NIeh9ab0oXFZjgccMFTgAW/fqSLi0VEaoDKjY87fGrq21XpZOZovSkjmvaEXocuLp5xN7h0gbeISHVSufFxHeMj6d0imlK3hw+WaL0pY856EELqQ8YqWPa66TQiIj5N5cYPHF4tXOtNGRTeEM56wN7+/t+Qn2U2j4iID1O58QPndT6y3tR3a7TelDG9b4C4LlCYDd+NN51GRMRnqdz4geAAJ0NOsdeb0ozFBjmccMHT9vaKd2DdTLN5RER8lMqNn7gmsTmWBfM37mVjZp7pOP6r+alw6u329hcjIC/TbB4RER+kcuMnmkWHcXb7xgC8u0ijN0adPR4ad4KCLLvgaIJFERGvUrnxI9cfWm/qk2U7KCjW7cjGBIbA314FZzBsmAVLXjWdSETEp6jc+JHTWsfQsmEYuYWlvL1AozdGxXY8sjTDrAcgc63ZPCIiPkTlxo84HBYjzmwNwItzNpF9sMRwIj+XeAu0ToLSQvj0RigtMp1IRMQnqNz4mct7NqNN43pkHyzh5bmbTMfxb5YFl7xorxye8au9uKaIiJwwlRs/43RY3D2oHQCvzd+iJRlMi4iFwc/Z2/Ofg/RVZvOIiPgAlRs/dE7HWHo2r09hiZvnUjaYjiMdBkP7i8BdCl+OBLfLdCIRkTpN5cYPWZbFPee1B+CDJWlsyco3nEi44CkIjoSdS2HJ/0ynERGp01Ru/FTiSQ05s10jXG4Pz8xaZzqORMZD0qElGVIehuwdZvOIiNRhKjd+7O5B7bEs+OqX3azamW06jvS6ARISoTgPZtytyf1ERKpI5caPdYyP5JJu8QBMnKl5VoxzOOyLix2BsG4GrJluOpGISJ2kcuPnks9pR6DT4qcNWfy8Kct0HGncAQaMsrdn3A35e43GERGpi1Ru/FzzhmH8vU9zAJ6cuQ6PToWYd9oYiGkLeRkw7TadnhIRqSSVG+GOs1oTGugkNe0As1dnmI4jgSFwxWuH1p76FhZMNp1IRKROUbkRGkeEcMOAlgA89e06XG6NFBgX1wXOm2Bvf/cQ7FhmNI6ISF2iciMA3Hz6yUSFBrIhM49pK3aajiMAvW+AjpeAuwQ+GQ4HD5hOJCJSJ6jcCABRoYHcOvBkACbNXk9RqWbJNc6yYPB/oX4LOLANvrxL19+IiBwHlRsp849+LWkcEczOAwd5f9F203EEILQ+XPm6fXv46i9gqWYvFhH5Kyo3UiY0yMldZ7cB4PnvN5JfVGo4kQDQtBckPWRvzxwLu1YYjSMiUtup3Eg5Q05JoEXDMPbmF/PavC2m48hhfUdAuwvBVQwfDYOD+00nEhGptVRupJxAp4Pkc9oC8MpPm8kuKDGcSAD7+ptLXzxy/c2023X9jYjIUajcyJ8M7hpPu9gIcgtLefnHTabjyGGh9eGqt+z5b9bNgJ//azqRiEitpHIjf+JwWPzrXHv05vX5W9mTW2Q4kZSJ7w7nP2Fvf/cwbJ1vNI6ISG2kciMVOqdjLN0S6nOwxMXkHzaajiO/12s4dLkKPC745AbIyzSdSESkVlG5kQpZlsWYQ6M37y3azs4DBw0nkjKWBRf9Bxq1h7x0+Hg4uHRnm4jIYSo3clQDWsdw6kkNKHa5eT5lg+k48nvB9eCqtyEoArbNg+/Gm04kIlJrqNzIUVmWxd2D2gHw8bIdbMnKN5xIymnU1r6DCmDBC7DqU7N5RERqiVpRbiZPnkzLli0JCQkhMTGRxYsXH/XYqVOnctpppxEdHU10dDRJSUnHPF5OTK8WDTizXSNcbg//mb3edBz5o44XQ/9R9vYXd0LmGqNxRERqA+Pl5sMPPyQ5OZnx48ezfPlyunXrxqBBg8jMrPgiyTlz5vD3v/+dH374gQULFpCQkMC5557Lzp1a7LG6/Otce/Rm+spdzN+YZTiN/MlZ46DVQCjJhw+uhcJs04lERIyyPB6zM4ElJiZyyimn8MILLwDgdrtJSEjgzjvv5N577/3L57tcLqKjo3nhhRcYOnToXx6fk5NDVFQU2dnZREZGnnB+fzH2s195f/F24iJDmDnqNOqHBZmOJL+XnwUvD4ScHdBmEFz9LjgDTacSEfGayvz9NjpyU1xczLJly0hKSirb53A4SEpKYsGCBcf1GgUFBZSUlNCgQYMKf15UVEROTk65h1TeuIs60ComnPScQu7/fBWGO7H8UXgMDDk0wd+Gb+1bxF2aXVpE/JPRcpOVlYXL5SI2Nrbc/tjYWNLT04/rNe655x7i4+PLFaTfmzBhAlFRUWWPhISEE87tj8KCAnh2SHcCHBZf/7qbz5brNGCt07QXXP0eOINgzXT49EbdIi4ifsn4NTcn4oknnuCDDz7g888/JyQkpMJjxo4dS3Z2dtkjLS2thlP6jm4J9RmVZK8aPn76b6TtKzCcSP6kTRIMeQccgbB6Gnx+swqOiPgdo+UmJiYGp9NJRkZGuf0ZGRnExcUd87lPP/00TzzxBLNmzaJr165HPS44OJjIyMhyD6m6285oTe8W0eQVlTL6w1RKXW7TkeSP2g6y16ByBNq3h0+7Ddwu06lERGqM0XITFBREr169SElJKdvndrtJSUmhb9++R33ek08+yb///W9mzpxJ7969ayKqHOJ0WPxnSHfqBQewdNt+XpqjhTVrpfYXwJWvgyMAfv0IZtytVcRFxG8YPy2VnJzM1KlTefPNN1mzZg233XYb+fn5DB8+HIChQ4cyduzYsuMnTpzIuHHjeO2112jZsiXp6emkp6eTl5dn6iP4nYQGYTx8cScAnk3ZwPLt+w0nkgp1GAx/exWwYOn/4OfnTScSEakRxsvNkCFDePrpp3nwwQfp3r07qampzJw5s+wi4+3bt7N79+6y41966SWKi4u54ooraNKkSdnj6aefNvUR/NLlPZsyuFs8LreHUR+kkluoO3NqpU6XwaDH7e3Z4+C3z83mERGpAcbnualpmufGe7IPlnDBcz+x88BBLu/RlElDupuOJBXxeOCbe2Dxy/at4sO+hOaJplOJiFRKnZnnRuq2qNBAnru6Ow4LPluxky9SdXt4rWRZcN4EaHcBuIrg/athr66VEhHfpXIjJ6R3ywbceZZ9e/gDn6/S7eG1lcNpX38T3wMO7oN3r7BnNRYR8UEqN3LC7jyrNb1aRJNbVMoo3R5eewWFwzUfQf3msG8zvHslFOlCfBHxPSo3csICnA6eHdKdiOAAlm3bz8s/bjYdSY6mXmO47jMIbQC7lsPH/9AyDSLic1RuxCsSGoQx/tDt4c99t4H1GbmGE8lRxbSBaz+GgFDYOBum36U5cETEp6jciNf8rWdTzmrfmGKXm7s/XqnTU7VZs95w5RtgOWHle5DyiOlEIiJeo3IjXmNZFo9f1oWIkABW7sjmlZ90eqpWa3ceDH7O3p43CRa8aDaPiIiXqNyIV8VFhfDgRR0BeHb2Bjbo9FTt1vN6OOsBe/vbsbB4qtk8IiJeoHIjXndFr2ac2a4RxS43Y3R6qvY7bQz0H2VvzxijgiMidZ7KjXidZVlMuLyrTk/VFZYFSQ9B/5H29zPGwJJXjUYSETkRKjdSLeKiQhh36PTUpFnrWbZtn+FEckyWBUkPQ7877e+//hcsfc1sJhGRKlK5kWpzZa9mXNi1CaVuD7e/u5w9uUWmI8mxWBac82/oe4f9/VejYfZ4cJWazSUiUkkqN1JtLMti4t+6cnKjcDJyirjr/RW6/qa2syw499Ej1+DMfxbevAhydplMJSJSKSo3Uq3qBQfw8vW9CA9ysmDzXp6etd50JPkrlgXnPGzPgxMUAdsXwJQBsDHFdDIRkeOiciPVrnXjCCZe0RWAKXM38e1v6YYTyXHpdBncMhfiukDBXnjnb/D9Y+B2mU4mInJMKjdSIy7qGs8N/VsBMOajlWzJyjecSI5Lw5Phn99B7xsAD/z45KEVxfeaTiYiclQqN1Jjxl7QnlNa2quH3/TWUnIKtWBjnRAYAhf9By6fCoFhsOl7eGUg7FxmOpmISIVUbqTGBDodTL62J3GRIWzMzGPUB6m43Fqwsc7oehXc+B00OAmy0+C18+zbxbXopojUMio3UqMaR4TwytBeBAc4+H5tJk99u850JKmM2E5w8xxofxG4iu3bxWeM0XU4IlKrqNxIjevarD5P/u4C4y9SdxpOJJUSEgVD3rFnNcayZzP+aCiUHDSdTEQEULkRQy7p3pTbzzgZgP/75BdWph0wG0gqx7JgwGi48nVwBsHar+CtS6FAM1GLiHkqN2LMmHPbkdShMUWlbm56aymb9uSZjiSV1ekyuP5zCI6CtIX2dTgH0kynEhE/p3IjxjgcFv8Z0p32cRFk5hYx5OWFrM/INR1LKqvlALhhJkQ2hax18OrZsH2h6VQi4sdUbsSoiJBA3r0xkQ5NIsnKK+LqVxby265s07GksmI7wj9nQeOOkJcBb1wIC1/SnVQiYoTKjRjXsF4w79+USLdmUezLL+bvrywkVdfg1D1RzeCfs6HT5eAuhZn3wqf/hCKdbhSRmqVyI7VC/bAg3r4xkV4toskpLOW6VxexeIsuTq1zguvBFa/BeU+AIwBWfWqfpspYbTqZiPgRlRupNSJDAnnrhj6celID8opKuf5/i5ildajqHsuCU2+DYV9BvTjYs9ZeePPrMbqbSkRqhMqN1CrhwQG8/o8+nN3evovq1neW8d6i7aZjSVW06Au3/GhP+OdxwZKp8N/u9rU4Li29ISLVR+VGap3QICcvX9+Lq3o3w+2B+z7/lWe/W49HF6fWPRGxcPW7MHQ6xHaGwmz7WpyX+sH2RabTiYiPUrmRWinA6WDi37py51mtAXj2uw3cP20VpS634WRSJScNtEdxBj8HYTGQtR5ePw++exhKi0ynExEfo3IjtZZlWfzr3Hb8+5JOWBa8t2g7/3h9CQcKik1Hk6pwOKHXP+DOZdDt7+Bxw7xJMPUsSF9lOp2I+BCVG6n1ru/bkpeu7UVooJN5G7O4ZPJ8TfZXl4XWh8umwFVvQ1hDyFgFr5wBc5+EkkLT6UTEB6jcSJ1wXuc4Pru9H82iQ9m2t4DLJs9n9uoM07HkRHS8GG5fCO0uAHcJ/PAYTO4Da77U5H8ickJUbqTO6NAkkul3DODUkxqQX+zi5reXMmnWOopLdR1OnVWvMVz9Hlz+KkQ0gQPb4MPr4K2LIeM30+lEpI6yPH52C0pOTg5RUVFkZ2cTGRlpOo5UQYnLzSNfrubthdsAaB8XwcS/daVbQn2zweTEFOXBvP/Az8+DqwgsB3S9Gk4fAw1PNp1ORAyrzN9vlRups75I3cnDX65mX34xDgtuPO0kRie1JTTIaTqanIj9W2HWA/bpKbBLTper7JIT08ZoNBExR+XmGFRufMvevCIe+Wo1X6TuAqBFwzAev6wL/VvHGE4mJ2zHUpg7ETbMsr+3HND5bzBgNMR2MptNRGqcys0xqNz4ppQ1Gdz/+SrSc+y7bS7v0ZT7L+xAw3rBhpPJCdu53L6Tav03R/a1PQ8GJEPzRHO5RKRGqdwcg8qN78otLOHpb9fx1sJteDxQPyyQ+87vwJW9m2FZlul4cqJ2pdrX5Kz+Ajj0n60W/aHfndDmXHseHRHxWSo3x6By4/tS0w4w9rNfWbM7B4DEVg146opuNG8YZjiZeEXWBpj/LKz80L6FHKB+c+j9T+g5FMIaGI0nItVD5eYYVG78Q6nLzevztzJp9noOlrgIC3Jy3wUduDaxuUZxfEX2Tlj0Eix/GwoP2PucwdDlCuh7B8R2NBpPRLxL5eYYVG78S9q+AsZ8vJJFW/YBcFqbGCb+rSvx9UMNJxOvKS6AVZ/C4lcg/Zcj+9ueB/1H2auTi0idp3JzDCo3/sft9vDGz1uZOHMtRaVuIkICuOPM1lzdpzlRoYGm44m3eDywY4k9T86aLym7LichERJvhXbnQ6BKrUhdpXJzDCo3/mvTnjz+9dFKUtMOABAe5OTK3gkM79+SFg3DzYYT78raCD//F1a+D65DC60G1YP2F0GXK+GkM8AZYDSiiFSOys0xqNz4t1KXm8+W7+TVeZtZn5EHgGXBuR1jGXl2WzrG658Jn5KbDounwi8fQfb2I/vDYqDrVdDjOs2ZI1JHqNwcg8qNAHg8HuZtzOLVn7Ywd/0ewC45g7vGk3xOW1rGaCTHp3g8kLYYfv0YfvsMCvYe+VmT7nbJ6XIFhEYbiygix6ZycwwqN/JHGzJyeS5lA1/9shsAp8Piqt4J3HV2a5pE6RoNn+MqgU3fw4p3YN03R24ndwbbK5X3HAotBoBD6wqL1CYqN8egciNHs2pnNk/PWsecdfZITqDT4tLuTbll4Em0bhxhOJ1Ui/ws+5TVincg83erkEe3hB7XQ7erIaqZsXgicoTKzTGo3MhfWbxlH0/PWsfiQ7ePA5zTMZZbB55Ez+bRmifHF3k8sGu5PWfOr59Ace6RnyWcaq9p1elSqNfYWEQRf6dycwwqN3K8lm3bz5S5m5i9OqNsX0KDUM5uH8tZ7RuTeFIDggM05b/PKc63l3hY8S5sm0/ZLeWWw17u4eSz4KSB9rU6WvJBpMao3ByDyo1U1sbMPF75cRPTUndRXOou2x8W5OT0No04r3McZ7ZvrDlzfFHOLvhtmj1J4M6l5X8WHAUt+0OrgXDymRDT1r4qXUSqhcrNMajcSFUVFJcyf+Nevl+bQcqaTDJzi8p+Fui06HdyDIM6xXFamxiaRYfq9JWv2bcFNsyCLT/C1p+gMLv8zyOa2PPnnHQGND8V6rdQ2RHxIpWbY1C5EW/weDys2pnD7NXpfLMqnQ2ZeeV+3igimN4tounVIpqeLaLpEBdJaJBOYfgMt8te6mHzXNg8B7YvgNLC8seEN4Zmp0Cz3vYjtrMW9RQ5ASo3x6ByI9VhY2Ye3/6WzndrMli1M5sSV/l/rRwWnNSoHp3iI+kUH0nvlg3o1qw+Tof+n71PKCmEtIV20dnyI+z+5cgt5r8X2dSeNDC2M8R1gSbdILqVbjsXOQ4qN8egciPVrbDExS87slm2bT/Ltu0jNe0AWXnFfzouOiyQgW0bcWb7xpzephHR4UEG0kq1KCmE3Sthx2J78sDdK+HAtoqPDap3pOzEtIWY1tCwjV2EVHpEyqjcHIPKjZiQmVPIb7ty+G1XNr/uzGbBpr3kFJaW/dyyoFVMOJ3jo+jSNIrOTaPo2CSSqDBdpOwzCnMg4zfIWAXpv9qPzNV/Pp11WEAoNDgJGrSyH9Gt7O8btlbxEb+kcnMMKjdSG5S63CzffoAf1mXyw9pM1qbnVnhcw/AgTmoUzkkx9TipUTgtGobRLDqMhAZhujvLF7hKYe8G+zRW5m/2gp97N8C+zeAuPfrzAkLtkhPT2h7tadQOGnWw9wVoBFB8k8rNMajcSG20N6+IVbtyWLUzm1932KM7Ow8cPOZzIkMCaN4wjIRDZSehQRgJ0aE0iw6lSVQo4cFa9brOcpXA/m12ydm/FfZvse/W2rf5UPGp4HoeAEcANDjZnlW5XmP7Ed4YIuLsffWbQ704jfpInaRycwwqN1JX5BeVsiUrn0178ti8J5/NWfls31fAjn0F7M3/8zU8fxQVGkiTqBDi64cSGxlMo4gQGkcEExsZQsN6QdQPDaR+WBCRIQEEOPXHrs5wldrX7+zdCFkbYM9a+5G5tvzMykfjCISophARD+ENIbyRvUp6eAwEhUNACASGQWAohETaZaheY3BqpFDMqnPlZvLkyTz11FOkp6fTrVs3nn/+efr06XPU4z/++GPGjRvH1q1badOmDRMnTuSCCy44rvdSuRFfkF9Uyo79B0nbV0Da/gK27ysgbd9BduwvYOeBg+QWHuOURgUiggOoHx5Ig7Ag6ocF0SA8iPphh74PD6JBWBDRYYE0iggmNiqEiOAAzeNT23g8kLMT9qyD3HTIz4S8TMjLsL8/kGb/3OOqwotbENbQnssnItYeCaoXd+hrrP2zsAb219BoFSGpFnWq3Hz44YcMHTqUKVOmkJiYyLPPPsvHH3/MunXraNz4z+u4/Pzzz5x++ulMmDCBiy66iPfee4+JEyeyfPlyOnfu/Jfvp3Ij/iC3sITd2YXsOnCQ3dmFZOYUkZFrf83MLWRvXjE5B0vILapcCTosLMhJXJQ9ElQvOICQQCchgU5CA52EBwcQFRpI/bBA++uhEaLo8ECiw4IICdR8P8a4SiF3N2Sn2aUnP+vQYw8U7IWSg1BSYH8tLYSD++3jjnX9T0WCIyGkPoQefkQf+j76yPfBkYdGiH43UhQQAgHBR746g8By2ktfWA77ynuVar9Vp8pNYmIip5xyCi+88AIAbrebhIQE7rzzTu69994/HT9kyBDy8/P56quvyvadeuqpdO/enSlTpvzl+6nciBxR4nKTc7CEAwdLOFBQwv78YvYXHH6UcKCgmH35xezPL2FfQTGZOYXl7vKqitBAJ1GhgYQFHSpEQU5CAh0EOR04HQ4CHBZOp4XTsnBY4LAsrEPbAU4HoYFOQoPsryGBTpwOC8ehn9vH/f55R77+kfWH17b483EWlP3cYVk4HBDodBx6WAQ6HTi8+Mf28Hsc/gx2ArMcuAko3E9gQSYBBek4CzJx5mfgzM/EmZ+Oo2APjoP7cRTuw1F4AIua/5PicQTicQaCIwiPMxCPIxAcTjxWgL3+l8OJp6wkWYADj8MJzmA8ziA8ASF4nMF4HAFHihQWWIeOcwTYX60A+7Wdh78G2e/tCLCvd/rde3rKCtnvXq+i7H88xmv/PFl4LKt8Bsth/6/zp32/L45//Jfg0O+r7LUqOKYCASFhxDRt46XPYqvM32+jVxwWFxezbNkyxo4dW7bP4XCQlJTEggULKnzOggULSE5OLrdv0KBBTJs2rcLji4qKKCo6Mk1+Tk7OiQcX8RGBTgcN6wXTsF7wcT+noLiUjJwi0rMLycwtpLDExcFiFwdL3BwscZFXWEr2wRKyDxaTfbDkUEkqYX9BMS63h4MlLg6WVOXUiJjnAOIOPSr6qZso8qhv5RNFPvWtPCLJJ9rKI4p8oqwj+yKtAoIpJpRiQigm1CommGKCKSHEOsoF00dhuUuw3CVAwQl/QvGOtQEdiHlgobH3N1pusrKycLlcxMbGltsfGxvL2rVrK3xOenp6hcenp6dXePyECRN4+OGHvRNYRAgLCqBVTACtYsIr9TyPx0NuUSn78+3SU3ioDB0sdlFY4qLE5cbl9lDq9pR99Xg8uD0ePB5we6C41E1hqf2cokNfXR4OHePB7QbXoeMPP9dVwUDC4QFrz6Hn2o+KQpf/ucvtocTlPvSwt91eGvz2eOz1xz2/ey/TPB7PoUz27+GvMzkooD4F1GfXUV/0D1+PclAgLoIowcKDAzeOQ1//yMKDEzdBlBL4u4cDNwG4cOImwHKVex0LDwG4CKKUYIoJooRgSnDiOvQ+9rGHX9tpuQjAjRMXgZQScOhrIKUEUYrz8PvgKntfq+x1PDisP+c+nN3O4y471lusQ5/TUfa53dhjMO5yn89Rdpy7wve3x2mO5LMq+N+gIkWOUK99lqrw+XtFx44dW26kJycnh4SEBIOJRPyTZVlEhgQSGaKLTUV8XX3D72+03MTExOB0OsnIyCi3PyMjg7i4ioc94+LiKnV8cHAwwcHHP+QuIiIidZvRyS2CgoLo1asXKSkpZfvcbjcpKSn07du3wuf07du33PEAs2fPPurxIiIi4l+Mn5ZKTk5m2LBh9O7dmz59+vDss8+Sn5/P8OHDARg6dChNmzZlwoQJAIwcOZKBAwfyzDPPcOGFF/LBBx+wdOlSXnnlFZMfQ0RERGoJ4+VmyJAh7NmzhwcffJD09HS6d+/OzJkzyy4a3r59O47fTRXer18/3nvvPR544AHuu+8+2rRpw7Rp045rjhsRERHxfcbnualpmudGRESk7qnM328tKCMiIiI+ReVGREREfIrKjYiIiPgUlRsRERHxKSo3IiIi4lNUbkRERMSnqNyIiIiIT1G5EREREZ+iciMiIiI+xfjyCzXt8ITMOTk5hpOIiIjI8Tr8d/t4Flbwu3KTm5sLQEJCguEkIiIiUlm5ublERUUd8xi/W1vK7Xaza9cuIiIisCyryq+Tk5NDQkICaWlpWqOqBuj3XbP0+65Z+n3XLP2+a5a3ft8ej4fc3Fzi4+PLLahdEb8buXE4HDRr1sxrrxcZGal/OWqQft81S7/vmqXfd83S77tmeeP3/VcjNofpgmIRERHxKSo3IiIi4lNUbqooODiY8ePHExwcbDqKX9Dvu2bp912z9PuuWfp91ywTv2+/u6BYREREfJtGbkRERMSnqNyIiIiIT1G5EREREZ+iciMiIiI+ReWmiiZPnkzLli0JCQkhMTGRxYsXm47kkyZMmMApp5xCREQEjRs35tJLL2XdunWmY/mNJ554AsuyGDVqlOkoPmvnzp1cd911NGzYkNDQULp06cLSpUtNx/JJLpeLcePG0apVK0JDQzn55JP597//fVxrFclf+/HHHxk8eDDx8fFYlsW0adPK/dzj8fDggw/SpEkTQkNDSUpKYsOGDdWSReWmCj788EOSk5MZP348y5cvp1u3bgwaNIjMzEzT0XzO3LlzGTFiBAsXLmT27NmUlJRw7rnnkp+fbzqaz1uyZAkvv/wyXbt2NR3FZ+3fv5/+/fsTGBjIN998w+rVq3nmmWeIjo42Hc0nTZw4kZdeeokXXniBNWvWMHHiRJ588kmef/5509F8Qn5+Pt26dWPy5MkV/vzJJ5/kv//9L1OmTGHRokWEh4czaNAgCgsLvR/GI5XWp08fz4gRI8q+d7lcnvj4eM+ECRMMpvIPmZmZHsAzd+5c01F8Wm5urqdNmzae2bNnewYOHOgZOXKk6Ug+6Z577vEMGDDAdAy/ceGFF3puuOGGcvsuv/xyz7XXXmsoke8CPJ9//nnZ92632xMXF+d56qmnyvYdOHDAExwc7Hn//fe9/v4auamk4uJili1bRlJSUtk+h8NBUlISCxYsMJjMP2RnZwPQoEEDw0l824gRI7jwwgvL/XMu3jd9+nR69+7NlVdeSePGjenRowdTp041Hctn9evXj5SUFNavXw/AypUrmTdvHueff77hZL5vy5YtpKenl/tvSlRUFImJidXyt9PvFs48UVlZWbhcLmJjY8vtj42NZe3atYZS+Qe3282oUaPo378/nTt3Nh3HZ33wwQcsX76cJUuWmI7i8zZv3sxLL71EcnIy9913H0uWLOGuu+4iKCiIYcOGmY7nc+69915ycnJo3749TqcTl8vFY489xrXXXms6ms9LT08HqPBv5+GfeZPKjdQZI0aMYNWqVcybN890FJ+VlpbGyJEjmT17NiEhIabj+Dy3203v3r15/PHHAejRowerVq1iypQpKjfV4KOPPuLdd9/lvffeo1OnTqSmpjJq1Cji4+P1+/YxOi1VSTExMTidTjIyMsrtz8jIIC4uzlAq33fHHXfw1Vdf8cMPP9CsWTPTcXzWsmXLyMzMpGfPngQEBBAQEMDcuXP573//S0BAAC6Xy3REn9KkSRM6duxYbl+HDh3Yvn27oUS+7e677+bee+/l6quvpkuXLlx//fWMHj2aCRMmmI7m8w7/faypv50qN5UUFBREr169SElJKdvndrtJSUmhb9++BpP5Jo/Hwx133MHnn3/O999/T6tWrUxH8mlnn302v/76K6mpqWWP3r17c+2115KamorT6TQd0af079//T1MbrF+/nhYtWhhK5NsKCgpwOMr/2XM6nbjdbkOJ/EerVq2Ii4sr97czJyeHRYsWVcvfTp2WqoLk5GSGDRtG79696dOnD88++yz5+fkMHz7cdDSfM2LECN577z2++OILIiIiys7NRkVFERoaajid74mIiPjT9Uzh4eE0bNhQ1zlVg9GjR9OvXz8ef/xxrrrqKhYvXswrr7zCK6+8YjqaTxo8eDCPPfYYzZs3p1OnTqxYsYJJkyZxww03mI7mE/Ly8ti4cWPZ91u2bCE1NZUGDRrQvHlzRo0axaOPPkqbNm1o1aoV48aNIz4+nksvvdT7Ybx+/5WfeP755z3Nmzf3BAUFefr06eNZuHCh6Ug+Cajw8frrr5uO5jd0K3j1+vLLLz2dO3f2BAcHe9q3b+955ZVXTEfyWTk5OZ6RI0d6mjdv7gkJCfGcdNJJnvvvv99TVFRkOppP+OGHHyr87/WwYcM8Ho99O/i4ceM8sbGxnuDgYM/ZZ5/tWbduXbVksTweTc0oIiIivkPX3IiIiIhPUbkRERERn6JyIyIiIj5F5UZERER8isqNiIiI+BSVGxEREfEpKjciIiLiU1RuRKTazZkzB8uyOHDgQI2+7xtvvEH9+vVP6DW2bt2KZVmkpqYe9RhTn09EKqZyIyJed8YZZzBq1CjTMUTET6nciEitVFxcbDqCiNRRKjci4lX/+Mc/mDt3Ls899xyWZWFZFlu3bgVg2bJl9O7dm7CwMPr161duReyHHnqI7t278+qrr9KqVStCQkIAOHDgADfeeCONGjUiMjKSs846i5UrV5Y9b+XKlZx55plEREQQGRlJr169WLp0ablM3377LR06dKBevXqcd9557N69u+xnbrebRx55hGbNmhEcHEz37t2ZOXPmMT/jjBkzaNu2LaGhoZx55plln09EageVGxHxqueee46+ffty0003sXv3bnbv3k1CQgIA999/P8888wxLly4lICDgT6sxb9y4kU8//ZTPPvus7BqXK6+8kszMTL755huWLVtGz549Ofvss9m3bx8A1157Lc2aNWPJkiUsW7aMe++9l8DAwLLXLCgo4Omnn+btt9/mxx9/ZPv27YwZM6Zc3meeeYann36aX375hUGDBnHxxRezYcOGCj9fWloal19+OYMHDyY1NZUbb7yRe++915u/QhE5UdWyHKeI+LU/riR+eLXg7777rmzf119/7QE8Bw8e9Hg8Hs/48eM9gYGBnszMzLJjfvrpJ09kZKSnsLCw3OuffPLJnpdfftnj8Xg8ERERnjfeeKPCHK+//roH8GzcuLFs3+TJkz2xsbFl38fHx3see+yxcs875ZRTPLfffrvH4/F4tmzZ4gE8K1as8Hg8Hs/YsWM9HTt2LHf8Pffc4wE8+/fvP9avRURqiEZuRKTGdO3atWy7SZMmAGRmZpbta9GiBY0aNSr7fuXKleTl5dGwYUPq1atX9tiyZQubNm0CIDk5mRtvvJGkpCSeeOKJsv2HhYWFcfLJJ5d738PvmZOTw65du+jfv3+55/Tv3581a9ZU+BnWrFlDYmJiuX19+/Y97t+BiFS/ANMBRMR//P50kWVZgH3Ny2Hh4eHljs/Ly6NJkybMmTPnT691+Bbvhx56iGuuuYavv/6ab775hvHjx/PBBx9w2WWX/ek9D7+vx+PxxscRkVpKIzci4nVBQUG4XK4Tfp2ePXuSnp5OQEAArVu3LveIiYkpO65t27aMHj2aWbNmcfnll/P6668f1+tHRkYSHx/P/Pnzy+2fP38+HTt2rPA5HTp0YPHixeX2LVy4sJKfTESqk8qNiHhdy5YtWbRoEVu3biUrK6vc6ExlJCUl0bdvXy699FJmzZrF1q1b+fnnn7n//vtZunQpBw8e5I477mDOnDls27aN+fPns2TJEjp06HDc73H33XczceJEPvzwQ9atW8e9995LamoqI0eOrPD4W2+9lQ0bNnD33Xezbt063nvvPd54440qfT4RqR4qNyLidWPGjMHpdNKxY0caNWrE9u3bq/Q6lmUxY8YMTj/9dIYPH07btm25+uqr2bZtG7GxsTidTvbu3cvQoUNp27YtV111Feeffz4PP/zwcb/HXXfdRXJyMv/617/o0qULM2fOZPr06bRp06bC45s3b86nn37KtGnT6NatG1OmTOHxxx+v0ucTkepheXTyWURERHyIRm5ERETEp6jciIiIiE9RuRERERGfonIjIiIiPkXlRkRERHyKyo2IiIj4FJUbERER8SkqNyIiIuJTVG5ERETEp6jciIiIiE9RuRERERGfonIjIiIiPuX/AUYL8ey48thMAAAAAElFTkSuQmCC", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "flip_flop_index_proportion_exceeding(fcst1, \"lead_day\", np.arange(0.1, 10, 0.1), days123=[1, 2, 3]).days123.plot()\n", + "flip_flop_index_proportion_exceeding(fcst2, \"lead_day\", np.arange(0.1, 10, 0.1), days123=[1, 2, 3]).days123.plot()" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "b5f8158f", + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "hide_input": false, + "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.12.0" + }, + "toc": { + "base_numbering": 1, + "nav_menu": {}, + "number_sections": true, + "sideBar": true, + "skip_h1_title": false, + "title_cell": "Table of Contents", + "title_sidebar": "Contents", + "toc_cell": false, + "toc_position": {}, + "toc_section_display": true, + "toc_window_display": false + }, + "varInspector": { + "cols": { + "lenName": 16, + "lenType": 16, + "lenVar": 40 + }, + "kernels_config": { + "python": { + "delete_cmd_postfix": "", + "delete_cmd_prefix": "del ", + "library": "var_list.py", + "varRefreshCmd": "print(var_dic_list())" + }, + "r": { + "delete_cmd_postfix": ") ", + "delete_cmd_prefix": "rm(", + "library": "var_list.r", + "varRefreshCmd": "cat(var_dic_list()) " + } + }, + "types_to_exclude": [ + "module", + "function", + "builtin_function_or_method", + "instance", + "_Feature" + ], + "window_display": false + } + }, + "nbformat": 4, + "nbformat_minor": 5 +}