Weather_Stations_Stat.ipynb

{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "provenance": []
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    },
    "language_info": {
      "name": "python"
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "source": [
        "# Modules"
      ],
      "metadata": {
        "id": "bYBIymv8X98d"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "import numpy as np\n",
        "import matplotlib as mpl\n",
        "import matplotlib.pyplot as plt\n",
        "#import random"
      ],
      "metadata": {
        "id": "2Vs0H_bIYCAY"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "source": [
        "# For load file \"weather.csv\" into Jupyter Notebook Colab  \n",
        "(not needed if not using Colab)"
      ],
      "metadata": {
        "id": "niHGNJDH8uKM"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "from google.colab import files\n",
        "uploaded = files.upload()\n",
        "! ls"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 91
        },
        "id": "ru9WHyRl8u4i",
        "outputId": "7233784f-289a-4e51-c966-1d74290f778c"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": [
              "<IPython.core.display.HTML object>"
            ],
            "text/html": [
              "\n",
              "     <input type=\"file\" id=\"files-38134645-f3e1-485f-89c8-fade22069b2a\" name=\"files[]\" multiple disabled\n",
              "        style=\"border:none\" />\n",
              "     <output id=\"result-38134645-f3e1-485f-89c8-fade22069b2a\">\n",
              "      Upload widget is only available when the cell has been executed in the\n",
              "      current browser session. Please rerun this cell to enable.\n",
              "      </output>\n",
              "      <script>// Copyright 2017 Google LLC\n",
              "//\n",
              "// Licensed under the Apache License, Version 2.0 (the \"License\");\n",
              "// you may not use this file except in compliance with the License.\n",
              "// You may obtain a copy of the License at\n",
              "//\n",
              "//      http://www.apache.org/licenses/LICENSE-2.0\n",
              "//\n",
              "// Unless required by applicable law or agreed to in writing, software\n",
              "// distributed under the License is distributed on an \"AS IS\" BASIS,\n",
              "// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n",
              "// See the License for the specific language governing permissions and\n",
              "// limitations under the License.\n",
              "\n",
              "/**\n",
              " * @fileoverview Helpers for google.colab Python module.\n",
              " */\n",
              "(function(scope) {\n",
              "function span(text, styleAttributes = {}) {\n",
              "  const element = document.createElement('span');\n",
              "  element.textContent = text;\n",
              "  for (const key of Object.keys(styleAttributes)) {\n",
              "    element.style[key] = styleAttributes[key];\n",
              "  }\n",
              "  return element;\n",
              "}\n",
              "\n",
              "// Max number of bytes which will be uploaded at a time.\n",
              "const MAX_PAYLOAD_SIZE = 100 * 1024;\n",
              "\n",
              "function _uploadFiles(inputId, outputId) {\n",
              "  const steps = uploadFilesStep(inputId, outputId);\n",
              "  const outputElement = document.getElementById(outputId);\n",
              "  // Cache steps on the outputElement to make it available for the next call\n",
              "  // to uploadFilesContinue from Python.\n",
              "  outputElement.steps = steps;\n",
              "\n",
              "  return _uploadFilesContinue(outputId);\n",
              "}\n",
              "\n",
              "// This is roughly an async generator (not supported in the browser yet),\n",
              "// where there are multiple asynchronous steps and the Python side is going\n",
              "// to poll for completion of each step.\n",
              "// This uses a Promise to block the python side on completion of each step,\n",
              "// then passes the result of the previous step as the input to the next step.\n",
              "function _uploadFilesContinue(outputId) {\n",
              "  const outputElement = document.getElementById(outputId);\n",
              "  const steps = outputElement.steps;\n",
              "\n",
              "  const next = steps.next(outputElement.lastPromiseValue);\n",
              "  return Promise.resolve(next.value.promise).then((value) => {\n",
              "    // Cache the last promise value to make it available to the next\n",
              "    // step of the generator.\n",
              "    outputElement.lastPromiseValue = value;\n",
              "    return next.value.response;\n",
              "  });\n",
              "}\n",
              "\n",
              "/**\n",
              " * Generator function which is called between each async step of the upload\n",
              " * process.\n",
              " * @param {string} inputId Element ID of the input file picker element.\n",
              " * @param {string} outputId Element ID of the output display.\n",
              " * @return {!Iterable<!Object>} Iterable of next steps.\n",
              " */\n",
              "function* uploadFilesStep(inputId, outputId) {\n",
              "  const inputElement = document.getElementById(inputId);\n",
              "  inputElement.disabled = false;\n",
              "\n",
              "  const outputElement = document.getElementById(outputId);\n",
              "  outputElement.innerHTML = '';\n",
              "\n",
              "  const pickedPromise = new Promise((resolve) => {\n",
              "    inputElement.addEventListener('change', (e) => {\n",
              "      resolve(e.target.files);\n",
              "    });\n",
              "  });\n",
              "\n",
              "  const cancel = document.createElement('button');\n",
              "  inputElement.parentElement.appendChild(cancel);\n",
              "  cancel.textContent = 'Cancel upload';\n",
              "  const cancelPromise = new Promise((resolve) => {\n",
              "    cancel.onclick = () => {\n",
              "      resolve(null);\n",
              "    };\n",
              "  });\n",
              "\n",
              "  // Wait for the user to pick the files.\n",
              "  const files = yield {\n",
              "    promise: Promise.race([pickedPromise, cancelPromise]),\n",
              "    response: {\n",
              "      action: 'starting',\n",
              "    }\n",
              "  };\n",
              "\n",
              "  cancel.remove();\n",
              "\n",
              "  // Disable the input element since further picks are not allowed.\n",
              "  inputElement.disabled = true;\n",
              "\n",
              "  if (!files) {\n",
              "    return {\n",
              "      response: {\n",
              "        action: 'complete',\n",
              "      }\n",
              "    };\n",
              "  }\n",
              "\n",
              "  for (const file of files) {\n",
              "    const li = document.createElement('li');\n",
              "    li.append(span(file.name, {fontWeight: 'bold'}));\n",
              "    li.append(span(\n",
              "        `(${file.type || 'n/a'}) - ${file.size} bytes, ` +\n",
              "        `last modified: ${\n",
              "            file.lastModifiedDate ? file.lastModifiedDate.toLocaleDateString() :\n",
              "                                    'n/a'} - `));\n",
              "    const percent = span('0% done');\n",
              "    li.appendChild(percent);\n",
              "\n",
              "    outputElement.appendChild(li);\n",
              "\n",
              "    const fileDataPromise = new Promise((resolve) => {\n",
              "      const reader = new FileReader();\n",
              "      reader.onload = (e) => {\n",
              "        resolve(e.target.result);\n",
              "      };\n",
              "      reader.readAsArrayBuffer(file);\n",
              "    });\n",
              "    // Wait for the data to be ready.\n",
              "    let fileData = yield {\n",
              "      promise: fileDataPromise,\n",
              "      response: {\n",
              "        action: 'continue',\n",
              "      }\n",
              "    };\n",
              "\n",
              "    // Use a chunked sending to avoid message size limits. See b/62115660.\n",
              "    let position = 0;\n",
              "    do {\n",
              "      const length = Math.min(fileData.byteLength - position, MAX_PAYLOAD_SIZE);\n",
              "      const chunk = new Uint8Array(fileData, position, length);\n",
              "      position += length;\n",
              "\n",
              "      const base64 = btoa(String.fromCharCode.apply(null, chunk));\n",
              "      yield {\n",
              "        response: {\n",
              "          action: 'append',\n",
              "          file: file.name,\n",
              "          data: base64,\n",
              "        },\n",
              "      };\n",
              "\n",
              "      let percentDone = fileData.byteLength === 0 ?\n",
              "          100 :\n",
              "          Math.round((position / fileData.byteLength) * 100);\n",
              "      percent.textContent = `${percentDone}% done`;\n",
              "\n",
              "    } while (position < fileData.byteLength);\n",
              "  }\n",
              "\n",
              "  // All done.\n",
              "  yield {\n",
              "    response: {\n",
              "      action: 'complete',\n",
              "    }\n",
              "  };\n",
              "}\n",
              "\n",
              "scope.google = scope.google || {};\n",
              "scope.google.colab = scope.google.colab || {};\n",
              "scope.google.colab._files = {\n",
              "  _uploadFiles,\n",
              "  _uploadFilesContinue,\n",
              "};\n",
              "})(self);\n",
              "</script> "
            ]
          },
          "metadata": {}
        },
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Saving weather.csv to weather.csv\n",
            "sample_data  weather.csv\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "# Load the data into appropriate NumPy arrays\n",
        "\n",
        "\n",
        "\n",
        "**BE CARE**  :  all following 3 notebook cells  (including that contain the vector \"head\") need to be run\n",
        "\n"
      ],
      "metadata": {
        "id": "zcymN_nK-p5b"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "a_n = np.genfromtxt('weather.csv', delimiter=',', names=True)  # with names=True the column headings have become the field names\n",
        "#print(f'{a_n = }')\n",
        "print(f'{a_n.shape = }')"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "2mD49T3H9tV2",
        "outputId": "d7af253b-22e4-47c5-cfd2-aeb135a33e72"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "a_n.shape = (151,)\n"
          ]
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "a = np.genfromtxt('weather.csv', delimiter=',')\n",
        "#print(f'{a = }')\n",
        "a = np.delete(a, 0, axis=0)  # This  clean the header row\n",
        "print(f'{a.shape = }')"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "NH0q9D2munzt",
        "outputId": "6cff7816-40ca-4adf-b044-b3a8704cbc04"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "a.shape = (151, 493)\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "The first three columns corresponds to the year, month and day in which the observations took place.\n",
        "\n",
        "The remaining columns are organized in (49) groups of 10 columns each:\n",
        "\n",
        "Each group represent a set of atmospheric observations collected at a given weather station in a given day. The following observations are reported (in order):\n",
        "- Maximum Temperature (°C)\n",
        "- Minimum Temperature (°C)\n",
        "- Mean Temperature (°C)\n",
        "- Maximum Relative Humidity (%)\n",
        "- Minimum Relative Humidity (%)\n",
        "- Mean Relative Humidity (%)\n",
        "- Maximum Pressure (hPa)\n",
        "- Minimum Pressure (hPa)\n",
        "- Mean Pressure (hPa)\n",
        "- Rainfall (mm)"
      ],
      "metadata": {
        "id": "EpeigKoAiF-d"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "head = [('YEAR', '<f8'), ('MONTH', '<f8'), ('DAY', '<f8'), ('WS_1_MAX_T', '<f8'), ('WS_1_MIN_T', '<f8'), ('WS_1_MEAN_T', '<f8'), ('WS_1_MAX_H', '<f8'), ('WS_1_MIN_H', '<f8'), ('WS_1_MEAN_H', '<f8'), ('WS_1_MAX_P', '<f8'), ('WS_1_MIN_P', '<f8'), ('WS_1_MEAN_P', '<f8'), ('WS_1_RAIN', '<f8'), ('WS_130_MAX_T', '<f8'), ('WS_130_MIN_T', '<f8'), ('WS_130_MEAN_T', '<f8'), ('WS_130_MAX_H', '<f8'), ('WS_130_MIN_H', '<f8'), ('WS_130_MEAN_H', '<f8'), ('WS_130_MAX_P', '<f8'), ('WS_130_MIN_P', '<f8'), ('WS_130_MEAN_P', '<f8'), ('WS_130_RAIN', '<f8'), ('WS_5_MAX_T', '<f8'), ('WS_5_MIN_T', '<f8'), ('WS_5_MEAN_T', '<f8'), ('WS_5_MAX_H', '<f8'), ('WS_5_MIN_H', '<f8'), ('WS_5_MEAN_H', '<f8'), ('WS_5_MAX_P', '<f8'), ('WS_5_MIN_P', '<f8'), ('WS_5_MEAN_P', '<f8'), ('WS_5_RAIN', '<f8'), ('WS_10_MAX_T', '<f8'), ('WS_10_MIN_T', '<f8'), ('WS_10_MEAN_T', '<f8'), ('WS_10_MAX_H', '<f8'), ('WS_10_MIN_H', '<f8'), ('WS_10_MEAN_H', '<f8'), ('WS_10_MAX_P', '<f8'), ('WS_10_MIN_P', '<f8'), ('WS_10_MEAN_P', '<f8'), ('WS_10_RAIN', '<f8'), ('WS_140_MAX_T', '<f8'), ('WS_140_MIN_T', '<f8'), ('WS_140_MEAN_T', '<f8'), ('WS_140_MAX_H', '<f8'), ('WS_140_MIN_H', '<f8'), ('WS_140_MEAN_H', '<f8'), ('WS_140_MAX_P', '<f8'), ('WS_140_MIN_P', '<f8'), ('WS_140_MEAN_P', '<f8'), ('WS_140_RAIN', '<f8'), ('WS_142_MAX_T', '<f8'), ('WS_142_MIN_T', '<f8'), ('WS_142_MEAN_T', '<f8'), ('WS_142_MAX_H', '<f8'), ('WS_142_MIN_H', '<f8'), ('WS_142_MEAN_H', '<f8'), ('WS_142_MAX_P', '<f8'), ('WS_142_MIN_P', '<f8'), ('WS_142_MEAN_P', '<f8'), ('WS_142_RAIN', '<f8'), ('WS_143_MAX_T', '<f8'), ('WS_143_MIN_T', '<f8'), ('WS_143_MEAN_T', '<f8'), ('WS_143_MAX_H', '<f8'), ('WS_143_MIN_H', '<f8'), ('WS_143_MEAN_H', '<f8'), ('WS_143_MAX_P', '<f8'), ('WS_143_MIN_P', '<f8'), ('WS_143_MEAN_P', '<f8'), ('WS_143_RAIN', '<f8'), ('WS_145_MAX_T', '<f8'), ('WS_145_MIN_T', '<f8'), ('WS_145_MEAN_T', '<f8'), ('WS_145_MAX_H', '<f8'), ('WS_145_MIN_H', '<f8'), ('WS_145_MEAN_H', '<f8'), ('WS_145_MAX_P', '<f8'), ('WS_145_MIN_P', '<f8'), ('WS_145_MEAN_P', '<f8'), ('WS_145_RAIN', '<f8'), ('WS_21_MAX_T', '<f8'), ('WS_21_MIN_T', '<f8'), ('WS_21_MEAN_T', '<f8'), ('WS_21_MAX_H', '<f8'), ('WS_21_MIN_H', '<f8'), ('WS_21_MEAN_H', '<f8'), ('WS_21_MAX_P', '<f8'), ('WS_21_MIN_P', '<f8'), ('WS_21_MEAN_P', '<f8'), ('WS_21_RAIN', '<f8'), ('WS_22_MAX_T', '<f8'), ('WS_22_MIN_T', '<f8'), ('WS_22_MEAN_T', '<f8'), ('WS_22_MAX_H', '<f8'), ('WS_22_MIN_H', '<f8'), ('WS_22_MEAN_H', '<f8'), ('WS_22_MAX_P', '<f8'), ('WS_22_MIN_P', '<f8'), ('WS_22_MEAN_P', '<f8'), ('WS_22_RAIN', '<f8'), ('WS_151_MAX_T', '<f8'), ('WS_151_MIN_T', '<f8'), ('WS_151_MEAN_T', '<f8'), ('WS_151_MAX_H', '<f8'), ('WS_151_MIN_H', '<f8'), ('WS_151_MEAN_H', '<f8'), ('WS_151_MAX_P', '<f8'), ('WS_151_MIN_P', '<f8'), ('WS_151_MEAN_P', '<f8'), ('WS_151_RAIN', '<f8'), ('WS_24_MAX_T', '<f8'), ('WS_24_MIN_T', '<f8'), ('WS_24_MEAN_T', '<f8'), ('WS_24_MAX_H', '<f8'), ('WS_24_MIN_H', '<f8'), ('WS_24_MEAN_H', '<f8'), ('WS_24_MAX_P', '<f8'), ('WS_24_MIN_P', '<f8'), ('WS_24_MEAN_P', '<f8'), ('WS_24_RAIN', '<f8'), ('WS_152_MAX_T', '<f8'), ('WS_152_MIN_T', '<f8'), ('WS_152_MEAN_T', '<f8'), ('WS_152_MAX_H', '<f8'), ('WS_152_MIN_H', '<f8'), ('WS_152_MEAN_H', '<f8'), ('WS_152_MAX_P', '<f8'), ('WS_152_MIN_P', '<f8'), ('WS_152_MEAN_P', '<f8'), ('WS_152_RAIN', '<f8'), ('WS_154_MAX_T', '<f8'), ('WS_154_MIN_T', '<f8'), ('WS_154_MEAN_T', '<f8'), ('WS_154_MAX_H', '<f8'), ('WS_154_MIN_H', '<f8'), ('WS_154_MEAN_H', '<f8'), ('WS_154_MAX_P', '<f8'), ('WS_154_MIN_P', '<f8'), ('WS_154_MEAN_P', '<f8'), ('WS_154_RAIN', '<f8'), ('WS_28_MAX_T', '<f8'), ('WS_28_MIN_T', '<f8'), ('WS_28_MEAN_T', '<f8'), ('WS_28_MAX_H', '<f8'), ('WS_28_MIN_H', '<f8'), ('WS_28_MEAN_H', '<f8'), ('WS_28_MAX_P', '<f8'), ('WS_28_MIN_P', '<f8'), ('WS_28_MEAN_P', '<f8'), ('WS_28_RAIN', '<f8'), ('WS_29_MAX_T', '<f8'), ('WS_29_MIN_T', '<f8'), ('WS_29_MEAN_T', '<f8'), ('WS_29_MAX_H', '<f8'), ('WS_29_MIN_H', '<f8'), ('WS_29_MEAN_H', '<f8'), ('WS_29_MAX_P', '<f8'), ('WS_29_MIN_P', '<f8'), ('WS_29_MEAN_P', '<f8'), ('WS_29_RAIN', '<f8'), ('WS_156_MAX_T', '<f8'), ('WS_156_MIN_T', '<f8'), ('WS_156_MEAN_T', '<f8'), ('WS_156_MAX_H', '<f8'), ('WS_156_MIN_H', '<f8'), ('WS_156_MEAN_H', '<f8'), ('WS_156_MAX_P', '<f8'), ('WS_156_MIN_P', '<f8'), ('WS_156_MEAN_P', '<f8'), ('WS_156_RAIN', '<f8'), ('WS_158_MAX_T', '<f8'), ('WS_158_MIN_T', '<f8'), ('WS_158_MEAN_T', '<f8'), ('WS_158_MAX_H', '<f8'), ('WS_158_MIN_H', '<f8'), ('WS_158_MEAN_H', '<f8'), ('WS_158_MAX_P', '<f8'), ('WS_158_MIN_P', '<f8'), ('WS_158_MEAN_P', '<f8'), ('WS_158_RAIN', '<f8'), ('WS_159_MAX_T', '<f8'), ('WS_159_MIN_T', '<f8'), ('WS_159_MEAN_T', '<f8'), ('WS_159_MAX_H', '<f8'), ('WS_159_MIN_H', '<f8'), ('WS_159_MEAN_H', '<f8'), ('WS_159_MAX_P', '<f8'), ('WS_159_MIN_P', '<f8'), ('WS_159_MEAN_P', '<f8'), ('WS_159_RAIN', '<f8'), ('WS_161_MAX_T', '<f8'), ('WS_161_MIN_T', '<f8'), ('WS_161_MEAN_T', '<f8'), ('WS_161_MAX_H', '<f8'), ('WS_161_MIN_H', '<f8'), ('WS_161_MEAN_H', '<f8'), ('WS_161_MAX_P', '<f8'), ('WS_161_MIN_P', '<f8'), ('WS_161_MEAN_P', '<f8'), ('WS_161_RAIN', '<f8'), ('WS_163_MAX_T', '<f8'), ('WS_163_MIN_T', '<f8'), ('WS_163_MEAN_T', '<f8'), ('WS_163_MAX_H', '<f8'), ('WS_163_MIN_H', '<f8'), ('WS_163_MEAN_H', '<f8'), ('WS_163_MAX_P', '<f8'), ('WS_163_MIN_P', '<f8'), ('WS_163_MEAN_P', '<f8'), ('WS_163_RAIN', '<f8'), ('WS_37_MAX_T', '<f8'), ('WS_37_MIN_T', '<f8'), ('WS_37_MEAN_T', '<f8'), ('WS_37_MAX_H', '<f8'), ('WS_37_MIN_H', '<f8'), ('WS_37_MEAN_H', '<f8'), ('WS_37_MAX_P', '<f8'), ('WS_37_MIN_P', '<f8'), ('WS_37_MEAN_P', '<f8'), ('WS_37_RAIN', '<f8'), ('WS_42_MAX_T', '<f8'), ('WS_42_MIN_T', '<f8'), ('WS_42_MEAN_T', '<f8'), ('WS_42_MAX_H', '<f8'), ('WS_42_MIN_H', '<f8'), ('WS_42_MEAN_H', '<f8'), ('WS_42_MAX_P', '<f8'), ('WS_42_MIN_P', '<f8'), ('WS_42_MEAN_P', '<f8'), ('WS_42_RAIN', '<f8'), ('WS_43_MAX_T', '<f8'), ('WS_43_MIN_T', '<f8'), ('WS_43_MEAN_T', '<f8'), ('WS_43_MAX_H', '<f8'), ('WS_43_MIN_H', '<f8'), ('WS_43_MEAN_H', '<f8'), ('WS_43_MAX_P', '<f8'), ('WS_43_MIN_P', '<f8'), ('WS_43_MEAN_P', '<f8'), ('WS_43_RAIN', '<f8'), ('WS_170_MAX_T', '<f8'), ('WS_170_MIN_T', '<f8'), ('WS_170_MEAN_T', '<f8'), ('WS_170_MAX_H', '<f8'), ('WS_170_MIN_H', '<f8'), ('WS_170_MEAN_H', '<f8'), ('WS_170_MAX_P', '<f8'), ('WS_170_MIN_P', '<f8'), ('WS_170_MEAN_P', '<f8'), ('WS_170_RAIN', '<f8'), ('WS_172_MAX_T', '<f8'), ('WS_172_MIN_T', '<f8'), ('WS_172_MEAN_T', '<f8'), ('WS_172_MAX_H', '<f8'), ('WS_172_MIN_H', '<f8'), ('WS_172_MEAN_H', '<f8'), ('WS_172_MAX_P', '<f8'), ('WS_172_MIN_P', '<f8'), ('WS_172_MEAN_P', '<f8'), ('WS_172_RAIN', '<f8'), ('WS_174_MAX_T', '<f8'), ('WS_174_MIN_T', '<f8'), ('WS_174_MEAN_T', '<f8'), ('WS_174_MAX_H', '<f8'), ('WS_174_MIN_H', '<f8'), ('WS_174_MEAN_H', '<f8'), ('WS_174_MAX_P', '<f8'), ('WS_174_MIN_P', '<f8'), ('WS_174_MEAN_P', '<f8'), ('WS_174_RAIN', '<f8'), ('WS_176_MAX_T', '<f8'), ('WS_176_MIN_T', '<f8'), ('WS_176_MEAN_T', '<f8'), ('WS_176_MAX_H', '<f8'), ('WS_176_MIN_H', '<f8'), ('WS_176_MEAN_H', '<f8'), ('WS_176_MAX_P', '<f8'), ('WS_176_MIN_P', '<f8'), ('WS_176_MEAN_P', '<f8'), ('WS_176_RAIN', '<f8'), ('WS_179_MAX_T', '<f8'), ('WS_179_MIN_T', '<f8'), ('WS_179_MEAN_T', '<f8'), ('WS_179_MAX_H', '<f8'), ('WS_179_MIN_H', '<f8'), ('WS_179_MEAN_H', '<f8'), ('WS_179_MAX_P', '<f8'), ('WS_179_MIN_P', '<f8'), ('WS_179_MEAN_P', '<f8'), ('WS_179_RAIN', '<f8'), ('WS_187_MAX_T', '<f8'), ('WS_187_MIN_T', '<f8'), ('WS_187_MEAN_T', '<f8'), ('WS_187_MAX_H', '<f8'), ('WS_187_MIN_H', '<f8'), ('WS_187_MEAN_H', '<f8'), ('WS_187_MAX_P', '<f8'), ('WS_187_MIN_P', '<f8'), ('WS_187_MEAN_P', '<f8'), ('WS_187_RAIN', '<f8'), ('WS_191_MAX_T', '<f8'), ('WS_191_MIN_T', '<f8'), ('WS_191_MEAN_T', '<f8'), ('WS_191_MAX_H', '<f8'), ('WS_191_MIN_H', '<f8'), ('WS_191_MEAN_H', '<f8'), ('WS_191_MAX_P', '<f8'), ('WS_191_MIN_P', '<f8'), ('WS_191_MEAN_P', '<f8'), ('WS_191_RAIN', '<f8'), ('WS_192_MAX_T', '<f8'), ('WS_192_MIN_T', '<f8'), ('WS_192_MEAN_T', '<f8'), ('WS_192_MAX_H', '<f8'), ('WS_192_MIN_H', '<f8'), ('WS_192_MEAN_H', '<f8'), ('WS_192_MAX_P', '<f8'), ('WS_192_MIN_P', '<f8'), ('WS_192_MEAN_P', '<f8'), ('WS_192_RAIN', '<f8'), ('WS_65_MAX_T', '<f8'), ('WS_65_MIN_T', '<f8'), ('WS_65_MEAN_T', '<f8'), ('WS_65_MAX_H', '<f8'), ('WS_65_MIN_H', '<f8'), ('WS_65_MEAN_H', '<f8'), ('WS_65_MAX_P', '<f8'), ('WS_65_MIN_P', '<f8'), ('WS_65_MEAN_P', '<f8'), ('WS_65_RAIN', '<f8'), ('WS_193_MAX_T', '<f8'), ('WS_193_MIN_T', '<f8'), ('WS_193_MEAN_T', '<f8'), ('WS_193_MAX_H', '<f8'), ('WS_193_MIN_H', '<f8'), ('WS_193_MEAN_H', '<f8'), ('WS_193_MAX_P', '<f8'), ('WS_193_MIN_P', '<f8'), ('WS_193_MEAN_P', '<f8'), ('WS_193_RAIN', '<f8'), ('WS_67_MAX_T', '<f8'), ('WS_67_MIN_T', '<f8'), ('WS_67_MEAN_T', '<f8'), ('WS_67_MAX_H', '<f8'), ('WS_67_MIN_H', '<f8'), ('WS_67_MEAN_H', '<f8'), ('WS_67_MAX_P', '<f8'), ('WS_67_MIN_P', '<f8'), ('WS_67_MEAN_P', '<f8'), ('WS_67_RAIN', '<f8'), ('WS_194_MAX_T', '<f8'), ('WS_194_MIN_T', '<f8'), ('WS_194_MEAN_T', '<f8'), ('WS_194_MAX_H', '<f8'), ('WS_194_MIN_H', '<f8'), ('WS_194_MEAN_H', '<f8'), ('WS_194_MAX_P', '<f8'), ('WS_194_MIN_P', '<f8'), ('WS_194_MEAN_P', '<f8'), ('WS_194_RAIN', '<f8'), ('WS_195_MAX_T', '<f8'), ('WS_195_MIN_T', '<f8'), ('WS_195_MEAN_T', '<f8'), ('WS_195_MAX_H', '<f8'), ('WS_195_MIN_H', '<f8'), ('WS_195_MEAN_H', '<f8'), ('WS_195_MAX_P', '<f8'), ('WS_195_MIN_P', '<f8'), ('WS_195_MEAN_P', '<f8'), ('WS_195_RAIN', '<f8'), ('WS_198_MAX_T', '<f8'), ('WS_198_MIN_T', '<f8'), ('WS_198_MEAN_T', '<f8'), ('WS_198_MAX_H', '<f8'), ('WS_198_MIN_H', '<f8'), ('WS_198_MEAN_H', '<f8'), ('WS_198_MAX_P', '<f8'), ('WS_198_MIN_P', '<f8'), ('WS_198_MEAN_P', '<f8'), ('WS_198_RAIN', '<f8'), ('WS_199_MAX_T', '<f8'), ('WS_199_MIN_T', '<f8'), ('WS_199_MEAN_T', '<f8'), ('WS_199_MAX_H', '<f8'), ('WS_199_MIN_H', '<f8'), ('WS_199_MEAN_H', '<f8'), ('WS_199_MAX_P', '<f8'), ('WS_199_MIN_P', '<f8'), ('WS_199_MEAN_P', '<f8'), ('WS_199_RAIN', '<f8'), ('WS_72_MAX_T', '<f8'), ('WS_72_MIN_T', '<f8'), ('WS_72_MEAN_T', '<f8'), ('WS_72_MAX_H', '<f8'), ('WS_72_MIN_H', '<f8'), ('WS_72_MEAN_H', '<f8'), ('WS_72_MAX_P', '<f8'), ('WS_72_MIN_P', '<f8'), ('WS_72_MEAN_P', '<f8'), ('WS_72_RAIN', '<f8'), ('WS_202_MAX_T', '<f8'), ('WS_202_MIN_T', '<f8'), ('WS_202_MEAN_T', '<f8'), ('WS_202_MAX_H', '<f8'), ('WS_202_MIN_H', '<f8'), ('WS_202_MEAN_H', '<f8'), ('WS_202_MAX_P', '<f8'), ('WS_202_MIN_P', '<f8'), ('WS_202_MEAN_P', '<f8'), ('WS_202_RAIN', '<f8'), ('WS_204_MAX_T', '<f8'), ('WS_204_MIN_T', '<f8'), ('WS_204_MEAN_T', '<f8'), ('WS_204_MAX_H', '<f8'), ('WS_204_MIN_H', '<f8'), ('WS_204_MEAN_H', '<f8'), ('WS_204_MAX_P', '<f8'), ('WS_204_MIN_P', '<f8'), ('WS_204_MEAN_P', '<f8'), ('WS_204_RAIN', '<f8'), ('WS_205_MAX_T', '<f8'), ('WS_205_MIN_T', '<f8'), ('WS_205_MEAN_T', '<f8'), ('WS_205_MAX_H', '<f8'), ('WS_205_MIN_H', '<f8'), ('WS_205_MEAN_H', '<f8'), ('WS_205_MAX_P', '<f8'), ('WS_205_MIN_P', '<f8'), ('WS_205_MEAN_P', '<f8'), ('WS_205_RAIN', '<f8'), ('WS_85_MAX_T', '<f8'), ('WS_85_MIN_T', '<f8'), ('WS_85_MEAN_T', '<f8'), ('WS_85_MAX_H', '<f8'), ('WS_85_MIN_H', '<f8'), ('WS_85_MEAN_H', '<f8'), ('WS_85_MAX_P', '<f8'), ('WS_85_MIN_P', '<f8'), ('WS_85_MEAN_P', '<f8'), ('WS_85_RAIN', '<f8'), ('WS_86_MAX_T', '<f8'), ('WS_86_MIN_T', '<f8'), ('WS_86_MEAN_T', '<f8'), ('WS_86_MAX_H', '<f8'), ('WS_86_MIN_H', '<f8'), ('WS_86_MEAN_H', '<f8'), ('WS_86_MAX_P', '<f8'), ('WS_86_MIN_P', '<f8'), ('WS_86_MEAN_P', '<f8'), ('WS_86_RAIN', '<f8'), ('WS_92_MAX_T', '<f8'), ('WS_92_MIN_T', '<f8'), ('WS_92_MEAN_T', '<f8'), ('WS_92_MAX_H', '<f8'), ('WS_92_MIN_H', '<f8'), ('WS_92_MEAN_H', '<f8'), ('WS_92_MAX_P', '<f8'), ('WS_92_MIN_P', '<f8'), ('WS_92_MEAN_P', '<f8'), ('WS_92_RAIN', '<f8'), ('WS_108_MAX_T', '<f8'), ('WS_108_MIN_T', '<f8'), ('WS_108_MEAN_T', '<f8'), ('WS_108_MAX_H', '<f8'), ('WS_108_MIN_H', '<f8'), ('WS_108_MEAN_H', '<f8'), ('WS_108_MAX_P', '<f8'), ('WS_108_MIN_P', '<f8'), ('WS_108_MEAN_P', '<f8'), ('WS_108_RAIN', '<f8'), ('WS_111_MAX_T', '<f8'), ('WS_111_MIN_T', '<f8'), ('WS_111_MEAN_T', '<f8'), ('WS_111_MAX_H', '<f8'), ('WS_111_MIN_H', '<f8'), ('WS_111_MEAN_H', '<f8'), ('WS_111_MAX_P', '<f8'), ('WS_111_MIN_P', '<f8'), ('WS_111_MEAN_P', '<f8'), ('WS_111_RAIN', '<f8'), ('WS_124_MAX_T', '<f8'), ('WS_124_MIN_T', '<f8'), ('WS_124_MEAN_T', '<f8'), ('WS_124_MAX_H', '<f8'), ('WS_124_MIN_H', '<f8'), ('WS_124_MEAN_H', '<f8'), ('WS_124_MAX_P', '<f8'), ('WS_124_MIN_P', '<f8'), ('WS_124_MEAN_P', '<f8'), ('WS_124_RAIN', '<f8')]\n",
        "print(len(head))"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "UXJ0Symv8p9D",
        "outputId": "cee2b424-f740-4336-f4fd-18c89c7b9757"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "493\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "# Perform 1\n",
        "\n",
        "Calculate the minimum, maximum and average rainfall and standard deviation of the rainfall over all the stations and all the days (a single set of numbers for the entire dataset)."
      ],
      "metadata": {
        "id": "k4ctMt9QAICn"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "def Perform_1(b=a):   \n",
        "   #print(b.shape)\n",
        "\n",
        "   # Tuple to build an appropriate mask for rainfall-data\n",
        "   ht = (0, 0, 0) \n",
        "   tt = (0, 0, 0,  0, 0, 0, 0, 0, 0, 1) * 49\n",
        "   rt = (ht + tt)\n",
        "    \n",
        "   ini = np.array(rt)\n",
        "   #print(ini.shape)\n",
        "   ini2 = np.zeros((151,1))\n",
        "   #print(ini2.shape)\n",
        "   ma = np.array(ini + ini2, dtype=bool) # broadcasting\n",
        "   #print(\"Mask shape : \", ma.shape)\n",
        "   #print(ma)\n",
        "\n",
        "   p1 = b[ma]  # Filter rainfall over all the stations and all the days \n",
        "   #print(f'{p1  = }')\n",
        "\n",
        "   #Descriptive statistics\n",
        "   print(f\"{np.min(p1) = } \")\n",
        "   print(f\"{np.max(p1) = } \")\n",
        "   print(f\"{np.mean(p1) = } \")\n",
        "   print(f\"{np.std(p1) = } \")\n",
        "   return np.min(p1), np.max(p1), np.mean(p1), np.std(p1)\n",
        "\n"
      ],
      "metadata": {
        "id": "ysDCwlKT7YDz"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "rainfall_min,rainfall_max,rainfall_mean,rainfall_std = Perform_1()"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "FK5h145h7nKC",
        "outputId": "ae9b1a79-1d58-490a-f01a-56143072f940"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "np.min(p1) = 0.0 \n",
            "np.max(p1) = 110.8 \n",
            "np.mean(p1) = 1.232098932288147 \n",
            "np.std(p1) = 5.0243051172827 \n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "# Perform 2\n",
        "\n",
        "Calculate the minimum, maximum and average rainfall and standard deviation of the rainfall at each station over all the days (a different set of numbers for each station)."
      ],
      "metadata": {
        "id": "fR6nEdcT-chd"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "def Perform_2(b=a_n, name_station = 'WS_1_RAIN'):    # Note : you have to give in input the name of station\n",
        "\n",
        "   p2 = b[name_station]\n",
        "   print(\"Data about rainfall of station : \",name_station[:4] )\n",
        "   #Descriptive statistics\n",
        "   print(f\"{np.min(p2) = } \")\n",
        "   print(f\"{np.max(p2) = } \")\n",
        "   print(f\"{np.mean(p2) = } \")\n",
        "   print(f\"{np.std(p2) = } \")\n",
        "\n"
      ],
      "metadata": {
        "id": "glZhEcq7_fXa"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "Perform_2()"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "1kTKYocpFBuU",
        "outputId": "988f27fd-aa3b-4797-aedd-e116b638c3a5"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Data about rainfall of station :  WS_1\n",
            "np.min(p2) = 0.0 \n",
            "np.max(p2) = 13.4 \n",
            "np.mean(p2) = 0.6344370860927152 \n",
            "np.std(p2) = 2.024716360075423 \n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "# Perform 3\n",
        "\n",
        "Calculate the minimum, maximum and average rainfall and standard deviation of the rainfall across all stations at each day (a different set of numbers for each day)."
      ],
      "metadata": {
        "id": "SihzQtmhFZtF"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "def Perform_3(b=a):   \n",
        "   #print(b.shape)\n",
        "\n",
        "   # Tuple to build an appropriate mask for rainfall-data\n",
        "   ht = (0, 0, 0) \n",
        "   tt = (0, 0, 0,  0, 0, 0, 0, 0, 0, 1) * 49\n",
        "   rt = (ht + tt)\n",
        "    \n",
        "   ini = np.array(rt)\n",
        "   #print(ini.shape)\n",
        "   ini2 = np.zeros((151,1))\n",
        "   #print(ini2.shape)\n",
        "   ma = np.array(ini + ini2, dtype=bool) # broadcasting\n",
        "   #print(\"Mask shape : \", ma.shape)\n",
        "   #print(ma)\n",
        "\n",
        "   p3 = b[ma]  # Filter rainfall over all the stations and all the days \n",
        "   p3n = (np.reshape(p3, (151, 49)))\n",
        "   #print(p3n)\n",
        "   \n",
        "   for i in range(151):\n",
        "     print(\"Day : \", round(a[i][0]),round(a[i][1]),round(a[i][2])) \n",
        "     print(\"-  Rainfall mini   : \",np.min(p3n[i]))\n",
        "     print(\"-  Rainfall max    : \",np.max(p3n[i]))\n",
        "     print(\"-  Rainfall mean   : \",np.mean(p3n[i]))\n",
        "     print(\"-  Rainfall st dev : \",np.std(p3n[i]))\n",
        "\n",
        "   "
      ],
      "metadata": {
        "id": "5SV0hMtkHd28"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "Perform_3()"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "0BxPu4FWLZ_S",
        "outputId": "3d9737e1-bcc9-4e5b-9016-c75f95188000"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Day :  2022 1 1\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  1.0\n",
            "-  Rainfall mean   :  0.07142857142857142\n",
            "-  Rainfall st dev :  0.16659862556700858\n",
            "Day :  2022 1 2\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.8\n",
            "-  Rainfall mean   :  0.06938775510204083\n",
            "-  Rainfall st dev :  0.15280250242722354\n",
            "Day :  2022 1 3\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.4\n",
            "-  Rainfall mean   :  0.05714285714285714\n",
            "-  Rainfall st dev :  0.10497813183356478\n",
            "Day :  2022 1 4\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  12.8\n",
            "-  Rainfall mean   :  0.9346938775510204\n",
            "-  Rainfall st dev :  2.1348948280178135\n",
            "Day :  2022 1 5\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  26.8\n",
            "-  Rainfall mean   :  1.883673469387755\n",
            "-  Rainfall st dev :  4.075629258477939\n",
            "Day :  2022 1 6\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  26.8\n",
            "-  Rainfall mean   :  0.6367346938775511\n",
            "-  Rainfall st dev :  3.7872877526601267\n",
            "Day :  2022 1 7\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  25.6\n",
            "-  Rainfall mean   :  0.536734693877551\n",
            "-  Rainfall st dev :  3.618021439087156\n",
            "Day :  2022 1 8\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  25.6\n",
            "-  Rainfall mean   :  0.7244897959183674\n",
            "-  Rainfall st dev :  3.684636088820541\n",
            "Day :  2022 1 9\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  6.6\n",
            "-  Rainfall mean   :  0.5306122448979592\n",
            "-  Rainfall st dev :  1.325758016771438\n",
            "Day :  2022 1 10\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.4\n",
            "-  Rainfall mean   :  0.0489795918367347\n",
            "-  Rainfall st dev :  0.10902620725871902\n",
            "Day :  2022 1 11\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  25.6\n",
            "-  Rainfall mean   :  0.5408163265306123\n",
            "-  Rainfall st dev :  3.6177520569900716\n",
            "Day :  2022 1 12\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  25.6\n",
            "-  Rainfall mean   :  0.5265306122448979\n",
            "-  Rainfall st dev :  3.619154005903956\n",
            "Day :  2022 1 13\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.004081632653061225\n",
            "-  Rainfall st dev :  0.028278380531736768\n",
            "Day :  2022 1 14\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.004081632653061225\n",
            "-  Rainfall st dev :  0.028278380531736768\n",
            "Day :  2022 1 15\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 1 16\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.3\n",
            "-  Rainfall mean   :  0.014285714285714285\n",
            "-  Rainfall st dev :  0.05714285714285714\n",
            "Day :  2022 1 17\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.3\n",
            "-  Rainfall mean   :  0.018367346938775512\n",
            "-  Rainfall st dev :  0.06283593599322716\n",
            "Day :  2022 1 18\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.4\n",
            "-  Rainfall mean   :  0.012244897959183675\n",
            "-  Rainfall st dev :  0.06270323059484578\n",
            "Day :  2022 1 19\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.8\n",
            "-  Rainfall mean   :  0.04897959183673469\n",
            "-  Rainfall st dev :  0.1540241817478629\n",
            "Day :  2022 1 20\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  10.0\n",
            "-  Rainfall mean   :  0.2897959183673469\n",
            "-  Rainfall st dev :  1.4602550654970765\n",
            "Day :  2022 1 21\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  10.0\n",
            "-  Rainfall mean   :  0.2183673469387755\n",
            "-  Rainfall st dev :  1.4130114706104244\n",
            "Day :  2022 1 22\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.3\n",
            "-  Rainfall mean   :  0.014285714285714285\n",
            "-  Rainfall st dev :  0.05714285714285714\n",
            "Day :  2022 1 23\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 1 24\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.028571428571428574\n",
            "-  Rainfall st dev :  0.06998542122237653\n",
            "Day :  2022 1 25\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.4\n",
            "-  Rainfall mean   :  0.036734693877551024\n",
            "-  Rainfall st dev :  0.08735075330217429\n",
            "Day :  2022 1 26\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.4\n",
            "-  Rainfall mean   :  0.04081632653061224\n",
            "-  Rainfall st dev :  0.09016621238520214\n",
            "Day :  2022 1 27\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.6\n",
            "-  Rainfall mean   :  0.036734693877551024\n",
            "-  Rainfall st dev :  0.10438132124526059\n",
            "Day :  2022 1 28\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.6\n",
            "-  Rainfall mean   :  0.036734693877551024\n",
            "-  Rainfall st dev :  0.10438132124526059\n",
            "Day :  2022 1 29\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.4\n",
            "-  Rainfall mean   :  0.0163265306122449\n",
            "-  Rainfall st dev :  0.07914579359047068\n",
            "Day :  2022 1 30\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.004081632653061225\n",
            "-  Rainfall st dev :  0.028278380531736768\n",
            "Day :  2022 1 31\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.02040816326530612\n",
            "-  Rainfall st dev :  0.060540395813025824\n",
            "Day :  2022 2 1\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.4\n",
            "-  Rainfall mean   :  0.012244897959183675\n",
            "-  Rainfall st dev :  0.06270323059484578\n",
            "Day :  2022 2 2\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  6.2\n",
            "-  Rainfall mean   :  0.24081632653061225\n",
            "-  Rainfall st dev :  1.0728406632454963\n",
            "Day :  2022 2 3\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  1.8\n",
            "-  Rainfall mean   :  0.04081632653061224\n",
            "-  Rainfall st dev :  0.255485426196389\n",
            "Day :  2022 2 4\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  4.8\n",
            "-  Rainfall mean   :  0.1836734693877551\n",
            "-  Rainfall st dev :  0.7404761123558247\n",
            "Day :  2022 2 5\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  1.8\n",
            "-  Rainfall mean   :  0.07346938775510205\n",
            "-  Rainfall st dev :  0.2790475858797912\n",
            "Day :  2022 2 6\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  1.4\n",
            "-  Rainfall mean   :  0.06938775510204082\n",
            "-  Rainfall st dev :  0.21112341069606666\n",
            "Day :  2022 2 7\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  25.8\n",
            "-  Rainfall mean   :  0.7755102040816326\n",
            "-  Rainfall st dev :  3.696304310526709\n",
            "Day :  2022 2 8\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  25.8\n",
            "-  Rainfall mean   :  0.5387755102040817\n",
            "-  Rainfall st dev :  3.6464585182221034\n",
            "Day :  2022 2 9\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 2 10\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  2.0\n",
            "-  Rainfall mean   :  0.04081632653061224\n",
            "-  Rainfall st dev :  0.2827838053173677\n",
            "Day :  2022 2 11\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  3.6\n",
            "-  Rainfall mean   :  0.17142857142857143\n",
            "-  Rainfall st dev :  0.7171371656006363\n",
            "Day :  2022 2 12\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  4.6\n",
            "-  Rainfall mean   :  0.6285714285714286\n",
            "-  Rainfall st dev :  1.2917319491951216\n",
            "Day :  2022 2 13\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  3.4\n",
            "-  Rainfall mean   :  0.07755102040816327\n",
            "-  Rainfall st dev :  0.4811827715554828\n",
            "Day :  2022 2 14\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  36.2\n",
            "-  Rainfall mean   :  3.510204081632653\n",
            "-  Rainfall st dev :  6.951324786318167\n",
            "Day :  2022 2 15\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  15.2\n",
            "-  Rainfall mean   :  3.048979591836735\n",
            "-  Rainfall st dev :  4.069153371753001\n",
            "Day :  2022 2 16\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  15.0\n",
            "-  Rainfall mean   :  0.8020408163265307\n",
            "-  Rainfall st dev :  2.4083604208467997\n",
            "Day :  2022 2 17\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  4.0\n",
            "-  Rainfall mean   :  0.16326530612244897\n",
            "-  Rainfall st dev :  0.6106850878738049\n",
            "Day :  2022 2 18\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  1.2\n",
            "-  Rainfall mean   :  0.1897959183673469\n",
            "-  Rainfall st dev :  0.3435973545069237\n",
            "Day :  2022 2 19\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  3.6\n",
            "-  Rainfall mean   :  0.22653061224489798\n",
            "-  Rainfall st dev :  0.6113803419576707\n",
            "Day :  2022 2 20\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  2.2\n",
            "-  Rainfall mean   :  0.08979591836734695\n",
            "-  Rainfall st dev :  0.31767323601666453\n",
            "Day :  2022 2 21\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  6.6\n",
            "-  Rainfall mean   :  0.25102040816326526\n",
            "-  Rainfall st dev :  1.123791405035866\n",
            "Day :  2022 2 22\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  4.6\n",
            "-  Rainfall mean   :  0.11020408163265306\n",
            "-  Rainfall st dev :  0.6528571109594189\n",
            "Day :  2022 2 23\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  1.0\n",
            "-  Rainfall mean   :  0.036734693877551024\n",
            "-  Rainfall st dev :  0.17922039510285837\n",
            "Day :  2022 2 24\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 2 25\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.004081632653061225\n",
            "-  Rainfall st dev :  0.02827838053173677\n",
            "Day :  2022 2 26\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.004081632653061225\n",
            "-  Rainfall st dev :  0.02827838053173677\n",
            "Day :  2022 2 27\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 2 28\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 1\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.004081632653061225\n",
            "-  Rainfall st dev :  0.028278380531736768\n",
            "Day :  2022 3 2\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.004081632653061225\n",
            "-  Rainfall st dev :  0.028278380531736768\n",
            "Day :  2022 3 3\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 4\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  1.0\n",
            "-  Rainfall mean   :  0.09183673469387756\n",
            "-  Rainfall st dev :  0.2038774488774233\n",
            "Day :  2022 3 5\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  2.1\n",
            "-  Rainfall mean   :  0.14081632653061224\n",
            "-  Rainfall st dev :  0.4275790259202547\n",
            "Day :  2022 3 6\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 7\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 8\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.8\n",
            "-  Rainfall mean   :  0.0163265306122449\n",
            "-  Rainfall st dev :  0.11311352212694707\n",
            "Day :  2022 3 9\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 10\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.8\n",
            "-  Rainfall mean   :  0.0163265306122449\n",
            "-  Rainfall st dev :  0.11311352212694709\n",
            "Day :  2022 3 11\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.8\n",
            "-  Rainfall mean   :  0.04489795918367347\n",
            "-  Rainfall st dev :  0.14716535818220366\n",
            "Day :  2022 3 12\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  5.6\n",
            "-  Rainfall mean   :  0.4877551020408163\n",
            "-  Rainfall st dev :  1.0807466529118153\n",
            "Day :  2022 3 13\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  8.4\n",
            "-  Rainfall mean   :  0.589795918367347\n",
            "-  Rainfall st dev :  1.5544859243373457\n",
            "Day :  2022 3 14\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  13.8\n",
            "-  Rainfall mean   :  2.5816326530612246\n",
            "-  Rainfall st dev :  2.8169801075408247\n",
            "Day :  2022 3 15\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  15.2\n",
            "-  Rainfall mean   :  2.66734693877551\n",
            "-  Rainfall st dev :  3.0486121134721285\n",
            "Day :  2022 3 16\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  8.6\n",
            "-  Rainfall mean   :  0.436734693877551\n",
            "-  Rainfall st dev :  1.5111414630832536\n",
            "Day :  2022 3 17\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  1.6\n",
            "-  Rainfall mean   :  0.17755102040816326\n",
            "-  Rainfall st dev :  0.36324237487815497\n",
            "Day :  2022 3 18\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  1.6\n",
            "-  Rainfall mean   :  0.0653061224489796\n",
            "-  Rainfall st dev :  0.2437387530911611\n",
            "Day :  2022 3 19\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 20\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 21\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 22\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 23\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 24\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 25\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 26\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 27\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 28\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 29\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 3 30\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  21.4\n",
            "-  Rainfall mean   :  5.653061224489796\n",
            "-  Rainfall st dev :  3.917678510881018\n",
            "Day :  2022 3 31\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  14.0\n",
            "-  Rainfall mean   :  3.010204081632653\n",
            "-  Rainfall st dev :  4.09081823550029\n",
            "Day :  2022 4 1\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  21.8\n",
            "-  Rainfall mean   :  6.795918367346939\n",
            "-  Rainfall st dev :  5.9865140680711395\n",
            "Day :  2022 4 2\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  21.8\n",
            "-  Rainfall mean   :  6.795918367346939\n",
            "-  Rainfall st dev :  5.9865140680711395\n",
            "Day :  2022 4 3\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  21.8\n",
            "-  Rainfall mean   :  6.795918367346939\n",
            "-  Rainfall st dev :  5.9865140680711395\n",
            "Day :  2022 4 4\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  21.8\n",
            "-  Rainfall mean   :  6.795918367346939\n",
            "-  Rainfall st dev :  5.9865140680711395\n",
            "Day :  2022 4 5\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.012244897959183675\n",
            "-  Rainfall st dev :  0.04794832703865604\n",
            "Day :  2022 4 6\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  2.0\n",
            "-  Rainfall mean   :  0.10204081632653061\n",
            "-  Rainfall st dev :  0.3639296836760614\n",
            "Day :  2022 4 7\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  14.2\n",
            "-  Rainfall mean   :  0.8367346938775511\n",
            "-  Rainfall st dev :  2.655352229031943\n",
            "Day :  2022 4 8\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  14.2\n",
            "-  Rainfall mean   :  0.8367346938775511\n",
            "-  Rainfall st dev :  2.655352229031943\n",
            "Day :  2022 4 9\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  14.2\n",
            "-  Rainfall mean   :  0.8367346938775511\n",
            "-  Rainfall st dev :  2.655352229031943\n",
            "Day :  2022 4 10\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  14.2\n",
            "-  Rainfall mean   :  0.8367346938775511\n",
            "-  Rainfall st dev :  2.655352229031943\n",
            "Day :  2022 4 11\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  1.2\n",
            "-  Rainfall mean   :  0.02857142857142857\n",
            "-  Rainfall st dev :  0.17142857142857146\n",
            "Day :  2022 4 12\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  1.2\n",
            "-  Rainfall mean   :  0.024489795918367346\n",
            "-  Rainfall st dev :  0.1696702831904206\n",
            "Day :  2022 4 13\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 4 14\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.004081632653061225\n",
            "-  Rainfall st dev :  0.02827838053173677\n",
            "Day :  2022 4 15\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.004081632653061225\n",
            "-  Rainfall st dev :  0.02827838053173677\n",
            "Day :  2022 4 16\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.3\n",
            "-  Rainfall mean   :  0.006122448979591836\n",
            "-  Rainfall st dev :  0.04241757079760515\n",
            "Day :  2022 4 17\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.3\n",
            "-  Rainfall mean   :  0.006122448979591836\n",
            "-  Rainfall st dev :  0.04241757079760515\n",
            "Day :  2022 4 18\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 4 19\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.004081632653061225\n",
            "-  Rainfall st dev :  0.02827838053173677\n",
            "Day :  2022 4 20\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  4.2\n",
            "-  Rainfall mean   :  0.8020408163265306\n",
            "-  Rainfall st dev :  1.0854303756459203\n",
            "Day :  2022 4 21\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  11.2\n",
            "-  Rainfall mean   :  1.522448979591837\n",
            "-  Rainfall st dev :  2.4088497810791525\n",
            "Day :  2022 4 22\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  12.4\n",
            "-  Rainfall mean   :  0.8387755102040818\n",
            "-  Rainfall st dev :  2.073232350262691\n",
            "Day :  2022 4 23\n",
            "-  Rainfall mini   :  0.8\n",
            "-  Rainfall max    :  110.8\n",
            "-  Rainfall mean   :  35.881632653061224\n",
            "-  Rainfall st dev :  26.46378097432211\n",
            "Day :  2022 4 24\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  102.0\n",
            "-  Rainfall mean   :  8.910204081632653\n",
            "-  Rainfall st dev :  17.337175403702844\n",
            "Day :  2022 4 25\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  15.0\n",
            "-  Rainfall mean   :  0.8755102040816326\n",
            "-  Rainfall st dev :  2.518221806398795\n",
            "Day :  2022 4 26\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  15.7\n",
            "-  Rainfall mean   :  0.6489795918367347\n",
            "-  Rainfall st dev :  2.3896991911270553\n",
            "Day :  2022 4 27\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  1.8\n",
            "-  Rainfall mean   :  0.06326530612244897\n",
            "-  Rainfall st dev :  0.27006717044364514\n",
            "Day :  2022 4 28\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.004081632653061225\n",
            "-  Rainfall st dev :  0.02827838053173677\n",
            "Day :  2022 4 29\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.0\n",
            "-  Rainfall mean   :  0.0\n",
            "-  Rainfall st dev :  0.0\n",
            "Day :  2022 4 30\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  5.2\n",
            "-  Rainfall mean   :  0.363265306122449\n",
            "-  Rainfall st dev :  0.8858741509070657\n",
            "Day :  2022 5 1\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  3.0\n",
            "-  Rainfall mean   :  0.2530612244897959\n",
            "-  Rainfall st dev :  0.6331334581613884\n",
            "Day :  2022 5 2\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  5.2\n",
            "-  Rainfall mean   :  0.7183673469387756\n",
            "-  Rainfall st dev :  1.3296728874069765\n",
            "Day :  2022 5 3\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  13.2\n",
            "-  Rainfall mean   :  3.0408163265306123\n",
            "-  Rainfall st dev :  3.3353345374142216\n",
            "Day :  2022 5 4\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  28.8\n",
            "-  Rainfall mean   :  4.969387755102041\n",
            "-  Rainfall st dev :  5.625459766484856\n",
            "Day :  2022 5 5\n",
            "-  Rainfall mini   :  0.8\n",
            "-  Rainfall max    :  37.0\n",
            "-  Rainfall mean   :  13.13673469387755\n",
            "-  Rainfall st dev :  7.731031926198363\n",
            "Day :  2022 5 6\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  28.2\n",
            "-  Rainfall mean   :  7.910204081632654\n",
            "-  Rainfall st dev :  7.501734354514137\n",
            "Day :  2022 5 7\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  27.2\n",
            "-  Rainfall mean   :  6.889795918367348\n",
            "-  Rainfall st dev :  6.731929550177628\n",
            "Day :  2022 5 8\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  17.2\n",
            "-  Rainfall mean   :  5.761224489795917\n",
            "-  Rainfall st dev :  4.221413648826929\n",
            "Day :  2022 5 9\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  12.4\n",
            "-  Rainfall mean   :  0.9489795918367349\n",
            "-  Rainfall st dev :  2.4528303059501666\n",
            "Day :  2022 5 10\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  5.4\n",
            "-  Rainfall mean   :  0.12653061224489798\n",
            "-  Rainfall st dev :  0.763123418002427\n",
            "Day :  2022 5 11\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.2\n",
            "-  Rainfall mean   :  0.004081632653061225\n",
            "-  Rainfall st dev :  0.028278380531736768\n",
            "Day :  2022 5 12\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  0.4\n",
            "-  Rainfall mean   :  0.00816326530612245\n",
            "-  Rainfall st dev :  0.05655676106347354\n",
            "Day :  2022 5 13\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  2.6\n",
            "-  Rainfall mean   :  0.05714285714285715\n",
            "-  Rainfall st dev :  0.3681171064778608\n",
            "Day :  2022 5 14\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  23.0\n",
            "-  Rainfall mean   :  2.5183673469387755\n",
            "-  Rainfall st dev :  4.457448408454754\n",
            "Day :  2022 5 15\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  8.0\n",
            "-  Rainfall mean   :  0.4959183673469388\n",
            "-  Rainfall st dev :  1.3480082334384098\n",
            "Day :  2022 5 16\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  3.0\n",
            "-  Rainfall mean   :  0.4224489795918367\n",
            "-  Rainfall st dev :  0.7287857685057453\n",
            "Day :  2022 5 17\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  1.2\n",
            "-  Rainfall mean   :  0.04081632653061224\n",
            "-  Rainfall st dev :  0.17603149920692102\n",
            "Day :  2022 5 18\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  15.4\n",
            "-  Rainfall mean   :  2.2816326530612248\n",
            "-  Rainfall st dev :  4.214875187315522\n",
            "Day :  2022 5 19\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  15.0\n",
            "-  Rainfall mean   :  1.006122448979592\n",
            "-  Rainfall st dev :  2.412755911203424\n",
            "Day :  2022 5 20\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  2.0\n",
            "-  Rainfall mean   :  0.09387755102040818\n",
            "-  Rainfall st dev :  0.34310001116191857\n",
            "Day :  2022 5 21\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  2.0\n",
            "-  Rainfall mean   :  0.05306122448979592\n",
            "-  Rainfall st dev :  0.28792150348650447\n",
            "Day :  2022 5 22\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  5.8\n",
            "-  Rainfall mean   :  0.4224489795918367\n",
            "-  Rainfall st dev :  0.9594354858345351\n",
            "Day :  2022 5 23\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  14.2\n",
            "-  Rainfall mean   :  1.430612244897959\n",
            "-  Rainfall st dev :  2.839387982462436\n",
            "Day :  2022 5 24\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  37.4\n",
            "-  Rainfall mean   :  6.671428571428571\n",
            "-  Rainfall st dev :  8.809897031991106\n",
            "Day :  2022 5 25\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  37.6\n",
            "-  Rainfall mean   :  3.6530612244897958\n",
            "-  Rainfall st dev :  7.769566532200646\n",
            "Day :  2022 5 26\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  8.0\n",
            "-  Rainfall mean   :  0.7081632653061225\n",
            "-  Rainfall st dev :  1.7142662778295628\n",
            "Day :  2022 5 27\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  8.0\n",
            "-  Rainfall mean   :  0.32244897959183677\n",
            "-  Rainfall st dev :  1.2505125979377416\n",
            "Day :  2022 5 28\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  22.9\n",
            "-  Rainfall mean   :  3.473469387755102\n",
            "-  Rainfall st dev :  5.8648999333880125\n",
            "Day :  2022 5 29\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  25.8\n",
            "-  Rainfall mean   :  3.010204081632653\n",
            "-  Rainfall st dev :  5.858997786344649\n",
            "Day :  2022 5 30\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  3.2\n",
            "-  Rainfall mean   :  0.1510204081632653\n",
            "-  Rainfall st dev :  0.5473954840058853\n",
            "Day :  2022 5 31\n",
            "-  Rainfall mini   :  0.0\n",
            "-  Rainfall max    :  3.2\n",
            "-  Rainfall mean   :  0.1510204081632653\n",
            "-  Rainfall st dev :  0.5473954840058853\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "# Perform 4\n",
        "\n",
        "Find the station that registered the highest temperature on each day (a station id for each day)."
      ],
      "metadata": {
        "id": "YJgszgkDL1rh"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "def Perform_4(b=a):   \n",
        "   #print(b.shape)\n",
        "\n",
        "   # Tuple to build an appropriate mask \n",
        "   ht = (0, 0, 0) \n",
        "   tt = (1, 0, 0,  0, 0, 0, 0, 0, 0, 0) * 49\n",
        "   rt = (ht + tt)\n",
        "    \n",
        "   ini = np.array(rt)\n",
        "   #print(ini.shape)\n",
        "   ini2 = np.zeros((151,1))\n",
        "   #print(ini2.shape)\n",
        "   ma = np.array(ini + ini2, dtype=bool) # broadcasting\n",
        "   #print(\"Mask shape : \", ma.shape)\n",
        "   #print(ma)\n",
        "\n",
        "   p4 = b[ma]  # Filter \n",
        "   p4n = (np.reshape(p4, (151, 49)))\n",
        "   #print(p4n)\n",
        "\n",
        "   for i in range(151):\n",
        "     print(\"Day : \", round(a[i][0]),round(a[i][1]),round(a[i][2])) \n",
        "     mm = np.max(p4n[i])\n",
        "     print(\"-  Temp max      : \",mm )\n",
        "     ft = (np.where(b[i] == mm))\n",
        "     print(\"In station id    : \",head[ft[0][0]][0][:6])\n",
        "     print(\"\")"
      ],
      "metadata": {
        "id": "gDwd0deMNuRR"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "Perform_4()"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "fkrQhVKGN1Au",
        "outputId": "cbe6e358-c66c-4751-ce94-b9e3c2dc7494"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Day :  2022 1 1\n",
            "-  Temp max      :  20.0\n",
            "In station id    :  WS_156\n",
            "\n",
            "Day :  2022 1 2\n",
            "-  Temp max      :  19.8\n",
            "In station id    :  WS_163\n",
            "\n",
            "Day :  2022 1 3\n",
            "-  Temp max      :  16.2\n",
            "In station id    :  WS_205\n",
            "\n",
            "Day :  2022 1 4\n",
            "-  Temp max      :  14.4\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 1 5\n",
            "-  Temp max      :  14.8\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 1 6\n",
            "-  Temp max      :  13.9\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 1 7\n",
            "-  Temp max      :  13.6\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 1 8\n",
            "-  Temp max      :  13.7\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 1 9\n",
            "-  Temp max      :  15.1\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 1 10\n",
            "-  Temp max      :  15.0\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 1 11\n",
            "-  Temp max      :  18.2\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 1 12\n",
            "-  Temp max      :  17.7\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 1 13\n",
            "-  Temp max      :  17.3\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 1 14\n",
            "-  Temp max      :  17.2\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 1 15\n",
            "-  Temp max      :  16.8\n",
            "In station id    :  WS_199\n",
            "\n",
            "Day :  2022 1 16\n",
            "-  Temp max      :  15.8\n",
            "In station id    :  WS_156\n",
            "\n",
            "Day :  2022 1 17\n",
            "-  Temp max      :  17.9\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 1 18\n",
            "-  Temp max      :  16.2\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 1 19\n",
            "-  Temp max      :  14.1\n",
            "In station id    :  WS_172\n",
            "\n",
            "Day :  2022 1 20\n",
            "-  Temp max      :  15.7\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 1 21\n",
            "-  Temp max      :  16.8\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 1 22\n",
            "-  Temp max      :  13.4\n",
            "In station id    :  WS_205\n",
            "\n",
            "Day :  2022 1 23\n",
            "-  Temp max      :  14.1\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 1 24\n",
            "-  Temp max      :  14.4\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 1 25\n",
            "-  Temp max      :  14.1\n",
            "In station id    :  WS_172\n",
            "\n",
            "Day :  2022 1 26\n",
            "-  Temp max      :  15.4\n",
            "In station id    :  WS_22_\n",
            "\n",
            "Day :  2022 1 27\n",
            "-  Temp max      :  15.2\n",
            "In station id    :  WS_21_\n",
            "\n",
            "Day :  2022 1 28\n",
            "-  Temp max      :  15.2\n",
            "In station id    :  WS_21_\n",
            "\n",
            "Day :  2022 1 29\n",
            "-  Temp max      :  15.8\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 1 30\n",
            "-  Temp max      :  23.2\n",
            "In station id    :  WS_156\n",
            "\n",
            "Day :  2022 1 31\n",
            "-  Temp max      :  17.6\n",
            "In station id    :  WS_163\n",
            "\n",
            "Day :  2022 2 1\n",
            "-  Temp max      :  17.7\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 2 2\n",
            "-  Temp max      :  22.2\n",
            "In station id    :  WS_199\n",
            "\n",
            "Day :  2022 2 3\n",
            "-  Temp max      :  20.6\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 2 4\n",
            "-  Temp max      :  14.9\n",
            "In station id    :  WS_21_\n",
            "\n",
            "Day :  2022 2 5\n",
            "-  Temp max      :  17.2\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 2 6\n",
            "-  Temp max      :  16.4\n",
            "In station id    :  WS_172\n",
            "\n",
            "Day :  2022 2 7\n",
            "-  Temp max      :  19.2\n",
            "In station id    :  WS_108\n",
            "\n",
            "Day :  2022 2 8\n",
            "-  Temp max      :  18.8\n",
            "In station id    :  WS_156\n",
            "\n",
            "Day :  2022 2 9\n",
            "-  Temp max      :  16.9\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 2 10\n",
            "-  Temp max      :  16.9\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 2 11\n",
            "-  Temp max      :  16.4\n",
            "In station id    :  WS_85_\n",
            "\n",
            "Day :  2022 2 12\n",
            "-  Temp max      :  16.6\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 2 13\n",
            "-  Temp max      :  14.2\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 2 14\n",
            "-  Temp max      :  12.2\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 2 15\n",
            "-  Temp max      :  13.0\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 2 16\n",
            "-  Temp max      :  15.1\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 2 17\n",
            "-  Temp max      :  24.0\n",
            "In station id    :  WS_29_\n",
            "\n",
            "Day :  2022 2 18\n",
            "-  Temp max      :  20.9\n",
            "In station id    :  WS_22_\n",
            "\n",
            "Day :  2022 2 19\n",
            "-  Temp max      :  17.1\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 2 20\n",
            "-  Temp max      :  19.1\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 2 21\n",
            "-  Temp max      :  18.3\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 2 22\n",
            "-  Temp max      :  19.8\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 2 23\n",
            "-  Temp max      :  22.2\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 2 24\n",
            "-  Temp max      :  21.8\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 2 25\n",
            "-  Temp max      :  18.2\n",
            "In station id    :  WS_199\n",
            "\n",
            "Day :  2022 2 26\n",
            "-  Temp max      :  17.4\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 2 27\n",
            "-  Temp max      :  17.2\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 2 28\n",
            "-  Temp max      :  16.8\n",
            "In station id    :  WS_199\n",
            "\n",
            "Day :  2022 3 1\n",
            "-  Temp max      :  16.5\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 3 2\n",
            "-  Temp max      :  14.6\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 3 3\n",
            "-  Temp max      :  17.4\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 3 4\n",
            "-  Temp max      :  17.2\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 3 5\n",
            "-  Temp max      :  14.9\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 3 6\n",
            "-  Temp max      :  13.9\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 3 7\n",
            "-  Temp max      :  14.4\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 3 8\n",
            "-  Temp max      :  12.7\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 3 9\n",
            "-  Temp max      :  15.9\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 3 10\n",
            "-  Temp max      :  16.8\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 3 11\n",
            "-  Temp max      :  17.2\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 3 12\n",
            "-  Temp max      :  15.0\n",
            "In station id    :  WS_65_\n",
            "\n",
            "Day :  2022 3 13\n",
            "-  Temp max      :  13.9\n",
            "In station id    :  WS_65_\n",
            "\n",
            "Day :  2022 3 14\n",
            "-  Temp max      :  14.7\n",
            "In station id    :  WS_21_\n",
            "\n",
            "Day :  2022 3 15\n",
            "-  Temp max      :  16.7\n",
            "In station id    :  WS_21_\n",
            "\n",
            "Day :  2022 3 16\n",
            "-  Temp max      :  17.4\n",
            "In station id    :  WS_21_\n",
            "\n",
            "Day :  2022 3 17\n",
            "-  Temp max      :  17.9\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 3 18\n",
            "-  Temp max      :  18.8\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 3 19\n",
            "-  Temp max      :  18.1\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 3 20\n",
            "-  Temp max      :  18.2\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 3 21\n",
            "-  Temp max      :  16.4\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 3 22\n",
            "-  Temp max      :  18.3\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 3 23\n",
            "-  Temp max      :  21.9\n",
            "In station id    :  WS_199\n",
            "\n",
            "Day :  2022 3 24\n",
            "-  Temp max      :  24.8\n",
            "In station id    :  WS_156\n",
            "\n",
            "Day :  2022 3 25\n",
            "-  Temp max      :  23.9\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 3 26\n",
            "-  Temp max      :  23.3\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 3 27\n",
            "-  Temp max      :  24.1\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 3 28\n",
            "-  Temp max      :  24.3\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 3 29\n",
            "-  Temp max      :  24.1\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 3 30\n",
            "-  Temp max      :  21.6\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 3 31\n",
            "-  Temp max      :  19.4\n",
            "In station id    :  WS_199\n",
            "\n",
            "Day :  2022 4 1\n",
            "-  Temp max      :  18.2\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 2\n",
            "-  Temp max      :  18.2\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 3\n",
            "-  Temp max      :  18.2\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 4\n",
            "-  Temp max      :  18.2\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 5\n",
            "-  Temp max      :  23.4\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 4 6\n",
            "-  Temp max      :  20.7\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 7\n",
            "-  Temp max      :  25.9\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 4 8\n",
            "-  Temp max      :  25.9\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 4 9\n",
            "-  Temp max      :  25.9\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 4 10\n",
            "-  Temp max      :  25.9\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 4 11\n",
            "-  Temp max      :  20.8\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 12\n",
            "-  Temp max      :  21.8\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 4 13\n",
            "-  Temp max      :  23.3\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 14\n",
            "-  Temp max      :  26.3\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 4 15\n",
            "-  Temp max      :  30.6\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 4 16\n",
            "-  Temp max      :  28.7\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 17\n",
            "-  Temp max      :  27.7\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 18\n",
            "-  Temp max      :  22.9\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 4 19\n",
            "-  Temp max      :  24.3\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 4 20\n",
            "-  Temp max      :  23.3\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 21\n",
            "-  Temp max      :  20.3\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 4 22\n",
            "-  Temp max      :  20.8\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 23\n",
            "-  Temp max      :  20.8\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 24\n",
            "-  Temp max      :  18.8\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 25\n",
            "-  Temp max      :  22.8\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 4 26\n",
            "-  Temp max      :  23.8\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 27\n",
            "-  Temp max      :  26.0\n",
            "In station id    :  WS_199\n",
            "\n",
            "Day :  2022 4 28\n",
            "-  Temp max      :  26.3\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 4 29\n",
            "-  Temp max      :  25.2\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 4 30\n",
            "-  Temp max      :  24.8\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 5 1\n",
            "-  Temp max      :  25.4\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 5 2\n",
            "-  Temp max      :  25.6\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 5 3\n",
            "-  Temp max      :  25.6\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 5 4\n",
            "-  Temp max      :  24.9\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 5 5\n",
            "-  Temp max      :  22.6\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 5 6\n",
            "-  Temp max      :  20.3\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 5 7\n",
            "-  Temp max      :  24.6\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 5 8\n",
            "-  Temp max      :  25.8\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 5 9\n",
            "-  Temp max      :  26.8\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 5 10\n",
            "-  Temp max      :  29.4\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 5 11\n",
            "-  Temp max      :  30.7\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 5 12\n",
            "-  Temp max      :  30.6\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 5 13\n",
            "-  Temp max      :  29.6\n",
            "In station id    :  WS_151\n",
            "\n",
            "Day :  2022 5 14\n",
            "-  Temp max      :  31.1\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 5 15\n",
            "-  Temp max      :  30.6\n",
            "In station id    :  WS_151\n",
            "\n",
            "Day :  2022 5 16\n",
            "-  Temp max      :  31.8\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 5 17\n",
            "-  Temp max      :  33.5\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 5 18\n",
            "-  Temp max      :  32.8\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 5 19\n",
            "-  Temp max      :  32.0\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 5 20\n",
            "-  Temp max      :  32.0\n",
            "In station id    :  WS_154\n",
            "\n",
            "Day :  2022 5 21\n",
            "-  Temp max      :  34.8\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 5 22\n",
            "-  Temp max      :  33.3\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 5 23\n",
            "-  Temp max      :  30.7\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 5 24\n",
            "-  Temp max      :  29.5\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 5 25\n",
            "-  Temp max      :  29.3\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 5 26\n",
            "-  Temp max      :  30.3\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 5 27\n",
            "-  Temp max      :  34.0\n",
            "In station id    :  WS_92_\n",
            "\n",
            "Day :  2022 5 28\n",
            "-  Temp max      :  33.3\n",
            "In station id    :  WS_111\n",
            "\n",
            "Day :  2022 5 29\n",
            "-  Temp max      :  32.9\n",
            "In station id    :  WS_28_\n",
            "\n",
            "Day :  2022 5 30\n",
            "-  Temp max      :  23.3\n",
            "In station id    :  WS_124\n",
            "\n",
            "Day :  2022 5 31\n",
            "-  Temp max      :  23.3\n",
            "In station id    :  WS_124\n",
            "\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "# Perform 5\n",
        "\n",
        "Find the day in which the lowest temperature has been recorded."
      ],
      "metadata": {
        "id": "sB66qQI-Alnh"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "def Perform_5(b=a_n, min_tem = 'WS_1_MIN_T'):    \n",
        "\n",
        "   p5 = b[min_tem]\n",
        "   #Descriptive statistics\n",
        "   mm = np.min(p5)\n",
        "   print(\" Temp Min : \", mm)\n",
        "   ft = (np.where(b[min_tem] == mm))\n",
        "   print(\" Recorded in station : \", head[ft[0][0]][0][:6])\n",
        "   print(\" Day : \", round(b[ft[0][0]][0]),round(b[ft[0][0]][1]),round(b[ft[0][0]][2])) \n"
      ],
      "metadata": {
        "id": "4t10O-6RBZE_"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "Perform_5()"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "0lGr67xRB9e8",
        "outputId": "551d9a74-adf7-40ec-a4b1-cd0671ab42eb"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            " Temp Min :  -3.2\n",
            " Recorded in station :  WS_130\n",
            " Day :  2022 1 22\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "# Perform 6\n",
        "Find the average relative humidity recorded in February for each station."
      ],
      "metadata": {
        "id": "J33fgzedDPaj"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "def Perform_6(b=a):    \n",
        "\n",
        "   # Tuple to build an appropriate mask \n",
        "   ht = (0, 0, 0) \n",
        "   tt = (0, 0, 0, 0, 0, 1, 0, 0, 0, 0) * 49\n",
        "   rt = (ht + tt)\n",
        "    \n",
        "   ini = np.array(rt)\n",
        "   #print(ini.shape)\n",
        "   ini2 = np.zeros((151,1))\n",
        "   #print(ini2.shape)\n",
        "   ma = np.array(ini + ini2, dtype=bool) # broadcasting\n",
        "   #print(\"Mask shape : \", ma.shape)\n",
        "   #print(ma)\n",
        "\n",
        "   p6 = b[ma]  # Filter \n",
        "   p6n = (np.reshape(p6, (151, 49)))\n",
        "   rt,sp6n,rt2 = np.vsplit(p6n, [31,59])  # we need only February : sp6n \n",
        "   #print(sp6n)\n",
        "   #print(sp6n.shape)\n",
        "   mea = np.mean(sp6n, axis=0)  # average relative humidity recorded for each station.\n",
        "   \n",
        "   c = 0\n",
        "   for i in range(8,490,10):\n",
        "     print(\"Station : \", head[i][0][:6] )\n",
        "     print(\"Average humidity recorded in February : \", mea[c],\"\\n\")\n",
        "     c += 1\n"
      ],
      "metadata": {
        "id": "XJZa4hiXDWkI"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "Perform_6()"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "fRxgU5G1e4XI",
        "outputId": "84363977-02ed-4f1d-fd71-4adc0545f5d5"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Station :  WS_1_M\n",
            "Average humidity recorded in February :  61.48214285714285 \n",
            "\n",
            "Station :  WS_130\n",
            "Average humidity recorded in February :  74.27142857142857 \n",
            "\n",
            "Station :  WS_5_M\n",
            "Average humidity recorded in February :  57.14999999999999 \n",
            "\n",
            "Station :  WS_10_\n",
            "Average humidity recorded in February :  59.76071428571429 \n",
            "\n",
            "Station :  WS_140\n",
            "Average humidity recorded in February :  60.271428571428565 \n",
            "\n",
            "Station :  WS_142\n",
            "Average humidity recorded in February :  72.52142857142857 \n",
            "\n",
            "Station :  WS_143\n",
            "Average humidity recorded in February :  75.45714285714287 \n",
            "\n",
            "Station :  WS_145\n",
            "Average humidity recorded in February :  47.66428571428572 \n",
            "\n",
            "Station :  WS_21_\n",
            "Average humidity recorded in February :  60.01071428571429 \n",
            "\n",
            "Station :  WS_22_\n",
            "Average humidity recorded in February :  28.582142857142863 \n",
            "\n",
            "Station :  WS_151\n",
            "Average humidity recorded in February :  75.7357142857143 \n",
            "\n",
            "Station :  WS_24_\n",
            "Average humidity recorded in February :  56.77142857142858 \n",
            "\n",
            "Station :  WS_152\n",
            "Average humidity recorded in February :  61.30714285714286 \n",
            "\n",
            "Station :  WS_154\n",
            "Average humidity recorded in February :  73.16428571428573 \n",
            "\n",
            "Station :  WS_28_\n",
            "Average humidity recorded in February :  70.16785714285713 \n",
            "\n",
            "Station :  WS_29_\n",
            "Average humidity recorded in February :  57.03928571428572 \n",
            "\n",
            "Station :  WS_156\n",
            "Average humidity recorded in February :  48.50000000000001 \n",
            "\n",
            "Station :  WS_158\n",
            "Average humidity recorded in February :  65.90357142857142 \n",
            "\n",
            "Station :  WS_159\n",
            "Average humidity recorded in February :  48.017857142857146 \n",
            "\n",
            "Station :  WS_161\n",
            "Average humidity recorded in February :  52.607142857142854 \n",
            "\n",
            "Station :  WS_163\n",
            "Average humidity recorded in February :  53.853571428571435 \n",
            "\n",
            "Station :  WS_37_\n",
            "Average humidity recorded in February :  55.164285714285704 \n",
            "\n",
            "Station :  WS_42_\n",
            "Average humidity recorded in February :  68.425 \n",
            "\n",
            "Station :  WS_43_\n",
            "Average humidity recorded in February :  52.314285714285695 \n",
            "\n",
            "Station :  WS_170\n",
            "Average humidity recorded in February :  62.117857142857154 \n",
            "\n",
            "Station :  WS_172\n",
            "Average humidity recorded in February :  70.35 \n",
            "\n",
            "Station :  WS_174\n",
            "Average humidity recorded in February :  75.42857142857144 \n",
            "\n",
            "Station :  WS_176\n",
            "Average humidity recorded in February :  53.63214285714286 \n",
            "\n",
            "Station :  WS_179\n",
            "Average humidity recorded in February :  62.092857142857135 \n",
            "\n",
            "Station :  WS_187\n",
            "Average humidity recorded in February :  57.77499999999999 \n",
            "\n",
            "Station :  WS_191\n",
            "Average humidity recorded in February :  64.64285714285715 \n",
            "\n",
            "Station :  WS_192\n",
            "Average humidity recorded in February :  47.4 \n",
            "\n",
            "Station :  WS_65_\n",
            "Average humidity recorded in February :  72.90714285714287 \n",
            "\n",
            "Station :  WS_193\n",
            "Average humidity recorded in February :  56.70714285714287 \n",
            "\n",
            "Station :  WS_67_\n",
            "Average humidity recorded in February :  69.57142857142858 \n",
            "\n",
            "Station :  WS_194\n",
            "Average humidity recorded in February :  55.86428571428571 \n",
            "\n",
            "Station :  WS_195\n",
            "Average humidity recorded in February :  65.11428571428573 \n",
            "\n",
            "Station :  WS_198\n",
            "Average humidity recorded in February :  55.55714285714287 \n",
            "\n",
            "Station :  WS_199\n",
            "Average humidity recorded in February :  58.042857142857144 \n",
            "\n",
            "Station :  WS_72_\n",
            "Average humidity recorded in February :  62.07142857142858 \n",
            "\n",
            "Station :  WS_202\n",
            "Average humidity recorded in February :  64.37142857142857 \n",
            "\n",
            "Station :  WS_204\n",
            "Average humidity recorded in February :  55.36785714285715 \n",
            "\n",
            "Station :  WS_205\n",
            "Average humidity recorded in February :  44.353571428571435 \n",
            "\n",
            "Station :  WS_85_\n",
            "Average humidity recorded in February :  66.79642857142856 \n",
            "\n",
            "Station :  WS_86_\n",
            "Average humidity recorded in February :  59.29642857142858 \n",
            "\n",
            "Station :  WS_92_\n",
            "Average humidity recorded in February :  70.26071428571429 \n",
            "\n",
            "Station :  WS_108\n",
            "Average humidity recorded in February :  57.35714285714285 \n",
            "\n",
            "Station :  WS_111\n",
            "Average humidity recorded in February :  53.625 \n",
            "\n",
            "Station :  WS_124\n",
            "Average humidity recorded in February :  66.54285714285716 \n",
            "\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "# Perform 7\n",
        "Calculate the average relative humidity for each month in the dataset."
      ],
      "metadata": {
        "id": "UsadFBGpDW8h"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "def Perform_7(b=a):    \n",
        "\n",
        "   # Tuple to build an appropriate mask \n",
        "   ht = (0, 0, 0) \n",
        "   tt = (0, 0, 0, 0, 0, 1, 0, 0, 0, 0) * 49\n",
        "   rt = (ht + tt)\n",
        "    \n",
        "   ini = np.array(rt)\n",
        "   #print(ini.shape)\n",
        "   ini2 = np.zeros((151,1))\n",
        "   #print(ini2.shape)\n",
        "   ma = np.array(ini + ini2, dtype=bool) # broadcasting\n",
        "   #print(\"Mask shape : \", ma.shape)\n",
        "   #print(ma)\n",
        "\n",
        "   p7 = b[ma]  # Filter \n",
        "   p7n = (np.reshape(p7, (151, 49)))\n",
        "   gen,feb,marz,apr,mag = np.vsplit(p7n, [31,59,90,120])  # split in month \n",
        "   #print(mag)\n",
        "   #print(mag.shape)\n",
        "   print(\"  Average relative humidity for each month\")\n",
        "   print(f\"January : {np.mean(gen)}  ,  February : {np.mean(feb)} , \\nMarch : {np.mean(marz)} , April  : {np.mean(apr)} , May  :  {np.mean(mag)}\")\n",
        "   #mea = np.mean(mag)\n",
        "   \n"
      ],
      "metadata": {
        "id": "5u7s38jXo9dR"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "Perform_7()"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "5LiAMNCHqwWd",
        "outputId": "9d120390-254a-43a0-8602-c878cc4f46e0"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "  Average relative humidity for each month\n",
            "January : 63.95385121790651  ,  February : 60.477332361516034 , \n",
            "March : 58.27432521395656 , April  : 64.35013605442177 , May  :  70.95602369980251\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "# Perform 8\n",
        "Compute the monthly average of daily temperature ranges across all stations. Visualize them as a bar plot.\n"
      ],
      "metadata": {
        "id": "9hpDTyF0DdqD"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "def Perform_8(b=a):    \n",
        "\n",
        "   # Tuple to build an appropriate mask \n",
        "   ht = (0, 0, 0) \n",
        "   tt = (0, 0, 1, 0, 0, 0, 0, 0, 0, 0) * 49\n",
        "   rt = (ht + tt)\n",
        "    \n",
        "   ini = np.array(rt)\n",
        "   #print(ini.shape)\n",
        "   ini2 = np.zeros((151,1))\n",
        "   #print(ini2.shape)\n",
        "   ma = np.array(ini + ini2, dtype=bool) # broadcasting\n",
        "   #print(\"Mask shape : \", ma.shape)\n",
        "   #print(ma)\n",
        "\n",
        "   p8 = b[ma]  # Filter \n",
        "   p8n = (np.reshape(p8, (151, 49)))\n",
        "   gen,feb,marz,apr,mag = np.vsplit(p8n, [31,59,90,120])  # split in month \n",
        "\n",
        "   # Monthly Average relative daily temperature ranges across all stations.\n",
        "   mea_gen = np.mean(gen, axis=1) \n",
        "   m_gen = np.mean(mea_gen)\n",
        "   mea_feb = np.mean(feb, axis=1) \n",
        "   m_feb = np.mean(mea_feb)\n",
        "   mea_marz = np.mean(marz, axis=1) \n",
        "   m_marz = np.mean(mea_marz)\n",
        "   mea_apr = np.mean(apr, axis=1) \n",
        "   m_apr = np.mean(mea_apr)\n",
        "   mea_mag = np.mean(mag, axis=1) \n",
        "   m_mag = np.mean(mea_mag)\n",
        "   #print(mea_gen)\n",
        "   #print(mea_gen.shape)\n",
        "   #print(m_gen)\n",
        "   print(\"Monthly average of daily temperature ranges across all stations.\\n\")\n",
        "   print(f\"  January : {mea_gen} ,\\nTot mean January : {m_gen},\\n\\n  February : {mea_feb},\\nTot mean February : {m_feb},\\n\") \n",
        "   print(f\" \\nMarch : {mea_marz} , \\nTot mean March  : {m_marz},\\n\\n April  : {mea_apr}, \\nTot mean April  : {m_apr} ,\\n\\n May  :  {mea_mag}, \\nTot mean May  : {m_mag}\\n \")\n",
        "\n",
        "\n",
        "   fig = plt.figure(figsize = (8, 4))\n",
        "   mon = [\"January\", \"February\", \"March\", \"April\", \"May\"]\n",
        "   val = [m_gen,m_feb,m_marz,m_apr,m_mag]\n",
        "   # creating the bar plot\n",
        "   plt.bar(mon, val, color ='green', width = 0.5)\n",
        " \n",
        "   plt.xlabel(\"Month\",color ='green')\n",
        "   plt.ylabel(\"Temperature\",color ='green')\n",
        "   plt.title(\"Monthly average of daily temperature ranges across all stations\",color ='green')\n",
        "   plt.show()\n"
      ],
      "metadata": {
        "id": "MDhZd8tgsAK2"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "Perform_8()"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 1000
        },
        "id": "Q4roE9ObvfFO",
        "outputId": "29abc4cc-6fbe-4c12-da3b-d7efa9efc098"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Monthly average of daily temperature ranges across all stations.\n",
            "\n",
            "  January : [8.5755102  6.14693878 5.51632653 5.71428571 5.67346939 2.54897959\n",
            " 0.60612245 0.71428571 1.85714286 1.9122449  2.23469388 1.03265306\n",
            " 1.52857143 4.5122449  5.8122449  4.38979592 5.22244898 4.95306122\n",
            " 2.35510204 3.08163265 2.16530612 2.2244898  3.25918367 2.15102041\n",
            " 0.76122449 1.67346939 2.4        2.4        6.39795918 7.41632653\n",
            " 5.7755102 ] ,\n",
            "Tot mean January : 3.5810401579986837,\n",
            "\n",
            "  February : [5.73877551 9.66938776 7.46938776 5.31020408 5.51836735 3.49387755\n",
            " 7.44285714 5.7244898  6.41632653 5.78163265 6.12040816 3.53061224\n",
            " 2.57346939 1.68163265 3.07346939 3.64897959 8.79387755 7.0122449\n",
            " 7.56734694 7.83061224 7.85714286 8.38571429 8.47346939 6.51020408\n",
            " 6.70612245 5.4        3.31632653 4.44081633],\n",
            "Tot mean February : 5.91027696793003,\n",
            "\n",
            " \n",
            "March : [ 3.31020408  3.99591837  6.36530612  5.89183673  4.24081633  2.91836735\n",
            "  2.54693878  1.89591837  3.76938776  5.76734694  4.81632653  3.96326531\n",
            "  4.21020408  5.6877551   7.05714286  8.69387755  8.60816327  9.1244898\n",
            "  8.47959184  6.94693878  6.70816327  6.81020408  9.03469388 10.95510204\n",
            " 10.94693878 11.20204082 11.79183673 12.23673469 11.08571429  7.98571429\n",
            "  9.16734694] , \n",
            "Tot mean March  : 6.974654377880184,\n",
            "\n",
            " April  : [ 6.40612245  6.40612245  6.40612245  6.40612245 12.60408163  9.59591837\n",
            " 11.13673469 11.13673469 11.13673469 11.13673469  8.73265306 11.04693878\n",
            " 13.2877551  14.96734694 16.48367347 16.55102041 12.76938776 11.43061224\n",
            " 12.17755102  9.56938776  9.61632653 11.51020408  9.91020408  9.26734694\n",
            " 11.72040816 12.32857143 13.75102041 15.06122449 14.80816327 13.76734694], \n",
            "Tot mean April  : 11.370952380952378 ,\n",
            "\n",
            " May  :  [13.98571429 14.23061224 13.10204082 13.4244898  11.14285714 11.27142857\n",
            " 13.43877551 13.81836735 16.08979592 17.53877551 18.81428571 19.08979592\n",
            " 19.05918367 19.74897959 19.63265306 19.47346939 20.68979592 21.37959184\n",
            " 20.5755102  20.95306122 21.75918367 21.38571429 20.36530612 18.26734694\n",
            " 17.35918367 19.68367347 21.9        21.89387755 14.38979592 14.26938776\n",
            " 14.26938776], \n",
            "Tot mean May  : 17.516194865042788\n",
            " \n"
          ]
        },
        {
          "output_type": "display_data",
          "data": {
            "text/plain": [
              "<Figure size 800x400 with 1 Axes>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAArwAAAGJCAYAAABo5eDAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABPXElEQVR4nO3deVRU9f/H8dcAMiiyqKBA4r4vuJZLlpqmopmaqZHlbt/6ptnX9Fe2KJhlu1ZuLSqWu2VqpaaRa2pukUtqbogbKi4gqIBwf394mBxZZGQQGJ+Pc+Yc5nM/9973nXvn8uLOZy4mwzAMAQAAAA7KKb8LAAAAAPISgRcAAAAOjcALAAAAh0bgBQAAgEMj8AIAAMChEXgBAADg0Ai8AAAAcGgEXgAAADg0Ai8AAAAcGoEXd50pzKQhy4fctl94ZLhMYSZFXYrK+6Jgs4PnD6rdt+3k9Z6XTGEmLdm/xKb510atlSnMpLVRa21ed+jaUJnCTFZtFSZWUL8l/WxeFgDHk9nvj1bhrdQqvFW+1ZRTpjCTQteG5ncZDofA60DS3+CmMJM2Rm/MMN0wDAVOCJQpzKTH5j6Wp7VsOr5JoWtDdenapTxdD/JP3yV9tfvsbr3zyDv6ttu3ahzQOL9LsptTl08pdG2oImMi87uUAof3Nu519jg/LD+4nFB7lxF4HZCbi5vm7p6boX3dsXU6EX9CZmdzntew6fgmha0L45eig7qaclWbT2zWwAYDNeSBIXom6BmV9Sx719b/5sNv6uobV/Ns+acun1LYujACbyZ4b+NeZ4/zw/KDyxW2LizTaVffuKo3H37zjpeNzBF4HVDHqh216O9Fup523ap97u65auTfSH7F/fKpMqRLTE7M7xJy5dyVc5IkbzfvfFm/i5OL3Fzc8mXdjqagHIt5VUdB2b7CLM1I07Xr1/K7jHuGm4ubXJxc8rsMh0PgdUAhdUJ0/sp5rT682tKWnJqs7/7+Tk/XfTrTeRKTE/XKL68ocEKgzOPMqj6puj7a9JEMw7Dqlz7+dsn+JaozpY7M48yqPaW2Vh5aaekTujZUI1ePlCRV/LSiZZjFrWNxs1tGZvou6SufD3yUkpqSYVq7b9up+qTq2c6/4dgG9VjUQ+UmlJN5nFmBEwL1v5X/09WUf68UfrTpI5nCTDp26ViG+Uf9Okqub7vq4tWLlrY/TvyhDrM7yOs9LxV7p5hahrfU79G/W82XPt7073N/6+nvn1aJ90uoxcwWkqRdZ3ap35J+qvRpJbmNc5PfR34asHSAzl85n2H9a6PWqvGXjeU2zk2VP6usL7Z/kelYVkmavWu2Gn3ZSEXfKaqS75fUU989peNxx7N9fdL9efpPBc8Jlud4TxV/t7jafNNGW05ssdqe8hPLS5JGrh4pU5hJFSZWyHaZJ+JPqOv8rnJ/112lPyyt/638n5KuJ2Xol5N9lF5DZtud7sjFIzKFmTRh84QM0zYd3yRTmEnzds/LdN61UWt1/1f3S5L6L+1vOX7DI8MtfWzZ7/+c/0fPLH5GXu95yfdDX73121syDEPH446ry/wu8hzvKb+P/PTxpo8z1GEKM2nBngV6PeJ1+X3kJ/d33fX4vMcz3Zd341jM7r0ddSkqw+uU7tYxidnVId358WuP91r6Mg5dOKR+S/rJ+z1veb3npf5L++tKyhWrvldTruqlFS/J5wMfeYz30OPzHtfJ+JOZjsE8GX9SA5YOUJmPyljOeTP+nJFhGz7/43PVnlJbxd4pphLvl1DjLxtn+ondzZJTkzV6zWg1+rKRvN7zkvu77npo5kNac3RNhr5pRpo+3fKp6k6tK7dxbvL90FcdZnfQ9lPbLX3Sz/Nzds1R7Sm1ZR5ntpyfb3d+kKSU1BSFrQ1T1c+rym2cm0p9UEotZrSw+p0UkxCj/kv7q+wnZWUeZ5b/x/7qMr/Lbb+zYcs5806tPrxaLWa0kPd73ir+bnFVn1Rdr0e8Lun254ecnMP6LemnydsmS5Jl/pvPZ5kdPzl53dOHNf4e/buG/zJcvh/6yv1dd3Vb0E3nEs9Z9d1+arvaz24vnw98VPSdoqr4aUUNWDrALq9fQcWfEA6ogncFNQtspnl75im4arAkacXBFYpLitNTdZ7SZ398ZtXfMAw9Pv9xrTm6RgMbDFR9v/r65fAvGrl6pE7Gn9SEDtahYWP0Ri3et1j/vf+/8nD10GdbP1P3hd0V/XK0ShUrpSdqPqF/zv+jeXvmaUL7CfIp5iNJ8i3mm+NlZObZoGf1zV/f6JfDv+ixav+OQY5JiNFvR3/TmJZjsn1dFv29SFdSruiFxi+oVLFS2npyqz7f+rlOXD6hRT0WSZJ61u6p/1v9f1q4d6FGPjjSav6Ffy9Uu8rtVKJoCUnSb0d/U/CcYDXyb6QxLcfIyeSkmZEz9cg3j2hD/w164L4HrObvsaiHqpasqncfeVeGbvwhsfrwah25eET96/eXX3E/7T23V1/u+FJ7z+3VloFbZDLdOAn+efpPdZjdQf4e/gprFaZUI1Vj14+1ek3TvbP+Hb215i31rN1TgxoM0rkr5/T51s/1cPjD+vM/f2Z7VXbv2b16aOZD8jR76v8e/D8VcSqiL3Z8oVbhrbSu3zo1KdtET9R8Qt5u3vrfL/9TSJ0QdazaUcVdi2e5zKspV9XmmzaKjovWSw+8pACPAH2761v9FvXbHe2jnKhUopIeDHxQc3bP0f+a/c9q2pxdc+Th6qEuNbpkOm9Nn5oa22qsRq8drecaPqeHyj8kSWoe2FyS7fu913e9VNOnpt5r855+Pvizxm0Yp5JFS+qLHV/okYqP6P2272vO7jkasXqE7r/vfj1c/mGr+d/Z8I5MJpNeffBVnU08q4l/TFTbb9sq8j+RKlqk6B3VdKfHYnbv7fSr/rbIrI7cHL+53b6b9VzUUxVLVNT4NuO18/ROff3n1ypdrLTef/R9S59+S/tp4d6FejboWTUt21Trjq1Tp7mdMtRzJuGMmk5vKpNMGnL/EPm6+2rFoRUauGyg4pPi9XLTlyVJX+34Si+tfElP1npSw5oM07Xr17TrzC79ceKPLC9WSFJ8Ury+3vm1QuqEaHDDwbqcdFnT/5yu9rPba+vgrarvV9/Sd+CygQqPDFdwlWANajBI19Oua0P0Bm05scVqHP5vR3/Twr0LNeSBIfIp5qMK3hVydH6QbvzRMH7jeA1qOEgP3PeA4pPitf3Udu08vVOPVn5UktR9YXftPbtXQx8YqgreFXQ28axWH1mt6LhoVfCukOW22rofbbX37F49Nu8xBZUJ0tjWY2V2NuvQhUP6/fiNPx5vd37IyTnsP43+o1OXT2n1kdX6ttu3OaopJ697uqErhqpE0RIa03KMoi5FaeKWiRriPEQLnlwgSTqbeFbtvm0nX3dfvdbiNXm7eSvqUpQW71ucq9euwDPgMGb+OdNQqIxtJ7cZk/6YZHi862FcSb5iGIZh9FjYw2gd3towDMMoP6G80WlOJ8t8S/YtMRQqY9y6cVbLe3Lhk4Yp1GQcOn/I0qZQGa5vu1q1/RXzl6FQGZ//8bml7cPfPzQUKuPoxaMZ6szpMtK3J30ZqWmpRtlPyhq9FvWyWt4nmz4xTKEm48iFI9m+Pumvxc3GbxhvmEJNxrFLxyxtzb5uZjT6opFVv60nthoKlfFN5DeGYRhGWlqaUfWzqkb7b9sbaWlpVuuoOLGi8eg3j1raxqwZYyhURsh3ITmqad7ueYZCZayPWm9p6zy3s1HsnWLGyfiTlraD5w8aLmNdDIX++zaOuhhlOIc5G++sf8dqmbvP7DZcxrpkaL9V1/ldDde3XY3DFw5b2k7FnzI83vUwHp75sKXt6MWjhkJlfPj7h9kuzzAMY+LmiYZCZSzcs9DSlpicaFT5rIqhUBlrjq6xtOd0H6W/pjcrP6G80feHvpbnX2z/wlCojH3n9lnakq8nGz4f+Fj1y8y2k9sMhcqY+edMq/Y72e/PLXvO0nY99bpR9pOyhinUZLy34T1L+8WrF42i44pa1bXm6BpDoTLu+/g+I/5avKV94Z6FhkJlfLrl0zuuKTfHYlbv7fRj4tbXzDBuvOfHrBlz2zpye/zaY/vSlzFgyQCrvt3mdzNKvV/K8nzHqR2GQmW8vOJlq379lvTLsL0Dlw40/D/yN2ITY636PvXdU4bXeC9LbV3mdTFqT66d7TZm5nrqdSPpepJV28WrF40yH5ax2o7fjvxmKFTGS8tfyrCMm48dhcpwCnMy9p7da9Unp+eHelPrWf2OudXFqxdzfP64VU73462/PwzDMFrObGm0nNky2+VP2DzBUKiMc4nnsuyT1fkhq/oyO4e9+POLGc5h6W49fnL6uqdvc9tv2lrtz/+t/J/hHOZsXLp6yTAMw/hh3w+WrHAvYUiDg+pZu6euXr+qn/75SZeTLuunf37K8grB8oPL5Wxy1ktNXrJqf6XZKzJkaMWhFVbtbSu1VeWSlS3Pg8oEydPsqSMXj+S4vjtZhpPJSb3r9tayA8t0OemypX3O7jlqHthcFUtUzHad6VfDpBtDOGKvxKp5YHMZMvTn6T8t03rV7qUdp3fo8IXDlrYFexfI7Gy2XBWMjInUwQsH9XTdp3X+6nnFXolV7JVYJaYkqk3FNlp/bL3SjDSr9T/f+Plsa7p2/Zpir8SqadmmkqSdp3dKklLTUvXrkV/VtUZXBXgEWPpXKVlFwVWCrZa3eN9ipRlp6lm7p6Wm2Cux8ivup6olq2pNVMaPONOlpqVq1eFV6lqjqyqVqGRp9/fw19N1n9bG6I2KT4rPcv6sLD+0XP7F/fVkrSctbcWKFNNzDZ/L0Den+ygnetbuKTcXN83ZNcfS9svhXxR7JVbPBD1j83ZId7bfBzUcZPnZ2clZjQMay5ChgQ0HWtq93bxV3ad6psd/n3p95GH2sDx/staT8i/ur+UHl99xTXd6LNrbrXXk5vjNbrmS7dt36zIeKveQzl89b3kPpH/E/9/7/2vVb+gDQ62eG4ah7/d9r87VOsuQYbVd7Su3V1xSnGX93m7eOhF/QttObsvRdqZzdnKWq7OrpBtDFi5cvaDradfVOKCxdsb8u23f7/teJpk0plXGT8NuvTLasnxL1fKtZXluy/nB281be8/t1cHzBzOtt6hLUbk6u2pt1FqrIWI5kdfHafonCEv3L83wvrG1vtyew6Q7Oy8/1/A5q/35ULmHlGqk6ljcjaF66dv40z8/ZTpE0FExpMFB+br7qm2ltpq7Z66upFxRqpFqFThudizumAI8Aqx+qUo3PrqRlGE8aznPchmWUcKthC5ey/mJ606X0adeH73/+/v6Yf8P6lOvjw7EHtCO0zs0rdO0264zOi5ao9eM1rIDyzKsJy4pzvJzj9o9NHzVcC3Yu0CvP/S6DMPQor8XKbhqsDzNnpKkgxdunMj7Lumb5frirsVZhj9IUkXvjIH8wtULClsbpvl75+ts4tlMazqbeFZXr19VlRJVMsxfpaR128ELB2XIUNXPq2ZaUxHnIlnWe+7KOV1JuaLqpTKOha7pU1NpRpqOxx1X7dK1s1xGZo5dOqYqJatk+IVa3SfjenK6j3LC281bnat11tw9c/X2I29LuvHH0X0e9+mRio/YtKx0d7Lfy3lZH+teZi+5ubhZhgPc3H7+asZxiFVLWu9Lk8mkKiWrWMY63s1j0d5urSM3x292y5Vs375b91v663fx6kV5mj117NIxOZmcMvyhfet78tyVc7p07ZK+3Pmlvtz5Zab1ptfz6oOv6tcjv+qBrx9QlZJV1K5SOz1d92k9WO7B22yxNCtylj7e/LH2x+5XStq/Iebm1+LwxcMK8AhQyaIlb7u8W19DW84PY1uPVZf5XVRtUjXVKV1HHSp30LP1nlVQmSBJktnFrPfbvq9XVr2iMh+VUdOyTfVYtcfUp16f236pOq+P0161e+nrnV9r0I+D9FrEa2pTsY2eqPmEnqz1pJxMt79GaM9zmHRn5+Xsjl3pxh8z3Wt2V9i6ME3YMkGtKrRS1+pd9XTdp2V2yfu7OOUXAq8De7rO0xr842DFJMQouEqw3b5R7+zknGm7ccsX3PJiGbV8a6mRfyPN3jVbfer10exds+Xq7KqetXtmO19qWqoe/fZRXbh6Qa8++Kpq+NSQu6u7TsafVL+l/az+kg/wCNBD5R7Swr0L9fpDr2vLiS2KjovW+23/HbuX3v/DRz+0Gh93s1vHtd78l3+6not6atPxTRrZfKTq+9VXcdfiSjPS1GFOhzu6upBmpMkkk1b0XpHpa5zdWNv8Zss+yqk+9fpo0d+LtOn4JtUtXVfLDizTfxv/N0e/uDJzJ/s9s/3gbMr9eyg3NeXFsWhS5mMnU9NSs5zn1jrsdfzaY/uyPEfJtn2Uvuxngp5R33qZ/1GSHgRr+tbUgSEH9NM/P2nl4ZX6ft/3mrJ9ikY/PFphrTO/hZV040t+/Zb2U9caXTWy+UiVdi8tZydnjd843uqTKltk9hrm1MPlH9bhlw5r6f6lWnVklb7+82tN2DJB0x6bZvnE4+WmL6tztc5asn+Jfjn8i95a85bGbxyv3/r8pgb+DbJctr3PmbcqWqSo1vdfrzVH1+jngz9r5aGVWrB3gR6p+IhWPbMqy+NCyptz2J243bFrMpn0Xc/vtOXEFv144Ef9cvgXDVg2QB9v/lhbBm0p0L8ncoPA68C61eym//z0H205scUyWD0z5b3K69cjv+py0mWrq7z7Y/ffmO5d3uZ1Z/XLzx761Ouj4b8M1+nLpzV3z1x1qtrJ6upVZnaf3a1/zv+jWV1nqU+9Ppb2m781fLNetXvpv8v/qwOxB7Rg7wIVK1JMnat1tkyvXOLGcAxPs6faVmp7R9tx8epFRRyNUFirMI1uOdrSfuvHgKXdS8vNxU2HLh7KsIxDF6zbKpeoLEOGKpaoqGqlqtlUj28xXxUrUkwHzh/IMG1/7H45mZwU6BVo0zKlG8fPnrN7ZBiG1VXeA7HW67F1H+VEhyod5FvMV3N2zVGTsk10JeWKnq337G3ny+r4tcd+t1X6Fdx0hmHo0IVDlpB0N49FKevXJv09eOv9edM/Rs2J3By/2bFl+3KqvHd5pRlpOnrxqKqW+veK9K3vSd9ivvJw9VBqWmqO9o+7q7t61emlXnV6KTk1WU8seELvbHhHox4aleWt+L77+ztVKlFJi3sutnqPjVlrPXShconK+uXQL7pw9UKOrvLeuh22nB9KFi2p/g36q3+D/kpITtDDMx9W6NpQqyE+lUtW1ivNX9ErzV/RwfMHVf+L+vp488ea/cTsTGvIi/2YGSeTk9pUaqM2ldrok/af6N0N7+qN397Qmqg1alupbZbvAVvOYTn9HZlX52VJalq2qZqWbap32ryjubvnqvfi3pq/Z77VPnIkjOF1YMVdi2tqp6kKbRlqFdZu1bFqR6UaqZq0dZJV+4QtE2SSKcM40Zxwd3WXlPGXnz2E1AmRyWTSsJXDdOTikRyNx0y/onbzFTTDMPTpH59m2r97re5yNjlr3p55WvT3Ij1W7THLNklSo4BGqlyisj7a9JESkhMyzH/rLWAyrckpY02SNHHLxAz92lZqqyX7l+jU5VOW9kMXDmUYX/1EzSfkbHJW2LqwDMs1DCPbW/c4OzmrXeV2Wrp/qdWtgc4knNHcPXPVolwLy5AOW3Ss0lGnLp/Sd39/Z2m7knIlw8e7tu6jnHBxclFInRAt/HuhwiPDVbd0XUtQzE5Wx6899rutvvnrG6sx69/9/Z1OJ5y2vC/v5rEoZf3aeJo95VPMR+uPrbdqn7Jtym3Xny43x292bNm+nGpfub2kjNv3+dbPM6y7e63u+n7f99pzdk+G5dy8f27dPldnV9XyrSVDRrZjLS3bd9PV5z9O/KHNxzdb9etes7sMGQpbm/Fq8e0+XbDl/HDrdhR3La4qJasoKfXGrQivpFzJcF/fyiUry8PVw9In2+2043681YWrFzK0pX9ykn4rxazeA7acw3L6OzIvzssXr17M8Breuo2OiCu8Dq5v/azH9aXrXL2zWldorTd+e0NRl6JUz6+eVh1epaUHlurlJi9bfbkspxr5N5IkvfHbG3qq9lMq4lxEnat1tgqNd8rX3VcdqnTQor8XydvNW52qZrwN0K1q+NRQ5RKVNWL1CJ28fFKeZk99v+/7LL8wUdq9tFpXbK1PNn+iy8mX1at2L6vpTiYnff341wqeE6zaU2qrf/3+us/jPp28fFJrotbI0+ypH0N+zLYmT7OnHi7/sD7Y9IFS0lJ0n8d9WnVklY5ePJqhb2jLUK06vEoPznhQLzR+QalpqZq0bZLqlK5j9d9+KpesrHGPjNOoiFGKuhSlrtW7ysPsoaMXj+qH/T/ouUbPaUTzEVnWNK71OMs9KP97/3/l4uSiL3Z8oaTrSfqg7QfZbk9WBjcarEnbJqnPkj7acXqH/Iv769td36pYkWJW/WzdRznVp14ffbb1M62JWmM1LCU7lUtUlrebt6ZtnyYPVw+5u7qryX1NVLFExVzvd1uVLFpSLWa2UP/6/XUm4Ywm/jFRVUpW0eBGgyXd/WMxu/f2oAaD9N7v72nQskFqHNBY64+t1z/n/8nxtub2+LXH9uVUo4BG6l6zuyb+MVHnr5633JYsfXtvvoL3Xpv3tOboGjX5uokGNxysWr61dOHqBe08vVO/HvlVF169EbLazW4nv+J+ejDwQZVxL6N9sfs0aeskdaraKcN3LG72WNXHtHjfYnVb0E2dqnbS0YtHNW3HNNXyrWX1R1Driq31bNCz+mzrZzp44aA6VLkxDGBD9Aa1rtBaQx4Yku025/T8UGtKLbWq0EqN/BupZNGS2n5qu777+zvL8v85/4/afNNGPWv1VC3fWnJxctEP+3/QmcQzeqr2U1muPy/2463Grhur9cfWq1PVTirvXV5nE89qyrYpKutZVi3K3binc1bnB1vOYenvo5dWvKT2ldvL2clZT9XJfNvtfV6e9dcsTdk2Rd1qdFPlkpV1Oemyvtr5lTzNnupYtaPNyyssCLyQk8lJy0KWafSa0Vqwd4FmRs5UBe8K+vDRD/VKs1fuaJn333e/3m79tqZtn6aVh1be+Ohv2FG7BF5J6hPURz/985N61uqZo0H2RZyL6MeQH/XSypc0fuN4ubm4qVuNbhrywBDVm1Yv03l61e6lX4/8Kg9Xj0xPAq0qtNLmgZv19vq3NWnrJCUkJ8ivuJ+alG2i/zT6T462Y+4TczV0xVBN3jZZhmGoXeV2WtF7hQI+CbDq1yigkVb0XqERq0borTVvKdAzUGNbjdW+2H2WoSfpXmvxmqqVqqYJWyZY/nVloFeg2lVup8erP55tPbVL19aG/hs0KmKUxm8crzQjTU3ua6LZ3WZnuNdjThUrUkwRfSI0dMVQfb71cxUrUky96/ZWcJVgdZjTwdLvTvZRTjQKaKTavrW1L3afetftnaN5ijgX0ayuszQqYpSe//l5XU+7rpldZqpiiYp22e+2eP2h17XrzC6N3zhel5Muq03FNprSaYrVHwx381jM7r09uuVonbtyTt/9/Z0W7l2o4KrBWtF7hUp/VDrH25ub49ce22eLb7p9I7/ifpq3Z55+2P+D2lZqqwVPLlD1SdWthh+UKV5GWwdv1dh1Y7V432JN2TZFpYqVUm3f2lZ/hP2n0X80Z/ccfbL5EyUkJ6isZ1m91OSl2/6b2X71+ykmIUZf7PhCvxz6RbV8a2l2t9la9PcirY1aa9V3ZpeZCioTpOl/TtfI1SPlZfZS44DGlvvIZien54eXHnhJy/5ZplWHVynpepLKe5fXuEfGaWTzG/c2D/QMVEidEEUcjdC3u76Vi5OLavjU0MInF6p7re7Z1pAX+/Fmj1d/XFGXojQjcoZir8TKp5iPWpZvqbBWYfJy85KU9fmhX/1+OT6HPVHzCQ19YKjm75mv2btmy5CRZeC193m5ZfmW2npyq+bvna8zCWfk5ealB+57QHOemHPbux0VZibjTr4lAeSzpfuXquuCrlrfb73lxt/3oq7zu964/c9Q+45hczQNvmigkkVLKqJPRH6XkmNro9aq9azWWtRjUZZ3WEHBExkTqQZfNNDsbrPVOyhnf2AByHuM4UWh9NXOr1SpRCXLR0z3glv/ve7B8we1/OBytSrfKn8KKiS2n9quyJhI9Qnqc/vOgA1ufU9KN8aTOpmcMvzHPAD5iyENKFTm75mvXWd26eeDP+vTDp/m+t9IFiaVPqukfvX6qVKJSjoWd0xTt0+Vq7Or/u/B/8vv0gqkPWf3aMepHfp488fyL+6vXnV63X4mwAYf/P6BdpzeodYVWsvFyUUrDq3QikMr9FzD5+74m/MA8gaBF4VKyPchKu5aXAMbDMzwH44cXYcqHTRvzzzFJMTI7GJWs7LN9G6bd61uiYR/fff3dxq7bqyq+1TXvO7zsrylE3Cnmgc21+ojq/X2+reVkJygcl7lFNoyVG88/EZ+lwbgFozhBQAAgENjDC8AAAAcGoEXAAAADo0xvJlIS0vTqVOn5OHhcU99KQoAAKCwMAxDly9fVkBAgJycsr+GS+DNxKlTpxQYyDdsAQAACrrjx4+rbNmy2fYh8GbCw+PGv3A8fvy4PD1t+x/VAAAAyHvx8fEKDAy05LbsEHgzkT6MwdPTk8ALAABQgOVk+ClfWgMAAIBDI/ACAADAoRF4AQAA4NAIvAAAAHBoBF4AAAA4NAIvAAAAHBqBFwAAAA6NwAsAAACHlq+Bd/369ercubMCAgJkMpm0ZMkSq+kmkynTx4cffpjlMkNDQzP0r1GjRh5vCQAAAAqqfA28iYmJqlevniZPnpzp9NOnT1s9ZsyYIZPJpO7du2e73Nq1a1vNt3HjxrwoHwAAAIVAvv5r4eDgYAUHB2c53c/Pz+r50qVL1bp1a1WqVCnb5bq4uGSYFwAAAPemfA28tjhz5ox+/vlnzZo167Z9Dx48qICAALm5ualZs2YaP368ypUrl2X/pKQkJSUlWZ7Hx8fbpWYAAHDvMYWZ8ruEfGOMMfK7hEwVmi+tzZo1Sx4eHnriiSey7dekSROFh4dr5cqVmjp1qo4ePaqHHnpIly9fznKe8ePHy8vLy/IIDAy0d/kAAADIJ4Um8M6YMUO9e/eWm5tbtv2Cg4PVo0cPBQUFqX379lq+fLkuXbqkhQsXZjnPqFGjFBcXZ3kcP37c3uUDAAAgnxSKIQ0bNmzQgQMHtGDBApvn9fb2VrVq1XTo0KEs+5jNZpnN5tyUCAAAgAKqUFzhnT59uho1aqR69erZPG9CQoIOHz4sf3//PKgMAAAABV2+Bt6EhARFRkYqMjJSknT06FFFRkYqOjra0ic+Pl6LFi3SoEGDMl1GmzZtNGnSJMvzESNGaN26dYqKitKmTZvUrVs3OTs7KyQkJE+3BQAAAAVTvg5p2L59u1q3bm15Pnz4cElS3759FR4eLkmaP3++DMPIMrAePnxYsbGxlucnTpxQSEiIzp8/L19fX7Vo0UJbtmyRr69v3m0IAAAACiyTYRgF8/4R+Sg+Pl5eXl6Ki4uTp6dnfpcDAAAKEW5LdnfYktcKxRheAAAA4E4ReAEAAODQCLwAAABwaAReAAAAODQCLwAAABwagRcAAAAOjcALAAAAh0bgBQAAgEMj8AIAAMChEXgBAADg0Ai8AAAAcGgEXgAAADg0Ai8AAAAcGoEXAAAADo3ACwAAAIdG4AUAAIBDI/ACAADAoRF4AQAA4NAIvAAAAHBoBF4AAAA4NAIvAAAAHBqBFwAAAA6NwAsAAACHRuAFAACAQyPwAgAAwKEReAEAAODQCLwAAABwaAReAAAAODQCLwAAABwagRcAAAAOjcALAAAAh5avgXf9+vXq3LmzAgICZDKZtGTJEqvp/fr1k8lksnp06NDhtsudPHmyKlSoIDc3NzVp0kRbt27Noy0AAABAQZevgTcxMVH16tXT5MmTs+zToUMHnT592vKYN29etstcsGCBhg8frjFjxmjnzp2qV6+e2rdvr7Nnz9q7fAAAABQCLvm58uDgYAUHB2fbx2w2y8/PL8fL/OSTTzR48GD1799fkjRt2jT9/PPPmjFjhl577bVc1QsAAIDCp8CP4V27dq1Kly6t6tWr64UXXtD58+ez7JucnKwdO3aobdu2ljYnJye1bdtWmzdvznK+pKQkxcfHWz0AAADgGAp04O3QoYO++eYbRURE6P3339e6desUHBys1NTUTPvHxsYqNTVVZcqUsWovU6aMYmJislzP+PHj5eXlZXkEBgbadTsAAACQf/J1SMPtPPXUU5af69atq6CgIFWuXFlr165VmzZt7LaeUaNGafjw4Zbn8fHxhF4AAAAHUaCv8N6qUqVK8vHx0aFDhzKd7uPjI2dnZ505c8aq/cyZM9mOAzabzfL09LR6AAAAwDEUqsB74sQJnT9/Xv7+/plOd3V1VaNGjRQREWFpS0tLU0REhJo1a3a3ygQAAEABkq+BNyEhQZGRkYqMjJQkHT16VJGRkYqOjlZCQoJGjhypLVu2KCoqShEREerSpYuqVKmi9u3bW5bRpk0bTZo0yfJ8+PDh+uqrrzRr1izt27dPL7zwghITEy13bQAAAMC9JV/H8G7fvl2tW7e2PE8fR9u3b19NnTpVu3bt0qxZs3Tp0iUFBASoXbt2evvtt2U2my3zHD58WLGxsZbnvXr10rlz5zR69GjFxMSofv36WrlyZYYvsgEAAODeYDIMw8jvIgqa+Ph4eXl5KS4ujvG8AADAJqYwU36XkG+MMXcvVtqS1wrVGF4AAADAVgReAAAAODQCLwAAABwagRcAAAAOjcALAAAAh0bgBQAAgEMj8AIAAMChEXgBAADg0Ai8AAAAcGgEXgAAADg0Ai8AAAAcGoEXAAAADo3ACwAAAIdG4AUAAIBDI/ACAADAoRF4AQAA4NAIvAAAAHBoBF4AAAA4NAIvAAAAHBqBFwAAAA6NwAsAAACHRuAFAACAQyPwAgAAwKEReAEAAODQCLwAAABwaAReAAAAODQCLwAAABwagRcAAAAOjcALAAAAh0bgBQAAgEMj8AIAAMCh5WvgXb9+vTp37qyAgACZTCYtWbLEMi0lJUWvvvqq6tatK3d3dwUEBKhPnz46depUtssMDQ2VyWSyetSoUSOPtwQAAAAFVb4G3sTERNWrV0+TJ0/OMO3KlSvauXOn3nrrLe3cuVOLFy/WgQMH9Pjjj992ubVr19bp06ctj40bN+ZF+QAAACgEXPJz5cHBwQoODs50mpeXl1avXm3VNmnSJD3wwAOKjo5WuXLlslyui4uL/Pz8clxHUlKSkpKSLM/j4+NzPC8AAAAKtkI1hjcuLk4mk0ne3t7Z9jt48KACAgJUqVIl9e7dW9HR0dn2Hz9+vLy8vCyPwMBAO1YNAACA/FRoAu+1a9f06quvKiQkRJ6enln2a9KkicLDw7Vy5UpNnTpVR48e1UMPPaTLly9nOc+oUaMUFxdneRw/fjwvNgEAAAD5IF+HNORUSkqKevbsKcMwNHXq1Gz73jxEIigoSE2aNFH58uW1cOFCDRw4MNN5zGazzGazXWsGAABAwVDgA2962D127Jh+++23bK/uZsbb21vVqlXToUOH8qhCAAAAFGQFekhDetg9ePCgfv31V5UqVcrmZSQkJOjw4cPy9/fPgwoBAABQ0OVr4E1ISFBkZKQiIyMlSUePHlVkZKSio6OVkpKiJ598Utu3b9ecOXOUmpqqmJgYxcTEKDk52bKMNm3aaNKkSZbnI0aM0Lp16xQVFaVNmzapW7ducnZ2VkhIyN3ePAAAABQA+TqkYfv27WrdurXl+fDhwyVJffv2VWhoqJYtWyZJql+/vtV8a9asUatWrSRJhw8fVmxsrGXaiRMnFBISovPnz8vX11ctWrTQli1b5Ovrm7cbAwAAgAIpXwNvq1atZBhGltOzm5YuKirK6vn8+fNzWxYAAAAcSIEewwsAAADkFoEXAAAADo3ACwAAAIdG4AUAAIBDI/ACAADAod1R4L2edl2/HvlVX2z/QpeTLkuSTl0+pYTkBLsWBwAAAOSWzbclO3bpmDrM6aDouGglXU/So5UflYfZQ+9vfF9JqUma9ti0vKgTAAAAuCM2X+EdtnKYGgc01sVXL6pokaKW9m41uyniaIRdiwMAAAByy+YrvBuiN2jTgE1ydXa1aq/gXUEn40/arTAAAADAHmy+wptmpCnVSM3QfiL+hDzMHnYpCgAAALAXmwNvu8rtNHHLRMtzk0xKSE7QmLVj1LFqR3vWBgAAAOSazUMaPnr0I3WY00G1JtfStevX9PTip3Xw/EH5FPPRvO7z8qJGAAAKLVOYKb9LyDfGGCO/SwAk3UHgDfQK1F/P/6UFexborzN/KSE5QQMbDFTvur2tvsQGAAAAFAQ2Bd6U1BTVmFxDP4X8pN5BvdVbvfOqLgAAAMAubBrDW8S5iK5dv5ZXtQAAAAB2Z/OX1l68/0W9//v7up52PS/qAQAAAOzK5jG8205tU8SRCK06vEp1y9SVexF3q+mLey22W3EAAABAbtkceL3dvNW9Vve8qAUAAACwO5sD78wuM/OiDgAAACBP2DyGFwAAAChMbL7CW/HTijIp65toHxl2JFcFAQAAAPZkc+B9ucnLVs9T0lL0Z8yfWnlopUY2H2mvugAAAAC7sDnwDms6LNP2yVsna/vp7bkuCAAAALAnu43hDa4arO///t5eiwMAAADswm6B97u/v1PJoiXttTgAAADALmwe0tDgiwZWX1ozZCgmIUbnEs9pSqcpdi0OAAAAyC2bA2+X6l2sAq+TyUm+7r5qVaGVavjUsGtxAAAAQG7ZHHhDW4XmQRkAAABA3rB5DK/zWGedTTybof38lfNyHutsl6IAAAAAe7E58BqGkWl7UmqSXJ1dc10QAAAAYE85Dryf/fGZPvvjM5lMJn2982vL88/++EwTNk/Qi8tftHkM7/r169W5c2cFBATIZDJpyZIlVtMNw9Do0aPl7++vokWLqm3btjp48OBtlzt58mRVqFBBbm5uatKkibZu3WpTXQAAAHAcOR7DO2HLBEk3Qui07dPk7PTv8AVXZ1dV8K6gaZ2m2bTyxMRE1atXTwMGDNATTzyRYfoHH3ygzz77TLNmzVLFihX11ltvqX379vr777/l5uaW6TIXLFig4cOHa9q0aWrSpIkmTpyo9u3b68CBAypdurRN9QEAAKDwMxlZjVHIQutZrbW452KVKFrCvoWYTPrhhx/UtWtXSTeCdUBAgF555RWNGDFCkhQXF6cyZcooPDxcTz31VKbLadKkie6//35NmjRJkpSWlqbAwEANHTpUr732Wo5qiY+Pl5eXl+Li4uTp6Zn7jQMA3LNMYabbd3JQxhibIobDYJ/fHbbkNZvH8K7pu8buYTczR48eVUxMjNq2bWtp8/LyUpMmTbR58+ZM50lOTtaOHTus5nFyclLbtm2znEeSkpKSFB8fb/UAAACAY7D5tmSSdCL+hJYdWKbouGglpyZbTfuk/Sd2KSwmJkaSVKZMGav2MmXKWKbdKjY2VqmpqZnOs3///izXNX78eIWFheWyYgAAABRENgfeiCMRenz+46pUopL2x+5XndJ1FHUpSoZhqKF/w7yoMc+NGjVKw4cPtzyPj49XYGBgPlYEAAAAe7F5SMOoiFEa0WyEdr+wW24ubvq+5/c6/r/jalmhpXrU6mG3wvz8/CRJZ86csWo/c+aMZdqtfHx85OzsbNM8kmQ2m+Xp6Wn1AAAAgGOwOfDui92nPvX6SJJcnFx0NeWqirsW19hWY/X+7+/brbCKFSvKz89PERERlrb4+Hj98ccfatasWabzuLq6qlGjRlbzpKWlKSIiIst5AAAA4NhsDrzuRdwt43b9i/vr8MXDlmmxV2JtWlZCQoIiIyMVGRkp6cYX1SIjIxUdHS2TyaSXX35Z48aN07Jly7R792716dNHAQEBljs5SFKbNm0sd2SQpOHDh+urr77SrFmztG/fPr3wwgtKTExU//79bd1UAAAAOACbx/A2LdtUG6M3qqZvTXWs2lGvrHpFu8/s1uL9i9W0bFOblrV9+3a1bt3a8jx9HG3fvn0VHh6u//u//1NiYqKee+45Xbp0SS1atNDKlSut7sF7+PBhxcb+G7R79eqlc+fOafTo0YqJiVH9+vW1cuXKDF9kAwAAwL3B5vvwHrl4RAnJCQoqE6TE5ES9suoVbTq+SVVLVdUn7T5Ree/yeVXrXcN9eAEA9sI9We897PO7w5a8ZtMV3tS0VJ2IP6GgMkGSJHdXd017zLb/rgYAAADcTTaN4XV2cla7b9vp4tWLeVUPAAAAYFc2f2mtTuk6OnLxSF7UAgAAANidzYF33CPjNGL1CP30z086ffm04pPirR4AAABAQWLzXRo6zukoSXp83uMymf4dlG0Yhkwmk1JHp9qvOgAAACCXbA68a/quyYs6AAAAgDxhc+BtWaFlXtQBAAAA5Ambx/BK0oZjG/TM4mfUfHpznYw/KUn69q9vtTF6o12LAwAAAHLL5sD7/d/fq/3s9irqUlQ7T+9UUmqSJCkuKU7vbnjX7gUCAAAAuWH7XRo2jNO0x6bpq8e/UhHnIpb2BwMf1M7TO+1aHAAAAJBbNgfeA7EH9HD5hzO0e7l56dK1S/aoCQAAALAbmwOvX3E/HbpwKEP7xuiNqlSikl2KAgAAAOzF5sA7uOFgDVs5TH+c+EMmmXTq8inN2TVHI1aN0AuNX8iLGgEAAIA7ZvNtyV5r8ZrSjDS1+aaNrqRc0cMzH5bZxawRzUZoaJOheVEjAAAAcMdsDrwmk0lvPPyGRj44UocuHFJCcoJq+dZScdfieVEfAAAAkCs2B950rs6u8nD1kIerB2EXAAAABZbNgfd62nWFrQ3TZ1s/U0JygiSpuGtxDX1gqMa0HGN1qzIAAAAgv9kceIcuH6rF+xfrg7YfqFlgM0nS5uObFbouVOevnNfUx6bavUgAAADgTtkceOfumav53ecruGqwpS2oTJACvQIV8n0IgRcAAAAFis23JTM7m1XBu0KG9oreFeXq7GqPmgAAAAC7sTnwDnlgiN5e/7aSridZ2pKuJ+mdDe9oyP1D7FocAAAAkFs2D2n4M+ZPRRyJUNkJZVWvTD1J0l9n/lJyarLaVGyjJxY8Yem7uNdi+1UKAAAA3AGbA6+3m7e61+pu1RboFWi3ggAAAAB7sjnwzuwyMy/qAAAAAPKEzWN4AQAAgMLE5iu856+c1+g1o7Umao3OJp5VmpFmNf3CqxfsVhwAAACQWzYH3md/eFaHLhzSwAYDVaZ4GZlkyou6AAAAALuwOfBuiN6gjf03qp5fvbyoBwAAALArm8fw1vCpoavXr+ZFLQAAAIDd2Rx4p3Scojd+e0Protbp/JXzik+Kt3oAAAAABckd3Yc3Pilej3zziFW7YRgymUxKHZ1qt+IAAACA3LI58PZe3FtFnIpo7hNz78qX1ipUqKBjx45laP/vf/+ryZMnZ2gPDw9X//79rdrMZrOuXbuWZzUCAACg4LI58O45u0d//udPVfepnhf1ZLBt2zalpv571XjPnj169NFH1aNHjyzn8fT01IEDByzPTSbuJAEAAHCvsjnwNg5orOPxx+9a4PX19bV6/t5776ly5cpq2bJllvOYTCb5+fnldWkAAAAoBGwOvEMfGKphK4dpZPORqlu6roo4F7GaHlQmyG7F3So5OVmzZ8/W8OHDs71qm5CQoPLlyystLU0NGzbUu+++q9q1a2fZPykpSUlJSZbn8fF8+Q4AAMBR2Bx4e33XS5I0YOkAS5vJZLorX1pbsmSJLl26pH79+mXZp3r16poxY4aCgoIUFxenjz76SM2bN9fevXtVtmzZTOcZP368wsLC8qhqAAAA5CeTYRiGLTMcu5TxC2Q3K+9dPlcFZad9+/ZydXXVjz/+mON5UlJSVLNmTYWEhOjtt9/OtE9mV3gDAwMVFxcnT0/PXNcNALh3mcLu3e+RGGNsihgOg31+d8THx8vLyytHec3mK7x5GWizc+zYMf36669avHixTfMVKVJEDRo00KFDh7LsYzabZTabc1siAAAACiCb//GEJH3717d6cMaDCvg4wHLFd+KWiVq6f6ldi7vZzJkzVbp0aXXq1Mmm+VJTU7V79275+/vnUWUAAAAoyGwOvFO3TdXwVcPVsUpHXbp2SanGjTG73m7emvjHRHvXJ0lKS0vTzJkz1bdvX7m4WF+U7tOnj0aNGmV5PnbsWK1atUpHjhzRzp079cwzz+jYsWMaNGhQntQGAACAgs3mwPv51s/1Veev9MbDb8jZydnS3jigsXaf2W3X4tL9+uuvio6O1oABAzJMi46O1unTpy3PL168qMGDB6tmzZrq2LGj4uPjtWnTJtWqVStPagMAAEDBZvMY3qOXjqqBX4MM7WZnsxJTEu1S1K3atWunrL5bt3btWqvnEyZM0IQJE/KkDgAAABQ+Nl/hrehdUZExkRnaVx5aqZo+Ne1REwAAAGA3Ob7CO3bdWI1oPkLDmw3Xi8tf1LXr12QYhrae3Kp5u+dp/Mbx+vrxr/OyVgAAAMBmOQ68YevC9Hzj5zWo4SAVdSmqN9e8qSspV/T0908rwCNAn3b4VE/VeSovawWAQu9evT/nvXo/VgAFQ44D781jaHsH9VbvoN66knJFCckJKu1eOk+KAwAAAHLLpi+tmWR9ZaJYkWIqVqSYXQsCAAAA7MmmwFttUrUMofdWF169kKuCAAAAAHuyKfCGtQqTl9krr2oBAAAA7M6mwPtUnacYrwsAAIBCJcf34TWZ7s1vFgMAAKBwy3Hgzeo/nQEAAAAFWY6HNKSNScvLOgAAAIA8YfO/FgYAAAAKEwIvAAAAHBqBFwAAAA6NwAsAAACHRuAFAACAQyPwAgAAwKEReAEAAODQCLwAAABwaAReAAAAOLQc/6c1APZlCjPldwn5xhjDvyoHANw9XOEFAACAQyPwAgAAwKEReAEAAODQCLwAAABwaAReAAAAODQCLwAAABwagRcAAAAOjcALAAAAh0bgBQAAgEMj8AIAAMChFejAGxoaKpPJZPWoUaNGtvMsWrRINWrUkJubm+rWravly5ffpWoBAABQEBXowCtJtWvX1unTpy2PjRs3Ztl306ZNCgkJ0cCBA/Xnn3+qa9eu6tq1q/bs2XMXKwYAAEBBUuADr4uLi/z8/CwPHx+fLPt++umn6tChg0aOHKmaNWvq7bffVsOGDTVp0qS7WDEAAAAKkgIfeA8ePKiAgABVqlRJvXv3VnR0dJZ9N2/erLZt21q1tW/fXps3b852HUlJSYqPj7d6AAAAwDEU6MDbpEkThYeHa+XKlZo6daqOHj2qhx56SJcvX860f0xMjMqUKWPVVqZMGcXExGS7nvHjx8vLy8vyCAwMtNs2AAAAIH8V6MAbHBysHj16KCgoSO3bt9fy5ct16dIlLVy40K7rGTVqlOLi4iyP48eP23X5AAAAyD8u+V2ALby9vVWtWjUdOnQo0+l+fn46c+aMVduZM2fk5+eX7XLNZrPMZrPd6gQAAEDBUaCv8N4qISFBhw8flr+/f6bTmzVrpoiICKu21atXq1mzZnejPAAAABRABTrwjhgxQuvWrVNUVJQ2bdqkbt26ydnZWSEhIZKkPn36aNSoUZb+w4YN08qVK/Xxxx9r//79Cg0N1fbt2zVkyJD82gQAAADkswI9pOHEiRMKCQnR+fPn5evrqxYtWmjLli3y9fWVJEVHR8vJ6d/M3rx5c82dO1dvvvmmXn/9dVWtWlVLlixRnTp18msTAAAAkM8KdOCdP39+ttPXrl2boa1Hjx7q0aNHHlUEAACAwqZAD2kAAAAAcovACwAAAIdG4AUAAIBDI/ACAADAoRF4AQAA4NAIvAAAAHBoBF4AAAA4NAIvAAAAHBqBFwAAAA6NwAsAAACHRuAFAACAQyPwAgAAwKEReAEAAODQCLwAAABwaAReAAAAODQCLwAAABwagRcAAAAOzSW/C8ANpjBTfpeQb4wxRn6XAAAAHBhXeAEAAODQCLwAAABwaAReAAAAODQCLwAAABwagRcAAAAOjcALAAAAh0bgBQAAgEMj8AIAAMChEXgBAADg0Ai8AAAAcGgEXgAAADg0Ai8AAAAcGoEXAAAADq1AB97x48fr/vvvl4eHh0qXLq2uXbvqwIED2c4THh4uk8lk9XBzc7tLFQMAAKCgKdCBd926dXrxxRe1ZcsWrV69WikpKWrXrp0SExOznc/T01OnT5+2PI4dO3aXKgYAAEBB45LfBWRn5cqVVs/Dw8NVunRp7dixQw8//HCW85lMJvn5+eV1eQAAACgECvQV3lvFxcVJkkqWLJltv4SEBJUvX16BgYHq0qWL9u7dm23/pKQkxcfHWz0AAADgGApN4E1LS9PLL7+sBx98UHXq1MmyX/Xq1TVjxgwtXbpUs2fPVlpampo3b64TJ05kOc/48ePl5eVleQQGBubFJgAAACAfFJrA++KLL2rPnj2aP39+tv2aNWumPn36qH79+mrZsqUWL14sX19fffHFF1nOM2rUKMXFxVkex48ft3f5AAAAyCcFegxvuiFDhuinn37S+vXrVbZsWZvmLVKkiBo0aKBDhw5l2cdsNstsNue2TAAAABRABfoKr2EYGjJkiH744Qf99ttvqlixos3LSE1N1e7du+Xv758HFQIAAKCgK9BXeF988UXNnTtXS5culYeHh2JiYiRJXl5eKlq0qCSpT58+uu+++zR+/HhJ0tixY9W0aVNVqVJFly5d0ocffqhjx45p0KBB+bYdAAAAyD8FOvBOnTpVktSqVSur9pkzZ6pfv36SpOjoaDk5/Xuh+uLFixo8eLBiYmJUokQJNWrUSJs2bVKtWrXuVtkAAAAoQAp04DUM47Z91q5da/V8woQJmjBhQh5VBAAAgMKmQI/hBQAAAHKLwAsAAACHRuAFAACAQyPwAgAAwKEReAEAAODQCLwAAABwaAReAAAAODQCLwAAABwagRcAAAAOjcALAAAAh0bgBQAAgEMj8AIAAMChEXgBAADg0Ai8AAAAcGgEXgAAADg0Ai8AAAAcGoEXAAAADo3ACwAAAIdG4AUAAIBDI/ACAADAoRF4AQAA4NAIvAAAAHBoBF4AAAA4NAIvAAAAHBqBFwAAAA6NwAsAAACHRuAFAACAQyPwAgAAwKEReAEAAODQCLwAAABwaAReAAAAOLRCEXgnT56sChUqyM3NTU2aNNHWrVuz7b9o0SLVqFFDbm5uqlu3rpYvX36XKgUAAEBBU+AD74IFCzR8+HCNGTNGO3fuVL169dS+fXudPXs20/6bNm1SSEiIBg4cqD///FNdu3ZV165dtWfPnrtcOQAAAAqCAh94P/nkEw0ePFj9+/dXrVq1NG3aNBUrVkwzZszItP+nn36qDh06aOTIkapZs6befvttNWzYUJMmTbrLlQMAAKAgcMnvArKTnJysHTt2aNSoUZY2JycntW3bVps3b850ns2bN2v48OFWbe3bt9eSJUuyXE9SUpKSkpIsz+Pi4iRJ8fHxuajeRtfu3qoKmrv6Ohck7PN70z2639nn96Z7dr+zz+/qugzDuG3fAh14Y2NjlZqaqjJlyli1lylTRvv37890npiYmEz7x8TEZLme8ePHKywsLEN7YGDgHVQNW3m955XfJeAuY5/fe9jn9yb2+70nP/b55cuX5eWV/XoLdOC9W0aNGmV1VTgtLU0XLlxQqVKlZDKZ8rGyuyM+Pl6BgYE6fvy4PD0987sc3AXs83sP+/zewz6/N91L+90wDF2+fFkBAQG37VugA6+Pj4+cnZ115swZq/YzZ87Iz88v03n8/Pxs6i9JZrNZZrPZqs3b2/vOii7EPD09Hf7NAWvs83sP+/zewz6/N90r+/12V3bTFegvrbm6uqpRo0aKiIiwtKWlpSkiIkLNmjXLdJ5mzZpZ9Zek1atXZ9kfAAAAjq1AX+GVpOHDh6tv375q3LixHnjgAU2cOFGJiYnq37+/JKlPnz667777NH78eEnSsGHD1LJlS3388cfq1KmT5s+fr+3bt+vLL7/Mz80AAABAPinwgbdXr146d+6cRo8erZiYGNWvX18rV660fDEtOjpaTk7/Xqhu3ry55s6dqzfffFOvv/66qlatqiVLlqhOnTr5tQkFntls1pgxYzIM64DjYp/fe9jn9x72+b2J/Z45k5GTezkAAAAAhVSBHsMLAAAA5BaBFwAAAA6NwAsAAACHRuAFCjGTyZTtv80GcotjzLH169dPXbt2tTxv1aqVXn755XyrB8grBN4C7taTERxPv379ZDKZMjwOHTqU36Uhn6UfG88//3yGaS+++KJMJpP69et39wtDvti8ebOcnZ3VqVMnuy3z008/VXh4uN2WB/viHGA/BF7YRUpKSn6XUKh16NBBp0+ftnpUrFgxT9aVl/uK48D+AgMDNX/+fF29etXSdu3aNc2dO1flypXL1bLZX4XL9OnTNXToUK1fv16nTp3K1bJSU1OVlpYmLy+ve/I/ixYmeXkOuJcQeAuRlStXqkWLFvL29lapUqX02GOP6fDhw5bpUVFRMplMWrx4sVq3bq1ixYqpXr162rx5s6VPaGio6tevb7XciRMnqkKFCpbn27Zt06OPPiofHx95eXmpZcuW2rlzp9U8JpNJU6dO1eOPPy53d3eNGzdOVapU0UcffWTVLzIykquVOWA2m+Xn52f1cHZ21tKlS9WwYUO5ubmpUqVKCgsL0/Xr163mPX36tIKDg1W0aFFVqlRJ3333nWVa+jGxYMECtWzZUm5ubpozZw7HQSHSsGFDBQYGavHixZa2xYsXq1y5cmrQoIGlLafnh1uPBUmaMWOGateuLbPZLH9/fw0ZMsSqhtjYWHXr1k3FihVT1apVtWzZsjzeatwqISFBCxYs0AsvvKBOnTpZXZVdu3atTCaTfv75ZwUFBcnNzU1NmzbVnj17LH3Cw8Pl7e2tZcuWqVatWjKbzYqOjuZTxELAXueARx55JMN7+9y5c3J1dc3wH2odEYG3EElMTNTw4cO1fft2RUREyMnJSd26dVNaWppVvzfeeEMjRoxQZGSkqlWrppCQkAwhKTuXL19W3759tXHjRm3ZskVVq1ZVx44ddfnyZat+oaGh6tatm3bv3q2BAwdqwIABmjlzplWfmTNn6uGHH1aVKlXufMPvURs2bFCfPn00bNgw/f333/riiy8UHh6ud955x6rfW2+9pe7du+uvv/5S79699dRTT2nfvn1WfV577TUNGzZM+/btU/v27XO0fo6DguPW13TGjBmW/zaZLqfnh1uPhalTp+rFF1/Uc889p927d2vZsmUZ9lNYWJh69uypXbt2qWPHjurdu7cuXLiQdxuMDBYuXKgaNWqoevXqeuaZZzRjxgzdehv9kSNH6uOPP9a2bdvk6+urzp07W13Fv3Llit5//319/fXX2rt3r0qXLn23NwN3yB7ngEGDBmnu3LlKSkqyzDN79mzdd999euSRR+7OhuQnAwVa3759jS5dumQ67dy5c4YkY/fu3YZhGMbRo0cNScbXX39t6bN3715DkrFv3z7DMAxjzJgxRr169ayWM2HCBKN8+fJZ1pCammp4eHgYP/74o6VNkvHyyy9b9Tt58qTh7Oxs/PHHH4ZhGEZycrLh4+NjhIeH53Rz70l9+/Y1nJ2dDXd3d8vjySefNNq0aWO8++67Vn2//fZbw9/f3/JckvH8889b9WnSpInxwgsvGIbx7zExceJEqz4cB4VD+vv/7NmzhtlsNqKiooyoqCjDzc3NOHfunNGlSxejb9++mc6b1fnh1mMhICDAeOONN7KsQZLx5ptvWp4nJCQYkowVK1bkfgORY82bN7fsu5SUFMPHx8dYs2aNYRiGsWbNGkOSMX/+fEv/8+fPG0WLFjUWLFhgGIZhzJw505BkREZGWi331t8xLVu2NIYNG5an24Kcs+c54OrVq0aJEiUsx4RhGEZQUJARGhp6NzYl33GFtxA5ePCgQkJCVKlSJXl6elo+fo6OjrbqFxQUZPnZ399fknT27Nkcr+fMmTMaPHiwqlatKi8vL3l6eiohISHDeho3bmz1PCAgQJ06ddKMGTMkST/++KOSkpLUo0ePHK/7XtW6dWtFRkZaHp999pn++usvjR07VsWLF7c8Bg8erNOnT+vKlSuWeZs1a2a1rGbNmmW4wnvrvsoJjoOCw9fX1/Ix9syZM9WpUyf5+PhY9cnp+eHm/XX27FmdOnVKbdq0yXb9N59T3N3d5enpadM5Bblz4MABbd26VSEhIZIkFxcX9erVS9OnT7fqd/O5oGTJkqpevbrVucDV1dVqX6LwsMc5wM3NTc8++6zl3Lxz507t2bPnnvnSm0t+F4Cc69y5s8qXL6+vvvpKAQEBSktLU506dZScnGzVr0iRIpafTSaTJFk+0nBycsrwMditX1zp27evzp8/r08//VTly5eX2WxWs2bNMqzH3d09Q42DBg3Ss88+qwkTJmjmzJnq1auXihUrducbfY9wd3fP8DFyQkKCwsLC9MQTT2To7+bmZvPyb8ZxUPgMGDDAMv5u8uTJGabn9Pxw8/4qWrRojtZ98zlFunFeuXWoBPLO9OnTdf36dQUEBFjaDMOQ2WzWpEmTcrycokWLWn4noPCxxzlg0KBBql+/vk6cOKGZM2fqkUceUfny5e/aNuQnAm8hcf78eR04cEBfffWVHnroIUnSxo0bbV6Or6+vYmJiZBiG5cQXGRlp1ef333/XlClT1LFjR0nS8ePHFRsbm6Pld+zYUe7u7po6dapWrlyp9evX21wjbmjYsKEOHDhw23GvW7ZsUZ8+faye3/xFhsxwHBQ+HTp0UHJyskwmU4Zx2Hd6fvDw8FCFChUUERGh1q1b50ndyJ3r16/rm2++0ccff6x27dpZTevatavmzZunGjVqSLrx3k//1v7Fixf1zz//qGbNmne9ZuQNe5wD6tatq8aNG+urr77S3LlzbfqDqbAj8BYSJUqUUKlSpfTll1/K399f0dHReu2112xeTqtWrXTu3Dl98MEHevLJJ7Vy5UqtWLFCnp6elj5Vq1bVt99+q8aNGys+Pl4jR47M8ZUgZ2dn9evXT6NGjVLVqlUzfNyOnBs9erQee+wxlStXTk8++aScnJz0119/ac+ePRo3bpyl36JFi9S4cWO1aNFCc+bM0datWzN81HkrjoPCx9nZ2fLxtLOzs9W03JwfQkND9fzzz6t06dIKDg7W5cuX9fvvv2vo0KF23wbY7qefftLFixc1cOBAeXl5WU3r3r27pk+frg8//FCSNHbsWJUqVUplypTRG2+8IR8fH+7A4EDsdQ4YNGiQhgwZInd3d3Xr1i3P6y4oGMNbwKWlpcnFxUVOTk6aP3++duzYoTp16uh///uf5SRni5o1a2rKlCmaPHmy6tWrp61bt2rEiBFWfaZPn66LFy+qYcOGevbZZ/XSSy/Z9G3egQMHKjk5OcM3SGGb9u3b66efftKqVat0//33q2nTppowYUKGj5/CwsI0f/58BQUF6ZtvvtG8efNUq1atbJfNcVA4eXp6Wv1Rki4354e+fftq4sSJmjJlimrXrq3HHntMBw8etHfpuEPTp09X27ZtM4Rd6Ubg3b59u3bt2iVJeu+99zRs2DA1atRIMTEx+vHHH+Xq6nq3S0Yessc5ICQkRC4uLgoJCbF5eFxhZjJuHciHAqVDhw6qUqVKofrYYcOGDWrTpo2OHz+uMmXK5Hc5yCccB8DdsXbtWrVu3VoXL17kn0jgtqKiolS5cmVt27ZNDRs2zO9y7hqu8BZQFy9e1E8//aS1a9eqbdu2+V1OjiQlJenEiRMKDQ1Vjx49CDn3KI4DACh4UlJSFBMTozfffFNNmza9p8KuROAtsAYMGKDnn39er7zyirp06ZLf5eTIvHnzVL58eV26dEkffPBBfpeDfMJxAAAFz++//y5/f39t27ZN06ZNy+9y7jqGNAAAAMChcYUXAAAADo3ACwAAAIdG4AUAAIBDI/ACAADAoRF4AQAA4NAIvACAbJnCTFqyf0l+lwEAd4zACwAFVL8l/WQKM+n5n57PMO3Fn1+UKcykfkv62W19oWtDVX9afbstDwAKCgIvABRggZ6Bmr9nvq6mXLW0Xbt+TXP3zFU5r3L5WBkAFB4EXgAowBr6N1SgV6AW71tsaVu8b7HKeZVTA78Glrak60l6acVLKv1habmNc1OLGS207eQ2y/S1UWtlCjMp4kiEGn/ZWMXeKabm05vrQOwBSVJ4ZLjC1oXprzN/yRRmkinMpPDIcMv8sVdi1W1BNxV7p5iqfl5Vyw4sy/uNBwA7IfACQAE3oP4AzYycaXk+488Z6l+/v1Wf/1v9f/p+3/ea1XWWdv5np6qUrKL2s9vrwtULVv3e+O0NfdzuY21/brtcnFw0YNkASVKv2r30SrNXVNu3tk6/clqnXzmtXrV7WeYLWxemnrV6atcLu9SxSkf1Xtw7w7IBoKAi8AJAAfdM0DPaGL1Rxy4d07FLx/T78d/1TNAzlumJyYmaun2qPnz0QwVXDVYt31r6qvNXKlqkqKbvnG61rHceeUctK7RULd9aeq3Fa9p0fJOuXb+mokWKqrhrcbk4ucivuJ/8ivupaJGilvn61eunkLohqlKyit5t864SkhO09eTWu/YaAEBuuOR3AQCA7Pm6+6pTtU4KjwyXIUOdqnaSTzEfy/TDFw8rJS1FDwY+aGkr4lxED9z3gPbF7rNaVlCZIMvP/sX9JUlnE8/edjzwzfO5u7rL0+yps4lnc7VdAHC3EHgBoBAYUH+AhqwYIkma3HHyHS+niHMRy88mk0mSlGak2TSfJJlkytF8AFAQMKQBAAqBDlU6KDk1WSmpKWpfub3VtMolKsvV2VW/H//d0paSmqJtJ7eplm+tHK/D1dlVqUaq3WoGgIKCK7wAUAg4Ozlr34v7LD/fzN3VXS80fkEjV49UyaIlVc6rnD74/QNdSbmigQ0G5ngdFbwr6OjFo4qMiVRZz7LycPWQ2cVs1+0AgPxA4AWAQsLT7JnltPfavqc0I03P/vCsLiddVuOAxvrlmV9UomiJHC+/e83uWrxvsVrPaq1L1y5pZpeZ6le/nx0qB4D8ZTIMw8jvIgAAAIC8whheAAAAODQCLwAAABwagRcAAAAOjcALAAAAh0bgBQAAgEMj8AIAAMChEXgBAADg0Ai8AAAAcGgEXgAAADg0Ai8AAAAcGoEXAAAADu3/AWP4ONXu99qAAAAAAElFTkSuQmCC\n"
          },
          "metadata": {}
        }
      ]
    }
  ]
}