Skip to content

Instantly share code, notes, and snippets.

@lkluft
Created May 16, 2019 12:53
Show Gist options
  • Save lkluft/74aaa9eb62eb29d19abe1236ce278d4a to your computer and use it in GitHub Desktop.
Save lkluft/74aaa9eb62eb29d19abe1236ce278d4a to your computer and use it in GitHub Desktop.
RRTMG with McICA
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"metadata": {
"ExecuteTime": {
"start_time": "2019-05-16T12:52:39.807763Z",
"end_time": "2019-05-16T12:52:42.046320Z"
},
"trusted": true
},
"cell_type": "code",
"source": "%matplotlib inline\nimport matplotlib.pyplot as plt\nimport typhon as ty\n\nimport konrad\n\n\nty.plots.styles.use()\n\n\ndef bulk_profile(plev, x=1, pbottom=900e2, ptop=750e2):\n tau = np.zeros_like(plev)\n tau[np.logical_and(plev < pbottom, plev > ptop)] = x\n \n return tau\n\n\nplev, _ = konrad.utils.get_pressure_grids()\natmosphere = konrad.atmosphere.Atmosphere(plev)\nsurface = konrad.surface.SurfaceHeatCapacity.from_atmosphere(atmosphere)\n\ncloud = konrad.cloud.DirectInputCloud(\n atmosphere['plev'].size,\n cloud_fraction=bulk_profile(plev, x=0.15, ptop=750e2, pbottom=900e2),\n lw_optical_thickness=0.02,\n sw_optical_thickness=0.02,\n)\n\nrrtmg = konrad.radiation.RRTMG(\n mcica=True,\n cloud_optical_properties='direct_input',\n)\n\nfig, ax = plt.subplots()\nnp.random.seed(1234)\nfor i in range(20):\n rrtmg.update_heatingrates(atmosphere, surface, cloud)\n lw, = ty.plots.profile_p(plev, rrtmg['lw_htngrt'][-1], c='tab:red', alpha=0.1)\n sw, = ty.plots.profile_p(plev, rrtmg['sw_htngrt'][-1], c='tab:blue', alpha=0.1)\nax.legend([lw, sw], ['Longwave', 'Shortwave'])\nax.set_ylim(top=100e2)\nax.set_xlim(-3, 1)\nax.set_xlabel('Heating [K/day]')",
"execution_count": 1,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 1,
"data": {
"text/plain": "Text(0.5,0,'Heating [K/day]')"
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/plain": "<Figure size 720x444.984 with 1 Axes>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAoEAAAGQCAYAAADRMbobAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xl8XHW9//HX55zZs8wkmba0lC600LCjFpDtXi4UAa/aq4AoCC4oCCoqAmpbQEAEUeSqbIIbIiCoPy9yQbysggXRyiIUutCF7stkncns53x/f5xJm4a0JJDkJJPP8/HIo8k533P6nskk88n3fL/fI8YYlFJKKaXU2GL5HUAppZRSSg0/LQKVUkoppcYgLQKVUkoppcYgLQKVUkoppcYgLQKVUkoppcYgLQKVUkoppcYgX4tAEZkpIj8RkZdExBGRJ/toY4nIpSKyTkSyIvKEiOzXR7v9ReTJSpt1IrJARGRYHohSSiml1Cjjd0/gfsD7gWWVj77MB+YB1wJzgTLwqIg0dDcQkSbgscq+ucA1leO+PmTJlVJKKaVGMfFzsWgRsYwxbuXz3wFJY8wxPfZHgc3A940xV1a21QNrgGuNMddWts0HLgKmGmM6K9suAy4EJhhjCsP3qJRSSimlRj5fewK7C8BdOAKoA+7rcUwn8DBwUo92JwJ/7i4AK+4D4pVzKKWUUkqpHvy+HPxWZgEOsLzX9iVAc692r/Vqs7xybDNKKaWUUmoHI70IbAAyxhin1/b2yr6e7Tp6Nqgck+nVTimllFJKAQG/A/RDX4MWpY/t/W3n7RA5BzgHYN99933P4sWL30lGpZRSwySdTtPy9W+Qf/VVjOMQbm6m8XvXkUgk/I42pFKpFO+/9XkAHvr8u0kmkz4nUj4alNVPRnpPYBtQJyJ2r+1xvN7Anu12+OmvHFPTq902xpjbjDGzjTGzo9HoIEZWSik1lOrq6ph2803UHn442DaFZctoOf8LtL/8st/RhtSKyvWucMDWAlANipFeBC4FbGBmr+3NeOMCe7brPfZvr8qxS1BKKVV1Jl/zHcZf9DXsRIJSSwtbvjmPlv+53+9YQ6Y1nSUWtJkUD/sdRVWJkV4ELgTSwKndG0SkDjgB+FOPdg8DJ1T2dTsFb5zgM8OQUymllA+aPvABxv/ohwQbG3A6O9nygx+w8ZZb/Y41JDZ0FrBtqA31vjim1Nvj9x1DYiJyioicAuwOjOv+WkRixpg88D1gnoicLyJzgN8DeeAnPU51K1ACfi8ic0TkPCoLTOsagUopVd0SU6fSePPNhKdPB2Nov/deNt50k9+xBt2KVCfpvIPb91B3pQbM74kh44Hf9trW/fV0YDVwNd4AyAV4M32fA+YYY9q6DzDGtFQKxB8DD+CNEbwG+O5QhldKKTUyNDQ00HDHL1lzydfJvvQSmf97hHVtbUxesMDvaINmc1sRgEhwpF/EU6OFr3cMGSlmz55tFi1a5HcMpcY8x3HIZrOICK7rUi6XcV0X27YREejowCmVMIAdCBAIhTDGUMpmMeUydHS++aSWhR0KErADuK5DuTONKZUQ28aORMCycLJZTGcHEokQqKnBAE5nGlMsQi4LjvOmc1rRKHY8jpPLUe7sxGSzIIIViyF2ALdYwGnvANsikEhAIIjT0Y7pykLl966by207pQgQDmMnEljBIKUtW3C7usAYJBrFCodxCwVMPg/FInZDIxIO4XZ24nSmwfXW3jfGQKm4/bzBEHYijhWNUs5kKKdaENdFwmEkGsWUSl4Ox8Gqr8eORXE70zidnZhej9tUziuBAHZtHVZ9HeVsDretFVwDxmDV1GCM8bI7DlY0gl0fx81mcTs6vO/TjmfFlMrec1dbi5VIYPJ5nLa2bc+7FYtiLAuTy2OKBSQUxm5IYBwXJ5Xyvk89X0fZLFT+n8h++zH9lpv78eob+U679WlWpXJMikeYMa6GaMBij8Yoe06o510TQjpZZGwZlNnBfvcEqj60tbVh2zaWZeG6Lo7j4DgOlmVhZbOI6+KUy5hyGSsUwg4GERFK+Tzu1q19L4pjCZYxBKMxDFDKZDCFAohgRyKIbeO2t+EWi2AHCNbWIIEA5c5O3B6/jHckSCSCCYexArb3RpjJeL/Mo1EkGKTc0em9adkWdjyOFQ7jdHbiZjLguDu8CW47aziMHa/HsmzKmbT3Ruo6SDiMFYvhlkq47R3eG0w8jl1fh0lncHq8wZjijqMAxA5g19Ui9fWUO9O4bW1eIRAOY9fU4Lqul6n7OW1sxM3ncdpavTeoCtPjzRXLwq6pxYrXY0olnNY23FwOsSzvjVAEUylOxBLsxkZwXJy2tje9aZlSaXvWSAQ7kQBjcDo6cLNZxBIk4j2nbj6HKRYRsbxzAuXWVu/72fuPOmPAuN6bZn09EgxsL1gcB4lEIBiCUhG3UEBcgxWvx4rGcFpbcfP5Pr7vYBzHe5y1NViRKE5XBjeb8wqWQAAJhXDLZUyhgLguVjSKVV/vfe9zuW0Fy5u/+YIdjSI1NTi5nJfTdRHL8l5rrus9d66LFQxixeNewVBptzMSCmHX1no/N5mM11YEKxIBEe9xui5YFoF4HOM4OJUiZqfnDASwuwvG9PYizAqFIBDAFIvbXo92fT1YFm46/abCage2jR2LIYEATiaz7XUhwSASDGJKJe+cxmDFYliRCE46vcPrZ6fPaTjsFbvFovd9sm0kHMYtl72CyXWxQiGs2lqcri7v9bQL234euwvTHt8n13GgVPK2BQLY9fW4+bz3vd9Fx4OEQlgtLd5rp6trx+8T4BYK275PJpvt1/epvG7dLh/HaDKuJsya1jzZUplnV7ZSdnd8LsMBi1jQpi4SIBELcnTzOPZK1lMbKDOpLqhFonoTLQKB/LJlLD3pJO8Xl+t6vwzDYQKJBpx02vslYwyIeL0RsP0XuQh2LOb94sxkcLPeX/kiAqEQxnW3/TK0QiHshgbcXA6n+41oJ6xwGKuuzuv1yGS8N2wRpLKcjVsoeNtsGzsex5TLb/2mFQxi19ZiXNf7/7v/yg6Hwbbf/KZl27h99AbscE7b9no+huJNKxLx3oz68abldnV5z8ku7PJNq8c5B/ymFYvhOs6b37SM8Yrq7uIi1bL9TWtXBUsg4BWpsMPrpK/iws3l+vl9qhSOgcD25xSQYgkJhbzvU6kExmBEoFjaXjDshLEsJF/AFQu3uP347sfb/bo33R+OUylKd3H1QcT7ObNtsO1tX2PbEAh4r/lg0Ps3HELCYa/ILRXBcbc9f71ZkQhWXR24Dq4x2352JBIBAQmHMIUiViSM3dCAMQYJ2Jhij9ep3WMwvoj3B0SlWLfq6jD5PGJbby7W7QCBpiYkFKLc0uK99no+5GBw++eVc0ow6P0BUCggeH8YWNGo95rM5xFjsCeM9157rW242Sy978IpwdC258NOxLFiNZTTadzOzm2P3YrFdugJtONx7Lo6KBb67AkEIBjyCvD6Oqy6etxsF05Hp/fzEwwSaGr0nuOuLu8PrdpagvEErnG3t+uDVVuDFY1hNSSQcplyR4f3XAUCBBoawA5g8l5PIJEIwfo4RsBpa8eSnSwJa9tYxx/f5/83Gt145qH86cU1bOgssGRjB2vacqRzDulCmULZ3fbRliuxpi3HlnSeeDjEslQGY7wicXxdhAMn1/OVY6YSj8f9fkjKZ1oEgvdmnenaYZNbKlPuo6eqr7cvJ53GlMs7tDUAvY51i0Xo6PB+se6iCACvyJNwGNO5/fKWMcbrGemZx3G8nqJQ6E0FYO9fi6ZU8i7TVArIbduLRe+Nv8ebsymVdnhz8k7oFcGmxxu5ofImVnmj3/4AXEy57G3rPq+IVxxWig4ALMt7wzdm23m7C2CJRr3etMrzKMEAEo1huQ6m7CCAVVODPS7p9bBkc14PVSi0rVjHsjAiWGGvJ8xONlHKZqGtDcTa1kPkGoPp6vJ6V2MxQhPGez2B7R079v5texwhr4Csq8OO13s9F61t3vettgarpqbSW5Hbdvkq2NSEhEOUWlq3FbY7nDYa9WqeaBS7ocHrCexMY7JdYNsEGhu3XWY0+TyEwgTj9UgggCkUtv0B8iaWhdXQ4L252zZuVxYn57W14nGsSNR7/RYKGIFgfT2WHaBUyOMa+n6tCt6lw7o6LNvGLZcpFwoYx8EOec8N8TiO42CMIRgMYlneOKZSqUQikdj2tVKq/046eMpO9z2zJMWKljaWb+6iPV9i9/oQLdkya9tzZEsOhbJLKlNg4etbeWltG985eT+aJ2rv4FimYwKBAxqbzB+OP94rmCq9FRKNEmxspJzP47a2ekVLpXcLp+z9RVoqe3+5NzYigNPWVuk5cr3xMd29hmVvvEt3DxMiuO3tuKUSxrK8njhjvHOWy2DbSE0NdiiEm07vuoer+9JrJIJJp3FLJa8XJhDAClT+cq70ThGJEIhGMbkcTh+XYXc4bXdx1N3DJgLBIHYg4BULlSJSwmHsWAxTKOB0X1Lb2Tm7e9gKBW+8UGWbhEJQKnnZ8YpKKxYD1/V6ON6iJ9KurwfHoZxO73A5EsfxCm9jvOw1NYhleT2WPXu4RHY4HwG7cuk0hOno8HoJbRurJuaNVcvmwHUgGEISCexwCJPp8p5Tx0GCPf+2qvRqVR67nUh4r5OuLggGsOrqEcvyCjjwLts3NBBMJLxxWek0ViSMVR8nuPsk7HHjYPwEmDlDL+0opQbkuRUpXt/Syc/+uopMwSEasrnmQwdxRLP+LhmFBmVMoBaBjNyJIZlMBsdxKJfLXi9ZKrV9Z7u3dLwEvDFEVrmMUyjgZDJIoYAVCiEITqd3SUUCXsFhMLjZLE46jdtrEP22nkzbxg6HkZoYTqYLd8tmjFW5fBQMeufM5bziKVZDoLYWk89T7i5uKnboQRPLy9DUiOnK4lQei1R6jUyxiJvLYUSQUJBAfT3GcSi3tG6/fNZ9qZFKD6fresVqPA7lMk5rq/c8BQJYwaDXS1oprCUQwKqv94r1zs5dXubEtr2CMRDw2nYXgaEQBrw8xni9mLW1Xk9Yj0usfapcIpZo1Lt0Wy7vuliPRHBzuZ2OyYPK5f1EgsC4cV6hXBlYHxg/ntAekwnU1eHWx5Hp02iaNm3n2ZRSY8ozS1Jc9r//ojNfJhKw+f5pB3Lo9HF+x1IDo0XgYBmpRaAaeqmehXVLCwBWMEgwEsF0duLkCzjFAlIoejNJXRcnn8ekO5FgkEBNDZYxlDMZnM5ObwxfD262UlgHA97kiFCIUlsbdM9mjcWwLItySysUixAIILEodjSKSWcotbRsm1Xp5nOYfMErmItFr+jrHueYyXjbehaW3T22lRmrgcokElwXKmMZ7doarGSSyJQpJE8/fTiecqXUCJBKpfjEnS+Tzpdpqg2z4AP7aCE4umgROFi0CFSjUevrr2Ol05hNm8ivXUt5zRqcbM4bl1go4LS2ej2+3bOgew4F6OPnPvHRjzLxS1/04ZEopfzQ0dHBV373Gmtauyi5hmvnHqyXhkcPXSJGqbGscWbvW2rvXGdnJ7J5M+VVq0g/9hiFxa96Yyi7x3CGwwQPOWSIkiqlRqJ4PM5/n7IPn/zVi2zoKDD/gZf4bfIgHW88huj0PKXGAKelhfZf38XWH/2I3D+fx1TGPYanTWPceZ9nn0cfIfnew/yOqZQaZvF4nGtP2Y9oyKar4HD23Yv9jqSGkV4ORi8Hq+q15e67yT79V/KrVnmTeVwXKxKh5qijiH/sNOpmzfI7olJqBHhmSYqL/vACZcdw+PQGfniGXhkY4fRysFLqzcrlMuu/8U3yL72E2z1TPBAgMGki0dmHEPnUJ2lqavI7plJqBDmiOcmp75rIPYs28Nzqdv6waDUfnj3N71hqiGkRqFSVaf3Nb8guWuTdVi8SITp7NvUfO43EQQf5HU0pNYJ99cT9eX5NB0u3dHHTkyu1CBwDdEygUlUmfPTRSE0NUlfHpJtuZMo139ECUCnVLzd8ZBbhgEVnvsx3HnjJ7zhqiGkRqFSViU+dSmy//bCAlh/+yO84SqlRJJlMcvy+3nqB/1jdQbbXrUpVddEiUKkqVDdnDm6hQO5f/+L1U07dcVFspZTahcs+dBCz90ywW0OM//nXZr/jqCGkRaBSVajhfcfT8KlPIoEApc2baTnjE6T+9pzfsZRSo8QpB09h+eZObv3LCpZs1D8iq5UWgUpVqQlnnklywQLvtnLZLC2XXqqFoFKqX47ddzcitkW+5PK9h5b7HUcNES0ClapiyeOOpemH/40dj+Pm87Rcfjmp51/wO5ZSahQ4+6hpACzelOHlNR3+hlFDQotApapcsrmZxu9/D7u21usR/OY3aWlp8TuWUmqE+/Dsaew5LkbAEm59SnsDq5EWgUqNAcnmZhpv+AGEQrjZLM6DD/odSSk1Csw/6QASsSCL13fw0Mtr/I6jBpkWgUqNEcnmZoKJBFYkQmDvvf2Oo5QaBQ6YEmdaMka25HDjoyv9jqMGmRaBSo0RbatW4ZbL2A0N1L3nPX7HUUqNEl88djoBW0h1Fbnh4Vf8jqMGkRaBSo0BLcuWsfVLX8Lt7MROJAgGg35HUkqNEs0Tkxw1oxGAh17ZQkeHThKpFloEKlXl2levpvWSS3A6OrHCIRJfON/vSEqpUWb+CTNJ1oYQC25+ep3fcdQg0SJQqSrXMm8+5XQGO15P05VX0aD3EVZKDVA8Hufso6bjuPDgKxt5ZokuIF0NfCsCReRUEfmjiKwXkYyIPC8in+jVxhKRS0VknYhkReQJEdmvj3PtLyJPVtqsE5EFIiLD92iUGpnan/s7pZYWrGiU8TfeSNOhh/gdSSk1Sp08eyqJSIBi2eXbf9KxgdXAz57AC4FO4CvAh4DHgTtF5Ms92swH5gHXAnOBMvCoiDR0NxCRJuCxyr65wDWV474+DI9BqRHNyecxuRziuiSmTfM7jlJqlLv65H23TRL58l3/8DuOeof8LAI/aIz5hDHmt8aYx40xFwH3AF8FEJEocDFwjTHmRmPMI8DJQBQ4t8d5Pg+EgI8YYx4xxtyEVwh+Q0TCw/mAlBppTLoTjMEpFOh4Rf9yV0q9M80Tk5z6rokAPLuqjTsX6iLSo5lvRaAxpq8BBS8A4yqfHwHUAff1OKYTeBg4qccxJwJ/ruzrdh8Qr5xDqTFLJu8Btg2Ow+YFl7L2sstpe/FFv2MppUaxr564P7PG1wBwy1OrdXzgKDbSJoYcAbxa+XwW4AC9/8xYAjT3+HoW8FqvNssrxzaj1BjWdPBBjL/ySiIHHICIkHvuObZ89UJWfvazpP72nN/xlFKj1A0fmUWyJkTZMXzzjy+xZKMWgqPRiCkCRWQO3pi+myubGoCMMcbp1bS9so8e7XZYtKhyTKZXO6XGpKZ/O5rpP/4RE674FsGpU3HLZQpLl7H1kktY8fHTWTtvPqm77yaV0l/iSqn+SSaT/OBj+xMN2eSKDl//7WKy2azfsdQABfwOACAi04G7gQeAX/bYZfpq3sf2/rbr+X+eA5wDMGXKlAGkVWp0qj/wQOpv+wmpRx+j42c/o7hu3baPzNNPwy230ppIEJ45k8isvQlNnQrvfS8NDfq3lFLqzZonJrnmQwdxyf+8SEeuzO3PrOHLc/QC3Ggixuy0ThqeACJJYCFeb94xxphsZfv5wI+AcM/eQBH5FnCeMWZC5estwE3GmCt6tLGBAvAlY8wtb5Vh9uzZZtGiRYP3oJQaBVKPPU7usccorluHs2ULTlcX2DZWNIoVCuF0dGAcB7shQXDqNMJTpxLefz845BCampr8jq+UGiEWrUpxx7NvsKkjxwETa7nsvw72O9JYMCjL4PnaEygiNcD/VnJ8oLsArFgK2MDMyufdmvHGBfZs1/tPj70qxy5BKdWn5HHHwnHHbvu6paUF+etfKa7fQGHVKgqLF+Ok0zidaczSpRSXLKHjj3/02ibihGc1E56xJxKNYjc1wcyZJJu1F0CpsWb29CSduTLz/vAyb7RkCdgvMe+Duij9aOBbESgiAeC3wAzgSGPMll5NFgJp4FTg25Vj6oATgO/2aPcw8DURqTPGpCvbTsHrWXxm6B6BUtWlqakJ5s7dYVt3YZh/fQWF116juH4dbqYLp6OT/KuvUnz1VZxcDlMqAbBVxOtJjEWx6uq9+xRPmkhk//0JTpmCmTqVeDzux8NTSg2hY/fdjcNeWMuzq9r447+20BBbynnHzfI7lnoLvl0OFpHbgM8BXwb+3nOfMeZvlTaXAt8ELgKWAZcABwD7GmPaKm2a8GYUvwRch9cLeD1wpTHm2v5k0cvBSvVfS0sL8uyzFF5/ndKWLRRXrsJpa8PN5aCv3ye2jR2LIcEgTmcniGDX12Mnk9j19UgohN3URGSvmQRmzSK+zz7D/6CUUoPirJ8uZNnmLhqiAT5x+DTOOHxPvyNVq0G5HOxnEbgamNrXPmOMVNpYwAK8BaEbgOeALxhjFvc61wHAj4HDgDbgFuDbpp8PTotApQZHaskSWLGC0spVlDdupNzaiikUwLYxmQylzZsx5TJi20g4DI6DWyh4B1sWVm0tVigE5TJGBLumBiseJ9DYSGDiRIJTp8Dee+tlZ6VGsEvufYEX17WSLRnOPHQy5x6rP69DYHQXgSOJFoFKDZ+WRf9E1q+juDVFeeMGSmvWepNQSiUkEMAUizjp9LZLzH2RYJBgUxNWQ4NXZAaDBMeNw540ifCkSTBuHGa3CST33nsYH5lSqtt5v/o7/1zTjiXCpw/fQwvBwadF4GDRIlCpkSOdTlN65RVKS5ZQWruOciqF096Om+7Ezea8y86WhV1Tg9g25Y4OcF2vJzEUwoBXGBoDwSB2UxOBhgbEuIBgx+ux6uMEd5+EPW067D5JexaVGgJn3raQpVu6sES48LjpfPS9M/2OVE20CBwsWgQqNbq0t7djr1mLm9pK1+LFlNeuxRSKXi9iVwZn61acbA5EkJoabNvG6erCFIt9n7AyoSUwYQLB3XbDzeexQkHs8eMJz5hJYPJkzPhxNEyfPrwPVKlR7lM/XcjyVI6GaJBPHzmNk2f3OQpMDZwWgYNFi0ClqlNHRwfyxhpMWyvZV1+jtHIlTjqN29GBm89t71k0xutJjEaxwuFtayRiWRAMIpaFyee39S5a0SjBceOoO+F9jPv4x/1+mEqNaJfc+wKL3mil4BjOOkzHCA4SLQIHixaBSo097e3tsGYNZu1auha/SnnFCpxCAbJZym1tuPn8Lo+XQIDQ7ruz56/vHKbESo1eZ/30GZZsymCJ8PUT9uLDs6f5HWm0G/2LRSul1GBLLVkC6zfgrF5FuaUVt60NAOM4lLZupbx5s9er57pIOIxbLkO57H0dCGDV14OI91H5I1kCASQUQiIRJBzCrqsnOGE89R/+sJ8PValR41efPYJTbnqaNW05fvDo60yoreWI5qTfscY8LQKVUiNed6+du2ULhddfx9ma8i7XYrzZxakUbi7nLTfTx9UNCYexolHcYnF7ARgIYEWj2JEIAFZtDYHxEwjvsQc0JAiEw5iJk2DmDJJJfbNS6p269bR9OOOOl2nLlbjsf//Fb5IH6s+Wz/RyMHo5WKmRJJ1O4z7/PJnnniO/+FXKmzdhiiUkFMItlbb12lmhEFZtLW42u+OlW8vCikSQSBg7VoOVSBDefRL2njMI1tZgEg2YCeNJztK7GSg13JZsTHHunS+RKzpMrA9z/wX/7nek0UrHBA4WLQKV8kdqyRKchQspvPYapfUbcCrLvUgoRDmXg0rPntg2Vn29t8A0YMVihCZNIjR1CtTVYQdDsMdk2Gsv7VlQaoR7ZPEGLv/jYsqO4fDpDfzwjEP8jjQa6ZhApdTokd6wgfwzz5D7xz8orVlLubUVN5t9Uzu7pobA+PHY9fUE43Gie++FdcIJNDY2+pBaKTXYjt9vEq+sbeWeRRt4dlUbTy9ew9H7TfE71pikRaBSatB1dHRQ/L//o+vZv1XW8PPu6uGm096SLN0q9xEOTBhPaOo0Iu86GPuoo2hoaPAvvFJqyH31xP1ZsTXLki0Z/vBKC/tMTpKMx/yONebo5WD0crBSgyX1wAOk//QnSuvW43R17TDrNpBMYsViSDRKcNIkIgcfBEccoZdvlRqjstks1z2yguffaGFcXZSfffowvyONJno5WCnlv8ymTXTecw+Zx5/AaW/37tIRDhNobCQ4YTzBadMJv/e9JI860u+oSqkRJBaL8ZkjZ/DIq5vZ1NnBTx5fogtJDzMtApVSA5ZKpSj99nfk/v53776+xSJuJoMEg4Rn7En8ox+l8fjj/Y6plBrhpiRjzJ4S59lVbfz67+v52HsmEo/H/Y41ZmgRqJTqt9TfniP9y19SWLoUUy57s3Zra4nMmEH4XQcTmjtXJ3AopQbk0hOm8/FfpkkXHa57ZDVXn3KQ35HGDC0ClVJvKbVkCen//iH5V1/dthizFYsRe/e7qD39dBoOOMDnhEqp0SqZTHLhnJn86MmV/OONFAuXbeDIvSf5HWtM0CJQKbVT7e3tpL9/PV0LF2LKZQBCkycT/9AHSX784z6nU0pVi5MOnsL9L23k+bUdXPPQMv5Xi8BhYfkdQCk1MrX8z/+Q+uxnyTz1FKZcxo7HGXfR15hxz91aACqlBt1X3jcD2xK2ZIrc8PArfscZE7QnUCm1g9QfH6D9zjsxhQLGcbATCepOeB8Tv/AFv6MppapY88Qkh05N8OyqNh54ZQufObxDJ4kMMe0JVEpts+7yy9l6/fWUNm3CGEN87oeYdPddWgAqpYbFpSdMJxELgmu4+el1fsepeloEKqVoX72a1Z87h/RfFwIQbm5m3E9vZ8I551BbW+tzOqXUWJFMJvnk4VMJBIRHXt3Eio0pvyNVNS0ClRrjUnffzZbPn0d+9WrsmhoaP/VJ9rz9NhomTPA7mlJqDDrj8D2JhWwyRYd597/md5yqpkWgUmPYmku+ztZbbsXp6iIYjzP+yiuZ8OlP+x1LKTXGfeE/9gJgVSrHHxat9jdMFdMiUKkxKL1xI6sv+DJdzz8PlkXkgP0Z/4ufkzhYF2lVSvnv+P0mMWt8DQC3/mWVz2mqlxaBSo0xrU8/zeYLLqC4ejUpc2m7AAAgAElEQVSBRIKm889n+s03U1dX53c0pZTaZv4HZxGwhbZcie888JLfcaqSFoFKjSEbf/Rjtiy4lPLWFHZtLbt952rGn/ZRv2MppdSbNE9MctQM7zaUf1q8lVRKJ4kMNt+KQBF5v4g8LSIpEcmLyDIRuUxEwj3aWCJyqYisE5GsiDwhIvv1ca79ReTJSpt1IrJARGR4H5FSI9vqC75M+29/i3FdgrtPovGmG6nbe2+/Yyml1E5d99H3UB8JUCi7XPrQ637HqTp+9gQ2AE8BZwMnAT8FvgH8oEeb+cA84FpgLlAGHhWRhu4GItIEPFbZNxe4pnLc14f+ISg18qWWL+f1U04l98ILANQcdigz7rqLhoaGtzhSKaX8d+q7vVvIvbA2zZI12hs4mHy7Y4gx5q5em54QkXrgSyLyRSACXAxcY4y5EUBEngPWAOfiFYYAnwdCwEeMMZ3AI5XC8BsicoMxpjAMD0epESn12OO0XHcdbjaLBAI0fOIMJpx9tt+xlFKq3849tpmnlrewqTPP/a9soXlK0u9IVWOkjQlsBYKVz48A6oD7undWiryH8XoOu50I/Lmyr9t9QLxyDqXGpE0//zlbr7wSN5vFqq0ledWVWgAqpUala05+F/vunmDRG63c+dQKv+NUDd+LQBGxRSQqIkcCFwA3GmMMMAtwgOW9DlkCNPf4ehbQezXJ5ZVjm1FqDNr4gxvouO+3IEJo8mSa7vwVyaOO8juWUkq9LVOSMQ6aUscbLVluWbiSJXonkUHhexEItABZ4K/AX/DGBYI3ZjBjjHF6tW+v7KNHu46eDSrHZHq1U6rqtbS0sPLMs0g/9hh2NEr8+DnMuOdukkm9fKKUGt0+e9TeNESDlB3D1Q8s9TtOVRgJReAxwNHAhcAHgVt67DN9tJc+tve33fadIueIyCIRWbR169YBBVZqJGpra6Pt/C9QWL0ayiUaz/4Mk+bP9zuWUkoNmotOnIUlwopUlt8vesPvOKOe70WgMeZFY8xfjTE34F0OPkdEZgBtQJ2I2L0OieP1BnZrAxI9G1SOqenVrvf/e5sxZrYxZva4ceMG46Eo5ZtisUj7xZdQ2rQJCYVo/MpXaPzAB/yOpZRSg+r4/SZx4v7jqA1b3PHMKtakOt76ILVTvheBvbxY+Xc6sBSwgZm92jTjjQvstpQ3j/3bq3LsEpQaAzZ88UsU167Fqokx7uKLSJ500lsfpJRSo9C35h5MQ02YTZ1FLrz3X37HGdVGWhF4ZOXfVcBCIA2c2r1TROqAE4A/9TjmYeCEyr5up+CNE3xmSNMqNQKs/OSnyC9dCiJMuOACmk480e9ISik1pC48vhlLhDVtOb2l3Dvg5x1Dfi8iX6/cOeR4EfkWcD1wrzFmhTEmD3wPmCci54vIHOD3QB74SY9T3QqUgN+LyBwROY/KAtO6RqCqdm9ccAGFlSsxQPJznyWhBaBSagw4bEaS9+3jDeV64OUt/H2Vju1/O3xbLBr4B3A6sCfeBI4VeMVbz4khV+NN8FiAN9P3OWCOMaatu4ExpqVSIP4YeABvjOA1wHeH4TEo5Zut995LbslSsG0azzqTppNP9juSUkoNmys/cjDLt/6VDe0Frnt4GT87PUQ8Hvc71qgi3pJ8Y9vs2bPNokWL/I6hVL+1/WMRLd/9Lm6hQN0xxzDxaxf6HUkppYZdR0cHZ9/9Cls6ckxKhPnN54/2O9JwkcE4yUgbE6iUeguptWvZeuUVOO3thJtnaQGolBqz4vE4l5y4N0XXsDKV45L7/ul3pFFFi0ClRpn2r16I096BRCIkFizwO45SSvnq0Onj+NAB4wF4clkLdy7sfaMxtTNaBCo1imz88Y2UNm9GAgEaL7lYx78opRQw74MHccBEb5GQ255+Q28r109aBCo1iqT//GcAYoceQvLf/s3nNEopNXJ8d+5eNESDFMouF/7mFb/jjApaBCo1Smy+/Xacjg4kECB28cV+x1FKqRElmUxy+X8eQMAWUl1FzvvV3/yONOJpEajUKNFx/x8BiL3n3SSTSZ/TKKXUyHNEc5JPHrYHAP9c08mtj2mP4K5oEajUKLD113dt7wX8xjf8jqOUUiPWucc2c/j0BkTgLys6WJPK+h1pxNIiUKlRIP344wBEDzpIewGVUuot/PCMQ/jIuyZTExC+df9LZLNaCPZFi0ClRrh0Oo3T0Y4Vi1F3+ul+x1FKqVHh7COn0dJV5pWNac799Qt+xxmRtAhUaoQrPPoolB2CkybReOghfsdRSqlRIRmPcdbhU7BEWLqli2sffMnvSCOOFoFKjXDpBx/E6eggOHmy31GUUmpU+fDsaRwzqxGAhxen+PuqrT4nGlm0CFRqBHNdl8LKVRjHIXrkEX7HUUqpUefaU97DgXvEiYZsrn7wNTo6OvyONGJoEajUCNb+4IMY18WqrSF54ol+x1FKqVHp+v9qRsRiY3ueL96ry8Z00yJQqREsu2gRVihEeN/9/I6ilFKjVjwe57yjpwGwdEsXP3l8ib+BRggtApUawQqvvgpAzSGzfU6ilFKj2wffPYX3TKkH4J5F6/WyMFoEKjVipRYvprh5C6ZUIqiXgpVS6h27du4+jK8NY1vCZQ8t9zuO77QIVGqEyj/wv2AMgYYGEomE33GUUmrUi8fjXDZ3X0IBYdHqdm5/YmxfFtYiUKkRqrh0KVgW4QMP9DuKUkpVjUOnj+OAiXFKjsuv/7GeVCrldyTfaBGo1AiU3rgR09mJXVdH/Rl6lxCllBpMlxw3hWjIJld0uOrPq/yO4xstApUagTK/+Q1ONkto6lTie+3ldxyllKoqyWSS02fvDsBzq9tZsnFs9gZqEajUCJR54kncXI7oAfv7HUUpparSucc2k6wJ4RrD1Q8s9TuOL7QIVGqEST3xBOW2NjCG4Ec/6nccpZSqWl99396At3bgHxat9jeMDwI72yEiK9/BeQ3wCWPMs+/gHEqNSem77wEgOHkyjY2NPqdRSqnqdfx+k/jVwlUs3dLFL595gw/PnuZ3pGG1q57AaUAeeGOAHxsqx0aHKLNSVSv1wgvkl3qXJeKnneZzGqWUqn7zPzgLS4SNnQX+tWqj33GG1U57Aiu+bYy5eyAnFJEksOXtR1Jq7ErfdLN3GXi33Uh+6IN+x1FKqarXPDHJe6bGWb45w8JVaQ6cPtHvSMNmVz2BS4G3c0+VUuXYroEcJCK7i0hGRIyI1PbYbonIpSKyTkSyIvKEiLzpRqoisr+IPFlps05EFoiIvI38Svkiteif5JctAyBx1pk+p1FKqbHjguP3ZkpDhCeXbSXVkfU7zrDZaRFojNnHGPPgQE9ojOmoHPvcAA/9HpDpY/t8YB5wLTAXKAOPikhDdwMRaQIeq+ybC1xTOe7rA82vlF/St966vRfwg9oLqJRSw2XWhARlI6xKdXHp/f/yO86wGRGzg0XkaOBE4Pu9tkeBi4FrjDE3GmMeAU7GG294bo+mnwdCwEeMMY8YY27CKwS/ISLh4XgMSr0T7YsXU1i1CiyLxKc/5XccpZQac45rTgLw0vr0mLmLiO9FoIjYwI+BK4Hez/oRQB1wX/cGY0wn8DBwUo92JwJ/ruzrdh8Qr5xDqRGt/bbbsSIRovvvR/L97/c7jlJKjTlnHrkX9ZEAZcfw6+fHxtSGt5oYsoPKJdizgcOABt5cRBpjzHEDzPB5IALcBJzRa98swAGW99q+BPiPXu0e7dVmeeXYZuCJAWZSathsvP12CkuWYNfX0XThhX7HUUqpMevYfZL83ytbeWlt51s3rgL9LgJFZCqwEJiEN2GkHmhlezGYYuCTQZqAq/DWFCz1MY+jAcgYY5xe29sr+3q222ESizHGEZFMr3ZKjSgtW7bQed99uPkCtcccQ92MGX5HUkqpMesTh03l6aUpVqa6eH5FinfPSPodaUgN5HLwt4EEcBywFyDAaXjF4DVAGjh6gP//1cBzxpiHdtHG9LFN+tje33beDpFzRGSRiCzaunVrv8IqNdgy374at1jCrq9n929+w+84Sik1pk1JxomFA+RKDnf8bZXfcYbcQIrA44DbjTFPsL2wEmNM1hgzH3gZ+G5/T1ZZ5uUzwBUikhCRBBCr7I5XJoW0AXWVcYM9xfF6A7u14RWoPc9vAzW92m1jjLnNGDPbGDN73Lhx/Y2t1KBpf/JJb2HoQICm88/zO45SSingqBnenZpeXF/9l4QHUgQ2Aa9UPi9V/u15V5BHgOMHcL69gCDwLF4R14Y3LhBgHd5kkaWADczsdWwz3rjAbksr23qf3+7VTqkRobW1ldQPf4QBat7zbpr+8z/9jqSUUgo4c/ZuBGwhV3R46OU1fscZUgMpArcC3TcyTePdUm5aj/0hBnaruL/iTe7o+dHdk/h+vHUDF1b+r1O7DxKROuAE4E89zvUwcEJlX7dT8MYJPjOATEoNi/avXUQplcKybeovvtjvOEoppSqSySQHTIpTE7J56Y2033GG1EBmBy8GDgJvCrCI/B04X0T+iFdMnsMAet2MMSngyZ7bRGRa5dOnjTGZyrbvAfNEpBVYBlyCV4D+pMehtwIXAL8XkevwegHnAVcaYwoDeIxKDblNt9xKYcUKEKHxC+eT0OEISik1ohw6vYHXt2Z4ecPbuXHa6DGQnsD7gcMrY/XAW9dvL2AVsKLy+VWDGw/wJo9cCywAHsDrcZxjjGnrbmCMaQHmVPY9AFyKN1ml32MUlRoO6cWL6Xj4YSQUov6E95HUy8BKKTXiHL/veDKFMiu3dlX1wtFiTJ+TZ/t3sMhs4HS89fj+YIwZlZdeZ8+ebRYtWuR3DFXl2traaPvSlyinM0RmzmDK9df7HUkppdROzPn+43Tmy5x9xBTOPbb3tAPfvWlNvbdjQItF92aMWQRo9aRUP7RfdDHFdesJNDXSMG+e33GUUkrtwtSGKC9vTPP3VW073Ke2mrxlESgijcCn8WbopoB7jDGvDnUwparJ2ssvJ79sGYjQ8NnPUtfU5HckpZRSu3Do9AZe3pjmjbac31GGzC7HBIrIZOBfwHXAucB84AUROWEYsilVFbbedReZJ54EoHbOcTSddNKuD1BKKeW7kw/07haSKThVOy7wrSaGXA5MxFu/74PA14As8KMhzqVUVUj97W+0/OznYAzhWbPY47LL/I6klFKqH5LJJNGQjWsMT6/O+B1nSLxVETgHuM8Yc4Ex5kFjzA3AhcBMEdlz6OMpNXqlUilar7wKUyoRaGpiwg//2+9ISimlBmB8bQiAZ1e2+JxkaLxVETiJXmv5AU/gzUqZNBSBlKoWHRd8GSedxopEaPjutdTU1PgdSSml1ABMT3p3s11XpeMC36oIDOLdsaOn7q9Dgx9Hqeqw/qqrKK5fj9g2TV+7kOSsWX5HUkopNUDH7zeOaNAiErD9jjIk+rNY9M4WEnz7CwwqVcW23HEHXYv+iRWJED/jdJInnuh3JKWUUm/DkdOT7Lt7gnA4SDab9TvOoOvPOoHXisg3e3xt4xWAPxWRrl5tjTHmoEFLp9Qos/nnP6f9rruxolHqP/Jhdvvc5/yOpJRS6m2KxWKsTXXRmi3y8uYODpse8zvSoHqrnsA1gAvU9fiIVbZbvbbXAfVDllSpEa7l/j/SesevcItFQvvsw27nVuvyokopNXa4xlB2Dc8srb7JIbvsCTTGTBumHEqNau3PP0/qJz8BIDR9GlO/d52/gZRSSg2Kppogqa4iq1uqb3JIf8YEKqV2oXXJElJXXYVxHEJTpzLh5pv9jqSUUmqQ7DmhFkuEbMn1O8qg0yJQqXcgtXEjLV+7iHJrG4FEguQNP6C2ttbvWEoppQbJgbsnaKgJUqjCIrA/E0O2EZHDgS8CewFNeOsF9mSMMTMGKZtSI177xZdQ7uzEikRovPIK4npPYKWUqip7j0/QmSuRK5b9jjLo+l0EishZwC+AErAMb3KIUmPWhu9fT3nTJiQYpOnSBTTqWoBKKVV1DpgSp+QYSo4hlUqRTCb9jjRoBtITOB9YCswxxmwYojxKjQqpe++la+FCrFiMxOkfJ/lv/+Z3JKWUUkMkHLAolF2WpaCKasABjQmcCtyiBaAa61JP/oWW23+K6eqi9tj/YPzHPuZ3JKWUUkMoEvTKpc3pTp+TDK6BFIHrgPBQBVFqNEht2EDLNdfgFgrYySSTvvIVvyMppZQaYpMbosQjAbJVNixwIEXgrcAZIlKdN9BTqh/av/wV3GwWKxYjfv33/Y6jlFJqGMSjQYqOYW1b7xuljW47HRMoIr0HOS0CTgb+LiI3AasAp/dxxpinBjWhUiPEuiuvpLRpE2LbNM2bR3LiRL8jKaWUGgalskuu5PBGaowUgcCTePcI7ql7SZif7mSfwbu3sFJVJbVhA5kn/wJA7Zw5JP9dJ4IopdRYURP2SpvWbMnnJINrV0Xgp4cthVIjXObb38aUStj19UxeMN/vOEoppYbRxHpvSkS+WF0LRu+0CDTG3DGcQZQaqVqee47cK4sBaDjrTJ/TKKWUGm7JukoRWGV3DRnQHUOUGos6f/ZzxLIITtmDcaed5nccpZRSw2y3RA0Aeae6isCdzg4WkSNEZMD3wBIRu3Js3TuLppT/2h56iPLWrdgNCRouucTvOEoppXywR2MIgNJYKQKBp4ET3sY5E5VjD9lVIxH5lIiYPj4+26ONJSKXisg6EcmKyBMisl8f59pfRJ6stFknIgtEpPd9jZUasJaf/xy3q4uao4+mcf/9/Y6jlFLKB80TvduElJ3ec2JHt11dDhagSUSmDPCcjWyfRdwfxwK5Hl+v6PH5fGAecDHeLesuAR4VkX2NMW0Ald7Kx4CXgbnA3sD3gTJw7QCzK7VN59atlDZvARFin9Z5UkopNZYFbKHsGJZsTG0rCke7txoT+N+Vj4EaSKn8D2NMpvdGEYniFX/XGGNurGx7DlgDnMv2Au/zQAj4iDGmE3ikUhh+Q0RuMMYU3kZ+pSi++CJYFlYsSqKhwe84SimlfJSsCdFVKBMLxvyOMmh2VQRe8Q7PvfIdHn8EUAfc173BGNMpIg8DJ7G9CDwR+HOlAOx2H17+I4An3mEONUY569dDIECgmu4WrpRS6m2ZOaGeXMkh5xSB6igEd7VEzDstAvtrRaXn7nXgemPM7ZXts/DuSLK8V/slwH/0+HoW8GivNssrxzajRaB6m8obNmAFAoQmTfI7ilJKKZ+tTqXZ3Fmgo6MLJiT8jjMoBnLv4MG2EbgUOBP4APA34DYR+XJlfwOQMcb0vjVde2UfPdp19GxQOSbTq51SA1J84w1MoYDopWCllBrzuvIOJcewdGvW7yiDxrd1Ao0xfwb+3GPTwyISAi4VkR93N+vjUOlje3/bbd8pcg5wDsCUKQOd+6LGgvLmzRjHwdLLwUopNebZljfnNVvo3Tc1evnZE9iX/wc0AdOANqBORHrfiziO1xvYrQ1vWZptKsfU9Gq3A2PMbcaY2caY2ePGjRuE6KrauPkCiBCavIffUZRSSvksFBBEwDLVswLdSCsCu59ZF29JGBuY2atNM964wG5LK9t62qty7BKUeptCe+wBkQil5cv8jqKUUspnDbEwjTUhYrGg31EGzUgrAk8BUnjLwCwE0sCp3TsrdyE5AfhTj2MeBk7odYeSU/DGCT4z1IFV9YoddSSBSITC66/7HUUppZTPskWHzlyJANVz1xDfikAR+Z2IXCQiJ4rIB0TkTuCjwOXGGNcYkwe+B8wTkfNFZA7weyAP/KTHqW4FSsDvRWSOiJyHt8D0tbpGoHonYv/+75hikdKyZXSl037HUUop5aN0vkTJMbTnyn5HGTRva2KIiMwEJgCvGGM63qr9TiwDPgfsgXcZeDFwhjHm7h5trq7sW4A30/c5YE733UIAjDEtlQLxx8ADeGMErwG++zZzKQVA3bRpbAqHKbe20nnffdScfbbfkZRSSqlBM6CewEqP3Qq8cXhPAe+pbB8vIq+LyCn9PZcxZp4xZpYxJmaMiVYmadzdq41rjLnSGDOp0uYYY8ziPs71cmVftNL2KmNMdd3gT/kitPvuAOSf0ZEFSik1lrmVBUdC1TMvpP9FoIgcA/wBaMW7G8e2p8EYswXvnr8fG+R8SvkqevTRABTXrfc5iVJKKT81RoPUhGx2a6iOu4XAwHoCLwNeAg4Dbupj/7PAuwcjlFIjRfzD/4UEg7jZLOtvuMHvOEoppXxkWVXUDcjAisDZwF3GmJ1Ni1kH7PbOIyk1coQjEWoOPxyAzIMP0bJeewSVUmos2pQuks6XqR2jS8TYwK5m2yaB4juLo9TIM+Hyy7Dr63ELBdJXXuV3HKWUUj4oOV4fmO2OzSLwNeDoXez/AN7lYqWqSigUouEznwEg99prtD7+uM+JlFJKDbeS400M2buK7iQ6kCLwZ8ApInJ2j+OMiMRE5EfA4cBtgx1QqZFg3MkfIbLvvkggQPuv7kQnnyul1NjiVn7vJ6vofvL9LgKNMbcA9wK3A8sBA9yDd2eOLwK/NMbcNRQhlRoJGhfMJ5BMUm5rY+svf+l3HKWUUsNkycYUAJaM3YkhGGM+AZwMPIZ3X95W4CHgVGOMrqSrqlp8jz2ofd/7MLkcbXffTWrxm5asVEopVYVaK7fFCNpjuAgEMMb8wRhzsjFmP2PMvsaYucaY3w9FOKVGmuSnPoldW4ubL9A6bx7tW7b4HUkppdQQ25zJABAO+Ha33SHRr0cjIrUi4ojIpUMdSKmRLBAIkPjW5Vg1NTjpDC2XfJ2uri6/YymllBpCq1NeERiyx2ARaIzJAO2AdnuoMS954IGMv+IKAnW1lDduZPOXLqBUKvkdSyml1BCpCUWxRNi9Mep3lEE1kJL2CeDfhyqIUqNJw2GHkvzilzCOQ2H5ctadd77fkZRSSg2R2pjFu6Y2cOi0Jr+jDKqBFIEXA0eJyBUiUj9UgZQaLRqOn0PjmZ8AEfLLlrFm3ny/IymllBoCr6xrZ9nGDkqm7HeUQTWQIvAxIAIsANpEZJOIrOz1sWJoYio1Mo3/5CdpOOUUJBgk+9JLbP7FL/yOpJRSapAt25whXXBoT+/qxmmjT2AAbdfgrQ2olOphtwu+BOEw6QcfpON3vwMDEz7zab9jKaWUGiTpvNcD2FQT9jnJ4Op3EWiMOWYIcyg1qu127jmUN20i89RTtN5xB6XWViZf9DW/YymllBoEpcrdQmZNrq7RcNU111kpH026dAGRWXuD65K+/37WXn2135GUUkoNgrAFiajNno01fkcZVFoEKjVILMti6k03ETvoQAAyD/+ZNd+6wudUSiml3omOjg6yJYPrClOScb/jDKp+F4Ei4lYWjN7VR3VNm1FqgESEqTfeSM2hhwDQ9dhjvPHNeT6nUkop9XYt3lgiV3IoONU3LWIgPYG/6uPjbuC5yv5/AXcOajqlRqkp119P7ZFHApD9619546KLfE6klFLq7VjR0gZALFh9F0/7/YiMMZ8yxny618eZxpgjgKOBPYCfDFlSpUaZPa69htpjvPXVs/9YxNr5C3xOpJRSaqDaurxlYaJjuQjcFWPMM8AvgOsG43xKVYs9rrqK+rkfQkIhci++yPrrvud3JKWUUgOQLxuiQZs9mqprUggM7sSQ5cC7B/F8SlWF3S+6iMRZZyGhEJlHHmHlZz9HLpfzO5ZSSql+yJUM4aCQrA35HWXQDWYReAyg72xK9WG3Mz9BfO6HMOUyhaVLWXf22aRbWvyOpZRS6i2s2pqmPVvGMa7fUQZdvxeLFpGzdrKrEZgDnAT8dDBCKVWNxn/qUxCN0nrb7ZQ3bWbLl7+Mc/m3SOw10+9oSimldiJT8Iq/iC0+Jxl8A+kJ/CXeuL9f9vr4AfA+4GfAhQP5z0UkICLfEJHlIlIQkXUickOvNpaIXFrZlxWRJ0Rkvz7Otb+IPFlps05EFohI9X3H1Kg2/rTTGHfFFYQSCUpbtrLlwq+Seuopv2MppZTaiWzRAarvlnEwsHsH/0cf2wzQCqwyxnS9jf//F8BxwBXAEmASsH+vNvOBecDFwFLgEuBREdnXGNMGICJNwGPAy8BcYG/g+0AZuPZt5Pr/7N15eBzVlfD/76nqvdXavRsvgI0xNgRsIAHbhB3COBiCMcw7MwkEwvImZHmBEGISmwxLkknCMJAAGZaZ/IYMDOAJ+x7ABAhxtmExYIMXbGzJstbea7m/P6olJNmANUi01Tqf59EjqupW9Wk9Rn10l3OVGjIN8w7HHj2K7ZddRnFbCy1X/gDv619nzMK/KXdoSiml+sk4QRI4ujpW5kgG30D2Dn52MF9YRE4AzgAOMMa8/gFt4gTJ3zXGmBtK534PbATO4/0E73wgApxqjOkEniglhpeJyM+MMYXBjF2pj6t2+nTMddex/evfwNm+nfYbbsB2XRpPWVTu0JRSSvXi+0GR6DFVVWWOZPB97IUhIjJHRI4VkYGmyGcDT39QAlhyGJAC7u4+UUryHiWYg9jtBOCx0rVudwM1pWcotdupmzSJcXfcTmzmvphQiLY77mDrL3VarVJK7U6c0k4h0xvLHMgQGMi2cReLyAP9zt0JvEyQlL0iImMG8NqHAm+JyA0i0lmay3efiIzv1WYfwCMoP9PbG8CMfu1W92uzpnTvDJTaTSVratjjhhtIzp+HKRRo/8//ZP3/uxjHccodmlJKjXgtLS34xmCJ0NhYeVngQHoCzyAYhgVARI4qnftPgnl74wjm6+2qscCXgE+VnnMWMAdY0WtBRx2QNsZ4/e5tL12jV7uO3g1K96T7teshIl8RkVUismrbtm0DCFupwRUKhZj47W+TOPwwjOuSe/llNpx1Nl2treUOTSmlRrQ/NxUBCFfgymAYWBI4haAHrtsiYAvwd8aYa4GbgIUDeJ6Uvk42xjxsjLkL+HvgEOCoXu12tmOz7OT8rrYLGhtzizFmrjFm7qhRowYQtlKDT0SYeMUV1O7/AVkAACAASURBVJ58MlgWzoYNbDn7bNreeqvcoSml1AjmIgI18YGsox0+BpIEJoFsr+OjgCeNMd1J1uvAhAE8rw14xRjTu2Lu80ARmNmrTUpE7H731hD0BvZ+Vm3vBqV7kv3aKbVbG/etb9L41a8i0Sje9laaL/o6Lc88U+6wlFJqRPJ8m9njqjhyxthyhzIkBpIEbgb2BxCRyQSJWu8Vw3XAQFbhriboqetPgO6y3G8CNtC/mu4M+vZKvsmOc/+mle59A6WGkVGLT6Nh+TLsVAo/k6HlBz9gy69+Ve6wlFJqxHl7cxfr2/I4hcqcpz2QJPAB4HwRuQG4hyDhe6jX9VnA+gE870Fgtoj0nmm5AAgDfy0d/w7oAhZ3NxCRFHA88Eiv+x4Fji9d63YawTzBFwYQk1K7hVGHH87oX/yC8JgxmKJD+223s+WXvyx3WEopNaL8aXMbnTmXpmxlVpobSBJ4JcFw7YUECd83jDFN0FPP7xTgtwN43i3AduABEVkoIn8L/IpgiPl5AGNMHvgxcLmIXCgixwD3Anng5l7PuglwgHtF5BgRuYCgwPS1WiNQDVe1kycx8d/uIDZzJmJZdD3wIO9dfz3vz8BQSik1lDryLgCjkuEyRzI0djkJNMa0GWOOJph7V22MublfkyOAqwbwvE6CeYVtBCuMbyTY9eP0fk2vIigKvZSgNzICHNO9W0jpWdsJ9i+OlNpcAVwD/HBX41FqdxRLJpn08xtJLVwYJIK/uZ+3v3QW6ba2j75ZKaXUx9K9ZdzYmniZIxkaor0KMHfuXLNq1apyh6HUh2q6/Xbaf/X/4TsO9vjxTPjFz0nW15c7LKWUqlin3PAsnQWXa07dn0Om7laVRAalZs1AikUfIiLn9jt3soi8IiKbReTqwQhIKbVzY846i/oLL8CKx/G2baPpO5eT7ez86BuVUkr9r3TlPDCwT32k3KEMiYHMCfw+8PnuAxGZBPyaoOhzB/BtETlrcMNTSvU26rTTaLzkYkJVVRTXrWPLN79FPpv96BuVUkoNyNtbWugsuGSLPjU1NeUOZ0gMJAk8gGC1brczCLojP2WMmQk8DnxlEGNTSu1Ew7HHUnvOOYgxFN96i01f/SrFYrHcYSmlVEV5pzX4vRoJDSRVGl4G8s4agK29jo8HnjPGbC4d309Qm08pNcRGfX4hNf/nbxHbxln7Npu/8x1c1y13WEopVTHWNwfTbZLh/vtVVI6BJIHtwBgAEYkCnwae63XdAJW5fEap3dDYL32J2v/zf7ASCYpvvkWzlo9RSqlBk3N8IiGLMTXRcocyZAaSBP4FOEdE5hCUYIkBj/W6PhVoGsTYlFIfYey551D/93+PFQ6TfvwJNl19tSaCSik1CDJFQzJiMb46Ue5QhsxAksAfAOOAlwkKMT9pjOldV+VvgN8PYmxKqV3Q+LdnUnXEEfiOQ/rRx3j3iivKHZJSSg17bzV30pZ1CYUq9w/rgRSLfgE4CPgG8CVgYfc1EWkgWBjyi0GOTyn1EUSEMRd9jcSnPw1A5rnn2PSz68oclVJKDW+t6WBhSNQelJJ8u6UBLXkxxrxljPkXY8y/G2OKvc5vN8Z80xjz3Ifdr5QaGpZlMeHK5cTnzgEDXfffz9YHHyx3WEopNWzlHR+APep1OLiHiEwVkXNE5LsiMqV0LiIik0SkMqspKjUM2LbNuB/8gMgee4Dr0nHDjXS8+mq5w1JKqWEp7wVJ4NjaZJkjGToDSgJF5IfAW8AtwJXAnqVLMeB14MJBjU4pNSDRqioaf3Aldk0NfiZD0/IrybW3lzsspZQadhyvuyewcvu3BrJt3HnAJcCNwHH02rfOGNNJUCdw4c7vVkp9Umr22ou6b3wDKxLB27qVrd//Pp7nlTsspZQaNlpaWnC9YEHIjHGNZY5m6AykJ/BCYIUx5hvAn3dy/X+AfQYlKqXUxzLqmKOpXrwYCYcpvr6a1vtWlDskpZQaNt5qCb6HKnhRCAwsCZwOPPEh17cBlZsuKzXMjD73HJLHHIMkEnTccw9tL71U7pCUUmpYaEqnAQjblbtlHAwsCcwDHzY7cjLBriJKqd2AbduMuehrRCZPxmtpoXnZcjr+8pdyh6WUUrs93wt6AGtjoTJHMrQGkgS+DJyyswsiEgP+HvjdYASllBockaoqGr/5DaxkEj+Tofm663Acp9xhKaXUbi2ZCDFjfBWHTq4tdyhDaiBJ4I+Bz4jIr4D9S+fGisjxwDPAROCfBjc8pdTHVTV1KrX/90IkEsHb1kL2FS0bo5RSH2ZTR4btnUWsiA4HA2CMeRK4ADgNeLJ0+lfAw8ABwLnGmBcHPUKl1MdWd9hhRPbcE3Fdun7723KHo5RSu7VXNnSwLV1k7dZMuUMZUgMa7DbG3CIi9wOLgRkEZWLWAHcbYzYPQXxKqUEQSqWIHnAAxXfeIffyyxQ7O4lUV5c7LKWU2i2l8y4AyYhd5kiG1i4lgSISBQ4Fthhj1gD/MqRRKaUGXdXhh5P+7//GbW6i49lnGbVQy3oqpdTOeCaoEZiM63AwgAc8BZw4hLEopYZQ1Yx9CI8bB65H12+fKXc4Sim124qFbRJhm8m18XKHMqR2KQk0xrjAVnrtEqKUGl7seJz4UUcisRjuli3km5rKHZJSSu2WsgUX24aqeKzcoQypgfRz/hdwuohUdt+oUhWs7qijCO+xB+K65F5/vdzhKKXUbmlbpkhX3iNU4V1fA0no/hVIAE+IyEIRmSEik/p/DVGcSqlBEB43jtCYMfidHXQ88iimNO9FKaXU+7KOD0AkVNn9XgNZHfwqYAiGhD/7Ie0qeymNUsOYHYkQHjOGbL5A8Z238YpFQtFoucNSSqndiu8HfyCPqaoqcyRDayBJ4JUESeCgEJFngCM+4PJhxpgXS0PP3wXOA+qB3wNfNca81u9Zs4AbgEOAVuAm4Cqj3RxK7cBtbwPfR+IJrFBlb4mklFL/GxFb8I3FxMZwuUMZUrv8CWCMWTbIr30h0L9Q2ZXAgcAfSsffBS4HLgHeBC4FnhSRmcaYNgARaSBYufwKcDIwnWDnEhe4dpBjVmpY8z0Pr70DADtVhWVrx71SSvUXj4SIhA01lZ0D7nKdwFHAnkCLMebtwXhhY0yfWekiEgHmAncZY1wRiRMkf9cYY24otfk9sJGgZ7A7wTsfiACnGmM6CeYsNgCXicjPjDGFwYhXqUpgCgVMVxcAoZqaMkejlFK7p7ZsERBqKvz35IfOeBQRS0RuArYALwBvicjzpaRwsJ0A1AG/Lh0fBqSAu7sblJK8R+lbr/AE4LHStW53AzWlZyilSkyxiMlkQARJ6Y4hSinVX0tLC45n8PzKn1H2Uctevgp8haBG4H0EQ66HATcPQSxnAJuBlaXjfQiKVK/p1+4Ngi3r6NVudb82a0r3zkAp1aPY0YlfKEDIJjx6KP6WU0qp4a3FCb5blb0wGPjo4eB/IEiwPm2M6QIQkV8CXxKRWmNM+2AEISIJYCFwS6/FHHVA2hjj9WveXrpGr3YdvRsYYzwRSfdr1/81v0KQ4DJpkla2USOD39WJKRZBLEK1teUORymldjutpYzClgovEshH9wTuA9zRnQCW/AtBGZjpgxjHQqCK94eCu+2sL1Z2cn5X271/gzG3GGPmGmPmjhqlPSJqZPDaO8DzsMIhQnX15Q5HKaV2OxmvCIBlVX4S+FE9gUngvX7n3ut1bbCcAaw1xqzqda4NSImI3a83sIagN7B3uz5dGiJil+IblJ5KpSpFccN6JB6HcJjQlMnlDkcppXY76Vyx3CF8YnZlxPuDet0GJUUWkRqChR79ewHfJOhx3Lvf+RkE8wJ7t+s/929a6d43UEoBkN24ka7fPgMiJBcsIDFZk0CllPogYe0JBOBzIjK213GCIBFcLCKf6tfWGGN+NsAYTgGi7JgE/g7oAhYD/wggIingeOCHvdo9Cvw/EUn1GrY+jWCe4AsDjEWpiuT7Pu2/+Q3Ohg1INErqiAVaKFoppXaiKh4BwBkBq4N35VPgb0tf/Z23k3MGGGgSeAbwV2NMnxW+xpi8iPwYuFxEWoG3CIpF5+m7Ovkm4CLgXhH5EUEv4OXAlVojUKlAdu3b5F58CT+fJzJtGom9+3ewK6WUAgiVBjpt7QnkyKF8cRFpBI4GrviAJlcRDDsvJVjp+3vgmO7dQgCMMdtF5BiCBSsPEMwRvIa+vYVKjVie65J+6km87duRmhrqFi0iFI+XOyyllNotjW+MM7khQX0sUu5QhtyHJoHGmGeH8sWNMS3AB27KYozxCbaSu/IjnvMK8NlBDU6pCpF/6y2yf/ozxrJIzDmI6nmHlzskpZTarVliga3DwUqpYczL52m7bwXuxo1YjY3UHH88oUSi3GEppdRuK1NwWdeSJhGu/L3VR0A9bKVGrvQf/0juT3/CKxSITJpE8lP913IppZTqLZ8NvntGewKVUsOM57q4ra0U332XtvsfwOvowKqvp+YE7QVUSqmPYn/gJLXKo0mgUsOYMQY/m8Xr7MJp3U5h/Xrya9ZSeOcdTEcHXiaDxGLEDzyQ5MEHlztcpZTa7c1oDFMVDeH7hkf+spETP1W5W8tqEqjUMGJcFy+bxevsxG1rw21txd3WgtPcjLtlC8WmJrzWVvxMBisaJTJtGvE5B1F7zDGEtRdQKaU+Uk1NDZ+Z2sBLG1r4+XPvaBKolPrkGd/Hdxz8TAa/vQMv3YXX2YmfzeGl0zjNzXitrXjpLvx0Gj+dwXddsG1CY8cSmT6d+lNPITF9OjICNkJXSqnBcumxUzj55haackUuvfuP/Oj0OeUOaUhoEqjUbsAYA66LXyjgpdP4XV34XV14+Tx+Lo9fKOB3deF2diKui/E9vLZ2vGIB4xskHMGqCWH5HnZDI7G996bq0EOI7bVXud+aUkoNOzU1NZy8/xh+veo9nn+7lTe2tDBjXGO5wxp0mgQqVQbG8zCO837S19mJyeXw83n8QhFTSgj9XA7juphsBt/zMZ4LImBZSCKBbQkmmYRsFj8UJjJ6NLG99iJxyMGEq6vL/TaVUmrY+uYJs3jqjRaa00V+8tgafvklTQKVUgNkjME4DjgOfi6H19WFn8ni57JBIugGCSG+Hyz0cIrgeuB52LU1GM/DLRYRHMSKILEoEg4jmSx+Joyfy2KNGk28oYHozH2Jz5yJFR5By9uUUmqIfPdz+3DFg2+wqb1QkYtENAlUapAZ1w2Su2IRk8/jpdOYbA4/n8M4DsbzMb4Pvh/M1QvZSDSKHY9hHAdna1OQFHoexvMIJRJYUyYHw7/bW8A3REaPQUI2hU2bMfk84XFjic+aRWTCBJ3/p5RSg+Qz08cxf69tPL66mZ8+9bYmgUqp9xnfDxK70tCun81icjlMPo+fzwfXfAO2DaXCoxKNYEUiQc9fvoCfSWNcDzuZxKqtQtrakFAIjI+EQvilBBIEq7oGweDnc0QmTSY+ejTiOEQnTcLW4V+llBp03/jsZJ5avY2OnFNxi0Q0CVTqf8F4Hl5HB8ZxgyHczk78bAY/k8EgSCiE2DZ2TRKDQSwLKxZDYjG81lbcjg5M0cEUCvjFIlYiAb5PqKYamTwZP5cL6v9tb0VSqWD4N5/HD9mYrjR+axt+bS2JuXN16FcppYZQTU0Np3wqWCTy3JrKWiSiSaBS/wt+JoNxXBDAdZFwCGwbu6EROxFH4vGexM7P5gCQSDjoHfR9rEQSE3YgESdSWxskeaEQ7rZtGNclPGoUhc2bkZoa7FgMqUoGC0iampGQTXjMBCJ77qkJoFJKfQK+ecIsHnq1mc68yx/eadMkUKmRyngefj4PgF1Xh9/eju/72DW12Ik4dkMDIoIpbd8GYFenkFAIr7UVEKx4DEkmkGgMKxEHY/BaWzGejxWJ4KXT4HpIOIRVW4Pb1IQpFLHiMcKTJhGdtAd2Q0MZfwpKKTWy2BU439oqdwBKDTd+LgcGrGgESnP+/GIRKxrFSiSCBNAYvM7OoF0sihWPB3MG8wWAUg9ekAwCwQ4g6QxWyAbAuB5WNIJdXY2z8V2cLVsxxSKRiROJTp1CqLFRF4AopZT6WDQJVGoAjO/jZ7MASCIRLARxikg4gliCxOMA+Ok0xnGRkI1VXY3xgpIvkkhgRSP4hSJWVRWI4KczeO0dWCEb33HwcnmsSBhqa3E3b8ZpasKKRYlMnkRkj0mE6urK+SNQSilVIXQ4WKkB8LNBL6BEwqUSLj5+oYBdlULiQS+gn88H7QTs6mpMPo/b2oqfL+AX8ghgJZJIOITJ5nDb2wGD77rBkG/Ixgf8piac1jasVBXRadOITp2KnUyW+SeglFKqUmgSqNQuMsZgckEvoJVM4nd1Bb2AoTAIWIk4xnWDYWDALvX0eV1dmGwWv+iA7/XMA/Tz+dIK42LwfMftSQC9d99FIlHis/Yj1NBAqLEROxYr11tXSqkRLxYWJA9ja6PlDmXQaBKo1C4yuVxpn94QGINxPUyhiCSTwUpgkffnAUYjSDweLPYoFvEcF1MsIpEwkkwEdQW7uoKhZQHfcbHDIXzXpfjeFrz2dqJ7TCR+yMFYpSFmpZRS5TN7Yh2j2jPUJRLlDmXQaBKo1C4wxvTMBbSSyaBEjOeCbSOWYCUS788DtC2s6urg2PXwcrlgTqDvYyerEMsKVhSn0/i+j3geEgnjFQrBDiDpNKGGBmIzZiDRyvmLUymlhrP32rJs6SwQs8sdyeDRhSFK7QKTz2M8H+leveu4mHweicWw4vGgZ69UD9Curu45Nk4REQtyWaxkEsIh/I5OvHQar1AIdhEJ2fhdaQob3+1JABNzDiIycSJi6f+iSim1O3i3Pcf2jMOa5nS5Qxk02hOo1EcwxuBnMgBILNYzFxCxghXBkQheRwcAdlUyGBZub8cYH4wJksRkEisWDXYWSafx0hkkGsV4Ll5nGmdrE36xQGT0aOKfOpDoHnsgdgX9uamUUsNcaefPiqLdDEp9BD+dDnoBI2FwXXzPwxQKQU3AWBy/q+v9eoDJYGcPTDCH0MvmIGQT3XtvALzOLvyudFDnLxHcW9i0GVPIEx0/nsTBBxOdMlkTQKWU2s10V2YNVdAIjfYEKvUheg/zSigUDPHm8kgsDrYFrhskiOEQVnU1AFYigVtsx+tK43d1QsgGY/ALBUyxCIkEYgl2qhZn02YsEcJTppI48FNExowp59tVSin1AeoTIXwxpCpoUqAmgUp9CC8dzP2QaKQ0LzDYL1hCIRDBOE6Q0NXUBEPGntezSMTr7ETiMdzOTtym5qC9beF3duC0bCM8ajThiRMITxhPdMoUwloEWimldlsu4LnljmJwlbVPU0TOEJE/iUhaRDaLyL+LyPh+bSwRuUJENolIVkR+KyL77eRZs0TkmVKbTSKyVHRfLfUx+LkcpuggtgW+wfgGU3SwYnHAgOsFBaFravBzefxMFr9QwO3oAKQ0fOxBLhckkAJeayvOe1sorluPs2kT4VGjSMyYoQmgUkrt5jqyLlnHoyvvlTuUQVO2JFBETgZ+DbwAnAx8G1gAPCQiveP6LnA5cG2pnQs8KSJ1vZ7VADxVunYycE3pvm8P/TtRlcj4Pn6pFxDbxjhOsNI3GgXPxZRmCNupVNC2tHDEGIPX2oYkE4TGjw/aFopIQz2mqwunZTteZyfh+noiU6cSqqvTOoBKKTWMuL5f7hAGTTmHg88E/miM+Wr3CRHpAO4H9gFWi0gcuAS4xhhzQ6nN74GNwHkEiSHA+UAEONUY0wk8UUoMLxORnxljCp/Um1KVwU+nMb4BKxjyNab0P70xGM9D7BBWMoGEQrhtbcG1kI2zfgNuVyehjqCAtLFDmIjB2bAh6FXEEJ04geje04juORW7goqOKqVUJeteHSyhylkYUs53EgY6+53rPu4exj0MSAF3dzcoJXmPAif2uu8E4LHStW53AzWlZyi1y0yxiJ/LYzDB//UG8DwIhfHzOSQSRSJhrHg8KA1T2kvYbWrC7ezEb23F7erCbWvH5PM4GzbgbNiI37KNyMSJxGbNJrb3XsG2ckoppYaFaMQiErKojVTOwpByJoG3AfNF5B9EpFpEpgM/AH5rjHm91GYfwAPW9Lv3DWBGr+N9gNX92qwp3TsDpQagZzEIEiSAvgd2CD/dhRWLISEbq7o62Pe3VDrG6eig2NSM29yEVZVCkGALuHffxd3ahIRDhCdPJrrnXkQnjA+2mVNKKTVsJCJhUrHKWk9btiTQGPMQ8CXgFqADeJOgd/DUXs3qgLQxpv8szPbStd7tOvo93wPS/dr1EJGviMgqEVm1bdu2j/FOVCXxczmM4xJkfyW2jSkUgiHgcBi7pgby+Z4t4ozv42x8F3fLFqxkVTBvMBaFXDZYLVxXS2z6dJKzZhMZMzrYOUQppdSw0pEp0JYpEtXh4I9PRI4EbgL+GTgSOAOoBVaISO++1p3V6JadnN/VdkFjY24xxsw1xswdNWrUQMNXFap7f2BKxUDFtjCuh8llkUQcK5UKysOEw8Hq32SSwtq1FN97D7EtQtXVGMDZuBEnkyUyaz+q5swhvs8M7Lo6TQCVUmqYKvoG30DaqZytQ8rZr/kT4H5jTM8KXhH5C8FQ78nAfUAbkBIRu19vYA1Bb2C3NoIEskcpkUz2a6fUB/JzOUyp7Evv/YH8TAaJxrCjUSQUwkunsRIJQnV1ZFevprDxXchmsadOxTgOhXfewdu+nfCUyYRTKUI1tVjJRLClnFJKqWHJ94PPhaQdKXMkg6ecfZozgL/0PmGMeRPIAXuVTr0J2MDeO7n3jV7Hb7Lj3L9ppXvfQKld0NMLCMFij5CNl81hnCJWIo4kk3jt7fiZLKZYpLBlC/nXV+O1bic8cQLGdci//TZeczP22LFEJ0/Grq7BSsR1EYhSSqndTjmTwA3AQb1PiMi+QBxYXzr1O6ALWNyrTQo4Hnik162PAseXrnU7jWCe4AuDHbiqPH4+H/QCWu/XFzfG4KfTWLEYdnU1JpPp2SLO8zxyf/wj3rZthEaPxiAU31mH29xMaPx44vtMJzplKnYijp1KfcgrK6WUGg48UznDwN3KORx8E/AzEXmPIKEbA3yPIAF8GMAYkxeRHwOXi0gr8BZwKZAHbu73rIuAe0XkRwS9gJcDV2qNQLUruos9A6VewFL9P2OChSC+j18oIpYgySS5l17CeW8LVqoKSSQovP02XksL4cmTSEzfh/CUKdjxIHlUSilVOepryh3B4ClnEng9UAQuICj23A48D3zHGNPrE5mrCGZpLSVY6ft74BhjTFt3A2PMdhE5BvgX4AGCOYLXAD/8BN6HGuZ6egEFKM358J0iplDAikWxEgm8zi4AJJUi/+qrFN5+J1g8UlVFYe3b+Nu3E540icS++xLeYxJ2LIqlCaBSSlWMCtoopEfZkkAT7Lv1i9LXh7XzgStLXx/W7hXgs4MVnxo5enoBRYIFISEbr3k7AHZ9fc/2cVYygdfWRn71anzXITxxD/zOTvyODkKT9iA+exaRiXsEhaRratCtq5VSqvI0hssdweCprKqHSg3QTnsBM1nwDVYi2NPX+AaJhBHLws1kIBzGqquHqirC9XWEJ0wg1NhAuL4BiYSxa2s1AVRKqQoTTBkXEhVU7F+TQDWida8IFsvCeD4Gg8mkQYTQmDGYTAaxLXzHIfs/r2DyeSQWo/jWGtx166g65mii06aB6yHhELb2ACqlVEXae2wS46BJoFKVwPh+sDuIvF9R3BSLGN9gp1LYiQQmFsPP5Sj8+c84mzbhey7uthb8YpHouHHB9m+uh4TsoAfQqpxK8koppQLZbJZcwSdsV86+waBJoBrBjOO8f+D5GM8FxwERTDiEu20bVl0dxc2b8ba34mYyGNfF5LKEa2qJHbA/oVQ1CEEPoCaASilVkbLZLBtbs4Qq7Pe8JoFqxDLFYvDdN4gIxvMwxSISDiOeh7FsvI4OJBLFc4r4jgPFIuEpU4hNn064oREJhbCrq4Ot5JRSSlWkFidYN2h2vhPtsKWfXGrE6ukJ9D2M74MxGN+A+Eg4gink8bZnMASJord1K3Z1ivje0wjV12PF41jJBFYsVtb3oZRSamh1lDaUskV7ApUa9rrnAxqnCOEwfi6HZdtgWVghGz+fw23vwGluorB+Pe6WrYhlEZo0Cau2BitZFawE1u3gdqqzs5Pm5mac3kPuakQKh8OMHj2aaq2bqYYxsQyJsE11rILqw6BJoBqhunsBje8jjoNYNl6xGAwJR8L4nZ14mQx+axvu1iZMPk94yhQSM2ZgJ5JYpZXAakednZ00NTUxYcIE4vG4rpYewYwx5HI5Nm/eDKCJoBq28o5LJAShyloXUta9g5UqG+M4GGOC7eBKtf9M0QHj43d2YlwPt7kJt1jETiSI7jeT5MFzCVXXYMWiwTzACpsgPFiam5uZMGECiURCE8ARTiSoqTZhwgSam5vLHY5S/2sdaYf2nEe66JU7lEGlPYFqRDKOE5SD8TxAwAkWifjpNBKLU3xvM4V16/GbmghPn0Zs6p5Exo9HEOyqKiQSKe8b2I05jkM8Hi93GGo3Eo/HdWqAGtbacwUApLLWhWhPoBp5jDFBElgoYBwHiUbBBAmgcT289jaKGzfiNm3FRMLYiQRWdQpBsOKxoDag+lDaA6h6038ParjLFoIeQNuqrH/LmgSqkcd1wYApFsDzsWJRJBbFuC5WMkmxZTvetm2IZRPeYxLRyVMIJauQcAgrlSp39EoppT5hdliwLSEWrqxJgZoEqpHHBP35xi3N7bBtnK1NGN/HS3fht7VhHAe7oYHI+PFY1SmsRAI7ldIejRFi2bJlNDY2ljsMpdRuIh4KU5sIM6oqWu5QBpUmgWrkKSVyxvcB8Lu6cFtaMI6D19aO09ICkQih+nrs6hR2sgorHkPClVUaQCml1K5pTRdpzzoUXF0YolRFME4RI+BiwHXxC/kgASwUtf/tmgAAIABJREFUsGtrg4LQsRgSj2Elk+UOVymlVJls68rj+Yac65Y7lEGlPYFq5BF5vzxMoRCUiPF9vI7O0jZxkaAgdCJBKJXCTiaRCts0XH0869atY9GiRVRXV5NKpVi0aBHr1q3r00ZE+Od//mcuu+wyGhsbGT16NF/72tcolrYr7PbMM8+w//77E4vFOOigg/jd735HY2Mjy5YtA+C2226juroaz3u/B2L06NGMGTOm59jzPFKpFLfffjsAq1ev5vTTT2fChAkkk0kOOOAA7rzzzj7xiwgPP/xwn1g8z2Ps2LEsX76859yrr77KSSedRCqVIpVKsXjxYrZu3frxfoBKDTP5Ug9gpMJKg1XWu1FqFxnjY3yD35UGx4FQCCxBigWs+jrsaBSJRrGqq3U1sOqjUChw9NFH8+abb3Lrrbdy2223sXr1ao4++ugdEryf/OQnbN26lf/4j//gkksu4Re/+AXXX399z/XNmzfzuc99jtGjR3PPPfdw4YUX8sUvfpFcLtfTZsGCBXR1dfHnP/8ZgDfeeIO2tjZaW1t54403APjzn/9MOp1m/vz5ALz33nvMmjWLW265hQceeIBTTjmFf/iHf+C+++4DYOrUqRxyyCHcddddfeJ99tlnaWpqYsmSJQCsXbuWww8/nGKxyK9+9Stuv/12XnnlFRYuXBj8IaXUCFFw/HKHMCR0OFiNPCLguph8Ptg2LmRjYjHC48djV9eA5yHhMKGqJKGaGi0Krfq4/fbb2bhxI2vWrGHq1KkAzJkzh+nTp3P77bdz3nnn9bSdMmUKd9xxBwDHH388zz//PPfccw8XX3wxANdddx2JRIIHHnigp7ZidXV1TxIGsPfeezNu3DhWrlzJ3LlzWblyJQceeCDGGFauXMmMGTNYuXIlY8eOZe+99wbg6KOP5uijjwaCkkgLFixgzZo13HbbbZx66qkALFmyhOXLl1MoFIhGg8nud911FwcccAAzZswAYPny5YwfP56HH36YcGlO7OzZs5k5cyYPP/wwJ5100pD8jJXaXUXsyvo8qKx3o9Su8jz8bAZsGwMU162jsPZtRMCuqyU8YQJ2QwOivYCDwmtvx2lqLsuX194+qO/l5ZdfZs6cOT0JIMCee+7JQQcdxMsvv9yn7bHHHtvneObMmWzatKnn+A9/+APHHntsn+Lan//853d4zXnz5vHcc88BsHLlShYsWMCCBQv6nOvuBQTIZDJceumlTJkyhXA4TDgc5s4772Tt2rU9bU4//XS6urp49NFHAXBdl/vuu69PAvrkk0+yaNEiRATXdXFdl7322ovJkyezatWqXf+hKTXMFUvTMSLhyqoQoUmgGpH8YhG3vR1ne2uwMjiXw/g+vufhptOIMdi1tVoSRu1gy5YtjB49eofzY8aM4b333utzrq6urs9xNBoln8/3HG/dupVRo0b1aROLxaiqqupzbsGCBTz//PM9vX/z589n/vz5rFy5EoDnn3++TxLYPfR80UUX8fjjj/OHP/yBM888s89rT5w4kcMOO6xnSPipp56ipaWlTxLY0tLCtdde25NIdn+tW7eOd999d5d+XkpVgqITTH+IVNj8cB0OViOPCF5XGlMs4he6sOIJvGKRSCqFWBZWIgkhG1s3ux80dm0tlfKrc9y4cbz22ms7nG9qamL27NkDetbYsWPZtm1bn3P5fJ50Ot3n3Pz582lpaeGJJ55gw4YNzJs3D2MMGzdu5PHHH2fbtm19ksAVK1bwta99jW9961s959ydrGpcsmQJl19+OblcjrvuuouDDz6YPffcs+d6fX09p556Kl/+8pd3uFfrKKqRpOhV5pxA7QlUI48IxveQUAivoxM/3YUdi2GKRdzWNvA9QrW15Y5S7aYOPfRQ/vjHP/ZZDbxu3Tr+9Kc/ccghhwzoWQcffDBPPPFEn4Ug999//w7tZs+eTW1tLVdddRUzZ86koaGBxsZG9t13X6666iqqq6vZf//9e9rncrmeOXwAHR0dPPbYYzs8d/HixWSzWVasWMGKFSs444wz+lw/+uijee2115gzZw5z587t8zVlypQBvVelhrOuYjAcXBOvrHqx2hOoRhwRIT5lKt727eTfWgOugz1pEn4mgykW8do78AGTzyO95mqpkaVYLHLPPffscH7RokX88Ic/5KSTTuLKK6/EGMP3vvc9Jk2axFlnnTWg1/jGN77BjTfeyMKFC/nmN7/J1q1bueaaa0gkEli9FiRZlsXhhx/OQw89xPnnn99zfv78+dx8882ceOKJfdofddRR3HjjjUybNo1UKsU111xDfX39Dit6x44dy4IFC7j44ovp6Ojg9NNP73N92bJlHHLIISxcuJCzzjqLhoYGNm/ezBNPPMGXvvQlPvvZzw7o/So1HK1a10J1IgQmykULJpY7nEGlPYFqRLIb6knMnk1o7FhMsYjJ5vAdB6+tDeN5eNu24ba1lTtMVUZdXV0sXrx4h681a9bw5JNPMm3aNM466yzOPvtspk+fzlNPPUUkEhnQa0yYMIGHHnqI5uZmTj31VK6//npuvfVWPM+jut90hO7h3t7Dvt3/PW/evD5tb7jhBg499FDOP/98vvrVr3L66afzhS98YacxLFmyhC1btnD44YczcWLfD7jp06fz0ksvEYvFOPfccznxxBNZtmwZsVisZyWyUpXu5mfW0pVzOWn/CRU3DUK01hPMnTvX6Eq3kcdtbaX9sccprF1LeM+peJs2YTyP6JQphMaMJdTQQGT8OEKjRukCkQFYvXo1++67b7nDGLa6F3k8/fTTHHnkkeUOZ9Dovws1HD3x2nt8d8Wr2JbwwP89nMaa3aZixKB8KJW1J1BEFonI/4hIQUTWici3dtLGEpErRGSTiGRF5Lcist9O2s0SkWdKbTaJyFLRT271YUpDvVYsSqSuHqu+HlMo4hUKABTXB2Vj3JYWTIVtFaR2H9/5znf4r//6L5599lluvvlmzjjjDPbff3+OOOKIcoem1Ij3ryvfBmDvxsTulAAOmrLNCRSRw4H7gNuAi4FDgR+KiG+Mua5X0+8ClwOXAG8ClwJPishMY0xb6VkNwFPAK8DJwHTgnwAXuPaTeUdquAnF44QaGvA7OpC6Oqy2Nqx4DCsSDeYHuh7u9u0Y1yU8flxQO1DnCKpBlsvl+Na3vkVTUxOpVIrjjjuOn/70p33m+CmlPnmvbOzg3bagrNK5C/YqczRDo5wLQ74PPG+MOad0/LiI1AHfE5GfG2OKIhInSP6uMcbcACAivwc2AufxfoJ3PhABTjXGdAJPlBLDy0TkZ8aYwif4vtQwEps6BSsUQowhXFODn21AQiGsaBSrtgY3k8U4RZwtW/HaOwhPnEBo9Gj9gFaD5rrrruO666776IZKqU/Ur158h6qIxZ6NVSyYMbbc4QyJcn6SHQA82e/c40Ad8JnS8WFACri7u0EpyXsUOLHXfScAj5WudbsbqCk9Q6mdcrZtw2luxs9kIB4jOnkKocYGJBJBolFCVUmsZBVWIo6fz1F8Zx2FN9/Ey2TKHbpSSqkh8saWFl5Yt52CC2ccXFkrgnsrZxIYA4r9znX32HXPHt4H8IA1/dq9AczodbwPsLpfmzWle2eg1E74hQJ+Zxd+JoPvuYjjIpYQamggNKoxKA9jh0AEO1lFeNw4sC3clu3kXnuNwubNmGL/f8JKKaWGu39+Yi1F1ycVs/nszAnlDmfIlDMJXAvM7Xfu4NL3+tL3OiBtjPH6tWsvXaNXu47eDUr3pPu16yEiXxGRVSKyqn/FfjVCeB5YApEIlmXhuw4ml8c4TjD3r1jECPj5HH6hgBWLEdlrL0L1dZiiQ3HdevJvv43X3o7xK7OavFJKjTQtLS38dXMXAIsPqtwEEMqbBN4ELBKRc0WkTkSOJ1ggAkEPXred1bCRnZzf1XZBY2NuMcbMNcbM7b93pxoZfMcJVv2KIJYFoTDG+IgxSDQabCFn2aUeQRtcFwuI7bMP0SlTsGIxvPYOvGwWU9Bpp0opVQmeXtuO5xtq4iG+OH9aucMZUuVMAm8jSAR/AbQCK4B/LF1rKn1vA1Ii0n/b0RqC3kB6teuzz1fpnmS/dkr18NJpTDYHxmDicfxsBhELqakB14VwkBRaloUVjUI4HMwVDIWQWBSJxbBiUcQOIeHK2kpIKaVGKjsUpjYeYmJt5ZWE6a9sSaAxxjPGfBUYBewPjAFeKF1+qfT9TcAG+pemn0EwL5Be7frP/ZtWuvcNlNoJt60NL5PBCofAN4jvQyyKnUjg5fPBfD/fYCWSiG0FyV4shjEmGAJ2HOxUComEkZDuwKiUUpXAcx0yRZ9c0Sl3KEOu7HUujDFtxphXjDFdwIXAC8aY7sTtd0AXsLi7vYikgOOBR3o95lHg+NK1bqcRzBN8AaX68YtF/O2tGMfBqkphMkEvoF1dHQwNOw6+E/wCkFgs6AG0BAmH8bNZTKEInotUVQW9hEoppSpCW9al6Pp05vsvR6g8ZUsCReTTInKxiBwjIqeKyH8RJG4XdLcxxuSBHwOXi8iFInIMcC+QB27u9bibAAe4t/S8CwgKTF+rNQLVzvjpNG5XF9g2ViSM19WFRMLY9fWYQiGYK+i6WNFoMPwbDkrGiEjP/sJWPI5YFhKLlfvtqCFwxx13MGfOHFKpFHV1dRx44IF861vvb2q0fv16RIQHH3xwyGIoFossW7aMv/zlL0P2GkqpkaucPYEOsAT4b+AOIAEcboz5n37triIoCr0UeICgKPQx3buFABhjtgPHlK49AFwBXAP8cGjfghqufMfBT6eRUAg/HEZcF4nFsRMJ/Hwev3soOJVCQsGUVInF8F0Xv6sr6EFMpbAikaDnUFWUq6++mnPOOYfjjz+e++67j3//93/n5JNP5v777/9E4ygWiyxfvlyTQKXKwPN3uq60opRtIpMx5o+8XxLmw9r5wJWlrw9r9wrw2UEJTlU8O5Uiuuee4BTxEfxEnPDYMT3LyU2hgIQj2IkE2HbPULDX1obxDdgWVjQWrBxWFefGG2/kggsu4Oqrr+45t3DhQr7//e9/YjEUdMW5UmWRiPZfi1q5tAtDjTi+71NYvx6/UECSVdi1tURmzCDU2AihEKFRo7Hr67Hr6rAiYQTpGQp229owroudrAoSw0ik3G9HDYH29nZ2VjpKRHY4l8lkOPfcc6mpqWHixIksX74cY/r2IDz99NMceuihxGIxxo4dy0UXXUQul+u5/swzzyAiPPbYYyxcuJBkMsm3v/1tUqlgmvNZZ52FiCAirF+/ngULFnDhhRf23P/QQw8hIlx66aU95+666y6i0Sj5fLD36b/9279x2GGHUVtby6hRo/ibv/kb1q5d29P++9//PmPHjsXvV/PywQcf7Hndbv/6r//KfvvtRzQaZfLkyfzoRz/alR+rUmo3o0mgGnG87a3kX32Vwvr1OO1tOOvWBXMALQvjOGAJdm1dUBS69KEvsViwICRfwLgOkqrqSQxV5ZkzZw7XX389d9xxB9u3b//QtpdeeinV1dXcc889/N3f/R3Lli3j3nvv7bn++uuvc8IJJzB27Fjuu+8+rrjiCm6//XbOOeecHZ51zjnnMGfOHO6//37OPPNMnn76aQCWLl3Kiy++yIsvvsi4ceNYsGABzz33XM99K1euJBaL7XDu4IMPJlaas7px40bOPvts7rvvPm699VaKxSILFiwgU9oC8YwzzqCpqYlnn322T0x33303hx56KFOmTAHgxz/+MRdeeCGnnHIKDz74IOeffz5Lly7lhhtuGMBPWCm1O9C6FmrEcdrb8NrakUQCv60d8jnwfcgX8Nra8DNprLp67HFjMUWnZyjYKSUDViyGJRaWLgjZZe3ZIgW3PLuqREMWtYmB9dj+/Oc/Z9GiRT09cPvuuy9f+MIXuPjii6muru7T9ogjjuAnP/kJAMceeyyPPPII99xzD6eddhoAV155JVOmTGHFihVYpfmjVVVVnHXWWSxdupR9992351lLlixh2bJlPcfpdBqAvfbai09/+tM95+fPn8/VV19Na2sr9fX1rFy5krPPPptf/vKXZLNZEokEK1eu5HOf+1zPPVdccUXPf3uex5FHHklDQwOPPPIIp512Gvvuuy+zZ8/mrrvu4sgjjwSCIenf/OY3LF++HIDOzk6WL1/O9773PZYuXdrznjOZDP/4j//IBRdcgG2PnKE0pYY77QlUI4rvebib38P4PiYaAddFokFtQL+Qx89m8PKFoIev9IEt0SjG94PagMVisFjEtnQouILNmjWL1157jQcffJCLLroIy7L4wQ9+wNy5c3sSs27HHntsn+OZM2eyadOmnuOXX36ZU045pScBBFi8eDHGGFatWtXn3hNPPHGX4jvssMOwLIuVK1eSz+dZtWoV559/PtXV1bz44ou0t7fz6quvMn/+/J57/vrXv/K5z32OxsZGQqEQ1dXVOI7TZ0h4yZIl3HvvvbiuC8AjjzxCV1cXixcHVbpefPFFMpkMp512Gq7r9nwdc8wxNDU19XnfSg1XVdGg+L+/8w3HKor2BKoRxc9k8NNprOpq7FQK8gVCDfVYiQTOtm146QwSj2M3NkCxCAQ9f25rG/gGLMGKxrB0QciADLQnbncQjUY56aSTOOmkkwC47bbb+PKXv8ytt97K17/+9Z52dXV1O9zXPQ8PYMuWLYwePbpPm0QiQVVVFe+9916f8/3bfZBUKsWnPvUpVq5cSW1tLclkklmzZjFv3ryexBCCZBGgo6ODE044gcmTJ3PjjTeyxx57EIlEOPLII/vEumTJEpYuXcrTTz/Ncccdx1133cX8+fOZMCHYP7WlpQWgT+9lb++++y6TJ0/epfeg1O5q5h611MZDjKup/N/zmgSqEcVtacFrbYVkAmIx7GSSyLRpiGVjVVVhJbsIVddgx+P4mSwSsjG2jbu9BeO6WKkUCLoqeAQ6++yzueSSS3jjjYFtQjRu3Diam5v7nMtms6TTacaPH9/n/EDmmM6fP5/nnnuO2tpa5s2bh4gwf/58HnroIfL5PLNnz6a2NthN86WXXmLr1q08++yzTJ8+vU8Mve29994cdNBB3HXXXcybN48HHnigz6KP+vp6IFgsMmbMmB1i2meffXY5fqV2V3E7wqS6GNXJyp/yo0mgGjGMMRTXr8fLZJBwGNOyHVNfB76Pn8thCgVCtXWEx44JFogAVjyO19YOng8Y7GRVT5FoVbmam5t36JXbtm0bHR0dO01+Psyhhx7KihUruOaaa3qGhO+55x5EhLlz537ovZHSlIPevXXd5s+fzw033EA4HOaUU04BYMGCBVxxxRV0dXX1GQruXokc7rXH9d13373T11yyZAnXXnstRx11FPl8vmduI8BnPvMZ4vE4W7du7ekhVarSJMLwelOGiJUtdyhDTpNANWK4nZ04TU0Y20aMwRQKweKOQgGnpQWvvZ3wuPFYySReRycIEI3ibdoc7BCSTAJgJSp/U/GRbvbs2Zx88skcd9xxjB49mg0bNvBP//RPJBIJvvjFLw7oWUuXLuXAAw/klFNO4fzzz2fdunVcdtllnHnmmR84rNotEokwdepU7r77bmbNmkUsFmP//fcnEokwf/58XNflhRde4Kc//SkABx54IJZlsWrVKi655JKe53z6058mkUhw3nnn9fRmXnfddSR28m/59NNP59vf/jaXXHIJRx55ZJ9kuLa2lmXLlnHRRRexYcMG5s2bh+/7vPXWW/z2t79lxYoVA/rZKLVbcrK4nsH1Kn9OoHZnqBHD3b4dP5PFlFYvSjhEqLERUyzitbcjoTB2bU3fXsCOTozjYDBYVSmsWBTR1Y8V73vf+x7r16/noosu4rjjjuOKK65gv/324+WXX2bq1KkDetZ+++3HI488wpYtW1i0aBHLli3ji1/8Ir/85S936f6bbrqJlpYWjjnmGA4++OCeeYSjRo1ixowZJBIJDjroIABCoRCf+cxnAJg3b17PM8aOHcuvf/1rNm7cyOc//3nuvPNO7r77buI7mdYwZcoUDj30ULZs2cIZZ5yxw/VLL72Um2++mYceeoiTTz6ZM888kzvvvJMFCxYM6Oei1O6qsbERqzQ1o3sebKWS/kVNR6K5c+ea/qv0VOXJvvUWhTffwi0WsApF7Po6YtOm4aXTuJs3Y9c3ED/oQPz2djBg19bgvPsuXi6PHYtiVdcQqq9Deg2pqR2tXr36I3u41Mij/y7UcHLkj5+i6Pr8+txDmNRYU+5wdmZQitRqT6AaEbx8HmfTJqxEglAiiRWPE2poxLguJpvFHjWayKQ9ENcFAxIJ4+fy+IUiGB+pSiGRsCaASik1AtQmIlTHw2ztKpY7lCGlSaAaEQpvv03x7bcpbNiAn8lAOISVSuGnM7gdHWDbhEaNwi9NwJdoFK91O8YYpLQQpHtOoFJKqcrmOC4dWYeOdO6jGw9jmgSqiud7HsVNm/AdF1wHk88H5WCEYC6gWISSiWCxiOMiluAXHfxcHowf7BMcsrG0OLT6/9u78yhJq/KO499fvd09vcxUL9M9zaIngwqMShQNeFwAwQSNIDIjGDUkYNRoNMZEI0piUPQkMQH3Y4yg0dEQF8wR3BAUcSQKCAORyYwzw8AwzNYos3ZP71315I/7FtRUV3dXVVdXdVU9n3PeU9W33lt179O3u59+3/fe1znXECZTMJU2Hhscr3ZTFpQnga6upSYmGNuwgfFdu5g6fJh0Og1mqLmZiYEBUsNHwh1DentJDQ9jGOmpKdL7Hic9MRHuDxxFPiPYOecaSJQIl9yNjKeq3JKF5UvEuLqUGhpiZPNmRn71AJO7djF14ACJJUsAsOFh0qkpoiWtRO3tJJa0kB4bIz04SOrwYUhE2OQkihJEK1agpgj5fYKLYmZFLXzs6ptPQHS1pmNJxL5h2D88We2mLChPAl3dMDPSY2OMPbyd4bvuZGLHo0wdOBBebG9HHR0kmpqIlvcStbaS6Ook0duHooiJh7djqRQkEiRal9Dc002ip4emZcvCNYGe0BSsubmZ0dHRvGvQucY0Ojp61ELVzi12T+9t59DIJJPp+v4HxpNAV/MsnWZ8925Gt2wlPTZKenCIiUd3kh4dJdHTTctxxxEtXYaamkgsaaGpr49o2TIslSI1NERKgihBU2eSRGcnTT09JDo6/BrAEq1YsYI9e/Zw/PHH0+YJdEMzM0ZHR9mzZ0/Rd1pxrpp6lrYwkTIGDk2/W0898STQ1azUxAQTO3cycv/9jG97iInHHyeKIhJ9vUQ9PTT19aFkkrYTTiDq6Q5H+sywkdFwRxAgWrKEaNlSou5uoq6ucEs4Xwx6XpLJJAB79+5lcrK+T6W4uTU3N9Pf3//EuHCuFrS1iNHJFLsP1vet4zwJdDVnaniYia1bGdm0iYmBx5jaty8s7RIliPr6SLS2kujoIFq6lKinm0RXJ0qlSA8OYek0AFFHO1FXV0j+li5Fra1+xKqMksmk/9F3ztWsp/Ymgb0MTfjEEOeqzsxIjYwwsmEDI/ffz+TOXaSGhkCQaO+guX8F0bIkS048kab2Nswg0dGBpaZIHThImrAAdFNnkkRXF01dXUTt7chP+TrnnMtxxlPCNc3jk+kqt2RheRLoFjVLp5nYv5/x3buZfOQRRh7YwNTAAEqnsfZ2Wvp6aepbQaKjg6ZkMtzWLYpgdIz00BAGRO1tRMlkOPKX7CTR1uqnfJ1zzs2ot7eXpkhMpYwtA/tYdWxvtZu0IDwJdItSemKCsd27GVm/nokHtzF16BBqbSU9NYk6Oog6kzT3H0Nzfz/Nfb2ouZmERGp0jHQqFZZ3SS578pRvR0dY889P+TrnnCtAe3PEYGqKzXuOeBLoXCVMjY4yvm1bmOyxfTtT+w+QHhlBLS1ECdHSs5yop4eoM0nU3U1zTw/pyUnSwyMYoCUtNHd3kejupikz0cOXpnDOOVekjpaIwbEpduw7Uu2mLBhPAl3VmRmpwUFGNm9m9L77GH90J6nDh7HxcdTcRNTTQ0t/P03LlxP1hFm8amtDE5OkDg9iQNPSjrC8S1cXUWZtv4TfEMc551xpkq1NDAyOM1DHt47zJNBVXXpoiCP33MPIveuZ3L8fRkdJJJNEHR00dXWFU77H9BN1dRO1tzF16HBY7gUj6u4Kp3w7O4mWLn3iriDOOefcfPR0tADD7B+aqHZTFoz8dj4gaQjYWu12LDK9wL5qN2IR8rjk53HJz+MyncckP49Lfh6X/FrN7JT5vokfCQy2mtlp1W7EYiJpvcdkOo9Lfh6X/Dwu03lM8vO45OdxyU/S+nK8j1805ZxzzjnXgDwJdM4555xrQJ4EBtdVuwGLkMckP49Lfh6X/Dwu03lM8vO45Odxya8scfGJIc4555xzDciPBDrnnHPONaCGSwIlXS/pIUkjkvZJulXS6QXWvUjSJkljkjZKWrPQ7a0EScdKukbSA5KOSNol6auSji+g7lpJlmd7RiXavpDmE5e4fr2Ol3dI+oGk/fH3+uwC69XtWIHS4xLXrcuxAiApIelKSbvj37s/lfTsAurVxXiRdIqkdXHfd0v6BxVw/0pJ3ZK+IulQvH1FUlcl2lwJpcRF0soZxsT1lWr3QpP0DEnXxn93UpLWFVivpPHScEkgEAEfBV4FvBloAX4i6emzVZL0UuAG4DbglfHjtySdsbDNrYjfAy4E/gu4ALgceAFwl6RlBdTfArwoZ9u1ME2tqJLjUufj5VKgB7i1hLr1OlagxLjU+VgB+ADw98C/EH6epoDbJHUXULemx4uk5cBPCH2+kPC35wPA+wuo/i3gLMLfqTcBZwLfWJiWVtY84wLwXo4eEx9agGZWy7OB84AH461QpY0XM2voDVgKjAHvmWO/HwO355TdDtxS7T6UIQZdQFNO2UmAAW+co+5aYH21+7AI41LP4yURP54Sx+LsAuvV7ViZZ1zqeay0AYPAB7PKksAw6JUfAAAMgUlEQVQh4Ip6Hy+ExOYgkMwq+2Dc/yWz1HtJPIbOyio7Ky57YbX7VcW4rIxj8Kpq92EBY5PIev7fwLoC6pQ8XhrxSGCuYWAcaJ5pB0mtQOa/9Ww3AC+TVNP3KjOzQ2Y2lVP2IDAC9FWnVdVXalwaYLykq92GxaiUuNT7WAFeDCwjq39mNgjcQjjqWe/+ELg17nPGDUAnITaz1fuNmd2RKYif/4b6iFupcal7Jf5+LXm8NGQSqKBJ0jHA1YRD0l+bpcrTCEni5pzyLXH50xakoVUk6TlAO/DrAnZ/lqRBSeOSfi7pzAVuXtUUGJeGGy9FaJixUqB6HysnAylgW075FmBVAfVrfbyczPTv7TZCTGbrf756UHjcFrtS45Lx5fh6ub2SPh7/M9XISh4vjXrbuL8GPhk/fxw438xmu84kc+3K4ZzyQzmv1wVJCeAzwMPAj+bY/X+BXxKSol7gbwnX+5xhZvcuaEMrrIi4NNR4KULDjJUi1PtY6QaOmFkqp/wQc/etHsZLNznfWzNLSTrC7P2fVi9WSNxqQalxGQf+jfD7d5BwyvMKwj9LdTOZqgQlj5eaTwIldQLHzrWfmW3J+vJrwM+BY4C/AG6WdKaZbZrrbXI/fobyqisxLhlXAy8EzjGzyTnqfzrnc38AbCJcCL7ofigrFZfM2+R+/AzlVTXPmBSswcZKsWpirEBJccnXB81Qnl2/psbLLErq/zzq1Yqi+2dmA8A7s4rWSXoMuFbSc83sgTK3sZaUNF5qPgkk/DL4cgH7PTH13Mx+C/wWQNItwEbCrKRLZ6h7MH7MnW7dGT8eYvEpOi4Akt4FvAd4g5ndVeyHmtmYpB8Cry62boVUIi61Nl5Kisl81etYKVKtjRUoLi4HgWWSopyjgZ0U2bcaGC/5HCTneyspAjqYvf8HyX/dcSfhOq9aV2pc8vk2cC1hNYdGTQJLHi81f02gma01M821zVJ/ipAEznbtzXZgkunn1lfF5dvn249yKyUukv6IcJr8fWb2zXl8vIBFOXmgQnGpqfEy35+heaqrsVKCmhorUHRcthKW5cpd228V4XqlYi3a8TKDrUz/3p5IiMls/c9XD0qP22JTalzyyYy1WhoX5VbyeKn5JHC+4gtKnw88MtM+ZjYG/Ax4bc5LFxOWdhhfuBZWhqRzgP8EPmtmH5vH+7QR1ji6r1xtq6ZS4tII46Uc6m2slKIBxsovgCGy+hevsfkK4IfFvFGNjpdbgFfkrCt6MeH6rTvnqHdM9lqR8fN+iozbIlVqXPK5OH6spXFRbqWPl9nWj6m3DTifcD3gJcDZwOuB/wFGgedl7fdGwnn0lVllLyXMXPpUXPeThFnFZ1S7X2WIyzMJP3y/IkzPf2HW9pSZ4kI41HwHYXHKlwGvA+4mrLv4/Gr3q1pxaYDxchrhF++7435/KP76tEYdK6XGpd7HSty/KwnLKr0D+APCRf0DQHe9jxdgOeF03I/ivr89jsUVOfvtANbmlN1GOBL8mnh7mLCsStX7Va24xD9TVwOr43ofIfz9vqHafSpjbNrj3xsXA3cRroPNfN1e7vFS9Q5XOLirgO8AewmzjB4lrE307Jz9pv2ijssvIsxUG4+/MWuq3acyxSXT33zbVTPFBWglXI+xK47JIcJ/HadXu0/VjEsDjJe1M8RkbaOOlVLjUu9jJe5bgrAQ8F7CH+x1c/3OrafxAvxu3OfROAZXAsrZZwfT/6j3AF8l/CN6OH7eVe3+VDMuhAM36+N4TBCWlbkKaKl2f8oYl5Wz/N1ZWe7xoriyc84555xrIA1/TaBzzjnnXCPyJNA555xzrgF5Euicc84514A8CXTOOeeca0CeBDrnnHPONSBPAp1zzjnnGpAngc65hiLpKkkmaWUV27AybkNm+2KFP3+tpIqsDyZpS1Y/d1TiM51zhfEk0Dk3b5LOjv/Iv3eWfUzS9yvUntWSrqrEZ83TjcCfAkclgZJ2SNqYu7OkpKQ74li+P8/rl0sal5RcuCYX7XJCH+vhnrfO1RVPAp1z9Wg14RZT+fwj0Ea4Y1C1bTCz683s7rl2lNRHuMPCS4C3mtm/5tltNfATMxssbzNLZ2bfM7PrCbcJc84tIk3VboBzzlWSmU0R7s1bMyQ9FfgxcALwBjO7Ic8+xxDua/32CjfPOVej/Eigc66qJJ0m6UZJ++JTmVslfUBSU85+L4ivZXtQ0oikIUm/kLQmZ791wGXx8+zr7t4Yl027JjCr7GRJ/yxpd9yWBySdl6fN7ZI+IWkgbsvdkn5/Ia61k3QS8HPgqcCr8yWAsQvjx+9m1W2VdI2kvZJGJd0j6eUzfE6h8f1MHKsT87zHsZKmJP1HSZ11zlWUHwl0zpVTu6TeQneOE6wbgYeAjwMHgBcBHwFOBV6btfsaYBVwA+FU7nJCsvdtSZeY2dfi/f6J8A/umYRr0TLuLKBJXwEmgY8BLcDfADdJOsnMdmTt9y3gPOAm4DbCEbobgUcK6XehJJ0K3Bq35Vwzm60Pa4C7zeyxrLKvE04Rfy9+n6cD356hnYXG91rgr4A3AX+X8x6XARHgSaBztcDMfPPNN9/mtQFnA1bA9v2sOq3AY8AdQFPO+7073v/srLKOPJ/bDmwFfp1Tvjb8esvb1qvi916Zp+z7gLLKT4/LP5pVdl5c9oWc982U5/3cnH1XxvteNcPrO4D9wKE4Rs+Z4/2SwDjw3qyyl8efsTZn39X52llkfO8E9ub5vj2Yu2/Wa+uAHdUeq7755tuTm58Ods6V03XAuTNsuc4F+oEvA12SejMbcHO8zxOnLs1sOPM8Ph27nJCk3A48s0wzYj9tZk+czjWze4EhIPvU5wXx4yeyK5rZzcDmMrQho4OQ3B0gJIKzOZ9wtPCmrLLV8eM12Tua2U2ExI6c8mLiex1wLPDKrDpnEeLkRwGdqxF+Otg5V07bzOy2fC9Iyi16Zvz4pVnerz+r/grCzN4LgRV59u0C5jsrdnuesgOEU6MZJwBpwinsXFt5sl/z9RDwWeBzwDpJ55jZTDNs1wAbzSy7TU+L2/lgnv03AydnFxQZ328CnwLeTDjVTPx8Avjq7N1yzi0WngQ656olkxVeDvxqhn32AihkkD8iJFifAe4FDgMp4M+AP6Y8E91Sc7Q1+/mCL7ZsZp+PJ5r8O08mgkcdFZS0hHBE7lM51adl3TO9Vmx8zWxU0vXA2+JZySPAxcB3zezxojvqnKsKTwKdc9WyLX4cnunoYZbnAM8FPmJmR63/J+ktefZfyATtEUJCdCLTT/+ePH33+TGza+NE8PM8mQgOZO1yLrCUMDEl28OE0+knAZtyXluV83Wx8YVwSvgvgUsJCWM7firYuZri1wQ656rlVuC3wBWSenJflNQmaVn8ZeYIXe4RrFMIp0JzHYlfn/a+ZZA5/fnunLacR/lOBR/FzK4D3kZI6H4q6bisl9cAO83s/pxq34kfL89p52qmJ6vFxhcz2wDcQ5gl/GZgJ+FoonOuRviRQOdcVZjZsKRLCZMZtkr6EuE6uC7CkarXEBKQdYQjbpuA90nKzFg9iZAYbQSen/P2dwPvBD4n6QeEZV9+aWblWMLlZkIC++fxJJbMEjFvBTYQjqqVnZl9IT4ieB0hETyHMGHkAsJSMLn73yrpe8BlcTJ8C2GJmEzMTsnavdj4ZlzHk7e8+7CZpefXS+dcJXkS6JyrmjhROR24AvgToA84SDiV+QlCUoWZpSSdT1i/7zLCzNmN8fPnMj1J+TrwPOD1hLUGE4Rr2+adBJqZSbqIsB7hGwjX420gJKzv4OiZxGVlZl+ME8EvAD8D3kKIWe6p4IzXESZ7XEI4bbwRuChu9xNJYAnxzfgG4fu0lDDL2zlXQ5S1GoJzzrl5kPR/QLOZ5V5zl7vfSkJCeg1wNTBmZkdK+LxPEhbE7jezmSa1LJh4UsoAcK+ZvWKGfboIBxy+AxxvZisr10Ln3Gz8mkDnnCuSpLY8ZecTjq79uIi3uhx4nOkzewu1GXhXNRLA2CVAN+EuIjO5m9DHF1ekRc65gvmRQOecK5KkjxJON/+UMDP2VMIEiUHgVDPbPUf9VuCMrKI9ZlbOhaYXlKQLgN8h3GnlN4Q7muRNRCW9iHB6GWDUzH5RkUY65+bkSaBzzhUpngl8BfAsoJOwoPTtwJU5CzbXJUk7gOOA+4C3mFnuEjTOuRrgSaBzzjnnXAPyawKdc8455xqQJ4HOOeeccw3Ik0DnnHPOuQbkSaBzzjnnXAPyJNA555xzrgF5Euicc84514D+H8iGB7S45MGfAAAAAElFTkSuQmCC\n"
},
"metadata": {}
}
]
}
],
"metadata": {
"kernelspec": {
"name": "python3",
"display_name": "Python 3",
"language": "python"
},
"language_info": {
"name": "python",
"version": "3.6.5",
"mimetype": "text/x-python",
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"pygments_lexer": "ipython3",
"nbconvert_exporter": "python",
"file_extension": ".py"
},
"toc": {
"toc_threshold": 6,
"toc_number_sections": true,
"toc_cell": false,
"toc_window_display": false
},
"gist": {
"id": "",
"data": {
"description": "RRTMG with McICA",
"public": true
}
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment