{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Describing ''Non-Fickian'' diffusion \n", "\n", "Dealing with multiplicative processes requires employing new rules of calculs, and taking sides on the \"Ito-Stratonovich\" dilemma!" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### 1) Some Conventional Actors" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "ExecuteTime": { "end_time": "2019-02-06T03:18:24.893042Z", "start_time": "2019-02-06T03:18:24.888006Z" } }, "outputs": [ { "data": { "text/html": [ "\n", " \n", " \n", " \n", "
\n" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%%html\n", "\n", " \n", " \n", " \n", "
" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### 2) Less Conventional actors" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "ExecuteTime": { "end_time": "2019-02-06T03:18:25.041435Z", "start_time": "2019-02-06T03:18:24.898333Z" } }, "outputs": [ { "data": { "text/html": [ "\n", " \n", " \n", "
\n" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%%html\n", "\n", " \n", " \n", "
" ] }, { "cell_type": "markdown", "metadata": { "ExecuteTime": { "end_time": "2019-02-06T02:58:41.084079Z", "start_time": "2019-02-06T02:58:41.077825Z" } }, "source": [ "#### 3) Richardson's dispersion as a backward-Ito SDE\n", "###### Source:\n", "\n", "Richardson, 1926, Atmospheric diffusion shown on a distance-neighbour graph, http://www.phys.sinica.edu.tw/jctsai/Jia2016/Reference/%5B12%5D_Richardson_1926ProcRoySoc_Atmospheric_diffusion_shown_on_a_distance-neighbour_graph.pdf" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "ExecuteTime": { "end_time": "2019-02-06T03:18:25.123301Z", "start_time": "2019-02-06T03:18:25.052522Z" } }, "outputs": [ { "data": { "text/html": [ "\n", " \n", " \n", "
\n" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%%html\n", "\n", " \n", " \n", "
" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "##### Richardson's equation interpreted as a \"1D Ito diffusion\" in terms of the algebraic separation x\n", "\n", "The diffusivity is $$K = \\epsilon x^{4/3}$$\n", "\n", "##### ITO - SDE: $$ dx = \\sqrt{2 K} dW \\;\\; \\rightarrow \\;\\; \\partial_t q = \\partial_{xx}( K(x) q)$$\n", "\n", "##### Backward-ITO SDE : $$ dx = \\sqrt{2 K} \\diamond dW = K_x dt + \\sqrt{2 K} dW \\;\\; \\rightarrow \\;\\; \\partial_t q = \\partial_x( K(x)\\partial_x q)$$" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "ExecuteTime": { "end_time": "2019-02-06T03:18:29.081759Z", "start_time": "2019-02-06T03:18:25.125726Z" } }, "outputs": [ { "data": { "text/plain": [ "(0.1, 100.0)" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABHgAAAI4CAYAAAARel4VAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3Xl8VOW9x/HvycKOGQJoRauYWHdrCdFqta5JtSquibhVxUribgUl2Palvb29lUQFK9pKRGtpLQKxWtEuMl5tXWo1iVfQalWiFsVWIEyUsGSZ5/5x5mSWzGQhM3Nm+bxfL15nmTNzfiZpM/PN8/weyxgjAAAAAAAApK8ctwsAAAAAAADA0BDwAAAAAAAApDkCHgAAAAAAgDRHwAMAAAAAAJDmCHgAAAAAAADSXJ7bBQAAAGQyy7I8kkollUhqNsZ4XS4JAABkIItl0gEAABLHsqwqSV5jTItlWauMMeVu1wQAADIPI3gAAAB2gmVZJZJKjTH1IecqJPkkeSS1GGOanccD17e4UiwAAMh49OABAAAYJMuyyiTdIjvIcc4VSSo3xniNMQ2Bx0NNl1STvCoBAEA2IeABAAAYpEAfnVURp8tkj97pERi144zsuV1SYVIKBAAAWSelp2hNmDDBTJ482e0yAABAimhqatpojJnodh0xeCRtCjlulVQUGNlTK3t6VrOijOIJ9OmpkqTRo0dPPeCAAxJfLQAASAsDff+T0gHP5MmT1djY6HYZAAAgRViW9ZHbNQxWYLpWQz/X1Euql6TS0lLD+x8AAOAY6PsfpmgBAADEhy/iuFA0VQYAAElCwAMAABAfyyUVhxx7jDHNbhUDAACyS0pP0QIAAEhFgVW0yiV5LMtqDqyc5bMsa0XgMcnuuwMAAJAUBDwAAACDFFhFyxvjPAAAQNIxRQsAAAAAACDNEfAAAAAAAACkOQIeAACAFGBZ1jTLsurb2trcLgUAAKQhAh4AAIAUYIxZaYypKigocLsUAACQhgh4AAAAAAAA0hwBDwAAAAAAQJoj4AEAAAAAAEhzBDwAAAAAAABpjoAHAAAAAAAgzRHwAADgkpaWFtXU1MiyLNXV1am+vl719fWqrq5WQ0NDz3UNDQ0qLi5WcXGxvF6vJMnr9Wrq1KkaN25c2LV1dXUaN26c6uvro75+XV2dampqwp4TTXNzc89za2pq1NLS0nNP57WKi4tVXl7e87rFxcVqaWnpeQ2fz6e6ujp5vV55vV7V19f3e18AAADsHMsY43YNMZWWlprGxka3ywAAIGF8Pp/GjRunyN/HlmVp1apVKisrk2SHPLfffruampp6rvF6vaqsrNTmzZvDntvQ0KCKioo+X7+mpkbjx4/XnDlzBlxbfX29zjvvPHk8HklSZWWlioqKVFtbK8kOhSSppKREPp9PJ510kp599tme653XaGpq0qJFiwb+RQr/ujQZY0p36slpgvc/AAAg1EDf/zCCBwCAFFRSUtITmEhSRUWFmpube42QKSws7BnVI9mhjxPu9KW2tla333572Ov1p7CwMCysiVZza2urJDv8qa6u7nV9VVWVGhsbw2qGzbKsaZZl1be1tbldCgAASEN5bhcAAEDCtLdL77yT3HsecIA0evSQXqK5uVk+n09VVVVh5ysqKrRo0aKeETOh55yRPj6fb8D3KSsrU0NDQ5+jeCLvP5DXlOygKbTOUKWlpWE1w2aMWSlpZWlp6cxBP9myYj+2aJHk/CzV10vV1X0VEdyfOlUKCRljXgcAAFICAQ8AACnAGdHi8/m0bNkyrVixotfol+rqalVXV6u2tlYtLS0qKSlRSUmJiouLe67pa4RNpMLCQr322mvx+Q8I4Yw8Kioqivp4cXGxmIKUXZyf75aWFvl8Ps2ZM0cNDQ2aOXOmysrKtGLFCvl8Pk2dOlUVFRWqra2V1+tVS0tLz8+pc666ulo1NTUqLCzs+d9KtNcPvbczUm3t2rUqLy9XS0uLqqqqot4DAIB0RcADAMhco0fboxDSQOholoqKCpWXl6u6ujps1ExZWZlaW1t7pmo5jxUVFamhoUEej2dQo2JaW1t1+OGHx+8/IsAJdpwQKtLatWtjhj/YSQMdUVNVFRzN05+Qfk9DVV1draamJpWVlam6urpnKmFra6vWrl0ryQ4na2trVVFRoZaWFtXW1mrVqlWS7J/Vuro6zZkzR2VlZT19nJyfo2iv7/xvIbRPVXFxsW655RaVlZX1eQ8AANIRPXgAAEhBlZWVqqmp6XX+vPPO07Jly8LOVVdXa9GiRYOaniUp7ENwPDlBU6w+O16vV9OnT4/7fZG6mpqa5PV61dDQoNbW1p4RNVVVVT0rq4UGgk5g6azAFtr/yePxaGoguHWuj/X6Tp8qh8fj6ekT1dc9AABIR4zgAQAgBa1atSrq6JfKykqVl5f3jHqQ7BE/NTU1UQOhWOrq6lRVVRX1HvGwYsUKTZ06VVVVVWEfmuvq6lRWVjagfj7IDM6KaitWrFBRUVGvaYFlZWWqr69XYWFh2M9FYWFhzAAyNLTp6/U9Ho88Ho98Pl/Pz2Ho6LG+7gEAQLoh4AEAxM/770v77ut2FWmjpaWlZ7nwurq6ng+gTU1NKiwsjLqUuBOOhH5ILSoqUkVFRa8PqtFe3+fzadOmTSouLu5zKkrkc0Pv2dzcLK/X2zNVrL6+XmVlZWE1eTweNTU16fbbb++ZBtba2qqioiKmwGQZZ3SM8/PhjDRzRpDV1NSovLw87Oe9oqJClZWVvV4nWhjT3+tPnz5dy5cvV2FhoZ599tmdugcAAOnAMim8CkJpaamhCSMApIl166Qzz4y96g4QB5ZlNRljSt2uI5Ey7f2Pz+fTzJkzNX369J4Qc8WKFaqsrOwJU5xphqGcELGkpEQ+n08VFRVqbm7WzJkzVVRUpFtuuaXnsb5ev7q6Wo2NjT1TsKqrq8NWe4u8BwAAqWag738IeAAA8VFaKvn9BDxIKAIeDIazStZ5550nj8ejlpYW1dTU9IRDAACkg4G+/6HJMgAgfrq73a4AAHqsWrVKZWVlYf13pk+f3tOEGQCATEIPHgBA/Kxe7XYFQNqyLGuapGn70scqbmpra3v6TxUWFvasoFU10KXiAQBIIwQ8AID4SOEpv0A6MMaslLSytLR0ptu1ZBKaegMAsgVTtAAAAAAAANJc3EfwWJblkfSBpNDJzcuMMXWWZVVI8knySGoxxtCJEwAyhWVJU6a4XQUAAACQlRIxgqdU0j7GmKnGmKmSbg+EO0WSyo0xXmNMg6RbEnBvAIBbjJG2bpXuv9/tSgAAAICsE/eAJxDg+CTJsqwSBUfylMkevdMj8LgizlVZltVoWVbjhg0b4l0eACCRPv1U+uUv3a4CAAAAyDqJ7sFTFjINyyNpU8hjrZKKIp9gjKk3xpQaY0onTpyY4PIAAHH1+efSq6+6XQUAAACQdRIW8FiWVabwPjwAgEy2Y4fbFQAAAABZK5EjeCoVHvD4Ih4vFAEQAGSOt95yuwIAAAAgayUy4ClVeKizXFJxyLGHVbQAAAAwUFOnTpXX63W7DAAAUlLcl0mP0OrsGGN8lmWtCEzdkqTaBN8bAJAsfr/bFQDIAitWrFBRUa8WjgAAQAkMeAJLpEee408uAJCJCHgAJAHhDgAAsSV6FS0AQDbo7ra3//M/0uTJrpYCpCvLsqZZllXf1tY22Ce682+AGhoaNG7cOFVWVkqSfD6fiouLVVNTI6/XK6/Xq/r6etXV1fU8x+v1qri4WPX19WpoaFBlZaWam5s1depU1dfXh13X3/O9Xq8aGhpUXV0dVldzc7Pq6up6Hvf5fD3Pde5bU1MzuO8FACDrrFu3Trfddps2btzodikJn6IFAMgGfr/9gW/aNGncOLerAdKSMWalpJWlpaUz3a4lnioqKtTa2qq1a9dKkjwej2pra1VRUaHi4mI1NTWprKxM1dXV8nq9Kisr6/nX1NSkRYsWqaioSCUlJSorKwt77erq6pjPr6ioUFNTk6qqqiRJt99+u1paWlRUVCSfz6eZM2eqqalJklRXV6f6+npVVFSotrZWq1atkiS1traqrq5Oc+bMSeJXDACQThYsWKAFCxbogw8+0JIlS1ythYAHADB0zgpaublSDoNDgaQyxu0K+lVVVaXi4mLV1taqpaVFJSUlkqSmpqaepsmtra1qaQkusOrxeFRcbK/P4Vw/fvz4sNft6/mS3ZTZUVhY2DNKZ/ny5SotLe15zAlw6urq5PF4el6zsLBwiP/lAIBMtmnTpp6RpbNmzXK5GgIeAEA8nH/+oKdtAMguZWVlqq+vV2FhoSoqKuTz+XTSSSf1NE5+7bXXej2nr4BlIM/fGYWFhb1GCgEAEM19992n9vZ2nXLssfraYYe5XQ49eAAAQ/TOO9LatXa4w+gdADHU1NSotrZWHo9Hkt3rprCwsKdxcmgPnIEYyvPPO++8Xtd5vV5VVFSosbGx13kAACK1t7frnnvukSTV/PWv0rnnulwRAQ8AYKic6SEEPAD6UFRU1NMfR7JH9Hg8HjU0NMjr9aqysrKnJ05zc7O8Xq+WLVum5ubmnnPLli3TihUr1NLSMqDnO9d6vV41NjZq0aJF8vl88ng8WrFiRU+j54aGBpWWlqqoqEi1tbVhzZcZzQMAiOahhx7Spk2b9HVJx0nS44+7XJFkmRSet11aWmoi/4oCAEgx770n7befvf/++5LXK0WsVgPEi2VZTcaY0v6vTF+8/wEAILV1dnbqK1/5ij766CM9LuksSfriC2nMmITcb6Dvf/hTKwBgaHJzg/v04AEAAECGW7ZsmT766CMdIOkMSfrv/05YuDMYBDwAgKEJDXWYogUAAIAMZoxRbW2tJGmOAqHKtde6WVIP3okDAIbG7w/ub9/uXh0AAABAgv3hD3/Qm2++qT0kXeScDCwg4DYCHgDA0IQGPJ2d7tUBAAAAJNi8efMkSbMkDZOkhgY3ywlDwAMAGJrQgMeY4KpaAAAAQAZ56aWX9OKLL8ojaaZz8sgjXawoHAEPAGBoQgOe0aNptAwAAICM5PTeuVbSWOfkHnu4VU4vBDwAgKFxAp5DD5WKi92tBUhjlmVNsyyrvq2tze1SAABAhLfeeksrV67UyJEjdb1z8k9/crOkXgh4AABD4wQ8xx3nbh1AmjPGrDTGVBUUFLhdCgAAiFBXVydJuvzEEzXROXnyya7VEw0BDwBE88wz0pIlbleRHh591N7efbe7dQAAAAAJ8K9//Uu//e1vlZubq9lPP+12OTER8ABANK2t0r//7XYV6eH11+0tvXcApICWlhZ5vV63ywAAZJD58+erq6tL06dP1z5uF9MHAh4AiMayWA1qoJyvEwEPgCRpaGhQc3Ozampq5PP5wh7zer0qKyuTz+dTQ0ODGhoaVF1d7VKlAIB0t2nTJj3wwAOSpJqaGper6RsBDwBEQ1gxcJdcIk2cyNcMQFI0NzerpaVFJSUl8nq9am1tjXqdM4qnoqJCHo9H9fX1ySwTAJAh7r33Xm3dulXf/va39dVNm4IP3HWXe0XFkOd2AQCQkhjBM3B5edKYMcHjK6+U+Gs5gDhqaWmRJBUVFamkpERFRUWqr69XdXW1ioqKeq5raGjQeeedJ8kOdkKfX15entyiAQBpr729Xffcc48kae6cOdIJJwQfnDXLpapiYwQPAMRCwDMwnZ3SD37gdhUAMlhtba2am5t7jj0ej6qqqrR27dqe8EeygxyPxxP2XK/Xq/LycpWVlSWtXgBAZnjwwQfV2tqqo446St8MDXdSFAEPAETDCJ7+GSN9/rnU1SVdeKHb1QBZy7IsV/4lk9NXR7KXqW1oaJAkjR8/vmffmbYVqrm5WYWFhaqqqgoLiAAA6E9nZ6fuvPNOSVLNDTco7Dffjh2u1NQfAh4AiIZ+Mv1bu1Y67TQ74MnPd7saABnIaaQsScuXL1dLS0tPT52GhgatXbtWVVVVksJDIOe5lZWVqqmp0dSpU2P26gEAIJqlS5dq3bp1OvDAAzXt+98PPrBwoTRsmHuF9YEePAAQDSN4+vf730svvih95ztSbq7b1QBZy6TB/1fV19ertrZWHo9H1dXVPaGMo6GhIaxnjqOkpESNjY2qqKgIe47Tdyfac0Kfu3bt2jj9FwAAsonf71ddXZ0kac6cOcqZMSP44LXXulRV/xjBAwDRWJb0m9+4XUVq6+y0t8Yw4glATC0tLfL5fFq7dq2ampokSZWVlT3Lm7e0tIT10Ym0atWqfhske73enubKAAAM1dNPP6233npLe06apAtDw52TT3avqAFgBA8ARGNZ0jvvuF1FauvqcrsCAGmgpaVFc+bM6TmuqqpSWVmZKisr1draqqKiIj3wwAMxn9/c3Nxvg2QaKAMA4qm2tlaSNPvrX9ewxx8PPvCnP7lU0cAQ8AAAds7o0VLEajU9/vMfabfdklsPgJQULXwpKirSqlWr+n2uz+frWRUrsscOAACJ8OKLL+qll17SuHHjdEVouPOrX7lX1AAxRQsAsHMmT5ZmzYr+mDN9CwCGwOPxqKysTF6vt6fvDgAAiTRv3jxJ0nWbN2uMc/LCC6VLLnGtpoFiBA8ARHP22W5XkPo6O6Vbb5Xuv9/tSgBkMGeYPAAAibZmzRo9/fTTGjlihK7bvj34QKw/aqYYRvAAAHZOXz14aLoMAACANOOsnHXF9u2aEPrA1Kmu1DNYBDwAEMupp7pdQWrz++1t5BLNo0YFHwMwYJZlTbMsq76trc3tUgAAyDofffSRli5dqtycHIWN10mjhUUIeAAgliOPdLuC1Ob8stu2Lfz8736XVr8IgVRhjFlpjKkqKChwuxQAALLOXXfdpe7ubl3g92uyc3LzZik318WqBoeABwBicVaB+vRTd+tIVU6IExnm5OVJTzyR/HoAAACAnbDxjDO0eOFCSdKc0AdirRibogh4AKA/06a5XUFq6u62t3/9a/j53Fzp5puTXw8AAAAwWJalhStXapuk0yQd6pxPw1VhWUULADA0Tz0VfpyXl1ZDWQEAAJCluru1RdLCwOFc57zfn5aLhjCCBwAirVkTfhzZRBh9y80Nju4BAAAAUtWiRVosabOkb0g6RkrbcEci4AGA3trbw4+bm92pI13l5RHwAAAAIOV1XHON7grsz33ySfsPu2ka7kgEPABg+/zz4DzbNJxvm1KYngUAAIA0sFTSx5IOlnTaaae5XM3QEfAAgCQVFEgXXGDvE/AMzRtvuF0BAAAA0Cf/K6+oNrA/5+67lZOT/vFI+v8XAEA8lJZK775r7xPwDM2JJ7pdAQAAANCnlUcdpbclfVnSBVdf7XY5cUHAAwCSdMQRwX36xwzNyJFuVwAAAADEZJYu1bzA/mxJ+fn5bpYTNwQ8ACBJkycHV8/q6gp/bMqUpJeTVjye8GN68AAAACBVnXOOXrjwQr0iqVDSFR9+6HJB8UPAAwCSvfKTIzLgSeNO+q4I/VoCAAAAqeTxx3t671wnafTee7tZTVwR8ACAFBx18sEHBDqDseuu0rPPhp9zvpaR5wEAAACXrZb0B0mjJF37zjsuVxNfBDwAIAVDiaIiafRoqbDQ3XrSRVmZVFISfs75Wr79dvLrAQAAAGL57W9VF9idOWWKJuy/v6vlxBsBDwBI4Y2Vu7vpIzNQjzzS+5wzRcvvT24tAAAAQB8+uOgiPSopT9KsJ55wu5y4I+ABAEm64Ybgfne39NRT7tWS7pxwjIAHAAAAqWL7dt0lqVvShaWl2muvvdyuKO4IeAAgUleXtPfekjFuV5KecgK/Wv7rv9ytAwAAAJCktjZ9NnKkHgwczlm82NVyEoWABwAiMUUrPnw+tysAAAAAJI9HCyVtlzRN0sGHHeZyQYlBwAMAkQh4BmbRIrcrAAAAAPpmjL6QdF/gcO6ll7pZTUIR8ABAqG98IzzgYZpWbK+/HvuxMWOSVwcAAAAQy0sv6QFJmyUdc/jh+sbDD7tcUOIQ8ACA0wx49GjJssIDnu7u4KpQGLgRI9yuAEg7lmVNsyyrvq2tze1SAADIGB3f/KbmB/bn3nabq7UkGgEPAHR02Nsf/zgY8DihzuLF0urV7tWWrlhBCxg0Y8xKY0xVQUGB26UAAJAZOjr0iKRPJB0i6dRTT3W5oMQi4AEAJ4wwRnr5ZXsVrdxc6bPPpKuukrZvl955x90aU9HZZ8d+rLs7eXUAAAAAUfiHD1dtYL/mhz+UZVmu1pNoBDwAsHatvTXGDnucKVrOyJ7cXOmCC9yrL1Wdckrsx/LzpcLC5NUCAAAARHhS0j8l7SVp+q23ulxN4tFYAkD2eOUVO3TYb7/w8x99ZG+rqoJTsywruD9sGCNSBuvb35aOOcbtKgAAAJClTFeX5gX2b/rxj5Wfn+9qPcnACB4A2eO556L309lzT+nyy6VddpFGjQqed1bQMkZasyY5NWaKJUvcrgAAAABZ7K/5+fq7pAmSvjtrltvlJEVCAh7LsjyWZVVZllVmWVZVyPmKwLkKy7JKEnFvAIjp00+jL3s+bJg0Y0bw2PkF4PTmKStLfG0AAAAA4uPXv+4ZvXOdpFGjR7tZTdIkagTPCmNMvTHGK6lckizLKpJUbozxGmMaJN2SoHsDQHQLF4YfP/OM9L//G2yqLNkjfNrb7X0n4Dn//OTVCAAAAGDn+f36v0su0Z8kjZZ0TRaNxI97wGNZVpkkn3NsjKkM7IadD1zbaxRPYORPo2VZjRs2bIh3eQCyXegInpNPlj75JNhUWZLuuy/4+P/+r73N8G77AAAAQMb42tdUF9idue++Gn/IIa6Wk0yJGMFTJNlBT2AqVkXgvEfSppDrWp1rQwVG/pQaY0onTpyYgPIAZLXIKVrDh4cHPL//ffCxv/0t+nMAAAAApKSWNWu0TPaKUrOefdbtcpIqEatoeSQpMD1LlmWtsiyrOQH3AYDBiwxr8vLsgGfECPt4l12Cj02YIG3cKDGaEAAAAEh9HR26S5Jf0iXf+Ia+vNdebleUVIkYwdMS+OfwSSpRxPQsSYUR1wFA4kUbjRPagyd0+cQ33rC39fWJrwsAAADAkPznzjv1UGB/zh13uFqLGxIR8HgVPvWqSFKzpOWSikPOe4wxjOwBkFyRAc+554ZP0crPl6680t6fNMneOs2WAQAAAKSse+bP13ZJZ0o68BvfcLucpIv7FC1jjM+yrGUhy6PfboxpkSTLslYEmjBLUm287w0A/Yo2guef/5S+9CV7Pz9f6uwMf5yABwAAAEhpn7/wgn6+yW77WzNsmMvVuCMRPXgUWAY92nlvIu4HAAN28cXSRReFn/ve9+zl0SUCHgAAACDdbNum+mOPlU/SsZKO+slP3K7IFYmYogUA6WXrVrvZstQ74Nl11/C+PAAAAABSyo4FC7QgsF8jSbNnu1iNexIyggcA0k5oD57QgOff/5a++EJ6+GFXygIAAADQh82b9Zsf/EDrJR0q6dvvvy/lZOdYluz8rwaASLECHsuyl04fPtydugAAAABEt22bugsL5ayXNfeAA2QVF/f5lExGwAMgO0U2W44V8AAAAABIPevWSaNG6feS/ilpsqTz1qxxtyaXEfAAyE7d3eHHBDyD4/fbo5sAAAAAN0ydKiNpXuDwptpa5eVldxcaAh4A2eP884P7Trp/6aXS7rsHmyxbVvRVs6Itr57NuruDoRgAAACQbBs26HlJr0makJ+vGdde63JB7iPgAZA9QputrV9vbx9+WPrKV4JhRU5O9DDniy8krzfhJaYNAh4AAAC4JTCS3Bm9c8Ott2rUqFHu1ZMiCHgAZIcdO6QJE4LHX/5ycL+zM7gUumVFD3ja2qSFCxNbYzoh4AEAAIAbAu9BX5f0jKTRo0fr6quvdrWkVEHAAyA7PP+89Otf2/vG2IHPiBH2cUeHNGyYvR8r4MnN7d23J5sNJOB5/317u21b4usBAABAdgi0U6gNHFZXV6uwsNC9elIIAQ+A7GBZ9r+DD5a6uuxQ5+WXwx+X7Cla0Xrw+P2MWAnV3R3sWxTLHXdIb7whHX98UkoCAABAhnvySUnSWkkrcnKUn5+vG2+80d2aUggBD4Ds4Pfb4c2cOXbAs2OHNHp07+tijeDp7g7v4ZPtfvrTgQVe993HyCcAAAAMndcrnXmmJOlOSX6/XxdffLH23HNPd+tKIXxaAZAdnGW98/KCAc/w4b2vi7WKFiN4wj3yyMC+HsawAhkAAACGbsYMSdK/Jf1SkmVZuvnmm10tKdUQ8ADIDsbYI3Dy8+2AJ7SxcqhYq2gxgifcQAOvjg6CMQAAAAzdxx9Lku6RtEPSmWeeqQMPPNDVklINn1YAZAdnilZenh3uxOohE2sEjzHBPj2wQ7L+evBIdsBDMAYAAIChePZZSdLnkn4+ZowkqaamxsWCUhPvugFkh8gpWrFWgYrVgwfhuroGPoKHgAcAAABDUVYmSVokqW3LFh1//PE68sgj3a0pBfGuG0BmmznT3v7jH8ERPF1d0pYtwYAiNNDJy5OWLIn+WgQ/QW1t0saN/V83fDhTtAAAADBk2yXND+wzeic6Ah4Ama252d5+//vBHjydnXaTtmjBQ06OtN9+vc8zPau3LVv6fnzKFPsfAQ8wIJZlTbMsq76trc3tUgAASB333itJ+rXsBsuHHXaYTj75ZFdLSlUEPACyR26uPULnppvsY6eHzJ139v9cevD0Fm0VslDHHENzamAQjDErjTFVBQUFbpcCAEDquO46dUu6I3BYU1Mji/flUfGuG0Bmsyzp7rvt/bw8+98TT9jHzsiSE07o/3V27JCGDUtMjelq1Ki+H7/nHjvgee456bHHklMTAAAAMs7jkt6TtM8++6iystLtclLWAJZAAYA08fvfS2eeGX6uqUlqabH38/PDV34azNShjg5pxIih15gJnF5EA1lFy1mRbO3axNUDAACAzHTJJTKSagOHN910k/IG8h40SzGCB0DmOOus6Oc3b5Z2390OJPLzg+cHM3XoqKOkL31paPVlipdesrcDGRp766321ueL/vgLL0hvvBGfugAAAJBZfv1r/a+kRkkTJ07UjBkz3K4opRHwAMgOt91mj9gJHbUzmLm7l10m7bln3MtKS86onMEEZB0d0c/xhbK1AAAgAElEQVS/+qr05ptDrwkAAACZJTBqfF7g8Hvf+55GjhzpXj1pgIAHQHbw++1Ah6XOh84JxgYT8HR2JqYWAAAAZKbWVjVJ8koaM2aMrrrqKrcrSnkEPACyw7Zt9rary906MoET8AxmBFSsrzurkwEAACDSb34jHXBAT++dK6+8UuPGjXO1pHRAwAMgO8yebW8nT3a1jIwykJ5ERx9tb/tbUh0AAACQ7D/Mfuc7em/jRj0mKV/29Cz0j4AHQOYoKen7ccuSJk6U7rgjOfVkqssvt7f9fb2lYEPmQw5JXD0AAADIDN3d0qhRkqQ7JfklXVJaqj322MPVstIFAQ+AzDHQ/jqD6R2D3rZsGdz13/++/cva8fTT0oYNwWOmaAEAAECSGhslSZ9KeliSJenmhx92r540w6ccAJlv/Hi3K8gsRx01uFDm1FODK29J0umnS2+9FTym8TUAAAAk6bXXJEk/k9Qh6exzztH+Bx/saknphIAHQGa7+WZpr73s/Z/+NHj+yCMH/1qxlvrONuXlduO7gYocMVVcLO2yi72fmxse/gAAACB7XXed2iT9InBYU1PjZjVph4AHQGbLybFDBEn61rfsrWVJV189+Ne6/vr41ZXOLEu68MKBXRttOctjjpHefddeOj03l5XNAAAAIG3cKEm6X9Lnkk444QQdccQRrpaUbgh4AGQ2y6Lnjpssy/4X2oPnkUfssKy1VVq/PvwxAAAAZKeJE7Vd0oLA4dy5c92sJi3xqQdAZojVx2X4cKmoyO4Dg/jYmabI11wjPf+89OGH9oidzk77fG0tAQ8AAAAkSUsk/UfSlClTVF5e7nY5aYeAB0Bm8PuDU7FC7b67tHBh+DSgSy6RzjgjebVlmp1tivznP0t/+5u97wQ8krRjx9BrAgAAQPoyRt2S6gKHNTU1slhpddDy3C4AAOKiq0vKy5M+/ljac8/wx3Jzw0eJsKpW8vzP/0hjxtj7Tzxhj6aSggHP3XdLY8e6UxsAAABSQ06OHpO0VlJxcbHOPfdctytKS4zgAZAZnIDnK1/padDWIy+PaUBu8Xjsr78kvfNOcMWsjg57qtewYfb+0qXu1QgAAAD3rF8vI6k2cHjTTTcpL4+xKDuDgAdAZujosIOE7dulN98Mfywnh4AnVYROzTLG/t4YI91xh3s1AQAAwD3XXCOvpGZJuw0bpssuu8zlgtIXAQ+AzPCDH0iNjfb+8OHhj+XmSlOnDv0ee+wx9NfIdh0dwX1nNM9g3XWXHeQBAAAgva1eLT3xRM/onRt+9CONGDHC1ZLSGQEPgMzQ3h5s1pufb2+NsUOEESOkBQtiP3egDj1Ueumlob9ONpszJ7i/swHPk08yIgsAACATHHaYGiU9K2lsXp6uuuoqtytKawQ8ADLDkiV2PxcpuMrTH/8ozZ8fv3u8+qq0ZUv8Xi+bXHSRvQ0NZnY24AEAAED6C7xnd0bvXHXOOfJ4PO7VkwHoXAQgc7S321snONi2LXxK0FB1dERfih07x++3R0QNdhUteioBAACkv1tv1buSHpM0LDdX37v7brcrSnsEPAAyjxPwdHYGR/XEw7Jl0re+xciTePH7pV//2l757L33Bv68nBy+BwAAAOnuJz/RHZKMpEsvv1y777672xWlPaZoAcg8zuiOjo5gP5542Gef4PQvDI7f3/t74YQ0gwl3JAIeAACAdLd+vdZLWiLJkr00OoaOgAdA5igutrfOh/8dO3qvqIWhMWbnpkdZlt3sOpTzffrxj6M/57rrop9/4w0CHgAAgHS2xx66W1KHpHMrKrTffvu5XVFGIOABkDmOPloqLQ1++P/tb+PbFNmy4vda6eq556Sf/GTwz3vwQWnffcPPOd+n3XaL/pwXXoh+fsMGevAAAACkqxdekE/S/YHDmpoaN6vJKAQ8ADJHTuD/0pzg4PnnpS++iN/r5+fTZLmzU9q+ffDPGzEiGJA1NNjb7m7pzjujX+/32yN1YnG+xwQ9AAAA6cPvl449Vr+Q9IWkk447TqWlpW5XlTEIeACkv/Xr7e2ECdIppwQ//J98srR8efzuU1wsnX9+/F4vXQ31l/BBB9lbv18aM0baurX3NV1dfb+G8z0+8URp5syh1QMAAIDEM0bKzdU2Sc56WXN/+EM3K8o4BDwA0t/mzfb26qulyy8Pfvi/6KKhhxGRHnkkvq+XbkaPlubP3/nnT50qHXigve98nz7+uPd1/QU8zz4rNTbaI7Ref33n6wEAAEByBBop/0rSZ5JKvvY1nXTSSa6WlGkIeABkjn32sadpOdN2urqkvDx3a8o0Q1l6/vPPw4+dgCeyYfLbb/c/9eq11+zr3n9f2rRp5+oBAABA8vzqV+qSdEfgcO73vy+LHpdxRcADIP2NGmU3VJbCl9Du7qZnTrwNZen599+Xli4NHkcLeFpbpXPOkZYt6/u17r3Xft4XX0gffrhz9QAAACB5zj5bj0lqkbTvvvvqnHPOcbuijMOftgGkP79fmjjR3s/NDQYGHR3Bxsvxctpp8X29dNPZufMBz5Qp0le+EjyODHg+/FB6+GHpnXek22+P/hqhI3tirbIFAACA1LJ9u8zixZoXOLz55puVyx9i444RPADSX+hIndARPNdcE/+lzadNk6ZPj+9rppOhBDyR34uODnvrfL/+9a/gY7F+4Xd2BvcffHDn6gAAAEByvf++Vkn6P0m7jR+vSy65xO2KMhIBD4D0FxnwrFuXuHstXZrdU4KG0oMncpWE00+3t8bYW2f59L50dUnDh+/c/QEAAOCOc87pGb1z4803a8SIEa6Wk6kIeACkv8iA57rrEnevv/zFXto7Ww2lB8/ZZ4cfO6NxnBE8v/pVMOxxtpFef13af/+duz8AAABc8ep77+k5SbtIuvLKK90uJ2MR8ABIf5EBTyLtt9/Oj2DJBEOZohWLE+aErrL1/vvRr33xRWn16t7nI1fiAgAAQGp48knVBnavuvRSFRQUuFpOJiPgAZD+QpdDT/RSi+++G/+AI53EM+DZay9729UVnK7Vn1gjiLZti09NAAAAiJ/OTv3zzDP1uKThkm746U/driijEfAASH+hI3j6m+ITD9kc8BgTn1FSU6ZIN9xg73d3S4ceGnx9x5gxvb+Pe+8dfSUzp2EzAAAAUsc//6k7JBlJlx53nHafNMntijIaAQ+A9Bct4GltTdz9fve7xL12thg/XnrySemqq8KXPv/v/w7uV1aGPyZJkybZK5lFirwOAAAArvvk+ee1RHbwcPPixW6Xk/ESEvBYllVrWVaRZVkey7KqQs5XWJZVFtiWJOLeALJQaMAzdqy9fe+9xNxrzz0T87rZ5s9/lnw+O5Brb7cDn0j77987uOnulkaN6n0tAQ8AAEDqaGuTLEsLrrtOnZIqiou17777ul1VxkvUCJ4SSaskPSBpuSRZllUkqdwY4zXGNEi6JUH3BpBtQgMepwFyoqbszJghHXNMYl47m+TkSC+9JE2dKh12mHTFFb2v2X333s2Tu7uDIZ7jrLMIeAAAAFKJx6PNkhYFDmuivddD3CUq4FlkjCk2xlQaY3yBc2WSfKEXRRvFY1lWlWVZjZZlNW7YsCFB5QHIKKEBj6OjQxo9Ov73ysuTdtkl/q+brYYPl269VYq2mkJubu/gpqtLKimRJkwInjv3XAIeAACAVPHnP0uSfi5pi6RySSVz57pZUdZIVMBTFDIVqyJwziNpU8g1rZKKIp9ojKk3xpQaY0onTpyYoPIAZJRoAc/27dJTT8X/XtFCh2yxcmX8X/P++4P7V10V/li0r3V3tx2yhf4BYPjw7P2eAAAApJrTTtM2ST8LHNZ4vW5Wk1USEvAYY+pCp2JZluVJxH0AQFL0gKe2NjGrXV10UfZO0TrjjPi+3p57Sp6QXw8PPxz+eE5O9IAn8ns9YgQBDwAAQKoYOVK/lLRBUmlpqU488US3K8oacQ94AqN25oScckbq+CIuLZTUEu/7A8hCkR/6J02Sjj46MQHP5MlSNo8u/PnP4/daX/ta+Pdtn33CH8/Njd6DJzLgGTmy93VACgksOlEW8f4IAIDM0tEhXXihurZs0Z2BUzU1NbIsy9WyskkiRvC0SAodg1VojGmW3Wy5OOS8J3AeAIYm8kP/qFH2FK1EBDzZrKREWr06fq8X+X075JDwx6NN0fr5z8Ofk5dnr8T1+efxqwuIs0A/whZJUZaLAwAgQ4waJS1dqhWSPpD0lUmTdPbZZ7tdVVbJi/cLGmOaA6N4imSP3KkJnPdZlrXCsqyywKW18b43gCy1ZIl05JHB464uaccOAp5U19XVO6xxOKN7IgOeF18MPmfSJPu6xx+3v98PPpj4moEQgcUiSo0x9SHnKmSPWvZIauGPWQCArNHdLaPgB/05//Vfyo0ceY2EinvAI0mB3jvRztNdCUD8PfaYtGhR8Li1lYAnHUQGPM6qZyeeKB1+eOyG1s5zcnLsqVmXXSatX5/wcoFQgT9YVUt6LeRckaRyY0x14HiFpEp3KgQAIIl+9ztJ0p8lvSFp91131Xe+8x1XS8pGiVpFCwCSKzQoeOop6Ze/tOcBI3U991z4qJ2xY+3tBRfYYY8T4ERyvtfDhtmPT5woMfwXSRb4o9WqiNNliug5GBjlAwBA5jJGOvdcSdK8wKkbb7pJw4cPd6+mLEXAAyAzhAY8znStHTvcqSUT+f124BJvod+3MWPs7RVX2AFPrBE8Th0vvBC96TLgHo+kTSHHzkITkh3+lARG+QAAkBksq+e92SuS/iKpoKBA1dXVrpaVrQh4AGSG0PDB2d9tN3dqyUSdnfEPzHJywsOZsWPtRs6SVFpqP9bVFf6cn/3MfiMh2Y38/P7wUUBAijLG1Btjyo0xUVcQtSyryrKsRsuyGjds2JDs8gAAGLznnw87dHrvXH311dpll12SXg4IeABkitDlF3Ny7KBgzz0Tcy9jEvO6qcxpXB1P3d1S6NDdAw8M7h97rFRYKG3eHP6cESOC+84IHyfg2bYtvvUBg+eLOC6UvXpWvwIBUKkxpnTixInxrwwAgHg74YSe3bclPSFp+PDhuuGGG1wrKdsR8ABIf4sW2aM5HE7YExr6xJNlRe8Nk8k6OxPf0+i008KPhw2z7xuL06PHCXi++CJxtQEDs1xSccixh1W0AAAZqb097PCOwFT7GTNmaDdG0buGgAcABqu/4CETNTRIH36Y+PuMGxfcz8/vO1RyRvA407yyLXSDqwKraJVLKg/syxjjk7TCsqyywLnavl4DAIC0FTI96+NXXtFvduxQTk6ObrrpJvdqQmKWSQeAjOYEPKHTi3bsCD/ONJ9/npz7PPNMcD8/v3eQFjo9zrLsY2cEz5132v+AJAisouWNcR4AgMx2+uk9uwuWL1dnZ6emT5+u4uLiPp6ERGMEDwAM1rBhwZElmzdLPp903HH2aJJMHUUyfry9LUnwis+hzbL7GymVlydVVQUDnldeSWxtAAAACJue1Tp2rBYtWiRJqqmpcasiBBDwAMhMiWyEHBrwNDRIS5faTYgfekhasiRx93XTZZfZ2wcfTN49BzJFa8YMpmgBAAAk0w9/2LP786uvVnt7u771rW9pypQpLhYFiYAHAAYvNOAJ3X/7bWn7dvfqSobQZc0TLdoUrWicmrq7E1sPkGCWZU2zLKu+ra3N7VIAAIjt7rslSVsl/Szwx7+5c+e6WBAcBDwAMFihI0tycoKjhRYsyOyluidMkA45JHn3izZFK9rKaM65rq7E1wQkkDFmpTGmqqCgwO1SAACI7uOPe3YfOvBAbdy4UUcccYSOP/5492pCDwIeABisyODhrbeCIU8mjyI5/vjELT0fTbQpWn1NvduxI7H1AAAAZLPzzpO+/GVJUqekO7dskWT33rGS+R4RMRHwAEh/K1cm936h07IsS1q82G62LGV2yBBooJc0eXnS73438OszOVwDAABw04IF0ooVPYfLJX20bp32339/nXXWWe7VhTAEPADS31NP9T73+uuJu19owOPo7rZXmJo0KXH3dVsy++9Idnj29NODux4AAADxN2tWz66RVHvwwZKkm2++WTk5xAqpgu8EAAxWtIDHsqSzzsrsPjDJDngGi4AHAAAg/mprg/sTJuiPTz6pNW+9pUmTJuniiy92ry70QsADAIMVK+DJz5eqqtypKRny8ty5b2hodt997tQAAACQjTo6pNAVsjZs0Lw77pAkzZo1S8OHD3epMERDwAMAgzVsmL1aVltbcNSIE/BkMrcCnhNPtLdbt0pr1kS/5vTTpYqK5NUEAACQyYyx39+GBjj77KOXX35ZL7zwgjwej6oy+Q+baYqAB0Bqe/TRvh9/+WVp112TU4sjP1867TSpqCh47phj7ODH0dUl/eMfyW8AnUhuTNG680472JHsr3GsHkcrV2Z2/yMAAIBk8np7n1u7VrWB6VrXXHONxo4dm+Si0B8CHgCpyxjpggv6vmbHDns0TTI5QY7PFxzBc8wxwRE8bW3SPfdI118v3XZbcmtLlPPPd6fHTegbh+5u6cc/jn1tX0uoA2nAsqxplmXVt7W1uV0KACDbfetb4cebN+sfb7+tJ598UiNGjND111/vTl3oEwEPgNT1i1/Y28WLY1/j90tffNH7/OmnJ6YmKRjw5OSEhwo5OdLee9uB02OPSRs3JnY1r2Q64QR37ltTI338sb2/erX9NY6FJstIc8aYlcaYqoKCArdLAQAgqKVF8nhUV1cnSbr88su1a7JH0GNACHgApK516+ytE/RE09kZ/Xwip0aFBjzd3eGP3XKLHToZw4iSePD5pP/8J3jMMpwAAADJs2SJtM8+WrdunR555BHl5ORo9uzZbleFGHinDCB1zZtnb/taetyNZcmdgOeqq6RNm4LnnUBnyxZ7P7QnD+KDgAcAACBx/P7wUdGlpZKk+fPnq6urS9OnT1dRaB9KpBTeKQNIfZ2d0adhSdKYMdLMmcmtx2k2vOuu0qxZwWbLzi/DG2+UXnmFMCKe/H57y9cUAAAgcc47L/x4//21adMm1dfXS5JqampcKAoDxTtlAKnv7belm26K/pjHI11xRXLrcThhQ+Ty6BMn2ttXX01uPZnMGalFwAMAAJA4jz0W3P/kEyknR/fdd5+2bt2qU045RYcddph7taFfvFMGkB5yc4PLZYfq7nbvQ39enr2NDHhOOSX5tcRDQ4PbFcTm9Doi4AEAAEiMNWuC+8uXS5Mmqb29Xffcc48kae7cuS4VhoHinTKA9LB+vfT8873P+/3BKVPJ5kzJcoIeR2Tj5XTh9DxKJQcdZG87OuxtfwEPja0BAAB2zoMPBvcrKyVJDz30kDZt2qSvf/3rOvbYY10qDANFwAMgPfj90YMTN0fwON58M/zYjcbP8ZCK4YgTog004DnzzMTWAwAAkIm6uqSf/czeD7RG6Ozs1J133inJHr1jhTZfRkoi4AGQmiLDhvb26AGP3+9OwHPzzcHGv6GBzrPPSn//e/LryVTOSmQ7dtjb/kZrffxxYusBAADINMuWhbccKC+XJD366KP617/+pQMOOEBnnHGGS8VhMAh4AKSm7m7pkkuk446zjzdvjj2Cx40pWqNHB4OdkpLg+eZmuym0o6IiuXUNVleX9MYbblcRW2TA01+Yt359YusBEsiyrGmWZdW3tbW5XQoAIJucf374cXm5/H6/amtrJUlz5sxRjtsj5jEgfJcApKauLumkk6Tvfc8+fv311BrBM2qUvXy7ZE8jMkaqrpba2uxeQd/5jnT22e5PH+vLJ5/YjasbG92uJLbvftfeDnSK1n/+ExxZBaQZY8xKY0xVQUGB26UAALLFhx/2PrYs/eEPf9Bbb72lPfbYQxdddJEblWEnpPAnDwBZravLbl581lnBRrvRetu41WR57Fjp88/tfWOkK6+099vb7W1enh38pHLAM22aHVI5QVUqcnrqDDTgkdK3BxIAAEAydXVJ++wTPDZG2ntvSeoZvTNr1iwNc0ZUI+Xl9X8JALjACXik4If6VGqyXFVlB0sHHSTdeqsdREnStm32Ni/PrivVm9GFBjyp2GTZCe8GOkVLSt9VzAAAAJLpxhujnn7ppZf04osvaty4cZo5c2aSi8JQpPCflgFktdCAx/mQn0pTtJx73nZb9JWbcnNTfwSPZI+MSeURPM73/kc/srd9fT2vvtreNjQktCQAAICMcO+9wf2QKe7O6J1rrrlGY8eOTXZVGIIU/+QBIGtFG8Hz5JO9r3OrybIj1r3z8uz+Np2dqTuipLNT2rgxvJdQLG6N7nG+vn/8o3TEEX0HPM4bk3ffTXxdAAAA6W6PPeztgw/2vA988803tXLlSo0cOVLXX3+9i8VhZxDwAEhNoQGP4/HHe1/n1gieaHbbLbiflyd9+qkdntx9t3s19eXNN+3pZU7Ak5MjPfVU9Gvdmmo2YkR4DQP5Xju9kQAAABDdf/4j/fvf9h/TKit7TtfV1UmSvvvd72rixIluVYedlCKfigAgQmjA01fnfrdH8IQGH/PnS5Mn2/u5ufYS3xMnSn/5iyulDUh3tx3wrFkjrV4tzZkT/Tq3RvDk50sLF9qrlv397/bKX/0h4AEAAOjNmOC/L33Jfh949NH24iGSPvroIy1dulS5ubmaPXu2y8ViZxDwAEhNoQHP7NlSSUn061JpBM+wYVJhoV2rM0Vr/PjU7nHj99v1NTba/Xjeflt64one17ndLNpZGWsgAc8XXyS2FgAAgHRjjP2eOSdHGj48eD7kvdX8+fPV1dWl888/X5OdP1oiraTIpyIAiBA5RSvWCBK3R/CEhkv772+POJHsmu66S9plF+lPf3KntoF45hk74AkNcN57r/d1bq6wlZMT7GMUOmUrkrNUPQEPAABAuGefDe6H/vFxwQJJ0saNG7V48WJJ0pxYI7qR8gh4AKSmaD14onFzBM+hh0pnnGHvr11rH//lL9I779i1H3+8Paon1UU2gd6yJfzYGHdH8MyYIZ1/vr3fV8Bz8cX2tr098TUBAACkkx/+sPe5V1+Vpk2TJN17773aunWrTj31VH31q19NcnGIFwIeAKkpMuC5667o13V3uxfwrFkTvHdRkb0dPtyempWXZ4/iGTnSndoGIzJI6+gIP+7sdDeoGjnS/jqPHRs+pDiSM3pq+/bk1AUAAJAOurrsXoaRDj9cktTe3q6FCxdKkmpqapJZGeKMgAdAaooMeE44Ifp1fr+7U7SimTLFDkQsK/Vqc4ROucrNDR+hE9kzqKPD/ZFI3d12s+2+vp6jR0sHHEDAg7RlWdY0y7Lq29ra3C4FAJBJAuGNJOmll6QVK6TW1p5TixcvVmtrq4466ih985vfdKFAxAsBD4DUlA5TtGKxrOBIk9Gj3a0llpUrg/vz5oU/FjmCJ1UCnv7CsoMPlurqpG3bklMTEGfGmJXGmKqCggK3SwEAZIrOTmnWrODxkUdKFRXSuHGBhzt1V2Ck/Ny5c2W5vbAGhiTFPhUBQMBAA57t21NzlMyXv2xvL7zQ3TpCVVXZ288/792IeNUqe/SL1DvgaW9PnYCnvzcdxkg7diSnJgAAgFS2aVP4e7h583r9YXTp0qVat26dDjroIJ1++ulJLhDxRsADIDUNNOCZOdO9ETz77x/7sXPPtbep9FeQBx6wtwUFvacxPfJIcH/r1vDHjj46NQKeffaRdt+97+v8/uTUAwAAkMo6O6UJE8LPRfTX8fv9qq2tlWSvnJWTaqPiMWh8BwGkptWrBxbwSO6N4Bk1Kvr5yFDn/vsTX8tgRY7SkaTPPrO3katQffbZwL8XifL3v0uVlVJZWd/XGZMeja0BAAASqbEx/Pijj3pd8vTTT+sf//iH9txzT11wwQVJKgyJRMADIPUYI91ww8BDhVT7a0NfS3mnCqdH0IknBs85zfa6usKv3bEjOOXMLYcdNrAg78wzpRtvTHw9AAAAqWzGjOD+mjXSXnuFPWyM0e233y5Jmj17toa5PVobcZFin4oAQNLDD9vb/gIeZzpOKk2DkqTnnnO7gticBsTDh9sNiVetCg7f9XrtbbQgpb+pUYn21lsDC3hyctwfbQQAAOAmv1/65z/t/aOPlg45pNclL774ov72t7+psLBQV1xxRZILRKIQ8ABIPU4D4P4+qHd3J76WvsSaLpSfn9w6BuOTT6SpU+0pWtdfbwciTz5pP1ZcbG9Dg5SmJns7Zkxy64z0r38NPLh57DF7+c9nnklsTQAAAKnE75eOPTb8vdyf/hT1Uqf3zrXXXqsxbr/PQ9wQ8ABIPQ89ZG8jP9Bfd134sdsBT12du/ffGc6S552dvYOoyZMln0864YTgOWflrVSYBjfQXkvPPCO9+64dZgEAAGSLN96QXngh/FyU8GbNmjV6+umnNXLkSF0X+f4aaS0F3rEDQIT33rO3kQFP5PBStwOedLRli90jyJhgaHPggcH9goLwIMUYady45NcZzWCaafOzAQAAss28eeHH778f9bK6wB8pr7jiCk2IXGkLaY1GBQBSj9NbJ9WnaKWbSy+1Ax6nwbLD45GmTIn+nP/7v9RZBWwwU98iG0UDAABkuuXLg/u/+U1w+n2IDz/8UEuXLlVubq5mz56dxOKQDAQ8AFJPZ6e9JeCJr/x86Ve/spc9d0I0x0EHRX+OMYmva6AGurrDN78p/eUv0t57J7YeAACAVPHii8H9V16Rvv71qJfNnz9f3d3duvjii7U375UyDlO0AKQeJ7iJDHi2bw8u5R16HQYmL0/64x+l5mbp1VfDH1uyJLjvhDpPPZW82uJp7lzpr38NrhgGAACQyR591P4Dl+Pww6NetmHDBi1evFiSNGfOnGRUhiQj4AGQuiIb+95wg3TaacFjAp7Byc+XNmyw9yMDnlDOsvNtbYmvaaBKSgZ+rdOrZ/v2xNQCJIhlWdMsy6pvS6X/7QEAUpsx0gUXhJ+LsTjGwoULtW3bNnArkDQAACAASURBVJ1++uk69NBDk1Acko2AB0Bq2bAhGDBE8/rr0po19j4Bz+A4PWyOPnpgPWpiNOZzRV8/E5GckV/77ZeYWoAEMcasNMZUFRQUuF0KACBdHHVUcH+XXWK+P96yZYvuvfdeSVJNTU0yKoMLCHgApIbPPrNXz9p1V2n0aGmvvaJft2OHdNll9n53d/gvtVSVKn1s8vOligp7Fa0dO2Jf195uL6f+ox9JBx+ctPLixgl4oiwLCgAAkDFWrpT+/vfg8WefxRy988ADD2jz5s06+uijdcwxxySpQCQbAQ+A1HDKKcERF1u2xPzlpPHjg4FJd7d0223JqS/dGWMHPJ2d9tfso49iXzt/vjRpkr1//PFJKa9fsQK/aPLypLKy8GaDAAAAmeaMM4L769f3Xik1oKOjQ/Pnz5fE6J1Ml9CAx7KsCsuyyiKPA9tBNFQAkHV22aX3ucMOk449Nnjc3R3st5LKBjO9KFH8fvuXfleXPUpq7NjY137yibRpk73/6KPJqa8/v/vdwK/Ny5N8PjvI+ve/E1cTAACAWzo6gvsHHSTtvnvMS3/729/q448/1sEHH6zTQvtZIuMkLOCxLMsjabokT+C4SFK5McZrjGmQdEui7g0gAzz/fO9zr70mnXxy8DhdAp6f/CS49LtburvtqVmdnfbX7NRTY18bGqJt2mRP2UoneXn2SmGSHVYBAABkihkz7D8eho7WWb065uV+v1+1tbWS7NE7ObFGySMjJPK7WyrptZDjMkm+0AsYxQOgR+Qol2HDel/jNAl2pEvA89Ofur9ktxPwPPOMPf2trwbVI0aEH8+aldja4i0vzx6xJMUcqgwAAJCWHn6497k+3g+vXLlS77zzjvbaay+df/75iasLKSEhAY9lWWXGGG/EaY+kTSHHrZKKojy3yrKsRsuyGjc4y/kCyHyRjYhj/aKyrOC17e3SqFGJrSsexoxJnYBHsr+2fQU848Ylp6ZkiBYUAgAApKPXX+997vHHY15ujNG8efMkSbNnz1Z+5B9LkXHiHvAEpmK17OzzjTH1xphSY0zpxIkT41gZgLTS1/BRZ7TP9u3pEfCMHClt3epuDZEBjzPCJZrIUS8TJiSurkQIDa/6+u8EAABIJ8cdF9z/xS/snoNnnRXz8hdeeEGvvPKKxo8fr+9+97tJKBBuS8QInhJJJZZlVUg6XNLhgdDHF3FdoYYQBAHIcAOZetXRkR4jNLZvl84+290auruDwU1/U7RuuCH8+OqrE1dXIoweHdzv678TAAAgXbS0SF98Ye/fdZd05ZVSQUGfT3FG71x33XUaHfr+CBkr7gGPMabB+Sc7wHnNGNMiabmk4pBLPcaY5njfH0Ca6uoKP441gseZnmVM+gQ8w4ZJn37qbg1dXXZvGqn/KVolgfZozjztU05JbG3xtv/+0pQp9j4BDwAASHft7VJxyEfp6up+n7J69Wr98Y9/1KhRo3TttdcmsDikkkSuolUiu7HydMuyiowxPkkrAsukl0mqTdS9AaShyB48/S0tfsop6RPwjBzpfg+ezs5gk+r+Ah5HQYFUXi4ddVRia0sE5+fnpZfcrQMAAGCoxowJ7hcUhI9WjsFZOWvmzJkaP358oipDiklYwGOMaTbGTDXGVAZG8CiwRHrPv0TdG0AaCg14+hluqvXr7dWg0iXgGT7cnqblps7O4NcqJ0f6/PO+r8/JsUf9PPNM4mtLBOfn6eqre4eHAAAA6aK1Nbg/bJjdd6cfH3zwgR599FHl5eVpVrqthoohSeQy6QAwcG+9Fdz/299iXxf6YT1dAp499pC+/nV3a+joCB/B88orfV+fl2eHQukqdARYR4d7dQAAAAxF6HvkAY4Iv+uuu+T3+3XhhRdqr732SlBhSEUEPABSywknSAceGPtxvz8Y6qRLwLP33tLFFyfvfps39z7nTNFasGBgDazz89M/4HF+jiL7OwEAAKSLf/zD3l56ad+rzAZ89tlnevDBByVJc+bMSWRlSEEEPABSxxFHSN5+Zm9ec01wdEZnZ7BxMIJOOqn3OScM+973BvTmQHl56R+M/PWv9pZGywAAIF29/LK9/epXB3T5woULtX37dk2bNk0HH3xwAgtDKiLgAZA6hg0bWPjgBDzGDOz6VNBf0+hEi2yy3J/nn7dHS6UrY4LhX3e3dMcd7tYDAAAwWL/4xf+zd9/hUZRdG8DvSQ8BQ5EqVapY6AgKKggiTUQQsYCCdAQBhYTX3mkiIEhRBBEUERRpIqKCCKJIEUFRIFKUIjW0tE3m++Nkvpnd7GZnk52t9++6cs3s7OzMQ0Kys2fOcw6wYIGsd+jgdveLFy9i+vTpAIDk5GQrR0YBKkg+GRFRyHCVFTJ5svvsHY2xVTrZW7bMeUFnY4AnIkKmjeUnKspUh4aAVbWqHsjKyQHGjAGuXPHrkIjcURSls6Ioc1JTU/09FCIi8rfPPpNmEZo6ddy+ZM6cOTh//jxatmyJW265xcLBUaBigIeIfKtVK+fbExKk25QZZ87IctAg74wplMye7TyQYaxXFBkJlC6d/3Hq1gVmzPD++Hxl2TI9g+eJJ2R5zz3+Gw+RCaqqrlRVdUCiu06CREQU2mw2oFs3/XGXLm5fkpGRgcmTJwMAkpKSrBoZBTgGeIjIt0xW/8/X5cuFP0aoUlXn08EcM3jcZT9FRABFi3p/fL6kBXg2bJDlN9/4bShEREREpvXqpa+PHQssXer2JYsWLcKxY8dwww03oIOJ6VwUmhjgIaLgMmGCv0cQnIxt0sNFVBRQqxaLLBMREVFw+eorff311902FcnJycGE3GvkpKQkKP6u/Uh+wwAPEflWYevmuJriRSIiwnlAIyvLvqV8MLSXLyxFkbbwWrHocuX8Ox4iIiIid44dA86dk3WtLIEbX3zxBf78809UqVIFDzzwgIWDo0DHAA8R+ZarOworV5p7faNGwGuveW88we6pp+wf5+Q475KVnGyfwaNNWwoHigJ07syi3ERERBT4VqyQZdu2QMmSbndXVRXjxo0DADz11FOIDreMbbLDAA8RBYbOnc3tpyhAXJy1Y7HCwIGm78J4ZMMGoHdvPXhhszkP8Pz1l32AJxwyeAD5/6LVJQrmtu9EREQUHgYPlqXJ5hAbN27Ezz//jKuvvhqPP/64hQOjYMAADxEFn2bN/D2CgsnM9P4xFQVISdELT9ts+jztf/6x3zdcgjrOmCksTURERORPp07p623bmnqJlr0zfPhwFClSxIpRURBhgIeIfMdbbc1vucU7xwl2WsCiSBHg5ptl3ZjB43jnJ1xTdrUMHgZ4iIiIKJAZ26HXru129127duGrr75CQkIChg4dauHAKFgwwENEvjN7thSOc1SQD95ma/aEstdeA7ZvB2JjgdOnZVt2tvMpWoDr+keh7uxZTtEiIiKiwLZ1K/Djj7LeurWpl4wfPx4AMGDAAJQ0Ua+HQh8DPETkW1lZkn56+TLQt69sy8522/4xjy+/9P7YrGY2kLV7N5CR4X6/fftkuWqV3jlr2zbJ6CGhBXWYwUNERESBKj0daN5cf7x6tduXpKSkYMmSJYiKisLIkSMtHBwFEwZ4iMi3Tp8GypQBNm2SQAbgujBwfsqWBZo29f74rPTQQ8DOne73e/FF4Phx9/vt2KGvZ2Xp6y1aeDy0kLVxoywzMoKzODcRERGFtqwsID5ef/zJJ6auWSZNmoScnBw88sgjqFSpkoUDpGDCAA8R+YaWYaIxThcqSAZP165AnTqFH5cvbdwINGzofr/IyLzfL2eKFtXXL1zQ12vW9HxsoUqrQ3T+vGQ2GYsXEhEREfmbMVunTx+gRw+3Lzl58iTef/99AMCYMWOsGhkFIY8CPIqiVDUuiYhMS0+3f2ycMvPvv54HeA4eBBYs8M7YfOGzz8zvu2cP8O237vczTjmKjZXlnXc633f6dPPnDyVaNlPr1sCAAcCSJf4dDwUMXtMQEVFAmDlTX58zx9RLpk2bhoyMDHTp0gXXXXedRQOjYGQqwKMoSgNFUcYBGK8oSj8AbawdFhGFHMf6JxERMsUoPV3aQHoa4DGT4RJIPGmRfvmy82LUjn75RV+32WR5//3SrSw723462BNPmD9/qElMBH77DShRIm+gkcIOr2mIiCigaNnFH31k6nr4woULmDFjBgAgOTnZypFREHL6P0hRlKqqqh7SHququhPATkVRquVuKu6DsRFRKHE1RevECeDIEc+nzjRoEFxtvz0J8PTsCbRs6dnxHb+/P//s2etD2bPPAqNHA3ffzQBPGOI1DRERBaypU/Ubcvfea+olc+bMQWpqKm677TY0a9bMwsFRMHKVwZPkbKOqqn8DOJN7cUREZF52tn3dHceW3SdPena8atWAm24q/Lh8xUxXLE1UlJ6R48zata6fGzRIlmfO2G+vXNn8+UONFghUVQZ4whOvaYiIKPBs2waMGKE/NhZadiEjIwOTJ08GwOwdcs5VgEfJTVt2ZrxVgyEiLwukLA6bTVJP33xTHp89K0ut+1OxYp4dT1GAunW9Nz6raQGe1q3d76so+Qd4nn3W/TFOngQee0x//Mwz7l8TqrR0Z0VhgCc88ZqGiIgCy7FjwLhx+uOtW029bOHChTh+/Dhuuukm3H333RYNjoKZ0wCPqqqDAHyqKMrTiqJcpShKa0VRliiKchaA+7LeRBQYBg709wh02dlAQgJQrpw8zsmRpTZ1qUgRz48ZTEWW09Jk2amT+31ffz3/AI9jPSNn+vWz76wVziIj9XUGeMIOr2mIiCigZGcD11yjN+DYuBG4+WYTL8vGhAkTAABJSUlQHLPhieAiwKMoSmsAjQA0BXAecodrlqqqJQE09t3wiChk2GzyQTsuTh5rQQotg8f4ITwU9ctNIND+/e5cuuT6uYwM4PRp+22xscA//8i69oZ/443682aCQqFKVYGOHWWdAZ6ww2saIiIKKH/9Zf/4lltMvWz58uX466+/ULVqVfQw0UqdwpOrKVrrIXPWF6uqGgFggPZE7px1IgoGgRTZz86WIE7XrvJYm7KUmSm1dAJprFYoXhxo2BAYMsTc/g8/7Pq5vXvtizbXqiVzuMeNk+9jRO6f9i+/1ANo4axcOb2LWEqKf8dC/sBrGiIiChzGa8G33zbVOUtVVYwfL7OKn376aUR52n2WwoarAM/9qqq2U1X1M+D/O05sVxSlm6Io9/lueERUKIEUNMnOljcwRZF6O48+KncssrL04E+o06Zp5Se/qVlGWoBnwgSgSRMpOj1jhmSrdOsmz6WmAleuFGysoaRrV+mgparAt9/q0wMpXPCahoiIAoOqAhs2yPpLL+k3oNz47rvvsG3bNpQuXRp9+vSxbnwU9FzV4FnmZFtq7vb/WT4qIvKOiIi87bP9RZuiBehjatVKpszs3Wvq7kXQ++MP9/v87SahYMYMWebexUGvXsC0afb7jB0ryzfe8Kw9e6jTgmeB8jtBPhFM1zSKonRWFGVOamqqv4dCRERWOHhQlsWKAc89Z/plWvbO8OHDUaQgdSspbLjK4MmP03ajRBRgnnlG2kMHyhQdY5aO9kG7ZUtgzhxZZycAobX0BiTbxNGBA7KcNUuW5coBJUvqzxcvrtfbiYnRp8KFcw0ejfa7YDZLisJBQF3TqKq6UlXVAYmJif4eChERWWH6dFmWKmU6037Hjh1Yt24dihYtiqFDh1o4OAoFHgd4VFX9xoqBEJGXrV0rwYJAyeDQpmgBeo2Y+Hjg+HGgeXOgdGn/jc1XGjVyv0+E4c+yswBPVJT+fRw82P65q68GPv4YqFxZHsfGys//4EFg8eKCjTmUaIEdBngoF69piIjIp+bPl6U2nd4ELXtnwIABKFGihAWDolBSkAweIgoWgZTB88svegZPr16yjIqSGjHFivlvXP7w1Veun9Pqw3TsCOzenff5yEg9wDNoUN7nVVXuCg0ZomfwpKYGTqDPn7TvLQM8RERE5GsXLsg1GQCMHm3qJQcOHMDSpUsRHR2NkSNHWjg4ChUM8BCFkvR0vb6IqgZWgKdvXz0wMWKELEuVkgBPONTfAfRpUjNnut5Ha3++ejWwcqXzfbSU3thYfdsXXwBnzuiPZ8yQAE+/fpIpNW5cwccdKhjgISIiIn955x1ZtmwJlC1r6iWTJk1CTk4OevXqhYoVK1o4OAoVDPAQhZIJE4A1a2RdUeQDfqAEeAA9g6d8eVlWqiSFh8MlwANIZs0XX7iuieNuKlVOjh7giYnRtzdoACQk2Ad9AGDLFuCtt8Lre+yK9j1ngIeIiIh8bflyWZqsO3nixAnMnz8fiqJgtMmMHyIGeIhCSUSE/YfXQMrgAfQAjzZ/WCsobCwsHMoqVACKFpV1V1Om6tSRFvKu5OTodXqMQZvmzaUmT6tW+jatwHJaWvh8j81ggIeIiIh86Z9/gJ9+kqxqLZPdjalTpyIjIwP33nsv6tSpY/EAKVQwwEMUyrZtC8wADwA0bKgHKMIlu2TFCj2LZNMm5/tce610QHOlenW9RboxWycyEpg40b4jQ0KCLBcuZIBH078/AzxERETkW1pjjIoVARNtzlNTU/FO7pSupKSAavhIAY4BHqJQowUQ0tKAI0cCq7iuMZAzb54ejAiXAI8x+PL77873sdkkGOOqdWZkpExxe+kloEwZfbuzKV+NG+vr4fI9dic6Wg/wPPWUf8dCREREoe/IEWDVKlnv3dvUS2bPno0LFy7gjjvuwM0332zh4CjU8IqfKJQYgwJafZZAyuCJj9fXb7pJXw+n4IMWiNGyaxzt3g3UrSsZTq7maKuqfTt1wP3P2Zg9Fc6iovQAz3ff+XcsREREFPqqVNHXn37a7e7p6el46623AADJyclWjYpCFDN4iEKV9oE+kDJ4HAsAa7SuAuGkWjXn20ePliyTzZuBypWd71OxogSBjNxNOwqkQJ8/GQM8rrKkiIiIiLxhyxZ9fdgwIC7O7Us+/PBDnDhxAvXr18ddd91l4eAoFDHAQxRqtA/yWoZHoHywb98euOoq58+ZmIscMsx0coqOlmCYq05bTZsC991nv83Vz1kryufqWOHEZpMLK9bgISIiIqv9+itw663646lT3b4kOzsbEyZMACC1dxTejCIPMcBDFEoiI4GePWVde0MIlABPly7OM3hatvT9WPztvffy/7loU9ZWrzZ/TFfZUdo0pHr1zB8rVA0fDnTuLN/7nBzpaEFERETkbZcuAfXr64937TKVOfz555/jwIEDqFatGrp3727hAClUMcBDFEqMnZICLYPHlQ0b/D0C36pdG6haNf+fi/ZzXLkSOHXK3HHvvVcKazuKjwe6dfN4mCEpMVFqH9lswJUrwH//yfasLCmA+Nln/h0fERERhYbRo/X1V14xdaNNVVWMGzcu9+WjERVONSrJaxjgIQolWmDgnXeA9HRZD/QAj2Ox4FDXvz/QoIH9NKHMTPspVMY39GnTzB1XUZzP605MBJYuLdhYQ5HWRctmk6luAPDFF8DOncC33/p3bERERBT8tmwBZs2S9XvuAZ591tTLvvnmG2zfvh1lypTBY489Zt34KKQxLEgUSrQAz7JlwN9/A//7X+AEeDiHWBcVZf9z6d8f6NNHD8oZAzyHDtm/Njvb8uGFtKgoCajZbEDx4rItM1MCjaxTRERERJ7atw84fx64+Wbg66+Bdu305+bNM32Y8ePHAwCefPJJxBs7zxJ5IMxunROFOC0bJiYGuHAB6No1cAI8/PCsi462/7lkZQFnzgCLFsljYzBs4UL71w4ZYv34QllUFNC2LfD553pANCtLAmf8P0pERESeOH0auO46oHlzuQ43BncaNwZKljR1mO3bt2P9+vUoVqwYhvBajwqBAR6iUDN9uv7BtXbtwAnwkE6bJmQUE5O3pf311/tuTOEiMlKWzzwjRayPHpWfxblz/F0hIiIizzjeiNNMmwZs22b6MFr2zsCBA1FcyzAmKgAGeIhCTVSUBAuAvJki/sQpWrrIyLwBHpsNWLLEftsdd9g/fvFFK0cVHrT/h1rx6tOn5Xdk4EDg+++Bs2f9NzYiIiIKHqmpwMiRzp8bOtT0Yfbv34+lS5ciJiYGI10dj8gkBniIQsmiRcDcuRLgmTlTAjwXLgB79/p7ZGSkKMCePfbbtIDPrbfq2xy7J2zcaO24woFjUW+bTQ+C/vUX8Nprvh8TERERBZf9+4EmTfTHX38N5OTIl6p61ERk0qRJUFUVvXr1QoUKFSwYLIUTBniIQsmmTZIOGh0tQZ7ISODPP4GJE/09MnL09tv2j7Ugg7GNpuOUrRIlrB1TOHDMJFMU+yy3jAzfjoeIiIiCy7x5QK1aEuQBgDZt5EtRPM5YP378OObPnw9FUTDa2FqdqIAY4CEKRdHRQGysrGdkOG+f7WssYOtaRARw8aKsaz83ABgxwn4/ZmIVnmNXiu7d9QDPXXfJnTciIiIiZ779Fujb137bunUFPtyUKVOQmZmJ++67D7Vr1y7k4IgY4CEKDZMm2T+OidEDBXPn6u23/eXPP1mDJz+bNgFjx8q6Vj8JkLtDs2bpj//6y7fjCkWOhQsPH9YDPNdcw0AkEREROZeeDtx5p/64QgXJnC/gNW5qaipm5V7nJSUleWOERAzwEIUEx5TOmBj7QIEH84At8fDD/ODsyDgVq2pV6eIEAOXL2+/neNHw7beWDiusPP+8LLX6R5GRzOAhIiKivNLS7LOABw4E/v1XWqEX0MyZM3HhwgW0bt0aTYz1fIgKgQEeolDTt69k7xjfhObNA264wX9j2r4d+PVX/50/EA0Zoq936SLLYsWAwYOd75+dLctWrawdVzh56SVgyhQpjAjI99ixsDURERGRdq2geeutQh0uPT0dU6ZMAcDsHfIuXskShYqcHMlIqFxZ6u4YU0gB/9dvOXDAv+cPVFlZUjMJkACDtu5o507fjSmcxMcDW7bIus0mWTxERERERsuW6esHDuSt6eehDz74ACdPnkSDBg3Qtm3bQg6OSMcMHqJQcNNN8uG0WjXJQMjK0qdlNWgQGEWWHe98kLjtNj2oU6tW3ilZkZHys83JYR0jb9m9G2jWzH7b/ffL95kZPERERGR04ACwYIGsv/suUL16oQ6XnZ2NibkdbpOSkqDw+o68iAEeolAQGSlBnago+dJqigASFPB3kWVyLSNDgnEPP+w8gBMbK+3SjZk+VDg33ABcfbX9NpuNGTxERESU10cf6eudOxf6cMuWLcPBgwdRvXp1dOvWrdDHIzJigIcoFERGShAnOjpvgEcrblyxon/GRu6pKrBwofMAT0yMBIEyMuwLZ1PBKQpQpIj9tuxs+b3RgmgsCk5ERETDhgEvvCDrn3wClC1bqMOpqopx48YBAEaPHo0oZg6TlzHAQxQKfvkFuHLFeYBH06KF78dFnnEWVNAyeCIjgU6dfD+mUGUM8FSqpGfwZGYC77xjP9eeiIiIws+xY8D06fpjx/qWBbB+/Xrs3LkTZcuWxaOPPlro4xE5YoCHKBg5a+U8bJgEd/79174gr5YV0rSpb8bmTK1aQPv2/jt/IMrMBFaulHUtsOMswKNl8MTFAeyy4D0JCfr64MF6gCcjA5g9G7hwwX9jo7ClKEpnRVHmpKam+nsoRER04436+tSpQKlShT6klr0zYsQIxAVCjUwKOQzwEAWjbt0A7QOAFhRYuVIyeDZutM8+2LFDCi0bP9D62qhRwJo1/jt/ILpwQb8rlJUlyx9+yLuflsEzdSqLLHuT8fdBq2GVliYBtawsCfSQ51SV09sKQVXVlaqqDkhMTPT3UIiIwpfW2OLsWXncuzcwfHihD7tt2zZ8++23KFasGAYNGlTo4xE5Y0mAR1GUNoqiNFQUpbuiKAMM27vnPtddUZSGVpybKCwsXw7s2SPr2dlAYqK8GUVHy2NHEYzlBpz0dAneAPpFg7MgXEyMBB4+/pjFsr3JOEVLUSQokZYG/PijrDPAUzCJifx7Q0REwSs7O2/DhRkzvHLo8ePHAwAGDx6M4sWLe+WYRI68fhWmKEpxAEmqqu5QVXUpgNm5268F0FZV1fW528d6+9xEYWXoUFlmZ+vTSVx1Wdq+3XvnTUvz3rHCWWoqcNVV7veLjdV/vqVLWzumcKIF03bvBvbvB155BThzRgI8+/YxmFZQFy/6ewREREQFt2WL/eP164GiRQt92D///BOfffYZYmJiMGLEiEIfj8gVrwd4VFU9r6pqWwDIzdKZk/tUGwDnjfs6y+JRFGWAoii/KIryy6lTp7w9PKLQ8euvwOrVQHKyPiUiO9v1NB5v1RQpUkSfUmTWWMZz87h0CTDzNy41Fci944Nrr7V2TOFEu3M2Ywawa5cUIT98WH+eGTxEREThp08ffT0ryyuFlQFg0qRJUFUVjz76KMqXL++VYxI5Y1kedW7w5gFVVQfmbioO4Ixhl7MA8nxaUVV1jqqqjVVVbVyad6uJnBs5UpZ79wJTpujbs7KARx6x31crEDd6tPfO36qV+X1VFTh3znvnDhU2G7BunRTEjo93vV/16kDJkr4bV7h4/HFZKgrQpo2sz5qlP89MNSIiovBx9CjwwAPAwYPy+KmnpHmJFxw7dgwLFiyAoigY7c3rcSInvPO/1glVVXcoipKiKMpBVVWrW3UeorCkVfHX5gSXLw8cPy4ZPN26AfXr6/tWqAD89pv3zt2woWfTV3Jy7D84kzC2ss8vmHDNNVIk+4MPrB9TONEy3dq0kSLgANC3rx74OXxYsni0OklEvpSWJtMHHdWoodeP+ucfvQCoo/h4oGZN/bGzY2muuUZ/TzlzRjoxunLTTfr6/v2u3wtKlAAqVpT1tDTgwAHXx6xRQw9ynzsnx4yJ0b+io1nXiYisdfEiULmy/bYuXbx2+ClTpiAzMxPdu3dHTePfZiILeD3Ak5u5UzK31s55RVGgKIo2PctYTaokgBRvn58oLGhtL5rLTAAAIABJREFUFY8ckWXLlsCSJXqBZeM0rS+/lIvjV1/17Rg1Nhs/JDtjDPDkJyqK04Ws9PXXzv9/Ll4sd+8aN3b92k6dgFWrrBtbMCpTBnjoIX+PIvj9/jtQr17e7b/8AjRqJOsvvwy8+67z1zdsaF97zdmxNLNnAwNy+2EsWwYMHOh6X2OHtJ49pUujM/37A3NyZ+j//nv+v0fGf9PYsTIeR5GRcoytW/Vt114rfx+jo+2DQTExErTt1En227xZfp+LFNG/4uP19R499ADSgQPyb4yLk2Nrx4+K0s9BRKFFVYF27ey3zZkjU7e94Ny5c5g5cyYAICkpySvHJMqPFRk8jSHTr4xSAPwCYLxhW3FVVV1cGRBRvozF3ooVkwvf1auBatXy7qsowNtve+/CNDNTsoXM+u033n11xhjgmTfP9X4REc47o5F3NGokv0Oam27Ssx3yy1S7fFl+55jlYy8yUg9AU8HFx0tmi7Ptmmuu0afgOnJ8rav9APspoCVL5r+v4zkyM50/d801+npcHHDDDa6PY/z/UqwYUK6cHDcrS5aZmfI3MCdH3y87G/j7b9fHfPhhff3XX4Hp053vpygyJUPTo4dMm3Vm8GDgnXdkfds24Lbb9OCPMRAUFQWsXAnUrSv7vvqq/K2IiLD/N2VmSrbt0qWyn6rKlNxixeTnUKKE/pWYCNx7r/6z+fFHYO1aWTcGuaKj5fvZv78+7tRUOabxfTg7W8YSG+u6bh9ROFBVoEkT+4D4Tz8BTZt67RQzZ87EpUuXcOedd6JxfsFuIi/xeoBHVdU5uYWSu0Nq7CSpqpoCAIqifJqbzQPYB3uIyBPGC7KLFyXTYORI13OFn3jC+V3RgtDas5vVr5+cn+wZAzzuumkxwOM72u9QyZLAlSuu99OCrDabfKAyZjaEs8hI1938yLy6dSWzJT8vvCBfZuQ3Rcuoe3f5MuOTT8ztd/315qcJT5woX46ys+3/ZkZEACkp9gET43qdOvq+LVsC06bJ77Pjl81m/35aqZK8p6any3M2mxxX+z3XZGXJPq6CwMa/2fv322ceGWlT4wA5b35Bq6pV9QDP1q2SweWMY4CnZUvJokpI0L8/WrDsiSfkBhAA/PmnFPSvUkW+KleWYNvVV8vfQy/VIiEKKP/7n31w5/x5Cah6SVpaGqZOnQoASE5O9tpxifJjyV9rVVXnuNi+3orzEYUdVQVeekm/uI+IsL87p9UUCQTR0Z533QoHxg8r11+f/74M8FgnMjLvtoYNpQbPjh3AXXfl/3qzH5zDxT//8IMgeV9kpP3vqqI4z1h15sYbzWclffGFuf2aNZMsPscgkLZeqZK+7/PPy7S37Gz72kKxsRJ00RQtKgGec+ekttK5c/pXaqp9FlTz5vr7v3ZOLcjl+DctI0PObeykqSjy3mxsZrJ/f/7ZpIcOSeAHAMaNA/btkwBVqVISVMrJkfPUri3ZRoAEwD77TG5inDkDnDghAaMOHaR2oNH583IcZgCSryxcKP+XNRkZXp+GOX/+fPz3339o1KgR7vRSNy4id3gVRhSsypVz/ZxjJy3Ae2nYTZpIerpZDPA4ZwzwlCiR/77MDrHOxo32j6dMAUaMkN+XsWOlzpU2LcOZ/O64hyOtLgpRKIuI0Itdu1O9unyZOWbVqvLlTrNm8mXGn39K8OfKFQkqxcToQSDjtLd69aQhwpEjEuA+fBg4dQo4fVoCTsZso3XrgO++c36+++/XAzzHjtlPlzOqU0cykbRsiUcekePWry8ZbAkJ8j559KiM7Y039DGvXi1/o6+9VgJKzgL1RPlJSwN69dIff/+914M7NpsNkyZNAiC1dxROhyQf4VUYUbC7+mpz+2VkSJq2VhegoB5/HGjb1vlz48cDjgXkYmJc12kIZ9nZQPv2UgTb3Zs+M3is43i3uGVLWWo/k5kz8w/wVK4sGT8kIiNdf6AjIv9wVSDamPlbqZLrAts5OfbvUy+8IB+OT5+Wr8xMOVZkpH2ntYgIqWuUmirZQmXKAH/9BaxfLxlAx47pAZ7y5eW9btu2vDeRjHWi9u8H7rlHf1ysmIx7xAj72k+A3Bzhh2pyRqthBcj/Nwtq4yxduhQpKSmoUaMG7rvvPq8fn8gVBniIglm1akDu3QG3/v1XUsW1go4FsW+fTFvROp44+vRT+wCPqspdEmbw5KWqUheiVy/3GTquOtWQNRQl/8Lgxp8Xg5f2IiPNB52JKDg4/j28/Xb5cqdqVee1mq5cAfbu1ad8AdIRbvJk4OefZTpYWpqct2JF4Lrr9P3OnZMOaVlZkp106JBcB02aJBnGP/+s79upkwSniheXjKuyZSUI9MMPcj2kHXfWLKllGBEhU+EeeECuc1q0YIAoVGm1xoYMsSS4o6oqxo+XcrOjR49GJLPMyIcY4CEKRooi89XbtgXM3hUYP77w2Ts//SQXRg0bAps2yR05YzFLjdZZ6NAhudgy3m0jcfPNQIUKcuHpznqWL/OpuXOBu+92/fyhQ/r6999LUCMnh93iAPlecIoWEeWnSBEJxjgqVgxwV6ekWTPpUqbZvh14/XXgq6/kZojNJn+DTp2SLA3jNDSjm2/WAzxLl+pF9bdu1Ytilysn03a1Tmu7d0sR8O7dJXjED+3B57ffJHNa+3/RqZMlp1m3bh127dqFcuXKoXfv3pacg8gVXoURBSNVlTTnBQs8646VkVG489psUlNHUYD335d26cY01wsXZFuXLhLY0VLCXU3pCmevvy7LuXPtaxs4U7++fJFv1Ktnn6WTnW1/IX/mjCzj46XYefPm8rtRkPn7Z8/K75O7OkzBQmsTTUTkC40aAcuWyd/plBQ946Z0aXm8Y4cEb9LSJDh/9Ki0wO7YUT/G6tV6oeqVK4EtW6Tg9okTEtDRAjxz5khh3oUL5QZNx46SST1woP00slOngF9/lWlvFStKtlLZsvYZS+R7aWkS2EtL07fldzOnEMblFm8eOXIk4lg4nHyMV2FEwSo+3j644s4HHwCvvFK4c2Zm6h9058/P+/z+/ZIpZLNJAOqPP2R706aFO28oq1DB/T4//1y4qXXkOWMR7MxM+X3TaFlX69ZJzZ6oKNl/6FCZZuCJTz6RoGm/foUfcyBo3ZqZTETke5GRQM2a9tu0lu/uxMbKFwA8+qh8TZ1q3zkMkOlbzz0nAZ6//9b/3r/+ugR5Jk2S94uaNaXukKPGjYHkZKBbN3l89KjUJdLOrapS06hUqfz/jtpscv64OPuObZS/bdvsgzuTJlkyBe+nn37Chg0bcNVVV2Ggq7pWRBbiVRhRMFIUeZMy3jFyp1EjvYBsQR0+LEVlXaU8A0CNGrLcvVvSYMk7mBXhW8bC1o6ZbzYb0LMncOut8jgrS6YHbN/u+XlcdaIJVrVrs2YFEQW/uDiZgm4M7g8aBLz8MnDwoExZnzhRrqsuXZLgNiCZnF26SKZIlSry97BKFWkV/8svMr0ekOuo7t0laDRligSVypaVgE+TJnogQlWBN98EVqyQrKJ+/aQrW61aQN+++thsNnbcdGfOHFlWrCgNFJ56ypLTaLV3hgwZgkStiDiRD/ETA1GwunLFfJtWALj+eim0XBjjx0sL1c8+04+ZlSUZCBqtfXtUFJCeXrjzkc74PSbrOWbwGGVlAU88IRfukyYBTz8ttbAK0k3r00+BNm0KN1YiIvIdRZHM5KZN5e//8ePSBUwzf74e6Nam+KalyVQwrW7iiRPAxYvAgQPAyJH2xz91Sg8s7dgh53DmyhWZGn/VVcCwYdIp9fhxqWWUkwP07y+ZRawVJMW5Fy2S9fHjgYcesuQ0+/btw/LlyxEbG4snn3zSknMQucMAD1EwUlV5YzfeWTLD2L60IDp3ltThESPkcUKC3M0yTv3SPhh37eq8iCIVTEKCv0cQXipXlotlwHkGjxZwM3Y4K+jdU2PKOBERBRdjcAewz2LUgivx8XrnJkCmZ//6q2SSrFoFNGgA9OkjU660m2iAvPePGCGNLeLi5IZAkybSxaxoUdnn8GGZNnbpkv04hg6VQNLLL5v/t6gq8M03Mh539QGDyfPP6+sWtiyfOHEiVFXFY489hnLaDU8iH2OAhygYKQpw7bX6m7tZkyZJSnFBdehgnzX088/AbbfZ7/Pgg3JhAMh8Z/IOFunzrZdekjugJ086D/BoU+aMHbWysjw/T6dOEkwKBevXS90K1hwgInIvOhoYPly+jHr10tfr1AHeeiv/41SpAuzaJbUW9+2TAMZff8n7mPHa4dAhmZ7kasr34cOSTfT55xJomjZNppsF+7Tb338Hpk+X9QkTLLue+ueff/Dhhx8iIiICT7vKuiLyAQZ4iILV5Mm+P6ezN3ljPZ6iReUO0oEDvhtTOLjppuC/wAo22p3WH35wHuDRimJmZMhdzjNngD17PD/P9dcDvXtLqn6wu3gx7x1kIiKyXvXqeTN1GjfW36sA4LHHpPlFjRpSw2fkSOD774HBgyXT6JlnJLgDSAHorl0lq2juXNn2009SM8jTm4v+pgV3AKl1ZJEpU6YgKysLPXr0QA2tHiWRH7DIMlG48UYRvnfeAXJbQNrVKpkwQZbGD6tffVX484U71t/xn5iYvDV4jBk86elAYYooakFRR0eOAP/9V/Dj+kN0tLSYJyIi/+vQAbjzTlm/ckVuQvz3n7SBnz9f/l4PGwYMGSLPa1PD6taVbGxF0W8wbdkCNG8O3HFH/o02HGVny1TmzZudP3/2rOvr0owMee6332QMiiJT0VJTzV/LLl8u0+AAmXpWpoz5sXvg3LlzmD17NgAgKSnJknMQmcUAD1G4cfywWhAlS+pvrleu6NudtfX05EKAnLv6an+PIHzFxubN4MnK0gM8e/fKnVPAXMt7R6VLO7/grFIFWLPG8+P50+LFwF13+XsURETkqEgRqePTtq3+OCEBuOUWCfAUKSLXc//+K7WBPvpIpn0NGiT7f/yxLLdvB779Nu/xbTbJ+nH01FPSxbVFCz1I07s38OGHsq61hF+0SAI3U6ZIC/hp02QqlWMZgF69gOLFpdGBO6oq59JonS8t8M477+DSpUto27YtGhak4QKRFzHAQxROOnXK+2G1ICIj9cDNe+/JUlXlQsAxm8EbAaVwx3bz/uMswGMssqztAxTs//rw4a4zdfzd8vbKFenCYtaRI/bt5YmIKHBcdx2wbp28t1y+LNmjmzfbZ15WqKDfwLjpJr1ezfTpwKuvyvqUKRKMAaRJwIoVUnuxcmU9iJOcLOcpUSLvOLZuzTud7JFHJHAzcqTUmNQ6UP3wg1xz1qljv3+5cnJNu3Gj1A5y5ocf9IzylBT76WpelJaWhqlTpwIAkpOT8+5gs0l3MyIfYYCHKJz06SMfVk+dKtxxOne2Lwo4fboUowXsM3oA52/u5BnW3/Gf2FjnbdK1C+CyZeXnM21awQI8+b3G3wGeCxekXoNZ6emWXUATEZGfPf64LFevlmCMVneuS5e8Nej275f3sEqVgE8+kelfrVvLjcDNm/VOVlrQCMhboP/RR6VodN268l50+LBkNHfsKFmuq1fLlLGqVfWyAUZTpsjynnuAatUK+693ad68eTh16hQaN26MVq1ayTj37gXGjJHrg+hoCZy1bg3MmyfbJk3y/3s8hSwGeIiCSWampMeuXVuw18fGyt2WDh0KN47YWKkdMnSoPN6zR1+vXVvfp2FDoGXLwp2LyJ9iYvJm8HTqpAd4WrQARo0C6tcvWBctZ4UYFy6Upb8v/g4ckACWWbVrS7o/ERGFnnLl7Dt8Xbok13916gA33yydu7Try/btZepV375Ajx5SsPmbb6QeTunSwBtvyI3BZ56RbKING4BZsyQTtEEDuWkyf75+TQlIhtCpU9JW3rgdAMaO1bN1Tp2SwNAXX8hjxy5lXmSz2TAxtzttcnIylEOHJOA0ZEjerrXffSffDwAYPVqymPz9Pk8hiQEeomBy7hzw7LN6lwNPxcZKWq23WkRqnQn++0/uznzwgaTqAvLmXtBxEgWK2Fjg4MG827UAz9KlcgcxOlrSsD2douSsleqmTbL094VfdLRedNOMtDQWWSYiCmWzZ8vXgQNAs2ZAkyYyjXzzZqBmTcmqUVWgX7/8jxMRodefK1IEuP12Wa9USYoyDxuW/+tvvlm6eo0eLY8rVpSaQa+9Jse99VZ5P+7aVS80bYElS5bg0KFDqFmzJu6tV0+mlwHSnczRZ59JkEzz4ovAW29JzaMSJSTIReQFDPAQBZPCfuDbskXmVHt7GoVWZLZBAwkedehgP42FKFjFxgIjRuiPL1yQpWNns4wM+T//4osFO49xfr72e+PvAI+xQ54ZLLBMRBTa4uOBAQP06z5AMlYiI30/lqZNpXvr4sXA6dNys7FjR3nu7FlZtm5t2elVVcX48eMBAGPGjEFkzZr6k2vXynt4drYsVVWCTceP27+3PvWUZACfPw+0asW6leQVDPAQBZPs7MK9iV6+LEtvB3i0ec6KIgGe9HTJevBWphCRv8TE2D8+dEiWRYrYb09JkaWxnoAntN9NQO9G5+8ADwsmExFRoHvgASncfP310hmsc2f9ub17LTvt2rVrsXv3bpSPjUUvY0OCHj2Adu1k3Vl32chIydYZNkzGa9Snj16oOinJ/9cBFJQY4CEKJu+/L28MBa1ro7WL9GbgZd06vXNWRIRkNmRmyvxnBngo2DkGeLSgh+P2W26RulSe0rrRGQtpR0RIW9ljxzw/njfZbHkDWURERIGmXj2pBVSzpn3HrVmzgOeeA3bvlulj7uzZIzdqhg2TrNxBg4CePZ3uOv6RRwAAIzMy8P+3TZs1k6LS7tx+u9QZ+vVXKTh9zz2y/aOP9H0mTJBOYKoq9Y28edMlLU0PJE2dap9V9Ouv7N4KyPdAUYAPP9S/V86+vNGd2MsY4CEKJu+8I39469cv2Ou1OwkFDbysWpV3W9u2+gfbiAj5QDh3rrx5saMOBTtFkQtEjasLrNq1gRde8Pz42u+ksUCzqsqdyPLlPT+eN2Vnu6+DQEREFCjS0vIWN758WQJAnTrJe/rYsfLelp4uz2/YIK/5/XfgxhslIDR9urRunz1bAjaKIvVy7rgDGDsWP65ejY1nzyIRwP/3/vrkE2DNGs/GW706sGyZFIRW1bzX93/9JdcJtWvL9G3HYMKxYzK2AQPka8wYvf6Pqrpu/mC8eTNihPzbU1Pleqd+fSm1MHiwnHfZMv142jWQMcBx3XX6sf74Q75HR4/qzwNSUNuTrpz+duWKXjBcuznuynPPWT8eD7FABlEwiY6W9otvvy2Rf09VrixLZymjZjgL8BiPFxEhX7VqAevX+2dONpGV8uuU5WnNGs2XX9q/9vJlufgyZvX4g83GOlpERBQcLl0CihWT9fLlpVTA999LHZ5vvpEsHkBaqvfuLfvfcguwc6dsHzPG/nhLltg/HjVKlhs3YnzupiEArnrrLema5ZjZWxA7dkhwISdHrqGN2UiAtKofNEgyhQDgmmtk+e67+j4TJ0rnMO1zwttvSyt77TOAM84aJGg3t7p3z3/M+/bJ9crBg9LSHrA/1/jxQHKy/vjmm2VKnS9lZ8sYzXz+uekm4LffzB974kT5mWgFtgMAM3iIgoljYVdPaXN9PZ3Tq3XLqlLF+fPGArEafjCkUOEYaHHVCrxhw4Idv0QJ+wCPzVb433VvsNkYpCUiosB35Yoe3AGABQukIHS7dvJ+unWr/Xt3aqpkp2jBHUeOU661bl0A/gDwBYDYmBg8eeyYZMB4I7gDyPVGQoL8W4oUyXvdvWiRlGmIjgY+/dT1cYw3gYcNcx3cWbIEePTRwo8bsC+8bWQM7gDS/czZDazISNnepo0EjTSffy7ZNFrhbGdUVbKEnBkzRj6TaMd/7jlZFitmn4n09NPyfXMW3DEGhurVs6/zBMi/XcsQ27LF9Th9hAEeomBSkBofjq66yvMAzwcfACtXSlqrM1qaK6v/Uygy/r4oCjBzpvP97rqrYNlxGRlyUfHhh3LXbsEC4O67CzZWb/K0qPuiRdaNhYiIyJW2bfX1KVMkSGAUHy+t3D/4QN6rhwyxDxi88ooEexYtkvfhixdlatLvv8s1wIQJEkRKTMSE3Gvxvo8/jrJWT6XetEnOn5kp61rGDgA89JD8uwBg6FDgvfdkelXp0u6P26GDHPf++4EZM+wDLhkZwIkTrgMmml27gOXLXT//8MP5v96xlo1Wk/Cbb/RpX+fPS42iL78EBg7UG10Y7dgh117lyjmvkeM4ZU9rhnHpkv32N98EnnzSfpvWAc3YDW3XLmDFCumUZqTVeOrXL/9/tw8wwEMUTAYNKvwx1qyxL+JmRny8FIBzl5VT0KlfRMEiMzP/O3XvvOP5MU+dkumPb70lqdSAfvHlT5s2mQ/wnDihz/snIiLyldRU+6wJxw/pRr17A199JUGeDRuAr7+WwMKzz0pw5KGH9GDHfffZ15eJj8fR337DoowMRERE4Omnn7bkn+NUdDTQooXUtjHSijBPny7Tt3bu1AMztWoB587p+370kfy7b78d6NZN356QIN8DLYAREwOULQuUKSPT2qpWleNoz2tf9erJ1C9j/cHBg/XnFy7Ut3/9tWy7csX8v1lRJMNZs3QpUK2adE07dUp+flOnSlMKb6tWDTh9Ov992rVzXrpizhzvj8dDnENBFEy0AMqMGQU/xq23yjInx3xARtvPmP5KFI7cBXgKolUruSh6+WXX6eL+MHmy1Cowo2NHa8dCRETkzPvv6+taRrkZt9/u8aneeustZGVloWfPnrjWHzVXtEwX7brcMQtF28d4g8i43q0b8O23eTOcXLnxRuDvv/Pf58UXgV69ZFp37dr2zzneqIqPl8DJ1VebO78zS5bkrY/kzsCBwHffSUe0l1/O+/zUqfaBwZQUc8ft2FH+jf36SYMZQAJxfsbb7UTBqH37wh/Dk3aL2h9oV60AGzQASpYs/JiIAp0VAZ6SJaU7niN/F1kGzNXSMnbWICIi8qXPP5flu+9a2r317NmzmJObnTHGsSCzLymKTCH74w/gtts8e21MjEwB93adzOrV8wZ3XClVSgJT99wj05pOnNCnc61dqxfDNurf3/1xa9TQp3k1aKAXov7jDykY/eefwEsv6RlGc+fKlLd9+6Qo9bJlwJ13Fqzt+Xvv6ccNAAzwEAUL7Y9GgwaSOlgYTz/t2QcyY1tEV6ZPNzfvlygYab9/7gI82dlA166eH/+LLwo2LqtFRbm/YLnvPucBKiIiIiv9/bc+ndg47cgCM2bMwOXLl9GuXTs0aNDA0nO5VbRo3g5bwSQhQa57OnSQ6WALF0pwpl07yRrKzpab2SVLAosXy7QnVQXS0vIeS5v+tX+/nr20Ywfw44+y7ur71Lcv8M8/emDqvvukA7C3b+L5AQM8RMHCm0VX69f3rKWz1hraVTpryZLAgw8ywEOhKTZWLyDuLsCTkZF/0UEzatSQZSDcCUpMlLte+dE6fiUlWT8eIiIizXvvybJzZ/t6LV525coVTMvtTJXE9zprGG8iR0RIzdAzZ6TmjiYuToIwgBSWVlXz083CCGvwEAWLLVskom3sFFBQf/wh81/btTO3/88/A//95zqAs25d4cdEFKhiY2Vef1SUBHDyC/Bo6cFmOcuKC5RAaVyc3GW7+WYpYulKixbSrrUwc+qJiIg80b693sno/vstPdX777+P06dPo2nTprjjjjssPRe5oU2j8vY0sxDC7wxRsLh0SVoHapHrwtizR5ZmAzxA/h862T2LQtlPP0mAMydHfv8efND1vp4GeIxZOps3SxH0H3+Ux/6uwTNlihRE1DL4XNGe58UWERH5QlqafZvq7t0tO1VWVhYmTZoEAEhOTobi7/dmColpVFbipzKiYGG2XbHZY7EoKpE558/L3O5z56SwYX6/i1rAxkxNmpwcPYhz//3ALbfI+u+/yzIqSp8a5g+KIgGetDTpMOGKNt2TgV4iIvKFe+7R19PSLP3Av2TJEhw+fBi1a9dGly5dLDsPkbfwaowoWHgzIFOrllSxNyMnxzvTwoiCVZs2QKVKwCuvAFeu5J9Zo2XwNGzo/rg2m16/Rmv5uXQpcN11sl6+PHD8uOvXO2uP6m3x8bL86Sfnz2/bJplHADN4iIjIejt36tnsM2bIdGKLqKqKcePGAZDOWRG8kUFBgP9LiYKJtz5ADR8uQR4zbLb8p6QQhbpOnYAKFWR9+/b899UCPLNmuT9uVpYe4NEYu4AULw5cuOD69Z62R/WElomkBXhcXUBv3QqsXClF4AcMsG48REREgH5DBAD69LH0VGvWrMGePXtQoUIFPKy18iYKcAzwEAUTb6WgxsSYn/px8SKnc1F4i4kBevc2t68nNXhOngR273b9fFSU+/o3VqlVC9i7V/87ceONQI8eeffTAlBffskMHiIist6mTbL84AP9JoRFxo8fDwAYNWoUYmNjLT0XkbcwwEMUTIx39wsjJsb8B8f337e/W0IUbjwJrGoBHjMtzjMzgebNXT8fHa3Xt3Fmzx7g5ZfNj80TBw4A06bpdXWqVZNOWY6efdaa84cYRVGKK4oyRlGU7oqimJi/R0REeSxapE8LNtbhscCWLVuwadMmFC9eHAOYoUpBhAEeomDirQ9T0dHAQw+Z2zc93bsFnomCjSf//x95BKhSxdy+2dlAmTKun88vg8dmk+fyq9HjDfXrAyNHOg80edoxLLwNADBHVdWlAB7w92CIiILOkSPyHgtIHcnixS09nZa9M3ToUBQrVszScxF5EwM8RMGiRAnvdanxJCMhK4vtCCm8edIS9dprgddeM7dvZmbeGjxGWgbP4cN5n8vIkKXVQZboaOCFF4CzZ/M+17y5TF1r1MjaMQQwRVEaKooywGFbd0VR2jhk6zRRVfV87vq1vh0lEVEIaNVKX8+vs6MX7N27FytWrEALdUfKAAAgAElEQVRcXByGDx9u6bmIvI0BHqJg8cYbQI0a3jlWRAQwc6a5fa++2tpirkTB5M473e9jthBjZmb+wVMtg6dqVX2bVg9LC/DMmWPuXIURFQUsXw6ULi2PtWye/fuBBQuA5GTrxxCAFEVpA2AsgOKGbdcCaKuq6vrcbJ2xhue0/ay97UxEFGouXwZSUmT9vffMv88W0MSJEwEAffv2RZn8Mm2JAhADPETBwEw9D0+ZzUqoVElPiSUKZ9OnA+vWee947rLjoqMlkNOkib6tdWtZagEeK2hZQbNnyzIqCkhMBPr3l6DUXXfJdm3qWvfu1o0lgKmquh7A1w6b2wA4b9yQm8WzDUDJ3E3nQURE5k2Zoq8/9pilpzpy5AgWLVqEyMhIPP3005aei8gKbHlBFAxstvyncoTquYkCxeOPAwkJ3psmCZjL4Fm1Cjh9Wt928aIsrQzw2GxA48Z6pk5UlNTiioqS82pZRGY78YWX4gDOGB6fhUzJmgNggKIo5wG84Y+BEREFpfPn9RqU/ftbXhdy8uTJsNlsePDBB1GtWjVLz0VkBQZ4iIKBuw+CVunZU7oUMMBD4a5hw/w7WhWEu9/r6GjpGGIM5miZd1YGeLKygF9+kfENGSJBrYwMGU96uh7kunTJujGEmNz6OxPy2ye3ls8AAKhcubIvhkVEFPhattTXp0+39FRnzpzBu+++CwBISkqy9FxEVuEULaJgkJXlnyDLgQP+OzdRIMmvo1VBuSuyHBUlGTsJCfq2s2eB1FQJuDRrlv/x16zJu81dUeaMDODoUSAuDujRQ7YpigS3YmOlDsKGDfqxYmPzP174cZx+VRJAipkXqqo6R1XVxqqqNi6t1TsiIgpXmZnA0KHAnj3yeO5cy292Tp8+HVeuXEH79u1Rr149S89FZBUGeIiCwcKFwA8/+P68qanyFcVkPwpzq1ZJNo1ZZoJBZmrwALLPmTPAO+8Ahw4Bw4dLIGbGDKBXL9evf/75vNsiI/UAjTM//ghcd51cVPfpo2+32aTg+sqV9vt/843rY4WnJQCqGx4XV1V1h78GQ0QUVGw2oG1bubEQGyvvewAwcCDQt6+lp758+TLefvttAMzeoeDGAA9RMNiyBbhwwbvHnDBBplvk58AB+QBn8XxnooDXpg1Qtqz5/Z94wv0+I0e6r8FTpQpw++3y+//dd7J9+XIJ8MTGAi1auH79b7/p9XKMevd2/Zq4OFl26QJUqKBvz8oCypSR4FKdOvr2665zfawQl9tFqy2Atrnr2lSsT3PbpLcBMN6fYyQiCiqzZgHr19tvi4oCJk+2/NRz587FmTNn0KxZM9zG7rEUxBjgIQoGt9/u/G58YaSkuC+S+txzwIoV5jtuEYWqZs0kwGHWiy+63+foUfcZPIcPSxctYye9++6T4Ky76VGZmUC/fs7P60y3bvp54uPtn7PZ9HbtLVrIfo88It21wlRuK/T7VVVtm9tRy7j9/7/8OUYioqAxYwYwbJj9thEjgHPngCJFLD11VlYW3nzzTQCSvaPwupeCGOddEAWDmBigYkXvHnPJEvcBnuLFvXtOomClKO7r1xiVL29uv/xq8GjFjKdOlXNv2yaPmzWTAsfFirk//vz5wLx55sZy6BCQlibrjhfTNhtQrpw+5kuXgPbtmd1HRESFN3OmnvmakCDlAXz4/rJ48WIcOXIEderUwT333OOz8xJZgRk8RIHuwAH5cOftN7r4eL0Tz/DhrvchIqBePeCVV7x7zJo182+7rt1B/OcfYO9e4ORJ/bmLF4GiRe33HzIk//Olpsqya9e8z2kdsdLTgRo18v7uq6qerZOQIFM3jx/P/3xERETurFhh//713ns+De7k5ORg/HiZTZuUlISI/N6XiYIA/wcTBboePaSOhrffcGJi9Aye3KJyeUREAJyHTCS/L55M0XKnTx/g6afN7Tt4sARlRo7Us+rS0/V6OZqff5bl559Lhp4jrebWHXfkfa5lSwkapaVJUNkxwJOWphdbj4+Xvx01apgbPxERkSvduunrw4YBPXv69PRr1qzB3r17UbFiRTz00EM+PTeRFRjgIQoGqur9AE9srJ7BU6mSfY0PI62wKxGZZ7Pl//z8+eaPVbWqBFhSUoBx42SbltU3f77UJzC67z7ggQeA5s3tt2sBXcfAECDZQvv3A927y2PHKVq//y7LxEQpNl2mjH2xZfIKRVE6K4oyJ1XLtiIiCmWbN+vvlzt2AJMm+XwI43LfV0eNGoUYi9uwE/kCAzxE4SomBnj/fVm/5hr58GcM8nz+uSyZqkrkuaFDvXesZ58FNm0CPvkkb8Hzbt2A8+dl3TGo1Ly5fRv1zEzJAnJm5059vVEj19Mzq1eXTJ7sbD2jh7xGVdWVqqoOSAzj4tVEFEY++EBfb9Ag/8YDFvjhhx+wefNmlChRAv379/fpuYmswk9uRIHKVUaNt8TGAhMnynpcnGQItGunP//KK8DGjdaOgShcmf0A36CBZOvMnCmPT5+Wpfb3oUkTmXq1YgXw66/660qWlBo/t96qb8vKAu691/05c3JcF3/WAkw2GwssExFRwfzyi3RkfPddebx3r1+GodXeeeKJJ1DUsa4dUZBigIcoUA0cCBw7Zt3xjXdJ4uLkQ6L24RGQO/off2zd+YlCmbv6OjVqmA/iDhigr+/bJ0EWrdtVbKz87jre9Tx7Nu/xMzNdB26qVJHlSy8By5fnff6GG2SpHfOVV9y3aSciInL07LNyc2LzZnk8YABQt67Ph7Fnzx6sWrUK8fHxGObYnp0oiDHAQxSo/vpL/0Dmrp5HQRg/nF24AGzfLnfu33hD316qlPfPSxQO3AVwoqLM/V5HREjRSU1MDDBoEDBqlDyOi5NaWiVK5H3t+fP2U7qyslynv2v1uF54AahcOe/zjsGczZud1/IhIiJyZc8e4LXX9MfdugGTJ/tlKBMmTAAAPP744yhdurRfxkBkBQZ4iALVxo2SQRMd7bqNeWEYP+jt2gXMni1TPP73P5nuAQD9+nn/vETh4NAh4NtvXT8fGSl1bNzZuNG+Hk5cnH3gSMu+0wooGzlm8WgBY2cZOlqHLXeMgSRm8BARkSdatdLXv/kGWLoUSEjw+TAOHz6Mjz76CJGRkXjqqad8fn4iKzHAQxTI5s2TO+U5Od4/tjHAU7kysHKl/njDBlmWLev98xKFgyNH7H+nHEVFSZq6OwkJ9pkyL71k/7yrAE+tWpKZp9XZAmSfmBjgyy/znicjA7jzTtfj0DpmrVunb2OAh4iIzHr1Vb0UwLZtQOvWfhvK5MmTkZ2djZ49e6Jq1ap+GweRFRjgIQpk+/dbd2wtwJOdLXU9NM88A7z1lqwzZZWoYOLjgalTXT//zz/AxYvmjqUFeN5/H3AsAqkFeLKy9P3q1pUpXGPGAAcP6vtmZTmvwaOqUtPnqqtcj2HhQlkau+qxyDIREbmjqkCnTsBzz8njMmWAxo39NpzTp0/j3dzizklJSX4bB5FVLOlxqihK99zVJgC2qaq61LD9PIDiAFJUVd1hxfmJQk63bt4/pvbh7OxZ++3Gu/JFinj/vEThwFWbcc0dd5gvol68uOuuVVqR5YgIIDkZuHwZ+P13yfypXh24/XZ9Xy2Dp317ueDW6vNoAd0aNcyNhyyjKEpnAJ1r8GdBRKFi2jRg9Wr98R9/+G8sAN5++22kpaWhQ4cOuPHGG/06FiIreD2DR1GUhgDOq6q6VFXVJADvKopSXFGUawG0VVV1fW7AZ6y3z00UUsqWBWbMkPVly7x//KuuAkaPlg+agwfr248f19fbtPH+eYnCQceO7luS9+1r/niOwR3t7qdWZFkL3qSnAw0b6vs/+KD+Gq0Gz6lTUvdAo9Uf6NRJXuuO2e5f5DFVVVeqqjogMTHR30MhIio4VZVurJUrAyNGyLakJMkkLVnSb8O6dOkSpk+fDgBITk722ziIrGTFFK1rAbQ1PD6bu60NJHvn/+UGg+CwbYCiKL8oivLLqVOnLBgeUZB4+WWp0zF9urkPXZ6KiJA79idPArfcItsSEoCZM2X9rrvyn7JBRK7dfbdkyrhy003A2ALe56hYEejRQ9ZjYyXAk5Ul07569wa6dwduuy3v67QuWjYbcO4csGaN/tyKFcCNNwKvv+7+/NnZ8u8jIiJytHGjXGM2bAgcPSrb+vcHxo2T61o/eu+993D27Fk0b94cLVq08OtYiKzi9QCPIXMHiqIUB1AydypWcQBnDLtqgR/H189RVbWxqqqN2bKOws78+XnvjkdFWXfHPCUFOHNG3owBYOBAYMoUmV7Sq5c15yQi+R2rWLFgr33kEcm+A2SalarqGTznzwPXXw9cc03e1/34o2TwxMbKVK6OHWV76dJAy5bSIatdO/fnz8qyzwwiIiIC5IbDHXfYbxs5Um5W+llmZiYm57ZkT05OhqJNUyYKMVYXWR4PoJHF5yAKHCkp+XfOcefDD+UDmsZmA06c0GtleJs2BWzbNrnT8uabMobateVDJBFZozBB2zfeyLstK0tqZqWl2d8hjYjQ27G/9ZYEgSIigAkT9H1efVXq/Jj1ySfAF18UbOxERBSabDb7mwQjRwJLlgCTJ9t3bvWTjz/+GEePHkXdunXRqVMnfw+HyDKWBXhyCyrPVlU1JXfTeYddSgJIAVEoOXoU+Pbbgr9+1y6ZKmGzSWbNP/8AL75orp1yQXz1lSyHDdM/cBYrJu2ViSg4rF4ttQ2io/U6O5rYWAn+aKKjgVmzJHAMyBRN7e+AWVu3SjCbiIhI89lnekb48uUS2Ln/fv+OKVdOTg4m5N7YGDNmDCIi2EiaQpcl/7sVRWkDYIeqqjsMBZaXAKhu2K04u2hRyImMlOBMQZ09C+zdC+zZA/zvf3o2T9eu3hmfI632TqtW8sEQkAAP3/iIrJOd7d0W4ytXyt+OuDj5+/Hvv/pzMTGSMq+JiACqVdMfnzkjF+WeuHAB+OWXwo2ZiIiC38SJkmWuKMADD8i2Z58FunTx77gcrFq1Cr///jsqVaqEBznFmEKcVV20ZgP4VFGU7QD+VlU1RVXV87nb2uQGgMZ7+9xEfvf66/p0iIJo2xaoWxdo0EAeX7nim4J0JUvKB0NApmeVK2f9OYnCVVaWfZaNt2h1ebSiloAEeLTg7dSp+jZNWprn5ylWzLsBKiIiCj7r1wNjxthva90aeO45/4zHBVVVMW7cOADAqFGjEBMA08WIrOT1T465WTnVXTy33tvnIwooX34JDBpU8NdrhU41GRl6VysrFS+u3+WfOBFg60gi63g7wJOYCKSm6kGXhAT9OWOAJy5OlsZzp6Z6bxxERBTa0tIkiLN1q/325GRpHpCUFBD1dox++OEH/PjjjyhZsiT69evn7+EQWc6/veqIQlFBP7j99JP9dIpPPpF0V2cFVb1JS6M1viFfd5215yQKZ94O8DgWbDYWTI6MlExAI20KZrlywJ13em8cREQU2l55xT64c++9wLx5nhXq9zEte2fYsGEoWrSon0dDZD0W2iDytoJ+cGvWTD6M5bZwxF13ydLqOyHLl+fd9uij1p6TKJxZMUVLa0t7/fX2HfBWrQImTcq7/wsvmGuJTj6lKEpnRVHmpDKziogCjapKUX/NwIHAwoUBHdzZvXs31qxZg/j4eDzxxBP+Hg6RTzDAQ+Qt2l10LUBTEDk5kskDyBtm1arSBYeIgs/Agc63ezvAEx2t/52IjbUPCo8aBdx9d97XWNmdjwpMVdWVqqoOSExM9PdQiIjs7dgB7N6td1udNct+SnAA0jpn9e/fH1dffbWfR0PkGwzwEHmLsRWxt1SqxAAPUTAzdrDSeDvA8/rreqtzRbF/LjERWLvW+evi42W5bp33xkJERKHpu+9kef/9EuQJcIcOHcLixYsRGRmJUaNG+Xs4RD7DAA+Rt2iFTL15jEuXAq5YHRGZNGoU8N9/ebd7O8DTrRvQqJGsL11q/1xMDDBnjvPXaX9b2rb13liIiCj07Nqld2ps08a/YzHpzTffRHZ2Nh566CFUqVLF38Mh8hkGeIi8JTNTr4Mxe7Z01PKUYxbQzp3WtFMmIuvFxwOVK+fd7u0AT0yMniZftar9c9p5LlzI+7oo9lkgIqJ8qKrUcWvQQB5HRQHt2/t3TCacOnUKc+fOBQCMcWzlThTiGOAh8paMDGDkSKltMWgQcPGi58fQAjydOsny8OG8Uy6IKDg4C6oA3g/wJCQA/fs7f07L0tm/H/jww7zPt27tvXEQEVFomTFDz9wBJNgTwEWVNW+//TbS0tLQqVMn3HDDDf4eDpFPMcBDVFg5OcDixTLNITpaAjyA3orYjB07ZKlN0Vq2TJbO7v4TUXC4/37n270d4ImIsO+cZRQZKedKT5cOW0ZFigDsKkJERM6kpADDhsl6+/ZyI/PJJ/07JhMuXbqE6dOnAwCSk5P9PBoi32OAh6iwbDZgyhTJ2LlyRd+enW3+GCdOyPLNN4ExY1h3hygUNGzofLsVbdLzU7QocPmynmKviY0Funb13TiIiCg4qCrQr5/+eMmSoLk2fffdd3Hu3DnceuutuPXWW/09HCKfY4CHqLBycuQu+ZEjgLEF45Qp5o+Rnq6vjx/vvbERkf9oXaoc+TrAk5AgBds53ZOIiMz46ivpmlWqFPDvv3KjIAhkZmbizTffBMDsHQpfDPAQFdby5foUq9tv19NXt241fwwtwFOhgnfHRkT+ExEBOLvA9EeA5/JluSNLRESUn82b9ULKTzwRVNemixYtwr///ovrr78eHTp08PdwiPyCAR6iwvrtN/sMnEcf9bztsFZcuW9f742LiPwvt4uHnUce8W2A588/gbNngbg4352TiIiCS0oKUK8e0KKFvi2IrktzcnIwYcIEAEBSUhIiPKmFSRRC+D+fqLC04M4zz8iyalXg8cdlPSXF3DHKlQNmzXI9pYOIgtOoUXkzZ06e9G2ABwDS0qSoMhERkaNz54AmTYDdu+VxhQryvhFEzT5WrFiBffv2oXLlyujZs6e/h0PkNwzwEJnlanpDWpos//pLliVKAK1ayfrixeaOnZEBdO4MJCUVboxEFFgqVtT/RhhFRfluDIsXA88/HzQFMsOZoiidFUWZk5qa6u+hEFEou3wZeOMNIDFR6rOVLCmZngDQoQOwdGlQZX2qqopx48YBAJ566ilE+/omClEAYYCHyIyzZ4EBA5y3FL50SZanT+vbtGKm69ebO77NJh++IiMLN04iCixFi8rfiAUL5HFOjix9efH5wAO+r/tDBaKq6kpVVQckJib6eyhEFKoyM4Hu3YH//Q+4cEHfXqIE8PffwOrVQPPm/htfAXz//ff46aefUKpUKTyuZdEThSkGeIjMKFcOeO89YMYMfdvOncC8edIeHQDeflt/7qqrZPndd+aOn5Xl2zv6ROQbRYvK34hp0+TxlSuy9Ecwl/UIiIjCV1oa8NxzQGwssHatvCfUry/PFS0KfPmllBkIQlr2zrBhw5CQkODn0RD5Fz9REpmhFUE2euEF2b52rTy+/nr9udhYz4/Pu+tEoSchAdizR5/iqQWE/VFvi23SiYjCk6pKJvrChfq2NWuAdu38NyYv+fXXX7F27VoUKVIETzjLtCcKM7ydR+TOzz/n3XbyJLBypR7cado07z5Vqpg/h83GDB6iUKSqwL336o8vXgS++cY/v+8M8BARhaeVK/XgTrNm0go9BII7ADB+/HgAQP/+/VGqVCk/j4bI/xjgIXJn5059XZuiVa+e/T4bN+Z93aZN5s/BDB6i0KQVqdyxA3j2WWDcOEmF97VBgxjgISIKR59+KrXYAGDiRODHH4FbbvHvmLwkJSUFn3zyCaKiojBq1Ch/D4coIDDAQ+SOcbrVwYNSa+fkSft9nHUaKFdO7pKYkZPD+hhEoahBA339yy+lbpc/2pUXK8a/MURE4URVgc8+Ax5+GEhPB+67Dxg50t+j8qo333wTOTk5ePjhh1E5iFq6E1mJc0KI3NGCNw0ayF347dvtn1+zxvnroqOBxx6zdGhEFOCMWTOqCrz0EnDddb4fxw03AOXL+/68RETkH488Anz0kaw/8ADw8cchlcn533//4f333wcAjBkzxs+jIQocDPAQmTFqlHTE2rAh73M1a7p+XX5vpBs2AMWL6x0MiCi0qSpQtqx/Omj17u37cxIRkX+sXq0Hd/r2lU6OIRTcAYBp06YhPT0d99xzD+rWrevv4RAFDOZrE7midb1JTwdefdX+jbFpU6BECeDFF4EaNQp2/FatgPXrCz1MIgoSIXZxTUREAeTcOWDsWJmG1aWLbHv+eWDuXOnoGEIuXryIGbl1MZOSkvw8GqLAwgAPkSsdO8oyM1Pq8EyaBNSpI9t+/lkK1L3wQsGO/fXXshw9WpZaMImIQs/YsbLMzOTvOhERed/y5UDJklLIf8oUIDsbeOghKe4fgubMmYPz58+jZcuWuCVECkYTeQsDPESuaIWUVVWKk7ZqBcTE6M8vXlzwY2dm2j/mnX2i0NWrlyz37pXuJURERN5y7Jj+PgMAnTpJ1s6CBSHZoTUjIwOTJ08GACQnJ/t5NESBhwEeIld27AByi7f9P2MRt8K0OtY6cw0eLEsGeIhCV+XK+sX3ggX+HQsREYWOrCygTx/g0iWgbVtZrlwpdXf8Ue/NBxYtWoRjx47hxhtvRPv27f09HKKAwwAPUX6O/197dx4fRX3/cfw9OTgDCQgIihwR8KiiciieoIJHkaKCUg8sagVEFFAgaFv7q9KWGEDwQoFWpFYreJZqRbCo9cSIgojiASJeyBVuyPX9/fHZ7SYhARJ2d/Z4PR+PeczszOzshxmymXzm+/18fyjfpeKqq6QePcoPfbwvaWn2y7ei4BOV4LDFdNsAElf9+tIZZ/gdBeKA53l9PM+bvmXLFr9DARCrCguld9+Vpk6VOneWXnnFumfNmJFwtXYqKikp0T333CPJau94PCAF9sIoWsC+pKTs3bpm0SLp3nsP7P3160s7dthoWWWVlloLoVmz7DW/oIDEVlwstWwpffut35Eghjnn5kma16VLlxv8jgVAjNi2TXr4YSk/3+4X58+XCgrK7zNtmtS6tT/xRdELL7yglStXqnXr1howYIDf4QAxiQQPsC8lJdI//ykNGVJ+/ahRB/b+WrWkf/1Luvrq8usLC20UruCwjrTgARLb4sXShRdKZ53ldyQAgHhQWCi98Ybdcy5fXn5bhw7WMnTdOunMM6XLLvMnxihyzik3N1eSNHr0aKWl8WcsUBl+MoCqDB9uLW9eeungjnP77XsneIqKrJsWLXeA5JCSInXtuvd3AQAAZTkn/fe/VqdxxQpb17Sp9Mc/Wqvwjh1t4I8ku4d87bXXtHjxYjVp0kTXXXed3+EAMYsED1CV9u2lm28+uGNkZlqXjGBCZ/Nma7lTWGite4Itd5LslzSQdOrVk372M7+jAADEsqIiK5D8+OOhdd27S3PmSM2a+RdXDAi23rnllltUr149n6MBYhcJHqAq4Rh9oF07m7dqJa1ebUXwSkstwRMstPzNN3TRAhJdUZEVXQcAIOiZZ6QHH5SOOko6/XRpwgTpk0+shuNtt0lDh0qHHhoalCNJffjhh5o/f77q16+vm266ye9wgJjG3SZQGefCk+Bp1crmP/5oRfIkacEC6YorrJltaakVxRs//uA/C0DsuvLKUMIXAJDc9uyRJk+WfvMbu+dctMgKKUtSnTrSwoVSt27+xhhDgiNnDR48WI0bN/Y5GiC2JXc6GKhK796hpyUjRhz88S64wBI6krR+vc3T00Pr/v3vg/8MALGre3drwQcASE7bt0t5eVKPHlKTJtIdd1hyZ+hQq892/PGW1Jk7l+ROGV999ZXmzJmj9PR03XrrrX6HA8Q8WvAAlfnpJ2tdI0ktWhz88dq1s2ROp07S7t22Li0t1KoHAAAAienHH6Wf/1z68MPQug4dpIkTpT59/IsrDkycOFGlpaW65ppr1LJlS7/DAWIeLXiAyqSnW82Mtm2lQYMO7lilpfZUpm9f6e67pV27bL3nhVrw3HXXwX0GAAAAYs/KldKpp1pyp107K5j8wQfSsmUkd/Zj3bp1evTRRyVJY8eO9TkaID7QggeoTK1aluDJyrLidgfD86wJ7ldf2TCX+flSr162rbjY5l27HtxnAAAAIHZ88YWUkyP9859SSYl08snSv/5l94I4IFOnTtWePXt08cUX65hjjvE7HCAu0IIHqMwbb1iCJ1zDl8+fb/OsLGnr1lCLnfbtLZnUoEF4PgcAAADRtXu3FU6W7KHehAlSx47Sc89Zcqd3b+k//yG5Uw1bt27VQw89JEnKycnxORogftCCB6jo5ZetwHJxsRSo2n/QXn/d+lo3a2Z1d+rUsfU33ST95S/h+QwAAABE16efWsvsjRuljAx7OBgcUOOKK6Sbb5ZOOSXphzqvrkceeURbtmxR9+7d1Y2i08AB45sGCLrzTptfeKFNRUXSueeG59ibNtmw6w0blk/wSKFhMQEAABA/xo+XTjpJ+u47a8WzYYMldxo0kJ55RnriCau/Q3KnWvbs2aN7771XEq13gOri2wYIeuml0HKwyHI4ffqpPdV58kk7ftDJJ4f3cwAAABA5S5dKv/iF9LvfWdesyy+XPvrICiqvXGlJnksv9TvKuPW3v/1NP/zwgzp27KgLLrjA73CAuEIXLSSfMWOkUaOkww4rv945+yXduXOoyHK45OaG6vBIdnwAAADEl7//Xfr1r63FTmqqNHu2dOWVfkeVMEpKSnRPoERCTk6OvHDVwwSSBC14kHwmTpTefFNat27vbT16WKIn3C14xo6VgtX/TzxROuKI8B0bAAAAkbVunTRypHT11Zbc+dWvbKQskjth9fzzz+uLL75QmzZtdPnll/sdDhB3SPAgOWVmWnHjCRPs9U8/SUuWSO++G9rngQfC+5nBkU4gCiMAACAASURBVBN4EgEAqITneX08z5u+ZcsWv0MBEPT119JFF0nNm0tTp0ppadJDD0mPPiq1bet3dAnFOacJgXvzMWPGKC2NziZAdZHgQXJxzuY7d9ov7Ntvt9c7dpTfb9cuG0UrnI47zuYkeAAAlXDOzXPODc7MzPQ7FCB5lZZaAqd7dyuQ3Lat9OKLNkBGjx7SokXSjTdyPxcBixYtUn5+vpo2baprr73W73CAuERaFMnlrbdsPmSIdZUKmjy5/H5nnCFNnx7ez+7Xz+Zdu4b3uAAAADg4RUXWjX/WLOnzz8tvO/106dlnpWbNfAktWQRb74wYMUJ169b1ORogPpHgQXLZsMHm69dLCxbY8vnnS6+8Yk9kzj7btqWnS4ccEpkYGBYdAAAgunbtkt5/35I07dpZV6v1662r1UcfWTf91att32bNpLw8qUkTac0aaeBAKSPD3/gT3AcffKAFCxYoIyNDw4YN8zscIG6R4EFyadxYWrFCOvZYqU8fad486b//tW3Nm9t861YrngcAAID4tnOn9NJLNuBFMIHTsKF1nV+xQiooCO3burX04IPSeefZwz5ETXDkrCFDhqhRo0Y+RwPEL2rwILns3CnVqycNGybVrm3rdu2yef360vjx0pYt0ssv+xcjAAAADs6CBdLFF1srnMsus+RO27ZSq1b2MO/tty2506uXDbzx3nvSypVS794kd6Lsyy+/1NNPP6309HSNGjXK73CAuEYLHsSvH36QWrSo3nu2bZMaNLCnM9nZtq51a6lRI+mww6xVz4wZ1hQXAAAA8aGkxLpcLV0qPfaY9OSToW1du9pw5jfdZMmbb7+1wTbq1pU6daJgss8mTpyo0tJSDRo0SIcffrjf4QBxjQQP4tff/y6NHGl9qI891mrbnHXWvt+zbJnUt68tr1pl8+xs67qVmip17Gg3BgzLCAAAENuKi6VHHpHmzLEaOoWFoW116ki//a00aJBUMWnQsqVN8N0PP/ygRx99VJ7nacyYMX6HA8Q9umghfhUUSHv22PKnn0r5+bbsXGg49IrGj5dq1Sq/bvVqe+oTlJlp3bUAAAAQe0pLpfnzpVNOkYYPl954w5I7TZtaa53f/EZavtzmtAiJaVOnTlVhYaEuvvhiHX300X6HA8Q9mikgfhUUSD/+KB15pL3escPmc+ZYV6xf//rAjxV8LwAAAGLPypXW9erDD6WPP5a++87Wt2ol/fnP0gUXWItsxI0tW7Zo2rRpkqScnByfowESAwkexKeiIum556yWTnDEqzvvlH73O+mtt0JJH+fK96vOzNz7WJ06SRMnRj5mAAAAVM05q6Ozfr0NjFG3rg1hPmOGtdIpq3VrafBg6ZZbGMI8Tj388MPaunWrevTooVNOOcXvcICEQIIH8SnYzSo9PdTHevx4W3f//dLkydYvOz091F3LOSkwBGM5zkndu0cnbgAAABjnbPTS77+3B3czZ1rx48rUry/98pfSRRdJ7dtLRx9t9RMRl3bv3q0pU6ZIksaNG+dzNEDiIMGD+FS7ttXfKSqy1y1b2sgIQR99JC1ebMvz5tnoWBdcYH21K3riicjHCwAAkMycs1Y4994rvfiiPYjzvL3rJmZmSs2bW0Jn927rdjVwoCV3Gjb0J3aE3ezZs/Xjjz/qxBNP1Hnnned3OEDCiEiCx/O8TpK6OOeml1nXX1KBpCxJq5xzSyLx2UgSe/ZIp50mvf22vf7uO+m446zQcr9+duPw7be27Re/kD74QNqwQbr7bumuu8ofi4JuAAAA4bFnj7Rkid2jvf22jWC6a5dNmzaV39c5617VvLl00klWP7FnTymFcWASWUlJifLy8iRZ7R2PYeqBsAl7gsfzvJ6Shkh6v8y6bEm9nHNDAq/nSros3J+NJLN+fWj57rutBk/XrjZc+jPPSFu3hrb//OeV198BAADAgdu2Tfr8c+nLL23avNmSOnv2SCtW2MO24CinFTVtKg0bJt14o9SsmY2GRTerpPPss8/qyy+/VHZ2tvr37+93OEBCCXuCxzm3MJDQySqzuqes9c7/eJ7XiVY8qLGMDGnkSOn99200haVL7WnPhReG9snPl7KzpVWrpHXrbLrlFv9iBgAAiAclJdLGjdJPP9m0erX07rvSO+9YEqdit6qKjj1WOv10a23dpYs9ZKtVSzrkECmtzJ8fJHeSjnNOEyZMkCSNHj1aaWlUDAHCKVo/UVmSNpZ5vUlStiQSPKiZq66SLrvMngItXy5dfLFUWBgqvvyf/0jnnCMNGmQte4Lq1vUlXAAAgJjknD0Ee++9ULeq/PzQKKUVpadLHTpYoeP27a1VTu3adg/WqpV06qlSo0bR/Tcgbrz66qtasmSJmjVrpkGDBvkdDpBwYi5l6nneYEmDJalVq1Y+R4OYtG2bdMYZdkMhWe2dL7+UduywVj2S1K2bzcv26T3pJHsvAABAonPO6hEuX24JnJIS6xJVWmoPxb74Qvr4Y9u+YcPe7z/kEOtG1ayZ1KKFdYM/9VS7n6pTJ/r/HiSEYOudESNGqC4PXoGwi1aCJ1hcOaixpFWV7RgozDxdkrp06bKf9p9IOsXF0o8/2o1GRfXrSx07Wlet4I1H167SAw9InTvbtv01KQYAAIgn27ZZF6pVq2xaudKSNsuXl69HuC+ZmZa4CXar6tbNRq8Cwig/P1+vvvqqGjRooGHDhvkdDpCQopXgmSMpt8zrLOrvoEbS06U//rF8rZ2ynnrKkkA33CB99pk1IT7//OjGCAAAktuuXVa7JlgDcPNmqzeTnh6a0tJCy54nFRTYKFMVp4ICu7cpKQnNS0qsFc7atZW3vgk65BDp+OOt61Ramj0ES0mxWNq0sW3HHSe1bFm+1TMQAbm59ufg0KFDlZWVtZ+9AdREpEbR6iUpy/O8Jc65hc65As/z5ga2SeWTPUD1fPGF1eCpzLZtNtSmJB11VPRiAgAA8W/3bmv1snOnJWl27qx8Cm7bscOKEQcTOcEpml3C69SR2ra1KTtbatfOkjbHHWfdq0jcIAZ8/vnneuaZZ1SrVi2NDJZUABB2ERlFS9LCKtYDNRe8WZo1S7rvvsr3OeYY6Q9/iFpIAAAgRpSWhobr3r7dWs2UnQoK9l5XcXtVhYWrKz3dkiuHHmpTo0YWX3GxVFRkU9ll56SsLOsWVXHKzLTjpaaGprQ0mw4/3B5spaSEJ24gQiZOnCjnnK655hoddthhfocDJKyYK7IMVOnmm0PLGRmV7/Pvf0t5edGJBwCAMPI8r4+kPu3btJG+/z4yH+KcJRqC3Xyqs1zT9+3vGIWFocTM7t2h5eqs273bEiYHKz3dEir16lU+1a279+vGjUOJnEMPtcROo0a0nAECfvjhBz322GPyPE9jxozxOxwgoZHgQfxo1szmV11V9U3T1q12cwYAQJxxzs2TNK+L592gww/3O5z4VLu2TfXrW5IlK8vmBzrVrUtiBgizKVOmqLCwUP369VOHDh38DgdIaCR4ED/at5dOOEF6/PGq92GULABAvEtPl5o0idzxU1NDhXYPZLk6+9bkfbVqhRIztWtbTZmyrw90Xa1aJGeAGFNQUKBp06ZJknJycnyOBkh8JHgQPwYPttEiAABIZB07Svn5fkcBAAft4Ycf1rZt23TOOeeoa9eufocDJDwqsiF+BIfxBAAAABDTdu3apSlTpkiSxo0b53M0QHIgwYP4MXy43xEAAAAAOACPPfaY1q1bp5NOOkk9e/b0OxwgKZDgAQAAAACETXFxsfICI9uOGzdOHvWxgKggwYPY989/St9843cUAAAAAA7AM888o1WrVunII49Uv379/A4HSBoUWUbse/11G+4UAAAAQExzzik3N1eSNGbMGKWmpvocEZA8aMGD2FdUJNFvFwAAAIh5CxYs0IcffqhDDz1Uv/rVr/wOB0gqJHgQ+5Yvt7lz/sYBAAAAYJ8mTJggSRo5cqTq1KnjczRAciHBg9h3+eU2z8nxNw4AAAAAVVq8eLEWLVqkhg0b6sYbb/Q7HCDpkOCJFbROqZrnSZdeKr32mt+RAAAAAKhCsPbO0KFDlZmZ6XM0QPIhwRMrhg9npKh9GTRIOuYYv6MAAAAAUImVK1fqueeeU61atTRy5Ei/wwGSEgmeWLFokbRr1/73+/pr6fbbIx5OTPnsM6lPH4k+vAAAAEBMysvLk3NOgwYNUosWLfwOB0hKJHhiwc6d0qefSu+8s/99N22S1q6NfEyxYudOKSvL7ygAAAAAVOG7777T7Nmz5XmeRo8e7Xc4QNIiwRMLXnnF5tdeu/99O3dOroTHKadI7dr5HQUAAACAKkyZMkVFRUXq37+/2rdv73c4QNIiwRMLatWSliyx5dLSqvcLFmLeti3yMcWK5cv3fU4AAAAA+Gbz5s16+OGHJUk5jHoL+IoETyx44QUpO9uWX3+96v2Ki20+e3bkY4olwX83AAAAgJgybdo0bd++XT179lTnzp39DgdIaiR4YsH06VKDBlJGxr6HS//+e2nsWOm446IXm9+OP97vCAAAAABUYteuXZoyZYokWu8AsSDN7wCSnnPS5MlSSor01ltWQHnqVGnEiL33bdNG+vzzUGufZNC0qdS4sd9RAAAAAKhg1qxZWr9+vTp37qxzzz3X73CApEcLHr99+WUoYXPIIdKyZdLIkVXvf8QR0vz5luhJNM5Ju3fbSGGS9I9/SAMGSH37+hsXAAAAgHKKi4uVl5cnSRo3bpw8z/M5IgAkePw2YICUmmrLGRnSHXdIV18tXXJJ+f1GjpRyc6U6daSlS0NFmePV738vrVkTer1li7ViuvRSS3Q5Z//e//s/30IEAAAAULm5c+dq9erVat++vS6p+LcLAF+Q4PHbkCHSRRfZcmam1KiR9Pe/S88/X370qKlTQ4mgvDzpiiv2Xa8nVnzxReVx3nWX9MYbodcbN9o8mPn/7W+ljz4KnRsAAAAAMcE5p9zcXEnSmDFjlBr8OwWAr0jwxJrNm6U777Tl77+3LktBe/bY/LzzbL5sWXRjq4kzzpA++KDybdu3S++/L/30k7Rhg6176SWb/+lPNj/55MjHCAAAAOCAzZ8/X0uXLlXz5s01cOBAv8MBEECCx28lJeVfP/igNGyYLU+aJN1/vxVf7tQpVHg5I0Pq1SvU6iWWNW0qTZwYev3TT6HlrVulxx6zOkTr19tIYhdcIF17rW0//XSpR4+ohgsAAABg34Ktd0aNGqU6der4HA2AIBI8fvr0U+mmm8qvGzbMRo2aOVOaMkVavdpawaxcKdWvH9pvwYL4SPDcfLN0zjnS669Lf/iDdOihtr57d2ncOEtopaRY4ufZZ6Ujjwwlvd58U2rXzr/YAQAAAJTz3nvv6bXXXlNmZqaGDh3qdzgAyiDB46f58ytfn5ZmNXZOO02aNs3W5efvvV9BQeRiC6eJE6Vvvw2N/LV7t/37JCsovX27FVTu2VO68ELp0Uf9ixUAAABAlYKtd2688UY1bNjQ52gAlEWCJxr69au80PCoUdZNqTL16knXXGPLnTpJRx9dfrtze3fvikXOWaHlPXukJ56QHnpIqlvXWu088oi1UNq+XWrb1vbv3du2MXoWAAAAEFM+/fRTPffcc6pdu7ZGBMtHAIgZJHgibdMm63q0YIG0a5d1y1q7NrQ9I2Pf7+/XT8rKqnzbjTeGL85I2LpVeu89G+K9oEA64QSrHSRJDRtKZ59t9YVuuWXvYeF///voxwsAAACgSnl5eZKkQYMGqXnz5j5HA6AiEjyRFkxo5OZaq5y//MXq0UhWRDk4LHhlli6VTj1VeuWVqvcpLAxfrFV57bWave+mm6RZs6QTT5QWLbKWOU2ahJJW7dvbUPBr19rw8AAAJDHP8/p4njd9y5YtfocCAHv59ttv9fjjjyslJUWjR4/2OxwAlSDBE2nOWR2d4DDgkyZJAwdK33xjXZP25YILpG3bpNTUyrffeafUoUN4463MyJE1e1/9+vZvT02V/vUv6cMPLbEzfHgooeN5lsQ666zwxQsAQBxyzs1zzg3OzMz0OxQA2Mu9996roqIi9e/fX+0YCAWISWl+B5DQSkosuZGSIi1bJmVnS6tWWW2dqVNt3b40ayadeWbV28eMsS5OkZZSwzzgaadJl11mNXcmTZJuuMHWn3lm+ZZLb7998DECAAAAiIhNmzZp+vTpkqScnByfowFQFRI8kZSWZsWRU1OlTz6RXn1VqlVL6tZNuvZa6T//2ff7u3Xb9/aMDDvm9u37r+VTU/PmWcub4mJLUK1Zs+9uZWXt3i3VqRN6PWyYzatqkQQAAAAg5jz00EPavn27zjvvPHXq1MnvcABUgS5akbZ6tZSZKR17rCViGjWy182a2Twc/vKX8BynMsuX23zFCquV07v33vucfnrl733mmVAyqGfPyMQHAAAAIGJ27typqVOnSqL1DhDrSPBESlGRJTWOOEL67jtbN2SIdNFFNoLU/PkH3hJmX2bNCg0xHgn160u//a20eLF09dXSunXltxcXWxer4JDuQe+8U744dKDiPgAAAID48eijj2rDhg3q2rWrzj77bL/DAbAPJHgiZdMm6brrpJdftmLJktS4sY2k1bChdP754fmc/v1tOPJwKy21+X//a0WR771XOvRQacsWKxwd9OCDNv/b30LrPM/q77RpE/64AAAAAERFcXGxJk6cKMla73jheEANIGJI8ETKhx9aQufII6Wjjiq/LS3NEj/hUK+edZ8Kp9JS6bDDrBXS4sXWnWzFCqlVK6llS+nbb0P7jhwpbd5shZh797a6O0Hjx4c3LgAAAABRM2fOHH399dfq0KGDLr74Yr/DAbAfJHgi5fLLbfjvSPM8S7qE044d1hVr61bpjjvsM845xxJWP/uZdb/q2NG6b0k29PnOnTZqWJs21tJHkq66KrxxAQAAAIgK55xyc3MlSWPGjFEqA6UAMY8ET6Tk5VlXrGiYMSN8xyoulr7+2pabNJFefNGWX33VEjwNGkgDBkgff2xJnRNPtO21a0tvvmmJoT/8IdTFCwAAAEDcefnll7Vs2TK1aNFCAwcO9DscAAeABE+k3H139D7rxhvDd6wZM6x1zuTJ9vrdd0PbateW0tOlww+3bmfPP29JnaAdO6Q775QuuSQ8BaQBAAAA+GLChAmSpFtvvVW1a9f2ORoAByLN7wAS1u9+F73P2rHDCh//9FOoe9T331sLooyM6h1r9GibB4c1/+ij0LZzzpHOOku6/XbrntW3b/n3zp4t9ehhNXsAAAAAxKV33nlHb7zxhjIzMzV48GC/wwFwgGjBEyllR5qKtKZNpSeflJo3D62bPFmaO9da0hQUlN9/7VqbP/uszb/5Rpo2zZYnTLARwI4/PlRsOcjzrAVPvXpWb6eigQNtWHgAAAAAcStYe+emm25Sw2iVnQBw0EjwRMKxx0p79kTv8372M+nPf7blzz6T/vpXadIkG6Zdkv70p9C+X3xho2GVlkp//KMle+64w94X3N6okS3vq5tVCv91AAAAgESzYsUKvfDCC6pdu7ZuueUWv8MBUA38lV4T8+ZJl1669/qdO0PzIUOiF0/79tLy5VYI+dhjpeuvD21bs0ZasMDm27dbcWTJRr5assSSPVlZUosWtr116+jFDQAAACCm5OXlSZKuu+46HRos/wAgLpDgqYnbbpOee85GynrkEVu3aZP0y19aq5fcXKlOnejFk5EhffWVdMopoa5hv/mNtdJp1UoqLLRtDRpI/fpJn3wibdsmjR9v+z74oL2/fXsbKQsAAABA0lm7dq0ef/xxpaSkaHSwNieAuEGCpyZGjJDOPVcaO1YaOtTq0bz7rrWSkawmTrRlZ1vh4xEjLMkzfnyoi9WaNTZ8uWQtj4491pZzcqRdu2xo9DFjpB9/tCQVAAAAgKQzefJkFRcX6/LLL1d2drbf4QCoJkbRqq7duy0hsnBhKIHyy19Kn39uSZVu3aQmTfyLb8qUvdc995x03HGWAAoWSfvHP6S0NJskqUMHadkyqW7d6MUKAAAAICZs3LhRM2bMkCTl5OT4HA2AmqAFT3VlZUmvvGLLpaXSl19K3btLffrY1LTpvosT+6FXL6uxU7YC/oABe+93/PHRiwkAAABAzHjwwQe1Y8cOnX/++TrxxBP9DgdADZDgORBbtkhFRbZcWirNnm3LnicdeaQNUd6ihX/xAQAAAEAN7dixQ/fdd58kady4cT5HA6CmSPAciOxsaeBAS/Q8/LB0yCHlt7/9tnTJJf7EBgAAAAAH4a9//as2btyok08+Wd27d/c7HAA1RIJnf/LzbYSsp54KDSdekXPSYYdFPzYAAAAAOAhFRUWaNGmSJGu948VauQkAB4wEz/4sXmyjUM2cKXXsKJ16qt8RAQAAAEBYPPXUU1qzZo2OOuoo9e3b1+9wABwEEjwVLV5sdXaCUlOlVq2k66+Xli61VjwAAAAAEOecc8rNzZUkjR07Vikp/HkIxDN+gis6+2xL6jhniZ677vI7IgAAAAAIu5deeknLly/X4YcfrquuusrvcAAcJBI8ktXY2bFDKiyU7rtP6tJFSkmxRM9jj/kdHQAAAACE3YQJEyRJo0aNUu3atX2OBsDBIsFTWmqjYmVkSN27S8cdJ731ljRhgvTJJ1LPnn5HCAAAAABh9dZbb+nNN99UVlaWBg8e7Hc4AMIgze8AfPPZZ9Yda9Mm6f33pVmzpA4dpK5drfVOTo7fEQIAAABARARr7wwfPlwNGjTwORoA4ZC8CZ6nnpKeflo64gibunTxOyIAAAAAiLjly5dr3rx5qlOnjm6++Wa/wwEQJlFN8Hie119SgaQsSaucc0ui+fnl/P73vn00AAAAAPglLy9PknT99derWbNmPkcDIFyiVoPH87xsSb2ccwudc09Luj1anw0AAAAAkL755hs98cQTSk1N1W233eZ3OADCKJpFlnvKWu/8j+d5nSru5HneYM/z8j3Py1+/fn3UggMAAACARDd58mQVFxdrwIABatu2rd/hAAijaCZ4siRtLPN6k6Tsijs556Y757o457o0bdo0asEBAAAAQCLbuHGjZsyYIUkaO3asz9EACDeGSQcAAACAJPDAAw9o586duvDCC3XCCSf4HQ6AMItmgqegwuvGklZF8fMBAAAAICnt2LFD9913nyRp3LhxPkcDIBKimeCZI+nIMq+zfB1FCwAAAACSxMyZM7Vp0yZ169ZNZ555pt/hAIiAqA2T7pwr8Dxvrud5PQOrcqP12QAAAACQrIqKijRp0iRJ1nrH8zyfIwIQCVFL8EiSc25hND8PAAAAAJLdk08+qbVr1+qYY45Rnz59/A4HQIRQZBkAAAAAElRpaanuueceSTZyVkoKfwICiYqfbgAAAABIUC+++KI++eQTtWzZUldeeaXf4QCIIBI8AAAAAJCgJkyYIEm69dZbVatWLZ+jARBJJHgAAAAAIAG9+eabevvtt9WoUSPdcMMNfocDIMJI8AAAAABAAgq23hk+fLgyMjJ8jgZApJHgAQAAAIAE8/HHH+vFF19U3bp1dfPNN/sdDoAoIMEDAAAAAAkmOHLW9ddfr6ZNm/ocDYBoIMEDAAAAAAlkzZo1evLJJ5WamqrbbrvN73AARAkJHgAAAABIIJMmTVJJSYmuuOIKtWnTxu9wAEQJCR4AAIAI8jwvy/O8np7njfU7FgCJb/369Zo5c6YkaexYvnaAZEKCBwAAIIKccwWSVkk6xO9YACS+Bx54QLt27VLv3r11/PHH+x0OgChK8zsAAACAWOV5XidJXZxz08us6y+pQFKWpFXOuSV+xQcAZW3fvl3333+/JCknJ8fnaABEGy14AAAAKuF5Xk9Jt8sSOcF12ZJ6OecWOueeDmwHgJgwc+ZMbd68WaeddprOOOMMv8MBEGW04AEAAKiEc25hIKGTVWZ1T1nrnf/xPK+Tc25JoGVPWQXOuYWRjhMAJKmwsFCTJk2SZK13PM/zOSIA0RbTCZ4PPvhgg+d5ayL4EU0kbYjg8VE5zrs/OO/+4dz7g/Puj0if99YRPPaByJK0sczrTZKyJS0JtOipSk9JnTzPy3bOraq40fO8wZIGB15u9zzvR0lbqjhWZiXbYu3/e2Ux+n3c6rz3QPfd33772l7VtqrWJ8M1jtb1rc7+0brGYb2+ffv2DcdhkuEa8zMcW8fke7pqB3b/45xL2klSvt8xJOPEeee8J9vEuee8J9OUaOddlnQZW+b12AqvH5HUPwKfO70622LtvO8rfr+OW533Hui++9uvutdxP+sT/hpH6/rG4jWOteubLNeYn+HYOibf0wc/UYMHAADgwBVUeN1YNkJWuM2r4bZYEakYD+a41Xnvge67v/1qch3j4fpKkYkzWte3OvtzjWPnmJG4xlzf2Dom39MHyQtkmpKS53n5zrkufseRbDjv/uC8+4dz7w/Ouz8S7bwHuk5lOefuCbzOkpTrnBsSeL3AOdfLzxgDcSTUecfeuMaJjeub+LjGiS8WrnFM1+CJgun73wURwHn3B+fdP5x7f3De/ZEw5z0wilYvSVme5y1xNnJWged5cwPbJCnXxxDLSpjzjipxjRMb1zfxcY0Tn+/XOKlb8AAAAAAAACQCavAAAAAAAADEuWTvogUAAJDQAnWDusiGc1/lnFvoc0gIozLXt5OkJVzfxFP2GgfrgSG+Ba7pYFmR/lXOuSU+h4Qw8+vnlhY8AAAAiS1Y8DFflgRAYrlc9gfiPZJy/A4G4eecK5AlAg7xOxaEzWDZUNtPSxrgdzAIP79+bmnBAyQJz/P6SyrgyR6AcPI8r5OkLs656WXW9ZcNJ54lnkwelHCcX+fcQs/zsiXlBEf/QmwI0/WdXuZYqyIYLmqA78jkUo3r3bVMq47s6EeKmor1n+mkTPDE0gVIJoHzLkldJb0fyFgjCgJNBAdIesrvWJJJ4LxfLrvhzi77iwCRExjd/+8MdQAABBxJREFUKEv2PS+SmpETONdDJL1fZl22pF5lhhGfK+kyfyKMb+E8v865VZ7n5XqeN9c5x/WIARH4+RkgWvDEFL4jk0t1r7fneVmBVh5ZPoSLGoiHn+mkS/DE2gVIFoFMZ7D1yNOe5232PG9h4EsNkddFZb6IEDVznXO9pP9915DgibBAUi27zBPtsZJI8ERImZYhZW9OeyqQXAvyPK8TD1Oqr7rnt8yDlKCCwDFyJT0SSPLwpDhGhOv6BvbpL+nPkhpXfD/8E85rjNhXnestuy8P/rzyMxsn4uG+J+kSPIqxC5BEsmUtd4K/pDYF1nHeI8zzvJ6BLyPqLkRRIMP/v+8anphHh3OuwPO8IYEEMvUK/JElaWOZ13zfh1eV53cfLWOfkpQd+OORFh6xrdrXN3Bdc2WtRZeIaxzravIzLNnfMJ08z8sO/H5DfKjqek+XNNjzvAJZchbxa1/3PVH/uU3GBA83nj4I/MJ6WvrfE/bGJNUiL5Bh5ibAH9lSue5C2s+NG8InR9IHnuflB1tQAcmszO9bWgIkoLL3WEhcgZaptAROEIFeDIyIluD8+LllFC34IVdSZ7+DSBKdZFnj/rIWVF1pnh81waTOwsDN9xDOfdR0knSDpMae533gdzBJqGJT88Yi0RxOnN/ExvVNfFzj5ML1TnwxdY2TMcETUxcg2QQSDY/QtDQ6nHNPByfZ//P3OfdRs0rlv1sKxPDEERf4jlkS+H/fWVJ+oBUVomeOpCPLvM6ixWZYcX4TG9c38XGNkwvXO/HF1DVOxi5ac2QtSIL4IYuSwB9ZSwJFHoPdtEg2REGg/k5PWQ2GJZz3qFgoG9EkiK6g0dFY5c/zXJHEj5jA93ovSVmB75aFgTpIc8sk1nL3cQjsA+c3sXF9Ex/XOLlwvRNfPFxjzznn5+f7ouLTXKrTR14gwTBXoRZU2c65Rj6GBERcoDVJ48DLTdTgiQ7P8wYHFjeJEUgAAACQJJIywQMAAAAAAJBIkrEGDwAAAAAAQEIhwQMAAAAAABDnSPAAAAAAAADEORI8AAAAAAAAcY4EDwAAAAAAQJwjwQMAAAAAABDnSPAAiHme52V7njfY7zgAAACihfsfANVFggdAPOgpKd/vIAAAAKKI+x8A1UKCB0BM8zyvk6QhkrI9z8vyOx4AAIBI4/4HQE2k+R0AAOyLc26J53mrnHNP+x0LAABANHD/A6AmaMEDIKYFnlpt8jsOAACAaOH+B0BNkOABEOu6SFoQaKoMAACQDLj/AVBtJHgAxLpVkhoHJgAAgGTA/Q+AavOcc37HAAAAAAAAgINACx4AAAAAAIA4R4IHAAAAAAAgzpHgAQAAAAAAiHMkeAAAAAAAAOIcCR4AAAAAAIA4R4IHAAAAAAAgzpHgAQAAAAAAiHP/D/RRfp1VD7wlAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "figheight=8\n", "def Euler_Maruyama(sfunc,z0,time=np.linspace(0,1,100),W=None,args=()):\n", " M,NT=len(z0),len(time)\n", " z=np.zeros((NT,M))\n", " if W is None:\n", " W=np.random.randn(NT,M)\n", " z[0,:]=z0\n", " for i in range(0,NT-1):\n", " tau = time[i+1]-time[i]\n", " drift,diffusion,dummy=sfunc(z[i,:],time[i],*args)\n", " z[i+1,:] = z[i,:] + drift*tau+ diffusion*(W[i+1,:]-W[i,:])\n", " return z\n", "\n", "def Richardson_BWDITO(z,t,epsilon):\n", " drift,diffusion,diffusionprime=(4/3)*2*epsilon*z**(1/3),np.sqrt(2*epsilon)*np.abs(z)**(2/3),(2/3)*np.sqrt(2*epsilon)*np.abs(z)**(-1/3)\n", " return drift,diffusion,diffusionprime\n", "\n", "def Richardson_FWDITO(z,t,epsilon):\n", " drift,diffusion,diffusionprime=(4/3)*2*epsilon*z**(1/3)*0,np.sqrt(2*epsilon)*np.abs(z)**(2/3),(2/3)*np.sqrt(2*epsilon)*np.abs(z)**(-1/3)\n", " return drift,diffusion,diffusionprime\n", "\n", "NT=2**13 #Timesteps\n", "M=500 #MC sampling\n", "\n", "time=np.linspace(0,10,NT)\n", "W,DW=np.zeros((NT,M)),np.random.randn(NT,M)\n", "for i in range(NT-1): W[i+1]=W[i]+np.sqrt(time[i+1]-time[i])*DW[i]\n", "z0=np.random.rand(M)\n", "\n", "fig,ax=subplots(1,2,figsize=(2*figheight,figheight),num='Samples')\n", "\n", "z2= Euler_Maruyama(Richardson_BWDITO,z0,time=time,W=W,args=(1,))\n", "#z3= Euler_Maruyama(Richardson_FWDITO,z0,time=time,W=W,args=(1,))\n", "\n", "\n", "a=ax[0]\n", "i=1\n", "t=time\n", "a.plot(t, z2[:,i],'r',lw=0.4,label='BWD ITO')\n", "#a.plot(t, z3[:,i],'b',lw=0.4,label='FWD ITO')\n", "\n", "#ztrue=z0[i]*np.exp((lam-0.5*mu**2)*t+W[::every,i]*mu)\n", "#a.plot(t, ztrue,'k',lw=1,label='True')\n", "a.legend()\n", "a.set_xlabel('$t$')\n", "a.set_ylabel('$X_t$')\n", "\n", "a=ax[1]\n", "a.plot(t, z2.mean(axis=1)/np.sqrt(z2.var(axis=1)),'r--',lw=2,label='average')\n", "a.plot(t, np.sqrt(z2.var(axis=1)),'r',lw=2,label='variance')\n", "\n", "#a.plot(t, z3.mean(axis=1)/np.sqrt(z3.var(axis=1)),'b--',lw=2,label='average')\n", "#a.plot(t, np.sqrt(z3.var(axis=1)),'b',lw=0.8,label='variance')\n", "\n", "tmp=np.sqrt(z2.var(axis=1))\n", "i=np.argmin(np.abs(t-1))\n", "a.plot(t, 0.7*tmp[i]*t**(3/2),'k',lw=2,label='$\\sim t^{3/2}$')\n", "\n", "\n", "#a.plot(t, ztrue,'k',lw=1,label='True')\n", "a.legend()\n", "a.set_xlabel('$t$')\n", "a.set_ylabel('$X_t$')\n", "a.set_xscale('log')\n", "a.set_yscale('log')\n", "a.set_ylim(1e-1,1e2)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "ExecuteTime": { "end_time": "2019-02-06T03:31:44.510459Z", "start_time": "2019-02-06T03:31:44.498258Z" } }, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "
\n" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%%html\n", "\n", "\n", "\n", "
\n" ] } ], "metadata": { "hide_input": false, "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.8" } }, "nbformat": 4, "nbformat_minor": 2 }