diff --git a/.gitignore b/.gitignore index f48659b52618b2a50cc5e8989af6a60ae864c32d..2195ecf344b5428a6e89ec23990deaca89ad7399 100644 --- a/.gitignore +++ b/.gitignore @@ -2,3 +2,4 @@ *.pyc binary_c_python_api.o output/* +david_results/* diff --git a/david_calculations/single_star_full.ipynb b/david_calculations/single_star_full.ipynb index 30c50c50e881ec48ddea6d58d863dbb22eda7dc2..080265bf547707ecacd5a8d5b0a4e950f4412fda 100644 --- a/david_calculations/single_star_full.ipynb +++ b/david_calculations/single_star_full.ipynb @@ -81,12 +81,12 @@ }, { "cell_type": "code", - "execution_count": 4, + "execution_count": 5, "metadata": {}, "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEGCAYAAAB7DNKzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deZgd5X3g+++v6qy9qxdtLbV2gcHIYGQwxI4J18Z4ggzBG8STZGzGGu4NXjJ37sR2PLYz9/EymXFu4oHYlwAx9tgQOw4YJsIwjjHYMQaJVRKLEEJLtyRaUku9d5/tN39Une7Tre7TZ+3T5/Tv8zz1VJ23qt56j86j+vW71FuiqhhjjDGzcSpdAGOMMQubBQpjjDFZWaAwxhiTlQUKY4wxWVmgMMYYk1Wg0gUoh/b2dl27dm2li2GMMVXj6aefPqmqHTPtq8lAsXbtWnbt2lXpYhhjTNUQkUOz7bOmJ2OMMVlZoDDGGJNVTQUKEdkmIrf39/dXuijGGFMzaipQqOqDqrq9ubm50kUxxpiaUVOBwhhjTOlZoDDGGJNVTQUK66MwxpjSq6nnKFT1QeDBrVu3fqKgDPb/M3TvgmAEAhEIhCEQ9dZBfx2ITC4zHSdS2i9ljDEVVlOBomg//RycfKW4PNxwRgDJCCQTgSYj4OQakKYcO0uerv2UxpjysLvLFApv2gbXfRsS45AYm7rEp31OjEN89Oxj4xn7E6NTjxsbmOX80eKKLm5uwScY9T/762DdtPVsaRmfAxGrORmziFigmM4JQLjBW+aTKiRjpQtI088fG4BEr3/eKMRH/GPHCivvTMHkrOCTLehM2xeq95Zgeh21YGTMAmGBYqEQ8f/qD0NkHp8DSaW8IBPPXEamrTMDy/Rjph0fG4bhk2fnkRzPs2CSETzqINTgf04HlYaMANMwmZ4ONNMXC0DGFMwCxWLnOJM303JKJSdrMNODSGwE4sNekMlc4iMQG/I/+9tjAzB43E/3A1M+zXbiQKgRwplLw+R2aHp6k5/eMO2cRi+oG7MI1FSgEJFtwLaNGzdWuihmOsctX5NeKjlZmzkr2GR+9oPO+ODUZWwABo5OTUPnvq4bygggTVMDTrjRqxlGmiHSMrmOtmSkN1uwMVWhpgJF0cNjTXVy3MmbcymkUl7gSQeN2LTAMj4E4wMZ+4cmt4d64dRr3v7RM5CKZ79WIDoZNKYHkUjL2fuiS6CuDaKtXnObMfOgpgKFMSXhOBm1nxWF56PqNbWN9XtBY6w/YznjL/1T9w/1wslXJ/dpavb8A1Goa/WCRl1rxnbbtO0lk9vhRuujMXmzQGFMuYhMjvBqXJ7/+ap+v0w6kJyB0dMwcgpG+mC0z1unt4/v9rdPM2vTmRP0aiX17d7SsAzql0JDh79eCvUdk2s3WNQ/gakNFiiMWahEJpvUmlflfl4q6QWXkT4vqEwElGnbQ73QvdNbx0dmzivaOi14LIXGZdDUCU0rvaVxpff8jqlZFiiMqTWOO9kURY4DO8aHYLgXhk74a39Jbw+fgJ5n/KAyfPb5dW1ewEgHj8xA0tTpBTrrU6laFiiMMZN9Mq3r5z52bAAGj3kjxSaWHj+tB3p2eTWW6eqXwpI1sGQttKzxtlv8z02dNg3NAma/jDEmP5Emb+k4Z/Zj4mMZwaQHzhyC04e89ZEnYc8/giYnjxfXq3UsWQvtm6B9s7du2+QFEaemJrquOgs+UIjIm4BPA+3AP6vqtypcJGPMXIIRaF3nLTNJxr0AcvrgZAA5fRD6XocXfugNL57Iqw7aNvrBww8gy86H1g1WC5knFflXFpG7gGuAXlV9c0b61cBfAy5wh6p+XVVfAm4WEQf4LmCBwphq5wa92sOStWfvU/WHCe/zl1e9dfdTsOfHTIzoCkRg6Ztg2Zth+QXesuz8+Z0CJ63vdW89W2CscpUKx98BbsW78QMgIi5wG/AeoBvYKSIPqOqLIvJ+4P8EvleBshpj5pOIN7KqcRmse+fUfbEROLUf3tjrDQd+Yze8/E/wbMatoWUNdL4VVr0NOrfCireUf1TWNy/01l+uzZemVSRQqOrjIrJ2WvIlwH5VPQAgIvcC1wIvquoDwAMi8k/AD2bKU0S2A9sBurq6ylRyY0xFhepgxRZv4UYvTdXrDzm+B46/4C3du2Dvfd5+J+DVNjq3wupLYM1vQXNnxb5CNVpIDXydwJGMz93ApSJyBXA9EAZ2zHayqt4uIseAbaFQ6OJyFtQYs4CITA7F3XzVZPrgcS9g9Ozy1s/9AHb+rbevdT2sfae/vAOaingCfxFYSIFiRqr6C+AXOR5rcz0ZYzyNy+FN13gLeA8ivrEHDv4KXv8l7L0fnrnb29e2ETZcCZve6wUOe4BwioUUKHqA1RmfV/lpObPZY40xs3Jcr79ixVvgsj/2AsfxF/zA8Tg88z146nZvDq3174JNV8Hm9879VHxshgcQK+DZw6d59vAZPvr2LsIBt6R5L6RAsRPYJCLr8ALEDcDv55OB1SiMMTlzXFh5kbdc/knv3SgHfwX7HoZXH4Z9P4V/wuvbOP86OO86aFl9dj7P3zPvRZ/JY/tO8Fc/e5U/uGxNyfOu1PDYe4ArgHYR6Qa+pKp3isgtwMN4w2PvUtW9eeZrNQpjTGGCUdj0Hm/R/+oNyX1lh9dE9cgXvGXV27yAcf51Xk0jlYLffNs7f9mbs+dfZqeHYzRHgwTd0j+cWKlRTzfOkr6DLB3WOeRrNQpjTPFEvCfPO86Bd/wJ9B3wAsbe++CRP/OWVZd4o69Oveqdk21K+HlwajhGa32oLHkvpKYnY4xZmFrXwzv/vbeces0LGHvvh513eJMhrrzQCyYVdHokxpK68kwLX1OBwpqejDFl17YBfvs/eMup17ynzP/XF73O8Qo63j/GOctL9JbHaWpqpi1VfVBVtzc3V+ARfmPM4tO2AVq6vEkNK9j0pKocPTPGyuZoWfKvqUBhjDEVIc7U2XDn2emROKPxJCtbLFDMSUS2icjt/f21Od+KMWaBcipbozh6ZhTAAkUurOnJGFMR4nhzTlVIjx8oVi2xQGGMMQuTSEU7s3tOW40iZ9b0ZIypCHE5PTLGV3e8VJHLHz0zSiTolG14bE0FCmt6MsZUhDgkEkluf/wAvzkww/vCy+xo/ygrW6KISFnyr6lAYYwxFSEODl5n9qfvfXaiz2C+9JwZo7NMzU5ggcIYY4rnuDgo79rcwUgsyR/d9RS9A2Pzdvme06MWKIwxZkETBwflzZ1N3PGHWzl6ZpRrb/sXnnq9r+yXHosnOTk0XraObKixQGGd2caYivCbnlwRLl3fxo9uvgzXET78/z/Bzd97msf2nWAsXp5RUcf7vZpLOQNFTc31ZLPHGmMqQfEDheP97X3+ymYe+ZPf5lu/eI3vPnGIn+49TtAVOluitNaHCDgOdWGXr/zeBUU3GaX7Q8rZ9FRTgaJYn4/EeGTkGQI/eDuCN3ogvZ5cycTIgunHZI44mJ425ViZ/ZiZ8pp+7Gx5z3q96WXLckxmXnOVbcZ/m3zKNi1vEcERZ6Icjjg4OJPpGcfMmp75OeOY6enTr5OZL8KUY9LnTllnljF9nH9Oegk4AW8t3tp1XFxxJ9fTts86x3HmPM4Vt2wjXUzu/rDvV8RXLuEdGW00daEA//dV5/B/XbGR3xw4xc6DfRzuG6F/NM5YPMkvXjnBL17p5aOXFveiIQsU82y3q4yjfGjj701JV7wnLlV1ynbmvpnSph87U9pMeWfL76zrTqzyLNtMZc3Ia9bvPcN1czlmpuuq6lmfU6S8taYm0zQ1e3rmZ/zP04/LSFd0yjEpUqBMHFON0gEjHVAccQg6QQJOgKATnHnbnX1/5jHZ8gm6QSJuhJAbIuJGCAfCs34OOLV9q3ku3geREO9yzm7Nj4ZcfufcpfzOuUsn0lIp5fwvPcyBE8W/RrXn9CgisKw5XHRes6ntXy9PovBedwl/esmfVroopkJyCTITgYaz09NLUpPekkqS0hQJTZBKeemJVGIyTVMkU8nJ4zPP8Y/LTE9vZ+6f6bhEKkFCE8STceKpOIlUgnhqcjuWjDEcG55x3/Tjpv8BU4iABAgHwoTdySUSiExsRwNR6oJ11AfrqQ/UT2zXBeuoC/jpwXrqAnWT+/z0hVSjyvXlco4jrGuv57UTQ0Vf80jfCMubIiV/T3YmCxTGZBARXCnff7hqlEwlZw0i48lxxpJjxJIxxhJjOX2eWBLevvHkOAMjAwzHhxmJjzCSGGE0kdtzCK64NIWaaAo30RxqpjHcSHOomaZQE83hZprDk9utkVbaom20RdqIBCJl+bdyZ6hRzGZ9Rz3Pd58p+pqH+kboaq0rOp9saipQ2IuLjCk91/GatOZTMpVkNDHqBY/ECCPxkYntiXVsmIHYgLeMD9Af66d/rJ/DA4fpH+9nMDY4a22oMdjoBQ0/cLRH22mPtrOsfhkr6lewsmElS+uWEnTymxIj4OReu1nfXs+O3ceIJVKEAoUPQD10aoT/I6NZqxxqKlDYqCdjaoPruDSEGmgINRScR0pTDMYGJ4JI31gfJ0dPcmr0lLce89b7Tu/jiaNPMBgfnHK+Iw5L65aysn4ly+uX09nQyerG1axpWsOapjW0RlrPavbKI06wtr2elMLhvhE2Li3sew6NJzg5NE5Xm9UojDEmb444E81Pq1k95/FjiTGODx/n6PBRjg0d4+jwUe/z0FGeP/E8jxx8hIQmJo5vDDbS1dRFR7RjIi2fisHa9noADp4cLjhQHD41AsAaCxTGGFN+kUCEtc1rWdu8dsb9iVSCY0PHODhwkEMDhzg0cIiDAwc5PnJ84phAHjWKdW1+oDhV+Minw33euWta6wvOIxcWKIwxJgcBJ8DqptWsblrNO3nnlH1//sMP8PDQy7iS+wixJfUhmqPBogLFQb9GUe6mp5qawsMYYyrBxSEp+dUowGt+ev1k4YFif+8QHY1hmqPleQ9FmgUKY4wpkiMOCSTn5yjSNnY08OobhT9L8crxQc5Z1ljw+bmyQGGMMUWarFHk93Dim1Y00js4zqmh8byvmUwpr/YOcs5yCxSIyHUi8rci8vciclWly2OMMdM5OCRFcCS/aWDOXd4EeDWDfB3uG2EsnqrdGoWI3CUivSKyZ1r61SLyiojsF5HPAqjq/ar6CeBm4COVKK8xxmTj+LdSR/KbSvzcFd5N/qUCAkU6uNRyjeI7wNWZCSLiArcB7wPOA24UkfMyDvmCv98YYxYUR7xbqZJfoGhvCNPeEGZvT/7v0NnT04/rCJtrtUahqo8D01/9dAmwX1UPqGoMuBe4Vjz/BXhIVZ+ZLU8R2S4iu0Rk14kTJ8pXeGOMmSZdoxDieZ97UVcLzx7Jf86nZ4+c5k0rGomGyj+9ykLqo+gEjmR87vbTPgm8G/igiNw828mqeruqblXVrR0dHbMdZowxJefg3awlzxoFwMVrlvD6yeG8OrSTKeX5I/1ctHpJ3tcrxEIKFDNS1W+q6sWqerOqfjvbsfYqVGNMJQjppqfEHEee7eI13s3+mcO51yr2vTHI0HiCC1e35H29QiykQNEDUyZkWeWn5UxVH1TV7c3NzSUtmDHGZDPZ9JR/jeKCzmZCAYcnXjuV8zm/fNVrXr98Y1ve1yvEQgoUO4FNIrJORELADcAD+WRgNQpjTCU4RdQoIkGXy9a38egrvTmf8/i+k2xe1sCK5vK9/jRTpYbH3gM8AZwjIt0icpOqJoBbgIeBl4AfqurefPK1GoUxphKK6aMAuPLcpbx+cjinN971j8Z56vU+3rV5/vpiKzXq6UZVXaGqQVVdpap3+uk7VHWzqm5Q1a/km6/VKIwxlSB4kzylCqhRAFx1/jIcgfufnbu1fcfuY8SSKa7ZsrKgaxViITU9Fc1qFMaYSii2RrGiOcq7Nnfww11HiCdnf7pbVbn3qcNsXNrAllXzd5+rqUBhNQpjTGUU9sBdpj+8bC1vDIxzz1OHZz3m0Vd6eb67n4//1rqz3q5XTjUVKKxGYYypBEf9QKGFB4orzungsvVt/LeHX+HADH0VJwbH+cJ9e1jXXs8HL15V8HUKYS8uMsaYIonf9FTIqKeJPET4iw9u4f23PcYH7vwRH7y4i4tXnIuIS/fpUe761eucHolz7/aLCeXzztUSqKlAISLbgG0bN26sdFGMMYuIlKDpCeDI6HM0bPw6Z2J93HsU7jkSIX7mbcT63sl5S1dx6+9fxFvm6SG7TNb0ZIwxRZKJpqfCaxR7Tu7hkz//JB31bfzXS7/If1pzI5d3XESk/de0nfuXfPQ9R7ioa36m7JiupmoUxhhTGV7TU6rAPoqUpvjPT/xnlkSWcOeGf82SH98C8WE+DHSv/22+2tLE1576GpFAhOs3XV/CcuempmoUxhhTCcXM9QTws0M/46W+l/j0+R9jyQOfgbb18Af3w+/8GauOv8g3n36IyxrX89Unv8qBMwdKWfSc1FSgsOGxxpjKKK6P4kf7fkRnQye/2/0KxAbhg38HG34H3vUf4VPPElhzOV996TeEJMA3nv5GKQuek5oKFNZHYYypCPWanpIFND0dGzrGk8ee5Nr123Ce+z6c+7vQvmnygEgTfOhu2p0wH0/V8Xj347x46sVSlTwncwYKEblMRG4TkRdE5ISIHBaRHSLyxyJid2RjjKHwzuyfH/k5ivK7kZUwchLe/IGzD6pvg8tv4cOvP0PECfGjfT8qtsB5yRooROQh4N/iTdR3NbAC7zWlXwAiwE9E5P3lLqQxxixkqVT6Vjr79Buz+Zeef6GrsYuuw7vADcHGd8984IUfpUmFq6KdPPT6Q8SSscILnKe5ahR/oKo3qeoDqnpUVROqOqSqz6jqN1T1CuDX81DOnFgfhTGmEtSfFFDQvM4bT46z8/hOLl95ORx5EjovhvAs78Bu7oSuy7iq7w2G48PsOr6r2GLnLGugUNWTc2WQyzHzxfoojDGVkCI971J+geL53ucZS47xW8svgWPPw6qt2U/YcCWXHnuFiBvmF92/KKishZir6WlQRAYylsHM9XwV0hhjFrLUxAN3+TU9PX/ieQAuSgUhGYNVl2Q/YcOVRFS5tGENvz46f405WR+4U9VZ6kDGGGPSUn5FIt9A8cLJF1jbtJbm3pe9hLlqFMsvADfMW5Mujw0d4tToKdqi5X8das7DY0XkLSJyi79sKWehjDGmmkzUKPJoelJVdp/YzQXtF8CJlyHSAo0rsp/kBmH5BVw04L1f+7kTzxVc5nzkFChE5NPA94Gl/vJ9EflkOQtmjDHVIl2jSKVyf47i+PBxTo2d4oKOC+DkPug4B3J5x0TnWznv6EsEnSDP9S6gQAHcBFyqql9U1S8Cbwc+Ub5iGWNM9UilvBt8PjWK1wdeB2BTyyYvULRvzu3E5VsIx4fZ3LiGl/tezrushcg1UAhMeTY96actKDY81hhTCan0A3ep3Psojg8fB2CFG4XhE16NIhd+QNkcamXf6X35FbRAuQaKvwOeFJEvi8ifA78B7ixfsQpjw2ONMZWQUvHXuTc9jcRHAGgY8AIG7bkGCm96j80aoG+sj5Oj5X9CIadpxlX1L0XkF8A7/KSPqeqzZSuVMcZUkaQKiDddeK7iqTgAwTPdXkLbhtxOrGuFunY2jXmBZt/pfbRH2/Mqb75y7czeAOxV1W8Cu4F3isj8v2bJGGMWoHSNIp8+iolAMejXKJo6c79g+yY2978BwKunX839vALl2vT0YyApIhuBbwOrgR+UrVTGGFNFkv7wWDT/QBHoPwL1SyEYyf2CbRtY0neIlnALhwYO5VPUguQaKFLqTYt4PXCrqv4/eBMEGmPMopeaiA95ND0l4wSdIDLQA82r8rtgcxcMHaerYRWHBw/nd24Bcg0UcRG5EfhD4H/6acHyFMkYY6pLYuLJ7PxqFEEnCGeOQMvq/C7oB5bVkTaODBzJ79wC5BooPgZcBnxFVV8XkXXA98pXrEkisl5E7hSRf5iP6xljTK76hmMMjMUnphkvKFD0d0NzYYGiy63n2PCxsk85nuuopxeBT2V8fl1EflXoRUXkLuAaoFdV35yRfjXw13hvKr9DVb+uqgeAmyxQGGMq6cxIjD09A+zu6Wd3zxle6O6n+/Qoq1ujXBZNQiS/uZ68QOFCYrTgQLEaF0XpHupmffP6/PLIQ9ZAISIu8GGgE/ipqu4RkWuAzwNR4KICr/sd4Fbgu9OudRvwHqAb2CkiD/hByhhj5s2poXF29/Sz9+gAe3r62XO0nyN9oxP7u1rreMvqFq7ZspJvP/Yaze4wNOU56ikZJ5h+bjnvPgo/UIx7NYmewZ7KBQq8h+pWA08B3xSRo8BW4LOqen+hF1XVx0Vk7bTkS4D9fg0CEbkXuBawQGGMKYvh8QSvnRhi3xtDvNo7yKtvDPHSsQGO9Y9NHLOmrY4tnS38/iVrOH9lE1tWNdNSF5rYf2YkxgvP+FNp5FujSMeV5jyGxgIEwtCwjGWjZwDoHenN7/w8zRUotgJbVDUlIhHgOLBBVU+VoSydQGavTDdwqYi0AV8BLhKRz6nq12Y6WUS2A9sBurq6ylA8Y0w1UlV6B8c5eHKYQ30jvNY7xL43Bnm1d4ju05O1hJDrsL6jnkvWtXJBZzPnr2zmvJVNNEezj9v57PvO5cZnfuZdK8/nKILp4xuW5//FmlfRPnACQXhj5I38z8/DXIEipn6jm6qOiciBMgWJWfnXuzmH424XkWPAtlAodHH5S2aMWShiiRTH+kc53DfCwVMjHD417K9HONQ3zFh88i/9UMBhQ0cDb+1awke2rmbTskY2LWtgTWsdATfnNy9MaKkL8Vc3vJUPPvN3+Xdmp8fV1hfwZHXDcoKnD9LW3lbxGsW5IvKCvy3ABv+zAKqqpXwvRQ9eM1faKj8tZ6r6IPDg1q1bbWZbY2pEKqWcHB7n2Jkxjp4ZpefMKMf6ve2j/vrk0PiUZ93CAYc1bXWsaavnnZvaWdNez5rWOta01dHZEi0oIGSzekkDUEBntiahrs17z0S+GpZC91Ms7bqY4yPH8z8/D3MFijeV9epT7QQ2+UNve4AbgN/PJwMR2QZs27hxYxmKZ4wppVgixcmhcU4MjtM76K297TFve2ic3gEvLZacegOOBl1WtkRY2RLl3HOWssLf7mqtY21bPUsbwzjO/E1wLU7+r0JNJBMEkgnvqexCNCyD4ZMsiy6lezivv6nzNlegOKxz1KVEROY6ZoZz7gGuANpFpBv4kqreKSK3AA/jDY+9S1X35pOv1SiMqZyxeJLTIzH6hmOcHo7TNxLj9LD/eWRynQ4MZ0biM+bTWh+ioyFMR2OYS9fV09EUprMlyormqBccmqO01AWRXF7yM0+8QZv59VGkSCHJODTkOeIpraEDUJaFmni695nC8sjRXIHiURH5MfATVZ14TlxEQngzyf4R8CjecNecqeqNs6TvAHbkk1cmq1EYUxxVZTiWZGA0Tv9onIHROANjicnPY3HOjMSn3PhPD8fpG44xGp99iu2WuiCtdSGW1IdY397Apeva6GgMs7TRCwjppb0hTLDEzULzQST/uZ5UFScZh6ZlhV20wTtvmRNhIDbAaGKUaCBaWF5zmCtQXA18HLjHbxI6A0Tw/uJ/BPirhTTduNUozGIXS6QYHk8w5C/D4wkGxxIMjM1+4x8YTWRsxzPmLZpZYzhAa0OIJXUhljZGOGdZE631QVrqQrTWe+mt9SFa64MsqQvRHA2WvE9goRHyf2d2SlO4yZjX11AIP1As9W/jvSO9rGlaU1hec8gaKFR1DPgb4G9EJAi0A6OqeqYspTFmERpPJBkaSzA8npxyg59+sx8eTzAcm9z29ienHBtLzN1GHg26NEUDNEWCNEWDtDeE2NBRT1M06KcFaJ7Y9tbNUS+9IRyo+Zt+ISb7KPJohU8lcTRVRKDwzmuIe898DMeHC8snBzlN4QGgqnHgWNlKUgLW9GTKRVUZi6cYiSUYiSUZjXs36NFYkpFYkpF4kpHxqftGYklvf8a+kXiS0ZgXFEZi3s09nszt5lIXcqkPB2gMB6gPB6gPu3S2RGkIuzREAtP2TW57N3tv3RgJEA64Zf7XWoQk/87sVHLc2yi0M9s/LzI+BMB4Or8yyDlQVANrelqcVJVYMsVYLMVYwrs5T6zjKcbiScbi3g18LJ7y194y4t/oR/0A4C2Z2/6+eDKf5mdcR6gLuf4SIBr0tpujQVY0RSZu+vXhAI2RAPUhl4ZIkIawl96QXvwAUB8K4M7jKB6Tn3QfRT5NT5pK4KAQLfAdcKE6CNYRjXk1idHE6BwnFK6mAoVZOJIpJZZIMZ5IMp5ITb1R+zfysYkb+uT+8Vlu6JmfveOmfs5v3J1HBOqCLtGQ99d5+mZeFwrQ3hCmLuTtqwu51GdsR/0AUB8KTGynz0vvD7nOghqVY8qrkOGxqklEATc057GzirYS8gPFeKLCNQoRqcfrm0iJyGbgXOAhvzlqwbCmJ0/6L+zxRIrxeMrbjns3bC9tcjvzZj57esZ2ln2xdP6JZM7NKdOJQCTg3WwjAYdIyJ38HHRYUhckHPRu6pGg46+9Jb0dDTlEAu5Z52YeGwk6djM3JZMeHks+NQpVb0rAQh62S6tbQmBsAICkzj7qrFi51igex3tP9hK80U47gY8AHy1XwQqxUJqeEumbtH/TnLiBTrnhJjNu4lNvxOPxJOMT6WfflCdvyJP5eDfqyTyK5QhEgi7hgEMo4BAOeNvh4OR2Qzgw677p6Zk37Ck39Yy1d77dvE31mXiOIo+qbUqTfqAookZR10ZwdACCkNBE4fnMIddAIao6IiI3AX+jqn8hIs+VrVQVoiiDY3H+y09fnvyredrNOjbDDXr6vuRc4wtzEA44/k12hhtvwKEpGsxIn/lmHc7cF5zcDs2Z7tjIFmPyIAV0ZquqN6i2mBpFtBV38AgEIZmqfI1CROQyvBrETX7aghs6UWzTUzKlDI4luOOXB6bebKfcrL0mjJZocOKmHHKdKTfc9A3bS3dnvJFn3pQnj/GuYU0ixlSZQjqzNeXVKJximp5acUfPQEP9gmh6+gzwOeA+Vd0rIuvxnsheUErR9BQOOLz6lX9Vwk6NruYAABInSURBVFIZY2qeOIhqXk1PqilEtejO7MBoP1Bf+RqFqj4GPJbx+QAZr0Y1xphFTSQ9pXbOp6TSNYqiOrNbcfGau+Kp8o0tynXU06PM0J2vqleWvETGGFNtxPFfappHjYJSjHpqw/UvuRCanv5DxnYE+ABQvi52Y4ypJn6gyLePwuvMLq7pyfWvuRCanp6elvQvIvJUGcpTFHuOwhhTUXnOHlt0Z3a0hcA81ChyGgMpIq0ZS7uIvBdoLlupCqSqD6rq9ubmBVc0Y0wt80cp5tWZTbozu4hAEW6aGH6aSFX+OYqn8RrfBK/J6XUmh8kaY8wi53dm53HGZGd2EU1PkSZcPzhVvI9CVdeVrQTGGFMDRPPro6AUD9yFmyZu4hWrUYjIlar6cxG5fqb9qvqP5SmWMcZUEREEJZ86RSrdSOMU8exyMIo4AVykok1P7wJ+DmybYZ8CFiiMMYb8n6NQVRwpcqocEb+fQio315Oqfslff6xsJSghG/VkjKkImb2PYiQ+wpPHnuTl0y9zfPg4KU0RdIKcTo2XZqqeiNehnUoVPxnobOZqevr32far6l+WtjjFWSizxxpjFqfMPorB2CC3v3A7f//K30+8VKg92o4rLvFUnFFNsrYU/c/hJgIMVHT22EZ/fQ7wNuAB//M2YME9R2GMMZUxtelp76m9fObRz9A70sv71r2P6zdez5aOLUQCkYkzUj/5JM6Jh4u/dKSZgA5Uro9CVf8cQEQeB96qqoP+5y8D/1S2UhljTDWZaHpSnut9jpt/djPNoWa+977vsaVjy4ynOKlEcUNj08JNuOO6IJ6jWAbEMj7H/DRjjDF+jeJYYohP/fxTtEfbueOqO1hev3z2U5Kx4obGpoUbCYxp5Z+jAL4LPCUi9/mfrwPuLk+RjDGmyog3OPafhw/REGzgv1/537MHCYBUvLjpO9IiTbinU5WvUajqV0Tkp8A7/KSPqeqzZSuVMcZUmWHHG+r6+Us/z7rmHJ5RTsZL1vQU1CTJSgcK8CYGFJEjeLPHIiJdqnq4bCUzxpiqITQmU0gwwjXrr8ntlFI1PfnTeCQS48XnNYtcJwV8v4i8ijfH02P++qGylWrqtetF5G4R+VsR+eh8XNMYY/Iiwk96jvLz1R/O/dmIZLxEfRRNBBQSybHi85pFro8F/r/A24F9/rxP7wZ+U+hFReQuEekVkT3T0q8WkVdEZL+IfNZPvh74B1X9BPD+Qq9pjDHlI3QkU4TzmY6jVE1PkSZclESi8oEirqqnAEdEHFV9FNhaxHW/A1ydmSAiLnAb8D7gPOBGETkPWAUc8Q8rX7f+ZEnKfwljTI3KY1LAVBycnFv/ZxesI6AsiD6KMyLSADwOfF9EeoHhQi+qqo+LyNppyZcA+/33cSMi9wLXAt14weI5sgQ2EdkObAfo6uoqtGjGGJO/dHNTPvOMJ2OlqVG4QQIoiWT53pmda43iWmAE+BPgp8BrzDxRYDE6maw5gBcgOvEmHvyAiHwLeHC2k1X1dlXdqqpbOzo6Slw0Y4zJJt0SkUekKFUfhRvy+igqOIUHAKqarj2kgLtFxAFuBL5froJNu3ZOkxLapIDGmIooZHK/UgUKJ4iLMlbGpqesNQoRaRKRz4nIrSJylXhuAQ4AHy5xWXqA1RmfV/lpxhhTHfKYZrx0TU8hv4+icu/M/h7ehIC7gX8LPAp8CLhOVa8tcVl2AptEZJ2IhIAbmJyEMCf2zmxjTGVUsukp6D1HUcGmp/WqegGAiNwBHAO6VLWocVgicg9wBdAuIt3Al1T1Tr+28jDgAnep6t5irmOMMfNiojM731FPJeqjoLw1irkCxUQ3uqomRaS72CDh53XjLOk7gB2F5mt9FMaYipBCahQlanpyAgTKXKOYq+npLSIy4C+DwJb0togMlK1UBbKmJ2NM1UgmStSZ7eACCa3QG+5UtYi3fs8/q1EYYyoq787sUgSKdI2icp3ZVcVqFMaYypntrdkzUPX6KErR9CSu10dRxhpFTQUKY4ypGJHcaxTpZx5K0ZntuH6NwgJFTkRkm4jc3t/fX+miGGMWnTweukv6LwwtRdOTuLha3j6KmgoU1vRkjKmsHGsUE4GiFKOeXG+uJ+ujMMaYBS6fpqek3/RUkhqFQ1MqxagmiJdpYsCaChTW9GSMqZw8OrNL2fTkuLQmvdpE31hf8fnNdImy5Foh1vRkjKmYfCYGLGnTU4DWpNc/YYHCGGMWukqMehKXtjLXKErweiVjjDGzNj2pwv6fwXM/gJ5dMHpm8s12bgluwY5b9hqFBQpjjCmFmTqzx/rhvpvhlR1QvxTWvgMalkJ8BBIx6Lq8BNd1yt5HUVOBoiRTeNgrs40xBZlWoxgfgrvfD2/sgau+Apdsh0AJ+iTOuqxQjxDC4dToqdLnT431UVhntjGmYjI7s1XhJ38Mx3fDR74Pl99SniAxcWmXVjfMqTELFMYYs7Clm5723gcv3g9X/hmcc3X5r+u4tDphG/VkjDELm1+jSIzDI/8JVlwIl396ni7tskSC9I+X5xkyCxTGGFMK6c7sp78DA93w7i+XZlRTLpwAUXEYTYyWJ/uy5GqMMYuOgCbh17dC12Ww/or5u7TjEMECRU5sCg9jTEXt/xn0H4ZL/11+T2oXS1wiOIwlin5T9YxqKlDYqCdjTMWIwKn90LAMzr1mfq/tuEQRxpIWKIwxZuEaH/DWF3yoNJP95UNcIgijiVE0n9ex5sgChTHGlNJ5183/NR2XKA4pTRFPlX6qcQsUxhhTSp0Xz/81xeFDwQ52XL+DgFP6kVY1NYWHMcZU1OarwanA399OgBZ1aGlcXZbsLVAYY0wpfPH0/I50yuS4k1OXl4EFCmOMKYVK1CTSxPWe4SiTBd9HISLrReROEfmHSpfFGGMWJMeFVKp82ZctZ0BE7hKRXhHZMy39ahF5RUT2i8hns+WhqgdU9aZyltMYY6qaOGWtUZS76ek7wK3Ad9MJIuICtwHvAbqBnSLyAOACX5t2/sdVtbfMZTTGmOrmBKq3j0JVHxeRtdOSLwH2q+oBABG5F7hWVb8GFPw4o4hsB7YDdHV1FZqNMcZUnzIHikr0UXQCRzI+d/tpMxKRNhH5NnCRiHxutuNU9XZV3aqqWzs6OkpXWmOMWejcICSrtEZRCqp6Crg5l2NL8SpUexOqMabqOAEo04SAUJkaRQ+Q+VTIKj+taDYpoDFmUXKDkIyVLftKBIqdwCYRWSciIeAG4IFSZGzTjBtjFiWnvE1P5R4eew/wBHCOiHSLyE2qmgBuAR4GXgJ+qKp7S3E9q1EYYxYlNwBlmAwwrdyjnm6cJX0HsKPU1ytFH4UxxlQdJwjJ8gWKBf9kdj6sRmGMWZTcUFlrFDUVKKyPwhizKLmB6u2jmG9WozDGLEpO0GoUxhhjsnCtjyJn1vRkjFmUglGIj5Qt+5oKFNb0ZIxZlEKN3gN3ifI8dFdTgcIYYxalcIO3jg2VJXsLFMYYU+1CfqAYHyxL9jUVKKyPwhizKFmNInfWR2GMWZRCjd563AKFMcaYmUzUKKzpyRhjzEwm+iisRjEn66MwxixK1keRO+ujMMYsStZHYYwxJivrozDGGJNVIOxNDGg1CmOMMbMKN1gfhTHGmCxCjVajMMYYk4XVKHJjw2ONMYtWqMHmespFscNjtcTlMcaYeROMQGK8LFnXVKAwxphFyw1DYqwsWVugMMaYWhAIey8vKgMLFMYYUwsCYWt6MsYYk4VrNQpjjDHZBEJlq1EEypJrCYnIdcDvAk3Anar6SIWLZIwxC48bhmQVNj2JyF0i0isie6alXy0ir4jIfhH5bLY8VPV+Vf0EcDPwkXKW1xhjqlYV1yi+A9wKfDedICIucBvwHqAb2CkiDwAu8LVp539cVXv97S/45xljjJnOLV9ndlkDhao+LiJrpyVfAuxX1QMAInIvcK2qfg24ZnoeIiLA14GHVPWZcpbXGGOqViAMKKSS4LglzboSndmdwJGMz91+2mw+Cbwb+KCI3DzbQSKyXUR2iciuEydOlKakxhhTLTrOhfN/zwsUJbbgO7NV9ZvAN3M47nYROQZsC4VCF5e/ZMYYs4Cc935vKYNK1Ch6gNUZn1f5aUWzV6EaY0zpVSJQ7AQ2icg6EQkBNwAPlCJjmz3WGGNKr9zDY+8BngDOEZFuEblJVRPALcDDwEvAD1V1bymuZzUKY4wpvXKPerpxlvQdwI5SX09EtgHbNm7cWOqsjTFm0aqpKTysRmGMMaVXU4HCGGNM6dVUoLDObGOMKb2aChTW9GSMMaUnqrX3pmgROQEcKvD0duBkCYuzkNTyd4Pa/n61/N2gtr9ftXy3NaraMdOOmgwUxRCRXaq6tdLlKIda/m5Q29+vlr8b1Pb3q4XvVlNNT8YYY0rPAoUxxpisLFCc7fZKF6CMavm7QW1/v1r+blDb36/qv5v1URhjjMnKahTGGGOyskBhjDEmKwsUPhG5WkReEZH9IvLZSpen1ETkoIjsFpHnRGRXpctTLBG5S0R6RWRPRlqriPwvEXnVXy+pZBkLNct3+7KI9Pi/33Mi8q8qWcZCichqEXlURF4Ukb0i8mk/vep/uyzfrep/O+ujAETEBfYB78F7NetO4EZVfbGiBSshETkIbFXVanjwZ04i8tvAEPBdVX2zn/YXQJ+qft0P9ktU9U8rWc5CzPLdvgwMqep/q2TZiiUiK4AVqvqMiDQCTwPXAf+GKv/tsny3D1Plv53VKDyXAPtV9YCqxoB7gWsrXCaThao+DvRNS74WuNvfvhvvP2nVmeW71QRVPaaqz/jbg3jvpOmkBn67LN+t6lmg8HQCRzI+d1MjP3AGBR4RkadFZHulC1Mmy1T1mL99HFhWycKUwS0i8oLfNFV1TTPTicha4CLgSWrst5v23aDKfzsLFIvHO1T1rcD7gD/2mzdqlnptqrXUrvotYANwIXAM+EZli1McEWkAfgx8RlUHMvdV+283w3er+t/OAoWnB1id8XmVn1YzVLXHX/cC9+E1t9WaN/x24nR7cW+Fy1MyqvqGqiZVNQX8LVX8+4lIEO9G+n1V/Uc/uSZ+u5m+Wy38dhYoPDuBTSKyTkRCwA3AAxUuU8mISL3fuYaI1ANXAXuyn1WVHgD+yN/+I+AnFSxLSaVvor7fo0p/PxER4E7gJVX9y4xdVf/bzfbdauG3s1FPPn/I2l8BLnCXqn6lwkUqGRFZj1eLAO896T+o9u8nIvcAV+BN4fwG8CXgfuCHQBfeNPMfVtWq6xSe5btdgdd0ocBB4N9ltOlXDRF5B/BLYDeQ8pM/j9eWX9W/XZbvdiNV/ttZoDDGGJOVNT0ZY4zJygKFMcaYrCxQGGOMycoChTHGmKwsUBhjjMnKAoVZ9ESkLWNmz+PTZvr8dZmueZGI3Olv/xsRURF5d8b+6/y0D+aZ770isqnU5TWLmwUKs+ip6ilVvVBVLwS+Dfx/6c+qenmZLvt54JsZn3fjPeiZdiPwfD4Z+rMgfwv4j0WXzpgMFiiMyUJEhvz1FSLymIj8REQOiMjXReSjIvKU/56PDf5xHSLyYxHZ6S+/NUOejcAWVc0MBL8ELhGRoD9X0EbgOf/4K0Xk/ozz3yMi96XLJyLfEJHngcv8fN4tIoHy/IuYxcgChTG5ewtwM/Am4A+Azap6CXAH8En/mL/Gq5G8DfiAv2+6rZw9jYMCPwPeizflduYUMo8C54pIh//5Y8Bd/nY98KSqvkVVf+XPJ7TfL6sxJWGBwpjc7fTfOTAOvAY84qfvBtb62+8GbhWR5/Bu9k1+DSHTCuDEDPnfi9f8dANwTzrRn031e8C/FpEWvJrDQ/7uJN4kdJl6gZV5fztjZmHVU2NyN56xncr4nGLy/5IDvF1Vx7LkMwpEpieq6lMicgEwoqr7vDnmJvwd8CAwBvxIVRN++piqJqdlFfGvYUxJWI3CmNJ6hMlmKETkwhmOeQmvD2Imn8Xr6J5CVY8CR4Ev4AWNbDZThTOUmoXLAoUxpfUpYKv/NrMX8fo0plDVl4Hm9NTv0/Y9pKqPzpL394EjqvrSbBcXkWXAqKoeL6z4xpzNZo81pgJE5E+AQVWdqbN7tnNuBZ5V1TvnyHcg2zHG5MtqFMZUxreY2ueRlYg8DWwB/scch54B7i6iXMacxWoUxhhjsrIahTHGmKwsUBhjjMnKAoUxxpisLFAYY4zJygKFMcaYrP437BUn7ZQo6w4AAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEGCAYAAAB7DNKzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deZwcZbnw/d/V2/TsySzZl0lIQlgSCIQEkAiisgghrALiAiI5vIrCefR9DoLPUc8LyDmix4OgvBGQxRgElfWE5SAgKlsWAiEsCQlZJntmkpnM3sv1/FE1PT2TWbp7uqdneq7v51Ofqr6r6q6701DX3EvdJaqKMcYY0xNPtgtgjDFmcLNAYYwxplcWKIwxxvTKAoUxxpheWaAwxhjTK1+2C5AJFRUVWlVVle1iGGPMkLFq1ap9qlrZ3b6cDBRVVVWsXLky28UwxpghQ0S29LTPmp6MMcb0ygKFMcaYXuVUoBCRhSKypK6uLttFMcaYnJFTgUJVn1bVxaWlpdkuijHG5IycChTGGGPSzwKFMcaYXuVUoLA+CmOMSb+ceo5CVZ8Gnp47d+41KWXw8V+geiX4g+ALgi8PfPnO2u+ufcGOpbvjRNL7pYwxJstyKlD023Pfh30f9S8Pb15cAIkLJLFAExdwEg1InY7tIU+v/ZTGmMywu0snCkcshPPvgXArhFs6L6Eun8OtEGo+9NhQ3P5wc+fjWup7OL+5f0UXb2LBx5/vfnbX/oIu657S4j77glZzMmYYsUDRlccHeUXOMpBUIdKWvoDU9fyWegjvcc9rhlCTe2xLauXtLpgcEnx6Czpd9gUKncXfvs63YGTMIGGBYrAQcf/qz4PgAD4HEo06QSYUvzR1WccHlq7HdDm+rREa9x2aR6Q1yYJJXPAogECR+7k9qBTFBZiijvT2QNN1sQBkTMosUAx3Hk/HzTSTopGOGkzXINLWBKFGJ8jEL6EmaGtwP7vbLfVwcJeb7gamZJrtxAOBYsiLX4o6tgNd00vc9KIu5xQ7Qd2YYSCnAoWILAQWTps2LdtFMV15vJlr0otGOmozhwSb+M9u0Gk92HlpqYf6HZ3T0L6v6w3EBZCSzgEnr9ipGQZLITiiY50/Ii691IKNGRJyKlD0e3isGZo83o6bczpEo07gaQ8abV0CS2sDtNbH7W/o2G7YAzUbnf3NByAa6v1avvyOoNE1iARHHLovfyQUlEN+mdPcZswAyKlAYUxaeDxxtZ+xqeej6jS1tdQ5QaOlLm454C51nfc37IF9Gzr2abTn/H35UFDmBI2Csrjt8i7bIzu284qtj8YkzQKFMZki0jHCq3hM8ueruv0y7YHkADTvh6YaaKqF5lpn3b69a627vZ8em848fqdWUljhLEWjoXAUFFW661FQWNmx9vr79U9gcoMFCmMGK5GOJrXSCYmfF404waWp1gkqsYDSZbthD1SvcNahpu7zyi/rEjxGQfFoKBkPJeOcpXic8/yOyVkWKIzJNR5vR1MUCQ7saG2Axj3QsNddu0v7duNe2L7aDSqNh55fUO4EjPbgER9ISsY7gc76VIYsCxTGmI4+mbKpfR/bUg8HdzojxWLLdjdtO2xf6dRYuiocBSMnw8gqGDHZ2R7hfi4Zb9PQDGL2yxhjkhMscZbKw3s+JtQSF0y2w4EtsH+Ls972Jrz3Z9BIx/HidWodI6ugYjpUzHDW5dOdIOLJqYmuh5xBHyhE5AjgeqAC+Iuq/jrLRTLG9MUfhLIpztKdSMgJIPs3dwSQ/Zuh9hN491FneHEsrwIon+YGDzeAjD4Kyg6zWsgAycq/sojcD5wL7FHVo+PSzwL+C/AC96rq7ar6AXCtiHiAhwALFMYMdV6/U3sYWXXoPlV3mPB6d9ngrKvfgvf+RGxEly8Io46A0UfDmFnOMvqogZ0Cp13tJ866p8A4xGUrHD8A3IVz4wdARLzA3cDngWpghYg8parvi8h5wP8DPJyFshpjBpKIM7KqeDRMWdB5X1sT1HwMu9c5w4F3r4UP/xvejrs1jJgM44+DCSfA+Lkw9pjMj8q681hn/aPcfGlaVgKFqr4qIlVdkucBH6vqJgAReQRYBLyvqk8BT4nIfwO/7y5PEVkMLAaYNGlShkpujMmqQAGMne0sXO6kqTr9Ibveg13vOkv1Slj3uLPf43NqG+PnwsR5MPlTUDo+a19hKBpMDXzjgW1xn6uB+SJyGnAhkAcs7+lkVV0iIjuBhYFA4PhMFtQYM4iIdAzFnXFGR/rBXU7A2L7SWa/5Paz4jbOvbCpULXCXU6CkH0/gDwODKVB0S1VfAV5J8Fib68kY4ygeA0ec6yzgPIi4+z3Y/Hf45G+w7glY/aCzr3waHHY6TD/TCRz2AGEngylQbAcmxn2e4KYlzGaPNcb0yON1+ivGHgMnfcsJHLvedQPHq7D6YXhriTOH1tRTYfoZMOPMvp+Kb+vmAcQseHvrft7eeoArTpxEns+b1rwHU6BYAUwXkSk4AeIy4EvJZGA1CmNMwjxeGDfHWU7+tvNulM1/h/XPw4bnYf1z8N84fRtHnQ9Hng8jJh6azzvLBrzo3fnr+r384sUNfOWkyWnPO1vDY5cBpwEVIlIN/FBV7xOR64DncYbH3q+q65LM12oUxpjU+PNh+uedRX/qDMn9aLnTRPXCD5xlwglOwDjqfKemEY3CG/c4548+uvf8M2x/Yxul+X783vQ/nJitUU+X95C+nF46rBPI12oUxpj+E3GePK88HE75Z6jd5ASMdY/DCzc7y4R5zuirmg3OOb1NCT8AahrbKCsMZCTvwdT0ZIwxg1PZVFjwv5ylZqMTMNY9ASvudSZDHHesE0yyaH9TGyMLMjMtfE4FCmt6MsZkXPlh8OnvOUvNRucp8//5V6dzPIt21bVw+Jg0veWxi5yaaUtVn1bVxaWlWXiE3xgz/JQfBiMmOZMaZrHpSVXZcaCFcaX5Gck/pwKFMcZkhXg6z4Y7wPY3hWgORRg3wgJFn0RkoYgsqavLzflWjDGDlCe7NYodB5oBLFAkwpqejDFZIR5nzqks2e4GigkjLVAYY8zgJJLVzuzt+61GkTBrejLGZIV42d/Uwm3LP8jK5XccaCbo92RseGxOBQprejLGZIV4CIcjLHl1E29s6uZ94Rm2o66ZcSPyEZGM5J9TgcIYY7JCPHhwOrOvf+TtWJ/BQNl+oIXxGWp2AgsUxhjTfx4vHpRTZ1TS1Bbha/e/xZ76lgG7/Pb9zRYojDFmUBMPHpSjx5dw71fnsuNAM4vu/gdvfVKb8Uu3hCLsa2jNWEc25FigsM5sY0xWuE1PXhHmTy3nsWtPwusRvvj/v861D6/ir+v30hLKzKioXXVOzSWTgSKn5nqy2WONMdmguIHC4/ztfdS4Ul7450/z61c28tDrW3hu3S78XmH8iHzKCgP4PB4K8rzcesGsfjcZtfeHZLLpKacCRX/dFGzjhabV+H5/IoIzeqB93bGS2MiCrsfEjzjomtbpWOn5mO7y6npsT3n3eL2uZevlmPi8+ipbt/82yZStS94igkc8sXJ4xIMHT0d63DE9psd/jjuma3rX68Tni9DpmPZzO63jy9h+nHtO++Lz+Jy1OGuvx4tXvB3rLtuHnOPx9HmcV7wZG+liEvfV2r8TGjeSU+LaaAoCPr57xuF887RpvLGphhWba9la20Rdc4iWUIRXPtrLKx/t4Yr5/XvRkAWKAbbWq7SiXDLtgk7pivPEpap22o7f111a12O7S+su797yO+S6sVWSZeuurHF59fi9u7luIsd0d11VPeRzlKiz1mhHmkZ7To//jPu563Fx6Yp2OiZKFJTYMUNRe8BoDyge8eD3+PF5fPg9/u63vT3vjz+mt3z8Xj9Bb5CAN0DQGyTPl9fjZ58nt281a0K1EAxwqufQ1vz8gJfPzBzFZ2aOiqVFo8pRP3yeTXv7/xrV7fubEYHRpXn9zqsnuf3rJUkUzvSO5F/m/Uu2i2KyJJEgEws0HJrevkQ04izRCFGNEtYw0aiTHo6GO9I0SiQa6Tg+/hz3uPj09u34/d0dF46GCWuYUCREKBoiHA0TinZst0XaaGxr7HZf1+O6/gGTCp/4yPPlkeftWIK+YGw735dPgb+AQn8hhb7C2HaBv4ACn5vuL6TAV9Cxz00fTDWqRF8u5/EIUyoK2bi3od/X3FbbxJiSYNrfkx3PAoUxcUQEr2Tuf7ihKBKN9BhEWiOttERaaIu00RJuSehzbAk7+1ojrdQ31dMYaqQp1ERTuInmcGLPIXjFS0mghJK8EkoDpRTnFVMaKKUkUEJpXimleR3bZcEyyvPLKQ+WE/QFM/Jv5e2mRtGTqZWFvFN9oN/X3FLbxKSygn7n05ucChT24iJj0s/rcZq0BlIkGqE53OwEj3ATTaGm2HZs3dZIfVu9s7TWU9dWR11LHVvrt1LXWsfBtoM91oaK/cVO0HADR0V+BRX5FYwuHM3YwrGMKxrHqIJR+D3JTYnh8yReu5laUcjytTtpC0cJ+FIfgLqlponPxjVrZUJOBQob9WRMbvB6vBQFiigKFKWcR1SjHGw7GAsitS217GveR01zjbNucdbr96/n9R2vczB0sNP5HvEwqmAU4wrHMaZwDOOLxjOxeCKTSyYzuWQyZcGyQ5q9kogTVFUUElXYWtvEtFGpfc+G1jD7GlqZVG41CmOMSZpHPLHmp4lM7PP4lnALuxp3saNxBzsbdrKjcYfzuWEH7+x9hxc2v0BYw7Hji/3FTCqZRGV+ZSwtmYpBVUUhAJv3NaYcKLbWNAEw2QKFMcZkXtAXpKq0iqrSqm73h6NhdjbsZHP9ZrbUb2FL/RY2129mV9Ou2DG+JGoUU8rdQFGT+sinrbXOuZPLClPOIxEWKIwxJgE+j4+JJROZWDKRBSzotO/Hj17E8w0f4pXER4iNLAxQmu/vV6DY7NYoMt30lFNTeBhjTDZ48RCR5GoU4DQ/fbIv9UDx8Z4GKovzKM3PzHso2lmgMMaYfvKIhzCS8HMU7aZVFrFhd+rPUny06yCHjy5O+fxEWaAwxph+6qhRJPdw4hFji9lzsJWahtakrxmJKhv2HOTwMRYoEJHzReQ3IvIHETkj2+UxxpiuPHiIiOCR5KaBmTmmBHBqBsnaWttESyiauzUKEblfRPaIyHtd0s8SkY9E5GMRuRFAVZ9Q1WuAa4FLs1FeY4zpjce9lXokuanEZ451bvIfpBAo2oNLLtcoHgDOik8QES9wN3A2cCRwuYgcGXfID9z9xhgzqHjEuZUqyQWKiqI8KoryWLc9+XfovLe9Dq9HmJGrNQpVfRXo+uqnecDHqrpJVduAR4BF4vh34FlVXd1TniKyWERWisjKvXv3Zq7wxhjTRXuNQgglfe6cSSN4e1vycz69vW0/R4wtJj+Q+elVBlMfxXhgW9znajft28DngItF5NqeTlbVJao6V1XnVlZW9nSYMcaknQfnZi1J1igAjp88kk/2NSbVoR2JKu9sq2POxJFJXy8VgylQdEtV71TV41X1WlW9p7dj7VWoxphsENqbnsJ9HHmo4yc7N/vVWxOvVazffZCG1jDHThyR9PVSMZgCxXboNCHLBDctYar6tKouLi0tTWvBjDGmNx1NT8nXKGaNLyXg8/D6xpqEz/nbBqd5/eRp5UlfLxWDKVCsAKaLyBQRCQCXAU8lk4HVKIwx2eDpR40i6Pdy0tRyXv5oT8LnvLp+HzNGFzG2NHOvP42XreGxy4DXgcNFpFpErlbVMHAd8DzwAfCoqq5LJl+rURhjsqE/fRQAp88cxSf7GhN6411dc4i3Pqnl1BkD1xebrVFPl6vqWFX1q+oEVb3PTV+uqjNU9TBVvTXZfK1GYYzJBsGZ5CmaQo0C4IyjRuMReOLtvlvbl6/dSVskyrmzx6V0rVQMpqanfrMahTEmG/pboxhbms+pMyp5dOU2QpGen+5WVR55ayvTRhUxe8LA3edyKlBYjcIYkx2pPXAX76snVbG7vpVlb23t8ZiXP9rDO9V1fP1TUw55u14m5VSgsBqFMSYbPOoGCk09UJx2eCUnTS3njuc/YlM3fRV7D7byg8ffY0pFIRcfPyHl66TCXlxkjDH9JG7TUyqjnmJ5iPAfF8/mvLv/ykX3PcbFx0/i+LEzEfFSvb+Z+//+CfubQjyy+HgCybxzNQ1yKlCIyEJg4bRp07JdFGPMMCJpaHoC2Na8hqJpt3OgrZZHdsCybUFCB06grXYBR46awF1fmsMxA/SQXTxrejLGmH6SWNNT6jWK9/a9x7df+jaVheX8dP6/8n8mX87JlXMIVrxG+cyfc8XntzFn0sBM2dFVTtUojDEmO5ymp2iKfRRRjfJvr/8bI4Mjue+wLzPyT9dBqJEvAtVTP81tI0r4yVs/IegLcuH0C9NY7sTkVI3CGGOyoT9zPQG8uOVFPqj9gOuPuoqRT90A5VPhK0/AZ25mwq73uXPVs5xUPJXb3ryNTQc2pbPoCcmpQGHDY40x2dG/PorH1j/G+KLxnFP9EbQdhIt/C4d9Bk793/Cdt/FNPpnbPniDgPj42aqfpbPgCcmpQGF9FMaYrFCn6SmSQtPTzoadvLnzTRZNXYhnzVKYeQ5UTO84IFgClzxIhSePr0cLeLX6Vd6veT9dJU9In4FCRE4SkbtF5F0R2SsiW0VkuYh8S0TsjmyMMaTemf3StpdQlHOC46BpHxx90aEHFZbDydfxxU9WE/QEeGz9Y/0tcFJ6DRQi8izwDZyJ+s4CxuK8pvQHQBB4UkTOy3QhjTFmMItG22+lPU+/0ZN/bP8Hk4onMWnrSvAGYNrnuj/w2CsoUeGM/PE8+8mztEXaUi9wkvqqUXxFVa9W1adUdYeqhlW1QVVXq+rPVPU04LUBKGdCrI/CGJMN6k4KKGhS57VGWlmxawUnjzsZtr0J44+HvB7egV06HiadxBm1u2kMNbJy18r+FjthvQYKVd3XVwaJHDNQrI/CGJMNUdrnXUouULyz5x1aIi18asw82PkOTJjb+wmHnc78nR8R9ObxSvUrKZU1FX01PR0Ukfq45WD8eqAKaYwxg1k09sBdck1P7+x9B4A5UT9E2mDCvN5POOx0gqrML5rMazsGrjGn1wfuVLWHOpAxxph2UbcikWygeHffu1SVVFG650Mnoa8axZhZ4M3juIiXvzZsoaa5hvL8zL8ONeHhsSJyjIhc5y6zM1koY4wZSmI1iiSanlSVtXvXMqtiFuz9EIIjoHhs7yd5/TBmFnPqnfdrr9m7JuUyJyOhQCEi1wNLgVHuslREvp3JghljzFDRXqOIRhN/jmJX4y5qWmqYVTkL9q2HysMhkXdMjD+OI3d8gN/jZ82eQRQogKuB+ar6r6r6r8CJwDWZK5Yxxgwd0ahzg0+mRvFJ/ScATB8x3QkUFTMSO3HMbPJCjcwonsyHtR8mXdZUJBooBDo9mx5x0wYVGx5rjMmGaPsDd9HE+yh2Ne4CYKw3Hxr3OjWKRLgBZUagjPX71ydX0BQlGih+C7wpIj8SkR8DbwD3Za5YqbHhscaYbIiquOvEm56aQk0AFNU7AYOKRAOFM73HDPVR21LLvubMP6GQ0DTjqvpzEXkFOMVNukpV385YqYwxZgiJqIA404UnKhQNAeA/UO0klB+W2IkFZVBQwfQWJ9Cs37+eivyKpMqbrEQ7sw8D1qnqncBaYIGIDPxrlowxZhBqr1Ek00cRCxQH3RpFyfjEL1gxnRl1uwHYsH9D4uelKNGmpz8BERGZBtwDTAR+n7FSGWPMEBJxh8eiyQcKX902KBwF/mDiFyw/jJG1WxiRN4It9VuSKWpKEg0UUXWmRbwQuEtV/1+cCQKNMWbYi8biQxJNT5EQfo8fqd8OpROSu2DpJGjYxaSiCWw9uDW5c1OQaKAIicjlwFeBZ9w0f2aKZIwxQ0s49mR2cjUKv8cPB7bBiInJXdANLBOD5Wyr35bcuSlINFBcBZwE3Kqqn4jIFODhzBWrg4hMFZH7ROSPA3E9Y4xJVG1jG/Utodg04ykFirpqKE0tUEzyFrKzcWfGpxxPdNTT+8B34j5/IiJ/T/WiInI/cC6wR1WPjks/C/gvnDeV36uqt6vqJuBqCxTGmGw60NTGe9vrWbu9jrXbD/BudR3V+5uZWJbPSfkRCCY315MTKLwQbk45UEzEi6JUN1QztXRqcnkkoddAISJe4IvAeOA5VX1PRM4FbgLygTkpXvcB4C7goS7Xuhv4PFANrBCRp9wgZYwxA6amoZW12+tYt6Oe97bX8d6OOrbVNsf2Tyor4JiJIzh39jju+etGSr2NUJLkqKdICH/7c8tJ91G4gaLVqUlsP7g9e4EC56G6icBbwJ0isgOYC9yoqk+kelFVfVVEqrokzwM+dmsQiMgjwCLAAoUxJiMaW8Ns3NvA+t0NbNhzkA27G/hgZz0761pix0wuL2D2+BF8ad5kjhpXwuwJpYwoCMT2H2hq493V7lQaydYo2uNKaRJDYwF8eVA0mtHNBwDY07QnufOT1FegmAvMVtWoiASBXcBhqlqTgbKMB+J7ZaqB+SJSDtwKzBGR76vqT7o7WUQWA4sBJk2alIHiGWOGIlVlz8FWNu9rZEttExv3NLB+90E27Gmgen9HLSHg9TC1spB5U8qYNb6Uo8aVcuS4Ekrzex+3c+PZM7l89YvOtZJ8jsLffnzRmOS/WOkEKur3Igi7m3Ynf34S+goUbeo2uqlqi4hsylCQ6JF7vWsTOG6JiOwEFgYCgeMzXzJjzGDRFo6ys66ZrbVNbK5pYmtNo7tuYkttIy2hjr/0Az4Ph1UWcdykkVw6dyLTRxczfXQRk8sK8HkTfvNCzIiCAL+47DguXv3b5Duz28fVFqbwZHXRGPz7N1NeUZ71GsVMEXnX3RbgMPezAKqq6XwvxXacZq52E9y0hKnq08DTc+fOtZltjckR0aiyr7GVnQda2HGgme0HmtlZ52zvcNf7Glo7PeuW5/MwubyAyeWFLJheweSKQiaXFTC5vIDxI/JTCgi9mTiyCEihM1sjUFDuvGciWUWjoPotRk06nl1Nu5I/Pwl9BYojMnr1zlYA092ht9uBy4AvJZOBiCwEFk6bNi0DxTPGpFNbOMq+hlb2Hmxlz0Fn7Wy3ONsNreypd9LaIp1vwPl+L+NGBBk3Ip+Zh49irLs9qayAqvJCRhXn4fEM3ATX4kn+VajhSBhfJOw8lZ2KotHQuI/R+aOobkzqb+qk9RUotmofdSkRkb6O6eacZcBpQIWIVAM/VNX7ROQ64Hmc4bH3q+q6ZPK1GoUx2dMSirC/qY3axjb2N4aobWpjf6P7ualj3R4YDjSFus2nrDBAZVEelcV5zJ9SSGVJHuNH5DO2NN8JDqX5jCjwI4m85GeAOIM2k+ujiBJFIiEoSnLEU7uiSkAZHShh1Z7VqeWRoL4Cxcsi8ifgSVWNPScuIgGcmWS/BryMM9w1Yap6eQ/py4HlyeQVz2oUxvSPqtLYFqG+OURdc4j65hD1LeGOzy0hDjSFOt349zeGqG1soznU8xTbIwr8lBUEGFkYYGpFEfOnlFNZnMeoYicgtC8VRXn409wsNBBEkp/rSVXxREJQMjq1ixY55432BKlvq6c53Ey+Lz+1vPrQV6A4C/g6sMxtEjoABHH+4n8B+MVgmm7cahRmuGsLR2lsDdPgLo2tYQ62hKlv6fnGX98cjtsOxc1b1L3iPB9lRQFGFgQYVRzk8NEllBX6GVEQoKzQSS8rDFBW6GdkQYDSfH/a+wQGGyH5d2ZHNYo30ub0NaTCDRSj3Nv4nqY9TC6ZnFpefeg1UKhqC/Ar4Fci4gcqgGZVPZCR0hgzDLWGIzS0hGlsjXS6wXe92Te2hmls69h29kc6HdsW7ruNPN/vpSTfR0nQT0m+n4qiAIdVFlKS73fTfJTGtp11ab6TXpTny/mbfio6+iiSaIWPRvBotB+BwjmvKOQ889EYakwtnwQkNIUHgKqGgJ0ZK0kaWNOTyRRVpSUUpaktTFNbhOaQc4NubovQ1BahKRShqbXzvqa2iLM/bl9TKEJzmxMUmtqcm3soktjNpSDgpTDPR3Gej8I8H4V5XsaPyKcoz0tR0NdlX8e2c7N31sVBH3k+b4b/tYYhSb4zOxppdTZS7cx2zwu2NgDQ2p5fBiQcKIYCa3oanlSVtkiUlrYoLWHn5hxbh6K0hCK0hJwbeEso6q6dpcm90Te7AcBZ4rfdfaFIMs3PeD1CQcDrLj7y/c52ab6fsSXB2E2/MM9HcdBHYcBLUdBPUZ6TXtS+uAGgMODDO4CjeExy2vsokml60mgYDwr5Kb4DLlAA/gLy25yaRHO4uY8TUpdTgcIMHpGo0haO0hqO0BqOdr5RuzfyltgNvWN/aw839PjPznGdPyc37s4hAgV+L/kB56/z9pt5QcBHRVEeBQFnX0HAS2Hcdr4bAAoDvth2+3nt+wNez6AalWMyK5XhsaoRRAFvoM9je5RfRsANFK3hLNcoRKQQp28iKiIzgJnAs25z1KBhTU+O9r+wW8NRWkNRZzvk3LCdtI7t+Jt5z+lx273sa2vPPxxJuDmlKxEI+pybbdDnIRjwdnz2exhZ4CfP79zUg36Pu3aW9u38gIegz3vIufHHBv0eu5mbtGkfHksyNQpVZ0rAVB62a1cwEl9LPQAR7XnUWX8lWqN4Fec92SNxRjutAC4FrshUwVIxWJqewu03afemGbuBdrrhRuJu4p1vxK2hCK2x9ENvyh035I58nBt1Rx795REI+r3k+TwEfB7yfM52nr9juyjP1+O+runxN+xON/W4tXO+3bzN0BN7jiKJqm1UI26g6EeNoqAcf3M9+CGs4dTz6UOigUJUtUlErgZ+par/ISJrMlaqLFGUgy0h/v25Dzv+au5ys27r5gbddV+kr/GFCcjzedybbDc3Xp+Hknx/XHr3N+u8+H3+ju1An+keG9liTBIkhc5sVXUG1fanRpFfhvfgNvBDJJr9GoWIyEk4NYir3bRBN3Siv01PkahysCXMvX/b1Plm2+lm7TRhjMj3x27KAa+n0w23/YbtpHu7vZHH35Q7jnGuYU0ixgwxqXRma9SpUXj60/RUhrf5ABQVDoqmpxuA7wOPq+o6EaxUPD0AAB0cSURBVJmK80T2oJKOpqc8n4cNt34hjaUyxuQ88SCqSTU9qUYR1X53Zvua64DC7NcoVPWvwF/jPm8i7tWoxhgzrIm0T6md8CnR9hpFvzqzy/DiNHeFopkbW5ToqKeX6aY7X1VPT3uJjDFmqBGP+1LTJGoUpGPUUzle95KDoenpe3HbQeAiIHNd7MYYM5S4gSLZPgqnM7t/TU9e95qDoelpVZekf4jIWxkoT7/YcxTGmKxKcvbYfndm54/ANwA1ioTGQIpIWdxSISJnAqUZK1WKVPVpVV1cWjroimaMyWXuKMWkOrNp78zuR6DIK4kNPw1Hs/8cxSqcxjfBaXL6hI5hssYYM8y5ndlJnNHRmd2PpqdgCV43OGW9j0JVp2SsBMYYkwNEk+ujIB0P3OWVxG7iWatRiMjpqvqSiFzY3X5V/XNmimWMMUOICIKSTJ0i2t5I4+nHs8v+fMTjw4tktenpVOAlYGE3+xSwQGGMMST/HIWq4pF+TpUj4vZTSPbmelLVH7rrqzJWgjSyUU/GmKyQnvsomkJNvLnzTT7c/yG7GncR1Sh+j5/90db0TNUTdDq0o9H+Twbak76anv5Xb/tV9efpLU7/DJbZY40xw1N8H8XBtoMseXcJf/joD7GXClXkV+AVL6FoiGaNUJWO/ue8EnzUZ3X22GJ3fThwAvCU+3khMOieozDGmOzo3PS0rmYdN7x8A3ua9nD2lLO5cNqFzK6cTdAXjJ0RffLbePY+3/9LB0vxaX32+ihU9ccAIvIqcJyqHnQ//wj474yVyhhjhpJY05OyZs8arn3xWkoDpTx89sPMrpzd7SmeaLh/Q2Pb5ZXgbdVB8RzFaKAt7nObm2aMMcatUewMN/Cdl75DRX4F955xL2MKx/R8SqStf0Nj2+UV42vR7D9HATwEvCUij7ufzwcezEyRjDFmiBFncOxfGrdQ5C/il6f/svcgARAN9W/6jnbBErz7o9mvUajqrSLyHHCKm3SVqr6dsVIZY8wQ0+hxhrreNP8mppQm8IxyJJS2pie/RohkO1CAMzGgiGzDmT0WEZmkqlszVjJjjBkyhOJIFPEHOXfquYmdkq6mJ3caj3C4tf959SDRSQHPE5ENOHM8/dVdP5uxUnW+dqGIPCgivxGRKwbimsYYkxQRnty+g5cmfjHxZyMioTT1UZTgUwhHWvqfVw8SfSzw/wNOBNa78z59Dngj1YuKyP0iskdE3uuSfpaIfCQiH4vIjW7yhcAfVfUa4LxUr2mMMZkjVEai5CUzHUe6mp6CJXhRwuHsB4qQqtYAHhHxqOrLwNx+XPcB4Kz4BBHxAncDZwNHApeLyJHABGCbe1jmuvU7SpL5SxhjclQSkwJGQ+BJuPW/Z/4CfMqg6KM4ICJFwKvAUhHZAzSmelFVfVVEqrokzwM+dt/HjYg8AiwCqnGCxRp6CWwishhYDDBp0qRUi2aMMclrb25KZp7xSFt6ahRePz6UcCRz78xOtEaxCGgC/hl4DthI9xMF9sd4OmoO4ASI8TgTD14kIr8Gnu7pZFVdoqpzVXVuZWVlmotmjDG9aW+JSCJSpKuPwhtw+iiyOIUHAKraXnuIAg+KiAe4HFiaqYJ1uXZCkxLapIDGmKxIZXK/dAUKjx8vSksGm556rVGISImIfF9E7hKRM8RxHbAJ+GKay7IdmBj3eYKbZowxQ0MS04ynr+kp4PZRZO+d2Q/jTAi4FvgG8DJwCXC+qi5Kc1lWANNFZIqIBIDL6JiEMCH2zmxjTHZks+nJ7zxHkcWmp6mqOgtARO4FdgKTVLVf47BEZBlwGlAhItXAD1X1Pre28jzgBe5X1XX9uY4xxgyIWGd2sqOe0tRHQWZrFH0Filg3uqpGRKS6v0HCzevyHtKXA8tTzdf6KIwxWSGp1CjS1PTk8eHLcI2ir6anY0Sk3l0OArPbt0WkPmOlSpE1PRljhoxIOE2d2R68QFiz9IY7Ve3HW78HntUojDFZlXRndjoCRXuNInud2UOK1SiMMdnT01uzu6Hq9FGko+lJvE4fRQZrFDkVKIwxJmtEEq9RtD/zkI7ObI/XrVFYoEiIiCwUkSV1dXXZLooxZthJ4qG7iPvC0HQ0PYkXr2a2jyKnAoU1PRljsivBGkUsUKRj1JPXmevJ+iiMMWaQS6bpKeI2PaWlRuGhJBqlWcOEMjQxYE4FCmt6MsZkTxKd2elsevJ4KYs4tYnaltr+59fdJTKSa5ZY05MxJmuSmRgwrU1PPsoiTv+EBQpjjBnssjHqSbyUZ7hGkYbXKxljjOmx6UkVPn4R1vwetq+E5gMdb7bzpuEW7PFmvEZhgcIYY9Khu87sljp4/Fr4aDkUjoKqU6BoFISaINwGk05Ow3U9Ge+jyKlAkZYpPOyV2caYlHSpUbQ2wIPnwe734IxbYd5i8KWhT+KQywqFCAE81DTXpD9/cqyPwjqzjTFZE9+ZrQpPfgt2rYVLl8LJ12UmSMQu7aXMm0dNiwUKY4wZ3NqbntY9Du8/AaffDIeflfnreryUefJs1JMxxgxubo0i3Aov/B8YeyycfP0AXdrLSPFT15qZZ8gsUBhjTDq0d2avegDqq+FzP0rPqKZEeHzki4fmcHNmss9IrsYYM+wIaAReuwsmnQRTTxu4S3s8BLFAkRCbwsMYk1Ufvwh1W2H+PyX3pHZ/iZcgHlrC/X5TdbdyKlDYqCdjTNaIQM3HUDQaZp47sNf2eMlHaIlYoDDGmMGrtd5Zz7okPZP9JUO8BBGaw81oMq9jTZAFCmOMSacjzx/4a3q85OMhqlFC0fRPNW6Bwhhj0mn88QN/TfFwib+S5Rcux+dJ/0irnJrCwwx+oVCI6upqWloy05ZqBrdgMMiECRPw+we4aWagzDgLPFn4+9vjY4R6GFE8MSPZW6AwA6q6upri4mKqqqqQgRwVYrJOVampqaG6upopU6Zkuzjp96/7B3akUzyPt2Pq8kxkn7GcjelGS0sL5eXlFiSGIRGhvLw8d2uTHk/2AoV4nWc4MmTQBwoRmSoi94nIH7NdFpMeFiSGL/vtM8TjhWg0c9lnLGdARO4XkT0i8l6X9LNE5CMR+VhEbuwtD1XdpKpXZ7KcxhgzpIlnSNcoHgA6TZ0oIl7gbuBs4EjgchE5UkRmicgzXZZRGS6fMX165ZVXOPdc5wGqp556ittvvz3LJTKmC48vo30UGe3MVtVXRaSqS/I84GNV3QQgIo8Ai1T1J0DKjzOKyGJgMcCkSZNSzcYMI6qKquJJYpTKeeedx3nnnZfBUhmTgqEcKHowHtgW97kamN/TwSJSDtwKzBGR77sB5RCqugRYAjB37tz0P5po0u7HT6/j/R31ac3zyHEl/HDhUT3u37x5M2eeeSbz589n1apVzJs3j7Vr19Lc3MzFF1/Mj3/8YwCee+45brjhBgoKCjjllFNi5z/wwAOsXLmSu+66iyuvvJJzzz2Xiy++GICioiIaGhrYuXMnl156KfX19YTDYX7961+zYMGCtH5PYzrx+iGSW4EiKapaA1ybyLHpeBWqdbXlvg0bNvDggw9y4oknUltbS1lZGZFIhM9+9rO8++67zJgxg2uuuYaXXnqJadOmcemllyaV/+9//3vOPPNMbr75ZiKRCE1NTRn6Jsa4PD7I0ISAkJ1AsR2IfypkgpvWb6r6NPD03Llzr0lHfiazevvLP5MmT57MiSeeCMCjjz7KkiVLCIfD7Ny5k/fff59oNMqUKVOYPn06AF/+8pdZsmRJwvmfcMIJfP3rXycUCnH++edz7LHHZuR7GBPj9XfMNZUB2RgeuwKYLiJTRCQAXAY8lY6MbZpxk4jCwkIAPvnkE+644w7+8pe/8O6773LOOeckNcbf5/MRdYckRqNR2traAPj0pz/Nq6++yvjx47nyyit56KGH0v8ljInnyWzTU6aHxy4DXgcOF5FqEblaVcPAdcDzwAfAo6q6Lh3Xs2nGTTLq6+spLCyktLSU3bt38+yzzwIwc+ZMNm/ezMaNGwFYtmxZt+dXVVWxatUqwBkNFQo5k7Ft2bKF0aNHc8011/CNb3yD1atXD8C3McOa1wcZmAywXaZHPV3eQ/pyYHm6r5eOPgozfBxzzDHMmTOHmTNnMnHiRD71qU8BznxES5Ys4ZxzzqGgoIAFCxZw8ODBQ86/5pprWLRoEccccwxnnXVWrKbyyiuv8NOf/hS/309RUZHVKEzmefwQyVygkEzMXZ5tc+fO1ZUrVyZ93hfuPZoxkSLu/6c3MlAqA/DBBx9wxBFHZLsYJovsv4EM+PM/wdbX4Ia1KWchIqtUdW53+wb9FB7JsD4KY8yw5PUN3T6KgWZ9FMaYYcnjz2gfRU4FCmOMGZa8me2jyKlAYU1PxphhyZ8Pocw92JlTgcKanowxw1KgGCJtEG7LSPY5FSiMMWZYyity1m0NGcneAoUZdqqrq1m0aBHTp0/nsMMO4/rrr489VW3MkBRwA0Xroc/7pENOBQrrozB9UVUuvPBCzj//fDZs2MD69etpaGjg5ptvznbRjEldhmsUg3722GTYpIBDzLM3wq7UHxDq1phZcHbPLxZ66aWXCAaDXHXVVQB4vV7+8z//kylTpjBlyhReeOEFGhsb2bBhA9/73vdoa2vj4YcfJi8vj+XLl1NWVsbGjRv51re+xd69eykoKOA3v/kNM2fOZOPGjVxxxRU0NjayaNEifvGLX9DQ0EBDQwOLFi1i//79hEIhbrnlFhYtWpTe722Gt0Cxs261pidj+m3dunUcf/zxndJKSkqYNGkS4XCY9957jz//+c+sWLGCm2++mYKCAt5++21OOumk2FQcixcv5pe//CWrVq3ijjvu4Jvf/CYA119/Pddffz1r165lwoQJsfyDwSCPP/44q1ev5uWXX+a73/0uuTgjgsmiWI0iM01POVWjMENML3/5Z8tnPvMZiouLKS4uprS0lIULFwIwa9Ys3n33XRoaGnjttde45JJLYue0trYC8Prrr/PEE08A8KUvfYnvfe97gNPcddNNN/Hqq6/i8XjYvn07u3fvZsyYMQP87UzOivVRWNNTn2xSQNOXI488kj/+8Y+d0urr69m6dSs+n4+8vLxYusfjiX32eDyEw2Gi0SgjRoxgzZo1CV9z6dKl7N27l1WrVuH3+6mqqkpqOnNj+mSjnhJnz1GYvnz2s5+lqakp1owUiUT47ne/y5VXXklBQUGf55eUlDBlyhQee+wxwKktvPPOOwCceOKJ/OlPfwLgkUceiZ1TV1fHqFGj8Pv9vPzyy2zZsiXdX8sMd9ZHYUz6iAiPP/44jz32GNOnT2fGjBkEg0Fuu+22hPNYunQp9913H8cccwxHHXUUTz75JAC/+MUv+PnPf87s2bP5+OOPaf+D5YorrmDlypXMmjWLhx56iJkzZ2bku5lhLMN9FDbNeBybZjzzcnmK6aamJvLz8xERHnnkEZYtWxYLIqZDLv83kFX/VgEnfQs+/+OUTu9tmvGc6qMwJptWrVrFddddh6oyYsQI7r///mwXyQwneUX2HIUxg92CBQti/RXGDLhAsfVRGGOM6UUGaxQ5FShsCg9jzLAVKLK5nhLR3+Gxudetb4wZNvxBCLdmJOucChTGGDNsefMgnJkHOS1QmGHH6/Vy7LHHxpbbb0//VCLJPJcR7wtf+AIHDhzo9Ziqqir27duXUv7pyOeee+6JPbD4wAMPsGPHjn6XxaSBL895eVEmss5IrsYMYvn5+UlNwZGK2267jZtuuinh41UVVWX58uUZLFV6XHvttbHtBx54gKOPPppx48ZlsUQGcAJFhpqeLFCYrPn3t/6dD2s/TGueM8tm8i/z/iWlc5977jluuOEGCgoKOOWUU9i0aRPPPPMMP/rRjygqKopN8nf00UfzzDPPUFVVxfnnn8+2bdtoaWnh+uuvZ/Hixdx44400Nzdz7LHHctRRR7F06VJ+/vOfx56r+MY3vsENN9zA5s2bOfPMM5k/fz6rVq1i+fLlnHrqqaxcuZKKiopu8+7JPffcw8aNG/npT38KODfwlStXctddd/G73/2OO++8k7a2NubPn8+vfvUrvF5vp/O7Kx/AQw89xB133IGIMHv2bB5++OHYv0dVVRUrV67kiiuuID8/n1tvvZXf/OY3sYkR/+d//odf/epXPP744yn9HiZJXqtRGJM27Tfxdt///vdZtGgR11xzDS+99BLTpk3j0ksvTSiv+++/n7KyMpqbmznhhBO46KKLuP3227nrrrtitZZVq1bx29/+ljfffBNVZf78+Zx66qmMHDmSDRs28OCDD3LiiScmlHd5eXm35bjooos46aSTYoHiD3/4AzfffDMffPABf/jDH/jHP/6B3+/nm9/8JkuXLuWrX/1q7NyeyhcIBLjlllt47bXXqKiooLa2ttM1L774Yu666y7uuOMO5s6di6ry3e9+l71791JZWclvf/tbvv71ryf072jSwBcYvjUKETkfOAcoAe5T1ReyXCSTJqn+5d9f3TU9rVmzhilTpjB9+nQAvvzlL7NkyZI+87rzzjtjfzFv27aNDRs2HHIz//vf/84FF1xAYWEhABdeeCF/+9vfOO+885g8eXK3QSLRvNtVVlYydepU3njjDaZPn86HH37Ipz71Ke6++25WrVrFCSecADhBctSoUQmVT0S45JJLqKioAKCsrKzXfwsR4Stf+Qq/+93vuOqqq3j99ddjfRlmAHjzIDIEA4WI3A+cC+xR1aPj0s8C/gvwAveqao+9iar6BPCEiIwE7gAsUJgB5fP5iEajsc/tU4S/8sorvPjii7z++usUFBRw2mmnJT19ePvNuatU8r7ssst49NFHmTlzJhdccAEigqryta99jZ/85CdJlStVV111FQsXLiQYDHLJJZfg8w36v0VzRwZrFJke9fQAcFZ8goh4gbuBs4EjgctF5EgRmSUiz3RZ4v/0+YF7njFpN3PmTDZv3szGjRsBWLZsWWxfVVUVq1evBmD16tV88skngDN9+MiRIykoKODDDz/kjTc6JpP0+/2EQiHAmdrjiSeeoKmpicbGRh5//HEWLFjQa3l6y7snF1xwAU8++STLli3jsssuA5xp1f/4xz+yZ88eAGpraw+Z5ryn8p1++uk89thj1NTUxM7tqri4mIMHOx7yGjduHOPGjeOWW26JvW7WDBDvEO3MVtVXRaSqS/I84GNV3QQgIo8Ai1T1Jzi1j05ERIDbgWdVdXUmy2uGh659FGeddRa33347S5Ys4ZxzzqGgoIAFCxbEboAXXXQRDz30EEcddRTz589nxowZsfPuuecejjjiCA4//PBOTUiLFy9m9uzZHHfccSxdupQrr7ySefPmAU5n8Zw5c9i8eXOPZewt756MHDmSI444gvfffz92rSOPPJJbbrmFM844g2g0it/v5+6772by5Mmx84477rhuywdw8803c+qpp+L1epkzZw4PPPBAp2teeeWVXHvtteTn5/P666+Tn5/PFVdcwd69e22G2IHmywMUohHwePs8PBkZn2bcDRTPtDc9icjFwFmq+g3381eA+ap6XQ/nfwf4GrACWKOq9/Rw3GJgMcCkSZOOT+XlMGffezRjbZrxjBoqU0y/8sor3HHHHTzzzDPZLsqQc9111zFnzhyuvvrqbvcPlf8Ghpz3n4J1f4YLljjNUEka0tOMq+qdwJ0JHLdERHYCCwOBwPGZL5kxpqvjjz+ewsJCfvazn2W7KMPPkec5SwZkI1BsBybGfZ7gpvWbqj4NPD137txr0pGfGb5OO+00TjvttGwXY8hZtWpVtotgMiAbU3isAKaLyBQRCQCXAU+lI2ObPXZoyMW3KprE2G8/NGU0UIjIMuB14HARqRaRq1U1DFwHPA98ADyqquvScb3+zh5rMi8YDFJTU2M3jGFIVampqSEYDGa7KCZJmR71dHkP6cuBtE9qIyILgYXTpk1Ld9YmTSZMmEB1dTV79+7NdlFMFgSDQSZMmJDtYpgkDfrO7GRYH8Xg5/f7mTJlSraLYYxJgk0zbowxplc5FSisM9sYY9IvpwKFdWYbY0z6ZfzJ7GwQkb1A8o9mOyqA/r8+bHDK5e8Guf39cvm7QW5/v6Hy3SaramV3O3IyUPSHiKzs6TH2oS6Xvxvk9vfL5e8Guf39cuG75VTTkzHGmPSzQGGMMaZXFigO1fdrzYauXP5ukNvfL5e/G+T29xvy3836KIwxxvTKahTGGGN6ZYHCGGNMryxQuETkLBH5SEQ+FpEbs12edBORzSKyVkTWiMjKbJenv0TkfhHZIyLvxaWVicj/iMgGdz0ym2VMVQ/f7Ucist39/daIyBeyWcZUichEEXlZRN4XkXUicr2bPuR/u16+25D/7ayPAhARL7Ae+DxQjfPOjMtV9f2sFiyNRGQzMFdVh8KDP30SkU8DDcBDca/Z/Q+gVlVvd4P9SFX9l2yWMxU9fLcfAQ2qekc2y9ZfIjIWGKuqq0WkGFgFnA9cyRD/7Xr5bl9kiP92VqNwzAM+VtVNqtoGPAIsynKZTC9U9VWgtkvyIuBBd/tBnP9Jh5wevltOUNWdqrra3T6I806a8eTAb9fLdxvyLFA4xgPb4j5XkyM/cBwFXhCRVSKyONuFyZDRqrrT3d4FjM5mYTLgOhF5122aGnJNM12JSBUwB3iTHPvtunw3GOK/nQWK4eMUVT0OOBv4ltu8kbPUaVPNpXbVXwOHAccCO4GfZbc4/SMiRcCfgBtUtT5+31D/7br5bkP+t7NA4dgOTIz7PMFNyxmqut1d7wEex2luyzW73Xbi9vbiPVkuT9qo6m5VjahqFPgNQ/j3ExE/zo10qar+2U3Oid+uu++WC7+dBQrHCmC6iEwRkQBwGfBUlsuUNiJS6HauISKFwBnAe72fNSQ9BXzN3f4a8GQWy5JW7TdR1wUM0d9PRAS4D/hAVX8et2vI/3Y9fbdc+O1s1JPLHbL2C8AL3K+qt2a5SGkjIlNxahHgvP7290P9+4nIMuA0nCmcdwM/BJ4AHgUm4Uwz/0VVHXKdwj18t9Nwmi4U2Az8U1yb/pAhIqcAfwPWAlE3+Sactvwh/dv18t0uZ4j/dhYojDHG9MqanowxxvTKAoUxxpheWaAwxhjTKwsUxhhjemWBwhhjTK8sUJhhT0TK42b23NVlps/XMnTNOSJyn7t9pYioiHwubv/5btrFSeb7iIhMT3d5zfBmgcIMe6pao6rHquqxwD3Af7Z/VtWTM3TZm4A74z6vxXnQs93lwDvJZOjOgvxr4H/3u3TGxLFAYUwvRKTBXZ8mIn8VkSdFZJOI3C4iV4jIW+57Pg5zj6sUkT+JyAp3+VQ3eRYDs1U1PhD8DZgnIn53rqBpwBr3+NNF5Im48z8vIo+3l09EfiYi7wAnufl8TkR8mfkXMcORBQpjEncMcC1wBPAVYIaqzgPuBb7tHvNfODWSE4CL3H1dzeXQaRwUeBE4E2fK7fgpZF4GZopIpfv5KuB+d7sQeFNVj1HVv7vzCX3sltWYtLBAYUziVrjvHGgFNgIvuOlrgSp3+3PAXSKyBudmX+LWEOKNBfZ2k/8jOM1PlwHL2hPd2VQfBr4sIiNwag7PursjOJPQxdsDjEv62xnTA6ueGpO41rjtaNznKB3/L3mAE1W1pZd8moFg10RVfUtEZgFNqrremWMu5rfA00AL8Jiqht30FlWNdMkq6F7DmLSwGoUx6fUCHc1QiMix3RzzAU4fRHduxOno7kRVdwA7gB/gBI3ezGAIzlBqBi8LFMak13eAue7bzN7H6dPoRFU/BErbp37vsu9ZVX25h7yXAttU9YOeLi4io4FmVd2VWvGNOZTNHmtMFojIPwMHVbW7zu6ezrkLeFtV7+sj3/rejjEmWVajMCY7fk3nPo9eicgqYDbwuz4OPQA82I9yGXMIq1EYY4zpldUojDHG9MoChTHGmF5ZoDDGGNMrCxTGGGN6ZYHCGGNMr/4vn2BzuzJGpsUAAAAASUVORK5CYII=\n", "text/plain": [ "<Figure size 432x288 with 1 Axes>" ] @@ -104,9 +104,17 @@ "plt.xlabel('Time (Myr)')\n", "plt.ylabel('Radius (Rsol)')\n", "plt.yscale('log')\n", + "plt.legend()\n", "plt.savefig(os.path.join(result_dir, 's'))\n", "plt.show()" ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] } ], "metadata": {