{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "1b5ab34d",
   "metadata": {},
   "source": [
    "## BSSE correction for H2O adsorption on ZnO(10-10)\n",
    "\n",
    "This example reproduces the BSSE analysis from the *Crystals and Surfaces* tutorial.\n",
    "It runs BAND counterpoise calculations for five basis sets and plots corrected vs. uncorrected adsorption energies.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9206f829",
   "metadata": {},
   "source": [
    "## Define the slab+molecule system\n",
    "\n",
    "The full structure is embedded below as an AMS `System` block, so this notebook does not depend on external files.\n",
    "Atoms 1, 2, and 3 are the adsorbate H2O fragment used in the counterpoise setup.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "bd7b85b4",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArwAAAEsCAYAAAAhNGCdAADVxklEQVR4AeydB5xeVZn/z7nl7X16S520SYEkkEAQJqB0paiJbVkr4BZXXdld/euaxG266666rqtiX+smKkVAOhlKICEJCUkmvU+feXu/7fx/550ZBBYUJMm8SZ7z4c7b733O9w6T3/vc33kexmgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQgTdOgL/xXdAeiAARIAJE4I8k8PK/weKP3A99jAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABInB6EJBZ0fHt9Ij4VaLcfMst+nvnT4pOCrNoZ0dd4FXeRk8TASJABIgAESACRIAInEUEXm4BOC2nvooxRQY+pzb8088vOydx/7uuSHz3yqW/kc+JVasqr8n7NIgAESACRODEEKA/rCeGI+2FCBCBU0NAXNUeC8mM6CcvaPWemkOevKMki+XpGmfRtqA/Wu9x+U/ekWjPRIAIEIGzmwAJ3rP7/NPsicBpQWA8I1rnd6+ZE208+K2rrji4pG5K1wWtzItVXlxup8VExmJdP5bhtRzWnS8bIlcsC9u0G1oZ87I1a8RpNJ/TBDuFSQSIwNlOgATv2f4bQPMnAqcBgfXjAtFiOndEbaPfGw1p6txGJ+iD0oU+PC2GImOV23LGHBmxI5yDhmlzx7K5IkRbR2vIezrMZxXOx4oVK9S12FatGrVnnBZngIIkAkTgrCWgnbUzp4kTASJw2hDoGou0YFlmulQW+WKZeTgv1nvqLMayp8s8HChzZVI4HP6FplksHs+WbeEumSYTtsO8imI2qbwihKt9QmukYF+3jq2r9kApPiJABIjAGAESvPSrQASIQLUSUFZ1dipz6+vFw4cO8du3bGFl295WMEyZEbXcTAm3hLSZCP5ZtmoVl1aAKpwIlxnb6bFYqMGj/uyq9tapF7TWN+0dSvzrXz4c/yLUrVOybGZbFtMEU2oMH/4mp6twGi+EhOkw8Y72yQsnxQKRWo8H6V0z9ZnHn38Oz1dee+GddIcIEAEiUEUEyNJQRSeDQiECROAlBJw1XV3WynXrbIhdc+yVRMG0mGGYAhlerd6tVxZ6revulmKrKgfSuvxAIiHipfJFzX5fx5RgIBpT1TYZrO04m0qYj2PbpouxcDDi7pDPr1tRnTYB2BcqnPekUv/dGvA+esX0lkfrPd7/ljGjukTVngMZHw0iQATObgIkeM/u80+zJwJVR0D6Q2VQSxtqlt96zozVX1y++I6vXXbe/WLtWhVPW3mIXcOwGAQiUx07Kt9bxUOsXl0RiVkhxKZMoeQUSoYom9YCGbPBWKmMDK9l2syNDGmdpuijc1lRnVNaMxpWrmzmcQ4cDSbkgKrmqzNYiooIEAEi8DsCJHh/x4LuEQEiUAUE1o8J3l2pzNs8mrrq8ultN8yMBt/y9Lr/kBp3Z8lyMlBbmluaBbiyWIZcNzRUvdnFMZFoWLZuWraiOIL7NXXUq2vbeskeFbzSXxZ1aWU5n2ofNrLWmI8C77HiYoL+Han2E0bxEQEiUJ2Xzei8EAEicPYS6BqbesG0jybyJSubL+FyP8/vPZIJ4yUDlga1DB+vhuJdtZpaPF1IlR2hIrPLuIMFapxjsV1lpEumbRm2yV1wIKfK1gz57AQKeC6rL6z6A80vTMfWZFZazkV3To+FdqO46ScRIAJnKwFatHa2nnmaNxGoXgKQfljN5TjbsKBLs1CZAQIxlORM+lvXl2wnWzKsAMqTMeGIdvneYSxsk7fVOLpHF3Mx03bihmUyLLiTWzNMCy5UOdhlOHaubNgRN9poeBSlMp8JmofM1Drr4Jl+teOPzwU1Gp6DTr/Eln5qy5w5i7EgFg3mcBIqi/Re7fP0PBEgAkRgogjQpaiJIk/HJQJE4FUJjPkTtCIElRRVyPAyf6W/BHSwcPYUx0p54Rt7RSCu6OioWsE7NCZ4Ic+fQ81dZkHwQqy3TZlVA9suE7DCcgPzgdWBRTgfX5z3qmxO4gvSZuH9+OI5S79y+fnnjh3nJVaR8bmgkFrcghVDindk2sOTQyEdb6zac3ASmdGuiQAROE0IkOA9TU4UhUkEzkICFoStMCESdUweIqtGMsjbthulyWSJAxZUlNLpwgUl1VS5QE1AKPoUbp6vKlIg8qLjKPD2YgGeYD7OXzW7erLmuWrMMz055PvsJ8+fs+fd82c8MysU+DXS6dJl8YrDtoUmz4s8Bx7GrdqIbHZHgwgQASJQvQRI8FbvuaHIiEA1E5CZv5dk/05gsKMLuhjrLltOuiAXqGHnsC9UFqhZjniugCoNUiQ6tpi1rKZyOb1q2/F2jYExYHetNJlA3BDwYteALf/+FqBz92IBmPRwSDG/UL79VFo01owJ3oF8OZovm5NUwawal6v53efNmooTjGpjq/7PvxO2EE6FP+LW0DwuPDZHuiECRIAIVCuB//OHrFoDpbiIABGoGgLQQZXL1yc7q2chK+rgcn9lgVpU0yqX+w0hEnnYHKQ1wMNFaJ4TrbTsrRo6rxKIKcSBsrQ0wMyLchNeI+KbhLfaluNkDMxHdlvzYB2e/PiKjnUnm+2Lo5TtjmVTj98cTGSc4WzOjrpd7il+/4XyTcvXr3/h34musU+VhdAMxMuQrUbDDBHSlPEvKWPvoBsiQASIQHUReOEPWXWFRdEQASJQpQQqYldm/W655ZaKODuJcSpFy1YNw2YK8rceJnUiCtfCxZuD5xXdyRxcTlfqvCwkn1998jLOcvdvZFTEKxbh7ZclyNBkwvYy7olqrFXuFOLRVbJHqzcEFD4RZcnGxeqOwUI5NZDOuV1cYUFdWy7jW/7SBYGVuZi2vdXAXGDPsDThhFzCPVO+l401pqjcpx9EgAgQgSoiQIK3ik4GhUIEqpnA2KVt8clPftL7+OOP/8/evXs/Mxbvif47IqDApLDOwL6wz0DbXZn9NC3nfHm8ksOeziO7ayLFi+oNwRpNR5EAxuZWr+Ad937o0sNrSxsABHxYVSqlycpIYsvMr2M5cjFbGFPRoN5PaYb386O2hiS6wW3vz+SZzDgHNfVCrAh0s7VrxwWxxDw+4hC9Mj0t3IJrQU2RcbN13eNTHX8b3RIBIkAEqoPAif6HqjpmRVEQASJwQglIsbtmzRrnuuuuCx44cOALHo/nuN/vz5/Qg7xoZ6tHxau8+p9GCTLGIXh9CnocYMilUlLwyooH8olWzX26LFzTUFINVgxL2gAQu/wJbwhnmyB6Ubas4kmeeQFjIdT2OqWe5PWjglckyubjxzM5lisZTkTXW9/X3lLHORerxny+Ml75TQRDw7mBSLeYikijij5eV3j0VfpJBIgAEagyAiR4q+yEUDhEoNoIjItdNCQI4LL8l9Ei9yFN076nKIrvZMdaFLZLliZjaO2Fy/2jJbtsO5tDEdgSmk94oBZLzJ5xsuN4g/uvCFvsowc+WYRtKXIRXskyFsr9mtC6RQhHR4p6LF2bG5LPntrRVdHdFR9v1/FsgaUKRSvm0n1zYqHzZCRzV6wY07kvxCWkhxfuDDSeEMwtqNvaC2ToDhEgAlVJgARvVZ4WCooIVAeBF4vdQqHwbxC8D9xzzz0PQvTWYXulS90nJPDusUSijYoMJTRrgH1BXj6f8clWhvYMbAe8sLmCZVWqN/ic0Vq8dZ2dLxdlJySWE7CTFwSv6TgF2BdUHVlRP1MqZdYM5rhl5lfWtNWR2I3pqMB26sf4udzeVyyP9GXyLq+isojmWi5DeYXOb4oUvPAjVzK8IdWhDO+pP2d0RCJABF4HARK8rwMWvZUInE0EXi52MfcH77vvvl/L5yF2T2q92PEGB1jQNYLWuxXfq4uJyNxijYY4lLzjqEU8r0Jz1yhK8TQ5LzqsCxzJ6UqTieBYe2HDVg4UpR8WFSlQ+NZXIwKtY/M5lQJefH7MxztSLu/pz+UrZd+8Gn8TYlFRqWH8fI+L932oLpGFeNdcEO9DZbFIxlw3VLVfOsaQ0g0RIAJnKwESvGfrmad5E4HfQ+CVxO699957Byoz6NLL+3s+ekJfKtvCi0xuJZPo5YoRH917wXTE3lKllJeAJcCpiK3hrq5xMXZCYziROys6QpE1d0ebTEg7MmMWdw7I6g1YzIbqDcwTZkpF8K57wS57IiN49X2tH/fxmtZz0taQLRtOSNVm/2NrrFH6eAFXCnCYSCojbaIUL7LuXIOlwStYJVs9+hL9JAJEgAhUHwESvNV3TigiIjChBF5N7MLDqzY1NY1n+k5qjF1jey8L+0ABmVwUZKhkPy2vMxkv2ShNls3D6sAgtjx8VGytgBg7qUH98Tt/oeoEXCD7ZMMGgSYThbEmEwCqF2EPkK16ZfWGBkUx/vhD/fGfBPOKmAXvR4/liyxdLJu1Lpe/PRKtVMdYt4JV/r0Yg4xFawLZ6lFLg58LsjT88ejpk0SACJwCAiR4TwFkOgQROF0I/D6xu27dulMidsdYVcQXdOH+yuV+KF4fZ56oUNvk60WorjyEsCxX5h8VumOJx+okvXo0O2ohM52W1SVGW/K+kBU1kfkVBsqS6cidu53ResMTMJPxzP2zPYVSHrYGHZUxWERVlstYXmZXEAZOjiyxpmJBIS1am4CzRYckAkTgdREgwfu6cNGbicAZTYBLu4KsxiAXqGGmD0obg8zsnmKx+2LIXmRCmWxyIIt4RZla8evC5PCMfF7WtUX6dO61jEWgeJHvHa+a9eJdVM/9ssNUaV9AvV0IdT6aybXtXRCPGTSf0N1IUmcZWyojrjv1cxnPRA8MG9aevnxRQWKdBTRV+ni1F/l4ZXgcHe/wOpqCQPiiWsZYhnc8Ny/fQoMIEAEiUD0ESPBWz7mgSIjARBKQGVLxvve9L1QlYnfsyjkqZFmOhW5r3I0A80xUSpBBOCLDi8oNEFxeJtRpyJlOJLzXemxYNFQp0mVLXmRF5QI8HsSj8liGV3p7IwK9HLBQbN+pF7xsOY6LmOy4bT15PFdkWdRADinqjL8OBCJjPl4G9PJ3JYOCEvuk4JVZ9tyYPWN5V9XaShAyDSJABM5mAiR4z+azT3MnAi8iADuDViqVPo4KDA/JzG5nZ6c2gZnd8cVR+wzHzqAWL6oBMFmTDM2/UC+W80R21PdqQx16JgVclYVeeImvhWiTghFb1fx9e6HMGmPPyWy1bNqAUmRzb0GTCWR0s6nKY4u5pCeZiSwUpY3XZMb6lM6ha8zHizrH648USixTKhu1uh64pCE8X3Jfh3hWjwpeJH9FGt2dEaUj467By3IgZBpEgAgQgeojcEr/mFbf9CkiIkAEQEBm7ER3d7c/m81OGxwc/I308nZ1dU1o1nRMOalFWyho2FCpbFCnKvlKvEI8n4dINKEcfYx7WwVvwvtldlJZCbEoBSO2cU8qnp7YMTTKmJmCj0hLg2wygcoG+ncYS78lFps+ye3SucNEBiqyVvdevaW5+bq5mI+cw6pTK3rHmT17vFTO9RWKWlBVlbCidEqCL651bFa6rUl7hpzL/1m0Jn+naBABIkAEqoYACd6qORUUCBGYWAIQu8LlchVbW1v9p7L02B+YtVNwHKuE0mTSK4pStQ14v4B/NzuE5+B7rbQXdgtVikMpcq1/VdWvHnC5fiy83krWF0K4av7Ooa6wSzaZEBC9Qc5y76+v//g1Hs9TfxOJ+jpcHsxF5XO9/uWLVPWuR5ubH3+4qWnRmlOb6f2dj9ey9kkfryP9xqp6Cbjz4fouMZ6tNlFyQmaquVy0Nibo5bmUGXZ5jqqJu4yLBhEgAmc3gar5h+DsPg00eyJQHQRgZ1DKZemQnfAhPj8qVNOQh91ocFCRfVjUNYM1Nflura9/zxKfH+1thZKHjXRROPyjQ83NX/+i3/u9K1z6x6e7XX9SEM6dEF+yM5u8zj6hGceuMZwQiYbM8Ar4kksOq58Ecd7OebQeer3N5+PRaJSxWNRhwaDV4PEsPUdR7r9n8uRGOYVVp0i4I0Mu/12whyx753HYGrLwSiPLe84nQyy6ch2zO8BSxoJZuNB8Au90UFlCSPuFFLpcZthxi2px1ZNhRzw0iAAROMsJkOA9y38BaPpEoMoJiCK6ulVEIlKNgvHIBxi770q3+2Pv8AfEVN2tNnm8fLLP1zjV5frLPw9HP9To81llCMRHHGfmXahxWyXCC3pQakPn2RKkIXy8ilc4aj1S1h5FEVG3m7kjEcYaoW3b2hTW0qKxcNhAHdy6uYbxRcxBrD5Foh0+XRwOq9Is8diRUhkNKMpGjapGL/LWnCOf729irjUQs4Om/WQeC/DkwsGC7czEhypWkv/n1v5yk0vvzrndX5EiGBv9OyPB0SACRGBCCchVwjSIABEgAlVLAPlmWBoqncicIuPnz+CchR3HDrldaq0/wJjfD1kFXVUo2MFslisKVx43jOKPhP2xX0G3rcCrUoxN8AQ5hB8iZ+lMRSSilBeELx4riJ3pELyVecgMr5xPschgUNZZPu8EDePddzY3/wvr6zsg7QJjGdSTOR0HcSEr7jx9tGwafYWyOi2A7LOmdK5irGtNPyt8pKFhrs+yFrYpXGCBm6hX9dZnm5q+vzOfbbxEiKunyZnazieOCtdTUwzjl1L4ymdOZtC0byJABIjA7yNAgvf30aHXiAARmDAC60czg44txLPwvV6OhV4CQlcEBRdul0sNeeBWkFnRGhQIcLk4BK/mDAw4/oSjRHz+3b/K5n4ks4tVILSkfrTkj0vD4QX1wmZFiF5dKGyq4qDIAcp6QSDGpEgc39DwQWga56rqeBXFbdp8Bl7dK/chxSNu2GrcXXNybAPSqCAjOTxgm4d7i6VZwu9lEUXtlMd7d1PTJ+Zw/q8X6rrejJbDNYqi+Lxupgjng/Mh3gswbYyUSnwPstibLSskY0XWmAYRIAJEYEIJkOCdUPx0cCJABF6NQA6iC+KOu4WdR/2rijD0ISlqoiCsDWHFdR0OXYjeYHD0VtOYggwvy2TEFE2b/cv6+mnrh4aOPQbRe+lodhG7O+VDCkextLU1doHj/GQJ51fPwuNGiMSworKpkK9Z02TxbJblAwHmz2QYM4zKxpHldSyLx2GsrjGLn3zE7y58K1/egv2lx2cBY6+yfj22LuasOYHid8zHa8Rt8cyxUmmWtC7Uq8rkZbW1n2xT1f+YCkkcVBSnwetVPDIjLc+FaVpuzMOBTL+nVC58xzK/9JDj/FCeQ8RM2d3xk0a3RIAITAgBErwTgp0OSgSIwB8goGyBhIJQYud7fC1+iMMCRJeHK6wNz+UhuNJw9IYhhCsCEQJYNnOANpTiyoHo8sc1rfGdjB0aPw6ErzYMsQkxhwSmTJSe9CEXd/E10Wio2bbvn6Oq57egskGTy6U0SaHu8VQCCBcKLJDLsf7BQdTntVlYikfMy8JzI4ahDOIzF4f8lzkBz2UzLKtnjWU9FLfMe57KGZv4mmIPdiLnIyfE13d2qsNdXW94jsjISvQsaYn1R8rm+1Nl0wkxpbldVb4cwgI1+I55yO1WPCEkcGWGXYpew9DskRHhHR7mSwLB3IplF/2rev/9MrbKvuT+aBABIkAEJooACd6JIk/HJQJE4NUISIHkXICs6MXIip7P+ZXt0HMeoWiNuofNQuGsOLKiIxCECYjGmBS7Lhe6UZQZz+chjC0+DOXb6Zhf2RkNrz1mFn97Tc7Yeylk5PgBH0NTDSkM4e+VftWTIn5RzUBfAxk4y+X66yaIXZ/jGD5dd4WlOJRe3XAYuWdYF5DVVYeGWCydZt3DwywAP68LwrdkGiyB9g59lqnkC6aY5vNaNUFvKzKsHzQd+4OTIkb2ett+dAhNIo6Z5Qf5QKKbdXW9MEchW0LDTLBi3R+eI4S5wjo7lbn19WJFR4dY3b1OWc3mqu516zYeNszyQLnsmuzS3XPxJSMOYA7iVpFRr2TYpeiVGzLTimHIDLs92TAaHtm27QNQu99eP2rBeCGu8XNAt0SACBCBU0mABO+ppE3HIgJE4A8RqGRF/6O5OTbVtn/bAaHYDPHaCo9unQ9ZURcWdyG7W4/L/WEI3kMjI6wEoRWVohdZ0RKypf0QXZBd1rxQcAkOtqTRNL50IGrvS1rmvYOWfcdbh5MbL/0jheEfCn7s9UpVgm5IwMWLF+tGX9+tOoy6yDrrKuJU5QK18cwoHnMpHBG3D3PiEO3bMQ80cmAo+WUnFUX9ZaG4J14qti9Cc4rZHjdrBwdsRm3AH2zR9evbhXP9lLJh76qNbExa1oP9BeuhFcf6NvJ1616wEUjxC3sHX9/V9X+sD7LJSKXuMl570fyMNQwzYGz3gCOG+0pGa7+iOksVleeRRXfhHKQQYxPeoECcy+y6rLwLt4kU8fKO8HDeKnewXP6gQQSIABGYYAIkeCf4BNDhiQAR+B0BZFyVNZBP84T4jxaIXS+yokF0wwjB38pisVG/LoSVSKWYG+KqCaJ3K+6jMQLTIMLKpingeeX7SkVtV6nA5vh9VovXw9uC+pzpnM+B6L3tUCzUjUoJDx4yjPX3HB96FMIwOx6BgPhbv349PLH/VxiOv2fslq+CfWA8I1p5bm43X72rQ0jxCNnHpvvdl+/fu/vvJgXD9TK9Wcal/RJihkNXNmoYHaMCETlQCHZkTeXnNqKurQ+yHo/VQ7a9+X6fbznLZqc8aztLGkvG9c3Z/KUzXHqoAwvFZvh9OI7XavKhNJsnuAw1fZeNGObq3Q3R7mGzvGmkbN35q8Ejj2GOMAePDqlL1y/vVOUcV0OmcsT77vPPmTszGnxnrcd1YdTr0gO6tiNr2ff+5tEtS5apak0E7ytYQjkfYh2xMdn5bhh+3SE8bpQCHhn3iujFczLTm8RB0nIOGOvlDxpEgAgQgQkmQIJ3gk8AHZ4IEIEXCChsxQrmeeqpNo8Q74CAdeDI1eWlc1V6XsezohCJqF5QKd0VkQu7ILC24zlkRW2ISfUIpNcv8oVvOqZ5WXuu0DHXrbMOfH6G32tP8vtEnc/fMVnlHZNt6xPzY5Gev7DMh+Ml645Hh1JPQfzFEU1FqL2a+IUoh1WA2WtelCV+YQa4M782vLg3X/r6gYJxYWPAz/LIeA4gvzsZrw1CzDYiXvRwZoqcg8zuQrSzUokVISItCMUGy3A2IZu6y1a+unew/1Ny30gZ70IL313Yzw8GLNa21bLPjxZKb21M59481aVPmuPzsJkBH5sW9NttmGNtJNzRpiodecv6wKzamJzjpqRp3LM/kXmE8/5jjI1aH/Dlgt2wcN5XZkSDfzG3LqI3oRpD2O1i9QHvZT98dvfHJwnOZiNONMhg02Af0eR5QMxuzKEVcxhBhr0HMdfheZnptXA+UpjLIOZ7hPOsjJ0GESACRKAaCJDgrYazQDEQASLAFi+GkFy3zlRratq5x+M3sFjLwDXyEtKeMkOKIl3IkY5tEGAo28VQdLdiwH0cSd8wV9Q0Gn5BbH3sQKn0TfkRZEiXdJnWjU2FcmdbOrtorteldcASMMPvt6dCHNaEQq1tmvqBvG19YFJdJH6DYT6VMow7tyVTD0H89mAfo+J3PCuKJ6TQfVN7e90FbXXvbwm4L6/1ejS/rsKCwLbnyuVvfujXD/990KVf+PbGGjPmdikb82W107HYNYi5D0XS9kMUyrnUwb7gkvOAQCzBe3zUsoRf18TfeXU+YFsKLBrzhhpjn9tdNh57dzK3CYdGGrUyjmuOfTwp2K+Ttu3H6+c8Vijd2JjKXNvqds2c7fNqs4IQvyE/mxIKGjWhYGuTS2s1Heftk2pihSemtm7LlI278qXSXV/LO7e2x0KfaA34RNjjtqLBIJ9SE2UHUxnx+M5DyhWaImRxNCxQY5osASez7DKjizkILK6rheDdk0iwXgjegMywYw5FeKj3CmHt5fwuGe36MYbyPg0iQASIwEQRIME7UeTpuESACLxAYIVcYAWx+8ELFq0sM2f1lqODTg5iN4V3xJEVTUAcNstMKARu5fI/7nNkGXPILnqQCPY4NnuW8U3PmvzDLDm4U+4Y+VML+nIDupptOILs4xHT7HiqbFzTlC3e2KKnlyIr6pqDjGZ70C+mh4JmbSBQ0+zSroMwvK6tNpa9bOrkR0eM4vqjucKDnB+DoXU0K7p87qwr59eEvz+nNtzcBmFZB1tBAIIXkV1Wtp2P3rqow+ztG7STpqWpuSL/esDLFru8iBO+gHIJApBXSqxJH3IAwlcgtgLm0oesaMEwuenRRZvXI2pc+hXou3ZFG+wCB2PR3SnLuveYadz/j33xzVsqjgE5S5aHIN0Ai8GGw7b9/w6Xzfanc8VrGxLpq5o8riXtfk9wVjDApocDbEo4ZNQFvN7pHveyoK4uu39g5J985bQW1BTbo2tKyO/T6mJRFpg2hd11x0PMLJaYC6JZfqnQ5KJAudhOLrST9hLEXrEx4DzEsG2GBxnzR4VkYWYVxbXdcb7x3319+9fiNKxENrwSKf0gAkSACEwgARK8EwifDk0EiACDi6Eidu2LZs24rSno/TfJ5FHbFkfReKFJ1qjF4wMQWGoyyWLIhsrsqMBjmRWFyBOtbl18V9F5yjJbCqr92aMNsYf3G4WH3p8sHbPFKGHoYu5SeDda+3b3Os6XeyF+N5WMpQ3p/PXNrtSlM3yeUAfEa3sowKaFg0a93xts9LivR53c61sLRfuppoZnLWH/73eHM715t+t/p4YDPOzSzUjAr9RHQqw2HGIuCNWHN+/w7kumvUeSWXZrJMBuaq5lAQjcnGmx3YWCc4+BhXWaW8lDwA9hHtIPq4xaMZTd8B7flU7ta8nwmQvcbnXOqEfXbIM/tzngnzNN4XOmWdZtHbBhwF/7UMI07vzlwb6u24VIj86SmWCD3YjdvYb5ZWytmzP5y2sSmbc2ul1Lpvg8rbMhfCdhjrNjIevheE6TFgqkzeEawWI6DRlblBjbM5Ji33jsGdaJLwM5dIMrQslC3kLO/g5mJdM+9uVDx0sjyMYfk1OB2D3sOE/dZRifWYXHELuVDLn8OA0iQASIwEQSIME7kfTp2ETgLCcg/bC/RDWBKa2Nl9T7Pf8W0FS7JeDngaBfKSWSbDkW/HegG9nmImfPIxPahEwiXKSyIQPLIyvaI2sf4BK65XXbDYFAS7OmvbvJsd89xQwVdkXMbSnbvqPXMO7jHCW7xqqPwRQhRXM3ysl2D5rODwZNs/W5fPHyaCr31hZ3csk0KQyR2ZwRDrLp0ZBRH/Br02qiF6CX2wWejIlCES7mVhXH53HrETS9qGtqYmnUCf7sL+5i31+/SVwV9LK72lu5rLKbypZYDSTfDtMQn7VspdfrYzMQfgL+ZNlZDe/h8Cmrw3hvP2M3by+bP9quaUseMws3NOWL105J5+bM9bjYHNgw4NG1pwT8dl043DpNUz8If+4Hp8WiPTeb1qa4Yd6zLZ166NO9iR7sqjJUzuX9H8TLptxCuzL5S+8fiC+cGQm+f1lL/RQRCDgLwwG0wOAsb6IMGjLpDmwi8YNHmSrZ6i422zJYGnEO4ItHDOz9WCAoy7/JhWmyBJyN81AA2AMQt48K0eOyrG9vHxj4Eo5rrsHccDumknGPBhEgAkRgAgmQ4J1A+HRoInA2E5CZ3eihQ4rYssXRFf0LARc8ucgi1tWElXfUhNiixAhbHAwxFWLLLuYZfKEsDqEFBynDAi4bhlZ1i22bDyWS+dm6GpkLb+ksLN5C5QIDHcC8M8KhZXD4Lusvl7/YHYtuHbKMruFS6Y4Vx4Y3mfCYjrOH+O2BKvtB0jTlFtqZKyxfn8jc0OiJXz7V721tD/nYrHDQNnFZf1Mqzy5sciuQ2YpcTKfBx6pPamPfWHcf+zmyoj+cO40vcuns18NJ9vPhNLsXl//dmuZcJJQMd/m+m1T4CijvycchM2EMYPIPMLKsuyB4v3Cwv38tkqlMWNZTyKg+dYixTx+yrEWPl43OlnzxxuZkeikqM7g6ZGWGoB8e3YBRF4Q/V4c/Vzhvn1IbyV4+pfXRwXJp/dFC6cE/65Y2jNEBd0gGmd+7Prx8SWuLx6NPCQdEW8DLpffYkJaRUhmWYZM9+PDj7Mkt21kNFvotKBfYdUjQ1uFLR3dZYwexSG0yvnT4MWcFQteA4IXVxIGqVoqO86ndpvkdFo+j1V1F5ZLYHYdPt0SACFQFARK8VXEaKAgicHYRQDEG6dmV3k57zpQpk1WFLZZFGWxNU57ddYAt6ulll8YiFShPFky2lqssjEdJLIzSIXwhONVBPEZm8Zatpnn/Vse5qKZsXtuUL1zensq2zoUdYJb05wZQlgyWgLaw63yUJTt/xDRv666v7Zb1ao8VSut/tO/4Y/cLkakcCD8gfuX9uzOWdXcmZ0X25Qoznsnk/uJNgr9/aizsXNRcozRDVMPnC5uCjYyozoaP9bLZ8SH2o3nT2J5CmX3+QA87VDaYS9ajVRQbmwbTbvgxs/xVrvs+73ecZQO2fRHmkMHCvK25WOwZ1t1tjFV/kBYApROa/vHRtXrPoh7vs0cN58tHDbPj6ZJxRXOmsLwpnrpsps8TnAN/8IxQEDaMgFHv8wUbPK7rpwpx/eRiyd7Y2LAxbZp3HM4Wuh7uTe7ONNZ8rcnj+dC0sJ81I4PdgAx2FHV9w5iPzNj+ePNO9pX1m5gG+8XHYH243udG+2AsxkNGt1wssN0ILIf3hWFlQF1hVkQJOGSZ2S7HSSz1+b73g0OHsjej2cbto3FTZhe8aBABIlA9BEjwVs+5oEiIwFlBAJlGJBy5/a7z5795VjT8VyiFtcyncW9dMMB3HR9g7n2H2EUxeGLRRnhzIs2+gFzsDiji2kKRhSAifarqcIXvHmHstr2Dgw9UsqK2/St4fX8Vt+3QTsO59OFS+YaWbG4JuoN1dHg9bDZKdrUHA3YbqhHURyIdkzWlY5ppfWJefW3P36MsWaJs3nH/YN9T3+jLSsvw+Eh98KLzaicFfedhc1qCPo46tZXr9KliGb5Whx04eJj13HUfe75viN2H7O+BPHKzqLzgRfYXWU9+QyajfS0QGLpE0z7DE4k+JZEQxxh7BAeQW2XwgQH2Tnxq3e8Wd9ldYy/hRmZKoZyZjUV43fhS0I3FbV/tgw1jS758eU0q+9ZWd2rJVJ+7dQ6yvu3ShhEJGY0BnzrZ614GN8KydnRsG3R7ep6zeWtlgZpLV2KoxtBSV8fCTfXsYN8gW3X3g+xX23azd0SC7OONtWwmPL2wR7NeMH+4WHL26B7ebFocXwTgSYaqxUnEMFPw7D7P2L//16FD6VX4vrDmd5UkRmdAP4kAESACVUKABG+VnAgKgwicDQRWobEDxK5z+bxZX5kWDn1iNqwLdRCk0+ti7OhwgoUPH2NvqY8xPwTj8/EUux2i61PveRtLZPNsW0+/UNBW2MV5/mfZ5PK+Z/fEx6oAVLKiyJCyXzKWgTa8K2ezu/batrq3bC7tKpRvbEpnr21zJ2eiMoM2WwrDUMCeFArateFga5uuoSyZ/YG2umj8XTPNp3CZ/q7+QvmB/yk71/pU5dvNyHQ2+D2iMRTkdairG0b22AdBuO1oD/vqI8+w9Uf62DEsRpPS1K1raH6B0lzw6aJ/Q3qbYXy9Mx7/JkIbhGiF/hwVsHgshawccimY8yKxO/rs6E+8VPHAVhpZ4L7U9gqEpfMPDOVvmfODODzIcdMKbYcN49FE9oYmT+LyyuI06UHGYrppsZDZGg4quw2n1efS0BFYUV2wJLhRUzfY3Mju2n2Qvf1fv8Wmwsf79alN7CrMTy6kS5cMtiufF3eUTOdBzaM262421bHtWnRQg+BFvp2pRYjdY7a9dm1//xcRkwKxCwg0iAARIALVSYAEb3WeF4qKCJxxBKSNAV3I7I6pk/5iciT0iQaf2/ZAOLbW1Sq9+QLf/Piz7CrYEEK4zL47mWG3Q3Rd+pZlbBY8tIZHY/ODbidVKCk7hhP7brYb0qyzQVv5u+YPNkSjHFJIKmPi18b9SskulCX79JGyOQslu65pSGZvbHbrS2cFvK7ZqFgwPRwUU8MoS+b31zTGXNdpnF335HAy7xtI+WOovAAPrgj7A0p9TYw1NjWwNHysX77vUfbfj21kRzK5SkbXhQw0yp9JscsgMO3l7cgh6+pHfvLM9jug8NnPhFARmIxHjvHb0Uev/acU9rI1sBwV8Yt5ilGR/4INI7g3V1jyZDx7Vb1n5NopAe+cJojY3ehDt7A+qpioumBhg52CKVhwt23/EVgXXOzvZ09lzXgOGWR2CHN6MlNwtoeCSrw5qFoD8T2HLXvGCPoih3BgNzZHiCHU2/06Fqj9I1S5FPI0iAARIAJVTYAEb1WfHgqOCJwxBNDdljmBQFNtQNf/IYDar6giwGPBkFKERl3/yFPsKghb2ZnsaCrLvpvJs4svPo/5oaQe33eYSbVoWradNky1O576xk82brfWYtHbK9CRWdGXiF88loJMZh+74Yft7i0bX8bW8WyuuLQhnr4eJbsuhdc3hIVpbAou6Z9bHzO3lyw/bBNCRwkDVVO5G5UZQuj0ZoQj7NrP/RvbsPdgReh6sECtZJiVdsFNEM6dU1vMhWjz0JsrbLvr8MBvFy5erG/ZssVe+ceL3FeYYuWpivh9BZGfxauP5G3zkcN58zOH88XzPbnyh6bW1364aNlC1vnNWzZLlbFA7TcPsHk9R9l7F8yotGaWNo0tKEm2EU08Nnhl9QZ9iKdy//WR6b3/8u29jTOSQlwID3WrqSgbIPqfZ4ODQ2NiVzKnQQSIABGoagIkeKv69FBwROCMISDFqeXx2TcgYxqV69N0NDtIYhHUc+s3sEux+KoNC6UOjiTZf/cPsz95x9VsdtDDfrrzoEgYpq3jknsZboaDyezDd27d8T/SGrES2eI/QKcifuV7pDDDjdwqflikR7sFxO9gGWXJymbr9kq92vTF9S7X5TOO97ceV12iORZmaCSBBWqCGXJP8OaW0IxhC6wMlYGMbgnbNCyua2+oEdNqo9yru1xP9w33bR8cev/R3qHSB44eVbYgITr6gZP2szLPF4vfTsx1OY77Bcaero1GbYjdj6TKhjMCwavDKpLFHCahCsM5aEqB2hhsbzLNnkhkGGtucFCCjW/vPvKFrt37voGIhzYewE82sBs/5PbCQHb5xVnrF56nO0SACBCBaiRAgrcazwrFRATOUAJI6k42pNrFiqeBfIkd2X+YXebV2ZTaKOtJpdmTqs6uf9d17JJJjWzXkWNo3cX5c8MpGB84S5eNXzzTO/AxqDuBtr+vl5AUhXL7P35YiEKpYH8QLxlyC+7OFa5taqj/UZ3l6LJhRA4Z3AyEbm88CftDmH3xndeyf7jjftYMgf6+c2axhRC739pxiD/ZM2wojP90fyb7uWI83od9KojyZItdHOYlQ86xsuita9T2oCqGMZLhSvZwphBoQSfjZqvE58M6UhurZyksstsyMMJ2wis945w5LAx/c3ZwhP9bvX7PkoNs6GNXtbu/fv8BFN4dtVCMHanC8VV8xy8Jhh4QASJABKqFAAneajkTFAcROEsIFLDafwSX1JXsALsCTRpm1kRYfyrD1hcN9ic3rWSNWBTWffCwKODSe94wR7YNJe/ROf9Ff3//gxIR0rT4ryJe3wixiiUAglQOJHwrpcDY44xlUQv4F0jqrk6UjVmpkuEk4BsezGSxIE3HUQX74OK57KqpLUwYqF1bKjndQ3E+ldnP3TmSeA/LZvfKHa7C/rDvUy125aFfPJxKHENDhxZPnXTf8trAu66sC9utfi9KMDhsTzzNNqBecBEiflZLI3MB6ba+ETbfLMPDa0SxI35x8ID19dE9Vni9eOd0nwgQASJwOhEgwXs6nS2KlQic5gQ0RTk6gsVoIpvhb2mKsVnREEsiy9iVzrNpi85h2XSaHT2aZvli0U4Xy+h4y/97ZGBgtVS4q1ZBRK6pCF2ZYTyRoyLmukb3qCH57EDYbkqUzFmD+ZLjg/VCw6I0Gyo4i5bGRwd15sHF/DyyvkO5vMjmi8pit5qQYvfbt9yi33r77VYViF05Gym6xYK66G2z3fyKG6Y2OvVYEJjMFpwdg3HlyUSO6dEIa44ErOOpjHo4V2ILUIh4ssqyzyQNaV8Qu9a94S8Wo1TpJxEgAkRgggmQ4J3gE0CHJwJnCYGKlaA/l7vjsrrwV97eWhOYHotgEZXBn0CtXRPteSE02XMoSyYrBZSxOO14tmAcTqV/DT7qosWLIXa3mKeAVUX8ovTZfxZM891HsnkO0Q0TBlNKsDekUJcWAhgpZqzsQpxHIRJnOhbjaB2M2JRbmprsW9949vlETFN+R3DmTwpHPab6pc1DKeXn+3ust09t1u4+1Mt+O5Rhqs/HGrJF1lMoa1nYNiY7trisrZ4nCgUzUyqlZRCrMReIZhpEgAgQgdOeAAne0/4U0gSIQPUTwAInLhdVXRhwv++6tjrfuS11drlk8qcGE/xpU7ApWPy1s29QuFWsTkOGdcSwlEIydfjo1u5uaEtnC9/yhxaonSgIUvCqxwcGNjc3Nv5T1rRX70/lWMlyGAQwS6EFrxuNL2Qx4SSsGaJUFu+I+tkhoWzE55z169fLv6lyHxM9KgLcVZfOOT01D9Z4XFdtH05qyUJx7W+ODReWNtV9AK2cRaJQNHaUrVVWMnnOVW1173FBJqOSg3rQsujfhok+g3R8IkAETigB+qN2QnHSzogAEXg5AbmaXy5wWtwQW2oK8bUf7Dlq+dwuDYWv2Ld7Eszl9bL8UFIEXHDqQnANG5ZYhGK484zyzm+issPalSvUlazShvjluz5Zj6W4VvsGBr7Q1NSUSJedz2XNbN3BbIH5dY1H3bqCOrRMFAr2zfURFkTQWVmAt/qG2LKFmR2t5vtUW/8roySO/bzn2PdRieKvL/a72A3tbfx4Km2+7bEtX5rkcf/KjTbCOqS66oh9z2QyGWTccT74ibaPVB8liogIEIGzgoBcrEGDCBABInDSCHSMXeK3mMhJ9QSxqK3de9T+7t5jzxWROc0VS6I7keGb+kdKO/tH2HmWKW5EownDtJ6TQdUNDUEGn/JRyShjodzXS8XikqhlZC8JefgVER/zlEoPP9M7+JZwLpeNqYpaRuY3Z5WrUvBKat09mcTG/vjq50YS35c+aJeuniMrE5umzXIlI4O3eHyKEvHjXwMNNg3HsuRzFlstv37QIAJEgAicGQRI8J4Z55FmQQSqlsCa0Uv8yvbB5C7YRN/pU9XHEiXjI3v6h+/609oA++qcSewvm8I9xwaH/v48u8zeGvCwDKojjBRLT8tJDXd1TVSWkcvsdCqVck3ijv6+xph9S2st+5um0PftdPqgV1UjXs4tC8ldjxDS0sCG6ycsVnn4Vxu8E1r2FsZ0uejPxdg0P5wjbnhFPJztwYdKQU1xebnCFFRvyCPD/mo7oueJABEgAqcrAbI0nK5njuImAqcXAelrRU3dxK9wKze2MBrcGkGnshi6mU1S+bNey1InuaUMYyoEb2Z7Iv+8fB9E50R5YvkvRzukTdNV1YsGFOZQoawajO1DWA3eimjEG5ApNUysXJNDGpWrb4guZGyXY1EdQhPI5hZ8Y+LWMiph64oQc714UXVs5mFORbyvX9+J93dNFPvqo0gREQEicFoToAzvaX36KHgicFoREDLTCEGr+BhrjGjq9BgqHuASujicLeyN6up59XjshXcUC6cO3ZPLpWCVldfVJyzjOHbggg/i1gPUHIvrcqWy/LvJ/YjMhTcYiH+oXK5YIKr4bPDVoxy9LoW1SbGuyGoYljWCmB3MT0O2mgnUPkZju1HxXsWTodCIABEgAq+XAAne10uM3k8EiMAfTQCZRgEB6xQYuvG6XIGIqlqlssH3ZvN7GjVtYQzte+VldQje7TiIxVZWspJ/9PFOxAddqlrrxwI7WAGYaZjCMe2CoihL5aI7XTDdtOz0cKrULY81gdnoPzhV/LGX+t2HomqtHohbBeIdC+62yg/6uaJ4UGzNgeBNm+ML8HC2aBABIkAEzhACJHjPkBNJ0yACpwOBTmRGZZy1Lm1xM6odQEgqhZKRfTKRMppdemtUUeCJRStfo/ycfN/6odH3y/uneozH6lbYItgAkJoWGhaoFd4zc9GekKL4ffjrqSMFLfVhWc2dLllRI6Cqxng2V2bXwXU2Mrw+v9TzlvQjs1E/ctfEZdZP9bmm4xEBInDmEyDBe+afY5ohEagaAvVj9oSIpnU2wr8L36himOahwzkj1ObS3CFcaketWytTZk/JoJd3TZh/9wVmOldNaQHQHcFty3HUdevsgKboPiSf8Zy0NFjb0i+8vVrvVL5oILhWiF2/T5ozkM3NGYa0NCiwZ3A37sjnigZZGqr1JFJcRIAI/PEESPD+8ezok0SACLxOAmtHF6C5YpraUatpzAUfad4wun1CtLXisVswDTVtkw+MjOwf27XMQE7ogEB0+WABUJHNRak0jlVcOjK7S+RzLjyn26J7O2NpBApZXLVZ0XHB24ZsrgfVGRwTmfS5mrYFcBsCWMQmPcooJyz67Kr3I0/o7wMdnAgQgdOTAAne0/O8UdRE4HQkAGNARRA21qpqGxatCQE7QG+hvAuL1c6vh3/XjYxp3jT3bmJM1uydUAG5fKw6hMrFEmkBcEPcKsKRFRosv6J6fJCQKvzGSIjKBWsTLsz/0C/EWIB5WZJM+pFl2HFH4JQoCwOYi+YIDcnq9HApV/V+5D80V3qdCBABIvByAiR4X06EHhMBInCyCFTadqE17+wml+5D0wYrXzbYs5nC0UZNmRuF4IUKYwXL3oYAnPVQYicrkNezX48yluGFGNdsJyE/62UCLSg44yjjlT89qhpUMrxuRZkZkOId9oyyYdlD+WI6qDBTVpyAR5mjQgPvzbFqrzjxek4fvZcIEAEiUCFQFf+g0LkgAkTgzCfQObZgDXV3lzXrKsqPcQ3VGOIbM3nRomqNMc7NsmGyhGFVKgesrxIkXqZYMsMrS5Lhkj9UOdNVzmbCFlDJ8HqEg4Q0Fth1VodAfyVsYF8ZsGHMkBUndCR0UZIsf9OxoUMxlbf48KoLgt607cwBZLCRDa4I5LGP0Q0RIAJE4LQnQIL3tD+FNAEicHoQWD5mEajR+IUNst4uKpTBR7ovbts1kzVVRZUAJVU2jFSuWFmwhlnBLjthY7xubbAibqH/pH1BOOIgIjICqGLrlZoQtgBDljaQo6vys6p/+FSeq2Rzpbi1bG02dC7KlC2s1BSGn1q3nb17GctiEpD41W/TqGrYFBwRIAJVRYAEb1WdDgqGCJyxBPgXRgUsrAzabHh4pbhCzVdre0CIya147BVMzZrm0H2G0SsprJ5gwYU/jtL26sJirnCl6QQEYUSIPXjOD/+ujgoTjoOsbxr1DvBcVY/x6hg6U2ZIcSu90ihJ1os5iAjaI/sh3mUjirxty2INNIgAESACZxwBErxn3CmlCRGBqiQwfom8rU5RGmOK6kj7wuGSkWjT1cvCCheyrW3SNJ/bzlgBShOJ1YnNMI4FbMAGYEr7gszmlm1bWhpm+jnz+gSz4XllHsveKImvR4Zaxr0Wm3xcTaNjjKXORLu0L8CvyxD/wS2ovesIZwqEbyWDXYatV85hfZX4pxEWDSJABIjACSFAgveEYKSdEAEi8PsIrMAlcpkuRdODjhZddYU5F8PlMpsttNv+M1Z77kLV5QRMhzVwPTUmdJGDlJ6BCRuVeHF02ZXMj+xzpW4tFqilIlhaF0BoHohGGxleGyW+8D6+BlYHBGyvxIbYq+pv62oEKEdA4QWZ4eUQ7ynT8n+qoWH5XJd7CjK8QrYVRuZ3j5zD8tE5TCT/Srz0gwgQASJwoghU1R/lEzUp2g8RIALVRWBoTLzWKMryFmlfQMaxZFhsnq67LvT7xexAUG30+MSVgeBNPc3N9/9rTU0AM5Cic6JE1/hxZaMGD+K1LQjCt6ns4RRjQZkhRc1gVjQt57l8OYM/pOJ7aKbRpWnfv8/luhEfdh5DIrVKzoKcC6JlXjfjbVLwmpblZG1x0UWq9tjfR2Puc3U3BLzivCkQ+vjR5ub//M9YLIgPYZYTxr9K0FEYRIAInCkESPCeKWeS5kEEqpjA+rEFaLWcLWyU5ccsmyvIg7YEAoLV1HDR2Ij2Bw2cBYN2i9t95QqX64tSNGJKE/Y3SipEjBLsC2gwwVRUNXBuN9hHV9XWfuUCSGDUrVWRETXfU1tbh0BbnhDi3h5V+aBLU359Z8B386WY5WZUdBjdzcT+BHE5HR8sDa3S0mDZDo9yrrSrigPerCEU4sFIRAkFAoFJLtfHVng8D323piaIt07kl46JhUZHJwJE4IwiMGH/mJxRFGkyRIAIvCqBztFMp1Lr8zXVqOp8eHhZybKUoKqxQDjMRW0t4y0tkIwtzKmrU5jLZTWo6p8/0NCwFFrTXjuBnliXymorFQxQoxYNMZQGxv/mfV7v+df7g6xNdSkzvH73ebp+709rau75XyF8dwphdKOAmVtVbv9NNHzzefDIjsWPqUzsiCEW1OA1UA5OLljDijt81/B4FH84zJj8wtHayipfPIJBo8nlWnKhy/UvCHpCv3RMLDE6OhEgAmcSARK8Z9LZpLkQgSojAO+u2oWEIoSTNcLYW2IKj0SR2zUMk2sakp8uFAXwIecYCDARDDIeQM8vj4f5kAWOKsrbpFiU/t9TPS2I9Mox3UxdJCsYoIWwIzPS9YpqR3TdCYWCzF1XyzSIddXvZ+8NBs/9XCjE7rJs13pUO+hGEQSh8tt/XR9bNebpnbBMKeaiopOEkoCsRSk1PypOSL8u1xSVeZDdZYibYR7IsDOOzYlEdDRPtlEX+cM/bG1tAQjpST7l5+BUn3M6HhEgAmc2gWrxmJ3ZlGl2ROAsJCBFkhRLH21puWApY+/Imea7PKYhYGJQhmXpWq+CS+sW0wyDCSxg4/Ld0jWKz8m0YpqxRWNiUcWzsvKBI98h33AqhhTrT8iSXTioimyoXKDm111q2I8lazXIl0qhKGNOpRgbGrL/n23zB8rGwfuMUrupcKesKfa5Ln312ua6Ft43fIsU7wj+lMxhFQTuXES3C9ua0S8crMblmhOEH1mHeEcNXm6pOivLZs/SYqLhnwIpfiHWuduNtmsaxwJDT7NlzQPr3nWj1hLoZhpEgAgQgdOTAAne0/O8UdREoJoJcAguKQWdDzc2fmUhY59YgEeTkbkNI5srIGqnwJo7iLJkx7JZNgViS4HoPZ5ADjKfY83wlw7hershnBm/rQ9eyIeyT49PVqxg6jqoLwhhqYlPlvjlsxD/7RCKQYhYWbJLhxAsYtGa5lGZinnAa8xYDKIXgpFL0ZjPYxVeSflFLPpsc2/vlx43rO+UObdNVTXO9XluXju1ka08PFARvWAjIEJl/CdycAh05c8R93LsW7If3/nKpqZFixRl1SzHeVM96NcIrrgUjVkQ8XF80ajFpudy+AQ+gvPA8Rj3HQP7G+Yc31XYA3XY7/j+6JYIEAEicDoSIMF7Op41ipkIVDEBKXaloHtLU9O3ZqvqrVPRniGmaSLk9aqa14sOB2htC1EVgNg9nkyyY5jLE0cOsx1DA+wc1BGY4vUrBQ4RjJqxkzT3hmdafRvyivLDQyX7br5uaHB86o91dmrLu7pQFezECF/EDeXKFMRuQeyatzU3X7gvnV5ZwxXh2I4iFV8BsZdxxwfRLoW7HBXBq6pIUHM24DjLFc7fW3AcZUPZ/LaJJw1VMRZiEdsvZraxlfuO3yI/I4/VvWIF7xga4quXL4cHopuvx/3XMZ/KlwqGdsZzoWJXrmM2vgfIrTKu8npbb3CpSzZ4vBc3wpqwiPPgZHyxaNN15pXZXIwmiNsEhO4xxD0Vj5V0mll4TikUWKZc5scxv7TDd2GWCvY7OtnKJ+kHESACROD0I0CC9/Q7ZxQxEahmAioEo91cX//mWkW5NeA4pltVNYhdrsuMaCQC366LCWRE5aX0VmR1d8IS8MWRERZGEjHCVfgZLGbgKvuAURYWWhA36cqyRpe2zO/R/2F9bMp9w8z+xYo9xx+5tKsLvgjG1q5YIe0ObOW6da8361vJinZAs65GzNCx8vPOBa2tsWWO85PFnF99UzTKavAk8rk85vKykVKJDUAkTkamWpWZXczBhGjXSyVeQoY0y0X4/9VF5v/jUPJ203GyG0riZ7bCXZaqWguCvpt/0jFt0s6c9Z41x44lmUxVY6zp6qrcjv+Q88Fc/o99QIrkuStkpnUFexdeB2fBukYzuU3Q4O+vj55Tx9mVUUV9U72iLCkpalBBrncqpCqy1HYMC9S8QVQbg+dYfumQ5yCG2BXM57eY10g2zdrhNA54kIVnKt/vONZR2zwguWwerTbxf2Iaj5luiQARIALVToAEb7WfIYqPCJxGBHBZnUkZ51GUv4LIkpfWuYLspyb9oVJoQUDKRWkcGywAjEN0RYtFtsDrZb9MpWyvYoiLNUWb7eLMcmn829kCWwgh2e64zEaPu6FW1z7oVl0ffGjRrOdzjP142LR/AXHYM47o1cTi+OvQfhzxKdCMEMijWVH5GsQj+1u/Ppe7/W8ZcKw/PVfRFk0Twm7SdaXOgz5rMiuKOMLIfiaRCe1BVrQG8T/R38dS6SSb7Q9wN2wCOcG87+TKk1c3194zSXd/Z+rx/o8+ni2szmC9m6ly87xI8Mr5UdeD/+aZtPJgU8PcGSH/ZRG3a4Ff144qTHly40j8XsxnSKxapfA1a8TaFUzZNdTJV49msp0K3NEf2scbYzMmafoFAU15a0hRzouq6qQ6xImvDKgR7LD/KBi2AWuIT1FUt6apPjCu8JdfPJDp5ZkMVtYJFsGXjr0Qvo+Xi+xP/V4UHC6zfYqbW46tXKtYPz+/KXrzef3Jp9ZC3gPbKfEgj58vuiUCRIAInCgCJHhPFEnaDxEgAlJM2qyuLoBU6SJTCKWES/pFcJFbCJu8Li43ZSw7KjONcqHaLGRMz4egegYiTBQ5s1VVzA/6OJSs+Z/Dqb4FPvfkuT6TTbc8rMXnFk0+z4JWj+vfQoL9/b0Xn3NXmvMfv/fxbY++WmZUXvpfvRxibU0lI2qPakam/m00uLRO12+sVZXOJkVZ9ABStlGhMGSm7YDLpUZQPaIiEqX3GIvW5AK1KLLRAnH+GPd/MDjA3gaR6BhFFLl1s354fv1chJo0/b1Zhb13R3vrpjtyxf23JzON9yazPAU+M6Kh84Zra3bPiYbds6JBFoONI+h2yb7EH2oMuPvOf8+1H4LYfQBfFiqinCGNKwX5iqbopA6v77yQwt8W1bSlYVWdUa9rmldma+HHRVtgMWRaPFkusf0lg21imnoBXjPxWVu+R4p2+cUDc6ncYj5cimBkqS/AF5DvlQrWTtSWuFhTOSbL95eKaDOnzwlr2vrfTG76m7cd7f8qdsWk8JWLCeV9GkSACBCB04UACd7T5UxRnETgNCEQ1HU3LudH0og3DTGbQmZ0GN7QCLK5bohEmd214N9V8djGawXIq+2WZR31+H5RyOcPP1Uof56rCheaal9RE1bOiQSe+PM9Rz/5dLbwnhaX/pY5QV9kpt9rt3hdoi3kD7WFgze5VX7TPVdduKPEnTXv+O3GX6/v7FTX4/gvZEZx6X9NF2PnB4M1VwVcb6pV1RuQDV0C0dhRp6kM1QvYbrQ23lQ0rMuZUODLVd2wXmjjWWlYMbgUvLKiATK7MTkPVJrYCvWul037zZqqzOQOd2ucrcsZYpLtsFlet1WnKUuujwXZJbEQe6ZY4rsNS93KdWdeLOJu9nuciM/rNETCPIbmxRwsagK+5pgne/+6m655x4of33ff385uuziqaheFVPUKxHpOncvlCyIGyFdmYBFdHt3qeotl1l8ssYO4PV42BgfLxnBJsHwpEl1QQj49h6x2FuchK6tM4DzIzPoLC9TknPAaKjKwtKJpT+H1Zq+LzXTr7OmC4L8ulu2LVUWd4ta/cu+MSRceylt/vrKvLy790+OWEmCmQQSIABGoegIkeKv+FFGAROC0ISCTtzym6/lSudw/xNgMbCLqONyPRWoarACTIbAeHR5iI8k4mwbh26Z72GF86s+Exb7gWGp7a80DocMDfY+m81/KCxa0IEaXRAJ/8uvzZvc9Ggj/CtfbF0TdrkiD28WFZak7h+Nia8+gMz0aErPgnTVsexUyo79iY/5emRl9X2OsY57PswSX/W8ISX+rS6uBiJQZVVZCPMOG6exDhYV7TZsZqA4szbxlfA5VFpgDIagiMyqzowL3OUSwtAMIvDYFj2+ERePhdFq1IDYNRREL/B4+Nezn3xtOFSeXTe85qHo72XTZbT4PuzgSlNlTlkB6O5HPCW9tRImhu1ljYz1zRcIg5TBvImkjKa72J1Pf+9b5HWlksicjw1rpUYxGb2hlbDrHS2VluFRmR/NFdqRYzPSWjJ6cafU2M55cqKnud7tc00KaNv2TnLvBnychaIcxz16cAy8y0zjSaBkyiGRZXSID4T6ChXrI+v5oSy77J7MMy9XgcfOLoyG2qndYzSmcLcXauzkB38qpUffC70enfAxi9wHgU1ZjV2Asvc80iAARIAJVTYAEb1WfHgqOCJx2BNRjR4+Wos3Nv+5VlL9rsG0H9gBFlh2Tf2w2I7P7z3297Gqvh5lmmSWxSKoHBRAmC0fX3K73HnTYe7fDBvCDTC7+63Qh/HAyK+IQkrGG+r9tDwfYDFgA6mEhgAUAjRMUli2V+eMHj6l3dx+wJ+cKTnPYvwGH4R+f0rCkye2+GJnRG5G9XNrsdqkyiym7i6HLm3OsZPBeZFyRFXWOQJ0bZcM4jMv7Xp8rnIFAzGAnSYjELERiBIu6uBS7ELkyu8sgEA28fgxlGva4PcVw0HloczK1DFnhOlNV7EVBn/Lxtobe2w73/c8zpdJbp3s8S2YVymyaFL/YGiF+m8slziA+PTPbmQudzkRdHZy3gqHCg+rO50Wt3xdhNZGIDSPtULmsJgol1guBeyhXYD3F0rGUYe33YFHcHE0z36W7G6f7fe0xRW2NqKoOk4hsHcyaMddj6Ko2CWEHMBc3vNLyNR3HzWVSrBZC2q25RBo2jD7bHupNJj8i9fGmUvk2ZM+tZeGAdmUsdPy7A/FoAkngEYjeeUHfjBq36/6fLZq5mm/dB607umjwlawk8jUaRIAIEIFqIUCCt1rOBMVBBM4MApCLjC0pFu85ruuf2qxpKDHGWRmX622IR2lh6OUKuw9ZVaaruIc+B6rD7oYgbLMdMcvnttFWeMlNtRH21too2wJRug/d2LjP48T8HuFHpYEY2hHXo9OZrirMKBZYayjA8o5Q1j67Y/jt6NP2LSxoQze0efUunblwBCn+CjKLWygpxwpFdjBfFLAB5NF8Qa1HR7TzNaV0kdfj3+Byuf8LijgOZTuMWfSbJvMjA9oOYa3g9mkIdRcWdrW7PPJ4Ti3yvB82CgeWeJQ1FwnH2pgr3pkXYqqlKs754UD7L89pv2hzJvf2z+3vCW8rlD41uVhe0V4o+6cHTV4DZTu/poZNASy5cExmm2XWWIpqR/7QdbEtneO7h+JqT6E0OFgs9wjb2TtZVfOdHr3pnIC/Lqqpi5GxjkgPrynniC8XR1DbeEhuyACXNN087vPrB5E51rB/acFI5PLsrvgIewusFws9LmZoLt6r6EJldsPjDbH7jtrKh29KJC6d5vUsbvRZ9gdbaj0Qw5/67Uj6L4Zse0ECC/kWRIJ8atC3+n8vmHdhr2l9EGK3fxWyvZTpPTP+B6ZZEIEzlQAJ3jP1zNK8iMApJjAmeuwvuVyz9HT6J6hBq8iuXVuROd0GK8Ah6LgZyJS+CzaAnyeTyv0QubABOIsDXqU9EmTfHUqW2soGbABeNsVys0nIhl6Iy+oXwQYQL2QVI+tikZZm1ohsqK++FqlKZCiRrVQGsHBsxiQ+wyjXxTT1/VJGl01LpMqGGCwUlCPIiiIzag+Xy1nddowWxHW17lJm+byeRvQ3dnMegOp0FlqWbDLBD0B7RoTgXmRFBaoyoDMZuzedYk8ND7PLA14WNw0pFBWOLnE36sqClKpu3T21+f6f50oPfjWZef+vk1lPCjr2Uk29ssnvfeYTb1r009ZoZK5fU1yoxsA98CcPpjKsgMVjQ+kMC/m8LIL8rgIBLyBIc4YpLMxhWNUPb+2P33tRJBC+LhqeNsWjXxxQtZawrIcGASsbYSRMy4kXS8pA2WBHy4Zy2DCLg5a9JW87T16B8r8eRb10q8ejG5iPhYx0ECXftmFejbBlTEP1jCDWpz1nmnySaSiTPK6rZ3rUrveYvs0PZ/Lzmn1utSHgq/vUlMYb7x1KLtyRK35nwLA+eAhZ8fOMsLW4Pnpls1vf8o1l53zgLzZsf1BqaiSTcUODCBABIlB9BEjwVt85oYiIwOlKQF4xR30D+/0ezifDB1u6vFx2v8cy2Q5UIXgWmdqNms7ium7XhMI7u5OJNqR/YzIjujTk55+c1ND7ucP9a58pli+d5taXTIUNYfqYDaAl4GP+RJwpznTmQpczARuALG8mF5F5YDnwuNLM7fPyvnRODBaL/HAmzw5nC1ambOTDKGAwDYaJt/h8ynSXHkETDEhlLv27Tty0Kt3GjpdNiEWDqV6/cdzjcQUxC8xELULscojeXRCUj0DKRR3UL4B4bFAcNiBs9suSLTpgUUBzh6veXRNkF0T94v5cie8tm/y3yZwTaKxvnRWN/N3UsI81Bf0MVgXmg5WATW5ifakcG4zHKzKxEX5ajkzswOAw6+3t4yNwQM8rlhsunT3l/RFdC8p4K1lceHgPlE3UAy7D2lBmhyFy+w1rd79lb0SFiN/AnrAV5+DIvar6DUWIKy9IpZynfT7+U1SbGETmeB54NWIh3oZCjs1ChvdNOhbsocLErwolDq+uOUdVZny0LjqjDl8QnssX1ZDPIxqj4av+7epLnmgO+d09mRyP5wpK/0iSPzgQt89rrGmK+d2//e7l553DH9q8c6ycGnl6T9f/gyluInAGEyDBewafXJoaETiVBOaOZ/eE6MbKK5aHJH0Mt3MUxWqB5eB9lsHfhlTqNmGlgx71G59RlB07C6WvonbtElmGbGkk0P4/C6ZffzyX/8Dnk4WZdaHAf+YCvlgSZbtSRpl7DYfVYX8mrAYaqhM4iskULOSCVRXa182Owhuwbs9hSFBRqOfcnOvWPfNjEa3NrYfg5VWhXZEVtZxefBa+WHa8ZCgHyoZ91DL3Ddjst/tN866Q11/wOc6D+ziPloSwMgCYQuZ3rtcrdqH72F2Fogr3LWwQqtKKbO9j2QJ/DML1Atyf4risFiz2+mB9VMlDXO6ybGVQU0QIvt6A16vUhMNKc02UeeFftiGuffgCoMeTbMehI+z5bA7NK4qoS1xkKIQm2y07Lr/XC5HL+5GpThbLXFZiOFYsK1i0Fu83rW1DlvXbPsN6HG+XIrdSJgwZVva/aHqc5vzP/dItAV4rYMdYAP3+JX9g5EFNy9e53dGhfD6wxTB5K9heGgqI/0nnh3+aLTRcis/INtDXRoPqQtgktrvRMKSuViyIBJY1IPO+pLWBuZChTiKWB3YfUh490uOc2xBVWnwyahpEgAgQgeolQIK3es8NRUYETisCY7VZ+WfQpeyvVHXwDsZuQpbx+vm2E5oHEdUuuJgCQXWeyuE84Lff0Va/56AjHv/kQLwdNWpjGZQpOE9R5g6HIo/eNrO9ODMSjMWQPfXpOjdhH4hDDI7gkn/fSIJNQqZShyXALuRZD0TjweN9zEhmRGfAZ80N+vyNLt2NhgtongD1a1jiIATugKxsAMF4qGQUei1r24Bt33HMsO4D5L3YKoIxg4xroLn5zViU9k1YMJZioVal01oYb5gMK8bBRKL8QL7kLsOKsRST6IyF+I/i6eSXh5L9WKzW0eFFrWCfm7UFvGK21807mMWzibjma6hhNWj44G9sQNu2IFNQKzeCTK+sEjHQN8hK+Rzzw/ahhAJOCVlnxKoMZXLsuKzEkC85sCzsixvWgz2m9VjRtp9EOEgNjw6ZVn8Hct3yUQfyxbug0duF+MftnH+uHyp/Orq8LbZt5UeFrL1V17/40QULv8uefPKaLWXrrpmmZaIShP7t6c1br+k+sjXFxGd7IXSHsIIOVRn4ZSjegFJxPNLSYNdHUT4t6Fc0MInAO33TYi9y3Yp99659X9jQN7JzFXy8qB9M2d3R00I/iQARqDICJHir7IRQOETgNCeAPCpj/2nbD+FGbq1HhXj3Bse5aQoWPc2zhTYbl+6nMZfZpKqzp3tds38yrVk8VTT4zrLBNziKmFETDbYGfMEwMqGN4RCLweOqozRWEVnFI4kkO9Dbx/px62AB2lBvP8sOjzAvROJcXVNCrY1BA6UNMrjsL8t29eFS/REsfEM2t6fHNJ/tN53fDJvmI4jrGDYYBUa3i1GoqwtiXD7s6+t7bjFjF/c3N1894jjLE4wtwPuOWZw/2VBX9+jeY8c+vT5XvNVEMzjYMZRbm2u1axuiXX/VfeTzzxXKtzbn9DfPwmX/WbAwTENliRhErR+xaG5YMLxoYAHBK9src1gx3PDS+oMBloinWF82z3qzBeVYJe5iz1DZ2jRoGPcMl81nENd+bBa2Ssyfh7hcjw0xgyZz1r20EQQ6d9h/j6YeOyBOv3kBE7EkGnicqygN8zj75hPd22e+ybb/Br2Q79tedl8zybLNGX7v1Yk3zft07MmdDz9mi+9AWE8fhn1jLkqs1cC77NTG1IaOWcwtW0PDBuFBiTne38/ePm86WzS54YErvnsHn7tiBYoQIxIaRIAIEIEqJECCtwpPCoVEBE53AmhBW8k4/pKxHnhhvwzR+JWEzd78vF1+d7NpXjPTsBo6UDWg3XJbbciEvinsVy7Rw3wYVoBEISc8aNRQG43w5vo65o9FKjgCWDiGDrnsucPH2Y7nnmchZBqDeNzm9QjZSzddNtTD6SzrLRSVQzmI3ZLR3WOYm3oN805kRR/DTqRDoSIYIRL5csQoRS4UusBtRUziZTmULRCIrK/vbtyX2wtjBPc4Yx+FzYE9ni3eanBul1QliOYYf3bX0nn56zfuvCJlmueiHNi6LUVzemMiLRbVRZW3wW+cRVmwELZKLV/4mh1UU5ALzyw4itceHxS7UzkLNorNA4XSFw3DWI9DvRCvDOCSUVEunRnOmlFxjmm84hBrMbd3O85aTOq5J4Xyg0HGLhqyuL0IC+Pmqdonn26tX/rOTP6/H8gUFkPs1jTCI/1Y3vgJTtqCgmVdtsd27e8djLuQCWdzkLEOIUPuhg9YyHrB8GPjNDFvNuv48gUUXTY/hyje+oqR0JNEgAgQgSohQIK3Sk4EhUEEziQCL8o4Qh8y1AKQgtJ+EPruwWO23XDMsq97tlz+QJteXiYF1UwsTpsCH2wjsrot8OsKZBU97VOZDzYCVocCYBBqCsqMBSC8GlCGzIyidi0uvacKJXt/Ii1Ld6mHsnl7sGRsHDDMh2BdeBDH3YjNllxlEO9EHEO42wWxiMdSLCKmVxzyNfkRZex2/LF8TgpkzoX4KDqXPfdUpvAtVKOwTVgc5sdCt31n6dzYzRt3fTgciRyt9XnaQxrHwrValoVw7IEVA03LWBS2DGRgWQbd2uK4n0Z9tGtmTlEWjyQ2fn5z9yWo54vquS+NFw9fLsrx1KuPMXuJ/Pu+vxNZ6i4h/i0h2CewcI0lFaU8X3Mv+8uayKxMIiUezeS1Fr/HXhgOzP/hknkfWmtwkBXudD7vPN7Tzy1ko6+f3MZQqxh5eQt96IAF5eVMW6hFvBO2kXNWdHYEUJ4Mdmyq1PDqZ4VeIQJEYCIJkOCdSPp0bCJw5hOABmKWFEK4lQJS3hmEYPpO2mHfSZvFZbvLxgcac8Ubp6TzNbNCPjbV7+UtqOYwCe91INBkFS60O0N+U3Y4g3cX2dIDqFCw5WivGDIddchyzL5c8Yv9+fwv8JFueQw5xkWjvMiO47/8sn/lPb/nh4y7IpZf6T14UUUr4G8XsN9n0rnbIXoZ7A3GgtrIh9Zdujj5/h1HD6IE2ZtnwAIgF471JDMIApUhkNUNY6GYXEFXQLY3BVvDcVSUWBbwQmUGi5+HlP7W4sX6rVu2WC/60vBKIbyW56SgV7rAHwL9k3Eh7r+/ZHx/UIhm1Bo25qtKzW2NNewRWCi2G7aadXsFPLrfez9sGCEX6kKgfNvhVDszIdKlsO2DdaQJE9fxxcOEdzqB2FM4DyijFjILjhvHyr2WoOg9RIAIEIGJIECCdyKo0zGJwNlH4AUBiTsV8YtbB1J2g23bG9Lh6PEDKv/C0XTRqsuXtfciowiLAksmU6xGFhHDAi+BmrtpiEQT2dIhVReP5S3m83rXl3yeT/cPjmyUO0Uqli+HGIXIq1z6PwGiEXt9xWFL7QfR+50SXt6czN5u4rjotGbNrY186lPtTUM/PBZn3fE0auaiuwU2C1nRDOIPQDBWUt6oFZw2DDGcyXPF72YjZSsK77B+y+bN1q0Qvq941Nf/ZCU7jZ2p3LYfMBg7b2NJ/GDIEVcOQ9AvRKnkN4UCyhC8xFl80ZiKxXXIsosY7nvwRWNWTYQdg03kCBYGyi51WfimZVk10yizdC7nJNHMI14yDyV9DWkhdkMjn7C4X/9M6RNEgAgQgd9DgATv74FDLxEBInBSCEgxZ0uBiuHCZgfd+iY/vKHtTbXKpJC/UoHh2HACtWkZK6Cdrxde3RLEYjyVZv2oXtDs94rbF81UXEbptsse27plRUeHa113t4V9SoH3alYFebwTOUyEB9OC850y5vBcMvtViN5giXNraXNt/U1T69hPjozw/amcbaKeQdG0FZkR9WPRFyS8g0YQCtoF80A6IzRPPVaf2TNQQC0I0Qj3wQm1BkjeMssufdX9iPcq1O9dnXKcVUMQ4QlUZJiFzm+hgUG0OQ6wlimTeA062WnI8pooaeaBQEfFYrZ3aIQlUD5NngvUWmboxuxkDFPFQru7uro2WutWrpT7f9Ws+IkET/siAkSACLxeAiR4Xy8xej8RIAInkgA0Iiy6qrozUyrnUbfWW+Nxi95UFjLMYQayuSOZLEN3MmbjvhS/R1HNQFZlmIHGCY/F41Pw+edWzJ1rQ/BKsXuqhxSSGtKf3xeOsnPbSOr+efXRKCZVbg96XEtqAvzhwbRqQFiOFA0WceuyJbLQFEUpoCTYvr6B7qv97gVC1DE063DmInpZkuEkDSlGZRdjfJEQq5O2s/HRTOGnfaYV6eRRMRtKm/X0seDS85iOjnYC3l1XNsPCyE63QqjvT2bZ7sERFoYAFlDnhiP0o+nc0OaR3NflLuHhnQj+JwkV7ZYIEIEzjUDFU3emTYrmQwSIwGlDQOovtaenpxci7Fe4PK4cSuXMXrQDPprKskMjSbYfImsfOpAdGE6II+mscyieZioWVMkMse7wc3Hj1A2hcuzEDZlRRijWJnTvvbLO605E3C63BfPuBXX+P8HiulvixfKuQ6mMszeetgfTqCWcSG/fOpS85ODgyP/K7K+N1Xw6hH9M1+UXgJM5pCiVdg9VWNb9bY0NJdhD+P/2DrHf9A2ztM/PPAE/KjF4mEBHNrmp8Ey7YHkIQ6yDP9sxkhI7Ehl100AiuTORftuze/bEV62q2FTkuaRBBIgAEahKApThrcrTQkERgbOKgBRhCrK8f23a9pzeYul8AxnEvGlytAZWwsg0aqjDCzXFB9FxrBe1dxdisZW8fh7RFdhSq2KYF1/SqXV1dT1bF3Rfhg5pHygZ1hO3Pbb111KJY4K/nhQJHrplztTgZc11vGdo5M6rn9j+jM+lXVIybQbxabsc4XMHRCtaSuzGR+THTpaAlIkOO1Zff1Gdz9vQ2T7FGc5mlSgy5pPbp7FsvsAiqB6hICuNlDozZPk0WBjQwY5tHUo5aHUMH4ro7s4UbpI1i1fh3K1ZU7GSYLc0iAARIALVSYAEb3WeF4qKCJxNBCrCDuIp3traehWqdP2kP1+6OgMh6M0WUaWBCzcWUOmcD/Xmiqu1kcQtRn10oSxLhiZqCyWo4fr6kyUOX/N5gNitVEX4zGPbtuNDn5QflGLwHmRTtzCWdoSz1yPE+T5o2bhpd+JltWyLZ8sQk+h47Hg590QUdyue371uVPDKXZyMUWHlUpS0tFXES2XXea2NbGpNmOWxUPDIwDBrQcY5gHbNFsrAjSTTLAHhmy4Z4i1Bt/jTaS1qPJna/+Y9B5779i2oKHH7lpOdlT4ZDGifRIAInGUEyNJwlp1wmi4RqFICUoQpsDYk+vv7r4E14IqRdOaJsFF2LkapssuD7oK3kF1+6Pjxb6EmbyKPTmpYScX8nNfI+WDRWkXEVcHcnFWYx2OdndraFSvUNUjuvnV0IZdl2E4aJbxQgtdmoUqtNVmNl6klzMNG62QNCdVaRZGiGRM6qTOpsPJ4PAdMx4kPw597JJ1zDo2k2HFUxTg0NMx2H+tl3UeOsb3He9mR4TgbyOXF8XSORxybR9DWWeXKnLc2Md8tt2+R8fKTGi3tnAgQASJwAghQhvcEQKRdEAEicEIISGuDFE8aRO9DuFWmR733f6S1zlRtyxsr5erDaKSAOmaegmlWGk/4kf0d+wxuqmNIkbumq0vO5SWjYNuuEhbbSaHuEqOd6PAGCF6Z4bUrXeQghE9FtlQyU48ePVpqbGz8baJsfrgvV7A8GnLpWLdWwuLABCwNXpQf41g4WDZMJ2fZyvojvcWZtuVBppf5FN48y9sQ4GywgJ3hUzSIABEgAtVNgDK81X1+KDoicLYRkGLMGWtNbMaRyR3Ml9DVy1G8mktfB7kIp8NWXIqHnQHeV8eZe22YRfiaNagSW52ZRnTCqOhBdCR7TgpeWDaQ0RUdX2tvd+POobLtlMumrboxAQjLGfKE1w11nmwNWfly4XK5/h6+6aHj+ZJ+IJk1Zab3cCIj5MLAw4kUtrRzOFtQnh9Ob9h4fOCj6BjHC1DAIVXzd3j8M2Ws61aMNhSR92kQASJABKqVAGV4q/XMUFxE4OwlIH5ZWeeF9VPoMmFaBve6vVjYZU8DkkdQ1SCN51GmzGYextWIGoMzIFG1tGQ7YxlcWYiRIuJ2ELdfYao/mZQ1iI9DzBcN03LrqEOsOU77KZqI/GKhHDt2rL+1oeFtZYv9GKJ3Zhax1cLi4MPCNChiuWwNCWln62DRusk0zfJIySyjnbM2KeRT/Jq6EC8/OSrOu05R2HQYIkAEiMAfR4AyvH8cN/oUESACJ4+AFGNy7MPl9Uy+bGluaEa/pkjBK4VjJgfhCBuA7WHMM9XFWuTzq6s0wytjkwML1DwF2AWkhxfzcc6d4pJ1cfUyZCWyrEzFLbpOnApLQyUe/JB6VukZHNzkK5WWAOea4UJp0+54utCDEmqNnPEOzoYefn731Tv27z+E9xYSppFOFEto9MxZRFVmje+IbokAESAC1U6ABG+1nyGKjwicvQTUnOWoJSz0UmBYCI32TJCCd0sOvgYUGHB8nHtbdK1VIppbpYIXuc+KgEc1hi0FxA1PhuVlLLAjw6UloFwSwiXtGSrKfvm4KkXwqRwV0XsomUzDN716oL9/aTYZ/8YlHoV/tr2Ff6Q56vnn2c0yJvlvRTpj2QcSyADLChnw8XbIQJcvXy73QYMIEAEiUNUESPBW9emh4IjAWUlAfH5UYGXQzWtfEQvU5EIv23YukDSgvlx5ZETLhsV0SMlWTS1VO6WKp4GxkQKErWFbAl3VNNRjqEXcJurb9kpfrxSRecuSNgG2vKtrPMt9KqYmBasMUdqIOSpgPIOycHwoXzQ4V8KNPr8UtvI9TtqyDw6j2530Inu4Mg1eay9bs0Y2shibIt5FgwgQASJQhQRI8FbhSaGQiAARqBCwsKArn0eGV7pJw7o6KgJtuydv2aUiKgd4obNwf7F8d11n1YsuNyo1MBNiUcNM3AqXmVPhCHHcwPNyjnAqV8qsyeflnE7hkMcz8Q+CvN01WCgZyXxRjWgqj6r6uTIOFbUY8pazcRhVGgqGaQcUpWlxY2wKlK78DAleCYkGESACVUuABG/VnhoKjAgQgYJjpyo2AAhCbotmZBTlQq8jEI4lVGpQdWitgMJjVU5KjKnBvSXLzhZNW4P3mJmOvVDGbTpCLyPDi/kxH7oMT+BcxkX28YRhDQwXiqqGyFH6rfKFQs4BfuPuQTSnwJcQO6yqrik+11QZ7zoSvBN42ujQRIAIvBYCJHhfCyV6DxEgAqeUwPpRSwNDRYataDHMoK8YxG3bm2IxqRVZwXEcaXVQR729Y3rylIb4eg42LiTTaDJhlNEyWYcJIMAVaWlgJdTklbV4ZYbXjZ+vZ8cn+L3jVpJCyjR3D6G7mlxMF1CVBZ2ojTw2iUPDhlmIF8uajyvwVavzZQx1nSe9jNoJnirtjggQgbONAAnes+2M03yJwGlEIG8LVxYWAFnZwM8VI2Hb8m9W3hRsJywNzIHvFR3MKt7e5V0TKhZ/L9UxsajnHceFmrsQ6li1xrkhP+Qw5znMgVnYkPWdidIHQTw9Ib7Y9WNfNFK2vXWwUJSZXCykU6bN8ftj8BrLcEeSlj2cKJUVucgOjTIqtXiXV0FrZxkcDSJABIjAqxEgwftqZOh5IkAEJoxA15iHFZUNtkrBa45WNgjG3OpsBCVQ2UDJQzgyZHjDSrUneCtpWxlkAfaFnhLaCMtsLuwNiyRg3B2BVxmC12YuWJWvQKk1+fxEjHHuYLu1H9UY8mgyEVG14MUtsYqwRUz5tI0uxLA1SIHuEnwR/A46W7tuIjPTE4GKjkkEiMBpRoAE72l2wihcInC2EBiTsUksSkNXNdPxorJBxKMF5PwLQiTyEI4y8wvR24IULyp9nfKFXq/1VIjVox7XIgTvcQjdSkUGj+B1q5BRjQuB1r0WPLwO8whhfp2xOOYu06nqaz3ACXzfuN94R3/JMOKFkhJCy+EwrzSZqKxMy9jO7iEIXumtRmmyyfMYku98YjLSJ3DetCsiQATOcAIkeM/wE0zTIwKnKYExFwDLZ23bKhm2IleruW0RlfNx0KZXZnilFQDm0lb04/VJkYgPcSkicStLOlTd37eiED5pxeBI60a4kluDqcz1ehv9aOUAYW+7GPc9Wl//kT+vqwtgPtYEzAGHrIzjcRML12STCblwTRtduCaT6VlhbxyE1QEeaiuiKIHLasNTxz6DV2kQASJABKqTQNX9g1CdmCgqIkAETjGB8UzjvqLtZJAB1bxSxXKxCHHIdmS+vIMmDhCPUIbGXIjDx7Cwaj3esgYiEsrLxgbDQ3WUy+oei8MQzgFUaZB2AKck7Nkfa2j41acikY9e4w8Ij8OVBt3jXe7xfHuNru98vrl55dgcTuXfafH50S8KWLhmdQ8Vxxau8d8tXINWPzBgmk6mbIggV/RazqbJ3w2q1CAp0CACRKBaCZzKP6TVyoDiIgJEoAoJjKUa1ZztqEX4eGXKNspVaF0msHjtsWHp7UX21w3fq+HT2y+F6JXbZ3Xlw/t0/S7h850DwSizvhP+d25oXPAyvr8IGwbaIjvI5rbOV5S3T+act7lcvC0Y5MFolPFg0K71+SbPUdX/3d7c/IFTLXrXj/FK2c7zA9LHa1hCLlxr97Ma+Q0C49CIZedRrUGHzYR5uXKefBIl4xAq42tHs+tVmWGXcdIgAkTg7CQw4f8QnJ3YadZEgAj8AQLjmcaMIca6rSG1qDvODNbQ4Mcl/ys6PF6RM2wuuKr+Sbh27aGWlr/7W6/3rzpV/bsz3O7rbOE8+AwaseE4FavDHzjeKXm5LOyA7BInvccewYQfmeiQprGIH/dqUaWspQURt6pOTQ2K9XrEJMa+el9jY92pnEPXGImMsPf2Vby6phlV1eBloSicI5WRxsK1gaRhoGuG9JIo50AHS7ErHQ9i5Wh2vZJhH3s/3RABIkAEJpwACd4JPwUUABEgAr+HgFV00G0N1gW0FhZpxi+8RVGeucbj+cS1/gBr1V1qs9vDp/p806bq+hc/EYl8bYEP5geVm1uEqH9U01rGRNiE/q2DiJS5UW4JfhAWDWlp4C4huIZsqEvXmebzIX0dZay+nkHQM6W2VmU+nx3R9XAz5x/GHOQ4VYvYpCUEhhD2XL9hWqi5y9FkQqlRtXPxLHK6rDjo2JuxqA0dnx3U4mVT8Jy0oBiYpPYvLu3zB9z6nUL3yffLiU8oexkDDSJABIgA/t7SIAJEgAhUL4GCYKmCFLxQiV6mTJ4FyRVyHDvmdquxAIo2SLHIucMKBdGUzbIM5+KZXEb/qWOv/i/L2igvscus40TOEKW7lC2IAZ3KNo+gIgNaCUvBy8pQgzlFYY6qMgW3DLcVSQkRzHSd40mBx5ch9i9iU8bmgrzqSa1IIcW5HGgyYSWHSuU6zetGRQblfDwns+Wsze0dMZHUPYbVhEJ1zdrV2PxkXNjP9OWyF16m8GV1aErhmNbsDSZbiLejt0ZlVuP7xVM0iAARIAKnlgAJ3lPLm45GBIjAaySwfjQz6FiMPZez7BuQTRRh4Yigw4TH7VGDUujGYqMbMpAMYheyyoYY1hJuz7P/lcl9AQpLXmafULG7CjGsQQfhq2KxUJuifHiqbYuyha5qQmHzZHClEhvG1lQojApeKXzlfXRkQ6c5jioVmCRz4a2GRIdb9mgn02SjDdyX4vdEj4qdZA0KMiQcsXcQgtdC2+MgU+bjQOoHmpv/+0LOb56LY8M/rTa43MyrKBchBX9RTgp3o2Qly2Xtadua+gRqCi9DRlgK3hMdJO2PCBABIvB6CJDgfT206L1EgAicMgI5iCQplLhliQSyorA0MEhcbmKRlw19y10uVN9F+d1gEHIQ96VQzGY5y+XYAt1V+8CCBh97flAmUaXaOhnC8A+yGBO7zvmNjefNVdWfdnI+czoy0A1I58ImwKZiDxa8sEPDw2wA9oD6YrFy/V9A8JYKBeUQ5l7v0hYnmuu391nWPT1F846rstmNl3YxfA+opHk5VospKJEgxa+c6gkZ60e/bFgpx9nWXzbeVAD/qMLrL6uv//50RfnTVnz5aNB1pQ1fOlT5xUNVbVYui0Aux7OcK3cWS31rLfvvlzKWlgwmiv8JgUE7IQJE4IwgQIL3jDiNNAkicMYR4LAAmBBKbJ7f7/EoKsuaNkf+kLXgojp8vSyDW59sPGGaUhYjVYqGDqPiSpSw9Otng6zuSsaOSDLy+fWjm7Pm1IlfGT6sudFwK+drZ3I+NYzGEjGXS68Zs2JIhapB5Dan06xvZITtLZdZEHORIjiLrQ8qf5pLFdFAaHaU89lTTPO2w7Hg82nb/PG+QvE+nsh1Q+wCAua4apXCutfwEyF+u+QOMVBzd0sfau7mDNPWudI0S+F/6kOjDC8y6iGPh6mRyGiGHfYShnmwoSE7mEwqy8KRXR/K933/tyR2KxzpBxEgAhNPgATvxJ8DioAIEIGXEuDICvIfR6OhK73eLy/l/AOzIVK9KMcQ09ysHfI1A5GbSqVYDqIrIMWu9LxCcHFkRjPIRhY5029zzB/9eTj4g83l7CO8xI7jEJUsL0QmX9/J1PWwBJxM8duBqLB/Y5bL9dcNqjrV7TimF2L3BSuGLEGGCg0sk6kI9qZEgh2CLWMvMr86hC5EO9/u2Oq64SRb6vWwOQG/3e73OpMC/gVImf5bU6n0L3tqajaNmMYPt2RSv+Vr1vSMYxQrVqjroHxXIvOL56Su/n2Dr1ixQlkx9o6V6yptgh3ky5nqsB29pmXFDVOdpyq8xbFFXtUUjmy6KpnL7G44PHoL/qJcVnkuJ+rK5cvubmqazfv79+DglOH9ffTpNSJABE4JARK8pwQzHYQIEIHXSmAVhB6EIrvK7V43V1Eub0NGsRmXz1thX+Bud2U3MXhew7Au9A4MsDyyomGIXgERXIToOg7Bi2641txo5BLUzbpkkhHMXG7b6+Omeec+o/AIT5aOsRdZAtZ3dkL8dp1I8QsFztRuiF0ISdeOJ5+8VXWwyg6tkQWEYsWKIYViKDRqxZCxQ+hybCHE/iQ2uDaUQexkh+188uFCyXrINN83NV88f5HHpaMKBesIBqxJAR+fHtGXTRdi2aRoONPZbK0/Xiz/6vGDx+/g69bB0Dw6Huvs1IbruwTEbyUTPP48bvlaCF0IXHud3F70wi2LF2vf3rJFVl7YNuKIzGDJiLWgTNpirvIRmWHHBusCCvNCT8ssOx6PFn/DJxTF0RUFPULEudjlHuxX8qBBBIgAEZhQAiR4JxQ/HZwIEIEXE0CWUV0DCTWtsfH9rZp2ecBxDL+uuyIQiFyW7ZLZRGRABbKiKm7rkeXdFo8zBZlSDcKrCM/rEBZ6Hc0XtAPlEpvl8xn1fn9ohq5dN00417UaRuH5qLkt6Vh3oMLAfXwg0c26ul7ww0rxO9wFcTiaDYaKe+WBF/hqvHdufb1Y0dEh1rP1ynB3PUTlOikqsRyNWarO5q67445/n1lbW1+GOMxBxebxYgHFa1+wYkjfsRSNeF1uUjY+gPljsdrRkhB/v2Vk5Ocy03rMtP8L24JnTOum1nzxmhmZbMdCCN95AT+bEw4YzZhjh0u/bpptXzenLvxPK03rob58/ufX7zn+2KXj80O9h/XLO8fFvZxbJd7FTU2+K+dO7ZgcCgQNy8l+7O6Hn7t9yxbzdrzho42N/76A8xBsFk7BcpRL8IVDVpdIQpTHIdBTOC8RnAOGLyAMFgyez6M6g4nycWgOwlgzdsHqSPBKDDSIABGYYAIkeCf4BNDhiQAReCmBpsWLfe7e3k97BAayohqELZowjGZEa2oq9gWZJRUQWR5YGNzI8G6FAMP7LQPv3+c4Q/+Tzj6AlrdvnpXNN8/zuNhsiMOZfp/R4PN6Z4VDyyAtl00uGV/cWRvbmLaNO3rzEL/HBl4QvzKiF2VGoUh/ZwtYsYKpXGZLx4TkWPTyPQwZXTX99NPhB3t6Phnk+t++v7nWdb+tih5UXJiDVshDmEt/qcwCEOweKXKlpWHMipG3LMfAPjqt8q6jhfLHfpjLocjBmPJGllbv6nq+bNt/c9BmnzloWkueKhofaE3nrp6VSLWeG/CxueEgmxkOGY3BUGurrn2w1TQ+uKu+rjteNu47kCvewfnBDYyNinu53/b29tClDeG/nh4OfLAl6JvU4PcwDQL8oQ/csCvvOP/6vYeeXdahqrfOxmmYghJprTgH4xn2OrCPQtweHxpicDywCOaBxYXMwPmIG4ZyxHF4XIiN8jjDL2InH9MgAkSACEwEARK8E0GdjkkEiMD/ISCzu7j8bbMtW3iosbFVluSSmVEshWImNtRjGM2EyqwohGNFLEoxjKc3oWsZLqNrqCqQynP7BlQ0eLoP+eDnHdH5UNm4oSVXuHy6K9M6D17Y2QEvmxHwW81+L5/sDi5TOVvWXCp/cXtDbGPStB48VMw99KG9faiE8DtxOC5+16Ekgrz8/7Z58xoWNEauqfe5L/HBbVHjdT+dKphPfWjdOiRo2T8vqQnfektzrdNXMu1gIqe+HVbat2kKG7AUdqAE/ysWqTVCqKO2bWXRXRFC8ahl8ZCmOqs9/nNSHtd9nw35ukZM61ebcsWHeFdXj5w+BsfiNEdbs2YDypVt2G3bod1lY/njueI7WxOZG+YG4sEFoQDriITs6ZGQUx+NdExSlY6WUum2LU01G1Jl447DqdRv/mmwkDon4r17Ziy8ZHLIxxoDflEf8AmvpvLWkH/uj7bt/ZEPtY9DimJjgZoakXaS8Qw7YpYWDA2L7Boxj+3IsLvh54VHmZVM05Gd5LYLMdLlcu2SmXDEXPkyUImefhABIkAEJogACd4JAk+HJQJE4HcEZGZUCskb589e5PW6v7SxZ9iDKgwioyg87gg2gsvlAYhCBQKrskAN2UUOwWhAZKH1sGPZpvM047/aIcQn2PDwgNyzxlka/tK7U4zdnYIw3GU7yx8tlm9szuQum+5OTZqLrO+soJ/NCPqtFvhh27zBZdPQNGFSObJ6c1397oxp3TtYLt7xnm0HXyR+17HrFs77UEcs9C+zaoL1rfh8VDZlUNW3FOEh3vBn77rv5+s3znGZlv2NY4O80bbVr4UDbKHuETpkX7ZU4r1wxqKqV6X2rlxwp4wKRdGPrGgJi74KXrfd6vX523Xt2mnCvrY1bGQvt607+63yL79/PL4ei9Mycn6r4JZdLURW5fxulG2TW+vzpfLlj6Zy72kbTlw6P+DX50dDbHYsbE6JhFhrKLhsssKWtcTC/9SpDiXdPm9DUFetACosQBjz1poo99fE2IEjx5yfPLdXTIMFA0N9xQw7rA1YoMa8yE67cLsN2V0/zoONLHtcUVzPO86au44eTeELzIQ3/ZCsaBABIkAESPDS7wARIAITSmBc7C6ePv2Gpmjwf2ZGQ8ENvcPiGPwMrZbJmuB5PVhCJhde0QaIK5e0AeCSehmL1o5DZLo0lX3HrSoly5o3aJn/cag2/MjmvPHAp4vF8awoxC/PMMe+OyfY3fts27evbJz7eL50Y2Mqe+0Ut2tOB7K+s0IQv+Gg2RYMsJaYf840TZmTKBu3ba6r250wjfsts/Trf40bC5oC3m9MwnujXo9VA8HbEAzwmnBA4Nq+ev/OvddszRRY/0iK3QqheytsAF7UDE5xxvcX8uwOtFgraEzJWCYfgnMYS9eEBpGI2sL6bsdJ3JtO90/LqnOxOI3NhT92dsBnNPn9wRZdu2mSbd00LRju+YRl/jZeKv3wqmPDm9agSrE8eatgefinxx/vwdq4HwwacjM7tmaL17Qk0jdO7h9eNj8SZB0Qv/PqYrbqcbsSTGmYqSlo8KZqXsQYQpk0/5TJ7DdbdrD3ffV7Sn0ozOYgawsfMYNNBPXhuGwuIRekjWbXZYUGeR6QYZfLCBG7ULDQLg+xi0WD37lrYOC/ZMYeXuiXL5Sb0N81OjgRIAJnLwESvGfvuaeZE4FqIKB0YNFXXWtre13A/f06ryfY5PNYU2qjWv+RXnY+6mJdBr/u7hJnO5AJTSKj6Ifokn7RArK+/ehaZghb8asuuzEYmDtVVedOduz3zCwbmbc5sjKDcedzhdxDH0+8RPwWMPENuPS+4aBtf+ZgqbzkyVzhxqZ4+tpJXtecORCxs+CHnR4O2m2RgKiNhudMd+lz9mVynwxmB1nEpQuPpoqQz6vV4DJ/XWsz2zM4wv7mJ3ezR7sP2peHfMrUoI/XoptayIDXQlhOj2OWbzNtq8cXCE7CwYehTGMQkRC8Cio36MO2HU9zfuUuw3pulyqu6zKt6yYVSpfPyWRbz/V5WAcWpyEbbTSEQ60tqnpz0jRu3l1T04353XcsX/7xe+HvHT+Z0n5x+eOPd6MzXTeqNnwZ27JnU7m3twwm3jqzZ6DdcntYxut3pkGjwjaCCNTR8miwLXzx3seQhS6zpnoP68OXiTbEOATuqHDBAsiqK7LBhxS7MsMO/hZeQ86dPQXrAtom7yhw/tXdAwPfl1YGeBnIyjB+UuiWCBCBCSdAgnfCTwEFQATOTgIys9uxa5e6Zs0aY0pby2eiHlcUzlwzGAzob5/cxOwjR9mb0UUtgMv/vFhgsqBrGsJLenlhA3Ckx/cAtruTqfx0hfsXIFM5C+JwZsBXrg8GQq26ep3pONdNNiKZy+qt9UnLvHNPovDQLYnEC5lf2AEsCLMNSLxuOGLbnz5SKi99OpO7onEkeXmLx710TiigzUB29JzaiLUxbyhJ+FpdKqrQqgrXIf58yMKK2lr2wX//Hju0az97d1OteihfZL9BhhfNJQT3cgOmBfd8Q/z6mZD/85Ns+/MHHeetfZzXQOwyZHfjAraLI5x/MdvXt08KRcW274Rn4c6dth3cadk3PlYy3tGWyS+f53WHFmBxWkfQb08LBkR9JNwxWeXw55Y/ua1e+o+NO3bkMvfAe7xn/DdKQPxqjz++AYvdNhzKFW5LKdpfLoz6/mFpfTTU5PcIFfySsIokiyVWn82xGyc3s+nDQ2y3abCLjBK7GhncQ8hQH5LWBWTY62RGHc/JBYPlbFYk8DgOp0Yv5zcc7evrwnFxtipiFzc0iAARIALVQ4AEb/WcC4qECJw1BFatYsqaNZUSXvbCWW3NhqFeI+RlcU3V9vQNseDzu9lVsTBzw67Qnc2zn2Dtk4OcIRalybJYlbwkMqJsvxD/8XSp/B9PK8oVdQYWp+Xzy2enMyG5OG2W3yf9uajM4AvNRMkuG2XJ2qLlzI4p9Y8lDeOuXZnCxj+TlRnGBsSvLcWvZYsNxwvl1djmPZPIXIYKBNfMrwld7olEeFskxKWiKyJ7W0SrYwXZThWdyBa4YLmALeJX8bTsCMd0XO7/Ofy652qa+wJd31UjxL8rx44dOs7YB1paW2M5y5qfwnFTmraD9fQkZAiroOPHsqIq7ADslxCSMPv+z4jc4M99zjCvbswXPzA5mV4yH9nl+VicNicSMCcFA4r0507B4ruGYugfNzc1PDBSMu7cFk+/ZLHbu89b8Pl5ddGPzYz6g7VerxJ06yhT7LACRKuVz7Ef/fwOrBfczuoh6r+kc3Y+5hPEXBWjzJ4Af1h0GawUlU5wHAvTUP/YzgqhHRPiKYjdx+RbYWHAGkCyMcjzSYMIEIHqIkCCt7rOB0VDBM54AjKzK8XuOxbOf/Os2tBftQS8F/lUJVbv9/E9I0mW3byTvRViK+R2sYPpLPuXssV+q7pYrFBEgwlFBFDtAN7ToazjfHZbf/934SqVwu0Hww77wbDltG4z7ctriuWLUbLr8naPq3WuH/5cWAJmhAJGo98XmurSr5/MxPX1kZK9samukhk9nCne92cHjr0gftGQQUVN3Z2PrXr/gV+t3zO3zqM7M2MhNebBYi3I7SxErg2h2H34GPvFD9ex+5Dd7YPQxRo7ZEAVYVioB8x58oOp1NcQ3pewyYpjMlTeMypwu3C/MqTXFYu7xJrfWQDGm0BAOzIFQlh8gbEezPE7A4bczI5nC+VrWjK5m6aOuBcsCI6WJJsVDRnNQb8+xRu+rlWI6xqj4eyyKU13j5TL6/4zbrx1ciT4EWR1WQjWhXosYqsFkwaI5gRKpH318c3syS272NXwDq9ARrsRgr1kWGx7LmfdZ1hqVtG5URZsCPOWmWlkhh3YGUQPvnTsQgkzGSfELqq1kditnFT6QQSIQNURIMFbdaeEAiICZy6B8QVqV5/TcfO0WPD2OTVh1gAbwiQsqLIgJA89+zy7Lgwhhuzs0VSG/RTp1Jv+9B3sXRBaGw4ehQgzUL+WGc8ks1ds27Zr+9rRhVHSK6pCGDpSGCIr+oM4xG/cskPby/byR/PFG5oTmcune92tc2EJGFucJsWv2h51LYNGXTapWPrixta6jWnTvuN4PncfxG53x9TWJT944PnfLG2qrW+DqGyAtzeGWMMQvRDo7K6d+9lfPPgU25/J4+hYvAUhXjZRixZxvGlqq9oS8P/ZL3fsXWetWqWwNWvGs7fy5FaE7NhZdn6PSETSlNkQwnJUPoMnZMvfbjg6uo+XjK9gW7oxW7ixNZ5+63SfZ7ZcnDYbLGdGw0Zd0B9c6Pe+7+GhxPv8uRQLaqqNRm1KFIvs6mMxVtfWwh7YtY994ad3sqnw5P5Tcx2bj7mhqzE7msuLZ7IFtlZwbb/bx+YgDlS6EFEIXD+CQRxaVlGU/Zb11V8PDq4fE+20QK1yqugHESAC1UiABG81nhWKiQicgQRkwwZZeqxj6qQVdT7P7U0+t+3VNNEYjagasrv3/+YRdh4um0vx2wcbw/dzRTb7zRehk5ifGajOMG3OZBTZLfKdQ4lEwOnf24ws7Ip166TYlcLQGhOGMouqQIAJWAIy0Gl3Jy2GDWXJSuXl67P5G5riqcunQfzCC4vFaQE2LRy0WkIBjtJdy5ChXdZSDH/p0fqaJ/+5NzW9Neir92mqGfB69TosUGttrGOof8tW/foB9q2nn2N5ZHVR9ICZELpl3Efs4n3nzOZ1oUApbdlP2Dv2qqvXrHlx9lae2YqQlXdex6h8RqpejMoccV9WaKj4cw8W7M8dLJSu3JDK3dA6mLh8RtDX2hENskX1MXZvsmCr6PQL23GlxJhLd7G6SW1sNxiv+uaP2YdrQuy6mZMZ2kqwFBasdacz1mOWo+2pq2U7S+aPikXjquc5bzgOsRvGAV3YUApuCBz+ccvAwNfHxC4tUAMXGkSACFQvARK81XtuKDIicCYR4OvWyVX7dQGvqn0t7NKECnUZ8nlUbyjIfnnvo2xhuchm1cXYCDKLP03l2PyLz0NJL84e33sQxlCkdU3Tzlu2vieZ/s63njhaWtU5BaV2K+LxxZyk8JJZUzmkPpTiV/phIX7Z3RnLlltob6F06ePp3A0tw8klk33uDixOQ1Y0xKZEgtacmgiPW+ISqEPmUrhAqzfdB09wLWrUltHa+Ko1X2PPHjiCYgVotoD3lMoG8+L+W+ZME1fOmGLCa+x6pmfo+z97dvtAxb4xKspx+BM2KnOEAq7MD3sFHVbG7d3wBt+9J2sF90DYPxVPvbO5b7jTDobDCxpqWBmZW/lBgVif697Hup9+lv3njDY2E3MvlU22I5FmGyCCn1B1LRuLmppL/8tjuw/cPrW+vqGg8GtyjE3qlZ9Hp2NUlbiXDQ4OyRhw7D86s7t27Vq1rq6OL1++HHtmKETH/+h9yR3QIAJEgAi8GgESvK9Ghp4nAkTgRBKQWUm7pp6/WddUlNZFLTHUOhDINj7y5LOso1RkC5rrWRyC64f9w2zWJUvZ5dOa2S93HWRDJcNBty8LZa9cR9O5w93Dmf9YBZsAqjv8IXEEPcZe4od9kfi9C+L5rn1WUd2XLy7dkMre2DAQv6bN7+1YGAuxbpuj1plPgUDkDnSlXDDn9/vZtv6RitiVtYAhwCt8Lp3exmY11DlRv1fZk8y6jqayz27oGfz02MI8qTFP1qjMT+5cCk/cjIv7LO7/OFE25TZ7sj+0qWTZwRxW4/VnCzyHhWnB4WF2Hjy80Zoo68lk2TNDSdbv94rAOXPM5LHhnx0YTnx5cHBwVyesC11DQ4PY3w+wvWSA5RtZoAZty8XKlStfcg7xHBwpoI3XXnIwekAEiAAReIMESPC+QYD0cSJABF4TASnIUIaAn4duXALlwkQOGccn0OhgnlFk58I/mi6U2MOwMSy5/kp23azJbAc8u1It7hxJK7qquhLF0r6+TOF9e48fT78Vnli89HpEkXzvK4lfKbg2wJqw4ZBZ/AxKip3/bDL7rlhtzV/OR8mxAmwKBdgV8ljA1RdPsnOnTmZ/de2b2X/e+wib21TP/vbixWxGNCi+snWvsvFQXxki/mdbBodvY5lMds2aigh9PTEilD96vHx+GgSpM3TO5IE96ETRky9hIRpn7nwGDS3cbFJrIwS7xTYOxNlmvFYzc5ronNrCjybSxSdr+M18525rbOGetE3Ic4eKcS8ZkuVLxOpLXv39D+T+KrU2/vmf//m6c889d+GkSZPgHbafO+ecc+6WYheDRO/vZ0ivEgEi8DoJkOB9ncDo7USACLwhAo4UkPC3skPH+9hF+At0fksDK8AW8DBKenXecA2b19bE9hw4BH+sKUaKZefZgfhGtO6942BPz3dw5DQ2ZQ3E3BuI4uXiUOmEqHt81BP7dMGyng4J/o5M2WhNl02RLJX5MLq6+eMJWR2Cffn6N7N3zJ7KWjw6s1Gaa9fAkOKUy3c+N5T5NMvG947FVRF1byDGN/JROT8TNg4mth9NtTXW71nqU5deF/U5aIih6tCvx1JZtmFghA16vGzJhYtZzKOL/lSW74mn913dO6SicxvDwj0pduWQ+xu/X3niDfyQapbNmjWr9sYbb/zJ0qVLr1yyZAkLoNObHIlE4oHt27f/Ce7GSfRWkNAPIkAEThABErwnCCTthggQgd9LQIom5uJ8Z9KweCKZ5lfVB9mS2hgzsQisazjJ9JntzKOrbPOefSxbKFqZYllDU7Xv9vf3f/RFe5aZ3Tcidl+0q8pdGZfdNfqsFKmaLDiL/mhb4yWjZaRQsv3wL7ghdNEdmBVQhzYw6GK1Lp31JdNsOJMTyXyJX+bTNv0KYvdjV13l/vr996NQw+vKPo8e/cT9rPgB5nYwvd1s+5v3zmib3hEOwBaisEKhLLYMxflT+HLhqqlhs/BlI43GE72Jsogji+1Kp9P3HzhQXvq+9ymsa4zKiYuLwYoiM7fszW9+808heq9oa2szIXYldzlENBq9csaMGT/Fe66Wwlg+J3/QIAJEgAi8UQLyHw8aRIAIEIGTTaDSgau3WHw44liDb6sPqQub6gTDoq/N8RTrC0dYDIvGtsHGcBDi93gyo+yOZ5xjiexdUDz8qvZ2NwKUwuhEit2Xz1mKK7l/2YBiXda0+NFs3unPFeBzzbHDsDTsHxxm+/qH2JbDx9mh4QQ7gJJkadQKtgolHZ9T3l4sysv8Ey3SZAk0wYej3zyQzP7jtkSmxut1i8OZvLrucB//fl+CDflDjKNT3L7BuNMN3htH0iyUL7AFzKlYF1ZjEidhyK56TigUWgw/9BWwf1iapkluaDgnNHh35X2rvr7+ioceemgxRK8jF7WdhDhol0SACJyFBCjDexaedJoyETjVBOAnlSLMXhbQzn37tIbgpVOaHQ4P786BEf5w3mT1/jDb1Tso3KidJcXaUNlUcsl0pnng2Cb5WBw4gBa9p0RISsGqXHzxwM+feKLx/cMl8y08mTUN29GlFUO24MUCOkeDzzRv22qiUFauh0zL6OxpfM4Zru9CmBM7wJpVqlQ4wisQzX2H+ni9x6Vu7Y9vfShZmKG5PQHfSJrvTmZZjcetlCyTzUUN4be11rFNOczkJA8I2VCpVBKZTEZJoV0xhC4adryQe5FdowVaNodkGKhycZKjod0TASJwthAgwXu2nGmaJxGYIAKrxjy3Cxtjc9A0bd29Rwd9M2oitocJ9VtHhlkCrtJJ6LDm1zXUilXYSNl02vDafLu8f35fNv0vAqWv0Nn2FIYvUEJNNDcr7zYt66f9RXFlxnRYRC7uQs1dnw51CCGZQee390UCrFFR2d68JW0MVTE6Rr8YcNTA+BzKSyxA9jT6P91H/nXbUOJrwWDw2ikKu3NROMrDKJfx8FDym04mc+Wl7W1TFWnkcJwDchLrurtPpnDnQ0NDvLe31zl06BCDjYG1tLQwZHsr9Zb7env5gQMHKsdH3eaqYEpBEAEicPoTIMF7+p9DmgEROB0ICFRnQBkrXtOHSgz/9dw+tcHvff5ooTxPdzkcTSEganneLYTnooCHfaguxHYlixsuxSXux5Z34u9U14laNPVaWElxzfv6+uK4vaqpqekWq1T8qsOF51jJ4Ch2cKA/X9xyuZu/q702JPLI/OYMo2r+lq4Zs33sGEofQvzndNTVeXcMDedkDjWbzWamBN3ae5uiVh18yIFC5q67U841eA0VcG0WFkpF8NZBkOLtJ3pIuwhH/eJdKHmWfv7550MQ4BaqM2jw8rJgICBK5TLfvmNH/3Nbtsg2zxwZXvkZGkSACBCBN0ygav5Iv+GZ0A6IABGoSgJjAkx5fjC5c35d5ANhl/7PQ4Xy3U/3Dh/5y5mTvnhJa4M4ks3H//q53Z+4NhT43srmiOKgisNQvrRVTmj9xMxKit5KO2AsmvvuJfXRf/no5DZvK1oT41L8l9/6xLGB2OSmd0EVKpZpZdOl4j4Z5i5khicm3Fc8qhStdvfwcA5ZXxcUpMxCIz+NZhngm4SV4GjRCrrhIpHd06TFpGRZ8u7JGhWmIyMj/Wg28Ynu7u4fSLE7PDwskOEVqL4rTbyaK5e75xvf+Eaf9O/iS9AfW/rsZM2B9ksEiMBpSoAE72l64ihsInCaEahk93YMp360rKbm15viieycoP+3NW6dxeDbNRS+TzVM1uZ2I++oWKliyUqVjO2VOXZ1TVSWb1y81ulcUfDHUoSRoj5ctryIq90Ln4UL9WKR4TUHSzwtY12N96yRd6pjyPil6GUQu+PC0fLgmcoKwFGBq7gZF7p8JzK8sHCMz1l+7GQMeS4ViNwf1tbWFnfs2PHZ48ePzw+Fw3xaYyNb/ad/KqaFQm8N/cVf1CC7m6DSZCfjFNA+icDZSeCFlQJn5/Rp1kSACJxCAgJLkNSn4/EsjukPu9S5MbS5lZnF4WK528XV85sgdz3I8qEGbvLh4sghGdvqU7NYTR7q5WP87+MsF2dhH1esEkp31WvOFtRLK3mR/9Xgey1byFJaVrUmD6SAfUH4osLaIi980hC4im1ahWPF8nFkeHEWMHAeypZTEcgvB3GCH1dELzK9/4vs+SL4eZc8u2nTalc+r/iKRVPzeJo6Z8/+a2R3BVu3bvwcnOAQaHdEgAicbQToj8nZdsZpvkRgAglgCZKA2pGiqrVW0xpiqmo7aNE7WChubXKpHXUoU+bGOwqGuXdLkmWh1ColtiYwZHloDSKxctlfwAogbDUd5HymDxleHbViLctOb8lkTMR6KsTiH4Wic+xTOhc1nkrcTBG2k322kO1HxtcthbucW8Ec65fMuv6o47yOD0nRK0uOWfF4/FlUZljzyNat++7r6tLN3l6rzuu97bcf+tAslGlwBNpIv4790luJABEgAq9IgP6QvCIWepIIEIGTRAByiwmU9prf5HK5Qpri5EqG83Q61d+k63OiaPCAThTosmZtxfGd9Z1sQv9GydSojMMD2S2zoI5tsf3prM+jKO0e6FsXhKKbOfuwui3LVmHh18RloyuB/qEfLiR2peBF0Vu5lV2Weg5aakRRMcOCcGcRzjbKfQx3nRIv8rjNAv1IOOuJx//pvk2b+N7ubtuH3435DQ3/r5LlxZeesS8TVfuF4g9xp9eJABGYeAIT+o/JxE+fIiACROBUEkCmsSJaalzaosaKfYHrhmH2Px3PsUZNa4oo3CwhyWgYxnoZ1ykSXr8XAdoJ66NZUcFM0xKDpmH4FcXwYiYoGstKpj1qZ1jze3dTFS+6VQ7Bi9QqMtNl01KgOHU3Z4r08Ap4eLljmxMQqPHOd75TZnt/+vTBg48/vHWrq9TXZyLL+65Hb7hhMV+zxkLI8ruH7L2mYnv5v1t4mQYRIAJE4PcTePkfjt//bnqVCBABIvAGCCwfK5kVVZWF9S6NeSrCy9yTN+ymSagf4GdcSZfKxjHDxDorNB6YwIxp55g4d3O+VNoX0A5Mt0w7+/zB/j0hhQfhNYb/2GZF05bVD6p61I9xVARfKMU7FttJ68g+VEkY9OI7iOpgeRjmgjrJEyIeZb1dHNgeymZX3bNzJ39+716BGN3zamu//L729tCPZs+uAWCUtat0wXPW4r6s4jC+qA23lcdVfRIoOCJABCaUQLUutJhQKHRwIkAETgoB/oVRwRuAd3cuNgYHL8uWzJ1hVV3UjMYDEMBq1jB718azx2UEEDgyszehw82Fq5LhhX0BC+zK9zBWex5jc1HdoBK/1xGbZIAV+0XXSW19/Edz6BjjqDN4eLEX6dlVbZHC3ZKs2CAFsINzUTz5VRpwtFcc9ueRucXvx/rN/f0P3bdr1+Wz3W47199/ycpcbk+borj3NzcP4JdhQ6+ifO3Snp6dbOXK8R3B+TBavgzCF0XWuPQH0yACRIAIvIQACd6X4KAHRIAInEQCMnsoBezMGk1tgeC1LMPUDpZKB1s09d21KqoHwCKQN6w9/YwV8MZKHdyTGM9r2rWbq4asyCBtAKZtyUvvEMBMk8JRVjawnerP8I5P1K0othTqCrK5JbRGxvOOFLzIXlcsDRnTqhRsGH//qbxdg+83OOecl0p/tuX48Z2/Tafddek0X6LrTXX4MoQufDH0Ie6oN4z3P93Sctv9n/+87+IpU97WMmlSIBwOH0R95G9A7D5CovdUnjU6FhE4fQiQ4D19zhVFSgROawKdELBdEFg+XVnQrOsKbAEOyo+Jp1K5RKOqzIooihCmxVEp4Gk50fV4P24mLFs3bgPAiqpaedlfinEbFRkQUzYAW4AUvNIGkDbNUZGIyVX7gO/YGbViyHJqlmS7sCJ4HaEUbKs05Jg9cg4TYSVZBZwc25uamm49F70w3Om00+p2qzGfT6h+P1S5JmDudsL5vL4omfyayOXY3MsuYyGIYYwF9fX1N6Km7y0Qvd8h0SuR0CACRODFBMjD+2IadJ8IEIGTTqCeq4ubdLVSbxfrv45uQ6vhVk2tiXJu5dF8omiaG2QQw6PZ4JMez6sdYNwGwJlYKEWivOwfYex53I2iTFnUh4pkJiob6MLZKvcx0fHKGF5l8NWjLIMK4zNdeJOs0hAWQsZd8SKrwlGRqi64ckZF8OJ5JFtP6VDW4MtNoLFxTpCxT8WQ6fXjS1HQ7WYuNKVgDQ2ctbYqrLlZY7GYQD1he8ovfmGjLJys9OA4NlLvquqgmcW3N23aNE/aGqToPaUzoIMRASJQ1QToD0JVnx4KjgicOQSWj2ZrlaimnFsP/y68r8iQOltzjhObjMc+LERKlc3SHsPZI2e969SLrleE7VVYpZSXYgtkeBE0YzEfHA5uGb5lYcGXPfKKH6yiJ/GHXsata2gWJzO6stmHR/A4njflYyxaYw6aTlh4z0SE3TmazWchxm71KeiMAbOIQO3jSksMLxrbBSGDo9HKJgIBVBOuUbV9+9Shn/xE2jKk9UWmeR2Px8P9fv87x+ZA/76NgaAbIkAExv7IEAgiQASIwEkmML5gLRjjyqwaWW8X2dEBwzjYrKrn10DjuHFZHfaA/Y+WSoNQZ5WM30mO6TXt3o1WxzLDK32vOcuSCdIiLrCPLfSyWdIUUjNKgc4fg/BaC/0oH1fhcGBpsCrBYi7IpMvVXaYLsOUiNgvPySHZn+rYu0a/DIEgP9dACEXcKyCjnpfeYmSjZdtjdPiQudxKaBC4jMPKcOg3v2EGythx+fs0OgTaFk8evz/+JN0SASJABE75HzZCTgSIwFlJQGoUOaY1qGo0Br9uomwwj8n+7Nu1te89T3M5fghglP8qHGCsjDfLVOr4Z0Y/eWp/8tWjWVGvjovpFcELseVynP3IN04LIDS3ENzAKrZhw8nK0Lrx/kshy1bK7OTExv5qpLDQjiuVigxgfcywg2263iznhkw7Q/KafwVVGyR7PDUh/zbYqNKRAjtsLAHegxCzqXyesWSSsREk0hMJxrNZZuP5NMTw0XKZDaZSrFgoMCxaY8eOHeOHDx8+io+z9evXT+TvjwyBBhEgAlVEoOL2r6J4KBQiQATOQAKdEFBdEFIQugtQkUEL4qL6iGGo7aoWWOSGDOOKNGKKJtW1tN/tu/8JRXkv6+lJSuEI1YKbUz+g+ORxveh00CrTurKE2hSVb0b2sVa2FEZlA61s2fF/KBa34I0K/pjav9GVmz1MzXPT/Jl8Dh+rLMQ69dG/5IgyDpm+7cACvLAX7gV0slPyDl/1X7W1rmkg7IbSnK17Qn0tLXf3WdZn+eDgRpmpluL9JXs6eQ+UFTjXjyhKN2TtmwbALQLGAcNgOgTuJGR3Q7oO64XDDIjbZKnEelCKLDdtGtu8eTOLuVzCxi/Rvn37xPbt238pw4TgHU0Hn7yYac9EgAicRgRI8J5GJ4tCJQKnK4H6MdEaVZRLYWGQdgDZy5a1+MKCh0JceD2wDDicZbN2Y6Fw5dJS6VsQuishGuW16lMlul6CV6YHcXwBv64pqzSgqgF7pmDe9HexmP9cuEbh6RVYuOb/TU3NNTwev/tdqvq1Eab81VRNYXfpnkt4ofRRKRqxD0xsYkT7SyYE5wKy0xyWBiHtC5MVxdsBEdnscsHcq2FtHrS5ab45yHnn5sbGj5w3MPAjKdoR+8kWjnwVAl2D87yoUOjt8fnYfsQibSMoXMesYpGlIXwj8vcGgleehwIqTGxFBYeuw4fZ5AceYAFUccihasPQ0NAtv/jFL3auWrVKWbNmzcmO+2V46SERIALVTIAEbzWfHYqNCJwBBCBmlLdhG5o8WS8ODlzYgHq7jmXD1IDVXxC7rLaW8XAYMgx+zXhc5QMDZqtlrdjc0nIp7+19DG+rdNg61ShwXDkciESGrCgroK2wx2FXfATiql56jpEZ1dF1eKbj3LVD0759/vDwDUVoyUuRQT1HV2+9IxJUbkxlb6kG0Ts+F2lf0GEWQXk15uGKQMkvrssFYYFAJZfuZLNWIJnU2ovF/17f0rKB9fYewGdPmuhFVldFOtaG2LXe53WvvlJYn3sE/oo7kRxHExLLxLEhbpU4BDqiFMi0CwtfllKcqw/a9qMbt251GvbtU1SXK53NZr8BS8MjJHblry0NIkAEXk6ABO/LidBjIkAEThgBKWhk5g6bAw/mnOWctWPRmoOFRooGE4AOOwOTNVZDIehdSCt4M2HG5EqxKDTb/gjE4uMIhp+wgF77juQxpU6cCRtASGakhS14BLFHkRENeVG7ARUB5CIqDy6xz1OUWz8fChX/PpPRGMSw6Qhzoa7efEdNhN0YT02o6F0M0bgZc8GEVB8mJD28BQhIVde5F9lUlPka3RC34vVqWBxmhS3LFzOMv8NnPiIF72vH9trfKX831o1m74OfDHi/uszj+dB0j8eZaVkcGV6lW9NZFrEO43enDrcQvIqKemN5xNlrWd+5f2DgFnm05IsOSWL3RTDoLhEgAi8hQIL3JTjoAREgAieKAISS9N/aH2pomHeeorwfnszLy6bB6+DXteTKeo+PlSFkPLgvcMmaqy/8OeL4JDeFaJQeUrkfbOpqiLY1J//y+vj0ZQUDjogiyPBqEIkOBg9qmhKS2VApEiHS5XCSSaEMDTmfcRzlGcv62QOFwrscRVFhGzAXetw3/6q5hr2jLz4RopdDVCoQlSbOA/Or6nKUJZM+WEegIoYs+eXIpg24hXjHLPGqPA8eD5SvIuBT7lw1ebKHHT1arkz0BP6QWW95bqe5XDNWel0/ucTnXdLocpUzqIKxrVjk9WXjv3fW1B3vFeIvhzlv8ePYOAcWF2J3jvOvHhgY+P4qzA2/DxwiV3R3d8spMtgY7BMYJu2KCBCBM4jAC//CnEFzoqkQASIwwQSkGIECcVY2N6+az9hnFyqKPgWZuToprDBMVWdJLALrSSXZZFSHlYvARrIZVsKK+1Z4eZO2LXJMTLk/FrqcJzKP4SOW/ByEr7JuBeMr1p1cX2wHBBnSmlJs52TEqLlbsQEIuXBKzkHaAKTohaRHtVjO8nm0Liu710Wjfe5c7j2Pm9ZaU+HC0VTjXJ/75l9ObmTvPDpwskWv9MJK4aesBnvJX2ZQ39XcPBNZ3n9pcpwr6qF1g4JrQUVlGTT5GMLir1aIXIYstYDg5fK+bUNXCo7FefVl0/RjPyVMX355wc0bHhxiV5Fi93KP56KL/O67L/b7YmFdL4+Uyu7Hszmrq1z+yyfK1rdhp2BNNTXfiLtc56YVRSs6zoAxMLAXEVQqeMj5yWggct9wULQDIkAEznwCJHjP/HNMMyQCp5TAuI3hgqamT8xUlNVTIbLQNtiKwgag+nyKzCSqyOo2YfV9PJ1mA1iMtGdwgD18/BjKCGhsdiCoOOgjkGFi2gy368FNrfXdWc7XDVjWL3h/fA+DipNDrFiBS+Lr2Mp1J2RR2wtiEbt2IKGg/JiyzBeYF0KuV4UgL2KxVBrJUBPSDyXJKvVh0RwByVBIYxX6GM9vt6x3aJz/Tdlx/uzJsvlNVA5gtgrRG/DcvLa9la080PN/RK8Uk3I++PHHCMqKgNyFj38BXwoQt9yHjJ99Oshqdvka3juDsX84V1HC08C5FYLdjVv5LmGU2XA8zkbwuBYWB6nbs8kEcxWKrACrBjKrRlYGduKGArHLpdh9Z8h764Ve33+dF/BrXlUp9+SL7ocy2cSTpnHd9rL1FN4nM8B8bzwuQ3hiPAQJ6p14DbeUyR2HQrdEgAi8JgIkeF8TJnoTESACr5GAvIRuh+rq2nHB/8sxx7E8UIQhj0dzy4VpNTWjl89RP1WgrmoNhJWJ+qn/jMzuEDK+XlncyyrLlmDsWLmMDCmK4OpaR7NLX+X3uD79+OzAAwmb/2pvInsHX7euoseg3fj6zk51eVeXDSEkBd9rHa8oFuWHb2pqunIh5/82k7H5zUh31giu1GtuNJmwWD/ibUOWF37SSsMDCzVitVKJFx00pkDJssV1dQ0bh4a+BdFbeqbMfoBcsWqpirkw5Lv5F3OmspW7D1dE7y2LFyvf3rzFwm4qMcus6vrly9VLMQ+E8KrzWAUhzjqZMreeCSn2pYAcm7D6qVhgVr2iXT5FVS5FxYg39atqzWRIaVgC7IjHo7plZhq+XSFjR33besQ+ODTEDiKzuzUxwlg+i7l5Hd3lVQ8KseMb/f1xBIKzMJphHzvO676RX4J+iRgQK/uTWPBby/3+W2f5fUJlorwvm3c/nMlv+lm28Cf4lrFfvvdFc5IaF98oKgNfM5iQv19jj+mGCBABIvCaCZDgfc2o6I1EgAj8IQKdECddyDD6FeWdAfhYcW3c4hC8ivSKygVS0vcK0SVwKb3iF0WWN4jtMrSP/RIWqt1bLNmXq151lktw1a3xn2QLbBqE8GwU6W1yu/R6t3ZdwKNeFw3U/tM9k2p/mxDsh3zb3g2sq6tieViLrO+KdetetQyYFItzIZAZVNXLxeLHA4GZ5+j82i2a+6ow52+eD1EIG4aY5HJxl/S5YtRBqGcgzvsh1Gtw/+nBfhZHVrQd1Q40RUeHMOH7FrMeTjTU/PQyt//b/NixFQ9nCz+CqvWipIN5Tihw808XtLOVzx+4hW3ZYt8OOffeN82PhllYZliTOMS4dQNHf0H0Vry4HZ2dfPWoqHdY1+jlfHx9CN5SX7OkQVevCin8WnSwm9mo6yrKv7F/LpjoEWfbHoUrHjznl/xle15YMSqdyZBdh52XNUD03gXR+1XYS/7U54VX1kYz5bITEmLubxuiK/lgci2EprIawa0ZsxHg7mseYK7hc3Je0T+vj/y8Mxi8stXrsUxUitieybufyOW//7/Z/CfwelaK3ZcJWhyaBC4Y0CACROANEiDB+wYB0seJABH4HYGucUGkKG82Idhki9gipFtBvqWSn0PNLogsaQOQoktaAqQR8/+3dyZwdlR1vj/nVNXdl769b+nsC0kgCUvAICSiAiM8RsQEnuhHZSTijIroe+PMPCXpGf3MzJP5qDOOqLiMwswIjQIODpsszSKQhZ0khM7a+3b77rdubef9TnW3Ig8xSNJ9If+S28u9VbfO+VZMfvdfv/P7N0JQvqeujj+ISCxRrsDjq8nVsTA/KaDL6wfH+1rKlTmro2G2ELaC9lDAbY6E2psioSsTQlx5z9mrfp2R/LahivvTTV1dfb8dzeRPEFx+RXTrBvhaOyeForJFKLF4RV1ybZsGsaiJC+ZALI5zoeVQJz4Jow1zLmthwwhMi3TlSi4UWQLdvgxkvnZBqF8/MsTeB7HuVEwWMDyedj2tTfCV8wPG32/n9pcOLu644aqh9O335oofzHieXoK1Y10qfuW3Vy9hz8SSzy9Jxi6vCxtL4sEAu3DBRbtNx7v1kv+481s4k71l/Xp9RWMjqrhdri8Cu7uV4GQfa03NWWgE3lOjGxemNLG2QTfaa3UNJhAchA8HZdvxHgHDZxzJ34n4LqU0sYhucnEaxuqnYigBj+vAId6RfcyW4Vr0oRnIneh+x/WwWIoCcJ3gjULqN9/e3vQO3jd8Dd6GTVkNjrjCOi12T01G1pwbi/7gHcn4mrpgsJKFX3dbNsceyZe2/qpQ6lQlXPw5eLXYVaekjQgQASJwVAiQ4D0qGOlNiAAReCUB9GQIZiF28WATEFajEKr1WBiVgMAS+Nks5CG2sujwBRsA7uXv91y+S4i9lXD45m253CegiFs8TbinJKP8/y6e0//5l/s+9fToxDkthv6hk5KxpsWW4zUWy7I5HBQLUol17aHgunjEu/a2Pzmjy5b6X71Y1iZaCwW+eScsA0pLoSLa2c3YxbXh9jXh2Np6oV2Y0Ph7G3S9vU55ipWD1vXkd8qWE4W1ADYMLYBmDGElEFVVFFYMrqrUsDN4GHMYc4jBErBTcqlZtufomljCPR7Bzf+flyrefKQ6LAgHo0XOPveP7fXuxabF7syVxDP5khgSmlfX0njlyvokW5CMsfpoiKXQeCOi6+tQyF73yJWXXDKYyV+0qeve9BRT7VMdDae3hcLvhbA9Fw0YVjcFAxG18Ex9iCg7rhyuWN6AaWr7IHR7S6b9EhrBFaMJoeIVinioeK8Cxu0nYqAyrcQuU4kMEt/xWi0E7/raWu1+VKvDZgVzFWxxKIhX0VyD88/duajthL1W5eObDo8NHqHoVckaqkrtvCeV2HRWMvqd05OJVMzQzeFiOfRoOjvxcKF41fai6VePsR9VcqcuNn0jAkTg2BAgwXtsuNK7EoHjlQC0C9QL54fGp1rEpiCu4hBZOiqJSyGunjt8mD010MsWBYNscTjKerH/u5GTdRF33bmpcM/ltv1nDxdK12el12FpwltXGz/jn5fP+8xzdfU/sj33fS3RcFN9OMjtiiXG01kseEu7LeGg7KhNxOfE41dkHPvmzru771Xj+KR6NNcunxMJntuAWK4o185pDhrxJASeklhoDewOW7aWxvieNy2+RwT0s1ANtXCcg+9YcDZZGcVYpXqonGCV0oDXVqFK+uFUit+azWg2RKKpCXliJMQbElHxk/GsmIcUBFWVnhsKstXxiDi9NsEOQyg/g7I2BLKbDAVkMhrmzbUp0ZBMIB1MSKdk4vngmY7r/ff/OXHRP6EH3QV1hrE2qesnNKIKjExghtdY2bLdHviGIW7Ffgj/3nKFW7YzXgO+qwR3TzECNf+C+vkYNOcEnhvDeSGIWRg2Bgh6P4LMT2YollgeiRi9sJ48axg3l0Phu35dLF4f0LWEputyYSQoxlz4sBk/b1E0tvOmJZGNm/Yefuz1rCNbQAcL6HxbyWWNtZ87rSb29RMTcRYQwjqUL4QeHM+8fH/evPxA2dy+GXTxBwZQaSMCRIAIHFsCJHiPLV96dyJwPBLgSDX4T5QnP3wADQSSIKBDcHloEZuDYPzL/j62SNkDDBfCrcTS3OABx2IrQ4ETClz78b+21g50V2znxmyBP5ItiLQQXm1z03nzoqHzYAFgDbEIwnGDLIzb+BVULV8cHNVu2v4CeySdkyubUjmma7kvzG8/LRUQl9Ub+llxTTu5LRjQQqiIuhCLJm75H4ZYHoRIfblkavvNysiIZb9UgS60agIXoCJqFFCdzEEYZvH+MYhhieo0x7Eci9bwJv41TaMy/VIgmE7E4o/vyGTOQLJtfUVw7xSI2/MbU+x7Q+MjD5cqYfweX4Z5z4+EWBsE8PlBA13byloRtohoeytrqK9jkcYGZJ8FuJHPCz406LXEo6ef0Vx3SxyZxZDcaswyXbHkcMlkBwslfqhY1hDjVcKYjDaEip1t6LGVwWgUFfBgDGZgFKu9mz3JD6Jy3YrjYxhzAIJXWQc4vLrCLLOUb2tgbADXpowPHLebxezJgcCDPJc7/YGi+VxAEwbXNbkgGtIHbdtF6+cWfLB46Kdrln4aNovvKgjKdjHlK4aKZsoarUHsqji3+J+1NX7jjNr4FfNjUQdBHWx3Jhf41VjmsZ+NZy/Crun1+GPxPRK7ChttRIAIzAABErwzAJlOQQSOIwK4R64KvPwh3PZ/do+un6RuoKPBhECmKysiySBhBNh/l0peEGLK1TWtKSDZi6jx7ciX5Jpo2G0OBVrPiEfY6ak4212x+PMCQV+hgKcqotEwKqg1SdFUk2RBRDnYENGNWGiFNFyz8+5HnDZDl6vi0X+vbUotaAwYiBfgzMIt/4LlyN5imfchcuvlYtk7WK70jlvO3oTnZU5By7ePBI05cU1f91nO9WFMYBxidxQisQ8iMQQbRh1+T+Pn/rER1oyzaY4jlS/5f9qlbWeGxedPTnuB7QXzF0jjXYBEBu+0RIz/05IO63sD6f91y+h4R0NG+xAE79LlsYi2EHObAxtDrZZnOirFRiTKpPIJK/uE+j2bFYFAwNORYtA7mpaDpYo4kC+KvkLJylt2IC4Znwex/65I2FgYMAQsGUlDMXa9QBaC/AAqy6NmRWjoVHYAnewOoJStQYEa4N8D8b5tYpy9H/Fvc1AxtowwG4SxAZVzvcbQ/2KH4J/aP6/1b981MPLFX2SK16IYndDVvpGQthtWjXS5oq1Ixb9z69rlp5bTQ1/4SHd3rlNd8C1bxNaHHhKdWDzYEYm0vK8xecfZ9TWnNUciVrliac+mM9ojE7nv3jWe/TQ+6jgfnPTr+gv0cDhtRIAIEIFjToAE7zFHTCcgAscHAWgqaBnG/pGxpBwevsPjfJ4TDMrnsPDrWbTjPQAx14FK39kQd7A7CCQysAqqt6ehe9lKVG5vnsjz+9I5fV0szBZAuM2NhGQ7bvkvheDKV0oih4NizY2oiNazWBMqovC9GvDR6oNDbGVzXfAra1cEgo4rDCFqVEV0DBaFMZzjECqi+4rl3OGSicKuc3CeLornGHp0eSK6oEET76zRdPTT5Vjw5bJmVG0PognvHMwkit8RN8YExO5BnOff+g6zeXC1zkd1ORkI80FI+3MF/5Mybvfvntdy+91l63t/PTZxORaurcxBgJ5Zm2j/1NzGf/j8snn/49sVJ1sv+HVt0ZCoDwV5Ef7ZMng0ApoJERpHtdmPCsO8Oaq6RiAgBjCmm14+zDiYNaAiflowYEMsGx1GgMchqvG0pvy7fVhoNoqKtfq+D48+x3224NgHFhtGoc8IXPY0MjLgqYbbWMq0ZfFtqFSfrOm8DaI5A3vzf1ueCJeLzAka1jwkYcCzu/XeuS0935rIT/ysUKqpYPHeudj3hGRMPJMrysfTOXdlU90n8m3zTrr+xNVfKhbtPbyzUzlTvHOaateenIjefkZDbUtNMGClS2bgidE02zGRv+ohCF71ZwRo/ei64+P/FTRLIkAEqoUACd5quRI0DiLwFiewFVW7TlTvvqLrl8NLun4Mt7bPhM90s2OzF3D7/FF4YJ9AdfdJiFwjkXhuzDT7HyiZF7jo2XCarnmXNqXEz9K5O34wnJ5YGA6+f46h1yyEDQBWBjYHVdEUqsMCVUojioqoypNVLX5RxTTyBeVN5cjd4oNDo6huWmx/rsBx6394tGxjPZwztNzQrE2RYMvimsRJKV10JGFvgI5lRfhXB2xbjlRsPgLRGYAQTKPiulei+Iz5OLjdr2wAuyF8f4Dq8KU4dxBC1ENS1hi8yXfbtncCFs41cf6BcxORD6yMh7M350oc1VDNxFKx+bWi1Q6xnee0NbBFNegEAUtDHJVndGdjGYjxCVSPe1GhXgA+IVR4y+g29/KBXnbw0GFW7B9m78MxiyHsGwxDBjmPwXLsZi1b9qOCO4wFagchcA/A2jDouM+gMcd9h6R790cddvASXb8raLvLT89ktO/H4/JJRKtloHlbcQ0ctHB+Eu8xPxRgLQGNDRVLX3q8WD73OddbvxYqepWUFkT9ok/ClvHeugS7u1DmO3MlHKexVDLGD+lBbaw25c5Nxtam0BjEq/Xy911x8X/d+OhTO+YHtM5TGmrjsI9UBguIHBsaSz82kfvQi7nSPbeAJ8SuugOgHrQRASJABGaUAAneGcVNJyMCb18CK/xlYPiC0AUT04RQ1J5HFRNBvA58pPyjji3+FF3Lntd0G6FkP35XQ8ON7QMDf3Vntvg5VBH5Oojej7TWv/svFzR//LSdPf8Qq098t5CMnZ2H6O3TBT8MgRhF9dVCdTQMUSsNpAlAAKtFXKoqWoQo/GFPX1+2Yu2JMJZdEwpoH6iLd7QZxplYpNYYxlhsCLoiPLw9qIgi2YAdgmWip2K7g667rd9x7x1h7KF4JPqVQ0Kc6XqeDduCyOMcSHNgKyGEb4OX1+IRvhapDHPCAfag64iHckV2ur84zXBaQ8HkZ1vqUMGGZxUi+gW4mWshaGvDIWXHEI01NawuHsVKLc6a8F4D6Ql2uH+AHR4YYi4W9eVGx5mF9Io4Xl9pGF6ksV4iZoxhrNo45j9QhucY3/tN+/CQbT/Q6zoPZ23vfuA+PP0nayNi1hCptjqNDxxLcCzyda1vxxN7bwsGi5jfUisS0Z4r5BOnYDHaHNg5vtZU03BKJrdhe8X6xEu2ffUOy165Fhm5K9E0ZH4krF3RVMtRNWbbbI+NBaNsDeY3LxbWmuNRrzYSkSGNI1hZfujja0/8UBoiHZG/zoGJfPCR4bHtdw9NXD5mWa9uJjE9VPpOBIgAEZgxAiR4Zww1nYgIvL0JbILAwgz5l1337j/XtPf8SsrPwqS5YaXnJlZgzT4Wqsm5krlrBAuEdO3rA0J2PrWo/aefGhzveTBbXJbGQquyxmO6jHb95TtPObiyIdnaEovyBLymOmqDOVQ0yzYqmqjiLkKFNgbBKFE9HodQ3I/b5mPjE/IdsXD/8rpke0PAOAvJBiH1FxySGFgGldh9SDQYQhX3QBme2Io1PuR6jw659u0DtrcNu+2avjpJxv7Uct1/P6hp52XwZCMe9XishFgdgBC9C75hR+Pe6RjU+toE/8FoJv+18UxhbTTSsjRsw44RkrBjeCfGo9oqCPV0AU2SExFei+ObsDgtiEVq6L7GAiqiDSI8faiX9eza54tcJDIwPZmQSFzwlLgdzuQZvMcMY3dRhX5qyHK6D1Uqt2E4z+ABPT7pI0HJVLXshZWXuSXX7Ulz3p/jvO05hDqggYbxN6Xi3I9Z5o/eZwTPR/MPns9k7nzGttctgKht1o3P3LJ03ncue+ng93Meu+lJp/IhCN9rVlTslafYLlvuht250bB2flRnE8ivMNGSuCbeyFob6kQdWmZoyEx2ikVvJea0C8Xbh5/bq+9IZ2+5o2/0KoxnYj3q8F2o/Kux0kYEiAARmC0CJHhnizydlwi8PQkopwBiAlxVdVSP9v1SXPaYZ1/c4brrVqD5wjLssQCr85Elm4CHdvO35zaz5y2H/bpUFs+aDks0xLBgLTEvhfivOiQyNCHSKgohKF2HjaD6eXh4mOVR5QzATzsyOMwyw6hbFkrww0pxYkv9GRxVYFRx5QBSDcYR19UPgboPC78GLWtXr+1sG7Dd25FJ+yjGNq4uAW6zq3vsfANkaDcK1AMDA+r585ubmz8yLsTFSGNI4C9KDhtCLpFI/Ly/UHgnUgw2u1zYZ2DR3afbGsS2ovntvzs4MLjdND7fVjSXL4+GtMXwIs9DNbcxZLCAaj8MS0AghmIohK/qnYzeDiwIOwbS1FgDhGNRCdtCye0vlLQ+pEf0l8p9QxV7G2LT7hyt2E9gTLvx8Dc15g9ivJg5x5hVBJiyCVh4cIjLng7HOW1UiC5YLc58L/Il0ppIrNC1q+/nztJz3EMf4C67+jnT3nEi8oxbY5FAyHX/Bm/wMblxo611df0wY3n//phT2rzbrHxpZclsPDlpyRNqYrwd/uUUos0iwYWsBtnEGrzUEnYJHdFsCQQUz21p0qyDg/fe8VzPpar1MrzDAuMjsetfNfpCBIjAbBIgwTub9OncROBtSADiS6p4quX4/rcMDbw87zooyOvGIXjhE/1Yi2X/yULLaD8Bt9sXO0EP1gBvIW73r6yJiTR4jHqWdBCblWiq4y1YoNbUgKYPsAUwiNwgqr0eH2NPvvACk7ASoH0x/K06C8QjDpwJAtmuYgTRXSqRAQvV3GHbfnLYcu49aFr34a2fxMNVyDFGdi3sCg/hAUE2LRinhZl6mQ0NDd2Ib+rxOxtevHHMdviDhdKVrgZTha5FT6qJdj7yjhMvO+vx51eMOd75e0xrc1O+/P65E3mGLFt+2qK5rE3ZMTDnIB5qgRqDYFcWCwEhvBtjvm/PAWnqhjbuuBNDZfMLo9n8z3DinDq5GhAEKd8wJcrxocKDsPXngpdeucktmFMnQ/iC521AlfprP/fk51TiRE4T1iojcP4DTt2zuUXR9yf29d72QiV8MRYIOmjgcel/njT/q7yray/eTI2ugmi5f0lb9i0PS/bXeyr2Z5fkih6afoiT2lDdhReb40OIn0mMDyYSlXcRDHIDczlr8bzma5xgOHHeeZXOzk4MmzYiQASIwOwTIME7+9eARkAE3nYEXiHGlFZTqViuYOzXWHP16z6Xxfts7+IdFfuStlLlvCXhQHBpxGLz4YNtxQK1ZuhPD7f7gx1tLJpMosdtPSwBMfQoLvsdzpJoltCBWDJTCUjP9RDXxYeKZf1gvsSQxDAOv+ujg5Z9O6qi23Du31gV1EBUVRRjUz0nPIhCJcZeS5DhZX9Tfz+qn/3fIeJV9VR1EFOJB5vTEK7dWYheVHrRoEKzDf2WH65dftUV23Z9t7G1tYAuZhePeI5nBIL8VHRrG83lWWw8zZohJwNInrCQw5vGIw+rhhsOy/1a0NGNwH0yrn9pdGj0aRVVfAmSGEZwzm6ME7+qsU6Lcn+Ar/Vlal4C40R2rnsNqtTP32PLb4wzHs8IUVmlBxdvt60df16X+Nmt+aK5NB7W22LhoCaCX/7o+lWfOTVVr33m9vszDhbu4f2H57a2/lR47mdHYB/ZZtpsSW0tK2DuebRZrkWF2s8mhmCvwLOcR6U+V7HmwxISgdgtYwy+tH+tcdJzRIAIEIGZJECCdyZp07mIwPFHQIlF1d5XbZr6gp/z0G4/KTreT/Y6zrK9FevSJ4rlS+fkikuQyiCW4tY5uqGxlfDeYjnapGRSbXSRLiBVG2BUESuo9N7TOygH4cftr9hsuGje2F82by277mM4xfjUeZRC5BtwZDdEKgby+6qiavfX2n5HXCqhPLUp0asp0TuBBW1ojnGlK7ilMoUXJmPf+emZJ5Uve2zP4/PnNLK59XV8KTqsIT+X9WKBmsT/8vAeRxDzZSGOLAv/cT+EMKrU8l9XdGi7hga/dPVzPU9vXM4CXbuY/YoPDtPnPtLvSiDjP0QRe94Py/D8Pibl9wFmTRpjXa2J8EWpxId1XXd3WY4W0ZBxXBP/8LtikQtq0VL4ro++f8hy5a179o9e98Xnn2cLUjV8cVO9XIK5MFyDwQl8IEGFFwv7WASRc06lzEYnsiyNJAukSLgTQfO1Pkgc6dhpPyJABIjAUSdAgveoI6U3JAJE4PcQ8G/Bq6ofXhdbIEJhedgjPa9zrK75W2jPu2/U9pJPD6blmW2NfAHsCqPw7IYRyxVAPBjiGVgZVcUsKqJSaN72isN7y+7TJSP4v8cL4/erc6o3vvb3WxXULkdjwxQYzAiTojcDMfvoRP5KDx2DXV24C2viP/7JOSd+4spnD28LasW1QXRfwxyFOgAN5Ng4BC4W7TFYXlUGrzyE2K8GJEY0GZp4nhsNeG+xccVGt2tXlzrPm9nU8Q6+qL/nn/I87927KtZ3Mp7clEb192TBxbtr4towBOxETYJ3NNextkgolQJvdLGrxZivbYoGNl5Xa2z+yhO79w0EywsShuYFcJxq6IH3YznkE0dgKXFxfYpoJgJvNR8smPv6tdq8akaBfF4Svm/mCtKxRIAIHDUCJHiPGkp6IyJABI6QgBJibufkznA6MG394sXZvXv3HmhPxlZ3xMKyo66Gj8CiEBwdQ+yYy2qUhxeL1lRFdAy30ZFY4P3tqsV6eizdtemx5+7/5qJFwat7epS4ez2rwhEO74h2U3NQGhbZE3Jzznb6HxvPbnWwXA85t3JeMnrD5R316TsGsFgNeyh9aEP0FiDaE8jhRZc5HO25piu1QxD182FxYFqQNekBlbygRKJfDT+ikfzhnVSlWr3fBMZ6KVoJP35XvvR11U0uK4Rchk4T87EQMFwTY80dbbIeGcdGQJduueTE0e4Z3eC+9fWnXn5itGwtjOllBz5d1YsDiRkQ71gsqFo8Kz+ySpbIWo4B//S3uh/Z7nQ1Nh7NOfzhWdIeRIAIEIHXIUCC93Xg0EtEgAgccwJK3PHu7m63sbFxJ6wJq1MB3RtHYkMQUWTqlnkBi9VisDAIiCoTldAMLAFjEFprDY5ObX71Ujuprc1lPT1+BfmYj/i3J1CiVxl8EVLmdeYr1sATo5nvwd7gRiJBcXZTTV3FQyeI4ayGtAKWQ7OHBjPEleDVYNBFFrF2MFc0xfi4DC5oC+rYJ1MxF+MtH20YGVHF6qO5KTbKisEx1m/kOd/dnSv9eMR2G/MNXC5FBIU81MeSp6zm4fY2xgIG1wv5ACIrnJa4terv1q189NMP7NwXFGIhWrzZJmrZBcvmSVwXmH9VZd3GXI0DmfwT/7n92f/YsoWJTZ1dVN09mleQ3osIEIE3RYAE75vCRwcTASJwFAj4whHe3OvKjvOx/bkiF7ABeIi08jN04QtVt81Vmy4senOHypYwCkUeRhewku0twPnd0dmrJqqxu/hiQMPeoCq48BV/ry4SdEquJ85sSLz4q6GJGyZM66+xoKv1cL4sUfFlSXh2oeV/9eJ4+iuncvcWy/WaDbxRSIhFR4Hn73sLnzO+aAB5z9z2tuJgLsdvHRjzTktG+bsRjxZC9zqpHogfQw4bC2QzgmdzsiMReSfE7aVFy7npsOctyyKfNxVCagb81EEdrZCFMCbK1mP7csWL8P42nyzf++f7fYOh54kAESACM0lA3U6kjQgQASIwmwT8W/iDg4N7oGk/n0fmQU+mIHom8s5+NF44MJFjB9I51oPHnomCtrN3iIeRbqBqqwHGfIG4cfny2RZX9pfPXq9j/DcsSUY3I8dBD2iIRHC9nw4PD38rmbCWmfnsj84Ma/yLC1v5VU2Jwo59+y8wM5lHK1JGbCx+02CQSHBuH+MLgSEyN9nQsBofMDreu2yBbGtqEBOhMFu4YhmzYR/h+IAh0KCDYzEgXAq87Li8aLkLBzaf+6JtV9CEze0czhcO9GULnuE40imZ5WeG0p/s3vXShr6+vrQ6AbbZvh6To6CvRIAIEIEpAlThpT8KRIAIVAMBdctdg+j955aWllK24n3VdL3GASxQg3BEVVTIMHwAZcf78VBv36qzFrSuRiIZiwuhuhhXxdbZ3e17ZX+xr/+Gc+a1PAMVL655cMeTi+AvfumlnnxCFy83o0y9Dl7ZkZwn/3ltbYRtS3MTotJCIoUSvLDTzogtI4jzlJC/O2Ha2ruXzGXza9ExDYvXDqKL3RxwDYeCzDVNNpbJsQlcAzTysD955wux8fHxMcDemgjo1sYF7V+9cmk7y+QL5Q33b/8+qifetZhz56QHuSquCQ2CCBABIjBNgATvNAn6TgSIwGwTUGJPQPR+f0FT0y+KlndFuGJeu7qpLtgRDYvBbP6739+776qWYPCfSo63xsZtdc+VS9fVsTjr7CygpFgNma/+HD53/47t0zB7Jr3F3JY8Dh8vy5kWusZJL93j71GAMH4JgvdULGJjiFVbo56FReNYVUj995VSHsAiuvJI2Qz254pSV/4FvOJC7BaQdxyDz9hD8kIOlfRsqSQHimXzpXLZ3IgPJYhK4/AkB9FNz0+dgPWEX7EoVid7CmNbUdmF4KWNCBABIlB1BMjSUHWXhAZEBKqfAAQTxwatdtQ39Z76geHhEVgB/iHl2qWLGxLi0/Ma2YWp6KBaDOUwr1SEcHRQFQ1xmZjnJJWV4FgJxD9mgt7GjRu1W/B4xcES1oWXEVHLpON6Qc7CVk2kA6+7judBB6P7GsRmSLJ6dczG5W86kuwVp/6dHxUnbXR0tISPB79MmzaHmHXQuIP1ZfNs/+g42zswxHb3D7K9QyPsUDrj9OfL/GC2eFf3rl2FkblzldXYsTzm+6sxF2YgXLjRMyw8X03X4HcmTb8QASJABKjCS38GiAAROGICu3btgrZUpVTciXec2BEf+MZ2hPTzq7X1OaQBjOGWerbisIrroHsX86KG5xXhecX5WQjBA4td95XC8o2d6Rjt3dXV5aISOr35QhC5tS+bqIoif9eNMB5K6awdO7wAIRwwEbmGZXpolcwr0wcdw+9qPDKg61txES88lC/XcPRLs11Px6I0Pl4y/KgxFTZmeZ5xKF+c2JvO/D2O4YVDh3zLBby+T2GhnSoJO1h+lwjI4BK8vp1tQWBDJwnfY3jt6K2JABH4IwlQhfePBEeHEYHjkABXQk7NGwuezgkGg/+lfka192hX9qbfb8L1JJqxuUxADDquPFWdDz28Hi+gUuo6rh3iPFlbH1uhnofArNq/z/xPCUhyqDieGjdTEWRJDcVqbFjYZlWUPQOvoUVFEk/pbOsxFY2qis57e3v3YVybSpaT3p8rGbvTOb4nnZUvYXHgs6MT8vHhtHx8YHz0hZH0ed279vWgus53/tafO64WuMECIYNYoBfXhRo3Q3e4qamq32gjAkSACFQPgar9B6J6ENFIiAAR2IKuWaAgcas+dsEFF3wN1coDd9xxxz14DgmuKt716G5bIcjwjg6qjlnTgi5E8FdMCIQy4Ef8V/T9ozYqvIy3GoYScG+FTVcVXlWZ1kEs4H8FVM62QfT6CQnwJC85g6mwBr8V8rEUj4qZ8kvfhwi4d6DS+6O+fHF8dzrPnx1Js9F8STa5juVmMhc9vvfAdmXRUNV1BXlqUDquDbppOEyDeE8J7ov3t8JFoDESASJwfBIgwXt8XneaNRE4YgJK7HaiRawSu6VS6WsQuE/ceeedtysRhDc5pmKzLN2A6ujFIa5wu9/2B+3yfNF2HdN2OaqLqJC6C494MrOz4/QHgj7LdU34dUUQ4zAda40ajtK6EJywOiAdATxXJGZskOraaQMDA3uHBgevqNX1JWO53KfbbNP74rwmd8uittDn21K1WyCMsVjt1ZtUlga4M5iBdImgRPcJ2ogAESACVUyA/pKq4otDQyMCs03g1WIX47kXYvdnSuxO2xuOxRh3TRUSccf86YqDPFh4XLG6a/E17UoTOs8jzaCAWC1dCceQxxarMTSsX38sK6JvZpqQ5f7WB09sCfYFzUBVNMpEnXq2wrxgGeIR3l5UfaWsNdiksJ885lh/VRYV9e+AvhsZusVc7t8mSmZhtFDUYbBgFcnP74QIf43Ob0IJXjVmVeGNaB5VeI/1laL3JwJE4E0RIMH7pvDRwUTg7UsAovZ3KruY6b2//OUvbzvWYlcRHZkSvFjQNVZWfl08IAZrVpTr1EJbUfQ8TVV+VbvhOiHK6phq3qZKvAasC1yJd4GqaJxP2gAqruiBgFfeXi8A7VgnY+1Tc5kpAa8qvc4HUe3Fd3vYsl8eQLc7Vb1NGNo6PKdteOgh37uNn6er1XuRLpG3bE8PQPCmK/JkvMYaRqr2Q4caHm1EgAgcxwRI8B7HF5+mTgReh4C/QA3i1rcxYL8ZE7uvHFPFleGSf7vfZWEurPHJF9FRWL5kQvAqb2/F83yxNdrdPS3GXvkWVfVz2ZPit00mlB0ZSpN7PcrbC4HponwdSjLhC14swpspweszmvqQYU3Y9qN9hRIrmJaX0PRlX2mvbVY+bcBV45muVmcd6dk2urDpEO9hyfxqtf9G9IUIEAEiUIUESPBW4UWhIRGB2SSABAZV6fMXqCnPLn6ecbHbPQUAFd6eEqq7SGrwYF+IOGFvLl5yIXLzJVgdmN+dbFJswWdarYIXdWhfLOYQaLEXIpFJpE6UpppMQPUaytKg0ihUekOTECrTdsY3MPfFLHg/dAi5vNlyxa4PBKKLalKnqcF0bZxMwZiCrOPDCEebYd/SEOWSLA0zfsXohESACLwRAiR43wgt2pcIvM0JoJInCoVC+cILL4xA7F6H6c642J1C7IsvlV07dbvfr36mpDZHvV6C6sLCNb9ZQ3RS6M5oNXRqjEf8beuk4HXQ3Syrmk+olWpImJiuitqo/EoLsWQqbyLoeX4axRG/+dHbUVkb1La9v2wWBwpFIyIEi2v8HPXkq+wKEpYGP2JNQ5U9KCX9W6Ig0UYEiEDVEqC/pKr20tDAiMDMEkBcFvQuL8fj8RjO/I943DNTnt3XmWlYCV4LvlcV5ZVimu/XRabBE340mV8tZSsuYKwGihf13pm1AbzOuF/zpYrHNGVfQN4uFq3xyUqu676IxWw5NJ8wgihS5xk7XR3cMPNzma5ED43Yzh50YBNYNMhShr4Gw+EbNmyYFsRqeBzVd6EqvCojOYRir3qSMdSJaSMCRIAIVCEBErxVeFFoSERgNgh0dHR4qOqmbNv+Om693zfLYnfqzjnLFB3PsSxEkAFKkUk/kQHCERVeJDdAcIWZ1BagZjobzN7oOSvS1RBNhtG6qiqqFuCpyrRT8aSnvL3KD5tA57U3+r5Ha/8NkwvX3AnbebQXtoZ8xWYJoS/bkkikOKLp1EUBejXmHAIl9sJp4lfZC1P2jA3dVWsrOVqI6H2IABF4ixIgwfsWvXA0bCJwtAnk8/mKaZphCN6fQ+z+YibSGF5nDtOLo/aq6mfZcZAGwFhYsEXqmArn6bzyvULxBiULtcUCber5rZNiTP1YVdtvYtYYQ8yai5g1j5mut+IL8CVjoHZJSpXgwDQlePlke+ENszAD1Gd97jnl4y2ZLGtWLOTz1q9piK9Sw+lCQsbWKZGO6m9W+Y4ZKrwBT07bM6Y/qGBv2ogAESAC1UOABG/1XAsaCRGYLQK+SEGurvXNb37zow888MBtGIg4ljm7RzLRKeWkYUGXUK13lRis46KkjsXt9KdQVWTww3oRxsOtkrer51dUqeCdSkBgtuRjytKgmkxEJNM7JiuqjiVlxq/8Yo5YkOeLekxnNsTj9Dl39JYrpcFSWY9jDWNSiPWK7yuzjm1k96r2wsqeEZa0aE3xoY0IEIHqJUCCt3qvDY2MCMw4gVNPPVUJSvX3wiv9mjM+jlec0Ct5nlPxs2v92+f16jXUFUN5VBaVSFT3/1s5N19xTNX+CKEeQGXXb9gQ4sxNQzRisFizxvaquaj0Bgj72RS8GIFvWRgdduy+fvh4PQwuqmlnY5x8tLFbTlerbYlcMswFgb0MdhNlc6CNCBABIlC1BEjwVu2loYERgVkhoIRLNYhdee2k8M5iMLuUDUCJQQxttU9FiENFV5rI6BUh1CTzUjtFPT8LC7384fyhL7AK+BtEojVd4UVLXvX3rxK8qtuaMJHSoKqlcTm1mM0/Yua/bJrkXhlzvMd7yybLwysd18SqaxIstamLucunxK3lyYASvCpxIggP8syPlM5IBIgAEThyAiR4j5wV7UkEjgcC07e0q2WusiI9V3VbUxkMNWg+4Q/MtvtL0jOLjqupTl8p5tWq5zf4L1blF58r/MhPliHcPce14UlO1EaDcGEwZiKRTD2v/LBBP2xi9uYAn65frUVniYcOlitYuGZZtZqeOjNct2pqVP6/G7giPYhZw/o7j8F2smQpY3G8Pt2gYvYmQGcmAkSACLwGARK8rwGFniICRKB6CJSxNk1Fk0kIqwCEIUamBJlWkNJVz6s2vchR80Va9Yz6/x/J1ABtE7m1yvuKyjSPS9+RAZUod1Zc10FSA7J62YL1jMWwP2Y2K/PyJv9h8B4/VLGtwbKp1Wg6T+lig5pVrr1dU989znsg4H3Ba0gvORcNNDBmDJk2IkAEiED1ESDBW33XhEZEBIgACDw0eWudOVJuU8JWiUTX85ZvVHG8jKWREvBiGVYHNKdguL3uZ9fi+Wq/tS5M3P1XVVHVVS3OmZ9fi7K1CXmrJ7GQLSLl/A70n8BcZmvzUL1V+vzAkG0f6IPgVXXbGk2coQbkhPp8xrZ0YWnAnpPi3VlIYne2rhedlwgQgSMg4PvHjmA/2oUIEAEiMCsEbC4t1XpXicQwPK9zJoWwLEJgIakB/gDJkqguzsrg3vBJdVSrpVTthQOCM8GFEu+Ykbgv7XnlDKwbYclGFzCWVdXd2aqYTvl4rXEpn0Baw9ICxlvD9RM+zFi0toeV1ZDRWdjDBw184ECZ2pO8cXaq0WootBEBIkAE/iABErx/EBHtQASIwGwSwOI0Q3VVU9oWNgDVdxdlRbQXRu8JND4QWPylFnuFEW2gqqKTHt/ZHPBrn9tTJVPJnF221LMl100GuYaEA+8kPH3rmOv+6qvqkZ34zdGdv/lp5n+Y9vFOwMd7oGJ9NFuxbVR42y5JRpZdnC3tVHNRiRMqO1ji2mjSkxHkCc/8SOmMRIAIEIEjI0CWhiPjRHsRASIwwwS6p26RV6R4qghhhR4TDrqTxVvD+mI1lIDn/UXZludMOM7Zlq5f1APBBSFW7R5Sx4RoR7oBV7nCKeXHwPZlTGcLqtcYvP7gVHKDen4WN9/Hi/L5rw/admWoUhEpoYkmTTtLjUn9w2Fz/jIWFDrwVlsGkslg7G1Tr22lSq/CQBsRIAJVRoAqvFV2QWg4RIAI/JaAqiRKnY8hi5fBBiCxOM1o9WRC7dHL2MGr3crBqzOV3x5Q/T+JMvywaoGaQEOHmOR+VXSrWgMGe2xn9Ywfo/S3vn7HG+0zK+08EEAMg7Yez35DvYJFawcF2iPjgwcWDcpAkrF5eHrPChK8Cg9tRIAIVBkBqvBW2QWh4RABIvAqAlKGSqiGqsxXA/XQlDuZU4uqqH4L0hogzNRDaeOq3uDDUGO00VVNC7mejoqo7knZUaWDlh+c7AJXGvO8HYcrFVYC/6jgpykfr1+Wdt0XbU92jtvOx01PngNLw5NqLvD/+i9X6bxoWESACBynBKjCe5xeeJo2EXgLEICOxea6fSXO8mXHCQSEcDKMrcazD2/Al3dN+XnxY7Vvai5K8OYhdD+Tcd2F8L/aJSnvnxo49HB1bV1THyLGGes+YNnvz9q2XS+05ksMY9lNtr0T1ZLR61136/X5/KsHPnndXv0s/U4EiAARmEUCJHhnET6dmggQgdclIJVChHoagO81buDWeVgIlpJyjjpqg/ryFtowF18I7vC8G/+0jKAD9Zjapl+b/r1Kvk/6eDXtkf2O645UbNaiCa2e87Mwvp3Xw14ygDltwC94qLkpWwaJXYCgjQgQgeojQIK3+q4JjYgIEIEpAlPqyYIN4K/GXS/kMJenhLgdVV+1x1vy1jlyhLU/R/V0AyaATYnEap3HJH7b3jOgaYO9lUrLqkjEqRHibIz7G5snx+52+tOgL0SACBABIkAEiAARIAJE4C1IQIlzNexmIX5xQygkZVOTHIyEd38GEXBQw9DqtBEBIkAE3hoEaNHaW+M60SiJwHFNYP1UXBdElr5lqgPbcQ1khiefY+y23Y5zU8G2rw5wsRGeBmeGh0CnIwJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBCBP5YAn4qAU80x6M7gH0uRjiMCRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEAEiQASIABEgAkSACBABIkAEiAARIAJEgAgQASJABIgAESACRIAIEIG3NIH/B7bHNqpQRRifAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<PIL.Image.Image image mode=RGBA size=700x300>"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from typing import Dict, List, Tuple\n",
    "\n",
    "import matplotlib.pyplot as plt\n",
    "import scm.plams as plams\n",
    "from scm.base import ChemicalSystem\n",
    "from scm.plams import Settings, view\n",
    "from scm.plams.recipes.counterpoise import CounterpoiseEnergyJob\n",
    "from scm.utils.conversions import chemsys_to_plams_molecule\n",
    "\n",
    "ab_system = ChemicalSystem(\n",
    "    \"\"\"\n",
    "System\n",
    "    Atoms\n",
    "        H   4.194124   2.321011   9.990591\n",
    "        O   3.277064   2.369351   9.662700\n",
    "        H   3.215647   3.302746   9.160403\n",
    "        O   3.121052   4.371977   8.152437\n",
    "        O   3.157940   9.599596   8.160532\n",
    "        O  -0.155188   4.377678   8.158283\n",
    "        O  -0.092407   9.610401   8.171480\n",
    "        Zn  3.135686   1.046635   8.068212\n",
    "        Zn  3.127411   6.251538   7.820969\n",
    "        Zn -0.123253   6.228628   7.789648\n",
    "        Zn -0.108077   0.986306   7.727521\n",
    "        Zn  1.540445   9.053723   7.248234\n",
    "        Zn  4.781682   9.052463   7.248636\n",
    "        Zn  1.454711   3.779398   7.208485\n",
    "        Zn  4.770754   3.793547   7.144487\n",
    "        O   1.537161   1.798390   7.074007\n",
    "        O   4.752079   1.794734   7.054472\n",
    "        O   1.506053   7.034856   7.081453\n",
    "        O   4.747135   7.034630   7.078047\n",
    "        O   1.504136   4.449837   5.256964\n",
    "        O   1.510233   9.664240   5.251216\n",
    "        O   4.745874   4.467809   5.219866\n",
    "        O   4.765413   9.642604   5.251871\n",
    "        Zn  1.527610   1.213182   5.119746\n",
    "        Zn  4.738978   1.206559   5.113718\n",
    "        Zn  1.504703   6.460518   5.092362\n",
    "        Zn  4.750778   6.477069   5.090083\n",
    "        Zn -0.110739   9.072159   4.273393\n",
    "        Zn  3.135224   9.070253   4.279334\n",
    "        Zn -0.117138   3.823484   4.274300\n",
    "        Zn  3.118895   3.832597   4.260451\n",
    "        O  -0.115211   1.822988   4.141734\n",
    "        O   3.131372   1.834356   4.108967\n",
    "        O  -0.117871   7.064012   4.120018\n",
    "        O   3.122850   7.064255   4.112580\n",
    "        O  -0.120166   4.440964   2.316795\n",
    "        O  -0.116391   9.670504   2.304567\n",
    "        O   3.130312   4.443558   2.298592\n",
    "        O   3.130933   9.674750   2.306875\n",
    "        Zn -0.116763   1.207088   2.147541\n",
    "        Zn  3.130704   1.213453   2.119764\n",
    "        Zn -0.119363   6.435899   2.118992\n",
    "        Zn  3.130940   6.437444   2.109042\n",
    "        Zn  4.758737   3.618851   1.589370\n",
    "        Zn  4.755815   8.858895   1.584735\n",
    "        Zn  1.508461   8.855081   1.581480\n",
    "        Zn  1.502167   3.625979   1.585735\n",
    "        O   4.763512   1.779351   1.204246\n",
    "        O   1.500332   1.785698   1.205474\n",
    "        O   4.757821   7.015026   1.196018\n",
    "        O   1.503517   7.011011   1.191741\n",
    "    End\n",
    "    Lattice\n",
    "        6.50000095   0.0   0.0\n",
    "        0.0   10.46000194   0.0\n",
    "    End\n",
    "End\n",
    "\"\"\"\n",
    ")\n",
    "\n",
    "view(ab_system, width=700, height=300, direction=\"tilt_x\", guess_bonds=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "58090102",
   "metadata": {},
   "source": [
    "## Run counterpoise jobs for each basis set\n",
    "\n",
    "This section can take significant time. For production use, run on appropriate compute resources.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "aa97a3ad",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-03-13T11:35:34.914708Z",
     "iopub.status.busy": "2026-03-13T11:35:34.914463Z",
     "iopub.status.idle": "2026-03-13T19:23:50.226699Z",
     "shell.execute_reply": "2026-03-13T19:23:50.225829Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[13.03|12:35:35] JOB DZ_counterpoise STARTED\n",
      "[13.03|12:35:35] JOB DZ_counterpoise RUNNING\n",
      "[13.03|12:35:35] JOB DZ_counterpoise/AB_geo_AB_basis_AB STARTED\n",
      "[13.03|12:35:35] JOB DZ_counterpoise/AB_geo_AB_basis_AB RUNNING\n",
      "[13.03|12:56:26] JOB DZ_counterpoise/AB_geo_AB_basis_AB FINISHED\n",
      "[13.03|12:56:26] JOB DZ_counterpoise/AB_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|12:56:26] JOB DZ_counterpoise/A_geo_AB_basis_AB STARTED\n",
      "[13.03|12:56:26] JOB DZ_counterpoise/A_geo_AB_basis_AB RUNNING\n",
      "[13.03|12:57:05] JOB DZ_counterpoise/A_geo_AB_basis_AB FINISHED\n",
      "[13.03|12:57:05] JOB DZ_counterpoise/A_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|12:57:05] JOB DZ_counterpoise/A_geo_AB_basis_A STARTED\n",
      "[13.03|12:57:05] JOB DZ_counterpoise/A_geo_AB_basis_A RUNNING\n",
      "[13.03|12:57:12] JOB DZ_counterpoise/A_geo_AB_basis_A FINISHED\n",
      "[13.03|12:57:12] JOB DZ_counterpoise/A_geo_AB_basis_A SUCCESSFUL\n",
      "[13.03|12:57:12] JOB DZ_counterpoise/B_geo_AB_basis_AB STARTED\n",
      "[13.03|12:57:12] JOB DZ_counterpoise/B_geo_AB_basis_AB RUNNING\n",
      "[13.03|12:58:14] JOB DZ_counterpoise/B_geo_AB_basis_AB FINISHED\n",
      "[13.03|12:58:14] JOB DZ_counterpoise/B_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|12:58:14] JOB DZ_counterpoise/B_geo_AB_basis_B STARTED\n",
      "[13.03|12:58:14] JOB DZ_counterpoise/B_geo_AB_basis_B RUNNING\n",
      "[13.03|12:59:11] JOB DZ_counterpoise/B_geo_AB_basis_B FINISHED\n",
      "[13.03|12:59:12] JOB DZ_counterpoise/B_geo_AB_basis_B SUCCESSFUL\n",
      "[13.03|12:59:12] JOB DZ_counterpoise/A_geo_A_basis_A STARTED\n",
      "[13.03|12:59:12] JOB DZ_counterpoise/A_geo_A_basis_A RUNNING\n",
      "[13.03|13:02:44] JOB DZ_counterpoise/A_geo_A_basis_A FINISHED\n",
      "[13.03|13:02:44] JOB DZ_counterpoise/A_geo_A_basis_A SUCCESSFUL\n",
      "[13.03|13:02:44] JOB DZ_counterpoise/B_geo_B_basis_B STARTED\n",
      "[13.03|13:02:44] JOB DZ_counterpoise/B_geo_B_basis_B RUNNING\n",
      "[13.03|13:25:31] JOB DZ_counterpoise/B_geo_B_basis_B FINISHED\n",
      "[13.03|13:25:31] JOB DZ_counterpoise/B_geo_B_basis_B SUCCESSFUL\n",
      "[13.03|13:25:31] JOB DZ_counterpoise FINISHED\n",
      "[13.03|13:25:31] JOB DZ_counterpoise SUCCESSFUL\n",
      "[13.03|13:25:31] JOB DZP_counterpoise STARTED\n",
      "[13.03|13:25:31] JOB DZP_counterpoise RUNNING\n",
      "[13.03|13:25:31] JOB DZP_counterpoise/AB_geo_AB_basis_AB STARTED\n",
      "[13.03|13:25:31] JOB DZP_counterpoise/AB_geo_AB_basis_AB RUNNING\n",
      "[13.03|13:51:00] JOB DZP_counterpoise/AB_geo_AB_basis_AB FINISHED\n",
      "[13.03|13:51:00] JOB DZP_counterpoise/AB_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|13:51:00] JOB DZP_counterpoise/A_geo_AB_basis_AB STARTED\n",
      "[13.03|13:51:00] JOB DZP_counterpoise/A_geo_AB_basis_AB RUNNING\n",
      "[13.03|13:51:50] JOB DZP_counterpoise/A_geo_AB_basis_AB FINISHED\n",
      "[13.03|13:51:50] JOB DZP_counterpoise/A_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|13:51:50] JOB DZP_counterpoise/A_geo_AB_basis_A STARTED\n",
      "[13.03|13:51:50] JOB DZP_counterpoise/A_geo_AB_basis_A RUNNING\n",
      "[13.03|13:51:57] JOB DZP_counterpoise/A_geo_AB_basis_A FINISHED\n",
      "[13.03|13:51:57] JOB DZP_counterpoise/A_geo_AB_basis_A SUCCESSFUL\n",
      "[13.03|13:51:57] JOB DZP_counterpoise/B_geo_AB_basis_AB STARTED\n",
      "[13.03|13:51:57] JOB DZP_counterpoise/B_geo_AB_basis_AB RUNNING\n",
      "[13.03|13:53:12] JOB DZP_counterpoise/B_geo_AB_basis_AB FINISHED\n",
      "[13.03|13:53:12] JOB DZP_counterpoise/B_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|13:53:12] JOB DZP_counterpoise/B_geo_AB_basis_B STARTED\n",
      "[13.03|13:53:12] JOB DZP_counterpoise/B_geo_AB_basis_B RUNNING\n",
      "[13.03|13:54:19] JOB DZP_counterpoise/B_geo_AB_basis_B FINISHED\n",
      "[13.03|13:54:19] JOB DZP_counterpoise/B_geo_AB_basis_B SUCCESSFUL\n",
      "[13.03|13:54:19] JOB DZP_counterpoise/A_geo_A_basis_A STARTED\n",
      "[13.03|13:54:19] JOB DZP_counterpoise/A_geo_A_basis_A RUNNING\n",
      "[13.03|13:58:34] JOB DZP_counterpoise/A_geo_A_basis_A FINISHED\n",
      "[13.03|13:58:34] JOB DZP_counterpoise/A_geo_A_basis_A SUCCESSFUL\n",
      "[13.03|13:58:34] JOB DZP_counterpoise/B_geo_B_basis_B STARTED\n",
      "[13.03|13:58:34] JOB DZP_counterpoise/B_geo_B_basis_B RUNNING\n",
      "[13.03|14:26:31] JOB DZP_counterpoise/B_geo_B_basis_B FINISHED\n",
      "[13.03|14:26:31] JOB DZP_counterpoise/B_geo_B_basis_B SUCCESSFUL\n",
      "[13.03|14:26:31] JOB DZP_counterpoise FINISHED\n",
      "[13.03|14:26:31] JOB DZP_counterpoise SUCCESSFUL\n",
      "[13.03|14:26:31] JOB TZP_counterpoise STARTED\n",
      "[13.03|14:26:31] JOB TZP_counterpoise RUNNING\n",
      "[13.03|14:26:31] JOB TZP_counterpoise/AB_geo_AB_basis_AB STARTED\n",
      "[13.03|14:26:31] JOB TZP_counterpoise/AB_geo_AB_basis_AB RUNNING\n",
      "[13.03|14:54:12] JOB TZP_counterpoise/AB_geo_AB_basis_AB FINISHED\n",
      "[13.03|14:54:12] JOB TZP_counterpoise/AB_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|14:54:12] JOB TZP_counterpoise/A_geo_AB_basis_AB STARTED\n",
      "[13.03|14:54:12] JOB TZP_counterpoise/A_geo_AB_basis_AB RUNNING\n",
      "[13.03|14:55:04] JOB TZP_counterpoise/A_geo_AB_basis_AB FINISHED\n",
      "[13.03|14:55:05] JOB TZP_counterpoise/A_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|14:55:05] JOB TZP_counterpoise/A_geo_AB_basis_A STARTED\n",
      "[13.03|14:55:05] JOB TZP_counterpoise/A_geo_AB_basis_A RUNNING\n",
      "[13.03|14:55:12] JOB TZP_counterpoise/A_geo_AB_basis_A FINISHED\n",
      "[13.03|14:55:12] JOB TZP_counterpoise/A_geo_AB_basis_A SUCCESSFUL\n",
      "[13.03|14:55:12] JOB TZP_counterpoise/B_geo_AB_basis_AB STARTED\n",
      "[13.03|14:55:12] JOB TZP_counterpoise/B_geo_AB_basis_AB RUNNING\n",
      "[13.03|14:56:33] JOB TZP_counterpoise/B_geo_AB_basis_AB FINISHED\n",
      "[13.03|14:56:33] JOB TZP_counterpoise/B_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|14:56:33] JOB TZP_counterpoise/B_geo_AB_basis_B STARTED\n",
      "[13.03|14:56:33] JOB TZP_counterpoise/B_geo_AB_basis_B RUNNING\n",
      "[13.03|14:57:46] JOB TZP_counterpoise/B_geo_AB_basis_B FINISHED\n",
      "[13.03|14:57:46] JOB TZP_counterpoise/B_geo_AB_basis_B SUCCESSFUL\n",
      "[13.03|14:57:46] JOB TZP_counterpoise/A_geo_A_basis_A STARTED\n",
      "[13.03|14:57:46] JOB TZP_counterpoise/A_geo_A_basis_A RUNNING\n",
      "[13.03|15:00:44] JOB TZP_counterpoise/A_geo_A_basis_A FINISHED\n",
      "[13.03|15:00:44] JOB TZP_counterpoise/A_geo_A_basis_A SUCCESSFUL\n",
      "[13.03|15:00:44] JOB TZP_counterpoise/B_geo_B_basis_B STARTED\n",
      "[13.03|15:00:44] JOB TZP_counterpoise/B_geo_B_basis_B RUNNING\n",
      "[13.03|15:30:10] JOB TZP_counterpoise/B_geo_B_basis_B FINISHED\n",
      "[13.03|15:30:10] JOB TZP_counterpoise/B_geo_B_basis_B SUCCESSFUL\n",
      "[13.03|15:30:10] JOB TZP_counterpoise FINISHED\n",
      "[13.03|15:30:10] JOB TZP_counterpoise SUCCESSFUL\n",
      "[13.03|15:30:10] JOB TZ2P_counterpoise STARTED\n",
      "[13.03|15:30:10] JOB TZ2P_counterpoise RUNNING\n",
      "[13.03|15:30:10] JOB TZ2P_counterpoise/AB_geo_AB_basis_AB STARTED\n",
      "[13.03|15:30:10] JOB TZ2P_counterpoise/AB_geo_AB_basis_AB RUNNING\n",
      "[13.03|16:07:32] JOB TZ2P_counterpoise/AB_geo_AB_basis_AB FINISHED\n",
      "[13.03|16:07:32] JOB TZ2P_counterpoise/AB_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|16:07:32] JOB TZ2P_counterpoise/A_geo_AB_basis_AB STARTED\n",
      "[13.03|16:07:32] JOB TZ2P_counterpoise/A_geo_AB_basis_AB RUNNING\n",
      "[13.03|16:09:30] JOB TZ2P_counterpoise/A_geo_AB_basis_AB FINISHED\n",
      "[13.03|16:09:30] JOB TZ2P_counterpoise/A_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|16:09:30] JOB TZ2P_counterpoise/A_geo_AB_basis_A STARTED\n",
      "[13.03|16:09:30] JOB TZ2P_counterpoise/A_geo_AB_basis_A RUNNING\n",
      "[13.03|16:09:37] JOB TZ2P_counterpoise/A_geo_AB_basis_A FINISHED\n",
      "[13.03|16:09:37] JOB TZ2P_counterpoise/A_geo_AB_basis_A SUCCESSFUL\n",
      "[13.03|16:09:37] JOB TZ2P_counterpoise/B_geo_AB_basis_AB STARTED\n",
      "[13.03|16:09:37] JOB TZ2P_counterpoise/B_geo_AB_basis_AB RUNNING\n",
      "[13.03|16:11:24] JOB TZ2P_counterpoise/B_geo_AB_basis_AB FINISHED\n",
      "[13.03|16:11:24] JOB TZ2P_counterpoise/B_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|16:11:24] JOB TZ2P_counterpoise/B_geo_AB_basis_B STARTED\n",
      "[13.03|16:11:24] JOB TZ2P_counterpoise/B_geo_AB_basis_B RUNNING\n",
      "[13.03|16:13:01] JOB TZ2P_counterpoise/B_geo_AB_basis_B FINISHED\n",
      "[13.03|16:13:01] JOB TZ2P_counterpoise/B_geo_AB_basis_B SUCCESSFUL\n",
      "[13.03|16:13:01] JOB TZ2P_counterpoise/A_geo_A_basis_A STARTED\n",
      "[13.03|16:13:01] JOB TZ2P_counterpoise/A_geo_A_basis_A RUNNING\n",
      "[13.03|16:19:13] JOB TZ2P_counterpoise/A_geo_A_basis_A FINISHED\n",
      "[13.03|16:19:13] JOB TZ2P_counterpoise/A_geo_A_basis_A SUCCESSFUL\n",
      "[13.03|16:19:13] JOB TZ2P_counterpoise/B_geo_B_basis_B STARTED\n",
      "[13.03|16:19:13] JOB TZ2P_counterpoise/B_geo_B_basis_B RUNNING\n",
      "[13.03|16:58:56] JOB TZ2P_counterpoise/B_geo_B_basis_B FINISHED\n",
      "[13.03|16:58:56] JOB TZ2P_counterpoise/B_geo_B_basis_B SUCCESSFUL\n",
      "[13.03|16:58:56] JOB TZ2P_counterpoise FINISHED\n",
      "[13.03|16:58:56] JOB TZ2P_counterpoise SUCCESSFUL\n",
      "[13.03|16:58:56] JOB QZ4P_counterpoise STARTED\n",
      "[13.03|16:58:56] JOB QZ4P_counterpoise RUNNING\n",
      "[13.03|16:58:56] JOB QZ4P_counterpoise/AB_geo_AB_basis_AB STARTED\n",
      "[13.03|16:58:56] JOB QZ4P_counterpoise/AB_geo_AB_basis_AB RUNNING\n",
      "[13.03|18:28:31] JOB QZ4P_counterpoise/AB_geo_AB_basis_AB FINISHED\n",
      "[13.03|18:28:31] JOB QZ4P_counterpoise/AB_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|18:28:31] JOB QZ4P_counterpoise/A_geo_AB_basis_AB STARTED\n",
      "[13.03|18:28:31] JOB QZ4P_counterpoise/A_geo_AB_basis_AB RUNNING\n",
      "[13.03|18:30:44] JOB QZ4P_counterpoise/A_geo_AB_basis_AB FINISHED\n",
      "[13.03|18:30:45] JOB QZ4P_counterpoise/A_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|18:30:45] JOB QZ4P_counterpoise/A_geo_AB_basis_A STARTED\n",
      "[13.03|18:30:45] JOB QZ4P_counterpoise/A_geo_AB_basis_A RUNNING\n",
      "[13.03|18:30:52] JOB QZ4P_counterpoise/A_geo_AB_basis_A FINISHED\n",
      "[13.03|18:30:52] JOB QZ4P_counterpoise/A_geo_AB_basis_A SUCCESSFUL\n",
      "[13.03|18:30:52] JOB QZ4P_counterpoise/B_geo_AB_basis_AB STARTED\n",
      "[13.03|18:30:52] JOB QZ4P_counterpoise/B_geo_AB_basis_AB RUNNING\n",
      "[13.03|18:34:38] JOB QZ4P_counterpoise/B_geo_AB_basis_AB FINISHED\n",
      "[13.03|18:34:38] JOB QZ4P_counterpoise/B_geo_AB_basis_AB SUCCESSFUL\n",
      "[13.03|18:34:38] JOB QZ4P_counterpoise/B_geo_AB_basis_B STARTED\n",
      "[13.03|18:34:38] JOB QZ4P_counterpoise/B_geo_AB_basis_B RUNNING\n",
      "[13.03|18:38:03] JOB QZ4P_counterpoise/B_geo_AB_basis_B FINISHED\n",
      "[13.03|18:38:03] JOB QZ4P_counterpoise/B_geo_AB_basis_B SUCCESSFUL\n",
      "[13.03|18:38:03] JOB QZ4P_counterpoise/A_geo_A_basis_A STARTED\n",
      "[13.03|18:38:03] JOB QZ4P_counterpoise/A_geo_A_basis_A RUNNING\n",
      "[13.03|18:40:15] JOB QZ4P_counterpoise/A_geo_A_basis_A FINISHED\n",
      "[13.03|18:40:15] JOB QZ4P_counterpoise/A_geo_A_basis_A SUCCESSFUL\n",
      "[13.03|18:40:15] JOB QZ4P_counterpoise/B_geo_B_basis_B STARTED\n",
      "[13.03|18:40:15] JOB QZ4P_counterpoise/B_geo_B_basis_B RUNNING\n",
      "[13.03|20:23:50] JOB QZ4P_counterpoise/B_geo_B_basis_B FINISHED\n",
      "[13.03|20:23:50] JOB QZ4P_counterpoise/B_geo_B_basis_B SUCCESSFUL\n",
      "[13.03|20:23:50] JOB QZ4P_counterpoise FINISHED\n",
      "[13.03|20:23:50] JOB QZ4P_counterpoise SUCCESSFUL\n"
     ]
    }
   ],
   "source": [
    "ab_mol = chemsys_to_plams_molecule(ab_system)\n",
    "ids_state_A: List[int] = [1, 2, 3]  # H2O fragment\n",
    "basis_sets: List[str] = [\"DZ\", \"DZP\", \"TZP\", \"TZ2P\", \"QZ4P\"]\n",
    "\n",
    "base_settings = Settings()\n",
    "base_settings.input.ams.Task = \"GeometryOptimization\"\n",
    "base_settings.input.band.KSpace.Quality = \"Basic\"\n",
    "base_settings.input.band.XC.GGA = \"PBE\"\n",
    "\n",
    "jobs: Dict[str, CounterpoiseEnergyJob] = {}\n",
    "for basis in basis_sets:\n",
    "    settings = base_settings.copy()\n",
    "    settings.input.band.Basis.Type = basis\n",
    "    job = CounterpoiseEnergyJob(\n",
    "        name=f\"{basis}_counterpoise\",\n",
    "        AB=ab_mol.copy(),\n",
    "        ids_state_A=ids_state_A,\n",
    "        settings_state_AB=settings,\n",
    "    )\n",
    "    job.run()\n",
    "    jobs[basis] = job"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "a74da19e",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-03-13T19:23:50.232460Z",
     "iopub.status.busy": "2026-03-13T19:23:50.232263Z",
     "iopub.status.idle": "2026-03-13T19:23:50.422974Z",
     "shell.execute_reply": "2026-03-13T19:23:50.422090Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Basis  Ebind_uncorrected  Ebind_corrected  Counterpoise_correction\n",
      "DZ                -37.83           -34.13                      3.70\n",
      "DZP               -27.05           -24.05                      2.99\n",
      "TZP               -24.50           -23.19                      1.31\n",
      "TZ2P              -23.74           -22.59                      1.14\n",
      "QZ4P              -22.43           -22.28                      0.14\n"
     ]
    }
   ],
   "source": [
    "rows: List[Tuple[str, float, float, float]] = []\n",
    "for basis in basis_sets:\n",
    "    energies = jobs[basis].results.get_all_energies(unit=\"kcal/mol\")\n",
    "    rows.append((basis, energies[\"Ebind_raw\"], energies[\"Ebind_cp\"], -energies[\"BSSE_tot\"]))\n",
    "\n",
    "print(\"Basis  Ebind_uncorrected  Ebind_corrected  Counterpoise_correction\")\n",
    "for basis, raw, corrected, cp_corr in rows:\n",
    "    print(f\"{basis:5s}  {raw:17.2f}  {corrected:15.2f}  {cp_corr:24.2f}\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "91018bea",
   "metadata": {},
   "source": [
    "## Plot adsorption energies\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "99809f23",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-03-13T19:23:50.427619Z",
     "iopub.status.busy": "2026-03-13T19:23:50.427452Z",
     "iopub.status.idle": "2026-03-13T19:23:50.621777Z",
     "shell.execute_reply": "2026-03-13T19:23:50.621110Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfAAAAEYCAYAAACju6QJAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAAsTAAALEwEAmpwYAABP0klEQVR4nO3dd3gVVfrA8e+bkJBCCSX0Kko3IATsgihiR6yruIoNsZe1gKzYVyyrrqgoa+GngKJIESsgYl2U3rvSe0kCJCHt/f0xk+QmpFyS3Dsp7+d57nPvnCnnvRnImzNz5hxRVYwxxhhTsYR4HYAxxhhjjp0lcGOMMaYCsgRujDHGVECWwI0xxpgKyBK4McYYUwFZAjfGGGMqIEvgxgSJiPQWka1exwEgIitEpLfXcZQHIvK8iNzvdRz5ici/ReQOr+Mw5ZclcFMpiMhGETk3X9kgEfnF/VxdRN4TkU0iclBEFovIBfm2jxGR0SKyU0SSRWSZiNwUzO8RCCIyVkSe9S1T1U6qOsejkIolIo+JyKF8r8MioiJyg5/HKPZ8ikgscAPwjrscLiKT3H9Pmv+PHHG8ICL73NcLIiJFxHC1iPzm1j+ngPVdRWSBu36BiHT1Wf0y8JiIhPvzfU3VYwncVBXVgC1AL6A28E/gUxFpBc4vbmAW0BI41d3mYWCkiDzoRcD+EJFqXscQCKr6L1Wt4fsCXgVWAp8Xt/8xnM9BwNeqmuJT9gtwPbCzgEMPBi4DugBxwCXA7UWEsh94DRhZSIzTgHFAHeD/gGnZCVtVdwCrgUuL+q6mClNVe9mrwr+AjcC5+coGAb8Usc9S4Ar38y3AbiA63zbXAIeAWoUc4z84fxgkAQuAM33WRQJjgQM4iedhYKvP+keBbcBBYA1wjlteHeeX/nb39RpQ3V3XG9jq7rsT+Min7DFgr/uzGOhuPxhIB9Lc7zE9/8/Lz/r+4f58dgA3FfEzbQJ8gZO41gO3+ax7EvgU+ND9ziuAeD/P74Xuz7hDvnP+kHseE4GJQMSxnE9gNnB9IXVuBXrnK/sNGOyzfAsw14/4bwXm5Cs7zz3/4lO2GTjfZ3k48IHX/7/sVT5f1gI3VZKINATa4iQRgL7AN6p6ON+mnwMROK24gswDugJ1gQnAZyIS4a57AmjjvvoBN/rU3w64G+ihqjXd9Rvd1cOBU9zjdgF64lwxyNbIra8lToLOLqsPNHXrGSMi7VR1DDAeeFGdluwlBXwHf+qr7R77FuBNEalTyM/jE5zE1wS4EviXiPTxWX+pu00MTqJ/o5Dj5HCvknyE88fAqnyrrwbOB1rjtIgHueX+ns8Tcf548lcnYInP8hK3rCQ6AUtV1Xc866X5jrcK55wYcxRL4KYymSoiCdkv4K2CNhKRMJyk9n+qutotro/TusxDVTNwWrX1CzqWqo5T1X2qmqGq/8ZpzbZzV18NPKeq+1V1C/C6z66Z7rYdRSRMVTeq6gZ33UDgaVXdrap7gKeAv/vsmwU8oapHNO+l38fdsh+Br9z6/VFcfenu+nRV/RqnBdsu/0FEpDlwOvCoqqaq6mLgXZx7zNl+UdWvVTUTJykXmZxEpDowCRivqhML2OR1Vd2uqvuB6Th/hID/5zMG52qAv2rgtPazJQI1iroPfgzHyj5eTZ/lg26MxhzFEripTC5T1ZjsF3Bn/g1EJAQncaThtICz7QUaF7B9NZxf9nsLqlBEHhKRVSKS6P7RUJvc5NAE5/J6tk3ZH1R1PXA/zmXl3SLyiYg08dlvU779mvgs71HV1HyhHMjX2sy/T1GKq2+fm/iyJeMkn4KOs19VfRPiJpyWezbf+8rJQEQx9/H/g/MHxD8KWZ//eNlx+Xs+D5A3YRbnEFDLZ7kWcEhVVUTe9ulw91gJjpV9PN+fX00g4RjiM1WIJXBTZbitpPeAhjj3vtN9Vs8CLhCR6Hy7XQEcAeYWcLwzgUdwWrp13D8aEoHs1tgOoLnPLi1891fVCap6Bs6lcAVecFdtd8t899vuu2sBX69Ovth99yluysHi6vPXdqCuiPgmxBY493mPmYj8Hefnf3W+c+UPf8/nUpxbKf5aQd6rBl3cMlR1iOZ2uvuXn8eKy9d6jyP3tg5AB/JesjcmhyVwU5WMxvmFeEm+S8/gtMq34tzDbiUiYSLSD+ey95Oqmv9SJzitowxgD1BNREaQt0X1KTBMROqISDPgnuwVItJORPq4l4hTgRScS+MAHwP/FJFYEakPjMDpqVycp9zHoM4ELgY+c8t3AccVsV9J68vDvU3wG/C8iESISBzOPfNjPpaIdMa5BTLQPe6x8vd8fo3zZIJv3dV9+jGEu98lO8l+CDwoIk3dKyb/wOmoWNj3CHWPVQ0IcY8V5q6eg3Mr5V63zuwrQrN9DtEL+OaYv72pEiyBmypBRFriPO7TFdjpc6lzIICqHgHOxbnk/TtOj+dXgOGq+lIhh/0O+BZYi3OpOJW8l8yfcsv/AmbgJJVs1XEeLdqLcxm4ATDMXfcsMB+ndbgMWOiWFWUnzuXg7Tj394f43N9/D+dee4KITC1g35LUV5hrgVZuHFNw7tXPKsFxHgSigcly9PPgxV6ePobz+SFwoYhE+pStwfmDqinOOU4h9wrFOzj32pcBy3H6GrxTRCh/d/cfDZzpfv6vG2MaziNpN+BcJr8Z5zZQGoCINAY6AlOL+76mapK8HSCNMRWNO9jIOFVt5nEoFZKI/AvYraqveR2LLxH5N7BBVQvsjGmMJXBjKjhL4MZUTXYJ3RhjjKmArAVujDHGVEDWAjfGGGMqoEoxEUL9+vW1VatWXodhjDHGlLkFCxbsVdXY/OWVIoG3atWK+fPnex2GMcYYU+ZEZFNB5XYJ3RhjjKmALIEbY4wxFZAnCVxEXhKR1SKyVESmiEiMW95XRBaIyDL3vU8xhzLGGGOqJK9a4DOBzqoahzMMZfYQkntxxqk+EWdO448K2d8YY4yp0jxJ4Ko6w2d6wrlAM7d8kapmz4K0Aoh0J3swxhhjjI/y0Av9ZmBiAeVXAAvdSQmOIiKDgcEALVq0KGgTY4wxJmimLtrGS9+tYXtCCk1iInm4XzsuO6lpwOoLWAIXkVlAowJWDVfVae42w3GmYxyfb99OOHMjn1fY8VV1DDAGID4+3oaTM8YY45mpi7YxbPIyUtIzAdiWkMKwycsAApbEA5bAVfXcotaLyCCcOYvPUZ/xXN15k6cAN6jqhkDFZ4wxxhyLjMwsUtIzSUnPJDUti+T0DFLSnOVnvlxJ38wfeST8U5rIXrZrfV7MuJqXvguveAm8KCJyPvAI0EtVk33KY3Dm1x2qqr96EZsxxpiKJytLc5JrdlIt6D05PZNUdzk5LZNUn/U5y/nWJadlkJqeRVpmVqH1XxryCyPD3iVK0gBoJnsZGfYuw5IAAvNAlVf3wN8AqgMzRQRgrqoOAe4GjgdGiMgId9vzVHW3N2EaY0z5Fex7riWlqhzJyCLZN6Fmf07PJCUtw33PcpNlQUnUbfmm+3x2E3JKWiZHMgpProUJDw0hIiyEyPBQosKrEREWSmRYCFHh1agTFe6Uh4USGR5KRFgoUeGhRIaFEuG+Zy9HhmTQZsJdRJGW5/hRksaw8M+A58voJ5mXJwlcVY8vpPxZ4Nkgh2OMMRVOWd1zVVXSMrPytlILbYkWkkTzJOOC349VaIgQ5SbLnMTpJs1GtcKccje5RpbkPSyUaqHH+CBWRhrsWw97VsHu1bnv+/8ECv6ODdl7zN/dX+WhF7oxxhg/paRlsj85jee+XnVUYkxJz+TxactZuSOJ5LQMUtKy3GTrtnDTs/K0drM/Zx1jN2ARcpOnT0KMDA+lfo1wt6wakeEhbnk1dxundevbgo0I82nN+iyHV/NwoNDMdCcp714Fe1bnvu9bD1nuE9ASAnWPgwYdoNMAmP8eJO876lBSu1nAwrQEbowxHlBVktMy2X84jYTkdPYnp5GQnMb+w2kcOJzGgTxl6TnrirtUfDA1g4/+tyk3Ofq0PGMiw2hcK8JZV0wSzd/q9T1OeGgI7u3Pii0rE/b/5bakV+Um6r3rICvd3UigbmuI7QDtL3LeG7SHeidAWETuseqfANPvhfSU3LKwSDhnBIFiCdwYY0pJVTl4JIOEw07SPZCcm4Sd9+yy9DyfC+sUJQIxkWHUiQqnTnQ4TWMi6NykFnWiw52yqDBe/G4N+w+nHbVv05gIfh16TqC/csWSlQkHNuZtTe9eDXvXQqbPUCMxLZ0W9QnnOe+x7aF+WwiPKr6OuKud9++fhsStULuZk7yzywPAErgxxvjIylIOpmYUm4jzt44zCrkOHRoiTjKOdhJvi7pRdGkWk7OcnZTrRocRExVO3ahwakWGERpSdAs3Iiw0zz1wgMiwUB7u175Mfx4VSlYWJG7Oe39690qnRZ3h0zKu3dxJzm16Q4OOzufYdhAeXbr6464OaMLOzxK4MabSysxSklLyJtvspLw/OS2nxZx9eToh2Vlf2D3haiHiJNlop3V8XP0adG/pJOK60eE567ITcZ2ocGpGVCOkmGRcEtkd1SpCL/Qyp+q0cnevytuhbM9aSD+cu12tpk5ybn2W896gg5Ooq9f0LvYyJD5jqFRY8fHxOn/+fK/DMMYEUEZmFgkp6Ucl4gPJvp/ztpYTUtIp7FdceGgIddxE7Fyq9v2cr3UcFU5MdBg1q1erHPd+KwpVSNp+dK/vPWsg7WDudjUaOfels+9Px7qJOjLGs9DLkogsUNX4/OXWAjfG+KUsnzlOz8zKe084XyJ2WsXpbqvYaR0npWYUerzq1UKoG52biDs2qZUnEee0jqPCiXGXo8JDLRmXF6pwaNfRvb53r4YjibnbRcc6Lemu1/q0qNtDVF3vYveQJXBjTLGKeub4ghMb5STbwpJyTiJ2L1sfPFJ4Mo4KD83TIm5RN8pNwAUn4uwBN0wFcWjP0b2+d6+C1ITcbSLrOsn5xCud9wYdnFZ1dD3Pwi6P7BK6MeYo2b2qdySksiMxhfs/WUxCSvpR2wlQ1G+QGtWr5b007XNZOqd1HJV979hJyhFhlowrhcP7chN1dmt6z6q8z0pHxOS2on3fo2OdrvgGsEvoxhiXqpKUmsHOxFS2J6awMzGVHQkp7EhMZWdSKtsTnLLDacWPnqXAw/3aOa1hn0RcJ8rpyOXpYBwmOFIO5Ls/7Sbtw3tyt6ley0nOvs9Rx3aAmo0sUZeCJXBjKhFVJSklgx1JTkLekZDKzsQUtiemOok60SlPzpecQwQa1IygUe0I2jasSa+2DWhcO4LGMRE0rh3BneMXsivpyFH1NY2J5K6zCxwZ2VQ2qYlO57GcS99uwj60M3eb8BpO57ET+vl0KusAtZpYog4AS+DGVBDZyTmn1eyTkLPfdxaRnBvHRNCukZOcm8Q4ybpx7Qga144ktmZ1wooYF3rYBR0Keea4XcC+r/HIkYO5idq3Q1nSttxtwqKcAU7a9Mnb+7tWMwixqy7BYgncmHJAVUlMST8qGW9PSGVnUop7Lzr1qLGvQwQa1nKScftGNTm7XYOcpJydoBvUrH7skzbkU6WfOa6s0g47iTp/r+/EzbnbVItwEnWrM/Lep45paYm6HLBObMYEWHZyzk7G2xNS87Sgsz8XlpzzJ+TGtSNzLm3H1ih9cjYV2NJPix+6Mz3FGTI0e1Sy7ISdsJmcLoih4U6ijm2f99J3nVYQYp0KvWad2IwJAFUlITm90MvZ2cup6XnHvA4NERrWrE6j2hF0aFKLPu0b0Kh2BE1ichO1JWdTpKWf5p08I3ELfHE3bP4dImvndig7sBHU/fcXEgb1joem3aDrwNxHtOq0hlBLBxWNnTFjCqGqHEhOz2klOx3Bci9n70wqOjk3jomkY5NanNO+AY1jIt3Wc+495+LGujYmR1aW8/jVoZ1wcJfz/u1jeWe+Asg4AvPfBQmFem2g0Ylw4tW5rep6bSA0zJvvYMqcJwlcRF4CLgHSgA3ATaqa4LO+BbASeFJVX/YiRlO5ZSfn7EemdiQ5j1LlebQqMfWoqRtDQ4RG7j3nTk1qcW6HBjSqHUmT2tmdwiw5m2OQme6MQJadlH0/+74f3p07D3WxBIbvgGrVAxq68Z5XLfCZwDBVzRCRF4BhwKM+618BvvEkMlNulHToTlVl/+E09/K122o+qtd2Kmn5knO1EMm559y5aW3O69SIRtn3oN0WdP0alpyNH9KSj07CBb37DmqSQyC6PtRo6LwadHTeazbK+z72YkjaevTutZtZ8q4iPEngqjrDZ3EucGX2gohcBvwFHMZUWYUN3amqnNU29qiE7PTYTnEvaxednOOaxdCvUwSNakW4j1M5Leh6lpxNUVSd4T7zJGH3dXBn3vcjSUfvH1ItNynXaQnNex6dlGs2ckYh8+cy97lP5L0HDhAW6XRkM1WC573QRWQ6MFFVx4lIDZzWeV/gIeBQYZfQRWQwMBigRYsW3Tdt2hSskE0QnD5yNtsSUorfkNzk7JuM8/TYdlvOgZjS0VQCWVmQvDdfEi6gtXxoN2SkHr1/WJRPEm7gzIxVs2G+90bO+N5l/eiVP73QTYUX9F7oIjILaFTAquGqOs3dZjiQAYx31z0JvKqqh4qbJUhVxwBjwHmMrIzCNuXE9iKS9xOXdMxJzI1jIqgfbcnZFCAjLV8L2U3C+VvLh3aDFjBsbETt3CTc/JS8ydi31Vy9pnejjMVdbQm7CgtYAlfVc4taLyKDgIuBczT3MsDJwJUi8iIQA2SJSKqqvhGoOE351CQmssAWeNOYSG46vbUHEZlyI+3w0Um4oPeU/QXsLM4l6uxk3LBzAa1l9zJ3WGTQv5oxx8KrXujnA48AvVQ1ObtcVc/02eZJnEvolryroDNOqMfEeXk76NjQnZWYqjMpRp6WcSE9stMOHr1/SJjbKm7oPNPc4pSCL2VHx9rzzqbS8Opf8htAdWCme6l8rqoO8SgWU84s35bIlEXbadewRs6UljZ0ZzlQkvutWZlweG8hPbF9H5vaBZlHT5ZCWHRu8m10Ihzf17nPnOcydiOIrGNDe5oqx6te6MVOX6SqTwYhFFPOJCSnMWTcAupHh/Px4FOpGx3udUgGChn16x7Y9yc06uQm490FP7+sWUcfLyImNwm3PLWAx6TcVnP1mkH9msZUJHYtyZQbWVnKg58uYVdSKp/ebsm7XJn1VAGjfqXCj8/nLkuIc4k6Owk3istNyr4JukZDCIsIbvzGVEJ+J3ARiQZSVQvqrmlM6b01Zz2zV+/mmf6dOKlFHa/DMUcOwboZsHJawQOGACAweI6TnKPq2/1lY4Ko0P9tIhIC/A0YCPQAjgDVRWQv8BXwjqquD0qUptL7ed0e/j1zLZd1bcL1p7T0OpyqKzUR1n7nJO31s5xWdnQshEc7vb/zq90MmnQNepjGmKJb4D8As3CGOV2u6tzIEpG6wNnACyIyRVXHBT5MU5ltT0jhvk8Wc0KDGvzr8hMpbgwAU8aS98Oab5yk/ecPkJkGNRtDtxuhY3+nR/fyz23UL2PKmaIS+Lmqmp6/UFX3A58Dn4uITWtjSiUtI4s7xy8kLSOL0dd3JyrcLsEGxeG9sPpLJ2n/9ZMzUUbt5tBzMHS4FJr1yNurO7u3uY36ZUy5UdRvy5pFtYRUdX9BCd6YY/HsVytZvCWB0QO70Sa2htfhVG4Hd8Kq6U7S3vSr0zu8Tms49W6npd3kpKJHFLNRv4wpV4pK4AsABQr6H63AcQGJyFQZUxdt48P/beK2M1tzwYmNvQ6nckrcCiu/gFVfwOa5gEL9tnDmP5yk3bCzd8OAGmNKpdAErqo2XqUJmDU7DzJs8jJ6tqrLI+e39zqcymX/X07CXvkFbJvvlDXsDL2HOUm7gf28jakM/LrhKCKXAme5i3NU9cvAhWQqu4Op6dwxbgE1IqrxxnUnERZqI2iV2t71sHKqk7h3LHHKGnd17lN36A/1ix07yRhTwRSbwEVkJM5jZNkzht0nIqep6mMBjcxUSqrKI5OWsml/MhNuPZkGtWxAjxJRhT2rnfvZK6fB7pVOebMe0PcZ6Hgp1GnlaYjGmMDypwV+IdDV5zGy/wMWAZbAzTF775e/+Gb5ToZf2IGTj6vndTgViyrsXJabtPetAwRanArnvwAdLnZ6hxtjqgR/n9mJAbLn5qsdmFBMZffHX/t5/pvVnN+pEbeeaV0s/KIK2xfmJu0DG50hS1udASffDh0ucUZBM8ZUOf4k8OeBRSLyA06P9LOAoQGNylQ6uw+mcteEhbSoG8VLV8XZYC1FycqCrfOchL3qC2fikJBq0LoXnPEgtL8Iout7HaUxxmPFJnBV/VhE5uDcBwd4VFV3BjQqU6lkZGZx94RFHErNYNwtJ1Mzwsb/OUpWJmz6zUnYq6bDwR0QGg5t+sDZj0G7C5wpM40xxuXvJfRYn+1PExFUdXKAYjKVzIvfreGPv/bz2jVdadfIpofMkZkOG392Hvda/SUc3gPVIuD4c6HjZdC2H0TU8jpKY0w55U8v9PeBOGAFkD2xrwKWwE2xvl2+gzE//cnfT2nJZSc19Toc72WkwZ9zYNU0WP0VpByAsGhoe57zjPbxfaG6jUhnjCmePy3wU1S1Y8AjMZXOn3sO8dBnS+nSPIZ/XtzB63C8k54CG2Y7Le0138CRRKheC9qe7ybtc5yJQYwx5hj4k8D/JyIdVXVlWVUqIi8BlwBpwAbgJlVNcNfFAe8AtXBa/D1UNbWs6jbBkZyWwR3jFhIWKowe2I3q1UK9Dim40g7DuplOR7R1MyDtEETEOL3GO14Kx/WGatW9jtIYU4H5k8A/xEniO3HmBBdAVTWuFPXOBIapaoaIvIAzZemjIlINGAf8XVWXiEg9wCZMqWBUleFTlrN290E+vLknTWKqSOsyNclJ1iunwrpZkJECUfXhxCudGb5anwWh1oHPGFM2/Eng7wF/B5aRew+8VFR1hs/iXOBK9/N5wFJVXeJut68s6jPBNe73zUxZtI1/9G3LmSfEFr9DRZZyANZ867S0N8yGzCNQoyGcdL3T0m5xGoTaFKnGmLLnz2+WPar6RQBjuBmY6H5uC6iIfIfT8/0TVX2xoJ1EZDAwGKBFixYBDM8ci8VbEnh6+grObhfLXWdX0vG3D++DNV85SfvPHyErHWo1hR63OPe0m/XMO5e2McYEgD8JfJGITACm41xCByj2MTIRmQUUNETUcFWd5m4zHMggd5z1asAZOM+cJwPfi8gCVf0+/0FUdQwwBiA+Pl79+B4mwPYfTuPOcQtoWCuCV6/pSkhIJRqs5dDu3Lm0N/4CmgkxLeGUO9y5tLtZ0jbGBJU/CTwSJ3Gf51NW7GNkqnpuUetFZBBwMXCOqmYn4K3AT6q6193ma6AbcFQCN+VLZpZy3yeL2Hs4jc+HnEZMVLjXIZVe0vbcpL3pN0Ch3vFwxv1O0m4UZ3NpG2M8U2gCF5FrgRmqelNZVyoi5wOPAL1UNdln1XfAIyIShdNDvRfwalnXb8ref2at5ed1exl5+Ymc2KwCD5efsNl53GvlNNj6h1MW2wF6PerOpd3BkrYxplwoqgXeAvhMRMJwWsDfAH/4tJZL4w2gOjDTHRN7rqoOUdUDIvIKMA+nlf+1qn5VBvWZAJq9ehevz17PVd2bcU2P5l6Hc+z2bcgdd3z7IqesURz0+aczl3ZsW2/jM8aYAkhx+VhEagLnAucDPYFVwLfAd6q6K+AR+iE+Pl7nz5/vdRhV0pb9yVw86heaxkQy+c7TiAirIM9771njzvD1Bexa5pQ17e487tXxUqh7nLfxGWOMy+0LFp+/3J/JTA4CU9wXItIRuADn+fB+ZRynqUBS0zO5Y/wCVJW3r+9evpO3KuxakdvS3rPaKW9+CvT7lzPASow9zWCMqTiKugferYj9flDVfwcgHlOBPDV9Bcu3JfHuDfG0qBfldThHU4Udi3Nb2vs3OHNptzwd4m+BDhdDrSZeR2mMMSVSVAu8qAStQJ8yjsVUIJ/O38LHf2zhrrPbcG7Hhl6HkysrC7YtcEZDW/WF0ylNQp1R0E67G9pfDDUaeB2lMcaUWqEJXFXPDmYgpuJYsT2Rx6cu5/Tj6/Fg33Zeh+PMpb3l99yW9sHtEBIGbc52eo+3uxCi6nodpTHGlCm/xngUkc5ARyAiu0xVPwxUUKb8SkxJ545xC6kTFc5//nYSoV4N1pKZAZt+de9pT4fDuyG0ujuX9pPOXNqRMd7EZowxQeDPfOBPAL1xEvjXOB3YfsHpxGaqkKws5R+fLmZ7QgoTbz+V+jUCPJvW0k/h+6chcSvUbgZnPwbRDXLn0k7eB2FRcEJf5xntE86D6jUDG5MxxpQT/rTArwS6AItU9SYRaYgzY5ipYkb/uIFZq3bz5CUd6d6yTmArW/opTL/XmUsbIHELTL3D+Rxe02lhd+zvtLjDy2EHOmOMCTB/EniKqmaJSIaI1AJ2AxVwtA5TGr+u38u/Z6zhki5NuPG0VoGv8Punc5O3r6j68MAKCIs4ep0xxlQh/iTw+SISA/wXWAAcAv4XyKBM+bIjMYV7P17EcbE1GHn5iUgwhhJN3FpwefI+S97GGIN/A7nc6X58W0S+BWqp6tLAhmXKi7SMLO4av5DU9Ezevr470dWDNLd1rSaQtO3o8trNglO/McaUc8XOfygiA0SkNoCqbgQ2i8hlAY7LlBP/+noVCzcn8OKVXTi+QY3gVKoK0bFHl4dFwjkjghODMcaUc/5MYPyEqiZmL6hqAvBEwCIy5cYXS7Yz9reN3HJGay6Kaxy8ihd84Iyg1vlKqN0cEOf9ktch7urgxWGMMeWYP9dDC0ryQbqOaryybtdBhn6+lPiWdRh6QfvgVbx7FXw7DNr0gcv/CyH+/I1pjDFVjz+/HeeLyCsi0sZ9vYLTmc1UUoeOZDBk3AKiwkN5c2A3wkKDlETTU2HSLc6z3Je9bcnbGGOK4M9vyHuANGAi8AlwBLgrkEEZ76gqj05ayl97DzPq2m40rBXEHt8zH4fdK+Cy0VCzHI2vbowx5ZA/l8IbqOpQ3wIR6QHMC0xIxkvv/7qRr5btYOgF7Tm1Tb3gVbzmG/hjDJxylzOymjHGmCL50wL/XESaZi+IyFnA+4ELyXhl3sb9PP/1Ks7r2JDbzzoueBUnbYepd0KjODjX+kcaY4w//EngtwNTRaSRiFwIjAIuLE2lIvKSiKwWkaUiMsUdKAYRCROR/xORZSKySkSGlaYe4789B49w1/iFNKsTyctXdwnOYC3gzCQ25XbISIUr34dqAR5f3RhjKoliE7iqzgPuBWYATwLnquqWUtY7E+isqnHAWiA7UV8FVFfVE4HuwO0i0qqUdZliZGRmcc/HC0lKTWf09d2pFREWvMp//Q/89RNc8CLUPyF49RpjTAVX6D1wEZkOqE9RFJAIvCciqOqlJa1UVWf4LM7FmTAFt75oEakGROJ0nksqaT3GPy/PWMvcP/fzytVd6NC4VvAq3roAfngOOg2Ak64PXr3GGFMJFNWJ7eUgxXAzTg93gElAf2AHzh8MD6jq/oJ2EpHBwGCAFi1aBCHMyum7FTt5+8cNDDy5BZd3C+IwpalJ8PnNULMJXPwaBOuSvTHGVBJFJfCfVFWLWI+ISGHbiMgsoFEBq4ar6jR3m+FABjDeXdcTyASaAHWAn0Vklqr+mf8gqjoGGAMQHx9fZJymYBv3HuahT5cQ16w2Iy7pGNzKv/oHJGyGm76ByJjg1m2MMZVAUQn8BxH5HJimqpuzC0UkHDgDuBH4ARhb0M6qem5RFYvIIOBi4ByfPwKuA75V1XRgt4j8CsQDRyVwUzopaZkMGbeA0FDhrYHdqF4tNHiVL/kEln0KvR+DFqcEr15jjKlEiurEdj5Oa/hjEdkuIitF5C9gHXAt8Jqqji1JpSJyPvAIcKmqJvus2gz0cbeJBk4BVpekDlM4VWX41GWs2XWQ167pSrM6UcGrfN8Gp/Xd8nQ466Hg1WuMMZVMoS1wVU0F3gLeEpEwoD6Q4k5mUlpvANWBme7jSnNVdQjwJvCBiKwABPjApi4texP+2Mzkhdu4/9wT6N2uQfAqzkiDz2+FkGpw+RgICWKr3xhjKhm/JiVxL2nvKKtKVfX4QsoP4TxKZgJkyZYEnvpiJb3axnJvnyA/tvXDs7B9IVz9kc3rbYwxpWSzRVQhBw6ncef4hcTWrM5r13QlJCSIPb83/OA88919EHQs8ROIxhhjXDYtaBWRmaXcN3Exew4e4bMhp1InOjx4lR/e64y2Vr8d9Hs+ePUaY0wlVmwLXETuEZE6wQjGBM6o2ev4ae0enry0E12axwSvYlVnnPOUBGeo1PAgdpgzxphKzJ9L6A2BeSLyqYicL0EbJNuUlTlrdvOf79dxRbdmXNuzeXAr/2MMrPsOznsGGnUObt3GGFOJ+TMW+j+BE4D3gEHAOhH5l4i0CXBspgxsPZDM/RMX065hTZ69rHPwJikB2LkMZvwT2p4PPQcHr15jjKkC/OrE5g60stN9ZeCMkjZJRF4MYGymlI5kZHLn+IVkZipvX9+dyPAgPraVlgyTbobIutD/TRsq1RhjylixndhE5D7gBmAv8C7wsKqmi0gIzqAujwQ2RFNST01fydKtiYz5e3da1Y8ObuXfDYO96+CGqRBdP7h1G2NMFeBPL/S6wOWqusm3UFWzROTiwIRlSmvSgq1M+H0zQ3q14bxOBQ1JH0Arp8GCsXD6/XBc7+DWbYwxVYQ/Cfw/ACJS16fsoKqmq+qqwIRlSmPl9iSGT1nGKcfV5aHz2ga38oQt8MU90KQb9PlncOs2xpgqxJ974AuBPcBanEvme4CNIrJQRLoHMjhz7BJT0rlj/AJqR4Yx6tpuVAsN4lg9WZkwebDzfuV7EBoWvLqNMaaK8ee3+0zgQlWtr6r1gAuAL4E7ccZKN+WEqvLQZ0vYdiCFtwZ2I7Zm9eAG8NPLsPk3uOjfUPe44NZtjDFVjD8J/BRV/S57QVVnAKeq6lycCUlMOfHOT38yc+UuHruwA/Gt6ha/Q1naPBd+HAlx10CXvwW3bmOMqYL8uQe+Q0QeBT5xl68BdolIKJAVsMjMMfnfhn28+O1qLoprzE2ntwpu5SkJzixjMS3gwpeDW7cxxlRR/rTArwOaAVOBKUBztywUuDpgkRm/7UpK5Z6PF9K6fjQvXBEX3MFaVGH6fXBwB1zxPkTUCl7dxhhThRXZAndb2f9R1YGFbLK+7EMyxyI9M4u7xi8kOS2Tj287hRrVgzw/zaKPYOVUOOcJaGZ9Go0xJliKbIGraibQUkSCOHWVORbPf72a+ZsOMPKKOE5oWDO4le9ZC988Cq17Oc98G2OMCRp/mmt/Ar+KyBfA4exCVX2lNBWLyDNAf5z76LuBQaq63Z0s5T/AhUCyW76wNHVVVl8u3c77v/7FoNNacWmXJsGtPOMIfH4zVIuAAe9AiE0tb4wxweTPb90NOI+NhQA1fV6l9ZKqxqlqV/f4I9zyC3AmTzkBGAyMLoO6Kp31uw/y6KSldGsRw2MXdgh+ALOedCYruewtqNU4+PUbY0wVV2wLXFWfAhCRKFVNLquKVTXJZzEaUPdzf+BDdwKVuSISIyKNVXVHWdVd0R0+ksGQcQuJCAvlzYHdCK8W5Nbvupkw9y1nhrF2FwS3bmOMMYAfLXAROVVEVgKr3eUuIlImA7iIyHMisgUYSG4LvCmwxWezrW5Z/n0Hi8h8EZm/Z8+esginQlBVhk5exp97DjHq2pNoXDsyuAEc3AVThkCDTtD3meDWbYwxJoc/TbfXgH7APgBVXQKc5c/BRWSWiCwv4NXfPdZwVW0OjAfuPpbAVXWMqsaranxsbOyx7Fqh/d9vG5m+ZDsP9WvHaccHeZavrCyYOgTSDsOV70NYRHDrN8YYk8OvZ45UdUu+Z4sz/dzvXD/jGA98DTwBbMN51jxbM7esyluw6QDPfrWKczs0ZMhZbYIfwNw3YcNsuPhVaNA++PUbY4zJ4U8LfIuInAaoiISJyENAqWchE5ETfBb7416iB74AbhDHKUCi3f+GvYeOcNf4hTSJieTfV3chJCSIg7UAbF8Es56CDpdA95uCW7cxxpij+NMCH4LzWFdTnJbwDOCuMqh7pIi0w3mMbJNbDzgt8QtxBolJBqp8tsjIzOLejxdxIDmNyXeeRu3IIM/ydeQQTLoFajSAS16HYI70ZowxpkD+9ELfi9PJrEyp6hWFlCtl8wdCpfHKzLX8tmEfL10ZR6cmtYMfwDePwIG/4MbpEBXkSVKMMcYUqNgELiKxwG1AK9/tVfXmwIVlss1cuYu35mzg2p7NuSq+efE7lLVlk2DxeDjrEWh1RvDrN8YYUyB/LqFPA34GZuFn5zVTNjbtO8yDny6mc9NaPHFJp+AHcGAjfPkAND8Zej0a/PqNMcYUyp8EHqWq9ts7yFLTMxkybiEhIowe2J2IsNDgBpCZ7kwRisDl/4XQIE+SYowxpkj+9EL/UkQuDHgkJoeq8s+py1m9M4nXrulK87pRwQ9izkjYOg8ueRXqtAx+/cYYY4rkTwK/DyeJp4pIkogcFJGkYvcyJTZx3hYmLdjKPX1O4Oz2DYIfwF8/w8//hpOuh84F9jU0xhjjMX96oQd5jsqqbdnWREZ8sYIzT6jPfeecUPwOZS15P0weDPXawAUvBr9+Y4wxfvFnLHQRketF5HF3ubmI9Ax8aFVPQnIad4xfQP3ocP7zt5MIDfZgLaow7W44vAeueA/Co4NbvzHGGL/5cwn9LeBU4Dp3+RDwZsAiqqKyspQHJi5mV1Iqb13fnbrR4cEPYv57sOYrOPdJaNI1+PUbY4zxmz9di09W1W4isghAVQ+IiAfZpXJ744f1/LBmD89c1pmuzWOCH8CulfDdcDj+XDjlzuDXb4wx5pj40wJPF5FQ3Pm63YFdsgIaVRXz09o9vDprLQNOasr1J7cIfgDpKfD5LVC9Flw2GkKCPL+4McaYY+bPb+rXgSlAAxF5DvgF+FdAo6pCtiWkcN8ni2jboCbPDeiMeDHO+IzHYfdKGDDaGe/cGGNMuedPL/TxIrIAOAcQ4DJVLfVsZAaOZGRy5/iFpGcqo6/vRlS4B4OlrP4K5v0XTr3buXxujDGmQvB3PvDV5E73acrIs1+uYsmWBN6+vhvHxdYIfgBJ22HaXdAoDs4ZEfz6jTHGlJjd7PTIlEVb+WjuJgafdRznd24c/ACyMp3nvTPS4MoPoFr14MdgjDGmxGyAaw+s3pnEsMnL6Nm6Lo/0a+dNEL++Bht/hv5vQv3jvYnBGGNMiVkLPMiSUtO5Y9xCakaE8cZ1J1Et1INTsGUezH4OOl0OXct8qndjjDFB4M9IbJeLyDoRSbSx0EtHVXn4syVs3p/Mm9d1o0HNiOAHkZroPDJWqylc/Cp40evdGGNMqfnT/HsRuFRVa6tqLVWtqaq1SlOpiDwjIktFZLGIzBCRJm75QLd8mYj8JiJdSlNPefPfn//kuxW7GHZBe3q2rhv8AFThq39A4la44l2IjAl+DMYYY8qEPwl8VwAeG3tJVeNUtSvwJZDdBfovoJeqngg8A4wp43o9M/fPfbzw7Rou6NyIW85o7U0QSz6BZZ9B72HQ4mRvYjDGGFMm/OnENl9EJgJTgSPZhao6uaSVqqrvJfho3FHeVPU3n/K5QLOS1lGe7E5K5e4Ji2hZN4oXr4zzZrCWfRvg64eg5Rlw5oPBr98YY0yZ8ieB1wKSgfN8yhQocQIHcEd1uwFIBM4uYJNbgG+K2H8wMBigRQsPhh/1U3pmFndNWMjhIxmMv/VkakaEBT+IjDSYdDOEVIPL34GQ0ODHYIwxpkz5MxLbTSU5sIjMAhoVsGq4qk5T1eHAcBEZBtwNPOGz79k4CfyMIuIag3uJPT4+XksSYzC8+O1q5m08wGvXdKVdI4+mVp/9DOxYDNeMg9qV4qKGMcZUecUmcBFpBowCTneLfgbuU9WtRe2nqv6Oyzke+Bo3gYtIHPAucIGq7vPzGOXS18t28N+f/+KGU1ty2UlNvQliw2z47XWIvxk6XOJNDMYYY8qcP53YPgC+AJq4r+luWYmJyAk+i/1xh2kVkRY4l+b/rqprS1OH1zbsOcQjk5bStXkM/7yoozdBHNoDU4ZAbHs47zlvYjDGGBMQ/twDj1VV34Q9VkTuL2W9I0WkHc60pJuAIW75CKAe8Jbb0StDVeNLWVfQJadlcMe4BYRXC+Gtgd0Ir+bBYC2qMO1OSEmAv0+B8Kjgx2CMMSZg/Eng+0TkeuBjd/laoFSXtlX1ikLKbwVuLc2xvaaqDP18Get2H+Kjm0+mSUykN4H8/jasmwEXvAQNO3kTgzHGmIDxp2l4M3A1sBPYAVwJlKhjW1Xw0dxNfLFkO//o25YzTqjvTRA7lsLMEdD2Auh5mzcxGGOMCSh/eqFvAi4NQiwV3sLNB3jmy5X0ad+AO3t7NEFI2mFnqNTIus5EJTZUqjHGVEqFJnAReURVXxSRUbgDrfhS1XsDGlkFs+/QEe4av5CGtSJ49equhIR4lDi/HQZ718EN0yC6njcxGGOMCbiiWuDZw6fOD0YgFVlmlnLfJ4vZdziNyXecRu0oDwZrAVgxFRb+H5zxABzXy5sYjDHGBEWhCVxVp7sfk1X1M991InJVQKOqYF6btZZf1u/lhStOpHPT2t4EkbAFpt8LTbvD2cO9icEYY0zQ+NOJbZifZVXS7NW7GDV7PVfHN+OaHh4N6ZqZAZNvg6wsZ5axUI+uABhjjAmaou6BXwBcCDQVkdd9VtUCMgIdWEWwZX8y93+ymE5NavF0/87eBfLzy7D5f3D5f6Hucd7FYYwxJmiKuge+Hef+96XAAp/yg8ADgQyqIkhNz2TIOOfHMnpgdyLCPJogZNNv8OMLEPc3iLvamxiMMcYEXVH3wJcAS0RkAiBAe5ze6GtUNS1I8ZVbT0xbwYrtSbx3Yzwt6nk0ylnKAfj8NohpCRe97E0MxhhjPOHPSGx9gXeADTiJvLWI3K6qhU71WdlNnLeZifO3cPfZx3NOh4beBKEK0++DQzvhlhlQ3aOZzowxxnjCnwT+CnC2qq4HEJE2wFcUMVd3ZbZ8WyKPT1vB6cfX44G+bb0LZOGHsHIanPuU0/PcGGNMleJPL/SD2cnb9SfOffAqJzE5nTvGL6BedDiv/+0kQr0arGXPGvh2KBzXG06z8XSMMaYq8qcFPl9EvgY+xbkHfhUwT0QuB1DVyQGMr9zIylIe/HQxOxNTmXj7qdSrUd2bQNJTYdItEBYJA96BEA9mOjPGGOM5fxJ4BLALyB7aaw8QCVyCk9CrRAIf/eMGvl+9m6cu7US3FnW8C2TWk7BrGVw7EWo28i4OY4wxnioygYtIKLBUVV8NUjzl0i/r9vLvGWu4tEsTbji1pXeBrJ0Bv4+Gk4dAu/O9i8MYY4znirz+qqqZOPN/V1k7ElO495NFHN+gBiOvOBHxanavgzth6h3Q8ESn45oxxpgqzZ9L6L+KyBvAROBwdqGqLgxYVOVEWkYWd45fyJH0TEZf352ocH9+XAGQlQVTbnemCr3yPQiL8CYOY4wx5YY/Gamr+/60T5kCfUpTsYg8A/QHsoDdwCBV3e6zvgfwP+BvqjqpNHWV1HNfrWTR5gTeGtiNNrE1vAjB8b9R8OccuPg1iG3nXRzGGGPKjWITuKqeHaC6X1LVxwFE5F5gBDDEXQ4FXgBmBKjuAk1dtI2XvlvD9oQUYqLCOJCczq1ntObCExsHM4y8ti2E75+GDpdA90HexWGMMaZcKfYZJBGpLSKviMh89/VvESn1nJmqmuSzGI3Tqs92D/A5Tss8KKYu2sawycvYlpCCAgeS0wkR6NDYwxHOjhyEz2+BGo3gktfBq/vvxhhjyh1/HiJ+H2fglqvdVxLwQVlULiLPicgWYCBOCxwRaQoMAEYXs+/g7D8q9uzZU+pYXvpuDSnpmXnKshRembmu1Mcusa8fgQMb4Yr/QlRd7+IwxhhT7viTwNuo6hOq+qf7egrwa85KEZklIssLePUHUNXhqtocGA/c7e72GvCoqmYVdWxVHaOq8aoaHxsb6084RdqekHJM5QG39DNYMgHOehhanuZNDMYYY8otfzqxpYjIGar6C4CInA74ldVU9Vw/4xgPfA08AcQDn7iPa9UHLhSRDFWd6uexSqRJTCTbCkjWTWIiA1ltwfb/BV8+AM1PgbMeCX79xhhjyj1/WuB3AG+KyEYR2QS8gdvZrDRE5ASfxf7AagBVba2qrVS1FTAJuDPQyRvg4X7tiMw3p3dkWCgP9wtyr+/MdPj8VpAQ59J5qEePrhljjCnX/OmFvhjoIiK13OWkovfw20gRaYfzGNkmyuCPgtK47KSmADm90JvERPJwv3Y55UEz53nYNh+uGgsxLYJbtzHGmApDVLXoDUTuw+m0dhD4L9ANGKqqQX3Eqyjx8fE6f/58r8MovT9/hA/7w0nXQ/83vI7GGGNMOSAiC1Q1Pn+5P5fQb3Zb3ecB9YC/AyPLOD5zeJ8z2lq94+GCF7yOxhhjTDnnzw3W7IePLwQ+VNUV4tmA4JWUKnxxNyTvg+smQni01xEZY4wp5/xpgS8QkRk4Cfw7EamJc9/alJV578Kar51JShp38ToaY4wxFYA/LfBbcMZD/1NVk0WkHnBTQKOqSnatgO+Gw/F94ZQ7vI7GGGNMBVFoAheRbvmKjrMr52UsPQUm3QIRteGy0TZUqjHGGL8V1QL/t/seAXQHluLcD48D5gOnBja0KuC74bBnFVw/GWqUfjQ5Y4wxVUeh98BV9Wx3JrIdQHd32NLuwEnAtmAFWGmt+hLmvwen3QPHn+N1NMYYYyoYfzqxtVPVZdkLqroc6BC4kKqAxG1Or/PGXaHPCK+jMcYYUwH504ltqYi8C4xzlwfiXE43JZGVCZMHQ0YaXPk+VAv3OiJjjDEVkD8J/Cac8dDvc5d/pJipPk0RfnkFNv0C/d+Cem28jsYYY0wFVewldFVNVdVXVXWAqg4AFgCvBD60SmjLH/DD89D5Cuh6ndfRGGOMqcD8mupKRE4CrgWuBv4CJgcyqEopNRE+vwVqN4WLX7VHxowxxpRKUc+Bt8VJ2tcCe4GJOJOfnB2k2CoPVWd+78RtcPN3znPfxhhjTCkU1QJfDfwMXKyq6wFE5IGgRFXZLJ4Ayz+HPv+E5j28jsYYY0wlUNQ98MtxngH/QUT+KyLnkDuxifHX3vXw9cPQ6kw440GvozHGGFNJFNoCV9WpwFQRiQb6A/cDDURkNDClPM0HXm5lpDn3vauFw4B3ICTU64iMMSWUnp7O1q1bSU1N9ToUU0lFRETQrFkzwsLC/Nq+2E5sqnoYmABMEJE6wFXAo0CJE7iIPIPzR0EWsBsYpKrb3XW9gdeAMGCvqvYqaT2em/007FgMf5vgdF4zxlRYW7dupWbNmrRq1QqbF8KUNVVl3759bN26ldatW/u1jz8jsflWcEBVx6hqacf+fElV41S1K/AlMAJARGKAt4BLVbUTzh8LFdP67+G3URB/C7S/yOtojDGllJqaSr169Sx5m4AQEerVq3dMV3iOKYGXFVVN8lmMBtT9fB0wWVU3u9vtDnZsZeLQHpgyBGI7QL/nvI7GGFNGLHmbQDrWf19+PQceCCLyHHADkAhkP5rWFggTkTlATeA/qvphIfsPBgYDtGjRIuDx+i0rC6be4Tz3fcNUCIv0OiJjjDGVUMBa4CIyS0SWF/DqD6Cqw1W1OTAeuNvdrRrO1KUXAf2Ax93n0Y/iXsqPV9X42NhyNBXn72/D+plOy7thJ6+jMcZ4ZOqibZw+cjath37F6SNnM3VR2UziuHPnTv72t7/Rpk0bunfvzoUXXsjatWvL5NgAc+bM4bfffiuz4xXltNNOC0o9gbJx40YmTJhwzPsNGjSISZMmlbr+gCVwVT1XVTsX8JqWb9PxwBXu563Ad6p6WFX3Aj8BXQIVY5nbsQRmPQHtLoIet3odjTHGI1MXbWPY5GVsS0hBgW0JKQybvKzUSVxVGTBgAL1792bDhg0sWLCA559/nl27dpVN4JQsgWdkZJSormD9oVCY/HEf6/coaQIvK55cQheRE1R1nbvYH2fQGIBpwBsiUg0IB04GXvUgxGOXdhgm3QJR9eDSUTZUqjGV2FPTV7Bye1Kh6xdtTiAtMytPWUp6Jo9MWsrHf2wucJ+OTWrxxCVFX7X74YcfCAsLY8iQITllXbp0QVV5+OGH+eabbxAR/vnPf3LNNdcwZ84cXn75Zb788ksA7r77buLj4xk0aBCtWrXixhtvZPr06aSnp/PZZ58RERHB22+/TWhoKOPGjWPUqFG0b9+eIUOGsHmzE/drr73G6aefzpNPPsmGDRv4888/adGiBf369WPKlCkkJiaybds2rr/+ep544gkAXnnlFd5//30Abr31Vu6//34AatSowaFDh9ixYwfXXHMNSUlJZGRkMHr0aM4880xmzJjBE088wZEjR2jTpg0ffPABNWrUOOrn8sILLzBu3DhCQkK44IILGDlyJIsXL2bIkCEkJyfTpk0b3n//ferUqUPv3r3p2rUrv/zyC9deey3Tp0/Ps9y7d28efPBBDh06RP369Rk7diyNGzdm/fr1DBkyhD179hAaGspnn33G0KFDWbVqFV27duXGG2/k3nvvZejQocyZM4cjR45w1113cfvtt6Oq3HPPPcycOZPmzZsTHl42s1B6dQ98pIi0w3mMbBMwBEBVV4nItzjTlWYB77rzj5d/3zwK+9bDjV9AdD2vozHGeCh/8i6u3F/Lly+ne/fuR5VPnjyZxYsXs2TJEvbu3UuPHj0466yzij1e/fr1WbhwIW+99RYvv/wy7777LkOGDKFGjRo89NBDAFx33XU88MADnHHGGWzevJl+/fqxatUqAFauXMkvv/xCZGQkY8eO5Y8//mD58uVERUXRo0cPLrroIkSEDz74gN9//x1V5eSTT6ZXr16cdNJJOXFMmDCBfv36MXz4cDIzM0lOTmbv3r08++yzzJo1i+joaF544QVeeeUVRowYkec7fPPNN0ybNo3ff/+dqKgo9u/fD8ANN9zAqFGj6NWrFyNGjOCpp57itddeAyAtLY358+cDMH369Jzl9PR0evXqxbRp04iNjWXixIkMHz6c999/n4EDBzJ06FAGDBhAamoqWVlZjBw5Ms8fSGPGjKF27drMmzePI0eOcPrpp3PeeeexaNEi1qxZw8qVK9m1axcdO3bk5ptvPsazfzRPEriqXlHEupeAl4IYTumtmAKLPnJGWmtd/H8aY0zFVlxL+fSRs9mWkHJUedOYSCbefmqZx5PdegwNDaVhw4b06tWLefPmUatWrSL3u/zyywHo3r07kycXPEfVrFmzWLlyZc5yUlIShw4dAuDSSy8lMjK3o27fvn2pV69ezrF/+eUXRIQBAwYQHR2dU/7zzz/nSeA9evTg5ptvJj09ncsuu4yuXbvy448/snLlSk4//XTASbqnnnr0z27WrFncdNNNREVFAVC3bl0SExNJSEigVy9nGJEbb7yRq67KfSr5mmuuyXOM7OU1a9awfPly+vbtC0BmZiaNGzfm4MGDbNu2jQEDBgDOgCsFmTFjBkuXLs25v52YmMi6dev46aefcs5PkyZN6NOnT4H7HyvPeqFXGgmb4Yv7oGk8nP2Y19EYY8qBh/u1Y9jkZaSkZ+aURYaF8nC/dqU6bqdOnY6p81O1atXIyspt9ed/xrh69eoAhIaGFnr/Nysri7lz5xaYtLKTcrb8j0H5+1jUWWedxU8//cRXX33FoEGDePDBB6lTpw59+/bl448/zrPt77//zu233w7A008/7dfxi4s7e1lV6dSpE//73//yrD948KBfx1VVRo0aRb9+/fKUf/311yWKsziePAdeaWRmwOe3gWbBFe9CqH/D3xljKrfLTmrK85efSNOYSASn5f385Sdy2UmlG5GxT58+HDlyhDFjxuSULV26lJiYGCZOnEhmZiZ79uzhp59+omfPnrRs2ZKVK1dy5MgREhIS+P7774uto2bNmnkS1nnnnceoUaNylhcvXlzovjNnzmT//v2kpKQwdepUTj/9dM4880ymTp1KcnIyhw8fZsqUKZx55pl59tu0aRMNGzbktttu49Zbb2XhwoWccsop/Prrr6xfvx6Aw4cPs3btWk4++WQWL17M4sWLufTSS+nbty8ffPABycnJAOzfv5/atWtTp04dfv75ZwA++uijnNZ4Udq1a8eePXtyEnh6ejorVqygZs2aNGvWjKlTpwJw5MgRkpOTj/pZ9evXj9GjR5Oeng7A2rVrOXz4MGeddVbO+dmxYwc//PBDsbH4w1rgpfHTi7BlLlz+LtT1b+g7Y0zVcNlJTUudsPMTEaZMmcL999/PCy+8QEREBK1ateK1117j0KFDdOnSBRHhxRdfpFGjRgBcffXVdO7cmdatW+e5bF2YSy65hCuvvJJp06YxatQoXn/9de666y7i4uLIyMjgrLPO4u233y5w3549e3LFFVewdetWrr/+euLj4wHnsamePXsCTie2/HHMmTOHl156ibCwMGrUqMGHH35IbGwsY8eO5dprr+XIkSMAPPvss7Rtm/fJ4vPPP5/FixcTHx9PeHg4F154If/617/4v//7v5xObMcddxwffPBBsd89PDycSZMmce+995KYmEhGRgb3338/nTp14qOPPuL2229nxIgRhIWF8dlnnxEXF0doaChdunRh0KBB3HfffWzcuJFu3bqhqsTGxjJ16lQGDBjA7Nmz6dixIy1atCjwVkBJiKoWv1U5Fx8fr9kdEoJm028w9iKIuwYGFPyP2RhTeaxatYoOHTp4HUa5NXbsWObPn88bb7zhdSgVWkH/zkRkgarG59/WLqGXRMoB59J5nVZwYcXqb2eMMaZysEvox0oVvrgXDu2CW2ZA9ZpeR2SMMZ4bNGgQgwYN8jqMKsVa4MdqwVhY9QWc8zg07eZ1NMYYY6ooS+DHYvdq+HYYHHc2nHqP19EYY4ypwiyB+ys9FT6/BcKjnE5rIfajM8YY4x27B+6vWU/AruVw3WdQs5HX0RhjjKnirBnpjzXfOtOEnnwHtD3P62iMMRXB0k/h1c7wZIzzvvTTUh9y48aNdO7cOU/Zk08+ycsvv1zqYwfC4sWLSzQKWe/evQn6o8EVkCXw4hzcCdPuhIYnQt+nvI7GGFMRLP0Upt8LiVsAdd6n31smSTyYMjMzi1wuTkkTuPGPXUIvSlYWTLkd0pLhyvegWnWvIzLGlAffDIWdywpfv3UeZB7JW5aeAtPuhgX/V/A+jU6EC0aWOKTevXtz8skn88MPP5CQkMB7773HmWeeSWZmJo8++ijffvstISEh3Hbbbdxzzz18//33PPTQQ2RkZNCjRw9Gjx5N9erVadWqFddccw0zZ87kkUceYejQoXmW69atW+AUn/PmzeO+++7j8OHDVK9enZkzZzJixAhSUlL45ZdfGDZsGBdffDH33HMPy5cvJz09nSeffJL+/fuTkpLCTTfdxJIlS2jfvj0pKUdPBGOOZgm8KL+9Dn/OgUteh9jSTUJgjKlC8ifv4srLSEZGBn/88Qdff/01Tz31FLNmzWLMmDFs3LiRxYsXU61aNfbv309qaiqDBg3i+++/p23bttxwww2MHj06Z57uevXqsXDhQgCGDh2as7x3714uv/zyo6b4zE7yEydOpEePHiQlJREVFcXTTz+dZ3S2xx57jD59+vD++++TkJBAz549Offcc3nnnXeIiopi1apVLF26lG7d7BFdf1gCL8y2BTD7GejYH7rd4HU0xpjypLiW8qud3cvn+dRuDjd9VeJqC5vdK7vcd3rQjRs3As50m0OGDKFaNefXfd26dVmyZAmtW7fOGVf8xhtv5M0338xJ4IVNtzl37twCp/hcs2YNjRs3pkePHgCFTmM6Y8YMvvjii5x79qmpqWzevJmffvqJe++9F4C4uDji4uKO7QdTRVkC97X0U/j+aUjc6jwmVr02XPIf8HNKPGOMAeCcEc4973SfS8FhkU55KdSrV48DBw7kKdu/fz+tWzuTKfkzPag/ippus6ApPpctK+J2gg9V5fPPP6ddO7uiWRasE1u2/J1OsjIh7TCsm+l1ZMaYiibuaufWW+3mgDjvl7zulJdCjRo1aNy4MbNnzwac5P3tt99yxhlnFLpP3759eeedd3IS+v79+2nXrh0bN27MmarT3+k2C5vis127duzYsYN58+YBzvzZGRkZBU63OWrUKLIn0Vq0aBHgzAc+YcIEAJYvX87SpUuP6edSVXmWwEXkGRFZKiKLRWSGiDRxy2uLyHQRWSIiK0TkpqAE9P3Tef9aBud+1fclmzDeGFPFxV0NDyyHJxOc91Im72wffvghzzzzDF27dqVPnz488cQTtGnTptDtb731Vlq0aEFcXBxdunRhwoQJRERE8MEHH3DVVVdx4oknEhISwpAhQ4qt23eKz7i4OE499VRWr15NeHg4EydO5J577qFLly707duX1NRUzj77bFauXEnXrl2ZOHEijz/+OOnp6cTFxdGpUycef/xxAO644w4OHTpEhw4dGDFiBN27dy+Tn1Vl59l0oiJSS1WT3M/3Ah1VdYiIPAbUVtVHRSQWWAM0UtW0wo5VJtOJPhkDFPSzEOc/oDGmSrPpRE0wVIjpRLOTtyua3OypQE1xemXUAPYDJb+Z46/azY6t3BhjjPGQp/fAReQ5EdkCDASye3e8AXQAtgPLgPtUNauAfQeLyHwRmb9nz57SB3POCKeTia8y6HRijDHGBEJAE7iIzBKR5QW8+gOo6nBVbQ6MB+52d+sHLAaaAF2BN0TkqGcSVHWMqsaranxsbGzpgw1QpxNjTOXh1S1HUzUc67+vgD5Gpqrn+rnpeOBr4AngJmCkOt9kvYj8BbQH/ghMlD7irraEbYwpUEREBPv27aNevXqFPo9tTEmpKvv27SMiIsLvfTx7DlxETlDVde5if2C1+3kzcA7ws4g0BNoBf3oQojHG5GjWrBlbt26lTG7ZGVOAiIgImjXzv9+VlwO5jBSRdkAWsAnIfobhGWCsiCwDBHhUVfd6FKMxxgAQFhaWM2CKMeWBZwlcVa8opHw7YHN2GmOMMUWwkdiMMcaYCsgSuDHGGFMBeTYSW1kSkT0499HLSn3A7ruXH3Y+yhc7H+WPnZPypazPR0tVPep56UqRwMuaiMwvaNg64w07H+WLnY/yx85J+RKs82GX0I0xxpgKyBK4McYYUwFZAi/YGK8DMHnY+Shf7HyUP3ZOypegnA+7B26MMcZUQNYCN8YYYyogS+DGGGNMBeTlWOieE5FMnDnHw4AM4EPgVVXNEpHngIt8No8C2gC1VfVQ0IOtpEp6DoArgZeAbUC4u89/gxl7ZSQi9YDv3cVGQCawB2degjAgzWfzOJzzs8p9rcE5Fz8Bd6pqVpDCrlRKeA52AKOBWu72z6nqRPd4c4DGQCpwCLhZVdcE/ItUIiLSDHgT6AiE4sye+Q9gBMXkCRHpAfwP+JuqThKRVpTV/xdVrbIv4JDP5wbALOCpQrYdDzzrdcyV7VXScwAMAt7w2W8P0NDr71OZXsCTwEOFrBsM/IhzFa8VsNwtr+b+Qrrc6/grw+sYzkFb4AS3vAlOQo9xl+cA8T77fOH196pIL5w/nP4AbnKXQ4H3gP8UsG2ePOFuOxsn4V/plpXZ/xe7hO5S1d04/7jvlnyT/YrI9cDxOP+ZTICU9By4+20AWgYhzCpPRNritDz+rvlaDaqaAfyGc65MgOQ/B6q6Vt3pmdWZEGo3cNTIXTjJws7NsekDpKrqBwCqmgk8ANwgIjWyNyrkd9Q9wOc45+Mopf3/Ygnch6r+ifMXU4PsMvdyx0hgoPvDNgFUknMgIscBxwHrgxRmlSUiYcAE4B+qurmA9VHAOTi3RUwA+HEOeuJcmt1QwO6XYOfmWHUCFvgWqGoSsBE38Rb0O0pEmgIDcG5tFKi0/1+q9D3w4ohIKDAOeFxVLTl4oJhzcI2InAEcAW5X1f1BD7DqeQZYoe79VR9tRGQxoMA0Vf0m6JFVHYWdA0SkMfARcGO+qyPjRSQFJ+ncE5Qoq4gifke9BjyqTn+e/LuVyf8XS+A+3JZcJrmXO/4J7Mi+dGIC7xjPwURVvTtowVVxItIbuALoVsDqDaraNZjxVEVFnQMRqQV8BQxX1bn5Vg9U1fkBD7ByWonTaTaH+7NuhNMRrbDfUfHAJ27yrg9cKCIZwGLK6P+LJXCXiMQCb+N0jFIROQWno1RBv6xMANg5KL9EpA7wAXCdqh70Op6qqKhzICLhwBTgQ1Wd5EV8ldj3wEgRuUFVP3Rb3P8G3gC6UMjvKFVtnf1ZRMYCX6rqVPdye5mo6gk80r2Mkf0I00fAK+66p3AeCfgh3+WPK1S1oHtLpmRKdA6CGaABYAhOv4TR+c7F88DvnkRU9RR1DsKAs4B6IjLILR+kqouDGWBl5DYmBgBvisjjOJ0DJ6rqcyLyHR7mCRtK1RhjjPGTiJwGfAwMUNWFnsZiCdwYY4ypeOwxMmOMMaYCsgRujDHGVECWwI0xxpgKyBK4McYYUwFZAjemEhORTBFZLCJLRGSh24O2JMd5V0Q6BiC+ViJyXVkf15iqwHqhG1OJicghVa3hfu4HPKaqvTwOK4c7sthDqnqxx6EYU+FYC9yYqqMWcABARGqIyPduq3yZiPR3y6NF5Cu3xb5cRK5xy+eISLyIhIrIWHfdMhF5IH8lInKVu36JiPzkloWKyEsiMk9ElorI7e7mI4Ez3asERx3LGFO4qj4SmzGVXfZIdxFAY5ypEQFScQaiSBKR+sBcEfkCOB/YrqoXAYhI7XzH6wo0VdXO7vqYAuocAfRT1W0+628BElW1h4hUB34VkRnAUKwFbkyJWAvcmMotRVW7qmp7nOT8oTvXugD/EpGlwCygKdAQZ1rDviLygoicqaqJ+Y73J3CciIwSkfOBpALq/BUYKyK34UwNC3AezvzJi3GGXq0HnFCm39SYKsYSuDFVhKr+D2dWpFhgoPve3Z0VaRcQoaprcSZmWAY8KyIj8h3jAM4EDnNwxuZ+t4B6huDM0NQcWCAi9XD+YLjH/WOiq6q2VtUZAfmixlQRdgndmCpCRNrjtIj3AbWB3aqaLiJnAy3dbZoA+1V1nIgkALfmO0Z9IE1VPxeRNTjzIOevp42q/g78LiIX4CTy74A7RGS2W2dbYBtwEKgZoK9sTKVmCdyYyi37Hjg4reAbVTVTRMYD00VkGTAfWO1ucyLwkohkAenAHfmO1xT4QESyr94NK6DOl0TkBLe+74ElwFKgFbDQvYS/B7jMLc8UkSXAWFV9tZTf15gqwx4jM8YYYyoguwdujDHGVECWwI0xxpgKyBK4McYYUwFZAjfGGGMqIEvgxhhjTAVkCdwYY4ypgCyBG2OMMRXQ/wMU9WCnv4aIbAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 504x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "labels = [r[0] for r in rows]\n",
    "uncorrected = [r[1] for r in rows]\n",
    "corrected = [r[2] for r in rows]\n",
    "\n",
    "x = list(range(len(labels)))\n",
    "fig, ax = plt.subplots(figsize=(7, 4))\n",
    "ax.plot(x, corrected, \"o-\", label=\"Counterpoise-corrected\")\n",
    "ax.plot(x, uncorrected, \"o-\", label=\"Uncorrected\")\n",
    "ax.set_title(\"H2O adsorption on ZnO(10-10)\")\n",
    "ax.set_xlabel(\"Basis set\")\n",
    "ax.set_ylabel(\"Adsorption energy (kcal/mol)\")\n",
    "ax.set_xticks(x)\n",
    "ax.set_xticklabels(labels)\n",
    "ax.legend()\n",
    "fig.tight_layout()\n",
    "ax;"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.16"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}