{ "cells": [ { "cell_type": "markdown", "id": "cc3ca078-f681-46f1-9035-44bd4ae01b18", "metadata": { "tags": [] }, "source": [ "## Initial and final (approximate) states\n", "\n", "By default, the NEB task performs geometry optimizations of the initial and final states before creating interpolated structures (images).\n", "\n", "But it is a good idea to optimize the molecules before starting the NEB, to make sure that they are in reasonable positions.\n", "\n", "## Create the systems\n", "\n", "Check out the GUI tutorial or the related \"building structures\" example to find out how to get the initial structures used here." ] }, { "cell_type": "code", "execution_count": 1, "id": "c05063dd-baf3-4ecb-97d4-a9cf53fa7ec4", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAJYAAACWCAYAAAA8AXHiAAAcuElEQVR4Ae2dCZQdVZnH7616S+9LurN0JyQEAiQkomyBDMMqgo5BcDQBRxRlGFEQjwuemTNzmCTIHJijRwygMKIHAQFJ8MyALEpkiXLgsCSsWUnEhCbdWTrdnU6v71XV/P71+sUmGzEG+uX1vefcrvdqe1X/+vX3ffe7t6qMccUp4BRwCjgFnAJOAaeAU8Ap4BRwCjgFnAJOAaeAU8Ap4BRwCjgFnAJOAaeAU8Ap4BRwCjgFnAJOAaeAU8Ap4BRwCjgFnAJOgYJVwBbskR2cB4aeC7x3H/qs0BgbvXue++YU2DcFBJS/51UjlkXD6p84sWcx3JJ9U2AW0CwMjJkd1NV9qjKZPPVY3z+EeTUmCPq3W/vwsuZm253bV4Q1s1iw4i/D6r/owF/OHFR1dWc3ptMnzIuiIz7ueePHGTPaWFtlIhxgFIXro6jvjt7eP9zY1nZ5R86yzQbE4i4OrP2/voqlwoaGU88xZvovfP/EhigaB1CjQs+rN55XweKkBS5VrNem1UGw5uKWllNeBLeit1w7BZr7r/Lw2nJODNXo0VNOsvbQ33re4Q1hWJX1vKoolar3ystrvKqqpFddbWxlpWGeySQSo470vKN+N3bsH47JuUPBVbylqE/ufbpssZWvrTXVWKZ7rT0Ed1caeF5ZIpmstGVllaaqypq6OmNGjjSmvt5Y1k2m0yabTNbVhuHYOydNmp8eOLai9RgOrL+OPkC4KmXMtWEq1fAtz6s9LIrSgbVJ3/eTJpUqMaWlCcAyZsSIGCqBZWpqjMFyJXzfZBOJQz/S1XXiZ3NW6ymC/OIsDqx9u64ApXSC8lE39x1/fJT0vPSnjUmohYeGETFVZHzfApcxJSXGlJUZU0GYla+lpUbLrLWsaEo/z3bs84yibSE6sPYO1iCg1JKLUp73rS+/+uq8pbi/D0VRhuUZNOxnL6q91L+wYlmarx5rUfkbMDc4acKEOdU5q1Wc+a3hnsfiImNAds2MM38OFurarPJTrABQ376YeOpb1k6elkxiksx9VKWncjWKukwYbjOZTI/p7S033cwO2FK1F94yGZCLmRN8WeKyEn6jeMtwBGsAmqkQJWh0fWOrgTVZOHClNX8eUAmof7/YmEbiqcOmJRIjcXdlge9vZSNZq22qQLINaNqBpxKIKlinke9pk0zmYOrpMUZV0IVhG7ttS4ThJv1w0ZZhBlY+Sx5DY8aNO7nU9y9Mr1tn27nCwJQv00Z73j98wZjDLvH9idN8v94kEkmC9Cw5qXbT19du0ula63lb2GArtRpgygGnBIASgBaY/v56Qchnj88h83tZ3hZpe2NaWoxpArXidIOcoBlGYOWgGj16xijfP++TWKGZYVh7fBBUVzY2LllDy255b+/qB1tbn/mQ7x/1dc87YlQiMQI4EjFQOTe32YbhO142u9xks23dqVRPWRhuII4qByyiduMDUcTnPuDbRkxVDlhJ1g2pPbjF7aG1HWj+zm+amhYC1lw+5yDXxSimMkzAykN1wWzfP+ZGz5vcaMwEruMIoKihVk0ntJ5eUnLIl0pLxwBDCW6sFFfY54VhB5Zmsw2CJi+KXmfZigej6JmfhuHUI4mrbjSmmcipBDfqA1TE8n6mXUDVzv5pGvK/q/nAhrXqweKtzUTRqp+zQOUvkX7ue9H8HQZg5aBqaDjzX7BCPzXmSC5yfdb3q2wiUW+TyQpcGrNCEwVBja2vn2G2bVvh9fc3EQd1YWXWActr1GX/a+0LN3nezC5rb/p8Z+f4L3veb6LS0nuAJQGchGy0DIGHbTrZYSXfFeTLNapTJ5vxvK2pbHbl91pa7l1hjI5r3iD3WzRMxSdS5GDloBo1auLZ9OH9NIrGcCHLAaoikUrVksysoNKjlzS+Ln1vb0RrDuNjJpnW1mVRf//CfmvfeszaJfOtvRBALvqGtWM+nUhMJOYqxZJdjavMROn0/wVhmGXDTtxmO0F5HevmwIoiPyI2S3jeNqBqur2l5XfX5qDSiIjiLcUMFhd4QThunC3NZqt+bG09diONi0p7iUQJScxyZcPjLLkSmipdXRZQ5LbSDIH5tN/dvXVRa+tvv+95//xdum8uSCTGA1RJ5Hk2lFsLgi1hR8cnEyUlW/2KiufZdzah4NxaQVVBVUaUED3q3oQVu66lZfnNpDEAcF7RukDpqFLMYHEBbZDJ1HyU1hnxUBIL4fkefo/+F7LgXmytlBkvL2dNAFB6IJu1slx+f38SUM46x/crZiYSUwAqMQBUN8H7tiibbfbDcDVx16ubOjs3/ySZ9B4rKYlOx31iHbejKz9v1NoM/khq4d6NG82mXCtweIzHKmawQEXApM611o8IouSSmJOVZaGGwMJ/FlVA5cFSl0wiYVmueSWpdHqqQAywOjaT6SB+2kDLcBWW6dWNUfTmbQTjt7LTjVuVdTDmBer340+7/JFbtkXt/gafcjGDhbuZg9e67RiCavro+gBNWW8lKbvjqryTMuJ9fXF3S/xZmXLFW2qw+X4IYD1hJrMNoJpoIa5k2WsbrV19Wxiui4ECFtY9DR0Xy73JTGnreA9MtVCVvRZ3TMU5vqscDGDtuBOBK5S/YO86iV2/zPGONw/7S8y8TGjrer2I/hXTSd0OGIKkg3xTK0nLcqwT41kwZLJY6npRV0yu+0XbdBDAv+4FwQq2W/GGMWtvCsO1D7Fg4yCggGaxTKGK4HIFBQoVLDuHxNBUQJrNf7v+5VWgyj7NfGo4bzcXcVbchFfHzLxgCcsPN+dN7e9bPTFb0hHZkG4Y24alqgKscqyUAvkEn0cDUylg2Rio7u6QZd3M34Lr3ApQj66OoqX/HYZN9xiMWw6oz+LWFgLRDqBY5MpgBQoOLIDyBA01tgInjxtXeq7vp0sirvP69W0cfDwfyBikkuuGEVALaAFCDtbDmMlmxrRKv+Hb6Wj853t7R6fWp39NpmkTbqoWUMqAJoXF8vkcAlM3Lo+hnraE+TjN/j6q8lDt/MSGsL6+9vzm5qaVYM3PnYFei1lreLm1wcDs6+eCAmsWsACU4EhcPWbMPx3ueeePCsMTaoOgkoZb9Fpj4xpyBo8tz2bvtps3r+XKe3PNHDMPCwVU5iP+tI+NMId8IWXrLyz3xqWiRDXRVTraEjxnerwNgW8qsTQpoMmFQtzkgBvcjuWq4Cdz84MgA2A9WWs7CdBX/7G5+c5V5J1SWCqisRzUTF15DwUKBixBtRCoTmtsPPaEKLr1GN8/6XAOvp46AjdVSU1z1wJmZnqjtd98o6HhGtvcfDM5IXMaQFVFNd+ttCM/VuHTJZOoYpRUKugyvTboD8LK7smJ3qrFNA3X01dnoSoErgy70lCXDn4BbNk7RpFKxx4YRtE6uF31HRYwPVpWkqkr+6oA5n3oy5wB9zdtzJgTTrT28RmM+R3PFR5Lm66Rtn91Mml95QXUNaIe3UyWfpLAPNkT3j9/y6jy8X7dzDJ/JEBVhBmvNOqKPNsRBOHWbHdic7jVtEYbmjIj1q9PpO3fhWEScBrwq40AVsf5VyOA+vSSgONjLfuwalvopnn7K83ND90+XBKaB5qCQrBYObhra6tHWbtgAlCVAVUFkXU1SaVq+lx89btwNwI+yxIAeUFXV+R3d4ZTvNSFJyZHk4cqj/o9P+yM+mxXpifaHPT6G8N2rznauLHda721PdjwY7PFbGloaPwfz8t8JYqa0LEdTkcAWQ3xVRnHkMSYBbi/NjqJW69sbl4KVKejT3GOPjjQIO28vyEHCxeoYD2YmEp9e5TnTSRIz6QZw5LGQpUDk6/B49xHFafHyZgr6eRjwbBdfkPQEUwp6zbLMiXe9mybacn2e01hp9cUbdm4xbbd1ha23orXG0gNRHjODZc3NDQ8wM99HYBOof+uDqDQRDJ46/npx4Ng8/yNG7veGOjPixsKO4vmvr+3AkMNllVcZcaPr01mMpencXWYjoSa/nHXi9yfrJX6XHTri1LkSjZpKEJvn7F9Xf6U1Dvmoc4+DFKX/1a0ranFdvy8Lez4CadOF4qMYZy8VIMAV2eAq3kR00Xjxo0bEQTBVM/rTQRBllvhW7kVPh5nzOK4P0/buLKfCgw1WGqeBZU9PVOwQqPxSyEmgg4XDfy1hjtA5f4IfbBUgkxVfS4A5iX57CWjSeluM7pk/ZYHe7r/tSPswBopEzoYqHflmgQXOzNRU1PTVqZ/pA4q8WgI1in+TuJBJ/2+fNSFHfKS9rxxgklD4XKdLrTL8HXbsUyh0uL9/bELjK2V+l80T3cmkIgKiZBm1USvANUdEANUAkeGLwZKIO1cZInYOKaPdWPQNGWbOD+lZa78jQoMNVgyLXJ7VfTmcVuCOl24PQGoOkBmK30r7boLYTtzOzpydRtrdXXhEklkktFkrtkcekT2xr8vB0ne7TFrr0XQad183R2Ee92BW7hnBfAnQ1rii0mQswywSAyo04VYHbAqgKYUS8WdoXFwVM1nRsvl4itg6+3tMW1Yrjas3BYT36EgQGR5XCkABYYarNjtRL7/JsNatrcyQo6IO6oAsBIsVgJwIoL1PiCrYVqOVWMQFN0wGdNNJVcVsY3dEEVPS8uROfdWALK6QxhqsGSxfFpq5JgaHuzxvC+8zRjfUmsTGivMIHRDUsl0MW3DWpWwMokmwzhg0wtUXdxbtYocwSue94gu5Zm52EkfXRliBYYaLJ2+4LJYp3lBMvmZVu5AeIsb80iNM5o8Mj1AtQ2QsGIG4OLB6Yy60w3tmXbPS63MZq+9v6Vl9SwAjVMX2qMrQ65AHDwP+VEMBN2NY8ZcQiT/Cx3PCHXp0K83CuhqqAKLCF0P8wQzE3WSmV+dzS56qKXlXKDyHFRSrXBKoYAlRRR4B42NjbOJx2/knqxGEqamhloLTIq7yMhbLFncciRo/9nzPT1Xz2lr65yXs3pxQ6BwpB3eR1JIYOlKxHCNHj2a3h3vCoYizKLvZgKp+HLFXLQJN7POc13G3LIpl0HXNjoHB5WUKKBSaGBJmhgufVDc9OzYsQ1ZY47Icv9Ua1/fm6atrUPLIAnbFRcHVU4H93cfFBAzAmyXMod4ipm7XbbLym7GkClQiBZrsBg6vnzVfOW9nIWSEq44BZwCTgGngFPAKeAUcAo4BZwCTgGngFPAKeAUcAo4BZwCTgGngFPAKeAUcAo4BZwCTgGngFPAKeAUcAo4BZwCTgGngFPAKeAUcAo4BZwCTgGngFPAKeAUcAo4BZwCTgGngFPAKeAUcAo4BZwCToEhV2DOnDneU089xXPnIl7HEld/wYIF7lEDQ35lDuIDmDVLjwnffRmArOAelVAIT/TbvWJubqwA4PBUMBtcdtllU6dPn/4lXnxw3JgxY6KSkpL1XV1d97FML0TQa4a0nnuuhePmvRWQ+2Mte/HFF8+55ZZb+pcsWRK1trbCUK7wZo1ow4YN991+++08ADHiUWJ6vr0rToG9KJB3fyeffPKcefPmRY8++mjwpz/9KQNMWQBS5RGtcdVbNh7Q+nzfo8vcy0+5RcNIgdhS1dfXH/mJT3yif/78+Zknn3wyzFsrHl0OQ3otRzzVCz2iF1544WLpw8eCgEsn4MrQKZB/9ld+Gh/J8ccf7yteIob6OI/MTGKleOtLv83qVS8Ulg2e6hqGlZWVn4lnFsgfB9bQXAhvwNUp2N5RNU916dKlcnPm9NNPv+Css86KjjnmGEvAbjo7O82WLVviFyjosLUORdfQA8JpM2fOLOOzzJiLtaTMMCp2AKj8KfvAU3P++efXMGOwC6v84he/+Nu77747ArJw69atebfHizp6I+CKgAx+4hLo8+OPP/6mwGJOQQTxLt2Qv8Tv/1RWJFq4cGEwY8aMaRUVFd8oKys7JZ1Oj0mlUhaQmn3ff7G7u/sXvKnjH6dMmXIuMVZ2xIgRiZqamh3uj/V52WzabN68We5RNWxrazME8G88/PDD3ezfnz17tt4rNKTFgfXByB9DBSRVkyZNuo546PLa2toUYMWQABavYkzUAszRa9euvYQclaHyoo5kgnVjqLBEO+DSIbO9efnll822bdtCIEu88sorv/5gTmXffsXFWPum09+yliUfZevq6ip5OcIirM9V5eXlsJTKAk+I5Yqqq6ujqqqqkJrFYsm/xa07Be20/Hb5bS0HRC0LyG2l6Or59U033XSf3GwhWCsdsLNYu1y2AzuDi+2RhwoOO+ywGwBoOu6uH0uUBKoEUBm5OU2pFgvmPfLII/EBKFBXlZsDyndZq4FWYchb0PwnnnhiAfHVV4GN9y0UTszuLNaB5WjnvfmKqQDjBNzcFaQOsoCVgisr9ydXiJUyvPks3u6Xv/ylWbFiBe/57DKbNm2K6zvvvGNaWlqMrFe+8Fk5LYsbXANUFwJU28CyuJmYX28op85ifQDqA9VX5LpkUfJV3/NwKRi/5pprzKpVq2LrRSvQrFu3LrZmirfk+oDIYPHi7bdv3x52dHQkWOdZltmnn35aea9ckusDOJ99+QkH1r6otH/ryC8pQCqlnik4qBryEoOSnwowuby33347jpv0kk9SCoYgPrZoWq+HN8qqFSiXicXjNdnZSOC99tprdwhU3G3BWKq8VA6svBLvzzRiNEIprmukgnBqnD1XBl0A9fHidKxPHEORWjDr16/nPerZ2AVqSmAeu0XirKyWCyxKAhA9rNuPcLNPK2CXu31/Dn//9+rA2n/t9mlLXJmvmEkg5avyT7JCmq/gXO7u0ksvNXfddZdpb283J554ohk5cqR57rnnDAlS5agSSi9oPbbdhBu8Geiu0+gHGga7Nhv36cje35UKpxnx/p7nUOxdmfT43T9A8lZpaemhwKGUgifLo/yUAnfFTfquqayaYimC+xgwQbdmzZq+5cuX/4jAvwe3uB6XqGbjJj7jBQt3/JWzWAceOUs3jb948eI4mB4/fvxVQFWpDDkWKm4NEifxItlcg1ytPcVUslSar5iKXJbmBVg4DwBfowvn3wYfptwfUBWc+xt8jM5iDVbjb/8sKxVfcGKraViiHwDGuWTcYyu0cuXKeCpLBWw7qlqFCuIFm6yWYAPELDUBVF9+9tln754wYUKSgF2wav8FF6zvLJ3LY+2syP59z3cu66InJk6cOBe3t4QRCefyttiATHtENaNGjVKHcaCAXfGVpnJ9xEyxxZLV0ndqhuUJclmvMO9eYim1AvvYt8AqeKgkoXOFUmHXIlA8LqwFhvhCHn300ZZYJ+5E3mn12EqpZYbb+yhW6r+Il06SVSIBGuDefFxaKGCwQBrlKZCCgZahxzyrmIrVIiyWfivEYiUBqw0IL1m9enU/82QADgqg8to4V5hXYmD6Xs13Bc30x3mAFAo+AYV1qia7/j2s0lXqfgGuLLDoVi3FSlmgSigPReJzPgAdSj1fPyegSJKGSpTiCj3lpOQGSUM8D1RfYzz7y6wmqAqy5adz2FNxYA1SZhBUyTPOOOM0LM7JxD+TAUEQrCT+eer+++9/dtAmhljqIizUPPr8jsRKhYqdsDwecARYHV/uDaieB6r/4PsT2pbO6EuB7jLqFFatoWp2O3U5v/FzgLqLz3J7ByVUHHf8WlxNh33J54SOO+642Vig/wSWqUoJKHc04KryeahFWJbrFi1atAorM591LmT9HVYKixORo5Lr8zdu3NgPWNfTr3cdAgsUhR47gm/6COuBa6rEZ5/LGIe1RZ8HykELlY7fWSxEyENFHHUjwfY3BYriI2pEjeGS1cKaMPEtAEQvvvhiK/DVy0qpVQcYsZUiIPfV1wdMvyNG+w7x1bIBUHa0GPm+I8c1sCw/0fXIA3VQxVT5E8hPHVi5i6xhLd+kFXcjLi1LjOQBjGqcxFSLTnHQQBogeOaZZ3y16AAwDs4JxEPlngjKFfBvIrF5HZbq5gGR32Wl8sIPTPMg6WucTN1p+UH7dTi3CuOWH88/COmHm4zV+QHwaFiLz9TKBQKY4IlHGcgqCSzcnE9eKlLeieprYB6QebrJASv1K4a4fAcaNlBleVTkAvdUZJXkGouu5E++6E5sLydkBx6mEacOcGER0FxNVTZbgbdVMK24SoG44MKKGSU5VfV58uTJVmkCQIo0XurPf/7zKurngOpz/K6g0j+sLNBB15rjmA9IGVZgKZZCtUjDdwmcy66//vpanokwDiM1kyBaJc4zSdkByGLABJkslqyYLJVuxaJlp7HmsmJrialOocX3KzZT7CT3tjcrxeLiL8PJFcYjAc4+++zGU089dR6B+jlHHHFEBVBoeG+13BxgxV0q6lYZSGAqqRlXwZWfr3l0sVhgs7jBWu77y9xzzz17CsiLn6LdnOGwsFgDlirkOQjTjz322Jd4astlwDD+wx/+8AhcXbVacbJQGh+VH86ijmEF5MqYKxeldTTVd6xTvAxLp218BuXlA/TdSDw8Zw0Hi2Xnzp0b/fCHP6yje+Y3WKlRtPYyxEvxuSueotVnBY1agbQO485gtfrUEtRyWS+5QRUBp2UCUOBhvbLKWw1PfPZ81kUPloawcPoBUH2DcVGjcGkZXFhSrT3FSBdddFHcCcxNpAYXqdgpjqcY9qsBdrH7E0yCTEWD9fRdMGGtPD4ve/XVVzvyubB4Jfen+Duhr7zySo2Iiw455JBPAYc+e7I4iqk09JeBdPHoTUEli6U4Sq5OWXdZJblAWSYF7bJe2k5wAVXEVDHWo3AUPc0NDUyd5Rr4pyrqBGm+7++8886bedRRRz1A/1/q0EMPtQJIaQSeN2Uee+yxeCiwLJJcnkDS2HNBR64zTjlomeIpgaUAXwlRqiUhuoks/IeAM98Vc1Bnyw+koS1asPKuCUv0JYC6A5cYEj9ZOo3jYSp5EclFxYDJEtG3p1uwwt///vd6mJkvAAWVWoQCSwE+YAVYNQ9rpvv6znnjjTcW5QHO79NNi3c8llILERn1I4mtfowbjJTsZEhLDJWszgAk8R0yco0vvfSSrFQEVJ5urWJ9PeElUP5qEFgaCuMDlW6Fv4LWoINqD/9FRZluwDrFiVAszlcJ0suwNgGWx5P1URFU+SLIlPB86623Aj7bqVOnrsXNPSpXB3AJ3F2Cvr8ELlIjOrNYtUVM/x6obneWKq/irtOibBUSSAfAI+tyOoBEpAQ08iCOoQZLkAcM8BTc+8RWbbQKZwPRUlzmdNzjDLar0TZYr3b29TjLVgzsoyDv5xs4tiGfFCNY8GIjWnXVtNomEQcpyI7v31NgrqB9p6LWnR5k1nznnXdewPpLsXgJ7rJ5gfVUdy4yd6pF2Xm888nu7/didIVqmVlA6cBSrcF1Rap6sAYdxYbRmXG6QC5QATsgBSzzSIDeyecX6aBODdy6pfSB/vEG19jFMs+lFRBhb6UYLZbOV1BksVyLSQccp9umyK6TLfDibhsCemXbY7CUs9JzE7hp4Sm2ETh5S5SfMsuVv1aBYgVLFsUSO90GOJfzaCAF8Flcnp57YIAtbvUxL0PmPPX6668vuuGGGxYpRVEoDy77ay9koa1fzGB5uL3VBONX0Kr7Gd0uat3p5odQ3Tm4Pg9XmGKdN3nk4tfUIpw7d26hXZ+D9nj+0u4+aE9hrwce9xOOHTv2JFp038MVnklQn5Ab5HMnYD3Aszu/yx5aqdLCZc73Kue+Lyx2sKSE4qY42CZfNQV3OAbI1PG8mvnvaAXKjnVyX91fp8C+KSBwdvknUoJzd/P3bZdurb0psIvYe1u5CJYNBkxuz6UNiuCiulNwCjgFnAJOAaeAU8Ap4BRwCjgFnAJOAaeAU8Ap4BRwCjgFnAJOAaeAU8Ap4BRwCjgFnAJOAaeAU8Ap4BRwChSIAv8PKKSx66H+wEcAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from scm.plams import view\n", "from scm.base import ChemicalSystem\n", "\n", "initial_sys = ChemicalSystem(\n", " \"\"\"\n", "System\n", " Atoms\n", " C 1.0065881512963437 -1.417955004240563 0.12885420257373528\n", " C -0.08873018725926754 -2.1350000996605933 -0.03448649693154524\n", " O -1.5081108856471217 1.002080375323799 -0.2045570136880532\n", " N -0.4803377947492119 1.514182931120784 -0.052917508314934095\n", " N 0.4992559331764258 2.001851603909277 0.09160944741412266\n", " H 1.3565548327663348 -1.1143852514314225 1.1071131753303691\n", " H 1.6148659535066119 -1.095826278535031 -0.7065218137603257\n", " H -0.6985763255773102 -2.454234856430313 0.8007371044735704\n", " H -0.4402537210929026 -2.435675244924285 -1.0129788309893883\n", " End\n", "End\"\"\"\n", ")\n", "\n", "view(initial_sys, width=150, height=150, direction=\"tilt_pca3\", guess_bonds=True)" ] }, { "cell_type": "code", "execution_count": 2, "id": "f5ba9744-4d09-40fa-b391-74ec6b926a69", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAJYAAACWCAYAAAA8AXHiAAAqbUlEQVR4Ae2dC5wdVZ3nT1Xde7v73n7cfqW78yRvIBGDQZ5qwIFFdwZGVxLUFR1liWLQBRx2PrM6nyS77u7M+tbxAX5EcRTZxM84zkcd5GESRqIIJMGEBBISHkm6k07S79fte6tqv79z701aCMgjdN/b1vnkn6pbdarq1P/86v86/3PamKhEHIg4EHEg4kDEgYgDEQciDkQciDgQcSDiQMSBiAMRByIORByIOBBxIOJAxIGIAxEHIg5EHIg4EHEg4kDEgYgDEQciDkQciDgQcSDiQMSBiAMRByIORByIOBBxIOJAxIGIAxEHIg5EHIg4EHEg4kDEgYgDEQciDkQciDgQcaDUOOCUWoNOXXtWu8YsGvN+j4fGrA1O3f2jO/0pcQAgrY4ZE44B1NjX1/HQe/HzY+tG+6+FAy/SAa/llhN17XIAs94vPn3q1E/PcJyZ84ypRVLFnTDMPtXe/r79xfOAC4nmRBLsBENO6d4kAVYeVPX1S+uqqv7sujCc+W7HaV7iOG1JY2pgWKUJQ38wDM1jjjP0nXT6r36wc+fOUWPWAcYVx8F4Sjn7J36zSQCsPKja2v78cmyqr3re4gVhOJVubYDqAsepNY5TxX4CCRVnmzW+3/54GO77aEfH2x/Mq0YnAhecOZUFZpdzKYJq2XWOM/tu1527IAjSOaRU4LrpMJFodisrG9yqqhQU57fheDznebMWue4bN7a1PXAd6hBQPQLirP11Kj807qV7yomwNt+pvHfJd1oZv2weVK2t867yvKXrw3AJAJmBdJrmeV6ricdbTEVF2iQSjolhzju8aQ7IjaIAR0aMz9YLgiPYX09ccPDg2x460VOyvTYWPjidlze5BnKgP1oEJNTrIuqeTMXKcVjPTU527o/eu6wqlCuwbLvb2qobjTljp+Oc12TMXDpzhhuLTQNUU01lZatJJhOmCi1YUXECWMPDxgwOWrLgyuX27Ovq2v6xTObokwDoqDEdQyfvQQHO8NyN0MkAJ8l0Ipwxffpynnx+FU5DMhZbOPDss+/u5/qCyrX3erlgPXlzSvwo33I5lmV8+ZtyYVi7yvMamsOwAlkUi7luzIji8QRgigEslGINpnulMR5XSGKhDm1h34NCz2uaE49fdI/vjw6F4VuOYeTvCcPMdsfpfToI/O3GHNxjzJ19YKr3hZxay6Ei4KydFmttven9rnvaXwZB8znYebXGJOPZbCozdeq2w9h3vzam73vt7c7m/L107eT0TMsUWJv8M880ia6u+AfCMKEvXx1k1Z3Lnuc5gMuxkkqgEsCkDgUslYI6NJkMBpbvhdXVibCqqiUZBCbp+2YG9HYII3/ABEHXcBiePxAEwa4w7NuHV/lUEDgPG7NvvzEPtxcB19Jy5Xmuu/Arnnf6ecbMwgttoD11PK2GbTIFpBvC0DsjDLuva2t74i7HufcGAHYsD8zJB65yBBayx/jHjjUscV13TkgMwXH05UvL5NgX+UiuEHKgvLQSsFS0FUmCiRwn5jjOMcfzKsJ4PB46Thygggv+D4JqF6oCZFWArhlgvo0tgBvm3Hyw9rYu3/9vOysqHulx3Xve6boL4kiqHJ6og/PAPeshAUtNNiH3gRo51/he369Z2to65/2HDjmPTEZwlSOw1En0fKxGgMKGoatHsXsyHB3leAZSxyusUAGdkFR5KaRr82RvZEJ+95q+vr2EJuoAWyV4TWkLAKrYJlCXccCR4DxVLeCqANwMyDQJbJ63FwC3cbJa4Y2YwBSPy3Gos5KzAGpHErMgLbPGTJ1Py9ZNn373OQcOmB5axDu8LAch3+wS/79sgYWYEWRgL9Z4gcJwCEANAaZ+k832oeqazPAw/UWRdJK0kfGOCrRqEZHDmRFAM8Tvn3q+L3Y08rtOJJAQnqjh2hQgSyJ9koCuim0lQIoj4aoBzwFu8i+AaqHrOJWuzsXjKey6ahwHxzoPY+w6+3zaEseBAFwzZ+dyvd8CUCvywdqica8Wl3cpW2BhrwAVqb+B4xSG/YCnD5AkkQyVdCLhKVMHyJzjwBKohoYwo7MhAJQN1RO6rg9IHvi3XA4H02ldGoZ1dWGYxjprdH2/nrr13Cf9fMABtjovlfqZF49L1VVyXs5DHGBVWOchhWVVXZ13HhTuyAKluJpEQXrGaYfvujOWT5t23yUHD166YTKNBJQtsBzH7Q/DHNKily7roatqkRopqApgJZBK6myHDsyyLzsnwTnHAiqTyXFskHq6rhdwHcgEwZe/HQS1W425gmOPD6O+TkOgLQQQ88MwfrbjNC4CcG1hWF+Ft+dms01IpRxA6aX+WZBsPQOwBGLXxs8U5lC4Q86DPgMBqwgwYmmSnKjWdJjL1XyYGwCs5WwmRylHYElMOY2Ntdu6urr2Oc7w3DA8itqqwUyvBCQSCQKVVJ8PeIbpVMSGoYdNDKAFdPAoJLU5yLGDeIaxyjA8918SiWMce/BQEHg7MNAPcs9Hfd/dCND+nzG39BizNQRw0wU4QHNBXd2c641Zwj18ICW16gMcOQ4BFOIk5IOzsrEKdpaVWJJa+g3YCHkgTo379ra21cmODnTyJLG1yhFYdKDxNIg8derUHyAB1hBBB2wpgBVHBHiAQ1a2z3YUGqSz0UkG8QGwgkAqUOeGqTPAsWfp4L900ukUJldKBnkrp0VSVx8KggzbRdT/yVEk3pM4B/uCYNrvqqpW+Z630QmCHKAWOAjpW8o7D3q+bLbCfayk0r7MOjkAklz5Msq+mzZmSiO/AdYanaFGeZdyBZaVWuiRr9PZq1BwBEnb6WB0EEWdGobyEOUdDtBx6CKTgGK2UwED+yOEFg4j0eJIrwQG+VEMcxtyQI3JA3RVl0dUAJ5WDHvTBCia2F6Ed3dNLnct9/6F7sMz0WvWgcg7D3pmNjuIqquzxrpAJFUor1A2nlSi9sMwIJ1nkGeFKcIiZ3CT/X+YnMiRMi3lCix90V5HR8fRtra2j4aOv96EffTUAUSBMEfP0dl0GlRUg3lg5YUBlrsZBjDtpqdnKlLlHry/aox4S4BMXqC8v2Qx5IAKSyjGJfByT/6FiwFUGurmWVKpcgRE8kh7AZBCFQpTJC2IFE+TxBKw5JnKeZA0pY2Aq3skm+0HVCrLy15a6S3KFVjqAGi96e2INzSm5jl+7TOh73QaN8xisg9iSfcjJ9LUkRqsghShJ+wlMeFkeO8O09+/vX10NKzHuyMI2kDHN1BPVA8g0iLq1lnAOU41QKsGLIpzKfSgWNU06j4NadsHpQBL3nkYHY1xvcDUyKB3CptLalhxrJDfPgCT89BFm/oA1tE+33+0AKzyV4MwomyBxfe/LjAVzpy6XNtXGkcucCqcS51jyftNX+wxTPchjKweXq+GrqxhW0Unx7mGsR93FDuqux1J8+GBgac2M+TTiGoigm9mIY3mhWHyzUiZGWFYMxWqE9BeCLga8DkjqK2tDysrHwOkXdxaHqmchwT1PRsIlcrNZocx1PNtkP2Xy8mhyEBSm320r5/4V8fuI0fOHEKa0sbJMbxTphJrmTog15Rd8A/NzoJUdbzKryZdpmboI2E7XntP1b8TIejD8R+i3hE63dpWDLWEuzHo1wdBxzcOHzad3EMpNANUeBaiSI32izxjDjfguZ2BlJlLcGsWgHsD4JkBAb7q+iBo8UZHhwg5yBY7BHAJe1nnwQVYsp/yzoMA5Lq9nK/gmBwLOQ85tnIe+jl3FInV8V1jPgfIdvLcyVHKEFg2DyuXcptWNDmzrkq7TX5NrNpz3QzqpTeoGG3z3KEF1wStv9lJ+KHB9zN+LJbFnsoe6ujo2k23CT0Uex86E3SdIHZXc+x/qM4RDGyh8oG8XSY7SmTqAdx8190+LQiSl6BuP4E683EekDyMflPQtgDnhPOQV8XySmWfyUCTczHkO84gAdgnd9fX/+IH7e12SEftmRTFMqKM3gQJRM8Yp3mW98btM50lTXWJGaYq3uAO5wK/d7Tba/cf+/He4GcywF6kHAcU93nRUuQLz7PAYyOjGvWLnXbiqtC0trZ+iczCG4lVZInUM6I9lWubqVIP1UJ55yEvsXQ7Xa7A7jA77cSx9lzW2fnz+wpAL4CeKmVeigwsi9dYjpRZz0ycKd70dXPN2cubYrP9moopnh/Gwr7MQNju7xl8xt+y5Ebz/mfW2izQTWNAYHtUv8ceezXvLZ4VAaftaGvr1NtI1bkOcAEMRfmnoNIaOVUHpSDrPLCVptOHMYx2PYJUO7iyo2PDtycbqHjJsjLeLahqvdrLpprpy5vcBj8dr/Lwtcxgtt8fDjpjA+Fzf9tr9uzbaW2VTUSKXpciYBYliwXZoUPtK5Fcg7iLN8pG02A4Bj11BKwa0JQEbLLzRrgWj9VkenAe/rqjY9ftkxFUvLT98rQt9eKsNqvVoTUtpuWb05yWMB1LOilFHf1hf9TvifUE++7fG/zm60WpNk4vpDaJ3EOHDt3E9lI05700Kue63dh9+z3HeYLDjwKqR2XEEws7zR0dbftwe7tAtTI+di4kFSdNKQvjHbC4a81av82b8vezTNvcRrfGr4tXeEGYDX3c+e5gf9+x4NmP0sfOejvxYVz7pwgur729/X6efH9LS8siPMQF2ezIfyLT4gOx2LlhMvlWJ5l8Wzg8XEkoYhtDAS5fxfu49rZxbex4PawcgGVVYMpLXTrDNH+8za3z04kKr4KcuMHRAX8w6Ip1hQf+Yb95au8yc3FsE7nw48W85z1H6lFGVHD48OHH2YoejsVmLE+lZiYcZzoSy2cKWoYQx+A7kVQ/KuRgUW3yFRmfpVwcpJXaF59p6r8+2zTh61c5NahABxWYyfXGDgf7t24PtnxeKhBQFW2fiXonPV8SDIAt00fbQ1rOEWVQBEEnUfwBZhFh2MYqFxfO4x1OzlJqwEIO2fxc660uo4PkBc7w6j89zzQtaPGq/PpYnNzNbJgdHTBHgkPZY0HntXTN6JnmTHWoqBQKALMeKZMwhn9HpB8PsJ+2DTmxWAZgVcyprl6Yzr+q4leTr0w4sOA2hrmJrSv44nBZnBZA3F8T8KmLx8+eZer++0yn2m+IV7hJWhxkh4KBoNs7Ehz+37vMU1uXIR2wwUrt6y8AJvtIECjPq5fY6CDAGg0qKlI1dXWzz8rDaf2E98HrAesJtbFQcgpVS31Yu+j86dOrLme6TAu5MB9/7rk+TjjT/fR3Tnfr4k2xCr8u5jlZAgsjuWFvf9Cx+zfh4/9rnL3AV9IH+mYAVLjV9yWxehij1BBOQ5BIVMUymYrzudmvli1rdjZteiW3LY+6EwYscV2gendLy5TFnvfBqcZc1hgEp9cxL4pworNl6tQDmzOjQ/f31Jzd7FaFTXFQFeaIWWXCw0EXVkvXKlicxQiWwSwJV2ql2KbHkVaoxO5qqUOGnhyABcDiS9XgTZs0q3rylYkAllSfQBVc09Lyifme95nTXXcKwDL1BKVrSFxKsZ8kZ+WsZCVzq3LhQ0OGHE1lBgxhsg95B82xb2z3n7uvhKWVkIJqtvbTQdKcdyO13kS+FiAaRh1qwoV3blvbX5COvELRVNhRkh+H3uNVlfHW7xZUa2H6Va2tty2Ix786h/GPBvJLahmZTXte0BiLhelEIqwki9OLxf03JEPnmnS7OTt+ICDC7u7yu7Y/4/f8LQFTF8O+1Oyq53eC+EsbR8mh7yKroT/AzsIzHCXsUNnsOAum5y9YLWBNqjKuEgubikCn8d8KqGZ63nWovmw1hnstWXB1+OENzDKo1EwDEaNpGlYLSZaqyIyY892jwbRYLNaQya3/556uPhbgUNutbVbCPWIBQ479r3y/54OoRFJkBlGDtUFlZbIimay6iLbvhgoALOE3eYVNGzeJJUN9PaCa3dKyqgVQ1TGFGVsqTrK3w3IspoZZnZWaK6VVPOrqyK9rQDeSyan96hoTVKS8GRVhuDxlbtowffq8NdwLvTJu7X+FfC1WR1pZbG33/d4AIk1mAGBlkFgxhR2WqOKyZcXqk2c7Xh3DUIsJqltbm6sd539K7SFusDT4D8GUIOBZoenCY4EFqAQsk05bsLFyGhP2YkG9G6tv9P2/o7vAVb7XSrg7aKOaaXZjH2Jr9TLVrJ9A6YgbjwcIZheJtSy2ceMaecaTqowXsKznxuS9K5KslCFAccAhM46v12XUjAl2Un/FGZ6aPizJVVub32pKsZaNSWi5xzBg0YarWfRgHuDyV5e21BKqxOMBMht2aYIsEot41jA581mAFZvf1laTRusLgLzO5CnjBSz72TLUf4UHW8VFWd2kUFpiigzZ3jQlP96Rt7EENEmxIum31idiJmiaeVrYZ29VN1xc2sBSEy2PmZH9mxPAGhKwZGdV19efdpYqLV++Yrz6Qo973ct4vUwwa9YsJYWfQQ6vk4XZ0IkpnkAtI7hpGotISeOi4m+dUwGInLXrFmG1vyl/sOT/p/HWD9mshUp8v8/aWXlgJRDW3gV6g87OMyeVxBoPr1AMw7nLos/M1Az/jcBpptaRnBsaltEzgwConwTzajxATxPvlGclUKnot9b+0QxPjrHncI1DxHGJTl/MUW1LuBS+CoMqHBjI5fqqfVYGZwoagdI4Ajlxtto+ZcrOYr0SfpWX37TxkljYE3E0nzOsaOAANFggAMJUBN/0AJwepsyEWiB0gBp9fXnSvpaHAWA56vRTX9dnHCYRUtZLHJR2kZBVGw/mciO781JrEHtghECpL7Nx8bx58yrq6+vdZctsRkSpv8/L4vZ4AEtfovvss89qFHb3COqMHRYbIqEEaSXqAVjHkEqHAdZRgDXaT3pvby8VIAAmsGU414NE62bwrYdroR1sNG2hHDpCfCY4mnscicUHMgjYhlGDWcIOFbP7+s6adtttt2U3bbK5ZOHy5ZrwURbvpS44aYmd9OipPyjG5vACH2Nyn1aQDdMcSEKVSKC4vENAIwRmAdkAIFPacQLQSf3lODfEsQFJtjB0D1HvMEu/sNH8LF1W6sUCi/mEG4aHu6/JZDQNfxhHt9KkUpnEW95S8c/J5Ps3jI7mto6MjPxk/fr1/YUX0kdTDu/3Av7rhcejWOagCn+iaHoPX6P02DHoKHQE8BwBNIcBTwcq7yCLG+xH/T0H7Wf/IMcOoQYPswbRMVbrYFnjoxtzufu4qZbCK3UbS+CQrxLOnz9/WmNj0jQ1jbq1te1EVHY7bW2dTCFLv7GhoenG5ubmO0hr3nnNNdf8HfVxg024evXq8eojHnfqyniqEfvVsojHJry7txFE92cyw7ONd5kCNSK10mxr2aagSkicdQAiOsSwMprWc/E7QeeuIPj0bR0df78OMx9g4T6WbKH5oZkxY0Y9oPlmY2Pjiurq6iCZrOQvZmjVv0pIi4fEAhZr4zVD1VdExnR3dz944MCBK++5554ugWvt2pLLN3tJpo+XKiw2grUT3E/CuM14hhUHABfukSdkjMJMTHRFEpXZwBJqhLUglvGwyOGTDztQgwBscFYQfFfSak1pqwkHQPAZOM7ChQt/XldXdz4eYA4nJqalJCsrU6z0l2KwQeCqUNhBK5bgAI+GGZYcbGpqughe3Y299Q7Y0AMdV4vUc1CXVpJxnki+AqylVcZTYunNZZT6zMH7EHOjvidAIblyGODeFLiT5ncNJNurAhKwaKCCqU47Eqs+kfA/nkh4e/r7f/zOw4eXb6DKJaU7EG3f9bTTTvsS0urG2tra0VQqlUBisS5pNWqwhoGFWrtfxVAW4+tYCVqJZoS/nDGIYzyAuenHWapp/fe///0VBYPerFu37gVA4joXoDkrVpTOn1IZb4klLHnMwbuD1fiG+US/nGUda4EG75BlZM1xo17AAlTMmnCcbo63MFj93poarxZwsUrHVXcxOfSSjo4vlyi4JE18VN/pSKhPABoWsXXjAg+SSaGXguGeOg4upBl+Cuoem1J1AEu8r68vB/iWX3311X921113/UqSCTI333xz09y5cxcJoEePHt3OsS6ep2sY3LALhOnnhJbxBpZe1oKLOXjrYPyvYehnOXbVAAt+DsE0vLzC2nuWUTZwWpfJOBeyuP9htlWojZrq6twsY770jalTd1zS3n5fCdpafBMajUrcAqmzmcDK4jcUgUYAE5DyXmHKSi/UoZVYAp0AhjoUyByOh9DN3OP+66+//j+cfvrpf40UfDMfZprYl8Zaj73nPe95kClnX1UdwMVm4lXjRABL/LXgOnbsWDv7H5k+ffrfIfbPDVl+WBH5QQjmBk48vuFIZ6fXwwKgDyWT6djAQOjx57xmILkaamvD08LwrtVTppy7vLPzaa5Vz5WChygA6f2q6GDZRyp/4NlxXIL4OMgEJgENUNggsPZ1THYX11LVWYrUuhOv8n2LFy828MvgPQZSpwCrkTpX8vvKJ5988svUvZn7ENmZ2HW2JgpYYraYr05w8X4Osv1JgdicKKrQ73krf5fJrK+NxXKVAwOxGOBqTSb9mel0Yy9/Qcnp7Lxw3fLlToidQX0wNqFFTQ6xI2ezbaSTVRgf1NJFWh5SSxzlScfGklpdPKe6FOYjBrquZebMme9DSoVI+QDD3sUZKIJVFRnQrgwXLFhw444dO3oA1VqumVC1WGycXmIiipgigKkdMnYF9LHkUSFOEtOPj+ZyX/7tyEhs59BQbjeR+U6WzctWVeXm1tef98158765Yv16f+My/VWwCS8CFn3rzIAqAA5/PycPIA1JiRg3tapOhvoQsToZ60WSjaXzItXVtQIXv31CEA4fIUKbP8JIEfggG6Kgjt49R2jj0xs2bDidZzPerxUCJ6ZM2IOf97pSYQIYI81/QPYY4JLLdNOzo6P3/W54OPbE0JC/F3B18wd1CF3n5qfTK78wb96HL2FIZEN+vO15tx+3nw7jfvowQsIGXer4IpiKgOK49fwEIDw/Rqz6WGC3R3ErS70MYwEeg81kwafrVDeZTHryInX897//vQWb3goAWSracNSJp9PpTxXe2AKwsD+um1IB1ku9dFGqOWEQvPeJkZF9W4aG3CcHBoJn+NL7XddL1tT4M+rqvrZ64cJzBK5w/KPVTiEcED711FOkY5g4IYZ3qtOL0qe4lVFeDCkIVAIUtqY5cuSIBdj27duNSBJM1wMSg/qz9piAKIMfr9reV+dVjz+kYO/JM9SfpLa5FxNDU9TGBl1firmv17kJQ/SreCGJekmw86vj8c0X1tYGS9Np94zmZqeNfRZid/d2dz+97dChN//jwYNd9IqNgb2K57ySS5SR4G7KDx4jPFMtBEOvpPM/2dDQsPiZZ54JOzs7HXlvnJPUsaS41djf8ghlsAt8Dz74oD0nMGFHWWAprCAQSSXKqNd9eI6VaLpW91ddjP2Q+zqAu2fVqlXTH3vsMcWcda0+znEtE2m8v9IXFajU3t8yGH39tsHBbwGwXII/peqxwmwTKnF6Q8Ns4mH8kUnnstVSiYUOf6UPehn1XSSUot8+oApmz57dgmS5ng7/GNsWdTwdzDCN50p1SUIVY1hIk+O3l/0kMOmcrjl48KDB+DbTpk2z6k8qUKpTdbivBY+k1u7du+UNmjP5a6A8m1yuKVaSSfVSN8DQr/jCF75wFnz4TcHOioB1nOsn35ENpgXTb+3MZN60pb9/JdPFcjHA5TLtpQZj/rSmpkvXnn32mtWbNq25denS+EcffVQDwKeqeADKCFCQ3P7FuPkfAhTXCFCSREgrFtPlj+jgkAgM8uKkqgQygUGSR0UgELAEHEkhgUgAE4ieeOIJCxZUJRlDg47OydYSQGV/CVCEHiwAeT7DQlX2nvxHyrNnAFYVbfrZ3XfffR7H9gpcPHdcQzHlpAqPM4+doi2x+cza2vPe3NjoL2xu9mY1NYWVRLmP9fXFdnV1Lf/Ctm0/Xgm4uufMCaSSLr744oC/wePQMQ5bohN2wusf+5plP7ljh1LmzJmzGMDcTIf+Z9RQQh2LSsrRqVLXjsCCtxcAiBDbyRMY6Fifep4kk+ykovoTqAQogU5bGfKSWkTUFUJwCS0EhC4EYlfnCI6aiy66yOAkWHXIoH6RL2O3ku4eAL33jDPOuBxg8fgIWGMZ9GL7FlicnF0Ri/3ujfX1DUuamsyClhanrb4eV9GEB3t6en9/9Oiy78kSfokio/tFAGYBJelUvJyo92UEJW9Azb1DgCoAxAIK+8cCCgnjQ64kjYAAPco5DTKfzTYHoDwApoi6tasKgVBrnOs51LPOCio0hq10DHA2ysZC+oU4BPo4zKJFiyyoaA8pOE3F5h3fShoKSIAzxGZb9K53vYs/LjW+UqucbKzjjGNHYl3SYV8ml3vvE31996aUOUCneRUVTlM6HTQ1NNTP9P0733XeeSu6U6nzalOpuTD7AqRCH4byY1y79bnnntsAcPrZN3hRxdQUB4Pck0FeAJVHtPvtdO4tAOoyGckCFPfgr3+R6xIEMaRTCAlQTn9/vydPj7DB79n/Kse/y0SSWlScbL/L5RWi1nKSaopHjZVYnJeB7ume0CNI1iu3bt16Kff9G+ovQjopJsafAM7HwrR9iRLQ1hjt/gvqPLlx40Z9jOOmDssVWOJn0Zi/r290dO3jPT2rWe8hi72lEV6PzzyUZ5YbHNzRjPdYjVcmVVOgd8u+ASgHsFW+s2vXri+S79S3cuXK+JgU4QSA+gAd81+RUmcphoSUYW2rCrnziqR7dLhsIAHKw6j2ZEsBqHuxjT4HgH5VaKNDWnYP++8AGJ8APJ8BNFNkLwlc0p4iqUIVAEg+o/+t/fv3f5GffdA/QXfR1lup/1dIMB+w6nk2BiaJpXcaWwRQFd2f/dqx58Zr/w9bNF5PPXXPyRvzYbjmyPDw4l3d3e/BS/Qhj8QnZ6s+fYY+klVVPpKBv3tTyQQGm2EgdePSWdNREauXLFlyBbbLdYBqC+qmlaGTKzj3Sb74xWwFKA0EK12FBf9zsp8koZhzMxADRAIUafr9dyOp/hFQ3Dvm9SRV9QHYniYF5mt4cOtow3+kaW+FZgEktYVbOz38/hn3/AVAPVy4h9AWhzJItyc4LvUaSsVKKmI32gFsGfBFYBaus86A4mO0LY/Y4olx2pY7sMQm23EET1fu7+ubhTJcyhLFwQ6yIGKplFuXBwV/u1lJddbIPv6FAxKbVIe99CY69pcXXnjhXXTu1RjmzQVACZAChUtdT8a4JAYgEqBidJwA9UM6WdJlhxpTMJTVmQFUtM+KDoIHGASa7xaIzUmLAKnrRfp45O0RFx4ygFMGvQWU3kfSTp6kYlmA38a6aKOCro5CHaj7Z7h+3MtkAJY6Te/RhY/9f0dratZ1EUOqwduqUpYAdoxsIgEFtWbjRXz9ltGoCrnyccChTIEm1OINSAXVKYYMPHUadQKOh0gTT4BCChzmmn9CctzBjSyg2AoM1vtjUwSUDo0tRemlump3EXCqIwCLVGfs9QKX7LqHacMRgNUIiKzHKCklla42K7Sh91ShvWqrg4ofYfjn5zqG0T/2njr0upbJACwxSMz3+JJXNhEsTKDu4gXbRcaxGC43X3aSUk30GwllI9cARNcrLTjknAxydTp/2SsfMhCgIK+gfg4DrG9xzTepU1RXRZC83I4TmF7S6ub82KL63r59+3pR0V8BNJ/FWxxFncYBlSNbTeOLei9JLBWka5Z2JnCIb7/jjjsOESoZ90yHyQAsdayPnXEumLgUXRSQWO9B1u4QTiShBCZJLXWAQCZgqVPAk1Ul7DMzORHjSzeASDaUBondgoe3A3B9lQ78V55VBJR4J0C/XEBR9VUXPUN/oOD/ENOaRTuuIw4X0EYfADmoPZf3so6FVDHgSwC+h++///6/kbdLyrKVeq/66a/iwskALKkPlasFIgAjJlqDVeCRutBWRrtIINOXrWMqcv+LLr+uF2GbqLOKIYOv0IE/oOqovYAOZqtnvBKpU7j0NW0sOBiAXonkGgRQN6Lq9KdWrH0FsDy9H6CSt3jX5s2bV/G0AbxdvehYlfuaGvFyL54MwBLDwYmj4QtrPGOAWyO2uJVUki1S/K19AUjHVQQyEb8DBS8pO/bu3fsp6m3gdBFA4pUkx3hIKB7zgqLGCiQOkusmxhPvQl3/F2yuBXwI8/k49nBuL5L1R9iARc90QkCllpc7sMQ4rWSThskLBBhIQUsLJIUGZHzLXhJJ9UnFqQhYxWPaqi7XKQDpMnzyNF/8vQBNY2wTJaFsO5/3X1HyeAxYP8Q5kcYja5Bc/cW6hWCv6hbrF0+N23ZCYhyn+u0EBkCQEDiKJEAVQSVA8SXzx8X7raFbMMTtvo4JbFKJqq/rqd9EG701a9aoqZJQE9ZBasBJirW5OC7Qu6htgUoZF/rtFSa3Tmiby11iiXkOU6H6mEiwF4C9CXBoySRnLKhke0hCqei4bCz9LoDIAqsALmV+hhx/iqo+BnL+Intlyf03ViVbyc0QVMk0styBJUbaRDuGSx5Hap2N5JF3xJLXeWNdxnvRgJeKLBrrOqbfApqOyeiVZGNfLvwW3VgZEdqWQZlQ6XQy/kwGYBU7/5cAi5VcLDCsRCoCiuNWOkkqKdRQ9Ap1HClXVJmKwruEGEZQj/8mZl1Mmg2D0SfjW3Tsj3BgMgBLKsElGv1T4jv7MLbnACCbvcm+lUoa7lAek2JZyo1SKEFhBwGraOgjuXxUoIZqfvnwww/vlr2CrTJW3fwRVkanx3JgMgCLddJt3vkAL7aKjIZ/BTSuVJu8PYFJqbvEfmy8RyBTXjkgtCqyAKwAUElaDWHk38J9lAxYcuplbMeV+n65e4VSg6Fyp5h+/sEf/vCH115++eWMRwd2SEben8YHizlUUo0CmrYCljxFwORDLCDYqyj7x7dt27ZHGaMFz6rU+69k21fOEsuCCsnSQKruneecc87lmrnC1mzZssWmlUjVEeM5vvCGQCVAoQqVQ2WNdWws1tO1U7JW7tmz5w6kXwzvqhgULdmOK/WGlavEsmtPwdw6ZqncTaru5ai7LJ6hf9999xlCDzaEoJBCMcVE44QCljw/gonKUlAelQfIHoIuBVTf5n42c7TUO60c2leWEqugqnwmCtzElPI3A5pRhmISAo9iOTLaJZFkT5G4ZweedU7HAVSOCQ5bAN3jSKufElL4BR2VFaigyFg/RagtR2DZ+XyEDVqRQKsgZXbataeU+KbsBalAxbEkrTDmrSoEWPL6tO7BvwOmtwtkqlcoEaiKnDhF23JUhbbNxKKWAJ4mgQPJw7ohvlV/n/3sZ+0yP4DouH0laQXolKsuMG2Gd84NN9yg5CUBSrZaJKlgwqksZSmxxABU2TnEnkJsJs3fcxX8xGYyb3jDG8znP/95G1K44IIL7KROeYZKjRHYuG6YxDcXFel/7WtfiwB1KtE05l7lCKzjzZcdpZCCAKU4FaEC86Mf/ciQ7WCuvfZaG79SqEFqUfEq1UFqteD5aYmfsn7340wo0Z1yZK41jADLI0gpR5ML8Abt0j5a+2Dp0qWaYm5jVwo1KMJeKDHZY+eff/71DC5vx8b6NuAa96nnxcZM9m052lg2kxIbaxuqrUszUZ5++unwoYceshJKUXbZVDLax4DK9iNSy6YfI7lupf4CwKVlfsqRByWPy3JkqiSWxyTQQwDji6hDl5koWpXYeoTaClAaaH5+4ZwM9Rzq0cF7/Az7+i2KyinmQDkCSyyQ1HKRSp9Dav2WMUGWJU0E8gxFsqeQRC/GKvuXHwDZFbfccgt/SdguqRiB68W49SqPlyuwhJoQW0nToP4cEO0BXMqjkodog6PyEl+kaM1Oxb6qWbHlTNUhqFqufHiRV5z4w+XMUKsSSYPpQko9Qw6VFuTQYPLx6ecavjlZYdDZrnugJL+ovD4cKEevcCwnBC5JIE1311pUdhhHcSs8QJvFoLWllNwn20tqUuAjPKElwLsZrN6pmzFEZB2CsTeO9l8bB8odWJK4Wp/ql0ihS1inQLOZrfEuIEkiMTZoh3QUy1KAlHoBpBVgHr799tsHJ2KW8GvrsvK4utyBpci5JNbtgOZTzLdrYuhGHmJxVRirFhV+UOQdiaX045B6Wu3uc1wbYF+V8oSJ8kDRSVo5GbwhO4DMBM6rANh6AYgAaQ7yFNNS5L0ALJzFIESCxVgc9tY777zzY2MWWzsJa6JDr4UDkwFYen8LLgKfHwFcX0I11mr1FQFLqTPKepBqlNHO8om3PvDAA0VQWe/ytTAwunbyc8CqNMC1AIP92ywS8ixDPT2MG/Yz77CT0MJP+X2Z2CBJxWayfFSTv2dL4A2P20uMGSYB1xSk1VTaVVdsGxItAlSRGdH2FXFA0ug4wIpXFqTUC44Xz0fbU8uByfz1Pv/dZE9FJeJAxIGIAxEHIg5EHIg4EHEg4kDEgYgDEQciDkQciDgQcSDiQMSBiAMRByIORByIOBBxIOJAxIGIAxEHIg5EHIg4EHEg4kDEgYgDEQciDkQciDgQcSDiQMSBiAMRByIORByIOBBxIOJAxIGIAxEHIg78KXPg/wP8f+1DKa05nwAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "final_sys = ChemicalSystem(\n", " \"\"\"\n", "System\n", " Atoms\n", " C 0.8176518265162523 0.006185364744877703 0.11652129166388352\n", " C -0.5531099061901967 -0.6363379504647356 -0.08529176454211394\n", " O -1.4164729952124286 0.48607450912407413 -0.19677333330677277\n", " N -0.6859030532374636 1.632537947036214 -0.08098563281192021\n", " N 0.5074709667559635 1.4283821326486141 0.08689807666792455\n", " H 1.2798721321383 -0.22791534308315498 1.0777663111036928\n", " H 1.5305100901334066 -0.2099050616627963 -0.6819652771198405\n", " H -0.866292408825807 -1.2485848049850283 0.7685478972017538\n", " H -0.6137266520780276 -1.2304367933580558 -1.0047175688566132\n", " End\n", "End\"\"\"\n", ")\n", "\n", "view(final_sys, width=150, height=150, direction=\"tilt_pca3\", guess_bonds=True)" ] }, { "cell_type": "markdown", "id": "3550a296-b915-41bd-8c35-a98976ce7ce3", "metadata": {}, "source": [ "## NEB and DFTB settings. Run the Job." ] }, { "cell_type": "code", "execution_count": 3, "id": "86411f51-575d-4595-b2f0-218eb448de12", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Properties\n", " PESPointCharacter Yes\n", "End\n", "\n", "Task NEB\n", "\n", "System\n", " Atoms\n", " C 1.0065881512963437 -1.417955004240563 0.12885420257373528\n", " C -0.08873018725926754 -2.1350000996605933 -0.03448649693154524\n", " O -1.5081108856471217 1.002080375323799 -0.2045570136880532\n", " N -0.4803377947492119 1.514182931120784 -0.052917508314934095\n", " N 0.4992559331764258 2.001851603909277 0.09160944741412266\n", " H 1.3565548327663348 -1.1143852514314225 1.1071131753303691\n", " H 1.6148659535066119 -1.095826278535031 -0.7065218137603257\n", " H -0.6985763255773102 -2.454234856430313 0.8007371044735704\n", " H -0.4402537210929026 -2.435675244924285 -1.0129788309893883\n", " End\n", "End\n", "System final\n", " Atoms\n", " C 0.8176518265162523 0.006185364744877703 0.11652129166388352\n", " C -0.5531099061901967 -0.6363379504647356 -0.08529176454211394\n", " O -1.4164729952124286 0.48607450912407413 -0.19677333330677277\n", " N -0.6859030532374636 1.632537947036214 -0.08098563281192021\n", " N 0.5074709667559635 1.4283821326486141 0.08689807666792455\n", " H 1.2798721321383 -0.22791534308315498 1.0777663111036928\n", " H 1.5305100901334066 -0.2099050616627963 -0.6819652771198405\n", " H -0.866292408825807 -1.2485848049850283 0.7685478972017538\n", " H -0.6137266520780276 -1.2304367933580558 -1.0047175688566132\n", " End\n", "End\n", "\n", "Engine DFTB\n", " Model GFN1-xTB\n", "EndEngine\n", "\n", "\n" ] } ], "source": [ "from scm.plams import Settings, AMSJob\n", "\n", "settings = Settings()\n", "\n", "# Input options for the AMS driver\n", "settings.input.ams.Task = \"NEB\"\n", "settings.input.ams.Properties.PESPointCharacter = \"Yes\"\n", "\n", "# Input options for the engine (DFTB in this case)\n", "settings.input.DFTB.Model = \"GFN1-xTB\"\n", "\n", "# You can pass multiple systems to an AMSJob by using a dictionary of molecules.\n", "# The key of the dictionary will be used as the header of the 'System' block\n", "molecules = {\"\": initial_sys, \"final\": final_sys}\n", "\n", "# Create and run the job\n", "job = AMSJob(settings=settings, molecule=molecules, name=\"NEB\")\n", "print(job.get_input())" ] }, { "cell_type": "code", "execution_count": 4, "id": "6b440e6c-7270-4081-a41d-132f8727737a", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[28.02|09:45:50] JOB NEB STARTED\n", "[28.02|09:45:50] JOB NEB RUNNING\n", "[28.02|09:45:51] JOB NEB FINISHED\n", "[28.02|09:45:51] JOB NEB SUCCESSFUL\n" ] } ], "source": [ "job.run();" ] }, { "cell_type": "code", "execution_count": 5, "id": "ac5d94e9-651f-4188-800e-9cf1ce357948", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Successful NEB calculation!\n" ] } ], "source": [ "# job = AMSJob.load_external(\"/path/to/ams.rkf\") # if loading results from disk\n", "assert job.ok(), \"Looks like NEB calculation failed?\"\n", "print(\"Successful NEB calculation!\")" ] }, { "cell_type": "markdown", "id": "dabbde26-8922-4673-a104-c84caab3f814", "metadata": {}, "source": [ "## NEB Results" ] }, { "cell_type": "code", "execution_count": 6, "id": "8d132792-ea30-45c5-a55f-04c0b26ef03f", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "['NEB_TS_final']\n" ] } ], "source": [ "print(job.results.engine_names())" ] }, { "cell_type": "code", "execution_count": 7, "id": "55175a0a-2584-4c57-90ba-2c49f928676f", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "pes_point_character='transition state'\n" ] } ], "source": [ "from pprint import pprint\n", "\n", "pes_point_character = job.results.readrkf(\"AMSResults\", \"PESPointCharacter\", file=\"NEB_TS_final\")\n", "print(f\"{pes_point_character=}\")" ] }, { "cell_type": "code", "execution_count": 8, "id": "4d6c6ffd-e1d6-4c05-832f-aa77e5cc4cb9", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "{'Climbing': True,\n", " 'Energies': [-17.333016762960952,\n", " -17.33294947958635,\n", " -17.33256809516803,\n", " -17.331527201901654,\n", " -17.328829238015103,\n", " -17.319753903625962,\n", " -17.307297193368917,\n", " -17.322575012605878,\n", " -17.364900430989824,\n", " -17.393606196413522],\n", " 'HighestIndex': 6,\n", " 'HistoryIndices': [309, 301, 302, 303, 304, 305, 306, 307, 308, 310],\n", " 'LeftBarrier': 0.025719569592034475,\n", " 'Molecules': [Molecule('C2H4N2O' at 0x7e3cea270370),\n", " Molecule('C2H4N2O' at 0x7e3cea270670),\n", " Molecule('C2H4N2O' at 0x7e3cea270790),\n", " Molecule('C2H4N2O' at 0x7e3cea270280),\n", " Molecule('C2H4N2O' at 0x7e3cea270b20),\n", " Molecule('C2H4N2O' at 0x7e3cea270dc0),\n", " Molecule('C2H4N2O' at 0x7e3cea6c9160),\n", " Molecule('C2H4N2O' at 0x7e3cf20742b0),\n", " Molecule('C2H4N2O' at 0x7e3ce9ffb0a0),\n", " Molecule('C2H4N2O' at 0x7e3ce9ffb310)],\n", " 'ReactionEnergy': -0.06058943345257006,\n", " 'RightBarrier': 0.08630900304460454,\n", " 'nImages': 8,\n", " 'nIterations': 216}\n" ] } ], "source": [ "neb_res = job.results.get_neb_results()\n", "pprint(neb_res) # energies in hartree" ] }, { "cell_type": "code", "execution_count": 9, "id": "120d755f-dcea-4170-8e39-0eb74ee34e72", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZcAAAEGCAYAAACpXNjrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuNSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/xnp5ZAAAACXBIWXMAAAsTAAALEwEAmpwYAAAszklEQVR4nO3deXxU9b3/8dcnewIhbGEJyL5viRrFpVJ36wa1rfveWttbbW1re++1VrvY9tfeVu2mrUptXahe26sF0VZxr1pF1BA2QQTZwo4hQCDbfH5/zAmMNAkDmcmZJO/nw/PgzPecM+edAeeTs32/5u6IiIgkUlrYAUREpONRcRERkYRTcRERkYRTcRERkYRTcRERkYTLCDtAKujdu7cPGTIk7BgiIu3K22+/vcXdC5tapuICDBkyhHnz5oUdQ0SkXTGzVc0t02kxERFJOBUXERFJOBUXERFJOBUXERFJOBUXERFJOBUXERFJOBUXERFJOBUXEUm6tR9V879vraYhoiE+Ogs9RCkiSTXvw21c+9DbbNtVy9urPuKnn5lEWpqFHUuSTEcuIpI0T7y7lkvue5OC3EyuOHYwj81by62zFqJBCjs+HbmISMJFIs6dzy3jNy8s59hhvfjdZUdQkJtJbmY697yygpyMdG4+eyxmOoLpqFRcRCShdtc28K2/zOepBeu5sPQwbvv0BLIyoidJ/vvMMdTUR5j+6kpyMtP51hmjQ04ryaLiIiIJs6lqD198cB7l67Zz81ljueaEoR87OjEzbj1nHDX1Dfz2xeVkZ6Tx1VNGhphYkkXFRUQSYlHFdq55YB7bd9dx7+WlnDaub5PrpaUZP/70RGrqItw+ZxnZmWlcO2V4G6eVZFNxEZFWm7N4Izc8+i4FuZn85cvHMr6ooMX109KM//ncJGrqI/zk6ffIzkjnyuOGtE1YaRMqLiJyyNyd6f9cyU/+voRJAwq474pS+nTLiWvbjPQ0fnlRCbUNEb43axHZGWlcdPSgJCeWtqJbkUXkkNTWR7jp8QX8+OklnDmhH49ee2zchaVRZnoav73kcD45qpCbnljAE++uTVJaaWsqLiJy0Cqra7ny/rk8+tYavnryCH578RHkZqUf0ntlZ6Rzz+VHcszQXtz42HyeKl+f4LQSBhUXETkoKzbv5Ly7X+ftVR9x54XF3Hj66FY/cZ+Tmc4frirliEE9uOHRd5mzeGOC0kpYVFxEJG6vf7CF8+5+ne276/jzFydz3uEDE/beeVkZ/PHqoxg/oIDrZrzDy8s2J+y9pe2puIhIXB6du5or/jCXPvnZzLzueEqH9Ez4PvJzMnnw6qMZ0acr1z44j9c/2JLwfUjbUHERkRY1RJyfPL2E/358AceN6M3/feU4DuuZl7T9FeRl8tAXjmZwrzyueWAe8z7clrR9SfKouIhIs3bV1POlh97m3ldWcOWxg7n/ylK65WQmfb+9umbz8DWT6dcth6v++Bbz11QmfZ+SWCouItKkisrdfO73/+KF9zbyg6nj+cG0CWSkt91XRp/8HGZ8cTI9umRyxf1zWVxR1Wb7ltZTcRGRfzN/TSXT7nqNtduquf+qo0J7er5/QS5/vuYYumSlc9kf3uT9jTtCySEHT8VFRD7mqfL1XHDPv8jOSOP/vnIcJ47uE2qew3rmMeOLx5CRZlwy/U1WbtkVah6Jj4qLiADRrlx+8/z7XPfnd5g4oICZ1x3PqL75YccCYGjvLsy4ZjKRiHPJfW+wZlt12JHkAFRcRISa+ga++dh8bp+zjPMOH8CML06mV9fssGN9zMi++Tz0hclU1zZw8X1vUFG5O+xI0gIVF5FObsvOGi65702eeHcd3zp9FHdcUEx2xqF15ZJs44q68dAXjmZ7dR2XTn+TTVV7wo4kzVBxEenElm3cwafveo2F67Zz1yVHcP3JI1N+6OFJA7vzp88fxcaqPVw6/U227qwJO5I0QcVFpJN6edlmPnv369TUR3jsS8dy9qT+YUeK25GDe3L/VUex5qNqLvvDXCqra8OOJPtRcRHphB54/UOu/uNcBvbMY+Z1x1N8WPewIx20Y4b14t7LS/lg006uvH8uVXvqwo4kMVRcRDqR+oYIt85cyPdmLeLkMX3565ePpah7btixDtmUUYXcfekRLKqo4uo/vsWumvqwI0lAxUWkk6jaU8fnH5jHg/9axbVThnHP5UfSJbv9D0Z76ri+/ObiwylbU8k1D8xjT11D2JGEkIqLmZ1vZovMLGJmpTHtl5pZWcwUMbOSJra/zczKg3WeNbOimO3LzWyBmb1uZsVt+GOJpKzVW6v57N2v8/ryLfzssxP5zlljSW/lGCyp5MyJ/bnjgmLeWLmVax96m5p6FZiwhXXkshD4DPBKbKO7z3D3EncvAS4HVrp7WRPb/9zdJwXrzQZuDdpXAp9094nAbcC9yYkv0n689eE2Pn33a2zaUcODXziaC4/qmOPUTysZwM8+M4lXlm3muhnvUtcQCTtSpxbKMbG7LwEOdMvjxcCjzWwf24NdF8CD9tdj2t8AEjeSkUg7U1PfwD0vr+C3LyxnQI9c7r/qKIb27hJ2rKS64KjDqKlv4JaZi/j6o2X86qKSNu1sU/ZJ5ROuFwLTmltoZj8GrgC2Ayc1scoXgL+3sP21wLUAgwZ1zN/kpPN6bfkWbvnbQlZs2cU5k/rzo09PoHteVtix2sTlxw6hpj7Cj55aQlZGGr84v7hDnQJsL5JWXMzsOaBfE4tudveZB9h2MlDt7gubW8fdbwZuNrObgOuB78VsfxLR4vKJFra/l+C0WWlpqbeUR6S92LRjDz9+agkzyyoY0iuPBz9/NFNGFYYdq81dc8Iwauoj/PyZpWRnpPGT8yaSpgLTppJWXNz91FZsfhHwSJzrzgCeJiguZjYJmA6c6e5bW5FBpN1oiDgz3lzFz59ZSk1dhBtOGcl/nDicnMzU7MalLVx30ghq6hr49QvLyc5I4/tTx6d87wMdScqdFjOzNOAC4IQW1hnp7u8HL6cB7wXtg4DHgcvdfVmys4qkggVrt3Pz3xZQvnY7nxjRmx9OG8+wwq5hx0oJ3zhtFHvqI9z7ygqyM9O56cwxKjBtJJTiYmbnAb8BCoGnzKzM3c8IFk8B1rj7iv22mQ783t3nAT81s9FABFgFfDlY7VagF3B38A+o3t1LEemAqvbUcfszS3nojVX06prNry8+nHMn9deXZwwz46Yzx1BT18C9r6wgJyONb54+OuxYnYK563JDaWmpz5s3L+wYInFxd54sX89tsxezdWcNlx8zmBvPGN0mY9u3V5GI850nFvDoW2v4r0+N4T9OHB52pA7BzN5u7hf4lDstJiLNW7llF7f8bSGvLt/CpIEF3H/lUUwcWBB2rJSXlmb8+LyJVFbXceecZVwyeRAFuSrGyaQbwEXagT11Ddw5Zxln/PIV5q+p5IfTxvPEV45XYTkI6WnGl08cTm1DhGcWbgg7ToenIxeRFPfKss3cOnMhH26tZlpJETefPZY++Tlhx2qXigcWMLhXHrPmV3DBUYeFHadDU3ERSVEbq/Zw2+zFzC5fz7BgDPnjR/QOO1a7ZmZMLS7irheXs2nHHhXpJNJpMZEU0xBx/vTaSk65/WWeXbyRb542ir9//QQVlgSZWlxExOGp8vVhR+nQdOQikkLmr6nk5r8tYOG6KqaMKuSHU8czpIP3B9bWRvbNZ0y/fGbNr+Dq44eGHafDUnERSQHbd9fxi2eW8vCbq+iTn81dlxzBWRP76ZmVJJlWMoCf/eM9Vm+tZlCvvLDjdEg6LSYSInfnb++u45TbX2bGm6u46rghPPfNT3K2HoZMqnOL+wPwZHlFyEk6Lh25iITkg807ueVvC3n9g60UH9adP119FBMG6NbitjCwRx5HDu7BrLIKrjtpRNhxOiQVF5E2tqeugbtfXM7vX15BTmYaP/r0BC4+epC6hW9j00qKuHXmIt7bUMWYft3CjtPh6LSYSBt6aekmTr/zFX79wnLOntSf5288kcuOGazCEoKzJvYnPc2YVaZTY8mg4iLSBjZs38N1M97hqj++RWa68ecvTubOC0sozM8OO1qn1btrNscN78WT5RWoj8XE02kxkSSqb4jw4L9WcfuzS6mPON8+YzRfPGEYWRn6vS4VTCsZwLf+Mp9311RyxKAeYcfpUFRcRJKgriHCO6s+4gdPLmbx+ipOGl3ID6ZO0G2vKeaM8X35zhNpzCqrUHFJMBUXkVb6aFctS9ZXsXh9FUvW7+C9DVW8v3EntQ0R+nXL4feXHcEZ4/XMSirKz8nk5NF9mF2+nu+ePZaMdB1RJoqKi0ic6hsifLh1F4vX72DJ+qq908aqmr3rFOZnM6ZfPlcfP4Sx/btx6ri+dM3W/2apbFpJEf9YtIE3VmzjEyPVxU6i6F+9SBO2V9exeH0V721oLCI7WLZxBzX1EQAy0owRfbpy3PDejO2fz9j+3Rjbvxu9u+oCfXtz0pg+dM3OYGbZOhWXBFJxkU6tIeJ8uHXX3qOQ94Kjkorte/au06tLFmP7d+OKYwczpl+0iIzo01UX5TuInMx0Th/fl38s2sCPzptAdkZ62JE6BBUX6TS2765j6YaPn9JaunEHe+qiRyPpacbwwi4cNbQnY/t3Y0y/fMb170Zhfraul3RwU4uLePyddby0dDNnjO8XdpwOQcVFOpxIxFm1rZr3ggLSeI1kXeXuvet0z8tkbL9uXHL04L2ntUb06UpOpn5r7YyOH9GbXl2ymDW/QsUlQVRcJGU1RJyq3XVU7q6jsrqWyt11bK/eN19ZXcf2/ZZ9VF3L9t11RIJn4tIMhhV25fBB3blk8iDGBddG+nbT0Yjsk5mexlkT+/PYvDXsrKnXTRgJoE9Qkq6uIRIUgaAQVO8rGHvb939dXUvVnvoW3zc/J4MeeVl0z8ukIDeTgT3y6J4bnT+sZy5j+3djVN98HY1IXKaWFPHQG6uYs3gD5x0+MOw47Z6KSyfi7tQ1OLUNEerqI9Q2RKiN+bOu8fXH2pzahoagzfcur4vZpiZYv64+QnVtA5W7gwISHFnsrGm+SKQZFORm0j0vi4LcTHp2yWJY7y57X3fPC6bcLAryMukerNstJ0PPJEhCHTmoBwO65zKrrELFJQFUXFqhfG0lf35zNe7gOO4QCeaD/3D34E+IBPPErB/b7sFCD7aN+L75xr6PGvcVicS8B9HrDB/7oo8tDjEFI5HMICs9jayMNLIz0shMTyM3K50eeVn065bD6H75dM/N2lsgGotI971FI4v87AzS1GmjpIC0NOOc4v784Z8r2barlp5dssKO1K6puLTCpqoaXnhvE2ZgGGnG3vP4ZuxtN4M0MwzAwIiuZwTtwXdrY1vjtvu2iWmP2caCNzOi54y7ZGeQlRH9ss9KT9v7xZ+Z3thme5fva9v35962YHl2xv5tRnZ6+t55HTlIRzO1uIh7Xl7B0wvWc9kxg8OO0661WFzMLAc4BzgBKAJ2AwuBp9x9UfLjpbZTx/Xl1HF9w44hIgkyLrhrcNb8ChWXVmr2V08z+wHwGnAs8CZwD/AYUA/81MzmmNmkNkkpItIGzIypxUXMXbmNiphb1+XgtXReY667H+nuN7r7n939OXef7e53uPu5wKWATkqKSIcytbgIgNnlGkSsNZotLu7+VOxrM8vbb/kmd5+XrGAiImEY0rsLkwYWMGu+iktrHPCKrJkdZ2aLgfeC18VmdnfSk4mIhGRqcREL11XxweadYUdpt+K53edO4AxgK4C7zwemJDOUiEiYzi0uwgxmleno5VDFdS+pu6/Zr6khCVlERFJC3245TB7akyfnV+x9xkwOTjzFZY2ZHQe4mWWa2beAJUnOJSISqmklA1ixZReLKqrCjtIuxVNcvgxcBwwA1gElwWsRkQ7rzAn9yEw3ZpatCztKu3TA4uLuW9z9Unfv6+593P0yd9/aFuFERMLSPS+LKSMLmV2+nkhEp8YOVjx3i40ys+fNbGHwepKZfTf50UREwjW1pIj12/fw1ofbwo7S7sRzWuw+4CagDsDdy4GLkhlKRCQVnDauL7mZ6czUMy8HLZ7ikufuc/dra3mgjQMws/PNbJGZRcysNKb9UjMri5kiZlbSxPa3mVl5sM6zZla03/KjzKzezD7Xmpwi0rnlZWVw6ri+/H3BeuoS3Kt4RxdPcdliZsOJ9uxO8IW9vpX7XQh8BnglttHdZ7h7ibuXAJcDK929rIntf+7uk4L1ZgO3Ni4ws3TgZ8CzrcwoIsK04iI+qq7j1fe3hB2lXYmnuFxHtNPKMWa2Dvg60TvIDpm7L3H3pQdY7WLg0Wa2j703sAtB4Qt8Ffg/YFNrMoqIAEwZVUhBbqbuGjtIB+pyPx34irufamZdgDR339E20bgQmNZCth8DVwDbgZOCtgHAecHro1p6czO7FrgWYNCgQYlJLCIdTlZGGmdO6Mes+RXsrm0gN0vDZsejxSMXd28APhHM7zqYwmJmz5nZwiamZgtGzLaTgWp3X9hCtpvd/TBgBnB90PxL4L/c/YAnR939XncvdffSwsLC+H4oEemUphYXUV3bwPPvbQw7SrsRz0iU75rZLOAvwK7GRnd/vKWN3P3UVuS6CHgkznVnAE8D3wNKgUeD0SB7A2eZWb27/60VWUSkk5s8rBd98rOZWVbBOZOKDryBxFVccoh2WnlyTJsDLRaXQ2VmacAFREe/bG6dke7+fvByGkGPze4+NGadPwGzVVhEpLXS04xzJhXx8Bur2L67joLczLAjpbx4LuhPd/erYyfgD63ZqZmdZ2ZriY5y+ZSZPROzeAqwxt1X7LfN9Jjbln8anGIrB04HbmhNHhGRA5laUkRtQ4RnFm4IO0q7YAfq8dPM3nH3Iw7U1p6Vlpb6vHka90xEmufunPiLlxjYI5cZ1xwTdpyUYGZvu3tpU8uaPS1mZscCxwGFZvbNmEXdAN0uISKdipkxtbiIu15czqYde+iTnxN2pJTW0mmxLKAr0QKUHzNVAXryXUQ6nanFRUQcnipv7XPkHV+zRy7u/rKZvQpMcvcftGEmEZGUNLJvPmP7d2PW/AquPn7ogTfoxOJ5zkX33YmIBKYWF/Hu6kpWb60OO0pKi+dusTIzm2Vml5vZZxqnpCcTEUlB5xb3B+DJcvWU3JJ4ikvscy7nBtM5yQwlIpKqBvbI48jBPZhVpuLSkgM+RBk81yIiIoFpJUXcOnMR722oYky/bmHHSUnxjESZY2bXmdndZnZ/49QW4UREUtFZE/uTnmY6emlBPKfFHgL6AWcALwMDgbbqGVlEJOX07prNccN7MWt+BQd6EL2ziqe4jHD3W4Bd7v4AcDYwObmxRERS27SSAaz9aDfvrK4MO0pKiqe41AV/VprZBKAA6JO8SCIiqe+M8X3Jykjjyfk6NdaUeIrLvWbWA7gFmAUsBv4nqalERFJcfk4mJ4/uw+zy9dQ3HHAIqU7ngMXF3ae7+0fu/rK7D3P3Pu7++7YIJyKSyqaVFLFlZw3/WrE17Cgp54C3IptZNvBZYEjs+u7+w+TFEhFJfSeN6UPX7AxmlVVwwkiNaBsrntNiM4kOyFVPdCTKxklEpFPLyUzn9PF9+ceiDdTUN4QdJ6XEMxLlQHf/VNKTiIi0Q9NKBvD4O+t4aelmzhjfL+w4KSOeI5fXzWxi0pOIiLRDxw/vRa8uWXqgcj8tDRa2APBgnavNbAVQAxjg7j6pbSKKiKSujPQ0zprYn8fmrWFnTT1ds+M5IdTxtfQpqHNKEZE4TC0p4qE3VjFn8QbOO3xg2HFSQkvFZau772xpYzPreqB1REQ6uiMH9WBA91xmllWouARauuYy08xuN7MpZtalsdHMhpnZF8zsGUAX+kWk00tLM84p7s+r729h267asOOkhGaLi7ufAjwPfAlYZGbbzWwr8DDRjiyvdPe/tk1MEZHUNrW4iPqI8/SC9WFHSQktXnly96eBp9soi4hIuzWufzdG9OnKrLIKLjtmcNhxQhfPrcgiInIAZsbU4iLmfriNisrdYccJnYqLiEiCTC0uAmB2uZ55UXEREUmQIb27UDywgJl6oDKuYY5vN7PxbRFGRKS9O7e4iEUVVXywuXM/pRHPkcsSomO6vGlmXzazgmSHEhFpr84tLsKMTt8dTLzjuRwPXEG02/1yM/uzmZ2U7HAiIu1N3245TB7ak1nzK3D3sOOEJq5rLmaWDowJpi3AfOCbZvZoErOJiLRL00oGsHLLLhauqwo7SmjiueZyJ7AUOAv4ibsf6e4/c/dzgcOTHVBEpL05c0I/MtONWfPXhR0lNPEcuZQDxe7+JXefu9+yo5OQSUSkXeuel8WUkYU8OX89kUjnPDUWT3GZD4w2syNipuFmluHu25MdUESkPZpaUsSGqj3M/XBb2FFCEc/AA3cDRxA9gjFgArAIKDCz/3D3Z5OYT0SkXTptXF9yM9OZNb+CY4b1CjtOm4vnyKUCONzdS939SKLXWVYApwH/k8xwIiLtVV5WBqeO68vTC9ZTWx8JO06bi6e4jHL3RY0v3H0xMMbdVyQvlohI+zetuIjK6jpeXb457ChtLp7istjMfmdmnwymu4O2bKAuyflERNqtKaMKKcjN7JQPVMZTXK4ElgNfD6YVwFVEC8shPUhpZueb2SIzi5hZaUz7pWZWFjNFzKykie1vM7PyYJ1nzawoZtmJQfsiM3v5UPKJiCRCVkYaZ07ox7OLN7K7tiHsOG2qxeISPDz5tLvf7u7nBdMv3L3a3SOtGOJ4IfAZ4JXYRnef4e4l7l4CXA6sdPeyJrb/ubtPCtabDdwa5O1O9AaEqe4+Hjj/EPOJiCTE1JIiqmsbeG7JxrCjtKkWi4u7NwCRRPcn5u5L3H3pAVa7GGiyBwB3j33stQvQeCP5JcDj7r46WG9Ta7OKiLTG5KG96JOfzaz5nevUWDy3Iu8EFpjZHGBXY6O7fy1pqaIuBKY1t9DMfky0v7Pt7Ds9NwrINLOXgHzgV+7+YDPbXwtcCzBo0KDEpRYRiZGeZpwzqYiH31jF9uo6CvIyw47UJuK55vI4cAvRU1hvx0wtMrPnzGxhE1OzBSNm28lAtbsvbG4dd7/Z3Q8DZgDXB80ZwJHA2cAZwC1mNqqZ7e8Nbq8uLSwsPFAkEZFDNrWkiNqGCP9YtD7sKG3mgEcu7v6AmeUCg+I4lRW73amtyHUR8Eic684Anga+B6wFtrr7LmCXmb0CFAPLWpFFRKRVigcWMLhXHrPmV3DhUZ3jTEk8HVeeC5QB/whel5jZrGQFMrM04AKaud4SrDMy5uU04L1gfibwCTPLMLM8YDLR8WhEREJjZkwtLuJfH2xl0449YcdpE/GcFvs+0Q4qKwGCu7eGtWanZnaema0FjgWeMrNnYhZPAdbs/5CmmU2PuW35p8EptnLgdOCGINsSokWwHJgLTG/p1JqISFuZWlxExOGp8s5xaiyeC/p17r7dzGLbWtWXgbs/ATzRzLKXgGOaaL8mZv6zLbz3z4GftyafiEiijeybz9j+3ZhZVsHVxw8NO07SxXPkssjMLgHSzWykmf0GeD3JuUREOpypxUWUralk9dbqsKMkXTzF5avAeKCG6EX2KqJP6ouIyEE4t7g/AE+Wd/xnXg5YXIKn8W9296OCW3dvdvfOcUVKRCSBBvbIo3RwD2aWdfwRKuO5W2yUmd0b9OH1QuPUFuFERDqaqSVFLNu4k/c2VB145XYsntNifwHeBb4LfDtmEhGRg3TWxP6kpxkzO3hPyfHcLVbv7r9LehIRkU6gd9dsThpdyCNzV/OVE4eTn9Mxu4OJ58jlSTP7ipn1N7OejVPSk4mIdFBfO2UkldV13P/qh2FHSZp4x3P5NtHbjxv7FZuXzFAiIh3ZpIHdOX1cX6b/cwWV1bVhx0mKeO4WG9rE1Kon9EVEOrtvnDaKHTX13PfPjjlifLPFxcz+M2b+/P2W/SSZoUREOrqx/btxzqT+/PG1D9m6sybsOAnX0pHLRTHzN+237FNJyCIi0ql8/dRR7Klr4PcvfxB2lIRrqbhYM/NNvRYRkYM0ok9XPn34AB781yo2VnWsZ9NbKi7ezHxTr0VE5BDccMpIGiLOXS8uDztKQrVUXIrNrMrMdgCTgvnG1xPbKJ+ISIc2uFcXzi89jEfmrmbtRx2nQ8tmi4u7p7t7N3fPd/eMYL7xdcd86kdEJARfPXkEhvHbFzrO0Us8z7mIiEgSFXXP5ZLJg/jL22v5cMuusOMkhIqLiEgK+MqJw8lMN379/PthR0kIFRcRkRTQp1sOVxw7hCfK1vH+xh1hx2k1FRcRkRTxpSnDyMtM55fPtf+jFxUXEZEU0atrNp//xFCeWrCexRXte7wXFRcRkRRyzQnD6JaTwR1zloUdpVVUXEREUkhBbiZfPGEYzy3ZSNmayrDjHDIVFxGRFHP1J4bSIy+zXR+9qLiIiKSYrtkZ/MeJw3ll2WbmrtwWdpxDouIiIpKCLj9mCIX52dz+7FLc2193jiouIiIpKDcrnetOHM6bK7fx+gdbw45z0FRcRERS1MWTB1FUkMMv2uHRi4qLiEiKys5I5/qTR/Lu6kpeXLop7DgHRcVFRCSFnV86kEE987j92WXt6uhFxUVEJIVlpqdxwykjWVRRxTOLNoQdJ24qLiIiKe7Thw9gWGEX7pizjIZI+zh6UXEREUlx6WnGN04dxbKNO5ldXhF2nLiouIiItANnT+zPmH75/PK596lviIQd54BUXERE2oG0NOMbp41i5ZZdPP7uurDjHJCKi4hIO3H6uL5MHFDAr557n9r61D56UXEREWknzIwbTx/FusrdPDZvTdhxWqTiIiLSjnxyVCFHDu7Bb154nz11DWHHaVYoxcXMzjezRWYWMbPSmPZLzawsZoqYWUkT299mZuXBOs+aWVHQXmBmT5rZ/OD9r27DH0tEJOkaj142VtUw483VYcdpVlhHLguBzwCvxDa6+wx3L3H3EuByYKW7lzWx/c/dfVKw3mzg1qD9OmCxuxcDJwK3m1lWUn4CEZGQHDe8N8cN78XvXlpOdW192HGaFEpxcfcl7r70AKtdDDzazPaxg0t3ARqfKnIg38wM6ApsA1LzkxcRaYUbTx/Flp21PPD6qrCjNCmVr7lcCDzS3EIz+7GZrQEuZd+Ry2+BsUAFsAC4wd2bvKXCzK41s3lmNm/z5s2JTS4ikmRHDu7JiaMLueeVD9ixpy7sOP8macXFzJ4zs4VNTNPi2HYyUO3uC5tbx91vdvfDgBnA9UHzGUAZUASUAL81s27NbH+vu5e6e2lhYeHB/XAiIingxtNGU1ldx/2vfhh2lH+TtOLi7qe6+4QmpplxbH4RLRy17GcG8Nlg/mrgcY9aDqwExhx8ehGR1DdxYAFnjO/L9H+uoLK6Nuw4H5Nyp8XMLA24gGautwTrjIx5OQ14L5hfDZwSrNMXGA2sSE5SEZHwfeO0UeysrefeV1Lrqy6sW5HPM7O1wLHAU2b2TMziKcAad1+x3zbTY25b/mlwiq0cOB24IWi/DTjOzBYAzwP/5e5bkvrDiIiEaEy/bpwzqYg/vf4hW3bWhB1nL2tPg88kS2lpqc+bNy/sGCIih+SDzTs57Y6X+fzxQ/nuOePabL9m9ra7lza1LOVOi4mIyMEZXtiV8w4fyENvrGJj1Z6w4wAqLiIiHcINp4ykIeLc9eLysKMAKi4iIh3CoF55XHDUYTwydzVrP6oOO46Ki4hIR3H9SSMwjN88H/7Ri4qLiEgHUdQ9l0smD+Kv76zlwy27Qs2i4iIi0oF85aThZKYbv3r+/VBzqLiIiHQgffJzuPLYIfytbB3vb9wRWg4VFxGRDuZLnxxOXmY6v3wuvKMXFRcRkQ6mZ5csvvCJoTy1YD2LKraHkkHFRUSkA/rCCcPolpPBnXOWhbJ/FRcRkQ6oIDeTa6cM47klmyhbU9nm+1dxERHpoK46fig9u2Rx+7MHGvg38VRcREQ6qK7ZGXz5k8P45/tbmLtyW5vuW8VFRKQDu/yYIRTmZ/OLZ5fSlr3gq7iIiHRguVnpXH/SCOau3MZry7e22X5VXEREOriLjj6MooKcNj16UXEREengsjPS+eopIylbU8kL721qk32quIiIdAKfO3Igg3rmccecZUQiyT96UXEREekEMtPT+PqpI1lUUcUzizYkfX8qLiIincS0kgEML+zCHXOW0ZDkoxcVFxGRTiI9zfjGaaN4f9NOZpdXJHVfKi4iIp3IWRP6M6ZfPnfOWUZ9QyRp+1FxERHpRNLSjG+eNooPt1bz+DvrkrefpL2ziIikpNPG9WXSwAJ+9fz71NYn5+hFxUVEpJMxM248fTTrKnfzv/PWJGUfKi4iIp3QlJG9mVpcRPfczKS8f0ZS3lVERFKamfHriw9P2vvryEVERBJOxUVERBJOxUVERBJOxUVERBJOxUVERBJOxUVERBJOxUVERBJOxUVERBLO2mo85VRmZpuBVYe4eW9gSwLjtHf6PD5On8c++iw+riN8HoPdvbCpBSourWRm89y9NOwcqUKfx8fp89hHn8XHdfTPQ6fFREQk4VRcREQk4VRcWu/esAOkGH0eH6fPYx99Fh/XoT8PXXMREZGE05GLiIgknIqLiIgknIpLK5jZp8xsqZktN7P/DjtPmMzsMDN70cwWm9kiM7sh7ExhM7N0M3vXzGaHnSVsZtbdzP5qZu+Z2RIzOzbsTGExs28E/48sNLNHzCwn7EzJoOJyiMwsHbgLOBMYB1xsZuPCTRWqeuBGdx8HHANc18k/D4AbgCVhh0gRvwL+4e5jgGI66ediZgOArwGl7j4BSAcuCjdVcqi4HLqjgeXuvsLda4FHgWkhZwqNu69393eC+R1EvzwGhJsqPGY2EDgbmB52lrCZWQEwBfgDgLvXuntlqKHClQHkmlkGkAdUhJwnKVRcDt0AYE3M67V04i/TWGY2BDgceDPkKGH6JfCfQCTkHKlgKLAZ+GNwmnC6mXUJO1QY3H0d8AtgNbAe2O7uz4abKjlUXCShzKwr8H/A1929Kuw8YTCzc4BN7v522FlSRAZwBPA7dz8c2AV0ymuUZtaD6BmOoUAR0MXMLgs3VXKouBy6dcBhMa8HBm2dlpllEi0sM9z98bDzhOh4YKqZfUj0dOnJZvZwuJFCtRZY6+6NR7J/JVpsOqNTgZXuvtnd64DHgeNCzpQUKi6H7i1gpJkNNbMsohflZoWcKTRmZkTPqS9x9zvCzhMmd7/J3Qe6+xCi/y5ecPcO+dtpPNx9A7DGzEYHTacAi0OMFKbVwDFmlhf8P3MKHfTmhoywA7RX7l5vZtcDzxC94+N+d18UcqwwHQ9cDiwws7Kg7Tvu/nR4kSSFfBWYEfwitgK4OuQ8oXD3N83sr8A7RO+wfJcO2g2Mun8REZGE02kxERFJOBUXERFJOBUXERFJOBUXERFJOBUXERFJOBUXSQlmtjPsDI2C3nuHtdG+hpjZwrbYV7C/l8ystIn26Qfb0eih/p2ZWamZ/fpQtg22/4WZnXyo20vb0HMuIjHMbDyQ7u4rWvk+Ge5en6BYSefu17ThvuYB81rxFr8B7gNeSEwiSQYduUhKMbMTzexlM5tpZivM7KdmdqmZzTWzBWY2PFjvXDN7M+gI8Tkz6xu0F5rZnGC8jOlmtsrMegfLLgvep8zM7gmGTdjfpcDMYP2pwbplwbg9K4P2I4OMb5vZM2bWP2h/ycx+aWbzgBvM7JQg3wIzu9/Mspv5sTPMbEYwzslfzSwveL9bzeytYNyPe4MnujGzrwXj5pSb2aNBW5dgH3ODfU4L2nPN7NHgvZ8Acpv53Pce0ZjZTjP7sZnNN7M3Yj7boWb2r+Dn+dF+2387yFpuZj8I2s4zs+ctqr+ZLTOzfsHf8exgne8HuV8K/r6/FvOetwSf+6sWHffkWwDuvgroZWb9mvt3JCnA3TVpCn0CdgZ/nghUAv2BbKL9tf0gWHYD8Mtgvgf7HgK+Brg9mP8tcFMw/ynAgd7AWOBJIDNYdjdwRRM5XgYmNtH+GHAdkAm8DhQG7RcS7Z0B4CXg7mA+h2iv2aOC1w8S7cxz//cdEmQ8Pnh9P/CtYL5nzHoPAecG8xVAdjDfPfjzJ8BljW3AMqAL8M2YfJOIPhVe2kSOlxrbgzyN+/of4LvB/KzGzyz4LBr/zk4n+pS5Ef2FdTYwJVj2MHB90HZxzN/x7GD++8HnmR38PW0NPuOjgLLgc8wH3m/8XILt7gM+G/a/W03NTzpykVT0lkfHh6kBPgAauyRfQPTLGKIdhT5jZguAbwPjg/ZPEO0sEnf/B/BR0H4KcCTwlkW7pzkFaOq6Sn+i3cPvZWb/Cex297uA0cAEYE7wPt8NsjT63+DP0UQ7KFwWvH6A6JgmTVnj7q8F8w8HPwPAScHR2QLg5JifsZxoVyqXES0WEP2C/+8g00tEv5QHBft8OPg8yoNtD6SWaDEAeJt9n/nxwCPB/EMx658eTO8S7dZkDDAyWPZV4Cagxt0foWlPuXuNu28BNgF9g33NdPc9Hh0f6Mn9ttlEtFdhSVG65iKpqCZmPhLzOsK+f7O/Ae5w91lmdiLR34BbYsAD7n7TAdbbTfSLObqR2anA+ewrDAYscvfmhund1WIIs8PY90X5e+AfRI8UYrlFh769m+jRxBoz+35MrrODPOcCN5vZxCDXZ9196X77aylOc+o8ODwAGvj490RT/UUZ8P/c/Z4mlg0k+vfW18zS3L2p8W1i/773319zcoj+XUmK0pGLtFcF7Bvi4MqY9teACwDM7HSip88Angc+Z2Z9gmU9zWxwE++7BBgRrDOY6FDW57t74xfZUqDQgjHgzSzTojcB7G8pMMTMRgSvLwdedvc17l4STL8Plg2yfWPKXwK8yr5CssWiY+R8LthfGnCYu78I/FfwOXQl2oHqV2OuyxwebP9K8J6Y2QSip8YO1WvsG5L30pj2Z4DPBzkxswFm1seiIy3eD1xM9HP95kHu61wzywne95z9lo8C2uwuOzl4Ki7SXn0f+IuZvQ1siWn/AXC6RW/vPR/YAOxw98VET2E9a2blwByip8D29xTRawIAVwG9gL8FF/Wf9uiQ1p8DfmZm84leF/i38TjcfQ/Rnn//EpzWihA9UmnKUuA6M1tCtBj+zqPDAN9H9Av0GaJDPEC0B+6Hg/d8F/h1sO5tRK9VlJvZouA1wO+ArsF7/5Doaa5DdUOQcwExo656dCTFPwP/Cpb9leh1ku8A/3T3V4kWlmvMbGw8O3L3t4he4ykH/k70lOh22Dtu0Ahad8eZJJl6RZYOJbgjq8GjQyIcS/SLuuQgts8FXiR6gb0hSTElDmbW1d13BnfPvQJc6+7vmNl5wBHufkvIEaUFuuYiHc0g4LHg9FEt8MWD2djdd5vZ94j+Zr46Cfkkfvda9MHOHKLXy94J2jOA28OLJfHQkYuIiCScrrmIiEjCqbiIiEjCqbiIiEjCqbiIiEjCqbiIiEjC/X8fK2KoSEi91QAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "\n", "neb_res = job.results.get_neb_results()\n", "fig, ax = plt.subplots()\n", "ax.plot(neb_res[\"Energies\"])\n", "ax.set_ylabel(\"Energy (hartree)\")\n", "ax.set_xlabel(\"Image (zero-based indexing)\");" ] }, { "cell_type": "code", "execution_count": 10, "id": "00f32c45-a0a1-499f-8b58-e2847022ba84", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA4sAAACoCAYAAABJw3tmAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuNSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/xnp5ZAAAACXBIWXMAAAsTAAALEwEAmpwYAAC7bUlEQVR4nOzdd5zdVZk/8M9zzrfdXqZmZpJJJ0AIKbRQpaMIKCisoGsHe+9iR8Vl117Qta/YVywUAVGKSpMmECCE9DKZfvv9lnOe3x/feyeh7K7+RGaSnPfrdclkwszce+fc732ec57zHGJmGIZhGIZhGIZhGMbuxHTfAcMwDMMwDMMwDGPmMcmiYRiGYRiGYRiG8RQmWTQMwzAMwzAMwzCewiSLhmEYhmEYhmEYxlOYZNEwDMMwDMMwDMN4CpMsGoZhGIZhGIZhGE9hkkXDMAzDMAzDMAzjKUyyaBiGYRiGYRiGYTyFSRYNwzAMwzAMwzCMpzDJomEYhmEYhmEYhvEUJlk0DMMwDMMwDMMwnsIki4ZhGIZhGIZhGMZTmGTRMAzDMAzDMAzDeAqTLBqGYRiGYRiGYRhPYZJFwzAMwzAMwzAM4ylMsmgYhmEYhmEYhmE8hUkWDcMwDMMwDMMwjKcwyaJhGIZhGIZhGIbxFCZZNAzDMAzDMAzDMJ7CJIuGYRiGYRiGYRjGU5hk0TAMwzAMwzAMw3gKkywahmEYhmEYhmEYT2GSRcMwDMMwDMMwDOMpTLJoGIZhGIZhGIZhPIVJFp9FRL1EdIkg+k964ud/LYiuN78LwzAMwzAMwzBmDGu678C+gugowTykAbBtr7I6Om7wPK+fKpWNPvOVQfz/fJOA/wLzzTzNd9cwDMMwDMMwjH0cMZu85J+N6AQh5X1ULC6fT7TqHKKBo4XomwMUBTONax3ezzx5dan0vT/5/jWV6b6/hmEYhmEYhmEYJln8JyM6VUh5Z6JYPP6llrX8nUIc1Af0O0QdEsiBKKUBO9K6WlNq43Xj45ddEgS3PMK8xfxiDMMwDMMwDMOYNiZZ/CciOl0IcZPX3X3kG4U49OPAUgeYraXsI8vqFFKmIYRgrYmjiEkpVlG06a9K/fZ14+MvvDeK0po5bX5BhmEYhmEYhmE860xTlX8SorME0fWyq6v/X4To/yTQawEZWFZWum5epFJp5HIC+TxRoQCRzRK5riDbnrvcsk79EHBaH7B9uh+GYRiGYRiGYRj7KJMs/hMQSSKSNHs2F4g6PgB0SMAjITyyrCR5nkeZjKBCgairC9TZCSoUgEyGpW0DQgyc3t39rdMs61SPaKP5HRmGYRiGYRiG8awzicgzjCgi4CIb+CVFUfZMITLzmR0NCCGEJMuy4LoWkklCNgvk80CxCOTzoHQa5LqAlDYxZ/+VGUXzKzIMwzAMwzAMYzqYTOQZQvRBIvqQZLYY+EYEnDhbqeJrAKEBJkATkYYQGpbFsG3AdQHP23VzHILjkBACIHKW9/R8Yr6Uy+n//umGYRiGYRiGYRjPLHPO4j+I6F0ELCfgMgYYRHfMEeKDFwFzXhNFP+iSsgaikIAAQBNAE1qH0NqFUgSlAGZAa4AIIAIRKSaSKebcfGZxJ9G3mPnVerofq2H8PYgWEfAcARwI4EQGKgDGCFgFoKiZPdO8yTAMwzD2cUSfEMB+BORbn+kGUGXgY8x8g4kVpplJFv8/Eb2RgFUEfAOA79j2SQcodcw5RD2vEWKwS8oihPgxgBIBdQA1MFehVBlhWEGzmUC9LqF1vHDYaADNJhBFDK19AA0GtCazrmjsgYg+JYB3MNHbCSAJkABAAGkAKh7bNwjglcz8mHkjMJ4xRI8QcAcB88F8zLMyyUb0HwIYJCABYAmAdQy8hZkfNWPb+IcQvV8AncT8TvXEz7+egDcRADAfaCaTjT0S0StascInGBAijhWmAl/F/BYmuk0AlzDz1eZ6Ok1Msvg0iE4g4AICfs7M13L8ubcK4GgC1jDQBeBKAr6bsKzTDmJe+nqtB14s5SJXyiykFIpICyCieDWlBKAErUuIogn4fgq1mgPmPBoNCQAIAqBeB3zfh1IlAJMEDBOzAuBM11NhGH83oqMk0T12InHgQDZ78WohBlYx5xcReUkguR6gu8Pw0VsnJy/cFIabK0T3E/PB5k3A+P9GdDEBTMyf1EIcQsyQREIDZRD9qwDeBeZlz2hATTTYCtYVhPg4AdJqBToagNK6pohuFMAyZu4y49v4uxClCLiEiD4viKSXSHzD6+5+ZRZIIQx9TqV2NBuNo8vMyif6KIAHmfnbZpwZewyi5wmiP1iW9dt0d/ebFmg9fxVRYR6QYWanRKTuq9WufKhSed6I1h/xAZjxPU3MOYu7id/8CyBaYxElk0TSy+fPSCcSFo2P31rz/VJd62aTmRPAKQdJufzVQP+LhJiXtKwUhLAVUSC0LlMUlXUq9VWRSDxEzAeBaB6AOZCyH7Y9C67bA8fphGVlAdiIIkazGcD3yxyG42BeXw/Dm18wNnbNTVr/XDEvM78oY0YjmkfAEZbn3T43lzvtXCkX/yvR4CJgAERFAGkACTBrMNfGgcqvRkc/9Tmlfv+YUtcEJmE0/l5EHgFHkRB32oCVcd356Xz+glnMiQGi7ES9/seNlcqVJUBXtf54APy3fiZKmoi6CZgrhXjYS6UO7M5kTlsADCxhzvUDmWEhsGl8/NIHff/RYSCsaH2NZj7MrP4YfxOiN5GU/ymIcp0dHRccIcTg2UDPoUB+FpAF4DKRNQTgZq03XFupfOuuev3MEa0vjKb7vhvG30JKWwiR7igWzz1DygUvJ+pbRTSQBDpBlAbggDkKtW6uB+q/rFS++JN6/dcPa/2cgPnnJlZ4lplksYVoNhH1WkKszScSq1ZlMmecRdR/CFAcANICSG4F+K7JyZtvbzabxwnR90Ih5mSlTEEIoYkCUqoKpUa11ptY68fGksk/17PZ0bnMvQQsANEcIpoFIXpgWZ2QsgAh0gBcaA2EYYgoqimgRFo/csfw8CXnKXXLNub9TZBhzGhEvQScYHnezSsLhVMuITroRGA+gB4WoksLUSAhUiBymFmQ1hDMgFLDa4C/fnRk5M3XhOGxDeavmwuS8TezbccSoru/UDjnGKJ5LxGibylQ7CTKJZjdEJDjzPwA8/hvyuWv3NBoXLVe61rwj/xMIkFSuq7rLtgvkzntxULMewFR30KiWU48KeK1JkT0duboRq3X/tfExMX3BMHlJeYj1P/5A4x9GtGpwrJu9bq6Dl0OrPoAsPAUIRZYQCeAPIgyABIgcgAQtK5UmUd/VSpd+pV6/bj7mK/2mX/2jF5Hia4i4HYC8gQkOe7B8CVmftxcr42/C1EnWdaEyOUOOMh1j3w/sOjFRAvA3AmgQwuRZ6J0a3xbglkQcwjmiXuD4O7LKpXPXuv7KytaX2bi4meRSRZbhDjMTSZHl2SzJ1wIHPxiIRZ3EvWDqACiLIAkAIFmcyOq1VEAAkIIzeyT1mUoNcJKbQ6YHx4FHlkHTPzMdVUlnXb+07IqLtAPYJCI+kDUDaADQmRBlEI8g0LQOmLmuiIaDojuef/Q0M3fjqKX1JkvMS8KY8YiGiCgYHne9qWFwtGfE2LlMczzGBggKWcJ2+6GbWdhWTaEiBs6RRE4CKCjCKTUSDUMH3p5pfKqqxqNhxWzay5Kxv+K6NXkuj+xOzsXHq318e8DFj1HiPlWvEUgD6I0EyUBOARIaF1TWk/cVC5/7AtB8Kebo+g1deYC/r6qphBE98HzfphKp1edbNuHvZ1o8ZHAIJg7IURRE+U4Xu20CLAovq5PTDDv+Nro6KsvD4KXDjG/yaz+GE+L6EKyrO+6nZ0rzyVa8UkhlvUzzwLQq6XsJCmLJGUKRDYAwUqBtQZprUip4btKpfd/tF4/8WbmR33gOZr5pP/vaynRfAKOISHqEOIWG3ASRMIBrBBAEEXbG8zfArAfM68y12zj/0RUIODFmDXrV4cBx31biBVLmAcY6IMQ3ULKTrKsPKR0IYQAM1hrcBiClQJpPb4NePAT5XLhp9XqYFnrrBl3zxKTLAIQ4hA7mdy6NJs9+sNEh5wFLGGiAbasWWTbRViWCyLBSgFRxNRsVrlW2yCCYCuiaAhab60DD23Set39RLU7hcCI1q9eRTTv+YnEY/PT6e8IogaAWQDayWKR4nKSJMd7EomBkIFKk3ntDba99R3bt/duCsMrTKJoPOOIfkjAXQAWgPlNf/dFgOgVBNgEvBjAt8lxru8tFud9UohDXgbsp4FBKeUA2XYfPK8LnufC8whWa5e078dNnep1qCAARdH2x1KpK89Zt+6Ah5mfw8xkLkzG0yIaFK47ZHV1rThBqRVfFOLABcz9RNSrhegmKYtkWWkSwmKAWGtwPCkB0np0E7D2U8PDH70pDLfstsK4eyexJ3cV2/V3246c7u7gGOaTPkK0dIB5gIlmkZTdZFmdFG9HkCAiaA3dCnJEFPnQeutXx8df9hnf/8Q25hPNdd14AqJBSiY3Uzo9+/lSHv4DohUpYI4mmiWlnEWW1Q3HycNxrKnrqFKA74N9vz3xtu2uycn3fLBev/teIKwBTQV8TQE5APcT0OqohxIDBOYP6yfeh2UETACog2iubVlDHYnEfrNTqeMWArPnE6VnAflRIto5OfmNR8Jw/SbmxvYw/KTP/EIzpo3/FdFzhG3fMauj47BfSXnoCuZ5CuizheiDZfXAdbvgOEk4zq5YIYrieKHZhAoCIIp2bvO8K966bVviet+PfOABzfwNEy/8k+3zDW6EeJtIpx/qTKeXvoZo1plAhyZKkW1npevmkEx6cF2ClPGFudFgSJkG8yz2/fvAfPMG5vv+AOBmZh0SveFE5s4XWlZ/QcoUMc9Hs9mA4/yUpdzBQBPMZQKKaCWLILI4Xm9pVIFtN0bRuo8OD799cxS9zlx8jWcM0fkEbCUh7ofr9iWYySISvmUNBUoNAMgx80t41/9/BDHf3mrw9EER7zksA/gcAQMJwC0AtxeJUjnPO+QIIv9F8f/jSiKHpHRg2y48z0ImQ0gmAafVq8n3AcsCmCGVQsTct7haXf12x3nFG33/3hCQ5uJvPAVRHwHHoaPjpmOVGvgS0az5zBkNZKQQGek4WThOBq5rwbJAAEgpIAjAzSZ0GHbOCQL/o0TvvksIVYsTQQIgEHeejj+Otf+kqb9LuTWh1M8Ok7J7FnOaiVJSyjQcJwvXTcJ1Ldh23NNPKYggABoN6GbTZaVmv6Gj4/JoYuLCDxJZDebjzPXdAAAQEQH/gr6+nYPVat+/EfWmgAQDniVEApaVhOsmkUxaSCTi6yhRHEg3GqBaDaJehwa6Dk0mP/TOIFh3u1ITjzPvaACfHwGaAcAVgEMgbACiDjQV0RU1QAP4vAJ8AC8mor+SbV+bchxxWDr90rOEWHA8c+8CIXpcoMCARwCQz//HMHN4u9brryyV/uMmomgb8GXFfLO5dhtPQXSUcJw7rUJh8D1CdK9kTirAtYk8SOnBthPwPA+JBE2NcSAe4/U6QASpNRRzz0CzeeaFzB8eBzasA+oVIrsKXEHABgCTYP60jn/mamK+7W8aj0RLiflBM3b/B/t0skh0ANn2q6xCIXFoGBbPi2ffHJLSFbadQDJpI5MhpNOAbcfJYrUqQMTQulMDx4la7ZFbyuU1t0h54UmWlTtDyr6cEAkSgjTgs9aTol5fyGF4cJBM/nmbbY+6zNUuYMQBUgy4Or5wB9t8f/JGrXdcXqmUHlHq9Rp43XQ/RcZeQohVludVOoH0XNc9vS+ZfM4As5cRwh6t1+/Y0WxeuRWob5DyhrLWDR1PZpxEAJhogUOUzDM7nYBdJDp5FnDiHCF6FhMllgCF2UL8tJ9ojR2vNgoIQRBCwLYlXJfgeUAqBXheHOQ0m/H0SBgCQQChFDOQP7NQeMulQ0M/fZzoSGIeNBduY0ocUB9Onndlj9YHflLKjvnMcUAdBxtPDKhdN07a2gG1lJC1Glip7lmelzxDKcauRPDJfz7dz2dY1iYI4RJzCoArhPBgWR5cN4Fk0kIyGY/x9uRiswlICQFANxoOMw+8tqPjozds3br9GiJNzMKMcQMAKJv9MarV/ouEKCwGPAXYFpEFIWzYtgPPc5BMAun0rjEWhlOJI2kNWa9LFqLvVNftPE2poMkcaOazh5lVnVlvZo7qROURrcUI81ATKG8BHo2Yv7gTiBTAZctqyHx++WrLOvpVwMKDieaAqAtERRYiy0QeiCQxW91a0xlEHc8pFL7wM6KLP1+vH/0o0UmK+XdmXBtTiIiYWc+ald2fKPUq5pQmsiRgQQgJKSUcx4Hriqkx3r5+B0E81rUGoiiOFZTKHG/bH3OUCtZovWUL8K21zN8eB6JhwJ8gOmgS+EAInExEPwIAML+Eic4lYATMf+An3r9FBBxNjnMi2fYOEU+IWGg2H1Fav5+BEpgv3afH9D6dLAINwfy+XBDMfr4QiQKzZCGEFELAsgQcJx64mUw8cJWKB28YEppNCMeZjWbzRcdKufz5tj2nKIRNADTQZK0nRBTt1MyPB8yPPaLUo3dEUfOPqVSUdl0cBPiziCYEs674frQxioJ76vXVtwTBQ6OtF9Y+PTCNZwaRINctZNLpxUelUueeLMSC1UDXQqAjS5S2AC9MJlePJ5P8MPPQzbXaD6+vVu9fx+wkme1eouNnEZ3VL8SSBUDHAqLMHCJ3UIhEnigLIRwIQbDtMsWBNgNQiI8TUBBCQ8p4FdG2d82IM8cf23YcTBOxEiLbqfXs1cDzHgeuFMCgaQZi7GYeMd+ue3t7zxEieyhzO6CWEMKCbdtwXRuJRBxsJJNxkBFFu2aplQIpJTmRsDgM/db3ZTz1astP+ZhoErb9J8T7ER2Kf66EZdmwbQnPI6RS8c9tTy46ztQeXRGGUMyFRBT1vD2TufTGSqXkmwlBAwCYWS9c6BTrdeelzC4TxWfNERGEIEhJsG2C6wKJxBPHNlGcNDabQBAIsizBiUSCtc54zAIADcZ7Z2n/OKwIoHVIzAvArOpaH6eYMcTsh8zBDte9W0h521KgvxPoZCAHKfNkWQVh2ylIKdpl1hyG4CjKZ6Io+apc7j9yRO/9eK227GGiUzTz9SaGMdoEAEXkvIDISzFL3UohQRTH1VLG103XjSdDEon487YdJ4qtUlQKApBStptI9BwvZeJ4pfqh1MoRrZtDWgdrtR7bzLxuPfMvtwDhCPDXHcDEJNHFVeBtGniciN4B5s+2qqbmCsuq5RznoYWJxKLBROL5BaDDApIol28sR9FdO4HwcSE+v5n5Vg10YF8se92nk0WiTey6ySIgVsd/VwREIIogRAQp40C3PYiFiN/8LQuwLBCRAlFqnuPMhxCkmeus9YRQaodmXl9lfuxB4KF7tB69jqj2Z9/vHldqByWTnGRGB5GWzGHV94OJINBR+0xHkygazwQhBGUyhY5EYtXLLeuIlwLzD4qPh+kiIeKuekSeA1i9zNSjdffqdLrvGN//3QNRlMwCuXlCePOJEt1CpDwhUhDCAhGIiAFoBiImmiQhtgMAxV3yQgABmH0APpgVtJbQeteEi9aYCtDj8j/NFH+wAsAPADkNz5gxs6U5nXaTmUz3SwEXgIgLR4laAQfBsugJwYZtxwE1EM9QN5tAsylISi0bjQ1QqrFbKSjv9idPLTEyawAStv0Que7DAPZDXLrHIGIIsWsypP2zHSce48Cu/bmWBRGGmoHMYYXChfPK5fMfJlpDzAeY6/0+jIgEM+t6PXUwIHrjz2rRnnQj2jXxFq/CTMUgINr1sWWhVdHBgnk7lNKayAaRDSFsSGkxYCN+3UgGMsQskvHeRZFmJjAP7+84YyAqAMiAKCUsKw3bzsB1k/A82Z7gg1Ig34dsNMDNpqOBWefkch9WQnz0zZXKZ8ZhzsQzdtGJRMLO57OnAADFR9G2bnHMHd+4NUGyK4HUetf4ljK+3gME267CsmwwZ0hr0a01dSmFg7Tupyha2FDq5B1a+1u1DtYzb9jJ/IdNzHdtAfwhoLqd6IVjUs5yU6nMCsc5+kTPO/ooop5FWhcLQFYCCaTTK+pEcitz7a9S/ud1WtPvarUFo0TvBvNl+9T43meTRSKXpDiK8skNOQXdD/gENAE0wFyH1nUoFSAMHYQhTQ3aMNwVBAAKlqUBVHUUTQiltjHzuirw+B1aP/wHYOLPRJUHmOeOM9cACERRxKXSZJWIau3vwsxMrUjZMJ4JRMdQKiVzmczCNwIHvB7o6wQ6BVGOLSsPyypASg9SxklZFAFhmPKiyD4xmTz+hGbThxAe4tUTpjg4VgzUwRxorZvQuk5a14UQG+F5QwR0AGig/RpiriGKqgjDDILAQqMhoPWu0pJ4szoQlwP6FH9tfQ4AAkwJqrFLfH18DjKZ9FwiWhZfM5V4crAhhGoFx5i67R5MSwkIwZDSBbAFvn8bdq2Ix7d4rm73xFEDELDtRxA3KmPEEyIR4g6pAeKgnqcS1/Zt9/shBIhIayGS6ShKrwbwMPC4AA4wK+j7tvi9n61VJCCYIwWEIp548wE0obUP5hBa24gb7cVfqNQTJ94ADSJA601oNNYS4EGIBIAEhHABJJjIYSEcInK5nUzuSihvhWWtAbCAAJeEcFtl1h6SSTm1lcCynlBmTQBEo2Fr5oFzMpk33lypXPl1ooMFszD7cg0wM/f29iaZMZ8oAhDsNr59MDehdRNKRVDKQRjGY6xVevqkca5BpEnrcfj+JIRIsxAJktKFZTkgsqF1IcFM85TCXK31MVHUxUotL2sdbtO6sV3r0S3M92/yvD+KbPaUk4iWHAwMpIh6QJQnoiyIEiCyE8yyqHWwLJv98LHMW48Q4vJvVyrFvxJ9Aswf2mfilD0+WfwXIvJaH48CuOpvWJWTdArNxSBtTte5fmotdH+XlkANQBXMVWhdRhRNwPezaDQcCOHA9wnM8UbbZhMIwwjMNTBPIAzvIeZ1VebH/sS89nqg/Cfm8gPAHJ+5zPEqS/Np71d7FdGsJhrPFKIVgsi2Uqnuc4HBNzJ3FYnSJEQalpUWjpOB5yXgugK2HX9NFIHiGWIbzLOJaISjqAIghNZNrXUTzDXSukpAmbQuMfM4mMdYyi0BUdljToCogvh1VIFS7ddRAlJa0DqBIBDtn4dGA2g2GVEUQOsKA1UAjdH4Hjlm8sTYTY9k/kPU3Vvotghe3Dk6pN0DauYGlAqglIsooqkSvXZw3Q424qROs+cl4fu3PxhF24bif2ENsACwDeA6pjrcMAtBRdsWL2Lutoia2DW52ITWDWjdRBQlEATxygu3WpYFwRMDHeaQ4z3vek78vZPT8WQaM0t7hdyLJyAau93aE9dVhGENQeCg2YzDtt0TtvbEm9YRgABSNpTWV49o7aeJUhJIuPFYS1F8pEwSRCkQJQjwQJSAZYWUSv0ORC4AF0QCQojWfjKJROKJZdZaY2qsRxEonkhPWVFUeH2hsPbHExPvKAGfn47n05iBLMsSSkV2fN2sY9ekch1a16BUBWGYhu/Hexi13rX3u9GIKzSiiKF1iPjcxSE0GmvAnCEhiixlnqXM0q7kMQnLcojIhuumBHM2pzVyWuv9leqD1l1Nx0FEVM4C3UxUgJQ5SFmElElIabXKrQlRlIBSiTlR5L02m333YuALH6hUqvcSvYKYv7tPxO57bLL4diL6HDNvRVxYAQB/YObPEdFGAF94muTrcDqc7uA7WPH13EuzRCHgeYnH9zud9YijxQQIkwAmoXUGUZRCs+m19gPkYdsumAV8n9FoRAjDKpQaZ6JhDfz2OqUe+jVQvgOoPQIsCpgnNFEPMYf/a5JoGM+keAVmENns8KAQ+78dKHYALhO5oj1L7HkeUimBVCoumyOKgw3HAcUlohYrlSTff4DCcJiZq2CeAPMYgBEA41uZhzcBYw8L0djiOD4BzjuBaoZ5EkR5Yk5B6wTC0AWRBWZGGObQbCYAWIgiIAg0ms0AUVSF1hMAxgBMPBaHH1sJWDytz6UxMxCtIsDXaZo7L53tfxlSE9wOpOOV7nawUUUYVuD7LmzbnSrRawfUcbABKBVB6ypbVo67u192+fbt77pL64rCruXECcQZ6NSEhRCiz7J6Tgb8DqDOQJ2AGpirrSCnjGbTg5RJMIupILq1z6aVNMYTjPH9bu6x773GM20lA4AKnT7pRADiCbfdJ6+VSiIIPDQaFoAMwtCaCqTbZc5BEEHrOjOXYduy1tXVdenIyFV1pUIphJwNWC5gz2O2JXNiAPBSRIkcs5sD0tJxAkdKn4ADEa/etMusear01XHilcX2+NZ66r2jVWbNmii3fyJx9vLx8XNvIrpdMB9hVhcNBEGghBABUQPM1SdNLCcRhgn4fny+otYpBIHV7irduo4ywjCEUpXW15dC5t8MRZHqJupyoqiLgE4QdVCc+BUgZZ6FyJCUKQiRIiFctm2XbJsI2JIkGmWgF0QpkjJNtp2F46RaXa2p/fOptSdY+37SiaL+k9Pp1ynb/s9XjY9/ZJjom8T8mr0+nt8j37BWE9FtzPyCVCo9r1jsLWjtekR4Q0fHyKgQYz9VSv2AiF7aSshWE4nZeDHfwXfwC+mFdI9csyAn5p2RCQafl/jLwP6jmbudCfkYJIYJlIHWHqLIBkDQOoTv11tBgI0oYoRhE2FY1cwlSLlzo9Z/uRgorQGWhMwPaaJOAgDmnXv9ADJmHAFsokSi64VCpBcDjiaSst2Mw7YtOI6YagLS7qrX7jjWOmaAwjANx7Hh+9cEWg9tBLb/lai2XojmxjAMh4HzLCGySWY9JwgSC8LwJuG6j7cSygwAF1o7iCIJIJ4NDMMapEwDcFqzdQpB0EQUVZm5JIDNkVIb7mTWTPQlBk6Y1ifSmF7xxMd3wHw3F8SC43M44O2e33NsOXkdQBWKA+oKmCtgTiGKdg+o8whDB1LSVEBbrzN8P4JSVWYuE3O5EUXNe5mf+xjzxgC4ywdOTcXHC1wfAKckADDwB1/rk2zme7YT1YvMZSaqTK2eC5FAELgQIh7rUZSElALM1Go8wgiCCFFUAXMJwDiAyZHpfG6NmYHoMAEM6yzt99p8Yv/nR/YaBkoSKAHIoj3pFkUufD8O15SK4PspCGFDa4EwjMd1GNahVImZJwmo7KzV1l2rVHkjc0RaI80MSYQiMwiQOUDYAKWIRNG27Tmum3sfMC/bKu2muMS6fYvLrNvzJ+0mZULsKrcmAsV7zz1LqfQhAG4Cxk2FiAEiIsdx6h0dhUeBoBeY0EBeAmlonQCRgyiy4fsSgEYU5WFZKRBZ0JpaYzxEENSgdUkDJSKqbdG66xKtr4ssS/cwu/sze0uBrvlEHZ1KdRJRF4BOEqIIog4WosBSZoVtV+F5fwaRRUASQrhkWQm4bgKJhEQyGe99b082to76EkTgRsNRwOBJnvfiVzHP+zTRKwT2gf25e2SyiESi96UdHacdmsmctJ+Ugz2AkwbgAxOPFQr3fLWj49fXEt31UiD6AB0jbos7FGCxXLAqSdkXL8KBx6TkrHkk8z2QOSi9mEr8GDPtACFJgA2lCMwRtG4gDMsgipNF5nbZXF0TjUvmzT/RuvIgcKoGbm+VlI7u9QPHmLH00qVLE6OjE8+NuzaC4mYdCkQaQvATOpO2k8VdbdgJts0QQupkcikRrbxuYuJH3yA6x2W2u7Xm+UKI04BV3UJ0Z4jsIlEyq1RXkvnTAMaIOQkiB8wSWjPCUEFrH1FURVwC5YBZQGsFpQJoXdNEFaH1hlsmJ39+P5EjmD9lZqP3YYfSoRTvSzxB9MjLXjeLlr0uKebu76JPBup+qlvbSfIkQFkAKSjlIS6d2z2gTkOIuCJkV7DRgFIljoON4eFKpcd1nPd/lXn0Wt8//ceu+6HvAvxL37/hCtf98CVEfrPZvP0TjvOJtzFf1g/cAGBSAJMAUmCOJxaJZOvnhgiCHITwwCxbe4FD+H4DSpXBPCmItgLY+tf4kf6VgROf7afXmAGIjhDMd+qCXPCGIg58TypcMqsUjZNvbSXJ4wClAXitfYpxRTSzQhQ1IWUGRO0xphFFYStZrDAwScC2h+r1dwdEI6cxVx9g/sBm4EUp5se3ALc1mKPWOXREzLCYg9lSZt7zxLLueB9ZXAbbRBR5CEMLvv/EMuv2nrJYe/8wMvHf3Wf7aTVmJuH7ftTX13cNER/HPAEgj7g02oXWNqJIgJmhdYAwbEDKVOvfJJRSiKIQUVSD1lVi3o4gKHQAF51H9Nw1WluPMqd+zbzzp8BbhG3LAa3tAWZ3EZBdoVTnXKDTi5v7FeA4jbjhMBag3U3bsmw4joVEIj4ur90gTet45b61pYGiCELrhGDuuqir62PfY/7IdqKXEPOP9uq4f49JFomIPM+zFqfTxxyTybxhpWUdPgfo7CFyikRIE8EG+OBE4jCf+QV9RL+8vDd92ae4Or5Kzj80xekLloqFJ2dkd5+Q+RxkWjfJopquadlMwBM91PC2M3gLE0iAVGtTeQ1EkwASAOz2vpMIqAmt169Np3/8XaILFPN1zDy+Vw8WY+ZjZl64cGEGwP67zxDvasQRN+OIJz52vekDu5pyxE06NMXlowcvI0q9T4ileaJMmsjJCiEyUtqSyNZCEBExAceJKLoVlnUTAw7ijnsMZgXmAMx1EKVBlACzDUBCaw3AbzXN2Qrg0c8HwZYq2nGGsU9aQkvEJCY5ZednzaE57+nG0hcnRV+vbRWFlBmk1FzUrHUMjBEQB9TMNrSWiCIGcwSlmgiC9nizoFQ8yRdFDWhdZqAM5pHbm82fPu55334O0cAW3yc4Tt9zgfJ638cVtp06DfHy5aelTDwnig4tSHmtAiZt5tFWsO5Aa9lqfBYH8kFQAVFcbh2voIdQqgmtqzr+uZs2ed4NfwUuAfPbzaTIPohotWC+XXfIBW/pxYHvTovZvZ7skjI8BDusTWCMguAh3mQjW3teFZh9KFUDUQqA1/o3bpVWN8FcAnMpIGr+IpM5MsmcvjgMaz+u1T765UTivT/SurrR9498s+NceRFQmQuc8yHgnKwQfz2I6Hty9/LuJ5ZZl9BsupAy1dqOEz+OJ5dZa11D3P+hblYUjd1pAGDmXzLjYqJmhnknE7kEWIiPdtGtzuntMR4nkvHEclydFI/xCQJKkPL4XCrVd5JShx2hlCwrxRWltjeYr1inlN6otfsgYN0E/LQB/DBhWfYiZvewZLL7ha47r9ja5kK7H91hWQTH2XVETburNVE8KRJ30wYFAZQQHbNte/XzAHwTOFsgDrj2WntEsjiXiE71PPuxfP7UVUJcvD9wUC/gdAtBPZZFecuCZVkgZiKtM24QZJ6XTL7+wcBeeIpYMrRS9D8vLTu7SKRTsDLcJIGKrqMWlfS4qohRHmk0Go2NcDhDYqKfWSlwk4iqxFwi5gyABANWq01eUzBvqyn1+Hu3bt22gbmfmXeYRNGYEZrNpgTYZQ5A1KR4prjRmiGuIYoaCEMHvi+mLpLxHsJdjUAARXGJX/cc1503R0qbiSQLAQKUACIw12QcvDS01g0VhkcjlXqYpJwgZq0BJeIENd6jEB9k3g6ABAClgYCBMcG84d9HR2+8GXg1M7/fvJb2US+kF9Ij/IieJeesnMuzL+0WCw9Pyu6kJbNCSxuBqsBpdmrheoIxzsQugZ4cUAdQqg4h4vEWr3DvqhRhrgLYRkFw901K/XWs0fhEHbAaQCQajY9vAsRGgNFofO5eIKwDDdVofOrBIJiYcBy3l6i2PzDEgCMACWaGUvHP1bqOKEohbrsWJ4taxwEQUNVE41KpTV/ftu3+CaLfCubTTLK4jzmSjhTMt+luOf/t/VjyzpyY3ePaBSkEg4KFmLD7ddPaAcGWaK9Yx+M6BHMTRO1A2gFzHMLFjW0ammiSms31Vq129WNEzR4pj5kDqJ1ApG3bOQYY8H0fnEyefD4QHFAqXfzhbPYHFwhx2LuY/5Ai2j61nywu7y63Gje1954DYZic6qLdPo4mCOK953GZ9SSAsb06cDb+Pq2O/8KyrPWdnZ3fBOFt4ElmtmV8EMATJ5a1TiNeoHHQvsYyh4ib4+yA1gsgZYGSSUsyywwzpbXWpFS/iKL+pVGEhtZUVapa17prROtzH9LafsS2738gkfjiWXHu2j6+I2rd2sfT7Kq2irtnP6WjNsWT6hYxW61k8fi9PmbZI5LFWQAeAJYdQPTWfuDgHGBlhBA520bW82AnEkRxpw5GGApqNGDX6+mFCs8/xe1WDSokSSZ0HYIqusblMMSoroshHg9HefgvkzzyTaWav8tx9mBAvZ9o8jDmJpjHARR1vNLhEZEkICCg9Ciw4ZMTE1t+F4ZKMX9yrx8oxsxHRHQoDhWJhX5UrUqOZ4eriGeIa9C6CqVKCIIUGg0HRElE0dQm7lZ3UrTKRutgrkJKhxzHBhBQHIQ3OO6OWmfmCmldFsyTgnkUwE4Aac5kxiIhKhZzoJmrrZK/AgFTyaJmBgORIJoQzKOXjY7e/tkwPKy+Lx52awAriegeZr6Sr+R5cvD8uZj9gQ4MLE5aHcKyUiIiQk1V0YgaOtBVQUECcOpgjDDF5daIV8xDRFE7oE6gvfoSr6BHAHxFVCal1l5fLv/+91ofEWh97QuAYAxQURT9+FSAJpm1Jkq9tVVsV1fqv98JiCSQeR6gvwyUwLyN431aCq2uwa0GNvGMePz2SmBWAJoRUJVar3swl/vRd4aGXsPAKwHsmIZn25guq2m1uI1v0/1y7ttnY/G7CqKvy7OyUggHATegAp8dzKZmfgfAFWbeivj85zhQjjuwpxGPsTiQjrdLRQqok9Y7NgfBTd8Mgt7HiX6xgPnX3wWCe4BcsVb7wDZA3A8w1eufnAtkdgDDIop+sVyIQrfnRfcB4Yq40ViWgBS0jif3gkC2GtoErTLrBADZmmSMEAR1RFEZzJNEtANKDd0TP+LxZ/9JNmasKIp0taq+nvbci2D7CcIoM2tqHaeBeGK5ttvEsoN4YpkBhBoocxAE3GgwEa0lKTMkZUYIkYSUCViWDde1HIBspTijdVZE0cELlcKKKEIT6NXMPy0S1REniD7aZ0IzN6e6aYehQBi2z4B+6vEdRO09vdRJ9DYCPr3Xxy0zPlmURCLhOLnubPa0AnB0ErAcIuFZFiVdF3YqBcpkgFQqPpE2DIFKBQyIjJ50ViSb/OemgK+rGI8UD+um2MaTegjj943rsW/Uuf5rjYMnEbzHD0bO3+Z53mO5XO58oH4u0JgLjFmti7JmxhYg+NHY2ND3osjfqPVrQ+av7fWDxJj5jqAjJDMrlzydH53bIbyKFXfVKyEObNPQOokw9CBEfAHWugjfT0II2W5sg0Yj3vui9TiYSwjDAGHIWusqa10R8Ur7BMUH1QwD2MnMY1uAyXuJJu4LAv9erTNpIv02ouYhwBAxT4I5w0CK4n1lUhApMNe2MY99dXR063+G4YGTzNeb2eh90EIiUQNYEqUXisE3L8bc1xVEd59n5aWUDkJEaER11KM6T/KwGNUbrmxGk99NON53gWaBeYcCIhmvYjcRT45kwLxrZjrWDqi3rdf6j5/W+mebgFdqoLKGWbWv45t33bO1o3GJ9CImenRUqc+I8fErftnR0XEgoN4AjGtmRUAg4vMXq4irUHZPFpnjMKMumTeXlXrkDY8/XhoBijDVKPuWI+gIcTvfrvvlnLfPxoJ3dYjeroSVsYS04HOTgrAe1XlUNhvV/0Qe6wj4DKMK5s0RUBNEVQGUEZdeJwA4DIjWRoImMY+Mav3op4Pg8Z8zH1JifngSGH0QcW0plPrNCWiVymn9xXcoReNAoJrNS/5C5HEisX2tZcmVwDgDWQCJ3cq70SoRbEKIamvlPD7zUakIUeS39kuWSOv1o8nknXfGj/oiDbz42X+yjRnocGK+Xeetxe9KN1clVG4b150HiRABHDBQB1GF4rGXQjzGLQYEMysQNUmpbVSr+QiCIQA9IOqCEB1aiCILkSMp0yRllqRMtY/OgG3bBEiPmV2ti0TUIeNJuiZ2HeNRh9ZVRFENQZCA7zsQYtd5pr4f38IwvprHyWUNQK23WBxNjI97jWl5Sp9FMz5ZXA3gbqLZKSmPTxBZEoAkIikELNsGeV6cKGaz8W7UIIj7dfkBwW9if3eMrqta2KYisRkV3s7V341g8stVXbuFwXVmDogswXweE/1LVK/XH200Gpek0+nLM5nMfK39PFEdzLo0MVFZFwTRODM342Vzkyga/7j3EIk0QEcAvGpXVy3aAFAK0Pv/L8esPJ+eL67iq/TtfLvqosHu+eKYV3nlwVeXxYPpmrsekicAygNIgjlulhAEAswaSvmtRgluKyCIu5OGYQ1aT0LrEfj+JkTRqIiPzNgJYLjBXHoYGHsYKN1GFDystX4Y4CrAkVKIRkaa5HneHfl8/iwgOhOoHgjUcwBJQFSZsQ5Q1wTB2E9LpXBdFJ0UMF9nEsV9wP1EdDdACwE+lpkPooNIYTZrUZq7kro+Osh9ZxRkR9qx0pKkRIMDNKKqrqiKGOYdNILNnxnlTZegwiESqX/1kL6CqJEFdmjmOgMVEa9iZ7BbQK3jDrtNwTyyXsp1bxsf33CbUs+PmCf/x9cWcz/v+ng/Br6tiL4zMj4+flmxWLQAfSEwwUAzYi4TUCCiLLWSRc0sGVCSqCGZRzYzb3zL6Ojjd2j9Nmb+uHnv2IccRoeJO/lOPVv2vWU+5r2zU3R3pay0JaTdShSbUU2PyUne8P1t/Je3Y0ijsyc3QSQ/KqD7wDvAXFJAjuNkMQnAEgCIyCfm8qPAuo+Vy2O/CYL9m8zfVkRCAys18FfNHDB2W8Ymol8hrg5URPJbQojvSdk5n5nfSNToZx7SRI4ARGs/2a4ya6Jya7+u3VptjPdSArWIqGQxb/re5s2/GyH6lmAeN2XWBoieL5hvV2nR+5IesfJfc5iNZOMIUcYCGnauZ8YoE+rEPLbbGPcASBFXbzQBlO+t19df6fvjFrN7CFHmYKDQq1SXVKoHQA+ALhB1shBFLUS+lThmSMq0iM9OlAD6AaxF+2z1uOKqPHV8R7NpgyiPKLJgtYpW4on09t7cJrSuIJ6UbG4fG/tBPd7u4z3tY99bzPhk8TbHgZPP9wvmZUwkFABFFN+EAEsZJ4nt9o5A/KdjA9KmXrvC/d4WfUXFv6lGlc+VtX8jg4PW2rcAgHZnsFZdtWbmoFKpbKtWq9vjbxi/r2vNmuh4AdwN5rJ5szf+IVkiOgSg37e69T4Jd7Q++AkRnfekhPG59FxxLV+rr+KrdBd1dfVh/9csokNfm6SBflvOst3IRs19HIwxIrT3brW7/DKYA0RRvK+LyAWzhFLcKsOoA6iwUo8gin4xwTz2EFHpj8DE7UrxIwBvjg8wjwv/49fSKTLuFAlmZm40GhuazeYXk8nkN7LZbJ4Is5k1BYG/o1wujUZR1GTmCDiJmK8zAcVe7DEi+j0gLmJWBzPzwbu1GX+AH+C5onDUMvR+bADdR+WsvOVYCamFQE03qB760aSuyC3YOrmTt7yqpMevZtbhC+hsd8Owe5znFXko97v7AppcLjBB8d7YEb3b6gsRQRA1AVR/H4brPzY6WrojDDlkbvzd13Bmjoho6/j4+MeKxeIWAG8gqs8CfDCPgTkJIAEiWxARmJsRc+1qYMPFIyPBo1o/XzP/h3nv2IesolViGMM8R/a8YSEG390junpSVkpaUqLBTfLDZlTVY3KUN39vk773dfPwatb6d5nEzpWHu2729p0dVz0MhVcTVXrjANdCqykIA2obc/Df5fLkN+v12sPM0MwPtLqya8b/3HRDY+p1+PxI6zV6dHTjxu7u7q8BfAnzpGa2qHUcRmsvWbvMOn7PaDcmaZd3AzVL68e3NJt3flbrdwB4IwOv/ic/u8ZMF3fcZW3bqQP6aNkXOzBopawsPMsh1z+auZHjseT1pOwxCApFXBFltW5gZrUBCH5QKk38V61W2xR/10AwV3JS7pjPTAcA9lFA+kAgcyDQmVKqWyrVgzDsBtAJITq1EAUWIk/JpBaO06D4TOgM2sd3RJEHovaReQphmIUQcVf3uNxaodn0EYZlaN3emztaiu/Prwg4b6++rs/oZJGISDqOnSLqVETpEHFVcxNAgxlNZmitIeP+5LtaFmkNYka8JVbg6FRYs4vheW/Z6I9/G2fSK/ErAAA/TZDOcVDO7Z//5PvzdF9jGH+vtxBRmZm7bNu5aNasFQuEOLWb+eACczIFwAFKDNy7TevrPu84a84D/K8T0S04ia7gG/S1fK0mIvsQOuSVS8WqD+ZpoM8SRUlWRikpIXkxjUbzuGKt1xYPE8iheJaYW4eD794oIU4W47bVIeIZ4ooUYsuvmB94m9aN7YiDDm69PkD0QgH8cur18TQlpExEqNVqtWq1WiWibcDUhAwh3ogA5hvM62kvdiqRuI5ZLwLUkkTCelFHR2+nUv1Joua15WDHw3U6eSV1vX+AivtnrBQsy6GINCpRBdXIV6Nckpt4xx3r5GMXccgPMrMu0LyBKrbvOM5558YE+hLBRHPD1uLP/s1h60Lm5koiPxu/kQsAqDDzPUDzvycny1c2GsF25v95pf5vEa/G0M5arfaZrq6un9Xr9TMBdQpQ35+5lgKgmbGJCHdr3fjx2Fj19jBU9fj949d7dUBhPNHBdLAoo8xEjdcsxsB7Zolib8ZKSktKNNmnIGxGFT0pR7D5uxv0Xy/SrEMA6KL9xIA9cHSB9u+nHe49W/I/P8SxvENYhPsTqSRzc5JZPVCt1tZUq42dwBc185v/5rG1exzD/KvWHDypkZGRH3R2dh5BpJ/PPKqZlSAKqF2u12r2B+apMmsAOorLrIdKzJtfMTn5lyHgFDA3zVg3AJxBRBbNEou/2UPzOrIyr5N2kpgVN8IRTaHtU9D3n9yx86cs7JVA1MUcshBcZ1YPVyq1u6vVxhBwGgNbCXiQgJcq5v8aJqIRgO8E/O8BPgFjA0JsWgLgSCC5mii9BMjP0rpLKNULolmwrEk4jt6tm3YSzC60thCG1Kq6ChAENQiRaE2yA1EUIQwbUKoG5hITjQEYuil+jA8QcN70PcXPAvoH3zf/qYhI2LadKRaLr88Jccl+AO1HRPOIaNC2MZhIYHY6jUw2C0qlQLYdJ4u1GnhyEiiVeazZwDql/D9qedT7hrbc910ALzUJnzGNziQS1wvBh+Tz8w9z3Y8vJDp7IZHTDVABQJqofToVV7QujzH/4KKxscv+5PubAYCI7NW05Owiet6eE12r3FaSqGVCBmSjrjUqKtJjtBk78zfGB3ShgLj8oouIOgDkAGRaTUBcMAuO3/gDDdQF84ZxpR44ZWTkjvuYD6E4MLiVW/ni33XRaCWH7cTy7/56Y8+VtSz5/M7ORfsL8fJe5tM6ieYXiNIZIByNuHF7OelWw6JXsNOwrCR8lpiMFI+oJm3lCWzgnd/fpLe/WbOuBBTQShx1SJ99wDVJ2fn+q5pfvuJA+YIPjfGmcKu6/SOe5dm9uVndvhMOCqE9Zh0y17ZVKtWttRpC4DACtjDz9mdk/BGRaFWoELC/KBQyiUxmyGk0BE9ODtfCsB7GkyOnCeAWjstkjb3dY0R0PSAeBHAjoBuUf9mB1PeRftE9J2sVpC1T8JkwGTbVsC7Lzbz9ew/rBy9cjuV8P/5qH0QnvqkLi759yeJL6e2Pv/dnO9QjP9+gb/wK8MQJ7F0TdyTiUE7/Q+OrfZ22bXt+Z2fnF4jouXEjsqQGOgRQpLjxzdSeXI7P722AeaTEPHLuzp13/575TGb+hYmxDBxBR8jb+XbVLQc/PQdL39Uh5iHl9kgpU6iFAZeCcd6u77t1I99wOjPXn7xA88Qx/tTYofX51teMaebikxd6UraN+UrhcMBeAaSPSCS8hdks0nHlyQCIBmjX/sciiPKQMts66ivemxt3tY6gdQCta8xcYeZ1gZT3LNyy5Ufb/qlP4gwxo5NFAHBd1y4Uiy90ge/OJfKWEGEuEc2WEv2Og9mJBHoSCaQ8DyTl1AmaulZDuV7n8TDEGqVK11jW/K9t3jwx3Y/H2Le9l4g+K4Q8tLv75JVEX10mxGA/M/cAulMIyguBtBCQRAytAa0FmGm71vf8rtF4709L+YLLxbd2iM7DU6JgkcxothIUkE01DZRVqEtRE5M8KcawI2ykNq61UuGB8SpLWgM9AuhEfCDu7u2pCURQ8bEwQxXmja8eHr7zN1qfHDH/xrzxG38zj4jeBOByITLndHW9c4mUb9qPqNCDeNRlAE4BSBLRZCT4nmaRHgs7UGOJahRip/KxDhOVDTz60W169KuatU9EDoG8Y+jU9Fxv2Q1N5nkPNP94whr+8+3AVNDwtNUiRMsIqBGw/hmfqNhtIuQJ33u31XM21Sh7v88QiXEAn3nS73pAZM87WPR+bA51LczLHNlWknwWmAwDNaQrYiNv/6+H9GMXvRwvV9/ib0Wr5LHnL3IO/y8Wzp+21O99yW24dvtuAbN80o9lPMOTb+3x7LpuX7FYfCcRXglQllkgfq9IAnAJIAbCVjDPfxRix7uGhkbu0foMZv6lGe8GBAmhWeucLJw5B/v/pF8stjJ2j3StHJpac9mv6CH9+PB2/eBJo/zIw4fQIfJu3P1030r/rWN8t+sxPd3XCUEylbKOSKe7riCi2UCWgQFB1EtxXFQAUR5Emd3OMo3LrQENrX0GaoqoLKNozdeGhr7yRqIXCeaf7/VjfkaXoQIAESnbsibCKJqYAHpHmZEBKKEUnDCEBBBqjYLvIyUEJDOiIEDD9zEZRTzJjCFgaFTrpiIiOdOzY2Ov9UEiukwIWtjZecwcoh8NEmWzzCpLJHNSWgXLQtq2IS0rPtxHayAMmYNA9+loZTfyVy+kbjsji2TJtGbL1QFBVHUTtajKE5HPE7oiRjASjWL7Xyq88/OqFl7dlex6CxF/CKjYQJOBSR2vNGYISHJ87iELwBfME5u13vme0dHHr9b6BJMoGn8Xm4ieB9CXhcg8v7v78kVSnjeHGfl4nIuilJQTQnhCsCCC52g6xZ2guY0mfldN8P2R4g2YXLeJxt4xqmvXa+awn/q9k8Uxv+m353V8Y/Kbh70o85oLtqltF4/CH3NoJy3HGXG7PDxlFRvAVAL5T7nu775t4Wk+b5o27eWyRHQoQDfFY0Bel0ymzsvlZieY++6p8rJlVHz1QsosLEiHHEsInwNuhKEa11XayiM/eFive+MoRoMe9Ni30q2L1mP9Tw4Qhx5TDkdfVhaVg1nxtqW01HqQH4zaY/yfqT2eiWjb+Pj4+/L5/E+FEBcwR8cD5QGglgISQog+EC1HGM6ujY3NfaNSb30IOIvMiqIBAIIErcEa7qf+nm7R+x+zaJaTkymkLAcRB1BBXdf0jmic1108yo88fBgdJv/Cf/mHx/eTr8dPTB49AKdys/n7+zMZ/AfAn4u7C29rfV1ERLvvzW2voEvE834RA74G6lLr9Y8r9ej7iT6yTySKwB6QLPq+j5GRkS3FYvHuCtHpQ8ycYmYHIKEUFDOaSqEkJRJEsLSG1hp+FHFZKV0ikhuBX/9s69bGh4jkJ8wbuDENHCJ6GxHJxYs7uycnv9RPlE0CKgVYKSmRtm0kPA/S8wDXjU9/1RrwfUK9IdGoqcM9dnYEHo+TyyGRqOkAlaiBSeXrUd2gIZ6IhjF83ySPfrHBtZ8xc0hENDIy8tmOjo6HhBDvJApXEI3ZQEkwOyByAAgNRGPM/g2+v/3fSiV/jVInMfM1+8RF0HjmPA+ga4RIHtnT85U5ROcVmaMMYGWFsPKWhZxtI2HbIMui9oSIFYY4wG5i0KnxvLrS94jKR78wUf8dgfQcmj13FMNb5tirr+61Cv92dvoF3/xNdPUrALxo10+9c+qa/j8lb4bxTHsTEVWJ8JBt597Y2XnKPK3/ZRaJ5UWgOyuRWJAQeFglOGCXXGmRQsB+FGFC17GZR3/4gN74lrNwaHQUjrSOkMsvH7QWvvAs76xTP1v67Bv3x/IfPMh3/wkAHuQHo2f7sbX2lgdDQ0O3Abgtl8t5Utpv9P3+d7nuKV3p9PHCthdBqbKXydy6uFSKHmP+avhs309jZlqB5WIlVopO6vzSXBoY7JQZztgJAWg0w7quq3GM8rofb9X3fedwOlzeyXf+U+Lyp3s/IKLG6OjoTzo6OuYQ4a1AqXVWrg/m+AgkogzFK+ntZBEAVGtv7vadWm966cjIbWWiDxD2ETM+WUT8i97CzNf6zMcOE6VdxCevaK0pYEZNa2SiKF4wZgYzc5OZawA9rvXI/czfdInoLXEHMMN41r0EoMuEEHNKpbcXhTjAA5QDWLYQcISA6ziwEgkgnQaSyThZZAaaTZCUYGZZ0GXeP1mhPzbSVI18HlORGtYNsYMrwQ6euHecJ75a4cpPmTkEnrCnyt+5c+cvHce5vlAoriCio4HmEqBRYOadRPxXrau3lErRo74fBXEjjhtMwG38XV5OJH4jJc3t6rqgi+j8HBAlADtJhIRlIek48DwPlEjEEyJSgloHfOq6jRTqOEMEcnmU+eDPROP2ZTx3wVJrzvf7rOe8+eP1H3/5qGjlwSO6/OAL6AVCQeEqXGX2vxrT4s1E4puWZR/S1bXsaOBji4DT5kuJTsS7wdMAH+SBj3HrtD4Q9GBT81pf8OaoobaJyR8+qkbeGXBUskl6AoIPECfc3iuz566rrf31AA446EG++09vpbfSF/gL0za+Wwlja0+YFwDBFtuulSwr6GKushDjbNuwfD91ACCuA26PgBdM1901ZohlRGI9gJzsfOtczDqtW+QoZyWEJUD1sKaaqkSj2PjXEm9421baSnfwHc/qAk6rQdno+Pj454rFohACrwUqSWafgTIDOcHc3pvrAJAMRCBqsNYj67Ueevnw8OY7ic4SzL/aZ3KKGZ8sti5Y9cmJievy+fwhdSEu2MpsaUD4zFxjRolo16ksABQzAiI9ShSsCYIP3jg6uu4UgC4ze0eMafJJgH8wOJh1Go2XJJhZEgkBQAgBYVmQlgXhukA7YXScOFm07XiFMQiAwKGFzjj/oSp4c8h6O9eam7i0ZYTK3yzr2peZOQB2JYntvVLtJhzMXAfwp9btKYgOI6JOMgG48fcSRLQZ4CsGBzsy9fr7coB2AOkAsIWAIyVcx4F4ugmRRgNCCLDWgqNI9Ut5wI+7O15yyY7eH/TK9PCwGv7OSsw/4vfRHa+c7sdpGO8jEl+X0lnW0fHClcBnFgox0MvMOWadk1IUhaC0lMISAkkwVrg1LPFqeKQh8EgUVRup2hWlyksqh4rBzCly6SV91HnzD8Obv3+C9qNhXRKbkR4/iA6yHuAHnvUVxSfb1cDpDQTUtgBXDwOTC5mrECKCZSUgpTefCDaz3Zzu+2tMryIVRRU2p6S7ejY63jxAhUTBSoqEJeCrBvthHeO8rTqETW8e4qHSaXSa+C1++6zHG61xPTQ+Pv6JXC73VynlhUIEBwKjDvOkw+wyYBMgOK5idUOl7O9H0V2fHB9PbQJOpn0pUQT2gGQRmEoYN0yWSpfn8vnuqhDP2czsVQGajDtIUgqAC2gRL6vQJDC5MQi+9ejo6HeOB+g3JlE0pklrq6zKzpp1sEU0mwBwXPsJBUATgaUEW1bc0de245WX+MRjwHVBjgOWNjqtAIu8nfhNGD76OJe/N6rr39PMO1s/5wlJ4u6e1FWsdRfA+F82gxvG34oBMRvQnfX6eUkhBi1mLYgEEQFCQEgJYduA58UTIqnUrgkRKQGlQEEA8gOwamhS7ivvor/8hMMFb96p/FdtQWN8Ha2jt+AtdI0pjzamSScRnUskeorFQ+cI8ZXZRLlMvB/Xykgp8paFrOPAsm3AsuL6N6WQCAJaYTd5RcT5LZH7b4Vzv/+yP/+8szyn0XlCTuRffhavuuA36i/fqXB1hl6DP8sAhrTumGAuE3OVAJ8ty4NlycVAygYON6/LfdwX8Fn+Er4yUOOhS+ZS56xOKyky0iLFIcKwqSp6LNyBTR/aptb/6VA6VN7Fd03btrBWwlgaHh7+XjKZvDaTyRzHzMcCwclEciFRl3Cc+RDiAPh+hz8+fvq/R9HBm4DPEvO/7nNjfY9IFoGphPEvVCq9M5vNXtSU8vQhoHuSOZlilq3jl31NVKswbyiH4ZfqIyM/ORfQl5sg2JhmRETdvb1HgIgiZhUCIgTQPjs0QHw4m2TelSQ+8RuABMFn4uXJiM6Q/sf/Y6zxCx3vS/wfk8Qne5r/x7w2jGcC27YtLaIzJDOYiBXiun9NBC0EtBBgKXdNiLSTxSiKP7YsQEqpCSpJtOAcZ84rbmg8/Mn7WN3a/iHX4BozXo1pczZA37KszoOl/PdeolwC0AnASgiBlGUh6XnxdgLPAxwH1Dr7mZtNcFUKNOp6NuPge37T8/6HKpMfOMATL9ocbP/KJj2ePRJHiVNoAV3Pj0/7iuLudk00do4x+8NalyOtqxbQICFSsCx7QMpT01qvH407Shp7u2uJ6OFWI7ERANcB3ABoJyhdoOJ7l9HAql6ZEnnLJUkalbCqq6pM27Dl54+qR794Cp0ipjNRbGvlFcTMwwB+BuBnUnofdt1Fb0+lTs3a9uFk2wOs9Xgqm/3hiomJxgZmOe33ezrsMckiMPWLfQzAxclk8hcgOjEUYnldyjxrrZVSGzRw+7jWV+cmJna+AIBJFI2ZgokyEeLjjRutW72157amFNJhiITvg9rNbQDA9+MS1CgCa40mAA3CSjcZ7U8l8cnWquW0PSjDQDwJMTAwkFZKLVJEiFqTIe2JEJ8ZITNY63ifYvvWnhhpX6YJUCC4IkJ/0s9+pcniv2g/62X86IwKoI19zwlE4g8Azers/Jcc0aEJQNmAZRHFpda2Ddt14xLrdDpOGIUAwhBUr8flHFoJ6Eh1k/yXi9yl45+q3f5hD+qczVwvTffj+7+dWmP+7y3MlZrW1RzQgJQhS2n3ZjJ9AxMTz9lM9D5m/vA+t+qyL5hPRO8DaFIIfM62bZfZseLjDHVHFAXbmIVHqXMXUuGsAZFOdliecAWoHtV1LaqJbdhx3+NyywfX0lq6ATOnJ0I7YUTrqCMiMaF1tc48mWEeI6Ik2zYLy0rNA0gAhX0y3tqjksW2RqPRrNfrfwLwp0Ii4brptNUIAlSr1YZSSgPAUUTiO6b01JhBBKACAFUAtdafVWaUtcZEGMJrNiGFgMMM8v34i4IAqNehfR9+FKGqNcoAtjPTfczqIrMyaMwcLgPFJoAmQO1JkbrW8YRIFCEZBHB8Pz4TV7Xec5vNqQkRrTV8gEMizHepcRqgf4q1Zowb0+4PANuOk3WZX+kRQcaRI4gIJCVEe9W8vfc8mYxLrKMIIIqTxnjyj7qsgLqdkQsO5r5v7Gw89tBxROLmGRyvtIJpDcgh5kqVuZJlrgspQ23bruO6xQWAuANYtE8G0nu704loO5HzKdueNc91+96Qycydo3VPPh4Xjauq1fWTjTDRrTMXzRfpnm7Lo5QUCFWTG2GTd/Lojk089E4OsfW1eC3NtLNnWwljq8eD3sZcKzFXe4AmpFRs2zYsy9qPiCzm3D7Z9XePSxZ323tFrb/7APz259oleX+aYYPRMIh5NCLCJIASgCyAFDMSSsENArRm6ZCNIjiWBQHEAXSzicD3UY4inmDGCHO4TeuSBnjWtD0aw3girbUGEMXNx3dNiFSYUVYKqTCE6/vISwlL6ziwZgbCEGg0oIMAgVLxhAgzJhgBI16ZNIzp1I4rir29cwWwhOINT0IBUERQrVJrFiJOEKWMy6otK04U26XXUsZ7eImwX6I6OlhMjizfADzGe8akH7PapnV1UutKH3MdQkRs2w6ktA8ggtXTM88ki3uRQSLqJpLbXLd3pW0f/fxs9vkDzEd0AX0dUnoZABYzL8tmxs9ydXR3NVUU7Fl5ywJziFro6zFVrj2Kbf/e1HRLP/rFzXzzTB8jO4BGmbkC5hqImiSlhGWJhUDWApbtk7nFHpcstj1dMw5zxpYxQzEz8+z+/tsUs18GnDFmThORxwxHqbj+AUCoNWpBgJSUsJkBpRBFEepBgIpSPAnQNuZtt/j+ukMAfZIpszZmACKidDrdTKXT2xtAcZKZS0SUYUaSGZ5SsMMQ1GxCMyMTBLAtC6I1xtsTIpUowgSzGGbGBNEjDKgHAT5ruh+gsa8jIqLOnp79IIQdMesAEAGAgBk+M5qtUmtLa5DW8co5UfznbqXWTEQBAAmadbBSne8Cdj44rQ/t77JN68Ykcxla10DkU9wRVR4EHGj19x9kOqLuJfqIaBaRPZpMHnRSNvvKA4Bz5wMds4ioSIQ8EaeI4BKRDXTMTjAWiQhrgxAjOsRY6PO4rvFGGv35I2ro2y/HKv4237UnJFrbtW5OtJNFoElSepDSmmXbByR8/2cl4FXTfR+fdXtssmgYe4p2+/GOjo5HHMe5wweOHY4bI0gLgGAGK4WIGU2tUQlDJISAzQzBjEgp+PGKC48D2Mz829tLpdG3Mbe7mhrGdBPVarWRzmZvioClY8ycR+ukKq1hAbsmRJRCzbKQbI1xaI0oDOMJkSjiSWbaCpTWSvnYK4jUlcz8wel8ZIbRIoVY0Np3zlNl1ohLratKoR6GcIIAVqMRf4Flxcliu9RaKSitqcGsQ6KEE0X9ANb8BcCiGXwt362iaxtzc0TrSnySNRpkWQ4sy57neUck7r77WxXgLdN8b41/1EIiSgDW9mRy6eps9qNziZ7bF588wDkhUJCS8pYlXCkhhGBigLXixVZI3c4EPdoI9B2BxGau3LyOxy7tBVd34C8zdnwDTxjj48zBqNZVpXWN4pXFAJYlezKZud1BcNEQ0WJiPnpGP55nmkkWDePZQWNjY7Wurq4vW7a9vMSc3tE6XoABRMzwtUadGSkiJOKDq0DM0PGsta4B2My88ZZ6/YezmRu/BPgT0/ygDKOFATAx/4qB11YAe5iZXSKyAAit40SRGY0oQkVKJIhgtca4ak2IVLTWJSKxWetb/lSvb/gMM75iVs+NGUIAEwpABaAqdpVZV7TGZBQhEQSwGw2kAMgoivfmah0ni80mVBiipjVXmKkMhCNaVxjg/VqdJfcAZa2DYa2rgdZVj7lJlpWCZdmzMpl6R6PxzmGipcR8gnnN7sFcgDamUj3LMpn3zCJ6XgeAFIC0EDLTOh7Gc10IxwGkjMeuUgTfR77ZwAoqcYIIfhDefG25suUzAN68B1zHW9vb6szhDuZak7mebjdysm3H87y+eYB4EPgLAyZZNAzjGdZaXRSVSuX6fD7/LUX0+jHA1cw6BIRPhAYzKsxIAfAA2AAI4IiZm0Q8BoyuIfqKXavdT8zqgT3g4mvsG9rju1wu/yWbzV6jgBcMA9pilkQ0VWLtt7r/ppWCh/gNSDBDxaV8uhaXWU/eKcR3uyYmyp9kxmnT/NgMA62tBHNmzXooAnQZoBKAHFr7zrWGE0WwgwAEIFIKSd+HJUQ8URKGiJpNNIIAk1HEZYBGmSceFmLLRwCs3AN6LLQC6ZA53K51pa5UzQMaECJk23aTjtO7EBCPAg8yYJLFPdVSIvEwkdWXTr+4U4izcgB7gPCIyJMSCduG63kQyWTcyMlx4qO9whBo1MFCwOEaHQwfeYfP2ITMD7aVyxs/TUTv3zNiFs0cbdW6VmWupoAGxcmiJSzL3g/A1XPnYsa/Xp9pJlk0jGdJK6CuViqVL2YyGRkSvXYc8AJA15hRBihH1E4W2QJYM1MIoMQ8soHoK49Vq1f0al17aM+46Br7Fq7ValUAn8lms4t8ogOGmCPNLCMi8gE0mJFpTYi4iN+ACIACdBPAGKA2aP1fd9Xrt76MOfr3PSCINvZ+7e6NmVxu3Xi5PNwAesaYdYZIeAAcZlhKAWGIiBkNpZCREm5r9VxHEYIwRD2KUNaaSwANMT90e7lceeOetZ2AAbVJ62pF62pB6zoRRdq2HWnb3oEArl2wIDSv2T0UEdF8gL1sti8jxIVpZscmggWQFAKWlLBtG/LJR8S0uv3CskCagSgiKKWLjENekkq9+ORK5bOl+DU008d567Wot2pdr2hd69a6AaIAtm3BsuQBAITvW/vcGDfJomE8u7hWq21l5s/mcrlxBVxQBhY1iXiCGWlmJAG4RCQAhIBuAA9sU+q7o/X6T6q12sgOrfe5C5Ux87Xbj9fr9XsKhcIntNaXNojmDsVl1KoKiDIRsgAlAbgAS4DBzCERlZj9LcB/P1KpfG1lvT7x72acGzMIEYlisTjquu5PFdFbxph1EnEFiGztvVVhiEBr1KMIKSGmJkRYKURao64UqszYScQbgP/eVq3WfwLgyGl9ZH8fZn5c6/qk1pU5WtdBFJBtO7Ase38Acr/9EuZ1u+cS6wEuJpPHukSLrHhSWrQ28wFCxOdA23a8ouh58epi6zxRMMdnQzebBN9HIoqQAP7l5ELh8lePj1f+a5of3N+AAYAZW7RuVJirFHf9DUAkYVliEdGgs2PHqftck26TLBrGs6i1iVoR0Vat1JelZf05m80+p8l8uA8sKTHnJcCCeZyIHvSB26vV6h+rvn+vCsMyaz3TZ+aMfVgrYVRjY2NXE1Ejl8u9rgmcOkokqswYY0YG0EkADkCCiCIi1IEdo8AV28rl7xfr9fXXaj3T26sb+x6emJiIFixY8M1atXpeg6h7iFnLeF4PqtUNtcGMKhGSRHABSABgRqQ1GoCqEtFjWt9xTan02xXM6gt70Op56/W9XetmWesqxcmiT1JakFIsJTrAvffe5/nTfT+Nf4iwgCMpThKVAkgBiHa7aSIIIaYSSLT35raPjRECEEJEABxg0cmpVP+bx8YeVkQkZ3hVVGuMb2P2J7SugbneOj7DgmWJuZ2dfcmdOxdVpvt+PttMsmgY06BVkjpBQfCHMAzvA9DleV5HNpfL6jDkSrVaqQbBsAZG94ui8nat1fgMv8gaBjA1IVInouuYeUMul/u9InpBDVjaALLjcTACCfhg3hQCt01WKr+uB8Fts8NwbII5mu7HYBhP1goisWPHjscymcynNHBZGbCZWSsiEQJoAqjtViHiIE4WGUAEqAZA27QeexT4t0N9f/vN2HP2Pu12XNmw1v6QUtVQ65oN+LAsF7ZtDebzg9nt2+eWpvWOGv8IBiCJeZECEAAUID7I3Ec8vputiQ9LKSCK4hVFreOPo2jXMTHMFABaECW7tB6wiB65AaDjZ3Ap6m5jvKR1UNK6prWuC8AnIVw4jlXIZAYGiOQw8E0wv3LGPpZnmkkWDWOa7LbKOA5gPIoibjabAsxQUcSqdd16cA+aeTaMtlZ16SNRFG2xbfuqfDbbB6IBH8jUtWYlxEijVlsfNho7oygaY+ZgG4CamRQxZqhWwhiGYfijYrE4WxO9qQK4EbOqA1QlojJAaQAJAHar1FoDaAI0wlzewHzpo5OTf3i11uHWPWyst7tFah1uVarqK1W3gKYQIsm27WRyub5BIrkVeDcxX2ret/Y83PqP7SM+GmbqiJhWc7JKFCEdhrCbTUgp48RQiHgNstGYOiImjFfS4RNRxGwDwPC0Pay/XWuMN7WOdihV95nrgrkphEhq27YdyyosIKJ7jz/+9wBMsmgYxrOEn5gMmvI7Y6/BzIqIKlEUVYMgeJwAmwHJcYgRsdYRx1sTmfewwNnYN7WqQibHx8c/WygUNkHK99aBgYAZpbjMWmUwlSwKAOQTocz80MYounSsVLr2nCCo7KHNmwgAM0eblKrVtK4l4wYgETuOa7luaiER/XnlSt5TjgIxngYRhT6ACuLjYaoAKlojpRQSUQTP9yGFQJoZMgh2JYu+DzQaiIIAVa1RiRv38diet61AA2qL1vW6UnWPuQEpI3YcS0hpLQFAzH3TfR+fVSZZNAzDMP5pWkkgt2Zsgyf9k0kQjT0OM0dENDJaKn3Xmj37z8Vq9fwAOCsE5leJxAjQPkM0ZKLHG8w/npic/HnT9x9/idbhD/bMRHE3vEXrWkPrGrSOA2nblrAsuQrAD1Ip17yu90wCcenomgh4TomZSwAyAJLM8FpHxEjfBzMjiiIkbRs2EUgpcBgi9H3UggAlpXgSoCGta48C6yJmLtMeMYfQmhDRm5Vq1LWuF7VuQogQti3hOHIeAMrnl+5TY9wki4ZhGMY/XTtpnO77YRjPhNaqeV2tX3//Tua1juN8KV8oLFbM80LmTggxpIAN1VLp8aDZHL1Ca/9f9vgkcapb5ONKNctKVUmpOoCALMuFZVkHAKDVq/ecvZjGU2gC7mTmN04CYpwZqVazJltryCiK9+AyoxFFSLWOiJFax0fERBHqUYSqUroE0HZg/RUjIyPvJ6LZe8bkYHuMb9S6WdO6LpgbrUZOLixLLgYglyzZEE7z/XxWmWTRMAzDMAzj79RueqOBRrPZ3LJzaGgbgJsxdV4bVGuShPaCRHGqAQgRbdW6ORFFNShVB7Pf6hYp58+atSrxhS9cVrn00o9P9901/n4aAIjod1rrDQHR3GFAJwBhAxCtI2KiKIqPiBECSaK4kRMzWOu466/WqAIYI6LtSn23oXX9DoBm7yGThe2uv0oFJaXqrHWDgEBYloSUWDAwsDp7xx03TxDRPlMdI6b7DhiGYRiGYeyJOKZ5F83MqvXn1L9P9/18htWUCncqVQ+1bnDcLTKCbduzUql5g41Go112buxB2hMbO3bsGGXgiyCiEqB2MGMYwAiAYWYMa42hKMKOMMT29i2KsEOp+MashwDxkNYbf18o/GgHoF6xByWKrQ/HtY5KWjdZ6wYAH1JGsG0rW6/PSt14440RM/MLX/jCfSKP2icepGEYhmEYxj/Tboli+7a3JYmgeDklYFZrlao1lGowsw8pI3Zdx3Hd9DwAOPvss02yuAdiZq21VsPDwz9grf/MRNYoEG1lxg4AQwCGmLGDGTu0xnalsK190xrbAGwGcD+z/qvW7/vd2rVj72Sm6h60Ajdv3jyZzWbrYdhY6/vVRhg2tVIhiBQSCfYOO8z6xPnnn/eC5z73uT1XXnmlBoCzzjprr86nTBmqYRiGYRiG8bdoLxpuDMN6U6l6Jm5yk9OOY0vbdpcSid8ce+zFJlncczGACRBdyMw/1kQHjjCHAbNsAFQDkAWQAuABsBGvPDERqsw0CvCxUuK5UbR8HvOvf7QHnSfamgxRixcvntPZ2TlIZNuFgk+ZzE6kUmPo6AjgedkLADrf81I7X/ayl/16YmLiY1ddddWOM888U/z617/eYx7r32OvzoQNwzAMwzCMZ0x7n9Z6rZsNpWqkVANAANtGuwEI9fbuW90i9yYcLy/y8PDwI0KI81nr25hITgDYCKjHAbUe0BsAvQlQG1t/f4AZGwE+PZMRZ+Vy1JNMvueyjo7VZxPxj4hmfL5xyimn0MqVKzFnzpzTOzo6rtpvv/3OXLFisTt3bpI6O0Nks4qy2SQlkymdTCa5UCj0dnd3X9jf33/7ueeee/YNN9zAp59++ox/nP8/9soHZRiGYRiGYfzTbI8ivx5FddK6CSAUUjIsCyuIVjrf//5ZJlncgzEzK6V469atD46OjZ3FzB8mYH1AhBEi2gTQYwCtIaKHiHAvgB1RxCu0xmwAJdtGRzotBpPJL7+iWOzWAKkZvI91//33p40bN9L27dv/tbe394vpdHqpbdtk2w5s2yPPy1AqlUM2m0M2m5W5XE6k02k4jqMymczsnp6eH5xxxhkvv+WWW3DKKafsdbmVKUM1DMMwjD3AOeecQxdccIFIJBLYf//9ccMNN/BrX/vavbLsyZiZdmsAsk2pYDiKGkuiqCHjjqgSti3nLFniZE466aSwo6ODf/nLX5rxuYdq/66DIBhLpVKfmjVr1nfqtdoJzHxoHZhXixecJglYU61W5w9E0ZnbPK+4xnGEcF3q9Dw1SLT/0cCnLpyYeF1ea5wOqGl+WE/RqqumQqFwTF9f30Wu686xLEs7jiMcxyHPc5FMJpBOp5FIJOA4DoiIoiiC7/uyVqsxEdnd3d1fOOmkkyZvueWWX51++uni6quv3mvGvkkWDcMwDGMGO/7446mnpwdXX301/fa3v7WEEIKIFDOrwcFBsWrVKu7o6DArOcazohVcN5mjTUo1QqWaxBwIITx2XTfl+7MG77zzqoqUMgSA/ytwXrFihbj44otx11138Re/+EUQEb7+9a9jzpw5NGfOHJ47d64Z29OsXq/TunXrtgP4Qev2BETUvUGIrowQp2YbDTvtumQ5jpXxvGhQqVdc2t9/zYu3bv3FFUTihTOw8VMymezMZDLnO46z0rIsYVkWWZZFtm3D8zykUilkMhmk02m4rgspJaIoQq1WIyklAAitdapQKHz83HPPve1f//Vfh1esWCHe/OY3Y2BgAKtWrWIAtHnzZrIsiw866KAZ9xz8b2gPalBkGIZhGPscx3HsuXPnFk8++eTFy5cvP7S7u3tWOp3e1mw27/3kJz/56KOPPjp21VVXRatXrzZv6MY/nZRSEhGI6H35/Ox39fcflO7vXyBTqRyHYYmEGPqd59EffV/fcu21v35ASllZvXp1dN111z0lQD799NPpxhtvtFzXTR155JHZE088MZHNZuknP/lJ5cEHH6wEQVD/7//+b3XCCSeYsT3N2itw2HWOKLc+FgA0iI7OSvmlIxKJA47IZuWyQgFz0mm2tcaGUmnsjlLpmDcNDa39ruPQxb4/I5Kl9mPKZrPndHd3X1IoFBZkMhlKp9Mik8kgk8kgn8+jUCigUCggl8vB8zxYloUwDFGtVlEqlTA5OcmlUkk1m83a5OTkJ/r7+z8XRZG4/PLLHcuyXCGEJCImokhr3fjoRz8avelNb5oRz8HfwqwsGoaxx/vxj38slFK45ZZb+Ec/+hFddtll3NfXR/39/bxy5UoTZBjT7sQTTxTVapWYGR0dHbxw4UK68847ccghh+ivfOUrTztGBwYGaHR0NL1ixYrnnHXWWW9YsmTJ0fPnz0/19vYimUwiCILmT37ykz+vWbPm8294wxtu+epXv1pdvXr1HhOAGHueQw89lHp6emzP8xakUqkDXTdtFwoQqdQYksmqsG2G6xZOIhIn+b4fveQlL/lzuVz++k033XTtq171qtK3v/3tqbF+2GGHWRs3bsyfeuqpxxx//PHnzps374ienp58LpejE088cUhrfetdd931w3e/+913n3nmmbWPfOQjM66EcV/SKkt98rWKAWgiEmC+tazU9x4KgvfmG42OtOMI13WpK5HQ/blc537A51aMjf3Lvx98cG0a7v7TYmZOJBKZZDJ5qBBiHu0CIoIQAnF5tQ3HceC6LpLJ5FSyqLWG7/twHAeWZQkiSlYqldO/9KUvfXvevHlzP/ShD50+d+7c1d3d3YVMJhNYlrV+ZGTkussvv/zPyWRy6FWvelU43c/B38KsLBqGscdasGABpdNpsXXrVi+RSBSWLl1aPOCAAxL33Xdfefv27ZNBEIyfe+654aWXXmoCaGNaHHXUUeLuu+9mKWVSSpkiIguAI4RoAAgcx6nv3LnTP/7448Uf/vCHqXF6yimn0Lnnnpv95je/ef4JJ5zw/pUrVw7Mnz8fAwMDulAosG3b7Rl9lEqlyccff/xDr33ta7/3u9/9rl4oFMx4N55xK1asoJGRkZRlWSf29PRc3NHRcXAqlbISCReu65LrenBdD1JKhTgOlwBQr9eDcrn85Ztuuuk/jj322KG7774bS5YsEXfcccfAscce+7rDDz/8tcuWLSvOmTMH3d3dyGQyaJX2oVQqDW3evPkzr3/963/wiU98YuL44483Y3sG2m3VMS2F+MZ+nnfmkdmsc3ChIBYVCpSxLFWt18UDpdKH/n3Llsv+86KLoud//vPT+rtsnwGTTCb3KxQKl6TT6Rfmcjmk02mRTqeRyWSQzWZRKBTQ2dmJjo4OFAoFJJNJ2LaNMAxRq9UwPj6O8fFxTExMoFwuc6lUGtZa/3jZsmXPO+iggxbOnTsXXV1dyGazsCwLjUajPjY2dv1HPvKRLxx33HF3/vCHP2z+9re/ndHJmEkWDcPYI82fP59yuZyrlBo88cQTT1i+fPmpfX19BxeLxXwmkxmybfsv27dvv/qTn/zkbcuWLdt+6aWX7hEzeMbeoR2IpFKpomVZc5YvX35EOp1e5rpup5Sy23XdxwAMhWH44B/+8If7iWjbokWLqrfeequeO3cuHXjggc7dd999ygknnPC5o446au5+++2H+fPni76+PvI8r/1jGHHDCFEqlSr33Xff+VLK6x977DH9yle+0gTVxjNm0aJF5Pu+E4bh2bNmzfpYoVBYmEqlOJVKUSqVomQyhWQyiUTCg23bEEJAKcVBEGjf9xFFUbhjx44f33bbbe89+OCDR2+++ebOxYsXv+m4445716pVq5z58+djcHCQCoUCWZYFtFasAIhqtdp47LHH3v2yl73se7feemvTTIbMTO1rHoADElJ+f1U6vfSQXM46sKOD5uRycJn1UKXiP9hovOATa9bc8KWODvHmsbFp+11KKWnWrFkin88fOjEx8el0On1MNpuldrKYTqenksVisYhisYhWMgnLstp7FlEqlTA+Po7JyUlUKhUulUro6urCEUccgXnz5unZs2ejUCjAtm0AADMTABodHf3Le97znvflcrk/ff7zn/en63n4W+x17V0Nw9g3ZDIZr9FoHHHsscdeunz58sv233//Mw888MDBAw44ILdw4cL95s2bd/7KlSu/8dWvfvWja9euXfyzn/3MlN0bzwoiItd13Ww2u/zAAw98y+GHH35FsVj8YqFQuLBQKJzd0dFxTC6Xe2U+n/9APp///llnnfWd44477tXbt2+f95KXvIQGBwfFnXfe2dfd3X1ePp+f5zgOHMeRiUSCLMsCM6M10UuIt5PoXC6XHRgYeNd73/veXLFYnN4nwNjrrFu3DpVKZVkul3ud53nzbNvWjuOQ4zgUl+YlkMmkkc/nUSwW0dHRgWKxSPl8XqTTaWHbtuzu7j575cqVb7j//vu97u7uQ7q6ul5TKBRc13UplUrJZDIppJTtFSoBwGJmTqVSif7+/vdddNFFi/70pz/RGWecMWOPYNiXtcpUCcCahlJffKTRGHm0XtcbqlXs9H34to1iLuf1e96nluZycwaWLJm2+5rJZMi2bVq9ejVOOeUUAGClFKIoQvvPMAwRBAF830ej0UCtVkOtVkOlUkG5XEalUkGtVkO9Xofv+wjDEFEUkdYaQgj4vs++70vHcaRlWZKZZWu1nQDorq6uVR/+8IfftHnz5jlr166V0/Zk/A1M8GQYxh7n+OOPt7Zt23bAkUce+d7BwcFTe3p6uFAo6FwuR67rkhCCAXAikUgNDAy8/LOf/Sze8Y53fODUU08dzmazZt+L8U9DROQ4TqK3t/e0np6etxYKhcPS6bTreR67rsue58F1XZZSQggBZraY+TDXdQ8+8sgjV91xxx2f7uvre7Szs7PfcZzjWl33SCmF9j6aJ2sFIJzJZI585zvfud95551325o1a2j+/PmmdMj4hxERWZaVchznRY7jHGxZlpBSkpSSbNue2sfVLttLJBJTKy/1ep2EEKy1tqIoSnV3d7/qOc95zrp77733yO7u7p7W3i8BAEI87fqFBBB1dnbOPvLII0874YQTHtm2bVuIp+6dM2YAZo73LwI/HgvDZQ/X66/POA55rits1xVZ19W9udyKC+bN+8DFa9a894bDDy9/9Y47ntXfZTKZpO7ubrFz587kjTfeONvzvEMty+qJoojaSeLuiWK9XodlWRBCQGuNZrMJy7KglML4+DjK5TKYGUqpdqJJtm2j2WzSpk2bkM/nkU6np0qrEb+kBABdLBZPvfDCC39+//33b/E8T8+ZM2dGjmuTLO5llixZQslkEvfcc8+MHHCG8XQGBwdJCIENGzY8ZdzOmzePlFLYvHkzA0BnZ6colUq5uXPnntHV1XViLpfTiURCpNNp4bpuO5huR9Tasizd09Nz/oc//OE/Xnnllf91/vnns23bpozJeMbNmTOHurq6hJTypHw+/7FUKrXEcRyybZtd16VkMiniUr34rK5W8MFBELBlWbbrumcffvjh7h133PF5z/MOmT17dlcmkyGtNWq1GiYnJ+E4zhO+voUAaNd1Hc/zDgdw5/3338/z5883EyPGP6RdWphIJJYmk8nVUsp0fHJL3ARECAHLsuC6LhKJBFKp1BPK9FrlqOT7PjebTSGEmOV53r+fc845Tj6fp0wmAyJCqVRCtVqFZVloT6S0fj7Q6r4phDgOwFcAmC0FMxsDCJj5s9t9f8UjtdrRKceB53lijuOIZCKhZyn18tP6+u75+l//+p/qec8TX7/mmn/ae3JPTw/t3LmTTzrpJDruuOOE4zg9vu/vN2/evJWdnZ3Huq67anh4uGd8fJyCIKB2oiilnBqPRARmRhAEcF0XjuOgVCph/fr1SCQSU5Mk2WwWWuupZDIMQ+zcuROzZ89GIpFAFEUAANu2SQiBTCbj9fb2rn7jG994/de+9rVwzpw5M/KabZLFvch+++1HHR0d6OzspMWLF0vXdZHNZtWLX/xivO1tbzPJozFjZTKZ9swd9fX12QDgOA42btwYLlu2bGo1xXEc6unpQa1W689kMqckEgnZClioPfP3JAKASiQSdjKZfMW3v/3tK0dGRibe9a53PbsP0NjrEREdeOCBSCaTBymlXp9IJBbatk22bQvbtsnzvKnVl3Q6PdV+XSlFzWYT1WqVarWa09nZ+dwTTjihL5/PF9PptJPP59Hd3U3tIKRWq0EpBdd1p74HESEMQ26VSXUCoHQ6ba75xj+Mmdm2bWlZ1iFEtKCdPDLz1Er303WNtCxrKmCu1+uwbZtaSaAUQnR3dnYikUiQEycRiKIIk5OTICK0Vt+nxjYABEFAQRDksevYBmOGYmZujZPtgdaXrK/XL886zvyE57GdSFBfNkvpbNYe0PodywqFB6tbt9526Jw5YsJxuNlsAgAsy0Imk8EDDzzw//W7Xrx4Ma1du5aFEDR37lxKJpNeKpVatHbt2oMWL168OpvNHu553v7JZDLhui5JKalSqUx1Nm1PWOyeKCql4Ps+LMuC53m4//77MT4+jlmzZgGIV8bb13fLsqYm9yzLwtq1a5FOp6G1huu6yGQy8DwPUkoEQbDfBRdckFyxYsX4M/ZLeIaZZHEv4jiOzcyLPM9b+epXv3pWR0cH7rvvvqE1a9Y89Pa3v/2RRqNRv/zyy81F1pgR5s+fT9u2bYPv+zwxMTHH87ylixYtmjN//vxOKSUSiQSn0+lxrfUGZr5/9erVOwBAay0sy+pSSi32fZ+azSYajQbK5fLUqott27uX64nW/q79TzvttMy73vWuGXtBNvZsa9eudQuFwrldXV2HWZZlSylFu1SvPTbT6TRyuRxSqdTU6kuj0aBWSaoIwzARhuERURTBcRzq6urC3Llz0d/fP9VNT2s91VzBcZz2ygzt2LEDY2NjjzOzrtfr0/10GHsJ13Utx3EGABSYmbTWYGZoraGUesqfURRBSolarYb169ejUqlAaw2t48WjarXKmzdvpv7+fsyePRt9fX3o6upCMplsr0S2E0y4rosgCFAqlTAxMaHa38OY2VoJo2Ctf09EX1hbrf5bynEc13Wl7bpUTKXQmc/PO3lw8OLrNm260PO8bbZtI4oiAgDLsth1XQDA2WefTb/4xS/+pth1cHCQNm3axI7j4CMf+QgVi8WeIAiWDw4OHtrV1XV4Mplc5nlef6ujKQkhmIiQTqc5kUjw5OQkNZvN9sTGVKLYvua2S1Bd18Vjjz02NXnXbDYRhiGUUigUCnBdF0IITE5Ool6vQwiBnp4edHZ2IpVKAUC71JV83+/IZDL2bmWqM45JFvcCz3ve82jevHlesVh80RlnnHHesmXLls6ePbtYLBbp7LPPHm82m4/eeuutV1577bW/fPe73z102WWXmYTRmFYLFy4kx3FQLBa7586de+qiRYtOy+fzB3qe12fbdtq2bZZSYvXq1TUAWxuNxv0PP/zwrxcvXnz9+Pi473mek06n3UKhgHw+D8/zEIYhKpUKwjCc2kPTvvgGQYBmsyk9zzNNvYxnXHu1xXXdha7rHiGlzLY/t/uqS3sVJZlMIp1OT3XHE0IgCAI0Gg1yHAfDw8PcbDaxYMEC9Pb2ore3Fx0dHf+vvTsPk7MqE4Z/n2d/au/qqt73ztLZmqwkgQRBQ0AiQUEYxQVHPx0F31kcHcR3ZvxcPpVxXBF4xVcdRkUHAdkJSyAQIIkhKyEhS3c6vXfXXs++nfP9keepqSw6qJBukvO7rr4gnV6qKqdPn/uc+9x3JTAMCtyYphkUWiClUgkdPnxYf+2113YTQnB3d/cUviLU2SKdTqPW1lakqiqDMWaCe1lBUOi67glFQIL0aMMw4Nlnn4WhoSGQJAnC4TBIkgQYY2QYBvK/NrS2tkJbWxvE4/HKKU5wgqPrOmiaBqqqQj6fh5GRkX6Msdff3496e3vpOmaa8+8vIgD4Zc62Vx5SlA+ERRGLksQyoggRUWTrk8l3rHDdG58bH/++w3E1siw3chzHI4QU27aH169fn3344YfNlStXoi1btpz233zevHnotddeIwAAnudBOBwW4/H4grvvvntpV1fXyng8Pl8UxVmhUCgkiiKwLIsQQsTvkYhM0yS6rgfrBWKaJjAMUzlBDzZBHMeppKUGqajHjh0D0zRJqVQCTdOIaZpM0H+RYRjQNK1STbWlpQXq6+srmyL+zxNpa2tr2LRp08Kvf/3rI42NjdYXvvCFaTe2abD4Nrd06VIGYywMDg5+aP369V+YMWPGjM7OTtLc3IwkSUIMw4Q8z2tet27dDFEUpeeee+4Xt956a+7mm2+edoOROjfMnDkTMQwDpVKpN5VKfaquru6ycDjcHgqFkCRJSJIkFKRuMAwjAkBSluV5y5cvX5LJZJb19/fvjEQi6xcuXCgsWrQIZs6cCU1NTRCNRoOCIZWdPlEUwXEcKBaLkMlkyOjoKB331FsGIdTLMExLVc8xFAR2Qbpe9T0YSZLg6NGjsHv37qCSHliWhbLZLOrp6amcQIbD4UqgWPW9QJZlsG0bjhw5gkdHR9GBAweefvLJJ49+9KMfRfPnz6dHMNSf7brrrkP33nsvsW0bdF3nMcZCEMidHCQGgR3P80AIAYZh4MCBA/DUU09BNBqFcDhcScG2bRtCoRDU19dDIpE4oRVBwM8sAQCAwcFBMj4+DsPDw96BAweewBg7mzZtgt7e3ql6aag/neJ53v83quuzIuXyEkkUMS/LqEkUGSkclpK1tR/rDoc7nHg8EZKkGkEQOIZhDIzxRD6f71u5cuULrutuAYBC9Rft6OhgBgYGsKqqAADQ1NTUJAjCis7OzlW1tbWLwuHwTFEUG2VZJjzPM37hu8o8axgGMQwDqaoKiqJoiqIMYoxrCSH1hmFg/xQdBZsjQbDIMAyYpgkNDQ3gui4ZHx8nw8PDSFEUpCgK1nUdKYqCPM+D9vZ2aGhogFQqBclkEsLhcGUTm2EYxLIsaW5urv/Yxz72ubvuuuvQpz71qf3/+Z//iT/60Y9Oq7UKDRbfxhBCqKWlBUzT7FqxYsXNDQ0NnYlEAqLRKBJFsbKxzbIshMPh9tWrV39icHDw9QsuuOApON6bi6LOuGKxCADQlEwmvxSPx98ty3JIkiRGFEUiyzI6qQAIwhgTx3E4y7LmIIRaWJYdj0QiLXPmzOFbWlqgsbERpVKpyikNAFR2prPZLKiqSkqlEgwNDQ089thjqmmazG233UYX0tSbKh6PM6FQqAUhlAhS9arfqk9jHMcBQghMTEzAr3/9a8hms+C3xwDLskCSJJBluXIiWX13CwBO+P9oNIoRQmj79u3HHn/88e/PmzdPef3116fVQoN6+5g/fz7at28fuffee8lFF10UjkQii2VZvszzvDWmaXK2bVcqRfI8XwkUg/vitm0DwzCwd+9eOHbsGNTW1oIoihCJRCAUCgEhBBYsWACiKFbSTKsDxWqyLIMkSWRwcBDt2bNnw2OPPfZib2+v97d/+7d0/n6bCO4vEkIOIIT+ZUBR7o2IYkSUZSJIEgrFYmhMkhpQKHRd/PiGMQn6dCKEIB6Pa21tbVeMj49vXLZs2c9N09wvSZK3fft2Ui6XSSwWk+Lx+KKZM2eubmtrWx2Px7tFUWyXZVmqOkUM1gTINE1iGAbRNA0pigKKokyqqrpF07SNtm3vC4fDiyVJutF13S5d17HneSioksrzfKU2QlBxPZ1OE1EU0ejoqJrNZsuapjWVSiVIJpPQ2toKsVgMJEmqvJ2mrgJiWZbU19ef/973vvd9n/vc5/q///3vmzDN1ug0WHwbW7RoEfOlL32J++xnP3uNJEldCCFMCGE4jkPVA9K/r0Wi0ejMCy+8cOV//Md//J5hmPzixYvpgoI6o1atWoVisZho2/aNkUjkUlEUwxzHBYsGJhQKQSwWq6Qs+Xe6gjQRxjTNWLlcjjY3NyO/eW4lnS9IzwOAyqlNPp+HY8eOkfHxcebgwYO/Onr0qBJUI6OoN0MqlUJ/93d/h8rlMrNhw4Z6z/Minuehk4NDPxUaDMMAjuMgEonAiy++CFu2bIFYLAYIIeB5HjDG0NHRAZ7ngW3bYNs2uK5bOZ08CWFZFiKRiD08PPyNvr6+369bt8773e9+R+f2c0RnZydqaWlBAAD79+8nuVzuT/6353kedXV1oYMHD+Lx8XGQZTne0NCwOhQKre/o6Fgai8XaVFWNDw0NMf6YREEBG3/hDAAAflpfpRqkaZowOTkJ1Xd26+vrK/dug7dg3j4NHI/HoVgs/v7xxx//xmWXXTbxgx/8gAaKbzNV7TSeU133+0cKhS9KosjicJjxeJ6YggAhWWZkWQZRFIMNhOAed8R13bmdnZ2ttbW1XWNjY187ePDgrubm5rqampqLIpHIpTU1NQsjkUirKIopWZYJx3FMkGUUnCKapol1XWc0TUPlctlRFKVfVdUndV1/wXXd1zzPGwEAHWP8KsMwk4IgfIYQstiyLNH/GoTneexnhSCGYVCwKei67oFYLPb9XC43rOv6J0dHRy8nhMjd3d3gf9wJxaCqBVkogiCgaDR6JcdxdyGEjKn4d/pjaLA4jS1fvhw5joNisRi0t7eToaEh1NbWho4dO0aee+45vGrVKvKDH/yAbWlpuXjmzJnQ0NCARFFEqqqCIAjVd7YQABCO49hwONyzd+/e2kwmUwBaUYw6gxBCwQngO6LR6HpBECLBxMvzPDq5AIh/AR1YlgVFUVA+nycDAwPIMAxUfeoS3HPxvwcAQCUVKhaL4Vwux2zbtm3jI488cv8111zj3HnnnXSxQf3F6urqUDKZhGw2SwCAnHfeeXWhUCjteR4bpOmdHCgGqXoYY5BlGQYGBqBcLkNwxwUhBOFwGGzbPqEBtKIoEIlEQBCEUx6HbduIEGIvW7bslccee8y6/vrrz/yLQZ1RDQ0NKBaLwaFDh4jfbqjyu/zyyy9H5XIZ7dy5k5im+Ud/x19yySWM31MOJxIJJp1OtySTyXe1tLSsr6mp6Q2FQo2CIMh+5V2Uz+dRsVgE0zRPmHuDjTp/UQ2GYUA6nYbzzjsPDh48CJZlgWmaYFkWpFKpyp91XQfDMCrpfdWCoqsMwzDd3d07isXigRUrVtC5+23KDxhd4nk/cnj+qpIo9o7zvCchxMqCQPzWKygUCqEgaPT7cyK/CEyUZdk1LMvyuVzuAAAsjMfjXbIs10mSJAaniACAgmsohmFg0zRB0zRGVVVUKpUMVVW3apq2wbKsFzDGAxjjAlS1YnFdN18sFh+MRCL7ZVlejxBajzGe43me6DhOpV0MIcQDgAOO4zyaz+d/CwAHAcDGGPcRQrYCwFcxxqLjOMR//PBHNqqRv2HSNXfu3PhVV101sXHjxrf03+NPRYPFacg/sifbjjcqrZ7sK///qU99Cj322GPMBRdccPMFF1ywZOHChdDR0cEkk0nwm9yCqqqVinsIIeL/0Iie53GHDx+Gyy677Mw/OeqcRQghkiRJ8Xj8apZlOxmGYYNmXUGvrqAASDgchmg0CoIgwJ49e+C1116DQqGAhoaGIBqNAsa40jjX8/5gtgYOh8NIEIRtTzzxxBfvuuuukZtuuolukFB/kUQiwXR1dZHJyUkSDoeZRCLR2tDQ8J5wOLy+tbV1KcaYz2azpzR2Du4pAhxP1SOEQENDAwAA6Lpe2eDgeR7K5TIUCgUoFouQy+UgHo+DJEmQTCZPSLe2bRsKhQLk83k8PDxMPM8jCxcupGP8LLZw4UI0Pj5OAAAaGhrCbW1tzaIoxhBCtuM4Y+Pj4/mjR496AADXX389uueee04ZD2vWrGGeeeYZ/Nxzz+FkMsm2trbOa2pqem8kErkiHo93SpKUFASB9zM2UDDHhkIhr1QqgWVZTJDeBwCVe+In96ebPXs21NfXQz6fh0wmQ1RVJYZhoHK5jMrlcmV8x2Ix4Hn+lIDRNE1QFAVM0wSM8Qljn3p7QQgxhBCcTqdnpFKpVimZJMDzDMfzIPA8Cvp0ntRaqBL4qarKEEIkWZYv6erqutC27ZAkSYjneeSvmcFxHDBNE0zT9HRdZzVNY/ziXyOapj2uadpjnuftxRjnAUD1HxohJx1tI4Q0RVF2a5p2FADuqampmclxXLPnefMAgEMIvUYIGcvlcq8RQrKO4+SDHwaO447OmDHjWU3TbtE0TQwKNJVKJVAUpVKt/WT+4+YFQZC3bt36Fv9r/OlosDjNLFmyBAUDd+7cuT3JZHIVz/NLRFGUWJZ1GYbZgzHe/Pvf//7Vnp6ef128ePH/mj17drStrQ01NTWhYOc5WEwHO9mu60I2m0WDg4OWruuuYUy7U27qLBZMpNFotInn+R6GYSQ4fuJd/TGV08KgX9GWLVvgxz/+cWUxEvQoUlUVVFWFcrl8wqZI9dcCAGJZForH49n7779/eM6cOd6hQ4foQvoc09raynR3d8OhQ4dOKHC0atUqxnVdtGPHDuw4zh8dF5IkoWuvvRb94he/wMViEdfU1Aitra0LOjo6ro1EIpfE4/EuSZJigiDwLMuioNS/ZVkn9OoCgMqi2jRNiMfj8K53vQv27dsHtm2DqqrAMAwoigITExNQW1sLoVAIghLypmlW2me4rguKopBisQhjY2NDDz30UHHDhg1v5UtJTQOjo6P8ggULFkQikasXLFiwPBqNtkmSJDIMgwkhhaampn1XXnnl05s2bXpEUZRiEBgihNDSpUuZ7du3e/6f2ZkzZ17Y1NT04Xg8fnE0Gk1JkhT1+yGiIEj0T8WxYRjItm2GYZig3UuwVEHV1SKDYDG4cyaKItTV1ZGamhpcKBSYiYkJyGQykM1mK8VtBEEAQggkEolKkRzLsqBYLKKRkRGYmJh43fM8q1gs0vn77YsghNiWlpZvSpKU4I6fTJ+wUSzLMoTD4cpVlOo+ncF4dBxHiMVivG3bCCGEgtNsf4yCn2rKlkolT1XVLYqiPGTb9lMY4zFCSBmOnyIiQsgfPKX2BzZBCJUIIQWO444CAA8AwbrFBADHdV3bT68NPgcYhsHHjh0rJRKJsWKxmCgUCiSfz6NgrPM8DzU1NZV76EFV60KhAJlMBgYGBqzp2B6GBovTyOLFi1Emk4G2trbOurq6z7a2tr43kUgkZVkWRVFEHMcRlmWvdV3XME1zVzweX51IJEKRSATF43FUXS0vWJyYpgmDg4NE0zQYGxvD+/fv36uq6sSBAwfopEudSQgAgGGYmQihZFXRD3Ryf67qe16PP/54pfS667ogyzIoigL5fB4KhUJlwSGKIiQSiUrASAgBwzAgn8+Tcrlc/Md//Efv61//Orriiium9EWgzpyOjg4kSRIMDQ1hAIBUKsWfd955YiqVYhRFcQ4fPmwWCgUMAHDZZZcxTz755Gl/Q69du5YxTRMDAGlvb5cFQbiovb39g9Fo9OJoNForiqIULLCDyqeSJEEulwO/BHvlTlew6REUB2FZFrq6uqC5uTnYPQdVVcmRI0fIyMgIRKNRJriPa1kWlEoliEajwPN8sIgihUKB6evre6xQKBSPHj2KVq9eTef2s0wqlULr1q2DzZs3N7S0tNyYTCY/GIvF6iVJkgRBYILq0SzLtjEMMwdj/J7169d/ZP369V995ZVXtt54443IT5vzUsetXbhw4cfi8fh5kUgkKgiCyPM8ZhiGrUonJX4qH1IUhSmXy4qiKM8ghPbxPP8B13Vn6rqOMcZMVQGQyrj2g0WCMSYYY891XWKa5nbHcZKZTGbG6OgohEIhJggOTdOERCJRqYKq6zpRVRWOHDlS3LNnzy7bts1LLrmEju23oeBUMZVKreB5fqU/PiopnUHAyPN8pbVQUAHa8zwINiiCvso8zyM/MATTND1N01hN05hyuQzlcnlc07SHVVW9H2O8mxCiAYABAAwA4CAQfCOPOwgoPc9z4XjBGSv4q+BrnBx0+mPZ8DxvqFAo9IyPj5NYLAayLFfu6mqaBtFoFFiWDbIAcbFYRGNjY0M7duwo7dixA73Rx3im0GBxmkAIIY7jIJ1Or25sbPxqOp1eGg6HJVmWid9SICidziOEYhjjNZIkcdU/aNUXZ4MLtTzPg6Io+NChQ3D06NFX7r///ucvu+wy5bbbbptWA5E6e82ePRstW7aMef311zEhRCaE8Cc3cQ52pP1ecyAIAmiaBoVCodIIlxACPM+DqqowPj4OtbW1lV8oGGOoq6urtM8IFtaTk5NoaGiof3Bw0KZNys8tAwMDpKGhQV66dOl5kUjk3UuXLl0lSVKTIAhcOp3W29vbX7v22mufNgzjyZdeemnsH/7hH9D3vve940clCDHr1q2DRx99FD/11FOY5/nYjBkz3p9Op98fjUYvCIfDkiiKPMdxJCh04DhOsMPtaZoGtm0jz/PQyf26qhfV1UGkn2ZHQqEQicfjqFAoQH9/v4cQYlzXRbquQ6FQqOxOY4xd0zSZvr6+/rvvvvve97///erNN99MPvrRj07hq0692cLhMFqxYgV65JFHWlpbW/8tmUxeGQ6HRT9lH8myjILKon6aM++6LtfY2HhxLBabu2bNmv/3nnvuufuxxx6bEYvFruvq6vpALBZrCYVCIZ7nGUEQCEIIYYxZP2WaGIZBdF1nFEWBUqmUVxTlCV3X/8N13V0AYNTU1OwVRfGbhJAZhmEQx3E827aRP64Ry7LErxaJ/GASNE27b2ho6EvJZPI6RVG+NTAwgAVBYACgcqoeVI4EAHAcx1MUhTl48ODjjz76aP8tt9wC3d3ddN3y9oQAAFiWvZphGJ4Q4hFCuOrsz+rMoiBwDMa04zjVmxCIZVkwTRMmJiaIoihsqVTCiqJsUVX1t5ZlPQIAk4QQkxDi+YEqgb+guuifEmD6z6XIMMzjqqpeOjw8XDklrToxh3A4fMqm35EjRx7dsWNH+be//S368pe//Oc+3LcEDRangeCOYk1NTVcikfhyLBZbJQgCkSSJqc7hDkpMu64LmUyGra+vB9M0UdDn6HRYliWapqGtW7eOvPjiiz/p7u7e/vzzz9MJl3rL1dfXM+FwGPr7+zEAuM3NzXWyLPcahhELesoF/z25AAghBARBgN7eXti7d2+lGqSqqsBxHIyOjlZOWE6egDmOA9u2iWEY0NfXp+zcufPZrVu36pdddhl+//vfP9UvC/UWmzVrFpo3bx688sorczs6Oj6fSCTWRaPRmCRJnCAIKFjMchw3hxDyXsMwJq688srv/fjHP779e9/7nvfe976XCXaLm5ub22praz+2dOnSq6PR6CxZlnme5xme54MFNgqKdfgLbCiXy4yqqpUdcf905YRUvWDhE9xhBADAGBPP81zHcXSWZR/zPC+czWbXO45DNE3zisUiqqmpgVAoRDiOQ47jsBMTE8rBgwe/nsvlDkqShMfGxujcfpZZtWoV2rx5c11bW9u3YrHY+/zxx/hFZ4IehshfHxBCCOOnOWOGYZowxt+dM2fOxQihyxKJRFQQBF4QBFKVaoqCVNOgGEi5XIZSqTSsqupvDcP4led5B+C/T1VQPp9/uKamZkCSpM8hhK5yXTfkeV7lRDG4z+ifpL9qWdaPxsfHf0EIsXO53G9SqdTF2Wz2ckKIY1kWo+s6UywWUTQaBVEUMUIIW5bFDA0N9b388st3XXXVVdkf/vCH5JZbbpnKfwrqzxdsws0HAKhuLRRsGlf/Ocgwqv67IGsjqO7vrwOy2Wz2IUVR7sEYvwIAJiHEQQixAID973VGczr9Qj5GNpvdUFtb+/5sNnshx3EuIYQLNkUSiQSEQqHgpNE1TZM5fPjw8IYNG+6+4YYbjO9973uYBovUKb71rW/BJz7xCVYUxc+GQqFl/q4K6+dwo2g0CvF4vJLfv23bNlIqlZCqqqhYLEI+n6+k4wW7cgDHd0P8RfNr+/bt+8revXsf2rVrF+0bQJ3WBRdcwMiyjMLhMMyePZsMDQ1BR0cH6uvrIw888ADxPO8NLUSvueYaZmRkhExMTGCGYVB9fX1XY2Pj1fX19R8SBGGO67p8sVhE1QVAgjLswQl5kKrX29sL73znO2HHjh0QDochHo8Dz/OQyWTIkSNHCACAbdtMMAH7DW+J53luuVxmDx069LOHH374tSVLluC+vj66kD7LdXZ2omXLlsGzzz57fnd397cSicRq/84fCRbXwaabv+PMC4LQHolE/v2v/uqvVvzv//2/P/vggw/m586du1KSpL/u7Oy8OhqNxiVJYnieJ356P/I8rzpIxKqqBkU7iKIorzIM80A4HF4PAL2GYRCMMRv06zpdqh4hBBNCsOM4uFgsbhwdHb3RMIx0Mpk0crnc+3VdZ8fGxiAej4Msy0Ga9Z7BwcGvjY6OPmnbtvHDH/6Qju+zDEKI4Xlebmho+KtQKHStXwSG9VtRoOCOV9D4nud5FGye+cU9AGMc6ezs/ICu64wkSV5Qqt/vO0f8vnOMqqqMoihuqVQ6oKrq3YZh/A4AjkHViUzlkiJCpFAo7ASAG+rq6hZzHHep53kXep4XB4AoQihHCBnEGD80OTn5JBxPAyRw/IRpTJKkz5umyWSz2csMw4BsNgupVAr7wSJyXZfN5/OHDh8+/Hf9/f1btm/f7txxxx1n/PWn3hyEELxo0SIOABLV10+q2wsFhcCC9kIsy4LruoAxBsMwwLKsStVoz/OIIAiora3ty8PDw7/EGKt+P0fGP3yZ6h6FhBAywDDMVx3HuWNsbKzbtm1XVVUoFosomMdZliWu6zLj4+Naf3//Fw8fPty/Zs0aPB3v5tJgcYoFlSCj0ejyVCp1AcuyIZZlGT+1FImiCKFQCGpqakAQBLjjjjvgxRdfROFwGDKZDExOTkIymYRIJAIcx0FTU1NlMaHrOoyNjYHnebtjsdiWNWvW4IULF6Ldu3dPu4FITZ3FixejbDYLg4ODf3QH7pprrmHuv//+035MNBpFK1asQE8//TS+//77cSKRkHt7e5eed955H+E47gpRFBtEUWQkSSKRSARlMhkYGxtDQZXIIJ062I22bRt0XQeO42D16tUwf/58ADheuKlYLMLOnTshk8kgz/OIqqpuoVCo3HdhGIYxDIMbHR198plnnrl98eLF+Ztuumn63Rin3nStra1ow4YNvV1dXd8Ih8OreZ4ngiCwwTwaCoWCanSIZdlgwQyWZaHa2tq/2rlzZ+eSJUsm6uvr3xOJRBDP88i/ywUAwFQVUyD+nSqmXC4z5XJZ0zTt95Zl/YwQ8gQAFARBuI9hmB8DwHLTNJHjOJ5lWSCKIgqCRYRQcELJuK7rapp2/+Tk5GcVRSkhhJRCofCpZDL5U13Xr1RVde74+DiDECq4rruxVCo9aNv2+Mc+9jH4+c9/Tuf0sxMhhNQLgvB3PM8zwakdy7KVhvahUAii0Sj4gVbl/h/Lssg/rYZ4PB58LdZPlyamaWL/rhcqlUquoihbdV2/0zTNpwEg539/BKcpBuKfniAA8CYmJn6PENoO/12wjAH/VMf/nqQqFZD4f3dAFMVP1NbW/k2xWPyAoijdIyMjjCAIwDDMpOM4G8bGxr6rado+mGZ3t6g/j39N6pQg8eTMIn8MgOd5lSsm/uZHpf2Kv9FMIpGIwjCMt2DBAhYA3DN9iviH+GPeyWazz6dSqU85jvONiYmJFYqiwOjoaCXd2r+/eHBkZORfh4eHH3Ucx/rOd74zLcc7DRanmD+oEABc6LcTAIDjQSRCqLL7nEwm4fHHH4dnn322suMyOTkJAwMDleNsx3GgWCxCJBIhCCEwDAOPjo5aL7/88vBHPvKRwtVXXx3cI6AoAADo7e1Fuq7D4OAgmTdvXk9jY+M7RVFcIYpiPcMwhOf5XYSQ3w8NDT330ksvFU/+fFEUmUsvvRQURcEAQMLhcKK7u/uqzs7ODwuCcIF8vDoTSJKERVEkHMchhBDU1NTAxMQECRY11YFisCAPfmkE/FNI4jgOaWtr023bHsxkMmlN01JBWqooiuC6bllRlN/19fV9TVGU/t7e3lM6mVNnH4QQI4piW0tLyz/KsvyOIEPDX1SjYFEdlGUPFiSGYRDLstgdO3Zg0zSXt7a2As/zhOM44u9SB2l9QaopVlWVLZVKSFGUvGEYjzmO81OM8SuEEN1/LGyhUNgfj8ffJ0nS5xFC/4/neQnTNFFQLTJIQSWEYM/zdlqW9Z3x8fFf+5+P/AW5Njo6+gwAPHO653zZZZchGiienRBCLMMwfDqdfi/HcR0IIeyfnFTudp1cECTYKA7uQpmmiXRdJ5IkMblcDlRVrRQEUVWVLRaLuq7rj5bL5Tswxi8FJzJBGt8bqBgJfjXIYI7FhBA3OL3038jpvo5lWaMA8OVwOPzdmpqaGaqqxhzHsXVd77dte/zkdgbU25c/n7kNDQ0lPxWfBJlFwYZx0DqFYZhKgBgUrHMcBwzDAE3TwLIs4rousm1bKxaL/Rhjc+/evdMiSKzmr+3d4eHhTbIsX51IJNYpinJxuVyewzAMYhim33GcjcVi8SHHccYvvfRS9NRTT03bMU+DxWlAkiRREIR2QkgcY8xU52gHOduEEMjn85VfEpZlga7rcOTIkeAXA1FVFTKZDITDYWAYBtu2jY8dO7Z306ZNT/f395s/+9nP4PHHH5/qp0tNE6tXr0bRaBQVi8XZq1at+nRbW9sHY7FYSpZlqKqut9bzPIhEIqM33HDDT375y1/e6nmeuXr1anbz5s2eZVlYEARuwYIF3SzLXjNnzpwP8Tw/T5ZlkCSJSJLk8TzPsizLEL8PkmVZxLIsTAhBhmEwLMsGi/JKAZDqapFVd1+IfxKETdN8XZKkv87n88Q0zbW6rp83OTmJCCF7yuXyJsuyXiOEWAAAe/funbYTMPWXGxkZQf/0T/+EPvnJTwoPPvjgCkEQ3ufPkwzLsqi6f2ckEqmUZQ823URRRIqiQCgUYmpra4kgCAQAGM/zUBAk6rqOdV0HRVEYv1/WkGma97qu+3MAOBBs+gVV//zCCqhYLGYA4Oa6urof8zy/FmP8Ds/zko5zvAc0IWQvxvi5bDb7NCGk0hg6WCj7AWOwGA9OZoKCOfjJJ5+kY/ssxTAMQQhJAHA5QoiAv/4M/r66kmR1Nclg3RAsvv2fAcIwDIyNjYGu66hUKo2oqvoLTdN+QQg5UPU1K+P3jT7O0wWCb7QgiP/9SgCwo/r9y5YtQy0tLWh4eJiO77MDAwAeQmg/IeSdlmURURRPuIZS3W7Fb81SKQpT1cIFTNMklmUhwzDy4+PjxTd6PWYqVB0GTYyNjf0UAH5aW1vLep4HxWKx8jOGEEJPPfXUtAt4q9FgcRqQJAkJgsAEedxB4Y/qHO58Pg89PT2AEALbtgEAQNM0EAQBdu/eDeVyGbLZrFdbWwuyLCNCCCoUCn07d+78oaZpz+/duxdouwwqcP7556NwOMwcPHjwyvb29q/G4/H54XAYRFH0gip7frAWVLRrAoAvf/zjH1971VVXXf/SSy8Ncxwn9fT0XDRv3rwbOI5bK0lSrX8njEiSRHieZxBCbLBLGPRB0jQNBX0SPc8jmqYdv7Tl32VwXbdyfzG40wUAQSDpWJaVGx0d/b+HDx/e5z+d105+fitWrKCniWex6667jrn00kvhk5/8JG5ubia/+tWviCRJ6VQq9R6O42R/oV3Jzjg5YAzS9l3Xrexkp1IpZJpm0ASaBKmmuq6jcrmMSqWSpWnaIcdx7nYc535CyCBAJQsEnbxArj55IYT0A8D/8d9OK1ion/z+07yPzuNnuWuvvZb593//d/bBBx+MHT58uNsvCMIEG2pBIZCTi4E4jlO58119MIcxJtFoFJmmefDYsWM/tizrHkLIJEDlZBCCe7Nn8nkG3y94DP/9bnqqeJYJehBuIIR81nEcZJomCbKKggyiYK0gyzKIoljZLA7Gd1Ct1z9p3D4yMlJYu3bttGszUa1qI5H1/+ghdLzn6djYGIyOjr4txjsNFqcYQgglk0kGIcQEk30QKAY53EGec0tLC/zN3/wNPP/88yAIAixevJhgjGHbtm1w8OBBMjo66vqLINNxnF2FQuE7/f39j031c6SmTl1dHfr2t7+NAI5vLjz33HNk586dEI1G2aNHj67r6Oj4aiQSmcfzPBYEgZEkifOr64EkSZVUOcdxiGmaXjKZXMlx3MaVK1fepWna+yVJWipJUlA4BPv9QBEhpNLQ2TAMbBgG0jQNyuWyUy6XBzRNu1+W5TmCIFxqmqaEMWb9Xl3kpDLZwcki9jwPO45jTE5OPuR53k+uvfZadN9991VSnfyn/AfTnqizw1VXXcU88cQT6IknnhCfeeaZyNKlS0XbtmHdunVzt2/ffnFwN6b692910PiHyrLX1tbCoUOHQNd1YhgGKIqCFEXxyuWyoev6Ntu27yKEbCSEFPyvWVlk/7HHW3W/i4Gq/lw+BvzxSscsdcstt6BvfvOb5OGHH0aPPvoowzCMFI/H036GEaoOEKtatVTudwcVeIP3W5YV3AvDLMuyCxcu/NWBAwe+D/DGUk3PlOnwGKi3jiAI5MILL2Q7Oztf6Ovr24cQWmCaJuE4DgWBYiKRgPb2duB5HnK5HJRKpRP6JmOMK+mrhmF4pVLpQf/+9lQ+tTfEDwZPLhQ17QPEajRYnEIcxyEAQIVCwfY8T49EIp5lWXxwohgcwwc7K4ZhQGdnJ3R0dIBlWZUGzsuWLYO+vr7B/fv33zY5OWlblrXbNM29nueVp/o5UlMjHA6jr33tayifz5ObbrqJ8Xu9kfb2djJ79my0efPm85qbmz8tiuJcjuOI32/rlOp6giAErSmQruucpmnYNM0unue/VVtbG6SaEr/4AhMUp6k6mQFVVd1SqaQpirJb07T/sm37EQAYl2U5gjH+LsMwV9u2HXddF9u2DX57gxMKgJDjlSJziqL8NJvN3qppWqWAwhS+zNQZ9vnPfx49/fTT4U9+8pPvnjNnzl+3t7cvSKfT4VgsRmpra/H27dsTJ5/AnHwSE7yxLHtCur8sy9DY2AgvvfRS0AJjQlXVJ2zb/ikAbK06BfmTF9knLxaqTHXVPmqa+PznP49uu+029kc/+lHkC1/4wqqenp5LyuXyBf/2b/8mO46DgmIg1RvKhmGAqqrAsmylhRBCCBzHAU3TgqbllUIiuq4LCCFuyZIlQAiZ/qts6qzwmc98Bn7yk58ghBBqbGy8U9O0HwVF7IKrVul0Gpqbm6Gmpgby+Tzs3bu3EjD6czjxN4yhXC4/s3///peXL19uP/nkk3Sj4QygweIUaW1tRT09PezBgwejyWTy3bW1tRcRQnh/1wRVn6r4OyqgaRrk83kAgOpSw8SyLBcAtmaz2e9Wf49FixahXbt20cX0Och1XfS1r30ttmTJkoYbbrhhVlNTkwwAxksvvdT385//XBFF8SpRFC/2C20wHMehoLpeOByGWCwG0Wi0ukx/0AgcZTIZiEQiOBQKAcMwLDreRw4MwwC/DDvRNA0URdEVRRlVVXWTrusPuK67lRBSBji+4M5kMgpC6NPJZHKLIAifwhh3W5YlOI4jBOknCCGbEKLbtr3dtu27MpnMo1P80lJTIB6Po1WrVqFnn3028cEPfvDf586d+7EFCxZAU1MTpNNp4Hketm7dCpZloeq73qcry67rOjAMUynDbhhGJZOjrq4OFi5cCLquw/j4+E1btmz5XfAY/pz7XBT1RlxxxRXMs88+y61du/aSZcuW/XNPT8+FPT09wPM8/Od//ic6evQoRCKRUwLFYJ0AcLwIiCAIlTtfQQXJoO2A/+YRQtxLLrmE/R8eEkW9KT7zmc+ge+65J37ttdeuOu+88z4yd+7c83/3u9/Bb3/7W9bviQzpdBqSyWSlWKNf8wBGR0chKMroeR72PA80TRvO5/PfdV138LnnnqPr2zOEBotTgGEYhmVZSdO0BV1dXR+LRqNXJxKJtCzLaGxsjJTL5Uo/LwCo5GoLglBJC8QYE38RRAzDKCmKcnc6nWYBgM1msy4hBNNA8dxz4403orvvvptdtWpV7zve8Y7PzJw5c317e3u6oaEBZFnGy5cvH1m+fPlr3/jGNzowxqK/2KjcGzi5HLssy4AQAp7nASEEuq6jeDwOGGMWAIjjOCRINdV1HTRNcxRFKZbL5aOapm2wLOshjPG+IE/fT90jwYLb/+/Pksnk7ziOW00ImQMAs4PnQwh5HWO8tVAovBxU2Xs75PdTb67rr78efvazn4UvvvjiO2bMmHFdU1MTjsfjJJlMsuFwGHRdB9u2TynJfrqy7EGmRlA8wTTNoEok6LoOsVgMOI4jnudNsiyLFi1axL3yyisOTZWj3grf+c530NNPP43OP//8dTNmzPj5rFmzYi0tLSSRSOBEIoF6enrYAwcOVDY3amtroa2tDQghUCwWQVGUSuZRML4xxtUbJMQ0TUbXdU/TtK0Ax9MCp/p5U2e/d73rXcwrr7ySeve73/313t7eG+bPn8/PmjUL/CsklcyOuro6SCQSEI1Gged5KBQKpFgsgqqqwDBMcA3FNQxjSFGUr42Pj2+xbduja4EzhwaLZ5C/WObD4XBTbW3tZdFo9FOxWGxRPB4noigSlmUhkUiQUqkEqqoiQgjx7zEiy7Iqd2zAvwDueR62bdstFApPuK77zHXXXYdvv/12uut9jmpsbES5XI5ZuXLliiVLltw1d+7cns7OTmhqaoJkMkkEQUCEkFZJkloBjt8DIIRU7nb5J3mVu11+42cAAJAkCQzDqJRon5iYAMdxwC8AQlRVLSuKMqmq6i5d1x+xbftZQsi4/3VPWwAk4J/YFADgYf/ttP5QARDq7Pa5z32Oueuuu0h7e/tf19XVXReNRl3/bi0Kh8PgeR7s2LED7rvvPmBZFmzbBkIIxONxaG5uBowxmKYJmqZV2mWYplkJFv1xHKTsEdM0kaqqZqlUKnmeR2644Qaarke9JRKJBCoWi+Sb3/xmbzKZvKupqSkqy7IXi8W4mpoaJjj1DrS0tMD8+fOhvr4eOI6DY8eOwdGjR8G2bRBF8YTWA0ENBNu2iW3bqFwuHx0dHX0ZAICm7lFvteuvvx7t2LEjunDhwm/39PR8sLm5mUskEl40GmXvuusuyOVyEIvFIJlMQkNDA8TjcQiHwwAAYFkWyWQyuFgsIo7jsOd5RdM0d6mqeluxWNzouq5JA8UziwaLfwGWZdEFF1yAqi/hbtq06bSTcE9PDyeKYkSW5UXpdPoTiURifSwWC4dCIeJfTEeqqhJVVU3XdT3HcSpFPxzHIUErg+qCH67rmoqivJjL5b48MjLivu9976MVIM9R6XQarVu3DjZs2DCvoaHh/7a2ts6MxWJePB5na2pqkCAIyHEcMj4+Tnbv3g25XA7V1tai6vtc1f9ffTITnMRgfHxox+NxcvDgQSiVSkTX9XyhUDhiGMbTlmU96bruDkKICQCVk/H/aVI/qSIegv9u6AxQ1SaABornpu9+97v4jjvuqJEk6R/C4TA5nphxvCUGwzAwOTkJt956K7zwwguQTqchHA7DrFmzYM6cOZBKpQBjDAMDA5DNZisnLkFBEOL39awuy27bNtJ1fWxycrIfAODuu++mixLqLTF//nzEsizT2dl5S21tbVIURSwIAhdU7t28eTM89NBDkEwmYebMmbBs2TJobGwEURTBtu1KY29d18GyrBPSUv1TduI4DrYsC5fL5Xv6+/u1iy++mPlD6xSKejN85jOfYTZs2ECi0eiHGhsbrwma0NfW1rKjo6Po4YcfBoZhQNM06OzshLq6OojFYhAKhcCyLFIul43JycmMoig6xnjQcZynNE37nWmaxwD+5zUF9eajweKfYfny5aijowNSqZRw8ODBsJ/ChzDGbiKRUObNmwcvvfQSTiaTyDRNuPTSS7nR0dFZTU1N74nFYh+KRqPzY7EY8U9tkGmaRFEUUiwWx7PZ7POe5+U4jltpmuYsx3FCgiCc0LCUEOJijEfK5fIzmqb9Wy6XO0bT885ts2fPRnfffTc3d+7czyWTyVmiKHqiKHJBkRrbtmH//v3o3nvvRZs2baoskj3v+EH0SWlL1XcUT6iuF5wuplIprGka5nn+4UKh8M26urq+0dHR6lYB+E8dj38gGKRj+hwWjCVRFHsYhmn1T0uY4M6W67qQz+dh27ZtwLIs1NfXw8qVK2Hu3LkQj8cBIQSWZUE0GoWxsTEwTfOEHp4AULkTXn0Ko2naE4cOHTKuvvpq5oEHHqALa+pNFxTu4nm+heO4d7AsC57nsdXVTvfu3QsAAHPnzoWLLroI2traQJZl8It8QDabhVKpBLZtg23blVYDAFBdEIRVFOXl8fHxOxsaGsimTZvonEq9pe6880782GOPJWRZ/kQoFBL9ugjAsiySJAkSiQRks1lACEFjYyMkk0mIxWIkFAqB4zhOqVTans1m77Isa9iyrAOmaeaCFhR0nTs1aLD4Jzr//PPR0NBQcmBgYObMmTPn1NbWdgiCgERRBNd1zUwmc+jo0aMHGhsbj6TTaXdkZCS9efPm81taWv46Go1eHo/HBVmWCcuyyK9ORkqlkpLNZvcXCoX7VVW9R5KkLM/z7+B5/gpCyDzP8+osy2L9XwJFx3H6bNvemM/n7yeEGPQH6NwWLDpEUWzmef69LMsSQghLCKksHCYmJuC2226DBx98EDiOA0mSgOM4qK2thaamJhBFsVIMIaiuF6TvMQxTOY0JPqa+vp7EYjGcz+dfO3r0aN/ixYtZhJBHjqOLa+pN5Xme4K8VKpUfi8Ui5HI5YBgGFi9eDIZhwNq1a2H27NmVdKagn2c+nwdVVSuBYdAEGgAq1VBd18WO4yBN08ZHR0fv9DyPLFu2bCqfNnV2QwBAIpFIJ8a4JtiwMAwDyuUylEolmDlzJixfvhwuvfRSmDFjBkQiEXBdFzRNg7GxMXLo0CHIZDIoaAnDsiwJEjr8sY5M09w9OTn5rwihTHt7+xQ/ZepsF2zwFYvFnmg02gYAyHVd5GfPQU9PD3z605+G22+/HViWha6uLkgmkxCJREAURYwQKtu2/SzLsr/hOA4Vi0V80lUWagrQYPENmjFjBurq6oLJycn5LS0t19TW1l4eDofnhUKhUNDAHCEEyWQyP2vWrJfGx8cfffnll480Njauj8ViV0aj0c5wOIwEQSCEEGQYBlEUxcvn80P5fP6JYrH4a9u2d2KMDQAAhNBGWZZfEAShlRDSDQA8AAAhZMSyrD7btlW600L5EAAQWZZne54X8zwPu67LWJYF5XIZZFmGwcFBeOaZZ0DXdUgmk9DW1gZz586Fzs5OSCQSwDAMDA8Pw9jYGJTL5cpJoyiKJ1TXMwyjkvJk2zb4aSJQLpffFo1lqbeX+vp6xDAMamxsjHV3d8OcOXPIjBkzoKWlBZLJZKUY09e//nUYHx+v9Af1PA9KpRIplUpw9OhRdPjwYVAUBTiOqxQKq0rpB4wx9jyPWJblFgqFO48cObL/n//5n9H27dvpmKbeUgzDQDQahbq6Oqivr4fa2lqIRCJg2zYsXrwYvvjFLwLLsiBJEgAA5HI5GB4ehl27dkFfXx8AAAnG9Ek1DTTLsrYXCoVvT05OvrRw4UKydetWOp6ptxoCAGBZtsfzPN6vkA66roNhGDA5OQlXXXUVZDIZkGUZVqxYAXV1dRCNRgEAiOd5CkKof+HChQzHccESl47bKUaDxTegsbERzZgxAw4dOrQkkUh8PpFIrItGo6FQKIRCoRAJhUJEFEVgWRYRQpKEkHW2bV/Y0dFxBCG0LJFIgCRJlUW3pmmkWCyWstnslkKh8IBhGBtc1x0J7ngFdF23EUL9ANBf/f4gSAz+/wy+FNQ01NzcjAAAJElKplIpaG5uhoaGBkgkEsCyLOi6DuFwGNasWQP79++Hjo4OWLJkSeXui+M4YJomhEIhcF0XyuUyuK4LpmkCz/OVgglBmp7f7BkpimJqmrbDNE3S09NDTxOpN1V3dzcKh8MkmUx2L1++/OMrV64ks2fPRu3t7dDc3AyRSKRyyrh8+XIwDANGRkYgn89DsVgkmUwGDhw4ALt378b5fB4JggCCICDbtoP+tcQPFgnGmNi27ei6/suRkZE7EELoXe9611S/BNRZbMmSJRCNRvmurq6VF1xwAVq0aBF0d3dDa2sr1NTUVIqLtbS0QD6fh3w+DxMTE5DJZGBwcJAMDQ3ZpVIJeJ4XTNMkVSmoZULIQUVRNtm2/V/lcvlVQoi3c+dOulagzgQCAMCyLLEsi5RKJSiVSqAoChw9ehT27dsHoijC1VdfXUlBDca6ZVnIcRxWkiRRkiTc29uLHn74YTpupwEaLL4BDMPAoUOHZiYSic/F4/ErJEmSeZ4HSZIgGo2iaDSKwuFwUG4dHMdB/f39NaFQaFk0GkXB+/3TRCeXy/XlcrmHFUW5z7btfQBgA5wY+L2BoiD0B4gCAACO43AqlapZsGDB5RdddBH09vairq4uaGxshHA4DAzDQE1NDfzLv/wL7N+/H2zbhkgkUumfaBgGFAoFmJiYAEVRKvcZg4qRVacvQcCIHcdBpVJpz/Dw8P7Zs2ejgwcP0vFIvalYloWhoSFp9uzZ/9TS0nJlKpWCSCSCEolEJVAMIIRAlmWor6+HQqFAhoeHYc+ePc727dvx+Pi4gBACP1jEHMcRf8caARy/s+i67rhhGPfl8/lv5/P5PJ1fqbfSrFmzUH19Pamvr79s4cKFn+rp6UF1dXUolUpBbW1t5RQR4HiadCqVAo7jYGhoiOTzedLa2qp3dna+kM/n93iel3Zdl/E/3NB1vc/zvD26ru/AGCvgnzROyROlzkUEAIBhmH22bduTk5NMQ0MD9PX1AcYYUqkUzJkzBxoaGqCpqQn8tfPxTyQERaPR5FVXXbX2kUce2QIAr0/lE6H+Gw0W/4jLL78cpdNpJpVKCZlM5kOhUGgNx3EhjuMYQRBAlmUUiUQgkUhAPB4HURRRsMiuqakBRVEIz/PBKQwUi8VyLpfbUCgU/ss0zS0Y4wmAPx740Ume+mO+9KUvoaGhIXbBggX/q7e39+ru7m6STqdROp2unCwCHF9wdHR0QH19PQwNDUE2mwW/+i4ZGxtDhw8fRkNDQ5V+no7jVIqABAfefhBJ/Lu2ZqlUun1kZET7+Mc/zgAtREO9idrb25ljx47hRCKxKBaLXROJRCCYdyVJOiVQDO7nhsNhEgqFYGBgQO3r69tYKpUOOo6zCiE03/O8qF8EBPkVIzEAjGKMtxUKhUccx3lS07TMVD1n6txh2zbs3Lkz2dPT8/na2tp2WZaRJEkQiUSQKIqnfHyQqlpfXw+hUAgMwxiwLOun5XL5ac/zop7nBcGiZZpmEWPs0msq1FQIxp0sy4ej0ejhfD6fPHbsGDM5OQkzZsxADQ0NUFtbCzU1NSdv+iEAgFAoFFq+fPm7TNM88sQTT9z65S9/Wf3KV75Cx/AUo8HiHzExMQGGYWDDMDoFQbiY47ga5jhgWRbxPA+yLEM0GoV4PA6yLIPneaCqKuru7obJyUnI5XJEVVW3VCrtz+Vyv1EU5WHHcQ4DAG0oSv1FgovkP/rRj+bX1tbeWFtbK4miiGRZRuFw+IQFdUAURUgmkzA+Pg5jY2Nw+PBh9Oqrr5LBwUHAGCNRFIMTmEoRkKAUe9W9LlQul38+NDT0+LJly5jt27fTFFTqTaXrOgAACIJwJc/zEfALdmCM0emmzSArn2VZEo/HLULIpmPHjv1/mUxmwDTNh2RZni+KYjvGuM113QgAZAghw4qivOp53uuGYRwGv2ULnZept1JnZyczMDCAE4nEClmWF/n1DhBCqJLJUS34M8/zpL6+njAMk33mmWceOnbs2PPhcFgdGRlRqz+eXlOhpgGk63q5oaHhPyzLWjgwMCB3d3dDPB4nPM+jIDPvdGsUAABJkmoWLly4ZuPGjc84jvP8zTffjG699VY6lqcQDRZPsnbtWvTUU0+R6667LrpkyZKGVColbtu2bd3Q0FAnAKCgkTnA8UmcYZhK0QRBECqFQeLxOLS3t0MmkwGO44qJROLbw8PDT8Tj8XImk/Gm+GlSZ4GgUmMoFLpaluUUwzAkWHQEAV4gOH1hGAYkSSKapsG2bdtG9+7dq42OjrY4jiPxPI8dx0GiKFbudfkniwSOL9aRZVnIMIz/Mzw8/G/Nzc06DRSpt0I2mw0WBjMdx2FM00RBkQRd16GmpuaUMe4jkiRpPM9vTyQSu1OpFH7xxRe3Oo6zi+O4MADEAUAAAI0QUrYsi6bpUWdUsVgEAACO4y5GCAlVbVuQZVmVU/KTEAAAnucxxnhXc3Pzz//+7/8+/9BDD8Gdd955ytil45maSoQQjBBCkiQ9EIvF5gLAZ4TjFxODNfTpxnilNzNCiIii2Dlz5sxlsixvXrRoEV1nTDEaLFY5//zzkWEYsfXr11+8bNmya5YuXTqjs7OT43m+/r777qv3PI8JBrrf8LbSDylYXFf3r6upqYE5c+YQXdfJyMiIxnFcKZ1On/oTQlF/hoGBAQAAYBhmsed5yLZtYllWpVJpKBQ63acRlmUJwzCv9vX1/XB4eLifYZj3YYyvs2273vM8sCwLeJ4nLMsiP1hE/r2uIcMw7spms3fruj7GMAxdkFBvKUIIr2kaUhQFyuUyFItFyGazEKT/Vy84CCFgWRaoqkp0XXfmzJnjNTQ0wObNmwkAmAghCwAKJ355uqimzqxisRiMuXbDMFhd15GqquBfVYGamprTzt3+RjSJx+MlVVXHLr/8cnL55Zef2QdPUX8CwzDygiB8P5VKfdB13Xq/XzMxDAMZhlHphXsyQgiDMY7V19c3btiwgf/oRz9qTsHDp6rQYNF3/vnno5GRkbq5c+d+YtWqVR9fsmRJU2dnp9jQ0ADpdBoQQsh1XWTbNpimCYIggGmaoGkalMtlYFkWbNsGjHGlvYB//4vYts25rtuEECLBRV6K+ktlMhkCAOB5XrCgDu7GQqFQgHA4DDzPVz4eIQSe5wV3agfWrFnzQqFQOLpv377XQ6HQZlEU1zmOc6HrunW2bUcQQg5CyCCE7CeEPJvL5R7xPO+Q4zhl2kuReoshACAIoclSqYQnJydRMpmEaDQKoVAIOI4D13UhGo0Cx3GVnqDFYhHGxsbsYrGoyLIMkUik8gX9wJAGh9RUQ3B8HOYVRcH5fB4VCoVKtVNZlqGhoQFkWa5khPjXW0ihUHAVRVHmz59Ps5OoaS3YiItEIiVCiBdsiAQbf7lcDqLR6AlztP954DgOqKrKWJbF9fX10QOWaYBGLgDAcRyKRCJiOp2+pLu7+3+1tLTUpdNpSKVSCAAgm82CrusgiiJEo1Fob28PJm8I+iu6rguiKFZ2tzVNA03TwLIsZFmW43neYKlUIhdeeCEd+NSbDRcKBZLNZiGZTEI8Hq8EirW1tZVeiZ7ngaZpMDk5CSMjI/iVV17xYrEY7u7unjxw4MCjDMNsjsfjKUEQooSQboxxzrbtQrlcniSEFF3XLcLxNBIaKFJnBMuyGwzD+NDw8DAXj8crPRKDsewXFqts0pVKJejv789MTk7uLxQKMDw8PNVPgaJOhgAAGIbZquv6h0dHR9lEIgGhUKiyngiK5AXFbvyNEDw2Nlbev3//Qcdx7Ntvvx3ddNNNdPODmta6u7uN4eHhsUKh0JjP5yGXy6FYLAahUAgEQYCGhoZK1XZCCNi2HWSQuMeOHbMmJiboemMaoMEiAFx55ZVIluW2/v7+D9TU1KQlSSKhUIgRBAHdd999cN999wFCCHp7e6G3txfS6TQQQmDPnj1QKBTAb4tRmdiDvnX+CSNSVVU1DOMgAEBLSwud3Kk3CwMAHsMwW0ql0tqRkRESi8UqF8c9zwNFUSAajQYn30RRFDI2Nma9+uqrw6+//vrk3r17CRw/wbEwxplCoTDp3xvYBceLMGGMMQEAhp7MUGcQQQihxsbGjRjj/ZlMZmFfX5/HMAzreR4KThGj0ShIkgSEEGJZFs7n8/aWLVt2t7a2bl++fDnccsstdLxS002w+H3GcZzc2NhYUzQaJYIgAEIIOY4DiqJAPB6v9Ge2LAvK5TIcO3ZsfHh4+PnXX38dXn/9dbjpppum9IlQ1B/jF+GzGxoadpfL5UWjo6MkFouxwRol2OQLNv2Cdl6lUgmPjIwUdu7cOdDW1uZ84AMfQL/5zW/oXD6FzvlgMR6Po+XLl8MLL7xQN2PGjOUsyyLXdcGyLGSaJmzduhVc14V169bBeeedB7FYDDDGoGkayLIMk5OTgDEGXdcrKX9+Pzpi2zY2TdPLZDKvj42NDXzoQx9C//Vf/0UHPPVmIQAAHMfd6zjOPw4NDYVkWcYsyzLBBkaQjspxHHieR3Rdx8PDw/1bt259rqOjw7juuuvQvffeS04KBAlCyIET73TRtCfqjAkqOo6Ojhaampr+1nXdB4eHh2sxxq5hGKyqqigovS4IAgEArGkaDA4OvrJnz54fAIDW0tIy1U+Dok7hj22GYZjRurq675qmeWtfXx8BAM5xHKJpGgo2QkRRJAgh4jiOl8/n1ZdeeunRXbt27UqlUrBp0ya6lqCmu2CM3uE4zoeHh4e5oFJ7cIpYLperN/3Atm1SKpXgyJEjB0dHR18oFAr0avk0cM4Hi0uXLkW7du2KXnLJJSsXLFhQs2zZMtTc3AzhcBhKpRKcd9550NHRATNnzoRQKAS2bUM+n4fx8XEYGRkBVVXBdV3gef6EBuau64LrutgwDE1V1Z+Njo66a9asof3oqDeNX3GMYRjmUH19/Vd1Xf/24cOHsed5nmEYrKIo4Kc3EY7jCMbYy+fzxr59+56WZXljuVwmDz/88GnHI001paZaEDAihLbW1dVd6zjOr4aGhhpLpRJMTExAMpkk4XAYcRyHHMchuVxu2+Dg4N9hjF+dOXMmPPDAA3SupaalYO6emJj4P/X19fM0Tfv4oUOHQFEUr1gsQjKZRNFoFHieJ4QQVlVVt7+/f0M2m/3BkiVL3Pvvv5+ObWraCzZGEEJ76uvrf2IYxk2HDx92Pc/jbNtGmqZBIpGobPohhIht297ExMT4Cy+8cJ9pmvtTqRQ899xzdLxPsdP2rDpX3HjjjehXv/pV/J3vfOc/XHTRRTcvWbKE7+rqQul0GgWnhJqmwfDwMKiqCrquQyaTgddff528+OKLMDk5iSRJAlEUT2gzQI7n7rmu65JsNvvjgYGBzy1btoxs3LiRLsCpN50/GTMNDQ3fB4AbJUmCdDoN9fX1JJlMBoUSkK7r6uDg4G9fe+21r3R0dAxFIhHYuHHjuTsBUG8L/viGdDrdwrLspzHG72MYZoYkSSzHcTZC6KhhGL8eHR39aSqVGuvq6iIvv/wyHdfUtOdvhvD19fU3IoS+yDBMXSQSgWg0CsEJjGVZuUKhcPvBgwd/4Lpu4X/+qhQ1ffjXWhBCiGtoaLgTIfTXoihCfX2919jYCMlkEkUiEeJnP7HFYrG0b9++nzmO86+NjY36E088QefyaeCcDRbf8573ME8//bTQ3d197Tvf+c6frV69mpkzZw7q7u6GUCiEghNCAABVVeHw4cPQ398P27ZtI5s2bYKxsTGQZdkTRRGJooh4nge/zx0hhLCu63qlUunh8fHxj334wx9Wb7/99nPzhabOGIQQ29jYuJYQ8gWE0GpBEFi/aIKLMd5XLpfvGBoa+g0hREMIgX8XkaKmPeRPrAzDMO973/uEl19+eSbGuE7TtDFN04bgeN9E/OEPfxj98pe/pOOaeltBCDHpdLqR47h3e553IUKoFQDyGOMXSqXS07ZtH6bZHtTbVdA/EQCYhoaGv/U3RtKyLFdSUBmG8SzLOpDNZr9z7NixX2OM7Sl90NQJztk01NmzZ8NTTz2VjMfjNyUSCVYURRwKhZAgCJVxHQSM4XAYEEJk9+7dxDRNBAAFz/PihmGwtm2DYRjBySLyy/6WbNv+96GhoW/QZGvqDCKjo6NPMAzzVGtr6xzDMOqKxSJgjDVFUV4jhKgAEPRNpOOSetsI5lFCCNx///0mALxa/fdLly5lAQBooEi9XU1OTo4AwP9FCP20+v2EELJ27Vpmih4WRf3FqmsfIIS+X19f/xsAuEJRlPZSqdSFEOp3HOfVQqHwLCEkN5WPlTq9czJYDKLBaDTawTDMQgAgrusyGOPTnrQihEg0GoX6+nqEECqvXLnygw888ICGMb7e87xZtm33AkCGEDJECHnGcZxf5/P50TP7rKhzXXAPhhDiAcC+k/++t7eX3bt3r0c3MKi3K3+MI/DbD4BfhOmVV16hBZiot63g1BAhxMJ/1zVA/vswPVWkziJofHx8DAB+erq/vOiii5gXXniBjvdp5pwMFsHPn+Z5vsNxHN7viYgURQHLsiq9jvwPJIQQYFmWDA0NTT711FN/f/DgwRcMwzAB4CUAgLa2tvjQ0JDqL9KDz2PoBE+daVWLDubk9+/du5cuqKm3PdrChTpbVa8hKOpsVLXhd8ppOSHEo4Hi9HSuBosAAMCyLKiqinK5HMnlcmRiYgIlEgngeb7SMxEAQNd1ks/nrf7+/r8fGBh48Ctf+YqTTqeZbDYLcHztUgKoLNARANCdQGpK0fFHURRFUdR042/40Y2Rt5FztsANy7JMfX39Uo7jnp45c2Zk4cKFaNasWaizsxNaW1uhtrYWBEEA0zRxqVSCAwcODA0PD38wHo///iMf+Qgd5BRFURRFURRFndXO2ZNFjDFMTk4eTafTz4+MjFwRjUYJx3GM67qgKArE43EQBAFc1yXlcpkcOHDg2bvvvnvk/PPPJx/5yEem+uFTFEVRFEVRFEW9pc7ZYBGO33kpEkK+oyjKvMOHD7d7nuepqspms1mIx+OE53nsui4eHBwceO211369du3aUdoMl6IoiqIoiqKoc8E5GywSQghCyM1ms1tTqdQ/FAqFbxw4cKA1k8nI6XQaRSIRzDCMUSqVDg8NDf2gr69vqyiK3sjICA0WKYqiKIqiKIo6652zdxYDCCGG4zgukUgsZBjmOp7n5wGAyDCMBgB78vn8fYIg7DMMwzMM49x+sSiKoiiKoiiKOmec88EiwPE+GoQQEolEmEQiIWYyGda2bUeSJNswDLJ69Wq0efNm+kJRFEVRFEVRFHXOoMGiz+/7EpT0BQCAhoYG5DgO5HI5+iJRFEVRFEVRFHVOocEiRVEURVEURVEUdQpmqh8ARVEURVEURVEUNf3QYJGiKIqiKIqiKIo6BQ0WKYqiKIqiKIqiqFPQYJGiKIqiKIqiKIo6BQ0WKYqiKIqiKIqiqFPQYJGiKIqiKIqiKIo6BQ0WKYqiKIqiKIqiqFPQYJGiKIqiKIqiKIo6BQ0WKYqiKIqiKIqiqFPQYJGiKIqiKIqiKIo6BQ0WKYqiKIqiKIqiqFPQYJGiKIqiKIqiKIo6BQ0WKYqiKIqiKIqiqFPQYJGiKIqiKIqiKIo6BQ0WKYqiKIqiKIqiqFPQYJGiKIqiKIqiKIo6BQ0WKYqiKIqiKIqiqFPQYJGiKIqiKIqiKIo6BQ0WKYqiKIqiKIqiqFPQYJGiKIqiKIqiKIo6BQ0WKYqiKIqiKIqiqFPQYJGiKIqiKIqiKIo6xf8PEPpd7HS+ifMAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "from scm.plams import view\n", "\n", "neb_res = job.results.get_neb_results()\n", "mols = neb_res[\"Molecules\"] # these are PLAMS Molecule, not ChemicalSystem\n", "indices = [0, 5, 6, 7, 9] # show only these\n", "fig, axes = plt.subplots(1, len(indices), figsize=(16, 16))\n", "for ax, i in zip(axes, indices):\n", " img = view(mols[i], width=160, height=160, direction=\"tilt_pca3\", guess_bonds=True)\n", " ax.imshow(img)\n", " ax.axis(\"off\")\n", "axes;" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.12" } }, "nbformat": 4, "nbformat_minor": 5 }