{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "d5cf2b3d-0f34-4d13-9ba3-dd2e3fb66337",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Figure 3: Plot 2-D histograms of WD and WDA for 96839 and 96253\n",
    "\n",
    "from docx import Document\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import matplotlib.ticker as mticker\n",
    "from itertools import chain\n",
    "from matplotlib.colors import LinearSegmentedColormap, Normalize\n",
    "import matplotlib.colors as mcolors\n",
    "from matplotlib.colors import BoundaryNorm\n",
    "\n",
    "# Load the Word document\n",
    "file_path = \"/Users/cdz/Desktop/Khafid_etal/BMKG_SurfaceStationData/CombinedStations.docx\"\n",
    "doc = Document(file_path)\n",
    "\n",
    "len(doc.tables)\n",
    "\n",
    "table = doc.tables[0]\n",
    "\n",
    "data = []\n",
    "for row in table.rows:\n",
    "    data.append([cell.text.strip() for cell in row.cells])\n",
    "\n",
    "df = pd.DataFrame(data)\n",
    "df.head()\n",
    "\n",
    "# Define column names\n",
    "df.columns = [\n",
    "    \"station_id\",\n",
    "    \"station_name\",\n",
    "    \"region\",\n",
    "    \"longitude\",\n",
    "    \"latitude\",\n",
    "    \"elevation_m\",\n",
    "    \"value\"\n",
    "]\n",
    "\n",
    "id = pd.to_numeric(df[\"station_id\"], errors=\"coerce\")\n",
    "nid = id.shape[0]\n",
    "\n",
    "y_labs = ['WD (degree)', 'WDA (degree)'] \n",
    "x_labs = \"Local Time (hour)\"\n",
    "labs = ['(a)', '(b)', '(c)', '(d)', '(e)', '(f)', '(g)', '(h)']\n",
    "\n",
    "cmap = LinearSegmentedColormap.from_list(\"white_to_dark\", [\"white\", \"black\"])\n",
    "norm = BoundaryNorm([0, 1, np.inf], ncolors=cmap.N, clip=True)\n",
    "\n",
    "count_std = np.zeros(6) \n",
    "\n",
    "y1 = [0,  0]\n",
    "y2 = [360, 360]\n",
    "\n",
    "vsh = ['WD', 'WDA']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "51ce5ba3-dcbb-4870-a512-b2c338c0e4f1",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1 WD WD (degree)\n",
      "96011\n",
      "96581\n",
      "2 WDA WDA (degree)\n",
      "96011\n",
      "96581\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmMAAAIACAYAAADOqGuJAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAACL50lEQVR4nO3dd3xUVfo/8M8Q0kiDJIT0GFoA6QliIlUgFGmrK7i6igKyfBGUXlQWK6Cu2AFBFgUV1BUQBYWoJJTISl0gsNSEQEgMJYUUUs/vD36ZzZBk5pnkTiaTfN6v17yUOyfnnLlz55kz957nHp1SSoGIiIiIrKKRtTtARERE1JBxMEZERERkRRyMEREREVkRB2NEREREVsTBGBEREZEVcTBGREREZEUcjBERERFZEQdjRERERFbEwRgRERGRFXEwdofr16/Dx8cHSUlJovKzZ8/Gs88+a9lOERFVony8io2NhU6nQ2ZmZpXlf/jhB3Tr1g2lpaW110kiMomDsTssWbIEI0aMwF133SUqP3fuXKxduxaJiYmW7RgR0R3MjVfDhw+HTqfDl19+admOEZFZOBgrJz8/H2vWrMHEiRPFf+Pj44Po6GisXLnSgj0jIjJUnXgFAE899RQ++OADC/WKiKqDg7FyfvzxRzRu3BiRkZEAgJKSEkyYMAGhoaFwdnZGWFgY3nvvvQp/N3LkSGzYsKG2u0tEDdid8arMvn370KVLFzg5OaFnz544fvy4wfMjR47E77//jgsXLtRmd4nICA7Gytm9ezciIiL0/y4tLUVgYCC+/vprnDx5En//+9/x/PPP4+uvvzb4u3vuuQeXLl3CxYsXa7vLRNRA3RmvysyZMwf/+Mc/cODAAfj4+GDkyJEoKirSPx8SEgIfHx/s2bOnNrtLREY0tnYH6pKkpCT4+/vr/21vb4+XX35Z/+/Q0FDEx8fj66+/xpgxY/TbAwIC9H8fEhJSex0mogbrznhVZtGiRRg0aBAA4LPPPkNgYCA2b95cIWZJk5SIyPI4GCsnPz8fTk5OBttWrlyJTz75BBcvXkR+fj4KCwvRtWtXgzLOzs4AgLy8vNrqKhE1cJXFKwAGly09PT0RFhaGU6dOGZRxdnZmvCKqQ3iZshxvb29kZGTo//31119jxowZGD9+PHbu3ImjR4/iqaeeQmFhocHf3bhxAwDQvHnzWu0vETVcd8YrY3Q6ncG/b9y4wXhFVIdwMFZOt27dcPLkSf2/9+zZg6ioKEyZMgXdunVD69atcf78+Qp/d+LECdjb2+Puu++uze4SUQN2Z7wqs3//fv3/Z2Rk4MyZM2jXrp1+261bt3D+/Hl069atVvpJRKZxMFbO4MGDkZCQoP+12bp1axw8eBA7duzAmTNnsHDhQhw4cKDC3+3Zswe9e/fWX64kIrK0O+NVmVdeeQW//PILTpw4gSeffBLe3t4YPXq0/vn9+/fD0dGxQhYmEVkPB2PldOrUCREREfpsycmTJ+PBBx/E2LFj0bNnT1y/fh1Tpkyp8HcbNmzA008/XdvdJaIG7M54VWbp0qV47rnnEB4ejtTUVGzduhUODg765zds2IDHHnsMTZo0qe0uE1EVdEopZe1O1CXbt2/H7NmzceLECTRqZHqsum3bNsyZMwfHjh1D48bMhyCi2mNuvLp69SratWuHgwcPIjQ0tBZ6SEQSHD3cYdiwYTh79ixSUlIQFBRksnxubi7Wrl3LgRgR1Tpz41ViYiKWL1/OgRhRHcMzY0RERERWxDljRERERFZk1cHYihUr0LlzZ7i7u8Pd3R2RkZH48ccfDcqcOnUKI0eOhIeHB9zc3HDvvfciOTlZ/3xBQQGmTZsGb29vuLi4YOTIkbh8+XJtvxQiqucYr4jIUqw6GAsMDMTSpUtx8OBBHDx4EPfffz9GjRqFhIQEAMD58+fRq1cvtGvXDrGxsfjPf/6DhQsXGtx1evr06di8eTM2btyIvXv3IicnB8OHD0dJSYm1XhYR1UOMV0RkKXVuzpinpyfeeustTJgwAY888gjs7e2xfv36SstmZWWhefPmWL9+PcaOHQsAuHLlCoKCgrB9+3YMHjy4NrtORA0M4xURaaHOpACWlJTgm2++QW5uLiIjI1FaWopt27Zh7ty5GDx4MI4cOYLQ0FAsWLBAfwPDQ4cOoaioCNHR0fp6/P390bFjR8THx1cZ3AoKClBQUKD/d2lpKW7cuAEvL68Ky4YQUe1RSuHmzZvw9/cX3arBWhiviEjTeKWs7NixY8rFxUXZ2dkpDw8PtW3bNqWUUqmpqQqAatKkiVq2bJk6cuSIWrJkidLpdCo2NlYppdQXX3yhHBwcKtQ5aNAgNWnSpCrbXLRokQLABx981NHHpUuXLBNwaojxig8++LjzoUW8svqZsbCwMBw9ehSZmZn49ttvMW7cOMTFxaFp06YAgFGjRmHGjBkAgK5duyI+Ph4rV65E3759q6xTKWX0F+OCBQswc+ZM/b+zsrIQHByMS5cuwd3dXZsXRkRmy87ORlBQENzc3KzdlUoxXhFRGS3jldUHYw4ODmjdujUAICIiAgcOHMB7772HDz74AI0bN0aHDh0Myrdv3x579+4FAPj6+qKwsBAZGRlo1qyZvkx6ejqioqKqbNPR0RGOjo4VtpdlSRGRddXVy2+MV0R0Jy3iVZ2blKGUQkFBARwcHNCjRw+cPn3a4PkzZ84gJCQEABAeHg57e3vExMTon09NTcWJEyeMBjciIi0wXhGRFqx6Zuz555/H0KFDERQUhJs3b2Ljxo2IjY3FTz/9BACYM2cOxo4diz59+qB///746aef8P333yM2NhYA4OHhgQkTJmDWrFnw8vKCp6cnZs+ejU6dOmHgwIFWfGVEVN8wXhGRxdR41lkNjB8/XoWEhCgHBwfVvHlzNWDAALVz506DMmvWrFGtW7dWTk5OqkuXLmrLli0Gz+fn56upU6cqT09P5ezsrIYPH66Sk5PN6kdWVpYCoLKysmr8moio+uryZ5HxiojK0/KzWOfuM2YN2dnZ8PDwwI0bN4zOwSgqKhLV5+DgYLKMlnNipHVJbiwpTc/Vqv/Sm11mZWWZLJObmyuqKycnx2QZyXuoNRcXF5NlpPtd8rF2dnY2Wab8DUuNadzY9En24uJik2Wys7Ph4+ODrKwszoeqQlm8SklJMbqP7OzsRPVpeQsRSV3SY1hSThqT8/LyTJYpf/uQmpQBbt9+xBTpZ0vyWZa+h5LPvDSOahkjJe+15DVKYgwA2Nvbmywj+fyUTeDXIl7VuTljRERERA0JB2NEREREVsTBGBEREZEVcTBGREREZEUcjBERERFZkdXvwF+X5OfnG82ykGaPSDN8JCRZJtKsKUmGjzQrR5IFKclsyc/PF7Unyaa884abVcnMzDRZRpLZqDVJhk9AQICoLslxI1nCw9XVVdSeJDNMkhWm5WenvsvMzDT6OdQyw05KEiNrO2MbAG7dumWyzPXr1zUpA8iOY+mNDJo0aWKyjPS7SavPKSD7PpHENGmbkvak34WVrWhRHZLMfCmeGSMiIiKyIg7GiIiIiKyIgzEiIiIiK+JgjIiIiMiKOBgjIiIisiJmU5bj4uJiNItOmmWiZRaQlhmQWvZLkrVS25mgYWFhorokmaDZ2dmiuiTZSdJsNS3fa0kWpKRf0vdHQpJZVVhYqFl79Z2Xl5fR9fCk+1LLDEgtSdY7lXxmAMDHx8dkGcnagi1atBC1J1nDUpppKIlX0vUkJftUGq8k8V2yJigg65dkf0nXOpYc85LYJ30PJXhmjIiIiMiKOBgjIiIisiIOxoiIiIisiIMxIiIiIiviYIyIiIjIiphNWc6tW7fEa3wZI8kykazbKK1LkokiJc2a0irz7+bNm6L2kpKSNKtLkul048YNzeqSroPm4eFhsoz0vZYcx4GBgSbLSNemlGQVSdbrk2Zf0e31Wo19DqWxTMusZy1jkSTGSLPnJPFWcnz+8ccfovYksUi6Lq+kTWnsk+wvydq9gCxeSda/ldbl6empSRlAdmcESXzUMsuYZ8aIiIiIrIiDMSIiIiIr4mCMiIiIyIo4GCMiIiKyIg7GiIiIiKyI2ZTlODg4GM1Akq6DJiFd00qS9SFdc1LSf2nWlKT/kgwmacZXSEiIyTLS9dkka/a1atVKVJckm1K6RqCxdVHLSNbCBGQZnJL2pOuxStazk2RycW1KuWbNmhldT1Ga2SjJCJPGGElGopZr6Upfo+RzI+m7NC5IsoKvXr0qqqtZs2Ymy0izvyVtNm/eXFSX5DUGBQWJ6mratKnJMl5eXibLSNYXBWTHg+T4kxwzUjwzRkRERGRF1TozdunSJSQlJSEvLw/NmzfH3XffLb6XEhFRbWK8IqK6Tnxm7OLFi1iwYAHuuusu3HXXXejbty+GDh2KiIgIeHh4YNCgQfjmm2/MupS3ZMkS9OjRA25ubvDx8cHo0aNx+vRpgzI6na7Sx1tvvaUvU1BQgGnTpsHb2xsuLi4YOXIkLl++LO4HEdUvjFdEZEtEg7HnnnsOnTp1wtmzZ/HKK68gISEBWVlZKCwsRFpaGrZv345evXph4cKF6Ny5Mw4cOCBqPC4uDs888wz279+PmJgYFBcXIzo62mDuT2pqqsHjn//8J3Q6HR566CF9menTp2Pz5s3YuHEj9u7di5ycHAwfPlx8d2Yiqj8Yr4jI1uiUYIbunDlzMHfuXNHEvu3btyMvLw9//vOfze7M1atX4ePjg7i4OPTp06fSMqNHj8bNmzfxyy+/ALi9JEjz5s2xfv16jB07FgBw5coVBAUFYfv27Rg8eLDJdrOzs+Hh4YEbN24YnQCo5QR+6STWujqBX9J/yQT+W7duidqTTJTXcgK/dJ9yAv9tWk3gz87ORvPmzZGVlSWejHunhhKv0tPTG/wEfmm/tJqMLZ2wreUEfsmySVpO4JceN5LXGBwcLKrLVifwZ2dnw8fHp0bxqoxor5c/xW7KsGHDqt2ZrKwsAFWvL/XHH39g27Zt+Oyzz/TbDh06hKKiIkRHR+u3+fv7o2PHjoiPj680uBUUFBh8iWZnZwMAcnJyjAYKLddtlH7JSbINpb+oJQeXNMtTMmiTBC7pmmppaWkmy+Tk5Ijqkuwv6fsjGYxJg5tkDTrpWm+SwKXljwvJgFOyH6SDc2MaSrwqLCw0ut+lx7BW68xK29RykCh9jZJykkHPtWvXRO1JPssZGRmiuv773/+aLHPhwgVRXefOnTNZRhJrAdn706JFC1FdDzzwgMkyPXv2NFlGOjiXvNeS70Lpd45EtbIpi4uL8fPPP+Pjjz/Wf5leuXKlRh1TSmHmzJno1asXOnbsWGmZzz77DG5ubnjwwQf129LS0uDg4FAh/bdFixZVHlRLliyBh4eH/iFNvyUi28N4RUR1ndnZlBcvXsSQIUOQnJyMgoICDBo0CG5ubnjzzTdx69YtrFy5slodmTp1Ko4dO4a9e/dWWeaf//wnHnvsMdEpRqVUlaPkBQsWYObMmfp/Z2dnM8AR1UOMV0RkC8w+M/bcc88hIiICGRkZBvNE/vSnP+nnRZhr2rRp2Lp1K3bt2oXAwMBKy+zZswenT5/GxIkTDbb7+vqisLCwwinf9PT0Kk+ROjo6wt3d3eBBRPUP4xUR2QKzB2N79+7Fiy++WGEuU0hICFJSUsyqSymFqVOnYtOmTfj1118RGhpaZdk1a9YgPDwcXbp0MdgeHh4Oe3t7xMTE6LelpqbixIkTiIqKMqs/RFS/MF4RkS0w+zJlaWlppROgL1++LJ5cXOaZZ57Bl19+ie+++w5ubm76ORMeHh4Gv2Kzs7PxzTff4O23365Qh4eHByZMmIBZs2bBy8sLnp6emD17Njp16oSBAwea+eqIqD5hvCIiW2D2YGzQoEF49913sWrVKgC3sxdycnKwaNEiszOTVqxYAQDo16+fwfa1a9fiySef1P9748aNUErhL3/5S6X1vPPOO2jcuDHGjBmD/Px8DBgwAJ9++qn4Ng1lXF1d4erqWuXz0rRmSbvSLKDazoCUkvRf0i/pF2KTJk1MlpFmZkpI18zU6hYfUtL3UHJrC0lWm3Q/aHWPLK3vjF+f45Wbm5vRz4/0PdEyLmhJEmOkGcGSOCq51Yu0PUldZdm4pkjW5U1MTBTVJbkdhTSOSrJBpbe2kGRdVjUloDzJ9wSgXeyTZm9KiO4zVt6VK1fQv39/2NnZ4ezZs4iIiMDZs2fh7e2N3bt3w8fHR7PO1Rbpfcbq6mBMqrYHY5Iy0lsZSL5YOBj7H1sdjGl53x6gfscrU/uIgzHtSWOM5JY30sGY5HunIQzG2rZta7JMbQ/GsrOz4e3tXXv3GSvP398fR48exYYNG3D48GGUlpZiwoQJeOyxx0Q3fiQiqi2MV0RkC6q1ULizszPGjx+P8ePHa90fIiJNMV4RUV1XrZu+rl+/Hr169YK/vz8uXrwI4PY8iO+++07TzhER1RTjFRHVdWYPxlasWIGZM2di6NChyMjI0M9LaNasGd59912t+0dEVG2MV0RkC8y+TPnBBx9g9erVGD16NJYuXarfHhERgdmzZ2vaudpWVFRkdLKkdE01ycRZ6Vpvkrq0XChc2i+tSBMZJBP9pftBkuUpTdaQ7C/pBGnJvpBOqJfUJTm2pJPAJeUk+13r468+xytTa1NKPw9aJqFIjnUtk5K0jKOS1yj9LEteo/SznJmZabKMdP6jZHK+9HsuICDAZBljdycwl+T9kcYrSYKTVklqUmZHvsTERHTr1q3CdkdHR+Tm5mrSKSIiLTBeEZEtMHswFhoaiqNHj1bY/uOPP6JDhw5a9ImISBOMV0RkC8y+TDlnzhw888wzuHXrFpRS+P3337FhwwYsWbIEn3zyiSX6SERULYxXRGQLzB6MPfXUUyguLsbcuXORl5eHRx99FAEBAXjvvffwyCOPWKKPRETVwnhFRLbArMFYcXExvvjiC4wYMQJPP/00rl27htLSUpu8izUR1W+MV0RkK8wajDVu3Bj/93//h1OnTgEAvL29LdIpaykoKDC6jIV06R5Jloy0LglpdpIkG0WanSTJuJEsvZGdnS1qLycnx2QZ6RqdTk5OJstIljMBaj8zTLrchyRzVpJRJNnvgHbLi0iOGan6Hq+Ki4uNfqalx6aW2d+S90+6hJEkQ1CazSbJgszPzzdZ5vr166L2JHVJlhMCZEsdnTt3TlTXjRs3TJaRfjd5eXmJyklIlkMyljlcJjQ0VNSe5P2RLG8k/Z6QMHsCf8+ePXHkyBHNOkBEZCmMV0RkC8yeMzZlyhTMmjULly9fRnh4eIXV6Tt37qxZ54iIaoLxiohsgdmDsbFjxwIAnn32Wf02nU4HpRR0Op34pmtERJbGeEVEtsDswZjk+jURUV3AeEVEtsDswVhISIgl+kFEpDnGKyKyBWYPxrZu3Vrpdp1OBycnJ7Ru3Vqc0VDXNGnSxGi2mnS9QkmWmqQMIF9rS0LLbEpJRpQkG0WaHSjJfpGu9SbJ/NNyvU8ts9qkr1FyrEqy1aRrEkpoufamVH2OV/b29uIMYmMknwdphp0krknXd9Qyy1OSiXfnfMLKSLM3JX338PAQ1SW5FUunTp1EdUliUWpqqqguSfzQ4vgs4+vra7KMdJ9K+i7ZV1qupWv2YGz06NH6ORfllZ+H0atXL2zZsgXNmjXTrKNEROZivCIiW2D2sC4mJgY9evRATEwMsrKykJWVhZiYGNxzzz344YcfsHv3bly/fh2zZ8+2RH+JiMQYr4jIFph9Zuy5557DqlWrEBUVpd82YMAAODk5YdKkSUhISMC7776L8ePHa9pRIiJzMV4RkS0w+8zY+fPnK50L5O7ujgsXLgAA2rRpg2vXrtW8d0RENcB4RUS2wOzBWHh4OObMmYOrV6/qt129ehVz585Fjx49AABnz55FYGCgdr0kIqoGxisisgVmX6Zcs2YNRo0ahcDAQAQFBUGn0yE5ORktW7bEd999B+D2enYLFy7UvLOWlp6ebjRrT7oOlZubm8ky0iw1ydqAkswQKUnWIiDLgpRkFN28eVPUniRr0dXVVVSXZH9Js/okmWHS91qSASldy1PSL8k6dZJMO0D2GiX7XboWplR9jlcZGRlG97tkDVZAlhEmzQiW1CWNo5K6tMxuzMrKMllGGq8k606W/4FQU9Kse0m/pLFPsh6m9BiUJM9IYkPLli1F7Um+5yTfcVquMW32YCwsLAynTp3Cjh07cObMGSil0K5dOwwaNEj/4Rk9erRmHSQiqi7GKyKyBWYPxoDbv5KGDBmCfv36wdHRUfyriYiotjFeEVFdZ/acsdLSUrz66qsICAiAq6urfrmRhQsXYs2aNZp3kIiouhiviMgWmD0Ye+211/Dpp5/izTffNLi23KlTJ3zyySeado6IqCYYr4jIFpg9GFu3bh1WrVqFxx57zGCScOfOnfHf//7XrLp2796NESNGwN/fHzqdDlu2bDF4PicnB1OnTkVgYCCcnZ3Rvn17rFixwqBMQUEBpk2bBm9vb7i4uGDkyJG4fPmyuS+LiOohxisisgVmzxlLSUlB69atK2wvLS0Vr91YJjc3F126dMFTTz2Fhx56qMLzM2bMwK5du/D555/jrrvuws6dOzFlyhT4+/tj1KhRAIDp06fj+++/x8aNG+Hl5YVZs2Zh+PDhOHTokHgNtDIODg5GM0mkr0+SgSZd00qytpc0O0kyV0aStShtU5LpJN0PkgxV6TpokqxLLecVSTMSJf2XHoOSfS/pl7Q9yfsoydLSci07oH7Hq8aNGxvdX1pmQGqZtSh9nVrGD63WRZVmrkvak65PLMkilO7T4OBgk2XS09NFdXXo0MFkGcmaoADQtGlTk2Ukn1fp8SD5/pJkiEuzayXMHozdfffd2LNnD0JCQgy2f/PNN+jWrZtZdQ0dOhRDhw6t8vnffvsN48aNQ79+/QAAkyZNwscff4yDBw9i1KhRyMrKwpo1a7B+/XoMHDgQAPD5558jKCgIP//8MwYPHlxpvQUFBQZvhvR2AURkWxiviMgWmH2ZctGiRZg6dSreeOMNlJaWYtOmTXj66aexePFi/P3vf9e0c7169cLWrVuRkpICpRR27dqFM2fO6IPWoUOHUFRUhOjoaP3f+Pv7o2PHjoiPj6+y3iVLlsDDw0P/CAoK0rTfRFQ3MF4RkS0wezA2YsQIfPXVV9i+fTt0Oh3+/ve/49SpU/j+++8xaNAgTTv3/vvvo0OHDggMDISDgwOGDBmC5cuXo1evXgCAtLQ0ODg4VLhhXIsWLZCWllZlvQsWLNAvGpyVlYVLly5p2m8iqhsYr4jIFlTrPmODBw+u8pS6lt5//33s378fW7duRUhICHbv3o0pU6bAz89Pf5q/Mkopo/MlHB0dxdfrici2MV4RUV1XrcFYbcjPz8fzzz+PzZs344EHHgBwOwPq6NGj+Mc//oGBAwfC19cXhYWFyMjIMPi1mZ6ejqioKLPb/OWXX4xO0JQsHwPIJs5Kg2tCQoLJMtLJm9L+S1y/ft1kGclkXumEy969e5ss07FjR1Fd/v7+JstIlucAZMt9SCeKS8pJJ9dKlnaRLMcibU+yLIhkgrSWy4vUJmvEq5SUFKPvs3TytCShRTqB39ykCGO8vb1NlpEuNSY5jv/44w+TZSTLCQGy2KdlckxmZqaoLsl3hSQ+ArLvAF9fX1FdkuWJJHMlpUviSery9PQ0WSY3N1fUnoRoMNasWTNxZo5WX/hFRUUoKiqqcCDa2dnpM/7Cw8Nhb2+PmJgYjBkzBgCQmpqKEydO4M0339SkH0RkWxiviMjWiAZj7777rv7/r1+/jtdeew2DBw9GZGQkgNtZRDt27DB7sd2cnByDxUYTExNx9OhReHp6Ijg4GH379sWcOXPg7OyMkJAQxMXFYd26dVi2bBkAwMPDAxMmTMCsWbPg5eUFT09PzJ49G506dTJ6WYCI6i/GKyKyNaLB2Lhx4/T//9BDD+GVV17B1KlT9dueffZZfPjhh/j5558xY8YMceMHDx5E//799f+eOXOmvr1PP/0UGzduxIIFC/DYY4/hxo0bCAkJweuvv47Jkyfr/+add95B48aNMWbMGOTn52PAgAH49NNPzb5nDxHVD4xXRGRrdEo6GeD/c3V1xdGjRyvcSPHs2bPo1q2b6AZ1dU12djY8PDzw8ccfc86YEOeM3cY5Y7dpOWfs73//O7KysuDu7i5q25j6HK/27t1rdI4M54z9D+eM3ablnDHJa5TOGZP0q67OGRsyZIgm8crsW1t4eXlh8+bNFbZv2bIFXl5eNeoMEZGWGK+IyBaYnU358ssvY8KECYiNjdXPwdi/fz9++uknm194t7Cw0OgIXXoGSpJhIc0ak4zOJb9QANmZCenSSpJfmpL90LJlS1F7kl9rfn5+orokS29If+VIjgnpEj+Ssw/S90erWyFIzvwBtzP7TNFy2R2p+hyvAgMDjR6n0qXNJGfypWfZJO+xNLlCqyXXANmSQpLPqYeHh6g9yedGEocA2ZI7kuXiAFm/pJ95LX/MSN5HyZlEyXccIPtsSI55qy6H9OSTT6J9+/Z4//33sWnTJiil0KFDB+zbtw89e/bUrGNERDXFeEVEtqBa9xnr2bMnvvjiC637QkSkOcYrIqrrRHPGzL2xmZY3QiMiMgfjFRHZGtFgrHXr1li8eDGuXLlSZRmlFGJiYjB06FC8//77mnWQiMgcjFdEZGtElyljY2Px4osv4uWXX0bXrl0REREBf39/ODk5ISMjAydPnsRvv/0Ge3t7LFiwAJMmTbJ0v4mIKsV4RUS2xqz7jF2+fBnffPMNdu/ejaSkJOTn58Pb2xvdunXD4MGDMWzYMPF9o+qSsvv2ODg4GM30kWQ2ltVnijSbsl27dibLJCYmiupycXExWUZy7ylAlkkjObSGDBkiai88PNxkGel9xiT3L5JmOhm7L10ZaSaaJENVmk0puafcL7/8YrJMcnKyqD3J/ZeSkpJMliktLcWlS5c0uW9PfY9XCQkJRrPotMy8lR7DkuNO2i9Jm9J7J0oy8SQZwXl5eaL2JHFBmpkpuVeXNGNbyyxPyWuUfr4k/ZK8Rml2o6RfkszMmzdvomPHjprEK7Mm8AcGBmLGjBlm3bWaiMgaGK+IyFbY3s9CIiIionqEgzEiIiIiK+JgjIiIiMiKOBgjIiIisqJq3YG/vlqzZg2aNGlS5fOS7BFAlkUoTWKVZB5JszgkGUXSNQ0bNzZ96EiyTyUZnoCsX9K1QyV1SdfPk2QBSbPHJNlVkvcQAP744w+TZZo1a2ayjDSLTpJ5lJ+fb7JMbm4uHnzwQVGbDV3Tpk2Nfvalx4r0WJdwdXU1WUZ6TEluxis5hgHZ2pSSmCZdt1ESY6RrO0rqkvZL8l0h2Q+ALLtRum6y5BgsLi42WUaSeQrI7mYg2adaZmObPRg7e/YsvvvuOyQlJUGn0yE0NBSjR48WL/hMRFRbGK+IyBaYNRhbsmQJ/v73v6O0tBQ+Pj5QSuHq1auYP38+Fi9ejNmzZ1uqn0REZmG8IiJbIT7HtmvXLrz44ot44YUXcO3aNaSmpiItLU0f3ObPn4/du3dbsq9ERCKMV0RkS8RnxlauXImJEyfipZdeMtju6emJV155BWlpaVixYgX69OmjdR+JiMzCeEVEtkR8Zuz333/H448/XuXzjz/+OPbv369Jp4iIaoLxiohsifjM2B9//IG77rqryudDQ0ORlpamRZ+spkePHpqs9SbJKJJknwGydS6NZYCWJ8nglGSsSOuS9EuaASnplzRDVbLvpRlmWu0HQJbhI31/WrRoYbJM8+bNTZaRZlaVlpaaLCPJ7pOuLWdKQ4hXubm5RrO5pMewJCNMehxIMiWln3lJm9JsNsnnRrImo7Q9yfq3kgxPQJbFL40xkuxGaRyV7Avpey2JDdK1IiUk39GS73tpZrCE+MzYrVu3jO4Me3t7TTtGRFRdjFdEZEvMyqb85JNPqhxRavWLlohIC4xXRGQrxIOx4OBgrF692mQZIiJrY7wiIlsiHowlJSVZsBtERNphvCIiWyKeM3bu3DlL9oOISDOMV0RkS8Rnxtq2bYuAgAD0799f/zCWrWSLSkpKjGabSDOKJJkh0kwnSTnp+l+SLBnppGbJmpKS/SDNKJJk5Uj3w/Xr10XlJCTHhDSjSLLvJev1SduUHFvSrF/JsSVZY0+aLWpKQ4hX169fN5rxJc2Kk6wxKMmWBWTHp7RfkmzXM2fOiOqSfOYl61xKvwMkmZmSbGZAtu9v3Lghqsvf399kGekazJL3WhqTJSRr20q/VzMyMkyWkex36feXhPjMWFxcHP72t7/hypUrmDp1Klq1aoXQ0FBMmDABn3/+OVJSUjTrVHk3b97E9OnTERISAmdnZ0RFReHAgQP655VSeOmll+Dv7w9nZ2f069cPCQkJFukLEdkGxisisiXiwVjv3r3x4osv4ueff0ZmZiZ27dqFp556ComJiZg0aRKCg4MRFhameQcnTpyImJgYrF+/HsePH0d0dDQGDhyoD6Zvvvkmli1bhg8//BAHDhyAr68vBg0axGwpogaM8YqIbIlOSc8ZVyI/Px979+7Fjh07sHr1auTk5Gh+WtLNzQ3fffcdHnjgAf32rl27Yvjw4Xj11Vfh7++P6dOnY968eQBu36itRYsWeOONN/C3v/1N1E52djY8PDyQkJBg9KavUpJLRNLLP5IbgWp5KUl6Y1vJZUot29PyMqXkdLcUL1PeptVlyps3b6Jjx47IysqCu7u7qG2p+havfvvtN6M3r+Rlyv+x5cuUkrqkl+e0vEwpOSbq6mVKyRQa6WXKqKgoTeKV+MwYcHtg8Ouvv2LhwoXo3bs3mjVrhmeffRY5OTlYsWIFkpOTa9SZOxUXF6OkpAROTk4G252dnbF3714kJiYiLS0N0dHR+uccHR3Rt29fxMfHV1lvQUEBsrOzDR5EVL8wXhGRrRBP4O/bty8OHDiAVq1aoU+fPpg2bRr69u0rWnalutzc3BAZGYlXX30V7du3R4sWLbBhwwb8+9//Rps2bfS/nO7sQ4sWLXDx4sUq612yZAlefvlli/WbiKyL8YqIbIl4MBYfHw8/Pz/0798f/fr1Q58+fUTrb9XU+vXrMX78eAQEBMDOzg7du3fHo48+isOHD+vL3HlqUill9HTlggULMHPmTP2/s7OzERQUBG9vb6OnGqWn6iWn/aWn6rW4bFpGcopd+holp4MlZaTtSU4rS/eph4eHyTJartGZl5cnqkuyvpz0eJBcHpBcjpBeptTqmJeusWdKQ4hXPj4+Ro8H6eVxyfEpeX8BwMvLy2QZ6fw4yfHi6+srqkuS9Sa5DCZ5fYAsrkniECD7zEs/N5JpIdLLbTWY4VSB5LtC0p702JLEdy0vsUqIL1NmZmZi1apVaNKkCd544w0EBASgU6dOmDp1Kv71r3/h6tWrFulgq1atEBcXh5ycHFy6dAm///47ioqKEBoaqv8g3jm3ID093egvYEdHR7i7uxs8iKj+YLwiIlsiHoy5uLhgyJAhWLp0Kf7973/j2rVrePPNN9GkSRO8+eabCAwMRMeOHS3WURcXF/j5+SEjIwM7duzAqFGj9AEuJiZGX66wsBBxcXGIioqyWF+IqG5jvCIiW2LWQuHlubi4wNPTE56enmjWrBkaN26MU6dOadk3AMCOHTuglEJYWBjOnTuHOXPmICwsDE899RR0Oh2mT5+OxYsXo02bNmjTpg0WL16MJk2a4NFHH9W8L0RkmxiviKguEw/GSktLcfDgQcTGxmLXrl3Yt28fcnNz9Xe5/uijj9C/f3/NO5iVlYUFCxbg8uXL8PT0xEMPPYTXX39dP4dh7ty5yM/Px5QpU5CRkYGePXti586dms61IiLbwnhFRLZEfJ8xd3d35Obmws/PD/369UO/fv3Qv39/tGrVytJ9tLiy+/b88ccfdW4Cv5aTJDmBX17OGhP4GzUyPWtAOinbVifwZ2dno3nz5jW+b09DiFfnz5+vcxP4JaSTrCX3BpPeV4oT+G+TTOCX3NcM4AT+svY6d+6syX3GxGfG3nrrLfTv3x9t27atUYN1WX5+vtGgo+WXvbQuLQdQkg+iNOhq9eGR9l3ywZCuE5aVlWWyjORmu4Dsi0z6hSH5MpB+wUpuYCkpI70BpGQgKaHVWm8NIV799ttvRr+EpfcjCwgIMFlGumbt6dOnTZa5du2aqC7JF+sHH3wgqksyWOndu7fJMh06dBC1J8ny9PHxEdUl+cwHBQWJ6pL8uJLER0D23ST9PpG815LBufR7NTMz02QZyU2qpT+0JcSDMendoYmIrI3xiohsiTY/Z4mIiIioWjgYIyIiIrIiDsaIiIiIrIiDMSIiIiIrqvZNX+ujJk2aGM26ka5VpWUauCQFV5phJ8kOkWajSNrUMptSkt0ozfyT7FMnJydRXdJyEpK1E6XHliRrStJ3aXta7VPJLUzotpYtW8LV1bXK56UZwZKsWilPT0+TZSRZw8DtZaJM8fPzE9UlyZ7r1KmTyTKS1wdoGxckt8CQrrsqicnGjqnyJN8BWsYrSUas9JiXvEbJ69Mq+xvgmTEiIiIiq+JgjIiIiMiKOBgjIiIisiIOxoiIiIisiIMxIiIiIitiNmU5ptbHkq6/p+Wi1pIMSMmak4BsfTlJe4Asa0WSjSJdt1GyBtjVq1dFdUn2l7SujIwMk2WkWZ4uLi6alAFkWV+SfSpdk1CSnSQ55qXZUHQ7K9HYZ0z6Wc7NzdWqS5qurxocHGyyTGJioqguyXqYycnJJstIF6KW7HvJ2o6ALF5Js9IlmZnS9RYlGdRSkv5LYq20T5LvaEmGp5b7gGfGiIiIiKyIgzEiIiIiK+JgjIiIiMiKOBgjIiIisiIOxoiIiIisiNmU5bi7u8Pd3b3K56XrbEkyNaQZRRKSrA9AlrkozW6U7Atpho+EJAtIsnYZIMv4kmRyAbL9JV3TVJKBJc3ekWZqmSJdY0+SaSzJppRmGRPQtm1bo/FKmgkrOYal70tgYKDJMlqu3duiRQtROcm+kGRK+vr6itqT1CXNjJZ8trRcC1PL7zkt15qV9Evad0mGu6Tv0uxaCZ4ZIyIiIrIiDsaIiIiIrIiDMSIiIiIr4mCMiIiIyIo4GCMiIiKyImZTlmNvb280G0OaUaRltpCENGtR0n/p+psSkrqkfdcqOxCQ7QdpFpCkX9LjQZLVJt0Pkja1yoAEZFmekoyv2v7s2LK8vDyjWdnSDDvJ2ofSdS4lx5Q0xkiOdWPZpOVlZmaaLOPv72+yjHQdT29vb5NlpJnRkmx5LbPgtYy10n5Jvge0vAOBJLtWst+13Fc8M0ZERERkRfVmMLZ8+XKEhobCyckJ4eHh2LNnj7W7RERUKcYrIiqvXgzGvvrqK0yfPh0vvPACjhw5gt69e2Po0KFITk62dteIiAwwXhHRnerFYGzZsmWYMGECJk6ciPbt2+Pdd99FUFAQVqxYYe2uEREZYLwiojvZ/AT+wsJCHDp0CPPnzzfYHh0djfj4+Er/pqCgwGDSalZWFgAgOzvbaFvSyczSSYta0XLZIS0nSUr2l5Z9l06IlUze1HICv5TkuJHuL60m8EtJ3mvJfi/7DNbXZZG0jFc5OTkm25L2yRTJJH9AdkxJY4zksyWdUC9ZtkYSP/Ly8kTtSV6jNSbwS/apNKZJjhvpUnCSuFbbMUESQ8uOKy36ZvODsWvXrqGkpKTCGmUtWrRAWlpapX+zZMkSvPzyyxW2S9cjJCLLun79umg9UlujZbwKDw+3SB+JyDxaxCubH4yVufOXgVKqyl8LCxYswMyZM/X/zszMREhICJKTk2vlCyA7OxtBQUG4dOmSODXb1tqs7+1Zo8363h5w+6xPcHAwPD09a6U9a2G8qltt1vf2rNFmfW8P0DZe2fxgzNvbG3Z2dhV+Vaanp1f49VnG0dGx0vvmeHh41NqbCNy+R05ttmeNNut7e9Zos763B2h7CbUuYbyq223W9/as0WZ9bw/QJl7ZfMRzcHBAeHg4YmJiDLbHxMQgKirKSr0iIqqI8YqIKmPzZ8YAYObMmXj88ccRERGByMhIrFq1CsnJyZg8ebK1u0ZEZIDxiojuVC8GY2PHjsX169fxyiuvIDU1FR07dsT27dsREhIi+ntHR0csWrRIvORHTdV2e9Zos763Z40263t71mqztjFe1b0263t71mizvrendZs6VV9zyImIiIhsgM3PGSMiIiKyZRyMEREREVkRB2NEREREVsTBGBEREZEVNfjB2PLlyxEaGgonJyeEh4djz549FmtrxYoV6Ny5s/6mdJGRkfjxxx8t1h4ApKSk4K9//Su8vLzQpEkTdO3aFYcOHbJYezdv3sT06dMREhICZ2dnREVF4cCBA5rVv3v3bowYMQL+/v7Q6XTYsmWL/rmioiLMmzcPnTp1gouLC/z9/fHEE0/gypUrFmkPAJ588knodDqDx7333lvt9iRt5uTkYOrUqQgMDISzszPat29fo0WmlyxZgh49esDNzQ0+Pj4YPXo0Tp8+bVBm06ZNGDx4MLy9vaHT6XD06FGLtlfe3/72N+h0Orz77rvVbrO+YLzSFuMV45UW7ZVX3XjVoAdjX331FaZPn44XXngBR44cQe/evTF06FAkJydbpL3AwEAsXboUBw8exMGDB3H//fdj1KhRSEhIsEh7GRkZuO+++2Bvb48ff/wRJ0+exNtvv42mTZtapD0AmDhxImJiYrB+/XocP34c0dHRGDhwIFJSUjSpPzc3F126dMGHH35Y4bm8vDwcPnwYCxcuxOHDh7Fp0yacOXMGI0eOtEh7ZYYMGYLU1FT9Y/v27dVuT9LmjBkz8NNPP+Hzzz/HqVOnMGPGDEybNg3fffddtdqLi4vDM888g/379yMmJgbFxcWIjo42WIQ5NzcX9913H5YuXVqtNsxtr8yWLVvw73//G/7+/jVu19YxXmmP8YrxSov2ytQoXqkG7J577lGTJ0822NauXTs1f/78WutDs2bN1CeffGKRuufNm6d69eplkbork5eXp+zs7NQPP/xgsL1Lly7qhRde0Lw9AGrz5s1Gy/z+++8KgLp48aJF2hs3bpwaNWpUjes2p827775bvfLKKwbbunfvrl588UVN2kxPT1cAVFxcXIXnEhMTFQB15MgRTdoy1t7ly5dVQECAOnHihAoJCVHvvPOOZm3aIsYrbTFeaY/xqvrxqsGeGSssLMShQ4cQHR1tsD06Ohrx8fEWb7+kpAQbN25Ebm4uIiMjLdLG1q1bERERgYcffhg+Pj7o1q0bVq9ebZG2AKC4uBglJSVwcnIy2O7s7Iy9e/darF1jsrKyoNPpLPrrOjY2Fj4+Pmjbti2efvpppKenW6wtAOjVqxe2bt2KlJQUKKWwa9cunDlzBoMHD9ak/qysLACotcW6K2uvtLQUjz/+OObMmYO77767VvpRlzFeaY/xivFKq/Y0iVeaDRdtTEpKigKg9u3bZ7D99ddfV23btrVYu8eOHVMuLi7Kzs5OeXh4qG3btlmsLUdHR+Xo6KgWLFigDh8+rFauXKmcnJzUZ599ZrE2IyMjVd++fVVKSooqLi5W69evVzqdziL7FCZ+aebn56vw8HD12GOPWay9jRs3qh9++EEdP35cbd26VXXp0kXdfffd6tatWxZrs6CgQD3xxBMKgGrcuLFycHBQ69at06S90tJSNWLEiCrPUGj9S7Oq9hYvXqwGDRqkSktLlVKqwZ8ZY7yyDMYrxist2tMiXtWL5ZBqQqfTGfxbKVVhm5bCwsJw9OhRZGZm4ttvv8W4ceMQFxeHDh06aN5WaWkpIiIisHjxYgBAt27dkJCQgBUrVuCJJ57QvD0AWL9+PcaPH4+AgADY2dmhe/fuePTRR3H48GGLtFeVoqIiPPLIIygtLcXy5cst1s7YsWP1/9+xY0dEREQgJCQE27Ztw4MPPmiRNt9//33s378fW7duRUhICHbv3o0pU6bAz88PAwcOrFHdU6dOxbFjx2rtzEBl7R06dAjvvfceDh8+bNHPoi1ivNIW4xXjVU3b0yxeaTJctEEFBQXKzs5Obdq0yWD7s88+q/r06VNr/RgwYICaNGmSReoODg5WEyZMMNi2fPly5e/vb5H2ysvJyVFXrlxRSik1ZswYNWzYMM3bQBW/NAsLC9Xo0aNV586d1bVr1yze3p1at26tli5dapE28/LylL29fYV5LhMmTFCDBw+uUVtTp05VgYGB6sKFC1WW0fKXZlXtvfPOO0qn0yk7Ozv9A4Bq1KiRCgkJqXG7tojxyrIYrxivqtueVvGqwc4Zc3BwQHh4OGJiYgy2x8TEICoqqtb6oZRCQUGBReq+7777KqTgnjlzRrwgcU24uLjAz88PGRkZ2LFjB0aNGmXxNoHbvzDHjBmDs2fP4ueff4aXl1ettFvm+vXruHTpEvz8/CxSf1FREYqKitCokeFH187ODqWlpdWqUymFqVOnYtOmTfj1118RGhqqRVer3d7jjz+OY8eO4ejRo/qHv78/5syZgx07dli0b3UV45VlMV4xXlW3Pa3iVYO+TDlz5kw8/vjjiIiIQGRkJFatWoXk5GRMnjzZIu09//zzGDp0KIKCgnDz5k1s3LgRsbGx+OmnnyzS3owZMxAVFYXFixdjzJgx+P3337Fq1SqsWrXKIu0BwI4dO6CUQlhYGM6dO4c5c+YgLCwMTz31lCb15+Tk4Ny5c/p/JyYm4ujRo/D09IS/vz/+/Oc/4/Dhw/jhhx9QUlKCtLQ0ALcnWzo4OGjanqenJ1566SU89NBD8PPzQ1JSEp5//nl4e3vjT3/6k0VeY3BwMPr27Ys5c+bA2dkZISEhiIuLw7p167Bs2bJqtffMM8/gyy+/xHfffQc3Nzf9PvPw8ICzszMA4MaNG0hOTtbfA6nsS9PX1xe+vr6atufl5VXhS8ne3h6+vr4ICwur1musDxivtMd4xXhV0/Y0i1c1Pndn4z766CMVEhKiHBwcVPfu3StNj9XK+PHj9W01b95cDRgwQO3cudNi7Sml1Pfff686duyoHB0dVbt27dSqVass2t5XX32lWrZsqRwcHJSvr6965plnVGZmpmb179q1SwGo8Bg3bpz+lHRlj127dmneXl5enoqOjlbNmzdX9vb2Kjg4WI0bN04lJydb7DUqpVRqaqp68sknlb+/v3JyclJhYWHq7bff1k8eNVdV+2zt2rX6MmvXrq20zKJFiyzS3p0a+gT+MoxX2mK8YrzSor07VSde6f5/Y0RERERkBQ12zhgRERFRXcDBGBEREZEVcTBGREREZEUcjBERERFZEQdjRERERFbEwRgRERGRFXEwRkRERGRFHIwRERERWREHY2RRL730Erp27VqjOnQ6HbZs2aJJf8x1/fp1+Pj4ICkpCQAQGxsLnU6HzMxMq/Tnww8/xMiRI63SNlF9x3ilLcYrOQ7G6rknn3wSo0ePtnY3KvXSSy9Bp9MZfSQlJSE1NRVDhw61Sh+XLFmCESNG4K677rJK+3d6+umnceDAAezdu9faXSHSHONVzTBe2S4OxshqZs+ejdTUVP0jMDAQr7zyisG2oKAg+Pr6wtHRsdb7l5+fjzVr1mDixIm13vadlFIoLi6Go6MjHn30UXzwwQfW7hJRg8J4Jcd4ZT4Oxhq4uLg43HPPPXB0dISfnx/mz5+P4uJi/fOlpaV444030Lp1azg6OiI4OBivv/66/vl58+ahbdu2aNKkCVq2bImFCxeiqKhI1Larqyt8fX31Dzs7O7i5uVXYVv60f1JSEnQ6Hb7++mv07t0bzs7O6NGjB86cOYMDBw4gIiICrq6uGDJkCK5evWrQ3tq1a9G+fXs4OTmhXbt2WL58udH+/fjjj2jcuDEiIyMrPHfo0CFERESgSZMmiIqKwunTpw2eX7FiBVq1agUHBweEhYVh/fr1+ufKXsPRo0f12zIzM6HT6RAbGwvgf5cXduzYgYiICDg6OmLPnj0AgJEjR2LLli3Iz88X7Wei+oLxqmqMVzbO7CXMyaaMGzdOjRo1qtLnLl++rJo0aaKmTJmiTp06pTZv3qy8vb0NVrafO3euatasmfr000/VuXPn1J49e9Tq1av1z7/66qtq3759KjExUW3dulW1aNFCvfHGG/rnFy1apLp06SLqa1Ur3QNQmzdvVkoplZiYqACodu3aqZ9++kmdPHlS3Xvvvap79+6qX79+au/everw4cOqdevWavLkyfo6Vq1apfz8/NS3336rLly4oL799lvl6empPv300yr789xzz6khQ4YYbNu1a5cCoHr27KliY2NVQkKC6t27t4qKitKX2bRpk7K3t1cfffSROn36tHr77beVnZ2d+vXXXw1ew5EjR/R/k5GRoQCoXbt2GbTTuXNntXPnTnXu3Dl17do1pZRSOTk5SqfTqdjYWNF+JbIVjFe3MV41PByM1XPGgtvzzz+vwsLCVGlpqX7bRx99pFxdXVVJSYnKzs5Wjo6OBsHMlDfffFOFh4fr/22p4PbJJ5/on9+wYYMCoH755Rf9tiVLlqiwsDD9v4OCgtSXX35pUO+rr76qIiMjq+zPqFGj1Pjx4w22lQWdn3/+Wb9t27ZtCoDKz89XSikVFRWlnn76aYO/e/jhh9WwYcMMXoMkuG3ZsqXSvpV94RDVJ4xXtzFeNTyNa+f8G9VFp06dQmRkJHQ6nX7bfffdh5ycHFy+fBlpaWkoKCjAgAEDqqzjX//6F959912cO3cOOTk5KC4uhru7u8X73rlzZ/3/t2jRAgDQqVMng23p6ekAgKtXr+LSpUuYMGECnn76aX2Z4uJieHh4VNlGfn4+nJycTLbv5+cHAEhPT0dwcDBOnTqFSZMmGZS/77778N5770lfnl5ERESl252dnZGXl2d2fUS2ivGK8ao+42CsAVNKGQS2sm3A7fRsZ2dno3+/f/9+PPLII3j55ZcxePBgeHh4YOPGjXj77bct1ucy9vb2+v8vew13bistLQUA/X9Xr16Nnj17GtRjZ2dXZRve3t7IyMgQt1/WTvltZcrv60aNGum3lalq3oqLi0ul22/cuIHmzZtX2Xei+obxivGqPuME/gasQ4cOiI+PN/iQxcfHw83NDQEBAWjTpg2cnZ3xyy+/VPr3+/btQ0hICF544QVERESgTZs2uHjxYm11X6xFixYICAjAhQsX0Lp1a4NHaGholX/XrVs3nDx50uz22rdvXyGVOz4+Hu3btwcAfVBKTU3VP19+cqwp58+fx61bt9CtWzez+0ZkqxivGK/qM54ZawCysrIqfHg8PT0xZcoUvPvuu5g2bRqmTp2K06dPY9GiRZg5cyYaNWoEJycnzJs3D3PnzoWDgwPuu+8+XL16FQkJCZgwYQJat26N5ORkbNy4ET169MC2bduwefNm67xIE1566SU8++yzcHd3x9ChQ1FQUICDBw8iIyMDM2fOrPRvBg8ejAULFiAjIwPNmjUTtzVnzhyMGTMG3bt3x4ABA/D9999j06ZN+PnnnwHcPmV/7733YunSpbjrrrtw7do1vPjii+L69+zZg5YtW6JVq1bivyGyFYxXjFcNkrUmq1HtGDdunAJQ4TFu3DillFKxsbGqR48eysHBQfn6+qp58+apoqIi/d+XlJSo1157TYWEhCh7e3sVHBysFi9erH9+zpw5ysvLS7m6uqqxY8eqd955R3l4eOift9SE2PKTScsmj2ZkZOi3rV271qAfSin1xRdfqK5duyoHBwfVrFkz1adPH7Vp0yajfbr33nvVypUrjbZ15MgRBUAlJibqty1fvly1bNlS2dvbq7Zt26p169YZ1FuWVeXs7Ky6du2qdu7cWemE2PLtlImOjlZLliwx2m8iW8R49T+MVw2LTqly53yJyMD27dsxe/ZsnDhxQj93wppOnDiBAQMG4MyZM0Yn8xJRw8N4Zbt4mZLIiGHDhuHs2bNISUlBUFCQtbuDK1euYN26dQxsRFQB45Xt4pkxIiIiIiuy/nlMIiIiogaMgzEiIiIiK+JgjIiIiMiKOBgjIiIisiIOxoiIiIisiIMxIiIiIiviYIyIiIjIijgYIyIiIrIiDsaIiIiIrIiDsTtcv34dPj4+SEpKMll29uzZePbZZy3fKSKiSpiKV7GxsdDpdMjMzAQA/PDDD+jWrRtKS0trr5NEZBIHY3dYsmQJRowYgbvuustk2blz52Lt2rVITEy0fMeIiO5gTrwCgOHDh0On0+HLL7+0bMeIyCwcjJWTn5+PNWvWYOLEiaLyPj4+iI6OxsqVKy3cMyIiQ+bGqzJPPfUUPvjgAwv1ioiqg4Oxcn788Uc0btwYkZGR+m0JCQl44IEH4O7uDjc3N/Tu3Rvnz5/XPz9y5Ehs2LDBGt0logassni1fft2tG3bFs7Ozujfv3+lly9HjhyJ33//HRcuXKjF3hKRMRyMlbN7925ERETo/52SkoI+ffrAyckJv/76Kw4dOoTx48ejuLhYX+aee+7BpUuXcPHiRWt0mYgaqDvj1aVLl/Dggw9i2LBhOHr0KCZOnIj58+dX+LuQkBD4+Phgz549tdldIjKisbU7UJckJSXB399f/++PPvoIHh4e2LhxI+zt7QEAbdu2NfibgIAA/d+GhITUXmeJqEG7M16tWLECLVu2xDvvvAOdToewsDAcP34cb7zxRoW/DQgIECUpEVHt4JmxcvLz8+Hk5KT/99GjR9G7d2/9QKwyzs7OAIC8vDyL94+IqMyd8erUqVO49957odPp9NvKX8Isz9nZmTGLqA7hYKwcb29vZGRk6P9dNtAy5saNGwCA5s2bW6xfRER3ujNeKaXEf3vjxg3GLKI6hIOxcrp164aTJ0/q/925c2fs2bMHRUVFVf7NiRMnYG9vj7vvvrs2ukhEBKBivOrQoQP2799vUObOfwPArVu3cP78eXTr1s3ifSQiGQ7Gyhk8eDASEhL0vzanTp2K7OxsPPLIIzh48CDOnj2L9evX4/Tp0/q/2bNnD3r37i06i0ZEpJU749XkyZNx/vx5zJw5E6dPn8aXX36JTz/9tMLf7d+/H46OjlVewiSi2sfBWDmdOnVCREQEvv76awCAl5cXfv31V+Tk5KBv374IDw/H6tWrDeaQbdiwAU8//bS1ukxEDdSd8So4OBjffvstvv/+e3Tp0gUrV67E4sWLK/zdhg0b8Nhjj6FJkya13WUiqoJOmTPRoAHYvn07Zs+ejRMnTqBRI+Nj1W3btmHOnDk4duwYGjdmYioR1S5z4hUAXL16Fe3atcPBgwcRGhpaCz0kIgmOIO4wbNgwnD17FikpKQgKCjJaNjc3F2vXruVAjIiswpx4BQCJiYlYvnw5B2JEdQzPjBERERFZkVXnjK1YsQKdO3eGu7s73N3dERkZiR9//NGgzKlTpzBy5Eh4eHjAzc0N9957L5KTk/XPFxQUYNq0afD29oaLiwtGjhyJy5cv1/ZLIaJ6jvGKiCzFqoOxwMBALF26FAcPHsTBgwdx//33Y9SoUUhISAAAnD9/Hr169UK7du0QGxuL//znP1i4cKHBjQ6nT5+OzZs3Y+PGjdi7dy9ycnIwfPhwlJSUWOtlEVE9xHhFRJZS5y5Tenp64q233sKECRPwyCOPwN7eHuvXr6+0bFZWFpo3b47169dj7NixAIArV64gKCgI27dvx+DBg2uz60TUwDBeEZEW6szM85KSEnzzzTfIzc1FZGQkSktLsW3bNsydOxeDBw/GkSNHEBoaigULFmD06NEAgEOHDqGoqAjR0dH6evz9/dGxY0fEx8dXGdwKCgpQUFCg/3dpaSlu3LgBLy8vg6VEiKh2KaVw8+ZN+Pv7i7IDrYXxiog0jVfKyo4dO6ZcXFyUnZ2d8vDwUNu2bVNKKZWamqoAqCZNmqhly5apI0eOqCVLliidTqdiY2OVUkp98cUXysHBoUKdgwYNUpMmTaqyzUWLFikAfPDBRx19XLp0yTIBp4YYr/jgg487H1rEK6ufGQsLC8PRo0eRmZmJb7/9FuPGjUNcXByaNm0KABg1ahRmzJgBAOjatSvi4+OxcuVK9O3bt8o6lVJGfzEuWLAAM2fO1P87KysLwcHBuHTpEtzd3bV5YURktuzsbAQFBcHNzc3aXakU4xURldEyXll9MObg4IDWrVsDACIiInDgwAG89957+OCDD9C4cWN06NDBoHz79u2xd+9eAICvry8KCwuRkZGBZs2a6cukp6cjKiqqyjYdHR3h6OhYYXtZlhQRWVddvfzGeEVEd9IiXtW5SRlKKRQUFMDBwQE9evQwWAcSAM6cOYOQkBAAQHh4OOzt7RETE6N/PjU1FSdOnDAa3IiItMB4RURasOqZseeffx5Dhw5FUFAQbt68iY0bNyI2NhY//fQTAGDOnDkYO3Ys+vTpg/79++Onn37C999/j9jYWACAh4cHJkyYgFmzZsHLywuenp6YPXs2OnXqhIEDB1rxlRFRfcN4RUQWU+NZZzUwfvx4FRISohwcHFTz5s3VgAED1M6dOw3KrFmzRrVu3Vo5OTmpLl26qC1bthg8n5+fr6ZOnao8PT2Vs7OzGj58uEpOTjarH1lZWQqAysrKqvFrIqLqq8ufRcYrIipPy89inbvPmDVkZ2fDw8MD69evR5MmTaos5+zsLKpPskulu718SntVpPNGioqKTJapbG5KZSTrcXp6epos4+LiImpP0i87OzvN6pLOASh/Q8+qSN5D4PYtC0yRvIcA8Mcff5gsk5GRYbJMYWGhqD0HBweTZfLz802Wyc3NxYMPPoisrCzOh6pCWbyKj4+Hq6trleVu3bolqq/8/LWaunbtmskykuMAuD2XzpQLFy6I6srMzDRZplOnTibLSGIaIIsLUh4eHibLeHt7i+qSfO8YO6bKk8Rbe3t7UV1ZWVkmy0iOG+kxn5eXZ7KM5PXl5OSgX79+msSrOjdnjIiIiKgh4WCMiIiIyIo4GCMiIiKyIg7GiIiIiKyIgzEiIiIiK7L6Hfjrkh49ehhd1kCaFSfJRpFknwG3M6dMMZYBWp4kk6a4uFizuiT9kmZASvolzVCV7HtpNqVW+wGQZQJJ358WLVqYLNO8eXOTZSRZs4B2maA3b94UtUe3s8uMfX6kmdG5ubladUmU8Sb9zAcHB5ssk5iYKKpLkuWZnJxssoz0+JTse+lnS/K9I/n8AbLMTEmmISCPRRKS/kuyv6V9ksQiyfeElvuAZ8aIiIiIrIiDMSIiIiIr4mCMiIiIyIo4GCMiIiKyIg7GiIiIiKyI2ZTllJSUoKSkpMrnpdkvkkwNabaepJyxPpcnyfyTrkUoWVNSsh9ycnJE7UkysKT74fr166JyEpJjQpo9Jtn30sw3SZuSY0ua9Ss5tiQZZlpmJ9V3bdu2NboenvSzLDkOpJnKgYGBJstI1yuUkGQNA7J9IcmU9PX1FbUnqUu6Lm+jRqbPmWi5Fqb0/ZHEd+lauhKSfkn7Lllnurazv3lmjIiIiMiKOBgjIiIisiIOxoiIiIisiIMxIiIiIiviYIyIiIjIiphNWY63t7fR7CTp+l+SjA5pdpKxtTLNJcn8k75GSQaWpIy0PUlmi3SfStZn03KNTulab5I1LKXHgySzVJJRJM2m1OqYl67jSbePK2NZs9IMO8nah9J1LiWZf5IygCxeGYvX5WVmZpos4+/vb7KMNJvZ29vbZBlpjNFyLV3J51R61wAJab8k3wPSrHQJSXatZL9rua94ZoyIiIjIijgYIyIiIrIiDsaIiIiIrMjsC55JSUnYs2cPkpKSkJeXh+bNm6Nbt26IjIzU9C7AREQ1xXhFRLZAPBj78ssv8f777+P333+Hj48PAgIC4OzsjBs3buD8+fNwcnLCY489hnnz5iEkJMSSfSYiMorxiohsiWgw1r17dzRq1AhPPvkkvv76awQHBxs8X1BQgN9++w0bN25EREQEli9fjocfftgiHbak/Px8o9km0mw9STlpXVpmQEqypqRre2m1np2075LsQOk6l1lZWSbL3Lp1S1SXJFNSmlGUn59vsow0o6hZs2aalJFkXALyDDlTpO+hMQ0lXrm6uhrNrpUeK1pmsEqOAy3X5XV1dRXVJemXJD5Ks4slx7GW2a7Sz5/kmMjIyBDVJd0XElpl4UozVCXfc5J9pWWGp2gw9uqrr+KBBx6o8nlHR0f069cP/fr1w2uvvYbExETNOkhEZA7GKyKyNaLBmLHAdidvb2/RPVaIiCyB8YqIbE21ri2cP38eL774Iv7yl78gPT0dAPDTTz8hISFB084REdUU4xUR1XVmD8bi4uLQqVMn/Pvf/8amTZv018aPHTuGRYsWmVXXkiVL0KNHD7i5ucHHxwejR4/G6dOnDcrodLpKH2+99Za+TEFBAaZNmwZvb2+4uLhg5MiRuHz5srkvjYjqGcYrIrIFZg/G5s+fj9deew0xMTEGE/j69++P3377zay64uLi8Mwzz2D//v2IiYlBcXExoqOjDZacSE1NNXj885//hE6nw0MPPaQvM336dGzevBkbN27E3r17kZOTg+HDh4smfRNR/cV4RUS2QKekaX3/n6urK44fP47Q0FC4ubnhP//5D1q2bImkpCS0a9dOnIVWmatXr8LHxwdxcXHo06dPpWVGjx6Nmzdv4pdffgFwOzOuefPmWL9+PcaOHQsAuHLlCoKCgrB9+3YMHjy4Qh0FBQUGGSrZ2dkICgpCenq60bXOpMFSmpEoIckO0XKtN+nhIMki0TKbUnJcSdavBICbN2+aLCN9r2tyvN9Jsk+lx5Yk00mSzSVtT3KcStrLzs6Gn58fsrKyxOsOGlOf41VCQoLRbEppXJBmzEpIjjtpv6RZlxKSDMHr169rUgbQdi1dSbarNLNR8hmU9ksSu6XxQ6vvCml2o3StVVNycnLQvXt3TeKV2WfGmjZtitTU1Arbjxw5goCAgBp1puyWA56enpU+/8cff2Dbtm2YMGGCftuhQ4dQVFSE6Oho/TZ/f3907NgR8fHxldazZMkSeHh46B9BQUE16jcR1U2MV0RkC8wejD366KOYN28e0tLSoNPpUFpain379mH27Nl44oknqt0RpRRmzpyJXr16oWPHjpWW+eyzz+Dm5oYHH3xQvy0tLQ0ODg4V7pnUokULpKWlVVrPggULkJWVpX9cunSp2v0morqL8YqIbIHZyyG9/vrrePLJJxEQEAClFDp06ICSkhI8+uijePHFF6vdkalTp+LYsWPYu3dvlWX++c9/4rHHHhOfaq3qNLejo6NmpymJqO5ivCIiW2D2YMze3h5ffPEFXnnlFRw5cgSlpaXo1q0b2rRpU+1OTJs2DVu3bsXu3bsRGBhYaZk9e/bg9OnT+Oqrrwy2+/r6orCwEBkZGQa/NtPT0xEVFVXtPhGR7WO8IiJbUO01TIKCgtCpUyc8+OCD1Q5sSilMnToVmzZtwq+//orQ0NAqy65Zswbh4eHo0qWLwfbw8HDY29sjJiZGvy01NRUnTpxgcCMiAIxXRFS3mX1mLC8vD9OmTcNnn30GADhz5gxatmyJZ599Fv7+/pg/f764rmeeeQZffvklvvvuO7i5uennTHh4eBhk+GRnZ+Obb77B22+/XaEODw8PTJgwAbNmzYKXlxc8PT0xe/ZsdOrUCQMHDjTrtZnKspNmAWmZSSO5PCFZuwwACgsLNWkPkGURSjJbpBlTkjUgr169KqpLsr+kdUmytKTZai4uLpqUAaqeVF6eZJ9KjhlAtkag5JjXMjsVqN/xysvLy2gGl/S90zIDUktarsvr4+NjsowkG65Fixai9rRcB1iS2V3+9irGSPapNF5JYrckxgCyfkn2lzQLXnLMa5ndLmH2J2zBggX4z3/+g9jYWIO5EAMHDqxwSt6UFStWICsrC/369YOfn5/+cWc9GzduhFIKf/nLXyqt55133sHo0aMxZswY3HfffWjSpAm+//57TRfxJCLbw3hFRLbA7PuMhYSE4KuvvsK9995rcN+ec+fOoXv37sjOzrZUXy0mOzsbHh4eOH/+vCb37dHyfl6SM1XS+2tJfslIz4xJfnVreWZMclzxzNj/SM6MSY5nLc+MScrcvHkToaGhmt1nrD7Hq7S0NJ4ZE54Zk8QZyVlZ6Zkenhm7rb6fGdMyXpn9CSu70eGdcnNzNb1JHxFRTTFeEZEtMHsw1qNHD2zbtk3/77KAtnr1akRGRmrXMyKiGmK8IiJbYPYE/iVLlmDIkCE4efIkiouL8d577yEhIQG//fYb4uLiLNFHIqJqYbwiIltg9mAsKioK8fHxeOutt9CqVSvs3LkT3bt3x2+//YZOnTpZoo+1xt3d3eh1X+k1fskcLi0n60rXJZNclpFeupHsC+l8DgkPDw+TZSRruAFAfn6+yTLBwcGiuiT7SzqPQTJfT7IGJCCbgyEhuWEpIJtTJJknaeYUVpPqc7zKysoy+hnTMi5I45VWxx0gix/Sz5bkcyOJ23/88YeoPcn6t5I4JG1T0h4g21+ZmZmiuiQx2dgcbHPrksyDlZQBZHGmqnsIlqflXEqzPjlFRUWYNGkSFi5cqE8VJyKqixiviMhWmDWss7e3x+bNmy3VFyIizTBeEZGtMPsc25/+9Cds2bLFAl0hItIW4xUR2QKzL/C3bt0ar776KuLj4xEeHl7hvkfPPvusZp0jIqoJxisisgVmD8Y++eQTNG3aFIcOHcKhQ4cMntPpdAxuRFRnMF4RkS0wezCWmJhoiX7UCfb29kazBKWZXlquVyUhzVqU9F/L7BBJXdK+a5mlJdkP0lUNtLpzNCDLapPuB0mbWmVAArJsNUlmptafnfocr5o1a2Y0+1t6rEiOA2mWteRzI40xWn4eJMeepO+tWrUStaflWrrNmjUzWebGjRuiuiRtNm/eXFSX5DUGBQWJ6mratKnJMl5eXibLSO+CLzkeJMef9HtCovbXuCAiIiIiPbNPN8ycObPS7TqdDk5OTmjdujVGjRolvt8HEZGlMF4RkS0wezB25MgRHD58GCUlJQgLC4NSCmfPnoWdnR3atWuH5cuXY9asWdi7dy86dOhgiT4TEYkwXhGRLTD7MuWoUaMwcOBAXLlyBYcOHcLhw4eRkpKCQYMG4S9/+QtSUlLQp08fzJgxwxL9JSISY7wiIlugU2auPxIQEICYmJgKvyITEhIQHR2NlJQUHD58GNHR0bh27ZqmnbWU7OxseHh4IDMz0+gEQOmukk521YqWyw5puUyTZH9p2XfpUkGFhYUmy2g5gV9KctxI95dWE/ilJO+1ZNJsdnY2mjZtiqysLPFkXGPqc7xKT09v8BP4pf3SajK2NC5oOYFfsmySlhP4pceN5DVKl5Wz1Qn82dnZ8PHx0SRemf1NkpWVhfT09ArB7erVq8jOzgZwe8dKvvDqmtLSUqNfdtIPomRNOC0HbNIBlGRdMi0HnFq+xqysLJNlcnNzRXXl5OSYLCNd109Ld94DqzLSfSoJ4M7OzibLSNemlAz+JDFB67hRn+NVYWGh0X5LP8uSAb50ACVpU8tBovQ1SspJPjPSAbtkPcmMjAxRXf/9739Nlrlw4YKornPnzpksk5aWJqpL8v60aNFCVNcDDzxgskzPnj1NltEyPkpimuS7RKpalynHjx+PzZs34/Lly0hJScHmzZsxYcIEjB49GgDw+++/o23btpp1koioOhiviMgWmH1m7OOPP8aMGTPwyCOP6C8LNW7cGOPGjcM777wDAGjXrh0++eQTbXtKRGQmxisisgVmD8ZcXV2xevVqvPPOO7hw4QKUUmjVqhVcXV31Zbp27aplH4mIqoXxiohsQbVn8KalpSE1NRVt27aFq6ur+No9EVFtY7wiorrM7MHY9evXMWDAALRt2xbDhg1DamoqAGDixImYNWuW5h0kIqouxisisgVmX6acMWMG7O3tkZycjPbt2+u3jx07FjNmzMDbb7+taQdrU35+vtEMCmmGnZbrVUmyQ6TZlFpmTUkyMyW3mpBktQCybMrTp0+L6srMzDRZRpLZqDVJ9k5AQICoLslx4+bmZrJM+ct5xkiyLrVcE1SqPscrNzc3o++h5DMKaHs7Gy1peWscyedB8pmXtiepSxLTACAkJMRkGekarJLbUdy8eVNUlyQbVHprC0nWZWBgoMkyTZo0EbUnyeit7bsimD0Y27lzJ3bs2FFhx7Rp0wYXL17UrGNERDXFeEVEtsDsy5S5ubmVjj6vXbsGR0dHTTpFRKQFxisisgVmD8b69OmDdevW6f+t0+lQWlqKt956C/3799e0c0RENcF4RUS2wOzLlG+99Rb69euHgwcPorCwEHPnzkVCQgJu3LiBffv2WaKPRETVwnhFRLbA7DNjHTp0wLFjx3DPPfdg0KBByM3NxYMPPogjR46gVatWZtW1e/dujBgxAv7+/tDpdNiyZYvB8zk5OZg6dSoCAwPh7OyM9u3bY8WKFQZlCgoKMG3aNHh7e8PFxQUjR47E5cuXzX1ZRFQPMV4RkS2o1irHvr6+ePnll2vceG5uLrp06YKnnnoKDz30UIXnZ8yYgV27duHzzz/HXXfdhZ07d2LKlCnw9/fHqFGjAADTp0/H999/j40bN8LLywuzZs3C8OHDcejQIbOzhFxcXIxmwVhjoXAtMyBrez3M2s4EDQsLE9UlyTIrW7fQFEkWoWQNSEDb91qSBSnpl5aZdlqtX2mu+hqvTK1NKf28S44pSWY0oF1ckKrt7G/peyR5jdLsfEn2tzTGSDIlpWuHSjK7pdnYEpL3R5pBLJkvKvm+1/J+haK9fuzYMXGFnTt3FpcdOnQohg4dWuXzv/32G8aNG4d+/foBACZNmoSPP/4YBw8exKhRo5CVlYU1a9Zg/fr1GDhwIADg888/R1BQEH7++WcMHjxY3Bciqh8Yr4jI1ogGY127doVOp4NSymDEXzYqLL9NOjKV6NWrF7Zu3Yrx48fD398fsbGxOHPmDN577z0AwKFDh1BUVITo6Gj93/j7+6Njx46Ij4+vMrgVFBSgoKBA/2/pWRAiqvsYr4jI1ojO8SYmJuLChQtITEzEt99+i9DQUCxfvhxHjx7F0aNHsXz5crRq1Qrffvutpp17//330aFDBwQGBsLBwQFDhgzB8uXL0atXLwC3lzhxcHBAs2bNDP6uRYsWSEtLq7LeJUuWwMPDQ/8ICgrStN9EZD2MV0Rka0RnxsrfAfjhhx/G+++/j2HDhum3de7cGUFBQVi4cCFGjx6tWefef/997N+/H1u3bkVISAh2796NKVOmwM/PT3+avzJ3/iK+04IFCzBz5kz9v7OzsxngiOoJxisisjVmT+A/fvw4QkNDK2wPDQ3FyZMnNekUcHuZnOeffx6bN2/GAw88AOB2ED169Cj+8Y9/YODAgfD19UVhYSEyMjIMfm2mp6cjKiqqyrodHR15w0eiBoDxiohsgdmDsfbt2+O1117DmjVr9JlkBQUFeO211wzWfqupoqIiFBUVVciWsbOz02edhYeHw97eHjExMRgzZgwAIDU1FSdOnMCbb75pdpu3bt0SZ7gYI8mkkWYnSeqSZr9ISLOTtMr8k66DlpSUpFld5effVOXGjRua1SX9IvXw8DBZRvpeS45jyVpv0mwoSaakZN1Jydp55qjP8aq4uNhoHJFmLUrmzUnjguT9k67vKMkQlGazSbIgJevkXr9+XdSepC7J2o6AbN3Jc+fOieqSxLVbt26J6vLy8hKVk5CsTSnJtK7sh1dlJO+Pu7u7yTKS+C9l9rf4ypUrMWLECAQFBaFLly4AgP/85z/Q6XT44YcfzKorJyfH4CBKTEzE0aNH4enpieDgYPTt2xdz5syBs7MzQkJCEBcXh3Xr1mHZsmUAbn95TZgwAbNmzYKXlxc8PT0xe/ZsdOrUyehlASJqGBiviMgWmD0Yu+eee5CYmIjPP/8c//3vf6GUwtixY/Hoo4+KVqov7+DBgwZLkpTNixg3bhw+/fRTbNy4EQsWLMBjjz2GGzduICQkBK+//jomT56s/5t33nkHjRs3xpgxY5Cfn48BAwbg008/1fT+SERkmxiviMgW6JSWdy2zUdnZ2fDw8EBaWpro1KQpvEwpr4uXKf+nvl+mlPS9bHJ6VlaWJp/F+qgsXqWmphrdR9LLlJJy0rgguZRUVy9TSm4ZwsuU/1Pblykl9wSUXqaUHFuS+JOdnY3AwEBN4pXoE/bbb7+JK8zNzUVCQkK1O0REVBOMV0Rka0SDsSeeeAKDBg3C119/jZycnErLnDx5Es8//zxat26Nw4cPa9pJIiIpxisisjWiy5RFRUX4+OOP8eGHH+L8+fNo27Yt/P394eTkhIyMDPz3v//VL8C7YMECdOzYsTb6rpmy0/43btwweqpRenpdQnraX3IaXno5QtJ/6dwVrdazk54Sl1wOzM3NFdUluZQi3aeSfknXW5TMYZKshQnILo1K2pNeBpKc9peuCdq8efManfZvKPHq2rVrmlzKlVw+ln5OJXVJY4zkeJF+TiWZvJJ+VTW4v5Ok75JLmYAsQ1WahSzZX6mpqaK6JJ95yfQFKV9fX5NlJFM9AFnfJfsqOzsb/v7+mlymNHvO2OHDh7Fnzx4kJSUhPz8f3t7e6NatG/r37w9PT88adcZaOBj7Hw7GbuNg7DZbHIyVV5/jFQdjHIyZUwbgYKxMXRyMmT3zu3v37ujevXuNGiUiqg2MV0RkC2SnZ4iIiIjIIjgYIyIiIrIiDsaIiIiIrEi7u4XWAzk5OUYnpmt5Q1TpxGjJzTslk0UB2YRE6YRLyWRXyaRZ6Y1a09LSTJbRcnKt9P2RTOCX3qj1jz/+MFnGzc1NVJfkhoxaJqRIkhS0nChOt28aaixJRprsIYlr0onykrqk6/lJ6pJ+TiWf+aysLJNlpPFKckPXq1eviuqSkH4HSPolXZ9ZcqNZ6THYrFkzk2Uk8b1ly5ai9iTJE02aNDFZRst4pdmZsevXr+Pdd9/VqjoiIothvCKiuqRGgzGlFHbs2IExY8bA398fr7/+ulb9IiLSFOMVEdVV1RqMJSUl4e9//ztCQkIwbNgwODk5Ydu2baJLSUREtYnxiojqOvFgrKCgABs2bMCAAQPQvn17nDhxAsuWLUOjRo0wf/58DBw4UHwzPyIiS2K8IiJbIp7AHxAQgA4dOuCvf/0r/vWvf+kn3P3lL3+xWOdqm6urK1xdXat8XjIhHZBNbpdOPK3tSfdSkv5L+iWdkC6ZTCmdXCshncSq1UoEUtL3UHIHfsmEeul+kE4gNkXSb4mGEK8aN25s9DOm5aR7LSfKS49hSZtarmQiOdYld26Xtic91iUT16X7NDg42GSZ9PR0UV0dOnQwWUa6+kjTpk1NlpF8/0qPB0kSiSRua/mdIz4zVlJSAp1OB51Ox1+URFSnMV4RkS0RD8ZSU1MxadIkbNiwAb6+vnjooYewefNm8a8vIqLawnhFRLZEPBhzcnLCY489hl9//RXHjx9H+/bt8eyzz6K4uBivv/46YmJiNLtUQURUE4xXRGRLqpVN2apVK7z22mu4ePEifvjhBxQUFGD48OFo0aKF1v0jIqoRxisiqutqdAf+Ro0aYdiwYRg2bBiuXr2K9evXa9UvIiJNMV4RUV1l9mDs7Nmz+O6775CUlASdToeWLVti9OjRCA0NxcyZMy3Rx1pTVFRkNGNDuqyN5PKHNOtDUpd0Hoxk+Rtpv7QizdKSLDsh3Q+SLE9p5qxkf2mZPSbNbpTUJTm2pJfyJOUk+13r468+x6uUlBSjS/hIM9mMZZCXkX5OpZ8bCW9vb5NlpJnKkgxByXJkkuWEANnnQcsYk5mZKapLEov8/f1FdV2/ft1kGV9fX1FdkuWJsrOzTZaRHMvSujw9PU2Wyc3NFbUnYdZgbMmSJfj73/+O0tJS+Pj4QCmFq1evYt68eVi8eDFmz56tWceIiGqC8YqIbIX4Z+iuXbvw4osv4oUXXsC1a9eQmpqKtLQ0XL16FfPnz8f8+fOxe/duS/aViEiE8YqIbIn4zNjKlSsxceJEvPTSSwbbPT098corryAtLQ0rVqxAnz59tO4jEZFZGK+IyJaIz4z9/vvvePzxx6t8/vHHH8f+/fs16RQRUU0wXhGRLREPxv744w/cddddVT4fGhrKhXeJqE5gvCIiWyK+THnr1i2jmVz29vbi7J26qqCgwOiaVZKMPkCW8SatS0KaRSjJPJJms0kyS/Py8kyWkWS1ALL12aRrdDo5OZksI1m7DJDvewnJvpes0QnIMmcla+NJ9jug3TqXkmNGoiHEq8DAQLi7u1f5vOQYAGTHsHRfSY5h6WdG8hmUZnlqFT88PDxE7UlijGQ9RkC2/qF0jV9JvyRlAMDLy0tUTkLyPkqyT6XZ5pLPhuSY13JtSrOyKT/55JMqU0e17NSd9S5cuBCbN29Geno6unXrhvfeew89evQAcPtNfPnll7Fq1SpkZGSgZ8+e+Oijj3D33XdbpD9EZBsYr4jIVogHY8HBwVi9erXJMlqbOHEiTpw4gfXr18Pf3x+ff/45Bg4ciJMnTyIgIABvvvkmli1bhk8//RRt27bFa6+9hkGDBuH06dPiXwtEVL8wXhGRLdEp6XleK8jPz4ebmxu+++47PPDAA/rtXbt2xfDhw/Hqq6/C398f06dPx7x58wDcPrXdokULvPHGG/jb3/4maic7OxseHh5ITk42etpfuqskp0q1vDkiL1PexsuU/6PVZUrpMa/VZcrs7GwEBwcjKyvL6GexLqrteJWUlNTgL1NKb/qamJhosozkNUrjAi9TmseWL1O2b99ek3glnsB/7ty5GjVUHcXFxSgpKalwcDg7O2Pv3r1ITExEWloaoqOj9c85Ojqib9++iI+Pr7LegoICZGdnGzyIqP5gvCIiWyIejLVt2xZBQUF44oknsHbtWiQlJVmwW7e5ubkhMjISr776Kq5cuYKSkhJ8/vnn+Pe//62/iSOACgv+tmjRwmim1JIlS+Dh4aF/BAUFWfR1EFHtYrwiIlsinjMWFxeHuLg4xMbGYurUqbh16xaCg4Nx//33o3///ujfvz8CAgI07+D69esxfvx4BAQEwM7ODt27d8ejjz6Kw4cP68vcedpbKWX0VPiCBQsM1qXLzs5GUFAQmjRpYvQykPTSouTyj6QMIF8bUELLy5SS07yS07bSy26Stcukp6gll9S0XO9TWpfkvZa+Rsmx6uzsbLKM9DKQhJZrb5rSEOJVbm6u0c+r9JKaJBZJL1PeuHHDZBlpvyRtStoDZJ+HlJQUk2Wk2b6Sz5Y0M1OyvqP0/dHy8qmWax1L+iWZhiJZLxOQ9UsSi7T8fhYPxnr37o3evXvjxRdfRFFREX777TfExsYiNjYWGzZsQEFBAVq3bo3Tp09r1jkAaNWqFeLi4pCbm4vs7Gz4+flh7NixCA0N1R+kaWlp8PPz0/9Nenp6hV+f5Tk6OooHQ0RkexiviMiWiC9Tlmdvb48+ffpgzpw5WLBgAaZMmQJXV1eLztNwcXGBn58fMjIysGPHDowaNUof4GJiYvTlCgsLERcXh6ioKIv1hYhsB+MVEdV1Zt1n7NatW4iPj8euXbsQGxuLAwcOIDQ0FH379sWKFSvQt29fzTu4Y8cOKKUQFhaGc+fOYc6cOQgLC8NTTz0FnU6H6dOnY/HixWjTpg3atGmDxYsXo0mTJnj00Uc17wsR2Q7GKyKyFeLBWN++fXHgwAG0atUKffr0wbRp09C3b1+jp9e1kJWVhQULFuDy5cvw9PTEQw89hNdff11//Xju3LnIz8/HlClT9DdR3LlzJ+/ZQ9SAMV4RkS0R32fM3t4efn5+GD16NPr164c+ffrA29vb0v2rFWX37bl27ZrRSedaTuCXsuUJ/JK6pBPEOYH/tvo+gT87Oxu+vr41vm9PQ4hXCQkJRgdxnMD/P/V9Ar+W91iUTuCXvMbansAvXVlDqwn8N2/eRMeOHTW5z5h4MJabm4s9e/YgNjYWu3btwtGjR9G2bVv07dsX/fr1Q9++fdG8efMadcZayoLbyZMnNQlukl+50i85yc1OJR8KKcmgB5BlQUoGF9IPj2TQU9XSN3eS7C/poMfOzs5kGel7LfnCkA7OJf2SfNlJBq6A7DVK9vvNmzfRunXrGge3hhCvUlJSNPnxqOWNiyUxUjqwy83NNVlGGq8kcVSycLz0hqiSAa70pqmSuqT9knympJ95yeBIGq8kx6AkxkgHpZK1oSX7tCyzWYvBmPgypYuLC4YMGYIhQ4YAuB009+7di127duHNN9/EY489hjZt2uDEiRM16hARUU0xXhGRLalWNiVwO9h5enrC09MTzZo1Q+PGjXHq1Ckt+0ZEpAnGKyKqy8RnxkpLS3Hw4EH9af99+/YhNzcXAQEB6N+/Pz766CP079/fkn0lIhJhvCIiWyIejDVt2hS5ubnw8/NDv379sGzZMvTv3x+tWrWyZP+IiMzGeEVEtkQ8GHvrrbfQv39/tG3b1pL9ISKqMcYrIrIl4sHY3/72N0v2o05wcHAwmkUnzU6SZKNIU34l2SHSLE9Jxooka1HapiRRV7ofJBmq0kwaSdallhlmWmYnSY9Byb6X9Evanlap4tL30JSGEK9MrU0pPYYl7530GJZkSkoyfaVtSuOHJBNPcksHaXuS26hIMjwBWRaydI1fSXaj8AYLon0hfa8lcUZ6qwkJyXeAlpnBEtWewE9ERERENcfBGBEREZEVcTBGREREZEUcjBERERFZEQdjRERERFYkzqZsCH755RejmSvSRWklWUzSxcQTEhJMlpFmrEj7L3H9+nWTZSSZO9LspN69e5ss07FjR1Fd/v7+Jss0a9ZMVJdk/TJpRqKkXHp6uqiurKwsk2WuXr2qWXuStd4k2VCSeui269evG834kmbFSTJYpVnWkvUkpf2SrBV55swZUV2SeCX5zEuzSiWZmdK1USX7XhrbJbFPutax5L2Wrk0pIVmHVJpBnJGRYbKMZL9LM2IleGaMiIiIyIo4GCMiIiKyIg7GiIiIiKyIgzEiIiIiK+JgjIiIiMiKmE1ZTmFhodHMRGnWoiTLRJo15unpabKMNGNFks0mXedSkmUn2Q8tW7YUtSfJAvLz8xPVJcl0cnd3F9UlOSak6y1Kssyk7480W9cUSbYoAKSkpJgsI8mclWbaEeDj42N0zVZpvMrLyzNZRnoMe3l5mSwjXT9Qciz4+vqK6pJkvUmy9SSvD5Bl4nl4eIjqkqzLK12bUhI/pLFPy8+qJAtS0p702JKsVaplJqgEz4wRERERWREHY0RERERWxMEYERERkRVxMEZERERkRRyMEREREVkRsynLmTVrltGsDklmIwBkZ2ebLCPNpmzXrp3JMomJiaK6XFxcTJaRrFcIyLKKJNkvQUFBovauXbtmsoxk/TlAvn6ZhGQdt8LCQlFdkgxVaTalZK26X375xWSZ5ORkUXuStd6SkpJMlpGugUjAb7/9ZjSLThKHACAgIMBkGekxfPr0aZNlJJ9lQJYZ98EHH4jqkmQbSta/7dChg6g9SZanj4+PqC5JVqw0jkoy6iXr2gKydTqlWbiS91oS36UZnpmZmSbLSDLSJZnIUjwzRkRERGRF9WYwtnz5coSGhsLJyQnh4eHYs2ePtbtERFQpxisiKq9eDMa++uorTJ8+HS+88AKOHDmC3r17Y+jQoeJLLEREtYXxiojuVC8GY8uWLcOECRMwceJEtG/fHu+++y6CgoKwYsUKa3eNiMgA4xUR3cnmJ/AXFhbi0KFDmD9/vsH26OhoxMfHV/o3BQUFBhOhyyYsmpr8J51cLJlEKJ1oKFmSQVqXlpOjtaqrqKhIVE6S8CBZfgmQLfEjXUpGsqyGdPKzZMkWaV2SiaWSZADp+yPZD5JjpqxMfV0WSct4Zeo9lizvA8g+N1p+TqVJKNJjXUJyPEmOYWnfJftB+v5IYpE09kneR8nEfGk5aV2S2Cd5jdK4IYmPku/esno0iVfKxqWkpCgAat++fQbbX3/9ddW2bdtK/2bRokUKAB988FFHH+fPn6+N8FHrGK/44KP+PbSIVzZ/ZqzMnbcrUEpVeQuDBQsWYObMmfp/Z2ZmIiQkBMnJyeLFW2siOzsbQUFBuHTpknhRVltrs763Z40263t7wO2zPsHBweLbyNgqxqu61WZ9b88abdb39gBt45XND8a8vb1hZ2eHtLQ0g+3p6elo0aJFpX/j6OhY6T1EPDw8au1NBAB3d/dabc8abdb39qzRZn1vDwAaNaoX01krYLyq223W9/as0WZ9bw/QJl7ZfMRzcHBAeHg4YmJiDLbHxMQgKirKSr0iIqqI8YqIKmPzZ8YAYObMmXj88ccRERGByMhIrFq1CsnJyZg8ebK1u0ZEZIDxiojuVC8GY2PHjsX169fxyiuvIDU1FR07dsT27dsREhIi+ntHR0csWrRItPyBFmq7PWu0Wd/bs0ab9b09a7VZ2xiv6l6b9b09a7RZ39vTuk2dUvU0h5yIiIjIBtj8nDEiIiIiW8bBGBEREZEVcTBGREREZEUcjBERERFZUYMfjC1fvhyhoaFwcnJCeHg49uzZY7G2VqxYgc6dO+tvShcZGYkff/zRYu0BQEpKCv7617/Cy8sLTZo0QdeuXXHo0CGLtXfz5k1Mnz4dISEhcHZ2RlRUFA4cOKBZ/bt378aIESPg7+8PnU6HLVu26J8rKirCvHnz0KlTJ7i4uMDf3x9PPPEErly5YpH2AODJJ5+ETqczeNx7773Vbk/SZk5ODqZOnYrAwEA4Ozujffv2NVpkesmSJejRowfc3Nzg4+OD0aNH4/Tp0wZlNm3ahMGDB8Pb2xs6nQ5Hjx61aHvl/e1vf4NOp8O7775b7TbrC8YrbTFeMV5p0V551Y1XDXow9tVXX2H69Ol44YUXcOTIEfTu3RtDhw5FcnKyRdoLDAzE0qVLcfDgQRw8eBD3338/Ro0ahYSEBIu0l5GRgfvuuw/29vb48ccfcfLkSbz99tto2rSpRdoDgIkTJyImJgbr16/H8ePHER0djYEDByIlJUWT+nNzc9GlSxd8+OGHFZ7Ly8vD4cOHsXDhQhw+fBibNm3CmTNnMHLkSIu0V2bIkCFITU3VP7Zv317t9iRtzpgxAz/99BM+//xznDp1CjNmzMC0adPw3XffVau9uLg4PPPMM9i/fz9iYmJQXFyM6Ohog4V5c3Nzcd9992Hp0qXVasPc9sps2bIF//73v+Hv71/jdm0d45X2GK8Yr7Ror0yN4lWNV7e0Yffcc4+aPHmywbZ27dqp+fPn11ofmjVrpj755BOL1D1v3jzVq1cvi9Rdmby8PGVnZ6d++OEHg+1dunRRL7zwgubtAVCbN282Wub3339XANTFixct0t64cePUqFGjaly3OW3efffd6pVXXjHY1r17d/Xiiy9q0mZ6eroCoOLi4io8l5iYqACoI0eOaNKWsfYuX76sAgIC1IkTJ1RISIh65513NGvTFjFeaYvxSnuMV9WPVw32zFhhYSEOHTqE6Ohog+3R0dGIj4+3ePslJSXYuHEjcnNzERkZaZE2tm7dioiICDz88MPw8fFBt27dsHr1aou0BQDFxcUoKSmBk5OTwXZnZ2fs3bvXYu0ak5WVBZ1OZ9Ff17GxsfDx8UHbtm3x9NNPIz093WJtAUCvXr2wdetWpKSkQCmFXbt24cyZMxg8eLAm9WdlZQFArS3WXVl7paWlePzxxzFnzhzcfffdtdKPuozxSnuMV4xXWrWnSbzSbLhoY1JSUhQAtW/fPoPtr7/+umrbtq3F2j127JhycXFRdnZ2ysPDQ23bts1ibTk6OipHR0e1YMECdfjwYbVy5Url5OSkPvvsM4u1GRkZqfr27atSUlJUcXGxWr9+vdLpdBbZpzDxSzM/P1+Fh4erxx57zGLtbdy4Uf3www/q+PHjauvWrapLly7q7rvvVrdu3bJYmwUFBeqJJ55QAFTjxo2Vg4ODWrdunSbtlZaWqhEjRlR5hkLrX5pVtbd48WI1aNAgVVpaqpRSDf7MGOOVZTBeMV5p0Z4W8apeLIdUEzqdzuDfSqkK27QUFhaGo0ePIjMzE99++y3GjRuHuLg4dOjQQfO2SktLERERgcWLFwMAunXrhoSEBKxYsQJPPPGE5u0BwPr16zF+/HgEBATAzs4O3bt3x6OPPorDhw9bpL2qFBUV4ZFHHkFpaSmWL19usXbGjh2r//+OHTsiIiICISEh2LZtGx588EGLtPn+++9j//792Lp1K0JCQrB7925MmTIFfn5+GDhwYI3qnjp1Ko4dO1ZrZwYqa+/QoUN47733cPjwYYt+Fm0R45W2GK8Yr2ranmbxSpPhog0qKChQdnZ2atOmTQbbn332WdWnT59a68eAAQPUpEmTLFJ3cHCwmjBhgsG25cuXK39/f4u0V15OTo66cuWKUkqpMWPGqGHDhmneBqr4pVlYWKhGjx6tOnfurK5du2bx9u7UunVrtXTpUou0mZeXp+zt7SvMc5kwYYIaPHhwjdqaOnWqCgwMVBcuXKiyjJa/NKtq75133lE6nU7Z2dnpHwBUo0aNVEhISI3btUWMV5bFeMV4Vd32tIpXDXbOmIODA8LDwxETE2OwPSYmBlFRUbXWD6UUCgoKLFL3fffdVyEF98yZM+IFiWvCxcUFfn5+yMjIwI4dOzBq1CiLtwnc/oU5ZswYnD17Fj///DO8vLxqpd0y169fx6VLl+Dn52eR+ouKilBUVIRGjQw/unZ2digtLa1WnUopTJ06FZs2bcKvv/6K0NBQLbpa7fYef/xxHDt2DEePHtU//P39MWfOHOzYscOifaurGK8si/GK8aq67WkVrxr0ZcqZM2fi8ccfR0REBCIjI7Fq1SokJydj8uTJFmnv+eefx9ChQxEUFISbN29i48aNiI2NxU8//WSR9mbMmIGoqCgsXrwYY8aMwe+//45Vq1Zh1apVFmkPAHbs2AGlFMLCwnDu3DnMmTMHYWFheOqppzSpPycnB+fOndP/OzExEUePHoWnpyf8/f3x5z//GYcPH8YPP/yAkpISpKWlAbg92dLBwUHT9jw9PfHSSy/hoYcegp+fH5KSkvD888/D29sbf/rTnyzyGoODg9G3b1/MmTMHzs7OCAkJQVxcHNatW4dly5ZVq71nnnkGX375Jb777ju4ubnp95mHhwecnZ0BADdu3EBycrL+HkhlX5q+vr7w9fXVtD0vL68KX0r29vbw9fVFWFhYtV5jfcB4pT3GK8armranWbyq8bk7G/fRRx+pkJAQ5eDgoLp3715peqxWxo8fr2+refPmasCAAWrnzp0Wa08ppb7//nvVsWNH5ejoqNq1a6dWrVpl0fa++uor1bJlS+Xg4KB8fX3VM888ozIzMzWrf9euXQpAhce4ceP0p6Qre+zatUvz9vLy8lR0dLRq3ry5sre3V8HBwWrcuHEqOTnZYq9RKaVSU1PVk08+qfz9/ZWTk5MKCwtTb7/9tn7yqLmq2mdr167Vl1m7dm2lZRYtWmSR9u7U0Cfwl2G80hbjFeOVFu3dqTrxSvf/GyMiIiIiK2iwc8aIiIiI6gIOxoiIiIisiIMxIiIiIiviYIyIiIjIijgYIyIiIrIiDsaIiIiIrIiDMSIiIiIr4mCMiIiIyIo4GCOLeumll9C1a9ca1aHT6bBlyxZN+mOu69evw8fHB0lJSQCA2NhY6HQ6ZGZmWqU/H374IUaOHGmVtonqO8YrbTFeyXEwVs89+eSTGD16tLW7UamXXnoJOp3O6CMpKQmpqakYOnSoVfq4ZMkSjBgxAnfddZdV2r/T008/jQMHDmDv3r3W7gqR5hivaobxynZxMEZWM3v2bKSmpuofgYGBeOWVVwy2BQUFwdfXF46OjrXev/z8fKxZswYTJ06s9bbvpJRCcXExHB0d8eijj+KDDz6wdpeIGhTGKznGK/NxMNbAxcXF4Z577oGjoyP8/Pwwf/58FBcX658vLS3FG2+8gdatW8PR0RHBwcF4/fXX9c/PmzcPbdu2RZMmTdCyZUssXLgQRUVForZdXV3h6+urf9jZ2cHNza3CtvKn/ZOSkqDT6fD111+jd+/ecHZ2Ro8ePXDmzBkcOHAAERERcHV1xZAhQ3D16lWD9tauXYv27dvDyckJ7dq1w/Lly43278cff0Tjxo0RGRlZ4blDhw4hIiICTZo0QVRUFE6fPm3w/IoVK9CqVSs4ODggLCwM69ev1z9X9hqOHj2q35aZmQmdTofY2FgA/7u8sGPHDkRERMDR0RF79uwBAIwcORJbtmxBfn6+aD8T1ReMV1VjvLJxZi9hTjZl3LhxatSoUZU+d/nyZdWkSRM1ZcoUderUKbV582bl7e1tsLL93LlzVbNmzdSnn36qzp07p/bs2aNWr16tf/7VV19V+/btU4mJiWrr1q2qRYsW6o033tA/v2jRItWlSxdRX6ta6R6A2rx5s1JKqcTERAVAtWvXTv3000/q5MmT6t5771Xdu3dX/fr1U3v37lWHDx9WrVu3VpMnT9bXsWrVKuXn56e+/fZbdeHCBfXtt98qT09P9emnn1bZn+eee04NGTLEYNuuXbsUANWzZ08VGxurEhISVO/evVVUVJS+zKZNm5S9vb366KOP1OnTp9Xbb7+t7Ozs1K+//mrwGo4cOaL/m4yMDAVA7dq1y6Cdzp07q507d6pz586pa9euKaWUysnJUTqdTsXGxor2K5GtYLy6jfGq4eFgrJ4zFtyef/55FRYWpkpLS/XbPvroI+Xq6qpKSkpUdna2cnR0NAhmprz55psqPDxc/29LBbdPPvlE//yGDRsUAPXLL7/oty1ZskSFhYXp/x0UFKS+/PJLg3pfffVVFRkZWWV/Ro0apcaPH2+wrSzo/Pzzz/pt27ZtUwBUfn6+UkqpqKgo9fTTTxv83cMPP6yGDRtm8BokwW3Lli2V9q3sC4eoPmG8uo3xquFpXDvn36guOnXqFCIjI6HT6fTb7rvvPuTk5ODy5ctIS0tDQUEBBgwYUGUd//rXv/Duu+/i3LlzyMnJQXFxMdzd3S3e986dO+v/v0WLFgCATp06GWxLT08HAFy9ehWXLl3ChAkT8PTTT+vLFBcXw8PDo8o28vPz4eTkZLJ9Pz8/AEB6ejqCg4Nx6tQpTJo0yaD8fffdh/fee0/68vQiIiIq3e7s7Iy8vDyz6yOyVYxXjFf1GQdjDZhSyiCwlW0DbqdnOzs7G/37/fv345FHHsHLL7+MwYMHw8PDAxs3bsTbb79tsT6Xsbe31/9/2Wu4c1tpaSkA6P+7evVq9OzZ06AeOzu7Ktvw9vZGRkaGuP2ydspvK1N+Xzdq1Ei/rUxV81ZcXFwq3X7jxg00b968yr4T1TeMV4xX9Rkn8DdgHTp0QHx8vMGHLD4+Hm5ubggICECbNm3g7OyMX375pdK/37dvH0JCQvDCCy8gIiICbdq0wcWLF2ur+2ItWrRAQEAALly4gNatWxs8QkNDq/y7bt264eTJk2a31759+wqp3PHx8Wjfvj0A6INSamqq/vnyk2NNOX/+PG7duoVu3bqZ3TciW8V4xXhVn/HMWAOQlZVV4cPj6emJKVOm4N1338W0adMwdepUnD59GosWLcLMmTPRqFEjODk5Yd68eZg7dy4cHBxw33334erVq0hISMCECRPQunVrJCcnY+PGjejRowe2bduGzZs3W+dFmvDSSy/h2Wefhbu7O4YOHYqCggIcPHgQGRkZmDlzZqV/M3jwYCxYsAAZGRlo1qyZuK05c+ZgzJgx6N69OwYMGIDvv/8emzZtws8//wzg9in7e++9F0uXLsVdd92Fa9eu4cUXXxTXv2fPHrRs2RKtWrUS/w2RrWC8YrxqkKw1WY1qx7hx4xSACo9x48YppZSKjY1VPXr0UA4ODsrX11fNmzdPFRUV6f++pKREvfbaayokJETZ29ur4OBgtXjxYv3zc+bMUV5eXsrV1VWNHTtWvfPOO8rDw0P/vKUmxJafTFo2eTQjI0O/be3atQb9UEqpL774QnXt2lU5ODioZs2aqT59+qhNmzYZ7dO9996rVq5cabStI0eOKAAqMTFRv2358uWqZcuWyt7eXrVt21atW7fOoN6yrCpnZ2fVtWtXtXPnzkonxJZvp0x0dLRasmSJ0X4T2SLGq/9hvGpYdEqVO+dLRAa2b9+O2bNn48SJE/q5E9Z04sQJDBgwAGfOnDE6mZeIGh7GK9vFy5RERgwbNgxnz55FSkoKgoKCrN0dXLlyBevWrWNgI6IKGK9sF8+MEREREVmR9c9jEhERETVgHIwRERERWREHY0RERERWxMEYERERkRVxMEZERERkRRyMEREREVkRB2NEREREVsSbvpLNuHXrFgoLC63dDSIii3NwcICTk5O1u0G1hIMxsgm3bt1CaGgo0tLSrN0VIiKL8/X1RWJiIgdkDQTvwE82ITs7Gx4eHrh06RLc3NyglDJ4AKjxtjuf06JOY/9v7raq6tKyzvLhQMs6q+qvVnVWZ39Iy5qzP8x9PyTHXHXeY1P7uCbHTVX7o6bHYk2PO2N9kuzn6rzHpvZzdY+bMllZWXB3dwfVfzwzRjbF3d3dYoOxul5nXewT62SddeXvrVWnJViqXqq7OIGfiIiIyIo4GCMiIiKyIg7GiIiIiKyIgzEiIiIiK+IEfrIp2dnZFpuke+dzWtRp7P+1mDhsiQnOZSw1abo6+8Oc/axlnebuD3PfD8kxV5332NQ+rslxU9X+qOmxWNPjzlifJPu5Ou+xqf1c3eOGGh4OxsgmKKXg6uqKoKAga3eFiMjiXF1dOTBrQDgYI5ug0+mQk5ODS5cu1Zv77mRnZyMoKIivqQ6rb68H4GuyBWWvR6fTWbsrVEs4GCOb4u7uXi+CbXl8TXVffXs9AF8TUV3CCfxEREREVsTBGBEREZEVcTBGNsHR0RGLFi2Co6OjtbuiGb6muq++vR6Ar8kW1LfXQ6ZxoXAiIiIiK+KZMSIiIiIr4mCMiIiIyIo4GCMiIiKyIg7GiIiIiKyIgzGyCcuXL0doaCicnJwQHh6OPXv2WLtL1fbSSy9Bp9MZPHx9fa3dLbHdu3djxIgR8Pf3h06nw5YtWwyeV0rhpZdegr+/P5ydndGvXz8kJCRYp7NCpl7Tk08+WeE9u/fee63TWYElS5agR48ecHNzg4+PD0aPHo3Tp08blLG190nymmztfVqxYgU6d+6sv1ltZGQkfvzxR/3ztvYeUfVxMEZ13ldffYXp06fjhRdewJEjR9C7d28MHToUycnJ1u5atd19991ITU3VP44fP27tLonl5uaiS5cu+PDDDyt9/s0338SyZcvw4Ycf4sCBA/D19cWgQYNw8+bNWu6pnKnXBABDhgwxeM+2b99eiz00T1xcHJ555hns378fMTExKC4uRnR0NHJzc/VlbO19krwmwLbep8DAQCxduhQHDx7EwYMHcf/992PUqFH6AZetvUdUA4qojrvnnnvU5MmTDba1a9dOzZ8/30o9qplFixapLl26WLsbmgCgNm/erP93aWmp8vX1VUuXLtVvu3XrlvLw8FArV660Qg/Nd+drUkqpcePGqVGjRlmlP1pIT09XAFRcXJxSqn68T3e+JqVs/31SSqlmzZqpTz75pF68RyTHM2NUpxUWFuLQoUOIjo422B4dHY34+Hgr9armzp49C39/f4SGhuKRRx7BhQsXrN0lTSQmJiItLc3g/XJ0dETfvn1t+v0CgNjYWPj4+KBt27Z4+umnkZ6ebu0uiWVlZQEAPD09AdSP9+nO11TGVt+nkpISbNy4Ebm5uYiMjKwX7xHJcTBGddq1a9dQUlKCFi1aGGxv0aIF0tLSrNSrmunZsyfWrVuHHTt2YPXq1UhLS0NUVBSuX79u7a7VWNl7Up/eLwAYOnQovvjiC/z66694++23ceDAAdx///0oKCiwdtdMUkph5syZ6NWrFzp27AjA9t+nyl4TYJvv0/Hjx+Hq6gpHR0dMnjwZmzdvRocOHWz+PSLzNLZ2B4gkdDqdwb+VUhW22YqhQ4fq/79Tp06IjIxEq1at8Nlnn2HmzJlW7Jl26tP7BQBjx47V/3/Hjh0RERGBkJAQbNu2DQ8++KAVe2ba1KlTcezYMezdu7fCc7b6PlX1mmzxfQoLC8PRo0eRmZmJb7/9FuPGjUNcXJz+eVt9j8g8PDNGdZq3tzfs7Owq/BJMT0+v8IvRVrm4uKBTp044e/astbtSY2VZofX5/QIAPz8/hISE1Pn3bNq0adi6dSt27dqFwMBA/XZbfp+qek2VsYX3ycHBAa1bt0ZERASWLFmCLl264L333rPp94jMx8EY1WkODg4IDw9HTEyMwfaYmBhERUVZqVfaKigowKlTp+Dn52ftrtRYaGgofH19Dd6vwsJCxMXF1Zv3CwCuX7+OS5cu1dn3TCmFqVOnYtOmTfj1118RGhpq8Lwtvk+mXlNl6vr7VBmlFAoKCmzyPaIasFrqAJHQxo0blb29vVqzZo06efKkmj59unJxcVFJSUnW7lq1zJo1S8XGxqoLFy6o/fv3q+HDhys3NzebeT03b95UR44cUUeOHFEA1LJly9SRI0fUxYsXlVJKLV26VHl4eKhNmzap48ePq7/85S/Kz89PZWdnW7nnVTP2mm7evKlmzZql4uPjVWJiotq1a5eKjIxUAQEBdfY1/d///Z/y8PBQsbGxKjU1Vf/Iy8vTl7G198nUa7LF92nBggVq9+7dKjExUR07dkw9//zzqlGjRmrnzp1KKdt7j6j6OBgjm/DRRx+pkJAQ5eDgoLp3726Qzm5rxo4dq/z8/JS9vb3y9/dXDz74oEpISLB2t8R27dqlAFR4jBs3Til1+7YJixYtUr6+vsrR0VH16dNHHT9+3LqdNsHYa8rLy1PR0dGqefPmyt7eXgUHB6tx48ap5ORka3e7SpW9FgBq7dq1+jK29j6Zek22+D6NHz9eH9eaN2+uBgwYoB+IKWV77xFVn04ppWrvPBwRERERlcc5Y0RERERWxMEYERERkRVxMEZERERkRRyMEREREVkRB2NEREREVsTBGBEREZEVcTBGREREZEUcjBERERFZEQdjREQALl26hGHDhsHd3R3t27fHd999V6HMN998Ay8vL1y9etUKPSSi+oqDMSIiAOPGjUNWVhb+9a9/4U9/+hPGjBmD8+fP65/PysrCc889h3/84x9o3ry5FXtKRPUNl0MiogYvLy8Prq6u2LdvHyIjIwEAbdu2xcyZMzF58mQAwOTJk3H69Gns2rXLml0lonqosbU7QERkbYWFhVBKwcXFRb/N1dUVt27dAgDEx8dj3bp1OHr0qJV6SET1GS9TElGD17RpU7Rr1w5vv/02MjIysGXLFvznP/9BVFQUioqKMGnSJCxYsABt27a1dleJqB7imTEiIgBr1qzBQw89BE9PTzRq1Agvvvgi7rnnHrz22mtQSmHevHnW7iIR1VOcM0ZE9P8VFRXhwoUL8PX1hYeHB86ePYsuXbpg586diIqKwssvv4y1a9eisLAQf/7zn/GPf/wDTk5O1u42Edk4DsaIiKowYMAAtG7dGh9//DHWrFmDl19+Gb/88gtcXV0xZMgQjB49Gi+//LK1u0lENo5zxoiIKvHpp5/i5MmTeOONNwAAP/74Ix5++GG0adMGfn5+mDBhArZv327lXhJRfcA5Y0REd7h27Rpmz56N5cuXo2nTpgAApRRyc3P1ZXJycsALC0SkBQ7GiIjuMHPmTPTs2RNjxozRbxs8eDBmzZqFyMhIuLq64v3338fEiROt2Esiqi84Z4yIqJxffvkFo0aNQkJCAkJCQvTbS0pK8MILL+Czzz5DYWEhHnroIbz33ntwdna2Ym+JqD7gYIyIiIjIijiBn4iIiMiKOBgjIiIisiIOxoiIiIisiIMxIiIiIiviYIyIiIjIijgYIyIiIrIiDsaIiIiIrIiDMSIiIiIr4mCMiIiIyIo4GCMiIiKyIg7GiIiIiKzo/wH0qDtT0hFQ3gAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 600x500 with 5 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# plot histograms\n",
    "\n",
    "id[0] = 96011\n",
    "id[1] = 96581\n",
    "\n",
    "vmin = 0\n",
    "vmax = 30\n",
    "\n",
    "fig, axs = plt.subplots(nrows=2, ncols=2, figsize=(6, 5), layout=\"constrained\")\n",
    "axs_flat = axs.flatten()\n",
    "axs[0, 0].set_xlim(0, 24)\n",
    "axs[0, 0].set_xticks(np.arange(0, 25, 3))\n",
    "\n",
    "ylimit1 = 0\n",
    "ylimit2 = 360\n",
    "dy = 360/25\n",
    "\n",
    "ybins = np.arange(ylimit1, ylimit2, dy)\n",
    "xbins = np.arange(0, 25, 1)     # hourly bins \n",
    "\n",
    "\n",
    "# loop through variables\n",
    "nv = 1\n",
    "nl = 0\n",
    "for v in vsh:\n",
    "\n",
    "    print(nv, v, y_labs[nv-1])\n",
    "\n",
    "    row = nv - 1\n",
    "    im = axs[row, col].pcolormesh(xbins, ybins, HW.T, cmap=cmap, vmin=vmin, vmax=vmax)\n",
    "    \n",
    "# loop through stations\n",
    "    ncol = 0\n",
    "    for n in range(2): \n",
    "        sid = id[n]\n",
    "        print(sid)\n",
    "            \n",
    "        data_in_path = f\"/Users/cdz/Desktop/Khafid_etal/Manuscript/2Dhist_stations_Data/hist2d_data_{sid}_{v}.npz\"\n",
    "        \n",
    "        data = np.load(data_in_path)\n",
    "\n",
    "        HW = data['H']\n",
    "\n",
    "        col = ncol\n",
    "\n",
    "#        print(\"nv = \", nv, \"mv = \", mv, \"row = \", row, \"col = \", col)\n",
    "        \n",
    "        axs[row, col].pcolormesh(xbins, ybins, HW.T, cmap=cmap)\n",
    "        axs[row, col].set_ylim(0, 361)\n",
    "        axs[row, col].set_yticks(np.arange(0,361,90))\n",
    "        axs[row, col].set_xticks(np.arange(0,25,3))\n",
    "        axs[row, col].set_ylabel(y_labs[nv-1])\n",
    "        axs[row, col].set_xlabel(x_labs)\n",
    "        axs[row, col].set_title(labs[nl], fontsize=10, loc='left')\n",
    "\n",
    "        ncol = ncol + 1\n",
    "        nl = nl+1\n",
    "\n",
    "    nv = nv + 1\n",
    "\n",
    "    \n",
    "cbar = fig.colorbar(im, ax=axs, orientation='horizontal', location='bottom',\n",
    "                    pad=0.02, shrink=0.6, aspect=40)\n",
    "cbar.set_label('%', fontsize=12)\n",
    "\n",
    "fig_path = f\"/Users/cdz/Desktop/Khafid_etal/Manuscript/Figures/Fig3.png\"\n",
    "fig.savefig(fig_path)   \n",
    "plt.show()\n",
    " "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "bedbd958-32ec-4ff4-9386-b1fd963a5a41",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "623f876b-139f-4224-ada2-c38ccc7558f5",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "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.11.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}