pymatgen phase diagram

pymatgenpd.ipynb

{
 "metadata": {
  "name": "matgen_pd"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Add your API Key\n",
      "----------------"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "api_key=\"your_key_here\""
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 13
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Imports\n",
      "-------"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from pymatgen.matproj.rest import MPRester\n",
      "from pymatgen.phasediagram.pdmaker import PhaseDiagram\n",
      "from pymatgen.phasediagram.plotter import PDPlotter"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 14
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Initialize the API\n",
      "------------------"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "mp_api=MPRester(api_key)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 15
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Grab all compounds with Ca C O\n",
      "------------------------------"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "entries = mp_api.get_entries_in_chemsys(['Ca', 'C', 'O'])"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 16
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Create a Phase Diagram\n",
      "----------------------"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "pd = PhaseDiagram(entries)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 17
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Plot the PD\n",
      "-----------"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "plotter = PDPlotter(pd)\n",
      "plotter.show()"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAkkAAAHECAYAAAAtRr6vAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XdclWXjP/DPARmCg+VGnEgGrhBzRO6dZZpppmWWaZkC\nIipOUNzI0MqGJfmk1aOZ/czcuVLRNPeCHA85EAEV2ePcvz/O1yuOoTLOOdcZn/fr9bxeXHA490fz\nwY/Xdd/XpVIURQERERERabGSHYCIiIjIGLEkEREREZWAJYmIiIioBCxJRERERCVgSSIiIiIqAUsS\nERERUQlYkoiIiIhKwJJEREREVAKWJCIiIqISsCQRERERlYAliYiM3tmzZxEUFIQWLVrAyckJ9vb2\n8PDwwODBg/Hjjz+CpysRkT6oeHYbERmroqIiTJ48GcuXL39iEerQoQPWr1+PunXrGjAdEZk7liQi\nMkqKomDIkCHYuHGj+JydnR1atmyJwsJCnD9/Hnl5eeJrrVq1wuHDh2Fvby8jLhGZIS63EZFRio2N\n1SpIb7zxBm7evIkjR47g+PHjuHr1Kvr37y++furUKSxevFhGVCIyU5xJIiKjk52dDQ8PD6SnpwMA\nunXrhl27dv3rdbm5uWjTpg0uXboEAHB3d0dSUpJBsxKR+aokOwAR0aM2bdokChIAzJ49u8TX2dvb\nY8aMGTh27Bi6dOmCzp07GyoiEVkAziQRkdEZO3YsvvzySwCAg4MDHjx4AJVKJTkVEVka3pNEREbn\n/Pnz4uMmTZqwIBGRFCxJRGR0ii+1Va9eXWISIrJkLElEZHQKCwvFx2q1WmISIrJkLElEZHScnZ3F\nx/fv35eYhIgsGUsSERkdT09P8fHly5dRVFT02NeeO3cO69atQ1ZWliGiEZEFYUkiIqPj7+8vPs7N\nzcWBAwce+9q4uDiMGDECNWrUwODBg5Gfn2+IiERkAbgFABEZnXv37qF+/fpidqhHjx7YsWPHv16X\nmpqKZs2a4d69ewAAPz8/HDlyxKBZich8cSaJiIyOk5MTgoODxXjXrl14++23te5PSkhIQL9+/URB\nAoApU6YYNCcRmTfOJBGRUVKr1ejduzd2794tPmdvb49WrVrhwYMHuHTpkta9SoMGDcKGDRtkRCUi\nM8WSRERGKz8/H2PHjsU333zz2NeoVCoMGzYMq1evhq2trQHTEZG5Y0kiIqN39OhRrFq1Cnv37sWt\nW7dQVFSEunXrokOHDnj33XfRpUsX2RGJyAyxJBERERGVgDduExEREZWAJYmIjNq2bdvg5NQJVla9\noFINgUr1OlSqIbCy6gUnp07Ytm2b7IhEZKa43EZERmvcuJn4/PPLAJYCcC/hFdcBhGDs2Cb47LMI\nw4YjIrPHkkRERklTkM4C2IgnT3qrAQzG8OENsHZtjGHCEZFF4HIbERmdbdu2/d8M0tMKEv7v6z9i\n3TobdOnylv7DEZHF4EwSERkdJ6dOuH//B5S8xPYkb+Ds2Znw9vbWRywisjAsSURkdKysekFR/n1W\n29Mlwtl5FNLTD+o8ExFZHi63EZHRUZTq5fxOT9y756jTLERkuViSiMgIqSrwveUtWERE2liSiMgI\nVeQugPs6S0FElo0liYiMjkpV3qKTCCenTJ1mISLLxZJEREanWrUsaDaKLAsFwGwcPrxaD4mIyBKx\nJBGR0fn++1kAQqDZKLI01ACmoHNnG3h5eekvGBFZFJYkIjI6ffr0wdixTQAMxtOL0sMdtwuwd+8a\n/YcjIovBfZKIyGi9+WYg1q2zgebstpIoAKbAweF7pKUlwt7e3oDpiMjccSaJiIzWyJF9AKyGZkYp\n4ZGvJgJ4A8B/kJ19HVFRUYaOR0RmjjNJRGSUCgoK0LJlS1y8ePH/PuMGzTEltQAUwtk5B+PHd0NE\nRAQAwNHREQkJCahbt66kxERkbliSiMgoxcbGIjAwEACgUqnw6I8qRVFQWFiINm3a4OzZswCAt956\nC998843BsxKReWJJIiKjc+fOHXh6euL+/cfvl/TwR9fu3bvRo0cP8fn4+Hg8//zzes9IROaP9yQR\nkdGZPXu2KEhubm5PfG337t3x6quvinFAQADU6tJuHUBE9HgsSURkVE6dOoUvvvhCjN3d3Z/6PZGR\nkbC1tQUAHDlyBGvXrtVbPiKyHCxJRGQ0FEVBYGCgmAnq3r07zp0799Tva9y4MYKDg8V46tSpyMzk\n8SREVDEsSURkNH766Sfs3bsXAGBtbY2BAweioKAAANCyZcsnfm9oaChq164NALh16xYWLVqk16xE\nZP5YkojIKOTm5mrNBo0fPx5nzpwR4379+j3x+6tWrapVjCIjI3H16lXdByUii8GSRERGISoqCteu\nXQMAuLq6Ys6cOfj111/F159WkgBg5MiR8PPzAwDk5eUhJCREL1mJyDJwCwAiku7GjRvw8vJCVlYW\nAODTTz+Fv78/WrRoAQCoXr06UlNTYWNjI77ncT+6Dh8+jI4dO4rxb7/9hq5du+oxPRGZK84kEZF0\noaGhoiC1aNECY8aM0ZpF6t27NypVqlSq9+rQoQNGjBghxoGBgSgqKtJtYCKyCCxJRCTVkSNH8J//\n/EeMY2JiUKlSpTIvtRW3aNEiODg4AABOnz6NVatW6SYsEVkULrcRkTRqtRodO3bEkSNHAACvvvoq\nNm7ciPv378PV1VXMACUnJ6NWrVpQqVTie5/2oysiIgKzZs0CoLnHKTExEc7Oznr6lRCROeJMEhFJ\ns3btWlGQbG1tERkZCQDYuXOnKEi+vr6oVatWmd87ODgYDRo0AACkpaVh7ty5OkpNRJaCJYmIpMjM\nzMTUqVPFODg4GI0bNwaACi21PVS5cmVRugDg448/xoULF8qZlogsEUsSEUmxcOFC3Lp1CwBQp04d\nhIaGAtAswW3dulW8rrwlCQAGDx6Mzp07AwAKCwsxadKkCiQmIkvDe5KIyOCuXr2K5s2bIy8vDwAQ\nFxeHt99+GwDw559/wtfXF4DmXqLbt2/D2toaAMp0T9JDJ0+exHPPPSdev2XLlgoVLyKyHJxJIiKD\nCwkJEQXJz88PI0eOFF8rvtTWp08fUZDKq3Xr1hgzZowYBwUFIT8/v0LvSUSWgSWJiAxqz549+PHH\nH8U4NjYWVlb//CjSxf1Ij4qIiED16tUBAAkJCfj444918r5EZN643EZEBlNYWAhfX1+cPn0aADBi\nxAitPZLS0tJQs2ZNqNVqqFQq3LlzB66uruLr5Vlueyg6Olrck1StWjUkJiaiZs2aFfnlEJGZ40wS\nERnMqlWrREFycHDQOpAWAHbs2AG1Wg0AaN++vVZBqqjx48fDy8sLAJCRkYGZM2fq7L2JyDyxJBGR\nQdy9e1ermISGhqJevXpar9HHUttDtra2iIqKEuNVq1bh5MmTOr0GEZkXliQiMoi5c+ciLS0NANCg\nQQMEBwdrfb2oqAjbtm0TY308gdavXz/07dsXgGa5LiAgoMzLdkRkOViSiEjvLly4oHWzdGRkJCpX\nrqz1mmPHjiE1NRUAUKtWLbRu3VovWaKiosRhufv378eGDRv0ch0iMn0sSUSkV4qiICgoCIWFhQCA\nzp07Y/Dgwf96XfGltr59+2o98aZLzzzzDCZMmCDGkydPRk5Ojl6uRUSmjSWJiPTq119/xfbt2wEA\nVlZWiImJ0XpKrfjrHtL3Zo+zZ8+Gm5sbACApKUnr+BIiooe4BQAR6U1+fj5atGiBhIQEAMD777+P\nzz///F+vu337NmrXrg0AsLa2RmpqKpycnP71uopsAfCozz//HOPGjQOgedLu0qVLcHd3r9B7EpF5\n4UwSEenNxx9/LApS9erVERERUeLrit+w3alTpxILkq699957aNmyJQAgOzsb06ZN0/s1ici0sCQR\nkV6kpKQgPDxcjOfMmYMaNWqU+FpDLrU9ZG1tjdjYWDFeu3YtDh06ZJBrE5FpYEkiIr2YOXMmMjIy\nAABeXl4YP358ia8rLCzEjh07xNiQh8926dIFr732mhgHBASIzSyJiFiSiEjnTpw4gVWrVolxdHQ0\nbG1tS3xtfHw87t27BwBwd3eHj4+PQTI+tHTpUtjZ2QHQbEOwZs0ag16fiIwXSxIR6ZSiKAgMDBQ3\nVvft21ds4FiSR5faSnryTZ8aNmyIyZMni3FoaCgePHhg0AxEZJxYkohIpzZs2ID9+/cDACpVqqR1\nFEhJZNyP9Khp06ahbt26AIDk5GQsWLBASg4iMi7cAoCIdCYnJwfPPPMMkpKSAABBQUFPLEk3btwQ\nj93b2NggLS0NVatWfezrdbkFwKO+/fZbjBw5EoDmnLfz58+jSZMmOr0GEZkWziQRkc5ERkaKguTm\n5obZs2c/8fVbt24VH7/44otPLEj6Nnz4cLRv3x6AZn+n4ktwRGSZWJKISCeuX7+ORYsWifH8+fOf\nut+RMSy1PWRlZaW1JcCmTZuwe/duiYmISDaWJCLSiWnTpiE7OxsA0KpVK7z77rtPfH1+fj527twp\nxrJLEgC0a9cOb731lhgHBgaKM+eIyPKwJBFRhR06dAhr164V45iYGFhbWz/xe37//XdkZmYCABo1\nagQvLy+9ZiythQsXwtHREQBw9uxZfPHFF5ITEZEsLElEVCFqtRoBAQFi/Nprr6FLly5P/T7Zj/4/\nTt26dTFjxgwxnjVrFtLT0yUmIiJZWJKIqELWrFmDY8eOAQDs7OywdOnSUn1f8Zu2jWGprbigoCA0\natQIAJCeno6wsDC5gYhICm4BQETl9uDBAzRr1gzJyckANEeRzJs376nfd+3aNVFC7O3tkZaWBgcH\nh6d+nz63AHjUTz/9hEGDBgHQnPN26tQpeHt76/WaRGRcOJNEROW2YMECUZDq1q2LqVOnlur7is8i\nde3atVQFydAGDhyIrl27AgCKiooQFBSk92JGRMaFJYmIyuXy5ctaG0UuXrwYVapUKdX3GtOj/4+j\nUqkQExMDKyvNj8mdO3fil19+kZyKiAyJy21EVC6vvvoqNm3aBABo3749Dh48KArFk+Tm5sLFxQU5\nOTkAgL/++qvUO1sbcrntoQ8//BArV64EADRt2hRnz54VB+ISkXnjTBIRldmuXbtEQQKA2NjYUhUk\nANi3b58oSM2aNTP6oz/mzp0rNsX866+/sHz5csmJiMhQWJKIqEwKCwsRGBgoxm+//TbatWtX6u83\nhaW24tzc3BAeHi7G8+bNw+3btyUmIiJDYUkiojL54osvcO7cOQCAo6MjFixYUKbvN7WSBAAffPAB\nmjdvDkDzRF/xfZSIyHzxniQiKrX09HR4enqKzRUXLFiA0NDQUn9/YmIimjVrBgBwcHBAenp6me7v\nkXFP0kPbt29Hnz59RI4//vgDvr6+Bs1ARIbFmSQiKrWwsDBRkBo1aoSgoKAyfX/xWaQePXqY1A3Q\nvXv3xksvvQRAU9ACAgK4JQCRmWNJIqJSOXfuHD799FMxXrZsGezt7cv0Hqa41FbcsmXLYGNjAwA4\nePAgfvjhB8mJiEifWJKI6KkURUFQUBCKiooAAN26dcPAgQPL9B5ZWVnYt2+fGPft21enGQ2hWbNm\nWufUhYSEIDs7W2IiItInliQieqpffvkFO3fuBABYWVkhJiamzAfS7tmzB3l5eQAAHx8feHh46Dyn\nIcycORM1atQAAFy/fr3UZ9URkelhSSKiJ8rLy8OkSZPEeOzYsWjRokWZ38fUl9oeql69utYTfYsX\nL0ZSUpLERESkLyxJRPREy5cvx19//QUAcHJywty5c8v8HoqiaJUkU1xqK+6dd95BmzZtAAA5OTml\nPrOOiEwLtwAgosdKTk5Gs2bN8ODBAwCanbUnTpxY5vc5f/48vL29AQBVq1ZFWlqauAG6LGRuAfCo\nAwcO4MUXXxTj/fv3w9/fX2IiItI1ziQR0WPNmDFDFKTmzZvjgw8+KNf7FJ9F6tWrV7kKkrHx9/fH\n0KFDxTgwMBBqtVpiIiLSNZYkIirR8ePHsXr1ajGOjo4ud7kxl/uRHrVkyRKxDcKff/6JuLg4uYGI\nSKdYkojoXx7dLPGll15C7969y/VeGRkZOHDggBg/3LXaHHh4eGDKlCliHBoaioyMDImJiEiXWJKI\n6F9++OEHHDx4EABgY2ODZcuWlfu9du3ahcLCQgBAmzZtULduXZ1kNBZTpkyBu7s7ACAlJQURERGS\nExGRrrAkEZGW7OxshISEiHFAQIA4b608zHWp7SFHR0csWbJEjGNiYpCYmCgxERHpCksSEWlZsmQJ\nrl+/DgCoWbMmZs6cWe73UhQFW7duFWNzLEkAMGzYMHTq1AkAUFBQgODgYMmJiEgXuAUAEQlJSUl4\n5plnkJOTAwD48ssv8d5775X7/U6dOoXWrVsDAFxcXJCSkgJra+tyv58xbQHwqGPHjsHPz0+Mt2/f\njl69eklMREQVxZkkIhKmTp0qClKbNm3wzjvvVOj9ii+19e7du0IFydi1bdtW6/crKCgIBQUFEhMR\nUUWxJBERAM3miN9//70Yx8bGVrjUmNMu26WxYMECVKlSBYBmA83PPvtMciIiqggutxERioqK4Ofn\nhxMnTgAAhg4dqlWYyuPu3btwc3ODWq2GSqVCcnIyatasWaH3NObltoeWLFkijilxcnJCYmIi3Nzc\nJKciovLgTBIRIS4uThQke3t7rae1ymvHjh1iB2o/P78KFyRTERAQgCZNmgAA7t27hzlz5khORETl\nxZJEZOEyMjIwffp0MZ4yZQo8PDwq/L7m/uj/49jZ2WntK/XZZ5/hzJkzEhMRUXmxJBFZuIiICKSk\npAAA3N3dtXaQLi+1Wm0Rj/4/zssvv4wePXoA0PxeBAYGGu3yIBE9Hu9JIrJgiYmJ8Pb2Fk9hrVu3\nDm+88UaF3/ePP/5Au3btAAA1atRAcnIyrKwq/m8yU7gn6aGzZ8+idevWKCoqAgD89NNPGDhwoORU\nRFQWnEkismDBwcGiIHXq1AnDhg3Tyfs++lSbLgqSqfHx8cEHH3wgxsHBwcjNzZWYiIjKyvJ+chER\nAM2N1Zs3bwagmaGJjY3VmqmpCEteaisuPDwcLi4uAIArV64gJiZGciIiKgsutxFZoIKCArRu3Rrn\nz58HALzzzjv4+uuvdfLed+7cQa1ataAoCqysrJCamgpnZ2edvLcpLbc99PHHH2PChAkAgCpVqiAh\nIQF16tSRnIqISoMzSUQW6LPPPhMFqUqVKliwYIHO3nv79u2iwHTs2FFnBclUjRs3Dt7e3gCAzMxM\nrScJici4sSQRWZjU1FTMnj1bjGfNmoXatWvr7P0tbZftp6lUqZLWMltcXBz++OMPiYmIqLRYkogs\nzJw5c3Dv3j0AQJMmTRAQEKCz9y4qKsK2bdvE2JLvRyquR48eeOWVV8R44sSJJrNcSGTJWJKILMiZ\nM2e0zhOLioqCnZ2dzt7/yJEjuHv3LgCgTp06aNWqlc7e29RFRkbC1tYWABAfH49169ZJTkRET8OS\nRGQhFEVBYGCgOCqkR48eGDBggE6v8egu27p6Ws4cNG3aFIGBgWI8depUZGVlSUxERE/DkkRkIX7+\n+Wf89ttvAABra2tER0frvMRY6lEkpTVjxgzUqlULAHDjxg0sXrxYciIiehJuAUBkAXJzc+Ht7Y0r\nV64AAD766COsWLFCp9e4efMm6tWrB0Bzs3JaWhqqVaum02uY4hYAj1q9ejVGjx4NQHOY8IULF9Cw\nYUO5oYioRJxJIrIAMTExoiC5uLggPDxc59cofsO2v7+/zguSuXj77bfh6+sLQFNedXFWHhHpB0sS\nkZm7desWIiIixHju3LliF2hd4i7bpWNlZYXly5eL8fr167Fv3z6JiYjocViSiMzc9OnTxQ3C3t7e\nGDt2rM6vUVBQgB07dogxS9KTdezYUesg4cDAQHEQLhEZD5YkIjP2xx9/IC4uToxjYmJQqVIlnV/n\n0KFDyMjIAAA0aNAAzZs31/k1zM3ixYtRuXJlAMDJkyd1diwMEekOSxKRmVIUBRMnThTjV155BT16\n9NDLtR7dZZuP/j9d/fr1MW3aNDGeMWOG2OSTiIwDSxKRmVq3bh3i4+MBALa2toiMjNTbtfjof/lM\nnjwZHh4eADQHA8+bN09yIiIqjlsAEJmhrKwseHl54caNGwA0GxcuWrRIL9dKSkpCgwYNAGjKWHp6\nOhwdHfVyLXPYAuBR//3vfzF06FAAmq0Tzp49Cy8vL8mpiAjgTBKRWVq8eLEoSLVq1cKMGTP0dq3i\nT7V16dJFbwXJXA0ZMgT+/v4AgMLCQgQHB0tOREQPsSQRmZlr165h6dKlYrxw4UJUrVpVb9fjUlvF\nqFQqxMTEiFmyLVu2aBVPIpKHy21EZub111/H+vXrAQC+vr44evQorKz08++hvLw8uLi4IDs7GwCQ\nkJAAT09PvVwLMM/ltofGjBmDVatWAQC8vLxw5swZ2NjYSE5FZNk4k0RkRvbt2ycKEgAsX75cbwUJ\nAPbv3y8KUtOmTfVakMxdRESE2KX80qVL+OSTTyQnIiKWJCIzUVRUhICAADEePnw4OnbsqNdrcpdt\n3alVqxZmz54txmFhYbhz547ERETEkkRkJr7++mucOnUKAFC5cmW9Pc1WHO9H0q0JEyaI2bj79+9r\nlSYiMjzek0RkBu7du4dmzZqJmYfw8HC9/wV7+fJlNG3aFICmlKWnp8Pe3l6v1zTne5Ie+uWXXzBg\nwAAAmnPe/vzzT7Rq1UpyKiLLxJkkIjMwb948UZA8PDwwefJkvV+z+FJbt27d9F6QLEX//v3Ru3dv\nAIBarUZgYKDZFkIiY8eSRGTiLl26pHWq/NKlS+Hg4KD363KpTT9UKhWio6NhbW0NANi7dy82btwo\nORWRZWJJIjJxkyZNQmFhIQDA398fQ4YM0fs1s7OzsWfPHjHu27ev3q9pSZo3b46PPvpIjCdPnozc\n3FyJiYgsE0sSkQnbunWrmNFRqVSIjY01yOGye/fuFX9pN2/eHI0aNdL7NS3NnDlz4OrqCkCzQWhU\nVJTkRESWhyWJyEQVFBQgKChIjN999120adPGINfmUpv+OTs7ax14u2DBAnHUDBEZBksSkYn65JNP\ncOnSJQBAtWrVEBERYZDrKoqCLVu2iDFLkv6MGTMGLVq0AKA5tDg0NFRyIiLLwi0AiEzQnTt34Onp\nifv37wMAIiMjDXYw6sWLF9G8eXMAQJUqVZCWlgZbW1uDXNsStgB41G+//Ybu3buL8eHDh9G+fXuJ\niYgsB2eSiEzQrFmzREHy9PTEhAkTDHbt4o/+9+zZ02AFyVJ169YNgwYNEuOAgACo1WqJiYgsB0sS\nkYk5deoUvvzySzGOiooyaFHh/UiGt3TpUvHf+OjRo1i7dq3kRESWgcttRCZEURR069YNe/fuBQD0\n7t0bW7duNcgTbQCQmZkJFxcXFBQUAACuX7+OevXqGeTagGUutz00ffp0LFy4EABQp04dJCQkoEqV\nKpJTEZk3ziQRmZCNGzeKgmRtbY3o6GiDFSQA2L17tyhILVu2NGhBsnShoaGoU6cOAODWrVuiMBGR\n/rAkEZmInJwcreNGPvroI3EDtaFwqU2eqlWrah1avGzZMly5ckViIiLzx5JEZCKioqJw7do1AICr\nqyvmzJlj0OsrisKSJNmIESPQrl07AEBeXh5CQkIkJyIyb7wnicgE3LhxA15eXsjKygIAfPrpp/jg\ngw8MmuHMmTNo2bIlAKB69epITU1FpUqVDJrBku9Jeig+Ph4dOnQQ499++w1du3aVmIjIfHEmicgE\nhIaGioLUokULjBkzxuAZis8i9e7d2+AFiTTat2+PESNGiHFgYKA4u4+IdIslicjIxcfH4z//+Y8Y\nx8TESCkoXGozHosWLYKDgwMA4PTp01i1apXkRETmicttREZMrVajQ4cOOHr0KABg0KBB+PHHHw2e\n4969e3Bzc0NRUREAIDk5GbVq1TJ4Di63/WP+/PmYOXMmAM09aomJiXB2dpacisi8cCaJyIh9++23\noiDZ2dlh6dKlUnLs2rVLFKS2bdtKKUikbdKkSWjYsCEAIC0tDeHh4XIDEZkhliQiI5WZmYlp06aJ\n8aRJk9C4cWMpWbjUZnwqV66sVZo/+eQTXLhwQWIiIvPDkkRkpBYuXIhbt24B0OywLOsEeLVarXVe\nG0uS8Rg8eDA6d+4MACgsLERQUJDFL0MS6RLvSSIyQleuXMGzzz6LvLw8AMA333yDt956S0qWP//8\nE76+vgA0977cvn0b1tbWUrLwnqR/O3nyJHx9fcWht7/88gv69+8vORWReeBMEpERCgkJEQWpXbt2\nWo98G1rxpbY+ffpIK0hUstatW2ttCREUFIT8/HyJiYjMB0sSkZHZs2cPNm7cKMaxsbGwspL3f1Xe\nj2T85s2bh+rVqwMAEhMTsWLFCsmJiMwDl9uIjEhhYSF8fX1x+vRpAJpjKIrvkWRoqampqFmzJhRF\ngUqlwp07d+Dq6iotD5fbHi86OhqTJk0CAFSrVg2JiYmoWbOm5FREpo0zSURGZNWqVaIgOTg4aB1o\nKsOOHTtEGWnfvr3UgkRPNn78eHh5eQEAMjIyxB5KRFR+LElERuLu3btaf7FNnz4d9erVk5iIS22m\nxNbWFtHR0WK8atUqnDhxQmIiItPHkkRkJMLDw5GWlgYAaNiwoVg6kaWoqAjbtm0TY5Yk49e3b1/x\n30lRFAQEBHBZkqgCWJKIjMCFCxfw8ccfi3FkZCQqV64sMRFw7NgxUdpq166N1q1bS81DpRMVFSXO\n9jtw4ADWr18vORGR6WJJIpJMURQEBQWJYz86d+6MQYMGSU6lvdTWt29fqU/YUel5eXlhwoQJYhwS\nEoKcnByJiYhMF3/qEUn266+/Yvv27QAAKysrxMTEaD3FJQvvRzJds2fPhpubGwAgKSkJkZGRkhMR\nmSZuAUAkUX5+Pnx8fJCYmAgAGDt2LD777DPJqYDbt2+jdu3aAABra2ukpqbCyclJcipuAVAWX3zx\nBcaOHQtAc87bpUuXUL9+fcmpiEwLZ5KIJFqxYoUoSNWrV8e8efMkJ9IofsN2p06djKIgUdm8++67\naNWqFQD4/3RMAAAgAElEQVQgJydH67BkIiodliQiSVJSUjB37lwxDgsLQ40aNSQm+geX2kyftbU1\nYmNjxXjdunU4dOiQxEREpocliUiSmTNnIiMjA4DmZtvx48dLTqRRWFgo7pECWJJMWefOnfHaa6+J\ncUBAgDgIl4iejiWJSIITJ05g1apVYhwdHQ0bGxuJif5x+PBh3L9/HwDg7u4OHx8fyYmoIpYuXQo7\nOzsAmm0d1qxZIzkRkelgSSIysEc3+evXrx/69u0rOdU/Hl1qM4Yn7aj8GjZsiJCQEDGeNm2amMEk\noidjSSIysPXr1+PAgQMAgEqVKiEqKkpyIm28H8n8TJs2TRxxc/v2bSxYsEByIiLTwC0AiAwoJycH\nzzzzDJKSkgAAkyZNwrJlyySn+seNGzfg7u4OALCxsUF6ejqqVKkiOdU/uAVA+a1duxYjRowAoDnn\n7fz582jSpInkVETGjTNJRAYUGRkpCpKbmxtmzZolOZG2rVu3io87d+5sVAWJKmb48OFo3749AM3+\nXJMnT5aciMj4sSQRGcjff/+NhQsXivH8+fONbv8hLrWZL5VKpbUlwKZNm7Br1y6JiYiMH0sSkYFM\nmzZNnKHVqlUrvPvuu5ITacvPz8fOnTvF2JhuJifdaNeuHd5++20xDgwMRGFhocRERMaNJYnIAA4e\nPIh169aJcWxsLKytrSUm+rfff/8dmZmZAIBGjRrBy8tLciLSh4ULF4pl1HPnzuHzzz+XnIjIeLEk\nEemZWq1GQECAGA8ZMgSdO3eWmKhkfPTfMtSpUwczZswQ49mzZyM9PV1iIiLjxZJEpGdr1qzB8ePH\nAQB2dnZYsmSJ5EQl4/1IliMwMBCNGjUCAKSnpyMsLExuICIjxS0AiPQoIyMDzZo1w+3btwFojiIx\nlkNsi7t69SoaN24MALC3t0daWhocHBwkp/o3bgGgOz/99BMGDRoEQHPO26lTp+Dt7S05FZFx4UwS\nkR4tWLBAFKR69eoZ7UnsxR/979q1q1EWJNKtgQMHolu3bgCAoqIiBAUFsXgSPYIliUhP/vrrL0RH\nR4vx4sWL4ejoKDHR43GpzfKoVCrExMTAykrz18DOnTuxefNmyamIjAuX24j0ZODAgfj5558BAB06\ndMDBgweN8mbo3NxcuLi4iO0JLl++LJbejA2X23Rv/Pjx+PTTTwEATZo0wblz58SBuESWjjNJRHqw\na9cuUZAAzSP/xliQAGDfvn2iIHl5eRltQSL9CA8PF5uaXr58GcuXL5eciMh4sCQR6VhhYSECAwPF\n+O2334afn5/ERE/GpTbL5ubmhvDwcDGeN28ekpOTJSYiMh4sSUQ69vnnn+PcuXMAgCpVqmgdRWKM\nipck7rJtmT744AM0b94cAPDgwQOtfZSILBnvSSLSobS0NHh6euLu3bsANLsbG+sTbQCQmJiIZs2a\nAQAcHByQnp5u1Pej8J4k/dm+fTv69OkDQPP7fPToUbRt21ZyKiK5OJNEpENhYWGiIDVu3Fhr2c0Y\nFZ9F6tGjh1EXJNKv3r17Y8CAAQA0BTQwMJBFlCweSxKRjpw7dw4rV64U42XLlsHe3l5ioqfj/UhU\n3LJly2BjYwNAc97gDz/8IDkRkVwsSUQ6oCgKgoKCUFRUBADo1q0bXnnlFcmpniwrKwt79+4VY96P\nRJ6enlrnDIaEhCA7O1tiIiK5WJKIdGDz5s3YuXMnAMDKygoxMTFG+8j/Q7/99hvy8/MBAD4+PvDw\n8JCciIzBzJkzUbNmTQDA9evXjfasQSJDYEkiqqC8vDxMmjRJjMeNG4cWLVpITFQ6XGqjklSvXh0L\nFiwQ48WLFyMpKUliIiJ5WJKIKig2NhaXL18GADg7O2Pu3LmSEz2doiha57WxJFFxo0aNwnPPPQdA\nsyP7lClTJCcikoNbABBVQHJyMpo1a4YHDx4AAJYvX44JEyZITvV058+fFye+V6tWDampqeKGXWPG\nLQAM58CBA3jxxRfFeP/+/fD395eYiMjwOJNEVAEzZswQBal58+YYN26c5ESlU3yprVevXiZRkMiw\n/P39MXToUDEOCAgQDyYQWQqWJKJyOnbsGFavXi3G0dHRJlM2uMs2lcaSJUvENhYnTpxAXFyc3EBE\nBsaSRFQOiqIgICBALPkMGDAAvXv3lpyqdDIyMnDgwAExfrjLMtGjPDw8MHXqVDGePn067t+/LzER\nkWGxJBGVw/fff49Dhw4BAGxsbLBs2TLJiUpv165dKCwsBAC0adMGdevWlZyIjNmUKVPg7u4OAEhJ\nSUFERITkRESGw5JEVEbZ2dlaT/sEBgbC09NTYqKy4aP/VBYODg5aeyXFxsYiMTFRYiIiw2FJIiqj\nJUuW4Pr16wCAmjVrYubMmZITlZ6iKCxJVGbDhg1Dp06dAAAFBQUIDg6WnIjIMLgFAFEZJCUlwcvL\nC7m5uQCAVatW4d1335WcqvROnjyJNm3aAABcXFyQkpICa2tryalKj1sAyHP8+HH4+fmJ3/dt27aZ\nzH14ROXFmSSiMpgyZYooSM899xxGjRolN1AZFZ9F6t27t0kVJJLL19cX77zzjhgHBQWhoKBAYiIi\n/WNJIiqlAwcOaJ2KHhsba3Ilg0ttVBHz589H1apVAQAXLlzAypUrJSci0i8utxGVQlFREfz8/HDi\nxAkAmns0vvvuO8mpyubu3btwc3ODWq2GSqXC7du3UaNGDdmxyoTLbfItWbJEbAvg5OSExMREuLm5\nSU5FpB+cSSIqhbi4OFGQKleujMWLF0tOVHY7duyAWq0GALRr187kChIZh4CAADRp0gQAcO/ePcyZ\nM0dyIiL9YUkieor79+9j+vTpYjxlyhR4eHhITFQ+3GWbdMHOzg5RUVFi/Nlnn+HMmTMSExHpD0sS\n0VNEREQgJSUFAODu7m6SJ6Kr1Wps3bpVjHk/ElXEgAED0LNnTwCaP1vFd58nMie8J4noCRITE+Ht\n7S2e4vnuu+8wbNgwyanK7o8//kC7du0AADVq1EBycjKsrEzv30i8J8l4nDt3Dq1atRKH3m7cuBGv\nvvqq5FREumV6PyWJDCg4OFgUpE6dOmmdim5KHl1qM8WCRMbF29sbH374oRhPnjxZbI9BZC74k5Lo\nMbZv347NmzcD0MxgxMbGas1kmBI++k/6EBYWBhcXFwDAlStXEBMTIzkRkW5xuY2oBAUFBWjVqhUu\nXLgAABg9ejS++uoryanKJyUlBbVr14aiKLCyskJqaiqcnZ1lxyoXLrcZn08++QQfffQRAMDR0REJ\nCQk8NJnMBmeSiEqwcuVKUZCqVq2K+fPnS05Uftu3bxeFomPHjiZbkMg4jR07Ft7e3gCArKwsrSdB\niUwdSxLRI1JTU7X2fpk1axZq164tMVHFcKmN9KlSpUqIjY0V42+++QZHjx6VmIhId1iSiB4xZ84c\n3Lt3DwDQtGlTTJw4UXKi8isqKsL27dvFmCWJ9KF79+4YOHCgGHNLADIXLElExZw5cwafffaZGC9b\ntgx2dnYSE1XMkSNHcPfuXQBA3bp10bJlS8mJyFxFRkbC1tYWABAfH49169ZJTkRUcSxJRP9HURQE\nBASIozt69uyJAQMGSE5VMY8++m+qT+eR8WvSpAmCgoLEeMqUKcjMzJSYiKjiWJKI/s+mTZuwZ88e\nAIC1tTWio6NNvlTwfiQypBkzZoj7927evGmSZxwSFcctAIgA5Obm4tlnn8XVq1cBABMmTMDy5csl\np6qYmzdvol69egA0N9empaWhWrVqklNVDLcAMH5xcXF45513AGjOebt48SIaNmwoNxRROXEmiQhA\nTEyMKEguLi4ICwuTG0gHtm3bJj729/c3+YJEpuGtt95C27ZtAQB5eXkmedYh0UMsSWTxbt68iYiI\nCDGeO3eu2EXYlHGpjWSwsrLS2hJg/fr12Ldvn8REROXHkkQWb/r06cjKygKgOY9q7NixkhNVXEFB\nAXbs2CHGLElkSB07dsTw4cPFOCAgQByES2RKWJLIoh09ehTffPONGMfGxqJSpUoSE+nGwYMH8eDB\nAwBAgwYN0Lx5c8mJyNIsXrwYDg4OAIBTp06Z7LE+ZNlYkshiPXzk/6GBAweie/fuEhPpzqNLbab+\nlB6ZHnd3d0ybNk2MZ8yYITZpJTIVLElksdatW4f4+HgAgK2tLSIjIyUn0p2tW7eKj7nURrJMnjwZ\nHh4eADTH/cybN09yIqKy4RYAZJEyMzPh5eWFmzdvAgCmTp2KRYsWSU6lG0lJSWjQoAEAzSPYaWlp\ncHR0lJxKN7gFgOn573//i6FDhwLQbEVx9uxZeHl5SU5FVDqcSSKLtHjxYlGQateujRkzZkhOpDvF\nZ5E6d+5sNgWJTNOQIUPg7+8PACgsLMSkSZMkJyIqPZYksjjXrl3D0qVLxXjhwoWoWrWqxES6xUf/\nyZioVCrExsaKWcBff/1V688okTHjchtZnCFDhmDDhg0AgLZt2+LIkSOwsjKPfy/k5eXBxcUF2dnZ\nAICEhAR4enpKTqU7XG4zXe+//z6+/PJLAICXlxfOnDkDGxsbyamInsw8/mYgKqV9+/aJggQAy5cv\nN5uCBAD79+8XBalp06ZmVZDItEVERIhd3y9duoRPPvlEciKipzOfvx2InqKoqEjrkf/hw4ejQ4cO\nEhPpHpfayFjVrFkTs2fPFuOwsDDcuXNHYiKip2NJIovx1Vdf4dSpUwAABwcHszyhnCWJjNmECRPE\n7Ob9+/cxa9YsyYmInoz3JJFFuHfvHjw9PZGamgpAcz6buf2A/uuvv8RfQJUrV0Z6ejrs7e0lp9It\n3pNk+rZs2YKXXnoJgOa/559//onWrVtLTkVUMs4kkUWYO3euKEgeHh6YPHmy5ES6V/zR/+7du5td\nQSLz0K9fP/Tp0weApugGBgay8JLRYkkis3fp0iWsWLFCjCMjI1G5cmWJifSDu2yTKVCpVIiKihJn\nJO7btw8bN26UnIqoZCxJZPYmTZqEwsJCAIC/vz9ee+01yYl0Lzs7G3v27BHjvn37SkxD9GTNmzfH\n+PHjxXjy5MnIycmRmIioZCxJZNaKb1z36KZ25mTv3r3Izc0FoPkLqGHDhnIDET3FnDlz4OrqCkCz\nwWtUVJTkRET/xpJEZis/P1/rCIT33nsPbdq0kZhIf/hUG5kaZ2dnREREiPGCBQtw48YNiYmI/o0l\niczWJ598gkuXLgEAqlWrpvUD2ZwoioItW7aIMUsSmYoxY8agZcuWADRLxtOmTZOciEgbtwAgs3Tn\nzh14enri/v37AIBly5aZ7cGaFy9eRPPmzQEAVapUQVpaGmxtbSWn0g9uAWB+9uzZg27duonx4cOH\n0b59e4mJiP7BmSQyS7NmzRIFqVmzZvjoo48kJ9Kf4kttPXv2NNuCROapa9euGDRokBgHBARArVZL\nTET0D5YkMjsnT57EF198IcZRUVFmXRx4PxKZuqVLl8LOzg4AcPToUXz77beSExFpsCSRWXl0c7o+\nffqYdXF48OAB9u/fL8Z89J9MUePGjREcHCzG06ZNw4MHDyQmItJgSSKz8uOPP2Lfvn0AgEqVKiEq\nKsosH/l/aPfu3SgoKAAAtGrVCvXq1ZOciKh8QkNDUadOHQDArVu3sHDhQsmJiFiSyIzk5OQgJCRE\njD/66CNxQ7O54i7bZC6qVKmideh0VFQUrly5IjEREUsSmZGoqChcu3YNAODq6orZs2fLDaRniqLw\nfiQyK2+++SbatWsHAMjLy9P6Rw+RDNwCgMzCjRs30KxZM2RnZwMAVq5ciXHjxklOpV9nzpwRe8w4\nOTnhzp074jwsc8UtAMxffHw8OnToIMa7d+/W2iKAyJA4k0RmYdq0aaIgtWzZEmPGjJGcSP+KzyL1\n6tXL7AsSWYb27dtj5MiRYhwYGCjOXiQyNJYkMnnx8fFajwzHxMTA2tpaYiLD4FIbmauFCxfC0dER\ngGbG9Msvv5SciCwVl9vIpKnVanTo0AFHjx4FAAwePBgbNmyQnEr/7t27Bzc3NxQVFQEAkpOTUatW\nLcmp9I/LbZZjwYIFmDFjBgDNPYaJiYlwdnaWnIosDWeSyKR9++23oiDZ2dlh6dKlkhMZxs6dO0VB\natu2rUUUJLIskyZNQsOGDQEAaWlpCA8PlxuILBJLEpmsBw8eaB2IGRwcjEaNGklMZDhcaiNzZ29v\nj8jISDH++OOPcf78eYmJyBKxJJHJWrhwIW7dugUAqFOnDkJDQyUnMgy1Ws39kcgiDBo0CF26dAEA\nFBUVISgoiMusZFC8J4lM0pUrV9C8eXPk5+cDANasWaP1RIw5O378ONq2bQsAcHNzQ3JyskXcqA7w\nniRLdOrUKTz33HPi0NvNmzfjpZdekpyKLAVnksgkhYSEiIL0/PPP480335ScyHCKzyL16dPHYgoS\nWaZWrVrh/fffF+NJkyaJ/+8T6RtLEpmc3377DRs3bhTj2NhYWFlZzh9l3o9Elmbu3LmoXr06ACAx\nMRErVqyQnIgsBZfbyCjFx8dj2LAw3LqlQlFRNQAqAAqsrDKgKNdQWHgRADBy5EisWbNGalZDSk1N\nRc2aNaEoCqysrJCSkgJXV1fZsQyGy22WKyYmBkFBQQCAatWqISEhgU91kt6xJJHRmTdvBWbPPgRg\nKQD3El5xHcAkVKq0E9eunbWok+/XrVsnlhY7dOiAQ4cOSU5kWCxJlqugoAAtW7bExYuafyC99957\n3GSS9M5y1ijIJGgK0m4Aa1FyQcL/ff57FBZ2RmzsfwwXzghwqY0slY2NDaKjo8X4q6++wokTJyQm\nIkvAmSQyGpqDLWOhKUil6e9qAFMxfDiwdq35byJZVFSEWrVqIS0tDYDmKbfnnntOcirD4kwSvfTS\nS9iyZQsAwN/fH/v27dP6c0GkSyxJZDQaNuyD//1vFR4/g/Q4w3Hlynyz30iy+OnotWvXxo0bNyzq\nhnWAJYmAS5cuwcfHRxx6+8MPP+D111+XnIrMlWX9hCWjduuWCmUvSAAQjh49PtR1HKNTfKmtb9++\nFleQiADAy8sLEydOFOOQkBBkZ2dLTETmjD9lyWhonmIrD0/8/bf5T7fzfiQijVmzZqFGjRoAgKSk\nJK3jS4h0iSWJjEj5i46iVNVhDuOTnJyM48ePAwCsra3Rs2dPyYmI5HFycsL8+fPFeNGiRfj7778l\nJiJzxZJERqT895ioVA90mMP4bN++XXz8wgsviI31iCzV6NGj0bp1awBATk6O1mHXRLrCkmQBzp49\ni6CgILRo0QJOTk6wt7eHh4cHBg8ejB9//LFUN8AePnwYEydORIsWLeDq6go7Ozu4u7ujV69eiImJ\nwYMHFS8p1tYZ5fzORNSvr67w9Y0Zl9qItFlbWyM2NlaM161bh4MHD0pMROaIT7eZsaKiIkyePBnL\nly9/YhHq0KED1q9fj7p16/7ra1lZWRg9ejTWr1//xGu5ublh5cqVGDx4cLnzlu/pNgXAm/jf/xbB\nw8Oj3Nc2ZoWFhXBzc8P9+/cBAKdPn0aLFi0kp5KDT7fRo15//XXx88nX1xdHjx7lQw2kM/yTZKYU\nRcHQoUMRGxsr/jKxs7ODn58f2rRpAzs7O/Haw4cPo3///sjNzdV6j7y8PPTq1UurIFlZWaFFixbw\n9/fX2uk6NTUVr7/+eoV2wP3++zAAIdDsf1QaagBT0LOnrdkWJEDz3+dhQXJ3d4ePj4/kRETGY8mS\nJbC3tweg2Tvsm2++kZyIzAlLkpmKjY3VOgT2jTfewM2bN3HkyBEcP34cV69eRf/+/cXXT506hcWL\nF2u9x5QpU3D48GEx7tmzJy5fvoxTp05h3759+Pvvv7FlyxZxfpKiKJgwYUK5d8Ft37495s7tCOA1\nPL0oqQEMBvAl9u//Hv/973/LdU1T8OhSGzfOI/pHw4YNERISIsahoaHIyCjv0j3RIxQyO1lZWYqr\nq6uiUqkUlUqldO/evcTX5eTkKM8884x4Xf369cXXLl++rFSqVEl8rWPHjkpBQUGJ73P69GnFzs5O\nvLZ///4Vyj9gwJsKMFYBlMf8T60AkxV7+7oKNOttCgBl0aJFilqtrtC1jVHLli3Fr3HTpk2y40hV\n/L830UOZmZlKvXr1xJ+NqVOnyo5EZoIzSWZo06ZNSE9PF+PZs2eX+Dp7e3vMmDEDEydOxMaNG3Hq\n1CnxtfXr16OoqEiMIyMjUalSpRLfp0WLFhg9erQY//rrr7hx40a5smdnZ+PUqQMAvgXQD0DCI69I\nBDAcw4cD58//jmeeeUZ8Zdq0aRg3bpzYidccXL9+HadPnwagObuqe/fukhMRGR9HR0etmfDo6Gj8\n9ddfEhOR2ZDd0kj33n//fTGr4+joWK7Zla5du4r3qFev3lNfv3PnTvF6lUqlrF27tjzRlfDwcPGv\nQVdXV6VBgx6KjU1fpVKl1xUbm75K48Z9lCtXrojXp6enK126dNGaYejdu7dy//79cl3f2HzxxRfi\n19WjRw/ZcaQDZ5LoMdRqtdKhQwfx5+OVV16RHYnMAGeSzND58+fFx02aNCnXPSxXr14VH7ds2fKp\nr3/0ZuJLly6V+Zp///03Fi1aJMYLFy7EtWs7kZ//KwoKfkB+/q+4fHmr1hltzs7O2LZtG0aMGCE+\nt337dvj7++P69etlzmBs+Og/UemoVCqtLQF+/vln7Ny5U2IiMgcsSWao+FJbeTcdvHPnjvjYxcXl\nqa9/9DUPT6ovi6lTpyInJwcA0Lp1a60lvCexs7PDmjVrtJYVT58+jeeffx4nT54scw5jkZ+fj127\ndokxSxLRk/n5+WHUqFFiHBQUZFbL72R4LElmqPgPBbW6fJssFn8PpRT70ZT3Og8dPHgQ3333nRjH\nxsbC2tq61N+vUqkQHh6OuLg4ce/UzZs34e/vj61bt1Yomyy///47MjMzAQCNGzdGs2bNJCciMn4L\nFixAlSpVAADnzp3D559/LjkRmTKWJDPk7OwsPn64v05ZFZ8ZSk1NferrH31NaWafHlKr1QgICBDj\nIUOG4MUXXyz19xf39ttvY/v27WIGLTMzEwMGDDDJH5R89J+o7OrUqYMZM2aI8axZs8o1s00EsCSZ\nJU9PT/Hx5cuXtZ5Se9S5c+ewbt06ZGVlaX2++KxFaZasij8ZB0DrqbOn+eabb8Thrfb29liyZEmp\nv7ck3bp1w6FDh9CgQQMAmp3Hx40bh6lTp1Z4xsuQipekvn37SkxCZFoCAwPRuHFjAMDdu3cRFhYm\nNxCZLJYkM+Tv7y8+zs3NxYEDBx772ri4OIwYMQI1atTA4MGDkZ+fDwBaMzl37tzBnj17nnjN4ps5\nqlQqre/PyclBeHg4xo8fjy5dumDYsGFISkoCAGRkZCA0NFS8NiQkBA0bNizdL/QJnn32WcTHx6Nt\n27bic0uWLMGwYcPEfU/G7OrVq7hw4QIATXHs0qWL3EBEJsTe3h7Lli0T45UrV+Ls2bMSE5HJkv14\nHene3bt3lSpVqojH8Xv27Fni6+7cuaM4OzuL17Vr10587dHNJNu0aaPk5OSU+D6HDh1SrK2txWv7\n9Omj9fXQ0FAlKSlJjEeNGqW4u7sraWlpypQpU8Qju/Xq1VMyMzN18Dvwj8zMTOXll1/WenS8Q4cO\nSkpKik6vo2uffPKJyNu3b1/ZcYwGuAUAlZJarVa6d++utYWGOW42S/rFnzRmas6cOVr7Fr311lvK\nvXv3xNcvXbqk+Pn5ab1mw4YNWu8xduxYra+/+OKLSkJCgvi6Wq1WNmzYoFW0bG1tlT/++EO8Jicn\nR3F0dFTCwsLE5xISEhSVSqVMmzZNsbW1FT/Eyru30tMUFhYqEydO1PoLtkmTJsqlS5f0cj1d6N+/\nv8i6YsUK2XGMBksSlcWZM2cUKysr8Wfm559/lh2JTAx/0pipoqIipUePHlolp3Llykr79u0Vb29v\nrVkilUqlDB48+F/vkZOT868iZWVlpfj4+CgvvPCCUqdOnX997eOPP9Z6j8zMTKVBgwZaxwTk5uYq\nKpVKadiwodbsjr7/lRcTE6OoVCpxTRcXF2X//v16vWZ5ZGdnK5UrVxY5L1++LDuS0WBJorL68MMP\ntf5xlJubKzsSmRD+pDFjeXl5yqhRo7SKzKP/s7KyUoYPH67k5eWV+B5ZWVnK8OHDn/geKpVKcXFx\nUb777rtS5Tp//rxWWQGgHD16VJe/9MfatGmTVgGxtbVV1q1bZ5Brl9bWrVtFPi8vL9lxjApLEpVV\namqq4uzsLP7cLF68WHYkMiG8cduM2draYvXq1YiPj8d7772Hpk2bwtHREfb29mjcuDHefPNN7N69\nG2vXroWtrW2J7+Hg4IC1a9fijz/+wEcffQQfHx84OzvD1tYWtWvXRteuXbFkyRJcuXIFw4YNK1Wu\nb7/9FlZW//zRGzVqFPz8/HTya36aV155Bfv27UOtWrUAaDZsHD58OBYsWFCq/aAMgbtsE+mOq6sr\nwsPDxXjevHlITk6WmIhMiUoxlr8ZyCIkJSXh2WefFVsOVKlSBQkJCahTp45Bc1y7dg39+vUTT5AB\nwLvvvouVK1fCxsbGoFke1bRpU1y+fBkAsHPnTvTo0UNqHmNSfK8o/uii0iooKEDr1q3FkU3vvPMO\nvv76a8mpyBRwJokMpqCgACNHjtT63MyZMw1ekACgYcOGOHToELp27So+99VXX6F///7l3oBTFxIT\nE0VBcnR01NrOgYjKx8bGBjExMWIcFxeHY8eOSUxEpoIliQwmKCgI1atXF7NI9evXR2BgoLQ8Tk5O\n2LZtG9566y3xuZ07d+KFF14Q+zgZWvGlth49esDOzk5KDiJz07NnTwwYMACAZhYyICCAs5H0VCxJ\nZBAxMTHw9vbWKgH+/v7SS4CtrS3i4uK07lk4e/Ys2rdvjz///NPgebjLNpH+LFu2TCynHzp0CN9/\n/73kRGTsWJJI73766Sds3rwZy5YtE0ekuLu7G81NySqVCrNnz8aaNWvED9Bbt27hxRdfxJYtWwyW\nI/ZP2G4AACAASURBVCsrC3v37hVjliQi3fL09NSavZ4yZcq/jmQiKo4lifQqLS0NI0eOxN69e8W9\nNgBw8+ZNeHl5SUz2byNHjsSOHTvg5OQEQFNaXn75ZXz66acGuf5vv/0mjoXx8fGBh4eHQa5LZElm\nzpyJmjVrAgCuX79e4bMiybyxJJFeubq6Ii0tDY0aNRKfGz9+PIqKirTOVTMWXbp0waFDh8T5cWq1\nGuPHj8fkyZP1fjguH/0n0r9q1aph4cKFYrxkyRJp9yCS8WNJIr2LjY0Vs0jOzs5a9/8Yo+bNmyM+\nPl5r76Zly5ZhyJAhyM7O1ss1FUVhSSIykFGjRuG5554DoDkEfMqUKZITkbHiPkmkV8nJyfD09ERm\nZiYAYPny5ZgwYYLkVKWTnZ2NN998E5s2bRKfe/755/H//t//E9P1unLu3Dn4+PgA0PxLNzU1Vfp+\nTcaI+ySRrvz+++9aW2zs37+fW27Qv3AmifRq+vTpoiA9++yzGDdunOREpefg4IANGzZo3eh55MgR\ntG/fHhcvXtTptYrPIvXq1YsFiUjPXnjhBa1TAgICAsSDJUQPsSSR3hw7dgyrV68W45iYGJP7y9/a\n2hrR0dFYsWKFOErl6tWr6NixI/bt26ez63CpjcjwFi9ejMqVKwMATpw4ofXzigjgchvpiaIoeOGF\nF3Do0CEAwMsvv4yff/5ZcqqK2bx5M4YNGybuS7KxscHq1avx5ptvVuh9MzIy4OrqisLCQgCaJ/9k\n7EJuCrjcRroWHh6OsLAwAEDNmjWRkJCA6tWryw1FRoMzSaQX33//vShINjY2iIyMlJyo4gYMGID9\n+/ejdu3aADTHrIwYMQIREREV+gt7165doiA999xzLEhEBhQSEoL69esDAFJSUhARESE5ERkTliTS\nuaysLK2nRQIDA+Hp6Skxke74+voiPj4e3t7e4nOzZs3C6NGjxR5HZcVdtonkcXBw0NorKTY2FgkJ\nCRITkTFhSSKdW7JkCa5fvw5AM309c+ZMyYl0q0GDBvj999/RvXt38bm4uDj069cP9+7dK9N78dF/\nIvmGDh2KTp06AdDMEAcHB0tORMaC9ySRTv3vf//DM888g9zcXADAV199hdGjR0tOpR/5+fkYO3Ys\n4uLixOe8vb2xZcsWNGjQoFTvcfLkSbRp0wYA4OLigpSUFFhbW+sjrlngPUmkL8ePH4efn5/4c7Vt\n2zb07t1bciqSjTNJpFNTp04VBcnX1xejRo2SG0iPbG1t8fXXX2PevHnic+fOnUP79u1x7NixUr1H\n8Vmk3r17syARSeLr66v1D7qgoCAUFBRITETGgCWJdObAgQP44YcfxDgmJkY8Nm+uVCoVZs6ciW+/\n/Ra2trYANBtodu7cGZs3b37q93Opjch4zJ8/H1WrVgUAXLhwAStXrpSciGTjchvpRFFREfz8/HDi\nxAkAwLBhw/Ddd99JTmVY+/btw6uvvoq7d+8CAKysrBATE4MJEyZg27ZtGDZsHjIyHKEo1QGoABQB\nuAzgFoAUpKSkoEaNGvJ+ASaAy22kb0uXLhUPnjg5OSExMRFubm6SU5EsLEmkE6tWrcKYMWMAAJUr\nV8bFixct8hT7ixcvol+/frh69ar4nI9PZ5w9WwfAUgDuJXzXdQAfYuzYlvjsMz5+/CQsSaRveXl5\n8PHxwV9//QUA+OCDD/Dpp59KTkWysCRRhd2/fx+enp64c+cOACAsLAxz5syRnEqelJQUvPzyyzhy\n5AgAJwCdAWzEk1e31QAGY/jwBli7NsYQMU0SSxIZwubNm/Hyyy8D0MwInzhxAi1btpScimRgSaIK\nCwkJEZtF1q9fHxcvXoSDg4PkVHJlZ2ejZ8+eOHTIA8BalO72PzWAqejc+Tb27l2j34AmiiWJDEFR\nFPTp0wc7duwAAHTt2hW7d+/W+vNHloEliSokISEBPj4+4imQ7777TuvQSEtWvXonZGT8gJKX2J7k\nDZw9O1Nrw0rSYEkiQzl//jxatmwpDr3duHEjXn31VcmpyNDM+9Ej0rvg4GBRkDp16oShQ4dKTmQ8\nHjxwRNkLEgDMhb//+7qOQ0Rl8Oyzz+LDDz8U4+DgYLG9CVkOliQqt23btuGXX34BoPkXfmxsLKej\ni9E8xVYenrh3z1GnWYio7MLCwuDi4gIAuHr1KqKjoyUnIkNjSaJyKSgoQFBQkBiPHj0avr6+EhMZ\no4oURp5CTv+/vfsOi+pY/wD+XToCCogVFUVBVMQSG0SxBSx4o4IF0Whi16tSLBgUFUFDsIEdEkvs\nUWxJUK/lWi9FiQUBAygmQVFRUZHe5vcHPyccWZCycJbl/TwPz7NndvbMu/vsObw7Z84MEZu+vr5g\nstg1a9YgOTlZxIhITaMkiVTKjh078McffwAAdHR0sGbNGpEjkkdVGTPzTmZREEIqb+bMmTA3NwdQ\ntHi3h4eHyBGRmkRJEqmwV69eCW7x9/T0RJMmTUSMSD5JJJVNdBKgq5su01gIIZWjoqICf/9/puX4\n6aefcPPmTREjIjWJkiRSYStWrOCr3bdr1w4LFiwQOSL5pKOTjqKJIiuCAViBsLA91RARIaQyBg8e\njFGjRvHtBQsWoLCwUMSISE2hJIlUSFRUFAIDA/n2xo0boa6uLmJE8ocxhl9//RVaWn8CcEPR/Efl\nUQhgCbp3z0T79u2rLT5CSMWtX7+er88YERGBQ4cOiRwRqQmUJJFyY4zBxcWF/4KytbXFiBEjRI5K\nvsTGxmLo0KH48ssv8ezZMwC/ARiNTydKRTNuAz/gwYMLOHXqVHWHSgipgLZt28LNzY1vu7u7Iz2d\nLosrOkqSSLmdOnUKly9fBgAoKytj06ZNdMv//0tNTcWCBQtgYWHBZ+kFAF1ddXTrlglgaRmvZgDc\noaISCuAdsrKyYG9vLxgHQQgRn4eHB5o2bQoASE5Oxvfffy9yRKS6UZJEyiU7OxsLFy7k23PnzkXH\njh1FjEg+5OfnY/v27TAxMcGWLVv47LxKSkqYM2cOEhIScPv2BfTv/xyAI4CEj/aQAGAC+vd/gZiY\n62jbti2Aol47V1dXLFiwgO+TECIuHR0dfPfdd3x73bp1+PPPP8ULiFQ/Rkg5rF27lqGoy4Pp6+uz\n169fM8YYu3//PnNxcWHm5uasQYMGTF1dnbVs2ZLZ29uz4OBgVlhYWC3xyKrd0NBQNn/+fGZubs70\n9fWZmpoaMzQ0ZDY2NmzTpk0sLS2t1NdevHiRmZub88/lw9/AgQPZvXv3StSPjo5menpWTCKxYRLJ\nGCaR2DA9PSsWHR3N66SkpDBLS0vB/v71r3+x9PT0in9ICqj450KIGAoKCliPHj3493DMmDFih0Sq\nEZ1pyCc9ffqUaWlp8ZPCtm3bWH5+PnNxcWFKSkpMIpGU+mdlZcWePn0qs1hk1W56ejobN25cmfuQ\nSCSsUaNGLDg4WPDahw8fslGjRpVIjtq0acOOHz9e5cQwMzOTjR07VrDvzz77jCUnJ1dpv4qAkiQi\nD0JDQwXfxcuXL4sdEqkmdKYhnzR58mR+MjA3N2e5ubnMwcFBkExoaGiwXr16se7duzMNDQ3Bc127\ndmVZWVlVjqOwsFAm7WZnZzMrKytBXWVlZWZhYcGsra1ZixYtBM8pKSmxoKAglpaWxpYuXcrU1NQE\nJ0gtLS22du1ambzHDwoKCtiSJUsE7bRq1UrQ61QXUZJE5MXEiRP5d7FLly4sPz9f7JBINaAzDSlT\nRESE4B/TxYsX2aZNmwRJhJOTE0tNTeWvefbsGRsxYoSgzqpVq6oci6zaXbBggaCOra0t+/PPPwV1\nzpw5w5o2bcrrqKqqsoYNG5boPZo8ebJMe8o+tnPnTqasrMzbq1+/Prt48WK1tSfvKEki8iIpKYnV\nq1ePfx8DAwPFDolUAzrTkFIVFBSw3r1785PAqFGjWEZGBmvYsCFPHgYPHiz1tVlZWczMzIzXa9my\nZZVikVW7jx49YioqKoLLcnl5eVL3FRUVVaLX6MNf7969WURERJXeU3mdOXOGaWtr87ZVVFTYnj17\naqRteUNJEpEnq1ev5t9HAwMD9ubNG7FDIjJGZxpSqv379/MTgJqaGnv48CE7ePCgoBfm6tWrZb7e\n2dmZnTx5UtDjw1jRpbMff/yRWVtbM319faasrMw0NTVZu3bt2LRp00r07MiqXV9fX8F+QkNDpe4j\nKSlJ0J3+4a9JkyZs//79rKCg4FMfn0zduXOHNW/eXBDL8uXLq21gvLyiJInIk8zMTNaqVSv+nXR1\ndRU7JCJjdKYhUr1//17wT3np0qWMMcZmzpzJEwwtLa1K/5OeNWtWmQOm69evz2JiYnh9WbU7cOBA\nvh9DQ8MSz2dmZrLVq1cLutGL/+3atatS7cpCUlISs7CwEMQzceJElp2dLVpMNY2SJCJvjh49Kujl\nffDggdghERmieZKIVL6+vkhOTgYANG3alK98HRsby+u0bdu2UpNJnjlzBkFBQXzb2NgY1tbWMDIy\n4mXv37+Hr68v35ZFuwDw+PFj/tjCwoI/Zozh2LFj6NChA1asWIHMzEzB6yQSCSQSCf76669KtSsL\nLVq0wPXr1zFkyBBedvDgQdja2iI1NVW0uAipy8aMGQNra2sARfOmFZ+Vm9R+lCSREv7880+sX7+e\nb/v6+kJHRwcABP+MGzRoUKn93717F4aGhpBIJJg/fz4ePnyIK1eu4PHjx5gyZQqvFxMTwx/Lol0A\nePnyJX+sr6/P4xkwYADGjRsnSII+zJ5dPCF7/fp1pduWhfr16+PXX3/FjBkzeNm1a9dgZWWFxMRE\nESMjpG6SSCQICAjg54mzZ8/izJkzIkdFZIWSJFLC4sWLkZOTAwDo2bMnvvrqK/5cfn4+f1zZVbA9\nPDyQlJSEN2/ewM/Pj5fn5OTwKf+Bot4kWbb78X6ysrIwa9YsdO/eHdeuXePlBgYGCAwMxO3bt9Gv\nX79Kt1VdVFVVERgYKOhpi4uLQ58+fRARESFiZITUTV27dsX06dP5tpubG3Jzc0WMiMgKJUlE4MqV\nKwgODubbAQEBUFL652uip6fHH797965KbTHGcPbsWSxZsgQDBgyAnp6eYC2k4gmNrNr90HsEAKdP\nn0ZQUBAYYwAAFRUVuLq6IiEhATNnzoSysjJevXpV6uvFJJFI4O7ujiNHjkBdXR1AUS/ZgAEDcPz4\ncZGjI6Tu8fHxQf369QEU/WjZtm2byBERWaAkiXAFBQVwdnbm2xMnToSlpaWgjomJCX/86NGjMtcV\ni4mJwaFDh5CRkSEoz8nJgYuLC5o2bQp7e3usX78e169fR4MGDWBubi51X7JoFyjqJWJFNywI9jFs\n2DDcv38fGzduhK6uLi+/d++e4PVmZmaltiuG8ePH49KlS2jYsCGAojX2xo4di40bN/LkjxBS/Ro3\nboyVK1fybS8vL6SkpIgYEZEJMUeNE/myc+dOfpdGvXr1WFJSUok6gYGBgrvQypqOf9GiRUwikTBN\nTU1mb2/PcnJyGGOMzZ07l7++U6dObN++ffyW/927d/Pn2rRpI7N2Hzx4wIYOHVribrUWLVqwkJCQ\nUvc1efJkwczbHz6TjIwMtmXLFubp6ckWLlzIRo0aJbgbr6bFx8ezdu3aCd7b3LlzS50DqrYC3d1G\n5FhOTg4zNTXl39GZM2eKHRKpIjrTEMYYY2/evGEGBgb84Pb29i61nra2Nk8cbGxspNZ7+fIl09PT\n4/V69erFGCs6iairq/Py33//XfC6NWvWSE2SKtvuZ599xlxcXJiKikqJBOlTS6aEhoYyZWVlvq+h\nQ4fy5+bPn8+GDRvGt52dnVnjxo3Zu3fvpO6rJrx8+ZJ9/vnngvdoZ2fH3r9/L1pMskZJEpF3ISEh\ngnPMnTt3xA6JVAGdaQhjjDFXV1d+YBsZGbHMzMxS665cuVLQqzN58mT29u1b/nxcXBzr2bOnoM6H\nRWKTk5MF5evXr+evu3XrFmvSpEmps2VXpl0dHZ0SiVGHDh0EdaytrVl8fDzfT2FhIQsODhYkW2pq\nauzWrVu8jouLC/vss8/49vbt25lEImE3b96sxKcvO1lZWWz8+PGC99ytW7dqXTqlJlGSRGqD4r3W\n/fv3r3OTvioSOtMQ9uDBA0FPy9GjR8usX1BQwL744gtBoqGpqcn69OnDOnXqJFj2QyKRMAcHB/7a\n/Px81qhRI8HzFhYWzMLCosSEkrq6upVu9+NeIwDM2tqa3blzh2VlZZVIppSUlJi5uTnr27cva9as\nWYnntm7dWuZn4ubmxnR0dOSi16agoIAtXbpU8N5btmzJoqKixA6tyihJIrVBbGys4Jx67NgxsUMi\nlURnGsKGDRsmSCTK86snJyeHff3112XOmq2kpMScnJz4WKQPPh5fVLy+q6ur4BLXixcvKtSutOTI\nyMiIHTt2TPC+MjIymJOTU5nxSyQSpq+vzw4fPlzmZ/Hu3TvWsmVLdujQoQp86tUvKCioxOK458+f\nFzusKqEkidQWLi4u5e6dJ/KLzjR1XFWvn0dERLAZM2YwExMTpq2tzTQ1NVnbtm3ZpEmTyhxcHRwc\nzPr06cO0tLRY/fr12eDBg9nZs2cZY4x1796dKSkpMSUlpVJX1v64XWljjurVq8e8vb3LPDlFRkay\n+fPns86dOzN9fX2mrq7OmjVrxgYNGsTWrVsnuJz3saysLObr68tGjRrF/Pz85LJL/dy5c4JLjioq\nKuzHH38UO6xKoySJ1BapqanlGudJ5BudaeqwnJwc1r59+1p7J0ZBQQHbv39/iYVfAbBJkyZJvTuv\nuri4uLD+/fuXOghcTHfv3mWGhoaCz2fZsmVymdR9CiVJpDb5+I7hJ0+eiB0SqSA609RhGzduFFyK\n+fjSljwLDw9nvXv3LpEc9ezZk4WGhtZ4PAkJCUwikTAvL68ab7s8njx5wrp27Sr4rCZMmFDrFsel\nJInUJvn5+YJFqSdNmiR2SKSCaDLJOiolJQVeXl58e+XKlWjcuLGIEZVPcnIypkyZUmIJjqZNm2Lv\n3r0IDw8vMQGmrD1//hzNmzfH6tWreVmrVq0AADdv3qzWtivL0NAQ165dw7Bhw3jZ4cOHYWNjI/p6\ndIQoKmVlZfj7+/PtAwcOICwsTMSISEVRklRHeXp68uU9TE1NMW/ePJEjKlt2djbWrl0LU1NT7Nu3\nj5erqalh6dKliI+Px5QpUwRLqFSX58+f4/nz53j79i0v+7BwrrGxcbW3X1k6Ojr45ZdfMHv2bF52\n/fp1WFlZ4dGjRyJGRojiGjhwIBwcHPi2s7NzldafJDVM7K4sUvPu3LkjuBOsrBmnxVZYWMiOHz/O\nWrduXeLS2ujRo9mjR49EiemLL74QtL1x40amp6fHEhMTazyeiiosLGR+fn6Cz9LAwECUy5QVBbrc\nRmqhxMREpq6uzr+7e/fuFTskUk4SxmiBp7qEMYaBAwfi6tWrAIrWLDtz5ozIUUkXFRUFFxcXXL58\nWVBubm4Of39/DB48WKTIgNevX8PHxwd5eXnIy8vDq1ev4OPjgw4dOogWU0UdO3YMX331FXJycgAA\n6urqOHDgAMaMGSNyZKWTSCT8MZ26SG2yfPlyrFmzBgDQrFkzxMXFQUdHR+SoyKdQklTHBAcHY+zY\nsQCKVr2/f/++3C3a+urVK3h6eiIoKEjQLa2vrw8fHx/MmDEDKioqIkaoOEJDQzFy5Ei8evWKl/n5\n+WHRokWChEReUJJEaqv09HSYmpri2bNnAIBvv/0Wa9euFTkq8ik0JqkOycrKwqJFi/j2vHnz5CpB\nysvLQ0BAAExMTLBz506eICkrK2PBggVISEjAnDlzKEGSISsrK4SHh8PU1JSXLVmyBHPnzkV+fr6I\nkRGiWLS1tfH999/z7Q0bNiAxMVHEiEh5UE9SHeLj4wNPT08AgIGBAeLj46GnpydyVEXOnTsHV1dX\n/PHHH4JyW1tbbNq0CR07dhQpsrrh9evXGDVqFG7cuMHLhg8fjiNHjsjVJQHqSSK1WWFhIaysrPid\nuaNHj8aJEydEjoqUhZKkOuLp06cwNTVFZmYmAGDnzp2YNWuWyFEB8fHxcHNzQ0hIiKC8Xbt22LRp\nE+zs7OTyso8iys7OxtSpU3H48GFe1rVrV/z2228wNDQUMbJ/UJJEaruIiAj06dOHb1+6dAmDBg0S\nMSJSFrrcVkcsXbqUJ0gWFhaYPn26qPG8e/cOixYtgrm5uSBB0tHRwbp16xAdHY0RI0ZQglSDNDQ0\ncODAASxbtoyX3b17F71790ZUVJSIkRGiOHr37o3JkyfzbRcXF7q0LceoJ6kOCAsLg5WVFd++fPky\nBgwYIEosBQUF2LNnDzw8PPjcQkBRD8HUqVOxZs0aNGnSRJTYyD927dqFWbNmoaCgAEBR8nrs2DEM\nGTJE1LioJ4koguTkZJiamiIjIwMAsH37dsyZM0fkqIg01JOk4AoLC+Hs7My3HRwcREuQrl27hh49\nemDGjBmCBKlv376IjIzEjz/+SAmSnJg2bRrOnj2L+vXrAwDev38POzs7/PDDDyJHRkjt17x5c3h4\nePBtT09PpKamihgRKQ0lSQpu//79uHXrFoCieXDWrVtX4zH89ddfGD9+PPr374+7d+/y8pYtW+LI\nkSO4du0aunfvXuNxkbLZ2Njgxo0baNmyJYCiXsCZM2fi22+/pRmDCakiNzc3tG7dGkDRjRPFl4ki\n8oMutymw9+/fw9TUFM+fPwcALFu2DD4+PjXWfkZGBvz8/ODn54fs7GxerqmpCXd3dyxevBj16tWr\nsXhI5SQnJ2PEiBG4c+cOLxs/fjz27t0LDQ2NGo2FLrcRRXLixAm+ZImysjKioqLoTl45Q0mSAvPw\n8MB3330HoKh7Ny4uDtra2tXeLmMMhw8fhru7O548eSJ4bsKECfj+++957wSpHdLT0+Ho6CgYZP/5\n55/j1KlTMDAwqLE4KEkiioQxhsGDB/NVBWxtbXHu3Dm6YUWOUJKkoBITE9GhQwfk5uYCAPbt24ev\nvvqq2tuNjIyEs7MzQkNDBeXdu3dHQEAA+vbtW+0xkOqRn58PZ2dnbN++nZe1a9cOZ86cgYmJSY3E\nQEkSUTRRUVHo1q0bv4T966+/YsSIESJHRT6gMUkKatGiRTxB6t27NyZOnFit7T1//hxTp05Fz549\nBQlS48aNsWvXLty6dYsSpFpORUUFW7duxYYNG3iy8vDhQ1haWuJ///ufyNERUjtZWFhg5syZfNvV\n1ZWfu4n4qCdJAV26dAlffPEF3w4PD0fv3r2rpa2cnBz4+/vDx8cH6enpvFxVVRUuLi5Yvnw5v0OK\nKI7jx49j0qRJfKyZuro69u3bh3HjxlVru9STRBTRy5cvYWpqirdv3wIA1q1bJ1hCioiHkiQFk5+f\nj27duiE6OhoAMHnyZPz0008yb4cxhl9++QULFy7Eo0ePBM99+eWXWL9+fY1dgiHiCA8Px5dffimY\nzsHX1xdLliyptjEVlCQRRRUQEAAXFxcARfOSJSQk0JQocoCSJAWzY8cOzJ07FwCgpaWF+Ph4NG/e\nXKZtREdHw8XFBZcuXRKUd+zYEf7+/rCxsZFpe0R+JSYmYvjw4YiLi+NlM2fOxLZt26plIWJKkoii\nysvLg4WFBV+/cvr06TQvmRygMUkKJDU1lS9gCxTd3SbLBOn169eYN28eunTpIkiQ9PT0sGXLFty7\nd48SpDrG2NgYoaGhsLa25mVBQUEYMWIE0tLSRIyMkNpFVVUVmzZt4tu7du3C7du3RYyIAJQkKRQv\nLy+8fv0aANC6dWu4ubnJZL/5+fnYunUrTExMsG3bNn4XhrKyMv79738jISEB8+bNq5aeAyL/9PX1\ncf78ecHNAf/5z3/Qr1+/ElNAEEJKN3ToUNjZ2QEo6il1dnamHlOR0eU2BREbGwsLCwu+1tbx48dh\nb29f5f1euHABrq6uiImJEZQPHjwY/v7+MDc3r3IbRDEwxrBy5Up4e3vzsubNmyMkJARdu3aVSRt0\nuY0ouvj4eHTq1IkvenvkyBGMHz9e5KjqLupJUgCMMbi6uvIEaeDAgRg9enSV9vnw4UOMHDkStra2\nggTJ2NgYp06dwoULFyhBIgISiQSrV6/G7t27ea9icnIy+vXrh7Nnz4ocHSG1g6mpqWC9zcWLFyMz\nM1PEiOo2SpIUQEhICM6fPw8AUFJSgr+/f6XvLkpLS4O7uzs6duyIX375hZdra2vD19cXsbGxGDly\nJM0IS0r1zTff4Ny5c3zqh/T0dPzrX/9CYGCgyJERUjt4enqiUaNGAICkpCSsX79e5IjqLrrcVsvl\n5uaiU6dOePjwIQBg9uzZ2LFjR4X3U1hYiL1798LDwwMvXrwQPPfNN99gzZo1aNasmUxiJnVDdHQ0\n7Ozs8Pfff/OyJUuW4LvvvoOSUuV+n9HlNlJX/PDDD3ySSU1NTcTFxdFyTiKgJKmWW79+PRYvXgwA\n0NXVRXx8PP8FUl7/+9//4OzsjN9//11QbmlpiYCAAPTs2VNm8ZK65dmzZxgxYoTgLp2xY8fip59+\ngqamZoX3R0kSqSsKCgrQo0cP3L17F0DRupeHDh0SOaq6h5KkWuzFixcwMTHB+/fvAQD+/v6Ca9mf\nkpSUBHd3dxw+fFhQbmhoCD8/P0yYMIEuq5EqS09Ph5OTE3799VdeZmlpidOnT1c4oackidQl165d\nQ//+/fn29evXaXmnGkZjkmqxZcuW8QTJzMyMTyL5KZmZmfDy8kL79u0FCZKGhgY8PT0RFxcHJycn\nSpCITGhra+PkyZOYP38+LwsLC4OlpSXi4+NFjIwQ+WZtbS1Y6sfFxYVPwUJqBvUk1VK3b99Gjx49\n+K/pc+fOYciQIWW+hjGGo0ePYvHixUhKShI8N27cOPj5+cHIyKjaYibE398fbm5u/Hurr6+P06dP\nl/vXMfUkkbrmr7/+gpmZGV8ncffu3fjmm29EjqruoJ6kWujjScbs7Ow+mSDdvn0b1tbWcHR0t8pA\nsgAAE5VJREFUFCRIXbt2xdWrV/Hzzz9TgkSqnYuLC44fP87HI6WmpmLw4MElLvkSQooYGRnxcacA\n8O2339Js9jWIkqRa6OjRo7hx4wYAQEVFBRs3biy17osXLzBjxgz06NGDvwYAGjVqhKCgIERGRgqW\nlCCkuo0ePRpXrlxB48aNARTdoenk5ITvvvuOeocIkcLd3R2GhoYAis7pa9asETmiuoMut9UymZmZ\nMDMz471BCxculDqHRm5uLjZv3ozVq1fzcUtAUVLl7OwMT09PNGjQoMbiJuRjjx8/xvDhw/mCngAw\nbdo07NixA6qqqlJfQ5fbSF116NAhvvSPqqoqYmJiYGJiInJUio+SpFrGy8sLq1atAlDUG5SQkCBI\ndhhjCAkJgZubGxISEgSvtbOzw4YNG9C+ffuaDJmQUr158wb29va4cuUKL7OxsUFwcDCfjLI4SpJI\nXcUYQ9++fREaGgoA+PLLL3H69GmRo1J8lCTVIklJSWjfvj2ysrIAFE02Nn36dP58bGwsXF1d+ezb\nH5iZmWHTpk0YOnRojcZLSHnk5uZi+vTp2L9/Py/r3LkzQkJCSkyeR0kSqcsiIyMF89adP38eNjY2\nIkak+ChJkjPh4eFwdFyFZ88kKCioD0ACgEFZOQ0qKo+RmRkHoGjAdWRkJJSVlfHmzRusWrUK27Zt\n4+u3AUCDBg3g5eWFuXPnlnr5ghB5wBiDl5cXvLy8eFmzZs3w22+/ISUlBY6O3khL0wJj/xwTEkka\n6tfPwJEjnvQDgNQZ33zzDfbu3QugMZSVjaCtbYai4cUM6uqZ6NBBB4GB7nTFQFYYkRurV29mgCMD\nkhjApPwlMWA4A1TZ1atXWV5eHtu+fTtr2LAhA8D/lJSU2OzZs1lKSorYb4mQCtm7dy9TUVEp9l1u\nWI5jwpHNmrVM7NAJqRFbtuxhgF2Zx4Sq6iS2f/9JsUNVCNSTJCe8vbdgxYpLAE6g7JsOCwEMx7hx\nrfHgQSju378veHbAgAEICAiAhYVFNUZLSPX573//C3t7e7x7VwhgIICT+PQx4QAnJyMcPOhfIzES\nIoYDB05h+vSDyMn5GZ86JtTUHLF58yjMmuVUU+EpJEqS5EB4eDgsLQMAHET5ZmUoBLAEwF4ArwEA\nrVu3xoYNGzB69GiaKZvUeoGBgZg9+woqdky4o3//F7hyZV+1xkaIGOLi4tC5sw/y8n5CeY8JZeVv\n4enZCitX/ru6w1NYlCTJgdath+Kvv34E0KKCrxwPTc1fsXz5cri5uUFDQ6M6wiOkxunqfo53735G\nxY+JCYiOXo5OnTpVR1iEiGbAgKm4enU1KnpMqKlNRlLSej4vGakYmkxSDjx7JkHF/xkAgA8aNfoc\nHh4elCARhZKWpoXKHROr0a/fTFmHQ4joHjx4j8ocE7m5nnBy8pB9QHWEitgBEPz/XWyVYYK//y6k\ny2tEAfX/dBWpTPDmjQodE0QBfVnJ15kgKuqdTCOpS6gnSS5U5YSuK7MoCFEMdEwQRVT5/xN5efVk\nGEfdQkmSXKjKsLC3MouCEMVAxwRRRJX/P6GqminDOOoWutwmB5SV01BsDsgKSICxsSoePaKx90Sx\nKCnZonK3lCRATy8Pqal0TBDF0qTJWKSkVOaVCejShdbprCzqSZIDzZoxAE8q+CoGYCUuXw6qhogI\nEVf9+hmo3DGxAmFhe6ohIkLE1aGDDipzTKipeePwYd/qCKlOoCRJDhw5sgrAYhTN9VIeRfMkOTkZ\nolWrVtUWFyFiOXLEE5U5Jvr3V6XlGIhCCgx0h6rqt6jIMaGsvBQeHr1hYGBQnaEpNEqS5ECfPn2w\nerUVgDH49AFQCGAMFi9uiIMH11V/cISIYOjQoZg1qy0AB5TvmHCAk1MeTSRJFFb79u2xe7cD1NQc\nUZ5jQk3NEdu2daGJJKuIkiQ54ek5H4sX9wLwbRm1GAB3LF7cC35+S2soMkLEsXOnD5ycjACU9V0v\nOiZoSRJSF0yaNAqbN4+CsnJZ8x4xKCt/+8klSaKjo+Hq6orOnTtDV1cXGhoaaNWqFRwcHHD8+HHQ\nPNNFaMbtcggLC8Phw4dx+fJlJCcnIz09HY0aNULHjh0xfPhwTJs2DTo6OjJpa+LExTh06AmA1QBM\nij2TAGAFnJxaUA8SqVMGDJiMq1dzAXij5DHhif791agHidQpXl7bsHZtOHJzV+DjY0JNbTU8PPqU\n2oNUUFCARYsWYfPmzWUmQpaWljh27BiaN28u2+BrGUqSypCRkYGpU6fi2LFjZdYzMDDAjh074ODg\nIJN2Hz9+jC++mIukJAkY04FE8h4tWzJcvLgdbdq0kUkbhNQmMTEx6NdvJt6+1QLQAMA76Opm4Pr1\nIFqChNRJKSkpcHLyQFTUO+Tl1YOqaiYsLBrg0KG1pS5BwhjD2LFjceLECV6mrq4OCwsL5OfnIzY2\nFjk5Ofy5Ll26ICwsrE6v6EBJUilycnIwaNAghIWF8TIlJSV06tQJurq6SExMxNOnT/lzEokEO3fu\nxIwZM8QIlxBCCCmTv78/3Nzc+PaECROwdetW6OnpAQCeP3+OGTNmICQkhNdZuXIlVq5cWeOxygtK\nkkrh7OyMLVu28G0bGxsEBQXByMiIl509exZTp07FixcvAABqamoICwtDt27dajxeQgghpDSZmZlo\n1aoVUlNTAQCDBg3CxYsXS9TLzs5Gt27dEBcXBwBo0aIF/v777xqNVZ7QwG0pEhMTsX37dr5taWmJ\nkJAQQYIEAMOGDcP58+ehpqYGAMjNzYWnp2eNxkpIbSDGINHCwkIEBwfDyckJ7dq1g46ODrS0tGBs\nbAxHR0ccOXKEBqeSOuPUqVM8QQKAFStWSK2noaGBZcuWYcGCBThx4gTu3bsneP727dtwcnJC69at\noa6uDlVVVTRt2hRDhgzB6dOnq/U9iIKREnx9fZlEIuF/oaGhZdafM2eOoP6TJ09qKFJC5Ft+fj5z\ncXFhSkpKgmPk4z8rKyv29OlTmbX78OFD1qVLlzLblEgkzNzcnMXExMisXULk1cyZM/n3XktLixUW\nFlZ4H5cvX2bq6uplHlPe3t7VEL14qCdJiv/85z/8cfPmzWFpaVlmfXt7e8H21atXqyUuQmoTxhjG\njx+PgIAA3mOjrq6Onj17olu3blBXV+d1w8LCYGdnh+zs7Cq3m5iYCEtLS0RFRfEyTU1N9OrVC1ZW\nVmjQ4J8lGmJiYmBlZYXbt29XuV1C5FlsbCx/3LZtW0gkFVswNz8/HxMnTkRubi4AoH79+ujbty+6\nd+8OZWVlXm/NmjXIysqSTdBygJIkKR4/fswfW1hYfLK+ubm5YPvDtVxC6rKAgADBXTQTJkxAcnIy\nIiIi8Pvvv+Px48ews7Pjz9+7dw/ff/99ldosLCzEmDFj8OrVK17m6uqKlJQUhIeH48aNG0hJScHG\njRuholK0dGVaWhrGjh2L9PT0KrVNiDwrfqmt+A+F8oqOjoa+vj7U1NTQpk0bPHr0CNeuXUNkZCQu\nXLjA6+Xm5iIhIUEmMcsDSpKkePnyJX+sr6//yfof13n9+rXMYyKkNsnMzISPjw/fHjRoEA4ePMjv\nogGApk2b4tixY4JlRHbt2lWldg8cOIC7d+/y7YULF2LDhg3Q0tLiZaqqqnBxcRHcmPH48WPs3Lmz\nSm0TIs/y8/P548LC8i5t8o+uXbvi/v37yMjIQGRkJBo2bMifK/4/kDGmUD84KEmSoviXiZVjYGdl\nvnCEKDKxBokePnyYP9bV1YW3t3epMc6aNQsdOnTg20FBtFg0UVzFf6C8e/eu0vtRUlJCXFwc1q1b\nhzFjxsDQ0BBdu3YV1Cn+P7S2UxE7AHmkr6+P58+fA4Cg2740H9cpT+8TIYqs+Li8evXqoV+/fqXW\nnTRpEiZNmlSi/MqVKxg6dCgfA/FBSkoKLly4gAsXLmD16tVYvnw5gKK5za5cucLrDRs27JOT4I0e\nPRoPHjwAADx69AjJycl1foZhophMTExw8+ZNAEXf9YKCAsFYouJiYmJw7949jBw5UtALe/r0aTg7\nOwumBNDQ0EDfvn1x48aN6n0DIqGeJClMTU354+Jd96X5+NevmZmZzGMipDYRY5Do8+fPS8wW/CnF\nxxMyxvDHH39UKE5CaoviP1Sys7Nx/fr1Uuvu3bsXkyZNQqNGjeDg4IDc3FyEhYXBwcEBf//9N9TV\n1bFq1SpEREQgLS1NoW9WoiRJCmtra/745cuXuHz5cpn1jx49yh9LJBL++qysLHh5eeHf//43BgwY\nAEdHxzo9KRepO2pykOjDhw8BCMcSApUbT1g8bkIUybhx4wS9QmvXrpVa79WrV3xsYHZ2Np48eQI1\nNTUEBQXxoSUzZszAihUr0LNnT6ioqODZs2fV/wZEQkmSFF9//bXg1+rChQtLvTU5LCwMBw8e5Nu2\ntrZo0aIFAMDb2xtTp07Ftm3bcOXKFWhqauLzzz+nEzFReDU5SPT9+/cl2vzw3KfQeEJSV+jq6mLh\nwoV8++LFi5gyZYpgfFJ8fDyGDx+Ot2/f8rIlS5YAAJKTk3lZREQE/5/47t07LFiwQNCWIo1JoiRJ\nCmNjY0ybNo1v3717F0OGDBHc1sgYw/Hjx2FnZ8dPtKqqqnygaHZ2NjZv3ozdu3fz13h4eODp06fY\ns2dPDb0TQsQhxiDR4m0C5bvLlMYTkrpkxYoVGDx4MN/ev38/mjVrBktLS5ibm6NTp06IjIzkz9vb\n2/OF24svJH3r1i20bt0aVlZWaNmypWCqDwD8h4sioCSpFP7+/ujRowffvn79OszMzNC5c2f069cP\nhoaGGDt2LM+4JRIJNm7cyF9TUFAAAwMDwaRarVq1AlA0aI4QRWZiYsIffxgkWpqYmBgcOnQIGRkZ\ngvLTp0+jTZs2sLKygru7O06cOIE3b96gb9++UvdjZGTElwgCUK4JIouPJ5RIJDSekCg0JSUlhISE\nYMqUKbwsOzsbERERiI2N5cepRCLBhAkTcOjQIV5v4cKFaNSoEd/+MPdYeno6unbtikGDBvHnFGqu\nQPEm+5Z/GRkZzMnJ6ZNLG+jr67PDhw9/cn+xsbFMIpGwdevW1UD0hIgnMDBQcIxcvny51LqLFi1i\nEomEaWpqMnt7e5aTk8NCQ0OZsrIyk0gkTENDg3l5ebGbN2+yvLw8VlhYKNj31atX+b769evHy7W1\ntVlaWlqp7RYUFLA2bdrw+m3btpXlR0CIXIuIiGAzZsxgJiYmTFtbm2lqarK2bduySZMmlXq8JiYm\nsgkTJrBGjRoxDQ0NZmZmxnx8fFh2djbbvHkzk0gkTElJiVlaWtbsm6lGlCSVQ2RkJJs/fz7r3Lkz\n09fXZ+rq6qxZs2Zs0KBBbN26dezt27fl2s/y5ctZkyZNWGpqajVHTIi43rx5w7S1tXkCYmNjI7Xe\ny5cvmZ6eHq/Xq1cvxhhjX3/9NS+bP3++4DVPnz4tNUnau3ev4LnZs2eXGqOfn5+grq+vrwzeOSFE\nkVCSVEP++usv1rBhQxYSEiJ2KITUiJUrVwqSkMmTJwt+UMTFxbGePXsK6gQHBzPGGLO1tRUkTllZ\nWYwxxt6+fcscHBwEr7l06RLfZ15eHuvcuXOJRKl4u1lZWczb21tQx8jIqMxeJ0JI3SRhrBy3gJAq\nycvLw7BhwzBnzhw+CI4QRVdYWIghQ4bg0qVLvExDQwNdunTB+/fvERcXJxirZG9vj+DgYACAm5sb\n/P39+XONGzeGsbExoqOjSyx5cPLkSYwcOZJvx8XFoW/fvoKB2+rq6rCwsICysnKJfWhra+PixYvo\n1auX7N48IUQhUJJUA+bNm4cRI0Zg6NChAIoGsrZt21bkqAipfrm5uZg1axZ++umnUutIJBI4Ojpi\nz549fOD106dP0b179xJzHwFF0wPo6+vjv//9LwDA19eX36b8QWJiIsaNG/fJwdvt27fHzz//XK6F\nrAkhdQ/d3VbN/P39YWdnxxOkvLw8HDlyROSoCKkZampq2LNnD8LDwzF9+nS0a9cOWlpa0NDQgLGx\nMSZOnIhLly7h4MGDgjvTDA0NER4eDkdHRxgYGEBdXR3t27eHt7c3wsLCeM+RRCLBqVOnSrRrbGyM\nyMhInDhxAk5OTjA2Noa2tjY0NTVhZGQEe3t77Nu3DzExMZQgEUJKRT1J1ejkyZPYunUrbG1t+cR2\nT58+haWlJZycnESOjhBCCCFloSSpmrx+/RpGRkbIysoSzPwrkUgQEREhmIOJEEIIIfKHkiRCCCGE\nECloTBIhhBBCiBSUJBFCCCGESEFJEiGEEEKIFJQkEUIIIYRIQUkSIYQQQogUlCQRQgghhEhBSRIh\nhBBCiBSUJBFCCCGESEFJEiGEEEKIFJQkEUIIIYRIQUkSIYQQQogUlCQRQgghhEhBSRIhhBBCiBSU\nJBFCCCGESEFJEiGEEEKIFJQkEUIIIYRIQUkSIYQQQogUlCQRQgghhEhBSRIhhBBCiBSUJBFCCCGE\nSEFJEiGEEEKIFJQkEUIIIYRIQUkSIYQQQogUlCQRQgghhEhBSRIhhBBCiBSUJBFCCCGESEFJEiGE\nEEKIFJQkEUIIIYRIQUkSIYQQQogUlCQRQgghhEhBSRIhhBBCiBSUJBFCCCGESEFJEiGEEEKIFJQk\nEUIIIYRIQUkSIYQQQogUlCQRQgghhEhBSRIhhBBCiBSUJBFCCCGESEFJEiGEEEKIFJQkEUIIIYRI\n8X9QAKDOHfiRhQAAAABJRU5ErkJggg==\n"
      }
     ],
     "prompt_number": 18
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from pymatgen.phasediagram.pdanalyzer import PDAnalyzer\n",
      "analyzer = PDAnalyzer(pd)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 19
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Get Energies\n",
      "------------"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "for e in entries:\n",
      "    print(\"%s: %f\" % (e.name, analyzer.get_e_above_hull(e)))"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "CaO: 0.000000\n",
        "C: 0.832340\n",
        "C: 0.006175\n",
        "C: 0.944071\n",
        "C: 0.528945\n",
        "C: 0.000000\n",
        "C: 0.003043\n",
        "C: 0.506401\n",
        "Ca: 0.000000\n",
        "CO: 0.692016\n",
        "CO2: 0.028775\n",
        "CaC2: 0.015081\n",
        "CaC2: 0.014315\n",
        "C: 0.028801\n",
        "C: 0.022568\n",
        "CO2: 0.000000\n",
        "CaCO3: 0.312190\n",
        "CaO: 0.065085\n",
        "C: 0.145186\n",
        "C: 0.004413\n",
        "CaC2: 0.022313\n",
        "C: 0.174406\n",
        "C: 0.029940\n",
        "C: 0.762539\n",
        "C: 0.135631\n",
        "Ca: 0.013330\n",
        "O2: 0.000000\n",
        "CaCO3: 0.008149\n",
        "Ca(CO2)2: 0.046302\n",
        "Ca: 0.013583\n",
        "CaCO3: 0.007359\n",
        "CO2: 0.025395\n",
        "CaCO3: 0.024281\n",
        "CaCO3: 0.013584\n",
        "C: 0.160479\n",
        "C: 0.002697\n",
        "C3O2: 0.513353\n",
        "CaCO3: 0.000000\n",
        "CaCO3: 0.002015\n",
        "Ca: 0.000508\n",
        "Ca: 0.396161\n",
        "Ca: 0.017145\n",
        "CO: 0.668938\n",
        "CaCO3: 0.169799\n",
        "CaC2: 0.762125\n"
       ]
      }
     ],
     "prompt_number": 20
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 20
    }
   ],
   "metadata": {}
  }
 ]
}