Skip to content

Instantly share code, notes, and snippets.

@JiaweiZhuang
Created March 28, 2019 16:20
Show Gist options
  • Save JiaweiZhuang/7a620ccea7d686891fed0221f386a50b to your computer and use it in GitHub Desktop.
Save JiaweiZhuang/7a620ccea7d686891fed0221f386a50b to your computer and use it in GitHub Desktop.
debug_xESMF_CMIP5
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Debugging https://github.com/JiaweiZhuang/xESMF/issues/49"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Preparation"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"%matplotlib inline\n",
"import matplotlib.pyplot as plt\n",
"import cartopy.crs as ccrs\n",
"import numpy as np\n",
"import xarray as xr\n",
"import xesmf as xe"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"# https://drive.google.com/file/d/1BCSJGoU8W1CSzXoqILlL-BFvYuXsaH41/view?usp=sharing\n",
"infile = './thetao_Omon_IPSL-CM5A-LR_historical_r1i1p1_1850-01-01.nc'"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.Dataset>\n",
"Dimensions: (bnds: 2, i: 182, j: 149, lev: 31, time: 1, vertices: 4)\n",
"Coordinates:\n",
" * time (time) float64 15.5\n",
" lon (j, i) float32 ...\n",
" lat (j, i) float32 ...\n",
" * i (i) int32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 ...\n",
" * j (j) int32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 ...\n",
" * lev (lev) float64 5.0 15.0 25.0 35.01 45.01 55.03 65.06 75.13 ...\n",
"Dimensions without coordinates: bnds, vertices\n",
"Data variables:\n",
" time_bnds (time, bnds) float64 ...\n",
" lon_bnds (j, i, vertices) float32 ...\n",
" lat_bnds (j, i, vertices) float32 ...\n",
" lev_bnds (lev, bnds) float64 ...\n",
" thetao (time, lev, j, i) float32 ...\n",
"Attributes:\n",
" CDI: Climate Data Interface version 1.9.2 (http://mpim...\n",
" Conventions: CF-1.4\n",
" history: Thu Mar 28 13:15:50 2019: cdo seldate,1850-01-01,...\n",
" source: IPSL-CM5A-LR (2010) : atmos : LMDZ4 (LMDZ4_v5, 96...\n",
" institution: IPSL (Institut Pierre Simon Laplace, Paris, France)\n",
" institute_id: IPSL\n",
" experiment_id: historical\n",
" model_id: IPSL-CM5A-LR\n",
" forcing: Nat,Ant,GHG,SA,Oz,LU,SS,Ds,BC,MD,OC,AA\n",
" parent_experiment_id: piControl\n",
" parent_experiment_rip: r1i1p1\n",
" branch_time: 1850.0\n",
" contact: ipsl-cmip5 _at_ ipsl.jussieu.fr Data manager : Se...\n",
" comment: This 20th century simulation include natural and ...\n",
" references: Model documentation and further reference availab...\n",
" initialization_method: 1\n",
" physics_version: 1\n",
" tracking_id: a73c6f24-c8fb-46c9-a9dc-a41507a05368\n",
" product: output\n",
" experiment: historical\n",
" frequency: mon\n",
" creation_date: 2011-07-04T20:11:18Z\n",
" project_id: CMIP5\n",
" table_id: Table Omon (31 January 2011) d2d6beec2b8fea5bbed3...\n",
" title: IPSL-CM5A-LR model output prepared for CMIP5 hist...\n",
" parent_experiment: pre-industrial control\n",
" modeling_realm: ocean\n",
" realization: 1\n",
" cmor_version: 2.7.1\n",
" CDO: Climate Data Operators version 1.9.2 (http://mpim..."
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ds = xr.open_dataset(infile, decode_times=False)\n",
"ds"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAVUAAACOCAYAAACBgpt2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsnXmcXEXV97/n3t5mn8xkkslkTyb7DlmAsO9hFRAERREVeERlR0QUEF4EQQVlEYFHfBRQdgxogABJMAECBEIWspNlJtvs+9Ld9573j7rd0z1bZkICeXz6N5/+9HRV3bpVde/93VPnnDolqkoKKaSQQgr7BtaX3YAUUkghhf8kpEg1hRRSSGEfIkWqKaSQQgr7EClSTSGFFFLYh0iRagoppJDCPkSKVFNIIYUU9iF6TaoiMkREGkTE3h8N+k+CiGwRkeO/7HakkEIKXxz2SKrtiUFVt6lqpqo6+7dpew8R+aGIfCgirSLy53Z5w0REvRdD7PPzburaIiLNCWVf3+8d2M8QkVu8MeiS8EXkMBF5X0TqRWSFiByekCcicpOIbBOROhH5u4hkJ+QPFJF/iEiViJSKyH+1q/t0EVnljec7IjI+Ie98EVknIrUiUiYi/9Ou7mEi8i8RqRaRXSLygIj4EvKnisgyEWnyvqcm5B0jIgu8urd00ufbRWSliERF5NZ2eQNEZK6I7PDGblgnxx8vIh+JSKOIlIjIeT3pc7s63vLqT+zTAhEp98b6ExE5s4tjH/eOLe4sP4UvBv+p0/8dwP8D/tRNmVzv5ZCpqrfvob7TE8qeuO+a+cVDREYCXwV2dlMmD5gL3APkAncDL4tIH6/It4BvArOBIiANuD+hiieAzUB/4FTglyJyjFf3KOBJ4L+8ul8G5iaQyBJgtqrmACMAH+ZaxvAQUAYMAKYCRwGXe3UHgH945+8D/A/wDy8doBFzT1zfRdc3Aj8G/tlJngu8CpzT2YEeST4F3ATkeG1b1sM+x+r4htff9rgSGKCq2cClwBMiMqDdsYcDI7voVwpfJFS1yw/wV8zN1Aw0YG64YYACPq/MQsxN/45X5mUgH3MT1QEfAMMS6hwLzAeqgHXAed214fN8vHb9uV1aUvt7UMcW4Pi9PH/8WMwL7CfAJqASeAbI8/JeBX7Y7thPgLP3w5jMA07prl/AacDqdmnrge96/z8HXJ+QdxjQAqQDmd74FiTkPwL81fv/h8A/E/Is7/46rpN2ZAJ/Af6VkLYGOCXh9z3AH73/TwS2A5KQvw04uV29xwNbuhmjJ4Bbu8jzef0b1i79KeD2Lo7ZY58xRLweOKS7+xOY6Y31zHZt+hiY7B1bvK/vm9Sn559uJVVV/aZ3U8Yktbu7KHo+RnIZiHlbvgs8DuR5D8EtACKSgSHUp4B+wAXAQyIyobNKReQhEanp4rOiu7b3AFu9qenjItJ3D2Wf9KZfr4vIlL083xXAVzCSVRFQDTzo5T2FGQsgLvUMpXOJiW7GpEZEftJVA0TkXCCsqv/aQ1vF+7RPm9hFvgBBYFRCevv87o5NzEdEDheRWqAeIxnel1D+d8D5IpIuIgOBOZiXEsAEYIV6TONhhZe+v3EIgKc+2CkiT3gSP/Sgz8AvgT8AuzqrXEReEZEWYClGkPkwIftq4G1V/bzPRAr7APtq+v+4qm5S1VqMJLRJVd9Q1SjwLDDNK3caRkJ4XFWjqvoR8DxmOtoBqnq5quZ28Zm8l22tAGZgSOtgIAsjVXeFb2Ck26HAAuA1Ecndi/NeBtykqqWq2grcCnzVmwK+CEwVkaEJ53zBK9cB3YxJrqre1dkxIpKJeXCv6kFb3wGKROQCEfGLyEWYl2W6lz8P+J6n38wBbvDS01W1HjOF/7mIhETkIAwxxo6dDxwlIkd70/KfAoGEfFR1sZrp/yCMJLoloW2LMCRZB5RiyOUlLy8TqG3Xl1rMNd7fGIQRLM7BvFwSVSLd9llEpmNUKfe3rzQGVT0N049TgNdU1fWOHYy5t27eD31KYS+wr0h1d8L/zZ38zvT+HwrMSpSsMARSuI/asUeoaoOqfuiR+m7M1OzERGNIu/JLVLVZVZtU9U6gBjhiL049FHgxod9rAAfo7xHRPzESP953d0S/N/gFZgq+eU8FVbUSOBO4BnMtTwbewJAYGL3k3zAS02rMy4aE/G8Aw4ESjPT1ZCxPVdcCFwEPYPS6fYFPE45NbMd2jBT6dwARsYDXgBeADO/YPsCvvEMagPbXMRsj8e5vNGOEi/Wq2oB5gZ0C3ffZ69NDwJWeENIlVDWiqvOAk0TkDC/5PuA2T6BJ4QBAT0h1X4axKgEWtZOsMlX1+50VFpGHJdlKn/hZvY/aFOtf++lud+V7WjYRJcCcdn0PecQBhqQuEJFDMVLOgq4q6mZMGkTkp10cdhxwhRiL+S5gMPCMiNzQWWFVXaSqM1Q1DyOBjQHe9/JcVb1FVYep6iAMsW73PqjqVlU9TVULVHUWRsf+fkLdz6nqRFXNx6iGhmJ0753BR5sBJs9r9wOq2uqR/+N45OW1Y7KIJF6fyV76/sYKunlWuulzNjAdeNq7LrFxKBWRrl7eiWNyHHBPwnUFeFdEvv75upPCXmNPSlfgPeDShN/D6Gio+l5CfpJxCGMU2Oj9nwVsxTykfu8zAxi3LxXFmJsuBNyJMbaFEto7C0MQFuZhfxpY0EU9QzDTsoBXx/VAOZDv5R9thrDLdmyhzVB1tTdWQ73fBcCZCWWDGD3rfODefTkeXv35mBlB7FMCnAtkdlF+mnd9sjHS0JKEvDzMQy3AeGBVu3tknHetA8CFGJVLouHqYMD2xuBp4KmEvG944y4Y4lmEUYXE8j/DGPx8GEv6i8CTXl7Au7+u9Mbzh97vgJdveddxjpceiuV5+X4v7SnvPg4BdkJ+CCMhq3cPhRLyvoPxeBiBmdY/g2ec667PXj8Tr8sMr/6BXn/Geu1N89p3IRAGDvKO79fueMXod9P29T2U+vTwWdtjATMN3IaZ9l7H5yBV7/cYzFS3HGMFfwuYuk87ZfSV2u5zq5d3gXfzN2KmYn8BChOOfRh42Pt/AkYCafTa+iYwPaHsN4F3umnHFpKt/9dgPB7qMV4Av2xX/r+9ts7Y7xe+nfU/sd/e779h9JG1Hgn0S8gb7fWjCUNO17Sr+yrv+jYCixPHzMtf7I1BFfBHICMh7w6MKqDR+34E7yXm5U/17rlqDFk/265t0zCuTM3AR8C0hLyjO7kvFibk/7mT/G8n5LfP03b9+oXX73LMy7xPT/rcro5hJD9f4zDGqXrMM/gBcFY31zVl/f+SP+JdiBT2AiLyGPCsqr72ZbclhRRSODCQItUUUkghhX2I/9QVVSmkkEIKXwpSpJpCCimksA+RItUUUkghhX2IFKmmkEIKKexDdBYRp0ukSYa20LS/2pJCCil8fmxV1WFfZgNOOiZDK6s6RgZdtqL1NVU9+Uto0heKXpFqC00cb52354IppJDCl4I33GeG7rnU/kV5VZQlrxZ1SE8v2rKnwEX/EegVqabwxUFVUVxcXCThD4TkVZgppHBgQVEiB24M+/2OFKm2g6pLlCgRwjRQQ8QLFFVDJRFacXFRFAsLhyiKi2BhY8e/bXyAMECGkUMeUaK4OKAuNj5EbBqoxU+QEGls0TU0Uk8rzQA4ODRQiyBYWBh6TbhJPddiC5sAIYKEmCKHE5DgFztYBxjE6viyUVe7TE9h/0CBVlKk+n8WqkoJG6jWMtLIZBvrAcjOzmby5Mnk5Q2ktraWq75+Kf369SM9PR3btmlsbCQrK4tgMEhraytr167l8ssvT6p7h27G5/ORlZVFWloalmVRW1uL4zjYTUEKRvThs88+i5dfsmQJruvi8/mYNm0awWAbSQ4YMIBdu5JDbbo4tNBIC41UTdiEz+ejpKSEyspKAAoZykRr1j4dK4coFjZRwvgJ9lpqbtEmWmjyXkoONj6y4xsK9A6dkWVXZVx1qaOKBq0lQpgc8siVAizpua02RcQ9gwKR/8OLig48UjVhIiF2syf+bp/Xm2pVqaOagCcdilfHSn2XMi/q3FlnnUXDogqqq6upq6tj8eLF8eN37tzJunXrAHjppZfYtm0bWVlZbN68mQULFrBmzZp42TFjxnD88cdz4YUXMmvWrCTiCYfDiAh+vz/eri1btlBWVsasWV0T4M6dbbufqCrV1dWUl5dz2WWXcfTRR/OLX/winh8gSIA0ciS/1+MEhvg266dEiRAgRAO1VFPWeWGFw+U0QpLeeX47LNZXOqQdwRkEJdSj43tCpO3hqsObznMAZJNHjuSxXktopZkRMoFMK5s0MgmS1u1Loqfnbi8dd0bG7ev6TyJsVSX8f5hUe7VMVUR0vxqqYqQJbSTaAwJVVSKEcXHwE6BJG1iK2Z8vj/4cZB1FiW5knX4UP8aHnygRAPwEidBKGpkUUEQzjZSbKHZkkssEmUEGOVhi0aA1vKem7sMPPzyJeGP45JNPmDx5b2No9x6qimVZjBkzhiVLlnB+v8v3fFA3aNQ63lUTTD+NDIplMnVazVbWxssUFhYybtw4Vi1Yz3iZgT++DVTP6m+lhSxye3VcZ+gp0bVoEyXuRnZrCS004iNAhI4xwI+1voptHXgbBfeUdN9wn0FVv1Sl+8TJAX3+nx1tUmOH7FymqtO/yLYk7L7QKVS1al+f88CSVNsT6B4ItUnr+Uw/ZRdb42kFDCIxrGWYFqq0jHwtZBAjCZDGiOMLWbVqFf3792fFihVxvamfAHXUUJMglTVQQz01tNLCcvffAOTSl+EynrolHfVGF110EeFwmJnWCdjYpJOJJTbznad7Oxo9Ql1dHccddxwAmZmZ5OfvnXSaiAzJ5kjOYLOuoYQNbNJVRIkwiGLKKCFMK7t27SJ9d1+mWLP3qv6MDrGk2+CqQyN1tNBMlAjNNFCpu4nQykH2UaRJRq/PGZJ0hlljGcAwUKWFZkrdjVSx29OTmxf6W+5zDKQYCwsffrIlj3wpxJZ2j0oPX/hhbSXs6cr9BAkQ6lwaThQoOsH/JslWgfCB4wK/jK5jICsmVOM+xYElqfa0HZagqmx2PmUTqwDoQwHVlHd73IcffsgjjzzCI488AoDf72fgwIFs2bKlQ9lQyNz8zc3N8TQLm3wKGSfT40ahqEbYoJ9QTTmKZ4jybiiHKE00kEs+42UG6ZK1z8m1vr6e7GxDUNM5mtzOtttKePjjD2c7QlCnc8NCizZRTTkBguTL/tugwVWHtfoRDdRSR+fCgyAcap9MRuebNLSrz2WHbqZBa6jRSlppJkqYEBn0H96XQCDA9OnTyc3N5YYbbmDgwIEsWLCAjz/+GMdxaGho4Hf/70Hqqaa/DGay7/A9nrP9lL/GLecD540O5fJlAJW6EwubIGlM9R1JpuR0XXEnhNsVqR4Ikur4yQF94pWO98rBQ0u+cEn1y8CBJakmQCyJ3zjt39KuurzpPAuYabyDQw3GOJNGBs00xssOGTKEkSNHsnXrVv72t79xxx130NjYyJIlSxARMjPNTi9+v59IJIKIIeyWlhYALMsiGAzS3NyMhY2FxW62MVBHYImNT/yMk67vk7C28pmu5gN9iwnM3HcDBOzevZsHH3ww/jtEO72mR5ydEan4vUvvSU2WFYz/1lYjuaurpJFFmrfFU1fE+3lRq1V8oIZ8QqRjYVMog8mTQjIlm3TJwcZCEZpppEWbCRJEUQQrfs02uCtooh5FqdAdAPzmN7/hkEMOYfjw4eTm5pKWltZlO4477jiOO+441qxZw/jx4+PprbSwKvoutvjJkCz6WP1o0SbK3FICEsTCJs8qpA/941Ko2JBuZTOEMdS6FdRpFYqR0gvtoThOlBotp5kGRDoSchIS7QnxpE48Gg4QPaaL0KIHHrV4W6yPwgQbB0BV397X5+l1z1VcLLF7Nf1w1EEASzrqqmI3R2euL13py7brJgBm2SeQQTZ9vgLPP/88QJxQs7KysCyLbdu2sW3bNgDuv/9+MjMzefJJs/1T7GEMBAJxUp01axY1NTUcf/zxPPDAAxx77LHMmTOHe+65h127drGbEnZrCWVsZyijCZFBBtldGjgCEmSsHESBFvGxvs1Eaxar3KU9HruuUFpayuDBg+O/Z3F8srGok6mp2Gb8rVxP0vP7wXXB512XGCF45bWqpo1Io1HE8qGOs0+nnlt0DRt1pSfdK/1kEMXWpLbpttcPRx02usvZ6W4mQjipjv4yBIcIFbqTgTKCB1/8LZZl0b9/f2bO7P2LbNy4cWzdupUNGzbw7rvvEgqFyM/Pp6amhldeeYXq6t2U7d7NyJEjOeWUU2hoaOCOO+5gov8wcqWAJq0nICEcogy0i9nmrIvX3awNbHE+xU+QAmsgOZJPUJKJPlGg6A0cut3i6guDIkT2glRFJAS8jdm1wQc8p6q3iMhwzD5leZjA499U1bCIBDFB5g/GBJH/mqpu6aLu72F2hBgELMfsjvAucGyvG7oH9Lrnq9ylTLYPSyLD9ghrC7vYRr3W0EQ9jdTh4pCmmbi45NEPvwRJkwwCbhAfARSXTM1p0zklGKocjVKnlUSJ0qh1KIqfIO87b6K4ZDyfPB3Mysrioosu4vrrr2fnzp089dRT3HrrrRQXF7Nu3TouvfRSNmzYwAcffEB6ejqhUIjc3Fy2b9/OCSecwO23305ubi6PPfYYgwYNYvjw4RxxxBG8+OKLzL37LZ5c9kfmzZvHz2+8mX5FBdTvaCZb+1AoQ+kjBZ2OW0yKXK3vIyK88MILnHXWWb0dfgBaWlqSCDWfQrIsTx/fhUQjPuNtYOV5G8F6BNueUOOIRJGsTKjz9swLBIiEm6iVClqsJurdKsrd7cbR2yO5EOlMl2Pw99BfViyhKdrgHZtGM41s0/WMYjIRDWNjx12e3ovOIyBBDvWdwpzbZ3HjjTfG67n4hq8ze/ZsHn74YS6++GLOPPPMHp2/OwwZMoQhQ4bE9dUxXH311Z2Wr66u5qGHHiIvL4/Ro0fT0FBGRkYGZRUVsAkWLFjAIYccQnV1NZWVlVRVVVFVVcWzzz7L0089Q7pkMWLyEBobG9mxcTdpksl4e2bbizJBWm3SeiKECRAEBBubeq3lI3fh5+73voCqENa9Mva1AseqaoOI+IHFIjIPs2PGvar6dxF5GPguZkPJ7wLVqlosIudjNoD8Whd1X4nZquY9VT1GRMZidmrY5+i1TvUIOY2QnWwoUFc9f8+NhLWZHWwhj37kSl8yyGHe7qdJT09n6dKlrFu3jsrKSlzXZfPmzVRVVVFZWcnqJetppA6HKJmSQ64UUKPlNGodABmSg6MR/BKkQWvilnsffvwEPSNTsi7uiSeeoLS0lLfeeovKykr69evHvHnzAJg6dSrLly83dfh83HnnnVx//fUAjBo1ivHjx7N8+XJKSkpwXRcQApipXguNDJbRLN2+kPz8fN59913uuusuXn31VQ6WY7ok1vXu8rgfLPR+uhaJRCgtLeXNN9/kkksuiadPlEMolCGxi5F8UEwF4JGoBPxYuTlgeZKsJW2E6vPFTgSuojVmg84mt56t4dXsiHxGppXDqRecyIABA7j44otJS0sjJyeHc845h4ULF5Kenk6gKYOgpJElOfTVAcZzwrao1UosbNLcDHy2IfkWbabC3UGF7iRTcokSRrDYpka6CxAkR/pS7u2PuGbNGrKzsxk4cCAAn332GcOHD+/VOO4vVFRUUF5ezrhx43p1XHNzM5988gnBYJDMzEyi0SiXX345CxcuJBgMEolE6NvX6MqLi4t55513APMSUxQHB1AyyKaWyi9dpzpqUrr+9h/FHdLPGLmyxzpVEUnHbEHzfcz2S4WqGvU2xrxVVU8Skde8/9/1tnrfhdkLrcODJSIfqOoMEVkOzFLVVhFZrqpT976nnaPXkupH6QvIbMqnj/QlU3LIJBcffsK0sl4/BmAQIxktU3nTfS7p2GOPPZZjj+1e2m5oaOCDDz5g6dKlTJs2jZkzZ+L3++O6TzCSWn19PSJCfn4+a9asYfXq1Sxbtoxf/crsVpyXl8cPfvADamsNMRxxxBEsXLgQn8/HmDFjkggtGo1SX19PNBrlkksuYdWqVbz88svk5eXxySef8N2pV/O+8wazfCdhYbMo+gIlup6ioiJGW9NY53zE4sWLefXVV1mmC8jT/gyTsWSSGzdoiSWMkslsc9YTJJ1/vPpCj8b7ePkq29jAFtYSJUKQNOPmxCT6MoDM9kabrvRv6iKBAKiaj+NAwG+I1PHKW5ZJ9/lwd5XhRMOsjX5IS241l/zoEq644goGDBjQaTsXLFhANBqlrq6ON954g0gkwtKlS3n55ZfZvaWCKTqb9xOMNgOcodRoBVEicUk3Nm32EzBSr+94/rb5QRYtWsTjjz/Om2++ycKFCzn33HPj9cyfP59LL720R2O5v9G3b984+fUGaWlpHHLIIUlpb731Ftu3bycvL49AIEB5uTHCLlq0iOLiYm655RZGjDCGa1Vl9+7dLF++nDlz5nz+jnxOGOf/TiXVviLyYcLvR1T1kcQCImJjLPbFwIOYvdxqtG377lLMpoh43yUAHuHWYja4rOjk3KUikgu8BMwXkWpgx150b4/otaT65z//mcsuu4zW1lZCoRAtLS2MZgqttFBDBbWewQjAcRws64t1rSgpKWHjxo2AWRW1fv16Ro4cycyZM2ltbaW+vp6PP/6Y9PR0cnJy8Pl87Nq1i6OOOiquF21tbeWOO+7g9ttvJ4DRjeVIPlOtw4kQZZ27DAsfu3QLY63prHE+YNOmTcyePZvi4mKmTJnCQw89RAEDk1yOXHF4y3meCdZMiuyRvB5+qkP7VZUZ9nHUUUUT9VTqLoKEGM8M457lD8SlzjWtSynRDfShHwGCRInQRD3NNDLdPpY8vwlqIQEjFUpeLlpZjeT3aZNKPYOV2oI0NINYaEUlzW4Dy5sXMOfc43n00UfJysqKt29q7onURHZxw6+/z2WXXYZttz1AN954I3fddRcAs2fPxu/3s3Dhwg79PO2006icFyWDbOZHn8Z13R7fK7HrdOqpp/LUU0/FvR9SiNsJvlRJdeSkDP3li+M7pJ8/6sPeSKqxnXJvBh5X1WIvfTDwL1Wd5G1Tf5Kqlnp5m4CZarYu767uo4Ac4FVVDXdXdm/Qa1LNIZ8GasmjH+XsIINsRjGZNSyjlWZy6UuN96I477zzePrp/eOfub+wefNmTh/9dbY4a2hUs/5eUQpkICNkAkvd15PKz7CP5/3ofNauXcv06dNpbDSGsjz6M15mJBmPxBJ26FbWOcs4OnA281v/3uH84+yD2eyuoYCBZFrZ5EoBmdIHO80YLGVAf7AN+by65s74ccNHnEjAn0FrtIEtm94AlJNGXGsesmCCg32Mt2wbCUdRv20kV0DqGtGGRna0bmBd41LuuPs2rr322jiJua7LyKyD2NL0Sby62Xnnsbiy7Ro7jsN5553HCy8kS+KffPIJzc3NTJw4kYyM3vuZJmLTJmOoHDly5Oeq5z8FO3fupKioiAyyaKT+SyfV4ZMy9dYXOi5++fbod3vlUiUit2B27L2BvZz+i0i2qtZ1tQgg0flfRH7bg2bVqeqt3RXo9fR/hnUczdpAOTsJahqlbGI5i5kohxCy0ljmLIyX3Zup0JeFcDjM73//e66//noKrSH4MUSk3kKCct1OrVYy0p5MvVtNQIK0agsfOG9QYA2Mu/AApJPJQdZRSfWLJYhtM5ARlLODj6KLOCbrfBbUtxHrb3/7Wza6K5liHU6e1a9NDxoMItlZEI2iQXPJon3SOTLv57y95HYASksW47oO2TmDGTbyeIoGzqQ1lIWvMYo4LhJtUwmo7RkZ0wPx9NaSz9ja8Am7mzaSZmVyUPZJXHfddfFjXnvtNc459QIanWoAhqVPYVj6lA76ddu2ueeee+LxDi688EKysrL26Qqz/0tkWllZyb/+9S9OOeWUThd2uK7L9u1G39xI/RfdvE6hKl1N/7uFiBQAEVWtEZE0zPb2vwIWAF/FeABcBPzDO2Su9/tdL/+tTvSpTwGnkbwIIPE70fn/HOC2PTTzOuDWbvvRW0n1BN/X4no7VZcWmqh1q8inkDV8SLWW8ea/32D69OmEQh3Xc+/atYvrr7+eO++8kx07dpCVldVrxf6+RktLCwPThyFYDLZHUSBGZdOiTdRoBSudJRTbk9norCCDbKJETDARhBaayCQHh2jcUDbVdwQFtIsnmWAwUnXZyCoqo9spTp9O0Eon29eXhZVPMNI3iYEMjxMqtt1GrmkhyDfBR8L9Mo2RCUOSGjM2CWgsXUDaXV5fU5REOcaKmgLvLL2HfunFvPLO/zB+/PgkFzHXdSkqKmJAyxSiGqGkaTWH5p0NwNzS3xONRrv1/0yhZ2htbWV8+sFUuWU4RBlsFbPR9wnhcJgPPviA6dOn09TURFpaGm+//TbnHXsh1W4ZfoJkSDYuDpW660uXVIdMzNYfP9dRIP3RuAXdSqoiMhn4H8DGzKmeUdXbRGQEbS5VHwMXeoamEPBXYBpQBZyvqp91XvueISLXqeqvP2+ZXkuqEgwiloW6LgKkO0HSnGw2O5+y2y0B4OyjzifbymeIPZa3Wp5GVRERnnzySS688EIAnnvyRVq0kaCVgVg2af5s/L50stMKKcqexKIND/S2aXuNxx9/HMFimn1kPNAKmKWNhTKEQmsIYgkD7OFsja5DcdjmGCv+hx9+yNVXX82qJWs5wXcB5bKTNZH3SQvmkOlkJRmNxO9DAgEEGGXNxN+ykm0tn9Lo1GCLD0tsiuyRJm6qneDqFPs/EEADyU777ZH0OAm4vjayBQgH/IgL/roIbsCCqINv0w5sfAzwD2fChAmsWbMGv99PcbGx4FqWxaRJk1j81jtMzTmB2mgZr5c9iotDIPAwAPfeey9XXXXV3l+E/wPYsWMHpaWl/NfsG7DFRlFq3UoyJYeXNv6FYcOGUSADGWlNoEy386n7ATGX3JNnnUlEW6jXGhQlXTIZao9ljHUwIUmPu8693tpRT/9Fwxiqeu+nqqorMATZPv0z6LhyRlVbgHPbp3cGEXlTVY/bQ9rfujl+jqrO2xOhwl7ko9UpAAAgAElEQVSQqpXXB6zkQMlueSWDdAy+Fh9Zdh41WsGG6Mdsiq5A5BnAWC0vvPBCBItcfz8KtIj87OFk+wpoDIYJSwut2kR1Yynvbv4T/XPeIOTPJpRdwNpNr+zXwMzRaJQq3U2pu4nB9qguy4VIZ6e7mQF2W3D16dOn8+KLLzKwYBDvRV+jRRsYYA3no5Y3Odh/LBlkgljmphdBIxHE70dcZXhgIsMDE3GjEWojZaRbWVixaZOrYFuGhD19KmkhotltPqCuz/K+BeyE8fF8h9UWT4o1yTs2/JvNq+aSnp7OoQMvJRj1zqVKrhSwrmYJRf6R7Iyal/3u3bvp168fH374IUuXLmVQxnhy0wcxU77C+1UvAXDppZfiui5nnHHGXo//fzp2797NQUWHUuZuJ10y8RPExTG+25LBRucThg0bBsBQazR1Ws1uLaFIhrNDN9PXGkg/axBBSSfP6k+tVpIrfeOLaazszDbD436xZ/cOxvn/wAhK40mz6RjPgz60rW3JhvbTSd4QkZNUdVu7Or6FmfL3KE5Ar6f/Jw+72khOXVlqXRdcpSFcyc6qFWxqNpGhAnYGYccYcY4e+F1CoT5oRgg3zY9EHNRvo7YFIoQjTVTUricabmL7rg/J7TuSLZ+9jc+3/5a+rVq1ikmTJsWDaAywhlMoQ5JWgYkl1LvVvB+Zz6GBU/l360vxvEgkwty5c5k9ezaFhYU88sgjXHbZZRzn+5qxjkubX6gEAnGLPJYVNxRpOALRaJsUOmgAboYhUSfdZ0g5Fh/UL0bShCTijE331TudawuiEA038f7cmwmEcgi31DK6z+GMyJ0FYSMKRZrr2dmwlk/D78X7dP7553P33Xdz6LiTCVhpjMk+jF2tn7Gycj6Ohrnyyiu577779s0F+A+E4zhMDBzCJnclRdZwRvomYbeT4MQSXHUQLCJuK+9FXyNdMhlhTWTerr9TUFDAcHsCo3xTAIyPceIsJq72Md+vlvzuS5/+D5yQq5c9fWSH9FsmvfxlRKm6ErgKQ6DbaSPVOuBRVX0goezpwG+Ak2NqBBG5Hvg2MKc92XaF3rOUz2eIwO6cVHfVbWDFjpdxtW3J3Ph+JzIgazzVzSXkZA4mEMxAAfUZEnUyPKOQbaGWYGVk0y/XjH3foQezYtljDB48izGjzmTh2zf3usk9wcSJE6moqCAajfLee+/x7bMuYa1+SB/ph42PIdZoEIvtzkaz2KBdyDq/388555wT/33OOedwzz338N6mefRzB1NkDSddshBscF20pRXJzID0NPD7cDJDqN9qMyLZFuoTXLvtoVEhLpG6PolP7WO6U3EVFfM7rle1oam+nLLVJkRhuMX47Q6aFsDZ6mKHDGn7AwEGyyQ2Vi0nrCbuwd///neef+YlQnYmRRnj2M121tcuYUjuNF5Z/Ke4dJVCR6xcuZLJkyeTJXkc7DuWLMnFdR3PPJJgNHQtLNuPBPyUh7eRTV92R7YAMC7jMAByff3BtrEyMyDd8yZxveXDcVI9YKJCHVCSqqr+DvidiPxIVe/fQ9mXRaQVeE1EzgS+AxwOHLknN61E9J5UPUJVO/kNiZpVVSt3vIKrUYr7H0VhznjS/NlIMIj6bfLpg9oWDmCFHVy/jfoE9Vnxqaza4pGJqbZi+6c01O/EjUZYueYpjjrJz6LXbuzYrn2AmHX1zDPPpNo9k/LychYvXszPfvYzqtaXURYpoV6rOSbtPN5q6t5VLD8/n/Xr1zN//nzmzp3LIw8+xiT7MPpaA9DWViQYNGPnT7gEIkZaB1yPYNUjTtdus9q7PkFtkqVTh2RjlTd+1TvXsfGdJ/jR5Zdx331LALNiDMA3ei3RDeOxsChr2URlZAthbcHn8xGNmpdixG1B1aWqpYSt9cs55/wzuPHGG5MCjqTQERs2bAAgmG/xXsU8+koRFbqD8fZMBsrwjsFtVAkGs3Ach5MG/ZCSxlWsbXqX4uBB9AsORUJByM5KPkli+MG9CN69v9CN8/+XBlW9X0QOA4aRwHuq+pd25V4XkUswMQiWAseoajO9QO9JNUao3hREbTMFkYgJmnL85J+YNE89oCK4PitODjFpNJpppr+uRxziWaHVJ6hlJC21IL1oOKwGK5RGc3MVH733AEedYrPoXz/uddN7i4KCAs466yz69+/PsYefwFDfODS6qedr20U48cQTOeKII1i9ejUfL1pJX7wVSTEVSsxH1HEg2vaMiE9wrbaXi+u3cD2NgWubsYlb9l3QRDN/gv65takaJ9zMX56eTyizgJaGciKhsfzzn89QVFQEvMedd97Jr3/9It/7/rfIy/su3/zmN2lpaSEvL49ly5YxderUpK1dUtgzzj77bBzH4d1332XUqFFMnz4dSqDGLWOgPTxOiGJJXB/aJziYT+repPC4epY8tYLp2afSN81bfpyZ0bbiLYZEQrUseq7I279wVWh1D6woVSLyV2AkJphKbBAVE5AlVqaaNnerdIyUul2MQUdVtduA1zH0PkqV3/YkqmR9jhv0g23IVX0WmrC23EhWHsn6jF4wRp7gkUTQiv8W15vC2hDKL2TaGT8lEm6CxlZW/fthPlz8W558ciDf+MY3etv8vcJhhx3GHx9/iFtuuYXf3nR30rr7RNTV1fHUU0/xqx8+QMAKEiKdareM7c5GXC8IsjpO28NkCTS3QCCA1RxBo6759ttmDCzBiipuwDKSfdCL2OQHJyREg96LygJxBcuJKVTB8m6bQFYfcorGkp0/nGi4mZaGcrasfJnCYVOpra7G7/dz99138+yzz3LQQQfF+xJbpdTdFi8pdA/Lspg9ezaVlZXxfcNqtQpVN+5loq5CayuSFsJvBcnx9+PZJ19iVsE5hEiD2MKN2As49u3N7PDsEGofSJKqEHUPLEkVmA6M7ywuQAL2iWN970nVsoxeL/FCivGTVJ+A30iiWCZfIm6bfjAmrcaINZFUfQl6QMub2or5307Lx+fmE66tApSm+t1ceOGFfPeyK9m6aTX9+/ffF2ORhDkDfgDAvJ0mXulFF13ERRdd1GX5X/7yl9x0001tCQ6kkUmeNQAbP663dcf86NMc6j+FrEYbbTFpVloahIJmigeIzzaSf9iT5kM2TsgHtrkfXBvEaSNOx1OXxMbPchTX061mDR1L1tCx+JuVhtodBJvzcKKt1O5eR58+fbjpppu49tpr99GopdAZVq1aRXFxMStWrKCROnbqVorEBIGJGywBRJje90xEbOPtogpRB2IxG6BN/WZZbavjDrAty830/8DR8XpYBRQCO7sqoLpv9tXufc8TCFUtc3FjpGmI0kilrs/CtQUnZOMELZyghRswqgPXNgTrev+rDy/NSGFGfwhuLN0WHL9g983noO/8mjGHXwxAa2NlfJ35vsKcwsuZU9i2x1OMXLuC4zhMmTKFm266iUmTJsV9NSf7j2BC8DCq3F1EaCXX7hc/psQ1sQliwU1U3bZAJ6rGE8JvG4nf542vRcKnzQsAABNEy+hZ7bYxixmzXL8QzrQomn0mrU1VRMONFI45grKysn1KqFVVVV5ELzhszm2Mm34hQ0YdGw/43R5HHHsbk2Z8l4MPv5IjT7t7n7XjQEJZWRmnHnMmtatNVLXzzjuP3IQoZho26SJGArW8oNuAIVDLaosk5rPb/rdos/6LN/3/guNsdAX1pv/tP3uCiAwWkQUiskZEVnuWe0Rkioi8KyIrReRlERNFSESGiUiziCz3Pg93U31f4FMReU1E5sY+7c7/fg/auMcyexNJNk6oMeW462uTTKHtoY+/rCxJIoREIrUcNb/tNtIAo4tVizYHCMz/KhbZIyfSr/IIGsu3csEFF/S6Cz3BvJ0PMmeIIcg5I66DoN9M1YF5W+6Nl9uxYwcrVqzgr3/9K+eccw7vvPMO9913Hysi/46XmRo8hlqtpMYxe1/lBPthZWcbKQQgFIT0tPgSVDc9gBuwcf1mCuUGLJyQheM3g+EEBSdAwkyhTbKP/cZNHEujUrGzsxj/zZupWPlvarZ+2q3vr6ry6KOPYlkW3/nOd/YY7OT+++/niiuuAKCwYAq7ytviAzQ1NSWtrisrK2PqjDPZXboM2xckEm5gyqHf77b+/41oaGhg/JDp5Et/+jGIrazhmWee4QRfcshPiXnTtL8e7T1sJObtkUCmXroeQNLq55j+R4FrVfUjEckClonIfOAx4DpVXSQi3wGuB37uHbOph+H7bu1BmUki8lE3+YKJgtUtej/9jxGqkDz9t9vcf2JO57GHPObiE5euYi4/Qpte0I7dLN4xMYJNqCOW7gQEKxQia8jYvYrs3h3m7Xqo7UdGGrRGOpQ5ceZtvP6+ce0aPHgw1dXVvPzyy2RkZqNucvT1obnT6Z8xkcKsDMZwCuqzkLCDhrzIUU2tcX9d1+9Z/oM2jt9CYyQasIiGJG6ocvxtLyFTSTKpqgI2cQ8BNxomWl+Dz5eGnZVFv0NOpHzF2x1INRwOU3zShUSbm8gsGsmGF38PwE9+8QCDvnIRy+/vPEDzvffeyzXXXMPMCZfw/upHkwg1PbuQky+8j/f/1bak+uyzz6apYTdDRh7Llg2vkZ5VSDC9T+cX5H8pjhn1Iz4qeZYsO5/R9jR2tpggMAUMZKOuYqRMapMXYtchGjUGTBFDqKpg2XEpNcnjJoFUVQRsOWAMVQpE92L6r6o78abnqlovImsw4f3GYKzxAPOB12gj1Z7WvUhEhgKjVPUNL15re+af2IOq9ri9wt7pVK2YtBpzRO/E/cdKJlU8yTM2C1ArVi72O5koXNtYtzsQq5g6Giu2kTdm3xJqh74G/RD08+ry2zhp2i1IwGcMcgmIRqMUFg0nLasfxdPPY9NHz5FTOIox6YeRlVaIpgfQ5ghOmqcfDVjGGT/NW8+vmbTm+IxOOf7ikLg+GQyJuv42N6kkQoW4vTLpofLUAZG6atb/8Xb8OXnkBux4gPDTTz8dy7IoKytj6IRphOuqENtHMKsP4YZamitKOfvss/nHP1+lsvQTmp+9j6FbV7B17uNJ/a+srOTHN9zE5BmXECwcxfSiG6mu2EDfwVOor9jCpuUvsvytexlzRAnr/m2OPfzww1m95WUyh4+BDa+RWzyFN1/aP25yXyQ+/vhjTj7yfByi1LfsZmifgxleeATlG99jtfse6XYu5c52yp3tDPWNxU8A8fvQaDSZYG3LkCmYKb/PqIIA1Lbbpv4Jx+iBRKoqhD+n9V9EhmGWrC7F6EPPwARSORcYnFB0uIh8jHHm/5mq/ptO4LlJXYqJHzASQ9YPA/FlqqrePk2fE73uuRvohFRjvqUxQ5TVjlQ9YkxUCagVI4cESdUj2thvt51KAMCKmrokZkHfj3h1eZt09drHv0jKU1WOOeYYFi1aBEB2n6FkTzqYscdMxwlBNdBYC75m8DUrEosD7VnoE6fy4cRwoK5ngBKSpPPOiDTuRdVuFZU5kecVkGWiSI245MesvecGXNeNS6gffPABx53/bSyfn7HfvgmntYnSN58l2lRHpLmBfy15HydsVB7NlTvZ9spfeOaZOZx77rls2LCBM8+5GREhGmlm47q5ZFYNp2Lrx1g+P1tW/4u8IZMpHD2bim3L2bF2Ybxp1113Hb/7w5+ortxE8fHfZeMb/01WVhb19fVJwcgPVKxevRpVZeLENsFm48aNHHzQdIbmTSc9vYCxGSeSkdkfJ91PBbu9Ui5FuZMY0VpMwMoAx0XSQojPZ/YLA/Pt98Wn/kanbset/eqzkvyRzbenXvsC+t4TKBB1O5VU9xikGkBEMoHngau8sH3fAX4vIjdjIlPFYqDuBIaoaqWIHAy8JCITVL3tQpLxA0z8gKUAqrpBRPp1Uu5zo/ekGrTj5KgJ0/1k53SS3qSJlvxEUnVtiZvK4gTqSybW9mQS8wzQgE1TcyXjb7qXT+/ofFq6v7B7926+973vxQm16NxvkTF6AvVBI1E6IXN7R7KNlV4cQWILYBzx3KW8/tht0jt4K6OisRVS7U6cQKLi0lZnIqm2I2In6iA+P44YNUZMN7pt2zYOPepocicfQt9ZxyDBDPxuBs1lpaZPJ51H8+5SKj54k5byHWjUHH/JrXeyfft2rrnmmqSmFfQJcdl3z+Cii56nqKiIkpISnnnmGbZt20ZW1uQkz4m+ffvy/uI3mTx5MsNOuTie3tDQcECT6ubNmzn6sK+zbZdZyluYM56stP7c+9i1nHPOOQzLm8HogcfjZBr9cUtugGiGjZXXB9+uEOkZhYye+Q2Cy4xAJMGgIdBgEPX2CtOAz0zlYyTqzQzdBD9vwLw025HrgUOq0hWpVuxpmaq3N9XzwJOq+gKAqq4FTvTyRwOneumtmH2tUNVlXpDq0cCHnVTd6m0WGDuPj/00ZL0m1Wi67elI20mZdjspVLzpe4LxJGn5ZEz6ijm3e9Nf9SWoCHzJRCzeNNcVyD3oEHbNf5GcaTP2tu97Bdd1KSw0e5pnn3wkfS86iWiDn1Y8B35LwectZFDAFSQq4MTcnWgbHBLGx4PEjvHI2BRqR7KuR6qJ46LeWLYj1cYdnyE+H5qWvKz2wQcfxE7LoODE00yVjiHpvKmHUbdpNXZmFs0rSwgVDqal3ETpGPqtK9gx9wlu/e8nABh51LfYtOgvjPreT1n/6B1J9Q8ePLhbz4LXX3+dUF4hmaMmkjFoJGLZ8XE9EKGqjBgxAstqU//47RAbdi1g0aLptLa2MnTkoTiZISI5ZqybCmwi0kqFW0rUaSEtux+1IwJki4kHq7bR2ahfEp4TAdUOL9RENRnt7hmTcQCRqkJ4LwxVnpP9fwNrVPW3Cen9VLVMjHPvzzDT9lj81SpVdbzwgKOArkL/LRKRnwJpInICcDnwcjdtGYTRvy4Qs2urT1Ube9KPXpNqJF3aphuJhNnOct+xlV75REk1USr1tZFsnFTtdseoIRKxIH3ceNI3rWLXolcY/cts1v/0i5FWE407/S4+llColRbbPAGuY5Y5SYwwVVBHcKMWOJ7jvmsc9ZPgESlgFtq4pu9JJKrmG8zx4kAs1kusXGzBBMTGUkkbPxa3pZnP7vw5ufP+gSBEmxpp2rSWrCkHs2Pe00Qb6tFwhEhFGeIPEKksZ+39P23rsz9A8fdvwpeVRf8TvkLps38CIDokh8zREymd9xR9j/mISEMNoeEjiJSX0VKyldL33yMvr/NFKEceeSTOLbez6e+/w2ltJpi3X2Zi+wwiwjPPPMN5550XT+tfdDAlVR9xzz33MKz/YayvfYdh401+7SBhy4t/pHHrerJnHIJVs5uCCbNpKgAnYEg3ft0S+cdT/4hnDhEXpJMdi2OzE3EhttjjgCHVriXVPWE28E1gpbdBH8BPgVEiEvNtfAGIKfaPBG4TkShmldR/JUbyb4efYHZfXQlchtmS5dHOCnrqhh9itlwZCQwFHsIEzd4jei+ppnnuUu1cnRJ1gJ2iPfHG/FATpdKYCiAmvdqx47ylnK5gRdvO0/fUM9j24G+pXryQYb+LsOXK/b90ddxzPyfvzEOJ7ignIyOKbblJJApmmR6A6wqOY3l86uI6npitCVYlVwyRxk4gEl9yGou7IZ4La7yIS1xaNQVjGQkPqGWCq1hYDL7mBpo3b6LixeSNGFt3bOfXN1zHoEGDsG2b0aNHU1tby6JFi7j++uspvu03hHeUImJj9cnCBdLHTaDfnLMpe/VFtv71fnKmHULG+Ik0bVqHlZFO2fN/N5HKgNLS0i5JdcaMGTRUlbFgwQKam5uZPXt2p+UOJJx77rkceeSRlJSUkJ1zLDlpg5k64mss/+xpWsJ11PW32OD/iIoXnkdbW7Ey0sg961R8rT40EqZhZj+iOUo0zZuCOgnGWIirdaxIgmon9hKOkWhM9ZPwwo2R8AFDqgqRvZBUVXUxncvhAL/rpPzzGFVBT/AjL7hKnEhF5EovrT2uIFn/ur43+tfeS6qZCX1ufxWlXXoi0ca4JGG67/oNmcZ+G0LVJGkWMFNqDMnElmQiCsEgBRd9k7K//IWqt9+g/6a1pE0ax5ZLru9tt3qEMfdezJZ75mL3yWTo5SeQGTQrovyWJ6mq4KjgeG/pSNTGVUW8fAsLFUVdQWPSqkV8Hyyv120n9AZPaVMLgJmK4rblx9Us0jZ2gFE7AOHqMmoXL+Kuu+7i6KOPZufOnUydOrXLKFOxbbxbqnbiH1SA5QvhRk1DNOqQMWkyhTnZVL/9BsGigeQcNJOcg2YSiTRSv3wZ6ZMn0bR6NUOHDu20/hgCgQAnnXRSt2UONIwdO5aamhoyho2hYpCPxqZhBKryaRydQWTzZspXrIiXdRubaV6+Cl+/vvS/8jJaB7qorfHrgiPJOzMoWBGjLmqvLyfhBStOAul6/1sHEKlC2zNwAOEiOhLztztJA2hpp39NUKztGXshqWKuXqL1GTq9oknO+5I8tXf9MVKNBVT2JNJEzrY0fqxJUMSW+JQHIDBsIEW33UDzp2uo+NPT9PvRxfS/5kJKf/U4fn+y+9PnwV//+lc23fo0439wOEUnjiXocwm7rbREfVjeQDiuhUZtojHDUXxlobe8FBfDem2DpR5jxsO2dqEvM99tOlZNkF46e5jEMSqa1u3bqXjuOQrO+xo//vGPexTs+5prrmH79u289cpLbCotIeuwQ8ibcypubT07HniATJQGv4/MqdOQ/rmEnQaccAuR6goG3fYzqv4xl6yjDiMnJ2eP5/rfhj/84Q/MmjWL3ev+TVbRMVjD+jPsmpuot3ZQ/7Ol8XJZJx1JaMJg0qZMQETMbEYiJvJfXFcu3gvW/BRHcCwL8Xc93cdJ1rdbThuxHiikqnDAkKqIXAB8HeN6lbiCKhvoKpzfEhH5MRASkWMwngOv9PSce2Gooo1QE63Rib8hTqKJkqqRTj0S9RupNEaqWOqRcIJ1q/0GS+BNl6UtzowAfkifPobG90ay6y6zVv/O3NHcfHPPY68OeegnZOS4WAEfO599nzRfM8EQMGwooaEFrL/iMQ6/93QGjM8BWmiKBmiMBOJTfYhJqhaJMYItUWL3lwWouEYiTygjogmLYxIJN8GDQkG9isyq1gRpV72B0OQHVlyh+dM1ZE6eSsaYcT3ePSEzM5M//vGPgPFDPfXUU/ng5luws7Ow83OpXL+Rf/7zn8ydO5cn/v02ZVueQCxBxUKbmkCEdWvXdlm/qjJ37lx27NjB17/+9f9V5GtZFieeeCK//v3DTG46Om4/CA7NI2PWRMLbyym68wpsH9g+h/bLyV23bfWTev/HXAON/cD1LmfyjCU+3Xc8SdZ7Y7vtSPZAgGqXLlVfBt7BuF71xQSgjqEeWNHpEfBjjE/rWuBKzGKDP/b0hHshqWpcxyMJU5IOhGqB60uYynuGE42Rqs9FfG3irlhtxJI0A26vqFXM3RefCoPYitgufS44iaYPVwJwyy238J3vfIdBgwZ125+NGzdy1VVXsX3hW2gkAq6iUYf0b8ymubSKqlfXoFEHX246kZ2VMN4QgNuJAjnW1phk6hfHkF87XWus3WCMvVaCcaszUlUVHFdwvRvVVcF1JE6yrmMsxuokOCtGhciOchpWreSXV1/Ntdcmu0D1FPn5+bz33nvs2rWLZ599lpNPPhnLsjj99NM5/fTTeRSzcaLrupSWljJmzBjGjxvH/PnzGT16dKd1RqNRvvKVr5Cens7KlSt56KGHOi13oKK4uBg30hr/7RQ3EwxEGXzDV71rHI3PXtyEaxh7+cVDqYoiliaogmLTzYR7wPUMFk7M00aNy1XsGFfjbnqZgWQPjy8P5n49EKCqW4GtwKEi0h+IuQutUdUOq6O8qf6fVPUi4A97c85eb6cy4p7f7DdS7aCTpXekClDxh6doeGdVvPjq1au7DKi8ZcsWZsyYQeapo+C4OWRkR2heW0rYzmTgpAwyAq1UNmWw7vI/kjZqAMGGSsZ+cwppeWnYgwZQGw4lEWXUsQk7dsKDxBdIqqCOhXr6uqr/+QeN7y+jz8knUfHiS/t1j6+29mrcD9YKBQg3NJmtZDrBGWecQSgUYsmSJRx55JH84he/YNmyZYwZMyYpBOGBhpaWFqZMmcKuwUXknXkaVoGnVw9EEdH4Ne6KVB2nzW1GvedG3YQy7IFUXU8qjZMqHqma35uvuvZL304lrbhIi3/bMTzmqjNv+8K3U4lBRM4Ffg0sxDDNEcD1qvpcJ2VfB05V1Y5r1HuAXkuqTpp2cPGIO58n6U8VfGq+aSM+25vuW7aLWBq/+ZLckLxqYlPcRAICbwoVv8kEBHwBM80q+OphtK7fSqTC7IE+YcIEKisrO7VCb9iwgZEjR5L2rUMJOw2k+8MwK4eA5QCNBGyH7EArzldG89mzK6itaGTXUrPH+tff+x5h16Y52qa3tS2XYBfT96TxEk16gBJJNfHYmF7KcY1aIeq5ZUUdiyh2PEaruBaua6GOIFGLmn/Op3XdBobcfjNbrv8pXxREhIULF3L00UcTGJjP4EtPpuCsQ/jktNs7lH3ssccYPnw4TU1NZGZmMm3aNJqamgCjt5wyZQp5eXnU1dXR0NBARkbGPo/zsDe44447KM2CAd+YgWQ2Y3svc8szRvqstuv4/9l7zyg5irP9+1fV3TOzSbur1SoHlJGEUCJLJoloMiYZMMkkkU14eIQBG2PAYIyJxn5ssMEmGEwUJtlkEAJMkLAQWQhlaRU2zc5Mh3o/VMeZWWlXJP3Pee9zdLTT09Nd1V111XXHcl0ZS3rS0RiG8MIxr4gz2CSoqvjC7MnQ7ho6OYOx5Qm8AGg3D3IIEI7VzUguBbZVSq0CghjXfwMloIqOdX1FCPEYEMamKqVu7sqNur9FdW2h5JgSRXZBqZBSIY1oAJmGh2F4GP7gM0QEqIEUe8+VAs+TZVd/141YGxBet2pUbyb/bQZvH34TTmuedGM1w7YZzq437M2jB9cfTSIAACAASURBVN+fuN+iRYsYOHAga4CMaftgCinDRQpFQ6odD0H9D4fR/uEyrKp+fPzIxwybPpgaM4fMKNqcNAXXrya1iaNaohK/dTyJ40lcGYCqxI6Z5pQSeJ5KxgMrgShInBXraH3+VfqePeNbBdRAdtllF2qnjsFe18aKe16iesIWZc/r3bs3J510Erfeeivvv/8+Q4cOZf78+QA8+uijzLz5ZxRacqRrUpgZi9bF65lw1o4M2Xskf99xQxXevhlRSjHknP1Yfs/LOOvaWXbNPQz9dcTGTMPFkArTH4d5x8DzRGgKAn/RNR3SluNfU2C7RkgaPARB/f5gPJTTUoJ5EX72BJ4rSxbm70oUoqx57DsWGQCqL2vovPTpanThlkr/X7ek26Das6ENiQrV1kACUIOIfdlNLaRr0xhp7SGXRUDqqejhez6YerFwDC+mOgfXhWigAbhKJuw3lnRJmR7T/ngkLx51N/mmNsYcsS3PnfUkB2VO5LF9/hyeu3DhQiZMmMAcczWmjC0AwiNlOPRMtTP375/w9NXant1/Ui+G7jqAwds1YgmPGitH2nBw/TY6npGYGJAEWt1fmXRukTQNOJ4EGfTdjyooegfCZ7YiiN8NHHxSsfaxx7nkJ+dz5ZWl7PDbkhduuYfJkyfT74c78cUv7mfnkVW8vG/prqu33HILW2yxBa2trQwZMoS99tqLCy+8kGeffZYJv96D2lGNZAwNQKvfXcLrP3uBNfNXwoNdb8vM/xzIgqcX8+m7bUw7fzJWhcktk+/p0m/HPvpzcgtX0Pb8W7R/ugo8lz5HTGXpH56lYbsh1FV20GFrTaW+IkvKcDF9m1jWsWjNZxLXq0nn6JnuoNLUxMRRkoJrhtlH8bkQHxeeioLpQ1OS/73rST1vlA7M+1oqgnxViZk0NiN5WgjxDHCf//lI4MlyJyqlulUBq1i6DaoDa5o1AJEESFmSqA4f/OM9Vi5ppm1JC1Mu3Im64Q14KmYXLFrRggHkxQZOMYsrPicA4uCzKTQbqBxYzb6zToD1zTSOqmf1/DV88uD7sE/Uvh49erBgwQIy0ysxhYenBCveWMzSR9+lbVWWyYcP45mrIwfhsnebAFj06jLMtWsZf8LWVFakQ1DNe2bJCu3E0mVsz8BTXqK9ARAnJlTwLzaxik3fUqoQdDF1FHluaRMdc+dz9jPPlbyLb1MmTpyot5558I9IU5JdvJbdnr+AF3b/Tcm5xamsq0/oSdszHax4byXWsP7hGOs3pR+7/2wHnjz7eaafvxW5FptdzxnHVVs/3Gk7Jh02lPce+iL83HN4LROOGt3p+Uop7rzzTs6aeR65Ne1ktuhLfmkTtZOH0njI9vTaaQRuNs/qR+Yw/vARNPRYw+K2OgCG16whLZ2w4n3eM2kyq8MSeKbw6JVpoyGVxfSDTB1lkHfN8DeOMnD9hTcQTwkcJUuOFYNuAKqbgygINcnuiBBiEHrPqL7oyNz/U0rdJISYgE5NrQa+AI4JiqYIIWaiM6Vc4Byl1DNl26TURUKIQ9H7Tgn/2o900o5/USZCTSm1V1f60W1QHVC5HqOLL2/qKaNo+qyZF342m6ePfZj9/7A3fSZvOL87Pqi8olU5OBZne8XALIUiJR1M4VHdaFDRt4a0zJEyXVa9tpDfzN+TlY6eCBeccCvjx49n1LTd2XJiho68ZO6NrzLy4NEsvPlNMimPn8/bHyH0lrv59R38do9/4RQ8/v2HhSx4eTVLF7Qx9Hv9GLJjX8YdviUOkY21w7VKTRyIsA/JPpXvT9AnIcAIIsANDaqBduB6EsPwcBsrkDWVLFmyhN69v7u0TyEEhx56KL+6/QYGHDyJ/172KJNuPJIfzjmV+3YoKUoUyiGvncnbFz5CzcjeDDl4PCnDCU0yFYZN72mNpK/fjo+eW8YHTy1hl7PLOyAve3s//n3jAuY/uZhTH9iFQWOqmP9GOw+d/Sr9RlbxU/PQsmD8pz/9iRtuuIH+00fx+f3vUFFvstPtp+GZKWzPwJI2VEr2ePB4qq0CvdLrwxJ3jalW0tKhzdVb4uQ9EydlhO+/wijQYLVTaUTms5xnYQo3XHhdFY0Nt0iDCedEwFCLxogh1GYDql+BqXarSLUQYixwFDAO6A/8WwgxqnhbFCHEwcAI4H2lVFfCYC6N/Z0BfoBfuKUr0m1QrTezJcfK2kbRwFc9KsNhN+7IfSe9SI9KhzpL7/YqfYAoMQcgE+p0OXFLGK4MmbIUirSMIiUs4WIJl4MvGsGDP3mT2098mzGHjmTLfQbRp08frr/+embMmEHzsGEsXryYnXbYgdN3O51Dbj6ELaf3S3jNq+tTDN2mnk9mr+F7xw3mlbu/BGDhK8tZ+MpyPnl+KUv/s4rDHzyIui1qMaWbaJunpDYNCAhTZHSeVWiLC/MDAOkzUYkG0GKwDVRB2486sHoLarcfxfLlnW7D863J2LFj6aEyVDWkGbLHMJbf/SJDLt+10/Pb2tp45+rnaf1kNbvddSTVGZeUjEA1JR0s4TJx7z7s9P06Ln9tBVbLupLrtLS08Pfz/oOZlpz82D409LVw8BixXT07HDeCV26ax4h7dyvbhnfeeYedd96Zh97Ucd59dhmFTJmYwqXStOmVace1Xf591lO0r8py7pPTaUi3AVBpFEhLO9ya2RIuVsYLx3mlUaBSakCNb99s4GH448MVMpwDgbiqlDTExRDBGPHIyIpOn++3K2KTQHUTilQfBNzvV6taKIT4FJ1e+nrYEiF+hwbd2cCVQojtlFIbtI0ppd4oOvSSEOKlrvbjW9lHtqZPBSf/c79v41adSl3/Ss65dzvenrWcF/7yMXPu/Ig5e45j3H6D2P3ScezEIZx88sk0NjaGG/i1NeXJrsnRa4uoHN2k/frjOYqJ3+9Lukea1UvykE4xdOf+ZFtc6kc1UNnruxvcuUUrv9ZMsk0V0zR5/vnnmTx1W/b88w94+ZxZfPnvz2HH5HnNzc1ccskl/PGuP1E1uI7d7z6SqgG1bCw/aNDWtSya26wDY3yZO3cuu+y9PcOn9WH/n46nQ1YQL9Q+9ZRRvPKHj3jznk/Z6vKhvP/I54lF85e//CUTJ06k3yFD2eqi6Zj9SiNGvnxxESvf0/6O5uVZ6L35lCrMe93anv6bk3gSSlK6VE8V6GqR6gHAnNjPlvjH4rIzMMGvZFUJvAJsEFSDPbB8kcAUCPaW37h0G1Tj6kt41xgTA+08CsQosrWWs70G4imJq2RMxUmaAYpX6bhYvp1Kh6x4ieMGirS0MSyP3Q7vxfcOa+TTt5p5cVYzfz32RWr6VvDq+l/w9Kd/ZO/zR/PwbL1NzV0/no3d4TBiam/2v25HPCUYdcBwRh0wnLxnstXooUCSUY/aa4ivurlI5eEoI3JcCW1PTZo4PMzYnj6mf048I2Vj2SlBaFfh82V0fLyMbbf9dsshdiZbbbUVRxzwA1749QvUj6hnyfOfh5tgHPvaCXz86Md8/PACevStYpdb9qf3qFqkKYE8UigqjEKodVQYNrVmlhqZQwrF+Gl1vPXgl9yx51R+PPo13njjDXbddVd2P20Q2/zYzx5z9dixlYWnBLmsR49+lfz3n0tZMX8ty5cvp3///mF7GxoaePLJJ9l///2pGt6b2gE9AUWvTDsp6dCQyjJ37mIAqntXUN/bImdorbBSFrCEg2foot5SeNjKCE1lUngYKHLKDMe1JVydaeiPd6kC/0Gkvbp+pkBo+qHIQRxoaJuL6h9Iee//RuupAt0pUl3uJsUPohCYA5RSWdG1gO35/nUEelVeCJTfl76MdB9UZQSqxeAViB5Mpep9XEJHU1FUQ9KeFP0NJGy5wb0NFFJ4EajG1CnwVWcRqWEAHpLx21cxcrs6uHIIACvXSP5y5jtct8eLjN65kakH9+LEXwxli1XXcMopp2AJlw4vFUYaBLtDekpghaaHYID795aKvKcS7XYDAA0WomAB8cGj2J4cP6ezqIKCa1Jp2qxv0jG09fWbz35Pd955J5dffjnXXnstVk2KvX+zK4vmrGLhs5/TMKYXE06dwoid+2JYBuBi+uMrLR3S0gnHW6WRp1IWqJJ5bGWwx7F9eP7Pi1m5sANGg23bTJw4kWmn1tPmmKD0NQzhYeHiKkm6TvKTZ3X1tr+d/jp33XUX//rod+x68nBqGlL8ZMyzjB8/nhkzZvDH1+6jcbImPb1SbbotRoHtDx/A3Ac/o21VB1dOeZLLX92dqvqUXrTxSEsdL27gJeZE8F4NFIZMxpQnxnuZvZ0S47mIlMTH1mbjb++cqW5UulOkGs1M41urDASWFV1ySyFE4G0WwHD/s9CXVluXacaw4sB/v6h1l2QTQDWPJVxyyiIjkoMj4w8WWxkJEIPI9uOGDpkkK41LMLACu2w5CcA6GLyWX4DSKALx4DjEmDQKWQS+Q3pJfn7/liz7tIMPnlvJvb9tQgjoPfk8xo0bB0DejT0u2+WL+W28fOmLOFmbsUdvxcTjxmKaIhz4FYbtT+wYqCpRdrIVs/34sXJ9dzwjWnB8bX/EHnW8+jPBzEe3YdCYGs4Y/ULZ33+bkkql+NWvfsUHH3zArFmzeP23bzPiwNFse/4OCLvAug9Wkt6zL5bQ4FntM7+0dMhIm0qZZIIZaetnlIbxO1Xz/r9Wwj56N4ZmsQDYKbx3jZHDRWB7ZriIpaVNRjiM37knl1yi43hfuGsJ5/9jexijf/dC2wM47QUqTT2ee6dasYRLWtqoBv3uavumaV6R5xfTnueaN3fB6KG0VuTfp6BMwoh9CBEvHXunAfEIiIVXBlCLxRIOLrKEtBibHCX9Dckm2FS7W6QazVrvFULcgHZUjQSKt5AeswmtfwMoTut7s8yxsrJJoNrDyNHupcMXm5E2Oc8KJ4CtDB3cjxca3eN/gwYXDZrJgeQWnQOEv3NjHtAAmALwDoArBNcYqIeAHgxEJAYSN8580UVNB4/KMGrL/mx7+CB+stNb8EgTu+3WiFTVeLZN2lK4BY+/HPBPXNujo1lPvHdvf5sFDy6gYas+KFex7JUvADji/3Zm6PZ6n3dPCUypEm0LnBbxyRG2KTgmvBL1TpsRZPh9ABaHXzKCW06Yy49vGsedYhonjXqVzUEuv/xyZs2aRe2AKj584ANEysSqTtP8SRO7nbcVlvSokAVqfPW5UubJSJu0v3BnpE1KuOF7AjjunAbOPegzHjlhNPffvzUNA3VcaGCiqpR5bGWGsa6gx4clXAIyf/oNI8m1udx93nuwYhiDnXN55bb5bH/pLpETShaolHlyLQV+MX0OB5w3lFk3Lgyv2bpoHb227qHV//j4jf0dvL+4JhUeCxxV8YiQGETGwwrTwk7co1LmyXrpxEL9nUsshb2b0q0i1Uqp+UKIB4AP0Gr6mcWefz/3v0vi10zth94dYDwR+e9BN5IAug2qPYwcGVHAig3UjLA1i/AnQDA4UsJNqO9xKcdS3ZC9BiBahsUWAXEwcOOMFEoZawDiyXslWXMgnifp0SDZ55SBPP3HJbzwwgvM3a+CtUu0I2D4Tr3oNSDNBfdP4afbv01DQwNz587lxRdfZOHChSxevJhPen7CY489xnuPLGLYlHptK/S9u5ZU0cRSyYkVmit8G1wgxZ/jEpg/UsJh32N70bef5M7zP+DyRyboHXs2A9lmm214//33+eyzz9h6660ZOnQot912G7NnzyYtC0ihqDWz4RiqMXIJDSTl/23ghe9x4BCTsZMruHHmcl7+58cATDh0HWMnaXDNCJucUkULrF6g1i/T77J1RQd7/bg/hZYcN5y1kMZhl3L03/Yg1a+e3imd6uwiKHS4nLed9ons+sO+PPunL8m36fnbMXs6cus3MYTCILbYFw1fD5nQ4PT5sUVUej7oRu+5oAw9j4LMQaGwhIPtFyLOiELi+81DRFiroDvS3SLV/m+uAq4q990myH7ASWgzQrzKTyvd2BK72wVVnv18S1JFq2LwkkuADVUCqsnVOA5ykXMqDrbFNtdy13CVSAxEiUcq1pZyYO0iyXlW4jrx34MewJ8vyHHZgckKYdufOIJpJ2/J9TvNoqlJJwT06tWrpJ1z5szhwgsvZFHLfzntbzvQRjWOMqg2cqHtLe+3Ic5aLOEihZcwr8jYBAwAOM7Wg3cSvIObLlqMlHDKZf04esL7JW37rsXzPPqOrGW/mePos91ALOHSx2oO2XmdkcUSTqJflg9YcRvzp/Nz/O9xi9lq+0rmvtbO8J0aufhWbSe3lYGrZNkF96azP+HNp3VI1s0vT6S6p8Wzf13JiEO2pCNTR4Us0MvSoJoRNl+8t557L/+IX/19OGu9Gs6aEkXdNPSz+O1LkxMgGVfRQS988fA6IPY+o7AoQ3gJe6ytDCzhJrS0jHAo+BpOMBYCjWefYQu+84Iq6cGDVL+Lzys5vuisC7+zgirdESHEEUqpBzb1991mqlW+ip9kUQqDpH1VhgMlGYdq+IPFRWjWEdgU0cBmxADbRYRMLi4JmyQCVyQHcPycENhFzHwghLZ5+ThqKxOD5Gqf8oFt/DjBo59N4D8vtOAJyex5aVZ+2Ey6WoNhOTANZIcdduC5555j1JQGbj/sFVL1lTSMrGPnI/owaKTWJmzMRMSCJVwywkYKFYJqSQSFD6KRs84L7XnBczjjsl7cetlKLj7sE9pv3oLhY9IsWW3RtDRHn0EpWlsUT9+zhv+81EbLGpuRk6o599YRHDfmrU7783XKXXfdRabGZPgODXR4gV08skkafh8DQLRwQ80nHr08YlyG+97ekteeamH2061YFtRIzUKzXrqE96xdUWDeqy384OSejByX5p7frODSg97nrN+PY8+TBrHctsAr0MdqDp1kq77s4I4LP2T1kgIvPLSW7x1hsdtxA3jl/mU4BcX6Jt1GWxnh/QJWndCYAvOFL8HCX2wGMHybvz6m/w7Gs4fU46OIfRtB5fjNRDZR/f9GxC/nd5dS6tiunK+UekAIsTc6vjUTO351V37fbVANPN0pXAqUD86XQoWpdwH4dsZYO2Or+rtkHn38WAg0SmIJO2QxG2pz0ARXSTLCIeWXUywoEw8Zrv6ekqHTLZjUO+5eia1MVjZl+fLFZoZnVtEVSafTfPruOp588kmEELz33ntcd/yVnHT1CCZObyCnTGpkLjw/I+2wf1YR+zSKJlogKeEknwnQWA+X3dyH5x5ZzwVHL2bg0BRffmbT2N9i5eIChinY66ie/OT6gfTqY3D9eYu55tgP+PH7kpNm9uUPVxQ7Ub8++eSTTzjnwtM4/vfbYkiFVLrtLjJ8j5qpF8gEoCqccIGMj50cBjnPYvJePbn2kQoaR9ZSJ9eG3xc7JOf/dz2/n7kU0xIcfHofqmokJ87sxwPXLuTS+7fy7acaVFPC5fP5Wa4+YSEHnTmQp/60jD/+cgXTj+vHuT/rxfFn1vL0rBw7H9EXQzjYnqSvud5vr1tirpFFmgYE0SvJMoG6/xsGyEwR+w4cfZuLlGxu+R2KH6PaKIRIKaU2+qD8hIE6dIzrn9EZVXM2+KOYdBtUgxdfEEYInKCBKvR8egY2BhYutg+agW0oAK7oe1+N8W1lrpKJrJRitTxQh1J+UHe7l6ZK5kNGENwnGMBeTA2LAH7Dom2UEUOKOu/y8r1L6dVo0mC2dfmZWZbFQQcdBOg6omPHjuWci49jyh71sTjGWPZUmWW+mImXk/g7cJVECNjj0HomTavhww8cRm/Xg+oKhesoHE+Q9msaG8LjklsH8taLbbhH1vP7y5fxizNW0qdPny73savy+eefM23PrTng3KEMGJes+J+wfZbprxTKjx6MIiUsXGxhgoDRW1cQhTBCpciX2O132quGo87qxUP/t4Z/37ua6lrJe/8czqovXkPYNpbU3nWnrcDvr1zKW8+1MOOXA7nk2Peov6qOax8fS6WhHVI9exl879jBuK7CVip0gkGSWW9IisFzY2C6IemGJe+bFQVdGK7ftnyB3iblcZLl/G4oc+40pdTWQoi5SqnLhBDX0fUNBrsPquu9itBeFNg7A6N5AGSB9x+gxdXsOVBXm12t9lo+8GW9dPgZNBMN1K52T2CrCMDW+78dYK0L/17h1DHIWqPVeWClU0t/a13IZnLKoq/RTE5ZIeC0e2kM4SVsqkEbdVscqkQhwQYLSnLt/6zhk7lZmpqaaGho6O6jC+XQQw/lRyccyVuPr2DkfpX+YqK/iyc6BKaWwGYcAL0+N4oFjkdYFGuAhvDo0dtgu95pwNW/8bsdeJUlil4Ngn1/UAPtHdw6U30jW5yccu0I/n7tIg44azA7HjGAtY5+Z7VmljojS6UskPEZV6ApFIfmxRma6//dQ0VMPw5iVSKfsKXr3whOOr8XhXaHN19sp2dvg9deew3HcVjw6hpG7prmP7OWc96FnzN4VIZ5734abpCoPLj32i9p7AVISXu7YuEXX9A4rIr3n17BX94YS0PfbNgOS3ixeGzf0VSGmYbfidIwwjgzh1gkTEyrM/xEg80GVNm81H9flvn/JFCzkXODAZUTQvRFlwncoqs36jao5lQqnORxdRmIgawRA1gz8X8AfjnPCuNdgRBcg++C6xl44OhUwGa3MnRerHZqwt+tFj3Ieim/ff51i67R6lUkjP8ZYYftD2x3wYSuFAUy0iGFG06A9+favPzIWt5///2vBKigC47k2j1uu3Ahe8/Lsu1RW1C7hb9NSwAAKpkjXhy/WpytA4HKGamYEpUwY6SEi+H3MTCrFKuRH31SoO8Ak0wmWbbuq8pvfvMbZt2+lP+5eyz9R9Ww1o0WuRqZCwP7g3cQz5DTnyN7ccjmVHA8iSZxs5SrZEk4X0q4nPnTXrz2TBsHHVPH3ofDry9cwaQpJqbRTtOCtfz4ggb++OvVYRprUED7y486mPtapL0ceEwtq5drh9ayZYqPnm1j5XKXPQ+pYcwoGSWGxEA1EEP49n7/0IZYbRxM4+PCRYYOuc3Fpqp3fN181H8ApdQVAEKIKqVU+0ZOf1IIUYfeKeA9tHJ7V1fvtQk2VUfHecaylGw0Mw3DPvzLSuEl1ICMtEPd2xOCGpljjVvtXyMIgbFCcA4cMAFg1hjaAbHGrQ7ZrBSKVi8TOnUy0qbVzURggkt7yIaT8YrxdmVEITGRU7hYwsPyJ/Csv7dRVS3CRICvIkIIvvzyS4455hhWzG3htqffobJniiMuG0n92D5Up4OQquTCFLETQ1fq8q8nERgoPDwMERwL1NFkSFKxUyQulnBY+FGBPv0NlFJf2xYsd999N1dcdTFXPbYVDf3S2EpHNtSYmhBolprXZhyieGMDlXA4BgAVApPQ2zoZMRNNAUmGIDtNhLVpCyHe+GPLEFz1x/5c8KOl3PToFgA89+dlHH5OP3bdI83Fp67izOP7Mmn8SgAuvfRSttulgjdfivLr9zmiB2dd2Y+FH+XxMs2890aWB36lbarPP5Xl5At68vi9LYwYbbLvodVM2EpiGOFWVHgKOjo8Olo9GnsbFJAYSuHGnrtEA3NgarNIRtUE2WJGfIvhzUA2N6YqhNgRnVhQDQz2ywmeppQ6o+g8CTyllFoPPCiEeAKoUEqt7eq9ug2qddJXb5QMHT9x1Tr4zkOHLQVgFzgNPBkxBks4NBhatQ9sowVlJIL9DeFhxeJ5q2SeNW51CJo1ooO1TvS5zsiy3q0MPcBVMk+LV+GzIDtx3XjwfHH4VwE/Z98fpxUVgvY2RTabpaqqqruPrUQGDRrEyy/rwjv5fJ7q6mpe+o1k7rxnsCpM9p0xhOknDkQIEfYlYPnFmThBQHkApLpPOrwnZHexv4NnoK8RmT0MFPseUsVDd7dw6vk9Ofm8OrYfEgW5b4rMnTuXH598Ar/+x3AGDwDQ79gSbmjeqJEdPmg6sbz4ILQoAtKImetre0r/i783W8mYxhSF0hWzQAPFiLFp9vhBLTf/dAV/eWEYV51qUFmzkqNn1DN6S5PtJq9ixOgUI8emmPVgO1feUJ8A1f+5uIp/zWrhmnO1U++NJ9cyYKBk6RKPLxa63HzVWv7wu3/w5ptvctVFj7C+9RP23r+CY46pYIuhJi6CJ/7Rzs9/2sLBP8hw3Y31eOiFIWK0vsMrcLIigHyiNKBLEIO9mYDqJgb/b6Ce6kR0FlUGHeR/hlLqTSHErugiK8EgfVgp9YtOLn8jsDc6Cwul1FwhxM4lTVfKE0LcBOzgf+4AulWpptug2uB7qvWLDF5sstSg6zMnoCRCIOdFt0wJl3alWWhj8FslQ5NASjihZx4iJ1Qd2YSa2NNsC4EhIwr0NKI41fjkNUKWph1m8RCWuO3K9Qe1DlzRv3n15TyXXVX7tQBqsaTTaWxbA77jOCxcuJDvHz6JD2evZZv9+jD1kN5hVIBuf2kGTjxaoJzYytCLHzFTgb+wxBfEdEZy81/7cOphKxg60mT7Mzq7YtfkkksuYf+zBjNqfJqMSI7NMGQsGCv++9XH9ELnxYLeIQLUQOLjMJCQ4cdszfFoZ5vIPHXSRb254tTFXH/RcvoMN7BS+sw/3deLBQscVixzeeqJHKPHWeyxb5rLYtU4L7tgHedf0ZPjz+jB0OEmny0oMGX7FH/9Uytvv2GzarnLEy+fyAGn9WPaKQb9m//DAw88wJEHX8tvf19PNi946y39DB59KEff/i2cfXY1mViRM08l+xyMASvG1iPZTEAV2MQEr87qqV4HXKGUekoI8X3/867+b15RSu3flYsrpRYXaV+dtfJfQoiDlFKPbUonNkH9V/4AVRHvKTPQA6ksarfng2FwTqUf3xpfcW0VgWAhblf071wpkiCeKlJr42FYroqcaqENmKjeq25LcjDmlBWym9AmZgpeebmAmvn1qcXlxDRNRo4cybzXV3Pfffdx7rnnMm6nWvr2S3r2y4WbRbY1LVYsSFx/b8bqM3hhQHqwiGXMdRgo9+G4AAAAIABJREFUGhoNzr20ntt+tZ7dv78FO2zxxSb1JZvN8uSTT3LHzNE0yPaSgPbANhgswvFN8eLHIXDwlCaUlITqIUNVuVxmn2a8kR22ypJc+6c+3HFzM/f9bg2vPg2zXykwfIRJoaDYa/9Kfn1rT138HI/Tz65m+QqPxYscDj6ikuFbSH5ycVApTjtP63oanPCD1fzhwT7cfG0rZo91HHh8AxMnTmTixImMGTOG4444LtHu2+9pYMYxa/j303muuKYH225jYRjCd16Vsuzy/d+M7JibgO8bqKeq0KmiALWUFk3piiwWQuwEKCFECjgHWNDJuWcBtUKIPJqlBsVXSmtBlpFuZ1R9tlhX7tfB+kkJpnNwPD4pyokuNh39HRSfDm2HnXhCk7yjnIe0NCOrs0wu3e4ok0V/VmH4VxDy1bK0nQt/tJTtz57EQ+cX17D95mTq1KlMObmNKdOqwj7YyijJ2IFk9plXVLilXBZWUJyk1dNOqWGpVeH7cmzF8Xt8wZW3NHLigZseszpu3DgOPbOZQ/aXCe3FwKPdi/apD1hzwFSlSIJfRrgloOoprSrHx0hOmbHFM0phjoeueUqE4B0H4N9fu4Z7f7++bD9+fl0d779ToGeDZF0LZNs8jj2lmuw6l9mv5DnujFpqanUfOzo85s+12W7HFFsP1pXDbrm7FwO3rOLn5zWRbXH45L95xmxTyeWXV7LVOBPLEvzntRy/ub6NxV+6NK32+MmFVZx+ZjXCjL9XwucQSPA8lILRg5d/5xlVmf6D1BanlhbY/+iK87ucUeXXU30Z2AoNrM+gwU0COymlFvnq/0PoalXL0LsDzO/ker3Qqa57+Nd5Fr39Somt1E8WKJHiugKdSbeZqh2G/ETHAruWFWMZgQ0sDl9G7DdGAJzh4NaTxUaAiuDZK/Lk2iTZWPA5biIorg1gKyPpMY3Z3XRbSrOxqmQ+4SwSfXsx8aLh/PuaN/jz93dEGoLjR87e2OP6SnLJ49vx7n/fZffRu9DqxkPYjJKyiJCMECi2G1t+5fx46muzC7VGNjJ7KAlCR0pgQFWdxadLvlod80mTJrHi86cwqCTjA2bOMzFEMhRMgyUJNhsfLzllQJH9VKcamyHTBmj1MuFzCuKVjVjKr6XckjRrwzcCnvm/9Rx83hCWLlV0LF7Lx++0s3KZw1P/aOXll22efyKLlOD5w+XZJ3M4ef3h7JOqELUWldKlsgp23skkI1xOO6UKFFx50VpWrWxizJRKjjilnt/870oW/CfLDw/M8tZ/+2BZgu2mZbhvaob2HBz+/VX89vp2evczOeTwynC+Bf/HF5bg8+ai/As6Vf+7VKS6TD3VXwI/UUo9JIQ4Au1w2gN4BxiilGrzzQKPoitVlZPRSqljiu4zFXit+EQ/WeAodAnAq4UQA4E+wNsb7Hhw3e4y1fcWDQRKt0GJHyvnVAiOB9PeKFpLbRUx10LM0aBZmfA/m74jS4ROGw2YSXYSV4MDtmsrMzYQk+cE9sjQo+qDT7xK0hK7J/Nb+vHWWQ+yyxF9mHpI4zcKqosWLWLrbUZx6MwRDNxzFJVGZPKwPZN4laqAeQXb2ATH4hIwv6DosyUdMn50QPCb/tY6esgOVjla03rq9kWs+LLA8/9Ys8n9mD9/PlOmTGH2/Abq0kE0gyCnDNY7KaQUoRPTxgjBryu1JIJFNgGqbkUsL173Ob44BgtmOD7941nfti9RoeO1Rmj/QUq41PjRJldf3swHC1wuuGkwFx/zJUef2ZPdD+rB9yqyLHNtHr0ny3PP57jqqlpsT9G7r0RKQXurx87TVnPdnX3YcbLkuWc6+OXPWlm1zOGX19Vy2FHadBCM/Vf+3cHN17Xw0YcO874cUKLxxZ9D8K4VMGnIku+cqVb0G6SGnlTKVBdcvXGm6tdTfQJ4JgjMF0I0A3VKKeWXB2xWSvUo89svgG2UUk1lvntHKTV5Y8f847eiAy12VkqNEUL09NvTperv3aYhBQw9iWNOhECKbXwGHoaK1K4UXgiytgrOidJYiz25nbZBmSHzzKkUOc+K2edkyOT0fTQo5z0L24vOCXauDNsaS/3UBU0UlUY+ZE55z2L3nh9Ree5YHrx4Lv13G9q1B7YJ0tHRwS6Hb8OORw5k630HsLyQfE26f8l4RYg2RdR9FMnQm6AfPgCZnostC1QaEcA0u5UYsYSM9hYXo3wmcpdl3LhxpFIpOrKKmrRuS62E+//czLVXt7F+fY5Xl40np6I4Y9BjJ6fMcEwFNmQ7lvQAybEQ/D6IUQ7s6/HwP0OoRC3glJ8FFdwn5TvwqhIpslFo3c+u7EFOSTyV5bkX6n0i0UxWCe67o52f/bwFgBdeXE+hUOAXV/XgmOMqefNFl7VrFbP+2sx2k3uy696VTN2rGstzwJAU/PcXmEZ2mF7FtrtWsrpJ7xYQl3jiTfzvzYWpwqY5qjqrp4pW7XcBXgR2Bz7xz+8LrPTBdju0aWBN0TV3RBfZbRRCxJG+B3SSZ6/NC5OFEO8CKKXW+nbYLkn3g//DivdlCuqWeash21B+oHkRmJbLe47bVJNeeaFtgypuKxUJ9d72J1nIbj39Oe9Fk8/xdIqtUwRK8UXBlC4tTiZsn+MZrLMrWTioL6LPUu77q81r7T/hX7v+tgtPrXvy4IMPUmh3mHTSeJpsvUhIFS1ieivsyCEDpbttBhIsHAW/cLnpRTbLDmGR98ywj61uBlL62a5alOWlR9dx81Odb+ncVamtraWl2aZXnW7j5587/PwyHTD/6quvwjD97LMq2u672GZsKzNRQCY4BiSccXHTTs6zKDbtBH2NxzVnRCGMLrCJ0pPjqcOW/w+0UaogJDUozrugmTlzCuyxcwV//lsbzz77LNOnT+fjjz9mzJgxTJmaIacky5d67LFLBY891MGO03PsGWTSyRSeF2VWJaIxJNQ0SnKetnH/950cU3ZIYyMTWlYg3zFBjWQTQ6rovJ7qKcBNfvX9HHCq/91hwAwhhIN2KB2lSlXvFDo21SSZSdXi/76c2H68qgIQQjTQjcTbboNqVqVwKd32oThOtdzxeB59cfB5wKrKObbKlQWMFxcJinEAIWjHmWo5UM17Zrgne/GW0aATFwKVGqDgmWQdi9Xt1fQ8ancWXvMgTmuOrZpmItMW8w7oLDyu+1IoFOg5sg7bSIPS7etwIxXXVka4MAT9CZhc0Kd42+OfzXDrF4UlXTrcVGgSSBsOhu1RbeSY+/w6xuzam+O2/e8m9+OTJf1pXu/R0rqKtNELtxn+/vcs//fHdq655hpmzpzJrX8/ihn/20BKuOSUR8HPobWViRdjYQH79IrU/7jdGCJNBXwzD8msP32OoBkdt2QJ10+PjSqCWcIFAzxPq/yG9BKxou3KpCZX4LBjm1m81Oayn/2OG2+8kZNO2pE999wTgC+++IJ0Os3br+UYPawCuyAYOzLN5O3P5n/PuI7pCyvwhEiM+SCcMO6gC8b+xx8XOOeoZdz08BDGTKwgvq8VgOsoPpwXpet+17IpoLqReqpTypx/K3DrRq75Eno31L90o2D1bWibbqMQ4grgCOCKLv5209NUtY2yCFg3AObl6kuW814XV+wvzv0O6qeW7m1VzNDigJn8zlYycTwApWCDveA7j+icvGOSd00M6VG79UBG/PZklt/5L9pOu5UBZ+wHB3Ta9W7LokWLqGysxPbbE2eTYfsQCXYaHItvQugoWVKgIyg3orfF9ossCx+wbXDTgqybYuHHBfqObWRTxXVd/u93bdx9ZzvHHFPJoEEmv76mhRefy/OrG+s4dL8zmTlzJjsfWB9me/WQOQwV1GfwLxRoNn5qbpyJeqGWEg3juCkgbjcvBuP4HmmWcEObdUY4pKWtt70xm/22pKgxs2Rj5qNJk1bRnlW89957TJgwgVNPPZW47LPPPsyePZvTTjuNB+7L0d7ezu23387ee+/NrFmzWDr3GvpPOtcvPpQcy/F44qCvg8aYDB6Z4sMPXQZunUmM/2yryy+P/4RP55VuH/9diNh0pvpNSlYI8WtKy/ntXnyiUupuIcTbaGcYwOFKqS6zizI6/P8vXZFU71qGXXIoA87anyU3Pc68efM2/qMuyhNPPEH/SZsOaF+H9BpWw4f/WsYeL5QWG96YtLe3s+e+1bz0fJ4//bUnMy/RWteqFR7HHFfJDlPTvPnmm4ycUMnQMZvLXvVdlzkv6ySG5x7px4QJEzo9b/Lkybzxxhucf/75nH766SGLPfLII7njjjvItncdeYQQ3P7MSPY4orR+b9PyAp/Oy3Lcxf3L/PK7EeGV/vuO5R7gQ2AomnV+AWyoeLCBzp4v0E2c7DZTXeNUhypX8eoftqZoDx6I7RraiWkg+K54J9Tk9UrfjO1v+RvZvzxcZYb3C9T/QGXekJSLaAjbKj3SsT3k80pQYdlsv08t3uz+zJs3j623LrcxY/dkzpw5fLToAyZtfVgYLgOaTZfbHDBqrwzZalzi7LtE3NLX7yiJKTz6HLItr//9ERr/uxx2614fzjvvPDxPcdff6klnBLafTlrXx2TBxx7tyuS+++5j+33q/ZRVP45UOAT5aob0EmzM8hzavbQu6u2Pg3aVDiM5Asl5VvTuY++7OCoiCOfKuya2klSb2g9RbeRIK5M2Nyoos2V6Be++VOC55/M88WQHubzi/j83svvBG4/flVJy4oknJo7NmDGDc845h0O+t5TfPjCYQcNi7NoPEYxHNPzipM+pqTeZ8auhuIaZaL+rJLXDK/njR9NQSnD3td9cHdxuyXcPosXSoJS6Qwhxbswk8FK5E4UQPwWOBh5BmyPuFULco5S6pis36jaoBqX7Niad7QQa906XO57cbkKVHCtn+4xPrGD3zHigvOMlbZAFzyxrR42uqWK1O5PfBc4Ax5VkCymajEpsI00u99XtWR0dHRx1xveZdsY4ZMYi70Web1O4yWcU25rDhLBwSLBwOD7IBqAaHAskNHGo5MLYXMiEi1mPA77Hf29+lclDZ/LO/l0aTwCMHTuW1+Z4tKUytHlQJWzyzQ533d7KvHnz2LLfljz4cAU3PT6MKpn3VfpkUoYlXO2c819ARtphZbSo1KT25LuxgsgeIgTTuAMzcEQFO+IGz85WklYnE2457lqCtLTC5IlF85o57qhF9B9icdKPLuaxWQcxceJETHPT43d79+7N/fffzxln9GTGgX/gz/8eQmNfHcHS7qX9DLBg/JpsMb4Hj922hLmvtHDT7G1whUk8ZDAe7bJZiAK5Ge1D6EsQoL1cCLEfOqJgYCfnHgtMUUplAYQQV6FjVL8ZUG2NreAhg/SBLS7BsSCUxyUJYmU3/tvIwNCb5rnYXjRZgvCnePZQm5sOJ0kAqnHHlOMZZQGms8+BOEqS82uA5m2Tgg31FVnyq1oZPHjwBtvemXR0dHDyzdvx8t+XseTdtYyc2osx+w+lHRLOKWmoEAg0o/bC7Y49JZFKakD1uxVEC0TRC1I7sjwZfnaVxHYjlqT8fgfn1Ow+mTUPvkhhVXO3+jRlyhQuusim1Q/NypLCUG1U1wiGDh3K66+/Tu8BFvX9K7EVtJLR3v0YvSmu2h+U8DMQENSVDdKFY8MmYVN1o9CzKDTODBcpAEcZFFyTrKOZqiEU7b69ucmu5rXHPwfg0wXNVFR8vaaKK664gttvv53br2rikpv64wrdkTVOdaK62IHnDGHwlJ6889RqWp0MyigFUV1PdTPx/rNZqPvF8kshRC1wAXALOqTqJ52cu4gkNprA5129UbdBtc1Nh4zC3EAwmiVcPCVCz3LAGsqqoUQOqHLihMwjUOEiQJZFbNYqepsdrhWCaDGLCx1VSpa0q4TFeRLXk9iOf42CifI082lfsp6RIztL5Ohcli1bxvhpW9Krn8X2xw7nh7/bEVGRpqAsTDcqNgOBsyqKeIhL8NxCdu2LFCqKvzRcZMyeYLsGbgCwbjKVN2/rYeHYKWRdA6s/SAbhb0juu+8+Lpp5HLvvGwHQGreKxhqPgUNTHHb2SD55p5UJ0xvDAuZ5ZVFQZljqMZD4VuQQxdoWb10ed25WykLCmw7guUkA6nAtqozgHH2trOPH74oK/3k5FDyT5Z93MP5H4752QAVobGxk+fLl7Lvvvjz3h72YdtpTZITNMq8+XBwqZQFbGQzfMcPwHRtR6PEflMMMK8Kp4ijx71D0mr9ZiBCiXim1Tin1hH+omY0btLLAfCHEM+je7AW8KoS4AUApVZrZEJNug2qHG010W0QqdrlQKFO6Jd7YOEAWx4kWi+1pEI4zTIi83YGkpEvG0Oy+3WdnxfZHT8kEO9XHOgF4n62p2L1dT+K4EtcNJqpA2RI372Kvz24SU7388ssZ9r3+TL9gK9qdNB2A7UjynknOtWizNdOTQuHIWNykH+5VzmYd9VeUPKcgaQPAMvyUUcdKfA/RthxCKtKjhpJb8GmX+uM4DkcffTQ/v6UPe+xfwbu5AfpeQte0Pe7qSp768wp22K+BKYdvwXLb3/xQGTQ7uoi4JUtBNPi73LbdxRLfsdQSrjaLxKhs3vMouCZVRiF2neiabT5jTXkGOddi6ZxlLJ2zTBek+wakb9++PPfccwwcOJDxx41Bpk2ancoQNIujAPR7La/RZeTm4fQTgHQ3G4j/SAixGpiNTkmdrZT6eCO/+af/L5Au708FmwCq/78kpWPxWlIN1RibkHq0du1aGqd8/duWfJ1SMXEcTbffzcAbZ7LkvA2blJqamqiurmbnvasQZfZaGjiqklOuGQZAs2vgbCZsZkOy1Zk7suCOt9jrxXN5dteyW89/ZenZsyd1dXUs+zzHwDHVm3ydrLuZxKluRiFVSqneQohR6KyqnYALhRCNaKB8TSl1XZnf3PFV7tltUF2VrylhSgEzDJhk8NnyolzzoAJVFJQvQ4cRgBkLtg9YV6C2F9xIbXc8mXC+BJIxo+1+K007kakVxKCWy4cP7l3wjNBeGsQxurH7uL4NMmCqyhWQN1j62hIqJnU/62jHK/dg3uw3OPL0PWkqVNNc0CyjxU5TcI0ES7Z8lpoyIuYmUTHmmrSfBv0ux8yDZ2C7OjbXkB5evLyiJwhMlablktlyAFXTprDunlmwAVBtb29n//33Z/IhdSjLxFUua50IIAzhkQ63dnFpczOsc7Sv3/ZkaAOPZ7BBcXypLsUY10KkiMYYQHW8RoLfr6xrxUw9RmK332L7erudQqIoSK0l9dx/W1IPz2fde0ujCp7fgAwYMIDVzXl6eilsZYT29A435cdVJ+s8BBI3Xf2/blPdhCLVAl156vtolf0EpdQ7xdf1menHwF+EEMP9889Fq/UloCqE2Ae4EhiCxkhBN0r/dR9UO6pQSiBENKlN6WEKD1NGjhQplH8sKDWnElk+gY0zBFDfORA3wLfYaRwvcqY4SuK4RsLbG0h7wd8yxXBptdPUWPmwLabUoBmAc/EgTPuAHLS/YJvkHDMBqkr5jpzAPqcEeGA3dyCru8cudv3rMbx342ymXrcP+doGWjoyNBe0fbGjYFFwDeIlW60iV6oQfoEUv73B/+X6FpcSG7ESKCXC56mUwHaMkGUKIFVh0+P701l+2W946qmn2HfffUuu67ouU4/eEq+34tiZo3HVMlwk6xyt3rc6GUzhUm1qldv2JB1uilYn7bdLhotAfMEszgrrbFGMg4wUipT/vEzpkpJOAjg9Jck6Fi12tCeaGdugr+AvrIb0/HGXpnaHUSybs7Sk32G7PI+JvzuB1fPXYpjQY9sRWD1rmLv/lZ3+Ji62bfPOO++w/1Z7s9ZRtDgZ1hb0gpN1UziepOAvMjnHjKJTiEBWCrX5gKr61opU74uuSjUS2B643f8/FL+G6k7AjsAgtMNpDtrDXwLAvtyKzqJ6n02wDncbVNd16IkiRbSVsiEVIjbJg+9M6cUGeAS6QMg44w4n02cHwQDPOWbCQ13wnSueF7HIAOxCIBBgGi7thcibaxluEevU7DUYhKadLI4M2nlTsJOPRwiFn/CDcgXSESjPwrO7Hs6yePFiZp/9MGPOmEZ65ECyjktzIcP6rGaqthM584IpkhelrymuXodO8BgIacaZ7FP8OSkfUBVJR47nSp0SExOZsqjadgKH/+wC+ra+zadHXOpfQ3Hgg0cw/7bZkC/ww5u2pdpoZb1XyaJCL77M6oU961iJFFnQC2TeCWIu9TtXMe1Af47aEGesxbbk+GchYjZhocEx3u9AOmwrfGamjN5/S0ca0/AQQuEFURA7bMmSPzxLsQz/+ZEs/8dbdMxfSMXgBtKjBuO1d7D01n9Ss91IhrWkuSTbl+OPPx7Lskp+H0jfY/YmteVg1qpa0rZDU6Ga9Xk9HtrtFI6SYXs78pb/7JOLZ2WmsNnspioA4XW/MZtQpPog4G4/33+OEKJOCNHPv04gr6LB8wbg0SBMaiOyBHhPKbVJRozue//b0yglkFIhpB+cLz1Mw0P6oBkyKeGRMt3wHEu6nTIpQ3rg+jsABPGkjonrRZ9txwhBR4XsquhaSiCkh5TBzALLb4MXY2Qq1gYp9KIQ/EYKhe0YOI6Mrq0EhuWGACRcgZJQ+OwLaveb2unzKhQK9OrVixkzZvBI6iO+uGMOw46aTP30CbTb0GoL2vJpcgV/91dHatOCIHy+8QUjEM8Tmi2jt05WSoTn62P+MwqfTblBrtk2seenVHRfBAjDw6oukBrQg6anXmTtA89R9/PbyC5Zh5MtYKRNphwzku+dNhrXkqx3XZYUevJ+ywBSvlqeSjm0FDK023qhy7smjhuBasCWPS/pwY6PEOX3KTg/3t4NiVKipH+o6FoI/XyFHw3g2RIhFUKC4Y+bda98SMcXqxl0zek0P/Mm2XmfgpSotiyNpxxIv4uPhVQat6BZfo8ftrDiF79n4TGXcQowZswYpk4tP0beffddmp99i0G/Op0Wp4lqo8CqjmrWZDVTbc+lcB2Jk/NBOWugTIXX1EbbG/+hx27ToMqgULEZuUcUyPIBI12qpwphkepJwBvAecAzQojr8YtU+6cNABbHfrbEPxYH1f5E9tTT/aIs7wCvA68rpcqFSv0PMEsI8SIQ2pOUUjeX7VWRdPtNuLYuFuxBOOodqSgIkIYPqj5zFYBhRCqqaUTAL4rUFUN4IWgEYS6Oq9X9YOI4jqFtmkpEjFEJPUlC8AOEjJiJABszATgBmwkmVnCuNGNqtCvxHBlNSKXBi2DrXVegPl2FvXQtX159Z6fPy7IsPM/juuuuo3bSEKqnbofcZxKtBd/W7IdpBW3xHIkKmeOGEEOEiolyJMJLYIffCQGO30dVdDkfa4qPCyKiqoSuF+7WFei550TcVU2kOtYy4KQdqB3RCzfvUNOYobEmR0HmwIUVqo7luR68+8VAanrodM7aihwtuUxojwY/b98PT/NcgfKk/r/TrY1F1M54P4o6XQKyCkTcXOTFxkn8cjLqtJIKPIGb1qBaudf3qFqdY+VtD1E1aRQDrzgJzzMg3QNhVeAUQGQlOAIlFIbRQK8Tf8SKG34Hrsu0adMYfNWVLLrk0pJeTTvxeGqm74ISfZjXpO25q1fUQdb3T/hqtGH779EReCmF12Gz/tEnWf/okwy6+So88yvWaPyapRObalNXKv93o0h1ucGSGAFKqRXAw/4/hBCVwEnoVNWhlC//dwU6WaCOb0P9V3YMaAJVEkASBi+LgMUqEME8EgrD8EJwC1TVAAwD80EcaF1PasblPybP1SxOA5B/nueDYzjhNFKETfRRQgkS5yhFNMF88YxANdZ9UI6MNGFP4MYmvJCK1Y8/TPXWk0ilOi+1KISg+tyTaL/6Fnqcejrp9WlcrxXb8TN8bBPPlTgF/W6VI6EgdduKQaSYuvntNxyBcEjsta7M4KH5P/UEwk0CZnjNOKiqIpwS4BbS5DMWPfY5BFFj015po1IFSMG6dsnH6yxqKzWA9q9p0WadgkFzk7Y153tY5FrTeuwAIuWGzxf89+cJ3X43um+JxN8xZRaKop/p70XSG93ZOpVUdvRzyPmmCHMAjceckDjdRCDaKd3fPvAzNG7B8Guuw8610vTAgzQ/9zyjChm8Qh61tg2zqobWz+aTnfs+fQ7/IbJFsu7jBoQDqXy8MZS8H+kKzHQjfX5wFCsfup/F5/yU+qMOpGHi+E469y2L2jT1H8Ii1Q8B9yilHvYPH492LAE8CPzJ/3sJ2k4ayECK9q/yA/6Dmqo7odnvp8AsylT996W3UqqkKlZXpfs6Q8ASlUDE1SpBtNpLFY2D2OT1fFtV8FnFQE0IhS2SKm7AJlWcXbgyvD/oQZ3Y8ryYrknl0y+iNUeBUCICzABoZNQWKD1H2EYEVlLhZFup3XrLso/p8ccf58iZ54FrkP9yCfVHH4yZS2HkoWNdBrdGO22cvKnBseDf3BEI/1+oQvnAED3LIpunLZA2GLnIQeBUCOJhwMLT38UZhJKUgFeCYfj3NLMCN+WnAXdI8laKjirtWBOuQLYZNJna1t5UX6MX1VYTs90v29dhYLYYBKGhdo1+jqGJ3RNh+4ip5QiFKmeuDoaQ8kFtA/M3GBuJR9aZjQF/+BQtNMpMstv4NUvAOhh/frtN1YM+ux7Isof+yqev/BRl++UEKyqpGDSU3vv+gJSd0byovajdlGk74Fn6NukGvV9cxdDhrLv/cWpKi+F/JyKUQmxCnGp3i1Sjt5s+SwhxP9pB1VxkTwUNoHPQcapXAm/6205vSJ4TQuyulHq+251gU0DVFlqdiq+gSuixGbPFFQ9MrW7KSJs2VBJAQU+iMp79OKMUjkgwlKAtcUaZAFlENHHKqIbxa4SHDRWlPsZZggdBlq5b4eK1ZzGt5JbV2WyWIRedQtsr/8FZ00zPIw6lqmYwRnU1xnows5BZZmH38L3Mts8iA0bp6n/SiQBSeCSVEKHtuUEbhQOyAFYHSH9LBScj8FKcqjFrAAAdGElEQVTRQhGwtcRC4qv3nZK3GHgF781rFigpUFbMS5+PVGx3ZQYjB6nmKADcqTAwcip8D7meAjeTfPnFoK8MUIYgrInij6mSNdMpvxAkPncDVKGMGaX4cxG4QinwxTWBiqpGRhxzPk6hg45li+j44lOMdCVV9QOoGD0aL7k5cJJZF/VJGYABuVWLaXrhaQCc1hYqho/CGrNVaWe+I/mWilQ/iQ6P+hQdUnVi8QWVUptS7u0UdDxrFl2lSuhLfUMhVUaL0ekAD/dK8dXt0oEWA8fgWJyJCJVkh8USTA4lSga2KGY9cfAIBmWRGlV63eBLUcKQlIE2NfiT3O5owW1t5d7DDk6c9+yzz9L0u3sB6HPpDCp7jiC9TkCHBh8Aq1WzP9CAmFgk4gATB9UiYFACwjKiCqQNZk5hFHz1MwduShcI0RcpBgGl30fxe4xZVsLnGlfnygBMcpEDo6D8hUEfrFijUDICSOlI3JRIWJGEpxLvxzMFyozGh5I+0ApKGbiX/FzOtNHZmOo0CmljZpHi0zvRelT0+LHMCnI5m9X/eZ6xY8fywWtPMGb8DTGNJOoD8ccdu7cS4LW3s+SuP+B1aEe229rKkNPOR3Z9x49vVhSbxFQ3oUi1As7s9o02LqX1Fbsh3QbVVItWK5WMAaIo8zlpdgU2PCj1uaUnqPgkJxhwIjlZYqAaTrIyqlrSwRG/b+mxYsbjWf6E9VXBwmfadHPllVfy7LNRuE21H7Nae+g+ZAYNxVopQjCVdswrGmiJTgSiur0aXKSbPFYMXMnnLRCewsgpDJ+pegYYpsBLieg5xkT3WZV8pwyRZLf+BCkBpSIGHzFthXQVwlExFVaDarDmGrYq2x7tGNNfBKAa3k4KPB9kPf8dhOBa1JbO7MXJG3ZyvLM+qqLrxkAu/n3CZi2Sv1ECaodsRaq2gQ8++IB0fW+s9th94u6KMguF8jxaF7zHytefInA0DDryNKq2GIF0DOh6iYZvXDYFVDcX+aq7qXYbVM129OA2SYKoEV+dRZcG7VdJZUvYsuIsLsaugvNK2CwkJ0y5Y0WTxPMLh/m7euCt1vvDz5gxI/zNM888w/5H/5C6A/ahdo/pyKzEyGvWCL7N00uyc+kopBuBrfCUz1IjIIuDU9C2hErot1MWFNLP/TQIgMgHIFNsdFEDfJU7plEUsUh97+SESdgWA1D1VOQkC1hxJyL862lAj4GqIRKTM1jQgj65KX1OZ+p6l2PhOztvA6BaApplrhmfBwHLFkow7IDT6Vi9hNrBYzAKsXkTe6fx8Sxd8JTLwqfuwM1luficM7jq2uto2HkvarYYjbQ7Z+LfiRTPtf/HJL6bKnA12rTwe+Cb2U3VzGkWVFY1Cz4HfxfbtuLMIVBzi9V0in5XbtAqkhkbRYM/boMUlAfVYoAtNyjjE0b4rEi1g1vIseaxR6mdtAOHHHJIeM6xxx6L0bOe2h2nIfMSo0NgdGg7KoCZ1wZ8JSM2KDwNrIHaLlzlM9ekOhz0TTciyVyD5yptD+F64XWVBGn4HnbTt4WGwCDCa8Wv7xkCGThm/OsmWEcn6nRoHgiYracSi0BiV/rgeOIdqASoCieIB46/NAFK4flOMyMtca0iU00MvIu1nEQfio5vqOpkMUAEWoKSnZMHfe/I0aa8iIVXVDZQMaRBA7VXtKgHuRfx5+xBbv1KWr/8kKr+w7nm+t8yYI8jqBs9CdmeJBGbgwg2P6YqhMgAP6Z0O5WTypz+7e6mKgsgJChHM1Z9MKmOeoavmsUHapGKHQyWEja2MSlipmXDaoocCOHnYpOBigHBRkA1YESm0jZBgHSqJnH+ZZddxrnnnou3ZC2p/gMx8vgsVF/c7FAYBT9v3YwuLlyF9EG1mOWVE+ElQSsAG61ye+E1g3br/yXKjByFoS4eY55CgTQEnplEGKE23J7Sc1WC0ernX/SA488f/b3wIlBVXrC4hAHJKKA9v4ZMbW+EkHr2eiLqC6A0gsc+Q0lISUmji4C4pE9FTQ/UdJEkFvrk+HnRIqbnhyoxCSCi5UbENKNi8lFV25/xx1xB26qFpOt6U1HfF9oDYrJ5ARhKJRfDzUP+it5OZW/gF8AxwIJOzv12d1M1bL36ejGPdTFT1bGj5Vf//6+9c4+Rq7rv+OfcOzM7+16vvesXGHv9CNQ2Nq9ggw0hNTKUNKEkaUKlVk1RSSI1/YNITQIkJEE0TaVW7R8kNAJKVQWUpxwRIkgMwkAAQYxlAtjGNrYxdliv9707uzNz7/n1j/s6986d8e6ymK0yX2m1M3fOPefcc8/5nt/rnGPaPquSo+EoSbs/5pxII9UYgUpUTgWpxsupiOE06iKu73zTivLwMNmWDpZeFK2Df/XVV/nyHV+n+4abaJy/FFXyHFNW2Xxmj/issmCHg01hORpV1lGaaq8vNgkYpBpU1RWCrZ/CgRYE/9sacS3ECkRkiFZBRG0llkJNZfuoBBHFHX8SJ9IkgQbPYpgHwoR+w3gRGUZojgivvvVTegdfZ83yG1i2eBPK1ViWQmxjj4akmUOpONGm8Y+SGDGfCZHqL+EkEJkEIhu2UhLafcWPjfZCDwnvwV9oEMvDrJomtKk35FpoWLoeywUpzjEiTWCmcarvI1aJyKeVUp8Qkf9RSj0MPGEmUEplRMThbJ+mqlw8VZPoqIukOqV8T3lsoCXGqZWUHpPp0yAmqZpSkFF2GlEmVVapdPycCToDgUMon22jPD5CeaA//P3666+n8+ItzN+4BcuPgisef4ejT/0U286xfMtnyNjzPNW6LFjFyH4RkK33OaUySWkkINREUk9KjCQ7MF6Lq8C1UcESHctClKokO6WiFUg1pDoxSEgl7aXhRFap8qc+U/A8ShlmHYnfJ0JP9xa62lbT1fEhrJLjLR6wiJGqbjBOMlCBXST4nILQ/mpI/QlUEB1R0pA0LTxiNuKw43n5k5dJ3oF5JVhEY/n3mEn8vmoKEeZClqqhYh8kZO6p/0THqQwppdYB7wLLE2leAi5OnKaqmOZpqtMnVfGIVRNtbOH6NkJzmbny/8JYaIeQFIN8klJmmqc5ZtfUvlfcINUKx5RBnGY5MXU0KDdp860B8W1dblaRJcP8nkvYt+PfueCqXuY5r3Py5Ek6rE4a10Tq/qlXXqDQ9zYA9nARq8UbWGIRxnnaEy6q5KJcn00qpDc8okyV7OLPyhmkA2XudqWURzTBQA4up96oEpKpilbK+eQcoppt0jSNmmRq/nfjHSIpjbflFtDW0OXlpV1UymyoyhlDElSRsyicDILrKc8InrQYmEz8axVtEhIxFdpazElmmge0qrCzVoazeWwpxn0BglcXMx1UqddcwBxU/3+glJoHfB1v0UAL8I1Emsi8LfI68PpMCpo+qeogHMYIWNceaYY21aB6YiyySpBqpSRD3P4ZlIcx+fsOnLj6n3DaBOmTpGo4flLVfON6DIZn2qzbkhVX0H/4ZfY/662Yu+eee7jjjjs4/Oj3Wbb2elo7l/GhtTexZv1N4T2WK+GkEZvJLcJ5VAUEapQdSnXVyNNMX8PL7s2EKuWa2XgqXUI1oMBYpKFRSWJNk1wNkhDtS8jJdOb+DK5EEmyQR/A/xckW3eca6rTyzAkGiaqwY6Y8FIYpQ0XXICJa74vRDkSdSbkSplPEJUpEUHF3XbzP+xfC8Cl/ZIYagdmfg2cL7ppDZAq+dD3HJFURCZa27gJ6qiTrUkpVPSolscqrKma0tY1YXoB+sJue0t5KqODdWn6nNAUJqyShlBsgGVAdBr6bpJEk4sCBYXZWg1TD68ZLrQwL8okyMahrGfwlLNu7p73lHC7/82/zzptPc+roS9z1rXsAGO49yET3ejqazzGIOC5VJwPdY2q8CLhSm0SDdObngGxM51XqPYZUGoQ6BeVbiXzNfJL22yBLfNun+X0KSDU9JIWbWBvEfog+JshROTqSSm1vKbVYVrShjnjSa9inQilWovx8CTJsQkW0J4EJQ8r32tKomsJz6Bre/1AyNfKK2aODfqpUaMXRWUIBJXxGHW+DMNJmLvFYmhnrDKixSfWPgGAn+A5gSEQ2+jtZ7QMO+L+9KCJfqJL3QrzwqCUicr1S6k+AzYld/m08CfY9TVPT9/6XvaWQokKfAm5OkTEM58HKGXMWtkoSSmohEmpYakwkRKQikRSqYp3MVJUDqc6otMS90d49lSqWibQwGhP2hEOmMcfKnu30rLmO0uQoWIpG8Y/wLuqKTh7O4FqwAsdUWaPKbmRLdQVVdsBxwfGDV7WGYgma/TOIMpl0ojVV+YBkk+SqjRFsqXjvSSGOKFogeIi4auy1jVXxeywP24pHK1Cp2keRAoZNuJbJw6hzfHKMXpzS3vMrpaM0tkIFZoHguWJ7KnrspAwzQNwtH9QnIN6EjTT23EY5VYapueAhyBeEwJ+oXSs+RrQkqxIL0ZsTEKbm7KxE6ibVIvKZIIFS6t/wDu8LcFhENk4h74eA/wbu8L+/CfwIb6+BAH8QkW/PpOImZuCo8ju/SlwzCDLsIxKt/7b8FTaBQyZIkOwMpmRX6XRKqMMQkkSq99n4nuaBVjGJobbUqhKDV2mNcr2YSbEUGZrABctxYx7rCvj1UGWvcOVocI17XB0R6oS3akAcB1zX+wNUW0u6JKoT5FbLFBDWp4Y67UuTFap9oFYH35NZqmRH8AnIbIQkqU5jDIa2zrCfGFJbbO8I8dQmy4okaUN1j9JH5gAvAsKPkQ3tpQkiC8unos0qhITACZUmmQf3BDZcSxnmAu+f5SSWfPvmI9O2akrEcwKmw3Rat1XdpPoNINhw5S/xNlWZLhaIyI+VUl/z83eUqjifYFYMKTMLqUrYm4DIe0pkHjDJTTnR8sUYkp0h5hWmUu0M8k0JSI/SVA7YtBAkkmSRQAWxEq+Pcvxg+2CQORpVMtYKKm/FT+QEicgrIFUc7ZXjBN9dKBSQUhk95m1bJE4ZqyEfSUWuhuT+mcryvc9GWTW890AktU5lQBr5ms+UnIxCx1DFJKejdxZ8TtpLzc9pEQ9++VUddSTuUSo0T4SX3ICQfEIOwwDNd6RjQoJnYzUl87jGYyLRGsbOZzXGayAYmJOQT7Taj7GJmQjMiU55Y09NwRZ+NlFFUp3pJtUBtgK9InLQuLbCD9IfAe4UkWerVGncjzcVP/9NxCVegD+tcu+0MH1SLbhx0oPKl2lVDlxv+aUYXvsqHSAZapOErxJWTVPre/I9G2RS0eXT1M5YXlWuA9jhSKqU6gKbqWNMkhOTUPK2hNOjY+jJIiROctDFSeycv0bWcby8MonXV0s6rdZeSkX256SEGUixFgT6aEiopu3WbEeto/djLkTQOixHiXgTQ7JOU5Gsq9qKfWgdS6Mcv75VQsBiZgD85wwk1eRKrZSJpBbBp9VXEu2VBs+cosLqQJy8xTBZxBY3zBVSDd5vJWa0SbXx083AI8b3PwDLRKRfKXUJsEMptTZxT4Db8Lz+K5VSvwW6gE/Fqy0DZ6rbVDCDtf/+pphpcYZBOyZJJEijidvMzP/JPE1UI7Pkb9OFEd94RmfQdH5zXESE08W3eXt4D32jB/nud7/Lb762G3teO0pZSBBCJYKUy6H9VJdKqVlaDXlUsBl2NltJqDOtb1L9T5KsnRiwWsW99GU3XdvQOlIBNd4kEuahfXOGyRQSr0tFZEB1c0PV/hFI7iY7JqIflGXF87R0SLoVsa3GuokzTgBVtARVhZzjdXKrtoOoYJtNUn+fM5iB+g9QZZNq/G3/bsLYrUpEivjHnYjIbqXUYWAN8DsS8O20V+M5vBRwQETKyXSzgenbVMtudXXaHHhQqYL6UhqkSJpnkmBr/XYG6SD1mlKG+utLMkline7AEaGsS7w19CJHBl8KL1vY/Pgrv2KetRB3MKFxBBKpOoP4YlvRJGAp73NQP629umuJD7YK73pK+1kp78iKypG0wZ2Y/EQliAkdNy6aeQf1MPpPySkwUR7CVlmabW/SqaynFUmLZlm1+ks4OaQMcIOkR5x+ynqSzsZzUdooR4SSW8CyM2Tt/LTsvn2jh7FVls4mY2N6S8WcehUmgXByi9tPYpJxIKXqxD0wtyRVZ/pbZtXYpBq8QPz9IvKOkb4LGPB3lerBO1X1rUSelwHHReRd3456CfBJ4JhS6puzJZ2amD6pTjpTU9+moobWkkBj9yQNV9Ulj1TyTJadJJypDJa0CSRUj736HBvczf6+pxA/w+985zs8ebuxW1itwxlr/GblcijzNM5kW2nx7ldW3LufVkbFcyQkOZOszf+BKg+IbSOZGpOAZXkqajZco4m4mXCBg0xMMlo+zeD4cd4Z+z2TzjiNditlPYlCcX7zZtozXeTtFopMYmGRtRp8ydmuLsWJZqj4Lo4uUtaTTDgjvF34PavbNmGpDJay0KJpsJv5+oOfY3h4mEcffZQ9u3ZRGNBoXBY1fwgsmCiP0Dd+OIxD7WxaRktmPm0N3ejSBO8WDuJImazVgKPLZO08JbdAURfQ4jLhehpoT+uliAga11fbLXLZJs5pv5B8rj2xcEIq+6LlzwlJyTStX88ZUqWa+n8mpG5SLSK/Aj5LXPUHbxepbyulHLzYlS+kkOR/4REySqmrgH8BvgRsBH5AwgQwG5iBpOpU2kzSBuuZ7DxaopVCZvq0PMzvwc5Atci0Wl7ViD5NtXTTCCgxEVjemrKxYh/5TDv7+nZy8cUXc+TIEXbs2MGVV17Jk7d/Nr3MacDbXMSwFQYSa2hKSfyvhXBQGlJvIJn6hBqGEgV21GTIVMKmKna0/r5QHuT0wAE6F5xPPteOZdlMjvTRP/YW4xOnGBo6yvj4KRry7dx884187GP/yI033ohlWYgIDzzwAI888givvPIY5QEoSgERjaVsslaeZa0bWb7gw1H9RCgVRzk89CKnCm9hi0WD3cy2v7iSlpZlnHvu9Rw4cACtNeVymWw2y4kTJ7jvvvtoamrikksu4YEHHqC5uZmdO3eya9cuMpkM55xzDtu2bWP58uWcOHGCZ555hr1797Jv3z4ymVZu+/g3Wbx4MUNDQ7S3t3Pq1CkWLVrE0qVLyeVyLFq0iD179vDkk0+Sy+VobGzEtm201nzjK//ModPPsaRtHS35Lhqz7eTsJiadMbqaV5CzG6P3pf02Dg6pFO29n1gctlT2zQ8UMiP1v9Ym1SLytynXfoZnKqgF2yDaz+A5x34G/Mwg7lmFkmm8CKWUXLfiyxGJBsQa2M9i9lPx7IamHU0ESdpUg7xtG9XQ4Ku5wRbxxmAHz46YNlPXUtNNj/h0Ol0aqULc614s82bvUxwtvMp5K5bR3d3Nfffdx5133snu3bt5/vnnuXXlP029zCqwcjms1lZo9Hcsy2RApwSVWnalk8yYMCqWkxrqdMkpMDL5LoNjx8lm8yyat46+4TeZLI8itoXjTCCiyVp5Whu7act147olRLsUdYGRch+j439gcPQY7fNXMtj/JqJdmlsXM1Ho569u/hQbNmxg8+bNXHbZZdi2XdMjXigUOHjwID09PbS0tHDy5EluueUWnn76aT7S8jeMOKc5xkHKToGRyV4+9/d/za233srGjRtre9rnAESE1157jRdeeIH9+/dz9OhRent7aW1tZeevn6Ih00xjrpMl7etY0NKDnWvAUt4a2PDJTDOaP94Q4fGD/4okz38/y2jPdssVnZ+suP74qft2T8VRNZtQSr0GbPRV//3ArSLyTPCbiMz6GTTTX1E1OoaUSl7sZADTNpm85ktPEqioJpSFyvpVyGW9P9uOS462bbhAPVtiqgc14a1PcwIoQIvLieHXyGWayVl5bCtDa647NeYwSdySsUPHTe/YId7p383pwiEee+wxbrvtNg4fPsx5553HL3/5S+69916uvfZalst7H+SqoQEachFh2r532DWI1bYhm0EClVspxpwBjvT+FkvZDE+8i9ZlOlqXUSqN4eoyuUwzE6UhmlodxsfHWbduHddt30h/fz8/+cl/sH37dtauvYzBwUHWrl1LNpulv7+fxx9/nBO9u2lvb8eyLBZ2d/ORCzZw+eW3cumll7Jy5UpEhL6+Pl5++WW2bt1KW9v0DqVrampiw4YN4fdSqcQTTzzB1VdfzZ/dtI677rqL22+/nfPPP59Vq1ZxwQUXvKc2PptQSrF+/XrWr688/XRycpLjx4/z3HPP8fOf/5zHf/V9tDhocVFY5DOttOUXcW7nxcxvWeHd5LhQLs+dBVUikSP2g8cjwC6l1GlgAngWQCm1isqQqlnBtCXV7etu92LQHNdb5QMwMYmemECK3rkhqQ2qLE8a9UlU5XKopsZI+kqSKXhqjuGdFaUo6nGO9+9mYedaMnae0dETNGbaUBpaG7sBcHUZ1dgUlW0BGsZGT/C7g/9L2Y0fpriqYzONViuOLpJVDZwo7KMjtxjbztKeX0LObsTONDBQOkk+304+28bRgZd553RkL21oaEBrzTPPPMOmTZu8fFetov3wEjrVwim3cazJfInd7uiAfEMkJWcykLGjbfz85wMo6GFODOxleOIkw4WTLF+4mb/7h2sZGBhg69atvP322yxZsgSAI0eOcOmll9LT08PSpUuxrGim0lrHvn/QGB8fZ8mSJYyOjnLFFVfw0EMPsWrVqg+6WmcVxWKR5557jm3btrGs63LOP3d7+JtyPevvr/fe/cFLqpku2dz6iYrrTww9cNYlVSCISV0M/FpExv1ra4AWEXlltsubtqQqDVkk5wff5z3niWpuxCq3Rh6/srcCyJRmVSBxBqFA2SxkM7g2OG6RTNZmrNyPQ5miM87YRC9jk30oZeGKQzbTiCBMloYZGj7KWyeeprV5MaPjwYm0iqbcPJRSjBf7WTh/HblsK20tS2hq6UJsmwKDIaF+9KMfZcWKFTz//PN0dMDy5fNpaWlhYGCAm9bfguM4DA0NsXfvXvr7+zl58iQXXnghtl1g376nOD10mp6eHtavX8/ChQtZuXIlW7ZsYdOmTRw6dIjvfe97DA4O0sRMDnMMHsm3aZZKIDok2aAdlWnWKJWYGO3j2QnPxLRjxw6uueaaaUuIAeYSoQI0NzczPPy+CBZzHtvP/RLH3P2cHPw9E+Uhli34MGvOu86Lugj2FrBr53FWEYQJzhGIyIsp1958v8qbPqma0cjB8roMkM+lxpiKeKs9dKCyB1InwtDESV5+w9s8Rimbxqb5ZDKN5HLNNDcvpGvBRu68/ePk83kGBwexbZuWlhaGhoYoFotYlsXw8DAiwtatW2lvb8d1XRYtWsSDDz6I67q88cYbvP76CyhXce55C2lovIg9e/ZQLBa5//77K+o7FRw8eJDR0VG01ti2TW9vLzt37uTuu+/m0KFD9Pf38/nPf55f/OIXfOuq/5xRGQDieB3THY06aEHGGGeEIhM43hREgTFGGcLFoZtz+MmuH3LVVVfNuNw65g5GRkb47bs/JN/Qwfld19AxbyUql60MdbOZO46q6sH/fxSYwdZ/2juWIxk0DR5xEo9hfen1+xkeO05zYxeOM0FbtgtHlxgtnqYh08y2bdt4+OGHmTdvHpm0gPYZ4qtf/eqs5WXi2LFjrFmzJvy+YcMGmpqauOGGG/jiF7/I6tWrWbNmDbYvVf5Gb6nI41rr0wBMSoFBTiNo2phHi2qPpRMRBuljnBHGGWWI05Qo0koHeRr57JdvorOzk5UrV3LRRRexevXqOe+kqWN6OHDgAJPuKOvn3UBLttMLC1MqDGkL9pFQJWfOkKpQxQT4R4Jp21Tfx7rUUUcd7x3HRGT5B1kBpdTjwIKUn06LyHVnuz5nG9Mi1TrqqKOOOmpjbnkj6qijjjr+n6NOqnXUUUcds4g6qdZRRx11zCLqpFpHHXXUMYuok2odddRRxyyiTqp11FFHHbOIOqnWUUcddcwi6qRaRx111DGLqJNqHXXUUccs4v8ADfUnyk4ng+gAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 864x144 with 2 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.figure(figsize=(12,2));\n",
"ax = plt.axes(projection=ccrs.PlateCarree());\n",
"ds['thetao'][0, 0, ::].plot.pcolormesh(ax=ax, x='lon', y='lat');\n",
"ax.coastlines();"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Check grid structure"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.DataArray 'lon' (j: 149, i: 182)>\n",
"array([[78.00004, 80. , 81.99996, ..., 76.00008, 78.00004, 80. ],\n",
" [78.00004, 80. , 81.99996, ..., 76.00008, 78.00004, 80. ],\n",
" [78.00004, 80. , 81.99996, ..., 76.00008, 78.00004, 80. ],\n",
" ...,\n",
" [79.77896, 80. , 80.22104, ..., 79.66438, 79.77896, 80. ],\n",
" [80. , 80. , 80. , ..., 80. , 80. , 80. ],\n",
" [80.22104, 80. , 79.77896, ..., 80.33562, 80.22104, 80. ]],\n",
" dtype=float32)\n",
"Coordinates:\n",
" lon (j, i) float32 78.00004 80.0 81.99996 83.99992 85.99989 ...\n",
" lat (j, i) float32 ...\n",
" * i (i) int32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 ...\n",
" * j (j) int32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 ...\n",
"Attributes:\n",
" standard_name: longitude\n",
" long_name: longitude coordinate\n",
" units: degrees_east\n",
" _CoordinateAxisType: Lon\n",
" bounds: lon_bnds"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ds['lon']"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.PathCollection at 0xd1c30ae10>"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsvXt0VFW6L/qbq6ryIpAHCeT9BkISBJJAi0AC2kceQaBH4xZUsK9gC7a979njDsXTArYS9rnqOafH2K1gK3ZvQRtb7NuABNBuJTxalBBAIQmvvJ8QTEIgz6q15v1jrjkz16pVlaBg78b6HI4Kq16r5vrm7/t9v++bcxFKKXzmM5/5zGd3vin/6BPwmc985jOffT/mA3yf+cxnPvuBmA/wfeYzn/nsB2I+wPeZz3zmsx+I+QDfZz7zmc9+IOYDfJ/5zGc++4GYD/B95jOf+ewHYj7A95nPfOazH4j5AN9nPvOZz34gZv9Hn4BsERERNCkp6R99Gj7zmc989k9lpaWlVymlkYO97r8U4CclJeHEiRP/6NPwmc985rN/KiOE1A7ldT5Jx2c+85nPfiDmA3yf+cxnPvuBmA/wfeYzn/nsB2I+wPeZz3zmsx+I+QDfZz7zmc9+IOYDfJ/5zGc++4GYD/B95jOf+ewHYj7AvwNN0zR8du4KNE37R5+Kz+5A8/nXP6/94AH/+3Le73OSHDx3BT/fVoKD567c9u8C7swxvB32z37+3D47dwU/33YCn31P/nU77U65JkO1Ow7wKaW4dOUGhnpz9uILV7Hm3VIUX7h6yz7z234PcGscUNM0UIrv7MRDPZeh/jZvNpQxvpnvuRXX7FZ/x2Dn/13PeajX67v6GNU0UEpBv6OPDnYOt2IuDDamN+u734df3U674wC/srULT//xJCpbuwZ9LaUUcaEB2PJINmaNjRj8M6/cGNLFtnKKWWMjsOXRHK/fA7g74D+SPQ91Mph/281MCv7ayis3Br1u+WNGYsOCDOSPGTno55r94FZMVPNn3Gpf+65+NtTr9X372Lf1LavXDPU6DtWvZo2NwJZHshEXGjAk3/g+/Oq2GqX0v8z/OTk59Luapmn04uXrVNO0QY9fvHydzvnNIXrx8nWvr+PHLrR0ur3e6jsuDvF15nOklFJVVemnFZepqqqUUko/rbhMxz6/TxyTn/Nkn5xppCnP7aWfnGn0+jr58+Tv8XQuQ/0NYlxbOj2+xvzaCxavNX+H1fXydC4e3zuEc/L0W83XdSjf6zYmQ3jvt/Uz+Xp5+3yXy/WdfOxm/etvZc0efcvlchnOUz5vq3MZ6nX05FffBQesjn9bv5I/x5Nf34wBOEGHgLH/cJCX/78VgM9tKGDh1aEsBl/TNHqxpZNeaOn0CnSDOdlgk9+T0/OJ+deyZq+TcigTUlVV+ptPztMxvyqynOjyOQw2STyd/1CA60JLJ704yHgaxsli/L2di/xdhnMa4gTVNI1+WnHZY1Ayj5Un/zGP77f1M2/nM9h5eBqnm/WxwfzLDPTmz/J2nkOdGxear9HZrx6kF1s6LcfPk199WxwYql99WnH5pnzKm8/ejP0gAd+bI1lNIm+RXVVVa6bW0mnpaGYn83Yug01+T4FDnkhjflVEf/PJ+ZuekN4+w9M5W42l/Fs9Acl3GUtP7/u252L+7KFMUD4x7/8/xZavs2Lans7dfF6efMX8O83fe/Hydbfz8Xgekr8PdZyG4mOD+RcnEjLQe/It81zwODdMc+FCSyed9epn9IKF78hg6smn5OszmE95O4+h+srNvu5m7QcJ+FZOPNiFHQxYzBfFytGsnMztXKTJ523yy+fFJ7AV8P/mk/M07X/svakJOdhktAIuAUhDPH9vzHKwsbQ6D/m4FTiaJ6Ll9fYgFw028S62dNJZr3xGPy1vsQwaVtfH02+QgcxTgDeDOD8/b2zUivHz75DPYbBrZh4Hbz42mH/J7xksGJnHydP1NAdoT1mQ+bd4uh7ego6V3wyFKGiaRj8tb6GzXvnMY+bhyee+K8v/QQK+FTB4BTEvwDDUCUeptZN5ZAgmB/PGSOXfcCsm5N/Kmmnq/9hL/8/HFV4Bw/z93ljZUILoUMeSf5Y34PJ0nS9YXGdPYCH7i6cJqmka/Vt5i9fg7g3o+fX3FCTNbFOMj/SZ3sBdBi7uG+bxks/BKhDK42z1mzz5mJV/eQL7wYKReQzlczYDv9mHzBmiVebk6fs9EUFv5+DmyxZE4kLzNTr9f35KLzRfc/MnK7D3FJBu1n6QgE/p0NiXGeTdLvBlLzLQTTjZUDINK/C0Al4rZ1FVlf7vj8/RlOf20r+WNYsx8DQhza+9GSf0xMo8ZUlm1moZKE1j6TUAWIHnINdWDsCeJinPMv5mYvFW7N5qvLwB/WDjNJgPmn+H+Zw8ZVdW5yADHj8uPycfk3+fld98/FU9TVq7l378Vb0YLzOZMI+VG9A3X3Mr2g42Hub5aAZXsz95C4aDXRP5PZbZxWXrOsK38Sdv9aGh2lAB/45py6TU2IZFAbz2cDbSRgUjNXIYKlu7Blq09Bao15ZNBgXwyx2nAEKQNiqYHdMHhxACAHhyeykqr9xg3wOAgrJHSnHwfKt4PyEEla1deHl/BZ6dm46UiCAcPN+Kp/94ElVXu8V5JI8MFM9XXe3GKwfOobK1S5wTAfD0H0/ikt6eBwCzx0Vi7bzxeHl/BQ6eb2UXT1EwMS4EhABnGjq8ttUdPHcFm4sv4enZqbg3fZQ495f3V2DtvPGYPS4ShBBxXB7DlIgg0Yb2yx2nQCnFs3PTkTwyEJWtXUiNHIaqq92GlsKUiCA8OzcdL++vQGVrl/idxReuMqZhMZbFF67ilQPnUHW1W1yD15ZNRtqoYIAQ/HLHKVS2donPppRattylRg4b+O4rN1DZ2oVZYyPw+iM5oPrY8jFMGxWM9QUZKNxbLq4x9yenShEfFij8oPLKDWz8qAxr56Zj9jh2NznzWKVGDkNdew8bJ+kYtyqp9XIwH0yJCMKzc8Zh40dl7PUAQID48CDx3fx9/JxfWzbZcA6afq2gj1Xxhavi3Chg8Le69h6kRA4TfvbZuSuoau3CXbEjQADsOd0Il8uFM03XQQGcaboOAHC5XNh1qhGaBkyIDUFlaxcOnrsifGvW2AgUX7g6ME7LJqOuvQevHDiHQxe/8Tgede09SI0cJsaDm5iPhMBhJ4bf3+/SDAD3+sPZmD0uUvgn/35KqcCDKsl/ZN/hx+VjB8+3AoDwSQqK2rZucX5W/kQpew2/dvI841iRNioYaaOCha/dLrtjAJ/3x3oDeg5UlFLLCQZAAAsHBBmUAHbx1y/IFJ/NL1hq5DBLJ3t5fwWemTPO4GCHLn7jEeT5pHvt4WwQAL94r1Q4mRXoz04fhadmpeG1g5UeVz5qmoavGq6BUmBCXCgIIV7BXv5NwABI8fED4HWy8gA3e1wkXn84ewC8pADgbSx5oPQEZjx4ABCvl3vYrQC+6mq3mFQyiBJCGIASGCau2cyTFoDlhOXXnB8b1Ad1++3SSW6gVtXaJc6NmsdLDz7PzhnnNl6y39XrwEr1scpLCzcEADPIHzzfivjQAKydm46Ne8vwsz+UID48CPOyorDn6xb8x2eVyIoZDgIgK2Y4NE3D/7PzDD4604IHJsYgITwI/9cfSrBxbznWzk1HXGgAivk4zRknft8rB87h2TnjEBvij2f1+SGPx7NzxjHQvdptMSfZtUiJCBLjweeuw07EtTH7z2+XThKfLxMXQRy8+I7V3LMCd0II4sMC4VSpuJ6VV26gcG851hdkIG1UsEey9X3YHQP4qTpIeptkAHM0EGLJXDlwcWerbO1C2qhgvLU8FwBzmqrWLhTuLUeVzjQ5gwEGnIwCgom+/nC2+HzOtvLHjBSALrOttXPTRSAA4DYRAXfQJ4RgQsxwdn4eGL6Z3VsxVZ4dvby/AmvnpiM+LNAj0PPxyx8zcnDw0idr8YWryB8zEmvnjUdKBANMT2Mpg6YBzKRxfe3hbIAQlhFIgZNP0uILV5EaOcwyazMDvCeWL5t50prHkAdMHvBnjY3w7oOAYezqO3oHBTV5vGS2X9XaZQB/Phb8u2Tg/9OJRhEAVkxLEv4WHxaItXPTUVhUjp9vLwUFsL4gAw4bQWNHLx6YGAVF8jGFAIqi4OC5Kyj6ugkL74rC/1qShfqOXthtwLqC8aAAnny3VIA/Z+6UUqydNx7xYYFY/d5JbNhT5hYIZSYuZzoAsH5BJiMcF64OjIfub3yuWgVDPsberoUn37GaewCGTBbkbNGT73wfi6/uGMAnOlsHTGm2lD5yoDIz19q2bkMqGR8eJJgpN55GyoxfnoCcoT4zZxzq2roFM+CgJE+6qqvdAGBg8r/ccQpxYYGCff3iPbbCcNbYCEvQf1ZnYJeu3ICiKCAEONvUaSnraPpWC1kxI0AIsWSqfLxefzh7YGIOMn7VJqZtBV4ys+eZDZd1xFjKco4+ec1Z0rNzxqFWGld+DgZJRDcesPl4mRli2qhgARqVrV2GIOBpuskASwhx/7cs7wEeJQw+VoAUPJdNRlxogDjvuvYerC/IEFmUccyoge2nRASJ6xkXFigkm+M17YxFS9+1/O5E/Offq7D87kQQAO98Xo3ldycKtq9RinXzx2NdwXhxbusKxmNjUTkSRwbj6XvTsPlQFb5uvGb0LQALJkThSCX7zg0LMoXfr5s/Hs8XjAelVLB6eSxWTk+GXSGo7+jFKwfOoa69B8/OGYdXDpzD8Zp2QRbiQgPgVNnV4XMiLjRAjEfllRt4cnupuN4gQLw+Htx/4kIDWCDUz818Lfh3e/IdPvf4fOSEUPYlSzNLTx58Z6irtr+L3TGADwxo6pbpo4VEUdvWLZxLlne4mRmXOa3mgKVJzC4xPEh8H2db8nE5IMhM/rWHs6EQIib8Wh0kPYF+YngQQIH6tm7MGhfpUdbRNE1orURRPDLVZ+eME9ovn5iDjZ8l8zeBV9XVbnH+PCPgk4WPJQDB6CmlAxNXv6avP5yNeD6ukjxm0HlNgYqPnzlrkycuD+qUUsN1ZYNlnKSG6y6xSXOKzgOTJx8UYKOPFbd/ff+0IcvjpCMuNADP3D8WGz9iLNjM9g/pgXJdQQbq23vw8v4KLL87Ee98Xi2uzcv7K7BiWhIeyo3Fz6anYNuxGtS19+C5eeOx7VgNatu6sWJaIjbtq0BhUQVACJ6dMw6FReUMnCjQ2NGLzKhh0CjgcjqhUUBVVSEVNl/rFaw1PiwQhXvL2e8nRHzu8mmJAviX350IANj+RS1WzUxGvO4zPPtYOzddD0gJ2PhRGerbe8T14L8dgCHocfLArw30+cSzi399/7QYbwPwE4Jn9XEfiu9wf06JYIAt/9vsO1a+4tF3pDl2u+w7Az4hZBwh5LT0fych5L8TQn5NCGmUjs+/FSfsycy6WHx4kJj8rz+c7SZRLL870TCxzGyeg9nL+ytYemuaaDz1f2t5Loj+HmBAhskfM9KtQMaBn09KzqzMEo7Z8axAPyVyGNYtyMDGonJUXe32KOvIcs7scZFu7F78W2JDhvQbwOuP5BgLbx6AnhsHL54OV7V2CbYqyzp8LHmKP2tshCiEmRkz+HXk11W/hlbFNXFddR3WHAQEe8NA5mYuKpsBXi4oG9gkBrR8Pm6vHDjncQyhj4tcWJWJB/dH7tOr3z3JvoRAsGCZ7ceGBujySIAA2SlJYeIzXt5fgfkTYrD9i1ocvtSGdz6vxrysaLz68XnEhQVi+bREvLCnDG8drsLK6cl4fn66AHqqUbRc68Gv5o3D+t1nUXyhFRTAibproAD2fNWMzcVVoABab/SjX9WgahrqdJ9qvtYrWP7KmcnYeqQKG/aUYV5WlAg6a+emY+uRKqx+75Qh4GmU4rl54zElMUxkXm8+miPGhYKivr3HEPTWF2SIMTX7DA8yPAuVC8LAAIjzuccx4cntpUzukuYe92d+HeRMn/uhXGsx+4pcW/i+tfzvDPiU0vOU0kmU0kkAcgB0A/iL/vRv+HOU0n3f9bu8GZdUZE3dPLFkoNp2rAZr5403PC+zeUiRmwBuE622rduQQnKmz4Hh0MVv3CQeXiB7bn4Gtn9RKwppXMKRC2aeQJ9LOZWtXQaW70n7k+Wc6qvdgt2nRg7DwfOtholiBi3OpqE7tQhUsAZ6M3jVt/e4BU+D7goYJq4cRDk744xZnsScnW07VoMV05IM2uussREC7LgOa8XerLRyuZhrkOtMBVLBJmHs3OHM1sxy5cyIgw0fz7mZUW5Az8fyy+o2UFDEhwfhreW5iA8LdGP7hACFe8tR39YjgHH1uyeFfLjiniTsOt2IZ+4fi9gQfzx2TxL2n20WfrnlYCUopVg0ORbbv6hFy/U+vX2KYvHkWLy0twKHL36Dy529IGCAETncHwCQGBYAAFiWG4s/lzbApVI8/5ez2FhUjoUTY7D1SJVg+duP1WLRxFhQUHxwoh7zsqLENXnz0Rw8PiPJcI0Ki8pRWFSB+o5erC/IwKsfnzf40/oFmYgLDUC/S0NsiL+oCXC/NfvM2rnp2HasBgAEUeCbsnHZSAb9uvYeQPfPuvYe9tzcdDH2hu4pOtAdxK+zXFvgz1kV2uWa2fdht1rSuQ9AJaW09hZ/7qCWGjnMwERlLZUDh6xFy6yfT7yX97N0lhd+eCSXtT3O9hPCAgXIAwykuCTzjK4Xyox+Xla0oUD226WTRArOmY5cMPME+glhgQLkZZbfcr3fq45v7jKRC31mVir//eS7peycKMWKaUli0jwzZxyg/3Y3lqozIq6xysFT6K7SuMWHB7GJpQdR/pxZyuGTmLOztbokAcDA9uUADgzorhxY+cQFIIq1Vfp14d0V8kSVi258MvOgKWdIde09A9mTriNzXV3OLF/eXyHGc9+ZJjfiwcdy27EarCvIEDqvzPbr9CBZ9003+l0a4kIDsH5BJuLDBkDw2bnpiNLBubmzF6vfO4nNxZVCamJMnuLB3HjsP9OMeVnR2PZ5DRZOisGGPWX4oKQeq/NTcPjiFazJT8HBC1exdEoMDlR8AwDQKNPvQwJtIArBv+TGIdBPwaKJsdh1qgE9Lg3N19icm5cVjf1nm/FgdjxAgF2nGsWcOF7bjrePsHpC4d5y8Vt+NW+c+PfaeeMRG+KPfhfLIAr3lqOhoxcOO0FJbYfQ5D0BP2ftnCjMy4oW10EOAHImDUKEPl984SoDZX3s+d9WxVpzAJCPO1W2a6rIpiU14nbr98CtB/ylAHZI/36aEPI1IeT3hJAwqzcQQn5OCDlBCDnR2tpq9ZIhGSEEqZHDDADPtVQ+qFxD58Yn19q56dh3pkmwf/k5DmJxOohxts+dwdxpw0GeyxoUwHPzxhtYFS9syoGAM531CzIMQG8G/br2HlFIq2ztQkJYIJwuih8lhQkdXy6qcWu+1msJVAAEE7UC/nUFGfjdo9moa+8RjBpgYPHkuyct5QgeWImiuAVPoo8bH2NeO2FgFTgo8PPAzAFansRr57HOEDmAm3XX1MhhokCLIbZkyuYJ/LkeLNcBuNy38aMyPHP/WDGGa+emD/ytA5QV49/8SI4YJ87248ICDcAOAA47K3oW7i1HQ3sv7DaCkroOFBaVY9O+c1g1Mxlbj1Rj5fQkOBSC5k7meytnJOOpWWnYf6YZ8yfEoOjrRtydEoFdJxvQ62SfzQWGjm4nWm/0o1FnvsCAtPOfn9dB1SjeL6nDwomx2H26EatmpIBqFJsPVmJeVhT2nWnC/Akx2H+2GWvyU7FqZgr2nWnCvKxovHW4cqCNkQD17T3YVFTB5A393zLAN1/rFRnW+oIMbDtWYyB0de09bj7Dj3OisO9Mk8gOtx6pwoppRtDn2ZZbANAJAjBAFnh9oa69B5qmoY7XGwDLv/nct1IjbrfdMsAnhPgBWAhgp35oC4BUAJMANAP431bvo5S+SSnNpZTmRkZGfuvvN2j4siyhp0upkcMMLZJmHV+eXN6AXsgwoEJ/ZsVGTYAUBWtN5FIO17MBCFkkf8xIPDc/Qzj9qx+fFwBkBnrxvaZCWn0b6/Zx2Agar/UhK5p1KVGNAZehYEuIG1Ctmz9edAwAFsCv64pcO56XFc0YNSFuaTj/fZSycXnm/rHi309uL4WmaaIoK4+bvPisvqPXwPi/rG5zS8v3nWnC2rnpbmyfZ2zyhOYszqy7pkYOE4GI6638+0WhVi7aevhb7rZItajdABCsn+vZMvngBIP7I2f8/Dgfq0d/lICtR6qwbv54NOjXoqSmAw478xc5AMSFBugAX4Xn56Xj+YLxguWPHhGAJ/JSsO3zGsybEI03DlVh65EqzJ8Qg31nmnB3agTe+7IWeWMjEeyvYEpyODYXV2JiXBj+eLweD+XGobS6XVyPkUF29gcBfpweibYuJ979ohY9erCw2wh6XRo+KKnH/AnROuhHY0txFbYercL8CdH4y6kGrC/IwOr8FGw9Ws18klL0SVlLXGgA+nTGvL6AyaGrZiRjU1EF4kID8Ny88dD0YjAHd7PPcF/ggP7cfBYo6jt6AcLGRs6kud9W6X353sgCvwbcv3jgl0mA/Df/vFljI1DV2iWk4e+jF/9WMvx5AE5SSi8DAKX0MqVUpZRqAN4CMPUWfpebyRq+1YTiqb7cLbHtWI3b5Fp+d8KQgL6+rUfoe/K/ubbPNdTl0xKZFqnLN7yAVnW1G/FhgXhOBye5iDsY0BNAsHwKCFmn6VqvYUy+brwmCrY/Sg4Xv1Nm95uKKgzALwMTl8d4JiLLD/Udvdh+rFZkKLwV88vqNvHZT24vRZ2ud3OWVt/e46bbc02UT1SiMynO3OS0fPMjOYYinCctlgfTFdOS8PL+ChHEvU3oWWMjRCCQ+9/lIpzcDgkwFs8zy1cOsN5wzuhr21ggWzWDAfDauekGRgcM1Epkxs9/Cx9LSimIwgLMYAGgvqMXmw9WCtB9YU85XvioHKtmJOPXH5Vjc3ElVkxLws4T9eh3qVg0KZYB4oRo/LWsGcH+dnxacRn5Y0fjo6+aMS8rCiXVVxESYAelFDfUAf/65EIHAKDHBRSdacG8zNHwtxM4VQ1biivx0+xYUKrBqWm6bh+ND07UA6BYODEGH5yoh0ujaO7sxZZids6UUrzwUTmcqmZg9gA7zgF+1HB/gADHa9pRWFSBwqJyrJiWiF2nGtx85jkTAZBZfHxYoKgRyJk0z6hq27oNZMHKL9JGBQsfyksLF77ESYD8d/6YkaLVtqq1y9Bd9M/Wh78MkpxDCImWnvsJgLO38LvcTNbwAfcJJVrlPDB+PrmmJIXj9Yez3VYocuA/Xt1uSK25lm/+N++emJoYBkqp6DqQpR0O7rw9jmv5Znbf3NlrAPrCfRUC/BvaexAfGgCniwpNnRvVKCgFMqOHo+FaHwMEfRKZQf7Vj8+7AZMZ7J8zyQ8iCOgrKrnuvL4gQxRk44UEEWj5by7fcJ2Wd1HEhwXiufkZQo81M7fXlk3GlKQwMUHNgM4zqKVT4iz1ex6Iatu6Ratqtb4+4sntpSg+3yq0fbkIx/8uPt8qmBmvhzxz/1hW05EyprVz0zElKRxEIZZ98uZslPvs3EzWybJoIiumrivIYB/qIQC8dbgSq2YkIz4sAC5Vg11/3uliCD16hD+cvCaBAT/ZWcqAeOeJeqgUsClAr4ti/5kmDA+w49ilq+jqZ6D9x5JG2EzzLoAAeamh6OhVsf/sZUxNHgkboehxqXjvizpc61Fx//jRcKoa/ni8FiAED+bEY/fpRrg0ijX5Kdh6pAo9ThV2HY2cqib+5sz+xQfY769r68aGPWXYtK9C+Ouv5o3D+gWZmJIYBqIQQVq4z3gCfcC6HZZ3Dq2cnjyQjepFWt4RVHzhqvAb2S/kpgS5o89MCHhDhCwJ/tP04RNCggD8NwD/n3T4FULIGULI1wBmA/i3W/FdXs5BMK1f7jglgE2WKDwx/uV3JxoKZID7CkVeuH3n82qsmpECRerhlbslAGBdQQYadEZf39GLDQsy8L8+uSCyhoSwQDf93puMs/VIFVbOTMYrH593A/+NReWoa++Bn13xmBI2d/aJbhxZwtlUVIGV05NFN8XoEQEgCjG0aMpgz9sM5boHP/6v758WIB0XGiAyB3lBzXp9fGUtnwM/B8C69h5DjzQPgjy74OdwvKbdbbLKDO39kgbRLSVrsrzYFh8eJDRYGfw1TRso1lod0wty/JimaW4Z0/qCDJHux+mrK998NEfqk2cAbeWbXDrjmjf/3QDw6sfnPQYAFwW2Hq3G8Zp22G0Ea/JTkRAe6Pb3kuw4bD1ShadmpeLff5IFp0ax43gtnKqGpbnxUAiTAl0a+629qobOPhX9esYgEXwAQC8FDlcypj8pfgQ++roFLgqoKoVT0xASaMfH5S1wqRSgwJLJsdh1ugGLJ8UiyM+O0SMCsGpGCgIcCl5ayGo7dhvBU7NS9SK6hhc+KgchBH52Bc3XemFXCFbNSBEA39LZx7p49I4eXpTlAO8J9GvbuvHM/WNRuLd8oB6k+whRCKJCAkSRli/84sRm40dlqP2my6NfyMdqvxlokOCfw/xwoLAPfD86/i0BfEppN6V0JKX0mnRsOaV0AqX0LkrpQkpp8634Li/nYFj8IpiWDhzyhDIzfrlvmW+BIPc2yzLPz6anGFsqLSSeF3aXiYUmhbrO+MyccajXv5MvepEdUqMUK6YlWso4BATRIwIE+D8+IwmF+yqgUQqniyI+NADrFmQgOsTfcmyiR/iJDhmzhENBBSt69ePzhhZNPmlW3JMkFsyIoqIE7lyz5wH0eE27geUuvzsBhUXlhiX0/FrIk3T53QmCUUEfA7kvmmdmz+qLcmQGZ2Bo88aL580SDqAXU/U9T6zAHwASLI5RSkXhjQdX3jFiBfLyGg3+OjmwciY6FLCXnzMHgG3HarAmPxXPz0/HloOVUDXg7b9Xo6SmAyAEbxyuQn1bD1x6YfV6r4o3DlXhcmcf7AqBS6VwaRR/0lm+ohCoGkVnnwpVpQi0A32D3O7WTwFK6ztBAWRFj8D1PhXX+zTMz4oCAQFA4aIUO/Tv31nagFUzkvHCnnL89uAluDSKls5ebNphD2ToAAAgAElEQVR3Di6VYsuhKjR09MKuMMYv6/eLJ7GgJ/uOFZvnY+gJ9F85cA7N13oFqJt9hPuAXN+p7+gVoO1NfpF9BcDA1gqKAoediByLEwF5A8bbaXfMSltDP3R4EF79+DzWFTBJQJZ3rFLof33/tGCYMqs/Xs2AKzcxTJ9krKL/zP1jDa8xSzo2BVg5M1nIOcdr21FYVI6Ne8s9Ot+mfRV463DVgCNKTH79ggzB7on+HyhwubMPDhtBSW07Nu4tQ9O1Psuxab7WJ4CLOypn928fqTawfD42MvBsPVKFwiKjhOMN3Lcdq8GqGSkCUKckhmHVjBS2hF7PvFjRPAEb9ba75+ZniEU2vFgrAzlfZs+7c8wTVy4cGzX8RIOEYzWxZfDnRVmiKIYCHZ/svPCWMHKYVMi1Bnk+6SkomjoGirdMPksx1JBuFuxfOXAOc7Oi8JeT9Si4KwbvflGLlmt98HMoKFycicenJ+Otw5V4MDuWdeZc64VNB3eNaujpV7G5uBJP5afCZmNtny6VgmoULnWgpbDbRYVS6A0s+vWA4KcARyrbEGgnGO6nMHavsdZcqrHAwjOI0cP94FQ1Aeibiyvx/Px0/HJ2qtC1Ny7Kgt1G8MJH5dA0DSvuSRoUxLnf7DrVgBXTkrDrVINo+eREhfuGLN2YfWT1uyeF1DhrbIRbcZ+Dt1zM58esfIWDvCci8E/VpfOPNs5Y5fQIYGDOQYYHAyu9VF6KLmuo6woyRGGMBwCunS6/O8Eg8XAJ54UHMvHuF3UAIXhiZgq2H6vF8/PHY1VeCvafbcbyaYnYyFeYzk3HX041YOX0JKxfMNC1s+3zGiyaFItX9GKSzO63HjU+bj5UiX6XhrYuJ4hoMZTGxtQX3KSzGs7u+SP03zVXaqPbd6YJ6wsysHJmsiXblzV7Du6P3ZOMbcdqhGxDCMH2L2rxxMwUg4Sz9Ug1nKomgq0nxsaX2XsqxvFNwOS9V7iGz4OIlYQjgz9fzMMntlyU48f4KuFZYyNEKyl/DqBuIM8n88rpyXj7aDVWzUgWhdmlU+L0GgXPIlnB8bF7koYM9jtL6tHroth1uhFz9Q6qVTNTQAC8UVyJHpeGD082YqWudT+YHQebjQifMO8dREHRr2ro7NOgKATD/QgCbUCvruMMQvIBMOD3t7FAAQAqhTUTpmxtgMNGYFcIluay/vyWzl78/u+1cKmsSEspxZp8FgB+teusoavICvTnZkUJaZXLMkQh4vfWtfcYfEOWbiiAdfPHo3BvuZBrOLGTV4nnjxn5rXzFykc4EeD76txuu2MAP21UMN5aMcUtcsqFNHMw8DSZLCeYFAAACOYvSzwrpiXhVV2eeeb+sahv62b7gUxjm1Vt+5z1sU9NDEO/rk0CgFNjKSzgno7OzYwySDsyw3epQNRwfxBKsDovBbu/asSyqQk4UM6WwBOFpYiKwi5zc2cfKCjeOFQp2P3CiTHYfqxWgHa/qmHniXrMmzCg3fOVklZs/7F7kpmsU8SKaADwlt7XTCnFRh1g+UpHGdTXFYzHmvxUIefwIGjF2DY/kiNWr3KJjUs9vHgryzhcw6/v6MWbj+Z4lHD4xOb+cljarIwX5Xghzqogx4/xBVC137i3vHLQyU0Mw3PzMwydGoJE6IVY/pr4sAA9I4yy9M9dpxrgohSFizJET/v8CTHYergKv9p1FhQUD+XGA2CafI9Lw/sl9brrEgT42bAkOw5bii9B1TclUzWKbifFMIfCvIwQ9JhF+yFYn/4eVfUcIlwaxdtHq/FgTjwcNgXvfVkHl8pY/uPTk/D07AGQf/toNR7MjkeQw4ZFE2M9yjVr56ZjZ0k9+lW2IOzNR3OELLP9C7bS95UD59x8g7/m1Y/Ps/mojwdn5Xx1My/kD8Uf+KNcsK1v6zH4yKoZydj+Ra0IGv9sXTr/UONFW3PkBCFC3uFaGWeU8kpHYOgBQP48LvHUtXVj65EqptvvLcevdp1FYVEFVtyThC3FVeLvt45Uob69BxsXZopOio0LM0EBoXHLgLfzZD36XZoo0nJWv+VwJSgoWjr7QAhB1Ah/qBrBw7mxum4K3BUbgqdnp2L0CH8B9KvzUkBAmP6vUew8US9krlc/Po8lOXEAYDguj8GqGSlYVzAQlLYdq8HxmnY4VQ2biytR395j+LvfpWLDnjLEhQYweUWScBLCg7D9i1o8dk8ypiSGwalq4rWP3ZNsmNB8qTvvpFo7N12w9dXvnnRL1d86XIlHf5QgWkF5z74Z/I199VQEa7koJxdo5UItlx3k7IkvzJL9JCE8CG8+miMyxfdLGsQiK3MnjpxNOlXWzrjWQvZZPCkOQQ4bLt/oF8G46OtGLJoUiyCHDQ/mxGPP6Ub0OlVsOVTJOndsBEunxCPAoWBJdhx2nqgXC6zYybMHjbJCrUulCLhJhHAQIEhv5elWIYq93DSVwm4jsClgK3JPN+C+9FHo7HXh/vGjAQL89uAl/P7vNViSHacXaJOx/2wzFk+Ow/6zzQLsBSmSWL2f/tte/fi8W2GfSzuyb3D/4dJNfFggHHaChJHDsL4gA5uKKsQ1tfIDb8V8c8E2PjzIIAPKBICvxv6n6NL5r2C8aMv6nlPEBDJr+LJ8IxdrvQWAoq8b8dg9yaK7ZsW0JMSHBeCF3WVMEwfE37mJYVg1MwUBdkXo+HZJ0+/XQQ0A1heMZ7sTwhr0d59uhNOlYUl2HF4+cA5fN+g1cQrB6mXgB4CS2g7sL2tBwV0M9H/72SU8v+usAHrOOt84VImfZsfBz66g6VoPNu4tw9ysKOw/04wHc+NFR4RZW5azmb+casCKaazDaUl2HBw2BfFhgViTn+r2d317DzYXVxokHN6Bw5n/Swsz4bApOF7Tjm3HasSE5r30VoDd1DFQLPeUqnMQ5xPcSsKZPS4S6xdkigVwZj1/MH3W/BrekRUfFoDV754UGQRvofQm4fBAsCQ7Dn52GyMFEtvn47J4chy2f14jfLTgrlgc0IFx35km5I0dhWB/O5bkxMHfzoDwz6WNcFGKd7+sZZ04OsrbFfabhzkIelwDssxQZBzZnHRA+vFTgC4d8DWVQtMorvepmJMZBQDYUVKHG70qDl24god/lIDDF69gyeQ42BVg4cQY7D7dCJUCUSMCBkjQhGjDuM2bEG1g9RsWZGL/2WbMzYoyADtgLfFw/5FlRU+kwHyNvWn4xueM3VoJ+t5ICiGixVfesfZ22h0D+PLCK1kflSfbimmJhsklF2s9BQA+kTjQPXZPMrYfq8Hx6nY47AS/XsgkHruNYPUspuVvO1aDn0yOEzr+z/NSsf0Y6yrgwM7lHDkAvPRABpyqNrCgqmA8CCHYdboBd6dEYHPxJcxMi8SWw5WiUMqBf/OhS3BqjFmoGvDR1y34qpFtXet0aYga4Q8XpdhSfAk/zWYs/sPSeiyaFKv3QQ9IOfvPNGPxpDhDe5ucQhcWleM/PrsIl0aFZr//bDM2LMhgC7J0vd6s3TtsDHSsNHre4rphQYbYwoHXDKYkhroVb3mbJ1t1mTLkVJ23W/Lrzm8I4k2jNT/Kr+OFWh4w+Gt4my6XejRNg1NlhUqXBkSN8HeTcMxa/v6zzVi/IEM/d8b2+YrVx6YnY9fpBiG9zb8rxvDvaamReL+kDnljI/Hn0kZolODDk0wGmjN+NDp7XLg/YzTsCoHdRjAnIwo3+jQoBBjup0Av6aBfA+w3KS33qkCgnb13hJ+C+VlR6OxToSgED0+Nx6cVl6FqrE1TUYCfTI7DoQtX0NWnYWdpPfpcFLtPs2yFapStO6EUy6clYWdJPZwaFStjd5bUAwSC1XM9f9epBrg0Y63CLPGsmsFqK7xpQWb9VqTA06OVhm9+5Po9zyj+UTr+HQP4/I5JZn1UHlCzPjo3c7ShWGsOALmJofjZ9GSxWKpeXw3JC2S8IPnKgXMM4I/V6nuLJ2HX6QYh9bx1uFJo8fFhgXgqP9VSzmnp7EOQw4bFk5jW2HJ94N+HLlzGmlmp+PT8ZdbPnB1rAH6Xi7rdvYNvsQAw/d6uEPS5NOw8UYc1+alwaRQflNRhfUEGnp7F/r3rVCMeu4d1NggwMaXQiyYy2eClhZmob+/BW0fY6k2+vH1eVhTeOlKJOn0xFJexNizIMKTkntYhyLr8loOVePK9gdXH5nZM/tpn54y7qVSdT3D+Hm8aradHvviKUip8jgcQ7nvN13otmV98eJBYr+FJy187bzxAKZPasuMAAiHxRI1grYE8COw61Yju/oGg8PmlK1g2NRGHz1+Bi1L8YlYKVuelABQ4dOEKCu6KwZGLV6FqQL9K8em5ywj2t4HorYE9LiDAxv53UcbW/bzMP38JSQJsgI2wLh2bTcEnFZcxXP/sj8tb0OeisBGCh6bEg4Bg1+lG5I0ZhWB/m8hGVs1IwYcnGmC3ESyaGIuNReX47acXDeDefK0Xfg4FD+bEC1bPN1vb8EAWXlzIgH3hxBhsKqow+IQ5EyQKk0Vl6YYTxcH8wJuGb67xyFKPWcfnN3C6nXbHAD7fj4WvZOPLz/mAWumju041uqXScgAw7y5opefz7RL2n2nC3KxoFO7V2a8KIfV0O1XsKKmD08Wi+tYjVVg0MQaFRYzlc9DfXHwJiybF4i+nGjAvKwrv/L0aeWNHYfepRvQ4NUQG+8NPUfBgThx2nmwQwL/50CWolKXlV2/0izHhC4IoIJh9gMMGVWc+vDDX0tmHt/9eDQBYOSMZo4L90OPU8MGJOmkPFA760dh1ugFP5DHWvmFPGbr7VfzHZxfF8vYPTtSju18VmQvX5uX6hDfQ59orB8uFd8W4dWOs0GslfBIDMMg4g6XqhonOuzQ86LCDPXKddkCuGWhN5b7HdXxZz+drEzz5Kn9++d2J2H26EU/lpwqJp7CoAkuyY0UQWF8wHr+8N1X8u+CuWBy+eAWLJ8fCbiN4rbgSbxyuAggwa9xo7DvThIlxoaCguN7rQm5iOBQFuNHHNOm8tHD0qoyt+9sYWx/wLHdTCBBkZ2DfqwI3nKweouqtmH52Ak2jcKoabArLtv58sgEgQN4YPRsZNwp/Lm3UWS5j5g/mMN1+0cRYBPopWJITh92nG7H8btZS+cRM1vn22D1JQto5XtMupDnO9p0qdfMJniVyFp6ga+yyHENB0djRPagfmGs65tcAsJR6ZB2f3+P6dtodA/j8Zga8I4dNYoLRIwIMMg3XR19cmIn1D2R67YZYeFeM2F3QDPaUUmz8qAzzJ0Rj/5kmzJOLZn42g9SzdEq84Vi/xhagyKC4cSHbYGlHSZ0OtvWYlhqBHcdrMWMM21Ru69EqLJoUK8D+wZw4fFBaj95+tpR+dX4KdpTUCaK/r6wVy3Lj4G+3odfJmP2SnHjYbQS/PXgJGtWwJj8VW49UYeFdMQCYtv/87jJo0NDXr+EvpxoMYP/BiXr0OAeKsg6bgqW5rIOCr370s9uwNJd1X1AALy3MBAXFhj1lhl5qT6D/zufVmJM52jCh5de8/nA2Ro8wArunlP1mUnXzhB/qI5+85m6ch3Jj3aRFWWJ8/J4kANbyzsa95Vh4VzQbQ8r62KNCAvHWcrbJn1PVsLO0AQ/qOj8IwdtHa/Cgngn88XgtuvtV8ZpAu4J7x42CXSH4W8VlPZC3gAAomBCF/WcvQ9OA+Vmj0dWv4dClNgxzKCjIHCW6bgDG5ANtTLIZGWTD/IxIBNmBHhXodjGwDwlQkD8mHF1OtnjLrhDcnTJS/D03Mxo3+ly4N30UKKX4W/llPDQlAX8raxEMfuvRavQ5New63Yh5E6Kx+3QjnpiZig9PNMCprwImCsGo4f7iRil+DgU/zY41XHMz25d94oXdZdhYVC6yMlkKlH3p7aOMDHm6/kOp6fAN1sw+x6W/4vOt+Pm2Ex7vqXyr7I4BfKIXTuTOiPUF47Fp3zlLDZ/or/XE9vn+3avy2F4fckCglOKF3WXoEQ4ZM1A0K2sR2/8W7i3H4kmxurOmICGctRU+mB2LQIcNS6fEGwq1T+WniOMgBH8tb8GyqQk4epHdrGPRxFjsPFnvBvaEEDhs7FI6nQOzkwIYEWiDn8MGAhhA32FjrXeUUvS4NPa5AFbnpeh90QkYHmjH4kmx+MspaeMrQgSYx4cF4omZydj9FQt0W49Woa6tG6tmJuNAWQsWT4oVTGvjwizLXmpv7XVcq7Vi9Zv0NlB550RvKTvPdCgovqxps0y7AQxJsx3skU9iT+k81/OjQwLgZ1cs5R0O6E/MTBa3ARS+LLH73aebDDo/B/gg4UfAztIG5I8bhT+dqMd96aNBFIIvq66i4K4oBPrZUVLThoK7otDVp+Lvl65i2ZQ4hAQq6HZq+HvlNyjIHIUQf4IgO0GfxsC91wXMy4rB/vJWdOvyT7ADCA1QMDMtEocvtmGYgyBimB13p4xE0ZnLCA2yY9lUVpxdOjUBn1Ww23H2qRQHzjaDKEQw+FUzUoS0w32BUgo/h4I1+alCe//1nnJsLCrXtzDJ1PfoGdDuv6xpw6aiCtHFJXxiBFuV/vg9SZZSoOxLGx7Iwu9/luOWpfGN18w1HXPQkDu0eHAxS3/cN28vv7+DAB8Y0PE5o6KUig2jPGn4ZrZvDgyyVsrb4wqLyuHUNNEDzcF+/5kmdt9OTRMBYWdpPW70qWyfk+p29GsUO0sbRCB4MDsWFMDzu8uw9Ug1Fk+KxZ6vmrA6LxmqrrcunBQLQgg+PFkvunY+KK1Hn5PCrrAWt5/q0o55Cfz2Y7VMyvFzB33B9DUNS7JZFkIIgUqBPV814SeT47CztAE9Tmnjq+xY7P6qEatmJOvyVLV4vrtfNQL61434UfJIvKLvJf7SA5mizc6y40Ja7u7nULBmVqooxJlZPQhE9saP80ljlnHMjG3rkSrWTqe/3pzqv3+i4aa0/MH0Wp6mE6mTw25jBMWTvLMmPwV+dpvQ8tkWx+xcObC/9IDkt1IQ2FnagCfyUrHnqybB+P9a3oKnZqXh8MUreDA7DirYDpdLsuPg0oCSagb6nb0q9p9twYofJSEk0IaOXhX7yq4AREGgn4IgO2P4w/1tmJ4aiuF+ipBxHHY7NEpQdPYyhvkpGB5ox7TUCAH2j/4oAX8ubUSfCzhwthkqgF/OToNNYVt8L8mJEwweAPz9WH8nZ+5vH63G+oIMQ6eNw07wwF3RHtn8loOV6HOpIsgDENmcn0NBdGigpX4v6/ibiipQ39bjlqU9/k4pNuwpc6vpmIOGe4dWgriBPF/bwzt3breOf8cAPqUDtxfkUROEuE0sKw1fZvsUMMhAvGDmZ1OEtMOLlrwHWu6Y4AtfAKBwcSbWzEpFsJ8NCyfG4J3Pq/UJOCDd7CxtwJLJsVAIxSI9CHT3qyAgomD74akG2BW28RXA2tlcKoVNAWNxYCl8T58KYgJ8lwrsLKllPc02BjS8cLskOw4EbMHNnq+asHhSLLYcqoRNIVg8KUZi9HFsQk6OFQGAa/YrZyZj0yLWTvnL2amGxTF3p0Zix/Fa/CglghWhO/vg77AhKiTQjcXLfdJbj1RhfUEGpiaFsz2FBmH1fH8gDujyvjyDMTa504e9j/Vpc6AeiobPO394ALFq0ZP90Cr7NMs7U5LCsX4BY/2FRWyLY75ilwcDEILCogocr24zSjw2BRSsd3xnaQPW5KWCEIKIYNYhtLO0Acty49n6DErhsBH0qRQlNW1YNjUeHT0uvPdlHf7tx2MQEmBDsJ8Ch43g4amJ6HYxhn+jT0V5cxeu92voVYER/gr++31pcNgIhvvbEBxgw5Skkdj7dYsB7EGAH4+PBED1MQH87DY8NSsNH55gmj4H94EVyikiEDRd6zH4wRMz2XP9qmbJ5nlGwIF/qPq9WccXDNzgEyr6Xeqg/uBpAZ5Z8gPg69IZqvG2THPUNE8sXlAza/ic7RtSZ53tm6WdgX+z5f6jh/ujp58FkJUzk2FXFLy4kO0M+faRaiyeHIddpxpQcBcr7t43bhQC7QpbaALg/RP1UDW2Ve2iSTGw2xRsOVTJ9PpSxuofzInHnq+axS6G/L0c/FUN0CjbvVC2Xg3o6mcMnIJiaW4869ZxathxvA5OjUIhwKKJMfigpB5OTRPAzhn9hydZmryjpE4EALsCLLyLrbwdWPwVILaPmD8hBscqW7FsaiKOVbaKm6csmhhr6LCQOy7kQmrztV5RuJVZvazZiq0BCNB0rcdNzpFZP5+oZsbGX88n30O5sSIgDNZ/PZg+O2tsBPPDsECDH8rZpzd5h7P+x+9JYr+Dam5afr+L7YkjSzyr8lLw9tEa1pUjCqBsE7XVeSlwqhQfnmzEgznx2HK4CrmJYfC3s8xu7OhhCAm0o6PHhYuXbyDI3w5CAIUomBA7AqGBdgQ7CEKC7MiKGY6QQDuG+ykYFmAHuxUj4G8nuC89Cnu/bsGIABv+7b4x+LCUbYc8e9wo7DjeAEIYgXmjuEr4qp9DMbB8rtNHjfAHIfq6k4OV6Fc14Qe5CaGG58xsfsMDWZiaNAD8Zv2eZ3PA0OQ8Q2CwMSn1ZvzBm+R3uwu2wB0E+FzOMUdNK1mH39jESsOPGhEAnjrLbJ9LO8bXsf1ANu2rAAXFugJ2nBCgpbMPhUXluN7nErLOrtONuFvvj56UEI4/nagXhbSHcuLg1DTsPMEYv1PTAVkH8w9K6jAjbSRsBHCqwI6SBvQ7NbZXCShUTQOBZmiPA1j6TanGdkRUWbHYqekdFJSCUEDTgB3Ha3Gj3wVK+WuMwC+zfDkA8A4dp6rpslQV5k2IZrr/hBgcvnAFP5tulHD4Yq3Fk1gHBt/EirP0dfPHi4kpd6+smmHUbK1kmsEYnGBsUgAwyzI8IPDMgWu1Q53Ahy5cxca95fjTiQYUFpXj8W2lBj80ZwGe5B0O9AAFAfD20RqPWr4s8UxJDNVX+7L3bTlUhTV5KbpHMEbv0ihGDrNjRIADe79uwYM58aAUeONwNf7tx2MQHGDHH4834L5xo3C9T0OfS8OZxk4E+tmgKPr2HoQBiN2mYHVeCn7ztwto73FBpQQHzjZjRIANQf42gLAsEmBSTmiQA7+cnYadJ9i6gNV5KYLJm4+tm8+29vBzsKyFgzv3g5LaDibNhAS4Ab/sCxsWZGJq0kBNx2EnWDkjWWRzQ5XzgIHA8NKiLLy0cGj+4AngZV+Aby+doZvcB7txbzmO63essZJ1zHrpiw9kGLTQRRNjpEJZBY7XtOHVj8/jxYWZWLcgA9D7zRdPisX2z2uExNNyvQ+b9lWg16myjpqJsQj2Y73FwX42LJoUK/qjyxrbsWxqIo5cbEWfS2OLYlwUvS4V735ZK5xh2ZQE+DtYAe294w1wuigcNmDZlDi214dGoaoUnX0aXBpx38aWANf7eYucBlVjIG+36TsngkKxAf0qcKNXxZzxbBUkKATw3zsuEru/akKfk+LD0ga4NIqnZ7EVtPKjXSFYOT1ZbOq180QdelwaRg/39yjhrJ2bbtgumLN6rtGaOzFkzdYs07y4cEDOMXdgDBoAAJhTd75Qy6zVDkWz73epeOsw24mRrWIlbuzPiv3L8k50SIAAep6RTkkMFZkB78Lh7Zpc4impaQfA9mdaLYCe2ZZDVWwzMh3c//t9aRgZ7EBEsJ8IBIQAw/wUhAY5MGZUEEID7bjW68LrxVW4b9wodPZp6Oh24mzjdbT3uHCtx4Wr1/vR0aNiuL+C//veVLF//VP5adhSzNpBn56dCoUQBPgpIgjYbWz8XRowerif2zG+sytfJEX06wWwgisnCjwoRIcGemTzsg9seCALS3PjTH4jyXmwknCMBIHfTWwwQDf7lhmTvi/9HriDAF8e1H6Xiq1Hqj3KOma91KyFyiy+36Xh9YOX8OyccSCEYNO+CmzYUyYWV82bYOzmWTQxRoD7vjNN+El2PA6caUb+uNHY93UT7kkbhWOVreIxb+woBPvZcF/6aAQ4FKgaRUe3C9OSI+CwKfjwZAP+JScOAQ6CYX4sZVQ1YH8ZW4ClUablD3ewjRP8zCSBsuXyoMD1PrYtraIAczJG44beb7csNwF+uvb6cXkLHHYblk2JByjbC2XH8XrkjYmEn53gp9mxDMAwoDcaH1knxdOz01g7pr5Skks4P82OxSZ9OwmiEKHJyqsdvbF63okh91VbabecwbHi2OABwFPqDkJAKZ+88DiBzWz9D4/lYv0DmVg6JQ4bHjCyP2/sHwSIDmWrkuXPSQgPFJkHBdtV0q6wMZ+SFI4NOmHp6nNhc3GlRzmHH+Pgzpnx7w5Xi0Cw5VAVnpqVBptC8LsjNULLDwmyY2ZaOEID7QZJh/8dGuSAv8OGmNBAsYCL7fHEQJwQwGG3YU1eqggCa/JS8UZxFdj2IL2w2xTDsaZrPWzLkoRQAc5vH62Gi1K2elwKCjwbs2Tz05PdZByezXG/eSg3FusXZOLFhYP35nOJzcovrIL7UCQ+Xx/+TVjaqGBx0wI+SWRZp9+loamj221CvfmoUQs1s/0l2bGG/UwW3hUNh76lsNy9wyUf3iGx+3Qj5k9gy93v1vvpM2JC3R65vPN+ST2mJI1ESIAd+WMjsf9sC3ISwoTM8+N01h/939JHwaVqaOt2IjM2BNf72FYKGthiF2K6oj0q0OVkzhjkILokRPFJ+WWMCLBD02+KodKBlrAlk2PxYSnbx8RuA5ZNicffyi+jz0Xx55ON6HNqeONQFXr7VbfHrUerkTdmFN49VoPL1/sR5M8K1gzEBwpsFGwrWq7JcjZvZvW5CaEAGJu7mU4MzuAGCwCeUnhg8NTd6pFnhfXtLJBZsb+uPhdeP3gJK6cnubF/MyDwz6n7pht9Lg3N13oAULx9tIad16IsEURKatphV9hOlFEj/KGq7HVr8lMAMPDlNxdZk58Ku01BS6e+T74OoMDA+/nfrL5gg8OmQFEUOGyK9d8KgcOu6PdtYFkk2wKZgfiWQ/DpCJQAACAASURBVFUD5yYxeSHFHKoCITDKM8VVcGqaGAfWCpyKQIdNEA05G1v5Tqklm5+aFGaQcbr6XZas/dd7ylG4r2JI2v5gfiFLO4NJfHXtPb7N027G+MVf/e5JwyThq9z4JLFi/ZuKKjx2Plj1Oa+ayfRFORiI4CB1SPBWRd4Sd6ruGyydkiAeT9ZexfwJ0Tjb0IaZYyKx72wzJieG4/CFVmQnhmHf2Rb09WvodWnYd7YFAPDRmcu43q8hwE5w6OI3CLKznS+7nBQBdsChwP2+ozZA1TR0OSn6nBqu9aro7tegUrb3eZ+L6fvX+1T0u5iG36s7o6bpy+FVCoVQ/DQ7Fn52gtnpkRjmb8O96QNL4oMDbMgbMwo7jtdiaspIA8jLfdTmpfG8K8OK1QuNNjTQY181b+E06/jypPUUANxSeP2R30PBW+pupdlvPVIFp8t7NmBm5jJAbNhdZgAEeaUm92Eu7ySEBwpf5YSFb1ZHCEGAQ8GqvBRMTQoDMJAVcNDl4M+3+mjR96f/tiaAXWLvW4qrJDZOYVdgDALS84QQdu9aLs+EBEDVWAbbdK1HgDNrzxzvVlQHGPg/fk8SNuwuM8hwsg+IgGHB2u02GPxC9gPz42B+we+Cx1UHbxJf/piReHZuOlIigr71+A/F7Lf1079nk7eqNWtlf3gsF/UdvYL1/2reOCnCQnQ+1LZ1i8nz0sJMduejsAD8fPtJLMmOxe7TTQLQd5Y2YOOiLFBIi15y4kUw2HWqEYsmx2LP6SZQSlkxq6wZhBDsP9sElwYUnWlGwYQY7P2qCUH+Npys/QZ5YyNx6EIrAh1sEyuWOgLzMyNRVNaK/LRwHLrUhkA723zqej+Fn8IWwwDGKG7DwA0s/BTG+PPTwnDoUjuGOwiGORQohBWGh/uzv0HZXuaKAoCyIrFDYVLSn0+ybR7eL2nAgrti8KcT9Vg6NUFffBWH3aca8NSsNHxQUgfFpuBKZy/87AoWTophoD6TSV/y4+aDlwysPjZkYFHM1iNVWL8gQ+xh/+JHZ6FqBF/WtOEPf6/FG49MxoYHspCXFo7okECAUrx84BweuycZbx6+ZHjt8/PTsWnfOcN7okYwnTx/bAQUJRN5aeFQlEzEhfgbUnfoO41a6v0YmNDrFrDurPwxIw2fl5cWDkIyEB8agP/82RTUd/QKBvj8/HQU7juHldMTQalm2I/dLBNxH37y3VNYOT3RAPRvHK5CdEiAyF43LMhA4b5ziBo+Dk5VE8x+y+EqgyQUFcKyU/78G4erMFSTmTyXcHgdQP4b0ANCfgreOjIA8qpGYbcTcTc3gJGnJ2YmIyE8CEH+bAU39xEu8VBKsWnfOawrGM/uNlfAbrcpyAEoKIXYFqGurRv/vv88np+fLgIGf1+e4VoZ/YLfsMfy2g/iF3FhgWzH3tAAce3yx4xETFiQ8Auena2akYTtX9QhITz7tmr5dwzDBzDk1JgzJj6Rhsr2B9g9a3ljbWjU0C0hL2lfPDkOe75qQt7YSGw5VIl7x41CkMOG+8ePhkNhmvtDOXE4VfcNlk2NR7C/gskJjOH72xUEOtjWBJoGXOtVUXzhGxAAJ+s6MMyhoNfF2BowcIs5wLiHmgq2R7n8mpKadhCwFY5dTsbwb/Rr0DTAbldwf8Zodj9Tjd2xKMChYE5mFIID7JidHokAB9su4mTtVbEkvtdJsePLWrgoC7z+DhvuHReJ14tZe6m51U6jmijUeeu84BqtuTjL2by8atasx5pfa8XezSm8laSzcfEE/OFnuV67dV5cmIn1CzKE1u6tY4f7opkBbj1SjadmpSHIz27Qf80yEZeFGNAzX5VlSu7PnAA1X+uxlHq4JJQQHghKqSTlDJ3ty0zdqXnfTJmAYEpimEGSMb+fUoqefhVvHakydNfwLO7to9Xidd5knI2LJ+ClhZliWwSZxfMNB70V47lfeKv1DEXS4RmbJ4mP3RSFkYVn56b/82yPTAipIYScIYScJoSc0I+FE0L+Sgi5qD+G3arvs7LUyGGGHmZP7XBmoDcXzfiiF+jdFvLEGq2vyGOFL6ZX8mCwJi/VcLs2WdJZOoUtKedL3KckR+B6r4q9Xzejq0/Dx2Ut6Omn2Pt1C3ISQ9Hn0pAUMQzXe1WoVIO/ne0tnhkdjK5+/T6gAFRpnxOerpnLPk7K2D3A9kLpdjHpp18DMqOD0euiyIgORpdTQ28/xcdlzRjmp8ClqrjRp2J81Aj86UQ9UkcNx47jDeh1UnxS1oJ+leBjfUn8j8dHwm4j+OnkWAHyf6u4jJHBfggfZhegzhfUDBTq+gS4A+6dF3LKLcs5bmBu0WbJX2suxMmszZzCi6CgL9I7fIl1aDV09GH1eydRWFSBw5faxHP8cdO+Cp1Y9Hpk/zydNzPAhLBAtyKvXNzdeqQK3f2qQd6xkoV4YJD9WWb2ZqlHLgQb/MX0Hm8ms3qrgqOcAQBgNw06Wi0Y/EBgH3g/7/bisgzgRY6BtYwj+8nvf5bjFqzjwwK9FuO5XxQWlQ96za2ee+WA7gtevkPTj/NrP2tsxG1feHWrJZ3ZlNKr0r+fA/AppfT/JYQ8p/977S3+TgADK20Li8oByi7W1qPVbqw/PjSA3Xhk3zk8r8s6ANy2rC3cV4F+l4bVeSl443AVpiSGITo0CFTTwFvenpqVhq1Hq/UbkVQBoCCKgj+fatSLVKxb4rcHK4WU89fyFqzJT2V7z2dFoehMC4L9AT87654I8rOhpLYDE2JG4ETdNfjbCBRC0Oui8LMBZ5tvICHUH3Ud7IYn/dIcc3kZH87u+Q2I+Nwsa74Bh8Ies2KG42zTdSjEhhv6G6YkhaL4wlU8cFc09nzdjBB/G0Ap7DYF96ZH4rNzV7AkOw6vHaxExHA/RAT7ISokAOHD7FAUBU/lpWBLcaW4d66Z1b99tBrrF7DbK/Jea9F5oafsg8k58WEBqGvrsWyzrG/vwa/3lMNhJ2juTMbvj9bgjUcmY/0C6xRelnYISUdeWjgwlz0+P5+t35CPmZ/LHzMSZIFRynlxYab4bHM6v2F3GQgBXnwgA5v2nQMkeYcHilUzU7D1aLVB3pFlIS4J8PsacICRZZvo0ECxsOsJfdyjhvuLQjDv2uFzSezZj+8GQA69R39zcSU0/bM5gx89IkNIMQ4bEZ09IoioGuzEKMmZ5RgrGUcO/lzG2bCn3DDObzwyGRsXT0B8WAAoxaDSjrdrbvUcnTPOIN9ZSXwtnSxAr5vPiGdieBDSRg//TuM9mN1uSWcRgHf0v98BsPh2fRFfafv49CSsKxhvaIfzxPoBC6CXFr0AEJpoSV2HWOLO5ZzRw/28LnDxs9tw9Ua/SI0fyonTbxZO0drVj2NVV/GLWWyTqKW5bG8buwIMD7Dh66ZOZEYHo19lzjvcX0G/yph5XUfft75wTHVkrB9gDuDUgDGRQTjbdB0AMGbUMBAAuYkhKKnpQFiQAwlhAYgJ8dc33GK30NtxvAEqZbfEGxnswC9mseC4Oj8FvztULcYBOribWb2nXmsABq3WSs7p6ncZVs16WiQVHxboXojj7P3iN4NKO0Mt0t3Mc/y+uZ4YIJd3eHFW9mVP8g7r+mEZBG+PlWUbeR7wjMEqAzAbDxw3U8z1pOt7Y/BcNpMLsUH+NjyVn4bNBwe6asxrIwB36c2cGVq1UNZ39LLMrb3Xq7Tzba/5YFLO4YvfWBb4b7fdSoZPAXxCCKEAfkcpfRPAaEppMwBQSpsJIaNu4fcZLCk8AFOTwrHlUCUC7HYQJYPdRWlu+qCs34rRy8Wv9QXjsbGoAiunJwkm77ApaLneB872ebGLb8K05VAVHsyJx+vFl/BUfip2lNTj/RP1IITo8k8K/nSiAQDQ3uPE9i/q4OdQQAiwIHM0dpxoRFnzDQTaCXpcFCmRgShv6YJTvwOR61t6h78CsThLAQsAaZFBuNjaDYB9dml9J6YkhaGkph1hQXbckxKOzYeq8fCP4vHHL+sRMdxPgPxDOXF4vbgKEcP98M0NJ1pv9OOKPtFZq12lWCwjs3q5WCv3Wkc9kAVQipf2lgEghu6Mf99/Hm88MhlPzEwV19Gsy658p1R04HA2Lxd1OaN7ds44xIUGWEo7quahEOetSHcTz8WGBuAZLwyQZwTx4YFY/e4pwfrdWbyx60dm8VafsWFPOQAqfBywZvMEEFsG8MDBs1m+pQAhpi0mpGNB/jasmpmCtw5Xio3QZONBwKlqIBoMTJwXbGeNjYCiZCEuxB+bCWAnpl54fUzNBXn3z2BjasXmCUlnxflbdF2H+pxG2ULINx/NMVz7263fA7cW8KdTSpt0UP8rIeTcUN5ECPk5gJ8DQEJCwrf+8j8er8P2L+sQEmDDxn9hjCYhPBspEUGIX54LAiAlchjuTo1ASkQQkiKCAUrdgoG5+CWkH4nJv320RuylI8s5PAjwAllmdDDCh/mJVYxOleIX+Sl4/RCbbG3d/Xj3i1o8lBuHHccbkJcQiiMXv8H+sy1YmhuDHSeakJsQgiNVHShv6UJG1DCUt3R9a7AnGAD7MH+C9j4KhwJc0sE+LzUUhys7EOxQcL6xHSEBNjw8JQ6vH6qBAiDE34bQIAf+JScWm4urETHcDxQUI4MdeFovqD2Vn4o/ldSBEtMeKCZZJ2qEP/wciqED43htu0Gq8QT8nros4kL83YC7vr1HAAF/5On971fkDEna4Y8vLcq6Jc95knL4+f16Tzn87AqeLxhvYP2yvMOLup4+/wX9N/56wXjDZl+8e4bf/MYKvAP9bHgiL0UUo81ylNVvNR5jYwkKz0GCEJEJ8GuraZqQeyggArYVWG9cPEH4hzlAO1UNW49UIzqEta1aXXtZ2rlV13Woz5mlHH7tE5bn3HZJ55YBPqW0SX+8Qgj5C4CpAC4TQqJ1dh8N4IrF+94E8CYA5ObmfuusplW/iTfVKBJGDoOiKKK9aYw0iPxY2qhgUErxpkUwUBQFrxw4h2fmjAOXJOQeaEII4kL98UaxqZ9ZZ0Gc7SuKArtNwe+OVIt2uG/+//bOPLqt6zr3vw2QojiJpEjKGihaoizJlmeLjtMklR1ncjzEzfBax22cvnQtPTd1m+Q1devKjUEn6nOctnkvr4kdJXXr+iWxUrue5HkemsQxJVuWRI2UZVESNYviIE7APe+PewGBIACCxHgOzrcWFsB9we+efQHs+5299z13YITjAyP84q293NA6n1+81UVdeSnV00t4decxrjnvDJ7afAhB8AGv7+7h08saeLrjKJ2HBij3u62V82aUsr93FID66XBsKPnx8eN27ADMrfZzoM/969L5M/jV+70AVE4rQYDPXDibn7cfwDcaYs+xUwhww6Xz+Nlv90Zy+w1V08aoewRCSjhvbjW/UK5iDrfaxQ/87g2nm2qn81Mv8P/w5V04MEbNxwZ+RFJS80nb6jz1FZ7WE76n7lVns/qprV4wPGdcsLgjJlgk23b7Nedwt8cd4fRmnX/1qaVJFaDfB1/5yIIx9YDYwDsmd//Ge4yGFG3XLUva3hkbdOfPrBgTmMJto+F2zjkzpsf1OVVbuDgbTrdEB+2m2uljZgIQPtl434WXxn4XEgXrv7v23HGBNXysomdwiZR3V8/QuM8n8n341NJJfear4py0x3FedTY/fX135I54sePLNjIS8EWkEvAppfq8158E7gQeB74M3OU9P5aJ/cXDVy9fwCvbj3DX55al3McqIpw1q4pdh/sjd5rpPNzPvBnT+MPL5jOvZhqrrnbX1w9c53btdPe4OdA/+d0WQirEvuMD9A8NM93vdrUMjgQ5dPJ018Od1y3jjic63KlzyOHnv93L718yjxe2HWZRQznVZSX8+6/fZ1qpjxnTS1h+Zg1PbT7EU5u6+YPWuTzYfoCa6aUADCv49JJ6nt56jP29o5EgXltewjGvCb+E+MXbpQ2ldBx1TxDNtdM50DdAmQ82dLnBfsVZdTy99aib5lFQV17C1eedwc/e2k/t9BLqKqbRO+Rw4wfm89D6/YhPuKCphtm1ZZF1S8LT93ChLhzcoztxwoH/jsc78PmE+25azreuGzt1j5fGif1hT6Tmw4XaVVefM64gFwly3rQ+XiEulWCRStqmqaaMH//hJXT1DI3ZTzwFGC/gJgsusao/tr0zPBMINytEp49E3FTm6qe2Rk4W34qaFYzGzArGLgGcmu3UiNvhNjvKj+jZVfxrIcZ+F2Lz8GOC9YnBuMclpc8+6rhzzbIxqb7w55ROuid80o4t5DbVlI1L5YTVv04pnTOAR7yWohLg50qpZ0TkLeCXIvInwF7gv2Vof+Pwv57dxabuPtau7+bO+Q0p/1/n4X7++F/fIhg63UMcDDn0DI5SW1FKic83zl49vZQfvryLoKP4pxc7OeXGUX7w0i5OnAryw1ffo6a8hDue6KDtumWMhhT/96VdjAQdeodDrNvUzcCIw/df2EWvp5ivXdLIsx2HOD4QpKa8hJODQeoryxDcG1VUlwp9o4oz6yuAYwB88hw3+HeeCHLOrHK2Hh6kaQbscWM4zVWw17tj2jlzaug46jZQvd3lXr792Yvn8OD6bgCWza7itV0n+PR5s3hw/QH8AkvOqGZmRSl+Hzy8YR+N1WVcvrSRl7YfRkQ4s6GKf//KB2lpqGBubQXfXtfBwd4hiEnn/PSN91i5omVsG15MB0Z46j5GzUMk8Aee2Dqu0yaZmo9uq/OJj9uvXTZOxXf3LuC+N/bEVW1j0z2ZT9vEmxF858kOBGHVNefABMGlqa48brol3kwg3jhi8/nRQXvMDVsmyNcntHnXqIy/qYi7j2jlHnssEqVZYm2JAnLCz34CxZ7JdM9Br1Y1J87JLswvUbPAu5/ZlvWLriBDAV8ptRu4MI79GPCxTOxjItywfC7/0d7FhxbWsqP7JF0nBt2bVcvpXt1YG8Ceo/0MDo96XxXXHgyGCCkYGQ3i+E8vVBC2h0JBRrwFDIZH3NRIcNTB8QtVpV7aweuAUEpFLmpxHLfP23EcKqf5EVQkRXPt+WfwO2c1cv+v3uN/fmIJ97y2mwvm10aC/xcvbeLB9n30Dp0+MS1sqCIc/C+YV8O2w4MsaqxmT6/XbXNGNXv73dclpe5HHb5KV4D6yrIIV4m/hBIffGDBTN587zgiwty6CmaUl0SWo/CJ0NJYyYKGKgRYNKsqciyvPHsWAtz+6CZExnZoDI0649rwwldQxkvjhNX8/LrySJEuWs0PBdWYgBFf0SVuq7v9mmUJ878ppXviTNMnm7aJNyOI19oZe3VuWBVG5+mTpVbcK8pDlJb4xu0/utCbSsonWaCLF4yTnRjiKveodso7ntg6oW+JWyvjf/aTVexJ0z0TfB/ipm3ifOa3RrVzZntZBTBoaYX9vcMMBRWrHnNXYoxV6IlU+2gwxInBINVlPvy+1FvP4kFEKCmJXckm/vvcfQklPoWEFH6/nxsva+KylnpaGir4nbMaaWmoYG1DFftODLJicT3L5tXyb/+1m7/46CIe2rCPi86cSf36fTgI58yrwb/xINOmlQFukI9+fc6cGmZVH+fGy5rZ0t2LiHBhcx1za7r50ysW8cVL53PJwgYuX1zPhxbPitQ1FjZWs6ixcsxJcnGcwpKIRE6m0X3XMLYNL1EaB8ar+XgBHK8QGB0wkim62B786OUV4uV/owNI9PQ+1T78idI2yWYE8QJ4WBXe+Zlz+d6z2xE5O/kJJcaGyJgi7GSCdqonlXRUerQt+gSfim/hk/KtV53Nt9dtiaSmJn+CSDw2NcnPPFnaJt7xDS+toV3RNt84eHIYBQQdB8FdEiDkOJF+9ZDjjLPB6avdfD6JLPur/D4YDeH3+yK2RHafDwi5zyUxJ4zI3XD8Akri/r/y+5DREOLzRWoKcLq4vGT2DJbMngHAjZc1R04In7m4iZaGChY2ViPAwoYKmhtm8JGWWha91IkAf37lIj7vqfkrljbykaWzaWmooGXWjEhAb5k1IxLQrzzb7ZqNV+ROFYJ4uXqvlU4SX0gT/QOMvnAqOoAHnbGBJX6HTmJFl6jdMpF6iw4ggsRNBaWq9jIyI/Bs0Wow2UU98YJWvPRCqrODqebwU1XpidJoiTpz4rVWrjhrpru8tlfHSPUEkcrnO9nPPJyWu/3aZWNO0Ak/85hZRrZhTMD/g9Z5/Pa941xz/hkIsPlAH+fPrfZ6id0vYqwtkX0y7x0NjvKvv97Hf/9gE6WlpWPG5PP5uNxbc0M5TsJ9dRwc4IolE9cd4p0QooNzOGB/86pzIraPL5sdeZ2saymTmF9XPk61hddNuXlFCz96uXNcB0Zs/3W0mh/zgxDcC7l+tsH9YV1zupsmnqJLtkjaZNXbZBV+MgWYaiE3bItWg7G54DF9+l7Hzp2fOTeSH04lvZBodjDVHH6qKj1RGu17z27n1k8t5e9i6g/x/I/tCEp28gufIFL9vG5PUNDN1HcEr6D+ky+1sigHN0AxJuD/cv0Bnth0kNd3HcMvqad04tkn895TI0H6hkNs6+6jomzs4Swt8dF23TLueLyDYMhJuK+TQ0EuaKrhY1HBWWd0nRgcl5aJvgQ+rP4jaRml4nZ0JFvVUMXrpkmg6Cat3hIU9ybbopesMJuKAoy2RffkxwuasR070TOC269ZNqU8fK5sidJoCXPucWyR3vvaiqmp80kWdCfzmSdqz42eEfz4S8vpPDIwLn2aaRgT8C9dUEvN9BJ8KEZGxxddJyrETvW94aLt8EjIXeckCgoVme4GQ07ifXkFXVMQLlR/5UNj8/aJLhRqqimjxDuvJlLzifKzE6UBJqvoJhtoJqOYowuzk50lRBdy4wXNVPv0J9tLn6jtcVqJj+hrFdLhm2yXTDL/p6rOp9SCmeJnfusEM4jmOndZ75UPrGdNlvP4xgR8n/iYVuIb015pkWPI6XVQwu2X8S6fj7eYWbJiWneS/Oy4NECS/GwqQW+ygSYTRdBk6ZhUZybx8sPx+vTTDWDRPebJZleTsiWZoU3V/4nSQqkWmdOdwUSn4uJdjPW9Z7fzzzdegoq6H3I2YUzAt8g/FjVWjuldj27Ni9dVE++CqUnnZ5NcPBVtS6WoN5VAM9kgPNl0TKozk3i2RKt0ppuCCReBP3BmXcJZWyrpm6l2yWRjtjaVFsyUP/OYE+WKs2bifHIJoWCQb35iMU4oxN7jp3LxE7UB3yJz2H1kYEybWap90mOWP0g3P5sgFxvvpJEs+KQaaCYbhFJNx4SPY6p54mycoBLZwkXgqdRI0u6SmWR9ZTKztUwU8ePZwqm4QysW8cBv3mdubQXfenwzh04OR+p54RsDZRvmBHyvoyBZ+2OqrZb5asvUHZEgncJ65Iny8OnmZxMpungnjWy36GVCFaZTG8jECSpR3jl2dpWr/U+2vpJuWigTn2/4+f5fvcdNH1pAU00Zd153Ln/7yGbc+0Z7S7X7fN66XdmDMQF/UWPlhO2PurdlFjrCKZ2mmrK465Gn0o6X9EcUZ0GzVAt98U4amQh+ydosM6EKM32BVGy9IN77M1FDSFQ3SbdonOjzzcRsLVsnyPB35K+vcr+v/+J3u67+/rPnjcnZ+3y+rK+nY0zA331kYML2R9uWmV3sPnqKu5/Zxjc/uSSl/Gy8PHyyXui9UfeAzUShLxPqbaptlrmoDaRSL5joIq90Al3cukmaRd5MXTCXj+/IqmvOocQH1104lzse30IwNLZAW1ri48z6yqx26fgDgUDWyCeLNWvWBFauXDml/z3WP8wv27sYDTmMBkMMhaBE3C9NyFFxbYnsk3nv4HDIXXY4pPD7JfLekLdy5JVLG3mu41DScQ2HFNeeP5tFs7J7WXW2UVdRyry6Cr69bguv7zzGmQ1V3P3MNppmVvKdJzsS2+or+c6TW2meWcEdj3fwyo4jLGyoZLVnW/3UNi5sqo1cnHXt+bNZ0FDN5YvrWdhYxZVLZ3HFkoaIrXlmJSsWu8tENM2sZMVZM+Pakm1Lx5ZrjkTHIa6toYrLlzRy7QWzWdDo2mLfX1NRyq0PbUr4eaRq+98v7hz3eSUdW4o23T7LiG1JA2Wlfu7/1R6GRkJjYkU4Xnzu4nnUV5VN8Esbj7a2tu5AILBmovcZo/At8g8R4YolDThejjdVVThhe98ki7f5zL9nOu+b6bGlsq+JLvJK52KoQj7mueD96eu7GRoNkd3my8QwRuEf7x/m4Q373duHOYoRB/cG4JFc+XhbIvtk3jsadAjh3mBkWunYwqvfJ3zs7Fm8uO3IhOO67oI52it8cJeb/trajZNShddfOJeLz5yZUNHFqtOCVG8Z5p1XV8H82ulcd8EcLmyuy9jYUhrvFFR3PFtLYzV/eNl8qstKuKCplhW59KFAeaeV+NnQdQKUe9+JaPh9wucvabIKPyXYLp2CwFRU4URFW13Umw61gVweh2zl1XX7LKNtD/zmfb7zexegHIfb/nPTmN9O+Naq2YRk+8quyaC1tVW1t7dP6X8dx+GlbYcLskvnlR1HJ+zS+fMrF1FSov/5N/w5hAu0r+06zoqzZvLKjqMp2V7deSzSbbG/Z4gVi+snzZGqzfLqyau7D5cvrueV7Uf4u8c2c/OKFubUTAfcePHRs2eNySqkChFZr5Rqneh9+kcYD7uPnnLbwILujUeODQxTX1Xm9rYCgyNB26WTA4Q/B6sKLW8+6xAFzYu7OmzQUfz4tfcii6VNK/FxZkNVVu96ZYzCV0qx63A/eAt3Rd/d6v1jA/zlf7zDySGH6mlCScyCaH2jaow9ni2RfWAoyIiCaUBl+fiAv/r6c/nbR7dEFk+Lt6/+UcWaP7qYT5w3d0q+FxLSVfjZUmp25mAOrwk+vLrzGE01ZXT1DNHsxSkRGXMXucmg6BS+iIxZ533JnJrojW5+LHIPptRgc/iTR7oK3+aTLW8h7Svb38d/eG4H/3xj9u9lG4YxCj8ZlFLsPNTHLmHdJQAAGr5JREFUm+8d40cv7yQYOr1tNOTQMxSittwfSf+E17X4+8+ez5kzK9h7/BSO49B9cgjAu6OTG7BDoRBPbDrEdeefgd/vR0RonlkROWO3NFbSeWQg7swjvK99PUNcsbRxSrm7QkOhKnyrNs3hNcGHcC7/vWODGVkDv+gUfjKICEtmz2DxGdV8sKUBok5yiW5uHp5edR4ZYPXT2zg1HBpXFwAi9YL2vT2U+twlmtfc1DrmjJ1w5pHEpisKVeFbtWkOrwk+3P3MNppn5k7ZRxC+WcVUH8B84GVgK7AF+JpnDwD7gXe8x9UTcS1fvlwVGkKhkHp+S7d6ZmOX+oenO9SzG7vUc5v2Rx6x9he2dKtQKJTvYecN4eP1wpZuNTo6ql7cekiNjo6mbJvs+/PNa4IPuvGa4MOLWw9lNE4A7SqFeJ0JhR8E/lIptUFEqoH1IvK8t+37Sql/yMA+8oawYp1I4Yftuai0FzKswtfPB914TfBBW4Uf+wAeAz6Bq/C/OZn/1V3hP/vuPnX/G53q+c0HUjp7O46jdh7qU6FQKPK842Cv2nmwd2q27pPqhY6DEQURDAYzwxvHlmhbMBgsWEVl1aYZvCb4oLPCj0BEFgAXA28CHwZuEZGbgHbcWcCJOP+zElgJ0NzcnMnhZASTUfgoONo/jIiw5qZWrjx7VlLuziMD3PLzDdzqnfHjdZFMxvatRzdzqG+IP7viLO55tZNvXete2Zcubzxbsm2Fqqis2jSD1wQftFf4QBWwHvic9/cZuEvM+IDVwH0Tceiu8J/Z2KW+++QmddtDG9TTG/eq5zbtV89vPqC2H+hRO7pPjlPJ27tPqhejFHm6quG5TfvVPz67VQ0ODqrvP7ddDQ8Pa6l8dOI1wQfdeE3wQWuFLyKlwMPAz5RS/+mdSA5Fbf8JsC4T+8o1Jqvwj/QN4wAPvnUAcJfGaKg+nd9PppzTVQ1t67ZxqG8In/i459VOGqvLeOA372unfHTiNcEH3XhN8CFfCj/tPnxxexnvB44rpb4eZZ+jlOr2Xn8DuEwpdUMyrmz14acD5V3BqxwnbvumimrrBHj/2ADdJ4eYXVPmfshx+vJ3Hz1FS0MFnUcGEMic7XA/e08Mcrl3FWm4zzcb+8qWD7rxmuCDbrwm+LD76KmM9N+Hkcs+/A8DXwI2icg7nu1vgS+KyEW4CyjuAf5HBvZlYWFhYTFFpL0efiAQ2BsIBNoCgcA9gUDgXu+xMxAIPOLZ7gkEAr8IBAL9E3Glsx5+ttB5ZICVD7Tzy/Z9/HJ9F890HOKJjd089s4BHnvnAA9t2B+xP/L2Ada27+OlHUd58t2DrHv3IOs2dfOs9z9PbT5Ic30lqx7ZxLy6Cv7qoY2se7c7Y7avP/gOa9u7ACHwxBbqq8q46+ltWdlXtnzQjdcEH3TjNcGHVY9s4ncWNTCzclpG4lSq6+FnvC0znYcJRdu71r2rbr7/1+qudRvVd9e9q+5+avOYi7KyWQiyRVtbQCwGXhN80LpoazKmWrT1cQwYX7TNZuHKFm1tAbEYeE3wQfu2zEw8rMK3Cl83XhN80I3XBB+swi9QWIVvjqKyatMMXhN8sArfKnyr8DXkNcEH3XhN8MEq/AKFVfjmKCqrNs3gNcEHq/CtwrcKX0NeE3zQjdcEH6zCL1BYhW+OorJq0wxeE3ywCt8qfKvwNeQ1wQfdeE3wwSr8AoVV+OYoKqs2zeA1wQer8K3CtwpfQ14TfNCN1wQfrMIvUFiFb46ismrTDF4TfLAK3yp8q/A15DXBB914TfAhXwo/7dUyM4l0VstUyl23/sTACLUVpXQeGaC2vCQtW+eRAY73D/P1tRt5tuMQb3Qe443O47zReYwXth3hpe3u47WdR3lzTw/PbzvM67tOsPPIIOv39rJ+70nW7+3hv3a7//PKjqMsbKhk9ZNbaZ5ZwR2Pd/DKjiMZs616dAuv7TyKAn70ym4aqsv4wYu7srKvbPmgG68JPujGa4IPq5/cyocX1XN8YIQTAyPUVU5La238VFfL9E30Bl0QXsZ45QPreWXHUW75+QZe2XE0LdstP9/A3hODKOUwzS989YpFVE3zc/PlLfhE4Re45aOLmDWjjD/76CJK/TCzws+NlzZxRlUJX72ihbqKEqb5hZsvb6HEJyhAodxn5bhn3gzZHMfBcRTK8ayOytq+LK9BPujGa4IPKPaeGIzEnc4jAwljW0aRyjQgV490Ujrh1Isu0zmb0jGD1wQfdOM1wYfYbemmdyi2lE7n4X6+tnajNtM5m9Ixg9cEH3TjNcGH2G2XL2lkZlXZlGNnqikdY7p0lFKMBEOIuKmNkaDjPadjywRHLnlDhBxFMBjEcRShUEhDH3TjNcEH3XhN8GHsNlekZx/G5PAtLCwsLJLDBnwLCwuLIkHWUzoichXwfwA/8FOl1F3Z2I9SilDIAREc5TDqODjKScuWCY5c8oZCjtsB4LjPTkg/H3TjNcEH3XhN8CF2W65SOlkN+CLiB34IfALYB7wlIo8rpToyva+uE4OcGAwiwOb9fRzrH2bz/r60bJngyCVv74gDwJt7elBA+94e7XzQjdcEH3TjNcGH2G1dJwZZMqcmQ9EwMbKt8D8A7FJK7QYQkQeB64GMB/ym2ulUT/OBCMvmVFJTXsqyOZVp2TLBkUvecj8MhmB5UxVv7+/n4nnVvLmnRysfdOM1wQfdeE3wIXZbU+30TIfEuMh2wJ8HdEX9vQ+4LBs76joxSJ+ncDcf6KPn1CibD/SlZcsERy55B0PusVi/rxeADV0ntfNBN14TfNCN1wQfYrd1nRhk6dxaso1sB/x41wqPSVaJyEpgJUBzc3NaO1OxL1R6tkxw5JwXCAbdL5ETzgvq5oNGvCb4oBuvCT6M25YjZDvg7wPmR/3dBByIfoNSag2wBqC1tXXK/s+vK2eGNz06b141tRWlnDevOi1bJjhyyRtO6Vy2oJZNhw7S2lxLe1efVj7oxmuCD7rxmuBD7Lb5deVTDX2TQrYD/lvAYhFZCOwHbgBuzMaO9vUM0TfiIEBH9wAnB0fp6B5Iy5YJjlzynk7p9APw9v4+7XzQjdcEH3TjNcGH2G37eoZYOjcTkTA5shrwlVJBEbkFeBa3LfM+pdSWbOxrfl05NeUliMB5c6upryzjvLnVadkywZFL3qpSoX9U0do8g7f397O8uYb1e3u18kE3XhN80I3XBB9it5mi8FFKPQU8le39iAilfh8i4PP5KPX7Is/p2DLBkTPeEj8yGqSkpAQf4Pf79fNBQ14TfNCN1wQforelszTyZGCvtLWwsLAoEtiAb2FhYVEksAHfwsLCokhgA76FhYVFkcCY9fARodQviAji81Fa4j2nY8sERw55S3yCAH6fH58PfKKfD9rxmuCDbrwm+BCzjRwVbfN+W8Poh73Fob3FoW68JvigG68JPuTrFofGKPzdR0+x+qmtCIL4lnH3M9vgqrPTsmWCI5e8beu2cahvCJ/4uOfVThqry3jgN+9r5YNuvCb4oBuvCT7EbjuzoYqzZlVlP1CmclbI1cMqfKvwdeM1wQfdeE3wIV8K397EvABuYGxvYq4vrwk+6MZrgg/2JuZpQgFKOYC4r1Hp2zLBkUNex3FwHIVyPKujtPNBO14TfNCN1wQf4mzLCVKZBuTqYVM6NqWjG68JPujGa4IP+UrpGKPwbdHWFm1tAbE4eE3wwRZtrcK3Cl9DXhN80I3XBB+swk8TVuFbhW/VZnHwmuCDVfhpKnzHcdSOg71q58FeFQqF1M5DfSoUCqVlywRHLnl3dJ9UL3QcjCiIYDConQ+68Zrgg268JvgQu81xnCnHPqVSV/h2LR0LCwuLIoE5ffhHBlj5QDvr3u2mub6SVY9sYl5dBX/10MYp2zLBkUverz/4DmvbuwAh8MQW6qvKuOvpbVr5oBuvCT7oxmuCD7HbPrK4kZmV06YcO1Ptw897Gif6YYu2tmirG68JPujGa4IPtmibJmzR1hZtbQGxOHhN8MEWba3CtwpfQ14TfNCN1wQfrMJPE1bhW4Vv1WZx8Jrgg5YKH/gesA14F3gEqPXsC4BB4B3vcW8qfFbhW4WvG68JPujGa4IPuir854HblFJBEfkucBvw1962TqXURWnypwyr8K3Ct2qzOHhN8EFLhR/9AD4L/EydVvibJ8thFb5V+LrxmuCDbrwm+KCrwo/GV4C1UX8vFJG3gV7gdqXU6/H+SURWAisBmpubp7xzq/Ctwrdqszh4TfChYBU+8AKwOc7j+qj3rMLN4Yv3dxlQ771eDnQBMybal11awS6toBuvCT7oxmuCD/laWmFCha+U+niy7SLyZeBa4GPejlFKDQPD3uv1ItIJLAHaJ3U2srCwsLDIGNJaWkFErgK+DXxKKXUyyt7Y1tY2HAgElIi0AN8E/j4QCAwm47NLK9ilFXTjNcEH3XhN8EHLpRWAXbjpmjHtl8DngS3ARmADcF0qfLZoa4u2uvGa4INuvCb4oGXRVil1VgL7w8DD6XBPFrZoa4u2toBYHLwm+FCwRdtcPqzCtwpfN14TfNCN1wQftFT4hQSr8K3Ct2qzOHhN8MEqfKvwrcLXkNcEH3TjNcEHq/DThFX4VuFbtVkcvCb4YBW+VfhW4WvIa4IPuvGa4INV+GnCKnyr8K3aLA5eE3ywCt8qfKvwNeQ1wQfdeE3wwSr8NGEVvlX4Vm0WB68JPliFbxW+Vfga8prgg268JviQL4Wf1lo6mUZaa+kc7udrazfyyo4jLGyoZPWTW2meWcEdj3dM2ZYJjlzyrnp0C6/tPIoCfvTKbhqqy/jBi7u08kE3XhN80I3XBB9it12+pJGZVWVTjp2prqVjTEpHAUo5gLivUenbMsGRQ17HcXAchXI8q6O080E7XhN80I3XBB/ibMsJUpkG5OphUzo2paMbrwk+6MZrgg/5SukYo/Bt0dYWbW0BsTh4TfDBFm2twrcKX0NeE3zQjdcEH6zCTxNW4VuFb9VmcfCa4INV+FbhW4WvIa8JPujGa4IP+VL4ti2zAFqybFumvrwm+KAbrwk+2LbMNKF7S1YmbLYtMw+8JvigG68JPsTZlhOkMg3I1cOmdGxKRzdeE3zQjdcEH/KV0jFG4duirS3a2gJicfCa4IOWRVsgAOwH3vEeV0dtuw3YBWwHPpUKn1X4VuHrxmuCD7rxmuBDvhR+WkXbtra2K4DXlFI3BQKBewOBwE4AEVnmnQwuBh4D1ra1tf0wEAgkTVXZoq0t2urGa4IPuvGa4INpRdvrgQeVUsPAeyKyC/gA8Oss7U/7gk0mbLZomwdeE3zQjdcEH+JsywlSmQYkeuCq+D3Au8B9QJ1n/2fgj6Le9y/AFybisykdm9LRjdcEH3TjNcGHgk3piMgLbW1tt7a1tX015rEfeBS4C7gXuBT4fCAQeKytre0aYF8gEHgXoK2t7XpgZyAQ2BqHf2VbW9uP29raVpaWls79xje+MaUTl03p2JSOTS8UB68JPhRsSkcp9fFUdigiPwHWeX/uA+ZHbW4CDiTgXwOsAWhtbZ3yzEb36VwmbDalY9MLRcFrgg9xtuUEqUwDEj2AOVGvv4Gbtwc4F9gIlAELgd2AfyI+m9KxKR3deE3wQTdeE3wo2JROMrS1td3T1tZ2Z1tb258CdcAtgUCgPxAIHGlra6sHfgrcCPyFUmrHRHy2S8emdHTjNcEH3XhN8KFgUzrJoJT6UpJtq4HV6fBPaizoPZ3LhM2mdGx6oSh4TfAhzracIJVpQK4eNqVjUzq68Zrgg268JvigZUon07ApHZvS0Y3XBB904zXBBy1TOoUEpRQjwRAibmpjJOh4z+nYMsGRS94QIUcRDAZxHEUoFNLQB914TfBBN14TfBi7zRXp2YcvJ3uxsLCwsMg7jFH4iFDqF0QE8fkoLfGe07FlgiOHvCU+QQC/z4/PBz7RzwfteE3wQTdeE3yI2YZIbuJkKon+XD1s0dYWbXXjNcEH3XhN8MEWbbFFW1u01Y/XBB904zXBB1u0TRPKFm1t0dYWEIuE1wQfbNHWwsLCwiKLsAHfwsLCokhgA76FhYVFkcAGfAsLC4sigQ34FhYWFkUCG/AtLCwsigQ24FtYWFgUCWzAt7CwsCgS2IBvYWFhUSQw6krbUMgBERzlMOo4OMpJy5YJjlzyhkKOexcdx312Qvr5oBuvCT7oxmuCD7HbcnWlrTEBv+vEICcGgwiweX8fx/qH2by/Ly1bJjhyyds74gDw5p4eFNC+t0c7H3TjNcEH3XhN8CF2W9eJQZbMqclQNEwMYwL+/LpyaspLEIHz5lZTX1nGeXOr07JlgiOXvFWlQv+oorV5Bm/v72d5cw3r9/Zq5YNuvCb4oBuvCT7EbptfV56TOJlWwBeRtcBS789aoEcpdZGILAC2Atu9bb9RSt2czr5SgVuQEBAij3RsmeDIJa9PBFD4fT7PJFnbl+U1xwfdeE3wYdy2HCGtgK+U+oPwaxH5R+Bk1OZOpdRF6fBPBjalY1M6Nr1QHLwm+KB1SkdEBPh94MpM8E0FTbXTqZ7mAxGWzamkpryUZXMq07JlgiOXvOV+GAzB8qYq3t7fz8XzqnlzT49WPujGa4IPuvGa4EPstqba6TmJk5nK4f8ucEgptTPKtlBE3gZ6gduVUq9naF9xsa9niL4RBwE6ugc4OThKR/dAWrZMcOSSdzDkHov1+/oBeHt/n3Y+6MZrgg+68ZrgQ+y2fT1DLJ2biUiYHBMGfBF5AZgdZ9MqpdRj3usvAr+I2tYNNCuljonIcuBRETlXKdUbh38lsBKgubl5suOPwCp8q/Ct2iwOXhN8KFiFr5T6eLLtIlICfA5YHvU/w8Cw93q9iHQCS4D2OPxrgDUAra2tU25G7ToxSJ+Xw958oI+eU6NsPtCXli0THLnkPa3w3fPqhq6T2vmgG68JPujGa4IPsdu6TgyydG4t2UYmUjofB7YppfaFDSLSCBxXSoVEpAVYDOzOwL4SQilF+GzhOA5Kec/p2DLBkUPeMILBIAChUEg7H7TjNcEH3XhN8CFmm04XXt3A2HQOwArgThEJAiHgZqXU8QzsKyHcurELn/gQcZ/TsWWCI5e8p+E1e/n080E3XhN80I3XBB9it0XHr2zCN/FbkkMp9cdKqXtjbA8rpc5VSl2olLpEKfVEuvuZCM31lTRUltBYVcqFzXXMrpnOhc11adkywZFL3poy9+P88FmN+IAPtjRo54NuvCb4oBuvCT7Ebmuur8x2iHShlCqYx/Lly9VU4TiO2nGwV+082KtCoZDaeahPhUKhtGyZ4Mgl747uk+qFjoNqdHRUvbj1kAoGg9r5oBuvCT7oxmuCD7HbHMeZcuxTSimgPZUYKypHuaNU0Nraqtrbx9V1LSwsLCySQETWK6VaJ3pf2ikdCwsLCws9YAO+hYWFRZHABnwLCwuLIoEN+BYWFhZFAhvwLSwsLIoENuBbWFhYFAlswLewsLAoEhRUH76IHAHeT4OiATiaoeFkE3acmYUu4wR9xmrHmVlke5xnKqUaJ3pTQQX8dCEi7alcfJBv2HFmFrqME/QZqx1nZlEo47QpHQsLC4sigQ34FhYWFkUC0wL+mnwPIEXYcWYWuowT9BmrHWdmURDjNCqHb2FhYWGRGKYpfAsLCwuLBDAi4IvIVSKyXUR2icjf5Hs80RCRPSKySUTeEZF2zzZTRJ4XkZ3ec12exnafiBwWkc1RtrhjExc/8I7xuyJySZ7HGRCR/d5xfUdEro7adps3zu0i8qkcjnO+iLwsIltFZIuIfM2zF9QxTTLOQjym00XktyKy0Rtrm2dfKCJvesd0rYhM8+xl3t+7vO0L8jzOfxOR96KO6UWePT+/p1QWzS/kB+AHOoEWYBqwEViW73FFjW8P0BBjuxv4G+/13wDfzdPYVgCXAJsnGhtwNfA0IMAHgTfzPM4A8M04713mfQfKgIXed8Ofo3HOAS7xXlcDO7zxFNQxTTLOQjymAlR5r0uBN71j9UvgBs9+L/Cn3uuvAvd6r28A1uZ5nP8GfCHO+/Py2Zug8D8A7FJK7VZKjQAPAtfneUwT4Xrgfu/1/cDv5WMQSqnXgNh7DSca2/XAvysXvwFqRWROHseZCNcDDyqlhpVS7wG7cL8jWYdSqlsptcF73QdsBeZRYMc0yTgTIZ/HVCml+r0/S72HAq4EHvLsscc0fKwfAj4mObhhbJJxJkJePnsTAv48oCvq730k//LmGgp4TkTWi8hKz3aGUqob3B8fMCtvoxuPRGMrxON8izcdvi8qLVYQ4/RSCRfjKr2CPaYx44QCPKYi4heRd4DDwPO4M4wepVQwzngiY/W2nwTq8zFOpVT4mK72jun3RaQsdpwecnJMTQj48c7ehdR69GGl1CXAp4E/E5EV+R7QFFFox/keYBFwEdAN/KNnz/s4RaQKeBj4ulKqN9lb49hyNtY44yzIY6qUCimlLgKacGcW5yQZT97GGjtOETkPuA04G7gUmAn8dT7HaULA3wfMj/q7CTiQp7GMg1LqgPd8GHgE9wt7KDx9854P52+E45BobAV1nJVSh7wfmAP8hNMphryOU0RKcYPoz5RS/+mZC+6YxhtnoR7TMJRSPcAruDnvWhEpiTOeyFi97TWkng7M9Div8tJnSik1DPwreT6mJgT8t4DFXtV+Gm6h5vE8jwkAEakUkerwa+CTwGbc8X3Ze9uXgcfyM8K4SDS2x4GbvO6CDwInw2mKfCAm3/lZ3OMK7jhv8Lo1FgKLgd/maEwC/AuwVSn1T1GbCuqYJhpngR7TRhGp9V6XAx/HrTm8DHzBe1vsMQ0f6y8ALymvSpqHcW6LOtELbp0h+pjm/veUi8pwth+4Fe8duLm9VfkeT9S4WnC7GzYCW8Jjw80pvgjs9J5n5ml8v8Cduo/iKo4/STQ23CnoD71jvAlozfM4H/DG8S7uj2dO1PtXeePcDnw6h+P8CO60/F3gHe9xdaEd0yTjLMRjegHwtjemzcC3PHsL7klnF/AfQJlnn+79vcvb3pLncb7kHdPNwP/jdCdPXj57e6WthYWFRZHAhJSOhYWFhUUKsAHfwsLCokhgA76FhYVFkcAGfAsLC4sigQ34FhYWFkUCG/AtLCwsigQ24FtYWFgUCWzAt7CwsCgS/H82csBbanry/AAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# multiple titles are stacked into a single 2D array\n",
"# just passinsg this will crash ESMPy \n",
"plt.scatter(ds['lon'], ds['lat'], s=0.2) "
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.PathCollection at 0xd1c403e10>"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAGxtJREFUeJztnXu0XFV9xz+/m4RQQIGQIIheEzBq4xsjQmsrPloeLk2ttI3W+u5dKqi4io+srFX31maVImpFQEiBIhpFxSIRCI9EqaVLHjeahJCQ5AIigbxoBRJZvXBnfv1jzoUhmXvn3DtzZ86+5/tZ66wzs8/JZ35zkuzZ891nzjF3RwghxOSnp9sFCCGE6Azq8IUQoiSowxdCiJKgDl8IIUqCOnwhhCgJ6vCFEKIkqMMXQoiSoA5fCCFKgjp8IYQoCVO7XUA9M2fO9NmzZ3e7DCGESIrVq1c/4u6zmu1XqA5/9uzZ9Pf3d7sMIYRICjN7IM9+inSEEKIkqMMXQoiSoA5fCCFKgjp8IYQoCerwhRCiJKjDF0KIkqAOXwghSsKk6PDdnYGde6hWq/ust+zYzcCO3eNqS9FRqVRYuWE7qzZsp1KpJP1e5JCjTI6BnXuY6FvOTgkhTOgLjIWlS5eGvr6+Mf+5e3f9njO+9yuOOvQAFl9917PWn71qLdeu20bvYQeOuS1Fx5SeHs760VquXbedWc/Zn7NX3JPse5FDjjI5Fl99FyccM5MZB+435j4wxrgthLC06Y7uXpjlda97nY+HSqXiqzbu8Keeemqf9c13b/OVd28bV1uKjsHBQf/Kig1+7ooNPjg4mPR7kUOOMjlWbdzhlUplXH0g0J+nj50cI/yde/jsVevonXEAS67b+Kz1F5dv4JbNu5gz88Axt6XomDa1h3Nv2sydD/yOWc/dn/NWDST7XuSQo0yOJddt5I0vnsmMg6aPuQ/MO8Iv1LV0xosDjjdeexWw8bWl6Kg6lSwG9Kqn/V7kkKNMjmw9obQjimnXokhHkY4ccpTVMSkiHTM7OcZ4XYzxzBjjASGEW0faV5GOIh055CirI/lIx8ymABcAfwZsBe40s+XuvqGdr+PuPDlUpVqtNlhXMGOcbek5KpXK05FOpVJJ+r3IIUe5HFV8gk/LnOjz8I8DBtz9Pnd/ErgSWDDBrymEEKIBEz1pexTwYN3zrcAb6ncwsz6gD6C3t3d8r2LGtKmG9fTsu55imI2zLUFHjz3zGd5jab8XOeQolWOqgdn4+sC8jHVidSwL8FfAJXXP/w745kj7a9JWk7ZyyFFWR/KTtjHGg4EFIYTvZs/fDuweaeJWk7aatJVDjrI6kp+0Be4E5prZHOAhYCHw3na/iGvSVpO2csiRvGPiJ20ntMN39yEzOwO4EZgCXObud0/kawohhGjMhP/S1t2vB66f0BfRpK0mbeWQI3VH6pO2Y100aatJWznkKKsj+UnbsaJJW03ayiFHWR2TYdK2IzS8YFGqF1DSxdPkkKOcjmw9obQjimnXokhHkY4ccpTV0YlIZ1KM8O975AnOueEeOPll+6yXXL8Rw7CeeWNuS9Gx7fgXcdEv7gOMIw45gO/c9kCy70UOOcrkOOeGe+idcSwvPvygiess2zEyb9eiEb5G+HLIUVaHRvg50QhfI3w55EjdoRF+TjTC1whfDjlSd2iEnxON8DXCl0OO1B0a4eekWq36lh27vVKp7LPevP1x37L98XG1pegYGhp6evQwNDSU9HuRQ44yObbs2O3VanVcfSA5R/gTfQMUIYQQBWFy/NJ21+8543u/4qhDD2Dx1Xc9a/3Zq9Zy7bpt9B524JjbUnRM6enhrB+t5dp125n1nP05e8U9yb4XOeQok2Px1XdxwjEzmXHgfmPuA/P+0rbrMU79okhHkY4ccpTVoUhHCCFE21CkU7CvdYp05JCjnI7CRzrU7ll7N1AF5u+1bREwAGwCTsrjU6SjSEcOOcrq6ESk0+p5+OuBvwQurm80s3nUbmf4cuD5wEoze4m7V1p8PSGEEOOkpUgnhPBICOF/YowfBG4KITwMEGP8KLDe3W8JITya3bx8Swhh62g+RTqKdOSQo6yOwkc6wwtwC3WRDnA+8L6655cCpzXzKNJRpCOHHGV1FCLSMbOVwBENNi1292tG+mONPltG8PcBfQC9vb3NyhFCCDFOmkY6IYQrQggXNlg2De/TINJ5NXBwCOHW7PnpwE8bRTohhNUhhKUhhKXLli1TpKNIRw45SulIOdJ5ObAWmA7MAe4DpjTzKNJRpCOHHGV1dCLSsayDHhdm9i7gm8As4FFgjbuflG1bDHwYGALOdPcVzXzz58/3/v7+cdcjhBBlxMxWu/v8Zvu19Etbd7/a3V/g7tPd/XnDnX22bYm7H+PuL83T2bdYBwM791CtVvdZb9mxm4Edu8fVlqKjUqmwcsN2Vm3YTqVSSfq9yCFHmRwDO/fQygA8D/qlbcFyPGX4cshRTkcyGX67FmX4yvDlkKOsjsJn+O1GGb4QQoydjmT4RcGV4SvDl0OOxB3K8HOiDF8ZvhxypO5Qhp8TZfjK8OWQI3WHMnwhhBBNKVWGL4QQojmTosN3Tdpq0lYOORJ3aNI2J5q01aStHHKk7tCkbU40aatJWznkSN2hSVshhBBNKdWkrSvDV4YvhxyJO5Th50QZvjJ8OeRI3aEMPyfK8JXhyyFH6o7CZ/hm9hXgHcCTwL3Ah9z90WzbIuAjQAX4lLvf2MynDF8IIcZOpzL8m4FXuPurgM3AouzF5wELqd3q8GTgQjOb0uJrCSGEaIFW73h1k7sPZU9vA16QPV4AXOnug+5+PzAAHNfKazWpQ5O2mrSVQ46kHUlN2sYYvwL8JISwLsa4ENgQQliXbXszsCOEsGE0hyZtNWkrhxxldRRi0hZYCaxvsCyo22cxcDU8PSdwAfC+uu2XAu8ewd8H9AP9vb2945qw0KStJm3lkCN1R+EnbQHM7APAx4C3uvsTWdui7MPkn7PnNwLB3X85mkuTtkIIMXY6MmlrZicDnwfeOdzZZywHFprZdDObA8wF7mjltUbDleErw5dDjsQdhc/wY4w3AAcAp8UYPxZjfE0I4doQwq4Y42HAJcB7qZ2WubmZTxm+Mnw55CiroxAZficX/fBKGb4ccpTVkUSG306U4QshxNgp1cXThBBCNGdSdPiuSVtN2sohR+KOwk/athtN2mrSVg45yurQpG1ONGmrSVs55EjdoUlbIYQQTSnVpK0rw1eGL4cciTuU4edEGb4yfDnkSN2hDD8nyvCV4cshR+oOZfhCCCGaUqoMXwghRHMmRYfvmrTVpK0cciTu0KRtTjRpq0lbOeRI3aFJ25xo0laTtnLIkbpDk7ZCCCGa0qk7Xn3ZzNaZ2Rozu8nMnp+1m5mdZ2YD2fZjW3mdZrgyfGX4csiRuKPwGX6McbW7/2sI4aIY46HAaSGEa2OMpwKnAMcDvwbODyFc0synDF8ZvhxylNWRVIYPLAK+lT2+GHhP3bZNwJHNHMrwleHLIUdZHUlk+Ga2BHg/8BjwZnffZWbXAme7+63ZPquAz7v7qAG9MnwhhBg7bcvwzWylma1vsCwAcPfF7v5CYBlwxvAfa6Bq+MliZn1m1m9m/bt27WpWTkNcGb4yfDnkSNxRiAw/hHBFCOHCBsum+v1ijA8C/xRCuDDG+CbgdyGE9dm2xcA3Qwh7GvhXhxCWhhCWLlu2TBm+Mnw55Cilo/AZPjC37vEngauyx28HVlAb6R8P3JHHpwxfGb4ccpTVUfgM38x+DLwUqAIPAB9z94fMzIDzgZOBJ4APeZP8HpThCyHEeOjIefju/m53f4W7v8rd3+HuD2Xt7u6nu/sx7v7KPJ19i3Uow1eGL4ccSTsKkeF3Ep2HrwxfDjnK6ih8ht/uRRm+Mnw55Ciro/AZfrtRhi+EEGOnVDdAcWX4yvDlkCNxhzL8nCjDV4YvhxypO5Th50QZvjJ8OeRI3aEMXwghRFOU4Sea4ynDl0OOcjqU4edEGb4yfDnkSN2hDD8nyvCV4cshR+oOZfhCCCGaogw/0RxPGb4ccpTToQw/J8rwleHLIUfqDmX4OVGGrwxfDjlSdyjDF0II0RRl+InmeMrw5ZCjnI5kMnwzOyvG+N8xxgtCCE+YmcUYz4sxfjPG+NEY4+0hhG3NPMrwleHLIUdZHUlk+MALgRup3eJwZtZ2Ks++p+3teVzK8JXhyyFHWR1JZPhmdhXwZeAaYL67P2JmFwO3uPv3s302ASe6+6ijfGX4QggxdjqS4ZvZO4GH3H3tXpuOAh6se741a2vk6DOzfjPr37Vr17jqcGX4yvDlkCNxRyEyfDNbGWP8XIzxE3stDwGLgQ+HEAZjjGcCS0MIT8QY/xa4NYTwW4AY4weAmxvl+CGE1SGEpSGEpcuWLVOGrwxfDjlK6Sh0hg+8EtgJ/CZbhoDfAkcAFwPvqdt3E3BkM6cyfGX4cshRVkcnMvxxRzrufpe7H+7us919NrXY5lh33w4sB95vNY4HHvMm+b0QQoiJpW2XVtgr0hkATgDOA04B+kIIDzdz6LRMRTpyyFFWR6EjnYlYFOko0pFDjrI6OhHp6NIKQgiROLq0QqKnZrXqqOi0TDnkSNJRiNMyO4kyfGX4cshRVocy/Jwow1eGL4ccqTuU4QshhGiKMvxEczxl+HLIUU6HMvycKMNXhi+HHKk7lOHnRBm+Mnw55Ejd0YkMf1JEOkIIIZqjSKdgX+sU6cghRzkdinRyokhHkY4ccqTuUKQjhBCibSjSKdjXOkU6cshRTocinZwo0lGkI4ccqTs6Eenol7ZCCJE4nbqJeTCzh8xsTbacWrdtkZkNmNkmMzuplddphrt+aTvcVtEvbeWQI0lH4X9pG2M8EfiFu78/hHBRCGELgJnNAwLwWuAa4AcxxgtCCKO+G2X4yvDlkKOsjsJn+NQ69bMatC8CFtU9vxE4oZlPGb4yfDnkKKujExl+O07LPMPM1pnZZWZ2aNZ2FPBg3T5bszYhhBBdommkY2YrY4yfizF+Yq/lIeAnwNnARcDrgXeHEK6JMb4d2BpCWAcQY1wAbAkhbGzg74sxXhxj7Js2bdrzP/OZz4z5TSjSUaQjhxypOwof6dQvwGxgvSvSUaQjhxxyTL5Ix8yOrHv6LmB99ng5sNDMppvZHGAucEcrryWEEKI1Wj1L51sxxi/FGD8OHAqcEULYE0LYFWM8DLgEeC/wKXff3Myns3QU6cghR1kdSUU67VgU6SjSkUOOsjoKH+kIIYRIB108rWBf6xTpyCFHOR2KdHKiSEeRjhxypO5QpCOEEKJtKNIp2Nc6RTpyyFFOhyKdnFQqFV+1cYc/9dRT+6yH443xtKXoGBwc9K+s2ODnrtjgg4ODSb8XOeQok2PVxh1eqVTG1QeSM9KZOuaPkgJy3yNPcM4N98DJL9tnveT6jRiG9cwbc1uKjm3Hv4iLfnEfYBxxyAF857YHkn0vcshRJsc5N9xD74xjefHhB01cZ9mOkXm7Fk3aatJWDjnK6tCkrRBCiLahSduCTdxo0lYOOcrp0KRtThTpKNKRQ47UHYp0hBBCtA1FOgX7WqdIRw45yulQpJMTRTqKdOSQI3WHIh0hhBBto+VIx8w+GWP8Tozx4zHG2SGEm7P2RTHGy2OMp8cYt4QQ7m3mUqSjSEcOOcrqKHykA7wZWAlMz54fnq3nAWuB6cAc4F5gSjOfLq2gSyvIIUdZHSlcWuHjwNnuPph9eOzM2hcAV2bt95vZAHAc8MsWX68hurSCLq0ghxypOwp/aQVgDRCB24H/BF6ftZ8PvK9uv0uB05r5NMLXCF8OOcrqKMQI38xWAkc02LQYmErt5uXHA68HfmhmRwPW6LNlBH8f0AfQ29vbrJyGaISvEb4ccqTuSGGEfwNwYt3ze4FZwCJgUV37jcAJzXw6LVOnZcohR1kdKZyW+RPgLQBm9hJgP+ARYDmw0Mymm9kcYC5wR4uvJYQQogVaOi0zxrgO+GCM8UvAu4BPuPv9IYRdMcbDgEuA9wKfcvfNzXw6LVOnZcohR1kdhT8ts92LJm01aSuHHGV1FGLSNgU0aatJWznkSN1R+Enbdi8a4WuEL4ccZXVohJ8TjfA1wpdDjtQdGuHnRKdl6rRMOeRI3ZHCaZlCCCESQTdAKdipWTotUw45yunQaZk50aStJm3lkCN1hyZtc6JJW03ayiFH6g5N2uZEI3yN8OWQI3WHRvg50QhfI3w55EjdoRF+TjTC1whfDjlSd3RihD85ztLZuYfPXrWO3hkHsOS6jc9af3H5Bm7ZvIs5Mw8cc1uKjmlTezj3ps3c+cDvmPXc/Tlv1UCy70UOOcrkWHLdRt744pnMOGj6mPvAvGfpTIpIxwHHG6+9Ctj42lJ0VJ1KdqsZr3ra70UOOcrkyNYTSjuimHYtinQU6cghR1kdnYh0JsUIX5O2mrSVQ47UHYWftAV+QO1G5muA3wBr6rYtAgaATcBJeXwa4WuEL4ccZXUUfoTv7n8z/NjMvgo8lj2eBywEXg48H1hpZi9x90orrzcSGuFrhC+HHKk7Cj/CH14AAx4E5vozo/uO3cRcI3yN8OWQI3VH4Uf4dfwJsMPdt2TPjwJuq9u+NWubEDTC1whfDjlSdxRihA+sBNY3WBbU7fMt4B/qnl8AvK/u+aXAu0fw9wH9QH9vb++4Pt00wtcIXw45Und0YoTf9IdXIYQrQggXNlg2AZjZVOAS4FMhhMcBYoyvBg4OIdyaPT8d+GkIYWsD/+oQwtIQwtJly5bph1f64ZUccpTSkcoPr94G3OPu9Z35cuB7ZvY1apO2c4E72vBaDWn4Y4dUf3yhH17JIUc5Hdl6QsnzNWC0Bbgc+FiD9sXAvdROyzwlj0uTtop05JCjrI5CRDrNCCH8JITQ36D9v0II54UQzg8hDORx6Vo6inTkkKOsjlQina7T8GtUql/rFOnIIUc5Hdl6QsnzNaBTiyIdRTpyyFFWRxKRTjtRpKNIRw45yupQpJMTd+fJoSrVarXBuoIZ42xLz1GpVJ6OdCqVStLvRQ45yuWoDp/wMmFMig4fM6ZNNaynZ9/1FMNsnG0JOnqs5+nD0mNpvxc55CiVY6qB2cT2la1k7u1elOErw5dDjrI6lOHnRBm+Mnw55EjdoQw/J64MXxm+HHIk71CGnw9l+Mrw5ZAjdYcy/Hwow1eGL4ccqTuU4edEGb4yfDnkSN2hDD8nrgxfGb4cciTvmPgMv6f5LkIIISYD6vCFEKIkqMMXQoiSMCkyfJ2W+UybTsuUQ45EHUU/LRN4DXAbsIbajciPy9oNOA8YANYBx+bxjfe0zGq16lt27PZKpbLPevP2x33L9sfH1ZaiY2ho6OlTvoaGhpJ+L3LIUSbHlh27vVqtjqsPJOdpmeYtzAqb2U3A1919hZmdCnzO3U/MHn8SOBV4A/ANd39DM9/8+fO9v3+fm2cJIYQYBTNb7e7zm+3XaobvwHOzxwcDD2ePFwBXZB8+twGHmNmRLb6WEEKIFmg1wz8TuNHMzqX24fFHWftRwIN1+23N2rbtLTCzPqAPoLe3t8VyhBBCjETTDt/MVgJHNNi0GHgr8Bl3/7GZ/TVwKfA2ahn+3jTMjtx9KbAUapFOzrqFEEKMkaYdvru/baRtZnYF8Ons6Y+AS7LHW4EX1u36Ap6Je4QQQnSBVjP8h4E3ZY/fAmzJHi8H3m81jgcec/d94hwhhBCdo9UM/++Bb5jZVOD/yLJ44HpqZ+gMAE8AH2rxdYQQQrRISx2+u98KvK5BuwOnt+IWQgjRXlo6D7/dmNku4IFx/vGZwCNtLGciSaXWVOqEdGpNpU5Ip9ZU6oSJq/VF7j6r2U6F6vBbwcz68/zwoAikUmsqdUI6taZSJ6RTayp1Qvdr1cXThBCiJKjDF0KIkjCZOvymt/cqEKnUmkqdkE6tqdQJ6dSaSp3Q5VonTYYvhBBidCbTCF8IIcQoTIoO38xONrNNZjZgZl/odj31mNlvzOwuM1tjZv1Z2wwzu9nMtmTrQ7tU22VmttPM1te1Nawt+9X0edkxXmdmxxag1mBmD2XHdk12We7hbYuyWjeZ2UkdrPOFZvZzM9toZneb2aez9kId11HqLOIx3d/M7jCztVmtMWufY2a3Z8f0B2a2X9Y+PXs+kG2f3eU6Lzez++uO6Wuy9s7/3ee5aH6RF2AKcC9wNLAfsBaY1+266ur7DTBzr7ZzgC9kj78A/EuXavtT4FhgfbPaqP1yegW1C+MdD9xegFoDcFaDfedl/w6mA3Oyfx9TOlTnkWQ3/AGeA2zO6inUcR2lziIeUwMOyh5PA27PjtUPgYVZ+0XAx7PHnwAuyh4vBH7Q5TovB05rsH/H/+4nwwj/OGDA3e9z9yeBK6ldj7/ILAC+nT3+NvAX3SjC3X8B/O9ezSPV1tV7HIxQ60gsAK5090F3v5/aJT6Om7Di6nD3be7+q+zxbmAjtUuDF+q4jlLnSHTzmLq778meTssWp3b9rquy9r2P6fCxvgp4q9lE3ztw1DpHouN/95Ohwx/p2vtFwYGbzGy11a79D/A8zy4ml60P71p1+zJSbUU9zmdkX4cvq4vGClFrFiW8ltpIr7DHda86oYDH1MymmNkaYCdwM7VvGI+6+1CDep6uNdv+GHBYN+p09+FjuiQ7pl83s+l715kx4cd0MnT4ua+93yX+2N2PBU4BTjezP+12QeOkiMf5W8Ax1O6tvA34atbe9VrN7CDgx8CZ7v74aLs2aOtYrQ3qLOQxdfeKu7+G2qXWjwP+cJR6ulbr3nWa2SuARcDLgNcDM4DPd6vOydDhF/ra++7+cLbeCVxN7R/rjuGvbtl6Z/cq3IeRaivccXb3Hdl/sCrwbzwTMXS1VjObRq0TXebu/5E1F+64NqqzqMd0GHd/FLiFWuZ9iNWu1Lt3PU/Xmm0/mPxxYLvrPDmLz9zdB4F/p4vHdDJ0+HcCc7MZ+/2oTdIs73JNAJjZgWb2nOHHwJ8D66nV94Fstw8A13SnwoaMVFvh7nGwV975LmrHFmq1LszO1pgDzAXu6FBNRu3Obxvd/Wt1mwp1XEeqs6DHdJaZHZI9/gNqd9XbCPwcOC3bbe9jOnysTwN+5tksaRfqvKfug96ozTPUH9PO/t1P9KxwJxZqs92bqeV6i7tdT11dR1M7s2EtcPdwbdTyxFXUbhizCpjRpfq+T+1r+1PURhsfGak2al8/L8iO8V3A/ALU+p2slnXU/vMcWbf/4qzWTcApHazzjdS+lq8D1mTLqUU7rqPUWcRj+irg11lN64F/zNqPpvahM0DtjnvTs/b9s+cD2faju1znz7Jjuh74Ls+cydPxv3v90lYIIUrCZIh0hBBC5EAdvhBClAR1+EIIURLU4QshRElQhy+EECVBHb4QQpQEdfhCCFES1OELIURJ+H/+cn7uATltzQAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# Get a more well-defined 2D mesh (subset of the full grid)\n",
"plt.scatter(ds['lon'][:80,:], ds['lat'][:80,], s=0.2)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Correctly build regridder"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [],
"source": [
"ds_subset = ds.isel(i=slice(0,80), j=slice(0,80))"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [],
"source": [
"ds_out = xe.util.grid_global(5, 4)"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Overwrite existing file: bilinear_80x80_45x72.nc \n",
" You can set reuse_weights=True to save computing time.\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"/Users/zhuangjw/Research/Computing/miniconda3/envs/sci/lib/python3.6/site-packages/xesmf/backend.py:36: UserWarning: Input array is not F_CONTIGUOUS. Will affect performance.\n",
" warnings.warn(\"Input array is not F_CONTIGUOUS. \"\n"
]
}
],
"source": [
"regridder = xe.Regridder(ds_subset, ds_out, 'bilinear') # now works"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"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.6.6"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment