Skip to content

Instantly share code, notes, and snippets.

@oroszl
Created July 11, 2017 18:41
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save oroszl/478a2813fa7cb13a93879b350868286c to your computer and use it in GitHub Desktop.
Save oroszl/478a2813fa7cb13a93879b350868286c to your computer and use it in GitHub Desktop.
pybiding example
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Populating the interactive namespace from numpy and matplotlib\n"
]
}
],
"source": [
"%pylab inline\n",
"import pybinding as pb\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"def mymodel(dim=22):\n",
" s1=matrix([[0,1],[1,0]])\n",
" s3=matrix([[1,0],[0,-1]])\n",
" lat=pb.Lattice(a1=[1,0,0],a2=[0,1,0],a3=[0,0,1])\n",
" lat.add_sublattices(('A', [0,0,0],-5*s1 ))\n",
" lat.add_hoppings(\n",
" ([1,0,0], 'A', 'A', -s1),\n",
" ([0,1,0], 'A', 'A', -s1),\n",
" ([0,0,1], 'A', 'A', -s1+1j*s3)\n",
" )\n",
" model = pb.Model(lat,\n",
" pb.primitive(a1=dim, a2=dim,a3=dim), \n",
" pb.force_double_precision())\n",
" posmin=model.system.positions.z.min()\n",
" posmax=model.system.positions.z.max()\n",
" return model,posmin,posmax"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEKCAYAAADjDHn2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXl8VOW5+L/PmT0QEpYgEkRQcUMsWFxara3dsHqrqbdW\nXHq1m7deva1L6cXWulBbsf663ksXba2916qotSlWe9W63VaLggZFVCpQRIJiIAkkZDKZmfP8/jjL\nnJmETIBMkkne7+eTDzPnvGfmJTnnfd5nF1XFYDAYDIbesAZ7AgaDwWAY+hhhYTAYDIaiGGFhMBgM\nhqIYYWEwGAyGohhhYTAYDIaiGGFhMBgMhqKUVFiIyKkislZE1onIwl7G/bOIqIjMDRy72r1urYjM\nK+U8DQaDwdA74VJ9sIiEgCXAx4DNwAoRWaaqrxaMqwS+CjwXOHYkMB+YCUwG/iwih6pqtlTzNRgM\nBsPuKaVmcRywTlU3qGoXcA9wZg/jvg3cDHQGjp0J3KOqKVX9B7DO/TyDwWAwDAKlFBa1wFuB95vd\nYz4icgxwgKo+tKfXutdfLCIrRWTlzJkzFTA/5sf8mB/zs2c/fWLQHNwiYgE/AK7a289Q1VtVda6q\nzk0kEv03OYPBYDDkUTKfBdAIHBB4P8U95lEJHAU8JSIAk4BlInJGH641GAwGwwBSSs1iBTBDRKaL\nSBTHYb3MO6mqO1R1gqpOU9VpwHLgDFVd6Y6bLyIxEZkOzACeL+FcDQaDwdALJdMsVDUjIpcBjwAh\n4HZVXSMii4CVqrqsl2vXiMi9wKtABrjUREIZDAbD4CHDpUT53LlzdeXKlYM9DYPBMExJp9Ns3ryZ\nzs7O4oOHIPF4nClTphCJRApPSV+uL6XPwmAwGIYNmzdvprKykmnTpuH6WcsGVWX79u1s3ryZ6dOn\n79VnmHIfBoPB0Ac6OzsZP3582QkKABFh/Pjx+6QVGWFhMBgMfaQcBYXHvs7dCAuDwWAwFMUIC4PB\nYDAUxQgLg8FgKAH1DY2cuPgJpi98iBMXP0F9Q//lFT/11FNcdNFF/fZ5fcFEQxkMBkM/U9/QyNUP\nrCaZdtLDGluTXP3AagDq5nQrc7fPbNmyhdNOO81/v3r1ajZs2MCBBx7Yb99hhIXBYDDsITc8uIZX\nt+zc7fmGTa10Ze28Y8l0lq/f/zJ3P7+px2uOnDyG6z45c6/mM3nyZFatWgXAkiVLePrpp/tVUIAR\nFgaDwdDvFAqKYsf7yvHHH08qlaK9vZ3m5mZmz54NwM0338y8efN45plnuO222/jrX/+6T9/TE0ZY\nGAwGwx5STAM4cfETNLYmux2vrU6w9F/ft9ff+9xzTo+4p556ijvuuIM77rjDP/f222/zhS98gWXL\nljF69Oi9/o7dYRzcBoPB0M8smHcYiUgo71giEmLBvMNK8n3pdJqzzz6bm2++mUMPPbQk32GEhcFg\nMPQzdXNquemsWdRWJxAcjeKms2aVxLkN8Oyzz7Jy5Uquu+46Zs+ezezZs9myZUu/focpJGgwGAx9\n4LXXXuOII44Y7GnsE7v5P/QptdtoFgaDwWAoihEWBoPBYCiKERYGg8FgKIoRFgaDwWAoihEWBoPB\nYChKSYWFiJwqImtFZJ2ILOzh/JdFZLWIrBKRv4rIke7xaSKSdI+vEpGfl3KeBoPBYOidkmVwi0gI\nWAJ8DNgMrBCRZar6amDYXar6c3f8GcAPgFPdc+tVdXap5mcwGAyGvlNKzeI4YJ2qblDVLuAe4Mzg\nAFUNVuIaBQyPpA+DwWAYZpRSWNQCbwXeb3aP5SEil4rIeuB7wFcCp6aLSIOIPC0iHyjhPA0Gg6Gs\nOO200/xM7aqqKn7zm9+U/DsHvZCgqi4BlojIecA1wIXA28BUVd0uIu8F6kVkZoEmgohcDFwMMHXq\n1AGeucFgGLH8aSG8s7p/P3PSLPjE4j4NffjhhwF44YUX+NznPkddXV3/zqUHSqlZNAIHBN5PcY/t\njnuAOgBVTanqdvf1C8B6oFt1LFW9VVXnqurcmpqafpu4wWAwDHW2bdvGZz/7We666y6qqqpK/n2l\n1CxWADNEZDqOkJgPnBccICIzVPUN9+3pwBvu8RqgWVWzInIQMAPYUMK5GgwGQ9/powbQnyxZsoTb\nbrsNgIceeogLL7yQa6+9lqOOOmpAvr9kwkJVMyJyGfAIEAJuV9U1IrIIWKmqy4DLROSjQBpowTFB\nAZwMLBKRNGADX1bV5lLN1WAwGIY6l156KZdeeikACxYs4Oijj2b+/PkD9v2m6qzBYDD0gaFUdVZE\nmDlzJuGws99ftGgRZ5xxRtHr9qXq7KA7uA0Gg8GwZwzGJt+U+zAYDAZDUYywMBgMBkNRjLAwGAyG\nPlLOPt59nbsRFgaDwdAH4vE427dvL0uBoaps376deDy+159hHNwGg8HQB6ZMmcLmzZtpamoa7Kns\nFfF4nClTpuz19UZYGAwGQx+IRCJMnz59sKcxaBgzlMFgMBiKYoSFwWAwGIpihIXBYDAYimKEhcFg\nMBiKYoSFwWAwGIpihIXBYDAYimKEhcFgMBiKYoSFwWAwGIpihIXBYDAYimKEhcFgMBiKYoSFwWAw\nGIpihIXBYDAYilJSYSEip4rIWhFZJyILezj/ZRFZLSKrROSvInJk4NzV7nVrRWReKedpMBgMht4p\nmbAQkRCwBPgEcCRwblAYuNylqrNUdTbwPeAH7rVHAvOBmcCpwE/dzzMYDAbDIFBKzeI4YJ2qblDV\nLuAe4MzgAFXdGXg7CvC6ipwJ3KOqKVX9B7DO/TyDwWAwDAKl7GdRC7wVeL8ZOL5wkIhcClwJRIEP\nB65dXnBtbQ/XXgxcDDB16tR+mbTBYDAYujPoDm5VXaKqBwP/AVyzh9feqqpzVXVuTU1NaSZoMBgM\nhpIKi0bggMD7Ke6x3XEPULeX1xoMBoOhhJRSWKwAZojIdBGJ4jislwUHiMiMwNvTgTfc18uA+SIS\nE5HpwAzg+RLO1WAwGAy9UDKfhapmROQy4BEgBNyuqmtEZBGwUlWXAZeJyEeBNNACXOheu0ZE7gVe\nBTLApaqaLdVcDQaDwdA7oqrFR5UBc+fO1ZUrVw72NAwGg6HckL4MGnQHt8FgMBiGPkZYGAwGg6Eo\nRlgYDAaDoShGWBgMBoOhKEZYGAwGg6EoRlgYDAaDoShGWBgMBoOhKEZYGAwGg6Eopaw6azCMeK6p\nX83dz71FVpWQCOcefwA31s0a7GkZDHuMERYGQ4m4pn41dy7f5L/PqvrvjcAwlBvGDGUwlIi7n3Pa\nuZwdeor3ytpuxw2GcsJoFgZDiciqUkkHt0RuJaMWh6Tu9I8bDOWG0SwMhhIREiFCBoCw2HnHDYZy\nwwgLg6FEnHv8AVh01yLOPf6AHkYbDEMbY4YyGErEjXWzCO96x2/pZaKhDOWMERYGQwm55IMH+cJi\n/U2nDe5kDIZ9wJihDIYSYmczgz0Fg6FfMMLCYCghapvIJ8PwoKTCQkROFZG1IrJORBb2cP5KEXlV\nRF4WkcdF5MDAuayIrHJ/lpVyngZDqTCt4w3DhZL5LEQkBCwBPgZsBlaIyDJVfTUwrAGYq6odInIJ\n8D3gHPdcUlVnl2p+BsNAYMxQhuFCKTWL44B1qrpBVbuAe4AzgwNU9UlV7XDfLgemlHA+BsOAYxsz\nlGGYUEphUQsE6xpsdo/tji8Afwq8j4vIShFZLiJ1PV0gIhe7Y1Y2NTXt+4wNhn5GbaNZGIYHQyJ0\nVkQuAOYCHwwcPlBVG0XkIOAJEVmtquuD16nqrcCtAHPnzjVbOMOQQ23jszAMD0qpWTQCwVTVKe6x\nPETko8A3gTNUNeUdV9VG998NwFPAnBLO1WAoCXbWCAvD8KCUwmIFMENEpotIFJgP5EU1icgc4Bc4\nguLdwPGxIhJzX08ATgSCjnGDoSxQtYsPMhjKgJKZoVQ1IyKXAY8AIeB2VV0jIouAlaq6DLgFGA3c\nJ05xtU2qegZwBPALEbFxBNrigigqg6EsMJqFYbhQUp+Fqj4MPFxw7NrA64/u5rpnAVNAx1D2GJ+F\nYbhgMrgNhhJizFCG4YIRFgZDCTGahWG4YISFwVBC1DaahWF4YISFwVBC7EBSnnF2G8oZIywMhlIS\n0Cyypk6UoYwxwsJgKCFBn4URFoZyxggLg6GEBH0WpgKtoZwxwsJgKCGqOQGRNT4LQxljhIXBUCLq\nGxq59alc7cvTfvgU9Q3dyqMZDGWBERYGQwmob2hkwX0v0ZXJaRbtnSkW3P+SERiGssQIC4OhBNzy\nyFrStmKR81mEsUlnlVseWTuIMzMY9o5ehYWIHCsikwLv/0VE/iAiPxGRcaWfnsFQnmxpTQJgkWuz\n4gkO75zBUE4U0yx+AXQBiMjJwGLgv4EduE2HDAZDdyZXJwAIBTQLcQWHd85gKCeKCYuQqja7r88B\nblXV36nqt4BDSjs1g6F8WTDvMCKW5GkWIWwiIWHBvMMGcWYGw95RrER5SETC6sT/fQS4eA+uNZQ5\n9Q2N3PLIWra0JqlKRBCB1o40k6sTLJh3GHVzemupPrLxfjd/+/1f/WNV8RALznyP+b31AXPvDT2K\nLfh3A0+LyDYgCfwFQEQOwTFFGYYp9Q2NXP3AapJpJzegNZn2zzW2Jrn6gdUA5qHthbo5tdS+WQur\nnPe3/csx1B5kfl/FMPfe0KRXM5Sqfge4CrgDOElVPZ3aAv69tFMzDCa3PLKWinQza2KfY468AcBZ\n1v/xIctZ+ZLprB/VU9/QyImLn2D6woc4cfETJjQ0gGouEc82FWjz2N19c8sja0mms3w5tIyZshGA\nq8L3cmfkO0D+vWcYOIqaklR1uYicAnzObX26RlWfLPnMDIPKltYkH7P+zihJcWm4ni+mF/CD6M8B\nmNZ5lz+mcBfY2JrkiqWruHzpKmqNySCvkOBI721RaFra1ZUhnXX2n0GNYUtrEsFmYeQeFnIP0zrv\n4t/D9c6HuEqGiSgbeIqFztaKyHPA9cBB7s/1IvK8iBRdAUTkVBFZKyLrRGRhD+evFJFXReRlEXlc\nRA4MnLtQRN5wfy7c4/+ZYZ+YXJ2gjQoAKqXnB1OBq+59iWQ6ywGylf3Z7h+H3AIwojWNQKc8HcHl\nPrxNRWNrEsUxLaWzylGygQSdgKMxXL50FZYIcScIc7eYiLKBp1g01H8BP1PVD6rqle7PB93jP+3t\nQhEJAUuATwBHAueKyJEFwxqAuap6NHA/8D332nHAdcDxwHHAdSIyds/+a4Z94ZTDa8hoCIAxdOx2\nXNa1TP4ldgVPxa4EoIZW5oeeAIzJIKhN2CO4xapnWhpDO58NPQooETL8MXYNd0dvzBubVaWCVK+f\nd8rhNSWcraEnigmLI1X1jsKDqvrfwOFFrj0OWKeqG1S1C7gHOLPgc55UVW8lWg5McV/PAx5T1WZV\nbQEeA04t8n2GfuTJ15uIiFOqolI6iJBfMfVQeYuvhZcCSg2tAMTEsRHcGv0BiyO/ZD+cqOvG1uTI\n9WcENYsR6LPw/BKNrtnoW+E7+XbkDt5nvUrCFQizrQ2As8m4IfxrImRIyO40C2dz8uTrTaWeuqGA\nYsKix/MiYgGhItfWAm8F3m92j+2OLwB/2pNrReRiEVkpIiubmszN059saU0SdQVEhAyVBdrFzyM/\n5LLwH5hIK6MKzFSHiCMQEuIsBjW0INgj0iwV1CxGms/CMz21tjYz2r1/4q4Q2I+Wbqam70R+xYXh\nxzjResUXJIXEXKeF8VkMPMWExR9F5DYRGeUdcF//HHi4vyYhIhcAc4Fb9uQ6Vb1VVeeq6tyaGqOW\n9ieTqxO+NmGhVEq+sMi4e4X9ZbsvVDw8H0clSappY0X8Ur4evgcYgWapoGYxwsxQnulpRezfeCX+\nRQBadTQAY6XNFxweE6UFgA6N5QmLEDkhG3MFjPFZDDzFhMXXcfIp3hSRF0TkRWAjsBP4WpFrG4ED\nAu+nuMfyEJGPAt8EzlDV1J5caygdC+YdRtxyFjdB/R2dRwuVANTKNqIF5zxGS9IXMp8K5ZLTRtSu\nUIOaxcgSFt7fuUJyC/8unEV+FJ3dtIextAMQlUzeNcF7L0aaRCRksuAHgV5DZ1U1DXxNRILlPdYH\n/Ay9sQKYISLTcRb6+cB5wQEiMgen/tSpqvpu4NQjwHcDTu2PA1f34TsN/Ujc8jQLmwhBE4rSqVEA\nqmQXUe25A9wYOoi7D/ooN+IFoCoRKc2EhyI6cjvlVVdEaOnI30h4WkJM0t3MUJ6PLEGKdMDKHQuM\nq5Au5r23dmSHYw8SRfMsRGQ8ziLvObRfE5G7VXV7b9epakZELsNZ+EPA7aq6RkQWAStVdRmO2Wk0\ncJ+bw7FJVc9Q1WYR+TaOwAFYFKhRZSgR9Q2NXL9sjZ8x++lQGizHDBV0cEfIoggAFaSISv4imNII\nMUkzmiQVrpAYLTlh4fypRwZBbWKkmaE601nCeSZK9QVEnC5/I5E769wYCVKEA0tTULOIkubO5Zu4\nc/kmQiKce/wB3Fg3q3T/CYNPr8JCRI4AnsBZ8BsAAY4FviEiH1bV13u7XlUfpsC3oarXBl5/tJdr\nbwduL/YfMPQPXrOetJ0rfBf0WQSFRZS0n0uRIJVnhhJsMljEgFGS7DEEsrWjZ7PVsCRghmKEmaGS\naZtqcibHUXT65qU4XX4ABDj3jUdCuhDN3YdelJ13nUdWlTuXbwIwAmMAKOaz+DbwVVW9SFV/rKo/\nUtULcUp9fKf00zMMFF6zng9aL/HH6DeYLm/7AkJQwpJb9CJk/HMVkspzcMfp8neIcbp6jGoZUc7J\noBlqhGkWkK8VOGbJnLAImpeiZALaameethr8jAgZzg49xe+i1zGOnQDc/VwwcNJQKooJi1mqem/h\nQVX9HXBUaaZkGAw8Z+T5oT9zlLWRo2V9QLOw8wRCNBAHX0FnnmaRyBMW6W6ahcDIck7aIzfPYmxF\nJE8rSEiKhCsgEtLlv4Z8v0SCrrz7LSgswmS5PPw73mu9wQw3RDsb0EIMpaOYsNi1l+cMZYa320/h\nOJ/HSAdR1xlpoX5kFEA04JysIJVnonLszc51ccmZGjLq3GoV0RBXLF01chL08pLyRoaD20vEa+lI\n5y30jqbp3DcxuvJCZ2Ok/U6Ccck3bVaGAsJCslS7UVPevRUaSU6wQaSYg3uiiFzZw3EBTGLDMGLB\nvMNYcN9LdHnCgg6/y5uFTSIUjHVP+zu/Csl3cCck5S8QB1eH2NnuvM5iEQkJu7pyBQdHRKnpvNDZ\n4b8DLiwsGfQx7FcBVZoGG6rCGeJdAWER2IBEyeZpFpVWTjuNkvHPeePPPT4YZW8oFcU0i9uAyh5+\nRgO/LO3UDANJ3ZxazjnuAL/15xjZ5YcyRship3MPdlxyPosEqW52aUucz2hta+OICU4IpCKks8ox\n8neOk9eAEZKgF/RT6PDP4L7hwTUk045fYTw78sxLV31oCjF1ouIk21lghkr77yNk8jSLcCYXqR8m\nS8T1n1XQiQXMPXBcKf9LBpdieRY37O6ciFze/9MxlJpgmejCrmNPvt7E+9wHtopdtOJk21qieZEr\nlVElqm7+hHTm7QKrpd1/HbZTvLV1e15hmAdi1wO5MueNwz1BbwT5LOobGmnpSHOwbOGWyK08bs3h\n9myupNsrm7by/mwnWG7orAT9FCnfDxYh429UwEnu9AgH8n0S0oWNE5zh3cO93d+GfaOYZtEbPZmn\nDEOYwjLRhbWaGluTvpbg+CwCD2ygNlQ21ZmnWQR9Fp49GVy7tLuTjEsaAv2ovfwLa5ibmyVohhrm\n0VCelniEOOGs+0l+/ae/vd7oC4QEXX5kFMDYwCYjGjBzAowOhN/mR94513sbjmL3t2Hf2BdhMcwf\n8+GHV6tnrrzOpSGnmYxnCqpvaETIRZ6MYVeeEBgdyMCekICY57MoyLOollzcQ5yuvGio8W6oI+Tq\nAA17M36eg3t4m6G8iLr93L9tCDs/kinbmUvKk648QTIucG9EJJN3XVBYBItWemYrIadRhNJtfDd8\nm1/4ckSYOgeIfREWw/0xH3Z4D/P9sUUsiNzrZ9c2tia55ZG1KLkEqDEFZcmDpoBkZzKnWRQ4uKsC\nO8Q4XVQEMrcnSy7pv7ceGcOKESQsvIg6b3GPUhANJV2+RhkjP3R2guSERSzgxIb8e68qEITpJfgp\nzkaosTXJxeE/cl74Sc4NPe6PG1G1yEpIsU55bSKys4efNmDyAM3R0E94D3NWHaVw/8Di7anyOc2i\nI99uXGAK8LSJwjyLoBkqLvl5FsHvG+MWGKwe7nWi8qrODu/91YJ5hyHkFvcq2ZWXZzGGDkJu8EOc\ndJ7PYnxAWASTPgFGBe69qoDmGkz49O7fGnYAuYKFMMKSQEtIMQd35UBNxFB6Fsw7jCuWrqKJaibR\nwhTZxlu6X94YT1gUFggM7u7idPkP/ShSBQ5u52HOqnTL4N5fcuW9xrg7xNZkmmkLH2JsRYTrPjlz\n2Dkjgz6L4RoNdU39au5+7i0/Oc7bWIylLW+zMFbaAEhq1K0N1UWnRohLOt8MRYZoQMhU7kazKCyN\nD7kAiy53aRtxSaAlZF/MUIYyo25OLYpT6A/yK8EeIW8yie1+qOMYdu3Wyejt7jo05sfHt6rT8mQs\nzoKwg1GuGapnzSK4QwRo6Uiz4P6Xhp8zMs8MNfw0i2vqV3Pn8k15WdTe4h4SpUZa/ePevdFMJRHJ\nUiW72M4YIKdZJDVKtCDzP3ifjpGgsEgzSzYw3tUmIOdI9e5XZZjn8QwgRliMMMZWRLAD1T09/hS7\nmj/HFvhmg5Cob18GJ0TW/wzX1PSuVgOOL2IXcdIa8oXATh3l2qh3p1l0ECbDBaHHqMFxiKazOvyc\nkcM8z8Kry3SC9SoftF4C8jcWk9y/eVKjvmbR4hosxrOTZvf1FNkGwDs6lohk8wTE7nwWUcnwYOwa\nHox90z/m3dve9cPezDmAFC1Rbhg+1Dc00t6ZwQ47ewTPZuyFxY6WTl/rgPzdfyVJMmoRFtt/6Nfp\nZKaxlYOtLaQ1TFJivs9iB6OYzHaykqJZRzNO2pkkzaQ0jIUyRnZxfuhxboj8hsNlE9dkvgA4tucj\nvvUnUhkbWyn7MtQyTBzcheXrKyIWsUjI1yjuid4IwImdP85b3PeXZro0xC7ijHPvm47QGFAYJzt5\nxZ4OwAzZTErDvKUTqZQkMekirSEikmU0nbRpgkpJ+r4uyGkqkwObEC8Pw9vc7OrKUN/QaLSLfsBo\nFiMIr7JssG/Ax6yVfstLcPwRO9VxCI6hg10aAxwnYwdxUhoJCAvnAZwi2+giTJKobzPORMcQp4tR\nkmKH21Vvf7bTQiU7qaCKXcwPPQnAUdY/8uaZTNt+SK1Xhvqa+tWl+JWUnjzNojzzLLzy9Z6gAOhI\n235jo6AZ6BOh5xlNkkYdD8AkmkkRpZMo4y1n8xEa7ZwbJ+00u/dGWGze0omkcMxQo+n0k0JH00EH\nzn3oaRZtmmCKNOXmGP0WNYG+3p5mMSy11UHCCIsRhBdCmDNDdfGZ0NN5Y2Kk/T7JVbKLdjeqpFKS\ndBGmnTjj3B3dBt3fj6zqIkJSY1S5O7+3U3FidBHTFC3u502SZlq0kp1awVHWPzjCcpK3DhSnSeIE\ndjAZxxxxjPyd+6LX+1pPuZahztcsylNYeJuM91uvcFfkRt88eaRsJEwmb9E+PfQcoyXJWzoRcDSL\nFBE6NcoYde6bNS05g0arjsZ276E3dT9ShImSYZR0+vfNaEnSqdE8M2eLjqbWNV0BzLbW8x5rg19p\nIBiya0Jn+wcjLEYQhSGEFZKffQ3ODs/b0VWxy7cpj5V2MoSwEU4OObv8nTqKLToBwNcsPHboKKKS\npVKSbHc/IypZmrWSnYxitrUBgPsyJzNW2qmkg6XRRTwb/woxuvif6E0ca/19GJShtv3FsFwzuL3F\n9rzQ47w/9CqnWs9znLzGw7Fv8K+hP/r5Mw9lj2OOtY4pso1NtiMsYpJ2hAVRP0vb0yYAdhH3I5c2\n6UQ6iVElu9hfmv0+7+OknXYSZAj5mkULlcQKOjTuJy2+jywYTWVCZ/sHIyxGEAvmHUbEEl9AxEkx\nUVrZpmNYbU/zx3mahSVKs47xj6c1zNuueQEgSYx3cNqkpwnTQdw/t5OKbp8HzkO+w42cyqrwF/to\nAKbKVg623gbgUNnsR2J5Zq1yLRcg6nQO9F6XI95iu9P9u02SZt5nvQrAyaGXfWHxlD3bv+Zdqn3/\nV6c6ZiiPlkBE/i5N+MUrt+pYOjTGRDeCarPmClu3UUGakF+k0ruHAJZmPgQ4VQFylWsdYREJiQmd\n7SdKKixE5FQRWSsi60RkYQ/nTxaRF0UkIyKfLjiXFZFV7s+yUs5zuHNN/WoOvvphLl+6irStfsRT\ngi7Gyw7+nD2GB7If8Me3BHZ+rYzyd8ZpQlyZvsQ/16ExfwFJaYRtWuWf2xYQMsHPa9HRviDZylg2\n6CQAPuRG0gBMkSZ/AfFMXiKUZVitIyycnXO5OrhPOdxZtL1NRo3s8KOcDpYt7C/b6dAYrwY2HG1a\nQYuroTomypywaA4Ii+aAhrCDUb5v4l2t5tbM6Xmfl3Z/j50aIemOA/jPbB1NWkUNO3wzlNeLJZ1V\nLl+6ioOvfrh8/V5DhJIJCxEJAUuATwBHAueKyJEFwzYBFwF39fARSVWd7f6cUap5Dne8OPhx2uwX\nXvPsuQlSjKGDnYzKe/haA7u2FFHfb5EmzEZ3cQfoIEabey5JzA+l7dIQTTrWH9cUECLNOD4LgEad\nwCbXtn1m6Bl/TFBYeM50WylLR6Voloy4ZdrLVLN48nXHJ+HlT9RIq69NTJCdvNd6g7d1HJtdkyQ4\nmqWnQaQI0xbQNN/RXEnxoOBo0wqSbkDFq/aBvuBwPs8xQ4Fz33lNugA6NcZOraBSOvwSIp5mUUMr\nCTrLP1BiCFBKzeI4YJ2qblDVLuAe4MzgAFXdqKovA+X5FJUBnmN4RfxS7os6Fec9YVEpSRLSRZsm\n8nZ+rQFNIK1hXxPIEPJ3yQAdxH3Nop04W10B0UWEVnICZ0tgEWlxfRbgCIudjGanVnCo1UirjqJD\nYxwgTX5lWYbDAAAgAElEQVSGeLDkeWNrsvwedlWy3u+sDIXFNfWr/VIaXnG+idLKJGn2NwdzrHU0\n6gR2BP7mTVrtC4J2TfgbkC4N5fkstgc00J1U+AKigxgZzd1rQc2ig7jfpAsgSZQ2ElSS9HOHomQI\nkWVF/N+4NfIDf2y5BkoMBUopLGqB4F9ms3usr8RFZKWILBeRup4GiMjF7piVTU1NPQ0Z8WRVqXZN\nObOsjUTIVfT0Kr/uZBSdgV1cS8DHkCZEm+Y0iyAdmtMsdmmCd3EWjwpSeTbldwJaRouOZpc6vg3P\nbLXF9YMst4/kbR3HEdab/vhRdBIm42eWl9vuUDRL1nvMykxYnH/b37hz+Sa/NlOlG+lWQyv7SzOP\nZ+eQVme3/7pOJehZatIq3/zYRgU7XJPUDrdAjIeXwQ2uZuHeh0miefdbGxVk3O/q0BhdAUGSIkqb\nVlAlu9xS+E7l2uMtp8nWB0Kv+GPLN1Bi8BnKDu4DVXUucB7wIxE5uHCAqt6qqnNVdW5Njenyujum\nyzv+6+Ot1/xOY54jsU0TeZFMrXnCImdCSAe7GOGYnjyhoMC7rlCwRP2IKiDvdQuVfjLgGzoFgJds\n5097X/ZktupYjpScsBhNktsjt/BU7Eq8Qse/Xb5pT38Fg0J9QyPN7Z2k3EXu4Ze3lE3v8fqGRp5Z\n38w0eZs18S9wdugpPzN7urWVMdLBm7ofv8+eBMCz9kwAtrraRlCz2Kmj/Huq0H8RDJjYwSi/JMzb\nOj7vfmvToBkq7puhUhomQ4g2KvJKi0RJ8xGrAcDPFTLsG6UUFo1AsDnuFPdYn1DVRvffDcBTwJz+\nnNxIojKQ9XqMvOG/nuiW2WgjZyuG/MW9izBtro8h7e7mbk7PB6CdBM3uzrBSkr5ZApzFwiPPoamV\n/DrzCb6dPp/7sycD8K3M5/in1I08br+XrYxlVKCeVLXs4uTQavaXZg4SJ1qqHPaGXiMe1PZ3xCHs\nsmnI4/mHPma9AMCFoUfz7iNwNMJrMp/nn1I3+pFQZ6Ru5NvpC9jCeF+z2EXcNz0JmhcZFdQy3tFx\nPJ19DwAPZD/QTbPwhMcuzQkL77N2aoV/P4NjhvKSPXsqOGjYc0opLFYAM0RkuohEgflAn6KaRGSs\niMTc1xOAE4FXSzbTYUxtdSKvb8CRARNP1NUw2qjI0yzaNU6Xu8A5mkW+Gepn2TOY1nkXNpavWVTT\nzkadxHp7f25Oz88Lo20JmBpaqWQbVfwqezpZ9+HvIsIrehAAWwPOzxYdzUGyxX8/3RUWMPQjo7xG\nUxa2vyMW1zVXDg15PD+Fl3AXIssYkr7mAM7iHvzbAWxlHL/KngaIH8hQSQcv2IcCsM6uZRdx/i87\ny98seHQQ53k9gmmdd7FBJ+cLC63Ic3B7uRmesGijwr+fwdEsvM1FRLJ+JFetybnYa0omLFQ1A1wG\nPAK8BtyrqmtEZJGInAEgIseKyGbgbOAXIrLGvfwIYKWIvAQ8CSxWVSMs9gAvXNZpleoIi6xKnqru\nsTNgKwb88gzgaBO+ZlFghoKcCak+exIdxPlI1/f5WdYJXjsrdT1npa4H4Mnse0hryC8DsTuCi9Hf\ndYqfewFwQCBTeKgvtl4im4X6i5wV0ImGelZxSBz/wyTXrzVNthKTNBvsXBubtxnX47UeT7s5NM/p\nEWzQyXy169+4Kn0JIPxL+mq+lv6yMy57NDu0otv1eWYoChzcbg6Hd6wtcH1Kw0Qlwxh2kXK14YSb\ndd7YmmTOokeH/GZjKFLSQoKq+jDwcMGxawOvV+CYpwqvexYoz8pxQwAvXNbD63u8lbHU4AiLlIb9\n+PadVKCac046ed1BzaLCfx1ExPFTTOvsKfIZXtRD/ddfSl9FCBsN7E9qqxOccngNT77e5O9kvVBa\nyE/egvxOe95DP1R7YEyuTtDYmsTC9jWooLAYilnFXmvSxoAg28/Np/DuofW6P+9zlfygFtjT33Kd\nTmFm56/Y5WqZf7BP6vF7L0x3S8Fyyd2TbZrw779dmtMssurcT22BZkc7GcUYOohJhk12DVOliQpS\n7HTNq145fDDly/eEoezgNuwlXnjguaHHWRpd5BdVe1erqRGn6NuOgF+irUCz6CTqh8h2Efajo7KB\nnZ4AP/zMbBKR7tpGT2QI59mnLzhhKs8s/DA31s3imYUf9s0DK20n2/ZF+xA/vwOcBaJQK2rpSHPV\nfUOzB8aCeYeRiISwUH+H7AmLRCQ05LKKPR9LY4HGs1/B73y1a3LarBPywleDf8sfnTObiOUs9E7H\nur7n3yciIX50zuxux1uo9H0/SeL+veRFmgXzOHZqhV+d9h1X+5ltrefx6FW+WdMUGNxzjLAYhnjh\ngd8K38nx1uuc4JZmeFfH+jvEYHx7O/nRUCmN+AtcmrA/VgLpMJOrE9TNqeWms2ZRW51AcHaXPzpn\nNhsXn86PzpntH09ELNy1g5AIF5wwtVvJcW/x3MFoTk99l/O6vkm75vweb2gtE2llijQRJ8UM2ez8\nX23lhgfXMNTwfjdhsf0FTbCprU5w01mzhtyO1vOxxOhiqmylhhaqaaOGVl4JZGZv1zGclPox/+ya\nF3uibk4tt5z9nrxeEhURi7EVkV7vk+DvxjODeWzTKv/3uCvgs/AqKAfNUNvIJYF6uT/fivwPB1tv\nc5z1un9uqJsChxqmn8UwJCRCVtU3enghssEIpXe1miNwTFVZQnRSYIbSEIgjLLwdZPDx9UpA1M2p\n7XHh293x3VE3p5ZvPPAyHWmbNToNwDdf7NQEb+t4Pmq9yF9jX/WvOT31HdbodL9U9lCjbk4trz4c\nAkKQhg8cMo6vf+7Dgz2tHvEWzgXhpXwx/Ke8c2uy0zjK2gg4i3KwZhM4DbUK6evff3fjTjhoLM+s\nz/WpaCNByN2sdGjOZ+EJkGAtsmDFAC9bPOrXQ8sFewxFU+BQxmgWw5ATDnJ2U96DMZY2OjVCeyBC\nyQtz9UqMB30JKSK+fThNyBcymwL9ur0SEP3Jd8862jdfgJPoB04+R5NW+/khHodIzvw0FE1RABY2\ntjX0M7ir3QXfEwpBXtUD/ddB3wA4hfqu++TMfp/Pxu2Fu34h4vrYgtFQtnvfBgVEU14It3Pca4o0\nTnK9vr0Nj6FvGM1iGOI9aF7JjHHSRpJYXpa2l20dtDt7dJEzQ3VpmOf0CL7YdRX/50a3QGlUeG+H\n6XVk84Rbh8bycjg8plvv+IVibnlk7ZAz7YBTSNB2a0MNZWHhJTZ3BhLmPNbYOWER3MHXVidYMO+w\nkvzevfvrpNSP/Da+Xr2nYFKeZ4YKCoigiTVYbh9gvFvNAEqz4RnOGM1iGOI8aLnIm0pJ0kmUjkDi\nnWfL7WkRzhLyo6E8R/ef7ffmCZZSqfB1c2oZFfOiXjzNIk4T3ecZTMJqbE0OSe3C8fNYjgY3RIVF\nfUOj3wVvYg+h1W8Gikd6vgERx6ldKgHt3V+bdaLvVPcqye7SeCC/whkXrEvVUzUCy9845TQL47PY\nM4ywGIZUV0QYQ362bVKjeQ/RcvtI/pB9P/8v8xn/2K8z8/h+2qkU75V2CBZ98yh1NI/3EL/l2sbH\nyC422s6CtV0rObfrm7xuH8CEwIMPsGAIRkaJ2qhYjrlkCAqL+oZGrrx3lf9+orTwh+z7+VTqBv9Y\nXv0mV7ModYklL5osiHf/dhLleftwHssew+2ZT7hnhV9n5nFj+vy8TU1h+LXXBxyMz2JPMWaoYUAw\nPl5wdIppBQtpiigpyWkWbVrBV9OX5Y25IXOh/3p/N75+nZ2/cyyl6cHDy1F4QQ/lefswbsuczgo9\njEez7+Xe7If4mz2TbTqGCbKDGbKZRp1AB3HStnL9sjVDyhwlKCqWYy4Zgm1Vv37/S9gKB0sjKSJM\nkJ2styfToDP4ReZ0tmkVNha/zHyCLiK+j6DUeH/DYN7H266zupp2OonxpfTX8q7x7t9z3N7uAP98\n0lFOLQmX8eSei8bWJNMWPsTYisiQzdcZShhhUeZ48fHJtKOiJ+gkRcRvGuSRJEqHHazJ091XEeRL\n6as4P/Q4TYEwxB+dM3tAHqgF8w7jiqWryBLiM13X+ccvTl/lv95GFXXWszwW+zrP2YdzTpeT6+mZ\nU4YKlmZRhCwWokOr+VF9QyNdWSVGF4/HFvjHve6HN2XO94/dmPls3rXBsNhS4UVKTVv4EADfy8yn\nRnbwuN17mTivEi7AD/66jfmBOoLjZSegJEiRJkyGsEnS6yPGDFXG1Dc0ctW9L5FMZ6minUXhX/Nq\n/PP8NPJjPhxq8COdAJIaKyjpkXvYQyLMmDgqL7b9L/bRfDl9BV7AbHUiMmAPUt2cWs4/YWqvqVzB\nrnzHW68zNrBjHEqmKEezCDmaxRAzQ3n5Kadaz+cdfzdQUr4nIpZw/Rn9HwG1O7yEzU26H+d0Xetn\nYnuMrYhw4sHj/Ps3aIYqLCMyVtqZLet5Lf55no19ha+EHgCcJL0r7l01pO6doYYRFmVKfUMjC+5/\nyU/AuzJ8H/8SfgyAeaGVfCr0V560Z/u1mDqJ5nUe84RFbXWC9TedxmNXfoj1N53Gj87pnpWdiIQG\ndHEAuLFuFj8MJGxVJyJ58fxBYQHw8dAL/uvrlw2dJD0LG0Rc883Qqpfr5aecGXo27/jWAmHh/e69\nxLlbzn7PgO7Ae/JfeJneGxefTsO1H+e3X3of6286jdrqhB9WC8597rUFvj97Mu0a58rwfYDjzL8y\ncj9j3GgrVbh86aqy6pcykBgzVJlyw4NrSGeVMbTzr+E/5hXZA9ifZu7TDzJNt4I4ZqhgGXIv2qmw\nvEPQVrylNcnkAfBR7I6eErbmLHqUlo60n6XbpFV0aIzTreUszZ4COKao+obGIWFSEGxUQtgydDQL\nz8cFMJoOTrRe4b7MyZwd/j8gX1gMlOmxN/bknmxsTTLDCprIhE6iVJBim1axTmvzsrgBTg2tYAwd\n/DLr9Pz+7fJNzD1w3KD/v4caRliUKd6u8HuR2zg1tMIvKe5hidKslb42kSKaFyPvmZcKyyrAnmdf\nDySt7v97g70/AC/aM2jUCZwXepwQWb9+1eVLV7HyzeZuZUUGGkud4omKIPbg+yw8jTSddbSc2dZ6\nYpKh3j6Rs3GERUvAzDNU7oO+3pMhkTzNAnIbow6NsVXHMttan3f+e5HbAPizfQwbdX+UoZu3M5gY\nM1SZM8vaAJDnj/Bo0Ur/eKdG/f4CQcqtzaQX7viiHsq/dl3OtzOfZZV9MHFJc7L1Mj+MLPH9F3cu\n3zToNmjBBrGwCTEUzFCeRvo+aw3XhX/D4eKUfHnNPpBzu77Jaanv5mXzlxtZ1TwNGhy/ETjadaGJ\nLchRstF/bXIwumM0izKlImLRkbb9B6GqoIsZODkS3oPTRoKdgcQlj3JrBuNFSinwiH0ckMvs/XX0\nFgDesKfw0+yZwODvEC0/z2JomKE8jfQ/I//JBNlJs46mScfQzBj+Zuf7pQYi4qm/qa1O0LQjP4HT\ny/JOEutVWCQCHRqryvD/XmrKdwsxQqlvaGT2DY/SkXYWnt4ihpp1jK+St2uix3DZcquP01Ok1D90\nUl6F2mA3wEKfzEAzNKOhlEo3aXOctPN3+4BuIwY64qm/OOXwmrw6UZArNtihMbbSi7AgJyxak2nT\nJKkAIyzKBE9IXL50FeHkNiJk+Eb4t37yXE80a6UTjQNub4juoqUc6+N4kVLezlex8naMB8o7/uue\nfDIDSTAaSoaAsBCcwpJe4yuABj0kb8xgRDz1F0++3uT3uki6da4ynrAg3qtm8ZXw7/mk9SxjaKeC\nTlo60ly+dJURGi7GDFUG5CfeKSvjl/B09mg+GHo5b9z/Zo/l5sx8now5yWutjPbFg1fu+9Opa3mH\nXGvTcrXNFiZs/TRzJqeHltPKaD5qveiPy6oOamSU5foshoJmUd/QiAK1sg2An2f+iWOttTyQ/UDe\nuGcWDs0y6n3Bu58/nbqWt92wca/pVEeBGWpa529ZF/ssYXH+LhNkJ/8Z/S/Aae50UuongGO6u/oB\nJ5y2HAVof1FSzUJEThWRtSKyTkS69U4UkZNF5EURyYjIpwvOXSgib7g/FxZeO5LwGtMcIW9yiuXU\n8SkUFOA48Fo155dIEfWrbnq+jZV6eF4/gnKvj+NVNP+dfTKfT3+dtfYBjJEORpETggvuH7yaUTkz\n1OBqFvUNjVy51Ll3PGHxYPZ9/HPXDWzQXF/tcvRTBPHu55V6OI049/l4tx7UZq0p0Cykx8AQgCmy\njfmhJ/z7KJnOjvjOeiUTFiISApYAnwCOBM4VkSMLhm0CLgLuKrh2HHAdcDxwHHCdiPSeVjqM8XZL\nf4pd7TtxeyKp0W5ObC8pL9j0xWMotvfcU847fmree283ebS1gZnyD8DJzh2sbnqeZmGLBQyesPjG\nAy9jo5xkreZA2QrQrYmRJZSlnyJITwl8KXUMKJu1Jq86LeBXr+2JxZFfck34Tv99uWrh/UUpzVDH\nAetUdQOAiNwDnAlut3dAVTe65wqfonnAY6ra7J5/DDgVuLuE8x2yeIX1ipEi6ucZPJN1HvqfZ84g\nSoZ7sx/yx4n7mYOVbNefeHkUdy53QkA94Xh39DsAHNr5G7qIDFo3PV9YIIOmWdQ3NNKRtvmQ9RJ3\nRL8HQLvG8xbO4VJMr6cChGd13cBc6+95ZUBWu61iOzXaa5TIBLdnPZS/Fr6vlFJY1AJvBd5vxtEU\n9vbabnexiFwMXAwwderUwtPDhgXzDmPBfS8VHeep1O/pvJUOr3EQcRZnzvPH1FYnytom3RM31s3i\nydebaGxN+pqFx2xZx/N6BAAnLn5iwAWkFzqrWKWv690DXv0wwO/FDp5W4aySYysiNFz78QGfW6nw\n/FknLn6CxtYka3Q6a7LT/fOzOn/pRwkme9EsIJfQNxy08H2lrKOhVPVWVZ2rqnNrasorBLSv1Dc0\ncv2yNaTt4guNF/2xg9F+W9Qgw/mG9/5fhaGR061cZFRja3LAe15YKEgIFXEbIQ0chfXDpgeixLYF\nusml0oOfWV4KejJJgdOTw4+Y2o3PwsMLu/V8FiM5KqqUwqIRCAZwT3GPlfraYYMXBdXXstue/dXT\nqgsLwN101qyyNzPsDu//5ZnhPKbLO5wTepJj5O8Afs+LgSIXDTXwJcq9bG3v/x+sHxb0V3g5O8ON\nujm13HTWLL8YZW11ggtOmOonogq9+yyAvP4dja1Jrn5g9YgVGKU0Q60AZojIdJyFfj5wXu+X+DwC\nfDfg1P44cHX/T3Fo40VBnW4tp1aK50MkiQ1Ic6KhzrzUYqbJOywI38vhsokvhx8EYFqnE0cxkD0v\nLDcayh4EM1RLR5r9aObmyG1k1CJJjDszH2Gd1ub5sIYzvdWUqm9oJPX73qO/RtHJ9yM/48eZs9ik\n+/kaxkh8vkomLFQ1IyKX4Sz8IeB2VV0jIouAlaq6TESOBX4PjAU+KSI3qOpMVW0WkW+T63G1yHN2\nD0e8KqBbWpNUJSKIOAXzvKVlSfQnffocKzr8/BF7glcCZa1OZa1O5Wx9mo+GGvzzgj3gdY/ETcrT\nQXJwH245rr+w2FSS5B+6P3dkT80bU+7hsntL3Zxanvx970vgR9z7J0GKf0tfDjgaxomLn+CUw2t4\n8vWmQa/OPFCU9MlR1YdV9VBVPVhVv+Meu1ZVl7mvV6jqFFUdparjVXVm4NrbVfUQ9+fXpZznYOKZ\nmhpbkyjOrrclICiKcUXXJTS59ufmVGjEqsgA3z3raD/vAmCjTso7X8Uu//VA/Z5CnhlKBrafhff/\nm8COvONvFYTLlmtZj/5CXW3vZ5lP8ku/n3d32jU/EqqxNcmdyzf5z+1IMFGVtYN7OOCZmg6Vtzha\nnNLJlXRwfujPfXKI/t7+APiF0qID7sAdStTNqeUHn5ntvy8UFsEwyKvuHZjfU0gUrNCAahbBCKjx\nki8s3tApee/LtaxHfxBscrRJJ/JwdvfBml4FhBpa+XToaf/4J6znGO3W2RruiXtGWAwyXqLPo7H/\nYFnsWwBcFv4934nczqdD/8cXQg/1+bOSxAbcgTvUqJtT6zsw1wcykwFqZAc1tFBBJ1lVrixxG007\n6zq0Xc1iIBzc9Q2NXHnvKsZrMwk6mSA7886/qfv5r2urEyNWUNQ3NPJbNzcHnHyL3pzdh8tbXBB6\njLujN/L/Ir+ginYOlHf4WfTH/CDyM3/ccE7cM7WhBoGgj8IS6dZTwovtvjb8P1RK8ZvvHR1Ljexw\nEowYWAfuUGTBvMO4fOkqnrOP4NeZebyt4/hG5G72Zzt3x7/DCvtQzu66HludzOZSLZi27QZeurWh\nZADMUN944GVCmuH5+KU8mz2SdxjHZp3A77Mnsd6enBfdM1zDqPvCLY+sRcFvCNZCpVtss2feF3qV\n94VyeSqj6CSMI/wPlc3+cWVw8nkGAqNZDDCej2JL6y7G0N5j86F2t0lRrIcSHT3x+a4F/GemjjU6\nrT+nWrZ4D6mNxQ2ZC7kn6zj9PxJyCgwea/3dH1vKsNFs1qnsKhLCJjQgZqiOtM17ZB0A7w+9ygR2\nsE2r+H7mM9TbJ+WNHW6L2Z7gaQBLMnUsTs/nafvovLpqxaiQnLDI9cFQqmgftv4LIywGmBseXEMy\nneXzof/lpfjF1NI9JLZCOgGISt/MFk2M5fuZz/iJRmMrRmZ0S5ALTshl9O+kgpRGOD7QezlBZ8nn\noLYjHNQa2KS8yYGy9QdbW/IS8DyCv5+RiFe64w2dws+zZ6BYtNG9k+TuGE0nMbfpVsLd1F0R/h0v\nxS+mmjaS6eyA+cUGCiMsBpD6hka/RtFJluNcO9zaRIRcb4GN8fP493D9bj/jpvS5fL7ra7s9HwkJ\n131y5Ea3eNxYNyuwIApNVOXZ72sCjt9SPdA5zULcpLzSmqG8/0eNtPrHamU72wqaAV1wwtRB700+\n2PSU3R0Mq76k66s8mD1ht9fXx66lPnotAHG3adJ5occBx0QFTnn84aRhGGExgHiREidaq31n2njZ\nmdehqxgr7MN4ynYifgrV5trqBLd8euRGtxRyY90sP5u9SfNbbU6kxX9dqgiWP73sLBKPvb6Nnaks\nbZ1dnLj4iZItHt7/IygsIL8EisCIFxTQPbu7sEnWn+zju1XlLcQSR/h7V3obkEOsLRzi+jGGU4SU\ncXAPIFtak0xmG7+N3uQfG09bnmbREzu1gjFuj+0UUWwsrui6hBd1hj9GKO+mNaWiKhGhNZn2W22+\no2OZJC1MlFY/7aGxNcmcRY/2a9XV+oZGbv7f1/h0yNmx2iqEJevbs6H/fAbX1K/m7ufe8v1fE6WV\nJq3yF6+tOs4fW21MlD7B7O7pbhOti7oW0KqVAH7ASDHsgrK1v4neDOQqBgyXCCmjWQwg1RWRvFh/\ngDAZor0Ii5ft6Ryd+qX/vsMtfPZ7+wO8GcgjGOnlk3eHt2FscReA5+3DAWdBjZL2k9ZaOtL92iTp\nlkfWkk47f1cbIYvlJOjRv7vNa+pXc+fyTWRVmYzT1GgiLXnJd+8EGv4MQuHbssB7fp6y57DKbTPr\nPWtPZGdzdOdtu71We6txzvB5No2wKDH1DY2cuPgJpi18iJaONJWuhuBhoURl96GuhZmjze6iF2Q4\nV5PdV1pdH9EqPRiAp7PvIa0haqSVG8J3sDJ+CWNxfBnprPbbIr6lNekLhwxObSgrEDrbX7vNu59z\nynl8yvoLz8a/wvutV6iRHbyrY3nRdha9F+xD/fE7RnhY9e7oyYfRgvOsCeon5fWEXWQZ9cqDlLvv\nwgiLEhIs5eExlva8MVdEfsfTsSt3+xmFERpeXPioaGhEVJPdV7xd3d3ZD3N+19U8YJ9EE1VMpJVz\nw08CcJL1ij++vxbxydUJX1jYWK5mkc073x94pqePhV5w/rVe8M1QX+z6GqenvssOt7Vuf37vcCPo\nw/DY4foEhe7VjINUSIqN8e41UuMBX+RwCKc1wqKEeKU8/sn6Gxvj51FBJ9XSXvzCAIXCQrGwBL7z\nqVn8Y/HpPLPww0ZQ9MKCeYe5RgLhGXsWikWTVuf1uQiW7q7qp6J6C+YdRiLsLOQZLLKEfOHRn5qg\nVwur2t2EHGO9wVhp512tppkxebk3wshOxCtG3Zxanln4YTYuPp0LTpjqJ8fubTJlNe18xHqBjfHz\nGEN72Tu7jbAoId4u9avhBwColW1+bHZfaXPNUF2a29nYWroInuFG3Zxazj9hap5V+V2tZm4gMW+K\nvOu/lt7Nz3v0vf/2wWlAULOw+10TjIWdR3iSm1vxHmsDAO8WNIES4PwTppqNRR958vUm3nEDA17S\ng/bqM6KS4RK3PP4hsgUob2e3iYYqIV7vbG9fYmH36szuiZTbN/j9qf8iLjm1tpxvuoHmxrpZzD1w\nHFfd63SNC4bR7tQKJkkujNbzcfQHJ04fC8/A/OOnY/39dRId/R+xlkzbgDJJWtihFVS5PrGmQG6F\n6XGy52xpTdLIgXwydSOv6oF79RlR0k7zqwDlbAY0mkUJ8ZxmngOsghQx6VsJD49d6jjWtlHFZp3o\nHy/nm24wqJtTi+3a9zcFfo8v2jOcMFoXBQ6++uG8iqR7i227GwMrhEoIq58KCdY3NDL7hkeZ5oZ7\njmEXFZLiz/Yx/phgwUBjqtxzvOdrtR7Uq7+iN2Kk/aAGC7vsA1GMsCghntPMC62blMgyY3zx2O2/\nZY/kNdvJPi6snAom+mlv8RaA19yd4uv2AbyjY/OEBThO4zuXb+L82/62T99nexncVsQRFuy7sPCq\nygaLRXqa0V+yR/vHCsuzG/aMnqKj3rKdcORbM6f7r3vjs3Mn+cJiUsIu+0AUIyxKTN2cWkIh56az\nO9t4t3lHkSvgs+mFnNv1TW5Jf4ZH7GPzzpnop73HE7DP2jO5Jf0ZrkpfwlbGMp4dhMhyitXApaFc\nqZVn1jfvU/SKV6LcClmohAj1QyHBGx5cg61QRTvXh+/gANnK/q6/4i2t4Stdl/LvXZft9W7Y4NBT\ndPQddCUAABmESURBVNS56Wv4dvoCvps5n+9kzi/6GY+8vNE3Q8UHoBZZqTHCosRcU7+arqyzu6ig\nk7AWt4lnCNNKJUuydf5DX52IsNFEP+0TdXNqsQTShFmSdar0Nmk1IVEOkHf5dfQWFkTuZRLb/Wv2\nJZDAzjp/a7EiqNU/moVXW+yTob9xUfhRfhxZwn6usNjKOJbZJ/Kg/X5/fGEZC0Pf8aKjPIGxWWv4\nVfY0wOkdUww7nTM5Zzt3ccXSVf1i3hwsSiosRORUEVkrIutEZGEP52MistQ9/5yITHOPTxORpIis\ncn9+Xsp5lgqvwYrn4B4lnXvs4AbH7DSSW1/2J+cdn19t9V3X2f2N8F3+sfdY6/3Xja3JvX7APc1C\nXM2i0Nm5J3h+Cn+OblfFY6x1fMhyuuJt1bHdrjv3+AP2+jsNDj2ZpDq0uLCIBnwWFZJCgd8u31S2\nuRYlExYiEgKWAJ8AjgTOFZEjC4Z9AWhR1UOAHwI3B86tV9XZ7s+XSzXPUuG1tlTwb5gaaWViIPKm\nLxizU//iVaP1dtxeeOTHQy+wws10dnwY6vcT2Vv/RVCzwAr7eRZ7StBP4c1porSyya6hRUdzWuh5\nmrSKdCC4MSRiqsv2E4VFB6sTEZK9dNXzmC5v+2XpvWKhCmXbybKUmsVxwDpV3aCqXcA9wJkFY84E\nfuO+vh/4iEj5681e5nZhY6Ovhn/PKaGX+vw51YmIMTuVgBvrZrH+JsecsFZzO++b0/PJqlAjrVwX\n/m/Wxi8i7GqCe+O/UNvVLKwwSH4G957g+SmOl9dYG7+IY+V1xstO3tAp3O02dnrdztcg1t90mhEU\n/YhnkvrH4tMR6ZsZ6huRuznCcsqxRAP5Va3JdFlqF6UUFrXAW4H3m91jPY5R1QywAxjvnpsuIg0i\n8rSIfKCE8+w3vDpQly9dRTKdZRLb+ULo4SJlxnbPSG+POhCkiPKzzCe5NXM6K/UwtlNFDTv4XPgR\nIL9l5p76L+yMI2isUBi1wlh76eD2/BSnhBoA+FToL0yQHWzXMfwycxqPZt/Lf2c/vlefbdhzWjrS\nfkh7X5kkLVwSWuZrGlfd+xLTFz5UVjWjhmpS3tvAVFXdLiLvBepFZKaq5nWfF5GLgYsBpk4d3M5f\nnjaRTGfwKtwvif6E91pv0Kl9KyFxZmoR8T3M8DbsPZY42fA3Z871jzVpVV4o7aGymVfdkhl7mgip\nbp6FhEIgob3WLDymiFNVdry0MZ6dbKOKZsZwcfqqvHHV/VSyxLB73mE8X09/iXZN8NPoT4qO/5fw\nYwCssA9lpR5O1t04lKJkfakopWbRCAR14ynusR7HiEgYqAK2q2pKVbcDqOoLwHrg0IJrUdVbVXWu\nqs6tqSke91xKvDpQf419laXRRQB++et4L1Vlg7ykh/CcHuG/N+1RS0uhsxucJkkzAtpE0Me0p4mQ\nXp6Fo1mE9tpn4S3+3lyOkDeJSLZbBzxwBKAJhigt3t/j3uwpNOqEPbo2Jk5W98b4+Vwe/h1QPg2S\nSiksVgAzRGS6iESB+cCygjHLgAvd158GnlBVFZEa10GOiBwEzAA2lHCu+4xTWVaZItv8Xs/F6tz3\nhmmPWnoKnd3gaBZTrVxhwaCW0drRtUcmA89nYYXCe6VZeGZNzxw5EWcu3vwKe2uPrYjwg8/MHvI7\n1HLn+jNmEnErOHb0Urq8J6rZRZVb9PFyt2YckFeZeqhSMjOUqmZE5DLgESAE3K6qa0RkEbBSVZcB\nvwL+R0TWAc04AgXgZGCRiKQBG/iyqjZ3/5ahQ0iEcZpbWK4P38E0a+tef5Zpjzow3Fg3ixvrZnHi\n4idobE3SREH71YCw2NWV5fKlq1j5ZnOfnMca0CywwoREUdtGrOJ7tPqGRhbc/xLprBckod0yzbfh\naBbViQirrjM+i4HCey6vuvelPa5IuyT6E25Oz+92vBzyYUrqs1DVh4GHC45dG3jdCZzdw3W/A35X\nyrn1N1nVvC54F4Uf7WW0wx+y7+db6YuoCz3D6ECGp61qBMUA4/kjggX43rJrqGEHM2UjHwut5EeZ\nfwaEO5dvYu6B44r+jWxfs4iA5cTpZ7MZwlbxsMsbHlxDOqtMYAdfDi/jtszpVEiKNfaBzLTedOfq\nCDbT0GjgqZtTyxVLV7GB/fnf7LH8V+ZMNul+vBz/UtFr/yNyT7djhZGTQxGTwb2PXFO/moOvduTh\nnvaqSBNmJ6P57+w8fprNRRWbIoEDj/c7D1akXaPTmCgtPBT7BpeHH+AQyZmg+mJjzmkWIV9YZDJ9\nW9i9CKgLwo/xxfCf+FbkTgAa3O53kMsRMX21B4fJ1QkyhPly+gpe0YPYyah9+rw5ix4d0pFRRljs\nA7n+xzajSDKOtj5dl1VH5Uxr9/o9pkjg4OD9zoPC4m0dl2f6OdbKCYi+2JidaHDHDCWWo8R7Tu++\ncpT8A4D3WU4iV4M9wz/XhiPgUun+qWZr2DN6yuzeU2J0MQrnXvr/7Z1/fBT1mcffz+xuNhEJUaRU\nYjEq4g8syGlFW2q1UtTaKiqKrfWlJ9ard9rz2vMutoporYf1ruirrX2dnG1VKmJbpdjU31p/UioI\nbaSIoGAhQgNqgJBks7vz3B8zs5lJNr8kyS7keb9evNj97ncmk8nsPPN9fnyevu4D39eYsdgNgv7H\nN8XvZ1XpzFwDmq44N3UzP82eDXi9mSvKErnKUKvWLhzBOX9HDwTgb+4ItmoFQ6XNKBzmN7CBHvqY\ng5WFEw+tLHpmLIK9HyqbAS9dFuCNUOe7YFZTevcFCo3e076yOyw6eGHqRv43c1a3+5gee5FVpTOZ\nIOuAvu0D39cUa53FHkHgZwwKuI5x1kc+fzV7NJ+O/TUytoN9cu0asxJn9tnjzDgUEVupYEbqRrZR\nzkRnXeSzMSFjkVVl0Yq6Tv92i1bUsXzFRo4Hrpj/Op91P+BEBy64+yWuOvNTXW53x5NrUCBBhtGh\nLn7guZ4ua72Ot1zTfCoGpk2sjP4tZ3v/pUh0yIas1woyOIwKPVSe5rwOwDmxV/lzxnMxFmtjM1tZ\nfASClEYgV5EJMF6i2b3XZ67gQ903MpaihLR6NvrAYaVmKIqIoK5lqR7F21qZExkEL1YQXlkAXNuJ\nimhQoJlq9XSctrdAi+8p2rajiesfqc3raggyoAIXV5VsIS4uT2ePy83ZzhD+4E7kPdry+60Ir/hI\nE6dKotmQj2Y/w72ZMyNjwcNA+wZcxVjZbcailwQ3guALHeS+AxzmbI7M3aVludS657MTaNAhvuCb\nt7Ko35kqugtiMHPTl8eRiLU9DQYBZIAX3fFUyjYmyDq+FX84pyA7P4+KaFCgGRThZXByUvMObqdF\nWEEGVCVbqY4vyD18/D57QmhW9Gk14YgV4RUJ4esgRYL7slNp1hKW+QKVTVrKLqLJK2Mc7wHkYNkS\nGQ8qu4vp/mDGohcESrLN6Syz4vdzVWxxlxlQjZTmbioPZk/j2NQ8UpTkbhzprBbdBTGYmTaxkjum\nT8j5ntept+qr1wrWuAfhiPLb5Cy+GV/E6c5rue3aq4gGboS4X4SXJZZ7QAgk6vO5GoIMqNmJ+/lG\n/DFuStyHq8Izrrey+LN7aGR+ZUUZd1xg9TjFQPAQGZAmzhJ3HEelfpFLx26ktFNNqQp2Acr9if9i\nqn9tNaezRaVQazGLbgh8yHUNzQjkSnAujz8BwNLUkZH5O7UsFxT1uld4N4mGkDsq8GUKmrsg7Atf\nHAQ+6KrqGlwcPpuaS4sm2U+imW6fctbwuDsJ6Cj4OKqijLqGZpK+zlcLJaR8fbASyYB2nR59qO/u\nKpdmNroj2Mk+TE3d3kFa4pXqz+/eL2v0GcFqMmhQGLiawYs9gZdpt6uTiu8h0sxwdnByrJaTY7VU\ntXj9VRqa00y85SkamtKMqijjutOPKNi9wlYWXdDe5aR4KYw/TNydmxPOlgEv3bINoRXvJrE9lIMd\npDyO8Iv49lTJ4r2ZIH6xUUeylQrWatsXNKVxxjvR+NSh19fk4hdBSmVgLFIkSPnXQZJ0h/Toi+ct\noaq6BvC6KR4iW3Lik6+pN+8t/UTEhWFxiuKifSr1DvbJvQ5WkzsYwi5t+xvuDL0eQkskJvZ/iTso\n92VBPmxKoxTeNWXGogu8p4UMZzpLc08HC0q+z3mxl3NzZsSej2yzWYdH3gd+76aQ/v2rWc/H/Laf\nphn8LKN4aB+/UBymp2Zxc/oS5me/wDjZENF6ctWLX9ywqDaXUlkez5BRh6FlpTgJ74ly5D5E0qMv\nnreEV95uy445St7FEeU76ZnclTmP29Idez1bnKL4CFKpg4SWsGFv9F83aimNoZVF2NuQlAyX+VmV\nAFNiK7g45iXRHCT1/IO8BRRWdNDcUJ2waEUddQ3NnOis5qcld3FP5ixuy9Ok/YuxPwHwge7L/tLI\nJo2q336j9VrOiv0xMl7HCCan7mJLqA1msabLDVaCm3ngggRYpkeyLHskZzuvMDP+OBNlLWfHlvBw\n9nO8oV48YcHSjdw6zTMGf1xaRutmT7fpjZdb4Bn4zylVHB1yIwSG4oLYHziAHWR9F+VL7id5xO3Y\nJrWywK4IIz9BGv0pqR920Iu6IX05te6hLNexkQypFr/b3lYdxgjZnruXBFSKJxj5cvJagJxrqq6h\nucu07f7CjEUeFq2o47pfeR3tghVFuC9zwN2Zs/nnuCekW6/7sb808pYeFJmzheHcm20rzhEBVToY\nFZP4KD6C+EVQqR+wQr18+AdK5lAmrXzKeZMzW72OwGGNH8m0kJIS9gHiJd7fN9va9lAQuBMq2Mkd\niXsAr6p/rVvJVqKGwhFMUbaIqfTjVNvZt8NnH1CeK8RtDAW4h/gu7DfdTzAi5rmkf5OdzPm+58Lp\nQqTwul9796eBvB7MDZWH2YtXkXZdpjqvMcQX+DtK/sZXY89G5r3ktimPfqBDAfire3Cn+y1LxLh4\n0ugOEgEm8VHcBFLmgVNqo47kocwpJEmT0gRHyCb2o60vV1V1DeNmPUHjrl2k/VhFPBkYixZuWFTL\nIdU1XLtwJdAmI5Lyg6JzM+dHfr5Jjxc/+aQ/Eo5EXJnQ5p7apuVs1I8BXjwK4H0dynuhJIZJzuqI\nm/sg2ZpzR6Wzys2PDWymlK0s2rFoRR0NzWmOkr9xT8nc3Hi5NHFb4t7I3HDR1nczM/mSu4RlOpYr\nWr+dC2gGhN0Hxx+8P3c8uYb3GpoLnuFg9IxAyvyQ6hoUqM5cyezMpYyVTSxO3sjJzl/4rTs5N39X\na5bt2kgq5l0HCX9l8dLqjcx/L9q0aJKzmhZNMCE1D/AKNwMqyhKsmGXy48VO2G0Z/l6HxxSv/8XN\n6Ut4zp1Ikyb5x/iTPJudyMz44yRJ0xgKeh/mbOZ2Z17u/cvJfwXa3FEfNqUH1B1lxiJEOFd6mOzq\ndv5Obct4WK8H8qPseQC5vPiAyoqySJpjB4kAY4/h4hNH51xSLSSp1UPYpuWcGltJC0m+GX+Eb6ev\n4k0dTZJWmtwEn5nzHGeMauZGYP0WL0ZxYex5vhZ7hmvS13CC8yYr3MMjRiLAAtl7Dp19r4OxoGfK\nz7NtVdw/yFzE8KCjJq25YHhPGciWrOaGChHkSq9NXsL8xG3dzt9JGQ9kpnBX5txO55iLae8icEkF\nKA4vuBOYFnuVOxM/YZzzLlf5caxS0qRIUNfQzGN/9VqiJiUNKP8e/xXjnfUsKpnFeGc9SzVar5OM\nO9w5w1xPexOdqdS+TznL3LFcnf5m5AG0MzaUfpXb4t6KYyCzo2xlESJojZqQnkk+N5PkxszlkTHB\n6y9QDEU0Rv8QuKSC2ohfZk6jSrawTYdRRopTnJU4uCRpzdXZBG7JUlo5UjbyMWngwcznOVj+ThOl\n3JeJuprW3BrVEDL2fNpn2MVE/IQIYXrrbAAOk57VUHw1/jzfyXiNlgaqJasZC58v/PAPAF32pFjn\njuI/0lcyRFr8HOlo8KosETOJ8UFE8GV/XcdyfuvNAHzZeZWTY7WMl3cYLjtzWW+NlJFRh/1kJ6c6\nXmB7bub8DllPEJW6NvYu2ruqAtd3s9+T5G2t5OrWa9ikI8jg8LvkDT3a78XzlvDLr5/UL8ccMOjd\nUItW1FFVXcPa+l3EyVAdX9Dp3CRpXtexvOSOp1Y76vSYoRhcfGVSR5nwl91jcFWY7NQySrZR5xdp\nZomxhf2plPeZEltOrVuV11A4grktBxFBAWe4Iv937kms1DGsivQu6cjUkD7ZK29/wJjv/L5fq7v7\n1ViIyBkiskZE1olIdZ7PkyKy0P98qYhUhT673h9fIyKn9/Wx3bColqpQ+uJQmngzeRkXxl/odJtw\nCX+AI3DnjGN5pfrzZigGGe3jFwAfUs4bWsVl8ScZJk2Riv73dDhTnOUc56zl8eykDvtLxh1LkR2E\nTJtYycqbpna4lrSb2/M9JXP5Wuxpyvz0/oyrXLtwJVXVNXml83eXfjMWIhIDfgKcCRwNfEVEjm43\nbSbwoaqOAeYCt/vbHg1cBIwDzgDu9vfXJ4SLrOJkuCb2CLWlVxCX/B3HmjTJZt2fWenLIuOW/27c\nOu2TbJhzFhvmtBVePpqdzAGyg1aNUeO2GYWFmVMZJk0APOaeGNnPhjlnsebWM+1aGsTcOu2T3Dnj\n2Mgq48eZc3gqe1zn2yR+zurSy5kkqyPjgfRMX9KfMYsTgHWq+g6AiDwEnAOEW8edQ663FL8Gfiwi\n4o8/pKopYL2IrPP3t6QvDixoh/pcybc4SLZSEgpoBxWW6/XjLHOP4DB5jwezp/GEe0JkH5b/brQn\niGE8kP0C+9LMK+4xbPILrwB+436WWDpLsybZqCNz4xajMAKCmMbEW57iw6Y0/52ZAcB3dT6CMslZ\nzXYdwuTYKuq1Itc0aWHyezRpkoXZU7g5cynQJj3TV4hq5yXlu7VjkenAGap6hf/+EmCSql4dmvOG\nP2eT//5tYBKeAfmjqs73x+8FHlfVX7f7GVcCV/pvjwB6lENW8vExnZvqnqCqmR3177vNO9/drf3s\nnRwAbCv0QRSC2LCRo2NlQ0fk/dBLeslLZnv9erd5R/cN3PdOBu310gUHOGXlbrx8RBXSk2bvndO6\nZd3yHkzbpqpndDdpj86GUtV7gHsK8bNFZJmqHl+In13M2HnJj52X/Nh56YiILMs2bS+6c9KfAe46\nIJwucpA/lneOiMSBYcD7PdzWMAzDGCD601i8BhwuIoeISAlewHpxuzmLgUv919OB59Tziy0GLvKz\npQ4BDgf+hGEYhlEQ+s0NpaoZEbkaeBKv2eDPVHWViNwCLFPVxcC9wAN+APsDPIOCP+9hvGB4BvgX\nVe1ZWfXAURD31x6AnZf82HnJj52XjhTlOem3ALdhGIax9zDoK7gNwzCM7jFjYRiGYXSLGYteIiIX\niMgqEXFF5Ph2n/WrRMmegojMFpE6EVnp//tioY+pUHQneTNYEZENIlLrXx/LCn08hUJEfiYi9X7N\nWTC2v4g8LSJr/f87iogVADMWvecN4DzgxfBgf0uU7IHMVdVj/X+/L/TBFIIeSt4MZk71r4+iqykY\nQH6Bd78IUw08q6qHA8/67wuOGYteoqqrVTVfpXhOokRV1wOBRIkxeMlJ3qhqKxBI3hgGAKr6Il4m\naJhzgPv81/cB0wb0oDrBjEXfUQlsDL3f5I8NVq4Wkb/4y+yiWEYXALsmOkeBp0RkuS/bY7QxUlU3\n+6+3ACO7mjxQ7NFyH/2FiDwDfDzPR99V1d8O9PEUI12dI+CnwPfwbgjfA/4HuDzPXGPwMllV60Tk\nY8DTIvKm/5RthFBVFZGiqG8wY5EHVZ3yETYbVBIlPT1HIjIP+F0/H06xMqiuid6gqnX+//Ui8iie\ny86MhcffReRAVd0sIgcC9YU+IDA3VF9iEiU+/gUecC5eUsBgpCeSN4MOERkiIkOD18BUBu81ko+w\nDNKlQFF4M2xl0UtE5FzgR8AIoEZEVqrq6XuIRMlA8QMRORbPDbUB+KfCHk5h6EzypsCHVQyMBB71\n1bfjwIOq+kRhD6kwiMgC4BTgABHZBNwEzAEeFpGZwLvAhYU7wjZM7sMwDMPoFnNDGYZhGN1ixsIw\nDMPoFjMWhmEYRreYsTAMwzC6xYyFYRiG0S1mLAyjF4hINqSmu7IzJVkRuVNETu5iP5f6aZPhsQNE\nZKtfq/OQiBze18dvGB8VS501jF4gIo2qum83c4YDNap6YhdzyoF3gNGq2uSPfQM4QVUvF5HPAV9T\n1a/34eEbxkfGVhaG0fecD+SKzETkOBF5wRfNe9KXctgBvAB8ObTdRUCw2ngJmCIiVjhrFAVmLAyj\nd5S1c0PNyDPnM8ByABFJ4FX8T1fV44CfAd/35y3AMxCIyChgLPAcgKq6eDL3E/rzlzGMnmJPLYbR\nO5pV9dhu5hwIbPVfHwEcg6esCp7sRyA/XYPXJKscT9LhN+0kYuqBUfiGxzAKiRkLw+h7moFS/7UA\nq1T1pPaTVLVZRJ7AE1u8CPhWuyml/r4Mo+CYG8ow+p7VwBj/9RpghIicBJ5bSkTGheYuwDMSI4El\n7fYzFlNjNYoEMxaG0Tvaxyzm5JlTg6ckit9OdTpwu4j8GVgJfDo092k8V9NCDaUmishIPJfXln76\nPQyjV1jqrGH0AyLyMvAlVW34iNv/G7BDVe/t2yMzjI+GrSwMo3/4NjB6N7ZvAO7ro2MxjN3GVhaG\nYRhGt9jKwjAMw+gWMxaGYRhGt5ixMAzDMLrFjIVhGIbRLWYsDMMwjG75f/h5rQ6vr71VAAAAAElF\nTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f6a0e412128>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"model,posmin,posmax=mymodel(23)\n",
"kpm=pb.kpm(model,silent=True)\n",
"ldosz = kpm.calc_ldos(energy=np.linspace(-12, 12, 1000), broadening=0.05,position=[0, 0,posmax])\n",
"ldosmz = kpm.calc_ldos(energy=np.linspace(-12, 12, 1000), broadening=0.05,position=[0, 0,posmin])\n",
"ldosz.plot(marker='o',label='+z')\n",
"ldosmz.plot(label='-z')\n",
"ylim(0,0.4);\n",
"legend();"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEKCAYAAADjDHn2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXmYHFW5uN+ve9ask5UkHUJCNkgIZCAEBEFAIYAQRhQM\nwhUURO+Fq140Gq4IiCBgFPjpRVlcUBEim0MggYDsBAMJTEIIIXtIZrLv26zd5/fHqequ7q6q7pnp\nnunpPu/zzNPdp6q6K6lT9Z1vF6UUBoPBYDD4EejsEzAYDAZD7mOEhcFgMBhSYoSFwWAwGFJihIXB\nYDAYUmKEhcFgMBhSYoSFwWAwGFKSVWEhIueKyAoRWS0iM3z2+7KIKBGZ5Bi70TpuhYhMyeZ5GgwG\ng8Gfomx9sYgEgfuBs4FaYKGIzFZKfZywX0/ge8C7jrFxwDRgPDAE+JeIjFFKhbN1vgaDwWDwJpua\nxWRgtVJqrVKqCZgFXOSy38+Bu4EGx9hFwCylVKNSah2w2vo+g8FgMHQC2RQWIWCj43OtNRZFRI4H\nDldKzWntsdbx14rIIhFZNH78eAWYvwz/VdfUqWNumaeGz5ijTr3rVVVdU6cANf3JJeqUO19xPeZ/\n/rFYnf7L1+LGbn72I1V520ud/u8xfx33N3/1DjV8xhz13rpd0bHlm/ep4TPmqBc/2py0/7JNe61t\nWxTouTfp9pfV8Blz1KTbX47OPfOX8b+06DQHt4gEgHuAH7T1O5RSDymlJimlJpWXl2fu5AxRqipD\nTDvxcLqVBJk/4yyqKmMy22uWRZRCJH6sOBigOZz2vDTkAY0t2mpcUpT8mHGrMmSP2XOnqjLE7y4/\nAYD7vloZN/cMHU82hUUdcLjj81BrzKYncAzwuoisB04GZltO7lTHGjqQsFIEE57+Iu43POhx597V\nNXX8Y+FGDjS2cOpdr1JdYy5lvlNdU8cPn1wCwLV/XRS95vY08ls2OOdOMKA/hU0Nu04naw5uYCEw\nWkRGoB/004Cv2RuVUnuB/vZnEXkd+KFSapGI1AOPicg9aAf3aOC9LJ6rwYdIRBEIJAgLBOVxyytA\nrKdCdU0dNz6zlPpmvcqs21PPjc8sBTArxTwl8Zpv298YveZHD+4FpNIsYnMtKiwikSyesSEdsqZZ\nKKVagOuBecBy4Aml1DIRuU1EpqY4dhnwBPAx8CJwnYmE6jzCSkVvWptEM5MTpVR0dThz3oroQ8Om\nvjnMzHkrMnyWhlzB75rHNAtvTcE5tYqiwiLTZ2loLdnULFBKzQXmJozd7LHvGQmf7wDuyNrJGdIm\nHIGAi3TwNENB9I7ftKfedR+vcUPXx++a27PIVbNwESD2vMsFzaK5uZna2loaGhpS75yDlJWVMXTo\nUIqLi9t0fFaFhSE/iEQUwQQdVMTH7qxiN/mQinLqXB4eQypMQEK+4nfN/XwWiQ5ucJqhMnySbaC2\ntpaePXsyfPjwOFNZV0Apxc6dO6mtrWXEiBFt+g5T7sOQEjcHN4inZhFxmKGmTxlLeXEwbnt5cZDp\nU8Zm/DwNuYH/Ndczw63pmj3iKixywMHd0NBAv379upygAO0H6tevX7u0IiMsDClxdXD7+iziwx/v\nvHgCFeVa9T2sVyl3XjzBOLfzGPua9yzThoshFWXRa57Oc1bIXQd3VxQUNu09d2OGMqTEzcGt8YqG\nUnE3fFVlCIXif/6xhFnXfoYR/btn6UwNuUJVZYjV2w7wu9dX886Mz0fH/R5XbtqGcXDnDkazMKQk\nHHHJsyBFnkXCU6EkqM0SzeauLxiawhFKi4Ku29wd3BaOuRPIMc2ikDHCwpCSiHI3Q3lZkd3Gi4P6\n+KYWc9MXCk0tkaTsbdsU4hb5FHVwO8a6smZRXVPHqXe9yogZczKejPr6669z1VVXZez70sGYoQwp\ncdcsvA0KSiWH2hZbD42mrnjXG9pEo5uwsF79/NVO23ouhc62ho5ORt20aRPnn39+9PPSpUtZu3Yt\nRxxxRMZ+wwgLQ0rCEXefhZuN2R5PNEOVWrG3zUazKBgaW8KUBBM1C/3qPnWSB2MO7s6PhnLys+eW\n8fGmfZ7bazbsSVoY1TeH+dFTH/L4extcjxk3pBe3XDi+TeczZMgQFi9eDMD999/PG2+8kVFBAUZY\nGNLATVikMkMlCosF63YC8NWHFhCqKGf6lLEmIiqPqa6p44Wlm6lvjnDqXa9Gr7etkfrmWTjGYqGz\n2T3fTOOlQbdXsz7ppJNobGzkwIED7Nq1i4kTJwJw9913M2XKFObPn8/DDz/M22+/3a7fccMIC0NK\nwgqX2lB+Du74aKjqmjoefGNt9LOpD5XfxEww+sHovN4nHNEHaEOeRY6ZoVJpAKfe9aprYmKoopx/\nfPszbf7dd9/VPeJef/11HnnkER555JHots2bN3P11Vcze/ZsevTo0ebf8MI4uA0piUQUwQRNQUS8\nzVDE3/Az562gsSVZJTf1ofKTdOqB+WsWscnTVR3cHZ2M2tzczCWXXMLdd9/NmDFjsvIbRlgYUuLl\ns/AisUS5qQ9VWPjWhop6uL2Pdy40uqqD205MDFWUI2iNIpvJqO+88w6LFi3illtuYeLEiUycOJFN\nmzZl9DeMGcqQkrBS7oUEPfZ3ligHUx+q0PCvDeUXRZdfSXlVlaGsCYczzjiDM844I/r5c5/7XNYL\nHBrNwpCSiIeD20taJEZDTZ8ylrKEEEpTHyp/SccE45pnYb06Z5pJyssdjLAwpMS1nwXirVkkmKGq\nKkP85IKjo5+zrZIbOhfbBGMnYjqvt2+JcjdpgXZy50IhwULHmKEMKYlEks1QvoUEUUnmhqnHhfhp\n9TJ+esE4rv5s20okG7oOVZUhHnpzLUMqyvnDlZOi4+m1VY2fO8GAdEkzVL5hNAuDL9U1dXy8eR9v\nrNyeVLLAOykvuWBcqZ3BbZLyCoZwREV9DjbRPIs0mx9V19TR3BLhgTfWmP7tnYzRLAye2PHyzVZG\nlDNe3sdl4VpIsNjO4DZLxIKhORKhKOiukbq2VU1ofmTPP3tPk5/TuWRVsxCRc0VkhYisFpEZLtu/\nIyJLRWSxiLwtIuOs8eEiUm+NLxaRB7J5ngZ3UvVS9m6rmmyGCgaEgBhhUUi4axYav6qzpn97bpI1\nzUJEgsD9wNlALbBQRGYrpT527PaYUuoBa/+pwD3Auda2NUqpidk6P0Nq/OPlvZ0WERczFEBJUcCY\noQqIlrCiKKkfr37x9VlYc8vk5+QW2dQsJgOrlVJrlVJNwCzgIucOSilnJa7u+M8hQwfjlQdhj7ua\nEvQGVwd4cTBgqs4WEC2RiKfPwk21SBxKNf8MHUs2hUUI2Oj4XGuNxSEi14nIGuCXwHcdm0aISI2I\nvCEip2XxPA0e+MXL+9aGSuiUZ1MSDBgzVAGhNYvWRdE59zH92705//zzo5navXv35i9/+UvWf7PT\nHdxKqfuB+0Xka8BNwJXAZmCYUmqniJwAVIvI+ARNBBG5FrgWYNiwYR185vmP7US84YnFRBRx1WKX\nb9nXKgc3aDNUc4tRHguFloiiKOC+Hk2n6qw9/37wxBLCSuVWteIXZsCWpZn9zkET4Ly70tp17ty5\nALz//vt84xvfoKqqKrPn4kI2NYs64HDH56HWmBezgCoApVSjUmqn9f59YA2QVB1LKfWQUmqSUmrS\ngAEDMnbihhhVlSH69SjlssnDmD/jrOiNKnjXKHcrUQ7GDFVotITdzFAaXwe345CqyhChPuVUTRwS\nN/8MsGPHDv7jP/6Dxx57jN69e2f997KpWSwERovICLSQmAZ8zbmDiIxWSq2yPn4RWGWNDwB2KaXC\nInIkMBpYi6FTcI1q8TMneNSSKg6KERYFRHMk2cEdbavqm5EdP3eKApJ7/SzS1AAyyf3338/DDz8M\nwJw5c7jyyiu5+eabOeaYYzrk97MmLJRSLSJyPTAPCAJ/UkotE5HbgEVKqdnA9SLyBaAZ2I02QQGc\nDtwmIs1ABPiOUmpXts7V4E9LOOLeKc9DtfBqalZSFDSd8goI39BZl/29BEggIKY2FHDddddx3XXX\nATB9+nSOPfZYpk2b1mG/n1WfhVJqLjA3Yexmx/vveRz3NPB0Ns/NkD5eN723g5uk0NrqmjpWb9vP\n8s374jqnGfKTf35QSzii+L/XVvPPmrpYpzyftqpuZiiwNIsca6va2fzqV79i/Pjx0U55t912G1On\nTs3qb3a6g9uQ+7REFEGXqBbP21epOEOCXya4ERj5R3VNHTf+M+b8dV7vM8cOBDzmjktbVdA9LYyw\niMffjJcdTG0oQ0r8avy44dYpz2TiFg4z562godmjM2JUs/B+2Lll/xth0fkYYWHwRSmlNQuXEMh0\nCwmaTNzCIq1OeS54+cCCuejgLkCMsDD4Yi/o3KKhvDvlxdeGMpm4hYXf9U6nn0WiPAnmkIO7M8w/\nmaK9526EhcGXFusmTW5+5O3gjkTAubvJxC0s/Doj+rdV1a+Ju+SKGaqsrIydO3d2SYGhlGLnzp2U\nlZW1+TuMg9vgS4ul/ydqFl72hOqaOlZu3c/HLlFPNz7zIfXNkdzKxDVknKrKEHvqm7h1tq4Z6rze\ns97bAMAdc5fzyDvrXedBUvMjkeiipTMZOnQotbW1bN++vbNPpU2UlZUxdOjQNh9vhIXBlxZrReeW\nZ5GIHfVkH5MY9fTqJ9tYWreX1354RtbO15AbfP6ow7h19sfM/MqxXDJJF3Korqnj1ueWRfdJnB9e\n6/WioNCYA2ViiouLGTGicLs8GjOUwRdb/S9OzMS1Xp0qeaqop6KgmEKCBUKLy7zxjZIiNpcSlVYT\nOpsbGGFh8MXTZ+GSXJUq6qk4EIiatQz5TUs4ed6kmh9eM0NHQ5l509kYYWHwxV7RefZSdoylinoK\nBiW64jTkNzHNIv2oOH8Hd+bP0dA6jLAw+DJ36WYAZjyzlFPvepXqGl042M2/nSrqqTiQG45KQ/ax\nNUhnfs70KWMpK3aPknLidHBX19Tx1qrt0TIx9vwzdDzGwW3wpLqmjl++GMuydjokbbSdWd/cdlTL\nD59cQkskuf9AUdCYoQoFe1HgbH5UVRmiJRzhh099COASFRc/N+yACdvPYcrEdC5GszB4MnPeChpb\n3B2SXtVDqypDDOvXjS8eOzip/0CR0SwKhhYP82W0odHZY5LmR6IZypSJyS2MsDB44uWQrNtTz8Nv\n6fYin/vla8mmAZWchQt6lWk0i8Iglp/j3s8i0XVVXVPHT6zig1f+6T2qa+pMmZgcwwgLgydeDkkB\n9jW0ALBpbwM3PrM0TmC4lSgH/eBoiagumQFraB1uZiiIZfY760DZ5qZdh5oB2La/kRufWUpFt2LX\n7zZlYjoHIywMnkyfMpaSouT8isRHfaJpQCWUKLexTRImZj7/8TJDuWkWXuYmpTBlYnIIIywMnlRV\nhvjmqcMBLSRCFeWesfBO04AivjaUjd1i04TP5j9eZqgoDu3Sy6y0t76ZOy+eQI9SHYcTqijnzosn\nGOd2J2GEhcGXE4f3BaD6ulOZP+MsQmlUkI0o5WGG0mMmizv/CXuYoUAvJNLNz6mqDPGNU4cjAm//\n+EwjKDoRIywMvtjRUKVWfHw6FWQT+1nY2A8OY4bKf5q9ClCiTVERh2aRak6VFgVQKvadhs4hq8JC\nRM4VkRUislpEZrhs/46ILBWRxSLytoiMc2y70TpuhYhMyeZ5GrxpsoVFkb6ZqypD3HnxBHqVadPA\n4N5lSaYBR+pFlOqaOu771yoAzr3vLZNclcdU19Rx87MfAXDFH99NutYBiS8TY8+pPpZDe2DP0rg5\nZc+9JqORdipZS8oTkSBwP3A2UAssFJHZSqmPHbs9ppR6wNp/KnAPcK4lNKYB44EhwL9EZIxSKt4L\nZsg6jS36v9zp6K6qDLF9fyN3zF3OS/9zOj3LkqNWErNwb3xmadSJuWVfg0muylMSr/XWfY1J11qQ\npNBZXXVW8T//WMI/vv0ZRvTvHt1mz73G5nDUf2HoeLKpWUwGViul1iqlmoBZwEXOHZRS+xwfuxMz\nZV4EzFJKNSql1gGrre8zdCDVNXX8Yu4nAHz59+/ErRCjhQRdjlNKmR7cBUpa11rcW6ja2obTclVd\nU8e9L68E4PzfGI20M8mmmA4BGx2fa4GTEncSkeuAG4AS4CzHsQsSjk1agorItcC1AMOGDcvISRs0\nSdrA3nhtwHZgKxfLQGI0lEmuKhzSudYBt/hrYuG0tlaajpZi6Dg63cGtlLpfKTUS+DFwUyuPfUgp\nNUkpNWnAgAHZOcECJdUK0S25yiaiVJwZyvTgLhzSudbaDOWmWcT3szAaaW6RTWFRBxzu+DzUGvNi\nFlDVxmMNGSbVCtEWBW6BTUrFV6U1PbgLh3SutSQ4uG0Sa0MZjTS3yKawWAiMFpERIlKCdljPdu4g\nIqMdH78IrLLezwamiUipiIwARgPvZfFcDQmkWiEGLNXCrXSHLvcR+2xHu/TtXgLAgIRoF0P+kBgt\nN8QlWi4gyQ5uiGmpAWvyGI00t8iasFBKtQDXA/OA5cATSqllInKbFfkEcL2ILBORxWi/xZXWscuA\nJ4CPgReB60wkVMeSaoWYSrNIjJ2tqgzx60uPA+DB/zjBCIo8pqoyxDWnHQnAmz9KTqTTLgs386W1\n3Zo6RiPNLbIah6aUmgvMTRi72fH+ez7H3gHckb2zM/hh3+AznvmQhuZIUu+BqIPbPR7KtTlSsVX6\nwSTl5T9ubVVtUpmhbM3Cnmt3vrCcrfsaqehWzK0XjjcLjU6i0x3chtzGq0CsWw9u5zHutaFMuY9C\noTmiKA6Ka9kXEXE1X9pOb+cRVZUhvnuWtlbvOdTMzHkrTPhsJ2GEhcEVO2zRLvdhdymLtlW1e3C7\nhkDGR0PZ2KUfTE+L/KclHPEsIijiGjkbG0vIs/j5nFgeb+I8NHQcRlgYXGlP6Gyig9vGrjprzFD5\nT3NYuRYRBG1mctVYVbyDG/Q8tNuq2pjw2c7BCAuDKylDZ6372TN01uVYU3W2cGiJRCgOemgW4Jpn\nEUvKi2HCZ3MHIywMrqQKW4w6uD2Sq1xLlFsrTdPPIv9pCSvXirNg+SxcxpWLZmHCZ3MHIywMruiw\nxfjp4RY66+rg9vhO24ZthEX+0xxW3pqFuC8yEkNnwYTP5hJGWBhcqaoMcdHEIdHPQRG+fEIoGrYY\nEG8HNyp+dWhTbGsWxgyV97REIp4+C8F/keHUSqsqQ3z5hFiobOI8NHQcRlgYXKmuqaO6ZlP0c1gp\nnn6/LhYNFfVZeNSGcnlOBE00VMHgZ4bycnAn1oYCPQ+ffj8W+ZQ4Dw0dhxEWBldmzltBQ4t3FEpU\ns3A51qX3EUDULGHMUPlPc9jHwS3ui4zEpDwwxQRzCSMsDK6kHw3lftO7hs7amkWkwMxQK1+Cbcs7\n+yw6lJZIitBZl3G3pDwTDZU7GGFhcCX9aKjkfRQe0VCWg7ugeik37IPHLoE/n9fZZ9KhNPsk5YHH\nIsN6NdFQuYkRFgZXpk8ZS0nCytA9GspDs3D5zhc/2gzAz5//mFPverUw7M7hZv1av7tzz6MDqa6p\n4911u1i8cY/rdRbP5kfJPgsTDZU7GGFh8MTpoOzTrTiu1HQqn0WitKiuqeOW55ZFPxdM2YYCK5Zs\nl4lp8igTA95mqMR+FhCLhrLHAoKJhuokjLAwJGHf8IccZRYSSy74+yxUUuhswZZtcOs7m8ek45D2\ndnDbPovY3LGjoezdIwoTDdVJGGFhSCKdGz5aGyrNch8F66iMFJZmkV4Pbq/QWXt7bMxEQ+UORlgY\nkkjvwa7vaC9HZaJ/u2AdlQWmWaTXgztFbSjH5CnYRUYOYoSFIYl0bnh/zSK5RHnBOioLTFikdZ09\nS5TbtaFiYwW7yMhBjLAwJDF9ythoaQ6b4qDE3fD+obPJmoXdm9keDlWUF0Yf7gJzcNvX2c7Wd7vO\nAY+GFm6ahZvwEeDMowZk/NwN/mRVWIjIuSKyQkRWi8gMl+03iMjHIvKhiLwiIkc4toVFZLH1Nzub\n52lwIfFmTvjs28/CI3S2qjJERbdi/uPkI5g/46z8FxTg3Wowj6mqDDGoVxkXHx9yvc5eZihcysRE\no6Gcu2Gc3J1B1oSFiASB+4HzgHHAZSIyLmG3GmCSUupY4Cngl45t9Uqpidbf1GydpyGZmfNW0JxQ\nkqM5opIiWiC5n0Wsvo979m5RMFBYGdwF5uC2aYlEoj3XE/FycEc8FhmvfbI9aUlinNwdTzY1i8nA\naqXUWqVUEzALuMi5g1LqNaXUIevjAmBoFs/HkCbpOBVjbVXjb2O3WHknxQEprEKCBeazsGnx6ZTn\nGTrrkflvnNy5QTaFRQjY6Phca415cTXwguNzmYgsEpEFIlLldoCIXGvts2j79u3tP2MDkGZEi5dm\nYW93XSNCMCiFVUiwQIWFXyFBcHdwR1S8c9vGOLlzg5xwcIvIFcAkYKZj+Ail1CTga8B9IjIy8Til\n1ENKqUlKqUkDBhiHV6bwch46x2MrwETNIrlkg5PiQKCw2qoWmIPbpiXSlhLl7ouMdOajIftkU1jU\nAYc7Pg+1xuIQkS8APwGmKqUa7XGlVJ31uhZ4HajM4rkaHLz2ibuW5hz3Cp2NaRbuFAWFsNEs8h5t\nhmpdpzxthkreP535aMg+2RQWC4HRIjJCREqAaUBcVJOIVAIPogXFNsd4HxEptd73B04FPs7iuRoc\ntMZnkezgtrZ7SItgIFBYVWcL1MHdHIkkhV/beETOepa2Nz6L3CBrwkIp1QJcD8wDlgNPKKWWicht\nImJHN80EegBPJoTIHg0sEpElwGvAXUopIyw6iNYl5SWYofCPhioOSmFFQxVg6Gw4olAKzxLl2gyV\nXk0xMD6LXCGrPgul1Fyl1Bil1Eil1B3W2M1KqdnW+y8opQ5LDJFVSr2jlJqglDrOev1jNs/TEI8u\nTx4/NdyycKF1mkV1TR3LN+/j9RXbC6dEeQGaoZ75oBaAe/+10r1EOcnzBrxDZws2+z/HyAkHtyG3\nqKoM8fXP6PxIwScLl+SkvKiwSLjt7Uq2tgnKlCjPT6pr6vjpsx9FP7tdZ/EpUe6mWdhZ4b3KigAY\n0rusMLL/c4yizj4BQ24ycVgFAC9+/3TGDuqZtD3W/Ch+PGaGih/3qx6a1zd9gWkWfqXo7evs5eCO\neKX+owXGzoNN/Pz5j3nhe6fTu1txxs/d4I/RLAxJVNfUcVO1Xh1+/U/vuq7+AwH32lAxzSKegnVS\nOh3cBeC/SC84wvu/wk2zAD0nf/vKKgDOue+N/NdIcxAjLAxx2OaiPYd0O9Ct+xpdzUX2LZ2YiRsN\nnU2zRHnv8jxfITo1izzXMqpr6jwf9vHBEeJaUyziUhvK/t4bn1nKnnr/OWnILkZYGOJIt9lMtOps\nwvG2eSHxoTF9yliKXZK0Dja15PdN7xQQeRxGaz/Qwy4qQ6IzWgTcAuK8fBamAVJuYISFIY50zUVe\nbVW98u2qKkP0KEt2kTWHVX7f9E4Hdx47u90e6KDnSaIzWvw0C5fvLlgTZo5hhIUhjnRj2qMrQI9S\n5m55Frst01Yidfl80zuFaR5rFl4PbjffhFforO6DYvIschUjLAxxpBvT7u2zUHHbnXiUCvIsDZIX\nRApDs/B7cCebMHFN4VYePguTZ5EbGGFhiMOOae9eom9Or4520TyLNJPyqmvqPE1UytqelxSIz8Kv\nqF+i5ujl4FYeVWftOTmwZykAfboVmzyLTsAIC0MSVZUhzp8wmMG9yzw72nn5LLwKCabyS+St36JA\noqH8ivolLhx0P4vk/SIuvdttqipDPPHtzwDw0wvGGUHRCRhhYXCloSVCWYLq78TTZWFHQyUsEVM5\nI/PWWek0PUVaOu88sozf9UvUPr1rQ3mbKoHofGxsyV+hm8v4CgsROVFEBjk+f11EnhWR34hI3+yf\nnqGzaGgOU1rkPT28OuXZK8bEez6VMzJvnZUFolm09vp51obyKlcM0fnY6BJ1Zcg+qTSLB4EmABE5\nHbgL+CuwF3gou6dm6EwamsO+moVdUDRxgfjCR5sA+Omzy+KKyE2fMtbTkS3W9rwkzsGdv8LC7/pV\nOBIvq2vqeG/dLhZv3JNUZNCrn4VNabElLIxm0SmkEhZBpdQu6/1XgYeUUk8rpX4KjMruqRk6i+qa\nOt5d635D27j1s6iuqeMXcz6JfnYWkauqDHH5ycOSBIYAl588LH9t0E5pmsfCoqoyxBUnD0saLw4I\nt04dD8QS9+yHfWKRQa9+FjYvLt0CwJ0vfFI4VYtziJTCQkTsTKrPA686tpkihHmIfUM3hd1vaJuY\nzyL2MJw5bwUNLe5F5ABur5rAvV+dSMgyWZQXB7n3qxO5vWpCtv45nU+BmKFAX9/h/bpRWhSIViue\neclx0YVAqkxs5ePgrq6p4yfV/tVsDdkl1QP/ceANEdkB1ANvAYjIKLQpypBnpFsd1q2tajqZtlWV\nIaoqQ5xz7xsc2b9H/moUNqowzFA23UqKOG10f/5w5YlJ21LND4W3g7tgqxbnEL6ahdWw6AfAI8Bn\nVcybGQD+O7unZugM0i+tYJuhYtIi3Uzb6po61u04yIvLtuS/OSFOs8jvqrPVNXWs2Lqffy3f5npd\nU80PPwe3KfnR+aQMnVVKLQD2AN8QketF5Eyl1Eql1AfZPz1DR5PuA//VT7YC8L1Zi6MPhulTxiZF\nUCVm2hZcE6QCScqLFhKMeF9Xv0zs6po6/vXxFtbtONgmQWPIPqlCZ0Mi8i5wK3Ck9XeriLwnIil1\nPxE5V0RWiMhqEZnhsv0GEflYRD4UkVdE5AjHtitFZJX1d2Wr/2WGVlNdU8fBxuRcALcH/q9fWhn9\nbD8YAH54Tmw/t+zvgqsgWiDRUOlcVzsTu9yKarLnB8CNzyylvtnbT+YmaAT/zHFDZkmlWfwf8Hul\n1OeUUjdYf5+zxn/nd6CIBIH7gfOAccBlIjIuYbcaYJJS6ljgKeCX1rF9gVuAk4DJwC0i0qd1/zRD\na6iuqWP6U0uiPQNs3EorzJy3Iil80X4wnHnUQAB+c1mla/Z3wZkTCsTBne51raoMcfa4QQzv1y06\nP9IVNF/hX4ZOAAAgAElEQVQ+IX4uKeDvCzZwU/XSzPwjDL6kEhbjlFKPJA4qpf4KHJXi2MnAaqXU\nWqVUEzALuCjhe15TSh2yPi4AhlrvpwAvK6V2KaV2Ay8D56b4PUM7+Nlzy6KmoURa88C3zRDBNJrg\npDPe5SkQYdGa6xoMSFzfi3QFjVtJEVtg5K0ZM4dIJSxct4tIAPDO2NKEgI2Oz7XWmBdXAy+05lgR\nuVZEFonIou3bvWvTGFLjVT7cbdzvwRAVFh4zq+AqiBaIsJg+ZSxlxf7+KpuASFzzo3QFjWcZdPK4\ntlgOkUpYPC8iD4tId3vAev8AMDdTJyEiVwCTgJmtOU4p9ZBSapJSatKAAcZ2mS0SV21+D4aIR6c8\nG9tu3aPUv6pt3lAgwqKqMsSM82LGBr/rGgwQXVRA+guIim7eLXjz1oyZQ6QSFj9C51N8KiLvi8gH\nwHpgH/DDFMfWAYc7Pg+1xuIQkS8APwGmKqUaW3OsIXNU+PTCvnX2srjPVZUhbr1wfPSz88EQ0yy8\nU3GrKkN887NHAvDWj87MX0EBBdPPAuBzY7S/6t6vHudZrRggGAjEmaHsBURJMN7xnXi8X02ovDVj\n5hCp8iyalVI/RD+4rwKuBI5QSv3Q8kP4sRAYLSIjRKQEmAbMdu4gIpXo+lNTlVLbHJvmAeeISB/L\nsX2ONWbIEnZJBjf21DcnaRdfOl7fyNOnjI17MIQ9qs4mUl4oFUQLKM/iUJOOpEvUEhJJ1CxAC4zx\noV6cNrq/q6CprqnjULP3XMlbM2YOkTLPQkT6AdcA37H+vmmN+aKUagGuRz/klwNPKKWWichtIjLV\n2m0m0AN4UkQWi8hs69hdwM/RAmchcJujRpUhC1RVhujjo+Yn2oSLrUqCD7yxhhEz5kRj4yPWQ6Ao\nhbDoZjVXsh8weUsBZXDXN+l/a3mJf2GIoEicsKiuqePUu16lZsMeFq7b5eqs9vNJVJQX57d2miP4\nXlURORpdD2oeOsxVgBOB/xWRs5RSn/gdr5SaS4JvQyl1s+P9F3yO/RPwp1T/AEPmOH/CIP7+7kbX\nbYk24dlLdHXZIY1rWcXQaGz81Z8dDnhHQ9nYq8/EkMm8o0B8FhC7lvZCwItAQKKLCjuZr6x5N4No\nZktLv2jOjlMA+PkkLjhucHtP3ZAGqTSLnwPfU0pdpZT6f0qp+5RSV6JLfdyR/dMzdCRzrKqebiTa\nhGfOW8EY2ci80hlcF6wG9MNi1kItbPzMUNU1ddz5wnIALv7dO/kd9lggwqK6po7vzVoMwH/9/QPf\naxqUWOisnWPxWukPWFCmKwi5JWn6+SSeX7K5vadvSINUwmKCUuqJxEGl1NPAMdk5JUNnUF1Txx6P\n8FlItgnX7alnpGjt4rjAmuj4jgPaleXl4LZXknZI7rb9jfld7iOS/8LCvqa7Duprvz3FNQ0GY2Yo\nW2OokIMAFNESN27jl6nt5lMzZJ5UwuJgG7cZuhj2Su6mor9xceDNuG2CwyQQibDouQfpTgNDZAcA\n9ZQCiimBhQztpk0RXqGzBVfuowA0i9Ze06jPYvsKTin7NG7bUNnOSKnjzLJVceOJCXm/KHqY8wLv\nxp2DIbukKlE+UERucBkXwCQ25BF6Jae4pkjnRT7TcDohtlNHf5Szx8DSJ5j0/o+4MngpPaQBgDKa\nOSWwjAdL7uXdAdP46qdTPTWLwiv34ezBnZ/+mdZe02gG9/2T+Tswlkei2/qxjz+XzKQXh6DpW1DS\nPfpd5TRQThO76MXXil7ja7zG8IbHfH/LkDlSaRYPAz1d/noAf8juqRk6kiEV5fRhf/TzeFnH/LLv\n8bXgq9FmRQBsWhx9W8EBAHpQzyTRhQUPbNIrvLdWuWfUm3If+Udrr2lAJO7/4lhZG33fRw7QS6wK\nQNtj8TNDKsp5tuSnfFD2HYpJjqDL2/mTQ6TKs/iZ1x+m+VFeMX3KWIYGYtHJI0Q7uz8XXMr0c8bA\n9pU6TyCs7dI9pZ5elp25XBoYLDsBqBAtQO55aaWrHbmwy33kZ55Fa6/payu20tthxR4Z2BR930di\nCxYO7YaGfbBvM9OnjGVMQM+n7sRrEXk9f3KIlHkWPriZpwxdkUiEqgkDuGLSwOhQCdoBfdzgblQ1\nz4H7T4SV86BeC5Q+7I9qFt1oJGT5Lyo4wLGyht8W/4bbn10S/zurXqbq6B7cefEEBvUq0/uXJ1e1\nzSsKoES5nYHdPbGES2g/bPkofudty/mvrT9jkOyODo2S2KLCnlMAi5avhnn/C/ccRdXImFlziLUw\n0b9Vxt0Xjc3f+ZNDtEdY+AfSG7oO1f8Jtw/gmIEl0aF7vjgEgEE9i2Kmp52r4JAWFn1lfzSCpZs0\n0k/2AVqzuLf4d1wQXMCwRofTceN78PevwL9upaoyxEs3nA7A9WeNyu8bvQDMUKAFxoXHDmFgz9JY\nBvbvToIHToWwI8ru+Rs4N7iQc4PvRYeOkFjxhr4OzeKtD1dCzd/0h+2xuTT3iiHR929c3pupz0+E\nNa9m4V9lcNIeYZGfOnUh8uEsACINThOAtXoLN0PEutn3b4lqFj2oj5oMutGgHZLolWEQ/VA8NrAW\n3pwJL94ItQv1d2zVdaZ2PXo1txY9wu1zlud3a9UCERbVNXW8sHg9f2z8AXfcfhPPvr8+tnHTYvjb\nxbD+bWjW8+R4iUU7DXZoCs73gYbdELQWMPsc82NPLHE0vOVj/ebfvu11DBkgVQb3ftyFggDGo5Rn\nBA5ujX04aN20DXshYNmj923SdmSguzTQF1tYNBIU/SAMioralI8q3gqv/lUfO/Hy6PfNXbCU82uf\n5aoiuLXlqrhOe3mnZRSAsLDzLAa3bGVC6XomtPyWc6pHcZHtxlj4B1jzCpRXwEEd+HBkIJZIZ4dg\nAwx1vi+tJ2rA2BnL5eFgTBNpaW6iFKApZr4yZIdUDu6eSqleLn89lVKpwm4NXYwix01o39Q07I1p\nGQ7NYoDso1yaaFDFlEkzveUQtap/dBvAKd1qY99naRQc2MrfX421by9FO8zzNtciTljkZ+isnWcx\nwBHz0i/smEvbrNX/zjVRDcEWCnWqH32toIh9qpyhEouiO+kwFdNqdzmExYHYd6uDpo9NR9EeM5Qh\nH3BkGBfXuwmLPVE/BbvXRc0I9mqwVsXSbTZEYg5ygGGHHKXNt+nyHtTvpuhAzKQwWHZyrKzh5MDH\n+Rkrn+8ObqU4df8L9OYAA2RPdDjk0BCiIbCbF5PIhshh0febVP84x/fQSF3s/2xnLLyWAzENOHJg\nR/Q8DNnFCItCpyX2gK7buC76/tBuy0zQsDcmOPZbY8XdYsdY2gTABhUTFu9GErruhhujbyd3i5kg\n+rKf2aU/ZVbJ7fkZK5/vZqida/hl8UP8tvi3ccLCGeFkh1vbLI0MB6BRFbOFPtFx51yiqDzOqe3U\nLPbuiIXa/vtDOxfDCItsY4RFodMUi3fvGdkXfR84ZK/YItDSAN0dWkPfkdG38cIitkp8Onwa70dG\ns6tiQuy4Ih0uO3VQbPU4SGK5HYcOHcg/R3dcifL8e6D98YX5AEwKrGSAxMxQYwKbaAmWxxYWZRU0\nFffinuavRLXRvXRnj+oRPWaTcnQ+6DcqJmSKyqMaLUDL3ljBy24tWkDt2ecIzjBkBSMsCh2HY7C3\nI8a9TBKKCg5y1I0MHR99W+uhWexUvbik6RYuqr8FulkPgcN0g6Ve+2IrxlESWyX2bd7C9KeW5JfA\nyGPN4qbqpXy0XJsXWwjQ3+GzGCW1NJb2g7IKPTD0RM4O/IHfhC+OCoi9qjt7iXZsZrNTWIQqY+8H\nHxf3u3aYtn6vhcS+fbsxZBcjLAqdptiKzc6bsIkoRyrNoGNj70MnRN86fRbO97tVTyIE2Li3CXpY\nGsfAcQD03reSRqUbLY0KxATDQNlDc1jll6PbqU3kmbB4/N2N0VDXRorjzFDDZBvr67vpCCiAbv3Y\nsFeX6diDFha76BmnWWynd+zLh58Wez90kuc59LUER2kkD/1dOYYRFoWOwwzVR+LDD7c5GyKOPAu6\n9YdJV0PPWLMZ52pwN7Eb334gBEWgxBrvNxIkaB3Xl/2qPE6zGIheHeaVozuPHdxhpRhsmRHLaWKA\n7GVVxNGwqKVnTLPo1i/qk9ptCYhDqjT6HmC7cgiLcRdBnxEw5ry4+bbPoYmALjwI0NMqamnIHkZY\nFCrNDbB/KzR7V5ov7jM09mHQBJi+Gi64J2ZWApoCMaf0fhVzfO9SPQGrJ7f9kOw5BLprs9UuerJT\n9Yr2xACiNu+8cnSrCASsKPM8qzobFIlqFj1El6xfqWLCor6kj0Oz6BOt37QbPTeKCLPXscDY6DBj\nUlQK362Br82CHrHx4jK9f9jSeovFauVKo47s2+Pe6dHQfrIqLETkXBFZISKrRWSGy/bTReQDEWkR\nka8kbAtbfbmjvbkNGeT1X8Cvx8RVkQU4qEqj7/uNcqj/5X3A7lFRMQyAnX2P54Aqi+7itD/vsR4I\nFeXFMUdl76FaO0ELk130pNThGxlomTH8Gt10OVQYAlZv8zzTLE4+sg9DHAEKfeUA69Wg6IN87JFH\n6uAIgIojqKoM0b0kyL8j2ne1Tg2OM0NFzZgleu5E59uAWJHA8u56225xaCGAoFj5+HS47xjYsCBj\n/0ZDjKwJCxEJAvcD5wHjgMtEZFzCbhuAq4DHXL6iXik10fqbmq3zLFjefRCAFe/Nixu2V32AfrCf\n8l04/1exGxegxwD4r3f52u5r2RLpFR0OE+TUhv/HyQ2/jY6JAP3HWMcdpo8Fdqle7FSxY7eqCgbK\nbi4I/JtpS74JB/Ik2cqpWeSLsFj7OvzhbJq3r2Ww7GSX44G/TfWJFo0bO3JkdGFhz4HiYIBaNYDz\nGu/kzpbL2ONYYDRSwhcbf8Hcs16I/71Bx8KES2Ha49HoKumRvKBoWvGKfvPp/Iz8Mw3xZDMLezKw\nWim1FkBEZgEXAR/bOyil1lvb8uQu6iKEW6IrvtK96yCgndkBURxQ5bESkSXd4cwb3b9j4FGsqF+T\nNFyX0BNrz6FmuPA+mPAV6D8KemnT1h7pFXWgH1KlbFQDGcge7i5+iO6qEZY8Dqd+NzP/3s5ERSCY\nZ8Lig79C7XtMC5fQJ3iAt8Pj+WxQJ2BuV70JiOXUrzgcKq+A0efAkIkA7K3XmuRydQRAnGYBsEwN\n547Xt3P+yY5BEfjyw/r9O3ohsvJAGZ9JOK1edjTfpuTkP0P7yaYZKgQ4DYi11li6lInIIhFZICJV\nbjuIyLXWPou2b8+TlWhHsH159O3wgM6G3WlpFI0Ux/Yr6UY6rIyEeDM8wXXbkIpyKO0JY8+L+87P\nHnMku9CaxW56sFEN4JTgx3QXK3lvwwLYtxm2LE37n5WTRML5oVlEwrDuTV1Y0jLzXBx8G4AlKpZ3\ns11VxI6pGAalPeCoL0aHEv1RtsN6ZvOl0bE6vwAHa/5sDcc04B2WhjosYD0DNtXkZU5LZ5PLDu4j\nlFKTgK8B94nIyMQdlFIPKaUmKaUmDRiQR3bubLNaq+sNKiYYdlk3XJywKI6PPEnEVkDOaZrJ15uT\nNRDXpjTHfBkGHcu4877DaROPBiCsAnwQGR3dZcthp+sqtf/8NjzwWdi9Ps1/WA6ilMNn0YUfYAv/\nAH+5EF79ua7vdPSF0U2vhGN5N/uL+vLu5N/AkEqdWJdAYqOkCAGGNzzG/eHYejDo0b8diJqh9qtu\n0fDrOAEFsHcj7FjZqn+eITXZNEPVAYc7Pg+1xtJCKVVnva4VkdeBSiDZ7mFoPds/gV4hVu/vzjFq\nJQ2qmANWEeFGlb5m4ffoC1WUM32KS1OaYSfDd94CoL5ch0QK8Ez4NMbLet6ikmsP68GgrW/COqtW\n1eYPYfZ34divQuXlrfmXdj5d3WfxzLdhwJhYlNH8/wfAK4OuYdWyInaFy6hRMaHwmcpjOOn8E4Er\nXb/Ong8z563w1CDCfkK1WM/TfXSjgWJKaWa76s3R1uZXwpV8Plij58wA0z0vk2RTs1gIjBaRESJS\nAkwD0opqEpE+IlJqve8PnIrD12FoI02HtBlh3yboFaLfYfrG3U859Ur3DYjTLHoNdfuWKCGPENdQ\nRXmsAY4Pt32oI1reiYznIOXMaLmWOS0nct8n8ZEubHwX1r0Bz/6X7/flJCrs8Fl0sdDZXet0r5NX\nboM9G2Lj3fpx8zth7mq6hIfCF6Icj5GXVqYuu1FVGWL+jLN8548nluDdpvqgLN12OzHNQmuoAjtX\nQ+OBvAtX7kyyJiyUUi3A9cA8YDnwhFJqmYjcJiJTAUTkRBGpBS4BHhQRu0zp0cAiEVkCvAbcpZQy\nwqI9rJ8PvxoND51hCYvBNJXoIm77VTfqdVeAOGExt7bU7ZuitLef9rK9pZzTeDc3t1wVN/72vvjq\ntax7M63vy0lUpOuGzu6L5cDEmQJDJ7Bpb3wS3EkN/8fkhvtblVDZpvkT1P+XAgSsJltOM9QW1ZdD\n5Yfpznm/PBIePD1WNdnQLrLqs1BKzVVKjVFKjVRK3WGN3ayUmm29X6iUGqqU6q6U6qeUGm+Nv6OU\nmqCUOs56/WM2z7MgeO9BXQdq60e6PWqvEG/WaXV/P91owNYsSvhO0/d5MzyBO17b6veN0d7LoYpy\nBEfv5TQbGA2pKGelOpxGSuLGB1b0hCufhy//UdcF2vJhbGNjF2ty05Ud3M5eEbvWwEnfgS/cChf+\nJslRvZW+bKNPqxIq2zR/Tvkuq2Q4c8InEbAModscwmInPVnX0BNq39OVjrd+BKteSvucDN6YBkaF\nwpaP4j/3HMyGhr3YikS9lYy3X5XzYmQyL0YmI3tTl1Coqgy1ubvd9CljufGZpdQ3x0wFgpWUN8KK\nrvroGdi8JHbQga06wqar0JVDZxMbC/UbBZO/BcD0KRGmP7WE5nDMv1AclLS1SptWz59+Izmn/hco\niAqLHY4yITtUbza19GK8U2HZshSOm9aq8zIkk8vRUIZMoJT+278Zhjki03sNIWKV7SihJapZ7Cfm\n1M522Y2qyhBfPiH+QaGAp9+vi1We7Z3gN1n/NjxyQVy3tJymKzq4V74Esy7XnRGdVBwR/znRD91B\nwV69y/UKxzZDOSsHbFV9OFhilaMp7qbL6e/f3LUj0XIEIyzynYc+B9X/qfsBDHA0JOo1hHMm6igW\nhRCxnIUHVExAtHaV2BZe+yQ5PyauxWqisHj+f2D9WzoxrCvQFYXFY5fAJ8/D8tkgjkdEn5iwmDlv\nBc2R+Adwc6RjKgbbkbVFaI3UWXLmUHFfjhplRdn3OEzXpjqwDe49Bl6/K+vnls8YYZHPNNdrE86S\nx/Vnu+wGQK8hnHTKGQDc1TKNg+gbrpSEPhZZxsshGh3vO0K/DrH6G9gRRdu7SBnzOAd3F1jdOs9x\nx0ro71gw9I5Fwqe8bllk9yE9R4uswg+bVH/+3vJ5VkVC3HHxcRw1ygrlLSqDst56cbGvFl6/M+vn\nls8Yn0U+4wx3BOjuSFzsORiKSjml9Gk2NTQyRmoB6O3oaTFz3oo2+yPSZUhFuWu8fdQENuZcOOk/\n4ZiLtfnJbs/aVZKuImEoLom9z3UOJAQ1dO8Pk2bq/IaSeBOl73XLIkGRuFyMrfThJy3fJCgB1lSG\nYKnl8C4qiZVIN7Qbo1nkM7s/jf9cXgEnfEM7Kou0Q/vU0VqAVIc/y/LIMP4SPie6e0esEr0qzEbH\ng8Vw3l1w+GQoixUeZMeqLrJSt81Q0jXMUIlCuOcgOOlaOP4/4oZTXrcsYguK6c3XMic8mQgBQAuQ\n6pq6WCfHz98cK5Fu05xHvVI6GCMs8pk9CcKirDdccC9cvwiA6po6Zi/WsfQ76M15TXfxqRoU3b0j\nVoluPgvPcftGH3aK7sOxrwu0X1URbfeXQNcSFoefpF8TndoWrbpuGcZO2nsyfAbXNX8/btuNzyyl\n+tNSuGUPjPpCsmaxa13Wzy9fMcIiHwk362ztxJpKZRXaO2h5CGfOW0Fji/sDrDXJde2hVbZv20E/\n4cvWTjVZOqsMosJdS1jU1ejeJUNP1J/7uAuLzvRZuCXz2USDI2wveKJmsWsNNOxNPtCQEiMs8pHH\nLoVfDE7WLBJuHL8buzXJde3BS3txHZ/6G7jkEaj8OpT2gie/AU98XZejyFVUBALB3BYWh3bpGlB/\n+xIsfhSOPAPOuBHOm6kLP7pQ0a3YdbwjtFE7mc+LuHmdqFl8MgfuGgZLn8rS2eUvRljkI2te1a87\nVunGMTYJN47XjR2qKO8QQQGttH0fNh7Gf0k7Li9/UofVfvwsvPVr2Losef9cINIFzFCL/qRrQK15\nFcZ+Ec69Syc+nnSt7mmSQHVNHQcaWpLG25KU11aqKkOeNaTi5nWZlbAXLNHtgFdYjZUW/SnLZ5h/\nGGGRb4QdN/H2T2Ihp6Afsg6mTxlLSTC+HHRHmZ9s2mz7HnYyXP0ynGrZrHO14U1X8FlsWaqd8F96\nEL76N+3U9sEtxwKge0lRhy0yIM3aUuUJ/TUadOtemrpY2ZgcwAiLfCMx9NGlp4BNVWWIr5ygk97a\nUtspE3iZwnwb4Nj0GKAjXorKYFuO1pnsCsJixyrtDD5umjaZpcDrmtld8DoK2xxlZ3QP7l2WPH8H\nWp2cB02AXo7xcLJmZPDHCIt8w1kpFLSD8psvwQ3Lk3atrqljzlJd0uGwXmXu/SeyjJ+N+6bqNLrk\nBYK6b0GumqFsB3cgR4VFJKzLefcfnXpfi24l7gLFy4+RbWxfdsQtlLrHQLjmFbj0r/HC4sCW5H0N\nvhhhkW98XB3/ueIIGHYS9BoSN1xdU8eNzyyNrga37GvQYYc1HRuOOn3KWLz6oj26YEN65zNwvBYW\nn74DLY0ZPb92k8sO7n2bYeU8nejYPz3T403VSznY5J5c2NFpL/Yc3mNldG/d1+g+h4dO0v4t5z1w\naGfuLjByFCMs8oWty3Rl2X//H5Q4qrJ6hD7OnLcirtorJNRk6iCqKkO+9ef+95kPfbZaHDYODm6D\nP58Hz3wrY+eWESI5GjrbuB9+dzLMukx/HuQdXWRTXVPHows2eG7vaDNUq+dwYp2xl36q+7wcyH5u\nSD5ghEU+EAnD70+BB07Vn8911MAp7+N6SGfGySdSXuw9DQ81p/GAdVbTXf58bsXRKwWSg5rFmtdi\nzt5gKRx2TMpDUi0kOtoM1eo53MNqqjXxchh1Nqx5BR45H56+OktnmF8YYZEPJCbfHXkm/PcH8PVn\nPQ9pVX5DlvFKDEyboZPgW6/CtMe1j2D9/MycWCbIVQf3+rehqBz+8x347/djPTd8SLWQ6GgzVKvn\n8BGf1RFfX7wHRsfK2rDujSycXf5hhEU+kFjPp9cQ6DdSJ1d5MH3KWIoTwmY7Mk7eiUsUZusJnQAj\nz9IVXmvfy8AXZggVtrLmA7lVSPDT+bre1mHjoeLw1PuTeiHR0WYot9DZaPMsNwIBHfFVXJZsng13\n7Ll3RbIqLETkXBFZISKrRWSGy/bTReQDEWkRka8kbLtSRFZZf1dm8zy7PInlutMIf3xy0Ya4LmcA\n4Yw8tVtPULxc3BAMSPpO9+Iy3YZ1Yy4JC6eDO0cKH9bv0T6uI05t1WGpigR2tFbq1Tzr0QUbUkfS\nJYaU71qb2ZPLQ7ImLEQkCNwPnAeMAy4TkXEJu20ArgIeSzi2L3ALcBIwGbhFRNyN74XOxoXwr1ug\npxXpMeyUlIfcVL2U+WuSm9hHFPzsuY6PELnsJO+Vbbi1DXUOPwnqPsidlWLUwR3MHTNU7SJA6Si5\nVuCXKNlZWunT79e6jqeMpOt7ZPznh86ElqYMnln+kU3NYjKwWim1VinVBMwCLnLuoJRar5T6EEi8\ni6YALyuldimldgMvA+dm8Vy7Lkuf1K8XPwjf+xCuSF3z5vF3N3pusxvLdCS3V03gipOHeW5vldP9\n8MnQUg9/vSjZl9MZqIjl4M6REuUrX4Inr9ICLDSpVYf6XYeZXzmuw3N0AOp9AiB8FxkicMMn8ONP\n4dhpuorxp29n4Qzzh2wKixDgfCrVWmMZO1ZErhWRRSKyaPv2Agt/azygmwG99yCMOF3/9TnCtZZP\nIuFcMYc4uL1qQnq1flIx+hw4bIK2yb/4vxk6u3aQSw7upkM68qdpP0z+tq7/1ApyoZZYa0hZBaDX\nYF0O5PyZ+vPfvgT/+ln2T6yL0qUd3Eqph5RSk5RSkwYMyH7TlZziw3/odpEQX/+pnVSUd04WLlhO\n90CC0z3QSvNGSTf4z7fhlO/Cyhdhfydn6uaSsFg+Gxr3wVVzdEOpVpJLQRE2AW93l68vLI6yXrE+\n6W/fkxsaaQ6STWFRBziN0UOtsWwfWxgsfkw/gMZ/CU7+r4x97a1Tx2fsu9pE4v2d5v2exPFX6kik\nxX9v7xm1j1zK4K55FPqMaLVjO45EpbSTldSvneRtvmyVBn35k1o7B/jwyXaeVX6STWGxEBgtIiNE\npASYBsxO89h5wDki0sdybJ9jjRlAR/vULYLzfqn7O6SoEppIH4/kqYry4k41J8yctyIpQqs53EoH\nt03/UTqu/oO/dW4UUsQROqs6MXR232adW3HctFgxpVbiVm22ubUBCBnm9qoJdPeqVdUaLXnkWXDl\nczD8NPjgL7kTIJFDZE1YKKVagOvRD/nlwBNKqWUicpuITAUQkRNFpBa4BHhQRJZZx+4Cfo4WOAuB\n26yxwqRxPxzcCQ374DeV8MezdZ3+4y5r9Vf59SLobK0i41nlx14Ku9d1bg2gXDFDLZ8NKK2JtpFc\nyvp3cseXJiSZLwEONrW0vtbZZ66DvRvh5/2h+jo9tmdD7oQ9dyJZ9VkopeYqpcYopUYqpe6wxm5W\nSs223i9USg1VSnVXSvVTSo13HPsnpdQo6+/P2TzPnOfXR8H/TdJFAu148BOvabWDEuAn/1yaE70I\n3Mh4VvnY8wCBT55v+0m1FzsaKhDs3AfOsn/qgosD2u5fyKWsfydVlSF6lCVnoDeHFTemU1vMyegp\n+sS77W8AAB+5SURBVP8JdNfATTVw3wSYOz0DZ9q16dIO7oKgfrdu1FK/C9a9Bd0HwI/WwZk/afVX\n+VUM7ejsWzfSambTGnoM1KG0r98Js/9bd63raKI9uDspdHbPBnjgNNjw73ZpFWA1yyqKf2R0dLMs\nL/Z4hHzXN0fSK3VvEwjANf+Cr1sW8w/+pl8XPtzOM+z6GGGR6+x29NFe9ZIuz9Ctb1pZ2on45Vd4\n9SjoSOyMXNugEBT48gmh9mk8J16jXz/4K6x9td3n2Go628E9/zew5UMo7q79Fe3Eae7p0624w5tl\neVHmU4zSb967UtINhp6o369+uR1nlV8YYZGL7NsEL94IzfWwxyEsGvb4dr5LhV90yCEPjaMjqa6p\n4+n366IBNmGlM3GPvHFO61aHTo69FL5box/WG97N2Lmmhf3/3Zk+i7pFuirv/3yUdg0oNy5/+N98\n/x+L4zTThnQqAncQfsUo25RXVNJNN0va4yjJ3nQItn0C/7q1IH0YRljkIu89BAt+Bwv/GK9ZAPQd\n2eav9Ys7z4Wp79afAHQZkrTq/XjR90itkXV0gUG7cGBnCYuWRt3j5PDJWhttI17lYTqj/4kXfmXN\n0s63SCSxJMieT3UPkLfv7fz8nU7ACItcxE4Q2rNBT9CAw3nXr+3Cwq8GU5tvqAySKqqm1eYEJ0Mn\n65pIHVn51RYOtrDo6KqzWz6CSDMMOb5dX/PYu94Njzo7EsrGb/76zXtf7HutqEy/7l5PdFm1f3Pb\nvrMLY4RFrvDR07DICvqyJ+L+zVqzGHh0bD/n+1Yy6Qjv1WWbb6gMkiqqpl1lSg6frAMFNi1u+3e0\nFuXULDqhkKBtbz8idXFJL6pr6nxX7Z3VdzsRr/kr+M97X+xWs4dZ0VHOyrT7N0O4Beb8AHasbtv3\ndzGMsMgVnvomPP99/d5u87i3Vk/QviPhqAugvC9UeGespsLLZFBeHOD2qtRtNbNNVqNqRn4eSnvB\no1+Cta9n73ec2MKhox3cSsFrv4DX79KNsOwOcW0glZmp0cVs2Bl4JecpUv8bPBkzRb9+5noo6Rlf\n+v7AVti8BBb+QUfaFQBGWOQakUis3eXu9doM1fdIuPSv8MOVvoemwstkkCuOyqxG1fQYAFe/pJ2W\nj321YxL1Es1QHeUUXfA7eONuHf10ySPt+qpUZqa02t52EF5BGm02lfUbCTdtg2Muhr7DY7XYQLfu\nrbf9OLng8cs+RljkGge2xHpI1++CSIsWFoEgBNun8udqUlW6tLvI4cCjdfy8BHTwQLaJc3B3UJ6F\nUvDv38GIz0HV73VV1XbQVeYGZGl+F5Xq1z7D4dDO2HjDXp3pDTGfRp5jhEUu4KxDU79bT8RgSWws\nMSqjjbhVDc2VpCqbbj7x8m0q35BIjwG6/o9zlZgtoppFB5qh9myAfbVw9IVtrgHlJFV3vM6sUpxI\nxpM6nfQZEXsfKNbdBg9ZmkVnF4jsIIywyAUa98feN+zTwsKugAntcmonUpSjSVU2v7j4WM+y020u\nKpjIsJN13/KDO9r/XX44zVCBoNYSs03tQv16+OSMfJ1vd7xA59cTc1JVGeLOiyfELTj8kvVaxWGO\nf2fFMH2P2haAlsbM/EaOY4RFZ7N7PTx7fezzyheg+ZBuDxo6QVfDbEeMvE11TR03PrM0rrNYrvgq\nnFRVhrjn0ome2zMSqjnsM/p1xVxobmj/93lhC4dA0IqGyrIz+OAOXVm2uFusvlE78e2Od0nndMdL\nhTM/b/ehZm58Zmn7NdKRZ+kAk8/9WJv2Nn0Q83tt+VA32gp3wGKgEzHCorNQSq9MHv0KrJgTG5//\n//Rr78PhW6/CFc+0+6eqa+r4wRNLkhLecimpyklVZSgzXfO8CFl5B7P/G/58XvYcz7bPIlCk/7KZ\nZ7FrrS549/6fdahsMLmwXlvoat3xZs5bQVM4fhFU3xzmB08saZ/A6DEQfrQWzvxf6DFIL/LWvKK3\nNR+CBffD+jdj2kYeYoRFZ7H0KbhrGOxc5b691xD92k67c3VNHdOfWuKZo5ArSVWJeNnKdxxobP8q\nsagUTvuBfr/pA22SygZOzSIQzK6wWGFppKDLbGeAm6qXerYmTeXL6Cy85nNYKaY/2U6BYd+L+zy+\n4/HL9D1dv6ftv5HDGGHRWdi2ZS8GZSbv4WfPLUtqKOQkV6NdvGzljS0RbnhicfsFxudvhu9aCXrZ\ncnYrp2aRZTPU+vk6YufWvdpk0k5uql7Kowu8M7f9fBmdid98bo4obp2dgZBpe6GRSItl0ty3qf2/\nkYMYYdFZ+GkM/7s5I34K0DZbP4b3y01h4afxRBSZuen7DNcmhQ0LsmOKiobOWj6LbDm4I2HY8I7u\nDJgh/Ep8QO5qpKnm855MlOIfNxW+4tNi5+C29v9GDmKERUeyfr72SSgFB31WZiXdOuyUFqzd3WG/\n1Rp6pwjJzMhNL6Ijo5Y+qTsQNh1s/3c6ifosgtnzWXz6DtzWV4dcD29Hb+0E/Ep8QO5qpB02n3v5\n+GsO7oD9W+HV2/Oq4GBmvGCG9HjkfP3abzSsdGkpftoPoKT13e/8qCgv9n2wtqveUhbpsLqGx39d\ndyDcvU5fk2Muztx3d4TP4v1H9GuvoTDm3Mx/vwe5lJvjJNV89uo/32qGTNTNpEp76l4pTt57CHau\nhjdn6hyqs3+Wmd/sZLKqWYjIuSKyQkRWi8gMl+2lIvIPa/u7IjLcGh8uIvUistj6eyCb59nhzLpM\nF7VL5NTvw2k3ZPSnLjhusO/2XKg264ZX5zMbv+S9VjHq8/CTrTrRavMS2LkmcyapbPosGg/o1eum\nxTD2fPj+0oyZLtPxB+ViJBSkns9fPNb/fkibolJdSsXN9LfxXd2dEWKle/KArAkLEQkC9wPnAeOA\ny0RkXMJuVwO7lVKjgHuBux3b1iilJlp/38nWeeYUpT0z/pWpHJG5UG3WjVRmjtLi5KJxbaa4DA4b\nB/Pvg98er81SmcDWLLLhs3j6Gvj1GNixAgYfp9uBZohU4dReSZO5QKr5nHHHfPd+/tvzKDIqm5rF\nZGC1UmqtUqoJmAVclLDPRcBfrPdPAZ8XydGlbntpOpR6nyz80/0ckVecPCwnqs26MX3K2LgWnomk\ncty3msHHxd5/8nxmvtPu+Z3pPAuldPKmzaBjM/O9Fqmc1187qe2Vj7PN7VUTuOJk7/PLuGO+PIU2\ndyB/nN3ZFBYhwNmtptYac91HKdUC7AVsUT1CRGpE5A0ROS2L55ld9m/RE6aDm6VU19Qx8WcvedbD\nDFWU56ygAG3mmHnJcb7ys82d89wY9YXY+51rMvOdUZ9Fhst9HNga/7kd/Src8OvHnssLDJvbqyZ4\nJnUqMjxvSnvFf04UHva12r4yu9UCOoBcjYbaDAxTSlUCNwCPiUivxJ1E5FoRWSQii7Zvz824b35/\nKvxqdHIiz6V/jdV/OvMn8I0XM/aT1TV13PDEYk/HdnFQctZB6aSqMsS9l07ES148umBD+/MtbI6e\nCpc/BSdeo52Ti/6s+0G0B5XhDO6DO2H2d2HZP/XnC+6Da15td2VZJzdVL43rs+2korw45wWFjZ9m\n+uiCDVz+8L8z80P9R0HVA1DWW3++bBYcd1ls+4Gtum/3/Sdqh3cXJpvCog5wGhCHWmOu+4hIEdAb\n2KmUalRK7QRQSr0PrAHGJP6AUuohpdQkpdSkAQNyM6OUQ1axuo0J/Z9Le0FLk34/7DNwxGcy9pM/\ne26Zb+hj95KinHVQJlJVGfLtFjD9yQx1vhOB0Wdrk05Lg25E9fqdsUZUbSGqWWTIwb36ZfjgL/Ci\nFSsy/LMw9IT2faeD6po630S8jIQrdxBVlSF6lHkHe85fsytzC42Jl0F36/lT2jNe22g6oGuQgdWW\nteuSTWGxEBgtIiNEpASYBsxO2Gc2cKX1/ivAq0opJSIDLAc5InIkMBpYS1fGTVic8WMo7Q2DM2tz\nTmXP39uFbnrAU7MAaI6kF72TNv0T1iRb22GysDSJSx58j9+/uZ7Gpqb2mUCcbT2RdnVNdONnz/kn\nOuZq5JwXqe6DjCR22pxzuzZB9TkiOVDFrtbQgflT2SBrwsLyQVwPzAOWA08opZaJyG0iMtXa7Y9A\nPxFZjTY32eG1pwMfishitOP7O0qpXXQ1nCGYn76jBYNNaU9dluHGDTEVtoPIlb7J6ZIqkDWjN32i\nsNjV9jXKX9/Rvo8mJbQQJEiERxdsaLvAcJ5LeUWsMU+GSPVwzdWcHC9SybaMakpjz4Mfr4OS7jFh\nUWaZBz99R7928ciorPoslFJzlVJjlFIjlVJ3WGM3K6VmW+8blFKXKKVGKaUmK6XWWuNPK6XGW2Gz\nxyulnsvmeWYNZ/Zm037o5YjxzkKYbLp0sXs+JRm96RNDIfdsgDd/BYsfb/VXvbVCOzfDBIgQoEgi\ngOLxdzf6H+hEKZjzQ903fNe62Hivoa0+Hz/S0c68nMa5SqfN8zLLDGUXA7VzLRr2dOmqtCaDO1u8\n8GN418ol7DkE9m+Kb4va47Cs/Gw6N31XM0P16VacctVbXVOXOT/MF3+tbdD/uhU+/TfUWibEY7/a\nqnyGgOWjCBOkRenjAqjWrdBrF8HCh/VfeV844Sp9bqPOTv870iCVCQpyN2u7rWQssTOR3partthh\nduo+ENa9qavSXvlcfHOzLkKuRkN1bfZtigkKgBO/qV+HOZzYGUyicpJOf4pcrevjxS0Xjk9qB5tI\nRk1RJ14D4y7SPoFah69pr39xvUS0JqE1i7B1qxURbp3tf7PDgV+/C/qOhLNugmEntepcUpFKGF9x\n8rAuExRhk6rla0YTO530tCwIkRZdGQB0pQCbl2/Ozu9mGSMsMsna1+Hvl8DqV+LHDz8JrnkFzv65\nDrO78DdZ+fnqmjrP/gM2udZzOx2qKkPM/Ip/zkVWInUSHcjb028UVV1T59AsAoTRD6YAEU4+sk/6\n5+A0PQH0HeG+XztIpY12hdwKN26dOj5lYmfGQmidDDwaTvwWVP0evrdYl8J35l9sqtHPikcuiPXx\n7gIYYdFeDu6AW3vD8ufh75fCqpd0fH6544FQ3geGTtJlJSZeBidc6f19baS6po7v/8M/jDRUUZ5z\nPbfTxc658OOkO17O7I/2G6Vf7cCE7Z/o6q6NLnW9Epg5bwUB3DWL9TvTyCLev1UXodu9TpswbPpn\nVtCnM2+6oqCAWGKnnyY3f82uzM+bQBC++CtdQqb3UC3g7VyYIZX69a8X6T4qSx6HJ78BD+R+3rER\nFu1lh9Xp7o27odgy7+yrhRGfi+1TlrmkKS9+/PSHvtu7FQeYP+OsLikobFKd+9b9TZx9z+uZ+8HR\nU/TrKf+tfUzbPoHfnQK/PyWl93TTnnqClrBoIRgVFgEiqUtO7P5U1316+RYdATX0RC0kuvWD/qPb\n/c+ySUdQpDLl5DpVlSF+felxvvts3d/OkOZ0sCOhRp8THxW5YxUse0b38c5xjLBoL4379evB7bH3\nAEdfCFc8rR1ZPQdl9RSqa+pobIn47nOo2X97VyFVielV2w5mLu9i4FFwwydw+g9hwFh9U+/fBHs+\n9W6taTGkopyg5bOIKKdmEUntM7J7Oy+4Xydy9R0B334Drl+U0fph6Ti1b506PmO/11mks0DyS0bM\nCMd9VfssJ39bJ1PaOMOhI7l9jxph0VZ2rdPmJ7v0wv7N8Rm6R5yq6w1d+Vx8FFSGual6acrVYT5x\ny4WpH17f/8fizK0Uew3WD+gBR8XaZgJs+cj3sOlTxlIiTs1C+yxKgyq1z8hZm6qlQQuL4vKMlSAH\nPW/SKcbYlTXR1jL+5hczm+DpZPBx8M0XdWi2s0nVJse9++x18Jvjs/P7GcCEzraVJbOs18fix6c9\npnMoemWobr4Plz/8b+avSc9B1tXNCTZVlSF+/PSHKTWpRxdsYN32A/z9WxkqozIg4QG/c7Xv7lWV\nIQLv94Ja7bNosYRF75I0fivxu/se2YoTTc3/b+/eo6OszwSOf5/cCASUm0SJdOFApBXbcltYxFV7\njhVsawWlgLp79JR62V23WrqcQsWVU8WCPUW71q7VSuttEasYWVPFFCsWFrrgJhZBQUBAwr3cKuQ6\nefaP3/smk8kkM+9kJjMhz+cczsy8729m3hnezPP+bs/vq0ve4ePDsVcFbCt7a2eTJbFX/ztdG2q8\n8EppkBw7yzUr7lkHa5c0bfd/S47vcTPBM4zVLIII1bvO6x2roa6VP7aisR0yhrqkvDLuQJElZ0dz\ngm/xDfGlR0lq/h+/Yzm3wLU5H//E9VtURV/Gs6S8kq17XV6pGnKpUResq6vPMG/F5pbHFf5af9nh\n2rZ9hZck5zPgahTxBIrCXnmdtmM7miBp1ZOWb6w1uflQfBX0Gxp9/6n98D+PNbVaZAgLFkEc3e6S\ny/333dF/JP72NuiVmsl24W5+an3cTU99euSyZPrIs6o5YcqoIh6d0fbIKN89yyuSMzyyaLRbtvTW\n16HvYNcM+c4iWDzETZyL8JNV28gO1QBQQx7VuCpFN+qoqgu1nA+zZjE8PNSlhji+G87/IlzzMFz6\nXeg5gPYqKa+k+IelcbXNZwF/uje5k/7SzV/nIp4en7oGGDK3NPWd3iOuj57q58xReGs+/PbW1L5/\nQBYsgqg55W5Pfto8z8vA0TD/sBsul2LjF5bFXaMAKP/3q8+qQOGbMqoo7vWU1+081v5RUnkFcNNy\nFzT6DHEdk+sfBxQ+eKVF8f0nqugmtYRUqCObGtyx5lPbuL+Zihdcn9fqH7nJXP2Gwfg74OoH2nfc\nNI16ineMw5I4A3Fn8+CUL/LJoq/HVVZJcirzaPJ6wA/2uCVxw53cl7r3bAcLFkGEr5sdvgJWz8Kk\nJ3WLNL9kM4PnlnLor7VxP6ez5fIJKp7Obt/Hh08zeG4pg5Nxxdh3iGuGqvVGv+3d0KLIwN7dyaeO\nGvIAaVaz8Pc3qvms6Qdir/fj1LeVJooA/HMm6ACIs/HiIlHrdh5j+Pw3UtfxLeLSAYWL0R+WLhYs\n4vXh67DzD02Pw9NARHZ+Jtn4hWUJDe3rbDO1g5oyqojiAQWBn/f8hr3tG/kSvozpxVPg8FbXnxWq\nbxz+OGfScAqy6qn2ahR+n0W+1DbNovfXMznwPmhD80VzCiOXq49fSXklX7jvjYTOmbOpU7s18dZI\nfTX1DclrzowmO2Kc0cZfNd3ftNStspcBLFi0ZcfqptXNlt8M63/efP+Eu9zt30wkmfwlUf0r4SC1\nCV9nzOWTiLLZVzJxaPAhpf7Il4Sap4ZcATnd3e3wa9zw1qPb4Tdfh19PhoYQU0YV8fn+OV7Ngsbb\n7lLPDWOKmDKoyq2g+MefuvQPAJfPgbyebjx+AlmJ/X6Je5ZXUJXAvJrOmtYjqHhyjUWzbucxBs8t\nTU1/hj8B9OqFzbe//j146ivu/qkDcCB9k/ds6Gxrdq+F56+H8f8EF01quX/qL12yuZE3QWHyRhrF\nO6yxNX165HL/tSO6RKDwvXDbBOaXbE7oStpvngLolpPF4hu+FPu7K+gH3//IrYB3Yo/btulp+NRr\njjryESX7e5Nz9Dh9vZqFX8PI0xpeea+Sm6rXcXH1CddHMWKqa4roNxTufr/lus5tSPRzR3p0xtk1\nCKIt/uf8yaptMXOpReP3Zzy/YW/yAuz0Z13z0/mXuAuIqrB+ydrPYNca93vUUA8L0pPmXPQsWdxg\n7NixumlTy1EpCdv2JiybEX3fuYPcbNrc/Ha9RUl5JQtWbklaErziAQWUzb4yKa/VGbU30EaaOLRv\n7HkaoTp4qAi8kU8AfPPnTHyriPtOP8RgOcjk2sVcKEdY2+1u/q3uDl4OXcFzPR7h7xs2Nj1n+Nfg\nxthrZpSUVzLnt/F3Vscjrs95loon5Um8BLg5GcFj91pXS23N/SeSOpOfthejbGQ1i2jeWQyHWpmh\ne+E4+E7iiceSdSUYqSv/wfvKZl/J/JLNvLBhb8zV9eLhNztEU9S7O3MmDXdXqYUjYP//ufk1+ytg\nfzn7T/QlP7e2cRRUY58FtYAyPPQxjJgCW0vcC/oJ5nA/YIle9QaRJW7+QVdoemqNX8uYt+LPCTXd\nhQuvcUCAmmqkwZe5UVKLW5mY97s5bhh38VXtOt6grGYRqfY0PDSw5fZbS2Hn2y63SxtzKUrKK5Ny\n4sXLgkR0QWa3t9d1WWtZkPssK4Y9xKyGl6H6JKMOzeep+nupJ5uZtfdRQBVb8mexuG4mr4YmsiH/\nX908iuqTHF/zC75xZgGVnNchxwt23kRTUl7J7JcqYs70TpaY/wefvOtmczfUub6LSMlrjrKaRWBn\njkUdM/+z+qk88sRJYAyUJbGpqx3sqrBtfj/Gsj99mvK1o19ruIzXai6DLVCd05vbs9cxvX4FY3O3\n83LIzeY/jRsq+4PcF+kv7o98akkV5ToCeDylxxcuW4Qbxw+y8yYKvwaQzKbhtrRVc23iBm/sjtbi\nvWuNq80mt0mqVSmtWYjIZOBnQDbwK1VdFLG/G/AsMAb4CzBDVXd7++YBs4AQ8F1VXdXWewWtWURe\neQ6SQzyY82uuyG452uDL1U9ykp5xv3Yq2RVhYjqqpjFatrOi24LGx34fBcAvc5cwKdudo0f0XC6t\neYy6Drpes/MmuJLySmYvryATcsF+M2sd/5HX8qLih3WzKAuN5gjNF9QK+P8dV7RJWbAQkWxgO/BV\nYB+wEbhRVbeGlfln4EuqeqeIzASmquoMEbkYWAaMAwYCvwcuUg1P69pckGAR/sPRjVoW5i5lWva7\nLcqd0h6MqXmiw/6gW5O0jjMDpDpwKDdlv00+tZSEJnKMXvh/i92oZYAcZ2b2H3gnNJKN+vkUHYPT\nVYbCplqyB6IkqgfVbM3/dtR9d9bew5sN45ptCxAw0h4sJgALVHWS93gegKr+OKzMKq/MehHJAQ4C\n5wFzw8uGl2vt/YIEi8FzSxkt23ku78cUSNMolo8aBnFYe1NPNk/UX8tRzmWXRum/6CB2NZh6qRpw\nkA4Jd6iauKViNFoQl2e9T43m8WDuUsobhjE9Z03jvpAKh+jDpTVN88F2x5feJO3BYhowWVW/4z3+\nR2C8qt4VVuYDr8w+7/FOYDywANigqs97258G3lDVlyPe43bgdu/hcCCuRZLzzh82ph0fLaVCVX89\nEjp5qDP/evUHjqb7IBKV07eoOCuve/wTHdLMzpf0yep+Tt/sXv0HSVZWxvb91h7c8V4cxY6q6uRY\nhTL2Q8ZDVZ8EnkzHe4vIJlUdm473zmT2vURn30t09r20lKnfSSrTfVQCg8IeX+hti1rGa4Y6F9fR\nHc9zjTHGdJBUBouNQLGIDBGRPGAmsDKizErgFu/+NOBtde1iK4GZItJNRIYAxcD/YowxJi1S1gyl\nqvUichewCjd0dqmqbhGRHwGbVHUl8DTwnIjsAI7hAgpeuZeArUA98C9tjYRKk7Q0f3UC9r1EZ99L\ndPa9tJSR38lZM4PbGGNM6liKcmOMMTFZsDDGGBOTBYuARORbIrJFRBpEZGzEvnkiskNEtolIlEUw\nugYRWSAilSJS4f37WrqPKV1EZLJ3PuwQkbnpPp5MISK7RWSzd35kRsK1NBCRpSJy2Jtz5m/rKyJl\nIvKxd9unrdfoKBYsgvsAuB5olh/ES1EyExgBTAZ+4aU86aoeUdWR3r/fpftg0sH7/38cuAa4GLjR\nO0+M8xXv/Mi4OQUd6De434twc4HVqloMrPYep50Fi4BU9UNVjTZT/DrgRVWtUdVPgB243Fam6xoH\n7FDVXapaC7yIO0+MAUBV38WNBA13HfCMd/8ZYEqHHlQrLFgkTxHwadjjfd62ruouEfmzV83OiGp0\nGtg50ToF3hKR97y0PaZJoaoe8O4fBFpfQKcDdep0H6kiIr8Hzo+y615Vfa2jjycTtfUdAf8JPID7\nQXgA+CkQPV2m6aouU9VKERkAlInIR95VtgmjqioiGTG/wYJFFKqayHqFXSpFSbzfkYg8Bbye4sPJ\nVF3qnAhCVSu928Mi8iquyc6ChXNIRC5Q1QMicgFwON0HBNYMlUyWosTjneC+qbhBAV1RPClvuhwR\nKRCRXv594Gq67jkSTXgapFuAjGjNsJpFQCIyFXgMt+5GqYhUqOqkTpKipKM8LCIjcc1Qu4E70ns4\n6dFayps0H1YmKAReFbccaA7wX6r6ZnoPKT1EZBlwJdBfRPYB9wOLgJdEZBawB5ieviNsYuk+jDHG\nxGTNUMYYY2KyYGGMMSYmCxbGGGNismBhjDEmJgsWxhhjYrJgYUwAIhIKy6Zb0VomWRF5VEQub+N1\nbvGGTYZv6y8iR7y5Oi+KSHGyj9+YRNnQWWMCEJHPVLVnjDL9gFJV/bs2ypwD7AI+p6pnvG13AuNU\n9dsicgXwD6p6WxIP35iEWc3CmOS7AWicZCYiY0RkjZc0b5WXyuEUsAa4Nux5MwG/tvFH4CoRsYmz\nJiNYsDAmmO4RzVAzopSZCLwHICK5uBn/01R1DLAUWOiVW4YLEIjIQOAi4G0AVW3Apbn/cio/jDHx\nsqsWY4KpUtWRMcpcABzx7g8HLsFlVgWX9sNPP12KWyTrHFxKh1ciUsQcBgbiBR5j0smChTHJVwXk\ne/cF2KKqEyILqWqViLyJS7Y4E5gdUSTfey1j0s6aoYxJvg+BYd79bcB5IjIBXLOUiIwIK7sMFyQK\ngfURr3MRlo3VZAgLFsYEE9lnsShKmVJcJlG85VSnAYtF5H2gArg0rGwZrqlpuYYNTRSRQlyT18EU\nfQ5jArGhs8akgIisBb6hqicSfP73gFOq+nRyj8yYxFjNwpjU+D7wuXY8/wTwTJKOxZh2s5qFMcaY\nmKxmYYwxJiYLFsYYY2KyYGGMMSYmCxbGGGNismBhjDEmpv8HjhEwdHH1RwAAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f6a0c323f60>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"model,posmin,posmax=mymodel(24)\n",
"kpm=pb.kpm(model,silent=True)\n",
"ldosz = kpm.calc_ldos(energy=np.linspace(-12, 12, 1000), broadening=0.05,position=[0, 0,posmax])\n",
"ldosmz = kpm.calc_ldos(energy=np.linspace(-12, 12, 1000), broadening=0.05,position=[0, 0,posmin])\n",
"ldosz.plot(marker='o',label='+z')\n",
"ldosmz.plot(label='-z')\n",
"ylim(0,0.4);\n",
"legend();"
]
}
],
"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.5.2"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment