manuscript_figure_S1_plot_nested_domains.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 79,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-07T01:26:50.676315Z",
     "start_time": "2018-06-07T01:26:50.594129Z"
    }
   },
   "outputs": [],
   "source": [
    "import xarray\n",
    "import matplotlib.pyplot as mp\n",
    "import cartopy\n",
    "import matplotlib.patches\n",
    "import matplotlib.patheffects"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-07T00:59:56.777111Z",
     "start_time": "2018-06-07T00:59:56.751928Z"
    }
   },
   "outputs": [],
   "source": [
    "root = '/beegfs/DATA/pritchard/blangenb/WRF_ARCHIVE/CHEYENNE/nested_FMAASO_10days/400/001/'"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "open example files for d01, d02, d03"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-07T01:07:35.443569Z",
     "start_time": "2018-06-07T01:07:35.407905Z"
    }
   },
   "outputs": [],
   "source": [
    "d01_filename = 'wrfout_d01_1979-01-01_00:00:00'\n",
    "d02_filename = 'wrfout_d02_1979-01-01_00:00:00'\n",
    "d03_filename = 'wrfout_d03_1979-01-01_00:00:00'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-07T01:07:50.591809Z",
     "start_time": "2018-06-07T01:07:47.160977Z"
    }
   },
   "outputs": [],
   "source": [
    "d01_ds = xarray.open_dataset(root+d01_filename)\n",
    "d02_ds = xarray.open_dataset(root+d02_filename)\n",
    "d03_ds = xarray.open_dataset(root+d03_filename)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-07T01:08:09.062429Z",
     "start_time": "2018-06-07T01:08:01.936979Z"
    }
   },
   "outputs": [],
   "source": [
    "XLAT_d01 = d01_ds['XLAT'].values\n",
    "XLONG_d01 = d01_ds['XLONG'].values\n",
    "\n",
    "XLAT_d02 = d02_ds['XLAT'].values\n",
    "XLONG_d02 = d02_ds['XLONG'].values\n",
    "\n",
    "XLAT_d03 = d03_ds['XLAT'].values\n",
    "XLONG_d03 = d03_ds['XLONG'].values"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-07T01:08:20.764895Z",
     "start_time": "2018-06-07T01:08:20.544058Z"
    }
   },
   "outputs": [],
   "source": [
    "latmin_d01 = XLAT_d01.min()\n",
    "latmax_d01 = XLAT_d01.max()\n",
    "\n",
    "latmin_d02 = XLAT_d02.min()\n",
    "latmax_d02 = XLAT_d02.max()\n",
    "\n",
    "latmin_d03 = XLAT_d03.min()\n",
    "latmax_d03 = XLAT_d03.max()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-07T01:10:02.875590Z",
     "start_time": "2018-06-07T01:10:02.649709Z"
    }
   },
   "outputs": [],
   "source": [
    "lonmin_d01 = XLONG_d01.min()+360.\n",
    "lonmax_d01 = XLONG_d01.max()+360.\n",
    "\n",
    "lonmin_d02 = XLONG_d02.min()+360.\n",
    "lonmax_d02 = XLONG_d02.max()+360.\n",
    "\n",
    "lonmin_d03 = XLONG_d03.min()+360.\n",
    "lonmax_d03 = XLONG_d03.max()+360."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-07T01:08:43.813750Z",
     "start_time": "2018-06-07T01:08:43.781716Z"
    }
   },
   "outputs": [],
   "source": [
    "d01_height = latmax_d01-latmin_d01\n",
    "d02_height = latmax_d02-latmin_d02\n",
    "d03_height = latmax_d03-latmin_d03"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-07T01:10:42.337716Z",
     "start_time": "2018-06-07T01:10:42.302622Z"
    }
   },
   "outputs": [],
   "source": [
    "d01_width = lonmax_d01-lonmin_d01\n",
    "d02_width = lonmax_d02-lonmin_d02\n",
    "d03_width = lonmax_d03-lonmin_d03"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 90,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-12T16:19:28.015089Z",
     "start_time": "2018-06-12T16:19:27.494726Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "-2.243599"
      ]
     },
     "execution_count": 90,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "latmin_d03"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 91,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-12T16:19:36.742477Z",
     "start_time": "2018-06-12T16:19:36.689260Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "12.03965"
      ]
     },
     "execution_count": 91,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "latmax_d03"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-07T01:10:14.367482Z",
     "start_time": "2018-06-07T01:10:14.318924Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "277.7049255371094"
      ]
     },
     "execution_count": 41,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "lonmin_d03"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-07T01:10:25.843791Z",
     "start_time": "2018-06-07T01:10:25.799454Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "313.6654510498047"
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "lonmax_d03"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 89,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-06-07T01:36:08.456386Z",
     "start_time": "2018-06-07T01:35:50.209545Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcwAAAFGCAYAAADn8K84AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xd4FHX+wPH3zvZkk02vpJKQQhKSAKFLb4LYu+dZztM79eye6Nk9vfM8bGf9eWc9FU8RQZHeayAEQgiEENJ7T3azdWZ/fwSBEEqCCaHM63l4HtiZnfnssrufmW/5fBUulwuZTCaTyWSnJvR3ADKZTCaTnQ/khCmTyWQyWTfICVMmk8lksm6QE6ZMJpPJZN0gJ0yZTCaTybpBTpgymUwmk3WDnDBlMplMJusGOWHKZDKZTNYNcsKUyWQymawbVD3Z2eBvcJnrzX0Vi0wmk8lk/aHE5XJFnm4nRU9K4ykUCteNX9/4a4KSyU6oKM/Bzve3ozaGEjDiWtR6D1Q6d3Qe3kSnTEFQKSnKWUvOv/+IrWIfxvThhGYMYMgEHQDmNpHcNWZKl2zHPXoELknE3lCIX3okLkmkraQJwwBvqtdmohDUTP1nDp6+IWccr91q5sCWH6jcMp/WvNUYU1KZ8XDYCfc1tUgsvns+hoQUUm+KIyRaQKns2riz8r1qRHsgUx756ozj+rWyl35IyQ8vMeEvo/EN7NsGqIYakVVPryHp9/8mbuRlfXqu3rb+owcxFWdjb6rA0VSFQqFA7R2AxtsHjY8Brbcad1897n5qLM1OGgvbaCusxFZTij4sCc9BY/FPGEdE0njcPH36++WcVdXFueR8/Qxt+esJmzmC9Es90Or65rPmtLvIWdtO0U85KN39iLvmWWKGXoqgFMhZ9TkFXzyMZ0I8TZkbcblcitMdT06YsnPCynerUXkMZ8ytryCcIJkAbP/hDVoKM1G7+6J086J65TvEXDeW2pwGWvP24BE3jvgrniAiaRySKLFn9Wc0F2VhqcrHWlOIo6UGz8SJhI65icRx15/0PD3149wx2OsLSblrLDFDVCdOhh9UUrdmHQCC3hP/sRlEZBiJTFChVHXsbzZJLP7DQibP24dPcFSvxHYmtnz9ApXL3iLulnEMGa/v9eObmiWyFjVRvS6TAZc+zMjrn+71c5xNkijR3tZAS20JLbXFmOrKsDZWYGssxd5Ugdrgi0d0Bv5xIwhPGING59bfIZ911UV7KNz8LY07f8BWX0LIlJEMm2NE7943idLUIrFraQtVa3aiDRxEzGWPETt8NoJSwG5tZ93bt2Eu2kb6PRlEJar56oav5IQpO3eJTonslRYa9jdjbzbRXlpI4l0fkjjuum4fY9u3f6Nu+wICRt5I4oSbMHgH0lxbyo4v5mKrLwaFAtFiwm/oHEbf/EKfvRZTcx17V35M1dp/g2TFJ3UQISmeRCapO105m9sktvy3jrq1qxEMfqjcjYiWVjxiYzFGGQEXNdvKmf3a7l5L5meqeM86st+5Fc/4KC65Ixid22l/S06rpUkke2ETNRu24502m5RrnsI/LL4XopWdayRRojx/KyVbF9Cw8wckSyveQwYTluHNwBQNKs2v/zydSHODxI7/1dKQuRNjynTiZz9MWPyITnGtmncjkrWYSQ9Go9N3xCEnTNk5q7ZCZOsH+Ug2gZAJd+DuF4aHfzjBkSkIKuUZH7e4IJvMuekETZ+Kb7QbRSsP0l6wl+AZ9zPud28d2a+xppiizB/xixpCRNK4jpjK9pG/4t80561GbQzALTgBj9AEfMIT8Q6KxmAMxOWSOLR7FTV71xM0eDzRqVM6xSuJEiV5Gyjf8RPNeauwVubjHh2PT2IAQfFuBIapMHgJlOY7yf5sH5LNRejUexEEJU0Ht2Eu3knMVU+TMOaqM39ze1F7axObPrgHc/E2Rv4pg9CBPRry0MXaT+to3FPNmMcW4Bca20tRys41NSV7yXz7VpzmOnzT4ojI8CIiXnmkJaUvLflbEYIuiuG/eQWvwIgu20v2bmT7q5dx+b9mHEmWICdM2Tlq+xIzh77bQOjMh8i4+gmUanWvHbswZx1ZL0xA7RuK2N6G15CZDL5qLsHRQ47sU1Wcx6anR+GZOJiWnGwueWYlgQljWP7KlYj2MmJnRGBtE2mpsNBeZcZSXYujuR7RYkKhVKMNDMNzUCgt+0sR2y0Ejb+DIZc9gLvRt0s85pYGirKXU5uzHFNxFvbaYlAq0QaEovf3wVRWir2qhEv+kUNQVHKvvQ+9bccPb1K+bB6z/j4Knf7Mf/QaqkVWPrmcqfNyMfqF9mKEsnNF9s8fcGj+XMJnj2H4Ze4n7J7oKyX7nWybt47Z7xSi1RtOuI/TbmPxA/Ek3ZpO/HDtkce7mzB/3SWjTNYD2xa1UbZsN2OeXdMpifWWgSnjCfyohva2BnTuXnj6BHfZR3JYkKxtSA4HmoAovKNSATBEpGJvbiVu6C9fIjfAFwgHwGGXsFnAYPzlByCQ0nwn+5as5ef738Z/3K2kXf1Ep3O6G31JmnAjTOi4yJRECVNzNfWl+2iqyEfjU0Czei0WU2Ovvxe9adjlD9C4fz0bP61jyj2BZ3wcL38FaqMPlQcyMfpd2YsRys4VJT/+g5S7J3ZKRmfLnv8VETb7sZMmSwCVRkv0Vc+S//0/iR/e84vUHifMqhInwRFynpX1TNZSM6VLshj77CoCwhL67DwGrwAMXgEn3R4aO5Sxf9tJec5qojMuQ6NzB8B/0Ajy5/8POHFCUGsE1JrOj4XHqQiPG0B9VTA7v93F2hemkvTbN1Hr3FHrDWj0BgxegWi0HQNnBKWAp29Ix+jctMm98nrPltF3v8fyx9JY9SGMuN7/mAuH7svdaANBR1zGnD6IUHYu0PiG4XRIZ/28NaUipgN7GP/oitPuq/P0RSGc2Z1vjzPf5tczufqN0Z0es5glyg44qD/UTnutjcAkI/EZWtRauS6CDHI3WylcsJUxT6/o02TZXSEDUwkZmNrpsdC4keyqOITDnoRa07PPrV+wksl/DGb9pzryvnoQ0WZDsluR2k3oBiQx66/rezP8fmHwCmDyK5ns+O9TLHlkMT5DUxBUAnC4FUsB7gE6Uqcb0J7ge+90SBz8Lpv4W974Vf3UsnObzj+atuoqoPdHV5+Kb5CAR+IQ1s27nslPLDxykXoidfs34xF98ovqU+lxwpz9t5FdHlv9VhHW6gY84y9B5xvBoWVLyftkP95pQ/BP9CJ8sBbfQPlLcjFyOiXyv8lm8F0fEhSV0t/hnJSbhxe6AYPJ22xlyISeD/tXKgUm3uFLRzNuh/wsG/kLq3p0nNaGStRad/QGY49j6GtGv1AmP/AJFQU7KctexpHxD66OO4q6nLUsXr6LmKuGkTxe16n/KntFOypjCINGyHeXFzJ9QAzm2uKzfl6VRsH0R6JZPq+IpXNHknzbm0SlTDjJvjrKV2xlSW0SURODSRih6/55ehrYscPkRVGiqdaFUq/CmDCBifd91LHhpmepryigYOM3VG7dwIH/rkep1WOIicYzwgOX5MLR7kS0SKjclPhEuxMySINvgJxULzT7tthR6rwZlHHuT0yPnv0ohQueJmlceq8MVujII90fVLdr+X84+PlD6ILjmPbC2nN2vl5obDqhsekn2PIXCrb/zL4vH6d4uYhXQiiSXUK0STTl5JL2p/n9Pl1G1rc8g6NpzmvjZF0bfUmlUTD90Uh2r7aQ/db1FA2ZxaT7/9Nlv5E3Potp5h/I3/g/cj/+Ky2VJ/osn+QcPQ2qplRk/9pmmvOrsVYWoVBp0QZGYxzWeZi4X2gsftc/BYDkFKkpzaVy7wZainciqPWo3LzQGz2xt9VTsjqLfZ/sBkGJe3gkhggvfKMNxKRrftWoPFn/a65oxxg//rz4ofQKjsFaWUR9VSqBA359vANi1ewuzcVpt6LSnPoqdstXz1Oz/t+MeHwa+xaXs2bejUx9fMF513wZO3wmA9OmsWf1Z7RWFaAxGFBp3QjO+BPRqZP6OzxZH/MKjqWgrhaI6ZfzK1UCQybpqdoeccprVYNXAENn30tE2lQ2vjyr28fv8bQSpcEX3xHXEjJ0FkED039VebFjSaJEY3UhVQcyaTqURVvhNqzVBxgwLYP0WR5y4jxPrf2kDkE7kjG/ebG/QzmlgzuXs/vtm4i5bgzpU9177bjfP5ZN4q1vEZM29aT75Kz8lIPzn2DSc+PxDRSwWSWWvpiLW+hwJvzxo/MuacouXk67jYV3BTPlr1PxDTr957aswEn2Z/k429sZ/1h6t55zOus+a6A5v5kZL6477YXqL765RtE38zCv+tLS7SB+rYoDO8j95jlMh7bJifM8ter/amjJr8Vz0Bh0vuEYQ2KIH3X1OZUE9q6fz/7//JGUuyf0+nD4NR/X0byvFq/BUwhIGEf86Cs6bT+Q+SN73r2VUX+eTFjs0Tmp5jaJFa/swS04nYn3f3JOvV8y2amsnHcTHiHNjL7G66T7mJoltn5VS+POPYRf8RSixUTlqncY/5dLCAjt+We9qsRJ2Z52Gg+00FZQxORXtvVorm+fJczrvu3+/r2lPH87e//3HKZDmQSMTiNosAeRgzW9Uq5L1rca60TKcu2YGqxYG+w07cln4DUvkzLl1v4O7Yit81+kYumbRF89mrSpbr062dra7iI/00JzWTv1mfkYB09j6I3Ps3/dV1Rv+BSnqY6k20eTkNE1UZtNEitf2YsuIIlJD3yOoJKnc8nOTN7Gb9n/6QMo9Z4odQZEWzuiuRndgMFc+tzyXj1X/tZF7P3o94TNGErCeANG76MJ0GaV2LXcTMniLfgMvYLht7x8ZBrYps//gql4JZfOje7R+UwtEj/96UeMabPwihnFwBGX97gW8wWVMH9Rnp/JoY3zac3fgKV8LyqDFyiVKISOP4JGg8rNDX2QkfF3+qM6C6WYZD2Tu9nKge/yuWzennOqX7M0bzO7PvoDSi2MuS+hT1bqMLVIrH3rIKYDuRhThhI5PojYdPUpP6fWdhfL/5aHxieeyQ9+Kd9pynqsoxXjtwz5wwT0BiU2i4vaAybKftpI2qPfE5M2pVfPJ4kSBdt/pHTdp7TkrsAQm4BCqcRSUYmjqRr3gSMYcus/CR00rNPz7FYzP94XS/q9oxmYrDnJ0bvat9VKwZIaZr286YxjviAT5rFsFhOt9eVIkogkOpFEJ3ZLG9a2Rg4ufIWQ0UFkXObR32HKjiOKEov+vJPIy54mZfJv+jucTiSnk02fz6UpZxGXvpjWJ0sOiaKE00GPjm21SCx7aS/uEWMYf/e759SFhuzc1lhTzJo/DyX9vonEpHYkofyddrJf/wn3+Am4BUaDywUKBRFjriMqeXyvnr+9tYn8zd8iKNX4RQ0hMDwRleZoa4rkdHJw51IqdvzE0BuepWDzAqo2vMfsFwZ3u6Vn3Wf1uEhh3B2vdTsuySlSuGsFZZvnEzBkBnlv33Bhl8bT6g0nXenA4BNM5quXkXjJlDOqSCLrO0qlQNLNKez64M8ERA+hrb6ClsoCLM3VeEemEJ48/oQl7c4GQaVizK1/Z+U/Cln7fgXTHzzx+pa/hlIpoOzhTaJOLzD58cEsf34lW796htG3vNTrcckuTPlrvsBrSPKRZAlQu9+MR8IQ9P52tMZDoICq9XsxBA/q9YTp5ulN2oy7Oj3W3tpIS30ZzdWH2P/5wyh1Wux1FbRNvYvU6XdRvuwtdq+1kD65e4PvWg/WEzVr9Ol3PKyhqpD1z09B0KrxS4sk/9P7u/3c8/YO83TWf/Qg9du+If7mkSSPPTuDlGTdt+qDjg5/bUAwOn8v1B4qzKXNmIsPoPIMwDBwBP7JU0kYd+2R8nVni629jUV3BTP7zdkYvM6dC67GOpFVT68j/rdv9WgZNNnFa9FDySTeMOiYGsldWcwSi+77mSnzcvHy7/2LxGMd3LmMXa9fh8rTG5W7BxFTBqIQFBz4Zgez3shD6+ZBad5mtr96GeOfm0zggFPf07U0iSx9aAkz/lV4ypKYv7DbLCx7cgx+6UGMvb5jUFJVsZO1T/zvwr7DPJ1LfvcGB1Kmsvfj+yjbGMToO8Px8pf7f84Vk+8OAI6vpxqA0xlDVZFIRV4zFevfpuDzB/EZdiVDrv0LvsEDz0psaq07Kk9/GmukcyphOmwuJIcdv8hzd2UT2blFcthoKLLCSRKmxSyx/OW9+GRc3efJ0mm3kvPubSTfM4XEEUfj2fJ9K4JKQ/aiN8m49knCE0dTNetRtrz9GbNeHHLKUpXbPqvGb8wt3UqWACV71uGSzIy6xvPIY8GR3U+D586vQR8YlDGLWa/novJMYMvH5f0djqwbVCqBsFg1Iy83ctmzcUx5eTpqfQVr52aQufB1Wht7VmruTGQteh1BJRAQfm59PcpyrSA62ffze9RXHOjvcGTngXFP/UzFiixK9ju7bOtIlrkI2kAybnmlz2OxW9uRLG2dkiWAwU+No6WG6o2fI4oOyvZtxSWokTCQtaSt0747fjbzw1O5LHmliGWvl9GSl0/GzS93P4b2FjTe3mc8Ev6CvcP8hUbnzvBbXmb5A4OwWsLkeZznGd8gJRNu96Ns3GR2f/0ly797FqW7D4aBwwkdeR1xI6/otTU1nXYbOxb8g4plbzD+mUnn3Gdl2AwDYYkz2bssj7VzR+ARP45ht83DN6R/qqrIzn2+wQNxC0+mYG0tBasVjPltIHr3js91faWEZHfgaMnn5/sGEnXdSwy7/ME+i0Wl1SM57Z0ey1lvYf/nGwm74i/U56zgh9v9UHv74zYgEFvJDiqdcXBFR+m6sgIHh77bSMKd7yM5HVjb6ombnYGbp3e3Y9C4GXG0mc74NVywfZjHW/rCDOyNRQSOiiFlssc51dQm6z7RKVFd5qI810zV5iIcTY14p1+O2uCLUqND0OhRadzQG/1x9w7Cd0A8Bi//0x730K7V7PnPfag9dWTcEUtg+LndfG9td5H5fRNVqzMZMOtRokZcgX/oIHmu5kXObm2ntb4MW3sb3kHRuHn6sOIf19KyZwUKpYbgSemMu9Gny/OqSpxseHEVGX/+kbCEUX0SmyRKLLjdhxGPTSUyQY25TWTRPd8z/MmlRA2ZxI9zR+MRoWPcrR3Nqz88sQtnfRmGmEQkhwNrdTnR1/2VtJl3n3EM7a1N/HR3MJe/f8WRCwfo/gLSF03ClJxOCrJ+pnTdp7TuX8/g28YyeLQ8GOh8V3bQScmOVpw2EckhIdldiHYRR5sFa00VuuAEZj7bdWK2027l0K5VVGb9RPOeZUiOdmKuHUHKcatsnOuqip1kzy/DXFKK2NaIJiAKt5AEoib/nphh0/s7PNlZIokSuWu+4OBXj4NSiVKjxWluI2zOXILiR9NQsoeDn97PmGcuIzhSiUrd9TO+a3U7BxfsYvqrO3Hz7JpUe0POys8o/O4vXPH6aFQqgSWvHMIr4UpGXPNnqov2sPPjB7DVHkAfEozL5UPyjS/TVlfSsYKP0Y/Q2KEnff115fupzt+GR0AEgVFDcDd2rBxkbmlg54K/01qwifaS3Wj9Q5g4Nx0vXzlhdsv+LT+Q99HvufaD82sRX1nPbPi6EXtbNBP+8G6nxyWnyMJ7wnFJdoLGpxMx1MiAGAHleV7owmKWqCkTqc5vp3xFDh6xY5nyyFf9HZasj1UczGb3Jw/haCkj5bY0YoZ0TCEpK3CQ/dkBnGYHE59bzY4v5tKcuwKxtQ61bzBRs4aQPq3zCPSV71bTWlCJd+qlJM3+Ez6Bkb0e77e/9WHGP6Zh9BHYtqiF8tUHCRp/O6FJ4wlLGEXB9p85tOJDfJMnMyB5AkGRSSc8TntbM4eyfqYmewkt+9aCaMc9Khp7UwvWqhJ8Rl6P/+CJVG1fiGQtIXpKBOEJ6hNONZQTJrDtm5exmhpx8w3H3laPm98A0mb8HpvFxOLfBRNx5XiGzTbIFYEuUMvfrsQQPosRVz/WZduBzJ/Y8/4dhM/MYOSVF16BC5tV4oc/rWLk3KUnvSqXnd9MzbVs/+IpGncsIGzWKIbP8kCl6fybv/X7Vqo2FTPj79uPLKpst7ZTcSCT7LdvZsCUNEZdfXTEqNMhUZDtoGRTNc529xO2zvxaP84dTcgoP4ZfasBsksjfYqHhQCsNW7cQNO1BbE2ltOSuRBcUirW6HLVXMIGjb2HonD9hNbewd/VnNGT/iLkoC/eoOHyTA4hI8yQoQnGkdaiuUiTrq0okp4igEhh1+wCMPif/nZcTJrDgrgG4DwpHEFTYm5oRVAHMfG4lABUFO8n+z31IlnpSfpvSo1JMsvPD9iUm6nNh2pOLTri9ujiXTc9ewlXvTzthE9X5buP8Zup2lOM79AoGpE0nYvCY/g5J1osW3huDe0QQo24P7VSv9ReWdgUL7/iSia/vO2GRl+qiPax/YhjXfnpll5sGm01i8UMbSL7nY2LSp51RfC7JhUJQUJSzlvKsJTjbm3EPiqVi5XsEjIjjkls6D9ZZ9WEtLqeIZ5g7oUnuhEYpcTolDmY7OPhzCdbaekRzE8akZEKG+TNwqBZ3Q+98b7ubMC/oEQLGxAlI9iLG3RXC5k+VuAUfbYINjU0n+KWN5Kz6hB2v/Ymgf83C3ePcHugh65lBI/Qc+t/P2K3tJ1yMOSgyCbVvGAd3OXp9lZJzwfArPMn1jqSldDM7XnsX0+3vMPiS6/s7LFkv8U6ZgWTPO2GyBFBrRASdAb2hIzG1Nlax4pEUUChQ6j1RKNXow2NO2MKm1QpEXzmMff99guiUSd0aTCaJEjUleyjaupCG7MVYSnahUKpQqDQETshA762mctUyoi9PJX1K12Ikk3/fdS6lSiUQP1xLbHoMB3dFEBSpwujbf7/TF3TCHHH7PLb++wGWPrkZ0dpO2OjHO20XlAKp0+6gfNX7lOc7iRsmJ8wLidFXicegRLZ/8xJjbj3xXC3/4VdTlrmG+OH9U46vL2m1AkOnGQADOYPGUbz0LTlhXkDSrn2K5Q8mUn/FNPyCu/52mZpcKJRqWhurMHgHUrhtMe5R0Yz8XSRWM1hNIp6nSD5DJuopW+dk1Ru3Epg6HZ+wRAIik4407R5v1bwbaMtfj3dKInFXRhCdEvNLmdqjtZOvO7NWDqVSOGW1orPlwmuHOobBK4Apj3yFV/IsxNZ6ds27iqpDu7vs55d2GbveW83ytyrYt82K037+NDvLTm3EHZFUr36fioPZJ9wePfoamrJzEEXpLEd2dnkGqBAtLf0dhqwXWc3NKAQlNmvXbTarxNrXcgia/AdCBqYCoNTosVRU4BIhNErJwGQN/iEnT5hKlcD4BxLQedVRtfV9Ml+bw5rXTlySMT/zR0wFm5jzxgSm/DGIhAwtWp2ATi/0ySIG/eXCeSWncMldbxE082FcDisbHk+lua6s0/aRNzzNhJczcR8wnbyvclnzYd9Xk5GdHb6BApGXjyHr3dtoqa/osr1y32Y0/kHn1VSSM2EzS4jtcsK8UNgsJra8dg3hs0cRGnU06TXWiWz5voWfnt6FLiiZUTe/2PF4TTGW5moku42t/y7s9nm8/JVMuM2fS5+Iwi9jMFqfruXznHYr+z59iLgbMy6o5HgiF3ST7C8ElYpL7vwnxRlzaCzdh5uHb5d9/EJj8bvuSdbUFqHzkcvoXUiGznSnvSGQlY8OIeqaF0ideQ+CUsBpt1K88EWS7xje3yH2mcoikZwFFbTk5eE/7k5sFhNavaG/w5L9Stv++zT22kOYa4LY9D8nCgXUZVdgrSzBmDKd2KteJG7kFUeWgsv86AEspTsYePUo4sd27c8/nZXvVtOcV8jwh7qW0Nv+3atofY0kju7/JtO+dlEkzF9EJo8n8jTL1+h8w2jen40oGi/4u46LhVIlMOE2P0pHTSLro7eo2vwlQSOvp3Ldx+iDAolLv3BGSJtNEuX5DuoKzDQXNGAuzEdw9wWXi5qlr5HlZmD0Tc/2d5iyX2nwpffiFZmKuaECa2M5ksNC9Ow7iMmYfcI+xvRb/8aWf1zFgU9/4sAnHd0P+uh4rng5rVvncw90o63Yi6CoFACa68oo3b2aqq3fYC7awdgnx18Uv5cX9LSSnnLabYiig6WPDyNs6iDSprnJczQvME6HROYPbTQdqCf8khASRmoviP/jwj129nyRh7WqGF1oIh7Rw7G1NWM6sJrYa0ZirrdTu72U6S9vQefuedrjyS48Zfu3kfWv36AQJHTBvkRPGHDC0eE2q4SlzYVCABQK6kqd7PkyF0nS4hYcg7koC7G9GUNMPL6JvqTN8Djvm2LleZhnYNnLc2jZuRhNyGCczRW4HFZ0YXF4xQUQGO9BbJqmy8Rgmaw/WdtdbPq8loasPcTc+CqDx9+ASqMj68d3KFvyKuPnjqS2xEHufzYx/sWN+IXG9nfIsrNMcopsX/B3ype8Rsx1445U98nbYiN/wV6MgwaAAJaqNizVlTibaxH0nuBy4XJJKPWeRF7xFxrzN+FyHGDwZaEER5z/lbGOJc/DPANa7xA8BqczcFoE9YcG0JRXi6X8EFVFe6heqWavMZDkO4ZfUE14sv5TvM9B7v+KcLSZES0mXE4HAC4AhYDa3QO10YjG6IbGqEbvrUFQKbE227C1OrC3WDEXF+ExaBzT/7kHg3cgAI1VRRR/+wyj505GoYDdb3yPZ+ocxONWipBd+CSnyLKXZmCp2INH/CAKvlxN4QJ31J7eONpaibnhb7SU5iCodAyYkIhfxGD8QuNQaY7eeUqiRPmBbRR++RiT/3riKSwXCzlhHmPUba+x4rnJ7P7Xz2gDwvAfFsvYF+Nx2CWyl5spWrCR3e9vpCA8jKRrIomM751lpWQXn5wNVvZ/upbIq1/ALzoVncEHlUZLq8XG9scHA6AyBCAJfoiKQFz401pWg+SwoTZG4xYUjG98EN5XDeq0uoQkSmz78G4Cxw9nwEAVX93QUUtWIdSy7pEkkh/6joQxV/XLa5b1D4/IoVirD6BQRTD+lU9QCEqaa4sJjk4/UqC8vjyfQ5k/4pk+A5VGiyRKVBzIpGjTfBp2fA+IhF065qJOliA3yZ6Qrb2N6qLd7Hz7FtxCg4iYEELccA2eqkYGAAAgAElEQVQ5qy0UfLMJfWQG1tIs3CMiSb4mkvC4i/tDJOseq0Vi3xYrlVuraC8pYejD3xEWPwqruZn2tnps7W34hsaxff7zVC2Zh0KjQ+Xphcsp4nI48RtzM8Ovf+6U6//tXvEJxYte5LK/D0etETiU6yAkWkXhbhs5b/2AZ8o01AY/mnYsxC16KH7pc0iafFufrU4hOz8c3Lmcna/MRNB7YkyehulQJkgOfIcmED3amwGxwgU9qEfuw+wF7a1N7Fv7BRWrP0Dr78HkByOoKZbYv6yGxp3Z6KNHYSvfhVtYOMNviyYwTE6csq5EUWLL/FYqVmzELTId3yEzEe0W6rd/h7ViHwq1DqWbx+EyYnoG3fAyUalTyFv9OY371tJWsBlnczUAiXe9T9L0E68HuG/jd+R99HuGPTiJ6KSurR8tjRIFme1YGuzETzRSecBG/pcbib72ZdJm/r5P3wPZucNus1BXspfG8v3oPf3xi0jEYAxg4Z3+jHh8OpV5JkISDYQNUl7QSfJYcsLsRXabhdV/vxJT/kbcYxPxSfQlINqdrLdWEXfHO1RtX4ygqGD6gwP6O1TZOcbplFj7YS1tRXVEzn6c2l1LaN71E54JSURNDCZmyNGBZKIokbXEROF/fyLh3i9Inngz0NHM2lRbhOSw4Rsah6DsfGEmiRLrP7yXxp0LSf39GGLTutfH3lAtsvyhb9FHpRMx4wFSJt3Suy9edk6xW838/Phw7LWFqH1D0AUG4jS1Y6spxyU6kRw2Eu+ayZAJPZ+neb6TB/30Io1Wz4xnlmJqqmHPik8o+eYJSgQVgs5Aw76NpFz7FBufHoOpOQSD18VxRSY7PYddYuWbpZiLKlCo1RR++xTB45IYecvMExaQzttso3jRVgZc/hSDxx2t+SooBXyDB578RC4XLfmbUBk8iU7pwVdaAcHTJyHaHVRt/kpOmBc4pVKDraaQOe9fcdxCE0m0NYnUVogERcop4VTkd6cHDN6BxI27jvKFLxB36xQszQ4a920n862bcYtIpXivnaQxuv4OU3YOMLeJrJ5XQHtJEbgkIqaOY/gcwymbuAoW5hJz02sMmXpbj861e/lHOBqKCZ8zDmu7C6VKQqkESQJrO1jNLgxeCtw9jp67NN9JW6OTcb/xxdTqYukjy2ioKjx1Ypad15RqNfoBiax4MQdj/ACCBnsQmaRB7y7g4a3E4ySrnsiOkhNmD/kERxE09V4aDu5gyj1BiKIni5/KwVpziIZDWpAT5kWvqtjJ5tczsdeWIrgZGfKHKd1aPswzNhxTzaEeny8gdjhFOg9KF2+idKGIS3LikkQABJ0Bpc4dp6kJXVA4QSMH4pJclP28FclmwtN/DmExKrxSktm/+jPG3Px8j88vO39M/+tGSnI3UJ27hgOL15Hzfh7hs8eRcbnHBTWvsq/ICfMMJEy9k3VPfYTjjmmoNQLJNyeS+fL3+ERduDVJZd2zZ5OFvP9biWQzowmMZOwjw7s9GCw0zUj+wiVUjb6a8vztDJ16O4LqxM+tOrSbpsoCEsdeQ2hsOlf9X9fC8sfa+NlfqFz0V4pK9qGPSGXYo9+z/W+XEhyuxGaVaMnLJ/6yeT1+vbLzi0bnTuywGcQOmwFATUkeme/cxpKcSiY9FtdrCzJfqOSEeQb8BsSh9h1ASZ5ITKrAwGQNfu9ec9KFXGUXB1GUOLh4H4ZBg/GIMDJ0jlePFiWPSdOQ81ExG58ei8tmorVoJzHjb8LDbwCi3YrW3YinTzANVYVseXkGWv9IEsde061jO9rqUPtF4DvsSkbf+gpVhdm4RCc/v3wAl0vELXwIA+LkC76LTWBEIrNe3syat29n5Su7mP5UEjo3uZrZycijZM/Qir9fjX+Sg6HT5ZUfZL1n/ReN1Gzai7Pp6Io5gt4DydKG26CxTHh8ASv/MhbP2ECasvMY/vgiwhNHn9G52lubOJS9DNFmIWrYDDx9LrxFtGXdIzlFVr95C7aGPKbNTUCnv7juNLs7Svbield6kS5gIKaaE6zcKpP9CsOv9CJgVCJKT3+UnkEoPQKQrGZUPuGo3IysfHIk3oPDmHrvAPzHpLL1mTHkrvv6jM7l5ulN0vgbGDLtdjlZXuQElZJJD3yOQuVL9s+t/R3OETVlTnLWW6gqcSI6+3+Rd7lJ9gy5+UfQtH9bf4chO4+JTglTi4u2ZheiA0SnC8npYtB4I2NumkTBTgcVWU0oNAJ6oxqdUYn71FR8glWsfLeahqxdhMyZS+zwWf39UmQXAEGlQmMMRKUz9XcoNDdIbPqgCPOhg7hFpLJ372r8J45n3G1BaLX9d58nJ8wzZG9rQG2Q3z7Zmfn+8SysFcUo9QaUBl8QndjrigDwHz+eKX8IISFDS0JGEGaTRFGOjZqcVop+KsLZ2oTf6JuY9vqXePoE9fMrkV1I7C3VuHmfvPTi2VKaa8PRJjL73WLUWne+vV5Je0UT3/9uK+4xg5n62KB+6Wu96Powb3lcnmcmk8lkvemrl0f06vF2LDVRu0tk+l9+AmD79/+kZuv/mPDED6z8y1iiLx9M2qTeq0gkV/qRyWQy2TmlqsTJrm8qaC+vwCc1lsSp3gQO6JqGdJ4qLJV7sVlMaPUGhl/5CFz5CACJv32DXa/Oxjf0WsLjzm4Kkwf9yGQymaxP1VeJbFvcyoYXV+ERNZO0P3yO5Ixl3XOr+Pb+jVQWiZ32TxqtwzMunJ8eSCD75w+RnM4j26JTp6AJjGH3lwWYTWd3INBF3ST7xauFZ3SM+ooC1s7N4Kr3p6PWyNccss5MrRLLn81E4xdDQNplRA6biW9IzAn3zVr8Lwq/eBiluxFB48bgO95lUIY8iEfWfUV71rFv/jPYag+i8fXHVluF2N6KoDs85c3lQqHWog2Ixi04HqXeC9HShOSwolCqj/xp2PYNSXeOI3Fk96qV3fjk0UGPz16ZSmNpO3ovDd7BGvxClBi8BGxWiS1fN1GzYTsecWMYNPsRolImHHme5BRZ/39/QlAdYOLt/l3OcXCXg9z5eag9I5j25I8Ih0tLig4Ha9+5k7aDGwi5JJFBYzzwDTzz32K5SbYP5S19H78R6XKylHXhdEisfm0fngnTGH/3u0e+4Ccz9LL7SJ52J20NFXj4hKDRXXwrRcjOTHtrE+vmXYulYi8Rlw4jbfq4w79JiVjMEpZ2V0cToqDAYXNRX+agqawCp60EXaAapUZAEiVcTheS6CL5d+NIGHFmpT3zF5bhNiCZ5sJaiuv2Ya8rBYUChUqNIXY0k17diU9g5JH9nXYbe1Z/RsK464gadwPZ7/0GUfTtUms5JlVNRGIyix7bSt6G+SRNuBHoqIs78f5PyN+6kIot37Dyx6V4D0lhxK3BfVpARk6YZ6Bh+7cM+2NGf4chOwfl77AjWhWM//2/Tpssf6HR6k96ByqTnYygFLDVlxE8IY2Myzw6bdO7C+jdO+/vH6KEM0yIpzP775md/i2JEqaWGtpbGwiKTOq0rabiIFtemYWjoQS1zkDiuOvJcfNnxZsVTLw3tMu0EbVGIOnWoez+4EHcfUOISh4PdLz+hDFXkTDmKtpbm9j2xVyWPvwVXqlDiJkUQGSiqtfX85Rvkc6AR8woKnLNJ93ucrmwNFgwV5lxmB30pNkbwFxj5sA3B9j/5X4OLjyIo93xa0OWnSWiw4XKzUhTbTFWc0t/hyO7gOncjSTd9hb12SX9HUoXglLA0ye4S7IE2P6vW/EfGkLQlHE0HcpCUApMe2ENksOLn5/bQ12l2OU5cUO1xN8ymp3zrmbzF8906tOEjiIcE//4PtPe2Id76FSyP8xi8VM5JzzWryHfYZ5A6c+vU7XxEzhJonNz2rGXtnHILYzo2dFHHndJLkpWlFCdWY2+yYabQkGzy4UQamDg5QPxivHq1vkr1lUwcHstAwWBraJIQ4iBoAx5vt35YGC6lvyvDrDqT7EAXDNf7PadpkzWUw6bGZXb2W/GNzX/isE2Lhc+0V5U721Hrev4jdXqDUz980Iyv/sba56ZR8x1YxkyWd/pDjHlEj2hcVPY+K+FLPzDfxC0evTB8Ux46Eu0bh132J6+IYy64Wmka+ay7X8vs/al9xn821Ekjjj9akHdISfME6je9Dl/trXiqTh5H3CrRsOrWTVHEqbkkMj7PI+AvCYeV6kI1R1t+siqtbHgvRwqUnyJvzEepebUbeySKDFQEBilVHJIkqg+B0pCyU7MYZdobYKmKifV+WZM5SYUwtEvuZwsZX1J6+6FvaUZUZR6vfnxZPJ32Mj5aBN3+biffucTCBp1I/nz5+E0tTDh5beOPC6olIy8/inKhkxh5zu3Ub5JR9otsYREC0dem2+gwKxnE6kpdVFdYOHAx4vZtfQDhs95qNPKPoJKxagbn2FtSx2Hlm0hccSgX/eiD5MT5nFEmxlXaw0DdKdu6++40T96B7rvv/uIymviNrUalUJBk8tFtSQRJAgMVSpJEQQ+2dPIIfdDxF4T26evQda3asqdbHojC1t557Ur/SfehV/yUPxSFPhFJhMaK6/+IetbEYPHkeOEoj1OYlI1fXKOukqR3KWNmCtbsFZV4nI6SfzdB/D9wz0+VsWBHXiFxFDtH4ck5uNmDOyyT1j8CELm5ZC1+E22vPoqYnsL2qAwRv0pneAIFSqVQGg0+IW6YW+/lKq17/PDincInvA7kqbfhcEr4MixXKIN74Suo2/PlJwwj2OpL8b/8J1lnijy3+Payn8huVwYon0AaC1pRcyp5zatFpVCwSqnk0UKCV24B9aSNq5EYIJKxS1qNc9ursQ8NgT3oDO7OusuSZSwNlhxWpzoffWoDeoj21yiC0u9BUe7A42nBp2PDsUp7qZlR23/qY2Dn//Y5XGvEdcz8d4P+yEi2cVMqVajC4qhrdEJ9H7CFEWJHV+W4XIFEzbxTvwikvAPjUOpVp9Rwtz54d3Y6w7hO2okDkME6+Zdx/SnliCoOqcipVpNxlWPknHVo9itZvYs/w8bXnwWY3ISg6YGEpWoRqsVGHWlEXHOMIr3Ojm4cinLHpxH2Jy5DJvzIIJKhb2xDP+4UzfHOu3dH2MiJ8zjWOuOJsxylwvVyBsJm3Z/p322f/MKonUPk+4MA6BiQwXTVCpUCgUHJYlF7kpS7x+KzkeHpc7Cd2/sJF7suNscq1SSs72agZedWYk+c7WZvM/ysNZ3rJSi1CoJHBZIdWY1kkNC0AoYI420Fbfi1e5ED9S5XGjjvIm5IoaGvAbK1pThbXbiDjS5XLiijUTPjsYz0vOM37eLxS/JMmjmI7gHRhMUP4rgyOQuX3iZ7Gxw2q2YCrYSfvu03j+2U2LthzVYq1uZ8MyCXqlbHDDyOkwlS5j8O2+cdi8WPrKJkryNneZmHk+jc2fonPsJT5vKwc3fs/31f6J6fNKRplpTs4vs97egUKtReRop/f55arZ+w6V/3YTKM4Cm0krgxP28NaUiqx//ptvxy9/y41jrigk7nDBrXC7cAmNwmJuQbO2oPf3RBwwkcsJvyXx1NvsyA0gYoaX5YDPJh/uttogiAyZEovPpaNLV++vxTwsgN7OWIEEgVhDYVNp2RrGZq8zseWc3NzkgTamm2eXizXYHdTl1XOGAS5RqllqdHMhr5F61mhBtx5WV6HKxprCN7/6+nVhB4CmViqDD2ySXi62lZr5/O5vI3yXjk+Dza9/CC5r/hIlojGMY85sX+zsUmYyyfVvQ+AbhG9T7cw+zl7djKmlmygsbcPPsnYLsnkExNO6xAKDSKMAFek+/bj3XPywe/+vnsl2jYetr/yDwklQuudmHzR+X4512ObFT7qR89yrKFv8NS8kubBYTcTPuY8ebN7C03saMhwd0OaalXQKF4qQDPI8nJ8zjHNskWyNJlC18gWCFAh1QCbQHxBA48gbib3ub/O+eJyomAfdWO36H+zzLJAnf4+7UjAONFGytZgoQLgiYyk24JBcKofvNoKZKE7nv7uY3DgXDlEoaXC7ecjgwXhaFttGKcmstGoWCQYLATJUKK7BHFPFRKAgVBKaoVBgUCoYKAi7goNQxkChaoWC0SoWXQsFHiwrlhHkcU7PE1vnVjPttEFXFIo07djPm2Tf7OyyZDFNTDXs+eQC/tKheP7bTKVG6dDeD7/ig15IlgFKjQ3J0TJMTRQlnSy1e/mE9OsbwKx9h4OirWf14OpVjpiKolDRs/w5jRAojrnmc4MSxNFcW4ObpjZvncHz+toOf7g4hZ70vKZfoOx0rLFZJdmgM1vKCbp1bTpjHsdYVEXA4YWYoldynVKI7pn+vsKmEFYv/Sk1QPI7aMtrrI/E/ZlRktctFWGDn23+dj47Gw1cw7goFHlYRW4sNnXf3JhGbKkyU/FTEbU4F6UoldZLEGw4H3lcMJHRcKAULjv5nxyuVZIoiX7mc6KK9MJWbuN7uZKxKxUilklJJ4gOHA/sAA6JdJKrOyv0aDYlKJZqadiz1FvR++lNEc3Eo3GOnbGcLw6/2pW7NOhas6Xg88f6vCIpK7t/gZDJg//r5aHz0jL6+97tSWptcOJpqCYxK6dXjat08sdVVs+ZjP2o3bEehdkPnbuzxcXwCIwm/4iky3/8Ps19KZucyIxVr/0PapfcQnjia8MTRR/Z18/Qm7s4PKPj2aSKSR3aqBKRUCqT/LpXNz3UvYcpj3o9jqS8+kgAvUaloB7JFkRxRxO5yMVAQuEej4ZKafHz8Ashd0XKkddzsciHpVah0na9D1G5q2o/5t16hQLR1f0Jt7ZYq7jicLGsliXl2Oz5XxxA6LrTLvlaXi89FJ4MfGkrKPSkMum4Qu6Sj01IWOJ14XzWQ9IfSGfboMPYpwXI4mQcIAtYma7fjupDt/k8WVT+vxN0gEDJzBgAxt79H0vgb+jkymaxDe0MZboGGPplO4uOvxCt9KNs+7vnAnlMJix9N0MR7MJWJjH1hI5P+nnXGx/IOG4ytrpLV71Vx6OuVBGRcfdJ9U6bcik/6Vaz5+x5s1s7T9ExNJx7YeSLyHeYxRLsFqbkK78P9ewscDtao3fCIH41kb8d6YBM3qpSMUCq5Uq1ml81GcWE5enXHHagTUKi6NrMKGgH7MW3kAh1FDrpDDfxerWaIUkmLy8U8u52A6wYRPDL4hPvnSxKGGC/cgztG4ao91LQePrfocnEAidHDOzrvFUoFanc1Fgvo6fjjtPZuZYzzUdaSRpzmjn7mhmqR2AmeVK12J33WPf0cmUx2VPjw2WT989+scQ1l9M2+aHW9lzi3/9ROa14+yX98steOCR1zLUff9Gy39jU11eC0W3HzCkCj1eO0W2muK6W1tozWulLKV3+Id3oSzbvymPluMQbvrlNUjjX29n+y8rVS1r5fyfQHj/ZnDh6lI6ebvSxywjyGoFSjMgbxfGsNSqAtIp30uz9F7d7Rht9edYAv5s0mxuXEV6FgvFJJDRLOwwlJA0j2rkUGJKeE+phmXRugVHevk36WSoX74ed6AF4KxbHTP7uod7lwCzraJKxQKI7sbgUEreq0hRMudhWbixFbGwBYOfdnXE4HEdfKg3xk55bI5PEoHv2BXR/cRXGuJ3HDeqeajeiUKF60hRFPLGFAXP/NJV713CSczVWIljYUKjUu0YnK0xe1ly+2uiokUyPth5TE3v7uaZMldBQRGfvHf/PT3aFYLSHo9B0XGPK0kjOkUKpIe3ojtqYKJIeVqJAEFMLR5OIWPAjfYVeRte0rpqlUDBQEcDppU3S88XqFAskq4hJdKJTHJMhmG7/0MoguF0240Hp178PtrlBQLkk0u1wkKZX8Rq3m1cWF+CT6oDV2PYYDEFQnvtJ0AcjTLU9rzkvplB1MYdPzi5EsbbjFjmHY5b3bNCWT9YaIpHEcCE+mrc4G/LqEabNKZH7fTGNuJSpjUJ8nS1NzLVZzMz6BAztV6QGwmluwVeQRMmsa8RONGLwEtDpQqQRMLRI/PVRO3N2f4BkUSUTiuG6f083DC314MgU77CSP6xhDYmqTE+YZaSnMpHLt/+ESnSi1bkTMfgKdb+cRXBrPwCN9fgAaYxClLdVH/h2oUGCuNmMINRx5zFxpJvRwv2iVy4U2wA1B3b3mk0JJ4h3RgUZ08ZwgECoITHaIbPuugMG3D5YLDvQRD6MCl6OjP3fCnxfIJe5k5yxrTSHOgcE01Ig014iU7WrFGKInKEZLcGT3fuKddhcr/1mIy+XDwMufJzp1yq+KqbGqiH0rPsLaWAYuCUGtJ27mvWj1BqoPZlG+/nPa8tcjaPRINgupD31HY0kOUcNn4R8Wj87dyIjnN3JozSese+57tIEh+CaHEpbqgVeAAJJI7Ig5ZzSCN2zyPRxa+k+Sx3UM3lvx5OpuP1dOmMeo2vAJ6bnLiBEEdksSRQZfBl79Qqd9Wot2EHw4+VW6XHhEDcVUuptaUy0BgkCMIFCZ39QpYbaWtJJxOLGVuVxH+he7Y6soEjgzEkuthcVZdVyrVnOpSkX23kbqd9fjn9p7ZZ9kR7kbFegGxKIyhHQqtSWTnWuc5mYqVpsoWViHQqXGb+RNVO88xIH/LmP0k9MIi+2o8iWKEvu22mksbkOhVCIIClySC5fkouVQAy6XF9OeWoRK0/VO1Wm3cmjXKporDmBtqsTeXMEtx2xfcGcQktWMoPdAFxRDe2kOvsPS8Yp0A4UCS2MdW16cjEKlRusfRNCIMMb/8VLcPQV+eDKHna8cLrygEPAPiwcgYvAYIgaPwWn/Fwe2LaI2dzWZby5C4+uLQq1Dozuz0fxJE27m0IIXWPluNcOuD8A7LYm6tWu69Vw5YR5D5xuGAAxRKglRKPjr5v/Smn45nlFDcYlOarZ+jbB/LemHBwVliiJ+qbMRVFr27lpEgCBwiVLJaytLMIQZ8IrxojarFvOuOoaqO8pWZYkivom+PYpLpVMRPTuaDXvqGSlKhAkCv1GreeO7Arxiu7cCiqxn1BqBqU+n8eP9i5GcYpcmI5nsXHH5vw4iKAUkpxNJcqLSdDQ1bv3mr+z64lsCnozD3Opi41v7cZpseKdehsvuAJcICiUKQYVnrDvDrnqsS7KsKNjJ3gUv0Zq7Cm3QANxC/dEaNXgP1MLRhjUmvTABnRu0NkjUFjmIGT4Dd8NxrTK/mX7C+IfcksCe/6lRuUcxdPZ9XbarNFoSx11L4rhrcdpfZ/fS/0Ol9zjyOntKqVYz+a+b2Dn/BZY99jWCpvslBeWEeYzAkTeycdV7zHC58BcEbhVFvn59DlZjEKKtnShbG/doNKgO9ysWatwZmjQVrXcIK7O+Z7RSSbgg8EfRxf+9l4NVpyTUJnG/Wo27QkGBJFFgUDE8pXuVLY6lNqiJnBPNf78p4HGNhmhBYKzVwd5FhSi18o95XzAYlQhaN5obyjqtFi+TnUt+6S4QVCqEY37S0y9/mFV561nw2/kABM14kDG/fbWjDuxJiA4H1SW5NJbupeHAVhq2fsmAGSMZdesMvPyP+53ZdPSvvoEdMbh7CN1uBv5FZKKaHKcL/5Tpp4wNQKXRMXTO/afcpzs8fYKZ8If3aL3uacwtdWx4PLVbz5MT5jH0/pEYMq7ji8z53K5Wk65UkioINFubUALGw3eWTS4XHzocRNz8KkqNHs+oYdSMuJ5PM7/hVrWaeKWSfwgCZhd4aDuaYgsliU8cDmJuSjzjUapBw4PYlVnNhrJ2xqtUXK5SsXtHLa1+fbOK+sVOFCVUnt40VhbICVN23tFo9aiOKQrQlLOM0rw5qHUGrG2NePiH4T8gvuPuVJTYu/4rDn79JAoBdMHBuId4MvmlqX1Sdu94gy4fyK43HyBu7HW9WlnodDx9Q/D0Den2/nLCPE70NS9SYG3jxd0/ccnhZbl8FQrsdCS9TaLITpeC4CueITDj2mOe9xKFLolnt33DFWo1qYKAx+EFpFc4nWzQC8TenIhf8unvLgWVwP7DhRKqJAnl4QFCCkHBoGsH8cNrWbicHQNeoxQKdtRZOCgICE4nRZLUaUCRoBZoBdY5ndgA4bi5WoJaYJvJgRsdU1J8ujkY6WKQ+YMJlygQEiMv0yU7/1QW7qJ170rm/N/16Nxg1woLu967FUGjQ6l3w95Qg2S3oQ8bjNjeiktsJ/n2Yb02PaUnIhLU7EKBzu3cXgBCTpjHUWr0xN/xAc35G1if9QM/5K3G3laLUuOG1i+c/2fvvsOiOLcHjn+3sLssvffeBAUFUVTs2DXGxCSaburNTS8mxvTEFGN67/emR5NomiaWqGjsvaAIqICA9N627+8P8tNwrSgwC7yf5/F5ZGd25mw9OzPve473gOkkDLgClbPXKfeLvvYN6lOv5+dfX+TLvJ3ITHoUTnZ49PUi+ZJw7BzOfrrh/wWNDKLQTk6ByYpSqyToH0nWwdeBsBviWJ9d3XKDDGIjXCjNr6fYaEGukhP4j6IGDv4OuEwOY11Vy4jP6LjWtWKjpkexfV8FACpnFa6RPfuaaGO9mSO7DJTsraIm4xAj5v2F1qlnPydC12Mxmdn10e0ET0nFwanlR3D/CQ70n/DPKRi9qSwzU5ZnApkPUYl2KM8wJa2j7V3VgNo30ubHCoiEeQauMcNwjWl5c1mt1vOevuEUmkT8vT+ia6pnxZxkrDV6KrNkFNZVkzLDE0eXc78h7b3siZwWecblngmeeP7PdVDvfqcfySmTyQgafebixm7RbrhFd94pEFt2cKue/R+vQRvcF4++k0i8+nPc/dq/sLUgdLQdv70JNDFgytlH5Ht4K/DwljZJpf+3iso9hxnwwA+SxnE+RMI8Dxcy11GjdWLya3spzcugrqKA4h2/suzeH7Dz8CXl7v4EhNv2L6mewmy2cHCznvJDdVRs3UfS/T8SnpgmdViCcEGa6mvY+ePzlKZ/xrAnxqGQ6IixLRqLqvAZdmOXmL5l+89mF6ZUaQiITiZ2yGWMvve/THj3CF4p17Lh2QatVbAAACAASURBVGWsfKvolCLAQufS6y2seCWPnJ+ykKv7k/rUGpEshS6prqqEDV8+xh/3RqEr30ra8+PwD+saP8ojJwZRvvkrVj3Ym5rSfKnDOStxhNmJHN18GHLtszRMvotNH/yL5fMOMObR3ieuMQidx2SysOX7KowNMiYu2IlKLVqaCV3TzmUfkLtoLm6J/Uh9bDQBXSRR/r+Y/mpi+iez9j8VpM8bh+/wWYSlTMUnpLfUoZ1CJEwJOLp6M+bhH1m14HJ2/17L0BniGmJnOpJhYN8XGcjtnEm+8z8iWQpd1t6V/yV/yTMMf3ocfiFd++t8+A3uZO92omDLCjY+9yZylQMylT1YzGhD+uIeM4yo1Ctwdve9oO2bDHr2/PEhVZnp6MvzUHsGY+8bg5N/zHlvo2s/w12YXKkgfNy/OfjF3ZivSO6QnnZCa2aThQ3fVlO6YQ/hV71A37E3ixqxQpfVVF9Dzhf3kvrUxC6fLAEUSjmxA9TEDvDHbPal6LAFs9mKXC6j9EgTlfu/JXfRXLRh/fEfej39xt9y3tsuytnJzvdnodQq8R8chItvJHXlBuqLd1O45o/z3k7Xf5a7sIjEcRz8Qkn2LiOxA07OfTKZLDTWW1t1Bhcujtls4c/3jtNc0szYV3a3abKyINiivD0rsQ+OIDCi+32NKxRygmNO/pgN6eUEk53QNYWQtVVH9vdP4ugRSGTy6cvt/VNTXRVbX55CyJRk+k90+MfBiR3ggF7nwZJZGecVV/d7prsQuUJO0IT7yfntI6KTYk68kDuWNpD7w0rc+qegdtUQN84Nn0DxUl2MbT830FRYzfgXN6O2dzz3HQTBxuWt+gD/1GCpw+hUGq2MvqPsMRsHkLnoccL7jTnn3M3MvxZh7+/PwClOp13elsbb4nyUxBLG3oSxtoGDmw0nbqs5Uov/5EdwCpuMhVjSn1rFph/rJYyyazObLBxbup5BD34vkqXQLZTkZdBcsJ/4Udpzr9wNJYy2R25nYsksN5Y9PoyNXz1BwaGtWMynzjwIiEul+Xgh5tMsaytx2CIxpUpNr+tfJfPzu4kekIpaLcc7wZ2SnbsY/8TvAOQkTiTzuwfhikSJo+2aSo5ZUTi44x0UK3UogtAuzEYDCntH1OqeecyjVMqZ+lwctdUxFGQaKd2/mR1vfsF2kxnXvhOJnXwvfuF9AagrP4a5uQGjXobiIn9fiIRpA3oNupTc5e+y87d6hlzhgleoirzfD51YHhCTwt6KEiqKzXj6df/rmk1NVrK2NFNb2IipwYzSXoHWSw0yGU2Vegw1BpwCHYhPczqvykkyGVgt5k6IXBA6h29oPMaachrrLT16WpqLmwKXIQr6DNFgNntReszK4U1FbJqXhkNYMhqvMCq3fk/yAxPRaNtegOZ/9dxn2sYk3vgahcs3U1tpxtNfgbGmhIJDWwHQOrvjNfR6di8plTjKi2M2WWisP/dpkb8+Pc6Rn/aBLAlt4CRkmkFU57pQfcQJZIk4h0+jJlfN7w+u5PBeA3qdha2/1fHz3L2s/qiUxobW+/AJlmE16SjNP9hRD00QOpXCzg6NfwyF2UapQ7EZCoUc/zAFw691Z+rbY/CIUyKTHSbtxTFE9jv/npdnI44wbYRvWAIeg65m4yc7GP7vMGKuH8X2V6chf/R3AiIT6X/l46x8OJHMbR7EDuz8bgIX6niumex11dRkFqIrOgIyOXYefrjEhGCo02OoqSN0XBT9Rp2cC+nXzx1DjZxhN718li3fz6HNv7DzrVnslIFD+EAiLn2Wwo0LWb3gEOMf63XiYr5CIcc9KYHstV/gM+ts2xSErmHvyv9iqCrAwz9K6lBsksZezqBp7d/5RCRMGzLw2hfY9PG/WfbAcqKuHkHE5QPZ8sI4kh9aQkifYfS5/RMyPr6ZqMRRKLtAG649a5vJ/mYdnqnX0+vqBwmOS8VOpSU/cwMlBzfi3ssfO40DWd/OoXinL86hLugqDVRs2Y6luZ6q0ryz9qHsNfhSvEK2YzYbT1yf7DVoGqsWTOfPN3IZ/3DYie4LMWnubHn9a0zXPHdKV3lB6CosZgt7ln9E7uKnGDJnVI+4RGNLRMK0IVonV8Y89B1FObvY+vIUgsYnEXvDUHa8Ph3FnGVED5zMwS+c2LeumaQxZ+9CIKWGWgvbfqykckcmA+f8TlDsoFbLw+JHEBY/4sTf4Unj2PXLm+iqS1HYa+lz570YGmtx9Tz3kHkP/9ZdXeRKJWmzF/HH3MHsSz/5PAVF2bHXz49tP7zIkGufbYdHKQidq7L4CJvfuBpzcyVD5ozqlvMvbZ14xm1QQFQSw59NZ90TqQx5dBQBY/qz9aUJqOdtJPne79g6fxIKZSp9R9rWkPLKUjO7F5dRuX0XLgnjGTN/23kVCNA4uDDkmqfbLQ6lSkPwhPsoTH+LpDEnR8YOvCWCtU++Rd3Ef19weS1BkEJB5ha2v3YZ/qOTSLk8XLK+lT2dSJg2yGTQoVSqwWrBTiUj5TIX/sgMImvNFwy98SWsc5ay/ZVLqTmWQECCC5H9zq8xdUc6lmVi6xvr8UiZwZjXvsDVW9oJ1bXH9uMU8j89QwMUKJ1cqa8sEglTsHn71nxNweqPMFYVYaovp9esMfQdKeoeS0kkTBtjMVtYck3LhyJo6kR8glquUfS5IpLtb31GQco0gnqloHkmnX1L5rPrvWVY7xxBVGL7jAK7EHqdha1vrCdy5sskjLlRsjj+yVhXhtbn1Le3c1Qo+Vt/JiCqvwRRCcL52fjlXMo2fkXMzIF4BPrj4SdrU0UaoWOIV8DGyBVyet3xOXIHVzwiNJhNLVMkwnqriJg+iJ3vXEtTXRVeQb1Iu+9z/Mfdw5G1Je1SxeJCHdqiQ+0TZTPJEiBizK0Ur9l1yhST6LHelKz7LxaT6ZT7GAyGU24ThM6Wf2AjJWs/Yey8ocQP0+AfphDJ0kaIV8EGJYy5kbhbPyZzUSY/3b+BtZ+X0dhgod8YLa5xoaycO5CSvJZiwb3H30ZTQSWrPyhB1yxN0jTqLKjdAyTZ95mExo/AOS6Nrd9VtLo9JEaJ0smJnJ0tHQr2pH/H8mfH8sczafx8jZrvr5Rz/MgeKUIWBMxGI3s//TfhVwzBxV18PdsamdVqPf+VZTLrVT+e//q26LpHIk78/+sFRySM5NwsZgu7//iAY78twGox0vvGQcQN0rDpxxrKdpZxySs7gJY2Pxvfv4W6Q+uIu2k4fYZozmv735q/68jwBaHbmqFu/+/B3L1ryfj6YZT2MPHxSJtv+Xf1Y1tP/P+7F1MkjOTifTfzO6xW6zlLAdn2K9LDyRVy+k+5i8s+yqf3rR9x8KsdrHqvmIFTXTCU51JWkAm0TEcZO2cxjtGpNFWLyh+C0NUcP7Kbna9PJ2ikPxMfs/1k2VOJV6WLiEm5hImv76d6935qKi24J/cl8/f3Wq1jqqtAX2vAZJTueqYgCG2398vZBE8ZQv/xjijElBGbJUbJdiEaBxe0IX0pOWogeboXKx9bSP7QqwnpnQpA0m3vs/Pjf/HzAwcIn9a/TcUNrlFcfcFx7VzRQNGmGia/sP6Ct9GRLGYLK1+YhMq1CY2rPVUZ1UxZsK3VOns3LMHOaiZu2JUSRSnYukX6iy/efSZNBfsZ/K9RHbZ9oX2InzJdjMVkABm4eikIu2wI+764/8Qyv/C+TJm/hb7/+oIjP2zl4BZdp8RUuC6XsAl3d8q+LoRcIaff9a9gNYdQmwt9b3qr1fKS3P1kvTmdgx/fQl3lcYmiFHqqiqJsrCYjji4dl5CF9iESZhfj2Xc8Bz5dy/qvq0gYraU5f29L41TTyfZV4YlpJNz9FRmfriNrp75D49E1WWkuPEJE8qQO2b7Z0Exz2RGay460/Fi4QL5h8aQ98BUTn15BUOzgVsuOH9yI2j8Cp5hYMv744GJDFoTzVlaQyZZ3b8YrNbnH9rbsSsQp2S4m5aoniBw6k20f/5tlj+9EExjG9tcvIytmBKPv/wb534MFIhPHYrnzS/Z+eBPlaQMYfIXTRQ8ksFqtFKwpoLG4EbmdnOAxwZSWKlF7h6G2dzxl3fq8nVTuW0Fz6WHkKg0azxB8h1yP5jynoFjMRnY/n4q/sggrUOEwgL4PLb2ox3A6CWNnUXkwnbrMVVgMtlujV+heasoL2DhvPN6psaRMd5c6HOE8iJ80XZCHfyTjn1xBQNo9GKsqibt+EE1Fu9ny3VOt1oseOJlRL2ylYncJf75z/KKLGzQUNWDenMOTMVkMMR6leHMxukYLSofWH3aL2UjOtw/S8NVkrpS9zCv9F/Jsry+4xPAimQtSyV+2ALO+6Zz7MzfXY284zl9zylh+XzlNJdkXFf+ZKFUaUm5+HXNdFTK5+EgIHc9itvDXgsvxTo1j+LUeqDXidGxXII4wuyi5Qk7y1Htx9o1g37vXkvLIWLa88BolqTPwDYs/sZ67Xxjj5q1j5ZMj+OurakbO8rig/VmtVo6tPsbNA5q4PKmZigY5W+ssGJosKOydTqxnMRnI/Ow2ksxL+eiBKuz/UbFvWiLcktrIc0vnk/7BJuLv+RGZomPeggZdE6V5+ynPz0AhV6Bx9kChVKPWOhEQNRC5snVbJGcPf4a/uh+X8ygWLwgX69DmJVgMNaTOjDj3yoLNEAmzi4seOJnC5MvY9vpvqDxDcPYIPGUdtb0jIx5byurHU9nqriJlqtNptnQqQ4OBkm0lNJU2UXuklliHOm67qrHVOvpGEwrNyUatJRu/It64lP/cXIVSAWsOqfkw3RGzFW5ObWRygo4Prqvm2k//Im/txwSOufPinoB/aKgpY9/Sdynb9A2GsqMA2IdGYefkisVgwmo2Y6yrZbvJjFfKVSRePgdHN+8T9/cN7dNusQjC2RSu/4rAUb3EfMsuRiTMbqD/Nc9z0DuchIl3onV2O+06zu5+DHl0KRueHkZg7zQCzqOXXv7yfOJrjzA2Tk/foQbi/E+tv1pf1Iy9XxLQcnRZuOpd3r+pHqUCvt6i5anVsYROewq5Us29Pz/D8Zr93Da8kacvqWPcRx/hP+o25IoL67airy6ifOcvWMwGassLyV37OY5hoXiHuVBn8cUks8NYlIdKqyU4MQClWolDgD/4+ZC1difpTwzC2SsQa2M5MqUav/5TCRl/Lyonz3PuWxAulMVsQV9dhJ321B+3gm0TCbMbcPEMYPDMJ8+5nk9IHAGTZrNvyVICHg455/omnYkpCTouT2o+4zp2jkr01S1TMerzd9PLuZD4QCNmC7y20om4O7/EISAOALV7AP/5NI3bhjcS7WMi0rGYpqJMHIMTzvORntRUepiD717BzN55uNhbKEPOtw5aFOXHuWd0A0TDB+mOeEa5cVNqI3F+R2jUy1ibpebHdQV4eTvg4tnIDYN209u/Jd70rFy+mfcVQZe9iO/gC5+XKghnYjLo2fjfh7CaG4gffn4lLAXbIRJmDxM19ErWLltAWVEgnKMlpL2XPfP/cOKD9JMjYJ+7tJbBESend0QOceaveT9RW/EMtYc3M/HvZdmlSurto08kSwCHgN4cMPhQWF1JoJuZ+AAjuwr3nzFhmo2ty/xZLGbKCjJxVCo48N6VLJiUwxX9mymulTPjI0+azErc1SYemVAPQGaJkndm1uCoOVn3c3KCjqGRBhZub+Kzh6px+EdN0LRYPTenNnLJuw/iEBCPU7A4RSu0r4y1X1N7aBWjHk0UTaC7IJEwexgP/0j8xt7Dru/WwANnXzdkXAhN/bywmlqSyuGfDlOvq2u1jl+IEs9B/di18Bl87aFfn5MJU+sX02pdmUyGY3BfMouzCXQzkxBoYH3BPuDaE+uUFxwia/Xn1Gato7kgg9jAk9M8LEYdm54ehr+PF69NzefypGaO18i58kNPLIMeZ+zo29k5J/zE+u9fW82eYyreXetIcY2CqwY0cceIRqYlNjOhTzNl9Qoe/8mJrFIlvXxNPHtpLWFeZm5KbeClN69Em3I9AfEjCYwZdGK6jiBcjMbyfFzjgnD1EO+nrki8aj2Qsa4Mufrcv5VkMhkOPg44BjjiGOCI0v7094ka4U7tgdWYGqvxcGyZunK0XInWJ/KUdZUObtQ0tbzt/FwsGOrKTyzbsuh51j8xBFPjdqIu8Wf8S+Na3Vcul+HtrOataXlcntRMQZWCy9/3RDb8eYLH33vKvnbnq7j66yiOD/wE9TW/M39dBDmlLY/BaJYx+R0vtvk+w+a6aL476MVH61qOpAdHGNBozejK/2THa5ex/LkJlOTuO20PTUFoC7O+gebSekyGrt31qacSCbOHKT66l4qt3xOQ3H4Tpf3D5FgMDTRUHMNV25IwG/RyFPbOp6yr1LpS29zytrNXWbEYTl4fLVrxNkOfHEvav3yIHajGwan121OBhQ+vrWJqPx15lQqmf+CBXdrLBIy67bRxfb9DS9CE2XgPuBznsGQcg/uSX9kyneTPTA3WkHEEjb+PYY/8htk+gLyKlmWejhbs7a2MucOXyCsHoy/NZPNLE/ntwXgyN/108U+Y0GMlX/kEFpMzK1/PlToU4QKIhNnD1JbkAlbCE1TnXPd8KZRyku4ajr6ygP8/cymTWcFymkIJViuyv+doG0wgU7aMkDUZ9Jgba/AOPP0EboUcPriumskJLfVxF27T0uA9Ev/hs84Y17Y8Fa4xQ0/8LVeq0Ztatn+4TIlTcF8AtM7uxF9yP5a/2+HZya0Ymizk7Dbg30vFqCdTGXBfKhofR/a/fjlVpXlYTGZK8w9QV1WM5e+CEDUluVQdP4y+ueH8njihx9E6uzP2sd9oKizh8J4LL/UoSENcw+xholOmkvNzFDuX1sEV7bfdiAQVxZ5ONOkrAXDWWDH9z/VOALOuHgfV3wmmSY7S3gWA0vwM7Fy8zjgQwl5lZVK8jrpmGc72Vu5Na+DnV9OpzkzHLXbkae/ToJPho3U97TKTGWSaM09n0dcbOfhjPobqYqwmI0oHNzR+UUTc8DbuPqGsXHAlNdt+RKZUIbPTYDWbkCkUyDUOWHRNuA+YjnuvofQbd/MZ9yH0TEqVBpV7AAadOC3b1YiE2cPIFXJSH1jE2scHtWvCBHAP1ZJdqiQpxIinowXDsZJT1mkqySFyQMu1wMwSO7R+MZgMOg4sfhHPgWcflbqv0I4b/+POFzdXkRBoZMEVNcxa9AhJc9eiUJ9aA9ZiBZnswk6iOHgFMeLJjWdc7hY1iMbcXWj9vUi9PYTy4xAYpUKjMVNZaubg6uMcWTiH0u0/4dV3PElTbLebi9D5LEY9Fzj9WJCQOCXbA7n7hTFo7h+tbrvYOrMAziHO7C1sOdXbL9hA7ZHtrZZbLWaaSzLp5duSMA8eV1JceJRlDyVibs4n8dLTHw0C6Ixw5YceyPpcz6OLXbBYYHi0gUvDc8hb+vJFx36+/v/064DLHiLp7i9pOHqYhlqIjFeg0bR0jPHwUTDsGncGPzwSq+U4h7+8n+XzJlFTdqzT4hRsm2ufMVQVikFkXY1ImD1UQFRSq783fldz0dt0CnZiZ35LwuzlY8KXIyyf3ZeKomya6qrZ9dUcEnzqcdRY0RlhV74dNUWbiJzWi4mPReLipjjjtk1mGU1WRyJnLiBXM4IvNmsBeOaSWky7PqEud8cFx20xW8jZ/HOr28wm4ynrHdr0M4uv0/Lz3dEsnzeBo+lfoHTyYdeXJ4vCN9ab2beuCV2TlbJcPf4pQQRPG4NCU8vap0cCYNA1sva921n18uX8+fo1/PX5o1QUZl1w/ELXo3Jwo/F4vdRhCG0kEqYAQPmOXBbft4m1/yknc7ue0gv49esc5MzhJgf2Ftghl8MrV9Tgryhnx3ND2TA7GkXGNzx3aS0AP+zQogr3ZOKj4SSMsD+vmpoymQyZTEbkVS+xYKU7pXVyXLRWnr+0hpxvZ19wv8yCQ5up2Lyw1W3NFcdY93HrqSqGplpkciUKjZy6vSsoX/sJ+uMH0ZWW8Of7JRTnm1j+2GYOfbOZX+5aRuH6coq3m6g62ET1ts2ETHkEgGOZm6jc9C3IizE1HKJk5VtsfvOaC4pd6JocPILQV5+5gpZgm8Q1TAGAqW8epCBrM8e2/kzOrxvRFWfjMyKZ1JluKO3O73eVTCEjZGwoc5c08MvdFaRGGvjt7nIWbrfH3s7KDYOb8HKycLxGzssrnQm/I/zcGz0NrW8ULkPv5plfX+aD66qZ0lfHT3v2sX/lOxe0vaBegykdeSvwzckbrVZKV76D+abXUNi1XGxKGHMjboG9yN/2Kwk3TcJiMmC1WFDYqcj89XX+emE1DmEpjJr9PQ3VJbj5hJ8oeGDQNaLStFxn9Q6JRxPYG2OtHTKFFsfoVHpd/sQFxS50TSZDM3LVmc+oCLZJJEwBALlSQUjvoYT0bpmGUVWax+a3ruPXRw4TO6Mv0QPObxqK7yBfMnOquew9mDetlsRgI09MPnnqadF2e55d6orX+Cgc/R3PsqWzCxx3Lytf+pn0rN2MjNHz4mW1jHr1XaDtfQXlCjlhSWNg98mE6egXycj3S09ZN6hXCkG9Uk653SOwFzsWPkufS+5Fbe+I2r510Yb/T5YAzu6+THl5a5vjFLqPhtJc7D3spQ5DaCORMIXTcvcJZeK89Rzc8D0HFz7NwYUGfMI8uOtbN+wUYLBA6FWnFiaQyWTEXhdL6U53pn9+lJHBdUT7mNCbZOwpsOOgwYOYO2POK1kq7ZXUyuyJe6ql6K1T2IATyxQqeyJnLuDmT65Do2g5Fauz2mNSKE6sb3SORvmP4glOYf25f9FSHv4RDBY7wm/pe2KZg38cK7/3IO4pNRYrOCec3Nf5cHT1YuQd77bpPkLPZDGZqdqzlIjJYVKHIrSRSJjCGckVcvqMmEnc0Ks4uudPjq56n6Lig9g5e+IUHYbV8/TNqGVyGb4DfPGM9yQ/o5KsKh1yBzn2kfYkxrqfd11WmUKGITCUGnVfBl/zDIq/52z+P9eYYSS9eAirSd8Sr50GmUKF+e/5nwqNEzL5ydNegaPvwHfIdWAxIVOqUahO/sLX+kaRNG8/FqMOkJ22SpEgtIf8A+ux6KrpndpL6lCENhIJUzgnuUJOZP9xRPYfh9lopCQ/g9y/viP9mU9x6d0bvyQPQvuqW41yzcs0kr26HF1pDWo3J/wHeOEdqWpTEfP1X9dQsWU3I+d/jPJ/ChCUFWSiq6/GMyiOsry9+EcPRGnX0i7pf9f9J6XmzEe2cjsNcjvRcknoWPbOnlhMJhSiW0mXIxKm0CYKOzsCIhMJiEyk8fI57F/xMQUbVnPoizU4hMfgFOpGbXYZ+opy/Eb/i6AhSdQUZnF05WKyFhUQf9NAYvqrz7kfs8nC8ZXpWA3NbH59Bs2FB1A6uaPx74XVqKf52D6Ubt4YK4qw8/BjR00lLomTGHzL2zi6+XTCMyEIF+bYnlVoA0Tz6K5IJEzhgjm4eDDoqrnAXPRN9Rxcv4j6wgOEjJtFr9TLUar+PlobOAUuf4jMDYvZ++EdNM4YQtKYUyvz/JNCKWf0S5eQt6sR5DIiBkTRXG+lptSIwk5OcK/gVsXZa6vN7Pi+gpUPJxI4ZTbJl9wvWnIJNqemvICCX19i4IMjpQ5FuAAiYQrtQq11InHCrWddJ3bodJx9wtg6fxIm/UAGTj77wB/vAAXeAf+4lugNgRGnf8u6uClI+5cPRzPc2f7KUyiUavpPuavNj0MQOkpZQSbbP/wX3qlJhMaKunhdkfgJLnSqgKgkUp9aTf6y3Wz8vqZdSvL9U3gfO+yDI7F38WrX7QrCxdr40iUYqo4yeKab1KEIF0gkTKHT+YT0ZsSz6ZTtOM7yl49SkGNCr7O0W/KUKeQ011W2y7YEob2MeHIVCntPdi8X7d+6KpEwBUl4+EUw4aXNOIWOZcvrG/np1p/58fYVrHqvmNpq8wVvt+CwCV1JCfFpN7RjtIJw8dz9wki87UOO/b4NXZNo7dUViWuYgmRUantSb3gBbngBgIqiHPb/8iq/3/05SmcPvAfHEz7YGWQyMn8vIWK0NxF9zl5xqDirGedew1tV1hEEWxEUO4idWheKjhiJiG+/Ju5C5xAJU7AZngFRjLrzI0y3vk3psYNkL3+Pbe+mY9E1YKotpXIThH4z46yF2k06M4qzzLUUBClZzBYM5cfwDesvdSjCBRAJU7A5SpW6Za7n3Z+euO3PN26gauNX1JRb8fA98321biqqsws7IcrWygsOYdA14heeKKazCGckV8jR+EdTkmcmoo94n3Q1ImEKXcLoez9nm180q598k6S7RxLZ9/Sns4p3laPxiAdgUe+2F2K/YL07b1e2bsYBcX3ubJwiUsjdcITQWJ/zamsn2A7xagldglwhZ9DMJ4i9+X12f7iJiuLTDwxqKigiMu2WTo5OEM5fwuWPUr3nEPvW6aQORWgjkTCFLiVu2FUEjLuHtc+kk5/Vusl1baUZc0MddhqtRNEJwrl5+EcSe/MHZH3xJz8/uoeCnLY3axekIU7JCl1OypVz2ahvZtPT8wj4asaJBtcGPVitVtTaUzuNXMxpwoaacja+O4uGozsw15XhmjQEU7OO0DHBxKfaYzKZWfZMJvryCuRKBchleCTFkDTNncYa0DdbyPj+CDKlB2lzf0WlvvA+iBVF2RzP3IJHaDwBkYkXvJ321qmnv7uB2CHTCO9XzIG137D55cconDiYgdOcxClaGycSptAl+fYeTtHPUFNhxdOv5TYvfwUeA/uxZ/F8Rt7xXrvsR99Uz6q5KbjFh5Fy0yi0jmDv0PpLTalUEDg8lOz/7CfipkvoN0Z74ovPxQ1AQdBjMfz22E7S35qFxiMQs74RY20xxvoKAkfdRt+xs05sb9UrV9F0bB8Wkx6LvhFzXXnr/bn5UOgTQ8Dz69rlMQrSUGudSJp8B/5xQ9jy2pUc9O5N/DDRLceWiYQpdEku3sEonD1xqVL+GgAAIABJREFU85ZxJMNAaC8lCqWcxMu8Wf34d2T1m9AuA3GOHdiA0lFL2u1n7oBiNlko3lYMgNXCaY8SlCoZcdckUJ5di9Vag9pBjn24CoXKm+zvn6Bk+084Rwyk6fghanYvJfnByTi4yVBpZKS/uAWZwgm1fxwNGb/j2j8RQ5X46HYXvmEJ+Kf9i7y139ArJQo7lTjKtFXiUyd0Sc4egcjVjix79hCN2XvJGz2atNt98PBRkHT3aHa9eyOcZeyPyaCjtqIQF88glKqWdmP66uPoKvIx6epRal1xDIqnLHsLzuFnbxe2b10zxppmoqc9TMXvX7C3wBu3CBeC04JbrReTpCImyf2U+4f3TeXg+kbqj69B5aRkyjuX4OhysrfowHtS2PTCcvRlOWC10JBbS++pD7Th2RJsXZ8xN/PXgdX88tAmhj0yGJ8gxbnvJHQ6kTCFLkml0TL5jQw2fzkXY20TFVt3s7QknPirwonsp0LxYBrw4yn3s5gtFGVvZfend6EvPoTV0Ixr7DAclOBqzCY6xISTg5WqWjnbjzrRrLfilBx+1li8Q9Xk1JbSuGsxXz5fh1JRx90vuuKf6o9Sc+6PmKOLgoGXOAOnXnsFCIpUMv2TSTTUwrGDBrK+3ULB3tVEpF55Pk+V0AVonVwZ/8QydvzyJutfWsC4F4a1asgu2AaRMIUuS6VxYMTtbwNg0Dezb/nHbH31OZidRlhc6/ZJh3etRFdbweHFz2I1NeGR3AdZnAem/DL8ZNncdVUD11/SiP0/eluXVlbx4yotby40UpurwSXM5bRxBEQocQyPwFJ2jLRBOjQqUCqAdpyOqFTKcfUA12EaIvoN57c7v0A/6xXUWqf224kgueRL76ex9Ch/vbeZiXMjxCAgGyNeDaFbUKntSb70PtTeERzbXXvK8v0f3kT+ipeJvjKO8c+nUne4CEORnFBVCWs/K+P2KxpR20H6djWL/7TnaKECHw8Ld81s4MO5FRz6+gBmw5mLwqc9EIpc2TkfJwcnOWrfENLfup6t37+IySDm83UnqbNeoSk/j+oyUQDC1ogjTKFbCRp7JzlfP8h3v1TzbebJ26e9NQylnZydy5tYdu9veAyagZe8kcevqsPdxUpWnpI7n3cjr9kNq0qNubKWK0fU8upDtaSl6BkaU0/hvgp8kk9/PVOtkcM5ZlZYzVas1pYvQZlchkwuw2KytPpbX6OnqawJmVKGg7cDdo4nj5T1dXqaSpqwWq30ujSK2gojFRk/8dvaT4ia+RJxQ68SZfm6AaVKjX1ALCW5Rjz9xGlZWyISptCt9B07i15Dr2DT53OA90/cvvbTYpz8nTDUmlB5hTHkhpfInJ/IhNSWo7PH33bBEB/DwBGByGQyzAYzP7++k+l7mxnc18BV45uYs7D0jAnzXI6tOUbhylxkfydMhYsGrY8j9dkVABjN4BLhhqy0ithwIwYj7MuzwykhgIBhARz99SiW4xXEhRtRyKHwiBKXSD/63RlGbk4QB799lMqsTQy75U2RNLsB9z7jyF/7K7Ep4eK0rA0Rr4TQ7ajtHRn179bzMJuK9OR+twz3CHv0RQfIWPEhYQFG1CpobJaxea8dJg9vZLKWw0SFSoFXghcbd7dc1IwJNaGrurBTn3kr8lDuz2bbNyXkLi9m/v01WM0Wag7XsGNhKbnLi4kNNXHnyGPsWFTCD69V8svblWz9ppRB2iPseW07d43IZdffyxa+UsmWr8u4rvdh9r6zk8jeCsY8k0LN/mVsXfTsxT15gk1InHofxnoz6Z+WYza1T2N14eKJhCn0CFPmb8Ej9Vpyfj5M6LWvEt5nGHIZrNyk5qfV9sidfNm7qJDSgpPXKZvKm/Dzavm7rlGGQt2202NWq5XcP3LRZB3mh9cq8Xa38PVSLXM/8Sb+X/1aXfN89q5a7prZwP5sO75fYU9xuRxXJyuvP1zDx09Xc8dVjRwvU/DrWg17DtnhqLXy0I0NzBhZS+H6QlzcFAy4M4mSdf9F13jqNVyha1FrnRj95EoaCpr4de5uMrfrMYnEKTlxSlboMdIe+PrE/026BvIybuahRaWAjJgZMzmSs48tH37JpS/EU7ylGHVxCZeMbAbgt3R73KLd2rS/vD/ycCnOZ9FrlXi4Wvj8Fy3zvvam7539sPdsXR6vf5yBWU+4s+GwG47+jjS8U8raT8sJ8jUzPlXHG1858vZiN5xD3ajLr+XpWeXMurSJ6WOaWDy/irBJYbj7ypEp5Sy9M4x+93xHZPL4i37OBOk4uHgw6YWNZKz9isxvXyLrRxlpjybgJKabSEYkTKFHUmociZzxcqvbXHqN4Jclz7DvExmeTSV8t6AKZwcrxRVyvl+pJeZe/zbtw600n0WvVeDmbOXTxQ68tNCLvncmonE/tfzZktVatlf4kvxIPHKFnMyvZWzZV0uQbzMlFXLeXOTKgEcHonJSUXGggpWb6pl1aRPBfmb0tXoAHJwUTH9rCId26Nn77jU03/g28aOuPbEPi9lCc0M1dmp7VKJAfZcgV8hJGHMjfUZdz1+f3c+qecsY8UgyHr6dnzSrys3k79VTW9hIc0UzV/9j2dJns/EfEkjS2O79vhIJUxAAs6GZ46s/IMTfmVtScrlrZgMqO8gtVHDTU+64DY84baI7k4RoIx89VYWrs5UfV9rz4ree9LsnEY3b6beRvk2NT7LPiQE7dlo7ahta/r9pjxq3aA9UTqoTyxqbWpY5OVgw6Vp3u+iVrMb+4TS2vfEQR354Cnv/GEx6HY1Z65HbabBarWiD43GOGU7wgEsIik0VA4VsnFwhZ9gtb7L1B09WP/EmIVNTSZ7kgKKDpzJVlZvJXFtP6ZZsjDUVOEUPQRs0AK/EcFj90on1gkY/SPa3j6Cv68fg6acvwNEdiIQp9Hi1R7aS8/X9jEvI4on/1hHgY8Zshk+XOPDal874TIgmeHDbji7/M68KB/uWEbHJfQwoTEasljPPqyssVaBN/MevcxknCh80NMlQOrTtoxoSo8T/7ZGUFVgozzdQX2Yh9Y7pePgoaG60kJdhoOTAbnZ/8D07DUY8ki8j5boX0DicvjiDID25Qs7gmU9RlHwJuz66jT/2wqj7I3Bwat+jTbPZQuYWPXlrCmk8mo1r4hR6z3qf8ITRyJX/eB/+I2H2GXUNfjEDWf/saDIDk4lN6Z5F5EXCFHq0yv0rKVp0Kx/PKSNtUMupzfU7Vcz70IXjCi9i7opC693200wO9lbeW+jIjPFNhPqbeejaWt75PpuEOxJOjMT9J6NJhuYMRwtWK8jONcnzNOxUcgIi5AREKIGTj8HeQU5siubvLzUvinLNZPyylxWPDmLQ7MX4hMS1eV9C5wmITMTnhY2s/+jfrHgqnaGzB+AdcGrSNJst5B0wkb+9hubSBnRlZZhqyrGaTai8A3HpFYydVo57qBNR/VXI5XB4l4HDK/Iw1jQRPHk2cXNmnHc1KQ//SEIve4qcZe8Sm9I930MiYQo9Vu2RrZR8fzML55fRN8ZIcbmcuW+5sibDmZjpkSTEe542uZ2PJ95x5ss/XMjKVfL23Bpum97Ir+lllGwrwS/Fr50fycUJCFPge08AO5a5suGZ4cTf+SXRAyZJHZZwFkqVmpF3fsqOn19j7ZMv4DGwP32neuLlr6C2ysLB9HpKNuVg1hvxGjSTgGF9cA+IwcU7BKWdiuLDuzi+bw0mfRNH/tjEgU8yQK5A7ROJz+BrSZpy94mmBG0RP/p68n+bz/alDQyY4ghAY70Fo85K0WEjxzaUUJ99iKs+mdDeT0mnEAlT6LFylzzD6/dX0jfGSFaekitne5N33EDKEwPwCri4j8biVVr6P9ifpR/uYeaeJob0M/Dq7BqmPngY91h31M5t/zLqSAqFnJSpTjh4DufAZ3dikn1KXPIYqcMSzkKukDNw+sNED72Kfb+9wZonvkCmtMNq1OGSMJ6Yq18jqv+klqbm/yOs7yjC+o468Xdd5XEsZhOu3sGnrNsWSpWawY/8woZnRlJ5KAaZUk5tRgZWK2h8IvBMvpz6rCzSP69i5KxTO/fYOpEwhR6psTgLZ91uJg7VYTDCjY+702gXCbJDZP2YS3Of/0loMvDs49mm07N2jnZEXh7Fo2828ucn5cSFm7h9ai3fLD5M75vaoVlnB4hNUVGRHc3Bt6/A+Yk/CYxOljok4RxcfUIYfuub6K99nobqEhxcvNp8LdrZo23X6M/GJ6Q3qU+toThrCxajnj7T38U3LOHE8poR1/Dn7AQ2OabSf4pzS1nJLkIkTKFHqj64hvGpzcjlcLxEwdSRzcDuv5dmnbJ+WZWCVUv86XtH3zbtx7O3JxlbfXl/YTP3X9/A/dfXs+z2Ysr3eV/8g+gACoWckbM82aAeRNaytwiM/krqkITzpLZ3RG0fKXUYAPiGxeMbFn/aZa4+IfR/+BcOfPsYmxvMjJzl0cnRXTiRMIUeqbk8l5iBLdMxQv3NPHZb/VnX333IjuXzz9yt5GwiL4vkvTequHR0M2EBZl59qIYZT+ZgbDZe0PY6Q/RwF9Kf/gOLydR6ZKQgtIOwhJG4+HzLmocTqZw4Dg+frlGMQXwShB7J3FzHjgMqaurP73RQcbninKXxFCoFHyxyRKkAo1mOTNEyYEjjpsF3VBj3zdczdnBLPVo/x2by6pR8tsQBe42Viho5IaqT21eoFfy1S02zXsauTDvkanmr/RRXyHnnW0dMZtpcsu98eAcosHP14OjeNUT2H9fu2xcEd59QZEoV1i5U8U8kTKFHCpn8MFu3R7DxaDOHV3yAY6gPrn5qHF1PPypWJpMRMc3rrNuMuyGORTnVAIRf7YziHwkwcHggRQoZ3+QZAFAmKukz2YElR1vqvrqO1WLvdbJcXsCwAA5vkZOVZ0bhrCBwkO+JZQ7+DmhHxPJNXkvyjbu+bSX7zpdnUhgFW34UCVPoEHWVx7Eadbh5X9hIdCmIhCn0SPbeEYRMfpiCzC3UL3mNSbMvft6YW7TbGevNyhQyAocHnnK7R+zpr9+onFSEjA05/bZkMvzbWEjhQoQPcmHji0vY+Ws8fSfcfkHTDAThTI7s+B3HqNgu1b6s60QqCB1gx1tX4xgZLXUYNikgXEHsjcMp+OMVDm3+SepwhG6mbOcvePezzcFvZyISptAjleTuY+X8y7Gamkm4KkzqcGxW/DAN2qAgjE11UocidDMNOVtxD1ZJHUabiIQp9DjpH/ybjfPScPBt4pLXRxEUJa5MnMn+DTrqc7IJ7DNc6lCEbsY+qA+1xbY7Uvx0RMIUehx9TTGhlw5i+HXuaLRdZ8BBZ8s9aCTz83UkP7QEr6BeUocjdCNmo5GmvF2E9utaRdrFT2uhRykryKT+0HpSbhgrdSg2z9Dc0i4lZ/m75G78HrWLDxpXH3yjB+IX3rYCDoLwT0f3/InKwxd3r64x//L/iSNMoUfJ+PlVfEcO6DITpaUU01/F4EfT8OjVjMohB0P1Woo3fsCuz+6mrqqEfau/oraiSOowhS4ob90XePXvemMHxBGm0GMUZm2nevdSUueMkDqULiMoUklQpOOJv4/sd2P3h9tZ/dgg7Fydyf78LlSeobBSuhiFrqUkdx91B/5k6E1dr7i/SJhCj7Hrw1vwGpyAb6g4sXKhNFo5xurjRF43kQGTHDEZenP0gAHYf2KdbT+9TtzIa3F085EuUMFm6Rqqkavt0bS9zazkxDeH0GPY+0Xj5KftUhOlbY1vqJxxb1zOgEktR51KlYzoxNYFDWqyFrP83miyt/8uRYiCjQuOG4baM5ztSxukDqXNxBGm0CMs6i2DL6SOomeYODuErJ2+7H3/Rjxe2IyHv2100BBsg1whp99Nb7F53mhih4/D1aPr/IAVCVMQhHYX01/NAQcn9i5egKmhlOaSHJwiB+GdMI7wxHFonbte82Ch/fhH9MNn5G2se+MPJj/VG6Wqa0zvEglTEIR2ZzJZQCHH3JSJV7wbbhPiKMks5tiqV8j8+FbsA2JxCO0PVgtWkx6LyYBTUB8Spz6Aqite3BLabMh1L7Di+V2kf1bKqNu9u8SlEpEwhW5rxgHrKbfpm+r54+EkwqbEkTTWQYKoegalUs7lrw1sdVtEvApwQ6+L4uh+I1V5+cgVMuR2cuQKGWW7f2HZ6o+IuHIefUZdj7wLfIEKF06uVDL8gYX8+eQw/vpSybAb3G0+aYqEKfQY5QWH2PTqdJyjA0WylJBaIyd2gBoG/E/3k0lOZO0M4uB3L5C/4h0G3PkffMMSpAlS6BQOLh6MemYN6c+NYd1/LYy4ydOmk6btRiYI7ej4kT2sf2oY/kPDGH2HmO5gq2L6q5k6Pwm/wX5sfC6Nw7tXSR2S0MGc3X0Z/cwaarPL2fBNjdThnJVImEKPcGDJC/inJTPoMmeb/gUrgEIpJ2WqE31uHc7et2ay708xvLm7c3T1ZuTjf1C2eT/7NzRLHc4ZiW8OoduzmMzU7PyN2JHOUocitEFciprBc0ZzeOEc1rx9EwWZm7GYLVKHJXQQZw9/Eu/9lszP11NaYJY6nNMSCVPo9uRKBS79JpG5rl7qUIQ2CoqyY+y8EaicitnxxhUsvk7L4lkerHhxKrrGWqnDE9pZWPwI/Cfcz85vjkkdymmJhCn0CAkznuH4qi3odeIIpatx9VIw4np3LntnGJd9Mo1xL49FJitny5ePSh2a0AESJt5J/aEMGmpt77MqEqbQI5Tl7kOucYCuMT9aOA2FQo5aI8fVQ07sJYHUZIgBQd2Rg4sHTrEj2b/G9krniYQp9AgWkwFzUwObvq0SR5ndQGCkAlN9BVXFuVKHInQAjXcE1dnlUodxCpEwhR6h37ibGf3KbppK7fjj2f3UVIqk2ZUplXKce8WRt2uF1KEIHaDfZQ/TmFvA/g06qUNpRSRMocdw9wll3KO/4tZ7Cn8++ZdNXiMRzp9HnAeV4rRst+TiGUDsrHc4sjRL6lBaEQlT6FHkSgVDZ72MY1Qqe1fUSR2OcBG0bnYY623vtJ3QPqIGTkZfWkhVue1MMREJU+iRoifdx/HVO8Wp2S6sPKse56hUqcMQOohSpcG5Txo5WxqlDuUEkTCFHimkdyrew24ifcFezCaRNLui2pxifHqPkDoMoQP5Jl9G+e4SqcM4QSRMoccacv2LgD0ZG/VShyK0ka7Zgq7oKMG9h0oditCBmquLMett5/MpEqbQY8kVcsKmPsLRpZmYRcm1LmX7L3U4RA5Cbe8odShCBypc+jJJN8VKHcYJImEKPVpc6pVYjFaydxqlDkU4Twe36Ches5NB//5M6lCEDmYf2IfyPHGEKQg2Qa5UEDJ5Njm/HZY6FOE8ZG7Tk/HpOhLv+RZ3vzCpwxE6WPCYOzi22nY+myJhCj1efNoNGCqrObzXIHUowlkc3KJj/8fpJN73PeGJaVKHI3SCuKFXYqqrJ2unbRxlioQp9HhKlZrgqXM4+EO2uJZpow5s0pHx6TqSHvhRJMseRK5U0uvGt9n70XqbKGkpEqYgAP3G346xvpms7eJapq0pzjeR8Vk6SQ8tIazvKKnDETpZ7JBpyO001FSIhCkINkFhZ0fklc9x6Id9Yl6mjdn1TT4BE+4nLF7MuezRLNK3GhIJUxD+FjdsBnI7Z/aubZY6FOFvWTv1NBcV0/+y2VKHIkjEoGvEWFOGh6/06Ur6CATBRsgVcmKvfokjS7bT3CiOMqVWW2Vh/3+3EjVzPiqNVupwBIls+vwRXOKTUKrEEaYg2JTI/uNwihnGpq9FUW8p6XUW0l/NwGPADPqMukbqcASJ7Fz6PtW7fyH1X8FShwKIhCkIpxh8y9tU7tzPkf1imokUivNN/HLPKlQeMQy9cYHU4QgSObx7FXmLn2LI7FRc3BRShwOIhCkIp3B08yH8yhfY//VBqUPpkQ6trsZz0DWMmf0DcqVtfFEKnctitnDg8/uIvmYY/mG28x4QCVMQTsNqMSPXqKUOo0eq3n+U0KEzkSvE11NPdezgX1iNjcQPt63PoHhHCsJpRA6aSnNBHpUlttO8tseQyVHa2dYXpdC5ctd9jXdKLAob+9FkW9EIgo1w9vDHe/iNbPowm8Z6MWK2M8lVKgx622kaLHS+upxNBCU6SR3GKUTCFIQzGHLDfDQ+SSx/fDOFR0xSh9NjyFUqTLomqcMQJNJQXYqxqgCvYNu5dvn/RMIUhDNQqjSMvuczgifPZdNLq1n7nwp0TVapw+r25HZ2mPQiYfZUB9Z+hUufBBwcbS892V5EgmBjkibfQdqC3eirnfjtoXXs+0snirR3IIvBgEIlrmH2VOXbfiRoiLfUYZyWSJiCcB5cvYMZO2cx8bd/Rs6P+/jjpSPiNG0HMJstNBfl4RuRJHUoggQsZgu644cIjLaTOpTTEglTENogasBEpryRgVvsNDbNX8sfr+RRlCtG0raXmgqQKexw9vCXOhRBAk31FWA24egqfRm801FKHYBgGxbpbfMNarOm/f2PCmAziJzZPjyBDwHxfuyRNFoXrGYjJhMobG/MjzjCFARBEGyDUqVG6epLVYltjhEQCVMQBEGwGUonT5pqbTNhilOyPdgMtZgi0Z4MukZ2/PASx/98n+BJqQyebnsTr23Z3rXNHPk1k/Hzt6JxcJE6HEEijqH9yV6RQWjvYFHpRxC6K5XGgSHXP8/gJ1dTsGKr1OF0ORU5dXgPvV4kyx4u9ZY30Ffo2LSoTupQTiESpiC0M7/QBCzN9aIJdRvpKhpx8gmXOgxBYiqNA8Pm/EzJX3vJ2qmXOpxWRMIUhHYmVypQeYVSVigSZlvoy8tx9YuUOgzBBrh6BxM48UFy08ukDqUVkTAFoQNovEKpKRYNqM/Xrj8bsRjNeIf0ljoUwUb0GnE1tRn7aKixnR+eImEKQgdQ2LtgbLadD7otKzpi4vDCv0i+fxEqjVbqcAQb4ezhj0vfCexZUSt1KCeIhCkIHUChccSoE9UMzqWx3syWt7YSPuNFgnqlSB2OYGN6TX6A4rW7MRlsY0S/SJiC0AHkKgfMepEwz2X3slq0IQNInHiH1KEINigodhBqn2i2/14vdSiASJiC0CHkdvaYDCJhnktDYQNeCeOkDkOwYYk3v03B0k1Ulkn/eRIJUxA6gKGmEK2baFF1No31ZhqPFeAeLAb6CGfmH9EPn5G3s+Gtg+gkHhcgEqYgdABd2RFc/ETCPJNtv9Wz7MHVOIYPJiAqWepwBBs35LoXsHMJZ+sP1ZLGIRKmIHQAfXk+HgGi8uTplBaYyftlE4MeW8GY2QtRimbRwjnIlQpC026lJqtE0jjEJ1oQ2plB14i5oRJXr57dokqvs3As00jxgQZqc8rRlRQBYDHq8Rp6AwFR/SWOUOhKHD0CMTVKO/hHJExBaGflhdnYefihVPbMEzi6Zgubv62kdN0mNH5ROEcPI3TCrfhEJCJTKMFixdUrWOowhS5GobTDajZJGoNImILQzqoLD6H28pY6jE7VWG8md5+B0ow6qvZk4hg5mLFvHcLVK0jq0IRuQqm0w2oSCVMQupW641lofXtGay+z2cKu5U0cXfIX9oF9cI0bzYAHnycodrDUoQndjKt3CObGGpoaLWgdpDl7IxKmILSzpuOHcIvq/iXeKkvNbP7oCMZ6PYOfWE1AVJLUIQndmFKlQRMQR0GWiZgklSQx9MyLLILQgZpLc3AP0EgdRoeqKDaz5qm1OIWnMenlHSJZCp3CKXwg5dkNku1fJExBaEcWswV9yWG8ghVSh9JhGuvNrF+wg4Dx9zF01stiWojQadyjBlF7uEqy/YuEKQjtqKG2pX+fo4vEgXSg9LdycYwcxoDpc6UORehhwpMn0nA4k8YGaSr+iIQpCO2ovrIIO1dPFIru+dGqrTbTkJPJiH9/hLybPkbBdjm6euEQmUL2Np0k+xfveEFoRw1Vx7Fz7r6HlyVHTWgCYsVpWEEyrr1GUpldJ8m+RcIUhHbUVF2CnYuD1GF0GGdPBfqKfCwm6TtHCD2Tb2wq9YcLJdm3SJiC0I50NaWonLvv0ZdXgBxLUw26ZtvoTyj0PEGxgzE1NXN4j7HT9y0SpiC0o/q8XWi9u2/CLC+yoHD2QuvkKnUoQg+lVGkIm/YEGd9loNd17uAfkTAFoZ2UFxyiLnMt8aO6b5Wfkhw92sA+Uoch9HB9x92CyiOaVa8eQa/vvKQpEqYgtJN9S17CZ9hAHJy658fKaLCQ/2cOfoNmSB2K0MPJlUrSZn+PTOnN2nc773pm9/xkC0Ins5gtVG1ZRL9L3KQOpUNUlplZ9vR+7FxC6DPqOqnDEQSUKjWjZ/9Aw9FjZO/Wd8o+RcIUhHYgV8iRq7vn6Fi93kL6i9twj7+UcXOXorD7v/buNTqq8l4D+DN7JslkkslkJmTIFXIFEgghIYjc5aq2yEXUVQurahVXVxc5VYq061i7vJ16RD3arqrVg3osCudwqVpFihQiBoEACZBAEnKDEJIh5DaTZK6Zved8YJ11yqomM5k9MyE8v4+szP/9z4fh2fvd+33fsFC3RAQACFdrkHHf8zi/4zxEd+CnZhmYRDJRxRjR2xmaHUgCqWyXGWrjRMxa+yIE1ejd8o9uTnmLfgKFQoOzXwd+MwMGJpFMhIgo2PtH1/pE00U3TCUnMGP9m9zZh0YkQSlg8rpX0Lj7OOzWwF6w8hdAJANLZyucV+swPic0xw4FgihKOPF+LVLv+TXiEjND3Q7R98ooWAzthHk4vjOwG7MzMIlk0HDsE8Tk5kGtUYS6FdmUfdoHj0uJGas2hroVoiHNePg1dHxbgUs1gdvQgIFJJANVRBQ8kifUbchCFCWUbu9B26FzmPnEdggqnjNPI1+scRyy1r6K8j+dhMMemKlZBiaRDAypk+Bo7wh1G35z2CUcfMuEjvLLWPjiERhTc0LdEpHX8pc+jMjkqSj/3BKQ+gxMIhnEj8uFq6MF7iC82h5Y6ijbAAAMdUlEQVQorU0i9v7qGKQBI5a+UIrY+NRQt0Tks5j0Ijj7AjMty8AkkkF4RDSgUMLad/NOy1b8uQ5jFzyOpZt3QRNjCHU7RD7rN19Dx8m/QKUOzGMEBiaRDK7UlSFMHw+d/uZcp3ixegCO9g7MWP3LULdC5BdXeyNyFwXmTFoGJpEMWsr3IXbyzbn0otMkouLdcoxb/hRU4epQt0M0bNGxRsTNWYvzB8wBqc/AJJJBT9V+JOUH5qrWF9dafds4wdTsRslzXyNx/mMoXF4coK6IgkdjzIDb5h7y70RRQmOVC5WHbV7X5vviRH7qbK2Do70eGfkZIe2j3yLh4KZdSLp7Keau1UM5xM48bpcHx9+sRMoPnsJt924KUpdEgaUdm47m/z6E6vx7kTvrxrNprf0STI1uXK2xov1YDQAV1AlZXtdmYBL56fwXf4BxVhEi1KGdsLFbr79w1PblV9hTGgeP6IYudzImr0xCYroS/RYgNk6AKEqoPeFC/ae1CNONQ9FKbkxAo8ek2Wvgsm1F1X9uQu3uMZi+vgCx8Qoc32ZCT0U51InZiEqbjvzHN2P8lAUQlAJ23ufdhiMMTCI/OKwWdBzdjkUvLA11K4hPUiKveBWqPziMqMxZKHxoCxqP7MKxV96E5OiDQhkGpUYLQR0JSEqkr3oaUxb8mHvE0qgiKAVMu/NRTJp7Hy4c/QTHXtoAiG7E3/EofvjOX6HRxg67NgOTyA+VB96HNmsC4pNGxtuxU2arkZa7FEfeu4Ljr92PFW+cg3vNZrgcdqijdDA1nYG1pw1ZhXfz5BEa1dRROuQvfRh5C9cBCkGWC0NeWhINkyRKMB16F5nLkkPdyg2iYwXkrU6G6OgHAKjC1dDE6CEoBSRnF2LCbcsZlnTLEFQq2WZRGJhEw9RSewySsw9Z+SPvQGV7nwjJaUXD6QOhboVo1Lilp2TXbb45183RCKKLAH57MtRdfLf4GGDHz4AdoW6EaHTgHSYREZEXGJhEREReuOWmZD/a0hjqFugm13D6AGp3bMTK3+WFupXvVVVqR+2OU1i2pRzRscZQt0M0KtxygUnkL01MPNzWvlC38Z3cbgn1FQOo+XMpCjfuYVgSyYiBSeSj2PhxGDB3wumQQr67DwA4HRLamyU0HOlB18lKqLRGTFj3OtLzFoS6NaJRhYFJ5CNNjAGRKbloOuNCzu3BOd1j36vNcHZZEB6jhVIdhoFeG1wWCwYsnZAcNoQbkmGYvhIzf70FyVkFQemJ6FbDwCQahpjsuTCdrwxaYCoEAerEQiQU3AmX1QKNPgHRccnQGcchOsbIjQiIgoCBSeSj3m4TOo5+hOkb5gdtTI1RA7dLh7yFa4M2JhHdKPQPYIhuMmf+8jIMBZOROTU8aGOaa69i7JTFQRuPiP4ZA5PIR30XTyG5QB/cQRUKuF324I5JRDfglCyRDyRRgv1KNcZmyr/h+oVyJ2r31EGbmYCJiwwYk6RAV5sHbRfscHaYYOsxyT4mEXmPgUnkA3PnZSgUCugM3h0464uqD8qQtOwJOLqv4NuXdkCymaHSJ0E9NgtTf74NWYXLZB+TiLzHwCTyQXtDBSJT0qAMwKHLYToD4sblIeveTZAefg0ejwRl2Mg7CYXoVsVnmEQ+6GqsQFSqHhUHrLhwygnRLclWO9Kox8XD29DZWo+jHz2NqpJtstUmIv8xMIl80Nd0DHEZWtR/fAjnPz6PvS/Uyla78IEECIoWlGzKg6X+b2jY9iT6e9plq09E/mFgEnlJEiVYmyogqBQIM6Qg7Z5fQaWOlK3+mEQlCu9PQMzkabA3NyIq8zZIkihbfSLyD59hEnmp+2ojJEc/Lh5sRMzEeeisLoF+okG2+tVlTpzbehiJi3+OucV7EaWLk602EfmPd5hEXjKMzcCU4u3Qpi1G1uJH0XuhFCl5UbLUtnSJqHr7K+Rv2I7Z655nWBKNQLzDJPKSrb8LvaYGZC38CdQaHUSbGUnp8lxztl92Q52Ug6zpXDpCNFIxMIm8VF/2V5gO/hFtX/0eHvcAtNmToFTJE5jmVifUCVmy1CKiwGBgEnnJ3HgKCQumYtYaHVrqRERq/du8wOmQ0Nl6fVnK1eMtSLqjWI42iShAGJhEXrJeKkfqmnFQKgWk5fh3Z3mhwoWq98ugUIQBShXC9UnIW7ROpk6JKBAYmERecLucsF+pRsqEDL9rVfzdisbdJzDpoT8gd979MnRHRMHAwCQagstpR8mrD0A7cQqitP4f1Hz5QB1yfvoWcmavlqE7IgoWBibRIBxWCw6+tAIqtQNLNqX7Xa/fIsF+uQ6C4H/wElFwcR0m0SDqT3wBa/0x2K92oOTNVhzdbYbD5hl2vWidgIKNq3Du3cdw5IPNcDl5xiXRzULh8Xj/41coFJ4Hdg//Pwuim5HLaYepoRyn314P0dqBO55dDGOyf3eI5g4RR95pQn/deagTJyBxwSMoWvELmTomIl/svE8Bj8cz5GvvnJIlGoK9rxuVH25EuF6Dec8sQWyc/xMzsfFKLP9NNhz2TJS81Yz+1gsydEpEgcQpWaIhVPzPc9AkReIHT2fLEpb/SHR70F/fhOwlj8lal4jkxztMokFYLV3oKtuJJS/dKduuPgBg6RFx5nMzrh4+AUPRaiRnF8pWm4gCg4FJNIi6o3ugzZ6EuLHyvdVqanbjm+cPQF9wD+a/eBTG1BzZahNR4DAwiQYhqMKhUMq7BKRiezOS73oSsx78rax1iSiw+AyTaBCRujEY6LfJVq/+tAu2ljZMX71JtppEFBy8wyQaRKR2DNzWfr/rWHpEVO3vRduhU8h68BWEqzUydEdEwcTAJBqEJjYebnMXBlwSwsKHNyHTeM6FU68fROzUu3DbU58jNed2mbskomBgYBINQm9MhyZjBko/7MSi9UafP+92Syh/owSTH38POXPuDUCHRBQsfIZJNAhBKWBO8X+hveQwLD2iz59XKAAhMgrRhiTZe+vtasOJT17HparDcDmsKH3vl/jbc8vQUnNc9rGIiIFJNKSL5fuhyciBTu/727JKpQBD/gRcKvtEtn56u9pw+N1ifPVkLsw1u3DmnUfw6UN6WNuOIDZHgxMv/xBlO38n23hEdB0Dk2gIgioMngE3RFEa1ucz5xpw7eutqDv5pd+9eCQP9m/IhLnqM4xfMQeOjm5oM1IxdcNyGCaOQW9TJzyShLisIr/HIqIbcfN1oiFIbhH7X14JpdCFKatSkZo99KN/t8uDym9suLi3GpLDBndPGwBg5rPfYPyUeX71U3fyS1T+/kcIMyQjc/UzMJV/BtvlSkRnzEBs9iyMz1+EMSkT/RqD6Fbi7ebrDEwiL/RbOlH6zi9gPbcPkePTMfX+DKTlhv3T31n7RJzcY0b7Nyeu/4NHgjJaD7e5HcYlxbj9wWcQodH63Y/DakF4RDQEFc/VJPIXTyshklGUNg626gOQbD2w1vTg+L/XoCa/CPMeTUJ07P8/2Ti7vxfdtd2AJCJtzXxMWxqNCLWAy3VuHN+yFZ3zfyTLvrHqKJ3fNYjIN3yGSeQFhaDA6q0mJC3fDCgEqJPTYD5dhi+KP8O+N9pQe8qBAZeE3hYrFEI4jHNnYuY9MYhQX/+JjZugQtqqOSh/+xG4Xc4QfxsiGg5OyRL5qKrkY9R9+C8wziuEo8cNy9lySC474PEAHgkqQwqWvTAbOsON16OiKGHfv9VDkqIwcfVvkDFtCQQlr1mJQo3PMIkCqK3xDM5u2wR7WzXSVxQhcWIkLpX3wdbhwKwfGxGt++4gdDoknP27FS37K6HUxCFp0XokZM/A2LQpUIWrg/wtiAhgYBIFRdPpg6jZ/Syc7fVIXFiAaXfpvjcs/1FXu4imchs6z3fD1mrCQFcbwo2ZWPCvXyJ27PggdE5E/4cv/RAFQUbBYmQULEZLbRlqv/gP7H1iHyISUqEeo0d0SjRyFsYgznjjm6xd7RJKnv8W0oALyUunIyYtBj21enigR7Q+IUTfhIiGwjtMIhn197Tj2qVzsJga0F1/DD3ln0KTmgG1MQ7hWgF9zd2wNl1A8t1PIn7SHDSXbocyIgph2jEoWr1JliUnROQbTskSjQC23h40V5ag/1ojHJZOaJMmIGfeAwxGohGEU7JEI4AmRo+cuTylhGg08PkOM4C9EBERhUKzx+NJG+qPfApMIiKiWxVXTRMREXmBgUlEROQFBiYREZEXGJhEREReYGASERF5gYFJRETkBQYmERGRFxiYREREXmBgEhEReeF/AZZRNnivrQMqAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 468x360 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fontsize = 24\n",
    "\n",
    "fig = mp.figure(figsize=(6.5,5))\n",
    "\n",
    "map_proj = cartopy.crs.PlateCarree()\n",
    "#map_proj = cartopy.crs.Mercator(central_longitude=-64.0, min_latitude=latmin_d01, max_latitude=latmax_d01, globe=None, latitude_true_scale=5.25)\n",
    "#data_proj = cartopy.crs.Mercator(central_longitude=0.0, min_latitude=-80.0, max_latitude=84.0, globe=None, latitude_true_scale=0.0)\n",
    "\n",
    "ax = fig.add_subplot(111, projection=map_proj)\n",
    "ax.set_extent([lonmin_d01-20,lonmax_d01+20,latmin_d01-20,latmax_d01+20], crs=map_proj)\n",
    "\n",
    "\n",
    "ax.add_feature(cartopy.feature.LAND, color='forestgreen', alpha=0.75)\n",
    "ax.add_feature(cartopy.feature.OCEAN, color='dodgerblue', alpha=0.75)\n",
    "\n",
    "\n",
    "ax.add_patch(matplotlib.patches.Rectangle(xy=[lonmin_d01, latmin_d01], width=d01_width, height=d01_height,\n",
    "                                facecolor='none',\n",
    "                                edgecolor='firebrick',\\\n",
    "                                linewidth=3,\\\n",
    "                                transform=map_proj, zorder=2) )\n",
    "\n",
    "\n",
    "ax.add_patch(matplotlib.patches.Rectangle(xy=[lonmin_d02, latmin_d02], width=d02_width, height=d02_height,\n",
    "                                facecolor='none',\n",
    "                                edgecolor='orange',\\\n",
    "                                linewidth=3,\\\n",
    "                                transform=map_proj, zorder=2) )\n",
    "\n",
    "\n",
    "ax.add_patch(matplotlib.patches.Rectangle(xy=[lonmin_d03, latmin_d03], width=d03_width, height=d03_height,\n",
    "                                facecolor='none',\n",
    "                                edgecolor='gold',\\\n",
    "                                linewidth=3,\\\n",
    "                                transform=map_proj, zorder=2) )\n",
    "\n",
    "text1 = ax.text(s='50 km', x=lonmin_d01+2.5, y=latmax_d01-2.5, ha='left', va='top', color='firebrick', fontsize=fontsize, transform=map_proj)\n",
    "text2 = ax.text(s='10 km', x=lonmin_d02+2.5, y=latmax_d02-2.5, ha='left', va='top', color='orange', fontsize=fontsize, transform=map_proj)\n",
    "text3 = ax.text(s='2 km', x=lonmin_d03+2.5, y=latmax_d03-2.5, ha='left', va='top', color='gold', fontsize=fontsize, transform=map_proj)\n",
    "\n",
    "text1.set_path_effects([matplotlib.patheffects.Stroke(linewidth=2, foreground='0.1'), matplotlib.patheffects.Normal()])\n",
    "text2.set_path_effects([matplotlib.patheffects.Stroke(linewidth=2, foreground='0.1'), matplotlib.patheffects.Normal()])\n",
    "text3.set_path_effects([matplotlib.patheffects.Stroke(linewidth=2, foreground='0.1'), matplotlib.patheffects.Normal()])\n",
    "\n",
    "ax.coastlines()\n",
    "\n",
    "fig.tight_layout()\n",
    "fig.savefig('manuscript_fig_S1_nested_domains.png', dpi=600, bbox_inches='tight')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.6.5"
  },
  "varInspector": {
   "cols": {
    "lenName": 16,
    "lenType": 16,
    "lenVar": 40
   },
   "kernels_config": {
    "python": {
     "delete_cmd_postfix": "",
     "delete_cmd_prefix": "del ",
     "library": "var_list.py",
     "varRefreshCmd": "print(var_dic_list())"
    },
    "r": {
     "delete_cmd_postfix": ") ",
     "delete_cmd_prefix": "rm(",
     "library": "var_list.r",
     "varRefreshCmd": "cat(var_dic_list()) "
    }
   },
   "types_to_exclude": [
    "module",
    "function",
    "builtin_function_or_method",
    "instance",
    "_Feature"
   ],
   "window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}