{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# MTSAT-1R/2 HRIT Reader Introduction"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Supported Files"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"MTSAT-1R JAMI and MTSAT-2 IMAGER data in [JMA HRIT format](https://www.jma.go.jp/jma/jma-eng/satellite/introduction/4_2HRIT.pdf). Both __segmented__ (data split into multiple files) and __non-segmented__ (data in one file) images, as well as __full disk__ and __half disk__ (northern/southern hemisphere) images are supported.\n",
"\n",
"Example filenames:\n",
"\n",
"- Segmented: `IMG_DK01IR1_201006150332_001`\n",
"- Non-Segmented: `IMG_DK02IR1_20140524` or `HRIT_MTSAT2_20140524_2201_DK02IR1`\n",
"\n",
"Full disk scans are identified by ``DK01`` , half disk scans by ``DK02/03`` (northern/southern hemisphere). \n",
"\n",
"Sample data can be found here: https://www.data.jma.go.jp/mscweb/en/operation/sample/hrit.html"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Calibration"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"You can either read the raw detector counts, or calibrate to Reflectance (VIS channel) or Brightness Temperature (IR channels). Radiance is not supported."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Examples"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The ``SatPy`` reader names are ``jami_hrit`` for JAMI@MTSAT-1R and ``mtsat2-imager_hrit`` for IMAGER@MTSAT-2. The following examples are for MTSAT-2, but you only have to change the reader name to adapt them to MTSAT-1R."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Mosaic"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import glob\n",
"import satpy\n",
"\n",
"channels = ['VIS', 'IR4', 'IR3', 'IR1', 'IR2'] # sorted by increasing wavelength\n",
"\n",
"# Create full disk scene of segmented files\n",
"filenames = glob.glob('IMG_DK01???_201006150332_???')\n",
"scene = satpy.Scene(filenames=filenames, reader='mtsat2-imager_hrit')\n",
"\n",
"# Save quicklook of each channel to file\n",
"for ch in channels:\n",
" scene.load([ch])\n",
" scene.save_dataset(ch, filename=ch+'.png')\n",
" del scene[ch]"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Then run the following command to create a mosaic from the `.png` files:\n",
"\n",
" montage VIS.png IR4.png IR3.png IR1.png IR2.png -geometry 256x256 -background black mtsat2_mosaic.jpg"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
""
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Half Disk Images"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Half disk images are supported, too! Information about the covered area are provided in the dataset's area definition:"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Area ID: NH\n",
"Description: Northern Hemisphere\n",
"Projection ID: geosmsg\n",
"Projection: {'a': '6378169.0', 'b': '6356583.8', 'h': '35785831.0', 'lon_0': '145.0', 'proj': 'geos', 'units': 'm'}\n",
"Number of columns: 2750\n",
"Number of rows: 1375\n",
"Area extent: (-5498000.088960204, -198000.00320373234, 5502000.089024927, 5302000.085788833)"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"scene_nh = satpy.Scene(filenames=['HRIT_MTSAT2_20140524_0101_DK02IR1'], reader='mtsat2-imger_hrit')\n",
"scene_nh.load(['IR1'], calibration='brightness_temperature')\n",
"scene_nh['IR1'].area"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Quicklook using cartopy"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDBkSEw8UHRofHh0aHBwgJC4nICIsIxwcKDcpLDAxNDQ0Hyc5PTgyPC4zNDL/2wBDAQkJCQwLDBgNDRgyIRwhMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjL/wAARCAC8AWQDASIAAhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEAAwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSExBhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElKU1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD3+iiuTuk1LU/GN5YQ63fadbW1hbTKlrHAdzySTqxYyxuekS4xjvQB1lFc/wD8I9qn/Q565/35sv8A5Ho/4R7VP+hz1z/vzZf/ACPQB0FFc/8A8I9qn/Q565/35sv/AJHo/wCEe1T/AKHPXP8AvzZf/I9AHQUVz/8Awj2qf9Dnrn/fmy/+R6P+Ee1T/oc9c/782X/yPQB0FFc//wAI9qn/AEOeuf8Afmy/+R6P+Ee1T/oc9c/782X/AMj0AdBRXP8A/CPap/0Oeuf9+bL/AOR6P+Ee1T/oc9c/782X/wAj0AdBRXP/APCPap/0Oeuf9+bL/wCR6P8AhHtU/wChz1z/AL82X/yPQBv005rA/wCEe1P/AKHLXP8AvzZf/I9Y2pumkki88d6ujKwVkWKzd1JGRlVtyRx3x6etFrkylGKvJ2O3P+eKbtz/AHa8buPFmrXT+Toet67MxA/e3EVkoU55yog5GO+4c/Tm9HqHiryx5vii+D99tta4H/kGs69WFDSo7Pt1OeWLpJ23+R6zmivIJ9Y8RWzAy+L7lV9DbWoJ/wDINVG8TeIgy48VXWz+L/R7Ukf+QaUZqSuvyZnLMaEd7/ce1bhRuFeKt4n1/wAwbPFF7t4zutrQH/0TQPFOvg8+Jb0j2htB/wC0a017Mn+1MN3f3HtWaM14uPE2vt/zM9+Mn/n2tD/7Rpl14l8SQg+V4punK9R9ntsnPt5H+fakpXly2D+08P3f3HtOD7flS7fpXia+KvEc8G6DxPfq/cSW1pgflDVRvEPjoy/L4tkVAPutY2xbP12f0qvXQf8AaWGvZv8AA90A/wB38qkH1rwuPxN40VNs3iu4aTk7ls7ZR7DHlH88/hVK6+I/iDSZm/tPW9R8uMYc28NqyliMgbxAdp6cMoOO3emaU8XSn8Gvoj6DorxLw38S9N1iMrqHjDX9IuFHzrPFYtHnsFcW/PHqB0Neh22lXV7F5tt441eaPON8aWLDP1Fv7ikdKlFu1zq6K5//AIR7VP8Aoc9c/wC/Nl/8j0f8I9qn/Q565/35sv8A5HoKOgorn/8AhHtU/wChz1z/AL82X/yPR/wj2qf9Dnrn/fmy/wDkegDoKK5//hHtU/6HPXP+/Nl/8j0f8I9qn/Q565/35sv/AJHoA6Ciuf8A+Ee1T/oc9c/782X/AMj0f8I9qn/Q565/35sv/kegDoKK5/8A4R7VP+hz1z/vzZf/ACPR/wAI9qn/AEOeuf8Afmy/+R6AOgorn/8AhHtU/wChz1z/AL82X/yPR/wj2qf9Dnrn/fmy/wDkegDoKK5//hHtU/6HPXP+/Nl/8j0nha6u7nRWa9unuZ4ry7tzM6qrOsVzJGpIQBc7UGcAUAdDRRRQAVz9n/yUPWf+wVYf+jbuugrn7P8A5KHrP/YKsP8A0bd0AdBRRRQAUUUUAFFFFABRRRQA0fSig/WuO8SeOtL0JXhhdby+VgpgjfAXkg7nwQpGDxyc445zVRi5OyJlKMVeTOolmjt4GlmdUjRSzOxwFA5JJ7CuUvfHmmx3EkWnW82oSpx5iDbFnOCN556DOQCDxz6cBqXi3W9fjEU7qkR6xxIQh6Hoc55AIznFRR28lrbCa/uSCxP33ycf7I5p1FGkveer2R5lfMd1S3Nm/wBc1a+jeKe7SOFvvRQkRpjGMHksQeeCcVyt2bYkB7uNtp5+VsfQDFQXmoJISsC7AcZY/eP+FLb6Hf3SiRbSQoTgEjBP51rCMaS560uX1sjz4xlOXPVu36tmxp2v6daQpDtkJ7usYAP61vC5t32lSGY4wO4z0yO1VdP8PWtmql082TgncBgH2rTMaRgDaxDtgAvnn8T7V83jcThHU5qKbb6/8OdUYSjHRJHJ6rboNSfdFK7McgiUY/LbUaRwqoBgwOvLnmtmWICV/P3f7BPOVzmsycKoZJFL7j27V6+HxTnCMeyR5lWcr/8AAIZLrTrbaJtnmMPlijLNI30Ucn8KaZNSuiTBYW1jCehut0jt/wABVgB+Jz7VnpottbiSfT7+5i3n98quGdj67mBb9ccVknXNRjuiDfyNbWkisnmRAGZN+19xHZfXHpVuFSWsX951U6akv3dm/NfodL9j1DJ2yafk9T9llH/tWmSaZeF1IWwkGPm3LImD7Yc1oQX6TcEbecZzxVvHFYyq1qbs9DnjWkndW+5HPjTbaDyI7qykvtQnLbUtZnQBQPvfM4Cgcc56mqk9mumSx6tJbS3N8mFMYn3iJCDkklcnAJ4GBnntmtW/WeGVLu1JN3b5MYPAkU/fQ/UY+hArClnuNbaS3tEMdq3+vkyVYhs7owCvXHX0zWlOEpNyk9DqjXlKKenn6dv63NVdRSaJJoSwBAdCpBBB5GeM1Bb2NnLeG6mv7wXWMB1VVwPTAHK+zZFPFmttBHCqFY0UKN2S2BxmhWSPaUcp7OtdipQcdDljWcbuC/D/ACLN34b0TUUkVIEt5J49jNHGFJGQcjrzkDrmoLfQtQ0SEf2BrU8M0LM0KSOPL56qRtGO/Ix1qcSSMFywwT95eauRSpK3lswWU8B9uN319/euaVNo0hjKj06/h/X9aHTaH8Qrt5Ps+pwESDl87cr05DJlccHAOD7iu50/WrDUm2WtwjTbd5h3DeBxzj05AyOM968caMh9kwDqc9f5f/rqjeaP5ebu1u7qCaBTJEN3mqjDkbVbJU8cbCprP4Xvp/XU7MNmWvLO9vv/AMn+Z9FDgUleK+A/iV4hltY4db093t1UBLhmIc42jHzctxnk5yerV61p+p2WoxGS2uEkx94D7y9RyDyOh61pdbHrRrU5ScE9UadFFFM2CiiigAooooAKKKKACue8Hf8AIFuf+wrqX/pbNXQ1z3g7/kC3P/YV1L/0tmoA6GiiigArn7P/AJKHrP8A2CrD/wBG3ddBXP2f/JQ9Z/7BVh/6Nu6AOgooooAKKKKACiiigCIZ2jAFVNQu7exs5bq6kSKGNdzu3Yf1+neob/VLXTYg88mN2diAZLEDPA/r06Vw2qXlxq85e6ZvKyDHArHaPTjueTzXPXxVOj8W/Y5q+IVNWjq/63K+veNb3UZprLSPMhiDAK6ZWWTHfP8AAvt94+2cVx0OkrDI29lLK23k8Zyen/163ppDah444Hct93yhwPqe1ZstldTOZrp1QFvu7gMd/Xn/AArWjXnJO3uxf3s8GvOtUu5N/oQXN5BZhktyJJMcykcD6CufvL2a6fczs2eMnvWxd2O7zG+0Iiq2wghzz6ZA607SdHgJErXGSvKkJxnt1/Ou6k6NCHtHq+/9foXQUKceeS17FWysI7X57lGaU4wv936iumsddeFTFNuQIhKtncVHqc/4+1Yl7DDazkmaeRmyBgAdc55z9aaI7aKDyhHM8kuHYeYAcdh938fyrGvh6WLinW1X9bChKfN7WU/67HW2ev2Ny/zmQHsNnX8s1tSPFgbHXPpXmyXkcEwU28iZ/vSD/wCJFX28RRiQyNjPdtxBx/IfhXg4zh7nqKdFtL7zqpY2pFNTs79Vc7CaKGfHnlRg5G7iq7JbQlhCq7iPmKtyP61ylv4zsA7fLOsQyFdI5JEyO2RwT/hTotUu9buTbR2sdrp7Z3u4KyyD/ZAI2/U8+wrCOVVqd1KT5V93/BCc4vVpJ92F/cQWU8lpAVe7dQZJHxstwc4Lf3j1wuMnHYc1xOsG0g1OG1hmMduIUh3cFm3FwTgg5BIXJ5/XjsdSs7XR7RPs9pG4Z9iQqpLyOfctyepJPpXKnU9808gFnDdCTZ5TAyEAHAXaG5Oc/d7nFfQYSEXHnU7vr/wwYaMb+7ol/Vy3YGSzSOFZZHCALliCSB611tjP9ot93O5eCCK5TSbjUbyeRLjS4kbcNpEbLx3yNx/pXX24FvhPIQOVBO3OOPx96vHOEopdTkxMYKd7q/kVtQby1R8ZGeahtpGklLEfKo+XA/nWlcJHNCyGND3XOev51nxZcFPIjBPGPn/U54rCjKMqVuqOSUI33/MXylniG5yxQ4wO3PH/AOunw6aBy4J+vWrkEUMf/LJB9Cf6mraFemw4+tYVcS43jB6FxpJ7szH05CCQCjn8QT61UIaI+XLwD0J7VutLFBEzzthQe5H+Fc9eapBc5XMijPYKea1wkqtVtbruOVC+zLrqzgDOJgvGeQw/x/nUSXBDBXHXoQc1QjkC9Lq445AMI4/ENmrhkSaEzLPGw4Em5Dw3rgA9a65UHaz1CVCTV0036l1MnDI2R19jV23ubm0kE1rM0UwBXcADkHsQeD+PpWFDJ5bHypoWH90yY/8AQsVrRBpFymG/3SG/lXDXozpu6FGM4vRM7rQvGcN/KllfqtpesQqDBKSnHVT/AA8g8H1ABJrr+K8VkiD8SKQexxyK6nw74qnt5FtdVmLxMQsdweCp6AMe49zznrxyNIVlLR6M9fC5g2/Z1vv/AMz0WioIpkmiSWN1eNwGVlOQwPQg+lT1qeuFFFFABRRRQAVz3g7/AJAtz/2FdS/9LZq6Gue8Hf8AIFuf+wrqX/pbNQB0NFFFABXP2f8AyUPWf+wVYf8Ao27roK5+z/5KHrP/AGCrD/0bd0AdBRRRQAUUUUAR/KQKzNW1a20m1E0xyzcIgPLH/D1NM1jWINH003cr7gvCoCAXY9AM/n9Aa83u7+XU7tr6/lO5+FjH8I7AD/Pf1NY1qnJHTc4cXjFR92PxP8PUvXN9Pe3onkbdv6lhgKOwUen+fWq16kjW7CN5A7f3SFx+NU7jWY7WLKoAx+6D1rFfULu/zI7sIc42g4z3xXBRwuJqz9q/dS79f69TyamIXK9b3+RspLb2aNlgzL9+QtxnnIz3Pt+tYtzdG8klkBEar1568jp+AH5U0wyXBW3RdwGPov8A9etvStDjs7lJJpDI6jeoxwPWuuTo4NOrOV5v+vl+hhBOtoV7ixumt4QIkSF4wSVwD2JPr2HrnvU6WJhtl2jge361uTSRk7nxkcAHtWPq96JbOWEBhvUrlWwQD79q8yhi8RX5YKNtdf66m1ajSp3bkcpdXsFxqKwKxkIDM0cf3to5OfTPA5xyRWNZ67LeXFzN5BYq3yi3BlB69eMADHfGe3rRIJdMs2s4kgklmYjEcYYvCv8AfGVAyeDzyVOOtdPomltfaCIJ41s+G3CA7Rzkcehx6f8A6/pKlZUKak9tjocIU6aja9/6/Iz2vzdaVcTyzFowpkZnQ7lKttH3sYJPHp+VKuiTiCG81OKEW5wUty295DtO3dg4A6cc8gHPGK0PBiMJpbKQma1td0VvvVQeuHLAADHcZyeTmtDTLaCXQZLdcSQR3UsNrkZyquwUcdhz+A968ytj5U6vsbWWl7ed/u2X4jUOSnJwZnSW8z6dEIYF+zpwjDj8P0/lU8EYj2BUdJeM89Cas/2toulWXkR31vNLEMMkT+Zhu+QucVQ0631HWmlubKSCG1Q/62WJm83PUqARwDwPU59KaxsXTlKWkU9PO/l/W5wywtRtRJNZt5dVn0/TGD7Jpy7eUTvVFVssD2GSo+jetVtF0W0NuUQLLHDdPHbzfeBX1BPUAllz321r6jpFvpuizhZ5JtT1BRZCSVsufMIUkAdlBLYAA+Umt6K1iiCRQgKEXaABgAdOleb9f5Y3p6LX/hzulhmqai3v/X/A9DGNrcW2Yo5EcEEhMbSPf5e1SixZQ8k7bmJBBH8P0rbEIVMDGe/vUEyEAApkVzxzGU3ZfN9zKWEUVdmQ8e3jr70zbg1bmTapdiFXuScCse51LJCWo3MTjcR/KvTw6nV0icU42ZeyqAFiB9ahfU4Y22g8EcEckn2H+NZqL5qyTXlwX24AA45PoKiifbdmdlBx90EZx/8AXrujg4u/NrYnmsWb0m6TDsyLydrfez7jtWZ9mRdvPXmtEqZUdzkMScZOSRxThYM7RkRNk9Rzj8a6adSNKPLexMaktkZhZkmTYTjHHan2rSLNuXaIWO11xjcK1pdJOwZQZHbIAqA2eAFXlh6jGPpVLFUpL3WX7Rx6amZeWstsx3qSuch16Edj7UtjemGQBmJTse4rRV3MMkJBLQ8geq5559uv51m3FoFfzEBKNyw9PpW6kprkmbc8JLll12OpjmZkG1zg+h4qVZA5KsARj+Lj+VYuju8cpt5GyrDMZb+Vbnl8gdPU14OJXsp8tyYyl0d0amh+IxoLlJ2kbT2OSMbjESeWXHUdyv4jnIPosMsNxAs0LrJE6hldTkMDyCCOoryNoskqQOR09aXwf4qk8P6u2iXT79PeXbGWYDyCT1BP8OTyPx9c9mFl7aLS3X5Hp4LFuP7uex7PRRRWh7AUUUUAFc94O/5Atz/2FdS/9LZq6Gue8Hf8gW5/7Cupf+ls1AHQ0UUUAFc/Z/8AJQ9Z/wCwVYf+jbuugrn7P/koes/9gqw/9G3dAHQUUUUAMyKq3N3Da2zzzsEiQZZj2qzwOa8x8c+IBMI7S2dWjTMm8f3hkZB9OcD8eoqXKzUVuznxNdUabl16HP6/rtxrOptN5h2I2IIyMBFz3GTyeM+/0GIIozsZ2cNIPmwwBIH+TUtjpu2MTSJ+8Y8LycexrUS2UXK7ShUH94DyT6Vy4nGU4+7Hp+h837OpVlzS3Zz66fLduC6MQTxx/n0/nVn7A0fyldu04Jx3rrAiopyPlUcD0rPeMXUjO5ZYvRe9cEM4nVb0tFG1TBciSvdjNPhjhh3sFLseOefapxcyy/MFZF9GxmqF39itFErysjZ4P3mPB46H19Kn0+/S6jG2CVR/eYVz1aMqiddRcvVbehtTukobLy/4Yr6pfi0h3HGemK5Oa5l1W4WJi8QlOwOD0z3wK6PxAZWIjhtmk45wvH51JoFsHt1kubKFWiP7twBk/X8zXsYRxwuFVZw1fpf7mRGlKdWzXoZOleHlmtmttxVrOU2hAXk4UHe3+9nd/wACHrW29oMxw28p2kKHYDhhk8k9yTjv61Lo8ck+p6ncyCOEs6wGFH3kbN2HbsGZWQ45+ULWg9ojTJnJ+bdjHX/PH5V5dfMK0p2qPbZfL9DrrYVyfMluYkPhK0tLZ/MmvLxSGb7OZAikkliPlC7uSfvE9awNUl1S70fT7nTmtI7GRzBHpcBZpGO0nDEFMEAEMmQAAck9K9HUYGNp/KsnTtAg069e6aS4uG3SGIS7cReY258YAySe5ycACuCljqnM5VHd/LXfTfT1O+FJRWxxlnoXirU4UEt/PawAfutuIFjABGPJVQcZxg7hwORziutTR7q9YyardsT/AA21nLJDEnuSpDOfrx6DudwknoppoP8Asmonja89kl6W/wAxOJnWegabYy+fFar9oxjznYySY9N7EnH41oHZGueAO5pd3baaazhQWbAAGSSelcjVWb953+a/zE0lqKelZeqataacAJm3SsOEXk/U+lZ+qeIUA+z6dcwo4OGkKFgPpgEGubWyUyl5buN3JyxO8kn3O2vfy7I5S/eYnRdur9e39bHHXrpKy/Qs3V/c6m+ShCZ+VM8f/XNBCWYGdrysMBF5H4n/AD6UOiCEhbuEMQFG0Pj8flpttYKo/wCPmJmJ9H/+Jr6VQhGPKtIroeS4t3k/z/4JGGAUAqCx5Zj61btreWfCog5PU/zq9aaIJCGllUoSOFBH8x1/lWslksa+XAuwZO5gece1cOJx8Ie7DV/h+Y44ac9WijZ2aJKyyFWcNhRz8o9+2a1RAsSZA/GnwWyxJhRz0zUyxk9cmvncTiJ1J3vod9HCuMdtSqIVY52gnrz2qs1moZyFHJyM81rhQq4A/SmsoZe341jDE1YvROxrLCXRhC1eCVZVkVd3DRkD5vbNV/7OiAcQMEVgeckj/wCt1rdmgDLx29DVVAwLRshGeMnk+x/z616lHF1JQb6+n9XOSeG5fda0MKG2e0niW5DEhxjHUZ6Hj8j9fpnozDyOaRYY5V2uAXXgNjp3/oPyqwEJA7fWuXF411Wm9Gvu+RVGha5VlhBU8DniuP1y3NvqaTABRKm4buxHUf59a7hkY8HoKxfEWnmfTpXGSyDcoHt/k11ZVi1TrpSej0+8c4Wd7HT+BPERmtY7G7kJXhbd37H+4T+WPxHoK9BO3FfO/hvU5LO9VQSFbBDA42kd/bHr14r27w/qZ1O0+ck3EWN/y7c5zg/jjn3B4HFezWvTrOnLrqv68j0MFVa/cyeq29DfooopnoBXPeDv+QLc/wDYV1L/ANLZq6Gue8Hf8gW5/wCwrqX/AKWzUAdDRRRQAVz9n/yUPWf+wVYf+jbuugrn7P8A5KHrP/YKsP8A0bd0AdBRRUTusSF3YKqjLMxwAPU0AYWu6kLC38uNh58gI6kFFwctx+Q6evOK8xFt9r1fzZ5ssPmCBTgL/CPy5+pNauqz3XiHWJZEdo7VTzn2PC4Oe3PHGTmsrWtwK6dp4JmdsSuvJBPqfU15/t1OvyRerT/7dXd+v37Hh4qq6z5re6n95pedBNEwtm89lO3CdMn3ptvYG3mWWZIw+Rg5J/yabZi20KwSAFnmIDP/ABYOKt2BfUTHPME2Dlcdf8//AFq8urWnSjNw0p66tav0Xb9BKMZzUftdlsvXcuzxs0DYJXI5zWHqmpR6davNI52IMAL1Jro7hsQnA6cCuWudMh1AxrMQ437io43DpzjHPX6Vz5TVUnes/dT6HRjfdnFJ+pyHna3r8gdJJljdsRqrEZH+c812mj6QNDsnlu7h2kYAsHfIU+g/z2rTt4ILCGPaqRxouAM4Cisi91P+1HmttmI0xgNjkmvVrY6tj37GlHlpLd9bXM61fkh73yX+ZctbSK4jkdiWaX75B4P0+metLqV1LapDY6bEn2mbJBZT5cSDG52xjPUADOST6ZIq3WtwaasVjaJ9rvGT91Crcntuc/wrnqT+p4rIuNWutOmjtYEjvb15kW6nk+VWd/m2Jj0XJ9FUDqTXJyVak3JrTon935f5DoxjGF1v/X9eZ0WjW39nW3k+b5krs0kspGPMdjknHb2HYADtWnbXEU0gbcpJBC49OD/UVXW3RUYzgEMMFeoAx3p0MDbGdTnLEg5yT2/pXl11SmpSb17mlOVVNGgepxSYzUEMr52uPxqC91eCyBBO9+yKeTXDDDVZz9nTV35HY69NR55OyLuKNveuYm8SO4cExxDOMAkt9emP/wBVZ76ndPJtMskm5s4ViPfHB/D1/GvWpZDiZazaicM8zpR+FNnT3mqWloMb/OlzgRQ/M59eP8a5rVtTuLoeRLtiQfegiYkg4/jbjP0H40kl3dCCKOFUtV+6FQYOMdcnn/PPrWNPcKsmDLuBIB284Ne9luVU6UuZq7Xz/LRfi/M4a2MnW92AvkhDhee2R3prJl9pOfRf61ozxRxWaeU++WRQVB4AHqc1c8O6RMxNy4b5h95hhf8A69enVxsKVJ1ZvRfictOE5ystyvb6YZAnyjI5bJrZtdIRSGddxzldxrbjs4Yhll5PQf5+tORdzHgBRwMCvk8Tnc6t1DRHp08t5WnPVldIMAEgcd/Sq8LvPfNEvyxR8/WtOUERnaB+NEUKxJhVGTycdzXlrF2hJtXb0Xl5nY8NeSS2W/8AkNEfbrT8ADpQq/OTkUsinjHSuVzu7XOhRsr2GkUm2lJZfvUoIbpRdpXCybsRlKp3du5dHTgLnPJ6H/Oa0cUmK3w+KlRmpIznRUlYzYzhApYlgcD1/GnGTbjOOTirWwLkFckcZxUVxFlCcYPX8a7niL1LNs5pU5wWj2BQCx28Z6+9ZeuyrDHGrhipPUHHqPT0rVgcbSnAbHFc74muD9kMeBksCfauvLVOeMjF7f1rsY4mS9j6nNwrbrfRiFriNi3yjYGHOeMgjA616Lpd4bO4truNty4w+0Y3qfvccfX64rgbPTDfXsZ6RhSWIIz97HSu5tkYWoRsllGMkda9nOsRGEqbT95EUJ88k12/E9LhljnhSaM5R1DKcYyDyKnxXHeFtRKztYSMu3BaLoO+WX365H412FdNGqqsFNHs0qiqR5kOrnvB3/IFuf8AsK6l/wCls1dDXPeDv+QLc/8AYV1L/wBLZq1NDoaKKKACufs/+Sh6z/2CrD/0bd10Fc/Z/wDJQ9Z/7BVh/wCjbugDdwK4zxTrcq6smh2yoyvbNPdyEkNGpYCNRj+9tkznsvvXZ14tF4gj1KXVNchikla9vnhgVlI3Rxjy49vscZPoWb0rjxTmqXLT+KWi/r0uYYmahTbZs3E4RhbQsFlK7i/ZB6msmfU7PTruHTLeVUuZFZgJOSxAJJOPxJ9ccd6o6xqc2nxQ6baTN/aVz+8mliUFo4wCS3PAyRhcg4HPas2206LT/EGllxJLJMkjF2kLeZIFwOSTgBTJwMd8AZxXl4fCJQcpfDr/ANvNdfTt955lSV3q7Pp5f117nQCwSaVZrsB0cYWMjoepP/6624ZPkIVDgcBRQ0irsHG/HT0qaBEUHH3q8zF4p1YpzW23b+vM0w9HklaLGx+bJGwnCn5uB1x+NY82t2kWqnTba4ikvlAZ4S6g4OOBkjk5zx6Vb1SS6u7yPSrKUwsYxNc3C4zFESQAv+0xDAHsAT1xXKeIIBJprWWiW0cKWl2AszNmSe4KldqDBLNlgS7EYK55wTTwlNOV5rR9Oi83/ktet9DrlQvZX179Tb1DU7O8ZLeGQ3dwmd1rZfvCT23N0Qe7EVVj0LWLlmaUQafCxx5cf7+XA/uucKM84yp7Z9ursrWKzsoYIYUhSNAojTovHSpth2nB59TzWUc1lSXJR0X/AAf67ilgqcnzSV2clY6TZ2mqNDZvuihTdO2S8ryHIyzdzx+Az0FGlabeHU5rm+hHlwNKlmEPO1nLlmH97G1cegP96uljt4bfz3H35W3SOx59Bn0AHA//AF02DUNOviBaX1rcE5I8qZWzgjPQ9iR+ddFXMpShom1or+i1endt/JChhGotp2uymukxreC9gaRXcYZGJZcHk8Z4Pv8AXilYSQTZa4aMBioVQMEdc45/HOK0wFjjZpZBgDJYnAAxz3+tZGq3EVhFJKLmQmUbgBzwCOhxwO3PHPQ1z4evVxNT2b16LS/yvbou5nXpRpQ51p8/66lfXNaS1iMKu3mMOGRc/r/+uuUefUJNyyuyFsFmJO4jtz/n+VNS4zcNJJGDIxwqsOpPTHv79ac891cvHCkBZ3yo2r1Oea+1weBhg4KEUu7bt/VjyKtWpVne2pV2xRybXIJz6muistIeeNPvshVmIQ8dOASccnP6fXEUXhqaG4geSHeikmX5s4OOOx4J9uBXWW+lWoYyrB5TYxhGYcYweOP8gGuDNM4pU4L2Ur+as/1+ZtRwcq0tfx0OYTwxNcSuBO0mMgZyI19s9Se/AxxVyDwRahgbmeSUjptAUD8Oa6pVWNFVQcAYA6mj94TwFUep5NfOVeIMdPSE+VfK/wDXp+J61PAUob6szk0LT1XYYABnONxy2OmfX6VpxxqkYVFAUdBVeXzIwSsm9vQnAz+RxUJluAivLs3DH7tD8ufXJrz5+3xKXPUv6t/k/wAzSLpUG7Qt6WLvl85POeKUoB04pIZVlTcvIpxrhfMnZnYlFq6G8MSMg4puNw4PHTIqKWJ1cyRltzDFVPLuluRKjggjB4yB/wDWrrp0IzV1NL17nLUrSi7OJejGAQfzoZto4oRt/J4Pp0prxbmyT06VnZc/vFa8vuggLIc8g0uwg+tEJJUjkc9PSn7TiiUmpNBGKcUxhJUE9aiVnY8jHvVgRnHzc0xlJ6cGqhOOwpwluM3AHJ6EU15FUDjr604gjG4d/wAKk2KR0zW05RXLJ/1YnllJaFAxghnK9+Koapp/2yAICFUegzz/AJxW24UkLjA9qiePcxVuFAzkcV14fGSpzU1ujjq4ZNNHJXFtHaRukYdrmSFvKaM5GB1A468elTeDr9poJbSSVpGjO5Gbup/w/rXQXFhHPFEJB5pjzgt3yMZI6Vyg8PXei3h1DRVWZkd5ZLVzt80NjhXP3SMH2Oa9pY6hisJOnUfvdG+6vb0X+ZhToSozSS/rqdYGltZ0ngOJI2Dryccdj7HofY139hdRXtpFcRnKyKG6jI9jjuOhrzmzvrbUbNLq1lSWCYbkkQ8EdPzByCK6rwtdFkmsnI/dHfHyM4bOePY85/2qMprSjJ0pnZQlyVOXo/zOsrnvB3/IFuf+wrqX/pbNXQ1z3g7/AJAtz/2FdS/9LZq947joaKKKACufs/8Akoes/wDYKsP/AEbd10Fc/Z/8lD1n/sFWH/o27oAj8Xa4mheGL695MqxMIVBILSEHaAfWvHWuU0fTYXmbzXtkVIlzzPKw3E9O+SSc8cntXY/FO7uZm0zSrdZGiMwnnIQFQFyy5PblfbqOtcJeW015ZXN3LLbpkmOGJZQRHlsOxOTyxTj0GR61m6PPNX66f5nlYu9SooP4Uzm7HU7iXULqWVy8su55ZBxuYgj8AOgHYYFd3d3dva6lolvuH7iSR2bdnKiNkP6sufSuc8OeHwb2FvtELureayqWyQOn8PritPW9PtZtWa8u7thFIgQRRj5WjHT8CfQgHvxW2NpUqtSNCOyXT7l+Fx4pwVW8PT1Oy024g1GCO8t5VmikB2unI9DzVyRD5RRRkkc9v1rF8JiJPD9qYWPly7nQFQv3mZ+gzjrW5OjPDtztGQc56YOf6V8RiIuNflurJmkKa5HYxr+4+ya5aTqwEU8Zs5pWztDE7oiT9d6/VwOM1etNA020v2vo7fNyzM+9nZgrN95lUnCk9yAM0XkENzpN5HII5IpIGVgnGQV7Ht7GnafdebpVndSOSZYI35G5jlQe3X8qVW7pcsHbo99fwN4aLnkzTXuc9/Wl9u1Zdxeuqq0cbnBBHUD8cdR+f4dahj1mBvkM0ZkIOQj/ACqPc446isFllWUeZa/exfXKSdmxnid86eNNhx51+6xYGRuQsqycgcYRjzVXTIYpvFjT28EG+1smtbqaJBgvvUhNwABICtkdsr0rL8R6+2kpBe77cD7UqFrgs2I3HzgAHqAAe468d6l03xTBpcUtrJAiwwtLLJ5UjLNtLbvMMMqq2DvBJyTknANepLBVKOHUI6tp9+r1/Cyt/wAAujNVffT0O0EWGOWJBOQMcD2FQ/2fZL0tIQMFfuDoeo/WlmlaNcjzHYjhVHJOCcdMDp1NU/t7vtUoYjIDtMrDhh2OOOgzwfbivHpUK0leE7Lyuv1X5BUnRjpJfgXRDEHBWMZHQ46f5wPyqRYlI6cZzzUcI8yJSXDdwQQQfQ1Ns9Gb86wqOzs5/ma04RavYfgClwKYEzyWb86Co9W/76rBxh/N+H/BN0O+lQzx7+ZJTGg6MrY6jHPam3U8VrEWYvnBOAecAZP1rFzqGpRs2Ej8xAVSRS+1Tn1xzxnnHbiu7CYPn/eOfLHu1+X/AAxy4ivCPuWu+3+ZqFIhOsqBZs8BwM7VwO/PJ49OvtUlxEJItvmBUwdwYZ7jnP4GsWLXdPgvpbMXM81xGWHz5UcdQMD+lbMNzFMAC22QjOwvzj1rfEYarRcZu9uja6dHb9dDKEqUrwbs30/4I8bbdUXIC5x+PQCrOeORz6ViX8DOqtBjchGxcZ/AD6VqxjKjLkn61y4jD01CNTnu3e+n/BNaFW8nC2iJcVGUYsD296dsx/E350u3/ab865UoraX4M6XFSDGe1NMec5J5/SlK/wC22frSYP8Afb9KEoraX4MTS6iogRcEk+5pQfbFN7cyN+n+FO2nH32/T/Chxi95fgxpLoJuPTFRb9zEA4YHBBFPKsBwx+nFIFyeWGfcCrjGC1uvx/yM5JvS4rLlT7ihVJQFeM80uG/vD8qau8LgMvBI6f8A16pRvTa5lo/P07FcqB/kXJNRvEs0ZjlUMjLgj1FSHcepWkMhXrt/OnGMl8LV/X/gGcknvsKAEj56AdTVdF+dj2PSpZNzREFRj61VhnVpnQEllxkY6VvRozcZSWvzRjVaUoor6Zo9po1vNb2akRSTPNtPOCxzgegHQCtGwuPsOpwzk4Qtscb9o2njJ9hnP4VRn1CKK4WEENIxOVBGQAMkn9KcJlmiDAjPQgHPPcV30oYmE41pJ9znnOLem6PTq5/wd/yBbn/sK6l/6WzVf0i7N9psExLFmTDFgBlhwTx7g1Q8Hf8AIEuf+wrqX/pbNX2EWpJNHemmro6KiiiqGFc/Z/8AJQ9Z/wCwVYf+jbuugrlHuorHxnr93cNsgg0eylkfBO1VkvCTgcngUAeR+M79bnx7qFw8ZEVikmUJLbmXCL+ewn8ageI2WlWtpH81wkSozbtwBwd2Pqcn8a5fX9Rb+07gPM7yzmFWkLEljlmYk98nA/Gum0u+hhe2mnwyxr8oAHLZwK6JU2vfWvL076Hg4iEmlU6N/wDB/U6TQNM+wWgeRP3znHqQOwP5/hRqMMMkA3qs6yP5hdS2MgYGDzgHII57+nIox+JFnuBDbl2LsGOw7S+BnGSOBkHJ44H1q94mmEOhXbIWecIyRsBx5jYQc9M7m/D0r5+qsRTxSqVt5f0vTUnl5o22fn5nPaPBe39vps1hepEbKCLyUdWaOWR42kaMgMvRTHzz0rr7XWGutGt7+bEXmruaMc7TyCoPchv5Gs3wWkC+GbUrAxuYXkLx7cMsmTlefYgDtjFYnihWsdVSKzS48m5SSVLcSFleTcpJUcheXYtzjv8AXnVOnicV7GatZvX5/rv952VJS5G4v0/W5p3HiMyWl3Y225rq4RoYFHzFpCpHbsOvTgAkmursbT7NZxwDhUVUXPUKAABx6VxngvRF0zzNSu5hJdsSrOycKO4UnoM8cehzXZTX6RwBwcliFBwcZPT6CubNqaVf2WGjp1fmKhKnCPvSvYZqEBeFw8m1MDJJCgY9eOlcfE8GnXTm0vILhzwUI+X06nqMntXQamk11BBbvwzqNzkZGcY7jjk/57cbrS3aB9E02IfaViD3UyMWECHONx7sR0A9SelejlUVCjyVJXv00tba/wDVuhy1Y+1qNxVvP8S499p+u2f2GS+dZpLh1guUXIRlAU4z1BBK8+p9qkl8Q+GdV05rLWbZbTFsfJaG4VnZQdiqhBDbiCeCMYJ7VycUcWm6ZZW4B2KrNz1J3kZP5Ve0Kz0NdXgl1I2drHbv9pjkaUKWwBtQqOBghmyeW7cFgO7MsthCmq13p27+Xrbt5nqUY+xkobo6fwPrE2pGaLzLq7hWIbZJ0BAZZGAww+XlfLbAJwc0XOn6tcalbz3eWyS72yknYucHBJIPQcA9xxRpGrL4fuzv042OjX8xeFpyTNvMe4s3zNwdpAHBGBXZwrbTwxSxopWQbkO3kgndn8eD9a8GeMqYKq6safuy0T36aq+yfe2u2pNakq+iZU0m7ubtHae2MEa8KN4PIJyP8+v56qsGAIOQecis+9lS3VQysd2QCDkg4PTqc9egJqHQr+1uLMRQygunLKWJIz15PXnNeRicO61OWJhG0brRbffr+JdCrySVGT1L1/diys5Ll/uxjJH6dqp6frtpqDusUqghS4GD0zjJyB7fn+NR6tPbXtvNp7PIszAFljHzgcH/AA+vOOa4uG+sreaVrdLi3Yx4WRl3Lu4O4nqCefz/AD9PLcohicNLni1Pp2tpb9TDEYuUan7t3X9f8A9BuYI7rKyqTGQQTnHcHHqPX8Kwv+EiFtqElhdWnlQxkhJd2QVHGcY69OPf82W2vTWtk63IR5PK3o3OH5xg/hj86q6l9nvtOhmW38tguQSRnBHI/QVvg8tcJuliY3heyaez7r8jmrYyLXtKe/X/ACMC71K3g1F7iwthFIxPzZJ3Z7nPIP41e/4SC60exgZU8yedsqznd/L6nn8O1ZbpGifvUGWGck4HX/8AXVi006G9ubJWKFGdgpRuenOA2RxjPI/mK+srUMP7Ne0V4re+t0k+pzQqR54yd/8AP+vyN/TvE738fk3lu8Fy/KlVKjb68npwe9a1vrANz5MksRdVIwudxI9cAgZIb3GPfiHRhO0b2d9ao6RN+6YLnK5PJz9M5HrTLi307Rblp4IJEuZSFXAdgfwztwM+2K+Vrwws60qEadn0s019+6Vt97anVzTS9qpWXz/p+Rpz69Z2s6xTsVDHAbHGc+nX8anTVLOS2F0k4aLCk7PmI3dAQOc1yF3pk1ve3HnuslqxLRoWztPU9uB8xrFs9b+xJJbpbsQSOd2Twe2CPTjr1960hw9h8RSUqEm3pfazv67MqOOxHNKNr2PUlmjkkKK6lwAWXPIz0yOoqQgGuJn8aPtC20A8xgMDaWOfT+XrVT7fdTQbtQ1iRPM+VovJYfL+GMA+o968+PDeJ3qPlXzb+5HS8wgltd/d+Z3KXdq1w8KzIZExuGenapdzFieNuOPWsLTbayhHmWtrvBXetxsA5xjA3HP9K14C20lyCOCDjB6V5mLw0KUnyN2Xe1/uW3o9TWjXlNLm/AfLKqJl3CLnGTVe4vIoF3/PJuzgRjcfyqlfWxvkcxu058zCbJAFj9c/r6miS2v1VSksCz7NrYBIGTyRxnrjvziuilhaKUXKevVPTt5N+un+ZhUr1W3yx077/rYv2l7BcYEchJxkBgVJH0NMuNRt7Z2V5F3bgMZ6E+tUBpA+1i6uJvNnznLDCjjHAz+dSS6bHJE/mnMuS6PnG0+35CtY4fB8799uLX4379V8l6B7bEcjSirlkX8UhJ3EKuct7/5NTQoGO8fme9Zdh4eghjT7RI0pUfdPTOST9f8A9frWiRI8ibZcRgD5VGAT7+1ZYiOHTcMPLTu1p/wb/cFJ1X71Vf1/XzLWQyjHQ1haox02Ca6GCwQhMLklifl6+hIrdxtUDpgcCud11s7Z5GVYonwqtk72wDyM9OM/hTymPNiFH7L38/L5/wCY8c0qfM1qjn3MtnbK6iRnkIkZ2YgsDgAcdcnJIxnrTZ7vVbGOMLZ+VFGPMAxyQcj5sHg89OtQW8d1d3cyWTeaIlJQcgZwEDYJyCM7vwFa6y3X2SdJ1KPCVkuGOD5vHzAgjj5QAMY4Gc19tVahJJqMu6e+u2n4bdemp48YqKu9/md/4G1KK/sZo0keQKyyLnooYY2jng5Vsj3+taPhAD+xrj/sKal/6WTV538LNTA8UXlg0qqrW5KruwJGDLjAPUgE/rXovg//AJAtx/2FNS/9LJq1lS9k+Q97C83skpdP6R0dFFFSdI2vPfFt3FYXfi2WZS6todlCAADzJLdxjqemWGfbNehV5L8S7uO3uNdt5A4e7stMiiIAIDLPdy8+2IyPqR9adJXml5g9jxzWtM3atBLks87oEXHQq4X/ANm/SuktNEVii3UgjjU8oThn4zge2OSeg3CtKzs7WYrf3hWRbXJTeMkZBJx9fYfT1EN3JLc2dxdQkoJd25FxnGeRnsMHnpycY5rp9tJzlCGm2vn2+48BYqpKMYp26f18jR02bThqBFlAhYEOHI2/KDycn1OMfjget3VNNGq2w0xJZI492+N4MK25WVh2xxx+VYnhazmvXlklk8pUUqg75PGR9M5rY0O7uLE2s13dyzxzO1reCZUzb3OQByqg7WII5znch714mZ8tGq+RuUkuvf8Ar/Lc0p4eTqNJ7P1KGjarNBAEktz5moMHikUgL5nlqki57Mro/pkZIPBrrJSFs54UOwouxXB5JK7iemB68f8A1q5ye3DeEtZiIjC6fNdSwleCrq5mQkf8CH1GK6aR44FUyALIY/MPGPQE/Xn9a8nETjNrlWt/8mr/AC0fzNqkHG8tlt/w39dDNdYdOaKFJZWkXjY3I549PeuX1q8uZNXSHzDFEw5G7BOMjJP+NXJb22sLxJZLi5leSTzDhh1OPlxg56npgelbuoJposzfzQpJ5MbMCegU5yctgYHU+n417FOawdSM6kXLmVr26/P9Thpxu21trZfhcW01e2sdKj+3XMKNbQZkCv2/hwO5wB+dVE1dr+4ksLmxNjf+Ss+wkSDy2woY8deoIxxtPOOazNO06URLqUsJTLiW3gkB2qADtdx3f26Lnpnky6teReILqwsbeGTzY3H2i+ijAaHKn5Y26hiSp46A8+3FUw0YV3VpK6u9ez1e3bfS51QnGpF0pvW39f19xma1ocsdsgjGzZFjJTIzuLfzP9aq+HbW6sbK0upLnyo4LiN5tgXYqFgGdxjLcZ+ZjxuLDFeianp0Wp2Rs5ZmiaQ5Up1OOorEtPCOpWMZNrqu0glguwjJx0+9jsK6/wC1aGIoOnXmoy6XTt67f8MaxjVg7X5rNflqiO6u9Q13xGkENpDJY2U/mpKg38iDjJyAd3nLgZGNjc11lgJzZwG5CibywGCjv/nHHajTtMtNMgaK0gjhEjtK4RcbmbqauYGeOtfI4vGRnD2UF7q2fX8+p3xpe9zsy76yubqOeJGURyBhh24IKEY6ZHIB4Pc/SsG28OajYarJIs4eB4+ZAed45BKnrzXYsivjeobBBAIzyOhqEiY3Rw48kp028hvUH861wua16UJU4Ws1rden4mNXCwe99Tktdu00/UxKlrM8rRiQXSSZOR7Zxgjjt1NYB1Q3F9cQxTb43yE3DgjsMH6Cu0v7F1lPkku2w5gWVkB5wvygEYwcE8DjPHNc3eeFrw3xlseCxyWcABT36c9+wr6vK8XhPZJVGk7bvy9dvK2nY8ypS5ptSXp6/wBf8EkmhY2pZ5I0dQSsTgKx9lHqOM8d/wA00a2mmuGjvd0cEab8PwB2GM/WrdtZ3NmkYh0mKJojhZnlJdx69B1J6H8RitS1gTUYSsjo8hx5+x8AdSMY4Pb/ABPWpr45U6Ulf3X1Vm191/xd7uxhHDrn5Urvz6/1oUm0CCRrdzte3ZCoL8kk/dz6kdKSHwlJDdRzJJ0ILo2Np65wMfp0rWcQJYoIJzKYXDJtIbLDr169av2kzXLG4G4RsowrDkEZz/SvHrZpjKdNyjLTVarz007tP8DtpYalKXI93bZ/f+K/ER/stqF3SiNVbHJ4BI7/AJjr7Ulzai4VgyLyUIcEgnacjOPf37028gTyGaRNzFiFA569Ouawnu/Em8iOKFlVtpwAB6euT+lcGFw0qyVSnUSa3cnb7v1udNWrGDcJR0fb9STXbm4tLiERW/mxvky7xwvGAB9ef8mvPNWWaC9mWNMxYDlE5ABA6V0txY6vfyH7Y7GPJcsVCqPc4+nT2q0mnoIy3lb5PK3K7Ajk8Dofavs8A6WBpRi2pO2tvw+78jzfrHJW50rp9DlLQxMTEWbj5hvGDg1eW4uLKZooHSSMqA8cmB7jkc5yOuc/lVLWLueS9jyI1YoFd4wAeRyKs2XmNcgL5kvlESfuSVZefUcfj9K9ea5oc0102/q36F1Iuyq9Gr23/wAv0PQNKnt/syQxkPKAA5iJaIHrwT/n1JrC8UeI2Cm20+UMAP3hRc479e/4elXbHUYXto2heOOZQR9nZQGA4P5+4459q5K4mF34g+zpp8cC7OYw2OSeoYdevHavmcvy+DxkqtWLfLd6239PytdeZcarcOVWsld+hp6X4m1ppSjxwSA/MWIKgDgc7fw7V29vPNLAWuLdVY8YjbeCPqQP5VhWGjWYZwDJFIj7Q6OAW5PIx9QPzGOK6aBCqfMckdTgCvKzuthHJexppP0a/Wx0YV1JybeiKd7byM8U0KDMeSwzz0P3R0z7/Wq1pbTIkX2i689yA5DrxwO2Rkc4OT/Wtk4BxmqWbdb1/LjT7SR8xIxuwOMn8a87C4qbg6dtk+n6vZehtUoRU+a+/wDX6Dpp/s8bySMioqknPJHPHHeszUvElpYwxPEyzGTldvTbz/X/AD6v1uxuNSsmgjnaHeyZ54AGd3QZOc9M447Vy174QuLWIy28jS7pcmNR0QDI9+3rzxXpZZg8DV5XiZ2d9uj9WY1ak1fldl6f1+R18es2z6abxnxEGIy3Xr0Nc74n1E3Fgvy4RgGBHIB57/TP+TVT7Ld3ugG1mneNFdXdt2Q6jIIGB2wOOeg9a07bw/bpZvZSOHlQq5GcHHOAf15/pXoUcLg8BV9rJ3ak9N7LTX8TllWnVile/f5efn2MnRHmgsJZoWBfaCwD446YOPQAn8ajh1S6RLglwTOWO08gHsR2P/1qbGosriPzGAUNllJx6EHOKoaiYYj5C5IDHAPPH1r340YVKkm1fms727HHFupOy6s1fB80WmeOtMuCMpJIbchTk5kBVT9MsPyNev8Ag8f8SS5P/UU1L/0tmryHwpbxHU9OuJ32iK7Rjk7RtU53Z9AeteweD+dDuP8AsKal/wCls1GImpTst0e3gZ80XF9GdHRRRWB3CV4/8Ugf7Xv2AJ2W1g3Hb/kIV7AeleUfEQp/bWorIjOrWunLtVtpOXve+DV0naaZlWfLTbOAlRRaxoZFYlRJsIxhQCeoPPOP0rqbLRS3h+G1ePMsxDSuThlyck5I7cDFZVjp9te6rh4JMeplUjapyOAvciu2iYbAQRivOznGSpqEKb1vzP8ATb+tDx8PCM1Zv/hyhHpyxXCQqSkSqW453ZJyPw9fpWdfBU1Gd5okuLDUQLWdFGNrhWKHng7lIUnOc7K6YDeVb5TjocVQ1PSUurCSBG2yeas8e0hMSIQw5weMqM8Hgmvn6eKU5fvpf1fV7X7HZDDuHwPT/gaHFm4vNPuLiDVE36XqcBtJruFtzZCFRIRgYYpyw5+7ntg7UWqs2gWF/cgmVoBG2AcMe59cErx7UfZo9T8OPcwQXDTCPz4oklAkWdM/KD0DBwRyCOvB6VLfvFqugi6RnmiHJ5ERypwwbcBtIO7IwOn4V6cVRlXirq17Pp0t+P6GVRTlQtLc87kvJ4dSSRIHumZmwkbDKqvBPP1FdLaKL1Jr3UnupCJQUtpZjtjQN8o2g7ScAHnPP0qO8bTrDVYRFDI4miaJnSUKpcENj7vUgHn/AGfxqpPqdtG6h7abLHOPOHrx/BX0k4fW1e1v8jB89VpU0lp3Oug1O3vme2hKkBejA8Dpz/ntTLSzTRfExhEhisr5Q6bh8r3AzlQf4TtC8d8HHesjS5YbK53x2khjk2nzJZgAAef7tatpJptzPf67LZpcRwyloryRFcJHGgO6PIyF3BumcnJzjFeBjsP7BuNN+41972+/8exeGjBt8zvsb8usaba3a2Vzf28Nww3LHJIFJGTjr9D+RqW71axsbU3NxcwpCBktvyT6AAdSfQVxqpY3DTR6nPqUd9ciOWSxhQOZMxr8w2oCwGCCx4BGOOBTPD+gaLrb3L29vqUEET7GFxFGiz5yDtcAsV4IOGAOa8Wtl1FNynLRb6fl/TPXSsdpa6pDf6ZDqULlLOSHzcyRkNt659uM/pVLRPEa6rHtu7KWwuSgmWGU53RsTtYN0PAGR2NarWyvA8DBfKZdhTbxtxjFcZY282n215c6mbyG4tsWcRndHV40+4yEKDznJ47VhhcLRr3gmr3Vt76/LXzJqzdODkdo86tFI0UoYpwdql9p69Bz3HFRpau6qbid5GAI4ygIPqucE/54rltO8VCW9n07TdNkuZ4+oXCIrdeScYXnqAT7Gtm21q5+2wWeo6dJZyzsUjkDLJEzBS20MCDnCk8qOlGIwNbCtwVl16Xt91++xFKUayUpfdqassZlBiYfI3UqxU4+o5qO30+C1ZmhVlLYySxOcf8A6z+JJ61YKsSCSuR0+X/69BZv7y+/H/1689Tmo8kZpJ7rX/I29lC/NLcz9R0+wuIdl2WICg/fOcDjP6/WqUFrp9nchIYXMpzlmPzJuJwB29f8mtzaWBGVPY5X/wCvUb2ytIsjJGzoflYryOMcV2UMZKEPZzqO3a7t+RzVcKpS5opfcQROhjjRUclju2gAYH49hkdKnjVYIliG7Cjqcn9fwqRVYDhgfcj/AOvS4Y91/wC+f/r1y1JqWl1b1f8AkbwpW16gsiugYEMp6EVmXt6FmZDFuKjcijl2IyOBg9+9Wrq2a5jCGQoAc/LwTVKWXyIYnVWnD/IXUZI45J9+MfWuvB0KXNzJpvtqvxtbuc2JnO1nou//AANxj28l3ukd3RJBsliVgcA4Hp/nnHuy505bSykl+1MAm0ozj7gHBA+v+H1qfTFjaEiOFkVwfmLFw4POc59+/wCFWLpJpphbyQI1uyn96DyG5GMduM8/X8ep1p06ypqaUVutNl8nr03vbTUwjQjOk5vVv13f9djyrULcT6gdy7Lgkk4kBxg9x710lnpyWOmSLIgWS4ACqflUEcZye/PT69azr+xttN1BkdLgyK+Qyzg5/wDHa7eC2tNTsLYzKFLkSbG4YlePbOM/yr6zM8cqVGlK/uPqvJXXyOb2cqqVOMtl1+44JrG6Fi0qM4CgyYychFYgE/l0z2rV0/RrLVLEmz8wXYx++kOVPPJwOePw7V2pVWb93Iu4LjcvYf5H6VNCGA6jJ9RXjYjiCpKn7toyv3f4rl1R0wwd3ZyMrT9KNuD50McyjOxpEBYAjkE9+47VtdVwOtIocryV59v/AK9L8/qv5f8A16+ZxNeeInzzkr/P/I76NCNKPLHYQnsQapxQSI2+Qrv3NtxnGMnrz9Ks7m3Y3An0xTvmO35l6+nsfenRvTbSa19ez8huEZtPsVoreQRDzpEkuMEb9uByew/AflRFGYiscj5Zhgc4H4CpnEium1QwzyQOB156/wCc1k3+qXEFxCiWhYCUK7swTHQZA5yPm/p1rooUquIk4wktdd0rb+Wnoc9RU6dm+nqRQ2zQay6pEGxtYll7HI4PqOT+OO9VNT8+BlW5cKXR4/MCjLoBnr26j5fUH8NaygmEjiSSQ7sFGfBI5OR0GO3Xnmm6jpn9optncqyHepB4Bx06dK9SniYQxMfayTSSTet/VfPdHC8M5UnyaXZxMwEOrRzNK5USDEjLtVgB15xjqvBznrzmq2t232S6COoL7AWGPu5J44PHGDWx/Z0UuoS2N3FKk3zOrGQEOCAM528DgfrTJdPEVsxljnaIZiYNMrBcjIA+X6Y98V9XTxEFOL5un332af3/APAOVe5JN7r+uxd8OiNtKmRJ/OIbJOwjbxjHP0/zmvTvB/8AyBZ/+wrqP/pbNXnujS+baSbbdogCAcsDuPqffGCa9D8If8gS4/7Cuo/+ls1cFJt4irfuvyPUwFrNr+tzpKKKK6j0BvWuc1jwho+v3X2vULed5/LRCYryaIFVLlcqjqCR5j4JGfmNdLSUCaurM5FPhx4ajOUgvlOMfLqt0OP+/lP/AOFf+H/+eepf+De7/wDjtdTRU8vdIS5V0OZHgDQR0Gpf+De7/wDjtJ/wgGg5zjU8+v8Aa93/APHa6fNGan2UP5V9y/yKsjjh8NvDIyBaXmGYucaldcsTkn/W9SSTmlT4c+HIYWhS3vlibJZF1O6Ckk5JI83uSSa7CijkXZfcibrscZJ8LfCEybJbC6kTP3X1G5I65/56evNRN8I/BLnLaTOx9TqFyf8A2pXc5ptaXa2Bcq2RxjfCnwa6bH025ZMY2nUbkjH/AH8qVvhz4Za3a2a2v2gZSjRHVLoqVIwQR5mMY7V1+aM0WXUFyrocfJ8M/DEoIktb1wV2kNqd0cj0/wBb0qWH4e+H7eFYYU1GONFCqiatdqqgdAAJeBXW0lZ+zjtyr7l/kVY5j/hA9D9dV/8ABvef/Hail+Hnh64jMc0WoyIeqvqt2QfwMtdXRR7KN/hX3L/IGkzif+FS+C/OSYaXcebGco/9o3O5T7HzOKnPw48NMyMba/LI25CdVusqfUfveDya7EUVTim7sLHLt4B0F+q6mfrq93/8dqpbfDHwtY7vslpeQbsZ8rUrpM4+kldjS1KpRWnKvuX+QWOZ/wCED0P/AKin/g3vP/jtNbwBoLjDDU2Gc4Or3Z/9q11FFHso/wAq+5f5Acz/AMIHof8A1FP/AAb3n/x2l/4QTQ/XVP8Awb3f/wAdrpqKXsYfyr7l/kByx8AaCTkrqWfU6td//HagT4b+GozlLe+U425Gp3Q49P8AWdK66ijkilsvuRNovoc1/wAIDoQ7an/4Nrv/AOO0HwHoZ6/2p/4N7z/47XS5ozT9jH+Vfcv8ijkG+G3hdpfNa0vDIDnedUus59c+ZUifD3w9EoEcOooB2XVrof8AtWupAHpS59hVXUlqieTyRzI+H+gr0GpD6atd/wDx2j/hAdC/6if/AINrv/47XUZozU+yh/KvuX+RVkcwPAWhKMKNTAznjV7v/wCO0v8AwgWh/wDUT/8ABvd//Ha6XNGaPYx/lX3L/IDmf+EB0H01T/wb3n/x2j/hAdB9NU/8G95/8drp6KPZR/lX3L/IDmf+ED0P11T/AMHF3/8AHaafAGhHqNSP11e7/wDjtdRSZo9lH+Vfcv8AIGr7nMDwBoI5A1MH/sL3f/x2j/hX+g+mpf8Ag3u//jtdPmjNHso/yr7l/kFkcnL8OvDcxXzYNQkK5279VujjPXH7ykHw78NiPyxDqATAG0ardY9Onm11tFVaNrWJtHscnH8OfDUSlY4dQQHqF1W6H/tWtzTNOtdHtEsrJHWBWd/3krSMWdi7Es5LMSzEkknrWhmkp21uCaWyJKKKKZR//9k=",
"text/plain": [
"