create_cascadedomain2d.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "outputs": [],
   "source": [
    "import pyvista as pv\n",
    "import numpy as np\n",
    "from ntrfc.cascade_case.domain.geoparas import DomainParameters\n",
    "from ntrfc.cascade_case.domain.domain import CascadeDomain2D\n",
    "from ntrfc.turbo.airfoil_generators.naca_airfoil_creator import naca"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "outputs": [],
   "source": [
    "pv.set_jupyter_backend(\"static\")"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "outputs": [],
   "source": [
    "xs,ys = naca(\"6510\",256)\n",
    "points = pv.PolyData(np.stack([xs,ys,np.zeros(len(xs))]).T)\n",
    "alpha = 1"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/many/miniconda3/envs/NTRfC/lib/python3.10/site-packages/pyvista/core/filters/poly_data.py:2772: PyvistaFutureWarning: The default value of the ``capping`` keyword argument will change in a future version to ``True`` to match the behavior of VTK. We recommend passing the keyword explicitly to prevent future surprises.\n",
      "  warnings.warn(\n"
     ]
    }
   ],
   "source": [
    "domainparas = DomainParameters()\n",
    "domainparas.generate_params_by_pointcloud(points,alpha=1)\n",
    "domainparas.xinlet=-3\n",
    "domainparas.xoutlet=4\n",
    "domainparas.pitch=2\n",
    "domainparas.blade_yshift=0.1"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "\u001B[0m\u001B[33m2022-12-19 12:14:12.538 (   6.212s) [         510C180]        vtkRenderer.cxx:1172  WARN| vtkOpenGLRenderer (0x455c5f0): Resetting view-up since view plane normal is parallel\u001B[0m\n"
     ]
    },
    {
     "data": {
      "text/plain": "<PIL.Image.Image image mode=RGB size=1024x768>",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABAAAAAMACAIAAAA12IJaAAA/7ElEQVR4nO3deZRcZZ0//ud2AolEICEsAgJBIWGpAkQdxGFAfgiigwt+AWkEYVRUZDA4AmIcBVEyAnFcDgMKDkfOaJd4lEEyMICHRYhgQGWpkiVmBAIJsnRYAiRk6fr9UUmn6e50V3c/3ffWva/X8YTqqlu3nq661X7eVZ/nuckRRxwRAACAYmhLewAAAMDYEQAAAKBABAAAACgQAQAAAApEAAAAgAIRAAAAoEAEAAAAKBABAAAACkQAAACAAhEAAACgQAQAAAAoEAEAAAAKRAAAAIACEQAAAKBABAAAACgQAQAAAApEAAAAgAIRAAAAoEAEAAAAKBABAAAACkQAAACAAhEAAACgQAQAAAAoEAEAAAAKRAAAAIACEQAAAKBABAAAACgQAQAAAApEAAAAgAIRAAAAoEAEAAAAKBABAAAACkQAAACAAhEAAACgQAQAAAAoEAEAAAAKRAAAAIACEQAAAKBABAAAACiQ8WkPAMihPbfb+96rHw5JkrQlISRhXFtIkpAkSVtbaEtCaAvjkhDaQltIkrbQloSksVkSQhKSpL7uQuOmepKEJITGrtpCSNrqSQhtoZ40bg2hLQlJCCHUk6TeFtZe2bhXEupJqLeFkIT62ruHegj1tqRxYd01ybptQkhCj80au133aUkS6knvf9deaNy6brN+t+zeLCQhJPUeN9XX3WXthZ6bNf5NknrjKUkaVzZ+y8Zzk4Qk6Wpc3/i1kqTe1niSknpI6uPWPkn1trauxjPX1tbV2HJc0tWWrH3axiVdSdLVloRxoastqbcl9XFJV1uot4WuJAnjkzVtSWgLXeOTriTU25L6+GRNW6i3hfr4ZE2S1NtCfVzSNS50jQv1JOnaKOlKQldbUh8futpC17gQkqRro9CVJF1toT4+qa+9b6i3hdCWhPEhtCVJWwhtIYxP2pKQtIVkXGgbl7T9/MZ/rFarY3kAA+SbAADEt9dBe5TK5UqlEkIIIVlX1CZJEtaWt0nj33UXwrofQ7Ku/F1/TX3d7T1uXb9Nj1vX36seemxWD6Gr8d+kcaH7pnr3Hbuved3e1v633ueafjZucvv+9t/jjvXeN4Xu56bea/tk3e8U1j1z3Y+cdF9ed991T+3aISfrntR6Ul8TQhLCmlBf/foNkvWBqBGdum9a+zK0hXHrrhm/Lt00wlp3tGkkrLU76b4p9AhKa2/tPkbWJr/1R8le5b1KpZIAABCRAACMirYt1rRtsSbKrpLBNxlMvc8FQgjrno9hPStJlFdmUH9e8H9z535vXZgEIAJzAID4KpXKcccdl/YoyI9yuZz2EADyQwAARkW1WlW0EUVHR0epVEp7FAD5IQAAo6JarSraACCDBABgVNRqNd8AEIWOMoC4BABgVGgBIiKHE0BEAgAwWjo6Otrb29MeBXmgowwgIgEAgKzTUQYQkQAAjJZaraZ1myi0AAFEJAAAo8XZW4lIRxlALAIAMIoUbQCQNQIAAC1ARxlALAIAMIqs4E4sOsoAYhEAgFFn+iZR6CgDiEIAAEZXR0eHFdwBIDsEAABag44ygCgEAGB0KdqIS0cZwAgJAMCocxYnYtFRBjByAgAw6qrVqqINADJCAABGXa1W8w0AUegoAxg5AQAYdVqAiMjhBDBCAgAwFqzgTiw6ygBGSAAAoJXoKAMYIQEAGAu1Wk3rNlFoAQIYIQEAGAvVajXtIZAfOsoARkIAAMaIog0AskAAAKDF6CgDGAkBABgjVnAnFh1lACMhAABjyvRNotBRBjBsAgAwdjo6OqzgDgDpEgAAaD06ygCGTQAAxo6ijbh0lAEMgwAAjClncSIWHWUAwyMAAGOqWq0q2gAgRQIAMKZqtZpvAIhCRxnA8AgAwJjSAkREDieAYRAAgLFmBXdi0VEGMAwCAACtSkcZwDAIAMBYq9VqWreJQgsQwDAIAMBYq1araQ+B/NBRBjBUAgCQAkUbAKRFAACghekoAxgqAQBIgRXciUVHGcBQCQBAakzfJAodZQBDIgAA6ejo6LCCOwCMPQEAgNamowxgSAQAIB2KNuLSUQbQJAEASI2zOBGLjjKA5gkAQGqq1aqiDQDGmAAApKZWq/kGgCh0lAE0TwAAUqMFiIgcTgBNEgCANFnBnVh0lAE0SQAAIA90lAE0SQAA0lSr1bRuE4UWIIAmCQBAmqrVatpDID90lAE0QwAAUqZoA4CxJAAAkBM6ygCaIQAAKbOCO7HoKANohgAAZILpm0ShowxgUAIAkL6Ojg4ruAPA2BAAAMgPHWUAgxIAgPQp2ohLRxnAAAQAIBOcxYlYdJQBDEwAADKhWq0q2oiiVqulPQSATBMAgEyo1Wq+ASCKarWqowxgAAIAkAlagIjI4QQwAAEAyAoruBOLjjKAAQgAAOSNjjKAAQgAQFbUajWt20ShBQhgAAIAkBXVajXtIZAfOsoANkQAADJE0QYAo00AACCHdJQBbIgAAGRIpVJRtBGFjjKADREAgMwxfZModJQB9EsAALKlo6PDCu4AMHoEAADySUcZQL8EACBbFG3EpaMMoBcBAMgcZ3EiFh1lAH0JAEDmVKtVRRsAjBIBAMicWq3mGwCi0FEG0JcAAGSOFiAicjgB9CIAAFlkBXdi0VEG0IsAAECe6SgD6EUAALKoVqtp3SYKLUAAvQgAQBZVq9W0h0B+6CgD6EkAADJK0QYAo0EAACDndJQB9CQAABllBXdi0VEG0JMAAGSa6ZtEoaMMoJsAAGRXR0eHFdwBIC4BAID801EG0E0AALJL0UZcOsoAggAAZJyzOBGLjjKABgEAyLRqtapoIxZhEiAIAEDG1Wo1RRtRVCoVxxJAEACAjNMCREQOJ4AgAADZZwV3YtFRBhAEAACKQ0cZQBAAgOyr1WoWAyUKLUAAQQAAsq9araY9BPJDRxmAAAC0AEUbAMQiAABQIDrKAAQAoAVUKhVFG1HoKAMQAICWYfomUegoAwpOAABaQ0dHhxXcAWDkBAAAikVHGVBwAgDQGhRtxKWjDCgsAQBoGc7iRCw6yoAiEwCAllGtVhVtxCJMAoUlAAAto1arKdqIolKpOJaAwhIAgJahBYiIHE5AYQkAQCuxgjux6CgDCksAAKCIdJQBhSUAAK2kVqtZDJQotAABhSUAAK2kWq2mPQTyQ0cZUEwCANBiFG0AMBICAAAFpaMMKCYBAGgxlUpF0UYUOsqAYhIAgJZk+iZR6CgDCkgAAFpPR0eHFdwBYHgEAACKS0cZUEACANB6FG3EpaMMKBQBAGhJzuJELDrKgKIRAICWVK1WFW3EIkwChSIAAC2pVqsp2oiiUqk4loBCEQCAlqQFiIgcTkChCABAq7KCO7HoKAMKRQAAoOh0lAGFIgAArapWq1kMlCi0AAGFIgAAraparaY9BPJDRxlQHAIA0MIUbQAwVAIAAOgoAwpEAABaWKVSUbQRhY4yoDgEAKDlmb5JFDrKgIIQAIDW1tHRYQV3AGieAAAAIegoAwpDAABam6KNuHSUAbknAAAtz1mciEVHGVAEAgDQ8qrVqqINAJokAAAtr1ar+QaAKHSUAUUgAAAtTwsQETmcgNwTAIA8sII7segoA3JPAACA9XSUAbknAAB5UKvVtG4ThRYgIPcEACAPqtVq2kMgP3SUAfkmAAA5oWgDgGYIAADwOjrKgHwTAICcsII7segoA/JNAAByxfRNotBRBuSYAADkR0dHhxXcAWBgAgAA9KajDMgxAQDID0UbcekoA3JJAAByxVmciEVHGZBXAgCQK9VqVdEGAAMQAIBcqdVqvgEgCh1lQF4JAECuaAEiIocTkEsCAJA3VnAnFh1lQC4JAADQPx1lQC4JAEDe1Go1rdtEoQUIyCUBAMibarWa9hDIDx1lQP4IAEAOKdoAYEMEAADYIB1lQP4IAEAOWcGdWHSUAfkjAAC5ZfomUegoA3JGAADyqaOjwwruANCXAAAAA9FRBuSMAADkk6KNiJwQAMgTAQDILUUbsVSrVR1lQG4IAEBuKdoAoC8BAMitWq3mGwCi0FEG5IkAAOSWFiAicjgBuSEAAHlmBXdi0VEG5IYAAACD01EG5IYAAORZrVbTuk0UWoCA3BAAgDyrVqtpD4H80FEG5IMAAOScog0AehIAAKApOsqAfBAAgJyzgjux6CgD8kEAAArB9E2i0FEG5IAAAORfR0eHFdwBoEEAAIBm6SgDckAAAPJP0UZcOsqAliYAAIXgLE7EoqMMaHUCAFAI1WpV0UYswiTQ0gQAoBBqtZqijSgqlYpjCWhpAgBQCFqAiMjhBLQ0AQAoCiu4E4uOMqClCQAAMDQ6yoCWJgAARVGr1SwGShRagICWJgAARVGtVtMeAvmhowxoXQIAUCCKNgAQAABgyHSUAa1LAAAKpFKpKNqIQkcZ0LoEAKBwTN8kCh1lQIsSAIBi6ejosII7AEUmAADAcOgoA1qUAAAUi6KNuHSUAS1HAAAKx1mciEVHGdCKBACgcKrVqqKNWIRJoOUIAEDh1Go1RRtRVCoVxxLQcgQAoHC0ABGRwwloOQIAUERWcCcWHWVAyxEAAGD4dJQBLUcAAIqoVqtZDJQotAABLUcAAIqoWq2mPQTyQ0cZ0FoEAKCgFG0AFJMAAAAjoqMMaC0CAFBQlUpF0UYUOsqA1iIAAIVm+iZR6CgDWogAABRXR0eHFdwBKBoBAABGSkcZ0EIEAKC4FG3EpaMMaAkCAFBozuJELDrKgFYhAACFVq1WFW3EIkwCLUEAAAqtVqsp2oiiUqk4loCWIAAAhaYFiIgcTkBLEACAorOCO7HoKANaggAAAHHoKANaggAAFF2tVrMYKFFoAQJaggAAFF21Wk17COSHjjIg+wQAAEUbAAUiAABANDrKgOwTAABCpVJRtBGFjjIg+wQAgLVM3yQKHWVAxgkAACGE0NHRYQV3AIpAAACAmHSUARknAACEoGgjNh1lQGYJAABrOYsTsegoA7JMAABYq1qtKtqIRZgEMksAAFirVqsp2oiiUqk4loDMEgAA1tICREQOJyCzBACA9azgTiw6yoDMEgAAID4dZUBmCQAA69VqNYuBEoUWICCzBACA9arVatpDID90lAHZJAAAvI6iDYB8EwAAYFToKAOySQAAeJ1KpaJoIwodZUA2CQAA/TB9kyh0lAEZJAAA9NbR0WEFdwDySgAAgNGiowzIIAEAoDdFG3HpKAMyRQAA6IezOBGLjjIgawQAgH5Uq1VFG7EIk0CmCAAA/ajVaoo2oqhUKo4lIFMEAIB+aAEiIocTkCkCAED/rOBOLDrKgEwRAABgdOkoAzJFAADoX61WsxgoUWgBAjJFAADoX7VaTXsI5IeOMiA7BACADVK0AZA/AgAAjDodZUB2CAAAG1SpVBRtRKGjDMgOAQBgEKZvEoWOMiAjBACAgXR0dFjBHYA8EQAAYCzoKAMyQgAAGIiijbh0lAGpEwAABuEsTsSiowzIAgEAYBDValXRRizCJJA6AQBgELVaTdFGFJVKxbEEpE4AABiEFiAicjgBqRMAAAZnBXdi0VEGpE4AAICxo6MMSJ0AADC4Wq1mMVCi0AIEpE4AABhctVpNewjkh44yIF0CAEBTFG0A5IMAAABjSkcZkC4BAKAplUpF0UYUOsqAdAkAAENg+iZR6CgDUiQAADSro6PDCu4AtDoBAADGmo4yIEUCAECzFG1E5IQAQFoEAIAhULQRS7Va1VEGpEIAABgCRRsRCZNAKgQAgCGo1WqKNqKoVCqOJSAVAgDAEGgBIiKHE5AKAQBgaKzgTiw6yoBUCAAAkA4dZUAqBACAoanVahYDJQotQEAqBACAoalWq2kPgfzQUQaMPQEAYMgUbQC0LgEAAFKjowwYewIAwJBVKhVFG1HoKAPGngAAMEymbxKFjjJgjAkAAMPR0dFhBXcAWpEAAABp0lEGjDEBAGA4FG1E5IQAwFgSAACGSdFGLNVqVUcZMGYEAIBhUrQRkTAJjBkBAGCYarWaoo0oKpWKYwkYM+PTHgBAq6pWq7Nnz057FClbsKD9nPCNXldWpk9PZTAtrdFR5rQAwBhIjjjiiLTHANCqGsu3VyqVtAeSmv/5n7n1kPS6sjp16qz9909lPK3LsQSMGS1AAAxTZ2f/XSvlzs72BQvGeDCtTkcZMGYEAIDhq9VqRV4MtLNzg3Ogj1uwQAYYEotKAWNGAAAYvoJ3bG/oGwCGp1qtNhqBAEaVAAAwIh0dHQUv2s4N5/R7/XELFpQ7O8d4MC2t4HkSGDMCAAAj8o1wbhLq/d5UEgCGwjQAYGwIAAAjomhr6Pd7gONMAxgK0wCAsSEAAIyIoq3hG+HcfjOAqcBD4nACxoAAADBSiraGm/Zf2fdK0wCGpFqtlkobXFsJIAoBAGCkFG3dOvqcA1gAAMgaAQCA0aULqHmmlABjQAAAGClFW0NnZ6kyfXp16tRe1/sSoHnayYAxIAAAjJSirScBYIQcTsBoEwAA4lC0NVT6TAMIMsBQmFICjDYBACCCjo4ORVu3vlOBnREMIDsEAACGaerUas8fOzs3+B2IbwCaZ0oJMNoEAAAiq5kGMALmAACjTQAAiMCntqHHNwB95wEHGQAgMwQAgAh8atuLaQAj4XACRpUAABBN0Yq2qVNraQ8hnywEBIwqAQAgjmq1OvhGedfdBWQaAEBmCQAAcfjUtienAwPILAEAgGHqtQxoCKGzc30E6ncqMM0wpxwYVQIAQByKtl58CTBsJgEDo0oAAIhG0TYwCwEBZIEAABBHMScB9+0C6lbpsxIozfMlADB6BAAAouleBQiAzBIAAKIp4Ke2A3wD0Jc5AM2zqBQwegQAgJgKXrT1+gbAQkAAGSQAAERTzGkAA+gVAHwDAJAFAgAAwzd1ai3tIeSTVWWB0SMAABCTecAAGScAAERTwE9tB54EXDMHYLgKOKEcGDMCAAAxdXYONA3aNACA1AkAANH41LYXqwABZJAAAMCIDOlUADRPngRGiQAAQEwmAQNknAAAEFMBP7Ud+BsAXUAAWSMAAEAWVavVgp9YGhglAgAAMQ3cAlSyChBA2gQAAAAoEAEAIKYCtm1MnVpLewgADIEAAEBkFgICyDIBAIARcR6AUVKr1Yq2ohQwNgQAAAAoEAEAAAAKRAAAILLOzmJNgwZoLQIAACPV/DSAmhMDA6RNAAAAgAIRAAAAoEAEAABGaoAWoHJn51iOBIBBCQAAROZEYABZJgAAAECBCAAAjBb9PwAZJAAAMHaqlgEFSJsAAMBITZ1aS3sIADRLAAAAgAIRAACIrHsVoJI5AADZIwAAMFp6TQI2AQAgCwQAAAAoEAEAIKZyuVyrmRFLBKVSqVrd4CmWAYZNAABgpKZO7b9O1QIEkEECAACjou9ZwGoCAEAGCAAAAFAgAgAAZJH5JMAoEQAAYiqXyyZuNrQvWNDrGnMAALJAAAAAgAIRAAAYCz7+B8gIAQAgGv0/PfVdBYghcTgBo0QAACC+vtW/bwAAMkIAAGAsOAnAkPj4Hxg9AgBANKVSSdHW0HcJIAAyQgAAYKQ6O8uDbqMFCCAjBACAaJy5qVuvOQCq/6HybRIwegQAACJ7T7gt7SEAsEECAACRHdQnAFSmT09jIC3Mt0nA6BEAAKKxckvDe8Jv0x4CABskAADEofrv1rcFyByAoXI4AaNHAAAgJtU/QMYJAACMVGdnKe0h5IqP/4FRJQAAxGHdxoZzwjd6XWMGMECmCAAAcVi2pcEaoCMnTAKjSgAAIJp+q39zAIZKmARGlQAAEEeR+7Y7O8uNC33PAKD6B8gaAQAggiJX/z31PQOAADAMDidgVAkAAETTtwWoJgAMUblcTnsIQM4JAAARFHzWZqMFyASAWIp8LAFjQAAAiMCszdDfAqCq/2EoeJgExoAAAMBocQaAYRAmgdEmAABEUORZm939P84AANASBACAkSpy9d+t3/4fLUDD4HACRpsAABCBio0oVP/AGBAAAEaqvb29yE3bnZ2lfvt/TAAAyCYBAICR6tv/EywBNCzt7e2VSiXtUQA5JwAAjFTB2zY6O8t9P/5X/QNklgAAMCIFr/7DBs7/pf9neBxOwBgQAABGxGmb9P/EovoHxoYAADAiTttU7uzsdY3qf3iESWBsCAAAI+JT2770/wyPMAmMDQEAYPhU//3yDcDwOJyAsSEAABCT6n94VP/AmBEAAIbPqu196f8ByLjxaQ8AgNbW8fqK3zcAwyNMAmNGAAAYpnK5HELQtuEj/yjK5fKsWbPSHgVQCFqAACBlJgAAY0kAABgmPRsAtCIBAABSJkwCY0kAABiOcrmsbQOAViQAAAyT6p8ozCYHxpgAADAcejYAaFECAACkSZgExpgAADAcJgAA0KIEAIAhU/0TkcMJGGMCAMCQlUolFRtRqP6BsScAAAxZuVyu1Wppj4I8MAEAGHsCAMCQ+dQWgNYlAAAMjeqfiBxOwNgTAACGplQqpT0EckL1D6RCAAAYmnK5rGmbKMwmB1IhAAAMjU9ticVsciAVAgDAEKj+icjhBKRCAACAFKj+gbQIAABDYNV2YjGbHEiLAAAwBOVyOe0hkBNmkwNpEQAAmtXo2dC2QRRagIC0CAAAMNZU/0CKBACAZpkAAEAOCAAAMNaESSBFAgBAU8rlsrYNYjGbHEiRAADQLNU/UZhNDqRLAABoip4NAPJBAACAMSVMAukSAACaYgIAAPkgAAAMTvVPLGaTA6kTAAAG197ermIjFscSkC4BAKAptVot7SGQByYAAKkTAACa4lNbAPJBAAAYhHM2EZEJAEDqBACAQZRKJRUbUaj+gSwQAAAGUS6XTQAgCrPJgSwQAAAG4VNbIhImgdQJAAADUf0TkcMJyAIBAGAgpVIp7SGQE6p/ICMEAICBlMtlq7YThdnkQEYIAAAD8aktsZhNDmSEAACwQap/InI4ARkhAADAqFP9A9khAABsUHt7uwkAAOSMAACwQeVyOe0hkBPCJJAdAgBA/xo9G9o2iEILEJAdAgAAjC7VP5ApAgBA//RsAJBLAgAAjC5hEsgUAQCgH+VyWdsGALkkAAD0T/VPFI21pBxOQHYIAAD90LMBQF4JAAAwioRJIGsEAIB+mAAAQF4JAAC9qf6JxWxyIIMEAIDe2tvbVWzE4lgCskYAAOhHrVZLewjkgQkAQAYJAAC96dkAIMcEAIDXUf0TkcMJyCABAOB1SqWSio0oVP9ANgkAAK9TLpdNACAKYRLIJgEA4HV8aksswiSQTQIAwHqqfyJyOAHZJAAAQHyqfyCzBACA9azaTiylUintIQD0TwAAWK9cLqc9BHKiXC4Lk0A2CQAAazV6NrRtEIUWICCzBAAAiEz1D2SZAACwlgkAABSBAAAAkQmTQJYJAAAhrJv+q22DKMwmB7JMAACAmMwmBzJOAAAIQc8GAIUhAABATMIkkHECAEAI1m0EoDAEAADVP9GUy2WHE5BxAgCAng1iUv0DGScAAEA0wiSQfQIAgBYgAApEAACKTvVPRA4nIPsEAKDoSqWSio0oVP9ASxAAgKIrl8u1Wi3tUZAH7e3tAgCQfQIAUHQ+tSUiYRLIPgEAKDTVPxE5nICWIAAAhVYqldIeAjmh+gdahQAAFFq5XLZqO1GYTQ60CgEAKDSf2hKL2eRAqxAAgOJS/RORwwloFQIAAIyU6h9oIQIAUFzt7e0mABCF2eRACxEAgOIql8tpD4GcMJscaCECAFBQjZ4NbRtEoQUIaCECAACMiOofaC0CAFBQJgAAUEwCAACMiDAJtBYBACiicrmsbYNYzCYHWosAABSU6p8ozCYHWo4AABSRng0ACksAAIDhEyaBliMAAEVkAgAAhSUAAIWj+icWs8mBViQAAIXT3t6uYiMWxxLQcgQAoIhqtVraQyAPTAAAWpEAABSOng0AikwAAIpF9U9EDiegFQkAQLGUSiUVG1Go/oEWJQAAxVIul00AIAqzyYEWJQAAxeJTWyISJoFWJAAABaL6JyKHE9CiBAAAGDLVP9C6BACgQKzaTixmkwOtSwAACsSntsRiNjnQugQAoChU/0TkcAJalwAAAEOj+gdamgAAFIUJAAAQBAAAGCphEmhpAgBQCOVyOYSgbYMoGocTQIsSAABgCBoTAIRJoHUJAEAh6NkAgAYBAACGQJgEWp0AABSCdRsBoEEAAPJP9U8sZpMDOSAAAPmnZwMAugkAANAsYRLIAQEAyD8tQADQTQAAck71T0QOJyAHBAAg50qlkoqNKFT/QD4IAEDOlcvlWq2W9ijIg/b2dgEAyAEBAMg5n9oSkTAJ5IAAAOSZ6p+IHE5APggAQJ6VSqW0h0BOqP6B3BAAgDwrl8tWbScKs8mB3BAAgDzzqS2xmE0O5IYAAOSW6p+IHE5AbggAADAI1T+QJwIAkFvt7e0mABCF2eRAnggAQG6Vy+W0h0BOmE0O5IkAAORTo2dD2wZRaAEC8kQAAICBqP6BnBEAgHwyAQAA+iUAAMBAhEkgZwQAIIfK5bK2DWIxmxzIGQEAyCfVP1GYTQ7kjwAA5JCeDQDYEAEAADZImATyRwAAcsgEAADYEAEAyBvVP7GYTQ7kkgAA5E17e7uKjVgcS0D+CABADtVqtbSHQB6YAADkkgAA5I2eDQAYgAAA5Irqn4gcTkAuCQBArpRKJRUbUaj+gbwSAIBcKZfLJgAQhdnkQF4JAECu+NSWiIRJIJcEACA/VP9E5HAC8koAAPKjVCqlPQRyQvUP5JgAAORHuVy2ajtRmE0O5JgAAOSHT22JxWxyIMcEACAnVP9E5HACckwAAIDXUf0D+SYAADnR3t5uAgAADEoAAIDXESaBfBMAgDwol8shBG0bRKEFCMg3AQAA1lP9A7knAAB5oGcDAJokAADAesIkkHsCANDyyuWytg0AaJIAAOSB6p8ozCYHikAAAFqeng0AaJ4AAABrCZNAEQgAQMszAQAAmicAAK1N9U9EDiegCAQAoLWVSiUVG1Go/oGCEACA1lYul2u1WtqjIA9MAAAKQgAAWptPbQFgSAQAoIWp/onI4QQUhAAAtDATAIhF9Q8UhwAAtDATAIhFmASKQwAAWphPbYlFmASKQwAAWpXqn4gcTkBxCAAAFJ3qHygUAQBoVVZtJ5ZSqZT2EADGjgAAtKpyuZz2EMiJcrksTALFIQAALanRs6Ftgyi0AAGFIgAAUGiqf6BoBACgJZkAAADDIwAAUGjCJFA0AgDQehrTf7VtEIXZ5EDRCAAAFJfZ5EABCQBA69GzAQDDJgAAUFzCJFBAAgDQeqzbCADDJgAALUb1TyzlctnhBBSQAAC0mPb2dhUbsTiWgAISAIDWU6vV0h4CeWACAFBMAgDQYvRsAMBICABAK1H9E5HDCSgmAQBoJaVSScVGFKp/oLAEAKCVlMtlEwCIwmxyoLAEAKCV+NSWiIRJoJgEAKBlqP6JyOEEFJYAALSMUqmU9hDICdU/UGQCANAyyuWyVduJwmxyoMgEAKBl+NSWWMwmB4pMAABag+qfiBxOQJEJAAAUi+ofKDgBAGgN7e3tJgAQhdnkQMEJAEBrKJfLaQ+BnDCbHCg4AQBoAY2eDW0bRKEFCCg4AQCAAlH9AwgAQAswAQAAYhEAACgQYRJAAACyrlwua9sAgFgEAKAFqP6JorGWlMMJKDgBAMg6PRsAEJEAAEBRCJMAQQAAss8EAACISAAAMk31TyxmkwM0CABAprW3t6vYiMWxBBAEACD7arVa2kMgD0wAAGgQAIBM07MBAHEJAEB2qf6JyOEE0CAAANlVKpVUbESh+gfoJgAA2VUul00AIAphEqCbAABkl09tiUWYBOgmAAAZpfonIocTQDcBAICcU/0D9CQAABll1XZiKZVKaQ8BIEMEACCjfGpLLOVyWZgE6CYAAFmk+icihxNATwIAAHmm+gfoRQAAssgEAAAYJePTHgAAjKKxCZPHHHPMCSec0Li8evXqlStXdnZ2Lly48Oabb37ggQfq9XrPjWfPnr1mzZqvfe1rA+/z8MMPP/XUU4877rhly5aN1riBQhIAgMwpl8shBG0bRDGWM4AvvPDCEML48eMnTpy49dZbv/Od7zz44IMXL1787W9/+7HHHhubMQAMSgAAILcaEwDGLEzecccdPX+88sorZ8yYMXPmzIsuuuhf//VfH3nkkcb1v/3tb7u6usZmSAB9mQMAZI4JAOTGI488csYZZ3R2dp5xxhkTJ05sXDl//vx77rkn3YEBReYbAAByKwth8tVXX/3+979/4YUXHnjggTfddFMI4ayzzuqeAzBlypRPf/rTb3vb21577bWbb765UqmsWbPmggsu2GOPPRp37+joCCH88pe/vPLKKxvXTJ069Zhjjtl111133HHH5cuXP/jgg1ddddWzzz7b0dHxgx/84De/+U1YN83gggsuOPHEE/fcc88pU6YsXLjwzjvv/N///d9ew9tss82OPvroGTNmTJs27eWXX3700Ufnzp1733339dpsyy23POaYY6ZPn/7mN7/5+eeff+CBB379618vWrRoFJ84YNQIAEDmlMvlWbNmpT0KiOahhx5atGjRP/zDPzQCQLckSb761a9uu+22119//eTJk4855pg1a9ZUKpWOjo7NNtvsbW9726GHHnrxxRe/+uqrTzzxROMuBx544Kmnnrpq1ao777zz1ltvnTp16r777nvRRRf99Kc/7fWgW2211b//+7+/+OKL8+bNe+211/bcc8/Pf/7ze+2110UXXdTdgLTffvvNnDmzq6vrrrvuuv322zfZZJN3vvOd3/zmN2+66aZLLrlkzZo1jc0OOuigU089ddmyZXfcccdNN9201VZb/f3f//3BBx/8wx/+sNdvBLQEAQDIFqu2E0umZpM/8sgj++yzT68rt9tuuxkzZlx66aXXX399COENb3hD44P/+++/P4QwadKkQw899M477+xeBWjHHXc8/fTT77333u9973vdV1555ZVHH330Jz/5yV4733777X/961//53/+Z2MNol/96leHHXbYaaeddvfdd996662NDc4444xqtfrd7363e2+/+MUvTjzxxKOOOmrZsmU/+clPQgjTpk2bOXPmbbfddskll6xevbqx2c9+9rNPfepT//zP//zMM8/0/boAyDgBAMiWLPRsQHQvvPDC5MmT+71pypQpjQsXXXTRwDs5+eSTly5dOmfOnOXLl3dfWa/Xf/GLX+y8884HHHBAz42XLVt25ZVX9lyB9KabbjriiCMOO+ywRgD4zGc+s2zZsjlz5rz66qs973jllVcuWrRo8uTJ48aNW7Nmzcknn7x48eKLL76458Tl1atXX3bZZbvsssvHP/5xAQBajgAAQD5lKkxuttlmL730Uq8rFy9efNdddx177LE77LDDddddV6vVep0xoKckSXbbbbe5c+f2rP67/eY3v+kVAJ566qlVq1b12uyxxx7be++9G3vbY489rrvuul7Vf0MjIYQQ2traZsyYMX/+/He96119N3v88cff9773bbrpps5UAK1FABiRxllaGpeffvrpT3/606PxKEceeeQJJ5xw7LHHrly5ctq0ad/73vfOP/98K0iQVyYAkEu77bbbU0891ff6Cy644H3ve9/73//+2bNnP/LIIxdddNHTTz/d7x622mqriRMnLlmypN9bFy9e3OuafnPCiy++OGHChBDClltuOXHixCeffHLgYW+55ZYTJkw48MADDzzwwA1tM2XKFAEAWosAMCIPPvjgv/3bv22xxRaf+cxn7rzzzlF6lPnz53/yk58slUr33nvvqaeeOm/ePNU/eWUCABFlJ0zus88+O+200w9/+MO+N61Zs+b666+//vrrS6XSF7/4xfPOO++zn/1svzt57rnnXnvttR122KHfW9/85jf3uqbfLxO6r3zuuedWrFix/fbbDzzyZ599dsWKFR0dHf/93/898JZAC3EegBFZtGjRnXfeuc8++zz11FN9V2CIZcmSJU888cQ+++xz2GGHbbvttpdddtkoPRCkrlQqCQBEkZ0wOWXKlC996UudnZ233XZbr5tmzJix//77Ny7XarWf/exn22233Zve9KZ+99PV1bVgwYIDDzxw0qRJfW897LDDhjSqer3+5z//+cADD9xkk0363loulw866KBx48bV6/UHH3ywcbnvZqeccsonPvGJIT0ukAUCwEgdcsghf/d3f/fd73535cqVo/cojf7LE0888Uc/+lHfLlLIjXK5XKvV0h4FedDe3p6FAFAqlc4555xJkybNmTPnlVde6XXr3nvvPWvWrO4P9TfffPMQQt/Nul1++eWTJ0/+0pe+tOmmm3ZfmSTJUUcd1R0kmnfZZZdtttlmZ5xxRs+9hRAOOeSQ888/v1QqNZYBvfzyy3fYYYfTTjtt44037vmgJ5100uGHH95YsAhoLVqARmTq1Kknn3zy1Vdf/fDDD4/qA/3+978/6qij5s+f3+s885Az2enZIAfGPky++93vDiGMHz9+woQJ22yzzV577bX77rs/++yzZ5999oIFC/puf8stt3z0ox8966yzbrzxxs033/wjH/nI7bffPkA//aOPPvr973//85///KWXXjpv3rxFixZNmTLl7W9/+1vf+tbLL7/8M5/5zJBGu2TJkgsvvHDmzJk//OEP582b9/jjj2+11ValUmm33XabN2/ej370o8ZmTz755He+850vfOELe+2117x585YsWTJ58uQDDzywcZ4BAQBakQAwIqeddtrSpUt/9rOfjfYDPf/88yGEX//616P9QJCi7PRskAOphMmvfOUrIYR6vf7SSy89//zzS5YsmT179j333NO9fH4vzz333KxZs0466aTjjz9++fLlN95446DdpLfddlu1Wj366KN33XXXQw45ZNWqVQ8//PBZZ531xBNP7L///o3/s2jePffcc8oppxx11FG77777QQcdtHz58scff/zrX//6vffe23OzO++88+GHHz766KNLpdL73//+Z5999uGHHz7nnHOeffbZIT0ckBHJEUcckfYYWtWhhx76hS984aqrrvrLX/4SQpg/f/7oPdY73vGOc84557jjjrPSAjnW3t7uGwCiKJfL7e3tjiWAfvkGYPje8573hBA+9rGPhRCee+65UQ0AO+2009KlS1X/5Fu5XM7Oqu20NLPJAQbgGwAgK+bOnfvBD34w7VGQB7Nnz65UKjIAQL8EgCFYsGDfzs7tOju3nTr1qenT/zh1aj+ndIllWdeKtgnjJ63yFQ2ZM2nSpJ///OeNy3/84x/PPffcgbf/7Gc/2+vvzJNPPnnKKad0/7jtttt+/OMfL5VKkydPfuihh2677babbrppgPOhNmnHHXc89thj3/rWt2622WYLFy78wx/+cO211/bc7cYbb9ze3v7ud7978uTJjz76aKVSyfF0xkGfjV4GfdUG3eFQHzGuVgmTQ303Zfxpb3XR/7h5OUZi5513PvbYY/fYY4+NNtrosccemzt37u9+97sBtvfuGBL1ZbMWLNh3wYK3Ny53dm57111HTJ/+x+nT/zR6j/jioVu/8tTyzR9dOeGFNaP3KDBUK1asOOecc8K6/rdBTZs27e67777uuut67qH78j777POv//qvjz766Lx58/bdd99HH330c5/73Lve9a7zzjtvJH+X3/Wud5155pl//etfb7jhhs7Ozl122eWEE07Yc88958yZ01ixN0mSr371q9OmTfvVr361dOnS/fff/5vf/ObXv/71++67b9gPmlmDPht9DfyqDbrDYTxiRC00m3xI76aMP+05EPePm5djJEql0rnnnrto0aKrr766s7Nzr732+vKXv3zNNddcccUV/W7v3TFUAkBTOju37a7+uy1Y8PbRCwAvdS1v22jcxvvvsGzP5U8v6pzwwprJi1ZLAmTBmjVr/vSnP4UQDjvssIkTJw66/U477fTzn/+8cZdeJk6cOHPmzPnz58+ZM+f888+/9NJLq9Xqb3/72wsuuOAjH/nIsM88uvHGG59yyil33333nDlzGguZ33777bfccst3vvOdI4444uqrrw4hlEqlfffd91/+5V8ak/jnzZu3ySabHHfccfkLAM08G30N8KoNusPhPWJEpVJptB8ilubfTdl/2nMg4h83L8dIJEny+c9/fuHChV//+tcb1fntt9/+0EMPzZw584477mj80e7Ju2MYBIARefDNH3zr554ejT1P/L9n1ixftdFmb9hoszdM3GbzEMKrLy1f/ezKiX96ftyil0fjEWE0bLHFFptuuumiRYv6vXWvvfbacsstr7jiinq93j0D+JFHHrnlllsOPvjgngFg1113PfLII3fZZZcpU6YsXrz41ltvve666za0tOJuu+22xRZbdHR0NP7QNzz22GO/+93vDjrooMbf+p122mn16tULFy7s3uDPf/7zRz/60Si/daY082z0MvCrNugOh/GIceVyNnn2n/aiGeHbZKyG2ZLe/OY377DDDl/72td6fjZ/8803f+pTn9p77737BgDvjmFwJuCmdHZu2+/1W7xjg+drjKtr1Zpxb9h45Yuvrt528M8kIDumTZsWQli0aNF22203Y8aMXh+q7bTTTsuWLevs7Gz0bHS3bSxcuHDHHXdMkqTx45FHHnnRRRdNmjTpmmuu+Y//+I9arXbCCSdceOGFvU5f2m2rrbYKISxZsqTX9UuWLJk6dWrj8lNPPTV+/Pg999yz+9ZSqfTYY4+N+DfOnGaejV4GftUG3eEwHjGuFmoBal72n/aiGeHbhAG0tbXdc889TzzxRK/rly9fPmXKlL7be3cMg28AmjJ9+p/6tgCFEDb6z4cmXDMqU4GfW/182/+bvrLz5TXLV65ZvnLSg68kC1/YzJxgWs20adPq9fp555230047hRDq9fr8+fMvueSSxumKnn/++Te+8Y0TJkzoda+pU6d2dnY25gCUSqV/+qd/uuyyy/7nf/6ncettt902d+7ciy66aObMmd/61rf6PuikSZNWrFjR85OehldeeeUNb3hD4/Kf/vSnBx54YNasWVdfffUzzzyz//7777HHHrNnz47622dCM89GLwO/aoPucBiPGFEuq/+Q+ae9gEb4NmEAjz/++Hnnndfryre85S1bb711z69tu3l3DINvAJo1ffof+14zqgsBvfb0i22/eWKjeU9v86un3/jQy1YEohVNmzYtSZL77rvvtNNOO/bYY7/xjW9sv/323/nOdzbZZJMQwp///Oeurq4PfOAD7e3t3T0bb3zjG9/73vc+8MADjR8/+clP3n333d3Vf8PTTz998cUX77fffo3/621SV1dX9+V6vX7ttddusskmJ5544plnnnnAAQc8+OCD//d//zfSX7h19Hw2ehn4VRvGDpvcgGHwtKdllN4m9GvixImnn376o48+escddzR/L++OAagpmzV9+p+mTn1qwYK3d3Zu+55dly5JHh3VJYC2Hz8l/L6I09LJmauuuurWW2+99957Gz/+8Y9/XLhw4SWXXHL88cdfdtllTz311K9+9atPfOITS5cu3WKLLSZNmjR9+vRPfepTkyZNuuqqq0IIG2200Vve8pb7779/77337rXnVatWrVixYt9993388ceHMbCDDjroS1/60tVXX33DDTe8+OKLM2bMOPnkk7/97W+feeaZr7766gh/61Y38KuW7tgG1TNMwuhp6bdJa9lkk02+/vWvb7755mefffaGpn4xVALAEEyd+tQPTrjjuF0271j4YmXhi2kPB1rA4sWLFy9e3POaF1988dZbb33729f21P30pz994xvf+IEPfOCMM84IIaxevXr8+PE//vGP//a3v4UQtt5663Hjxh111FFHHXVUv/ufMmXKP/7jP37uc5/rvua+++675557uucP9NTWtvY7z/Hjx3/605++4YYbfvKTn3Tf69xzz7300ksPP/zw/E0IG/jZ6GvQV23QHQ71ESPK5Qzghiw/7QU08rcJzdhss83OO++8zTfffNasWU89tcG2C++OoRIAgNHylre8Zeedd7755pt7Xf/ss882pmSFEOr1+rx583bcccfLL7981apVZ5xxxqpVq6699trGrc8888yaNWsuvPDCO++8c0OPsvnmmz/yyCPdPy5fvny33XabMGHCRhtttGrVqp5bTpo06eWXXw4hbLPNNpMnT77rrrt6jWrhwoW77rrrCH7jLHrllVcGfjZ6GfRVG3SHQ33EiHrNJs+TLD/tBTTytwnNmDp16je/+c1x48Z9+ctffuaZZza0mXfHMBQ3+gCj7U1vetPpp5++zTbb9Lr+rW99a8/1dtrb2zs6Ov7617/ut99+O+644w9+8IPuU4CtWrXqr3/96wEHHNB35zNmzDj77LO33HLLF198cWEPixcvbvz/xA477NDrLttvv/1zzz0XQnjppZdCCH3nfm2yySYvvpi3L/cGfTZ6GfRVG3SHQ31EmuFpz5SRv00Y1Jve9KYLLrhg9erVA1f/wbtjWAQAIKaNN9740EMP3XrrrUMIf/jDH/72t7999rOf7fk16+67737QQQfdcsstve643Xbbtbe3/+IXv+i1rvaPf/zjAw444AMf+EDPK7fccssvfvGLEydO7Pdv94MPPrh06dJjjz123Lhx3VfuvPPO++23X+ObhGXLlv3lL3/54Ac/2PNL4b322mvatGn333//8H/5TBr02QhDfNUG3WEzjzhKcjYBoOfrkuWnvSDivk0Y2A477HDBBRcsXbr0K1/5ygsvvNB3A++OEUqOOOKItMfQStp32dwcAGg4++yzJ06ceO655/a8cvLkyf/1X//1rW99a/78+SGE6dOnn3vuuU8//fSNN9740ksv7b777kccccQf/vCH2bNnd3/MP3fu3A996EPnn3/+5ptvPnPmzL5zvD784Q+fdNJJ1Wr197///QsvvLDLLrscfvjhTz/99Ne+9rUNfXs7Y8aMc84558knn7zjjjs6Ozt33XXXww8//Mknnzz33HNfeeWVEMKOO+44Z86cJ5988sYbb2xMAv7Qhz50zz33fPvb347+RKVu0GdjqK/aoDscdINRMnv27Eql0ootQM28mzL7tOdPlD9uXo5ha1T/K1euvPTSS1esWNHzpmeeeaYxE8C7Y4QEgKERAKBbM/8fGULYeuutP/axj+27776TJk167LHHbrnllhtvvLG7+i+Xy+3t7b/97W9PPfXUM888s2c3f0/Tpk076qijdtllly222OLxxx//3e9+d+211w68gtsOO+xw/PHH77nnnhMmTHj88cer1WpHR0fPBtAtt9zy+OOP32OPPaZMmbJo0aJbbrnl+uuv7x5Yzgz8bAz1VRt0h81sEF25XJ49e/YHP/jB0XuI0dPkuymDT3suRfnjFrwcw3XIIYecfvrp/d70y1/+8sorrwzeHSMmAAyNAABxzZ49u1qt5qltg7Q0wuSsWbPSHghA1pkDAKSsVqulPQTyIGcTAABGjwAApKmxbmPaowCAAhEAgNSo/onI4QTQJAEASE2pVFKxEYXqH6B5zgQMpKZcLg/atL3iuokhhNBYsD9JQkiSJAnJ2svrLoQkSULjf+s2C23rrg+hcTkkbT0u99iyscPujbtPDpAkobHId8+dNPS4vr7u7muv6d5J0mNX6+5SX7ur0H339TtZv+W6HSZJPYSQhO5drd0s6XH3EEII9bb1u+25TX39ozRu6h7nurus3W3osdveO1n3ECEkG9jh6y733mE/l3vsuZ40d5cN76Tx4xMvhut/+Jv6Tv9fr43XD36wXdV73XGAwW/omkF2Uh9oDD2vef2jvP4XqQ90xx4P1OOY7P4F633Htm6beu/r+16z7u5Jnw2S113usdfG8bv21vq6K8P64yvpvVk/16zbVdLn+h7vv3rPO7b1uGNbz2tef7nxivXY/7p3Z4/NGsYna/ruqvvu6y8n9cYbq3u3bUk9hHpbUk+67x4aO+lav8PQlfT4BRu7agtrx5OEevcde14OIYxbu5OudbeGJOlKui93j2fdAZWE+vika91O6ut2Ug8hjAu9d9WWdCUhJI2xrf0V1v4NSNaNJFk7krWjauxkXM9fMNTXj2f9r7bueVv7XK17TkIIIYzr8Qdy3NrnMEnW/RVsaAvJuCRJQmgL6//utyVJCElbSBobhBDGJW0hhA+fvFugP74BAFLjU1sAGHsCAJCO4VT/9cY/9e7Lg2w60O09d7KBjQc9J8DaDZo9dcBA273usfrfsOmH6f/5qffZZGg7GXiHA+2qz78b2GSgXW1gJ82OYYjjaXaHTTxkvdd/RzSeQQ/IZnc1uEF21dQTv/YgGmQ8zQ23x642fPQ2tav+xtP312nuqe57eQMbDe2128DNzR5LTTzOCLdbN9ihHGQbfO83+We2mX2P8KAvFgEAAAAKRAAA0mHVdgBIhQAApMMEAABIhQAApED1DwBpEQBG5OCDD547d+7ee+8dQpg4ceLFF198zjnnJD1XQQMAgCwRAEbk9ttvf/bZZz/84Q+HEE4//fTx48fPmTOn3sSMdSg4EwAAIC0CwIisWbPmmmuuecc73vGFL3zhbW9727e+9a1XXnkl7UEBAMAGCQAjddNNN7388svvfe9758yZ8+STT6Y9HGgB5XI5hGAOAOSbfljILAFgpLbZZpuNNtqoq6vrscceS3ssAAAwCAFgRDbddNOvfvWr999//6pVq4488si0hwOtwQQAAEiRADB8bW1tZ5111muvvTZnzpwbbrjhsMMO22yzzdIeFAAADEQAGL6TTjppl112mT179ooVK6655ppx48Z96EMfSntQkHXlctlJAMiJrPW4Z208QFaNT3sALeyKK6644oorGpc7Ozu1AEGTVP8AkCLfAABjygQAAEiXAAAAAAUiAABjygQAAEiXOQDN+uiapcmq5du98Nptf3hm5d+W/r/VXUsmTL4rmZT2uKCVqP4BIHUCQLMmbD75uG23CdtMDuPHhd3eHF54edaCV8PKtIcFLaVUKv3s4qu6lo4LIYSQhKSxbkmSJGHtCibJuuu7/w31tEYLALkkADSr8nLbcSGENV1h/LiwpissW15dacU1GJoHfvvgvVc/FJKNkrYkhCSMawtJEpIkaWsLbUkIbaEtCUlbaAtJkoS2tpAkISxNe9QAkCsCwBB0PPXacVuGEEJYvaa60geTMGR/XnL/xu8a2l1WPTxx/Q9rvyTo/q5gAIPl86TnTjawcdLcTppefX2gvxqve6z+d9jsH52k/+en592TZLBd9d3JwDsc6Ffr829/mwyyqw3spNkxDHE8Q9vhYI9WH3hnQxhPUm9mREM7KpvYST+7GsITP9i7pHH3wZ7sZP2uNvy+bGpX/Y2n76/T1Hj6udzPeBqXmtvhgLtKelwefFQDPs4It1s32KG8RTb43h/0mRl8yEmv/9IEAWAInluxMrz6WpiwUVi5utLZ5VCDMTDxH1ekPYRCqvf4t6ErnYFA6tb0+DcbGuXHuB7/wtBYBWgIfrN64nNLl4UQqq926f8BAKAVCQBDc9Or40IIty5+Ke2BAADAcAgAQ1Pp7Aqvraq+nKVvAgEAoGkCwJDNemjZ39ompD0KAAAYDgFgyKovWPwfAIBWJQAAAECBCAAAAFAgAgAAABSIAAAAAAUiAAAAQIEIAAAAUCACAAAAFIgAAAAABSIAAABAgQgAAABQIAIAAAAUiAAAAAAFIgAAAECBCAAAAFAgAgAAABSIAAAAAAUiAAAAQIEIAAAAUCACAAAAFIgAAAAABSIAAABAgQgAAABQIAIAAAAUiAAAAAAFIgAAAECBCAAAAFAgAgAAABSIAAAAAAUiAAAAQIEIAAAAUCACAAAAFIgAAAAABSIAAABAgfz/gEWtxQu9piMAAAAASUVORK5CYII=\n"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "domain2d = CascadeDomain2D()\n",
    "domain2d.from_cascade_parameters(domainparas)\n",
    "domain2d.plot()"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "\u001B[0m\u001B[33m2022-12-19 01:52:13.289 (  43.991s) [        B92B1180]        vtkRenderer.cxx:1172  WARN| vtkOpenGLRenderer (0x4209f70): Resetting view-up since view plane normal is parallel\u001B[0m\n"
     ]
    },
    {
     "data": {
      "text/plain": "<PIL.Image.Image image mode=RGB size=1024x768>",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABAAAAAMACAIAAAA12IJaAAA/7ElEQVR4nO3deZRcZZ0//ud2AolEICEsAgJBIWGpAkQdxGFAfgiigwt+AWkEYVRUZDA4AmIcBVEyAnFcDgMKDkfOaJd4lEEyMICHRYhgQGWpkiVmBAIJsnRYAiRk6fr9UUmn6e50V3c/3ffWva/X8YTqqlu3nq661X7eVZ/nuckRRxwRAACAYmhLewAAAMDYEQAAAKBABAAAACgQAQAAAApEAAAAgAIRAAAAoEAEAAAAKBABAAAACkQAAACAAhEAAACgQAQAAAAoEAEAAAAKRAAAAIACEQAAAKBABAAAACgQAQAAAApEAAAAgAIRAAAAoEAEAAAAKBABAAAACkQAAACAAhEAAACgQAQAAAAoEAEAAAAKRAAAAIACEQAAAKBABAAAACgQAQAAAApEAAAAgAIRAAAAoEAEAAAAKBABAAAACkQAAACAAhEAAACgQAQAAAAoEAEAAAAKRAAAAIACEQAAAKBABAAAACiQ8WkPAMihPbfb+96rHw5JkrQlISRhXFtIkpAkSVtbaEtCaAvjkhDaQltIkrbQloSksVkSQhKSpL7uQuOmepKEJITGrtpCSNrqSQhtoZ40bg2hLQlJCCHUk6TeFtZe2bhXEupJqLeFkIT62ruHegj1tqRxYd01ybptQkhCj80au133aUkS6knvf9deaNy6brN+t+zeLCQhJPUeN9XX3WXthZ6bNf5NknrjKUkaVzZ+y8Zzk4Qk6Wpc3/i1kqTe1niSknpI6uPWPkn1trauxjPX1tbV2HJc0tWWrH3axiVdSdLVloRxoastqbcl9XFJV1uot4WuJAnjkzVtSWgLXeOTriTU25L6+GRNW6i3hfr4ZE2S1NtCfVzSNS50jQv1JOnaKOlKQldbUh8futpC17gQkqRro9CVJF1toT4+qa+9b6i3hdCWhPEhtCVJWwhtIYxP2pKQtIVkXGgbl7T9/MZ/rFarY3kAA+SbAADEt9dBe5TK5UqlEkIIIVlX1CZJEtaWt0nj33UXwrofQ7Ku/F1/TX3d7T1uXb9Nj1vX36seemxWD6Gr8d+kcaH7pnr3Hbuved3e1v633ueafjZucvv+9t/jjvXeN4Xu56bea/tk3e8U1j1z3Y+cdF9ed991T+3aISfrntR6Ul8TQhLCmlBf/foNkvWBqBGdum9a+zK0hXHrrhm/Lt00wlp3tGkkrLU76b4p9AhKa2/tPkbWJr/1R8le5b1KpZIAABCRAACMirYt1rRtsSbKrpLBNxlMvc8FQgjrno9hPStJlFdmUH9e8H9z535vXZgEIAJzAID4KpXKcccdl/YoyI9yuZz2EADyQwAARkW1WlW0EUVHR0epVEp7FAD5IQAAo6JarSraACCDBABgVNRqNd8AEIWOMoC4BABgVGgBIiKHE0BEAgAwWjo6Otrb29MeBXmgowwgIgEAgKzTUQYQkQAAjJZaraZ1myi0AAFEJAAAo8XZW4lIRxlALAIAMIoUbQCQNQIAAC1ARxlALAIAMIqs4E4sOsoAYhEAgFFn+iZR6CgDiEIAAEZXR0eHFdwBIDsEAABag44ygCgEAGB0KdqIS0cZwAgJAMCocxYnYtFRBjByAgAw6qrVqqINADJCAABGXa1W8w0AUegoAxg5AQAYdVqAiMjhBDBCAgAwFqzgTiw6ygBGSAAAoJXoKAMYIQEAGAu1Wk3rNlFoAQIYIQEAGAvVajXtIZAfOsoARkIAAMaIog0AskAAAKDF6CgDGAkBABgjVnAnFh1lACMhAABjyvRNotBRBjBsAgAwdjo6OqzgDgDpEgAAaD06ygCGTQAAxo6ijbh0lAEMgwAAjClncSIWHWUAwyMAAGOqWq0q2gAgRQIAMKZqtZpvAIhCRxnA8AgAwJjSAkREDieAYRAAgLFmBXdi0VEGMAwCAACtSkcZwDAIAMBYq9VqWreJQgsQwDAIAMBYq1araQ+B/NBRBjBUAgCQAkUbAKRFAACghekoAxgqAQBIgRXciUVHGcBQCQBAakzfJAodZQBDIgAA6ejo6LCCOwCMPQEAgNamowxgSAQAIB2KNuLSUQbQJAEASI2zOBGLjjKA5gkAQGqq1aqiDQDGmAAApKZWq/kGgCh0lAE0TwAAUqMFiIgcTgBNEgCANFnBnVh0lAE0SQAAIA90lAE0SQAA0lSr1bRuE4UWIIAmCQBAmqrVatpDID90lAE0QwAAUqZoA4CxJAAAkBM6ygCaIQAAKbOCO7HoKANohgAAZILpm0ShowxgUAIAkL6Ojg4ruAPA2BAAAMgPHWUAgxIAgPQp2ohLRxnAAAQAIBOcxYlYdJQBDEwAADKhWq0q2oiiVqulPQSATBMAgEyo1Wq+ASCKarWqowxgAAIAkAlagIjI4QQwAAEAyAoruBOLjjKAAQgAAOSNjjKAAQgAQFbUajWt20ShBQhgAAIAkBXVajXtIZAfOsoANkQAADJE0QYAo00AACCHdJQBbIgAAGRIpVJRtBGFjjKADREAgMwxfZModJQB9EsAALKlo6PDCu4AMHoEAADySUcZQL8EACBbFG3EpaMMoBcBAMgcZ3EiFh1lAH0JAEDmVKtVRRsAjBIBAMicWq3mGwCi0FEG0JcAAGSOFiAicjgB9CIAAFlkBXdi0VEG0IsAAECe6SgD6EUAALKoVqtp3SYKLUAAvQgAQBZVq9W0h0B+6CgD6EkAADJK0QYAo0EAACDndJQB9CQAABllBXdi0VEG0JMAAGSa6ZtEoaMMoJsAAGRXR0eHFdwBIC4BAID801EG0E0AALJL0UZcOsoAggAAZJyzOBGLjjKABgEAyLRqtapoIxZhEiAIAEDG1Wo1RRtRVCoVxxJAEACAjNMCREQOJ4AgAADZZwV3YtFRBhAEAACKQ0cZQBAAgOyr1WoWAyUKLUAAQQAAsq9araY9BPJDRxmAAAC0AEUbAMQiAABQIDrKAAQAoAVUKhVFG1HoKAMQAICWYfomUegoAwpOAABaQ0dHhxXcAWDkBAAAikVHGVBwAgDQGhRtxKWjDCgsAQBoGc7iRCw6yoAiEwCAllGtVhVtxCJMAoUlAAAto1arKdqIolKpOJaAwhIAgJahBYiIHE5AYQkAQCuxgjux6CgDCksAAKCIdJQBhSUAAK2kVqtZDJQotAABhSUAAK2kWq2mPQTyQ0cZUEwCANBiFG0AMBICAAAFpaMMKCYBAGgxlUpF0UYUOsqAYhIAgJZk+iZR6CgDCkgAAFpPR0eHFdwBYHgEAACKS0cZUEACANB6FG3EpaMMKBQBAGhJzuJELDrKgKIRAICWVK1WFW3EIkwChSIAAC2pVqsp2oiiUqk4loBCEQCAlqQFiIgcTkChCABAq7KCO7HoKAMKRQAAoOh0lAGFIgAArapWq1kMlCi0AAGFIgAAraparaY9BPJDRxlQHAIA0MIUbQAwVAIAAOgoAwpEAABaWKVSUbQRhY4yoDgEAKDlmb5JFDrKgIIQAIDW1tHRYQV3AGieAAAAIegoAwpDAABam6KNuHSUAbknAAAtz1mciEVHGVAEAgDQ8qrVqqINAJokAAAtr1ar+QaAKHSUAUUgAAAtTwsQETmcgNwTAIA8sII7segoA3JPAACA9XSUAbknAAB5UKvVtG4ThRYgIPcEACAPqtVq2kMgP3SUAfkmAAA5oWgDgGYIAADwOjrKgHwTAICcsII7segoA/JNAAByxfRNotBRBuSYAADkR0dHhxXcAWBgAgAA9KajDMgxAQDID0UbcekoA3JJAAByxVmciEVHGZBXAgCQK9VqVdEGAAMQAIBcqdVqvgEgCh1lQF4JAECuaAEiIocTkEsCAJA3VnAnFh1lQC4JAADQPx1lQC4JAEDe1Go1rdtEoQUIyCUBAMibarWa9hDIDx1lQP4IAEAOKdoAYEMEAADYIB1lQP4IAEAOWcGdWHSUAfkjAAC5ZfomUegoA3JGAADyqaOjwwruANCXAAAAA9FRBuSMAADkk6KNiJwQAMgTAQDILUUbsVSrVR1lQG4IAEBuKdoAoC8BAMitWq3mGwCi0FEG5IkAAOSWFiAicjgBuSEAAHlmBXdi0VEG5IYAAACD01EG5IYAAORZrVbTuk0UWoCA3BAAgDyrVqtpD4H80FEG5IMAAOScog0AehIAAKApOsqAfBAAgJyzgjux6CgD8kEAAArB9E2i0FEG5IAAAORfR0eHFdwBoEEAAIBm6SgDckAAAPJP0UZcOsqAliYAAIXgLE7EoqMMaHUCAFAI1WpV0UYswiTQ0gQAoBBqtZqijSgqlYpjCWhpAgBQCFqAiMjhBLQ0AQAoCiu4E4uOMqClCQAAMDQ6yoCWJgAARVGr1SwGShRagICWJgAARVGtVtMeAvmhowxoXQIAUCCKNgAQAABgyHSUAa1LAAAKpFKpKNqIQkcZ0LoEAKBwTN8kCh1lQIsSAIBi6ejosII7AEUmAADAcOgoA1qUAAAUi6KNuHSUAS1HAAAKx1mciEVHGdCKBACgcKrVqqKNWIRJoOUIAEDh1Go1RRtRVCoVxxLQcgQAoHC0ABGRwwloOQIAUERWcCcWHWVAyxEAAGD4dJQBLUcAAIqoVqtZDJQotAABLUcAAIqoWq2mPQTyQ0cZ0FoEAKCgFG0AFJMAAAAjoqMMaC0CAFBQlUpF0UYUOsqA1iIAAIVm+iZR6CgDWogAABRXR0eHFdwBKBoBAABGSkcZ0EIEAKC4FG3EpaMMaAkCAFBozuJELDrKgFYhAACFVq1WFW3EIkwCLUEAAAqtVqsp2oiiUqk4loCWIAAAhaYFiIgcTkBLEACAorOCO7HoKANaggAAAHHoKANaggAAFF2tVrMYKFFoAQJaggAAFF21Wk17COSHjjIg+wQAAEUbAAUiAABANDrKgOwTAABCpVJRtBGFjjIg+wQAgLVM3yQKHWVAxgkAACGE0NHRYQV3AIpAAACAmHSUARknAACEoGgjNh1lQGYJAABrOYsTsegoA7JMAABYq1qtKtqIRZgEMksAAFirVqsp2oiiUqk4loDMEgAA1tICREQOJyCzBACA9azgTiw6yoDMEgAAID4dZUBmCQAA69VqNYuBEoUWICCzBACA9arVatpDID90lAHZJAAAvI6iDYB8EwAAYFToKAOySQAAeJ1KpaJoIwodZUA2CQAA/TB9kyh0lAEZJAAA9NbR0WEFdwDySgAAgNGiowzIIAEAoDdFG3HpKAMyRQAA6IezOBGLjjIgawQAgH5Uq1VFG7EIk0CmCAAA/ajVaoo2oqhUKo4lIFMEAIB+aAEiIocTkCkCAED/rOBOLDrKgEwRAABgdOkoAzJFAADoX61WsxgoUWgBAjJFAADoX7VaTXsI5IeOMiA7BACADVK0AZA/AgAAjDodZUB2CAAAG1SpVBRtRKGjDMgOAQBgEKZvEoWOMiAjBACAgXR0dFjBHYA8EQAAYCzoKAMyQgAAGIiijbh0lAGpEwAABuEsTsSiowzIAgEAYBDValXRRizCJJA6AQBgELVaTdFGFJVKxbEEpE4AABiEFiAicjgBqRMAAAZnBXdi0VEGpE4AAICxo6MMSJ0AADC4Wq1mMVCi0AIEpE4AABhctVpNewjkh44yIF0CAEBTFG0A5IMAAABjSkcZkC4BAKAplUpF0UYUOsqAdAkAAENg+iZR6CgDUiQAADSro6PDCu4AtDoBAADGmo4yIEUCAECzFG1E5IQAQFoEAIAhULQRS7Va1VEGpEIAABgCRRsRCZNAKgQAgCGo1WqKNqKoVCqOJSAVAgDAEGgBIiKHE5AKAQBgaKzgTiw6yoBUCAAAkA4dZUAqBACAoanVahYDJQotQEAqBACAoalWq2kPgfzQUQaMPQEAYMgUbQC0LgEAAFKjowwYewIAwJBVKhVFG1HoKAPGngAAMEymbxKFjjJgjAkAAMPR0dFhBXcAWpEAAABp0lEGjDEBAGA4FG1E5IQAwFgSAACGSdFGLNVqVUcZMGYEAIBhUrQRkTAJjBkBAGCYarWaoo0oKpWKYwkYM+PTHgBAq6pWq7Nnz057FClbsKD9nPCNXldWpk9PZTAtrdFR5rQAwBhIjjjiiLTHANCqGsu3VyqVtAeSmv/5n7n1kPS6sjp16qz9909lPK3LsQSMGS1AAAxTZ2f/XSvlzs72BQvGeDCtTkcZMGYEAIDhq9VqRV4MtLNzg3Ogj1uwQAYYEotKAWNGAAAYvoJ3bG/oGwCGp1qtNhqBAEaVAAAwIh0dHQUv2s4N5/R7/XELFpQ7O8d4MC2t4HkSGDMCAAAj8o1wbhLq/d5UEgCGwjQAYGwIAAAjomhr6Pd7gONMAxgK0wCAsSEAAIyIoq3hG+HcfjOAqcBD4nACxoAAADBSiraGm/Zf2fdK0wCGpFqtlkobXFsJIAoBAGCkFG3dOvqcA1gAAMgaAQCA0aULqHmmlABjQAAAGClFW0NnZ6kyfXp16tRe1/sSoHnayYAxIAAAjJSirScBYIQcTsBoEwAA4lC0NVT6TAMIMsBQmFICjDYBACCCjo4ORVu3vlOBnREMIDsEAACGaerUas8fOzs3+B2IbwCaZ0oJMNoEAAAiq5kGMALmAACjTQAAiMCntqHHNwB95wEHGQAgMwQAgAh8atuLaQAj4XACRpUAABBN0Yq2qVNraQ8hnywEBIwqAQAgjmq1OvhGedfdBWQaAEBmCQAAcfjUtienAwPILAEAgGHqtQxoCKGzc30E6ncqMM0wpxwYVQIAQByKtl58CTBsJgEDo0oAAIhG0TYwCwEBZIEAABBHMScB9+0C6lbpsxIozfMlADB6BAAAouleBQiAzBIAAKIp4Ke2A3wD0Jc5AM2zqBQwegQAgJgKXrT1+gbAQkAAGSQAAERTzGkAA+gVAHwDAJAFAgAAwzd1ai3tIeSTVWWB0SMAABCTecAAGScAAERTwE9tB54EXDMHYLgKOKEcGDMCAAAxdXYONA3aNACA1AkAANH41LYXqwABZJAAAMCIDOlUADRPngRGiQAAQEwmAQNknAAAEFMBP7Ud+BsAXUAAWSMAAEAWVavVgp9YGhglAgAAMQ3cAlSyChBA2gQAAAAoEAEAIKYCtm1MnVpLewgADIEAAEBkFgICyDIBAIARcR6AUVKr1Yq2ohQwNgQAAAAoEAEAAAAKRAAAILLOzmJNgwZoLQIAACPV/DSAmhMDA6RNAAAAgAIRAAAAoEAEAABGaoAWoHJn51iOBIBBCQAAROZEYABZJgAAAECBCAAAjBb9PwAZJAAAMHaqlgEFSJsAAMBITZ1aS3sIADRLAAAAgAIRAACIrHsVoJI5AADZIwAAMFp6TQI2AQAgCwQAAAAoEAEAIKZyuVyrmRFLBKVSqVrd4CmWAYZNAABgpKZO7b9O1QIEkEECAACjou9ZwGoCAEAGCAAAAFAgAgAAZJH5JMAoEQAAYiqXyyZuNrQvWNDrGnMAALJAAAAAgAIRAAAYCz7+B8gIAQAgGv0/PfVdBYghcTgBo0QAACC+vtW/bwAAMkIAAGAsOAnAkPj4Hxg9AgBANKVSSdHW0HcJIAAyQgAAYKQ6O8uDbqMFCCAjBACAaJy5qVuvOQCq/6HybRIwegQAACJ7T7gt7SEAsEECAACRHdQnAFSmT09jIC3Mt0nA6BEAAKKxckvDe8Jv0x4CABskAADEofrv1rcFyByAoXI4AaNHAAAgJtU/QMYJAACMVGdnKe0h5IqP/4FRJQAAxGHdxoZzwjd6XWMGMECmCAAAcVi2pcEaoCMnTAKjSgAAIJp+q39zAIZKmARGlQAAEEeR+7Y7O8uNC33PAKD6B8gaAQAggiJX/z31PQOAADAMDidgVAkAAETTtwWoJgAMUblcTnsIQM4JAAARFHzWZqMFyASAWIp8LAFjQAAAiMCszdDfAqCq/2EoeJgExoAAAMBocQaAYRAmgdEmAABEUORZm939P84AANASBACAkSpy9d+t3/4fLUDD4HACRpsAABCBio0oVP/AGBAAAEaqvb29yE3bnZ2lfvt/TAAAyCYBAICR6tv/EywBNCzt7e2VSiXtUQA5JwAAjFTB2zY6O8t9P/5X/QNklgAAMCIFr/7DBs7/pf9neBxOwBgQAABGxGmb9P/EovoHxoYAADAiTttU7uzsdY3qf3iESWBsCAAAI+JT2770/wyPMAmMDQEAYPhU//3yDcDwOJyAsSEAABCT6n94VP/AmBEAAIbPqu196f8ByLjxaQ8AgNbW8fqK3zcAwyNMAmNGAAAYpnK5HELQtuEj/yjK5fKsWbPSHgVQCFqAACBlJgAAY0kAABgmPRsAtCIBAABSJkwCY0kAABiOcrmsbQOAViQAAAyT6p8ozCYHxpgAADAcejYAaFECAACkSZgExpgAADAcJgAA0KIEAIAhU/0TkcMJGGMCAMCQlUolFRtRqP6BsScAAAxZuVyu1Wppj4I8MAEAGHsCAMCQ+dQWgNYlAAAMjeqfiBxOwNgTAACGplQqpT0EckL1D6RCAAAYmnK5rGmbKMwmB1IhAAAMjU9ticVsciAVAgDAEKj+icjhBKRCAACAFKj+gbQIAABDYNV2YjGbHEiLAAAwBOVyOe0hkBNmkwNpEQAAmtXo2dC2QRRagIC0CAAAMNZU/0CKBACAZpkAAEAOCAAAMNaESSBFAgBAU8rlsrYNYjGbHEiRAADQLNU/UZhNDqRLAABoip4NAPJBAACAMSVMAukSAACaYgIAAPkgAAAMTvVPLGaTA6kTAAAG197ermIjFscSkC4BAKAptVot7SGQByYAAKkTAACa4lNbAPJBAAAYhHM2EZEJAEDqBACAQZRKJRUbUaj+gSwQAAAGUS6XTQAgCrPJgSwQAAAG4VNbIhImgdQJAAADUf0TkcMJyAIBAGAgpVIp7SGQE6p/ICMEAICBlMtlq7YThdnkQEYIAAAD8aktsZhNDmSEAACwQap/InI4ARkhAADAqFP9A9khAABsUHt7uwkAAOSMAACwQeVyOe0hkBPCJJAdAgBA/xo9G9o2iEILEJAdAgAAjC7VP5ApAgBA//RsAJBLAgAAjC5hEsgUAQCgH+VyWdsGALkkAAD0T/VPFI21pBxOQHYIAAD90LMBQF4JAAAwioRJIGsEAIB+mAAAQF4JAAC9qf6JxWxyIIMEAIDe2tvbVWzE4lgCskYAAOhHrVZLewjkgQkAQAYJAAC96dkAIMcEAIDXUf0TkcMJyCABAOB1SqWSio0oVP9ANgkAAK9TLpdNACAKYRLIJgEA4HV8aksswiSQTQIAwHqqfyJyOAHZJAAAQHyqfyCzBACA9azaTiylUintIQD0TwAAWK9cLqc9BHKiXC4Lk0A2CQAAazV6NrRtEIUWICCzBAAAiEz1D2SZAACwlgkAABSBAAAAkQmTQJYJAAAhrJv+q22DKMwmB7JMAACAmMwmBzJOAAAIQc8GAIUhAABATMIkkHECAEAI1m0EoDAEAADVP9GUy2WHE5BxAgCAng1iUv0DGScAAEA0wiSQfQIAgBYgAApEAACKTvVPRA4nIPsEAKDoSqWSio0oVP9ASxAAgKIrl8u1Wi3tUZAH7e3tAgCQfQIAUHQ+tSUiYRLIPgEAKDTVPxE5nICWIAAAhVYqldIeAjmh+gdahQAAFFq5XLZqO1GYTQ60CgEAKDSf2hKL2eRAqxAAgOJS/RORwwloFQIAAIyU6h9oIQIAUFzt7e0mABCF2eRACxEAgOIql8tpD4GcMJscaCECAFBQjZ4NbRtEoQUIaCECAACMiOofaC0CAFBQJgAAUEwCAACMiDAJtBYBACiicrmsbYNYzCYHWosAABSU6p8ozCYHWo4AABSRng0ACksAAIDhEyaBliMAAEVkAgAAhSUAAIWj+icWs8mBViQAAIXT3t6uYiMWxxLQcgQAoIhqtVraQyAPTAAAWpEAABSOng0AikwAAIpF9U9EDiegFQkAQLGUSiUVG1Go/oEWJQAAxVIul00AIAqzyYEWJQAAxeJTWyISJoFWJAAABaL6JyKHE9CiBAAAGDLVP9C6BACgQKzaTixmkwOtSwAACsSntsRiNjnQugQAoChU/0TkcAJalwAAAEOj+gdamgAAFIUJAAAQBAAAGCphEmhpAgBQCOVyOYSgbYMoGocTQIsSAABgCBoTAIRJoHUJAEAh6NkAgAYBAACGQJgEWp0AABSCdRsBoEEAAPJP9U8sZpMDOSAAAPmnZwMAugkAANAsYRLIAQEAyD8tQADQTQAAck71T0QOJyAHBAAg50qlkoqNKFT/QD4IAEDOlcvlWq2W9ijIg/b2dgEAyAEBAMg5n9oSkTAJ5IAAAOSZ6p+IHE5APggAQJ6VSqW0h0BOqP6B3BAAgDwrl8tWbScKs8mB3BAAgDzzqS2xmE0O5IYAAOSW6p+IHE5AbggAADAI1T+QJwIAkFvt7e0mABCF2eRAnggAQG6Vy+W0h0BOmE0O5IkAAORTo2dD2wZRaAEC8kQAAICBqP6BnBEAgHwyAQAA+iUAAMBAhEkgZwQAIIfK5bK2DWIxmxzIGQEAyCfVP1GYTQ7kjwAA5JCeDQDYEAEAADZImATyRwAAcsgEAADYEAEAyBvVP7GYTQ7kkgAA5E17e7uKjVgcS0D+CABADtVqtbSHQB6YAADkkgAA5I2eDQAYgAAA5Irqn4gcTkAuCQBArpRKJRUbUaj+gbwSAIBcKZfLJgAQhdnkQF4JAECu+NSWiIRJIJcEACA/VP9E5HAC8koAAPKjVCqlPQRyQvUP5JgAAORHuVy2ajtRmE0O5JgAAOSHT22JxWxyIMcEACAnVP9E5HACckwAAIDXUf0D+SYAADnR3t5uAgAADEoAAIDXESaBfBMAgDwol8shBG0bRKEFCMg3AQAA1lP9A7knAAB5oGcDAJokAADAesIkkHsCANDyyuWytg0AaJIAAOSB6p8ozCYHikAAAFqeng0AaJ4AAABrCZNAEQgAQMszAQAAmicAAK1N9U9EDiegCAQAoLWVSiUVG1Go/oGCEACA1lYul2u1WtqjIA9MAAAKQgAAWptPbQFgSAQAoIWp/onI4QQUhAAAtDATAIhF9Q8UhwAAtDATAIhFmASKQwAAWphPbYlFmASKQwAAWpXqn4gcTkBxCAAAFJ3qHygUAQBoVVZtJ5ZSqZT2EADGjgAAtKpyuZz2EMiJcrksTALFIQAALanRs6Ftgyi0AAGFIgAAUGiqf6BoBACgJZkAAADDIwAAUGjCJFA0AgDQehrTf7VtEIXZ5EDRCAAAFJfZ5EABCQBA69GzAQDDJgAAUFzCJFBAAgDQeqzbCADDJgAALUb1TyzlctnhBBSQAAC0mPb2dhUbsTiWgAISAIDWU6vV0h4CeWACAFBMAgDQYvRsAMBICABAK1H9E5HDCSgmAQBoJaVSScVGFKp/oLAEAKCVlMtlEwCIwmxyoLAEAKCV+NSWiIRJoJgEAKBlqP6JyOEEFJYAALSMUqmU9hDICdU/UGQCANAyyuWyVduJwmxyoMgEAKBl+NSWWMwmB4pMAABag+qfiBxOQJEJAAAUi+ofKDgBAGgN7e3tJgAQhdnkQMEJAEBrKJfLaQ+BnDCbHCg4AQBoAY2eDW0bRKEFCCg4AQCAAlH9AwgAQAswAQAAYhEAACgQYRJAAACyrlwua9sAgFgEAKAFqP6JorGWlMMJKDgBAMg6PRsAEJEAAEBRCJMAQQAAss8EAACISAAAMk31TyxmkwM0CABAprW3t6vYiMWxBBAEACD7arVa2kMgD0wAAGgQAIBM07MBAHEJAEB2qf6JyOEE0CAAANlVKpVUbESh+gfoJgAA2VUul00AIAphEqCbAABkl09tiUWYBOgmAAAZpfonIocTQDcBAICcU/0D9CQAABll1XZiKZVKaQ8BIEMEACCjfGpLLOVyWZgE6CYAAFmk+icihxNATwIAAHmm+gfoRQAAssgEAAAYJePTHgAAjKKxCZPHHHPMCSec0Li8evXqlStXdnZ2Lly48Oabb37ggQfq9XrPjWfPnr1mzZqvfe1rA+/z8MMPP/XUU4877rhly5aN1riBQhIAgMwpl8shBG0bRDGWM4AvvPDCEML48eMnTpy49dZbv/Od7zz44IMXL1787W9/+7HHHhubMQAMSgAAILcaEwDGLEzecccdPX+88sorZ8yYMXPmzIsuuuhf//VfH3nkkcb1v/3tb7u6usZmSAB9mQMAZI4JAOTGI488csYZZ3R2dp5xxhkTJ05sXDl//vx77rkn3YEBReYbAAByKwth8tVXX/3+979/4YUXHnjggTfddFMI4ayzzuqeAzBlypRPf/rTb3vb21577bWbb765UqmsWbPmggsu2GOPPRp37+joCCH88pe/vPLKKxvXTJ069Zhjjtl111133HHH5cuXP/jgg1ddddWzzz7b0dHxgx/84De/+U1YN83gggsuOPHEE/fcc88pU6YsXLjwzjvv/N///d9ew9tss82OPvroGTNmTJs27eWXX3700Ufnzp1733339dpsyy23POaYY6ZPn/7mN7/5+eeff+CBB379618vWrRoFJ84YNQIAEDmlMvlWbNmpT0KiOahhx5atGjRP/zDPzQCQLckSb761a9uu+22119//eTJk4855pg1a9ZUKpWOjo7NNtvsbW9726GHHnrxxRe/+uqrTzzxROMuBx544Kmnnrpq1ao777zz1ltvnTp16r777nvRRRf99Kc/7fWgW2211b//+7+/+OKL8+bNe+211/bcc8/Pf/7ze+2110UXXdTdgLTffvvNnDmzq6vrrrvuuv322zfZZJN3vvOd3/zmN2+66aZLLrlkzZo1jc0OOuigU089ddmyZXfcccdNN9201VZb/f3f//3BBx/8wx/+sNdvBLQEAQDIFqu2E0umZpM/8sgj++yzT68rt9tuuxkzZlx66aXXX399COENb3hD44P/+++/P4QwadKkQw899M477+xeBWjHHXc8/fTT77333u9973vdV1555ZVHH330Jz/5yV4733777X/961//53/+Z2MNol/96leHHXbYaaeddvfdd996662NDc4444xqtfrd7363e2+/+MUvTjzxxKOOOmrZsmU/+clPQgjTpk2bOXPmbbfddskll6xevbqx2c9+9rNPfepT//zP//zMM8/0/boAyDgBAMiWLPRsQHQvvPDC5MmT+71pypQpjQsXXXTRwDs5+eSTly5dOmfOnOXLl3dfWa/Xf/GLX+y8884HHHBAz42XLVt25ZVX9lyB9KabbjriiCMOO+ywRgD4zGc+s2zZsjlz5rz66qs973jllVcuWrRo8uTJ48aNW7Nmzcknn7x48eKLL76458Tl1atXX3bZZbvsssvHP/5xAQBajgAAQD5lKkxuttlmL730Uq8rFy9efNdddx177LE77LDDddddV6vVep0xoKckSXbbbbe5c+f2rP67/eY3v+kVAJ566qlVq1b12uyxxx7be++9G3vbY489rrvuul7Vf0MjIYQQ2traZsyYMX/+/He96119N3v88cff9773bbrpps5UAK1FABiRxllaGpeffvrpT3/606PxKEceeeQJJ5xw7LHHrly5ctq0ad/73vfOP/98K0iQVyYAkEu77bbbU0891ff6Cy644H3ve9/73//+2bNnP/LIIxdddNHTTz/d7x622mqriRMnLlmypN9bFy9e3OuafnPCiy++OGHChBDClltuOXHixCeffHLgYW+55ZYTJkw48MADDzzwwA1tM2XKFAEAWosAMCIPPvjgv/3bv22xxRaf+cxn7rzzzlF6lPnz53/yk58slUr33nvvqaeeOm/ePNU/eWUCABFlJ0zus88+O+200w9/+MO+N61Zs+b666+//vrrS6XSF7/4xfPOO++zn/1svzt57rnnXnvttR122KHfW9/85jf3uqbfLxO6r3zuuedWrFix/fbbDzzyZ599dsWKFR0dHf/93/898JZAC3EegBFZtGjRnXfeuc8++zz11FN9V2CIZcmSJU888cQ+++xz2GGHbbvttpdddtkoPRCkrlQqCQBEkZ0wOWXKlC996UudnZ233XZbr5tmzJix//77Ny7XarWf/exn22233Zve9KZ+99PV1bVgwYIDDzxw0qRJfW897LDDhjSqer3+5z//+cADD9xkk0363loulw866KBx48bV6/UHH3ywcbnvZqeccsonPvGJIT0ukAUCwEgdcsghf/d3f/fd73535cqVo/cojf7LE0888Uc/+lHfLlLIjXK5XKvV0h4FedDe3p6FAFAqlc4555xJkybNmTPnlVde6XXr3nvvPWvWrO4P9TfffPMQQt/Nul1++eWTJ0/+0pe+tOmmm3ZfmSTJUUcd1R0kmnfZZZdtttlmZ5xxRs+9hRAOOeSQ888/v1QqNZYBvfzyy3fYYYfTTjtt44037vmgJ5100uGHH95YsAhoLVqARmTq1Kknn3zy1Vdf/fDDD4/qA/3+978/6qij5s+f3+s885Az2enZIAfGPky++93vDiGMHz9+woQJ22yzzV577bX77rs/++yzZ5999oIFC/puf8stt3z0ox8966yzbrzxxs033/wjH/nI7bffPkA//aOPPvr973//85///KWXXjpv3rxFixZNmTLl7W9/+1vf+tbLL7/8M5/5zJBGu2TJkgsvvHDmzJk//OEP582b9/jjj2+11ValUmm33XabN2/ej370o8ZmTz755He+850vfOELe+2117x585YsWTJ58uQDDzywcZ4BAQBakQAwIqeddtrSpUt/9rOfjfYDPf/88yGEX//616P9QJCi7PRskAOphMmvfOUrIYR6vf7SSy89//zzS5YsmT179j333NO9fH4vzz333KxZs0466aTjjz9++fLlN95446DdpLfddlu1Wj366KN33XXXQw45ZNWqVQ8//PBZZ531xBNP7L///o3/s2jePffcc8oppxx11FG77777QQcdtHz58scff/zrX//6vffe23OzO++88+GHHz766KNLpdL73//+Z5999uGHHz7nnHOeffbZIT0ckBHJEUcckfYYWtWhhx76hS984aqrrvrLX/4SQpg/f/7oPdY73vGOc84557jjjrPSAjnW3t7uGwCiKJfL7e3tjiWAfvkGYPje8573hBA+9rGPhRCee+65UQ0AO+2009KlS1X/5Fu5XM7Oqu20NLPJAQbgGwAgK+bOnfvBD34w7VGQB7Nnz65UKjIAQL8EgCFYsGDfzs7tOju3nTr1qenT/zh1aj+ndIllWdeKtgnjJ63yFQ2ZM2nSpJ///OeNy3/84x/PPffcgbf/7Gc/2+vvzJNPPnnKKad0/7jtttt+/OMfL5VKkydPfuihh2677babbrppgPOhNmnHHXc89thj3/rWt2622WYLFy78wx/+cO211/bc7cYbb9ze3v7ud7978uTJjz76aKVSyfF0xkGfjV4GfdUG3eFQHzGuVgmTQ303Zfxpb3XR/7h5OUZi5513PvbYY/fYY4+NNtrosccemzt37u9+97sBtvfuGBL1ZbMWLNh3wYK3Ny53dm57111HTJ/+x+nT/zR6j/jioVu/8tTyzR9dOeGFNaP3KDBUK1asOOecc8K6/rdBTZs27e67777uuut67qH78j777POv//qvjz766Lx58/bdd99HH330c5/73Lve9a7zzjtvJH+X3/Wud5155pl//etfb7jhhs7Ozl122eWEE07Yc88958yZ01ixN0mSr371q9OmTfvVr361dOnS/fff/5vf/ObXv/71++67b9gPmlmDPht9DfyqDbrDYTxiRC00m3xI76aMP+05EPePm5djJEql0rnnnrto0aKrr766s7Nzr732+vKXv3zNNddcccUV/W7v3TFUAkBTOju37a7+uy1Y8PbRCwAvdS1v22jcxvvvsGzP5U8v6pzwwprJi1ZLAmTBmjVr/vSnP4UQDjvssIkTJw66/U477fTzn/+8cZdeJk6cOHPmzPnz58+ZM+f888+/9NJLq9Xqb3/72wsuuOAjH/nIsM88uvHGG59yyil33333nDlzGguZ33777bfccst3vvOdI4444uqrrw4hlEqlfffd91/+5V8ak/jnzZu3ySabHHfccfkLAM08G30N8KoNusPhPWJEpVJptB8ilubfTdl/2nMg4h83L8dIJEny+c9/fuHChV//+tcb1fntt9/+0EMPzZw584477mj80e7Ju2MYBIARefDNH3zr554ejT1P/L9n1ixftdFmb9hoszdM3GbzEMKrLy1f/ezKiX96ftyil0fjEWE0bLHFFptuuumiRYv6vXWvvfbacsstr7jiinq93j0D+JFHHrnlllsOPvjgngFg1113PfLII3fZZZcpU6YsXrz41ltvve666za0tOJuu+22xRZbdHR0NP7QNzz22GO/+93vDjrooMbf+p122mn16tULFy7s3uDPf/7zRz/60Si/daY082z0MvCrNugOh/GIceVyNnn2n/aiGeHbZKyG2ZLe/OY377DDDl/72td6fjZ/8803f+pTn9p77737BgDvjmFwJuCmdHZu2+/1W7xjg+drjKtr1Zpxb9h45Yuvrt528M8kIDumTZsWQli0aNF22203Y8aMXh+q7bTTTsuWLevs7Gz0bHS3bSxcuHDHHXdMkqTx45FHHnnRRRdNmjTpmmuu+Y//+I9arXbCCSdceOGFvU5f2m2rrbYKISxZsqTX9UuWLJk6dWrj8lNPPTV+/Pg999yz+9ZSqfTYY4+N+DfOnGaejV4GftUG3eEwHjGuFmoBal72n/aiGeHbhAG0tbXdc889TzzxRK/rly9fPmXKlL7be3cMg28AmjJ9+p/6tgCFEDb6z4cmXDMqU4GfW/182/+bvrLz5TXLV65ZvnLSg68kC1/YzJxgWs20adPq9fp555230047hRDq9fr8+fMvueSSxumKnn/++Te+8Y0TJkzoda+pU6d2dnY25gCUSqV/+qd/uuyyy/7nf/6ncettt902d+7ciy66aObMmd/61rf6PuikSZNWrFjR85OehldeeeUNb3hD4/Kf/vSnBx54YNasWVdfffUzzzyz//7777HHHrNnz47622dCM89GLwO/aoPucBiPGFEuq/+Q+ae9gEb4NmEAjz/++Hnnndfryre85S1bb711z69tu3l3DINvAJo1ffof+14zqgsBvfb0i22/eWKjeU9v86un3/jQy1YEohVNmzYtSZL77rvvtNNOO/bYY7/xjW9sv/323/nOdzbZZJMQwp///Oeurq4PfOAD7e3t3T0bb3zjG9/73vc+8MADjR8/+clP3n333d3Vf8PTTz998cUX77fffo3/621SV1dX9+V6vX7ttddusskmJ5544plnnnnAAQc8+OCD//d//zfSX7h19Hw2ehn4VRvGDpvcgGHwtKdllN4m9GvixImnn376o48+escddzR/L++OAagpmzV9+p+mTn1qwYK3d3Zu+55dly5JHh3VJYC2Hz8l/L6I09LJmauuuurWW2+99957Gz/+8Y9/XLhw4SWXXHL88cdfdtllTz311K9+9atPfOITS5cu3WKLLSZNmjR9+vRPfepTkyZNuuqqq0IIG2200Vve8pb7779/77337rXnVatWrVixYt9993388ceHMbCDDjroS1/60tVXX33DDTe8+OKLM2bMOPnkk7/97W+feeaZr7766gh/61Y38KuW7tgG1TNMwuhp6bdJa9lkk02+/vWvb7755mefffaGpn4xVALAEEyd+tQPTrjjuF0271j4YmXhi2kPB1rA4sWLFy9e3POaF1988dZbb33729f21P30pz994xvf+IEPfOCMM84IIaxevXr8+PE//vGP//a3v4UQtt5663Hjxh111FFHHXVUv/ufMmXKP/7jP37uc5/rvua+++675557uucP9NTWtvY7z/Hjx3/605++4YYbfvKTn3Tf69xzz7300ksPP/zw/E0IG/jZ6GvQV23QHQ71ESPK5Qzghiw/7QU08rcJzdhss83OO++8zTfffNasWU89tcG2C++OoRIAgNHylre8Zeedd7755pt7Xf/ss882pmSFEOr1+rx583bcccfLL7981apVZ5xxxqpVq6699trGrc8888yaNWsuvPDCO++8c0OPsvnmmz/yyCPdPy5fvny33XabMGHCRhtttGrVqp5bTpo06eWXXw4hbLPNNpMnT77rrrt6jWrhwoW77rrrCH7jLHrllVcGfjZ6GfRVG3SHQ33EiHrNJs+TLD/tBTTytwnNmDp16je/+c1x48Z9+ctffuaZZza0mXfHMBQ3+gCj7U1vetPpp5++zTbb9Lr+rW99a8/1dtrb2zs6Ov7617/ut99+O+644w9+8IPuU4CtWrXqr3/96wEHHNB35zNmzDj77LO33HLLF198cWEPixcvbvz/xA477NDrLttvv/1zzz0XQnjppZdCCH3nfm2yySYvvpi3L/cGfTZ6GfRVG3SHQ31EmuFpz5SRv00Y1Jve9KYLLrhg9erVA1f/wbtjWAQAIKaNN9740EMP3XrrrUMIf/jDH/72t7999rOf7fk16+67737QQQfdcsstve643Xbbtbe3/+IXv+i1rvaPf/zjAw444AMf+EDPK7fccssvfvGLEydO7Pdv94MPPrh06dJjjz123Lhx3VfuvPPO++23X+ObhGXLlv3lL3/54Ac/2PNL4b322mvatGn333//8H/5TBr02QhDfNUG3WEzjzhKcjYBoOfrkuWnvSDivk0Y2A477HDBBRcsXbr0K1/5ygsvvNB3A++OEUqOOOKItMfQStp32dwcAGg4++yzJ06ceO655/a8cvLkyf/1X//1rW99a/78+SGE6dOnn3vuuU8//fSNN9740ksv7b777kccccQf/vCH2bNnd3/MP3fu3A996EPnn3/+5ptvPnPmzL5zvD784Q+fdNJJ1Wr197///QsvvLDLLrscfvjhTz/99Ne+9rUNfXs7Y8aMc84558knn7zjjjs6Ozt33XXXww8//Mknnzz33HNfeeWVEMKOO+44Z86cJ5988sYbb2xMAv7Qhz50zz33fPvb347+RKVu0GdjqK/aoDscdINRMnv27Eql0ootQM28mzL7tOdPlD9uXo5ha1T/K1euvPTSS1esWNHzpmeeeaYxE8C7Y4QEgKERAKBbM/8fGULYeuutP/axj+27776TJk167LHHbrnllhtvvLG7+i+Xy+3t7b/97W9PPfXUM888s2c3f0/Tpk076qijdtllly222OLxxx//3e9+d+211w68gtsOO+xw/PHH77nnnhMmTHj88cer1WpHR0fPBtAtt9zy+OOP32OPPaZMmbJo0aJbbrnl+uuv7x5Yzgz8bAz1VRt0h81sEF25XJ49e/YHP/jB0XuI0dPkuymDT3suRfnjFrwcw3XIIYecfvrp/d70y1/+8sorrwzeHSMmAAyNAABxzZ49u1qt5qltg7Q0wuSsWbPSHghA1pkDAKSsVqulPQTyIGcTAABGjwAApKmxbmPaowCAAhEAgNSo/onI4QTQJAEASE2pVFKxEYXqH6B5zgQMpKZcLg/atL3iuokhhNBYsD9JQkiSJAnJ2svrLoQkSULjf+s2C23rrg+hcTkkbT0u99iyscPujbtPDpAkobHId8+dNPS4vr7u7muv6d5J0mNX6+5SX7ur0H339TtZv+W6HSZJPYSQhO5drd0s6XH3EEII9bb1u+25TX39ozRu6h7nurus3W3osdveO1n3ECEkG9jh6y733mE/l3vsuZ40d5cN76Tx4xMvhut/+Jv6Tv9fr43XD36wXdV73XGAwW/omkF2Uh9oDD2vef2jvP4XqQ90xx4P1OOY7P4F633Htm6beu/r+16z7u5Jnw2S113usdfG8bv21vq6K8P64yvpvVk/16zbVdLn+h7vv3rPO7b1uGNbz2tef7nxivXY/7p3Z4/NGsYna/ruqvvu6y8n9cYbq3u3bUk9hHpbUk+67x4aO+lav8PQlfT4BRu7agtrx5OEevcde14OIYxbu5OudbeGJOlKui93j2fdAZWE+vika91O6ut2Ug8hjAu9d9WWdCUhJI2xrf0V1v4NSNaNJFk7krWjauxkXM9fMNTXj2f9r7bueVv7XK17TkIIIYzr8Qdy3NrnMEnW/RVsaAvJuCRJQmgL6//utyVJCElbSBobhBDGJW0hhA+fvFugP74BAFLjU1sAGHsCAJCO4VT/9cY/9e7Lg2w60O09d7KBjQc9J8DaDZo9dcBA273usfrfsOmH6f/5qffZZGg7GXiHA+2qz78b2GSgXW1gJ82OYYjjaXaHTTxkvdd/RzSeQQ/IZnc1uEF21dQTv/YgGmQ8zQ23x642fPQ2tav+xtP312nuqe57eQMbDe2128DNzR5LTTzOCLdbN9ihHGQbfO83+We2mX2P8KAvFgEAAAAKRAAA0mHVdgBIhQAApMMEAABIhQAApED1DwBpEQBG5OCDD547d+7ee+8dQpg4ceLFF198zjnnJD1XQQMAgCwRAEbk9ttvf/bZZz/84Q+HEE4//fTx48fPmTOn3sSMdSg4EwAAIC0CwIisWbPmmmuuecc73vGFL3zhbW9727e+9a1XXnkl7UEBAMAGCQAjddNNN7388svvfe9758yZ8+STT6Y9HGgB5XI5hGAOAOSbfljILAFgpLbZZpuNNtqoq6vrscceS3ssAAAwCAFgRDbddNOvfvWr999//6pVq4488si0hwOtwQQAAEiRADB8bW1tZ5111muvvTZnzpwbbrjhsMMO22yzzdIeFAAADEQAGL6TTjppl112mT179ooVK6655ppx48Z96EMfSntQkHXlctlJAMiJrPW4Z208QFaNT3sALeyKK6644oorGpc7Ozu1AEGTVP8AkCLfAABjygQAAEiXAAAAAAUiAABjygQAAEiXOQDN+uiapcmq5du98Nptf3hm5d+W/r/VXUsmTL4rmZT2uKCVqP4BIHUCQLMmbD75uG23CdtMDuPHhd3eHF54edaCV8PKtIcFLaVUKv3s4qu6lo4LIYSQhKSxbkmSJGHtCibJuuu7/w31tEYLALkkADSr8nLbcSGENV1h/LiwpissW15dacU1GJoHfvvgvVc/FJKNkrYkhCSMawtJEpIkaWsLbUkIbaEtCUlbaAtJkoS2tpAkISxNe9QAkCsCwBB0PPXacVuGEEJYvaa60geTMGR/XnL/xu8a2l1WPTxx/Q9rvyTo/q5gAIPl86TnTjawcdLcTppefX2gvxqve6z+d9jsH52k/+en592TZLBd9d3JwDsc6Ffr829/mwyyqw3spNkxDHE8Q9vhYI9WH3hnQxhPUm9mREM7KpvYST+7GsITP9i7pHH3wZ7sZP2uNvy+bGpX/Y2n76/T1Hj6udzPeBqXmtvhgLtKelwefFQDPs4It1s32KG8RTb43h/0mRl8yEmv/9IEAWAInluxMrz6WpiwUVi5utLZ5VCDMTDxH1ekPYRCqvf4t6ErnYFA6tb0+DcbGuXHuB7/wtBYBWgIfrN64nNLl4UQqq926f8BAKAVCQBDc9Or40IIty5+Ke2BAADAcAgAQ1Pp7Aqvraq+nKVvAgEAoGkCwJDNemjZ39ompD0KAAAYDgFgyKovWPwfAIBWJQAAAECBCAAAAFAgAgAAABSIAAAAAAUiAAAAQIEIAAAAUCACAAAAFIgAAAAABSIAAABAgQgAAABQIAIAAAAUiAAAAAAFIgAAAECBCAAAAFAgAgAAABSIAAAAAAUiAAAAQIEIAAAAUCACAAAAFIgAAAAABSIAAABAgQgAAABQIAIAAAAUiAAAAAAFIgAAAECBCAAAAFAgAgAAABSIAAAAAAUiAAAAQIEIAAAAUCACAAAAFIgAAAAABSIAAABAgfz/gEWtxQu9piMAAAAASUVORK5CYII=\n"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}