{
"cells": [
{
"cell_type": "markdown",
"id": "4b1b36d9",
"metadata": {},
"source": [
"# Wave turbulence in gongs"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "2cc6b20c",
"metadata": {},
"outputs": [],
"source": [
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"from scipy import signal as sig\n",
"from scipy import fft as fft\n",
"#import scipy as sc\n",
"import scipy.io.wavfile as wav\n",
"#%matplotlib inline\n",
"#%matplotlib notebook \n",
"#for jupiter notebook\n",
"%matplotlib widget"
]
},
{
"cell_type": "markdown",
"id": "9796e975",
"metadata": {},
"source": [
"# Acoustic measurement\n",
"\n",
"A microphone is available to record the sound emitted by the gong. Use the Audacity sofware to record the sound. Save the data in .wav file format that python can read using the scipy library.\n",
"\n",
"A quasi-stationary forcing can be achieved by hitting reapeatedly the gong with constant impact strength.\n",
"\n",
"#### 1 Perform a series of experiments with different forcing strengths (and possibly different gongs). Document precisely for dataset. Make long records."
]
},
{
"cell_type": "markdown",
"id": "c4feddd1",
"metadata": {},
"source": [
" "
]
},
{
"cell_type": "markdown",
"id": "d669b204",
"metadata": {},
"source": [
"#### 2. load one audio record (use the function `wav.read` from `scipy` to read .wav files). There are two channels but we will use only one (whichever). Plot the signal"
]
},
{
"cell_type": "code",
"execution_count": 17,
"id": "febe1854-031d-4cc8-8da8-35fadc046a8b",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"D:\\M2TMA\\20260915\\TP Gong 150926\\stationar regime small gong\n"
]
}
],
"source": [
"cd D:\\M2TMA\\20260915\\TP Gong 150926\\stationar regime small gong"
]
},
{
"cell_type": "code",
"execution_count": 26,
"id": "259a631a",
"metadata": {},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "979ad48439514f1aa14025857dda6af0",
"version_major": 2,
"version_minor": 0
},
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjExLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvctoD+AAAAAlwSFlzAAAPYQAAD2EBqD+naQAATA1JREFUeJzt3Qd4VFX6x/GXGnov0nuRKkjvTZSq2NZFRNFVsbHK/lkFV0VFwd5dCyiuoq51RUB6kSJdpPdeQwmEmlDm/5yLM0yGmcy9k7l37sz5fp4nJJncSS43k8wv55z3Pdk8Ho9HAAAAoI3ssT4BAAAAOIsACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACaIQACAABohgAIAACgGQIgAACAZgiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACaIQACAABohgAIAACgGQIgAACAZgiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACaIQACAABohgAIAACgGQIgAACAZgiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACaIQACAABohgAIAACgGQIgAACAZgiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACayRnrE4hnFy5ckL1790rBggUlW7ZssT4dAABggsfjkePHj0vZsmUle3Y9x8IIgFmgwl+FChWi990AAACO2bVrl5QvX17LK04AzAI18ud9ABUqVCha3xMAAGCj1NRUYwDH+zyuIwJgFninfVX4IwACABBfsmm8fEvPiW8AAACNEQABAAA0QwAEAADQDAEQAABAMwRAAAAAzRAAAQAANEMABAAA0AwBEAAAQDMEQAAAAM0QAAEAADRDAAQAANAMARAAAEAzBEBoZ8vBEzJy0jo5dCIt1qcCAEBM5IzNlwVip/tbcyXt3AXZcOC4jB3QjG8FAEA7jABCOyr8KSt2HTV9n/MXPLI75ZSNZwUAgHMYAQRMqDZskvG6V8Oy8s5fG3HNAABxjRFAaO3Y6bOWjv/5j72mj/19Z4qs3nMsw20ej0fSzp2XeKTO/cIFT6xPAwAQBQRAxD0VTCJx9NRZafjsVBk9d6ul+20/dFKe/ml1plPCqWfOSp/3F0jPd+YZ08deg7/5Q2r9a7LsOhJ/08kPjlsu7V+dJWfOxmeABQBcQgBEXNt68IS0GDlDPpm3LeLPMWLiOkvH3/bRQvnPbzvk7rFLQh5z5ES6723/APjj73uM158t2C7x5pfV+2XXkdPy68aDsT4VAEAWEQCRqbPnL8iirYddO205/Oe1ciA1TZ6bsFZ6vTMvS4UaM9YdMKZtw9mfesZ4vfHACTmdHv66pJ+/WHRidjTTPzCaMXP9Adl44Ljv/m9M2yiz1idb+hwAAL0QADWjAoKVKdNnf14jf/looQz9fpW4kf+atFV7jsnw8Wsi+jw7D5+Sez5bakzbWnHl05PDHjNu4Q7Tn6/K0ElGwcnm5IuBLpzlO1Pk7rFLpesbvxohfcqaA/LWjE0yIJPRSQAACIAaUWHp5g9+sxQOvli403j9w59Tl2a/zh1jFsnQH5wPjcnH04wQuGT7EUv323P0tG3ndPC4uYbTJ9LO+d5Wgc6MNXtTfW8v3Z4ie03+P9Rxv6zaZ7qoQ40EAwASBwFQIyv3HJNlO1Jk9oaDtj6hq68zd9Mh+WrxxfBolpr6zGqV6crdx2Tsgu1yywe/Wbrf0VOX1uxFg///YvS8baYKTSat2ud7W12GlbuPyoBPFxvTu+rarNuXGrUq3FajZsoD45bLA+OWmTp+2toDEX2diSv3yecWRkABAM6gD6BG3py+0fd2NhHZcfik8cR+e/NKkjd3jrD3X7bjiFxdqVjY485fsB4uVcDp9NpsScqZXaY82k6yZVNnGJ7Jw8JSYchOqtDkb22rWrpP73fnG6/X7TsuXeqU8o3GfjuwpTStfPH7kOZXkasC94SVl0KkGWrKeMSEtfJkjyszvebpfzbPtuLIyXR56MuL17Vl1WJSvVRBy58DAGAPRgA1Erj3bafX5hjB5LWpG0zd/6Z/ZxxVC7WXbiRdWfYdOy07Dp8yCitUpexHv24xwoN/QcSB1DPGtKXVIgmnqOvx25bDEbelyazoxBv+FP/RzTF+1c+ZhT//kcPAUUQ1QmklOM5Yl2xMsavRvcz8tOLSsoHDflXRAIDYIwBqzBukFltcL6eoKc0mI6bLe7M2hz32j11HZeSkdXLSb41bZlSvvBcnrTcChqrM9a6ja/7iDGOkbtyiS1OK4frpZTWMWbl/25dmyV8/XhjxdGkk9h27WJEcjArMszYky4SVe6X+8CnG20qw/5FaGmDWf5fuMkKoCuhmi1XOnndnaAcAXREAkaXeea9MCT96eP178+XDX7ca7UmsOv3nFGfTF6b7bvPvQ7f9cOgAeO78Benx9jy5Z+wSSTlpfQRKFZOoHoNm7+s919ku6ZOnAvOAT5fIw1/+LifTzxtvR9vUP8Pugi2HLgvj/tm535hFcird3B8AAAD7EQBhiZqaDcd/rGf2n6NOysbkExFdbTXtGwnVomXtvlSZsT5ZGj0/TbYdOmnp/qqYRPUYDFfEMPKXdfLuzE1hRw79bzsXZBo72oUowezPZMQwEi9P3iD/XbJT+n68SNq+PCvTY5dsNz/KCACwFwFQU2qqMhJqajYc/+xzl9+o0+87UqTvxwtl/f5LrUvsLNAInHb0r7LNzIdztvimS5VwdSYfztkqr07NOLo5M0gjZjUK57UzyNT1+7PDh+usUlvYBaMaWqtdVSLxeIgekYER1+wSAACA/QiAmvIfjVGtU5q9MF3mRGnq8uc/9ga9/XjaOVmw5bDcMWax6c91/My5oLc9/t1Kuf7deWKHkb+sz/J0qWpKHeiRPytivQHxqf9lDGNmdhXxH03cHMGI6uEQ09lqXZ8qCnpx0jrbqq3VXsLxuAcyACQiAiB8DZTv/MR8MMvMt8t2RaUxsvKvgJCkLNp2xAgsf+y+PGS5QaisNGtDxoCtppbVDiTJx61Py365eKd0eX2ORNtHv26VhVsPi13+Z6GhOADAPgRAXNbrTxVPWLE6yGhXOCq8qEbH8WCZif2BI9XulVnS7IUZlq/5B3Psmy622sA7lGDrIF8zUQik7hethtcAgOAIgBox09FE9fozs87PX2BxhZmvo6Yv7/bbku6ci9uEqJ1TrLLafUZVELv3CjhLtZdRTcHP+DW5BgBEFwEQl/lk/qXmwsEELuaPtDHzoRPpvlEi1SbGbqplzZQ1+8WNth9yz9o4tY2f1YrpaJq0ar/R3ke14QFU4VLgelkAWUcAhGWTV2cMUa9H0N/PS+3+ofyyOrJWL1bd/7m5vW+zIpKCiV7vzotouzU7qC3cOr46O9anIV8v2WXsOOI2qp/lY/9dYalJuOqByIimdao10n9+22Gsl42kl6dZ6nvplp8/wCkEQFj2j2//yPB+YEsTK+OBwfrhIbysbHBy/MxZ115itXuJv7V7rbUMcsJbMzYZ2xUu32luDasKFnWeniINhk+1fW3jrPXJxt7OVteUupX/7MKFKG+xGLj7kOobqrakNEt9L/t/sliG/rDStvPShbqWJ2gT5TgCIGKqz3vzZXeKe6Y/48XuFPNPVIG2OjS9G+r5+puloavE42nP4LRz5tYo7j168XuVfv6C0QrJTgPGLjH2dm738qzLwjRCU4Fehc2vFpkvgFq995ixK9FXizPvehCMGg1mRPiSOz5ZJPWemUKbKIcRABGU2uPVyhTXoRPmW7v4U0+Iqim1jX/cI0J2jVb98zt3jJgkp56Re/+zVOZuOhjxaJtV0Wq1FM7eY2cyNGFH9AXrUWqGGp2t/dRkY3/uRBmpPXb6rNz87wXyn9+2W7qf+uNfbQwwf/PF1lPfL99t0xkiGAIgQv5AT1lzcZ9Xs1Mokdp15LTx9eAMs0sUv4vgl/HM9Qd8W/d5XF7XrHpMTlt7wFJjcvX/81IjbVat2OVc66N1+9w3fe4UNe1u5Q/YSHw6f3tE1/rtGZt8OxW59feeKvR74vuVGfZdD7dF6NIdKfL0T9YKt9q8NEuue3NuhGeJrCIAIqQ/LPTpU78ovFN7uvdwm7fpkIwPsRtKrC3bkWJ63eVvW6w3hL577FJj6z5VSOJ2+yLYF1n9/6LxPbBDpFv5udn/Vlj/OTp8Is2YTmz6wgx5beoG01P13pHTSNbSdnvLfIh5e+Zm39tuXff2zszNRhGWWuNoxsm0847uaY7oIAAiJKtBzju1d17z+VzVwsTN/i+giCczZ89fMIoKpq9LthxGPvrV+giZk4Jt1+eEm/69wLb1f4nm+QlrLd/npxV7jfWWalmKCjL+I3XhfLdsd4ZRXrtd++av4jZq5DQrjeavem5aVM8H9slp4+dGHHth4jrTVY6BNM9/rudtvWPGf5fsimiqU02rqubWmT3JZDPRL8fuaTxdv6862WJxz+xP5m2XTrVLO9LB4MxZ82sAJ6/eJ6fSz8uNjcuLnX6zcSvIcAJ7zMJejAAiqEjDHxKLt4LVqszCnzL0h1USbZE2JLdbsLNyqtG2me0W1Siv1YpUda2PnXLn+rXnIhg19Ddv8yFTrZKiNZVv5o8cNRsz8IvlxlrrBVsOyRcLd9hWQOI/neu0j+duk9Pp7ADkFAKgRrbHcHcHxB/VhPn92fbsOazWF0WTWoNa5+nJ8oODVYRZGZwM12hbBSxVTKMqlbOi97vzwx6jWsaoHnhWQuCtH/4mDZ+bKrM2WK+EjgdWCoO8xkQwUq6oRtfh+PdA7PvxIqOA6aso/wy5ZdRdPRb3R7A+F9YRADVykr+s4FCvwUjN25yx6nD4z2tloYkpKbWvdNq5C1mqRneTz3/bbhTTdH5tjiPFMOo5f+OB46bv4x39GvDpEsttdNSTe6dXZ8vouda2f1Sh2KnRoUiqtdV6xdkRBOKPItwG8w+bKspVT8RY+9+K2J+DDgiAAFwj2BPvbR8ttPQ53NZgd/th6yPvszZcDFV2N46OhlenbLB2/NQNRjPyERPXWbpfsxdnGKNDVnvNOSkRei/+ErDVZ6KO8sLhAHj69GmZMWOGfPHFF7Jo0aKQx61evVo+++wz+fnnn437xPoYAPHjvVmX2mw4QU2Z7Tl6OuTU2VvTL/Z9ixXbA7HFza83+RVlqPWHVlntNadaMlkZ3XSjMzHcp1iN8lr1t8+WRvS9RYIGwK+++krq1KkjL774okyePFluuOEGad++vZw4kbFC68knn5SWLVvK+PHjZdiwYVKvXj3ZuXNnzI4B4JxIlx/5V2Sq1h9OB87Wo2bKGzEOeqG47Yn4/IVL52MtOkZGLQ3o+ob72q1Y8ca0jUFb1sxYl3nLmlW7j8k9Y5fY1nsylOnrDsikPxvCw70cC4DFixeXxYsX+0YA16xZIxs3bpSRI0f6jlm4cKEREFUg+/7772X58uVSqlQp+fvf/x6TYxD/TzbQtxm3U16dujHDDg+B3FmbHBuLtx2R1XsSe4eSyk9MlOHjrY1ShjN/c/DH8z2fZd6YvNe782TG+mSj96TTTafDjTwnyjZ48cyxANi1a1cpWbKk7/1ixYpJ06ZNZd26S+tAvv76a7nyyiulY8eOxvu5cuWS+++/XyZMmOAbKXTyGFj3rsOjL4hfT/4Y/VYwXv3GhF5iEqkUB3c38a/6dJsDWahMVtXDOhi7wH3rFA8ej2y/dmXD/syn0IMVaoXrS/nzSnfulqSTmBWBHD16VObPny+NGjXy3aZGBdU0sT/1/rlz54zRQqePCZSWliapqakZXpDRBH6oYdK4RZcvtwiVe5weLVCjF4E74fivXYsGNYUXapuyFX59OPcds29t8s9+WxZmMzkhG9isNyvTuG6IuU5Pj0YiNYI9g7MS+AKt35/5c936IAExXAupoy7tI6mTmATA8+fPS//+/aVw4cIyaNAg3+0qUBUpUiTDsUWLFvV9zOljAqnpanXO3pcKFSpE9P8H4N6gcOz0Wan91GTp/d48U/vOhrMjRBWw2pLvnRnBR8z9q39bjpwpu1NO2dKf7WOLrViUhVuPXBaWzez9HCzszvmz2jmW1LV1Oyt7FHuxqwZcFwAvXLgg99xzjyxZssQoBlFByitv3rxy/HjGvyS8YSxfvnyOHxNo6NChcuzYMd/Lrl32NOIEdOWJYtRTm9JHsjG9t6+dWqsWbieUGSb2SM5spGP2RnMtNtq8NEuG/bhaIlF/+FT5fOGOoAFy5W7r+yGv25d62ehP4+enGXvvZqbFizMuu+142lktdoPJjKogd8K4heEbTruFE8VBcDgAql9A9957r0yZMkVmzZolNWvWzPDxGjVqyLZtGbupe9+vVq2a48cESkpKkkKFCmV4AWC/Gk/+YnnrMbUpvXoxW5ikpgIfHLcsQ+hrNWqmbD3ojjXBXy2OvEPBU/9bHTaszjEZRkNZHmYqNSWKU36hetMt2ha6aXif9+dLl9fnSKfXZpv6wyBwCYCdnFjiEMme3nZy8TJXbWR3OvxNmjRJZs6cKbVr177smN69e8vSpUszrMEbN26ctG7d2qgidvoYWLflINvNwR73fb40Q0HG5uTjUV0/pSolJ63aLy9OWp/h9jkbYz9NGQ3h1jBGe3s+O6lAa3UE8PedR2Vz8gnZevCkqanvN6cHXwtu1iYLvQeDraFLBKpNTahwGziSDOfldOoLPfXUUzJmzBijzcqcOXOMF6V06dLSp08f4+2ePXsaL9dee60xTbxq1SqZOnWqMVro5eQxANwzKnDQb4qx0fPTjNfTB7eT6qUKOnFqCc/Orf9CBTOzhSfRZmbA7e2Zm6X3VeWkeqkCEX2Nrm/+Kiue7moEoOIFksQu//t9j9zQqJy4kbdNzQ8PtpLGFS+usfc65bIde3Tk2Ahg+fLljTYrZ86ckRUrVvheNm261DsrW7Zs8uOPP8rzzz8vhw4dkrp16xrBrFmzZjE5BoB7pJ29cNloQt+PF5kq6rBrwXxWR4ncZuXuoxFPgX4wJ3TV56fzg08/un2bscymlM38MdPw2aly9YjpYdf53TFmUcQ7tjz63xW2T1dndWbndhM/p0jgEcCBAweaOi5HjhzSr18/48UNxwCwj5qSMzvCop5EO7w6W+Y93sl3W7KJVhedXpsj20f1yPSYus9MEaeqM91QYKNG5IKtjXz06xXyr55XGm8//OXv8totDaVb/TIhP4+/5TuPGvsCn0w/J8/0qpvhY6H2/f1pxV555eaGkjvn5WMRdlQ9x4raKWbbyO7GwEOo9ZGvT9sow7pfvPaR7ISTO7t9o6mqyfngazKu2bfidLBwmzjf3rgVsz6AAHDdmxm36Np55FTMpilDSaAc4tP9rblGq5tAWw+dlLvHLjVeTqWflwfGLbf0ed+dtVk+nb89bPV0YD/EYPsl74ujcG3G2jBr3j761XpLnihu0QwNEQABxIz/Hr6RLoZPPh6/QWHv0dic+wYLBQofRxBMrGwJOezHVfLG9I3yx65L08+KG3L32XPRq849d97Z/9GB1Og1gnbayF/Wyz+/+yPWp5HwCIAA4lqzFy7vLxcvzDRQNrvOzi4vTAo+fRttZtZqeqn1cgs2H5K7Pl0su8KMGmfF8J/XSrwa+sPKy25TlfOqqj3S9YZO+mbp7lifQsIjAAJAGOnnLoTcts1Jo35ZL4u3HZFxi3YY4cepJ3JVfKO+9q8OtsTJbA3goRPp0nf0Ipm94aD8/evfRXfB9o4OVrjR5fVf5c5PFhu70CRS03dEhgAIIK5MXLnP0a+nRqZq/usXuWPMYnEDNTX25I+rjfDz/qzgW8lFey3jt8t2G6OP/T8xdw3en7VFlu90Zo/deJnqVCEtXGHLTyv2RPS5X5u6wdLxExz+GYqUnaO7IAACiDMPfWmtMCGr3ppxeYFCKE6MEm4/fOlJce0+ZxoI77FYfPPfpbvkxvcX2HY+8eje/yyVk+mZPz7eDFIMY8bHc+3f5eObGDQK/8uHvzn+NXXCCCAARMkT36/iWkbYyiUwaJu9q1N76WaVmraO59Y2//z+8jWFip3/p3hqsxSPCIAAtBh9cWJLtx9/j2wK7/CJNOn1zjz5fOEOsUvqmejtxWuWmT6N/nsxR2rZjiOWjje7jWC0nXW4EtgJb0yLrBl6HGfhhEEABJDwpq09YCx+d4JqgfLy5PWyYMsh0/dRu0Ws2nMs5B630fDhnK1RWbBv5YnbqVZ0d32yxNLxx89c2vklVHPmSKqTw1m0NfOdRbYdOin/+OYPOZUe2c40VkRr5E5tmYf4RAAEEFNOjkx51+jZOW1V48lf5P3ZW0xtUxcNE1buNb0lnQqYkRQa7PPrV3g6zDq2WDge4VZ+Ztz6wW+ONmf+fvlu+WC2/e1+qgydZLkNERILARBATL1noZI1q9QuFcqKgKbD8Wr6ugPGlm1miwfUFPPfv14hbnXs1KU/BsbMs17YcDLtfEybZkfL/lTra98iCeaj50Z395HMqO0HI6nqTSGk2oYACCCmjp50bgTQu9XcmbPR2+Ehqyo/MVHc7mzAji1mvTJlg6XK6Kd+Wi2fzNsm7V+ZJWMXXAzrVny9ZKfpkbjVe45Jnacny6CvzPURVH80bIpBGDTryqcnW1p2oDi5DG/gF8uk7cuzfG2czA7C3/D+fHtPTGMEQAAxpVqG3D56oSPrnhCZn//YG1FoUP0Dh49fY/r48X/slecmrJUdfq1urDbsNmPh1iPS8515xn7HZhw8niY3vDdfrnnjV2OdXqTsLnxQyw7sbA7+4LhlEY8aqnW4yuh51u4f6WMB4REAAcTc/M2H5bMF9lXAeiXHSdNgN7O6fvKrxc73j4u2vX6tZjq+Olsmr94vbvXo1ytkyXZzVdHbDp609P2ctGq/jJiYta0B1ej77ztT5DxlwDFHAATgCqcdGAFUa+YmrYqPXRDgXmorvkh8EeH9rJi8Zr/cYrJwRR37RoTNpyO1bl+q9Hl/gW9EELFDAASglbdnbGIfUsRspNuMb5buFid/HpTfthyWNXuPOfZ1EXs5Y30CAABk1fYsrM1zo33HThsFKne1quLI1/rrxwst32/LwRO2nA+cQQAE4A5mmqRpwM5F/NHg1g0cOrw6WxLJMz+tkSXbU4wXt+317HX4BH0E4xlTwAD049YUIyKfzLfe/w6JJ5o7kADBEAABuMIvDhVnqOLDWOyLG8muG4CbjbHY0gXuQgAE4Aqbkp1ZT6R2dqADReR0vHaLt5lrq6KbKWuo5I1nBEAAgKVmzbp5YVLWet8lIhq3xz8CIAAALrPI5aOOdZ6eEutTQBYRAAFo54Fxyx35Oqq3GuB2Xy4yt4cyEgsBEABsEklvNcBpP/y+x9LxyccpVEoEBEAAACy4oGMljJ9XJm+I9SkgCgiAAADYsKVbovp2mXNb1cE+BEAAcJHzmo8uAXAGARAAXGT2+uRYnwIADRAAAcBF9h5jgT3g78IFRsXtQAAEAACu9d1y1hzagQAIAABca+FWvYtu7EIABAAA7sUMsC0IgAAAAJohAAIAAGiGAAgAAFyLGWB7EAABAAA0QwAEAACu5WF3HFsQAAEAADRDAAQAAK7FGkB7EAABAAA0QwAEAACudeRkeqxPISERAAEAgGvN3XQo1qeQkAiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACaIQACAABoxvEAeOHCBVm9erXs3Lkz5DHp6emyYcMGSU5OdsUxAAAAicSxAHjq1Cl56aWXpEaNGtKsWTMZNmxY0OO+//57KVOmjHTu3FkqVqwoffr0Me4bq2MAAAASjWMBcM+ePXLkyBGZNm2atGvXLugxW7dulb59+8qzzz4ru3fvlh07dsiKFSvk8ccfj8kxAAAAicixAKhG/tQIYNWqVUMeM3bsWClevLg89NBDxvulS5eWhx9+2LhdTdU6fQwAAEAiclURyLJly6Rp06aSLVs2320tW7aUEydOyKZNmxw/BgAAIBG5KgAePnzYGJXz533/0KFDjh8TKC0tTVJTUzO8AAAAxBtXBcBcuXIZIcuf9331MaePCTRy5EgpXLiw76VChQpZ+N8CAADEhqsCoApU+/bty3Db3r17fR9z+phAQ4cOlWPHjvledu3alYX/LQAAQGy4KgB26tRJFixYIEePHvXdNmHCBKlevbovlDl5TKCkpCQpVKhQhhcAAIB442gA/OOPP4xWK8ePH5eUlBTj7bVr1/o+3q9fP6lSpYrcdNNNMn36dHnjjTfkww8/lBEjRsTkGAAAgESUzePxeJz6Yk2aNJFz585luE21X5kyZYrv/QMHDsjzzz8vS5culWLFisl9990nN9xwQ4b7OHlMZlQRiFoLqKaD42E0sPITE2N9CgAAWLZ9VI+oXrXUOHv+jvsAmGji7QFEAAQAxCMCYIKvAQQAAID9CIAAAACaIQACAABohgAIAACgGQIgAACAZgiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACaIQACAABohgAIAACgGQIgAACAZgiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACaIQACAABohgAIAACgGQIgAACAZgiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACaIQACAABohgAIAACgGQIgAACAZgiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACaIQACAABohgAIAACgGQIgAACAZgiAAAAAmtE6AK5cuVI6duwoBQsWlEqVKskLL7wgHo8n1qcFAABgK20D4JEjR6RLly5So0YN2bhxo4wePVpefvllee2112J9agAAALbKKZoaO3aspKWlybvvviu5c+eWMmXKyKBBg4wAOHjwYMmePbbZ2DsSmS1bNt/72w6dlGL5c8u/52yRD+dslUmD2krtKwoaH085lS7r9x+X/Ek55an/rZYi+XLJc9fXk0rF8smQ71bKn58GAABA3wA4f/58adWqlRH+vDp37iwjRoyQrVu3SvXq1WN2bit3H5Xe784Pe1z3t+dm+vGOr86O4lkBAIBEoe0U8L59+6RUqVIZbitZsqTxev/+/UHvo0YMU1NTM7zYwUz4AwAAiJS2ATAY77RvqEKQkSNHSuHChX0vFSpUcPgMAQAAsk7bAFi6dGk5ePBghtuSk5N9Hwtm6NChcuzYMd/Lrl27bDm3bSO72/J5AQAAtF4DqNb/qbYvZ8+elVy5chm3zZw505gWDrX+LykpyXixmyr82D6qR8iPqxHKlFNnjYIQKyo/MTEKZwcAAOKdtiOAAwYMkBw5chgVvykpKTJ37lx5++235e9//3vMK4DNBESr4Q8AAEB0HwEsUaKETJkyRR5++GFj1K9o0aLy0EMPyRNPPBHrUwMAALCVtgFQadKkiSxcuDDWpwEAAOAod891AgAAIOoIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACaIQACAABohgAIAACgGQIgAACAZgiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACaIQACAABohgAIAACgGQIgAACAZgiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACaIQACAABohgAIAACgGQIgAACAZgiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACaIQACAABohgAIAACgGQIgAACAZgiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgAAKAZAiAAAIBmCIAAAACacTQAJicny6hRo6RNmzYyfPjwoMccOXJEhgwZIh07dpSbbrpJJk+eHNNjAAAAEo1jAXDDhg1y9dVXS0pKiqSlpcnmzZsvOyY9PV06dOggixYtksGDB0uTJk2kZ8+e8uOPP8bkGAAAgESU06kvVKlSJdmyZYvkzp1brrvuuqDHjBs3zgiK+/btk2LFikmvXr1kx44dMmzYMOnTp4/jxwAAACQix0YA8+TJY4S/zMyYMUNatWplBDKv66+/XtavXy979uxx/BgAAIBE5KoiEDUCV7Zs2Qy3ed9XH3P6mEBq6jo1NTXDCwAAgDZTwCtXrpT77rsv02P69+8vDz74oOnPefbsWUlKSspwW968eX0fc/qYQCNHjpRnn33W9P8HAABkTfmiF5+b4ZIAWKVKFXnzzTczPaZcuXKWPqeajlWVuf4OHz5svC5evLjjxwQaOnSoUTDipUYAK1SoYOn/CAAAzLuiUB4ul5sCYMGCBaVFixZRPZnGjRvLZ599Jh6PR7Jly2bcpqp08+fPLzVq1HD8mEBqxDBw1BAAACDeuGoNoJoy3r9/v4wZM8Y3Ivfuu+9Kv379fMHLyWMAAAASUTaPGgJziGoAfe7cOaPSNleuXFKtWjUpVaqUjB8/3neMas/ywAMPGLerFi1t27aVb7/91hhxjMUxmVFTwIULF5Zjx45JoUKFxO0qPzEx1qcAAIAlTSsXlW8HtorqVUuNs+fvuA+Aaoo18Mup0bZGjRpluO3UqVOyceNGY51exYoVg34uJ49JlAcQARAAEG8IgHHeCFpp3ry5qePy5csnV111lWuOAQAASCSuWgMIAAAA+xEAAQAANEMABAAAruVcpYJeCIAAAACaIQACAABohgAIAABcq3KJ/LE+hYREAAQAAK5VPH/uWJ9CQiIAAgAA98oW6xNITARAAAAAzRAAAQAANEMABAAA0AwBEAAAQDMEQAAA4F7sBGILAiAAAIBmCIAAAACaIQACAAD3og+gLQiAAAAAmiEAAgAAaIYACAAAoBkCIAAAgGYIgAAAAJohAAIAAGiGAAgALtK6evFYnwLgLuwEYgsCIAC4SMVi+WJ9CgA0QAAEAADuRSNoWxAAAQCAJXXKFOKKxTkCIAAAsCQbo3JxjwAIAC6SL3fOWJ8CNFOlRP5YnwJigAAIAC5y89XlJVE0qVRUfnqodaxPAyEM6lxDRvdvIjMGt+caaYg/NQHARcoVzSuJNE3YsEKRWJ8GQihdKEm61CnN9dEUI4AAEt6aZ6/N8P7Ux9qJWxXKkyvWpwCX6NeiYqxPAQmMAAgg4eVPyjjZkTuHM7/6qpZkbRUiN6z7lfJEt9pcQtiCAAgANhn/cBuurUNeubmBZHeoMrVycfubdb/Yp75REDSwfTXbvxb0RAAEAJvQKcM5tzSpINkc6E3y1b0tZOpj9hdNdK3r7rV5VxTK49jXurt1Fce+lk4IgACAhGB3/JvwSBtpWa245M6ZPeH+eOhRv4zpY6+pU1pevLG+OKW0g2FTJwRAALBJou1hf4tLW9TkiGDu95ledSzfp165wpKorAyefty/iaVQ1qxyschOCrYiAAJIaLc3v1hJ+VBH59dS5c+dQxJtmtWNrMa/l26qLwMsTiv+zHrOqNk+qkf0PhkiRgAEkNBe6OPcVFUgJ9akOeGeNlWkY62SRmNnK66te4Wl43s0uDQN2b5mSdP36/tnyLfrcq8a3lXql0+80b+G5S/1aGxdvURMzwXOIwAC0EIuh1q/JGLft6d61pFPBzST7BanWu9qVdnS8QVy55S2NUpIiQJJcl+7qqbu065mSXmyx5USaaWtGQXjrDdj3bKFTB1XrWQB39u3Nqkg79/eOOx97mhRyff2kGtrSff65kJ+uSIXG5zTGsk92AkEgBbUlN+kVfukZ4Oyjny9d/s2snyfpJyJNWWcM4LQ/Z+7m8kFj8iyHSmmju9Qs6TvumUzJoM9lkYOa5cpKDe+vyDkMQ92iK82LGMHNJXGlYpKg+FTwx6b12+JglpH2d1EIUg+v/s81LG68bryExPD3u/Le5vL6LnbTAd72I8ACEALhfPm8rXv2HH4pO1fz6mg6Vaf3tU04mnzHBYGGjOMKJm838WgGF7XOqXln9fFVyPmDrVKmR7VdVKl4vnl+RvqOfo1kTn3z4kAAOJOx9rmgkhWvNCnnqW1glZlj/Kiwq/vayFuEeupWNVKJtwUM+xFAAQQc/e2rSLd6lkrGIBIsfy5ZflT18h7fcOv3YoHV/sVmXg84adyb29eKUOhTZ0y5ta+xYrdJUGtqhWPWXCzKmeI9aRq/SecQQAEEHN9GpWXf/e7Ouqft3HFS1WO8cxb5Rro8etqGSHQv3o2nt2cxT6Dbh89srOY5J/X1ZJ3HfhDoInNPf3qlXN3iE8kBEDEnFubyyL+fRWlKbfrryorA1pbq2g1WyX73cCWl1VKBipbOHjT3Ty5rBWNVIvxtF+4rcWsVhkHymll8eCfnGrUo0bO6piszo3Egx2qG38MxLuKxfJJzdKXqpNhHwIggIQVrarat25rJAPbR78aVE2DqREVVfmqeu31s3kEq0i++A8I8UqN1sKcL+91z1rJREYABAATTCxJs6xZlWK+XnaqKtOJPWbtEmqUMlL+LUpCFYBYUat0wahNvUfCv+ceMqfWAQ7qdLHFDOxDGxgAiJFoLai3c7ozyWQozequJ//qmbGZc/1M9t29sVE5owDETnP/2VEqFMsXtc9n9foUyuu+5tMjb3RuV50CeYgnduMKI6hGFYsYDT/nbz7MFYLtckWwdisRxMNWcRtGdLN8n+eur5vlvomZXZtI1vpZFc3wZ8XvT13j2p1r/hLFvaDDPfTvaFFZfttyWK6pQ3cAuxAAEdRjXWoaC5abjJjOFYLtqpdKjOkx77ZYd7asbGqxe6Lqb+L/j+CKOlTIYaVlTulCSUYFs5N/r6glAGr7QdiHAAgg5uJhJMyMG64qJ13rXmF6yy6zShWKzvq66+pdIUt3pBhP6AdS00Qn/usri+Rz3/Sq00pbeEzNf7yT8TMayc+p2mLO6V1HYA4BECEVSHLm4ZEIrQvgrNw5sstNV5czdoGYs/GgfLV4l21fS7VqifaavivLFJKqFooC1Jq3dftSpUXV4nL/58skK/8XNfKo9oo1M7p/u4UiiCol8sueo6fFre5vX1XenblZNh44Ll2udHbtpdv6E876vw627+msTH60rVQvWSDi+8NeBEAEpf7Qs9pjzO1BE4ljYIdqMviamsbb19UrY1sAVFt3+e9OkZlSBZNsG8lUT6DP9Lq4ru6RTtXlnZmbwxZKhPo83hFKdb7JxzMfBcxr4XdA5RL55MEOzaV4BDs5dLZ52zj1+dUaQ6f3Z76vXVU5kHrGVXvgFs+f2wjrTqh9BU2d3czRZ97t27fLxIkTZd++fVKjRg259dZbJW/ejI1P1fY/P/zwgyxdulSKFSsmt912m1SoUCFmx+hKjTRYpYb51TUdMXGdLeeExFQ1giejBztEvydfoL+1qRLRz4Hd/tG1ltGe5NDxdEujiGPubJLh/S/+1ly6vvFrVM+tVfUSIa/l6Hnbgn5szpAOltdD5g/xR2OoAF6peGwaYA/rnrGy2Q5NKhW1tBNMviRn/rCH+zk2Lvv666/LtddeK2vXrpU8efLIe++9J1deeaURBv2pUPjoo49K9uzZZf78+VK3bl35/fffY3aMriKpQFO/ege0rmLL+SBxjX+kTaZ7BAfqWqd0lkenKxQNHziqmShM+fHBVhF9/asqZG2LujKF80r98tZG/zoHTHvWjKAvXmYalg/9f3qyR+ggpMKf1ZHTQZ1qWDo+kX33QCt+78LdI4DXXHONDBo0SHLmvPglhwwZYgSuUaNGyVtvvWXcNmnSJPn+++9l9erVUqfOxUWjPXv2NELanDlzHD8GEtGCXyBaSwDua1dNPp4bfOQoKwK3HHu2d115YdI6ST93wXdbjYAA6JHLO0E3qmhuejhQ7jhpe9P7qvBTptMHt5flO1Pkpsaht3QMFfD+2qxiRNPmTlXKOlmg4v/Ys9P1Dcs58nXgfo6NANavX98X/pSkpCSpVauW7Nmzx3fbTz/9JI0bN/YFMqV///4yd+5cOXLkiOPH6Kr2FdEdGQAiVbJgkjHd2b+l/Yvo29UokaVN78uG2Mc3q/4bpf2MrVrwRCdpkMmonn8Ln1ubVMjyPr7RUDRK1b12r0kMtGp4V+OPEKvFGd51sGZ93L+JDOrsjtHTbCHaktuw4Q5CiFlpzo4dO2T27NnSrl07320bN26UatUyru1R76t1ZZs3b3b8mEBpaWmSmpqa4SURhduCCXDSi33qy3PX27+I3n9h/IYR11m+/9+CTFeHEm7U6/+6Xnpiz5fbmYmawHV4dgVaO7UJsQYxFn5+OPTShmB7Vt/ZqrLl4gyrYU5VqcfzdoOIroh/s2zbts03dRtKly5djGnVQCdOnJCbbrpJ6tWrJwMHDvTdfurUKSlYMOPoU6FChXwfc/qYQCNHjpRnn31WEl2kf8cnSCs3aOj6q8oaFbJqGq5Hg7LGE7JZ/7yulqzYeVSuM9n/zwz/J+lgU89WPWCiaOa9vo2l17vzTI1I/uWjhRIrH95xdciPRTKdbOV7bbZAqWGFIpbXaLppBmj9/uNR+TxI0ACopnArV868P1bRopevj1HhqlevXpKeni4zZ86U3LkvreVQgezo0aMZjk9JSfF9zOljAg0dOlQGDx7se1+NAFI1DMS/Ivku/h56NoKRxgc7RH/T+r7NK8mLk9ZH5XP95+5m0q5mybDHVSqRT7rVu0J+Wb1fOmUyBdo8xpXR10YxaCtXloleUBnarbbc397+CvVItA1Y4hBKsyrFohIAVXEKEjQAli1b1iiYsOL06dNG+Dt06JAR/kqUyPiAVOvx1O3+1q9fLzly5JCaNWs6fkyw0KteADivRIHccuhEunQx0WxZjT6oxr/xKpq9Mc2EP6VQnlzy6i0NpVv9MtKxlrn7JIJo9m40e62d6vd3+GS67/13/trIsa+95Mku9HeNA44tBjhz5owR/pKTk43gVbLk5T8oqi3LmjVrZN68i9MQ586dk48//li6devmG5Vz8hgAF1ldN6QKCMy4v535kDb1sfbGaNbNmVSbek1+tJ30aRT+OFzeX693w7LGvq92965Tbm2SON8jNYWudnhxS5HP63+5Kugotx1yBhQAuaAeCG5qA6NGC2fMmCF9+/aVF154wXe7mkb2jiS2adNGHnvsMWPdoHpRI3IHDhzI0JbFyWNgTaQ/8znipCWGzlQj3/dnbzF9vNkCgqEWGuWqLQPdNMJipaJ0xvrkDLf1ami+cW+kKhTLK+1qZH69vA2arTQSzqpapQsaO6wcOZlueY/jrS92FycMy6RvYSiRXMNRN9a3bUpdbZPoFLVjzPG0c459PcRZAOzRo4fUrl37stvLlClzWcNoFRLV7hy9e/eW7t27S4ECBWJ2DOxVs3QBY5/Mlydv4FIniG/ubymJKNJp2ceuqZkhAKpRoqsrWWsxE4lfh3QMO735RLfaxtZwDSs4V7CgTkkV3JgJf+WK5M2wv7BTrWaqWdhhJStua2Z+r2WnZI9gSrxNjRLG2lErcvKHvz4BUE3/mtWkSRPjxS3HwB5qQ/bRAdtTwX3UFmLLdlwskDKjdCHn18l+eW9z6fvxoojuqwofzIh0WjTw+bSDxfV1JSLYW/fi1w3/RK6CmFr071YqmPoHwEgVzMN+42bliiCYXX9VOcsB0KmQjdBoCITL2LWhPaLLyvq5SF1bt/RlW4iF44lBJ9dW1az3fxv/cGu5p00VefO2jGuloq1ohGuv1DSpatzrxn58d/7ZmFs16Y4H99rws6JmLxJRgSTrf+hEMjDLs0zs8WcREKcK5bV3ob4/O4oCpjx6qQl8LKhdLszsdJFVkQa4FjFut5KZp3rWkT6Ny0u9suaLHtTU97p9qXJj48i2IlMFQFa927eR9GwQfju7SAy+ppYM/GKZOG1g+2qybMcRudumfdfvaVtFFm49LL9tPSx2oiF17DECiKhh5DBx5U+K/u4wtTRtFJsIIx9q6viqCkWM12b99/4W8vk9zeSeNpGNxkVSAJRZP8OsyhejHZPKFckj3w5sZbTrsWut61f3tZDW1e39A+T2FpVCFpTAGQRAOGb64NiO+ABukEPTHhmqz2DbGiUd/f+b3UZP7ebi9dz1dU3dp1455wpnYtEup1u9iwGzbGFrldpZKaq6rWkFqRCwJSHsQwCEJUWysNl6uSL8YEerIbLSp1FkU2nR2LA90O9PXWP7uSBzdfx60N1kolciMu7mooqd7mtXVW5vHnxkKlhbIqe8fHNDcVrfZhXls7ubycRBbR17bhl1UwNHvhYuYg0gLMmWhV+Iuox8DOpUXd6eudmWzz3hkTZSt2whOXP2guSN0RRUqAbCOlAVzgdS08SNVHXx2n2pxtsjb6wf69OJO6rYyWrBUyJTLXec7CUI5zECiMu0rmZ97Uf/lpWkR5g1KfG66Pe3oeZ2tfAa3LWWEdTsoKad1FpLN4U/b+sI1a5Erd9xY9WqHRW93z9gvt9hvxb2V8sO6lzD6Os37bF2cfuzpoObr45sdPauVpWlSon8cqMLR3cDd0CJ5PdTNLc/hDlccVzmoU7WNrdXe64O7Wauc37XOqVl6toDcXXVyxRO3EATLSqUqqB8weMxFT7sWlfkJCvNnEfcUF8WbT0im5JP2LbrRp5cOYwKUbh/Z5jvlu22fL/hveuKx+OxVGyn2jhNWXPAdK/LSKl1ez8/3EZW7z1mNO82u/ZS9aBcvO2I8fbTPevYeo64HH8m4jJJOa399RYY/gZfUzPksY0qXtwDNB6oBcnL/tXF0n2ql4qsN5iaPldfz05qBMGMlyNch5MrR3ZTj50G5QvLF39rLvFIhaxIqbVUi4d1ltpXWN8vVleNKrj798Wbf7nK+IPn07uaurLTwuu3XiXv395YXr3F/jWE9csXlr82q2ipWvs1v/NqXd16L09kDSOAiLrM1vqpxrsvTV5vvP1gR3ePVtQpW0iKW9yFIZJVjuufv06ScmaX0XO3ZXpcVvfBVWsHw1FTuLc0KS///H6l8b4dPcHHPxx6ejxPruzG+sZI/LVZBflq8S5bRzvUE+l9ny+VQZ1qWL6vCgpW977V3V2tK0tSruzSysKylC5XlpLp65KNNjV2u6FROenVsKxr1zertbndbWoXE62RQ/X7L/38BW3WEbsJVxyOUk+C20Z2l5Pp51nzYXFUKavPMUkmv04s+zn+vXNN3x8IkUyRXVevjDSrbN/WZmqEd+Y/Otj2+XH5qHL/luZGrr1eu/Uq+fmPvY4FH7eGv3j6/ZeVkXVEjilgjUS6p6i/SHf5evnmBhkCRjQW/Lpxf89KxfNZ3qnArGd6metPFkrHWiWNl8y274pVc1uvMllYG6imn1XVotsKZOCswnlzSb8WlRxt0wLEIwIgpFrJ/FKyYNbDYWZubRL99W0z/tE+7DHFA54EnjfZ5DVSL/7ZfsPsX7RNLYxWqQrArMiZPbt8OqCZPNalZsjw9UkU1zLVDrLTh+q1lpmeNhVIAAAyIgBqJfj43Yx/dJAaERYvxFKpgpmPFn19XwtZOKyz/PL3tnJNndIy+dG20sPCvqA1SlnbquzetlV85xRpMYgTgo2QqRHZ+Y93koZRXDdVreTl1yBcnzUrW4sBACLnvjk0OOqDflcbr0fcUE/u/HSxPNDeWguYO0Ls5xhr39zf0mgx4O1R9XH/iyNPR06mh7xPxWL5pHiB3MaC85qlC0pLi/0QndwaKhLe9izBpt+XPdXFaPwaqGrJrI06ZkUMlyICtqkd0DMPiBUCoEb8R37mDOkglYpfenKvWrKAzP2ntYbH0we3l6ompiVvbBx+y7L721WVD3/daurrfnTH1WFH2BpVtD6SVb5oXvny3haW76e2S1q2/Yj0MjG6qMLXibRzvvedWq+2+MnOmS5WD2zf8t3AljJlzX55uKP1atdQRpncneKLe5rLM+NXsy0UEpJayvHDg62kRH57l90A4RAANfLKzQ3lto8WGm/7h79ImZnmVFtnvXBD+Cf+f3StJb9tPSxta5SQW66uIA99uVzW7L24rVWgrnXDt/nIEWL4KBpFDh/0ayyPf79Kjp0+a7yvCg/Mbpm09F9d5FT6eVm6/Yicv+CRQnnM7a38+T3NQvb2G7tgu+Xp8lqlC8qGA8dDHt+kcjHjJVpWDe8qBU3+X9vUKGEsSwASVeM46oeKxEUA1EiLqsVleK86UtTB6rgxdzY1Ncqlpif9+8Op7Y7W7F2bobLPG7jMCDadGa44w2OixDln9mxGq5GDJ9Llqf+tNnUuz/auK3M3HZJnr6/ra3lgJsT6a1sjeMCMtHhnymPtpPITE00fn9XZWLPhDwDgDFZca+au1lXk+qvCT8n6TwUGVpKG0/nKUr63a5SOrBjCGNka0DRDA2k3Kea3J6yZYojRdzYxtkhyknfKtXfD4FPT/+pxpW8P2XD6RLD/6J0mdx4BADiPEUBkSk0D9qhfRiau2me8/8otDeTJH1bJoyFaiShqq6tZ/9fBCElWt5XzUuvVOtS6FCTV1PDr0zbG/LvliaMChduaVTTCfqgR2L+1rWrsYlDKxChiZn0bi+YLPrqninDUWqcKRa31RgQA2I8AiLBeurmBLwCqtWMLhna2vWed18KhnSXt3HljzaIqDihVKHxYUS1f4tF19a6QFyaty3KrFX/hpt9LZ3FrsrtbV5E7WoauBGetEwC4EwEQrnaF384QqjjAjNxx2ktO7YvpNeTaWvLq1A1h1yVeWzd2YXfe4x2lPKN7ABCX4vOZEgiiRdVixlZsQ7vXtu36eKc77d7/s1/zSlLGxOhcsH175z/RyWj7Em2fDmhq/L+fu74u4Q8A4hwjgIhrai/dh7/83SgaeaZXnaCByKzaZULv/PHlvc3lpV/Wywt9LhZWdKxVShqUL2y8RNMfT3c1prwLh1hXF8z/Hmotv+9MkWd/vlg1bVexifo/b3mxuy2fGwDgLAIg4lrPBmWNl6x6qGM1ebBD6F1QWlUrIT/5takJbFsTLReD38XwN+qmBtL/k8Xy+HWZj2heVaGINCxfWOZtOuRoix8AQPwiACKsPH9uIaYUtdD+JJ4Muda+aeNItatZUjaMuM5UJbUa+Rxz16W2OQAAZIYAiLBy5sguc//ZUc5d8Di2dRkuirSNDgAAmSEAwnKFaqKZ/GhbWz+/3QUjAABYRQCE1qqVzG80rrbDvW2ryObkE9K8SvT21AUAIBoIgNDWE91qG42M7fJkjzq2fW4AALKCPoDQltpeTlXzAgCgG0YAoZ3pg9vJnqNnpG7Z6PbwAwAgXhAAoZ3qpQoaLwAA6Ir5LwAAAM0QAAEAADRDAAQAANAMARAAAEAzBEAAAADNEAABAAA0QwAEAADQDAEQAABAMwRAAAAAzRAAAQAANEMABAAA0AwBEAAAQDMEQAAAAM3kjPUJxDOPx2O8Tk1NjfWpAAAAk1L/fN72Po/riACYBcePHzdeV6hQIVrfDwAA4ODzeOHChbW83tk8OsffLLpw4YLs3btXChYsKNmyZYv6XycqWO7atUsKFSoU1c+tM64r1zXe8JjlusabeHjMejweI/yVLVtWsmfXczUcI4BZoB405cuXFzupHx63/gDFM64r1zXe8JjlusYbtz9mC2s68uelZ+wFAADQGAEQAABAMwRAl0pKSpJnnnnGeA2uq9vxeOXaxhses1xb3VEEAgAAoBlGAAEAADRDAAQAANAMARAAAEAz9AGMkT/++EMmTZpkvN29e3dp2LChLffRjWrs+c0338jOnTulZs2acsstt0ju3Lkzvc+yZcvk119/lfPnz0uzZs2kXbt2jp1vPFmwYIHMmDFD8uTJIzfccIPUqFHD9H3feust2bdvnwwbNszVfcFi4dChQ8ZjNjk5WRo0aGBcWzONaRcuXCizZs2SAgUKGI/zK664wpHzjSdTpkwxrpPq93bzzTeH7dt65swZ+fHHH2Xz5s2SP39+ad++vVx99dWOnW+8OHr0qIwbN85o9Pz4449L0aJFw95n9+7d8t1338mxY8ekRYsWcu211zpyrgiNEcAYGDNmjDRv3ly2bt0q27ZtM0LHJ598EvX76ObgwYPSuHFj+fDDD41f5M8//7zxC1y9Hcr1118vAwcONH6RqV1d+vTpI7fffrvW+0MGM2LECOMX9oEDB2TFihVSv359mThxoqn7qsfpc889Jy+99BL7ZgfYsmWL1KtXz3hiPHnypDz66KNGAMzs8ad2IBowYID06NHD+H7s2bNHunbtaoQWXKJ+rvv162cEjpkzZ0qdOnVkyZIlmQZx9bgeNWqUpKeny8aNG6Vt27by7LPPcln9jBw50riW6udf/Uyr6xuOuu7qPuoPFnW8+r6o7w9iTG0FB+ccPXrUU6BAAc/rr7/uu+3ll1/2FCxY0HPs2LGo3UdHjzzyiKdWrVqeU6dOGe8nJyd7ihQpkuG6BVqwYEGG95ctW6aeeT3Tp0+3/XzjxZYtWzw5cuTwfPPNN77bBg0a5KlQoYLn/Pnzmd533bp1nrJly3rGjBljXNddu3Y5cMbxo0+fPp62bdv6ruOmTZsuu9aB3nrrLU/+/Pk9Gzdu9N125MgRz/79+x0553gwf/584/Hm//N94403epo1axbyPh988IEnT548nuPHj2f4PauudbjHuU7mzZvnOXnypGfu3LnGNd62bVvY+zRt2tRz6623+t7/9ddfL/v+wHmMADps2rRpxl/6/fv399121113GVOX6mPRuo+OfvjhB7n11lslb968xvslS5Y0RknU7aG0bNkyw/tqWj1HjhzGiCAuGj9+vDEdpkZHvdQIlLpGmY2opKWlyV/+8hd55ZVXpGLFilzOAGqUSY2iqNEQ75Rv9erVpU2bNpk+Zt955x3jPv5T8GoKrnTp0lxjv98F6vr4/3yrx+zixYuNEdNgihUrZoy8njt3LsP3SF1bXfeKDaZ169aSL18+08er661+T9x5552+29TIarVq1TJ9nMN+rAF0mJpWKFKkiBQvXtx3mwoqao2K+li07qObU6dOGb9o1BOoP/W+lZD81VdfGVNsgcFQZ+oxVqlSJcmZ89KvC+91Vh9TSxOCGTx4sDGl1rdvX5k+fbpj5xsvtm/fbgSMYI/Z33//Peh9UlJSjKne4cOHG8F8+fLlUq5cOSOclyhRwqEzdz/1uAx2Xb0fU9cs0E033WRc986dOxshR00Jq2v97bffOnbeiWjDhg3G62DfD56/Yos/axymRvKCLYJXYU59LFr30Y33OgReJyvXaM2aNfLQQw/JkCFDpFatWracZzwK9vhThQdqpDTUtVUL6VXB0nvvvefQWerxmPWut3rttdeMkUD1x4oqIFGjXaFCo46CPWbVdfV+LNQfkSqsnD17VgoWLGg8xlXhklpzjax9L7L6uxn2YATQYeqXSrBFs6qqSn0sWvfRjfc6BF4ns9dI/SV6zTXXGCMpapEzMl7bwOuqlh+oqulQ1/aJJ54wqrC913LHjh3Ga/V+t27dpGfPntpf4kges97b1bSk/8i2GrVSFda//PKL9tc11GNWXVf/axjo5ZdflkWLFhm/C7xTnJ999pkxddypUyem2KPwOPevVOf5K/YYAXTYlVdeafwgqJYPXvv37zeqI9XHonUf3ah1f2qa0jvd4KXeD3eNNm3aJB07dpQuXboYFaus98lIXT/vdKX/dfV+LBhVzaquqVq6oF7UGkLvX/3eNZq6q1y5stFSx8pjVk3zquUfgdPu6n2qgC9R1y/YdVVq164dcgZAteHxX9+m2pWotazqdwQi473egd8PFbR5/oqxGBSeaE1VmKnK1BEjRvhuGz58uHGbf/XZ008/7Zk5c6al++ju//7v/zxVqlTxpKamGu/v3r3bqJR+7733fMdMmTLF89xzz/neV1WXqkr1jjvuoNIvhJ07d3py5crlGTt2rO+2v/3tb56qVat6Lly4YLx/+vRpz+OPP+5ZunRp0M8xbdo0qoCDuO2224zK1PT0dOP9lStXerJnz+756aeffMeoimBV+etf7d66dWvftVev1fuqohgXLVmyJEM1v7pG3bt3NyquvVTVtHrMqt8BypNPPukpXbp0hs4K//73v42qbCqsL5dZFfDo0aMz/L5o06aNp1evXr7H7OTJk437hvp9AWcQAGPg66+/NtoN3HLLLcZLUlKScZs/1Xrg+eeft3Qf3aWkpHgaNGjgqVOnjue+++7zVKxY0dOlSxdPWlqa7xj1S7548eK+9ytVqmSExCFDhhhPBt4Xb/jGRW+++aYnb968nn79+hlPpOrx6X+N1LVXv9A//fTToJeMABg6XKvHYJMmTYxQXapUKc/tt9+e4Zg777zT07BhQ9/7hw8fNh7nqrWGCoPqtWrJs3nzZh6uftTPdOHChT0DBgzwtGvXzvi5VwHba9WqVcZjVj02vY/hevXqeSpXrux56KGHjN+zuXPn9owaNYrrGvCzrH5Hqsepun4DBw687I+/zp07e3r06OF7X113df3bt29vfD8KFSpkfH8QW9nUP7EehdS1AazqUq+oBruqJN7f66+/bkw/tGrVyvR9cLH1yM8//2zsBKIKOdR6M/8p3blz5xqL5QcNGuRb96MW0gdSu4H4X3uIrFq1SmbPni1JSUnGGr6yZctmuO5vvPGG0XZHVf4GUmsAVYW1KrJRC+yRcT2lquhVSzxUGyK13syfahWjPqbWonmp6fgJEyYYj3PVYkc9zplav5xa0/fbb78ZSw969+6doZOCqvIdPXq0/PWvfzWWjyhqXataW6mmfNXjVP0OUGtZccm8efOMl0D+P/tff/210TVA7b7idfjwYeNxrpYzqS4LoboHwDkEQAAAAM1QBAIAAKAZAiAAAIBmCIAAAACaIQACAABohgAIAACgGQIgAACAZtgLGAAAuIpqUax6Mi5fvlx69eoldevWjfhzqT6k6nNly5bN6FfovyexzhgBBAAArjFjxgyjkb9q1D906FCjeX+kXnvtNalTp45Mnz5dVqxYIddcc40sXrw4qucbrxgBBAAArqF2bpk0aZJUr17dGLUL5fTp08ZOObt375aqVavKddddJ7lz5/Z9XI36DRkyRH755Rdj9ywlNTVVDhw44Mj/w+3YCQQAALiSCoCff/659OvXL8Ptars+NZqngl+9evVk6dKlcvToUZkzZ46ULFnSOKZ79+7GNpVqRBGXYwoYAADElTvvvNPYx3nmzJny9ttvy4IFC4wwOHz4cN8xah/oDh06GFPIb731lrEf+cGDB2N63m7CFDAAAIgbaspXhbvmzZvLq6++ahSMqBc1WqhuV9T7akRw8uTJ8s033xhBcPz48fLAAw/ITz/9JO3btxfdEQABAEDc2LNnj/FaTe8eOnTId7uqFG7Tpo3xtgqD+fLlk/3798uaNWskT548xu233367PPLII7Jy5UrRHQEQAADEjRIlShivb775ZunUqVPI41Qlce3atX3hT1Ejf99++61vxFBnrAEEAABxo1q1alK/fn155ZVX5Pz5877bz549K6tWrfK9rwKi6iN47tw5322LFi2SmjVrah/+FEYAAQCAa6jGzapgw2vChAnGur+GDRtKt27djNu++OILo+1Ls2bNjBYvKSkpMnv2bBk8eLARDpVHH33UWPfXunVr4xg1FTx16lRjDSBoAwMAAFxk8+bNMnr06MtuV2Hvxhtv9L1/4sQJ+fHHH43jy5YtK127dpUqVapkuI8aFfzuu+9kw4YNUqZMGendu7fxGgRAAAAA7bAGEAAAQDMEQAAAAM0QAAEAADRDAAQAANAMARAAAEAzBEAAAADNEAABAAA0QwAEAADQDAEQAABAMwRAAAAAzRAAAQAANEMABAAA0AwBEAAAQPTy/+WYCE2SHFafAAAAAElFTkSuQmCC",
"text/html": [
"\n",
" \n",
"
\n",
" Figure\n",
"
\n",
"
\n",
"
\n",
"
\n",
" Figure\n",
"
\n",
"
\n",
"
\n",
"
\n",
" Figure\n",
"
\n",
"
\n",
"
\n",
"
\n",
" Figure 24\n",
"
\n",
"
\n",
"