Documentation/KLayout-SiEPIC-Tools Python Programming.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "224a47cb",
   "metadata": {},
   "source": [
    "# Creating layouts by programming KLayout in Python\n",
    "This tutorial is based on an open source application for creating layouts (KLayout), and open source set of helper scripts in Python (SiEPIC-Tools), and an open source Process Design Kit (PDK) with technology files (SiEPIC-EBeam-PDK).\n",
    "\n",
    "Note this notebook is for educational purposes. Copy and paste the necessary code segments into your own Python project, using your favourite development environment, e.g., VSCode, Spyder, Jupyter notebook, KLayout's embedded Python, etc.\n",
    "\n",
    "### Software installation:\n",
    "\n",
    "***KLayout***\n",
    "\n",
    "> pip install klayout\n",
    "\n",
    "***SiEPIC-Tools***\n",
    "\n",
    "> pip install SiEPIC\n",
    "\n",
    "***SiEPIC-EBeam PDK***\n",
    "\n",
    "> pip install siepic_ebeam_pdk\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "7a1042dd",
   "metadata": {},
   "outputs": [],
   "source": [
    "!pip install klayout --upgrade\n",
    "!pip install SiEPIC  --upgrade\n",
    "!pip install siepic_ebeam_pdk  --upgrade"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f5ca9186",
   "metadata": {},
   "source": [
    "# Tutorial\n",
    "\n",
    "### Create a new layout\n",
    " - technology: EBeam\n",
    " - top cell name: Top\n",
    " - floor plan dimensions: 605,000 x 410,000 nm = 605 x 410 µm\n",
    " - database units (1 unit = 0.001 micron = 1 nm, configured in the EBeam technology)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "3ce45d0b",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "KLayout SiEPIC-Tools version 0.5.11\n",
      "SiEPIC-EBeam-PDK Python module: siepic_ebeam_pdk, KLayout technology: EBeam\n",
      "No libraries associated to EBeam technology\n",
      "Loaded technology libraries: ['Basic', 'EBeam', 'EBeam-ANT', 'EBeam-Dream', 'EBeam-SiN', 'EBeam_Beta']\n"
     ]
    }
   ],
   "source": [
    "import pya # KLayout Python API\n",
    "import SiEPIC  # import module for SiEPIC-Tools, helper functions for KLayout\n",
    "import siepic_ebeam_pdk  # import module for the SiEPIC-EBeam-PDK technology\n",
    "from SiEPIC.utils.layout import new_layout, floorplan\n",
    "\n",
    "tech_name, top_cell_name = 'EBeam', 'Top'\n",
    "topcell, ly = new_layout(tech_name, top_cell_name, GUI=True, overwrite = True)\n",
    "floorplan(topcell, 605e3, 410e3)\n",
    "dbu = ly.dbu"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1c264a9f",
   "metadata": {},
   "source": [
    "### Create a blank cell in the layout\n",
    "- name: cell\n",
    "- Instantiate it as a subcell of the top cell at position (0,0) database units\n",
    "- record the instance in a variable (inst_cell) so we can manipulate or query it later if necessary"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "fb55cd8f",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pya  # import module for KLayout\n",
    "\n",
    "subcell = ly.create_cell('cell')\n",
    "t = pya.Trans(pya.Trans.R0, 0,0)\n",
    "inst_cell = topcell.insert(pya.CellInstArray(subcell.cell_index(), t))"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1cf4ae2d",
   "metadata": {},
   "source": [
    "### Place cells from the PDK\n",
    "- use SiEPIC create_cell2 function, which is an enhanced version (error checking) of pya.Layout.create_cell\n",
    "- Instantiate it in the subcell at position (40 µm, 15 µm) or (40e3, 15e3) database units\n",
    "- record the instance in a variable (inst_gc1) so we can manipulate or query it later if necessary"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "1fa6c5cc",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAyAAAAIfCAIAAAAHSS/zAAAgAElEQVR4nO2d65HjOLBmoY5rWI0V15ySvNm1oseM9Yb7gz1sFl7EG5nAOVERo0KT+ABSJXyTmQJfx3EYAAAAAGjHr9kDAAAAAFiNcoP1+XwajgNFFFFEEcXxilNEUURxB8UXKUIAAACAtpAiBAAAAGgMBgsAAACgMb+u/OI90Zj1evyJc9V3OFHpsBWdqHTYXKj1TlQ6bEUnKh02F6r+RGqwAAAAABpDihAAAACgMWzTgCKKKKK4r+IUURRR3EGRFCEAAABAY0gRAgAAADTmf2YPoAGvf/5E8I7f33NHAgAAAGCWqcE6fn8fv79f/3wus9VbsRMooogiisuLoojiDoor1GC9/vncY1cEtAAAAGAu6g2W5a6sfzpf4LQAAABgJCvUYIW4fBVOCwAAAEaiuwYrEr66c1ZoXUVa6XVaEuaIIooooriYKIoo7qCoO0WYaLC8JxoCWgAAANAHxQar2F3dezhf4LQAAACgISvXYD1CkRYAAAD0QHcNVitCRVorzRFFFFFEUYgoiijuoLh1ijDeuSGgBQAAAEVoNVhd3dVd5XyB0wIAAIB0tq7BeoQiLQAAACiAGqwkxeKdtIoVh4EiiijurDhFFEUUd1BUmSIckx98HMP5YvpIAAAAQBqkCAuxsofYLAAAALjQlyIcGb5KmWPb1KGu+CeKKKKoXXGKKIoo7qCoL0UoIT8YgdQhAAAAkCJsDF88BAAAAGURLOHhKy8UaQEAAOyGvhosdYpZRVpK54giiigqVZwiiiKKOygSwRoNqUMAAIDl0WSwFnBXd3BaAAAAq0KR+zQohwcAAFgVTTVY5t/RgmPm+OM5PC9NCWYUUURRu+IUURRR3EFRTYpwsfxgBAJaAAAA2iFFKA5ShwAAANrBYMkFpwUAAKAUHTVYZ35QV/K1oeKPIq3qxx2mKHYFRRRRlKM4RRRFFHdQ1FGDtU8BVgoEtAAAAIRDilAfpA4BAACEg8FSDE4LAABAJr+u/OI90Zj1uveJV35wirqKE4/f3+8vcxVp1SsKn+8CJyodNhdqvROVDlvRiUqHzYWqP1FBDRYFWLkQ0AIAAJgLKcIFIXUIAAAwF+nbNNzDV2MU72hXtPZ38G7xoH2OKKKIojpRFFHcQVF6ipD8YFuIaQEAAAxAtMHCXfXjdFpcXgAAgB5Qg7Upp7UioAUAANAD6TVYKHbF2t+h+XN4vCx/VVFEUZHiFFEUUdxBkRQh/ICYFgAAQD1yDRbuai4UaQEAABRDDRb4oUgLAACgGGqwUIwppuyk1VaxByiiiKIoURRR3EFRaIqQ/KBYiGkBAAA8QooQ8uA5PAAAAI9gsKAQy2lhswAAAC4k1mCF8oO6kq/7KOYWaWmcI4oorqo4RRRFFHdQlFiDRQGWakgdAgAAkCKExlCkBQAAgMGCXlCkBQAA2yKuBiuSH9SVfEXxwirSWnKOKKKoVHGKKIoo7qAorgaLAqzlIXUIAADLQ4oQRkORFgAALI+sFGE8fKUrNojio2KP5/DEFYeBIopaFKeIoojiDoqyUoTkBzeHcngAAFgDUoQgiNNakToEAADtYLBAHBRpAQCAdmTVYKGI4p1WRVqS54giinMVp4iiiOIOioJqsCjAgkco0gIAABWQIgRNUKQFAAAqkBLBInwFBeC0AABAJtRgoahYMbFIS/UcUURxPVEUUdxBkQgWLAVFWgAAIAFqsGApKNICAAAJYLBgQdhJCwAA5iKiBisxP6gr+YqiBMWzQuv9ZXo/7tBi7auK4kqKU0RRRHEHRRE1WBRgwTCIaQEAwABIEcJeWNlDbBYAAPRgfgSL8BVMhIAWAAD04NeVX7wnGrNejz9xrvoOJyoddsGJbpEWV1jUiUqHrehEpcNWdKLSYXOh6k8kggXwA2JaAABQDzVYAD+gSAsAAOqZvE1DVviqiWIWKO6smPIQnraKNaCIoiJRFFHcQXFyipD8IGiB1CEAAKRDihAgCXaHBwCAdDBYAHngtAAA4JGZNVi5+UFdyVcUl1c8i7SuOq0BihFQRFGRKIoo7qA4swaLAixYCQJaAABwQYoQoA2kDgEA4AKDBdAYnBYAAEzeBwtFFBdWjBRpdVIcCYprKE4RRRHFHRSpwQIYBAEtAIB9mGawcFewLTgtAIDloQYLYDQUaQEALA81WCiiOE3RKtKqqdMSO0cUhStOEUURxR0U56QIyQ8CeDk9Fn8dAADaIUUIIIjTWpE6BADQDgYLQBwUaQEAaIcaLBRRlKuYXqSld44ozlWcIooiijsoUoMFoAmKtAAAVDDBYOGuACohdQgAIBxShCiiqE/RTR2uN0cUFxZFEcUdFIlgAawAMS0AAFGMNli4K4CuUKQFACABtmkAWAp20gIAkAA1WCiiuKBiw4fwJCoOAMVlRFFEcQdFUoQAW0BMCwBgJKQIAbaA3eEBAEYyNIJF+ApADpTDAwD0gxosFFHcVLGySEvFHFGUKYoiijsoEsECAGNIHQIANIUaLAAwhiItAICmYLAA4Ac4LQCAesbVYNXnB3UlX1FEUbtifDOtNeaI4hRRFFHcQXFcDRYFWADaIaYFAJAIKUIASMXKHmKzAABCDIpgEb4CWA8CWgAAIX5d+cV7ojHrdcrB5l/jthcrNumEE8Wq73Ci0mFbjX+LtF5/i7S4ULpOVDpsRScqHTYXqv5EIlgA0AxiWgAAJ9RgAUAz2OIBAOAEgwUA7cFpAcDmjNgHq1V+MF2xFSiiiGKlYnwzrR6KDdlBcYooiijuoDiiBosCLAC4IKYFADtAihAAhkL2EAB2YNyjcupBEUUUV1LMzR5qnKN8xSmiKKK4gyIpQgCQAjEtAFgGUoQAIAUexQMAy9A9gkX4CgDKIKAFAHqhBgtFFFEUqmgXab0WnON0xSmiKKK4gyIRLABQAzEtANACNVgAoAa2eAAALWCwAEAfOC0AEE7fGqy2+UFdyVcUUURxgGLDR/EkKnaFGiwUUVxGsW8NFgVYADAYYloAIAFShACwFGQPAUACGCwAWBOcFgBMhH2wUEQRxcUVrTqtAYrFUIOFIorLKFKDBQB7QUALAAbQ0WDhrgBAMjgtAOgHNVgAsCkUaQFAP6jBQhFFFHdXDG2mtdIcRYmiiOIOiqQIAQBsiGkBQCWkCAEAbMgeAkAlpAhRRBFFFIOKXR/F41Xs2r8QURRR3EGRFCEAQAbEtAAghV4GC3cFAGuD0wKACNRgAQCUQJ0WAESgBgtFFFFEsUqx+FE8xYptEXthUURRtSIpQgCAlhDQAgBDihAAoC2kDgHAYLAAADqB0wLYGU01WK9X4Y9AlRC6EswooohiCumbaVGDhSKKyyh2qcHqUYD1epn7SK1f4+2hgyOnZ/V/thuTqgIAm0NMC2B5dKQIXVtzHH6v420PHRxRyer/ak93cgCwM2QPAZZHgcEKZd/aeqxKD2fwWACQD04LYFV+XfnFe6Ix67XVeOUHc0/0vj79yuld3GPOdvdEb7vViaXiVX+/P9c/3U+8+reGdKUXi+cr5ESlw1Z0otJhc6H6nXgWab2/zFWnxRVe40Slw+ZC1Z/YvgarYQFWTd1VqN1trKm7oh4LADpBTAtANQpShBdXbs5NGnrTiImlV6Hj732G0pSRXCEAQA1W9hCbBaALTQbLZHqmy40VeKyCLx4CAPTgtFYEtAB0oW8frDK3lM7dXaXM8aptb8L4q4oiiiiqUEzfTKuhaCdQRHEHRek1WMb8yNnlDtatr4rUYEVa0oWKewAAyIKYFoBkpKcIc3N8vfuJgKkCgJGwxQOAZBSkCCtrxu+PsonUqnvJnWP9Y3N0xT9RRBFFCYqV2UMt00QRRV2K0lOE1nf6eoSI3Kp2YlEAoBpiWgDTkZ4ivOi998GAHCIAwBjIHgJMp3EEq+1jnsdUjruxK5wWACwGTgtgMOURrM/n8/292h+q+4ydwXMcr5hbruH9dLY6iX+CJ84xMrCUFSJxXm5XWXMJHdzkPja5NbkdPt7fxzlmXcAs7oqt7m8WblgopJKCt5N+nwCRmNYOH3QoojhekQhWm20a9FKw3lSakkTiA8vyPVn9NDFYTWhyawp6i1+TXBvXI17S8P5mMcBgjYSYFkBXyr9FOIbEuqvr63spP1BJwx0OiyUiB2QNb8BcRhKazkrXZKW5zKXfzqUAYOQbrPTa9uP483N/ff8xPIC5HV0/ixM79x5WMLDF1hV3Oitdk5XmIgecFkAPWu6D1TY/eFH2rJusflrtg1XPeMViij+F43OsiU8UD6n5ijL3Pt6n0++aTHhYlpj7248p75xLdJjT2uGjFUUUW9ZgNTdY6ftgpddRhY7cdh+s9EIQ99O2bUlKXCJ0QGQY8ZGExiy2BqvJrYn00+T+Zt2CXMbf39zOx3fSG+q0AIqRniK8yHU8reJecNH7EzZl+Uz5DmP84MR/TUFOXCR9LpEj0ztpUulVT+/7C4bsIUAFagxWAXgs1YQWSCELZ9dQjTRSvosXDy6CdnBaALkoeBZhIqEvDNZ4LF3pXi2kzLFtcKL5VX10V+vdx7jHWtVdbfIJkCXaxGntcGFRRFFNDZb7a3oP8foq9sG6/1pWctS8JKVss6UxlTHDYlfjb03Z3Ie5KzmVT3I6kQB1WgBepD+LsInX4ek30IqtMoMurp1aNXYF6fDcQwAv0g1WK2+Ex6qHqgux7mrkrQmFrO4HDBsMiAKnBXBHQQ1Wq9pz9sFK4aqrcH9aSei6qpErEFk8esxxwK2pR92CGr+kut6rckQf67R2uLAooqimBiu+D1ZKD/GWyn2wXIumJVpWvDzf7/XCNVghBjiJ8bcm8eDIwLpeFjn3lxqsXIhpwYZITxFe9N5DoT6HaHnBtRm/JxYA6IXsIWxIsxRh8/DVeFwPVxYbrHFp01OEidRkpnTFeEPEpy/nUTlQwxrvVVGiZ+rw/WVC2cNO7HArUZSm2CxF2MNgNdmmwZiHTlwVU52LVFRQ3+N5LF1zN14GbNMw/uty429NfYrwUaIGOYk5OZ0sADEtWBU1KUJTmsXLPasgF6nITtXw+PWxtfFOX0jgduSteRQquCY7v6/A+z8DEv6sACqRbrAs75Lllq4jyzxWlsomHmsYuUXuY5DsscaQuNFo6Jq0vWtEPhbD++UMbi7oRfo2Dfd40ql4tqT8uM4sVzddpVUNvpYarBrkJNELdjpI9A1L3sdQIjJyTdK3k7i+2H9/ff+JD8yrUr+ThZz36mKij4rNH30ocI4oLq+oowYrFB+S2a4rmtWkRkfpNg1exZROxtRjjb81kYMftwHL2v4gEuKKtLvDK3DGcsqnqMHKgpgWqEN0ivAyKyHXktUe8T1N+o+0w2Lslit8dFcRvA80LHBXie0hy7Xw3dkENnoAdYg2WKadx3r0Pb093AKoq0SOr6mV0xFV8t9kJGV50lzpVu7qcXheB5zVD8gEpwVaKDdYn8/n+/vvG73fW/yyMu/35/X6LvBAib7HPfKcYysPl8L9qk5h8FJtxu500JaIx+pxH8ffGi8h+xL/p0u9rbvyamXNscl7VU4nZYz/2GmimOW0lM4RRdWKbWqwBoTfm9djPQplteuNXRWvuDVhDG8/9TVYNSPpunFUGU1uTXE/6fc3Lne/sP0ygHG6vlfndgJ3iGmBKKSnCC/uX9PLemZz2b5Z6f1H2hdG7OdXWeZO7HQKcOdScE1SLkjWRgwDMoMhVrq5EIfsIYhCjcEy+VZJvpBShH9m5foJ4dPJIjSXsi/c5R7j/a6faeeurva2c4H1wGmBBH5dezzcN3vIej3+xLnqi52Yy/U5NWY/kqxhXwfkfsctZS5lI0l53fXW5F4TqxPvFkTuRft8PnFX5A6poP39ZQ8ykSW3JTvh4zT++u9+Wq8/+2mxYIk9UemwIweoqcECgJGEdhrrt8PC9HY2ptoBYlowDAwWAPyga5W6/HZs1ibgtKA3og1Wp+Jxd/OqASoAKgh9l1OI+xnZfr1mAV4bnBZ0QsGzCK+f9/tjtVw/6e1xFVcxt/+ISgrji0VQRPHkXm51/Hz8n9eFPNZdNW+/5jg3ptWVKeViGt+uDRWP1s89fFTsBIrSFKVHsNzRVe6D9bj9epN9tvRuiwV7EgnYCIwtzXJXRDg2gZgWNKGBwepXgNXj2cwpz7ep18VggRbiBkKIy5nbjsfaFpwW1CA9RXgvkPr7hfPAXqC57SGu4605tuo/gq74J4raFWvclasoyhU1aXc3uWiYQvJCilCUops9TL/7WuaIYj9FBRGstnGsHvEqt50IFgjnMSojx+XIaSeUBYawFiQj3WA9UvCVQO8pVhlWPRgsEAvuqqad9RVOeCdAHNGPynn0KK+XJ1ZU4GzcyBb2CFYFd1XZnn4ArI0b2uT9AD84Snm/3+cL8/Uu7qRM8eIc/n0SBRO6TnFfuIq9QRHFrorm6339hA5IbD8V6/tJb7/PcYyu5zPn6+/nXvxKljH+nTNFdFXF+1ti1TmimI7oFOGz9H9hJ/dFQSfXayqoYD2s4lxRMSGl7ZRkQQhiWmCEpwjTaeWKGn4fEEAsEtzJAu3nr+dSSroQ7pA9BLOMwTLtvBEeC9aDuquu7XgsCIHT2pr6xOT0Gqx4yyORTnSle1FE0UvDuqt7e6Q+aZn29DqzVvVY1GAtr9ijdC+u2A8U4yxSgxVpGdMJgEyIXQ1rJ0QBWfCGWR4Mln8fLAwWLADuanz7+YIlE9LhbbMqtQZrbtlBE2/kfnMQgwXa4TuDs9pZLKEY3jyrUZmYHFaAdYRrsNgHC0UULcVOdVeu4rb7YEXaa4psqMFC8Sh9C+ma4w6K5Qbrz/kDDZZHPeyNCjo5fI4NQB0D3BXt8fZ+hcywFbyRVLOIwToqvJF1ijFJPwAywV3JaWddhFbgtDTCPlj+frw1WPf94gEEQlW7nPb4AQBZsJ+WSoqtmZwarHjLI95TzsaUqq+26EowoyhKcXDsin2wHtsLog7UYKGYiPv3vt4ctSuukyIMtRR08thOihCkQWZQZjuZHRgAbzOZYLBihim3HWAWuCux7Sx+MAzebKJQX4PV8BnPVz/3EivKrUA+8bIMgfVJW7UDDINSLVkUWzM5NVjsg4XizoqzYlfsg5Xenh5XoAYLxbaKXWNaQuYoVnGRFCH7YMG2kBlU0U7uBubCO3A8ixiso8IbeffBAlBB5ENTgqug3WpheYPp4LSGsUKK8P5rk28RkiJEUYXi/YPSfVRO6JRW7aF8RD/dcyuKfv277dcch1W7kyJEcZhivdOSP8e5iutEsEItYzoBGA+xK3XtBA9AILwtO1FlJabfDwwWbEvoM1Ggq6D93s5KBmLBabUFg+VJLGKwQDi4K73trGEgH96lTaiqwRp89b01WPXeKLLXg650L4r7KLqffdrrk1Lal6kzi6xe1GChKEox7rTWmGM/xRUiWP32wQIQiPcjT0Jshvas9umfnwDpENMqYAWD5X1R0Mn1GoMFksFdLdDOcgVK4a2bziIG66jwRuyDBYpwP90EugfaU9pZpUA1OK1Hfn0+f55YdL3IfT3+RG8nx/HjuYHF6u/3xxjzetX+SLhQba8wJzY8saYT89/zxT6fz/X8O+vgeLvbc2L79eL1z+f9ZayD4+3XM9HuKlnti833ajEVbwb+2OWfqHTYjwccv7/fX+b4/f365/P65/N6caGc1zXubLp17bRNQ6SfxHZiYNAD638Zp8dgaK9sJwAAK0FMywKDVeulvO0YLOgB7mqxdhYkWBKqDE8wWLGiq+J2DBY05/6BJcol0F7TziIEa7Pzm7zcCCy8D9ZdscxjFRssXZt8oDhSMe6u2AerR/uA5y1ayw/7YKG4qmITpyV8jhbqI1juT5NO6n8AGvIYcqddb/u2/38Pe7JPTEu3wQLYBNzVwu2bLDYAFss7LQwWgAJwVwu3L7/MAMRZ9U9Adw0WiijuoBj59Fm7Psli7Tqz85+owUJxZ8VHp6VrjkSwAKTz6K5o196+5P++AxSzxh8FBgtANKEPGoEugfbi9jWWE4DmqP7T0J0i7PT1vVlfG1xpLhHSI66dviU6RqXVZUx0V+dVHekS7vdxjK7nE0CkWyprP280KUIUUfTyx2kZTXNUHMHy7jhVv6oldjhAt0mfg/fxaos7/jLrcz/Ra7BCZ6VLJHZbgGuwJLgB2pu3q/4/dYBhKPpL0Wqw7qtX6HVlt499NvRY9Wt85HS9HquJh7auQ7oTShQKXec2AUjno0SIG6C9R7uWZQNAAvKdlkqDFVkjK11OmRd5PKxMt6FjU+qxmngXy3+HUoSP56Yc0NDr/+kEd7VZu+TVAkAmYp1W+SIwqwarkyfwruU/qkz6eCxLN0Uxvc+UgY33WFk1WN7X7q/xTuJXNT7TiFDkOkcespTF/VNj5/qkk+XrzGatE2vX7qC4j2KPv6CNarAe18LiGFLk19z2rLW/ST/WWZXjt9onxrGaXJ/Hs8pyqY/Xua3BkuNyaO/aLvN/xAF0ISSmpclgpeRfioug3cbEH28nKcN4jDPVzCWUCCv+mULEDIUOeOwk3mfcuz/ed+v4euIfE7Sv1y5hVQBYhrlOq3wRmOiuQi2hxtyes/rxrrWPp6fEmSrnUuDVckM1vUm5NSmXOn582/hl+sAewV1t2I7BAujBlBR8hcEyA0dpjiOwD9ZjS0rnoSU2MfnqWqL4MCKH3evM6s1iSj/2I1Ya2dZ0xQiJ3uXxase9Zvq7KL29SQ1Wlrtavj7p2KPObMr/cEurpEERxa6KWX9lc2qwhn0ExF1CQc7osf/63Fykh/R4SZNo3KxOmpAS5/P+6v5TPHn62H9Bez+DJccN0N6pnSAWwAB6//+MdIOVWB/zeEyKRMPcXNby3DDHtJ7BiqjnpvZCvXXKCVZetNBfvkA3QHvbdgwWwGA6/d2JNlgFa2GTuqVWuTn3mCbGK30YxZ209QrFPN7WxMt+f+29RD08VhOD5TaGDqZ9mXYMFsAs2jqtCoPVvwbLWqKu+iTj+3ZbjTdyf+6KBUONGB33yPtP5T5YBd4oVNnW0Cs8Kobw3l8T/R6f98rn7oMVeo95z/KOp74Gy/0jpz5pkzozarBQRHG64uW0lq3BSglIWP80Pdbi/fUIDDt+VnEES0LOtAnx+1t2GVO872PQy3uY9bqtwZLjcmjv3U6WEEAONX+J0g1W8WI5mNyBpeSqOnmjULjOq9vEKxSTchlD7aEJFsz38SolnpLFfYmVsOrTPrIdgwUghJUN1nE8eIXHlgGkxI1SolOP800cSWhUiSMRGxeMt6TEmVL6LPO+bd+K9xiGnFWf9mHtGCwAIUwyWANrsIxTn2QdEGmpISX5mpK0Ch3mHtN8H6zpXsGlZh+siDEN/evxdPHTVSL3xWqsrMEqcFeb1Cd17V9C+3UfRxqshStpUESxkhqrIz2CFUpj3Q+wWroOIzG/Fk/MRYSs4wuGmtgiPy545NzfdI9VM9/4WzGxkxTiVTi0L99OBAtACMumCJ/H4LiZTgYrt71sGPVmsSCnJtxgHcn3N9FjdZpvW6+Pu9q8HYMFIIRlU4QWkRTh4+Jao5gVt6gZg2mxTUNWS0oOsWwkEQpShIn3N3TYpdhvvpZ7a5UidNtDx3uvalc3cG4o0K9/t/3vfRylaz9Iqr8uKUIUURSluGyK8HkMTmyjXwQr0WNVLastlvnHThIP6B0XTKHg/qbnCiM911jbJl7fu77KjLXQ3qOdCBaAEHZPEd5/7ZoiTFmzpxusx/GkO7BOccF0yu5vfNid5tvW67vrq4RVn/Zh7RgsACFgsGItbVW6eizz9P2+lB5CL0LdPuYQy0bShOL7mz79o9F823p9wy6je7djsACEQA1WrKVe0aR9i7DeY501WL3jKAU5xOk1WJGWx3NN2qNy7r/WG6yGNVjUJ13t+9SZUYOFIopCFKnBirWMHENlTs06vVMc5TFOJuGqhnSzRuK9HZ3m2/aiGXYZ3b5dwgcswEa8/4/3p+aP8ZdRzus1eQDH8XcM99et+qzsp2v/A6gZ5zVNRfO98/rnc/z+pn3DdrcRANSh3mBJWDs7eax+nCN8veZfukcqr+dIj9W8f2mrPu0T2wFgHN//a77/t76bX5/Pn/9bul4kvjb/ZhzsfV18otXJtXaOUX88wPJG5Ve4aHiPisdh3u/PcZjz5/32HPN6jbhQiSem39/4Aa7Hus+3fqjv9+fqp9OFOldf98TXP5/3l29e4fZrFb+rPLZbQ0pv944/t32H+XrHn/K61ccpJzY/Uemwt7tQ5v8ZY05f9afR8VjZisX5SgklAmK/71ZZg1XciVt0VVxO1OR7kZV0+n6f94AB39/M6I19sPZupwYLYDRn0ZXzuuaPcR2DdQRqulM6iXw3sHgkBQjxRm29Qg09vt8XqUbv/f3NjN7YB2vvdgwWwGiuwvb7r3V/jP9jFqK4zsZK6k2sTHq9GlRfnddBTj+t6F1H1fa6NURaPRDtY9oBYAKf/+t/nU9Fkfu/Nbol3LObAhWbrKn3Op4a0tf4+Bx7eBrJ97HVfO+VbZUkrsreeqOCftLbras6QPdenzRG915fNUZ3/F/HFFEUUdSiWEVxOE1CBLvVDkaVubmrBxX7YOV2UjaSJvS4v4r2wYqEpmlfvl3CBywA1PwlqjdYPeqWsvrxnlhfg9XQG2X14z14osFqfn9VGKwjusTSvnw7BgtACMsaLG/5uVuN7vqSx1MsFetFelchOxIyOlm9lZ1SdhlD1zPS7QC897esE29voT7r42Rlnfw4PbDECnQDtDdvx2ABCGGSwTIjDNadyPPd0tvjLa5ibv9xlcfj7ee7PR2faA4i/fif71bnaeIUPIswMt+UTu73MeR9vbpZKsfP+9jDYEVW5X2e09evf7d91vMWBxusHZ4ohyKKZSz7LMJcr5DS/thnDw93b6zv/96efoka9j8M977URLC8v3p7ro+TlXXyo0N2arj1idcAACAASURBVNi4nQgWgBBWThE2b291zKz+rfasU1rpDsObM63vJHRAZZzMet3WYElY9Wkf1o7BAhDCyinC+yp1jw0We4XIsueukW40sjjek5iEauiBEj3Wn5REWhyuCQUpwkhLSifxd06nOFlDg5WyKv+Y4xA34PnrkOFOGrZfCbthuqQIUURRlOLiKcIBOTXvuW3jPeleoUlQquCUrnVXZbQyWFl9NomTeVvyOvxviZXjNmgf1k4EC0AIK6cIe/x0Ep0yl066QmgytsceOqk0MViiVn3aR7ZjsAAksKzBAtiWiMES6AZob9uOwQIQwrI1WBY7pHtRRPHCu8rKcQOz6pP69e+2T6kzG++ulP6BoIhiP66/xBpFIlgAQnEXWjnuivZO7YSvACbS9g8QgwUgFMNODfu1464AxtPpf2wwWABCMZk7NdC+QDsGC2AYvQPGmmqwOn3PLkLDdG/Xrwq2Uhn8XcV+3Vo/8fs45qKV9HBLGAUP+E8lVJ/kHh/qxx28wPqkwe3j68ymGKyFK2lQRNEly1dtUYPlLsah5Tn9yJHUjD9ycOSwrP7j7XF110mk95/Vc0H/R+fxp7dfv+Zdn6edGjzvq5/uytseOp726e1EsAB6MKXAscp6DBto/Zo33WPVe4WUKYzxCo//1GpeYzxipfcNtUfuRcb79r/wldt++OYVd1eRfrLaPSOnvVE7BgugIeNN1Q/1qpOHDPpPtkK5x+qxZqcc0FU3YlxE9V/QVfP3W5n3vT4aQrGoWe3WIGlv0j53MQBYBiF/SuWm4/0eMfq2Oa9cj9W8Biv0T17FMTm1lOctxrvt57FaeevreYudvJTb/n7bYacy73v/mPhzNayMkvl7pNXoPb5hO/tgNW+ftSQoraRBEUWLHr5qTg3W0T+C1SMeMyuO1corRA5u1b/bnhueKWuPqHSNY7WNk4WOKfZ298+LFPdDzEl1+/T/5wZQh5B4lYsyg3X8txzW/EzB9Qo9ZuH+U/3lCpkSb0t6b+l9FvRfMBdLq7mQVzfp5/oWYfg7gF4hK7f4o8OQkPf6+IrAjqaugvZD8CIBIBP5fzIqDVbu6e46PR6vV0g8Mf06eL1CD5rcmhQLYv2aNZ3EQbrvkx4XLcVy+U/MfOqz8cW0LNcV6ccap1sBltgP7bntwpcKACEo+mMpX0kG1GDZn/WZ6+tRvXa2rcGyfvWOxKsY8VihNTt9vmVzrLE+V0WU1c+jwfK2xDkHFq9sK/a+EUKVbQXeN/Jp4q1PsmJaWe7KOuXuutzj2QerbbupfvBZMdIqaVBE0eW/WL6mOSqIYNXEFXqsnWVUeoXEONaw+bret6CH0ItIt5XvgZQ+C3x81jAK3s9eg5XulrKOv06hrmtku6L/KQcYieo/DQUG63rRZHFVarCOwODj7qHffMd431YGq3ecLHEY1usagxXLNPnmm7XqXzEwjykU6U4WaNe7hAD0QLWvuihfScakCO+v6w2WtyVCvxRhaCQpD3WJdFIw38oUYYH3vW+aYP24Ko8TTBnqfdOEMQbrcfOLrPez9Vnj/bu7/h5d71tSD+T9XsLP4/8mJUe5EOuqDtD1p+w76F73lxQhijsrPvoqXXPUEcHy/lrQQ1knTWjlFawT6w1WGfXeN12lbS5SQgSrQOX60HmOhVjznR2bof2xfY3/WQcoZtU/AQzWIHp4BbeTKQZrgErDXGSnOFnWMApU4h9AkXopr4oEV0H7/df1lhaAR1b1VRfrG6wxViBlJKEXBZ1crx8TT2sYrCMw39xOvD008b65wyhQCX0S2au1+dHu0ZXhKmi/fl17jQGw2Oc9X76SDK7BOnyPH0nsoXjt7FSDFfEKjzVY1q9e95A131bbNKSTrlhTt3TvIXGbhhrva+Gtwap0pd6PJLc+yc0MWjGtUOe57cPqky4inzmd2gfUmVm3lRosFFdVbOKrhM/RQnoEq+b/+70+Y3oE6/q1R8azh1dIGckYlVyhxPdPkzhZymC8Whk9OB9P/oCWr6qdeiyx7Zv8rzxsyz7xKhfpBuvwrUOm0RNRRtKkOirdYB23+eZeqMSRJHqFxJuVboZye3v0gl7v2/yieb1vrtD9c6pgdfdLpD2BJ6V/2nPbd157YGGuN/bm723RBut4ckWt2gcw3mBFJBK90eNIQiYmVyt0l5vfd+8BBbqV9y59/CGDFV/FPf3kt/uvlSR3skb75isQrASmyuLX5/Mxxhhjrhe5r7ue+H5/Xi9zHOb18hxztrsnetvPfoTP9/FE73xDr8+r5x5wb7+TO2zrOj+emHsfTWC+oePj87XeAJGhXv0/zrfsnlrX4XFe1uvXP5/j97d7wNVuXQe3/c9fxD+f4/e3O9+z3b2Yr38+7y/P8FLavePPbQ/N97Hdmq811Eh77/m6nYz/eJmrvsOJSoedfvDrn8/rn8/r9Tl+fx+/v71/HQKHPeLEGnc2xqg+hg36xWka4kn3lEZBIvNKjNCEessdyfU6MYKVMh63pUkcK2uEj+09oo/edn9sKe3JOX+uXl37378+ebEf7e38vz6ohvfwI9INVkHOK9Q+12M9eqP0Trwv4n229QrpuqF/yk3VVXqs+BwjF6FhrrnA+wbPynlyTn079Vg92lmcQCm8ddMRbbDK1qRIu4QarLZ1PK6TGFOXVmmwvO3xlrZ1V+6RreJhcZXc9ojBSqx2px5LZjtLFOgCX1VAueOYsg/W1V7zk84a+2A1vCCRkUSOTJR+tFy5U7iuRuJVvb+ovGiJ+2AVC7nV7qHn9NkdPn2X8LDeP2HXxT5Yxe2RtYp9sFAUpRj3VWvMsZ+i9AjW46IbP919XeYn6nn0RgWdeLsNvWhLusGK9+C9R1ktx1Pc7nFsid63kpT5ZvSW/OSc3Pb7OO8urVX/tB+Er0A8xKuaIN1gHXVr0pi1M3ckoZb6TobN1zIuNXOJeMH0i/YYx0oZyWNXNTT3vqGPv3o3YLmrw3pfzXYn2ttZukAsvDnbUruSdL0T1lJU70g6rZ0FI/G21HcybL7WHamcS+j+Zl200GRzDVbKKQU0977ej8ImbsD8V9XuMYWz3Yn2dhYwEAhvy06Uf8afickxBuso9QqVa2enGqzISB6rhR4dVYpKumJkJPcXBVc1xQs+zvfx+MSzso5PwVuDZf1ar3L/WOxUn+TNGJ4MeE6f1R6qM+un27zO7HEZowYLxWGK9b5K/hznKuqIYHl/LeihrJMm5HqFUCfWifUGq4wx3rcgR1zgscZHsBqqGKfa3VQ/8cbrBafHfhZoJ04AEuB9OAwM1iAKvIK3B+u1BINVppKVEMzKReZ6rCbeN2VUiePJ69Z9ArT5+a/eSxF1A56hCnAn2ttZ1WAuvAPHU/4ZPzhF6P6a2ENNJz1ShHGvEFJMie4cRfOtTBEmqriKuRVXWblI6zrcN/jwHuzVquExRdhE5U8/1rZY5m+7X/exTsjY7QUurUd7/DOnR/v9Plb2n7i2kSJEsa1iV18lZI5iFbeIYLVdO8so9gqRI0N2asB8x3jfmrql9DhWWZwsl65xMtdgRdxSqAfvwIhpNWkneACD4S0nAekGq3JNGrN2Zo3k+jU9GOOe7u3zGDjf+nhMihesVPFe5IZxsiy6et/7h6lx6qWsmJb3dOv4epdGu9XCUgcDwFeJQrrBOtwP99LAT7+1M4UyrxBZjx8tV9f5JnqXyOnu6+YG6+gfJ8saSehFm/6jj7dzXdf9ROtIt92NaT32Q/u9hTUPusJ7TCbln/Eja7DOF7mPyvF2mLWq9d6mIXcKj2uz1ys0+fEKlU0qq7fHlkcsxUSV5hfN630bqpzbNJjbE29+qD89IcfqzcoMWu2X4r1//8Vv2u6pbBPjokLtuSsfNVgoJuL+P9V6c9SuqCCC5b4OHfPYXrA8t6LJ+I8nj5juSCqvm7v0dtKN2JRELPeWeInqxx8feY/r9vcD19jtoeNT2r3u6scYaA+0E1eA5vCmUoRog3U09VgT3dURNUYNPeIYr2A1dur/ak95D0SwxtnJg7rtj2+/Hu8H12N5Vv2cdtxVcTsLITSEt5NGpBuso5HHmuuurgH09ohjvILbXhV3SQjI1dxHS8trdLrEkxKuYQ9v+sNjVX73DXdV2s5yCE3gjaSact8xoAbr4vygD+2fdD/G217srrrug+U9ILJ/UsosxniFrOG57ff7mB6QSzk+1MmPaqGAR2l73d5vz37oDft3f7UqouL1WN7+U+qu3ONHtivaB6t4UaQGC8Wj1FfpmuMOig0COyM9lvs6dMzjwYOpjMMlzqKHV4gckHhYbnvogDKvnJ6qa3vd0gffw/tW1mMJdFeK2ok6QBm8cxZDjcE6/vvcz/0RQtngc2dRqVI8gE6TanJDH3sYc9EG35q/Buv+PcGi7xJ6O6c91M4aCbnwnlkVTQYLALK4f3CLciGrtrNSQjq8W5an3GD93ZVq1PtDV/IVRRQlKMY/xKfUJxW3C98Hq8l6SQ3W8or9fJWcOaJ4QgQLYHFCH+ii3In2dqIREId3yIZgsADWx/1wF+VOtLezdkII3hs7g8EC2ALqsTq1s4KCC+8KOHTVYI38BtxJevK15stu9VPI6kqIUI16LrrS9p0Ur098Ue4kq13aPliPl7QAarBUK97fD6vOEcV0XsdxmDpe/3yO39+VnTyrvMw10vvr0DEp7aF/jR+f0klWD9fBZeN3j6y5Pola3sMq70vBZYdcXv98zhfW32zor5j2UPt1JY1zMWE3Qn9WsDn/M3sASVhLb8hD5LZLo9U4a65PzRhWvS8rcfz+PheDu2OQ41q0tLOgguFtAE8oMFhN1uwBa3lXb9Sqn4YeS8t9AYu7x7pa3MPkuBlp7Syrm8MbABL5VXzm5/N5PqiaK+flKl7tFrntp4qXrDlG+k/n8/k06ScyHqv9muO9PWUMxdc/8T42ZMx7VZHi48KQ4jbuimNcjjvHwe7q9c+n9+I6/p0zRVSd4nnrz5/j9/f501WxABSlKTaowTLdyrBCtUGJ63FiPZBbt1RTgxWvggqd647ZmmNWDVao8wJCuu58m0AoazDeQiI5sSJR7RRdbQjBKqhBQYrw5HIt6VXYWS6nJld1P7HAbXhNVcPcWW6yL8trXvPFG2nkShea//yEEDcjqh1rtRv4KmiC9BThnff7k+45rvhKOnev4E1meXGHVGM1psQ/LVNYkCvMVcw7oRoUI9yTHXcnccfrQj6fz2D3c81xmO45x6t9wHJLinCiYm4SsF6xLShKU1STInR/TezBmIdOQiqJhd7G5zDKgk9j5lugnn7RQC8hJyEhhjSrnUjGDnCXoRNqUoTFPObarAOuXx9PDFmxuW6jOLcYug6wCefqcv+C4c4ZQ9KCy4Ovgt5IN1hj9j4o8FghL1I54Ll7PeCxwKrK8h4z3f10bbdmzeq7GPgqGMavK794TzRmve564ll35ZKr7tYShQ44289f3fbzxWU7rD0O3PYm8y24wtf4y060roP34Me9LWS+o4ScKHnYVvXJfW+Ce92VdWK83VUpaH9/GasxsT20o4TVbpzAlbfnlNfS7qlY9WEnvl5/6qveX3/f4VxhUScqHXbkAB01WAV7H1g9XL+aaFF2Sj1W27orbw/D5pvSSaSlVdwOxGIZDpkxpybt91+JbSwD8SqYiA6DZXp6hfgpIdPTo+5qynzLOqn3gqCIuP8Q6JbS27FWS4KvAglIr8G6KNhfqkbosR6rd/xm2Hwr0TJOqMGqTLpK4I0wt5Te7paXsRIvAL4KRNFmHyy3MlQ1oTokd7+r3O1MH7lf1TF4FXP3D6tX7AqK9bg7A91rsywk74PlDjuy49F691GIaFvFlP2rtM8RRY2KbVKEpk+WUOy+UDVDqlR87MGUztc6hX2wIILXWkmOaRGyWg/iVSAcNSnCYlpVCGmpNJo135Sgl/yrB4m4W3GagOu6/mlifRWFViuBrwItSDdYorzRAI8lapzp/URiXVc7LMmj0xr83cOIyWM9Vg2+CtQh/VmErWqoy/px59i7pju071cu6eOM38dIP9720PF316UriY5iomKo/OXaX8oqfipwUfddqR77t0ZVtipPv6qriqYrptRXtVVsBYoo6qjBqtkXyiW9Bit0cKfvD+rdB+uxh1AjrE36d18e67fSRftFODrtRwMuxKtgAaSnCE9q4kauM8jVdW1B71zh3Lhdbv/jr08luj67sxZ1KzcXaYyceB2ZWxge3woh/l3jUL1UCvH9F1JIuUSRc7NO8RaEUYNv6u4CgEB0GCwzb7+luIforSuf3BziXO4bOLm/asfdjODyNPd273xDByRenMiXCu99WmYipXDK6tD1cPW3z3vdQgOo1HL7sS7O1dhESAW6/ocHIJ1yg/X5fL6/F/97OOc40mONv6pexaz4U268au4cQ+4ha9V3ow72Wf8a8/UQQ3IbQyGfFCP4en2Owz+1AsOUgntVvVcsFJuJfLMvZHE+n8/lsbxeOSIduuyuobyu+fH7+5zjdUy8mr4gAON2MuVzdcqf5PvfP6/H+CohH60o7qV4tMN8vRv2dhyHNbqCwRpT2Il7WOjEhpdw4nzdU+4nPl6Npu+jjoTeovf267X5eodex8/KbcxS8f4kTi2FkFzuuY//5L0mj/14RxW6bim3LDKwlBeP/xofvzXa3AuukeUnCHChJkVYTEE9UFZdkbR6o/F7NCxMbuwnPWN1D02ZQJ4oPp7HztOPr1dM6ad5Os8kXLeaXJsbJ/OO4S7k3tBQz+7w1k4RkgeEDZG+TcP4Wu/r+3Hm5xzv7TX9P6rXkzWeyH3sVAcm6nu2nVaylDnedxBId07eH2PM++vHGtbDcFiIuo9eUi5vvPTe/Fsol3hDXfXFtmkIbbIg/82DIooNFBtGwzqlCK8xFme7vL0lHh9/HTq3jPHzfezk8KUd3UYtKcKjOq9XdlZi4+NhcR4HnNJb1sGRcx+7KksReo/scXe8B5T1E5lCKMOrPYNGHhDgUJEibJvzimCpXLqh9q7jHDDfyn6U5hDvMYPr/6e9jSl4z0rvKuXIUKIqcmSZUPzEJqeHsEaeHnC68oPp20wkNj7m+OJl9Y/7WVhTMEukBckDAtxpttGo6bALn/W8YWOyl/OUPTPHtKcwcr7pPcS7Ld4TdQHYdlI13L5W4KsAvEivwbroun6H6qsivqFTPdZ1VYf5lcQ6sxAF8x2fRH+9VKXtUeypaFWjD1DMQlENVqi+qp9iDSiiOF5RTQTL/bWgh8eW0Ov0/ouH6p7Vab4p5z6OpH6oC6A3BBLfF7THibmEDFCT7yE26SdRyELpG+YO8SqARFoaLNN6yalMeHmzV4leocyXlI0zPpLcHowp6eQxCdh8qACgCHwVQC4KitzvpFuf67DcCqHex/emYDzxferlTA0ABoOvAihGeg3Wvb7nVDwX/sefO8V1UYlzrKy7uo/Z6sc7ndzxPF6ryPiz4oX1P53QlbZHEcXpojX1VWWKvUERxfGKolOEJvpluoLv93m7qv/mndt5QSqzyXhSasgS2+Oibe9LJBcJAMMgXgXQEOkpwitelb7vVCRHFlq/myztMvfBKrtuKQOo6T+lHQDGgK8C6IF0g2Uaeaz4Eq7FG7XqJ95ukjODbT0WAIwEXwXQFek1WCdXWMX7T4/t6fEYi4I5VnqF+z5YTTzH4/Vx55hl7FKuf2J7P3Sl7VFEsStufdWS00QRxemKjWuwTM/NgYoNR3yKTdZ7t5OsGqxH0muwKhkmVKMOALkQrwIYjIIU4UWn1bdJtzWdNJzXMIOCEwJQAb4KYBaaDBYAAKSArwKYjo4aLBRRRBFFFB8p279K3TRRRFGFoqYaLAAAcCFeBSAQUoQAACrBVwFIhhQhiiiiiKImxebPsZE5TRRR1K5IihAAQAHEqwB0QYoQAEA0p7XCVwHoon0EyxDEAgCohpAVgGp+XfnFe6Ix6/X4E+eq73Ci0mErOlHpsLlQA068SqzeX+YqseIK6z1R6bC5UPUnEsECAJgP8SqAxaAGCwBgGvgqgFXBYAEAjAZfBbA8XfbBOn5/Xx8fDYkodgJFFFFEsSGPW1iNn+MUURRR3EGxSw2WoQwLAOA/iFcBbAgpQgCALuCrAHYGgwUA0BJ8FQAYnkWIIoooothEsckjAqnBQhHFZRSpwQIAqIJH2QCASy+DZfBYAKCB4k8qUoEAEKE8RQiwBj22FFGKnEtRP5J+c2mSCgSA5aEGC0V9irlr56NiSofeY87G+z/9WX1ff58oF/8xtwXb2+Kqe0U/n4duXZXQ7BIvb+59jKiHjs9SjNyg+Fleh3SeeFcc46uowUIRxXUUj26Yr3e/zmFnrLfW9ev5wny9z5/r9f3I+8H3I90D7h26oq6idzCR0Xp7jrR4pxk6y9tJaBbuTCNy1uUNTT8k4Y7ncZyRsXn7t44MaXmFvG+PyFvl8S0HANtCihDW5HqcwPXCGw26xyGuYMZ1fPozCeJHWmGbSCAnEhdxJdyWgm5T5K6wjTf05f7TY8Dsfq75GRxyZ2Q1xudyj0idLy4J7zivg627b12Kez8mEPdy33IAsDOkCFFcUNFaYs2/D4eZ25J5LcP3ZTUXa45WRun6tWGa6VRslbq6hud+P+6SuBSvf7Ku4f0nJGTViXvd1fWv11UNOZh4ufrl4UIHWxb8+P39/vL864V1WdzGAkgRoojiMorlBuv7e3RpJ4ooFvL157+Ji71b5JTCfdVPmWMoIBT5NRLEiiuWBVTiZ1mKltOKRLCscJQVYYqM4fGqxkNHoSCccXzeNYy74ph69gl/HXt8CKCI4njFjju538MAAG2JRz68eOMx1r/eD0gvx3YlQqN9PMsKHbmjjScK72dFOnnkfqLlhxIv9WO395ZIbi73LlgfO96LY5zM4z2P7B1b6FcAgBAd98EybIUFsA29/9gb9l/sOwEA0qEGC0U1isUxg1ZzTKzdbqiYTtlDXR4nkqjYL5xzr39vssHHFYhKcVfC76NeURRR3EGRCBYALA4hKwAYT8caLACAieCrAGAifSNYhiAWAIwFXwUAEqAGC0UUUVxB8V6t1WpLBWlzXEYURRR3UCSCBQCKIV4FADLpbrAMHgsAWoOvAgDhUOQOAGrAVwGAFqjBQhFFFKUr/qmverWsr0ph7as6URRFFHdQJEUIAEIhXgUAeiFFCACywFcBwAJgsABACpGncQMA6GJEDdb9qfU16Eq+oogiiolYW1gNUExkB8UpoiiiuIPiiBosQxkWADiQCgSAhSFFCABDwVcBwA6wTQOKKKI4QrH+UTby56hRcYooiijuoEiKEAA6QrwKAPZkkMEyeCyAncBXAcDmUIMFAM3AVwEAnPy68ov3RGPW6/EnzlXf4USlw1Z0otJhhw54vf7WV72//rorLpT8E5UOW9GJSofNhao/kRQhABRCvAoAIAQpQgDIhi3XAQDijItgGYJYAMohZAUAkAj7YKGIIooPeLewWmyO2ypOEUURxR0UiWABgB/iVQAAxVCDBQA/wFcBANSDwQIAY/BVAABNGVqDdfz+vj7ExyhWgiKKyysWPyJQ0RxRlCaKIoo7KA6twTKUYQHIgHgVAEBXSBECbAS+CgBgDBgsgC1ga1AAgJGwDxaKKK6saJVYDVDsB4rLiKKI4g6K1GABLAipQACAuYw2WAaPBdANfBUAgBCowQJQD74KAEAa1GChiKJWxXt91ftrtLta9arupjhFFEUUd1AkRQigDOJVAADymWCwDB4LIB98FQCAIkgRooiidMWUR9lonyOKsxSniKKI4g6KRLAAhELICgBAL3yLEEAW+CoAgAXAYAGIAF8FALASc2qwjt/f13IyRrEMFFHsrWg9yqbYXUmeI4qSFaeIoojiDopzarAMZViwN8SrAADWZprBMngs2A98FQDAJlCDBdAdfBUAwG6wDxaKKHZUbFJflaU4DBTXUJwiiiKKOyiSIgRoDyErAIDNIUUI0Ax8FQAAnGCwAGrBVwEAgMXMGqzc3bB0JV9RXF4xsoXVMnNEcXnFKaIooriD4swaLEMZFiiEeBUAADwy2WAZPBYoAV8FAADpUIMFEANfBQAABbAPFoooeqh8RKCKOaKI4ixRFFHcQZEUIcAPzpAV70kAAKiBFCGAMaQCAQCgKb+u8Nc9Dpb1evyJc9V3OFHpsAtO/JMKfP1NBXKFRZ2odNiKTlQ6bEUnKh02F6r+xPkpQkOWEIZDvAoAALpCihA2Al8FAABjwGDB+uCrAABgMCK2aUh8Zk5DxURQVK14bbXw/jIFWy3UsPBVRXExxSmiKKK4g6KIGixDGRa0g3gVAABMhxQhrANbWAEAgBAwWKAeQlYAACANETVYJq0MS1fyFcXeiulPs9E7RxRRXFIURRR3UJRSg2Uow4I0iFcBAIB8SBGCDvBVAACgCEERLEMQCxzwVQAAoBEpNVgoongnvb6qlWJbUERRi+IUURRR3EGRCBYIgngVAACsgSyDZfBYu8IWVgAAsBIUucNMCFkBAMCSUIOF4gTFtiVWKYpd+0cRRb2KU0RRRHEHRVKEMA7iVQAAsAmkCKE7+CoAANgNcSnCyDNzdMUGUQzlAVeaI4ooalecIooiijsoiksRGrKEyiFeBQAAINFgGTyWQvBVAAAAF9RgQS1sYQUAAGAhrgYLRS2KVonVAMViUEQRRVGiKKK4gyIpQsiDVCAAAMAjQg2WwWMJA18FAACQDjVYEANfBQAAUAA1WCh6FLs+ykbIHFFEEcVZoiiiuIMiKUL4C/EqAACAJsg1WAaPNQp8FQAAQFuowdoatrACAADogegaLOu5hLqSr5IV7yVW768Bgj9Y9aqiiKJGxSmiKKK4g6LoFKEhS9gUUoEAAABjIEW4PvgqAACAwWCwlgVfBQAAMAvRNVjmZxmWruTrLMXcLaw0zhFFFFFULYoiijsoSq/BMpRhpUG8CgAAQA6kgCyoaQAAAzhJREFUCHWDrwIAABDIryv8dY+DZb0ef+JcdSEnvv75vF5/84BcYV0nKh02F2q9E5UOW9GJSofNhao/UUGK0JAlvEHICgAAQD6kCHWArwIAAFAEBks0+CoAAACNSN+m4eTcrGGk4sksxdytFuoVR4IiiijKUZwiiiKKOyjqqMEyyWVYn8/n+1trsId4FQAAwBoslSKc8j9/9eCrAAAAFqM8RTiY+5buXtJjV5YPu3/Z8uR6XTrYJEbmAQEAAGAkOmqw/vBvd4Xv7+/v7+/Tq50veqhEfJWuBDOKKKKoXXGKKIoobqF46MF8vUP/9H6/rRcRrGPcc7N6S8d8vc+fhn0CAACAQDTVYJ1ZQjeVJrywnRIrAACA3dBksCJY+9NL8Fv4KgAAgG3RVIP1+Xy8pe7fN85fs/psOcTq0nVlCWYUUURRueIUURRR3EFRzT5YF/ENsVLCV/frdZW03891X6SM6nxBvAoAAABWM1gpNMwh4qsAAADARV8NVqjUfST4KgAAAIigrAZrumLvrUElzBFFFFHcR3GKKIoo7qCoL0VoWmQJCxTPF4SsAAAA4BGVBsuM8lj4KgAAACiAFKGHa6uF95cZ/JTAha8qiiiiKFBxiiiKKO6gWBXBmvJZcPJ+m/e7fZ/WixMJ25YCAACAIvR9i/C0O+9/O2y1cHg6nGgiAQAAQCn6DNZJ/WYN1FcBAABAJzTVYDUh61E2utK9KKKIIooqRFFEcQdFfTVY98xgVhDrDFkVPByQGiwAAADIQmuKMB1SgQAAADAYfQbrR9js32AUzf1KoNKcJgAAAOjjKOX9fhef21DRfL2tX8+ffopdQRFFFFFcXhRFFHdQ1LqT+8VZhkUeEAAAAOSgL0XoMv7RhAAAAAAR1EewAAAAAKShaR8sFFFEEUUUFxBFEcUdFIlgAQAAADSmPIIFAAAAAF4wWAAAAACNoQYLRRRRRHFfxSmiKKK4gyI1WAAAAACNIUUIAAAA0JhfV/jrHgfLej3+xLnqO5yodNiKTlQ6bC7UeicqHbaiE5UOmwtVfyIpQgAAAIDGkCIEAAAAaAwGCwAAAKAxbNOAIooooriv4hRRFFHcQfH/A5RMB2qIbdLFAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 271,
       "width": 400
      }
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAyAAAAIeCAIAAADMFfxWAAAUtUlEQVR4nO3dUXLbOBYFUNjVC0tWJmllWZrmQz1qh5IcmnzMBYxzKjXlKIH6zszPbQIPfLterw0AgDrv6QAAAN+NggUAUEzBAgAopmABABRTsAAAim0vWJfLpTDHXyZ8ivApwqcInyJ8xLjJW3X4N9c0AADUskUIAFBMwQIAKPbPhjWXt/IYAAD9On3xRJUnWAAAxRQsAIBiChYAQDEFCwCg2JZD7gtfPfYFANC5nSN9nmABABRTsAAAik1asLwsKUX4FOFThE8RPmLc5K2HdxEudiWdwQIAvpmdbWfSJ1gAAMdRsAAAiilYAADFFCwAgGLvj2fmn56iX3O0fuXCbZ8MtFCGnQt7yLB5oQw7F8qwc6EMOxf2kGHzQhl2Llw5RbhyoSlCAIAlU4QAAH1RsAAAiilYAADFFCwAgGIKFgBAsUkLlrdRpgifInyK8CnCR4ybvFWHd00DAMCSaxoAAPqiYAEAFFOwAACKKVgAAMUmLVjGHFKETxE+RfgU4SPGTd5MEQIAHM0UIQBAXxQsAIBiChYAQDEFCwCg2KQFy5hDivApwqcInyJ8xLjJmylCAICjmSIEAOiLggUAUEzBAgAopmABABRTsAAAik1asMyRpgifInyK8CnCR4ybvLmmAQDgaK5pAADoi4IFAFBMwQIAKKZgAQAUm7RgGXNIET5F+BThU4SPGDd5M0UIAHA0U4QAAH1RsAAAiilYAADFFCwAgGKTFixjDinCpwifInyK8BHjJm+mCAEAjmaKEACgLwoWAEAxBQsAoJiCBQBQbNKCZcwhRfgU4VOETxE+YtzkzRQhAMDRTBECAPRFwQIAKKZgAQAUU7AAAIopWAAAxSYtWOZIU4RPET5F+BThI8ZN3lzTAABwNNc0AAD0RcECACimYAEAFFOwAACKTVqwjDmkCJ8ifIrwKcJHjJu8mSIEADiaKUIAgL4oWAAAxRQsAIBiChYAQLH3xzPzT0/Rrzlav3Lhtk8GWijDzoU9ZNi8UIadC2XYuVCGnQt7yLB5oQw7F66cIly50BQhAMCSKUIAgL4oWAAAxRQsAIBiChYAQDEFCwCg2KQFy9soU4RPET5F+BThI8ZN3qrDu6YBAGDJNQ0AAH1RsAAAiilYAADFFCwAgGKTFixjDinCpwifInyK8BHjJm+mCAEAjmaKEACgLwoWAEAxBQsAoJiCBQBQbNKCZcwhRfgU4VOETxE+YtzkzRQhAMDRTBECAPRFwQIAKKZgAQAUU7AAAIopWAAAxSYtWOZIU4RPET5F+BThI8ZN3lzTAABwNNc0AAD0RcECACimYAEAFFOwAACKTVqwjDmkCJ8ifIrwKcJHjJu8mSIEADiaKUIAgL4oWAAAxRQsAIBiChYAQLFJC5YxhxThU4RPET5F+IhxkzdThAAARzNFCADQFwULAKCYggUAUEzBAgAoNmnBMuaQInyK8CnCpwgfMW7yZooQAOBopggBAPqiYAEAFFOwAACKKVgAAMUULACAYpMWLHOkKcKnCJ8ifIrwEeMmb65pAAA4mmsaAAD6omABABRTsAAAiilYAADFJi1YxhxShE8RPkX4FOEjxk3eTBECABzNFCEAQF8ULACAYgoWAEAxBQsAoNj745n5p6fo1xytX7lw2ycDLZRh58IeMmxeKMPOhTLsXCjDzoU9ZNi8UIadC1dOEa5caIoQAGDJFCEAQF8ULACAYgoWAEAxBQsAoJiCBQBQbNKC5W2UKcKnCJ8ifIrwEeMmb9XhXdMAALDkmgYAgL4oWAAAxRQsAIBiChYAQLFJC5YxhxThU4RPET5F+IhxkzdThAAARzNFCADQFwULAKCYggUAUEzBAgAoNmnBMuaQInyK8CnCpwgfMW7yZooQAOBopggBAPqiYAEAFFOwAACKKVgAAMUmLVjGHFKETxE+RfgU4SPGTd5MEQIAHM0UIQBAXxQsAIBiChYAQDEFCwCgmIIFAFBs0oJljjRF+BThU4RPET5i3OTNNQ0AAEdzTQMAQF8ULACAYgoWAEAxBQsAoNikBcuYQ4rwKcKnCJ8ifMS4yZspQgCAo5kiBADoi4IFAFBMwQIAKKZgAQAUm7RgGXNIET5F+BThU4SPGDd5M0UIAHA0U4QAAH1RsAAAiilYAADFFCwAgGIKFgBAsUkLljnSFOFThE8RPkX4iHGTN9c0AAAczTUNAAB9UbAAAIopWAAAxRQsAIBikxYsYw4pwqcInyJ8ivAR4yZvpggBAI5mihAAoC8KFgBAMQULAKCYggUAUOz98cz801P0a47Wr1y47ZOBFsqwc2EPGTYvlGHnQhl2LpRh58IeMmxeKMPOhSunCFcuNEUIALBkihAAoC8KFgBAMQULAKCYggUAUEzBAgAoNmnB8jbKFOFThE8RPkX4iHGTt+rwrmkAAFhyTQMAQF8ULACAYgoWAEAxBQsAoNikBcuYQ4rwKcKnCJ8ifMS4yZspQgCAo5kiBADoi4IFAFBMwQIAKKZgAQAUm7RgGXNIET5F+BThU4SPGDd5M0UIAHA0U4QAAH1RsAAAiilYAADFFCwAgGKTFixjDinCpwifInyK8BHjJm+mCAEAjmaKEACgLwoWAEAxBQsAoJiCBQBQTMECACg2acEyR5oifIrwKcKnCB8xbvLmmgYAgKO5pgEAoC8KFgBAMQULAKCYggUAUGzSgmXMIUX4FOFThE8RPmLc5M0UIQDA0UwRAgD0RcECACimYAEAFFOwAACKTVqwjDmkCJ8ifIrwKcJHjJu8mSIEADiaKUIAgL4oWAAAxRQsAIBiChYAQDEFCwCg2KQFyxxpivApwqcInyJ8xLjJm2saAACO5poGAIC+KFgAAMUULACAYgoWAECxSQuWMYcU4VOETxE+RfiIcZM3U4QAAEczRQgA0BcFCwCgmIIFAFBMwQIAKPb+eGb+6Sn6NUfrVy7c9slAC2XYubCHDJsXyrBzoQw7F8qwc2EPGTYvlGHnwpVThCsXmiIEAFgyRQgA0BcFCwCgmIIFAFBMwQIAKDZpwfKypBThU4RPET5F+Ihxk7fq8KYIAQCWTBECAPRFwQIAKKZgAQAUU7AAAIopWAAAxSYtWOZIU4RPET5F+BThI8ZN3lzTAABwNNc0AAD0RcECACimYAEAFFOwAACKTVqwjDmkCJ8ifIrwKcJHjJu8mSIEADiaKUIAgL4oWAAAxRQsAIBiChYAQLFJC5YxhxThU4RPET5F+IhxkzdThAAARzNFCADQFwULAKCYggUAUEzBAgAopmABABSbtGCZI00RPkX4FOFThI8YN3lzTQMAwNFc0wAA0BcFCwCgmIIFAFBMwQIAKDZpwTLmkCJ8ivApwqcIHzFu8maKEADgaKYIAQD6omABABRTsAAAiilYAADFJi1YxhxShE8RPkX4FOEjxk3eTBECABzNFCEAQF8ULACAYgoWAEAxBQsAoJiCBQBQbNKCZY40RfgU4VOETxE+YtzkzTUNAABHc00DAEBfFCwAgGIKFgBAsYKC9fbz8vZz4ENtAAC1ap5gXX+dbjVrlKZlzCFF+BThU4RPET5i3OSttynCc7tcr6ffvvH/Hev667fPAQBGsXOK8J/CKDf3XqVpAQBzqi9Yd4umpWYBAJPYVbDObdVu5a1aeaAFAExi7xOsc1vblmwdAgCTeH88M//0FP2ao/UrF55/tOuv08fBw80Z+lwow86FPWTYvFCGnQtl2LlQhp0Le8iweaEMOxeunCJcuXDXFOG5Xc7ttP9dhJ5pAQBdyU8Rvn1I8PW21prdQwDge9lesD4/4f7qxtHbzuAnn9y3Du+/3ZwQACBi5xRhQft5WsU80wIAxnXgPVjb3BrV28/LrVFpWgDAcGreRfjoNif4+Kuta0iPf+f+DSUvPfSypBThU4RPET5F+Ihxk7fq8NunCG8jhK2184c/Wvllrx5HLQ5jbfsSAICd8lOEHz2dKHz7PWL78bJdPR6B/4TdQwCgTweewXp7+8KtDR8PXX11+89LDwGArhx1BuuPPtage7sq+c6dx7MAAHY6dopwuTn48EfXa1m7AgDoRFnBuu0GftKonippV186Gg8AcLTYFmH7cWk/nrSr9Rt89/saFtdArGGONEX4FOFThE8RPmLc5K3zaxrWej1IeP/5VVsyMwgAHK2vaxqu1ye7hI8V7u3n/Yf/5v4+b1d6FQAwisNflfP0AdniLobH1z8//pFeBQCMYufLni9/fN/z7Tasx8PsT69UcHcoAPANbC9Y53Y6t0t7eEb16pHV0y95dQZLrwIAxvU3pgj/OBj4OA94dCRjDinCpwifInyK8BHjJm/9TBG2/28RfvVc/W//eC+3AQD609cU4Uq2AgGAb2xvwTq3y+lP59zv9CoAYAY7pwj/Pef+Ob0KAJjKgVuEehUAMKf6KcK/PxK4gTGHFOFThE8RPkX4iHGTt66mCFtr53ZpP367mb3PRgUAsF54ivB2DOvxonYAgGkVnMFSrQAAPvobN7kDAExFwQIAKKZgAQAUm7RgmSNNET5F+BThU4SPGDd56+2ahvb1wUUAgM7tbDuTPsECADiOggUAUEzBAgAopmABABSbtGAZc0gRPkX4FOFThI8YN3kzRQgAcDRThAAAfVGwAACKKVgAAMUULACAYpMWLGMOKcKnCJ8ifIrwEeMmb6YIAQCOZooQAKAvChYAQDEFCwCgmIIFAFBMwQIAKDZpwTJHmiJ8ivApwqcIHzFu8uaaBgCAo7mmAQCgLwoWAEAxBQsAoJiCBQBQbNKCZcwhRfgU4VOETxE+YtzkzRQhAMDRTBECAPRFwQIAKKZgAQAUU7AAAIq9P56Zf3qKfs3R+pULt30y0EIZdi7sIcPmhTLsXCjDzoUy7FzYQ4bNC2XYuXDlFOHKhaYIAQCWTBECAPRFwQIAKKZgAQAUU7AAAIpNWrC8LClF+BThU4RPET5i3OStOrwpQgCAJVOEAAB9UbAAAIopWAAAxRQsAIBiChYAQLFJC5Y50hThU4RPET5F+IhxkzfXNAAAHM01DQAAfVGwAACKKVgAAMUULACAYpMWLGMOKcKnCJ8ifIrwEeMmb6YIAQCOZooQAKAvChYAQDEFCwCgmIIFAFBs0oJlzCFF+BThU4RPET5i3OTNFCEAwNFMEQIA9EXBAgAopmABABRTsAAAiilYAADFJi1Y5khThE8RPkX4FOEjxk3eXNMAAHA01zQAAPRFwQIAKKZgAQAUU7AAAIpNWrCMOaQInyJ8ivApwkeMm7yZIgQAOJopQgCAvihYAADFFCwAgGIKFgBAsUkLljGHFOFThE8RPkX4iHGTN1OEAABHM0UIANAXBQsAoJiCBQBQTMECACg2acEy5pAifIrwKcKnCB8xbvJmihAA4GimCAEA+qJgAQAUU7AAAIopWAAAxRQsAIBikxYsc6QpwqcInyJ8ivAR4yZvrmkAADiaaxoAAPqiYAEAFFOwAACKKVgAAMXeH8/MPz1Fv+Zo/cqF2z4ZaKEMOxf2kGHzQhl2LpRh50IZdi7sIcPmhTLsXLhyinDlQlOEAABLpggBAPqiYAEAFFOwAACKKVgAAMUmLVhelpQifIrwKcKnCB8xbvJWHd4UIQDAkilCAIC+KFgAAMUULACAYgoWAEAxBQsAoNikBcscaYrwKcKnCJ8ifMS4yZtrGgAAjuaaBgCAvihYAADFFCwAgGIKFgBAsUkLljGHFOFThE8RPkX4iHGTN1OEAABHM0UIANAXBQsAoJiCBQBQTMECACg2acEy5pAifIrwKcKnCB8xbvJmihAA4Gj5KcK3n5e3nwM3VgCAWv+UfMv11+nesa6/TiXfCQAwqL0F69wut0Z171WaFgAwuZonWB9pWgDA5OoL1t2iaalZAMAkdh1yP7dVZ9uvv063Q1r9HIc3R5oifIrwKcKnCB8xbvLW1TUN53Y5t9NXBxdtHQIAndt5TcOBW4SvOKQFAHxvgYJ1p2kBAN9SQcF6+/AM7ev7ja1pWgDA97K9YH1+wv3VYfaPV5I+/eR+HP7+280JAQAidk4RFrSfp1XsNnh43OyhMYcU4VOETxE+RfiIcZO3fqYIbyOErbXzhz+6f9n6J1ivLJ5d9fNM63K5nE7LDLf/V9Z/DgD0rNMpwk9q0JqO9bi8h3Nar7rtvXItuterzwGA723XFuF+963A26/75580sKN3Dz9xOp2ePqO6f3g6ne4l7NXnryz+wsfvubn/vO+/BABwuOInWE8nCt9+f8jWfjx/BPX287J+A7H18Uzr71g8BvM8DAA6d+ATrGWv+vwv/7zc6tGGkrR4pvXV5SU+2T38y0kAgLjYFuFiQ7Dk4dPtS9Z0rPLe8+qRkkdNC0M3TuFThE8RPmXc8OMmb9Xhj73J/ZOHWLc/ul7L2tWX6D0pQ/8vL3yK8CnCp4wbftzkrTp8WcG6nbj60rZgK3p2tbindOe3fdXtAPvj6ahXn68x9L8BAAC5dxH+uLRnfWj9Iaq/f7b942Rf+1B170OCixb16vM/fv+aqUMAoFvFF42u9XqQ8P7zq9r0XWcGzQYCQD/6umj0en2yS/hY4d5+3n/47znW5+3qu/YqAOD7OXyK8N+zWQ8veP7428X1CosBw9uvx8tI9xh6A074FOFThE8RPmXc8OMmb129i7C1dm6n84c/WvkE679//Ivi5XkVABAU2yI8t9OtYy3605cKW7dvdAYA2OxvTBH+8S4GvQoA+E62bxG2+yDhefue5fn833/yJUYOAeA4fU0RrnRvVNlqNW5HGfoUIQB8e3sL1rldzuv/8nn5AwDA97PrmobbRaN//mvn5S/2GPrxlfApwqcInyJ8xLjJWz/XNNyc2+VVYer/edXQW4TjhgeA/nV3Bqv/XgUAcKiCgrXY+BuoVw39JBMA6NbeLcL24Ur3qkwAAF0JbBGqVgAAH70/bpM93TizmwYAsLI4FWwRAgB8b1/dItx1DxYAAI+2PMECAOATnmABABRTsAAAim0vWEPPFQqfInyK8CnCpwgfMW7y1sO7CAEA+IQtQgCAYgoWAEAxBQsAoJiCBQBQTMECACjmmobxCJ8ifIrwKcKnjBt+3OStOvz/AEhgdLVwn/MqAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 271,
       "width": 400
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "from SiEPIC.utils import create_cell2\n",
    "cell = create_cell2(ly, 'GC_TE_1310_8degOxide_BB', 'EBeam')\n",
    "t = pya.Trans(pya.Trans.R0, 40e3,15e3)\n",
    "inst_gc1 = topcell.insert(pya.CellInstArray(cell.cell_index(), t))\n",
    "\n",
    "t = pya.Trans(pya.Trans.R0, 40e3,15e3+127e3)\n",
    "inst_gc2 = topcell.insert(pya.CellInstArray(cell.cell_index(), t))\n",
    "\n",
    "cell.plot(width = 400)\n",
    "topcell.plot(width = 400)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "502aa6ce",
   "metadata": {},
   "source": [
    "## Components and their pins\n",
    "- Query a cell to get the SiEPIC Component description\n",
    "- Get the pins associated with this Component\n",
    "- Alternative: Query a cell to get a list of Pins"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "1725d0b8",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Display component details:\n",
      "<class 'SiEPIC.core.Component'>\n",
      "- basic_name: GC_TE_1310_8degOxide_BB, component: GC_TE_1310_8degOxide_BB-0 / GC_TE_1310_8degOxide_BB; transformation: r0 *1 0,0; center position: -19.985,0.084; number of pins: 2; optical pins: [['opt1', '0,0', None]]; electrical pins: []; optical IO pins: [['GC_TE_1310_8degOxide_BB', '-20400,0', None]]; has compact model: False; params: .\n",
      "\n",
      "Display pin details:\n",
      "- pin_name opt1: component_idx 0, pin_type 1, rotation: 0.0, net: None, (0,0), path: (-10,0;10,0) w=350 bx=0 ex=0 r=false\n",
      "- pin_name GC_TE_1310_8degOxide_BB: component_idx 0, pin_type 0, rotation: 0, net: None, (-20400,0), path: None\n",
      "\n",
      "Display pin names:\n",
      " -  opt1\n"
     ]
    }
   ],
   "source": [
    "\n",
    "print('Display component details:')\n",
    "component = cell.find_components()[0]  # find component within a cell, return the first one\n",
    "print(type(component))\n",
    "component.display()\n",
    "\n",
    "print('\\nDisplay pin details:')\n",
    "pins = cell.find_pins_component(component)\n",
    "for p in pins:\n",
    "    p.display()\n",
    "    \n",
    "# or only pins\n",
    "print('\\nDisplay pin names:')\n",
    "pins, _ = cell.find_pins()\n",
    "for p in pins:\n",
    "    if p.type == 1:  # optical pin type, defined in SiEPIC._globals\n",
    "        print(' - ', p.pin_name)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "166e3def",
   "metadata": {},
   "source": [
    "## Waveguides\n",
    "- Find out what types of waveguides are defined in the PDK\n",
    "- Connect two components using a waveguide\n",
    "- Specify the waveguide type\n",
    "- Specify the path the waveguide should take\n",
    "- Return the wavelength details (length)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "d01e8e3e",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Waveguide types:\n",
      "  Strip TE 1550 nm, w=500 nm\n",
      "  Strip TE 1310 nm, w=410 nm\n",
      "  Strip TE 1310 nm, w=350 nm\n",
      "  Strip TM 1550 nm, w=500 nm\n",
      "  Strip TE-TM 1550, w=450 nm\n",
      "  Multimode Strip TE 1550 nm, w=2000 nm\n",
      "  Multimode Strip TE 1550 nm, w=3000 nm\n",
      "  Slot TE 1550 nm, w=500 nm, gap=100nm\n",
      "  eskid TE 1550\n",
      "  Rib (90 nm slab) TE 1550 nm, w=500 nm\n",
      "  Rib (90 nm slab) TE 1310 nm, w=350 nm\n",
      "  Si routing TE 1550 nm (compound waveguide)\n",
      "  Si routing TE 1310 nm (compound waveguide)\n",
      "  SiN Strip TE 895 nm, w=450 nm\n",
      "  SiN Strip TE 1550 nm, w=750 nm\n",
      "  SiN Strip TE 1550 nm, w=800 nm\n",
      "  SiN Strip TE 1550 nm, w=1000 nm\n",
      "  SiN Strip TM 1550 nm, w=1000 nm\n",
      "  SiN Strip TE 1310 nm, w=750 nm\n",
      "  SiN Strip TM 1310 nm, w=750 nm\n",
      "  Multimode SiN Strip TE 1550 nm, w=3000 nm\n",
      "  SiN routing TE 1550 nm (compound waveguide)\n",
      "190.793\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABkAAAARLCAIAAABTNvgaAAAgAElEQVR4nOzd0ZFc55Vm0R8TY8yYIVkx4w0Bb7qtUJvR3mAekAKg4mXifiRRudVcKxAhHjCB3EW9nbh58nz+p48fP37+jvF3j50S45OxU2K8P3ZKjE/GTonx/tgpMT4ZOyXG+2OnxPhk7JQY74+dEuOTsVNivD92Sn49fvj8+fMBAAAAgKr/9eoAAAAAAHjGAgsAAACANAssAAAAANJesMD69OnT+7/pDzWrTjVM1aQZpmrSDGtWnWqYqkkzTNWkGaZq0gxTNWmGqZo0w1RNmmHNqlMN+/TpkyPuAAAAAKT5CCEAAAAAaRZYAAAAAKRZYAEAAACQ5oj7Q7PqVMNUTZphqibNsGbVqYapmjTDVE2aYaomzTBVk2aYqkkzTNWkGdasOtUwR9wBAAAAqPMRQgAAAADSLLAAAAAASLPAAgAAACDNEfeHZtWphqmaNMNUTZphzapTDVM1aYapmjTDVE2aYaomzTBVk2aYqkkzrFl1qmGOuAMAAABQ5yOEAAAAAKRZYAEAAACQZoEFAAAAQJoj7g/NqlMNUzVphqmaNMOaVacapmrSDFM1aYapmjTDVE2aYaomzTBVk2ZYs+pUwxxxBwAAAKDORwgBAAAASPvf7/lmnz6857sBAAAA8LP88o4f6vMEFgAAAABpFlgAAAAAdL33EXcfIQQAAAD4n8FHCAEAAADg4V2PuL/xnos6AAAAAP6IF360zhNYAAAAAKRZYD18+vTp1QnXmmGqJs0wVZNmWLPqVMNUTZphqibNMFWTZpiqSTNM1aQZpmrSDGtWnWrYi4+4+wghAAAAwL+LFy52PIEFAAAAQJoFFgAAAABpFlgAAAAApFlgPTSvlJ1qmKpJM0zVpBnWrDrVMFWTZpiqSTNM1aQZpmrSDFM1aYapmjTDmlWnGuaIOwAAAAC3OOIOAAAAANcssAAAAABIs8ACAAAAIM0C66F5pexUw1RNmmGqJs2wZtWphqmaNMNUTZphqibNMFWTZpiqSTNM1aQZ1qw61TBH3AEAAAC4xRF3AAAAALhmgQUAAABAmgUWAAAAAGkWWA/NK2WnGqZq0gxTNWmGNatONUzVpBmmatIMUzVphqmaNMNUTZphqibNsGbVqYY54g4AAADALY64AwAAAMA1CywAAAAA0iywAAAAAEj7tsB6c6br541vvNv7vucPGCkxPhk7Jcb7Y6fE+GTslBjvj50S45OxU2K8P3ZKjE/GTonx/tgpMT4ZOyXG++OdF7yq0xF3AAAAAH7MEXcAAAAAuGaBBQAAAECaBRYAAAAAaRZYD8/Pkr1QM0zVpBmmatIMa1adapiqSTNM1aQZpmrSDFM1aYapmjTDVE2aYc2qUw1zxB0AAACAWxxxBwAAAIBrFlgAAAAApFlgAQAAAJBmgfXQvFJ2qmGqJs0wVZNmWLPqVMNUTZphqibNMFWTZpiqSTNM1aQZpmrSDGtWnWqYI+4AAAAA3OKIOwAAAABcs8ACAAAAIM0CCwAAAIA0C6yH5pWyUw1TNWmGqZo0w5pVpxqmatIMUzVphqmaNMNUTZphqibNMFWTZliz6lTDHHEHAAAA4BZH3AEAAADgmgUWAAAAAGkWWAAAAACkWWA9NK+UnWqYqkkzTNWkGdasOtUwVZNmmKpJM0zVpBmmatIMUzVphqmaNMOaVaca5og7AAAAALc44g4AAAAA1yywAAAAAEizwAIAAAAgzQLroXml7FTDVE2aYaomzbBm1amGqZo0w1RNmmGqJs0wVZNmmKpJM0zVpBnWrDrVMEfcAQAAALjFEXcAAAAAuGaBBQAAAECaBRYAAAAAaRZYD80rZacapmrSDFM1aYY1q041TNWkGaZq0gxTNWmGqZo0w1RNmmGqJs2wZtWphjniDgAAAMAtjrgDAAAAwDULLAAAAADSLLAAAAAASLPAemheKTvVMFWTZpiqSTOsWXWqYaomzTBVk2aYqkkzTNWkGaZq0gxTNWmGNatONcwRdwAAAABuccQdAAAAAK5ZYAEAAACQZoEFAAAAQJoF1kPzStmphqmaNMNUTZphzapTDVM1aYapmjTDVE2aYaomzTBVk2aYqkkzrFl1qmGOuAMAAABwiyPuAAAAAHDNAgsAAACANAssAAAAANIssB6aV8pONUzVpBmmatIMa1adapiqSTNM1aQZpmrSDFM1aYapmjTDVE2aYc2qUw1zxB0AAACAWxxxBwAAAIBrFlgAAAAApFlgAQAAAJD2bYH15kzXzxvfeLf3fc8fMFJifDJ2Soz3x06J8cnYKTHeHzslxidjp8R4f+yUGJ+MnRLj/bFTYnwydkqM98c7L3hVpyPuAAAAAPyYI+4AAAAAcM0CCwAAAIA0CywAAAAA0iywHp6fJXuhZpiqSTNM1aQZ1qw61TBVk2aYqkkzTNWkGaZq0gxTNWmGqZo0w5pVpxrmiDsAAAAAtzjiDgAAAADXLLAAAAAASLPAAgAAACDNAuuheaXsVMNUTZphqibNsGbVqYapmjTDVE2aYaomzTBVk2aYqkkzTNWkGdasOtUwR9wBAAAAuMURdwAAAAC4ZoEFAAAAQJoFFgAAAABpFlgPzStlpxqmatIMUzVphjWrTjVM1aQZpmrSDFM1aYapmjTDVE2aYaomzbBm1amGOeIOAAAAwC2OuAMAAADANQssAAAAANIssAAAAABIs8B6aF4pO9UwVZNmmKpJM6xZdaphqibNMFWTZpiqSTNM1aQZpmrSDFM1aYY1q041zBF3AAAAAG5xxB0AAAAArllgAQAAAJBmgQUAAABAmgUWAAAAAGkWWA/NM/unGqZq0gxTNWmGNatONUzVpBmmatIMUzVphqmaNMNUTZphqibNsGbVqYb5FkIAAAAAbvEthAAAAABwzQILAAAAgDQLLAAAAADSLLAemlfKTjVM1aQZpmrSDGtWnWqYqkkzTNWkGaZq0gxTNWmGqZo0w1RNmmHNqlMNc8QdAAAAgFsccQcAAACAaxZYAAAAAKRZYAEAAACQZoH10LxSdqphqibNMFWTZliz6lTDVE2aYaomzTBVk2aYqkkzTNWkGaZq0gxrVp1qmCPuAAAAANziiDsAAAAAXLPAAgAAACDNAgsAAACANAush+aVslMNUzVphqmaNMOaVacapmrSDFM1aYapmjTDVE2aYaomzTBVk2ZYs+pUwxxxBwAAAOAWR9wBAAAA4JoFFgAAAABpFlgAAAAApFlgPTSvlJ1qmKpJM0zVpBnWrDrVMFWTZpiqSTNM1aQZpmrSDFM1aYapmjTDmlWnGuaIOwAAAAC3OOIOAAAAANcssAAAAABIs8ACAAAAIO3bAuvNma6fN77xbu/7nj9gpMT4ZOyUGO+PnRLjk7FTYrw/dkqMT8ZOifH+2CkxPhk7Jcb7Y6fE+GTslBjvj3de8KpOR9wBAAAA+DFH3AEAAADgmgUWAAAAAGkWWAAAAACkWWA9PD9L9kLNMFWTZpiqSTOsWXWqYaomzTBVk2aYqkkzTNWkGaZq0gxTNWmGNatONcwRdwAAAABuccQdAAAAAK5ZYAEAAACQZoEFAAAAQJoF1kPzStmphqmaNMNUTZphzapTDVM1aYapmjTDVE2aYaomzTBVk2aYqkkzrFl1qmGOuAMAAABwiyPuAAAAAHDNAgsAAACANAssAAAAANIssB6aV8pONUzVpBmmatIMa1adapiqSTNM1aQZpmrSDFM1aYapmjTDVE2aYc2qUw1zxB0AAACAWxxxBwAAAIBrFlgAAAAApFlgAQAAAJBmgfXQvFJ2qmGqJs0wVZNmWLPqVMNUTZphqibNMFWTZpiqSTNM1aQZpmrSDGtWnWqYI+4AAAAA3OKIOwAAAABcs8ACAAAAIM0CCwAAAIA0C6yH5pWyUw1TNWmGqZo0w5pVpxqmatIMUzVphqmaNMNUTZphqibNMFWTZliz6lTDHHEHAAAA4BZH3AEAAADgmgUWAAAAAGkWWAAAAACkWWA9NK+UnWqYqkkzTNWkGdasOtUwVZNmmKpJM0zVpBmmatIMUzVphqmaNMOaVaca5og7AAAAALc44g4AAAAA1yywAAAAAEizwAIAAAAgzQLroXml7FTDVE2aYaomzbBm1amGqZo0w1RNmmGqJs0wVZNmmKpJM0zVpBnWrDrVMEfcAQAAALjFEXcAAAAAuGaBBQAAAECaBRYAAAAAaRZYD80rZacapmrSDFM1aYY1q041TNWkGaZq0gxTNWmGqZo0w1RNmmGqJs2wZtWphjniDgAAAMAtjrgDAAAAwDULLAAAAADSLLAAAAAASLPAemheKTvVMFWTZpiqSTOsWXWqYaomzTBVk2aYqkkzTNWkGaZq0gxTNWmGNatONcwRdwAAAABuccQdAAAAAK5ZYAEAAACQZoEFAAAAQNq3BdabM10/b3zj3d73PX/ASInxydgpMd4fOyXGJ2OnxHh/7JQYn4ydEuP9sVNifDJ2Soz3x06J8cnYKTHeH++84FWdjrgDAAAA8GOOuAMAAADANQssAAAAANIssAAAAABIs8B6eH6W7IWaYaomzTBVk2ZYs+pUw1RNmmGqJs0wVZNmmKpJM0zVpBmmatIMa1adapgj7gAAAADc4og7AAAAAFyzwAIAAAAgzQILAAAAgDQLrIfmlbJTDVM1aYapmjTDmlWnGqZq0gxTNWmGqZo0w1RNmmGqJs0wVZNmWLPqVMMccQcAAADgFkfcAQAAAOCaBRYAAAAAaRZYAAAAAKRZYD00r5SdapiqSTNM1aQZ1qw61TBVk2aYqkkzTNWkGaZq0gxTNWmGqZo0w5pVpxrmiDsAAAAAtzjiDgAAAADXLLAAAAAASLPAAgAAACDNAuuheaXsVMNUTZphqibNsGbVqYapmjTDVE2aYaomzTBVk2aYqkkzTNWkGdasOtUwR9wBAAAAuMURdwAAAAC4ZoEFAAAAQJoFFgAAAABpFlgPzStlpxqmatIMUzVphjWrTjVM1aQZpmrSDFM1aYapmjTDVE2aYaomzbBm1amGOeIOAAAAwC2OuAMAAADANQssAAAAANIssAAAAABIs8B6aF4pO9UwVZNmmKpJM6xZdaphqibNMFWTZpiqSTNM1aQZpmrSDFM1aYY1q041zBF3AAAAAG5xxB0AAAAArllgAQAAAJBmgQUAAABAmgXWQ/NK2amGqZo0w1RNmmHNqlMNUzVphqmaNMNUTZphqibNMFWTZpiqSTOsWXWqYY64AwAAAHCLI+4AAAAAcM0CCwAAAIA0CywAAAAA0iywHppXyk41TNWkGaZq0gxrVp1qmKpJM0zVpBmmatIMUzVphqmaNMNUTZphzapTDXPEHQAAAIBbHHEHAAAAgGsWWAAAAACkWWABAAAAkGaB9dC8UnaqYaomzTBVk2ZYs+pUw1RNmmGqJs0wVZNmmKpJM0zVpBmmatIMa1adapgj7gAAAADc4og7AAAAAFyzwAIAAAAgzQILAAAAgLRvC6w3Z7p+3vjGu73ve/6AkRLjk7FTYrw/dkqMT8ZOifH+2CkxPhk7Jcb7Y6fE+GTslBjvj50S45OxU2K8P955was6HXEHAAAA4McccQcAAACAaxZYAAAAAKRZYAEAAACQZoH18Pws2Qs1w1RNmmGqJs2wZtWphqmaNMNUTZphqibNMFWTZpiqSTNM1aQZ1qw61TBH3AEAAAC4xRF3AAAAALhmgQUAAABAmgUWAAAAAGkWWA/NK2WnGqZq0gxTNWmGNatONUzVpBmmatIMUzVphqmaNMNUTZphqibNsGbVqYY54g4AAADALY64AwAAAMA1CywAAAAA0iywAAAAAEizwHpoXik71TBVk2aYqkkzrFl1qmGqJs0wVZNmmKpJM0zVpBmmatIMUzVphjWrTjXMEXcAAAAAbnHEHQAAAACuWWABAAAAkGaBBQAAAECaBRYAAAAAaRZYD80z+6capmrSDFM1aYY1q041TNWkGaZq0gxTNWmGqZo0w1RNmmGqJs2wZtWphvkWQgAAAABu8S2EAAAAAHDNAgsAAACANAssAAAAANIssB6aV8pONUzVpBmmatIMa1adapiqSTNM1aQZpmrSDFM1aYapmjTDVE2aYc2qUw1zxB0AAACAWxxxBwAAAIBrFlgAAAAApFlgAQAAAJBmgfXQvFJ2qmGqJs0wVZNmWLPqVMNUTZphqibNMFWTZpiqSTNM1aQZpmrSDGtWnWqYI+4AAAAA3OKIOwAAAABcs8ACAAAAIM0CCwAAAIA0C6yH5pWyUw1TNWmGqZo0w5pVpxqmatIMUzVphqmaNMNUTZphqibNMFWTZliz6lTDHHEHAAAA4BZH3AEAAADgmgUWAAAAAGkWWAAAAACkWWA9NK+UnWqYqkkzTNWkGdasOtUwVZNmmKpJM0zVpBmmatIMUzVphqmaNMOaVaca5og7AAAAALc44g4AAAAA1yywAAAAAEizwAIAAAAgzQLroXml7FTDVE2aYaomzbBm1amGqZo0w1RNmmGqJs0wVZNmmKpJM0zVpBnWrDrVMEfcAQAAALjFEXcAAAAAuGaBBQAAAECaBRYAAAAAad8WWG/OdP288Y13e9/3/AEjJcYnY6fEeH/slBifjJ0S4/2xU2J8MnZKjPfHTonxydgpMd4fOyXGJ2OnxHh/vPOCV3U64g4AAADAjzniDgAAAADXLLAAAAAASLPAAgAAACDNAuvh+VmyF2qGqZo0w1RNmmHNqlMNUzVphqmaNMNUTZphqibNMFWTZpiqSTOsWXWqYY64AwAAAHCLI+4AAAAAcM0CCwAAAIA0CywAAAAA0iywHppXyk41TNWkGaZq0gxrVp1qmKpJM0zVpBmmatIMUzVphqmaNMNUTZphzapTDXPEHQAAAIBbHHEHAAAAgGsWWAAAAACkWWABAAAAkGaB9dC8UnaqYaomzTBVk2ZYs+pUw1RNmmGqJs0wVZNmmKpJM0zVpBmmatIMa1adapgj7gAAAADc4og7AAAAAFyzwAIAAAAgzQILAAAAgDQLrIfmlbJTDVM1aYapmjTDmlWnGqZq0gxTNWmGqZo0w1RNmmGqJs0wVZNmWLPqVMMccQcAAADgFkfcAQAAAOCaBRYAAAAAaRZYAAAAAKRZYD00r5SdapiqSTNM1aQZ1qw61TBVk2aYqkkzTNWkGaZq0gxTNWmGqZo0w5pVpxrmiDsAAAAAtzjiDgAAAADXLLAAAAAASLPAAgAAACDNAuuheaXsVMNUTZphqibNsGbVqYapmjTDVE2aYaomzTBVk2aYqkkzTNWkGdasOtUwR9wBAAAAuMURdwAAAAC4ZoEFAAAAQJoFFgAAAABpFlgPzStlpxqmatIMUzVphjWrTjVM1aQZpmrSDFM1aYapmjTDVE2aYaomzbBm1amGOeIOAAAAwC2OuAMAAADANQssAAAAANIssAAAAABIs8B6aF4pO9UwVZNmmKpJM6xZdaphqibNMFWTZpiqSTNM1aQZpmrSDFM1aYY1q041zBF3AAAAAG5xxB0AAAAArllgAQAAAJBmgQUAAABAmgXWQ/NK2amGqZo0w1RNmmHNqlMNUzVphqmaNMNUTZphqibNMFWTZpiqSTOsWXWqYY64AwAAAHCLI+4AAAAAcM0CCwAAAIA0CywAAAAA0r4tsN6c6fp54xvv9r7v+QNGSoxPxk6J8f7YKTE+GTslxvtjp8T4ZOyUGO+PnRLjk7FTYrw/dkqMT8ZOifH+eOcFr+p0xB0AAACAH3PEHQAAAACuWWABAAAAkGaBBQAAAECaBdbD87NkL9QMUzVphqmaNMOaVacapmrSDFM1aYapmjTDVE2aYaomzTBVk2ZYs+pUwxxxBwAAAOAWR9wBAAAA4JoFFgAAAABpFlgAAAAApFlgPTSvlJ1qmKpJM0zVpBnWrDrVMFWTZpiqSTNM1aQZpmrSDFM1aYapmjTDmlWnGuaIOwAAAAC3OOIOAAAAANcssAAAAABIs8ACAAAAIM0C66F5pexUw1RNmmGqJs2wZtWphqmaNMNUTZphqibNMFWTZpiqSTNM1aQZ1qw61TBH3AEAAAC4xRF3AAAAALhmgQUAAABAmgUWAAAAAGkWWA/NK2WnGqZq0gxTNWmGNatONUzVpBmmatIMUzVphqmaNMNUTZphqibNsGbVqYY54g4AAADALY64AwAAAMA1CywAAAAA0iywAAAAAEizwHpoXik71TBVk2aYqkkzrFl1qmGqJs0wVZNmmKpJM0zVpBmmatIMUzVphjWrTjXMEXcAAAAAbnHEHQAAAACuWWABAAAAkGaBBQAAAECaBdZD80rZqYapmjTDVE2aYc2qUw1TNWmGqZo0w1RNmmGqJs0wVZNmmKpJM6xZdaphjrgDAAAAcIsj7gAAAABwzQILAAAAgDQLLAAAAADSLLAemlfKTjVM1aQZpmrSDGtWnWqYqkkzTNWkGaZq0gxTNWmGqZo0w1RNmmHNqlMNc8QdAAAAgFsccQcAAACAaxZYAAAAAKRZYAEAAACQZoH10LxSdqphqibNMFWTZliz6lTDVE2aYaomzTBVk2aYqkkzTNWkGaZq0gxrVp1qmCPuAAAAANziiDsAAAAAXLPAAgAAACDNAgsAAACANAush+aVslMNUzVphqmaNMOaVacapmrSDFM1aYapmjTDVE2aYaomzTBVk2ZYs+pUwxxxBwAAAOAWR9wBAAAA4JoFFgAAAABpFlgAAAAApH1bYL050/Xzxjfe7X3f8weMlBifjJ0S4/2xU2J8MnZKjPfHTonxydgpMd4fOyXGJ2OnxHh/7JQYn4ydEuP98c4LXtXpiDsAAAAAP+aIOwAAAABcs8ACAAAAIM0CCwAAAIA0C6yH52fJXqgZpmrSDFM1aYY1q041TNWkGaZq0gxTNWmGqZo0w1RNmmGqJs2wZtWphjniDgAAAMAtjrgDAAAAwDULLAAAAADSLLAAAAAASLPAAgAAACDNAuuheWb/VMNUTZphqibNsGbVqYapmjTDVE2aYaomzTBVk2aYqkkzTNWkGdasOtUw30IIAAAAwC2+hRAAAAAArllgAQAAAJBmgQUAAABAmgXWQ/NK2amGqZo0w1RNmmHNqlMNUzVphqmaNMNUTZphqibNMFWTZpiqSTOsWXWqYY64AwAAAHCLI+4AAAAAcM0CCwAAAIA0CywAAAAA0iywHppXyk41TNWkGaZq0gxrVp1qmKpJM0zVpBmmatIMUzVphqmaNMNUTZphzapTDXPEHQAAAIBbHHEHAAAAgGsWWAAAAACkWWABAAAAkGaB9dC8UnaqYaomzTBVk2ZYs+pUw1RNmmGqJs0wVZNmmKpJM0zVpBmmatIMa1adapgj7gAAAADc4og7AAAAAFyzwAIAAAAgzQILAAAAgDQLrIfmlbJTDVM1aYapmjTDmlWnGqZq0gxTNWmGqZo0w1RNmmGqJs0wVZNmWLPqVMMccQcAAADgFkfcAQAAAOCaBRYAAAAAaRZYAAAAAKRZYD00r5SdapiqSTNM1aQZ1qw61TBVk2aYqkkzTNWkGaZq0gxTNWmGqZo0w5pVpxrmiDsAAAAAtzjiDgAAAADXLLAAAAAASLPAAgAAACDNAuuheaXsVMNUTZphqibNsGbVqYapmjTDVE2aYaomzTBVk2aYqkkzTNWkGdasOtUwR9wBAAAAuMURdwAAAAC4ZoEFAAAAQJoFFgAAAABpFlgPzStlpxqmatIMUzVphjWrTjVM1aQZpmrSDFM1aYapmjTDVE2aYaomzbBm1amGOeIOAAAAwC2OuAMAAADANQssAAAAANIssAAAAABI+7bAenOm6+eNb7zb+77nDxgpMT4ZOyXG+2OnxPhk7JQY74+dEuOTsVNivD92SoxPxk6J8f7YKTE+GTslxvvjnRe8qtMRdwAAAAB+zBF3AAAAALhmgQUAAABAmgUWAAAAAGkWWA/Pz5K9UDNM1aQZpmrSDGtWnWqYqkkzTNWkGaZq0gxTNWmGqZo0w1RNmmHNqlMNc8QdAAAAgFsccQcAAACAaxZYAAAAAKRZYAEAAACQZoH10LxSdqphqibNMFWTZliz6lTDVE2aYaomzTBVk2aYqkkzTNWkGaZq0gxrVp1qmCPuAAAAANziiDsAAAAAXLPAAgAAACDNAgsAAACANAush+aVslMNUzVphqmaNMOaVacapmrSDFM1aYapmjTDVE2aYaomzTBVk2ZYs+pUwxxxBwAAAOAWR9wBAAAA4JoFFgAAAABpFlgAAAAApFlgPTSvlJ1qmKpJM0zVpBnWrDrVMFWTZpiqSTNM1aQZpmrSDFM1aYapmjTDmlWnGuaIOwAAAAC3OOIOAAAAANcssAAAAABIs8ACAAAAIM0C66F5pexUw1RNmmGqJs2wZtWphqmaNMNUTZphqibNMFWTZpiqSTNM1aQZ1qw61TBH3AEAAAC4xRF3AAAAALhmgQUAAABAmgUWAAAAAGkWWA/NK2WnGqZq0gxTNWmGNatONUzVpBmmatIMUzVphqmaNMNUTZphqibNsGbVqYY54g4AAADALY64AwAAAMA1CywAAAAA0iywAAAAAEizwHpoXik71TBVk2aYqkkzrFl1qmGqJs0wVZNmmKpJM0zVpBmmatIMUzVphjWrTjXMEXcAAAAAbnHEHQAAAACuWWABAAAAkGaBBQAAAECaBdZD80rZqYapmjTDVE2aYc2qUw1TNWmGqZo0w1RNmmGqJs0wVZNmmKpJM6xZdaphjrgDAAAAcIsj7gAAAABwzQILAAAAgDQLLAAAAADSLLAemlfKTjVM1aQZpmrSDGtWnWqYqkkzTNWkGaZq0gxTNWmGqZo0w1RNmmHNqlMNc8QdAAAAgFsccQcAAACAaxZYAAAAAKRZYAEAAACQ9m2B9eZM188b33i3933PHzBSYnwydkqM98dOifHJ2Ckx3h87JcYnY6fEeH/slBifjJ0S4/2xU2J8MnZKjPfHOy94Vess6HQAACAASURBVKcj7gAAAAD8mCPuAAAAAHDNAgsAAACANAssAAAAANIssB6enyV7oWaYqkkzTNWkGdasOtUwVZNmmKpJM0zVpBmmatIMUzVphqmaNMOaVaca5og7AAAAALc44g4AAAAA1yywAAAAAEizwAIAAAAgzQLroXml7FTDVE2aYaomzbBm1amGqZo0w1RNmmGqJs0wVZNmmKpJM0zVpBnWrDrVMEfcAQAAALjFEXcAAAAAuGaBBQAAAECaBRYAAAAAaRZYD80rZacapmrSDFM1aYY1q041TNWkGaZq0gxTNWmGqZo0w1RNmmGqJs2wZtWphjniDgAAAMAtjrgDAAAAwDULLAAAAADSLLAAAAAASLPAemheKTvVMFWTZpiqSTOsWXWqYaomzTBVk2aYqkkzTNWkGaZq0gxTNWmGNatONcwRdwAAAABuccQdAAAAAK5ZYAEAAACQZoEFAAAAQJoF1kPzStmphqmaNMNUTZphzapTDVM1aYapmjTDVE2aYaomzTBVk2aYqkkzrFl1qmGOuAMAAABwiyPuAAAAAHDNAgsAAACANAssAAAAANIssB6aV8pONUzVpBmmatIMa1adapiqSTNM1aQZpmrSDFM1aYapmjTDVE2aYc2qUw1zxB0AAACAWxxxBwAAAIBrFlgAAAAApFlgAQAAAJBmgfXQvFJ2qmGqJs0wVZNmWLPqVMNUTZphqibNMFWTZpiqSTNM1aQZpmrSDGtWnWqYI+4AAAAA3OKIOwAAAABcs8ACAAAAIM0CCwAAAIA0C6yH5pWyUw1TNWmGqZo0w5pVpxqmatIMUzVphqmaNMNUTZphqibNMFWTZliz6lTDHHEHAAAA4BZH3AEAAADgmgUWAAAAAGkWWAAAAACkWWA9NK+UnWqYqkkzTNWkGdasOtUwVZNmmKpJM0zVpBmmatIMUzVphqmaNMOaVaca5og7AAAAALc44g4AAAAA1yywAAAAAEizwAIAAAAg7dsC682Zrp83vvFu7/ueP2CkxPhk7JQY74+dEuOTsVNivD92SoxPxk6J8f7YKTE+GTslxvtjp8T4ZOyUGO+Pd17wqk5H3AEAAAD4MUfcAQAAAOCaBRYAAAAAaRZYAAAAAKRZYAEAAACQZoH18Pyu/gs1w1RNmmGqJs2wZtWphqmaNMNUTZphqibNMFWTZpiqSTNM1aQZ1qw61TDfQggAAADALb6FEAAAAACuWWABAAAAkGaBBQAAAECaBdZD80rZqYapmjTDVE2aYc2qUw1TNWmGqZo0w1RNmmGqJs0wVZNmmKpJM6xZdaphjrgDAAAAcIsj7gAAAABwzQILAAAAgDQLLAAAAADSLLAemlfKTjVM1aQZpmrSDGtWnWqYqkkzTNWkGaZq0gxTNWmGqZo0w1RNmmHNqlMNc8QdAAAAgFsccQcAAACAaxZYAAAAAKRZYAEAAACQZoH10LxSdqphqibNMFWTZliz6lTDVE2aYaomzTBVk2aYqkkzTNWkGaZq0gxrVp1qmCPuAAAAANziiDsAAAAAXLPAAgAAACDNAgsAAACANAush+aVslMNUzVphqmaNMOaVacapmrSDFM1aYapmjTDVE2aYaomzTBVk2ZYs+pUwxxxBwAAAOAWR9wBAAAA4JoFFgAAAABpFlgAAAAApFlgPTSvlJ1qmKpJM0zVpBnWrDrVMFWTZpiqSTNM1aQZpmrSDFM1aYapmjTDmlWnGuaIOwAAAAC3OOIOAAAAANcssAAAAABIs8ACAAAAIM0C66F5pexUw1RNmmGqJs2wZtWphqmaNMNUTZphqibNMFWTZpiqSTNM1aQZ1qw61TBH3AEAAAC4xRF3AAAAALhmgQUAAABAmgUWAAAAAGkWWA/NK2WnGqZq0gxTNWmGNatONUzVpBmmatIMUzVphqmaNMNUTZphqibNsGbVqYY54g4AAADALY64AwAAAMA1CywAAAAA0iywAAAAAEizwHpoXik71TBVk2aYqkkzrFl1qmGqJs0wVZNmmKpJM0zVpBmmatIMUzVphjWrTjXMEXcAAAAAbnHEHQAAAACuWWABAAAAkGaBBQAAAEDatwXWmzNdP298493e9z1/wEiJ8cnYKTHeHzslxidjp8R4f+yUGJ+MnRLj/bFTYnwydkqM98dOifHJ2Ckx3h/vvOBVnY64AwAAAPBjjrgDAAAAwDULLAAAAADSLLAAAAAASLPAenh+luyFmmGqJs0wVZNmWLPqVMNUTZphqibNMFWTZpiqSTNM1aQZpmrSDGtWnWqYI+4AAAAA3OKIOwAAAABcs8ACAAAAIM0CCwAAAIA0C6yH5pWyUw1TNWmGqZo0w5pVpxqmatIMUzVphqmaNMNUTZphqibNMFWTZliz6lTDHHEHAAAA4BZH3AEAAADgmgUWAAAAAGkWWAAAAACkWWA9NK+UnWqYqkkzTNWkGdasOtUwVZNmmKpJM0zVpBmmatIMUzVphqmaNMOaVaca5og7AAAAALc44g4AAAAA1yywAAAAAEizwAIAAAAgzQLroXml7FTDVE2aYaomzbBm1amGqZo0w1RNmmGqJs0wVZNmmKpJM0zVpBnWrDrVMEfcAQAAALjFEXcAAAAAuGaBBQAAAECaBRYAAAAAaRZYD80rZacapmrSDFM1aYY1q041TNWkGaZq0gxTNWmGqZo0w1RNmmGqJs2wZtWphjniDgAAAMAtjrgDAAAAwDULLAAAAADSLLAAAAAASLPAemheKTvVMFWTZpiqSTOsWXWqYaomzTBVk2aYqkkzTNWkGaZq0gxTNWmGNatONcwRdwAAAABuccQdAAAAAK5ZYAEAAACQZoEFAAAAQJoF1kPzStmphqmaNMNUTZphzapTDVM1aYapmjTDVE2aYaomzTBVk2aYqkkzrFl1qmGOuAMAAABwiyPuAAAAAHDNAgsAAACANAssAAAAANIssB6aV8pONUzVpBmmatIMa1adapiqSTNM1aQZpmrSDFM1aYapmjTDVE2aYc2qUw1zxB0AAACAWxxxBwAAAIBrFlgAAAAApFlgAQAAAJBmgfXQvFJ2qmGqJs0wVZNmWLPqVMNUTZphqibNMFWTZpiqSTNM1aQZpmrSDGtWnWqYI+4AAAAA3OKIOwAAAABcs8ACAAAAIM0CCwAAAIC0bwusN2e6ft74xru973v+gJES45OxU2K8P3ZKjE/GTonx/tgpMT4ZOyXG+2OnxPhk7JQY74+dEuOTsVNivD/eecGrOh1xBwAAAODHHHEHAAAAgGsWWAAAAACkWWABAAAAkGaB9fD8LNkLNcNUTZphqibNsGbVqYapmjTDVE2aYaomzTBVk2aYqkkzTNWkGdasOtUwR9wBAAAAuMURdwAAAAC4ZoEFAAAAQJoFFgAAAABpFlgPzStlpxqmatIMUzVphjWrTjVM1aQZpmrSDFM1aYapmjTDVE2aYaomzbBm1amGVY64f/j747/O53/88m4xAAAAANz3wiPu//v93uqpL6urL5ssaywAAAAAvqossL74fo11bLIAAAAAKCywPvz905tF1dfRJgsAAACA9BH3z//45cuvD3//9OXXz3uv5pWyUw1TNWmGqZo0w5pVpxqmatIMUzVphqmaNMNUTZphqibNMFWTZliz6lTDXn/EfT165UgWAAAAwEv8pY+4T9soR7IAAAAA/mpev8D6HRzJAgAAAPjreOUC64/ftLLJAgAAAPgf78VH3P+sfdMfP/fevFJ2qmGqJs0wVZNmWLPqVMNUTZphqibNMFWTZpiqSTNM1aQZpmrSDGtWnWrYi4+4f/zbp5/3wJRz7wAAAAB/or/0EfefxLl3AAAAgP8ZXrbA+nje45k0R7IAAAAA/t298gms99wl2WQBAAAA/Jt68RH39/db596bV8pONUzVpBmmatIMa1adapiqSTNM1aQZpmrSDFM1aYapmjTDVE2aYc2qUw175RH3j+fT58+JZ6CcewcAAAD4IUfcX8m5dwAAAICy1yyw3ueC+8SRLAAAAICmF34LYXQ9ZJMFAAAAkPLiI+4fPjz79Z5+faXst869v7Ps+bRXJ1xoVp1qmKpJM6xZdaphqibNMFWTZpiqSTNM1aQZpmrSDFM1aYY1q0417GVH3D+eTx/PL798fmypPn/+l3XV5+9+P8UzWQAAAMBf1l/9iPuvH7b6rdXVm8egvjwe9eRvfnOg/Xe85utbfHka602GTRYAAADAz5ZYYP3an/IE1p3t0u/bQNlkAQAAALyb4gLrz7p+9eHvn364V/qDl61ssgAAAAB+thcfcf/+GavPn9+OP/QOh9W/vsXzO+4/7+J79nzaqxMuNKtONUzVpBnWrDrVMFWTZpiqSTNM1aQZpmrSDFM1aYapmjTDmlWnGvaaI+4fz6dzzvdH3H/LT0r79SGt+//2d7yLZ7IAAACA/wH+ckfcP57HTufLiurNdxH+1JXa9/up31otfX8b/s7nEH/rL3nzjjZZAAAAAL9D9AbWz9hh3Xm06uu66vuvJvzdO6w3b2STBQAAAPA7FBdY7+D59ur7l/0pO6w3b2qTBQAAAHDfi4+4/5YPH/607yJ8/IU/+uTgD1dUf+Jd9unie/Z82qsTLjSrTjVM1aQZ1qw61TBVk2aYqkkzTNWkGaZq0gxTNWmGqZo0w5pVpxr2siPuX25gfb319eYG1ld/Ytrzh55+a3v19fe//PGf98DUnctcAAAAAC/0lzvift+d3dZvvebry/7g9uodfP9GPmAIAAAA8L0XL7AuPyf4dRV1uZD6rY8W/vAjh398e/U+Ky2nsgAAAAC+13oC68t+6skq6v5hrCerrm9/2+3t1ddr7u/p66cX3+Gt7cgAAACArBcvsB4f8fvnpunXnxP8fRurO8ezCp8c/LFP//nlfz//7f/8pHf48PH/nXM+fPj4/W++52U0AAAAgOde/C2E3z9v9etd1ZN/9Uff91dbqi9n9p9stf7kgoYPH//f54//8eXX9Aez30rw6oRrzTBVk2ZYs+pUw1RNmmGqJs0wVZNmmKpJM0zVpBmmatIMa1adatjrv4Xw/iNXTx6w+vUR968Pdn1/xP3721K/4x7Wz/4iwksfPnw6fzvnnPNff+Zf+/njbz7P9eWBLE9gAQAAAG/4FsJvvm5O3myyvh8vl1yXL3j8w9++/+dP579++brb+vZHbmyv3s238+2ff8q+7PH3/9djfLLPAgAAAHi5xALr5oNX519vZr25n3X5iNa//LV/+/YI1dvt2L1nsn7241fv9rWDb/7+Dx8+fd1hff74H18ewgIAAACISCyw7p+4+v7JqR8+ovXtGav/+uXL9urLs1dv/85XX3N/t73Vb/n8+ZcPf//08+7EAwAAAPwR3464vznT9fPG5/7Inz3/eh7r4s/+7dObf/t1S/Xmxd9vr379+cEfRt75D/Lh75++/Pr4t/P5H798fbt3+z/i+/HzP3758PG/v44f/v7pJRnvOXZKjPfHTonxydgpMd4fOyXGJ2OnxHh/7JQYn4ydEuP9sVNifDJ2Soz3xzsveFVn7oj7V+s3Dz75HOLDPx/C+vb6G89e/emPR738eavf8v1DWB8+/vdPOr8FAAAA/Jv66x5xf/Ldgr/1sj/qb5++7LCmTw7+8WVTdm8FAAAAEJe4gfXV19tVb/7hua+brx89TPbPLxP87pr724Cnnxz8HeytAAAAAP6gFyywPp5fvn6K8NeX13+qz//45etG6ddPWn39nTerq9+xe/p33Vv98n/Pp/98dQQAAADAv/hfP37Jz/fmEaq33yF49conv/OD9/puo/T1jPqXX19/5/sXTxuor3/Vlz/4p2yvnt9L+xN9v937oXermjSrTjVM1aQZ1qw61TBVk2aYqkkzTNWkGaZq0gxTNWmGqZo0w5pVpxr2miPu55933H/51UcFf30S64cfKvzd+T9c1txcP71ZeP3Omob/397dHbd1ZY0W3bz1BdNh2FF0Z0MyGzsKOwxnw/sAmAKgAwiLP8CUOEapurglEJxS+2nVPguve9wtcQcAAABOfN0l7juvM6nDKdXDw8Zw6vvHDA9fdm4v+6aTq1jf/+YPmn/ShwQBAAAAfjatTyF85xqsN8+SzK0AAAAAsu5/A+vwmtX3nyd49ScMfrrd6MrcCgAAAODGEkvcC85tKTvZy37jqhVen3bvhA3NqlUNUzXSDGtWrWqYqpFmmKqRZpiqkWaYqpFmmKqRZpiqkWZYs2pVw+65xH2tVV4T/jUfFbTEHQAAADjnyy1xf1qPuxnWRxmtb7/8PrsvvtTcCgAAAKDs/juwPsQ7503mVgAAAABZv8gA623MrQAAAAD6vuIS95O97LvpVXNL2aqGqRpphqkaaYY1q1Y1TNVIM0zVSDNM1UgzTNVIM0zVSDNM1UgzrFm1qmF3W+K+dnvcf7vp1Sf3rX7IEncAAADgnC+3xH19wh73c8ytAAAAAH5qv+wOLHMrAAAAgF/DLzjA2o2uzK0AAAAAfg13XuL+ek/qQ97qcDX79NubW8pWNUzVSDNM1UgzrFm1qmGqRpphqkaaYapGmmGqRpphqkaaYapGmmHNqlUNu+cS97XW48vazZve854eFfxAlrgDAAAA53zFJe7vZG4FAAAA8EX8ZAMscysAAACAr+bOA6yXvx6veYrQ3AoAAADgy7rzEvfLTvayf+r0qrmlbFXDVI00w1SNNMOaVasapmqkGaZqpBmmaqQZpmqkGaZqpBmmaqQZ1qxa1bD7L3Ffa2OPu/tW92KJOwAAAHDOl17ifvgU4W50ZW4FAAAAwKv7D7B2jK4AAAAA2JQYYJlbAQAAAHBOeon7LTW3lK1qmKqRZpiqkWZYs2pVw1SNNMNUjTTDVI00w1SNNMNUjTTDVI00w5pVqxqWWOJOhyXuAAAAwDl3HOy4gQUAAABAmgEWAAAAAGkGWAAAAACkGWDtNbeUrWqYqpFmmKqRZlizalXDVI00w1SNNMNUjTTDVI00w1SNNMNUjTTDmlWrGmaJO0cscQcAAADOscQdAAAAALYZYAEAAACQZoAFAAAAQJoB1l5zS9mqhqkaaYapGmmGNatWNUzVSDNM1UgzTNVIM0zVSDNM1UgzTNVIM6xZtaphlrhzxBJ3AAAA4BxL3AEAAABgmwEWAAAAAGkGWAAAAACkGWDtNbeUrWqYqpFmmKqRZlizalXDVI00w1SNNMNUjTTDVI00w1SNNMNUjTTDmlWrGmaJO0cscQcAAADOscQdAAAAALYZYAEAAACQZoAFAAAAQNq3AdbJmq7PO5642c+95V8wUuJ44dgpcbz+2ClxvHDslDhef+yUOF44dkocrz92ShwvHDsljtcfOyWOF46dEsfrj9e84F6dlrjzjSXuAAAAwDmWuAMAAADANgMsAAAAANIMsAAAAABIM8Dau7yW7I6aYapGmmGqRpphzapVDVM10gxTNdIMUzXSDFM10gxTNdIMUzXSDGtWrWqYJe4cscQdAAAAOMcSdwAAAADYZoAFAAAAQJoBFgAAAABpBlh7zS1lqxqmaqQZpmqkGdasWtUwVSPNMFUjzTBVI80wVSPNMFUjzTBVI82wZtWqhlnizhFL3AEAAIBzLHEHAAAAgG0GWAAAAACkGWABAAAAkGaAtdfcUraqYapGmmGqRpphzapVDVM10gxTNdIMUzXSDFM10gxTNdIMUzXSDGtWrWqYJe4cscQdAAAAOMcSdwAAAADYZoAFAAAAQJoBFgAAAABpBlh7zS1lqxqmaqQZpmqkGdasWtUwVSPNMFUjzTBVI80wVSPNMFUjzTBVI82wZtWqhlnizhFL3AEAAIBzLHEHAAAAgG0GWAAAAACkGWABAAAAkGaAtdfcUraqYapGmmGqRpphzapVDVM10gxTNdIMUzXSDFM10gxTNdIMUzXSDGtWrWqYJe4cscQdAAAAOMcSdwAAAADYZoAFAAAAQJoBFgAAAABpBlh7zS1lqxqmaqQZpmqkGdasWtUwVSPNMFUjzTBVI80wVSPNMFUjzTBVI82wZtWqhlnizhFL3AEAAIBzLHEHAAAAgG0GWAAAAACkGWABAAAAkGaAtdfcUraqYapGmmGqRpphzapVDVM10gxTNdIMUzXSDFM10gxTNdIMUzXSDGtWrWqYJe4cscQdAAAAOMcSdwAAAADYZoAFAAAAQJoBFgAAAABpBlh7zS1lqxqmaqQZpmqkGdasWtUwVSPNMFUjzTBVI80wVSPNMFUjzTBVI82wZtWqhlnizhFL3AEAAIBzLHEHAAAAgG0GWAAAAACkGWABAAAAkGaAtdfcUraqYapGmmGqRpphzapVDVM10gxTNdIMUzXSDFM10gxTNdIMUzXSDGtWrWqYJe4cscQdAAAAOMcSdwAAAADYZoAFAAAAQJoBFgAAAABp3wZYJ2u6Pu944mY/95Z/wUiJ44Vjp8Tx+mOnxPHCsVPieP2xU+J44dgpcbz+2ClxvHDslDhef+yUOF44dkocrz9e84J7dVrizjeWuAMAAADnWOIOAAAAANsMsAAAAABIM8ACAAAAIM0Aa+/yWrI7aoapGmmGqRpphjWrVjVM1UgzTNVIM0zVSDNM1UgzTNVIM0zVSDOsWbWqYZa4c8QSdwAAAOAcS9wBAAAAYJsBFgAAAABpBlgAAAAApBlg7TW3lK1qmKqRZpiqkWZYs2pVw1SNNMNUjTTDVI00w1SNNMNUjTTDVI00w5pVqxpmiTtHLHEHAAAAzrHEHQAAAAC2GWABAAAAkGaABQAAAECaAdZec0vZqoapGmmGqRpphjWrVjVM1UgzTNVIM0zVSDNM1UgzTNVIM0zVSDOsWbWqYZa4c8QSdwAAAOAcS9wBAAAAYJsBFgAAAABpBlgAAAAApBlg7TW3lK1qmKqRZpiqkWZYs2pVw1SNNMNUjTTDVI00w1SNNMNUjTTDVI00w5pVqxpmiTtHLHEHAAAAzrHEHQAAAAC2GWABAAAAkGaABQAAAECaAdZec0vZqoapGmmGqRpphjWrVjVM1UgzTNVIM0zVSDNM1UgzTNVIM0zVSDOsWbWqYZa4c8QSdwAAAOAcS9wBAAAAYJsBFgAAAABpBlgAAAAApBlg7TW3lK1qmKqRZpiqkWZYs2pVw1SNNMNUjTTDVI00w1SNNMNUjTTDVI00w5pVqxpmiTtHLHEHAAAAzrHEHQAAAAC2GWABAAAAkGaABQAAAECaAdZec0vZqoapGmmGqRpphjWrVjVM1UgzTNVIM0zVSDNM1UgzTNVIM0zVSDOsWbWqYZa4c8QSdwAAAOAcS9wBAAAAYJsBFgAAAABpBlgAAAAApBlgAQAAAJBmgLXXXLO/qmGqRpphqkaaYc2qVQ1TNdIMUzXSDFM10gxTNdIMUzXSDFM10gxrVq1qmE8h5IhPIQQAAADO8SmEAAAAALDNAAsAAACANAMsAAAAANIMsPaaW8pWNUzVSDNM1UgzrFm1qmGqRpphqkaaYapGmmGqRpphqkaaYapGmmHNqlUNs8SdI5a4AwAAAOdY4g4AAAAA2wywAAAAAEgzwAIAAAAg7dsA62RN1+cdT9zs597yLxgpcbxw7JQ4Xn/slDheOHZKHK8/dkocLxw7JY7XHzsljheOnRLH64+dEscLx06J4/XHa15wr05L3PnGEncAAADgHEvcAQAAAGCbARYAAAAAaQZYAAAAAKQZYO1dXkt2R80wVSPNMFUjzbBm1aqGqRpphqkaaYapGmmGqRpphqkaaYapGmmGNatWNcwSd45Y4g4AAACcY4k7AAAAAGwzwAIAAAAgzQALAAAAgDQDrL3mlrJVDVM10gxTNdIMa1atapiqkWaYqpFmmKqRZpiqkWaYqpFmmKqRZlizalXDLHHniCXuAAAAwDmWuAMAAADANgMsAAAAANIMsAAAAABIM8Daa24pW9UwVSPNMFUjzbBm1aqGqRpphqkaaYapGmmGqRpphqkaaYapGmmGNatWNcwSd45Y4g4AAACcY4k7AAAAAGwzwAIAAAAgzQALAAAAgDQDrL3mlrJVDVM10gxTNdIMa1atapiqkWaYqpFmmKqRZpiqkWaYqpFmmKqRZlizalXDLHHniCXuAAAAwDmWuAMAAADANgMsAAAAANIMsAAAAABIM8Daa24pW9UwVSPNMFUjzbBm1aqGqRpphqkaaYapGmmGqRpphqkaaYapGmmGNatWNcwSd45Y4g4AAACcY4k7AAAAAGwzwAIAAAAgzQALAAAAgDQDrL3mlrJVDVM10gxTNdIMa1atapiqkWaYqpFmmKqRZpiqkWaYqpFmmKqRZlizalXDLHHniCXuAAAAwDmWuAMAAADANgMsAAAAANIMsAAAAABIM8Daa24pW9UwVSPNMFUjzbBm1aqGqRpphqkaaYapGmmGqRpphqkaaYapGmmGNatWNcwSd45Y4g4AAACcY4k7AAAAAGwzwAIAAAAgzQALAAAAgDQDrL3mlrJVDVM10gxTNdIMa1atapiqkWaYqpFmmKqRZpiqkWaYqpFmmKqRZlizalXDLHHniCXuAAAAwDmWuAMAAADANgMsAAAAANIMsAAAAABIM8Daa24pW9UwVSPNMFUjzbBm1aqGqRpphqkaaYapGmmGqRpphqkaaYapGmmGNatWNcwSd45Y4g4AAACcY4k7AAAAAGwzwAIAAAAgzQALAAAAgLRvA6yTNV2fdzxxs597y79gpMTxwrFT4nj9sVPieOHYKXG8/tgpcbxw7JQ4Xn/slDheOHZKHK8/dkocLxw7JY7XH695wb06LXHnG0vcAQAAgHMscQcAAACAbQZYAAAAAKQZYAEAAACQZoC1d3kt2R01w1SNNMNUjTTDmlWrGqZqpBmmaqQZpmqkGaZqpBmmaqQZpmqkGdasWtUwS9w5Yok7AAAAcI4l7gAAAACwzQALAAAAgDQDLAAAAADSDLD2mlvKVjVM1UgzTNVIM6xZtaphqkaaYEnutgAAEnJJREFUYapGmmGqRpphqkaaYapGmmGqRpphzapVDbPEnSOWuAMAAADnWOIOAAAAANsMsAAAAABIM8ACAAAAIM0Aa6+5pWxVw1SNNMNUjTTDmlWrGqZqpBmmaqQZpmqkGaZqpBmmaqQZpmqkGdasWtUwS9w5Yok7AAAAcI4l7gAAAACwzQALAAAAgDQDLAAAAADSDLD2mlvKVjVM1UgzTNVIM6xZtaphqkaaYapGmmGqRpphqkaaYapGmmGqRpphzapVDbPEnSOWuAMAAADnWOIOAAAAANsMsAAAAABIM8ACAAAAIM0Aa6+5pWxVw1SNNMNUjTTDmlWrGqZqpBmmaqQZpmqkGaZqpBmmaqQZpmqkGdasWtUwS9w5Yok7AAAAcI4l7gAAAACwzQALAAAAgDQDLAAAAADSDLD2mlvKVjVM1UgzTNVIM6xZtaphqkaaYapGmmGqRpphqkaaYapGmmGqRpphzapVDbPEnSOWuAMAAADnWOIOAAAAANsMsAAAAABIM8ACAAAAIM0Aa6+5pWxVw1SNNMNUjTTDmlWrGqZqpBmmaqQZpmqkGaZqpBmmaqQZpmqkGdasWtUwS9w5Yok7AAAAcI4l7gAAAACwzQALAAAAgDQDLAAAAADSDLD2mlvKVjVM1UgzTNVIM6xZtaphqkaaYapGmmGqRpphqkaaYapGmmGqRpphzapVDbPEnSOWuAMAAADnWOIOAAAAANsMsAAAAABIM8ACAAAAIM0Aa6+5pWxVw1SNNMNUjTTDmlWrGqZqpBmmaqQZpmqkGaZqpBmmaqQZpmqkGdasWtUwS9w5Yok7AAAAcI4l7gAAAACwzQALAAAAgDQDLAAAAADSvg2wTtZ0fd7xxM1+7i3/gpESxwvHTonj9cdOieOFY6fE8fpjp8TxwrFT4nj9sVPieOHYKXG8/tgpcbxw7JQ4Xn+85gX36rTEnW8scQcAAADOscQdAAAAALYZYAEAAACQZoAFAAAAQJoB1t7ltWR31AxTNdIMUzXSDGtWrWqYqpFmmKqRZpiqkWaYqpFmmKqRZpiqkWZYs2pVwyxx54gl7gAAAMA5lrgDAAAAwDYDLAAAAADSDLAAAAAASDPA2mtuKVvVMFUjzTBVI82wZtWqhqkaaYapGmmGqRpphqkaaYapGmmGqRpphjWrVjXMEneOWOIOAAAAnGOJOwAAAABsM8ACAAAAIM0ACwAAAIA0A6y95payVQ1TNdIMUzXSDGtWrWqYqpFmmKqRZpiqkWaYqpFmmKqRZpiqkWZYs2pVwyxx54gl7gAAAMA5lrgDAAAAwDYDLAAAAADSDLAAAAAASDPA2mtuKVvVMFUjzTBVI82wZtWqhqkaaYapGmmGqRpphqkaaYapGmmGqRpphjWrVjXMEneOWOIOAAAAnGOJOwAAAABsM8ACAAAAIM0ACwAAAIA0A6y95payVQ1TNdIMUzXSDGtWrWqYqpFmmKqRZpiqkWaYqpFmmKqRZpiqkWZYs2pVwyxx54gl7gAAAMA5lrgDAAAAwDYDLAAAAADS/u/eAWut9fD7/gHLl788tgYAAADAkcQAa/07utpNsoyxAAAAAHjVeoTw5a/Hl78eH35/3v265Y9urtlf1TBVI80wVSPNsGbVqoapGmmGqRpphqkaaYapGmmGqRpphqkaaYY1q1Y1LPEphA+/P5+7cuXRwhvzKYQAAADAOXf8FMLKI4SbXudWJlkAAAAAX9adB1hXPid4MskyxgIAAAD4Ou5/A2s0jTrc9T79XgAAAAB+Rq0l7lfa7Xr/2HXvzS1lqxqmaqQZpmqkGdasWtUwVSPNMFUjzTBVI80wVSPNMFUjzTBVI82wZtWqht15ifvTbx/2PKA7WR/CEncAAADgnK+7xP2j5k3WvQMAAAD8qu6/A+tjWfcOAAAA8Iv51QZYr6x7BwAAAPg13G2J+9O6xVaw69e9N7eUrWqYqpFmmKqRZlizalXDVI00w1SNNMNUjTTDVI00w1SNNMNUjTTDmlWrGnbPJe5P6/kua8LdybrAEncAAADgnK+7xP32rHsHAAAA+Ll8uQHWK+veAQAAAH4KX3eA9cq6dwAAAICy+yxxv80G95Gn39aV695vLLs+7d4JG5pVqxqmaqQZ1qxa1TBVI80wVSPNMFUjzTBVI80wVSPNMFUjzbBm1aqG3W2J+9N6flqPt9z19QZf8E6WJe4AAADAOV93ifvDw6U/veFsbeunW/cOAAAAEJDYgfXycjTJOjnenUkWAAAAwB0lBljfj6vO3b06WU21W1l14Z1PFrS/4TWvP2K3IeskwyQLAAAA4LN9W+J+sqbr847XeHjYT7Uuf+/lP3367cc/6G0TqNON7w+nGTf7x3R887FT4nj9sVPieOHYKXG8/tgpcbxw7JQ4Xn/slDheOHZKHK8/dkocLxw7JY7XH695wb0677zEfboD6zNuYF12eAPrwst+jTtZlrgDAAAA53zdJe6H665246qT42XvmUxd6fVHXB5jeboQAAAA4JPcYYD1tL5NnQ5vYJ3cxnp42JhhvW0qdPJd31/juv5Pr/9ZJlkAAAAAH+I+N7Ce1n6mc3jr6vU21qc+1Hg4n7pwnerw4tU7p2YmWQAAAADv8f9+/JKbu7wY6+1v+/vzhenVbj3Y67hq8zMH32b3Vt82vg/f7fK+tHtRNdIMUzXSDGtWrWqYqpFmmKqRZpiqkWaYqpFmmKqRZpiqkWZYs2pVw+6zxH23wX0d7Po6uYG1PucS1g/vXm1etvqM+1PZO1mWuAMAAADnfN0l7ud8+LOEb5tejV5wPU8XAgAAAFwvOsD6PKPp1evvf9LHHW4+pWiYBQAAAHCoPsD6/unCncPLWede8/qyyxedfji9uoHDH+RaFgAAAMChOy9xf3g4O3VaW9vcLzxU+MPV75enVCdbyi7vw/pUp0vfH57ftv39U2WXut07YVszTNVIM6xZtaphqkaaYapGmmGqRpphqkaaYapGmmGqRpphzapVDbv/EvfDqdP+ttTWHOqaG1ib33J4A2tzIDXd5n77W1GHMSczrA+PscQdAAAAOOfrLnE/GVp9P6X64b2qzRdfM5QrPDn4Y89/roO50u5/dx6e/vfw8HT42lvOIgEAAABu5s4DrMP51PezqnODrQ/4ucPpVerxvZ2Xpz/WWg9P/7t3CAAAAMDnuvMOrBMXrhAd/tE7bxq9YXp1rztZD0//vH6x+3XygpenP3aTLAAAAIBfVWuAtdZ6edn/OnFyV+vw17kXHH7x7bt+f15/P27c9ro4vbrxDazXre0vL49Pv62Xvx5fXva/Hv7+59ww65ayS93unbCtGaZqpBnWrFrVMFUjzTBVI80wVSPNMFUjzTBVI80wVSPNsGbVqoZVlri/HHzx6uSPdg53Zp3sz9rcgbX/nd++/evvplSv+933L/vR3avbbHB/HZNd84MeHp7XWi9P+61YD0//e///lZa4AwAAAOd83SXuO9evuDocPJ181/frtL7Nc/5+3M+w5nevri17h9Hc6tXLy63vhQEAAADcRWKA9YF+8MGFvz2vv4+GRNdMrz5pTvS2udXpmzz98+0SVurzEwEAAAA+SH2AdeGjCU/sRldnP9bw9RLW4VtNplcfNRv6wDc8XM718vTHfbdiAQAAAHySOy9xv/KzBT9yT9e/Y6yT6dXrlrLNqdb7h03f9rL/9bj7deU3Zten3TthQ7NqVcNUjTTDmlWrGqZqpBmmaqQZpmqkGaZqpBmmaqQZpmqkGdasWtWw+y9xX1tb2C98sen6v8Hhw4CjJwffPMD68Atc37//bu36+ojN65a4AwAAAOd8rSXuT+vxdYZ19om/z3H4zN33N61OPnPw8LumP+iz51YAAAAAX0diB9bJTavTzxDceuXh78x+1vEM6+RP3zm6us/c6vG/6/nP2/04AAAAgNtKDLBubPOa1bmX/dCHPGkIAAAAwDl3XuK+83qp6vDu1eYThZc/lPCHM6lDr8vUD6dOm7+53fzdUvZPml5l16fdO2FDs2pVw1SNNMOaVasapmqkGaZqpBmmaqQZpmqkGaZqpBmmaqQZ1qxa1bD7LHFf/+5xv/0S9/erLbfa7+16/nNZ4g4AAAB8pq+1xP3E4fjp9esLX9zLbnQVmVsBAAAAfB33H2DF1a5cAQAAAHw1dxtgPa3nx9UdCZlbAQAAAETcZ4n700ePrkbr2zfttpR9v5r9I+o+IKxG1UgzTNVIM6xZtaphqkaaYapGmmGqRpphqkaaYapGmmGqRpphzapVDbvbEve11tN67qwJ/3nvW1niDgAAANzGl17ifkc/79wKAAAA4Ov4igMscysAAACAn8g9B1j7x99u+ON2X5hbAQAAAPxE7rPEfX3CHvdzrtzL3txStqphqkaaYapGmmHNqlUNUzXSDFM10gxTNdIMUzXSDFM10gxTNdIMa1atatg9l7ivtZ5++8QbWF/kvpUl7gAAAMBtWOL+kXajq197bgUAAADwddx5gPWBa7C+yJUrAAAAgK/mngOsl78eX6dOb2ZuBQAAAPBru9sS93e6cjX79ZpbylY1TNVIM0zVSDOsWbWqYapGmmGqRpphqkaaYapGmmGqRpphqkaaYc2qVQ278xL3x5fxI4TuW52wxB0AAAC4ja+7xH33FOEPp1HmVgAAAABfVvpTCM2tAAAAACgOsMytAAAAAHh1/yXuh59F+LF72UeaW8pWNUzVSDNM1UgzrFm1qmGqRpphqkaaYapGmmGqRpphqkaaYapGmmHNqlUNu/8S9x1Xrt7MEncAAADgNr7uEvcdcysAAAAAzrnnAOtkbsfbPD+stf67+/L9/6TPf3/YWwEAAAB8iPvvwAIAAACAC74NsE7WdH3eEQAAAICf3c1GSXde4s77Pa3np/X4/dfvfLf3vxUAAADwa7vlEnePEAIAAACQdtMbWHyGh9+fX377z1rr4emfl5e3X5t6+P355a/H9fzn+98KAAAA4AO5gbWxouv5+dneLgAAAICIOwywUrOh15jDLx4fHx8fHyOd12e8PP3n4ff3NT//ee0LG/84J5pVqxqmaqQZ1qxa1TBVI80wVSPNMFUjzTBVI80wVSPNMFUjzbBm1aqG3XqJe81uVrX738PfOXnBD9/h++Pr/9+vg7DL7/Mer48QrrUe/v7n5a83/qDd+zw8/e/l6Q+PEAIAAAAd/3fvgF/W6yTr5ItPshs8rb/f+z4vT3/s3woAAACg4evuwLowUWrel7vgdeT05qcID69xrbXW438/qg0AAADgnb7oAOvyfahPvSr14U6fGfx7TWdY+88fBAAAAEj6ukvcn//1+vW9i7ZdGfbw9L/dFy9P/3n57T8Pvz/vfp19/b8v2N+9+nd9++v7fEjVjTWrVjVM1UgzrFm1qmGqRpphqkaaYapGmmGqRpphqkaaYapGmmHNqlUN++pL3HcuLGL/4VWskwXw6/zqq0/dgfXw8LDW2lxc9fD3P5vfsr9y9d3HDu4GWP6rAAAAADoscT81/dzA7z9w8F72s6fjMdbRZqtDz3+ug+VZ6+q7VwAAAAA35gbWu3z2ZwteaXcD69DoYwS/H135rwIAAADocAPrV3Ayb3r4/fnhafvJwe1vf/rDxw4CAAAAWV93ifuJZtV6U9jLX48vL4Nfb5heNf+5mlWrGqZqpBnWrFrVMFUjzTBVI80wVSPNMFUjzTBVI80wVSPNsGbVqoZZ4g4AAABA3R1uYAEAAADA9QywAAAAAEh71xL35oORZYWPLAQAAAD4ufyan0L4hjnR8/PzZ0+X3jbvu0HYG6gaaYapGmmGNatWNUzVSDNM1UgzTNVIM0zVSDNM1UgzTNVIM6xZtaph713inr2BFfy3XtX/CAAAAADi7MACAAAAIM0ACwAAAIA0AywAAAAA0n7NJe7Z5Vxv0NycpWqkGaZqpBnWrFrVMFUjzTBVI80wVSPNMFUjzTBVI80wVSPNsGbVqoa9d4k7AAAAAHw2jxACAAAAkGaABQAAAECaARYAAAAAaXcYYDU3rDerVjVM1UgzTNVIM6xZtaphqkaaYapGmmGqRpphqkaaYapGmmGqRpphzapVDbPEHQAAAIA6jxACAAAAkGaABQAAAECaARYAAAAAaZa47zWrVjVM1UgzTNVIM6xZtaphqkaaYapGmmGqRpphqkaaYapGmmGqRpphzapVDbPEHQAAAIA6jxACAAAAkGaABQAAAECaARYAAAAAaZa47zWrVjVM1UgzTNVIM6xZtaphqkaaYapGmmGqRpphqkaaYapGmmGqRpphzapVDbPEHQAAAIC6/w/PxJYKDpgTsAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 549,
       "width": 800
      }
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Message from klive: {\"version\": \"0.2.2\", \"type\": \"open\", \"file\": \"/var/folders/2_/hfctv7kn7h965wdl4jgrm8t80000gn/T/tmp9z98hcki/Top.gds\"}\n"
     ]
    }
   ],
   "source": [
    "from SiEPIC.utils import load_Waveguides_by_Tech\n",
    "waveguide_types = load_Waveguides_by_Tech(tech_name)\n",
    "\n",
    "print('Waveguide types:')\n",
    "for w in waveguide_types:\n",
    "    print('  %s' % w['name'])\n",
    "wg_type = waveguide_types[0]['name']\n",
    "\n",
    "from SiEPIC.scripts import connect_pins_with_waveguide\n",
    "wg = connect_pins_with_waveguide(inst_gc1, 'opt1', inst_gc2, 'opt1', \n",
    "                                    turtle_A=[5, -90, \n",
    "                                              20, 90, 10, 90, 1, 90],\n",
    "                                    waveguide_type=wg_type,\n",
    "                                )\n",
    "\n",
    "from SiEPIC.utils import waveguide_length\n",
    "print(waveguide_length(wg))\n",
    "\n",
    "topcell.plot() # in the browser\n",
    "topcell.show() # in KLayout GUI\n",
    "        "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9648c84c",
   "metadata": {},
   "source": [
    "## Developers\n",
    "- Install the tools and PDKs by creating a copy using GitHub Desktop, then installing as a symbolic link in pip."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ff650ee7",
   "metadata": {},
   "outputs": [],
   "source": [
    "!python -m pip install -e '/Users/lukasc/Documents/GitHub/SiEPIC-Tools/klayout_dot_config/python'\n",
    "!python -m pip install -e '/Users/lukasc/Documents/GitHub/SiEPIC_EBeam_PDK/klayout'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "36268e1c",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.10.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}