{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "### Calculation of DOS\n" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "/home/aksenov/simanrc.py was read\n" ] } ], "source": [ "from siman import header\n", "from siman.SSHTools import SSHTools\n", "from siman.calc_manage import add_loop, res_loop\n", "from siman.database import write_database, read_database\n", "from siman.set_functions import read_vasp_sets\n", "from siman.dos_functions import plot_dos\n", "%matplotlib inline\n", "read_database()\n", "header.PATH2PROJECT = 'topologic' # path to project relative to your home folder on cluster\n", "header.PATH2POTENTIALS = '/home/aksenov/scientific_projects/PAW_PBE_VASP' #path to VASP POTENTIALS\n", "header.varset['static'].potdir = {83:'Bi_pv', 34:'Se'} #subfolders with required potentials\n", "header.ssh_object = SSHTools()\n", "header.ssh_object.setup(user=\"d.aksenov\",host=\"10.30.16.168\",pkey=\"/home/aksenov/.ssh/id_rsa\")" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "Attention! You have chosen to override set dos\n", " \n", "\n", "Warning! You did not change NSW in dos set\n", " \n", "\n" ] }, { "data": { "text/plain": [ "{'aks': ,\n", " 'dos': ,\n", " 'ion_relax': ,\n", " 'opt': ,\n", " 'static': }" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#Creating new set 'dos' with parameters for DOS calculations\n", "dos_pack = {'NSW':0, 'LORBIT':12, \n", " 'ISMEAR':-5, 'LAECHG':'.TRUE.', \n", " 'EMIN':-10, 'EMAX':14, \n", " 'NEDOS':2000, 'KSPACING':0.15, \n", " 'savefile':'d'}\n", "read_vasp_sets([('dos', 'static',dos_pack, 'override')]) #new set 'dos' from 'static'" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "-- You provided the following geo file explicitly Bi2Se3_mp-541837_computed.POSCAR ; Version of file does not matter, I use *curver*= 1 as a new version \n", "\n", "-- check_kpoints(): Kpoint mesh is: [12, 12, 12] \n", "\n", "-- check_kpoints(): The actual k-spacings are [ 0.15 0.15 0.15] \n", "\n", "-- Attention! ngkpt = [12, 12, 12] is adopted from struct_des which you provided for it Bi2Se3 and kspacing = 0.15 \n", "\n", " \n", " Calculation ('Bi2Se3', 'dos', 1) successfully created\n", " \n", " \n", "\n" ] }, { "data": { "text/plain": [ "'Bi2Se3'" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "add_loop('Bi2Se3', 'dos', 1, input_geo_file = 'Bi2Se3_mp-541837_computed.POSCAR', run = 0) # calculting DOS" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "-- name | energy(eV)| Vector lenghts (A) | Stresses (MPa) | N MD, N SCF \n", "-- Bi2Se3.dos.1 | -20.1559 | 10.27;10.27;10.27 | -687,-688,-649 | 1,13, 13 \n", "\n", "Database has been successfully updated\n", "\n" ] } ], "source": [ "res_loop('Bi2Se3', ['dos'], [1]) #reading results\n", "write_database() # writing database" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "-- ------Start plot_dos()----- \n", "\n", "-- Bi2Se3.dos.1 e_fermi 4.2271 \n", "\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXm4HFW19t91kpwkZCAESEICAUIShmAioEAIwglwQQhe\nBQUU5foBioAigxAFriThovJJnIIKCHwKXLnoBWRwICBwBCXKEMJMmAJJIHNCBjKdpNf3x+6idlVX\ndVd3VXVVd7+/5znPqarurtqnT1W9tYa9lqgqCCGEkLzRlvUACCGEkCAoUIQQQnIJBYoQQkguoUAR\nQgjJJRQoQgghuYQCRQghJJd0T2vHIjIYwFUAxqnqgcVtkwEMBrAYwMcAXKGqc9MaAyGEkMYlTQtq\nAoB7fdv6qOq3VPUaAHcBuCbF4xNCCGlgUhMoVb0bwFrftim+Y3teJ4QQQhxSc/GVQ0TaAXwZwLlZ\nHJ8QQkj+qbtAiUgPAL8EcKmqzivzPtZgIoSQJkJVpZr31yOL78MBicg2AG4A8CNVnSMiJ5b7oKry\nJ8LPlClTMh9DI/zwe4r+w+sv2g/Pqeg/tZBmFt9hAE4DMERELgPwYwC3ARgD4BciIgC2AXB3WmMg\nhBDSuKQmUKr6GIDHfJs/m9bxCCGENBecqNsEdHR0ZD2EhoDfE0kanlPpIrX6BtNGRDSvYyOk2RGR\nmuMGhARRPKdylyRBCCGEVA0FihBCSC6hQBFCCMklFChCCCG5hAJFCCEkl1CgCCGE5BIKFCGEkFxC\ngSKEEJJLKFCEEEJyCQWKEEJILqFAEUIIySUUKEIIIbmEAkUIISSXUKAIIYTkEgoUIYSQXEKBIoQQ\nkksoUIQQQnIJBYoQQkguoUARQgjJJRQoQgghuYQClQD33QfMmAGsXJn1SAghpHnonvUAGp3XXgOu\nvNIsL1oE/OAH2Y6HEEKaBVpQMbnnHnf5oYeyGwchhDQbFKiYDBiQ9QgIIaQ5oUDFpDudpIQQkgoU\nqJiIZD0CQghpTihQMaFAEUJIOlCgYkKBIoSQdKBAJYxq1iMghJDmIDWBEpHBInKjiDxpbespIteK\nyHdE5CYRGZXW8evF1q3e9S1bshkHIYQ0G2laUBMA3OvbdgGAd1T1agA/BXBzisevC35B6urKZhyE\nENJspCZQqno3gLW+zZMAzCq+/iKAsSLSN60x1AO/QPktKkIIIbVR7xjUIHhFa21xW8Pit5gKhWzG\nQQghzUa9p5kuAdDPWu8PYGnYm6dOnfrhckdHBzo6OtIaV834LSgKFCGEAJ2dnejs7Iy1D9EU085E\n5HAA01X148X1yQAKqjpdRD4C4OeqenjIZzXNsSXF1VcDd97prs+cCWy/fXbjISQJRASNcP2RxqF4\nTlU1MSfNLL7DAJwGYIiIXCYiPQHMALCriFwO4EIAZ6Z1/HrBGBQhhKRDai4+VX0MwGMBL52X1jGz\ngDEoQghJB07UjQljUIQQkg4UqJj4LSi6+AghJBkoUDGhBUUIIelAgYoJkyQIISQdKFAxoQVFCCHp\nQIGKCbP4CCEkHShQMaEFRQgh6UCBigktKEIISQcKVEyYJEEIIelAgYoJLShCCEkHClRMaEERQkg6\nUKBiQguKEELSgQIVE2bxEUJIOlCgYkIXHyGEpAMFKiZ+gWKPN0IISQYKVExYzZwQQtKBAhUDVcag\nCCEkLShQMSgUSl16tKAIISQZKFAx8FtPAC0oQghJCgpUDPzxJ4ACRQghSUGBikGQQNHFRwghyUCB\nigFdfIQQkh4UqBjQgiKEkPSgQMUgyILiRF1CCEkGClQMggSKFhQhhCQDBSoGjEERQkh6UKBiwDRz\nQghJDwpUDGhBEUJIelCgYkALihBC0oMCFQNm8RFCSHpQoGLAeVCEEJIe3bM4qIhcDGBXACsAjAJw\nhqpuymIscaAFRQgh6VF3gRKRwQAuVdXti+v3ADgRwP/UeyxxYQyKEELSIwsX33oAm0Skf3G9L4CX\nMhhHbJjFRwgh6VF3C0pV14rIZAC/F5H3ACwA8Ea9x5EEFChCCEmPLFx84wBcAuCjqqoiMh3AFQC+\n43/v1KlTP1zu6OhAR0dHnUYZDcagCCEkmM7OTnR2dsbaR0WBEpG9AXwFwD4AegOYD+AuVb23xmMO\nA7BC9cNb+SIAuwS90RaoPMIsPkIICcZvVEybNq3qfZQVKBE5GcBJAGYC+CuALgADAXSIyCRVPavq\nIwIPADhWRK4BsBrAGAAX1LCfzKEFRQgh6REqUCLSBgCqelLAy78XkbEiMkZVq0pwUNUCgPOqG2Y+\nYRYfIYSkR6hAqWpBRHqVef35dIbUODBJghBC0qNSmvk3ReTHIjK+LqNpMGhBEUJIelRKkvgqgFcA\nnCIi/wfACwBuU9XVaQ+sEaAFRQgh6VFWoFT12eLiLQBuEZExAG4QkU2q+uXUR5dzaEERQkh6RKok\nISLdROQEAD8GcCyAjamOqkHYsKF0GwWKEEKSoVKa+QwAawCcAeA9ADcC+KyqrqvD2HLPuoBvgWnm\nhBCSDJViUP8HwO0APqWqz6Q/nMYiSKA4UZcQQpKhkkCdo6q/rctIGpAgFx8tKEIISYZKSRK/BQAR\nGQjgPwFsAdAJ4E1VnZv66HLO5s2l2xiDIoSQZIjabuNqAM8BKAD4J5qkEkRcNgakilCgCCEkGaIK\n1FxVvQXAKlVdCdMio+XZFNADmAJFCCHJEFWgxojITgBURLYFMCLFMTUMQS4+xqAIISQZovaD+jWA\np2AqmZ8L4POpjaiBCLKgmMVHCCHJUGke1H6q+qyqPg5gZxHZQVWX12lsuSdIoGhBEUJIMlSyoM4Q\nkbEAHgXwP8zcc1FlDIoQQtKkbAxKVc8DMBHAPwBMFpFHRGSyiAyvy+hyTFAdPoACRfLL5s3A7NnB\nsVNC8kjFGFSxweBDAB4SkR4wtfiuFpG/qur/S3uAeSXIegIoUCS/nH8+8NRTwK67AnfcAfTokfWI\nCClPpCQJERkJYC1My/cRAC5V1XfSHFjeCXsKpUCRPLJ+vREnAHjnHWDePGD06GzHREgloqaZfxdA\nb5hq5kMBTEltRA0CLSjSSCxd6l1///1sxkFINUQVqGcBvAtgL1WdDODV9IbUGFCgSCPxwQfedQoU\naQSiCtRYADMAPCgivcGJuqECxTRzkkf8ZbnCzl9C8kQ1tfheBfADAAfDJE20NGExKE7UJXnEL0hh\nWaiE5IlISRKq+pqI9AJwCIC5MO6+liaoUCxAC4rkE1pQpBGJ2vL9EgA/A/AfAEbBWFQtzaxZ7nLv\n3u4yY1Akj/gFihYUaQSiuvj6qupEAC+r6qMAWj7Eettt7rLduJACRfIILSjSiEQVqG7F344Dq28K\nY2kY/CI0dmz4a4TUG1Vg9WrvNlpQpBGJKlBbReQBAMeLyJ0APqj0gWZm7Vrv+pe/7C4zBkWypFAA\nzj0XOPJIYPp0d3sjWVCPPAKcdhrwu99lPRKSNVGTJKaIyNEw6ebPqWpLZ/HZc0h23hkYONBdZxYf\nyZLZs92KEXfcAXzta0C/fo1lQU2fbiYWv/IK8IlPAEOHZj0ikhVRkyROVtUHVXU6gCUick3K48o1\nq1a5y9ttB3Tr5q7TgiJZMn++d905V/0WU14tqE2bvFUv3n47s6GQHBDVxbeXs6CqzwOQdIbTGNgW\n1IABgFjfBi0okiWLF3vXHXd0o1hQ/goXfnc6aS3KCpSInC8i8wBcKCJvicg8EXkd0TvxNiXr17vL\nffsCbda3SAuKZIm/5l6YQOXZgrJZsyabcZB8UFZoVPVnAH4mIiep6v8mdVARGQ3gCwA2ADgMwFRV\nfTqp/afNli3uco8eXoFiFh/JEvvhCXAFyn/jz2tPKP84aUG1NlGTJDziJCInqOofajmgiLQB+LGq\nHl9cvwXAlvKfyhe2e6R7dwoUyQ9+S6nRXHx+4fSny5PWImo/qCNh6vDtABN/6g+gJoEC8HGzSzkP\nQB8AK1T1xhr3lQm2BdW9uzcGRYEiWeIXonXrgrdvyekjIS0oYhM1lvQFAMcAOAumJ9S3YhxzV5iC\ns6eo6joRuU1ENqnqrTH2WVfsp88ePbxZfBQokiVhAuW/8TeKQOU1VkbqQ1SBmquqq0Sku6p2ich2\nMY65BsCrqlq8dPB3AB0ASgRq6tSpHy53dHSgo6MjxmGTwy9QtKBIXrDLbgHhFlSjuPjCijKT/NPZ\n2YnOzs5Y+4gqUIeLyDMAeonITTAFY2vlXwC2FxFRVYWxqF4LeqMtUHnC7+JjDIrUwpo1wFVXmdTq\nqVOTmZDabC4+ClTj4jcqpk2bVvU+os6DOgXALJg41PMAvlr1kYqo6ioA34bJDvwuTFzrJ7XuLwv8\nFhTTzEkt3H67Keszezbw058ms88wC6pR+kHRgiI2US2oSar6ewAQkU4Ygbqk1oOq6r0A7q3181lT\nTqA4UZdE5eGH3eVHHklmn41uQVGgiE2tlSRaGvtptFcvWlCkNvxWjX8OU7WoRk8zz6tA0cVHbGqt\nJNGjPsPLJ/ZF07MnY1CkelS9NR2B+HXngjLe1q0zx/JbJnkVKP84mcXX2pQVKFX9maruDuAsVR2h\nqrur6ihVvaBO48slFCgSl/XrS+NFb74Zb59B1sa6dcHxprwKFC0oYhPVxfesiAwWkYEicoGI7Jrq\nqHKOfdH06sU0c1I9K1eWbps3L94+/YIHAB98EHyTz2uShF+ggv4m0jpEFajvAugNM0l3KIApqY2o\nAbAvmt692W6DVE9QhYQ33oi3zyAhKhRKK4QD+bWg/C6+vAopqQ+RLSgA7wLYS1UnA3g1vSHlH/vm\n0q8f222Q6gkSqJdfjveAE+YOW768dFteBcpvQRUK9Eq0MlHTzMcCmAHgQRHpDWBEekPKP/bNhe02\nSC0EtZF4/33zs12VdVqefx548EFgl12CX1+xonRbV5c5VyVnnd2Cqqxv3mxc6aT1iCpQVwM4FsCN\nMHX0Wrrlu9+CYpIE8bNxI/DYY8C++wZXiAgrgrpgQXUCtWoV8I1vlE9RDxIowJyrtns6D1CgiE2o\nQBXbYuyuqm+q6mtwyxE9Wny9F4AdVXVB+sPMD6ru5EfAWFC2KFGgiCpw1lnGZbfttsCdd5aKTphA\nLVwIjB0b/VidnZXnTwW5+ABjReVNoBopoYOkT2gMSlULAM4UkS+JyLbOdhHpJiLjAdwMICD82txs\n2OCKUK9ebFhISnn7bSNOgOln9Pjjpe8pJ1DVsGhR5feEWVB5jEOFWVCkNamUJPGfMMVcXxSRtSKy\nCkaUzgMwWVVbrluLbT3162d+U6CIzQsveNeXLCl9jx2D2nNPd7lagQpKV/fT6AJFC6p1qTRRt6Cq\n31PVXQAMAzAaQH9VPVVV363LCHOGP0EC4Dwo4mXxYu96kAvOftDZZx93uVqBsvcTRjkXX94IqhyR\nx3GS+hA1SQKqGpB31Hr4EyQArwUF5DM7itQPfwmjDz4ofY9tQY0e7S77xa0SUTrOhllZjWJB0cXX\nukSdB0WKBFlQAFPNiYvfpRYkUPZ5tMce7gPNsmXVWQxhFtT224ePxyGPlklQ5Ygkx6kK/OpXwLe/\nDbz3XnL7JekQ2YIihiALCiidrOu3qkjr4LdYglx89nk0cCCw447A0qXmBrp0KTBsWLRjhQnUkCGu\nMIU9MOXRgkpboJ57zggUYCxdZ5nkk6puoyIyRkQ+mtZgGoGgJAmA5Y6Ii1+gKllQ/foZQXGIkpkX\ntB+bnXaq/NlGEagkXXzPW82CZs/OpxVJXCq12/iNiPyquHwqgL8DuEVELq/H4PJImEAxUaJ12bwZ\nuP564NZbzf++kkCpemNQ/fp5BaUagQqzoBpRoIL6WQHJiog/CWXZsuT2TZKnkotvo6qeXVz+FoAT\nVLVTRG5IeVy5JczFx1Tz1uXnPzft2wEzN84vSH7B2rDBrdnYsyfQ3l6bQG3eHG5d2BZZGHkTqM2b\ng6+dJC0of0bjsmXBlT5IPqjk4nsPAERkKIAhqtppb29FoiRJUKBaizvvdJdvv73Uxbt8uVcM7Cw/\np8KELVBRM/nKpZhHsaDy5t4Ka62RpED5vzNaUPmmkkANF5G+MBNz77W275jekPINBYr4sW+gQfOY\nVL1P7rZADRhgftdiQZUTqBEjgqc6tLe7y3mzoMIEKslx+ov0hs0RI/mgkkD9EcA8AJMAfA8AROSv\noEABYAyKRLdC7GoSQRZULUkS5QRq0CBg8ODS7TtaV26jCFSSFpRfoIKqypP8UDYGpar3ALjHt+2o\nVEeUc6Jk8VGgWoeoLqJ33zV1G3/zG7dOHxDs4luyxJxDlaYqlBOo9nZg551L3YXDhpmxAPlz8YX9\nPWm6+ChQ+abiPCgR2QbAZwGMBLAQwB2tWIPPIYoFxTTz1iGozl4QTz8N3HQTMH++d/ugQeZ3796m\n8vnq1UY4VqzwWjtBVCpzNGyYOa5Dv35et3TeLKiwlPmkhHTLltI5aRSofFMpzXxXAP8AcASAdgCH\nAPiHiIS0Rmt+omTxsatu67B0afhr3a3Hv/vuKxUnwOvaqzYOFSZQH/uY+e2f7Dt0qLHiHPJmQYWJ\nRVIWVND3tXp1Mvsm6VApBvUtAMeo6umqeqmqng7g3wBckP7Q8kdQLygHljpqTcpZUAceaMoYlcOO\nE1WbyWc/LB14oKlI0asXcNppZtu4cd73H364VzRbzYIKEkBaUPmmkkCtUFXPM6KqLgHQkv/WTZvc\ni7q93ZsRxSy+1qSckPTtC3zxi+U/P3Kku2wLVJSq5vYNfdw44P77gZkzgQkTzLb99wcmTjTLvXsD\nxx+fb4GyxcKO6SYlUEECSAsq31QSqDBboCVthLAUc4AC1aqUs6AGDQImTfJWK7c57DCvKO22m7v8\nwAPh59EHHwAXXWRiWg79+plJv336uNtEgO9/H/jlL4G77sq/i8+e0GzH35wJvK+/Hs/dFyRQUarB\nk+yoJFADipN0P0REhgEYkN6Q8ktY/AmgQLUq5WJQgwYZS+CHPzQi1dYGfPObwF/+Alx3HXD11d73\nH3WUsXQA4K23TGHTIO65B3jsMe+2/v2D39ujh3H/OckYebag7Ljb8OHu8saNwNlnA1/4AvCVr9R+\nfYW5+OiSzy+VsviuAfCAiMwFsAjAUACjAByb9sDySFiKOcAYVKtSLpnBSYDYeWe3woST7RmUode/\nP3DEEcCf/mTWX3oJ2G+/0vf95Cel2/znYxi2BZU3gXrXaoG6227Ak0+a5UcecR8OX34ZePFFYOzY\n6vcflCSxdauxJG13PckPlTrqLgJwMEwViSUwc6LGF7e3HOUsKH+7DdL8bNgAvP++We7e3RtPAkrX\nozSx3Hdfd/mNN0pff/jh4M9FFai8WlDr1wNvv+2u2+5OvxuunNVajrCEiKBq8yQfRGm3sQnAMgBL\ni78DmjK3BmEZfEB67Tbefx946CFmG+URO0Fi8GBgzBh3vXdvYzlVi/2ZIOvsd78L/tyAiE53W6Dy\nFIN6+GHXdTdqFLDDDuHvDesQXImwayioXxfJB5XmQW0L4CkAUwEcBWAagKdFJHYMSkR6ichzIvLD\nuPuqF/aTll+g0ih1pAqcey5w6aUmKE7yhS0gO+0EfPazrqvoC1+orWmlPS8qKEPQtjIchg4Fdt01\n2v7z6uKzC+4ec4w32cNPrQIVNm+MFlR+qRSD+g6Ay1X1AWeDiBwN4FIA34557KsAzI65j7pin+D+\nCyiNJInly4HXXjPLc+aYJ8CwYDipH3/+s0lU6NnT3bbTTsA++wB3321uoHvvXdu+bYEKKnlkn4O3\n3grMmgUceaTXgi9HHl18ixaZeBtgBPSEE7zxKD+0oFqHSgIltjgBgKo+KCKHxzmoiHwJpvnhOABl\nnpXyRb0F6j1fU5NlyyhQWbNhA3DVVaXpzk66+JAh0XoxhWGXPNqyxVvyqKvLPW5bmxHBffapbv95\ndPG98oq7/NGPmr+/nEt75Upg7lzg738387xGjIh2nLCUcgpUfqnkhAj719VsFIvI3gD2KhaibSjK\nxaDSECj/U2RYtWdSP95+O3guTlQXWxTCSh75XcxRki785NGCevVVd9mxPMs9iM2bB5x1lknV//rX\no19vtkDZFTzo4ssvlSyoY4r9oPwcDOD7NR7zBAAbReTbAA4F0ENEvqmqM2rcX92wT2S/BZVGDMpf\nTYBPetljt8qwifoUH4UhQ9yb9uLFbkp1uQekqOTdgtpzT/O7XFbiO++4y8uWmf/J9ttXPo5tlQ0Z\n4k6y5nWVXyoJ1GYEW0s1n9qq+qGwiUgvAH3CxGnq1KkfLnd0dKCjo6PWwyZCuSSJNNptvPiid50W\nVPYElcYZODBZgQqryZeEQOUtSWL9em/FdSfNvq3NuDajtDNZvjyaQNkW1JAh7kRoWlDp0NnZic7O\nzlj7qCRQk1X1Kf9GETkg1lHNPk4E8AkA7SJyiqqWJNDaApUHyt0g0mi34cyxceCTXvYExUY+/eno\nSQpRCMvkS9qCyoNA3XOPa8kNHWp+HI491iSCAMAnPwn89a/BY16+3LW8wlANd/HxukoHv1Exbdq0\nqvdRqWFhiTgVtz9T9ZFK93E3gLvj7qeelHPxpdFuw58WSwsqe4IsqBNPTPYY9bKgsnbxFQrAHXe4\n61/+svdB7ytfMYLa1QWcfjowe3bwJN0obds3bHA9G716mUQMB1pQ+aViw0LiEjVJIikLyn/h8Ekv\ne/wW1FVXeQUlCewySPYNudksqFdfdTNV+/c3hXVtttnGzAN02G672gXK/r/172/27UCByi81TCVs\nXeptQW3y1eygQGWPbUFNnWpcT0njFHYFmlug5sxxlydMMJZNOQYODN6+YkXptkLBzCF0rhlboPr1\n816/vK7yCwWqCqLOg0rKgvLfQDZudJefeMJUxn7gAZA6YguU7SZKku23d11dq1a550GzJUnYWaqV\nYkhAuEAFWVBTpgCnngqcfLKJPfkLPdsWFAUqv9DFFxH/JEn/014a86D8MQJHoNasAS64wBznmWeA\njo7KT58kGZIQiUp0727cWStXmoedFStMUL/ZLCg7vhbFTRomUH6338KFpqWJc4xHHvE+TNDF1zjQ\ngopIpUmSSQuUaqmr0EmSePRR9xibNtVe3ZlUj/20Xa5eXFzsLDNnwnazJUlUK1DbbVd5P4B5aPOv\n+5uN0sXXGFCgIlJuDhSQ/ETdoKdbx4LyX4BRgsQkGeybmdNcMA3seVVXXmm6yTabBWVXyYhSHips\nrtPy5d7qHk5dP4f5870CRQuqcaBARaRc/AlIvt1GOYGyg8tAtMmMJBnsVP80LSi7xt7ChcCXvmRq\nzzk0ciWJDRuAzk5XNNrbw60jm3Lvsb0I/gnu8+eXZvGVs6Def9+kvwf14yL1hTGoiJTL4AOSb1gY\nJlAbNgQXkSX1oV4W1KGHAtdc465v3eo9rxrVgioUTP285593tw0ZEq2uYFgMCjDXxM47G0vT6QDg\nsGaNt65l377hSRKFgmkv/8Yb5v97993B3Y9JfaAFFZFK7pWks/iCnm43bAhOqaVA1YdCwU39F/G2\n20iaYcOA444Lf71RBeqNN7ziBHi755ajnEA59fm+973g1+1jDhhQKlCOW37BAtdy2rABePzxaGNL\nikIh2YanjQ4FKiLVCFSaMagggWIMqj7Y7r3evWtrSFgNU6cCF18c/FqjJknMnVu6bcKEaJ/1u/j2\n2MNdnjfPuPls957dlde2oHbYwfzvbAvY+d/6Owg49frqwYsvmmaNJ5/MxCcHClREKsWg6iFQ69fT\ngsoS282bpnvPoa0NOOWU4GPVGv/K0oJ64gnAX45t0CDg+OOjfd4WV8D0jnJ46y1v24799w/fr5Mx\naH+Hzv/WX62+ngJ1ww3m+PPmATfdVL/j5hkKVETykMW3alWwtUQLqj7UK0HCRsSbcu5sq1Ugs0qS\nWLECmDzZu238eODGG6tzlZ5zjvk9aJARbwe/QO21FzB8eOnnu3d3C9IGCZS/W+/CheGNDpNm1ix3\n+cEH63PMvMMkiYhUcvEl3W4jLAa1YEHp9qDW4CR56pUg4WfQINMo0aFPn9r/11lVkpg501sJ5eij\nTbyo2qaLp59uuugOGWImp/fsaeKCK1d63Xt77BEsUMOHu9dqFIECTHzLaQOSFv57RlBTzFaEt7SI\n2K0vgpqp1cOCAoJ9+Js2eeeUkHSwLSg7yJ42/iyyOBUssopB2ZmnkybVJk6AEeYRI8z339bm7WRs\nu+N23BHYfffSY3zsY+6y/T2WEyj74SAt/K11Nm/m/CyAAhUZp/sm4C3m6ZB0DCrs5jF7dvD2N9+M\nf0xSniwtKJs4AtXe7p6rXV31s6Ls2On48bWJUxB2IoR9Q99hB1PeaNw47/vHj3eXo8SggPoIVFAc\nOSje3GpQoCJiC1TQrPek08wr3Tja2oAjjnDXOakwfWyBalQLyh+/qleZH/tmG6X7bVRsgbJxjnHO\nOe61uccewCGHuO+xBcpx4QcJ1Ouvxx9nJYKy9vxWVSvCGFQECgWvi6KSQCUxUbeS++Xgg83T4COP\nmPUg1x9JlmYQKMDcmB2LYf16U1khbexEnrQFqq3NTUk/4ADgt78186A6Oryx4iALKqhj8pNPmnp/\nUcox1UqQBRXUHLPVoAUVgffecydoDhwYfEHX24I65hhgzBh3/dlnOcEvbbISKL+Lb8CAePvL2oIK\ns3pqIWhfAwZ4r8dRo4DPfrZUGCsJlJNk0dUF/PrXyYw3jCALigJFgYqEHd8ZOTL4PWnGoPwpze3t\n5mlw5Eg3YWPlSjPP5OyzzeROf7t4Ep+8JEmUq6gQhXpX8t640RWA7t2Dk4xqJUigogqgP0miUPBe\nN3Za/IMPJteINAgKVDAUqAjYAmXPXrdJc6Lufvt5rbaJE91U4/32c7effz7w9NOmEOdtt8UfA/GS\nVZKE/8k/rouv3haUP/6UVIKEsz8/UQXKb0GtW+d6Ifr0AQ46yLVe164FXnkl3ljLQRdfMBSoCNjx\nHbsNgk2aaea9ewMXXWQsp112Ac47z33NTpu1cRq2kco8/7wR9Ur/t6wsKDtuAsSvARjk2kqTtBIk\ngGAxinoUeibfAAAZOklEQVQMf5KEvyWHCPDxj7vbnnyytjFGgUkSwVCgKqDqnV9hx31s0my30aOH\nKdvS2WmqK9vB2kMPDf78okXJdfZtZmbPBs4807hFf/e78u/NyoICgM98xvzu1g048sh4+7K7yz71\nVLx9RcFOkEgy/gQEi1EtAvXBB16LxXFDHnigu83fhy1JaEEFQ4GqwLvvuhdYnz7hMaik223YMSin\nPE17e6l7ZPjwYKtOlTX6onDdde4DxS23lH9vpZYraTJ5MvDjHwN33RU/m8wuujpzZrx9RSFNC6pX\nr9L/Ra0CZVtQjojb9f5eeimdh76NG4OzBylQFKiKPPqou7zffuElZtLM4uteYTLAiScGb58/P/44\nmh27+2qlmoZZClR7O3DYYabnUVwOPthdXr06+OaYJGkKFFBqlUU9hj9Jwv4eHAtq6FA3a3LdunSu\nqbDK5RQoClRZVL1FG+2JsX7SzOKrJFAnnwyccUapJcWS/eXZurW05lm5h4ssBSpJbLcVUNpiImnq\nLVBOtfJK+C0of9ddwHgs7Dp8/m69SbB4sbtsTymgQFGgyvL8827mTo8eJrU7jDSz+CoJVFsbcO65\nwO9/D5x6qrs9qK4YcbGrgziUSxqwY1CNLFCAN3bp79CcNGkLlH+ftQjUunXBAgWkL1B2KaWxY93l\n99/n3EYKVBluvdVdPu648jPu07Sg/H1wymHHFyhQ5Vm4sHRbOTdfs1hQgDeOlXasMm2B8k9kjjpP\nzM7E9FtQdiLJRz7iLqchUG+95S7vuad7vW/e7K0A34pQoEJ48UXgb39z17/4xfLvz9KCsrFvABSo\n8gQJlF2L7c03TXbfb35jnmQrNa1sJOybeNrnif2dxp1kHMSxx7rX3+c+F70VSXu7+QGMu9d2idsP\no/vs4y6/9lryomFnB44Z4xXHVnfzsRZfAKomY8rhyCPD5z85JD0PihZU+lQSqO9+19yQOjuB0aOz\nmweVBvZ5kvZ8G/s8TEOg9twTuP1247L1x9cq0aePG4e0XZ22SPTrB+y2m3HFFQqmMaKd3ReHTZtc\nF19bm3Hxbbuta8mvXp1uDcC8QwsqgJkzTfwJMNaLPTE2jKQbFtKCSp+g5ABHoAoFI04O99zjxgMG\nDmz85pD1sqDWr3ctjvb29IR95EhgwoTqHuYAryUcJlBAenGoBQvc82qnnUzaPC0olwa/zJJn1Spg\n+nR3/eSTo6X22hZUGhN1o2IXEg1qHUBcyllQfqvCnm7QDE+09bK0/e69JMscJYEtUHaszB9vTisO\n9c477rLTfNG+hlu9mgQFykIVuOYa96QYPBj42teifTbpdhu1WlD205c98ZB4UQ2e0+LcrP3ibj90\nNJtApXkTTNu9F5ewWGK9LCg7QWK33cxv+3tq9aaFdY9BicgIAFcBeAbALgBWqOp/1XscQfzhD955\nT5ddFj0YnuZE3WosKKcVdqFgYiabN7uBYOLyzjvemJKDc7Mud2NoBoGql4vPFnpbFPNCWOsSv0CN\nHGlqIG7aZOYtrVwZX3DnzgVuuMFddywoe15XpcnjzU4WFtRAAP+jqj9S1QsAfF5E9qv0obT5y1+A\nq6921//9341POypZTtS1EfG6J2hFBRNWTNe5oZa7aUedZ5Nn+vVzz9l16yo3yKyVtFPM47LLLqXb\nevYsLcjbrZtJxnBIorL5jTd61x2BsturUKDqjKo+rar3W5sEQB1qKofz5z8DU6a4wjJ6tLcXTBTy\nIlCAV6DSLmPTqNipvWec4S47wlROoMJarjQSdtdZIL14Zd4tKEcUbPzzqhz23ttdfvnleMddu9Zk\nh9o455VtQbV6Pc1MY1Ai8hkAM1X1tYpvToGlS4FLLgGuuMIVlT32AH7xC5NNUw1JC5Rdgqfa9gpM\nlKiMPXvfrqoQxYJqBoEC6iNQeY9BOV1zbaII1A03AJ//vHcyfzXcdZd3/aij3P8HXXwumc2DEpEO\nAB1FN18gU6dO/XC5o6MDHeVqDVXJnDmmwZ9dHWCPPYDrr6/tSS/pGJQtUNXGkGwXQVA5n1bn/ffd\nWFPPnsBee3lfKxTCBWrAgHzeaGvB/jvSEij7BptHC2rUKDdm6zB4cPB77Qm7APDGG8CMGaaItJ3l\nFwXbAjvlFDMh3KFZBKqzsxOdfjOxSjIRKBGZBOBQVb1ARHYCsKuq/tP/PlugkkQV+P73veJ0/PHA\nhReWBkejknS7jTgCZQfxKVCl2Nl7u+5qvt++fU0splAw7pewG8P+++cvVbpW6jHfxj7/8phc0qeP\nER47M8/JpvOz227Gs+KvJNHZWb1A2Q9AEyd6z6nttnNFc82axk108hsV06ZNq3ofdXfxicj+AO4A\ncJCIPArgHgCj6zmGefO86Z3XXgtMnVq7OAHJNyzctMldrvbktJ8A7UrJxGCXtBk61Pz2zwsKE3a7\nw2qjY8cq6yFQYZZJ1pxwgnfdToawaWsLfs1u2RKVcskjbW3eba2cap5FksRsVe2nqkeo6kRVPUhV\na/Tk1sazz7rLEycC48fH32fSpY7ixKBoQZXHDjw77lB/PCYoON3eHr+bbZ6wH8jSSKbZutW1REW8\nruc8cdxx7jyn4cOBAw4If+9RR5VuW7Cg+mNWym5sFjdfXFqyFp/dwn2/hBLck06SiGNB2QJFC6qU\nSgL13nvuDbt7d+Dyy4HHHwc+9anmiT8B6QvUsmXutTBwYH7dVD16mJTv2bNNPLJcgtTnPgfMmgU8\n8YS7bdkyM28xarbtxo1u65bu3d3miDY77uimslOgWow33nCX/YHPWkkziy+OQC1aFH8szYZ9wTtP\nqvZT7OzZ7vJOOxlh+tSn6jO2epK2i89+OMqre8+hRw/goIOivW/GDLP8yU+ac6lQMG5jx11cCX9m\nY1BMk6nmhpYrdVQoeFOMK1Upj0qe0sy32859CvQ3YiPBFpQdGLcTj6LedBqRtC0o++EojwkScbEn\nbFfzIBhl8rL9faXd8TjPtJxALVzo3vx32KF8E8JqSDoGFcfFJ+ItcFuLj7yZCRIoO9XcvlkPG1af\nMWVB2haU/SAYVLGh0QkTqKVLvb3D/NgCFeYytq/foKLGrULLCZSdvZeU9QQk327DriRRrQUF8AQv\nh+3icwRq3LjgSZtRKtk3KmlXHJk3z11O8lrLC7aV47Tq+PnPgUmTTKm0sGoTtosvzILi9WtoOYF6\n8013OcmLJul2G3EsKIAneBgbNrhPtz16uDfpbt3MhEk/ScUo80jaLj5boHbfPfn9Z419jb37rvl7\nne7La9aY2p5B94IoLj7/vpO4pzQiLS1QI0cmt98k221s3eruQ8RrnUWFLr5g/NaT/WBx3HHe4HTf\nvtVPwGwk7Oyx1auTvQlu2eKdEB02+bWR8T8E+ssevfxy8Bwp2x0Ylnrfv7/7/9m4sXUz+VpOoOwM\nvrQEKu6F7s/gq6Vyge3zpwXlYsefbDECzA3hV78ytdja24Ef/rA292qj0N4O9O5tlrdudVOfk2DB\nAvcha8iQ9DrpZoktUM89Z4pO+/njH0u32Z17y7mQbavz1VerH18z0FICtXmzt4Nlki6+JLP44mTw\nOdDFF0xQ/Mlm+HBzo3niCeDAA+s3rqxIK1Gi2d17gEmdt+c+2V4Ph5kzvdcz4LUsyyXh2E0SH3oI\nuOOO4CabzUxLCdQ777gn0dChyT7VJSlQceNPgHlqdVyDy5cHN+drRYIy+FqZtOJQrSBQ3bp5sz8d\nbrjBnZ6wdi3wt7+5ry1f7pba6tmzvAVlC9Sf/wxMnw588Yut5bJvKYFKK/4EpGdB1SpQ3bp55/C0\n8lwKm7w30Ks3aWXy2dmyzSpQQGnpo4suMgWFjz/e3ea4+V55BfjWt9zte+9dPr58yCGlD9EbNgD3\n3x/8/makpQTKjj8l3dMnrRhUnBgI3XylREnxbSXSqmje7CnmDqecApx4otvk9NRTzXZboGbNMqWy\nzjzTmzQxdmz5ffftC5x0Uun2f/0r/rgbhZYqdWTPSxg1Ktl9J9luIwkXH0CBCsIWqDz2J6o3abj4\n/NVamtmC6tEDuOyy0u1DhxpLavZs831ceGHpe4IKz/o5+2xjZT32mPuA/frr1dX+a2RaxoIqFLw9\nX5JOH06y3UZSFpSdyddKfuty2I35mqnwa62kkSSxaJF7Dg8cmFy1lkajXP3G8eO9HXrD6NEDOPdc\nkyDh1DPcvLl1kiVaRqBefdVNo91xx+RrgyVpQSURgwJoQQVBC8pLGjGoVok/VeLII0sro48YYaYy\n/PSn1U8fscMS9nfczLSMQNmZNAcemHxX1DTnQdUKBcqLqteCokCl4+JrlfhTJbbZprR/2NlnG9df\nLZPv7e+yVQSqab2YmzebdM9HHjHrtotrwoTkj5e3NHPAO8di0aLW8VuHsX69K/69e7uTVFuZNFx8\nrZBiHpUzzzRxqMWLTS+piRNr35ctUPZ33Mw07e3ql78E/vu/S7cPGwZ0dCR/vLylmQMmfjVokJl3\nUSgYkWrGqtJRsVPMaT0Z0sjio0C5DB8O3HuveTjq2zfevuzv0p4y08w0pYuvUAieKzB6NHDttel0\n9kyy3UZSSRJAaTmWVoYJEqUkHYNSTa9jQKPS1hZfnADvdzl/vvGINDtNaUG9+KL3afAnPzE+3wMP\nTM/FlWS7jY0b3eW4AjVunNshdto042o4+eTWzKzydzIlybv4li1zk5H69eP3nCR9+rgekS1bjKWa\n9HSZvNGUFtTrr7vLxx0HfOITZlZ2mvGXtCyouNbesceaVFXAPN1efz1w9NHAV78KXHedd75Ks8ME\niVL8SRJxE3z87r2kk5FaHbv80VNPZTeOetGUAmWXs69XzCVvxWIdRowwTdRs//WWLcCzzwI332wC\nt+ecA9x2G3DffUbcm7X3DC2oUnr2dB+Curq81nstMMU8XQ46yF3+5z+zG0e9aEoXn13O3m7LnCZJ\nppknlcXncMABZqLfn/4E3HlnaafPp57yPo2NG2fcos3mBuQcqGAGDHALmK5ZEy+7kSnm6XLwwe7y\nM8+Yh9k0Yup5oektKLtgapok2bAwSQvKoVs304b61luBBx803T4nTvSO2+G554AZM5I5bp5gHb5g\nkoxD2dlltKCSZ9gwkxkImAfZF17Idjxp05QClbUFleQ8qDQa5g0caOqAXXONyXY8/3xT5PKQQ9z3\n3Hef+1TdLNCCCiapTL5CAXjtNXe92QP4WfHRj7rLc+dmN4560HQuvs2b3fkubW0m66UeOIkIgPHl\nxyFpF185Bg8GTjvNXf/qV018qlAAHngA+I//SPf49YStNoJJqprE/Plu37GBA0s7FpNkGD3aXbYT\nwpqRprOgFi92lwcPrq2kSC3YNbfiBprTcPFF5dOfdpdnzqzvsdOGSRLBJOXis9uS7703M/jSwhao\nZregmk6gsnDvAckKVNouvnJMnOhabXPnNk8aeleX6W4KGMvathpanaSqSdh9iqJU6ia1YTdbnTev\nuSfsNrVA1StBAkhPoOqdodOnD3Dooe56s3TvtOdADRgQnBzSqiQRg5o502SJOtjnEEmW/v3dbgxd\nXc3zEBlE012mdtXuegpU9+7uTW/r1nhPNc6TPpBMiZRqsfvY3H9/czyhMf4Uji1QtpBHQdVkhl5+\nuZsc9JGPAGPGJDc+UkqruPmaTqDsRl5OOmY9EEnOirKfYrNwRR1yiBvgXrkSeOKJ+o8haWyBYvzJ\ni90bzZ6iUQlVU9vSnpIwYoSZwsD4U7rstZe7/Mor2Y0jbTIRKBE5UkR+ISJXiMgVSe47K4EC0hGo\nKJNlOzs7az9YAN26Accf767/8Y+J7j4TFi4E1q7tBFBfy7oRqKVvWKEA/OAHxnpy2H9/4Kab3M6v\nrUDS115U7BjfSy9lMoS6UHeBEpHeAK4HcL6qXglgrIjE6JLi0tXlvcDsC68eJCFQqtkLFABMmuQu\nP/54cs3ssmL+fFegWrnlSBBDhrjZrosXe61NP4WCmRz6jW8Ad9/tbj/8cFNSq9mqj1QiK4Had183\npPDCC8CcOZkMI3WymAc1HsDbqupENv4BYBKAR2vZ2dlnu+nDy5a5KdpDh9b/YrEF6jOf8bZojsq6\ndW7Mp709uzImu+8O7LOPKYvU1WX+Hv+8Fqekk13aKWhb3PdX+myU99s33d12Kz1WK9PebuJGzk3u\npJPM75NPNr+7usx11dVlKpX7H76OOw6YMqV+UzqImWh+6KHAY4+Z9bPOMg/k9nzMvHDzzbXH0rMQ\nqEEArDQArAGwX607e/ttYPny0u3771/rHmtn++29xTLjNhWrtwXo54QT3Lp9a9Y0vhUFmJvxxz6W\n9Sjyx6c/7QqU83+u1FZcBDj9dPOQyKzI+nPRRWZS/dq1xrK1wxt5Ik5lHdE6l64WkSMAXKaqRxXX\nLwQwTFUv9r2vSWtqE0JIa6KqVaXPZGFBzQIwXER6qGoXgAkAfuF/U7V/CCGEkOai7hYUYLL4AJwE\nYCmALlX9r7oPghBCSK7JRKAIIYSQSuSumrmICICvAvgvABNV9WXrtXkAnJZo76rqaQG7aAkqfE9f\nhEk82QLgLVX9VTajzB8iMgXA4dam76nqw1mNJ28UvRsnAlgCAMWpICQAEZkFoFi/HVtV9d+yHE+e\nEJHBAK4CME5VDyxu6wlgOoB3AYwE8H9VtWw99twJFICxAP4F4IOA137NC+ZDAr8nERkG4GJV3a+4\n/qSIPKyqMXMKmwZV1SOyHkQeseYo7q2qW0TkThGZqKo1TQFpAf7C+1EoEwDcC2Ccte0CAO+o6nQR\n2RfAzQAOK7eT3AmUqj4HfGgh+DlMRC4G0B/m5JhV18HliDLf0zEAnrbWZwE4FsDP6zS0vCMichmA\nTQC6AbhWVTdU+EyrkOgcxRZgrIhcAmAbAE+p6p+zHlBeUNW7ReRw3+ZJAC4tvv6iiIwVkb6qui5s\nP5kIlIg8ADMf6sNNABTAFaparrDOd1T16eKT3mwRmaSqFWZrNC41fk9B88zq1LYxH5T73gD8HuYm\nvEFEzgFwLYCv1H+UuSTROYotwNXF+1EbgMdEZI2q/j3rQeUY//m1trgtXwKlqp+s8XNPF39vEJE5\nMGZk0wpUjd/TUgB2DYv+AJq876aXKr63RwBcXPFdrcNSmPPFoX9xGwnAuh8VRORxABMBUKDCWQKg\nn7Ve8fxqmPnfInKEiBxjbRoJgHGVUmYCOMBaHw/gLxmNJXeIyA+t1dHgOWTz4RzF4voEAH8q8/6W\nRUT2FJEzrE2jwHMpCDsE8SeY+xFE5CMA5pRz7wE5jEGJyAAAX4dR17NE5HZVfRJGaaeIyH4AhgG4\nS1WboBFEbYR9T6r6rohMF5GfwGTx3cgECQ9dxe9mGYB9AZyb8Xhyg+32FJGlAJ5jgkQoawAcJyI7\nAdgWwHxVvT3jMeUGETkMwGkAhhRjvj8CMAPANSJyOYyX58yK++E8KEIIIXmkYVx8hBBCWgsKFCGE\nkFxCgSKEEJJLKFCEEEJyCQWKEEJILqFAEdJAFKsW1P2zhGQBT1hCAIjIx0XkURH5h4hcISJTROSn\nIvLNrMfmICKjYSpEh70+QUSeF5EHRWRQcdt2IvInEfkFgI8U56AQ0hBwHhQhRYqtOPqo6uTi+kAA\ne+ahKHHR+vkHgH8rN/teRL4A4Ouqeqi17TpVPae4fCGANap6c9pjJiQutKAICUBEhgA4T1VnicjB\nIvKsiNwsItcVraz9i+/rV9z+bRG5VkSOLm6/WkTeKVpiD4jIj0VkNxG5T0SmicglIrJKRL4kIr8T\nkTkisq+I7F1skfLvviFNBLDIEScRGSoit4jIxSJyvTMeAHcB2FNExhTftxuAt6393A9WzyANAi0o\nQooULajjATwOoA9MU8wrrdfaVfVyEfkcgE+o6vki8n0Yi+RqEekF4FUAI4oFRNcDGAJTtXksgMsB\n/K+q/q+I7AHgIVUdISL9AMwBsA9MC5CpjhVnje2bAEap6nnF9dsB3Keqd4jIrgD+oKqOaE4H0KM4\nvqkArlPVJcXXegNYqaq90/kWCUmO3NXiIyRjHlXVySLSHcDOvtdeK/5eBrcq81gAy0VkMkxhzOcB\nDASwHMASVV1TfN9zRavm+8X1D6vwq+paEfkjgC8B6AngloBx9YSpregwFsASERkO4wlZYr12I4An\nROQ/AezgiFORLgDdRaRNVQvlvghCsoYCRUgAxY6yS8TbUTbI3TAHwGJV/TkAiMipAFaE7PYlmArq\nc+BtiQKYhpK/B/CEqv4y4LMLYIrb2sd92OkLJiILrbHPFZEXAdwG4DrffgbCuAopTiT3UKAIASAi\nB8C0n+5RrL4sMJ1SNxVv/ocB2FdEZsJUaR5bjPtcDeCHxc/0BPCeqqqInAmgv4hcoKo/LR7mOwBm\nFFsNLIQleKr6uoi8C+CBkCE+UPy8wyUArixaZdsD8Ff2vwHAVao607d9PEycipDcwxgUIXVCREYC\nWKqqa0RkFwA3q+rRItKuqptF5FonxhTy+S/DZBkGWVhRjj8AxqI6W1VX1/RHEFJHaEERUj92BfBd\nEXkBxsV3UXH7j4oJFY+X+7Cq3iIi+5Z7TwW2A3C6qm6MsQ9C6gYtKEIIIbmE86AIIYTkEgoUIYSQ\nXEKBIoQQkksoUIQQQnIJBYoQQkguoUARQgjJJf8fJZ7U2zeFr94AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "{'name': 'Bi2Se3.dos.1'}" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "plot_dos(header.db['Bi2Se3', 'dos', 1], dostype = 'total', up =1) #plotting total DOS; files saved in DOS folder" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "['Bi', 'Bi', 'Se', 'Se', 'Se']" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "header.calc['Bi2Se3', 'dos', 1].end.get_elements() # just showing list of elements in cell" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "-- ------Start plot_dos()----- \n", "\n", "-- Bi2Se3.dos.1 e_fermi 4.2271 \n", "\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAARgAAAGoCAYAAACdRPr5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXec3GW1/99nZne2peymJ5QAIREMggZECPxgKQrKtaGg\ngoCgoGCDC4RiIQGRzhUQEAtcUCxchRsvSJHICoQAhkDoEEklbZNstpdpz++P75TvzE6f+U4979cr\nrzzf/szszGfOOc95ziPGGBRFUZzAVeoOKIpSvajAKIriGCowiqI4hgqMoiiOoQKjKIpjqMAoiuIY\njgmMiEwVkV+JyItpzjtSRHwi8kGn+qIoSmlw0oI5DFic6gQRmQycDGxwsB+KopQIxwTGGPMg0Jfs\nuIgIcDVwOSBO9UNRlNJRyhjMpcCvjDE9JeyDoigOUleKh4pIAzAXCIjI0cB44Osi8rAx5qkE5+t8\nBkUpA4wxWXkbxbBgIh0SkWYRmWSMGTHGfNUYc70x5jqgB/hNInEJY4ypqH9XXHFFyftQK/3Wz0fx\n3udscXIU6QjgNGCaiFweslq+BlxpO6dORH4AjAPOEZF9nOqPoijFxzEXyRjzNPB03O474s7xYwV6\nr3aqH4qilA5NtHOI9vb2UnchJyq135VGrbzPkqtvVUxExFRCP5XSICI5xwiUzAm9z2UX5HWMPfbY\nAxGpin977LFHqd9ORSk4FW3BVNMvVzW9lmKj711xqDkLRlGU8kYFRlEUx1CBURTFMVRgisT3vvc9\n2trauO+++0rdFUUpGiWZi1SL3Hrrrbz22mul7oaiFBW1YBRFcQy1YBxk+fLlnHvuuTQ1NXHQQQfp\nUKpSc6gF4xA+n48TTzyRiy++mKeffpozzjiDF154odTdUpSiUtUWzEEHFfZ+y5dnfu6yZcvYtm0b\nJ510EgAHHHAAc+bMKWyHFKXMUQvGITZv3kxraytWZVCLCRMmlLBHilJ8VGAcYvr06XR3dxMMBiP7\nduzYUcIeKUrxqWoXKRuXptAceuihTJkyhT/96U985Stf4ZVXXuGtt94qXYcUpQToZEcHeemll/jm\nN79JU1MTc+fO5c0332T79u1cf/31/Md//EfMueX+WsoZfe+KQy6THVVgyoRqei3FRt+74qCzqRVF\nKStUYBRFcQwVGEVRHEMFRlEUx1CBURTFMVRgFEVxDBUYRVEcQwVGURTHUIFRFMUxVGAURXEMFRhF\nURxDBcYhrrrqKqZPn87ZZ5/N5z//eQ4++GA+/elP09XVVequKUrR0MmODnLmmWeydOlSXnnlFZqb\nm/nmN79Jf38/999//6hzy/21lDP63hWHXCY7VnU9mIN+WdiamcvPyb7AzAknnEBzczMAp512Gsce\neyy/+93vYirdKUq1oi6Sw7S1tUXaEydOxOfzsX379hL2SFGKhwqMw9hjLtu2baO+vp5JkyaVsEeK\nUjyq2kXKxaUpNE888QQDAwO0tLTw29/+lpNOOkndI6Vm0CCvg5x55pmMGTOGLVu2sG7dOqZMmcJ9\n992XcHWBcn8t5Yy+d8VBg7xlyOTJk7nttttK3Q1FKQkag1EUxTHUgnGIq666iscff5zGxkZ22203\nzjzzzFJ3SVGKjsZgyoRqei3FRt+74qCrCiiKUlaowCiK4hgqMIqiOIYKjKIojqECoyiKY6jAKIri\nGI7mwYjIVOAnwAHGmIMTHD8D+BiwGpgH3GaMWeZknxRFKR5OJ9odBiwGDkhyfAZwvjHGKyIHA78G\n9ne4T4qiFAlHBcYY86CIHJni+DW2TTfQ52R/iokxhm9/+9u8/vrruFwuZs+eza233kpTU1Opu6Yo\nRaOcYjDfAy4odScKxWOPPca6det4+umn6ejooKuri23btpW6W4pSVMpiLpKI3ADcbYx5Mdk5Cxcu\njLTb29tpb29Pe9+DDipwyczlmdeXaWtr47XXXuPJJ5/kmGOO4Q9/+AMej6eg/VEUJ+no6KCjoyOv\nezg+FynkIt1ojPloaLsZaDbGbBcRF/BfwJ+NMc+IyInGmAcT3COnuUilFBiwrJibbrqJd999l3PO\nOYfLL788abEpnU+TO/reFYeyqwcjIkcApwHTRORy4Cbga8B+wHnADcBngQ+J9c3bCxglMJVIb28v\nRx55JMcffzxr1qzhuOOOY9ddd+WMM84oddcUpWg4HeR9Gng6bvcdtuMXAhc69fxsLY5C8tBDD9Hd\n3c33v/999txzT3bddVcCgUDJ+qMopUDLNTjEu+++y4UXXsjQ0BB9fX0ccMAB3HHHHdTVJdb0cn4t\n5Y6+d8UhFxdJBaZMqKbXUmz0vSsOWg9GUZSyQgVGURTHUIFRFMUxVGAURXEMFRhFURxDBUZRFMdQ\ngVEUxTHKYrJjrsycObNqFpKfOXNmqbugKAWnohPtSo3X62X+/PmR7f3335+77767hD2qTTTRrjho\nol2Rif9Q61wjRYlFBSYPgsFgym1FqXVUYPIgXlDUglGUWFRg8kAtGEVJjQpMHsTHYFRgFCUWFZg8\nUAtGUVKjApMHKjCKkhoVmDxQF0lRUqMCkwdqwShKalRg8iBeUDSbVFFiUYHJA7VgFCU1KjB5oDEY\nRUmNCkweqAWjKKlRgckDjcEoSmpUYPJAXSRFSY0KTB7ET25UgVGUWFRg8kAtGEVJjQpMHmiQV1FS\nowKTB/EWjAZ5FSUWFZg80IJTipIaFZg80GFqRUmNCkweaJBXUVKjApMH8S6RWjCKEosKTB4ksmBU\nZBQligpMHiRyiVRgFCWKCkweJBITFRhFiaICkweJxEQDvYoSRQUmD1RgFCU1KjB5oAKjKKlRgckD\njcFUD0ET5MbnbuSM/z2DlVtWlro7VYMKTIFRC6YyWbp+KX98/Y+80fkGd710V6m7UzWowOSBukjV\nw+udr0faL258sYQ9qS5UYPJABaZ6aPG0lLoLVYkKTB6owFQPjXWNMdtBo3/HQqACkwca5K0efAFf\nzLY34C1RT6qLOidvLiJTgZ8ABxhjDk5wXICfAn3ATOBuY8wLTvapkKgFUz2MBEZit/0jo6waJXsc\nFRjgMGAxcECS4ycDY40xl4lIG/C8iOxjKsQMUIGpHuItmHjBUXLDURfJGPMglnWSjBOAZaFzdwJD\nwFwn+1RIVGCqh3hBGfYPl6gn1YXTFkw6phArQH2hfRVBIoF56qmnAPjiF7+Ix+MpdpeUHPEH/THb\nGoMpDKUWmE5grG17XGjfKBYuXBhpt7e3097e7mS/cubmm28GYOPGjVx88cUl7o2SKfGjRiN+dZE6\nOjro6OjI6x7FEhiJNESagWZjzHbgEeD/AfeLyASgAXgj0Q3sAlMupAoV/elPf1KBqSDiBSbeoqlF\n4n/IFy1alPU9HI3BiMgRwGnANBG5XEQagK8BV4ZOeQDoE5EfA9cBp1dKgBfSD0lX0EupeQLB2PKn\nKjCFwVELxhjzNPB03O47bMcNcJmTfXCSdALS3d1NW1tbkXqj5IMh9m+pAlMYNNEuD9IJzM6dO4vU\nEyVf1IJxBhUYB+nq6ip1F5QMUQvGGVRg8iCdBdPXlyoFSCkn1IJxBhWYPEiXVKdLyVYO8aNIvqAv\nyZlKNqjA5EE6C8bv11/BSkGHqZ1BBcZB1IKpHFRgnEEFJg/UgqkeAkZjME6gApMH6QRGLZjKIf5v\nqQJTGFRg8kAtmOpBLRhnUIHJAxWY6kEtGGdQgXEQFZjKQS0YZ1CByQO1YKqHUXkwAc2DKQQqMHmg\nQd7qQYepnUEFJg/Ugqke4gVGly0pDCoweaAWTPUQLyjxMRklN1Rg8kAtmOohXmC0WFhhUIHJA7Vg\nqge1YJxBBSYP1IKpHjQG4wwqMHmgAlM95CMw7/e+z13L7+LNbW8WulsVjwqMg6iLVDnEu0TZCMzl\nSy7nVyt+xfmPna/5M3GowOSBWjDVQ/zfMlOBMcZELJeuoS429m0seN8qGRWYPFCBqR7iLZj4EprJ\niF9iVkefYlGByQMVmOohVwumZ6QnZnvIP1SwPlUDKjB5oMPU1UOuMZj4NazjLZpaRwXGQdSCqRxy\nHUWKn7M05FMLxo4KTB6oBVM95JpoN0pg1EWKQQUmDzQGUz2oBeMMKjB5kG5dJBWYyiF+1ChngVEL\nJgYVGAdRF6lyKJQFo0HeWFRg8kBdpOoh11GkeIEZ9A0WrE/VgApMHiQSmDlz5kTaasFUDvEuUqaJ\ndhqDSY0KTB7YBebUU0/lmWee4aqrrorsUwumchhVD4bMMnLjhUhjMLGowOSBXWBEhKamJurq6iL7\nVGAqh1HD1DlaMBqDiUUFpkCICAButzuyTwWmcogXilxjMJkKU62gApMH8RYMoBZMhVKoUSSthBeL\nCkweJAry2i0YDfJWDqMEhtwERpc7iUUFJg/UgqkOjDE5WzC5lnmoFVRg8iCRwKgFU3kkEpNcXSS1\nYGJRgckDu8C4XNZbqRZM5ZEobqIxmMJQl+4EEdkX+AbwQaAJWA/8xRiz2OG+lT2JYjD19fWRttfr\nHXVcKT/UgnGOlAIjIicDJwGPA08CPmAC0C4iJxhjznG+i+VLshiMiFh+fTBIMBiMWDdKeZIobpJr\nHowKTCxJBUZEXADGmJMSHH5ARPYXkbnGmDcc612Zk0hgRIT6+vqI9eL1emlsbCxJ/5TMSOTW5JrJ\nqwITS1KBMcYERSTpN8MY86ozXapMwgID4PF4VGAqiETuUK4WjMZgYklnu39PRG4WkUOL0psKI5EF\nA5bAhNE4TPmTSGAytWDURUpNuiDv2cBbwJdE5GvAa8BvjTE9Ka+qEZKVa1CBqSwSiULOFozmwcSQ\nUmCMMS+HmvcC94rIXOAuERkxxpzheO/KnEwsGJ9PV/ordwppwaiLFEvaYWoAEXEDnwG+BRwC/DHD\n644BTgS2Ahhjrow7vgdwI/Ai8BHg98aY/8us66UnmcDoUHVlkU8MJl5Q1EWKJd0w9a1AL3AWsAn4\nFfAFY0x/uhuLSBPwC2BfY4xfRP4sIkcZY56ynbYAeMYYc4uIfBh4AKgYgbFjF5iGhoZIWwWm/Ekk\nJjlbMOoixZAuyPs1YBLwaWPMQcaYuzIRlxCHAmuNMeG/wFLghLhztgCTQ+0pwPIM710WqAVTHSRy\nazQPpjCkc5HONcbcn+O9pwB9tu1eLDfIzn8BD4nITcBHgauoIDQGUx1oJq9zpAvy3g8gIhOAHwJ+\noAN4zxjzTpp7dwLjbNvjQvvs/DfwK2PMn0RkErBKRPY0xnTH32zhwoWRdnt7O+3t7Wke7zz2ZUuS\nCczIyEhR+6RkTyJrRQUGOjo66OjoyOseGQV5gWuxXJx9geeBK4HvpLlmGbC7iNQbY3zAYcDtItIG\n+I0xfcCuwObQ+d1AgCRum11gyhG1YCqXQlow1TSKFP9DvmjRoqzvkanAvGOMuVdELjHGdInIhnQX\nGGOGRORc4DYR6QRWGmOeEpHrgB3A9cB/At8XkfnAXsDlxpiurF9FidAYTHWQyOrIVWAyva5WyFRg\n5orIdMCIyHgsMUiLMWYJsCRu3yW29lIsy6gi0US76sAXHG1lak3ewpCpwNwD/AtrJvV5wJcd61EF\noVMFqgNvYPTfKB8LxhgT83moZdLlwXzEGPOyMeYZYFcRmWSM2V6kvpU9KjDVgS9QOAsGrBwaQQUG\n0ufBnCUi/xSRhSLyARWX5GiQt3JJZMH4g376RvoSnB1LPrVkaoGUAmOM+S5wFFacZIGI/ENEFojI\n7kXpXZmjFkx1kCgGA/Dx336cf6z5R8pr/Sb3AHEtkLbUmjEmaIz5uzHm68BxwNvAtSJyluO9K3M0\nyFsdJHKRwLJiFvx9QcprE87ErqKh6nzJdLLj3lhZuT6sEaTLjDHrnOxYJaDD1NVBIhcpU/Ip9VAL\nZFos9kdYBb9vBmYAVzjWowpCpwpUB8lcpEzIJ4emFshUYF4GNgL7GGMWYLlJig17YW91kSqLZC5S\nJqiLlJpMBWZ/4FbgiVAZhowS7aod+1wkOyowlYW6SM6RqcBci2W1XINVcOrvjvWogtAYTHWQzkVK\nFswHdZHSkVGQ1xjzbmiFgfnAO1juUs2TTGC04FRlkc6CGQmM0FiXeGUIdZFSk5EFIyIXA7cApwOz\nsSwaxYaWa6hcBn2DKY8P+4eTHlMXKTWZukhjjDFHAW+GSl6OqtdSiySzYJqbmyPtoaGhovZJyZ4h\nX+q/Uarj6iKlJlOBcYf+D3+jxjjQl4ojmW/e0tISaff3Z1phVCkVA76BlMdTuVDqIqUm09nUARF5\nDGgWkYOBFQ72qWLIxIIZHExtfiulJ50FMxJI7OYaY/KqhlcLZBrkvUJEPoE1XL3SGKOjSCQXGLsF\nMzCQ+tdRKT3pYjDJ8mSSWSoag4mSaZD3ZGPME8aYG4GtInKDw/2qOJJZMAMDAymHOZXSM+hPLTDv\n975Pz/DoxUyT1d9VFylKpjGYfcKN0KL3WuyC1Hkw4ZGkYDCoI0llTjoL5gf/+AGf+v2nWN+zPmZ/\nMoFRFylKSoERke+LyBrgAhFZLSJrRGQVmcduqppUlondTdI4THmTLgYDMOIf4YE3HojZl9SCURcp\nQrp6MLcYY/YEzjHG7GWM2dMYM9sYc36R+lfW2KcK2OciQaybpCNJ5U26UaQwG3tj80vVRUpPRi6S\nMeZ/7Nsi8nlnulO5xNdgbWtri7T/+c9/Frs7SoYETZARv+XCigjnffS8pOe2NbXFbKuLlJ5Mg7zH\niMiLYTcJ+LXD/aoIksVgACZMmBBp33LLLUknRiqlxR5/aapr4pQPncL3PvY9XDL6q1Hnio0MqMCk\nJ9Mg71ewqtndBcwBdBSJ1DGYxsbYuSu9vb1Od0fJAfvo0LiGcTTWNXL6AafzhX2/MOrc+OHqZAJT\nTas75kumAvOOMWYnUBdapbEt3QW1QCoLprW1NWa7p2f0MKdSerqHo7Ne7C5Qvbt+1LnxGb1qwaQn\nU4E5UkSOBhpF5NdYJRtqnlQC8+EPfzhmWwWmPLELTGtD9Eehwd0w6tz4sg4qMOnJdLj5S0AQa13q\nb2At+6rYiBeYo48+OmZbJz2WJzuGdkTarY1RgfG4PaPOVRcpezK1YE4wxgwZYwaBDuBs57pUOaSy\nYOrq6jjiiCMi25psV55s6Ikusz5j7IxIO6HAqAWTNblm8iqkDvJCbOGp4eHkNUWU0mG3YKaNmRZp\nJxIYjcFkT66ZvKMjYDVIKgsGYkeSVGDKkx2DUYGZ2Dwx0s7HgtFM3igpYzDGmFuAW0TkpPhkO0UF\nphro80aXhx3fMD7SzicGo5m8UTJetkREporIBBE5X0RmOtqrCkQFpjKxC8zYhrGRtrpIhUEXXsuD\ndBaMxmDKH/sC92M80UKN6iIVhkyHqe0Lr31NRFIv2Fsj2NP/1YKpTGIsGE9qCybsIr2z/R0eeOMB\ndg7vTHhPdZGiZCowuvBaAtSCqWy8AW9koqPb5Y5ZmiSZi+QP+rng8QvoHOhMel91kaJku/DatejC\nawlJJDD2Bdj8fk2+Kjf6vdEyGmM9Y2OXnkniIq3rXpdSXEBdJDtJLRgRcQF7GmPeM8a8C7wbOvRU\n6HgjMNkYsyHZPaqddBaMCkx5Y4+/2AO8kNxF2tK/Je191YKJktSCMcYEga+LyFdFJDJ+JyJuETkU\n+A01vj5SukQ7u8D4fLkvsK44Q7L4CyR3kTIRGI3BREnnIv0QmAm8LiJ9IrITS1S+CywwxvSlvLrK\nSWfB1NVFDUQVmPIj2QgSJLZCgibI9sHtae+rLlKUdIl2QeBq4GoRGQc0ANuNlskHsnORVGDKD3sM\nJl5gWjwt8acD0DXUlfa+6iJFyTTIizGm1xizTcUlMYkExr5GtQpM+ZHKRZoxdgZf2PcLNNc3x+xP\nNjRtR12kKBkLjDIaDfJWNjGjSHFBXoDL/t9ldHytI6YQlX3ukh17zEZdpCgqMHmQzpjTGEx5k8pF\nCuMSF/Wu6A9FMgvGnkOjLlKUrARGROaKyIfTn1kbaAymsokZpvaMtmDC2MtnJovB2AVGXaQo6co1\n/LeI/DLUPgV4FrhXRH5QjM6VOyowlU2yiY7x2C0Yu9VjJ0Zg1EWKkM6CGTbGnBNqXwh83hhzALC7\ns92qPFRgKo9MXCRInBMTT0NddFqIukhR0gnMJgARmQFMM8Z02PfXOhrkrWxS5cHYSbTCQDzqIiUm\n3WTH3UVkDFZi3WLb/smZ3FxEjgFOBLYCGGOuTHDOdwED7AmMN8Z8I5N7lwMa5K1sUg1T27G7SMmw\nX68uUpR0AvMwsAbYDHwSQESeBBKP1dkIzbr+BbCvMcYvIn8WkaOMMU/ZzvkqsNMY87vQ9n65vYzS\noDGYyqaQLpJdYNRFipIuk/d/gf+N23dshvc+FFhrjAn7BkuBEwhNlgxxKvBoyIqZRoUtSasCU9mk\ny4MJk5EFY7teXaQoaevBiEgz8AVgb+B94I8ZzkGaAtjP6wU+EnfOTGCcMeYnIjIbeExE9qnEbGEV\nmMoiEAxE1qUWkVEZu3YyicGoBZOYdMPUM7Esj6MBDzAfWCoiu2Vw705gnG17XGifnV7gBQBjzKrQ\nOZncuyzQIG/lMuAbiLRb6lsSLnYfJhMXaXJLNCypMZgo6SyYC4HjjDERYRCRqcCC0LFULMMKEteH\n1rM+DLhdRNoAf8gKepJQdbzQZEoXkHA+/MKFCyPt9vZ22tvb0zzeebKdTW2MSXieUnxS1YKJJxMX\naVxD9Le0Wlykjo4OOjo68rpHOoHZYRcXAGPMVhHpTXdjY8yQiJwL3CYincBKY8xTInIdVpD4+tC/\n60TkUmAWcIYxxpvofnaBKRfS1eR1uVy43W4CgQDGGILBIG63u5hdVJKQ6QgSJLZg2pra2DlkTRuY\nv9v8GAuoWlyk+B/yRYsWZX2PdAKTLBaSUYzEGLMEWBK37xJbuxc4N5N7lTvJLJP6+noCAesXzefz\nqcCUCZmOIEFigTl6j6Opc9Xx3s73OP+Q81nbvTZyTF2kKOkEplVEZhhjIol1IrIL0JrimpohnYsE\nlsCEC377fL6YlQaU0pFpkh3EZumGaa5v5vuHfD+yvb5nfaRdLS5SIUgnMDdgjey8g5ULMwOYTSgn\nptbJZLBLR5LKk3xdJHvmLlCVLlIhSJcHs1lEDsHKxp0JvAg8aIwZKkbnyp1MLBjN5i1P8nWR4gXG\nLVHXVwUmSibrIo0A27BWdtwW2lbIXmB0qLp8yDTJDkaLSaJ9bldUYJKt+FiLpBSY0GoC/8ASlfVY\nVsy1InK0MaamVxSIRwWmssi0Fgyoi5QP6SyYS4EfGGMeC+8QkU8AlwGXJL2qRlALpnJRF6k4pCvX\nIHZxATDGPAHoN4XMgrwqMOVJpsWmABrco0eR1EXKjHQCM5hk/0CS/TWFWjCVSzbD1Ooi5U46F+m4\nUD2YeA4BfupAfyoKFZjKpdcbTUa3p/knIlEeTPw+dZESk05gvCS2VnS8NQ4VmMrCXrx7QtOElOdm\n4iLZLRhNtIuSTmAWGGP+Fb9TRA50qD8VhVowlUnQBOkejg6CtjamTkzPKMhri8HoVIEoKWMwicQl\ntP8lZ7pTWWiQtzLp9/ZHRKDF05K2HENCF8mtLlIm6MJreaAWTGUSngUN0NbYluJMi2yDvOoiRVGB\nyQMVmMrEvjpjuvgLZO8iqQUTRQWmQKjAVA6dA9ESR5lYMImmCqQaRdIYTBQVmDxQC6Yyea/rvUh7\nZuvMtOfHC4zH7RlVYlNdpMSowOSBBnkrk1VdqyLt2RNmpz2/qa4pZjtR/V51kRKjApMHWq6hMtk2\nuC3S3m18+hrz8X9bb2B0VVd1kRKjApMHdoFxuRK/lWrBlB89wz2R9viG8RldYxeZRBaKukiJUYHJ\nA43BVCY9IzaBacxMYNLlyqiLlBgVmDzQGEzl4Q/6GfBas19c4ko70TGMkHq5GXWREqMCkweZFv0O\nowJTenpHopMcxzaMTbngmp10Vom6SIlRgckDdZEqj1ziLwCfmPWJSHv+bvNHHVcXKTGZ1ORVMkAF\npjKwWzCZxl8ALjjkAt7c9iZ93j6+c/B3Rh2PsWDURYqgApMHasFUHjEB3iwsmPGN43ngpAcImmBC\ntyq+4JQuE2yhLlIeaJC38rCXaUhXaCoRyWI2LnHFCIrJbPHTqkcFJg/Ugqk87C5Sujow2aJu0mhU\nYPJABabysAd5c7FgUqE1YUajAlMgVGAqg1xjMJkQU9VOh6oBFZi8UAum8hjyRVc9bvG0FPTe6iKN\nRgUmDzTIW3nYJyrWu+pTnJk96iKNRgUmD3Q2deVhF5hEtXbzId5FUpFRgckLdZEqDyctGLuLdM/L\n93D43Ydz3bPXFfQZlYYKTIFQgakMHLVgbC7SH17/A96Al/95839iSnTWGioweaAWTOXhDToYg7G5\nSHbW7FxT0OdUEioweZBJkFdnU5cXvkA0DlZoCyZZlq99FclaQwUmD9SCqTxGAiORtpOjSHbsuTe1\nhgpMHqjAVB5OWjB1rsRzh+3Zw7WGCkyBUIGpDJy0YJIJjH2CZa2hApMHwWA0z0EFpjIohQXT5+0r\n6HMqCRWYPFAXqfIohQXT7+0v6HMqCRWYPFCBqSyMMTEWTLqVArJFBWY0KjAFQgWm/LGn77vElTRv\nJVeSCcyAb6Cgz6kkVGDyQC2YysKexVto6wXUgkmECkwe6GzqzAiaYFmMpKjAFB8VmDzQ2dTp8QV8\nnLn4TI6971juW3lfSftSKoEZ8A5k9GNUjTgqMCJyjIjcLiI/FpEfpzjvVBEJikizk/0pNOoipee5\nDc/xRucbANz6wq0lLWFQKoEJmiBD/qGEx6odxwRGRJqAXwDfN8ZcCewvIkclOG8fYF+o7DLsyQTG\n7bYVIQoGa+6XbHP/5pjtrf1bS9ST0gkM1K6b5KQFcyiw1hgT/tleCpxgPyEkQhcDCyHN4r9lSCYW\njIjEiEytWTHhdaDDbOjdUKKelFZg4t+HWsFJgZkC2FMYe0P77FwNXGkToYoiU2uklt2k+CzWnUM7\nS9ST0gpMrbpITq7s2AnY14UYF9oHgIjsCrQCX5Loz/9/isjfjDErHOxXwcjEggGrZMPIiJVBWmsC\nE+8alHKK9lNsAAAgAElEQVQ0yckkO0gtMMP+4YI/rxJwUmCWAbuLSL0xxgccBtwuIm2A3xjzPnBW\n+GQRuQa42RgzmOhmCxcujLTb29tpb293sOuZkanA1LQFMxJnwQyXzoKxTxNQgUlPR0cHHR0ded3D\nMYExxgyJyLnAbSLSCaw0xjwlItcBO4DrAURkEvBNrCDvAhG5yxizOf5+doEpR1RgYMXmFfzv2//L\nsXsdyxEzjwDUggljXy6lUoj/IV+0aFHW93DSgsEYswRYErfvkrjt7VixmKud7IsTqAUTJRAMcNET\nF9E70suTq5/kkVMeoa2pjX5frMDYl24tNmrBFB9NtMsDDfJG2di3MSIe3oCXl7e8DIx2kUpZfKmk\nFkyNBnlVYHIkXlxq3YJZ2702Zvv93veB0S5Sr7c2LZhKdJEKgQpMjqjAxLJjcEfMdjjWEi8wtWrB\nqIuk5IzLlfptrAWBiY+tDPuH8Qa8MbknUNoC2CW1YNRFUrIhk3KZYWpBYOKFY9g/nDA9fsA7gD9Y\nmvfAaQumsa4x6TG1YJSsyFVgqnVGdXwwN5nAJDq3WDhtwTTXJ5+rq3ORlJxRFylzCybRucXCaQtG\nBWY0KjA5oi5SLPExmH5vf9IvValyYZws+A0qMIlQgckRFZgob257k3d2vBOzr3u4O6krVKqRpGJb\nMBOaJkTaKjBKztSyi/TAGw9w+kOnjxKTrqGusrNgnJ5NHS8wu4zbJdJWgVGyIhsLpr4+ao5Xm8Bc\nv/T6hPt7RnqSxlpKNR+p2AIzfcz0SFsFRskKdZGs+UfJMMawrntdZNueI9I11OVov5JhF5h6t/Mx\nmBljZ0Ta/d7+mqtmCCowBaFWXaQt/VtSHl/dvTrSnjVhVqRdKoHxBW3LxroLu2wsJI7BhC2lWq3L\nqwKTI2rBROcbJeO9rvci7VltUYHZMbQj0emO47QF45LYr9MYzxjGNoyNbNeim6QCkyMqMOmFYtAX\nrR02d/LcSLscXCQnYjAAp3zoFABaPC2079HOGM+YyLFaFBhH68HUCrXqImUTrJ03fV6kXc0C852D\nv8OB0w9k7wl7M65hnApMqTtQqagFk7nAHDHzCGa2zoxsdw11YYxJ+74VmmIIjMft4cg9joxsj6mv\nbYFRFylHVGCSrxBgD3ZOaJrAovZFeNweWjwtAPiD/lGrDRQDe5DXKYGJp9YtGBWYAlCrLlKiAt7f\nmPcNbjn+FuZOmcu86fN44KQHIoHOiU0TI+dtH9xetH6GGfE7O1UgEbUuMOoi5YhaMLEp/7uN342D\nZxzMqR86lbENY7n3c/eOOn9C0wTW96wHLDdpr7a9itZXUAumFKjA5IgKTOys6OuPvZ7ZE2enPN8+\nN6cUgd5ixGDiqXWBURepANSqi2QXmPGN49Oeb3eR4ktsFgOnJzsmQgVGyQm7BZOOahQYY0yMizS+\nIb3AlNKCMcY4nmiXCE20U3LCLjC1aMEM+YcipS8b6hpoqEufej+x2WbBFDmb1x5/cbvco7JunUIt\nGCVv0glMNc6m/uPrf4y0M7FeINZFKrYFY3ePnJiHlAwVGCUnatlFWtu9ljv+dUdke1LzpIyuK6WL\nVAr3CFRgVGByxD71vtZcpJVbVsZsf2DiBzK6zu4iFTsPphQjSKACowKTI7UsMJ0DnTHbX/nQVzK6\nzm7B7BzeWdT6KDEWTJGS7EAFRgUmR2rZRbKn+Z/30fMyTphrrGuMTCPwBXxFnS5QiiQ7GC0wtVZ0\nSgUmR2rZgrHX37UHbjOhVHGYUrlIda66yAhbLRadUoHJkWwExj6K5PV6U5xZGdiLdtvzPDKh1gQG\nattNUoHJkWxcpIaG6LDoyMhIijMrA7trM9aTncCUKptXBaY0qMDkSDYWTLUJjN2CGdcwLqtry8GC\nKWaQF1RglBzIVWCGhyt/EfQYCyYPFylRuQeniEm0yyDruJC01LdE2kM+jcEoGVDLLpL9S5JqudRE\ntDW1RdqZWjCrV6/mW9/6FpdccknOAl1KC6aprinSttcprgVUYHKkll2kYX/0S27/8mRCjAWTpCJe\nPDfffDPLly9nyZIlLFmyJKvnhbGvS11sC8YuwiowSkZkIzCNjY2RdqWPIgVNMGINiEjWAdNcYjDP\nP/98pP30009n9bwwdqursa4xxZmFRwVGyZpsXCSPJ/olrPQYTPwXNdvC3a2NrZF2MWMw9vyTbN26\nfLE/T/NglIzINtEufE4wGKzoZDv7FyRb9wjyD/JmI+x27MKYS7/zoak++rwB70BRn11qVGByJBuB\ngVg3qZLjMPb4Sy6uxriGcZFaLH0jfTGjO4mIF5RAIPl62KmIEcb64gpMzCiSWjBKJmT7S1otgd58\nYxkucWXlJg0NxX4hc3UxSxmDsQuaxmCUjMjWgqmWXJhCWAL2oep0i7fFC0z8dqbk69rlgwZ5lazJ\nR2Aq2YLJZ4g6zITGzEeSCiUwMf0usoukAqNkjbpIuQuM3YLZ2r819fMKJDD2L7ZaMMVDBSZHNMib\nuyWw94S9I+2XNr+U+nlx7mTOLpKvdEFeHaZWssYuMJnkglSNBePPP1h64PQDI+13d7yb+nmFCvKW\nMAZjf54OUysZkc3KjgDNzdFfsf7+yp1RW4jRmFkTZkXaa7vXRpY/Sfi8AglMKWMwLZ7aHaZ2dOlY\nETkGOBHYCmCMuTLu+AJgKrAFOAj4sTHmnUI9v7u7m/vuu4/GxkbOOOOMGCsiX7J1kcaMiU7ZHxio\n3F+xQlgCYzxjmNIyhc6BTvxBP5v6NrH7+N0TPy9OYMKJivYqgZlQLhZMrcVgHBMYEWkCfgHsa4zx\ni8ifReQoY8xTttNajDEXhs4/GbgB+Eyh+nDjjTfy2GOPAdYclh/84Afss88+Wae3JyJbF6mlJfor\nVskWTKEsgWljpkWKh28f3J6xwIBlxdgFOxPKJgbjG8IYU5DPYCXgpIt0KLDWGBO2f5cCJ9hPMMZc\nEdeXglaBfuGFFyLtt99+m9NOO43TTz+dNWvWANDT08Ppp5/OggUL2LRpU1b3zlZg7F+IShaYQqXc\nT26eHGmnWsIkkUuUbQzLGBNjORQ70c7tckcmhQZNMGZmd7XjpMBMIVYwekP7RiEiHuAM4IeFenhv\nby87d47OEn3rrbf45je/yTvvvMMxxxzDm2++yT/+8Q/uvPPOrO6fzdKxUD0Ck+9UgTD2xdq2DWxL\nel4iCyZbgfEFfQSN9fdyu9xFrwcDtTtU7WQMphOw11McF9oXg4jUA3cAlxlj1iS72cKFCyPt9vZ2\n2tvbUz5869bY/Irjjz+eJ554gmAwSFdXF6eeemrM8UcffZSrrroq5T3t1IIF4w14R5VjKNScnskt\nUQtm22B2ApNtoDc+ObAU7klzfXMka3nINwTF9dJyoqOjg46Ojrzu4aTALAN2F5F6Y4wPOAy4XUTa\nAL8xpk9EmoGfAzcYY94SkRONMQ8mupldYDJh27boh/bggw/mJz/5CSeeeCLnnXdeQWYz5yMwfX3F\nWw8oF4ImyAWPXcDzG5/n+x/7Pqd86JTIsUJZMJnWhSmEBVPK+Eui5w74KiPIH/9DvmjRoqzv4ZiL\nZIwZAs4FbhORK4GVoQDvpaH9AL8F5mMJz1PAJYV6vl1gJk2yzPF58+Zx1VVXJRxNcrlcWWXnZusi\nTZ4c/cV+7733Mn5OKXhx44ss3bCUQDDAzctujrgXUDiBaWuMZvOmqmxXEIEp4QhSmFqty+voMLUx\nZgmwJG7fJbb2F5x6dmdn1BubMiUa+vn4xz/OvHnzuOuuu3j44YcjFeaCwSD9/f2MG5dZlfxsLZi5\nc+fi8Xjwer2sW7eOzs7OmH6VE+u618Vsb+jZwMzWmUABBcY2XSDVjOpCBHlLOZM6TK0OVVdtop3d\ngon/Ik+cOJHLL7+cZ599lhkzZkT2d3enntlrJ1uB8Xg8HHDAAZHtZ599NuNnFZsdQ7HrFa3vWR9p\nO+EipRKYRDlDlWjB5BPkDQahUlecrVqBsYtFW1tbwnNcLhetrdHaJD09PRnfP1sXCYjxZ5988smM\nn1Vs4oeN7esgFWI2NcSVzhzamXTN5kQCk22QN59VEApFrgLz73/DZz4Dn/wkrF7tRM+cpSYEZvz4\n8UnPsx/LRmCytWAAjjrqqEh7+fLlo0a6yoX4FRd7RqLvS6EsmKa6psgIlTfgTZpCn2jELWuBKWE1\nuzC5Tnj8059gyxbYvh3uuMOJnjlL1QqMXSzsVko8doFx0kUCy1WbN28eYFlAd999d8bPKybbhzKz\nYPIRGBHJaAmTRAJTkTGYHKvaLVsWbec5YlwSakJgUlkwxXSRAL7+9a9H2osXL2b79uRZrKUi3oLp\nHenl2fXP8pk/fCamAl2+6wtlEuitlhiMfRQpm+Vj419qjjXPS0ZVCowxJmOBKaaLBFZOzv777w+A\n3+/nL3/5S8bXFoOgCY7KS+ke7uanz/yUTX3R6RT17vq8v6z2oepEuTDBYLAgAlPKmdRhJjZPjLTD\nc7DS4fdD/EeyDH+PUlKVAjM0NITPZ1Wr93g8McWe4rFbMF1dmS/GnqvAiAhf+cpXItuLFy/OeSkO\nJ+ge7o7JewF4e/vbo74UbY1teWfE2gUmUW3eHTt2JHxvshWYUlazCzNtzLRIe0v/loyu2bJltMWS\nYPZLWVOVAmOPpaSKv4A1ZB3miSeeyPgZubpIYI0mhS2nzs5Oli9fntX1TpJo4qF9mDqM/Rc5V+wu\n0kubXho1kpQsCJ7PKFKpLJjpY6ZH2is2r2BDz4a01ySaf5uFkV0WVKXA9PZGg5Kp3COIFZhsLIlc\nLRiA+vp6jj/++Mj2r3/967JZUjbVzGY7dusjV+z3eGTVI/zk6Z/EHE8mMPnEYEo1TL3b+N1i5l/9\n9Z2/pr3m/fdH7yvzWSajqHmBmTNnTqTt9/uT5mPEk4/AAHzuc5+LtFesWMFFF11UFis+Ziow9jyW\nXLFbMACL31kc457ZBcbu5mZrwdjLVJZKYFzi4sR9Toxsr+palfYatWDKFLuLlC71v6mpKfLh9fv9\nDA5mNoSYj4sEMHv2bM4555zI9nPPPcef//znrO9TaOJHkJJRiOHeRFaQPdi7ZUs0VrH77tGCVJVo\nwQAcvMvBkXa69aAgscDYfjsrgqoUGLsFky4GE39Oprkw+VowAGeffTaf/exnI9v2AlmlIhMLxiUu\nvvjBL+b9rERV7OzLmNgtmJkzZ0ba2QpMOVgwAOMbo9a0PbcoGRs3jt6nAlMG2IebM5m8mMtQdSEE\nRkT40pe+FNlet25dirOLQ7qF0OZMnMPPjv8ZcybOSXleJsxsncmnZn8qZp99hMU+YXW33XaLtPOx\nYOz5KMVmXEP0s2jPjk7E+vXwxhuj96vAlAGZ5sAkOidTgcl2VYFk2L84GzduLPmQdZ83dRTx/EPO\nZ/5u8wv2vEXtizhhdrSS6taBqNViFxi7BZN1DMZWf6VUo0gA4xuin7O+kb5R6QB2rr462g66hvHX\nWZa1CkwZYM+OTTbR0U6+FkwuMZgwTU1NESsrEAhkNV3BCexZpi4Z/brs6f2FQETYs23PyHbYRQoG\ngzEz4vOyYHzlYcG4Xe6IixY0wZRTBl4KrUfn82xj9YEns+rQY+ma8QcVmHLAHhycNm1aijMt7G5U\nMV2kMOGCWEDJpw7YBWbXcbuOOl5ogQGY0hItpxEun9nV1RWx5lpbW2MqAmYdgykTCwYym1Vtf3m9\nUx7DNd6K9m7f/R56eiurbkNVCow9OJiJwNitnESFwhNRKBcJYnNxyklg9mjdI+aYiMSY+YXCLjDh\njOH4gmG5Lr1rjCkbCwZGL2GSiA22HLzhMW8TLlkUqO9iR4UlwlSdwAQCgZTFphJh/4Lv2JHZMG2h\nXCSItWAyfb5TpBKYiU0TcbvcBX9mJgJjL3OaTQxmJDASiXV43B5H+p8NdoFJVpv3zTej7Ukzt+Kx\n1V3f6c1sHlO5UHUC09nZGbEuJkyYgMfjSXNF/gJTSAumlALjD/ojEwNd4mKvtr1ijk8dM9WR59rX\nSNo2uA1jzCiBsS9cNzAwkHFCpN0NKeUQdaI+JLNgXn892m4c34P996sv2Eloml1FUHUCs3bt2kjb\nnpyVilzyYArpIk2YYKuwn8WEy0Jjzxdp8bQwb/q8mOPzps2Lv6QgNNU3MbZhLAC+gI/u4e6EFkz4\nx8Lr9WbsJpWzwCSzYOw14aXZ+jy6Q4aX37O9oqYLVI3A9Pb2ctZZZ/Hd7343sm+vvfZKcUUUewwm\n0y94IV2kXGvSFBq7ezTGM4YZY2dw5ofPxOP2MGfiHE7d/9QUV+eH3YrpHOgcJTAiktNoXzkLTDIL\nJhwGNBhGsF5nWGAC9d0VNZLk6KoCxeThhx/m1Vdfjdk3a9asjK6dNm0aIoIxho0bNxIIBHC7U/vq\nhXSRcq1JU2jiBQbg2wd/m28f/G3Hnz2lZQqrd1pFZxMJDFijfeH4Wk9PD1OnpnfZyllgklkwYYEJ\nuvtxuUMrUoYtmPqdlDiTISuqxoJ55ZVXYrY9Hk/a1R/D2HNR4otVJaOQLpJdYHpL+PMUIzD12S0u\nny/xQ9WJBCYXV9YuMKUeQYL0FkwwGJ0xHajvjghLXV143062JV8Is+yoGoGxD+/Onz+fu+++O6Nf\nuDDZxkEK6SLZ83BKmWiXyIIpFnaB2dq/NUZgwovW2UcE7blOqYgpNlXiHBiIFZj3Owc4/3w4//zo\nxEZ7CeKG8dEfOrvAdFbQQFLVuEj2X/7zzz8/4/hLmLa2NtassZbGzlZg8rVg7L/MZWPBlFBg1mza\nwJo11mjWHns0R5Ls7DlNuQhMuVkwDz08hD+0PJbfDz//eexUAM/4nYRD2fX1ofMqTGCqxoKxr/c8\nduzYrK/PNtBaSBfJ3t/e3t6SzUcqpcDYg7xL/rme7m7o7obe3qmR9zcXgbG/pnKIwYRLdo6MwIYt\n0RjM889bBabsAlM3JvpDF7VgulVgSkG+AmPPs8ikJkwhXaT6+nqam0NzVEJL2JaCcrBg/H5Y9V60\nEMrWrVHLxi4wmzdvzui+fSPRz4V9NnOpaPFYn7OeHgi6Yz9nzz8fW7HONSaaExW2YAL1O8nwpZcF\nVSEwXq83UnKyvr4+4eL26YhP5EpHIV0kiI3D9JUo0aGUAjNtzDT6+2HVKvD6owJTVzeVsN6HYzGQ\neUKifXZ4ONemVBgDb7/WzLZtsGPHaIH517+sQt9h3DYLJpwvGnQP8eqbI9x+O5RJldWUVIXA2OMW\nY8eOzekLH7YgIDOBKaSLBLECU6o4TCkFZlzDeLo2WoH24NAgRsI/GFMjc3NyGUWyF3Ya6ymtwNx5\nJ/zqjubI0iPxArN8OaxcGd2uHxfrIkVzYXZyzz1w771O9zh/qkJg8nWPoLQuEsT2uxYtmK1bwXRZ\neUtm0E/QbYU36+unRiq7xecLZTJdoNQu0hNPwGc/C4ccAnffDa5A9Ics6B7koosgnOfZ0wN/tdUC\nb5oUG2wJz/f011uJMkuWONr1gqACE0JdpNIKzKuvQsNgaORv0IdxWQJTVzclUl2/oaEhMqs6EAhk\n9ENQShfpnnvg8sut0pfheu52gZm8ywAnnQQHHjj62l13heH6qKs4Y+wMwlPWAiGBee89yLCEdMlQ\ngQmRrYsUCAQi7XRZv5kQP5JUCvp9UYEJByOLxcqV0DBgWTD1AT9NY8MWzDRWvb+D+1bex2tbX8va\nTSqVi/Tqq3D77bH7XC447phm9tkH5syBveYM4XbDQQeNPu+Sy3yR2jgiwgcnf5CWFthrL/A3WJaN\nMVZpzXKmKvJgSuEi2QWmEC6SWjDQMDgLYwxm0Mew8RJsqMe0BPlV51l4/m8jrWM9zGyMzjzv6elh\nl112SXlfuwVTTBfJvlD97Nlw6aXW/353Mx33WvGU8FSBI4+EG26A8EfqlFNg9307MSstK3ly8+TI\nzPaGBpi033r6QsHg9ethn32K9aqypyosmPggby5k6yLZg7yFsGBqOci7YQO89RY0DOwFIwEkaBCB\nQPNO1hz6ZfpkIzt2wOp1XjYPRguCpbNgjDH0DEdzmorpIr39drR99tlwwAHQ3Byb7DfkG8IYw+TJ\ncNllMGECHH44nHMOvLczOqV6xtgZzBwfrUnM+KjZsiH9ApElpSosGHveSCarCCRCXaTSCExvr/UF\nBHAHxjJmeBwDgAjImPqYc42BrUODTA39LqZ7n/q9/XgD1mhUc31zURPt7MPNtnrluF1uPG4P3oCX\noAkyEhihsa6Rz30ObGvxsXJLdDhp7uS5MUu8DDasibQTLW1STqgFE6LULlKpR5H8QX9k8p1LXEX7\nMi5ZAvYqoQdN3DfSdo8bXSys3xV9b9JlXIdjGACTmielOLOwGGONioWJL6qYSV3e1zujVaf2n7o/\ne7btGanG18M6AnXWZ77ck+6qQmBKMYpUaAum1DEYezB0jGdMwhUFnMBeHvLEE+Gz+8yjtamVcQ3j\naJkwenKit85LwFjvfTqBsa+xlInA7NwJ775rCUQ+9PVFC3c3N8OYOGMwpmSDd/RnLWiCvLX9rcj2\nflP2o7GukX0nWeJb74ENcy8gUNejFkwxKNQoUviDVYoYTKldpFLli9hjFcceCxvXb2T6mOnsMm4X\nLvrkxbQOzcPlj35DA243Xr/l9qQTmFU7ous/x9cXjue99+CEE6wA6003Zf867Njdo6lTLXfPjn2E\nLpEFs2bnmsj+Sc2TItMojt7zaAA89TA0biVb9r6OrVujweFypOoEJtcYzP894uGtfw+xZm2A4eGR\nGAslEdU2imS3YIolMD4f/Pvf0e05c4L87W9/i2z/x1GH8/o1v+SJUzpo3fEJa6fHzaDXMg/SCfEb\n26JLI+4zKfVQy+LF0dT7P/4R8llkM5V7BLG1dhKt8Pjylpcj7f2m7BfJs/ryfl/miJlHIALuOuid\n/HeG67aUdX2YqhAY+wdtTLw9mgE7Brs4+4VDGBm7nm73u/T0mrRWTLUFeWMS0oqUL7J6NZEC1jNm\nwD333BI5VldXxz777ENbGxx6KOw2dg/rgMdF/4j1657KggkEA7y48cXI9gFTD0jZl/hlwZ94IvPX\nEY9dYBKVJLKXprC7cWGWb1oeaR80I5ok43F7uPm4m5k3fR6eegBD38R/RmrJlCNVITD5WjA/eOxq\ny8yst96OnqGBtIFeJ4ep+/r6Mq6aXyjsw7n2RdoLwWOPwU9/Cg8/HBvfeCsaZmD8+Be5//77I9uf\n+tSnYlaEmDnOWv1RGuro8/bQNdSVsm7Pyq0rI6NiU8dMHbVCgp0VK2ILbUNsbChb0gnM9LHTI+3X\ntr426m/9WudrkfZHpn1k1PXHzTouMru6f8IzKjBOYx+mzjYG8+BbD/KP1f8EQOqst2PYP1x0C6au\nro6mJiuoGQwGM4oDFRKnLJg//xl++EN48EFYuBAuvDCafRoWmJ07H2Lp0u9Erpk/fz6XXnppzH0+\nMvWjuALNSFMdxsC2gW1s2Zq8JszS9Usj7cN2OyzpdI433oCLLordZ4yPv//9DyxcuGhUnedMSCcw\ndrF76O2H+MIDX4hYW1v6t0SWz22oa2DvCXuPun7+bvOpD2nv4PiXWb2ufKdVV3weTHz9lGxcJG/A\ny10v3RWdzxGyYLyB7ASmEDEYsMp2bgwNC+zYsSMndy9XCh2DGR62Ju7ddBMYE6C39wmGht7gd79b\nw/331zNr1il4vQfT3/8cW7ZcTXiFmZaWFq644opR61nN3q2V6Y8u5P19LsZgjbRs3LoRX5JFgp7d\n8Gykffjuh8cce/RR+P3vrRorb78djb2MGwc9PUE2bryMvr4Obr8dHnroUa67biHHHXdcxnPO7AWh\nEsVgjpx5JBObJ7Jj0Co5sb5nPd/523dYcNgCnngv6pvtN3m/hAvFTR87nV3H7saO7RswrhFeWPMa\nkGBCUxlQ8QLT398fcVdaWlqysiYeXfUoOwZ3RARG6t0YIMAI27enFphCu0hg1TsJC8y2bduYac/Q\ncpiYOTtZZLxu2ACPP25NuuvtHaCnpxe/fyKvvOKht9fPwMCLdHbewshIrA+ycuWztLZ+nsFBK6AZ\nMt645pprYhaiCzNjBozbcTQTtp5EZ+M14Pcz5B+KWcUzzJb+LbzXZT2v3l3PR2d8NHLs1VfhRz8a\n/TpaW+Gaa/o5//xF9PV1ANZrWr3az0UX/ZD777+fa6+9Nu3UBIgVmEQrF7d4WrjzhDu5a/ldLHt/\nGYO+QYImyLXPXhtz3glzTkj6jEN3/ygr11hpvK92vYgKjEPY15K2r2+UjrD1Mjxs5SxMWv91Nnmu\nBMC4vLz17naOPjr59YV2kSB2Cdlir1EdE4PJYP3pvj7LCrjzzjfZtOlOhoffIRCIxkTc7vEEg0MY\n48XthlmzrBnFO3bAwIAVi+nufgiwLIe6OheLFy9m+vTpCZ+3xx7W/029+2Ma6jB+KzGwM0H9yGfW\nPRNpHzj9wEix70AArrlm9L2nToVrrx3kppu+w9DQ64wbFy1dGQzCmjWwc+dbfOMbF7B48e9Srhaa\nLskuzF5te3Hdx69jx+AOvvPod2KG1AHOOfAcPj3n00mf84kPfoy7nn4QY2CzvEh//7mj8m3KgZoU\nmLXda/n1il+zqbeTTZvA7ZvAxPfPoNP9c7xYn6xX3nk35T2ccJHsFdsS/TI7SddQVBwmNE1IcSb8\n9a9rOeecH9Pf/y7G+BOeEwj0UF9v1TppbYUxY5r4/Oc/z7Rp03j00Sd56aVXCQatL+SYMfDFL34x\nqbiAZQnMmAHergOguR5/9zCDrsGEcRi7m3HEzCMi7T/+Ed55Z5Censfw+bbQ3r4Lhx12GPvtt4mf\n/ewWXn/9dVwu2GUXOO20k1m69AzWr/8dXV1/ors7yDPPrOarX/07EyeewPTp8NWvWrObvV5Yuxam\nT7cEyZ5k15JmUvrE5oncecKdnPvIuazasQqXuLj08Es5cd8TU153yO4H4fEIIyOGobFv8Orb/cw/\nqKHUwUUAAAvrSURBVPwUxlGBEZFjgBOBrQDGmCvjjjcANwIbgb2B64wxq+Lvk4pMBKZ7uJuXNr3E\nS5tfYvmm5azeuRqfzzLvR0Zg6rpv0OhuZvL4Xdi4Zi0A7254L+G9wqRzkTo6OjJelylMKS2YHUNW\nPKDvnT4mft5yUUZGrFm/dXXWL/obb8DLL6/j8su/js8XtXgaGy2RqKtzMX78BAYGuqivD9LYCFOn\nTmX+/PmcffbZkWVHvvzlL/PXv/6V6667Dp/PR2NjI2eccUbK/onAN74Bi66cgbthHL5AH/5AkDfX\nxg73bOzdGMkjcYmLY/c6FoCVKwNcddUjbN58B37/dqZMsUaKEo0WLViwgJNPPpkVK+Diiy/kvffa\n2LbtDoyBv//9v9lrr+MBFw8/LBx8sJX929VlzXQ+5JDofWbMGJ1kF8b++WhtbOWez97Dsg3LmDVh\nVsy8o2SMbxzPLg0fYPXI2yBBHn1lBfMPOiLtdcXGMYERkSbgF8C+xhi/iPxZRI4yxjxlO+18YJ0x\n5kYR2Q/4DZDVu2QXmNbWVoZ8Q2zu38y67nWs2LyC5zcs551tq/D5rJyLkREYGrL+AbR0zWfi1i9w\nxVVw0+I92LjCGn14f1vqTKtwDeD+kQEuf/QGRp71UBccy5aufoaGg3QvW8nMT+5HXdMQLs8QMydO\n5+MfOJzDdpvPpOZJNNY10lDXQL2rPhI8LKUFs33QErSet/t4+vGJXP/kCCtW7KC+fjtjxnSyZs3L\n9PU9j9cbfV/cbsu9OPLIj7BgwcXstdde1NXVEQwG2b59O263O2E8xeVy8bnPfY5Zs2bx6KOP8qlP\nfSqjNaw+/WlYulT499LZ+NiIzwf/eiOaM2KM4a6X7rLaQcPclrks+8c7PP74Gn7/+78wNGT1vaHB\nmrmciPPOO4+TTz4ZgHnzrED1Y499kbPO+m/8/kG83jW8/fbHADd1dW2sXr0rHs8MoI76+mk88sh+\neDy7AnVMnTrMm28OU1dXR2trK62trRH3Kv4HqLGukaP2PGpUf7xeL93d3bz//vv8+9//xu/3s+ee\nezJr1iw+OuOjrO61UqGfWvU8WX51ioI4lW8hIkcDlxljPh7avgDYxRhzke2cp0PnLA1tdwO7GmP6\n4+5l4vv5278u4VsLzsK3sw8IYsTQMHcq7nnjMADGMr8TrQAiQQ9jug5n3LbjGN91FFdd6eL44+GC\nG27g1hsXRs7zjJ9g3QhrnWAiLUOwbxgTDGAI4PnknrjaGmOe4V22Cc+hM+LeE6uYkEj0l01w4Qo2\n4jaNBLeO0Pv38LCoUDfeZl+baA/sySTR98VE/wufa2z7MaH3ZPTf2+DHYLl8ge4RxpiDCQYTr5sc\nxu2GH/3oPD772cOZPXt2Qar6ZcKmTfCxry5i0xvXWzsEglsHaZg1lWDQT9D4IGgwviCe/t1wBWJ9\nFLcbPvzhSXzmM8eyatUqVqxYgcfj4ROf+ASnnXZa0vW0rrjiLm699Vf4fFbMyOdLX01ul12sc+14\nPB7q6+tZv349c+bMoa6uDrfbjTGGQCCA3+8nGAwSCATwer2MhP2tBAwFA7y9czXisj5Q9dJi/R1S\n/ily+zsd1f4F/vbLn2GMyeoGTrpIUwB7znsvEJ81FH9OX2hf2nU7egcGGN5pcyMM+MdDwpCAcdPU\nN5fm7oNo6TmQpt79qaOJD38Yzr0GPhLq1Sf/32HcemP0Mm9P+gXYgMjwdjqMSTRvJAgMAoOYCUEC\nbh8ELBEI9CT/cDmF+N0EJbG4NDRAfX0du+zyAS699AxOPjlFFNwhZsyA0z/5Ja79900wEojopy9u\neoXLP26UuDQ3t3Dhhafxve+dEinP4fV6cblc1NWl/ipcdtkZ7Nz5Js89F82vGRiwRMbjgbFjrVnh\nXV3W37mtbbS4hJ8XFo5MV0ZIRqPLhXsYAkFrqN6Lc/kw27pyW4zJaQvmcmPMsaHtRBbMP0PnhC2Y\nntA5oywYRzqpKEpWlJMFswzYXUTqjTE+4DDgdhFpA/zGmD7gEeBQYKmIfAh4JV5cIPsXpShKeeCY\nBQORUaSTgE7AZ4y5SkSuA3YYY64XkUbgBmALMAv4qTHm38nvqChKJeGowCiKUtuUdaKdWEMTZwNX\nAUcZY960HVsDhIuTbjTGnFaCLo4iTZ9PxQp0+4HVxphflqaXyRGRK4AjbbuuNsaU5RJf6fKsyhUR\nWQaEI+mB8EhrOSEiU4GfAAcYYw4O7cs6b62sBQbYH3gBSDQx6J4y/UAl7LOI7AJcZIz5SGj7RRFZ\nYoxJndFXfIwxpvjDQ1mSYZ5VufJomX527RwGLAbshXSyzlsra4ExxqyEiFUQzxEichEwDusPtqyo\nnUtCij4fByy3bS8DPgn8vEhdyxQRkcuBEcAN3GaMSZ0UUxoOBdaa6FyFpcAJQCUIzP4icjHQDPzL\nGPO3dBcUG2PMgyJyZNzuE4DLQsdfF5H9RWRMooGZMCUXGBF5DCv3JbILK7vhx8aYh1NceqkxZnno\nl2yFiJxgjFntZF8jHcytz4nygpJMhXOWVP0HHsD64g6JyLnAbcA3it/LtGSSZ1WuXBv67LqAp0Wk\n1xjzbNqrSk/WeWslFxhjzPE5Xrc89P+QiLyCZdIVRWBy7HMn1khZmHFAVvOuCkUW/f8HcFHas0pD\nJ9Z7GGZcaF/ZY/vsBkXkGeAooBIEZitgr+WR9j2vyIp2InK0iBxn27U3UG6xjHgeJ7Zox6HAoyXq\nS1JE5Hrb5hzK932N5FmFtg/Dyqsqa0TkAyJylm3XbMr3PYbYuQXhvDVS5a3ZKbkFkwoRaQW+jaWU\n54jI740xL2Kp5hUi8hFgF+AvxpjnStjVCMn6bIzZKCI3ish/YY0i/aoMA7wAvlAftwH7AeeVuD8J\nsbtwItIJrKyQAG8v8CkRmQ6MB9YbY35f4j6NQkSOAE4DpoVicjcBtwI3iMgPsKzxr6e9j+bBKIri\nFBXpIimKUhmowCiK4hgqMIqiOIYKjKIojqECoyiKY6jAKGVJKMu16NcqhUX/EDWGiHxURJ4SkaUi\n8mMRuUJEfiYi3yt138KIyBysmbzJjh8mIq+KyBMiMiW0r01EHhGR24EPhXI1lBKjeTA1SKgkQ4sx\nZkFoewLwgXKYMBqyPpYCH0+VJSoiXwG+bYw53LbvTmPMuaH2BUCvMeY3TvdZSY5aMDWOiEwDvmuM\nWSYih4jIyyLyGxG5M2TlzAudNza0/xIRuU1EPhHaf62IrAtZQo+JyM0isoeI/FVEFonIxSKyU0S+\nKiJ/EpFXRGQ/Edk3VLLiM3FdOgrYHBYXEZkhIveKyEUi8otwf4C/AB8Qkbmh8/YA1tru83+UaRZy\nTWGM0X819g+4AvgXcDNwF9YscPuxq0PtLwK3hNo/xZrBDtCI9WV2hbYHsaZGCFb9kAeAk0LHZmEV\n1wJrotx7QANWqYLrE/Tte1glIsLbvwe+HGrPBFbYjt1o699CYKrtWBMwVOr3utb/lfVcJMVRnjLG\nLBCROmDXuGPhdXO3EZ09uz+wXUQWYAnJq8AEYDuw1RgTWs2ZlSGr4qeh7cgMd2NMn4g8DHwVS2Tu\nTdCvBqy5WmH2B7aKyO5YFrdt5Wd+BTwnIj8EJhlj7Md8QJ2IuIwxCVbHUoqBCkyNY6xqcFvjqsEl\nCsy9AmwxxvwcQEROAZIt7PMG1kzsV4gtUQFWga0HgOeMMXckuHYD1iRL+3OXmFCdHRF539b3d0Tk\ndeC3wJ1x95mA5WqpuJQQFZgaQ0QOxCpzWB+aJStY7spI6Mt7BLCfiDyONZt2/1Dc41rg+tA1DcAm\nY4wRka8D40TkfGPMz0KPuRS4NTSl/32wL0ppVonIRuCxJF18LHR9mIuBK0NW0UQgftb8XcBPjDGP\nx+0/FCtOo5QQHUVSCo6I7A10GmN6RWQ34DfGmE+IiMcY4xWR24wx301x/RlYo1yJLJxMnt+KZdF8\nyxjTk9OLUAqCWjCKE8wEfiQir2G5SP8Z2n+TiAwCz6S62Bhzr1hFpXOlDTjTGDOcxz2UAqAWjKIo\njqF5MIqiOIYKjKIojqECoyiKY6jAKIriGCowiqI4hgqMoiiO8f8BIL+0oU1OMpsAAAAASUVORK5C\nYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "{'name': 'Bi2Se3.dos.1'}" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "plot_dos(header.calc['Bi2Se3', 'dos', 1], iatom = 1, dostype = 'partial', orbitals = ['s', 'p', 'd']) #plotting for atom 1" ] } ], "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.6.8" } }, "nbformat": 4, "nbformat_minor": 1 }