{ "cells": [ { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0 Surface(lbl='Obj', profile=Spherical(c=0.0), interact_mode=transmit)\n", "1 Surface(profile=Spherical(c=0.023413720440177945), interact_mode=transmit)\n", "2 Surface(profile=Spherical(c=0.005118231139318251), interact_mode=transmit)\n", "3 Surface(profile=Spherical(c=0.04878048780487805), interact_mode=transmit)\n", "4 Surface(profile=Spherical(c=0.0), interact_mode=transmit)\n", "5 Surface(profile=Spherical(c=0.06693440428380187), interact_mode=transmit)\n", "6 Surface(profile=Spherical(c=0.0), interact_mode=transmit)\n", "7 Surface(profile=Spherical(c=-0.06693440428380187), interact_mode=transmit)\n", "8 Surface(profile=Spherical(c=0.0), interact_mode=transmit)\n", "9 Surface(profile=Spherical(c=-0.04878048780487805), interact_mode=transmit)\n", "10 Surface(profile=Spherical(c=0.0), interact_mode=transmit)\n", "11 Surface(profile=Spherical(c=-0.023413720440177945), interact_mode=transmit)\n", "12 Surface(lbl='Img', profile=Spherical(c=0.0), interact_mode=transmit)\n", "0 Gap(t=10000000000.0, medium=Air())\n", "1 Gap(t=3.99, medium=Glass(nd=1.7343, vd=28.19, mat='', cat=''))\n", "2 Gap(t=0.2, medium=Air())\n", "3 Gap(t=7.18, medium=Glass(nd=1.67133, vd=41.64, mat='', cat=''))\n", "4 Gap(t=1.29, medium=Glass(nd=1.7919, vd=25.55, mat='', cat=''))\n", "5 Gap(t=5.35, medium=Air())\n", "6 Gap(t=7.61, medium=Air())\n", "7 Gap(t=1, medium=Glass(nd=1.65222, vd=33.6, mat='', cat=''))\n", "8 Gap(t=5.22, medium=Glass(nd=1.79227, vd=47.15, mat='', cat=''))\n", "9 Gap(t=0.2, medium=Air())\n", "10 Gap(t=3.69, medium=Glass(nd=1.79227, vd=47.15, mat='', cat=''))\n", "11 Gap(t=37.32, medium=Air())\n", " r t medium mode zdr sd\n", " Obj: 0.000000 1.00000e+10 air 1 4.1421e+09\n", " 1: 42.710000 3.99000 734.282 1 20.051\n", " 2: 195.380000 0.200000 air 1 20.041\n", " 3: 20.500000 7.18000 671.416 1 15.123\n", " 4: 0.000000 1.29000 792.256 1 14.636\n", " 5: 14.940000 5.35000 air 1 10.221\n", "Stop: 0.000000 7.61000 air 1 9.2442\n", " 7: -14.940000 1.00000 652.336 1 10.533\n", " 8: 0.000000 5.22000 792.471 1 14.644\n", " 9: -20.500000 0.200000 air 1 14.265\n", " 10: 0.000000 3.69000 792.471 1 17.200\n", " 11: -42.710000 37.3200 air 1 17.217\n", " Img: 0.000000 0.00000 1 21.882\n", "efl 52.02\n", "ffl -23.89\n", "pp1 28.13\n", "bfl 37.36\n", "ppk 14.66\n", "f/# 2\n", "m 0.0006947\n", "red -1.922e+08\n", "obj_dist 1e+10\n", "obj_ang 22.5\n", "enp_dist 20.21\n", "enp_radius 13\n", "na obj 1.3e-09\n", "n obj 1\n", "img_dist 37.32\n", "img_ht 21.55\n", "exp_dist -24\n", "exp_radius 15.34\n", "na img -0.2425\n", "n img 1\n", "optical invariant 5.386\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAb4AAAEuCAYAAADx63eqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdd3Rc13nv/e9UDKag9947wQKCJNiLRIqkGqku2ZIs24ns2K8TO/ZKshzHdnJzHd+bemM7thNbvRdKFCV2igQ7QJAA0XvvHVMwM2fOvH+AlEiJBQQGfX/W4uIiMXPOBgbAb/Y5ez+Pom9EdiMIgiAIC4RypgcgCIIgCNNJBJ8gCIKwoIjgEwRBEBYUEXyCIAjCgiKCTxAEQVhQRPAJgiAIC4oIPkEQBGFBEcEnCIIgLCgi+ARBEIQFRQSfIAiCsKCI4BMEQRAWFBF8giAIwoIigk8QBEFYUETwCYIgCAuKCD5BEARhQRHBJwiCICwoIvgEQRCEBUUEnyAIgrCgiOATBEEQFhQRfIIgCMKCIoJPEARBWFBE8AmCIAgLigg+QRAEYUERwScIgiAsKCL4BEEQhAVFBJ8gCIKwoIjgEwRBEBYUEXyCIAjCgiKCTxAEQVhQRPAJgiAIC4oIPkEQBGFBEcEnCIIgLCgi+ARBEIQFRQSfIAiCsKCI4BMEQRAWFBF8giAIwoIigk8QBEGYF6yjo1wqL+WV99+95ePU0zQeQRAEQfCI0dFRqurrKa+porm9DbPFgiS5cOMe1/NF8AmCIAiz0qjDQW1DPaVVlTS1tV4XcEqFEi+1Bq1ag5dKi1OyogDcQKB/wC2PK4JPEARBmFGjDgd1jQ1jAdfayojF/KWACzD4kpGYSEJoNB0DPZQ2V9M93I9SocAuOQGIDA/n2Ucex6jX3/J8ir4ReXxzQ0EQBEGYBIfDQW1TI6VVVTS2NjNivjbgFHiptfgbfIgLiSIjJomkiFi8NF4MWUbILyvgclM1vSMD+OlNBJn8aeptxy45iImI4tmHH0F/m8C7SgSfIAiC4FEOh4O6piYuV1fS1NrC8IgZySXhdo8FnFatJUDvQ2xIJJkxySRFjgXctfpHBjlRep6y5hr6LUMEGPzIiknBoNNxuOQMDslBbFQ0Tz/8KHqd7o7GJ4JPEARBmBCHJFHf1EBZZRUNrS0Mj4zgvDbgVGMzuNiQSDKjk0iMiMPby+umx+se6OVEWQHlLbUM2cwEm/zJjktlbcZyShoq+ejCMRySk4SYWL66+2F0dxh4V4ngm4TKlg58dFoiggNneiiCIAhTxiFJNDQ1UVZTRX1z01jASdcGnAY/gw+xwRFkRCeTHBmLt5f3uI7d0dfFidICKtrqGRm1EOoTyOL4dNZmLMPgbeB4yTk+LjqO0+UkKS6ep3Y/jE6rndTnIxa3TMJv33qbwKhUfvzYtpkeiiAIwqRJkkRdcxPl1dXUNzcxNDKMJEnIVwJOo9Lgp/dhSUwaGTHJJEXEYdCNL+Cu1dTVzqnyAirbG7DabYT7BbM+Yzl56Tl4e3khyzKflpxl/6V8JJeLlIQEnnpwN9pJBt5VIvgm6ERZIxokEsJDZnoogiAId0SSJBpaWyirqqKuuZHh4RGckhPZ7UZxdQanN7E4OnXsHlxELAZvw6TOWdfRzKnyQqrbGxmVHET5h3B39mpWpC3BSzMWaLIsc6joJIdKTiG5XKQnJ/P4A7vQqj0bVSL4JuiD4yfw8/ElIzp0pociCIJwQ5Is09jcTFl1FXVNjQwND+OQnLjdbhQo0Ko1+HqbWBSVQkZsMikeCLirZFmmuq2B0xUXqetswuFyEh0Yzo6cjeSmZKO5JsxkWeZgUT5HLp/BJctkpqby+H0PoPZw4F0lgm8CLjW0Y1JLSDKkRYsZnyAIM0uSZZpamymvqqamqYGhoeHPZ3CfBZyRrKhkMmKSSImIx6j3TMBdS5ZlyptrOFN5ifquFlxuF7FBkTy48m6WJmehViq/9PhPCo7zafk5XG6Z7PR0Htlx35QF3lUi+Cbg5f1H+c6GWP77RBO6KX6BBEEQrpJkmebWVspqqqhtaGBweAin85qAU6nx1ZvIiEwkMzqFlKg4THrjlI5JlmWK6ys4V11MY3crAHEhUTy6djuL49NRKr9cElqWZfadP8bxigJkt8ySjEwe2nnfl4Jxqojf2neour0XvcJOWpgvslsx08MRBGEekmSZltZWymuqqWlsYGBo8LqA06jU+HobSQ9PJDMmmZTIOHwMpukbnyRRVFdGQU0JTb3tqBQqEkOj+crGB8mISbph2F39vD46d5iTlReQ3W6WZWWza/uOaQu8q0Tw3aHf7z3IsyujARAbQQRBmAxZlmlub6O8poaahjoGBodwOB3XBZyPt5H08AQyopNIiUrAdxoD7lpOSeJ8VTEX6kpp6etAq9KQGBbL1+96hNSohFs+V5JlPjh9kDPVF3HjZvniJTy49Z6bBuRUE8F3B1p7h9DKNpbFjhVAlUVXJ0EQxkGWZVra2ymvraa6vp6BwcEbBJyB1LA40mOSSI1KwM/gM9PDxu50cK7yIhfqymjr78Zb60VyeCw7czeRGB5z2+dLksR7Zw5yrqYYgFXLlnHvlrtnLPCuEsF3B/7j/f08sjQcALvLhQtxqVMQhM/JskxbZydlNdVU19fSPzCIw+lEdsufBZzJ20BKWNzYDC4yHn+T70wP+zo2u50zFRcoqi+nc7AHg5c3qRHx7M67h9jQiHEdwylJvHvqEwrqSkEBa5bnsn3T5hkPvKtE8I1T77AZlXOEtSmJADT0mFFrbl56RxCE+UuWZdq7OimvqaGqrpa+wQEcjmsDToVRZyA5NJb0qERSoxMIMPnN9LBvymKzcLL8AsUNlXQN9+GjM5AWlcBTG+4jPHD8W7acksRbJ/ZR1FgGCgXrVq5k2/qNsybwrhLBN07//PZ+ti0K++zfDd0W9FOwHFgQhNlDlmU6ursor6mmqr6Ovv4B7A7HlYADtUqNSWcgKTiG9OgkUqPiCfTxn+lhj8uwZYQTZYWUNlXRc6XjQUZUEs/e9RAhfndWhtEpOXjj+D4uNVWgUCjYlLeGLWvXzbrAu0oE3ziMjDpwjw6yLXPJZ//XMWQjxD96BkclCIKnyLJMZ0835TXVVNbV0TfQj91+fcAZvQwkBEWRHp1EWlQ8gb63bnY6G/WPDJJfWkBpcw39lkECDL5kxaTwJ1m5E5qR2p0O3ji+l+LmKpQKJVvWrmPLmnVTMHLPEsE3Dr9842M2pgVd93+9Zgcp6aI4tSDMJbIs09XTQ3ltNZV1tfT19zN6bcAp1Rh1euIDo65coown2Hdu/5x3D/aRX1pAeWstg9YRgkz+LI5PZ13m8gmvELU77bx67ENKW2pQqpRsW7+RjXmrPTzyqSOC7zZGJQmHuZcHli657v+HrRKLYsJu8ixBEGaSLMt09fVQUVNLVV0NPX39jNrtyG4ZAI1SjUGnJzYgkvToRFKjEgjy8Z+1l+buVEdfNyfKzlPZWs/wlY4HK5KzWZuRM6mSZDa7nVeO7aGirQ6VSsnOzVtYu2KlB0c+PUTw3cYv3zzAisQvX7N3uGSigmbXaixh7rhcUY7W4EtqTORMD2XO6+rtoaKmhoraGnr6+q4LuLEZnDcxAeGkRSWQGpVIiG/AvAm4azV3t3OyrICq9gYsdhthfsGsychhTfqycbcIuhmb3cZLR96nqqMRlUrFfXdvY3VOjmcGPgNE8N2CJEmMDHTw5PYlX/qYqNoiTMbhT49Sa9Xzqx88N9NDmTN6+voor6mmoraG7r7esYCTrwacCqOXnmj/MNKiEkmLjifEN2heBty16juaOVl+gZqORmwOO5EBIWzKziMvbelnHQ8mwzJq46Uj71HT2YRarWLXPfewYsmyyQ98hongu4V/33uCrMgbXwMXe/iEyTB4qVEO2+keNBPiN7W1FOeanv4+KmpqKK+pobuvB7vdjuuagDN46Yn2DSU1MpG06ARC/ed/wF2rqrWe0xVF1HZ83vHgnqXrWZG6+LqOB5Nhtlp48eh71HU1o1ZreGj7TpYvXuyRY88GIvhuoaO9kR88mnnDj4kZnzBZAbj5j3cP8g9f3z3TQ5kRfQP9lNfUUF5bTXdPL6P20esDTutNlG8oqVHxpEcnEuofvKAC7qqxjge1nK28SF13Cy7ZRWxQBA+svItliZke7WQwYjXzwuF3aehpRaPR8Oh9D7A0M8tjx58tRPDdRGlzJz5aBV4q1Q0/LoJPmCx/rYaWgc6ZHsaU6xscoKKmemwG19uDbfT6gNNrdUT4BpMWFU9qdCLh/iELMuCuJcsyxQ0VnKsqprGnDbfbTVxIJI+u3s6ihHSPF3UesozwwuF3aeptQ6vR8uSDu1mUlu7Rc8wmIvhu4qUDx3ly6Y3L84w6XLgVIviEyQtWOHjnZDEPr537l5EGBgfH7sHV1dDZ3X3DgAvzDSI1LZ70qCTCA0XAXUuSZS7WlnK+uoTm3jaUChUJIVF8ZcP9ZMQkT8nXatAyzB8PvUNzXwdeWi1P7XqYrNRUj59nthHBdxMKh5mVCYk3/FjroBWlevI3jgXhiVhv/r2gYE4F3+DQEBV1tZRXV9HR043NZkOWZdyA6mrA+QSRnBpHWlQiUUFhIuBuwilJFNaUUFBz+ZqOBzF8bcvDpETGT9nXrX9kkD8efpfW/k50Xl4889CjpCcnT8m5ZiMRfDfwwpECogNuXoezdcCKbpLLgwUBIEivQes0z8pFLoPDw1TW1VJWXUlndzdW2ygu2QVcCTiNjhDfAFKS40iLSiIyOGza+6rNRU7JwdmKS1yoK6V1oBudWktyeBw7l28gMSJuSs/dN9TPHw6/S/tgN95eOp577AlS4m/dUmg+EsF3A5cry/nZA0k3/Xj30Cgm4+wtOCvMLXm+M7vIZXhkhIq6Gsqqqujo7sJqG0WWXZ/N4Lw1XgT7BJCcNNZRIDI4XATcHbra8eBifQUdg93or3Q82JW3jdjQqd/L2T3QywtH3qNjqAe9zptvPvkUiTFxU37e2UoE3xd0D5rxVkkE6m8+4+szOwj2n/leWcL8cFe0iRPlU7/IZdhsprK2hrKaKjo6u7DYbJ8HnEKJt1Y3FnCJcaTHJBIVHCECbhIsozZOlhVS0lhJ51AvJp2B9MgEHt+wk8jA6an61NnfzQtH36drqBeDtzfPP/VV4qJv30dvvhPB9wW/2nuMdcm3rs03aJNYlCRmfILnBCscvHOqhIfXZE/6WMNmM1V1tZRVV9He1YnVasP1xYAz+bMycTHp0YlEB4V7dEn8QjZsGeFkeSEljWMdD3y9jWREJ/HM5l2E+Afd/gAe0tHXzR+PvEvPSD9GvZ5vfeUZYqOipu38s534bv8C61A3O+659b4Vy6hEelTINI3oxkrrGulqbwTA1+TD8nlQTWEheyLWm38rOH9HwWe2mKmsq6Osuoq2zo4rASfjxo1KoUSn9SLYGEBuXDYZMYnEBEeIgJsCAyNDnCgtoLS5mn7LEP56H7Jik/lmZu60tyhq6+nkhaPv0WsewGQw8mdPP0t0hCiL90Xip+Aahy7VEGRUob7J3r2rHC6ZML+Z68W3/0Ilx/OPkRviTYtVSafZzNsdAfzT9rgZG5MwOUF6DV6OGy9yMVutVNfVUlpdRWtnB1arFZfrmoDTeBFk8ic3Lov06ERiQqI8VsFDuLGeoT5OlhZS2lIz1vHA6Mfi+DTWZeZOuOPBZLR0t/Pi0ffoswzhYzTx3WefIzIsfNrHMVeIn45rHDh7nu9tjL3t42QUM/bO+URpPZ+eOMbP1kVi8tIA8NvCEQoc8oyMR/CclT4y//raO2SEGWntaMdyg4ALNPmRE51BekwScaEi4KZT50A3J0oLqWipZXjUQohPACuSs1mTnoNxhppSN3a18dLR9xiwDuNr8uF7X/sG4aHj75i+UImfmitGJQm1bCMp5Pbv1mT3NAzoBorq2/no0AH+Nu/z0APoGzHjYzvOH149Oi3jsI46+c7XvzVlx88/d5qTxaUMuTS3f/Adcstu3OZeVO5rX0Q37d5JvPi9hz1+vttp7Oqior2bfpsDpxtU+hAaRo0EGv1YmppGRkyKCLgZ1NLTTn7ZBara6q50PAhiTfoy1mTkTLrjwWTUdzTz8rE9DNpG8Pfx5QfffJ7gwLndN3A6iZ+mK3619wRZEeO7ROGegXJlla09vP7BB/z1yjACDdcHgtstow+PA+P0jKu7ro43int4fHHwlBy/p6+HZ9fFERc6NR2uf3+0EkdnK8/EjO1JK+q08D8jw1NyrhspuHSJlq5uQiw9nOiTMam1pAeFkR4ayUudNnJyt7I1deFsJp5tGjpbOFl+ger2BmyOUSL9Q9m0KI+89CV4aW6+2ns61LQ18urxDxiymQn08+OHT3+LQP+51wl+pongu6K9vZHvP3TjgtRfJE9zZ4amngH+++13+cHyUMJMuhs+RnLDt3KmZwP077qlKQu96fDNzWn880d23mvvYnfE9E3fT1+4wCdHDyO5XPhqdGgDFvFcwvXbYnYGK3m5pFAE3zSrbm3gdEURNZ1NOCQH0QFhbFu6jpWp2WhmQZWm6tZ6Xj3+IcOjFoL8A/jhM98m0G96F87MJyL4gPLWbkwa0GlvvajlquksUN0xMMy/v/IW310SSKyfqBbjKT+4dzH/8G4BR7v6meqNKSfOnuVA/qfILpn0yES+sulBXjv4Olbzl/fIJRh1uLunb/a5UMmyTEVLHWcqi6jvakFyuYgJiuD+3M3kJGXNmtWv5U21vJ6/F7PdSkhgIN/+2tfw8xUNsCdrdry6M+zF/Z/yyOIbF6T+olGHC3maClQPmG38nxff4OuZfqQGza5yVvPBXz24jL9/+zxJ0uiUHP/IqXyOnjqF7JbJjk7liY33jutSmS8Sp+ubWZ0gNhp7kizLXG6s4mzVJRq6W8c6HgRH8vDqe8hOyJhVm/VLG6t44+Q+LHYbocHBfOeRr+PnI4pmeIoIPgD7CGuSx1evrrnfgmoaLn2MjDr4+X+/wtMpJpaGi2/4qaBWqfjR7lx+9spxFEgeOaYsyxw48Sn558+B282SuAweX7/jji6XbQ3w4v2L50TweYAky1yqLeN8dTFNve0oFQriQ6J5at19ZMalzLri2cX1Fbx16mOsjlHCQkL43iPfxMc0/dsj5rsFH3ynKpvw14/vEidAS78NnW5qLzmOShI/+e3LPJxgYGW0qBAzlQxaNesyIjlc0kJ1ey8pEROrriHLMh8dOcTZoiIAchOzeGjN9gmtxsz0M/BWw+CExiGAJEkU1FymsPYyLb0dqFVqEkOjeXbLblIjE2Zd2AEU1Zby7ukDWJ2jRIaG8f1HH8NoEFd5psqCD769pwt5KGP8VVj6zKMY9VN3jV2SJP761y+zPcqLTfHi5vV0MOpUqFxKfv3aW/zkT54hwGf8e7JkWWbPgf0UllxCgYK8lKU8uOquSd8jMrmdlLR3kR0h9mSNh1NycLaymKK6Ulr7u/BSa0kOj+WbWx8jOTJupod3UwVVJbx/7iA2p53o8Ai+/8hjGPX6mR7WvDflwZd//hwtHW1TfZoJcw93EGGIwWK1YdDffiY3aHPia5y6d2J/9dvX2BSq4t4UsSdnOimAh4Jc/Pz3r/CL730d3W2CS5Jl3t23l0vlZSgVStanL2fnyrs8dp9og6+WN8+fJvvBXR453nxkd9o5VVbEpYZy2ge60Wu9SYmI54FVdxMXOrvrUp6rvMSec4cYlRzERkbxzEOPoBeBN22mPPiOnTyJYpqX/4+XS3ZjULs4XFA+1ilaASqFAq1KibdOg7/RQHhwAKEB/gT6+qBSKRkZlQgPm5rg+6c3D5Coc/BwhqitNxOWhRjotZv561+9xL9+77kbPkaSJN7c+yGl1ZWoFEq2ZOaxPXeDxy+f5QYa2dvQ69FjzgeWURunyi9Q3FBxpeOBnrTIRB5bu5PI4OnpeDAZZ8qK+KDwCHbJQXx0DF996BH0uhtvURKmzpQHX3RIBM9ve2KqTzMh/3H8JGFefTyxOA63282IzUr30CAdg/30jgzT02+hpXsQSZaR3W6UCgVmq5KLnT2Y28rJTs8gPTkFnXbyi13O17Qy3FnP326Om/wnJkzY1mgjA3Uj/Og3r/LLbz312f87JInX97xHZW0taqWKbYvXcffSNVN6v8jglqjt6SMpeGHP/kesZk6WFVLcWPlZx4P06ES+unkXYdPY8WAy8i8X8FHRMZySk4TYOL666yF0IvBmzIK+x9fU08w370oBQKFQ4KM34KM3kBT+5RmXQ5LoGx7iV2eqCfXyorm5laq6OmS3jFKhQKPREOjvT2ZyCtnpmXdUPmhUknhj3z5+vGL2v2NdCB5LNPH/Knr49z3H+NMda3j1vXepbaxHo9Jwb85GNmavmpYFEnlGFa+cO8NP7713ys812wxahsm/XMDl5ir6zGMdDzJjkvlm5nICfedOpZJPS87xSdFxnC4nSfEJPLXrIY+8URYmZ0EHn0bhxOg1vnqQWrWa8IBAZJU3f37fU58tXrDZbVS2NlDSUEFzbwfHTp/h0Ml8ANRKFUajgfjoGJZkZpEYE3vDRQ8//Z932BzhfdOqLML0+1aynn8sO89PK0/jpdawa8VW1i3KndYxrAk0crh56hvUzhZ9Q/2cKCukrLmaAesIgUY/FsWksm5RLn6GubOlR5ZljhWf4UDxSSSXi9SERJ58cBdaEXizxoINvk8qaggz3XndPdl9fWcGby9vliZmsDQx4/PHyDItPR2UNFRQ3d5IRXUNl8rKcONGqVCi89ISEhRMdlo61f0SPtIQu9PFnq3Zoq6zk9P1rUQrNbT6pvKLp789I+NQqZToZAftQyNE+M7PvVydA73kl56nvLWOYZuZEJ8AchKzWJu5HJN+bi3nl2WZI5dOc6jkFJLLRXpyMo8/sAvtLKkCI3xuwb4ix6tKeG7p+Kq1XGs8lR2VSiWxoZHEhl5/yXTQMkxZYw2lzdV09Haz9/BhvHCSqrNT1gAZcVEoFLNvj9FCYbbZ+KS4HKsMKX6BrE3M4IWWQf716FH+YvPmGRlTjl7FC2dO8zf3bJuR80+Ftp5OTpQXUNlaj9luJcw3iNWpS1mdkYNhivfITgVZljlQdIKjl8/ikmUyU9N4/L77Z03ZM+HLFuwr45IspIck3fHz3JNYoepn8GFNZg5rMnMA+P6bL7M7Uk1vXwunapo5VdNEiFHHuswUAkU9vmkjy6CVR3mnqBk/rZYnMrMxXmk58/UYf35aW0ltz+IZWWRyV4iJky2t035eT2vsaiW/rJCa9gasjlEi/EPYmLWS1RlLZ7zjwUTJsswnBcf5tPwcLrdMdno6j+y4TwTeHLAgX6GmgUEMmokFmKcKVP/y0CEyjLAyMhQiQ8HtprqjnfMt9bx9rgStUkFaeBDLU5PRqsdfWUa4M61dPVyqrkcrj3J3fArxQV/eMP5smIH/PPgB//bUjbc4TCWVSonW5cAmSXjPsV+oNW2NnKq4QG1HE3bJQVRAKFuXrGVl2uJZ0fFgomRZ5qNzRzlRWYjsllmamcXuHffOqlqfwq3NrZ8kD3mtsIAV0RNbGeaJJjbnmlsY7mvir1ZeUx9UoSAlIpKUiEjsTidn66qo6OjiclsPvl5q8tISiQ0df4UZ4dZcbjf7TxXQ3D1IgMFI16D2hqEHEGf0In1gkF8cOMBfbZv+S44RKjcfFJfxeM7iaT/3nZBlmarWek5XFFHX1YzkchEdFM59uZvISVo052dCkiyz9+xhTlVdQHa7Wb4omwfu2SECbw6a29+JE9Q32MU9eePrvfdF7kkmn12SeOPkAX6y+Oab1L00GjakZbGBLDoG+zlVV8X+4iqUVBEbYGJNZioG77l3L2S2aO/p5ePTYyXGduWsor53kIq2+ls+56loP/6+rm5GyohtCNTzekPVrAw+WZYpbazmbNVF6q90PIgNjuShvHtYnDi7Oh5MlCRJ7Dl7mLPVF3HjJnfJUh64e9usrPkpjM+CCz67JKFVutCoJvZNO9kmtH+3dw/bI3wI1I/vvka4XwAP5+ThkmUuNTdQ0t7Cy/mF6FUKlsZHkRUvVoOOl8vt5sCpQpq6B4gPCmL74lyUKiX1veMrCP2NCCO/OrqX//v4s9M6e0ky6ZD6RqbtfLcjyTLFdeWcqy6mqacNhUJBQkgUT6y7l0VxqfMmECRJ4r3TBzhXWwJAXk4OOzffNW8+v4VswQXfO8WXiQuYWE08u1OGSQTfa4VF+MnD3BMbd8fPVSmV5MQlkhOXyIjVSn5tBWdrWzlT24IDPWppFJifS949oaO3n49PFQEKdi1bQWTQnXeQD/fWkq218r8PHORvd+7w/CBvQe1yYJckvGbocqEkSVyovUxBzWWaeztQq1Qkhsbw9KZdpEcnzqswcEoS7578hIL6UlDA2txc7tm4eV59jgvdggu+S43V/OWaic2Susw21BO8KS9JEhdqivhlbtyEnn8tk17PjuwccLup7e7kxepGwkdK+P0eN9EhvuRlZ+BrHH+HgflMdrs5eOYCDZ39xAUGsWPJ2Cxvoh6N9OMf6hopbG5lecz0FUIOU8l8VFbNQ4szbv9gD3FKDs5VFnOhrozW/k60c6TjwUQ5JYk3T+zjYmMZKBSsX5XH1nXrReDNQwsu+JSMEuE7sRlfl2V0wkuv/19+Pkv8vCd8ifWGFAqSQsPRN7TT5ZfLbv9uSlubee3gKbRqFUmRwazISsd7nNVpZgtdQDR/OFPMzx+cXGmqzr4B9p26gNut4P5lucQEeWZx0PORJv7t+McseeK5abvkuc7fm3dqy6Y8+OxOO6fLL3Kxvoz2gW68tTpSIuK5f+UW4sOip/TcM8UpOXj9+D6KmypQKBRsXr2Gu9ZtmOlhCVNoQQVffV8fRs3Eg6ff6kCnvfNFJXZJoq2zju97YLZ3UwoFy5NSWZ6Uis1u53xdFTVtHZQ3HcNbqyYrIYrFqUkeCV6zwxNrW29ux6plFNa1888Hq/jB1tQ7fr7sdnPo7AXqO/qJCQhk59IVqDz4hsOkdpMgD/F377zIkohwTHoDJm8jJm8DBp0Og7cRk84bby9vjzMO0oIAACAASURBVM0W0nz1OAeGPXKsL7LZbZwsu8ClKx0PjF56UiMTeGTtdqKD77zIw1xhdzp4/dO9lLRUoVQouWvdejavXjvTwxKmwYIKvn1lFaSGTLzm39CoA5PXnZdR+r+Hj7AmyIDCk7O9W/D28mJDRjYbMrIZspo5U13JpZoWCqsaMXlryUlNICUuGqXizu9X7itpp3506ktJ/eSpe/nR79/k5TONfDUvbtzP6+ob4KNTRbhluG/pcmKDJ7cC0+V20zXcT31vD11WM2anE4fsQqvSYJdc1LTZcbldOF0uXG4XLtmFyy3jkuXPStSpFMqxv5Vjf9RKNc7REXxGrewZdqBSKFAplKiUClRKFWqFAo1SiUqpQqtSoVaN/a1yqGnpafdIGJmtFvLLCylprKR7uB8fbyPpUYl8dfMDhPnP720zdqedV459QFlLLUqVku0bNrN+1aqZHpYwjRZU8LX0tvPo6olfrjHbnfjq7uwyqdnuYHCghYdWJdz+wVPAV2/kniXLAejs7+dMbSUnLlVx/FIlASZvVi1KIzp0fAs9KjuG+Likh3/89jNTOeTP/OPXHuJH//Uq4b6d3JVx684VstvN4XNF1LX3ERMQwM6lKyc0yxuyWWns66Z5aIAhpwOb5ESlUOKj1RHq7ceK8GCifQNRqdS83dGPJSiF76648SzB6ZKw2K1Y7LaxPw4bZruFEZuF4sZiNC4l/lonksuFhIzLJeOQXbhkGZfbjYyMLLuR3W5cuDEqvfnXD19Eq1IRYPAlPTqJ1RnLCPTxH9/nZhkhv/Q8JU3V9JkH8debyIxJ4Rt3PzKnOh5MlM1u55Vj71PRVo9KpWLnlrtYm7tipoclzIAFFXxut2PC9/cArJJMuPeddVD45cH9bAufHZXlwwIC2LViNW63m/quDgrqqtl3+iJKhYLwABOrszMI9LvxWIetDn59tI7vPPbwbbuTe4pareYfvvEEP/7dS4T56siK9Lvh43r6B/nw5AVkGe5dspy4kPHN8lwumZa+Hmo62tBJg/zhQjMuZAwqLYHeRnL8I4jzC8Z4kzc7D4X48tPWSlyuVahUX/6aaFRq/PQ++Om//DXt6W/DolKwOWL8l84rhi3sdeh5PD2Tc42XuVhXxqfl59GqNAQZ/VgUl8qqtCX4GD5f3ds3PMCJ0gLKmmsYsA4TaPQlKyaFdZm5+JsWRlk8m93GS0fep6qjEbVKxf1bt5G3LGemhyXMoAUVfBqFPKnnjzpd+N5BEd1eiwW7uYutWYmTOq+nKRQKEsMiSAyLQHbJlLY0cqmpgbeOnkWrUhIbHkBeVjoG/djnKrlc/HxvOXev2UBC6PTODPQ6Ld9/8hH+7bU3+fH9GYT6fP7GQ3a7OXL+IrXtvUT5+nPfslWo1Def5dmdEnWdbdR3ddJrNmN1OtCoVEiSGlnpxf1JSwg2+qEcZ6FwpUrFIq2SX184e9NZnyel+xh4q3WIzKgUMqPG+kg6XQ5KWqopbCzjVMUFDhTn46XSolVrkGQXNucowaYAliVmsC4zd851PJgMi83CS0ffp6azCbVaze57dpC7ZMlMD0uYBRZM8JW0d2HSTq7mpV1yEWgY/4zxnw9+wq6Y2X0JSalSkh2XQHZcAg6nxIX6asrb23ipNR+dRkVaTBiHGlwEhyewddmdLzTxhKggXx7dvoNffPQx//vhxei0KnoHh/jwRCEu2c2ORctICAv/0vOGrVZq2ltp7O1h0Gpl1OXEW60lSG8iJyqJxOBITDodp+sbqO+pItR056/V7WZ9nqZyOXC5pM/OpVFpyYnLIsI/hFM1FynrqGXANoxCoRi7z+h2Yxm10djZSqDJn2VJGXO6TuZ4mK0WXjz6HnVdzajVGh7aeS/LF82+qjfCzFkwwXeospyMsMld2rFLMiHG8QVf08AgSvsgqyNm5t7eRGg1avJSM8hLzcBsG+VcbQX5lV0MDku4uzv4q1+cGdu+r1CgQIFSoUClUqFSqdBqNOh0Ogx6PV0uIyvj/PH39cVkMOJn8sHX13dSnadXJEfR3L2Sn+wp4J44N1VtPUT6+nHfslWo1Srcbjc9w0NUtbXQ1j/AsN2GJMsYtV6EGP1YnJRAXGC4xzeAK1Uqsr1U/KrwDP/fynUePfaNhKjcHG1u5O74JFr6OjldV0R5Zz0jdguxAeFsTl1JbvwidFe23YyMmiloKKW4tYoPzx3hrdMfY9QaiAoIITdtCVkxyXO+huZVI1YzLxx+l4aeVjQaLY/d9yBLMidWmlCY3+bHd/w4dA528Y3syYWQJMsE6cd3qfNXR/bzdELQpM43k4zeOrTGKLrQ8i/PPYksy1hGrQxbzAzbzAxbzQyYhxk0DzFks2C1WxkdtdMz0kOvNMSJjkpktxu3W0b+QoFTBQoUirHgVKqUqFVqNBo1Xl46jN56fExGAnz98Pfzw8dowsfHiL/Jl8URJs6f6uN8nZtHcxajVqk5VnqJzuFhzHY7CgX46vRE+ASwPiSSCN9AVNOw+Xh3sA8/ba3C5cqb0lmf2y2TrpI4UHyYEyUfY3WOEhcQwbaMNeTGZ6JRffmNhUlnZHP6Kjanj61a7DcPcq6+hNL2Gt488RGvyBI+OgMxwZHkpS0hOTJ+zm3YHrKM8MfD79DU246XVsuTux5iUWraTA9LmMUWTPAp3E4CDJPr++WGcb07LuvsRucykx08d5eFD1odvH25hb9/8ClgrLmuSW/EpDdy8/LaY75/spFfro277v9kWcZmtzFsNTNkNTNsMTNgGQvOYZsFy6gVm91On6Wfzq4uJNfYtoAvBqcKBX1WJZ9cLsFbo8HfW09qcCTJoVEE6E0T2qIxWVM56xsZtVLb30nrcC9dNjMqpRJZH8X9K7axJDodzR0GbYDRj+3Z69mevR6AzqFeztZdoqKznv85/A5u3PjqjCSExbAmI+dLzZRnk0HLMH889A7NfR3otFqe3v0wGSkzczlemFsWRPBJkoR6kgtbANzj7MX3wolD/FnK9Fbw97R/OVnJ/cvWY/TyzP0gpVKJwduAwdtAeOD4vzYnSwv5oOAwOrWWZ1c9iJ/Bh4ruTg7UFfJPG3M9MjZP8NSszylJNA510TTYS4dlEJtLIkhnINoUwOroVAINvvxDcz+5cYs8Mu4w3yAeXHYXDzL25qS1v4MzDcVUdzVxsbEcBQr89T4kR8SxJiOH8MCZfzPXPzLIHw+9S+tAJzovHV975FFSE5NneljCHLIggu9McysB+sn/Ah9PvZIzjc2YGCXB78uLLeaKo3WdjKJnU/LMrUbtHxnkVx+9wqBthE3Judy/+PMiwaE+QVzs7eOFknqeneTla0+5Ouv7z4IzfG/V+Gd9brdMt3mQuv4u2syD9DusmNRehBt8WRedRrRvECrl9YuyVG4no5IDnYcXqSiVSmKCIokJGpvlybJMXXczZxtKqGyp42zNJdRKNQEGXzJixvYQBphuvMVkKvQN9fOHw+/QPtiDt07H1x9/kuS4+Gk7vzB/LIjgy6+tITd88j+g4wm+d8/n85fpczf0nC6ZD8vb+NkDX5mR88uyzJsn9lFYf5kwYxA/v/+7+Oq+vAT/L1Zt4Af7X6V1xEqUaeJ7Mz1pd7APP225/azPJcv0WYeQBzr4Q/cgChSE6X1IDwwnISAMvfbWe0WDVUpONDezNSHJ05/CdZRKJclhcSSHxY2N2yVR3lHH+YbLXKgt5VjZObxUGgKNfmTHp7EqbcmUbJfoGujlxSPv0jHUi17nzTeffIrEmDiPn0dYOBZE8PUP97Bh5eRvdrtv05LILkmoZRsRxrkbfP+SX8nyxGUeu8R5J2paGvjD0XeR3TJPr7iPnNtczvvT3Hv4z4KP+cXmZdM0wltTqlRk62496zvfUMJHJZ9iHx0hTG3g/th4ggy+KMa5dxAgVafhYkfjlAffF6lUahZFpbIoauw+ml1ycLm1isKGMk6WF7L/0gl0ai9CfAJZlpTJ8uRsvL0mfl+9o6+bF4++R9dwHwZvPc9/5WniouZnoWxhei2I4FMpXHhPcg8f3P4e3x/PnifdZ3ILaGZSSccAnVb4y5zp3eRrd9r5zb7XaO7rIDs8mWfW7BrXoo2UwECCfGJ48XI9zyyaHZc8dwf78PP2OuD64Ctvr+P9i4cYtlu4Ky2Ptu5GrBYFwcY7L3q+zM/AsZ5uD4144rzUWpbHLWL5lTcoNucoRY3lFDWXc6Aonz3nDmHQehPmF0xuyiKWJKaPaw9hW28nLx59j56RAYx6A9/+6rPERM7eRTbC3DPvg88uSaiVk1/YAre/1FnTVss/Lrl1TcnZ7MULDXx/68PTes6jl87wcdFxDFpvfnjXc0QH3tnX7/urNvD9/a+yJW52XPJUqlToZQfNQ4PE+PrR3NfOO0UHaR/qZk3CUnZkb8BLreXF7sYJn8Og1oDL6rlBe4i3Rsea5GWsSR6bgQ+NmilsuExxSxXvnz3EG6c+xqj1JiogjBXpS8iMTUF9zdaJlp52XjzyPn2WQUwGI9955jmiwufu1RNh9pr3wXe4up4w453V15yIXosFb7cDb+3c/JL+19kaIoLjifCdni7uPQN9/PqTVxkZtXB3Wh47szdO+Fh/mruN31zYz//auNRzA5yEdSYv/lh4nFC1jeruJpZGpfGn6x/F4OW5YFa6nB471lTx1RnZkp7HlvQ8YGwP4dkrewhfP74XSZbw0RkJNvnTMdzLyKgFX6OJ7z77HJE3qMQjCJ4yN39L34HCxlq2xI2vev3t3OpS5+/yT5IXPD2h4WldI6NU9Fj5l0c3Tvm5ZFnmlWN7KG6sJNI3hB9tfQ7TDRav3InUwCAkDJR0D5Ad4pnXeqIGrCMM9zUwJCkIjUzlb+75JkETKIV2O1pc9FjNBM+h2psBRj92ZK9nx5U9hEVNFbxdtJ/q7qbPHmOz2zlw/FM25a0hPiZmpoYqzHPzPviGrf2sjvFM2aJbXTDtH2zn3uVz88b7787Xsmv5+ik/T1ljNS8f/wDc8FzebhbHeK66xp+t3MJvz304Y8HXOdJPQVsdnbZh4kyBGHVBbMq6a0pCDyBKreRoYyOPZWRNyfGnUmNvGy+f/YAeyyDhfsH8+f3P4G/yo6GrlZNlBdS0NfG7114GFHh5aUmOT2Bz3hrCQ+f23lhh9pj3wadWuD3SdVx2yXCTVZ2VXT34KOVpazTrSYNWB702mXUJcVN2Dpvdxq/3vUrbQDfLItP4at79Hi/tFWXywYmest5BMoOmb29ZY38nFzobGLDbSPEPZUvCIvRaHUGDI7xWfJqfbb5/Ss6bbtDxaVcTzKHgq+tq5pXze+mzDhHpH8qPd377un2AieExJIaPzfJkWaaqtZ78sgJq6uq5XFmBQqFAr/MmIyWFjavyCPSf3QXghdlr3gefSjGe3Xe3Z5VcKG5SDuvlsyfZFT2zl9gm6vcFtaxPm7reZAcKT3Co5BQ+Xkb+Zts3CPObusofz+du4n8K9k35vT7ZLVPZ00pxdzOjLomMgHDuT0lAo9Z89pilfiY+bO6dsjEsMun5oHNgyo7vSTVdTbx6bi/9tmGiA8P4s/u/ip/h1j0qlUol6TFJpMeMbdmQJImSxkpOVxRRUlpGQfElFAoFJoOBJRlZrF2xEh/j3LnsK8yseR18NocDtdIzwWdxuFAqb/zlctoGWHFlk+9cYnfKtI44+P6iDI8fu6Ovm9988hpWxyg7M9dzd+Yaj5/ji+J8/XGgo6p/hNQAz99vdckuijsbuNzbhhIFi4OjyQyN+VJllat8kLjY1cHSUM8v1FCrVChlyePH9aTKjgZeL/iIAdsIcUERfPeBZ/A1TOx1UavVLEvKYlnS2AzX7nRQVFvKmcqLnLlQyInzZ1EqlPj6mMhdvJRVy3LQ66Z+UZswN83r4LvY2omvzjMbsa1O5w0r/R+priNYO/2FkT3h9wW1ZEV7tm2LJMuMdlzg/1RVEecfwbd2PI73bSqReNLXl27ilUv7+fsNntuLaJccFLXVUzHQgUGjJS88ieSg8NtuOt/hZ+Cd4jMs3brbY2O5lkqenSs7y9treb3gY4ZHzcQFR/HnD37N4xVdvDRa8tKXkZc+tnXCMmrjXOUlCmpLOHoyn4MnPkWlVBLgH0BeTg7Ls5egnSftl4TJm9ffCWVd7YSbPPNL1+Z0oVR8+Z39gZIC/iwp2CPnmE6yS6a2z8I/b1nusWNeqivn9fyPUOHF8+seJSNieiuLwNim9lHZi/qBERL8JzfrM49aKWivo3aohyCdgbvjsoj2G/9rnepjwN7cN6kx3Ipe4aa6r5eUwNnR/upyazVvXviEEbuFhJBofrDr6xj1hmk5t0HnzeYleWxeMrZ1YsgywqnyQi7VV7Dv0CE+PHgAtUpFSFAQ61euIist47o9hMLCMq+Dr32wn61xnnmnaXO6UKmuDz5JklBIFuJ9Z75i/Z16ubiR6JA4jxzLYrPwq32v0DXUx4qYLEp0i8mIiPXIsSfi2WUbeLH4ED9bP7Gu2/3WYc631dJsHiDS4MeDyUsJNk7sHm4gLvKbm1gX4/mvR5xWRX5zw4wHX3FzJW8V7cfssJIUGssPN38Dg/f0BN7N+BpM7MjdxI7cTQD0DPVxsqyQy03VvLX3Q9748AM0ajWRYRFsWLmK1MTEOdeHUJi4eR18ZusImSGeqd5uk1yov3Av5/WLJSQYNDd5xux2qX2QX+y6b9LH+ejcEY6Vncff28SPdzxPsCmA0kvNHhjhxKUHBmOR1LSNWIm8g2ou7cN9FLTX0W0zE+8TyGNpK/H1ntwbp3sDTbxafn5qgk+n5fjQ1C2guZ2ipnLeKTqAxWkjJSyOp7fswtvrzkuwTYdg30B2rd7GrtXbAGjr6eR42Xkq2+p56d23cePGS6MlLiaaTXlrRE3QeW5eB5/bLeGv90ztTLvkQv2FJfglDRX8OGvu7S36oLwVX1P4uJrq3kxLTzu/3f8GdsnBruzNbExb6cERTt5dSUt5u7KYP89Nv+1jey3DeDn7+biuibSAMLYmZeOt9swl8jiDN1Lv1IRTvF7HR/3DU3LsWyloLOW9i4ewOkdJj0jgqU0PTqoY9UyIDA7jyY1jW01kWaahs5n8sgvUtDTx27qXAAU6Ly9SEhLYuHo14cFz7+dcuLl5HXye2soAYJNk1KrrF8oo3A78vefejO9scy9/ue2xCT1XkiT+59BbVLU3khgUzZ+ufxSdZvb90tuakMyRmnO3fExTXxenGsppHbLgUup4ZtFSNFOwACJI6Z6Sy51+Wg2yNOLRY97K+YYS3rt4GJtkJyMyka9sfgCvWfja3ymlUkliRByJEXHAWBBWttSRX15AVW0dxRXln+0hzExNZcOqPAL95ub2JWHMPA8+zxSnhrEZ37UdA+yShMbtueNPF6dLZtSlxE9/5zOa81XFvHN6PxqVmu9seIKUsNndBFSh8qZpyEKs7/X3m9oGezlZX8aAzcLi8Dii/XS8P1A9JaEHkO2t4XRzzZRc7lTcsp6QZ5yuucgHl49ilxxkRafwxMb78NJMf9uq6aJUKsmITSYjdqyru1OSuNxQyamKC1y6XMr5SxdRKhQYDQaWZi5i3YoVGA1iD+FcMm+DT5IklB6c8TldLryuaalyuLqeUN3kWx1Nt4+q2gi+w+7ww5YR/nPfK/SZB1gdv4RHcu6ZEwsBtiYt5Z2qC/xgxdg+xc6hfk7Wl9FjGSYrPJYHs9fgpVZzur5hSseR62/kSGfXlBxb6XZNyXFlWeZUbREflR7HLjnIjk3liQ33jqut0HyjUatZlpzFsuSrewjtXKgp42zVRU4XFHD83BmUCiV+vj6sWLyEVUtz0Ik9hLPavA2+ko4efLw8dxnS7nKju2ZGUNhQw91Bc+9dXkFLP9/dsmtcj5VlmT1nDnKq6iJBBj9+sv1bBJrmziWezXEJ7K86Tc/IECfrS+kYGSAjJJp7s1ah00zfJWqdSo1SnppLkkrZs8EnyzL5NYXsKz2Bw+VkaVwGj67fviAD72a8NF6szljG6owrewhtFs5WFVNYc5lD+fnsPz62hzDQP4DVOctZlr1Y7CGcZebtq1HS0U6Yh/bwATglGZ3u81+Ww+Z+VqbOrdYpskvG6oQQ0+33tzV0tfLfB97E6ZJ4dOm2z3qszSU2xyiu0UFeLTpBVlgU2zKWo5+he1Iq2YnLJXm8RqkCGbtLwmuSx5VlmWNV59lffhLJJZGTkMlDa7dP2eXf+cTgbWDLktVsWbIagEHLMKfKLnCpoYK9hw6y5+B+1CoVocEhrFuxiqy0NLGHcIbN2+/q1v4e1kd7bi+RQ5bxuaYWowbJI8Wvp9OB2g78TbdeneaUHPxu/5vUdbWQFhLH19c9fN0l3rniQlMpb184iI/OhF0bwF2pE9vT5ylBSjje0szmOM92ijcqFFT29rI4dGINkGVZ5kjFGQ5WnEZyu8hNymZ33tZJrfhd6PwMPuxcsYmdK8b2EHYP9nGytIDLLTW8tXcPb3zoRqPWEBURzoaVq0iJF3sIp9u8/e4esg6TGeK5xQQOlxu9diz4Bq2jaKdhUYGnnWrs5RsbHrjpx0+WFvJBwWF0ai1/sfmrxAfPvb1MFruVl89+SG1PC/ct2siG1Fx+dOCVmR4WiwxenGuu9XjwBaoUVPXdefDJssyh8lMcqjyDyy2zKnkJD6zeKmYiUyDEL5Dda+9hN/cAY1uB8i8XUNnRwItvf76HMD4mhk15a4iNiprhEc9/8zb43LKTEA92XpdcMvor3dX3lVcQY5xbsyDZJWNxQrS/75c+NjAyxH9+9DKDthE2Ji/ngcVb5uQ70ILGy7xTdJBwnyB+vONP8NOPfa4atXHKCleP1zJfAwc7uz1+3FCthvqh8ZdFc8kyB0rzOVJ9Drfbzeq0pdy78i4ReNMoOjiCJzePvQGVZZm6jmbyywqoa2nmN3UvokCBTudFakIim/LWEBo890oiznbzNvg8uaITwCnLGLRjYVfW2sDTUbduqzLbHG/qxmi4vrSVLMu8nf8x5+tKCDUG8fN7v4Ovfu51kbfYrbx4+gMa+lu5P3sz65Kvb7N0T/JS9tacJ3Wl57tQjNfYAhfPF5WO1mkpMt++PZFLlvm45Dif1p7H7YZ16cvZvmKTCLwZplQqSY6MIzkyDhgr8l7VXEt+eSGV1TVcKi9DoVBg8PYmMzWdjStX4e83ff0m56t5G3ye3MMH4HK5MV4JPod9hIzAuXUZ8GRDD48s3/rZv2vaGvnD4XeQ3TJfyb2P3PhFMzi6iTtXX8y7Fw8R6RfKj7c/f8PgXh0Zyf5q+wyM7noq2emRhSjXitd7YR2x3vTjLpfE3pLjnKgtBGBDxgq2526YkzP6hUCtVJIZl0JmXAowtoewpL6CU5VFXCwp4dzFCyiv9iHMWsS6XLGHcCLmbfApPdwpSJJlfK7szVEzNXunppLZIZMZFoLdaee/9r1OU187i8KTeHbN7us25s8VI6NmXjr9AY0DHTyYvfmWq05VKjWSPPO/6ENUCk40N3J3vOe6VmhUKrjBlgaXS2LPpaOcqr+IQqFg06KVbFu2XgTeHKNRq8lJWUROytgbU7vTTkH1Zc5VFXP6fAHHz55BpVDi6+vDisVLWbV0mdhDOA5z7zfeOEiShBLPXuqU3W50aiUNfQPoPdTcdjpJbgVHi8/wcdFxDBpvfnjXc0QHTmwl4Ew7XXuRPcVHiPYL48c7n8dXd/t3vG7UjDicmLQzV2IuS6/lfHOdR4MPQHHNJnanS2JP0WFON15CpVBy9+I1bFmyWgTePOGl8WJt5nLWZo61E7PYLJypvEhhbSmH8k+w//gxVEolQQEB5OXksmxRtthDeAPz8isybLej9vCUT3aDVq1hX9lFUn1mZwX6G3G73dR39yBZevj4Qit3p+WxM3vjTA9rQoZGzbx4ag+tg508uHgLq5OWjvu5oaYQTrf1sC0+YgpHeGvLfY0c6Ozx+HGVbhdOl4P3LhzmbFMJKqWK7UvXszF7lQi8ec7gbeCupWu5a+laYGyh2qnyIi41lvPhwQPsOfAJapWasJAQ1q1YSWaq2EMI8zT4us1WtCrPlhNz40anUtLc3cbDqbPz5rLb7aZ/ZIS67naae3sYso5idzmxOkGj8uGnO/4E0zhmR7PRyZoLfFByjLiACP525/N3/HmsjknkTOPZGQ0+jUqFUnZ49JiDNjPy6AA/fPdf0KhU7MzZxMbs2dUpQ5g+/iZf7l25iXtXXtlDONDL8bLzlDfX8uaHe5DdbrRqDVEREWxYmUdyfPyCfHM0L4Ov32JDq/bsi+l2g1qtRnY5iPDgNonJGLSYqe3qwOZ04D9Uwm8OmwEF3moNAQYTOVGRxAeF84fLzTyalTcnQ2/QOsyLZ/bQNtjNI8u2siI+e0LHWREWzkcVM7/ARS1Lk3q+S3ZRP9BJWU8bPTYzLrcL2TuUB/K2s25RrodGKcwXIf5BPLJ2x2f/bu5u50RpAdXt9bzw9puAG61WS0JMLJvy1hATGTlzg51G8zL4BkZH0WmmpoC0coY2rtscdhp6uqjt7aLfYsHmGruvo1NrUCq02L2jeD49Ej/9l6vV9I86WB42t8qrAZyoLmDv5U9JCIzi7+79Ngav8TeV/SKVSo1rFixwUeFmxHFnATwyaqWkq4mGoR7Mkh2VQkmA1ot10fEkBIbwi8YhEXrCuMSERPCVG+0hbG7m17UvoECBt86L1MQkNq5aPW/3EM7L4BseteGlnprgU3h40cyNSC4XLf291HZ30Gkewep0ILndeKnU+HnpyQyJISkogkC9CYVCwT9dLMDsFXTD0AOQ3EqP14icSoPWIV44vYeO4V4ezbmH3DjPbLWQUWNxShg0M/e10CkUtAwN3fIxsttFU38Pl3ua6baZkdwuvJUqYnz8yI7MwO8LXeEVgM3hwFs7t4oqCDPrRnsIK5pqOFlRSEVVNRfLSq/sT8Tr1gAAIABJREFUIdSTlZbGhhXzZw/h3PlteAfMdjsGT1/qvPK3p1eL4nbTPTz0/7d33+FxVWfix78zd5qmqPdmFctylbsl947BNh0SBwLpJJvefqmb3U22JZtNNrvpZRdCIIFAQgcbMO5VttwkW5IlW713TZ977/z+kDEQN5U7mpF1Ps+T54nlmXuOhXTfOeee932pbm+hub+XAb8fv6pg1OmxmyxkRydRkJhGenQC0ij24v2KiqxOnPZJuyuP8HL5XqYmZfEPW/5uTKu8v5XkSOJwcxfrc8J3mtWu19HsvLJrusvn4UxHA7W9HQzKPiSdjliTmaUZ2RQkpV73g4sB6PX4ROATxsSg1zMnt5A5uYXAUA7hqQtnOVhZRtmp0xwuO45epyfabmP+7CKWL16C3ard7+d4ukkDn5ckW2j+aWNd8Q163NS0t1HX00mP14NPUdABNqOJZFsMCzNSyY5LJkqjwtBHWzuJt8Zrcq1Q6nH28ejB5+l09vCBxbexcMpszcdYmD6F8paysAY+h15Hu8tFEAj4+nmlqoI29wD+S6u6TEcMm9KmE2cffgUdgw563G7SYyZe1R0hchkNBhZNK2LRtKHn6r6Aj6NVpzlSfYp9Rw6z69AB9Do9cbExLJk3n+L5C7FMkA9fN2Xg8wX82EwhCnzB4Qc+fyBw+blcp8uJR5YJMlS+KiHKxqL0fPLjU4iJCt2hk8Z+FxnRkf3Aeue5w7xWsY9pSVP49JptRJlCc3hoenwih+q0Lxs2EtF6HaebynH11pLk6cVl1LMoLZPClPRRFxIw6qDX49V4poLwXmajmZWzF7Ny9tDzZKfbxcHKE5TVlvP6nj28tustJL1EUnw8yxYtYsHsoojt8hGZsxojv+wPSaKyT5aRrlEDNKiqNPd2U93RSutAP65AgEBQxaSXiDFbKUxIZ2pCGsmOOPQ6jcvKXEe/XyY9Vrv2TFrqdvby6IHn6HL38eCSrczPnhHS8VKsVnzy+FfdUYMKNZ0tnOxspM3vp8+WSqYjDpM9h3uzx/4J2aQDp0/bNAlBuBG71cYtC1Zwy4J35RBWHONkfSXP79jBX7cP5RCmpSSzakkxs6ZNj5jUiZsy8PnkAA4Nnw0N0VHf04dZ0kMwSI/Lyfn2Zup7e+j3+fCrCgadHrvRTJo9nqmJaWTGJIa9HJjTL5MWYbX8VFVlZ+Vhtp/dz4yUXD6z7gGijKFPEZEkA8Hg+HzoUFSF6q5mznQ20uPzIOl0pERZKc7IZQfpxJpcuJ3ajGXWDT3XFoRwinPEsLVkPVtL1gPQ1tvFvvKjVDTW8McXnicYDGIyGslKz2BVcQkFOeHLIbwpA5+iyMRaNN5rVhX+vO9ldK4Wfrm/FoAog5Ekm4OilFxy41OwhmiLbizcgQCZ0Ve2IgqXzsEeHjv4HD3uAR4uvoO5WdPHdXw1hIFPURWqOps509lAr9+LpNORarVxx5TZpEQPnYZzyQovtfpAwx9Pkw5cI0yREIRQS41L5P6Vm7n/0p/r21vYW3GU8y0XeezP7+QQ5mdPYe2yFWSlj19xiZsy8MlKgDjL6Lc6g8Eg3QP9VLY00dDdxaDXC26ZgaCdPGscDxQsJC7Kjm4ctyxHy6+oJEaFPyCrqsqb5w6y49xBZqXm87l1H8RiNI/7PLRORlFUhXMdTZzpaqT/UrBLs9pZk1dAov3KDxw2g0RQ1TZIWfQ6BsVWpxDhpqSk81DKXcDQ/aCmpZ59FaVcqG/k5zWPXsohtDB9aj6rS5aRkhi6HMKbMvCpQXVEh1v8skx9RxtVrc10Dg7gCQxV17AaTaQ64iiZMoP/q2gnPzafLKWZ+AnUs04N6sKew9c+0MVjB5+n3+vkIyV3MTtzWtjmosWKT1FkznY2Ut7ZTH/Ag0GnJ91mZ13eNBLtN+7TqENlKPtOGxa9jo6AWPEJE4der2daZi7TMnOBSzmEdUM5hBWV1ZSVX8ohtFqZUziD1cUlxMZot3N1UwY+PUH00rX3jnudg1S3NFHX1UG/x4NfUTBJEjEWGzOTsylIySLB5njPIRTd2U4GvE4SJljn9fF6pnU1qqqyo2I/b1YdZnZaAV9c/xDmMKzy3m20Kz5FkanoaKS8q5mBgAeDTiLDZmd9/jQShhHs3k0XVAHtciujJB1+sdUpTGAGvZ45eYXMyXs7h9DPidpKDlaWcfzUKQ6VHbuUQ2hnwZwili1aPKYcwpsy8L2bIqs0dHdQ1dpEW38fHn8AlSBWg4lkewzzMwrITUjDYrzx1qgn4CElKrw37pEKVwOltv4uHj34V5w+Nx9dejezMgrCNJP3Co5gpRWQZSo6Gjjb1cyA7MWok8iw27mloJC4Maz6tQ58VoMevxzeNA1B0JLRYGJJYRFLCodyCD0+H6XVJzlafZq9hw/x1sH96PV64mJiKJ63gCXzF4woh/CmC3wBWUYXcPLMoX30edz4ZAWjXk+0JYqpCekUJmeSHD26lAJvwEtGVOQcFBmOUB7muOp4qoq3+yw/rDnNnPRpfGDJZswaJeNrQQ3qUBQV6Ro7AgE5wJmOes52teCUfRh1Epl2B5uyZ1xRKmy0RpILOhxWSY8sAp9wE4sym1k1p5hVc4Y6jwy6nRw8d4KyCxXs2L2bV3ftvJRDmHAph3DOdXMIb5rA19rdwdP7XqGxpw3JoMcsxbE8J4e8pDRsGpy2DAaHburWMNZ5HI3xXPG19nfw6IHnkT1ePrPsXqan54/j6MMjSQb6/QHi37Vyl1UFSfHyxOm9OGUfJp1EliOazVmziYnSviSTTuNC5zaDhOwdW9cHQZhIHFY7mxauZNPClQD0DPZxoOIYJy6e47ntr/HX7a/y/W98+5rv13UPqiG9N/7zf/0XLp87lENcwSyBMSaR1mCiZtcc8AXRKX7i9QH8ushZwQA4gwYSLDL+q5yRd3r8+BjfruOGYABZF75O59ej10GyJUDgXd8rvxIE5wA6DQ+cXI/bEIUuKBMdcBLQj33rXEWHVzIiSxNrG14QQumpL953zb8L+fIlMymVT236gKbXHHAN8ud9r1LVehEdOuakTeOeRRuJudRv7uuvP8UP1mpT0f9t39p9EiVo5st5ySRHQHrAu/3gRCnO5EX8y+wra3J+fddp/n3jtX8Axqq5t53HDj6HV/bzweKtFKbm8c2TDfz7vOyQjTkWX331cZbFuWgb7Eev05EZH4fJEM320z38w8zMcZnDdy8OMj3WiNtp5SMaVG4p63HxZjCO791xpwazE4TIpaoqlY217Dt7jLrOZrwBHzqdDqslilmFhawuWUpCbNwNrzOh9u3OXKji+aNv0OsaIMZi4955G1mWP/+K7P9QdcyzGKNo9/ojLvBdT6jaKKmqykundrGn5hgLs2byvsWbMEqRtRJ+N1VVefTAX1HdTcjBaLbMX0h2YjIA+6pqx3UuIzlgMxxuWcFojswVtiCMhaqqXGxvYl95KTVt9bj9HkCHxWymcGo+a0qWkZaSMuLrRnzg8wV8vHDoTY5fqEBWFfITMvm7Fe8jNTb52m8K0YEOu9lKh6cfbvyBImKEIvA19bTx2KHn8Csyn1q1jWkpUzQfQ0sdA138ZOcf8CkBsKZy/9KF4Z2QTtsyTR41iMkgAp9wc2jtbmd3+VHONdXi9LoJEsRsNJGTncXapcvJycwa8xgRG/jq21v48/5XaOvvxGIws7ZgCZtmLxveqkKnI6CoGK+TyzdSOiDGYqfH3anZNceD/hpFtUdDUVVePLmTfbVlLM6ezX2LNoW9FumN7Kw4yIsVe0iLTeILdzzMd154OtxTIqhxxR+PomKJwHJ5gjAcnf3d7K84xpn6avo9g6jBIEaDgYzUdO4p3kJhfr7mNT0j6q4lqyq7Th5gd/lRPAEfadFJfGb1A0xLyRnRdSS9gcG/ObmnhSS7g87+8a/uPxZ6HXhlP5YxphQ0dLfw+0PPI6sKn1n9AfKTI/MZ3tt8AR//89aTNPW3s2XBatbNWwZAuKvMKYqKqtOj5XlbX5AJ0wdNEAZcgxw4V8aJ2gp6XP0oQRWDJJGSlMSmDWuZPX0mhhAXr46IwNcz2MfTe1+htr0BSS+xIHMGd81fP+ru2wbJQK/Hr3HgC5Juc1AbhrY2Y2ExGmgecJIfP7pmtIqq8vyJNzlw4QTFU+Zw38Jbwl4C7Uaq2y7y6/3PYJQMfPPeR0iMSbj8d/qwpfQP6Q7IGAwmQLtKKz51qJGxIEQil9fDkaqTHDt/ms7B3qH0Ib2e+Lh4tizdyKK58zCNc9++sN7BjlWf4eVjuxjwOImLiubBxVtZOGXWmJe1JsnIYED7hN7s6Bhc8sTKl7IbDDQ7B0YV+Oq6mvn9oecJAp9b+wC5iWPfWw8lVVV5qvQ1DtedYlbmVD6y8b4rfpbCveJr9wYwGe1oG/iCpFhEKoMQGXwBP2U15RyqPEFbfxcBRUav0xMTHc26FSspWbAQqyW8W/PjHvg8Pg9/ObCD0w1VqKpKYXIOX1r3EAl27U6MmCVTSAJfYpQFf2jTHjVntxjpcLlG9B5FkflL2RscrjvN0py53LNgQ8Sv8vrdg/zozccY9Ln4yLp7mZNbeMVrnD4/UpgDX5dPxhplBXlAs2v6AYdJBD4hPGRZ5kx9FQfOltHU04pPDqDX6bDbbCxduIgVS4qJtkdWT9Bxu5vVtNTx7P7tdAz2YDNZuG3GCtZNLw7JDdVsMOH0ezS/riQZND+KHmpxZiOtruHfZC92NvL7wy+gQ8fn136QnMSMEM5OGwdrTvDnsu3E22L43gNfIMocddXX1ff2EWXUrkbmaPQFVKIT7dCv3TX9QYidYDVkhYlLVVWqmi6wr6KUi3+TS1c0axZrSpaSEDe6RyvjJaSBT5ZlOvu6+ebj/4lf9pMVk8qX138o5DdTq8mEJzAYkmurhKdj8GjNTYqntLz9hq9TFJlnj7/OkfozLM+bz13zNyCFqTvycAUUmV/tfoqa7gbWzCzm9kudn6+lrLGJJHt4t1gGlCCJdgduDQNfIAixUVcP9oIwVqqqUt/exN6KUs631eP2eQGwmM1Mm5rH2pLlo8qlC6eQBL7Wznb+9MLzdHZ1IeklluXM4/Z5a8et8ajVaMEVCM0hlIm24suPdxCQL173NbUdDTx++EUMej1fXPcQ2Qnj1wl5tBq72/jp7icA+PIdHyUjMfWG76lsbeCB2Qk3fF0oOYMw0xHNBQ2vqQSDxFtFOoOgndbudvZcyqUbvJRLZzKayMnMZO3S5eRmR/ap7hvRLPCpqsq+I4fZefAA/oCf2Cg7D62+k66Obm4rWq3VMMPiMJvp9Gh7COXtQxETLfABGHRXr2UTUGSePbaD0oZyVk1dwO1z10f8Kg/ghRM7eav6CLnJWXxqy4PDPvrs8zuZnRzeZHuPClPitQ18chCizWKrUxi97oFe9peXcqq+in6PEzWoXsqlS+Ou4s1Mz5+qeS5dOI058A0MDvLUSy9wsbEBHTAtNYf3rdxCnGOofc8bHXvGOsSIxZotNIQo7UDr5OPxYNQH6ff5iHnXzbG6vZ4njryISTLw5fUfIjP+xiumcHP53PzXm7+ny9XH+5ZvoXj63BG936BXr9ugeDz4gpChYSdpGPqZvF4LFkH4WwOuQQ5VlnG8toIe5zu5dMmJiWxcv5aiGaHPpQunUf+2nD53lhde34HL4ybKaGHrgjWsKSqJiE8FMWYzPkX7ip2KIqOiR1XViPh3DldWjJ03LtRw34xZBBQ/fy7dwfHGs6wuWMTtRWsnxL/lZEMlvz/yPHazlX/c9jkc1pGdEutzezFHwD9T0ekxax6kJt6HMWF8eXwejlSd4mj1aToHe97JpYuNY/P6DSwumotpEhVBGNFvoNfv5/kdr3Lm3DlUVSUjLoVPbnw/WUmR9Uwo0WrDp2i74tMDMuCIiuFEzwALE2M1vX4oLUiN5+ULDVTFRfHEkZexGEx8ZcOHyYiL/AfSiqry6IG/cLrlPEvyi9i2ZuuorvPquSqyYm0az27kFL32NTUn4i6EEFoB2U9ZTQWHKk/S0tfxrlw6B2uXr2DpwkVhz6ULp2EFvvqmJv788ot09/VilAwsLZjP1uJ1mCO0WkSCxUJA661OnQ5ZUVmQkU9p8/EJFfhmxtt58vgZfr3vNOumFbN5zqoJscrrGOjiv3b+Ab8S4DO3PkB+Rs6or1XRUsf7ZoS3unibx4/OoP3pS1UX3hQNIfxkWeZMQzUHK47T2NOGT/YP5dJZbZQsWMDKJSVEOxzhnmbEuG7g27F7F/uPHSUgy8TbYvj4hvuZNWXaeM1t1KKMJhSN88z1Oh1eWWZDbi7/ev6QthcPkX6Pi6N1VZzvbkHCyJdu+QhpMdfpahFB3qg4yMsVe0iPTeLzdzyMcYy1Rn2+Qealhvck2qEeF1kpszW/rgh8k4+qqlQ3X2Bv+THqOpvwBHzo0GGNsjBn5gxWlywlKT68J5gj2XUD376jR5iRkc/9K24b8TOVcNP69KVBr8MV8JNgtSIT2Tea1v4ejtZV0tTfQ1ZcIvcULeenZXXE2yM7qRSGikv/984naB7oYMuCNaybt1ST60bCwZZaX5BtM6Zrft2gxm2OhMhU397M7vIj1LTW4bqUS2c2myjIz2Pt0mWkp0T+AbVIcd3Al52awYNr74zYLc3rUTXuyWfQ63HJQ2XQFL0Bj6IQJUVQAAwGqWpvpKyxhj6vm8KkdB4uXIfDMlToOz26g9drz3P7tBlhnui1VbVd5DfXKC49Fk6vD5M+/KXm3OiZnpKk6TV7fQH0kujFdzNq7e5g76VcugGv63Iu3ZRLuXR5EzyXLpyuG/gWzpvLzlMH2bxozThNRztKCAKf0+cHINGRwqH2XtalJ2o6xmgoqopPDhDbX8Z+p8rstCnMz5x6RbXzuwqy+O/jZyMy8Kmqyh+PvsLR+jPMzirgwxvu1fQZ5OPHyihIjNbseqMVioMtNU4fVmvkH1ISbqx7oJf9Fcc4XV9Fn3vwci5dekoqdxbfyvSpBRPi2fxEcP3AV1TEb8oeH6+5aEoJ6jVtRmvU6/BcWvGtzCng0Lk9YQ18sqJwvLmGU+0N+BQTOus0vjR/KtI1TvhlOKygal+/dKz63AP86I3HcPrdfGz9fczK0f4Z8oW2Wv55g/ZbjCPR6vGjM46uzdb1NHkCJKWF99COMDqDbicHz5VRVltBt7Pvci5dUkIi965ZydwZs0R+Zojc8LuakJTAucYaZmRNHY/5aMZqsnGh30VhvDYnmQySHveljg9LUtN4+VR42hMFFJljjTWc7mzAYbKwPm82zzY04TQlXDPovS3Vbuatugusy8kbp9le38GaMv5ctoMEeyzfu+9LRIWg+ojT58ekk7Gbw7sdeKjHRXay9gdbOhWYFh/+nQfhxjw+H0erT3K0+hQdA3+bS7eexUXzJlUuXTjdMPDdtnYtT//lhQkX+BJtsZzv6dcs8Jn1El5lKNhJkmHcD7j4lQBHG6op72wmzhLFpvy55MRf2uJqaBrWNe4qyOLXp0+FPfAFFJlf7v4Ttd2NrJ1VzNbi6xeXHovHjhxlVqq2lVJG44IvyAMztd9m7leCzMkQhxoi0VAu3TkOVZbR2teB/1IuXbTDzpply1i2YBFWq/a7AMKN3TDwWa1WFL3KoNs5oU525scl09hbodn1jBKXtzoBkMz0eP3EW0L7Cc0j+znaUM3ZrmYSo+xsmTafrNjRHZDIi3OgyjUaz3Bk6rtb+PnuPxIEvnrnx0lLCG16xcX2Oh65ZWZIxxgON3qmJWu/MvMFYUrcxMkpvZnJqkpFXTX7zx6jsbsVn+xHdymXbsn8+axYUkJsdPifNQvDTGBfvXQZr5Xt4X0rtoR6PppZmJZGWdMxza5nkqTLW50AU+Iz2dPRzN3ZaZqN8W5uv4/DDVVU9bSSbHVwR+EiMjQ45ZhkNXKouZGlGePfTf25E2+yu/ooeSnZfHLzAyGvBej0+bFIClGm8J++DcXBFoCgyOELG1VVqWmuY+/ZUi60N+EJeNGhIyrKwuwZM1hTspSkBJFLF4mGFfgK8nLZuXfvhKpRmWp34NMwi90s6enzvfNcb0NeIc8eO8/dmo0wxOn3cri+kvM9baTaYrl7xmJSHdrl391ekMXjFcfHNfC5fG5+/Mbv6Xb38f7lW1gywuLSo/W7Q4eZmx7+1VCT04feFJotLZG8Pr7q21vYW36E6tY6XD4PEMRsMlOQl8OapcvJSA3NB2FBW8M+MpSbM4Vj58+wpHB8blpa0DKXzyRJ+JR3At+0hAScGlZFG/C6OVRfRW1fO5mOOO6bVUKSXfub9ozEmBv259PSiYZzPH7kBRyjLC49Fo2dDXx206xxG+9aDvd5mJKWH5JrqyJ5PaTaejvYe6aUs401DPhcBINBTEYjUzIyWV2ylLzsKRNmMSC8Y9iBb/2KFTz25FMTKvBpmctnliRk5b0nORXd2I8a93ucHGyo5GJfF9nRCWybs5x4a2hr6iVYDRxvaWZhekbIxlBUlf/b/yxnWmsoKZjL+1aN7zZ5n9uLxaBgNob/pnTep/KpOXM0v66iqGKrU2M9g33srzjGqbrKv8mlS+H2JZuYWTBNBLqbwLDv3AaDAVOUibbeLlLjJsbxaSWoR1FUJA1y+UwGPYG/CXxmczSnegaYGz/yB9Y97kEO1VdRP9BFTmwiDxStIDZqfFZDH5iZw8/KDrIw/f6QXF/xO/n7F35CQJH57OYPkpc2/hUmfnf4IAvSw1+iTVFU3HoTWXHanyxt9PgxmsWpwLFwul0cOFdG2YUKugf7UIIKkl4iKSGBe1atZN4skUt3MxrRf9GNa1ez4629fGjdPaGaj6bMBjONTjc5MWMPKFEGCVn1vedr2+YU89zx10YU+LpcAxxsOEfTQC95scl8cO5Koi3j2y4nJ8aOWeejrr+XnBhtk59frziAs24/OXGJfO72D2EM002jtauBLy8J/+7Emx2DJCaEpuv7RZePaJtIZRgJj89HafVJjlafpn2g+3IuXVxMLLeuXcviefOxiFy6m96I7kppSSn0e50EZDlsN7SRSLDFUdUzoEngs0gSshp4z9dmJiXzqKIb1qGf9sE+DjVU0uLsoyA+hYfnr8Zu0r5FzXBtm5XDr0t38+8btDmeM1Rc+g80D3SiJi/ly1vXaHLd0fjdocNMS7RpVrVnLI65Fb60aXlIrn3Rq1KYlxmSa98sArLMidqzHKwso7W3A78SuJxLt3rpMpYvFLl0k9GIo1fRrJnsLj/Cxnmh+WXWUm5cEvX91ZpcK8pgQFav7OqeEJPGzrZuNqZfPbeuZaCbQw1VdLgHKYxP5ZaC+VhN2lcoGak5SXHo1It0u90kjPEXv6q1lt8c+Asmg5Fv3ftJ/vXMgEazHJ1zjVX8cHNRWOcAEFBUvJKZ5BD1QWuXgzw8bWIVlgg1WVU5V3eefWdLaXhPLp2VRfPmsmpJCbEx4S9oIITXiANf8YIF/PbxJyZE4JuXnMqf2k5pcq0oo4SqXnmM82MLlvLTPc9eEfia+rs43FBNp2eQmQnpbC5cSJQx/AHv3bZMy+R/ju7iu2tGd/BEVVX+eORljjaUU5Q9jYfX33Np5Ru+wPf7I8fIS4iKiNXeq+0DpKcUhOz6AZ1ErHXydtGGS7l0LfXsqyiltr0Rb8AL6IiyWJg1fTprly4TuXTCFUYc+PR6PbFxMdS2NpAfhkMLI5EbF487oE3OgdUgIQevvFaS1Y4n+E6boobeDg41VtPrdTM7KYOtMxZhGWMT1VBZmZnMi1VluPx+bCN8rvF2cWmX38PHNtwXMQ2KT9dX8INbw7/aAzjtCfLNLaH7gBiqpPhI19DRwp7yo1S3XHxPLl1e7hTWLV1OZlp6uKcoRLhRPai7bd06nn/x1YgPfABqUJtP/g6T8aorPoCClKn8sboCq7+DAb+HOUmZ3D2rGNME6JO2OieFn5bu4RvLNw77PQfOl/HMiR0k2OP42n2fDElx6dH44/ET5MRZIiKFwSsr+CVLyFZkXd4AOkNkfN9Dra23i33lR6loPD/Ul+5SLl12Riari0vIn5IjUgyEERlV4It2OPCpAVxeDzZL+A5oDIdWuXxGSQ/BKyvB1HU10dl2ikF3J3dnpbIwYypGKfIP/rxta34Gu+rKUBQZ6QbzDigyv9j1Ry70NLFuzlK2LF47TrMcnmO1p/n+Ju3z5UbjhbYB8jK078bwtuO9LpLjQ5MUH269g/0cqDjGyfpKel0Dl3Pp0pIv5dJNE7l0wtiM+g69bPFiXj+xj7uX3qLlfLSnN9Lh9pKsxSdv3TuBb9Dr5OnS1zjbdoElU2ZzvNvBtOTUCRX03rYwLZ6fHjvIF4tXXfM19d3N/Gz3n9Dpxqe49Eg9c+I0WTHmiKjLCVDlg+8uKw7Z9Wt8KoumREZ7qbFyul0cqjzBsdryob506qVcuvgE7lpxKwtmzxG5dIKmRv3TNGt6IfsPH9FyLiExJTaDfY0d3Fs49m1ZPUMVSXZU7OOtqqNkxaXwjVs/QbIjHkvVWf5Ud4LPTA9NzlYovb8wm6/tKqPduYAU+5WpH8+deIPd1aXkp2TzyDgUlx6NQzUn+NeN4S9PBjAYUAgYoogKYT5YrwIr8ibezxoMpb4crTrNkepTtPd3IasKer2euJgYNq1ew5L5C0QunRBSY/oYlZGRysnas8zLD3/bl2vZmFvIM+U7QIPAF5TdfPflnwPwcPHtFGW909X7rsKZfL3m6JjHCAdJ0vOhonx+fGg7P9h43+Wvu3xufvTGY/S4+9m2YiuLCyPj0Mjf+vm+/UyNiwp7s9m3Pd86wKycBSEdQ9EbMU+QVVBAljl14SwHKsto6XlvLt2qkqUsW7QYu8ilE8bRmH5zNqxew5NPPRPRga8gIR6X/8r8u+Fq7e/hXFs9db2dBL0+SgqWsmnWSqSrrHokk4PzA04KoidO38K3FSXHkVTbzHMYkZ/gAAAeYUlEQVRVFdxdOIsT9Wd5/OiLOMw2/mnb57Fbx7e6zHA19vZT31bND2+bH+6pXHYhAP+yZGFIx1D0kRv0ZFXlXP159p87RkNnC95LuXQ2q5WFc+eyuljk0gnhNabfHovJhN4k0T3QS0K0tqWvtBQYwclONRikqa+LqrYG6ns7kYMqWTEJrMqfTZW7k81zVl/zvXfNKuEv5Tv5xuyJF/gAvrhkJl/beYz21nLOtteGpbj0SP1k58s8sngq+gjI2wPo8PpRTLaQrsbqnF6MIWpzNBqqqlLb2sC+iqPUtjfi8b+TSzdz+jRWlywjJXF0zZMFIRTG/Nu5YdUqth/cx4Or79BiPiFhkCzU97uYEnP1VYsSDNLQ3U5lRyONvV0AZMclsn7aPLITUpB0QydDdbrO646zJD2DZ05oO/fxNOB1kmXq5GyXyuc2f5DcCE9X+cW+/eTEGpmZEjmrh8dbnNy9fHNIxyjr95CVGt6KLY2dLew5U0p1ywWcl3LpTCYT+dlTWLtsBVnpIpdOiFxjDnxZ6em86uxFVtWIPPQAkJ+YycGmdqbE5F7+mqKq1Ha1cr6jicb+bgx6iSmxSWyeuYj02ET0uivTIIaTGBHrSOH1lk5uuUYJs0ikBoMcqC3nRMtF8uLj6MZKk1sh98ZvDZuW/kHqImyLs8Prp1+KYmVeTkjHqfMFuWtq6CrCXE1Hbxd7K0qpaDjPgNeJeimXLis9g/uKSyjIyRUpBsKEocl+zPTCaRw8e4xVs5docTnNbcmfzm9Kz6OoKtUdTZzvaKZ5oAezwUhOXDJ3zi4hLebGLWx0uht3dP9iyVq++/oTEybw9btd/PX0AdwBP7fOmUd+ajor3T7+cecelk7Jithj5D964yUeWZwfMVucMLTauyfEqz0AV1BHUXpKSMfocw2wv/w4J+vOXs6lM0gG0lKS2bJkA7OmTReBTpiwNLmrrSwu5n//8GREBj5VVWnursc12MGvDryK3WQhJz6ZktyZJDtGtkU2nBWfzWQiMS6bvza0c092aG9OY3WsvppD9VUk2+1sW74K06UgF2c1U5wVx3/s3Mm3Nm0a1znJssygx8WA28mgx0m/y0nPYB997gGcHjcun4d2l4skmx5JjsYv2y/PO5zaPOOz2oPQlCpzeS7l0tWU0/WeXLp47lyxifmziyLi+ywIWtDkJ1mv12OLttHQ0UJ2cmTs7Q96nbx59hClDRUEgyqKwcIHFpQQbxtLpfwbr/gAvlCyhm+/9jh3ZSZF5KdiT8DHX08eoNfjYk3hDGZm51zxmgfn5fCt7Sd57VwVt80oHPa1A7JMMODmYlsjg24XA65Bepx99LkHcXrcuH0ePAEffjmArMjIqkowqKIGgwTf9f3VATqdDtCh1+uRJP1QKyy9kaBOJcvu561zZ1AqTmOSJOKsNvJTUpmWnoXdMv6Fm//QOsg9y7eGfByPrKBqUArPF/BRWn2GI1WnaOvvQlZl9Do9cbEx3LJqNcXz5mMJw/dREMaDZh/hblu7jldefYOPbQxNV+/hqmq7yJvnDlHb1UBadBJ3zV3Hoimz+cc9r9DuUYkfw6l8SQde2X/DotNmyUBuyjSevNjMQ/kZox8wBCpa6thVc4aYKAsfXrXuunU2v7NuDt96bTftrefRKT6cXhdunxdvwI9fCRBQZBRVGQpa7wpcChI/vaCi0w0FL71ej6SXMBgNmE0mbNFWUqx2YmOiiYuNJcYRTYzDQawjGrvDcd1nxV/62WN8fG06C7KHtqbdHi/VDc2cb2yl9GINB2uqkXR6HFEWchKSmJGRTUL08BsFj0abx8+AZGXlOCSU7+xwkpo48oMtAVnm9IVzHKgso7mn/XIuncNuY0VxMSsWLcJum5inkQVhpDQLfHGxsbhlL76AD/M4t9/xyX72VJVy6MJJBn1u5qQX8LVbPkpqzDtltWamZHOopY4ZiaM/AWgySHS43GTH3LiqxKcWlPC1155gm6pijIBVn1+WeeH0QdqcfSzJzWfR1Guv4nqdgxysPktjTy85FoUzjQMYdWAyGTGbzdhsNlJtduKiY4iLjSXWEU20w0FsdDR2m43PvlTHL+7UvpzW13/7NMU5tstBD8AaZWFeYT7zCofqVgZkhbqWNs5dbORcaxOnmurRocNqMpIeG8f09GwyE5OuenhptMZrtQdw1qPw8Jp5N3ydrKpUNdSw7+wx6jub38mli7KyYG4Rq5eUEBcbOw4zFoTIo+mm/eL583jjxH62Llmv5WWvqbm3ndfPHqS89TyxFgcrpi5g1bSFGKUrA9NteQX85ED5mMYzSxLdXg/ZMTe+YUiSgaLM2fxvTTWfmhbetIDarha2nyvDYjTwweWrib5KMrrb6+NITSU1HW34ZBmb2cjCGdnMnZrPrqoOnjvdy48//dCwxluVo33j1e89+QpZ9gAfWjb9uq8zGiQKsjMoyB5aaavBIK2d3Zyta6C5s4eajg7UYBCLwUCyI5rC9AxUefQFDsZztQfg0RmYmnRlf7l3culKqW1veFcunZkZhdNYU7KMlKSJceBKEEJN08A3f84cfn38cS0veQVFVTly4RT7ao7TPtjFtOQcPr1qG/nJ1w8usRYrPnVsn/ItBolej2fYr//QvEV85eVyPAGZKOP4HwxQVJVXyo9Q19dJUUYWK2bMvvTcbIgsy5RdrKG8uRG334/ZIDEtO4XFMwsxm955jrRhZiq1nW7+/rHn+JcP333DcbfN1fYG+8tXD6D3dvCVu0feeUGv05GRnEhGcuLlr/UODHLuYgN1rV28da6cfo9KVMDJM6eayYqOY2pSKgk2x3u+V9fyh1Yn9ywfnyR/r6wgS+/spjR3trHnzBEqWy7i9LkJEsRsMpGXPYU1JcuYkpk5LvMShIlG87txckoSFfXVmjcm7XP383rFQcoazyHpJZbkzOazax/AZh5+BQs5qEdRVKRRHoGPMhjo93lH9J4VUxfxq/Mn+dLMnFGNOVrNfV28VHEUvU7H+4tXkHjpOZcaDHKusY6yugv0e70YJT3ZyXEsnTsTh/XaLaY+uTqP775Yzs9f2ctntly7i4PWXjpaQVP9Ob5/n3bthuKiHSybO4tlc4f+/PqxSnb2XMSAm3M97ZzpbgfApNMTb4kiLz6J3IRkrKb3buGHarUXDAaRVQVZUfArCqqi4A8qHOhyYTHE890//ZQBz+BQLp3BSGZ6OvcVL6UgV+TSCcJwaB74Nq1dy9PP/FWzwFfRfJ43Kw9T19NMdmwq7190G3MzC0f1C263xHCis5dFqVduFQ1HlFHPoM83ovfcPX02X6k9QX8gQIwx9EWUlWCQnZVlVHY2My05hQ1FQ8WS6zraOVJTSbfLhQ5IibezoXgOyQnDLzX37S0z+MYzp9l5Mo3180KfQH3yYgv7jh7k3+6Zg0EKXbsho0FC1pvYMn8oEgaDKj2DTs63ttPcP8jh5jr2N9Uh6cAqGUi1RTM1MYmneiW2LtlIa3cHTq8bl9eN0+OizzXAoNuF0+fB4/PiDXjxyQH8cgBFVZBVBQI+Ej29/K7Df/ksa/DS/97+wtvrzbcXnl3mOHzmAIkpydy6YT2zp0+P2KIRghDJNA98VosFVYJ+1yAxo0wd8AS87Ko8wuGLp/EGfMzLnM6DS7aQ6Lhxkvn1zEvLo7Tl7KgDX7TZyAWne8Tvu23GUn5ZdZBvzA5t/7RO5wDPnz6ArKrcvbAYgyTxYukhWvr7CRIkzm5mU3EROaNMfjZIEn9/xyy+99xbIQ98rb0D/OGlV/jGlkJslvHdJtbp9CRER7/nNKiiyDR2dVPT3knTYA/nB3pxW5N4+cAzvAzo0F0+xarT6ZH0eiTp0klWowmbzUZylBWHzU5sdAynTpYSlOzckqFikCSMBgmDZMBokJAkCUl3ZUD7p7MefvjVz4zjd0IQbk4huaOsWbaM7WV7eP/KkZ10q+9u4fWKA5xrv0CiLY6N05eyLH/eDTuDD9em3Fz+pb501O+PtZhwdY9sqxNgQ95UtlcdocPjJTlK+9yoIGDytPCnsgukx8RiNZl4+UQpflXBYTGxvGgqM/KmaHKSMdZqYvrYPn/ckFeW+Y8nnuHjK6eQGRcZxZglyUBOSgo5KUMfGv69YpD1mzazrmh0HwBqq8txByRS44e3kvX9zfM9QRBGLySBLz8nhzd270FVb3xaTlFkDtSeYH9NGV2uPmam5vGFdQ8xJUH7RHiLwYRP0Y36OV+ixYI70DOqsd83dw2/Pv0G3ynStriw2+9lwOcm0dOIUa+nua+XKJPEnKkZzC8swGiIjI7kI/HM9jdYXxjLgikhjrCjtL1hEL89cdRBbzTeanYTk3hzdFwXhHAL2R5Sfl4OR6pOXfPvuwZ72HH2AKeaqogyminOLWLt9GKijKGtFuGISuBIazfLMkd+8jDFZsEvj+wZ39uWpGfw7GkLlX2DTI8d+3H/YDDIgbpKytovoqpR+HUmZmXYKJkzHet1ktInApvez8Is7VMitNDvk9k/aODfv/T+cR233BnkrltD2+NPECaLkAW+tcuW8+iTT1EYm3X5a6qqcqqxkl3VR2nsayM3PpMPldzJrIzx++R8Z+F8Xq3aNarAF2MyIivyqMf+8oot/OeuP/P9+QWjTmrv97g41FBJTW8HCkHsBh1GVaIjdg7rFkVmsLiZ/KLGw6Z1m7CMc91Kt87AvNzIKAcoCBNdyH57DQYDFqsZl8+Ly+fmzXOHKK0rR1YVFmbN5KPL7yHWGtpSUlczNyWFP54aXfCSJP2wOjRcS6rdwZL8xfzo7MkRHXRRVJXTrRc53d7IgN/DUMgMsmVuIdkpyfzrWydHPSdh+J6/6ESJTePWhddPoteaT1YIiOd7gqCZkH5svWXtGn73+BMcrDtJqiOBrXPWsDh3DlK4j2DrLLQ6PaTZr523di1jnfn7Zhbx7dZaXm/p4pb0xOu+tmWgmyMN52lx9WGWJNJs0Tj9HqKMOu5bsQSL8cal0wRtdLkDHPMY+eEn7xn3sfe2unAkRHJ3REGYWEIa+FISkwgS5IvrHiI7BIdVRmt+xjReqrnAI6M4kq8bZoeG6/mHNbfzzdeeYFF8NPGW9wYvT8DHsaYaqrrb8CkB0mwO7po9j5qONk61t5CfGM3GBXPHPAdhZH51wcvdm28PS3/CUwOwdZ14vicIWgn5b7HdZqemoyGiAt/dhTP53lunR/XesWx1vs0sGXho8Sb+89h2/m1+AcFgkPNdzRxvuUi3x4nDZGZuegZzs3IIqirPnDhCn9fDhln5TM2InO/jZPF07SBSQiarQpyHeS1OjCwpEOXHBEErIQ980VY71Z11rJtREuqhhs0sGfCrEgFFxTjCtAatavrPiI/HFhXHD4+XYld60QHZMbHcNnMODutQ7lpHfx/PnzmOpIcHVy7CFjXyrVlhbNqcPk77zPzn390ZnvFdPhST+O8uCFoKeeDT6XQkxMRR29Fww0LS4yklOpWd9W3cmjeyFdRotjo7B3uoar1AVXsdrf2d9Ptc+GQ/ZoMRf1DP1owprM/Peac2FVB64TylzfVkxFjZsmQ+uqtU8hBC77d1fj54z71h2eIEeLHRS1HRorCMLQg3q3H5bb5lwQr+um9HRAW+989ewKPHXhl54LvOVmdAkanvauZsay0Xu5rocvXh9LsJBoNEGc3ERDnISclkZk4BhZm5mI1mOgYH+Y+XnmJlbhCTpENRVZ4/eZR2l5Pl07KZnTM+7W6EKz1+fpCo1FwWTQ3fNmOrYuRza+aHbXxBuBmNS+CzRdmQgwoun3tE3RRCKcsRgysw8tWbpAOfIuMNeDnfepFzbRdp7m2n1zuAJ+DDKBmwmaJIdMRSXDiXWVMKyEhIvWZR7WSHg+IZS/jByeN8bnoiz548SpAg7182n1i76IgdLnX9XqpkCz95+I6wzaHfGyBgtIRttSkIN6tx+41aMXMxb547zJ3z1o3XkDdkkGyc7xmkIP7qid9+WabHNUCHs5/OwT56PS5UzyDf+MuPUIMKFqMZh8VGWlwSq7KWMCM7H4d15MHq/vlz+X9Vx/ntsVJyHEbuXLoo/Ckfk5hXlnm0UeaRB+4L6zyea/RQMFWc5hQErYU88CmqwqDbSWpCIm+c2E+/x4leqxMiY7Q0s4CXzp/mzvwkOgb66HEP0O9x4w748SkBZFXFJEmYJANWk4noqCiiohzMmjqPu+Yt1OSTuCzL/PTlPyC7WqkzJLJtTpYIemGkKAo/rPSwdPkaZmYmh3UuDX6J7966NKxzEISbUcgDX3NXG//wp/++/OfvvPg/aNAkQDMWA+yUm7EYjdjMJlLjYkiOiSE5Oo44mx3935z6rDt8nqDBpknQa+3u4L9f/j2yqvCZhz9MuzvIj195iR+szkIKYf854dp+VOkie+YC7l0e3lxJd0DGbzCPe2k0QZgMQv5b9Y9f/url/+/3+/nDn57h05sfDPWww/aVZ5/kgRUziDINL9BEWwx0u5xjHveNE/t5rWwvCXGxfPFjj2AwGMgCTk6fy3f2neKfV4rgN95+VjmAMX0qn946fh3mr+WFBjeZWdp1nRcE4R3juqdmMpmQzAa6+0fX2icUshOzef5c47BfH2sx0e9xjXo8WZb50XP/y6tle1hZUsJXP/np96weP75pGfHZM/jugeZRjyGMnKyDXlsaf//gyHpIhkqNV+ITty0P9zQE4aY07g+TNq5Zxatle8d72Gt6aPECTrf2Dfv1KXYLLt/Iu7DD0Lbvt5/4Ma19nXzuwx9l85qrH/T5/J1rMCXn8r39DaMaRxiZUw19BCSJ739yfFsNXUtAUfBJZqJtoW3RJQiT1bgHvozUNHo8A8jy6Nv7aCneZsMj6/AFbtw0FyAr1obX5xnxONtL9/CjF/6XmNhovvvlr5KRmnbd13/t/lsIRGfy/YNNIx5LGL7dlW10dPbhNMSFeyqXvdLgJDFVlCgThFAJy/HBWdML2VtRGo6hryo3JY+nztQN67WJUSYUdfhBOyDL/PCvv2XH6f2sW76Cr3ziU8M+GPOdB7fQZ07kx4dF8AuF8uY+dh2tZEP0IDrNitGNXblbzyc2h/85oyDcrMIS+JYvXkx50/lwDH1Vn1hWQnnbwLBeq5f06IdZtqyhs4VvP/EjOga6+cJHPs7GlatHPLfvffhuWnSx/Ly0ZcTvFa6tudfNU2+d5uMZXgyRE/NQFAWvZCE5VhQvEIRQCUvg0+v1RMc6qGuPjJWM2WBAwcKF7sFhvX44HRpeKd3FT158lIT4eP7pK18jLSVl1PP7t4/fT03Axu/KWkd9DeEdA24/v3yljAeSXcSZIyjqAW80u4lOTA33NAThpha2JKHb1q7jpVd28PGN7wvXFN7jjvnFPHXqAN9aN/uGr73eii8g+/mv5x+jtb+TjStXs375Ck3m94NPbuOrv3iCJ8+08+Cc9wZRNRhE9QxSfnHsaRbD4fT4Qj5GS3c/ijK8564joahBHtt1jnlSL30ulT4XNLsUzccZreOD8Ml7xTanIIRS2AJfbEwMXsWHx+cjymwO1zQuW5mXwwvH96Aq6hVJ639Lf40VX317Cz9/9Q8EdfCljz1CSlKSpnP8/iPb+Novn8RyrpN7Z7xzbavJwoleiZ+c0HS4a9MX8tSpTrbN1fbf97bieYt4cvcxGs95Q3J9WY7hVTka3hW/u2z5IRlrJPp9Mh6DlSlJkXPQRhBuRmEtC7Fk4UJeP7mPO4s3hHMal8U4Uthe08rmwozrvu5qYfHFw2+yq+II6ckpfPrDH8UQgrJjBoOB//i7B/naLx4nytjN5qkJANijzMTHz+Dnd4anUarWMtIz+NoD1/9vcDN6qs5N0exl4Z6GINz0wloUcu7MmdR3Rc6hjb9buZK9Fztv+Dq9Hjx+PwC+gJ/vP/MrdlUc4bY16/j8Rz8ekqD3NoPBwL99+mHebPbz69JaflNai9Ov/ZagMP6aFSMf2bgk3NMQhJte2AsBpqYmc+ZiJXNyp4d7KiTabHgVA71uH3HWa2+/RhkkLvT0ERX08svtf0Sn1/GVT3yKpISEcZmnxWDg+5/9KE+d6mRv3dCBHIdJFLaeyA62OjE6Rn8AShCE4Qt74Nu0Zg1//PNfIiLwASzOm80TJ6v43LLCa74mOsrIa8f30NxWTWZqGp9++MPX7LcXStvmJoXsOZswvnb1BHnk4Y3hnoYgTAphXyZYLBYw6OhzDS+PLtS2LZhLfe/VK7MEg0Gqmhvx9jfQ2Haeres38tkPfzQsQU+4eQwdaokiLyU+3FMRhEkhIu7Y61at5NXje8I9jcv0xmiON3Vf/rPX5+fN02X8eud2dp49Q5TFiD1jNisWi+cxwtgNHWqZF+5pCMKkEfatToCczCy2v/kWqqpGxOrpQyUrePrwdlIsCnsrz9LjdmE2SCwozGbe9ALK6nqoPTU+OXPCza9ZMfJlcahFEMZNRAQ+gIKp+RyqLGP5zEVhnUdA9nO2pozBvlaeL2sgMdrGvasXk5zwTm5VdoINv+/Gpz8F4UbEoRZBGH8RE/hWL1vG/z3+RNgCX2NnC3/e/xrNPe3o9TqCtmRMiTbuXzX1itcm2IwE1cip9iFMXOJQiyCMv4gJfAa9HqvdSnN3GxkJ41OrUFVVdp8+ws4zB3H7vdiirGy78y7mzpiJV5b57m8eu/pcJema1VsEYbhkNSgOtQhCGERM4APYtG4tr7++i4+svy+k4/S5Bnh678tUt9YTJEhuZjbb7riTaIfj8mssBgN+nYWajkGmJjuuuMZwOzQIwrW0BwIULQzv1r4gTEYRFfiS4hMY9LkJyH6MBpPm1z9RU8GLR3fS5xnEZDRyy6rVrCouueaBmttXLefxA3v43t1XFq7WR1ZRf2ECcukNfGRjcbinIQiTTkQFPoCFc4t489Qhbls48t51V+Px+Xj+0A7KLp5FURWSEhL5wrb3D6tN0JpZuby0ew+dA16Soi3v+Tux1SmMhS+ox20Zn0o/giC8V8QFvgVFRfz25BPcNsbr1LU18vT+V2nv70LSSywoKmLrhlswDbP7+dvWLF7ML3cf4x/umPWer0sE8coylhFeTxAAPvexR3jqlDgZLAjhEHF3bb1eT0JiHFVNFyjMHFm3AVmWefPEAfacPYpX9uOw2/nQ/duYnj/6ljO3L5nF7qNHGHD7iba+s/1qNUtUNnYwLzd91NcWJjdRbk4QwiPiAh/ArevW8+xfXhh24Ovu7+FPe1/hQmcjOqAgN4/7t96B3WrVZD6Liubz813lfHPLjMtfi7WZONfUKQKfIAjCBBORgc9utSLrVJxuF3ar7aqvUVWV0upTvHJ8D4NeFxaTmTs3bKJk4ULN5/OBVfP5ys9O4PLL2ExD37Jkh4nazu4bvFMQBEGINBEZ+ABWLi3htRN7uH/55vd83eVx8ezBHZypr0YNqqSlpPDIww+RFB/agwLTC2byy101fHXTUBeJjLgoTrRHRmFtQRAEYfgiNvAV5ueza+/+y/U7q5pqefbADrqcvRgMBpYvXswta9aGtOnru31iUwlf/Vk5PkXBLElMS3HgOdk0LmMLgiAI2onYwAeQMyWL321/mtqORvxKgLjoGB554CHysrPDMp/M7AJ+u/sCn11fQGqMlaASCMs8BEEQhNGL6MC3btUqfviLn1E4dSr33rZ5qHdfGH1283K+8avHkBVFlC0TBEGYoCI68JkMBr79+S+GexqXGQwGkpKz+f3Bej62Mg9JlC0TBEGYcMLf/G6C+co96zjVOHSoRSdWfIIgCBOOCHwjZDAYsMel8scj9UiiXqcgCMKEIwLfKHz9/bdypLYXgx66BkQndkEQhIlEBL5RsBgMmOyJKAEfp+rawz0dQRAEYQRE4Bulr2/bjNftprZVFBoWBEGYSETgGyWHxUTQlozP5wn3VARBEIQR0HUPquJooiAIgjBpiBWfIAiCMKmIwCcIgiBMKiLwCYIgCJOKCHyCIAjCpCICnyAIgjCpiMAnCIIgTCoi8AmCIAiTigh8giAIwqQiAp8gCIIwqYjAJwiCIEwqIvAJgiAIk4oIfIIgCMKkIgKfIAiCMKmIwCcIgiBMKiLwCYIgCJOKCHyCIAjCpCICnyAIgjCpiMAnCIIgTCoi8AmCIAiTigh8giAIwqQiAp8gCIIwqYjAJwiCIEwqIvAJgiAIk4oIfIIgCMKkIgKfIAiCMKmIwCcIgiBMKiLwCYIgCJOKCHyCIAjCpCICnyAIgjCpiMAnCIIgTCoi8AmCIAiTigh8giAIwqQiAp8gCIIwqYjAJwiCIEwqIvAJgiAIk4oIfIIgCMKkIgKfIAiCMKmIwCcIgiBMKiLwCYIgCJOKCHyCIAjCpCICnyAIgjCpiMAnCIIgTCoi8AmCIAiTigh8giAIwqTy/wG3F8fjLGX+0AAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "isdark = False\n", "from rayoptics.environment import *\n", "from rayoptics.elem.elements import Element\n", "from rayoptics.raytr.trace import apply_paraxial_vignetting\n", "opm = OpticalModel()\n", "sm = opm.seq_model\n", "osp = opm.optical_spec\n", "pm = opm.parax_model\n", "osp.pupil = PupilSpec(osp, key=['image', 'f/#'], value=2)\n", "osp.field_of_view = FieldSpec(osp, key=['object', 'angle'], flds=[0., 22.5])\n", "osp.spectral_region = WvlSpec([(486.1327, 0.5), (587.5618, 1.0), (656.2725, 0.5)], ref_wl=1)\n", "opm.system_spec.title = 'US 4'\n", "opm.system_spec.dimensions = 'MM'\n", "opm.radius_mode = True\n", "sm.gaps[0].thi=1e10\n", "sm.add_surface([42.71,3.99,1.7343,28.19])\n", "sm.ifcs[-1].max_aperture = 14.47\n", "sm.add_surface([195.38,0.2])\n", "sm.ifcs[-1].max_aperture = 13.53\n", "sm.add_surface([20.5,7.18,1.67133,41.64])\n", "sm.ifcs[-1].max_aperture = 12.01\n", "sm.add_surface([0,1.29,1.7919,25.55])\n", "sm.ifcs[-1].max_aperture = 10.745\n", "sm.add_surface([14.94,5.35])\n", "sm.ifcs[-1].max_aperture = 9.195\n", "sm.add_surface([0,7.61])\n", "sm.set_stop()\n", "sm.ifcs[-1].max_aperture = 9.0295\n", "sm.add_surface([-14.94,1,1.65222,33.6])\n", "sm.ifcs[-1].max_aperture = 8.75\n", "sm.add_surface([0,5.22,1.79227,47.15])\n", "sm.ifcs[-1].max_aperture = 9.635\n", "sm.add_surface([-20.5,0.2])\n", "sm.ifcs[-1].max_aperture = 10.19\n", "sm.add_surface([0,3.69,1.79227,47.15])\n", "sm.ifcs[-1].max_aperture = 11.48\n", "sm.add_surface([-42.71,37.32])\n", "sm.ifcs[-1].max_aperture = 11.985\n", "sm.list_surfaces()\n", "sm.list_gaps()\n", "opm.update_model()\n", "sm.list_model()\n", "pm.first_order_data()\n", "apply_paraxial_vignetting(opm)\n", "layout_plt = plt.figure(FigureClass=InteractiveLayout, opt_model=opm, do_draw_rays=True, do_paraxial_layout=False,\n", " is_dark=isdark).plot()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAKYAAAEYCAYAAADWP43GAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO2da3hTVdr3/zttaVKaQxsQyqktJJVGaDkWZQo4OIJgoaKIHMZHUBCvPqKvfV7nmdFxZt5BveYa5+o1Dsp4QGCG0wgollZAVOSkQNAWCqTQtKXlUE6mNE1p0tJmvx/SpDvJTrKTnTQ76fr5oTZZa+17rf5Zh3utdW/KYLLSIBAEhijcBhAIbBBhEgQJESZBkBBhEgQJESZBkBBhEgRJbDAKuVBbg5Jv9oO20piYPQYPPjDZ6fvaS5dQ+u1+XL95E4vy52H0yEzHdz+dqcCBH44CAKZPzsX40VnBMIkQ4fAWptVqRfH+fXhu4WLIpTK8t3E9MtVqDOjX35FGIZPhyUfn4PCJE055W81mfHv0CF5c+ixAAe9tsOVNEEv4mkWIcHgP5ZevNUCZlAylIgmxMTHI1mig01c5pUlWKJByzwBQFOX0edXFWqjS05EgkSBBLIEqPR1VtbV8TSJEAbx7zGaTCXKp1PG7XCrD5YarfuSVMfJK0WwysaY9caoM2lPlAICnH3sCffr04WG18GhuMUGWKPWdMMKwWq1ITEz0Ox9vYbLvZ1Ksn7rn5b4bOmnMOEwaMw4A0NraCqvVyjlvJPD39evw+xdfDrcZQUckCmxQ5j2Uy6VSGBm9nNHUDJmU278QuVQGo6mZkdcEmTT6eg2C//AW5pCUQTA0NqKxqQkdnZ04rdNBo8rglDcjfTj0F2vRajGj1WKG/mItMtKH8zWJEAXwHspjRCLMnTET6z/dBittxYSsbAzo3x/7Dx/CkJQUaNQZuHytAZs+3wmzxYLz1Xp8ffQwCpevRIJEgumTc/H+xg0AgId+MQUJErIiJwBUJB57i8Y55pvvvRu1c8yEhAT/84XAFgKBN0SYBEFChEkQJESYDDrWFuHunGnoWFsUblN6BCHXlwiTAb23GLBabT97AUKuLxEmA2pWPiAS2X72AoRc36Ace4sWYgsKoc29An39aqgrriAnS3hDXLDQVhRCP3QD1G8tQ05WYbjNcYP0mC7o6zeApjuhr98QblNCitDrSYTpgjp1GSgqBurUZeE2JaQIvZ69bijvWFsEem8xqFn5iC2wDWFrtAaU6k3IU0uxKqcoqodwOzlZ3fV0rr8SAHs79SRR0WP64/ZgW4mW6k2w0rafvRG2+vuzYg+F2ykqhOlPI7KtRPPUUogo289w0LG2CK/u+yxs/kS2+vuzYg+F2ykqDnH4O+xoKwqhr98AdeoyQQzbd+dMA6xWQCRCXMmhcJvjd/t4a/9efYijLPcKtj9fgbLcK5zSC21FSs3Kh5WiBONP9Ld9/G1/LkSkMM/fMmON1uD43VdDrtEaMHNLnSOP0FakZblXsH3F6aD+Yfng2j6u7ecKW/vb81xvuRuQDREpTMB5ou5LaK6T+5ysIiyZc1sQwzhg+8NSFC2YHty1fXwtDtna356nydwRkA0RK0zmRN2X0MK9uPGFOnUZaJoSTA/uiq/2Y2t/ex6FJDCPZFQsfsJJsPx9fE+wh9vv6IlAFz9BEaavEDEdHR3YXrobV69fR4JEgkX585CsUKCxqQlF6z5E/+RkAMCwQYMx75HZPp8nJGEGY0W9RmvA7gtGzL1X7nBwh8OOUBC2Vbk9RMyyBQvxyoqVOKU7hxs/33JKc7LiFCRiMV59oQC5E3Ow7+ABx3dKRRJefnYFXn52BSdRCo1gnNAp1ZsASsTLwS/kk0KB0CMhYnR6PcZ1BcsaNTIT1fV1oOmIm0GwEltQiPK3xuPToauhrQhsCM1TSwHaymsOHAqXTTjpkRAxzSYTFF2hYGJEIojj49FqNgMAGo1NeHf9Oojj4zFj6jSkDx3G+hxmiJiFefno368fX9ODBtNdMn3yR37n/930RBi1b+F3018Pmw2horW1NaB8PRIihjUUDAXIEhPx24IX0VeSgCvXr2HTZzvwyvKVEMfHuyV3DRHT0tLC1/SgoU5d5tgp4WMXn7zBsiHYCDpEjFwqQ1NXKJhOqxWWtjYkiCWIjY1FX4ltYjxkYAqSFUn4uZHdiStE7E7kE5bXw+4XzckqQuOAarxx5iWPjvBIokdCxGhUapSdqQAAnD1fiRGpaaAoCi2tdxyra0PTbRhuNyJZkcTXpB5DaKeShGYPAFhvXAsoX4+EiJmQPQbbS4rxzgdrIZGIsSh/HgDg4qXL+ProIYgoEUQiCo/NnCX4EDFMf2HehGcc5xiFQJ5a6rBHKH5Na1NTQPmIg91P/PUXchWILwe7v0ITil/z7rUaJKmz/c4XsVuSgcL1UKundNSsfNAUoM+8xck9FKyziv6Uo60ohD7zFmgKbn5NvvX3l2ZJs+9ELPQ6YVr37AKsVttPL3gSQmxBIXasPIuy3CucDl0Ey/HtTzn6+g0oy72CHSvPuvWuXAXOtZ18IYnv7zsRC71OmNUaA6wUjWqN95Wrt57Rn2NzwXJ8+1OOJ/u89aSucG0nXyT2ZfdL+yKq55hsl6z8OZ29pSQJNN0JiorBkjm3WcvzhWsZnvA1x+RaDhNXe/0pg62dAql/rz7B7gk294mvI3LMQ7GuPU8g7hi23svX/I3t+0AON7va688BYLZ26kl3VFQLM5BzmMzGd/3jBFIe2x/Y1zyP7ftADje72uvvAWBf5YWSqB7KAyGQ4crf8lxdPx1ri9C55wvEzH7M8Tuba6gnbAs2YT2P2dMI6TymL2ZuqYOVBkQU8NWSNADu8zdXn6OneTBbWUKn188xhRbr0W7PX6q2ug1/rpe3XG9Jerpcl6eWouDHjdi97RnB1TPY9kSNMIUW69Fuz6iyr/DVkjSnodJ1EeJ6S9LTQmdVjhJ5Nd9BRAuvnsG2J2qEeSunH6wUjVs5wjin6c0h7roIcb0l6W2hI7ST6qFq94gUJttd5e/GHMCO5yvw3ZgDLDm4E6yhydvJdrZ77sxbkp7cONqKQnw6dDXK3xrP+2BGsOrprd3XaA04f8scULkRKUy2u8rBCmLAdSuO7Q/r+pmnuaKrm+aE5XX8u7kcJyyvs35vx7U8LjbwracvvLU7H39nRAqT7a5yIH4+tp6J61Yc29zK/lnnnmLM3FIHs3gx6x+N6Q/8YPd4KK6l479kY6C4lo5/7BqHYbI4Vn+hXQRm8WLM3FKHzj2ebfA152Orp6+IG2x4a3c+/s6IFObAxLiglMPWM7X8+hHsXHkOLb9+xGtetr10+2fVmbcwMf4t7Lj5KusfbVWO0uHu6Uvr8ciODDz1URYe2ZGBJFE16o133RZMQLcIdtx8FRPj30I1y7431zkfWz2DvbOzKkeJkf0DO1/bq/2Y/jqYuew92z+z0jEwplS7lWt3np8fpcCp+7/DzB0ZUNwWgwIFGjSakiz46skqDD88DRMuNCNmtvv5yzVaA+TXVBBR7vveTJsaB1Tzql8wIA52DwTzJLerg5vpCD9heR2lehOevOcdSCxbPR4SacubBhFthZWi0aywQN4lSjs0aBiTLJA1iSGiKVgpEeJL3Q/62p9tFi/GjpuvOsTEtOmNMy8FzSEfaDsSYXqgPW8KKBqgKaBP6ZGAnmf/o5wdNxO/zVjM2qNw3ZUp+7/3477zsTAqLE49pevPpiQL5E1inBvZgXF/O+6xPG/PtfeAf6nailFlX/H6xxloO4Z15+dCbQ3+9tE/8c4Ha3Hw2A9u33d0dGDrF5/jnQ/W4v1/bUAj4x7Id8e+xzsfrMXfPvonqmpr/H62rwl7MM4V2hcUmp/2YPXof7AOc74OOKwouYrf71yBeMoIAG6iBOD0u+K2GAAQTxnx+50rsKLkKmu53p67KkeJ1aP/Ac1Pe3g7wX21YyALJ2+ENUTMjZ9v4bROh1eWP49nFyzCF/v3+T139DVh57qY8YZty5BGTabB46n1VTlKfGn8HC+sfpzVVVNnvItFJ76HWqeEiLYJkdlDsv0uoimodUq8sbESs7792K3MjrVFeGH14/jS+LnHOaG+fgNqMm2i4uOU99WOwV44hTVEjE5fhWyNBrGxsUhWKKBMSsblaw1+PT+QEHm+cP3XH1tQiNNv56B86nWvflI2V035H9+GJW8aXv7pX1Drkp16RW89pu3/bYhoCnn6b2HJm4byP77t9XmuqFOXoXzqdZx+O8ftpJI/PVxPh3oMa4iYZpMJwwYNZuSVotnE/i+OGSIm/+GZWLvpX47vpgAwaoE3tXxrY+OQfA5AibD7ghFG7eauT1MB/AkX64H9h99lzferIWkYe7kW1SNv4tCWX+GE5TV89NMeiGgKD+u/dkrL7CUBuPWcdpifx9C26cQzn9yPSeK3MW3kTagq+6F8SBq+eY/dJk92s9eRH2x/hz+89EpAZYU1RAyXvHZcQ8QEEkuyYfUTUJ64AcOkARj0xmce08ntbpMMOVbluD+HWc5nMz9yuFjuf/FlbNqdBAqdkOI6rmM4qjMNUOuUoGigU2RFjFWETpEVIlrU/Z1Lj/npygqMPzIYKp0SNGi09emAuW8HFLfFMCosuC4eDqn0J5RP6UTZlOt4eu4RnGS4ep746nmf9fRVR3/aKxSENUSMXCp1fG7La3LL6y/ehijliRsQ0RSUJ254LcPuAF+Vo2Qtj1mO69zqQvt80DQA2oo3T70Etc4296NAIcYqcvy0zx/ZWPBRFlRdgqVAIb491rFYUtwW481TLwF0J2iawoX2+QCc53hs9XStB7OOnvDVXsFe8DAJa4gYjSoDp3U6dHR0oLGpCYbGRgxNGcTLHm+TcMOkAbBSNAyTBvAqz16O5Z4E7N72DAp+3OiYWw0Y9hfQiMH4o4OQrWt0mkua+8SAhq3nZM4tXYdw5uKIBnBT0R1kjAKFbF0jFnychXFHB2PAsL8AcD6rabknwa2egSxOfLVXKO8ABcWPeb6mGqXffO0IETN9cq5TiJi7HR3YXlKMhhs3HCFilF0xig78cBQ/VpyGSCTCnIcexr0jVD6f582P2ZPXD+y+PSsFzFm4CTSAXMnbGB77Hzz1UZbT/NDOTUU87mlq9zBhcYYGcHKkAuob1xzuIzsUKFgpGn9//j5UU39AnfEuSv/zNEQefI2h2NXhUiZxsPtBIHMntj+C7jcaqM7ZwnSfHJmEP41dg7+fng1VZTIo2nmVbafbJWT7DXBO1y1kWy7m4sg5v+1zvaYR/yd7DyaJ38JLX5Tb5qFJFvTf/KNHu33RsbYI1j27UK0xoOXXj/CKYtfrr1Yw8TX34TrXZFJSZRu2Sqq6hy0rZfM0UKAw4cJtPC0dC1VlsmMotkM7BmV7PhrwMpTbJWlL1y1c196XAgVVZTKelo7FvXGfOuah8ttiR93Z7PYFvbcYFA2M0CVzfndSsIlKYfqa+/g711yjNTgkQcO2T637jQYZZ2O7ezsaWPjRfTAqLIBIhDr5UHRSIpwcmWQ7AcTYatzxfAXs3geK5T87258/h9tJ7S6CpdBJ2cq3UiIYFRY89fF9mLnDeV7//J/z0bD6CSe7uYrIfkqqRtPI+d1JwSYqXwvNDMfniraiEImmSvSDEvcof+G1HPswaA8XP0n8FjLiduJCnRULdKOdhlS7oORNEsSVHsI7JVdRZ7SdtP9r4lZofvoKe1S/xJ6HVqDOeBcJ6o3I03/jlBfo7jVPjkzCv5vLcfipOFxqvou/VNnK+FL1S/xzwlKkyePw8ZzBkOdNgYh23uIE4BgRns4ai6q783HC8jpKqkychvXYgkI0GL7HiBM0DMX1QJZ/bRwMet0cc0tJEuZ/eJ9tuPURos9+QAKwHZIoOmWbP1Zn2noeVWU/xMy2xfqkv9wFUBSorrvhng5X/H7nCmTE7UTV3fl4Y2Ol27AP2Ib61UszHenenN+9Helarv2ACYamApfrQc3Kx4ma7zDhwm0YFbaDINWZjSgcs8epLr5OG7Xl5XadbqIRX3qUQ8uyE+gcM+p6TF+TfXXqMtRo9kGlU0LEsnfMzP/kPe+gT+tWtCcsxvKH30fbNiVENKDWKUFTgG78LIzt2ubTdkV/U6dOQA489yiLTnwPVeV9qM78HlWZcbhXZ/ue2dtVZZqc0mF+d37Xcstyr0A/tALq1PHIyfq37TPtM1DsmOLYAlVVKvHVO2lY9/V/O+qzRvtnr+1kmDSge4EYQDvzJep6TE89lbaiEImbuwQ5e57H41/MXvJp6VhMODoQIyqViJn9OICuuzI07XDXnH47BzlZRdD9RoMRumTUaBqh+avOo+3M42P/WXEWCz+6z20o/8/z57Dw41GcjpmxPVdbUYjs17QQ0bbVv+hRW6/euedz1GQa8GPudWw2ldu+89J7eludcz3mR1blXXg6TKCv34AROpsbh3nowXV1ac9PAai6Ox8jKm2ngei9xYgtKESf0iMQPfq4Y8GTuHkfANgET1NQ6bqjym0pSXK7ISmaPQ8QiSCaPQ/tCYtZ69CesNgpHRPXcl2f27G2CNmvadGssMBK0RA9avtHSO8thoimMKJSiaq78x2itLcT2yrb2+o81HGMokaY9oYFwLrVZhvCG93uyLiuLu1bdXMypDjZ9jp042e73eOOLSi0zSfRLQhXIXl6VfKj8sfxwRufI7agEMsffh96jcHhCqJBQ68xYPnD7yO2oBBxJYfwzwnPOAnGtVzX59oFKGsSO50ost9H142fjZNtr2NuhtSpndhW2d5W58x7S6FwG0XNUM4lRhAbzLmS6+EHT9cXANtxtsyfvkLl+JkY+/9ecyuXzYnvycYLF9c58l24+5TXxY6vzQFPdq3RGnDj0m9xb5+duDfN1h6eDqN4mjOytaevIb3XL37YFhv6+g0Ye3ggRlRq0TG7iHVeuSpH6fhDtP3Z2fHenf80klUb8SGWOtL+NmMxrOrFtj8Iiz39tT8DNGX76cXGnKwivF7xkuP3uRnOQ6NrHrZymXiya8jWNXi+uhI1mQNRTm1ATlaR00bDqjeUXhcxHWuLkL1Hi8TMgSifusEpNGMo3EZRM5SznZZRpy5zmiP62q1wdbwz88+u/s6p8X3NsdhCuXg60ZMmjwNoGmnyOLfvXPP4ChHjya7Z1d855pj2YdnXRgOzvZhzVOawzuWUUiBE7FD+7vFbnNwVzNt9j8ofZx12vLk+/LkdyMeFEsj7yv15nrd6eCqHOUx/afycUzu4ltXrVuVct8Tsi4jYgkJHb/LkPe84rWy9leVPUP6efjOZP8/zVg/Xcuwr/yfvecfR+zLbMVg2eSNihck2ZNn3sNvzprBeCLMPOxLLVow9PBDZr2nRsbbI67Bsn2fa0zKf5eoO6slQ0J6ex2aX3YU09vBA1kMZzHKYaSWWrV6H6YbVT6AtLxcNq5/walMgROxQzuZg57rdyHRA+9qW9JSWq0OdK4EM5Wyw2WWPWMzcEPCEP2m5bFv2uqGcDU++SldysopsOzldiwhPznC2tHZcHdtCgc0u+4IpZvbjrEJj1t9XWiaB3AjgSlT1mHb8XYRsKUnqcgvZth65RKsI9ktEg9Vj+muX7cUEtq3K8qnXOb1DyJ/2DYsfs9VsxtbiXbhtbEKSXIHFj81Dgtg9utdPZypw4AdbVz99ci7Gd90x/3DLJpjutCAu1mbGc08tRmLfvnxMAuA8AeciTJtbyDZcd+wpxj8nPOMzX2xBIRDGt9p6wh+71mgNWLGnGDFdbqCWp2dxyudv+wYCr6H84PEfoEpNw6srC6BKTcOhY8fc0rSazfj26BH8938tw38/swzfHj2CVkt3lNmFc/Lx8rMr8PKzKziLsuXOJa+LHH8n4PbhupMSYY/ql6wrylCf2A4lnmwv1ZuwR/VLdFLchm473trXHji2YfUT2FKShJY7lwKymZcwdfoqjBs9GgAwbvRonNNfcEtTdbEWqvR0JEgkSBBLoEpPR1VtLZ/Hwtx2i/VAhh17zJ7kGyrHvNGXsGILCvHxHz7HhxOXsja4EF9SzxVPtueppfhw4lJ8/IfPPQ77ru2mrShE8g2VxxhO9uggyhM3QNOdMLfdckvDBV5DecudO5Al2v6IskQpWu60uqWxReqQOX53jbaxY08pRBSFUfeOxPTJuaAo9vuDzEgcj8+chBpNI1SVSsTPnY++ie530ZnbkVT+GpQmzHH8cX43PZX1Gb+bnojfTWev6zxNMnbpGjFPk4xElucFg1CV68l2b/W1U6qvc7TbS6e3Iru4e1ty+uSP3NLfmTsf7bs/w+37+4OizgX89l2fwly3bQtMd+64fT5j6jROD2CNwtHFwrn5kEtlaGtrw+Zdn6Hs7BnH/NMV10gc8i5XCA2gpaXFLT1z3ti+eyfy3ljsmLCzpffFC2NleGGs7R9YIPm5EKpy+djO3AtvX73TsS3Z8vQs9rJWrELcilUYCGAJbIufQPApzOWLlnj8bu/BA2huMUGWKEVziwmJfd1XX3KpDLWX6h2/G00mDB+W6vgOAOLj45GtuQ9XrjV4FKa/5GQVoWN29wqVeViDwB1mu3XMyge9txgxs/KRkxXahR+vOaZGlYGyM2cAAGVnzkCjznBLk5E+HPqLtWi1mNFqMUN/sRYZ6cPRabXiTqtt6O/s7MT5mmoM6M+92+eyGOG6jQZ4Ptgb7fhTby7tGaxFIq855rQHHsDWL3bhZMUpKGRyLHnMdv3gyrUGHC8vw/zZeUiQSDB9ci7e32jbCnvoF1OQIJGgvb0d6z/dhk6rFVbaClVqOnKyx3J+drBdFlyOyEUbno6y8SFYfxdewuwrScAKlqF+SMogzGfEIJqYPQYTs8c4penTpw9WLXsu4Gf7ew7Q150f5pyU3lssSB9lsGEeZWPzYQYSkSNY5zMjdkvS33OAnu782PG09RjN+Np+5BKRw5Vgnc+MWGH6C5d9dH/mpNGAr/pyicgRMtt69GlhJCerCPhrUfc+r9ZAVuleWKM1oFT+OPL+4Ht7NhREbI/JdvaSy4owkndwehKu7cRs8461RWjPmwLdbzS8vRsRK0y2OSOXxuzpw7yRCtd2YrZ5IHNST0SsMNnmjFwaM1SXp6INru3EbPNgzkmj8jxmJBKs85hCg5xgJ0QVRJgEQUKESRAkRJgEQUKESRAkRJgEQUKESRAkRJgEQUKESRAkRJgEQUKESRAkPRIiZv2n23Cp4SrShgzF0iefcnze2NSEbcW70GoxY/CAgVgwJx+xMTF8TCJECSEPEQMAUyfdjwV5c90+33vwAHIn5uDVlQWQiMX48fQpPuYQooiQh4gBAFVaOuL7xDt9RtM0aurrMGpkZlf+LJzTV/ExhxBFhDxEjCdazWZI4sWI6YrUIJfKnELHuMIMEbN47jxIxGKPaSOR/1m+MuCoFUKmvb09NGEI+YaI8QRr6Bj2sEUAnEPErNn4CVYtDfzqrxD5ZPu2qKsTEHi9Qh4ixhN9JQkwt1nQabUiRiSC0dTs6H0JhJCHiPEERVEYMSwVZ89XduWvgEat5mMOIYqI+d/X/vinQDMPGjgQh44fw4EfjqLVbMbcX81AXFwcrlxrwP4jhxxC/WDzv3GivAy3Gg04Xl6GAf36QZmUjMEpKfjywDc4dPwYxPHxmDH1Qc7zrCEDUwI1W7BEY52AwOoVkXd+CNFP9C0DCVEBESZBkEREiBi+W59C4kJtDUq+2Q/aSmNi9hg8+MBkp+87OjqwvXQ3rl6/jgSJBIvy5yFZoQiTtdzwVafaS5dQ+u1+XL95E4vy52F016aKNyKix+S79SkUrFYrivfvw7IFC/HKipU4pTuHGz87B88/WXEKErEYr75QgNyJOdh38ECYrOUGlzopZDI8+egcZGtGcS43IoTJZ+tTSFy+1gBlUjKUiiTExsQgW6OBzmUbVqfXY1xXuO9RIzNRXV8Hmhbu+pRLnZIVCqTcM8Djix/YiAhh8tn6FBK2N3h0byKwbcM2m0xQdMWmjxGJII6PR6vZDKHCpU6BIJg5Zqi2PoUEe79HuaTxb6s23HCpUyAIRpih2voUEnKpFEZGb2I0NUMmTXRJI0OTqRlymQydVissbW2sCz2hwKVOgRARQzmfrU8hMSRlEAyNjWhsakJHZydO63TQqJzrolGpUXamAgBw9nwlRqSm+TU362m41CkQImLn5465FVu/2IWmZqPj7RgJEonT2zEA29bnLYMBbXfbkSCRYP6sR5ExfESYrXfmfE01Sr/5GlbaiglZ2Zg+ORf7Dx/CkJQUaNQZuNvRge0lxWi4cQMSiRiL8udBqUgKt9le8VWny9casOnznTBbLIiLiUViYl8ULl/ptcyIECah9xERQzmh90GESRAkRJgEQUKESRAkRJgEQUKESRAkRJgEQUKESRAkRJgEQUKESRAkRJgEQUKESRAkQTmPyecy0k9nKnDgh6MAgOmTczG+61oBoXfDW5j2y0jPLVwMuVSG9zauR6ZajQH9+jvS2C8jHT5xwilvq9mMb48ewYtLnwUo4L0NtrxCPhhL6Bl4D+V8LiNVXayFKj0dCRIJEsQSqNLTUVVby9ckQhTAu8dku4x0ueGqH3lljLxSjxeZmPExn37sCfTp04eH1cLDfnUk2rBarUhM9P+qBW9h8rmMxHrxygPM+JjR+L7yv69fF7XvKw8oH98H87mMJJfKYDQ1M/KaIJNGX69B8B/ewuRzGSkjfTj0F2vRajGj1WKG/mItMtKH8zWJEAXwHspjRCLMnTET6z/d5riMNKB/f4+Xkc5X6/H10cMoXL4SCRIJpk/OxfsbNwAAHvrFFCRIyIqcEKGX0aJxjvnme+9G7RwzkJcDkJ0fgiAhwiQIEiJMgiAhwiQIEiJMgiAhwiQIEiJMgiAhwiQIEiJMgiAhwiQIEiJMgiAhwiQIEiJMgiAhwiQIEiJMgiAhwiQIEiJMgiAhwiQIkh4JEePpHdyNTU0oWvch+icnAwCGDRqMeY/MDoZJhAinR0LEMN/BfVp3DvsOHsDixx4HACgVSXj52RV8zSBEGT0SIibS3sFNCD89EiLG2zu4G41NeHf9Oojj4zFj6jSkD+J2PxwAAAtlSURBVB3G+hxmiJiFefno368fX9MFRyChVIROa2tg75bvkRAxnt7BLUtMxG8LXkRfSQKuXL+GTZ/twCvLV0IcH++W3DVETEtLC1/TBUc01knQIWLs7+AG4PQO7tjYWPSV2O4cDxmYgmRFEn5uNPA1iRAF9EiIGE/v4G5pveMIXGBoug3D7UYkC/wVyISeoUdCxEzIHoPtJcV454O1jndwA8DFS5fx9dFDEFEiiEQUHps5i4SIIQAgIWIEAwkR45IvBLYQ/GSN1oBD8jlYoyXzaztEmAKgVG8CKJHtJwEAEaYgyFNLAdpq+0kAQIQpCFblKDHNWIJVOcpwmyIYiDBDwBqtATO31DnNGdk+6+myIgkizBBQqjfBSsNpzsj2WU+XFUkQYYaAPLUUIgpOc0a2z3q6rEgiKOcxCb5ZlaMMaA4ZaL5Ih/SYISDUQ61r+dE45yTCDAH+DLXaikJs3i1H6tA3sHm3DNqKQr/Lj8Y5J9mSDAPaikJcqPsE444MglrnPkzrNQaUTbmKe9OWIyeryGd5a7QGlOpNyFNLBTfsB7olSYTZQ2grCqGv34DxR4dh+NnunpRieb0h8/xqdZdIMziKVGiQvXKBo6/fgLGHB2L4WSkoxn9sML9X6ZQYd2Qw9PUbetji8EKEGSS8LUC0FYUYe3ggVDqlQ4w0h/8AOMQ59vBA1vlnNC58ACLMoOFpAaKtKERV3To3UQLOPaPrf67pVDolqurWuYkzGhc+ABFm0PC0Eq+qW4dxRwa7padAOfWMduyfsQ3z444MdhNntDrbyeInxGzaLcdTH4526gXtonTtGe0wRemarlpjQPnU61gy53YP1YAfgS5+whqJAwC+O/Y9fjx9GpSIwtxfzUDG8BHBMEkQaCsKMe7IIMfvbKIE2FfmzPTMdCqdEmem3wmd0QKB91Buj8SxbMFCvLJiJU7pzuHGz7ec0jAjceROzMG+gwcAADd+voXTOh1eWf48nl2wCF/s3xcxPSEXqurWQc2YWwaKa496924Tr/IigbBG4tDpq5Ct0SA2NhbJCgWUScm4fK2Br0mChtlrckkXCNGwUg9rJI5mkwnDBg1m5JWi2cS+umRG4sh/eCbWbvoXX9NDTupQ59+5CM1VkGwCpWnb5TVPHJLPASgRdl8wwqjd7J/RQeYPL70SUL6wRuLgkteOaySOSLhRWPLdVhiTLJDfFjs+cxOZl56RbYFkTLJAJOrjtf5y+xZlhhyrcoTfTmzwFqY/kTjkMplTJA65VOqI0GHLa3LLG8kYTZWQ387yuPoG3IXKtihiClR+Wwyabvf63Gg4KhfWSBwaVQZO63To6OhAY1MTDI2NGJoyiO0xEQr7PjiX+aMnP2dvIayROAb074+szEwUrfsQIpEI+TNmBhyESYhkpD0HvWav08rcdf5o/8x1UeS6l85nMRSJEAd7iGFzsAPsArXjKS0NGnqNAXeenhUxJ43I6aIw48lFc2/ac9BrDKwi87Rn7vqdPT2bKKPBNcQGEWaQ8HSYIierCOVTGlDNIk42fyZzyHbbipzS4NZTkkMcBK94O0yRkfYcyqZcdfvcLkpPP5nYDgs/59dzIxkyx+whtBWFSNy8D6pzyayLGE+LG3tv2RJB80om5GpFhFD9+mQMO+U8UHlbCEWyKAEizIjCfhkNNI1xRwc7XUjTa26jLPcyQIFcRguBPSFF6ML0VyjaikJcuPgJ7k1/LiAhztxSBysNiCjgqyVpQahB8CDuIgHh70o5J6sI9Zf/zHm4di0/GhdARJghgE0owYz25lr+qhwlvlqSJrhhnA9kKO8hfA23nmKwC3mY5gIZygVOMKO99QZItLcegkR78w/SYxIECREmQZAQYRIECREmQZAQYRIECREmQZDwche1ms3YWrwLt41NSJIrsPixeUgQu78996czFTjww1EAwPTJuRjfFfzgwy2bYLrTgrhYmxnPPbUYiX378jGJECXwEubB4z9AlZqGBx+YjIPHfsChY8cw65fTndK0ms349ugRvLj0WYAC3tuwHplqtUPAC+fkY0hU3YwkBANeQ7lOX4Vxo0cDAMaNHo1z+gtuaaou1kKVno4EiQQJYglU6emoqq3l81hCL4BXj9ly5w5kibatMlmiFC13Wt3S2ELIyBy/u4aB2bGnFCKKwqh7R2L65FxQFPsVVWaImIV5+ejfrx8f0wVJYmL0BHuw09rqrgku+BTmum1bYLrjHvZuxtRpnB7g7cL+wrn5kEtlaGtrw+Zdn6Hs7BnH/NMV1xAxLS0tnJ4fSURjnQKNE+BTmMsXLfH43d6DB9DcYoIsUYrmFhMS+7qfIpFLZai9VO/43WgyYfiwVMd3ABAfH49szX24cq3BozAJvQtec0yNKgNlZ84AAMrOnIFGneGWJiN9OPQXa9FqMaPVYob+Yi0y0oej02rFna5uvrOzE+drqjGgf38+5hCiCF5zzGkPPICtX+zCyYpTUMjkWPLY4wCAK9cacLy8DPNn5yFBIsH0ybl4f6PtdSAP/WIKEiQStLe3Y/2n29BptcJKW6FKTUdO9lj+NSJEBeSgsEDwdFA40iEHhQlRBREmQZAQYRIECREmQZAQYRIECREmQZAQYRIECREmQZAQYRIECREmQZAQYRIECREmQZAQYRIECREmQZAQYRIECREmQZAQYRIECREmQZD0SIiY9Z9uw6WGq0gbMhRLn3zK8XljUxO2Fe9Cq8WMwQMGYsGcfMTGxPAxiRAl8Oox7SFiXl1ZAFVqGg4dO8aabuqk+7Egb67b53sPHkDuxBy8urIAErEYP54+xcccQhQR8hAxAKBKS0d8n3inz2iaRk19HUaNzOzKn4Vz+io+5hCiiJCHiPFEq9kMSbwYMV2RGuRSmVPoGFeYIWIWz50HiVjMw3Lh8T/LVwYctULItLe3B3RLMuQhYjzBGjqGPWwRAOcQMWs2foJVS91fkRzJfLJ9W9TVCQi8XiEPEeOJvpIEmNss6LRaESMSwWhqdvS+BELIQ8R4gqIojBiWirPnK7vyV0CjVvMxhxBFxPzva3/8U6CZBw0ciEPHj+HAD0fRajZj7q9mIC4uDleuNWD/kUMOoX6w+d84UV6GW40GHC8vw4B+/aBMSsbglBR8eeAbHDp+DOL4eMyY+iDnedaQgSmBmi1YorFOQGD1isgQMYToJ/qWgYSogAiTIEgi4iWnfLc+hcSF2hqUfLMftJXGxOwxePCByU7fd3R0YHvpbly9fh0JEgkW5c9DskIRJmu54atOtZcuofTb/bh+8yYW5c/D6K5NFW9ERI/Jd+tTKFitVhTv34dlCxbilRUrcUp3Djd+vuWU5mTFKUjEYrz6QgFyJ+Zg38EDYbKWG1zqpJDJ8OSjc5CtGcW53IgQJp+tTyFx+VoDlEnJUCqSEBsTg2yNBjqXbVidXo9xXeG+R43MRHV9HWhauOtTLnVKViiQcs8Ajy9+YCMihMln61NI2N7g0b2JwLYN22wyQdEVmz5GJII4Ph6tZnOP2ukPXOoUCIKZY4Zq61NIsPd7lEsa/7Zqww2XOgWCYIQZqq1PISGXSmFk9CZGUzNk0kSXNDI0mZohl8nQabXC0tbGutATClzqFAgRMZTz2foUEkNSBsHQ2IjGpiZ0dHbitE4Hjcq5LhqVGmVnKgAAZ89XYkRqml9zs56GS50CISJ2fu6YW7H1i11oajY63o6RIJE4vR0DsG193jIY0Ha3HQkSCebPehQZw0eE2XpnztdUo/Sbr2GlrZiQlY3pk3Ox//AhDElJgUadgbsdHdheUoyGGzcgkYixKH8elIqkcJvtFV91unytAZs+3wmzxYK4mFgkJvZF4fKVXsuMCGESeh8RMZQTeh9EmARBQoRJECREmARBQoRJECREmARBQoRJECT/H0DAOoyr5FpLAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "spot_plt = plt.figure(FigureClass=SpotDiagramFigure, opt_model=opm, scale_type=Fit.User_Scale, \n", " user_scale_value=0.1, is_dark=isdark).plot()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }