Skip to content

Instantly share code, notes, and snippets.

@naomi-henderson
Created August 7, 2018 17:53
Show Gist options
  • Save naomi-henderson/425bd1aadeb398f7eb2a7aaf67e84487 to your computer and use it in GitHub Desktop.
Save naomi-henderson/425bd1aadeb398f7eb2a7aaf67e84487 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"import xarray as xr\n",
"\n",
"\"\"\" Using climatological Clouds and Windspeed in the\n",
"surface heat flux formulas\n",
"\"\"\"\n",
"rhoa = 1.225; cE = 0.00125; L = 2.5e6\n",
"eps = 0.97; sigma = 5.67e-8\n",
"r = 0.8; ps = 1000; es0 = 6.11\n",
"delta = 1.0\n",
"f2 = 0.05\n",
"# 'a' should decrease when deep convection happens above 28 degC\n",
"# a = Ts-T0;a[a>28] = 40;a[a<=28] = 80;a = 0.01*a\n",
"a = 0.6"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"# basic parameters\n",
"T0 = 273.15\n",
"f1bar = 0.39\n",
"Ubar = 5.0\n",
"Tsbar = T0 + 25\n",
"Cbar = 0.6"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"# Find linearization of Q_LH (latent heating)\n",
"const1 = rhoa * cE * L * (1-r)\n",
"\n",
"def f_es(T):\n",
" return es0*np.exp(17.67*(T-T0)/(T-T0+243.5))\n",
"\n",
"def f_qs(T):\n",
" return 0.622*f_es(T)/ps\n",
"\n",
"def f_dqsdT(T):\n",
" return f_qs(T)*(17.67*243.5)/(T-T0+243.5)**2\n",
"\n",
"def f_QLH(T,U):\n",
" return const1 * U * f_qs(T)\n",
"\n",
"def f_dQLHdU(T,U):\n",
" return const1 * f_qs(T)\n",
"\n",
"def f_dQLHdT(T,U):\n",
" return const1 * U * f_dqsdT(T)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"# Find linearization of Q_LW (longwave)\n",
"const2 = eps*sigma\n",
"\n",
"def f_qa(T):\n",
" return r * f_qs(T)\n",
"\n",
"def f_Ta(T):\n",
" return T - delta\n",
"\n",
"def f_ebar(T):\n",
" qa = r * f_qs(T)\n",
" return r * f_qs(T) * ps / 0.622\n",
"\n",
"def f_QLW1(T,C,f):\n",
" Ta = f_Ta(T)\n",
" return const2*(1 - a*C**2)*Ta**4*(f - f2*np.sqrt(f_ebar(T)))\n",
"\n",
"def f_QLW2(T):\n",
" return 4*eps*sigma*T**3*(T - f_Ta(T))\n",
"\n",
"def f_QLW(T,f):\n",
" return f_QLW1(T,f) + f_QLW2(T)\n",
"\n",
"def f_dQLWdf(T,C):\n",
" return const2 *(1 - a*C**2) * T**4\n",
"\n",
"def f_dQLWdT(T,C,f):\n",
" ebar = f_ebar(T)\n",
" qs = f_qs(T)\n",
" dqsdT = f_dqsdT(T)\n",
" return const2 * (1 - a*C**2) * T**3*(4*f - f2*np.sqrt(ebar)*(4 + T*dqsdT/2/qs)) + 12*eps*sigma*T**2*delta\n"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAdkAAADkCAYAAADU1F7FAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXu0JVV1r78fDwUFVCISGjBAROVohAYkJm2MChrkSiP3\ncokmGFpFoiE+WhMFuQHRkGuiAcnQQdICAQVFbCC0b0ExRq6C3TTvxogICrQiUWzwAQHn/aOqzqlT\np967nvvMb4wzzt71XLVX1frVnGuuuWRmOI7jOI7TPJv1XQDHcRzHmVZcZB3HcRynJVxkHcdxHKcl\nXGQdx3EcpyVcZB3HcRynJVxkHcdxHKclXGQdx3EcpyVcZB3HcRynJVxkHacjJH1a0oPh31f6Lo/j\nOO0jz/jkOI7jOO3glqzjtIyko2MW7K8kPRr7fr+kx/ZdRsdx2sEtWcfpEEknAfuZ2WF9l8VxnPZx\nS9ZxumUf4Pq+C+E4Tje4yDpOt+wDXNd3IRzH6QYXWcfpCEnbAbtRQ2QlPUnSeZK+Jumbkl7ReAEd\nx2mcLfougOMsIvYGHgC+V2PfTwKnmdkXJG0LXC3pK2a2qdESOo7TKG7JOk537APcYLFoQ0lbSvqo\npK9I+oakHZM7SToIuNvMvgBgZg8ANwLP6KrgjuPUwy1Zx+mOtP7Y5wIPmNmLJcnSw/2XA2sSy34D\n+HkLZXQcp0HcknWc7tibhSK7DnhA0ieBl2fs93jgp9EXSU8AngJsaKOQjuM0h4+TdZwekfQ4M/uF\npO2BT5vZsnD5LmZ2V/j5DcBvmdkJkrYAzgXWmNlFvRXccZxSuMg6To9IuhDYlaDr5m/M7EuhkF5p\nZn8QbvMY4J+BpwGbA+eY2dl9ldlxnPK4yDrOwJB0ALC3mX2k77I4jjMZvfbJSnqipNWSbpW0QdLv\nSdpe0uWSvhP+f1KfZXScrjGza1xgHad5JG0l6RpJ10u6WdIp4fJSuiPpYEnflnSbpOPLnLPvwKcz\ngC+Y2TMJgkI2AMcDXzazPYEvh98dx3EcZ1IeAl5sZnsTRPsfLOl5lNAdSZsDHwZeBswAr5I0U3TC\n3kQ2jJB8AXA2gJk9bGb3A4cB54WbnQd4ZhvHcRxnYizgwfDrluGfUU53DgBuM7Pbzexh4MJwv1z6\nHCe7O/Bj4F8l7U0wlOEtwI5mtjHc5ofAgsH5juM4zuJD0p8C/ytnkzvM7G0Fx9icQG+eBnzYzK6W\nVEZ3dgZ+EPt+F/C7RWXuU2S3APYF3hRe5BkkTHQzM0mpkVmSjgWODb5tuR88OfNE++23pLAw69bd\nU7bclY7rDI86dT1G8u7PSX+Dsdz7TdZ11jUvlvupGhvvM7MdWjjwwa+Gw38jY+UnAEm3MqsNAKwy\ns1XRFzN7FNhH0hOBSyU9O36MPN2pQ58iexdwl5ldHX5fTSCyP5K0k5ltlLQTcG/azuGPtgpAWmLw\n55knWrv2ZML+7XDfk+etj6+rwtq1Jxdv5AyKunU9Rtatm/89uu+b+A2iYyefpaHRZH1nPe+L6Z4q\nz7vvbOvI2wJPzFi3GfO1IQ8zu1/SlcDBlNOduwmG20XsEi7Lpbc+WTP7IfADSVH+1QOBWwjSxx0d\nLjsauGzScyUfgjoPhdnJC/6ccdFGYzime0E6pfHfIDrmUIWmybrJusa0e2As98RiQ9IOoQWLpK2B\nlwC3Uk53vgXsKWn3cOz6K1mY7nQBfecufhNwQVjg24HXEAj/RZJeB9wJHNlj+QB/YKaBoYrANCGd\nMshnxezkTuo/S2j93hsUOwHnhf2ymwEXmdlnJH2DFN2RtAQ4y8wOMbNHJP0l8EXmksLcXHTCXkXW\nzK4D9k9ZdWDb544ahKIHYIiNhlONtixYZzyk1ZeLX7N09TIzCWZ2A7A0Zfl/kaI7ZnYPcEjs++eA\nz1U5Z9/jZHtl6DeEMzldCKzfR3NM82/hL1b5THPdT8KiFlnHcZpnyH20cbrqT19M4ryYrrUsfffJ\nDhK/UZyyjEFMnHziz3tecFPb55gGpvna6rKoRNbFc/HRZuPmDcr0MYZ+xaHgv1U5Fo272AXW8Xug\nW8biNu6aabkP27L6p41FIbJe6U5EE/eCi0c1/LdaPHhbu5CpdRd7ZTtpeIPvFNHU+NbFaOn587WQ\nqRTZab6JHcfphrbaEReixcXUiawLrOM4feIi6sSZmj7ZMeWQdfrBG7/+8H7shXh7tTiYCpEdy7Rb\nTn94A+90RVnx9HtycTB17mLHcYbLUCcRcMbD0Y+BZ2aYh5c+1G1ZyuAiOwHJN9G0wdneoDjOfFxo\nncXEVLiL69DW3Jplljnd4nXgOE5fLEqRjTe6dRtgb7jHgdeT0wduqTsRU+8ujjeyTeXa9IbbcZyI\ntPbARdaJmDqRzRPArHVV+ojqCKz3QXWPvwg5XZDXpjgOTKHI1iVp8bZ1fBfb9vCGbRz4M+AsJnoX\nWUmbA2uBu83s5ZK2Bz4J7AbcARxpZj/ttkxueY4NF1jH6Ye0tlJ6d/cFGShDCHx6C7Ah9v144Mtm\ntifw5fB7KdpuaL0hd5zJmYbnyDNYBbgxUkyvlqykXYD/AZwKvC1cfBjwwvDzecBXgXfmHWfdunta\nu+HLHrepmTucavjv7XSN33NOFfp2F38QeAewbWzZjma2Mfz8Q2DHzktFOZdx2np/s2uGvIbMX2ic\nvvB7zqlKb+5iSS8H7jWzdVnbmJkBlrH/sZLWSloLv2ipjP5A9UHR7+714jhOHSTtKulKSbdIulnS\nW8Ll+0j6pqTrQl05IOcYm0taL+kzZc7ZZ5/sMmC5pDuAC4EXSzof+JGknQDC//em7Wxmq8xsfzPb\nHx7XVZlj53eLtU9caMfNGOtvjGV2FvAI8HYzmwGeBxwnaQb4B+AUM9sHOCn8nkUyjiiX3kTWzE4w\ns13MbDfglcBXzOwoYA1wdLjZ0cBlPRXRGTj+ouOMg5Wxv+libC8eZrbRzK4NPz9AIJY7E3hMtws3\newJwT9r+sTiis8qes+8+2TTeB1wk6XXAncCRPZcnlTJ9hk49yva5ju0Bd+YzpqFyC++1uGCenrIs\ni2ib03O3GhNRPY7teZS0G7AUuBp4K/BFSR8gMD5/P2O3tDiiXAYhsmb2VYIoYszsv4ADmz1+tzfA\nmBqPxcIYGwFnoBwae7Y/vYlptFCr0uWztc0rYLvt09dtvjqI1wGOjS1eZWar4ttJ2ga4GHirmW2S\n9LfASjO7WNKRwNnAQYl9ZuOIJL2wbHmHME62VSKx61r0vEEfDv7CM1xG95wcmriXDt0ufbtCptN9\nPATi8TrhX1JgtyQQ2AvM7JJw8dFA9PlTQFrgU1YcUS6DsGSbwBtSJ4vRNeSLjKF7frScOXE9JrYi\n6pWLhPbTm2ocPSm00+NGHiKSRGClbjCz02Kr7gH+kMCj+mLgO8l9zewE4ITwOC8E/iqMI8plKkR2\nv/2WlNrOXYZOgDdsTjm0PPblmMTKY6gQ/lKWIut2uPdq/EVpwGkVlwGvBm6UdF247F3A64EzJG0B\n/IrQ3SxpCXCWmR1S94RTIbLDYiVDfhAWB1XdcF5nfTPECTS0nHnCutfy9bOfN6xZunCHLNdxLQs3\ni7R7u/97d0j1loeZfR1Qxur9Ura/B1ggsPE4oiIWnci2Z82unPfZrG5fjRNRva5idZBs8AobumE2\nXouNobiO4wIbF9daVOm3bcTlDF3eu0OoryGz6ES2CuUbeQ9gaItKdZDXmFUW3fCYC3DhnXa0hnkC\newSrZ9et5oh2T15GkCvdu+3ery6wxSxKkZ3fd5DfgNfJk+tWbJeEjUnUOKUFpqRx6HYNWg1lcHEe\nG0mBTSXZT9t4H20KlV4Yp29c7thYlCIbp6yl5EFTQyRhvSYbvCJqD7/IoFRjV4bF3SB24TLWmvnf\nLRbgFLmHCwU2jar3YBFlRDt+H2feg82LrVux5Vj0IgvpApo1w06Z7epwEifOfn4PpzZyzKnn0O3y\nG7VWoj9zaKzvrY/uh+kU9qSYprHX8vXMMBfMtNfyOYE9nEsBuJTDF+yTGvzUNFUt5cLhRG7Zdo2C\niW7Gzf77729r167tuxgTEYls9FAv5ZY+izMY0r0HK+cJbDIwpZPGrwxdCjw0HMU6jEa47ktskbjm\nBTMlBRbmi2xb/bKV79sy91fuPTFZHefVjaR1weQtzSLpYzcfyVEzGRmf9lsN635sWdHDveCW7ABZ\netWGYDSXkyDW/3pMA1GfbVPVdTipKNd1f6c2xOO2eMrcG0l3cFw84wIbfY+E9ghWtyK0aWXOFd4y\n8Qe5lm39oWvuKi6Pi2yDTOrynfdgXxW+jC0bv6ehGVbOC27Ka0SLXHl9iXOhpTJJf94kAp3bpzc+\nsZ2h2Eqtug4WCm0RTQhx4djciOjeyRPbhl6mXGCr4SLbEuuZqeTyTb45z/LGhOfjzMUjunOu4pWp\n1mvU0KU1Zk0JaVpjWrfxLFumWu7uKgJdFHUNGWLbvdBKp7BXLCrpFqr9NnlimPbMRSK69KpgutD1\ny/ZasE+yf7bOuZOUuaei+6e22DZg1brAVsf7ZBsk2a8K1fpW1zMz+3ADkJZ6OpYpM9kAVD3f0JHC\nxiAhsE0IX63I0QZo0s3YWN9zluimNsYdC+2hJ7PXmrl6L+MhiizZqgIbcSmHcziXsvSqDanPWHLb\npil7jxTWf97LVGZfbX79lhVZ75Odwy3Zhsl7eLugqgU9VOYJLMVjFpsWzax6nLRRLSpnFREusoxL\ni3CW9ZPqYuzHfRz9bidxYqHQpv3GWfUZf6mNBHXeS3KB0JZ53qveM3kemjiFlm1edH2NceJuxdbD\nLdmGCQYDxB7eCn2qCyxZmG/NJuZ7yHv4xy60s4nZj5kT2GSDVlfwun4RatLamcQSLiW6g7Jqg26C\nvdbMf8HKE9l4XEREXn0veN4yKLJoq1LlnihT57Wt2pL1WlVgW7VkP8ZRM7unr9/vtbDu28OyZF1k\nGyZVKKsEL12VuD9yRBbmP/zJfqSxCm18+EWWwEakNVZ9exPK0oT41hXd3Ea5UoMM7YhtaDFftt0C\nL8bhXJp5b8dFtug+KCuwWTQlvGXvg4nFtlS9pteli2x93F3cNucD56t8wNIymy+00X5J8SV4yNMa\nijkrumJZB8AM69krln0nEti8QJS6TNrIxqnT4OYF3pSlbv90rqsxy82YG6HapNAmUmXGKKrvWvdD\nwYtsFln3T9V7ISvhRZIybuTcyPquk7M4QI8iK2lX4KPAjoABq8zsDEnbA58EdgPuAI40s5/2Vc6q\nzD54aUFLZUmzfJfNiW3ph/gqjWIIUNqQi3hCgHhjViYYJU6TQlr1PEMR3lqNcl4/LeT01UI9wU1k\nuSoY89ta7MH5VBLaNNL6ectQRWxrC20aBf2z3hc7GX1aso8AbzezayVtC6yTdDmwAviymb1P0vHA\n8cA7eyxndSYR2AIqN9yTCG1y+FDEBMOI0vrM0pqL2cxXJUSyKyGtSlvCO4noZg13qmz9FA4HqUmF\nXNStCm2cCUQ3ywOTR5lhQrWF1q3ZzhlMn6yky4APhX8vNLONknYCvmpmz8jbd0h9sgCbthbbrQg/\nnzu3fLtfTi5OyQa3ksBUEdssgY1TUWzTBDZJMrKzdaq+EE1o5cTpuk8vIq9xrtWnF1EnrWOaxZoy\nUXpav3yQGW3+PRgFHmaRG1hYhZr3QZU6L1OvRV0DmfWZrMucftk6lqz3yc4xiD5ZSbsBS4GrgR3N\nbGO46ocE7uTRERfXtigtQlFD0nQf7RvL9zWvZ4bDKR+o1JjANu1VKBi7XIXkNdYV3SaHGxVatVCc\nwq8OCas1LelInKz7I27ZxgW3MXFN7l+x7qtYtk1YtJkkLdraUz86RfRuyUraBvh34FQzu0TS/Wb2\nxNj6n5rZk1L2OxY4FuCpT33qfnfeeWdnZS5i09YLX6RO/1Xw/+Qav3dWxGSmEGU1IGUtzzJWbIVj\nFlkXSVKvq0xwSotu+tI0ZOm2PWSkkWEhEXXcjxmu4CxxzbzvzyfzHmxVZONMUOdl6nlSi7aSdyJF\naOvMj+2W7By9WrKStgQuBi4ws0vCxT+StFPMXXxv2r5mtgpYBYG7uJMCl2S7X9o8oY0Etg5ZLtbK\nAluWKgLbIJUs1yGIaRY1LZwkTVm5WZSNVIWWcy6zMKlG0nLNFVjmXmqT3TFLuSV9SN2AKBPIN6lF\nW6nP3S3axukzuljA2cAGMzsttmoNcDTwvvD/ZT0Ub2IioU0K7KatVapvdu4tfO7han38Zw8Cm9sA\nDllMi2hIbCOqRlWXpYrYRkySzrHMFHMRufEHKfdG/MU2esYC93HL9/WEEclNCW0elYZsFc5J61Sh\nN3expOcD/wHcCPw6XPwugn7Zi4CnAncSDOH5Sd6xhhb4FJEmsiu3KhcAledirW3FFrl2a4jspnOr\nX08py2LMAptGjwFTTWcXapIyKRDzXL1ZsQ/z7snovj4zZVkbtBgUVVSXbeQ9thIT3ydxd/EcvVmy\nZvZ1sl8xD+yyLG2ycquF7uIia7ZqH2Zp8gKVJmh00iyIJJHrrpAWxDWzIV7R/LkyaWD8ZURVq7bu\nzDFNC27VhP2FffM5zLsnV5TbpzEmSG4xqUXbaN5jmLrhPln5GcJ1bwKOAx4FPmtm70jZfyXBr2ME\nBuJrzCy3Q3AQ0cWLhZVbFWwwO4dsTk7iNvpiawps1QjqpdySmrkKaExcq5Qpb9tWGuaRCG1Elihm\nNeBVJ2monFM4cY9Uvf+y+m5bpWK3QZN9tFBebOPEhbeOFTtwsvIz7AgcBuxtZg9JekpyR0k7A28G\nZszsl5IuAl4JnJt3QhfZDigUV5gnPlkPWp0AjnnjdFck1qVEQWdt2yo1BbbNYVLRsRv/HRrsq60j\ntHHq9vHVnfGodi7hCV7AMusv8uh0FYdQ4QWryT7aqi9KkBTehqZTHAjh8NCN4ecHJG0AdgZeD7zP\nzB4K16UG3BJo5taS/ht4HHBP0TldZFskGWU8u3xFYkGKdddEtG2qCJUc21okMo0JXImyL3g5aOrc\nJdh0botWLTQWhdxHJqmqx8+iTvBbZSs2tn1qd82ZVl5wk3VW9iWgJ6FNo6ylO80k8jO8H/gDSacC\nvwL+ysy+Fd/ezO6W9AHg+8AvgS+Z2ZeKzuMi2zGzDXbygW4wMCaLOoJRxyXXqTuuA1oTWmjMhTyJ\n2EZMmtSiavR74Ytki8Fvs67jFRV3TKuro2hNaCG/TsvmO86i2CsxPEt2w7578NDMY1PX/WKb783L\noRCyKhzyOUuYn+Fi4K1mtknSFsD2wPOA5wIXSdrDYpHBkp5E4FLeHbgf+JSko8wst/ZdZFsmEpw8\n1yyQatnMTlnX4Jv+pnNhuzPzizIpyWutKrrJsrcqciVpzX0MjbuQ4zQx7KfpoWNdiWtRXc1bP6nb\nOK3usq6jpX7aJJN6JoJjZs/dO1TiORTSyMjPcBdwSSiq10j6NfBk4MexXQ8CvmdmPw6Pcwnw+xTc\nsS6yHbHdLyu4oyjZODbkSuuNGo1p39fWidhCowFScdoYa1unHJl0JLC1qeJ5iLZrQGzrjJPOezlK\nCnDmtiOZyassOfkZ/g14EXClpKcDjwHuS+z+feB5kh5H4C4+ECgcO1pKZMOTnkmQV/jZkp4DLDez\nvy2zv1OS8GGbRGBLEeuX3W5F+8JVxYWca4GvaKpEk5PXZ9wIDc4EE6epmYHqnqspUifgWJG/fWMv\nSFVd/EXu5JJi20SXQEQZ78Rs/U2X0C4DXg3cKOm6cNm7gHOAcyTdBDwMHG1mJmkJcJaZHWJmV0ta\nDVxLEKW8nhyLOaKsJfsR4K+BfwEwsxskfRxwkW2Y5ANUdaxgFcHsUrTiMxN10f/cJWm/eWtRyXFa\nFN4smop6b4oqv3PTdVLqWY0o029bUrzrzldbhSGnopyEgvwMC359M7sHOCT2/WSg0rREZUX2cWZ2\nTWBpz/JIlRM55Ef1lp2MvSGB7StHMVDbzR1Zs327jMvQuqUL2fdCiy8xU9H41o0OjrPMFoYELVuY\nSKby71WjvzaLKgJcWM4KM2458ykrsvdJ+m2CLBdIOoJwrJHTDF1asENk7OXPoxNLN06eaEyZF6Ey\nEwjsrLs5p+tjwSTy0fSSWUlYsmgg6rzx6SJdaGtRVmSPI/A9P1PS3cD38Me1Nao+HJMK1Fisw2ki\n+Xt35rovEpW+n+oqw2HqHLsv6lxTw5NM1GLacoj3QCmRNbPbgYMkPR7YzMweaLdYi4tSyfNHHkmc\nJyJNpUIcM51bu1n04ILuhCpDbNL2Oz+sj7qWXJVkF0laCoIrfT5nInJFVtLbMpYDkAiBdtoibVqv\nc+sfLi1Sty9rdlpFswkGMbFBRAvDizIpGvrSBFWPfRTDibBty8It85ukzWjk5FJkyW4b/n8GQRaM\nKF30ocA1bRVqsVK2H7YtYWpKaE//Vcl8zU4tOgmqyqJLizbvXGVFsq4Fm7bvpENZktZs/Ph9uZPr\nvsh4/2xpckXWzE4BkPQ1YN/ITSzp3cBnWy9dWR5cNxdYkLxpBn4jRK7irqM22xh3Gp/SLym0aedz\nK3ZEDM1dPLTyTMKkFnsdsW3CSxAX2qqBXYuIsoFPOxIM0I14OFw2LNJunOjNcUjunpBCgW3Zik0O\nzm89KUXLx69Dcq7fsgzBUm/Vip0mEUtSFFzV5rWnvfSf35BAFfWnN+1+PwoX1xKUFdmPEuRzjNKE\nvAI4r50itcjsfK3DENuqOYnbEqkhil+b1BXWrGN0LbitiOs0i2oaQ7reSYKiyuCBTL1SNrr4VEmf\nB/4gXPQaM1s4229f/ITx3Uh5b4Bju5aQlVs1I2Bt0VbZuuqDblRcOxaZpjMTTUVSDGdRUDZ38VMJ\nkiVfGl9mZt9vq2CSDgbOADYnyB35vrbONQYWg7WZNf/upAxZ+MvQmLg2LKx9TTZQ5tyjE+Gmrdky\ndT3Sl/mxUdZd/FnCbE/A1gTz6X0beFYbhZK0OfBh4CUEUxB9S9IaM7slf8+RMIVWbERdazbKotO0\n0C5qgW1AVNsS0qYmiM9KdN9UuQcj1mWSdFSt7yG5zCtwBQeyAzukrrufczouTTFl3cW/E/8uaV/g\nL1opUcABwG1hEgwkXUgwWW66yG7POG6YRRIkUEVo01LUNSG0Lq7VmUSYmhLNrs5bdo7c9cv2mhPa\ntmM58qzZKM1il+OVnUaoNZ+smV0r6XebLkyMnYEfxL7fBbR5vmEw8qxOccr0UWblgG3DZdw2gxgb\nXKPRrSOsdYVtNUfU2q9pjmB14TXERXj9sr0W5iTugwbyGWcxyQvWYKz9gVK2Tzae+WkzYF/gnlZK\nVBJJxwLHAjx1eIOJqjNyN3Ed0uaZbUpg+wjCalJoa41jbrERhmriWldQN6xZMLdN46xenr8+KcLv\n4dSWS9Q9TXYDVJrybxFS1pLdNvb5EYI+2oubL84sdwO7xr7vEi6bxcxWEU6YO7P/1rZ+2e4LDpJa\n2X0N31kkruLFTiTsTYhtrUnGa06VVqbRjay7MmJ7BKuBamLbqsCeFf4/Jv88ey1fP1vm6Bo6pSgA\nqoEXqaoCm1ffaW734PgutBEyKxYdSf/bzD5VtKyxQklbAP8JHEggrt8C/sTMbk7bfmb/re2CtQtF\nNoteXD8TBDuN0V1clrgl20VkcSR+XVq5TbqRa/XVVmyUm2yEk3Quumclvh9Tbre9ls+NULxl4cyx\n3ZMU3g6Etmq3QFxw99WGdWa2f+WCFSDpY2+4+fVH7TCTHvj0kf3OYeO6jYOyaMqK7LVmtm/RskYL\nJh0CfJBgCM85Zpbpsxm8yBZZsSVcxYtBaNvui80KqOrardyU6HYhuDAc0YUJhPcs8gU2Lsax7azA\ntdwZDYsstCO0Ee/V37nIhuSKrKSXAYcARwKfjK3aDpgxswPaLV45luy/kx2z9rWlIwahI6Et6yIu\n2R87zULbNmUs5j4ikl1w56jTj9uYizlp8YbYmvTlvdCDNQv1hNZFdo4ikd0b2Ad4D3BSbNUDwJVm\n9tN2i1eOSGQjisQ2umlaD2goI7IVA55caMsRn1Eobx7QpOD2OfSnCcGtPfRnsQtuJLKf3jRvsdl2\n9Y/ZNg3EebQltC6ycxTNwnM9cL2kC8zskY7KNDHRDZEmtn2N52uKeCPqgjtHmrjMW1Zyaq4+U0M2\nkQ+59jR4NcZfxgMLhxA4FfWjdhGhPAiWdTcDzuFcOvq2sy+KJm2/yMyOBNZLWtBCmdlzWitZA1zK\n4ZVcyGPDBbeikKRMOJ3WT5sUuD5Et4ko5VrRyVBr6rSlV20obRV1JbZxCoX3GBIu49MBkMDs5NLn\n7oVl5qMXBkzREJ63hP9f3nZB2iIutMmH+iRO7H8MXJl0aSVoKr/tmMS68wnLe6BJsYXhDAeCatbR\nEayeKJnFJFaudIoL7ZQgaVeCWeV2JEgVvMrMzgjnSH898ONw03eZ2efK7Ft0zs3yVprZxvDjX5jZ\nnfE/2k2r2CiXcviwXR0DSo+23Yr5f4uBouvtPZMTgdhGf5Ow6dwaL1LnU+lFsEoygiqepk7GrR46\n4D7YHhl0+1mNR4C3m9kM8DzgOEkz4brTzWyf8O9zFffNpGwyipcA70wse1nKsl74KU9iNUf0M3i8\nKeJCO6DsT3WFNu6qzHJbDsJqLjnzSVxo+86L3ERmqVrWbQXLtqr7uItGfNH01TbMFAlsZDhuDD8/\nIGkDQRrfSfbNHapS1Cf7RgKLdQ9JN8RWbQtcVaZgXZJ0J41WdNMasQEJbxniDXdWI96b6OYEQcVf\nCrJIE7hFlcKxZNahKkLbJvXEdSVRvyyMxGVckabHyc5vf/+uRom6RdJuwFLgamAZ8CZJfwasJbBY\nM0fPJPbNpciS/TjweeD/AsfHlj9gZj8pOnjf9JoerWmiRm1kYluFvIY+Tfhq9zVCrgUbHauK6Hct\nvE27sCsHSbWcJ7kJSotrxhjZaaZIYMuKax+TPnyJP+Kx7JG67n4umZfXPmRVmIZ3FknbEKQGfquZ\nbZJ0JvBegr7W9wL/CLyWFJL7FpW3aAjPz4CfAa8KD/4UYCtgG0nbtDlpe5OM3pUcZ6Bu5b6pHUkL\nmb9pGasHIwjNAAAVqklEQVQ2j7JCWFWMh9BHXJY+rFkX13yaENihzKiURjyvfRqStiQQyQvM7JJw\nnx/F1n8E+EzZfYsoOwvPocBpwBLgXuC3CDJAtzJpe1V+df/jMh+sKKowS2hnWN9ebtK2I/4mtSRG\nItJVhK6y6zP5GyY8BpMKbRnGJJp9UtSwty2ug3YZl2xnJhHYssI65L5vSQLOBjaY2Wmx5TvFAn0P\nB24qu28RZQOf/pYgmuoKM1sq6UUM3lkUsGHN0nlCC5MPB5gaytZgT2I8iMAo6rmPnYCm+v3yntdS\njXpZYf10ofdvlLQtrkMW1gTLgFcDN0q6Llz2LuBVkvYhcBffAfw5gKQlwFlmdkjWvhmRyLOUFdn/\nNrP/krSZpM3M7EpJH6xyZa2T9RDFprZKiq1TkjqvUxMK8ySCVmsu1iQp45ddbLtnYuu1jLhWENZB\nWrMFVmyewE4iriMS1lnM7OtA2g+WKpRmdg9B/v68fXMpK7L3h529XwMukHQv8POqJ2uN+3LWJeaR\nTGaD6WUaq4wUduuX7TU9Ex73HSFdNC/nBExrpq0+xkW31sgXiesUWKzrCYZo5rVgTQtsIy81i4yy\nInsY8CuCuPY/BZ5AMGnAcIkeomhweTjV1YKbpO+prEYQpdkYHUVIx2fcmYgS2biKhGnoItymsLbl\noqwtrqWF9fTiTUZA1u/fuLi6sOZSSmTNLG61ntdSWerzs0ezH6C42GbMGTkkylizTUdrdm49lxDb\nugFHjQlsxIRpL6sOS2qbxkS1pxfDWo19obhWE9WhuYyreMCyBLayuLqwlqYoGcUDBB3BC1YBZoOe\nBypB/EELBbeTuSLLRBgnrNmuhzy0cb5SD32B2HYR2QslGqmWhk01YQX3kvqyAYFt1Jqqbb3Wt1iH\nJrRJqjzTjf3eU+CCb4OicbLbdlWQ5ih6cFbGbobhvCMMITNOkySvp5SINe1GjvfLFlikpa2BLIFp\nwQU+yNzRJQW27v3caFBiZqM/He7gONH9W7UPtpLAurjWomyf7MD5EeUfnPh2HbyJlhm/tgj6ZEuJ\nbsmI3sjCLeUaPjPmSUj5nePlij7Xcp/n1eFIxiNn0vD9WScXbnNWVXMCK50CDGcqPBfYYTIlIjtA\n8sQ1LiZn2myUYNMMIbF31iwrmYKWYXEmrbpKfa9ZLvtlxlJY8PunNVYT9Vt3aP02woSi2pkV26PA\nxhm667gKlYbluMCWoheRlfR+4FDgYeC7wGvM7P5w3QnA64BHgTeb2Rf7KGPrHEXQ+NdkCAJahmQ5\nk6Kb6qZtaI7deUyYfauscFQS4ynzYDSdcD6imYZ/+lzEZZm4ragc5LR4f+s0+rJkLwdOMLNHJP09\ncALwznBuvlcSpGtcAlwh6elm9mhP5azFemZaGX3blbA21S+WlsYyuoa42FYW2owZdArJeKlJK1Nd\nioRmasZBxyj7AtLo/Vup4W+/0e/amj2JE4HDC+/ZKtHE/jLTDr2IrJl9Kfb1mzBb44cBF5rZQ8D3\nJN0GHAB8o+Mi1qYN12+TjVOX2a7yph5MClum0MJ8sa0rsCVoUmyzmDTArW+Rrlv+onvYs7DV49Ic\noW2k3RjYy8wYGUKf7GuBT4afdyYQ3Yi7yJhQd/50Rk+odeKm3z7jAls2WjVPlOs+JG00WFVTqCUz\na0H61IPxRqKwn7YlgX0PpwKRddDcS00bYj22KPRO8hJHLNI+wqr3q7/QdEtrIivpCuA3U1adaGaX\nhducCDwCXFD1+PHpjKQltVvfNt088QaxigVSf6LkyWgyF2nyWHHRTYpt8m0816ptkfdw6qzQNkHT\n7v02Lewm6ev+TcetqziN/N49v8zc/pW94LYM4+Rnj+u2MCVoTWTN7KC89ZJWAC8HDjSzSCTvBnaN\nbbZLuKxVuuhPWSAcYf9g0kXZxVyOjYhpkRspkVErfs606QcXuGonCAqbhDJCW/T7tzV3cdfBbq26\nIZkgN/EizzakNbDX8oVeoTymMdn/WOgruvhg4B3AH5rZL2Kr1gAflxTNXbsncE0PRWyFeZZtYl0b\n4jrxgzNJY5bcNya6yekHk+7jyH3bF+/hVGaYs7zj5StTB5O+BGWdry3xzqJpUa/yu3ijvxBVzFA3\n8fyvi/xlpin66pP9EPBY4PJgHly+aWZvMLObJV0E3ELgRj5ubJHFSbICE07ixFJi0lnD1PYDFZsN\nCeZPP5i8xiHMPHELS2eFtus+rKzzdVWOScV8knLWnuUl1YU5Ha5iRZOYpORbH4RXC5iW37oN+oou\nflrOulOhZ1OmJnFBjVsBbQYm1HpI6grqJH0x8dmQYJ7YRlbthjVLsb5nRSpJG5ZWWrBY/FxZ65um\nj8CYiSZen+KApzyBrYv/1t0yhOjiqaKsoK7mCFazHmLu0ondO2nUEdQ2HqSMqQdhuK7BuDUb0WZZ\ni47d5e/UpqBXvo5aL4Xjt6w0wQtn7Xul1mQL4/+t28RFtiHSgmWy+tOSYtpoUEJnolrlwVq58FzR\n1IPHML/MA7Nk45Z2Jk242vMslR6maBzEi88imHi9Lo3XT+3f2gW2CBfZBphJWKRQXUgjWn/Lr9ww\nNfEQxY+xcq4cyTl+B0qlhOl1KXu8vn6vtsW90ft43A2/lPD6NPnb++/cOS6yDVKnLzXeH1lI64IK\n7T84p7NAaKPPwJCmH8ykbD2U/f0PLbjm5HGKts/at8p+SYbwMtTw5OtN0usEAU3XTan71gW2LC6y\nLVJGOBuzkgYpqEXnXTk+l19a3Ux6DVX3r3u+Jn/rSQS7LKXLO6UNftzbE7dmmxTVXjxbiwsX2ZAm\nE1K0GvE7VQ9FzKoNv0vDmZ+zkInf+FfmrKtKE/VcoTyDeDka8r3dML32Ty+i37kFXGRj1BHaIPvK\n3Pfa4/xgwgdlrA/CwMvd6pylQ7v2ScrT5AtDFkP7vXqiFUFdHL+tpF2BjwI7AgasMrMzJL2XYIKa\nXwP3AivM7J6U/Z9I0Co8O9z/tWaWO4GNi2wDLBDWKu4c7/8YH1Oc+KA+i/36m0E6JfwUi1soxH/7\nCjwCvN3MrpW0LbBO0uXA+83sbwAkvRk4CXhDyv5nAF8wsyMkPQYoTJbsIjsBWs7CyL/alqo/KBFd\nz81ZH68zpy2SXSnx5U5dzGwjsDH8/ICkDcDOZnZLbLPHE1ip85D0BOAFwIpw/4eBh4vO6SKboGwD\nnzpQPBJYf/vsjLk3/4W0ItSzdev15zRH3n0c4Pdb00jajSCN/NXh91OBPwN+BrwoZZfdgR8D/ypp\nb2Ad8BYz+3neeVxkJyVuufqA7U4obpDmtptEaGfHKy7A69OZoz2vi99nqVwMbJux7r7kXONA0O+6\nKr6ZpG3CI73VzDYBmNmJwImSTgD+EkhW7BbAvsCbzOxqSWcAxwN/k1dcF9kaSJsWDmFwC6cTygps\nfPv6jWBsqNGCZY7TDFXvaSef+FzjaUjakkBgLzCzS1I2uQD4HAtF9i7gLjO7Ovy+mkBkc9msTKGd\nAOmU9AfCBbYT6jZGdfbzhs9xpg8F076dDWwws9Niy/eMbXYYcGtyXzP7IfADSc8IFx1IMGNcLm7J\nplCqgZ3nGnZxbZrIAm1K7KLj1LNqvX4dZ0pYBrwauFHSdeGydwGvC8Xz18CdhJHFkpYAZ5nZIeG2\nbwIuCCOLbwdeU3RCF9mJ8Qa4LdqwJsu4j92KddrE76/+MLOvA0pZ9bmM7e8BDol9vw7Yv8o5XWRr\n4cI6ZsYzRMhxnLHjfbKO4zg18Zc1pwgXWcdxHMdpiV7dxZLeDnwA2MHM7guXnQC8DngUeLOZfbHH\nIjpTiveLOY7TBb2JbJio+aXA92PLZoBXAs8ClgBXSHq6mT3aTykdx3Ecpz59WrKnA+8ALostOwy4\n0MweAr4n6TbgACB3lgPHcZwh4h4Tp5c+WUmHAXeb2fWJVTsDP4h9vytc5jiO4zijozVLVtIVwG+m\nrDqRYPDvSyc8fiw/5RMmOZTjOE5lfLy1U4bWRNbMDkpbLul3CGYzuD7IcMUuwLWSDgDuBnaNbb5L\nuCzt+LP5KaUlC6YlchzHcZy+6bxP1sxuBJ4SfZd0B7C/md0naQ3wcUmnEQQ+7Qlc03UZHcdx6uIW\nrBNnUBmfzOxmSRcRJF1+BDjOI4sdx3GcsdK7yJrZbonvpwKn9lMax3EcZ9B89UEga67n4dlknvHJ\ncRzHcVrCRdZxOsBz3E4/mfNNO4saF1nH6QBvfB1nceIi6ziOUxH3TDhlcZF1HMdxnJZwkXUcx2kA\nt26dNHofwuM4jjN2vM/dycItWcdxnIq4qDplcZF1HMdxnJZwkXUcx6lBZM26Vevk4SLrOI7jLAok\n7SrpSkm3SLpZ0lvC5e+XdKukGyRdKumJGfsfLOnbkm6TdHyZc7rIOo7jOIuFR4C3m9kM8DzgOEkz\nwOXAs83sOcB/Aickd5S0OfBh4GXADPCqcN9cXGQdx3Fq4q7icWFmG83s2vDzA8AGYGcz+5KZPRJu\n9k2CucyTHADcZma3m9nDwIXAYUXndJF1HMdxRoTVXDcfSbsBS4GrE6teC3w+ZZedgR/Evt8VLsvF\nRdZxHMcZC9fBVzJW3Qpsj6RjJa2N/R2b3FLSNsDFwFvNbFNs+YkELuULmiqwJ6NwHMdxxsI/wwc/\nAA8Cj0msWg3sj5mtAlZlHUDSlgQCe4GZXRJbvgJ4OXCgmaWZxHcDu8a+7xIuy8UtWcdxHGcUmNnP\nYW9gXWLNvcDPMftX5e0vScDZwAYzOy22/GDgHcByM/tFxu7fAvaUtLukxwCvBNYUldlF1nEcxxkR\nl28D64GHY8u+BvxhmZ2XAa8GXizpuvDvEOBDwLbA5eGyfwaQtETS5wDCwKi/BL5IEDB1kZndXHRC\npVvF7SPpTcBxwKPAZ83sHeHyE4DXhcvfbGZfLD7WEoM/b7O4juM4Tmnevc7M9m/r6NJLLbARf4/A\niv08ZrfnWrF90UufrKQXEYQ+721mD0l6Srh8hsAEfxawBLhC0tPN7NE+yuk4juMMkcu3gac8CPtR\nwYrthb7cxW8E3mdmDwGY2b3h8sOAC83sITP7HnAbwdgkx3EcxwHifbNfoExfbJ/0JbJPB/5A0tWS\n/l3Sc8PltcYhOY7jOIuNy7cJkjMN14qFFt3Fkq4AfjNl1YnhebcnSGv1XOAiSXtUPP6xQDT+6SF4\n900TFHcoPBm4r+9CNIBfx3CYhmsAv44hUeYafqvtQgTWLIO1YCNaE1kzOyhrnaQ3ApeEY5GukfRr\ngoorPQ4pPhZK0to2O9m7wq9jWEzDdUzDNYBfx5CYhmvokr7cxf8GvAhA0tMJRhXfRzDm6JWSHitp\nd2BP4Jqeyug4juM4E9FXxqdzgHMk3UQw2Ono0Kq9WdJFwC0Eqa2O88hix3EcZ6z0IrLhDAZHZaw7\nFTi14iEzU2iNDL+OYTEN1zEN1wB+HUNiGq6hM3pLRuE4juM4046nVXQcx3Gclhi9yEo6WNK3Jd0m\n6fi+y1MFSXdIujHMlbk2XLa9pMslfSf8/6S+yxlH0jmS7g3706NlmWWWdEJYN9+W9Ef9lHohGdfx\nbkl3J3KaRusGdx2SdpV0paRbJN0s6S3h8lHVR851jK0+tpJ0jaTrw+s4JVw+tvrIuo5R1cdgMLPR\n/gGbA98F9iCIUL4emOm7XBXKfwfw5MSyfwCODz8fD/x93+VMlO8FwL7ATUVlBmbCOnkssHtYV5v3\nfQ051/Fu4K9Sth3kdQA7AfuGn7clGJk/M7b6yLmOsdWHgG3Cz1sSTAb+vBHWR9Z1jKo+hvI3dkv2\nAOA2M7vdgmCqCwlSM46Zw4Dzws/nAa/osSwLMLOvAT9JLM4q82DTZGZcRxaDvA4z22hm14afHyCY\nGWRnRlYfOdeRxVCvw8zswfDrluGfMb76yLqOLAZ5HUNh7CI79jSMRjAJwrowgxXAjma2Mfz8Q2DH\nfopWiawyj7F+3iTphtCdHLn1Bn8dknYDlhJYHaOtj8R1wMjqQ9Lmkq4jmBrmcjMbZX1kXAeMrD6G\nwNhFduw838z2AV4GHCfpBfGVFvhiRhX+PcYyxziToOthH2Aj8I/9FqcckrYBLgbeamab4uvGVB8p\n1zG6+jCzR8NnehfgAEnPTqwfRX1kXMfo6mMIjF1kS6dhHCJmdnf4/17gUgIXy48k7QQQ/r83+wiD\nIavMo6ofM/tR2Lj8GvgIcy6vwV6HpC0JhOkCM7skXDy6+ki7jjHWR4SZ3Q9cCRzMCOsjIn4dY66P\nPhm7yH4L2FPS7pIeQzAX7Zqey1QKSY+XtG30GXgpcBNB+Y8ONzsauKyfElYiq8yjSpMZNYQhhxPU\nBwz0OiQJOBvYYGanxVaNqj6yrmOE9bGDpCeGn7cGXgLcyvjqI/U6xlYfg6HvyKtJ/4BDCKIRvwuc\n2Hd5KpR7D4KIvOuBm6OyA78BfBn4DnAFsH3fZU2U+xMErqL/Juh7eV1emQlmXfou8G3gZX2Xv+A6\nPgbcCNxA0HDsNOTrAJ5P4Hq8Abgu/DtkbPWRcx1jq4/nAOvD8t4EnBQuH1t9ZF3HqOpjKH+e8clx\nHMdxWmLs7mLHcRzHGSwuso7jOI7TEi6yjuM4jtMSLrKO4ziO0xIuso7jOI7TEi6yjlMRSQ8Wb1X5\nmMsVziIl6RWSZmoc46uS9m+6bI7j1MdF1nEGgJmtMbP3hV9fQTCzieM4I8dF1nFqooD3S7pJwbzA\nfxwuf2FoVa6WdKukC8KsRkg6JFy2TtI/SfpMuHyFpA9J+n1gOfD+cM7O345bqJKeLOmO8PPWki6U\ntEHSpcDWsbK9VNI3JF0r6VNhXmDHcTpmi74L4Dgj5n8SJEvfG3gy8C1JXwvXLQWeBdwDXAUsk7QW\n+BfgBWb2PUmfSB7QzP6fpDXAZ8xsNUCoz2m8EfiFme0l6TnAteH2Twb+D3CQmf1c0juBtwHvaeKi\nHccpj4us49Tn+cAnzOxRgiTw/w48F9gEXGNmdwGEU4btBjwI3G7BnJsQpHY8dsFRy/MC4J8AzOwG\nSTeEy59H4G6+KhToxwDfmOA8juPUxEXWcdrhodjnR5nsWXuEua6drUpsL4I5QF81wTkdx2kA75N1\nnPr8B/DH4QTXOxBYlnmzj3wb2COcmBzgjzO2ewDYNvb9DmC/8PMRseVfA/4EIJzv8znh8m8SuKef\nFq57vKSnl7gex3EaxkXWcepzKcGMJNcDXwHeYWY/zNrYzH4J/AXwBUnrCMT0ZymbXgj8taT1kn4b\n+ADwRknrCfp+I84EtpG0gaC/dV14nh8DK4BPhC7kbwDPnORCHceph8/C4zgdImkbM3swjDb+MPAd\nMzu973I5jtMObsk6Tre8PgyEuhl4AkG0seM4U4pbso7jOI7TEm7JOo7jOE5LuMg6juM4Tku4yDqO\n4zhOS7jIOo7jOE5LuMg6juM4Tku4yDqO4zhOS/x/7FwU6sqFOb4AAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fe232290630>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"#dclim = xr.open_dataset('data/clim.nc').sel(Y=slice(-20,20),X=slice(120,280))\n",
"dclim = xr.open_dataset('data/clim60.nc')\n",
"\n",
"plt.figure(figsize=(8,3))\n",
"dp=dclim.sst.plot.contourf(levels=11,cmap=\"jet\",vmin=22,vmax=30,add_colorbar=0)\n",
"#,vmin=-2,vmax=2,add_colorbar=0)\n",
"plt.title(r'$\\bar T_{s,O}$')\n",
"plt.ylabel('latitude');plt.xlabel('longitude')\n",
"cbar = plt.colorbar(dp)\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [],
"source": [
"# set Q'_LW + Q'_LH = 0, solve for Ts' (assuming U'=0)\n",
"# Q'_LW = ALW(Tsbar,Cbar,f1bar)* Tsprime + BLW(Tsbar,Cbar) * f1prime\n",
"# Q'_LH = ALH(Tsbar,Ubar) * Tsprime \n",
"Tsb = 1.0*dclim.sst + T0\n",
"Ub = 1.0*dclim.wnsp\n",
"Cb = 1.0*dclim.cld\n",
"f1p = -0.003\n",
"\n",
"ALH0 = f_dQLHdT(Tsb,Ubar)\n",
"ALW0 = f_dQLWdT(Tsb,Cbar,f1bar)\n",
"BLW0 = f_dQLWdf(Tsb,Cbar)\n",
"dTse0 = - BLW0 * f1p / (ALH0 + ALW0)\n",
"dclim['dTse0'] = dTse0 \n",
"\n",
"ALH1 = f_dQLHdT(Tsb,Ub)\n",
"ALW1 = f_dQLWdT(Tsb,Cbar,f1bar)\n",
"BLW1 = f_dQLWdf(Tsb,Cbar)\n",
"dTse1 = - BLW1 * f1p / (ALH1 + ALW1)\n",
"dclim['dTse1'] = dTse1 \n",
"\n",
"ALH2 = f_dQLHdT(Tsb,Ubar)\n",
"ALW2 = f_dQLWdT(Tsb,Cb,f1bar)\n",
"BLW2 = f_dQLWdf(Tsb,Cb)\n",
"dTse2 = - BLW2 * f1p / (ALH2 + ALW2)\n",
"dclim['dTse2'] = dTse2 \n",
"\n",
"ALH = f_dQLHdT(Tsb,Ub)\n",
"ALW = f_dQLWdT(Tsb,Cb,f1bar)\n",
"BLW = f_dQLWdf(Tsb,Cb)\n",
"dTse = - BLW * f1p / (ALH + ALW)\n",
"dclim['dTse'] = dTse \n",
"dclim['ALH'] = ALH\n",
"dclim['ALW'] = ALW\n",
"dclim['BLW'] = BLW\n",
"dclim['QLW'] = ALW + BLW * f1p / dTse\n",
"dclim.to_netcdf('Q.nc')"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjQAAAGoCAYAAABRx3+mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsvXucJGV1//8+LAiss8slg6MQyYIRYb3gAl4S1EiiCfZE\nDUlPowm6a0zM4KLRb37xEr+xd0gweImoLzasJvG76wtvPRM3QaZDvEezQYVlAZUVNUKQiwOjXHa4\nCAvn90dV9TxdXdfuquqqnuf9es1ruuv6VHX3eT51nvOcI6qKxWKxWCwWS5U5YNgNsFgsFovFYhkU\nK2gsFovFYrFUHitoLBaLxWKxVB4raCwWi8VisVQeK2gsFovFYrFUHitoLBaLxWKxVB4raCwWi8Vi\nsVQeK2gsFovFYrFUHitoLLkgIp8XkSX37ytF7WuxWKpHFX/zVWzzqCM2U/BoIiK/C3xVVe+P2e5p\nwGeBpwDvUtWPFNE+i8VSHspuL0Tk74AFVf1QTsf/NvA6Vf1eHse3FIP10FQMEdloPBU8JCKPGu/v\nEZGD3U0/BByY4JBvwzFka7IyTiLyAhG5V0TEt/zLIvL2BPuvEZH3iMiPRGSfiNwkIheLyFFZtM9i\nWSmMgr1wf/evBT6axflC+ABwftKNrY0qJ1bQVAxV3aGqY6o6BrwHuNx7r6qHq+ovRORE4Buqem+C\nQ/4K0NdTiYiEGcANwLXa6/57NnBNzDEPB74BnAi8TFXXAC8EDnLbarFYEjIi9mIT0FbVB/s5b0Iu\nA84QkSfGbWhtVHmxgqbaPBu4LmD5b5LgacYd9z0DuNh9YjvBXX6SiHzNfYL7noi8wtjnZvep6Xrg\n/hAjtQHY4zvXOuBI//IALgJ+DtRV9YcAqnqrqv6Zql4dd00WiyWUqtqLlwH/6Vv/PhH5V+P9+12P\nzuNC2j7meqeeZCx7hojcISJrVPUhYDfwO3H3AWujSosVNNXm2cC1AcvXqOo343ZW1d/EedI4z31i\n+4GIHAR8HvgC8ATgTcAn3bFzj1cDk8Dhqro/pF1+T8wG4FZVXQxrj4g8GXgNztj8Y3Htt1gsqaiq\nvXgmcKNv/XtxPCobRGQaOBP4fVV9OKTtS8D3gVOMxRcC71HVfe77vcDJ4XfA2qiyYwVNRRGRtcA6\nAgyUqr53gEM/HxgDLlTVh1X1K8DlOEbJ4yOq+pMgF7Br4J5OryfmFGKGm4CXAHep6pV9t95isfRQ\ncXtxOLDPXKmqP8PxlOwA3gnUEgyZXeUeFxF5EbCebs/UPvdcUVgbVWKsoKkuJ+P8AG/K+LhHAz/x\nPX38L3CM8f4nEfuf5P7f61v+XByDEsUEcEuSRlosllRU2V7cDawJ2HcPjvfmnaoadQ6PjqAB3gf8\ntc+jswa4J+YY1kaVGCtoqsuzgesDAukG5XbgySJifjeOBW4z3kedcwJ40HQti8g48Bs4gXdR3AIc\n4zu3xWIZnCrbi+uBE8ydROSZwCU4Hpo/TtjWq4BTROQPgEOAT/nWn0RwjJGJtVElxn4o1aVnPFxE\nDhKRT4jIV0TkShGZ6OO43wIeAN7mHu/FwMuBzyTc/zvAwSLypyJyqBs4+Gmc2RXXu+3cLiLbA/a9\n3P1/oesiR0ROEJEPu0bOYrH0R2XtBdDGETheu4/BiduZBt4IPNM9r7c+zL5cBzwR+Hscr44a+xwC\nnAp8MeY41kaVGCtoqsvJ9I6HPwfY5wbv/bqqLqQ9qOuCfTnOzIJF4B+A16rq9xPu/1OgAbwZZybA\nFSxPvfR4MrArYN/7cGZcnAD8UETuAXYC90cFE1sslliqbC8+AdRcwbMWR+B8UFUvU9UHgPcDFxjb\nh9mXX+AIqJtV9d99q18OfE1Vb486jrVR5cZmCh4h3CRZM8BxwKWq+vkhN6kHd1rldcCzVPWRYbfH\nYlmpVMFeeIjIe4A74zIFR9kXd92PgIZ/VpeIfAt4vap+N+44lvJiBc0IISKrVfUBETkS+Lyqnj7s\nNlkslnKy0uyFiFwAHK+qr47d2FJJkqS6tlSHj7t5Eg4E/nrYjbFYLKVmRdgLETkF+CpOcPFZQ26O\nJUesh8ZisVgsFkvlsUHBFovFYrFYKo8VNBbLiCAiZ4p0Vyy2WCyWLBGRJ7rDeKVjJGJofklEjw0x\n46tOOZVHr9kduu+qU07NqVXl4Pbd4ddeBBMputcsP4uk1330qcvnzOJepbneKK5VFlX1qKTb/7GI\nrgKmnLexrRCRM4EPA6uAf1LVCwO2eTHwIZwqwouq+htJ9x1V1owfrE9Zt58Df/wYj9wNd7nLve+R\n9x06+nic/LbQ2c78ro0acb+do1i+Vx5HHwEc4bzef8Tys/W9HAbAYQRXMjjwbjcp8d1wu3uPjz71\n1J42eMd/5Mfd5876c/DOexRwkHFN3O2+dtvI8b22xvuBe/vtP+KAzvV7HMa9nWv+2RFHdJb7t0vD\nsUay4+t2P5bK1myAO24GRORQt6hnJEXampGIodlwgOh/Hhy8bu2Dyn2H9tr3tQ9W/7qjmCnJg/pb\nD0m3fZafS9J70DR+A4Pet7TXG8ZhD7FbVU9Luv1TRPTXgK8Bt8EBURlhRWQV8APgpcCtOBlUX62q\nNxjbHA78N3Cmqt4iIk9Q1TuT7DvKHHfakbrl6pcySZvxuSUWpmCbu66pyowIzVnn/cIUTEwBrZVp\na6bd/xNTMDPbu745C4v1sa5l89R6tpukHXj88bklaMGC79gTU0Bj+djedtDdjmaGfZ93D6a983s0\n6GrjhM/WNFVZcPedmILFltPmoPvg0Y5Yl4QabTbOtTrvZSq5rRGR44+F/3kq8DigrRppMIu2NSPh\noYnCipl4PENcBu47VDL7fPoRKoOIm6zETFpE5IVPB56KU5L4NngF8G8RuzwX+JGq/tjd/zPAKwHT\nUPwh8DlVvQVAVe9Mse/IM0+NyXqbiaklmt7ChtCccjrSR2WJbUBzBYmZ5pRvZcP97/ad3m9rRqQj\ndkzCOnFzuSluFutjUIeJuaWuZYv+Y9Wdv0naNOeWmJlabntWosY7zoJIRzRN49R1iNtnwrO/s9Cc\nW+oReR6DChlwxEyYQEzCBvif9Ti1LT6Ok2E5xktTqK0ZqqARkZtxCqY9CuxX1dPcnAifxakMezNO\nEqS7w45hGYw8xUy/HXyWosYj7DqjDFqaezMsMePylhe6L14IXA1vIVrQHEN3wcBbgef5tjkBOEhE\nvoZTtO/DqvqJhPuWjqxsjTcM0vHQeE/eRmc+3lgidcrdEWBm1hU1jYhtTG+E0XFHeSRM/NtN0max\nPpZo/3lqUIfN2u4IzixFDTjipOkedxvArCNstpHQI9SCcZaYrLcT35O0jHsCsBW9nR8ROXgC52lJ\ncAzEAnxIREzvzsdU9WPG+0JtTRk8NGf4Uka/A/iyql4oIu9w37897UFH3QsThfkkFEfWYiarjj0v\nUQPprjlK1AxZxJgc8a6DYa3bzIse4ngRudpY7zcySTgQp7bNbwGHAleKyDejdyk9A9uaA3/8GBvn\nWixMOSl2p+keRvCYgGXPzQjj/03NzMLm1rJQGZ9bggZMz9IZWpkm2dBKEvrZf54aG2dbnYCzrEUN\ndNuNbTHbetuDc48moONR8q4vC+9MBqxaB5zrvhHgG/BfqhrkcEtDZramDILGzyuBF7uvd+CEBaQS\nNCtZzAyTEnXwA5PnENzaTd3v79ue+SkWY8bEb8OpU+Pxy3RXRwbnaehnqno/cL+IfB2nHtCtCfat\nCn3ZmgXDGxMkZlY648awyWJ9bNkjgOPByUrMDMJifYzNCltlKX7jPkns4W0423jevoVZJwI2Lzqe\nsbq7IMd7QMG2ZtiCRoEvicijwEfdp8gJVb3DXf9TQoYhReQNwBtg+YqtkCmePEVMHl4aiHf9JhUz\nFz2U/Pr9IiZ2XZLHuv65CniqiByHYyBehTOObfJvwMUiciBO/N/zgIuA7yfYt4xkYmt+2V3mDS1Y\n4lmsjzExu9R5XYSQ8Xs0ar64keU2pBx3SUmcrVkQ6QxXTkw5YmZiyqnyaVKjnYmXxn+MQeJpElKo\nrRm2oHmBqt4mIk8AvigiXRVaVVVFJPAb4Rqkj4Ezy8mKmW6CgqGzoGgvTF6iJg+iREuZUNX9InIe\n8B84D4MfV9Xvici0u36bqu4VkStw0sU/hjNl0ivc17PvUC4kHZnYmtOeIjrxP9X4PhaF+QAQNvbg\neQWyFjNJO3lvO7+wac7CzNTyNWQ99JSIxrIna8J9Px4VHPzAZODy2ur5RKdrU+uZ6ZQXRduaoQoa\nVb3N/X+niOzEiWpeEJEnqeodIvIk4M7IgzBg/pJzAzr+S0bDYHni46LYTAHR+680yjTrKy9UtQ3d\n1l1Vt/nevx94f5J9y05WtoYj+rc1aYPSq0bPLCcf/YiZrGNHvM7cY0e9QXO21Zn5NAwWpmBi1hEw\n4/Q//NN+YDKxqDHPnTdF2pqhCRoReTxOvox97uvfBs4HLgM2Ahe6/6NmawRzrjiixBQrQSIlSMyY\n+1eUfrwzK1W8mJj3Lep+eALR3KYq3pmVSF62Jo1AMacq5zuaWCwdz0YGYqaowFe/qFmsjzHNUqGf\ni/+702wB9f5mfnkkFTPe1O3F+hjbBhBQZWSYHpoJYKebqf1A4FOqeoWIXAW0ROT1wP8SOQnQwC9O\ndo32E3YYQWLG7HiDOuMyECYIcgiYHRgrZCpHprYmynsXNGOmM1V5FidkYzZ6+6rgFzOdKey+u1gm\nMRPEPDU2TrW6Ppc8CRyiayRLpjfocNNWzmOrLLE10dbVY2iCxk2Wc3LA8p/hTN/KFk/gnO4ajzDv\njEfFvDRJvTJlEDJJRUBRYibu3nnt9bencx3nuP8vza5NluzIw9Z4uUWCCIvHCBvWqJqo6Umm16CT\nMdkfQ2N20iWZegz0emmKoksEGtmMdwRMz/bal5WYAbpmnPXDQeMwcYbzes0PgOsGOlzmDDsoOBuW\ndi93KiZByyDYe+MfoqoSIe1euykbUTAsL0RW7e/33EmWdX3Hgr5vaUSOKaC3pfsurv0jWPs4b99U\nu1pS4NXgmZh18sxExV6MYhyWX8wstpanZXtp//0zmcokZEzCAoXzYsaXhycox4wnWmqr5wOXm6SN\nlzGv0xTQW0boezoagiYNfjFjdiJV8MhEiJc0y6vC2k3k6y07V+KFyiAMKnIspWPa9zqpfqySFwai\nBZlZh8n0zngddVlFTCgNYLYYb9kV2ogVK2HeGEgvZFYSK0/QeGTVYVlGgyK/D/5zeQJnlzhDois0\n/qsKHIVXcHI5ZqTpK75oipw0gqesBF1Dp/BjY4kZd9k2YF2ImInqoKPIsvMO8sS03fZO1tuQY4Cs\nFxS+SsfYHDOEFEbcvTCvL8775Am3UfMirlxB4+HvPE4v+VNUmqGxrDrpMngU8vDS7JJyCFuzDQWJ\nGRE5E/gwTn6Hf1LVC0O2ew5wJfAqVZ1zl70V+BOcZHXfAV4XU6BuZLiLgMrRs8uixuu0vOyrE7Mw\n7Xowema2lNxj43V4E1PQbCwPr03jChnffWhOwdkxnoc05C1mvOVtau7wj5OXJS8vzcQUnM3Fqe5L\nknsQdG1h1+sIUee7WZSYKdLWjI6g8QuRlfaUm2fHXFSnHyecchp62nP6SaHrNuzam/n5ho2IrAK2\nAi/FSS9+lYhcpqo3BGz3XuALxrJjgDcD61X1QRFp4WTw3F5Q84dO2Cwmcxhms7qp/4FVCs0AAVAl\nzI5wYqpX1E0DO1qNgUUMZD+kEhcj4633kuxljVfDarE1lulQ0jCCmtNStK0ZDUEzNkBiPUt5SDJb\nKANRs4f1zosIIdPZ1t2mMGHTEea5CvLnAj9yZ/8gIp/BqWt0g2+7NwH/AjzHt/xA4FAReQRYDdye\nZ2PLxNGnBtuahVk6RQWhezbJYn2MxdZYZUWNV1uo6U07bzgeG7NqwI5Wo++4mTxjQtJ0+kHei0E9\nNabYPZuLB7rWLATMeMOoq+VeW85BwYXamtEQNEGMciyCOezk854EeRsq52U4h8yHufawng09v6EU\n+8eIn0zucXHDnccAPzHe34pTP6WD+3R0FnAGhpFxywd8ALgFeBD4gqp+gRVMU7XzFD7ecKpLex1/\nZxnLMTdVw4uNoU6nCjTAfL03XsbrsE1PRFWCWOepsVnbAxWs9GJlzKKlzVk4u7491XHy8r504r76\nFWrjLGdruhy4jjNE5C3GFh9zS4V4FGprRlfQwHIHYQobs9Pw56aJwHuqH6RTzJQEQibJOj+lET9R\n3hoveDaCjhcmZllW+O9x4ffxbODx7uttjIvI1cZav5FJwoeAt6vqY2I8wYnIEThPWMcB9wCzInKO\nqpYh0mpodDqwhjip7KeAurNsYUo6QbWxHUlDmJktZ3xN18ycGG9MWURMP/lm0g49mdOx/cuA0JpM\nRQ8ZTdJmK+H1tvrkq6r6+gGPkZmtGW1B4xHW+fmWJ3mKH/RJPxNy9DylET9JGLhj7yN+J0/hkrgN\naYeqsvXOLKrqaRHrb2O5SD04RaRv821zGvAZ18CMAzUR2Q8cBNykqncBiMjngF+nHKHjw6elPSW7\neypzN8QZoooQLUMtlhhD5aZkJ8CsOr1YH6M55QwPxgYIN5bFjBc8DYBPEBUtXsKEXHMKaBX6nSrU\n1qwMQZMArxNMIljezbsAOJ8Lcm9XHFkLkKzJbQhslyyLhmELzAj2nH5Sebxey1wFPFVEjsMxLq8C\n/tDcQFWP816LyHbgclX9VxF5HvB8EVmN4wb+LcD0BlliWJh1hnGaQSsbbsI+d2igDMJmGnhzhUVM\nnJdmMmCdGfMU+Rk0YHNrrKdswWZt86gsOYkGB2t+KvxetJ7rLlbMQMG25oBs275yOIud7GF9z1/e\nmOcou5gJI892l8E742fP6SdFX3PBqQJUdT9wHvAfwF6gparfE5FpEYn0SKvqt4A54BqcaZQHAGmH\nsyyETJuta2BFqRkRxxvQEGZEOjE7RbTvCnUalMUMpmHh5ZtJw2Kre6jI/3nNiDO8GJRIcJ6aU7sr\nWSXCgQm7PnN5P9XOB6VoW5NI0IjICSLyZRH5rvv+WSLyf5PsW0U8D4xJUi9AIaImroOsAFm2v4wi\nJojAax5S3iNVbavqCar6FFW9wF22TVV78sCp6iYvL4T7vqmqJ6rqM1T1Nar6i6zatRJsTZJEe/7q\n1Z3aP60xmrNuocu5/ESNJ5jWaf+zl8pIUMcf1dE3YwK5zc/SO7Z3fDP5YNg2/WIeJzaeyfXSbOSz\nA52zX4q0NUmHnP4R+Evgo+5JrheRTwF/m/SiRonzuYA97AxdHzdsZQqmqGGrqnTU/TLIlOiqCzpL\nKCvK1gTGaNQVWsuzZWbEiblZ5Q1t1GHjXMuZSVWPPrZH3JBVT9I/Nx5kR8rrqQr+4Zh5aoFDT36C\nPq9JNzlf4LCWkVHaLzyCClGGrR+ESdpD8c4Mg6SCZrWqflu6v/T7c2jP0NjADexhPTs5C3BEh19s\npInVCNrfW55ku5VEGmETJmSSxD55ny04Q4aWUjLytsab5h3pqfGCi+ekk7Cvq5hhHTa2WoljbKK2\n84sZrybTZi5ervjsDje1H5gsdAaTf5gr6NxbOQ/oHvpJUszRXx4gTNT4Sz+IzKC6HAG1VZaoaZub\nZTkxzzptME+NjY0WEziiwhMWYUNDWRE0ZX5Y3pmiSSpoFkXkKTjphxGROnBHbq0aEmaHF0WU5yTq\nGEFDWd5yv6gZxDuT9Dr8DLuTH9Tr4txHR9T475//nvjfD/vaLR1WhK2B5Y4yaiaNN2zh7wjb1Jhs\ntTslFvz4sxeH5VYRmWELy6Ub/DlTOp3+6vncRE2a2ByzGrUfz1OS5Dzm/qZ3xRONprCZmILp2WVR\ns4UtiIC6n9mMSJeY8fbfzMVO3p56m/G5pc5n2RGJGYiYJPdufG4p0pOXhkfGD2CxvhqA+3/yMPBw\nNgfOiKSCZjNOMM6JInIbcBPlqIKTKedzQZfoSCs0+hEzg9CvcElyvKp18KZnDeB81+OWdn/I99rz\nyjx83fNPYGztKvdd6WZVpWFF2JqJuMKADS9RX3hStnlqNKdaPQn7OmJmCnbUG+7rVuw05M06xg5q\nXaLB9H6wuttT45FU3GQZVGy2cTMXpz6PXxjFeWsmpgDjPi+Lmt65ak23ZpN5XE9QJBExg9ynwM8i\nZmhylEgkaNy0xS8RkccDB6jqvnybVU7Sek3SDCf5p4JviOiQsxYyUeeogrAJuh/O/ezvPhVx7ctT\nzi0mK9XWmEJjRqQTFByWlK1DA7bNLk8BN5O8nd3aDjid9GJruWCmeS7VJjOyhZlZR0maHpmtnMd4\nY6lnto+fYc5+GvTcfm+TPxZmsT7GOEvQ6h162sIWZmRL1/stOPdyvWzqeLt6YmfCPEU5FPOsQr2n\nLIkUNCLyf0KWA6CqH8yhTaUjiZDJ2ztThIiJOm9ZhU2e92UnZ+V63WXOn1M01tY4eIJkYdaZLt02\nPCYepldho1Eh2qM55dRWMvG8OUG1pLyq2uMNp+P2jj8+5+Rhac4t0a73Dj2NCmGixu+lCYt58oTO\nFrbQVHWH8bY42YY1WqhkdS/DxMwk7eU6Y0Mirtq2iLwS+BvgMZx4ubeo6n+56w4H/gl4Bs4w9B+r\n6pVh54rz0Kxx/z8Np8bCZe77lwPfTnFNlSFo2ImcOzaTMiXtM8m7c09LVkIm7Jq845dd0I0QK87W\n+DFFySodixQzEP70vdgaC0ywttgaY1qWOnE7Hp0hKF/OFC9j7mJ9rCvuw1+vKagzrbro6RI19Tbj\nLHUlPDQxh6NmRNgScsyl7eOMbco+zV7Ud2J8bokZQpI4FkDCattfBi5TVRWRZ+EMkp3orvswcIWq\n1kXkcTgFKkOJzEOjqjOqOoOTrvgUVf0LVf0L4FTg2D6urxKkFRN5eAnyiLkZlJ2cNTRPkb8deXOW\ne7V5nnOY0/JF5EwRuVFEfiQi7whY/0oRuV5ErhWRq0XkBca6w0VkTkS+LyJ7ReTXBm3PSrQ1TdWu\neJbmrBPH0pyFzVycaIZPGjrJ3kLwD2/NU2OxNcb43JIzTdxHbfV8aJui1pWNuHbOU3PuTcP5jJpT\nxl9MnpqtnMfFqzdTWz3PXW9cw772UexrH8VNq4/r+vO2GaTtppiZpL08tZ+QJI7F0Km2raoPA161\n7Q6quqTa+SE8nuUJAYcBLwL+2d3uYVW9J+pkSYOCJ+gOZ36Y4XqxSkNcR1fl4aYwhjlDKO29iBMl\naTxPeXhrhjFtv+inppSseFvjBfJG4XVejghxYmO8ekI7QqYF12g7XpfZJRamlodQvArR/uBaD69Q\nY3O2lSq4FYKrb5eBOPHg93R4nhqgJ8B2nCWmXWEzMdXtxdkqSzSnWtw822IrywkTx1nyFcBssVXb\nbHZ/SXH3y99+v5gZn1vqyoEDvdPNM2LgattO2+Qs4O+AJwDexR8H3AX8PxE5GdgN/Lmq3h/WmKSC\n5hPAt0XEs+S/x+jmXAJ6vTRBifTyFhplEzJhhLUza6GT1/3YyVmJcgaZ21d8CKrz1AQgIt5TU0fQ\nqKo5zzfoqWmTu13WczdXnK0xvTTznA1E51/xOi8vvsOrDu3Vh4qavjxPjZuner0tXnkD//HNLLcz\nU46o2VFvJAo27UocFyEgihA7SbwfcdfkT07XNQvKEw69t7Z7mMqopdQ0JpvNiLBVltiq5zn5f0KC\nheMCf/1ixh/EPOjg090cyTy/CcAN3Ahcl0W1bVR1J7BTRF6EE0/zEhx9cgrwJlX9loh8GHgH8Ndh\nx0k6y+kCEfl34IXuotep6p5BLqBKjHrG3rxIKkDi4ljywDxnWHCuP57KZFAvlbl/Fl6ay3k5B3Ow\n+27vuIiYRdyG+tSUhpVua/xxM1FDCh6L9TE2q5Ni35ut5McLEH3UjaHxaE45MTfzxvE3Npzg4a1B\nDXSnAIdl1TU7/SiBkFTshO4fIYL6GrrpYzaQFzS8WHeG5eKISn7oBWZvlSU+O7upk0jRjFmKa7f/\nMwkKYo6tHp49Saptd1DVr4vI8SIyjmOXbnVrOoFT16lneNwkkaARkWNxiobuNJep6i1J9s+bB/he\nV5bfIlz4VfGelJk4ETCoFyQs5qef43r7hH3uSYaj4mbCOSU1MhHPi6p62qAHyeqpKQ1ltzW3796d\nqqRAGs72OaKivDJBeBWidwSIio1zLd8Qh4MngLzjPipOEGkoDcMLEMBGn4siaNr5fEiZgDSJ5soW\nn5NE1CQd8nE+pyWmaTmzpELuS5CYCRpq6jm+K2oKiquJrbYtIr8K/I87vH0KcDDwM/f9T0Tkaap6\nI0617cipoUmHnOZxXc7AoThPaTcCT0+4f2ripnpFkeVMoaAOxoqZ9AxjiCaLcy5/h8K9NSaDZGmO\nqg+WMYU+NaWkUFuT1s4cfTw0T10eRnBmtWzJJDah5tUEihAyQZ2mWRaBVneSOE98LASIGQ+vEzSH\nRrycKlHbJ2G8tdQze6oTi+Iyb7Q3jMoUyGw4NbDM++Pd17ghH1Mcz7ilMTbHZD/ul6KChFV1v4h4\n1bZXAR/3qm2767cBfwC8VkQeAR4EzjaChN8EfNKN1fsx8Lqo8yUdcnqm+d5VUW9MflnpSBi02OEe\nDi9MZFgxk5yKx5n0EDUEFUSJr7/Qp6Y0FGlr0tqZDg26ssZmwQ43dsYjVMj4hUTDWeeJBNP7YYoZ\n//BDl2BpBU9HDmJhygl89TwA23CCkSOZdac2u4yzLMoW62NdXqewIopZlgsIwp8pOIxID5lXXXvO\nJ+IMQZNU9nZyA80tUasnFzVeIsAJgqWT+Tl7AcpbMv4u+1HVNnTfOLPStqq+F3hvyL7XAom9zUk9\nNP6TXCMiPWPuGRIbtGhyOPeEdB79e2jK5pnxPAVlnM5tUuJOvBCqcP1FPzUN2NY8bU0qO2PiBeJm\nwSJrIGCYyD+MEEiLjqjxi4M0nhSTpqqTc8XIguuxDXrEXJCoMQXUNMsCaGKK5TY1loWaJwa89g+j\nOnSUmIkSMlFDTd53pF8v3swUbNZ4QWMmAexkNw6gaQitznBgSJ2vKpI0hsbM4nkAzhj67bm0yCE2\naFFE3gBWgkBiAAAgAElEQVS8AeCJx/alyypH2cVMFUmbrTfKS1MFMeNR5FNTGgq2NUmDozu25tjx\nnFpCd9r9qCEmMz7CFAg9nZi7PLKqN76ndqP8Qhq8cwR5a8x1C6a3xhVjQFfxRqBTmbrMdH02McJx\nkCFJb/ZTUIHSms+7ZYqalUhkYj2DNcbfwTjj3K+M3CNnVPVjqnqaqp52xFHZCpqyeWeqRJXuky09\nUEpKbWuOWusuNDqwsFiTJCyyZnk2S9Rwh/FkPTG1/JeEIJERFiMzI8LCAPEV29z2TdP7F0eSmUJ5\n0nY/hyBvyDy1HoGVt2gwY2q2yhI3S6vnO+Jvq9dO8y9Lfs6RnXt0Qwln/yZVAjeoapejUUR89Ucz\nJVXQYhh7WJ+60yrrFO1+vDNpYz6yospVu5NgBptbr1nmFGlrMrEz0N902E4Ryoin+9C4DA9j+KYH\nd5kXTxGcl6SXOI9OHDNuNl3T85NE0JQJv1CoBQyFeVO2wecd832eaWJn/PgDhT0PjN9TE0YtxNMV\nNYRWZZIKmnfSa1CClmVFbNBiUvoRNX6q5HUw8aYCD7PTHWZW4SA27Nq7/Ob0wY9Xltimf+WVrMJ7\nYnzPUNsyIEXamtR2Zv8RBwCPQQOmZwfv/Gksd0iml8brhGKHEOrBi/3bjM8tMdFwA0UDgoC9LMNR\nwcHNWZZT6fu282JEvKGq0JpHJvEJkXsYJDg4SQK/sOngQQKnS9wYM7c64sa9V1kktINlcTPplj1L\ncg+CtjHbPmrCJq7a9stwotWOEZGPGKvW4lTFzIWwoMW8zhdFVcWMybBFjUmpij3uEjg9mzwijrDJ\nLI/MimMYtqZfO+PlHJnw+Y368dL4hUpQXEQSYmcK1btfb/Yl2fMH+/qDnj0Pi5OEr8Y695xb3YBS\nr9MOy2/SJWb6EDJ+vCnuWZM0YV/UuT1xM85S555mWXZgnH1sxBE23nBlmkSD5owu57vSR+R4ChJU\n2/4j4O2AAPuAc1X1OmP9KuBq4DZV/d2oc8V5aG53D/QKnIygHvuAtya6mj4JClocBkEdb5VETlmE\njJ9SCRtLGRiKrUlrZ+7lMOCRzvtpBvfS+DvnxHWS+hA+XcfXXmFj4h828h+jTY2atrlZljtEv5gJ\nK96YRfxJWlFTWz0/UJmFqHIUvSLBETbTOPc2Ky+NiTddP+4ehFVFj43byoCEqRFuAn5DVe92H2w+\nRndw/p8De3EebiKJFDSuSrpORD6pqrl5ZKpGXNbYMpJ35eh+Kbou0p7TT+oedsrQS2Ppn6rZms7U\n2IwGwvrxOAzioTCHubbqeUwHiJqe2VQuZhyHI2pwRM1ct5jxvDo76o2uWVtJxExeM5wGFTUmQULB\nL2wmtZ3btOhOeYSE19N+YHIYGZaT1I37b2P7b+LEsuFu/8s4ZVcuAMwZkIHEDTm1VLUB7BGRHquv\nqs+KO8Eo029HXLSgKLsXZASKPVaOIt3ACdtTGVtjTo/1DzulZasssVnbnafteWqc94BTRSlJ59Nv\n52weezMXL+dmxumMdwAbp1o91+bNRNpIi404hSrb1EBr7KDNOvc4G+danezEpvAyEwAOiyxFDQQL\nBc/7YQ7pZDnstJgijiaqre38h5wSpUYweD3w78b7DwFvA/eCY4gbcvpz9//ABsuyTNLOu0yelKrT\n5ZXJiarEzxTtBk5IpWzNPDUm6+1MvDTjjSV3BtISzakW7Zbz1J1nFeqwYQiTxdYYTV+uFXOm1DRw\npltvqFOZuwFn17fTrjvDUaYnqOv8AYkEE7c9AzHkXXeW3powUeMNS2Y17LQgwsQU7Gg5wUhZC7Qo\nfv7YEZ1zPfzwTwDOEJG3GJv4C+EmRkTOwBE0L3Df/y5wp6ruFpEXJzlG3JDTHe7LN6rq230nfy/O\nE9zQWb30UGCHtef0k1IfawM3dDqmYcd59DNMlCQAuGxDZnne30RCJoNhp6qIGZdC3cBJqIqtcWJo\nHDxRQ0hW1n6YmYX1somaNjqemjxJOgwRVD5hwg0A3uwKkoVZYBa26nls5uKe9tdWzweKkWHWaery\nVuQgCto5DDttw5lhZ5JEoOU03PRVVX19xPpEqRFE5FnAPwEvU9WfuYtPB14hIjXgEGCtiFyqqueE\nnSzptO2X0mtQXhawrFRs2LW3L1FTRs5iZyoBkkQkhG1TpNDJS8yk9siUJJamX0/SD9onw+qOs2Rc\nRK42Vvufmgp1A6ek9LamtwN2XPZpZziZM4L8YmGSNuf12T6Ape3dKY3HNi0m2i/QW9KKDnz2hqG8\nKt3TstQ1hOWRxCuUBf1O7Q6dsj2g0NnMxXx2ahMzs8tB01lUaQ+ckl2yKuQkqxt3LPA54DWq+gNv\nuaq+EydlA66H5v+LEjMQH0NzLk5huONF5Hpj1RpgV8ILGiqDJNZLKyLyJqw9/qri53NBV+XmtDEq\naUVG0ntUlKer7+GlIYuaDIfFFlU1k9IEWbiBE56nkramTY31fcQgeB3bNMEzieapDTyUkFTEmMTN\nfPF3xE1VaHQHA2/DFRWr04sB89z9eG7ymLWThUjY0WqAZBurktfU9SxJWDfu3cAvAf8gzu9if7/2\nK85D8ymcJ7O/A95hLN+nqj/v54RlJmjYoGzBqv7hIr+Y8TCHzvwCJ682jQQl8dTkTKFu4IRU1tYM\nMnV7YhZwg2c9j800cIW7vl9R04+Y8dOJE2otOYn3YrwK3jYzImyca9GuV8KDUBibdayTt6efnEXe\nfpC87EUZSFA37k+AP4k5xteAr8WdKy6G5l7gXuDVACLyBBwjNiYiY6p6S9wJhk7CDqpiMRCJvEem\nd8oUOJYVT6Fu4CRU2das0jGmZSlxJ+WVPFhsjeHIju74Cn9nVWTQp4mZrC9SzLScdRPGopkpN5Zm\ndbrpwoN6HMrosTC9RlG5fdKw2BrrO49MGe9RViQqTikiLxeRH+LMfPhP4Ga6x9TLy6XAuf0XWysz\nYd6ZlUwRs5lM9rA+E6FYZLvdPC+eG3gv0PLcwJ4rmG438LW+mJzcqKKtMacgp6lS7Z+63Nm30dvp\n5OHZqK2eDzzuJG02zjnDI6E1pCJoqjLt7lujvaK9Mj0MkCW5E38zNdi09+USnEPPW5s5SYOC/xZ4\nPvAlVd3gjqsP/FRmWcbr0OKCmAepS+Xtaz012VDl+1ikGzglpbY1h3Fvz7I2tcCcLUmYp8Y0ra4h\nqx314F4v7YyctELC70kZZ1+yWlEBTKiyyJrAYoplJGnnPsh1zFNjY585X0yh7HlnsqBGO1WLHls4\naDng/KrHZ9KGLEkqaB5R1Z+JyAEicoCqflVEPpRry7Lg0mE3oJskT+H+bUZlllYRDOzl8L4vCYpW\nZilmivYqlZzS25rAYMwGiQXNNE6n5LFKx0DcOJVZODtBZ5Wn12N8LpspxuPsg4awsREcT1NF+h3i\n8fZzsiSnu7+mmGnOwtlcnGr/qAKUw0xwmAdJBc09IjIGfB34pIjcCdyfX7OKJY8n7aw6qaSem5VM\nJvfaEzOXJAjU2yVsYPDPxAqZQCpha/yiJmlHNSNCcxZ2+BLONWdbLNbHUndWZWfGzU0TNI27TORZ\n18j7rvRzfHM23BXa6BG7g7Z7RVXbNngl8BBOkbg/Ag4Dzs+rUZlxDrl7aYrqlDrnSeA9GOj4JSFM\nLOTSzhy+I4Xfzx3AQcWeMidKbWvMISdT1GzmYj47u6mrQnUYi/XuIYM2tb49GIN0ZkXE6VSJnvuR\ncWffNkpmpGWVhg8zxYmaTrHMERMvQSQSNKpqPiHtyKkt+XKuJHv6DqBsnX2WlPXaCmvXgGKmrPev\nqlTB1kwaFZVNUeN5aUJnO811e2f6IctOtsg8JlXImeLHLDSZFWZtp6QkGWZKUkpiWMImQd24E4H/\nB5wCvEtVP+AufzLwCZwJdIqTIPTDUeeKS6y3j2BnoQCqqlnVccmPGC+Nf7ip9B1U0DT0Xb6ZFQnz\nqJT+WvOmHzFj3OsVf/8ypGq2JkzURDEz5eQi6bdjzyVpXFqhYVbUbtGZsj3KZCFsvPs8SZukBS0W\n3On9YUHiYcSJmyKFTcK6cT8H3gz8nm/3/cBfqOo1IrIG2C0iX/Tt20XktG1VXaOqawP+1pTNwPSF\n2zlt2LW381c5/GImbJlBZa81SwYUM5ZsqYqtuZfDmKfWE0zpVVZuhg05NZyEeZsHiJEpysPR6UBN\n8TInMCcs1sdcT5RTu2khwTT1aco/yykJWV5DXL4i777uaDUGOu/yBO3eYwR9j3OgUzdOVR8GvLpx\nHVT1TlW9CnjEt/wOVb3Gfb0PJ8XEMVEnSxpDMxqcKz0TQHPv2LOIz/BPWvV3rEHnuNTdxj/MtpI7\n5ZLNehsmRbqBVxKLLSfJnicGFutOdtjmFLxZtw98/MD6PXnEetSdY85zdicnzWJ9bLkDrMPGqZZT\nkLIhLMz2Jt/zAlpX6RijwiDemhrtxDPIJqaynZ4NucUIxVXbTls3LhARWQdsAL4Vtd1oCJqfE95Z\nnUN+WSyK6iC985zje5+EAWKHYlkpAsFMzJj0u1Tie1O0G3iUMKttm5idxaS2MQsPrHNnp6TJmJuG\nNJ1e0lpJ5pN72JDHjlbDGbZoLLENaIYcK4lXyt+5lt2jkyRmJYi4oPEFESZmHTGTt/dk+RpSxPTc\nybJtux2Ir7Y9MO6sx38B3qKq90VtOxqCJopLyU7QDLuTyvv8cccvYNbYUOlX+FXvnnTcwAAi4rmB\nO6JEVe8E7hSRrgxuqnoHcIf7ep+IeG7gFSFoILqzDXPve4nw8hI1SUkqFKK2MzvweWpMttogS05n\n7C9eOQVnJ2zbpOvB8M8CKzt5BQ9X6R5EkKhuXBgichCOmPmkqn4ubvvRFzSQTtRUr3OKx1/6oV9h\nMor3Ji1lvQdXdD24jPtKFQzVDbzS8Wf1DcvyW5Vp0/5pwvPUaM62oOV4GLysx9O4VaYTHnPjnHMM\n6qM5KyouIHiRNc5QU310huhIUDcuDHFKb/8zsFdVP5hkn5EQNI8uwn3bnddrN4VsVNaOaBjYe5EM\nVwjGfrcyxjvfACyq6mmDtyScNG7glYRfrJj1kuLKFeRZgDJrseQfctlRbzBZb7MKaLpxImk8LTXa\nLNbHGE+ZRbeMpB6OargPnK2xrgzSic+XQxmMrFDV/SLi1Y1bBXzcqxvnrt8mIk8ErgbWAo+5MTnr\ngWcBrwG+IyLXuof8K7dsSyAjIWhMgjqDNB1RUOflP2a/HVsGHVUkRXW4KwKfmMmCvD//FBTqBl4p\nxHUs/dRWyoo8vUJeB+55brx6RUliQGq0maTNo7LEo9BXPawyE+S16ZrhNCedgpX9xsxECWb/umEI\nmwR1436KY4P8/BdO2obEjJygCaKfjiRqH3Pd2k3l6agGFXNhx8xCKJXlHkWx9hL3RYiYCboXVbiu\nAAp1A68Usu4saqvnc/XaQLBHKRNaMFmPHjYyxYzXu212Z1JVbbhpUPyB1/1cf5R4KeK7VAZWhKDJ\nk7J3aH7x1c++ST1UZb8XSUkqZqtK0W5gS/8U3RFl8STvDR3FeRwmabNVqj/E1DduXp+sZzOFfXZV\nidEahKEIGhHZAvwpcJe7qGMQReSdwOuBR4E3q+p/FNm2ix5Ktt1bD8m3HXmQVWc8Cp16GPcdujLy\n9BTpBh4mZbY1/eDvlPISO1HDF5H7ed6Vevw5/DlZBi0LUXayqmI+VPY9AlcueG+G2pQghumhuchL\n1uUhIutx3N9PB44GviQiJ6jqo1EHWtDkQiQr0p6vigJoJbH2wXxy9RT9vbQEkpmtGQZRQsJcl0bc\nBB0zTdByFk/7O+oNNuvyVO1RFjNA2hJOlj4o25DTK4HPqOovgJtE5Ec4OTOuHG6zBieoY3vrIVYY\njQJWtFSSkbM1ScRNUnFkkpsXKIU3Z+SYi3ZyVnHaehkYpqB5k4i8Fmec/i9U9W6c/BjfNLa5lZDa\nDSLyBuANQEjuzvLTT0do7mPFTXYkGWoqt3DZCtgvRAiZ2ZrVx/5Szk3NhizjJVZC7EXetKmxsdHq\nmsUVFz/TT6K+MoqgBGVWxF1fAx4ANnk1nOL29ZOboBGRLwFPDFj1LuAS4G9wasH8DfD3wB+nOb6b\nKOxjAEeLjH7J1wCsd6c4yi1mVjZF2pojTztuRdoay+DsqDeYptUpb5AHnghKWvMrbwGUsMzKy4Cn\nun/Pw/nNPi/hvl3kJmhU9SVJthORfwQud98OlB/DEk1Up5yX2Ak7Z9hwW1A7/NsVKcyskCk/Rdqa\nw7g31ZNzGZ+YLcOhTY2NbkZlLyNyXiT9jtZo5x3aE1tmxX3/CVVV4JsicriIPAlYl2DfLoY1y+lJ\nbk0YgLOA77qvLwM+JSIfxAnUeyrw7SE0ccXhddxxgsJcP0hnH7ZvkmNmKTKSCChLdRm2rYl6Ys7j\nPB5JazEl2d6SHTvqDTbSYmEKJhpLTLbauRehzJBH4B6cSYGrcJwmPC+DattB2xyTcN8uhhVD8z4R\neTaOG/hm4M8A3FwYLRwFth/YXMZZB6NMXGc+ap39qF1PUooc1x4ypbA1WRYuzON8Udv3W8077b5Z\nEHYdZRFsbWpM1ttMTC11MgRXBVV9ROR04E3A8bjOzreUqY8eiqBR1ddErLsAuKDA5lgsK4qix7WH\nSda25jDuZdLXaVboCbsvBhFjRQicJO0bttAy2czF1FrxbfZ/z8rBf0/AkxbgDuDXUP2XODGTZGg3\nbJuDEuzbRdmmbVsslr64FuPnvCpm40LHtUedcnY8/ZOnQDOFRV/p/TO614O2Y9BzmlTp+6Oqdzpe\nmh8AVybRD0nKrFwGnOfakucB96rqHSJyV4J9uxgJQXMHLG2BG4fdjowYBxaH3YgMGJXrgOFcy6+k\n2Pav4IpnGe9PEJGrjfdDHdceJa7b/djSUbI0CrYm4jtdVAa49OcJ2WPA32cx15vgLMOymWlsDfDf\na4C1SYaakpRZwclYXgN+hDO8/bqofaPONxKCBrhRVU8bdiOyQESuHoVrGZXrgPJfi6p+k+6cKpb8\nGAlbU/bvdBpG5Vqqch2qugQkruOQoMyKApuT7hvFqAgai8WSnELHtS0Wi6UIDhh2AywWS+F0xrVF\n5HE4Y9OX+ba5DHitODwfd1w74b4Wi8VSOKPioflY/CaVYVSuZVSuA0brWgof1x4xRuW7MCrXAaNz\nLaNyHUNDnOEri8VisVgslupih5wsFovFYrFUHitoLBaLxWKxVJ7KCxoROVNEbhSRH4nIO4bdnihE\n5OMicqeIfNdYdqSIfFFEfuj+P8JY9073um4Ukd8ZTqt7EZEni8hXReQGEfmeiPy5u7yK13KIiHxb\nRK5zr2XGXV65a7HkR5XsDFhbU7ZrsXamIFS1sn84QYn/g1NY4nHAdcD6Ybcror0vAk4Bvmssex/w\nDvf1O4D3uq/Xu9dzMHCce52rhn0NbtueBJzivl6DkzZyfUWvRYAx9/VBwLeA51fxWuxfbt+RStkZ\nt83W1pToWqydKeav6h6aTgp3VX0Y8NKwlxJV/Trwc9/iVwI73Nc7gN8zln9GVX+hqjfhzDZ5biEN\njUFV71C3UKGq7gP24mSQreK1qDqJosAxNAfhFDKs3LVYcqNSdgasrSnbtVg7UwxVFzRh6dmrxIQ6\n+T0AfgpMuK8rcW0isg7YgPPEUclrEZFVInItcCfwRVWt7LVYcmFUPvNKf6erbmusncmfqguakUId\nX2Nl5tGLyBjwLzgl5O8z11XpWlT1UVV9Nk7W2+eKyDN86ytzLRZLEqr2nR4FW2PtTP5UXdAkSeFe\ndhbcKsa4/+90l5f62kTkIBwD80lV/Zy7uJLX4qGq9wBfBc6k4tdiyZRR+cwr+Z0eNVtj7Ux+VF3Q\njEIa9suAje7rjcC/GctfJSIHi1M+/anAt4fQvh5ERIB/Bvaq6geNVVW8lqNE5HD39aHAS4HvU8Fr\nseTGKNgZqOB3elRsjbUzBTHsqORB/3DSs/8AJwr8XcNuT0xbPw3cATyCMyb6euCXgC8DPwS+BBxp\nbP8u97puBF427PYb7XoBjmv0euBa969W0Wt5FrDHvZbvAu92l1fuWuxfrt+TytgZt73W1pToWqyd\nKebPlj6wWCwWi8VSeao+5GSxWCwWi8ViBY3FYrFYLJbqYwWNxWKxWCyWymMFjcVisVgslspjBY3F\nYrFYLJbKYwXNCCIiS/FbpT7mK7wqwyLyeyKyvo9jfE1ETsu6bRaLZThYW2MpE1bQWBKhqpep6oXu\n29/DqQZrsVgsmWJtjaVfrKAZYcTh/SLyXRH5joic7S5/sfsEMyci3xeRT7oZORGRmrtst4h8REQu\nd5dvEpGLReTXgVcA7xeRa0XkKebTkIiMi8jN7utDReQzIrJXRHYChxpt+20RuVJErhGRWbdWi8Vi\nqSDW1ljKwIHDboAlV34feDZwMjAOXCUiX3fXbQCeDtwO7AJOF5GrgY8CL1LVm0Tk0/4Dqup/i8hl\nwOWqOgfg2qcgzgUeUNWTRORZwDXu9uPA/wVeoqr3i8jbgf8DnJ/FRVsslsKxtsYydKygGW1eAHxa\nVR/FKYL2n8BzgPuAb6vqrQDilLRfBywBP1bVm9z9Pw28YYDzvwj4CICqXi8i17vLn4/jRt7lGqjH\nAVcOcB6LxTJcrK2xDB0raFYuvzBeP8pg34X9LA9fHpJgewG+qKqvHuCcFoulGlhbYykEG0Mz2nwD\nOFtEVonIUThPMVEVW28EjheRde77s0O22wesMd7fDJzqvq4by78O/CGAiDwDp0AbwDdx3M6/6q57\nvIickOB6LBZLObG2xjJ0rKAZbXbiVHe9DvgK8DZV/WnYxqr6IPBG4AoR2Y1jTO4N2PQzwF+KyB4R\neQrwAeBcEdmDM37ucQkwJiJ7ccasd7vnuQvYBHzadQ1fCZw4yIVaLJahYm2NZejYatuWLkRkTFWX\n3JkIW4EfqupFw26XxWIZLaytsWSN9dBY/PypG7j3PeAwnJkIFovFkjXW1lgyxXpoLBaLxWKxVB7r\nobFYLBaLxVJ5rKCxWCwWi8VSeaygsVgsFovFUnmsoLFYLBaLxVJ5rKCxWCwWi8VSeaygsVgsFovF\nUnmsoLFYLBaLxVJ5rKCxWCwWi8VSeaygsVgsFovFUnmsoLFYLBaLxVJ5rKCxWCwWi8VSeaygsVgs\nFovFUnmsoLFEIiKfF5El9+8rRe1rsViqRxV/81VssyUYK2hKjoj8rog8PsF2TxORa0Vkn4i8Oavz\nq+rLVXXM/fvNLPcVkb8Tkbdk1daA439bRJ6e1/EtlrJRdntRxt98FdtsCcYKmiEhIhuNp4KHRORR\n4/09InKwu+mHgAMTHPJtwFdVdY2qfiSjNr5ARO4VEfEt/7KIvD3B/mtE5D0i8iPXcN4kIheLyFEi\nchTwWuCjWbQ1hA8A5+d4fIulEEbBXpTxNx9lo9z1pWuzJRwraIaEqu7wngqA9wCXG08Jh6vqL0Tk\nROAbqnpvgkP+CvC9ftoiImEGcANwraqqb/mzgWtijnk48A3gROBlqroGeCFwkNvWTUBbVR/sp80J\nuQw4Q0SemOM5LJbcGRF7sYkS/eYT2CgoWZst0VhBUw6eDVwXsPw3SfBk4I77ngFc7D6xneAuP0lE\nvuY+wX1PRF5h7HOz+9R0PXB/iJHaAOzxnWsdcKR/eQAXAT8H6qr6QwBVvVVV/0xVrwZeBvyn79jv\nE5F/Nd6/3326e1zIdY+5T6pPMpY9Q0TuEJE1qvoQsBv4nZi2WixVoqr2omy/+Tgbhb/NadubQ5st\nEVhBUw6eDVwbsHyNqn4zbmd33PcbwHnuE9sPROQg4PPAF4AnAG8CPikiTzN2fTUwCRyuqvtD2uX3\nxGwAblXVxbD2iMiTgdcA71LVx0I2eyZwo2/Ze3GeVDaIyDRwJvD7qvpw0AFUdQn4PnCKsfhC4D2q\nus99vxc4OaytFksFqaq9KM1vPqGNIqDNqdqbZZst8VhBM2REZC2wjgADparvHeDQzwfGgAtV9WFV\n/QpwOY5R8viIqv4kyJ3qGrin0+uJOYWY4SbgJcBdqnplxDaHA/vMBar6M5ynph3AO4FaAvf5VW6b\nEJEXAevpfkrd557LYqk8FbcXZfrNJ7FRPW3us71ZtdkSgxU0w+dknC/zTRkf92jgJ76nj/8FjjHe\n/yRi/5Pc/3t9y5+L8+OMYgK4JWabu4E1Acv34DwVvVNVo9rn0TEUwPuAv/Y9La0B7klwHIulClTZ\nXpTpN5/ERkFwm9O2F6ydKgQraIbPs4HrAwLpBuV24MkiYn7GxwK3Ge+jzjkBPGi6lkVkHPgNnCC2\nKG4BjvGd28/1wAnmAhF5JnAJztPPH8ecw+Mq4BQR+QPgEOBTvvUnERxvYLFUkSrbizL95pPYqJ42\n99nerNpsicEKmuHTMx4uIgeJyCdE5CsicqWITPRx3G8BDwBvc4/3YuDlwGcS7v8d4GAR+VMROdQN\nHPw0zuyK6912bheR7QH7Xu7+v9B1kSMiJ4jIh10jB9DGMXa464/BGcOfBt4IPNNts7c+7FzXAU8E\n/h7niUmNfQ4BTgW+mPCaLZayU1l7wZB+8yHHSWKjuto8QHtTt9nSH1bQDJ+T6R0Pfw6wzw3e+3VV\nXUh7UNed+XKcKP1F4B+A16rq9xPu/1OgAbwZZybAFSxPvfR4MrArYN/7cGZcnAD8UETuAXYC9xvB\nxJ8Aaq7xW4tjOD6oqpep6gPA+4ELEpzrFzjG9GZV/Xff6pcDX1PV25Ncs8VSAapsL4b1m+85TkIb\nZbb5sH7b22ebLX0g2XsuLYMiTpKsGeA44FJV/fyQm9SDO03xOuBZqvpIn8d4D3Cnqn6o33O5634E\nNPwzPETkW8DrVfW7/bTPYqkCVbAXHkX/5ouyU3HnsXaqGKygKSEislpVHxCRI4HPq+rpw25TWRGR\nC4DjVfXVsRtbLCPISrMXVfzNV7HNVSRJimxL8XzczZNwIPDXw25MGRGRU4Cv4gTtnTXk5lgsw2RF\n2PoJYkIAACAASURBVIsq/uar2OYqYz00FovFYrFYKo8NCrZYLBaLxVJ5RmLIafxw0XVeWa+fA8ee\nCrfsdv4Dj16zO3TfVaecmn8Dh8jtu51rn5CYDQ1WjQcsPDJg2c8THMz7LPzH+rnv/Vi2n0PUZ25i\nfv5B+3Tuhf/6/ddurH9g7JBE5w5i9dJDAOy+kUVVPSrpfiJy1Fq48z54iqr+uO8GWCI5eHyNPn7d\nOMf/+GZnwfGncvvu3Rx9qvM98n5vTzj1AA68+zEnLRtw+910thlVvGsP4ujjgbvhkbvhLmPZ/iOW\nn6nv5TDu5bDA/Q/j3p733v195G446NRTe87vHf/AHz/G7XcbyzP+HB5xz3vQ8TEbHmGc98dGW49w\n2mle/72PHcZhBzjXfJi7BuAWju3sdu9jwfcqCd6xAe7efXMqW/NyEf0h8H3VFL1KMYzEkNNpJ4pe\n/c/um0t9Ky9R7ju0976vfbD61x1H0HV7rN0E920PXt7FOQlP5r/vUUQd8/TsPpeo6/dYuwm4ZPmc\ngd+VTQS32btm37o9p5/Us2kcG3Z1J1iVF7BbVU9Luv8LRPQ+nCeUaxIYGhE5E/gwsAr4J1W9MGCb\nFwMfwqk+vKiqv5F031HlyNOO05devYWtnMf43BILU7ANJzHJhCoLIkxM4Uxgbrk7tUbb1sxI+Ndt\nGpz7AczMdq/brGPMU+u8bxuvw6jR7ryepM14Y8l5Y95v9/2OeqNru4VZ57PyaGbU9y2ILH8HjGtc\nrI8BMD7ntrHF8nehISzMLt8bgMWWcz+8+9B+YBKA2ur5zjbesiCWtvc+iY5t6i25V1vtnGWSNvPU\n2CStxLZGRA5+Ajy0FngR8M8lszUj4aGJ5Fzp6qSDOvFRI6wj94uVrveXKJxr7Od20mbn7O90uzA7\n9SBxk1QY7ZJMRU0UPeKNbqHbdR8vBc7xiZWIuSQ7I+L/zmJn8kYmQESOeiLwp8B25/3xUV4aEVkF\nbAVeCtwKXCUil6nqDcY2h+PkIjlTVW8RkSck3XclsJmLqdXbTGqbptepzjliZrE1ttzRriAx0/SE\nnIvXoS8Cj4pzP5qqzIjQnIIdXsdtipoHJrs6cD89oqcFG+dcJdNwzjnvO26bGrQcMbSZNo/KUpew\nGZQJVfDugyGqxlnqamd3u5Wu7IcNYXxuicl6u+ca04oY//ogUTPpCsNJQyAmoQYPPQiciFPFNI6i\nbc3oCRqv4wzxGPifyFcKXZ33OSzfH+9eeGLGEB5+T0OY56FH6ESIlyjvRec4GYkavxcuUOil/C54\nbQy7jighE7RNRuLmL34N5xHmRcAn4a+AP4nY/rnAjzzRIyKfAV4JmIbiD4HPqeotAKp6Z4p9Vwzj\nc0udjhSA+vK6hVnoJ2VvFZk2XnfuBTBPrdNpemIGHNEHvWLG/O/HL3Ta1KjR7nhi4jw8bdf7UdM2\nzUbL8aZl5KUxvT0L7jV2PntP5EWI24VZYBYmZpeouaKmtnp+IDHj384TNjXXMwPpBI2IrJoAXodT\nu+FgZ9nzVPVbEbsVamuGKmhE5GacQmuPAvtV9TQ3l8JncSrK3oyTiOjusGOkYgUKmS7O8b2+FEc8\n+IZO/J211wGHdb7e9mEenKRDMF3HycFTs/bB7uHH+7bD2kuitweWxZ7rpfFIIl7i2MlZnMVO9px+\nUrQHLJrnPs198avOv5eKyNXG+o+p6seM98fQXWjwVuB5vmOeABwkIl/DKZz3YVX9RMJ9S0dWtuZY\nbum468cbS52hjnGWWKyPGcMLumLEDLhDObPQBObr3cJinhobp5ZdFNPusja1yA7bT5D3JkgQxbIa\naMGZ0uoIrKyGn4COSDKFTZy47QirOWGSXi9NUvGShrTeGeDgdcBb3DcCNOENIrLV2GaotqYMHpoz\nfKmm3wF8WVUvFJF3uO/fnvqoK1i8+IdN/EMrnc7T9NSEkKbT9o4bJ2CCjpn1MEwYax/0Da0lwT8c\nF8Cc+WgO1Jljjjp15tI2sS82Hwxr3SbOPMRimvibEA7EqS/zW8ChwJUi8s3oXUpPJrZm41wLWm4n\n1dKOV2Yc8H0NRh5PCHjCYMGIIfE65RptR/jNCc1ZJ7YlSbxMEGkEUOQxVuOU2hT/WFB2eDFVC7Px\n25qMzy17abLGLwjH2YcjTfrmq6r6+kEOQIa2pgyCxs8rgRe7r3cAXyONoCko/qIKRAXF9ngE/PEh\n9OeBiBIzSY83oLcinjiBEiNeknicPIFjCp2ixE0CbsOpO+Pxy3RXVQbnaehnqno/cL+IfB2njtCt\nCfatCqltzYF3P+bEQzRcMWMJpMtzQg3qdIJQ/UGvw6JNjc9OtXqClbOkEyge5gGaK3aiUEfMdTg7\n71MWamuGnYdGgS+JyG4ReYO7bEJV73Bf/5QQT52IvEFErhaRq++6B0fIWDETS09nfA6JxUwWQyxh\nBB57lzh/eXCJ9u3F84stv3fGT525RGKmn9lRfXIV8FQROc6tMfMq4DLfNv8GvEBEDhSR1Tiu3r0J\n9y0jmdiahf3ixELUra3xMAODzVk7Jm1qbObizjDT0MWM2wYvlidqptagRMbptJyYI+8Plmdn5YV5\n7/sYdkpLobZm2B6aF6jqbW5U8xdFpKuyq6qqiAR+G9xxuo8BnHbaaf1bF3+HeSmjMVx1bu9Qk5+g\nDjROtKQNaE0jgrx4kqETNOPLGJorUHjkgqruF5HzgP/AiSX+uKp+T0Sm3fXbVHWviFyBk7L9MZwp\nk17Bv559h3Ih6Ri6rTHjNfKI3Rg2/llOJkULGH/MSdBMH3Bmq23V89gqS87sq6I/j4YzxNSZYl5f\nnk5tMrZpcaA4mrCZTp3ZeDlRtK0ZqqBR1dvc/3eKyE6cqOYFEXmSqt4hIk8C7ow8SBTnitMZWc9N\n7LTpfrwvccHCg3h0ch92Sso5vv8uSa8tzTBTkderqm3ofjxT1W2+9+8H3p9k37KTt62JFSgNZ/00\n+XoDiqZz3a6YMWc4eaQRM0lzqSTd11wXJWo+O7Up16GnUNwQnsl6m81c3FmcVVBwWB4aj7TxPf1Q\npK0ZmqARkccDB6jqPvf1bwPn47iUNgIXuv//LfFBg2IfLgUulWCviyd4TM6h0FwomeO/BxFCJqsh\npCyHonq8NJcSme8lU8zvQ8j0df+1BgUDp+Esdi4LmTTJCS2JydrWRAmSsKf8mVlozgItmHY7kW0R\n21cB7z5M4+bemVtyZnn1ERgdJ0Y8gjroQTwXnZlTDaAgQdOVu8f9Tnic98DW3h36IEy81VbP89m5\nTZ2g9izz8ZSBYXpoJoCd4ny4BwKfUtUrROQqoCUirwf+l1Anpo+4mSveen/eFX8n4nVkVRQ1QxAz\nebFh1164NH5qdWb47l3aAOl+An47ws0KmbzJ1tbE4BcpM0bStaBOpGqixuyQp4FV2uuV8fB7Zwad\nfpzH9OWi8bInL7bGnOSCdSKnsCe55iRerNrqebZyHjMhcU6jwNAEjZss5+SA5T/Dmb6VnFt2w3NS\nbB8lfszOpSjPwKCECBn/sI3XSScRM0kCXUcGL47qnHT3CPq/D56Y8T6f+7b7EgFuSzcssfaPYO3j\nvH37atLIkqWt8eoFNaein3A9kWJ2/kMZ0sgYv3fBy8wbFFxqdtDDFiKxHX5+s7e7mBFhGtjRavSI\nmH7uUdLhOG+YyYuZMQX0lhRDoAeNw8QZzus1PwCuS7xrIQw7KLh43M7fK4FQtczBobWGPHzJ8aI8\nDXGiJYq8piObw03eZ3TfoZJf7a0IMRN1f5Jes3k9YTFHK6EcxygyMUXkMIV/aKpKXhiIbr+XTwaW\ng0vNWkReR112IRNVZiFrvPt5hTYGukdJRYyfWrXC3vpi5QkaPxUSM7GEZPqF7IRMELkkkCt4GGbP\n6SclEjKQTMwEBUrHzeDKVbhZMsNL878w67xeKQ4xU+B4yfFqbm2kGaDZWKLdGr6Y6avDd2No8hz+\n84bnNrtiJs396UeYRdXFCvIgjgIjIWgeTfn9NZ+IezLpllzgeOn7w6pix00pzlrM5MY5sJZivDRR\nYqYfT8wgJKkQbhkuXpr/MJpTy8NLpuCpqsfGqyLtxV5MA1dAJ7gUz1PVIHVH3S/9ein8eJ2+423K\nb9zJq8J+NhenukeDeJjivE9FiZmkFbNF5DnAlcCrVHXOXfZWnLp0CnwHeJ2qPhR2rpEQND1cslyz\nJy4XC7jxCwm2KwudtqbIj5KlkMkzfqarthF7u0RN1uxhPRhiBgYM7k2JFz+zdlPxw05FGplRxnvK\n9bw2E+705SZuaYRZaLaqGT/jXdvErDcd24m/mJiCjY3uDLvNKXcI6oH825WVmPF3+J4QzcNLsw3Y\n3BoLDPxNcz0DD5EVOJsLklfMdrd7L0YRbxE5BngzsF5VHxSRFk5yve1h5xsJQbPqlFPhkqsD1yUV\nK52YmiJm1GSBb3rxINWfwyhb4G9HpGbgqdnDeqD/+5NnAsC1Dyrk+PRUtJEZJY4+9VSaV/famm24\nQ1H+eVJumYQm1RQ1ALScgpPTtNjG8hRts2Pc0WrkPuU4C8IEQZsak602iCPashI15tR2M89M2mDe\nrFgwZjh515cmKLgPklbMfhPwL/RO7zkQOFREHsEp2nB71MlGQtAEYVZWHrmgy4TDS4N6HkqBb+gp\nCzwxE0aQWMk6i7GZe8a8toJiaAo1MqOO58nYBky7HYYncBZmKfSJOE/erNvdwNI2O+oN1rlf1TY1\nPmB4HvIUJP2SRBjMU6Np1HbKStQ0p+Ds1vYu70y/QiUusNefkM+/vTf8maEH6gwReYvxPnW1bfch\n6SzgDAxb42b2/gBwC/Ag8AVV/QIRjKyggeXOwRM2Fz20/EHed6hwkeskf+shQ2le/5yuA3sY8iSq\n40/S3q5hJ7cq+Fp3XT/CxhQwG3r67GTeln7ETNrMvxmKmXERMd0IQzUyKwFT1HhMRARdViV+pukW\nV+zMylntKzw55JpMQQzi1VhsjTEtS30HepseGa+Gk1cvKk27+p2R5O3nBWybTNJmYB/aOMteyMuB\n6zKptv0h4O2q+pgYvxcROQLnQes44B5gVkTOUdXQKSMjLWg8TG+Nucxz6yfqSPyJ+YZInIchiFLU\nSCJ4GnMUpqjhUjfeJCZAOPL+7BI2sOzdysXrkgT3J5mZ5+ls4PHu620squppAx4xMyOzUvBESleu\nFt+yOCETtO8w8QuyMAFT5PTnPJmnxqQmH3paMETsNG4JCNzhxQAxO8yp0+NzvTloCiBJte3TgM+4\ndmYcqInIfuAg4CZVvQtARD4H/DoRc2BXhKABR7Q0fcvMDzZxfMa5IWUUhkCUICiLgIniLHaGXoM5\nxOMXNRAtajbs2ltIAcm+ay8Zw0wXPeR4CAuerl2okVlpBHUYQZmDk4ibMoia5hSsCRAyoyJi/MxT\noznbSpRRd2IKmkYNK2+wrdlYYmZ2eYp/kULGP+wEblXt1rLgKpBOxWwcG/Mq4A/NDVT1OO+1iGwH\nLlfVfxWR5wHPdytwP4iTBDM4WNZlxQiapCSaHmzW/DHNeIGFMM1pxmZ8TBWETFLMRHRm1uNOkPcQ\nvWZZFJL0hjyHkHumUCNjCSZIsPRkFx7mkJVbUHNHq3f20qiKGY/F+hjNqaXOrCeALWxBdfmx2Fu+\nrtVbMWOy1abpZuUtMqIoSMyA651pAPVibU2SatsR+35LROaAa4D9wB7cqvdhWEHjI/H07aDn0UsB\n8q0B9W7eBaMQ7NsvvqnqPfWqdnV3AEGlH/xEiRNzn0yqYfu+N8N4Ai/ayFiypwjvzcwsbNYxjvPN\nXqqimEkyVGbGncz7Zj2BI2jo8fP37mvuP0jBzjQECRmvTZ2yFAWLGY8k1baN5Zt875uE3fQADkiy\nkYicICJfFpHvuu+fJSL/N+lJKkdQrSfTCzAEj8C7eZcrZsIpi3dmw669nb9+mKPelTfH89R0xEWM\niIliz+kndY6TpI39XoN3HvN8HmWYdaeqbVU9QVWfoqoXuMu2BRkaVd3k5aBx3zdV9URVfYaqvkZV\nf5FVu1aarQnywIQJFS9Go+kOZeSRGG1GxDluQ2hOOVONqyhgPNoPTEYGLscFNfuHaILueVAdK3C8\nIjNTTkLCsG0Goe2VmQjxypjtWCkk9dD8I/CXwEcBVPV6EfkU8Ld5Naxo1j6o8RW7kwgZz0sDgZ6a\noIDVoJk3QKyAyZN+41AG9WKYQsYcTgua+dSFWx09qs3mNfXTTn+xT/+6tFRudl0xjLytSZJyPijA\n2Evat1gfY2JqqZOq39zeJI1Y8lfQXmyNdeVNKSv+2JR+ZmCZJQJ6ZgcFJKITmekaenpUlqhpu5Ns\nsDm1PLPJY3xuicl6m/kY8RHZzn73bQGt4cdiFUEiDw2wWlW/7Vu2P+vGlIqop35P+Pg71RjCZt8E\nLR+mmCkTkZ6aADZwQ6hAzALT4xLmgYmjSlmph8CKszVJPC3NWaeT7ASftsY6AadBiMx0Xkdt553f\n80Q0Z518M8c9cFPHu+EXBmWbqr1xrsXGuRaTrr8C0rfRvE6/cAi7f5443Oa2oZO/ZtYRMDvqDTbr\nGM1ZJ6HdeGOJSXMIKKgdhtfF/5fmGkwWMsyD9Mj4ASzWne/h/Sc/LrsDZ0RSD82iiDwFJ9U5IlIH\n7sitVcPiEo0fvojz4iQUOVknasviWFnEiBQxu8h/Hq/dnjA0BU1Ye6K8LWnPn3R783xZlzu47vkn\nMLZ2lfsug1if4bEibE0SL423frOOsSOgQ9s423LqKAXst8V9vU4bXAFMSyswt4p3jh2tBs25llN0\n0u0Ua6vnu7wVZibgqMKHWRLUQfvjXibry3EvSb0YS9vH+04AuIUtzMgWN0jY+Rz9M6K84p0Ak/W2\n41FzY4eDvDRpvS9lE5VlIamg2YwT+HeiiNwG3ERq/0QF2eUL8I0SMzF3IygRnilq3s27OJ8L+m5q\nFgLJ7HSLEiZ+vCGmsNpTQdfpb6s/kV4/eXtywwhqtp6aQFakrfEH+ZoJ2sKGKRbrY0wTngRunS5X\nxJ6YBaZ6h6g8YXWmtNihja6OtUabjXMtFutjwZ1wDqImSUftz7jrDY0l2dcsChkkaoKuaWIKpmd7\nq6p7wsbDFKkbGy1oLQuVyVb3cFNeAsbfdi+ZXuKo2oqTSNC4KdJfIiKPBw5Q1X35NquE+MWMr5aS\nH6eTXc8GbuiZmWRiChxT1JzPBamHncxj9StuihQyUXl0omZvpRFve1gfun0SL02W96Mnn05J8hmV\niZVua/yxLKt07P9n793j7Lrquv/3txek00nT0ilDRWoLcgtIaalAKCJoUXrGcnmeyQQ1mEEQhkwQ\n5fF5gJ/KNCoIwsPtlSkBfTCBcMnJSLUwhypFEKUB2hKK0FpEWkuhTDulhJwE7G39/tiXs/Y6a++9\n9j57n8tkvV+veWVmn31Z+5yTtT77e7Uufg2CBXLrpmYc46EfO6/GmQ+Pa9EIs2zs3aSjhXhhqUlr\nOjimMbbMFC1WNsHxGV9T3ZpThl4tDUWOt3W4TrPU6NaVCdpMQm4LC5EdsWUsKKzXBOXuMoqvXfI9\nMT+DBq2hCAjOa4QrIi8A/gx4gMC9/PtKqX8NXzsV+GvgiQRW299RSh1Iu1amoBGR16ZsB0Ap9Q63\nWxphbFYZrZeSbUEclHVDpwpxUxdVtGvQa9TkEe1TtZsvbUy2a0ekBjUf4/i5pjuW5njlHpjbLWa6\njztbzfC80PVUNO27RYPG2HLqYpvWp2jQrhGbiNHJcjslRE0zXRjMEVhvLo3ljHGenPegjvdrkBWJ\ndRwb4X4GuEIppUTkSQTK+3Hha+8GrlRKTYvIgwh6x6WSZ6FZF/77WIJ+LleEf18CmIF7a4PIxZQT\nS6MHhqbVOrmcF3GecVyaKyWySGzgIAA3dB1ZnjQBUXRx71ffqLz3SKeIsNH3j44xxadL3Zq885rb\n011kG2oNYM6in09Njhxzc40tkyl+bT9x3IxpBcnLwlmO+i6FNMaWY9fT8Wo8UdZfH8fq9HjiuKgN\nwLKe+TOWv+jWLWLyRIoraWLGzHqCsFDeUpuVTUnXU17AtU6aFatBq3wGk+Xc0WcVVQceMLmNcJVS\nulo8mU783HrgWcBsuN89wD1ZF8sUNEqpHeGJPw+cH5l/ReRSYHSLExTBrAisPVnHi9iF3YdFr+lu\no7SF2vbaBg4mtkSvp7li8l63MSyNLbPel6z99HstKmyiY9IsKGlipsx7lmUtOzgAUdPvpyYX/FwT\nEGUbRUGleeIgyHoK1oM5gniPzcbiGJ8j/JTmVYv7tQaMO0RSY3WWabB1qRmPCci01FRNVeJFp2gw\nsG6pOV7BguHGWcloJ7BvZpbV5jgTYSuEuK7NDHEg8S3ABpqcrWacBGMWupiZWGrHGU4DbKOR2wgX\nQEReBPwF8FAgegPOAe4E/kZEzgWuA16jlDqSdjHXoOBJksronnDb2sWsIbO38xSV9+RuW/RcF+0I\nfbEuIoSy9u21qnDRexgEZS02UezSQTZ0WdmqHt8QuAD7+tRUkGNurtEXmlXWMbHUThUzNlfCwv5g\ncdyFvQx/RNQte4oWk5tgvjnOYmituVLNJLKYIqI4mkD8NGmojvspi14FzyCFjO3eEhaU6WShvMn9\nbRYiS4gmVCCIpZmTNnESfaIeTNIy9zxpBv/LCooaWyDwxFIbmklrUhWi5m4ewjK/DMAN3ARc/xwR\n+X1tl/crpQpXDldKXQ5cLiLPIrAMX0SgT84HXh1WKH838HrgT9LO4ypoPgh8WUSimfiFwJ6igx5Z\nvtCpOxMtdlkuBp00EXDjFR2X0uOff7Dr9bTjehEVoyBIimD2sdIxhY3+uZjvg+7eM60llzsGZtve\n27yx6dfoJcPNwoSI6P2VzEmmr09NBTmm55rlMIA3b0HTF9TV6XEW9of6c6kZHE/68feHC+yctJlX\nQV2bKIA4WhwX2Q6E1W6N67q4R7IEj8tinSY+soRO2TTsiFyRZrRGgOD9WJ0eh+lORd4opiYiLQvN\nzGrbRZBeH4kaGHwMkgOfVUq9LON1l0a4MUqpz4vII0VkgmBeuk0p9aXw5SUCQZOKa5bTm0TkU8Av\nhpteqpTqXoUHxPf4buzaWWK60viTmJQgzmghy+ocbaKLmehvm6jx5JPlhgJ3gVnV9dNey7OORRlu\nZQsqfpJL+Cl+KvzrxlWl1AWlTqRR1VNTwWsO91xz3XWJmJc6zPiJ1GnLImsrzBa5gyL3kElUU+YW\nrR7NLgBpc4Pa3bVPt60mYGKpTWM6P+A0S/TE8SklFuteRYuNIhlaWfc1NR1YRqI074jo/b6US1PT\np6NMs+AzafIemkFQeI8uKFNc6XFTdbTOsJDbCFdEfg74z9C9fT7wU8Bd4d/fEZHHKqVuImiEm+mf\ndxI0InIWQdPQy/VtSqlbC9xY36g0sNYSHGwTLr26JsqKmmmWhtLy0q+mmUVcbmXQLSc2seH6/rsE\nOvexOnRfn5qKMApzjb5IVBmbsIfN8e+2RTbRaFAjsha0aNCa7nZVRUJl60wzYW3RzwuBUNJja6xj\nDIVTWrXbaCx5WTZZbqsqrBL96j/VJXCmw1gbfVvo1stqbgnJIPFIbDZUK9MFZQsG1pncFFxxZb9d\n2NSNYyPc/wn8tojcC/wY2KxU/J/q1cCHw1i9bwMvzbqeq8tpmdCHDpxEYHa+CXiC850VJC8Lw0bd\nC7sZO1PkeqZVJm0fF1EzDB22XcdQR7yILh5dY42qIM2CklcM0NxvCOjrU1NB+jrXFJ1nHvqU47jy\n2mkWZjql7s3ePmVYZR1TkEi3NheoODZCP256nCk6RdtsYsaFKVpxPE0aC5tgVRuTrcbJVI71Jkvw\ntAy3V7/JEmFOrQdS9ymWaqRbTyIXny0QO+t9Wp0eZ4LO5zM5AwtN4u8sdAKUL62wNYKNvG7bSqm3\nAm9NOfargLO12dXl9PP63+EEt831IkVxzMKony90BwLXnRlkihpdCKnnB9YnM3akaitNkYU3T7BU\n0U6h65zYU6tt8TImdYqetPdt2NLj+/3UVHBsfZtreplnVpvjLISZK1WRSI82iC0i+to4E4iKqL8T\ndDKQzCf3vG7PuphZUCosEndpYp+V/YC0g6rD0ViMGOS8Qm5baSbGW4ZMd1ZF9VfM96tIt2w9Wyyq\n1FuUSNRMzLRpNFuZFi39njPHOQML2ucVfw45QnaUcLXQJFBKfUVEuoIIKyQ3C8NGlQt7kE6bT9bi\nX9VY1PODfyNXWnTuqp74qxQwYBExe+37VcF5lr5F5jZdjJrWlKJuSd1KU0SouAq7gxc+nhdxOX9W\naFTF6edTUy/UPNeUmmfiRcNSA6YUSwLGIm9bmLrEQigook7OkHTldC3uKePVF93I7WErEhd9Oeb0\nNOX9nS7gXdgMEzN20RNZmkxsqeRViBZXgeJUade4z60zHdFmCsWi7NgP8xmB2ImaMxou4zbF8FrA\nqdu2iLxW+/lDEfkI8L0ax2XLwni4MaZXiMi1InKtunM1XlymWapkoXddgOp44jczoHQhUwfDGIPj\nTAGx9CIuT3xPyhK9X3mWlPO+cGP840odFq1Ros9zTe48E44pnmvuulNblMKFLK1CrAurrLP2Skrr\n35RH1mIfZEN1b89zNZnsMn52hGnjZoNGZiw/KeOyUfY9cMH13FFnaWfCe5xYaidERRkxEx1zv7R5\nj8yyb2bWWcy5jrloa4Qf8JC4lcMNw9QjL8TVQrNO+/0+Aj/331Y/HHfCNNT3A2y44CQFyUyjaZbY\nQPEn8INs6FpUbO6mfouAoNBexwKVtyAXDTLOyxaKcGmp0NUSQmvI2G9sNYOicZdNlQ6+U8H36mDK\ne+AkTPZizZ47xkXNUM81FzxKlG0RKBMcHMVJRKnTJuaiG/UVKkOcZjzdYmFTt6tsYVMoTHoMFt2x\nqdODKmuxzFtw6xQzedewWW9s443vL630TwVVehPVpENrDXTcblEqedd3JUont413DeMqaG5QPuy4\nOQAAIABJREFUSiX+C4jIJqoxuNoolIWRRrAoV5/CbcusMQVAL4LHJkSi86cJDdOV8vjnH4zFXBkL\nj4twyipiNwyiph89taK6NYkO366CpMJ+Tn/HCzieaAJ7c3Un7j/9nGsKzzP3nXYcQTeIasjqph2h\nL1ppQbd6ppN+nPVczVYiqBncYoIWNAtMlyAyMnTmAcIg1EjguNIPMVPF9bMCoCeW2rGLLytd2xVd\nLE9pmj+ylpiWGxeB02HtCB1XQfMGuicU27aqyM3CcOWNPRYtq2JRdMlwisgSMxE2AaFbpiIxsoFA\n1NzAeaXdVi4F40zXSzQ+q6iBaoWNRRSYn1ma8Kqy7YBN2HhK0c+5ptQ8szo9Hi9YC3QW9zJWmrSF\n3hQjvWXZdO/XajbYNzMbj31CEzNzwOR+uwspco+dHW66RQIzRHTvUTDrorRjARQLnJAygsWpmJ+D\nO6bXnklZ50oE5053MtLy0rWLMsFhYHNsScrr5m2z4MTjHJImllWR1237YqABPFxE3qO9dAqBObgW\n0rIw8o6rKjukqIixWTHSApSLuIFchIz+mi5qoD63mEsxO13UAH2z1qSJGXNcw0pn/MeW22kQc02Z\neeYQ64F7O6KmRyaW2kT/ndIWpUS6bkbXrKL1WzY3d7NvZhZIWlx2QaeUv06TrrE2VEfUQLfLKnKn\nFRUxZYSHzVLRT2w1aSK/U9W9lLayD5aEqen8ys2215MVj+vtYOnQCPe3gNcR9II4DLxKKXW99vrx\nwLXAd5VSv551rTwLzffCEz2foMR5xGHgD5zupiS2LIwssjtKV1pWvjRlxYxtEY4Egpm2bNZoeWPs\nqupQV7XcrGaRqdYaqEzcZIkZfZv+fg6iOaSNfrjHhpyBzDVF5xkIu1DXsHBGXZfLFJYre8zm5u6g\nzYFjYLCZUdWiwaIKatjohQEjllNcIlVaSopQRWdryC90F93zvBovHHTdD5LvQX2CxrE0ws3ALyml\n7g4fbN5PshXLawie8E7Ju15et+3rgetF5MNKqdosMr3yQ061bu/laTyr47Jreni0X9EKwK5ixoa+\nf5rIq8uCkyZsMq01UJvFJq0dxbBZagYhZvr51OTCqMw1h1gfLwZRkO7c/vR+PXmsbIKt+ztpvrlP\n2xnCxdbnKK9VQOvoFPNj0FBGrZqltt3J14QJAlETFfSbZyeoMNddBQt51IIhKgw3KAFTB1mfQfRa\nY2xZu+d6rDS27LghxKUR7tXa/l8kiGUj3P9nCPrIvQl4bd7F8lxOTaXUDHBQRLo+CaXUk/IuMGjO\n+8KNcGGJYyxk9WtyKeqWh4uLyRxblvAyjzfHVlfbBDNIussqYlproJ74mj5jy5DLPcby2QWfU30u\np34/NTmOaaTmGl3U9BJUuYuorktwjn37Z9k8vbuypoSmyEkTOHpwaRRAqmdDzREW1gN27e+Md45m\nXNd561IgzObZmXBN6fdStI9T0arBsSusJteT67jj/cZgQx0WkLB2kXPclEPrhJLkddt2aoSr8TLg\nU9rf7wL+D8nsx1TyXE6vCf/t+Qms37yIy3tOfzXFQlrJfX2xNoVNlmgoUgsl617013SxYC6WNutN\nv1oHOIkaKGat6SFLKBpPMAaBC3t7chqxYOC+PjU5MhJzzaEH1seVeKPF5JYeF6yEdWdTOVGT1YXa\n3C9N1OgiQA96NvsA6ewCFkJrTJCyHfQfitOKjXsoKtR0i4d5fGZH7xRh029LUevoFA01k4gz6pUd\nIswBVyr3e0kTkkU/jx88cFp8zD33fAfyu207IyLPIRA0zwz//nXgDqXUdSLybJdzZBbWU0rdHv66\nTSn1X/oPNbY+qIKyYsaMp8g7T5rrIq2IW1T4L0/M6OctVJhNG79e2M0s8KaPzza2XgrPuR6b6mrZ\nQr5YMV4flhiUouOoqebMRFQILvx5hfG6U0E5jbSnpsryl0dtrokm9V4WyDT3w45NBHEtgyZMu04T\nM/H4m8RhGLuA51W4eEe0jk6VWox1d5cucvSfIpSxaLRoxBlf/WoKmUZjbLnfvbKcSiOIyJOAvwZe\noJS6K9x8IfB8EbkF+BjwyyKS+bjrmrb9XAJ/us7Flm1DQbxQ7wXeW53PshcqaS9gfpRpi/6Fytol\nHLqtTlkxN1XW1onO35U6nRZXA4UsMGlWKRuRS+xyXgQXhsd+oXcrTUSq9SkF22dStPXBN1vnwljs\n/VlVSlXSmqCKp6aCjMxcY/ZMKltYz0bR2Ijx2VVnK00VRPc6R3dNml10gm97cTE5CRbjM4COVSav\nJo+5f+K8WfsXvKfW0Sn2NGdAqo+lSRNkeULb1uyyJlwa4Z4FfBx4iVLqm9F2pdQbCEo2EM41f6iU\nylwV8mJoXkXwdPRIEfma9tI64AuON9RXdDHzo91wynsLHJwiAvRzm5VybQtXVuZR3vlN4vPbdGlU\naTZtIY6OeW9S4JguKtv1bWNOE2VFhE5aQG6msCmATRzY7mXYRE2fKfrUdLHlqakBPBg4RUT25k00\neYziXNMLkZiZI90CUnTRyRM1eQHCNszxmQvxpFKgCbOoDs3WmSY03Tpo92lhLYRNFHXtU0DYbD+6\nyE7VXbenLlxcbf0QNY6NcN8InA5cJsF36b6yD2R5vZw+AlwCXBH+G/08pdcJrA6ihVJfSH50UrUm\nPtfmjKZ7x3Ye8yeVLCNb2gJ8oQqETI6FKisA2hXTpZZHlsCrw3WU5lqLnGuX86LgunuBV1X3fSly\nL30WP/FTU9gx+8UE/8djsp6alFI/o5Q6OzzunyqaC0ZqrtFZMAvQ5bCiCYBJ7dhogdPPV9Q9oIuW\n8dnVxE+vpC3AUUE9/fUd+/tXtK3LFZWTVVXETejimnL9jFo0UttcuBIJ4cn9OTsapN1HP9xPSqmW\nUuoxSqlHKaXeFG7bFTXDVUq9XCl1mlLqyeFPl5hRSn3OJZsyL4bmkFLqFqXUb4S+7B8TxLSPhxPe\nUHAqP+zeGE6B7/yJo9/StM7sJRYSpeNxNGHjLF6KUMSdlmF5yIqzcaVo3E30Ttjoh6gB+5h/tLv4\nubMCgocltkcnTIuOnppuBJrRU1P05ETyqemrInJtzWMaibnGxp5m0MzHNT5iF4Hlw2wUGR8/k/5E\n7RIDUVTAmIJgYqnNyqZkMHARa0K078RM2ykWyHQz2eJl8ih6TJnYpzxR4yIOqkizngP2TM+Ujt8q\nEzs0Krh2275ERP6DIJXzn4FbSAYJDpyuBauX9F/92AvVwBYlU0j9aDcdq0uZ2KAS7pRKBZiFQYsa\nnarFTLyP473000rTz6emIgz7XLP+uEOZrxcRNavT49bGhraeO11xIpq4qSrQM1ogVzYF45u0WF5c\n0UWNy+JZlaWgiLApK2qyBEFdQbc7ROLv1vEqP117Kky/N38SY12DosZJ0AB/Djwd+KZS6hzgVwhS\nOdc+FbogeqHMYmulZIxIncImzVpz8MLHlxI2roUHe+EgGwqlag+jpWZIGfq5xiYi5jL2txHtbxMv\nZZ/iyy6k5nFliwR20QzmmqjInss4BiFsytIvYWOKZP37UVSUmcKmqKh5YOVE2rsnaO+e4J5rTi50\nbD9wFTT3hkGBx4nIcUqpzwKVZFEMBVnBwIb3vqp+UQNFc6cVZRSsNXn7m+KsaM2isjVnvKhxYqjn\nmvV0W2haNBKxMFnsEGFhU/CUvUwjqLQbxlXMAWer8q4EsFtu0n7i8de48O8Im1wWWTjrEDZp91hF\nFeM6hY0uZubVuPX7YRt/nihea00pI1zTtn8oIuPA54EPi8gdwJH6hlWeSsz2QxKCWGumTA/3WCRz\nqyh5WVARVbwviSDyvXDKLIVceWldvLOo7TPdA5xY/WkHwNDPNbZ+QKvNced+SMwkF5xlGpyt4Eos\nXbZrFBu2LJcWDd5Tc7PCQWEW6Uu8pr3vZV0xWZlRvQq0eRVWYMaeqm4rJKh/x2wCJmpdsZZwFTQv\nAH5C0CTut4D1wJ/WNaheSFswiqTJZVXbrYvMxXFLUF++MrdTBYJtEKJGZ5isHbb3IWv8wzT2IWQk\n5hpT1CzTYGFTM7BIZMw1UUBnhIuASVsM9cWrVytD1mK/1sirNly0Zo1tnyorEi/sh83sTIw77fNK\nq5A8SOHi0DfuccDfAOcDf6SUenu4/RHAB4FJggSB9yul3p11LSdBo5TSn5D2ON7HwDllFv5gd5Dp\nlMZBNoBlgXFxW7iQ1johjbwu1T86STjlx/XULyi60NYtasxr1UFknXHmC8J55L9X/Rr/WmNU5xog\nCPDNSKddEckM6CwiZrr2q3ARndzU6ds0SPpRJ8W1lUK8v836YhE5aZ/HFC0WXQc3E7gn90zPpLvM\nUsZvXntQAcCOfeN+APwe8ELj8PuA/6WU+oqIrAOuE5FPG8cmyCusd5i49VjyJUAppSppTFc1By98\nPOfRWaj+4MEp+xWMhUgLXHWliAAwu1Sfx432LoBpQcuOrpNerQV1ihqdOgRCYTETUuY9MwsyepKM\n2lxjW6yCAN8Ut9NM8P90PnzSrjIrqWoi99mKSFA4L2JJm2uanQrBLpbvUei2bYqGIvVlipL3nq0U\nyGiKx5ExflchVgMufePuAO4QkcQNhO1Qbg9/PywiNxK0aCknaJRSTh0uh5FYBITNDn+0m2DxjxZ6\nxyftiEEEA+fGaGRlYOn3moLLvWd18x4keZ9HEdEQufJOmc3ZUQse7+X7kNexvR/00wzswqjMNYdY\nb10com2R20knCgS+UnWetE1Rk/UEXocbKOt682Owb9NsZ6MmZFang2DmqekWc9LOzIiKAlrn1Tjn\nlOj2POgKwlV3qG7QCpp+OjC5CZjpdjUVIW/8ne9xTzFTVXfbtiIiZwPnAV/K2s81hmao+SGn2heJ\nC4N/EtaNVwls6XYF2RaTOhaeMhaNuDw/oespWljTYmEiq4NN1IT3X4R+xxNVQdG4FsholfGqnPe7\nIpJjrq8mTb/NwMcKLRrQhIWZZiwCVqfHWdgf/Lv9aNLZUOeC3cu5495DAGxmKkwtj+IwooVwSrVY\nWAosObtItzpEwaxrhV5cgys52XArIkzuD6xkyzRq+Y70dM476Kwv3wMq7LadRpgk8LfA7yulfpS1\n75oQNGlEQuXghY/v9OuhszDrC8iwp2NHoiaLONYmtEoB3aImXJTL1ndxOa5fbqii2ISrHmjtEpt0\n8MLHD+W9FaSvZuC1xKEH1mcuCC2maDWNQN/GFBwtf81BWCryGjtG6c6N6RZb9zdhkz0Yeo5AFade\nR7NANcIaKcsVpFL3G5dYnCLxM6P4HqTg1DcuDRE5kUDMfFgp9fG8/de0oAFDtFw48gtRJokGj7qo\n+YJ0fq/ZyhAJhlERNafMVpg9NjxMGK0KBmoGPtYoK0DMxpJV9F8qi6u7pUUDpgNXGwAzErvc5ki6\n2bKu1RhbZipyybj3uh1K0sRNYOHKdu9Mbgqsea0S1pmqG5NWRG637TQk6FT5/4AblVLvcDlmzQua\ntYRr3ZM4MypFvFTpOhrFYFfdcpcabA3WGKVhFGkAXJmwxK6W7VbrShEz8FrigbtOoL17omuByFpM\nymK7TtHj8yh6fttiHbnati4FHbajKshFAlp1pipOex4kTjFQYcD4anOcCUsrjDSKfOfq+H664NJt\nW0QeBlxLUJ3kgTAmZwPwJOAlwL+JyFfDU/5/SqnUaOZjQtAsaZK/SANF/fgyx6XRi/WiSLp4oh6P\nJm5cXUdZBeBcso6ibcMkApxFl0XMHLzw8ZnHD9N95tBXM/BapF8LRN3X6cUqlLAgjEFrOrDWJOrj\nOFgZGmPLcRNLWyuIfpL3fpcRmJ2aMcH7Yo01CoWMXqMoi0EJlDKEAqRlbNul/f59gjnI5F8Jshyd\nOSYETa9iRD/etkgPq2Ui4YKqkKIL9zC5n6xWLtOSlZI9licEh+k+c+irGdgzOugLZRlxo7uOXBbn\nSMzcH1ZZjlpC9JOylo5K3DhLAjOd7DHAyd00Prs6UqKmXxwTgqYqzP4/I7J4dVlrBpGtlJVFVocl\np7Qlq6LrDvN3o99mYE85Br1o9bJ4u1haGmPL7FuaZWVTp5hfFAxrLujRWKqKBanifa1K3OhixtMb\nXtA4kuVSKXu+qha9qlsFpC3wVQshc8zm32Xen56tZX3qrj5o4dNPM7CnOIMMCLZhEwC2MbaOTnEO\nbsGsDVrsCNOY5/bDZEp14qrERzTeuuKdoODntiRdYsYldmiglpnD98KBleiPwY0jhYEIGhG5FPhd\n4M5wU/yEJyJvAF4G3A/8nlLqHwYxxqLYRIBr+4Q895WZXh6JIX1fF1FTZD+TPFdLHRQROL2MIQgO\ndss6TnM7DbNF5lhmrc01tgVzUAtcloXCJaC5RYNFFbinVqfH2WOxzrjemykobMf1433Kuu8pWkyU\naCnhXUvuDNJC886o+miEiGwg8Oc/Afhp4CoReYxS6v5BDNBcJG0F5rLcFa4F6fIWY1O4RItn2UXU\nJoYSZNS7qSp2qMx50sY8rDFMWaSN3wujWhj6uUYnbUHM2z4s7qmsbTqt2dCa08C5Tk9ehtkwLP6u\n2Wm6W66r95LWx6ofGXVrhWFzOb0A+JhS6r+Bm0XkWwRFwA64HFx0YXN94jdFS9GYiyrjV9KsFq7W\nGdvvVYyjCOa1e+nTVFlPJ8S5r1OVlZPT3G5/VvhMi0BK0zKPjZ7mmrowF68iLoxhEDZFyFv40+5j\nFO7Pdm/LNNjabMI0rLIu3lY0RX3UPud+MkhB82oR+W2CwMP/pZS6m6Dg1xe1fW4Lt3UhIq8AXgHw\nsLNOKLWwpcWxROeqOjso65xlF8eyi3+WoMlyS/UjGNrF+tJXq0xKPR9d2JSxsoyiZWlEqWyu4SFn\nVT64quNl0ha8vIVwEAtlkcDatNeHdWE3761Fg62EbTHC2Jle6u2UsdwM4jN26Bsn4euRrW5WKfUV\nl2NNahM0InIV8DDLS38EvJfgIVSF//5f4HeKnD+sfPp+gA0XnOTWWtrCsbqoFOld5XJs0WsPk2ul\nrkwnnSKWNE8x+jnXHH/+k9WYowAZ9EJb1HWlv15k7FXF9WQFHhepgjvo991Ge/cErdmpuN5Mnpgp\nU206673Ksvy15wtfyhnHvnEXA48Of55G8H/2aY7HJqhN0CilLnLZT0T+Cvhk+GdPBb+qpA7rjOv1\nqswmMu/DtY5KXvBwkffHds1hW9ijqsEJt1OFbSKG7X7XEsM615RZjIeFPFHjak2pSlyUcb1FDErg\nmO9Be/cEm7ftpkGLrUtNpqa7+1ZV0berxu+Ugnu0P+8Jt2WS2zcu/PuDSikFfFFEThWRM4GzHY5N\nMKgspzPDJncALwK+Hv5+BfAREXkHQaDeo4Ev551vrP2TwgLEtsj2W8SkkeWW6nWMLsXhiozJ9Zom\nQ9u5W++B5UjZ96Zqy1AR+mkGHiRVzzVlGFbRkkcV4y7qJqrjvSoyhoTlIiVzqugY9f2DVgj5JZui\noOCu3lnGtrR96kAp9WORRwF/CkwAHwA4s4K+cbZ9Hu54bIJBxdD8pYg8mUDd3QK8EiAs7tUkUGD3\nAfN1ZR0Mi3gpQlVj7ledGdcx1HXdollpMTU38IwYoJjpqxl4wFQ616w/7hDPzevNQ/86ZNv6BOVd\nu18LYBbDENhaNGanKrG1/ehi0CbCgv7Z2D5bc5vLPnn9pLLbZZo8C/hn4PHAz6LU994OvD37mP4x\nEEGjlHpJxmtvAt7Ux+EMFrML9l76tqCamAKgXwtuEbdY1nG97tdFQUtNKrrVx/bZ9qkTukZfzcCD\nZFBzTW5Dwhopeu2s/esWO1VbZIZBrLng1LQyg6hflktQcZXfRaX+RkQepeCzwF1nOhzi4tpN2+dE\nh2MTiLI1yhoxROQwcNOgx1ERE8Bo2qiTrJX7gMHcy88qpc5w2VFE/pzAohLxIwjzQgMSZmARmQae\np5R6efj3S4CnKaW2a/t8EniLUupfw78/A7yOQNBkHruWWUNzjf//OXwM6j6c5xoAEXkKcK5S6gMO\n+54AfBP4FQIxcg3wm0qpb2j7TAHbCdzbTwPeo5R6qsuxJsNWh6YsNymlLhj0IKpARK5dC/eyVu4D\nhv9elFJ/DPzxoMdxjLAm5pph/04XYa3cy6jch1LqOuA6x31z+8YRtGBpAN8iiNd7adaxWddbK4LG\n4/G401czsMfjOXZx6BunAGvyuO3YLI4rOUaPxzO6XAM8WkTOEZEHEbQAuMLY5wrgtyXg6cChMFvI\n5ViPx+PpO2vFQvP+/F1GhrVyL2vlPmBt3UvfzcBrjLXyXVgr9wFr517Wyn0MjDURFOzxeDwej+fY\nxrucPB6Px+PxjDxe0Hg8Ho/H4xl5Rl7QiMjzROQmEfmWiLx+0OPJQkQ+ICJ3iMjXtW0PEZFPi8h/\nhP+epr32hvC+bhKRXxvMqLsRkUeIyGdF5AYR+YaIvCbcPor38mAR+bKIXB/ey45w+8jdi6c+Rmme\nAT/XDNu9+HmmTyilRvaHICjxP4FHAg8Crgc2DHpcGeN9FnA+8HVt218Crw9/fz3w1vD3DeH9/BRw\nTnifxw/6HsKxnQmcH/6+jqD40YYRvRcBxsPfTwS+BDx9FO/F/9T2HRmpeSYcs59rhuhe/DzTn59R\nt9DEJdyVUvcAURn2oUQp9XngB8bmFwB7wt/3AC/Utn9MKfXfSqmbCbJNntqXgeaglLpdhY0KlVKH\ngRsJGomN4r0opVQ7/PPE8EcxgvfiqY2RmmfAzzXDdi9+nukPoy5o0rp0jhKTqtMN+PvAZPj7SNyb\niJwNnEfwxDGS9yIix4vIV4E7gE8rpUb2Xjy1sFY+85H+To/6XOPnmfoZdUGzplCBrXFk8uhFZBz4\nW+D3lVI/0l8bpXtRSt2vlHoyQdXbp4rIE43XR+ZePB4XRu07vRbmGj/P1M+oCxqXEu7DzooEXYwJ\n/70j3D7U9yYiJxJMMB9WSn083DyS9xKhlPohQRvZ5zHi9+KplLXymY/kd3qtzTV+nqmPURc0a6EM\n+xXA1vD3rcDfa9tfLCI/JSLnAI8GvjyA8XUhIgL8P+BGpdQ7tJdG8V7OEJFTw99PAp4L/DsjeC+e\n2lgL8wyM4Hd6rcw1fp7pE4OOSu71h6A8+zcJosD/aNDjyRnrR4HbgXsJfKIvA04HPgP8B3AV8BBt\n/z8K7+sm4OJBj18b1zMJTKNfA74a/jRG9F6eBBwM7+XrwBvD7SN3L/6n1u/JyMwz4Xj9XDNE9+Ln\nmf78+NYHHo/H4/F4Rp5Rdzl5PB6Px+PxeEHj8Xg8Ho9n9PGCxuPxeDwez8jjBY3H4/F4PJ6Rxwsa\nj8fj8Xg8I48XNGsQEWnn71X4nM+PugyLyAtFZEOJc3xORC6oemwej2cw+LnGM0x4QeNxQil1hVLq\nLeGfLyToBuvxeDyV4ucaT1m8oFnDSMDbROTrIvJvIrI53P7s8AlmSUT+XUQ+HFbkREQa4bbrROQ9\nIvLJcPusiOwUkWcAzwfeJiJfFZFH6U9DIjIhIreEv58kIh8TkRtF5HLgJG1svyoiB0TkKyKyP+zV\n4vF4RhA/13iGgRMGPQBPrfwP4MnAucAEcI2IfD587TzgCcD3gC8AF4rItcD7gGcppW4WkY+aJ1RK\nXS0iVwCfVEotAYTzk41XAUeVUo8XkScBXwn3nwD+GLhIKXVERF4HvBb40ypu2uPx9B0/13gGjhc0\na5tnAh9VSt1P0ATtn4FfAH4EfFkpdRuABC3tzwbawLeVUjeHx38UeEUP138W8B4ApdTXRORr4fan\nE5iRvxBOUA8CDvRwHY/HM1j8XOMZOF7QHLv8t/b7/fT2XbiPjvvywQ77C/BppdRv9HBNj8czGvi5\nxtMXfAzN2uZfgM0icryInEHwFJPVsfUm4JEicnb49+aU/Q4D67S/bwGeEv4+rW3/PPCbACLyRIIG\nbQBfJDA7/1z42ski8hiH+/F4PMOJn2s8A8cLmrXN5QTdXa8H/gn4P0qp76ftrJT6MbANuFJEriOY\nTA5Zdv0Y8L9F5KCIPAp4O/AqETlI4D+PeC8wLiI3EvisrwuvcycwC3w0NA0fAB7Xy416PJ6B4uca\nz8Dx3bY9CURkXCnVDjMRFoH/UEq9c9Dj8ng8aws/13iqxltoPCa/GwbufQNYT5CJ4PF4PFXj5xpP\npXgLjcfj8Xg8npHHW2g8Ho/H4/GMPF7QeDwej8fjGXm8oPF4PB6PxzPyeEHj8Xg8Ho9n5PGCxuPx\neDwez8jjBY3H4/F4PJ6Rxwsaj8fj8Xg8I48XNB6Px+PxeEYeL2g8Ho/H4/GMPF7QeDwej8fjGXm8\noPF4PB6PxzPyeEFzDCEinxCRdvjzT4Mej8fjWZv4ucYzCHxzSk+MiPwFsKKUeldN5/8y8FKl1Dfq\nOL/H4xkN/FzjqQNvoTmGEJF1IvJmEfmWiBwWkZtFZKeInCEiZwC/DbyvxiG8HfjTGs/v8XgGjIg8\nU0QOiYgY2z8jIq/zc42nLrygOUYQkVOBfwEeB1yslFoH/CJwIvCzwCzQUkr9uMZhXAE8R0QeVuM1\nPB7PYDkP+KrqNv8/GfgKfq7x1IQXNMcO7wR+AEwrpf4DQCl1m1LqlUqpa4GLgX/WDxCRvxSRv9P+\nflv4lPUg2wVEZFxE7heRM7VtTxSR20VknVLqJ8B1wK9Vf3sej2dIOA84qG8QkbOBh4Tb/VzjqYUT\nBj0AT/2IyCOAlwC/qJR6IGW3nwduMra9Ffi2iJwHPA14HvBMpdQ9thMopdoi8u/A+cByuPktwJuV\nUofDv28Ezi19Mx6PZ9h5MmDGxpwH3KaUWhURP9d4asELmmODi4A7lVIHMvY5FTisb1BK3SUi7wT2\nAOsJJphDOde6hnCSEZFnARuA/6G9fhg403agx+MZbUTkROAJGBYagjnhK+Hvfq7x1IJ3OR0bTAK3\n5uxzN7DOsv0ggfXmDUqp7zhcK5pkAP4S+BPjKWsd8EOH83g8ntHj8eG/Nxrbn0owN4Cc20gXAAAg\nAElEQVSfazw14QXNscGtwMNFJOvz/hrwGH1DaBp+L8FT0+84Xusa4HwR+Z/Ag4GPGK8/Hrje8Vwe\nj2e0mAR+rJS6L9ogIhPALxEE6oKfazw14QXNscEnw3/fIiKnAIjIY0Tk3eFkA9AimHQIX3848Alg\nDtgG/LyIPFt7fbeI7LZc63rgYcD/JXjSUtoxDwaeAny6ovvyeDzDxb8BPyUivysiJ4nIY4CPAp9U\nSn0t3MfPNZ5a8ILmGEAp9SPglwmeiv5DRH4IXA4cUUqthrt9EGiEk9ApBJPOO5RSVyiljgJvA96k\nnfYRwBcs1/pvgkntFqXUp4yXLwE+p5T6XoW35/F4hgSl1PeBGeD3CLIqr6STqh3h5xpPLfhKwZ4Y\nEXkzcEde9c4wlfJ64ElKqXstr30LmFFKfdF47UvAy5RSX6925B6PZ5Twc42nDryg8VSKiLwJeKRS\n6jcGPRaPx7N28XONx8S7nDyVICLni8gh4FnAqwc9Ho/Hszbxc40nDW+h8XiOQUTkA8CvE5j9n2h5\nXYB3Aw3gKDCrlPqKuZ/H4/EMC95C4/Ecm+wmqMaaxsXAo8OfVxCk1Ho8Hs/Q4gWNx3MMopT6PEEW\nShovAD6oAr4InKr3zfF4PJ5hwwsaj8dj4+GAXq31tnCbx+M5hhGRnxeRFw96HDbWRC+niTFRZ6+3\nvPDQvg9leLmj4P53lrjGGSWOSUP/7O5I2W6+Vhdl3oteOAOu+z6rSinnd/TnRNRR7e/b4RvAT7RN\n71dKvb+yMR6jjImoUwm+6ic+Eu467TRufeAsHvjmCXDkXs5+ync4/e674W6475HHcYj1wc8D63ng\noKL/X6Z+cAbHnSeJLeuPC9owhXcPEL8Xd995OueecZATvh30yb33bvd3Rf8PceJpwGnBZ3DLnY/k\nuNPv46zjbmU9hzjh7mQP3vtOOy5xfW7ttTzNGXDyiXA6nHbGXfF9mtfVr38rZ3HogWChWn/cocQx\n9357cN+M2yk418DXfgCIyFVaHbNUROQPgJcDiqBu0EvDbujR65XF660JQXP2erjWVix7e9+HMpzs\nLLj/+wjaxxUlqhLxyozzZmEeF31+5vj1z7XovZUhb9wGO+4K/l04veB1tPuXN/NfRQ49mjycS+En\nSqkLCo5A57sEBc0ifibcdkxzKp33ee7bMLn/bvZMP5ftRxcBuGbsHCaWgCbsaU6z/egi7Ysm4MAK\nsMvxKnMF9h0CFhf4wLbNAEzRsuwwBsA5R79J++T7GT9yPF9tncHKpuDVMnc6B0FHqLthctPdrF57\nExNLbWjCyv7kvpObAB6ID1i9dpUzLjsM8ztKXFnjSPBz98MXuPeqVS4dm2eKVjAOg9XpMea5lNbR\nKdq7J6KhAzCzbQ/7ZmbZsb/rsNqYA45X49wvbR6G+1zzchF1AkGZ5zsDDSZZ+4dVoH8P2KCU+rGI\nNIEXE8TwRejxek8jiNd7WoHbiVkTgiYV22JnWyS3W/YdlBgyx1U15uL8yvTXooVZx2mRLigAuo6L\nxpQmVvohYnReSeKebO+LDXO/hdOT2woLnv5yBbBdRD5GMLkcUkrdPuAxDZy0x9jDLeOVGQwxA9lC\nZc7y966M1/PokyBaXGC3Jmbsi/k4AO2T7wd2sXPsBmj2dtno7uYgEALSfd2Y/ca7J23UfkEWVe+i\nBuDADtonzzG7uI/x2VUa08s0aCXej2UatI5OATA+Gxg12ruDrjPNy7bSaLZYmGn2TdRM7geW2rAf\n2OR+3D8Dv0JgKP8rgj5dDlaaE4CTROReAnVrmsfieD3giyJyqoicWWa+GaigEZFbCFq83w/cp5S6\nQEQeAuwDzgZuIagCeXfaOQrjuhjuZG1YePT7zRIaZUVIHbyPdCvPCOMqhPqBiHwUeDYwISK3AQvA\niQBKqV0E5egbBJVYjwIvHcxIq6HKucYmLVanx+PFK17An6GLmayj866UIU42Tgb/dl2nT+yF1rYG\njcgyYxEqE7Tj92RGnc3WpUtZ2V+N5HI9R5c0bMLu5mZm2VeNqGEXzEN7HpobF2hdNUVrbJnF6e1M\nLLVp0UjsHYmZiNnL9kET5qRZuxRd2ET8Oa02x4EMMaghIif9NBBlBjwB+H7QH/DJ2m4J17ZS6rsi\n8naCBsk/Bv5RKfWPxqnT4vVGS9CEPMdQeK8HPqOUeouIvD78+3V9H9VaEDMuWITMUCy8a1TURAza\nOpNXXTV8Wprv03D6Rc9zzYmPhMmnBAvBPDtpXrYVngFs0Xa6LPy3jMiIBErEgRWsIsjcb+PkYETN\ngR00BZob9zF+1SqHOSN2+UxGT/7NQNTMqCvZNzNbmZgpyy5gbj9sXWqy9aFNhEtzjtDff3PkFqtb\naLFpbtxK4+oWU9OtwNV00QRdRN+bvTDLPu5ULeakXe/7M0MsaObZSbLNVibyCDrvhgCfCfplvTz1\nAJHTCCww5wA/BPaLyBal1N4yQ89jGASNyQsInhwhaCX/OeoUNNuNf/vE/jMvYdPtn+h+YcBCKk/M\n7LhrQIux6So036e63VCa28l0HXlGlsJzzV2nncae5nO7Y2MOpB3hYJExxckWYG/4L5Pd++fS63Lo\nYBkyObCD9u6FxKaEsGnCvmZ/40TyWNkUjm3jAhxwtdK4WtiC78T2o4u0xpYN12PIxsmOu25bYOGa\nZyfvcRcYpVjRXEz7lmZ79f7lcRFws1LqTgAR+TjBI4AuaCqL1xu0oFHAVSJyP/C+0FQ1qfnOvk/K\n/2gReQVBwS/OOqUfQ7Wz/8xL4t+tAiXnmGGiaGzIQIVN2vY+xNf0KmYGbZ05RqlkrvmZsyRwD8yv\nAC6LoE0UpCyKWyy/b7HtSHI5yL1eWQoEJ29cCGJDHhq8iXpgbvT7MIU57yK4u5X9sFttZlbeQ/oI\no71Jik+bOzHx+g7aJy+GgsH2XVnoxNmEqmKxuZ37i95MQbS76TmWyYFbgaeLyBiBy+lXgGuNfSqL\n1xu0oHlm6GN7KPBpEfl3/UWllBIRa2+GcEJ6P8AF54pytmz0uOBliZFUq4vBpts/wf4zL0k9l6sw\nKkRKwG8v1oZarTVl3U22AO8hwIuYgVPJXHP8+U9WgZgpuzxnPOHHVpkSVOZuMhdlB1GzcYGZq/cA\nYSxRs83kpuEUMja2LjVBkSNqLMRuvvA92jjJ+FWrNMaWQxdTvuVnmQZbm82ONSsM1J3bVN/7Fn0D\nJzfBnuYMSH2qRin1JRFZAr4C3AccBN4vInPh65XG6w1U0Cilvhv+e4eIXA48FViJIpzDyqTVVhrp\n0aVjig1TlER/1yJKipKysOsCpldrQ22iRo+hKfqZVWmtccj8qpztwJuLHfLTx8OClmp/qXeJJahq\nrnng4B2UWmo2LiTiJYCOCDmwEiyOW4IUXggyX3TGZ1e7AkljKhEzRbOoOsfNXL2HRbazbAS+jgxN\nmJpuwcZ9Ga7DhF2jQ0LUJBm/ajVww2UI4BYNtmpmktgNVjOTm4AZwmDles00SqkFgqQDnV3a65XF\n6w2sUrCInCwi66LfgV8Fvk5gfor+N28F/n4wI6yXLMGjW29slpws605MH60UQxtP0ms8UsnMr4XT\nu39yWcMB0INmGOaa8atWmdm2JxAskbDZOEkwr++CAytByi8tGrTi1F4g8TvQseQkrDnheUqhLcZm\nPI/5usH4keNpmOnaMyWH0Ud0KwUQjD/XOqa9v/r7r2WatXdPxLVmdo7NBzEytvc0FLito1Psac4w\nuanzCe6oMWja/CQb1rpBo8sgLTSTwOVBkUBOAD6ilLpSRK4BmiLyMoKCP87/PWyLfF2WkixBYbum\nablxiaMpFWszIHdL3cHC5nvh/LkWdUH1kL5e+v7LWqI8rlQ714QultbRqbi2iitRLRLAsKzson3R\nAq2rG8lFZi+0meiImr2d7cWK9aVhiQ3pwnKNjQuMX7XKzrCY3J7pmTg1Wbc4DLu7qTCRNc1EyzKL\nLGlRsUWT8SPz7BzbHLi6WliLAVaFTYpO7if+vILvY/2BNP1iYIJGKfVt4FzL9rsIAoecufvE9ew/\n81nW1+oQOUXFjOuxeccNhRurX2iLvO09KxSMnZcRFYkeBzFTqWjTrTJd93kMfdY1U+VcA8CW8Ml2\nDJqLW2FejzOxZwi1d0/Qmg2e3DtixODADlpH56PCuh32Qnuv5m6qRMiY17Ys1GkFAUNBFxWPi2gd\nnWJxbHutC3SvpNqamjC+e5V2WceH+Xnqn5n+2uICN4+tS1Q1rlr0pd1jZI3aPL2bptxCUHrJLrrS\nOHECJp8T/L7um8D1pYZYG4MOCh4IZTKTXEg7V5qoKipubK6nxDXTLBF9KJo3yKDXnsReXVaRIi6k\nFNFWJyLyPIL+KccDf62Ueovx+nqCpfcsgnni7Uqpv+nrIIeRk3+6yx3U3jsJW8JYiUQ9mDlm1JVB\nPMz8Dtrzc+RlRrUvmqB11VTmPsVaKOQdowmxeNxJ8ZJcrOfsYoYGO8fmmZhpD2UwsLnIR4t7V82c\nXHYBC+lZZonvgPEObFzgzm3rmJhpx9eu6j0yXWipzERWwmJCZlRYE4LmLk7nQ7yEl/Chwse6LCRp\ni+WHeAlA5nWzLERlRI3LuPKwiQ+zJH9WXEzt4iXFOpP1fhcSNVXWsKkw9sXl+1QFInI8wYz2XIKq\nnNeIyBVKqRu03eaBG5RSl4jIGcBNIvJhpdQ9tQ5u2Am/+5G5vr17ojvdessk7A1iIBu0AsvMPDgt\nXwd2BNkxWwiDSclJFc4iJZA1dd8UFieJF/Hw+mbPptg6k3+2vmF1t6Qs+Cv7YWdznlk2uJ1cD/A2\nhaA1sHiO3VdvZmKpOsGn3198XyMQv1Qna0LQlCFaPKCzgOjb9O2FLQBa8be84N+IIrE1dbieTJEy\nFGnGO2HTdnfRV8pSU1bMVCVkwjVA/+4Fv9fqcnoq8K3QFUNY/+EFgC5oFLAu7IQ7DvyAIO3SAx0x\ng9abZ+9kx3ozG7iZpmjRGlumWaSmy4EdyXTpeIEs07Syx2VTy7xiGzSfsbVrl3l20r5ogonXtgfq\nanISMPqC3ySRXj5FlOmUU1/IGjycl9o+GcTM9BCuYhUwEY5CZnV6PEiM7hMi8liC9iIRjwTeqJR6\nl7bPbxEUtBSC9iSvUkqVcmataUFje+I1RUvaNvN426L6Ej6UeiyQXtU2A9Nqoy/Qw1qQr1IsQiHN\nkmWzyvXkfjKaUPYr86imz3VCRPQCVnqPFVvvFLO77U6CLKDvAeuAzUqpB+oY6EhxF1rV1zAIlMlg\n0b86EDqNseVg321BnZEGLVDQvGwh7BvkIkxsrxcRJ+by3puwiQKVm2ylQSssmU9c0n/31Zv7Zh3I\nsjlZLTBp4wpbAESiZmKpHaRan1z26hlsobCYSXOTAYl7ivpkudIYW6ZZqDpyeZRSNwFPhtgy/F3g\ncmO3m4FfUkrdLSIXE9R88t22s4RJrvjQ+MQNm7hkQ+dRw1w4zfNEr8eZTNvDBdV4+tcX26w4nrQF\nuWt7lmAadIE5Uxy4HuOAac2oXNSUJUu4ZnV+r4IzgN/R/n4zq0qpC3o4468BXwV+GXgUQTG6f1FK\n/aiHc44+R75nWQg6JTZiMWMQu572LljqltQUFlpR40o962r8qsACpVuoZq7eE1gfKqawcImwCBjb\noj+x1E6IGoCdY/PMUrDAnkPhwd3bNgeNNTTygndjjPux3YteA8h0B6btNwB+BfhPpdR/6RuVUldr\nf36RoPVBKdaUoMmiiJiJ/jVFjfP5w3aksbChI2BKpx+b7LT8nrZAlhEXZXllyu8VUESU9oUigmS4\nUrJdeqe8FHhLWPTqWyJyM/A44Mv9GeKoMAd7A8vF+OxqLGji9GwzY2kLcMDMiKqayFVVwfn3Ats6\ni+TU2DnMs5PG2DKNba2u0v0RvV7ZOUg3xypkLv7Rgj5FK35tQus2PUWLGdWgKRlWtCIVncOMsCla\n8VgnteaQMSn3kSdebNhe10VO87KtfbHOWHgx8NGcfV4GfKrsBdaUoInEiI4uSvL2te2TdnwW8aJ7\nZgFri07Zp/ksi8yoFm7bSRyLtP/MS6yiJs1KA+mVnYexMSiU/84V5Brg0SJyDoGQeTHwm8Y+txI8\nUf2LiEwCjwW+XffAhp8z6FquD6wAk7QJUrOhU4sk+jubEs0gnajgfBbrzqL2H8UUM5VVue3BfeUi\nAmzCRmf8yPEpricjy8nWzynej06Kvzmu6cxbsI45omXZXqhAXiJDaw5yu41rTND5bD4JXM9zROT3\ntT1013aMiDwIeD7whrRTi8hzCATNM90HlGRNCJpDPzktVaC4CJcszAUmOp/LovMhXhJba1IxrStD\n2Ido4DiImjRy41MGLGL08enWwTpRSt0nItuBfyBI2/6AUuobRn+VPwN2i8i/EQTrvU4ptZp60mOF\nk0+EJ6UUoZtfCVK4NbprkehPxqalRuujNJgnaAu7OlYYjcTiPIDMmqJxI2mYrpro78bYMs3FKObJ\n4EAUBxWhfW5x9lMQoxIFiRd19dhES97+aaJG//xa042kCNs4mdHuwYnPKqVe5rDfxcBXlFJW/6eI\nPAn4a+DisD5UKdaEoHHm7Snb/zD7MFPUFHmCzk0nNxfUtNiXnRQXPUNgcbCyRkRbkWDerCy2QbjR\nlFItSM6AoZCJfv8eQYsAj87pdPdk0jEL1KV1ZM4qjrcFt6aQfWGOiaVLU1+tSlj0SrRgp41HFysR\nkWiIhID5+sy2PTT3pgXP7oqrJXe5GSHIcrtogfZuaG3riI2iQiUeqx5sbkEXM7p7KbqnxPuycbIa\nd2QxfoMUd5OInAV8HHiJUuqbvVxk7QqaNPGSt2+KuOnlqTktHTwVF1GTd/wwU8YaZVhpTIrWIep3\nxph5Pf07UbdFxtMnUivtljxX4sm52qwlV2bUlbCU3KYvjlUGmWYFs2ahW4+yhM2UIVwatMLOWfZ7\naNCiddVUTtaTJia0eKnW0amwOWXQ26krlsogIYZs18gINo+wvX+6dWbP9ExwnS0QdAPpD2H/tOei\nBT4Y1uA3EjwqXBa2J7mvbELD2hU0NdJLnEPp2jaQL1QchEwVC3ntLRiiAOaCcT9liysOEi9mRpBE\nH6WUXkimmLFabVKekmPLT2SlmYu7cXf2Wci28lTA+JF53iNngDbV1SVmbOcrK3CATBdZwrUUipoI\n3RISb091Pa0QeGzpOldjbDkQNWZj0QK4WmQg+V6Z9x5hFW5b6NXllItS6ghxScp4m24Nfjnw8iqu\n5SRoROQxwHuBSaXUE0N/1/OVUn9exSAqJ886oxdwf6n2e47rSScSNS4xNWmLrClqMkVOSTFThyWi\nzDlTg3CzrDTvIylqNCsNZw7GVVMlXsx0M/RzzV3ArbqQ0ETNFpJuqMi1pIuRvdpr0Hk9EjkHLNWB\nw3MnFsdZaF/Uc/xDCnNhdeBVJjcFYsBcJG1iJsud4hq02ouIyUIf/+r0eMJakyVqstlFUFC7g3ku\nfbsNXQBlYTvefK8Sn5GRRbXaNCxWWW7TEcbVQvNXwP8mfHZWSn1NRD4CDMckU5aX5u+Shilg0qw2\neRYD575SWQu/RcwMWxG+OpprVpnGbascXRe1CLGHkvwevLn6S/SJ4Z5rjtxL0ioyZ7fO6NvM1ggu\nIkS34mzsVCDWF77mlq3Vx9rEXbQ3B120mzPWGjO2mJRCmTbGuXohTr+2WSb0oYfBy7pbKrqPSIhE\nlpVuF5D9fW7vnoBt9nE1xpaDFGlI38cigLLeR9t71XXf2j2v7A8yz+bZmagdtFZxFTRjSqkvh/6t\nCF8GHRJWGlPUFFkYMy0zBdO4h03MRJQSNSlWmoi0thV5ZO1fhbhJSxc3SS0rUOqqa4IRm2t2aQGW\nFotMGKsQt0fYPWEIoF3ZAZqLk53WAxiL3bag/o1hJOgQihOIqhunZE4tLgRNNsOKv3rcxdaM8rY2\nUVOELCGT5jaJiERMIktpWsvGsg3bEDfRvlPTdlEDdLqjA1ZRM7+SCPqdvWwf47Or7BwLP5RQyGw/\nuhjU7rHcs4v1xST1/bHd9wyViZl7J45jdToICDrynXuA4Wrt5ipoVkXkUQT9XRCRaeD22kbVK39I\nttvJZpnJczdF53N0S2UthE7tDHQhY4qaGsRMPy0UuWTFzqTce9Ex11WkL6vuUPTZRNfuQ72ZUWS0\n5hog2XNpzq3oWiKOJqWJ5MZAzKQ9sTdoBYLGaj2Yi8UM4JQ5NX7VauZCmhZo6ypqehEwaftPTVti\nRwwLhU5cJyfaJxQ2pqhxZxeto/MwFgb2zq/Q2HZl/GqXOyslvTrtvSn6vphELsP2ZRNrzr1kw1XQ\nzBP0V3iciHyXoPeCa63E2ln/4Ls5ZG6MhIdLtpMmUmw1Z9L2Kbow2Ra70h2iLfRimdEtHaMQXNtr\np3LIFzVF34Pag6WPDYZ6rsnGcD/FwcNB0b14mxWLqAmLsunxHmamzvjsatjJOzmOGXVlYkvQdiEl\n5mYvMEvHqkDSAtKikSimZ2IbnwuFFusUQ9FEMzjHyv6os1Y2kcBJCBuLqIkFigPtkxfD4S3CxoX4\nM9MDjxtjra7A3dR2DK5ooiy6j9SaQHpXcD1ma43hJGjCrrwXhelXxymlDtc7rOLorp8EBQJ9ixAt\ndi5P+qUXuoLp11U0sqxLyBR2qTmes2h38rrccS6fsc1KMyhE5HnAuwnSNP5aKfUWyz7PBt4FnAis\nKqV+qe5xDf9cc2e5w3Qhk5qdpImaxUl2b9ucWARTxUIixXuO8SPH2/dLs9Ic2EF790JXnMcyDbYf\nXQx+NxZkGy5ixmnBdmgLZVpeijaTiGJLTKwWpyjNOU0ELC7AfPiZHlhh+9FFWqF7SW8NoY95BZjc\n1I4FSGY9n7z3IxIzFlanx4NGovrYE/cxDLWOqiNT0IjIa1O2A6CUekcNYypNnqUkLaukiOn/kg37\nuxZ9c3GKXh/kE7trDEdV5x8URS01ZoG7sjE4tnO67mteu9frFyXsertIUBviNuAaEblCKXWDts+p\nwGXA85RSt4rIQ2se00jNNZmYZfHNIOFUOingkZjJKxjXGFsO3U7ExzbG9qS7ThYn7TE38ysJQRNZ\nZtoXTTBz9Z54W5Wp1F30IGT0v6sqGdcYWw4sW4nYE0MQblwIC/BtDUXlLtonB7fSUCTETJdTcH8Y\nYWUEKydw6flpdN6OPqfoMwwCkwdXeTqcS/4aeCKBK/l3lFJdtkIR+QWCd/HFSqkl83UX8iw068J/\nHwv8AnBF+PcljGCTuipjFtKe+jPFTBFLRMXF8YZFgCQoY5mxNOLMqsIbbc/qaB61U4hwcbn18n7a\nRFifLTZPBb4VWkMQkY8BLwBu0Pb5TeDjSqlbAZRSd9Q8pjUw1+xKyWCay3kqDjOZrooymcLu1Uac\nhylqWjQy3SJmddrx2dVwcbZYadIE1xbiwnANzXphChvXbs+9YIoZG65iJrbOGJYNWzzQ+Oxq6DKM\nUu2TtYeCZo/d1psWjTio2sUOYkuPz8QQMtARotFn3jo6NQyxM+8GrlRKTYc9nbrKDIYPWW8F/rGX\nC2UKGqXUjvBinwfOj8y/InIp4JKo3xdO565CrhKXAFgzNibKYLLtH9VFyaTo4p3W6mBYqgAPun1B\nirDJtdZkCCLoFje2fZzOb1zDdj6bqImoQH5OiMi12t9607iHA9/RXrsNeJpx/GOAE0XkcwRi491K\nqQ/2Piw7ozLXlMNhOdOihBbZnvlkvqwtWEmMQGANPXOnU/q+M664I7RBJIKabIVt6W0CzPFB9cJm\nclO3qClijelyMbm4e3SMzLUYPQtK68XVOjoViBTaLGBxc1mu3yVqItHSNP7WMMUMkMzSGmC8jIis\nB54FzAIopdJSo14N/C3Bw0xpXIOCJ41B3EM/aydXTJk4kSLWncTi59quwIZL+4OiwqLoWAYtXPIw\n3o9CrjYzi2x7D5aXEq0qqnQL3n3ievaf+SxtyydWy5YPDzkBeApBx+2TgAMi8sVee604sKbmGjd2\nwfwc7b2BlWZips2O/cFibS7CupBIpuLOxe4mk0Q68uyUZm0Iqg2PHzmeRc5JlsjXrT97SVhq8jBF\nTyRs9IXbaokwF2uLqNPfD5vFJrPbt3F+W9Vjs8BeRMe6RXeI+l666w4dgPbJ93PG4uEgZb8B7WdP\nhLE2IbeFx7wjEJRRzI21rk6OkInGnhAy0djqJa/b9jkEgWd/IyLnAtcBrwmrBwMgIg8HXgQ8hz4J\nmg8CXxaRy8O/Xwh0/88ZEKfdeyh1IXJdJPLcFuAWG9MlZvR/dd5n/F2wzH9poTHsAqUMGcIh9bPa\nafm9jPjMez8LWNf66Bb8LvAI7e+fCbfp3AbcFU48R0LLyblA3YJmqOea2ghdTodbZ7Cyv9vyEMVG\n6E/gidiOxcl4u1OV2yiuJxRBy6F7xLRWxE/4B3bR3rtA66oMF5fZMoCkmypCX7RN8hZxnS6Vm7G/\n7XoJcWg2rnTMcAK6K0THv4dCdT783URzTzZljtaRKRibtwpAG+b4rbVzKm6PcTcPYZlfBuAGbgKu\nz+u2fQJwPvBqpdSXROTdwOuBP9H2eRfwOqXUA0b9qcK4Zjm9SUQ+BfxiuOmlSqmDPV25T2QtEpF4\n0ffJc1tYBUuRJ3tTyOjbi4iaXqw0a5Eq3wPXQoa9uBEHyzXAo0XkHAIh82KCmBmdvwd2isgJwIMI\nXFLvrHtgozzXlCeoXdMYW05YJSY3ATM5Yia0niQK96VUpU1YXNgFLCSsDaZlpn1RJJiiOjsrtHen\nG8tas90ioEWwTRdZabVYIFlXJpdp911dWjUUEjEmkagp7d7ZRfvkObYfWQQtfT7CTNePsLqYahIz\nJbkNuE0p9aXw7yUCQaNzAfCxUMxMAA0RuU8p9XdFL+bay+ksYBW4XN8WBQwOBWmLSwGTv05a5lLq\ntatYqIqKmogeUp+PaWxNMOt+HwtWfa4DpdR9IrId+AeCtO0PKKW+oXfAVUrdKKGrAHMAACAASURB\nVCJXAl8DHiBI7f563WMb+rnm5J+GJ1XcGDJsbaAv8jYxk+dOaIwt09iWrIHSJYLmV4gzXg7sgC0L\nQal/jAXRVlXYrLFj0N6rWYwMt0yUiRUJL5vQ0ccL5Vsp5JHaXVuPM7KQcDulsdEeOOzWcyvIkJrd\nuA+2kCiqmCZkgO7vBRQQVnPApY77Fkcp9X0R+Y6IPFYpdROBC/sGY59zot9FZDfwyTJiBtxdTsuE\nlTsJ/OnnADcBTyhzURdc6mQ4UUBsZAVpQg8ugTSrTNq+ZUTNsFLHYj0I8WZzHY2wkFRKtSC5Yugd\ncMO/3wa8rZ/jos9zTWXzTGkM59JMJ5XXFDPJMvzdBO6mlj2m4qIJuwjbawgRF6FmpqXbiBZwI3uq\nvXciIXaabE3tSN0i231WVvBkVQOO3WYpwiZT1OiuJ6c0/RQO7IAD0GQhEYhtI1XkxsXzbGHTVSa3\nO/Fq4MNhhtO3gZfqD09VXsjV5fTz+t8icj6phs3ecamTUYi8QNoyi24/FzLTkjA8rotC4yga9OrU\neXzQgmKERc0w0s+5pqd5xpIpVBytSzdabMk03UIG7EIiXDjbuycYn12NewZ1uSDyhIqLkDmwYrE0\nZB1jpquH49VF2Ra6BIIucHRhYbPm9CJq9PN0vZ4jbPrCXmBbwft0ElULyT9r6djeQSn1VQK3ko71\ni6OUmu3lWq4WGvOiXxERM82zSlzqZGRjEwFpZAiEgbQB0K00Raw7vVChQMoSLoVqreSkwm+6/RPZ\n465LaJjfl4reuzo6ko86Nc81vc8zvRK6m3aOzadnq2RZRMLWCkDcXiEutOcSS1Eo5kOvx+si5Iyu\n5Ob1LOIGkgLHFDdVipoIs15PadKsZ2Xjag7soHV0PnHPad3BEzFUtkYhg69F0xdcY2j0Kp7HEUQt\nf6+WEQXk1skQkVcArwA46+HU9qTsJGaGxVpShhJjd7W02MRLWrVmG59gU2a6/IfOzBGbb0oJ4tYp\n697TRE2VVZiPdVHT57nGpR5PYq7hQWdVO4It3QGzuVYZk/i1yZTtLhSxNOVYZaz7ZYibyJIQ3au2\nIEeWp35gs9okPptQZFrdTpmCobwVr33RBFwd/B4FBnfV9xnrjC/tvYp7idmo2ULTT1wtNOu03+8j\n8HP/bfXDcSfMdX8/wAXnSuBz10XNK0laN7IEj4N1JnNhrNMFVGbB1e877/gSQaouKe42ioiZtGNM\ngROJpkLCU/+O9ILj5+5qlYruIXhfi4mauzjduM7IiqKhnmvkZy9QObsXQg8Gzqz8GweXuggbG5HQ\nsDULiF6vLsg5GI/jOZ1bQ5SnSJG/5TSrj1aDR8/oiq0iNlHTa1G7AztoyhzNjVsTtWqg+54aYxmW\npm3d7rPcAGcLP+Ah8flv4njg+sLnqBNXQXODUiqxkojIJqC6XgJJXOpkxNx94nqI+m2bokYnQ9SY\ni3Oqq8nVClSXu6OoJapMkHGBarcRNnFjZoqlNhAtQFSx2YXcGJyKPqO0VP+irQxGoct5H+jnXFNo\nnomJFq9e42hSCuGBZgXQ3QdmAbdShPETphsqvpcIe3uG7gXaIpAOzFnGmROEqouaMBXdhnOdHQ3X\npplR3Ze0/RMdtHXxEImF2WC/pFBI6XBeiKClRvMZCwlrTVSzx0x1T0vxbowlg6EjUdbdrX10cRU0\nb6B7QrFtqwqXOhkxp917KLnBXPRLxDtUtrBUFQOTdw9Z1qleMqccrVppVW/7IWoqy0YrSRVixguZ\nmH7ONYXmGYDjTr+PB6q6urZomy4Om/sg021Q+NqTsHchmZ69OKm5UxaS6dsbJ4PxbglERydmYwHm\n8xofOmTUZAg10+1kipoycTS2SsUufZTM4oNxTIsmFhLWGyaqc+lE1prFrezetjnuE6Wj30NuTZ/Q\n4uTS/3JUyOu2fTHQAB4uIu/RXjqFwBxcC2l1MgqdpMbsm9qxWZY0zPHGC3hdosaGpSeSOaYsUQPl\nXFD68WnXtVKgYq8TFbkYe+obtYYYxFxTdp6JmxVWFGhpW5DNp2mgK9erlwycqA0CzHdSug9Ae2+2\nBSgSFrGg2BamFydEjW6l0S1CISWsTHWImrKYwqdLOBgxLe35Ct157IK9C2yfXYTppKUmb5z6WPXA\n4roRkVuAw8D9wH22liwi8myCisEnAqtKqV8qc608C833gGuB5xP0YIg4DPxBmQu6YquTsWYoWRHY\nqZt0lsutDjKEgh5rEy3cVQobZ2wxT2XJEZdFGKRlybX+ioj8AsEz5ouVUks1Dmkgc03ReWb9cYe4\nN/pjCwQuheoWrCKxHpHboyyxaLo67Bg9v2KJZ9ECeEPrTCQkYgERixqzo7hWiXhLeuBvvD3FzRRh\nEzX6eGIrV0on8ELNMvPMFlqbhTTXD2OhqDlyPO2T3S+dy4EV2hdNsv2qRabGzsnfXyNyrenZUn2y\n0DxHKWWNWBaRU4HLgOcppW4VkYeWvUhet+3rgetF5MNKqdosMj1zB7m1ZYbGClORmDHjLRJtHAZV\npyXHagPUKmxKCQJzrFnv19oRM071V8L93gr8Y91jGpm5ho51o0xQpUlkaWjRYJHtweLYBGaST9NV\nd642aWxr0dpmaa2gF8jbG7hQWrNJ68gi21nctj1hQTrn6M20Lwqq3Ua0jk4l3GaZDR8tpBbgs1hr\n4vtKed+6ulrruKzw+j6huLEKm8rdOh333U5Li4TUCxkCTBc1Q8BvAh+PqoErpe4oe6I8l1NTKTUD\nHBSRruh+pdSTyl64dmpKq3XGjOMpai3JETP6v6nCxjhPgjqFjpH9YwuaNd1Q4CZsTHdTbuxJ0fu0\nBQxXnMGWK2bqF6Gu9VdeTZBh1FMHXBdGZa5Zz6FgkRwj2bm6TJyEJhAAlsda1rgIXcy4xHm4LGpA\n5qIWiw3TPRSOWRcvpgVkYqnNzdPnMH/1zqSgGCNxXFaacVFMa028XbPaOAvDGYopkGhfTdhE7228\ny8Z99nYSJRm/KuV9cxy7bYwlyeu2DUHl76tE5H7gfZbXHwOcKCKfI8hyfLdS6oNlBpPncnpN+O+v\nlzn5KJOacVI0y6hsb6aQNDFjbjPHqh9nffrvc4VbV1EDSdESiRvXzKbKyBExQ2PxS2dCRK7V/tYn\nGpc6Tw8HXgQ8hz4IGkZkrlkfZlNGooZtobAp41I4sANYiONxZrfsY2pbi4lmIFqKWhBW8v6L7A97\nRIVM0O1+6IqpMNORw1ib5vzzaG7cyvhVq7GI2H50kcOtM+KxL06n/yeKisK5kiZYbPtkoYuaKLjX\n+j4XFTUQW9b0c0bWmpmr99DsrZF0SHfbgi5rU06X8gg9s8uVHzxwWvw+33PPdwDyum0DPFMp9d3Q\nlfRpEfl3pdTntddPAJ5C0OfpJOCAiHxRKfXNQoMj3+V0e/jrNqXU6/TXROStwOu6jxouXKrWFs4w\nScmi6rkoWs4impU5Y75ms9qAIW4qTl/uOq8DNheUTpaQKWWdMWv0lKghVKWrqSoO/eQ007K1agu+\nK8C7gNcppR4Iu+DWyijNNXodkCjDZZYN5U6mVfplfoUzOMydzXXJBcqysNrES14Uz5x23OSm4Lym\nqMkfq3alA0GWU5z+e9EEm6/eHbvOzAXTZg3Qg3ltriJbgbs0oqq5efualppUYRMJgx6tNQD7mKW5\nqBwywtww7zHThZZBmWOKopT6bvjvHSJyOYGFWBc0twF3KaWOAEdE5PPAuUC1gkbjuXRPKBdbtg0l\nRVNo9eNSF6AUK0pXLEueWMhYTPXz6vegL1xpC36awOlySaWNoQbrTVZBvjxhk7a/ee5CRJlfNRRG\ntN3PEKVmu9RfuQD4WChmJoCGiNxXtgtuAYZ6rjnh7iBpO1oIIhfRbNmOxQmXzi6YJxA129ZZF5sy\nQkbfz5ZAbQoPa9fptPYJe4HZTu2V5jO20ri64zqL7sG0BOhWofi6FgtRXs8lkzQxY2ZB2QKFKxU2\nJk2Yae6hWWnGUwGL04AQkZOB45RSh8PffxX4U2O3vwd2isgJwIMIrMXvLHO9vBiaVxF4PR8pIl/T\nXloHfKHMBfvG9mDhdFkg8wqalVksY2uNTdT0sHia8SW2eBObyMkMIrbRS2BxjkBIK0QHyUW/iBAt\nbR0zqknnvi8WXITakJFbf0UpFadPiMhu4JN1iplRmmviBSOxwJVYqDYudFw6cXzFnHNcSaVtikPc\nUnk1WRRm3ES/wy5aNJiaTqYT29wbUxaB0aDF9qOLMNYd0JsW4GuON81KY0vt1q1SpjCIxh2ju3Ly\nxI3mfoqObdCidWSK9snViJrofsx7SHMjDUjoTAKXhw9GJwAfUUpdqXfbVkrdKCJXAl8DHiDIuvx6\nmYvlWWg+AnwK+Avg9dr2w0qpH5S5YF8wFlOX6rJlrThlxuOKufC7ZgGltQvICiLWr9lFGdeUg6hJ\nG0NEmtUmTSTU3QfJNs4qrtfv/k1p9Vf0SaZvg+kwWnONsaAVS82d6xSqS2FZFwVl4jlKYLqdrD2L\nIjGjW5aMVO/mZVtpbLMHOGddG4L7tmbvZGBacTLjbFJSu6Nrm+OpRNxo52mMLdPcuLW3YnuLnffa\nvB+b2zD3XmokTDw417J9l/H324C39Xq9vBiaQwQ9BX4DIAzqeTAwLiLjUZrVKJDXEyiLsgtOVQtV\n1WIryyKVOeaiVpuCLRRcrDaFsFnHzIKDjrjEzbiINCe2A28uf7gLtvoraUJGKTVb72hGf665eewc\nzth42C2TZeNknKUStzfQ0qP1dghmbMTkpo7bKa07k428Wr1dWS9jWmG4vZGoSTmL2RJhfoXts4td\nVXV10mrD2OJsiqQWm0IlzdqUl9pdWty4oH/eRdm4kEiFj0i7T9OCA0krzjC5p6rAtdv2JcA7gJ8m\nqPrys8CNwBPqG9rwMAyVW3up0VJLplDRLCmHGjVlm15mUmScO+FDbwrF45nZn3ved6KWezkGGNW5\nZmKpXWqhilsIhLEVejpuInhVs9JEmUqmsMljMmcK0TNzth9dDMRWlOHk0hxTG037ogn2XD3j1Bcp\nj6zU4rxzuAgc2zZbvI0+lkFYOyBI1y4SV2SLSdLFZHAfa0fUuAYF/znwdOAqpdR5IvIcnEohDYAe\ngzujVOIiTRBHkYE1QnSIJ6pdDESp9CmWmrz3par4Go+V0ZlrdJowvnvVodHfHGzpDl5tLm6FvQs0\nxvZ0LcKxlcZwc+QJlFTSLAqR62Q6tBzNG1WDo8rILhzYwexl+9i9bXOmqOmVohYdmzUmS+QUjblx\nGWtrdqpcQ8iNC+wc25xwz9linmzp67ZqymUKNj6wcmLHDXlNleWPq8FV0NyrlLpLRI4TkeOUUp8V\nkXfVOrIilC6U3EFfxF7Ch8jNwnSwOFRFNJ4qWwSkLdp9t0ZlvI+VunBsVpocYWO7btn3xwsbZ4Z6\nrrnvtOOAB2KLiZ511Ggu54dTbJxMuAwatIJFZRuwLai6G6Evnl0ugqJujjy0gU8stZOdvcvKydD1\nVLQ8f6/Ygo2zyLLiZPWIsrly0khYv3JHZDLHjLqSBklxmFY/SG9wGt+H9ntUmbrD2mlP6Spofigi\n4wS54x8WkTuAI/UNqwQFytK7VJuNtlkXsLIpzXqgrGM6d1ZG0ECpOq07J96mVqEVCZvt2QHTCUoK\nWpuwqeTevg+8vffTDAFDPdccYj2r0/fGC5Qe07JvaZbWkTtpn7yYfoItyQU0WqDitgchaVVcixZC\nyyLTVWIW1Esj0wW1i/ZFCyxfXX/rhjTKuKvKNrlMEzhmtldj2kH4aowfOb4jfA1s6e9ptXsiUaMX\nKQz2PfYEzQuAnxA0ifstYD3dueRDhasQKJyunbGQZy5MO41/89AWzE23fwLOpMtK4+ISK2LVGYZY\nIaDWtgOpbC9pPdGFje2zzehG7rEy1HPNIdazzFODDCTaybWgSc5iNcfMtj0JlwEYaeAVWF7SFmpz\nQSxbjC0/jkbf1x4kbV43T6iVGWfeOXsNQI7Ick8lrhF+MRrTLZquhRgXk24m6K7tY1qk0ixONoFT\npFrzKOAkaMIKfhHdIdbDQMoTs+4+MuNiMq0wWdcwyD2+F2uGkf5cNK7nkg37u0SNLuJKLaz/f3vn\nH2xHWd7xz7cpP00U5dI0iJYIWIsMRUoRkaFO64+QKNBOuEFbm1iUCUlaO4xVKB0T2rGDSktliEZK\nI+nIEG5SnWYgIwMVCmNASPglPwykhLYIhgRbyE0cJOnTP9733Ozdu+ecPefuOXt27/OZOZM9u3t2\n3+eezfN+z/u87/OUUfAyTb+ETl6atbMq7R8QquBrGh1EQ9QkZ5XMZWOo2wMTktFN3zONW9YvApqM\nwLQRM5MtJJg1/yMX6ZGa9JLtCUxcd9VOkDRWWRU5wbbdtbIET7cjOhuZ23rlVDeDIO9bzo1LFrBw\nffaHs0Rhnu+1k+zLRRCL3G4GfmJmH00dexPhCXs7QY9cbWbf6vZe7RLr7SYUlppwCDAze2O3Ny6U\nZG3ORAdywbLYWc+KOxKiOLeY6aD6cs+4Ltoyq7NyDY1zs0RQ1yME/RQzrSiqHR18h2N/syLuPSh/\nxwGhKr5mrDgliVwx6VUiWeGa9y1n++Ezxjq2RskByD9SkleEtBM+ba8zAjeOLGDRt285sO+Pwoqs\ncVWt9+apYbWYeRzf7iSg/0uIW90vr9i5rcnoTBY71nUykXsxw5vWNBUzWTSzZ1z1byauMOtDwOmz\nhJWKWf+HlwJPmtnHJB0FbJV0k5n9opsbtctDM6Obiw4MsdOYIGwihc6PadOGIq7TEDUNWomSdbM+\nlj/9ft5QST864XStpT7RLtxUqJgZACTNAb5GSKx3g5ldlTr+h4RyAwJ2A5eY2aO9ak9VfE2jOGWS\npCCZx8aQu2Vpch7NYnZumsHQ8PgOp1WHmrdz73SUIe89Fq4fYdmdK8dWtAwvWTM2j6MxmnLb4RtZ\ntPKWmOk4Oww1fc+0cJ90r9ltaK2o3rfN/fOKnU4mIOcTM4vH8hStZFkh9g6NjE6wt5H0sMg5WVlI\nOgaYB3wJuDTjFANmKKQSng78DNjX9f3Msn4UVYvTZsk2/0mbk3rVQXeaj6WLpG7jOvYORhTStaAm\njErlafdkSiBk0Y39SYoWOTnnzvRb0Ohv2dJJcUkNnWacmyiu/S01/XwcAn6aUDfpeUIphI+b2ZOJ\nc84EnjKz/5F0DrDCzN6bdb2pxGnHyTZ/OWw3S062a/50boshiDERMNxCoLTr3HvdiSeunyxeuWuk\nda6VpJ0jZy6cWLwS2GlXd277IM5RzSHA2pYciHOkFsy/kRE9l3HmYlgZVsGNTRJv8bdoW109gyxR\npXX5fI2kw5lx1h5OujfseOEa+M9LVzM+E/D1ZnZ94jPrCdm/ZwCfywg5zQA2AO+K5ywws65jYXkn\nBQ82O2neUTY6wMl0ROlrJzvVfo5aNAoptiNjhVRXYqYIJitgJnu9Po7yZFK0/cVwOrAtpiVH0lrC\nZNwxQWNmmxLn308oYDnlef1ZDoSNRrI76qGRURYywsLhkbFzW3Y+/Up31cF9dqyDmc2ESKKS9EJG\nQkhjEylRA6xcztD6FQeulzVCkapOPRk67eA7yuPTTmQNZwvbrLpfI19fCKQnS8eRmUW7moqZZvZ1\nVKekw2SMObjLzC7KOiDpo8BLZrZF0geafP4jwCPA7wLHAXdIutfMXu2mMfUQNK3oJoQxmJ3QgerQ\n7UivkGp2PA9tajKNkbdt/Sbru0wJ0guWtV4Rl1sElvncvAyMn0o3JCkxZDPul9Nbgf9OHHueUOG2\nGRcR6ixNeXYCV+btFFIdUBkFsibD4lT7xwRAqqMdYpSV85ex8c55oVDlfaFzvnHJAhg+0BEnO+QJ\nYiKOXuyaP50hRtlRXMqtpnQzwtGMrHSDE+ZWEUdx1BB9iSco1vaae/ht48TMZKqrt6JPz+L7gXMl\nzSWUMXmjpG+bWTKz0aeAqyyEirZJ2k4YrXmgmxvWX9AkKarD6bbzLuL+ydGadmSt/Cp6ZOabqe0c\nyepKJ/03bCFqcomZwbR1Vychq2bETL0XAWdNvkn1omoCpVPS9qUFDkRhMhJGpbaPzOa2TQfmkSwc\nHmkqGjLFzQgsnX8dt7BoUu0shXUTBW5a5OxYFyashRYvHpe4MDlHKXl+4+wqYmaXA5cDxBGaz6XE\nDMB/Ab8H3CtpJvDrwLPd3nNqCZoi6WTkpxejF63CYGmKzmrcrgMfzA5+IqnvJRmiq3memJ8Ab0u8\nPybuG4ekk4EbgHPM7OU+tc0ZULI61nEiR6PMiUMLyYSD7UiGo7JWC1WlQ88jAOORcdXWx82ZgbEM\n1FWxu1MkLYaxYrh/A9wo6UeEBQhfMLNdrT7fitoLmitTbnj5kT34fJ4OvNedfF6B0+kk5gaNz1RF\nrOQlQ9S0Pb/6PAicIGk2QchcCHwieYKktwPfAT5pZk/3v4lOFZjQiRPFyTCQo1NOnr9m/nBYpjyS\nfe2qkWx/VniysQw+PWemjmLGzO4G7o7bqxL7XwA+XNR9aito0kIkuX/5kQeOpwVKs8/1gmZtKIQs\ngZMlZNrNMaHJ5/LetyrknStUA8xsn6RlwO2EUfDVZvZE6pfTF4Ejga+HFZXsKyKE5dSbVYSRiZnD\nHXTKceVPY5lyHTv0GGQKJObLrGTZuFVgdbS9n9RC0LywvzMhkjy3nwKmGck29ETcwOSFRlWFShEU\nNO8o63su6/kzs40wfnw/9cvp08Cn+90up/qsAsg5qXfmOpgxdyfXsZSh9aO17tBXEbLIJefMDA2P\njgvN1dX2flGKoJG0AvgMYdEAwF9GB4ukywmTEPcDf2Zmt5fRxilNen5Q3cVMwXOc8oiUQRDSUwH3\nNYPNjgtgu80e69jr3qG3Sr5XCdt3v55Ymr+71KZkUeYIzTVmNq4+sKQTCfH8dwNHA3dKeqeZ7S+j\ngWmKHj1Jhr4GlrqLmSQdrgob+O/OaVA5XzNVWAWg/pY8KIvFwPF7D2SQnsfGKSHi+skvld2AFOcB\na83sNTPbDmwjJAGrLVlzeLyjHBDahJp6Fh7sihcIyboaL6cNU87XOOWyCrju8KUAY3NnXMwUS5mC\n5k8lPSZptaQ3x31ZCb/emvVhSRdL2ixp894eNrLRafWy88q69iCLmikhunLOmxksUeM0oRK+xqk/\nczTC8JI1DA2PjiVorDuSpkl6WNKtGcck6VpJ2+L/0VMnc6+ehZwk3Qn8asahK4BvENafW/z374B2\n1ZjGEbOeXg9wtNSTglT9EDOtaKzISr7PopP2tZuAnGflVRF/j7yCqG9/++Q8mmZiZiqF3ypEHXyN\nMzVYRRidKTJLcQVoVW37HOCE+Hov4f9r1zXjeiZozOyDec6T9I9AQ7nlSvjVL9KCopc0m0/Ty5GQ\npH3p+/Rq5VWn9vTzO8icHNxGxNR+pKoC1MHXOFOHlRotspbSQJOj2vZ5wD/H0gf3SzpC0iwze7Gb\n+5UScpI0K/H294HH4/YG4EJJh8SkXyfQZU2HIuj3yEyz+/V6tKQVvZgIPdBhmnQphzYMtC1OZXyN\nM7WYeUGhBSL7xWvwCrAnvn0M4ORGODa+Lk595h+AzwP/1+SauUO/eShrldNXJJ1CGAZ+jvi7OCb3\nGiFU/d0HLM2z6uDoabD8iObHk2GUdr+oy+6guhE1g3T9ybZjIOgwtNQsOWOnuWaS11nR45EfSXOA\nrxES691gZleljisenwvsBRaZ2UO9bVVPKNTXOE4RVHH+jJntD25jMaHk0hYIpQoyw7A5q20XSimC\nxsw+2eLYlwjDU4WR7CjSc1IGumN1Kkn6mRq0Z0zSNGAl8CHCL6IHJW0wsycTpxUa2y6Lfvsax6k3\nt78BfmVPcBu/g9n2VnPK8lTbLjT0O2jLtvvKoHU0jtMnTge2mdmzZvYLYC0hlp1kLLZtZvcDR6TC\nN47jTDHMbC+8B/g5sKalfjCzy83sGDM7lpDz6fsZ1bY3AH8cVzudAbzS7fwZqEnpgy37GdXLbC27\nHQUxBHRdbXSAqIsdUI4tv9bZ6S/eDiuGEjsOlbQ58f76uFoHsuPW6dGXZrHtrp1NHXgRRldQC1/j\n/z8Hj7Ls6MjXmH1Pk7lZqmbcRkJYexshtP2pyVy7FoIG2FqXwnmSNtfBlrrYAdWwxczmlN2GKUIt\nfE0Vnum81MWWutiRRYtq20YscVUEUzrk5DhTlDxxa1/W7DhOpXBB4zhTjweBEyTNlnQwIb69IXVO\nobFtx3GcXlOXkNP17U+pDHWxpS52QL1swcz2SVoG3E5Ytr06LmPuWWy7RtTlWaiLHVAfW+piR2mo\nyRJyx3Ecx3GcyuAhJ8dxHMdxKo8LGsdxHMdxKk/lBY2kOZK2xvLjl5XdnlZIWi3pJUmPJ/a9RdId\nkp6J/745cezyaNdWSR8pp9UTkfQ2SXdJelLSE5I+G/dX0ZZDJT0g6dFoy5Vxf+VscXpHlfwMuK8Z\nNFvcz/QJM6vsizCh8T+AdwAHA48CJ5bdrhbtPRs4FXg8se8rwGVx+zLgy3H7xGjPIcDsaOe0sm2I\nbZsFnBq3ZwBPx/ZW0RYB0+P2QcAPgTOqaIu/evaMVMrPxDa7rxkgW9zP9OdV9RGaPCncBwYzuwf4\nWWr3ecCauL0GOD+xf62ZvWZm2wmrTU7vS0PbYGYvWixUaGa7gacIWWSraIuZ2Wh8e1B8GRW0xekZ\nlfIz4L5m0GxxP9Mfqi5oCi09XhIz7UB+j58CM+N2JWyTdCyhuMcPqagtkqZJegR4CbjDzCpri9MT\n6vKdV/qZrrqvcT/Te6ouaGqFhbHGyqyjlzQd+Bfgz83s1eSxKtliZvvN7BRCNtzTJZ2UOl4ZWxwn\nD1V7puvga9zP9J6qC5o6pGff0ahiHP99Ke4faNskHURwMDeZ2Xfi7kra0sDM/he4C5hDxW1xCqUu\n33kln+m6+Rr3M72j6oImTwr3QWcDsDBuLwT+NbH/QkmHSJoNnAA8UEL751nCrAAAA0JJREFUJiBJ\nwD8BT5nZ3ycOVdGWoyQdEbcPAz4E/JgK2uL0jDr4GajgM10XX+N+pk+UPSt5si9CevanCbPAryi7\nPW3aejPwIvA6ISZ6EXAk8G/AM8CdwFsS518R7doKnFN2+xPtOoswNPoY8Eh8za2oLScDD0dbHge+\nGPdXzhZ/9fQ5qYyfie11XzNAtrif6c/LSx84juM4jlN5qh5ychzHcRzHcUHjOI7jOE71cUHjOI7j\nOE7lcUHjOI7jOE7lcUHjOI7jOE7lcUFTQySNtj+r42ue26gyLOl8SSd2cY27JZ1WdNscxykH9zXO\nIOGCxsmFmW0ws6vi2/MJ1WAdx3EKxX2N0y0uaGqMAl+V9LikH0laEPd/IP6CWS/px5Juihk5kTQ3\n7tsi6VpJt8b9iyRdJ+lM4Fzgq5IekXRc8teQpCFJz8XtwyStlfSUpO8ChyXa9mFJ90l6SNK6WKvF\ncZwK4r7GGQR+uewGOD3lD4BTgN8EhoAHJd0Tj70HeDfwAvAD4P2SNgPfBM42s+2Sbk5f0Mw2SdoA\n3Gpm6wGif8riEmCvmf2GpJOBh+L5Q8BfAR80sz2SvgBcCvx1EUY7jtN33Nc4peOCpt6cBdxsZvsJ\nRdD+Hfht4FXgATN7HkChpP2xwCjwrJltj5+/Gbh4Evc/G7gWwMwek/RY3H8GYRj5B9FBHQzcN4n7\nOI5TLu5rnNJxQTN1eS2xvZ/JPQv7OBC+PDTH+QLuMLOPT+KejuNUA/c1Tl/wOTT15l5ggaRpko4i\n/IppVbF1K/AOScfG9wuanLcbmJF4/xzwW3F7fmL/PcAnACSdRCjQBnA/Ydj5+HjsDZLemcMex3EG\nE/c1Tum4oKk33yVUd30U+D7weTP7abOTzeznwBLge5K2EJzJKxmnrgX+QtLDko4DrgYukfQwIX7e\n4BvAdElPEWLWW+J9dgKLgJvj0PB9wLsmY6jjOKXivsYpHa+27YxD0nQzG40rEVYCz5jZNWW3y3Gc\neuG+xikaH6Fx0nwmTtx7AngTYSWC4zhO0bivcQrFR2gcx3Ecx6k8PkLjOI7jOE7lcUHjOI7jOE7l\ncUHjOI7jOE7lcUHjOI7jOE7lcUHjOI7jOE7l+X9uSuuvJ8AWWAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fe228c87198>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.figure(figsize=(8,6))\n",
"plt.subplot(321)\n",
"dp=dclim.dTse0.plot.contourf(levels=11,cmap=\"jet\",vmin = 0.0, vmax = 0.6, add_colorbar=0)\n",
"#,vmin=-2,vmax=2,add_colorbar=0)\n",
"plt.title(r'$T^{\\,\\prime}_s$ for $\\bar U,\\bar C$')\n",
"plt.ylabel('latitude');plt.xlabel('longitude')\n",
"cbar = plt.colorbar(dp)\n",
"plt.subplot(322)\n",
"dp=dclim.dTse1.plot.contourf(levels=11,cmap=\"jet\",vmin = 0.0, vmax = 0.6, add_colorbar=0)\n",
"#,vmin=-2,vmax=2,add_colorbar=0)\n",
"plt.title(r'$T^{\\,\\prime}_s$ for $\\bar U(x,y), \\bar C$')\n",
"plt.ylabel('latitude');plt.xlabel('longitude')\n",
"cbar = plt.colorbar(dp)\n",
"plt.subplot(323)\n",
"dp=dclim.dTse2.plot.contourf(levels=11,cmap=\"jet\",vmin = 0.0, vmax = 0.6, add_colorbar=0)\n",
"#,vmin=-2,vmax=2,add_colorbar=0)\n",
"plt.title(r'$T^{\\,\\prime}_s$ for $\\bar U, \\bar C(x,y)$')\n",
"plt.ylabel('latitude');plt.xlabel('longitude')\n",
"cbar = plt.colorbar(dp)\n",
"plt.subplot(324)\n",
"dp=dclim.dTse.plot.contourf(levels=11,cmap=\"jet\",vmin = 0.0, vmax = 0.6, add_colorbar=0)\n",
"#,vmin=-2,vmax=2,add_colorbar=0)\n",
"plt.title(r'$T^{\\,\\prime}_s$ for $\\bar U(x,y), \\bar C(x,y)$')\n",
"plt.ylabel('latitude');plt.xlabel('longitude')\n",
"cbar = plt.colorbar(dp)\n",
"plt.subplot(325)\n",
"dp=dclim.cld.plot.contourf(levels=11,cmap=\"jet\",vmin = 0.0, vmax = 1.0, add_colorbar=0)\n",
"#,vmin=-2,vmax=2,add_colorbar=0)\n",
"plt.title(r'$\\bar C(x,y)$')\n",
"plt.ylabel('latitude');plt.xlabel('longitude')\n",
"cbar = plt.colorbar(dp)\n",
"plt.subplot(326)\n",
"dp=dclim.wnsp.plot.contourf(levels=11,cmap=\"jet\",vmin = 4.0, vmax = 8.0, add_colorbar=0)\n",
"#,vmin=-2,vmax=2,add_colorbar=0)\n",
"plt.title(r'$\\bar U(x,y)$')\n",
"plt.ylabel('latitude');plt.xlabel('longitude')\n",
"cbar = plt.colorbar(dp)\n",
"plt.tight_layout()\n",
"plt.savefig('Tsp4.eps', format='eps', dpi=1000)\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [],
"source": [
"# Now, save the dQdf and dQdT terms for using in INC:\n",
"dQdT = ALH + ALW\n",
"dQdf = BLW\n",
"\n",
"# Define the new Dataset\n",
"dQ = xr.Dataset({'X': ('X', dclim.X), 'Y': ('Y', dclim.Y), 'dQdT': (['Y','X'], dQdT), 'dQdf': (['Y','X'], dQdf)})\n",
"dQ.X.attrs = dclim.X.attrs; dQ.Y.attrs = dclim.Y.attrs\n",
"dQ.dQdT.attrs = [('units','W/m^2/K')]\n",
"dQ.dQdf.attrs = [('units','W/m^2')]\n",
"dQ['ALH'] = ALH\n",
"dQ['ALW'] = ALW\n",
"dQ['BLW'] = BLW\n",
"dQ['dTse'] = dTse \n",
"dQ.to_netcdf('dQ.nc')"
]
}
],
"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.3"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment