{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Definition and plotting of a CHS section" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[,\n", " ]" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAOUAAADFCAYAAACrZTUgAAAAP3RFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMS5wb3N0MSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8kixA/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABI8ElEQVR4nO19eXRb1bnvJ8my5VGeLTuOh3jMSCCQrOTShldySUo6916GsgrNo0lL6fS4pZAuCAXubZi6Lq+lXbS9C/L6CqTQewu9LaQFSrm3EFKGl8nzGE+Sbdmx5dmWtN8fP7b3OUfnHB3Jli0n57eWli3pzNq//X37Gy2MMUYmTJiIG1iX+wJMmDAhh0lKEybiDCYpTZiIM5ikNGEizmCS0oSJOINJShMm4gwmKU2YiDMkLPcFLAeCwSD19fVReno6WSyW5b4cExcoGGM0NjZGRUVFZLUal38XJSn7+vpo9erVy30ZJi4SdHd3U3FxseHtL0pSpqenExEeVkZGxjJfjYkLFT6fj1avXj0/3ozioiQlV1kzMjJMUpqIOSJdIpmGHhMm4gxLQsqf/OQnVFZWRg6Hg7Zt20Z/+9vfNLc9cuQIWSwW2cvhcMi2YYzRoUOHqLCwkJKTk2nXrl3U0tIS69swYWJJEHNS/vrXv6Y77riD7rvvPvrggw/okksuod27d9PAwIDmPhkZGeR2u+df586dk33/yCOP0I9+9CN68skn6cSJE5Samkq7d++m6enpWN+OCROxB4sxtm7dym6//fb594FAgBUVFbHDhw+rbv/0008zp9OpebxgMMhcLhd79NFH5z8bGRlhSUlJ7LnnnjN0TaOjo4yI2OjoqLGbMGEiCkQ7zmJq6JmdnaX333+fDh48OP+Z1WqlXbt20fHjxzX3Gx8fp9LSUgoGg3TZZZfRD37wA1q/fj0REXV0dJDH46Fdu3bNb+90Omnbtm10/PhxuuGGG0KONzMzQzMzM/PvfT7fYtzeBQXGiHw+Iq8X/3NwG8XICFFdHdGOHaH7WixEeXlE6eliexPRI6ak9Hq9FAgEqKCgQPZ5QUEBNTY2qu5TU1NDTz31FG3atIlGR0fpscceox07dlBdXR0VFxeTx+OZP4bymPw7JQ4fPkz333//ItzRykYwSNTZSXTmDNG6daEEcjqJysuJtPzcnKyVlfLPAwGQWW1FMjBAVFhIVFZmEtYo4s4lsn37dtq+ffv8+x07dtDatWvpZz/7GT344INRHfPgwYN0xx13zL/n/qOLAVNTRN3dIITFAnIMDRElJRGVlER2rMxMSMSWFqKqKvG5zUZUUICXGlatImptxf+MEZWW4vwm1BFTUubm5pLNZqP+/n7Z5/39/eRyuQwdw26306WXXkqtH/6qfL/+/n4qLCyUHXPz5s2qx0hKSqKki2gUeL1E58/jf4eDqLpafMcYUXY2SNHXR1RUFNmxnU6QW0lMPSQmim0ZI+rqIpqZwXFycnA9JgRian1NTEykLVu20Ouvvz7/WTAYpNdff10mDfUQCATozJkz8wQsLy8nl8slO6bP56MTJ04YPuaFBsaIzp0DUVpaoH5WVeGlVAjOnwcJCgqw3+Cg8fNYLFCBMzKIXC6i5ubIr9VigaSsrsb1BYM4TksLJLpZMYpib309evQoS0pKYkeOHGH19fXswIEDLDMzk3k8HsYYY1/84hfZ3XffPb/9/fffz/74xz+ytrY29v7777MbbriBORwOVldXN7/NQw89xDIzM9lLL73ETp8+zT796U+z8vJyNjU1ZeiaLhTrq8fD2EsvMdbczNj0tLF9WlsZCwbF+44OxoaHje07O4vtOXw+xhob9fdpaTF2bMYYm5zEvTz7LGODg8b3i1fEpfWViOj666+nwcFBOnToEHk8Htq8eTMdO3Zs3lDT1dUli6A/f/487d+/nzweD2VlZdGWLVvo7bffpnXr1s1v893vfpcmJibowIEDNDIyQldeeSUdO3YsJMjgQsXwMAwomZlQT42qkRxSg0tZGaSU3U6Ulqa/n91O5PeL9+npWC82NRHV1ER2DWqYmCCamyO65BK8b2yEkcjpXPixVxIsjF18CoPP5yOn00mjo6MrKvZ1fJyop4coKwvqZ0cHVFGfL1RN1UJbG1FFRejnDQ2wvIab11pbQ62vExNQPWtrjW2vhN8PFTYnB/clXa/29eH+SkqIUlL0jxNviHacmbGvKwAzM5AaIyMY+NzKGQhAikxMLPwca9eCGIFA5PumpmJSaGiIfN/ubkwUa9eqW2+LinDPXi+ewdxc5OdYaTBJGccIBDAQe3sxMLVS8oqKsI0R6OlFGzfCh6m3jcMBN4sSqakw4Bgl5vg4UX09JpWaGqFS+/1ECSqLqpISbNfVBXU5GDR2npUIk5RxCMagwrW1YSCuWaO+HR/IGRlEY2MLP6/FAmKePau9zapV2hNASgrWqPX12vvzexsaQgCDUqvr7tZWxS0WqN5VVeL5XIiIu+CBix2dnUTT01iHqUkMKaQSzeUi8njwdyGw2eCuqK8HaZQIF5WTnIy1aV0d0YeRkfPgUT+VlfBdqiEQCH/fVismq9lZaBIpKZEHQsQzTEkZJwgEiE6fhqSoqgo/MJXIzMSaczHAo32izYZLToZE4xLX7wfJAwEQXYuQkcJmwyTxy18SffDBhaPSmqSMA4yNYS22cSPRlVdCyszO6u8TDIbGqObnq8efcvh8oeqiFtLSEFLX2WlseyW4q+bXv8YaUMuQEw3m5nDMjg5YbP/n/8TatK6OaHJycc6xnDBJuczo70dUzYYNIj510yZYQvXSQ/v64MOTIjsbPkwtDA4S5eYav7bMTBhwenrkn+fmGosE6u4m+tjHIHmbmsJPNOEwM4PjdHVBxa6sxCRUVARVfuNGrHe93oWdZ7lhknIZ0dGBwaRmyNmwgai9XdvdMT0NNVGJnBztQcmYdgaIFvLyoCZKw5czM4lGR7X3YQyqq8uF/SsrsQbs7sYaUJJFJ8P0tHqg+tQUyNjXh+NUVGDy6uyExVeKqiocv6srsvuMJ5ikXAYwBlUrJ0ffMLNuHaRUJOmfeXmLLykKCzHQh4bCbzs3h7VxbS1UYG6M4pbT2lpIs8bGUE2gtxfWXY6JCZBxYABkLC8X3wUC2D81NfQaVq1CtFE0ftN4gGl9XWLwQID1640Zc2pqEBUTDEJCGUFWFtTYxcy+KCmBC8Ju116Xjo1BGm7apG+l5ZoBtzTzaB0uyX0+IrcbhNMK32tuVo8g4sjKwjFPnoTWEanhbDmxgi515eP8eaiBPLbTKCorocoGApCu4VBQAAmz2ClRFRWYUBISQB6psWlwEGSSulG0AgE4ysrw99w5qKgDAzhmWpp+LK3PF1rlwGYLPV9SEp51fT18nyslotJUX5cIvb2QJHqzux7WrEEUjNFUq4yMxXORSFFbiwnC5YJUJMLf2dnQmFq3O9QYpYbcXEw4J07AgBMux7OnJzS6qbg41CBFBOKuXw/VW6MwRdzBJOUS4NgxDNyZGZHzeO4cJMPYmDzzQg+lpVD3jAyuwsLYDcL163EPMzNQI1NS5GtBjpkZ/QD3uTms+4aGoBZ/6lPhz61mdSYKzWBRorwcBD15Mvw5lhum+hpjdHVBjbz8cqHqMYYBOTUFY4bXKwaU2lpMaiwhIvq//1ckMOsVVEhPj8w3aRRzc7D8PvQQ0eHDxlRqKXiondUKyWuxQC2229X9r9L9RkYir5bAkZEBjUVpUIo3mKSMIbxerHUuvxzGGl6Ww2JBVEtiYmS5grOzGMwf/zjI6HZDGuXkqPsfV63CYM/IAPmjTX3y+yHZ+eRgt+Pav/ENELSxEZLISMWVjg7cR2Ulng0RDD68sFZPj3bIXGtr5LmjUnR1EV12Gf6ePw9jUDzCJGWMMDmJH15am2Yh8Hgg9davx+AsKQFRampEWlNSktxtQAQL5sQEDExKn54WAgEcm4et2WzyKnejo3ilp4u1XXs7yFZVJcgmRW8vrl8tZ1Oq5moFGExP47h2u7F7UGJuThiBSkpEzGw8lm4ySRkDMAb3wcaN4jOeXRGp2sQYBlBurrwAVmIiBtnkJL7LzcXA5ZU716zBNqtXC4upGlmIQL5z50QupdUKAqttzxjWxxs2iAp1/HzBICS5zSaMPoODWDMWFanfu9HiXe3t6gHyUlgsuAe1625vlz+/2lqiU6fCu2+WAyYpYwC1DIm0NOM5jxyjo1DnamvFQAsEhMQqL5dnczgc2JYxDMK5OcTDJifLAxB4RTmeMMyLWRnx5TU0II6V7ycFz96Ym8MzeOUVon379C3OPl94Unq9xsIDi4vxjJXqbzAoQhilWL8e9xOO7EsNk5SLjLY27YLGTieIZmQd2dEBkijJ3dUlH3RqKVs8eoYIFt6pKaJ/+zesATlKSiLP1ujrgw9US+JyJCTgnDt3QnXWIpTRWNyBAWPESUpSV3/b2tRLkiQkQHp3dISq/csJ0yWyiOjrg1FFLfSLSBBIDzMzkDL5+erJvsp8w+xsrF211qz5+dhm0yZE+XDCRkrIyUkQjVtaZ2a0j3HmDIwoW7Zg3amVAjY0FErKjAx5XC03AkULxuTahRJOJzQMRWniZYVJykXCyAjUtrw8/e145Ika3G5IwvXrtYmthupq7Rqso6OQVjt3QoUuKUG0T1NTZPV42tvlwQFer7orpKFBROPYbJgQCgsx0SivS81Vk58vAiQCAZB/IQWzOjq0KzdwFBYiMGMxqjcsBkxSLgJmZ0EoI9bNiorQMhY8QN1uj87kb7Nh4CoH1fQ0pDdXzVavxnXW1uI8bW36WRsczc2h1zU1FUqW1las65QWzbQ07H/ypLDoejzh15LNzXLjTDSYmzOmFVRUyI1dywmTlAsEt45y40c4KI0N58/DWFNbG359pSVdiEA4aZhZMCiSi/l509NBJq7OVVfjvH19uIfx8dDjDgxAxQvnOuBBEunp6t8nJiIO9cwZ3Ee443GXSzSWUa7Kd3drFxtTw/r1oRJ9OWCScoF4+WUMgtbW0FdbG9ZxyjIVJSUYxG1tUC3Xrw9vPCGCWpefr/29tKrdmTNylwyHmqpbXg5yjoyAnDxRem4Ok0a4igFuN0gnDYBXW+Py3iFPPQUVnj+nvr7QZ9TbGxmhOKSFvSYnjS8DJibwe1ityx+KF3ft1X/xi1/QRz7yEcrKyqKsrCzatWtXyPZf+tKXQlqw79mzJ9a3EYLpaUgihwPWPeWLq42dnfjBpYPw5z+HZJqe1q8wEAmcTqiwdXUgmZpxw2JBCphahYLiYuzn94OczzwTvvL50BC218sL9XoxETQ3g7if+IT8OTmdWPvx5/PTn+LZtLZCzY0k8CI5WcQH600m4+MwQPFzjo7iWjZsgFpuNB45Foi5S4S3V3/yySdp27Zt9Pjjj9Pu3bupqamJ8lWm/b/85S9044030o4dO8jhcNDDDz9M11xzDdXV1dEqifd5z5499PTTT8+/X46uWp2dGMQtLeppSlYrBqEyhWpoiOjLX8ZgKysTwepSrFplzMAxMCD3QY6PY4ArXSlSFBZCZdZK7UpPhzXSZsO9MQbSKVVnnw/SVZkdMjwMycklck6OWBtKI2s4UlNxDJ8PKvj11+N8JSWYZNRKSebk6IfJjY7KJ4qxMbnlm69z1VBVBdU/2oyehSLmbQu2bdtGV1xxBT3xxBNEhK5bq1evpm984xt09913h90/EAhQVlYWPfHEE3TzzTcTESTlyMgIvfjii1Fd02K0LRgbww9fXIxB29xsvJ9GYyN+8NZW9aByxiBNpUWgAgGU2ODd/vivlp8v93u2tAgVTjrolO0KJidBPKV/rr8f95aZiUFcWgqSut343GKBNbOhAeSqrYWKK61K8OtfE111FdHf/V3ovbe2inIeUrS1gcgOByaz4WF9o5fXq56aVlQE1T0vT260SUszlkbG0dkJUi+kPU1cti3g7dWlrdCNtFeXYnJykubm5ihbMa3/5S9/ofz8fKqpqaHbbruNhnRqVczMzJDP55O9FgrpmodHixgpcSh1JVRWqksBiwXE4u3seMu42lpIh8pK8bkyEMFigYTNydGvU8NVNKmznYfNVVaCaBs2YHIgwoCurgahOjqIjh4FKVtbcRyuiubnowX72rXq98aYnJCjo5Daq1ZhghoexrWHM/Dk5oYuFyoqoHUcPYq/KSniOUVCSCJoMNFW8lsoYkpKvfbqWq3QlbjrrruoqKhIRuw9e/bQL3/5S3r99dfp4YcfpjfffJM+/vGPU0DDnn348GFyOp3zr4V2cR4eDlWdKiuN1Un1euW+zLy88I5rjweDfcMGWE+NDJbsbJBGr+RkZSVIxevqFBWFN+pYrZAen/+8IIP0fnp6QK7sbHze1KR9rNZWUa3A4ZBbl6PR3wYG4Fa6/Xb4ZTMyhE82mnk4PX2ZfJcxaMs3j97eXkZE7O2335Z9fuedd7KtW7eG3f/w4cMsKyuLnTp1Sne7trY2RkTstddeU/1+enqajY6Ozr+6u7sX1J+yoUH988ZGee9HJQYG1Psunj2rvd/sbGgPyJERxs6cCd3H52PM7ZZ/1tbG2Ogo+lIqMTnJ2JtvMvb886HH4n0lJyYYO3dO/l1jI2N9fYyNjck/7+hgbHxcfq7xccZ4a9HBQcaGhhg7fx73rOypKb3P5ubQ69VDe7u4d7WemL29OH5zM2N+v/Hjav3WRhCX/SkX0l79scceo4ceeohee+012rRpk+62a9asodzcXGptbaWrr7465PvFbK/e36/tluCSR2stNDSkbjyoroaUVXOUNzeHxn06nZjFz56FmsV9g/39oTGea9ZAPbRYRIYHNypxq3F2NgxEUh8jl1QpKXLrcFcX1PbUVHnLuqkpqNhKFwQ34pw+DcsoYzim0hA1OaleMjMcGDNWg4cHKgSDiE7itYDCBTBkZ+N3izSReyGIy/bqjzzyCD344IN07Ngxuvzyy8Oep6enh4aGhuZbsMcSvD25Gmw2/NhqqtfAgHYInt0OY4/Sec+d32rrK6sVfsjBQRhh1MBJaLfDvdHQAPVU2np9chJqMa+zQxQaK5qUJIip5fvTC2dLSsJ3P/0p7kfN/9jdLQ+0T0sLrzrOzMCoU1NjvLqC1YqJqLoakxB31Wipt9KwvyVD9MLZGCJtr/7QQw+xxMRE9pvf/Ia53e7519iHutLY2Bj7zne+w44fP846OjrYa6+9xi677DJWVVXFpg32GI9Wrejuhoqoh9lZdVXRiBok6SDPJiagehrBwABUs8ZG7NPcLF4TE9imtZWx+nqoslJIVb2TJ6HGdnTgPqRobGSss1McjzGhYnZ3QyXlUN7/8DBjTU0415kzUGmlmJsLVTmDQX0VdnhY+5lG0tKdo7cX16im3g4Phy4LjCDacRZzUjLG2I9//GNWUlLCEhMT2datW9k777wz/93OnTvZLbfcMv++tLSUEVHI67777mOMMTY5OcmuueYalpeXx+x2OystLWX79++fJ7kRRPuwjK4vlNu53fhhw2FkBAOcMay5tDA7izVUS4t4/fWvjH3966Gk4+BEaWtjrL9ffC4dwH4/SKM2qBsbQWopJiZAVOX9SknZ14dtpOdqbmbM65UfW21NrUWunp7Qda6R/YwgEMD+TU0gKkc0a8tox5nZXt2g3tPRgfWHkaXpzIw8EJz7JY2goQHqZlERQr94GpP0V7LboQJKnfDNzdgmMRFqtDLRV9rmvK8PbgxeYFnqvxwdhUvh6qvlavPgINZ80qrnRAgz/NrX5NfCj9nZiXWry4XzdXeLZ9LdjessLNT28aq1Zm9pgeVbL07YSEt3IxgbE0uD9HTcQySG+2j9lGaSs0FIW48XFurHVEqTbY30jPT5YKSxWDBQv/99ogMHsGaMZHCVlIgA8rNnMdB5TRunE872zEwQfmgI92S1gojcdZKQgHMr17HKqB2Or38d5OPPht9rYyP+51XdlTVgV6/GOX/1K6Ibbwx/b8EgwgcrKhaWyhUJ0tPxYozonXeI3niD6Hvfi/15TVIagNdLdMUVMPAwhgHGZ1AuNWw2DDROAl7YanpaDFTG4MeThtQxhh++shJE8PuJvvlNWPt4CF24lDDeooDHfTqdMHy0tOCvywXJ0tYmSJKTA4nZ1YXzR1slzmqVG3jcblhDq6rkbRaUNWAHBjAxlJfDGkoEKV9aGjohTExAU+GdyZYKPT2ilXxeHtHnPgfpqZUJs1gwSWkAw8Py8pBqZvRAAD8ir3tDBGvjP/yDiJThkTp6oVt9fUhxamtDVMz0NJzf3GqoNii9XnF92dmipk11Nf6vr5enlnV3Y6DX1KhnkiwEhYVEn/wkBnNjIwbwqlVi8uKFtPLycP72dqEN8GLVRJhgLBZRsHrDhsW9Ti14PMLqy+OPR0YwOVZWQrswSRkHMDI722xCos3MgFR///cgUyCAwWpkWTE7C/W3tlbkadbUQILydaNe63UeRcPXXLm5IGpDA/b3+7EeXWBQU1gkJ+MefD7cR2MjnkNurlhfd3TItYCkJDG5eL1E//7vkJ6f+1xsr9XrhauLCC4j6XIjGMRkyyeFpZDUJinDYGoqsqDkzk4MfL72mZkBqfr6oNppFU5WgmeYcKmXkIDj8PKVgQBU5ImJUMc29yvy67Zasf8VVyx9ZfCMDLxSUjDhSO/d7w+t49rfD8mUmEi0Zw+uOzkZE4rFAtKEm9yMEGdkRKyjc3O11ff6enmgA59kjeS/RgszyTkM1JrJqIFXEMjPhyQ7d07UTvX7ofJyUjU2ikBvPRQUYOBIrZ0Wi2jCOjBA9PzzkEbSbcrKcH4izPSnTuEelrNUf0kJJNDp06FtGHhJzKYmkLSmRqjXRFgjV1eDOJOTIGg0naHHx0WwwOwsjlldrR0M0tmJ65aSvLwcEj6WMCWlAejNvLzsRkZGaDgcrxwnNfnn5eHF1Tq1quZS1NTolxv5+McxkNvaxGDn0ojnVsZLweG0NEidM2dA0MREEbBeXi5q6XR3Q71Wu2aXCy/GQBqen8kb+CgxNSW6g6WlGa/5MzyM4yqlstW68Gr34WCSUge8iK8W+vrgTqipkYelSeuZau3P1Tpe1dxqVVehbDb45ZQFiT0eENrlwsCW+vl4rmF3N5EkuSYuYLNhkvjP/4REV1pUJydBtHC1cS0WMZnNzYmE6owM+XuHI/LiW34/lhp6ieKxhKm+6kCr5ujkJFTV5GRIMGXZDWV8bG6udstzXtW8ooLo/feJ3nortJK6y4W1Fp+hR0bkrpaCAnn8Kw9YjzdCSvHJT+KvUuq0t4cvCamEzSaCJt57j+jIERx39WrtZkF6CFc1fdUq9V6YiwVTUupAWfiYCO4Nm037R1NTbXJy5BZRNfBs/+3b8Veam1lUBOI2NWGS6O+XS0ZeJYA75xsaFt/VEQts2ABVlicBGSUk7wLGk8p575P+fnQ4W7dOlIycnkagh9EiXO3t2qowR0pK5C0oIoFJSoPwemFYqajQD7XTypZQtiOXYngY0q+yEvsr67/29kI6ezxEv/890Re+EHqswkKhTvOej/EObrRqboaBzG4Pjdbhyc+Dg2LCS0jA5CS1gPb14Xm4XCIYgbtbJiaEOltcrB0R5PXifGlpi3aLUcEkpQZ4NTTuH8zJMdbPQq2AFhGIKnWUc3i9MMhwIpeVheZWcuf75CS6HdvtUK2V5UeOH4daGGlLguVESgpU/ZdeQhkRaScv3lx2chKVDrQmmr4+fFdYCI1DSarUVPE8u7pg/ElOlqu2c3MgvtH6vZmZsetxaZJSA2Nj+KHcbvxQRiSPXmNWmy20+nZ/P/yY0nWr1nnOnoW6Z7NB5Vu3Ti4pgkEcXxratlKQm0u0bVtoxNLgIJq88hb0ar0ppTWEiKB16KXVciJOTgrLb3ExJsxIooby8kRw/GLDJKUG3nyT6LrrIot86e3Vt/RlZoqgcLdbBAAowZvccAtkfT2Oy0m4YQP8fZdcIvaRVgFYiaiokPcrmZrCs6qqgsai5hbixi1p2KPRNgV2O54nY0SPP46JQavyAa/awP/n6OqKzTM3SamBv/s7kGxubvFKQRQUQBWemMCPq2V8cLmwndMpjDvSdazFggHKe1POzmI9pdVZaiXAboc05G6olhZhALJYQicqtxvbGg2ImJ4W/koiEHLNGnx+001QRf3+yGq9xspfaZJSA3Y71Kn+ftHrQ2/Qu93hU7T8frg8qqqwbWsrftjUVBBfaUBqa4MqpqYSJyYKtWtmZvkKBy8mqqpAxkAgdP2+ahWsyk4n1vvhCDk+LtaaRHi2agH9vHA1j7riy4TlhElKFUxMCFWmoEAEeWdlaRNvbEx7LTM0BINOQgLRlVfKVR7GoKoNDIiULosF+Xvd3UT/+I/6ZSJ5H4+VYG0NB6sVhEtJgSFLqjY6nVAxT56E1HQ4QGDlff/xjyBYRkb4oAEeOcSRkoJ9Imm7rqxjuxgwSamCgQH5j2W1Ql0cGkKibXW13Oig1pSUx3P6/bAucr+iskAxL54slYZNTUSf/jQGaCAgciLVwLNGLhTs2AHHvDT0kDGEJX7wAdG77xLdcANIWlgoJ8T588goGR4OX6Vudhaqq9JmkJgIScmJqacd5eaqN75dKFbwKiR2UAsaIBJukY4OeUSH1Ok9MQFSNTfDiFNTE76RrBR1dRgoTqeYDBIS8Lk0V1OKC0FKcvC1pRLd3egM/b/+lwgIUN43D6DgPls96JUMsdlE2z6tZ04kqjksNkxSRgiLRZQnrK+H6smrETQ1YbauqZGX4jACxqCaVVXJpSZjmIn5ZCA1VjAW++Do5cb4ONZ569ZBY9Aqx3nunLBkZ2fD5aHVzYz7oKWkzs2Vl5K0WEDMpiZRfUCJWE2GJilVYGSgO50YKH/5C9ELL8BhXVMTXU9Fv1+oS1JzvrSfBZ8MnE5MBpOTobVSLxTk58PA5nbj78aNcjWytFTeJ8Xvx3pcWjeJF8ZWIhjExKm0qGdmiiJlUmzYgN9AraFurGCSMkq43bAG/o//AWPM4KD+rKqFqSkc55JLQtcvaqocTxHr60MEz0K6QsUrnE6iP/0J60itgl3S5rBaPtqystDcR2VGjRGsXYvfOxaqqhpMQ0+EGBkR7g9p63K+PunqgtqUnx8+umZsDGtTveBxqdSemxODLCEB0S4XKnbsADl5YD63xFqtYs3tdsOqnZWlrkqmpcEQ4/NhMhsawrbR+HOrqmA7CARi38LAJKUK1H7gmRkQwukMjSyRbs/VSY8Hs3Jmpnonq/PnYSUMF2s5Pk709tswFiUkYHDw86mpZxcScnNDLZu8QNnsLFTbH/2I6I478CwyMvCcpL9HaSnWpOvX4zdZSI7kmjWYdP1+8ZvGYl0Zd+3ViYheeOEFqq2tJYfDQRs3bqSXX35Z9j1jjA4dOkSFhYWUnJxMu3btohYjfeiiAGMgV28vHPRG25W4XFCTeHa9dA3ES0dy1SwQwPe81bf0ZbGI2qNlZReWpVUPWuv6uTlMkLy79He+A8NYZSVU+fb20GdotRI98YSxkLhw9oSSEpCSW3djYmiLvBh7ZDh69ChLTExkTz31FKurq2P79+9nmZmZrF9aO1+Ct956i9lsNvbII4+w+vp6ds899zC73c7OnDkzv81DDz3EnE4ne/HFF9mpU6fYpz71KVZeXs6mpqYMXVO4cvK8hwUvyR+udZqRMvkTEyjPf++9jP35z/J+H62toW3hlMeem0PbgK6u0Ou8EKF8pufOoZWAtF0BbyXQ2irvcaKG//f/0KtlclJ7m9FRxox2v/B4cC16v33c9hLZunUru/322+ffBwIBVlRUxA4fPqy6/XXXXcf27t0r+2zbtm3sK1/5CmOMsWAwyFwuF3v00Ufnvx8ZGWFJSUnsueeeM3RNeg9raoqx997DD65sRKMFo70rurpw7BMnQFAjv5WSeD4fyBmun8ZKR0sL7rWxEWRUkqmvT95U6PTp8MdjDL1UFqsPyeAgY6++qt1bNFpSxl179ePHj8u2JyLavXv3/PYdHR3k8Xhk2zidTtq2bZvmMSNprz44CFdDaal+a4JI0dUF/9qWLVj78IptStU2HNLTYUQ6ftxYX5OVipkZxAkXF8MVpMzgGB2VG9Kqq7W7RkvLhK5ZA2MP9zFHC8ZgqBsc1PaHRou4a6/u8Xh0t+d/IzlmJO3VV62Cb8rjEcWPw4FnMGihsxODijeb5ZEghYUgZ26uaAOuLJsoXUMOD8N9EgiIEvoXKpKSUPd1aEgUcubo7w81niUl4bmqtarv7ZUHr6enw63kdsuDMYxgdBTjoqUFx7ziiuia3erhovBTHjx4kEZHR+df3Tq/BDeXl5djVjVSozUvT7uxaEcHBoE01K6wEKTnSEkBOaurMYCamsTxGIOltqEBRo61a0HulZymFQl4qGJbm3AHaWX8u1z4Tq0erJqBbM0aELmuTl/a8QLYzc2QrrxebKwqPMT0p42mvbrL5dLdnv+N5JhJSUmUkZEhexlBQgJIkJICUmhpvVoW0bY2/Ohqfi016crLJvKSla+/jq5U58/jOtRcKxcqpM9UGs30xBP4HZSlUDh4gTGjyMiAm6S3V0hNrh2dPy+kYmkprkE5xFacSySa9urbt2+XbU9E9Oqrr85vX15eTi6XS7aNz+ejEydO6LZsjwRKlTUzE6Tw+YT6GA6trSCjVvXtoqLQtuiMYX3Z3Ay1bcsWon37oJq1tODV2Xnhx7tKIa1qPjGB4tOXXYbn29ysXlVuzZrQbJxwqKgA6QcHiX73Oxx7bk5IRa3eLbH4LWIePHDHHXfQLbfcQpdffjlt3bqVHn/8cZqYmKB9+/YREdHNN99Mq1atosOHDxMR0be+9S3auXMn/fCHP6S9e/fS0aNH6b333qOf//znRERksVjo29/+Nv3zP/8zVVVVUXl5Od17771UVFREn/nMZ2J6L8XF+BGam2E40GpR19wMqRauoDCXlsPDIhuhpEQYJdrb8V46IGZnRXL02bMXVtqWFKdOYfKTFr0iwr0nJIjPOGmJsETIysI+Dgeea7gaOlNT8oyfQADLAyMRWbFCzEl5/fXX0+DgIB06dIg8Hg9t3ryZjh07Nm+o6erqIqtkgbRjxw569tln6Z577qHvfe97VFVVRS+++CJtkKSDf/e736WJiQk6cOAAjYyM0JVXXknHjh0jxxIEgloswnLKe4dIo04aGyEF1TTksbFQQ8SzzxJ96UvqOXlqKWSJiRhoAwMXRrUBNTAGaTczox37yiFtRTAwIAhaWgppOTUlVzlHRuTrf4dDXpGAMRyvs1Ne8HopYbZXV2FPJO25+/sxI/MM9OJiqLlq5va0tNCUodZWBJ5bLPot0YmwjmpshFpcUIDzTE4uz8CJJTo7YQxLTMTEV1wstA5lO3g18GXA9DTR//7fyMHkcDqFFVwNzc2C5B4PjqFWJZ9Db6yY7dWXCQUFGAT/+q+ozm2zGe9wxf1nxcWYwXneIFccpNOl2w1VV1ruMiMDk8KFRkreo5MIRphz56Cmhnum/f3CTcQYlgT/+I94xkaqEvI8Sw6XC79LY+PSaiUmKVVgRHdgDCb62Vn8kPfcg4EzPg4rXmVl+B6GHo+QjpmZIleyqEishXgDWpfrwlVXlVBaNEtLhb+SJ47zRq/SbfPy5KQ6cwapdVyNDedP9PlCJ7jMTJD61KnQvE7e8WuxYZJSBenp2r3t/X6RnVFeLo+q6e9H0HMwiIEQDOK9lk9R2XzUYoFk6OmBSsxN9MrKbsEgBgR/DQzoq2QrCVwq8mB86QQZDBL9n/9DdO21WNfrBZg3Noq8yYqKhVWpcziQgM4Na5zcakEMiwGTlCrIy8PgkJLS54MKyTsqq/mn+ADibe2CQZEPqEdOJfx+ELOrSxgspLBYIDHsdkiQyckLh5T9/aK8pxpuuQXWVY8HZFGTfr29MJxJnfu82atWL9CBAf1aShYLJGVzs7DyTk8vfjQPkUlKVdhswjnt8WBdweNVI4HVin2CQSFdpfmQUgQCwtVRVobXZZeJpkJqaG8n2rwZa02eyLuSMTyMSSgrS7vMI0/ZcrlAspkZGGb4hDc2hiWFcv2ZmorfVes5jYwY62NZXS2MSLHCRRKsFTnOnhWdlmtrw5csJApVtzisVvyYvKszj6m1WCDlGhshDauqcC7u2UlPR/ABr9MjRW8vvk9Ph6Eolq3ZlgoDA1AHExOhxp86pb++Ly/Hc21uBkGDQWg4WtKwpCSy4H8tlJSA4KdOLfxYajBJqYH8fJBGWiQ5HAoLQ6N0pLDZIDkrKohOnEAI3fHjIhNCTb3lznCpg3toCANQqm7l5ODzlYr+frkKnpAAdfHUKXlInVJyWq2YyOx2WMATE9VLVHLwKuxSDA0ZL/Fx/jz2Hx2NrM9MJDDVVw3k5gpDQmcnVKLERP3s/9TU8MHr/f1QlcrL0cOiogIqslqCS1YWBktentgmLQ1qnrLKejBI9Ne/IntkJeKtt9DmTwreiv30aRCUN+SRwuMBQbKzUcC6ogJSMxCAlqP0/SYl4XcaHhYhkEND+qrr2Jj4fTIz8ewDAXUNZjFgklIDJSWihCNXh2ZmxCybk6M+u2oRtqcHUUD5+WJtOjYm+iqqYXhYXoenoQED7tprxecWiygZUlm5MrtvNTRgLdndjcEubVWQn49Kf6dPy63Q3d1Q/QsKhBtjdBTPghfGnp4WET7Z2SJqqqgImSFaBbeI5OF3aWmhz7SzUz+oYCEwSakBhyNUbU1KEjOq1yvWhtKuWNKZnDH8eDMzUFEjrQmbnY3X+DgGyGWXIXhAzTJ59iykitcLaW1kDRwPOHcOpEpJwT1K721kREyCKSnoYj0zA+IWF4eqj0qCORzi9xoelndz5tkkBQUixpVXC7RYsK/e5BYMhvdDRwuTlFGCV1pjDAOLq7dEwtrKGKRstHl3k5MwTKSmisABtQp2HR1i1s7NxT68D2Y8gzc94oESfr+8E3ZmpriH8XHcp8MB7UXNFcErxqtJPz7BEYH8k5NYH3o8sNSOjOB30rKOLyVMQ48OsrPDG08sFhCiuhpq6NtvEz34oKiYHg0hp6eh0g0Ogox6BgWfD8YOaUvxkhLRJTpewUkhdV1UV6tPOt3deBZlZUjd6uzEUkCJ9HRjlcyLiyGRnU6il18WGpDRaoG9vbHVRExS6iAnB2qPEfD15qZNRNdfD5Jw94dWVQLlAJidBRl5S3e11DBl8xoeYMDBGK5jfBwDTqcc0bLh/Hmi//xPqIstLaLVn8UCrYDHrwaDUMvT00EYvjRYuxbEUJZDycoy9nt1dOB3qa0l2r8ffsuzZ43lyRJhQlnM+k1KmOrrIsDtxgDZsEHUbhkcVF/P5OaGJj77/djPbsdA0ZutU1OF26WpSZxjakr44LgF8rOfxd94yiTp7cXgv+46+GYrK+Vrx/R0qJf5+aJ4ssUC6Sj1P3KpGggIFTcxUb9LlseDCaGsTKi/eXmYDNevh784NzeyLmmxgEnKMEhKgjqplqrJGCRbfr7cpJ6SAiLw9Y10PTM4KAg6Nob9ExLCk1GJoSHM8D4f/nc41COOysuhykrb9S0XWlogzbgVlEs+6drR5yP6wx9QdWHzZrGvdM3OUVkJYodrJTAyAu1C2mqCIzMT15Wfj+94ZfvqavXfg0d3xRImKcNg9WrMyEpL3PnzmPXXrhVWOKnllQ8YpaU0L0/MxC++KCSBVvkK6TEdDpA7GCT6298gEXNywof/FRTAGPXCC0R796q3a48lxsdhOf3oR/UbrPKInOuvlwdh+P3als6KCkhRZTAFEZYU7e2iKZIRuFx4xnV1WBYokxIGB2PvcjJJGQZqs2VbGySoMutAWsrQbsdAUWaCSLFxY/iEXSmmpzGTP/cc0d//vfZ20usYHsaaLTub6NJLMagmJzFZRNI/Mxrw/itpaZB8FgtUxMJCUVDM6wVRfT6o3xs2hD7zjg7951RWhn15VQdesoUXPosUvJvzuXPQQmLlj9SCaegxCMawbjt7FoNKzec4NSWXQmohXQvB1BQG+qFDQmXmA7C5GdKbCKSbnYUaNjcH1ZiHsJWW4v25c/heqyrcQhAIgHy8/wp/VoWFeD85iXMnJAiNY3hYnZBEotuWHkpKcK8vv4xzV1QYn/DS0tRr6JaWYuI4exbPM1b5k0qYpDSAsjKEgfX2YuAYVf8sFsy6Wq6JSAqx9PZCDVSmjeXmioprc3Mg57FjkKo1NSLfb3ZWSEbeus9qJTp6VFTK04sZDYfZWXGco0ehSSjXsPx+eRHqsTHE/87NhUqjhAR8PjkZPj2quxv3PTkJ67fdHhl5CgvVizgTQX3ljWPfeWdppKapvhqA3Y6BEU3lOF7QeSFVA1paoH6GC5rOzwexPvc5qIXSekB9fUKl5Qaq1auJtm7FfQWDGHh+v2jbkJCAYyonofFxqMR8DVhSInIgLRaQPTtbft+5uVAFpUaeoSG4JIaGIDlXrRL+1tJSrAeDQfW41L4+4ZMsLhbXaLOBZKdPY72+WFE31dVQsZcisMAkpUFUVMD4oBWnOjOjHSjA1SO1SgZ6CAZhcJBmu4eDzyfUVT4wExJANt4Zur4eay27XUgIqxUTiNeLdWdlJbbt7wcBmppEWc3UVPy1WOC+yM0NLeXodEJanj4NSZOZibV4bi6eVVOTIM3EBCRnT49YD6elhUbn9PeLAtaFhdoO/JQUrNeVRbf0EE5r6eqKbn0aDUxSGkRmJowsWqTs7tZ2ORQXh0rLcINgfBxSSGudpQavV0hTux2qXH29CC6YmIA03LhRe+3W34/r5OlMXLpy1baiQuzb3k505ZUgmFp9VYcD119XJyyWPh+It2mT2I7Xv+Vrz+5ubMPdGNyFVFBgvPwGL63S2Wms6BZXl9WMX4zh2SkzTmIFc00ZAfSSZMMZI5TRJnqxqQMDeHF3iRrUEqq9XrlbgA9MjwcqsNutf0xeXlEZycRjaysq5K4bHqeq/FwKqxWTwLlzOL7XGypxlL1VVq+G5Gxvx3XwZGYjEk+JsjJMDuFaGZSWav+27e2RWckXCpOUESAlBTNmNJVyCwpEOBkRBpvaGvHcOait4Rz9Lpd8IPt82uoxL+QljY9VYnQUaqmapJiZwcDm383NhfpPuRpKpP580tJEdQA1OByhtXJ37iTasQPPpKUl+hYBOTkg5+nT6s1/iDB5qB0/GMTkE6tmPqrXEqsDDw8P00033UQZGRmUmZlJt956K43rRAsPDw/TN77xDaqpqaHk5GQqKSmhb37zmzSq6IJjsVhCXkePHo3VbYSABwVEg4ICuS9NKrG48SUz01hIXGqqPChbasjhOH8equOqVQgaCAblFQyIxDX09oa6eYJBSAlpeBuXisqwt7IykEcNXV2QqNdcg3Xg2bOhMblakopbisvLsQSItHUdR1ISVOa2NuPxzETLk58aM1LedNNNVFdXR6+++ir9/ve/p//6r/+iAwcOaG7f19dHfX199Nhjj9HZs2fpyJEjdOzYMbr11ltDtn366afJ7XbPv2LdQ0QKux1+OOmMb7ScRFaW8CVKMTsrDDrRqGjSpqhEIHhTE1wEvKNUWRkIkZgYSh5p6hcHJ9nMjLyMJk+oHh4OlarFxaGkaW/HBFJQgGP292OdOTISmhFis4m1q7K6HA8EyMwUGTTRYO1aTGZ6yxCOuTlc05K3HYyo77NB1NfXMyJi77777vxnr7zyCrNYLKy3t9fwcZ5//nmWmJjI5ubm5j8jIvbb3/52QdcXbdtrjkAAbb85lC3Q9TA2hjbrvI33+fPyY0UCfoyGBvHZ0BBjZ88yNjsrPmtqku83NMRYayv+P32asfZ29eP/9re4tuZm8Wppwd8jR3BcZevz5mbG6urE/yMj8mNKn9XUFLbl2wSD4lmEe6YDA2gzPzYm/5zfVzgMD+Pc0tbo09Py1usNDdqt040g2nEWE+vr8ePHKTMzky6//PL5z3bt2kVWq5VOnDhBn/3sZw0dh/dgSFB4gm+//Xb68pe/TGvWrKGvfvWrtG/fPrLomChnZmZoRuLB12uvbgRWK2ZutSDpcEhLgwqZkACV0++PvHSlFIGAWA81N0OKrV+vv092NiRAUxPU0J07sS8vLcLR3Ez0iU+EOuK56l1Zif97ekTJRZsNUt9igcFGbx3rcCAmtadHuGH48cOBxxB3dUE6V1VFFjCQlYVrO3MGRiSHAxoBX3NOTsKGsBwJzzEhpcfjoXxFdeCEhATKzs7WbIGuhNfrpQcffDBE5X3ggQfoYx/7GKWkpNCf/vQn+trXvkbj4+P0zW9+U/NYhw8fpvvvvz/yG9HBmjUY1NEEBeTkIFTu1luF4z5aFamlBYOzvh7kVg7MwcHQQO1AAOohb7EwNRWacT8xgT4c3d2hJRul5UY4+Tj8fqKf/ATPx+PBX+m9paeH1l4tLoaaXF8PNTeSNXtJCUgsrVlkFNxt1NSE30QaLN/VtXxtIiIi5d13300PP/yw7jYNDQ0LuiAiSLK9e/fSunXr6Pvf/77su3vvvXf+/0svvZQmJibo0Ucf1SXlwYMH6Y477pAdf/UC6wNaLJhJefErNYyPq1epS04m+sxnkJo0MiKyHLTqxmqdf3QUA3nXLm3pODIiDBW8iLDVCoON1QojkNUaeg+9vZAgaq6EqSltv197O9Ftt8HN4XCAYIxhXblqFYxYLS2h6U9JSZCavb0Iady6FRMCV3DUijIrn8fYGNEPf0h0333y1nZWK4IptKzTNTWYaHh1upGR6Nb2i4WISPlP//RP9KUvfUl3mzVr1pDL5aIBqf2fiPx+Pw0PD2u2QOcYGxujPXv2UHp6Ov32t78le5hUhm3bttGDDz5IMzMzlCS1SEiQlJSk+d1CUFKCONPt20XguXTwpKbKne0cTU2QIDabqPUTKYJBxGJyK6xaGY3UVEhEPoBXrw6NDHK5YNXMzFR3h/CSKEYMWe3tGPzJyTjv6tViQpiYEEEAXi+ubXJS/TjcMJSaiuWBEek3PIxncugQpLXUFxoM4nq04lul8Hjgz/3IR8JvGytERMq8vDzKM5CWvX37dhoZGaH333+ftmzZQkREf/7znykYDNK2bds09/P5fLR7925KSkqi3/3ud4aawJ48eZKysrJiQjojKC4GIXlqklGsXo21VLQCe2YG0ig/Xzsmt7cX6itfSWgFnNfUwE2xcSPeS5O68/KEesehJs27uyH9eECE0h/IOzJPTGAyqa7Wr3Oj1ZZeDdzNw6sRpKbCys2jjKxWY24mnw/3sewFx6K3Leljz5497NJLL2UnTpxgf/3rX1lVVRW78cYb57/v6elhNTU17MSJE4wxWKq2bdvGNm7cyFpbW5nb7Z5/+f1+xhhjv/vd79gvfvELdubMGdbS0sJ++tOfspSUFHbo0KGIrm2h1lclOjtDrYxa8HoZGxzE/5FYbZWor8dfboFVg9Kq29srrKgej/y78XHG2trwf1OT3OrY0sKYxAAect1uN15SqF3X5KS4bqnFWAm9e1KiqQnPVLnv2bPGj8Hh9TLW3R35flqIdpzFjJRDQ0PsxhtvZGlpaSwjI4Pt27ePjUns1x0dHYyI2BtvvMEYY+yNN95gRKT66ujoYIzBrbJ582aWlpbGUlNT2SWXXMKefPJJFggEIrq2xSYlYxhkMzPht4vWlSLF4CBcAoxpD+DJSUwWWhgZEQRtawMJu7vhLlG6UIJB+WfS6x4chItHCeV1TU/LiaIkut6+avD74c6RumOmpsS1TE4ad4/w7ZX3vVDEHSnjGbEgJWMYJHqQ+uEYw6D5UAmICFIp09Ki7kuLxPc5Oyv8j//2b+pShg/Y8XFIXMbgY+XSVQkpsWZnGTt1Sv59MKguLWdnGftwDtbE2BiuUXnfbW3y59nRgesNh2Aw9PoWA9GOMzP2dRFRW4toEy2ohaZF2o+is1NeUlKtqZDfH5mLxW4XuZilpbCENjTIK8MVFWHtxtPXJiZgPAkXoxsIwG/J16scUuuoFN3d+utsvWB9ZdVyvdA/Kerqwvt2lxImKRcRdjsGrxbR5ubkwQbSPphG4PfDgCK1oKpZMdUKdmlhcBDW17ExWCzLyrBvbS1cKI2NOGdamgjGn5vDYA/XzzEYhHP+kkvUjWCVlaFWY712AEaD9aVQC/2Toq0N9xyrFgTRwMynXGQ4nUKKSOMnjJS1CIeWlvAObS55wlmCeTHj3Fz1Y1oscOcwBlfH7CwsskNDIKo0H1LrOk6d0iYkEQIdjBRA5sH6q1ZF7j/MyIB0V8uVdLvhu9SLOloOmKSMAYqKhGmeV9Lu6VGXLHp1ZaUYHcWADEe2cLl/vPtxUVH4xF+iUHL+678S3X9/+KibV18l+tjHwgdDlJZCAqpVgyfCZNDYKColRIPqapBaWmZydBTPYakr1RmBScoYobISqtv69frrO626skq43cbCvqQNcjiCQZGPWF4uz/qIBD4fnPPl5frqnt9PtHs3tmtvF1XM1aoGJCdrtyrnRZTDSeVwsFigtfT34xrm5jBJxtM6UgqTlDHEhg0gZn6+dhmLxWgok5wM9djrlRtJZmch0axWkD6c8YeH+inBVVHeuLW5WT+Ivr1dxNLyte358whCYAxrQuna2uXCpGOziTjdvj6Qx2gR5XDIzUVIYn4+/i6U6LGEScoYgg/KP/8Z2RbRgDFIKD1Vs6hIlO/n8bi9vRj4kbRDkBZxloJLfKkE1mo5p/VdVhZejEGFnpvDWo6vEzkpq6qgOTidkXW20qu8wFFTg6r0e/cuf7s7PZikjDFSUoiuugoDW6sIltOpXbOnpSW8ldNiQcxmURGkkbSfZSTgZT+kOHsWg1m6nluzBsRSs4KGk+rSTstjYyLgPSUF+05P4/tIu1oNDOivpXllwD17lra0RzQwXSJLgLQ0qGFnzoTWoSGCSqWWST85iQEUzlzv8aCosdsN9TXSjtFa4NUQlGtQu127u9X4uPFSmunpIHx1NST94cN4FtG2mdPL1uG+0oVawJcCpqRcIthsWMc0NyO420jWRbicvr4+qG0FBURf/arI6OjthZqYlRV9W7f6ehhqtKzCq1aFqrsTE5E1DxoYgIbAj/eDH0BFbmzEtRstJxnuHGNjocEL8QyTlEuM6mpY/pSRK0rXwcCAdkpXTw9mf2lB4qEhqH7l5ULVHB4WKWWJiUg1M7KWamyEi0KPYGlpIKXyusJVUZiehgvEYsGEUV0N7cHthiGKT1jcMJSRoV1rNxw6OyHll7I85GLAJOUyoLgYhGlqEoO4oABqKE8xGhoKrY967pxIMJaqqLzIVEaG3B8n7Ys5MyOiZxiDf05tbdXcjOMbUSFzc0XXLD1/JGO49tlZkERJ3I4OkeDMy3Bww9DoKCaJtDTjavlCgg3iASYplwnZ2Rh8vLWA0wk/msuFQSqNkeUO/5ISdenl8wnXQ16eegmQpCThC+Uk4fVb8/MhuQYGYIwyuibMycHEkpsbes1EmHi8XvzPY2qVcLvFRMStyNIQQacTr/FxkDMlRb9S+WIEGyw3TFIuIxwOrHV4sx0iEMXvxwBuaxMdqQzkexORqNeTm6tfXV0ayTIwQPSb32D9OjgoqhVkZYF4eiovr7nD43r9fvgpibC/nuWYMaipXLLrnSctTd5GLylJ3AOX0iMj0Dbi2QdpBCYplxkWCwZlRwfIwCNyeMmQaMz33NdntIhwfj5C4oJBsf5iDKqjWjgdJ4HDgfXeu+/CqtncLHyNRtau0RQeS0mB+js9DYlot2Ni8PtFc6OVDpOUcQIeuvbcc0Tf/nb0oXBEGKiJibCGGlkbStey3KJqscBvqlcaY3oaBqs33yS6+ebInP0jIyBTtFX8HA4hOX/1K6LPfjZ6S3O8wfRTxhFKSoi++124Qoz2ztDahgd6G8HoqCCgz2e8op7DAbXy1lvVK7/rQa3NQiQIBiFp+/rQ4/JCISSRScq4g8UC9a+8XBRL1oNWaByR8CXqQRmeFmlLeN7pa80asZYMB2WithSZmeEJ3tYmenzwRrUXEkxSxikSEqCeFRRg7aRFrulp7SgVpxOOc72cRbdbrnYmJED9VYs8UkLqa01OBjmkTYfUMDsLo5CWWs07PquBJ13zVnlL3uNjiXCB3taFg+RkkDM9HQNSUU43LGpqRL1VJXgDGyXKy8NLaMbgppAmCBvZr7U18jb1Hg/uPScHzyLeY1cXCpOUKwQZGRiQNhsGqNE1nMUS2rCWo61NO9olN1e/s5VWizi9xro8yMGousmrHCQl4d6jjYldaTBJucLApcXMDAaskV6LLldodfBwZUPy8oTjX4nZWdHkSIn0dFybshgzY+pBDWrwerGeDgZxr2qt2y9kWBiLtj/uyoXP5yOn0znf1Wslo69PrOPS07UrgU9NQQ3kUTft7Qhb01MFlftwNDSEz9Osr5cnKDc3y+NylejtFfeRkWGsonm8I9pxZvopVzikRhqfT1hObTZEvHBjSHIypBWvB6SsrKeG5GQRYcSl4ugoSBNOBeXVBAoLsfZMSpITMhCQt1tftWphLpILCXHTXp2I6Kqrrgppnf7Vr35Vtk1XVxft3buXUlJSKD8/n+68807yazXJuMiQkYF1XlUV1nbt7SBpczOc7Dw52eMxnhaldJG43cbIk52NSSIQwBqztBTkbGnBi9fAra7G62JZLxpBzCTlTTfdRG63m1599VWam5ujffv20YEDB+jZZ5/V3W///v30wAMPzL9PkURgBwIB2rt3L7lcLnr77bfJ7XbTzTffTHa7nX7wgx/E6lZWJBIS5FbO7m6oiD4f0fHjRDfdZOw4FouojDAzIy+bGQ4lJUS/+AXRFVeAiKmpxkP/LmbEZE3Z0NBA69ato3fffXe+m/OxY8fo2muvpZ6eHirSiMe66qqraPPmzfT444+rfv/KK6/QJz7xCerr66OCD6f6J598ku666y4aHBykRIO28gtpTRkNpqdhCeXGGKXRh79PSIBETUkRld/XroXEGxiA6ikdPcr9k5NB4gvdhaGFuFpTLqS9+jPPPEO/+tWvyOVy0Sc/+Um6995756Xl8ePHaePGjfOEJCLavXs33XbbbVRXV0eXXnqp6jEXu736SofDoZ/+xDE3B/L19UFSOhxC4pWURNbO3IRxxFV79S984QtUWlpKRUVFdPr0abrrrruoqamJ/uM//mP+uAWKxRB/r3fcWLRXvxhgt4v1Y6QOfxPRI67aqx84cGD+/40bN1JhYSFdffXV1NbWRhULqOkQi/bqJkzECnHXXl0K3vW5tbWVKioqyOVy0d/+9jfZNv0fesX1jhur9uomTMQCcdVeXYmTJ08SEVHhh5WTtm/fTv/yL/9CAwMD8+rxq6++ShkZGbRusUppmzCx3Fj8VplApO3VW1tb2QMPPMDee+891tHRwV566SW2Zs0a9tGPfnR+H7/fzzZs2MCuueYadvLkSXbs2DGWl5fHDh48GNG1xapprAkTUkQ7zmJmP3vmmWfo61//Ol199dVktVrp85//PP3oRz+a/35ubo6amppo8sPmiomJifTaa6/R448/ThMTE7R69Wr6/Oc/T/fcc8/8PjabjX7/+9/TbbfdRtu3b6fU1FS65ZZbZH5NI2Af2uwvdiusidiCjy8Wodfxoox97enpMQ09JpYM3d3dVBxB2fqLkpTBYJD6+vooPT2dLJIgTm6V7e7uviiDCrRgPhd1hHsujDEaGxujoqIiskaQkX1Run+tVqvuzJWRkWEOPhWYz0Udes/FGUU1aDOf0oSJOINJShMm4gwmKSVISkqi++67zww0UMB8LuqI1XO5KA09JkzEM0xJacJEnMEkpQkTcQaTlCZMxBlMUpowEWcwSWnCRJzhoidlrKrurTT85Cc/obKyMnI4HLRt27aQvFUlXnjhBaqtrSWHw0EbN26kl19+eYmudGkRyXM5cuRIyLhwGO32K8XiJqusPOzZs4ddcskl7J133mH//d//zSorK2UpZmrYuXMn279/P3O73fOvlZwGdvToUZaYmMieeuopVldXx/bv388yMzNZf3+/6vZvvfUWs9ls7JFHHmH19fXsnnvuYXa7nZ05c2aJrzy2iPS5PP300ywjI0M2LjweT8TnvahJWV9fz4iIvfvuu/OfvfLKK8xisbDe3l7N/Xbu3Mm+9a1vLcEVLg22bt3Kbr/99vn3gUCAFRUVscOHD6tuf91117G9e/fKPtu2bRv7yle+EtPrXGpE+lyefvpp5nQ6F3zei1p9DVd1Tw/PPPMM5ebm0oYNG+jgwYPzeaErDbOzs/T+++/Trl275j+zWq20a9cuOn78uOo+x48fl21PhKqCWtuvRETzXIiIxsfHqbS0lFavXk2f/vSnqa6uLuJzX5RZIhyxqrq3kuD1eikQCKhWCWxsbFTdR6uqoN4zW2mI5rnU1NTQU089RZs2baLR0VF67LHHaMeOHVRXVxdRPuUFScqVWnXPxMrG9u3bafv27fPvd+zYQWvXrqWf/exn9OCDDxo+zgVJyuWuureSkJubSzabbb4qIEd/f7/mM3C5XBFtvxIRzXNRwm6306WXXkqtra0RnfuCXFPm5eVRbW2t7isxMVFWdY9jMarurSQkJibSli1b6PXXX5//LBgM0uuvvy6b9aXYvn27bHsiVBXU2n4lIprnokQgEKAzZ85EPi4WbCpa4YhF1b2VhqNHj7KkpCR25MgRVl9fzw4cOMAyMzPnzflf/OIX2d133z2//VtvvcUSEhLYY489xhoaGth99913wbpEInku999/P/vjH//I2tra2Pvvv89uuOEG5nA4WF1dXUTnvehJOTQ0xG688UaWlpbGMjIy2L59+9jY2Nj89x0dHYyI2BtvvMEYY6yrq4t99KMfZdnZ2SwpKYlVVlayO++8c0X7KRlj7Mc//jErKSlhiYmJbOvWreydd96Z/27nzp3slltukW3//PPPs+rqapaYmMjWr1/P/vCHPyzxFS8NInku3/72t+e3LSgoYNdeey374IMPIj6nmU9pwkSc4YJcU5owsZJhktKEiTiDSUoTJuIMJilNmIgzmKQ0YSLOYJLShIk4g0lKEybiDCYpTZiIM5ikNGEizmCS0oSJOINJShMm4gz/H9JgJRO/gLJXAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from sigmaepsilon.mesh.domains.section import LineSection\n", "from sigmaepsilon.mesh.plotting.mpl import triplot_mpl_mesh\n", "import matplotlib.pyplot as plt\n", "\n", "section = LineSection(\"CHS\", d=1.0, t=0.3, n=32, mesh_params=dict(n_max=20))\n", "triobj = section.trimesh(order=1).to_triobj()\n", "fig, ax = plt.subplots(figsize=(4, 2))\n", "triplot_mpl_mesh(triobj, fig=fig, ax=ax, lw=0.1)" ] } ], "metadata": { "kernelspec": { "display_name": "sigmaepsilon-mesh-HBhBZswl-py3.10", "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.11" } }, "nbformat": 4, "nbformat_minor": 2 }