Error loading file viewer.
Select countries/regions from HEALPix
The snippet can be accessed without any authentication.
Authored by
Lukas Kluft
This notebook shows how to read Natural Earth shapefiles and use the information to select country data from HEALPix.
healpix_polygon.ipynb 61.01 KiB
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "ba18b705-0366-4c78-a6e7-f113e771aaa1",
"metadata": {
"execution": {
"iopub.execute_input": "2025-03-06T15:48:04.348637Z",
"iopub.status.busy": "2025-03-06T15:48:04.348281Z",
"iopub.status.idle": "2025-03-06T15:48:12.102594Z",
"shell.execute_reply": "2025-03-06T15:48:12.102328Z",
"shell.execute_reply.started": "2025-03-06T15:48:04.348612Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.image.AxesImage at 0x143da8260>"
]
},
"execution_count": 1,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAUoAAAF9CAYAAABiViEQAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAArbVJREFUeJzsnXecHHX9/5+fmdm+e/1yl0tvhBB6kRIMVZo0AQU7qHxpokSwASKgiCBNmgoIKFKkKAJSpIjUBCGUUBJCSC+X63fbd2Y+vz9m+87M4U8hyd08Hw8euZt972c+e+y+9vP5vJuQUko8PDw8PBxRNvUEPDw8PDZ3PKH08PDwGAZPKD08PDyGwRNKDw8Pj2HwhNLDw8NjGDyh9PDw8BgGTyg9PDw8hkH7qIbpdJpsNvtxzsXDw8PjE8Xv9xMMBoe1+0hCmU6nmTJlChs2bPivJ+bh4eGxudDe3s7y5cuHFcuPJJTZbJYNGzawevVq6urq/icT9Khk4cKF7Lffftxzzz0ceuihm3o6Hh4jnsHBQSZMmEA2m/3fCGWBuro6Tyg/Jq655hq22morjjvuON555x3effddHn74Yc4++2x23nnnTT09D49Rjfgoud6Dg4PU19czMDDgCeXHwEsvvcScOXO4+eabeeKJJ7j//vsrHr/tttv49Kc/zbRp0zbRDD08Rh7/ia79RytKj4+H9evXA3DZZZexcuVKLr30Um6//XZaW1tJJBKcdNJJRKNR1q5d631ReXhsArzwoM2AY489lsMPP5wPPviAk08+mR/96EcsXryY559/noULF3LccceRyWRQFO9/l4fHpsBbUW4mXH311UyfPp0LL7yw4rppmjz44IPFVaWHh8cnj7dE2UyYPn06V199NfX19RXXFUXhwAMP5Oabb+aHP/whg4ODH2k8r8yoh8f/Dk8otwCuvvpqzjzzTC6//HJOOeWUYe3PPvtsWltb+d3vfvcJzM7DY+TjCeUWQF1dHcuXLwcgFAq52nZ3d/OXv/yFnp4eLrjggk9ieh4eIx7vjHILIJFI8MgjjxCLxbjuuusA+OEPf8hdd91FMBhk/PjxtLa20tnZyXPPPVd83tSpUzfVlD08RhTeinILYMaMGRxxxBHMmDGDcDhMf38/f/jDH1izZg1HH3007e3t9Pb20t3dXXzOnnvuyQ033LAJZ+3hMXLwVpRbCN/97nc58MAD2Xbbbent7SWTyfDCCy8wZ84cpJQ89dRTfPnLX2bHHXfk0UcfZezYsZt6yh4eIwZvRbmFcMABB/Dkk0+yxx57cMwxx/Daa68xZ84c7rrrLrbddlsOOuggpkyZwpNPPumJpIfH/xgvhXEL5oEHHuC4447jkEMO4Qc/+AH77rsvQohNPS0Pjy0CL4VxFJDJZPjGN77Bsccey3333fcfC+Tq1avRdZ0pU6Z8TDP08Bg5eFvvLZQXX3yRwcFBLrjgAleRXLFiBbNmzeInP/lJMQjdMAwmTpzIXnvt9UlN18Nji8YTyi2U+fPnAzB79mxXu1NOOYXFixfz85//nFtuuQWAiy++GIC2traPd5IeHiME74xyC2T9+vXMnDmTgw8+mPvuu8/WZuXKlQQCAcaOHctOO+3Ebrvtxk033cT48eNZs2YNAB9++KG39fYYtfwnuuatKLdALr74Yvx+PzfffHPNY2vWrOHLX/4ykydPLnq/d999d2688UbuvvtudtxxRwAuuugiTyQ9PD4inlBugTz00EOceOKJNDQ0VFw3DINp06Zx1113FYtrnHTSSVx33XWoqspuu+1GV1cXU6dO5Qc/+MEmmLmHx5aJJ5RbGFJKIpEIv/3tbznxxBMZGBgArCro++yzD9lsllgsRl9fH6Zpcuutt6JpGj/96U+ZMWMGixcvLqY+enh4fDQ8odzCEELw73//m+9///vcfffd3HbbbWzYsIGTTz4Z0zS56qqrWLJkCUKIojd84cKFXHzxxfzkJz9hxYoV7L777pv4VXh4bFl4zpwtmN133x1N01i2bBkDAwOsXLmSMWPG1NjdfvvtnHTSSaRSKW8l6eGRx3PmjBJOPfVUXn75ZSZMmMBrr71mK5IAEydOBOCdd975JKfn4TFi8DJztmBOOukk9txzTyZNmlRRp1JKyTvvvMOyZcvYa6+92HvvvQG4++672WWXXTbVdD08tlg8odzC2XrrrSt+/93vfsdll11WLPTb1NTE9ddfD5RWlh4eHv8Z3tZ7BPHMM89w6qmnsnz5cubMmQNAb28vV1xxBQCHHnroppyeh8cWi7eiHEHsvPPOfOtb32Lx4sW89dZbAKiqytKlSznkkENsK54vWrSIv//97zQ2NjJ58mQOPvjgT3raHh6bPZ7Xe4SyatUqJk2axG233caJJ55Y8/jSpUs555xzeOihh4rXtt12WxYtWvQJztLDY9PhlVnzYOLEiRxyyCGcdtppPPfcc8yaNYtAIEB3dzdvvPEGjz32GOFwmLa2Nnp7ewkEAtx+++2betoeHpsl3opyBDM0NMSvf/1r7rrrLtatW0cmk6GpqYmtt96alpYW/vGPf2CaJt/5znc466yzaG5u3tRT9vD4xPhPdM0TylHIAw88wFe+8hV23XVX7rnnHsaNG7epp+Th8Ynjbb09Kli5ciU33XQTQgiCwSAXXHABX/jCF7j99tsrMnX++c9/csstt3DmmWeyxx57bMIZe3hsXnhCOQr4zW9+w2WXXUZHRwfr1q1j7ty53HXXXShKKTrs8ccf57DDDkNKSUtLiyeUHh5leEI5wnn55ZdZtGgRqqqydu1aBgYGiEajFSLZ19fHGWecwb777suGDRswTXMTztjDY/PDE8oRzmmnncabb75Z/L1Qp/L999/nvvvuY8GCBTz99NMkk0nq6+t57733OPvsszfVdD08Nks8oRzBrFixoiiS5557Lv39/Tz88MNcc801LFy4kGg0yoQJE0gmk4DlJf/DH/7AV7/61U05bQ+PzQ4vhXEE8/e//7348y9+8QsaGxv52te+RkdHB1deeSWHHnooK1euZNy4cTzxxBO8//77fO1rX/N6g3t4VOGtKEcwxxxzDPX19UyaNIlEIsHQ0BANDQ0899xznHfeeTQ2NnLeeedx2mmn0djYuKmn6+Gx2eIJ5Qhm7NixfOUrXyn+vnHjRg499FDeeecdzj77bH784x8TjUY34Qw9PLYMvK33KOH5559n6tSprFmzhgULFnDJJZdsEpE84+LrP/F7enj8t3iZOaOAF154gVNPPRUpJc888wxtbW2f+By+eu41ALzTbTK7ReGOX5z1ic/Bw6McLzPHo8hf/vIXjj32WACeeuqpT1wkFy5cyLEnfZvgVnMINVp9xt/t8eI0PbYsPKEcoeRyOd588010XS9e+yS82blcjgceeIDFixdz/e1/pmflYgDUxW+w7Vcvrgh09/DYUvCEcoRy4IEH8txzzxV/nzZtGrNnz/5Y7/n6669z8skn89prr+EP1+Gva2Hip79AonsVPe/NJ9WzlkjrhI91Dh4eHweeUI5Auru7ee6557j44os59thji4Hl/+2KctZPL+X+LxzJqlWrWLlyZfG/wu9r1qwh1NDG7M99F39LqT9P6qX1APhCMQBa/7IEfvdfTcXD4xPFE8oRSCAQAGD69Olss802/9VYAwMDTP/ClxhY+G9y3V1se/G5gNViYvz48UyaNInJkyczd+5c7pu/nrrx2yCDasUYRjYFgFA0uOtxVsohHnroIQ4++ODiXD08Nmc8oRyBxGIxYrEY69at+4+eZ5omjz32GO+++y6Dg4Ns3LiRW/50J2YmTWyHnWna90Du+85pTJo0iY6ODjRNI5lMss/pv2XJOmiY2AKAUTVu3biZ9C59lbfv/CkACipHHXUUs2fP5tFHH/W6Q3ps9nhCOUIplFT7qAwMDHDIIYcwf/58otEoKQRKMEhsx51pmLMPvoZGBDBr1ixaWlro7u4G4NBz/lgzlpoFww9K1gAJM98MEg3PxUCnSRtLVG1gyOhhRd8rnH/++fzxj7VjeHhsTnhCOUL5T4XyoosuYtHbb3P3Pfdw/gdrbM8zJbDXTX/ksRM+xwkX/9V1PDWRpfmxD4q/jwvMqHg8qjYxceJEenp6PvIcPTw2FZ5QjlA6OjpYvXr1sHbpdJr777+fG264gdAB+/CTZWsdnT5q1vrXVSTz6QvlIlnxsJSY4/t4fcVL9L3SxaOPPjrsHDdnTNNESomqWueyCxda1Zp23nkHx+csfC1vs4uLzatvWDa77uhs8xHu5fG/wRPKEUpHRwcLFiywfcw0Tf785z9z22238eyzz5LL5QhtN5v6z+wHKfvx1DSENg5zUylRdOeHM2aSD/wvs+btFey5555cffXf2H333T/aC9rMWLZsGRdddBH3338/4XCYww8/imM+dyy/vvJ1hICn/mUvXs889RyXf/dxhBA89vYOSCk5/PDDGRwcZM899+SSSy7h+X++yOVfvMuy6bYPD3jmqee47AdPI4Tg8dc8ofy48YRyhFLYekspa1aI2223HYsXL2arWbOIHnkYwa2m4+9otx6sEkotYf0bdNkhC9NaRopqL04ZE/eLcuej91AfruOxxx7j4IMP3mLLub3xxhvstddeNDc3c9555/HC8y9zxx1/5O677mXGtM/S1job0zR5/vnn+dnPfsbGjRvZZadPcchBn+W2X/y7Yqz333+fRx99lP3335+rrrqK1199A+WNSfZHH1JyzNGfJ5PN0Ld8LHXhsZ/USx71eEI5Quno6CCZTBbzWcH6oN1www28++67AKROPRG3DFffEAT63O8jTOkqkAA7HD+Oux65lWkzpvLUU08xZsyY/+CVbH48++yzpFIpfvSjnxMIBMkmYY9dt+H9Dx7mncV/5v1lESZO/B1r165lm222we8L8Yc/3saf776fxvBEVEVDNzMcdtjrvPzyy9TV1XHUQZ9n4wKdp/75JJMC2zI2MB2f4uett94inU7T1dXFVVdexzP/fAK/FiVnJNl52pdoqZu2qf8cowJPKEcoHR0dAKxdu5ahoSFuuOEGHnvssWLF89Zvfc3xuf7+/L+DzuMPt4pMb93GCcdMI5vLctOd16MoCo8//vh/LJJSSj744AM6OjqIRCL/0XM/LqZOnQrAXX98m2CgAYBgoJ7ttvkyQ/G1bOxexDHH7UZmyMfkidN45s8r2WvazqzoeYV4pgtTGviUAIFAgJ222Y29dp3Lo5cvZFxgpnU8kVrIyszbAOyww73F+/q0MJPH7Mn0jv2Zv/hm1vW+RUvdNK649BbO+fG3PvG/w2jCE8oRSnu7tZW+9957ueGGG4rhPAD1h36G0Pbb2j5PKu4CCZBsg+gad5vePZs59ZxvcM0117BkyRLeeOONongXOO+C33DJxac5jiGlZNrErVi+5gPOPfdcLrnkElu78795MUeecggvvvgie+yxB3vuuWeNzQXnXMvFV3zHdc4//OVvuexHp7ranHvxb9nw4TIAMpmholCClUtfFxtPYNIUYkqA917u5MOXVwIQDbSwbcdhFWNtO2sc/37uTV5dvKp4bWpoJyYFt2PI6MGUJtucMJv3n+9HUwIE/fXFLblEoil+pJQ8+vRGzvmx67Q9/ks8oRyh5HI5wAr7AVAUhVA4TOi4o4jssmONvZrOF6twOTZMt1r/6iHAQSgTHdYAesT697XXXmOnnXaqyDP/zvd+DcCi99KO9/reCT+lq28Dy9dY3vOlS5daz1m0iPnz5+P3+/nwww+586a7WbthLZfc+lM0TUPXdQ455BBOO+00Ojo6CIVC3HDZ/Sx/PQdXOL+2Uy+4jtfW5pwNgG//4Frmv7GexS/eQcBfRyRSuzpOt/jRQwrPLEniVO3TiFrZSC8/32/7AVSFRnCHnQD4oCtCLFSZvWSYORLpbmJNE0nUSdSw9zH+uPH+wiOUWbNmcc8993DCCScAcPzxx3PVVVexx/1/qrFVUwpimMpnqTYwhsk2jE8Q6EFLIAt6O3fuXP70pz9xwgkncPTRR/P4U6tYsSa/bbcR5QULFnD0IZ9nYLAfUxoERJivn/wVHn30Ua699lrOOeec4pdA0BcinIvRxnh++9ivOfDAA3nwwQe56KKLOOqoo4pjhvyN7LPV6YC1Sr300kt5/fXX+eIXv8gdT79DMBxlaY+0nQ/AsV/9CRvWLmXJkvfoWf82Pn+Enbc9CU0t/UFSbX4ADL/1hZOr99uOZdQFkZpq+1hxrB0nYNQFHR9XhEok2saGda+xYd1rtHbsAPzAdUyP/w6vcO8IRUrJrrvuSldXF7fccgsHHXQQAJOvLy2r1KT1oRZV74DY8tLPybxj1fRV2jS9U3rS0CRLYQx/SWkEsPjCeZimyRVXXMEvL7uOvt41TN/6CMZNmmPZCPjnYz8sPuegWV/gmSV/IaI00qS1Y2IyPrAVP/nzKRx22GFIKdl7770Z2z2dnsUDKEItjvMP476K1370PmexYU0/r664i46G7ZjdcSgt43184wdzOeCAA2hubqanpwehakzb70TqJ8wi1d9JZuUC5mw7CZ/PRyKR4OG/P0liqAcQROrH0jZpN8ZM3I3wQPnfKICp1aps4xulUAG9PpT/o1eWmdPeXVkaZ9dJ1t86XCmykZXx4s/ptrD1t1YMEoMbeOOF62jt2J6Na9/E4z/DK9zrQTweZ+HChdx4441FkSxHSyrF4HAnEuNAui9+GJyiYLq8ixRF4YVXTLbf7du8+uI1xIc21NikUin2mnIwb3a+SIPWyi7RQ1BFadDfz3uYVCrFF6aehPKhSt/6eFEkAaSEY8Z/jb+s+SOfn3sBAJneKD41jSl1+hIreXP1gySX9XDXPy4m2jyW7T73Hbr7h1g9/68sf/4uOnY6mDWvPIQWCNMe0nln8SoU1UesZSYTZk2nrmUqPn+4bNZZEh3WilKqzucVekP+OYp7KFRi9ynIgPvHMdUehvy9FDT0rBW7NWHG/q7P8/jv8YRyhBKLxZgxYwbvvPNOxXUtMXzh3MR461/pYjowTXG1kcDu51xFw/tZMlnrw51K9ZKIb2CwfyWq6kfV/NQ1XEt8YAMSkwmBWWwd2r1CBAGGOjMc23YSmcEcYH+OOLQuxZGfOp9ssqT+0UALszs+S+fge2T0ONHgGNr32J/IxB0ZMjQCsWYm7vV5Ft17Mavn/5XmGbsxcc/jCKYEs1uH+xsFkcOEgeqNYfvzhfJx9rLCe6Tm/MdOdeS9/VVD9WxcTCDUSCTW7j4Rj/8aTyhHMPF4vBhDCTDjjl+gdwi0dc6HjbkxlhD51/scbcDaig8XP5k0BaFs6dO99bbHER/agGFkMY0shp5BoZHxuck0qm3EtCbbcUQ4jEEYBtc73kubMI7cxgREwxXXxzfuwPhGK3MlPqMBKSCjloTYH6mnY6dD8IXraJ7xKRRTMBCSBHLOh7b9U62PTWSDs026SSHdGKb1TYdUJ0DrTaD1QnyGS6tgAaZfoOQql/9SSno736O5bWue+evZzs/3+J/gCeUIpauri/Xr17PDDpZI6LpOLr+60cdl0Na6e2ayY3OOYimMfIiK6iyWhRCjnm38NL9rJYm3tm9Pa/v21gMSlDIx8r/5obWHrrwTIhJmOLQJhSWwg4GwRDL/I8F+g3RDSSzH7viZ/OsqPSFTpxIYrH1x/dNKH5nEWMUSy6r7ppsKEQSCrh3DtLyZrDkH1noTxZ+jS/sssXRZfJo+gcjJoklyqJNMqo/msbPQNO9j/HHjNTAZoWzYYJ0FTpo0iXg8zvR7flXxuD4uU/nBFKXVZIHs2FzNh7cgkgWqzzCFrI3D7NnGX7lFrxJJgOwOU6Gin06tSGrjOiq3skKURLJwaShZIbhSEUWRdGO41TFUimSBRHvVR8hG7Lp3qBJ7G/9pdKlNClTVWNInkFiryRWLH8fnj/DKMzcMM2uP/wWeUI5QfD5rNfjue++x7YPX2drkWuIkP3if3EAfuVb7s79se+l6tUgWKIilkFbaox29W/stp4+NSBbvtf0UpKaBcF5Jah1ji4KqjR9nayPiqaIYJabX29oE+42aVV41mTq1eA7p5rApvH6pQLrR/iNl+vLPNyVaX9J+3nrZhBxuJ32CbHqA3o3vcd653ycYdA4j8vjf4YUHjVCklOyzzz68vGwx7eefghINYfQNkvlgFZmlq0i9sRh9gxW+EpmzC83f/LzjWL4uH0rG3SkhNUmwa/giF61vOAeZF/Av70Lo7ks8vbsbbZh0SKMpSnKK+/s1E1Nw3fMCWtpkaIK7+z/YZ5ILu687mt5NEdgQd7WJb+W+BQdASrS+R+nu7ubtt9+2NXl/6QdsNWP6MAONbv4TXfOEcgTz1ltvsfOndsPIZB1tAjMm0XrWSSghl5WJBH+Xhsg6C0GmTSeyUkVNuay8FGhcmsU35OwEEaZEyRpoq7sdbaSmIusjKH3OomOGfaSntFixnU7TlqDokkydahsHCWD6Id0srK25wyfF9EvSrZLwOoGath9HmBI1DQ1v9+OL29eiM0Iaycl1YIJwuJnhE2SaFI5oG+Liiy5kaGiIaLQ2B2ivE67gmnP251O77mw/aY//SNe8rfcIZvvtt2f96jWEm0pe1eB2M6j/3AE0n/J5xl1zHm3nnoYSHn77lm3VMYP2AmcGrOuJSQZ61P4DLlUrs6d7Wz8Zh+2pMCyRBJBB+8wW6dOQ9ZHiz3YYsQDpKVb/HjXrpG6WSAK2ThsAI2iJZGH+dis9I2iJJECyQ6KHa+8nDEskAfq3bSDbYP/akpPzH1YFpMOyMt2oICXct9javj/wwAPWPAyDG2+8kYaGBmJ1dSx66joO//z/8Y1v/R8PP/ww/f39tuN5fDS8FeUo4Lvf/S533H8vdZeVFYUwBdK0+TBWvxuqvbV9WjGjB8AImeh1lQIaWq/gGyjZmFptZk/DBzlCXSWBEoZZc3ap9A6hJDKlqQR8yGiowkYkMohs6RxVbwiRHddQYWP6RKXTyQTFqHxh2ahSTD8EK58921D79xFVi8HkuNovj9iHVdk3ydqPWGzZEMGu0mvL1QdId1RVR6paWepBhUxdXrilZOXTdzCwfBG77LwzGzduZNWqVXx67n6s6A2QineRHuwiNbQRI5dGVVV23313TjjhBL71rW8RClX+HUcjXmaORwWzZs1isLOLcE8/WnODs0iCtWoqfDZtvkL1Rh1UFXVIxYiY6NFaoUiNNTE1CPQotiIJ0D/dh+mHyFoDoZsoeu04ZlMMFAVlKIUM+ZE2K18ZCYAiEOksuZYIubbaN7ySk5hSIPPv9mqRBPDHTXJhS4zAXiTBOj4o5MXnYvZrjGyDxN+fz6DR7W2GpsWKQpltDpIZY+O8UoDCvcIK2WhZiqgQTNzvS/RPeA0ztZrDDz+caGMHzy0NMK4s/lxKSTbVx6zGLlLxHubNm8cll1zCj3/8Y84880wUxdtUfhS8v9Io4JhjjsEwDA5ak3UXyXJc9hl6nYEeM2xFskCm1STbIG1FssDgRB+5mGorkgXMhghmfcRWJItTDfmRoYCtSBZQdGnFIbr4iHxJE1/SJOu2uFAssczWS3J19n+kTJMk0yxRchLF+XiYxKQImTFhe5HMI4UgG60UyeJUVI2mmbtj7ngcHdN24bmltbGxQggC4SaWZ2fyl7/8hSVLlnD44Yczb948jj76aAYGBmqe41GLJ5SjgJaWFqvkWDLLL9u2c7UVqolQhyklBBiR4W2cVlwVNuFhksmB3l0a6ZvtXrTXqA+hJl0a9gDB9XEC6929zrmIQA6Tl51uM8g2u79+PQjZBlcTkh0hEhPct8CJdoVUq/vH1AgK7l/a72pjKoKzf3Yj06ZN45ZbbuGRRx7h+eefZ7fddmPVqlWuz/XwhHJUoCgKRx55JP/85z854eDPcvP02sK2BYQiEapE+FyEoLAidQlEVPLe31y9s00hBTA1xjlLqHdWkFSrRnKss40Z9CE11TE+EyCwYQg1nkUbzBBYZx/smWxVyYWtdEHNPtSRVLuBHpNITWIG7V+bmhYoOuSikHJKw87/CQ2/sF0tAsQ7FLJ1AsMPOQcbsM6AcxHhKqhSgReWpDj5R1Yt0MMOO4xXX32VXC7HnDlz6OzsdHyuhyeUo4YDDjiA9957j7feeouD9p7LX3Y+pMZGaCWhEYqDWJqicltuI5ZKWikGp5s+SbbBJhNlrVkUIz2kkhxbu7Xu3j5EprF0jN71qdo9sRn0IctKl2nx2sD5wNpBlGTpupqotUm0acVamlAV/J0nNc7AKPNqS1XaRgKUb+/1ECQ7akwqX4Mm8vGcleQipfmYPsjFasVSD5Wu6UFBsq12nPLQp7fW6Bz/3asAmDZtGnvttRdr1qwpFkb2sMcTylHCUUcdxdixY7ntttsA2GX7HXh6n1KQuaKZiCrRE4pE+MuEoFoki4ali0pKqTkHlJok21gaJ7baRM1U2hh+hcS4klh27RwmF63clmfrNDr3KmXamCF/hUiCFYrjGyodDAZXD6BkarfkoZWVZ3OGTcSOFq98sUag9sVXp3DaxZHWFDy2CzZQIVNXei2Dk+wED7J1ZcIYFjX3N/yCxNjyiANRc78VPZIjTruC6667jrvuuotf/epX7L333rWT8ijiCeUowefzMXfuXF5++eXitelTrCZZimY6bqOFkCh+01kki4YSJelcKV2qkGkyqVtpojh0XDA1hfj4EJ27RYre52r0sMqGuQ3WmE5niSb4Bi2xFDl7743IGoQ+7AdgaJx98IcwwTeUT4Wc5Hz+aYStF60mhePrj08sDOo4TDEffmCKWlEEucImL4x6RDiWuCusIE1frUgWePv1BXznO9/le9/7Huecc47zpDwATyhHFXPnzuXVV1/FNEuf5lBd2vWsEfJp08P4ZURWGTZ3GsU6d3O9lypcz+PAWn2Wx1faDwSBjQ4HjXmEYVr53G5TkpBplO4FjAVocffXL9Vh/8wAdG/rXggZLCEcrhZmqsX5Ra1a+AgrFtzPySd/i1/96leOdh4lPKEcRWy99dYYhsG//vUvALb9y0UIIQmG3ZtqAYjg8OV1TL+7EgQ7lWHfcQNTh39Ljv/bepAS30bn0BYlo4MEM+IcVtS3i1WdNzDgPO/4JECAr1+1/7KQENigWSLoEnZVl2+v4dabKD5RgCIwXBzhhaIj6jAp81II0g02557pOBvee5aO2Qdw0003eXGUHxHvrzSK2G+//Rg/fjz33nsv2WwWo7jEkQTD9gF/5TGXImjYbuXKc8CdxDLUqRSfOjRJqR1HVIqkXSMzYUomPLgOUUgmk5R+Lg2DktarHE7VA4miSDoiID658rm+gdplZbBTK5oIQNiIZd2HVMzHViyrnmYrllV/WiexLL+eblAqVp+5lKW0F8/7ov2TPWzxhHIUsW7dOtasWcOOO+7IDn//Rc3jlWJpH5guApViaVcow/TLko20RLKaoYlKccsrFRiYUmtTLpZKzmT8Q7UVzrWNgwgz39VRSkS69izRDAcr6lj27dxSYxMYkBVb4/ikGpPi6ykwXOdKN5vy61KxOlhWUy6WQoLPJgRUTZfNSdqLZ6beEstcJsEHz9/Otttuyxe+8IXhJ+9RxBPKUcTixYsBuFy3L80FBbEUSBcREAFrG+5WTcj0WZ/e0EZnm/gE6zxucLKzjREABIz7e21TsgJa1yDCMBEZ5+MBMxwAITD9zvfyD1piKV3OCH0DKsK0QoACG52cQAIQCB1iK5zHAuteifEuB47CElXNJU5ezVhCWh1JUMDIZujLruP9J68n4oeHHnqIWCzmPjGPCrxc71FE4cNhZpzPJKUUee+Nm3tWIHLDf8c6FfotZ3BKbTiRze3QwxqaS+aNMpjEjIYQNnncBfTGEENT3DN8tJSk32Z1V2HTr6JmnD3KAFKRhLvc/46mHxId7itTqUj8/S4e/vwd1Iz9Eer6l/9O18Kni7//6sYbmTJlivMNPWzxhHIU0dJibTnbPuhnYIp9yohpCnxBnVxGc8wJl1kF/KYllsNsLU0fjuFAsvxwz0Hf1IwkskEyuH0LTfPtV5Vm0Ed2rBWM7uu1P7gzQirJ8WEUA0yHWr2GH9ItCoE+HPPUhVGKlXQqvWb6JNlGg2wTNLyn4BuqNcpFYHCa9aK1BCg2XypqWhJZa/2shyWmU5X1wtFDWXk2KSXdbz1XFMlXXnmF5uZmTyT/P/G23qOIqVOncsghh/Dh/c+w87raD50sO/X3BXQU1Ua9ysRT+kxQ7OovCpRMvgpPvcSwcfBIpVRVSDoIl5ayRLJAtqnWg21E/EWRBJA2W2s9qpEcXyo8odiIuxEUpMtCagrVfypel14ZUG5XzNcIWCJZoH+WSaax9vUXRBJAj8jiUUU5BZEE0JL2lYgqCvyWlWXb+OqTrH/hb/iiDXz262ew2267MXXqVMQw7XM97PGEcpRx++23MzAwQMO7GzkiVfL8Sikwq1aQml+3gtELGAJZteWWPgllgip0gZKtHCdXB3pZXrRUqIkVlAoV70Z/XBLeWCkM8a0aSLeVBE+vC5IbU3nWpkf9mMGSdzpX7yM1ttaFXO640SOCdJPNiq6q06xd9fLqY4Ncfe05wuCMSmVON9eYoIcqv1B8g7WiqKVALWtbK6SsXdGb0P/+QjpfeZy2Tx3M2RddxiO3X197Q4//CE8oRxltbW185zvf4cYbb+SbuxwIgGkjkgU0n4GqmUhDQer2bxepSaQqLZHM2Y+jRwvbx1qRLI4jAMXyQAd77PfiySl1pMZH0RtC6M32541G2IcZsm6SHmMfRynKkpEy9fZz9sUFWr6rrFu5tIJY6S4VlZId1s1SbZLEOIc2D/kvk0CfJNRlP07Bqy2krFnNStNg/YJHWfXkn2icsRPfPfEELv3e/7lM3OOj4gnlKOTEE08EoKenB1Uxh61PKVQJ+jBbNoGjSBbIxZxrOBaQwtl7WyA1PkpqnLvX1ghp9M12b3EhTEi3DOO4SQpEjppVcvU4hl+6CmVivEmyHZJtrrdDyUGg12XOAr6918QKkZTSZHDVYt6/9yo2vv4MY/f4LJde8GN+fNqJ7jfz+Mh4zpxRSH29VVjiwjtuQD1sKn6fQTZnn6NnmgrSFChBHTPt/HYRpsAMmMWzSTukYq2atISz6PgSoAcEWspZUIVhnSmiqPji9i7z/hkaqWYV/1DlOWc5ibGCbMxqHKY6rBj1oOVokWrtNrtoE5JIzaqa5NRXSEkrZJosUXVKmfQPWCFFmUYI9NU+LgT86NjZfO7wg6iPPMZ5tz9G3wevM/DhIvTkIOH2yWx13Fn84msHc+zhB9vfxOP/C08oRyHjx49n33335aUXX2TaYVNRFRO/T5LN1b4dCqtNRZOIkI6RsnnL5FebUgNDmKhpu8o3JbHK1Zn4BmttfEMlD3k2KvDHbZwXeikbx/AL7Aqo987ykanPO5NiwlYoh8YL9HDJe20EaleyerisOpDAViz1SCkPXNEFIqVghCrFUk2Vys4J0wq/qs4dD/SVznb1kEAqkqDVTRgjl2HFy/cyuG4JX/hTFsMoTcIXqadh+o40TN2ecPtkFEV4Ivkx4AnlKGVJvBO9rN2rqkgCfp1MtvSWMIyqEmaqRA3rGMmyt40uEOU5cqrVcExNlZf6khVebdMH2XoTf1kDMv9gZeMuUxNkY+AfKivhppOP8SyRrdfwD5Se2L2dj1ykct5DEwWxVaXnDU4SGIHKVW31Kq9CJIt/gCqbqKx5njAEapVYVseTCinAsOpZFl9b1dbeCAhSrRDo1Fn27O3Eu1dy6v+dzHbbzkZVVVRVZcKECTS3tvGNm56yxhVw2ykH4fG/xxPKUchBf74IsdMEMje/R/ztlUS3tfL1FFESS8NQsCtRIxSJGslhJHyQEwi7uB6lJJbVIllAapBtNPH3KfgH7Le1pirI1EFgUNqKpGUDmQaNQL9O105+9EDtzfSAYGAy1K+QDEwRVvkxG/QgaHlniVO1IKmVBN1pC10ujOUdKytsynwxwW57G0OTfPjG/Qx1LuOOO+7gy1+yz8/+6w/qOeZXD/CX7x/LlMlOuZce/w2eM2eUMeeBC+gM6zR/djdC08ey8c8vVDyuCOkokgWEACVo2ItkcSDQY/ZFNApIxcpOccvMKWakuHRVlgpkGjVbkSxg+gQD0xRHkQRAWBXJc1FnE7DEcrh+QGpCQU0otkUyCiiGILTROTNpaMlb9L/9Ctdfd62jSAJMmTyJp392oieSHyOeUI4ynv7suQgEQhE0f3ZX4m8uJ7Ou5GZNZ31omktOHWAaYvgalqr9StLGklTb8IUaXQUOGJykDustl2L4QhZSYdjam0ZADltb0mjUMRrdm51JRZJuqZ2QlJKBd15j/eP3ctBBB3H66ae73wxoamoa1sbj/x9PKEcZwWCQhZ+9EAVBw5xZ+FrrWHPdI2S7BkhlSicxms8m7aQKUecWXGhhhpyVqZja56KBwf7hS/QMTlKLC2AnsRQm1suRzivYiq20zW0l+XYQhWk7lEvTm0o5m3pTzqH1gyzap8eUBpKGzrqH7mDdw3fyhaOP4q677rKfrMcniieUoxBFUXj9iIsJ+gNMmHcUicVrWPzN6/jwW79C7yt1KNR8tS0izGrHRF2WakGVVamPdmJZnf+caq9cgQpkjUjWrCqFYHDS8O1uP0o5NNvzxvJSaGB1XayeQtXYemNtYrvemKs4yaj++0BJLHtffYHBxW8x+eivcffdd9PcbJPG4/GJ4wnlKObfR11M+1aT2erGUwEw01mEvzLgRivrp1MtkgVEXa5oYycCkBfLfPELuyIRYGWtoFjhP4F+h3HyYikVweBE+7evminpmWM9yLJVpWsriMJ9HVrTliNt8t4LGAUBdVk9p8eYGJkU/miM5X/9w/CT8vjE8IRylPPisReT1S1x7Pj+Cag2rRNURTqKZAER1l2FAkBq9vGT5eRCOUKd7q0pTJ8gVVt7twIlg2vQOlhiqbgfI1qFJoYpAydM6wvCaBj+TNLtbySlxNywkm2meRV+Njc8ofRgYtx6GwitdhtrGoLMkL+islANugK9fpSU+9spskrD597vi0l3rKLhpdWocefzTzVt0P58P/4hZwXzJU3a/p3Gl3Ded6tZSWyliZpxFi9Fh2CXdM31lprEcDiLLBkJzFBt3GU54rHH6fvwfS688EKXgTw2BZ5QevD4mZcghCC7vqfmsWzcangts4pzTnhfabvutNWtLE9mP05oTamWZP0ra21ttIRO7APrHHXMywO2NoFBg+ZFllen9fUMgYFaQdXSksg6SyBjq6Tt6lPNQqDXuh5eV1tNCED6TfTm/ApYcfD0m6JYDV767cWy4cWXWP7U01xxxRUcddRRtq/LY9PhBZx7EAqFmDFjBv4P+ylfOBnVJdVyCvhMRPn2MVO5ChVpBQJmRcC2lhAEO0sXtBToQYEsS2uMfpCg7R+dFWP5NybIjilVCPIN5IisTlTYhDdkSbb7i78HewwaPqhc/jUvyrJu71KpNV9CEuqqFMboGkn/jLJA8Qz4q7ozhjoh3Qp6fkpm0MSor9puK7KyB7pZqgZvZrP0PfAwqbfeRWg+lHAIJRhCplMsX72GefPm8b3vfQ+PzQ9vRekBwPbbb09LTuXYzBgA9KxKLlmbSS1zCrKQ2phSYbD2u1ZkFEQ+/9s3WCmSBbR0qdpQ/aLBGpEEiL69kcCaQQACPZkakQRoeitOdJW1Eg136jUiWSCy3hI0/1CtSBYoOJC0VK1IFgiWlT+rEckCZV8k5S0zEv9+g/jz8wnvsgPhnbfDP34caiTMzEmTuP3227nyyiu9wrqbKd6K0gOACRMm8Pbbb/OLL53KX/96IXrKpl9sHqkLRFqDpHNojsgKQhtBSTvbqBkIbZS0PN/taBN5vwejKUpovc2+N0/D4iSKIQhvdD6PrF+WQ8lJhOk8n1CXRBjDO2/8/ZAeN4yRApTlb0vTZOCxpwhMn0LTcUcC1i799Lp2vv/lr7iP5bHJ8VaUHgAIITBNkx/ffQPKMN5rkVZdRRIgvF4QWqegJZzHCm8wia4zSU92zyrxJU30OmfhBoitzKDk3MUr1CXRXBw3AHUrdcIbh/Feq7imeAIoKaWm5JxMpxE+HzKXQwDnd0z1RHILwRNKDwDeffddEn54KLwRgGDMPsVFSakIlyK2BQr1FN3SCoN91gpQb3UuwpudZMUBGWEfuQb7Qrwy760XuomStRdLPaJi+hWUrERL2a88AwM6WtIk0G8SXWcvlqkxkI2BSCqoA/YbMiWhQkYBQyDyxwtCUWj+2vGkl37Ihqt/y0+axvPNo45xfN0emxeeUHqQSqV49tln0XcrFVXQ/AahOhuVG6aKOUBsZaWNb6h2FRdbWSloqa1rS39npo7BqC85YcygTb1MTQWl3KNeK4K5qIqpld7qim6tUssJ9umo6bKeNXGT2OrKeM7kWKv8WvFeaQW1r3JOypAG5f3OTYpiGd5+NmPPPp36RIprL/0lyeQwsVIemw2eUHrw3HPPkU6niew0reK66jMI15dCdpTE8OmCdR+W+swUn5erdI7ULTfwJSvF04gFSW7bUfw9vVU7ZrR2u51tLimV9FWKZPF+mdJqMBfTkGrt27z8HDLYo6Nk7Zp5la4lx1nFfWvGySpoPZbTSxnU7FtmmBRX4YFJE3j66ad57733uOqqq2ptPTZLPKH04Ju/vhCtuY7AhNaaxxTNJNKYQolrMEx2Tv1SUWx+VY0wwN8vafhAR0vbnxOaIR+JHcaTntWBDNrVLgfpU8i2RpA+zar3ZncvU6KkdXJ1WqlMmw2+hEmoO2fbBrb4mj60vOiG/XQsCn8Xt79PPsTy+SO/yHbbbccZZ5zBr371K7q7nR1ZHpsPnlCOcvb6ywV0LfqA6M7THENTMlkNpdVBAfOENijkIu5C6ktI9ID7W65vmwDdO7g3BUu2+RjYyr1opJDD30uYYNr0AS9H0UGqwxw3SPCv87vbAKF1KuPHjQPg3HPPxTAMbrjhhmGf57Hp8YRyFLPb335Cf18v2TXdRHeabmuTTvvyhRydxwmvV6yaGC42vqQs2mRj7g3IpCro3tZZeKSw7NzEMj4xAoKalg/lFHvvhJzn07u1NY/oGpcCvDnAxLVpWrBLQZqw/XnWdrulpYWvfOUrXH755XznO99hw4YNjs/12PR4QjlK2fnhn5BTJUMLPwRF0LDTxBqbdMpXEQajtLuvKgEyje6iWqA6ukYK6NuqzDHiMEaytWQjBfTPrBXL+KRIxfPtxFJUVUyvEUsBvbMCpe29hOhqm/PQqvoddmIZ7CqNnUMw+7yrMQyDyy67jK9//evcdtttnHbaaTXP89h88IRyFJJMJjHzZdEGFy4jsvV4eh5fyPsnXoksCIgEO7VS2tOlwrXSWk1Wk2moFMtqxw1ALqoUxVIqVnvZarq385dsqBTJcgrnkFIISyTtKK916dBWoiCWUoHere3jNsvFslokC5SLZblIFjCB7X56LT6fj3333Zd4PM7g4KD9YB6bBV5mzihkr2d+AYDUDeJvrsBMZUm8uxqA7JIVBLeeQirtvPVV2tKYG4OE1jp/z2YaBP5BWdFFsZpcVEFLmgxMc34b9mzrp+m9LOlGZ5uBGRHqPkyQ7Ag72hh+gaJLV8cNAAL6ZroHt4uPUHpNSwiMkFtJNdj1F79j+w8WMGXKFG655Rb3AT02Kd6KchSTWLwWM1WZGx2a1Eo4kEVRndMBJSAbcph+94ZfiXFguPg4TBWGJmrD9qjp3s6P1Nz38xt3i2EO87U/OFHDGKb3TjY2fAhUZJ0YvkalYb+arMCAf7y9glgsxpQpXg3KzRlPKEch0wYtr/LQwmXFa1pDhG1uP5NQg6VsjTH7YGgpwcxaYpKaYV+AQiqgR/Jb+6kC3cGJnWm03n5WK1p7m4L4JdoFpt9e5My8iCbbVMcmZIXrQxOcuzXqAcHQBM0qg+agp1JYfckDA8Kx6K+iQ6BfWJXWXSqsa0lQVI2uri57I4/NBk8oRyEPfOl8dhqIMfDykuK1+j23ItJmne8ZyQzdT75F+rW3ap5r5ipXXEa0VgkKIllgaHKt6hhVomcnOtUOn+QYgRGsvGj6Bdl6UWaj1NgYAUGytfRWj4/TyIUr3/q5sEJ8XKWjqFospQKyLJ7SPyiKFZAKqFmBf6CyXFv16lPRrQpFAJmBHhLhFp589nk8Nl+8M8pRyqUHfJ3b1/wIgOhOU+l54g30rgEynQOkV5cK+M68ZwY5aaUR2hXuTU/KEVijoQ1YAmo6BGZn60UxO0cPCHIx+wyWwle3FFTUtCyQahFE1+TDegKCXLR2nFSzQnStpU56SJBurF0PJNpVwl3gHzLJxhSSrbU3K0ZFybxI2nxa/EOQiwiMgETNCHzxWhs1a/1dTM1yAJXnvyuaDyOT5pwHX+XcwSGOP/Kw2gE8NjneinKU8o9//KP48+RzjgBTMvDqh0WRDI5vYrtb/o+6eoNgIItpCkzd/u2SGa+TazYw/WA4NOFKjIV0k0APOYgkoBj5/jOKvUgWyEXz49iIZLlNLmovkgWSrSqZenuRLFAQbDuRLOBLQKgHW5EsoOSslWR1kZDIuGkk131ILhHnF88s4bd/us95EI9NhreiHKU8/vjjADSOaSaSD0UUqsKUcw4nunUHgfaGom0wkMM0FFK6szc4266jhiVqv0uNSlPmt6HOAqcHLVFy61GT6LBEx2/fCcKyGWsJZHXeeTlGEOLjVMe0ywKKLjGGyc4Jb5Ak2ofJTIrL4nlqgYaZO7PhhYcZWPoGLTvuTVYfrtuZx6bAW1GOQgzD4O677wYgst9U4m+tpOUz2yENk1xfvEIkDVPBlIJgKEsk6q4oRp2B3mz/QY+uloQ3SvxxSaDf3sORi4Eelhh+aVuAAqzzTyMIuShkGp3d5XrI+i/nkulo+PP/hZxt1KxEyVBRWaiawpFCwKEqOljPV7LWeOVooSih9okk1i7jB/tM5zsnftF5Mh6bDG9FOQpZsGBB8eeh99az5s//BsDXFKXlgO2Kj1kiWXqeP5AjEXfPwzYjJlT1KKtbIfEPlpUwS0qEaZJuqvye1su27aZPgqgtslEuoHrYagEb7KlcpWXLylsaQUABX1U8d3m5NCOfgKRVOfrVjCw6YpSstSLWw5X3CvSZaPk5aikJpix68wtoKYnIf38IHVQpK7KFImMnI9a8y5c/dwQemyfeinIU8uXv/1/x54G31gAw8St7sMN1X0Krs5ZXepVIFojVD19DMddWmbJSLpIFqisIZRpqbUyt8lqurtbGCEKqtXQ9Ww9mVeym4beuF9Aj1MRcmlqpaRhUimQBoVNRsT3YUxLJAloGgr2lFbOWLIlkcRyDYotcIeCEOdsTHxqqeW0emw+eUI5CNqwotYKddOJeTPq69V+4LUZQsz7VDll+aJpBXYPLwR8gA5LcWEssG5a4tILI97jJNElHb7ketp6frZeOAeUFYcw01Apg0SY/fi5qBbrb2qjW41q6ViQLCMM6awx1magO56hqFkLdJr7EMONkJHefciivLJjPzjvvbG/osVngbb1HGVJKMp3WPnTG2Qcx9tDtKh4XQqII91QZVTVpbRmkq7vO+T4+iZoSaBnnDB9Fl8QnmahJZyeIVCA5SUcdcs+YGZpm4HNxJBXu51afEqD1zSzxsW7FJ62OjrnocPeyBDUXdnltBhxx4jyW//sVnn32WdfxPDYt3opylLF69WqkYbLNRUfWiCRA1lARQhLQnHP0/KqBopqMaXV2O4eWWsu8/mnOgrJhT4FUQY86C3O2VUeqEr3B2RtsREykBtlm5zmHOyVICPa6fwkoWUndaodqF3mEAf7BYXIYsYqGuKEnhlj+6pP84Ac/YO7cucOO57Hp8IRylPHggw+iaRoN246veSxrlETNEstacfJrRql6kCJpG1MrlqElAcgHp0sF+qfXiuX6OaIUKylAj9WqSnaMXhFJ5CiWovSvo1gWFrYuYtn+7/xe2pTUrbIXy2LTMQn+geHF0u/QhVKYkOpdjzQNHnnfCwna3PGEchQhpeTGG2/kmGOO4cVjbkAta11QLpIFqgue++1WmUJSnqgdWlIb1yMF9JWJpV16INiLZY1Ng14hjDUplAKyLUbFnMLrq8aVEOyRpWEktL+SocJ7JSV1q7Llt7LtzPhRxbL85Rbyv0ON7Si+ACtf+is7fe0X6F4M5WaLJ5SjiJdffpklS5Zw+umnA/DMQb8imFVsRbJA0KfjoGtF2toGEEKWVm125KuSmypsmOMyWv5m2THOoqHXW2JpRJxvmG0x8/UyncU30CNRdJO2VzL2RTkkxFZaYhlxaF8LH00sfQmrwnt5kQxfuI7pB51Msncdq+Y/yB6nXjvsOB6bBk8oRxHPP/880WiUvffeu3jtDzufjayuPlGFX9PxuZxZArS0DBJa7t43pn+mQucc9znqdaarSAIgBUbYTZVBG9KJrnAfR82YTLrPvQWDVIVrsDkAApTc8KvhYL9Zc24ZbZtC27b70Pfh60jT/TV5bDo8oRxFzJ8/n0996lOoqrWCfH/ZUo58/RpMKRzFUkqBdF1PgmEo9A2Fyc1OQMD+wy79kvRWaWjMguogKqpE1ucQAcP5nSkFSlKx2r86DOPvzjL9V8tpffgDtCH7GB4tYTDmX50AKDnnL4FUixUY4pZ7rucro7uKZf6hwGDl38fIpuhfuQh/pIGx7v3SPDYhnlCOEkzT5OWXX2b33XcHYMEbr/GVpbeWHrcRSksk3dF1lf54KQcwt3USGa4UHjNkkpleFpndkAWtSlA1E1lXcqAIf61YCtMSyeLvWWG73Z967crizy2PL8ffU5t62frixuLPgdX9Fb3Ai/Mu69CYqVdtiwfL6k+QwxZeLRNRxSj9vPL5P5ON97Hf507kkWvPsXmyx+aAJ5SjhAULFtDZ2cnBBx/M4y88w5kb7q+xKV9VfhSRzOY0BhK1KY36jFTxZzNqkJ2SqbGhPgf+ksrJWK1QCb9RsfoUqdq3q8hVimV0SW0wfNMzKwmsK10PdtXOJ7B2ADVZEmojoJBurFxGZuqUisLAUgUjWDknRZcVYinMSpEE8Mclqi6Jb1xB/8pFzD3sOB67/Vc1c/LYfPCEcpRw7733MnbsWFYnurkg/qStjSlFcRs+/IkbDCWdqwmZTZboZCe6lAGKWTbS73w2J3zWY0J3aTubE2BA/cJBxt+5ztam8cU1hFYMEF6XovH1Xlsb/4ZB1KEMekgh02C/187GFIyAQGoCw6FvuKJLq6+OiWOPHl9C0rfkVeoamnjqgdsdX5vH5oGXmTMKME2T++67j+OOO45rxGuu345SCnRTQVPdnTfprHv2ijEhg/lR+taGdcdzzQJKRgEXoQRofNeg7cFOV5v61zYiYu4Hgck2SE7UCPQ5f1Wkm1SUnLTtLllAS5kgZcX2vZqIjBLPDN8C2GPT460oRwEvv/wya9euJdHiR+ZXjU7kDJWsrpLTnb0XqYyfVMZXE2dZjqmr5Mborh51IUBEdYTiLDjKoGY5bwzncZrfytG8SMcc0+g8IUWx/ss6r3B7dq2jd+c60i3O9zKCAiOAa9FgsDzqalaiuqRwtgYmk0ql+MwBR7mO5bHp8YRyFHDvvfcSaazj9Z2t8B0nsczqKrphvSUswazdcCTSAdJZ67oQ0lYsTV3FLAibxFYshQLkBVKoEmHjCVf6fYiEJdjCAMVGLMf8O0fDUmv1azbFMMY2105IUS2RBGROh4y9WA5sU1ptJh2K8JZ3lczW2X98tFRJHJ084f7eLPXR8Uwbvx/PPPsIe+95iBdwvhnjCeUIxzRNbv7jbYTnzkQolc6acrK6imFW1Yc0FDK5kljGUwGyVc3FhJCW6OUxcmUiWbxZZb8dociaRGihSJQysVR6fbXOGwOUqi14bFXlEYGsj2BMaK18XlUhDKnrkK506Kw/sKni91xMkBhX+bzqWpR6EDINVbUnk2bNuaSaqlxVBrozaElLFCd3zGXq+P14af6THH744XhsnnhCOcLZ/6KTSfUP0bj3rJrHjLx4ZXStRiRLNtb1wWTQcTtecP0YWc22AVmB4mNOJopEyYcNiYzDW9MsieX4Zxza5UZCGJPbrV9U+zlLw4CUdT649rMtpDpsvPdhQXyCda9cRNiWcDP8FAsQ+xJGRehP8WUZEi1pCXpwYxo1XRJ3IQRjW3ZE00K8s2hlzXM9Ng88oRzBHPzni3h34Zv4WmJEth5nYyEwTIHpIm4AmmJiGO5vFUWRjjUsK+7och6ZnxJql3uGDybUfyAI9Dmf/8mgn8G9p9QmrJfbmCbrDh5Dptn5fkZQ0DvbPeDc1CA+TnHs4Q3WKcMdlx+Lki0ZGWaO5WufY/5bN6AqPqa2exXON1c8oRyhzL3/p6wPZOl7cQmNe8+q2HZXIoZNTxRAS517sd5cxueeEJ4fyBxGcNUVIZSU+0B1y92bjwEkxgUwQ+6e+bVfmIIe1WpaQJSTi5qYARia4tafR2D6BD2z7O8nBNz7x28wdeoU7nr4NADSmX7mv3UjH679Jx1jdmaP7U/l/se/7/6iPDYZnlCOUAYCJol3VqH3xYnttY2rrQBHsawLWdtTRUha6+37sebSvtJq0lmPizi1vVU/DEGhR41DSci6DylWDe/b2r4rWHx8sFigt3/uZFubNV+aihnMLxMdNDAXM4vzdlpR6iFRstEEPbPtV6etra3Ff//2zDms7XqdrJ5gj+3OYOakQ/FpoaKNx+aHF0c5wpBSssvDP0WaJutvf5rAxFbCM8ehGwJNdQnsxiqjli2cQwqoD6ZrbMbUx9k4UPIOZ9M2qyhBpfjYesaV4nkkMi+SVSg5KlpE1C+rHFcK6J0Voum9UiZQfELVWaOq0L/vVBqe/bB4afVXp9VsybVEZc8cO+0cnC6p+6DqxVT7raq+AxQB/3j4bP7xj3/w7LPPks1mEUKwpvPfxMJjiYRaUBTBo8+fa3NHj80FTyhHENlslk898XOkkHQ/tIDUsvVMu+xEhFvAYxV+zSBnKNQFbdIO8xTE0lYkCxTE0jXWUkFRTVuRLFAQy/oPnMfpnRWicXHK2UBA/35TaXh2eb6Mm/2kCmIpAT1m/6UyOF0SWyYQstYLXqB7Oz/Nb2dR0fneGTuw9957s2DBAsaNG0ckEsE0TY77/JGsXdqO6hM8/vxPnOfusVngCeUIYt+Hf470S/TBJBvueo6Wwz9FZOtSJXPdVNAUA/fDROkqkgVcRbIwkqEgqotfVNs4ebfL0NJWNJFbNbi+WSF8cenafqHvgCnEO4Y5I01Bus103I6DJZbh9e7Om87pWXjyTo444nz22GMPnnjiCT7zmc/8R19aHpsP3hnlCOKlYy+mOaXQ+/hCkJIxx5fqTgohifizVi8cFzUJaIZr5g5YYUWKNkyCoqz61wYlI6h/JYDi4ktSc9C8yCTQ76JKAgy/qOkTXo6pQf90xern7TBxqUqSEwxMv3QUZSnADJvEp+k17XRL40Dnkud4f8UK/va3v/HSSy9x0EEHeSK5BeMJ5QjjmS9cRGB9nMDYJrS6cPF62F/yjgRUw6pIXkX5NSexzBkKnX1W90XhJJbSWk1C/l8bPVGTgrrXLMdHaIO980ZLQdM7JYG0a/0qFUski/O2cboYARiYVnqr6zY7fdMnSY4vq2bkr520VCyRLJCYYolqxTgaDAX7GHjueSI77si8ZUt5Z8mS2ht6bFF4QjkC+cU35pFeuZG1Nz1B37/eJvnOcnKDlWd4/iqxFELiryqEYVZVEcrmVLr6YxU2QjMrRbdMJIuXqsRSG1SIvVHpHQ5tBLVsx++LQ+OSqoyWQROlLMtPKlSUPQMrU6bcAaSHYXBy7du83OliBCSpjtoVa7VjxgzV2iQmVf7N0mNMhKKAlKjhEBI46rG/1zzPY8vCO6McgXzxi19k+fLl/PSyS+h55N/WRUXQ/OmZTD51f/xNltfar1qOGwCfg0dcSgFCksn66B0K29oIVYJhZd5Ui2RxHENBqCa+XoXIEvvzzWA3ZJoswaxbYT8f/5BJLirQgwLTppAuWHUj/UMmahaGJjhkHAWt+wgD0u0Or90nIScQJkiHbTaAHjXR4gpG0LJRYzFCM7cis2ZN0eZvTz3FUQce6DiGx+aNt6IcgQghOO+88+j4xgGli6ak51+Lee2LN2KkStHafs1wFMkCiWTAUSQLKCkVEXf/3g2s04gudncCBXqtDolu+OKSnIPHuUA2plgxji4YAdAj7veSPokSy2EGnO1SY030Fp1MU8kmOGUy6Q+XM/Dc8yTff5/vzn+J2x75i+u9PDZfvBXlCOVX9/2eDXc9R8Pu0wmNb2Rg4QoynQP4x9RRiA4XIt8NR0jHcmhDiSBD8SBCsa/wU6TXhyLst6cFIsusA0Tp8q4z/YJUVCBVk0C//f3SzQIjBCigOpVzFJCNCtQUlq0NRlBihPLdIx2WDFo0hxbNIpM+9KRD5o3fIDvRhLgK+S+L2Jw5pD5YRu/fHrL+3prKj755BN2pON///Nec/wAemyWeUI5ALvjzjdzxxlPkugYZc/5x1M0Yw+T/qxWw8p7VTmI5FLcCuAsFLezEUvTlBUSCGlcxorVel8jykpdFmDa9ZgAjIIrni+lGhUB/7TjJNoVMg8jbOwilKIUSaSnyMY+VJnpYFleJwhRIZI1YanVZtLDlZfKFcwhFkotXnq0qAQNRqNAezff5GdRQfD7av/VNzFwO0+im/8ln2XjzQ1wVC9ObTXHpl0+xmbjH5oonlCOQh6Nr2fjXV4huN5HwtHZ0EySiwllT7fW2E8u+gUiFjTQtBSqPjRTdfkiWiWBSQTXBqCvdK/a+ir+7rClYISGn3KESrK3Ok2xVCHeVeZk7FLKx6tJnVORqS5sm5Gq6Uij1iIlZlWkoTIEsK9jha8igBivrQ2pBHSEk2SGrBYYS1BG+qi+OsGFVwOi3FF8dI9FCDbSefCRmMs2Gq+/h9rMV+vQ0v/36d/HYMvDOKEcYX7/rUhLvrSG5ZC2tR+9evG6YgqyRL4Ir7BvQCkoC2tMXJWUTVC4lyHyuttgYqBDJ4jhpBbXful73rlYhkkUbk2L8pBGyL2GWiwoS7dZzhybUiiRYWTu5vJ5LBccYSd9QfsxYrUgWyc/H35SuEckCaqDsy6ZaJAsEy75IQpa9UBTavvN5wjvOYOOND/B4Yg3Hew3Fthg8oRxBHHPPxbxVH6frb68Q6GiibrfpFY8bpiCV87kGigugLx4mk3XebEgJ5oAf0s5vH5FV0Po0fP0ud5OQarePfSygh4TlmHFx3kgNMjaFzctRcuAfdD8fFVIgDIHid6+mFGhIo4QdqnYUaM4iWisznISqMua0YxB+H53X389rAZe0S4/NCk8oRxBX7fV1Mp399L+0hNajdqsprWZKZdjak92DkWI7CEfiGhgCI+zuLScnGJrsMswEq/Bttt7d82yqEOwexkaTZJpcTfDt1Uv9pH5Xm9apPdSF3Rt+TW7tYXpHl6tNU8sQTc211ZbUaJi2M44l/d4Kvr7M+/htKXj/p0YQkydOYt9XJWo4QOMB21U8ZpqloO94xn7v2TUQLVY09zlsPYn7Sp6Sj/DusdtSFx8r+IBUyDbYC2FgIO9wMSDUZW9jBArjSNIOK0tt/x5EyET1mTRM6bO1GTOjG9VnoCom9RH71d7kpl401cSnGcwYu9HWprlhCFU1URSTpvraOp7p91cjhGDOXnvZT9Zjs8MTyhHE4OAgd955J/NO/TZasNRz26xq8yClYChT25O7Om3RViyrtMqIDLOqBAan1V4bmlI1rM07MdCfD90pTrDWxghW9d9RJOmWyklqn+muOE9UVEnjtMre3m1bdaGUxZMqQtJQJZZTmntQyhw+qmIys6OyRW5L41CFjRCyQiwFEHjuHb75zW+y//77174gj80STyhHELfeeivJZJKzzjqL1w77GapbcQsJQ+mSWHZWpSYWqBDLIfvak0bUdHSigLUAHZhO3lsEg1PsKwGVB2wH+qRtjnhoY+liIRPGbk7WvxLfwd0Vzc+KJmX3b5vZZduiQghJY8RyqU9t6bGtziaEZOtxGwBobRp0zKFvbogjBCz4zOl0dnay00472c/dY7PEE8oRgmEY/PrXv+b4449n3LhxCCF44cDza1aT1QymA44iWcAX1O1Fsvz+EXexBEssB6fgapdplI6B5gBIawvuKJJ50m0S7aAeV5um6b20b73RVtyKCJjW0u06DoBwq7mWpy2S4Pbbb8c0TQ4++OBh7T02H7w4yhHCgw8+yIoVK7j//vsB6OnpYd8Xfw0IFAchkFIQT1qrStUpjVFCdjAAfhORdRHd/MpSjdtXC0IBPWZ5k339qq2NkFZVoWQbhDdK2602UHQQ2VUTKuDftp+sruLXDBddlqhCYrisvFUhaQ8N0ZmOOVZUEkh2mriGtfF6R5toWuH2nU9j4uETOeaYY5g6darz5D02O7wV5QjhqquuYu7cueyyyy4sX7kyL5IWdrpllomkE9IUlkgWfvc7KFeZNhhRe5uCSALkGgzbd56aLA2UHCNs+9TIsoWtUx8bUTbPrK7adocUQhLyWccKqsMXiU8xmNlgOWzagkOoNqtGgSSkWaFC46IDaDY2rXEfTx/1CxobGwE47LDDvNqUWxieUI4AFixYwEsvvcS8efN4451FHPnm7ysel1Xl0gxTIVElkoahVAiKNAS5IRvveLWo2H3eq9McbdIec/WVy0G7nWuqtSq8KWCFFJVTLZYiZODberDiWtZQMcumIIQkqFU6qqrFMqjqzKiv3HKPCcbxlVUZVkRJJAuMjQ7iL7OZNBDkkWN/BljONoAxY8bgsWXhCeUI4Oqrr2batGlExzTx1WX32dqUi2UyZR8eZJolsazOaS6O45NWip4LRlgWy5JJTaI7rDILNsIA1aFFrZnXcyMECbvW5JTEUkR1fDOGbG1y+awkRZg1IlmgIJYRLcvUOvvzzZZAgoCiowqToGofdN4WGSKo6szui3HvFy4sXr/yyitRFIVtt93W/oV4bLZ4Z5RbOMuXL+f+++/n6quv5qzup9y74UhRjJN0wjTde8GAJXAKJjLnnFJjhiRKTmK4lCfTYyZaSqD1OI+TbhRoGZNM4zBb1Tod32T7droFVMWsKU5cTUjLMSlmH2dZoMGfYkh3P7YYE07wufZ9AXjrrbe4+eabuf7667n22muZMmWK63M9Nj88odzCueyyy2hqakKMbUDKbnDI4wbQdZWsS2oigMwqGBnV1TOt+E0Uv4GZBJlxEDnNtDJ3cgoYDoMpEr3eROigDdiPY0Qk6bESNQVK1sGZ0phDGZPGNBUUxV7lg5pOUNOR4FhSDiDqy5IyfIQcVosAOamgCRPdLvgzjyoM7g09yILz53PJJZdQX1/PpZdeyplnnun4HI/NF08ot2DWrFnDbbfdxv7HHcFvfFY/VymFbbhLLqeSy1n/u6vbbheQGRWz4Nl2aDWrBAwUnyVGSljHVCQyVfU28pnIvA0+E4QCenVJH1k8u8w1G6CC1lsplnrMRI9aNkbIavilZmonpYyxUg6t9Ex7sSz0DHIrKRfSrILG8VwAUwoiWrbGxsiLo6YYYGIrlpow0RQD05D88rJfctppp3HNNdfg9ztV4/DY3PHOKLdgrrjiChSfxvLDK4tfVItANqsVRbJAtUyY6TKRLA5U+asSLIlk+bWKp/jLRLKAZlqCWfxd1jh4cg0GudbKs8OCSBbnWCi0W37/tsq8bNMUGFW56hF/5eqwvEpS0caXJVR2dpnU/cSrtte6VNDL/raaYlQ4d8DylGv5a0IVhNsjvLL4355IbuF4QrmFsnHjRq7/zY3EjtgNNVz5gS7fXmayGrpuv60tfOTNlIbMDdPvOqSjOPToVqJlQuTUx7tcGB2cQXqs9Nxco/04pl+ih63nKx0pREPtFlnKklhGA1kCNs6bcrGM+jME1FqblO5jMGcVLtalYhtvqQqz6OX2K0ZFCFGiM05mIE18TJIL/nqR7evx2DLwhHILZb9vfR6pCJo++ynbxyWQSvkxHESygABk9ba4eiDsK5sXx/CZKLEsMujuKMFftbK0IdORI9tsumbeSJ9EmZhExBwKd2CJZSyYcXXeKAgubzq2Ipynmqzh42T11GGD0s8OnohS5QVb9Ps3UPwq231jJ5Z3rHZ8vsfmjyeUWyC9vb0sfmo+zZ/dBTXm0BAGKzZyuLJqyuogvl73o2otkquNn6wiFM0QaXCvr+iPZPFHa8/9yjH9YETcBTfclkCLuteDnNLUQ0uotnJPAYHgNxO/wV677O5q88u2M9hj5z04r/EcR5tzG89mTxubzECa6NgogXp3D7nH5o8nlFsgu552DNI0GfO5T+HX7EUllSydiTmJpbIyCFmB0MHfbS+WWszqFVNwgtgRDGURikRRTaL19rUc/fmeM0KRzmJpiqITySkLKNIeR/FZPck1h7TL6S3dBDQr1rHBpjiuQHDnNt9lm622BuDGaT+2tblh2o+YlbeZtdUsTMMk2ZWssZ211SwAZk6fSaonRe/7Pax9cTX9H/QRbA5ZgjvB23pvyXhe7y2MPe/9MaseXUDzwTvha4hg579OJQK2Xu1ylOWhmhJm/i6NbJlDRYtlKyrmCEAKWVH6JxTKVgioUEyiDSni/aWVrj+cq3CeCGGJZbY8qL0gkmXIqvzySMdQZZWfvFiWFxreqrWrYgusCJPGYJK+dJiVd7zEur+9TnhCE9/seJNQKEQoFCIcDtOuqjwTfxU1oKL4NbSgj2sDv+bRRx9lyZIl1liKQjabpWW7MTRMbUCaMD07hV26dmH9+vVs3LgRwyh9cYXbIux46q5cv/WV1NfX1/5P8Nhi8IRyC+JTD5/PYGcXRjxN46dnFa/7NYNs/iwymbDf5pmmKNZJVD4MOcQHlX7U6rKOfXXyOTeEQvZbciEksYYUQ/0hAhH71aMQkkAsQ2Yo4BxnSaVY2pVCs8TSQC/LvKmds6QpmOSDvgT6QIp95u5MIBAglUrR2dlJKpUilUohEglWDWzAyOTQcgp3Ni7jkEMO4dRTT0VRFAzDIBKJ8JNbfkrf0l6mhCfjH+tn1113ZezYsbS3tzN27FhaWlq4suc6wq0Rbt3lN57HewTgCeUWwq4Pn4+hSJIfWoViAxNayh6V+DWD/oGw/ZPzmIZAWxF0tfFv1DCnJ9376ghJMOR+1ugW+F5OIJYh0+8+JxkwibW4ZN0I8GkG05rcyqpJtv/GTvQ88R5z5szhhz/8ob1Vsee58+xPOumkYW32knsNa+Ox5eCdUW4hvHrEzwnogp7HFhLZdiJatFxcBMm0H3/A2QsMIAf95JqGqRspQE/4XLfuQkgmN/YN26Rs2sSNtsVuy43qImnXepBCkUyY2O06jgJMaex1LHFWuJlR30DrgVtxzTXXkE7bn6UKIYYVt/+VjceWgyeUWxAvHXExyaXraZgzs+J6Ml3a2vkDur2Ale1Ic40OYqlAtt3yJusJn215MlU1mZnvFTOpsc+2RJkiJGNi1gpwfFN/RWuEAkJQbOLV0jGAYhN+pKgm4ydYq8T6UNq+UJGQTG4sb+tQayURJHLW32jS8bvS2dnJ7bffbjOah4c9nlBuQZimiV+rrDRul4rnC+iVqzRTIKtKpuUazcr/+ypk2ypDbvSkr2J8n2Ywo62y++CEhn58ZSmDqpC0RivDcsY1DlR4qBVFEqvqdNjcPoBaFqyu+QzGja/sa1MfSlcUIfYpJpMaKgtYWKvK0pxNBMlc6bWHJzTQdsBMzjzzTK6//no8PD4KnlBuQVx00UUYhsHWk2cAVqWflENHRV+hN7UhkA5tHHINljBJnyQ7xj4uUU9qSFMQ9OeYNsa+JcK4+gECqoFPMWmJ2scujm0YxK8ZqIokGsrY2jS1DaL5DfwBnbEd9hV86oIZ1PwKdUJ9v61NuVimcrV/n21+8Bk6jtyOM888kxtuuMF2DA+Pcjyh3EKYP38+v/jFL7jwwgt56kc3IBGks+59bDAEMu5uI6I6uRb3s82IkmNyS6+rzbi6AZoitTGG5Yxv7Kcl5l4KrbF1iPEd7r1uYoGMY03JAqYU5Ez7rCRFU5lxxlwmHrcz3/72t5k3bx4bNmxwHc9jdOMJ5RZAZ2cnX/7yl9l111350Y9+xD3/eBhpClSHkmIARlbFTKtW3KMDaiiHOjGJ5nfOhGmujzN9fKdt2E05dgUiqon6MrRF7AvrFghpOSJ+d4+6phg0hZOuaYVZU2UoG3S0EUIw/lv7037SZ7jxNzcyYcIETj/99Io4SA+PAp5QbsZ0dXVxzjnnMHXqVAYGBrj77ru56e/3cmnmFcA6M7RrCqanNXKFlaQqkQ5FKPz11hbYH8mi2fXwBia3W9ttTZi2PWOA4vWwlrUtLgFUVAOvD9h7nCO+LPVB6zG7QhYAPtUo5m87rRgzhlY8l9QdbOJZPwOZEM1H78G0W8/is186jt/85jfcc889tvYeoxsvjnIzpbOzk913353e3l6+973vMW/ePJqamvjdO7dV2PlUA4EsBlyDda5YgSKtzMCyVEYtVrlq84dyCAG5stqS41orzwkLgmiU1WBUqwQ0qOYQSNJGacsfUnMEy3rLNAaTKMKkL12K+4z5MxUryZCWQ0GS0kvj+NVSCbPia80X0S2QMjQyus/VZjATIJEtnV0q0QDLjt2KQzYewhlnnMHSpUs5+uijmT17Nj7fMMcbHqMCIaVdEEglg4OD1NfXMzAwQF1d3Scxr1HNypUrOeCAA0gmk7z44ovF1gHn3HUtT8a6bJ9jmAo5XSWX8FkVyu2QAmEIfPUZtJCD8yarkU34mNzeTXO9/XmiKZWi+DhtybOmRkr3ua4y49kA3akI9YE0IZ/9fLKGSjLnx6/qaA5HDaow8SkmSd1H1rD/7lfyNv3pIKmcvfjN3RBCvLmGP/7xj8TjcYLBIF/96lf53e9+58VEjkD+E13ztt6bGel0miOPPBLTNHnppZeKIvnNP13uKJJg9YPJxv3OIgkgJP7GtKNIAmh+na0nrHcUSbBEpzUQdz239Cs6IS3nKJJg1YHsiA46iiRYq8iQlnMUSbBWuDlTdRRJsMQ9Z6qOIgnwXHuK1J6T2bhxIy+88AJf//rXufnmm1m2bJnjczxGB55QbmK+eM/P+OLdVjvTRYsWcc011/DWW28RjUZpb28H4IS7f87rjQOu42RzmuM5Y4EJ47tpb+l3tZnR3E1j1N173exPElJzNPud7Uwp0IRRsU2vRhUmYS1LWHV33qiKWRE/aYchBX4XUQbI6Gqxl7cTz6ZNvvD73zBnzhw+97nPAV4aood3RrlJOeK+n7K+PkduIMnUz+7BisdeKeYaL1q0iCeffJKbUq+xrj6HwMqKqW5zAJDJaUgpUDSTQH2GzEBtYYzJk7rw50WirW6IzsFYjc3WrRuLK0CfYtg6S1oCieJ533BebsDK7nHQmcJlRUjCWpakXhvzaJSdqypCuqYpKkIS0HQyeu3bOmuoSEBRTEL+HCmb0KrkhiimrrBIMdnnuivYbeUGGhoavK6JHp5QfpJIKVm8eDEtLS188flf0y0TdP7hFTY+9CoA076yHy/86i4uueQSmpqauCr5MoOxkjDYiWVBJIs2qiTQkCHTXxLLqVM6KzJjfIrJ2PpB1g+UzmVmj9lQsb0VWK0NslViqVVtt1sDcboy0Ypr1WKmm0rN1rl6HMUmu9ww7a66I5AEtRzpModOQSSL9xKSsD9HskwsE+tiyDJRXq5neeH66/jBvO+hKN7Ga7TjCeUnyJ133slXv/pVAAJjG9AHU5i6wZgjdqX92N3R6sMct+Banr/2Wub87XzSvtrVU7lYprMadss1oUgCjWkyfVbhDLsCt6qQdDQMsq7fEkunM8BysWwL1sZAqkLSFhyiM22tUJ1WfOViWS2SBaJaptjQS3foPz7cqrKajOHQL6hsKx9fG6uosQmg9/RiZrLcnR7g0o98N4+RiieUnyD77bdf8Wc1GqRhzta0Hb0bvqbSiizlk+z+9/PQXaJSBAXxcxYMISDYlGZCk30qIFiruHENAzQHnVsmgCWWjX7n0msCS0TXp9w9h7qpEBzmHDGqZejPOre3gI8mlkEtx0DGvXxbJJClc1mL7WNK2JpDMul+XusxOvD2FJ8g48aNo7Ozk9jsCaSWd9G83+wKkSzi0N6hnOywTcMkM8Z0uaf6CcmUWA91PvsAcLBWXhPCfUS1DPbVfq2rfdnwsCI43ONgxVyODQ26inLUl6HOb19NqICmmDSHhhe5+gn2TjK1vg6tpYnBfz7PhGvPQ9eHn7vHyMUTyk+YMWPG0PXaUoJt9ay947maxwtpiZpD90BJSSQt54x9mbPp5QUsbBRFCMnkaG/F73bjjA/1F3+PaLXe6YJIFnASw48qkgXaQ4O2nu6Ir1RQI+a3r2VZ/nKbQw4rYVnKla+fMFAzjhCCpi8eR279BtZdcB0TLzuHoSH39EuPkYsnlJuAQCDAxWf9mKGFyzGzJQGpzt22xLL0AZZArmol6fcZFe0YNMVkWmtllZ+AqlfYqIrJpDKRBGvLWy5MqjAZVyaSBaoFpVwkSzaVv38UkbQTvLbgUEXWj51NzFc5b0Ht37EplKywkVKQzlWeOtWNH0SUn+UKSXCXKYy95CzUuiidv7yZra46l/Ve8YxRiSeUm4iDDjoIM6sTeHWD7Ye7gHUWKZGyViQL+DWrK6Ff1ZnSYl95J6AaRaGZELE/t4yUiWVHyH5LGlazxUBzp4ITAUUviuVHCUFUhCSo2IvpmGC8aGO3ogVrK164l9PfsTGYQhUSKQWZnP3RfF3HEIpmWimfYWtFr9bHaPvRyWhjmum87GZ2ufUyFi9dOvyL8hhReEK5iZg6dSoAn/FNwT9MyTBNMV0Dt8EqkDFumKB0v2oQUN37YUe0jGs2DZS2yAM5Z6dLQLFSDgMOAlick6J/BBuDsINIFtAU56IdBRqCKUeRLFA/bhAZqjz2UCIhQjtujczkkLrBUa/+mQVvvOo6jsfIwhPKTcSiRYsAeLpjiWM1nQISYfWpdknj86kG8WyAnE1AeoGgmmN8ZMDVCeJXDMYGBhAuEYwF0fa7BJyrwiSmZVzFK6jmaA3Ehw1cby4LcrdDIKnzpWsKZlRjSkHAJTNHVUyCvhzBYO2XSWbJckI7bo1/UgetLYOc03Wf6708RhaeUG4iPvzwQwAi4xsAHFdx5cHkfgdvuF8zCPutD3ciZ9+uNqJlaQlYYUCKQ4fEgKITUy3RVpC2YqmbSrEyUNRh9akpZlFE/YphG6MZ0bI0+635xLT0R1rF2gmqgiz27QmrOUfxLoQTNYWShPz2q+qCiDY1xglXdZlUG2Lk1nYSLkvbvOuJv7rO2WPk4AnlJuLuR+5DiwbQwlbaXtSfJaRVfoDt+uFUB48HfHpNUYl0VQpfTEvTWJWXXe1RDik5ompli4bqvJicqZIxKwM8I2q2otakppj4RKVY+USlWMa0NA2+yvlE1EzFOHYEFb1GCKtfR1DN1Wzlq+W+IZgiGqh8rdVRBg0NCaLR0ko/euBeGD19bLjij/S/vhI9neU3ynxu+NsfXefsMTLwAs43AT/9w+X849EnaN9nWsX1sC+HIiSJnN9WJKFUfzJnqAR9OQI2q8y07kMiCGk56n0pYj77HjUKVnPGiJolqNiLlILERJA1VXKm/dslrOZIGz58iuG4RfYJAz3/vewUt2k5iqRtzneBgKIjkGRMzbFQRkC1mqulDauTpN36ORawHFeD6SA+zcBnE45VF0uhKJLBwRBtOzbR9NOjWXHTs7zzw3tBEdTvMIH+L+9F3/1xzj/udMc5e2z5eEL5CXP+/Zdx9z/uxcwZbPXNPWoeD2o6ad2H7pJ5oqkmPtW0/XAXyOgakyK9xWrgdgghadKcM27AOv9TwFEkC0S0DOYwDie/YhDR7EW7QFDJkcRZKAvjCOGc5liw0YTJYM45Oyfit2plJlyEORpJ0xRKYqgCdp9G425TSXy4kaF319H5+FssvuhBHr61mdBDF3D2kRe7zttjy8Xben+CnHX/hfw7/AGr/vY2E4/ejkBzpMYmmfOhCOlaf1FT3EUSYHpdF1FfBk2420XU7LBeZ4EkrLoLXL0vTb0v5WozPtTHGP+gq42CHHYcoGZ7X40qzGHLrvkVg6g/Q0PA+X4NgRSN0QR1eYebUATR6W2MPXInJp34afShNG29S+ma9A6/eOTsYeftsWXiCeUnxMl//TEftPXw9jVWNs7kz+9YY5PI+Ys9XhQhHcVwuNqMAHX+dGkcB1Fp8VmZJqowHbfeBTRhOq4GQ2o2X2XcqDl7LDA+2EdYyeITJnWa/dZbERIEruOU4+QA0hSzGDPa4LcXQb+qF+NBncZpCiaLj/ltvkw6n3ybQL2fmTuHMXImL7/1HHsesz0XXHABt956KwMD7uFaHlsOnlB+QqxvibPu6fdZ9+T7bDtvH8IttRktRtVW0m5LPNxKEmB24/rKcYTEJ3SklLz7wPs88f1/0Wj0VjhGrKBvd7FUMYlWiVzhXLGAppg0+SuLbEwM9VSsWlVM6rVKAasWf7txqlGQNVk/1kq8fD4GTYFK0Q2oeo2jqiVUe69qL3tjsMoBNXUMmYEsvz/2KX572OM8ct6rvP/6cm666Sa+9a1vMW3aNG677TY+QrcVj80c74zyE+DzT8wjm0ux6Mp/MXb/GXR8ZiugsgrOUNY+rMenGuTypcIsR4472zettd229y0f4O/ff4n+5dYqx6dtW2OjCElIzZEynEsXKUhiWpohPUhEzdiquSIkzf4EPdkIk8PdqHb56EgatCT9ethxhVw+jhMCSTDvTKoWydI4Js2BBD2ZSPE51ajCpDWUoCtl2bTaCKcqTJpCCXrzNsee2cba3eaw9J/rUP0q2x01iTFb1aPEg5wx+XLOPfdcvvGNb7Bx40Z++MMfOr4Gj80fTyg/Zj73zHfJhU36X+5Ej2fY+v/2rGgtoAhJfzqIU8k0Aa4OmXJ2bF6D1HXWvNZNsidNsjdNqjdD99J+Vjy/rmh36n2fRvM71Gn8CKVyBXJYp4wQkumRjY7ee2scaPUN0aPbVFAqG6fFH6c762KDFRakuziThJC0BBPEHeJMwRLUMeF4fkwH8UbSHEqwVbQTVUim7NXGlL3aKmzMaJrr153NH/5wD+PGjePcc8/lhBNOYNKkSY739ti88YTyY+b+fa7icy+eRfdrq/HVBQi1V7ZgWDNQD0A04Jyi19lvPaetwal6jWTXltVk41nuO/mfdC/pt7Xa60sTOOS7M5Dh2m1/OSE1S8pw9gTX5bfN/XrEUVYFlqgMtwYWQtLii9OdcxbClOkjomVI6M4iJxGo+RW6s9RLor60q1gWso7cnGAtgQS9uQitfucGbG3Rbn79zkGcccZtXHrppTz//POeUG7BeEL5MaOqKg/NvY7olxpo//Q0hGIJh5SwdrC+aBfP+GvEUgIb+0vC2tkfY0zDUIX0CCHZtXkVAI/+6OWiSDZPr6d9u2YC9QGiY0Ls/ZkIjR2F3OwsSdNZCAXW2WPSQSwL92/UEvTp4Zo4RYGkUbPO81Rh5otnVAumLDsDtLKA7OId3Y4BCphUrtCdinUUiPoyNmIpKlakulTzMaGVlYlK552CrmyMFn/V/w+g3Z934igmlz97NAAvvfQSX/nKV4Z9LR6bJ55QfkIcOfcw/v7qU4CVTrdusLYaeDzjJxLIFj945SJJ2bXW+iHyelsUydfuWMLKl6wSYEdf92kmzhlX/IgrSBoDlRWDwoq7WEKtWAoBdWqlE6Yxf85YEKvC2WM5qpAYEopiaeOJb/bF6c1FS6InrZVkOREtQ9KoDMY3bcRV4FRiuMxG5Buf5Z9ht223epcbRfsmm66T3dkYrf7SSr8oknnGTA7hCyh0dTm3GvbY/PG83p8QBx10EINLu9CTWVuRLJDIWMJU3n2wmq4BS0ALcYIfPreO569+E4Cv3X8Ik+eMRclLiCZMxleJZIGwkh32TDKsZhFCoghZI5IFGrQkCiYKZo1IFlCFBCHzHnj7bW2TL46KibQRyer5gL1IgnMuezkRLVMsAed2tlkQdzuRLGDkP0bVBUCklNx9/nsEIiozv7nO7qkeWwieUH5CtLa2ApCJD++YMaWge9D5zA6gYyDADo1r2fBOLw+d9QIAh1yyO01TSyIcUrJ0+PtdxwkNExIEEFGyjiJZoEFL0aC52xTOLd3HSaJL9zYXw/UBh/wW3KbdbjkRLWO73S9Hl8qwGUe92QgBRWeMr/IM+bVHOnnn2W6Ov2gWrZPcz4U9Nm88ofyEiESskJKNfb5hi9kmsz7Xkmo7JgLce9zPiKxo4+/ffwmAjh1bmHnIxKJNWM0yOdg9bI1GcBcvVZiE1YxtwHXFGMIsBnDbYUiFtOkjafpda2tmpUaTLzFs5k1kmBJuWUOjJxN2Tb3UTdXqXe6i3aYUpHQfacN5nHotyaRAN2rZnHvXpvjLL5aw02FtbHdAK2KV1xt8S8YTyk+IZ1950frBMK2t4TBiGQzkbNvM7puMct+JPwBAvhpjaIO1JTzwJ7sWw45iapoJAavVQ0jJDpvGGFRytqKjCZOQUlq9OQlh+XUnm4wsCU3K4Wy0fOXW6EsOm1oZVrO2RTjSho+ejLWC68uEyNiIXNbUinngJgK7mHBTimIM68ZUjKRNa8xmX5xtI2ut+eT/1n3r0vzm5DcIxTSOPW8rfB/uwHcOvsn1tXhs3nhC+Qlw+6P3c+1jdxPoaCIwvhngI4llwJ9DK6sO9LlsCzd/fR4Ab775JpdffjkAu+67E01TrC13g5akI9BfMU5QyeET7qITEHqN6ASVyi2uTxg1gmp3xll9ze4MsHqlZ0iFXNWWuzp7x45Q1TY8qfvpy1RWXh/IhkiViVzG8NV4vWWVWBpmSSQLdKejFc9r9w+wdbgyC2r9oh5+eeR8cmmD027ZibbeAzj9s1cO+zo8Nm88ofyYueovt3PVwAv0v7CYlkN3qgk2t2uYVU7Ap+P3GXxTTuDyL58GgK7rfOtb32LmzJl8+ctfZmBtnNlr5wLQ5lB0IqDow26f/UIvbnn9DsKqlYmlwH7+QpSK/uakSlbWrugyUiObF0tdqo7nkpGPcB5ZeF3xXICBrH21oKH86jGl+xyrBRXOK3VTKebcV9ObCTOQDTIh0Mu00MaaxzW/Qi5jsO+XO5hpfIFvHf7TYefvsfnjCeXHyEX3XcudwaX0PrUIBDQfsF2NjTKcexbwaTp1QctrahgGZ5xxBq+++irr1q1j7733ZunSpRww+zD22Hik6zh+oRMYxnnjEwY+4S6qpZAZZ5EXQpI2tZpVYjlZqbE+U+96ZhlRM67B72B9CRhSYcglkBygNx35aGO5lG8DaPYnmBi0b+I2cVaEYEQll8jxlcPOch3HY8vBE8qPkWebV2KuWcvGB1+hYc7WaPW1nk9VsWpLurFLy2o2TPk3hmHwpS99iZtuss67+vr62H///Wlra2Pq1KncetXw1bZ9wqjZUtfamMUcdDsK4udmkzQC5KTmatOdjdKXi7Ah4xwutT5TT8ZUHYPfAeJ6AFMKx4LAABlDYzAbJJVzD2D3KQYxv3t6ZntwkK5cbYxrLmNyx4XLSA4aTAt/1nUMjy0LL+D8Y+KU+y6mb9Uq3v7J3wiMbWDcifvV2Kj5cmBCSHz5quXV7Nm2gpZgHAkcccq+PH7fi8XHfvKTn7DVVlvx1FNPsd122xGPx9lx6FLeiP3EdW5ujbpKCExZW9UnZ2pl8YsCU4oaB07CCBS30tJhnM5MrNhWImEEWJduoCPYX2GzJt1IMp+ZkzVVpAzY5pgXtvB+1aBepG2334VVYs5UkTlB2Ff7ZVFoh+tXdOr8kkGbcWbXW2eSSTNAZ04h/s4a3nt5gLVLk6xZkqBrdYZ5Z5/J5b+8qua5Hlsu3oryY+ArD5zP4vpe3v3Z34nNGsvO1x5P47jKbaFaVjMRQFFkTfOwuWM/oCVYyide+OIiJs4eA8DZZ5/NhRdeCMC2227Ld7/7XT744ANmz9yJT5vXDTvHsOK+arKoXA1mK0TSnrgRrDlvlHlBLae6907K9LE63VhxLVmVvpiTCkN6pXgNVrXM9SkGjVWFeNNV3mrdVEhkK1eo0apwI59i1BT03aFhTUVR4Vcf6eQXxy/iiVvX0teZZc4xY7j//ge46opr0TRvDTKS8ITyf8wxD/2ItY0pjHSO3ECKjs9uhxbyoyqy2NBKqxLJAkKUxPKAce9T76/cSvqCKnWTfZz/9Je44oorUJTS/7558+YxODjID3/4Q2bP2p7PRty34YqQVpm0YSiE7GRNzTE4uzysx+m8sfy5a9MNtjbZMk/48mSLrY0hBYN5sRzIhWw96qowiznZKd1nu/03pEI8X9ou5svYlnpThUljfpxdmlYVV5ypwSwv3rqUB89/ndmHdPCD5w/mR3dux5Vn/osjDz/adt4eWzaeUP6P6YtazhIt5EerCxJfVsrxVQU8ucfZ4OIEUQQ8s9f3CGu1W8NQnY/0YA5aE9y08KCKx4aGhkin0/zxj5ZAThg/iZji3i9cIImq7jZgCeFwGSymVBjUQ642hlRYk2p0dd4sT7awLNlKzsXGlIL+XNh1HEVI6v1p1zJvphS0BYdc0zgVIfls+yICio5pSF6560N+fdCTPP3r99jti1M4+pKdkcLH0a0PMWG8Vx1opOIJ5f8YkRcUoSo07TaZngXLi9cf3evnNDU18a99LnR87rP7/JSmpibOm3FPzePhpgCJ3vyqNJrg92/vSy6X47LLLmOXXXZhm2224dFHHy0bTxJT3GMRM6Y2bBaMQe05ZDV9etjVww3w3lA7/Vl3MV0+1MTKoUZXm6Tutw3+LkeXyrChV1s1bBy2rcbejUvxCx1jxXpu+/oLPP7Lt9n+8Amc9Y/PcPD3t0Xzqfxoxj00NTa5juOxZeMJ5f+Yp/f7FSK/itEiAfShNALB0/v9imDQ2jIGg0FUm5WOIqmwEVWFMVSfwtDGNJMD1ip15ZI+Jm/TwLnnnst3v/tdXnvtNebOnYtpmtz5vtXhUQAxJWW7asqYGoVzSCexLK8naSeWEujTI/lCacJRLN8dbM+3jsXW2SKBFUNNxXFWxxupPiMFKlrZOqUVlm/HO6ID2OUozmzYWBw9YrN6B/h04/soSPrWp7nihIWkB7OcePscDjt/e+ra8oJvlv6feYxcPKH8HyOE4JkDrsCvKyRWdONvDNP58CIuvPBCfvzjH3Prrbeyfv16XjrkUvx6SQhET4pHPvV9TNMSIyklp7bewKoF/bx8+wfc8+35vPHXlex/fCsCeOXRbi794iL8QYUXXniByy+/nGAwSCqV4p5leyDKRM0Sy8ottuVMqRSiarG0K7pbLpYlkaykXCylFLw72F4zzkA2WMyEMaVg5VDtimx1vKFs6yxs+31Xi6XdmWVHdLBidWm30oxo2YpMKVWU3FaP3rCSQFjlR/fswKSdy+aZ1bhwhwdrxvIYeQj5ETofDQ4OUl9fz8DAAHV1zjFvHpXM/MreLHvgFcyMTnt7O36/n9WrVyOlZO+99+b000/nt+pCPrj9adY99jpSShRFoampiUQiQSpV2jaP266R9g6F5g4/y9+Ks/S1IXY6sIlTrt4Kn1/lqJa/AfBEz6GO8xkyg0hEjce5nJy0OtwM593WpUq/7l4RR0GyZGiMq02dL8PKuPtWe1xkwLUoBVBsMuZWMm1DIoYAptd3O9okDR+aMNmzflnx2k8PXsB2+zZz3I+nA9Cl1yGHQly4559d5+SxefOf6JonlB8zuq5jmiZ+v7Ua6u3t5aGHHuJPf/oTTz/9NHV1daTTaX7+c+v88vTTTyebdQ8I3+nAJjqmhzj4G+OI1FsCYkiFgJJDw/ksUQhJjx7FcDlLFEKS+QhhQH25SLEOoxPL4i2uAedgCfJAxv3cUhMmDYGUa4E2TTHRhOkqlIWWu25l0/yKznaRNcXf0wmds3d9ka9dOpPdj7ZWxkZPM2fOecB1zh6bP/+JrnnBXh8z1fF0TU1NnHjiiZx44oksWLCAG264ga9//esccMABPPjgg2SzWY455hjmzJnD2LFjGUoM8Hbk9zxzy0reeaaLE38+jU9/vrKZlS5V+owwGNCmDdrmaQshCQidDl8/G3L15Gzyr8FKc/SrOnEj4CKEgnotxaBNzGRxTqZCLN9ywak1gyEVkrrfsQ0EWPGMLUGrI2LWIf86oBh0hPoBWJ2yX536hUFb0MqDH9KDth7zsJplq/AGDBSU/GnpYLcVxdA41golMtZN4Mz977C9h8fIxRPKTcjuu+/O7rvvXvz90EMPZdy4cSiKwrx583j51X8xP3ork2ngxF/X88ErfczeudZx0GeUtsCdeh2t2hBBUZnTHSgTz3bfABv1GJmqcmfl549RNUPCDNSs0KQURVGrU9MMGcEaB07WVOnPWXOK+jIkdD96Vf60LhVS+TNHVZEYJjViGVD1ip7c5e19C4TUHO3BUiGQoJIjXXW0EFRytAZKgfsxLV2cX/m18iIXJgIViT9kzfutZ3qY1ro7Zx7mlUsbjXjOnM2IQCDAFVdcwf333891113H/OglxccUVbDVnk3kAuGKajzVogDQpccq+uHYeavHaEMV2TmKMGucORElU3GtXCQLxNR0xQo2Y2g1IhTRsvjKChHnTLUokgVUpbISUUjLVYgkWFvn8nCeqJatEEmA1kC8ogJ6WM1WiGSB8qIfjb6EbSUgiaBhTIADvjGe1x8Z4IDJZ9XYeIwOPKHczDjhhBP43ve+x7x581j9tn3JtIIIpkw/Q6Z9aEqPYbWSUIWJ3yH0p1mLE1VTqMJwDA8KK1kCQrcVyQJRNUNQyZEyfAw4BJ2HtSx+1SBraqQdOiuq+d48ES1Lg98+/lMTJqqQ1PnStAbs2/c2+xPEtAwxLUOzP2E/HzVbrKQ0yaESkEnhNQsG+hKoqnucqMfIxRPKzZDLLruMMVNDPHLlUkebrlyMuOlWVkzgW/aVGgGMD+isWJzCNK2VWaOaHLZIRkxND5vB47cp/FtNnS9NTHMfpyWQoDXk3C8b8imKLs2+ABp8ScYH7ZuqFQipuYoOinY8e+c6nr51Nddccw0zZ850tfUYuXhCuRmiaRrf+Or/8cGCPt5/qXa1M6CH6NZjxA2nQGfBzHVf5fOHfIfginnFq4YhueDLH3D24Ut4/E+lEBm37B2f0FExiLqkQ+byhXfdKpL7FR2/MGgOJItOFTsaAkn8ik5Qda6bWe9L5YPXnT3qMTVNWM0QcxH4ei1Fiy/u6JBa8nIfD1yylO9973t85zvfcRzHY+TjCeVmys9/ciWtY1pYOr9yVdSnR0jmV5Jp08dQzVZXsGvfmRx94JcAOOyAr9G84WdsWJXhlgvXsHKxJRx6TlKIDAspOerV2hWaX+ioZeFGjZr9NrbcgzzGYYXmK1ttNvjtxastVBJQn2J/FNDoTxbPIKW0j/es11LFFhEhNWsr4I1aknC+LqddGwoQLPy9n5133plf/epXtnPxGD14QrmZIoSgLlbPBGPv4rWeXLTGeZORGgP5wO94T5boc4exdNEqrr32Wt5++22klOwz53D+dlkH/7jbWp3uNDfGHy5dx33XdxbHCQidBrVSCJWqmMygyNKsVQphdZuHsJqlPTBQ+TybquoTIpVfAB3hgZpUwurCIM3+RO1YVcGVjVqypop7QMnRWNZvvNkXrxnHlErZaxFsve5s/r1gISeeeGJFlSaP0YkXHrQZk81mGTd2Ip/VzuUPqfMc86hzUmX+Y/08/ONXMHL/rHhs8uTJnHrqqfzoh+ex3bbbc+WVV2MYlrpM3Kpy6+4XBk1qnF4jSkDYb339QqfVN0hXro6M1Gyr8wQVnXHBftamGwipOdt2uGE1x+RoDyvizYyP9Nu2nlCFSUTLktD9tAbijt0kreLBkiZfwtHGlx+/1TfkeJYqpQAhOEi9npsf/T26rnP88cfb2nqMLjyh3Ezp7+9n9erVzJw5k5kzZvHt9ddx9cazbG3fuG8ZT/3idSZPmsSKFSuL1w888EDa29v56U9/SiaTKRZv2Lg6S/skP3sc3FAzliZMmtUhEi6OIg2Dif4e3k/X5nAX8AmDnWKreD/pbBNQDD7VtILuXNTRRhEm40N9riXVwDo+GK4t76RAj22js3K2Dqxj5pRZvPPOOwBeJpoH4G29N0tWrFjB/vvvTygUYu5cq7tix9hxNXbJvgz/+NlrPHXJ6+yy8y68/vobnHDCCYRC1rnlU089xeuvv843vvEN5s2bx5e+3cIPrmrn+FOb2Lg6y5svDGGXwaoK6eoECSo6qjCYGVzvaNPqszKEZkXWOdq0+IfwKwZjA87OnSaf5ZV3a3YWUbMIIV2cWxBVMmjCqOhTbodP6Dz8+g4IIQiHw7Z/H4/Rh5frvZmxbt069thjDzRN4y9/+QuzZ8/mpZde4oMPPmDNmjX87b07iHenGepM0vleP9KQ1NfXM3/+fLbeemsAEokE559/Ptdccw0A0WiUeDzOuMk+zr++g1BE4eJvb+DDd1J0TAmw7+ea2OfoRlo6KoPAB22Ep/psTyJYnOqouDbGN1hR1k0ieCdRKfTVYTkSq5FYOU2+RI2rJmPWrgijVX107EKZqqsnJW1WzLuEliOE5N6b+/jNL3qYNWsWCxYsIBarbSTmseXjFcXYgrn++us588wzOeOMM5g4cSLXXHMN69evRwhBS0sLY8eOZTDcSbglyLjtmrnju4/S1tZmO9bf/vY3jj/+eLbeSeVzJzVw4Sml1Z2qgVG1SFMUME0IRxV+/cQsGtt8DJWJpZ1TBmBxqqMYqtPmG7C1AXg7MT7/k6TVbx8ruS4vlgJLKO0oF8tqkSxeVzMUPD1Old5LYinZNby8eH2w3+CuG/v46x8HaG+bwNNPP8306dPtX5THFosnlFswK1as4LTTTuPFF19E13WOOuoozj77bHbYYQd8vpLH+4cLvshlu9897HhX3D2Ri09bS32Tyje+34qUoPkFG1bnWPJminffyLBhZe129Df/2oYx46wV5pAZdHTuFFiSHouUgjE+5230u0lLUFt8zgHlhVWlk0gWyJiao0gWiCgZ6lT3Cu8p6WeX0HLbxzasyfGDr6/DzIzhueeeY+rUqa5jeWxZeEI5AjAMyzHx36bNPfT6nvT2d3L2CauYMTvIxbeMQ5RVqJUI1vdoLH0riTQlW+8aJRKrvWehTqUTCpI+I+LqdBHCyihytUE6etPL0YR9g7YCEivkZ7y/13nOQjLO1+vaW6inU2feCWuRehvPP/88kydPdp2Xx5bDf6JrnjNnM0VV1f9JbvGRO71MXWgK3/lZGwv+meDZh0tngyYKcTNIrFFj533q2GW/ehuRlPiETljJ2Ib5AKhI6pQUk3zdjjnjVpm3HOP9vc42SAKKTp2adu1lY6C4imlBJAHHrBtNGIzzWSLqVDNTAGPbBNfd3Y5UN7DvvvuycuVKW1uPkY0nlKOAY/Z4mlkzd2P73UP882Fra2xIxTUECCzhKhe1kJKtEUsNsyIFcoKvh0CVh1rJ18Is0OHvqwkKV4SseF5MTVf02QZLAMtrZGakViNyUla2z92Qq68JCfIrOmN9/cXfEzJgu8r1Yb321naN6+7qwDB09ttvP1atWlVj6zGy8YRylPCFfR5kq+2DrHg/iy5VknJ4kbSLSwwp2WJao0/Y54CP03qL6YGKkLaFhNt9A8VQHUWYtjZRNVMMDpcITJu3a7ZKLO1sNubqSMvS+W6bVutwSkp/Re3N6i+Etg6VX9+jksmk2G+//Vi7dm3NGB4jF08oRxGtsV3ZsDpH75D7lt5JJAsElRx+oRMRzs6Udq2fBjVpK4AFxvgGiahpxzJwABE1g08xXFtTZKWGIRXXc83uXIykGSCiOM85Jf3kpIqCtG2p0T5O5bo/+xga6ueUU05xHMdj5OEJ5Sjiq8ddBEDnGueg6149Qr/h3jRMSlGxlXaiw+de5gxgqr9rWJtmLU6TQ0GOAkElN2wwuYrhGkhv2UhX4W4fr3HaBQH+/ve/c9ppp9Hf3+86nsfIwBPKUcSkSZMA6HltP0ebLr2ODbkGNuQaHG3MfP9uHeeVaeFss0l1DgVqVBM0qQkm+O0L54IV4hNRMrRoznUj69Uk9WqKqJpxrJs51tfPWH8/Wak6pjHGRIaYkkbFtK0KD5CWGnscWs/pF7bzhz/cxrHHHus4L4+RgyeUo4j29nZ8Ph/jWrYjvOoHNY+vLxPHbj3Gmmxtr+3y80BdKuRsxNInjKIjpkFN0qrWxla2qEPU551AHVo/U/21rRiiarpildhuE8zeqCWIqaXtdMRmVTnB30NrWXynXXGReiVVsS1XkTXOpJT0kZUaQgiO+Goj51zVyjPPPMO//vWvmvE8RhaeUI4iFEVh7NixrF27lkP2+yJjun5ZfGxNtolBo7K2Zb8RYWW2tfi7KZWaYrmGVMiV1VbxC71GYGJqmvYyB8oYbZBo1VlhqzbEVoFS7nidmqppkFanphjnL23nW7S4rTA2lJVUmxLosq2jWV6urlFJErIJqC936CTz55fl7PmZKFO3CXLRRRfVPNdjZOFVDxplNDY20tdnic3eexzMm4uaeEY717ZJGcCQEWRZpo0p/i7HgHNDCiQaUcU5/jGsZBin9WHifL7ZqCaZHVzD2lyTozMpomSYGOghbfoczxIDQqdJS9CkxSsaqFXMGYW06Wec1ucY1wkUe4Xb9QsXQvClbzfz89P/yfPPP8+nP/3p4mPJZJLnnnuON998kyVLlpDL5Tj99NPZc889He/lsfniCeUoo1woAf6w/jZax9mLZIHuXJRBI8j2odXORsLyGrt5wn1CtwxdSEk/MTVFqqqVbjmNagJTFfTpzuXZfMJgiq+LTsM548JAMGQGaVKdHUUCSVI6/332/EyMKVsH+OJJhzP/2XeIRCLMnz+fiy++mPnz5xONRtlmm2145ZVX+POf/0wikahIRfXYMvC23qOIRCLBG2+8UUzDO/Of3yLZPMTqdKPjc+JGADO/+nJCFRKBxEQQd+gKaZ1tuoskWFt5BVmMw7SjEL5TXW29nO2Cq9CEUbHlryaqpDFQ6DEijjYDZhCJqDlOKM5FEXztkq1IJtJMmDCBpqYmDjvsMPr7+3nwwQfp6+tjwYIFHHjggXzqU5/yRHILxVtRjiLuuusuBgYGOOWUUzj9pZMw66wtsInCqnQTE4OVedFxI1BxJvlacgq7hCsLSNhlz9gzvEiWr/4EkrCSqSmHVtGLHEmLNkS3XlkGbYfQymJQvIpJh9bPOr2hwiampIszMlHoNqK0VHnoB8xQxetRhVmTwbNeb2DitoIL/7YDbzzdiyoCnHXEvWyzzTbFnPqFCxfy1FNP8cc//nHYv4HH5olXFGMUscsuuzB+/Hg6zqtHqvb/2ycGe4oCMmRTj1Ig2S3yYfF3p5VWeaEJc5jq5BLo1OttHytPswwrGVu5lUC3XgdIdgqtrKiFWbIRRbEsF8lqq9a8WPab9v3JodRMbZ1uvxKXuso3t3mu+PvJJ5/ME088wYcffoimeWuTzQWvKIZHDRs2bGDhwoUce+yxjiIJsCrdjMReJMESnFcS0wBnkQSrNBv8dyIJFEN2Ig4iCYW1qmTn0ApbkbRsrLVxnaNIWlYbjZirSIKVlukkkgBCM/j9+1ZTuBUrVnDnnXdy8skneyK5BeMJ5SjhmWeeAWD5lFscqwCBVVRiWXLMMCXVTIbM0LCl0NKmz3UcE+gy3KuHSwRp6XPt4Q0wyd/taiMR1Cu1LXkr5yPo1Otd0yUB1rkE4xfHSlor4bPOOovm5mbmzZs3zDM8Nmc8oRwlPPzwQ7TPiBBrCTA93Gm7GpSyVJZsY6bOtvyYJgxmhjYAMGCGHUuUFSrvONWxNBB050XSKQvGRNCbd7T0GhFbARNCFs8tLRv7t3Qy74yyWvDay/dG3dp+rc01OdbMLFyvbuVbTq6/npN3fJquri4eeughLrjgAqJRZw+9x+aPJ5SjgFcXvsJDj/+FmXs1F69NC1VmwpSLZIGubOVqLyByzAh2VlwbNMMVoiIoiWSBarHUUegxKoWjWiwNFPqqvNF9RqSizBpASGSrbMI1NtXiaa2oK8WyOk5ynd5YE2CuS5UNZccE5WKZHNIxDUl2w1hO+dSjANxyyy0AHHHEEXhs2XiHJiOc195cwPeu+Tx6xmSfkyZVPBZUsqRNP6YUjiuonKniUwzCSpZJgW5bmyEzRINqNQLTsA/ezkm1GCDe6xCOowizeKbpVJij3whTp6bwYTiuRPuNMPVKCk0Y+cDy2pAcBZmXOUFWqvQatSu+DXoDY7RBAiJHVmp06bXHBAomd/xsBc/8aQPBsI/jjjme9t4H2X777fnlL3/Jt7/9bdrbnVv2emwZeCvKEc6fV57Fi3et4TOnTaF+TGWozcRgLwq1IS/l9OYiZEzNUSQL5KTmKJIFTASZYfpqK8Ictvd2IdWyOsWxnAEzhADHjCOwxDItNVuRLLBRryMnVVuRBFjxdpxn/rSB/Y/cnh9+/1zeeOMNPve5zzFt2jSGhoY4/fTTXV+Lx5aBt6IcwUgpefiCTprGh5j7tYm2NgFFt6oB2bSBLaDk+9i4lVbTpYIUwtHrDJCWvuJZoRMq5rCOm+EanQGERM7VaQWQlIFhS8r5hEHCDJJL5Fj2dpL3Xxti5TsJ2iYF6d2Q5a1/9TNhyhgev/9VfD4fF154IUuWLGH58uVMmTKFmTNnDjtXj80fTyhHKIZh8Jvf/IbFby5n3k+/heZfUWMzqFsrs5CSQ8G+Z/YY/yDt/kESRhBTyRCyaVkbElnCSrYocHZimUO1Vp35oG3DRgw1TDRh0KjGUTAZtAnTCSlZwvlzyZxUbfO0w0qWaD6VMiB021Vs3Azajl+OX+gEhE42Y/Ljz75Dz/ocmk8wfacYCx7tIVKnccgR+3LzdX+uyLiZOXOmJ5AjDE8oRyhf//rXufPOOznppJO46sKbue+J37N64j3Fx/v1MEN6KVYyoOQQyIqtakegv6K1bMq0MnXK0wvDSqai8o6dWGalVlG7UhUmSFHhdPEJo5hNA1aNSUWY9JedZ1bfKyc1JKKiinpUyRSFFCyhFKLydQ2Y4WH7BQVErnimauQkfRtzHPHNMRw5bzKarzBvwYlbPeU6jsfIwDujHKHkcpagvPXWW1x55ZXUiwnM3GC1L+jJRStEsoBf0QmplshMCPba9t9Om37ieZGJKmnb8mRWfxtLMDNVIllAFaV2C9UiWSCmpGnOZ8o43UuXKpl80Yo6JV0hksXXhVGxEh5OJIGKykSKKlAUQes4P5FAYZ6Cseu+N+w4HiMDTyhHKPfccw9PP/00DQ0NXHDBBRx88MG8/coyurIxkobzOaFPGEwJddGoOQdnZ00fmjCGaQdhBW8bLlXQFWH1ArcTyQJhJUOjmnC9lyEV0tLn6tzxYXnuq8OS7KhuKfH8Q70YumTHufVomAQVnRn9P+fgfY8bdiyPkYEnlCMUIQT7778/Tz31FENDQ5x11lmc/9PzGej7f+2deXiU5b33P88zazYyhJCEJQkgmylCUBDRgoAErVKRU49tFUGgImrRI/B6FHi1r1hUjgsUq3Kwlla0KsdCXQlwNBEEAVlCUDCAhpCQlSSTZDL7c79/DJlkMktEgUhyf64r18XM/J77/g1X+HIvvyXyzXSc3oHSRv2KBH19xEsbgDJPPHZhDIqFbI1LM1DaqmBFa0yKJ2KzM4Buugbcbfw6KwiS9MHV1lvSumCw5hW8t6acKydZSEk3AQo/c/4311w5IeI4ko6FFMpOgKqquEaV4nF6KN5bGba4bhe93X/ZEi7XO1Ff579AaYiwhW0K8QnVdqGJptYTrgghPE1B8AqEFctkvRWD4kWghGxNAfgvj8yKmz6G0KFOFl0jerwBXSj35VgpO+FkypxkQGGk4R0u+9nwsP5KOiZSKDsBD2+fjnmwhiFaT12pb0sdSix1Ld7TUKgLIZYtUx+1M4VvW1PiDiwYURWiwG6Ju2tAJtBJd7cgG19GT/PyNtRCN0VfG+CTgCCxbH3DblA89DMEZiZ11dkCMm0UBAbFw86Pa0kdYKb/0FgmdN1En/R+IbyQdHSkUHYCvLG+kmd6o4q7sfmsr6VYWgzBZ5ICBWuLPjpJIZp7CQgQy2J3QsgUwpYB28XuhKCqQppQKXIn+l/7VqLB0tgyHKiHvjak4Ld8J1QYEgQKfledLeRRwsmjDnZtquHqmxKY3Gsb3bolBtlIOgcyPKiDs3Dvb0AHtioH9loXXfsEZpioiqCL3h4xxNvqjWKguYxwxSQEUKeZsXqjwwaLCxQqPXFnwnTC2AiFE65EehpqI34ng+KlextnjS50ES+JAAYYy8NWN3c5NZ6fX0RyqpG/PFlAVFTkmEtJx0auKDs4z17hi50sOeA7l+uZGbgqUhVBg9cc9mpG4EsDPNiYGiFjRqHWG9NmRk25O/6MmIZDYFI9nPbGRrwq0ilaxLRDAIvqIC5COwmAONVJH0N1kNdHDzaycMpRSo47eWD5JcTERL6QknR8pFB2AnSNRsq+qiYuJYq4pOaVUctta0OI80gBOFtctOQ3poawUdpMAxQoVLibK0jXeaNCCKFPJJuoDiOWLbfM1d7YkDUxW26j41RXSPlu+V66oRoVgcup8fp/lfLIrccwmlWeenMoC6d+FfG7SToHUig7AcuvXof1qJ2Efs1iFepsr2VtSYESIJJNuFukA2qoWNsQSYBKd3BBiTpvVEB9yZYi2US1N7ZVDcpgn2u0mAC/VQTxrVrUxrYSS+XMewH+fFXOwluO8t5rVfz2P1L4w+qrWPCrAxG/l6TzIIWyk2CoiSUxMdFXUafWgfVUcItWm2bCK3xZNaHyvgEO23v6OzLWeds+twvVD7uJphAkJUy4EtAiDlME3Mq3pFaLwStUdAi6hOnjHau6UBGoiCCR3LXdwcxbK4g1a6z8Vz9umDKSaVd/EuFbSTobUig7CQkJCVTsr+PtWZ/w0vj3eXt2bki7Rs0Utmp5EwWOFGo9ba8kvUKlOkLvbfAJqTFiho/vBjycSDbhFAbiwohkE11UF11aiaStQeM/769m5GgTf9uQRLr5MqaNfS/iOJLOhxTKTsLy5cspOVVC8f7TACT0C11f0aR4iNE5Q1blacKl6ThgS48olm6ho0EztVnqzKaZwq5em6jToiLaRCsuUnR1EWfSIzArXsytvtfGt23YGjQeX96V0uOTuGn8hoi+SDonUig7CXnHd6N5NZIGxDHyt32ZuGBIkE2U6iZObwcgWucMudJzanr/xU+hI5Eqd/CK0SX0/rqTBtUTsVsjgNUbjb2NOpVWLbRNnOogSVcP+MKLQomlAQ1TC4E0tAgbytns4JrxZhw1d3DzpNci+iDpvMg4yk7A2g9fZOOhV9EZVWb+fQymmKb0Qq8/xTBadRGtC9y6RulcqFpziTK71xAkVsXOBDxCR4rRF4zu1PQ4ROAlUFM6YKRK6vVeM5pQiNGF3z7Xa2Y0FH8L23i1ka6qPcBGCAWU5kAlAxrGVkJtVDQQ4BIKBYdd3DE7jik3PBt2XolECmUH56X3llHedxfl31hJHtDFL5LgC9x2Cx2xOgfmEAV5oblOZZU7NuQtOECZKx6PUEk0NIRt9aBXfDnUrRuYtcSmmdBQiNM52rZRHUEi2URLsWwtkk0YFY3iIo06qyDZMi3sfBIJyK13h6e87y6EEHy3u4qeQyxBn+sVLaxINpGgt4UVySYcmqHNfjgOzdjmeaRdM4bczre2cUUQXPCJZVSEc9bGRo1lj1mJijIz/z+ejjiWRCKFsoNzg3Mh5d/UYT1lZ9D4wG6AXqGiILB5w1cBilMdGBUPA8/08g6HUfVGHAd8l0CRtt9NNo2aKSBAvTXRqtOXG+5JCGsTo3rQI0L+gn973MNvf1XNl7tUNmzYSFxc6IstiaQJKZQdnNEjxtCv8AYAEtKaU/E8QvUHc3uELmRmDjQXoYjTOciIPhXSJv7MBVCkcZr680D4zoguTe/vwR3OJkZ1+iujh2tUFntGJJto/Uv++tpGDn/lYceOnVx//fUhx5BIWiKFshMwfdoMkpOTWT//SzSPT2Ra52V7hUq9JzCA3KILrCgUpbq4LKY40EbfGHCrHWocqyc6YCWpCQV7KyF0avpWWThwymUJeB2rOoIKXRxzdw94Had60EUIFDpW4GH9PzzccccdDBs2LKydRNISKZSdgB49evD+++9TdsTKNzllYbNlNBT/yq+rrjFkDKRR8ZAZW+SzMdhCp0K2GKfWEx0ygF0IBfuZlhROTR+yoIZH6Ch2+mpbxqmOkP54hI4CdxIAXdTwLWp931ihtOgZXC4XU6dODWknkYRCCmUnYeTIkfTp04fiA9URq/wIFPLqUiO2etChMSSmOGKtIHFGLEMVrWi2gYN1vSL6o6Fy2hMb0R9NqFRqsW3ULgKTAndOu5OuXbvyzjvvtGEtkTQjhbKToGkaDYRugdCSb22JaELhk5pLw9q4hQ4VQW9jdUQBEyhh204IoKAhCU0oHKlPDjtGtOpCCIVvnd3DzpSga0AIONqqsnqoWU+V9mHp0qW88847zJ8/H7s9dIiRRNISKZSdALvdzpNfT8HZ6MYUo8N8JjayJQKfSPpfC4LODCG4B04vY02b84cSy6MNSQGvj9QnB60+A31U+M6ZFOxTq5dH3V2DVqgKYDjzlqJqTL75SZ555hleeuklZs6c2ab/EokUyg5O1ekq/uv4bTjq3TRUOonv4Ts7NKkevxBpQuE7W3Cbg0+rBwcIY7hGYWoblcRbEy5E6JsW4qkgAnpxN1Ho7O5vNaEoggQ1uPf4MbfF3wLC15Qs8HNFgd/cvoKlS59g/fr1VFRUBI0hkbRECmUH56VS34qpaH81QkDa8ObYwyaxLGwMbuzVxGc1A7FrhojdFHsaa9tsJ9u0qnQLHcdDiHITbs23rQ8lkk2ccCaiKhpd1eBScU1867agKMEi2Yxg3ITn0DSNnJyciL5LJFIoOzj6Ul+Q+be7TtMl2Uy39MC2Bu7vEQSeWz2Ianfkdgg9jLWY2sjwASixWyJ+XtjYrc2e4QIw0LYwt25g1prUFD29eujJz8+PaCeRSKHs4CzKWoOpsB9HPyun36hEFKV5iVXrjuKgtRf1blPYkCGnV0+lPZbN5RmUOeND2qiKwKR4SDVWEx2mJqQmfMWAe5itmHWhBVWvaETrXRQ5Emj0hg4mF/hCgg7a08K2oFAVgVHxUOyNxh5mJWxWdCTpookyK+Tn7w1pI5E0IYWyE5DeMJqqQhtX3Jrmf6/SGcvXdT38rxvcJtxaoKjYPQZOO5pXkjmVAznRapuuU7SAcmw9jbVBRS28Qg3I8U421ROjDyyga1A1ovTNAlritPgroDehoQQU1Tjs6EWVJzD9sLU/pd4oGlrll0crehJ1vrPaCT+P4rPcbD7+QJZYk4RHCmUnYNWqVUyaNMl/PlnqiA+6dQaweZqLVtjcRmqcwSu2ndX9+KbBF86jU7SQBX6TDVZ/Vo9HqLi04FVdorGBOINv9WlUvSFXmWXOLtSeCVzXUEIeERx1plDqsfhfh/KnQjNTdyYTSAESdM0C/PjCBGpqNU4de5ScrS8HPSuRgBTKDs+uXbv48ssvmTx5Mn2LbqSoMYHvbOEvb+weA3VuM1ZX+H44+2vT0IcRySYSDfX0MVYGrVJbkmCw0d3UgEkXvhVEpSuO0+6YiOeohc7unHR3i3ihVKWZqNGMdFEDt/QpSXoSE1RKyzyMHfICn25eFnYMSedFCmUH5rPPPmPChAmMGDGCGTNmUOcpJ0oXudf1JTGVpEZFjo3MSvyK+FZ54KFIM55u06ab0YbFEHmsUw4LJQ5LRJsqTxzF7vD/AYCvHYRHBF8U9UjWU1bpE1m9sXdkhyWdEimUHRSbzcaMGTO44oor2Lp1K3/57Alq+n1JsqmeS+NCl0wbHFdOsrmeXlG1DLUUh7QBuCyuBJPiJlYNX2A3QeeLbxwecyKsjfFMi9qYCOJt9xpwayrVrvD9eaJUFwbFS5UnjmPO0Fk+KToHcaobm3Bz2hvod3J3HaXlXvZ89zxjxk0PO4+k8yKFsoNy9OhRCgsLmTdvHlFRUdj6HPF/1s1oY2h8SdAzCcbmuMQkcwMjEoJF7tc99vj/bFQ8dAlRZTxRV49Z8Z05JunruDL22yAbk+oJ6KzY3RgcON7oMeLRmn9Fvw0RfxmtutC3qF5U6w0OY+qpsxPd4oLHLjxUepv9Tu1l4FBBMqNGTw569nyyefNmVq1adUHnlPwwpFB2UIYNG8bw4cNZu3YtT31yV9DnXfQOLrcU+V9nhlhBWox2rkr8zv96Ru8d9DIHbsv1ijdgG56krwtqStZVZ+OauKP+12adJyit0ah6SDbV+V/bPEa8rVIabV5jwCVUjM4ZsnHZIUfz9jlV1xjUeRHAJZreUxgyYgkFR49TXBx+FX2u2bhxI9dffz0PPPDABZtT8sORQtlBURSF9PR0Dp3chTcpdDGMaJ2bK7sWMqJrEdFhtr+xeidjuh9jTlou3QzBqz7wVROy6Gwk663owwSCx6oOxnU5glkXnGfehF7RSDFZafCE7y3u0PQcaUghRucMW1LNoRnJs6eRrrdhCNMzRwClnkYqyCE9vR8AjY1tn7ueC9auXesv87Z9+/YLMqfkxyGbi3VQGhoa+ODj9xh336CIdpfHFAJw0J4W1uZ3CXvQKRr7nN3D2gwzVqAABW5LWJte+lpuid/Pv6zDw9oMMJfTz1xJdtXPwtpMSDiMXtGo8YTPFrox9kibGT5eBL1SevorCCUmhk+t/LGsW7eOiooKtmzZwqZNmwB4/fXXueaaa87bnJJzh1xRdkAqKirIGNMHISBjUq+wdkOjitApWsS+2zMtezEqXnQIRpgqQ9oMM1WgUwSqIhgUpppQd10DKgIDXv4tPnQmzCXmClRFYFC83JgYOq3wxsR8DGc6OnbVh17hAv6VZPhvptAj2XcGO3bsWAC/gJ0t33zzDTqdDkVRGD9+PCtXriQ/Px9x5oa9rq6OO++8kwULFvjnePLJJ5k2TXZ/vFhQhAgRL9GKuro64uPjsVqtdOkSvumTpP05cuQII8dmIjS4bcVIUof5gsxbxjwqimBoVFHQs/sb+zTbIJhlCRa0L51J/nWagk8kQ/rhaq4N2V3XEFSwTQDvWkcEzHeJuSLI5qPKoQHz3dT9YNBc1Z7mro0Kgl91CRbZlisCr1dH3/TA88hRo0ZhMpnIzc1FURQ2bNjAmjVrqKysxO1243K5cLvd/p+Wr1tv2Y1GIy6Xi9TUVGbOnMmsWbNYtWoVK1asIDY2lsGDB/PFF18E+Si5sJyNrsmtdwfCZrPx1M7ZNFQ6mfbfo/0iCb6qPQbFi6poXBZ1MuTzw6MLyWtMR1UEM+L3hbQZYapgn6s7XqGEFUloKkqhoCih65f7yp958ZwpAtzPHDxWkzB+UDkUvaJxQ+KhkHMl6Buo9sSgVwS3xIW20fCJpaMxml69DrJu3Tq2bdtGdXU1JpOJqqoqvv32W1RV5dprryU3N5exY8dy+eWXYzQaMRgMQT8t39fr9RQXF/OLX/yCzMxMPv/8c9avX8/zzz/P0qVLmTBhAl6vF6vVSl5eHocPH+bSS8MXR5b8tJBC2UGoqCjnpfK7Serv+5+x6rsG+l0VeKboFjpGRIePawQYFn2Cvm3caVxurGR7xUBICi+UAw21HKjsR1pCaFEGuKXLfjbVDyPFWB1xvi46J2MTvolo09NQz9VRhRFtqqq6kbd/AUuWDOTUqVMMHTqU5ORk7HY7J082+2mxWHj99de54447AoqInA0TJ05k4sSJPPfcc7zzzju8+uqr/s/69OnDuHHjmDt3LnfddRd9+/b9QXNILhzyjLKD8NeqGUTpXOxcewy9WcfgCSlBNmbVhcFf0jY0tvKeDBmST1lZelibsrJ0fn15LmWlEWxK07ll+OdUl4e/JDIq8FC3fOLU8CmMMaqb5ambiIlQwq2L6uU3cVXow/w611k1/t9jKjdMqmDWrFmMGTOGY8eOkZeXx+bNm9m2bRuVlZW8+OKLVFRUsHHjRqZNm/aDRbIlsbGxzJo1ix07dlBYWEhVVRVbtmzhlltu4YUXXmDQoEEsXboUt7vtEnWS9kOeUXYA9ufvZbNuMbWldp69aTtXz7iEcb8P3NbF6hxcZm5eNbVu6iqE4NM3YrGfyqSiogKr1Yq98SAZP2tkyM8MTJ5sJipKobBwKNf8PNv/3GfbryetT+C5YVHhUMa2sPlo242kpAVuic0q9NU3C9vmxu6cblVazaI6+XlUcy/x7fae1GqmAJtEnYfJMc0r0jJvIy7hu8Jx2AV/fc3GSy868HoNzJkzh+nTp3P55ZeH/Hu80NhsNpYtW8YzzzxDRkYGs2fPZuzYsWRmZp4TkZZE5mx0TQrlRU7uzi18Ef8cHpfGmll7qCmxM//9n6NFNYfOdNXZGGw+FfSsBxWBgqNR4+m769n3xUkuu+wyevXqRXx8PI2NjeTl5XDyZD1duypMnDiI+fPXMnLkSFTVJ3KapvHCyomYo/Yz+mojpyvHk3XdW0Fzvbf1dnoP8MUMxqiQpg9e/eXYu1Hq8VX2SdLZudIcnGq525FChddXsKO33sXE6Nogm0qvnX1fO/i/i+rIO+Bl7tz7WLRoET169Aiy/Smwd+9elixZwtatW/F4PAwbNozFixdz6623SsE8j0ih7CR88Ml6vkr5CwDvLTvMrrdPMudvV5KeacGl6fEIHd31dfQ3lfufsds0st+ppSDfjtB85cuO5+mwVrt48803+eUvfxk0z9q/3s6//pXPlq3fYbPZSEpKYvTo0VRXV3Ps2DFKS0sBXy+aESNGcsMNNxAXF0dRUREnT56kqKiI6upq6uoqGXGlYPrUGK7INJLWW4+qBgrBF46uuIVGZphQJIADzu6YFD1jouqCPjtZ4uax5dW8vr6e9D46/v63TxkzZswP+vu90DgcDnJyclixYgXZ2dmMHTuWlStXkpmZ2d6udUikUHYSnvn6BgAOfFjKWw8f5JbHMrjq16lnPlVIOH4dqT97jaP5dg7saORUoYsdW+tprPcycGgUBoOCpzGZ9N6X8sc//rHNf5A2m40vv/ySjz76iAMHDpCUlETv3r256aab6Nu3L5s3b2bTpk1s2bIFr9dLamoqaWlppKamkpiYiKqqvPP2Xyg46hNWS7zKtNtiuG92HH3SDH6/t+VNI/Oyv0TwROHrQ3czdtgb/ndqrV6eXlXDn1610iVW5fGFCaQN+R2/nHxxlk3Lzs7moYce4siRI9x99908+eSTdO8ePuBfcvZIoewkvLVpDSfS3uX5m7eTmB7DnX9qOttSUHKvZPvH+6msLuKLz/OIilHpkWbgslHR/NvsbiT1NOI6dh83T7r/nPslhIi4ZTy+vzeHDrvZvtPBa280UGvVmHJjNA/dF4/H8CjjJ85j65Y/MTTj6RBPKxz8+j+ZmPUAOVtfZNSAFby81sofV1bjcAjmz7Ww8L4E9hf+B+Mm/v6cf7cLidvt5uWXX+bxxx9HCMHjjz/O/fffj9EYuk2G5Ow4K10T3wOr1SoAYbVav4+55AKy9u1XBCCM0Trxfz4eIx77fIKYfMsNIioqSvTv319kZGSIe+6bJj4oGCw+Pn7pmZ8Mkf3pm+3m82c5b4jGU+mi8VS6qDqWKlY+nSAG9NMLRVHE/PnzRV1dnRBCiBPHUkR5cY8WPz1FTs46/zgbN24UPVK6C1VF3D2tiyg+0Ed4SweIbTlvtNdXOy9UVlaKe++9V6iqKgYNGiQ++uij9napQ3A2uiaF8iKnpKRE9OzZQ+gMijCYVRETEy0sFotYsmSJaGho8Ntt37kpQCjbm11ffOwXy8ZTfcSOzz8Qzz77rMCXkCMAodPpRHq6TmQOM4jhmQYxePBAkZWVJdasWSPuueceAYgbb7xR5OekCW9pf+EtHSB2f/Fxe3+180ZeXp4YP368/3sfOXKkvV26qJFC2QnJ3faJuOGhgeLRRx8VxcXFIW0OHNwlPj6eIQ4c3HWBvQtN/sFdovFUH5Hfwp8vvvjCL5SDBw8W1103TkydYhajR48KENGoqCixevVqoWmayD+4S3hLBwSM01HRNE28++67om/fvkKv14uHHnpI1NTUtLdbFyVSKDspx7492qbN8W8LLoAn35/vQvjzwQcfiFmzZgUIY9OPoiiioKAgYLUcbpyOjN1uF8uWLRMxMTEiMTFRrF69Wni93vZ266JCCqXkoqe6ulo89NBDIiUlRQCiV69e4umnnxa1tbXt7dpPipKSEjF9+nQBiNtuu000Nja2t0sXDWeja/LWWyLpAGzYsIFp06YxZswYPvroI39CQEuEEOzbt4+cnBwmT57MoEGDKCkp4brrriMjI4N//vOf7eB5+3E2uiZzvSWSDsDUqVPZuHEj2dnZrFixIujzdevWYbFYGDFiBAsXLmTw4MEoikLv3r355ptvOHQodNUliQ8plBJJByErK4sHH3yQxx57jKqq5vYf27dv584776SuzpfJZLFYAp6bOnUqR44cQRIeKZQSSQdiyZIlALzwwgv+9z744AMArrrqKoqKiqipqUEIgdfrpb6+nn/+858ht+qSZuTfjkTSgUhMTOTee+9l1apV1NT42nI4nU4AHnzwQVJTU/22qqoSGxsbchxJIFIoJZIOxoIFC3C73Tz88MMIIcjIyACgoKCgnT27eJEVziWSDkZKSgp//vOfmT17Nk6nk4oKXyX6bdu2tbNnFy9SKCWSDsisWbOw2+0sXrwYq9UKELKEnuT7IbfeEkkHZO/evezYscPfMnf8+PHMnj27nb26eJFCKZF0MN59911Gjx7Nnj17WLBgAbt37+Z///d/iYmJafthSUjk1lsi6WA88sgjTJo0iQ0bNmAwGNp+QNImckUpkXQwXC4XqampUiTPIVIoJZIOxt13382aNWvYuXNne7vSYZBFMSSSDobH42H06NEIIdizZ4/s5BgGWRRDIunE6PV6/vCHP7B3715Z7OIcIYVSIumATJw4kbi4ON588832dqVDIG+9JZIOiMlkYu7cuTz99NNs2rSJlJQUunTpQnx8PP3792fu3LnyGO0skGeUEkkHRQjBm2++SW5uLlVVVdTV1WG1Wjl06BCpqamsX7+eYcOGtbeb7Ybs6y2RSMJy7Ngx/v3f/50jR46wevVqpk+f3t4utQvyMkcikYSlf//+7Ny5k1tvvZXf/e53lJSUtLdLP3mkUEoknQSv18vq1at5+eWX0el0PPbYY7jdbvLy8trbtZ888jJHIukkbN++nblz5wKwZs0aZs6cCUDfvn3ZvXs3ffv2pXv37u3p4k8WeUYpkXQShBA89dRTLF68OOTnt912G2+//fYF9qr9kGeUEokkCEVRWLRoEevWrfO/N2vWLP+fly9f3h5uXRRIoZRIOhm33347EydOBOC1116jd+/enDhxgvT09Hb27PzjdDopKiriwIED5Obmfu/npFBKJJ0MRVF47rnnAJg7dy6HDx8mLS2tnb06v9TU1PD444+TlJREeno6w4cP5+abb/7ez8szSomkk/LrX/+aHTt2UFBQQFRUVHu7c04pKSlh48aNFBQUUFxczNatW3G73dx7771MmjSJhIQEjEYjmZmZMuBcIpGE5+jRo2RkZPDEE0/w6KOPtrc754w33njDf6M/YMAAevTowdVXX839999PcnKy3+5sdE2GB0kknZQBAwbw4IMPsnjxYrp168acOXPa26Ufzd69e7nrrru4/fbbWbVq1Tlb2EmhlEg6McuXL8flcnHPPfdgMpmYMWNGe7v0g2lsbOSOO+5g6NChrFmzBqPReM7GlkIpkXRiVFVl5cqVnD59mkceeYTbb7/9om0hsXbtWgoKCvjqq6/OqUiCvPWWSDo9iqLw8MMPU1ZWxqefftre7vwghBD8/e9/59prr+XSSy895+NLoZRIJAwdOpT4+Hjef/99vF5ve7tz1uTm5rJr1y4WLlx4XsaXQimRSFAUhXnz5vHiiy8yYMAAnn/+eerr69vbrTbRNI3c3FxmzZrFyJEj+cUvfnFe5pFCKZFIAFi6dCm7du3i6quv5pFHHuH666/H7Xa3t1sRuf/++xk3bhwWi4W33noLVT0/kiaFUiKR+LnyyitZt24dubm57NmzhyVLlrS3S2ERQrBjxw4A9uzZQ79+/c7bXFIoJRJJEKNHj+aJJ57gueeeIz8/v73dCclbb73FwYMHef3119HpdOd1LpmZI5FIQuJyuRg6dCgWi4VPPvmE6Ojo9nbJT3l5ORkZGWRlZfHWW2/9oDFkmTWJRPKjMRqN/O1vfyM/P5+pU6fy6aefYrPZAJ/IHDx4kLy8vHa5JZ8/fz46nY5Vq1ZdkPmkUEokkrCMGjWKjRs3sm/fPiZMmIDFYiElJYX4+HiGDRtGZmYmw4cPZ+/evefNh3vvvReTycTgwYNZvHgxNTU1/M///A8LFy68YBXZZWaORCKJSFZWFuXl5Xz99dds376d0tJSBg4cyCWXXEJ9fT2PPvoo11xzDa+88gp33XXXOZ07OzubV155haysLNLS0li5ciUbN27E5XIxcODAczpXJKRQSiSSNlFVlSFDhjBkyJCgz8aMGcPvf/97Zs6cye7du1mxYsU5SyHMyckhJSWF7OxsFEUhKyuL3/zmNwD06NHjnMzxfZBbb4lE8qMwm82sWbOG1atX8+qrrzJlyhQ8Hs85GdtisVBWVkZpaSkAU6ZMCfjsQiGFUiKR/GgURWHOnDl8+OGHbNmyhTlz5uBwOH70uL179/aPDz5RXrlyJcOHD6dPnz4/evzvixRKiURyzsjKyuK1115j3bp1dO3alWuvvZZnn32WhoYGvF4vVVVVZzXexIkT0el0/OMf//C/98ADD7Bv3z5MJtO5dj8sUiglEsk5Zfr06Rw6dIhly5aRkJDAokWLSExMJCkpiaSkJO655x4KCwsjjuHxeNi+fTtPPvkkXq+XpKSkC+N8GGTAuUQiOa8UFRXx7rvvUltbS3R0NMuXL8dqtfLb3/6W1atX+wPZKysryc7O5sMPPyQ7O5uamhqSk5OZN28eixYt8m+/zxVno2tSKCUSyQXFZrPx17/+lQULFjBv3jxGjRrFn//8Zz777DOEEFxxxRXceOON3HTTTYwYMeK8pSeet545dXV1P8oxiUQiAd/2/MSJEzz77LMAXHHFFbz44otkZWUFNABrygQ6H5yNnn2vFaXD4aBv376UlZX9KMckEonkp0RKSgrfffcdZrM5ot33EkrwiaXL5TonzkkkEslPAaPR2KZIwlkIpUQikXRWZHiQRCKRtIEUSolEImkDKZQSiUTSBlIoJRKJpA2kUEokEkkbSKGUSCSSNpBCKZFIJG3w/wFPazm6yVaSCwAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 400x500 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"import cartopy.crs as ccrs\n",
"import cartopy.feature as cf\n",
"import easygems.healpix as egh\n",
"import intake\n",
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"import shapefile\n",
"import shapely\n",
"from cartopy.io import shapereader\n",
"\n",
"\n",
"def get_shapes(resolution=\"10m\"):\n",
" \"\"\"Return a dictionary that maps ISO A3 country codes to shapely geometries.\"\"\"\n",
" shpfilename = shapereader.natural_earth(resolution, \"cultural\", \"admin_0_countries\")\n",
"\n",
" records = shapefile.Reader(shpfilename).records()\n",
" shapes = shapefile.Reader(shpfilename).shapes()\n",
"\n",
" return {r[\"ISO_A3_EH\"]: shapely.geometry.shape(s) for r, s in zip(records, shapes)}\n",
"\n",
"\n",
"# Load ICON online catalog\n",
"cat = intake.open_catalog(\"https://data.nextgems-h2020.eu/online.yaml\")\n",
"ds = cat.ICON.ngc4008(zoom=9, chunks={}).to_dask().pipe(egh.attach_coords)\n",
"\n",
"# Select regional polygon (e.g. Portugal) from HEALPix grid\n",
"shapes = get_shapes()\n",
"polygon = shapes[\"PRT\"]\n",
"\n",
"icell = np.where(\n",
" shapely.contains_xy(polygon, ((ds.lon.values - 180) % 360) - 180, ds.lat.values)\n",
")[0]\n",
"\n",
"# Plot regional cutout\n",
"fig, ax = plt.subplots(figsize=(4, 5), subplot_kw={\"projection\": ccrs.PlateCarree()})\n",
"ax.set_extent([-11, -5, 36, 43])\n",
"ax.coastlines()\n",
"ax.add_feature(cf.BORDERS)\n",
"egh.healpix_show(ds[\"tas\"].sel(cell=icell, time=\"2025\").mean(\"time\"))"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.12.8"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Please register or sign in to comment