ramanActiveMode.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "outputs": [],
   "source": [],
   "metadata": {
    "collapsed": false,
    "ExecuteTime": {
     "end_time": "2024-04-05T02:48:32.294981Z",
     "start_time": "2024-04-05T02:48:32.292618Z"
    }
   },
   "id": "initial_id",
   "execution_count": 34
  },
  {
   "cell_type": "code",
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np"
   ],
   "metadata": {
    "collapsed": false,
    "ExecuteTime": {
     "end_time": "2024-04-05T02:48:32.304938Z",
     "start_time": "2024-04-05T02:48:32.302491Z"
    }
   },
   "id": "a5d455d7fc2115d",
   "execution_count": 35
  },
  {
   "cell_type": "code",
   "outputs": [],
   "source": [
    "# Original format data provided by the user\n",
    "data = \"\"\"\n",
    "     freq (  1 -  1) =        -56.9  [cm-1]   --> A               I+R\n",
    "     freq (  2 -  2) =        -27.2  [cm-1]   --> A               I+R\n",
    "     freq (  3 -  3) =        -25.9  [cm-1]   --> E               I+R\n",
    "     freq (  4 -  4) =        -25.9  [cm-1]   --> E*              I+R\n",
    "     freq (  5 -  5) =         -5.6  [cm-1]   --> A               I+R\n",
    "     freq (  6 -  6) =         27.8  [cm-1]   --> E               I+R\n",
    "     freq (  7 -  7) =         27.8  [cm-1]   --> E*              I+R\n",
    "     freq (  8 -  8) =         38.9  [cm-1]   --> E               I+R\n",
    "     freq (  9 -  9) =         38.9  [cm-1]   --> E*              I+R\n",
    "     freq ( 10 - 10) =         47.6  [cm-1]   --> A               I+R\n",
    "     freq ( 11 - 11) =         47.9  [cm-1]   --> E               I+R\n",
    "     freq ( 12 - 12) =         47.9  [cm-1]   --> E*              I+R\n",
    "     freq ( 13 - 13) =         50.4  [cm-1]   --> A               I+R\n",
    "     freq ( 14 - 14) =         50.4  [cm-1]   --> E               I+R\n",
    "     freq ( 15 - 15) =         50.4  [cm-1]   --> E*              I+R\n",
    "     freq ( 16 - 16) =         56.6  [cm-1]   --> A               I+R\n",
    "     freq ( 17 - 17) =         57.3  [cm-1]   --> E               I+R\n",
    "     freq ( 18 - 18) =         57.3  [cm-1]   --> E*              I+R\n",
    "     freq ( 19 - 19) =         60.0  [cm-1]   --> A               I+R\n",
    "     freq ( 20 - 20) =         62.9  [cm-1]   --> E               I+R\n",
    "     freq ( 21 - 21) =         62.9  [cm-1]   --> E*              I+R\n",
    "     freq ( 22 - 22) =         67.6  [cm-1]   --> E               I+R\n",
    "     freq ( 23 - 23) =         67.6  [cm-1]   --> E*              I+R\n",
    "     freq ( 24 - 24) =         68.8  [cm-1]   --> A               I+R\n",
    "     freq ( 25 - 25) =         76.4  [cm-1]   --> A               I+R\n",
    "     freq ( 26 - 26) =         76.6  [cm-1]   --> E               I+R\n",
    "     freq ( 27 - 27) =         76.6  [cm-1]   --> E*              I+R\n",
    "     freq ( 28 - 28) =         83.3  [cm-1]   --> E               I+R\n",
    "     freq ( 29 - 29) =         83.3  [cm-1]   --> E*              I+R\n",
    "     freq ( 30 - 30) =         85.9  [cm-1]   --> A               I+R\n",
    "     freq ( 31 - 31) =         97.8  [cm-1]   --> E               I+R\n",
    "     freq ( 32 - 32) =         97.8  [cm-1]   --> E*              I+R\n",
    "     freq ( 33 - 33) =        106.4  [cm-1]   --> A               I+R\n",
    "     freq ( 34 - 34) =        109.7  [cm-1]   --> E               I+R\n",
    "     freq ( 35 - 35) =        109.7  [cm-1]   --> E*              I+R\n",
    "     freq ( 36 - 36) =        115.4  [cm-1]   --> A               I+R\n",
    "     freq ( 37 - 37) =        119.3  [cm-1]   --> A               I+R\n",
    "     freq ( 38 - 38) =        121.7  [cm-1]   --> E               I+R\n",
    "     freq ( 39 - 39) =        121.7  [cm-1]   --> E*              I+R\n",
    "     freq ( 40 - 40) =        132.0  [cm-1]   --> A               I+R\n",
    "     freq ( 41 - 41) =        134.9  [cm-1]   --> E               I+R\n",
    "     freq ( 42 - 42) =        134.9  [cm-1]   --> E*              I+R\n",
    "     freq ( 43 - 43) =        135.6  [cm-1]   --> A               I+R\n",
    "     freq ( 44 - 44) =        145.1  [cm-1]   --> E               I+R\n",
    "     freq ( 45 - 45) =        145.1  [cm-1]   --> E*              I+R\n",
    "     freq ( 46 - 46) =        148.9  [cm-1]   --> A               I+R\n",
    "     freq ( 47 - 47) =        162.1  [cm-1]   --> E               I+R\n",
    "     freq ( 48 - 48) =        162.1  [cm-1]   --> E*              I+R\n",
    "     freq ( 49 - 49) =        163.3  [cm-1]   --> A               I+R\n",
    "     freq ( 50 - 50) =        168.6  [cm-1]   --> E               I+R\n",
    "     freq ( 51 - 51) =        168.6  [cm-1]   --> E*              I+R\n",
    "     freq ( 52 - 52) =        172.0  [cm-1]   --> A               I+R\n",
    "     freq ( 53 - 53) =        173.9  [cm-1]   --> E               I+R\n",
    "     freq ( 54 - 54) =        173.9  [cm-1]   --> E*              I+R\n",
    "     freq ( 55 - 55) =        182.4  [cm-1]   --> A               I+R\n",
    "     freq ( 56 - 56) =        189.1  [cm-1]   --> E               I+R\n",
    "     freq ( 57 - 57) =        189.1  [cm-1]   --> E*              I+R\n",
    "     freq ( 58 - 58) =        194.8  [cm-1]   --> A               I+R\n",
    "     freq ( 59 - 59) =        206.4  [cm-1]   --> A               I+R\n",
    "     freq ( 60 - 60) =        211.4  [cm-1]   --> E               I+R\n",
    "     freq ( 61 - 61) =        211.4  [cm-1]   --> E*              I+R\n",
    "     freq ( 62 - 62) =        240.9  [cm-1]   --> E               I+R\n",
    "     freq ( 63 - 63) =        240.9  [cm-1]   --> E*              I+R\n",
    "     freq ( 64 - 64) =        244.0  [cm-1]   --> A               I+R\n",
    "     freq ( 65 - 65) =        251.3  [cm-1]   --> E               I+R\n",
    "     freq ( 66 - 66) =        251.3  [cm-1]   --> E*              I+R\n",
    "     freq ( 67 - 67) =        253.0  [cm-1]   --> A               I+R\n",
    "     freq ( 68 - 68) =        319.9  [cm-1]   --> E               I+R\n",
    "     freq ( 69 - 69) =        319.9  [cm-1]   --> E*              I+R\n",
    "     freq ( 70 - 70) =        371.6  [cm-1]   --> A               I+R\n",
    "     freq ( 71 - 71) =        376.8  [cm-1]   --> E               I+R\n",
    "     freq ( 72 - 72) =        376.8  [cm-1]   --> E*              I+R\n",
    "     freq ( 73 - 73) =        386.1  [cm-1]   --> E               I+R\n",
    "     freq ( 74 - 74) =        386.1  [cm-1]   --> E*              I+R\n",
    "     freq ( 75 - 75) =        404.1  [cm-1]   --> A               I+R\n",
    "     freq ( 76 - 76) =        418.0  [cm-1]   --> A               I+R\n",
    "     freq ( 77 - 77) =        430.2  [cm-1]   --> E               I+R\n",
    "     freq ( 78 - 78) =        430.2  [cm-1]   --> E*              I+R\n",
    "     freq ( 79 - 79) =        433.0  [cm-1]   --> E               I+R\n",
    "     freq ( 80 - 80) =        433.0  [cm-1]   --> E*              I+R\n",
    "     freq ( 81 - 81) =        441.1  [cm-1]   --> A               I+R\n",
    "     freq ( 82 - 82) =        515.5  [cm-1]   --> E               I+R\n",
    "     freq ( 83 - 83) =        515.5  [cm-1]   --> E*              I+R\n",
    "     freq ( 84 - 84) =        518.7  [cm-1]   --> A               I+R\n",
    "     freq ( 85 - 85) =        522.3  [cm-1]   --> A               I+R\n",
    "     freq ( 86 - 86) =        523.6  [cm-1]   --> E               I+R\n",
    "     freq ( 87 - 87) =        523.6  [cm-1]   --> E*              I+R\n",
    "     freq ( 88 - 88) =        537.3  [cm-1]   --> E               I+R\n",
    "     freq ( 89 - 89) =        537.3  [cm-1]   --> E*              I+R\n",
    "     freq ( 90 - 90) =        540.7  [cm-1]   --> A               I+R\n",
    "     freq ( 91 - 91) =        541.7  [cm-1]   --> E               I+R\n",
    "     freq ( 92 - 92) =        541.7  [cm-1]   --> E*              I+R\n",
    "     freq ( 93 - 93) =        546.2  [cm-1]   --> A               I+R\n",
    "     freq ( 94 - 94) =        552.9  [cm-1]   --> A               I+R\n",
    "     freq ( 95 - 95) =        555.0  [cm-1]   --> E               I+R\n",
    "     freq ( 96 - 96) =        555.0  [cm-1]   --> E*              I+R\n",
    "     freq ( 97 - 97) =        558.4  [cm-1]   --> A               I+R\n",
    "     freq ( 98 - 98) =        560.9  [cm-1]   --> E               I+R\n",
    "     freq ( 99 - 99) =        560.9  [cm-1]   --> E*              I+R\n",
    "     freq (100 -100) =        891.1  [cm-1]   --> A               I+R\n",
    "     freq (101 -101) =        895.2  [cm-1]   --> E               I+R\n",
    "     freq (102 -102) =        895.2  [cm-1]   --> E*              I+R\n",
    "     freq (103 -103) =        906.2  [cm-1]   --> E               I+R\n",
    "     freq (104 -104) =        906.2  [cm-1]   --> E*              I+R\n",
    "     freq (105 -105) =        910.7  [cm-1]   --> A               I+R\n",
    "     freq (106 -106) =        912.1  [cm-1]   --> E               I+R\n",
    "     freq (107 -107) =        912.1  [cm-1]   --> E*              I+R\n",
    "     freq (108 -108) =        927.3  [cm-1]   --> A               I+R\n",
    "     freq (109 -109) =        930.7  [cm-1]   --> E               I+R\n",
    "     freq (110 -110) =        930.7  [cm-1]   --> E*              I+R\n",
    "     freq (111 -111) =        945.5  [cm-1]   --> A               I+R\n",
    "     freq (112 -112) =        949.6  [cm-1]   --> E               I+R\n",
    "     freq (113 -113) =        949.6  [cm-1]   --> E*              I+R\n",
    "     freq (114 -114) =        956.4  [cm-1]   --> E               I+R\n",
    "     freq (115 -115) =        956.4  [cm-1]   --> E*              I+R\n",
    "     freq (116 -116) =        958.5  [cm-1]   --> A               I+R\n",
    "     freq (117 -117) =        973.4  [cm-1]   --> A               I+R\n",
    "     freq (118 -118) =        975.5  [cm-1]   --> E               I+R\n",
    "     freq (119 -119) =        975.5  [cm-1]   --> E*              I+R\n",
    "     freq (120 -120) =        994.5  [cm-1]   --> A               I+R\n",
    "     freq (121 -121) =       1002.3  [cm-1]   --> E               I+R\n",
    "     freq (122 -122) =       1002.3  [cm-1]   --> E*              I+R\n",
    "     freq (123 -123) =       1034.2  [cm-1]   --> A               I+R\n",
    "\"\"\""
   ],
   "metadata": {
    "collapsed": false,
    "ExecuteTime": {
     "end_time": "2024-04-05T02:48:32.310195Z",
     "start_time": "2024-04-05T02:48:32.306245Z"
    }
   },
   "id": "4fd58fb668e1e04f",
   "execution_count": 36
  },
  {
   "cell_type": "code",
   "outputs": [],
   "source": [
    "# Split the data into lines for processing\n",
    "lines = data.strip().split(\"\\n\")\n",
    "\n",
    "# Initialize lists to hold the frequencies for A1g and Eg modes\n",
    "a1g_freqs = []\n",
    "eg_freqs = []\n",
    "a_freqs = []\n",
    "e_freqs = []\n",
    "es_freqs = []\n",
    "\n",
    "# Process each line to extract frequencies and symmetries\n",
    "for line in lines:\n",
    "    parts = line.split()\n",
    "    frequency = float(parts[-5])  # Extract the frequency\n",
    "    mode_symmetry = parts[-2]     # Extract the mode symmetry (e.g., A_1g, E_g)\n",
    "    activity = parts[-1]         # Extract the activity label (R for Raman-active, I for Infrared-active)\n",
    "\n",
    "    if activity == \"R\":\n",
    "        if mode_symmetry == \"A_1g\":\n",
    "            a1g_freqs.append(frequency)\n",
    "        elif mode_symmetry == \"E_g\":\n",
    "            eg_freqs.append(frequency)\n",
    "            \n",
    "    elif activity == \"I+R\":\n",
    "        if mode_symmetry == \"A\":\n",
    "            a_freqs.append(frequency)\n",
    "        elif mode_symmetry == \"E\":\n",
    "            e_freqs.append(frequency)\n",
    "        elif mode_symmetry == \"E*\":\n",
    "            es_freqs.append(frequency)\n",
    "\n"
   ],
   "metadata": {
    "collapsed": false,
    "ExecuteTime": {
     "end_time": "2024-04-05T02:48:32.314157Z",
     "start_time": "2024-04-05T02:48:32.311117Z"
    }
   },
   "id": "2212522e8ea2f8a8",
   "execution_count": 37
  },
  {
   "cell_type": "code",
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "eg_freqs: \n",
      "[]\n",
      "a1g_freqs: \n",
      "[]\n",
      "e_freqs: \n",
      "[-25.9, 27.8, 38.9, 47.9, 50.4, 57.3, 62.9, 67.6, 76.6, 83.3, 97.8, 109.7, 121.7, 134.9, 145.1, 162.1, 168.6, 173.9, 189.1, 211.4, 240.9, 251.3, 319.9, 376.8, 386.1, 430.2, 433.0, 515.5, 523.6, 537.3, 541.7, 555.0, 560.9, 895.2, 906.2, 912.1, 930.7, 949.6, 956.4, 975.5, 1002.3]\n",
      "a_freqs: \n",
      "[-56.9, -27.2, -5.6, 47.6, 50.4, 56.6, 60.0, 68.8, 76.4, 85.9, 106.4, 115.4, 119.3, 132.0, 135.6, 148.9, 163.3, 172.0, 182.4, 194.8, 206.4, 244.0, 253.0, 371.6, 404.1, 418.0, 441.1, 518.7, 522.3, 540.7, 546.2, 552.9, 558.4, 891.1, 910.7, 927.3, 945.5, 958.5, 973.4, 994.5, 1034.2]\n",
      "es_freqs: \n",
      "[-25.9, 27.8, 38.9, 47.9, 50.4, 57.3, 62.9, 67.6, 76.6, 83.3, 97.8, 109.7, 121.7, 134.9, 145.1, 162.1, 168.6, 173.9, 189.1, 211.4, 240.9, 251.3, 319.9, 376.8, 386.1, 430.2, 433.0, 515.5, 523.6, 537.3, 541.7, 555.0, 560.9, 895.2, 906.2, 912.1, 930.7, 949.6, 956.4, 975.5, 1002.3]\n"
     ]
    }
   ],
   "source": [
    "print('eg_freqs: ')\n",
    "print(eg_freqs)\n",
    "print('a1g_freqs: ')\n",
    "print(a1g_freqs)\n",
    "print('e_freqs: ')\n",
    "print(e_freqs)\n",
    "print('a_freqs: ')\n",
    "print(a_freqs)\n",
    "print('es_freqs: ')\n",
    "print(es_freqs)"
   ],
   "metadata": {
    "collapsed": false,
    "ExecuteTime": {
     "end_time": "2024-04-05T02:48:32.317482Z",
     "start_time": "2024-04-05T02:48:32.315407Z"
    }
   },
   "id": "59337e5366782035",
   "execution_count": 38
  },
  {
   "cell_type": "code",
   "outputs": [
    {
     "data": {
      "text/plain": "<Figure size 2000x400 with 1 Axes>",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABiIAAAGRCAYAAADo26vuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB0Z0lEQVR4nO3dd3gVVf7H8c/ckpuekEZCr4IFaYoVYQUBC4ooLv5YBctiQRF3bVhxXRTbytpdC1iwoqBiQUSRIoKggEiHBJAWWnq75fz+iLnmkgSSkCEI79fz3AfuzJk53zkzc+7c+83MsYwxRgAAAAAAAAAAADZw1HcAAAAAAAAAAADgyEUiAgAAAAAAAAAA2IZEBAAAAAAAAAAAsA2JCAAAAAAAAAAAYBsSEQAAAAAAAAAAwDYkIgAAAAAAAAAAgG1IRAAAAAAAAAAAANuQiAAAAAAAAAAAALYhEQEAAAAAAAAAAGxDIgIAAKCeTZw4UZZlKSMj47CLo2fPnurZs+chj6W+6q2JHTt26NJLL1ViYqIsy9L48ePrO6QjQkZGhizL0sSJE+s7lCq1aNFCw4YNq+8wamXMmDGyLEu7du2q71AOOcuyNGbMmPoOAwAA4KhEIgIAAKCOXXjhhYqMjFRubm6VZYYMGaKwsDDt3r37EEZ2eFmxYoXGjBlT7wmY2rr11ls1ffp0jR49Wm+++ab69etXZVnLsip9paamHsKIUR2ff/65LMtSo0aNFAgE6jucWnv44Yc1derUeqv/b3/7m8LDw7VmzZoK88aNGyfLsjRt2rR6iAwAAAD1wVXfAQAAABxphgwZok8//VRTpkzRlVdeWWF+QUGBPv74Y/Xr10+JiYm64oorNHjwYHk8nnqIdv+++uor29a9YsUKPfjgg+rZs6datGhxyOqtK998840uuugi3XbbbdUqf84551Q4HiIiIuwI7U+tefPmKiwslNvtrpf6J02apBYtWigjI0PffPONevfuXS9xHKyHH35Yl156qQYMGFAv9f/nP//R559/ruuvv17ffPNNcHp6err+9a9/6ZJLLtEFF1xwSGMqLCyUy8VXYAAAgPrAVRgAAEAdu/DCCxUTE6O333670kTExx9/rPz8fA0ZMkSS5HQ65XQ6D3WY1RIWFnZU1VsTmZmZio+Pr3b5Y445Rn/729+qVdYYo6KioqMyUWFZlsLDw+ul7vz8fH388cd65JFHNGHCBE2aNOlPm4iobykpKXr00Uc1fPhwvf766xo6dKgk6cYbb5Tb7dZ///vfQx5TfR1XAAAA4NFMAAAAdS4iIkIDBw7UzJkzlZmZWWH+22+/rZiYGF144YWSKh+bYdGiRerbt6+SkpIUERGhli1b6uqrrw7OnzVrlizL0qxZs0LWXdnz9ZctW6Zhw4apVatWCg8PV2pqqq6++upqPRZq37EaWrRoUeVjhspi2bhxo2688Ua1a9dOERERSkxM1KBBg0K2b+LEiRo0aJAk6S9/+UuFdVQ2RkRmZqauueYaNWzYUOHh4erYsaNef/31Srf/iSee0P/+9z+1bt1aHo9HJ598sn788ccDbq8kbdiwQYMGDVJCQoIiIyN16qmn6rPPPguJ3bIsGWP03HPPBWM/GC1atNAFF1yg6dOn66STTlJERIReeuklSVJWVpZGjRqlpk2byuPxqE2bNnr00UcrPDYoKytLw4YNU1xcnOLj4zV06FAtWbKkwvFQ1fgbw4YNq3BnSiAQ0Pjx43X88ccrPDxcDRs21HXXXae9e/dWGv/cuXPVrVs3hYeHq1WrVnrjjTcq1JOVlaVbb71VLVq0kMfjUZMmTXTllVcGxyyoaoyIVatW6dJLL1VCQoLCw8N10kkn6ZNPPgkp4/V69eCDD6pt27YKDw9XYmKizjzzTM2YMWN/zR80ZcoUFRYWatCgQRo8eLA++ugjFRUVVWvZAx030h/n7fvvv6+xY8eqSZMmCg8PV69evbRu3boK63zuuefUqlUrRUREqFu3bpozZ061xk+xLEv5+fl6/fXXg8fnvmNalB0v8fHxiouL01VXXaWCgoIK63rrrbfUtWtXRUREKCEhQYMHD9bmzZur1SbXXnutzjjjDN12223avXu33n33XX355Zf697//rcaNGx9w+SeeeEKnn366EhMTFRERoa5du2ry5MkhZSZMmCDLsvTaa6+FTH/44YdlWZY+//zzkHYpP0ZEbm6uRo0aFTwWU1JSdM455+inn36q1vYBAACg+rgjAgAAwAZDhgzR66+/rvfff1833XRTcPqePXs0ffp0XX755VX+tXtmZqb69Omj5ORk3XXXXYqPj1dGRoY++uijWsUyY8YMbdiwQVdddZVSU1P166+/6n//+59+/fVX/fDDDzX6EX38+PHKy8sLmfbUU09pyZIlSkxMlCT9+OOP+v777zV48GA1adJEGRkZeuGFF9SzZ0+tWLFCkZGROuusszRy5Eg9/fTTuvvuu3XsscdKUvDffRUWFqpnz55at26dbrrpJrVs2VIffPCBhg0bpqysLN1yyy0h5d9++23l5ubquuuuk2VZeuyxxzRw4EBt2LBhv4/82bFjh04//XQVFBRo5MiRSkxM1Ouvv64LL7xQkydP1sUXX6yzzjpLb775pq644opKH7dUlaKiogoDBMfExAQfybV69Wpdfvnluu666/T3v/9d7dq1U0FBgXr06KEtW7bouuuuU7NmzfT9999r9OjR2rZtW3CAbGOMLrroIs2dO1fXX3+9jj32WE2ZMiX4V+i1dd1112nixIm66qqrNHLkSKWnp+vZZ5/Vzz//rHnz5oW05bp163TppZfqmmuu0dChQ/Xaa69p2LBh6tq1q44//nhJUl5enrp3766VK1fq6quvVpcuXbRr1y598skn+u2335SUlFRpHL/++qvOOOMMNW7cWHfddZeioqL0/vvva8CAAfrwww918cUXSyodiPmRRx7Rtddeq27duiknJ0eLFi3STz/9pHPOOeeA2ztp0iT95S9/UWpqqgYPHqy77rpLn376aTBpVpXqHDfljRs3Tg6HQ7fddpuys7P12GOPaciQIVqwYEGwzAsvvKCbbrpJ3bt316233qqMjAwNGDBADRo0UJMmTfYbz5tvvhlsg+HDh0uSWrduHVLmsssuU8uWLfXII4/op59+0iuvvBK8i6HM2LFjdd999+myyy7Ttddeq507d+qZZ57RWWedpZ9//vmAdwRZlqWXXnpJnTt31g033KA5c+bopJNO0ogRI/a7XJn//ve/uvDCCzVkyBCVlJTo3Xff1aBBgzRt2jSdf/75kqSrrrpKH330kf7xj3/onHPOUdOmTfXLL7/owQcf1DXXXKPzzjuvyvVff/31mjx5sm666SYdd9xx2r17t+bOnauVK1eqS5cu1YoRAAAA1WQAAABQ53w+n0lLSzOnnXZayPQXX3zRSDLTp08PTpswYYKRZNLT040xxkyZMsVIMj/++GOV6//222+NJPPtt9+GTE9PTzeSzIQJE4LTCgoKKiz/zjvvGElm9uzZVcZhjDE9evQwPXr0qDKO999/30gy//rXv/Zb3/z5840k88YbbwSnffDBB5VuQ2X1jh8/3kgyb731VnBaSUmJOe2000x0dLTJyckJ2f7ExESzZ8+eYNmPP/7YSDKffvppldtijDGjRo0yksycOXOC03Jzc03Lli1NixYtjN/vD06XZEaMGLHf9ZUvW9mrbD81b97cSDJffvllyHIPPfSQiYqKMmvWrAmZftdddxmn02k2bdpkjDFm6tSpRpJ57LHHgmV8Pp/p3r17heOhqn06dOhQ07x58+D7OXPmGElm0qRJIeW+/PLLCtPL4i9/PGVmZhqPx2P++c9/Bqfdf//9RpL56KOPKtQfCASMMZUfw7169TIdOnQwRUVFIeVPP/1007Zt2+C0jh07mvPPP7/Cuqtjx44dxuVymZdffjk47fTTTzcXXXRRhbLNmzc3Q4cODb6v7nFTdt4ee+yxpri4OFj2v//9r5FkfvnlF2OMMcXFxSYxMdGcfPLJxuv1BstNnDjRSNrvOVkmKioqJMYyDzzwgJFkrr766pDpF198sUlMTAy+z8jIME6n04wdOzak3C+//GJcLleF6fszevRoI8k4nU6zePHiai+3b19SUlJiTjjhBHP22WeHTN+2bZtJSEgw55xzjikuLjadO3c2zZo1M9nZ2SHlJJkHHngg+D4uLq7a5zAAAAAODo9mAgAAsIHT6dTgwYM1f/78kEcSvf3222rYsKF69epV5bJlf2U8bdo0eb3eg46l/J0XZX+Vf+qpp0rSQT2CZMWKFbr66qt10UUX6d577620Pq/Xq927d6tNmzaKj4+vdX2ff/65UlNTdfnllwenud1ujRw5Unl5efruu+9Cyv/1r39VgwYNgu+7d+8uqfTxOQeqp1u3bjrzzDOD06KjozV8+HBlZGRoxYoVtYpfki666CLNmDEj5NW3b9/g/JYtW4a8l6QPPvhA3bt3V4MGDbRr167gq3fv3vL7/Zo9e3YwbpfLpRtuuCG4rNPp1M0331zreD/44APFxcXpnHPOCam7a9euio6O1rfffhtS/rjjjgu2syQlJyerXbt2IW3+4YcfqmPHjhXuEJBU5Z05e/bs0TfffKPLLrtMubm5wTh2796tvn37au3atdqyZYuk0nPn119/1dq1a2u8ve+++64cDocuueSS4LTLL79cX3zxRYVHUe2rpsfNVVddFTIOyr7H56JFi7R79279/e9/DxlceciQISHH9cG4/vrrQ953795du3fvVk5OjiTpo48+UiAQ0GWXXRay/1NTU9W2bdsK+39/yu50adSokU444YRqL1e+L9m7d6+ys7PVvXv3Cv1IamqqnnvuOc2YMUPdu3fXkiVL9Nprryk2Nna/64+Pj9eCBQu0devWascEAACA2iERAQAAYJOywajffvttSdJvv/2mOXPmaPDgwfsdnLpHjx665JJL9OCDDyopKUkXXXSRJkyYoOLi4lrFsWfPHt1yyy1q2LChIiIilJycrJYtW0qSsrOza7XOnJwcDRw4UI0bN9Ybb7wR8iNyYWGh7r///uCYBklJSUpOTlZWVlat69u4caPatm0rhyP08rXsUU4bN24Mmd6sWbOQ92U/3h7oB+WNGzeqXbt2FaZXVU9NNGnSRL179w55paWlBeeX7ZPy1q5dqy+//FLJyckhr7IBlMvGINm4caPS0tIUHR0dsnxl21Jda9euVXZ2tlJSUirUn5eXV2H8k33bXCpt9/Jtvn79+hr9EC2VPvLJGKP77ruvQhwPPPCApD/a4V//+peysrJ0zDHHqEOHDrr99tu1bNmyatXz1ltvqVu3btq9e7fWrVundevWqXPnziopKdEHH3yw32Vretwc6PgsK9+mTZuQci6Xq8I4HrV1oBjWrl0rY4zatm1bod1XrlxZ6fg3ldm8ebMeeOABnXDCCdq8ebMee+yxkPl79uzR9u3bg6/yfcS0adN06qmnKjw8XAkJCUpOTtYLL7xQaT8yePBgnX/++Vq4cKH+/ve/7zfZW+axxx7T8uXL1bRpU3Xr1k1jxow5YLISAAAAtcMYEQAAADbp2rWr2rdvr3feeUd333233nnnHRljggmKqliWpcmTJ+uHH37Qp59+qunTp+vqq6/Wk08+qR9++EHR0dFV/vW43++vMO2yyy7T999/r9tvv12dOnVSdHS0AoGA+vXrV2HA4+oaNmyYtm7dqoULF1b4q+Obb75ZEyZM0KhRo3TaaacpLi5OlmVp8ODBta6vpqpK9BhjDkn9tVHZmCGBQEDnnHOO7rjjjkqXOeaYY2pcT9lA2/va99gJBAJKSUnRpEmTKl1PcnJyyHu72rzsmLntttsq3DFSpuwH+7POOkvr16/Xxx9/rK+++kqvvPKKnnrqKb344ou69tprq6xj7dq1wcHM27ZtW2H+pEmTgmMt1IXD4fg8UAyBQECWZemLL76otOy+Sa+qlI2R88UXX+gf//iHxo4dq//7v/9Tq1atJEkDBw4MuaNp6NChmjhxoubMmaMLL7xQZ511lp5//nmlpaXJ7XZrwoQJweRuebt379aiRYskld6tFQgEKiQu93XZZZepe/fumjJlir766is9/vjjevTRR/XRRx/p3HPPrdb2AQAAoHpIRAAAANhoyJAhuu+++7Rs2TK9/fbbatu2rU4++eRqLXvqqafq1FNP1dixY/X2229ryJAhevfdd3XttdcG/3o5KysrZJl9//J67969mjlzph588EHdf//9wem1eXRNmXHjxmnq1Kn66KOP1L59+wrzJ0+erKFDh+rJJ58MTisqKqoQa00GyW7evLmWLVtW4cfFVatWBefXhebNm2v16tUVptd1PdXVunVr5eXlBe+AqErz5s01c+ZM5eXlhfxAXNm2NGjQoNK/+t732GndurW+/vprnXHGGVUOrF5TrVu31vLly2u0TNkP1m63+4DtIEkJCQm66qqrdNVVVykvL09nnXWWxowZs99ExKRJk+R2u/Xmm29W+NF97ty5evrpp7Vp06ZK7/qQ6v64KSu/bt06/eUvfwlO9/l8ysjI0IknnnjAddTk/KpM69atZYxRy5Yta5XwkqQpU6bok08+0VNPPaUmTZpo/Pjxmj59ukaMGKEvvvhCkvTkk0+G3DXTqFEjSaWP8QoPD9f06dODA7pL0oQJEyqta8SIEcrNzdUjjzyi0aNHa/z48frHP/5xwBjT0tJ044036sYbb1RmZqa6dOmisWPHkogAAACoYzyaCQAAwEZldz/cf//9WrJkyQHvhpBKkwf7/mV0p06dJCn4eKbmzZvL6XQGxwgo8/zzz4e8L/tRdd/1jR8/vtrbUN7XX3+te++9V/fcc48GDBhQaRmn01mhvmeeeabCX9xHRUVJqphMqcx5552n7du367333gtO8/l8euaZZxQdHa0ePXrUbEP2U8/ChQs1f/784LT8/Hz973//U4sWLXTcccfVST3Vddlll2n+/PmaPn16hXlZWVny+XySSuP2+Xx64YUXgvP9fr+eeeaZCsu1bt1aq1at0s6dO4PTli5dqnnz5lWo2+/366GHHqqwDp/PV639tq9LLrlES5cu1ZQpUyrMq+pugJSUFPXs2VMvvfSStm3bVmF++e3YvXt3yLzo6Gi1adPmgI81mzRpkrp3766//vWvuvTSS0Net99+uyTpnXfeqXL5uj5uTjrpJCUmJurll18O7uOyOA/0eLEyUVFRtdpHZQYOHCin06kHH3ywwr4xxlRo633l5uZq5MiR6ty5c3CskkaNGumhhx7Sl19+GXzcVdeuXUMeV1bWVk6nU5ZlhfQbGRkZmjp1aoW6Jk+erPfee0/jxo3TXXfdpcGDB+vee+/VmjVrqozP7/dXeMRTSkqKGjVqVOvH4AEAAKBq3BEBAABgo5YtW+r000/Xxx9/LEnVSkS8/vrrev7553XxxRerdevWys3N1csvv6zY2Fidd955kqS4uDgNGjRIzzzzjCzLUuvWrTVt2rQKz22PjY3VWWedpccee0xer1eNGzfWV199pfT09Fptz+WXX67k5GS1bdtWb731Vsi8c845Rw0bNtQFF1ygN998U3FxcTruuOM0f/58ff3110pMTAwp36lTJzmdTj366KPKzs6Wx+PR2WefrZSUlAr1Dh8+XC+99JKGDRumxYsXq0WLFpo8ebLmzZun8ePHKyYmplbbs6+77rpL77zzjs4991yNHDlSCQkJev3115Wenq4PP/zwgI96qWu33367PvnkE11wwQUaNmyYunbtqvz8fP3yyy+aPHmyMjIylJSUpP79++uMM87QXXfdpYyMDB133HH66KOPKn2W/tVXX63//Oc/6tu3r6655hplZmbqxRdf1PHHHx8cqFgqHavkuuuu0yOPPKIlS5aoT58+crvdWrt2rT744AP997//1aWXXlrj7Zk8ebIGDRqkq6++Wl27dtWePXv0ySef6MUXX1THjh0rXe65557TmWeeqQ4dOujvf/+7WrVqpR07dmj+/Pn67bfftHTpUkmlA2b37NlTXbt2VUJCghYtWqTJkycHHw9UmQULFmjdunVVlmncuLG6dOmiSZMm6c4776y0TF0fN2FhYRozZoxuvvlmnX322brsssuUkZGhiRMnqnXr1tW626Fr1676+uuv9Z///EeNGjVSy5Ytdcopp1Q7htatW+vf//63Ro8erYyMDA0YMEAxMTFKT0/XlClTNHz4cN12221VLn/vvfdq69at+uijj0LuMhkxYoRef/11jRo1Sv369avy3D3//PP1n//8R/369dP//d//KTMzU88995zatGkTMu5HZmambrjhBv3lL38J7sNnn31W3377rYYNG6a5c+dW2v65ublq0qSJLr30UnXs2FHR0dH6+uuv9eOPP4bczQUAAIA6YgAAAGCr5557zkgy3bp1q3T+hAkTjCSTnp5ujDHmp59+Mpdffrlp1qyZ8Xg8JiUlxVxwwQVm0aJFIcvt3LnTXHLJJSYyMtI0aNDAXHfddWb58uVGkpkwYUKw3G+//WYuvvhiEx8fb+Li4sygQYPM1q1bjSTzwAMPVBmHMcb06NHD9OjRI/heUpWvb7/91hhjzN69e81VV11lkpKSTHR0tOnbt69ZtWqVad68uRk6dGjINrz88sumVatWxul0hqxj33qNMWbHjh3B9YaFhZkOHTqEbKcxxqSnpxtJ5vHHH6/Qzvtub1XWr19vLr30UhMfH2/Cw8NNt27dzLRp0ypd34gRIw64vuqUbd68uTn//PMrnZebm2tGjx5t2rRpY8LCwkxSUpI5/fTTzRNPPGFKSkqC5Xbv3m2uuOIKExsba+Li4swVV1xhfv755wrHgzHGvPXWW6ZVq1YmLCzMdOrUyUyfPt0MHTrUNG/evEL9//vf/0zXrl1NRESEiYmJMR06dDB33HGH2bp16wHjr2w/7t6929x0002mcePGJiwszDRp0sQMHTrU7Nq1yxjzxz7cN+b169ebK6+80qSmphq3220aN25sLrjgAjN58uRgmX//+9+mW7duJj4+3kRERJj27dubsWPHhrTTvm6++WYjyaxfv77KMmPGjDGSzNKlS4Pbu++xXJ3j5ttvvzWSzAcffBAyvaptfvrpp03z5s2Nx+Mx3bp1M/PmzTNdu3Y1/fr1qzLWMqtWrTJnnXWWiYiIMJKC8T7wwANGktm5c2dI+crOf2OM+fDDD82ZZ55poqKiTFRUlGnfvr0ZMWKEWb16dZV1L1q0yDidTnPTTTdVOn/hwoXG4XCYkSNH7ncbXn31VdO2bVvj8XhM+/btzYQJE4Lxlxk4cKCJiYkxGRkZIct+/PHHRpJ59NFHg9PK9wHFxcXm9ttvNx07djQxMTEmKirKdOzY0Tz//PP7jQkAAAC1YxlzGI/YBwAAAKDWMjIy1LJlS02YMEHDhg2r73BwkAKBgJKTkzVw4EC9/PLL9R0OAAAAUG2MEQEAAAAAh5mioqIKYzO88cYb2rNnj3r27Fk/QQEAAAC1xBgRAAAAAHCY+eGHH3Trrbdq0KBBSkxM1E8//aRXX31VJ5xwggYNGlTf4QEAAAA1QiICAAAAAA4zLVq0UNOmTfX0009rz549SkhI0JVXXqlx48YpLCysvsMDAAAAaoQxIgAAAAAAAAAAgG0YIwIAAAAAAAAAANiGRAQAAAAAAAAAALBNtcaICAQC2rp1q2JiYmRZlt0xAQAAAAAAAACAw5gxRrm5uWrUqJEcjv3f81CtRMTWrVvVtGnTOgkOAAAAAAAAAAAcGTZv3qwmTZrst0y1EhExMTHBFcbGxh58ZDiqeL1effXVV+rTp4/cbnd9hwOgHtAPAJDoCwDQDwCgHwBAP3AkycnJUdOmTYP5g/2pViKi7HFMsbGxJCJQY16vV5GRkYqNjaVzAY5S9AMAJPoCAPQDAOgHANAPHImqM5wDg1UDAAAAAAAAAADbkIgAAAAAAAAAAAC2IREBAAAAAAAAAABsU60xIgAAAAAAAAAAhy+/3y+v11vfYRyQ1+uVy+VSUVGR/H5/fYeD/XC73XI6nXWyLhIRAAAAAAAAAPAnZYzR9u3blZWVVd+hVIsxRqmpqdq8eXO1BjlG/YqPj1dqaupB7ysSEQAAAAAAAADwJ1WWhEhJSVFkZORh/+N+IBBQXl6eoqOj5XAwcsDhyhijgoICZWZmSpLS0tIOan0kIgAAAAAAAADgT8jv9weTEImJifUdTrUEAgGVlJQoPDycRMRhLiIiQpKUmZmplJSUg3pME3saAAAAAAAAAP6EysaEiIyMrOdIcKQqO7YOdvwREhEAAAAAAAAA8Cd2uD+OCX9edXVskYgAAAAAAAAAAAC2IREBAAAAAAAAAMBhKiMjQ5ZlacmSJfUdSq2RiAAAAAAAAAAAHDI33nijnE6nLMsKefXr16/KZcaMGVOhvGVZ+vrrrw9h5PWjadOm2rZtm0444YT6DqXWXPUdAAAAAAAAAADg0NuVV6xiX6DK+R6XQ0nRHlvq7tu3ryZOnBhan2f/dR1//PEVEg8JCQkVypWUlCgsLOygYzxcOJ1Opaam1ncYB4U7IgAAAAAAAADgKLMrr1hjPvlVoz9cVuVrzCe/aldesS31ezwepaamhrwaNGiw32VcLleFZcLCwjRs2DANGDBAY8eOVaNGjdSuXTtJ0ubNm3XZZZcpPj5eCQkJuuiii5SRkRFcn9/v1z/+8Q/Fx8crMTFRd9xxh4YOHaoBAwYEy7Ro0ULjx48PiaNTp04aM2ZM8H1WVpauvfZaJScnKzY2VmeffbaWLl0anD9mzBh16tRJb775plq0aKG4uDgNHjxYubm5wTKBQECPPfaY2rRpI4/Ho2bNmmns2LGSKn800/Lly3XuuecqOjpaDRs21BVXXKFdu3YF50+ePFkdOnRQRESEEhMT1bt3b+Xn5x9ot9iGRAQAAAAAAAAAHGWKfQHlFHrlcTsVHxlW4eVxO5VT6N3vHROHk5kzZ2r16tWaMWOGpk2bJq/Xq759+yomJkZz5szRvHnzFB0drX79+qmkpESS9OSTT2rixIl67bXXNHfuXO3Zs0dTpkypcd2DBg1SZmamvvjiCy1evFhdunRRr169tGfPnmCZ9evXa+rUqZo2bZqmTZum7777TuPGjQvOHz16tMaNG6f77rtPK1as0Ntvv62GDRtWWl9WVpbOPvtsde7cWYsWLdKXX36pHTt26LLLLpMkbdu2TZdffrmuvvpqrVy5UrNmzdLAgQNljKnxttUVHs0EAAAAAAAAAEepCLdTUZ7KfyYu9vptq/ezzz5TdHR0yLS7775bd999d5XL/PLLLyHLHHfccVq4cKEkKSoqSq+88krwkUxvvfWWAoGAXnnlFVmWJUmaMGGC4uPjNWvWLPXp00fjx4/X6NGjNXDgQEnSiy++qOnTp9doO+bOnauFCxcqMzMz+GipJ554QlOnTtXkyZM1fPhwSaV3PEycOFExMTGSpCuuuEIzZ87U2LFjlZubq//+97969tlnNXToUElS69atdeaZZ1Za57PPPqvOnTvr4YcfDk577bXX1LRpU61Zs0Z5eXny+XwaOHCgmjdvLknq0KFDjbarrpGIAAAAAAAAAAAcUj179tSLL74YMq2y8R7Ka9eunT755JPg+/JjSnTo0CFkXIilS5dq3bp1wR/+yxQVFWn9+vXKzs7Wtm3bdMoppwTnuVwunXTSSTW6c2Dp0qXKy8tTYmJiyPTCwkKtX78++L5FixYhsaSlpSkzM1OStHLlShUXF6tXr17VrvPbb7+tkMiRSu+86NOnj3r16qUOHTqob9++6tOnjy699NIDPvrKTiQiAAAAAAAAAACHVFRUlNq0aVOjZcLCwqpcJioqKuR9Xl6eunbtqkmTJlUom5ycXO06HQ5HhcSE1+sNqSctLU2zZs2qsGx8fHzw/263O2SeZVkKBEofexUREVHteMrq7N+/vx599NEK89LS0uR0OjVjxgx9//33+uqrr/TMM8/onnvu0YIFC9SyZcsa1VVXGCMCAAAAAAAAAHBE6dKli9auXauUlBS1adMm5BUXF6e4uDilpaVpwYIFwWV8Pp8WL14csp7k5GRt27Yt+D4nJ0fp6ekh9Wzfvl0ul6tCPUlJSdWKtW3btoqIiNDMmTOrvW2//vqrWrRoUaHOsoSMZVk644wz9OCDD+rnn39WWFhYrca/qCskIgAAAAAAAAAAh1RxcbG2b98e8tq1a1edrX/IkCFKSkrSRRddpDlz5ig9PV2zZs3SyJEj9dtvv0mSbrnlFo0bN05Tp07VqlWrdOONNyorKytkPWeffbbefPNNzZkzR7/88ouGDh0qp9MZnN+7d2+ddtppGjBggL766itlZGTo+++/1z333KNFixZVK9bw8HDdeeeduuOOO/TGG29o/fr1+uGHH/Tqq69WWn7EiBHas2ePLr/8cv34449av369pk+frquuukp+v18LFizQww8/rEWLFmnTpk366KOPtHPnTh177LG1a8w6wKOZAAAAAAAAAOAoVVjFgNRVTa8r06dPV1paWsi0du3aadWqVXWy/sjISM2ePVt33nmnBg4cqNzcXDVu3Fi9evVSbGysJOmf//yntm3bpqFDh8rhcOjqq6/WxRdfrOzs7OB6Ro8erfT0dF1wwQWKi4vTQw89FHJHhGVZ+vzzz3XPPffoqquu0s6dO5WamqqzzjpLDRs2rHa89913n1wul+6//35t3bpVaWlpuv766yst26hRI82bN0933nmn+vTpo+LiYjVv3lz9+vWTw+FQbGysZs+erfHjxysnJ0fNmzfXk08+qXPPPbeWrXnwLFONkTdycnIUFxen7Ozs4E4Cqsvr9erzzz/XeeedV+FZaACODvQDACT6AgD0AwDoB4C6VlRUpPT0dLVs2VLh4eE1WnZXXrHGfPKrcgq9VZaJjXBrzIXHKynaU2WZmgoEAsrJyVFsbKwcjsPvgT3Dhg1TVlaWpk6dWt+hHBb2d4zVJG/AHREAAAAAAAAAcJRJivZozIXHq9gXqLKMx+Wo0yQEjl4kIgAAAAAAAADgKESSAYcKiQgAAAAAAAAAACRNnDixvkM4Ih1+D+ECAAAAAAAAAABHDBIRAAAAAAAAAADANiQiAAAAAAAAAACAbUhEAAAAAAAAAAAA25CIAAAAAAAAAAAAtiERAQAAAAAAAAAAbEMiAgAAAAAAAACAw9jEiRMVHx9f32HUGokIAAAAAAAAAMAhN3/+fDmdTp1//vnVKt+zZ09ZllXh5fP5bI60/v31r3/VmjVr6juMWiMRAQAAAAAAAABHs0BA2rVO2rK49N9A4JBU++qrr+rmm2/W7NmztXXr1mot8/e//13btm0LeblcrgrlSkpK6jrcehUREaGUlJT6DqPWSEQAAAAAAAAAwNFq21Jp+t3S57dL0+/9/d+7S6fbKC8vT++9955uuOEGnX/++Zo4cWK1louMjFRqamrIS5JatGihhx56SFdeeaViY2M1fPhwSdLcuXPVvXt3RUREqGnTpho5cqTy8/OD68vMzFT//v0VERGhli1batKkSWrRooXGjx8vScrIyJBlWVqyZElwmaysLFmWpVmzZgWnLV++XOeee66io6PVsGFDXXHFFdq1a1dwfs+ePTVy5EjdcccdSkhIUGpqqsaMGROybVlZWbruuuvUsGFDhYeH64QTTtC0adMkVf5opo8//lhdunRReHi4WrVqpQcffDB4d4gxRmPGjFGzZs3k8XjUqFEjjRw5slptbAcSEQAAAAAAAABwNNq2VPruUWnbz1JkAymxdem/237+fbp9yYj3339f7du3V7t27fS3v/1Nr732mowxB7XOJ554Qh07dtTPP/+s++67T+vXr1e/fv10ySWXaNmyZXrvvfc0d+5c3XTTTcFlhg0bps2bN+vbb7/V5MmT9fzzzyszM7NG9WZlZenss89W586dtWjRIn355ZfasWOHLrvsspByr7/+uqKiorRgwQI99thj+te//qUZM2ZIkgKBgM4991zNmzdPb731llasWKFx48bJ6XRWWuecOXN05ZVX6pZbbtGKFSv00ksvaeLEiRo7dqwk6cMPP9RTTz2ll156SWvXrtXUqVPVoUOHGm1XXap4zwoAAAAAAAAA4MgWCEhL3pEKdktJ7STLKp3uiSl9v2u1tPRdqWEHyVH3f88+YcIE/e1vf5Mk9evXT9nZ2fruu+/Us2fP/S73/PPP65VXXgm+v+666/Tkk09Kks4++2z985//DM679tprNWTIEI0aNUqS1LZtWz399NPq0aOHXnjhBW3atElffPGFFi5cqJNPPllS6eOijj322Bpty7PPPqvOnTvr4YcfDk577bXX1LRpU61Zs0bHHHOMJOnEE0/UAw88EIzl2Wef1cyZM3XOOefo66+/1sKFC7Vy5cpg+VatWlVZ54MPPqi77rpLQ4cODZZ96KGHdMcdd+iBBx7Qpk2blJqaqt69e8vtdqtZs2bq1q1bjbarLpGIAAAAAAAAAICjzZ4N0s7VUmzjP5IQZSyrdHrmqtJySW3qtOq1a9dq4cKFmjJliiTJ5XLpr3/9q1599dUDJiKGDBmie+65J/i+/OOKTjrppJCyS5cu1bJlyzRp0qTgNGOMAoGA0tPTtWbNGrlcLnXt2jU4v3379hUegXQgS5cu1bfffqvo6OgK89avXx+SiCgvLS0tePfFkiVL1KRJk2DZ6tQ5b9684B0QkuT3+1VUVKSCggINGjRI48ePV6tWrdSvXz+dd9556t+/f6XjaRwKJCIAAAAAAAAA4GhTnC35iiR3ZOXz3RGSb1tpuTr25ptvyufzqVGjRsFpxhh5PB49++yziouLq3LZuLg4tWlTeWIkKioq5H1eXp6uu+66SsdGaNasmdasWXPAWB2/3w1S/rFRXq+3Qj39+/fXo48+WmH5tLS04P/dbnfIPMuyFPh9YPCIiIgDxrJvnQ8++KAGDhxYYV54eLiaNm2q1atX6+uvv9aMGTN044036vHHH9d3331XIY5DgUQEAAAAAAAAABxtPHGSK1zyFpQ+jmlf3sLS+Z6qkwK14fP59N577+mJJ55Q3759Q+YNGDBA77zzjq6//vo6qatLly5asWJFlYmL9u3by+fzafHixcFHM61evVpZWVnBMsnJyZKkbdu2qXPnzpIUMnB1WT0ffvihWrRoUes7Dk488UT99ttvIY9y2p8uXbpo9erVVW6bVJrc6N+/v/r3768RI0aoffv2+uWXX9SlS5daxXgwSEQAAAAAAAAAwNEmoZWU3K50YOryY0RIkjFSzhapUZfScnVo2rRpysrK0tVXX60GDRqEzLvkkkv06quv1lki4s4779Spp56qm266Sddee62ioqK0YsUKzZgxQ88++6zatWunfv366brrrtMLL7wgl8ulUaNGhdydEBERoVNPPVXjxo1Ty5YtlZmZqXvvvTeknhEjRujll1/W5ZdfrjvuuEMJCQlat26d3n33Xb3yyitVDjhdXo8ePXTWWWfpkksu0X/+8x+1adNGq1atkmVZ6tevX4Xy999/vy644AI1a9ZMl156qRwOh5YuXarly5fr3//+tyZOnCi/369TTjlFkZGReuuttxQREaHmzZsffMPWQt2PMgIAAAAAAAAAOLw5HFKny6XIxNKBqYtzpYCv9N9dq6XIJKnj4DofqPq1115Tjx49Kn380iWXXKJFixZp2bJldVLXiSeeqO+++05r1qxR9+7d1blzZ91///0hj4SaMGGCGjVqpB49emjgwIEaPny4UlJSKsTs8/nUtWtXjRo1Sv/+979D5jdq1Ejz5s2T3+9Xnz591KFDB40aNUrx8fHBRztVx4cffqiTTz5Zl19+uY477jjdcccd8vv9lZbt27evpk2bpq+++konn3yyTj31VD311FPBREN8fLxefvllnXHGGTrxxBP19ddf69NPP1ViYmK146lLlin/cKsq5OTkKC4uTtnZ2YqNjT0UceEI4vV69fnnn+u8886rl+ePAah/9AMAJPoCAPQDAOgHgLpWVFSk9PR0tWzZUuHh4bVbybal0pJ3Sgeu9hWVPo4ppX1pEiKtY90GLCkQCCgnJ0exsbE1+pH+UGrRooVGjRqlUaNG1Xco9W5/x1hN8gY8mgkAAAAAAAAAjlZpHaWGHaQ9G0oHpvbElT6O6TBNEuDPiUQEAAAAAAAAABzNHA4pqepBj4GDRSICAAAAAAAAAIDfZWRk1HcIRxzurwEAAAAAAAAAALYhEQEAAAAAAAAAAGxDIgIAAAAAAAAAANiGRAQAAAAAAAAAALANiQgAAAAAAAAAAGAbEhEAAAAAAAAAAMA2JCIAAAAAAAAAAPiTsixLU6dOre8w9otEBAAAAAAAAADgkLnxxhvldDplWVbIq1+/frVa38SJEyusy7IsvfLKK3Uc+eFp27ZtOvfcc+s7jP1y1XcAAAAAAAAAAID6ZYzR1vytahTVSJZl2V5f3759NXHixJBpHo+n1uuLjY3V6tWrQ6bFxcVVKFdSUqKwsLBa13M4Sk1Nre8QDog7IgAAAAAAAADgKLds1zL9d/F/9cuuXw5JfR6PR6mpqSGvBg0aSCpNiowZM0bNmjWTx+NRo0aNNHLkyP2uz7KsCuuLiIjQmDFj1KlTJ73yyitq2bKlwsPDJUlZWVm69tprlZycrNjYWJ199tlaunRpyDrHjRunhg0bKiYmRtdcc43uuusuderUKTi/Z8+eGjVqVMgyAwYM0LBhw4Lvi4uLddttt6lx48aKiorSKaecolmzZgXnT5w4UfHx8Zo+fbqOPfZYRUdHq1+/ftq2bVvIel977TUdf/zx8ng8SktL00033RSy7eUfzbR582Zddtllio+PV0JCgi666CJlZGQE58+aNUvdunVTVFSU4uPjdcYZZ2jjxo37bd+DRSICAAAAAAAAAI5iARPQjIwZWpu1VjM2zlDABOo1ng8//FBPPfWUXnrpJa1du1ZTp05Vhw4dar2+devW6cMPP9RHH32kJUuWSJIGDRqkzMxMffHFF1q8eLG6dOmiXr16ac+ePZKk999/X2PGjNHDDz+sRYsWKS0tTc8//3yN677ppps0f/58vfvuu1q2bJkGDRqkfv36ae3atcEyBQUFeuKJJ/Tmm29q9uzZ2rRpk2677bbg/BdeeEEjRozQ8OHD9csvv+iTTz5RmzZtKq3P6/Wqb9++iomJ0Zw5czRv3rxgcqOkpEQ+n08DBgxQjx49tGzZMs2fP1/Dhw+3/S4YHs0EAAAAAAAAAEexX3b9opV7VqpxdGOt2L1Cy3ct14nJJ9pa52effabo6OiQaXfffbfuvvtubdq0Sampqerdu7fcbreaNWumbt267Xd92dnZIeuLjo7W9u3bJZU+jumNN95QcnKyJGnu3LlauHChMjMzg4+DeuKJJzR16lRNnjxZw4cP1/jx43XNNdfommuukST9+9//1tdff62ioqJqb+OmTZs0YcIEbdq0SY0aNZIk3Xbbbfryyy81YcIEPfzww5JKkwcvvviiWrduLak0efGvf/0ruJ5///vf+uc//6lbbrklOO3kk0+utM733ntPgUBAr7zySjC5MGHCBMXHx2vWrFk66aSTlJ2drQsuuCBY37HHHlvtbaotEhEAAAAAAAAAcJQquxvCF/ApOTJZG7I2aMbGGToh6QQ5LPseqNOzZ0+9+OKLIdMSEhIkld6tMH78eLVq1Ur9+vXTeeedp/79+8vlqvrn7JiYGP3000/B9w7HH7E3b948mISQpKVLlyovL0+JiYkh6ygsLNT69eslSStXrtT1118fMv+0007Tt99+W+1t/OWXX+T3+3XMMceETC8uLg6pOzIyMpgUkKS0tDRlZmZKkjIzM7V161b16tWrWnUuXbpU69atU0xMTMj0oqIirV+/Xn369NGwYcPUt29fnXPOOerdu7cuu+wypaWlVXu7aoNEBAAAAAAAAAAcpcruhkiNKh3wuGFUw0NyV0RUVFSVjxdq2rSpVq9era+//lozZszQjTfeqMcff1zfffed3G53pcs4HI4q1xcVFRXyPi8vT2lpaSFjNZSJj4+v9jY4HA4ZY0Kmeb3ekHqcTqcWL14sp9MZUq783Rv7bpNlWcH1RkREVDuesjq7du2qSZMmVZhXloyZMGGCRo4cqS+//FLvvfee7r33Xs2YMUOnnnpqjeqqCcaIAAAAAAAAAICjUPm7ISLdkZKkKHeUfAFfvY8VERERof79++vpp5/WrFmzNH/+fP3yS90MpN2lSxdt375dLpdLbdq0CXklJSVJKn1c0YIFC0KW++GHH0LeJycnhwwq7ff7tXz58uD7zp07y+/3KzMzs0I9qamp1Yo1JiZGLVq00MyZM6u9bWvXrlVKSkqFOuPi4kJiGz16tL7//nudcMIJevvtt6u1/toiEQEAAAAAAAAAR6F974YoU/6uCLsUFxdr+/btIa9du3ZJkiZOnKhXX31Vy5cv14YNG/TWW28pIiJCzZs3r5O6e/furdNOO00DBgzQV199pYyMDH3//fe65557tGjRIknSLbfcotdee00TJkzQmjVr9MADD+jXX38NWc/ZZ5+tzz77TJ999plWrVqlG264QVlZWcH5xxxzjIYMGaIrr7xSH330kdLT07Vw4UI98sgj+uyzz6od75gxY/Tkk0/q6aef1tq1a/XTTz/pmWeeqbTskCFDlJSUpIsuukhz5sxRenq6Zs2apZEjR+q3335Tenq6Ro8erfnz52vjxo366quvtHbtWtvHieDRTAAAAAAAAABwlCm7G8Ib8CrcFR5y90OEK0LegNfWsSKmT59eYVyCdu3aadWqVYqPj9e4ceP0j3/8Q36/Xx06dNCnn35aYUyH2rIsS59//rnuueceXXXVVdq5c6dSU1N11llnqWHDhpKkv/71r1q/fr3uuOMOFRUV6ZJLLtENN9yg6dOnB9dz9dVXa+nSpbryyivlcrl066236i9/+UtIXRMmTAgONr1lyxYlJSXp1FNP1QUXXFDteIcOHaqioiI99dRTuu2225SUlKRLL7200rKRkZGaPXu27rzzTg0cOFC5ublq3LixevXqpdjYWBUWFmrVqlV6/fXXtXv3bqWlpWnEiBG67rrratGS1WeZfR9iVYmcnBzFxcUpOztbsbGxtgaEI4/X69Xnn3+u8847r8pnuAE4stEPAJDoCwDQDwCgHwDqWlFRkdLT09WyZUuFh4fXaNmM7AyN/2m89hbtrbJMg/AGGtVllFrEtTjISP8QCASUk5Oj2NjYkAGl/wzGjBmjqVOnasmSJfUdyiGzv2OsJnkD7ogAAAAAAAAAgKNMk5gm+vuJf5c/4K+yjNPhVJOYJocwKhypSEQAAAAAAAAAwFHG5XDp+MTj6zsMHCX+XPe+AAAAAAAAAABQD8aMGXNUPZapLpGIAAAAAAAAAAAAtiERAQAAAAAAAAB/YsaY+g4BR6i6OrZIRAAAAAAAAADAn5Db7ZYkFRQU1HMkOFKVHVtlx1ptMVg1AAAAAAAAAPwJOZ1OxcfHKzMzU5IUGRkpy7LqOar9CwQCKikpUVFRkRwO/k7+cGWMUUFBgTIzMxUfHy+n03lQ6yMRAQAAAAAAAAB/UqmpqZIUTEYc7owxKiwsVERExGGfNIEUHx8fPMYOBokIAAAAAAAAAPiTsixLaWlpSklJkdfrre9wDsjr9Wr27Nk666yzDvpxP7CX2+0+6DshypCIAAAAAAAAAIA/OafTWWc/GtvJ6XTK5/MpPDycRMRRhIdwAQAAAAAAAAAA25CIAAAAAAAAAAAAtiERAQAAAAAAAAAAbEMiAgAAAAAAAAAA2IZEBAAAAAAAAAAAsA2JCAAAAAAAAAAAYBsSEQAAAAAAAAAAwDYkIgAAAAAAAAAAgG1IRAAAAAAAAAAAANuQiAAAAAAAAAAAALYhEQEAAAAAAAAAAGxDIgIAAAAAAAAAANiGRAQAAAAAAAAAALANiQgAAAAAAAAAAGAbEhEAAAAAAAAAAMA2JCIAAAAAAAAAAIBtSEQAAAAAAAAAAADbkIgAAAAAAAAAAAC2IREBAAAAAAAAAABsQyICAAAAAAAAAADYhkQEAAAAAAAAAACwDYkIAAAAAAAAAABgGxIRAAAAAAAAAADANiQiAAAAAAAAAACAbUhEAAAAAAAAAAAA25CIAAAAAAAAAAAAtiERAQAAAAAAAAAAbEMiAgAAAAAAAAAA2IZEBAAAAAAAAAAAsA2JCAAAAAAAAAAAYBsSEQAAAAAAAAAAwDYkIgAAAAAAAAAAgG1IRAAAAAAAAAAAANuQiAAAAAAAAAAAALYhEQEAAAAAAAAAAGxDIgIAAAAAAAAAANiGRAQAAAAAAAAAALANiQgAAAAAAAAAAGAbEhEAAAAAAAAAAMA2JCIAAAAAAAAAAIBtSEQAAAAAAAAAAADbkIgAAAAAAAAAAAC2IREBAAAAAAAAAABsQyICAAAAAAAAAADYhkQEAAAAAAAAAACwDYkIAAAAAAAAAABgGxIRAAAAAAAAAADANiQiAAAAAAAAAACAbUhEAAAAAAAAAAAA25CIAAAAAAAAAAAAtiERAQAAAAAAAAAAbEMiAgAAAAAAAAAA2IZEBAAAAAAAAAAAsA2JCAAAAAAAAAAAYBsSEQAAAAAAAAAAwDYkIgAAAAAAAAAAgG1IRAAAAAAAAAAAANuQiAAAAAAAAAAAALYhEQEAAAAAAAAAAGxDIgIAAAAAAAAAANiGRAQAAAAAAAAAALANiQgAAAAAAAAAAGAbEhEAAAAAAAAAAMA2JCIAAAAAAAAAAIBtSEQAAAAAAAAAAADbkIgAAAAAAAAAAAC2IREBAAAAAAAAAABsQyICAAAAAAAAAADYhkQEAAAAAAAAAACwDYkIAAAAAAAAAABgGxIRAAAAAAAAAADANiQiAAAAAAAAAACAbVz1HcCfViAg7dkgFWdLnjgpoZXkqP+8zq68YhV7fXJlpUtF2dquEiU2PFmWwylJ8rgcSor21G2lh2lblLcrr1jFvoAkyRijHYXb1DAiTZZlSaq6XWq6XPn2z83arc2+QsUkdpEnd5P25G9UamxLRTduJ1kOeVylbVRhf6V0lTtnk/Ky98jvjlFUo2NkZAXrzirwqtgXkMflUHykWzsKtqhxcZEcJbkynlg5k1orKSai4jaYgJx7NygzN10pMS3lb9AqNI59tzO8odzZGbKKc+SKjFdCk/YH3K+1audyx88ef4QKY5pLlqPq5SPdFY83SXt+WyVfQZaMJ1beuBbaUbSjWvu4fAxZe3ep2BUtX3xLGWO0e8ePSlWYFB4nX3xLedyuCuvY7zabgCJyNyrBWVjjc6PCeivZz7Icldcto/yta+T05io6LqF67bGf87ja+7W2fUEly+36/Tg/YJ2HWJVtISNXVro8vjzFN0g6cPvV4vxC1Wrbxx/R6uC8OiTtWlfXELVcT/lt3JNXrC15W5UUnlo6r2i7Gkc3UsLv23igz7DD9RroSFTdvnhXWJr8u9NlFecoEBajLZ5wNYxsXONjuKr6svKLpT0bFBHIV1Rsg9D1V+caIBCo+vqlNp/lAADUl5p+r6/F7y/5RSWKK9x80J/rdR13Xa5fqpvfZvYXz0H9dlLLba/sWmh3fokkaVt2kRxO7+H/HY7r/jpBIqI2ti2Vlrwj7Vwt+YokV7iU3E7qdLmU1rHewtqVV6xX3p+qTnu+UGPvRqW7CjUt2quTiltps2eANoa1UWyEW2MuPL7uTubDtC3K25VXrDGf/KqcQq8kKV/rlemYqZRAb0Wp9EfsytqlpsuVb/+0ko1aZeXp69gSXZhjVCKXvogN6C85kWrgbKv5MX20PeoYNS1eq27ZXwX31xdRRRqQL7Us8SjXa6lEYVoX2VKfxR2vZWGr1cB7trZmpsgfMHI5LHVJ/EkFjs80MDegY0skr+XWzoiW6tz/RjVofVJwGxpkr9QZ+TNUoHX6KqpAffMjFaE2mhd1jrZHHSNJKvYGgtuZr6kamGt0WuEuuY1XxuVRePsuijz5b1Xu11q1c7njp6SkQDv2+LTR0VTzos7RijCrwvLHWRkambRIkVnr/jjeIhNU5Atox6ZNsnzF8lpufR+RqCkxDkXp4v3u4/IxlGxfqZ2Ze1QUcCnHEad0d5FmR+xU/zy3WvgitMXdXEsSztW1lw0IrmN/23xcSUBn5M9Q88BmRSe4FBYWWe1zY9/1xvi+U46mhezn7Z4Wmh9zjtY5WofUfWJJe52fvVwNCtLlUYnywiL0Q+QB2mM/5/GumPbV2q8PnRJQgzWTa94XVFJ3QXwbvbLrJK0wLap/LB0CVe3vE0vaq3/ur2rs3ahwh0+RKQkKSz22yvarzfmFqtW2jz+i1cF5dUjata6uIWq5nvLbWOILaE32cvmiZ8uRc5YkKRA7W+68Hmobf7zCnI79foYdrtdAR6Lq9sWeGJcKsnJULJf8cmtFmDQlxqFY9Veuq3QfV+cYrqq+4wqPUY+dS9XCbFaJvFoS5dLU2NL1JwQaHfgaYNtSFfz4lnas+qnC9Utr38k1/ywHAKC+1PB7fW1+f2lUuEad9nypY11b5ZG31p/rdR13Xa7f4y77Y9WD+21mf/FUde1R59cXB7hO3pVXrEc+X6mzo6QxnyxXtkk/vL/Dcd1fZ0hE1NS2pdJ3j0oFu6XYxpI7UvIWSNt+lrI3ST3urLeDMLB1qXpnTlSsyVFWWLJmRRltcnvVMLBBA/Mn6VPXFVpX2CqYXT1oh3FblFfsCyin0CuP26lwl6Vd/iXyBraq0LVEjZzHqMhnlFPordAuNV2ufPvvdifp+7ACbXd79UtksbxyaZM7QvMjpZsK1io1d6dmF52h04vmKF65ygpL1ryIYm1xl2iJp1hdinL0m1qpRGFq5V+r6MBayQpXkWeJpHPkcjjU1qxXnPlYvznz9X1UrJIdjeX2F6lxwWqFz3tcirxbxRHHqEH2Sl2cO0kxgWz9L86h9DCHvpel4dlrlZKbqXd9g5UR1lZxEW6FuywVlXwry7dWm90OdfI31V7jkctfpFbbl0jfbalyv9a4nfc5frzhqdq7e6taedeqYX6mVkc0ltf6Y/nk/DXqvfsNOX0+KaFp6fGWs1VaP1PugJFXbZQV1ljhVrE2u9NlmYBcrm8UF1b1Pi4fgy88Vb9ZHiVZu9WheJFmR0ToN3eUvo2J1eD8OLUqWaukzB0KbG0pHdNtv9vssr7RwPwtivJma5cjSd7YRgpzFFf73Ci/3jb+9doTmKoMd+h+blK4Sv2Kt+nLpKHaGXWMdvmXSP50xQRWqpVfylC8ChWuVLd//+1xgPM40GWUcgrNfvdrg+yVCp/3ieTLqrovSDqu4oZWUbdz+xL13rNcJYlXBretOues3Srb32Vt3tpnaY8rWbtMmBp7IhRWRfvV9vxC1Wrbxx+x6ui8sr1ddyyX5j5+8NcQB3EtUn4b3S4pULJMCtsuR/QySVIgbLv8EcsU7ekgy3Lu9zPscL0GOhJVpy8u8hepWe4GJXlzVeiI0Kbw9poflas9zmw19k5RM0czrXO2rtYxXFV9cVqpdsannUrSdnn0fdRO7XFmq3nxexqSFylHwF/1NYAkffeonDk7tddEyxeWFLx+cRifGvg2qJU/rPqf5QAA1Jcafq+vze8v3XwO9cp7W47AHuU7U5XrjtD8iC21+lyv67jrcv3ZBSWSZR30bzP7i6eya486v76oxnVyccQxyi3ySlFSbLhL2w/n73Bc99cp7iGpiUCgNANWsFtKaid5YiSHs/TfpHal05e+W1quHmKLXPG+YgI52ulpoXXhDm12FSspEK5VnjDtcOxV9/yvZZk6iu1wbosqRLid8ro3qsDaqChnigqsjfK6NynC7Tz45fZtf4+lnZ58NfIFtDA8QssiHGri92uzx69FUQ0V58/WhfmTFRPIDpbf6spVY5/Rck+E1oZZaqodKrAi9GNUin5zFqqpz6tCa6Msz2Z5HEZd3V/oN1eBGpgoZbh92uwpkc8VrS2upnIU7fm9/X06I3+G4gI5+jE6RRlun5IDHmW4fVoUnaK4QI7OKpwpywQU4XbK58qQ37dCjfzSKk+Y1oc75HS6VeSIlK/BMdXar9Vtr8qOnyJHpHZ6Wmi7c6/8/hWKciarwNoonytD3fNnKCaQUxqHJ0ayHFLuVsnhlhwuNfDvktPp0rpwh1Z7wtTIL/n9K+RzZVS+j/eJwYTFyMihhMBuLfd4tCrMpcY+vza7irQ+3KGdnhaKCeQocuUHFbY/dJuT5fet0A7HXu30tFCRI7LW50aky1Ja4Uf6zVlYYT//5mymmECOuud/LZ8rQwXaqKY+rzY7C7UoqqEKrEgZy7H/9jAHPo8jV34QPD4q26+RLktn5M8oPeZq0hfspw/xNTgmdNtqeM7aLdgW5ds8OkU+V7SM5ZQJq7z96uL8QtVq28cfUWw4r2xr12WTD/4aoo6uRUrPz00yYVvk9CcpEL5eJnyDXP4kmbDf5HNvrtZn2OF8DXQkqrIvdkapoW+bLOPXXkeCXDLa69iqDLdPDUyUNjsL1ahwiiJd1kHV95urQPM8SSpyROq38BJtDPOXrt+Vrz3WNu30NK/8GmDJO9KSt0v/CKLBMSpyRMrpdP/+ee1WM2+JNrsKtCgqpXqf5QAA1Jcafq+vze8vhcpQo8IpijU52qDG8rqitMlTcnCf6zbEXZfrr/PvjuXiqXjtUcfXF9W9Ti7326TPvfnw/Q7HdX+dIxFRE3s2lN6GE9tYsvbp5CyrdHrmqtJy9RCba88a7XYmKWBJi11ZCsgozrjll9G3kS6llKSroW9rndV32LZFFYwJaJt3oYwJKNwRL2P82uZdIHOA5Ey1ltun/Ze5d0nyK9bvUJZL2u1wKMXnlSyffvJkK88VpST/TuVbUQpY0lLXLhnjV7xxyG9ZmhkZrmjlKVIF+smTrWKHS6m+YrlMgXwRS9XIbNXWyN9UbLkVZ8Lklwnuc1mW/NGNpMxV8myao8bejdrtTtJid3bIMbHYna3d7iQ19W1Sqn+rjAlob/FMuUyhohQZuk6p2vu1uu1V1fETsKRZES65AoWKN24Z49fe4plKKcnQbmfSH+WLc6SibCksUoGwaEUG8hQWyAvGHKVIOQKF2ls8s/J9XEkMESZf4YFczYwKl89yKMXvk5GvdJ2WtNuZJNfu1RW2v/w2xxu3HKZQ30a6FNj3WqiG50aid7NWWxkqdrgq3c+7nUlKKUn/fb8VKNVXrCKHSz+Fle7rgIx+Cqu6PVxZ6Qc8j127V6uhb2uV+zXRu1mNvRtLj7n99QV7Mw7Y/uWXK79tNT1n7VbWFuXbvOz8Cqqk/eri/ELVatvHH1FsOK9sa9ddaw/+GqKOrkWMCWhn4EdJAVmBGAUcOfI7smUFYiTLr0z/j9X+DONcPnSq6ovDAnmKCuQp4I6ULEtFjnDNDyuSkU9xJkxFDpdWW+lK9G4+qPqKLbdWROXJp4B+jcyVX0YJAUt+SbPCHQoz+aErKDs2tvwsbV0ScvwEfv9MMPIrKWDkl9Ey965qfZYDAFBvavi9vja/vzhNgVZb6drtSpJkBT8zS3/bqOXnug1x1+X66/y7YyXxlG/HOr2+qOZ1sisrPTh5u3fx4fsdjuv+OkcioiaKs0ufBeaOrHy+O6J0fnH2oY1LkoqzZfmLVWx5lOEs0CZnoeIDbklSg4BbG1xebXQVKiKQf4AVVb++w7YtqpBjNig7kK5wR4IkKdzRQNmBDOWY9INfrkL7FyvFF9BeZ+kXUp9lKcspNfBZ2uQs1HqXT075FZBTGc4CbXQWKcUfkJFDSX5pTZiltR5LWz0F2uwsUnzAI4fxKyYQpYD7NxVGLtLmMK/iAqXPymsQcGujs1AbnQWSJOOKlHxFcuTvkNt4tcZtKhwTG52FWus2cpsSRQbylWM2aG9gs5L8UsByVVinpGrt1+q2V1XHT4azQBtcXiX5JafxKdzRQHsDm7XRVahiq9yzAf1eKeCXHC7JcsmhgDaVO/YDlkvJfmlvYHPl+7iSGFzGp1VuaY3boaSAQ5aMEv2OYDuUWB5Z/uIK219+m53Gp2S/Ss+58m1XgzYsk+9fp3RnieKr2M/FlkcZrkLtDWxWjImWw/jVIODRJlehtoYVaYun6Pfjp/L2sIpzDngeW/5iRfx+fFS2X/P9635/VuUB+oKSnNDpB+hDSsptW03PWbuVtUX5Nq9wrkgV2q8uzi9UrbZ9/BHFhvPKtnb118E1RB1di+SYDco1GbICcTKOPEkBSQEZR54sf7xyTfU/w2oUPw5KVX3xJmeBHArIWKVPoP01zKk1YZaS/KVfexoEPEp3lijfv+6g6osLhGlrWLGWROVoW1ixGgRKr0VS/EZrwhzaVNU1gDe/9Hb+csdP2bV7WYwpvoA2uoqr9VkOAEC9qeH3+tr8/hJjopTuLNEad+kP0xtdhfv8tlGLz3Ub4q7L9df5d8dK4smoyW8nNVHN62Sr+I/fB3ICGYfvdziu++sciYia8MSVDkjireSLhSR5C0vne+IObVyS5ImTcXrkNkXBv5oKV+mtTOFyKiC/5ob7VeCIqLP6Dtu2qIRR6V87GhOQywqXJLmsCBnjL52uKv4StLrL7dP+PkvyGEuZLslS6Ym2zeWQ27jkl9FCT4G8ckryBctHGEuWAoo0kt8y+jbSrcWRJfLLKNIYBSynLEeUjOXT9qhV8spSpCnNiIfLKb9M8C/hLV+B5ApXIKqhii2Xlrr3VDgm/DJa4t6jYsutfEekdgZ+lN+SPHLJYXwV1inpgPu1Ju1V2fFT9ld/Rn555JLfcsllRchvSXPD/XKboj8KO92lt8QFfJLxySeHfvAUBrfTYXy/r0OV7+NKYiixnJoV6VZAAUWagIys0nX83g4uUyTj9IRs/77b7Ldc8silgPyhbVemmueGUUDrnOvlt6rez25TpLnhfvktyeGIVMByKtKU/iXl8qic4F9oVtUexhN7wPPYOD0qcERUuV/XO9er2HKVHnNVrEOucCks9oDtX56r3LbV5Jy1W/n9vW+bV9jf+7TfwZ5fqFpt+/gjTh2fV7a2q7MOriHq4FokuI0KSCZMAeduSZYsWfI7d8syYTLyVeszrMbxo9b21xcvCC+UTw5ZxqeApHnhRn5ZCvt9aLxIY+S1LK13rq/2MVx5fZLfCujH6CwFrMDv19sORRrJZ0kLfr8mCeEtlNxRfzxbWH9c+wRkysVoySdV67McAIB6U8Pv9bX5/cXhiJLXsvRL2F4FFNCS378/lf22UZvP9bqOuy7XX9YGdfrdcZ94yl971Pn1RTWvk40nNliPMebw/Q7HdX+dIxFREwmtSkdFz9kimX2+WBhTOj2lfWm5eojNl3CM9jp2aHO57PDvwSnV59fqMLfS3SV1Vt9h2xaVKFC6cs0fWdYy4Y4GyjUZKlDGwS23T/vHBTza5XQpxyF5AlJEIKAsp1NZDksN/C5tchVpfkSStrlzg+V9lltO45eMUUOvXz+Hh2lduFcN/C6FBYpU6IxRkSNKlglTkStbfhMhjymSfv8gahBwa5Or9C8BnXlbpZT2Km7WXT9EJGmLM1fxgdCx6RsEXNrizNX88CSlu4qVazLkdKSp0BmjsEDpeoPrdBVWa7/WpL0qO342uQq12Vmohj5/cHslyelI0+owl/Y6dvxR3hMrhcdJJQVylORpSXiU0t2+3499E2wzpyOt8n1cSQyrwxxaGeZWQ59fzoBPPssd/AuETc4C7XXskC+xXcj277vNRY4oFTpjlOrzl96h4Sr8o84anBsFSlemtVuRig3uj/L7ebOrQHudmVod5pLTkRasNyxQpAZ+lzI8hdoYXqAGAVeV7eGLb3nA89iX2E7p7pIq9+sOa7cWRCSXHnP76wsatDhg+5dfbq9jR3Db9q1zf+es3crv75A2D7j+OFekStvvYM8vVK22ffwRp47PK1vbNantwV9D1MG1SNk2epQgOfJK74IwYZIJk3Hkyzjy5FFitT7Dahw/am1/ffEGl08/h0fJ4S3QarfROndACQGnAtYfn4dRitUOa3e1j+HK6vMEiuQOOLTH7ZUrUPqVym+5ZCQl+402uHyVXwM07iw16hQ8fko/00vveChb3lhSXMBTrc9yAADqTQ2/19fm95ciR5SiFKstrjwtic7RZlf537tq97le13HX1fqLrI0qsjLq/rvjPvGUv/ao8+uLal4n++JbqlAbJUkeR4OQYofVdziu++sciYiacDikTpdLkYnSrtVScW7pX2MX55a+j0ySOg4uLVcPseUde6lmRzrkMIUKD/gVMH4pUKIwX64shWmPM1Z7HT8rUBfPWjuc22IfARPQXsdi+Y1Pxu+W1+cPvow/TH7jVZZjcYV2qdFy+7S/x/i12VU6RoBbfjkklchSlqtIMf58FVtOvRHbRrMjyvZXQPmOCPlU+herYUba5bSU5/Aqxp+vIsut7e7GKgn4FXDkSApoqzNMRXIpzJdbup+NFDA+rQ/bIl94gtRxsAKWQ19Gx6nYcirWVyAFSoLHRayvQMWWU19ExyjH9VPpdgbCtN3dWEWWW2G+XIUH/PIZv5aH7ZZj7/73a03ba9/jJ+D3ldZjCmWpNA6vP1C6fCBMe5xxpe279/fjzfilmEZSwKuA8WlmZLh8JqDwgF9hvtxgm5lAFft4nxhMcY6WerJVZLkVbgIKyKjQ8sj/+zodplCzIx3Ka39JcPsr3WZ/QNvdjWUpTA5TqOVhuxXw1+zcKL/eLHez4P4I3c+/aXaEpb3OOJlAWLDeIsutGH+e8p1e5Tp88uyvPaQDnsd57S/RHsdP+9mvPn0RHSefp0HN+oL99CGOvas1O9Lxx7ZV85y1W4X9Xa7NY30FChiflnmyZIpzKm+/gzi/ULXa9vFHpDo8r2xv1xMvPfhriIO8Fim/jQHjkt+5W0ZGAWMpYCwZ+UunGdcBP8MO12ugI9EB+2L5NSMqSl7LqUURRfJalgJWlEzAG/w8zHI3K/2MrcYxXFV9hZZTxc4CGRnlOr0yv19z5zuiZKyoqq8BOl0udfo/KTJRjr2rtTxst3zGH7x+yXdEKd8RVf3P8qOhbwMAHJ5q+r2+Nr+/+AOl34nl0PKYnfIbn8KMgr911fRzva7jrsv1ZzsXKce5qO6/O5aLp7Jrjzq9vqjmdXJAUpa1RJJk/K7D9zsc1/11znXgIgiR1lHqcWfpqOk7V0u+baW34TTqUnrwpXWst9C2xSZoeWRjuXxbtMNRKEtGxrJU7AxXjjNexYpQhHOXdhVtUVMdc/AVHsZtUd6uoi0yzt0K+KUcVT5Yd6CSdqnpcuXb/zdnoXIdfvksh4qcjt//lt0o1+XTakekCh0J8oU5lelPVqR/V3B/ZTo8suRWkWXJbxkZy2h1WLgKHQ1UbPLkD+yWnKXPnit2FWmJK0YJ/nx5TOnyDmMpwxOj9V2u0IlpHbVr2xplhfkV8KWoRFnBcuWPC5+zQFKBAoHS7cyRtCesgWL9peXdxijb5dfG5I5qf+rwKvdrjdt5n+Nne0mOsl2FKrQi9UtYvIpNnuTPK7d0hJZHRmljfIra527743hr01sbivO1M3eN3MavHQ7vH8d8uXVUto/Lx7Bt+1LlO0qfU7gyrPS2Oqd8spRb2l5WpJZHNta22ASlHmCb/2hDS9muQm3PXauYsNhqnxvl17tNeSH7o2w/rwmPUJYzUYVyq9C/NaTeCP9u+a0iWVZAvzmLpP21xwHO422KOOB+zQrzaX2XK3Xi5rlV9wVeb8UFq6h7Y3I7Lc/aoUIFgtu232PpEKlsf4eeL0XKc+RoW8lOta2i/Wp7fqFqte3jj1h1dF7Z3q4NT6iba4iDuBYpv415ZuPvn69+yVHu7lFntnIDm+RU2H4/ww7Xa6AjUXX64r1OnxbEJWuDcuWzLO2wimRkhV4fmOodw1XVlxkWqxKVyJJfBS6vfvOb4OetJMX6s/Z/DdDjTm384SVlZy6ocP0iqWaf5UdD3wYAODzV+Ht9zX9/yZHkdCYooN1KCvi1q+w7ei0+1+s67rpcv99dut0H+9vM/uKp6tqjTq8vqnGdvGvbGhnnHklSrn+7SlTx94LD5jqH6/46ZRmz770lFeXk5CguLk7Z2dmKjY09UPGjQyBQOip6cXbps8ASWtV7BswX8GnBluUq8hbLmbtNljdfxh0lf0yaym5+CXe7dUrjE+Ry1GEO6gBt4fV69fnnn+u8886T2+3ez4rs8Ue7VPJD6O8qa5eaLle+/ZXzmzbsXqciK0y+8AS5ivbI4SuSyxOllo1PktPhltvpkCWpxO8L3V/RqVLuFm3YvU4BV0SwvCT5jV9r96SrxO9XmNOptgkt5bSskOXdDZrplCYn7hOTV1Kg0uPijzj2zTT/UT4sIk4nt+sjlyusztu57PjxFe7RjzmZKohIVlU3a4W73Tol7Ti5sjaFHG++gE8/rv5KJYXZFY75/dZdLgbfrrVavGWJCh2e35dXhfYKd3tC1nHgbQ4osnCnTo5NkSsiodr9ROXrDd1/jrjGsixHJftN8huv0n9bJMtfpDYJrWVimxy4Pao4j2u0X+Wosi/Ybz+wT92++GZasG1FzY8lm+2/LUr3T0SgWF0bd5Irqe0B2q9m5xeqVuu+50h3kOeVXe1aoS+oq+upWqyn/DaW/3wt74/PWud+P8MOp+vBI111++KOaR20uLhQ3qzNlV4TS9U7hquqz2/8Wrt7gwJ5O+SRt5LP2wNfA/h8JVVev9Tqs/xPor6/GwCof/QDR6Cafq+v4e8vfuPXxqyNau10yekrqvXnel3HXZfrr6vfZvZnf9ceZers+mI/18m+gE/zNy3TnqWbFXVcqixnxTY57K5zuO6vUk3yBiQiYDsuMgDQDwCQ6AsA0A8AoB8AQD9wJKlJ3oDUDQAAAAAAAAAAsA2JCAAAAAAAAAAAYBsSEQAAAAAAAAAAwDYkIgAAAAAAAAAAgG1IRAAAAAAAAAAAANuQiAAAAAAAAAAAALYhEQEAAAAAAAAAAGxDIgIAAAAAAAAAANiGRAQAAAAAAAAAALANiQgAAAAAAAAAAGAbEhEAAAAAAAAAAMA2JCIAAAAAAAAAAIBtSEQAAAAAAAAAAADbkIgAAAAAAAAAAAC2IREBAAAAAAAAAABsQyICAAAAAAAAAADYhkQEAAAAAAAAAACwDYkIAAAAAAAAAABgGxIRAAAAAAAAAADANiQiAAAAAAAAAACAbUhEAAAAAAAAAAAA25CIAAAAAAAAAAAAtiERAQAAAAAAAAAAbEMiAgAAAAAAAAAA2IZEBAAAAAAAAAAAsA2JCAAAAAAAAAAAYBsSEQAAAAAAAAAAwDYkIgAAAAAAAAAAgG1IRAAAAAAAAAAAANuQiAAAAAAAAAAAALYhEQEAAAAAAAAAAGxDIgIAAAAAAAAAANiGRAQAAAAAAAAAALANiQgAAAAAAAAAAGAbEhEAAAAAAAAAAMA2JCIAAAAAAAAAAIBtSEQAAAAAAAAAAADbkIgAAAAAAAAAAAC2IREBAAAAAAAAAABsQyICAAAAAAAAAADYhkQEAAAAAAAAAACwDYkIAAAAAAAAAABgGxIRAAAAAAAAAADANiQiAAAAAAAAAACAbUhEAAAAAAAAAAAA25CIAAAAAAAAAAAAtiERAQAAAAAAAAAAbEMiAgAAAAAAAAAA2IZEBAAAAAAAAAAAsA2JCAAAAAAAAAAAYBsSEQAAAAAAAAAAwDYkIgAAAAAAAAAAgG1IRAAAAAAAAAAAANuQiAAAAAAAAAAAALYhEQEAAAAAAAAAAGxDIgIAAAAAAAAAANiGRAQAAAAAAAAAALANiQgAAAAAAAAAAGAbEhEAAAAAAAAAAMA2JCIAAAAAAAAAAIBtSEQAAAAAAAAAAADbkIgAAAAAAAAAAAC2IREBAAAAAAAAAABsQyICAAAAAAAAAADYhkQEAAAAAAAAAACwDYkIAAAAAAAAAABgGxIRAAAAAAAAAADANiQiAAAAAAAAAACAbUhEAAAAAAAAAAAA25CIAAAAAAAAAAAAtiERAQAAAAAAAAAAbEMiAgAAAAAAAAAA2IZEBAAAAAAAAAAAsA2JCAAAAAAAAAAAYBsSEQAAAAAAAAAAwDYkIgAAAAAAAAAAgG1IRAAAAAAAAAAAANuQiAAAAAAAAAAAALYhEQEAAAAAAAAAAGxDIgIAAAAAAAAAANiGRAQAAAAAAAAAALANiQgAAAAAAAAAAGAbEhEAAAAAAAAAAMA2JCIAAAAAAAAAAIBtSEQAAAAAAAAAAADbkIgAAAAAAAAAAAC2IREBAAAAAAAAAABsQyICAAAAAAAAAADYhkQEAAAAAAAAAACwDYkIAAAAAAAAAABgGxIRAAAAAAAAAADANiQiAAAAAAAAAACAbUhEAAAAAAAAAAAA25CIAAAAAAAAAAAAtiERAQAAAAAAAAAAbEMiAgAAAAAAAAAA2IZEBAAAAAAAAAAAsA2JCAAAAAAAAAAAYBsSEQAAAAAAAAAAwDYkIgAAAAAAAAAAgG1IRAAAAAAAAAAAANuQiAAAAAAAAAAAALYhEQEAAAAAAAAAAGxDIgIAAAAAAAAAANiGRAQAAAAAAAAAALANiQgAAAAAAAAAAGAbEhEAAAAAAAAAAMA2JCIAAAAAAAAAAIBtSEQAAAAAAAAAAADbkIgAAAAAAAAAAAC2IREBAAAAAAAAAABsQyICAAAAAAAAAADYhkQEAAAAAAAAAACwDYkIAAAAAAAAAABgGxIRAAAAAAAAAADANiQiAAAAAAAAAACAbUhEAAAAAAAAAAAA25CIAAAAAAAAAAAAtiERAQAAAAAAAAAAbEMiAgAAAAAAAAAA2IZEBAAAAAAAAAAAsA2JCAAAAAAAAAAAYBsSEQAAAAAAAAAAwDYkIgAAAAAAAAAAgG1IRAAAAAAAAAAAANuQiAAAAAAAAAAAALYhEQEAAAAAAAAAAGxDIgIAAAAAAAAAANiGRAQAAAAAAAAAALANiQgAAAAAAAAAAGAbEhEAAAAAAAAAAMA2JCIAAAAAAAAAAIBtSEQAAAAAAAAAAADbkIgAAAAAAAAAAAC2IREBAAAAAAAAAABsQyICAAAAAAAAAADYhkQEAAAAAAAAAACwDYkIAAAAAAAAAABgGxIRAAAAAAAAAADANiQiAAAAAAAAAACAbUhEAAAAAAAAAAAA25CIAAAAAAAAAAAAtiERAQAAAAAAAAAAbEMiAgAAAAAAAAAA2IZEBAAAAAAAAAAAsA2JCAAAAAAAAAAAYBtXdQoZYyRJOTk5tgaDI5PX61VBQYFycnLkdrvrOxwA9YB+AIBEXwCAfgAA/QAA+oEjSVm+oCx/sD/VSkTk5uZKkpo2bXoQYQEAAAAAAAAAgCNJbm6u4uLi9lvGMtVIVwQCAW3dulUxMTGyLKvOAsTRIScnR02bNtXmzZsVGxtb3+EAqAf0AwAk+gIA9AMA6AcA0A8cSYwxys3NVaNGjeRw7H8UiGrdEeFwONSkSZM6CQ5Hr9jYWDoX4ChHPwBAoi8AQD8AgH4AAP3AkeJAd0KUYbBqAAAAAAAAAABgGxIRAAAAAAAAAADANiQiYDuPx6MHHnhAHo+nvkMBUE/oBwBI9AUA6AcA0A8AoB84WlVrsGoAAAAAAAAAAIDa4I4IAAAAAAAAAABgGxIRAAAAAAAAAADANiQiAAAAAAAAAACAbUhEAAAAAAAAAAAA25CIgK2ee+45tWjRQuHh4TrllFO0cOHC+g4JQB155JFHdPLJJysmJkYpKSkaMGCAVq9eHVKmqKhII0aMUGJioqKjo3XJJZdox44dIWU2bdqk888/X5GRkUpJSdHtt98un893KDcFQB0ZN26cLMvSqFGjgtPoB4Cjw5YtW/S3v/1NiYmJioiIUIcOHbRo0aLgfGOM7r//fqWlpSkiIkK9e/fW2rVrQ9axZ88eDRkyRLGxsYqPj9c111yjvLy8Q70pAGrB7/frvvvuU8uWLRUREaHWrVvroYcekjEmWIZ+ADiyzJ49W/3791ejRo1kWZamTp0aMr+uzvlly5ape/fuCg8PV9OmTfXYY4/ZvWmwCYkI2Oa9997TP/7xDz3wwAP66aef1LFjR/Xt21eZmZn1HRqAOvDdd99pxIgR+uGHHzRjxgx5vV716dNH+fn5wTK33nqrPv30U33wwQf67rvvtHXrVg0cODA43+/36/zzz1dJSYm+//57vf7665o4caLuv//++tgkAAfhxx9/1EsvvaQTTzwxZDr9AHDk27t3r8444wy53W598cUXWrFihZ588kk1aNAgWOaxxx7T008/rRdffFELFixQVFSU+vbtq6KiomCZIUOG6Ndff9WMGTM0bdo0zZ49W8OHD6+PTQJQQ48++qheeOEFPfvss1q5cqUeffRRPfbYY3rmmWeCZegHgCNLfn6+OnbsqOeee67S+XVxzufk5KhPnz5q3ry5Fi9erMcff1xjxozR//73P9u3DzYwgE26detmRowYEXzv9/tNo0aNzCOPPFKPUQGwS2ZmppFkvvvuO2OMMVlZWcbtdpsPPvggWGblypVGkpk/f74xxpjPP//cOBwOs3379mCZF154wcTGxpri4uJDuwEAai03N9e0bdvWzJgxw/To0cPccsstxhj6AeBoceedd5ozzzyzyvmBQMCkpqaaxx9/PDgtKyvLeDwe88477xhjjFmxYoWRZH788cdgmS+++MJYlmW2bNliX/AA6sT5559vrr766pBpAwcONEOGDDHG0A8ARzpJZsqUKcH3dXXOP//886ZBgwYh3wvuvPNO065dO5u3CHbgjgjYoqSkRIsXL1bv3r2D0xwOh3r37q358+fXY2QA7JKdnS1JSkhIkCQtXrxYXq83pB9o3769mjVrFuwH5s+frw4dOqhhw4bBMn379lVOTo5+/fXXQxg9gIMxYsQInX/++SHnu0Q/ABwtPvnkE5100kkaNGiQUlJS1LlzZ7388svB+enp6dq+fXtIXxAXF6dTTjklpC+Ij4/XSSedFCzTu3dvORwOLViw4NBtDIBaOf300zVz5kytWbNGkrR06VLNnTtX5557riT6AeBoU1fn/Pz583XWWWcpLCwsWKZv375avXq19u7de4i2BnXFVd8B4Mi0a9cu+f3+kB8VJKlhw4ZatWpVPUUFwC6BQECjRo3SGWecoRNOOEGStH37doWFhSk+Pj6kbMOGDbV9+/Zgmcr6ibJ5AA5/7777rn766Sf9+OOPFebRDwBHhw0bNuiFF17QP/7xD91999368ccfNXLkSIWFhWno0KHBc7myc718X5CSkhIy3+VyKSEhgb4A+BO46667lJOTo/bt28vpdMrv92vs2LEaMmSIJNEPAEeZujrnt2/frpYtW1ZYR9m88o+BxOGPRAQA4KCNGDFCy5cv19y5c+s7FACH0ObNm3XLLbdoxowZCg8Pr+9wANSTQCCgk046SQ8//LAkqXPnzlq+fLlefPFFDR06tJ6jA3AovP/++5o0aZLefvttHX/88VqyZIlGjRqlRo0a0Q8AACQxWDVskpSUJKfTqR07doRM37Fjh1JTU+spKgB2uOmmmzRt2jR9++23atKkSXB6amqqSkpKlJWVFVK+fD+QmppaaT9RNg/A4W3x4sXKzMxUly5d5HK55HK59N133+npp5+Wy+VSw4YN6QeAo0BaWpqOO+64kGnHHnusNm3aJOmPc3l/3w1SU1OVmZkZMt/n82nPnj30BcCfwO2336677rpLgwcPVocOHXTFFVfo1ltv1SOPPCKJfgA42tTVOc93hSMLiQjYIiwsTF27dtXMmTOD0wKBgGbOnKnTTjutHiMDUFeMMbrppps0ZcoUffPNNxVul+zatavcbndIP7B69Wpt2rQp2A+cdtpp+uWXX0IuPmbMmKHY2NgKP2gAOPz06tVLv/zyi5YsWRJ8nXTSSRoyZEjw//QDwJHvjDPO0OrVq0OmrVmzRs2bN5cktWzZUqmpqSF9QU5OjhYsWBDSF2RlZWnx4sXBMt98840CgYBOOeWUQ7AVAA5GQUGBHI7Qn5icTqcCgYAk+gHgaFNX5/xpp52m2bNny+v1BsvMmDFD7dq147FMf0b1PVo2jlzvvvuu8Xg8ZuLEiWbFihVm+PDhJj4+3mzfvr2+QwNQB2644QYTFxdnZs2aZbZt2xZ8FRQUBMtcf/31plmzZuabb74xixYtMqeddpo57bTTgvN9Pp854YQTTJ8+fcySJUvMl19+aZKTk83o0aPrY5MA1IEePXqYW265JfiefgA48i1cuNC4XC4zduxYs3btWjNp0iQTGRlp3nrrrWCZcePGmfj4ePPxxx+bZcuWmYsuusi0bNnSFBYWBsv069fPdO7c2SxYsMDMnTvXtG3b1lx++eX1sUkAamjo0KGmcePGZtq0aSY9Pd189NFHJikpydxxxx3BMvQDwJElNzfX/Pzzz+bnn382ksx//vMf8/PPP5uNGzcaY+rmnM/KyjINGzY0V1xxhVm+fLl59913TWRkpHnppZcO+fbi4JGIgK2eeeYZ06xZMxMWFma6detmfvjhh/oOCUAdkVTpa8KECcEyhYWF5sYbbzQNGjQwkZGR5uKLLzbbtm0LWU9GRoY599xzTUREhElKSjL//Oc/jdfrPcRbA6Cu7JuIoB8Ajg6ffvqpOeGEE4zH4zHt27c3//vf/0LmBwIBc99995mGDRsaj8djevXqZVavXh1SZvfu3ebyyy830dHRJjY21lx11VUmNzf3UG4GgFrKyckxt9xyi2nWrJkJDw83rVq1Mvfcc48pLi4OlqEfAI4s3377baW/CQwdOtQYU3fn/NKlS82ZZ55pPB6Pady4sRk3btyh2kTUMcsYY+rnXgwAAAAAAAAAAHCkY4wIAAAAAAAAAABgGxIRAAAAAAAAAADANiQiAAAAAAAAAACAbUhEAAAAAAAAAAAA25CIAAAAAAAAAAAAtiERAQAAAAAAAAAAbEMiAgAAAAAAAAAA2IZEBAAAAAAAAAAAsA2JCAAAAAAAAAAAYBsSEQAAAACOWhdffLEaNGigSy+9tL5DAQAAAI5YJCIAAAAAHLVuueUWvfHGG/UdBgAAAHBEIxEBAAAA4JDbvXu3UlJSlJGRUa9x9OzZUzExMRWmDx48WE8++WQ9RAQAAAAceUhEAAAAAAdh2LBhsiyrwmvdunX1HdphbezYsbrooovUokWL+g6lUvfee6/Gjh2r7Ozs+g4FAAAA+NNz1XcAAAAAwJ9dv379NGHChJBpycnJFcqVlJQoLCzsUIV12CooKNCrr76q6dOn215Xp06d5PP5Kkz/6quv1KhRoyqXO+GEE9S6dWu99dZbGjFihJ0hAgAAAEc87ogAAAAADpLH41FqamrIy+l0qmfPnrrppps0atQoJSUlqW/fvpKkQCCgRx55RC1btlRERIQ6duyoyZMnh6wzPz9fV155paKjo5WWlqYnn3xSPXv21KhRo4JlWrRoofHjx4cs16lTJ40ZM6ba9fTs2VMjR47UHXfcoYSEBKWmpgaXLxMIBPTYY4+pTZs28ng8atasmcaOHas33nhDiYmJKi4uDik/YMAAXXHFFVW21+effy6Px6NTTz31gHWUj/Pmm2/WqFGj1KBBAzVs2FAvv/yy8vPzddVVVykmJkZt2rTRF198EVLXkiVLtHz58gqv/SUhyvTv31/vvvvuAcsBAAAA2D8SEQAAAICNXn/9dYWFhWnevHl68cUXJUmPPPKI3njjDb344ov69ddfdeutt+pvf/ubvvvuu+Byt99+u7777jt9/PHH+uqrrzRr1iz99NNPNaq7OvWUxRgVFaUFCxboscce07/+9S/NmDEjOH/06NEaN26c7rvvPq1YsUJvv/22GjZsqEGDBsnv9+uTTz4Jls3MzNRnn32mq6++usq45syZo65du4ZMq6qOfeNMSkrSwoULdfPNN+uGG27QoEGDdPrpp+unn35Snz59dMUVV6igoKBG7VSVbt26aeHChRUSLQAAAABqxjLGmPoOAgAAAPizGjZsmN566y2Fh4cHp5177rn64IMP1LNnT+Xk5IQkEIqLi5WQkKCvv/5ap512WnD6tddeq4KCAr399tvKy8tTYmKi3nrrLQ0aNEiStGfPHjVp0kTDhw8P3gXRokULjRo1KuQuiU6dOmnAgAEaPXr0AeuRSu808Pv9mjNnTrBMt27ddPbZZ2vcuHHKzc1VcnKynn32WV177bUVtv/GG29URkaGPv/8c0nSf/7zHz333HNat26dLMuqtM0GDBigxMREvfrqq5J0wDoqi9Pv9ysuLk4DBw7UG2+8IUnavn270tLSNH/+/JC7Lfand+/eWrp0qfLz85WQkKAPPvgg2F7Lli1Tx44dlZGRoebNm1drfQAAAAAqYowIAAAA4CD95S9/0QsvvBB8HxUVFfz/vn/5v27dOhUUFOicc84JmV5SUqLOnTtLktavX6+SkhKdcsopwfkJCQlq165dtWOqTj1lTjzxxJD3aWlpyszMlCStXLlSxcXF6tWrV6X1/P3vf9fJJ5+sLVu2qHHjxpo4cWJwAO+qFBYWhiRuDlRHZXE6nU4lJiaqQ4cOwWlld1CUxV4dX3/9dZXzIiIiJKnO7rAAAAAAjlYkIgAAAICDFBUVpTZt2lQ5r7y8vDxJ0meffabGjRuHzPN4PDWq1+FwaN8bnL1eb43rcbvdIe8ty1IgEJD0x4/xVencubM6duyoN954Q3369NGvv/6qzz77bL/LJCUlae/evcH3B6pjf3GWn1aW/CiL/WDt2bNHUuUDjwMAAACoPsaIAAAAAA6h4447Th6PR5s2bVKbNm1CXk2bNpUktW7dWm63WwsWLAgut3fvXq1ZsyZkXcnJydq2bVvwfU5OjtLT06tdT3W0bdtWERERmjlzZpVlrr32Wk2cOFETJkxQ7969D7j+zp07a8WKFTWqoz4sX75cTZo0UVJSUn2HAgAAAPypcUcEAAAAcAjFxMTotttu06233qpAIKAzzzxT2dnZmjdvnmJjYzV06FBFR0frmmuu0e23367ExESlpKTonnvukcMR+ndEZ599tiZOnKj+/fsrPj5e999/v5xOZ7XrqY7w8HDdeeeduuOOOxQWFqYzzjhDO3fu1K+//qprrrlGkvR///d/uu222/Tyyy8Hx2vYn759+2r06NHau3evGjRoUK066sOcOXPUp0+feqsfAAAAOFKQiAAAAAAOsYceekjJycl65JFHtGHDBsXHx6tLly66++67g2Uef/xx5eXlqX///oqJidE///lPZWdnh6xn9OjRSk9P1wUXXKC4uDg99NBDwTsiqltPddx3331yuVy6//77tXXrVqWlpen6668Pzo+Li9Mll1yizz77TAMGDDjg+jp06KAuXbro/fff13XXXVetOg61oqIiTZ06VV9++WW9xQAAAAAcKSyz70NlAQAAAByWevbsqU6dOmn8+PH1HUoFvXr10vHHH6+nn366WuU/++wz3X777Vq+fHmFOz0OBy+88IKmTJmir776qr5DAQAAAP70uCMCAAAAQK3t3btXs2bN0qxZs/T8889Xe7nzzz9fa9eu1ZYtW2o0ZsWh4na79cwzz9R3GAAAAMARgUQEAAAAgFrr3Lmz9u7dq0cffVTt2rWr0bKjRo2yJ6g6cO2119Z3CAAAAMARg0czAQAAAAAAAAAA2xx+D2MFAAAAAAAAAABHDBIRAAAAAAAAAADANiQiAAAAAAAAAACAbUhEAAAAAAAAAAAA25CIAAAAAAAAAAAAtiERAQAAAAAAAAAAbEMiAgAAAAAAAAAA2IZEBAAAAAAAAAAAsA2JCAAAAAAAAAAAYBsSEQAAAAAAAAAAwDYkIgAAAAAAAAAAgG1IRAAAAAAAAAAAANuQiAAAAAAAAAAAALb5f+UyM7XzU6joAAAAAElFTkSuQmCC"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Common y-value for all frequencies\n",
    "y_value = 100\n",
    "\n",
    "plt.figure(figsize=(20, 4))\n",
    "\n",
    "# Scatter plot for each set of frequencies\n",
    "plt.scatter(e_freqs, [y_value] * len(e_freqs), label='E Frequencies', alpha=0.6, marker='s')\n",
    "plt.scatter(a_freqs, [y_value] * len(a_freqs), label='A Frequencies', alpha=0.6, marker='o')\n",
    "plt.scatter(es_freqs, [y_value] * len(es_freqs), label='Es Frequencies', alpha=0.6, marker='^')\n",
    "\n",
    "plt.xlabel('Frequency (cm$^{-1}$)')\n",
    "plt.yticks([])  # Hide y-axis ticks as they're irrelevant in this context\n",
    "plt.title('Visualization of Frequencies Along the X-axis')\n",
    "plt.legend()\n",
    "plt.grid(True, axis='x')\n",
    "\n",
    "plt.show()"
   ],
   "metadata": {
    "collapsed": false,
    "ExecuteTime": {
     "end_time": "2024-04-05T02:50:19.826827Z",
     "start_time": "2024-04-05T02:50:19.662062Z"
    }
   },
   "id": "6fbb3b81a0bdd4a0",
   "execution_count": 40
  },
  {
   "cell_type": "code",
   "outputs": [],
   "source": [],
   "metadata": {
    "collapsed": false,
    "ExecuteTime": {
     "end_time": "2024-04-05T02:48:32.474498Z",
     "start_time": "2024-04-05T02:48:32.473288Z"
    }
   },
   "id": "dda613550f2b543f",
   "execution_count": 39
  },
  {
   "cell_type": "code",
   "outputs": [],
   "source": [],
   "metadata": {
    "collapsed": false,
    "ExecuteTime": {
     "end_time": "2024-04-05T02:48:32.476816Z",
     "start_time": "2024-04-05T02:48:32.475627Z"
    }
   },
   "id": "d3dc075ae0816dbd",
   "execution_count": 39
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}