Created
November 12, 2016 04:03
-
-
Save amjames/ce93100b611813786a6bc894e385d4a7 to your computer and use it in GitHub Desktop.
Psi4SureIsNeat.ipynb
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
{ | |
"cells": [ | |
{ | |
"metadata": { | |
"collapsed": false, | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "## init cell\n%matplotlib inline\nimport sys\nimport os\npsi4_rel_path = \"stage/usr/local/lib\"\nsys.path.append(psi4_rel_path)", | |
"execution_count": 1, | |
"outputs": [] | |
}, | |
{ | |
"metadata": { | |
"collapsed": true, | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "## this is beautiful\nimport psi4", | |
"execution_count": 2, | |
"outputs": [] | |
}, | |
{ | |
"metadata": { | |
"collapsed": false, | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "# create molecule/set options\nh2 = psi4.core.Molecule.create_molecule_from_string(\"\"\"\n H\n H 1 R\n\"\"\")\npsi4.core.set_global_option(\"BASIS\",\"cc-pvtz\")\npsi4.core.set_global_option(\"d_convergence\", 12)\npsi4.core.set_global_option(\"e_convergence\", 12)\npsi4.core.set_global_option(\"r_convergence\", 12)", | |
"execution_count": 3, | |
"outputs": [ | |
{ | |
"output_type": "execute_result", | |
"data": { | |
"text/plain": "True" | |
}, | |
"execution_count": 3, | |
"metadata": {} | |
} | |
] | |
}, | |
{ | |
"metadata": { | |
"collapsed": false, | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "# do scan\nenergies = []\nrvals = [0.65, 0.7, 0.75, 0.8, 0.85]\nfor r in rvals:\n h2.R = r\n energies.append(psi4.driver.energy('fci',molecule=h2))\n\n# Since h2 is the active molecule it will be used by default in diatomic.anharmonicity\n# However, if you need to provide the routine a molecule pass it as the third parameter:\n# phys_const = diatomic.anharmonicity(rvals, energies, h2)\n", | |
"execution_count": 4, | |
"outputs": [] | |
}, | |
{ | |
"metadata": { | |
"collapsed": false, | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "psi4.diatomic.anharmonicity(rvals,energies,plot_fit='screen',mol=h2)", | |
"execution_count": 5, | |
"outputs": [ | |
{ | |
"output_type": "display_data", | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAioAAAFyCAYAAADbBfG0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3XdYU9cbB/DvCShLUJChIjLdi+KsioILV93VSquirdWq\ntY46aq2KdfZnXa04WhUUnFi3VkXrqFuw1hZXBUHFBQo4QBnn90fITS4JEDDJZbyf58ljzsm597w3\nCDk59wzGOQchhBBCSHEkkzoAQgghhJC8UEOFEEIIIcUWNVQIIYQQUmxRQ4UQQgghxRY1VAghhBBS\nbFFDhRBCCCHFFjVUCCGEEFJsUUOFEEIIIcWWsdQBlFSMscoA/ADcBZAubTSEEEJIiWIKwAXAYc55\nUn4FqaFSdH4AwqQOghBCCCnBPgawOb8C1FApursAEBoairp16xq04gkTJmDp0qUaX3v4EOjRQ5ke\nPBgYP95AgRWBX6gfEl8lAgDs/rTD75t+lzgi/cjvZ1bSldZrK63XBZTeaysu18U50LMnkJAgT9ev\nD2zc+G7n1Oe13Xl2BwN2DBDSw98bjjHNx+ilLoXr16/jk08+AXI+S/NDDZWiSweAunXrwsvLy6AV\nV6xYMd86W7QALlyQPz99GggJARgzUHCF1PJmS+y/tR8AkCZLM/h7aSgF/cxKMn1e25gxwMqVufIO\njMHK7is1H6BD9DMreYrLdUVFKRspADB0KPCuYenz2sKPhQPVlOlJ/Sahnl09vdSlQYFDJ2gwbSnU\nv7/y+d278l+a4sqrivIX7+Xbl3j59qWE0ZDi5t49DXmpGjIJKUbCw8Xpfv2kiUMbnHNsvqa88+JZ\nxdOQjRStUEOlFMr9S5H7l6Y4ea/qe6L01UdXJYqEEELeHefiv7nvvQe4uUkXT0HO3juLuJQ4Ie3f\nwF/CaDSjhkop5Ooq7mYMD5f/8hRHXlXFXZlXHl2RKBJCCHl3//wD3L6tTKv2cBdHG68qB88wMAxq\nOEjCaDSjhkoJNGhQwf+RVH85/vsPuHZNjwG9AycrJ1Q2qyxPNASiHhbj+1TvQJufWUlVWq+ttF4X\nUHqvrThcV+4ebF01VPRxbWkZadj27zYh7ePig+pW1XVez7uiwbQlkDb/Yfv1A6ZPV6Z37AAaNdJj\nUEXEGINXVS8cjTlKDZUSSvXa4uPjkZiYqLNzJyerj7FKjk1GlAEGXtWuXdsg9UihtF5bcbiuTZuU\nzz08gJcvdTNOUB/XduS/I0iJTRHS7Wq102kdtra2qFGjxjufh/Hiek+gmGOMeQGIjIyMLBajzDVp\n3Bj4+2/589q1gevXi+fsn2kR07DozCIAgLHMGC++eQFTY1OJoyKFFR8fj7p16+L169dSh0IIKQbM\nzc1x/fp1jY2VqKgoNGnSBACacM7zbR1Rj0op1r+/sqFy86b8eePG0sakieo4lczsTFx7fA3NHJtJ\nGBEpisTERLx+/VqStYUIIcWLYp2UxMTEd+5VoYZKKTZwIDBzpjK9bVvxbKg0rdZUlL6UcIkaKiWY\nFGsLEUJKLxpMW4rVqgV4eirT27YVz9k/rpVcYWNmI6QvJ1yWMBpCCCHFCTVUSrmBA5XPY2KAyEjp\nYskLY0zUq3Ip4ZKE0RBCCClOqKFSyg0YIE5v3y5NHAVpWlXZUIl+Go1Xb19JGA0hhJDighoqpZyb\nG9BMZbjH9u3F8/aP6piUbJ5NC78RQggBQA2VMkH19k9cnHLDwuIk94BaGqdCSPElk8kwZ84cnZ3v\n3r17MDMzw7lz53R2Tn0ICAiAq6urweqbPXs2ZLLCf0y///77mDp1qtblX716hc8++wxVq1aFTCbD\nxIkTERcXB5lMho3vuu2zDlBDpQzIfftn2zbN5aTkaOmIKhWqCGkap0KKm5CQEMhksjwXxPLx8UGj\n4riqoh4wxsB0uCjTnDlz0LJlS7z//vs6O6c+6Pq6ASAtLQ2BgYE4deqUxvqK0lCZOnUqVq5ciSdP\nnmhVft68edi4cSPGjBmD0NBQDB48WKhf1aFDhxAYGFjoeN4VNVTKACcnoFUrZXrHDiA7W7p4NMk9\noJZ6VEhxlN+HlK4/wIqztLQ0fPvttzo5V2JiIjZu3IgvvvhCJ+craV6/fo3AwECcOHFC7bXvvvuu\nSAso9urVC1ZWVggKCtKq/B9//IGWLVtixowZ8Pf3x3vvvQdnZ2ekpaUJjRYAOHjwoE570rRFDZUy\nQvX2z4MHwJkz0sWSF9UBtbeSbiE5PVnCaAiRVnp6utQh5Kl8+fJF+qavyaZNm1CuXDn06NFDJ+eT\nWlZWFjIyMrQun9/q8DKZDOXLly90DIwx9O/fX+vbNk+ePEGlSpXU8suXLy9qgEu1kj01VMqI/v3F\ny+cXx9s/uRd5K637/pCyY8OGDejQoQMcHBxgamqK+vXrY/Xq1WrlXFxc0LNnTxw5cgTNmjWDmZkZ\n1q5dC0D+YTVu3DiEh4ejfv36MDc3R6tWrfDPP/8AANasWYOaNWvCzMwMvr6+iI+PVzv/jh070LRp\nU5ibm8POzg6DBw9GQkKCqExAQAAsLS2RkJCA3r17w9LSEvb29pg8ebLaB5SmMSoJCQn49NNP4ejo\nCFNTU7i5uWH06NHIzMzM9z3as2cPWrRoAXNzc7XX1q5dCw8PD5ibm6Nly5b4888/4ePjg/bt2wtl\ngoODIZPJ1K775MmTkMlkolsqf/75JwYMGABnZ2eYmpqiRo0amDhxosZG4e7du9GgQQOYmZmhUaNG\n2L17t1oZxTiOJUuWYPny5fDw8ICpqSmuX7+OjIwMzJw5E02bNkWlSpVQoUIFtG3bVtRzEhcXB3t7\nezDGhPEoqu9tXmNUQkND0aJFC1hYWMDGxgbt2rVDRESEqEynTp0QFxeHq1ev5vHOK9+ju3fvYv/+\n/ZDJZDAyMkJ8fLzaGJVhw4YJPTSKOI2MjPI8ty7RyrRlRLVqgLc3oPidDQ8Hli8HDPT/TCuaBtS2\nd22fR2lCpJGSkoKkpCRRHudc47fo1atXo0GDBujVqxeMjY2xb98+jB49Gpxz0a0Oxhhu3LgBf39/\njBw5Ep9//jlq164tvH7q1Cns3bsXY8aMAQDMnz8fPXr0wJQpU7Bq1SqMGTMGz58/x6JFizB8+HDR\nh1ZwcDCGDx+OFi1aYOHChXj8+DGWLVuGs2fP4sqVK7CyshJiyM7Ohp+fH1q2bIkff/wRERERWLJk\nCTw8PDBy5Mg835OHDx+iWbNmSE1NxciRI1G7dm08ePAA4eHheP36tVBHbpmZmbh06RJGjx6t9tq6\ndeswatQotGnTBhMmTEBMTAx69uwJGxsb0ZLs+Y0byZ2/Y8cOpKWlYfTo0ahcuTIuXryIn376CQ8e\nPMA2lW9vR44cQf/+/dGgQQMsXLgQSUlJGDZsGKpX17yz8Pr16/HmzRuMHDkSJiYmsLGxQWpqKtav\nX49Bgwbh888/x4sXL7Bu3Tp06dIFFy9eRKNGjWBnZ4fVq1dj1KhR6Nu3L/r27QsAwlgnTdcWGBiI\nwMBAtG7dGt9//z3Kly+PCxcu4Pjx4+jYsaNQrkmTJuCc48yZM2icx5Lk9erVQ2hoKMaPHw8nJydM\nmjQJAGBnZ6c2vmXUqFFISEhAREQEwsLCDNu7wjmnRxEeALwA8MjISF5SrFzJuXxysvxx7JjUEamr\nsbQGx2xwzAbvv72/1OGQQoiMjOTa/k4kJ3N++rRhHsnJurm+4OBgzhjL99GwYUPRMenp6Wrn6dKl\nC/fw8BDlubi4cJlMxo8ePapWnjHGzczMeHx8vJC3du1azhjj1apV469evRLyp0+fzmUyGY+Li+Oc\nc56RkcEdHBx448aN+Zs3b4RyBw4c4IwxPnv2bCEvICCAy2QyPm/ePFH9Xl5evFmzZmoxBQYGCukh\nQ4ZwY2NjHhUVpf7G5ePOnTucMcZXrlwpylfE3aRJE56RkSHk//rrr5wxxn19fYW84OBg0TUrnDhx\ngstkMn7y5EkhT9PPY+HChdzIyIjfu3dPyPP09OSOjo78xYsXQl5ERARnjHFXV1ch7+7du5wxxitV\nqsSTkpJE583OzhbFzjnnKSkpvEqVKvyzzz4T8hITE9XeT4XZs2dzmUwmpP/77z9uZGTE+/fX7m+j\niYkJHzNmTIHlXFxc+AcffCDKU1xbSEiIkDd27FhRPPkp6O+B4nUAXryAz1vqUSlD+vUDvvxSOZB2\n2zagfTHrsGharSniU+RduDSgtvS6dk3ew2cIp08Dbdro5lyMMQQFBaFmzZpqr02cOBHZuUapm5iY\nCM9TU1ORkZGBtm3b4siRI3jx4gUsLS2F111dXUXfiFV17NgRTk5OQrpFixYAgP79+4tumSjyY2Ji\nUKNGDVy+fBlPnjzBnDlzRGMdunXrhjp16uDAgQOYNWuWqK7cPSfe3t4IDQ3V/IZA/mV3z5496Nmz\nJ9577708y2mi6JmytrYW5Svinjt3LoyNlR9TQ4cOxddff12oOlSp/jxev36NtLQ0vP/++8jOzsaV\nK1dQvXp1PHr0CFevXsX06dNRoUIFoXyHDh1Qr149jYNb+/fvDxsbG1EeY0yInXOO5ORkZGVloWnT\npnnOHCvIrl27wDnHTNVN3PJhbW2NxMTEItVVnFBDpQxxcAB8fYFjx+TpnTuBn38GypWTNi5Vzao1\nw2/XfwMA3E2+i6evnsLOwk7iqAhRatasmcZNF62trdVuCZ05cwazZs3C+fPnRR9wjDGkpKSoNVTy\notpIAYCKFSsCgNqtiIoVK4JzjufPnwOQj4FgjKFWrVpq56xTpw7O5BpVb2pqisqVK6tdl+J8mjx9\n+hSpqamoX79+nmUKwnPdRlDE7eHhIco3NjaGm5tbkeu5d+8evvvuO+zbt090TYqfh6JuAGp1A0Dt\n2rVx5Yr6YpQuLi4a6wsJCcGSJUtw48YN0a3Bol5DTEwMZDKZ1ruTc85LxWw0GkxbxqjO/klKAo4f\nly4WTXKPU4l8WAw3JyJECzExMejYsSOePXuGpUuX4uDBg4iIiMCECRMAQK33xczMLM9z5TVoMa/8\n3B/82jLU4EgFRaMov4ZQQfL6IM7KyhKls7Oz0bFjRxw6dAjffPMN9uzZg4iICISEhIBzrvbzKAxN\nP7vQ0FAMGzYMNWvWxPr163H48GFERESgffv271RXYSQnJ8PW1tYgdelTiepRYYxNB9AdgCeAN5xz\nmwIOAWOsD4BRAJoAsAHgyTn/O1eZEQD8IR93YgmgEuc8VcfhFwt9+wKjRwOKgfjbtgF+ftLGpKpJ\n1Sai9KUHl9DFo4tE0RB9adhQfkvGUHVJYd++fXj79i327dsHR0dHIf+YokvTAJydncE5x82bN+Hj\n4yN67ebNm3B2dn7nOuzs7GBlZSXMQiqMGjVqwMzMDLGxsaJ8Rdy3b98WxZ2ZmYnY2Fh4qmwLr7ht\nlJycLBpke/fuXdE5r127htu3b2PTpk34+OOPhfzcs2UU78nt27fV4r1586bW17Zz5064u7sjPDxc\nlJ/7tk1hejzc3d2RnZ2N6OjoAhcXTEhIwNu3b7XufdGGVL0zJaqhAqAcgO0AzgEYruUxFgBOA9gG\n4Jc8ypgBOJTzWPCOMRZrlSsDHTsCv/8uT+/aBaxeDRRhqr5eWJtZw8PGA/89+w8AcPkhjVMpjSpW\n1N24keJK0Tuh+u05JSUFwcHBBouhadOmsLe3x+rVqzF8+HCUy7nPe+jQIVy/fh2zZ89+5zoYY+jd\nuzfCwsIQFRWl8bZYXoyNjdG0aVNcviz+PW/atKkwI2bYsGHCWI8NGzYgOVm8vpK7uzs45zh16pTw\n4Z2dnS1M71bQ9PMAgGXLlok+gKtUqQJPT0+EhIRg2rRpwu25o0ePIjo6Os/bPLlp6p26cOECzp07\nJ2ogKsYY5b4uTXr37o2pU6dizpw52LFjR74Nh8jISDDG0Ep1tc93ZGFhAUA+3iqvmVz6UKIaKpzz\nQABgjA0txDGhOcc4A9D4U+Wcr8gp004HYRZ7AwcqGyrJycDhw8AHH0gbk6pm1ZoJDZVLD2gpfVJ8\nFOaWSufOnYWFzEaOHIkXL17g119/hYODAx49eqTHKJWMjY2FKctt27bFoEGD8OjRI6xYsQJubm4Y\nP368TuqZP38+jh49irZt2+Lzzz9H3bp1kZCQgPDwcJw5cybfD7VevXphxowZePnypTB41djYGHPn\nzsWoUaPg6+uLgQMHIjY2Fhs2bIC7u7vo+Hr16qFly5aYNm0akpKSYGNjg61bt6o1SOrUqQN3d3dM\nmjQJ9+/fh5WVFXbu3KmxgbBgwQL06NEDrVu3xvDhw5GUlISff/4ZDRo0wMuXL7V6T3r06IHffvsN\nvXv3Rvfu3RETE4M1a9agfv36onOYmpqiXr162LZtG2rWrAkbGxs0aNBA45gfd3d3fPvtt5g7dy68\nvb3Rt29fmJiY4NKlS3B0dMS8efOEskeOHEGNGjVEvU/vSjHl+csvv4Sfnx+MjIwwUHU8gZ7QGJUy\nqHdvcQ/K5s3SxaKJ6jiVhy8fIuFFQj6lCTGcgrq+VV+vVasWdu7cCZlMhsmTJ2Pt2rUYNWoUxo0b\np/G4/NYC0fRafvmqhg4dim3btiEjIwPTpk3DL7/8gn79+uH06dNqDQht1yPJXXe1atVw4cIFfPjh\nh9i8eTO++uorhIaGon379hoXclM1ePBgZGZmYu/evaL8ESNGICgoCA8fPsSUKVPw559/Yt++fRrX\nMtm8eTNat26NRYsWYcGCBejQoQMWLlwoKmNsbIz9+/fjvffew8KFCzFnzhzUrl1b4+qtfn5+2LFj\nB7KzszF9+nTs3r0bwcHBaNKkSYHvhUJAQAAWLFiAv//+G1999RWOHj2KsLAwjedYt24dHB0dMXHi\nRPj7+2Pnzp2i86sKDAzE+vXrkZ6ejhkzZmDWrFmIj49Hhw4dhDKcc/z2228YOlS77/Ta/l/q27cv\nxo0bh8OHD2PIkCHw9/fX6vzvrKD5y8XxAWAogGeFPMYZQDaARvmUaQcgC4CVFucrceuoqOrTR7me\nipkZ5yrLBUju1N1TwloqmA2+6/ouqUMiWijMOiqEqPr000+5t7e3VmV9fHxE66gQdbt27eIWFhb8\n0aNHksVQqtZRYYwtAJDfftQcQF3O+S0DhVQoEyZMEKYKKgwaNAiDBg2SKCLt+PvLx6cAQFoasGcP\noDK+TFJeVb0gYzJkc3nX7YX7F9C7Tm+JoyKE6MusWbNQu3ZtnDt3rtjvoFwS/PDDD/jyyy/h4OAg\ndSgAgC1btmDLli2iPMV0cG1I3lABsBjAhgLKxBgikKJYunRpoQaPFRfduwNWVkBqztymsLDi01Cx\nKG+BhvYNcfWxfI+KCw8uSBwRIUSfnJycirRLMNHs7NmzUocgounLe1RUFJo0aZLHEWKSj1HhnCdx\nzm8V8Mh/V6tCVKej85R4ZmbyqcoKR44AT59KF09uLau3FJ5fSriErOysfEoTQsqS0rCIGdGe5A2V\nwmCMOTHGGkM+3sSIMdY452GhUuYGY6yXSto655j6kM/6qZNzjINKGYecMjVzyjTKKSNe17mUUR0H\nlZUF7NghXSy5tXBsITx/+fYlop9GSxgNIaS4+OOPPwy6Fg2RXolqqACYAyAKwCwAFXKeR0G+mJtC\nTQCqg0Z6ArgCYB/kPSpbco5R3dBiVE6ZNTllTuaUKUaTdnWvfXv5svoKYWHSxZKbao8KAJy/f16i\nSAghhEipRDVUOOfDOOdGGh6nVMoYcc43qqRDOOcyDcfMUSkTmEcZ9XlrpYiREfDRR8r02bNArgUi\nJVPbtjYqmijbmzROhRBCyqYS1VAhupd7AO3WrdLEkZuMydDcsbmQph4VQggpm6ihUsY1bQqobhIa\nFiZfXaU4UB2nEv00GqlvSuX2S4QQQvJBDZUyjjHxoNp//wWuXZMuHlWq41Q4OC2nTwghZRA1VAhy\nr4JcXJbUV731A9A4FUIIKYuooUJQuzaguu7O5s1Arv28JGFnYQd3a+UGZDROhRBCyh5qqBAA4l6V\ne/eAM2eki0VVi+rKcSoXHlwo1O61hJQEAQEBsLS0lDoMycTFxUEmk2ncHLCoLl68CBMTE9y7d09n\n5zS0kydPQiaT4dSpUwUXVvHNN9+gZcuWBRdU8b///Q/u7u4wNjYWVlp3cXHB8OHDC3UefaGGCgEg\nn6asuthjcbn909JR+Qv35NUT3E2+K10whOQICgqCTCbTyb40+e2cXFbo+vpnzJiBjz/+GE5OTjo9\nrz6sWrUKISEhGl8ryvsyfvx4XL16Ffv379eq/JEjRzB16lR4e3sjODgY8+fPBwDIZDJR/devX0dg\nYCDi4+MLHdO7Kg57/ZBioFo1wNcXOH5cnt6+HVi+HChfXtq4VHtUAHmviqu1q0TRECK3efNmuLq6\n4uLFi4iJiYGbm5vUIZVYzs7OSEtLQ7ly5XRyvr/++gsRERE4f75k3CoOCgqCnZ0dhg4dKspv164d\n0tLSUL6Qf4QdHBzQq1cvLF68GD169Ciw/B9//AEjIyOsW7cORkZGQv7Nmzchkyn7MqKjoxEYGAhf\nX1/UqFGjUDG9K+pRIQLV2z/Pnsn3/5GaZxVPmBiZCGkap0KkFhsbi7Nnz2LJkiWwtbVFWHFa0jkP\nnHO8efNG6jDyVL58eZ31qmzYsAHOzs5o3rx5wYUllJaWVmCZwjZSFAYMGIA///wTd+/eLbDs48eP\nYWZmJmqkAEC5cuVEeZxzyXr+qKFCBP36iXtQisPf3/JG5fFe1feENM38IVILCwuDjY0Nunfvjv79\n+2tsqCjGXSxZsgS//PILPDw8YGpqiubNm+Py5csaz5uQkIDevXvD0tIS9vb2mDx5stqYrMWLF6N1\n69awtbWFubk5mjZtip07d6qdSyaTYdy4cdi8eTMaNGgAU1NTHD58WBRXUFAQ3N3dYWFhAT8/Pzx4\n8AAA8P3338PJyQnm5ubo3bs3kpOT1c4fFBQknNfR0RFjx45FSkqKqIyPjw8aNWqE69evw9fXFxYW\nFqhevTr+97//aXyvco9RuXnzJgYMGAB7e3uYm5ujTp06mDFjhsb3TtWePXvQvn17ja/NnTsXTk5O\nsLCwQIcOHRAdHQ1XV1fRWIzZs2eLehIUgoODIZPJRLc+9u7dix49esDR0RGmpqbw8PDA3LlzkZ1r\nNoLivYiKikLbtm1hYWGB6dOnw9XVFf/++y9OnDgBmUwGmUwmxJ7XGJULFy6gW7dusLGxQYUKFdC4\ncWOsWLFCVKZjx47gnGPPnj35vlcymQwhISF49eoVZDIZjIyMhJ+D6hiVkJAQDBgwQLgWRdnCjp8p\nKrr1QwSVKgHduwO7dsnTu3cDqamAlZW0cbV0bCn0pFx5eAVvMt/AxNikgKMI0Y/NmzejX79+MDY2\nxqBBg7B69WpERkZq3LI+LCwML1++xKhRo8AYw6JFi9CvXz/ExMSIvq1mZmbCz88PLVu2xI8//oiI\niAgsWbIEHh4eGDlSuS3ZihUr0KtXL3zyySd4+/Yttm7digEDBmD//v3o2rWrqO5jx45h+/btGDt2\nLGxtbeHi4iK8FhoaioyMDIwbNw7Pnj3DokWL8OGHH6J9+/Y4efIkpk2bhv/++w8rVqzA119/jV9/\n/VU4dvbs2ZgzZw46d+6M0aNH4+bNmwgKCsLly5dx5swZ4boYY3j27Bm6du2Kvn374qOPPkJ4eDim\nTZuGRo0awc/PL8/3+O+//4a3tzdMTEwwcuRIODs7486dO9i/fz/mzp2b53EJCQmIj48XBoSq+u67\n7zBv3jz06NEDXbt2RVRUFDp37oyMjAxRubzGDGnKDw4OhqWlJSZNmoQKFSrg+PHjmDlzJl68eIFF\nixaJjk1MTES3bt3w0UcfYciQIXBwcICvry/Gjh0LS0tLzJgxA5xzOKhswJa7vqNHj+KDDz5AtWrV\nMH78eFSpUgXXr1/HgQMHMG7cOKGclZUV3N3dcebMGXz11Vd5vl+hoaFYs2YNLl26hHXr1oFzjlat\nWqnV3bZtW4wbNw4//fQTZsyYgTp16gAA6tatm+e5dYpzTo8iPAB4AeCRkZG8NPntN87la9PKH+vW\nSR0R51uvbeWYDeFx4f4FqUMiGkRGRnJtfyeS05L56bjTBnkkpyXr7BovX77MGWP8+PHjQp6TkxOf\nMGGCqNzdu3c5Y4zb2dnxlJQUIX/v3r1cJpPxAwcOCHkBAQFcJpPxefPmic7h5eXFmzVrJspLT08X\npTMzM3nDhg15x44dRfmMMW5sbMxv3LihMS4HBwf+4sULIX/69OmcMcbfe+89npWVJeT7+/tzU1NT\n/vbtW84550+fPuUmJia8a9euovOuXLmSy2QyHhwcLOT5+PhwmUzGw8LChLy3b9/yqlWr8g8//FAt\nppCQECGvbdu2vGLFivz+/fu8MI4dO8YZY6L3VzXunj17ivK//fZbzhjjw4YNE/Jmz57NZTKZ2rmD\ng4O5TCbjcXFxQl7unwfnnI8aNYpXqFBBeM84V74Xv/zyi1r5Bg0acF9fX7X8EydOcJlMxk+ePMk5\n5zwrK4u7urpyNzc3npqamtdbIPDz8+P169cvsFxAQAC3tLRUy3dxcRG9L+Hh4aJ4ClLQ3wPF6wC8\neAGft9SjQkS6dQNsbORjVABg0yZA6hlquQfUnr9/Xm0xOFKyXHtyDd4bvA1S1+lhp9GmRhudnCss\nLAxVqlSBj4+PkDdw4ECEhYXhxx9/VPsG/NFHH8FKpUvS29sbnHPExMSonVu150RRNjQ0VJRnYqLs\nSUxOTkZmZia8vb2xVcMmXT4+Pqhdu7bG6xgwYAAqVKggpFu0kP+ODR48WHTbo0WLFti6dSsePHgA\nFxcXREREICMjA+PHjxedb8SIEZg+fToOHDggGhRaoUIF+KsMfitXrhyaN2+u8foVEhMTcfr0aUyY\nMAGOjo55ltMkKSkJjDFYW1uL8hVxf/nll6L88ePHC7NcikL15/Hy5Uu8efMGbdq0wdq1a3Hjxg00\nbNhQVDYgIKDIdV25cgV3797F8uXLtZrObm1tjb/++qvI9RUnNEaFiJiYiHdUPnECiIuTLBwAgHNF\nZ1SpUEX3oPXKAAAgAElEQVRIn713VsJoSFmVnZ2Nbdu2wdfXFzExMbhz5w7u3LmD5s2b49GjRzh2\n7JjaMbmnx1aqVAkA8Pz5c1G+qakpKleuLMqztrZWK7d//368//77MDMzg42NDezt7bFq1Sq18SEA\nRLd6CoqrYkX5TuXVq1fXmK+IIy7nj0GtWrVE5cqVKwc3NzfhdYXc58vrulQpGjH169fPs0xBeK6x\nPYq4PFQ3NgNga2ur1qgpjOjoaPTp0weVKlWClZUV7OzsMHjwYABQ+5k4OjrC2LjofQN37twBY0zr\n94VLOPhV16ihQtQMGSJO5/pSZ3CMMbRyaiWkz9wrJqvRkTLl+PHjePjwIbZu3YqaNWsKj4EDB4Ix\npnFQbe6ZFAq5P0jzKqfq9OnT6NWrF8zNzbFq1SocOnQIERER8Pf317gQopmZWZ7nyqs+bePVlq7P\nV5DKlSuDc55vQ6ggeX24Z2VlidIpKSlo27Ytrl27hrlz52L//v2IiIgQxqbkHlCb389DH54/fw5b\nW1uD1qkvdOuHqGneHKhVC7h1S57euBGYPl28IJyhtXZqjd+u/wYAuJ96H/Ep8ahR0bBz+YnuNLRv\niNPDThusLl0IDQ2Fg4MDgoKC1D5od+7ciV27dmH16tWi2wG69Ntvv8HMzAyHDx8WfTNft26dXurT\nxNnZGYB8Ro5qj01GRgZiY2PRqVOnd65DsSbNP//8U+hjFYM8Y2NjRfmKuG/fvi2KOzExUa1Ro+hh\nSU1NFd22yz3V98SJE3j+/Dn27NmD1q1bC/l37twpVMza9nq4u7uDc45//vknz1lNqmJjY+Hp6Vmo\nWPIjZe8MNVSIGsbkvSqKmYC3bgEXLwItWuR/nD6p9qgAwJn4M6jRkBoqJVVF04o6GzdiCOnp6di1\naxcGDhyIPn36qL1etWpVbNmyBXv37sWHH36olxiMjIzAGENmZqbQULl7926BU1B1qWPHjihXrhxW\nrFghmrXz66+/IjU1VasFxgpia2uLtm3bYv369ZgwYUKhVpetVq0anJyc1KaAd+zYEcbGxvjpp59E\njamlS5eqnUPRIDh16pRwPa9evVKbPm1kZATOuajn5O3btwgKCtI6XgCwsLDQOAU8Ny8vL7i6umLZ\nsmUYOnSocFtOk9TUVNy5cwdjxowpVCwFxck51ypWXaOGCtHok0+UDRVA3qsiZUPFq6oXTI1NkZ6Z\nDkA+TmVQw0HSBUTKlD179uDFixfo2bOnxtdbtmwJOzs7hIWF6a2h0r17dyxZsgR+fn7w9/fH48eP\nERQUhJo1a+Lvv//WS52A+DaNra0tvvnmG8yZMwddunRBz549cePGDaxatQrNmzfHxx9/rJM6V6xY\nAW9vb3h5eeHzzz+Hq6srYmNjcfDgQVy5ciXfY3v16oXdu3eL8mxtbfH1119j4cKF6NGjB7p164Yr\nV67g999/h52dnahs586dUaNGDQwfPhyTJ0+GTCbDhg0bYG9vL9o7qFWrVrC2tsaQIUOEqcGhoaGF\n7nlo0qQJVq9ejXnz5sHDwwP29vbw9fUFIH7vGWNYtWoVevbsCU9PTwwbNgxVq1bFjRs3EB0djUOH\nDglljx49CgB5/n8tCk9PTxgZGWHRokVITk6GiYkJOnToYJDbSzRGhWjk7Ay0a6dMb90KSLmwZXmj\n8mhWrZmQpnEqxJA2b94Mc3NzdOzYUePrjDF0794dv//+u3AroTDrceT14aaa7+vri/Xr1+Px48eY\nMGECtm3bhh9++AG9e/fWqo6CXtMmBgCYNWsWfv75Z9y7dw8TJ05EeHg4Ro0ahcOHD6uNSdH2nLnT\njRo1wvnz59GuXTusXr0aX331FXbt2qXxWnMbPnw4Hjx4gLNnxYPu582bh8DAQPz111+YMmUKYmNj\nceTIEZibm4vKGRsbY/fu3fDw8MDMmTPx888/4/PPP1frnbCxscGBAwdQrVo1fPfdd0Ij8ocfftDq\nmhVmzpyJbt264X//+x/8/f3x/fff53lM586d8ccff6B27dpYsmQJJk2ahOPHj6s1SMLDw9GmTRu4\numq33Yg2/08dHBywZs0aPHnyBJ999hn8/f0RHR2t1fnfWUHzl+lRttZRUbVunXhNld9+kzaeaUen\nCWupyAJlPDW94LUEiOEUZh0VQvSpQ4cOfPDgwVqVzb1eSEn38OFDbmZmxvft2ydpHLpcR4V6VEie\n+vcHTE2VaR3uwl4kquNUsnk2LadPCNFo/vz52L59u+hWTVmxfPlyNG7cWCfjhYoLaqiQPFlZAarj\nBg8cABITpYtH04BaQgjJrXnz5khPTy/UQNzSYsGCBTh37pzUYegUNVRIvlTXVMnIALZtky6WyuaV\nUce2jpA+e58WfiOEvJv8xvOQ4oEaKiRfHTsCVZSLwkp++6e1k3K9gnP3ziErOyuf0oQQkr+YmBiD\nrkVDCo8aKiRfxsaA6ozDixeBmzeli0f19s+Lty/wz5PCLwpFCCGk5KCGCilQ7iX1N22SJg5A3KMC\n0DRlQggp7aihQgrUqBHQuLEyHRICZEl0x6VW5VqwNVcuMEQNFUIIKd2ooUK0otqrcv8+oGGjWIPI\nvUEh7aRMCCGlW4lqqDDGpjPGzjDGXjHGnml5TB/G2GHGWCJjLJsx1ijX69aMsRWMsRuMsdeMsTjG\n2HLGmFVe5yyLPvlEPl5FYcMG6WJpVV3ZULmbfBcJLxKkC4YQQohelaiGCoByALYDWFWIYywAnAYw\nBfJV8HKrBqAqgIkA6gMYCqALgF/fKdJSxt4e+OADZXrXLuCZVk1F3WtdI9c4FVpPhRBCSq0S1VDh\nnAdyzpcDuFaIY0I553MBHAOgNlmec/4v5/xDzvlBznks5/wEgG8BfMAYK1Hvj74NG6Z8/uYNsGWL\nNHE0rdYU5Y3KC2kap0IIIaUXfRBrVglAKuc8u8CSZUjXruI1VaS6/WNqbIomVZsI6dPxp6UJhBBC\niN5RQyUXxpgtgBkA1kgdS3FjbAwMHqxMR0YCV69KE0tb57bC878e/YWU9BRpAiGEEKJXkjdUGGML\ncga55vXIYozVMlAslgAOAPgHQKA2x0yYMAE9e/YUPbZIdU/EAFRv/wDS9aqoNlSyeTbN/iF6FxIS\nAplMhqioKKlDAQCkpaUhMDAQp06d0qr8yZMnIZPJhIepqSmqVKkCX19fLFiwAInvsJHX9evXERgY\niPj4+CKfg5ReW7ZsUfucnDBhgtbHGxdcRO8WAyjo4y5G30EwxioAOAwgGUBfzrlWK4UsXboUXl5e\neo2tOKlbF3j/fUCx51VoKPDDD0D58vkfp2utnVpDxmTIzrk7dzLuJLrW7GrYIIjecc6xb98+bNiw\nATdu3ECFChXg5+eH0aNHo1q1agaPpzjtCfP69WsEBgaCMYa2bdsWfECO8ePHo2nTpsjKysLTp09x\n9uxZzJ49G0uWLMH27dvh6+tb6Fiio6MRGBgIX19f1KhRo9DHk9Jt0KBBGDRokCgvKioKTZo0yeMI\nMckbKpzzJABJhqpOU2ZOT8phAGkAenLO3xoonhJp2DBlQyUpCdi3D+jXz7AxVDStCM8qnoh6KP92\neypOu2+VpHiIjY3F2bNnkZGRgSZNmqBhw4ZqZdLT0zFw4EDs3btXlH/58mUsX74cO3fuROfOnQ0V\ncrHDucY/ZwVq06YN+vbtK6QnTpyIa9euoVOnTujfvz+io6Ph4OBQ6FiKUyOOlC6S3/opDMaYE2Os\nMQBnAEaMscY5DwuVMjcYY71U0tY5x9SHfNZPnZxjHHJetwRwFIA5gM8AVGKMOeQ8StT7YygDBwJm\nZsq0ZLd/aii/RV5KuITXGa+lCYRo7dGjR+jTpw/c3NzwySefYNiwYWjUqBFat26N6OhoUdmvv/5a\nrZGi8PLlS/Tp0wd3794V5XPOcezYMYwdOxYfffQRJk+ejCtXrujrchAQEABLS0skJCSgd+/esLS0\nhL29PSZPnixqSMTFxUEmk2HJkiVYtmwZXFxcYG5uDh8fH/z777+ic/r4+KB9+/Ya63J1dRXOZ29v\nD8YYZs+eLdzOmTNnTpGuo2HDhli2bBmeP3+On3/+WciPj4/H6NGjUadOHZibm8PW1hYDBgxAXFyc\nUCYkJAQDBgwQYpfJZDAyMhJuSe3duxc9evSAo6MjTE1N4eHhgblz5yI7m+YqEO2UtA/iOQCiAMwC\nUCHneRQA1f6jmgAqqqR7ArgCYB/kPSpbco4ZmfO6F4BmABoC+A9AAoCHOf9W19N1lGhWVsCHHyrT\nhw4BCRKsuaY6TiUzOxPn7583fBBEa8nJyWjXrh12796t9trZs2fRpk0b3Lp1CwDw+PFjrF27Nt/z\nvX79GsuWLRPSiYmJ8PX1RceOHbFy5Ups27YNixcvhpeXFwYPHow3b97o9oIgvxWUnZ0NPz8/2NnZ\n4ccff4SPjw+WLFmiMf6QkBD89NNPGDt2LKZPn45///0XHTp0wNOnT0XnzKsuxWt2dnZYvXo1OOfo\n27cvQkNDERoaKuopKaz+/fvDzMwMR44cEfIuXbqE8+fPY9CgQfjpp5/wxRdf4NixY/D19UV6ejoA\noF27dhg3bhwAYMaMGQgNDcWmTZtQt25dAEBwcDAsLS0xadIkrFixAk2bNsXMmTPxzTffFDlWUsZw\nzulRhAfkDRweGRnJy6I//uAcUD4WLDB8DE9ePuGYDeEx8/hMwwdBBJGRkTy/34kpU6ZwyL8s5Pno\n3r0755zzjRs3FlgWAHd3d+ecc56ZmclbtmyZb9lhw4a90/UFBwdzmUwmur6AgAAuk8n4vHnzRGW9\nvLx4s2bNhPTdu3c5Y4xbWFjwhw8fCvkXL17kjDE+adIkIc/Hx4f7+vqq1R8QEMBdXV2FdGJiImeM\n8cDAQK3iP3HiBGeM8Z07d+ZZxtPTk1euXFlIp6enq5W5cOECZ4zx0NBQIS88PJzLZDJ+8uRJtfKa\nzjFq1CheoUIF/vbtW61iJyVPQX8PFK8D8OIFfN6WtB4VUky0bQu4uSnTGzbImyyGZGdhh3p29YT0\nqXgap1JcZWdnY/369QWWO3jwIB4+fIjk5GStzvv8+XMAwO7du3H+fP49aooBufowcuRIUdrb2xsx\nMepzAPr06YMqKosRNWvWDC1atMDBgwf1EldhVahQAS9evBDSJiYmwvPMzEw8e/YMbm5uqFSpktaz\nn1TP8fLlSyQlJaFNmzZ4/fq13n4epHSRfDAtKZlkMiAgAJg5U56+dUveWPH0NGwc9czbIRrysQ3n\n4s/j/KU3KG9kUsBRRB+uX8/7tadPn2o1/ZVzjuvXr8NNtRWcD3d3dwDA9u3btSq/fft2zFT8p9UR\nU1NTVK5cWZRnbW0tNKJUeXh4qOXVqlULO3bs0GlMRfXy5UtYWloK6fT0dMyfPx/BwcF48OCBMO6G\nMYaUFO3WLoqOjsa3336LP/74A6mpqUJ+Yc5ByjZqqJAiGzoUmDVL2ZPy6acSBNGgLdBfvvXTm+x0\nvN//EhDfRoJASH7Mzc21LmthYYE2bdrAwcEBjx8/zrfs4JwVCBO0HCSlbbnCMDIy0un58hqjkpWl\n1YoJRZaZmYlbt26JZmCNHTsWISEhmDBhAlq2bImKFSuCMYaBAwdqNRg2JSUFbdu2RaVKlTB37ly4\nubnB1NQUkZGRmDZtGg2oJVqhhgopsho1gE6dAMXYOzMz4PBhwMIi/+N06UmaN7pGKNNjfjiF4TWp\noSKF69flu2xrYmlpiRYtWuDChQv5nqNy5crw9PRE+fLlsXTpUvj7++dZtnHjxhgxYgQAiG6n5Efb\ncvpy+/Zttbxbt27BxcVFSFtbWyM2NlatnOpMG0D3a7rs2LEDaWlp6NKli5C3c+dOBAQE4IcffhDy\n3rx5o3ZrLq9YTpw4gefPn2PPnj1o3Vq5meidO3d0Gjsp3aihQt7Jp58qGyppafIPq88/N2QEjnCP\ndMed5/I/fP9lnIKX13RDBkC0NGXKFPQrYMGdr776ShjTMGjQIGRmZmLcuHFqH4xdunTBxo0bYWpq\nCkA+YyU8PLzAGPr371/E6HVj9+7dSEhIEBaru3jxIi5cuICJEycKZdzd3XHo0CEkJSUJt5SuXr2K\nM2fOiBZTU/RSaTueJz9Xr17F+PHjUblyZYwePVrINzIyUuv1WLFihVrvjoWFBTjnarEYGRmBcy46\nx9u3bxEUFPTOMZOygxoq5J306gXY2gKK4Qdr1xq6oQK0c24nNFTO3DuDzOxMGMvov3Zx07dvX0yf\nPh3z58/X+Hq/fv3UpqwOHjwYffv2xZ49e3Djxg1YWlqic+fOaNy4sdq5vby88h3gOWjQIDRo0OCd\nroG/44hxDw8PtGnTBl988QXS09OxfPly2NnZYfLkyUKZ4cOHY8mSJejcuTM+/fRTPH78GGvWrEGD\nBg1EYzxMTU1Rr149bNu2DTVr1oSNjQ0aNGiA+vXr5xvDqVOnkJaWhqysLCQlJeHMmTPYu3cvrK2t\nsWvXLtjb2wtle/TogU2bNsHKygr16tXDuXPncOzYMdja2orO6enpCSMjIyxatAjJyckwMTFBhw4d\n0KpVK1hbW2PIkCHCFObQ0FBaHI4UTkHTguhB05ML8vXX4qnKhn5Lgq8Ei6YpX7x/0bABEM55wdMR\nFY4ePcp79OjBK1asyM3Nzbm3tzffvHkzz8rKeqf6ExISePPmzTVOTe7bty9/9erVO50/r+nJVlZW\namVnz57NjYyMhLRievKPP/7Ily5dyp2dnbmZmRn38fHh165dUzt+8+bN3MPDg5uamnIvLy9+9OhR\nHhAQwN3c3ETlzp8/z5s1a8ZNTU25TCbLd6ryiRMnuEwmEx4mJibcwcGB+/j48IULF/LExES1Y1JS\nUvinn37K7e3tuZWVFe/WrRu/desWd3V15cOHDxeVXbduHffw8ODlypUTTVU+d+4cb9WqFbewsODV\nq1fn33zzDT969Gie05lJ6aDL6cmMG3pOaSnBGPMCEBkZGVmm9vrR5NYtoHZtZXrUKGDVKsPVH/s8\nFm4rlLNEFndajEmtJhkuAAJAuXeHlL8TWVlZOHjwIHbt2oWnT5/CyckJ/v7+aN26taTf4uPi4uDq\n6orFixeLbvMQUloV9PdAZa+fJpzzfOe6U/84eWe1agHt2gEnT8rTYWHA//4HVKhgmPpdKrmgulV1\n3E+9D0C+ngo1VMomIyMjfPDBB/jggw+kDoUQoiO04BvRCdVxKS9eANu2Ga5uxhjaObcT0qfiTiEr\nW79TOQkhhBgGNVSITvTtC9jYKNO//GLY+n1cfITnyenJ+OvRX4YNgJACqO7VQwjRHjVUiE6YmgJD\nhijTFy4AV68arv72ruLdZo/HHjdc5YQUwNnZGVlZWZgwYYLUoRBS4lBDhehMztpbAkP2qrhZu8Gl\nkouQPhZ7zHCVE0II0RtqqBCdqVcPaKOyKGxoKPD6teHqb++i7FU5HX8ab7PeGq5yQgghekENFaJT\nqr0qKSmAIfdaU7398zrjNS4+uGi4ygkhhOgFNVSITn34IVCpkjK9dq3h6vZ19RWlaZwKIYSUfLSO\nCtEpMzNg8GDgp5/k6bNngX//BQpY1VsnqllWQx3bOriReAOAvKEys91M/VdMRK5fvy51CIQQieny\n7wA1VIjOjRihbKgA8kG1y5YZpu72Lu2FhsrZe2dx5t4ZmBmbGabyMu5h1kOYmZvhk7y2UCaElCnm\n5uZq+0IVBTVUiM41bAi0bAmcPy9Ph4QA8+cDOZu96lV71/YIuizfmTUjOwNt1rcp4AiiU58D0OUA\n6j+nAm0W6fCEpNSLawMcXq5MN18OeG6ULp48THx/Ij5u9LHUYeiVra2taMfvoqKGCtGLzz9XNlSS\nk+Ur1Q4bpv96fVx8wMDAId/DapjnMIxtPlb/FRM9GiB1AKQEGTfYHWdQEQBQrnw2Dq2bBOvKX0ka\n072Ue+i3vR+yuHzF7Pp29THHfw4syltIGldJQQ0VohcDBwKTJgHPn8vTQUGGaahUNq8MzyqeuPLo\nCgAg+mk0vKqW7U0jCSkrYmOBsypj6AcOkKFDg0bSBZRj6fmlQiMFAJb4LaFGSiHQrB+iF+bm4obJ\n5cvApUuGqVt1mvKlhEtISU8xTMWEEEmtWQNwrkyPHi1dLArRT6MR9neYkPau4Y1Obp0kjKjkoYYK\n0ZtRo8TpoCDD1KvaUMnm2Tgdf9owFRNCJPPmDbBunTLduLF8rJzUZp2YJdyKBoC57efSnk+FRA0V\nojc1awKdOyvTW7cCSUn6r9e7hjeMmJGQpvVUCCn9wsOBxERlevRoQOr2wJWHVxAeHS6kO7l1Qlvn\nthJGVDJRQ4XolWrXa3o6EBys/zotTSzR3LG5kKaGCiGln2qPraUl4O8vXSwK3xz7RpSe236uRJGU\nbNRQIXrVvTvg5KRMr1oFZGfrv17V2z9XH1/Fk1dP9F8pIUQSV6/KF5dUGDoUqFBBungA4FjMMRy+\nc1hI96rdS/QFimiPGipEr4yNgZEjlek7d4CjR/Vfb0e3jqJ0REyE/islhEhi1SpxOvf4OEPL5tmY\nGjFVSMuYDPM7zJcwopKtRDVUGGPTGWNnGGOvGGPPtDymD2PsMGMskTGWzRhTm6vGGFvNGPuPMfaa\nMfaEMbabMVZb91dQNn36KVCunDJtiEG1rZxawaKccvrfkTtH9F8p0bkxYzTkHdCQScqs1FT5Tu0K\n7doZZsuO/Oz4dwciH0YK6eGew1HPrp6EEZVsJaqhAqAcgO0AVhVUUIUFgNMApgAqQ6/FLgMIAFAH\nQGcADMBhRkOzdaJKFaBfP2V6/34gLk6/dZY3Ki/apPDInSPgPK8fPymu7t3TkJeqIZOUWZs2Aa9e\nKdNffCFdLADwNustph+fLqRNjU0x22e2dAGVAiWqocI5D+ScLwdwrRDHhHLO5wI4BnkDRFOZXznn\nf3LO4znnfwGYAcAJgIsOwiYQfzPOzjbMrsqd3ZRTjh6+fIh/nvyj/0oJIQbDufi2j4MD0KePdPEA\nwNrItYh5HiOkx7cYD0crRwkjKvlKVEPFEBhjFgCGA4gBQF/ddKR1a/keQAq//ipf90CfOrt3FqXp\n9g8hpcvp0/Ld2RVGjADKl5cunhdvXmDOyTlC2sbMBlPbTM3nCKINaqjkYIx9wRh7AeAFAD8AnTnn\nmRKHVWowJp6q/OQJ8Ntv+q2zVuVacK7oLKRVR+ATQko+1V3aZTL5HmNS+vHcj3j6+qmQ/tb7W1Qy\nrSRhRKWD5A0VxtiCnEGueT2yGGO1DBBKKABPAG0B3AKwgzEmYdu89Pn4Y/n6BgorV+q3PsaYqFfl\nVNwppGWk6bdSQohBxMcDu3Yp0z17ipdCMLRHLx9h8dnFQrpGxRoY3awYrOFfChSHTQkXA9hQQJmY\nAl5/Z5xzRW/KHcbYBQDPAfQBsC2/4yZMmICKFSuK8gYNGoRBgwbpK9QSy9ISGDJE2UA5cwaIigK8\n9LhnYGf3zvgl6hcAwJusNzgdf1rtlhAhpORZtQrIUu7zh6+k3SAZ35/8Hq8ylKN6v/f9HqbGphJG\nVHxs2bIFW7ZsEeWlpGi/B5vkDRXOeRIAAyysLq9Oy3IyyAfemhRUcOnSpfDS5ydtKTNmjLgnZcUK\n/a5W28G1A2RMhmwuX2Xu8H+HqaFCSAmXliYekN+woXxaslRuJN7Amsg1ynjsG+Ljhh9LF1Axo+nL\ne1RUFJo0aaLV8ZLf+ikMxpgTY6wxAGcARoyxxjkPC5UyNxhjvVTS1jnH1Ie88VEn5xiHnNddGWPT\nGGNeOedvBWAHgNcADhry+sqCunWBTiobh27ZAjx+rL/6rM2sRatBHomhAbUliaaOyUENqLeyrNu8\nGXimspLWl19Ku6/P10e+RhZXdu8s7LgQRjKjfI4ghVGiGioA5gCIAjALQIWc51EAVJtlNQGo3ovp\nCeAKgH2Q96hsyTlGsV5qOgBvAAcA3M55PQVAK865yhZXRFdUu2jfvtX/VGXVacr/PPkHCS8S9Fsh\n0RmNDZWG1FApyziX98QqWFvLx79J5eidozhw+4CQ7uTWCV09ukoXUClUohoqnPNhnHMjDY9TKmWM\nOOcbVdIhnHOZhmPm5Lz+kHPenXNelXNuyjl35pwP5pzfluIay4KuXQEPD2U6KEjeYNEXPw8/UZqm\nKRNScp06Bfz9tzI9YgRgbi5NLJnZmZh4ZKKQljEZlvgtAa0VqlslqqFCSgeZTN5Vq/DoEbBjh/7q\na+7YHFYmVkKaGiqElFyqvSkymXjZA0NbF7VOtJDkCK8RaGDfQLqASilqqBBJBASIpyqr/vHRNWOZ\nMTq4dhDSR2OOCoNrCSElR1wcsHu3Mt27N+DsnHd5fUp9k4rv/vhOSFuZWGGO75x8jiBFRQ0VIgkr\nK2DYMGX64kXg/Hn91efnrrz9k/g6EVEPo/RXGSFEL4KC5FtwKIwbJ10s80/PV1vczd7CXrqASrFC\nN1QYYyaMsbaMscGMsZGMsb6MMVd9BEdKt9wj9Zcv119ducepHLxNE7oIKUlevwZ++UWZbtQIaNtW\nmlhin8di6fmlQtq1kiu+aiHxQi6lmNYNFcZYa8bYdgDJAI4DWAbgO8hXdP2PMXabMTaZMWaZ33kI\nUfDwALp1U6bDw4EHD/RTl0slF9S1rSukqaFCSMkSFgY8f65Mjxsn3ZTkqRFT8TZLOQPgh04/wMS4\nwGW3SBFp1VBhjO2FfIXWuwA6A7DknFfmnFfnnJtDPiV4LoAOAG4xxjrleTJCVKhOVc7MFO+Eqmvd\na3YXnl98cBFPXz3NpzQpDuLjNeSlaMgkpVruKck2NoC/vzSx/Bn/J3ZEK0f/e9fwRr+6/aQJpozQ\ntkflAABXzvkUzvlpzrlowxTOeUzONOAukDdWaKQi0UrHjvJF4BTWrAHS0/VTV7eayu4bDo7f//td\nPxURnRk7VkPeQQ2ZpFSLiAD+UU6uwYgRgJmZ4ePIys7CuEPigTE0HVn/tGqocM7XcM4ztCwbzTk/\n9hZlm10AACAASURBVG5hkbKCMfGAuMREeRevPrSp0QaW5ZV3JlUXaSKEFF9LliifGxvLt+KQwi9R\nv+DKoytCekjjIWharak0wZQhRZ71wxgrzxirzhirofrQZXCkbBg8GKikshP60qXyrl5dK2dUTrTP\nz+E7h5GZnan7igghOvPvv8DvKp2fAwZIs0ty0uskfHv8WyFtWd4SizouMnwgZVBRZv3UZIydBpAG\nIA5AbM7jbs6/hBSKhQUwcqQy/e+/wOHD+qlLdZxKcnoyzt07p5+KCCE6sXSpOD1xouZy+vbt8W/x\nLE25wVCgTyCqVKgiTTBlTFF6VIIhH4PSA/I9drxyHu/l/EtIoX35pbxLV+HHH/VTT9ea4j046PYP\nIcXX48fApk3KdLt2gJYb7upUZEIk1kYqNyWrZ1cPY5vTWClDKUpDxRPASM75Ic75X5zzq6oPXQdI\nygZHR/EGdBERwF9/6b6eKhWqoElV5V86mqZMSPG1cqV4H7BJkwwfQzbPxpeHvgSH8n70ii4rUM6o\nnOGDKaOK0lCJBmCr60AIyf1HSHUAnS6p3v659uQaTXclpBhKS5OvRKtQqxbQvXve5fVl09VNOHdf\neYv4w3ofooNbh3yOILqm7ToqVooHgKkAfmCM+TDGKqu+lvM6IUXSuLF8urLCli3A/fu6r0d1mjIA\nHLp9SPeVEELeycaNQFKSMj1hgnwTQkNKSU/BlIgpQtq8nDkWd15s2CCI1j0qyQCe5zyOAmgJ4BiA\nJyr5ijKEFNnXXyufZ2YCP/2k+zqaOTaDnbmdkKZxKsXXzz9ryOumIZOUKtnZ4kG0lSsDQ4YYPo7Z\nJ2bjyasnQnp6m+moUZEmtxqaccFFAAC+eo2CkBydOwMNGigXd1qzBpgxQ7zT8ruSMRm61uyKjVc3\nAgAiYiKQlpEGs3ISrCBF8lVDw2cCfVCUfgcPAjdvKtNffAGYmxs2hr8f/42fLiq/Kblbu2NSKwkG\nyRCtF3w7qe1D3wGT0o0x8fTDlBRg3Trd16M6TiUtMw0RMRG6r4QQUiSq49PKlzf8Am/ZPBsj949E\nFs8S8pZ1WQZTY1PDBkIAFG5TwimMMTOVdGvGmIlK2pIxFqT5aEK05+8PVFFZnmDZMvltIF3q4tEF\n5WTKUft7bu7RbQWEkCK5cgX44w9l+uOPxX8PDOGXyF9w/v55Id2rdi/0qNXDsEEQQWGGJi0AoNoB\nfwiAo0raHMBIEPKOTEzk66ooxMUBv/2m2zqsTKzQ3rW9kN53ax+ysrPyOYIQYgiLc41VNfQCb49f\nPsa0Y9OEtEU5C6zouiKfI4i+FaahknvXJdqFiejNqFHie9KLF+t+Wf1etXsJz5+8eoILDy7otgJC\nSKHExABbtyrTfn7yMWuGNOnIJCSnJwvpQJ9AGhclMQNP9iJEOzY2wPDhyvSlS8BJHY+A6lm7pyi9\n5wbd/iFESj/+KJ/xozB1qmHrj4iJQNg15a6ojR0a46uWXxk2CKKGGiqk2Mq9bsLChbo9v6OVo2jn\nUxqnQoh0Hj8G1q9Xpps3B3x8DFd/emY6Rh8YLaQZGNb0WANjmbaTY4m+FLah8hljbBxjbBzkU5sD\nVNKf6T48Upa5ucl3SlU4fBiIitJtHaq3f24m3cSNxBu6rYC8ky1bNORd05BJSrwVK4D0dGV62jT5\nLEBDWfjnQtx+dltIj2wyEi2qtzBcACRPhWmoxAMYAWBCzuMRgMEq6c9yyhCiM9OmidO67lVRbagA\ndPunuNHYUPmHGiqlTWqqeLn82rWBXr3yLq9rNxNvYsGfC4S0g4UDFnRckM8RxJC0bqhwzl04564F\nPfQZLCl7GjcGuqpseBweDty+nXf5wmpg3wCulZT/ben2DyGGt3YtkKwcv4qpUw23XH42z8aoA6Pw\nNku5++FSv6WoZFrJMAGQAtEYFVLsffON8jnnwA8/6O7cjDFRr8r5++fx+OVj3VVACMnXmzfiBd4c\nHeVrpxjKr1G/4sTdE0K6o1tHfNTgI8MFQAqk7aaEWv/UGGNOjLHWRQ+JEDFvb6C1yv+okBDgwQPd\nnb9XHWVDhYNj3619ujs5ISRfmzYBDx8q0xMnylejNYT7qfcx+ehkIW1mbIbV3VeDGXJwDCmQtj0q\nXzDGruesTls394uMsYqMsW6Msc0AogBU1mmUpMxTHauSkSHesOxdtanRBjZmNkJ6141dujs5ISRP\nWVniHlJra2DECMPUzTnHFwe+QOqbVCFvbvu5cLdxN0wARGva7vXTDsBUAJ0A/MMYS2WM3WaMXWOM\n3QeQBGA95INpG3DO9+ojWMbYdMbYGcbYK8bYMy2P6cMYO8wYS2SMZTPGGhVQ/lBOuZ75lSOG1b27\neOGnNWuAZ1r9DyiYscxYtDz20TtHkZKeopuTE0LytGuXeMzZ2LG63YA0P9v+3Yb9t/YL6eaOzfFV\nC1ozpTgqzGDavZzzTgAcAAwB8DOAMACzAbQAUI1zPo1zrs8b/OUAbAewqhDHWAA4DWAKgHzXNmWM\nTQCQVVA5YniMiXtVXr4EVq7U3fn71+0vPM/IzhD9ASOE6B7nwKJFyrSZmXjrDH16+uopvjykrKyc\nrBzW9VwHI5mRYQIghVLowbSc80TO+W7O+XLO+ULO+a+c80jOeXbBR78bznkg53w5gGuFOCaUcz4X\nwDHks+w/Y8wT8mnWw/MrR6QzcCDgqjKvbPly4NUr3Zy7k3snWJZXfpULvx6umxOTd+LkpCHPSkMm\nKXGOHgUuX1amP/0UsLMzTN3jD49H4utEIf2t97doYG/gtfqJ1mjWD4CcXaHDAIzmnD+ROh6imbEx\nMFk57g1JScD/2bvv8CiKN4Dj3zeF3pGqFOlNQFB+KgKCdOkIAoqIYEEQxYINC4JdRBGwK01CEwWU\nKhqaAkJUQHpRQbpCCCWQkPn9sZcryaXc5S53l7yf59mHzNzM7ruk3Hu7szOffeabfeeLyEfnmp3t\n5aV7l3L20lnf7Fx5zd1Vs0m3+fBSmgoIY+Dllx3liAh4/PHsOfa3u79l5taZ9nK90vV4ptkz6fRQ\ngaaJimU8sNYYo9f7g9w990Dp0o7ym29ajzf6gvPtn/jEeBbvWeybHSulXERHw7p1jnL//lC5sv+P\ne+biGR789kF7OUzC+KzLZ+QJz6bHjJRXAr6IgYi8hjVQNy0GqG2M2e2n43cBWgENvek/YsQIihYt\n6lLXt29f+vbt64PoVEr581uPLyaPV/nnH2t9kCFDsr7vdtXaUSCyAOcTzgMwb/s8etftnUEvpZSn\nxoxxfB0WBs8+mz3HfXL5k/wT55jb4NH/PUqTK5tkz8FzsaioKKJSTDMdG5v5BxbEmMCOGxWRkmT8\nOPN+Y0yiU58BwHhjTIl0+qQ8TiXgANDQGLPFqX488DCuA2jDgSRgtTGmVRr7awRs3rx5M40aNcps\nGMoH4uKsT1/JT/1UqAB79/pm7oXec3szd/tcAApEFuDEkycoEFkg6ztWSgGwdq01N1Kyu+6y5lLx\nt2V7l9H+y/b2cpXiVdjy4BYK5ino/4OrVGJiYmjcuDFAY2NMuqu4eXzrR0RaehuYO8aYf40xuzPY\nEjPeU+YO56buNaA+0MBpA3gEGOij4yofKlzY9X72wYPWJHC+cHsdx+2f8wnnWbp3qW92rJQCXK+m\niMBzz/n/mKcunGLQwkEudZ90/kSTlBDhzRiVpSKyT0RGiUi2Dr+3zXrbAKgEhItIA9tW0KnNThHp\n6lQubutTF+tpnlq2PmUAjDHHjTHbnTdb14PGmL+y7+yUJ4YNg2JOS3G8+qo1EVxWdazekXwR+ezl\nr3Z8lfWdKqUA2LABli93lHv3hlq1/H/cR5Y+4nLL5+EmD9PqarcXy1UQ8iZRuRJrDpXbgf22ydR6\ni0h2jEZ6GWvm2xeBQravY4DGTm2qA86DRroAvwKLsK6oRNn6PJDOcXQelSBXpAiMGOEo//mnby4f\nF8pTiA7VHKsgLtq1iPjE+HR6KKUyy/lqCsCoUf4/5tc7vmb6Fscfh+olqvN6ax8vw678ytt5VMYb\nYxpiTfS2G5gMHBaRCbarF35hjBlojAl3s612ahNujJnmVJ5qjAlz0+dl90ex78Mvs+sq3xk+HJzH\nMb/yCiT64Cah8+2fuEtxLN+3PJ3WSqnM2LwZvvvOUe7Rw3W2aX84ce4ED3zr+EwaJmFM7TZVx52F\nmCw9nmwbAPMa1hWWQliTpW0WkTUiUtcH8SmVpmLF4BGnGa/374cvv8z6fjvV6OTyuOKsbbOyvlPl\nlaFD3dR956ZSBb2xY13Lzz/v3+MZY3jwuwc5cf6EvW7kTSO5scKN/j2w8jmvEhURiRSR20VkMfAX\n0A4YhjW9fjVb3VyfRalUGh591HVtEF9cVSmStwgdq3e0lxfsWsC5Sz6aAld55OBBN3Vn3FSqoLZl\nC3zzjaPcuTM09GpCiMyL2hbF/B3z7eV6pevx0i0v+fegyi88nkdFRN4H+mINTJ0OjDTGbHNqck5E\nngAO+yZEpdJWvLh1C+iVV6zynj3w+uvQsWP6/TLyv4J9+QbrL+v5hPNMWLaIdlf2yWK0ylOxsRAT\nk3GdCm4jR7qW/X015XDcYYYudlx5iwiLYFq3aeSNyOvfAyu/8HgeFRFZCXwKzDfGuJ0TVEQigKbG\nmFVZDzE46TwqwePff6FiRTh/3oc7jTwPT5SBvLZp9Hd2gVkLfHgA5bW+XSBKh5CFqhtvhJ9+8t/+\nk0wSHb/syLJ9y+x1L9/yMs+38HN2pDziyTwqHl9RMcbcmok2iUCOTVJUcClZ0pqZdtw4R92YMVm9\nqlKAUTFdWfKPNeglovYSVvx8iiJ5imcpVuWZESNg/PgUdRth/BOBiUd57rHHYJXTu8Frr/n3eBM3\nTnRJUq4rfx1P3/x0Oj1UsPPm1k+XNF4yQDyw1xhzIEtRKeWhp56CDz90rKY8dapVFxnp/T6HFurL\nkigrUUk0CezLM59BjQZl0Ev5UtGikPKCZdFdqetUcPrlF9ckpVMnaNHCf8fbemwrI1c47jPli8jH\ntG7TiAzPwh8CFXDeDKb9Bvja9m/KbRmwV0RWiYh+9FTZplQp1yeA9u7N+my1baq2oUR+xyoNUdui\n0mmtlEop5ViUl9OcFCLr4hPj6Te/HxcvO0YkvNP2HWqXqu2/g6ps4c2ihK2A14HngI22uibAWNt2\nGvgIeBvQj58q2zzxBEyaZA22BOuPYv/+kNfL8XN5wvNwe+3b+TjmYwB+/PNHlu9bzhUFrvBRxCoj\nsfFViDmyP0VdLDFHdDRtsPttY0GWLatpL7fqeApT9gAxR/xzvLd/epttxx3PdXSq0YkHr3swnR4q\nVHgzmHYr8IAx5qcU9U2Bj40xdUWkNfC5Maai70INLjqYNjiNHev6KW7CBHj4Ye/3F/1nNC2n+nR5\nK+WJrX3gGp3HJiRNWQl/Jk9TnwRD6kOZP7Ll0CXzl2T70O2ULlg6W46nPOfXwbRY86SccVN/Bqhi\n+3oPoB87VbZ75BF47z04edIqv/IK3HsvFPRy7bFmFZtRrlA5jpy1PgZeU/oapnSb4ptgVSY9GegA\nlIc2ri3EkD9r2Mvtup3m1eenpdPDe6cunOKOeXfw74V/7XWfdvlUk5QcxJtEZTPwlojcbYw5ASAi\npYA3gV9sbaoDOiuTynaFC8PTT1u3gQCOHbNuB6WcxyGzwsPC6VOvD+PXW4+ebD2+laJ5i1K1RFUf\nRaxUzmIMPPyuoxwWBu+9XoKa5Uqk3cnrYxm6zurqkqQMbzKcbrW6+fxYKnC8GUw7GLgaOCQie0Vk\nL3AIqGx7Dazp9Me6766Ufz30EJQr5yi/8YZj3Io3+l3Tz6XsvMCZUsrVsmWu86T07w81a6bdPis+\n3PQhi3Yvspfrla7HG23e8M/BVMB4syjhTqAO0BWYYNu6AHWNMbttbb4xxuhfcxUQ+fO7rsr633/w\n7rtpt89I43KNqX2F48mBab9PI8kkZSFCpXKmpCTX372ICHjhBf8ca8uxLTy2/DF7OW94Xmb2mEm+\niHz+OaAKGI8SFdsaPyuBqsaYpcaYCbZtmTH6l1sFj8GDoVIlR3ncOGsGW2+ICAMaDLCXD5w+wNq/\n12YxQqVynnnzrFWSkw0cCFWqpN3eW2cvnaX33N7EJ8bb695o/QbXlLnG9wdTAedRomKMSQDq+ykW\npXwmTx548UVHOS4O3nzT+/3dWf9OBLGXp/3un4GBSoWqhAR47jlHOV8+/11NGbp4KLv+3WUvd6rR\nieH/G+6fg6mA82aMygx0fhQVAvr3hxqOBw+YMAEOHfJuX1cVuYrWVVrby3P+mMP5BF8uLqRUaPv0\nU2uixWTDh8NVV/n+OFN+m+LyQaFCkQpM6ToFEUmnlwpl3iQqEcAQEdkkIh+JyDvOm68DVMpbERGu\nM2HGx7teZfGU8+2fuEtxfLPzm3RaK1/4+283dbFuKlVAnT0Lo0c7ysWKWU/f+dr2E9tdVkUOl3Ci\nekZRskBJ3x9MBQ1vEpV6QAwQB9QArnXaGvouNKWyrlcvuO46R3nKFNi2Lc3m6epeuzuF8xS2l/X2\nj/8NG+ambrGbShVQ775rTQWQ7JlnoLiPF1E5n3Ce3nN7u1zJHNtqLE0rNvXtgVTQ8eapn5bpbK0y\n3oNS2ScszHVsSlKS95/0CkQWoFedXvbyiv0rOBx3OIsRKhXaTp50/R278sqszQadlkeWPMIfJxwz\n27ar2o6RTb2cIEmFFG+uqAAgItVEpJ2I5LeV9QahCkotW0KHDo7yd9+5rujqibsb3G3/OskkMWPL\njCxGp1Roe+UVa7B6stGjrSkCfGnm1pl8+uun9nK5QuWY1n0aYeL1W5gKIR5/l0WkpO0R5d3AYiB5\naq3PRGScL4NTyldefx2cU+knn7Rm0PRUs0rNqFyssr085bcpeLpellI5xZ9/wuTJjnLt2jBgQJrN\nvbL9xHbuX3S/vRwmYczsOVOnyM9FvElHxwMJQEXA+bGH2UB7XwSllK/Vrw93Oy6G8MsvMHeu5/sJ\nkzCXQbU7Tu7gp4M/pdNDqZzrhRfg0iVH+dVXrUHsvhJ3MY6ec3pyLuGc45jNX+CWyrf47iAq6HmT\nqLQFnjLGpHzQcw9QyU17pYLCmDGQN6+j/Oyzrn9kM2tgw4Euc6p8EvOJD6JTKrRs3gzTneYfv/FG\n6NrVd/s3xnDvwnvZeXKnva5t1baMaj4qnV4qJ/ImUSmI65WUZCWAi1kLRyn/qVDBWl052b598PHH\nnu+nUrFKtK3a1l6e88ccTsef9kGESoUGY+Dxx13r3njD9fZqVo1fP5552+fZyxWLVuTLHl8SHhbu\nu4OokOBNorIGcLqIjhGRMGAk8KNPolLKT55+2vWxydGj4bQXOcZ9je6zf30h8QJRW6N8EJ1SoeGb\nb1wHpPfoAc2a+W7/a/5aw8gVjid68oTnYV6veVxR4ArfHUSFDG8SlZHA/SKyBMgDvAlsA5oDT/kw\nNqV8rnhx12m+T560nlrwVOeanV0G8+ntH/+YONFNXUc3lSrbXLpkDUZPFhmZteUpUjoSd4Te83pz\n2Vy2101oP4Hrr7zedwdRIcWbeVS2YU30thZYgHUraD5wrTFmn2/DU8r3hg1zXSjtvfdcp/7OjDzh\nebinwT328q9Hf2Xz4c1pd1BeqVjRTV1RN5Uq20yaZN02TTZ8OFSt6pt9J1xO4I55d3D07FF73YAG\nA7i/8f3p9FI5nVcPoRtjYo0xrxhjehtjOhpjRhljjvg6uJRE5FkRWSci50Tkv0z26S4iy0TkpIgk\niUiqRRVFJNr2WvJ2WUQmu9ufCn1588LbbzvKCQnwxBOe72dwo8EuZb2qonK6f/91XZaiZEkY5cOx\nrU9//zRr/l5jL9cvU5/Jt03WdXxyOa8SFREpJiJtReQuEbnbefN1gClEAnOADzzoUxBrXM1IIK0J\nLwzwMVAGKIs1N4xOeZiDdetmTQSXbMECWLnSs31UL1nd5THJmVtncvbSWd8EqFQQSjmma/Roa10f\nX5ixZQbvrHcsF1c0b1Hm955PgcgCvjmAClkeP/EuIp2BL4FCwBlc3/wN4LcFUIwxo20xZHpKIWPM\nDFufSkB6afl5Y8yJrEWoQoUIjB8P117rmPhtxAj49VcI9+ChgsHXDib6z2jAWqhw9rbZDGqki4ur\nnGfnztSTuz3wgG/2venwJgYvdL1COa37NKqW8NE9JRXSvLmiMg74HChkjClmjCnutJXwcXzZ6U4R\nOSEiW0Xk1eSlAVTO1aABDHb627h1K3z2mWf76FmnJ8XzOR4jmrxpss5Uq3KkJ5+Ey47xrbz9tm8m\ndzt69ijdZnXj4mXH7BajbxlNl5pdsr5zlSN4k6hcCUwwxribSyVUfQncBdwCvAr0B6an10HlDGPG\nQGHHgsiMGgWxsZnvny8iHwMbDrSXY47EsOGfDT6MUKnAW7oUvv3WUW7TxnX9LG9dTLxIzzk9+Sfu\nH3td91rddVI35cKbRGUZcJ2vAhCR11IMZE25XRaRGr46njvGmE+NMSuMMX8YY6KwEpXuInK1P4+r\nAq9MGXj+eUf5xAkYO9azfQy5fojLTLUTN+rjsyrnuHjRerInWVgYjBuX9cndjDEMWzzMZQmKeqXr\nMbXbVF1sULkQTy9Ti8gg4AXgC2Ar1ro/dsaYhR7uryRQMoNm+40xiU59BgDjPbnVZBujcgBoaIzZ\nkkHbAsBZoJ0xZkUabRoBm5s3b07RokVdXuvbty99+/bNbGgqwC5ehDp1YP9+qxwZad0Gqlkz8/u4\nbeZtLN6z2OofFsnBEQcpU6iMH6LNXaKiIOWvUtTWKPpeo79f2eWNN6yJEpM9/DBMmJD1/U7+ZTJD\nFw+1l0vkL8Ev9/1CleJV0umlQlFUVBRRUa6TYsbGxrJ69WqAxsaYmPT6e5OoJKXzsjHG+H1+4ywk\nKvux5nvJKFFpCqwGGtjmjXHXphGwefPmzTRq1Cjzwaug9PXX1uyaydq2tS53Z/ZT4+I9i7lt5m32\n8tiWY3mu+XPp9FCZ0aULLEzx0adLVBcW9vXo85Dy0qFDUKsWnLOtCXjFFbB7t+vszt6I/jOaNtPb\nkJhkff4Ml3CW3bWMW6vcmsWIVaiIiYmhcePGkIlExZsJ38LS2fyapIhIBRFpgLX4YbiINLBtBZ3a\n7BSRrk7l4rY+dbGe+qll61PG9noVERklIo1EpJKIdAGmAqvSSlJUztOtm3XfPdny5TB/fub7t6/W\n3uWT4AebPrD/EVYqVD35pCNJAXj99awnKX+e/pNec3u5/H6MaztOkxSVplC7EfgyEAO8iPV4dIxt\na+zUpjrgfC+mC/ArsAjr8ekoW5/kB+suAa2xxt7sAN4C5tr6qVxCBN5/37rtk+zRR13/SKcnTMIY\ner3jMvY/cf+wYOcCH0epVPaJjoZZsxzl66+HgQPTbJ4psfGxdJrZiZPnT9rr7ml4D8P/NzydXiq3\ny3SiIiKLRaSoU/lpESnmVC4pItt9HaAzY8xAY0y4m221U5twY8w0p/LU5Ks9KbaXba8fMsbcYowp\nZYwpYIypaYx5xhijM3flMjVrus5Qe+iQZwNrBzYcSP4Ix1PtE3/RQbUqNCUkWGNRkolYU+eHZeGj\nbWJSIr3n9eaPE3/Y65pc2YQPbvtAZ55V6fLkx64dkNep/CzgPEYkAvBg+KFSwee556BCBUd53Dhr\noqvMKJ6/OHfVv8tejv4zmq3Htvo4QqX8b/Jk2OZ043vQIOuKireMMQxfMpzl+5bb6yoUqcCCPgvI\nF5EvC5Gq3MCTRCVlyqspsMpxChaEd991lJM/WWZ2zLnz7R+A8evH+zA6pfzv6FF44QVHuVgxePXV\nrO1zwoYJfLDJsfJJoTyF+Lbft5QtVDZrO1a5QqiNUVHK77p3h3btHOXvv4d58zLXt0HZBrSs7FhE\n6MutX3Ikzu/rdSrlMyNGwJkzjvLYsVCqlPf7W7RrESOWjbCXwySM2bfPpn6ZVOvDKuWWJ4mKIfWi\nfjpXuMpxRKx5IpwH1o4YAXFxmev/+I2P27++dPmSTgCXBc634ex1RdxUKp9YutR1AG3Dhllbz+e3\no7/R96u+GKe3infbvUvH6h2zEKXKbTy99TNFROaLyHwgH/ChU/lzv0SoVADUqGE9mpnsn39cZ7BN\nT4fqHah9RW17+YNNH3DuUiYfH1IuJk1yU3ebm0qVZefPw0MPOcphYfDxx96v53M47jCdozpzLsHx\nsz/s+mE8/L+H0+mlVGqeJCpTgeNArG2bARx2Kh/HjysnK5Xdnn0WKlVylCdMgI0bM+4XJmE8duNj\n9vKp+FNM+W2K7wNUyodefhkOHHCUhw3zfgDt2Utn6RzVmUNnDtnrOlTrwPj2OmZLec7jmWmVRWem\nzR2WLIGOTlep69eHTZtcbwu5E58YT6V3K3H83HEAqhSvwu5huwkP8/vEzUp5bOtWaNQIEm1zsF15\nJWzfDkWKeL6vhMsJdI7qzLJ9y+x19UrXY9296yiS14sdqhzJrzPTKpWbdOjgutbMli3WI8sZyReR\nj2HXD7OX95/az4JdOgGcCj5JSXD//Y4kBazJD71JUowx3LfoPpckpWyhsnzb91tNUpTXNFFRKgPv\nvus6bfjo0bB3b8b9hlw/xGUCuLd/ehu9gqmCzccfw/r1jnKXLtaSEt54/sfnmfr7VHu5cJ7CLLlz\nCZWKVUqnl1Lp00RFqQyULu16FSU+Hh58MOO5Va4ocAUDGgywl38+9DOr/1qdTg+lsteRI64rIxcs\naF1N8Wai2A83fcgra16xlyPCIviq91c0LNvQB5Gq3EwTFaUy4Z57oKVjehRWroRpmRg6/vhNjxMm\njl+zsWs8mJNfKT8yxkq4Y2MddWPHQsWKnu9rwc4FDF3sOtnh510+p03VNmn0UCrzNFFRKhNE4KOP\nIK/TIhKPPQbHj6ffr1qJavSp18de/n7/92w4tMFPUSqVeVFRsHCho9yokfWkj6d+Pvgzfb7qElwQ\n8gAAIABJREFUQ5JJste9fuvr9G/Q3wdRKqWJilKZVr06vPiio/zff9a8ExndAnr25mddys6Xx1X6\nhg51U/edm0rlkWPHXBcdjIyEL77wfM6UnSd30imqE/GJ8fa6YdcPY2TTkT6KVClNVJTyyBNPWI8o\nJ/vqK5g9O/0+dUvXpUftHvbyot2L+O3ob36KMGc5eNBN3Rk3lcojQ4daiXay555z/bnOjL9O/0Wb\n6W3474JjRz1q9+Dd9u/qasjKpzRRUcoDkZEwZYrrJ8+hQ62F3NLzXLPnXMqvrsniKm9KeWnuXCvB\nTtagATzzjGf7OHb2GG2mt3GZ0O3mijczo/sMnStI+ZwmKkp56NprrU+gyf77L+OngBqVa+Syvsm8\n7fPYcWKHH6NUKrUTJ1xvp0VEWLd88uTJ/D5Ox5+m3Yx27Plvj72ufpn6LOyzkPyR+dPpqZR3NFFR\nygvPPmst2JZswQKYOTP9Ps5XVQxGx6qobDd8uJWsJHv6aSvxzqzzCefpNLMTvx/73V5XrUQ1lt21\njOL5i6fTUynvaaKilBfy5El9C+jhh615KdJyU4WbaFnZ8YzzzK0z+eP4H/4LUikns2a5roxcty6M\nGpX5/pcuX6LnnJ6sO7jOXndl4StZ0X8FZQuV9WGkSrnSREUpLzVoAC+84CifOmVNRZ7eLaDRt4y2\nf20wvBD9QtqNlfKRgwdhyBBHOTzcuuXj/Lh9ei4nXab/1/1Zunepva5k/pIs77+cysUq+zZYpVLQ\nREWpLHj6aWv+iWTffmtNSZ6WZpWa0a5qO3t5/o75bD682Y8RqtwuKQkGDoTTpx11zz+f+ZWRjTEM\n+W4Ic/6YY68rlKcQS+5cQp1SdXwcrVKpaaKiVBZERsLUqa6DEUeMgJ070+4ztpXr7LSjfvTg+nsu\n47wgpL2unptKlaYJE6yZlJM1aWKNscoMYwzDlwznk5hP7HV5w/OysM9Crr8yk5mOUlmkiYpSWVSv\nHrz+uqN84QLceSdcuuS+/XXlr6N7re728tK9S1n791o/Rxma3CYq12iikll//OG6lk+BAjB9upVg\nZ8QYw+PLH2fiLxPtdeESzuzbZ9Py6pbp9FTKtzRRUcoHHnkE2jgtaxIT4zqLbUpjWo5BcEyK9dwP\nz+nKysqnLl2Cu+6CixcddePGQY0aGfc1xvD0908zfv14e12YhDGjxwy61urqh2iVSpsmKkr5QFiY\n9RRQyZKOujfegFWr3LevW7ou/a7pZy+v/ms1y/Yt82+QKlcZNQp+c5oAuWNHeOCBzPV94ccXePOn\nN+1lQZjabarLulVKZRdNVJTykfLl4RPHrXyMgf79raeB3HnplpcIF8csnk+ueJLLSZf9HKXKDZYs\ngbfecpSvuAI++8xaXDMjY1aNSbXK92ddPuOu+nf5OEqlMkcTFaV8qHt3GDzYUT54EO67z/0jy9VK\nVOP+xvfby9uOb+PzXz/PhihVTnb4MNx9t2vdp59C2UxMdfL62tdTPTL/UaePGHjtQB9GqJRnNFFR\nysfGj7dWWk721VcwaZL7ti/d8hJF8haxl0f9OIq4i3F+jlDlVJcvW+NSTp501A0fDl0zMazkzXVv\n8sxK10V/JnaY6JJMKxUImqgo5WOFCkFUlOsjy48/Dps2pW5bumBpnr3Z8azo8XPHeX3t66kbKpUJ\nr7wCP/7oKDdqBG++mXb7ZGNWjeGp759yqXun7TsMbTI0jR5KZZ+QSlRE5FkRWSci50Tkv4x7gIh0\nF5FlInJSRJJExO1i5iJyo4isFJGzIhIrItEiksl5G5Vy1bix9YRFskuXoHdv10m3kj1ywyNUKlrJ\nXn5n/Tv8Hft3NkQZ/P5289+g/zfurVoFox0TH1OokDVlfnqzzxpjGPXDqFS3e95o/QYjbhzhp0iV\n8kxIJSpAJDAH+MCDPgWBNcBIwO3znyJyI7AEWApcZ9smAklZCVblbkOHwu23O8oHDsCgQanHq+SL\nyMfrrR1XUeIT43l2ZSZn5Mrhhg1zU7fYTWUud+IE9OtnzUKb7KOPXG9BpmSMYeSKkakWx3yn7TuM\nbDrST5Eq5bmQSlSMMaONMe8BWz3oM8MYMxZYCaQ15v0d4F1jzFvGmJ3GmD3GmHnGmAQfhK1yKRFr\nEGOVKo66+fNh4sTUbe+oewc3XHWDvfzl1i/56eBP2RClCnWJidCnjzWINtm991qJS1qMMTy69FHe\n/vltl/qJHSbqlRQVdEIqUfEHESkF/A84abutdNR226dpoGNToa9oUZgzJ/V4lfXrXduJCOPajnOp\nG/LdEBKTErMhShXKnnsOfvjBUa5Tx5o2Py1JJokh3w1hwkZHI0H4uNPHOiZFBaVcn6gAyZ93XwQ+\nAtoBMcBKEakasKhUjtG4MbzzjqOckAA9esCRI67tbqpwE/3r97eXtxzbwqSNaTwupBTWE2XOg2UL\nF7au2hUs6L59YlIigxYO4qPNH9nrwiSML7p+wX2N7/NztEp5JyLQAYjIa8BT6TQxQG1jzG4/hZCc\nrH1ojJlm+/oxEbkVuBd4Lr3OI0aMoGjRoi51ffv2pa+7RUpUrvXQQ7BmDcyebZWPHIFevaxPws5X\nW95q8xYLdy0k9mIsAM//+Dy96/amXOFyAYhaBbMdO+Cee1zrpk6FmjXdt49PjKfvV335Zuc39rpw\nCWd69+m6fpLyq6ioKKKiolzqYmNjM90/4IkK8DbwRQZt9vvx+Mmfa3ekqN8BVMyo8/jx42nUqJHP\ng1I5i4g1M+j27bDVNsJq3Tp49FGYPNnRrkyhMrzS6hWGLbEGjMZdiuPx5Y8zs+fMAEStglVcnHVV\n7uxZR93TT1sTDrptfzGObrO78cMBxz2iiLAIZvWcRc86Pf0crcrt3H14j4mJoXHjxpnqH/BbP8aY\nf40xuzPYfHWjPtVTP8aYP4HDQMrPITWAv3x0XKUoWBC++QaKF3fUffCBlcA4e/C6B2lUzpH8Rm2L\ncnmDUblbUhIMGAA7dzrqWreGsWPdtz95/iStprVy+RnKH5GfhX0WapKiQkLAExVPiEgFEWkAVALC\nRaSBbSvo1GaniHR1Khe39amL9dRPLVufMk67fgsYLiI9RaSqiIzBSlxSvIUolTVVqlhzW4Q5/eY9\n9JDr4NrwsHA+uO0Dl9WVh3w3hAsJF7IxUhWsXngBvv7aUa5Y0ZpgMDw8dduDsQdp9kUzNh12zDZY\nNG9RVvRfQYfqHbIhWqWyLqQSFeBlrIGuLwKFbF/HAM7Xj6oDzoNGugC/AouwrqhE2frY1xG1PfL8\nGtZjyr8BLYHWxpgD/joRlXu1bQuvvuooX7pkXbI/eNBR1+TKJi5Tl+/+dzejVznN5pVLuHuUe2JH\nN5W5xIwZ1uyzyfLmtQbUXnFF6ra7Tu6i6edN2XnScemlTMEyrB64mqYV9aFGFTrEuFstTWVIRBoB\nmzdv3qxjVJTHjLHmvpgzx1FXvz6sXWs9uQFw6sIp6kyuw9GzRwHr6YwNgzdwXfnrAhCxCrSff4Zb\nbrES22RRUdbPUUrrD62nc1RnTp53LPpzdbGrWdF/BVVL6MOMKvCcxqg0NsbEpNc21K6oKJUjiMDn\nn8O11zrqtmyx3nQSbSOyiucvzoe3fWh/Pckkce+Ce7l0+RIqd/nrL+jWzTVJeeEF90nK1zu+puXU\nli5JSr3S9Vh771pNUlRI0kRFqQApWBAWLYLy5R11ixfDY485yl1rdeWOunfYy1uPb+W1Na9lY5Qq\n0OLioHNnOH7cUde7N7z4Yuq27294n55zehKfGG+vu6nCTay6ZxXlC5dP3UGpEBAMjycrlWtdeSV8\n+y00awbnzll1778P+fI5Pi3fd+X7LNuzktOXrE/IY1aPpVpSZ2oX01uOOV1CAowY4XikHeC66+CL\nL1wHZCeZJEauGMm4n11nN+5Ruwczus8gf2T+bIpYKd/TMSpe0jEqypcWLYKuXVMvWGhXLwpud1q8\n5VRlmD8DEvUNKMcyQGwluFDSXlX8iktELd9FqbKOZcguJl7kxegXWbF/hUv3R/73COPajiM8zM3j\nQEoFmCdjVPSKilJBoHNna5r9EU7rweXNC5MmWeNYjOnDU5vns/LIPOvF4n/CoJsDEqvKJt+/Amud\nVtHOc4ZTPVvQfuFvGXZ97IbHGNduXIbtlAoFmqgoFSQeeQT27HHMVHvxorWA4erVUL++MLv2h1zz\nwTqOnHUsEvRu+3dpVrFZgCJW/jL781K8ubaCvRwRmcSEacf4X3PH1E57/9vLiGUjOBznWDY5T3ge\nXm75Mg83eThb41XKnzRRUSpIiFir3h47Zs2NARAbC+3bW9PtX311SaZ1n0ab6W3sfV5d8ypbHtxC\nmUJl0thraIuKgpTLZkVtjcrRa9PMmwdvPe9aN3VKGP3uqG4vf7v7WwYtHMTZS4459IvnK87Cvgu5\nuaJeaVM5iz71o1QQCQ+HL7+EVq0cdUeOWJPEHTsGrau05rEbHI8FHT93nH7z+3E56XIAovW/FOuY\nWXXb3FTmEIsXQ79+rmOV3nrLqgMwxvDWurfoEtXFJUmpWbImGwZv0CRF5UiaqCgVZPLmtaZIdx6j\nvXcvdOgAp07Bq7e+Sv0y9e2v/XDgB16Kfin7A1U+9eOP0LOn9aRPskcftW7/gTVoduCCgYz8fiTG\nadmytlXbsn7weqqXrI5SOZEmKkoFoSJFYMkSqO703vPrr9CuHcSfy8vcXnMpnKew/bWxa8ayZM+S\nAESqfOHnn60B1fGO6U8YMADGjbNuCR6OO0yraa2Y+vtUl37Dmwznu37fUSxfsWyOWKnso4mKUkGq\ndGlYvtx1QrhffrHGrJSNrMFnXVzXzOwzpw93P3w3PXr0YOjQoURHR6PTDwS/X3+1rpYlz6MD0KsX\nfPqpNVfKmr/W0Pjjxvx08Cf76xFhEXzU6SPe6/AeEWE61FDlbJqoKBXEKleGlSuhjNNY2fXrrTe2\n9hV7MbzJcHv9mcQzTE+Yztfffs3kyZNp2bIlHTt25MyZM9kfuMqUzZvh1lutQdPJOnWyFh8MDze8\nu/5dWk5taV/vCaBE/hIsv2u5y6KVSuVkmqgoFeRq1YIffoBSpRx1P/0Et90Gz9/wJkXjnBYLLwd0\nA8QqLl26lDvuuEOvrAShDRusJOXUKUfdrbfC3LmQwDn6ze/HiGUjuGwcA6Ublm3IL/f9QsurWwYg\nYqUCQxMVpUJAnTpWsnLFFY66NWugxc3xxH5iIM6pcV2guaO4dOlSVq9enV2hqkxYuxbatDEuV1LK\nl9/D44+v5q+zu7jhsxuYtW2WS5+7G9zNunvXUaV4lWyOVqnA0kRFqRBRrx58/z2UKOGo2769KJxZ\nDbOKQ6JT45ZYCYvNnDlzsitMn6pQwU1dETeVISQ6Gtq3N8TFiVPt9xw+3JCOT7eg7oS6bDu+zf5K\nZFgkkztOZkrXKRSILJDt8SoVaJqoKBVCGjRIfRsIGsA/62FhCdfG3YFK1pdHjx4lFE2a5KbuNjeV\nIWLBAujQwXDunHOSsgTydIJu56EHXA533OopX7g8q+5ZxZDrhyAiqXeoVC6giYpSIaZBA+u2z1VX\nOdfWgC2/wprijqoIoC9QGsqWLZu9QapUPvkEevSA+HjnhGMhlO0C91+Ehq7tGxVrRMz9MdxY4cZs\njVOpYKOJilIhqGZNa5xD+fIXnGorwg+7YJvTQJZ8wF3QvEvzlLtQ2cQYePlluP9+SEpyfmU2XN8d\nBieC07eMJGAlNNndJMcujaCUJzRRUSpEVaoEv/ySj0KFDjgqTSn4eh/sL+2oKwIv7nnR5RFXlT0S\nE2HoUHjxxRQvFB4Ld/WB25JcV1yLBaYAa+D40ePZFqdSwUwTFaVCWPnywtatJSle3DH4kstFYPYe\nOOJIVnb9u4tbp93K8XP65pddTp+25kT54APX+qpdn4CHnodqKTrsAD4E/raKertOKYsmKkqFuMqV\ni/DPP3Vp2dIpCblYBL78ncIJjkdZt5/Yzq3TbuXEuRMBiDJ32bMHbrgBli1z1EnBf7nujTvYd+04\nyO/UOAFYDMwGnO7k9e7dO3uCVSrIaaKiVA6QP7/w/felGTnSqfJsWeIm/kje85XtVduOb6P19NZ6\nZcWPVq6E//0Pdu1y1OWrv4hiz9Rj04UUj4n/A3wEbHStbt++Pc2b67gipUATFaVyjLAweOMN61ZD\neLitMrYiFz/+kbC4ivZ2W45t4ebPb+av038FJtAcyhhrEcF27Zxmmy10lAIDehPfowunEh1jhCLC\nIqh2sBp8Bpx03U/79u2ZPXu2Po6slI0mKkrlMA8+aC1mWLKkreJ0ZZI+i4ZYx0Rpe/7bQ9PPm/LH\n8T8CEmNmDR3qpu47N5UBduoUdO8OTzwBly8DGGj0KRGP1Ob81XNd2ta+ojY/D/qZ3Z/sJvqHaB56\n6CF69OjBQw89RHR0NIsXL6ZIkSIBOQ+lgpEuu6lUDtSqlbXgXffu1uq8nL4avlgN/dtAyb0A/BP3\nD82+aMbCvgu5ueLNgQ04DQcPuqk746YygDZtgt694UDyw1cld0HnB6DyKpfJgsMkjBE3jGBMyzHk\nj7QGqbRo0YIWLVpke8xKhRK9oqJUDlWpkjXXyp132ipOV4bP18IRx8xip+JP0WpqK6b8NiUQIYa0\npCQYPx6aNrUlKXnPQJuR8NA1UHmVS9try17LxsEbebvt2/YkRSmVOZqoKJWDFSgA06fD++9DnjzA\nuTIwJRr+dAzUTEhKYOCCgTyx/AkuJ11Oc1/K4a+/rJWOH3sMLiUkQYNpMKwmNH0LwhPs7fJH5OfN\n1m+y8b6NNC7fOIARKxW6QipREZFnRWSdiJwTkf8y2ae7iCwTkZMikiQi9VO8XslWf9n2r/PW0z9n\nolT2EYFhw2DjRmtGWy4WhRnLYEs/l3bjfh5Hhy876BNB6TAGpk6F+vWtxQUp/wvcezN0HwCFXSfU\na1OlDVuHbOXJpk8SEaZ32ZXyVkglKkAkMAf4IKOGTgoCa4CRgHHz+t9AWaCc7d+ywItAHLAkK8Eq\nFUwaNLDGrQwaBCTmg/kzYOUrLm1W7F9Bww8bEv1ndEBiDGb79llP9NxzD5yJ3A29esP9TaDCzy7t\nKherzPze81l21zKqlqgamGCVykFCKs03xowGEJEBHvSZYetTCUj1vJ8xxgAuHyFFpDsw2xhzPksB\nKxVkChaETz+13nAfekg4ueZZOFEbevSHPOcAOHL2CLdOu5VRzUYxqvkoIsMjAxx19jHGsGrVKubO\nncuRI0coV64c3br1ZuPG5owdK8RHHoZOo6HRZxDmepssf0R+nrn5GZ646Qkdh6KUD4VUopIdRKQx\n1jqmQwIdi1L+0qsX3HILDB8Os2Z1h49/sa4QlLGm4k8ySby8+mUW7FrAlG5TaFi2Yfo7zAHOnDnD\nHXfcwdKlS51qOzJ5cmkodBSavw3XfwCRF1L17V23N2+1eYuKRSumek0plTWhdusnOwwCthtjNgQ6\nEKX8qVQpiIqCr7+GshG14dMNsHmwS5vfj/3O9Z9cz/M/PM+FhNRv0P7Wt6+bunpuKrPIGJMiSbkG\nWA5FJ0PHifDo1XDTO6mSlJaVW7Jh8AZm3z5bkxSl/CTgiYqIvOZmEKvzdllEamRTLPmAvsCn2XE8\npYJBt26wfTsMe6AAYd99AnNnwfkS9tcTkxIZu2YstSbWZt72eRhjMMYQHR3N0KFD6dGjB0OHDiU6\nOhrrTqrvuE1UrvF9orJq1SpbklIdmAplpkLXKBheDZpMhoiLLu0blm3I0juXsvLulTS5sonP41FK\nOYiv/7B4HIBISaBkBs32G2PscyfZxqiMN8aUSKdPyuNUAg4ADY0xW9Jo0x/4BLjSGPNvBvtrBGxu\n3rw5RYsWdXmtb9++9HX3F1apILdli3U7aNWmY3DbQ1Bnfqo2N5Zthvk+gfXz1qd6LXn691CbWbVf\nv9FEza4ENQrDDe+nmgclWYnEEky6YxK96/YmTAL+OU+pkBAVFUVUVJRLXWxsLKtXrwZobIyJSa9/\nwBMVb2QhUdkPXJtOovIjcMIYk+GypcmJyubNm2nUqFFmw1Aq6BkD8+bBM88a9uWdC+0fhcJHUjfc\nnReiL8Jh1+r27duzePHioF+rxhjrEePXJh1ixfGp0PgTKJbG+keHgdXQvU535n+VOnlTSnkmJiaG\nxo0bQyYSlZD6SCAiFUSkAVAJCBeRBratoFObnSLS1alc3NanLtZTP7Vsfcqk2Hc1oDnWFRWlci0R\na7Dtju3CJ4/2pvz8XbD2KUjM49qwxkW4H7i7JNQQ+zN1S5cuTf6kFJRiY2HiRxeo1CmKVp+3Y0W9\ninDrKPdJygFgBvAxsBPKlS2XzdEqpUIqUQFeBmKw5jkpZPs6BnCe8rE64HwvpgvwK7AIax6VKFuf\nB1LseyDwtzFmhV8iVyrEREbC4MGwb3th3uv8OuW+/gO2u5kDscq/0M/Aw8WhaVkoAnPmzMn+gNOR\nmAjzv4uj6f2zKfHAHTz8dykONukH1ZaDpLiqnBAGm4HJwFRgr+Ol3r0zvNiqlPKxkLz1Ewz01o/K\nbRIS4Kab3mXTwfJwyyyo/bX7hkYocKwuw1oOZ0jrjlQucWWWjutubpNevXrRokWLdG8vnTlj+HTB\ndqI2fM/vZ1eQUOH7VINineU9V5SLP8daSYqbB5xC5ZaWUqHAk1s/Oo+KUipTIiOhSZM9bNo0Ama3\nhDKvwQ3b4ZrZEHHJ0VAM58tu480d9/PmDihyviHXF29P+zo30+fmG7iqREZj5x3cz20CkydPTjVw\n98iJeGZF/8Z3v23k95MbOFnkB2v+k1JYm7tzMoXoXqM3w5reQ/1i9enTpw9LLyxN1S75WJqkKJX9\n9IqKl/SKisqNoqOjadmypVNNSSh4OzSOgGsXQ/EDGe4jz9nqlOEari5Um3plalO/UkVqXVWOWleV\no3SxgiTnAsYYWrYcwKpV0x2d8xaEsiVASkORghSrWpJ8FcL5l10kFN0B4YnuD+ok/HIBmpTowJAW\nPeletzOF8hSyv2aMYfXq1cyZM4ejR49StmxZevfuTfPmzTVJUcqHPLmioomKlzRRUbmRMYaOHTum\nusJhqQYVbiSy8WUSaiyFAplaN9RVQn7kcn7kcj5IzEPSNx9Dj7sgMh4i4l2v3Hgg76VyNC7Whvtb\ndKPXte0oEFnAq/0opXxDb/0opfxCRJg9e7bb2zGwl/Z1qzF76mxOxxbkkyUb+OaPxey8vJjEK35L\nPWjVncgLmMgLjtVDIy9AYc9Xc464WIqKYTfQqnJrHmzbmkYVausVEaVClCYqSimPFClShMWLF6d7\ni6RIERhz/02M4SaMGcv2fWeYs24jP+z5iV1nfuVU+A4Si+xNtbCfx5LCCTtdiXJ56lK1aG2aXn0d\nfZs3od5VFTUxUSqH0ERFKeUxEaFFixa0aNEiE22hbrUijK7WmtG0ttfHnb/Emm372XLgCH/9d4TD\nZ45w8ty/nL90kUQTz6Gjf3E6tgRsaAGJBhKS4EIYNNgKS8vD2ZMQe5wHH2jPpPcm+fN0lVIBpImK\nUiogChfIQ8cmtejYpJbb16Ojo2k57T9YkmI6+yrA36fsRZ3bRKmcLdQmfFNK5RItWrSgVKnS6bZp\n3749zZs3z6aIlFKBoImKUiooiQjXXdeY9u3bu31d5zZRKnfQWz9KqaAVERHJggWuA3djKscwJXqK\nzm2iVC6hiYpSKmhNnJh64O7fsX9TsWjFAEemlMoueutHKRW0KrrJRzRJUSp30URFKaWUUkFLExWl\nlFJKBS1NVJRSSikVtDRRUUoppVTQ0kRFKaWUUkFLExWllFJKBS1NVJRSQSsqyk3dVjeVSqkcSxMV\npVTQcpuobNNERancRBMVpZRSSgUtTVSUUkopFbQ0UVFKKaVU0NJERSmllFJBSxMVpZRSSgUtTVSU\nUkopFbQ0UVFKBa0KFdzUFXFTqZTKsTRRUUoFrUmT3NTd5qZSKZVjhVSiIiLPisg6ETknIv9lsk93\nEVkmIidFJElE6rtpU0ZEpovIERE5KyKbRaSH789AKaWUUp4IqUQFiATmAB940KcgsAYYCZg02kwH\nqgOdgHrAfGCOiDTwPlSllFJKZVVEoAPwhDFmNICIDPCgzwxbn0qApNHsRuBBY8xmW/kVERkBNAZ+\n9z5ipZRSSmVFqF1R8Zd1wB0iUlwsfYC8QHRgw1JKKaVyt5C6ouJHdwCzgX+BROAc0N0Ysz+gUSml\nlFK5XMATFRF5DXgqnSYGqG2M2e3HMMYCRYFWWMlKN2CuiNxsjPkjvY4jRoygaNGiLnV9+/alb9++\n/opVKaWUChlRUVFEpVgKPTY2NtP9xZi0xpdmDxEpCZTMoNl+Y0yiU58BwHhjTAkPjlMJOAA0NMZs\ncaqvAuwF6hpjdjjVrwD2GGMeSmN/jYDNmzdvplGjRpkNQymllMr1YmJiaNy4MUBjY0xMem0DPkbF\nGPOvMWZ3BltixnvK3OHc1BWw1V9OUX+ZIPj/USo3GzrUTd13biqVUjlWSL0Ri0gF2yPDlYBwEWlg\n2wo6tdkpIl2dysVtfepiPfVTy9anjK3JTmAf8LGIXC8iVUTkcaA18HV2nZtSKrWDB93UnXFTqZTK\nsUIqUQFeBmKAF4FCtq9jsB4jTlYda7xJsi7Ar8AirCsnUbY+DwDYrtZ0AE4AC7EeR74LuNsYs8yP\n5+K1lPf6coqcel6g5xaKcup5Qc49t5x6XpCzzy0jIZWoGGMGGmPC3WyrndqEG2OmOZWnGmPC3PR5\n2anNPmNML2NMOWNMYWPMtcaYmdl9fpmVU39gc+p5gZ5bKMqp5wU599xy6nlBzj63jIRUoqKUUkqp\n3EUTFaWUUkoFLU1UlFJKKRW0Aj7hWwjLB7Bjx46M2vlcbGwsMTHpPnYeknLqeYGem7dOn4aUuz59\n4HS2/F/q9yz05NTzgpx3bk7vnfkyahvwCd9ClYj0A74MdBxKKaVUCLszo4dXNFHxkm0hZrgzAAAO\nqElEQVRG3XbAn0B8YKNRSimlQko+oDKwzBjzb3oNNVFRSimlVNDSwbRKKaWUClqaqCillFIqaGmi\nopRSSqmgpYmKUkoppYKWJipBSESGisgBEbkgIutF5PpM9msqIgkiEpQP23tyXiLSQkSSUmyXRaR0\ndsacWZ5+z0Qkj4i8IiJ/iki8iOwXkXuyKdxM8/B79oXT98n5+7Y1O2POLC++Z3eKyG8ick5EDovI\nZyJSIrvi9YQX5zZURLaLyHkR2SEi/bMr1swSkWYislBE/rH9XHXJRJ9bRGSz7Xdst4gMyI5YPeHp\neYlIWRH5UkR22X7X3smuWANFE5UgIyJ3AOOwVoi+Fms152UickUG/YoCU4Hv/R6kF7w8L4O1GnZZ\n21bOGHPc37F6ystzmwu0BAYCNYC+wC4/h+oRL85rOLbvk+3fq4D/gDn+j9Yznp6biDTF+v36BKgD\n3A40AT7OloA94MW5DQFeAV7AOreXgEkiclu2BJx5BYHfgIew/jakS0QqA98CK4EGwHvApyLSxn8h\nesWj8wLyAseBMbZ+OZ8xRrcg2oD1wHtOZQEOASMz6BcFjMb64xQT6PPI6nkBLYDLQJFAx+6Hc2uP\n9QZeLNCx+/K83PTvBiQCFQJ9Lj74nj0O7ElRNwz4O9Dn4oNzWwe8kaLubWB1oM8lnXNMArpk0OYN\nYEuKuihgcaDjz8p5pWj/I/BOoOP296ZXVIKIiEQCjbE+AQBgrJ/G74Eb0+k3ELgaK1EJOt6eF9Yf\n2N9sl9mXi8hN/o3Uc16eW2dgE/CUiByyXcJ9S0QynEo6u2The+bsXuB7Y8xB30foPS/P7Weggoh0\nsO2jDNAL+M6/0XrGy3PLS+pJK+OBJiIS7o84s8kNpL7CvIzM//yqIKGJSnC5AggHjqWoP4Z1KT0V\nEakOvIo1DXGSf8PzmsfnBRwBHgB6Aj2Ag0C0iDT0V5Be8ubcqgDNgLpYVx0ewbqVMMlPMXrDm/Oy\nE5FyQAesWyXBxuNzM8b8BNwFzBaRS1g/n6ewrqoEE2++b8uAwSLSCEBErgMGAZG2/YWqsrj/fygi\nInkDEI/ykiYqIUxEwrDWG3rRGLMvuTqAIfmMMWa3MeYTY8yvxpj1xphBwE/AiEDH5gNhWJd4+xlj\nNhljlgKPAQNy0B/Qe7DeyBcEOA6fEJE6WGMcXgIaYS2fcTXwUQDD8pUxwBLgZxFJAL4GptheC9YP\nPyoX0UQluJzEGpdRJkV9GeCom/aFgeuAibanfRKA54GGInJJRG7xZ7Ae8PS80rIRqOaroHzEm3M7\nAvxjjDnrVLcDK8m8yucReier37OBwDRjTKKvA/MBb87taWCdMeYdY8w2Y8wKrMGP99puAwULj8/N\nGBNvjBkMFAAqARWBv4A4Y8wJP8bqb0dx//9wxhhzMQDxKC9pohJEjDEJwGbg1uQ6ERFb+Sc3Xc4A\n9YCGWKPaGwAfAjttX2/wc8iZ4sV5paUh1pt80PDy3NYB5UWkgFNdTaxPr4f8FKpHsvI9syXIVYHP\n/Bii17w8twJYA4OdJWE9pRE0VzGz8n0zxlw2xhy2jWnpAyzyZ6zZ4Gec/h9s2trqVSgJ9Ghe3Vw3\noDdwHrgbqIV1aflfoJTt9deAqen0D9anfjw6L6xxG12w3vDqAu8CCcAtgT4XH5xbQaxPrLOB2kBz\nrEeTPwz0ufjiZxGYDvwU6Ph9/D0bAFwEHsS65dMU6wpf0J2nF+dWHbgT62plE2AWcAKoGOhzSXFe\nBbE+gDXEShIftZUrpHFelYE4rKd/amJdAbsEtA70uWTlvGx1ye1/sf2+NQBqB/pc/PZ/FOgAdHPz\nTbF+of4ELmBl/9c5vfYF8EM6fYMyUfH0vIAngT3AOdsfzZVA80Cfg6++Z1hzpywDzmIlLW8CeQN9\nHj44ryK2c7o30LH74dyGAltt53cIa16VcoE+j6yeG1YyE2M7r1PAfKB6oM/BzTm1sL2RX06xfZ7O\n96w51hWmC7a/J/0DfR4+Oi937fcH+lz8tYntpJVSSimlgo6OUVFKKaVU0NJERSmllFJBSxMVpZRS\nSgUtTVSUUkopFbQ0UVFKKaVU0NJERSmllFJBSxMVpZRSSgUtTVSUUkopFbQ0UVEqlxORAyIyPJ3X\nvxCR+dkZU1pEpJKIJIlIfS/63ioi221r3ygPiMjPItI90HGo3EkTFaUCzJYIJDltJ0VkiYhcE+jY\nAimdBMnb6bTfAF42KabjFpF8IvKfiBwXkUgv9+01EWlh+74Xye5je2As1v+fUtlOExWlgsMSrCXo\nywKtsFbqDfXVa/3F4ysiInIzUAVrHZuUemKt4bMT6Ja10LwiZGIV5kAkUU6WAIVFpEMAY1C5lCYq\nSgWHi8aYE8aY48aYLcDrQAURKZncQETqichKETlvu+rykYgUdHr9CxH5WkQeF5HDtjYTRSTcqU0p\nEVlk28c+EennaaBieUZE9tv286uI9HR6PfkKQSsR+UVEzonIOhGpnmI/o0TkmIicFpEPReRVEfnV\n9tqLWCsWd7Xt67KINHfqXlVEfrDt+zcRuSGDsO8AVhhjLrl5bRAww7YNdnO+SSIySETm2463W0Q6\np2jTxVZ/XkSWi0h/56skIlJRRBbartycFZGtItJeRCoBP9h2c8p2np/b+vwoIu+LyHgROQEstdVX\nEJEFIhInIrEiMltESjvF8qLtezJQRP6ytZsoImEiMlJEjtj+359NcQ4v2drHi8ghEXk3+TVjTBKw\nGOiTwf+zUj6niYpSQUZECgH9gT3GmH9tdQWwVlv+F2gM3A60Bt5P0b0l1pWDW4C7gXtsW7KpwJVY\nK7bejrXKbikPQ3wWuAu4H6gDjAemi0izFO3GAiNs8SYCnzud4522/TwJXAf8Y4sl+bbM28AcrDfn\nMkA54KcU+34Ta3n73cBMEUnv71kzYFPKShGpCtwAzAbmAs1EpIKb/i8As4BrsN6wvxSRYrZ9XG3r\nO98Wz6fAq7jeopoM5AFuBuoBT2GtVvw31hUdgOq283zEqd/dwEXgJuBBERFgIVDMdk6tsb7fs1LE\nWxVoD7TDSi4GA98B5bFWFH4KGCsi19vO4XbgUeA+oBrWlaWtKfa50XZMpbJXoJdv1k233L5hLeOe\nAMTZtiTgENDQqc19wEkgn1NdB6wEoJTTfvaDtSq6rW42MNP2dQ3bvhs5vV7TVjc8g/jm277Og/UG\n+78UbT4BZti+boG17PwtKWK9DOSxlX8G3kuxjzVAjLvjOtVVssV7j1Ndbdu+a6RzDqeAO93UjwW+\ncirPB15I0SYJeMmpXMBW19ZWfh34PUWfMbaYitjKvwPPpxFbC+e2TvU/AptS1LUBLgHlU5x/EtDY\nVn7R9nNUwKnNEmBfin3tAEbavh5hK4en83/YGUgI9O+Lbrlv0ysqSgWHH4D6WJ/Ir8e6erLU6dN9\nLaw3w3inPuuwrorWdKr7wxjj/En+CJB8W6AW1htNTPKLxphdwGkP4qyG9Ua9wnZLIU5E4rCuAFVJ\n0db5E/kR27/JsdQEfknRfqMHcaTctzjt2538gPP/HbYrMAOwbvkkmwkMTO94xpjzwBmn49Ug43OZ\nADwvImttt1gyO1B6c4pyLeCgMeawUzw7sL6HtZ3a/WmLM9kxYHuKfR1zOoe5WN/XAyLysYh0c75l\naHMBCBORvJmMXSmf0ERFqeBwzhhzwBiz3xizGesKSkHbv55ISFE2+Pb3vJDt345YSVXyVgfolU4s\nycmTr2LxdN8ngeIp6tph3QabLSIJIpIARAEVReTWdI6XfMxMn4sx5jPgamAa1q2fTSIyNBNdz2X2\nGCm4izfNczDGHMJKuIYA54FJwKoUyUoJrJ/Ti17GpJRXNFFRKngZrCsBYF2WbyAi+Z1evxnrlsGu\nTO5vJxAhIo2TK0SkJtZ4h8zajjVmopItqXLe/vFgP7uwrhw5S1m+BKT8VA/ePZ78K1Yy5WwQVmLS\nENeka5bttczahTXOxlmTlI2MMf8YYz42xtwOjMORhCYP8HV3rintwBpkfWVyhYjUwfoe/uFBzKkY\nYy4aY74zxjyKNdbpJqwxOcnqYf0/KpWtIgIdgFIKgLwiUsb2dXHgYaxL8cmPKH8JvARMFZHRWJfs\nJwDTjDEnMnMAY8xuEVkGfCwiQ7CSnPFYn6AzxRhzVkTeBsbbPm2vBYoCTYFYY8x0W1N3j9o6170P\nfCIim7EGyfbBuvW1z6nNn0BbEamBNYg4Np19Z2QZ1sBUawcipbDGXHQyxrjcEhGR6cDXIlLMGJOZ\n22IfASNE5HXgM+BarFtKYEuqRGQ81jiR3VhXJlriuBXzl61dZxFZDFwwxri9kmKM+V5EtmEN5h0B\nRGJd/fjRGON1EiEiA7ASpQ1YPw/9bf/+5dSsGbDc22Mo5a3/t3P/PDZFURjGnx3NlPgC8wE0WqKi\nUUmUEkLFiFYhkamEjEhUOo1IKESjISEEUUj4CBKq8aeh01mK90jujBEzl3vuDs8vudU9Z99zbnNW\n9nrXcUdF6sNBYHX4vGSY7Kmq5wBV9ZW0KnaS/MMd4BEpaLbiBJmweQrcJQ/ZT1tZoKqWSVj0HHnY\nPiCtoLeTh2106sQat8lkzBWSw1gEbrA2R3Kd7Fa8Hq5x72bW/oVbwK6JEeljJHD6ZINjH5OH9NFN\n3ss7MkF1mIRmTwEXh69/tEm2AdfI/3Wf7G6dGc5fJQHYFeADP09yrXeIhIOfkcLhDdONDU/e1xey\nw/NiuIf9pIj7DDDs4OwhAWdpVG1t7k6S5qO19hB4X1XHf3vwdOtfJpM1p2ex/rrfOg+crKrFWf/W\nGIbdou1VtTTva9H/x9aPpNENWZsl0pL5BhwBDpD3gszKJRIW/euGVtor0qLaB5wlrbl/xUeSq5FG\n546KpNG11hZI/mY3sEBaPBeq6t5cL2xKrbWr5O23O8hL3G4CK5U3ukr6AxYqkiSpW4ZpJUlStyxU\nJElStyxUJElStyxUJElStyxUJElStyxUJElStyxUJElStyxUJElStyxUJElSt74DE2PfSz4klJ0A\nAAAASUVORK5CYII=\n", | |
"text/plain": "<matplotlib.figure.Figure at 0x106034470>" | |
}, | |
"metadata": {} | |
}, | |
{ | |
"output_type": "execute_result", | |
"data": { | |
"text/plain": "{'B0': 59.029031623822917,\n 'Be': 60.669383302996749,\n 'De': 0.045872631987149484,\n 'ZPVE(anharmonic)': 2170.2152781463119,\n 'ZPVE(harmonic)': 2170.2152781463119,\n 'ae': 3.2807033583476657,\n 'nu': 4123.5266903628035,\n 'r0': 0.7528142730476254,\n 're': 0.74256740791500009,\n 'we': 4412.731844935897,\n 'wexe': 144.60257728654668}" | |
}, | |
"execution_count": 5, | |
"metadata": {} | |
} | |
] | |
}, | |
{ | |
"metadata": { | |
"collapsed": true, | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "", | |
"execution_count": null, | |
"outputs": [] | |
} | |
], | |
"metadata": { | |
"_draft": { | |
"nbviewer_url": "https://gist.github.com/7eb983115833b1f79aeb7c2038bc53f4" | |
}, | |
"anaconda-cloud": {}, | |
"gist": { | |
"id": "7eb983115833b1f79aeb7c2038bc53f4", | |
"data": { | |
"description": "Psi4SureIsNeat.ipynb", | |
"public": true | |
} | |
}, | |
"language_info": { | |
"pygments_lexer": "ipython3", | |
"mimetype": "text/x-python", | |
"version": "3.5.2", | |
"nbconvert_exporter": "python", | |
"codemirror_mode": { | |
"name": "ipython", | |
"version": 3 | |
}, | |
"name": "python", | |
"file_extension": ".py" | |
}, | |
"kernelspec": { | |
"name": "python3", | |
"display_name": "Python [default]", | |
"language": "python" | |
} | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 1 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment