{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "51a8a104-4b79-4c02-9337-942edabafc8f",
   "metadata": {},
   "source": [
    "---\n",
    "**© 2026 MaLab @ Purdue University. All rights reserved.**\n",
    "\n",
    "**Authors:** Botao Du, Qihao Guo, Ruichao Ma\n",
    "\n",
    "**Group:** [MaLab @ Purdue](https://www.ma-quantumlab.com/)\n",
    "\n",
    "This notebook and its contents are the intellectual property of the authors.\n",
    "\n",
    "---"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d1fe891c-38cb-4ca8-9e12-12cb437aabd4",
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "from numpy import *\n",
    "from matplotlib.pyplot import *\n",
    "\n",
    "import math as math\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "51d0ee8c-044b-4a58-8349-be1793fa35e9",
   "metadata": {},
   "source": [
    "# Fig.3, Load Data"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "773e9d86-baa1-4780-800f-5659cdeb8536",
   "metadata": {},
   "source": [
    "Load decay rate from Fig_3_data.npz\n",
    "- theta_list: relative phases between the sideband to Q1 and the sideband to Q2, unit ($\\times 2\\pi$), shape: (index_theta,) = (26,).\n",
    "- Gamma_Engineered_all_removed: Engineered T1 decay rates from NQ single-excitation plus eigenstate, unit (2$\\pi\\times$ KHz). shape: (index_state, index_theta,) = (4, 26,).\n",
    "- Gamma_Intrinsic_all: Intrinsic T1 decay rates from NQ single-excitation minus eigenstate, unit (2$\\pi\\times$ KHz).shape: (index_theta,) = (4,)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "1168df80-ff8f-45b5-9326-daf6351ed28a",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Dataset for Fig.3, Programmable Superradiance in an Interacting Qubit Array, MaLab, arXiv: 2605.12442\n",
      "['theta_list', 'Gamma_Engineered_all_removed', 'Gamma_Intrinsic_all', 'note']\n"
     ]
    }
   ],
   "source": [
    "infile = 'data/Fig_3_data.npz'\n",
    "exp_data = np.load(infile)\n",
    "print(exp_data['note'])\n",
    "print(exp_data.files)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "70d8f32e-daf3-4bbf-aed2-172dca04c044",
   "metadata": {},
   "outputs": [],
   "source": [
    "theta_list = exp_data['theta_list']\n",
    "Gamma_Engineered_all_removed = exp_data['Gamma_Engineered_all_removed']\n",
    "Gamma_Intrinsic_all = exp_data['Gamma_Intrinsic_all']"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "96e6d6fa-b5bb-4e76-883a-d75d7d303f39",
   "metadata": {},
   "source": [
    "# Plot Fig.3"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "5399476c-f582-47bc-8a1e-fd225d3896f9",
   "metadata": {},
   "outputs": [],
   "source": [
    "order_2q = [1, 2]\n",
    "order_3q = [1, 2, 0]\n",
    "order_4q = [3, 1, 2, 0]\n",
    "\n",
    "# ---------- in phase ----------\n",
    "temp = Gamma_Engineered_all_removed[:,:,0]\n",
    "eng_in_phase_1q = [15]   # 15 kHz calibrated singel qubit engineered decay\n",
    "eng_in_phase_2q = temp[0][order_2q]\n",
    "eng_in_phase_3q = temp[1][order_3q]\n",
    "eng_in_phase_4q = temp[2][order_4q]\n",
    "\n",
    "eng_in_phase_data_list = [\n",
    "    eng_in_phase_1q, eng_in_phase_2q, eng_in_phase_3q, eng_in_phase_4q]\n",
    "\n",
    "eng_in_phase_k = [\n",
    "    [d[k] for d in eng_in_phase_data_list if len(d) > k]\n",
    "    for k in range(max(len(d) for d in eng_in_phase_data_list))\n",
    "]\n",
    "eng_in_phase_k1, eng_in_phase_k2, eng_in_phase_k3, eng_in_phase_k4 = eng_in_phase_k"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "c15f54d3-5620-4a30-b195-13c156a51d27",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmkAAAHbCAYAAACQmw0xAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuNSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/xnp5ZAAAACXBIWXMAACZzAAAmcwHzbHUKAAB9uElEQVR4nO3dd1xV9f8H8Ndhb0RwIqICThRXThTce+e23FpZVlY2HJmaad8sM03L3LnLvfcgtwg4E0VQEQd7IxfO7w9+98SVdYF77zlXX8/Hg4dnfD6f877X0ref8xmCKIoiFMjExASCIKBLly7Yv3+/3OEQERERGZSZIR/Wvn37Yte5ePFiseoJgoBjx44V+zlERERESiIYsidN3TumL6IoQhAEZGVl6e0ZRERERIZg0J40ICeRIiIiIqLCGTRJGzlypCEfR0RERGS0DPq6k4iIiIi0YyJ3AERERESUF5M0IiIiIgVikkZERESkQIpI0u7fv6+ztr799ludtUVEREQkF0UkaZ06dcLTp09L3c4333yDmTNn6iAiIiIiInkpIkkLCwtD586dER8fX+I2vv76a3zzzTe6C4qIiIhIRopI0gDg+vXr6NGjB1JTU4tdd/r06Zg7d64eoiIiIiKShyKSNE9PTwDA+fPn0b9/f6hUKq3rfvnll/juu++knQzeffddvcRIREREZEiKSNKOHDmCypUrS8fDhw/XavuoqVOn4vvvv5fKvv/++1iyZIleYyUiIiIyBEUkae7u7jh06BDKli0LAPjrr7/wzjvvFFrnk08+wcKFC6UE7cMPP8TixYv1HisRERGRISgiSQOAunXrYt++fbC1tQUA/PHHH/jiiy/yLfvxxx9j0aJFUoI2ZcoU/PTTTwaLlYiIiEjfFJOkAUCzZs2wY8cOWFhYAAD+97//4fvvv9coM3nyZCxevFhK0D777DP88MMPBo+ViIiISJ8UucH633//jSFDhiArKwuCIOC3337DuHHj8N577+G3336TErQvvvgC8+bNkzlaIiIiIt1TZJIGACtXrsT48eMhCAJMTEzQrl07HD16VLo/bdo0zJkzR8YIiYiIiPRHsUkakPO68/PPP4cgCAAg9aBNnz4ds2fPljM0IiIiIr1S1Ji0l3322WeYOnUqRFGEKIoQBAGzZs1igkZERESvPIP1pJ0+fbrEdT/77DNcvnwZ3bp1w9SpU4ss37Zt2xI/i4iIiEgJDJakmZiYSK8t9UkQhGLtWEBERESkRGaGfmBJc0JBEEpcl4iIiMjYGCxJq1q1qkF60oiIiIheBYqe3UlERET0ulL07E4iUhZ/f39pljXlGDt2LARBwJ49ewotd+jQIfTu3RsVKlSAjY0NGjRogKVLlxooSv1r27YtTE1NcePGDblDIXplMEkjksmoUaMgCAL8/f3lDoVK6MaNG1i7di0aNWqEXr165VsmKysLY8eORdeuXbFnzx4kJycjKysL165dw/vvv4+PP/7YwFEX7fnz5/jkk0/g5eUFa2truLi4oHPnzti5c2eBdWbMmIHs7OwC91wmouJjkkZEVEJfffUVsrKyMG3atALLjBo1CqtWrUKtWrVw4sQJpKSkICkpCV999RUA4Oeff0ZQUJCBIi7ajRs34O3tjR9//BF3796Fubk54uPjceTIEfTr1w8ffvhhvvU6deqE5s2bY+/evThz5oyBoyZ6NTFJIyIqgdDQUOzevRvly5dH37598y2zcuVK/Pnnn6hWrRpOnz4t9ZpaWFjg22+/RePGjSGKIjZs2GC4wAuRkZGB3r1749mzZ/D29kZQUBASExORmJiIuXPnQhAELF68GKtXr863/vjx4wEAP/74oyHDJnplMUkjIiqBlStXAgAGDRoEU1PTPPeTkpLw+eefAwCWLl2K8uXL5ynToUMHAMDly5f1GKn2fv/9d4SFhcHGxgb79u2Dj48PAMDGxgbTpk3De++9ByBna77MzMw89fv37w8LCwvs3bsXT58+NWjsRK8iJmlECpR7gH5WVhYWLVoEHx8f2NjYoGzZsujZs2eJ/2LP3XZ6ejq+/vpr1K5dG9bW1ihfvjyGDh2KO3fuFNlOceOKjo7Gr7/+ij59+qB27dqwt7eHra0t6tatiylTpuDx48cFPiszMxO//PILWrVqhTJlysDc3BwVKlRAgwYN8O677xb6ei00NBTvvvsuatasCRsbG9jb26Nhw4b45ptvkJCQUPQXlo/s7GysXbsWADB48OB8yyxfvhwxMTFo1KgRunfvnm+ZsmXLAgAiIyNLFIeu/fnnnwCAoUOHomrVqnnuT506FYIg4PHjxzhx4kSe+05OTujUqRNUKhXWrVun93iJXnkiEcli5MiRIgDRz88vzz0/Pz8RgDht2jSxS5cuIgDR3NxctLOzEwGIAEQrKyvx7NmzxX6uuu0vvvhCbNGihQhAtLCwEB0cHKS2bWxsxFOnTuk0rk8++UQqY2ZmJpYtW1Y0MzOTrpUrV04MDg7OUy8zM1Ns166dVE4QBLFMmTIadfv06ZPvZ/3jjz9Ec3Nzjc+V+7x69eri3bt3i/0dXr16VfreMjIy8i1Tp04dEYC4ZMmSAtuZM2eOCECsVatWsWPQtaSkJFEQBBGAuG3btgLL1a1bVwQgTp06Nd/78+bNEwGIHTt21FeoRK8N9qQRKdjSpUtx6dIlbNmyBcnJyUhKSkJwcDC8vb2Rnp5e4CBubSxbtgzXr1/H+vXrkZycjISEBFy5cgUNGjRAamoqBg0ahLi4OJ3FVbVqVcybNw8hISFIS0tDTEwM0tPTceXKFXTr1g3Pnz/HsGHD8uwssmnTJpw4cQI2NjZYv349UlNTERcXh/T0dDx48ADLli1D48aN8zxv//79GD9+PMzNzTFnzhxERUUhJSUFqamp+Oeff9C0aVPcv38f/fv3R3Z2drG+O3XPXf369WFhYZHn/o0bN3Dr1i0IgoB+/foV2I76+7W3ty/W8/Xh1q1b0nfv7e1dYDn1vZs3b+Z7v1mzZgCAc+fOcYs+otKSO0skel1p05MGQDxz5kye+5cvX5buR0REFOu5udvetGlTnvtRUVGio6OjCECcM2eOQeLKyMiQemhOnjypce/dd98VAYjvvPOO1u2pVCqxRo0aIgBxx44d+ZaJiYkRK1WqJAIQ//77b63bFkVRHD58uAhAHD16dL73lyxZIvX6VahQocAfa2trEYA4ePDgIp+5evVq6bst7k9+/429bNeuXVL5hISEAst99NFHIgCxSZMm+d6Pjo6W2gkKCiryuURUMPakESlYmzZt4Ovrm+d6kyZNUKVKFQDA9evXS9R2tWrVMGTIkDzXK1asiLFjxwIA/vrrL4PEZWFhgU6dOgEA/vnnH417Dg4OAICoqCit2zt16hTCwsJQs2bNAmdeli1bFt26dQMAHD58WOu2c8dSrly5fO+fPXsWQM5exU+fPi3wJy0tDQBQu3btIp9pbW2NChUqlOhHPfatMMnJydKxjY1NgeXU95KSkvK9X7ZsWZiY5PzVUpzfMyLKy+AbrBOR9t54440C77m6uuLRo0cFvpIsip+fX6H3fvzxR1y/fh0vXrzI80qvpHHdvn0bS5YswenTpxEeHo7k5OQ8rzdfnkDQrVs3LFiwALt27ULv3r0xcuRI+Pn5wcXFpcAY1ElSREQEKlasWGA5dWLy4MGDAsvkJzo6GkDOQPn83Lp1C0DOUhQFLVYbFRWFypUrAwBatWpV5DMHDx5c4CQFJREEAWXKlEFsbKz0PRFRyTBJI1KwwsYqWVlZAUC+SyFoQ50g5MfV1RVAzgzOuLg4VKhQodRxbd68GW+//bZ03cTEBI6OjrC0tASQkzClpKQgJSVFo56fnx/mzJmDb775Bnv27JG2X6pZsya6d++OiRMn5umJUvfgZGRkaLUURGpqapFlcktPTwcAKfaXhYWFAYDUq5if48ePA8jpRdQmSdM3Ozs76Tg1NVXqwXyZ+rvS5r8BdU8hEZWMIl53RkZGFnvgbm6PHj3K869xIlKO58+fY/z48cjMzMTgwYNx+fJlpKenIy4uDk+ePMGTJ0+kHqf8/l+ePn067t69i/nz56Nr165wcHDAnTt3sGjRInh7e+P333/XKJ+VlQUA6NGjB0RRLPLn5MmTxfo86teH8fHx+d5XvwosKNEBgB07dgAAunfvrpEgySV30l7Ycijqe5UqVSqwTGxsLADA2dlZR9ERvZ5kT9Ju376N5s2bY+DAgSWaCXTjxg288cYbGDp0qPQHMxEVTZu/iE1NTQt8pVccBw4cQHJyMurWrYuNGzeiSZMmMDc31yhTVI+Xu7s7Pv/8cxw4cACxsbE4efIk/P39kZWVhcmTJ+Phw4dSWXXPX3FfY2pL/ao1v1e62dnZUqIpCEK+9Z8/f459+/YBAEaPHq3VM7ds2YKKFSuW6Kd///5Ftl+7dm0p3sI2SVffq1u3br7309PTpZ7Gwl5JE1HRZE3SYmNj4efnh6ioKOzcuRP9+vXDixcvtK4fHByM9u3b4+nTp9i2bRsmTpyox2iJXi2FLQB76tQpADnLLeS3xERxPXr0CADQoEEDaVB5bqIoSq//tGFqago/Pz/s3bsXlpaWyMjIwLlz56T76teHN27cQHh4eOmCz0edOnUAAPfv389zz8TEREoSCxo4/8MPPyA9PR0+Pj4Fbsz+srS0tEInIRT2o+7ZKoydnZ20fMbBgwfzLfPo0SNp6Q31bgkvy/2daDMhgogKJmuSVrZsWcycOVM6379/P3r16iX9K6wwly9fRocOHfD8+XMAOV31X3zxhd5iJXrVhIWFYdu2bXmuP3v2DKtWrQIADBw4UCfPcnR0BJAz4zO/15krVqzAvXv38q1b2D/cLCwspC2ZMjIypOvt27eHu7s7srOz8fHHHxc6nCIzM1NjZqM21DNbC9pdoWnTpgAg9Zbl9s8//2DRokUwMTHBL7/8UmBv28tGjRql1avb0rzOHT58OICctely90yqff/99xBFEZUrV0a7du3ybePixYsAAC8vr3y3wiIi7cn+unPSpElYvny5dH706FF069Ytz+Dh3M6fP49OnTpJ/zpUb17s6emp93iJXhWOjo4YO3YsNm7cKA01CAoKQpcuXRAXF4fy5ctLezWWVocOHSAIAq5fv47JkydLY7kSExPxv//9D5MmTSpw/NLbb7+N0aNH49ChQ0hMTJSuR0REYOTIkUhNTYWFhYVG0mBubo5ff/0VJiYm2LlzJ7p164bz589LyVp2djZu3ryJ+fPno2bNmggKCirW52ndujVMTEzw5MkTqZcwtxEjRgDIWcJk2bJlyM7Ohkqlwvr169G9e3e8ePECs2fPRps2bYr1XH2bMGECatSogZSUFPTs2RMhISEAcnrx5s+fjyVLlgAA5s6dm+d1tdqlS5cAAG3btjVM0ESvMgOsxaaVNWvWiKampqKJiYloYmIitm7dWkxMTMxT7vTp06KDg4MoCIIoCILo6elZ7MU8iZRAm8Vsv/766wLrq8usXr26WM/NvS1U8+bNRQCipaVlsbaFKklckydP1lhgtUyZMqKJiYkIQOzSpYs4bdo0EYA4cuRIjXp9+vTJsyWUra2tdM3ExET8/fff841l48aN0oKx6s/p7OyssTUUADEgIKBY36EoimKnTp1EAOLy5cvz3MvOzpa2zQIgWltbixYWFiIA0dTUVJw9e3axn2co169fF8uXLy/F7uDgIJqamkrnH3zwQaH1q1WrJgIQDx8+bKCIiV5dsvekqY0cORLr16+XXl2cO3cOHTp00BiYe/z4cXTv3l2aOVWrVi2cOnUq342AiahwlpaWOHnyJGbOnAl3d3e8ePEC5cqVw5AhQxAYGKjznpCff/4Zy5cvR8OGDWFpaYmsrCw0atQIixYtwr59+2Bmlv+KQPPnz8f333+Pbt26wcPDA5mZmcjMzET16tUxcuRIXLp0CePHj8+37tChQxEaGoovvvhCem58fDzs7e3RvHlzTJkyBQEBAWjdunWxP4/6mRs3bsxzTxAE7Nq1CzNnzoSHhweys7Ph7OyMoUOH4vz585gxY0axn2co9erVw7Vr1/Dxxx/D09MTGRkZcHR0RMeOHbFjxw4sXry4wLrnz59HeHg4qlWrho4dOxowaqJXkyCKylq74u+//8awYcOgUqkgiiJ8fHxw5MgRXL58GQMGDJDW3alXrx6OHj2aZ/0mIiqcv78/Tp06ha+//hqzZs2SOxyj9eLFC7i7u+Pp06cIDw/nPxYBTJ48Gb/88gvmzp2LadOmyR0OkdFTTE+a2oABA/DXX3/BwsICgiAgJCQELVq0QL9+/aQEzcfHBydOnGCCRkSysbCwwPTp0yGKIhYuXCh3OLKLiYnB6tWrUa5cOUyePFnucIheCYpL0gCgV69e2LFjh7Sad1hYGDIyMiAIApo0aYLjx49z/R0ikt2ECRPg4eGBFStWaLWzwats0aJFSE5OxvTp0wvdjYCItKfIJA0Aunbtig8++ACiKEpT1E1MTLBmzRqdLK5JRFRa5ubmWL16NaZOnYqIiAi5w5GVs7MzZs+ejXfeeUfuUIheGYrdu3Pt2rX48ccfIQiClKhlZWWhb9++OH78ONzc3OQOkYgIbdq0UdxSGnL46KOP5A6B6JWjuIkDQM7Clu+++y6ys7MhCAIaN26M0NBQJCUlQRRFVKtWDcePH0e1atXkDpWIiIhILxT3unPp0qVSggbkzEQ7efIkjhw5gjJlykAQBISHh6Nt27a4e/euzNESERER6YeikrSffvoJkydPlhK0Ll26YN++fbC1tcUbb7yBo0ePwtnZGYIg4NGjR/D398ft27dljpqIiIhI9xSTpM2fPx+ffvqptK9f7969sWvXLlhZWUllGjVqhOPHj6N8+fIQBAGPHz+Gv78/bty4IVfYRERERHqhiCRt9uzZmDZtmpSgDRw4UFor7WXe3t44efIkKlasCEEQ8OzZM/j7+yM4ONjQYRMRERHpjexJ2vTp0/HNN99ICdqIESOwadOmAreIAXK2gzp9+jSqVKkCQRAQExOD9u3bSxv7EhERERk7WZO02NhYbNq0SUrQxo4di7Vr18LEpOiwPDw8cPr0aWmGZ1xcXL576BEREREZI9mX4IiMjETnzp3Rrl07LFmypET1O3TogG7duuGnn37SQ4REREREhid7kgYASUlJpdpGpLT1iYiIiJRGEUkaEREREWmSfeIAEREREeXFJI2IiIhIgZikERERESkQkzQiIiIiBSp4xVgiLdja2iI9PR2mpqYoX7683OEQERmlZ8+eISsrC1ZWVkhJSZE7HFIIzu6kUjE1NUV2drbcYRARvRJMTEyQlZUldxikEOxJo1JRJ2kmJiaoVKmS3OEQERmlqKgoZGdnw9TUVO5QSEGYpFGplC9fHpGRkahUqRIePXokdzhEREapSpUqiIyM5LAR0sCJA0REREQKpIgk7cSJE0hNTZU7DCIiIiLFUESS1qFDB5QpUwbvvPOO3KEQERERKYIikjQAyMrKwoMHD+QOg4iIiEgRFJOkATmzW86dO4fk5ORi13327BmWLl2Kc+fOgauKEBERkbFT1OzOkJAQ+Pr6AgDc3d1Rv359NGjQAA0aNED9+vVRq1YtCIKQb9309HR88MEHMDExwalTp9C6dWtDhk5ERESkU4pK0nL3gEVERCAiIgJ79+6VrllaWqJu3boayVvdunXh5OSEbdu2SW0EBQUxSSMiIiKjpqgkzdbWFoIgIDk5Od9Xlunp6bh69SquXr1aaDuJiYn6CpGIiIjIIBSVpLVp0wb79+9HWFgYQkJCEBISguDgYISEhCAsLAyiKBY43iz3a9AqVaoYKmQiIiIivVBUkqZWo0YN1KhRA3379pWupaam4tq1a1LSFhISgmvXriEhIQHAf69Kvby80K9fPznCJiIiItIZRSZp+bGxsUHz5s3RvHlzjeuPHz/G/fv3kZCQgDJlyuCNN96Aubm5TFESERER6YbRJGkFqVy5MipXrix3GEREREQ6pah10oiIiIgoB5M0IiIiIgVSVJKWmpqKjIwMucMgIiIikp2ixqSdOXMGdnZ28PDwkHYZUP9ao0YNucMjIiIiMhhFJWmiKCIrKwuhoaEIDQ3F33//Ld2ztbWFt7e3RuLWoEEDlClTRr6AiYiIiPREUUmaWn4L1iYnJ+PChQu4cOGCxnVXV9c8iZu3t7ehQiUiIiLSC0UlaZ07d8bPP/8sLVSr/jU8PLzA3QYiIyMRGRmJgwcPAsjZeUClUhk6dCIiIiKdUkSS5ujoiISEBJiYmKBWrVqoVasWBg4cKN1PTk7GtWvXNBK3a9euIT4+vsBtooiIiIiMmSKStLi4ONy6dQtPnz7N976dnR1atmyJli1balx/+PChRuIWEhKCf//91xAhExEREemVIpI0AKhTpw7q1KlTrDpubm5wc3ND9+7dpWuZmZm6Do2IiIjI4BS1TpoucN9OIiIiehW8ckkaERER0auASRoRERGRAjFJIyIiIlIgWScOZGVl4fTp07h9+zYSExNRvnx5tGjRotgTCADgp59+wrVr1yAIAlauXKmHaImIiIgMR7Yk7bfffsPMmTMRHR2d556Pjw/mzp2rMWuzKIcPH8ahQ4eYpBEREdErQZbXnRMmTMB7772H58+f51mMVhRFBAUFoVevXhg+fDiSkpLkCJGIiIhIVgZP0n7++Wf88ccfEEURgiAAgLTlkzphEwQBoihi8+bNaNy4Me7du2foMImIiIhkZdAkLT4+HtOmTYMgCBAEARUqVMBvv/2Gx48fIz09Hf/++y/mzZuHChUqSHXu3buHNm3a4Pr164YMlYiIiEhWBk3S1q5di9TUVACAq6srLl68iPHjx6NixYqwsLCAl5cXvvjiC9y+fRtjxoyRetuePHkCf39/XL582ZDhEhEREcnGoEnasWPHpOOlS5eiSpUq+ZZzcHDAH3/8gT/++AMWFhYQBAGxsbHo0KEDAgICDBUuERERkWwMmqQFBwcDAMqWLYtevXoVWX7MmDE4fPgwHBwcIAgCkpKS0LVrV41kj4iIiOhVZNAkLTY2FoIgFGsdtDZt2uDkyZMoV64cBEFAamoqevXqhX379ukxUiIiIiJ5GTRJy8zMBACYmpoWq56Pjw9OnTqFypUrQxAEpKenY8CAAdi+fbs+wiQiIiKSnUGTNCcnJwDA06dPi123Vq1aOH36NNzd3SEIAl68eIEhQ4Zg48aNug6TiIiISHYGTdJq1qwJURRx9+5dJCcnF7t+9erVcfr0aXh5eUEQBKhUKowcORKrV6/WQ7RERERE8jFoktawYUMAOXt2njhxokRtVKlSBadOnUK9evWktsaPH49Lly7pKkwiIiIi2Rk0SWvXrp10vHbt2hK3U6FCBZw6dQpNmjQBAGRnZyMuLq7U8REREREphcGTNAsLCwDAnj17EBkZWeK2nJyccPz4cbRu3VpX4REREREphkGTNEdHR3Tq1AmiKCIzMxOzZ88uVXv29vY4fPgwOnTokGejdiIiIiJjZmboB/7www8YNWoUAMDc3LzU7VlbW2Pfvn34+eefpS2niIjo9SaKIv65G4O9IY8RnZwBFztL9GxQGa09nSEIgtzhEWlFENkFRaVQpUoVREZGwtXVFY8ePZI7HCIiJKRm4p0/r+BcWEyee608nLFseBM42pS+k0CX+Gcp5cegrzuJiIj0SRTFAhM0ADh7LwbvbrjCITJkFJikERHRK+OfuzEFJmhqZ+/F4Oy9wssQKQGTNCIiemXsDXmsVbk9wdqVI5ITkzQiInplRCdn6LQckZwMlqTNnj0bKSkphnqcJCUlpdRLfRARkXFwsbPUaTkiORksSZs1axY8PDwwf/58xMfH6/158fHx+O6771CjRg188803en8eERHJr2eDylqV6+WjXTkiORn0defz588xbdo0VK1aFe+++65e9tu8ePEiJk6ciKpVq2L69Ol4/vy5zp9BRETK1NrTGa08nAst08qj6DJESmCwJO3MmTNo2LAhRFFEcnIyfv/9d7Ro0QI1a9bE1KlTcfz4caSlpRW73dTUVBw5cgSffvopPD090bJlS/zxxx9ITk6GKIpo1KgRzpw5o4dPRERESiMIApYNb1JgEqZeJ40L2pIxMOhitqIoYuPGjZgzZw7u3LmTE0Cu/1HMzMxQt25d1KtXDx4eHqhcuTKcnJxgbW0NURSRnp6OuLg4REZG4u7du7h+/Tpu376NrKwsjWcAgJeXF2bOnIlhw4bxf0Y94gKMRKREoiji7L0Y7An+b8eBXj6V0cpDmTsO8M9Syo8sOw5kZ2dj27ZtWLx4Mc6dO6cZUDH/53k5/ObNm+Ojjz7CwIEDYWLCyav6xj9YiIhKj3+WUn4MvncnAJiYmGDw4MEYPHgwbty4gQ0bNuDvv/9GaGhovqtAqxO3gvJJDw8P9O/fHyNGjED9+vX1GjsRERGRIShq787w8HCcOnUKV65cwc2bNxEREYHo6Ghp6Q5bW1u4uLjA3d0dderUQZMmTdC2bVvUqFFD5shfX/zXHxFR6fHPUsqPLD1pBalWrRqqVauGkSNHyh0KERERkaw4aIuIiIhIgZikERERESmQol53EhERkfZEUcSLFy80lqIi/TE1NYWFhYXBlnFhkkZERGRksrKy8Pz5cyQmJjJBMzBTU1M4ODigXLlyMDU11euzmKQREREZkaysLDx8+LBEu/RQ6WVlZSEuLg7p6elwc3PTa6LGMWlERERG5Pnz50zQFCAtLU3v+4OzJ42IiMhIiKKIxMRE6VwQBJQvXx5WVlYyRvX6SE9Px7Nnz6TF9ZOSklChQgW9jVFjkkZERGQkXp4kUL58eZQtW1bGiF4vNjY2AICnT58CAFQqFV68eAFLS0u9PI+vO4mIiIzEy5ME2INmeC8nZPqcuMEkjYiIiEhLhlp+A2CSRkRERKRITNKIiIiIFIhJGhEREZECKSJJ++uvv6BSqeQOg4iIiEgxFJGkDRo0CK6urvjss89w+/ZtucMhIiIiI3TkyBEIgoAJEybIHYpOKCJJA4Do6Gj8+OOPqFevHnx9fbF27VquqExERERau3r1KgCgSZMmMkeiG4pI0hwcHCCKovRz7tw5jBkzBpUqVcK7776Ly5cvyx0iERERKVxgYCAAJmk69eTJE6xduxZt27YFAClZS0xMxO+//47mzZujYcOGWLp0KeLj4+UNloiIiBTp6tWrMDc3R/369eUORScUkaRZWVnhrbfewsmTJ3Hnzh18/vnnqFSpEoD/ErZr165h8uTJqFy5slSWiIiI9EcURQSERuOLv0Mwbu0lfPF3CAJCo6W9K5UkOTkZoaGh8Pb2zrMrwJYtW2BnZ4eKFSsiICBApgiLTxFJWm6enp747rvv8ODBA+zatQu9e/eGmZmZlKylp6dj48aN6NChA7y8vDB//nw8efJE7rCJiIheKQmpmRi24gJGrLyAzZce4uitZ9h86SFGrLyA4X9cQEJqptwhaggKCoIoihqvOlUqFaZMmYIhQ4agXr16uHLlCnx9fWWMsngUl6SpmZqaolevXti5cycePHiAefPmwcvLC8B/vWthYWGYNm0aqlatir59+2Lv3r3Izs6WOXIiIiLjJooi3vnzCs6FxeR7/+y9GLy74YqietTU49EaN24MIGcT9A4dOuCnn37C6NGjcfr0abi6usoZYrEpNknLrWLFivjiiy/w77//4uTJkxgxYgSsra2lZE2lUmHPnj3o06cPqlatiunTp+P+/ftyh01ERGSU/rkbU2CCpnb2XgzO3iu8jCHlntl57tw5NG7cGOfOncOSJUuwatUqjVegly9fxttvvw1PT08IgoBZs2bJFHXhjCJJy61t27ZYt24doqKisHDhQlhZWUEQBClhe/z4Mb777jt4eXmhR48eOHXqlNwhExERGZW9IY+1KrcnWLtyhhAYGAgzMzOcOXMGfn5+UKlUOHbsGCZNmpSn7D///IPz58/D19cXjo6OMkSrHaNL0gAgODgY06ZNw9y5c5GRkQFAc1d6URSRnZ2NgwcPon379ujZsyeioqLkCpeIiMioRCdn6LScvmVkZODmzZtQqVT49NNP4e7ujitXrqBNmzb5lv/ggw9w584drFmzBmXKlDFssMVgNElaYmIili1bhqZNm6Jx48b49ddfER8fL/WgOTg44L333sNvv/0Gf39/AP+NXdu/fz9atWqF6OhoeT8EERGREXCxsyy6UDHK6du1a9egUqnQpUsXVKxYEQ8ePEBYWFiB5U1MjCP9UXyUp0+fxttvv41KlSrh/fffx9WrVzUWvm3evDlWrlyJx48fY8mSJRg/fjyOHz+Of//9F2PGjJF+Ix48eIC5c+fK/GmIiIiUr2eDylqV6+WjXTl9U49H69atG3bt2gUTExP0798f9+7dkzmy0lFkkvbkyRPMnz8fNWvWRLt27bBhwwakpaVJiZmjoyMmTZqE4OBgnDt3DqNHj4a1tbVGG15eXvjjjz+wdetW6dqePXsM/VGIiIiMTmtPZ7TycC60TCuPossYinpmp4+PD5o1a4aVK1ciJiYGvXr1QkJCgszRlZxikrTs7Gzs3r1bmqE5bdo03Lt3T0rMAKBFixZYtWoVHj9+jF9++UWrFYX79+8Pb29viKKIR48e6ftjEBERGT1BELBseJMCk7BWHs5YNryJxnhwOal70nx8fAAAw4YNw7Rp03Dr1i0MGjQIWVlZcoZXYmZyBwAAX375JdatWyctSpt73RVHR0e89dZbmDBhAry9vUvUfvXq1XH9+nWoVCqdxEtERPSqc7Qxx4ZxzXH2Xgz2BD9GdHIGXOws0cunMlp5OCsmQcvKykJISAjc3Nzg5OQkXZ8zZw5u3ryJHTt24MMPP8SSJUtkjLJkFJGkLViwQFpGQ61ly5aYMGECBg8eDCsrq1K1bywDBImIiJREEAS09nRBa08XuUMp0O3bt5GWlib1oqkJgoD169fD19cXS5cuRe3atfH+++/LFGXJKCJJA3J6z8qUKYMRI0Zg4sSJqFevns7anjt3Lj766COdtfeykydPol27dsWud+vWLdSuXVursmfPnsWaNWsQEBAgvbatUqUKfH19MWrUKLRq1arYzyciIjJ29erVK3DnA1tbW+lVqDFSRJLWqlUrTJgwAYMGDSp1r1l+dJnwGVpKSgomT56MVatW5bl369Yt3Lp1CytWrMCYMWOwePFi2NrayhAlERGR8Xj+/Lm02H1qaipu3ryJv/76C7a2tujWrZvM0f1HEUmaMe1IXxQrKyv4+flpVdbOzq7Q+1lZWejfvz8OHz4sXbO2tka9evVgZmaGmzdvIjExEQCwatUqREZGYt++fTA1NS35ByAiInrF3bhxAwMHDpTOt23bhm3btsHd3R3h4eHyBfYSRSRpr5IKFSrg4MGDOmlrxowZGgna+PHjMX/+fJQtWxZATi/b/PnzpfXfDh06hJkzZ+Lbb7/VyfOJiIheRf7+/oraHL4gHFGvUJGRkfjpp5+k87feegu///67lKABOe/a58yZg+nTp0vXfvrpJzx+rJy91IiIiKhkmKQp1OLFi5Geng4AsLGxwaJFiwosO2PGDLi5uQEA0tLS8PPPPxsiRCIiItIjxb/uTE9PR0JCgrSRujaqVq2qx4gMY/v27dLxoEGDNHrQXmZhYYHRo0dj9uzZAIAdO3ZgwYIFeo+RiIiI9EdxSVpKSgrWrVuHXbt24fLly4iLiytWfUEQjH7R2n///Rd3796Vzrt27VpknW7duklJWmhoKO7cuYOaNWvqLUYiIiLSL0Ulafv27cOYMWMQHR0tXTOGgX26FhwcrHHesmXLIus0btwYFhYWePHihdQGkzQiIiLjpZgxaXv37kXfvn0RHR0t7depTtAEQShw+4nC7skhPj4egwYNQrVq1WBtbQ17e3tUr14dffv2xZIlS6QlMwpz69Yt6djCwkIab1aYl8vlboOIiIiMjyKStJSUFIwaNUraALVFixY4duwYkpOT0aVLFylZy87ORmJiIm7duoVVq1ahbdu20r3x48cjPT1d9k1UExISsG3bNkRERCA9PR3JyckIDw/Hrl278MEHH6Bq1ar45ZdfCm0jIiJCOq5SpYrWSWjusXilWeflxx9/RJUqVbT6iYqKKvFziIiIqGCKeN25atUqxMbGQhAENG/eHCdOnIClpWW+Ze3s7FCrVi3UqlULo0aNwrZt2zB69Gj88ccfuHPnDo4dOyb7Xp3VqlWDq6srLC0tER0djZs3b0rj5BISEjB58mQEBQVh5cqV+dbP3dvm6Oio9XMdHByk46SkpBJGn/P8yMjIEtcnIiKi0lNEkpZ7wdYffvihwAQtPwMHDoSDgwO6deuG06dP45tvvsE333yjjzALZGJigo4dO2L06NHo2rVrnpmYycnJ+PPPPzFjxgxpvN2qVatQs2ZNfP7553naS0lJkY6Ls02WtbV1vm0Ul4ODA1xdXbUqGxUVhezs7BI/i4iIiPKniNedISEhAAAXF5dCNwovaBJBly5d0L9/f4iiiKVLlxp8dmfbtm1x5MgRDBs2LN+lMuzs7PDOO+8gMDAQ1apVk67Pnj0bT58+zVM+MzNTOjYz0z6Pzl1WPYGgJKZMmYJHjx5p9VOpUqUSP4eIiIgKpogkLSYmBoIg5DsbMXfikZaWVmAb/fr1AwDExcXh5MmTOo9RF9zc3LB582bpPDU1Nd9XnjY2NtKxekFbbeQuy43WiYiIjJsikjR1z1fu13Vq9vb20vGTJ08KbENXg+b1rXnz5vD395fOjxw5kqdM7o3XC0tMX5aamppvG0SkbKIoIiA0Gl/8HYJxay/hi79DEBAa/VouQURE/1HEmDQnJyc8e/Ys38Hu5cqVk47v3LmDGjVq5NtG7jFYuddZUyJ/f3+pt+/OnTt57ru4uEjHxZk9mTuJdXZ2LnmARGQwCamZeOfPKzgXFqNxffOlh2jl4Yxlw5vA0cZcpuiISE6K6Enz8vKCKIr59oA1aNBAOs6v10ntxIkT0nHuWY5KlHscV34JZa1ataTjmJgYjR6ywjx8+FA6rl27dikiJCJDEEUx3wRN7ey9GLy74Qp71IheU4pI0ho3bgwAePbsWZ5Xmh06dJDWCVu5ciXCwsLy1L9+/TqWLVsmnfv4+Ogx2tLLnXTlHn+mVqdOHY3zoKCgItuMjIzE8+fPC2yDiJTnn7sxBSZoamfvxeDsvcLLEFGOI0eOQBAETJgwQe5QdEIRSVqHDh2k4/3792vcq1atGjp27AhRFJGYmIhmzZph7ty5OHDgAA4cOIDp06fD19cXycnJEAQB1atX12obJTndvHlTOi5fvnye+82aNdNYhiQgIKDINs+cOSMdW1lZoVmzZqWMkoj0bW/IY63K7QnWrhzR6+7q1asAgCZNmsgciW4oIknr1KkT7O3tIYoi1qxZk+f+zz//DBsbGwiCgNjYWHz99dfo2bMnevbsie+++05a/FUQBPz888+yL2ZbmLS0NOzevVs6z2/JETs7O43EdcOGDUW2m7tMhw4dOLuTyAhEJ2fotBzR6y4wMBAAkzSdsrKywm+//YY5c+aga9eueWY01q5dG/v375cmEeTe21M9VsPa2hpr1qxBjx49DB5/ccyYMUNjbbS+ffvmW27UqFHScUhICPbs2VNgm4GBgThw4EC+dYlIuVzstFu4W9tyRK+7q1evwtzcHPXr15c7FJ1QxOxOABgyZEih99u2bYvQ0FCsXLkSR48exYMHD5CZmYlKlSrBz88PEyZMkGVh1cOHD+PQoUP4+OOPUaVKlQLLZWZmYsaMGVi4cKF0rXHjxujdu3e+5d988034+PggODgYADBx4kR4eXnlmRAQFRWFESNGSHuWNmzYEAMGDCjtxyIiA+jZoDI2X3pYZLlePpUNEA3Rf5JeJCE0LrTIcl5OXrC3sC+ynCEkJycjNDQUDRs2zLNz0ZYtWzB27FjY2dnhr7/+gq+vr0xRFo9ikjRt2Nvb46OPPsJHH30kdyiS1NRU/Pjjj1i0aBFat24NPz8/eHt7w8XFBRYWFoiOjsbFixexYcMGjdmXZcuWxcaNGwvcPF0QBPzxxx9o27Yt0tLSEBUVhebNm+Pdd99F27ZtYWZmhosXL2LJkiVSz5y1tTVWrFih9YbsRCSv1p7OaOXhXOjEgFYeOWWIDCk0LhQjD44sstzarmvRuEJjA0RUtKCgIIiiqPGqU6VSYerUqfjpp5/QrFkzbN++XettD5XAqJI0JcvOzsaZM2c0BvAXxMvLC1u2bNFYaiM/TZs2xYYNGzB8+HCkpaUhMTERCxYswIIFC/KUtba2xoYNG9C0adMSfwYiMixBELBseBO8u+FKvomaep00/sOLqGjq8WjqFSOePn2KQYMG4fTp0xg9ejSWLVtWrL3BlUARY9KMWe3atdG3b184OTkVWbZatWr4/vvvcfXqVTRq1Eir9vv164crV65oLEWSmyAI6NixIwIDA6WtsYjIeDjamGPDuObYMK45hrzhho51ymPIG27SNS5kS6Sd3DM7z507h8aNG+PcuXNYsmQJVq1apZGgLVu2DO3bt4eLiwucnZ3RvXt3XL9+Xa7QC8SetFKqXbs2duzYAQC4d+8ebt26hUePHiE+Ph5ZWVlwcHBA+fLl8cYbbxS4W0JR6tSpg6NHj+Lhw4c4e/YsIiMjAQCurq5o1aoV3NzcdPZ5iMjwBEFAa08XtPZ0KbowEeUrMDAQZmZmOHPmDL788ks4OTnh2LFjaNOmTZ6yCxYsQI8ePTBlyhSYmJjgp59+Qtu2bXHjxg1ZxrcXRBFJ2osXL1C7dm3Ex8ejTJkyuH79er6LvBZk4cKF+PbbbwEA3333HSZOnKivUAvl4eEBDw8PvbXv5uaGwYMH6619IiIiY5SRkYGbN29CpVLh008/haenJ06cOFHghL4rV65obJ/o5+cHd3d3rFu3Dp9//rmhwi6SIl537t69G+Hh4UhISMCQIUOKlaABwNixY5Geno74+HgsX75cT1ESERGREl27dg0qlQpdunRBxYoV8eDBg3x3KFJ7eX9rW1tbeHh45Ls9pZwUkaTlXuNr2LBhxa5fpkwZdO/eHUDOmmLF2ZSciIiIjJt6PFq3bt2wa9cumJiYoH///rh3755W9RMTE3Hjxg14enrqM8xiU0SSpp6RUbZsWXh7e5eoDT8/P+n4ypUrOomLiIiIlE+dR/j4+KBZs2ZYuXIlYmJi0KtXLyQkJBRZ/+uvv4aZmRnefvttfYdaLIpI0sLDwyEIAry8vErcRs2aNaXj+/fv6yIsIiIiMgLqnjQfHx8AOW/lpk2bhlu3bmHQoEHSgu/52bVrF37++Wf88ssv0s5GSqGIiQMpKSkAUKr9Ju3s7KTjpKSkUsdERET0uvJy8sLarmu1Kie3rKwshISEwM3NTWM5rDlz5uDmzZvYsWMHPvzwQyxZsiRP3cDAQAwfPhxTp07FW2+9ZciwtaKIJK1MmTKIiYlBTEzBq24XJXfd3AkbERERFY+9hb1idhIoyu3bt5GWlib1oqkJgoD169fD19cXS5cuRe3atfH+++9L9yMiItCjRw906dIF3333naHD1ooikrSKFSsiOjoaN2/eREpKSol61C5cuCAdV6hQQZfhERERkULVq1cPoijme8/W1lZ6FZpbQkICunfvDjc3N6xfv16xu3ooYkxa69atAeRsQr52bdHdqy/LyMjAn3/+KZ23aNFCZ7ERERHRq2XAgAEIDw/HV199hZCQEJw/fx7nz5/HzZs35Q5NgyKStN69e0vHM2fOxL///lus+p988gkePnwIQRDg7e0Nd3d3XYdIREREr4hjx44hNTUV/fr1Q8uWLaWf9957T+7QNCgiSevWrZu0IWpsbCz8/f2xd+/eIuvFx8dj5MiRWLZsmXRtxowZeouTiIiIjJ8oivn+nDx5Uu7QNChiTBoA/Pbbb/Dz80NaWhqePn2KPn36oEGDBujfvz+aNGmCcuXKwdLSEgkJCQgNDcXp06fx999/IzU1FaIoQhAE9O3bF2+++abcH4WIiIio1BSTpDVp0gSbN2/GkCFDkJaWBlEUERISgpCQkALrqJMzAGjXrh02btxoqHCJiIiI9EoRrzvVevbsiYsXL6JJkyYA/uuOzH2c+wcArK2tMWvWLBw+fBiWlpayxU5ERESkS4rpSVOrW7cuLl68iBMnTmDTpk04c+YM7ty5o1HG3t4eLVq0QKdOnTBmzBiULVtWpmiJiIiI9ENxSZpau3bt0K5dOwA5vWhxcXHIyMhAmTJlYG1tLXN0RERERPql2CQtN0EQ2FtGRERErxVFjUkjIiIiohyK7kn7999/ERgYiOjoaCQkJCA7OxszZ86UOywiIiIivVNckpaUlITFixdj+fLlePz4cZ77+SVpQ4YMwYMHDyAIArZu3QpXV1dDhEpERESkN4p63XnhwgX4+Phg5syZePz4cb5LbuSnVatW0r5b69atM2DERERERPqhmCQtMDAQnTt3RkREhJSUeXp6om/fvqhcuXKhdUeOHAlzc3MAwN9//22IcImIiIj0ShFJmkqlwtChQ5GUlAQAaNiwIc6fP487d+5g+/btqF+/fqH1HR0d0a5dO4iiiKCgIMTGxhoibCIiIiK9UUSStn79eoSGhkIQBDRq1AgBAQFo1qxZsdpo2bIlAEjbSREREREZM0UkaTt37pSOly9fDhsbm2K34e3tLR3fvXtXF2ERERERyUYRSVpQUBAAwN3dHU2bNi1RG7kXu42Pj9dBVERERGRMjhw5AkEQMGHCBLlD0QlFJGnPnz+HIAioXr16idswM/tvNRGVSqWLsIiIiMiIXL16FQDQpEkTmSPRDUUkaVZWVgCAjIyMErcRHR0tHXMLKSIiIh1LVf6kvMDAQABM0nSqQoUKEEURoaGhJW7jwoUL0rGbm5suwiIiIiIAeJEKrO2d86uCXb16Febm5kWuCmEsFJGkqWdmRkdHIyAgoNj1VSoVNm3aBAAwNTWFr6+vTuMjIiJ6rV1ZDTy9lvOrQiUnJyM0NBTe3t6wtLTUuLdlyxbY2dmhYsWKJcoz5KKIJK1Pnz7S8dSpU5GVlVWs+vPmzcPDhw8hCAI6dOgAe3t7XYdIRET0enqRCgQsyjkOWKTY3rSgoCCIoqjxqlOlUmHKlCkYMmQI6tWrhytXrhhVR45ikrSGDRsCyHlt+eabbyIxMbHIeqIoYt68eZg9e7Z0bcaMGfoKk4iI6PVzZTWQ8iznOOWZYnvT1OPRGjduDAB4+vQpOnTogJ9++gmjR4/G6dOnjW5vb0UkaQCwYsUKaX203bt3w9PTE1OnTsWBAweknQiAnPfNBw8exKxZs1CzZk3MmDED2dnZEAQBkyZNQqtWreT6CERERK+W3L1oagrtTcs9s/PcuXNo3Lgxzp07hyVLlmDVqlUar0DXrl2Lpk2bokyZMrCzs0OTJk2wdetWuUIvkFnRRQyjSZMm2LZtGwYNGoTk5GTExMRg4cKFWLhwoVRGFMU866ipN17v378/Fi1aZMiQiYiIXm25e9HU1L1pLSfJE1MBAgMDYWZmhjNnzuDLL7+Ek5MTjh07hjZt2uQpGxcXh759+6Jhw4awsrLCzp07MXjwYNjY2KBnz54yRJ8/QVRnOQpx9+5djB07FmfOnJGuCYIgHb8crp2dHaZPn46pU6caLEb6T5UqVRAZGQlXV1c8evRI7nCIiIyStn+WpqamIiIiQjp3d3cv0S49WnmRCvzskzdJAwDb8sCHwYCFnp5dTBkZGbCzs5PWSfX09MSJEydQpUoVrdvw9fWFm5ubNBGxIIb8PVDM6041T09PnDp1CufPn8ekSZPQoEEDmJiYQBRFKUGzsbFB+/btsWDBAkRERDBBIyIi0rX8etHUFDY27dq1a1CpVOjSpQsqVqyIBw8eICwsrFhtODs748WLF3qKsGQU87rzZc2aNdPYZD0hIQEpKSlwdHSEra2tjJERERG94vIbi/aygEVAk9GK6E1Tj0fr1q0bZs+eDT8/P/Tv3x8XLlyAh4dHgfVUKhVSU1Nx4MABHDlyBNu3bzdUyFpRXE9aQRwdHVG5cmUmaERERPpWWC+amoJ609QzO318fNCsWTOsXLkSMTEx6NWrFxISEvKt8+TJE5ibm8PR0RHDhw/HokWL0LVrV0OGXSSjSdKIiIjIALTpRVNTyExPdU+aj48PAGDYsGGYNm0abt26hUGDBuW7/qqLiwsuXbqE48eP45NPPsH777+vuJ40xb7uJCIiIhlkZwKD1hWvvIyysrIQEhICNzc3ODk5SdfnzJmDmzdvYseOHfjwww+xZMkSjXpmZmbSihHt2rVDbGwspk2bhv79+xs0/sIwSSMiIqL/WDkC7i3ljkJrt2/fRlpamtSLpiYIAtavXw9fX18sXboUtWvXxvvvv19gOw0bNsT69ev1HW6xGCxJq1GjhkGeIwgC7t27Z5BnERERkbzq1auXZ3kuNVtbW+lVaFHOnTuHatWq6TCy0jNYkhYeHq6x3llB8vuiX65XUBlRFLV6BhGVnCiK+OduDPaGPEZ0cgZc7CzRs0FltPZ05v9/RGQU2rVrhwEDBqB27dpIT0/Hnj17sGHDBqxatUru0DQY9HWntuvmvrx4bUH1Clvkloh0LyE1E+/8eQXnwmI0rm++9BCtPJyxbHgTONqYyxQdEZF2fHx88Msvv+DBgwews7ND3bp1sW/fPnTv3l3u0DQYLEk7ceJEkWUePXqEjz/+GDExMRBFEfXr10e3bt1Qr149ODs7w9LSEklJSbh//z4uXryIffv2ITk5GYIgYMiQIZg4caIBPgnR60kUxXwTNLWz92Lw7oYr2DCuOXvUiEjRFi1aZBRbSRosSfPz8yv0/tWrV/HRRx8hNjYWderUwbJly/Ldbyu3lJQUzJs3DwsWLMCWLVtgaWmJ1auVsWYL0avmn7sxBSZoamfvxeDsvRi09nQxUFRERK8uRayTlpCQgP79+yMmJkbatb6oBA3IGRD47bffYv369RBFEevWrcOPP/5ogIiJXj97Qx5rVW5PsHbliIiocIpI0tasWYOIiAgIgoBVq1bB3t6+WPWHDh2Kfv36QRRFzJ8/X9pglYh0Jzo5Q6fliIiocIpI0rZu3QoA8PLyQv369UvUxsCBAwEAMTExOH78uM5iI6IcLnaWOi1HRESFU0SSFhYWBkEQ4OrqWuI2cte9f/++LsIiolx6NqisVblePtqVIyKiwikiSYuPjwcAPHtWxGauhchdt6DNVImo5Fp7OqOVh3OhZVp5FF2GiIi0o4gkrVKlShBFETdv3kRYWFiJ2ti1a5d0XLFiRV2FRkT/TxAELBvepMAkTL1OGpffICLSDUXs3enr64vw8HAAwIQJE3DgwAGYm2u/IObRo0exceNG6bx169a6DpGIADjamGPDuOY4ey8Ge4L/23Ggl09ltPLgjgNERLqkiCRt/Pjx+PPPPwHkLHrbsWNHrFixAjVr1iy0XnZ2NpYtW4bPPvsM2dnZEAQB/v7+8PDwMETYRK8lQRDQ2tOFa6EREemZIpK0Nm3aYPz48VixYgUEQUBAQAC8vb3Rrl07dOnSBd7e3nB2doaFhQWSkpIQHh6OixcvYvv27YiMjJS2hLKxscHy5ctl/jREREREpaeIJA0Afv31VyQlJWHz5s0QBAEqlQpHjx7F0aNHC6yTe0N1BwcH7N27F15eXoYKmYiIiEhvFDFxAABMTU2xceNGrFixAs7OOQOTc2+urj7Ofa7+tUePHggODoavr688wRMRERHpmGKSNLWxY8fi4cOH+PPPP/Hmm2/C3d1dSsjUrKys0LJlS0ydOhXXrl3Dnj174O7uLlPERERERLqnmNeduVlaWmLYsGEYNmwYAEClUiE+Ph4vXryAg4MD7OzsZI6QiIiISL8U15OWHzMzM7i4uKBy5cpM0IiIiChfR44cgSAImDBhgtyh6IRRJGlERERERbl69SoAoEmTJjJHohtM0oiIiOiVEBgYCIBJGhEREb0mMp8+xYOJE5H59KncoRTq6tWrMDc3R/369eUORSeYpBEREVGBMp8+RcTbbyPl1GlEvP22YhO15ORkhIaGwtvbG5aWlhr3tmzZAjs7O1SsWBEBAQEyRVh8TNKIiIgoX+oELTPiQc55xAPFJmpBQUEQRVHjVadKpcKUKVMwZMgQ1KtXD1euXDGqNVWZpBEREVEeLydo0nWFJmrq8WiNGzcGADx9+hQdOnTATz/9hNGjR+P06dNwdXWVM8RiY5JGREREGgpK0KT7CkzUcs/sPHfuHBo3boxz585hyZIlWLVqVZ5XoGrXr1+HmZkZqlWrZsBotcMkjYiIiCRFJWhSOYUlaoGBgTAzM8OZM2fg5+cHlUqFY8eOYdKkSYXWmzJlirQdpdIwSSMiIiJJ1MyZRSZoapkRDxA1c6aeIypaRkYGbt68CZVKhU8//RTu7u64cuUK2rRpU2i9ffv24c6dOxg9erSBIi0eJmlEREQkqTR7Nszdq2pV1ty9KirNnq3niIp27do1qFQqdOnSBRUrVsSDBw8QFhZWaB11Qvftt9/CysrKQJEWD5M0IiIikphXqAD3deuKTNTM3avmlKtQwUCRFUw9Hq1bt27YtWsXTExM0L9/f9y7d6/AOsuWLYOtra20T7gSMUkjIiIiDUUlakpK0ID/Znb6+PigWbNmWLlyJWJiYtCrVy8kJCTkKR8XF4dvvvkGP/zwAwRBMHS4WmOSRkRERHkUlKgpLUED/utJ8/HxAQAMGzYM06ZNw61btzBo0CBkZWVplJ89ezZatmwJf39/Q4daLEzSiIiIKF8vJ2pKTNCysrIQEhICNzc3ODk5SdfnzJmDfv364fDhw/jwww+l6+Hh4Vi2bBm++uorxMfHIz4+Hunp6cjOzkZ8fDxevHghx8fIl5mhHlSjRg2DPEcQhELfQRMREZH21Ila1MyZOZMKFJSgAcDt27eRlpYm9aKpCYKA9evXw9fXF0uXLkXt2rXx/vvvIzw8HBkZGWjVqlWetpycnLBixQqMGzfOUOEXymBJWnh4uFbvfUVRzHPt5XoFlRFFUdHvlomIiIyReYUKqPrbb3KHka969erlmxcAgK2trfQqVK1hw4Y4ceKExrU1a9bg4MGD2Lx5M2rVqqW3WIvLYEkakH9ylZ/ciZYoigXWe7kcERERUWHKlCmTZyzayZMnYWVlpbgxagZL0l7OWvPz6NEjfPzxx4iJiYEoiqhfvz66deuGevXqwdnZGZaWlkhKSsL9+/dx8eJF7Nu3D8nJyRAEAUOGDMHEiRMN8EmIiIiI9M9gSZqfn1+h969evYqPPvoIsbGxqFOnDpYtW1bkSsEpKSmYN28eFixYgC1btsDS0hKrV6/WZdhERET0ips1axZmzZoldxh5KGJ2Z0JCAvr374+YmBhpQ9SiEjQg513zt99+i/Xr10MURaxbtw4//vijASImIiJjkPn0KR5MnKiY/SWJikMRSdqaNWsQEREBQRCwatUq2NvbF6v+0KFD0a9fP4iiiPnz50OlUukpUiIiMhbqjcJTTp1W1EbgRNpSRJK2detWAICXlxfq169fojYGDhwIAIiJicHx48d1FhsRERkfdYKm3ig8M+IBEzUyOopI0sLCwiAIAlxdXUvcRu669+/f10VYRERkhF5O0KTrTNTIyCgiSYuPjwcAPHv2rMRt5K6b3z5dRET06isoQZPuM1EjI6KIJK1SpUoQRRE3b95EWFhYidrYtWuXdFyxYkVdhUZEREaiqARNKsdEjYyEIpI0X19f6XjChAnIzMwsVv2jR49i48aN0nnr1q11FhsRERmHqJkzi0zQ1DIjHiBq5kw9R0SvIkMunq+IJG38+PHS8YkTJ9CxY0fcuXOnyHrZ2dlYunQpevfujezsbAiCAH9/f3h4eOgzXCIiUqBKs2dLG4EXxdy9KirNnq3niHTP1NRU4zw9PV2mSF5fGRkZGucv/57okkG3hSpImzZtMH78eKxYsQKCICAgIADe3t5o164dunTpAm9vbzg7O8PCwgJJSUkIDw/HxYsXsX37dkRGRkpZrY2NDZYvXy7zpyEiIjmoNwIv6pWnuXtVuK9bp7iNwrVhYWEBU1NTZGVlAfhvPLalpSX3rtYzURSRkZGhMQbezMwMFhYWenumIpI0APj111+RlJSEzZs3QxAEqFQqHD16FEePHi2wTu4N1R0cHLB37154eXkZKmQiIlKYohI1Y07QgJw9qx0cHBAXFwcg5+/BpxxbJxt7e3u9JseKeN0J5HQXbty4EStWrICzszMAzc3V1ce5z9W/9ujRA8HBwRpj24iI6PWkTtRefvVp7AmaWrly5WBtbS13GK89a2trlCtXTq/PUEySpjZ27Fg8fPgQf/75J9588024u7vnGaRnZWWFli1bYurUqbh27Rr27NkDd3d3mSImIiKleTlRe1USNCCnU8PNzQ1OTk56HQ9F+TMzM4OTkxPc3Nz0/v0LoiGnKZSQSqVCfHw8Xrx4AQcHB9jZ2ckdEv2/KlWqIDIyEq6urnj06JHc4RARach8+hRRM2fmTCpQcIJW0j9LRVHEixcvpDFqpF+mpqawsLAw2Pg/xYxJK4yZmRlcXFzkDoOIiIyMeYUKqPrbb3KHoTeCIMDS0lLuMEhPFPe6k4iIiIiYpBEREREpkqJfd/77778IDAxEdHQ0EhISkJ2djZlcIZqIiIheA4pL0pKSkrB48WIsX74cjx8/znM/vyRtyJAhePDgAQRBwNatW+Hq6mqIUImIdC81FrApK3cURKQAikrSLly4gKFDhyIiIgJA3v2xCppN0apVK2zduhWCIGDdunX48ssv9R4rGQ9RFPHP3RjsDXmM6OQMuNhZomeDymjt6cwVuklZXqQCa3sDYw8DFjZyR0NEMlPMmLTAwEB07twZERER0qK1np6e6Nu3LypXrlxo3ZEjR8Lc3BwA8PfffxsiXDISCamZGLbiAkasvIDNlx7i6K1n2HzpIUasvIDhf1xAQmqm3CES/efKauDptZxfiei1p4ieNJVKhaFDhyIpKQmCIKBhw4ZYvnw5mjVrBgDo1q1bvq8+1RwdHdGuXTscPnwYQUFBiI2NRdmyfF3wuhNFEe/8eQXnwmLyvX/2Xgze3XAFG8Y1Z48aye9FKhCwKOc4YBHQZDR700og6UUSQuNCiyzn5eQFewt7A0REVHKKSNLWr1+P0NBQCIKARo0a4fTp07CxKd4fTi1btsThw4chiiJCQkLg7++vn2DJaPxzN6bABE3t7L0YnL0Xg9aeXIePZHZlNZDy/xs3pzzLOW85Sd6YjFBoXChGHhxZZLm1XdeicYXGBoiIqOQU8bpz586d0vHy5cuLnaABgLe3t3R89+5dXYRFRm5vSMG9r7ntCdauHJHe5O5FUwtYlHOdiF5bikjSgoKCAADu7u5o2rRpidrI/XozPj5eB1GRsYtOztBpOSK9yd2LpqbuTSOi15YikrTnz59DEARUr169xG2Ymf335lalUukiLDJyLnbabZWibTkivcivF02NvWlErzVFJGlWVlYAgIyMkvdoREdHS8ecNEAA0LNB4bOC1Xr5aFeOSC/y60VTY28a0WtNEUlahQoVIIoiQkOLnpFTkAsXLkjHbm5uugiLjFxrT2e08nAutEwrj6LLEOlNYb1oauxNI3ptKSJJa9myJYCc3rCAgIBi11epVNi0aRMAwNTUFL6+vjqNj4yTIAhYNrxJgUlYKw9nLBvehMtvkHwK60VTY28a0WtLEUtw9OnTB2vWrAEATJ06FWfOnIGpqanW9efNm4eHDx9CEAR06NAB9vZc+4ZyONqYY8O45jh7LwZ7gv/bcaCXT2W08uCOAyQjbXrR1LhuGtFrSTFJWsOGDREcHIwLFy7gzTffxNq1a+Hg4FBoPVEU8d1332H27NnStRkzZug7XDIygiCgtacL10IjZcnOBAatK155InqtKCJJA4AVK1bA398fqamp2L17Nzw9PTFq1Ci0a9cOSUlJUrmrV6/i6dOnOH/+PDZs2ICwsDCIoghBEDBp0iS0atVKxk9BRKQlK0fAvaXcUbxyvJy8sLbrWq3KESmdIL68i7mMDh48iEGDBiE5OTnPayh1mAVdHzBgALZs2QITE0UMs3ttVKlSBZGRkXB1dcWjR4/kDoeIyCjxz1LKj6Iymq5duyIwMBBt2rSRNlnPnZwJgqBxXRRF2NnZYf78+di2bRsTNCIiInplKOZ1p5qnpydOnTqFixcvYv369Thz5gxu3LiBrKwsqYyNjQ1atGiBLl26YNy4cXBycpIxYiIiIiLdU1ySptasWTM0a9ZMOk9ISEBKSgocHR1ha2srY2RERERE+qeIJO306dPSsY+PDxwdHfOUcXR0zPc6ERER0atIEUmav78/BEGAlZUVHj9+LHc4RERERLJTxEh7CwsLiKKI2rVrs7eMiIiICApJ0ipWrAhBEFCmTBm5QyEiIiJSBEUkaV5eXhBFkWvDEBEREf0/RSRpAwcOBADcvXsXoaGhMkdDREREJD9FJGnDhg1DtWrVAACTJ09Gdna2vAERERERyUwRszvt7OywZcsWdOrUCYcPH0afPn3w+++/o1KlSnKHRkSkN0kvkhAaV/TbAy8nL9hb2BsgIiJSEkUkaep10hYuXIgpU6Zg//79qFGjBnr06AFfX194eHjA3t5e622f2rZtq89wiYh0IjQuFCMPjiyy3Nqua9G4QmMDRERESqKIJE29TpqaKIrIyMjAjh07sGPHjmK1JQgCVCqVrkMkIiJjlRoL2JSVOwqiYlPEmDS1lzdTV18r7g8REREA4EUqsLZ3zq9ERkYRPWlVq1bV6EkjIiLSiSurgafXcn5tOUnuaIiKRRFJWnh4uNwhEBHRq+ZFKhCwKOc4YBHQZDRgYSNnRETFoogkjYiUjzMRyehcWQ2kPMs5TnnG3jQyOkzSiEgrnIlIRiV3L5oae9PIyChq4gAREZFO5O5FU1P3phEZCUX3pP37778IDAxEdHQ0EhISkJ2djZkzZ8odFhERKVl+vWhq7E0jI6K4JC0pKQmLFy/G8uXL8fjx4zz380vShgwZggcPHkAQBGzduhWurq6GCJWIqFS8nLywtutarcpRMeTXi6bGsWlkRBSVpF24cAFDhw5FREQEAORZ86ygZTpatWqFrVu3QhAErFu3Dl9++aXeYyUiKi17C3uO39O1wnrR1NibRkZCMWPSAgMD0blzZ0REREiL0np6eqJv376oXLlyoXVHjhwJc3NzAMDff/9tiHCJiEiJCutFU+PYNDISikjSVCoVhg4diqSkJABAw4YNcf78edy5cwfbt29H/fr1C63v6OiIdu3aQRRFBAUFITY21hBhExGRkmjTi6YWsIi7EJDiKeJ15/r16xEaGgpBENCoUSOcPn0aNjbF64Zu2bIlDh8+DFEUERISAn9/f/0ES0REypSdCQxaV7zyRAqmiCRt586d0vHy5cuLnaABgLe3t3R89+5dJmlERK8bK0fAvaXcURDpjCKStKCgIACAu7s7mjZtWqI2ypYtKx3Hx8frICoiyo0zEYmIDEsRSdrz588hCAKqV69e4jbMzP77KCqVShdhEVEunIlIRGRYipg4YGVlBQDIyMgocRvR0dHSce5eNSIiIiJjpIgkrUKFChBFEaGhRW/eXJALFy5Ix25ubroIi4iIiEg2ikjSWrbMGegZHR2NgICAYtdXqVTYtGkTAMDU1BS+vr46jU8Jzp49iwkTJqBu3bpwcHCAg4MD6tatiwkTJuDs2bNyh0dEREQ6pogkrU+fPtLx1KlTkZWVVaz68+bNw8OHDyEIAjp06AB7e3tdhyiblJQUjB07Fq1bt8aKFStw69YtJCUlISkpCbdu3cKKFSvQunVrjB07FikpKXKHS0RERDqimCStYcOGAHJeW7755ptITEwssp4oipg3bx5mz54tXZsxY4a+wjS4rKws9O/fH6tWrZKuWVtbo2nTpmjRogUcHByk66tWrcKAAQOKneASERGRMikiSQOAFStWSOuj7d69G56enpg6dSoOHDgg7UQAAFevXsXBgwcxa9Ys1KxZEzNmzEB2djYEQcCkSZPQqlUruT6Czs2YMQOHDx+WzsePH49Hjx7h0qVLOHfuHB4/fozp06dL9w8dOpTvBvRERERkfATx5V3MZXTw4EEMGjQIycnJeTZTV4dZ0PUBAwZgy5YtMDFRTN5ZKpGRkfD09ER6ejoA4K233sK6dfmvpD1jxgzMnTsXQE5P2927d4vc71RXqlSpgsjISLi6uuLRo0cGeSYR0auGf5ZSfhSV0XTt2hWBgYFo06aNtMl67uRMEASN66Iows7ODvPnz8e2bdtemQQNABYvXiwlaDY2Nli0aFGBZWfMmCHNaE1LS8PPP/9siBCJiIhIjxSX1Xh6euLUqVM4f/48Jk2ahAYNGsDExEQjYbOxsUH79u2xYMECREREYOrUqTJHrXvbt2+XjgcNGlTo2m8WFhYYPXq0dL5jxw69xkZERET6p4gdB/LTrFkzNGvWTDpPSEhASkoKHB0dYWtrK2Nk+vfvv//i7t270nnXrl2LrNOtWzdpAkVoaCju3LmDmjVr6i1GIiIi0i/F9aQVxNHREZUrV37lEzQACA4O1jhXryNXmMaNG8PCwqLANoiIiMi4GE2S9jq5deuWdGxhYaHVDgovl8vdBhERERkfJmkKFBERIR1XqVIlz4zWglStWlU6Dg8P13VYREREZECKHZMG5IzNCgwMRHR0NBISEpCdnf1arAOWeyFfR0dHrevlXtw299pyxfXjjz/ixx9/1KpsZGQkACAqKgpVqlQp8TOJiF5nUVFRAIBnz57JHAkpieKStKSkJCxevBjLly/H48eP89zPL0kbMmQIHjx4AEEQsHXrVri6uhoiVL3Jvb2TlZWV1vWsra3zbaO4EhMTpeRLW9nZ2cWuQ0REmrhrDOWmqCTtwoULGDp0qPS67+V1dgt67deqVSts3boVgiBg3bp1+PLLL/Ueqz5lZmZKx2Zm2v8W5S774sWLEj/fwcFB60RXnUibmZmhfPnyJX6mIURFRSE7OxsmJiaoVKmS3OEYPX6fusfvVPeM5Tt99uwZsrKyivUPc3r1KSZJCwwMROfOnZGcnCwlZ56envD29sbFixfz7VVTGzlyJD777DOoVCr8/fffRp+kqbfHAiAtaKuN3GVLMwt2ypQpmDJlSonrK5V6Re9KlSpxRW8d4Pepe/xOdY/fKRkzRUwcUKlUGDp0qDSOqmHDhjh//jzu3LmD7du3o379+oXWd3R0RLt27SCKIoKCghAbG2uIsPXGzs5OOk5LS9O6Xmpqar5tEBERkfFRRJK2fv16hIaGQhAENGrUCAEBARoL2WpDvZaYKIoICQnRR5gG4+LiIh2rB5Nq48mTJ9Kxs7OzTmMiIiIiw1JEkrZz507pePny5Rqv+7Tl7e0tHederd8Y1apVSzqOiYnR6CErzMOHD6Xj2rVr6zwuIiIiMhxFJGlBQUEAAHd3dzRt2rREbeTe2zI+Pl4HUcmnTp06Gufq76cwkZGReP78eYFtEBERkXFRRJL2/PlzCIKA6tWrl7iN3DMbVSqVLsKSTbNmzWBpaSmdBwQEFFnnzJkz0rGVlVWxXxcTERGRsigiSVNPOc7IyChxG9HR0dJx7l41Y2RnZ4cOHTpI5xs2bCiyTu4yHTp0eC32OCUiInqVKSJJq1ChAkRRRGhoaInbuHDhgnSszV6XSjdq1CjpOCQkBHv27CmwbGBgIA4cOJBvXSIiIjJOikjS1DMzo6OjtXq19zKVSoVNmzYBAExNTeHr66vT+OTw5ptvwsfHRzqfOHEibt++nadcVFQURowYIa1S3bBhQwwYMMBgcRIREZF+KCJJ69Onj3Q8derUYm+LMW/ePDx8+BCCIKBDhw6wt7fXdYgGJwgC/vjjD2mrp6ioKDRv3hxffPEF9u/fj8OHD2Pu3Llo1KgRbt26BSBnW6gVK1ZovSE7ERERKZcidhzo06cPGjZsiODgYFy4cAFvvvkm1q5dq7FheH5EUcR3332H2bNnS9dmzJih73ANpmnTptiwYQOGDx+OtLQ0JCYmYsGCBViwYEGestbW1tiwYUOJZ8e+DqZMmYLExMQi/7si7fD71D1+p7rH75SMmSC+vEGmTK5cuQJ/f39pTTBnZ2eMGjUK7dq1w7fffouzZ89CEARcvnwZT58+xfnz57FhwwaEhYVBFEUIgoBJkyZh8eLFMn8S3bt16xY++OADHD9+PN/9TDt06IBffvmFa6MRERG9QhSTpAHAwYMHMWjQICQnJ+d5ZacOs6DrAwYMwJYtW2Bioog3uHrx8OFDnD17FpGRkQAAV1dXtGrV6pWYKEFERESaFJWkATm7BYwdO1Zj3a/cidnL4drZ2WH69OmYOnWqwWIkIiIi0jfFJWlqFy9exPr163HmzBncuHFDYzKBjY0NWrRogS5dumDcuHFwcnKSMVIiIiIi3VNskvayhIQEpKSkwNHRkQu1EhER0SvPaJI0IiIiotfJqzvKnoiIiMiIMUkjIiIiUiAmaUREREQKpIgdB3K7d+8e1q9fj3/++Qe3b99GfHy8tMCtNgRBgEql0mOERERERPqnmCRNpVJhypQpWLZsGbKzs6XrnNdAREREryPFJGmDBw/Gzp07mZQRERERQSFJ2o4dO7Bjxw5pZ4HKlStj3Lhx8PX1haurK2xsbGSOkIiIiMiwFJGkrV27Vjpu3749du3axQVrqVSeP3+Oy5cv49KlS9KvT548ke6vXr0ao0aNki9AIxMfH48TJ07gxIkTCAoKwp07dxAXFwdzc3OULVsWPj4+6NChA0aOHMkdQLSQmZmJwMBAnD17FiEhIbh9+zYiIiIQHx+PzMxMODg4wM3NDU2bNsXAgQPRqVOnV3pfYn0LDw+Ht7c3UlJSpGtff/01Zs2aJV9QRNoQFaBy5cqiIAiiubm5GBkZKXc4ZMSioqJEd3d3EUChP6tXr5Y7VKNw69YtsWfPnqKFhUWR3ykA0cbGRvzpp5/E7OxsuUNXtE8//VSr71P907BhQzEwMFDusI1W586d83ynX3/9tdxhERVJET1pMTExEAQBDRo0QOXKleUOh4xYeno6IiIi5A7jlXH9+nXs3btX45qpqSk8PT1RoUIFZGVl4datW4iNjQUApKam4uOPP8b169exYsUKaQgDaRJfGntra2sLDw8PODk5QRAEPHnyBHfu3JEmUQUFBaFt27Y4ePAgWrduLUfIRuvPP//E4cOH5Q6DqEQU0X/u4uICAChTpoy8gdArpVy5cujatSumT5+OXbt2yR2OUTMzM0Pfvn2xc+dOxMbG4vbt2zh16hQCAgIQHR2NnTt3wtXVVSq/cuVKLF++XMaIlc3a2ho9e/bE77//jtu3byM5ORnBwcE4efIkTpw4gVu3buHJkyeYNm0aTE1NAQDJyckYNmyYxis7Klx0dDQ+/vhjAECdOnXYCUBGRxE9abVq1cLjx48RGRkpdyhk5MqWLYtt27bhjTfegLu7u9zhGD1zc3OMGzcOM2bMQNWqVfMtIwgC+vTpg8aNG6NZs2bS2L+ZM2di3LhxMDc3N2TIRmHOnDlFlilXrhzmzp2LatWqYfz48QCABw8eYOvWrRg9erS+Q3wlfPzxx4iOjgYALF++HG+//bbMEREVjyJ60kaMGAEAuHPnDsLCwmSOhoyZg4MD3nzzTSZoOtKnTx+sWLGiwAQtNzc3N3zzzTfSeXR0NE6fPq3P8F4L48aNg4eHh3R+8uRJ+YIxIocPH8aff/4JABg9ejTatm0rc0RExaeIJG348OHw8fEBAHz66acyR0NEJdWrVy+N89u3b8sUyaulcePG0nHuWcqUv9TUVLzzzjsAcobT/O9//5M5IqKSUUSSZmFhgb/++guVKlXCrl278PbbbyM5OVnusIiomMqWLatxnpiYKFMkr5bcW905ODjIGIlxmDFjBu7fvw8A+OGHH+Ds7CxzREQlo4gxaepXIgsXLsT777+PDRs2YN++fRgyZAhatmyJihUrwsLCQuv22K1NJI+XZ9aWL19epkheHZmZmTh37px03rJlSxmjUb4rV67g559/BgD4+/tj5MiRMkdEVHKKSNL8/f01puqLooi4uDgsX7682DPEuME6kXy2b9+ucd6iRQuZInl1TJs2TXrFWbZsWS7CXAiVSoVx48YhKysLFhYWWLZsmdwhEZWKIpI0NVEUIQiClLC9vJYQESlXQkKC1IMBAA0aNEC9evVkjMg4qVQqPH/+HBcuXMCvv/6KI0eOAACsrKywcePGPK+U6T8LFy5EUFAQAODzzz9H7dq15Q2IqJQUkaRVrVqVi14SGblPPvlEY1D73LlzZYzGuLi4uCAmJqbA+x07dsTChQvRoEEDA0ZlXO7duyfNLvby8sJXX30lc0REpaeIJC08PFzuEIioFFatWoWVK1dK54MHD84z05NKxtfXF++//z7q168vdyiKNnHiRKSlpQEAfv31V1hZWckcEVHpKSJJIyLjdebMGbz33nvSefXq1fHbb7/JGJHx6dChAxISEgAAGRkZGttCBQQEICAgAM2aNcOWLVtQrVo1eYNVoNWrV+PYsWMAcpZ06tixo8wREekGkzQiKrHg4GD06tULGRkZAHJmcx48eBCOjo4yR2ZctmzZkudabGws/vjjD8yePRspKSm4ePEi/Pz8cOnSJc6azeXZs2fS+ppOTk748ccfZY6ISHcUsU4aERmff//9F507d5Z6gJycnHD48GHUrFlT5sheDWXLlsXUqVNx5swZ2NvbA8jZFuqTTz6ROTJlmTx5MmJjYwEA8+fPZwJLrxQmaURUbPfv30fHjh3x7NkzAICdnR0OHDgg7RxCutOoUSN8+eWX0vnmzZulpOR1d+7cOakXsmXLltIep0SvCkE00DoX69at0zjPvdHty/dKi5voUn5yzyBevXo115sqoUePHqFt27bSiu7W1tY4cOAA/Pz8ZI7s1XXv3j14enpK54cOHULnzp1ljEgZdu7ciX79+umkrfv373O8HymOwcakjRo1SvpLUhAEjUQq973SerltItKdp0+fomPHjlKCZmlpiZ07dzJB0zM3NzeN8+joaJkiISJDMvjEgYI67rhwLZGyxcbGolOnTvj3338BAObm5ti6dSt7dAxAPe5PrUyZMvIEojCWlpbF2pczLi4O2dnZAHJ6gG1sbKR7pqamOo+PqLQMlqS1bdu2wN6ywu4RkfwSExPRpUsXXLt2DUDOX2gbNmxA7969ZY7s9XDmzBmNcw8PD5kiUZZu3boVq1exWrVq0v6yU6dOxaxZs/QUGZFuGCxJO3nyZInuEZG8UlNT0aNHD1y+fBkAYGJigtWrV2PgwIEyR/Z6ePHihcbuDR4eHqhVq5aMERGRoXB2JxEVKCMjA3369EFAQACAnDGfv//+O9566y2ZIzNeR44cwWeffYbIyMgiy0ZFRaFXr164evWqdO2LL77QZ3hEpCAGm91JZCjjx4/H+vXr81xXL7gKAGZmZvmOQUlPT9drbMbm+++/x+effy6dOzk5oVmzZlrX79SpE9f1eol6RqIgCGjVqhXatGmD+vXro1y5crCxsUFycjLCwsJw5swZ7Nq1C6mpqVLd3r17Y+fOnRweUkK5X3d+/fXXfN1JiscdB+iVk5mZqZGQ5UelUkGlUhkoIuOVO0EAcgZeHzp0SOv6FStW1HVIrwxRFPHPP//gn3/+0ar86NGjsXz5ciZoRK8Rvu4kIjKgpk2bYsqUKahbt26RCZeFhQUGDBiAU6dOYdWqVbCwsDBQlESkBHzdSUQkk/j4eAQHByMsLAzR0dHIyMiAra0tnJycUKdOHfj4+MDKykruMIlIJkzSiIiIiBRIEWPSatSoUeo2LCws4OjoCBcXF/j4+KBVq1bo1q0bFygkIiIio6SInjQTExNpbEZ+4bw8bkObMgBQvnx5fPXVV/jggw90FCkRERGRYSgmScutsIStuOUEQUCPHj2wY8cO9qoRERGR0VBEknbq1CkAwKNHj/Dxxx8jOjoagiCgbdu28Pf3R506deDo6AgTExPExcXh9u3bOHXqFE6dOoXs7GyUK1cOP/74I1xcXPDs2TNcunQJmzdvlrYLEQQB7733Hn755Rc5PyYRERGR1hSRpAFAUFAQOnXqhNjYWLRu3Rq///47ateuXWid27dvY+LEiThz5gycnZ1x5MgRNGzYEEDOoqSTJk3C6tWrAeQsXnrz5k14enrq+6MQERERlZoi1klLSkrCgAEDEBsbC39/fxw9erTIBA0AateujaNHj6Jdu3aIiYnBgAEDkJiYCACwsrLCypUr0bVrVwBAVlYW1q5dq9fPQURERKQrikjS1q5di/v37wMAli9fXqwFG83NzbFs2TIAQHh4ONasWaNxf86cOdIxN3InIiIiY6GIJG3z5s0AcnrGvLy8il2/Zs2aqFOnDgBgy5YtGveaNGmC8uXLQxRFhIWFlT5YIiIiIgNQRJIWFhYGQRBQqVKlErdRsWLFAhMxDw8PAEBsbGyJ2yciIiIyJEUkaXFxcQCA58+fl7gN9UzO+Pj4PPdsbW0BFL2kBxEREZFSKCJJq1ChAkRRxI0bN/Dw4cNi13/06BGuX78OQRBQvnz5PPcTEhIAAA4ODqWOlYiIiMgQFJGktWzZEkBOT9fkyZOLVVddR91Lpm4rt9DQUAiCADc3t9IHS0RERGQAikjSRo8eLR3v3r0b/fr1Q1RUVJH1nj59igEDBmDXrl3StTFjxmiUuX79uvQK1MfHRzcBExEREemZIjZY79y5MwYOHIht27ZBEATs3r0bhw8fRs+ePdGuXTvUrl0bZcqUAZDz6vL27ds4ceIE9u7di7S0NIiiCEEQMGDAAHTu3Fmj7W3btknH7du3N+THIiIiIioxxew4kJ6ejr59++Lw4cMQBEFKvAqjLiOKIjp16oRdu3bByspKo0y/fv3w5MkTmJiY4ODBg7C3t9fnxyAiIiLSCcUkaUDOrgALFy7E7NmzkZqaKl1XJ2IvHwOAtbU1ZsyYgc8++4wbqBMREdErQ1FJmlpCQgI2bNiA/fv349KlS3mW5nBxcUHTpk3RvXt3jBgxQnoVSkRERPSqUGSS9rKMjAxp8L+jo2OeV5pERERErxqjSNKIiIiIXjeKWIKDiIiIiDQxSSN6hfn7+0MQBAiCgDVr1sgdjtF53b+/f/75R/r8s2bNkjucUklKSkK5cuUgCAJ8fX3lDodIK0zS6JUzatQo6S8Wf39/ucMhMkrZ2dnSDjAuLi745JNPiqxz/vx56f899c+IESO0fuaaNWukenZ2dsjOzi5x/C+zt7fHF198ASAn+dyyZYvO2ibSF8Ulaffu3cOsWbPQqVMnuLm5wd7eHqamplr/mJkpYn1eogJVq1ZN+ovo5MmTeq9HVBKbNm1CYGAgAOCjjz7Sao3Jq1ev5rm2b98+qFQqrZ6Zu36DBg1gYqLbv6Lee+89uLi4AAC++uorreMikotikjSVSoXJkyejdu3amDNnDo4fP47Hjx8jJSUFoigW64eIiEouKysLs2fPBgDY2tri3Xff1apefklafHy81v+oyF2/YcOGWtUpDmtra0yaNAkAEBYWhrVr1+r8GUS6pJgkbfDgwVi6dCmysrKYcBHpyMmTJ6X/l0aNGiV3OEbndf3+/vrrL9y5cwcAMGzYMJQtW1arekFBQdKxnZ2ddLxjx44i64qiiODgYOm8UaNGWkZbPO+++670xmX+/Pn8e4YUTRHvBnfs2IEdO3ZI20BVrlwZ48aNg6+vL1xdXWFjYyNzhEREr4+ff/5ZOh47dqxWdVQqFa5duyadT5o0CQsWLAAA7N69G0uWLCl0q7+wsDAkJiZK5/pK0ipUqIAePXpg165duHv3Lvbv348ePXro5VlEpaWIJC13l3P79u2xa9cu2NrayhgREdHr6fr16zh37hwAwMPDA82bN9eq3u3bt5Geni6djxkzBps2bcKDBw/w6NEjXL58GW+88UaB9XO/6jQzM4O3t3cJP0HRhg0bhl27dgEAfvvtNyZppFiKeN156dIlAICpqSnWrVvHBI0UJS0tDbt378aHH36INm3aoGLFirC0tIStrS2qVq2Knj17YvHixUhOTi6wjfDwcGnQf0REhHS9Xbt2eWbD5Z6VWtJ6aiVZQkIURezfvx/vvPMO6tevj3LlysHc3BxlypRBw4YNMWbMGGzbtg0ZGRlatbVz506MHj0atWvXhpOTE6ytrVG1alX07t0bq1atQmZmplZxFVdGRgbWr1+P/v37o0aNGrCzs4OZmRkcHBzg5eWF7t27Y9asWTh//nyBr7wK+/5y/94U90ebWcdyfXcbNmyQjvv27at1vZdfdXp6eqJPnz7StaJeeeauX7t2bb3uLNO9e3dYWloCAA4ePIiYmBi9PYuoVEQFsLS0FE1MTMQmTZrIHQq9AkaOHCkCEAGIfn5+pWpr48aNor29vdReYT9ly5YVd+/enW879+/f16qNl+MuaT01Pz8/6d7q1auL/LyXL18WmzRpotWz3N3di2yrUaNGRbbj5eUlXr58WYvfDe0FBQWJNWvW1Pp727JlS77tFPb9Fff3prDfp5fJ+d15eXlJ7R85ckTrelOmTJHqtWrVShRFUTx69Kh0rU6dOoXW7969u1R2xIgRpfoM2mjfvr30vJUrV+r9eUQloYjXnS4uLoiKiuJG6aQ49+7dQ1JSknRevnx5VKtWDfb29khLS8OdO3cQHR0NAIiNjUXfvn2xc+dO9OrVS6Mda2trdOnSBQBw6tQp6bXQG2+8ke+g7AYNGpSqXkns3bsXgwcPRmpqqnTN3NwcdevWhbOzM1JSUhAaGorY2FgAkPbTzc/Bgwfx5ptvIiUlRbrm4uICLy8vWFpa4v79+1LPYGhoKNq1a4dDhw6hZcuWJY5fLSoqCh06dNDoHSlTpgxq1aoFe3t7pKamIioqCuHh4VIPWknW48r9e6ON3L9/hZHzu7t//z5CQ0MB5LxyLE6buXvC1DMz/fz84OTkhLi4ONy6dQt37txBzZo1i6yvr/Foufn5+eH48eMAgEOHDmHMmDF6fyZRscmdJYpizr9oBEEQa9euLXco9ArQZU/a3LlzxdatW4u///67GBkZmW+ZM2fOiM2bN5ee6eLiIiYmJhbYpru7u1T2xIkTWsdSknra9qTduHFDtLGxkco6OTmJixcvzvdzBAcHi59//nmBPWmhoaGinZ2d1FazZs3EkydPitnZ2RrlLly4oNFb5O7uLsbHx2v1uQrz/vvvS21WrlxZ3Lt3r5iVlZWnXFxcnLhp0yaxU6dO4ubNm/Ntq7g9kQVZunSpRg/YsmXL8i0n93e3evVqqc0GDRoUq27ZsmWluitWrJCujxgxQro+f/78fOs+ffpU4/s5fvx4qT6HNvbv3y89r3z58np/HlFJKCJJW7VqlSgIgmhiYiLeu3dP7nDIyOkySUtOTtaqXFpamtiiRQvpuUuWLCmwrBKTtDfeeEPjL6wbN24U2XZSUlK+19u2bSu11atXL/HFixcFtpGcnCw2bNhQKj979uwin1uUGjVqSO0dPXpUqzoqlSrf67pI0o4dOyaamZlJ7UyaNKnAsnJ/d7kT3KFDh2pdLzw8XCPJunTpknRv27Zt0vUWLVrkW//gwYMa9WNjY0v9WYry8OFDjWc+ePBA788kKi5FTBwYPnw4fHx8AACffvqpzNEQ/UfbSSxWVlb49ttvpXP1zDFjcPz4cWnyDgD8/vvvqFu3bpH1cq+DpXbhwgWcPn0aAODs7Ix169bB3Ny8wDZsbW2xfPly6Xz58uWlXrfq0aNH0nHr1q21qmNqalqqZxbk3r17GDhwoLSyfceOHbFo0aJ8yyrhu7t586Z07OHhoXW93K8qX56Z2bVrV2kSwIULFxAVFVVofXd3dzg5ORUj6pKpUqWKNHkA0PzsREqhiCTNwsICf/31FypVqoRdu3bh7bffLnSmHJES5V6q4PLlyzJGUjybNm2Sjhs0aKAxI6+41q9fLx2PHDlSq3GmzZs3h6enJwDg8ePHuH37domfD0BjVmBISEip2iqNxMRE9OrVSxrD5+Xlha1btxa4dZ0SvrvcM4grV66sdb3cy2fUqlVL4/fAzs4OHTp0AACIopjvP2D0vdNAQXJ/xtyfnUgpFDFxQP2vx4ULF+L999/Hhg0bsG/fPgwZMgQtW7ZExYoVYWFhoXV7bdu21Veo9Bq7f/8+jh07hpCQEDx//hxJSUkF7v0XFxeH1NRUo1iI+cyZM9Jx//79ddZW+/btta7n7e2Nu3fvAgACAwNRp06dEsfQpEkTnDhxAgAwYsQIrF69WuseNV3Jzs7GkCFDcOvWLQA5Exf27NlTaA+REr479SQYAMXqzSpq0H/fvn2xb98+AMDOnTvxzjvvaNzPnaQZYtKAmpOTE+7fvw8AeP78ucGeS6QtRSRp6rWI1ERRRFxcHJYvX67Rna8NQRC4aS7p1O3bt/Hhhx/iyJEjxXqdlJCQoPgkLTs7W9r+B8hJcEpKFEXcuHFDOp83bx5++eUXrermXqk+d6JQEh9++KGUpIWGhsLX1xdeXl7o2rUr2rZtizZt2qBChQqlekZRpk6digMHDgDIeZW6ZcsW1KpVq8DySvnucs8otba21rpeUT1hvXr1gomJCbKzs3HixAkkJibCwcFBeqY6yQQMm6Tl/oy5PzuRUigiSVMTRVFa7FF9TiSn06dPo1u3bhrLUmhLm8Ve5RYXF6fx/1m5cuVK3FZCQgKysrKk87Nnz5a4ndLo06cP5s2bh+nTp0tLa4SGhiI0NFRKfOrXr4+BAwdi3LhxqFSpUqme97K1a9di4cKF0vlPP/2Ezp07F1pHKd9dbtr++RsbG4sHDx5I5/klaRUqVECLFi1w9uxZvHjxAvv27cPQoUMBAMHBwRpLoBjydSf/jiGlU0SSVrVq1UL3dCOSQ2JiIgYOHCglaPb29hgzZgw6d+6MmjVromLFirC2ttYYdG5s/x2/nEjmHkhdXLrqiSjJmmUv+/LLL9G9e3f88MMP2LlzZ54xrteuXcO1a9cwf/58zJ49G5988kmpnwnkJFcTJkyQzidMmIAPPvigyHpK+e5sbW2lRE+bNd0AzVedQME9YX379pWSz507d0pJWu76ZcuWRdWqVYsZdcmlpaVJx9zphpRIEUlaeHi43CEQ5bFq1So8e/YMQM7YlQsXLsDLy6vA8rkXvTUWLw9OL01PjKOjo8b5pUuX0LRp0xK3V1o+Pj5Yv349MjMzcenSJZw8eRKnTp3C6dOnpQQkNTUVn376KUxMTPDxxx+X6nkPHjxAv3798OLFCwA5wziWLFmiVV2lfHcuLi7SfwNxcXFa1cn9qtPNzS3fRZaBnB7OqVOnAgAOHDiAjIwMWFpayjYeDdD8jKXpRSbSF0XM7iRSoiNHjkjHH374YaEJGgBERkbqOySds7Gxgb29vXSuXm2+JOzs7DTG4KkTXLmZm5ujVatW+Oqrr3Do0CFER0fjt99+g4uLi1Tm66+/LlVvVkpKCnr37i19Zg8PD/z111+FLqGRm1K+O3d3d+lY2/+etZ2ZWbNmTWlSQ1JSEo4dO1as+kDO68nffvsNTZs2ha2tLZydndGxY0ccPHgQAFCtWjUIgoCTJ09qFfvjx4+l49yfnUgpmKQRFSD3OBttejXOnTunVbsmJv/9b1ecMTElrVeUFi1aSMfqmda6aOv8+fOlaktfbG1tMWHCBGzbtk26lpSUVOJ4RVHE22+/jeDgYACAg4MDdu/eDWdn52K1o4TvLvfM0NyD+QtTnO2ccm/YvmPHDqhUKo0JE4XVz8rKwuDBg/HOO+/gypUryMjIQHZ2No4dO4bu3btrPdFC7dGjR1KvJwCt1gYkMjQmaUQFyMzMLFb5NWvWaFUu99iX3GNi9FWvKJ06dZKO//77b2ldr5LIvZflhg0bNAbDK42/v7/Ga8anT5+WqJ2ZM2di+/btAHJmcm7evLlEf+Er4bvLPbv3+vXrRZZPT0/XWJutqJ6w3Gvw7dmzBzdu3NAY+1ZY/f/973/Ytm0bBEHAnDlzEBcXh7i4OERGRmLIkCH45JNPirWMRu5ZseXKlYObm5vWdYkMxaBJ2rp167Bu3Tqpm1uX5s6di/79+2PAgAE6b5teT7ln/f3zzz+Flt22bZvWvVAVK1aUjrXtrShNvaKMGTNGetWWmpqKDz/8sMRtjRs3TtqJICwsDN99951OYtRWcXoYMzIyNBLxgsZSFWbLli2YO3eudP7999+jW7duxW4HkP+7AzTXmLx9+3aRi4pfu3ZNI5ksKklr1qyZtIDs06dP8euvv0r3rK2tUbt27XzrpaamSt/H559/junTp0uv6StXrow///wTbdq0KdYs7NwLTvv5+Wldj8iQDJqkjRo1CqNHj9aYnl6QGjVqoEaNGhgzZoxWbf/zzz/YuXMndu7cWcooiXLk/oN7yZIlGv/yzu3QoUMYNWqU1u02btxYOl69erXWg/VLWq8ozs7O0oBuAPjzzz8xadKkQmf3xcXF4YcffshzvWzZspg+fbp0PnPmTHzzzTdF9krGx8dj8eLFGDx4cAk+wX8iIiLQpk0b7Nmzp8j1Er/++mvpL3ULCwuN143aCAwMxOjRo6XzMWPGYMqUKcUP+v/J/d0BOWPp1NtBZWVlFfmPk9zjycqUKYPq1asXWl4QBPTu3Vs6X716tXRcv379ArfnOnz4MBITE2Fubo7PPvssz30TExN89dVXhT77ZadOnZKOc/diEimJImZ35ic8PByCIODJkydyh0JG7PTp0xpb1Gjj33//hbu7OyZMmIDvv/8eaWlpSEpKQqtWrfDuu++iXbt2sLW1xYMHD7Bjxw7pVde4cePwxx9/FNn+0KFD8b///Q+iKCIoKAiurq5o3LgxnJycpCU8vL29NXpoSlNPG9OnT8fp06dx/PhxAMCvv/6KPXv2YMSIEWjevDmcnZ2RnJyMO3fu4NSpU9i/fz8sLS3z3Wt36tSpuHjxIrZv3w5RFDFr1iysXLkSQ4cORfPmzeHi4gKVSoXY2Fhcv34d586dw4kTJ5CZmamxtVZJBQQEICAgAC4uLujRowfeeOMNVK9eHY6OjkhLS8OtW7ewadMmjTGEkyZN0mobptx2796t8dr5/v376Nq1q1Z1GzRogO+//z7Pdbm/OyBn14n//e9/AHL+AVJYApM7SVPvv1yUPn36SIuU505AC+uFu3LlCoCcRK6gHk9fX1+YmZlptZh5SkqKlICamZmVais0Ir0y5G7ugiCIJiYmYrdu3XRaVhRFsWvXrlIder2NHDlSBFDin/v370ttrVu3ThQEocg6bdq0EdPS0gps52XTpk0rtD0/Pz+d1PPz85PurV69utDvLTU1Vezbt6/W35Ojo2OBbWVmZoqTJk0q9nffvHnzQmMsyv3794v9zF69eonp6en5tlfY9/f111+X+L+xgn5/5fzu1IKDg6U2q1evXmjZFi1aSGU//PBDrdrPyMgQHRwc8sS/bNmyAutMnDhR+r0qTMWKFUUA4okTJwott23bNum53bt31ypuIjlw4gBRId566y3s2rWrwNc4Tk5OmDZtGo4fP16sHru5c+fi+PHjGDFiBGrVqgU7OzutFsItaT1tWFtbY/v27diyZQu8vb0LLdugQQPMmzevwPtmZmZYsmQJzpw5g86dOxf4GgvIeQXWqFEjzJkzR2PGZUlUqFAB33//Pdq0aVPkfr+1atXCihUrsGvXrlIt4qtrcn13ag0aNJB65e7fv1/gTNPs7GyNDey13SnAwsIi33F7htxpYOPGjdLxxIkTDfZcouISRNFw+2KYmJhAEAR06dIF+/fv11lZAOjWrRsOHToEQRAUPaOMjJNKpcK5c+cQHByMxMREuLi4oFq1avD39y8yGTBW6r+gnz59iuTkZNjb26N69epo0qQJXF1di9VWQkICAgIC8PDhQ8TGxsLMzAxlypSBp6cnGjRooLFmma6kp6cjODgYoaGhePLkCdLS0mBra4uKFSuiUaNGpdqI3JDk+O42b94s7Qgwbtw4rFixQufPKI4ZM2Zg7ty5aNKkicaA/9xevHgBW1tbqFQqnDhxAv7+/vmWe/bsGapUqYLMzEzUqFEDoaGhGsvbECmJYsekESmJmZkZ2rRpgzZt2sgdisFUr169yIHg2nJ0dESPHj100pa2rKys0Lx5c52N1ZKLHN/dwIEDMXPmTISGhmLTpk1YsGBBiWa/6op6aZCQkBDExcXByckpT5mAgACtxqP99ttv0li4zz//nAkaKRr/6yQiIg2mpqaYOXMmgJxB9suWLZM1ns6dO8PR0RGZmZn5zioWRVGrJUvS09OlrbqqVaumMTuXSImYpBERUR7Dhg2TdgBYtGiRrHvT2tjY4IsvvgAAfPfdd/juu++kNdyioqLw9ttv49SpUxpba+Vn2bJl0pZb8+bN03rbLiK5MEkjIqI8TExMpK2WoqOjtVrfUp8+++wzvPnmmxBFEV999RXKlCmDsmXLwtXVFRs2bMDChQsL3SQ9KSlJmuzSunVracwdkZIxSSMiony1bt0aoihKa7bJydTUFFu3bsXy5cvRuHFjacJO+/btsW/fPnzwwQeF1re3t8fz588hiiICAgIMETJRqXHiABERGQVBEDBx4kQum0GvDfakERERESmQLD1pFy9eRPv27XVaNveiikRERETGTpYkLS4uTmNz2/yoV1HXpiwRERHRq8bgSZoBNzggIiIiMloGTdJGjhxpyMcREdFrJDw8XO4QiHTKoHt3EhEREZF2OLuTiIiISIGYpBEREREpEJM0IiIiIgVikkZERESkQEzSiIiIiBSISRoRERGRAjFJIyIiIlIgJmlERERECsQkjYiIiEiBmKQRERERKRCTNCIiIiIFYpJGREREpEBM0oiIiIgUiEkaERERkQIxSSMiIiJSICZpRERERArEJI2IiIhIgZikERERESkQkzQiIiIiBWKSRkRERKRATNKIiIiIFIhJGhEREZECMUkjIiIiUiAmaUREREQKxCSNiIiISIH+D1yN8mQaal1+AAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 437.5x350 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "ONE_COL_WIDTH = 3.5\n",
    "fig, ax = plt.subplots(figsize=(0.5*ONE_COL_WIDTH, 0.4*ONE_COL_WIDTH), dpi = 250)\n",
    "\n",
    "labels = [r'$k_1$', r'$k_2$', r'$k_3$', r'$k_4$' ]            # E3 > E2 > E1 > E0\n",
    "colors        = ['C0', 'C2', 'C1', 'C3']\n",
    "markers       = ['.', 's', '^', 'D']\n",
    "marker_sizes  = [4, 1.5, 1.5, 1.5]\n",
    "\n",
    "plot_order = [0, 1, 2, 3]   # last index is plotted on top\n",
    "\n",
    "data_list = [eng_in_phase_k1, eng_in_phase_k2, eng_in_phase_k3, eng_in_phase_k4]\n",
    "\n",
    "for ii in plot_order:\n",
    "    d = data_list[ii]\n",
    "    x = array(range(4-ii)) + (ii+1)\n",
    "\n",
    "    ax.plot(x, d,\n",
    "            linestyle='None',\n",
    "            marker=markers[ii],\n",
    "            markersize=marker_sizes[ii],\n",
    "            color=colors[ii],\n",
    "            label=labels[ii])\n",
    "    \n",
    "ylabel(r'Engineered decay ($2\\pi\\times$kHz)')\n",
    "xlabel(r'Lattice size ($N_\\mathrm{q}$)')\n",
    "xticks([1, 2, 3, 4])\n",
    "ylim(0,80)\n",
    "xlim(0.5,4.5)\n",
    "\n",
    "title(r'In phase ($\\theta=0$)', fontsize=7)\n",
    "\n",
    "ax.legend(fontsize=6,loc='center left', bbox_to_anchor=(1.0, 0.5), ncol=1)\n",
    "# ax.legend(fontsize=6,loc='upper center', ncol=4)\n",
    "\n",
    "# fig.savefig(r\"figure_out\\scaling_eng_in_phase.pdf\", transparent=True)\n",
    "\n",
    "show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ba0d7f53-b066-4ec2-915a-5a17d9c6da64",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.19"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
