diff --git a/arc/reaction/reaction.py b/arc/reaction/reaction.py
index 39f6752bd3..74fe9744f2 100644
--- a/arc/reaction/reaction.py
+++ b/arc/reaction/reaction.py
@@ -550,13 +550,16 @@ def determine_family(self,
                          ):
         """
         Determine the RMG reaction family.
-        Populates the .family, and .family_own_reverse attributes.
 
         Args:
             rmg_family_set (str, optional): The RMG family set to use.
             consider_rmg_families (bool, optional): Whether to consider RMG's families in addition to ARC's.
             consider_arc_families (bool, optional): Whether to consider ARC's families in addition to RMG's.
             discover_own_reverse_rxns_in_reverse (bool, optional): Whether to discover own reverse reactions in reverse.
+
+        Returns:
+            tuple[str | None, bool | None]: The reaction family label,
+                and whether the family's template also represents its own reverse.
         """
         if rmg_family_set == 'default' and consider_rmg_families and consider_arc_families and not discover_own_reverse_rxns_in_reverse:
             # these are the default values, don't bother generating a new product_dicts list, use the property
diff --git a/ipython/Tools/External symmetry and optical isomers.ipynb b/ipython/Tools/External symmetry and optical isomers.ipynb
index f5da0b3752..2fdca13a23 100644
--- a/ipython/Tools/External symmetry and optical isomers.ipynb	
+++ b/ipython/Tools/External symmetry and optical isomers.ipynb	
@@ -19,28 +19,14 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stderr",
-     "output_type": "stream",
-     "text": [
-      "RDKit WARNING: [13:53:02] Enabling RDKit 2019.09.3 jupyter extensions\n"
-     ]
-    }
-   ],
-   "source": [
-    "from arc.molecule import Molecule\n",
-    "from arc.species import ARCSpecies\n",
-    "from IPython.display import display\n",
-    "from arc.plotter import show_sticks\n",
-    "from arc.species.converter import molecules_from_xyz, check_xyz_dict"
-   ]
+   "outputs": [],
+   "source": "from arc.molecule import Molecule\nfrom arc.species import ARCSpecies\nfrom IPython.display import display\nfrom arc.plotter import show_sticks\nfrom arc.species.converter import check_xyz_dict\nfrom arc.species.perceive import perceive_molecule_from_xyz"
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -90,100 +76,10 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stderr",
-     "output_type": "stream",
-     "text": [
-      "WARNING:pint.util:Could not resolve planks_constant: UndefinedUnitError()\n",
-      "WARNING:pint.util:Could not resolve plank_constant: UndefinedUnitError()\n"
-     ]
-    },
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "SMILES:  [O]OC(c1cc(ccc1OC)C(C)C)NC(C)C\n",
-      "External symmetry:  1\n",
-      "Optical isomers:  2\n"
-     ]
-    },
-    {
-     "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAJkAAABpCAIAAACmrBb0AAAABmJLR0QA/wD/AP+gvaeTAAAPeklEQVR4nO2de1AT1x7HDw4oD6UEQhJe8h6oSOAqWiw+GqQVKmop0ooSn3Pj6HTiUO3EuVwbR4aWWtuJM9YRW0ujIGlsqaUU8AHcFkaFKiBalFRBoCUkoCRgEkgg5/6xECMghs1m8/wMfyQnu+f3C989m9/Z/f3O2kEIgQ2LYJaxHbCBGeatZVRUlFQqNbYXpoKd+Z5j1Wq1k5OTTCazt7c3ti8mgRmPS7FYTCAQbEJqMGMthUKhl5eXsb0wIWxaWg42LS0H89bS29vb2F6YEOatpW1camPT0nKwaWk52LS0HMz1ug+E0NnZub+/39HR0di+mArmOi77+/udnZ1tQmpjrlraTrCTMSctlUplamrqmTNngJaW+fn5KSkpSqXS2N4ZH7PRUqVSvf/++0NDQxkZGQAALy+vffv2AQAyMjLs7Ozeeeed4eFhY/toZMwj9lGpVO+9955SqSwuLp4zZ87kTxGZf/rpp8mfWg9mMC4nCllUBKKjgaMjIJMBgwGePHFwcPj+++8dHR1TUlKseXSagZYFBQWzZs0aG3N5eWDvXnDwIBCJwP/+B7q6QEICGB5G5HRxcSkoKDC2v8YDmhEKBSQQ4Jkzz1pkMkihwLw84/lkQpjBuHzGzZugvx+kpT1rcXYG69eDK1eM55MJYVZa9vUBFxcwb95zjV5eoK/PSA6ZFmalJZEIZDIwOPhco1AIiEQjOWRamJWWixcDNzdw4cKzFrkclJSA1auN55MJYVZJbE5OICcH7N8PXFxAYiIQCkFmJiCRwPbtxvbMNDB28DVzCgoglQpnz4ZEIty1C/b1Ic39/f2GsKZSqQoLC728vNhstlqtNoQJrDCP6z66sHv3boVCcfLkyblz52LS4eDg4DfffMPhcIKCguh0+tdff00gEPLz88lkMib9Y4+xDybMkMvlTCYzMDDw2rVrenbV09PDZrOJRGJycvL169eRRpVKxWazKRTKL7/8orezBsFytET48ccfyWQym80eGRlBsXtzczODwSASiQwGo7W1FWlUq9UdHR3I6+vXrwcHBzMYDJlMhpnTGGFpWkIIOzs733jjDRqN1tXVpfteNTU1ycnJFAqFzWY/fvwYaRwaGvr2228jIiLWrFmj2VIqlWZkZCxYsKCpqQlj1/XDArWEEKrVag6H4+npWVRUNP2Wo6OjJSUlS5cuDQ0N5XA4crkcaZdKpRwOx9fXNy4urqSkZHLUw+fzKRQKh8MxnYDIMrVE+OOPP0JDQ+l0+uDg4ORPBwYGOBzO/PnzJ6jV3t7OYrE8PT3pdHpzc/M0/be3t8fFxb311lvd3d0G+QIzxJK1hOMBUVhY2K1btzSNQqGQzWZ7eHhohzYQwoaGBjqd7u7uzmQyOzs7dekfCYi8vLxKS0ux936GWLiWCD/88AOJRGKz2Y2NjUhow2QyHz16hHyqVquvXLmSnJyMTCJRzFNNJCCyCi0hhB0dHb6+vo6Ojrm5uRq1hoeHuVxuREREdHQ0l8tVKpWo+5dIJJs3b46IiLh9+zZGLs8Ya9ESQshisY4cOYK8lkgkHA7Hx8fnRaENOrhcrqenp7ECIrO6tq4fIpHIz88PAFBeXh4UFHT79u1Lly7V1tauW7fOzs4OExNbt26tr6/n8/lJSUk9PT2Y9DkD8D98jMWaNWsqKioghFKpVCQSGc6QUqk8ePCgj49Pfn6+4axMxoq0pFKpeP6YnThxwt/fHzdz0KrOsTinulMolJiYGI1pHCxai5YqlWpgYMDDwwM3i5pDp7GxMTU1FQeL1qKlSCTy9PScNQu/76vRsquri4hLFou1aIl/LZHGIm6mbVoa3KJlaikQCHg8Hp4WNdi0xJLW1tb4+PjR0VHcLGpj0xIzBALBm2++mZOTs2XLFnwsTgBnLUdGRiQSCRLyWJSWAoEgISEhOzt727ZtOJibEpy11A6bcVtTakzLvXv33r9/3xAGuru74+Pjjxw5sm3btoULgasrePx47CMeD4xPpsHChUC7QovHA+HhWLqBs5Yac2q1uq+vj0Qi4WB0TMuYmJiVK1eePn0acwNz5849d+7c9vF0ZAcH8NlnmBt5OcbSsq+v75VXXnFwcMDB6JiWO3furKur++6771JSUh5rBg4WuLq60mg0zdvMTHDqFMD5FgKegwMB/8AHaP9eBgYG/v7771FRUZGRkeXl5QayFx4OUlLAJ58YqPupwXNwIBhZSwCAvb394cOHCwoKdu/evW/fPgMtzsFmg/x80NU1sX3PHkAkjv0xGFha7O7utvgJCZgyjo2Pj29ubhaLxa+//nprayvmJoOCwObN4MiRie2ffAKamsb+cnOxtGgNk0vwojmJm5tbUVERk8lcvnz58ePH9TEAIbx69WpLS4t243//C4qKwMOHz21JIABf37E/d3d9bE7EqrVE2Lp1a01NzdmzZ999910UAZFSqTx79mxkZORHH300oV7Hzw/s3Am+/HLG7qIAQlhaWiqTyfAwNo7JaQkACA8Pv3HjBpVKpVKpFRUVOnY6MDBw/PjxoKCg06dPf/rppw0NDfPnz5+wzX/+A4aG0HisO8jBFB0dXVdXV1FRcezYMbVabViTAAAA1Gp1b28vUg6G6ylBx/yDq1ev+vr6MpnM4eHhaTZra2tjMpkEAoFOp9+5c0f/vAd0IDnpfn5+mjS7zs7OVatWxcfH//3334a2jlz0QV4HBAS0tbUZ2iLCDPJ9ent7169fHxMToymA0ubWrVuapO8ZFeVoUKvVMplMIpGoVCoUuyNo56TX1dVpfzQ6OsrhcEgkEo/HQ92/LjQ1NVGpVOR1Zmbm9Ec/hsw4d4vL5RKJRA6Hg7zVJH0HBARopxGjpqKiwtfXF0XRXVNTE51OJxAITCZTU2KHsGvXrt9++w15XV9fHxISQqfTnz59qqerU9Ld3f3222/PmzdPLBYbov9pQJOH9+eff0ZFRaWkpBw9ejQ8PDwmJobH4+kzmCbQ09OTlJQUGxv78OFDXbbXLrd78uTJ5A0uXbrk7e2dlZWFODkwMMBgMMLDwxsaGrDyGUJ49+7dHTt2uLu779mzZ/v27T4+PpcvX8aw/5eCMqdyaGgoKiqKTCZXV1dj6s8YmqK7goKCF22DVBBERkZSqdS8vDyFQjFNh2KxeN26dUuWLBEIBEjLhQsXkCKT0dFRPb1FDiakhrdvfP2E6upqPz8/JpM5NDSkZ/86gj4/dv/+/ceOHcPQlcncuXMnMjIyLS1twql7cmijY4dcLpdMJueNr7nW0dGxcuXK1atX//PPPyjcQ2o3X3vttZCQEO3aTQ0SiSQ9PT0iImL62j/9QQ4g9Fpu3ry5sLAQO3+mBim6CwgIqK2thRB2d3ez2WxPT8+0tLT6+noUHba0tERHR6empiLFzyMjI7m5ud7e3j///LPunQwODnI4HH9//7i4OD6fP/1PuybCwLzIZGRk5MKFC7GxsXFxcVAfLWk0WlVVFXaOTUdxcTGJRAoODvbw8Dhw4AC6OFmDQqFgsVj+/v6agKiurk7HgGjKZSleSltb27JlyzZs2NDb26uP5xqGhoa4XG5YWNjixYu5XC5yMKHXMjw8/N69e5h4pgvnzp3z9vaWSqVYdXj58mUfHx8mk4mU6kmlUgaD8eqrrzY2Nk65/ZTLUuiOpuq2rKxMH7fFYjFyZkpOTr5y5Yr2R+i1dHV1lUgk+rg1I86fP5+eno5tn2KxODk5ecmSJX/99RfSgqxCkJubqx0QaYc2mmUp0FFVVYU6IHrw4AGTyXRzc6PT6S0tLZM3QKmlTCZzdnZGty86vvjiiw8//BDzbtVqdV5ennZA9OjRoxUrViQkJHR1dU25LIWe9Pf3b9q0aeHChboHRDU1NWlpaRQKhcViTbM2AkotHzx4EBQUhG5fdBw4cODzzz83UOfIjHnjxo3I9FSlUn388cezZ8+m0Wi//vqrIQpjJ1xymRIkTl62bFlwcDCHw3lpAT1KLWtqapDYCTe2bNkyzVxTf+Ry+QcffBAYGHjjxg2kZcLFI8y5f//+okWLUlJSNFNSDQqFYnJo81JQasnn8zdu3IhuX3TQaLTKykpDWyktLT158qShrWhQKpVsNtvb27u8vBxpEYlEmtAGmYbpDso1R414d9egrF271tAmtHFwcDh8+PDy5ct37NhBo9Hs7e0vXryYnp5+7dq1kJCQmfZm09L4JCQkNDY2BgUFJSYmCgQC1AV+KPPWcf7PKhQKpVLp5uaGm0WcIRKJFAolJydHn0pN89BSKBRSKBTczBkF/f+lZqOlpZ5gEZ4+fQoA0HMVY5uWJgEmXxCNliqVanBw0CjrOFgqRtOyp6eHRCLhsI5DYWHh4OAg0PqqxcXFYrHY0HbxB5PMejR64DNKsrKyjh49iqxuprEoEAhoNJpIJDK0dZzB5F+KZn6Jg5aHDh0qLS2trKxEwoH8/Hyk/eDBg3Z2dqtWraqurraks67Fanno0KGSkpLKysopJ1ssFgsAQKPRLElOoVAYGRmpZydozrEikchws73s7OznhJz05FIAAIvF2rRpU0JCwsDAgIHcwBlshgeKK8Lt7e3okp104eLFi8/uG5w6Bd3cYFERlEhgSwtMTIT/+hccv4vL4/GmXEjdHFmwYMHdu3f17MSE16m0pieXEggEPTMWoEmvU2k1Ty4dHh5WKBQEAkHPfkxYS6t5cqlQKCSTyfqvLW1wLauqqmpra7Vbbt68WVZW9vI9rebJpVjNCwyrZVlZWXp6upOTk3aji4vLrl27+Hz+lLvI5fKxV1bz5FLM5niY/HRPSXl5OZlMnjIhuKWlxcfHZ0L+DvLAnrVr1z5r+uor6OYGeTwokcB792BiIqRS4bR1I+bIiRMn9u7dq38/htJyGiERtOXUfmDPxJz0Fzy51JLIysrKzs7Wvx+DaFlTU0MikTQJbfD8eRgVBefMgSQS/Pe/4Xjw3dzc7O7uvnTpUtQP7LEMdu7ceUZ76oUWg2j58OHDZ4m8L57vDw8PR0RE5Obm6vPAHgsgKSlJz8IEBANfK3jZfN90HjhoRGJiYl5UxDIjsH9eNJ/Pl8vlY4sZ1taCFSvAwMBz08Tdu8GTJ88FqNZNa2traGio/veDsZ+TMJnMPXv2jB0iVjPf14ewsDBMbuyjzI+dhrKyMqVSOXYVQzPf15bTEuf7pgD243LRokWxsbFjb6xmvm8KYD8un8PJCeTkgP37gYsLSEwEQiHIzAQkEhhfGtgGhmAf+0xBYSE4ehTcvw9cXcGGDeCzzwCOOXzWAy5a2sAFE77nZWOG2LS0HP4PTxvULK95WTgAAAAASUVORK5CYII=\n",
-      "text/plain": [
-       "Molecule(smiles=\"[O]OC(c1cc(ccc1OC)C(C)C)NC(C)C\")"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    },
-    {
-     "data": {
-      "application/3dmoljs_load.v0": "<div id=\"3dmolviewer_15931954791471632\"  style=\"position: relative; width: 400px; height: 400px\">\n        <p id=\"3dmolwarning_15931954791471632\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n        </div>\n<script>\n\nvar loadScriptAsync = function(uri){\n  return new Promise((resolve, reject) => {\n    var tag = document.createElement('script');\n    tag.src = uri;\n    tag.async = true;\n    tag.onload = () => {\n      resolve();\n    };\n  var firstScriptTag = document.getElementsByTagName('script')[0];\n  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n});\n};\n\nif(typeof $3Dmolpromise === 'undefined') {\n$3Dmolpromise = null;\n  $3Dmolpromise = loadScriptAsync('https://3dmol.csb.pitt.edu/build/3Dmol.js');\n}\n\nvar viewer_15931954791471632 = null;\nvar warn = document.getElementById(\"3dmolwarning_15931954791471632\");\nif(warn) {\n    warn.parentNode.removeChild(warn);\n}\n$3Dmolpromise.then(function() {\nviewer_15931954791471632 = $3Dmol.createViewer($(\"#3dmolviewer_15931954791471632\"),{backgroundColor:\"white\"});\n\tviewer_15931954791471632.addModel(\"\\n     RDKit          3D\\n\\n 40 40  0  0  0  0  0  0  0  0999 V2000\\n   -1.2234    2.3879   -0.3134 O   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.1480   -0.1944    2.3919 O   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.2458   -0.4645    3.0320 O   0  0  0  0  0  1  0  0  0  0  0  0\\n   -1.7224   -1.4210    0.4870 N   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.4775   -0.8093   -0.0681 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.8681   -1.5928   -0.9629 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.3923   -0.0964    0.9409 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.8991   -1.0677    1.3827 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.2789   -2.1263   -0.8273 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.2076   -3.0540   -1.2315 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.8957   -0.6518   -1.6013 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.2112    3.6664   -0.9107 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    2.2407    0.0642   -0.1407 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.1181    0.4201    0.3389 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.0655    1.6772   -0.2774 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.0265   -0.3608    0.3959 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    2.2687    1.3137   -0.7474 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.1360    2.1193   -0.8224 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    4.2890   -0.2602   -0.5594 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.8888   -1.3820   -1.4044 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.2001    0.6353    0.8529 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    4.8237   -1.6513    1.4178 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    4.0633   -0.1294    1.9182 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.1290   -1.6289    1.9207 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.0030   -1.9430   -1.8689 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    4.1963   -2.7220   -0.8140 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    2.4858   -2.7261   -0.3701 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -3.1825   -3.3114   -0.8023 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.2564   -3.2473   -2.3058 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4556   -3.7105   -0.7894 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.6229    0.3977   -1.4624 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -3.8865   -0.8091   -1.1605 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.9707   -0.8388   -2.6762 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.2250    4.0514   -0.8151 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5186    4.3405   -0.3952 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.9445    3.6090   -1.9718 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.5758   -1.7118    0.9554 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.9441   -1.3333    0.8707 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.2057    3.0841   -1.3071 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.1963    1.6777   -1.1788 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1 12  1  0\\n  1 15  1  0\\n  2  7  1  0\\n  2  3  1  0\\n  4  6  1  0\\n  4  7  1  0\\n  4 37  1  0\\n  5  8  1  0\\n  5  9  1  0\\n  5 13  1  0\\n  5 19  1  0\\n  6 10  1  0\\n  6 11  1  0\\n  6 20  1  0\\n  7 14  1  0\\n  7 21  1  0\\n  8 22  1  0\\n  8 23  1  0\\n  8 24  1  0\\n  9 25  1  0\\n  9 26  1  0\\n  9 27  1  0\\n 10 28  1  0\\n 10 29  1  0\\n 10 30  1  0\\n 11 31  1  0\\n 11 32  1  0\\n 11 33  1  0\\n 12 34  1  0\\n 12 35  1  0\\n 12 36  1  0\\n 13 17  2  0\\n 13 16  1  0\\n 14 16  2  0\\n 14 15  1  0\\n 15 18  2  0\\n 16 38  1  0\\n 17 18  1  0\\n 17 40  1  0\\n 18 39  1  0\\nM  RAD  1   3   2\\nM  END\\n\",\"sdf\");\n\tviewer_15931954791471632.setStyle({\"stick\": {}});\n\tviewer_15931954791471632.zoomTo();\nviewer_15931954791471632.render();\n});\n</script>",
-      "text/html": [
-       "<div id=\"3dmolviewer_15931954791471632\"  style=\"position: relative; width: 400px; height: 400px\">\n",
-       "        <p id=\"3dmolwarning_15931954791471632\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n",
-       "        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n",
-       "        </div>\n",
-       "<script>\n",
-       "\n",
-       "var loadScriptAsync = function(uri){\n",
-       "  return new Promise((resolve, reject) => {\n",
-       "    var tag = document.createElement('script');\n",
-       "    tag.src = uri;\n",
-       "    tag.async = true;\n",
-       "    tag.onload = () => {\n",
-       "      resolve();\n",
-       "    };\n",
-       "  var firstScriptTag = document.getElementsByTagName('script')[0];\n",
-       "  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n",
-       "});\n",
-       "};\n",
-       "\n",
-       "if(typeof $3Dmolpromise === 'undefined') {\n",
-       "$3Dmolpromise = null;\n",
-       "  $3Dmolpromise = loadScriptAsync('https://3dmol.csb.pitt.edu/build/3Dmol.js');\n",
-       "}\n",
-       "\n",
-       "var viewer_15931954791471632 = null;\n",
-       "var warn = document.getElementById(\"3dmolwarning_15931954791471632\");\n",
-       "if(warn) {\n",
-       "    warn.parentNode.removeChild(warn);\n",
-       "}\n",
-       "$3Dmolpromise.then(function() {\n",
-       "viewer_15931954791471632 = $3Dmol.createViewer($(\"#3dmolviewer_15931954791471632\"),{backgroundColor:\"white\"});\n",
-       "\tviewer_15931954791471632.addModel(\"\\n     RDKit          3D\\n\\n 40 40  0  0  0  0  0  0  0  0999 V2000\\n   -1.2234    2.3879   -0.3134 O   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.1480   -0.1944    2.3919 O   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.2458   -0.4645    3.0320 O   0  0  0  0  0  1  0  0  0  0  0  0\\n   -1.7224   -1.4210    0.4870 N   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.4775   -0.8093   -0.0681 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.8681   -1.5928   -0.9629 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.3923   -0.0964    0.9409 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.8991   -1.0677    1.3827 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.2789   -2.1263   -0.8273 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.2076   -3.0540   -1.2315 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.8957   -0.6518   -1.6013 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.2112    3.6664   -0.9107 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    2.2407    0.0642   -0.1407 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.1181    0.4201    0.3389 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.0655    1.6772   -0.2774 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.0265   -0.3608    0.3959 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    2.2687    1.3137   -0.7474 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.1360    2.1193   -0.8224 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    4.2890   -0.2602   -0.5594 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.8888   -1.3820   -1.4044 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.2001    0.6353    0.8529 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    4.8237   -1.6513    1.4178 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    4.0633   -0.1294    1.9182 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.1290   -1.6289    1.9207 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.0030   -1.9430   -1.8689 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    4.1963   -2.7220   -0.8140 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    2.4858   -2.7261   -0.3701 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -3.1825   -3.3114   -0.8023 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.2564   -3.2473   -2.3058 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4556   -3.7105   -0.7894 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.6229    0.3977   -1.4624 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -3.8865   -0.8091   -1.1605 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.9707   -0.8388   -2.6762 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.2250    4.0514   -0.8151 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5186    4.3405   -0.3952 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.9445    3.6090   -1.9718 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.5758   -1.7118    0.9554 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.9441   -1.3333    0.8707 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.2057    3.0841   -1.3071 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.1963    1.6777   -1.1788 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1 12  1  0\\n  1 15  1  0\\n  2  7  1  0\\n  2  3  1  0\\n  4  6  1  0\\n  4  7  1  0\\n  4 37  1  0\\n  5  8  1  0\\n  5  9  1  0\\n  5 13  1  0\\n  5 19  1  0\\n  6 10  1  0\\n  6 11  1  0\\n  6 20  1  0\\n  7 14  1  0\\n  7 21  1  0\\n  8 22  1  0\\n  8 23  1  0\\n  8 24  1  0\\n  9 25  1  0\\n  9 26  1  0\\n  9 27  1  0\\n 10 28  1  0\\n 10 29  1  0\\n 10 30  1  0\\n 11 31  1  0\\n 11 32  1  0\\n 11 33  1  0\\n 12 34  1  0\\n 12 35  1  0\\n 12 36  1  0\\n 13 17  2  0\\n 13 16  1  0\\n 14 16  2  0\\n 14 15  1  0\\n 15 18  2  0\\n 16 38  1  0\\n 17 18  1  0\\n 17 40  1  0\\n 18 39  1  0\\nM  RAD  1   3   2\\nM  END\\n\",\"sdf\");\n",
-       "\tviewer_15931954791471632.setStyle({\"stick\": {}});\n",
-       "\tviewer_15931954791471632.zoomTo();\n",
-       "viewer_15931954791471632.render();\n",
-       "});\n",
-       "</script>"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    }
-   ],
-   "source": [
-    "\n",
-    "if smiles is None:\n",
-    "    mol = molecules_from_xyz(check_xyz_dict(xyz))[1]\n",
-    "    smiles = mol.to_smiles()\n",
-    "\n",
-    "rmg_mol = Molecule(smiles=smiles)\n",
-    "arc_spc = ARCSpecies(label='test', smiles=smiles, xyz=xyz)\n",
-    "arc_spc.determine_symmetry()\n",
-    "\n",
-    "print('SMILES: ', smiles)\n",
-    "print('External symmetry: ', arc_spc.external_symmetry)\n",
-    "print('Optical isomers: ', arc_spc.optical_isomers)\n",
-    "\n",
-    "display(rmg_mol)\n",
-    "success = show_sticks(xyz=xyz)"
-   ]
+   "outputs": [],
+   "source": "\nif smiles is None:\n    mol = perceive_molecule_from_xyz(check_xyz_dict(xyz))\n    smiles = mol.to_smiles()\n\nrmg_mol = Molecule(smiles=smiles)\narc_spc = ARCSpecies(label='test', smiles=smiles, xyz=xyz)\narc_spc.determine_symmetry()\n\nprint('SMILES: ', smiles)\nprint('External symmetry: ', arc_spc.external_symmetry)\nprint('Optical isomers: ', arc_spc.optical_isomers)\n\ndisplay(rmg_mol)\nsuccess = show_sticks(xyz=xyz)"
   },
   {
    "cell_type": "code",
diff --git a/ipython/Tools/Perceive_xyz_(xyz_to_smiles).ipynb b/ipython/Tools/Perceive_xyz_(xyz_to_smiles).ipynb
index 5eda953f57..dc73f256cd 100644
--- a/ipython/Tools/Perceive_xyz_(xyz_to_smiles).ipynb
+++ b/ipython/Tools/Perceive_xyz_(xyz_to_smiles).ipynb
@@ -20,7 +20,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -40,40 +40,12 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": null,
    "metadata": {
     "pycharm": {}
    },
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "SMILES: O=CN(C)C\n"
-     ]
-    },
-    {
-     "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAEQAAAAmCAIAAADyTaq0AAAABmJLR0QA/wD/AP+gvaeTAAAEIklEQVRYhe2YX0haXxzAzww2agku7ZKgmMNClJoRPWQkYT30j8EeLsSgXmz25H0LX0bCRsNHhRYzGCsYc1REXHqLwLInGQZRSzB7mXIp7sqcpmn6/T3cfmXarHVvTIafp+v3fM/33M/heu459wEAoH8F3t8aeHFx0e/3c1vzXCYej3NbtzDJZJKmaa1WOzo6GolEuCp7LqPRaIaHh/f397mq+zsymcz09LRcLpdKpT6fL5lM1tXV2e32dDrNQXUAAIBoNGqxWIRCocViicfjcD94PB6tVqtSqXAcp2maCXq9Xp1O19TUtLq6yrI+yv7h9/txHFcoFLOzsyzr5hAKhYxGo0QicTgc4XDYZDJhGGa321OpFJNAkmRtbW1/f//e3t6dR0H5oZWVlYaGBr1ev7m5eW2fSCSC43goFLrNACcnJ1arVSgUEgQRiUQu4j6fr6+vr76+fmlpKTsTwzCz2ZydCQBTU1MTExPZEbfbbTKZ0un0DTIAkEqlHA6HWCw2Go0HBwc5Jm1tbUajMZPJ3GhCkqRcLi8w38vLyyqVqqura3t7m4kEg8HBwUGJRDIzM8MMMTk5KZPJAoFAdsdYLKbX64eGhrJ9rpdhODw8JAiiurraZrMxz0M0Gu3o6DAYDDlTks/GxoZOp9NoNC6Xq3BmMpm02WwikYggiHA4zAQ9Hk9ra2tLS8vY2JhMJtvd3c3vGIvFOjs7BwYGLp7VQjIMOzs7PT09SqVyYWHhNiY0TRMEgWGYzWY7Ozu7sf7vemUymZGREalUeq3JtT43yzCQJFlZWalQKBgTtRr4fPh/QQKnE5qbAQDW1taYOT46Orpl5Wy8Xm97ezuzsrlcritP15cv8OwZPHoEGAavXsHPnxc+er1+fHz8D2QAgKKo+fl55lqthqoqGB3NlXG73VtbW3fQyIZZ2Xp7ey+Wb/jwAQQCcDohHIbv36G7G5qaIJFgGk9PT4+Pj/9MJhu1Gt6+BT4fKOqKDFcwK1swGAQAiMfhyRP4+PGyORaDmhpwOHJ63X1vplSiFy/Qu3ccvLjzKS8vN5vNNTU1CCH07Rs6OkI4ftlcUYGeP0fLyzm9WG00LRb06RP68YNNjUKUlZUhhBBNo8ePEZ9/pU0sRjSdk89K5ulT9PIlevOGTY1bIBKhWAz9+nUlSFFIJMpJZHsEeP0aOZ0oEGBZpiDNzUggQHNzl5GTE0SSqLMzN/Nuf1C1Gubmzq9NJqiq4ngByOX9exAI4OtXCIdhZwe6u6GxEfI2xBzIUBRUVNyzDAB8/gyNjfDwIYhEYDBcvuOyeAClYzMDRVHr6+tc3Qp7WMl4vV6r1crVrbDnr33QuA9KMsVKSaZYKckUKyWZYqUkU6yUZIoVVjI8Ho/HK6LpYHU4SyQSwWBQoVBweENs+KdOmv8Brc8HwMxCO/cAAAAASUVORK5CYII=\n",
-      "text/plain": [
-       "Molecule(smiles=\"O=CN(C)C\")"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    }
-   ],
-   "source": [
-    "from arc.species.converter import molecules_from_xyz, check_xyz_dict\n",
-    "from IPython.display import display\n",
-    "\n",
-    "\n",
-    "s_mol, b_mol = molecules_from_xyz(check_xyz_dict(coords))\n",
-    "\n",
-    "print(f'SMILES: {b_mol.to_smiles()}')\n",
-    "\n",
-    "display(b_mol)"
-   ]
+   "outputs": [],
+   "source": "from arc.species.converter import check_xyz_dict\nfrom arc.species.perceive import perceive_molecule_from_xyz\nfrom IPython.display import display\n\n\nmol = perceive_molecule_from_xyz(check_xyz_dict(coords))\n\nprint(f'SMILES: {mol.to_smiles()}')\n\ndisplay(mol)"
   },
   {
    "cell_type": "code",
diff --git a/ipython/Tools/Plot xyz.ipynb b/ipython/Tools/Plot xyz.ipynb
index ec8ff3e30a..52d25795f3 100644
--- a/ipython/Tools/Plot xyz.ipynb	
+++ b/ipython/Tools/Plot xyz.ipynb	
@@ -17,7 +17,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -32,239 +32,48 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 8,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
-   "source": [
-    "from arc.plotter import show_sticks, plot_3d_mol_as_scatter\n",
-    "from arc.species.converter import check_xyz_dict, molecules_from_xyz\n",
-    "from arc.species.converter import xyz_to_pybel_mol, pybel_to_inchi\n",
-    "from IPython.display import display\n",
-    "%matplotlib inline\n",
-    "\n",
-    "xyz = '/home/alongd/Code/ARC/ipython/Tools/ts1016.log'\n",
-    "\n",
-    "xyz = check_xyz_dict(xyz)"
-   ]
+   "source": "from arc.plotter import show_sticks, plot_3d_mol_as_scatter\nfrom arc.species.converter import check_xyz_dict\nfrom arc.species.converter import xyz_to_pybel_mol, pybel_to_inchi\nfrom arc.species.perceive import perceive_molecule_from_xyz\nfrom IPython.display import display\n%matplotlib inline\n\nxyz = 'ts1016.log'\n\nxyz = check_xyz_dict(xyz)"
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
-   "source": [
-    "path = '/home/alongd/Code/runs/T3/35/iteration_0/ARC/calcs/Species/C2H4_0/freq_a6950/output.out'"
-   ]
+   "source": "path = '/path/to/ESS/output/file/output.out'"
   },
   {
    "cell_type": "code",
-   "execution_count": 9,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "data": {
-      "application/3dmoljs_load.v0": "<div id=\"3dmolviewer_15931965379759881\"  style=\"position: relative; width: 400px; height: 400px\">\n        <p id=\"3dmolwarning_15931965379759881\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n        </div>\n<script>\n\nvar loadScriptAsync = function(uri){\n  return new Promise((resolve, reject) => {\n    var tag = document.createElement('script');\n    tag.src = uri;\n    tag.async = true;\n    tag.onload = () => {\n      resolve();\n    };\n  var firstScriptTag = document.getElementsByTagName('script')[0];\n  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n});\n};\n\nif(typeof $3Dmolpromise === 'undefined') {\n$3Dmolpromise = null;\n  $3Dmolpromise = loadScriptAsync('https://3dmol.csb.pitt.edu/build/3Dmol.js');\n}\n\nvar viewer_15931965379759881 = null;\nvar warn = document.getElementById(\"3dmolwarning_15931965379759881\");\nif(warn) {\n    warn.parentNode.removeChild(warn);\n}\n$3Dmolpromise.then(function() {\nviewer_15931965379759881 = $3Dmol.createViewer($(\"#3dmolviewer_15931965379759881\"),{backgroundColor:\"white\"});\n\tviewer_15931965379759881.addModel(\"\\n     RDKit          3D\\n\\n 22 21  0  0  0  0  0  0  0  0999 V2000\\n    2.6783   -0.0344   -0.2127 N   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.2400   -0.1946    0.1514 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3809    0.2797   -1.0317 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.0395   -1.2749    0.3399 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.9535    0.3566    1.0732 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.1181    0.0557   -0.8257 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.5613    1.3553   -1.2292 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.7051   -0.2600   -1.9496 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.6615    0.7865    0.3640 N   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.7024    0.2992   -1.7419 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.9662   -0.0688    1.4732 C   0  0  0  0  0  3  0  0  0  0  0  0\\n   -2.7781    1.7066    0.0055 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.0886   -0.5618    1.9005 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.6565    0.3225    2.2149 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.5558   -1.1768    0.9836 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -3.0537    2.3105    0.8860 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -3.6835    1.1912   -0.3555 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.4396    2.4098   -0.7729 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.2197   -1.4020   -0.5946 O   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.4537   -1.7208   -0.2269 O   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.2705   -0.4938    0.4719 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    2.9396    0.9463   -0.2247 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1  2  1  0\\n  1 21  1  0\\n  1 22  1  0\\n  2  3  1  0\\n  2  4  1  0\\n  2  5  1  0\\n  3  6  1  0\\n  3  7  1  0\\n  3  8  1  0\\n  6  9  1  0\\n  6 19  1  0\\n  6 10  1  0\\n  9 11  1  0\\n  9 12  1  0\\n 11 13  1  0\\n 11 14  1  0\\n 12 16  1  0\\n 12 17  1  0\\n 12 18  1  0\\n 15 20  1  0\\n 19 20  1  0\\nM  RAD  1  11   2\\nM  END\\n\",\"sdf\");\n\tviewer_15931965379759881.setStyle({\"stick\": {}});\n\tviewer_15931965379759881.zoomTo();\nviewer_15931965379759881.render();\n});\n</script>",
-      "text/html": [
-       "<div id=\"3dmolviewer_15931965379759881\"  style=\"position: relative; width: 400px; height: 400px\">\n",
-       "        <p id=\"3dmolwarning_15931965379759881\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n",
-       "        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n",
-       "        </div>\n",
-       "<script>\n",
-       "\n",
-       "var loadScriptAsync = function(uri){\n",
-       "  return new Promise((resolve, reject) => {\n",
-       "    var tag = document.createElement('script');\n",
-       "    tag.src = uri;\n",
-       "    tag.async = true;\n",
-       "    tag.onload = () => {\n",
-       "      resolve();\n",
-       "    };\n",
-       "  var firstScriptTag = document.getElementsByTagName('script')[0];\n",
-       "  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n",
-       "});\n",
-       "};\n",
-       "\n",
-       "if(typeof $3Dmolpromise === 'undefined') {\n",
-       "$3Dmolpromise = null;\n",
-       "  $3Dmolpromise = loadScriptAsync('https://3dmol.csb.pitt.edu/build/3Dmol.js');\n",
-       "}\n",
-       "\n",
-       "var viewer_15931965379759881 = null;\n",
-       "var warn = document.getElementById(\"3dmolwarning_15931965379759881\");\n",
-       "if(warn) {\n",
-       "    warn.parentNode.removeChild(warn);\n",
-       "}\n",
-       "$3Dmolpromise.then(function() {\n",
-       "viewer_15931965379759881 = $3Dmol.createViewer($(\"#3dmolviewer_15931965379759881\"),{backgroundColor:\"white\"});\n",
-       "\tviewer_15931965379759881.addModel(\"\\n     RDKit          3D\\n\\n 22 21  0  0  0  0  0  0  0  0999 V2000\\n    2.6783   -0.0344   -0.2127 N   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.2400   -0.1946    0.1514 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3809    0.2797   -1.0317 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.0395   -1.2749    0.3399 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.9535    0.3566    1.0732 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.1181    0.0557   -0.8257 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.5613    1.3553   -1.2292 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.7051   -0.2600   -1.9496 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.6615    0.7865    0.3640 N   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.7024    0.2992   -1.7419 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.9662   -0.0688    1.4732 C   0  0  0  0  0  3  0  0  0  0  0  0\\n   -2.7781    1.7066    0.0055 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.0886   -0.5618    1.9005 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.6565    0.3225    2.2149 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.5558   -1.1768    0.9836 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -3.0537    2.3105    0.8860 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -3.6835    1.1912   -0.3555 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.4396    2.4098   -0.7729 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.2197   -1.4020   -0.5946 O   0  0  0  0  0  0  0  0  0  0  0  0\\n   -2.4537   -1.7208   -0.2269 O   0  0  0  0  0  0  0  0  0  0  0  0\\n    3.2705   -0.4938    0.4719 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    2.9396    0.9463   -0.2247 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1  2  1  0\\n  1 21  1  0\\n  1 22  1  0\\n  2  3  1  0\\n  2  4  1  0\\n  2  5  1  0\\n  3  6  1  0\\n  3  7  1  0\\n  3  8  1  0\\n  6  9  1  0\\n  6 19  1  0\\n  6 10  1  0\\n  9 11  1  0\\n  9 12  1  0\\n 11 13  1  0\\n 11 14  1  0\\n 12 16  1  0\\n 12 17  1  0\\n 12 18  1  0\\n 15 20  1  0\\n 19 20  1  0\\nM  RAD  1  11   2\\nM  END\\n\",\"sdf\");\n",
-       "\tviewer_15931965379759881.setStyle({\"stick\": {}});\n",
-       "\tviewer_15931965379759881.zoomTo();\n",
-       "viewer_15931965379759881.render();\n",
-       "});\n",
-       "</script>"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    },
-    {
-     "data": {
-      "text/plain": [
-       "True"
-      ]
-     },
-     "execution_count": 9,
-     "metadata": {},
-     "output_type": "execute_result"
-    }
-   ],
+   "outputs": [],
    "source": [
     "show_sticks(xyz)"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "[O]N=C\n",
-      "\n",
-      "\n",
-      "InChI=1/CH2NO/c1-2-3/h1H2\n",
-      "\n",
-      "\n",
-      "multiplicity 2\n",
-      "1 O u1 p2 c0 {2,S}\n",
-      "2 N u0 p1 c0 {1,S} {3,D}\n",
-      "3 C u0 p0 c0 {2,D} {4,S} {5,S}\n",
-      "4 H u0 p0 c0 {3,S}\n",
-      "5 H u0 p0 c0 {3,S}\n",
-      "\n",
-      "\n",
-      "\n"
-     ]
-    },
-    {
-     "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAEYAAAAcCAIAAABTWV5CAAAABmJLR0QA/wD/AP+gvaeTAAAEdUlEQVRYhd1YbUhabRh+lLY4lWQlx9MXlAUuhY0QighqZIFpjEYMYiCMja0gqGSBPxoUg8XZn2URlD9GRDGREy36+BFWZFQYTXlRi+hHH64yY5o4dTPK5/1xSu0TKY+8vNev+76e59zXc/Gcc+7nHBqEEPy/QI+ynsVimZqaolQiqpa2t7dLS0t3d3eplYHRws7ODofD6e7uplooSpai5gdGx5LVas3Kyurp6YEQ8vmQwYC/fp0NqVRQIIiwXDSeJTqd3tnZ2dDQQKYPHoDPn6mUo7D2OVAUra6uDqQyGejrAwcHVMlF+yUOAHj0CDx/Djo6qKofPUt+vz8Qt7WB/n7w8yclQlGyND4+bjabAymHA16+BB8/UqIVQ0nVEBgMhubmZrfbrdVqQ/kPH0BeHsjKirwihbtkt9ubmpoqKytrampWVlYYDEboaGYmeP0afPkSZKRSqdFovKlaV1eXTqcLZQiCGBkZuWZqhJsChBDC4+NjhULBYrEaGxudTidJbm1tQQj5fEgQZ9OsVhgXF+xLAwMDbDZbKpUeHh5eKtje3s7j8Ww2WyhpMpkwDFOr1ZcmR96SRqPh8Xjl5eWrq6sks7e3J5VKMzIyLBbL7dc6HA65XM5ms3Ec9/l8JInjOJfL3d/fvzp/bW0tPT19aGgolIywJbfbXVZWNjExQaZerxfHcRRF5XK5y+UKs8j6+rpYLOZyuZOTk7f4IXHV1QVLfD5/cHAwkKpUKi6Xe1Ot4eHhwsJCBEFYLJZYLF5aWiL509NTMhgbG8vOzq6qqtrc3AzTTCg0Gk1qampOTs7BwcEZ9e0bfPIExsZCFIVv30K7naSNRiOGYfPz8/eypFQqGQyGUqnc3d09OjoiCKKuri4wajAYSkpK8vPztVrtHcwEsLm5aT9fN+zrg0wmVKmg0wnX1qBIBPPz4d+/5ODOzo7D4bi7Ja/Xy2QylUplKOn3+wMxiqIKheLk5OQ+fi7gzx+YlAS/fg0yHg/EMHhxDSTu0pf0er3T6aytrQ0laTRaIN7Y2EhMTLxD5Rvx4wc4OgIvXgSZuDjw7BnQaMC7dyTR2tqamZlZX19/eZcSEhJSzsFgMK7dpdHR0fj4+IjtQDj4/h1eVWxrg0+fkqHX642Njc3Ly4NXPy46Ojr+OQeO4yS5sLBQVFRUXFwsk8kAACkpKR6Px+VyRXIfbgeLBTwe8Pv3BdJqBSwWGSIIYjabZ2dnwdXTQ1JSUsY5kpOTSZLD4czNzS0uLlosFpPJJBAImEymWq0OvRBS+qdJIABMJiCIIOP1grExIBQGiNzcXAzDQJhnvLS0NDKIiYmh0+kIguA43tLSQqPRJBIJgiAzMzPT09O9vb2RtBEKBAGfPoH370F8PBCJgNUKZDKAouDVq2smh96ct7/xlpeXRSJRICUIoqCggOxLEokk0JcoxNAQfPwYPnwIWSz45k3wc/8iwj092Gy24uLiYNf7DyMsSz6fr6KiQq/XU72aiCAsS/39/RiGCYVCoVCo0+moXtM98S9BaguW+MSkZwAAAABJRU5ErkJggg==\n",
-      "text/plain": [
-       "Molecule(smiles=\"[O]N=C\")"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    }
-   ],
-   "source": [
-    "# Analyze graph\n",
-    "\n",
-    "s_mol, b_mol = molecules_from_xyz(xyz)\n",
-    "mol = b_mol or s_mol\n",
-    "print(mol.to_smiles())\n",
-    "print('\\n')\n",
-    "print(pybel_to_inchi(xyz_to_pybel_mol(xyz)))\n",
-    "print('\\n')\n",
-    "print(mol.to_adjacency_list())\n",
-    "print('\\n')\n",
-    "display(mol)"
-   ]
+   "outputs": [],
+   "source": "# Analyze graph\n\nmol = perceive_molecule_from_xyz(xyz)\nprint(mol.to_smiles())\nprint('\\n')\nprint(pybel_to_inchi(xyz_to_pybel_mol(xyz)))\nprint('\\n')\nprint(mol.to_adjacency_list())\nprint('\\n')\ndisplay(mol)"
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "data": {
-      "application/3dmoljs_load.v0": "<div id=\"3dmolviewer_1593194001914734\"  style=\"position: relative; width: 400px; height: 400px\">\n        <p id=\"3dmolwarning_1593194001914734\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n        </div>\n<script>\n\nvar loadScriptAsync = function(uri){\n  return new Promise((resolve, reject) => {\n    var tag = document.createElement('script');\n    tag.src = uri;\n    tag.async = true;\n    tag.onload = () => {\n      resolve();\n    };\n  var firstScriptTag = document.getElementsByTagName('script')[0];\n  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n});\n};\n\nif(typeof $3Dmolpromise === 'undefined') {\n$3Dmolpromise = null;\n  $3Dmolpromise = loadScriptAsync('https://3dmol.csb.pitt.edu/build/3Dmol.js');\n}\n\nvar viewer_1593194001914734 = null;\nvar warn = document.getElementById(\"3dmolwarning_1593194001914734\");\nif(warn) {\n    warn.parentNode.removeChild(warn);\n}\n$3Dmolpromise.then(function() {\nviewer_1593194001914734 = $3Dmol.createViewer($(\"#3dmolviewer_1593194001914734\"),{backgroundColor:\"white\"});\n\tviewer_1593194001914734.addModel(\"\\n     RDKit          3D\\n\\n  5  4  0  0  0  0  0  0  0  0999 V2000\\n    1.1746   -0.1531    0.0000 O   0  0  0  0  0  1  0  0  0  0  0  0\\n    0.0630    0.3515    0.0000 N   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.1271   -0.1133    0.0000 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.9380    0.6017    0.0000 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.2977   -1.1874    0.0000 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1  2  1  0\\n  2  3  2  3\\n  3  4  1  0\\n  3  5  1  0\\nM  RAD  1   1   2\\nM  END\\n\",\"sdf\");\n\tviewer_1593194001914734.setStyle({\"stick\": {}});\n\tviewer_1593194001914734.zoomTo();\nviewer_1593194001914734.render();\n});\n</script>",
-      "text/html": [
-       "<div id=\"3dmolviewer_1593194001914734\"  style=\"position: relative; width: 400px; height: 400px\">\n",
-       "        <p id=\"3dmolwarning_1593194001914734\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n",
-       "        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n",
-       "        </div>\n",
-       "<script>\n",
-       "\n",
-       "var loadScriptAsync = function(uri){\n",
-       "  return new Promise((resolve, reject) => {\n",
-       "    var tag = document.createElement('script');\n",
-       "    tag.src = uri;\n",
-       "    tag.async = true;\n",
-       "    tag.onload = () => {\n",
-       "      resolve();\n",
-       "    };\n",
-       "  var firstScriptTag = document.getElementsByTagName('script')[0];\n",
-       "  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n",
-       "});\n",
-       "};\n",
-       "\n",
-       "if(typeof $3Dmolpromise === 'undefined') {\n",
-       "$3Dmolpromise = null;\n",
-       "  $3Dmolpromise = loadScriptAsync('https://3dmol.csb.pitt.edu/build/3Dmol.js');\n",
-       "}\n",
-       "\n",
-       "var viewer_1593194001914734 = null;\n",
-       "var warn = document.getElementById(\"3dmolwarning_1593194001914734\");\n",
-       "if(warn) {\n",
-       "    warn.parentNode.removeChild(warn);\n",
-       "}\n",
-       "$3Dmolpromise.then(function() {\n",
-       "viewer_1593194001914734 = $3Dmol.createViewer($(\"#3dmolviewer_1593194001914734\"),{backgroundColor:\"white\"});\n",
-       "\tviewer_1593194001914734.addModel(\"\\n     RDKit          3D\\n\\n  5  4  0  0  0  0  0  0  0  0999 V2000\\n    1.1746   -0.1531    0.0000 O   0  0  0  0  0  1  0  0  0  0  0  0\\n    0.0630    0.3515    0.0000 N   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.1271   -0.1133    0.0000 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.9380    0.6017    0.0000 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.2977   -1.1874    0.0000 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1  2  1  0\\n  2  3  2  3\\n  3  4  1  0\\n  3  5  1  0\\nM  RAD  1   1   2\\nM  END\\n\",\"sdf\");\n",
-       "\tviewer_1593194001914734.setStyle({\"stick\": {}});\n",
-       "\tviewer_1593194001914734.zoomTo();\n",
-       "viewer_1593194001914734.render();\n",
-       "});\n",
-       "</script>"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    },
-    {
-     "data": {
-      "text/plain": [
-       "True"
-      ]
-     },
-     "execution_count": 4,
-     "metadata": {},
-     "output_type": "execute_result"
-    }
-   ],
+   "outputs": [],
    "source": [
     "show_sticks(xyz)"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stderr",
-     "output_type": "stream",
-     "text": [
-      "WARNING:pint.util:Could not resolve planks_constant: UndefinedUnitError()\n",
-      "WARNING:pint.util:Could not resolve plank_constant: UndefinedUnitError()\n"
-     ]
-    },
-    {
-     "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAb4AAAEuCAYAAADx63eqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXSU5b0H8O/7zpKZyQRCzEJiFJBFhbJchQK1t/UUFVoFq5bi0V6LPbW1LlwpUq0919N6W1wQrUu9Xm09aHu6aBVBQLHF7SJHMWhIZRESCFvWIdusmXnnfe8fDxOJEDJJZuZ9Z57v55w5kORl3l8CzHeeXTEMA0RERLJQzS6AiIgokxh8REQkFQYfERFJhcFHRERSYfAREZFUGHxERCQVez9f51oHIiLKRkpfX2CLj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIwLg9Xp7fbx69WrcdtttJ123du1aTJkyBdOmTcP06dOxZcuWTJVIRCliN7sAomwyZ84cLFiwAIqioKamBt/97nexZ88es8siogFg8BENwIktw2AwCEVRTKyGiAaDwUcEIBwOY9q0aT0ft7W1YcGCBae8ds2aNfj5z3+OlpYWbNiwIVMlElGKKIZhnO7rp/0iUa7wer0IBAI9H69evRpVVVV48skn+/wz7733Hu677z7885//zESJRDQwfXbHsMVH0ohGowiFQggGg/D7/YjFYjAMAzabDYZhwOfzwePxwO12J/V8X/va11BXVwefz4fi4uI0V09EqcLgo5xmGAaCwSBaWlrQ1tYGAFAUBU6nE6qq9lxjGAbq6+thGAY8Hg8CgQBO1RtSW1uLsWPHQlEUfPzxx4hGozjjjDMy+j0R0dAw+ChnxWIxHD58GD6fDw6HA16vt8/JKIqioKCgAIBoGSaCMhAI9JrQ8vLLL+OFF16Aw+GA2+3G3/72N05wIcoyHOOjnNTZ2Ym6ujoYhoH8/PxBhVN3dzcikQgqKipw5plnMuCIskuf/2EZfJRz2tvbUVtbC7fbDYfDMaTn0nUdfr8fJSUlGD16NMOPKHv0+Z+VO7dQTvH7/aitrYXH4xly6AGAqqoYNmwYfD4fDh8+nIIKichsDD7KGZqmoa6uDi6XC3Z76oavE+N/TU1N6OzsTNnzEpE5GHyUM44cOYJ4PA6n05ny51YUBR6PBwcOHICmaSl/fiLKHAYf5YRQKITW1lbk5+en7R4OhwOapqG1tTVt9yCi9GPwUU7w+Xyw2Wxpn3zi8XjQ1NQEXdfTeh8iSh8GH2W9eDyO1tbWpHdcGQqbzQZN09DV1ZX2exFRenABO2W9cDgMwzB6dmJJN5vNBr/fj8LCwozcLxcZBtDcDBw8CPh8QHc3EA4Dmgbk5QEul3iceSYwahQwbJjZFVMuYfBR1guHwxm9n9PpZItvgDo6gM8+A/btA3bsEL9Go+JrmgYoCqCq4lddFw9FAex2IB4HioqAiROBKVOAc84Bzj0XSMFqFZIUg4+yXiAQGPTyhUgkgmAwiFgsBk3TYBgG7HY77HY78vLyUFBQcNK4od1uRyAQgK7rGWtlZiPDAPbsAdatA9599/PP5ecDhYXJB5dhAJEIUFUFbNkiAjI/H7jqKuDSS4GSkvR9D5SbuHMLZb19+/YhHA4jLy/vtNdFIhE0NzejubkZ9fX1OHToEDo6OmCz2XpdZxgGFEXpCcHKykqMGTMGFRUVKC0tRWFhIQKBAKZNm5bS9YK5IhQSAfXSS8DhwyLgSkqAL/yYhyQcBo4dE7+fPRu48kpg8mQRikTHccsyyl179+5FJBI5ZfBFo1Hs3LkTW7ZsQXNzM1RVha7rcDgccLlccDqdp50JqmkaIpEIIpFIT+vO4XDgggsuwDXXXINRo0al7fvKNpoGrFkD/OlPYsxu2DDxSOdEW10HWltFt2l5ObB0qegOJQKDj3LZ/v370dXV1WtWZ2trK6qqqrBt2zbEYjG43W54PJ6ULHeIxWKIRCJobGzEtGnTsHDhQkyfPl3q1t/Bg8DKlWLsrqxMTFDJtI4OoLNTtP5uvBHweDJfA1kKg49yV3NzMw4fPoyCggIcOHAAmzdvRn19fc8+m6kOpER4+v1+tLW1obu7G8OHD8e1116LK6+8Mi07x1hVopX33HMi7IqL09vC6088DjQ0iDp+9jO2/iTH4KPc5ff7UVNTg61bt+LDDz9EXl4evF5v2iaeqKqKWCyG7u7uns+FQiEcO3YMo0aNwj333IPx48en5d5WcugQ8NBDwN69opvRSnl/YuvvBz8AMrDEk6yHwUe5a/v27XjmmWfQ1taGoqKitM+0tNlsCIfDp9yz89ixY4hEIvje976H6667Lmdbf7t2AffcI2ZclpSY28rrSzwOHD0qlkHcdx9w/JxhkgeDj3JPMBjEs88+i7Vr16KiogKFhYUZ2UpMVVUEAoE+vx6LxdDU1JSzrb+PPwbuvVeMoQ0fbnY1p2cYoutz1CjgN78R6wFJGgw+yi1NTU2488470djYiPLycjgcDuTn5yMej6f1vqqqIhqNIppYfX0aidbf8uXLMXfu3LTWlSmffCJaeoWFgNdrdjXJa2wERo4EVq2yflhTyjD4KHccOnQIP/3pTxEMBlFaWtrz+cQ5fOlq9SmK0tPa6+f/TY/E2sHbb78dV199dVaf4L57N7B8uQi8bOw2bGgAxowBHnggu0KbBo3BR7mhvr4ed9xxBzRNwxlnnNHra4qiwOv1Qtf1pINpIOx2O0Kh0IDP44tGo2hsbMRNN92E66+/PuV1ZUJDA3DrrWIxera2mBLdnlOmACtWcLG7BPoMPv7VU9Y4evQoli1bdsrQA8SOK6FQKC2TW2w2G6LR6KAOoXU6naioqMCzzz6LV199NeW1pZuuA488IiaLZDL0fv97BR98sKzn45qah7F9+y9Puq6jYw/WrZuN557LQ03Nw30+n6IAFRVijHLDhnRUTNmCwUdZoa2tDXfeeSfC4fApQy8hHo8jHA6ftA3ZUNhstp5F64PlcDgwcuRIPPbYY9i8eXPKasuE9euBmhrghF7ljLDZ8lBf/woiEd9pr8vLK8Ls2Y9j8uQ7+31ORREL7P/3f0Xrj+TE4CPLMwwDDz/8MHw+H0qS2JFY07Se8BtK609RlJ6W3lBCLyEvLw8lJSVYuXIlGrLkVffoUeCZZ0RYZHp4UlHsOO+8H+Ff/3r0tNe53aUoKZkBVU1u12uXS3RzPvKIaM2SfOTdY4myxubNm7F161acddZZSf8ZTdMQCAR6TXhJdtxPURSEw+Fe+3t6vV7MmDEDw044GC4YDGLLli09oThhwgScf/75p31ut9uNzs5OrFy5EqtWrbL06Q66Djz6qAgJl8ucGiZOvBWvvDIFU6f+LKXPW1oqWrEbNgDz56f0qSkLMPjI0nw+Hx599FGUlJQMOCQMw0A4HIbdbofD4ejZuswwjJNmfiZmbAKAruvYsGEDxo4di6985SsARFdrOBzuFXyKomD69Ok444wzEIvFsH79epSXl/d7QG1ZWRmqq6uxfv16LFiwYEDfUya9/roIh7PPNq8Gp3MYxo+/ATt3Pg6bLXXbr5zY5Tl9uth5huRh3bebJD3DMPDII48gGo3CM4QdhxNdn8FgEJFIBLFYrCfoEg9d1xGJRBAKhVBbWwtFUTBhwoSe5ygqKkJZWVmv5/V4PD3jjQ6HA8OHD0coFOq3HkVRUFpait/97neW7fKMxYDnn7fGriyTJt2Bzz77AzQtmNLndblEq3bt2pQ+LWUBBh9ZVqKLc+TIkSl5Pl3Xe/bYDAaDCAQCPY9wOIxYLIZ4PI6Ojo7TTqA5lUAggLa2NhQXFyd1vet43+HKlSszstvMQFVVAV1d1jjhwOUqwpgx38Vnn/0h5c9dUgJs3CjOECR5MPjIkqLRKJ588kkUFxdbehwMEFuUvfPOO5gxY8aA9uYsKyvDjh07sG3btjRWNzh//7s1Qi9h8uRlfc7uDIWa8Oc/V+LTTx9BdfWv8ec/VyIa7UrqeZ1OcZbf//1fKqslq+MYH1nStm3b4Pf7UVlZmfF7FxYW4uDBg0ldq+s63nnnHZxzzjkDPpRWURS43W689NJLmDVr1mBKTYtDh4BPPwUGMJcoLRYv/nw/VI+nDDfeeOpmmcczEtddd2TQ9xk+HHjxReCyy8zv1qXMsPZbaZLWX//61yGN6w3FyJEjoes69u7d2/M5n8+HpqamXtcZhoGtW7eisLAQEydOHNS9ioqKUF1djSNHBv/CnWqvvw7YbPKEQEGBWLaxe7fZlVCmMPjIcg4cOIDdu3djxIgRptxfURRcfPHFaGxsxCuvvIK1a9dix44dJwVxS0sL9u/fj8bGRrz22mt47bXXBhxgqqpCURRssMhWIpGIGPNKYrlkzlAUsRXb+vVmV0KZwr06yXKeeOIJrFu3DhUVFWaXkhHRaBSdnZ14+eWX4Tb5xNQ9e4Bly8TWXjLp7gbCYdHlKUtLVwLcq5OyQzAYxMaNG5PaoSVXOJ1OdHd3Y8uWLWaXgvp6sSenbJxOIBAA2trMroQygcFHlrJz505omgaHI7ntp3KF2+3GW2+9ZXYZ+PRTEQKyURSxQ02Sc5ooyzH4yFJqa2vTfpisFXm9XuzatSstxykNxM6d8p5VF48D+/ebXQVlAoOPLKW6uhr5+flml5FxTqcTwWAQx44dM62GUEicVG6l9XuZ5PGILdoo9zH4yDIMw8CePXvglbDJoSgKFEVBfX29aTUcOiTXMoYv8nrFkgaTG92UAVzATpbh8/kQCoVMW8ZgNl3XUVdXh+nTp5ty//Z2uV/0nU6g+XA34m9ugb3mY2DbNqCzU0z3jEbFuwK3WzwmTAC+/GVg6lRg/HjAzpfSbMK/LbKMAwcOQJG1uQExwWXHjh1YtGiRKfcfxOHy2a+rU0xlbWqCcqwNancJYjf8EPbI8XcBitK7CZx4Z/DOO8Af/yi+pmki/C66CJgzB5g7l0FocfzbIctobGyEJuWrr+DxeLDfxNkV0agkLT5dB44cEVNYfb7j37QBKCqgAFFNhTvZqa2JcPzsMzEz6E9/Ei3Cn/wEuPFGcfYRWQ7H+MgywuGw1C0+m82WkpPeB0vXczz4wiFgRzXwt78B770LtLZ8vo5BFYObCgB9oC+Lqiq2fnG7xQ+wqwt46CHg/POBRYuA99/P8R9s9mGLjywjEolY/iSGdFJVFdFo1LT7Oxw5OrFF10XrbseO3q071XbypVDgUIbY65Do5tR1sfHpW28BU6YAzz0HDHAjc0oPeV9lyHK6u7ulbvGpqmpqV6/TmYPB19EOvLZOtPR6xuz6ftkzoMCJFL35UFUgL0/cd/t2YMYMceS7Bc9flA2DjywjLy/P9AXcZtJ1HXYTJ0W4XDkUfLouFuW9th7o6ASg9PvNxQ0VdsRhR4rffCR2wdY04Be/AC65hFvEmIzBR5bhcrkseRp5pui6PqCDbFOtsjJHhqIC/t6tvCS7z4PwYIztIFQlTT8E2/Gu1Y8/Fq2/559Pz32oXww+sgyXyyV9iy8vL8+0+5eUiFaficOMQ9fRPqBW3okChheTbWk+lO/E1t9Pfwo88ECOvNvILgw+sozy8nJTu/rMFgqFMGbMGNPuryjAueeKUwqyks8HbNgoknsQk6TiUHGevTYNhZ1CovW3ciVw990Mvwxj8JFljB49WuoWXygUwtSpU02tYepUsWdn1mlvBzZtEi2pQc4MVqFjlHo4xYWd7obH6/z974F7783cfYnBR9ZRWloKl8uFWCxmdimmsNlsGDdunKk1jB2bhRNcurqAN14fUujFDfHnKtWGVFbWv0S9Tz0FrFqV2XtLjMFHlqEoCs477zwEsravbWh0Xcfo0aNNrWH0aNHrljUNb00DNr0BRGODDj0A8BtejLXVw66YcCRWou4VK4D16zN/fwkx+MhSpk2bhmAwaHYZGReNRuHxeFBcXGxqHSUlYsORrDmJvOojIBwZUugBgB9eXOF4M6lrj8TjuLKrC+Pb2jC2rQ3/GQgg2sc7hecjEYxva8P4tjY8f7pdeVRVHAh4881Z9MPPXgw+spRx48ZJuXuL3+/H+eefb4kF/AsXZsk4X3MzsG8fMMQfWdSww6nE8FXntn6vNQwDV/v9+LbTiX1FRdg7YgQChoFfnOIH1qbr+FUohA8LC7GtsBC/CoXQfrrlOg4HEAwCt98+lG+HkiDfKwxZ2qRJk6CqqnTjfOFwGBdffLHZZQAApk8HCgrEaTyWpWnAu+8CuoGhJl+rUYxvOv6JfKX/tH8rFoMLwI0uFwDApih4ND8fz0UiCH2h1bcpFsOlTieKVBUjVBWXOp14o79/1zYb8MYb7PJMMwYfWUpBQQHmzZsHn89ndikZE41G4XQ68fWvf93sUgCIhsfVVwMmHgbfv6qPgMjQuzgNQ2xK/U3nW0ldvzMex4VfWHIzTFVxtqqiNt57fPBoPI6zTqivUlVxNN7PGKKisMszAxh8ZDkLFixAPB6XZmmDz+fD5ZdfDo/HY3YpPS65RPxqyY10WltT0sUJAG1GIc637cVoW3LLGPpqX57q86f615tUyYkuz2XLkqqJBo7BR5ZzzjnnYMKECejo6DC7lLQzDAO6ruOKK64wu5ReSkqAf/93MYxmOf+qOZ7IQ0s+wwAC8OK7eeuS/jOTbDZUfWEj8S5dx2Fdx1hb79MeKm02HD7hncMRXUeF7eQTIU5JVYG1a4GmpqRro+Qx+MhyFEXBokWLpJjd2dbWhkmTJpm+jOFUbrpJnLBjqYku4RBw9OiQuzgBoMkoxUz7dsyyVyX9Z+Y4HAgBeOH4DM24YWBZMIjFeXnwfGFi0lyHA29Go2jXdbTrOt6MRjHX4UjuRonvb/XqpGuj5DH4JOD1ent9vHr1atx22219Xv/RRx/BZrPh73//e7pL69Ps2bPh8XgQstSrbmoZhoFQKIRFixaZXcoplZSICYatrRZa17d33/E+xKG19kKGC3Yljv90PzugBfuKomBNQQFe6u7G+LY2TGhvh0tRsCI//6Rri1QV/+XxYEZHB2Z0dOBejwdFAwlswwD+538AySZ6ZQKDj3qJx+O46667MHfuXFPryMvLw6233orW1tacHetraWnBxIkTMWvWLLNL6dM3vgHMnGmRHjddB3btwqlHz5JnGECLUYLbXc+iWB34BJKzbDa8Nnw49hUVoa6oCE94vcjrIz1/4HKhtqgItUVFPTNBk2a3i6m1mzYNuEY6PQYf9fLEE0/gmmuuQWlpqdmlYO7cufjyl7+MZksONA1Nd3c3dF3H3XffDVuy4z4mUBRgyRKLdHkeOQLEtdMeJJuMJqMUs+xV+IZjS4oKS6NIBPjtb82uIucw+CQQDocxbdq0nse9fWyIe/ToUaxZswY333xzhis8NUVRsGzZMqiqirClF5UNjGEYaG5uxo9//GNUVlaaXU6/El2eLS1ipr1pdu8ecgF+Ix+OQXRxmsbpBD75BDicwc2zJcDgk4Db7UZ1dXXP47777jvldXfccQcefPBBS7VAysrKsGTJErS0tORMl2dLSwsmTZqEb3/722aXkrRvfAO45hrx+mvKEgfDAHytQ5rUEjTc6DSG4ZeehwbVxWkKRRHN7e3bza4kp8h7+FmO0jQNzc3N8Pv90HUdNputZ8p8f1uBVVVV4dprrwUg1pZt3LgRdrvd9BfoefPm4e2330Z1dTXKy8tNrWWoIpEIdF3HXXfdZak3GP1RFOBHPxJn9b35JnDWWSmZWJm8YHBISxjCRh6OGUX4lftBTLXvTG1t6RYKAVVVQBa9UbI6Bl8O8Pv9qKqqQk1NDbq6umC323vt+ahpGlasWIHi4mJceOGF0L6wDinhwIEDPb9fvHgxrrjiCtNDDxBdnsuXL8ctt9wCn89n+kbOgxWNRtHS0oK77roLZ511ltnlDJiqAnfcISYZvvUWcPbZGQy/tjYxtjeIVn/QcOOYUYS73Y9hlvPjNBSXZqoKvP++2VXkFAZfFvP7/diwYQNqa2uhKEpPoEWj0V7XGYaBeDyO5uZm/OMf/8D27dsRi8UQi8XgSHZdkclKSkqwatUqLFmyBO3t7RgxYoTZJQ2IpmlobGzErbfeinnz5pldzqDZ7cDPfgZ4PGI7ycpK8bm08/nExBZ1YK3kLsOLLqMAv3I/mJ2hB4gf8M6dIvSzYmDS+pR+xk1yY1AlxxiGgR07duD111+HpmnQBzHoYrfb4XK5sHDhQpx99tlpqDI99u3bh6VLl0JRFBQVFZldTlJisRgaGhpwww034MYbb7TECQxDpevACy8Af/mL2NA67e9DXt8ItCQ/xqcbChqNMniUMO71rMIU+640F5gB27YBY8aYXUU26fM/Gie3ZBnDMLBu3Tps3LgR0Wh0UKEHiBZIIBDAH//4R1RVJb9zhdnGjx+Pxx57DA6HIys2su7u7kZjYyNuuummnAk9QOTP4sXAo48C+fli0ksfPeip0dYOqMn97AKGB4eMSlzk2Ibfe5fmRuipKlBTY3YVOYPBl0UMw8Crr76KnTt3puzYHk3T8Oabb2ZV+I0dOxaPP/44CgoKcPTo0UGHf7q1t7ejtbUVS5cuxfXXX58zoXeiiROBp58WMz4bGoD29jTcxDCAWBT9TWzRDQVH9ZGIwYlfuh/CPe7folDtSkNBJtA0ix+XkV0YfFnkgw8+wO7du1N+Vl0sFsOmTZtw8ODBlD5vOlVWVuKZZ57BJZdcgqNHj8Lv95tdUg9N03DkyBEUFBTg8ccfx4IFC8wuKa3cbrGvZ6L1d+gQ0NmZwm3ODP20Y1u6oaBZL8YRo+J4K+8OXOT8KLeGw+JxoLvb7CpyBsf4ssSxY8fw9NNP9zkjMxW8Xi+WLFmSNRNeANEK3rZtGx588EF0dXWhvLzc1BPc29vbEQgEsHDhQixevBhut9u0WswQiQDvvQe8+KLYS9puFwvgh7RyIxoF/vqXk3ZsSSxRABTMsn+Eq/Jex2TbrtwKvIRIBLjvPmDpUrMrySZ9/kvgrM4s8corryCe5m0zIpEINm/enFWzDhVFwcyZM7F69Wo89dRT2LRpEwoKClBYWJjRrsVIJAKfz4eRI0fi/vvvx6RJkzJ2bytxuYDLLgMuvRTYswdYt+74Qem6mATj9Q5+FqhhAN1wwm94EYELBUoQ/5H3Ei51vosSld2AlDy2+LJAQ0MDVq9enfIuzlOx2+1Yvnw5nE5n2u+VaoZhoKqqCn/4wx+wb98+qKqKkpKStLVgDcNAR0cHgsEg8vPz8Z3vfAcLFy6UrpXXn7Y24O23ga1bxfmxmvZ5N6jXK7pKVVU8FEV8LR4Xj2AQCAfjsL31T8RVOwqVLky0fYZLnO9hhv0TOJR0zqixkFgMuP9+4Cc/MbuSbNLnO18GXxZ4+eWXsXPnzoxs2eVwOHDZZZdh+vTpab9XuhiGgbq6Oqxbtw6bNm2CpmkYPnw4vF5vSlqBsVgMra2t0HUd5557Lq699lrMmjUrK98sZJquA42NwMGDwN69wL/+JbpEo1HxiMfFAeROpwjEceOAqVMMjLluNka5WzBc9edmV2Z/dB1YuRL4wQ/MriSbMPiyla7rWLFiRdq7OU9UVlZmmY2qh8rv9+Ptt9/Giy++iMbGRqiqCl3X4Xa74fV64XQ6TxuG8XgcgUAAgUCg5zqbzYZvfvObmD9/PsaMGZOTszUtp7xcJKOJ47emUlXgmWeAK680u5JswjG+bOXz+WCz2TIafD6fL6m9PbNBQUEBFixYgPnz56O9vR319fWoq6tDdXU1du/ejdbWVthsNiiK0tOi/uLvx40bhylTpuC8887DmDFjcOaZZ2bVPps5YfJk4MMPRVNQRroOTJlidhU5g8FncQ0NDRk/lcBms+HYsWMoKSnJ6H3TKbHTS1FRES644AIsXLgQhmGgvb0d7e3tPVu46boOh8MBh8MBt9uN8vJyhpwVfPWrwAcfmF2FOXRdTIsdPdrsSnIGg8/i2traMjKp5USKoqC9vT2ngu9UTgxDsrgLLhCDfhbdrCCtNE209tilnjLZ35eV4zIdegnpXC9INGBTp6Z5TzQL03XgoovMriKnMPgszp6Rre9Pxu49spTKSjHdU8YWn8cDZPEsayti8FncsGHDMh5+hmGgoKAgo/ckOi1FAWbMEOvZZGIYoqV74YVmV5JTGHwWZ8bkCk3TUFpamtF7EvXr5pvF1jAyiUaBf/s3ceQ9pQyDz+LKysoyPs43fPhw07pYifo0d66Y4CLTWJ/LJY69p5Ri8Fmcw+HI6EGxdrsdU7heiKzIbgduuUWe2Y2aJoI+i/bOzRYMvixw0UUXZWw7LMMwsnq7MspxixeL4JNhkouiiL052fuScgy+LDB27Fi4MjC2oaoqxo0bB6/Xm/Z7EQ1KWZk4+iHXJ7kkgn3xYlPLyFUMviygKAquuuqqtJ+TZ7fbcfnll6f1HkRDtnx57i9t0HWxL+fIkWZXkpMYfFli9OjRmDx5ctomnTgcDnzrW9/iMgayvgsvBG64IYVHvFtMLCaOsl+1yuxKchaDL4vMmzcPJSUlKV/e4HA48KUvfYmTWih7/OY3QHGxmO6fSwxD7PrnaJMAAAUoSURBVMv59NMAt9JLGwZfFnE4HPj+97+P0tLSlLX8HA4HJk6ciPnz5/N4HcoeHg+werWY+JFLXZ7xuJjFecUVZleS03geXxaKxWJ44403UFNTM+g9NRVFgc1mw5w5czBz5kyGHmWnpUuB558XraRsF4uJI+lratjaSw0eRJuL6uvrsWbNGkQiEUST7PJRFAV2ux2lpaW4+uqreTIBZbdQSJzc0NiY3Wf1JVqtzz/P1l7qMPhylWEYqKurw/vvv49Dhw7B4XBA1/Veu70k1gDG43FMnDgRs2fPRnl5uVklE6XW/v3A178OdHVlZ/glQu+ee4Bly8ytJbcw+GSgaRqam5vR0NAAv98PTdPgcDgwYsQIVFRUoLi4OCdOVSc6ya5dwJw5ogWYTeGXCL1bbgH++7/NrSX3MPiIKMd98glw+eVAMJgd4ZcIvR/+EHjgAXm2YsscBh8RSWD3bjErsqNDLHK3apjE4+LXu+8WC/KtWmd2Y/ARkSQOHQIWLgRqaz9fF2cVifP1nE7g4YfFQnxKlz6DjwM+RJRbzj4beP994Oc/B1RVLHK3wi4viVbehRcCH33E0DMRW3xElLv27AG+/31zW38ntvJ+/WsxpsdJZpnArk4ikpSmAY89Btx/v5hQYhiZOepH18WidKcTmDoVeO45YNSo9N+XEhh8RCS55maxzdlTTwHhMBCJiFBK9cQSTRPPqSjiCKUlS4DZszmBJfMYfEREAEQwbdokWoHbt4uPE61AVR14QOm6GL+Lx8WJ6W63WJe3eLE4P5DMwuAjIjpJbS3w6qvAli1AdTXQ2Qnk5YkWoWF83nJLSLxe6roIOEAE57nnilbdJZcAl13GU9OtgcFHRNQvnw/YsUMsht+2TQRhJCIedrs4FcLlAsaPB2bOFGN348dba8kEJTD4iIhIKlzHR0REBDD4iIhIMgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKpMPiIiEgqDD4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpMLgIyIiqTD4iIhIKgw+IiKSCoOPiIikwuAjIiKp2Pv5upKRKoiIiDKELT4iIpIKg4+IiKTC4CMiIqkw+IiISCoMPiIikgqDj4iIpPL/T9pjc6hmaRYAAAAASUVORK5CYII=\n",
-      "text/plain": [
-       "<Figure size 432x288 with 1 Axes>"
-      ]
-     },
-     "metadata": {
-      "needs_background": "light"
-     },
-     "output_type": "display_data"
-    }
-   ],
+   "outputs": [],
    "source": [
     "plot_3d_mol_as_scatter(xyz)"
    ]
diff --git a/ipython/Tools/TS Guesses.ipynb b/ipython/Tools/TS Guesses.ipynb
index 17b9f88738..af6ce4a3eb 100644
--- a/ipython/Tools/TS Guesses.ipynb	
+++ b/ipython/Tools/TS Guesses.ipynb	
@@ -12,26 +12,15 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": null,
    "id": "5bae17f3",
    "metadata": {},
    "outputs": [],
-   "source": [
-    "from arc.imports import settings\n",
-    "from arc.job.adapters.common import all_families_ts_adapters, ts_adapters_by_rmg_family\n",
-    "from arc.plotter import draw_structure\n",
-    "from arc.reaction import ARCReaction\n",
-    "from arc.rmgdb import make_rmg_database_object, load_families_only\n",
-    "from arc.species import ARCSpecies\n",
-    "from arc.species.converter import xyz_to_str\n",
-    "from arc.scheduler import Scheduler\n",
-    "\n",
-    "%matplotlib inline"
-   ]
+   "source": "from arc.imports import settings\nfrom arc.job.adapters.common import all_families_ts_adapters, ts_adapters_by_rmg_family\nfrom arc.level import Level\nfrom arc.plotter import draw_structure\nfrom arc.reaction import ARCReaction\nfrom arc.species import ARCSpecies\nfrom arc.species.converter import xyz_to_str\nfrom arc.scheduler import Scheduler\n\n%matplotlib inline"
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": null,
    "id": "c2be78fd",
    "metadata": {},
    "outputs": [],
@@ -45,7 +34,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": null,
    "id": "1571f9b1",
    "metadata": {},
    "outputs": [],
@@ -56,18 +45,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
-   "id": "826ede9f",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "rmgdb = make_rmg_database_object()\n",
-    "load_families_only(rmgdb)"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": null,
    "id": "9a2a6dc9",
    "metadata": {},
    "outputs": [],
@@ -78,432 +56,38 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": null,
    "id": "348b3b0a",
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "Considering the following reaction:\n",
-      "\n",
-      "[OH] + CCC <=> O + [CH2]CC\n",
-      "\n",
-      "Reaction family: H_Abstraction\n"
-     ]
-    }
-   ],
-   "source": [
-    "r_species, p_species = list(), list()\n",
-    "for r_smiles in [r1_smiles, r2_smiles]:\n",
-    "    if r_smiles is not None:\n",
-    "        r_species.append(generate_species_from_smiles(r_smiles))\n",
-    "for p_smiles in [p1_smiles, p2_smiles]:\n",
-    "    if p_smiles is not None:\n",
-    "        p_species.append(generate_species_from_smiles(p_smiles))\n",
-    "\n",
-    "rxn = ARCReaction(r_species=r_species, p_species=p_species)\n",
-    "rxn.ts_species = ARCSpecies(label='TS', is_ts=True)\n",
-    "rxn.determine_family(rmg_database=rmgdb)\n",
-    "\n",
-    "print(f'Considering the following reaction:\\n\\n{rxn.label}\\n\\nReaction family: {rxn.family.label}')"
-   ]
+   "outputs": [],
+   "source": "r_species, p_species = list(), list()\nfor r_smiles in [r1_smiles, r2_smiles]:\n    if r_smiles is not None:\n        r_species.append(generate_species_from_smiles(r_smiles))\nfor p_smiles in [p1_smiles, p2_smiles]:\n    if p_smiles is not None:\n        p_species.append(generate_species_from_smiles(p_smiles))\n\nrxn = ARCReaction(r_species=r_species, p_species=p_species)\nrxn.ts_species = ARCSpecies(label='TS', is_ts=True)\nrxn.family, rxn.family_own_reverse = rxn.determine_family()\n\nprint(f'Considering the following reaction:\\n\\n{rxn.label}\\n\\nReaction family: {rxn.family}')"
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": null,
    "id": "2db720a0",
    "metadata": {
     "scrolled": false
    },
    "outputs": [],
-   "source": [
-    "sched = Scheduler(project='TS_search_1',\n",
-    "                  species_list=rxn.r_species + rxn.p_species,\n",
-    "                  rmg_database=rmgdb,\n",
-    "                  conformer_level='wb97xd/Def2TZVP',\n",
-    "                  opt_level='wb97xd/Def2TZVP',\n",
-    "                  sp_level='wb97xd/Def2TZVP',\n",
-    "                  freq_level='wb97xd/Def2TZVP',\n",
-    "                  ess_settings=None,\n",
-    "                  testing=True,\n",
-    "                  project_directory='/home/alon/Dropbox (MIT)/runs/TS_search_1',\n",
-    "                 )\n",
-    "sched.job_dict = {'TS': {'tsg': {}}}"
-   ]
+   "source": "level = Level(repr='wb97xd/Def2TZVP')\n\nsched = Scheduler(project='TS_search_1',\n                  species_list=rxn.r_species + rxn.p_species,\n                  conformer_opt_level=level,\n                  conformer_sp_level=level,\n                  opt_level=level,\n                  sp_level=level,\n                  freq_level=level,\n                  ess_settings=None,\n                  testing=True,\n                  project_directory='TS_search_1',\n                 )\nsched.job_dict = {'TS': {'tsg': {}}}"
   },
   {
    "cell_type": "code",
-   "execution_count": 8,
+   "execution_count": null,
    "id": "0ff84291",
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "Considering the following TS Search methods: ['heuristics', 'autotst', 'gcn', 'kinbot']\n",
-      "dihedral_increment: 90\n"
-     ]
-    },
-    {
-     "name": "stderr",
-     "output_type": "stream",
-     "text": [
-      "/bin/bash: -c: line 0: syntax error near unexpected token `('\n",
-      "/bin/bash: -c: line 0: `source ~/.bashrc; /home/alon/anaconda3/envs/tst_env/bin/python /mnt/Data/Code/ARC/arc/job/adapters/ts/scripts/autotst_script.py [OH]+CCC_O+[CH2]CC /home/alon/Dropbox (MIT)/runs/TS_search_1/calcs/TS_guesses/TS/tsg1/ts_results.yml'\n",
-      "WARNING:arc:AutoTST subprocess did not give a successful return code for ARCReaction(label=\"[OH] + CCC <=> O + [CH2]CC\", rmg_reaction=\"[OH] + CCC <=> O + [CH2]CC\", multiplicity=2, charge=0) in the forward direction.\n",
-      "Got return code: 1\n",
-      "stdout: None\n",
-      "stderr: None\n",
-      "/bin/bash: -c: line 0: syntax error near unexpected token `('\n",
-      "/bin/bash: -c: line 0: `source ~/.bashrc; /home/alon/anaconda3/envs/tst_env/bin/python /mnt/Data/Code/ARC/arc/job/adapters/ts/scripts/autotst_script.py O+[CH2]CC_[OH]+CCC /home/alon/Dropbox (MIT)/runs/TS_search_1/calcs/TS_guesses/TS/tsg1/ts_results.yml'\n",
-      "WARNING:arc:AutoTST subprocess did not give a successful return code for ARCReaction(label=\"[OH] + CCC <=> O + [CH2]CC\", rmg_reaction=\"[OH] + CCC <=> O + [CH2]CC\", multiplicity=2, charge=0) in the reverse direction.\n",
-      "Got return code: 1\n",
-      "stdout: None\n",
-      "stderr: None\n"
-     ]
-    },
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "Got 4 TS Guesses (before clustering).\n",
-      "Got 4 TS Guesses (after clustering).\n"
-     ]
-    }
-   ],
-   "source": [
-    "tsg_index = 0\n",
-    "ts_adapters = settings['ts_adapters']\n",
-    "ts_adapters = [a.lower() for a in ts_adapters]\n",
-    "\n",
-    "print(f'Considering the following TS Search methods: {ts_adapters}')\n",
-    "\n",
-    "for method in ts_adapters:\n",
-    "    if method in all_families_ts_adapters or \\\n",
-    "            (rxn.family is not None\n",
-    "             and rxn.family.label in list(ts_adapters_by_rmg_family.keys())\n",
-    "             and method in ts_adapters_by_rmg_family[rxn.family.label]):\n",
-    "        sched.run_job(job_type='tsg',\n",
-    "                      job_adapter=method,\n",
-    "                      reactions=[rxn],\n",
-    "                      tsg=tsg_index,\n",
-    "                      dihedral_increment=dihedral_increment,\n",
-    "                     )\n",
-    "        tsg_index += 1\n",
-    "\n",
-    "print(f'Got {len(rxn.ts_species.ts_guesses)} TS Guesses (before clustering).')\n",
-    "\n",
-    "rxn.ts_species.cluster_tsgs()\n",
-    "        \n",
-    "print(f'Got {len(rxn.ts_species.ts_guesses)} TS Guesses (after clustering).')"
-   ]
+   "outputs": [],
+   "source": "tsg_index = 0\nts_adapters = settings['ts_adapters']\nts_adapters = [a.lower() for a in ts_adapters]\n\nprint(f'Considering the following TS Search methods: {ts_adapters}')\n\nfor method in ts_adapters:\n    if method in all_families_ts_adapters or \\\n            (rxn.family is not None\n             and rxn.family in list(ts_adapters_by_rmg_family.keys())\n             and method in ts_adapters_by_rmg_family[rxn.family]):\n        sched.run_job(job_type='tsg',\n                      job_adapter=method,\n                      reactions=[rxn],\n                      tsg=tsg_index,\n                      dihedral_increment=dihedral_increment,\n                     )\n        tsg_index += 1\n\nprint(f'Got {len(rxn.ts_species.ts_guesses)} TS Guesses (before clustering).')\n\nrxn.ts_species.cluster_tsgs()\n        \nprint(f'Got {len(rxn.ts_species.ts_guesses)} TS Guesses (after clustering).')"
   },
   {
    "cell_type": "code",
-   "execution_count": 9,
+   "execution_count": null,
    "id": "91452d76",
    "metadata": {
     "scrolled": false
    },
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "heuristics 0 (success: True):\n",
-      "coordinates:\n",
-      "O       1.43737500    1.58112542    0.11736861\n",
-      "H       1.91964767    0.94498355   -0.43190146\n",
-      "C      -0.56973933    0.33273535    0.85620513\n",
-      "C      -0.56973933   -1.07932962    0.29545433\n",
-      "C      -0.56973933   -1.07932962   -1.22387717\n",
-      "H       0.49495302    0.99495540    0.46428247\n",
-      "H      -0.56974009    0.30632916    1.95037139\n",
-      "H      -1.45698281    0.88458558    0.52960174\n",
-      "H      -1.45215391   -1.61664362    0.66021023\n",
-      "H       0.31267551   -1.61664346    0.66021017\n",
-      "H       0.31750400   -0.57210949   -1.61622422\n",
-      "H      -1.45698284   -0.57210981   -1.61622424\n",
-      "H      -0.56973903   -2.10599218   -1.60316723\n"
-     ]
-    },
-    {
-     "data": {
-      "application/3dmoljs_load.v0": "<div id=\"3dmolviewer_16580360278266184\"  style=\"position: relative; width: 400px; height: 400px\">\n        <p id=\"3dmolwarning_16580360278266184\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n        </div>\n<script>\n\nvar loadScriptAsync = function(uri){\n  return new Promise((resolve, reject) => {\n    var tag = document.createElement('script');\n    tag.src = uri;\n    tag.async = true;\n    tag.onload = () => {\n      resolve();\n    };\n  var firstScriptTag = document.getElementsByTagName('script')[0];\n  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n});\n};\n\nif(typeof $3Dmolpromise === 'undefined') {\n$3Dmolpromise = null;\n  $3Dmolpromise = loadScriptAsync('https://3dmol.org/build/3Dmol.js');\n}\n\nvar viewer_16580360278266184 = null;\nvar warn = document.getElementById(\"3dmolwarning_16580360278266184\");\nif(warn) {\n    warn.parentNode.removeChild(warn);\n}\n$3Dmolpromise.then(function() {\nviewer_16580360278266184 = $3Dmol.createViewer($(\"#3dmolviewer_16580360278266184\"),{backgroundColor:\"white\"});\nviewer_16580360278266184.zoomTo();\n\tviewer_16580360278266184.addModel(\"\\n     RDKit          3D\\n\\n 13 11  0  0  0  0  0  0  0  0999 V2000\\n    1.4374    1.5811    0.1174 O   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.9196    0.9450   -0.4319 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5697    0.3327    0.8562 C   0  0  0  0  0  3  0  0  0  0  0  0\\n   -0.5697   -1.0793    0.2955 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5697   -1.0793   -1.2239 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.4950    0.9950    0.4643 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5697    0.3063    1.9504 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4570    0.8846    0.5296 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4522   -1.6166    0.6602 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3127   -1.6166    0.6602 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3175   -0.5721   -1.6162 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4570   -0.5721   -1.6162 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5697   -2.1060   -1.6032 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1  2  1  0\\n  1  6  1  0\\n  3  4  1  0\\n  3  7  1  0\\n  3  8  1  0\\n  4  5  1  0\\n  4  9  1  0\\n  4 10  1  0\\n  5 11  1  0\\n  5 12  1  0\\n  5 13  1  0\\nM  RAD  1   3   2\\nM  END\\n\",\"sdf\");\n\tviewer_16580360278266184.setStyle({\"stick\": {}});\n\tviewer_16580360278266184.zoomTo();\nviewer_16580360278266184.render();\n});\n</script>",
-      "text/html": [
-       "<div id=\"3dmolviewer_16580360278266184\"  style=\"position: relative; width: 400px; height: 400px\">\n",
-       "        <p id=\"3dmolwarning_16580360278266184\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n",
-       "        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n",
-       "        </div>\n",
-       "<script>\n",
-       "\n",
-       "var loadScriptAsync = function(uri){\n",
-       "  return new Promise((resolve, reject) => {\n",
-       "    var tag = document.createElement('script');\n",
-       "    tag.src = uri;\n",
-       "    tag.async = true;\n",
-       "    tag.onload = () => {\n",
-       "      resolve();\n",
-       "    };\n",
-       "  var firstScriptTag = document.getElementsByTagName('script')[0];\n",
-       "  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n",
-       "});\n",
-       "};\n",
-       "\n",
-       "if(typeof $3Dmolpromise === 'undefined') {\n",
-       "$3Dmolpromise = null;\n",
-       "  $3Dmolpromise = loadScriptAsync('https://3dmol.org/build/3Dmol.js');\n",
-       "}\n",
-       "\n",
-       "var viewer_16580360278266184 = null;\n",
-       "var warn = document.getElementById(\"3dmolwarning_16580360278266184\");\n",
-       "if(warn) {\n",
-       "    warn.parentNode.removeChild(warn);\n",
-       "}\n",
-       "$3Dmolpromise.then(function() {\n",
-       "viewer_16580360278266184 = $3Dmol.createViewer($(\"#3dmolviewer_16580360278266184\"),{backgroundColor:\"white\"});\n",
-       "viewer_16580360278266184.zoomTo();\n",
-       "\tviewer_16580360278266184.addModel(\"\\n     RDKit          3D\\n\\n 13 11  0  0  0  0  0  0  0  0999 V2000\\n    1.4374    1.5811    0.1174 O   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.9196    0.9450   -0.4319 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5697    0.3327    0.8562 C   0  0  0  0  0  3  0  0  0  0  0  0\\n   -0.5697   -1.0793    0.2955 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5697   -1.0793   -1.2239 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.4950    0.9950    0.4643 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5697    0.3063    1.9504 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4570    0.8846    0.5296 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4522   -1.6166    0.6602 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3127   -1.6166    0.6602 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3175   -0.5721   -1.6162 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4570   -0.5721   -1.6162 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5697   -2.1060   -1.6032 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1  2  1  0\\n  1  6  1  0\\n  3  4  1  0\\n  3  7  1  0\\n  3  8  1  0\\n  4  5  1  0\\n  4  9  1  0\\n  4 10  1  0\\n  5 11  1  0\\n  5 12  1  0\\n  5 13  1  0\\nM  RAD  1   3   2\\nM  END\\n\",\"sdf\");\n",
-       "\tviewer_16580360278266184.setStyle({\"stick\": {}});\n",
-       "\tviewer_16580360278266184.zoomTo();\n",
-       "viewer_16580360278266184.render();\n",
-       "});\n",
-       "</script>"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    },
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "\n",
-      "\n",
-      "\n",
-      "\n",
-      "heuristics 1 (success: True):\n",
-      "coordinates:\n",
-      "O       1.44990567    1.56370672    0.12197742\n",
-      "H       1.17292977    1.98298507   -0.70654599\n",
-      "C      -0.55720866    0.31531664    0.86081394\n",
-      "C      -0.55720866   -1.09674832    0.30006314\n",
-      "C      -0.55720866   -1.09674832   -1.21926836\n",
-      "H       0.50748369    0.97753670    0.46889128\n",
-      "H      -0.55720941    0.28891046    1.95498020\n",
-      "H      -1.44445213    0.86716687    0.53421055\n",
-      "H      -1.43962324   -1.63406232    0.66481904\n",
-      "H       0.32520618   -1.63406216    0.66481898\n",
-      "H       0.33003468   -0.58952819   -1.61161541\n",
-      "H      -1.44445217   -0.58952851   -1.61161543\n",
-      "H      -0.55720836   -2.12341088   -1.59855842\n"
-     ]
-    },
-    {
-     "data": {
-      "application/3dmoljs_load.v0": "<div id=\"3dmolviewer_16580360278298295\"  style=\"position: relative; width: 400px; height: 400px\">\n        <p id=\"3dmolwarning_16580360278298295\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n        </div>\n<script>\n\nvar loadScriptAsync = function(uri){\n  return new Promise((resolve, reject) => {\n    var tag = document.createElement('script');\n    tag.src = uri;\n    tag.async = true;\n    tag.onload = () => {\n      resolve();\n    };\n  var firstScriptTag = document.getElementsByTagName('script')[0];\n  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n});\n};\n\nif(typeof $3Dmolpromise === 'undefined') {\n$3Dmolpromise = null;\n  $3Dmolpromise = loadScriptAsync('https://3dmol.org/build/3Dmol.js');\n}\n\nvar viewer_16580360278298295 = null;\nvar warn = document.getElementById(\"3dmolwarning_16580360278298295\");\nif(warn) {\n    warn.parentNode.removeChild(warn);\n}\n$3Dmolpromise.then(function() {\nviewer_16580360278298295 = $3Dmol.createViewer($(\"#3dmolviewer_16580360278298295\"),{backgroundColor:\"white\"});\nviewer_16580360278298295.zoomTo();\n\tviewer_16580360278298295.addModel(\"\\n     RDKit          3D\\n\\n 13 11  0  0  0  0  0  0  0  0999 V2000\\n    1.4499    1.5637    0.1220 O   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.1729    1.9830   -0.7065 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5572    0.3153    0.8608 C   0  0  0  0  0  3  0  0  0  0  0  0\\n   -0.5572   -1.0967    0.3001 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5572   -1.0967   -1.2193 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.5075    0.9775    0.4689 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5572    0.2889    1.9550 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4445    0.8672    0.5342 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4396   -1.6341    0.6648 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3252   -1.6341    0.6648 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3300   -0.5895   -1.6116 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4445   -0.5895   -1.6116 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5572   -2.1234   -1.5986 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1  2  1  0\\n  1  6  1  0\\n  3  4  1  0\\n  3  7  1  0\\n  3  8  1  0\\n  4  5  1  0\\n  4  9  1  0\\n  4 10  1  0\\n  5 11  1  0\\n  5 12  1  0\\n  5 13  1  0\\nM  RAD  1   3   2\\nM  END\\n\",\"sdf\");\n\tviewer_16580360278298295.setStyle({\"stick\": {}});\n\tviewer_16580360278298295.zoomTo();\nviewer_16580360278298295.render();\n});\n</script>",
-      "text/html": [
-       "<div id=\"3dmolviewer_16580360278298295\"  style=\"position: relative; width: 400px; height: 400px\">\n",
-       "        <p id=\"3dmolwarning_16580360278298295\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n",
-       "        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n",
-       "        </div>\n",
-       "<script>\n",
-       "\n",
-       "var loadScriptAsync = function(uri){\n",
-       "  return new Promise((resolve, reject) => {\n",
-       "    var tag = document.createElement('script');\n",
-       "    tag.src = uri;\n",
-       "    tag.async = true;\n",
-       "    tag.onload = () => {\n",
-       "      resolve();\n",
-       "    };\n",
-       "  var firstScriptTag = document.getElementsByTagName('script')[0];\n",
-       "  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n",
-       "});\n",
-       "};\n",
-       "\n",
-       "if(typeof $3Dmolpromise === 'undefined') {\n",
-       "$3Dmolpromise = null;\n",
-       "  $3Dmolpromise = loadScriptAsync('https://3dmol.org/build/3Dmol.js');\n",
-       "}\n",
-       "\n",
-       "var viewer_16580360278298295 = null;\n",
-       "var warn = document.getElementById(\"3dmolwarning_16580360278298295\");\n",
-       "if(warn) {\n",
-       "    warn.parentNode.removeChild(warn);\n",
-       "}\n",
-       "$3Dmolpromise.then(function() {\n",
-       "viewer_16580360278298295 = $3Dmol.createViewer($(\"#3dmolviewer_16580360278298295\"),{backgroundColor:\"white\"});\n",
-       "viewer_16580360278298295.zoomTo();\n",
-       "\tviewer_16580360278298295.addModel(\"\\n     RDKit          3D\\n\\n 13 11  0  0  0  0  0  0  0  0999 V2000\\n    1.4499    1.5637    0.1220 O   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.1729    1.9830   -0.7065 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5572    0.3153    0.8608 C   0  0  0  0  0  3  0  0  0  0  0  0\\n   -0.5572   -1.0967    0.3001 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5572   -1.0967   -1.2193 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.5075    0.9775    0.4689 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5572    0.2889    1.9550 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4445    0.8672    0.5342 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4396   -1.6341    0.6648 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3252   -1.6341    0.6648 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3300   -0.5895   -1.6116 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4445   -0.5895   -1.6116 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5572   -2.1234   -1.5986 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1  2  1  0\\n  1  6  1  0\\n  3  4  1  0\\n  3  7  1  0\\n  3  8  1  0\\n  4  5  1  0\\n  4  9  1  0\\n  4 10  1  0\\n  5 11  1  0\\n  5 12  1  0\\n  5 13  1  0\\nM  RAD  1   3   2\\nM  END\\n\",\"sdf\");\n",
-       "\tviewer_16580360278298295.setStyle({\"stick\": {}});\n",
-       "\tviewer_16580360278298295.zoomTo();\n",
-       "viewer_16580360278298295.render();\n",
-       "});\n",
-       "</script>"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    },
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "\n",
-      "\n",
-      "\n",
-      "\n",
-      "heuristics 2 (success: True):\n",
-      "coordinates:\n",
-      "O       1.44703222    1.55623294    0.10154324\n",
-      "H       1.34416191    2.42835646    0.51115195\n",
-      "C      -0.56008211    0.30784287    0.84037976\n",
-      "C      -0.56008211   -1.10422209    0.27962895\n",
-      "C      -0.56008211   -1.10422209   -1.23970255\n",
-      "H       0.50461024    0.97006292    0.44845710\n",
-      "H      -0.56008286    0.28143668    1.93454601\n",
-      "H      -1.44732558    0.85969310    0.51377636\n",
-      "H      -1.44249668   -1.64153610    0.64438485\n",
-      "H       0.32233273   -1.64153593    0.64438479\n",
-      "H       0.32716123   -0.59700196   -1.63204960\n",
-      "H      -1.44732561   -0.59700228   -1.63204962\n",
-      "H      -0.56008181   -2.13088465   -1.61899260\n"
-     ]
-    },
-    {
-     "data": {
-      "application/3dmoljs_load.v0": "<div id=\"3dmolviewer_16580360278336098\"  style=\"position: relative; width: 400px; height: 400px\">\n        <p id=\"3dmolwarning_16580360278336098\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n        </div>\n<script>\n\nvar loadScriptAsync = function(uri){\n  return new Promise((resolve, reject) => {\n    var tag = document.createElement('script');\n    tag.src = uri;\n    tag.async = true;\n    tag.onload = () => {\n      resolve();\n    };\n  var firstScriptTag = document.getElementsByTagName('script')[0];\n  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n});\n};\n\nif(typeof $3Dmolpromise === 'undefined') {\n$3Dmolpromise = null;\n  $3Dmolpromise = loadScriptAsync('https://3dmol.org/build/3Dmol.js');\n}\n\nvar viewer_16580360278336098 = null;\nvar warn = document.getElementById(\"3dmolwarning_16580360278336098\");\nif(warn) {\n    warn.parentNode.removeChild(warn);\n}\n$3Dmolpromise.then(function() {\nviewer_16580360278336098 = $3Dmol.createViewer($(\"#3dmolviewer_16580360278336098\"),{backgroundColor:\"white\"});\nviewer_16580360278336098.zoomTo();\n\tviewer_16580360278336098.addModel(\"\\n     RDKit          3D\\n\\n 13 11  0  0  0  0  0  0  0  0999 V2000\\n    1.4470    1.5562    0.1015 O   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.3442    2.4284    0.5112 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5601    0.3078    0.8404 C   0  0  0  0  0  3  0  0  0  0  0  0\\n   -0.5601   -1.1042    0.2796 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5601   -1.1042   -1.2397 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.5046    0.9701    0.4485 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5601    0.2814    1.9345 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4473    0.8597    0.5138 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4425   -1.6415    0.6444 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3223   -1.6415    0.6444 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3272   -0.5970   -1.6320 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4473   -0.5970   -1.6320 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5601   -2.1309   -1.6190 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1  2  1  0\\n  1  6  1  0\\n  3  4  1  0\\n  3  7  1  0\\n  3  8  1  0\\n  4  5  1  0\\n  4  9  1  0\\n  4 10  1  0\\n  5 11  1  0\\n  5 12  1  0\\n  5 13  1  0\\nM  RAD  1   3   2\\nM  END\\n\",\"sdf\");\n\tviewer_16580360278336098.setStyle({\"stick\": {}});\n\tviewer_16580360278336098.zoomTo();\nviewer_16580360278336098.render();\n});\n</script>",
-      "text/html": [
-       "<div id=\"3dmolviewer_16580360278336098\"  style=\"position: relative; width: 400px; height: 400px\">\n",
-       "        <p id=\"3dmolwarning_16580360278336098\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n",
-       "        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n",
-       "        </div>\n",
-       "<script>\n",
-       "\n",
-       "var loadScriptAsync = function(uri){\n",
-       "  return new Promise((resolve, reject) => {\n",
-       "    var tag = document.createElement('script');\n",
-       "    tag.src = uri;\n",
-       "    tag.async = true;\n",
-       "    tag.onload = () => {\n",
-       "      resolve();\n",
-       "    };\n",
-       "  var firstScriptTag = document.getElementsByTagName('script')[0];\n",
-       "  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n",
-       "});\n",
-       "};\n",
-       "\n",
-       "if(typeof $3Dmolpromise === 'undefined') {\n",
-       "$3Dmolpromise = null;\n",
-       "  $3Dmolpromise = loadScriptAsync('https://3dmol.org/build/3Dmol.js');\n",
-       "}\n",
-       "\n",
-       "var viewer_16580360278336098 = null;\n",
-       "var warn = document.getElementById(\"3dmolwarning_16580360278336098\");\n",
-       "if(warn) {\n",
-       "    warn.parentNode.removeChild(warn);\n",
-       "}\n",
-       "$3Dmolpromise.then(function() {\n",
-       "viewer_16580360278336098 = $3Dmol.createViewer($(\"#3dmolviewer_16580360278336098\"),{backgroundColor:\"white\"});\n",
-       "viewer_16580360278336098.zoomTo();\n",
-       "\tviewer_16580360278336098.addModel(\"\\n     RDKit          3D\\n\\n 13 11  0  0  0  0  0  0  0  0999 V2000\\n    1.4470    1.5562    0.1015 O   0  0  0  0  0  0  0  0  0  0  0  0\\n    1.3442    2.4284    0.5112 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5601    0.3078    0.8404 C   0  0  0  0  0  3  0  0  0  0  0  0\\n   -0.5601   -1.1042    0.2796 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5601   -1.1042   -1.2397 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.5046    0.9701    0.4485 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5601    0.2814    1.9345 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4473    0.8597    0.5138 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4425   -1.6415    0.6444 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3223   -1.6415    0.6444 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3272   -0.5970   -1.6320 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4473   -0.5970   -1.6320 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5601   -2.1309   -1.6190 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1  2  1  0\\n  1  6  1  0\\n  3  4  1  0\\n  3  7  1  0\\n  3  8  1  0\\n  4  5  1  0\\n  4  9  1  0\\n  4 10  1  0\\n  5 11  1  0\\n  5 12  1  0\\n  5 13  1  0\\nM  RAD  1   3   2\\nM  END\\n\",\"sdf\");\n",
-       "\tviewer_16580360278336098.setStyle({\"stick\": {}});\n",
-       "\tviewer_16580360278336098.zoomTo();\n",
-       "viewer_16580360278336098.render();\n",
-       "});\n",
-       "</script>"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    },
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "\n",
-      "\n",
-      "\n",
-      "\n",
-      "heuristics 3 (success: True):\n",
-      "coordinates:\n",
-      "O       1.43450155    1.57365164    0.09693443\n",
-      "H       2.09087982    1.39035494    0.78579648\n",
-      "C      -0.57261278    0.32526157    0.83577095\n",
-      "C      -0.57261278   -1.08680339    0.27502014\n",
-      "C      -0.57261278   -1.08680339   -1.24431136\n",
-      "H       0.49207957    0.98748162    0.44384829\n",
-      "H      -0.57261353    0.29885538    1.92993720\n",
-      "H      -1.45985625    0.87711180    0.50916755\n",
-      "H      -1.45502736   -1.62411740    0.63977604\n",
-      "H       0.30980206   -1.62411723    0.63977598\n",
-      "H       0.31463056   -0.57958326   -1.63665841\n",
-      "H      -1.45985628   -0.57958358   -1.63665843\n",
-      "H      -0.57261248   -2.11346595   -1.62360141\n"
-     ]
-    },
-    {
-     "data": {
-      "application/3dmoljs_load.v0": "<div id=\"3dmolviewer_16580360278368294\"  style=\"position: relative; width: 400px; height: 400px\">\n        <p id=\"3dmolwarning_16580360278368294\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n        </div>\n<script>\n\nvar loadScriptAsync = function(uri){\n  return new Promise((resolve, reject) => {\n    var tag = document.createElement('script');\n    tag.src = uri;\n    tag.async = true;\n    tag.onload = () => {\n      resolve();\n    };\n  var firstScriptTag = document.getElementsByTagName('script')[0];\n  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n});\n};\n\nif(typeof $3Dmolpromise === 'undefined') {\n$3Dmolpromise = null;\n  $3Dmolpromise = loadScriptAsync('https://3dmol.org/build/3Dmol.js');\n}\n\nvar viewer_16580360278368294 = null;\nvar warn = document.getElementById(\"3dmolwarning_16580360278368294\");\nif(warn) {\n    warn.parentNode.removeChild(warn);\n}\n$3Dmolpromise.then(function() {\nviewer_16580360278368294 = $3Dmol.createViewer($(\"#3dmolviewer_16580360278368294\"),{backgroundColor:\"white\"});\nviewer_16580360278368294.zoomTo();\n\tviewer_16580360278368294.addModel(\"\\n     RDKit          3D\\n\\n 13 11  0  0  0  0  0  0  0  0999 V2000\\n    1.4345    1.5737    0.0969 O   0  0  0  0  0  0  0  0  0  0  0  0\\n    2.0909    1.3904    0.7858 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5726    0.3253    0.8358 C   0  0  0  0  0  3  0  0  0  0  0  0\\n   -0.5726   -1.0868    0.2750 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5726   -1.0868   -1.2443 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.4921    0.9875    0.4438 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5726    0.2989    1.9299 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4599    0.8771    0.5092 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4550   -1.6241    0.6398 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3098   -1.6241    0.6398 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3146   -0.5796   -1.6367 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4599   -0.5796   -1.6367 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5726   -2.1135   -1.6236 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1  2  1  0\\n  1  6  1  0\\n  3  4  1  0\\n  3  7  1  0\\n  3  8  1  0\\n  4  5  1  0\\n  4  9  1  0\\n  4 10  1  0\\n  5 11  1  0\\n  5 12  1  0\\n  5 13  1  0\\nM  RAD  1   3   2\\nM  END\\n\",\"sdf\");\n\tviewer_16580360278368294.setStyle({\"stick\": {}});\n\tviewer_16580360278368294.zoomTo();\nviewer_16580360278368294.render();\n});\n</script>",
-      "text/html": [
-       "<div id=\"3dmolviewer_16580360278368294\"  style=\"position: relative; width: 400px; height: 400px\">\n",
-       "        <p id=\"3dmolwarning_16580360278368294\" style=\"background-color:#ffcccc;color:black\">You appear to be running in JupyterLab (or JavaScript failed to load for some other reason).  You need to install the 3dmol extension: <br>\n",
-       "        <tt>jupyter labextension install jupyterlab_3dmol</tt></p>\n",
-       "        </div>\n",
-       "<script>\n",
-       "\n",
-       "var loadScriptAsync = function(uri){\n",
-       "  return new Promise((resolve, reject) => {\n",
-       "    var tag = document.createElement('script');\n",
-       "    tag.src = uri;\n",
-       "    tag.async = true;\n",
-       "    tag.onload = () => {\n",
-       "      resolve();\n",
-       "    };\n",
-       "  var firstScriptTag = document.getElementsByTagName('script')[0];\n",
-       "  firstScriptTag.parentNode.insertBefore(tag, firstScriptTag);\n",
-       "});\n",
-       "};\n",
-       "\n",
-       "if(typeof $3Dmolpromise === 'undefined') {\n",
-       "$3Dmolpromise = null;\n",
-       "  $3Dmolpromise = loadScriptAsync('https://3dmol.org/build/3Dmol.js');\n",
-       "}\n",
-       "\n",
-       "var viewer_16580360278368294 = null;\n",
-       "var warn = document.getElementById(\"3dmolwarning_16580360278368294\");\n",
-       "if(warn) {\n",
-       "    warn.parentNode.removeChild(warn);\n",
-       "}\n",
-       "$3Dmolpromise.then(function() {\n",
-       "viewer_16580360278368294 = $3Dmol.createViewer($(\"#3dmolviewer_16580360278368294\"),{backgroundColor:\"white\"});\n",
-       "viewer_16580360278368294.zoomTo();\n",
-       "\tviewer_16580360278368294.addModel(\"\\n     RDKit          3D\\n\\n 13 11  0  0  0  0  0  0  0  0999 V2000\\n    1.4345    1.5737    0.0969 O   0  0  0  0  0  0  0  0  0  0  0  0\\n    2.0909    1.3904    0.7858 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5726    0.3253    0.8358 C   0  0  0  0  0  3  0  0  0  0  0  0\\n   -0.5726   -1.0868    0.2750 C   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5726   -1.0868   -1.2443 C   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.4921    0.9875    0.4438 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5726    0.2989    1.9299 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4599    0.8771    0.5092 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4550   -1.6241    0.6398 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3098   -1.6241    0.6398 H   0  0  0  0  0  0  0  0  0  0  0  0\\n    0.3146   -0.5796   -1.6367 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -1.4599   -0.5796   -1.6367 H   0  0  0  0  0  0  0  0  0  0  0  0\\n   -0.5726   -2.1135   -1.6236 H   0  0  0  0  0  0  0  0  0  0  0  0\\n  1  2  1  0\\n  1  6  1  0\\n  3  4  1  0\\n  3  7  1  0\\n  3  8  1  0\\n  4  5  1  0\\n  4  9  1  0\\n  4 10  1  0\\n  5 11  1  0\\n  5 12  1  0\\n  5 13  1  0\\nM  RAD  1   3   2\\nM  END\\n\",\"sdf\");\n",
-       "\tviewer_16580360278368294.setStyle({\"stick\": {}});\n",
-       "\tviewer_16580360278368294.zoomTo();\n",
-       "viewer_16580360278368294.render();\n",
-       "});\n",
-       "</script>"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    },
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "\n",
-      "\n",
-      "\n",
-      "\n"
-     ]
-    }
-   ],
+   "outputs": [],
    "source": [
     "for tsg in rxn.ts_species.ts_guesses:\n",
     "    print(f'{tsg.method} {tsg.method_index} (success: {tsg.success}):')\n",
diff --git a/ipython/Tools/Visualize 1D Torsion scan.ipynb b/ipython/Tools/Visualize 1D Torsion scan.ipynb
index 71c7a94058..e393c7ade2 100644
--- a/ipython/Tools/Visualize 1D Torsion scan.ipynb	
+++ b/ipython/Tools/Visualize 1D Torsion scan.ipynb	
@@ -17,7 +17,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -30,847 +30,17 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stderr",
-     "output_type": "stream",
-     "text": [
-      "WARNING:root:No kyotocabinet available\n"
-     ]
-    }
-   ],
-   "source": [
-    "from arc.species.species import determine_rotor_symmetry\n",
-    "from arc.parser import parse_1d_scan_energies\n",
-    "from arc.common import read_yaml_file\n",
-    "from arc.plotter import plot_1d_rotor_scan\n",
-    "%matplotlib notebook"
-   ]
+   "outputs": [],
+   "source": "from arc.species.species import determine_rotor_symmetry\nfrom arc.parser.parser import parse_1d_scan_energies\nfrom arc.common import read_yaml_file\nfrom arc.plotter import plot_1d_rotor_scan\n%matplotlib notebook"
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "data": {
-      "application/javascript": [
-       "/* Put everything inside the global mpl namespace */\n",
-       "window.mpl = {};\n",
-       "\n",
-       "\n",
-       "mpl.get_websocket_type = function() {\n",
-       "    if (typeof(WebSocket) !== 'undefined') {\n",
-       "        return WebSocket;\n",
-       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
-       "        return MozWebSocket;\n",
-       "    } else {\n",
-       "        alert('Your browser does not have WebSocket support. ' +\n",
-       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
-       "              'Firefox 4 and 5 are also supported but you ' +\n",
-       "              'have to enable WebSockets in about:config.');\n",
-       "    };\n",
-       "}\n",
-       "\n",
-       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
-       "    this.id = figure_id;\n",
-       "\n",
-       "    this.ws = websocket;\n",
-       "\n",
-       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
-       "\n",
-       "    if (!this.supports_binary) {\n",
-       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
-       "        if (warnings) {\n",
-       "            warnings.style.display = 'block';\n",
-       "            warnings.textContent = (\n",
-       "                \"This browser does not support binary websocket messages. \" +\n",
-       "                    \"Performance may be slow.\");\n",
-       "        }\n",
-       "    }\n",
-       "\n",
-       "    this.imageObj = new Image();\n",
-       "\n",
-       "    this.context = undefined;\n",
-       "    this.message = undefined;\n",
-       "    this.canvas = undefined;\n",
-       "    this.rubberband_canvas = undefined;\n",
-       "    this.rubberband_context = undefined;\n",
-       "    this.format_dropdown = undefined;\n",
-       "\n",
-       "    this.image_mode = 'full';\n",
-       "\n",
-       "    this.root = $('<div/>');\n",
-       "    this._root_extra_style(this.root)\n",
-       "    this.root.attr('style', 'display: inline-block');\n",
-       "\n",
-       "    $(parent_element).append(this.root);\n",
-       "\n",
-       "    this._init_header(this);\n",
-       "    this._init_canvas(this);\n",
-       "    this._init_toolbar(this);\n",
-       "\n",
-       "    var fig = this;\n",
-       "\n",
-       "    this.waiting = false;\n",
-       "\n",
-       "    this.ws.onopen =  function () {\n",
-       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
-       "            fig.send_message(\"send_image_mode\", {});\n",
-       "            if (mpl.ratio != 1) {\n",
-       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
-       "            }\n",
-       "            fig.send_message(\"refresh\", {});\n",
-       "        }\n",
-       "\n",
-       "    this.imageObj.onload = function() {\n",
-       "            if (fig.image_mode == 'full') {\n",
-       "                // Full images could contain transparency (where diff images\n",
-       "                // almost always do), so we need to clear the canvas so that\n",
-       "                // there is no ghosting.\n",
-       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
-       "            }\n",
-       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
-       "        };\n",
-       "\n",
-       "    this.imageObj.onunload = function() {\n",
-       "        fig.ws.close();\n",
-       "    }\n",
-       "\n",
-       "    this.ws.onmessage = this._make_on_message_function(this);\n",
-       "\n",
-       "    this.ondownload = ondownload;\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype._init_header = function() {\n",
-       "    var titlebar = $(\n",
-       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
-       "        'ui-helper-clearfix\"/>');\n",
-       "    var titletext = $(\n",
-       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
-       "        'text-align: center; padding: 3px;\"/>');\n",
-       "    titlebar.append(titletext)\n",
-       "    this.root.append(titlebar);\n",
-       "    this.header = titletext[0];\n",
-       "}\n",
-       "\n",
-       "\n",
-       "\n",
-       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
-       "\n",
-       "}\n",
-       "\n",
-       "\n",
-       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
-       "\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype._init_canvas = function() {\n",
-       "    var fig = this;\n",
-       "\n",
-       "    var canvas_div = $('<div/>');\n",
-       "\n",
-       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
-       "\n",
-       "    function canvas_keyboard_event(event) {\n",
-       "        return fig.key_event(event, event['data']);\n",
-       "    }\n",
-       "\n",
-       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
-       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
-       "    this.canvas_div = canvas_div\n",
-       "    this._canvas_extra_style(canvas_div)\n",
-       "    this.root.append(canvas_div);\n",
-       "\n",
-       "    var canvas = $('<canvas/>');\n",
-       "    canvas.addClass('mpl-canvas');\n",
-       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
-       "\n",
-       "    this.canvas = canvas[0];\n",
-       "    this.context = canvas[0].getContext(\"2d\");\n",
-       "\n",
-       "    var backingStore = this.context.backingStorePixelRatio ||\n",
-       "\tthis.context.webkitBackingStorePixelRatio ||\n",
-       "\tthis.context.mozBackingStorePixelRatio ||\n",
-       "\tthis.context.msBackingStorePixelRatio ||\n",
-       "\tthis.context.oBackingStorePixelRatio ||\n",
-       "\tthis.context.backingStorePixelRatio || 1;\n",
-       "\n",
-       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
-       "\n",
-       "    var rubberband = $('<canvas/>');\n",
-       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
-       "\n",
-       "    var pass_mouse_events = true;\n",
-       "\n",
-       "    canvas_div.resizable({\n",
-       "        start: function(event, ui) {\n",
-       "            pass_mouse_events = false;\n",
-       "        },\n",
-       "        resize: function(event, ui) {\n",
-       "            fig.request_resize(ui.size.width, ui.size.height);\n",
-       "        },\n",
-       "        stop: function(event, ui) {\n",
-       "            pass_mouse_events = true;\n",
-       "            fig.request_resize(ui.size.width, ui.size.height);\n",
-       "        },\n",
-       "    });\n",
-       "\n",
-       "    function mouse_event_fn(event) {\n",
-       "        if (pass_mouse_events)\n",
-       "            return fig.mouse_event(event, event['data']);\n",
-       "    }\n",
-       "\n",
-       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
-       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
-       "    // Throttle sequential mouse events to 1 every 20ms.\n",
-       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
-       "\n",
-       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
-       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
-       "\n",
-       "    canvas_div.on(\"wheel\", function (event) {\n",
-       "        event = event.originalEvent;\n",
-       "        event['data'] = 'scroll'\n",
-       "        if (event.deltaY < 0) {\n",
-       "            event.step = 1;\n",
-       "        } else {\n",
-       "            event.step = -1;\n",
-       "        }\n",
-       "        mouse_event_fn(event);\n",
-       "    });\n",
-       "\n",
-       "    canvas_div.append(canvas);\n",
-       "    canvas_div.append(rubberband);\n",
-       "\n",
-       "    this.rubberband = rubberband;\n",
-       "    this.rubberband_canvas = rubberband[0];\n",
-       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
-       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
-       "\n",
-       "    this._resize_canvas = function(width, height) {\n",
-       "        // Keep the size of the canvas, canvas container, and rubber band\n",
-       "        // canvas in synch.\n",
-       "        canvas_div.css('width', width)\n",
-       "        canvas_div.css('height', height)\n",
-       "\n",
-       "        canvas.attr('width', width * mpl.ratio);\n",
-       "        canvas.attr('height', height * mpl.ratio);\n",
-       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
-       "\n",
-       "        rubberband.attr('width', width);\n",
-       "        rubberband.attr('height', height);\n",
-       "    }\n",
-       "\n",
-       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
-       "    // upon first draw.\n",
-       "    this._resize_canvas(600, 600);\n",
-       "\n",
-       "    // Disable right mouse context menu.\n",
-       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
-       "        return false;\n",
-       "    });\n",
-       "\n",
-       "    function set_focus () {\n",
-       "        canvas.focus();\n",
-       "        canvas_div.focus();\n",
-       "    }\n",
-       "\n",
-       "    window.setTimeout(set_focus, 100);\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype._init_toolbar = function() {\n",
-       "    var fig = this;\n",
-       "\n",
-       "    var nav_element = $('<div/>');\n",
-       "    nav_element.attr('style', 'width: 100%');\n",
-       "    this.root.append(nav_element);\n",
-       "\n",
-       "    // Define a callback function for later on.\n",
-       "    function toolbar_event(event) {\n",
-       "        return fig.toolbar_button_onclick(event['data']);\n",
-       "    }\n",
-       "    function toolbar_mouse_event(event) {\n",
-       "        return fig.toolbar_button_onmouseover(event['data']);\n",
-       "    }\n",
-       "\n",
-       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
-       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
-       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
-       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
-       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
-       "\n",
-       "        if (!name) {\n",
-       "            // put a spacer in here.\n",
-       "            continue;\n",
-       "        }\n",
-       "        var button = $('<button/>');\n",
-       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
-       "                        'ui-button-icon-only');\n",
-       "        button.attr('role', 'button');\n",
-       "        button.attr('aria-disabled', 'false');\n",
-       "        button.click(method_name, toolbar_event);\n",
-       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
-       "\n",
-       "        var icon_img = $('<span/>');\n",
-       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
-       "        icon_img.addClass(image);\n",
-       "        icon_img.addClass('ui-corner-all');\n",
-       "\n",
-       "        var tooltip_span = $('<span/>');\n",
-       "        tooltip_span.addClass('ui-button-text');\n",
-       "        tooltip_span.html(tooltip);\n",
-       "\n",
-       "        button.append(icon_img);\n",
-       "        button.append(tooltip_span);\n",
-       "\n",
-       "        nav_element.append(button);\n",
-       "    }\n",
-       "\n",
-       "    var fmt_picker_span = $('<span/>');\n",
-       "\n",
-       "    var fmt_picker = $('<select/>');\n",
-       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
-       "    fmt_picker_span.append(fmt_picker);\n",
-       "    nav_element.append(fmt_picker_span);\n",
-       "    this.format_dropdown = fmt_picker[0];\n",
-       "\n",
-       "    for (var ind in mpl.extensions) {\n",
-       "        var fmt = mpl.extensions[ind];\n",
-       "        var option = $(\n",
-       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
-       "        fmt_picker.append(option);\n",
-       "    }\n",
-       "\n",
-       "    // Add hover states to the ui-buttons\n",
-       "    $( \".ui-button\" ).hover(\n",
-       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
-       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
-       "    );\n",
-       "\n",
-       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
-       "    nav_element.append(status_bar);\n",
-       "    this.message = status_bar[0];\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
-       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
-       "    // which will in turn request a refresh of the image.\n",
-       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.send_message = function(type, properties) {\n",
-       "    properties['type'] = type;\n",
-       "    properties['figure_id'] = this.id;\n",
-       "    this.ws.send(JSON.stringify(properties));\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.send_draw_message = function() {\n",
-       "    if (!this.waiting) {\n",
-       "        this.waiting = true;\n",
-       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
-       "    }\n",
-       "}\n",
-       "\n",
-       "\n",
-       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
-       "    var format_dropdown = fig.format_dropdown;\n",
-       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
-       "    fig.ondownload(fig, format);\n",
-       "}\n",
-       "\n",
-       "\n",
-       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
-       "    var size = msg['size'];\n",
-       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
-       "        fig._resize_canvas(size[0], size[1]);\n",
-       "        fig.send_message(\"refresh\", {});\n",
-       "    };\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
-       "    var x0 = msg['x0'] / mpl.ratio;\n",
-       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
-       "    var x1 = msg['x1'] / mpl.ratio;\n",
-       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
-       "    x0 = Math.floor(x0) + 0.5;\n",
-       "    y0 = Math.floor(y0) + 0.5;\n",
-       "    x1 = Math.floor(x1) + 0.5;\n",
-       "    y1 = Math.floor(y1) + 0.5;\n",
-       "    var min_x = Math.min(x0, x1);\n",
-       "    var min_y = Math.min(y0, y1);\n",
-       "    var width = Math.abs(x1 - x0);\n",
-       "    var height = Math.abs(y1 - y0);\n",
-       "\n",
-       "    fig.rubberband_context.clearRect(\n",
-       "        0, 0, fig.canvas.width / mpl.ratio, fig.canvas.height / mpl.ratio);\n",
-       "\n",
-       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
-       "    // Updates the figure title.\n",
-       "    fig.header.textContent = msg['label'];\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
-       "    var cursor = msg['cursor'];\n",
-       "    switch(cursor)\n",
-       "    {\n",
-       "    case 0:\n",
-       "        cursor = 'pointer';\n",
-       "        break;\n",
-       "    case 1:\n",
-       "        cursor = 'default';\n",
-       "        break;\n",
-       "    case 2:\n",
-       "        cursor = 'crosshair';\n",
-       "        break;\n",
-       "    case 3:\n",
-       "        cursor = 'move';\n",
-       "        break;\n",
-       "    }\n",
-       "    fig.rubberband_canvas.style.cursor = cursor;\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
-       "    fig.message.textContent = msg['message'];\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
-       "    // Request the server to send over a new figure.\n",
-       "    fig.send_draw_message();\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
-       "    fig.image_mode = msg['mode'];\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.updated_canvas_event = function() {\n",
-       "    // Called whenever the canvas gets updated.\n",
-       "    this.send_message(\"ack\", {});\n",
-       "}\n",
-       "\n",
-       "// A function to construct a web socket function for onmessage handling.\n",
-       "// Called in the figure constructor.\n",
-       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
-       "    return function socket_on_message(evt) {\n",
-       "        if (evt.data instanceof Blob) {\n",
-       "            /* FIXME: We get \"Resource interpreted as Image but\n",
-       "             * transferred with MIME type text/plain:\" errors on\n",
-       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
-       "             * to be part of the websocket stream */\n",
-       "            evt.data.type = \"image/png\";\n",
-       "\n",
-       "            /* Free the memory for the previous frames */\n",
-       "            if (fig.imageObj.src) {\n",
-       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
-       "                    fig.imageObj.src);\n",
-       "            }\n",
-       "\n",
-       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
-       "                evt.data);\n",
-       "            fig.updated_canvas_event();\n",
-       "            fig.waiting = false;\n",
-       "            return;\n",
-       "        }\n",
-       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
-       "            fig.imageObj.src = evt.data;\n",
-       "            fig.updated_canvas_event();\n",
-       "            fig.waiting = false;\n",
-       "            return;\n",
-       "        }\n",
-       "\n",
-       "        var msg = JSON.parse(evt.data);\n",
-       "        var msg_type = msg['type'];\n",
-       "\n",
-       "        // Call the  \"handle_{type}\" callback, which takes\n",
-       "        // the figure and JSON message as its only arguments.\n",
-       "        try {\n",
-       "            var callback = fig[\"handle_\" + msg_type];\n",
-       "        } catch (e) {\n",
-       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
-       "            return;\n",
-       "        }\n",
-       "\n",
-       "        if (callback) {\n",
-       "            try {\n",
-       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
-       "                callback(fig, msg);\n",
-       "            } catch (e) {\n",
-       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
-       "            }\n",
-       "        }\n",
-       "    };\n",
-       "}\n",
-       "\n",
-       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
-       "mpl.findpos = function(e) {\n",
-       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
-       "    var targ;\n",
-       "    if (!e)\n",
-       "        e = window.event;\n",
-       "    if (e.target)\n",
-       "        targ = e.target;\n",
-       "    else if (e.srcElement)\n",
-       "        targ = e.srcElement;\n",
-       "    if (targ.nodeType == 3) // defeat Safari bug\n",
-       "        targ = targ.parentNode;\n",
-       "\n",
-       "    // jQuery normalizes the pageX and pageY\n",
-       "    // pageX,Y are the mouse positions relative to the document\n",
-       "    // offset() returns the position of the element relative to the document\n",
-       "    var x = e.pageX - $(targ).offset().left;\n",
-       "    var y = e.pageY - $(targ).offset().top;\n",
-       "\n",
-       "    return {\"x\": x, \"y\": y};\n",
-       "};\n",
-       "\n",
-       "/*\n",
-       " * return a copy of an object with only non-object keys\n",
-       " * we need this to avoid circular references\n",
-       " * http://stackoverflow.com/a/24161582/3208463\n",
-       " */\n",
-       "function simpleKeys (original) {\n",
-       "  return Object.keys(original).reduce(function (obj, key) {\n",
-       "    if (typeof original[key] !== 'object')\n",
-       "        obj[key] = original[key]\n",
-       "    return obj;\n",
-       "  }, {});\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
-       "    var canvas_pos = mpl.findpos(event)\n",
-       "\n",
-       "    if (name === 'button_press')\n",
-       "    {\n",
-       "        this.canvas.focus();\n",
-       "        this.canvas_div.focus();\n",
-       "    }\n",
-       "\n",
-       "    var x = canvas_pos.x * mpl.ratio;\n",
-       "    var y = canvas_pos.y * mpl.ratio;\n",
-       "\n",
-       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
-       "                             step: event.step,\n",
-       "                             guiEvent: simpleKeys(event)});\n",
-       "\n",
-       "    /* This prevents the web browser from automatically changing to\n",
-       "     * the text insertion cursor when the button is pressed.  We want\n",
-       "     * to control all of the cursor setting manually through the\n",
-       "     * 'cursor' event from matplotlib */\n",
-       "    event.preventDefault();\n",
-       "    return false;\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
-       "    // Handle any extra behaviour associated with a key event\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.key_event = function(event, name) {\n",
-       "\n",
-       "    // Prevent repeat events\n",
-       "    if (name == 'key_press')\n",
-       "    {\n",
-       "        if (event.which === this._key)\n",
-       "            return;\n",
-       "        else\n",
-       "            this._key = event.which;\n",
-       "    }\n",
-       "    if (name == 'key_release')\n",
-       "        this._key = null;\n",
-       "\n",
-       "    var value = '';\n",
-       "    if (event.ctrlKey && event.which != 17)\n",
-       "        value += \"ctrl+\";\n",
-       "    if (event.altKey && event.which != 18)\n",
-       "        value += \"alt+\";\n",
-       "    if (event.shiftKey && event.which != 16)\n",
-       "        value += \"shift+\";\n",
-       "\n",
-       "    value += 'k';\n",
-       "    value += event.which.toString();\n",
-       "\n",
-       "    this._key_event_extra(event, name);\n",
-       "\n",
-       "    this.send_message(name, {key: value,\n",
-       "                             guiEvent: simpleKeys(event)});\n",
-       "    return false;\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
-       "    if (name == 'download') {\n",
-       "        this.handle_save(this, null);\n",
-       "    } else {\n",
-       "        this.send_message(\"toolbar_button\", {name: name});\n",
-       "    }\n",
-       "};\n",
-       "\n",
-       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
-       "    this.message.textContent = tooltip;\n",
-       "};\n",
-       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
-       "\n",
-       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
-       "\n",
-       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
-       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
-       "    // object with the appropriate methods. Currently this is a non binary\n",
-       "    // socket, so there is still some room for performance tuning.\n",
-       "    var ws = {};\n",
-       "\n",
-       "    ws.close = function() {\n",
-       "        comm.close()\n",
-       "    };\n",
-       "    ws.send = function(m) {\n",
-       "        //console.log('sending', m);\n",
-       "        comm.send(m);\n",
-       "    };\n",
-       "    // Register the callback with on_msg.\n",
-       "    comm.on_msg(function(msg) {\n",
-       "        //console.log('receiving', msg['content']['data'], msg);\n",
-       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
-       "        ws.onmessage(msg['content']['data'])\n",
-       "    });\n",
-       "    return ws;\n",
-       "}\n",
-       "\n",
-       "mpl.mpl_figure_comm = function(comm, msg) {\n",
-       "    // This is the function which gets called when the mpl process\n",
-       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
-       "\n",
-       "    var id = msg.content.data.id;\n",
-       "    // Get hold of the div created by the display call when the Comm\n",
-       "    // socket was opened in Python.\n",
-       "    var element = $(\"#\" + id);\n",
-       "    var ws_proxy = comm_websocket_adapter(comm)\n",
-       "\n",
-       "    function ondownload(figure, format) {\n",
-       "        window.open(figure.imageObj.src);\n",
-       "    }\n",
-       "\n",
-       "    var fig = new mpl.figure(id, ws_proxy,\n",
-       "                           ondownload,\n",
-       "                           element.get(0));\n",
-       "\n",
-       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
-       "    // web socket which is closed, not our websocket->open comm proxy.\n",
-       "    ws_proxy.onopen();\n",
-       "\n",
-       "    fig.parent_element = element.get(0);\n",
-       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
-       "    if (!fig.cell_info) {\n",
-       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
-       "        return;\n",
-       "    }\n",
-       "\n",
-       "    var output_index = fig.cell_info[2]\n",
-       "    var cell = fig.cell_info[0];\n",
-       "\n",
-       "};\n",
-       "\n",
-       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
-       "    var width = fig.canvas.width/mpl.ratio\n",
-       "    fig.root.unbind('remove')\n",
-       "\n",
-       "    // Update the output cell to use the data from the current canvas.\n",
-       "    fig.push_to_output();\n",
-       "    var dataURL = fig.canvas.toDataURL();\n",
-       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
-       "    // the notebook keyboard shortcuts fail.\n",
-       "    IPython.keyboard_manager.enable()\n",
-       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
-       "    fig.close_ws(fig, msg);\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
-       "    fig.send_message('closing', msg);\n",
-       "    // fig.ws.close()\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
-       "    // Turn the data on the canvas into data in the output cell.\n",
-       "    var width = this.canvas.width/mpl.ratio\n",
-       "    var dataURL = this.canvas.toDataURL();\n",
-       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.updated_canvas_event = function() {\n",
-       "    // Tell IPython that the notebook contents must change.\n",
-       "    IPython.notebook.set_dirty(true);\n",
-       "    this.send_message(\"ack\", {});\n",
-       "    var fig = this;\n",
-       "    // Wait a second, then push the new image to the DOM so\n",
-       "    // that it is saved nicely (might be nice to debounce this).\n",
-       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype._init_toolbar = function() {\n",
-       "    var fig = this;\n",
-       "\n",
-       "    var nav_element = $('<div/>');\n",
-       "    nav_element.attr('style', 'width: 100%');\n",
-       "    this.root.append(nav_element);\n",
-       "\n",
-       "    // Define a callback function for later on.\n",
-       "    function toolbar_event(event) {\n",
-       "        return fig.toolbar_button_onclick(event['data']);\n",
-       "    }\n",
-       "    function toolbar_mouse_event(event) {\n",
-       "        return fig.toolbar_button_onmouseover(event['data']);\n",
-       "    }\n",
-       "\n",
-       "    for(var toolbar_ind in mpl.toolbar_items){\n",
-       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
-       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
-       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
-       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
-       "\n",
-       "        if (!name) { continue; };\n",
-       "\n",
-       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
-       "        button.click(method_name, toolbar_event);\n",
-       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
-       "        nav_element.append(button);\n",
-       "    }\n",
-       "\n",
-       "    // Add the status bar.\n",
-       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
-       "    nav_element.append(status_bar);\n",
-       "    this.message = status_bar[0];\n",
-       "\n",
-       "    // Add the close button to the window.\n",
-       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
-       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
-       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
-       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
-       "    buttongrp.append(button);\n",
-       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
-       "    titlebar.prepend(buttongrp);\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype._root_extra_style = function(el){\n",
-       "    var fig = this\n",
-       "    el.on(\"remove\", function(){\n",
-       "\tfig.close_ws(fig, {});\n",
-       "    });\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
-       "    // this is important to make the div 'focusable\n",
-       "    el.attr('tabindex', 0)\n",
-       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
-       "    // off when our div gets focus\n",
-       "\n",
-       "    // location in version 3\n",
-       "    if (IPython.notebook.keyboard_manager) {\n",
-       "        IPython.notebook.keyboard_manager.register_events(el);\n",
-       "    }\n",
-       "    else {\n",
-       "        // location in version 2\n",
-       "        IPython.keyboard_manager.register_events(el);\n",
-       "    }\n",
-       "\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
-       "    var manager = IPython.notebook.keyboard_manager;\n",
-       "    if (!manager)\n",
-       "        manager = IPython.keyboard_manager;\n",
-       "\n",
-       "    // Check for shift+enter\n",
-       "    if (event.shiftKey && event.which == 13) {\n",
-       "        this.canvas_div.blur();\n",
-       "        event.shiftKey = false;\n",
-       "        // Send a \"J\" for go to next cell\n",
-       "        event.which = 74;\n",
-       "        event.keyCode = 74;\n",
-       "        manager.command_mode();\n",
-       "        manager.handle_keydown(event);\n",
-       "    }\n",
-       "}\n",
-       "\n",
-       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
-       "    fig.ondownload(fig, null);\n",
-       "}\n",
-       "\n",
-       "\n",
-       "mpl.find_output_cell = function(html_output) {\n",
-       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
-       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
-       "    // IPython event is triggered only after the cells have been serialised, which for\n",
-       "    // our purposes (turning an active figure into a static one), is too late.\n",
-       "    var cells = IPython.notebook.get_cells();\n",
-       "    var ncells = cells.length;\n",
-       "    for (var i=0; i<ncells; i++) {\n",
-       "        var cell = cells[i];\n",
-       "        if (cell.cell_type === 'code'){\n",
-       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
-       "                var data = cell.output_area.outputs[j];\n",
-       "                if (data.data) {\n",
-       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
-       "                    data = data.data;\n",
-       "                }\n",
-       "                if (data['text/html'] == html_output) {\n",
-       "                    return [cell, data, j];\n",
-       "                }\n",
-       "            }\n",
-       "        }\n",
-       "    }\n",
-       "}\n",
-       "\n",
-       "// Register the function which deals with the matplotlib target/channel.\n",
-       "// The kernel may be null if the page has been refreshed.\n",
-       "if (IPython.notebook.kernel != null) {\n",
-       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
-       "}\n"
-      ],
-      "text/plain": [
-       "<IPython.core.display.Javascript object>"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    },
-    {
-     "data": {
-      "text/html": [
-       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAeAAAAFoCAYAAACPNyggAAAgAElEQVR4nOzdeXgT1foH8EkopYVCQRG06E9QCUIrooKiXEAuFhAKwhW4oKLobWm5COWyg2hxoVUQXFAUUEQBi1SwiqAsactS1padlqVluu97uieZ7++PtmNjUugyyZlJ3s/zzONjMpnzThryzcycOYcDIYQQQmyOY10AIYQQ4ogogAkhhBAGKIAJIYQQBiiACSGEEAYogAkhhBAGKIAJIYQQBiiACSGEEAYogAkhhBAGKIAJIYQQBiiACSGEEAYogAkhhBAGKIAJIYQQBiiACSGEEAYogAkhhBAGKIAJIYQQBiiACSGEEAYogG+jsLAQ4eHhiI2NxeXLl2mhhRZaaHHgJTY2FuHh4SgsLGxxvlAA30Z4eDg4jqOFFlpooYUWcQkPD29xvig2gDdt2gSO49CuXTuz52JjYzF8+HC0a9cO7u7umDBhAhITE5vVTmxsrPhms/7lRQsttNBCC9ul7qAsNja2pTGmzABOS0uDu7s7PDw8zAI4Pj4e7du3x+DBg7F3717s2rULnp6e8PDwQE5OTpPbunz5MjiOw+XLl6UqnxBCiEJJmQmKDGAfHx+MHTsWr732mlkAT5o0CZ07d0ZxcbH4WFJSElq3bo1FixY1uS0KYEIIIXUcOoC3bt2K9u3bIzU11SyA9Xo9XF1d4e/vb/a6ESNGoGfPnk1ujwKYEEJIHYcN4OzsbNx555348ssvAcAsgK9evQqO48Tn61uwYAFUKhUqKiqa1CYFMCGEkDoOG8AvvvginnnmGQiCAMA8gKOjo8FxHEJDQ81eGxwcDI7jkJGR0eD2s7OzG7zgTgFMCCHEIQP4559/hrOzM65cuSI+1lAA79ixw+z1dQGcmZnZYBtBQUENdjmnACaEEOJwAazT6dC1a1fMnz8fhYWF4jJ16lS0a9cOhYWFKC0tbfEpaDoCJoQQcisOF8A8z9/2pugXXnhB7IQVEBBgto2RI0dSJyxCCCEt4nABXFFRgcjISLNl5MiRcHFxQWRkJC5dugQAmDx5Mrp06YKSkhLx9cnJyXB2dsbixYub3DYFsPwlxqUjJHAbAsasQUjgNiTGpbMuiRBipxwugBti6T7g+Ph4uLm5YciQIdi3bx92794NLy8vGojDTiXGpWOc51KMemihuIzzWkYhTAixCgrgWpYCGABiYmIwfPhwtG3bFh06dMD48eORkJDQrDYogOUtOHCLSfjWLSGB21iXRgixQxTANkQBLF9GoQqvjZiLUQ8tMAvgmT5rWZdHCLFDFMA2RAEsT4JgxKXsuZjrO9XiEfA7sz5kXSIhxA5RANsQBbA8Xc8PhpbXIPRQf/h4zjcJX58+8xB66Elkl+5jXSYhxM5QANsQBbD8JBd9Cy2vgZbX4ETaKFy7fBUhgdsw02ct3pu9Dj9pB0HLaxDBeyK3LIp1uYQQO0IBbEMUwPKSqftNDN+jKf9AuT7NbJ388mhE8J7Q8hpEJj2CgopTDColhNgjCmAbogCWj/rBGpX0OHRVVxtcN6fsICL43rXr9kNx5QUbVkoIsVcUwDZEASwPJZVXEJXUTzy1XFB+8ravydT9Ci3fC1peg8PJA24Z2IQQ0hgUwDZEAcxeeXUKjiY/U3vquReySvc2+rVpJaHiKeufI57D+3M20IhZhJBmowC2IQpgdhLj0rFyzma8OuJNzPWditBD/ZFS/H2Tt5Nc9E1Nb+k+82jELEJIi1AA2xAFMBuJcekY51U3xGTNQBtjPRc2OzDf+u87NGIWIaTFKIBtiAKYjZDAbZIGZsCYNRa3RyNmEUKaggLYhiiA2agJTOmGmJQ60AkhjokC2IYogNlYOec7SQOz5pT2MpNT2nQNmBDSVBTANkQBzMbJmO8l7zSVGJeO5f99F9NGzMBc36m4dOG4hBUTQhwBBbANUQCzcS7zDYQe6o95ftMxU8LbhsqrU8TbkhILPpGgUkKII6EAtiEKYNurMuSLo1hdzVsh+fbPpE+EltfgeKo3BEGQfPuEEPtFAWxDFMC2l1ayQzxKLag4Lfn2k4s2i9svrrwk+fYJIfaLAtiGKIBt72zmq7WTLQyCIBgl336FPlMM4Bv5NG8wIaTxKIBtiALYtqoMudDyD0PLa3At732rtROT8RK0vAbHUoZaJeQJIfZJEQFcXl6Omzdv4sqVK8jOzrZWM1ZHAWxbqcXbxKPToopYm7RTWBFjtXYIIfZFtgGclpaGoKAg9O/fH05OTlCr1eJy1113YeLEiQgPD4fRqJwjDgpg24rNeNkmR6ZVhrx6R9rvWa0dQoh9kV0AZ2RkYPr06WjdujU6duyI0aNHY/ny5fj888+xceNGfPjhh5gxYwb69u0LlUqF7t27IzQ0VIqmrY4C2HYq9Vni9IHXbXBt9mzm69DyGhxJfhpGQW/19gghyie7AHZzc8M///lP7Nq1C9XV1bdcNzExEcuXL8cdd9yBjz/+WIrmrYoC2HZSirfU6518wertpZeEie3ll9OgHISQ25NdAGu12ia/pqioCOfOnZOieauiALadmIx/Q8trEJ36T5vcn1ttKEIE7wktr0Fc7ltWb48QonyyC2B7RgFsGxX6jHq3Bq22WbsXsvyh5TU4nDwARqHKZu0SQpSJAtiGKIBtI7noWzGASypt915n6n4T280ti7RZu4QQZZJdAL/77ruNXt57T1k9TimAbeN0+otMhofUG0sRmdQXWl6DyzkLbNYuIUSZZBfAKpWq0YtarZaiSZuhALa++hMkJBQ0b77flriUPQdaXoOopH4wGCts3j4hRDlkF8D2jALY+pKKNogBrKu6avP2s0v3i+1nl/5p8/YJIcpBAWxDFMDWdyp9PLS8BifSRjGZnchgrEBU0mPQ8hpczJ5t8/YJIcqhiAC+ceMGNmzYgODgYGzcuBE3btywVlNWRQFsXWXVfL35eT9nVseVnIXQ8hpEJj0CvVHHrA5CiLzJOoAFQcCsWbPQqlUrk2u/rVq1wuzZyju6oAC2Lr5wvRjApVUJzOrILYsS68jU/cqsDkKIvMk6gNesWQOVSoX//ve/OH36NNLS0nD69GkxlNeutX0nm5agALauk2k+0PIanEzzYVqHUajC4eQnoeU1OJ/lz7QWQoh8yTqAe/fujTlz5lh8bs6cOejdu7fUTVoVBbD1lFYliEedNwu/ZF0O4nOXQ8trEMF7otpQyLocQogMyTqA27RpgwMHDlh87sCBA2jTpo3UTVoVBbD1JBZ8LgZwWTXPuhzklx8X60kv2cm6HEKIDMk6gLt06YJNmzZZfG7Tpk246667pG7SqiiArUMQBJxIex5aXoNT6eNZlwMAEAQDjiY/Ay2vwdnM6azLIYTIkKwD+OWXX0bXrl0RE2M6yfm5c+fg4eGBV155ReomrYoC2Dp0VVfFo02+cAPrckQRJ1Zgru9UvOLtj5VzvkViXDrrkgghMiLrAE5JSUG3bt2gVqvh5eUFb29veHl5Qa1W495770VqaqrUTVoVBbD0EuPS8dZ/38Yr3jMw13cqrlyMZV0SgJq6xnkuwaiHFmLUQwsw6qGFGOe1jEKYECKSdQADQG5uLpYsWYKnnnoKDz30EJ566iksW7YMubm51mjOqiiApZUYl45xXstkGXIhgdtq6zJdQgK3sS6NECITsg9ge0IBLC05h1zAmDUWa5vpo6xb5wgh1kMBbEMUwNKSc8jJ+ccBIUQeZB/Av/zyCyZOnIgBAwbgkUceMVn69u1rjSathgJYWnIOOUunx8d6LpLF6XFCiDzIOoBXrVoFlUqFLl264JlnnsGzzz5rtigJBbC0Eq6kwcdzniyvAQM1IRwSuA2vjXgTc32nYt9RP9YlEUJkRNYB3L17d/znP/+BwWCQetNMUABLq6z6JkIP9cdc36nwff5thARuk0341heXu0ycI9goVLMuhxAiE7IO4Pbt20Or1Uq9WWYogKWVVvKTeP9vSVU863IalKkLF+ssqjzPuhxCiEzIOoBHjRqFL774QurNMkMBLK3LOfOh5TU4nDwAgmBkXU6DKvQZYgAnFclnoBBCCFuyDuCrV6+id+/e+PXXX1FVVSX15m2OAlg6giDgWMoQaHkNLmQFsC7ntqJT/1k7O5Iv61IIITIh6wA2GAyYPXs21Go1WrVqhfbt25ssHTp0kLpJq6IAlk55dYp4VJlctJl1ObcVl7uk9jrwYzAKetblEEJkQNYBPG/ePKhUKjz22GN46aWXMH36dLNFSSiApZNe8rMYwMWVl1iXc1sZut316r3AuhxCiAzIOoDvuOMOLFmyROrNMkMBLJ0rOYtrjygfhyDIv5d8eXVqvSP2b1iXQwiRAVkHsLu7O/WCJhZFpw6rvaaqnHtro1PqavZnXQohRAZkHcCTJk3CypUrpd4sMxTA0qjQp9frVbyRdTmNdiVnkaKO2gkh1iXrAL548SJ69+6NTz/9FDdu3EB+fr7Z0lTnzp3D6NGjcd9998HFxQWdOnXCwIEDsXXrVrN1Y2NjMXz4cLRr1w7u7u6YMGECEhMTm70/FMDSyND9osj7atNLwv66b7mSPgOEODpZB7BKpYJKpYJarW5waarIyEj4+/tj69atiIiIwJ49ezBlyhRwHIf3339fXC8+Ph7t27fH4MGDsXfvXuzatQuenp7w8PBATk5Os/aHAlgaSh1ZyrTn9nesyyGEMCbrAA4KCsKKFStuuUjlqaeewn333Sf+/6RJk9C5c2cUFxeLjyUlJaF169ZYtGhRs9qgAJbG8VRvaHkNzmW+wbqUJlHavcuEEOuSdQDb0pgxY9CjRw8AgF6vh6urK/z9zTvLjBgxAj179mxWGxTALVepzxKPIvnCr1iX02RKGb2LEGJ9sgvgfv364YMPPkBcXJwUm2uQ0WiEXq9HTk4OvvzySzg5OeHrr78GUDMCF8dx+PLLL81et2DBAqhUKlRUVDS5TQrglsvU7REDuLAihnU5TZZeslMR41cTQqxPdgH8v//9D927d4darUbv3r2xfPlynDt3TopNm/D39wfHceA4Ds7Ozli/fr34XHR0NDiOQ2hoqNnrgoODwXEcMjIybrn97OxsXL582WQJDw+nAG6h+Ny3oeU1iEzqC6OgvOFJy6p5MYBTirewLocQwpDsArjO6dOnsXjxYvTs2RMqlQoPPPAAFi5ciFOnTkmy/eTkZJw5cwZ79+5FQEAA1Go1Vq9eDeCvAN6xY4fZ6+oCODMz85bbDwoKEgP+7wsFcPOdSBsFLa/B2cxXWZfSLIIg4GjKIGh5DS5mz2JdDiGEIdkGcH3nz5/H22+/jT59+kClUuHee+/FnDlzcOTIEQiCIEkbAQEBcHJyQk5OjiSnoOkIWHpVhlzx6PFm4TrW5TTb5Zx50PIaHEl+kq4DE+LAFBHA9cXFxeH9999Hv379oFKpcPfdd0uy3c2bN4PjOJw8eVLshBUQYN5TdeTIkdQJi5Hs0n1iABdUSHMmhIW0klBxP3RVV1mXQwhhRHEBXF9CQgJWrVolybamTZsGtVot3uM7efJkdOnSBSUlJeI6ycnJcHZ2xuLFi5vVBgVwy1zNe7fm+i/vBYOxknU5zVZanSgGcGqx+QAwhBDHoOgAbg4/Pz/Mnz8fP/30E6KiovDzzz/j3//+NziOw8KFC8X14uPj4ebmhiFDhmDfvn3YvXs3vLy8aCAOhk6m+UDLaxCb8TLrUlpEEAQcTX6m9jrwbNblEEIYkV0Ajx07tlHLuHHjMGnSJCxZsqRJvaQ3b96MwYMHo3PnznByckLHjh0xdOhQi0NRxsTEYPjw4Wjbti06dOiA8ePHIyEhodn7RgHcfFWGfPGoMbHgU9bltNil7MDa68ADJevHQAhRFtkF8P3334/u3bs3aunSpQtatWoFJycnhIWFSdG8VVEAN19O6QExgPPLj7Mup8VSi7eL+1NadYN1OYQQBmQXwE2Vl5cHHx8fPPLIIyyabxIK4Oa7lvcBtLwGEbwnDMZy1uW0WGnVjXrXgbezLocQwoCsA9hovPUtGnWzIf35559o3bq11M1LjgK4+U6lvwAtr0FMxr9ZlyIJQRBwJHkgtLwGl7IDWZdDCGFA1gHs59fwZOvFxcXo378/ACA3Nxdbtsh/VCEK4OapNhRBy/eCltcgoWAN63IkczF7NrS8BkeTn6HrwIQ4IFkHcNu2bREcHGz2eFlZGQYNGgQPDw+pm7QqCuDmyS3Tiqdr88qPsC5HMqnFW/+6Dlzd/HmmCSHKJOsA3r17N5ycnEzGZK6srMRzzz2Hzp07W33CBqlRADfP9fwPa6//9obeqGNdjmR0VVfFAE4rMR93nBBi32QdwACwdu1auLi44OjRozAYDBg7dizc3d0RGxtrjeasigK4eU6n/wtaXoMz6RNZlyIpQTDiSPKT0PIaXM6Zx7ocQoiNyT6AAeDNN9/EnXfeCR8fH7Rr1w7Hjh2zVlNWRQHcdHqjDlr+YWh5DW7kf8S6HMldzJ5Vcx045R90HZgQB6OIABYEAePGjYOLiwsOHjxorWasjgK46fLKosTTtLllkazLkVxK8RZx/8qqedblEEJsSHYB/Mgjj1hcevbsiQ4dOpg81rdvXymatBkK4Ka7kb+6NqAeht5YcvsXKExJVbwYwOklO1mXQwixIdkF8NChQ/Hss882elESCuCmO5MxGVpeg9PpE1iXYhWCYMTh5AG114EXsC6HEGJDsgtge0YB3DQGYxki+D7Q8hpczze/Hc1e/HH0v5jrOxWvjpiFkMBtSIxLZ10SIcQGKIBtiAK48RLj0vHum5/iFe8ZmOs7FTFnd7MuySoS49Ix1nMRRj20EKMeWoBRDy3EOK9lFMKEOADZBfDJkyeb/Jry8nJFhBoFcOMkxqVjnNeyv4XSUrsMpZDAbbX7abqEBG5jXRohxMpkF8Bt2rTBuHHjGtXbOSsrC6tWrcLdd9+NDz/8UIrmrYoCuHEcKZQCxqyxuK8zfdayLo0QYmWyC+CEhASMHz8eKpUKHh4emDJlClatWoUtW7Zgx44dWL9+PebPn49BgwahdevW6NKlC7766qvbTtwgBxTAjeNIoeRIPzYIIaZkF8B14uLiMHv2bDz44INQqVQmS9u2bTF8+HBs3rwZlZWVUjZrVRTAjeNIoeRIp9sJIaZkG8D15eTk4Pz58zhx4gRu3LiB6upqazVlVRTAjeNoHZMS49Kx4s1PMW1EXYez31iXRAixAUUEsL2gAG68P4/Nqb01578OcWuO3lhSb8pF+zvVTggxRwFsQxTAjSMIAo6mDKqdrH4O63Js5mSaD7S8BmczX2VdCiHEBiiAbYgCuHHK9Wni8IzJRd+xLsdm4nPfhpbXICqpH4yCnnU5hBArowC2IQrgxsnU7REDuKjyPOtybCajZJe43yWVyprrmhDSdBTANkQB3DjX8t6HltcgkveCUahiXY7NlFXfFAM4tdj+enwTQkxRANsQBXDjnE6fAC2vQUzGFNal2JQgCDiS/CRNzECIg5B1AH/88ccoKCiQerPMUADfXs0EDL2h5TW4kS//0c2kdj7LH1peg+jU4axLIYRYmawDuFWrVmjbti18fX1x/rzyrwVSAN9eQcUp8TRsdul+1uXYHF/4lbj/VYZc1uUQQqxI1gGclJSEhQsXonPnzlCr1Rg0aBBCQ0Oh1yuzhygF8O3xhRvEAKrUZ7Mux+YKyk+K+59Tevvx0AkhyiXrAK5TWVmJb7/9Fk888QRUKhXuueceBAUFISMjw1pNWgUF8O1dEE/BDmNdChOmp+BXsS6HEGJFigjg+o4fP46hQ4dCrVbD2dkZU6ZMUczpaQrgW6vphPRUbSekeazLYeavTmhTWZdCCLEiRQXwwYMH8cILL6BVq1bo3LkzXn/9dXh4eMDJyQmbNm2ydvMtRgF8a2XVfL3bcLayLoeZq3nv1tyGlfSIQ92GRYijkX0Al5SU4LPPPkOvXr2gUqnwyCOPYNOmTaioqAAAVFdX44033kC3bt2s0bykKIBvLUP3S72BKBz3PcrU/Sa+D8WVF1iXQwixElkHcEBAANq3b49WrVrhhRdeQEREhMX1jh49CpVKJXXzkqMAvrW6oRgjkx516KEYy6tTxABOKd7CuhxCiJXIOoA7duyIefPm4ebNm7dcLzc3F1u2yP+LigL41k6ljYWW1yA28xXWpTBlOhlFIOtyCCFWIusALi0tlXqTTFEAN0xv1NWbjm8N63KYu5g9C1peg2MpQ1mXQgixElkHsL2hAG5Yfnm0eNo1t8zypQZHklz0jfh+VOgzWZdDCLECWQfwsGHDGlyGDx+OF198ER9++CFycnKkbtoqKIAbdrNwXb0RoPJZl8NcYUVMvRHB9rEuhxBiBbIO4KFDh8LDwwMqlQo9evTAwIED0b17d6hUKnTr1g39+vVDmzZt0LVrV1y/fl3q5iVHAdywc1n/gZbX4ETqCNalyILBWIkI3hNaXoPr+StZl0MIsQJZB/DBgwdx//334+TJkyaPnzhxAv/3f/+Hffv24ebNm3jwwQcxefJkqZuXHAWwZYJgxOGkJ6DlNbiSs5h1ObJxJn0StLwGZ9InsS6FEGIFsg7gAQMG4Ntvv7X43DfffIN+/foBAL7++mvcddddUjcvOQpgy0qrboinW9NKdrAuRzau5wdDy2sQwXvCYKxkXQ4hRGKyDmAXFxccOHDA4nP79++Hi4sLACAqKgrOzs5SNy85CmDL0kt2igGsq7rKuhzZyC7dJ74vhRUxrMshhEhM1gF8//33Y+bMmRafCwgIwP/93/8BAH777Td06dJF6uYlRwFsWVzuUmh5DaKSHocgGFiXIxsV+kwxgJOK5D/UKiGkaWQdwCtXroRKpcJLL72EX375BcePH8cvv/yCKVOmQK1WIzg4GACwZMkSPPfcc1I3LzkKYMtOpI2CltfgbObrrEuRnWMpQ6HlNbiYPYt1KYQQick6gAFg6dKlcHFxgVqthlqthkqlgouLC5YtWyauc+HCBSQkJFijeUlRAJurNhSKR3mJBZ+zLkd2LmXPhZbX4GjyMxAEgXU5hBAJyT6AAaCwsBB//PEHtm3bhj/++AMFBQXWasqqKIDN5ZZFiQGcV36EdTmyk1L8vfj+lFensC6HECIh2QZweXk5nnnmGRw8eFDKzTJFAWwuseCT2oDphWpDMetyZKe48oIYwJm6X1mXQwiRkGwDGADat2+PyMhIqTfLDAWwubOZr0LLa3AybQzrUmTJKFQjMqkvtLwGV/PeZV0OIURCsg7gESNG4KOPPpJ6s8xQAJsSBAOikvpBy2sQl/sW63JkKybjJWh5DU6lj2ddCiFEQrIO4EuXLuGBBx7A999/D51OJ/XmbY4C2FRJVbx4ejW95GfW5cjWjfzVtQNy9IbeaF8zhBHiyGQdwG5ubnB2dhZ7QLu5uaF9+/bi0qFDB6mbtCoKYFOpxdvFAC6tkn8vdlZyyg6K71NB+QnW5RBCJCLrAH7ttdcwffr0Wy5KQgFs6krOQmh5DQ4nPwlBMLIuR7aqDHliAPOFX7EuhxAiEVkHsL2hADZ1PPU5aHkNzmf5sS5F9v56r2awLoUQIhEKYBuiAP6L6VHdetblyJ7p2QIakIMQeyD7AI6Pj8eUKVNw9913o3Xr1oiNjQUArFixAhEREU3enlarxeuvv45evXqhbdu28PDwwLhx4xATYz7YfWxsLIYPH4527drB3d0dEyZMQGJiYrP3hQL4L/Wva+aXH2ddjuyZXC+vbv5nkBAiH7IO4HPnzsHNzQ333HOPOP5zXQAvWLCgWXMAT5w4EcOGDcP69esRFRWFsLAwDBw4EE5OTtBqteJ68fHxaN++PQYPHoy9e/di165d8PT0hIeHB3Jycpq1PxTAf6nr2avlH6aevY1Qv8d4Rsku1uUQQiQg6wAeOXIknnnmGZSVlUGv10OlUokBvHPnTvTo0aPJ28zOzjZ7TKfToWvXrhg+fLj42KRJk9C5c2cUF/81OlNSUhJat26NRYsWNWNvKIDr++ve1hdYl6II9e+Zjs9dzrocQogEZB3Abm5uCA8PBwAYDAaTAD58+DBcXV0la2vYsGHQaDQAAL1eD1dXV/j7+5utN2LECPTs2bNZbVAA16gZ3emR2tGdVrAuRzH2Hv0P5vpOxasjZiMkcBsS49JZl0QIaQFZB3C7du2wb98+AOYBHB4ejo4dO0rSTlFRkXiNFwCuXr0KjuPw5Zdfmq27YMECqFQqVFRUNLkdCmAgMS4d783+Aq94z8Bc36k4FbOddUmKkBiXjrGeCzHqoYUY9dACjHpoIcZ5LaMQJkTBZB3ATz/9NPz8am5R+XsABwQEmJwybomXX34ZTk5OYkes6OhocByH0NBQs3WDg4PBcRwyMjJuuc3s7GxcvnzZZAkPD3foAE6MS8c4r2WmIeK5lEKkEUICt9W+b6ZLSOA21qURQppJ1gG8Y8cOqNVqzJ07F7GxsVCpVPjll1+wevVqODs747fffmtxG8uXLwfHcVi3bp34WF0A79ixw2z9ugDOzMy85XaDgoLAcZzFxVEDmEKk+QLGrLH43s30Wcu6NEJIM8k6gAFg5cqVaN26NdRqNVQqFdRqNZydnRESEtLiba9YsQIcx2HlypUmj0txCpqOgM1RiDQf/XghxP7IPoABIDU1FZs2bcLKlSuxYcMGJCUltXibdeG7YoV5J6C6TlgBAQFmz40cOZI6YTUThUjzWTx970Wn7wlRMkUEsNTee+89cByH5csbvp1j8uTJ6NKlC0pKSsTHkpOT4ezsjMWLFzerXUcP4JqORIupI1EzJcalI2jWR5g2oqYD24XzUaxLIoS0gCICODs7G6dPn8bhw4fNlqb6+OOPwXEcRo0ahRMnTpgtdeLj4+Hm5oYhQ4Zg37592L17N7y8vGggjhbSnngbc32nYtqIAAQH/kDh20Tl1anigBzJRd+xLocQ0gKyDuCMjAwMHz5cnI5QpVKJ14Hr/ttUQ2ZyAQgAACAASURBVIcObbBzFMeZ7kJMTAyGDx+Otm3bokOHDhg/fjwSEpo/bR4FMHAyzQdaXoPYzFdYl6JIgiDgWMpgaHkNLma/ybocQkgLyDqA//Wvf+GOO+7A6tWrsX//fkRFRZktSuLoAVxtKBKP3hILPmFdjmJdyp4LLa/B0eRnaGIGQhRM1gF85513YvPmzVJvlhlHD+DcskgxgPPKj7AuR7FSi7eK72NZdcs7JBJC2JB1AHfu3BkHDx6UerPMOHoAJxR8XG8CBh3rchSrpDKOJmYgxA7IOoBnzpyJwMBAqTfLjKMHcEzGVJqAQQI1EzM8Bi2vQVzuMtblEEKaSdYBHBkZiYceegizZ8/GgQMHEBsba7YoiSMHsMFYiQjeE1peg2t577EuR/HOZb4BLa/BidSRrEshhDSTrAO4rtdzXY/n+ktze0Gz5MgBXFgRI542zSrdy7ocxbtZ+IX4flYZ8lmXQwhpBlkH8JYtW267KIkjB3BS0QYxMCr0tx5Hm9xeQfkJ8f3MKbWffhKEOBJZB7C9ceQAPp/lDy2vQXTqMNal2AWDsQwRfB9oeQ2u53/IuhxCSDMoJoCvXr2KY8eOobS01JrNWJWjBrAgGHE4eQC0vAaXcxawLsdunEmfCC2vwZn0SaxLIYQ0g+wD+Pvvv0e3bt3Ea791Ha8mTZqEjRs3WqNJq3HUANZVXRdPl6aVmM+xTJrnen4ItLwGEbwnDMZy1uUQQppI1gG8c+dOqFQqjB07FuvXr4dKpRIDOCQkBM8995zUTVqVowZwWkmoGMC6quusy7EbOaUHxPe1oOIU63IIIU0k6wB+7LHH8MYbbwAADAaDSQCHh4fDw8ND6iatylED+HLOAmh5DQ4nD4AgGFmXYzeqDHliAPOF61mXQwhpIlkHsIuLCw4cOADAPICPHDmCNm3aSN2kVTlqAEenDoOW1+B8lj/rUuzO8dQR0PIanMv6D+tSCCFNJOsAvuOOOxAWFgbAPIBDQ0Nx9913S92kVTliAFfoM8WjtKSiDazLsTtxuUug5TWISnocgmBgXQ4hpAlkHcBjx47F0KFDIQiCWQCPHDkSU6ZMkbpJq3LEAM4q3SsGcGFFDOty7E56SZj4/pZUxbMuhxDSBLIO4DNnzqBNmzbo378/1q5dC7Vajbfeegs+Pj5wdXXFpUuXpG7SqhwxgK/lvQctr0Ek7wWjUMW6HLtTVn1TDODU4m2syyGENIGsAxgAIiIi8PDDD5sMS6nRaBAZGWmN5qzKEQP4VPoL0PIaxGS8xLoUuyQIAo4kD4SW1+BSzv9Yl0MIaQLZB3CdhIQEREdH49q1a9ZsxqocLYD1xhJo+Yeh5TVIKFjDuhy7dTF7FrS8BsdShrAuhRDSBIoJYHvgaAGcV35EPD2aWxbFuhy7lVy0ud442+msyyGENBIFsA05WgAnFKytDYZeqDYUsS7HbhVXXhADOFP3G+tyCCGNRAFsQ44WwLEZr0DLa3AyzYd1KXbNKFQjMqkvtLwGV/OCWJdDCGkkCmAbcqQANgpViEx6hELBRmIz6ccOIUpDAWxDjhTARRXn6LSoDSUWfEKn+wlRGApgG3KkAE4u+oY6BtkQdXgjRHlkHcB79uzBunXrLD73xRdfYO/evVI3aVWOFMAXsmbW3hozlHUpDkFv1NEtX4QojKwDeODAgQgJCbH43EcffYRBgwZJ3aRVOUoA1wwO8SS0vAaXc+axLsdh0KAnhCiLrAO4Y8eO2L9/v8XnDhw4gE6dOkndpFU5SgCXViXQ8IgM0LCfhCiLrAPY1dUVe/bssfjcnj174OLiInWTVuUoAZxeslMMYB1NEGAz9Se+KKo4y7ocQshtyDqA+/XrB39/y3PI+vv7o2/fvlI3aVWOEsBXchZDy2twOOkJmiLPhkynftzEuhxCyG3IOoA///xztGrVCu+88w6ysrIAAFlZWQgKCkKrVq3w2WefSd2kVTlKAB9PfQ5aXoPzWb6sS3E40SnDoOU1uJAVwLoUQshtyDqABUHAtGnToFKpoFar0bp1a6jVaqhUKrz66qtSN2d1jhDAlfps8SiML/yKdTkO53LOfGh5DY4kPwlBMLIuhxByC7IO4DpHjhzBsmXL4Ofnh2XLluHo0aPWasqqHCGAs0v/EAO4oOI063IcTlrxj+L7X1qVwLocIqHEuHSEBG5DwJg1CAnchsQ4y/fXS70esR5FBLC9cIQAvpb3AbS8BhG8JwzGCtblOBxd1TUxgNNKfmJdDpFIYlw6xnktw6iHForLOM8lOHvuD+SUHkSG7hekFm/D0VOfYaznQpP1xnouwokz3yBT9ytyyg4iv/w4LpyPwjjPJabb81pGIWxjFMA25AgBfDp9ArS8Bmcy/s26FIckCEbsjBiKub5T8frIhXRkowANHYkahWqUVMUjvWQnls5cbBKWdctc36niDy4tr8Fc36ktWi9o1iqUVieKly/oKNm6ZBfAPXr0wPnz5wEA3bt3R48ePRpcHnjgASmatBl7D+DrV25gru9LeMV7BpbNfIv+sTKQGJeOsZ4Lar9QF9CRjcxZOrId67kQvx15RZzMRMtr8Ir3DPHv+deyANNGzDAJVqnWi0rqhz1HXjU7mqbPkrRkF8DTp0/HzZs3AQCvvfYapk+ffstFSew5gBPj0jHOczF98TMWErjN4pFNSCANiCJHwYE/NOqI9X9+0yyu9+7sz1FSeQVl1UmoMuQ1uL2Vc75FWfVNFFdeQkH5Sax4c81t223oKJk+S9KRXQDbM3sOYPril4eAMZa/WGf6rGVdGqnHKFQhrSQUr46YZfFI9LWRgbiW9x4ydLtRWnUDCVdSza8BW/iBa/FacaPXW4qL548greQnxOe+jddGzGmgtjnILt1PvewlQAFsQ/YcwDVf/H//x0pf/LZGP4TkzShUIa34RxxLGdrko8y667EzfdY2qndzS9dr6LNUd5R8Ms0HWbrfabCdFpB9AAuCgFOnTmHnzp34/vvvzRYlsecADg7cQl/8MmB6ZEOXAliq34EpOPAHHD+zCcdShpicWt4V+TzGei6S5XVWy9enF2FnxLMm+3AidQROxmxFcOBW6qzVRLIO4GvXrqFXr17i4Bt/X9RqtdRNWpU9B3DM2d3w6TNPll8kjiYxLh1vz1qJaSNmYK7vVFy6EM26JIdjKbx8+sxD6KH+0PIaRKcOQ1rJTzAKVY0+YmXBUm16ow5JRRtwJHkgtLwGoYf607/9ZpJ1AD/33HO4//77ERYWhmvXriEpKclsURJ7DuCreSsQeqg/5vq+jIAxH8vui8TRlFXz4hFKctFm1uU4nIZO386f8QbSS3bCKFSzLrHFDMYyJBd9h3l+0+nsVzPJOoA7dOiAn36yn8EE7DmAj6d6Q8trcDZTeUOE2iNBEBCdWjMu9LnMN1iX43D8R39kuU/EmDWsS5NcwOiPLQTwAgSM+Zh1abIn6wD28PDAH3/8IfVmmbHXAC6vTqk3C88G1uWQWvG5b9fMD5z0CAzGStblOIxKfRYWzPB1mKPCho72F/nPhsFYxro8WZN1AL/11luYNm2a1Jtlxl4DOK0kVAzgksorrMshtbJL/xT/LvnldB3YFir06Tie+pxDXRe91fXumIyp0Bt1rEuULVkHcGhoKB544AGMHTsWGzduxK5du8wWJbHXAL6Y/WbtDDwD6d5AGak2FEHLPwwtr8GN/FWsy7F75dUp4ml/La+B9sQ7su1cJbX6nbWCA7/HvqN+4vtwOv1FVBsKWZcoS7IOYEs9n6kXtLwYBT0OJz0BLa/B5Zx5rMshf3MmfRK0vAan0l9gXYpdK6tOMrnF6Hr+SgiCwLosZoxCFS5kBYjvx6m0sagy5LEuS3ZkHcBRUVG3XZTEHgO4qCJW/EeWodvNuhzyN4kFn4p/nypDLuty7FJpVQKOpgwS3+cb+asdOnzrGIVqXMqe+9f9wmmjUKnPYl2WrMg6gO2NPQZwYsHn4j+wSn0263LI3xRWnBH/Ppm6X1mXY3d0VddwJPlp8T1OLPiUwrceQTDgSu4S8f3ZHTkWH8zZRAN21FJEAJeUlGD//v348ccfceDAAZSUlFirKauyxwA+k/FvcVg6Ij9GoRpRSY9By2twJWcR63LsQt31zhmjgzHP7zVxcI2bhV+yLk2WBMGIq3lBDtUxrbFkH8CrV6+Gm5ubyWhYbm5uWLNGeffT2VsAVxuKxU4+1/M/ZF0OaUDdtbijKYPo6KyFLA316dNnHo6dXse6NFkTBAHLZr7lMLdmNZasA/j777+HSqXC6NGjsXPnThw7dgw7d+7EmDFjoFar8cMPP0jdpFXZWwBnl+4XTy3llR9lXQ5pQGrxNvHvpKu6xrocRaPJLpqPZuoyJ+sA7tevH15++WWLz7388st47LHHpG7SquwtgMWBHngvGIwVrMshDaBhKaVDIdJ89OPFnKwD2MXFpcGRsP744w+4uLg0a7slJSVYuHAhvL290blzZ3Ach6CgIIvrxsbGYvjw4WjXrh3c3d0xYcIEJCYmNqtdewrg+kMdns18nXU55BYEQUB0Cg1LKYX3Zq+jEGmmhk7fx1+6yLo0ZmQdwJ06dcL27dstPrd9+3Z06tSpWdvleR7u7u4YMmQIfH19Gwzg+Ph4tG/fHoMHD8bevXuxa9cueHp6wsPDAzk5OU1u154CuKw6qd5R1TesyyG3EZ+7vHZYyr40LGUzGYUq/BI1gToStUBdBza/0e9hru9UhB7qj/NZ/g7bN0HWATxu3Dj07NkT6emmH+7MzEz06tULL7zQvMEFBEEQ/+C5ubkNBvCkSZPQuXNnFBcXi48lJSWhdevWWLSo6T1K7SmATa8rXmVdDrmN7NI/6g1LeZx1OYpUd0916KH+WP7fdx1ihCtrupyzoN4tcr+xLocJWQfwpUuX0L59e7i6usLHxwd+fn7w8fFB27Zt4e7ujitXWj7ucEMBrNfr4erqCn9/f7PXjBgxAj179mxyW/YUwNSzVlmqDYU0LGUL6KriEcH3EUcVs4fpBFmrNhSI91AfSX7SIUfKknUAA8C1a9cwZcoUdO3aFU5OTujatSteeuklXLsmTW/OhgL46tWr4DgOX35pfm/fggULoFKpUFHRtI5H9hLANfeW9qN7SxXmTPpEGpayGYyCHqfT/wUtr0EE3xsllcr+9ysn9c/MXMoOZF2Ozck2gA0GA65du4bCQusO4t1QAEdHR4PjOISGhpq9Jjg4GBzHISMjo8HtZmdn4/LlyyZLeHi4XQQwja6kTIkFn9QbltLxjjaaK7noG/F9SyigOW6ldjF7tvj+5pQeYF2OTck2gPV6PVq1aoV9+/ZJuVkztwvgHTt2mL2mLoAzMzMb3G5QUBA4jrO4KD2ATb/IaXxhpTD94eSY19yaqqw6CZFJj9SMZZw6gm63s4JKfQ4OJw8QL2k50sxJsg1gAOjWrRt+//13qTdrwlqnoO35CJhOZSqT6aWDxazLkT1BMCI24xXxR0thxRnWJdmtDN0v4vvsSJ9NWQfw/PnzMXnyZKk3a+J2nbACAgLMXjNy5EiH7YRV05mnF3XmUagLWf7Uea6R0kp2iKFwNW8F63LsmiAIOJ/l+9fIemXKmumuuWQdwFu3bsW9996LYcOGYd26dfj555+xa9cuk6WlbnUb0uTJk9GlSxeTyR+Sk5Ph7OyMxYub/ivNHgI4u3Qf3c6iYKnFW+vdPnaddTmyVaHPRFTS49DyGhxLGQq9Uce6JLtXoc8QJw5xlPdc1gFcN/lCQ4tarW72tvft24ewsDBs3rwZHMdh0qRJCAsLQ1hYGMrKygDUDMTh5uaGIUOGYN++fdi9eze8vLwceiCOuNy3xAEdjEIV63JIE5VV36w3gMp3rMuRpZqjMX/xfcp1kKMxOag76xB6qD+Wzlxi99MWyjqAIyMjERUVdculue6///4GO0nxPC+uFxMTg+HDh6Nt27bo0KEDxo8fj4SEhGa1qfQAFgQBx1KGQstrcD7Ll3U5pBlq/obP1gxLmfUf1uXITmJcOoJmrcIr3jMw13cqDkQvYF2SQxEEAb8ffd1hRhuTdQDbG6UHcGlVgnhUkEJHT4pV/ywGDUv5l5qxipeajFU8zmupXX7xy9n7czY6zHjbsg7gHj164Pz58xafu3TpEnr06CF1k1al9ABOKf5eDODSqhusyyHNRNfxLaPZeuTBkWacknUAq1QqnDp1yuJzZ86cadE1YBaUHsDns2bUdpAYTD1oFcy0J/tq1uXIxozRH4pHvvb+xS9njvRDSPYBfPr0aYvP/fTTT3B3d5e6SatScgAbhSpEJj0KLa9BXO4S1uWQFjqd/mLtvdzjWZciG0sCzL/07fWLX84sTVs41nORXV4KkF0Ab9myBcOGDcOwYcOgUqnwxBNPiP9ftwwcOBCurq4YNWqUFE3ajJIDuKD8pHjaMqt0L+tySAslFKylYSnrKa1KQOihJx2m84/c1U1b+NrI2ZjrOxU7I4ai2lB8+xcqjOwCeP369fDy8oKXlxfUajUefPBB8f/rlieeeALTp09HcnKyFE3ajJIDOKFgTe0Xdi9UGfJZl0NaqKDiNEIP9cdc36nwfX6FXd/q0RiXsufW3v4yAO/P/pqmGpSJnLKD4g/FxILPWJcjOdkFcH3du3dvsBOWEik1gBPj0rFwRgBe8Z6BhTMC6EvJDty4kkRHe7V0VVfFL/nLOfNYl0PqEQRBnIkqKukxVBsKWJckKVkHsL1RYgBbvjXDMb+o7YkjdXS5nYvZs2oD+GGUVd9kXQ75m7zyI+IPJHvrNCjrAN6zZw/WrVtn8bkvvvgCe/cq61qkEgOYvqjtkyPd6nErxZWXxC936lwoT4IgICZjinjveqW+6aMQypWsA3jgwIEICQmx+NxHH32EQYMGSd2kVSkxgGu+qOnWDHtDP6xqnM/yg5bXIILvg/LqFNblkAbU7wR6Le8D1uVIRtYB3LFjR+zfv9/icwcOHECnTp2kbtKqlBjAK+dspi9qO2TpVg9HG/WpqOKs+KUen/s263LIbZzNfK3mKJj3QoW+4bnYlUTWAezq6oo9e/ZYfG7Pnj1wcXGRukmrUmIAHz+ziTrr2Km6cY+njagZ9zjm7G7WJdmUPX6h27P6P5iu5r3DuhxJyDqA+/XrB39/f4vP+fv7o2/fvlI3aVVKDOAz6ZMQeqg/5vu9gZl2PjOJIzIYyxCV1A9aXoNL2XNYl2Mzpqc032ddDmkke7tkIOsA/vzzz9GqVSu88847yMrKAgBkZWUhKCgIrVq1wmefKeu+MKUFcGl1ovgllVBA13zt1ZWcheKRoD0OdvB3NZ16ptplpx57V7/T3BU76DQn6wAWBAHTpk0T5/5t3bo11Go1VCoVXn31VambszqlBXBiwSd/Tb5Qnci6HGIl+eXHxL9zeslO1uVYXV750Xq3tXzEuhzSRBez3xRvG1P695KsA7jOkSNHsHTpUvj5+WHZsmU4evSotZqyKiUFsCAYEZ06DFpeg9PpL7Iuh1iRIBhwNGUQtLwGsZmvsC7HqmoGdnixdmCHfjSqmwLpqq6Jk4kofeAURQSwvVBSABdWnPlr7t/iH1iXQ6zsen6I+Peu0NvvNf7cMm29oQ0/YV0OaabLOfPEoXF1VVdZl9NsigjgP//8E0uWLIGvr684/vPp06eRk6OsazdKCuD43OViZwc6SrB/JZVXxGBKKtrAuhyrEAQjTqW/AC2vweGkJ1BtKGJdEmmmsuqb0PIPI/RQfywO+B8CFNpBVNYBXFZWhueee068BqxWqxEbGwsAmDx5MubPny91k1allAA2GCtxOOkJaHkNLmRZ7oVO7IsgCDiZNhpaXoOTaWPsbr7nxLh0BL25Gq9419xydeQUHf0q3cHjixV/i6SsA3ju3Lno2LEjdu/eDZ1OB5VKJQbwN998g0cffVTqJq1KKQGcXbpPPBrKLt3HuhxiI3zh1+LfvaQyjnU5kqFBR+zTB3M2KX6QIFkHcLdu3cSxoA0Gg0kA00hY1nMhy188TWcwVrIuh9hIuT5NDODr+R+yLkcyNOymfbKH8cxlHcDOzs6IiIgAYB7Ahw4dopGwrKDKkI8Ivk/t4PRvsS6H2FhsxivQ8hocTRkEQTCwLkcSAWM+pvHM7ZA9/LCSdQD36NFDHGzj7wG8du1a9O7dW+omrUoJAZxavFU8CiqoOM26HGJj6SU7xb9/fnk063Ik8c6sYMV/URNzFi8teC5W1KUFWQfw/Pnz0aVLF5w9e1YM4LNnzyIpKQn33Xcf3nlHWeOBKiGAz6RPhJbXIDplGATByLocYmPVhmJE8J41Iw3lLGZdTosZjBXYFTlS8Z11iGWJcelYOWczpo0IwFzfqfjt8BRFdSCUdQCXlJSgX79+cHJyQr9+/aBWq/Hoo4/Czc0NTzzxBMrLy6Vu0qrkHsD1h56keyQd18Xs2eJAFQajsv6N/V1y0WZoeQ1CD/XHO7M+xEyftYq8XYXcWmLBp+J3V07pQdblNJqsAxgAysvLERISgn/84x/o2bMnnn76aaxcuRJlZWXWaM6q5B7ACQVrxQ9xWfVN1uUQRnJKD4qfg0yd5dnIlEBv1OFI8pPQ8hqcSB0Jo6BnXRKxEr2xBIeTB9TeRuejmP4Lsg9geyLnABYEI6JTaoaePJM+kXU5hCGjUCV+mZ3P8mNdTrPdLFxHt9M5EL5wQ70fjr+yLqdRKIBtSM4BXFBxWvzwphZvZV0OYexq3ju1I6H1RpUhj3U5TVZtKEBU0mPQ8hqcSh9P/RkcgMFYhqPJz0DLa3A89TkYhWrWJd2W7AL49ddfb/TyxhtvSNGkzcg5gONy36KhJ4mosCJG0WOB38hfJdafVxbFuhxiIynFP4h/97SSn1iXc1uyC+D7778f3bt3b9TSo0cPKZq0GbkGMA09Sf5OEARxNiylXZKo1GchMqkvtLwGMRlTFdUrlrSMUajCsZRnoeU1OJYyRPYDCckugO2ZXAPYdOjJP1iXQ2Si/nzQSuqUdzVvhVh3YcUZ1uUQG0svCfvr7E3Rd6zLuSUKYBuSawCfp6EniQWlVQn1bkv7jHU5jVJenSKO5HYu6z+syyEMGAU9TqSOgJbX4EjyQOiNpaxLapDsAvjw4cPQ6XS3XS83NxfffvutFE3ajBwDuP7Qk/G5y1mXQ2TmdPq/agZmSR2uiFO5V3IWiT8aiisvsS6HMJKl+138HPCFX7Eup0GyC2C1Wo1Tp06J/280GtG6dWucPXvWZL2TJ09CrVZL0aTNyC2AE+PS8fasD8Qp2s6fV84N7MQ2Uoq3IPRQf8z1nYoZz38g60EsdFXXoeV7QctrcDF7NutyCEOCYMSptLGyn/tZdgGsUqlMAvjvY0DXoQBuGctTtNHwfMTU1ctXFDOM48XsN2uPeh5GadUN1uUQxnLLtOJRcEKBPCfeoAC2ITkFsD3MJEKsTymfk+LKi+KXrT2MYU1aThAEnEmfJA6rWmXIZV2SGQpgG5JTANfMpUlTtJFbU8qcq+cy36i9j90T5dUprMshMpFfflwcC3zZzGUIGLNGVpdRKIBtSE4BHPTmx4o4siFsyf0IODEuHe/O/lzsx6A9oawZ0oj17T36hmwvo8gygH/44QfExsYiNjYWp0+fhkqlwvbt28XHYmNjsWXLFgrgZjIYy7Er8nnZfiiJfFjqK+DjOR8JV1JZl9bAfLBL6TNMTLw3+wvZ/oiUZQCr1WqT5VaPKYlcArhugIWaKdqCaYo2ckuJcekICdyG/zy/DHN9pyL0UH+klexgXZbsj86JPMj5cpvsAnjLli1NWpREDgFcVp0kTrh+On2CYqbtIuwZjOXiMH9Hkp9kfmvHjNEfyvaLlciHnH+oyS6A7ZkcAvh81gyxt2hRxTlmdRBlyi79Q/z8XMt7n1kdRqEKC2cEyPaLlciHpUsVYz0XIOFKGuvSKIBtiXUA55ZFiF+ecblLmdRAlE0QBMRmThOnKtRVXWNSx438jxB6qD/1YyCNUncZ5fVRC8TLKCnF37MuiwLYllgGsMFYiejU4TUjwyT3V+Qcr0QedFVXoeUfhpbX4GzmqzYfojKv/Kj4Q/K3w/9GcOAP1I+BNIreWILolJpZviJ5L+iqrjKthwLYhlgG8M3CL8QvrdTirTZvn9iXa3nv1ZtBa7/N2q0y5IqTrkcl9VPULE1EHmrmuq75AXkybQwMxgpmtVAA2xCrAC6vTkVk0iPQ8hqcShsLo6C3afvE/lQbCnEk+cmaiRpShtnkS0wQjDiX9R8x+DN0u63eJrFPiQWf1+vL8B6zOiiAbYhVAF/MnkXzoxLJpRX/KH6ubhaus3p7yUWbxfYu58xTxOxMRJ6Mgh4xGf8WP0+5ZRFM6qAAtiEWAZxXfqTel9YCm7VL7J8gGHAqbVzN9bSkvqjQW+/6a3HlJfH2uejUf0JvvP2UpYTcSnl1CqKSHq+9re4pVOpzbF4DBbAN2TqAjUIVjtdOTB2V9Bgq9dk2aZc4jsKKM+IPvEvZgVZpQ2/U4Xiqd23P6z4oqjxvlXaI48nU/Sp+fs9lvgFBMNq0fQpgG7J1APOFG8QPV3LRZpu0SRzPpZz/iXMG+z3/vuS9ka/kLBY/x0lFGyTbLiEAcDlngfj5Sin6zrZtUwDbjq0CODEuHR/M+QbTvAMw13cqwqMmwChUW7VN4rjiLl2odz+uNHNL19236ft8kHjf5tnM6TY/QiH2T2/UITr1n7U/Il/GjNHBNruljQLYhmwRwDWjviz92wD1S+j+SGI1Ug/1Z3ECiD7zEH/pksSVE1Ljwnmt5D8iG4MC+BZ0Oh0CAwNxzz33oE2bNnj00UcRGhra7O1ZO4AFQUDQm6tpeD5iU5bnDF4Av9HvNqunspzH7iX2idVnjgL4Fry9vdGxY0d8/fXXiIiIgK+vLziOw/bt25u1PWsGcFHlecRkTMEr3jNogHpiUw19ec31nYqzma9BV3W9scZwUQAAHPpJREFUUdsRBAF5ZYfx2ohA+gwTm2roR+T0kYEoKD9ptXYpgBuwd+9ecByHH3/80eRxb29veHh4wGBo+ixC1gjg8uoUXMqeK3YimOs7lY4eiE2ZnjJeKJ4yDj3UXxwz+nr+SlQbii2+XhCMyCk9iNPpE+gzTJi41Y/IuiFXiypiJW+XArgBvr6+cHNzg15vOmrUjz/+CI7jEB0d3eRttvTNruuYEjBmDVbO+Q6RJ1aK90bWfNF5IvLEynrXgGmAemIbdZ/NujGZr12+iqt574pD/tXca/k0TsZ8L36GQwK34nTsTziZ5iOuo+U1+En7D4z1XEifYWIzln5EjvVciJ+0g0w+m+ezfHHh/OF6n+GWddaiAG7AwIEDMWDAALPH696wDRuafjtES97s2x1lXMr5H8qrU8V1638Z0hcXYUVXFY/YjJeh5TUWZy+q/xk+nNwfNwvXodpQSJ9hYnOWPnNVhlxcz1+JSN6rwc9wS34cUgA3oGfPnhg5cqTZ4xkZGeA4DsHBwbd8fXZ2Ni5fvmyyhIeHN/vNbugUyUL//9LABETWBEFAlu53zPN73eJneJ7fa0gq2kCjWxHZqtBn4mpeEOb6vizp5REK4Ab07NkTo0aNMnu8LoBDQkJu+fqgoCBwHGdxac6bbbmTwELMHEMdU4gyBIy23EM/YMzHrEsjpFH8R39o+Xu4mR0EKYAb0NJT0LY6AqaOKUQp6DNMlE7qzzAFcAP8/PwsdsIKDQ1l0gnL0jVg6phClIQ+w0TppP4MUwA3YN++feA4Djt27DB5fNSoUcxuQ6KOKUTp6DNMlE7KzzAF8C14e3ujU6dO2LhxIyIiIuDn5weO47BtG/vTDYQQQpSNAvgWdDod5syZg7vvvhvOzs7o27evrIeiJIQQohwUwDZEAUwIIaQOBbANUQATQgipQwFsQxTAhBBC6lAA2xAFMCGEkDoUwDZEAUwIIaQOBbANUQATQgipQwFsQ7GxseA4DuHh4WbDVNJCCy200OJYS93wxLGxLZ9rmAL4NurebFpooYUWWmipW8LDw1ucLxTAt1FYWIjw8HDExsa2+BeTIxxF077a50L7ap8L7WvTl9jYWISHh6OwsLDF+UIBbAOXLzvOdWTaV/tE+2qfaF/ZogC2ATn+4a2F9tU+0b7aJ9pXtiiAbUCOf3hroX21T7Sv9on2lS0KYBuQ4x/eWmhf7RPtq32ifWWLAtgGsrOzERQUhOzsbNalWB3tq32ifbVPtK9sUQATQgghDFAAE0IIIQxQABNCCCEMUAATQgghDFAAW5FOp0NgYCDuuecetGnTBo8++ihCQ0NZlyWJo0eP4vnnn0fHjh3h4uKChx56CO+9957JOgcPHsTAgQPh6uqKO++8E6+99pqsOkDUV1JSgoULF8Lb2xudO3cGx3EICgoyWcdgMGDNmjUYOXIkunXrBldXVzz88MNYvHhxg6PifP755+jVqxecnZ3RvXt3rFixAtXV1TbYo4Y1Zl8BQBAEbNy4EY8//jjat2+PO+64A0OGDMHvv/9ucbty3FetVovXX38dvXr1Qtu2beHh4YFx48YhJiamwdcIgoDBgweD4zjMmjXL4jpy3Ndz585h9OjRuO++++Di4oJOnTph4MCB2Lp1q9m6sbGxGD58ONq1awd3d3dMmDABiYmJFrer9H2trq7GmjVr4OXlBRcXF7i7u+Ppp59GdHS02XorVqzA/fffD2dnZ/Tq1Quff/65VfeDAtiKvL290bFjR3z99deIiIiAr68vOI7D9u3bWZfWItu3b4darcaUKVPw22+/ISIiAps2bcK7774rrhMVFQUnJye88MILOHDgALZt24Zu3brBy8sLlZWVDKu3jOd5uLu7Y8iQIeLf6e+hpNPp0L59e8yYMQNhYWGIjIzEmjVr0KlTJ/Tp0wfl5eUm63/wwQdQqVRYunQpIiMjsWrVKjg7O8PPz8+Ge2auMfsKAG+//TY4jkNAQAAOHDiA3377Dd7e3uA4Drt27TJZV677OnHiRAwbNgzr169HVFQUwsLCMHDgQDg5OUGr1Vp8zbp163DPPfc0GMBy3dfIyEj4+/tj69atiIiIwJ49ezBlyhRwHIf3339fXC8+Ph7t27fH4MGDsXfvXuzatQuenp7w8PBATk6OyTaVvq8GgwFjxoyBu7s7Vq5cicjISPz+++949913ceDAAZNt+vr6ok2bNli1ahUiIyOxZMkSqFQqrFy50mr7QQFsJXv37gXHcfjxxx9NHvf29oaHhwcMBgOjylomLS0N7dq1w8yZM2+53oABA9CnTx/o9XrxsejoaHAch/Xr11u7zCYTBAGCIAAAcnNzGzwCzsvLM3ttWFgYOI4z+fWdl5cHFxcXzJgxw2TdlStXQqVS4cqVK9LvRCM1Zl8BoFu3bvjHP/5h8lhFRQXc3d0xbtw48TE576ulMy46nQ5du3bF8OHDzZ7jeR5ubm7YvXu3xQCW87425KmnnsJ9990n/v+kSZPQuXNnFBcXi48lJSWhdevWWLRokfiYPezrJ598ArVajRMnTtzydZcvX4ZKpUJwcLDJ435+fnB1dUV+fr5V6qUAthJfX1+4ubmZBBAA/Pjjj+A4zuz0h1KsWLECHMchKSmpwXXS0tLAcRxCQkLMntNoNPD29rZmiS12q1CyJDk5GRzHmfzj3bZtGziOM/uHn5GRAY7jrPqruiluta/du3fHmDFjzB7v2rUrJk+eLP6/Uva1vmHDhkGj0Zg97u3tjQkTJgCAxQBW4r6OGTMGPXr0AADo9Xq4urrC39/fbL0RI0agZ8+e4v8rfV+Bms/ws88+e9vXffDBB+A4DpmZmSaPHz9+3KpnLSmArWTgwIEYMGCA2eN1o7Fs2LCBQVUt989//hN33HEH/vzzTzz66KNo1aoV7rrrLvj7+4u/qP/8809wHIe9e/eavX7ixIm45557bF12kzQ1gL/77jtwHIdff/1VfGzJkiXgOA6lpaVm63fu3BlTp06VqtwWudW+fvLJJ2jVqhW++eYbFBQUICMjA//73//g4uKCY8eOiespZV/rFBUVidc969u0aRPc3d2Rnp4OwHIAK2FfjUYj9Ho9cnJy8OWXX8LJyQlff/01AODq1avgOA5ffvml2esWLFgAlUqFiooKAMrf15SUFHAch9mzZ2Pp0qXo0qULWrVqhT59+mDLli0m25kyZQruuusus+2XlpaC4zgsXbrUKvVTAFtJz549MXLkSLPH6349/v1Uh1L06tULLi4uaN++PYKDg8XrQq6urhg0aBAEQcD27dst/nIGgBkzZsDZ2ZlB5Y3XlABOS0tD165d0b9/fxiNRvFxPz8/tGnTxuJrNBoNRowYIVW5LXK7ff3666/Rpk0bcQ7UO+64AwcPHjRZRyn7Wufll1+Gk5OTSUestLQ0uLu7m/wwthTASthXf39/8e/l7Oxscsmn7jKQpc6gwcHB4DgOGRkZAJS/rydOnADHcejQoQP69OmDnTt3Yv/+/Zg4cSI4jsPGjRvFdb29vdGrVy+LbTg7O5udhpcKBbCV9OzZE6NGjTJ7vC6ALZ2eVYKePXtarP/TTz8Fx3E4ePCgGMAnT540e/2MGTMa/EctF40N4Pz8fPTt2xddunQx60Hq5+cHFxcXi6/TaDQWf5yxcKt93bx5M9q0aYP58+fj0KFD2LdvH6ZMmYK2bdvizz//FNdTyr4CwPLly8FxHNatW2fyuI+PD4YMGSJeGwcaDmC572tycjLOnDmDvXv3IiAgAGq1GqtXrwbwVwDv2LHD7HV1AVx3GtZe9tXZ2dnkkpkgCHj88cdx7733io95e3vj4YcfttiGs7OzxVP2UqAAthJ7PQU9cOBAcByHs2fPmjx+7do1cByHjz76yCFOQRcUFODxxx/HnXfeiQsXLpg9X3f6rqyszOw5uZy+Axre14KCAri6ulrsBTx06FB0795d/H+l7Gtd/4W/X7sMCwuDk5MTTp48icLCQnHhOA5+fn4oLCwUb7tRyr7WFxAQACcnJ+Tk5DTrFLTS97Vv375m6y1duhQcx4md9OgUtJ3x8/Oz2AkrNDRU0Z2wZsyYYTGA6z7sq1evFjthffjhh2av79Wrl+I7YdWFb6dOnczehzoNnQXIzMyUVQeWhva17vTd36+VAcD8+fPBcRx0Oh0AZexrXfiuWLHC7LmgoCDxNGZDyy+//AJAGfv6d5s3bxZrruuEFRAQYLbeyJEjTTph2cO+tm3b1mIA1/24yM3NBVDTs9tSJ6y6fwfUCUth9u3bZ/FUz6hRoxR9G9L+/fst/uNbu3YtOI7D0aNHAQBPPvkkvLy8TPaz7sP81Vdf2bTmprpVANeFb8eOHXHmzJkGt5Gfnw8XFxezL7qQkBBZ3cLR0L7W9ez+e/2CIGDQoEHo1KmTeLpW7vv63nvvgeM4LF++3OLzPM8jMjLSbOE4DuPHj0dkZKT4RS33fbVk2rRpUKvV4j2+kydPRpcuXVBSUiKuk5ycDGdnZyxevFh8zB72derUqWjdujV4nhfXEQQB/fr1w4MPPig+Vncb0t8PGvz9/ek2JKXy9vZGp06dsHHjRkRERMDPzw8cx2Hbtm2sS2uRsWPHok2bNnj//fdx8OBBhISEwMXFBT4+PuI6kZGRcHJywoQJE8Trwvfdd59sB+IAan40hYWFib+iJ02ahLCwMISFhaGsrAzl5eUYMGAAVCoVPvvsM5w4ccJkSUhIMNle3SAG/9/evQdFVfZxAH9c98Z1BZbrwi4qVwtImbGdTEnwgpqlMYiOkyAl1JS3sAy7ICHYhFCmko0kaqSpaHJnuShIoanNmElYNJp5qyklRwRhge/7B3M2D7vsrhrvru/7+8zsHzw85znnObuz3z3POec5a9asQX19PbKzsyGRSCw+iQFguq8A8Nxzz0EgEGD58uXQaDQoKSlBTEyM3mQHgPX2dcOGDWCMITo6Wu/9MnVvqKmJOKytr0uWLEFKSgr27t2L+vp6FBUVIS4uDowxvP7667p6LS0tsLe3x6RJk1BRUYGDBw/i0UcfNToRx8Pa119++QUjRoxAYGAg9uzZg/LycsydOxfDhg3D/v37eW1yE3FkZ2ejvr4ea9asoYk4Hma3bt3CsmXL4OHhAbFYjNDQ0P+JqSg7OjqwevVq+Pj4QCgUQqlUIjU1VS9Yq6uroVarIZVK4ezsjEWLFlntVJQAoFKpBh2CvHDhAi5cuGB0mDI+Pl6vzY0bNyIgIABisRhKpRJpaWkWn8YPMN1XoH/SjezsbISGhuqmolSr1SgsLORdrMSxxr5GREQYfc+MGSyAAevs6/bt2zFx4kTI5XIIhUKMGDECERERBqdnPHXqFKKiomBrawtHR0fMmTNH7wck52Hv6w8//IBZs2bBwcEBUqkUarUapaWlevW6u7uRlpYGpVIJsViMgIAAmoqSEEII+V9EAUwIIYRYAAUwIYQQYgEUwIQQQogFUAATQgghFkABTAghhFgABTAhhBBiARTAhBBCiAVQABNCCCEWQAFMCCGEWAAFMLEqBQUFvKkCJRIJ3N3d8dRTTyErK8vgVJbc02zuplKpMGvWrCHfXm56yoKCgn+1XW4/3D2JvLUZqr63tbXBxcXFrGlb4+PjoVKp/tX1W7uffvoJIpEI3333naU3hTwgCmBiVbjgKSgowLFjx3D06FEUFRVhxYoVkMlkcHZ2Rk1NDW+ZS5cu6U2sTwE89Iaq7ytWrEBISIjB+aYH+n8MYABISEjApEmTLL0Z5AFRABOrwgWPoUf9Xbx4ET4+PnBwcMDvv/9utB1rC2BDDzU35v81gK9fvw4bGxts3brVrPqWCuCOjg6zfiAMlVOnTj3UzxUn/SiAiVUxFsAAsG/fPjDGkJ6eriszNgRdWVmJsWPHQiqVIjAwEJ999plem9euXUNSUhIUCgVEIhF8fX2xdu1aaLVaXr0rV64gNjYW9vb2cHR0xLx583TPOL47hOLj42FnZ4czZ85g6tSpsLe3h1qtBtD/hKhnnnkGCoUCEokEo0ePRlJSku55swP3g6kAbm1tRUJCAvz8/GBjYwMvLy88/fTTOHPmDK8e93zb3bt3Y82aNfD09ISDgwOioqJw7tw5Xt2+vj5kZmZCqVRCIpEgPDwc1dXViIiIQEREhK7eYAH8888/Y8GCBXB1dYVYLEZQUBA2b95stB+cnJwciEQitLW16f2voKBA91SeoKAg7Ny502AAd3V1ISMjA4GBgRCLxZDL5UhISNB71N6dO3fw2muvwd3dHTY2Npg4cSJOnToFlUrFe7IV915oNBosXrwYcrkcjDF0dnbeU39v3ryJlJQU+Pr6QiQSwcvLC8uXL0d7ezuv3r59+zB+/Hg4OjrCxsYGI0eOxOLFi/XaCw4OxvPPP2/WfiXWiQKYWBVTAdze3o7hw4cjKipKVzZYAHt7e2PMmDHYtWsXNBoNYmNjwRhDQ0ODrt61a9fg4+MDlUqFTz/9FLW1tcjIyIBEIkFCQoKuXkdHB4KDgyGTybBp0yZoNBosW7YMSqXSYABzQb5+/XrU1dVBo9EAAD755BOsX78eJSUlaGhowM6dOxEWFobAwEDeI97MDeCGhgakpKSgqKgIDQ0N+OqrrzBnzhzY2NjwgpULYF9fXyxcuBDl5eXYs2cPlEol/P390dPTo6ubmpoKxhiSkpJQVVWFbdu2QalUwtPT02QANzc3QyaTISQkBLt27UJ1dTVSUlIgEAiwdu1ao30BgMjISIwfP16vnNsfzz77LEpLS1FYWAg/Pz/de8fp7e1FdHQ07OzskJ6ejpqaGuTn50OhUGDMmDHo6OjQ1V2wYAEEAgHefPNNVFdX46OPPoKPjw9kMpnBAFYoFEhKSkJlZSWKiorQ09Njdn9v376Nxx57DHK5HLm5uaitrcXGjRshk8kQGRmpO5puamrCsGHDMH/+fFRUVODw4cMoKCgwGLQvv/wy5HK5RY/EyYOhACZWxVQAA4C7uzuCg4N1fw8WwFKpFBcvXtSVdXZ2wtnZGcnJybqy5ORk2Nvb8+oB/zzIvbm5GUB/cDLGUFxczKu3ZMkSgwHMGMP27duN9rWvrw9arRYXL17Ua/t+h6B7enrQ3d0Nf39/rFy5UlfOBfDMmTN59bkRBe4c+o0bNyCRSBAXF8erxx3pmwrg6dOnw9vbGzdv3uQt/+qrr0IqleLGjRtGt9/W1hYvvfQSr6y3txdeXl4YN24cL2x+/fVXiEQiXgDv2bMHjDEcOHCA18bJkyfBGENeXh6A/h8KjDGsXr2aV49b3lAAL1q0SG97ze3v+vXrIRAI9D7XRUVFYIyhoqICwD+fu7///tvYbgIAbNu2DYwxtLS0mKxLrBMFMLEq5gSwm5ubWQHMDfveTa1WIzo6Wve3QqHA7NmzodVqeS/uC5r7wp43bx4cHBz02uOCzVAAD/xSBoA//vgDycnJ8Pb2hkAg4F3x/f777+vtB1MBrNVqkZmZieDgYIhEIl57d/eT286B51bPnTsHxhi+/PJLAEBFRYXBAAMAX19fowHc2dkJoVCIpUuX6u1Prl0uaAxpa2sDYwzvvvsur/zHH38EYwwbNmzQWyYiIoIXwAsXLsSIESPQ3d2ttw0eHh6YN28eACAvLw+MMb0ribVaLYRCocEAHvjj6176O2HCBISGhurVu3XrFoYNG4Y33ngDQP+IBmMM06ZNw969e3H58uVB91dxcTEYY6itrR20DrFuFMDEqvybQ9CGLsIaeB5TKBTyQmvg67333gMAREVFwc/PT6+9lpYWgwFsa2urV7e3txdhYWFwdXXFxx9/jCNHjuDEiRM4fvw4GGNIS0vT2w+mAnjp0qUQCARITU1FVVUVvv32W5w8eRJhYWG8fnIBvH//ft7yA0P0888/B2MMX3/9td66Hn/8caMBfPnyZaP7kjGGXbt2DdqXq1evgjGGrKwsXnljYyMYYygsLNRbJi4ujhfAU6ZMMbr+yMhIAEBGRgYYYwYDzt3d3WAAnzhxglfvXvrr5+dntF5iYqKu3UOHDiEyMhISiQSMMTzyyCPYvXu33nZqNBowxlBaWjroPiXWjQKYWBVTAbx3714wxpCRkaEre5AA9vDwwLRp03Dy5EmDrytXrgC49yNgOzs7vbrff/89GGPYsWMHr7y1tfW+A9jJyYl3rpqjUCjuK4CNHQGrVCqjAXz79m0MHz4cCQkJg+7Pv/76a9C+dHV1gTGGVatW8crv5Qh4/vz5cHFxGXT93Hnx+zkCHviZvJf+qtVqhISEDFrP0Pt8584d1NfXY/r06WCMoampifd/brh84C145OFBAUysijm3IclkMt4VrQ8SwC+++CK8vLxMnpu813PAhgL4zJkzYIzpTTCxatWq+w7ggee0AaCsrEzvfK25AXz9+nVIJBLdUC3H3HPAU6ZMQVhYGLq6uoxu92BGjRqFuXPn8sp6e3vh6emJ8PBwk+eACwsLwRjD8ePHja7n7NmzYIzphn45xs4BG/pMmtvfdevWwdbWFufPnzdaz5DTp0+DMYYtW7bwyjMzMyEQCEx+don1ogAmVmXgRByNjY04cOAAbyKOw4cP85Z5kAC+evUqVCoVgoKCkJeXh7q6OpSXl2PLli2YNWsWLl26BKD/aCcgIAAymQybN2+GRqPB8uXLB70K2lAAd3d3Y/To0VCpVNi9ezeqqqrwyiuvICAg4L4DeNGiRZBIJPjwww9RV1eHDz74AK6urvD29r6vAAb+uQo6OTkZVVVVyM/Ph4+PDzw9PTF58mSjyzY3N8PJyQnjx49HQUEBjhw5gpKSEuTm5vKWHUxiYiI8PT31yvPz88FY/1XQZWVlg14F3dPTgxkzZsDZ2Rnp6emorKxEbW0tduzYgfj4eBw8eFBXd8GCBRg+fDhSU1NRU1PDuwr67tt+jAWwuf1tb2/H2LFj4e3tjZycHNTU1ECj0WDbtm2IjY3V/WB45513sHjxYhQWFqK+vh6HDh3C5MmTIRKJcPbsWd66Z8+ejXHjxpncp8R6UQATq8J92XEvsVgMNzc3REREICsrS+9eTuDBAhgA/vzzTyxbtgwjR46ESCSCs7MzwsPD8dZbb/Hu0bx8+TJiYmJgb28PBwcHxMTEoKmpyewABvqHU6dOnQoHBwc4OTkhNjYWv/32230HcFtbG1544QW4ubnB1tYWTz75JBobG/X6eS8B3NfXh3Xr1sHb2xtisRihoaEoKytDWFgY7+h0sPuAL1y4gMTERN191a6urnjiiSewbt06o30BgLq6OoPnW4H+EPb394dYLEZAQAC2b99u8D5grVaLDRs2ICwsDFKpFPb29ggKCkJycjJaW1t19bj7gN3c3CCVSqFWq3Hs2DHIZDLeFeSmTouY29/29na8/fbbuvuTuduXVq5cqZtYpqysDDNmzIBCodB99mfOnInGxkZeW7du3YKtrS1ycnJM7lNivSiACSEmnT9/HmKxGJmZmUO+rpCQEL1bkf5bvvnmGzDG8MUXX1hk/ebKz8+HnZ0dDT8/5CiACSE8p0+fxurVq1FcXIwjR45g69atGDVqFORyuckpQP8NlZWVkEqluuH/oVJdXY309HSUlZWhrq4Oubm5cHNzg7+/v26WK2uk1Wrh7+9v1ogCsW4UwIQQntbWVkRFRUEul0MoFMLFxQUxMTF6U1YOpU2bNuHo0aNDuo7jx49jwoQJcHJyglAohIeHB+Lj43H16tUhXe+DOn/+PNLT0636RwIxDwUwIYQQYgEUwIQQQogFUAATQgghFkABTAghhFgABTAhhBBiARTAhBBCiAVQABNCCCEWQAFMCCGEWAAFMCGEEGIBFMCEEEKIBVAAE0IIIRZAAUwIIYRYAAUwIYQQYgEUwIQQQogFUAATQgghFkABTAghhFgABTAhhBBiARTAhBBCiAX8B23P2Jc7Rly0AAAAAElFTkSuQmCC\" width=\"480\">"
-      ],
-      "text/plain": [
-       "<IPython.core.display.HTML object>"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    },
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "\n",
-      "symmetry: 1, maximum barrier hieght: 39.10 kJ/mol\n"
-     ]
-    }
-   ],
-   "source": [
-    "if path.endswith('.yml'):\n",
-    "    results = read_yaml_file(path)\n",
-    "    units = 'degrees'\n",
-    "    energies, angles = None, None\n",
-    "else:\n",
-    "    results = None\n",
-    "    energies, angles = parse_1d_scan_energies(path=path)\n",
-    "\n",
-    "plot_1d_rotor_scan(angles=angles, energies=energies, results=results,\n",
-    "                   path=folder_path_for_saving, pivots=pivots, units=units)\n",
-    "\n",
-    "symmetry, max_e = determine_rotor_symmetry(rotor_path=path, label='label', pivots=pivots)\n",
-    "\n",
-    "print(f'\\nsymmetry: {symmetry}, maximum barrier hieght: {max_e:.2f} kJ/mol')"
-   ]
+   "outputs": [],
+   "source": "if path.endswith('.yml'):\n    results = read_yaml_file(path)\n    units = 'degrees'\n    energies, angles = None, None\nelse:\n    results = None\n    energies, angles = parse_1d_scan_energies(log_file_path=path)\n\nplot_1d_rotor_scan(angles=angles, energies=energies, results=results,\n                   path=folder_path_for_saving, pivots=pivots, units=units)\n\nsymmetry, max_e = determine_rotor_symmetry(rotor_path=path, label='label', pivots=pivots)\n\nprint(f'\\nsymmetry: {symmetry}, maximum barrier hieght: {max_e:.2f} kJ/mol')"
   },
   {
    "cell_type": "code",