Skip to content
Permalink
master
Switch branches/tags

Name already in use

A tag already exists with the provided branch name. Many Git commands accept both tag and branch names, so creating this branch may cause unexpected behavior. Are you sure you want to create this branch?
Go to file
 
 
Cannot retrieve contributors at this time
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# A26"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"array([ 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130, 140,\n",
" 150])"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"winds = np.arange(20, 160, 10)\n",
"#winds = np.array([20, 50, 100, 150])\n",
"winds"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"p_amb = 100 # Ambient air pressure in kPa\n",
"rho = 1.225 # Air density in kg m^-3"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"array([ 0.49 , 1.1025, 1.96 , 3.0625, 4.41 , 6.0025, 7.84 ,\n",
" 9.9225, 12.25 , 14.8225, 17.64 , 20.7025, 24.01 , 27.5625])"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Eqn. 15.44\n",
"delta_p_max = rho*(winds**2) / 1000. # Convert to kPa\n",
"delta_p_max"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0.5, 1.0, 'Max. tornadic wind speed vs. pressure deficit')"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEcCAYAAADKlrO6AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAA7cUlEQVR4nO3dd5hU5dnH8e9vC713kLI0QVFQioAVLIk1GivYRYPGJJrEvNEYkxhLooklmmisWBE0KvZesCOw9KoILL0vsIALW+73j3MWx3V3mVl2dmZ27891zbVz+j1nzp57zvOc8zwyM5xzztVeaYkOwDnnXGJ5InDOuVrOE4FzztVyngicc66W80TgnHO1nCcC55yr5TwRuL0maZikFRHDcyUNi9O2jpC0cC+WN0k9qjKmqiRpoqTLEh1HTRbLMSDp55LWStomqWX4t9selukczpdeNRHHnyeCSpC0VNIuSa1KjZ8RHmRZcd7+xZI+jec29oaZ9TGziXFa9ydm1ise63YukqRM4C7gR2bWyMw2hn8XV7ScmS0L5ysK15P0yd0TQeUtAUaWDEg6EKifuHCiJykj0TG41KNAQs8Z1XzstgXqAXOrcZsJ4Ymg8p4CLowYvgh4MnIGSSdJmi5pq6Tlkm6MmHaOpMWSmoTDJ0haI6l1RRuVtB/wADA0vPzcHI5vKulJSesl5Ui6oeSfNryC+EzS3ZI2ATdKelzSfZJel5Qn6UtJ3SO2c08Y81ZJ2ZKOiJhWP1w+V9I8YFCpGJdKOjZ8ny7peknfhNvJltSpjM/1hKRrwvf7hFdWV4bDPSRtCk9EpYuhlkr6naRZkrZIelZSvYjp/ydptaRVkkbtYd9eHH4neZKWSDqv1P77d7iNBZKOiViuqaRHw+2slHRLZLGApFGS5of7621JXSKmHReub4uk/wAqJ7YOkr6V1CJi3MGSNkjKDPfRR+F6Nkh6tqLPWuozV/TZJkq6VdJnwA6gm6Tekt4Nv5OFks6OmP9ESfPCfbhS0u/C8a0kvSZpc7jcJxHH5/eKasJj65bw/TBJKyRdK2kN8JikNEnXhcfURknPRe6XMj5juceApLqS7pC0TEER0APh8b0vUFIEuVnSB6VjDee7U8H/2xZJn4bjssL5MiTdChwB/EfB/+t/ovleqp2Z+SvGF7AUOJbgQNkPSAeWA10AA7LC+YYBBxIk3L7AWuC0iPWMBR4HWgKrgJOj3P7FwKelxj0JvAw0BrKAr4BLI+YvBH4FZBBcuTwObAIOCceNBcZHrO/8MK4M4BpgDVAvnHYb8AnQAugEzAFWlN4/4fv/A2YDvQhOcv2AlmV8plHAq+H7c4FvgGcjpr0csU9Lb2sy0CGMZz5wRTjt+HCfHwA0BJ4Jv58eZWy/IbAV6BUOtwf6lNp/vwEygXOALUCLcPpLwIPhOtqE8VweTjsNWERwnGQANwCfh9Nahds8M1zvb8LtXFbO9/4B8LOI4X8CD4TvxwF/JDjW6gGHx3AsVfTZJgLLgD5h/E0JjvVLwuH+wIaIfbUaOCJ83xzoH77/O8EPmMzwdQSgcNr3vhOCY/OWiO+7ELgdqEtw7P4amAR0DMc9CIwr5/NVeAwA/wJeITh2GgOvAn8Pp2WF82ZErC9y2fvC/bMPwTng0DCe7y0XzlPmd5osr4QHkIovvksEN4QH+PHAu+E/xu5EUMZy/wLujhhuFv6TzQYejGH7FxORCMKDcCewf8S4y4GJEfMvK7WOx4FHIoZPBBZUsM1coF/4fjFwfMS00ZSfCBYCp0bxmboDmwlOZA+E8a8Ipz0B/DZ8P6yMbZ0fMfwPvjs5jgFui5i2LxUngs3AGUD9Mvb3KsITVzhuMnABQfHBzshlCIoMPwzfv0mYkMPhNIJf1l0IrignRUwTsKK8kwZwGfBBxLzLgSPD4SeBh4COMR7L5X628P1E4KaIaecAn5Rax4PAX8L3y8LvrkmpeW4i+KFS1r7fUyLYRfgjJBw3HzgmYrg9UEDECTtiWrnHQLgPtwPdI6YPBZaE77MoJxGE3+O3hP8Tpbb5veVIgUTgRUN75ymCX68XU6pYCEDSYEkfKiiu2QJcQfArEAAz2wz8j+DXyp17EUcroA6QEzEuh+CXSonlZSy3JuL9DqBRROzXhMUZWxQUPzWNiL1DqfVFbre0TgS/7itkZt8A24CDCH4tvgasktQLOAr4qILFy/scUcdpZtsJTnJXAKsVFJn1jphlpYX/1RHr6kBwQs8Ml9kc7qsHCa4MCKffEzFtE8EJaJ/S8YXrL+t7KvE8QZFgB+BIgpPNJ+G034frnazgrq0Ki8FKKe+zlYiMqQswuOTzhJ/pPKBdOP0Mgh8VOWFR1dBw/D8JrozeUVD8dl0M8a03s/xSMUyI2P58oIggKZdW0THQGmgAZEes661w/J60Irjy2uOxnQo8EewFM8shqDQ+EXixjFmeIbjs7GRmTQl+6e4uA5Z0EEGxxzjg3lg2XWp4A8Evoi4R4zoDKytYplwK6gOuBc4GmptZM4LigpLYVxOc4CO3VZ7lBL/2o/ERQTFJHTNbGQ5fSFDEMCPKdUSKJU7M7G0zO47gF+YC4OGIyftIiiy/70zwS3o5wRVBKzNrFr6amFmfcL7lBMVEzSJe9c3s89Lxhev/Qf1JRHybgXcIvpdzCYpDLJy2xsx+ZmYdCH6R36/ob5Mt77Pt3nTE++XAR6U+TyMz+3kYxxQzO5UgEb4EPBeOzzOza8ysG3AK8NuIuogdBCfkEu34vtLH7nLghFIx1AuPmdIqOgY2EPyq7xOxnqZm1og92wDkE92xHfX/XqJ4Ith7lwJHh78oS2sMbDKzfEmHEPzzAqCgQvNp4HqC8tZ9FFaORmEt0FFSHQALblN7DrhVUmMFlZG/DddfGY0JymXXAxmS/gw0iZj+HPAHSc0ldSSoeyjPI8DNknoq0FdSy3Lm/Qj4JfBxODwxXPen4WeM1XPAxZL2l9QA+Et5M0pqK+knkhoSnNi3EfzKLNEGuEpBxexZBGX+b5jZaoKT852SmoQVmd0lHRUu9wDBvuoTbqdpuDzA60AfSacruBvmKn54EiztGYLkeEb4viT+s8LvAoJiPCsVf0XK/GzlzPsasK+kC8L5MyUNkrSfpDqSzpPU1MwKCOo/Sm6hPFlBhbYixpfENwM4V8GNBccTXAFW5AGCY71LuO7Wkk4tZ95yjwEzKyZI9ndLahOuax9JP97D9kuWHQPcpaAiP13SUEl1y5h9LVDhsweJ5olgL5nZN2Y2tZzJVwI3ScoD/kz46yj0d4Ky7v+a2U6CytlbJPWE3Q9lnVfOej8guKVtjaQN4bhfEZR3LgY+JThJjKnkx3qboGz7K4JL6Xy+f3n913D8EoKT4FMVrOsugs/9DsEJ4FHKv832I4IkVJIIPiX4pfhxOfNXyMzeJKiX+YCgWOKDCmZPI6gUX0VQfHMUwfdX4kugJ8EvwVuBM81sYzjtQoKiuXkEJ+HnCa4qMLMJBBWd4yVtJahYPyGctgE4i6DyfWO4/s/28LFeCedba2YzI8YPAr6UtC2c52ozWwJ7PJb29Nm+x8zygB8BIwj21Rq+q8iFoN5kafhZryA4rgnX/x5Bgv0CuN++e9bkaoKrhM0ExUwv7WEf3BN+xnfC/61JwOBy4t3TMXBtOH5SGPN7BDc2RON3BPV7UwiOmdsp+5x6D3CmgrvGYrnyrzYltfbOuXJIupigsu/wRMdS1WryZ3PR8ysC55yr5TwROOdcLedFQ845V8v5FYFzztVyngicc66WS8lWKFu1amVZWVmVWnb79u00bNiwagOqJh57YqRq7KkaN3js8ZKdnb3BzH7w5HRKJoKsrCymTi3v1v2KTZw4kWHDhlVtQNXEY0+MVI09VeMGjz1eJJXZzIoXDTnnXC3nicA552o5TwTOOVfLeSJwzrlazhOBc87Vcp4InHMuRWTn5HLfh4vIzsmt0vWm5O2jzjlX22Qv3cSIhydRVGzUyUhj7GVDGNCleZWs268InHMuBdz+9gIKioxig4LCYiYtLrPLiErxKwLnnEtyj366hMlLcklPE5iRmZHGkG7ldfQXO08EzjmXxF6esZKbX5vH8X3acekRXZm8ZBNDurWssmIh8ETgnHNJ65Ov1/O7/83kkK4t+NeIg6iXmc6grBZVvh2vI3DOuSQ0e8UWrngqm+6tG/HwhQOpl5ket215InDOuSSTs3E7lzw+mWYN6vDEqENoWj8zrtvzROCcc0lkfd5OLnh0MkXFxpOXHkLbJvXivk2vI3DOuSSxbWchlzw+mfV5O3nmZ4Pp3rpRtWzXE4FzziWBnYVFXP7UVOavzuORiwZycOequytoT7xoyDnnEqy42LjmuZl8tmgjt5/Rl+G92lTr9j0ROOdcApkZN78+j9dmreba43tz5oCO1R6DJwLnnEugBz5azGOfLeWSw7K44qhuCYnBE4FzziXI89kruP2tBZzSrwN/Oml/JCUkDk8EzjmXAB8uWMe1L8zisB4tueOsvqSlJSYJgCcC55yrdtOX5XLl2Gn0bteYB84fQN2M+D01HA1PBM45V42+Wb+NUY9PoXXjujx+ySE0rhffp4ajUW2JQFInSR9Kmi9prqSrw/E3SlopaUb4OrG6YnLOueq0dms+Fz46mfQ08eSoQ2jduG6iQwKq94GyQuAaM5smqTGQLendcNrdZnZHNcbinHPVasu3BVw0ZjKbd+xi/OihZLVqmOiQdos5EUhqCOSbWVEsy5nZamB1+D5P0nxgn1i375xzqSa/oIifPTmVb9ZvY8zFgziwY9NEh/Q9MrOKZ5DSgBHAecAgYCdQF1gPvAE8ZGZfx7RRKQv4GDgA+C1wMbAVmEpw1fCDnpkljQZGA7Rt23bA+PHjY9nkbtu2baNRo+ppv6OqeeyJkaqxp2rcULNiLzbjvhk7yV5bxBV96zKkQ+Ja9hk+fHi2mQ38wQQzq/AFfAT8CegLpEWMbwGcAbwAnL+n9UQs1wjIBk4Ph9sC6QT1FbcCY/a0jgEDBlhlffjhh5VeNtE89sRI1dhTNW6zmhN7cXGx/XHCLOty7Wv2yCeLExdUCJhqZZxTo0lNx5pZQRkJZFOYBF6QFFW1dzjfC8BYM3sxXM/aiOkPA69Fsy7nnEt2//5gEU9PWsblR3Xj0sO7Jjqccu0xEUQmAUnNgZ5AvYjpH5eVKEpT8Mjco8B8M7srYnx7C+oPAH4KzIk+fOecS07jJi/jrne/4vT++3Dd8b0THU6Foi6sknQZcDXQEZgBDAG+AI6OchWHARcAsyXNCMddD4yUdBBgwFLg8mhjcs65ZPTO3DX8ccJshvVqze1n9E1Y0xHRiqXW4mqCyuJJZjZcUm/gr9EubGafAmXtjTdiiME555JWdk4uj8/ZyWerp3Fgx2bcf15/MtOT/7ndWBJBvpnlS0JSXTNbIKlX3CJzzrkUkp2Ty7kPT2JnYTECfn1sTxrUSY2+v2JJVSskNQNeAt6V9DKwKh5BOedcqnlzzmp2FhYDIMG8VVsTHFH0ok5XZvbT8O2Nkj4EmgJvxSUq55xLIcs37WDCtJVAUP5dJyONId1aJjaoGOwxEUiqB1wB9ABmA4+a2UfxDsw551LBqs3fMvLhSRQWG3ee1Y8vZsxj5LGDGNCl+voc3lvRXBE8ARQAnwAnAPsTVBw751yttm5rPuc+PIktOwoY+7PB9O3YjJZ5i1IqCUB0iWB/MzsQQNKjwOT4huScc8lvw7adnPvIl6zL28lTlwZJIFVFU1m8+2ExMyuMYyzOOZcScrfv4vxHvmRF7g4euzi1ioHKEs0VQT9JW/nuGYD6EcNmZk3iFp1zziWZLd8WcMGYL1m8YTtjLhrE4BSqFC5PNE1MJLYPNeecSxLbdhZy0ZjJLFyTx0MXDOTwnq0SHVKViPo5AkkDyhh3StWG45xzyWnHrkJGPTaF2Su38J9z+zO8d5tEh1RlYnmg7GFJB5YMSBoJ3FD1ITnnXHLJLyjisiemMjVnE/865yB+3KddokOqUrE8/3wm8Lyk84DDgQuBH8UlKuecSxI7C4u44ulsvli8kTvP6scp/TokOqQqF8uTxYsljSBoYmI58CMz+zZegTnnXKIVFBXzy2emM3Hhem47/UBO798x0SHFRTRPFs8maCK6RAuCHsW+lISZ9Y1XcM45lyiFRcX8evwM3p23lptO7cOIQzonOqS4ieaK4OS4R+Gcc0mkqNj43f9m8vrs1dxw0n5cODQr0SHFVTSJYFnY12W5JGlP8zjnXCooLjauf3E2L81Yxf/9uBeXHdEt0SHFXTR3DX0o6VeSvnddJKmOpKMlPQFcFJ/wnHOu+pgZf3llLs9OXc5VR/fgF8N7JDqkahHNFcHxwChgnKSuwGaCPovTgXeAu81sRrwCdM656mBm3PL6fJ6alMPlR3bjN8ftm+iQqk00TxbnA/cD90vKBFoB35rZ5jjH5pxz1cLM+OfbC3n00yVcfGgW153QO+n7Ga5KMfWjZmYFwOo4xeKccwlx7/uLuH/iN5w7uDN/OWX/WpUEILYni51zrsb578RvuPu9rzhzQEduOfWAWpcEwBOBc64WG/PpEm5/awGn9OvA7Wf0JS2t9iUBqGQikOQtkjrnUtrTk3K46bV5HN+nHXed3Y/0WpoEoPJXBA9JagAg6cgqjMc55+LuuanLueGlORzTuw33jjyYzPTaXTgSU2VxhD8Dj0oqBGYAH1dZRM45F0f/eu8r/vXe1/Tr2JT7zutPnYzanQSg8lcENwMLCdogeq7qwnHOufi5652F/Ou9rwFYuCaPuau2Jjii5FDZK4Lfm9kGSQ2Be4DLqjAm55yrck9PyuHeDxbtHi4oKmbS4o0p399wVajsFcEvJb1JkASmR7OApE6SPpQ0X9JcSVeH41tIelfS1+Ff/1acc1XqkU8Wc8NLcxjQpTn1MtJIF2RmpDGkBvQ3XBUqe0XQDJgEPAn8JsplCoFrzGyapMZAtqR3gYuB983sNknXAdcB11YyLuec283MuPf9Rdz93lecdGB77j7nIGav3MKkxRsZ0q2lXw2EKpsINoXLrgvf75GZrSZ8KtnM8iTNB/YBTgWGhbM9AUzEE4Fzbi+ZGbe9tYAHP1rMGf07cvsZB5KRnsaALs09AZRSqURgZjdJ6gDcC8yJdXlJWcDBwJdA2zBJYGarJdWcHqGdcwlRXGzc+Opcnvwih/OHdOamnxxQax8Wi4Yq042ApP2BLsBsM1sR47KNgI+AW83sRUmbzaxZxPRcM/tBupY0GhgN0LZt2wHjx4+POW6Abdu20ahRo0otm2gee2KkauypGjfsXezFZoyZs4tPVxZyfFYm5/TKrNZmI5J5vw8fPjzbzAb+YIKZxfwC/gecDdwEPBHDcpnA28BvI8YtBNqH79sDC/e0ngEDBlhlffjhh5VeNtE89sRI1dhTNW6zyse+q7DIrhybbV2ufc3ufnehFRcXV21gUUjm/Q5MtTLOqZWtI3jXzJ4jhmcIFKTkR4H5ZnZXxKRXCDq2uS38+3IlY3LO1WL5BUX88plpvDd/Hdef2JvRR3ZPdEgpo7KJ4FBJxwMb+eGJvTyHARcAsyXNCMddT5AAnpN0KbAMOKuSMTnnaqkduwoZ/WQ2ny7awM2nHcAFQ7okOqSUUtlEMMfM7pCUAfSJZgEz+xQor6DumErG4Zyr5bbmFzDqsSlMW5bLHWf148wBHRMdUsqpbCI4WdJO4G0zm1mVATnnXLRyt+/ioscmM2/VVv49sj8n9W2f6JBSUmWfLD4H+Bo4XdLDVRiPc85FZV1ePiMemsSCNXk8eMEATwJ7YY9XBJIuAu4kSBqvAb8ws7XAW+HLOeeq1arN33L+I1+yZms+j108iMN6tEp0SCktmiuCPwHHAb2BHOBvcY3IOecqkLNxO2c98AXr83by1KWHeBKoAtHUEWw1s5KG5f4k6ct4BuScc+VZtC6P8x75kl2FxYwbPYQD9mma6JBqhGgSQfvwqd75wAKCh8Kcc65azV21hQsfnUxamhg/eii92jVOdEg1RjSJ4C9AX+A84ECgkaQ3gJnALDMbF8f4nHOO6ctyuWjMZBrVzWDsz4bQtVXDRIdUo+wxEZjZQ5HDkjoSJIb+wPmAJwLnXNxMWryRSx+fQqvGdRl72WA6Nm+Q6JBqnJieI5B0EDASGAGsIahAds65uJi4cB2XP5VN5xYNGHvZYNo0qZfokGqkaG4f3ZfgxH8usJ2gfaGjzGyppCVxjs85V0u9NWcNvxo3jX3bNubJUYfQslHdRIdUY0VzRbAAmAKcaWal+x6IvQ1r55yrQHZOLvfPyGfq2mwO6tSMxy45hKb1/R6VeIrmOYIzgKXAu5KeknSKJP9WnHNVLjsnl3Me/ILJa4oAuOZHvTwJVIM9JgIzm2Bm5wA9CJ4kvhxYIekxoEmc43PO1RLFxcbtby2gsDgoaBAwY/nmhMZUW0Td1pCZbTezsWZ2MrAfQef1s+MWmXOu1thZWMSvn53B5CWbSE8TaUBmRhpDurVMdGi1QmX7LN4EPBi+nHOu0rbsKGD0U1P5cskmrj2+N4dkNWfc+1MZeewg72S+mlS2GWrnnNtrK3J3cMljU1i6cTv3jDiIUw/aB4C87nU8CVQjTwTOuYSYs3ILox6fwrcFRTw5ajBDu3sxUKJ4InDOVbuPvlrPlU9n07R+Ji/8/FD2bevtBiVS1JXFkm6PZpxzzlXkuSnLGfX4FLq0bMiEXxzmSSAJxNJD2XFljDuhqgJxztVsZsZd737F71+YxWE9WvHcFUNp601GJIVompj4OXAl0F3SLL7rgL4x8FkcY3PO1RAFRcX84cXZPJ+9grMHduTWnx5IZnple8p1VS2aOoKxwJsEPZNdR5AIDMgzs9w4xuacqwHy8gu4cuw0Pvl6A78+tidXH9MTSXte0FWbaBLBG2Z2uKSfACdHjJckMzN/utg5V6Y1W/K5+LHJLFq3jX+e2ZezBnZKdEiuDNH0R3B4+LdR/MNxztUUC9fkcfFjk8nLL2TMxYM4ct/WiQ7JlcNvH3XOVbnPF23g8qezqZ+ZzrOXD6FPB+9bOJnFcvvoE5KaRQw3lzQmLlE551LWhOkruOixybRvWo8JvzjMk0AKiOWKoK+ZbS4ZMLNcSQdXfUjOuVRkZtw/8Rv++fZChnZryQMXDPAmpFNELIkgTVLzkjuFJLWIcXnnXA1VWFTMn1+ZyzNfLuO0gzpw+5l9qZuRnuiwXJRiOZHfCXwu6XmC20fPBm6NduGwGOlkYJ2ZHRCOuxH4GbA+nO16M3sjhpiccwm2fWchvxo3nQ8WrOPKYd353Y96kZbmt4emkqgTgZk9KSkbGE7wLMHpZjYvhm09DvwHeLLU+LvN7I4Y1uOcSxLr8vK59PGpzF21hVt/egDnDe6S6JBcJcRUtGNmc4G5ldmQmX0sKasyyzrnkkt2Ti6vz17NqzNXsi2/iIcvHMgx+7VNdFiukmRWcf/zkj4NHyjL4/ud1QuI6YGyMBG8Vqpo6GJgKzAVuKa8p5UljQZGA7Rt23bA+PHjo93s92zbto1GjVLzkQiPPTFSNfZ4xb0ot4jbpuRTWBwMX3JAHY7qWLWVwqm6zyG5Yx8+fHi2mQ0sPT6aK4Il4d8/mdm/qjQq+C9wM0GCuZmgHmJUWTOa2UPAQwADBw60YcOGVWqDEydOpLLLJprHnhipGnu84h731FQKi/MBSBO02qcrw4b1qNJtpOo+h9SMPZrnCPpL6gJcEj470CLytTcbN7O1ZlZkZsXAw8Ahe7M+51z87Cos5voJs3l77lrSBOmCOt6vcI0QzRXBg8BbQDdgWqlpFo6vFEntzWx1OPhTYE5l1+Wci591eflc+fQ0pubkcsVR3TlmvzZMXrKJId1aepeSNUA0bQ3dC9wr6b9m9vPKbkjSOGAY0ErSCuAvwDBJBxEklKXA5ZVdv3MuPqYvy+WKp7PZ+m0h/x55MKf06wDAoKy9KhBwSSSWu4Z+Iel8oKuZ3SypM9DOzCZHs7CZjSxj9KMxbN85V82em7KcG16aQ5smdXnh54eyfwdvbLgmiiUR3AcUA0cTVOzmAS8Ag+IQl3MugQqKirn5tXk8+UUOh/doxb9HHkzzhnUSHZaLk1gSwWAz6y9pOuxua8iPDOdqmPV5O/nF2GlMXrqJ0Ud24/c/7kWG9yZWo8WSCAokpRM+SyCpNcEVgnOuhpi5fDNXPJ1N7o5d3DPiIE49aJ9Eh+SqQSyJ4F5gAtBG0q3AmcANcYnKOVftns9ewfUTZtO6UV2ev+JQDtjHm4+uLWJpa2hs2NbQMQRPFZ9mZvPjFplzrloUFBVz6+vzefzzpQzt1pL7zutPC68PqFVibWtoAbAgTrE456rZhm1BfcCXSzYx6rCuXH9ib68PqIX2mAgk/bai6WZ2V9WF45yrLrNXbOHyp6aycfsu7jq7H6f375jokFyCRHNF0Dj824vgVtFXwuFTgI/jEZRzLr5enLaCP7w4m5YN6/D8FYdyYEevD6jNonmy+K8Akt4B+ptZXjh8I/C/uEbnnKtShUXF/O2NBYz5bAmDu7bgvvP606pR3USH5RIsljqCzsCuiOFdQFaVRuOci5tN23fxy2em8fk3G7n40Cz+eNJ+ZHp9gCO2RPAUMFnSBIJnCX4KPBGXqJxzVWrOyi1c/lQ267ft5J9n9uWsgZ0SHZJLIrHcPnqrpDeBI8JRl5jZ9PiE5ZyrKi/PWMm1L8yiWf06/O/yofTr1CzRIbkkE+vto9P4YVPUzrkkk52TyyuLdvHi6mm8MnM1g7Kac/95A2jd2OsD3A/FlAicc8kvOyeXcx+exM7CYmA1x/dpx70jD6ZOhtcHuLL5keFcDfPslGVhEgiaADiwY1NPAq5CUR8dCpwv6c/hcGdJ3rWkc0mioKiY299awHNTVyCCJFA307uSdHsWS9HQ/XzXH8FNeH8EziWNZRt38Kvx05m5fDMjBnXiJ/068OJH0xh57CDvStLtkfdH4FyKe3nGSv44YQ4S/Ofcgzm5b9CV5K4VdTwJuKh4fwTOpajtOwv588tzeWHaCvp3bsY9Iw6mU4sGiQ7LpSDvj8C5FDR7xRauGj+dpRu3c9XRPbjqmJ7eaqirtKgSgSQRNDDn/RE4l0DFxcaYz5Zw+1sLaNmwLs9cNoSh3b0y2O2dqBKBmZmkl8xsAN4fgXMJsT5vJ7/730w++mo9x+3fln+c0dc7lHdVIpaioUmSBpnZlLhF45wr08dfree3z81ka34BN5/ah/OHdCG4UHdu78WSCIYDV0haCmwnKB4yM+sbj8Ccc7CrsJg73lnIQx8vZt+2jXj6skPo3a5JosNyNUwsieCEuEXhnPuBJRu2c/X46cxasYXzBnfmhpP2p36d9ESH5WqgWBLBReWMv6kqAnHOfefFaSv400tzyEhP44HzB3D8Ae0SHZKrwWJJBNsj3tcDTgb8riHnqlBefgF/fnkuE6av5JCsFvxrxEF0aFY/0WG5Gi6W/gjujByWdAff9V+8R5LGECSPdWZ2QDiuBfAsQU9nS4GzzSw32nU6V5PMXL6Zq8ZPZ/mmHfzm2H355dE9SE/zCmEXf3vzBEoDoFsM8z8OHF9q3HXA+2bWE3g/HHauVikuNh746BvO+O/nFBYZz14+lKuP7elJwFWbqK8IJM0mbF4CSAdaE0P9gJl9LCmr1OhTgWHh+yeAicC10a7TuVSWnZPL+/PX8tmiDcxcsYUTDmjHbaf3pWmDzESH5mqZWOoITo54XwisNbPCvdx+WzNbDWBmqyW12cv1OZcSsnNyGfHQFxQUBb+trjiqG9ce39ufDXAJITPb81xVtbHgiuC1iDqCzWbWLGJ6rpmV2VyipNHAaIC2bdsOGD9+fKVi2LZtG40aNarUsonmsSdGVceet8v455RvWZYX/O8JOKNnJid3r9qnhH2fJ0Yyxz58+PBsMxtYenwsRUNnAW+ZWZ6kG4D+wC1hP8aVtVZS+/BqoD2wrrwZzewh4CGAgQMH2rBhwyq1wYkTJ1LZZRPNY0+MqordzHh11mpufGUuW741MtKEmZGZkRaXfgN8nydGKsYeS9HQn8zsf5IOB34M3AH8Fxi8F9t/heD5hNvCvy/vxbqcS1prtuRzw0uzeW/+Ovp1asa4M4awbWchkxZvZEi3lt5vgEuoWBJBUfj3JOC/ZvaypBujXVjSOIKK4VaSVgB/IUgAz0m6FFgGnBVDPM4lPTNj/JTl/O31+RQUF3PDSftxyWFdd98R5AnAJYNYEsFKSQ8CxwK3S6pLDLefmtnIciYdE0MMzqWMnI3bue6F2XyxeCNDu7XktjMOpEvLhokOy7kfiCURnE3wHMAdZrZZUjvg/+ITlnOpq6jYeOyzJdzxzkIy09L4++kHMmJQJ78jyCWtWBLBSZRRWRyfsJxLTQvX5PH7F2Yxc/lmjt2vDbecdiDtmtZLdFjOVSjRlcXO1Qi7Cou578NF3D9xEY3rZXLvyIM5pW97vwpwKaHaKoudq6lmLN/M75+fyVdrt3HaQR348yl9aOE9h7kUUpnK4uOoRGWxczXNt7uKuPOdhYz5bAltm9RjzMUDObp320SH5VzM9qayuD1eWexqqc+/2cB1L8xm2aYdnDe4M9ed0JvG9byNIJeaYkkE3wINgZEEjc1lApvjEJNzSWtrfgF/f2M+4yYvJ6tlA8aPHsKQbi0THZZzeyWWRHA/UAwcTZAI8oAXgEFxiMu5pPPuvLXc8NJs1uft5PIju/HrY/f1riNdjRBLIhhsZv0lTQcws1xJXiPmarTsnFye/2oXD389ic++2Ujvdo15+MKB9O3YLNGhOVdlYkkEBZLSCfskkNSa4ArBuRope+kmRjw8KWwqeiMjB3Xir6ceQJ0Mv0fC1SyxHNH3AhOANpJuBT4F/haXqJxLsNkrtnDV+Om7+wtIE3Rs0cCTgKuRoroiUPBUzMdANkHbQAJOMzPvvN7VKOu25vPPtxfy/LQVNK6bQWa6KCoy6mSkeaWwq7GiSgRmZpJeMrMBwII4x+RctcsvKOLRT5dw/4eL2FVUzOgjuvGLo3vw9dptjHtvSlz6C3AuWcRSRzBJ0iAzmxK3aJyrZmbGW3PW8Lc357N807cct39b/njifmS1CloJHdClOXnd63gScDVaLIlgOHCFpKXAdoLiITOzvvEIzLl4m7NyCze/No8vl2yiV9vGjL1sMIf1aJXosJyrdrEkghPiFoVz1Wh93k7ufGchz05dTvMGdbjltAMYMagTGeleEexqp1gSwVrgSuBwgltIPyVofdS5lLCzsIjHPlvKfz5YRH5BEZce1pVfHdOTpvW9aQhXu8WSCJ4keJr43+HwSOApvHtJl+TMjHfmreVvb8wnZ+MOjt2vDdefuB/dWjdKdGjOJYVYEkEvM+sXMfyhpJlVHZBzVWn+6q3c9Oo8vli8kZ5tGvHkqEM4ct/WiQ7LuaQSSyKYLmmImU0CkDQY+Cw+YTm3dzZs28md73zFs1OW0aR+Jjed2odzD+ns9QDOlSGmtoaACyUtC4c7A/MlzcbvHnJJYldhMU98vpR73/+abwuKuOjQLK4+pifNGnizWM6VJ5ZEcHzconBuL5kZ781fx62vz2Ppxh0M69WaG07anx5tvB7AuT2JOhGYWU48A3GuMrJzcnlt5iqyc3KZtXIL3Vs35LFLBjG8V5tEh+ZcyojlisC5pPLGrNX8atx0iixoGG7UYVn84cT9yPR6AOdi4onApZylG7bznw8X8cK0FYQ5gHRBy0Z1PQk4VwmeCFzKWLJhO//5YBEvzVhJRpo48YB2vDd/HYVFxWR666DOVdpeJQJJ7cxsTVUF41xZFq/ftjsB1MlI4+JDs7j8yG60aVKP7JxcJi3eyJBuLb1hOOcqaW+vCB4FTqqKQJwrbdG6bfzng695ZeYq6mSkMeqwrow+qhttGtfbPc+ALs09ATi3l/YqEZhZlSSBsEXTPKAIKDSzgVWxXpeaFq3L498fLOKVmauol5HOZUd042dHdKN147qJDs25GinqRCBpfzObV2rcMDObWEWxDDezDVW0LpeCvl6bx70fLOK1Wauon5nO6CODBNCqkScA5+IpliuC5yQ9BfwDqBf+HQgMjUdgrvZYuCaPez/4mjdmr6ZBZjpXHNWdyw7vSktPAM5VC1nJ/Xd7mlFqCNwODAAaA2OB282seK+DkJYAuQTNWz9oZg+VMc9oYDRA27ZtB4wfP75S29q2bRuNGqXm06Y1LfblecW88s0upqwpol46HNslkx9nZdK4jhIUZdlSdb+natzgscfL8OHDs8sqeo/liqAA+BaoT3BFsKQqkkDoMDNbJakN8K6kBWb2ceQMYXJ4CGDgwIE2bNiwSm1o4sSJVHbZRKspsc9fvZV73/+aN+esoVHdDH45vAeXHt6V5g2Tsz2gVN3vqRo3eOzVLZZEMAV4GRgEtAQelHSmmZ25t0GY2arw7zpJE4BDgI8rXsqliuycXF77ZhcbGi/n3XlreXvuWhrXzeCqo3sw6vCu3iCccwkWSyK41Mymhu/XAKdKumBvAwiLnNLMLC98/yPgpr1dr0sO2Tm5jHx4ErsKi3n+61k0qJPOVcf05NLDutK0gfcM5lwyiKXRuamSmgM9CYqGAKqiIbq2wARJJfE8Y2ZvVcF6XQLtKizmzTmruf2tBewqDEoQBVx2eFd+e9y+iQ3OOfc9sdw+ehlwNdARmAEMAb4Ajt6bAMxsMdBvjzO6lLBmSz7PfJnDM5OXs2HbTto1qUtGmiguNupkpnGUtwrqXNKJpWjoaoL6gUlmNlxSb+Cv8QnLpRIz48slm3jyi6W8PXctxWYc3asNFwztwpE9WzN9+WbGvTeFkccO8qeAnUtCsSSCfDPLl4Skuma2QFKvuEXmkt72nYVMmL6Sp77IYeHaPJrWz+Syw7ty3uAudG7ZYPd8A7o0J697HU8CziWpWBLBCknNgJcIbvHMBVbFIyiX3L5Zv42nvsjhhewV5O0spE+HJvzjjL6c0q8D9eukJzo851yMYqks/mn49kZJHwJNAK/UrSWKio0PFqzjyS+W8snXG8hMFycd2J4LhmbRv3Mzwsp+51wK2mMikPRKeZOAnwE/qdKIXFLZtH0Xz05ZztOTcli5+VvaNanHNcfty4hDOnsjcM7VENFcEQwFlgPjgC8JEoCr4Wat2MwTn+fw6qxV7CosZmi3ltxw0n4ct39bMrwXMOdqlGgSQTvgOGAkcC7wOjDOzObGMzBXfUo6d+nfuRmrt+TzxBc5zFy+mQZ10jlnYCcuGNqFfds2TnSYzrk42WMiMLMigrqAtyTVJUgIEyXdZGb/jneALr6yc3I5N3zyt6T5wW6tG3LjKftz+oCONKnnT/86V9NFVVkcJoCTCJJAFnAv8GL8wnLxtnzTDl6btZrHPlvCzsLv2g48e2BHbj+jr1f+OleLRFNZ/ARwAPAm8FczmxP3qFxcrN2az+uzVvPqrFVMX7YZgJ5tGrFp+y7MjMyMNM4Z1NmTgHO1TDRXBBcA24F9gasiThICzMyaxCk2VwU2bd/Fm3NW8+rMVXy5ZBNmsF/7Jvz++F6c0rcDnVo08A7gnavloqkj8FtEUszW/ALembuWV2eu4tNFGygqNrq1bshVR/fklH7t6dHm+xW/3gG8c7XbXnVe75LHjl2FvD9/Ha/OXMXEhevZVVTMPs3q87MjunFKv/bs376JF/k458rkiSCF7Sws4qOF63l11mrem7eWbwuKaNO4LucN6cwp/TpwcCd/4tc5t2eeCFJESS9fDbps5NuCYl6duYq3564hL7+Q5g0y+Wn/fTilbwcO6dqC9DQ/+TvnoueJIAV88c0GLhozhV1FxTz/9SQAGtfN4Ed92nFKv/Yc1qMVmf60r3OukjwRJKHiYmPe6q188vUGPlu0gS8Wb6So2HZPP+2gDtx2Rl/qZXpLn865veeJIEks37SDTxdt4NNFG/h80QZydxQA0KttY07o04535q2hsCjo5euCoVmeBJxzVcYTQYJs2VHA599s2H3yz9m4A4C2TeoyvHcbjujZisO6t6JNk6B76OycXO/lyzkXF54IqsnOwiKyc3L5bNEGPv16A7NXbqHYoGGddIZ2b8nFh2ZxRM9WdG/dqMw7fbyXL+dcvHgiiJPiYmPBmjw+W7SBTxZtYPKSjeQXFJOeJg7u1IxfHd2TI3q2ol+nZl7R65xLKE8EVcTMeGfeWiZMX8m2/EIWrNnKhm27AOjRphEjBnXm8B6tGNytBY29RU/nXBLxRFAJxcXGko3bmbtqK3NXbWHeqq3MXL6ZrfmFu+c5ct9W/KTfPhzeoxXtmtZLYLTOOVcxTwR7sLOwiK/WbGPe6i3hiX8r81dvZceuIgAy08W+bRvTpWUD5qzcigHpgsFdW3LmgI6JDd4556LgiSBCXn4B88KTfcmv/UXrtlEY3sPfqG4G+7dvwtkDO7F/hyb06dCEnm0aUycjjeycXM57ZBIFhcVkZqQxpFvLBH8a55yLTq1KBCXNNDTumkun5vV3n+znrtrKvNVbd9/CCdCqUV36dGjC0b3b0KdDU/p0aELnFg1IK6f5hgFdmjP2siHenLNzLuXUmkSQnZPLiIe+oKDIeP7rz783rUvLBvTp8P1f+m0ax16u7805O+dSUVIkAknHA/cA6cAjZnZbVW9j0uKNFBZ910zD0b3bcPmR3divQxPvl9c5V6sl/AZ2SenAfcAJwP7ASEn7V/V2hnRrSd3MNNKAeplp/GJ4DwZ3a+lJwDlX6yU8EQCHAIvMbLGZ7QLGA6dW9UZKyvBP75nJ2MuGeBGOc86FkqFoaB9gecTwCmBwPDbkzTQ459wPycz2PFc8A5DOAn5sZpeFwxcAh5jZr0rNNxoYDdC2bdsB48ePr9T2tm3bRqNGjfYu6ATx2BMjVWNP1bjBY4+X4cOHZ5vZwNLjk+GKYAXQKWK4I7Cq9Exm9hDwEMDAgQNt2LBhldrYxIkTqeyyieaxJ0aqxp6qcYPHXt2SoY5gCtBTUldJdYARwCsJjsk552qNhF8RmFmhpF8CbxPcPjrGzOYmOCznnKs1Ep4IAMzsDeCNRMfhnHO1UTIUDTnnnEughN81VBmS1gM5lVy8FbChCsOpTh57YqRq7KkaN3js8dLFzFqXHpmSiWBvSJpa1u1TqcBjT4xUjT1V4waPvbp50ZBzztVyngicc66Wq42J4KFEB7AXPPbESNXYUzVu8NirVa2rI3DOOfd9tfGKwDnnXARPBM45V8vV6EQgqZOkDyXNlzRX0tXh+BaS3pX0dfg3KdullpQuabqk18LhVIm7maTnJS0I9/3QFIr9N+GxMkfSOEn1kjV2SWMkrZM0J2JcubFK+oOkRZIWSvpxYqLeHUtZsf8zPGZmSZogqVnEtKSOPWLa7ySZpFYR45Im9vLU6EQAFALXmNl+wBDgF2HvZ9cB75tZT+D9cDgZXQ3MjxhOlbjvAd4ys95AP4LPkPSxS9oHuAoYaGYHELR9NYLkjf1x4PhS48qMNTzuRwB9wmXuD3sHTJTH+WHs7wIHmFlf4CvgD5AysSOpE3AcsCxiXLLFXqYanQjMbLWZTQvf5xGckPYh6AHtiXC2J4DTEhJgBSR1BE4CHokYnQpxNwGOBB4FMLNdZraZFIg9lAHUl5QBNCBoEj0pYzezj4FNpUaXF+upwHgz22lmS4BFBL0DJkRZsZvZO2ZWGA5OImiSHlIg9tDdwO+ByDtwkir28tToRBBJUhZwMPAl0NbMVkOQLIA2CQytPP8iOKiKI8alQtzdgPXAY2Gx1iOSGpICsZvZSuAOgl90q4EtZvYOKRB7hPJiLasnwH2qObZYjALeDN8nfeySfgKsNLOZpSYlfexQSxKBpEbAC8CvzWxrouPZE0knA+vMLDvRsVRCBtAf+K+ZHQxsJ3mKUioUlqefCnQFOgANJZ2f2KiqjMoYl5T3jkv6I0Gx7tiSUWXMljSxS2oA/BH4c1mTyxiXNLGXqPGJQFImQRIYa2YvhqPXSmofTm8PrEtUfOU4DPiJpKXAeOBoSU+T/HFD8ItnhZl9GQ4/T5AYUiH2Y4ElZrbezAqAF4FDSY3YS5QXa1Q9ASaapIuAk4Hz7LuHnJI99u4EPx5mhv+zHYFpktqR/LEDNTwRSBJBWfV8M7srYtIrwEXh+4uAl6s7toqY2R/MrKOZZRFUNH1gZueT5HEDmNkaYLmkXuGoY4B5pEDsBEVCQyQ1CI+dYwjqlVIh9hLlxfoKMEJSXUldgZ7A5ATEVy5JxwPXAj8xsx0Rk5I6djObbWZtzCwr/J9dAfQP/xeSOvbdzKzGvoDDCS7DZgEzwteJQEuCOyq+Dv+2SHSsFXyGYcBr4fuUiBs4CJga7veXgOYpFPtfgQXAHOApoG6yxg6MI6jLKCA4+VxaUawExRffAAuBE5Iw9kUE5ekl/6sPpErspaYvBVolY+zlvbyJCeecq+VqdNGQc865PfNE4JxztZwnAuecq+U8ETjnXC3nicA552o5TwTOOVfLeSJwzlU7Sd0kPSrp+UTH4jwRuAhhO+pPRQxnSFpf0h9CNcfyeYzz3yjpd/GKZ29jkHR5uH+Pihj3y3DcsdUXZfkk1Zf0UXU0k2xmi83s0lLbryPp47DlV1eNPBG4SNuBAyTVD4ePA1YmIhAzOzQR242jvgRPWu8Huxsqu5SgpdbZCYwr0ijgRTMrqqoVSjpQ0mulXmW23mpmuwiehj6nqrbvouOJwJX2JkE/CAAjCR6nB0DSS5KyFfTgNTocNyjsUaqepIbhtAMq2oCk30u6Knx/t6QPwvfHhI3rIWmbpCwFPZw9HK73nYgkhaQ/hr0+vQf0KmM7DSW9Lmmmgh7HzgnHZynoCeuJMPbnwxMzks6XNFnSDEkPlvw6rmB8hTFEODDcl73D4auA/wHFZra2ov1Vjc4jbJsoYh89Eu67sZKOlfSZgt7PftCmfln724J2eE4u9aqo0b6XwjhcNfJE4EobT9BIVj2CX7FfRkwbZWYDgIHAVZJamtkUgoa1bgH+ATxtZj/owq+Uj4EjwvcDgUYKWok9HPik1Lw9gfvMrA+wGTgDQNIAggb5DgZOBwaVsZ3jgVVm1s+CHsfeipjWC3jIgt6wtgJXStqP4NfoYWZ2EFAEnFfB+GhiKLEf8BzQW1LTcH2fE7RplHCS6gDdzGxpxOgeBL3N9SVIYOcSfEe/A64vYzUV7e/S22sp6QHgYEl/iJg0h4r3o4sDL4tz32NmsxR04jMSeKPU5Ksk/TR834ngJL0RuAmYAuQT/NLdk2xggKTGwE5gGkFCOKKM5ZeY2YyI5bLC90cAEyxspVLSK2VsZzZwh6TbCRrui0wyy83ss/D90+F284EBwJSg8VHqEzTj3KSc8S2iiKGkC8ONZrY4LBb5PfBvYF+C4iIkNTez3LKWryatCBJtpCVmNhtA0lyCLjBN0my++x4iVbS/v8fMNgJXlDG+SNIuSY0t6FXQVQO/InBleYWgp67IYqFhBO31DzWzfsB0oF44uQXQCGgcMa5cFrT1vxS4hOBX8SfAcIJ23eeXmn1nxPsivv/jpcIWE83sK4IT+Gzg75IiOw4pvawRdCLyhJkdFL56mdmNFYzfYwyhvnxXD5BH8Mv5CYLiopLxd0exnnj6lh9+d5H7vjhiuJgyfkTuYX/Hoi5BUnbVxBOBK8sY4KaSX4OhpkCume2Q1BsYEjHtIeBPBD1K3R7lNj4mKGL4mCARXAHMsOibw/0Y+Gl4p0tj4JTSM0jqAOwws6cJElv/iMmdJQ0N348EPiWoqDyzpDJTUgtJXSoYv8cYQpEn/H8CvwwrZA8EZiloh7+3wjuOJP1F0j2S/iqpo6QPJf1G0rOlPl87BXf5/D6s77g8rO84oIz1tJb0WLi+MWFR3G7h1Uh6WCRYKXvY39GuoyVQ0jGQqyaeCNwPmNkKM7un1Oi3gAxJs4CbCToXR9KFQKGZPQPcBgySdHQ47Y3w5FCWT4D2wBdhZWk+P6wfqCjGacCzBO3Wv1DOsgcCkyXNIGgT/paIafOBi8LP04Kga815wA3AO+H4d4H2FYyPJoaSOOaEcb9mZl+E4/cn6LRnA0Hdyh2S9gEyCYpphgD9gJfM7G6C7hsjHUxwl88/CBL1wwQV0F1Kr8fM1hN0vHMncFU5J9p3COoAKqui/R2t4fywSNLFW6I7RPCXv6r7RVC+PSfRcUTEMwo4Inw/hqCIrVv4/nqgdzjt6VLL/YGgEjcTeDQcdzPQuYz1NCLoaGdsBXEcDDyV4H3xItAr0d9JbXt5ZbFzibcBuEzSBmAuQZFZS4J6mIOBryS1AtaUWq4H8BVBHURJ3UqWmS0LK3cj13MvwVXN2ZKGmdnE0kGY2fSwGCrdqvBZgmiFdy69ZGYLq3vbtZ33UOacc7Wc1xE451wt54nAOedqOU8EzjlXy3kicM65Ws4TgXPO1XKeCJxzrpbzROCcc7WcJwLnnKvlPBE451wt9/8DMxXFNQxMQwAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(winds, delta_p_max, '.-')\n",
"plt.grid()\n",
"plt.xlabel(r'Max. wind speed $M_{tan\\ max}$ (m s$^{-1}$)')\n",
"plt.ylabel(r'Max. pressure deficit $\\Delta P_{max}$ (kPa)')\n",
"plt.title('Max. tornadic wind speed vs. pressure deficit')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# A27\n",
"## Parts a, e, and j"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABC8UlEQVR4nO3dd3xUZfb48c9JI4XQOwFDE6QJUgVRrCD2CrZVsf5W13V13VV37XVXd117w17gq4hdUFRiwwIq0iMRBEInlBSSkHJ+fzw3ZAgpk5DJzGTO+/W6r7n9npuB58x97r3PI6qKMcaYyBUV7ACMMcYElyUCY4yJcJYIjDEmwlkiMMaYCGeJwBhjIpwlAmOMiXCWCIwBRGSsiGT6ua6IyAsisl1EfhCRMSKS7sd254nIJ/sfbXCJyO0i8qo33lVEckUkOthxmbqzRGCqJCK/i0i+9x99k1f4NQ12XCHgMOBYIEVVh6vqV6rau6aNVPU1VT2ubFpEVER6BjLQQFPVNaraVFVLgh2LqTtLBKYmJ6lqU+AQYBjwz4oriEhMg0cV3BgOAH5X1bwGPGZAhMJ3Z4LPEoHxi6quA2YC/WHPr9mrRGQFsMKbd6KILBCRHSIyV0QGlm0vIn8XkXUikiMi6SJytDd/uIjMF5Fs76rjv978fapqvCuUY7zx20Vkuoi8KiLZwEUi0lxEnhORDd6x7q6qykJEEkTkRa96Zykuyfku7yQib4nIFhFZJSLXePMvAaYAh3pXSndUjFVEuojIDG/bLBF5zJt/kYh87Y1/6a3+i7efiRWO38T7O/b3mdfWu0JrJyJtROQDb51tIvKViNT4/7ksVu/72Ai8ICItvX1t8f4eH4hIis823UTkC++7mw208VmW6v1biKn4Hfl8T2XVSPHe95XlxT1PRNrXFLMJPPs1YPwiIl2ACcAMn9mnAiOAfBE5BHgeOAmYD5wPvCcivYFU4GpgmKquF5FUoKyAfhh4WFVf8aqd+uO/U4CzgD8ATYCpwCagJ5AEfACsBZ6uZNvbgB7ekIRLcmXnGgW8D7wLnAOkAJ+KSLqqPiciJcClqnqYt/5Yn22jveN+DlwAlABDKx5cVQ8XEQUOVtWMSpYXisgM7/j/8GafDXyhqptF5D4gE2jrLRsJ+NteTAegFe7KJgpIBF7w9h+N+x4fw32/AK8D3wLH4b7vD72/TW1dCDQHugCFwCAgvw77MfXMrghMTd4RkR3A18AXwL0+y+5T1W2qmg9cBjytqt+raomqvoT7zz4SVxg2AfqKSKyq/q6qv3n7KAJ6ikgbVc1V1e9qEdu3qvqOqpYCzYDjgWtVNU9VNwMPAZOq2PZs4B4v/rXAIz7LhgFtVfVOVd2tqiuBZ6vZl6/hQCfgBi+OAlX9uhbn5Ot1XCIoc643D9zfrSNwgKoWefcp/E0EpcBtqlqoqvmqmqWqb6nqLlXNAe4BjgB3Mxj397jFW/9LXJKsiyKgNdDT+zfyo6pm13Ffph5ZIjA1OVVVW6jqAar6R6/QL7PWZ/wA4Hrvkn+Hlzy6AJ28X7zXArcDm0Vkmoh08ra7BDgQWO5VFZxYi9gqHj8W2OBz/KeBdlVs26nC9qsr7KtThXO5GfCnGqMLsFpVi/07hWp9DiSIyAgROQD3C/ptb9kDQAbwiYisFJEba7HfLapaUDYhIoki8rSIrPaq2b4EWnhXN52A7RXuh6ymbl4BPgamich6Efm3iMTWcV+mHlkiMPvD9xfoWtwv7BY+Q6KqTgVQ1de9qpQDvO3+5c1foarn4ArsfwHTRSQJyMNVWQB7qlzasreKxy8E2vgcv5mq9qsi9g24QrtM1wr7WlXhXJJVdYIff5O1QFeph5uw3pXOG7irgnOBD7xf7Khqjqper6rdcdVx14l338WfXVeYvh7oDYxQ1WbA4d58wf2dWnrfSZmuVG2v7w1XDVV2PkWqeoeq9gVGASfiqvVMkFkiMPXlWeBK79eriEiSiJwgIski0ltEjhKRJkABrl64BEBEzheRtl6ht8PbVwnwKxDv7SMW97RSk6oOrqobgE+A/4hIMxGJEpEeInJEFZu8Adzk3ShNAf7ks+wHINu7oZogItEi0l9EhlW+q738gCs87/f+BvEiMrqKdTcB3WvY3+vAROA8yquFym7M9xQRAbJxf7O6PsKZjPtOdohIK9z9EwBUdTXuns8dIhInIofhEk9VFgCTRCRWRIYCZ/rEfKSIDPCSejauqsgeOw0BlghMvVDV+bj7BI8B23HVFhd5i5sA9wNbgY24X/83e8vGA0tEJBd343iSV6++E/gj7gmddbhfmjW98PUHIA5Y6sUwHVePXpk7cFUcq3AJ5BWfcynBFXaDvOVbvTia13B83217Amu8mCdWsfrtwEte9dPZVezve9y5d8LnhjbQC/gUyMXdyH1CVdMARGSmiNyM//4HJODO8ztgVoXl5+JuEm/DJYmXq9nXLbgb8Ntxf+PXfZZ1wH0n2cAy3D2nV2sRpwkQsY5pjDEmstkVgTHGRDhLBMYYE+EsERhjTISzRGCMMREu7JqYaNOmjaamptZp27y8PJKSkmpeMYTZOYSGcD+HcI8f7Bxq68cff9yqqhXfxQHCMBGkpqYyf/78Om2blpbG2LFj6zegBmbnEBrC/RzCPX6wc6gtEanyjXCrGjLGmAhnicAYYyKcJQJjjIlwlgiMMSbCWSIwxpgIF9BEICLjxXVLmFFZe+leK5WPeMsXer1cGWOMaUABSwReU7OP43qN6gucIyJ9K6x2PK4VxV7A5cCTgYrHGGNM5QL5HsFwIMPr5g8RmYbrY3apzzqnAC97Xex9JyItRKSj17a8MSYIVKG0FEpK9v4sG8qWVzdUXKds2nd+ZeO1+VywoBW5uW7cdyhbXtVQ03J/hrK/U22XV5y/cmVXvvmm8uWVjR92GBx3XP1/54FMBJ3ZuyvATFyb5jWt0xnXscceInI57oqB9u3bk5aWVqeAcnNz67xtqLBzaBiqUFQkFBdHUVQUxe7dQlFR1J4hOzuaRYt+oqQkiuJioahIqh0vKRFKSoTiYql0vKqhtLTiJ3tNl427ArV82hV45du4dXzHR1NaWoyqW6ds/dJSCfafvhYGBjuAelBTv0R7O+ecNcTFraz3KAKZCCr7F1Wx8wN/1kFVnwGeARg6dKjW9U08exMxNOzPORQXQ24u5OS4T9+hbF5eHuTnVz7s2rXvvIICKCzce9i9u37PuaKYmPIhNnbvad8hOtoNZeNxceXjvssrG6Kiqp63fn0mBxyQste8qKiax0XKp6saRNxQcf3qxitbVt2nCPz8848MHTpkz7TvfkRASkuIKt6991BStOdTiovcZ1H5tJQUEVW0242XLS8p9pYV75kvJcVQVD6P4iKkuNjNLylGiqv/pNiNZ2/bRrOmSXumKS1x/8hLyj+luLh8Xp9rYOyt9f/vsd73WC6TvfuETQHW12Ed04gUF8P27bGkp8P27f4N2dnlhXxhof/HEoGEhKqHdu3Kx5s02XeIi6t6fnr6IoYMGUBsrCvI4+Kocdy3cJcg//BOS8tg7NiU+t9xScm+WdZ3qGlexYxczdB/+3aaxsa6rF3ZUFpa/+dXUVTU3tm8uqxeSZbX+DiiWrbYN7tX9hkTAwcfHJDTCGQimAf0EpFuuK4GJ+G6vPP1HnC1d/9gBLDT7g+Ep127YMMG2LjRfVY1bN0KqlV14QtJSdCyZfnQrRs0awbJydC0aflnxcF3flKSK9zj4gJX4KalZRGWF2aqUFhI7I4dsGpV+SWU72dl8/Ly3JdcdklVNl5xujaZujIxMeUZNz6+8kyclAStWpGflETTzp3dF+3PUFOmrmzat4D3xktFKCkpcUNREcUFBZQUFFBSWOiG3bspLRv3mVc2Xuozvn7tWtq2alU+b/duSouK3Gdenvv0mdf7oIMCUiEWsESgqsUicjXwMRANPK+qS0TkSm/5U8BHwARc/7a7gIsDFY/ZP6qweTP8+iukp5d/rlgB69a5X+0VxcRA+/bQsSMccACMHOnGt29fwfDhvfYq8MuGuLiGP7ewkp+/7yXTtm2wc6f7EvwZioqoOhVXEB/vsmtioiuAExNdlm3WDDp0cOOJieXzfccTEtz28fF7j1c2r6zgj/L/QcYllVQxqiqlRUUU79pF8a5dFOXl7Rkvzs93nwUFFOfnU7Jz557pkrJ5+fluPd95hYXlBX1BAaXFxX7HuA8Rops0IToujqi4OIpU2bFlC1GxsUTFxREdF0d0fDxxzZq56dhYops0Icpbv3W/fnU/djUC2vqoqn6EK+x95z3lM67AVYGMwdROXp4r5MuGskL/119dWVOmSRPo2RP69IFjjnEFfMWhdevK/1+npa1j7NheDXdSoUgVduyALVvcsHlz+fiWLZCVVXk9WU2/uOPjXSHtOxxwwD7zVmzcSK/Bg8svoXwvp8rGExNdNm8AJYWF7N62zRXcubkU5eZSlJfnPqsY37puHbOefLK80Pc+tRYFtURFER0fT0xCAtEJCcQkJBDjTce2aeMK7fh4Ypo0ceMJCUQ3aUJMfPyeZb6fMfHxrgAvW99nPCoujqjYWMTnMjVU7vmFXTPUpn5lZcGXX8IXX0BaGixcWP64GkDXrtC7N5x/vvs88EA3dO3qqi5NBXl57hJp3TpYv778c9OmvQv6LVvcDZPKJCe7LFp2mdS3b/l4q1b7Xka1bAnNm7tC3s9LqnVpafSq5wKopLCQwu3bKdyxg8IdO9idnc3u7GyKcnL2jO/Ozma3N12Uk8PunTvZnZNDqR9356NiYohNTia2aVNiEhOhpIQmLVvStHNnYhITiUlMJDYpac94TGIisT7jMV5BX1bwxyQkEBUXt1fBHKksEUSYrVtdwZ+W5gr/hQvd/IQEGDUKbrsN+vVzhX7Pnm6+8eTmwsqVsGoVnT7/HGbPLi/0ywp838umMsnJro6sbVtITYXhw92479Cunfts08b9qg8yLS2lcOdOCrZupSAri4KsrD0FfOH27ezesYPCnTv3KvhL8vOr3J9ERxOXnExss2bEJScT17w5SR06ENesGbHJycR582OaNiU2KYnYpk3d4DMeXSHJhcqv6cbAEkEjl5EBv/ziCv60NFi82M1PTHQF/913wxFHwLBhrronopWWujvav/3mCvyVK/ce37x5z6oHgqs26dgROnWCgw5ydWSdO7vpzp3Lx5OTg3ZKFRXl5bFr40YKly9nVW7uXgV9QVZW+fS2bWhJSaX7iGvWjCYtWxLXogUJ7drR4sADadKiBU1atqRJ8+ZuWfPmbvAK/ZjERPvlHcIsETQyJSUwbx58+CE8/7z7kQqu4D/sMDjnHBg7FoYOjeAbs4WFsGwZLFrkhqVLXYG/atXedfBRUa4OrHt3OOUU9+kNc9esYdSpp9bq5mYgqSpF2dns2rTJDRs37jWe740X5ebu2eZb7zMqNpb41q2Jb92ahPbtadWvn5tu02bP/PhWrWjSqpW7iWl1go2OJYJGYMcO+PhjV/jPnOmqf6KioH9/GDAArr/eFf6xscGOtIGVlsLq1eUFftmQnu4yJrhs2Lu3+0V/4omuoO/Rw3127Vplttydm9vgSaA4P5/cdevIXbvWDZmZbli7lrwNG/atmhEhoU0bEtq3J7lbN9qPHEli+/Yktm9P+oYNjDz6aBLatCE2Odl+rUc4SwRhSNX9iP3wQzd8840r11q1gvHj4YQTYNw4d78xYqi6R5u++spdEi1a5OrBcnLK10lNhYED4bTTXIYcMAB69QqpDLk7J4edv/1G7po15HiFfd7ateSsXUvB1q17rRuTmEjTLl1o1r07HUePJrFDBzd4hX1C27ZEVXFuv6el0bx77Zo3MI2XJYIwsWkTzJ0LL73Ui4sucj90wZVrf/ubK/xHjoygJ3lKStzNj6++cne/v/66vA6/ZUv3h7nwwvICv3//kKqr352dzc7ffmNnRob79Ib8TZvKVxIhsX17mnbpQqcxY2japQtNU1L2fDZp2dJ+yZt6YYkgBJWUwJIlruAvG377zS2Lj+/AccfBzTfDhAmQEoBWAkJSQQH88IMr+L/6yv1Ryn7tp6a6S6AxY9zQu3fw23DwFOfns335cnauWMGO334j2yv887ds2bNOdEICzbt3p8OIETTv0YNmPXrQLDWVpM6d93lSxphAsEQQAnJy4PvvXRXP3Lnw3Xflb+q2awejR8OVV7rPnJyvOe64I4IbcENZsgSmT3ePac6bV94SXP/+7sWGsoI/RLJhaXExO1esIGvxYrIWLSJr8WJ2ZmTsefomJiGBZj160GHUKJr37EnzHj1o3qMHSZ06ISFy09lEJksEDUgVMjNd/X7ZUFadXVrqfsQOGADnnuse7Rw92rW14/vjNi1tn8ZZG5clS+DNN92wdKk7+eHD4c9/doX+6NHuZkiQqSo5q1fvKfCzFi1i+/LllBQUABDXvDmt+/en89ixtO7fn5Z9+pDYoYMV+CYkWSIIgJIS+P338sJ+2bLyT5+n92jdGgYPhltucQX/iBHuBdGIs3QpvPHG3oX/4YfDY4/BGWe4Nm2CTEtL2b58ORvnzmXTDz+w8eefeX/XLgCi4+Np1bcvPc8+m9YDBtC6f3+aduli9fcmbFgiqCVV1yxDVa1r/vorLF++9+PonTq5VgIuvth99u3rnlZs2zZ45xF0S5e6gv+NN0K28M9dt46N337Lxm+/ZdN331G4YwcAzXv1ImHwYPofdxytBwygeY8eRDVQmzzGBELE/OstLIStW+NYuXLv5s8rGy/73LXLPa3jW9Bv3AhFRfvuPznZvWTasycce6wr6MsK/Ij8lV+Z7dvhyScZ9swz7rEn38L/9NPdHzCIdmdnu1/7c+ey4dtvyV2zBoCEdu3odMQRdBg1ig4jRpDQti1paWn0tOYNTCMRMYngnXdg0qRRtd6udevy1jT79Km8lc2OHV2jjaYKW7bA//7nCvzsbIoGDgyJwl9V2b58OWtnz2bjt9+ybfFitLSUmMRE2g0bRu/zzqPDyJE069HDqnlMoxYxiWDoUPjLX9I5+ODeVTaTXtl4xDyXHwgbNsB//gNPPukus846C26+mQXbtwe1sbCcNWtY/dFH/P7hh2SvXIlER9N6wAD6XXEFHQ49lDYDB1b5IpYxjVHEJIIePeDkkzcwdmzvYIfS+K1dC//+Nzz7rGtq+dxz4aabXD0ZuNbvGlj+1q2smTWL3z/8kCyvydV2Q4fS54IL6HLccTRp0aLBYzImVERMIjANYOVKuO8+eOklN33hhXDjjS4LB0FRbi5rP/2U3z/8kE3ffYeWltKid28GXX89Bxx/PElBvidhTKiwRGD23/LlLgG89pprmvnyy127F127NngoJbt3s/6rr1j94YesS0ujpLCQpJQU+l56KQeccAItevZs8JiMCXWWCEzdbdniXvSaNs3dVPnzn11Tp506NXgoxbt2kfHWWyx/8UV2bdxIk1at6HHGGRxwwgm0Ofhgu9lrTDUsEZi6+eormDTJvVTx97/DddcF5cWI3dnZ/Dp1KumvvELh9u20GzaMYbfcQsfRo+2GrzF+skRgaqe0FO6/370O3aMHfPQRHHxwg4dRkJXF8ldeYcXUqRTl5tLpiCPod+mltD3kkAaPxZhwZ4nA+G/zZrjgAvjkE9fV2dNPN3jTznnr17PsxRf5bfp0Snbvpuu4cfS79FJalj2RZIypNUsExj9ffukK/6wslwAuu6xBm3rO/v13lk6Zwqr33weg28kn0/eSS2iWmtpgMRjTWFkiMNUrLXVPBN16q2s/o4GrgnasWMHip55izccfE92kCb0mTeKgCy8kKQg3pI1prCwRmKoFsSqoZPdulk6ZwuKnnyYmPp5+l11G7/PPJz6i+t80pmFYIjCV++ILV/hv2wbPPAOXXtpgVUHbly3j23/8gx3p6aSedBJDbrzR3vw1JoAsEZi9lZbCvffCbbe5qqCZMxusKsj3KqBJixYc/uijpBx1VIMc25hIZonAlNu50zUMN3t2g1cF7XUVcOKJDLnpJrsKMKaBWCIwTkmJaxxuzpwGrQra6yqgeXO7CjAmCCwRGOe229wTQU884R4NbQDbly/n25tvtqsAY4IsoD1pi8h4EUkXkQwRubGS5eeJyEJvmCsiDf+KqoHp0+Gee9xVwJVXBvxwWlLCoieeYNbEiRRs3crhjz7KqH/9y5KAMUESsCsCEYkGHgeOBTKBeSLynqou9VltFXCEqm4XkeOBZ4ARgYrJVGLRIrjoIhg50vUaFuDqoO3Ll7Pl3/9mQ2amXQUYEyICWTU0HMhQ1ZUAIjINOAXYkwhUda7P+t8BKQGMx1S0bRuceio0awZvvQVNmgT0cBvmzuXLP/0JjY21ewHGhBBR1cDsWORMYLyqXupNXwCMUNWrq1j/r0CfsvUrLLscuBygffv2Q6ZNm1anmHJzc2natGmdtg0V9XUOUlLCgBtvpMUvv7Dgf/8ju2/feoiuagWLF7Pt2WeJad+e+MmTadahQ0CPF2jh/m8p3OMHO4faOvLII39U1aGVLlTVgAzAWcAUn+kLgEerWPdIYBnQuqb9DhkyROtqzpw5dd42VNTbOdxwgyqoPvts/eyvGms+/VSnDhyoM886Swu2b7fvIQSEe/yqdg61BczXKsrVQN4szgS6+EynAOsrriQiA4EpwCmqmhXAeEyZqVPhgQfgj390N4gDaPWsWXx93XW07NePo6ZMsfsBxoSgQCaCeUAvEekmInHAJOA93xVEpCswA7hAVX8NYCymzM8/wyWXwJgx8NBDAT3UqvffZ+4NN9Bm4ECOevZZ4po1C+jxjDF1E7CbxapaLCJXAx8D0cDzqrpERK70lj8F3Aq0Bp7wuhIs1qrqsMz+27LF3Rxu3RrefBPi4gJ2qN9mzOD7W2+l/bBhHPH448QkJgbsWMaY/RPQF8pU9SPgowrznvIZvxQIbN2EcYqKYOJE2LQJvv4a2rcP2KFW/N//Me/OO+k4ejRjHnmEmPj4gB3LGLP/7M3iSHHDDa75iJdegqGBu+ha/sor/HT//XQ64gjGPPQQ0QF+JNUYs/8sEUSCl16Chx+Ga6+FP/whYIdZ+txzLPjvf+lyzDGMeuABogNY9WSMqT+WCBq7efPgiivgqKPck0IBsujJJ1n02GMccPzxHHrffUTFxgbsWMaY+mWJoDH75Rc47jjo0AH+7/8gpv6/blVl4aOPsuTpp0k9+WRG3n03UdHR9X4cY0zg+FUyiMhQYAzQCcgHFgOfquq2AMZm9ocqnHkm7NgBr78ObdoE4BDKgv/+l2XPP0+PM85g2G23WRIwJgxV+x6BiFwkIj8BNwEJQDqwGTgMmC0iL3nvAphQM2UKZGTATTfB8ccH5BCrP/qIZc8/T6+JExl+++2WBIwJUzVdESQBo1U1v7KFIjII6AWsqee4zP5YtQquuw6OPhruvjsgh8jfsoX599xD64MPZsg//oFEBbRFc2NMAFWbCFT18RqWL6jXaMz+Ky11zUpHRcHzz7vPeqaqzLvrLorz8+2egDGNgL/3COKBS4B+wJ63g1R1coDiMnX1v//Bl1/CCy9A18DU2q3+6CMyP/uMQddfT/Pu3QNyDGNMw/H35+IrQAdgHPAFrgG5nEAFZepo6VK4+WY4+WS48MKAHMK3SqhPgI5hjGlY/iaCnqp6C5Cnqi8BJwADAheWqbWiIlf4Jye7zucD0NOYVQkZ0zj5+2B5kfe5Q0T6AxuB1IBEZOrmvvtg/nzX/3CA2hFa/eGHZH72GYP/+lerEjKmEfE3ETwjIi2BW3BNSTf1xk0o+PFHuOsuOO88OOOMgBwif8sW5t97L60PPpjeAWymwhjT8GpMBCJyKtACGK6qHwP2UzCUFBS49oPat4dHHw3IIVSVeXfeaVVCxjRS1SYCEXkC96TQXOAuERmuqnc1SGTGP7fc4m4Sz5oFLVsG5BCrP/yQzM8/tyohYxqpmq4IDgcOVtUSEUkEvgIsEYSKr76C//zHNSo3blxADlFWJdRm0CCrEjKmkarpqaHdqloCoKq7gPp/FMXUTW6ue3GsWzd48MGAHMKqhIyJDDVdEfQRkYXeuAA9vGkBVFUHBjQ6U7UbbnBNSXzxBTRtGpBD/O5TJdSsW7eAHMMYE3w1JYKDGiQKUzsffwxPPQV//avrhD4A8rds4cd77rEqIWMiQLVVQ6q6WlVXA33Lxn3mBaZJS1OtmJwcmDwZ+vVzj4wGgKrywx13UFJYaFVCxkQAf98svkVEjiqbEJG/A6cEJiRTnV6PPAKbN7vuJwPUKfzvH37IujlzGHjNNVYlZEwE8PeFspOBD0TkBmA80MebZxrS9Om0//RTuP12GDIkIIfYq0roggsCcgxjTGjxKxGo6lYRORn4FPgROFNVNaCRmb3Nnw8XXkh2nz40u/nmgBzCqoSMiUw1vVCWA/gW+HG4N4vPFBFV1WaBDM54Vq2CE06Adu1YfM89jApQx/Abvv6adXPmMPiGG6xKyJgIUlPHNMkNFYipQlaW62qyqAg++ojdmzYF7FC/Tp1KfJs29D7vvIAdwxgTemrqszi1huUiIin1GpEpV1AAp5zirgjefRcOCtzTvHkbNrDhq6/ocfrpRAXoisMYE5pqukfwgIhEAe/i7g1swfVQ1hM4EjgauA3IDGSQEam01DUm9803MG1awN4XKPPbW2+hqvQ488yAHscYE3pqqho6S0T6AucBk4GOwC5gGfARcI+qFgQ8ykj0t7/Bm2/CAw/AxIkBPVRpcTG/vfUWHUePpmnnzgE9ljEm9NT41JCqLgX+0QCxmDKPPuoak7v6arj++oAfbv2XX5K/eTND//nPgB/LGBN6/H2hzDSUt9+GP/8ZTj3VdUQfgC4nK1rxxhsktGtH5yOOCPixjDGhxxJBKPn2Wzj3XBg+HF57DRrgOf7cdevY8PXX7iZxjL/vFxpjGpOAJgIRGS8i6SKSISI3VrPeMBEpEZHIvVO5YgWcdBJ07gzvvw+JiQ1y2N+mT0dE6BGgLi6NMaHPr0QgIp/5M6/C8mjgcVzjdH2Bc7wbz5Wt9y/gY39iaZQ2b3bvCojAzJnQtm2DHLa0qIjfZsyg45gxJHXq1CDHNMaEnpreLI4HEoE2Xuf1ZRXWzYCaSo7hQIaqrvT2NQ3XUN3SCuv9CXgLGFa70BuJXbvclcC6dTBnDvTq1WCHXpeWRsHWrfQ866wGO6YxJvTUVCl8BXAtrtD/yWd+Nu7XfnU6A2t9pjOBEb4riEhn4DTgKKpJBCJyOXA5QPv27UlLS6vh0JXLzc2t87YBUVJC/9tuo/W8eSy54w62FhRADfHV5zlkPf00US1asKK0lIwG/LuE3PdQB+F+DuEeP9g51CtVrXEA/uTPehW2OQuY4jN9AfBohXXeBEZ64y/iGrOrdr9DhgzRupozZ06dt613RUWql12mCqqPPOL3ZvV1Djlr1uhrffvqwscfr5f91UZIfQ91FO7nEO7xq9o51BYwX6soV2uqGjpKVT8H1onI6ZUkkRnVbJ4JdPGZTgHWV1hnKDBN3COSbYAJIlKsqu9UF1fYy8pyL4l99hncfDP86U8NHkLG9OlIVBQ9Tt/nazXGRJiaqoaOAD4HTqpkmQLVJYJ5QC8R6QasAyYB5+61A9U9TVyKyIvAB40+CSxZAiefDJmZ8PzzcPHFDR5Cye7drHz7bTodcQSJHTo0+PGNMaGlpiYmbvM+a11aqWqxiFyNexooGnheVZeIyJXe8qfqEG94e+cduOAC19l8WhocemhQwlg3Zw4FWVn0PPvsoBzfGBNa/H189F4RaeEz3VJE7q5pO1X9SFUPVNUeqnqPN++pypKAql6kqtNrEXv4KC2FO++E005zLYjOnx+0JADuTeLEjh3pOHp00GIwxoQOf18oO15Vd5RNqOp2YEJAImpscnPh7LPhttvc1cCXX7qXxoIkZ/VqNn33HT3PPNN6IDPGAP73WRwtIk1UtRBARBKAJoELq5FYtcr1J7BkiWtE7i9/aZC2g6qT8eabSHQ03e0msTHG428ieBX4TERewN0kngy8FLCoGoM5c+Css6CkxL0tfNxxwY7I3SR+5x06H3kkie3aBTscY0yI8Lfz+n+LyCJcRzQC3KWqkdskRHVU4fHH4dpr4cAD4b33oGfPYEcFwNpPP6Vw+3Z7k9gYsxe/m5tU1ZnAzADGEv4KC+Gqq+C551yzEa++Cs2aBTuqPTLeeIOklBQ6jhoV7FCMMSGkpj6Lv/Y+c0Qk22fIEZHshgkxDKi6l8PGjHFJ4J//dI+KhlASyF61is3z5tHzzDORKGt93BhTrqYrgj8AqGpyA8QSflRh1iy46y7Xl0DHjq57yRDs9zfjzTeRmBi6n3pqsEMxxoSYmn4avgk1NzkdcVTh3XddBzITJriWQ594AlauDMkkUFJYyMp33iHlqKNIaKAmro0x4aOmK4IoEbkNOFBErqu4UFX/G5iwQlRpKbz1Ftx9NyxcCN27w5Qp7v2AuLhgR1elNbNns3vnTnrZm8TGmErUdEUwCSjAJYzkSobIUFzsuo7s39+9HFZYCC+/DOnpcMklIZ0EwN0kbtqlC+1HjKh5ZWNMxKmpraF04F8istB7aiiyFBW5J3/uvRcyMlwimDbNVf+EyVu5OzMy2PLjjwy67jq7SWyMqZS/j4/+JCLPAZ1U9Xivy8lDVfW5AMbWsHbtguXLYelSWLbMff7wA6xfD4MHw4wZ7i3hMCtMM958kyi7SWyMqYa/ieBF4AXgH970r8D/AeGXCHbuLC/oyz6XLoXVq91NYICYGNdl5KhRcNFF7oZwkJuGqIvdeXmsfO89Uo49lvjWrYMdjjEmRPmbCNqo6hsichPsaWK6JIBx1b/33+fQyZNh69byeU2aQJ8+riXQyZOhb1839OgR8vX+/vj5X/+iKDubriHQvIUxJnT5mwjyRKQ1rp0hRGQksDNgUQVC585sHzKEDkceWV7gp6aGTV1/XWxbvpzohARSjjkm2KEYY0KYv4ngOuA9oIeIfAO0BULvgfnqHHIIy2+8kQ5jxwY7kgZRsH07O9LTOeiii4gKs/saxpiG5W+jcz+JyBFAb1yjc+mqWhTQyMx+yZw9Gy0u5oDx44MdijEmxNXUeX1VjdYfKCI1dV5vgmj1rFkkp6bSok+fYIdijAlxNV0RlHVa3w4YhevIHuBIII3qO683QZK/ZQub582j3xVXIGH4tJMxpmHV9ELZxQAi8gHQV1U3eNMdgccDH56pizWffIKWllq1kDHGL/7eRUwtSwKeTcCBAYjH1IM1M2fS4sADaR4iHeIYY0Kbv08NpYnIx8BU3COkk4A5AYvK1Fnehg1s+flnBl5zTbBDMcaECX+fGrpaRE4DDvdmPaOqbwcuLFNXaz52PYhatZAxxl+16arybcAK/xC3euZMWvXrR/IBBwQ7FGNMmLA3jRqR3LVr2bZ4MV3tasAYUwuWCBqR1bNmAVYtZIypHb8TgYgkiEjvQAZj9s/qmTNpM2gQSZ06BTsUY0wY8SsRiMhJwAJgljc9SETeC2BcppZ2rlzJjvR0qxYyxtSav1cEtwPDgR0AqroASA1EQKZu1syaBSJ0HTcu2KEYY8KMv4mgWFXDq9npCKKqrJ45k3ZDh5LYrl2wwzHGhBl/E8FiETkXiBaRXiLyKDA3gHGZWtjx669kr1zJAccfH+xQjDFhyN9E8CegH1AIvI7rlObamjYSkfEiki4iGSJyYxXrjBWRBSKyRES+8DMe42PNrFlIdDRdjj022KEYY8JQjS+UiUg08J6qHkN5n8U18rZ7HDgWyATmich7qrrUZ50WwBPAeFVdIyJWr1FLZdVC7UeMIL5Vq2CHY4wJQzVeEahqCbBLRJrXct/DgQxVXamqu4FpwCkV1jkXmKGqa7xjba7lMSLetiVLyF271qqFjDF15m8TEwXAIhGZDeSVzVTV6lo26wys9ZnOBEZUWOdAIFZE0oBk4GFVfbnijkTkcuBygPbt25OWluZn2HvLzc2t87ahouI57JwxA6KjWR0fz9owObfG+D2Em3CPH+wc6pO/ieBDb6iNynpE0UqOPwQ4GkgAvhWR71T11702Un0GeAZg6NChOraO/Q6npaVR121Dhe85aGkp7951F50OO4yxEyYEN7BaaGzfQzgK9/jBzqE++dv66Et12Hcm0MVnOgVYX8k6W1U1D8gTkS+Bg4FfMTXa+ssv7Nq4kYOvvTbYoRhjwphfiUBEVrHvr3lUtXs1m80DeolIN2Adrg+Dcyus8y7wmIjEAHG4qqOH/InJuCYlouLiSDnyyGCHYowJY/5WDQ31GY8HzgKqfURFVYtF5GrgYyAaeF5Vl4jIld7yp1R1mYjMAhYCpcAUVV1c25OIRKUlJaz5+GM6H344sU2bBjscY0wY87dqKKvCrP+JyNfArTVs9xHwUYV5T1WYfgB4wJ84TLktP/5IwdatdLWnhYwx+8nfqqFDfCajcFcIyQGJyPhl9cyZxCQk0Pnww2te2RhjquFv1dB/fMaLgd+Bs+s9GuOX0qIi1s6eTeexY4lJTAx2OMaYMOdv1ZDdjQwhG7//nsLt261ayBhTL/ztj+DPItJMnCki8pOIHBfo4Ezl1syaRWzTpnQ67LBgh2KMaQT8bXRusqpmA8cB7YCLgfsDFpWpkhYVsfbTT0k56iiimzQJdjjGmEbA30RQ9pbwBOAFVf2Fyt8cNgFWuHw5RTk5HBBGbxIbY0Kbv4ngRxH5BJcIPhaRZNxz/6aB5f/4I3HNm9Nh5Mhgh2KMaST8fWroEmAQsFJVd4lIK1z1kGlAxQUFFCxcSPcTTyQqNjbY4RhjGgl/rwgOBdJVdYeInA/8E9c5jWlA67/8Ei0stGohY0y98jcRPInrk+Bg4G/AamCf5qJNYK167z2ikpNpN3RozSsbY4yfatN5veI6lnlYVR/G3ixuUDmrV7MuLY3E0aOJivG3Rs8YY2rmb4mSIyI3ARcAY7xuKK2SugGlv/46UdHRJFmTEsaYeubvFcFEXMf1k1V1I673MWsoroHszslh5YwZdJ0wgejmte0x1BhjqudXIvAK/7eAsjeYtgJvByoos7ff3nqL4l276HPBBcEOxRjTCPnbxMRlwHTgaW9WZ+CdAMVkfJQWF/Pra6/RbuhQWvXtG+xwjDGNkL9VQ1cBo4FsAFVdgWtqwgRY5uefk7d+Pb3tasAYEyD+JoJCVd1dNuF1LblP15Wm/qW/8gpJKSl0tu4ojTEB4m8i+EJEbgYSRORY4E3g/cCFZQCyFi9my08/0fu884iKjg52OMaYRsrfx0f/DlwKLAKuwHU/OSVQQRkn/ZVXiElKosfppwc7FNNIFRUVkZmZSUFBQbBDqbXmzZuzbNmyYIexXwJxDvHx8aSkpBBbi2ZoakwEIhIFLFTV/sCz+xGfqYVdmzezetYsDjznHOuc3gRMZmYmycnJpKamIhJeDQrn5OSQnBze77XW9zmoKllZWWRmZtKtWze/t6uxakhVS4FfRKTr/gRoamfF1KloSQm9zzsv2KGYRqygoIDWrVuHXRIwlRMRWrduXesrPH+rhjoCS0TkByCvbKaqnlyroxm/FOfnk/HGG6QcdRRNu3QJdjimkbMk0LjU5fv0NxHcUes9mzr7/f33Kdyxgz5/+EOwQzHGRIBqE4GIxANXAj1xN4qfU9XihggsUqkqy195hZYHHUTbIUOCHY4xJgLUdI/gJWAoLgkcD/wn4BFFuI1z55K9ciV9/vAHu2Q3xjSImhJBX1U9X1WfBs4ExjRATBFt+csvE9+mDV3Hjw92KMY0qLfffhsRYfny5XvmPf3003To0IFBgwbRvXt3XnzxxeAF6Jk1axa9e/emZ8+e3H///VWuN3nyZNq1a0f//v0rXf77778zYsSISpdFR0czaNAg+vfvz1lnncWuXbvqJfaq1JQIispGrEoo8Hb+9hsbvv6aXpMmER0XF+xwjGlQU6dOZejQoUybNm3PvIULF3L77bezYMECpk+fzvXXXx/ECKGkpISrrrqKmTNnsnTpUqZOncrSpUsrXfeiiy5i1qxZdTpOQkICCxYsYPHixcTFxfHUU0/tT9g1qikRHCwi2d6QAwwsGxeR7IBGFoHSX32VqLg4ek2cGOxQjGlQubm5fPHFFzz33HNMnTp1z/xFixZx0EEHAZCSkkJJSUmwQgTghx9+oGfPnnTv3p24uDgmTZrEu+++W+m6hx9+OK1atfJrvytXrmTw4MHMmzdvn2VjxowhIyMDgFNPPZUhQ4bQr18/nnnmmbqfSAXV3ixWVWvXoIEU7tjBqvfeo9tJJxHv5z8eY+rTtdfCggX1u89Bg+B//6t5vXfeeYdjjjmGgQMHkpSUxE8//cQhhxzCokWL6NOnD6rKI488woknnliv8Y0ZM4acnJx95j/44IMcc8wx+8xft24dXXwe6U5JSeH777/frxjS09OZNGkSL7zwAoMGDdprWXFxMTNnzmS8V1X8/PPP06pVK/Lz8xk2bBhnnHEGrVu33q/jg/+Pj5oAy3jzTUoKCqyVURORpk6dyuWXXw7A2WefzdSpU2nbti25ubmMGzeO2NhYhg8fzuOPP16vx/3qq69qtb7rsXdv+/NQx9atWznllFN466236Nev3575+fn5e5LCmDFjuOSSSwB45JFHePtt1xXM2rVrWbFiRegnAhEZDzwMRANTVPX+CsubA68CXb1YHlTVFwIZUygqLSri16lT6XDoobTo1SvY4ZgI5c8v90DIysrihx9+YMaMGQBMnDiRI444grFjx3L00UfvU8++ceNGJk6cyDHHHENGRgajRo1i9uzZ3H777fTv35877riDbdu20aJFC66++mr+9re/cdddd3Hrrbfy9NNP79UGT22vCFJSUli7du2e6czMTDp16lTnc2/WrBldunThm2++2SsRlN0j8JWWlsann37Kt99+S2JiImPHjq23NqIClgi8fo0fB44FMoF5IvKeqvreWbkKWKqqJ4lIWyBdRF7zbfI6Eqz55BPyN21i+G23BTsUYxrc9OnTmTBhAk2auA4Qu3XrRocOHfjpp584+OCD91n/559/5vTTT2fy5MlccMEFXHbZZbRo0YLVq1fTsmVLioqKaNGiBd999x133HEHXbt25frrr+e5557bpyG22l4RDBs2jBUrVrBq1So6d+7MtGnTeP311+t87nFxcbzzzjuMGzeOpk2bcu6551a57s6dO2nZsiWJiYksX76c7777rs7HrcjfZqjrYjiQoaorvYJ9GnBKhXUUSBZ3bdUU2AZE1NNJqsryl18mOTWVTmPs6VwTeaZOncr7779PamrqnmHZsmXMnTuXgQMH7rP+ggULGDduHEVFRbRu3ZqoqCgWL17MgAEDuOWWW/j73//OhRdeSOfOncnNzWXlypXExMTQtB4ab4yJieGxxx5j3LhxHHTQQZx99tl7/ZKfMGEC69evB+Ccc87h0EMPJT09nZSUFJ577rlK95mUlMQHH3zAQw89VOWNZ4Dx48dTXFzMwIEDueWWWxg5cuR+n08ZqazOq152LHImMF5VL/WmLwBGqOrVPuskA+8BfYBkYKKqfljJvi4HLgdo3779EN/Hy2ojNze3Xv4x1KfdK1ey9T//ofnEiSQdfniN64fiOdSWnUPwlcXfvHlzevbsGexwauWqq67i0UcfZeHChXz55Zdcc801XHbZZTz77LM8+uijZGdns23bNnr16sXChQu58cYbefvttznkkEMYE2I/tkpKSogOQF8jGRkZ7Ny5c695Rx555I+qOrTSDVQ1IANwFu6+QNn0BcCjFdY5E3gIEFwzFquAZtXtd8iQIVpXc+bMqfO2gfLVX/6ib4wcqUV5eX6tH4rnUFt2DsFXFv/SpUuDG8h+yM7ODnYI+y1Q51DZ9wrM1yrK1UBWDWUCvk1npgDrK6xzMTDDizPDSwR9AhhTSMlbv561s2fT88wziUlMDHY4xpgIFchEMA/oJSLdRCQOmISrBvK1BjgaQETaA72BlQGMKaT8+vrrIMKB1dwgMsaYQAvYU0OqWiwiVwMf4x4ffV5Vl4jIld7yp4C7gBdFZBGueujvqro1UDGFktzMTH6dOpWu48eT1LFjsMMxxkSwgL5HoKof4fo39p33lM/4euC4QMYQilSV+ffei4gw6C9/CXY4xpgIF8iqIVOFzM8+Y/0XXzDg6qvtasAYE3SWCBpYUV4eP953Hy0OPND6IzbGhARra6iBLXriCXZt3MjoBx8kqsJbjsYYEwx2RdCAtqenk/7KK/Q44wzaDh4c7HCMMQawRNBgtLSUeXfeSVyzZgy67rpgh2OMMXtYImggv82YwdYFCxh8/fU0adEi2OEYE3IaW1eVqampDBgwgEGDBjF06L4tO4RTV5WmHhRs28aC//6XtkOG0O3UU4MdjjEhqbF1VQkwZ84cFixYwPz582t1nFDrqtLUgwX/+Q9FeXkMu/XW/erEwpjGqjF2VVkbId1Vpdl/m3/8kZXvvEPfSy6hRZi18mgiTBD7qmyMXVWKCMcddxwiwhVXXLGnB7aKrKvKRq5k927m3XknSZ060f/KK4MdjjEhqzF2VfnNN9/QqVMnNm/ezLHHHkufPn04vEJT8xHRVWWkS3/5ZXZmZHD4Y49Z66Im9AWpr8rG2lVl2fx27dpx2mmn8cMPP+yTCBp9V5WRLm/9ehY99RQpRx1FypFHBjscY0JWY+yqMi8vj9LSUpKTk8nLy+OTTz7h1ltv3We9UOmq0hJBgMy/914Ahtx0U5AjMSa0TZ06lYULF5KamrpnXlZWFq1bt+b888/fZ/0FCxZw2mmn7dNV5WWXXcYtt9zCww8/zJYtW1i7dm1Au6osKSlh8uTJ+3RVOWXKFAoKCjjttNMAV89/7rnn7qnnr6isq8pjjz2WpKQkTjmlYo++zvjx43nqqacYOHAgvXv3rteuKi0RBEDm55+zbs4cBl1/PUlVXDYaY5y0tLRarZ+RkcGBBx7It99+u+eJot9//52uXbvSr18/HnzwQbKyshg8eDDXXHMNd999N2+88QZpaWmMHTt2v+OdMGECEyZMqHTZRx+VN7b8yy+/VLuf1NTUPTeaW7RosdcTQ7m5ufus36RJE2bOnFmXkGtkiaCeFeXlMf+ee2jeqxd9Lrgg2OEY0+iUdQI/cOBARo8eDcArr7wCUOV7BjfccEPDBBem7D2Cerb4ySfZtXEjw265xRqVM8aEBUsE9WjHihUsf+UVup9+Ou2GDAl2OMYY4xdLBPWkpLCQH26/ndikJGtUzhgTVuweQT0oLijgq2uuYeuCBYz697+Jb9ky2CEZY4zfLBHsp+Jdu/ji6qvZ9MMPjLjzTlJPOCHYIRljTK1YItgPRXl5fPHHP7Llp58Yec89dK/i+V9jjAlllgjqaHdODmlXXknWokUcev/9diVgjAlblgjqYPfOncy54gq2LVvG6AceoOu4ccEOyRhj6swSQS0V7tjB55ddxs5ff2XMQw+RctRRwQ7JGGP2iz0+WgsF27bx2cUXszMjgzGPPmpJwJh6FC5dVU6ePJl27drRv3//aterqUtL66oyDOVv2cJnF11EzurVHPH443Su0JysMWb/hENXlQAXXXTRPk1jV1TbLi0rsq4qQ9CuzZv57OKLyV2/nrFPPknHUaOCHZIxjUq4dFUJcPjhh9OqVatq16ltl5bWVWWIy9uwgc8mT6Zg61aOfPppazrCNFo/3ncf29PT63WfLXv39qsp9nDpqtJftenS0rqqDHG569bx2cUXs3vnTo6aMoU2lXSSYYzZf+HSVaW//O3S0rqqDGE7V65k7ezZrJg2jeKCAo567jla13BjyJhwF6xOlMKpq0p/+dulpXVVGUJUlR2//sra2bNZ+8kn7PztNwDaDBrEsH/+k5ZeHaUxpv6FU1eV/vK3S8tG31WliDwPnAhsVtV9fk6Lu056GJgA7AIuUtWfAhVPRarKtsWLWTt7NmtmzyZ3zRokKoq2Q4YwZOJEuhxzDInt2zdUOMZErHDqqhLgnHPOIS0tja1bt5KSksIdd9yxp+qmrKvKTp06Vdulpa/G3lXli8BjwMtVLD8e6OUNI4Anvc+A0dJStvz0E2tmzybz00/JW78eiY6m/YgR9L34YlKOPpr4eqhvM8b4L9y6qvR9qqki364qq+vSEiKkq0pV/VJEUqtZ5RTgZXV3Vb4TkRYi0lFVNwQinnVffMGmf/6TDTt3EhUbS4dRoxjwxz/S+cgjadKiRSAOaYwJAOuqsv4F8x5BZ2Ctz3SmN2+fRCAilwOXA7Rv377WvyAAitatIyolhWannUZ8//6QkMAaYE2FGzKhLjc3t07nH0rsHIKvLP7mzZtXerM0HJSUlIRt7GUCdQ4FBQW1+vcZzESw77NUsO8zV4CqPgM8AzB06FCt6+VdWufO9XJpGEz1dXkbTHYOwVcW/7Jly0hOTg52OHWSk5MTtrGXCdQ5xMfHM3jwYL/XD+abxZlAF5/pFGB9kGIxxpiIFcxE8B7wB3FGAjsDdX/AGFO1yl5+MuGrLt9nIB8fnQqMBdqISCZwGxALoKpPAR/hHh3NwD0+enGgYjHGVC4+Pn7Po5qVvflqwouqkpWVRXx8fK22C+RTQ+fUsFyBqwJ1fGNMzVJSUsjMzGTLli3BDqXWCgoKal3ghZpAnEN8fDwpKSm12sbeLDYmgsXGxtKtW7dgh1EnaWlptbohGopC5RysGWpjjIlwlgiMMSbCWSIwxpgIJ+H26JiIbAFW13HzNsDWegwnGOwcQkO4n0O4xw92DrV1gKq2rWxB2CWC/SEi81V1aLDj2B92DqEh3M8h3OMHO4f6ZFVDxhgT4SwRGGNMhIu0RPBMsAOoB3YOoSHczyHc4wc7h3oTUfcIjDHG7CvSrgiMMcZUYInAGGMiXMQkAhEZLyLpIpIhIjcGO56aiEgXEZkjIstEZImI/Nmb30pEZovICu+zZbBjrYmIRIvIzyLygTcdVufgdaM6XUSWe9/HoeF0DiLyF+/f0GIRmSoi8aEev4g8LyKbRWSxz7wqYxaRm7z/2+kiMi44Ue+tinN4wPt3tFBE3haRFj7LgnYOEZEIRCQaeBw4HugLnCMifYMbVY2KgetV9SBgJHCVF/ONwGeq2gv4zJsOdX8GlvlMh9s5PAzMUtU+wMG4cwmLcxCRzsA1wFBV7Q9EA5MI/fhfBMZXmFdpzN7/i0lAP2+bJ7z/88H2Ivuew2ygv6oOBH4FboLgn0NEJAJgOJChqitVdTcwDTglyDFVS1U3qOpP3ngOrvDpjIv7JW+1l4BTgxKgn0QkBTgBmOIzO2zOQUSaAYcDzwGo6m5V3UEYnQOuleEEEYkBEnE9AYZ0/Kr6JbCtwuyqYj4FmKaqhaq6CtfHyfCGiLM6lZ2Dqn6iqsXe5He4nhkhyOcQKYmgM7DWZzrTmxcWRCQVGAx8D7Qv68nN+2wXxND88T/gb0Cpz7xwOofuwBbgBa96a4qIJBEm56Cq64AHgTXABlxPgJ8QJvFXUFXM4fr/ezIw0xsP6jlESiKorOulsHhuVkSaAm8B16pqdrDjqQ0RORHYrKo/BjuW/RADHAI8qaqDgTxCrxqlSl49+ilAN6ATkCQi5wc3qnoXdv+/ReQfuOrf18pmVbJag51DpCSCTKCLz3QK7vI4pIlILC4JvKaqM7zZm0Sko7e8I7A5WPH5YTRwsoj8jquOO0pEXiW8ziETyFTV773p6bjEEC7ncAywSlW3qGoRMAMYRfjE76uqmMPq/7eIXAicCJyn5S9yBfUcIiURzAN6iUg3EYnD3ZR5L8gxVUtcB7LPActU9b8+i94DLvTGLwTebejY/KWqN6lqiqqm4v7mn6vq+YTXOWwE1opIb2/W0cBSwucc1gAjRSTR+zd1NO5+U7jE76uqmN8DJolIExHpBvQCfghCfDUSkfHA34GTVXWXz6LgnoOqRsQATMDdpf8N+Eew4/Ej3sNwl4YLgQXeMAFojXtiYoX32SrYsfp5PmOBD7zxsDoHYBAw3/su3gFahtM5AHcAy4HFwCtAk1CPH5iKu6dRhPu1fEl1MQP/8P5vpwPHBzv+as4hA3cvoOz/9FOhcA7WxIQxxkS4SKkaMsYYUwVLBMYYE+EsERhjTISzRGCMMRHOEoExxkQ4SwQmIohIiYgs8FrgfN+31Uc/t08TkaHe+Ee13b6KfQ4WkSk1r7nXNtNEpNf+HtsYX5YITKTIV9VB6lrg3AZcVdcdqeoEdQ3P7a+bgUdruc2TuLabjKk3lghMJPoWr0EvERkuInO9BuXmlr1BLCIJ3q/vhSLyf0BC2cYi8ruItBGR1Aptzf9VRG73xq8RkaXe9tMqBiAiycBAVf3Fm75dRF4SkU+8/Z8uIv8WkUUiMstrbgTgK+AYryVRY+qFJQITUbw23o+mvImR5cDh6hqUuxW415v//4Bd6tqNvwcYUstD3QgM9ra/spLlQ3Fv+vrqgWuy+xTgVWCOqg4A8r35qGop7u3Ug2sZjzFVskRgIkWCiCwAsoBWuA5CAJoDb3q/7B/CdQwCrg+CVwFUdSGueYnaWAi85rX0WVzJ8o645q19zVTXMNwiXAcys7z5i4BUn/U241oSNaZeWCIwkSJfVQcBBwBxlN8juAv3y7s/cBIQ77NNTe2vFLP3/yHfbU/A9Yo3BPixkqqc/ArrAxTCnl/9RVre/ksprjls3+Pk1xCbMX6zRGAiiqruxHXd+Fev3r05sM5bfJHPql8C5wGISH9gYCW72wS0E5HWItIE17QwIhIFdFHVObgbuy2AphW2XQb0rONpHAgsqeO2xuzDEoGJOKr6M/ALrmnsfwP3icg3uOqYMk8CTUVkIa4w36dJYK8a505cz3Ef4O434O3nVRFZBPwMPFTxKSNVXQ40924a+01E2uOubjbUZjtjqmOtjxoTJCLyFyBHVf1+l8DbJltVnwtcZCbS2BWBMcHzJN59gVrYQXkH7sbUC7siMMaYCGdXBMYYE+EsERhjTISzRGCMMRHOEoExxkQ4SwTGGBPh/j/E31DHiXkbRwAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"# Import some needed libraries\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"\n",
"# Pressure deficits in kPa\n",
"R0 = 25. # Core radius\n",
"R_inner = np.arange(0, R0 + 5, 5) # A set of inner radii\n",
"R_outer = np.arange(R0, 130, 5) # A set of outer radii\n",
"Delta_p_max = np.array([0.1, 0.5, 1.0]) # Pressure deficits (in kPa) from parts a, e, and j\n",
"colors = [\"blue\", \"red\", \"brown\"]\n",
"# Plot the pressure deficit\n",
"for dpmax, col in zip(Delta_p_max, colors):\n",
" # Rankine Combined Vortex model\n",
" # Plot inner core pressure (Eqn. 15.41)\n",
" RCV_p_inner = dpmax * ( 1 - 0.5 * (R_inner / R0) ** 2 )\n",
" plt.plot(R_inner, RCV_p_inner, color = col, label = r\"$\\Delta P_{max} =$ %0.1f kPa\"%dpmax)\n",
" # Plot outer region pressure (Eqn. 15.43)\n",
" RCV_p_outer = (dpmax / 2) * ( R0 / R_outer) ** 2\n",
" plt.plot(R_outer, RCV_p_outer, color = col)\n",
"\n",
"plt.legend()\n",
"plt.xlabel('Radius (m)')\n",
"plt.ylabel('Pressure deficit (kPa)')\n",
"plt.title('Pressure deficit vs. radius')\n",
"plt.grid()\n",
"plt.gca().invert_yaxis()"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"array([ 9.03507903, 20.20305089, 28.57142857])"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"rho = 1.225 # Standard air density in kg per cubic meter at sea level\n",
"# Use Eqn. (15.45)\n",
"M_tan_max = np.sqrt(Delta_p_max * 1000. / rho) # Remember to multiply Delta_P_max by 1000 to convert from kPa to Pa.\n",
"M_tan_max # meters per second"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0, 0.5, 'Wind speed (m s$^{-1}$)')"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEWCAYAAABmE+CbAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABbMUlEQVR4nO2dd3iUxdbAfyc9tIQaIAFCE6kCARQRAaliARsqFqzIFbvftV4V9dquDRAVUcQONpqoKJoEUZo06YhSTEIvgQQIpJzvj9kkSwrsJtnsbjK/55ln920zZ7a8552ZU0RVsVgsFoslwNsCWCwWi8U3sArBYrFYLIBVCBaLxWJxYBWCxWKxWACrECwWi8XiwCoEi8VisQBWIVj8GBFJF5FmLpwXKyIqIkFl3P5NIvJrKevoKSKbykomX0REeotIstP2OhHp7T2JLMVhFUIlxXEzzS05InLMafs6b8tXEBFJFJHbnPepajVV3eItmcoCVV2gqq1yt0Vkm4j086ZMnkZV26pqorflsBSmTJ+YLP6DqlbLfS8i24DbVPUn70lk8QdEJEhVs7wth8Uz2BGC5SREpJuILBKRVBHZKSITRCTE6biKyCgR2SwiB0XkTRERx7FAEXlVRPaJyFYRuct5qkZEIkRksqPeFBH5r4gEOo7dJCK/isgrjnq3isiFjmPPAT2BCY4RzAQnWVo43l8kIitF5LCIJInIGBf7+4iIfFVg3zgRGX86mYuo61wR+V1EDjlez3U6VktEpojIDkf/Zjr2502niMjHQGPgG0c/HxKRb0Xk7gLtrBaRoUW0P1dE7iqw7w8RuVwMr4vIHod8q0WknYufkYrIaBHZDGx2+oySHJ/3chHp6XR+uIh84OjneqBrgfryRkGO8/7rdKzg9NLDjs89TUQ2iUhfV2S2lBBVtaWSF2Ab0M/xPg44BzN6jAU2APc5navAHCASc/PaCwxyHBsFrAdigJrAT47zgxzHZwLvAFWBesBS4A7HsZuATOB2IBD4F7ADEMfxRMwohgKytHC87w20xzzkdAB2A0Mdx2Kd5ShQRxPgKFDDsR0I7ATOcVHmXx3vawEHgRscn921ju3ajuPfAp87PpdgoJeT3MlFfReO7WHAEqfts4D9QEgRfbkR+M1puw2QCoQCA4Hlju9NgNZAAxd/HwrMc/Qx3LHveqC2o68PAruAMMexF4EFjvMbAWuL6yPwAfBfp2N5nwfQCkgCGjp9j829/X+pyMWOECwnoarLVXWxqmap6jbMzbBXgdNeVNVUVf0HSAA6OvYPA8aparKqHsTcGAAQkSjgQoxyOaKqe4DXgWuc6t2uqu+qajbwIdAAiHJR7kRVXaOqOaq6GphahNxFXbcdWAEMdey6ADiqqotdlDmXi4DNqvqx47ObCmwELhGRBo56RqnqQVXNVNX5rvQLmAW0FJGWju0bgM9V9UQR584AOopIE8f2dcB0VT2OUbbVgTMxSnaDqu50UQaAF1T1gKoeA1DVT1R1v6Ovr2KUTu5ayDDgOcf5ScB4N9pxJttRbxsRCVbVbar6dwnrsriAVQiWkxCRM0RkjojsEpHDwPNAnQKn7XJ6fxTIXY9oiHmiy8X5fRPMk/FOx3RUKkbZ1CuqXlU96nhbDRcQkbNFJEFE9orIIcxopaDcxfEZ5okeYLhj21WZc2kIbC+wbzsQjXlKPuBQkm7huJl/AVwvIgEOOT8u5tw0zEgkV2FdA3zqOBYPTADeBHaLyCQRqeGGKM7fJSLyoIhscEw/pQIR5H/eBX8HBT8Xl1DVv4D7gDHAHhGZJiINS1KXxTWsQrAU5G3Mk21LVa0BPIaZYnCFnZjpolwaOb1PAo4DdVQ10lFqqGpbF+s+XVjez4DZQCNVjQAmuiH3l0BvEYkBLiNfIbgj8w6MAnGmMZDiqKeWiES6IEtR/fwQ87TfFzN6WXSK66cC14pIdyAcM4IzFauOV9U4oC1wBvBvF+QpJJdjveBhzEigpqpGAofI/7x3cvJ33/gU9R4Bqjht1z+pUdXPVPU8zGerwEtuyGxxE6sQLAWpDhwG0kXkTMxcvqt8AdwrItGOm9/DuQcc0xM/Aq+KSA0RCRCR5iJy2mkdB7uBU/kcVMc8hWeISDfMk75LqOpezBrFFGCrqm4ogczfAWeIyHARCRKRqzFz+HMc9XwPvCUiNUUkWETOd7WfDgWQA7xKMaODAnI0AZ7BTC3lAIhIV8coKhhzE87ATMmUhOpAFmb9KEhEngScRxtfAI86+hoD3F1EHbmsAgY7Ft3rY0YEOGRuJSIXiEioQ95jpZDZ4gJWIVgK8n+Ym2ka8C5mIdRV3sXcQFcDKzE3pyzy/8Q3AiGYheeDwFeYdQJXGAdc6bBcKWpO+k7gGRFJA57E3JTc4TOgH/mjg1xckllV9wMXYxZY9wMPARer6j7HKTdg5vE3AntwuvEV4AXgP44pqv9z2v8RZtH8k1N1wjHFNL2IvtTAfD8HMVM4+4FXAETkMRH5/lT1FuAHjIL701FXBidPET3t2L8V83s4lRL7GPgDs9D8Iyf/3kIx61D7MNOJ9TAjVouHyLXgsFjKHDFmoxNVteBUisVNRORGYKRj+sRi8Qh2hGApMxz254MdUybRwFMYyxdLKRCRKpgR0CRvy2Kp2JS7QhCRMBFZ6nCYWSciTzv21xKReWIcnuaJSM3yls1SagQzXXAQM2W0ATN9YykhIjIQM1e/m8LTWRZLmVLuU0YiIkBVVU13LHD9CtwLXI5ZFHxRRB7BWC88fKq6LBaLxVJ2lPsIQQ3pjs1gR1FgCMa8Dsfr0PKWzWKxWCozXgluJyYWzHKgBfCmqi4Rkahcz0lV3SkiRTn/ICIjgZEA4eHhcY0aNSrqtNOSk5NDQIB/L6H4ex/8XX6wffAF/F1+KP8+/Pnnn/tUtW6hA96Mm4GJq5IAtANSCxw7eLrr4+LitKQkJCSU+Fpfwd/74O/yq9o++AL+Lr9q+fcBWKa+FstIVVMxDkGDMO70DQAcr3u8J5nFYrFUPrxhZVQ314VfRMIxDjQbMWEHRjhOG4EJ6mWxWCyWcsIbawgNgA8d6wgBwBeqOkdEFgFfiMitwD/AVV6QzWKxWCot5a4Q1IQm7lTE/v2Y4F0Wi8ULZGZmkpycTEZGhrdFcYuIiAg2bNjgbTFKhaf6EBYWRkxMDMHBwS6db1NoWiwWAJKTk6levTqxsbEYdyH/IC0tjerVq3tbjFLhiT6oKvv37yc5OZmmTZu6dI1/22pZLJYyIyMjg9q1a/uVMrAUj4hQu3Ztt0Z8ViFYLJY8rDKoWLj7fVqFUAnIPnGCxDvvZM3bb3tbFIvF4sNYhVAJyD5+nB3z5/P31197WxSLxeLDWIVQCdj5228AHN21iyM73cmrbrFYKhNWIVQCkuPjCalRA1TZNmeOt8WxWE7LjBkzEBE2btyYt++dd96hfv36dOzYkWbNmvHBBx94T0AHc+fOpVWrVrRo0YIXX3yx2PNuueUW6tWrR7t27Yo8vn379mKPBQYG0rFjR9q1a8dVV13F0aNHy0T2orAKoYKTk5nJjgULiLngAurGxbF11qzcWFEWi88ydepUunTpwrRp0/L2rV69mjFjxrBq1Sq++uorHnzwQS9KCNnZ2YwePZrvv/+e9evXM3XqVNavX1/kuTfddBNz584tUTvh4eGsWrWKtWvXEhISwsSJE0sj9imxCqGCs2f5cjIPHyb6ggtoNmQIh7duZf+aNd4Wy2IplvT0dObPn8/kyZOZOnVq3v41a9bQunVrAGJiYsjOzi6uinJh6dKltGjRgmbNmhESEsI111zDrFlFR9w5//zzqVWrlkv1btmyhU6dOvH7778XOtazZ0/++usvAIYOHUpcXBxt27Zl0qSySaZnHdMqOMkJCQSGhtKge3dysrNZ9txzbJ01izodOnhbNIsPc999sGpV2dbZsSOMHXv682bOnEm/fv3o0KEDVatWZcWKFXTu3Jk1a9Zw5plnoqqMHz+eiy++uEzl69mzJ2lpaYX2v/LKK/Tr16/Q/pSUFJzD78fExLBkyZJSybBp0yauueYapkyZQseOHU86lpWVxffff8+gQYMAeP/996lVqxbHjh2ja9euXHHFFdSuXbtU7VuFUIFRVVLi44k65xyCqlQBIKZvX7Z/9x2dH36YwJAQL0tosRRm6tSpjBw5EoBhw4YxdepU6tatS3p6OgMHDiQ4OJhu3brx5ptvlmm7CxYscOv8oqZeS+PHsXfvXoYMGcLXX39N27Zt8/YfO3YsTzn07NmTW2+9FYDx48czY4ZJWZ6UlMTmzZutQrAUT+qmTRzZsYO2d9yRt6/pkCFs/+47dsyfT6P+/b0oncWXceVJ3hPs37+fpUuXMn36dACuvvpqevXqRe/evenbt2+hefhdu3Zx5ZVXcumll7Ju3TrOPfdc5s2bx5gxY2jXrh1PP/00Bw4cIDIykrvuuouHHnqIZ599lieffJJ33nnnpBg/7o4QYmJiSEpKyttOTk6mYcOGJe57REQEjRo14rfffjtJIeSuITiTmJjITz/9xKJFi6hSpQq9e/cukxhUViFUYJLj40GE6N698/bV796d8Lp12TJrllUIFp/jq6++YvDgwYSGhgLQtGlT6tevz4oVKzjrrLMKnb9y5UouueQSHnroIYYOHcrtt99OZGQk27dvp2bNmmRmZhIZGcnixYt5+umnady4MQ8++CCTJ08uFPDN3RFC165d2bx5M1u3biU6Oppp06bx2WeflbjvISEhzJw5k4EDB1KtWjWGDx9e7LmHDh2iZs2aVKlShY0bN7J48eISt+uMXVSuwCQnJFDnrLMIr1Mnb19AYCCxF1/MjgULyDhwwIvSWSyFmTp1Kt988w2xsbF5ZcOGDSxcuJAORax7rVq1in79+pGZmUnt2rUJCAhg7dq1tG/fnieeeIKHH36YESNGEB0dTXp6Olu2bCEoKIhq1aqVWtagoCAmTJjAwIEDad26NcOGDTvpyX7w4MHs2LEDgGuvvZbu3buzadMmYmJimDx5cpF1Vq1alTlz5vD6668Xu0ANMGjQILKysujQoQNPPPEE55xzTqn7A3aEUGE5umsXB9evp+P99xc61nTIEDZMmcL2776j1fXXe0E6i6VoEhMT3Tr/r7/+4s4772TdunV5Fkjbtm2jcePGtG3blldeeYX9+/fTqVMn7rnnHv773//yxRdfkJiYSG+nkXNJGTx4MIMHDy7y2HfffZf33tlaqiiaNGnC2rVrAYiMjDzJwig9Pb3Q+aGhoXz//fclEfmUWIVQQUlOSAAg+oILCh2LbNmSmm3asGXWLKsQLH7N5MmTSUtLo2PHjnkLrx9//DFAsX4K//73v8tLPL/DThlVUJLj46nepAk1iomD3vTSSzm4fj2pmzeXs2QWi8VXsQqhAnIiLY09S5cSc8EFxZrBxV50ERIUxNZTzFNaLJbKhVUIFZCdv/5KTlYW0X36FHtOWK1aNOzZk21z5pDjZY9Pi8XiG1iFUAFJjo8ntFYt6hTwdCxI00sv5djevexatKh8BLNYLD6NVQgVjNxgdtG9ehEQGHjKc6N79yakRg07bWSxWACrECoce5YvJzMtjZgirIsKEhgSQpMLLyT555/JLMK0zWKxVC6sQqhgJMfHExgaSn0XHVWaDh1K9vHj/PPDDx6WzGKx+DpWIVQgVJXk+Hjqd++eF8zudNRu354aTZuydfZsD0tnsVh8HasQKhCpGzdydOdOl6aLchERml56KXuWLSM9OdmD0lksFl/HKoQKRHJCAojQsFcvt66LvfhiELGjBIvPUNFSaMbGxtK+fXs6duxIly5dCh2vtCk0RaSRiCSIyAYRWSci9zr2jxGRFBFZ5ShFBwixFEtRwexcoWrDhkR168bW2bNtek2LT1DRUmgCJCQksGrVKpYtW+ZWOxU9hWYW8KCqtgbOAUaLSBvHsddVtaOjfFd8FZaCHNm5k4Pr1xNzCme0U9F0yBDSk5LYu2JFGUtmsbhHRUyh6Q6VKoWmqu4Edjrep4nIBiC6vOWoaKScIpidKzTq149lzz7L1tmzqRcXV5aiWfwRL+bQrIgpNEWEAQMGICLccccdeRnhClKpU2iKSCzQCVgC9ADuEpEbgWWYUcTBIq4ZCYwEiIqKcjtcbi7p6eklvtZXcO7D/q+/JrBePVb+8w/880+J6gvu0IEtc+ZwtEcPpBzSa1a078Bfye1DRERE3g0x9MQJAsr4CTznxAmOF3HDLcjHH3/MzTffTFpaGkOGDOHDDz8kPDyc9PR0+vfvT1BQEHFxcbz44oukpaWRnZ1d5I3cXZzDVRekqPqPHj1KZmZm3rFjx46dtO3MDz/8QIMGDfLSZDZu3JgePXrkHc/JyWHPnj1ccsklfPLJJzRv3vykenNzQXTv3p1hw4aRlpbGyy+/zJw5cwCTQnPVqlV069atUNsZGRmu/0ZV1SsFqAYsBy53bEcBgZhprOeA909XR1xcnJaUhISEEl/rK+T24fjhwzq1Qwdd8fLLpapv56JF+mmbNrr122/LQLrTU5G+A38mtw/r16/3riCqum/fPq1Tp45mZGSoquqWLVu0UaNGOmfOHB04cGCh83fu3Kk9evTQl156SW+88UadOHGiXnHFFbpmzRpVVR0zZozec889+uSTT+qePXv0pptu0qSkJL355pv1xIkTJ9V13nnn6VlnnVWozJs3r0hZFy5cqAMGDMjbfv755/X5558/bR+feuopfbnAf3XNmjXasmVL7devn77zzjsnHatatWqhOhISErRHjx565MgRVVXt1atXsb/For5XYJkWcU/1yghBRIKBr4FPVXU6gKrudjr+LjDHG7L5I7nB7NwxNy2KqG7dqFK/Pltnzya2mKQfFosnqYgpNI8cOUJOTg7Vq1fnyJEj/Pjjjzz55JOFzvOFFJrlrhDExGOeDGxQ1dec9jdQs74AcBmwtrxl81dyg9nVLuIP4w4SEEDTSy9l/XvvcWzvXsLr1i0jCS0W15g6dSqrV68mNjY2b9/+/fupXbs21xeRzGnVqlUMGDCgUArN22+/nSeeeIJx48axd+9ekpKSPJpCMzs7m1tuuaVQCs333nuPjIwMLrvsMsCsAwwfPjxvHaAguSk0+/fvT9WqVRkyZEiR5w0aNIiJEyfSoUMHWrVq5dcpNHsANwBrRGSVY99jwLUi0hFQYBtwhxdk8ztyg9k16tfvtMHsXKHppZeybtIktn37La1vuqn0AlosblBRU2j+8ccfp6yn0qbQVNVfgaKytlgz0xKwZ9kyl4PZuUKNpk2p3aEDW2fN4swRI4pNsGOx+AI2hWbZYj2V/Zzk+HgCw8Ko3717mdXZbMgQUv/8k4NOXqIWi6XiYxWCH6POwezCw8us3saDBhEQHMzWmTPLrE6LxeL7WIXgx2QlJ3N0164ymy7KJTQykkb9+vH3jBlkHCzkCmKxWCooViH4MRmrV4MI0W4Gs3OFtqNGkXX0KBvef7/M67ZYLL6JVQh+TMaaNdQ56yzCSumuXhSRLVoQe/HF/PnZZxzbu7fM67dYLL6HVQh+ypEdO8hMSirz6SJn2v/rX+RkZrLu3Xc91obFYvEdrELwU5Idwew8qRCqN2lCs6FD+euLLziyY4fH2rFYLL6BVQh+SkpiIkFRUdRo2tSj7bQbNQqAtR6MwW6xWHwDqxD8kBNpaexZupQwRwRET1K1YUNaDBvGlpkzObx9u8fbs1gs3sMqBD9kx4IF5GRlEdq+fbm01/b22wkIDmbtW2+VS3sWi7+k0LzllluoV69esekvczldqs1Km0LTUnpSHMHsQjw8XZRLeN26nDF8ONu+/ZbUzZvLpU1L5cYfUmgC3HTTTcydO/eU57ibarMg5ZlC06sJcizuk33ihAlmN2AAGQHlp8/b3Hormz//nDUTJtBz3Lhya9fiHZa/8AIHN20q0zprtmpF3KOPnva83BSa8+bN46qrrmLMmDGASaE5bNgwwDdSaAKcf/75bNu27ZTnOKfaBPJSbbZp06bI87ds2cIVV1zBpEmT6Nq160nHevbsyerVqwGTQjMpKYmMjAzuvffeYrOwuYMdIfgZe5YtIzM9vcS5k0tKaGQkZ44YQdJPP3Fg3bpybdtSuSgqhSZQLik0c4PkOZeffvqpVPUWlWozJSWlyHM3bdrEFVdcwZQpUwopg9wUmu0dU8Xvv/8+y5cvZ9myZYwfP579+/eXSk6wIwS/wzmY3V/F5G/1FGfeeCN/fvopf7zxBn2s1VGFxpUneU8xderUvKfdYcOGMXXqVOrWrUt6ejoDBw4kODiYbt268eabb5Zpu+4myHEVk6DsZIqKIpybXvPrr78+Ka/CsWPH8iK59uzZk1tvvRWA8ePHM2PGDMCk0Ny8ebN/51S2uIeqkpKQQINzzy3TYHauElK9Om1uuYVVr7/O3hUrqNu5c7nLYKnY7N+/n6VLlzJ9+nQArr76anr16kXv3r3p27dvofn6Xbt2ceWVV3LppZeybt06zj33XObNm8eYMWNo164dTz/9NAcOHCAyMpK77rqLhx56iGeffZYnn3ySd95556SsaT179iwyH/Irr7xCv379StynmJgYkpKS8raTk5Np2LBhofMiIiJo1KgRv/3220kKIXcNwZnExER++uknFi1aRJUqVejduzcZGRklljEXqxD8iIMbN3J01y7a33WX12Q4Y/hwNn70EX+MH0/fKVNsvgRLmeJPKTRdxdVUm76QQtOuIfgRyfHxJpjd+ed7TYagKlVoO3Ike37/nd3lPGVlqfhMnTqVb775htjY2LyyYcMGFi5cSIci/G5WrVpFv379CqXQbN++PU888QQPP/wwI0aMIDo6usxTaAJce+21dO/enU2bNhETE8PkyZPzjg0ePJgdO3aclGqzdevWDBs27KQRgDO5KTRff/11Zs2aVWy7gwYNIisriw4dOvDEE0/4dQpNSwlJiY+nbseOHglm5w4thg1jw5Qp/DFuHFFnn21HCZYyw99SaE6dOrXYY84pNE+VahMqcQpNS8k4smMHBzdupKMP2F4HhoTQbtQolo4Zw47584kugz+WxVISbArNssVOGfkJ5RHMzh2aDR1KtUaN+GP8eDQnx9viWCyWMsAqBD8hJSGBGs2aUSM21tuiABAQHEz70aNJ3bSJpHnzvC2OxWIpA0qkEESkqogElrUwlqI5cfgwu3//vdyd0U5Hk8GDiWjenNUTJpDjA16jFouldLikEEQkQESGi8i3IrIH2AjsFJF1IvKyiLT0rJiVmx0LFqBZWUT7mEIICAyk/V13cXjLFrbNmeNtcSxlQFFOVBb/xd3v09URQgLQHHgUqK+qjVS1HtATWAy8KCLXu9WyxWWSExIIq12b2uUQ7tpdGvXrR83WrVnz1ltknzjhbXEspSAsLIz9+/dbpVBBUFX2799PWFiYy9e4amXUT1Uzi2jwAPA18LWIBBe+zFJask+cYOeCBTQeOJCAQN+bpZOAADrccw/z//UvtsycSUtH8DGL/xETE0NycjJ7/SyHdkZGhls3PV/EU30ICwsjJibG5fNdUghFKYOSnAMgIo2Aj4D6QA4wSVXHiUgt4HMgFtgGDFPVg67UWZHZ8/vvZKan+9x0kTMNe/akTseOrH3rLZoMGkRIjRreFslSAoKDg2laTiHVy5LExEQ6derkbTFKha/0we1FZRF5uJRtZgEPqmpr4BxgtIi0AR4BflbVlsDPju1KT3JCAoHh4dTv3r3klWRkQNeucN99ZSaXMyJC3KOPknHgAMuef94jbVgsFs9zWoUgIl84lS+B20rToKruVNUVjvdpwAYgGhgCfOg47UNgaGnaqQicFMyuNMPJ7dth2TL46CPIdGkg5za127Wj3ahRbPvmG/754QePtGGxWDyLKyOEw6o6zFGuAkoXHNwJEYkFOgFLgChV3QlGaQD1yqodf+Xghg0c3bWr9NNFb74JgYFw8CBMmVI2whVB29tvp1bbtix9+mmO+dk8tMViATmdRYGINFXVrSJSR1X3iUgtx2Jy6RoWqQbMB55T1ekikqqqkU7HD6pqzSKuGwmMBIiKiopzTrHnDunp6WUW4MpTHJ4zh/S5c4l64QUCq1cvdNyVPgSmp9N92DD29ehB+I4dhO3Zw5JPPyUnJMQjMmfu2sXeF18k9IwzqPWvf50yzpE/fAenw/bB+/i7/FD+fejTp89yVe1S6ICqulSA2a6e60JdwcAPwANO+zYBDRzvGwCbTldPXFyclpSEhIQSX1tefHvZZfrj9dcXe9ylPrz2miqoLlumGh9v3r/2WtkJWQQbP/lEP23TRjd//vkpz/OH7+B02D54H3+XX7X8+wAs0yLuqe4sKpdJSEsxj4yTgQ2q+prTodnACMf7EUDxsV8rAekpKaRu2lS62EXZ2TB+PJx3HsTFQZ8+0K8fPP88FJEIpKw449priTrnHFb873+k/fOPx9qxWCxlizsKoay8VXoANwAXiMgqRxkMvAj0F5HNQH/HdqUlxREGOLo0CuGbb2DbtpOti557Dvbtg7FjSyHdqZGAAM7573+RoCAWP/aYDWthsfgJ5T5CUNVfVVVUtYOqdnSU71R1v6r2VdWWjtdSr1P4M8nx8SaYXZMmJa9k7Fho3BiGDMnf160bDB0Kr7wCBzz3EVdt0IAujz/O3pUr2fD++x5rx2KxlB3uKATvZd2uZJw4fJg9y5aVbrpo1SqYPx/uvhuCCvgfPvusmTL63/9KJefpiL34YhoNGMCaCRM4uHGjR9uyWCylx2WFoKprPSmIJZ8dv/xS+mB248ZBlSpw662Fj7VrB8OHm/WFnTtL3sZpEBG6PvkkIRERLHr0URvryGLxcdzyVBaRq0SkuuP9f0Rkuoh09oxolZfcYHZ1ShrMbs8e+OwzuOkmqFnIctcwZoxxUnvuuZKK6RJhNWty9rPPkvrnn6x+4w2PtmWxWEqHu6ErnlDVNBE5DxiI8Sh+u+zFqrxknzjBjgULiO7dGwkoYf6iiRPhxAm4557iz2nRwoweJk0yC88eJLpXL5pfeSUbpkxhz/LlHm3LYrGUHHfvOLnmIhcBb6vqLMAzHk6VlN1Ll5J15EjJ1w+OH4e334YLL4RWrU597hNPQEAAPP10ydpyg84PPUS1mBgWPfYYmUeOeLw9i8XiPu4qhBQReQcYBnwnIqElqMNyClIcweyizjmnZBV88QXs2gX33nv6c6OjYfRoE+Now4aSteciwVWr0v355zmSksKKl17yaFsWi6VkuHszH4bxMB6kqqlALeDfZS1UZUVVSUlMpEGPHiULZqdqTE1bt4YBA1y75pFHzOLzk0+6356b1O3cmTa33MLfX39NckKCx9uzWCzu4ZZCUNWjqjpdVTc7tneq6o+eEa3ycXD9eo7u2kVM794lq+C332DFCrN2cIoYQidRty488AB89ZW51sO0v+suIs84g6VPPUW2B72lLRaL+9jpHh8iOT4eCQigYUkVwrhxxqrohhvcu+6BB6BWLfjPf0rWrhsEhoTQ/cUXOXH4MIemTbPpGi0WH8IqBB8iOT6eOp06EVacqeip2L4dpk+H22+HqlXduzYiwkwdff89LFjgfttuUrNVKzrccw8Zq1ax7p13PN6exWJxDasQfIT05GRS//yz5NZFb75ppolGjy7Z9aNHQ/368NhjZi3Cw7S++WbCu3Vj9RtvsGXGDI+3Z7FYTo+7jmldRGSGiKwQkdUiskZEVntKuMpEbjC7mJJ4Jx85Au++C5dfbmIXlYQqVYwZ6q+/QjlkPBMRIq+7jvrdu7NkzBh2/Pqrx9u0WCynxt0RwqfAFOAK4BLgYserpZQkx8cT0bw51UsSzO6jjyA1tfQ5k2+7DWJj4fHHy2WUIEFB9Bw7lojmzfn1vvs4sH69x9u0WCzF465C2Kuqs1V1q6puzy0ekawSceLQIfYsW1ay2EU5OSYmUZcu0L176QQJCTFOaitWmPWIciC4WjV6T5xISGQkiaNGkZ6cXC7tWiyWwrirEJ4SkfdE5FoRuTy3eESySkTKggVodnbJ1g9+/BE2bjSjA1dNTU/FddcZP4b//Mck2CkHqtSrR5+JE8nOzCRx1CiOp6aWS7sWi+Vk3FUINwMdgUGYqaLcaSNLKUiJjyesTh1qt2/v/sXjxkGDBnDVVWUjTGCgCY+9cSN88knZ1OkCES1a0GvCBNJTUpg/ejRZGRnl1rbFYjG4qxDOUtUuqjpCVW92lFs8IlkloTTB7Kr88w/MnQt33mmme8qKyy83U1APPwy7d5ddvaehXlwc5774Ivv++IOFDz1kM61ZLOWMuwphsYi08YgklZTdS5eSdfRoiaaLor/+GkJD4Y47ylYoEXj/fTh0CEaMMOsU5UTjgQPp/PDDJP/8M8tfeME6rlks5Yi7CuE8YJWIbLJmp2VDSnw8geHh1Hc3mN3Bg9T/8Ucz51+3btkL1r49vP66MUF97bWyr/8UnHnDDZx5001snjqVDZMnl2vbFktlJuj0p5zEII9IUUlRVZIdwewCQ0Pdu/i99wjMyHAtqmlJueMO+OknePRR6NULunb1XFsF6PTggxzbvZtVr79OeP36NL3YLlVZLJ7GpRGCiDFfcTY1LWh2mnuOxXUOrFvHsd273XdGy8qCCRM42KkTlDSrmiuIGIe3hg3hmmvg8GHPtVWw6YAAznn+eep17cqSxx9n1+LF5da2xVJZcXXKKEFE7haRk9xgRSRERC4QkQ+BEWUvXsUmL5hdr17uXThzJvzzD8mXl4PFb82aJh3n9u0walS5OKzlEhgSwvnjx1O9aVN+ueceDm7cWG5tWyyVEVcVwiBMtrSpIrJDRNaLyBZgM3At8LqqfuAhGSssKQkJ1O3c2f1gdmPHQtOm7C+tI5qr9OhhHNamToUPPiifNh2E1KhB77ffJrhaNRJHjSLtn3/KtX2LpTLhkkJQ1QxVfUtVewBNgL5AZ1Vtoqq3q+oqTwpZEckNZue2d/KyZSbvwT33GJ+B8uKRR6BPH7jrLo9nVytI1QYNjOPaiRPMu/56G+LCYvEQbkc7VdVMR2KcVA/IU2nIzRjmtrnpuHFQrRrcfLMHpDoFgYHGUa1KFbOeUM6OY5FnnEH/Tz4hICSEn0aMYNeiReXavsVSGbDhr71ESm4wO3eik+7cCZ9/DrfcYnIYlDcNG8KHH8Lq1fDv8s+cGtGsGQM++4xqMTEkjhrFtm+/LXcZLJaKTLkrBBF5X0T2iMhap31jRCRFRFY5yuDylqs8OZ6ayp7ly4l2d3Tw9tvGwujuuz0jmCsMHgz33w8TJpjF7XKmSr169PvwQ+p07MjChx5i40cflbsMFktFxSU/BBF54FTHVdUdz6UPgAlAwX/y66r6ihv1+C07ShLMLiMDJk6Eiy+GFi08J5wrvPAC/PKLGanExUGjRuXafEiNGvSZNImFDz3Eipde4tjevXR84AGs5bPFUjpcHSFUd5QuwL+AaEcZBbgVykJVfwEOuHNNRSMlIYHwunWp3a6d6xdNnQp793rWEc1VQkNh2jTIzIThw82opZwJDA2lx2uv0fLqq9nw/vssfuwxcjIzy10Oi6UiIe7EihGRH4ErVDXNsV0d+FJV3fJgFpFYYI6qtnNsjwFuAg4Dy4AHVfVgMdeOBEYCREVFxU2bNs2dpvNIT0+nWrVqJbq2NGhmJrsefpjwLl2IHD7cxYuULrffDjk5LJs8OS/Mtbf6kEvUvHm0fv55tt14I9tKsMhdFvKrKulz55I2Zw6hbdpQ87bbCHDX67sUePs7KAv8vQ/+Lj+Ufx/69OmzXFW7FDqgqi4XYCMQ6rQdCmx0pw7HdbHAWqftKCAQM2J5DnjflXri4uK0pCQkJJT42tKQ8ssv+mmbNpqcmOj6RQkJqqD67rsFdieUqWwlYsQI1YAAI6OblKX8m7/4Qj9r107nXn21HjtwoMzqPR0+8R2UEn/vg7/Lr1r+fQCWaRH3VHcXlT8GljoWgZ8CllB4LcBtVHW3qmarag7wLtCttHX6KskJCQS5G8xu3DioXdsEsvM1JkwwaxrXXQf79nlNjBZXXUXPceNI/fNP5l1/PekpKV6TxWLxV9xSCKr6HCZJzkEgFbhZVZ8vrRAi0sBp8zJgbXHn+jOak0NKQoJ7wey2bIFZs0ygufBwzwpYEqpVM+sJ+/aZUNleWE/IJeaCC7jgvfc4fuAAPw4fbkNdWCxu4pZCcASwawNEqOo4YL+IuPU0LyJTgUVAKxFJFpFbgf85hdLuA9zvTp3+woF16zi2Z4975qYTJhinsDvv9JxgpaVTJxNO47vvjFLwYmKbup070+/jjwkICuKnESPYsWCB12SxWPwNd6eM3gK6Y+IXAaQBb7pTgapeq6oNVDVYVWNUdbKq3qCq7VW1g6peqqo73ZTLL0hOSEACA4l2NZhdWhpMnmzSY0ZHe1a40vKvf8Hzz5tAeLfdVq5JdQoS2aIF/T/5hKoNG5I4ahQrX33VWiBZLC7grkI4W1VHAxkAaiyByjB3Y8UmN5hdaGSkaxd88IEJOX3ffR6Uqgx59FEYM8bIPWqUV5VC1QYNGDB1Ki2GDWPD++8zb8QIu65gsZwGdxVCpogEAgogInUB7/3r/Yi8YHa9e7t2QU4OjB8P55wD3fxojf3JJ41iePddE4DPiykwg8LC6PbUU5z32msc/vtvvr/ySpLmzfOaPBaLr+OuQhgPzACiROQ54Feg1IvKlYHk+HjAjWB2330Hf/3lP6ODXETguefgwQfhzTfNq5fzIjceOJALv/qK6o0bs+C++/j92WfJPn7cqzJZLL6IWyk0VfVTEVmOCX8twFBVLd9YyH5Kcnw8ES1auB7MbuxYs25QHklwyhoRePllOHHC5GUOCTHhLrwYWqJao0b0//hjVo8fz4YpU9i7ciXnvfoqNZo29ZpMFouvURIroy5AbVWdABxx18qoMnI8NZW9K1a4PjpYuxZ+/tnkHggO9qxwnkLE+E+MGgUvvWTWFrxMYEgInf7v/+j19tsc272buVddxRYvBOizWHyVcrcyqozkBrNzORnOuHHG5+D22z0rmKcRMdNGt9wCzzxjppJ8gOjzz+fC6dOp1a4dix9/nIWPPkrmkSPeFsti8TrWyqgcSI6Pdz2Y3b59JhHNDTcY72R/JyAAJk2C66+H//zHTCX5AFWiorhg8mTa3Xkn2+fMYe6wYRws50xwFouvYa2MPEz2iRPs/PVXonv3RgJc+LgnTTKhru+5x/PClReBgTBlClx9NTz0kBkB+QABgYF0GD2aCyZPJuvoUX4YPpwNH3xAjhe9rS0Wb1JSK6N61srINXYvWULW0aOueSdnZpoplv79oW1bzwtXngQFwccfm0Xy++4zyX58hKhu3bjw669pcO65rHz5ZeYOG8a+P/7wtlgWS7njbiyjT4GHgBeAnRgroy89IVhFITk+3gSzO/vs05/81VewY4dv5DzwBMHBJq/DJZfAnXdS34dSYIbVqsX5EybQc+xYjqem8uN117F0zBiOp6Z6WzSLpdxwy+wUQFU3YsJgW06D5uSQkphIg/POcy2Y3bhx0LIlXHih54XzFiEh8OWXMHQoZ77yigmO98QTZq3By4gIjfr3p/6557LmzTfZ9MknJP38M53+7/9oeumlNiObpcLjrtlpmIg8ICLTReRrEblfRMI8JZy/kxvMziVz08WLYckSs3bgAzdHjxIaCjNmsGvgQGOOOnQoHDrkbanyCK5alc4PPcSgL76gWqNGLH7sMX6++WYO/f23t0WzWDyKu3eej4C2wBuYvMitMTkSLEWQHB+PBAbS8PzzT3/yuHEQEQE33eRxuXyCsDA2PvwwvPEGfP+9Cc/hY1Y+Nc88kwGffEK3p54iddMmvr/8claNHUvWsWPeFs1i8QjuKoRWqnqrqiY4ykjgDE8IVhFIdjWYXXKymUa59VYzhVJZEDHOdz//DKmpRin4mKOYBATQYtgwLp4zhyaDB7P+3Xf5dsgQUubP97ZoFkuZ465CWCkieam+RORs4LeyFalikJ6UxKHNm4lxxRntrbdMvJ+77vK8YL7I+efD8uXQpg1cdplZU/BipNSiCKtdm+4vvEDfKVMIDA1l/p138su995J9sMjU3xaLX+K2YxqwUES2icg2TKKbXk7JbSwOXA5md/QovPMODBkClTmuTkwMzJ9vvJr/+19jieSDFj65Jqpn3XsvOxcsYPfTT7Pi5ZfJsIrBUgFwVyEMApoCvRylKTAYuBi4pGxF82+SExKIaNmSao0anfrETz+FAwf8L6qpJwgLg/feMz4K8+ZB166wbp23pSpEYEgIbUeO5KJvviG8c2c2ffQRswcMYPUbb3AiLc3b4lksJcZdhdANOKCq24EbgNcxge62O/ZZcApmd7rpIlWzmNyxI/TsWS6y+TwiJiBeQgKkp8PZZxv/DB+kWnQ0NW+8kcEzZ9LgvPNYO3EiswcMYN2775J19Ki3xbNY3MZdhfCEqqaJyHnAQOBDwHdcTn2EHb/84lowu59/Nk/A993n1dDQPkmPHmZdoX17k0L00Ue9mqv5VEQ0b07P119n0JdfUqdjR/4YO5bZgwax6dNPyT5xwtviWSwu465CyP1HXgS8raqzsMHtCuFyMLuxY6FePbjmmnKRy+9o2BASE2HkSHjxRbjoIti929tSFUutNm3o/fbb9P/4Y2o0a8by55/nm8GD+fvrr218JItf4K5CSBGRd4BhwHciElqCOio02cePm2B2ffqcOpjd5s3w7bcmOb0rXsyVldBQs+g+aZKZRmrd2gTK83IWtlNRt3Nn+k6ZQp933yWsdm2WPPkk3156Kdu++w71Mespi8UZd2/mw4AfgEGqmgrUAv5d1kL5M7uWLCHr2LHTWxeNH29i+4waVT6C+Tu33w5//GGC/t1yiwkA+Ndf3paqWESEBueey8Bp0zj/jTcICAlh4b//zXeXX86WmTPtVJLFJ3E3uN1RVZ2uqpsd2ztV9UfPiOafpCQkEFSlClGnCmaXmmqecq+9FurXLzfZ/J4zzzSmqRMnwu+/m/WFl14yUWJ9FBEh5oILGDx9Ouf+738ALH78cWb178/ad96xwfMsPoWd7ilDNCeHlIQEE8wu5BRLK++/D0eOVNyopp4kIADuuAPWrzdBAB95xHg4L1/ubclOiQQEEHvRRQyeMYM+775LzVatWD1+PDP79uX3Z57h8LZt3hbRYrEKoSzZv3Ytx/buPbW5aXa2id/Tsyd07lx+wlU0oqNh+nT4+muz0NytGzz4oFG0PkzuVFKfSZMYPHMmsRddxN/TpzPnoouYP3o0u5cuRX14fcRSsXEp/LWIPHCq46r6WtmI49+kJCScPpjd7NmwbRu88kq5yVWhufxyuOACM1J47TWjJCZOhIEDvS3ZaYls2ZKzn3mGDvfcw+bPP2fz1Kn8fPPN1GzdmjNHjKDxwIGnHmlaLGWMqyOE6o7SBfgXEO0oo4A27jQoIu+LyB4RWeu0r5aIzBORzY7Xmu7U6Sskx8efPpjduHHQpIkJVWEpGyIjjRL45RdjlTRokMlJvW+ftyVzifA6degwejRDfvqJbmPGkH38OIseeYTZgwax7t13OeYn/bD4Py4pBFV9WlWfBuoAnVX1QVV9EIgDYtxs8wNMCAxnHgF+VtWWwM+Obb8i7Z9/OPTXX6e2Llq1yiyK3nWXSSlpKVt69jSf8RNPwOefm0Xot98GP7HoCQoLo8VVV3HRrFn0njiRiGbN+GPsWGb27csv99xDSmKi9WeweBR31xAaA87/rhNArDsVqOovwIECu4dgvJ5xvA51Uy6vk5KQAJwmmN24cVC1qglzbfEMYWHwzDOwYoXxWbjzTqMYPvrIZz2dCyIBATTs2ZML3nuPi2bP5swbbmDfqlXMHz2aWf37s2rsWNK220gxlrJH3FnAEpHHMb4IMwAFLgM+V9UX3GpUJBaYo6rtHNupqhrpdPygqhY5bSQiI4GRAFFRUXHTpk1zp+k80tPTqVaGuQf2jR1LzpEj1Hv88SKPBx84QPdrrmHn4MFsLqNAdmXdh/LG4/KrUmvpUppOnkz1zZs50qQJW2++mX3nn19moULK6zvQ7Gwy1qzh6KJFHF+3DlQJadmSKt27E9apEwGlWGuwvyPvU9596NOnz3JV7VLogKq6VYDOwL2O0snd6x11xAJrnbZTCxw/6Eo9cXFxWlISEhJKfG1BMg4e1M/atdNV48YVf9LTT6uC6saNZdZuWfbBG5Sb/NnZql9+qXrmmeY7iItTnTtXNSen1FV74zs4smuXrn3nHZ01cKB+2qaNfnH22brk6ad135o1mlOCPtnfkfcp7z4Ay7SIe6rbZqequkJVxznKyhIqqILsFpEGAI7XPWVUb7mQMn8+mpNTvLnp8eMmCc6FF0KrVuUrnMX4Llx5JaxZYxwC9+0zC8+9esGCBd6Wzm2qREXRduRILvnuO/pOmUJ0795snTmTH66+mu8vv5z1kyeTnpLibTEtfohbCkFEQkVkuIg8JiJP5pYykGM2MMLxfgQwqwzqLDdS4uMJr1ePWm3bFn3CF18YW3mb88C7BAWZnNWbNsGECSae1PnnG0Xt445tRSEBAUR168a5L77IZYmJdH3iCQJDQ1n12mvMHjCAH665hg0ffMCRnTu9LarFT3B3hDALswCcBRxxKi4jIlMxmdZaiUiyiNwKvAj0F5HNQH/Htl+Qffw4O3/7rfhgdqomqmnr1ib+jsX7hIbC6NHw998m9MXSpdClS/4owg8JqVGDltdcw8Bp07j0hx/oeP/95GRns/Lll5nVrx8/Xncdmz75hKN7/GrwbSln3LV9jFHVgiajbqGq1xZzqG9p6vUWuxYvNsHsipsu+u03Y/EycaJXcx48+ihcfLFJM2BxUKUKPPSQCYXx2mumfP21cXS7914Tbjsw0NtSuk21mBja3HYbbW67jcPbt/PP3Ln8M3cuy194geUvvki9uDgaDxpEo/79Ca9Tx9viWnwId0cIC0WkvUck8VNSEhIIqlq1+GB2Y8dCzZrGUcpLfPaZSSfQrx8sW+Y1MXyXiAh4+mnYvt18UH/+aRwHW7UypsKHD3tbwhJTo0kT2t1xB4NnzOCib76h/Z13cvzgQZb997/M7NOHn2+5hSOJiRzZscPbolp8AHcVwnnAchHZJCKrRWSNiKz2hGD+gObkkJyQQMPigtlt3w4zZpgEL1WqlL+AwJw5Ztr8rLMgKspMmc+Y4RVRfJ9ateDhh2HrVrPuExVl1n1iYszr3397W8JSEdGsGe3vvJOLZs9m8MyZtB05kmN79nDoyy+Z1b8/3w4dyqqxY9m7ciU5fuKzYSlb3J0yutAjUvgp+9euJWPfvuJTZU6YYKaJRo8uX8Ec/PADXHGFUQY//QQZGebB94or4H//M7HgbObOIggKMmk7r7rKDKnGjTNWYuPHm3m3++6DPn38+sOLbNmSyJYt6XD33fz0+edEHzvGjvnz2fD++6x/911Ca9akYc+eRPfuTYMePQj2czt/i2u4pRBU1bpHOpEcH2+C2fXsWfhgejq89565+zZqVO6yxcfD0KHQpo1RDBERpiQkwIgR8O9/m/wyb7xh8vRYiqFLF/j4Y6NB337brAX17WtyMdxzDwFe+G7LmqCoKFr37k3rm27ixOHD7PztN1ISE0lJTGTr7NlIUBD1unQhulcvonv1onqTJt4W2eIhXI12+quqniciaRgPZXF+VdUaHpTRZ0mJj6deXFzRwew++sgkwvGCqemCBXDJJdCiBcybZ2ZCcgkPh2nTzLEXXsifHYmIKHcx/YsGDUxIjMceg6lTzajh9tvpXq0aXHcd3HgjdO/u16MGMNZKTS68kCYXXkhOVhb7/viDlMREdsyfz4qXXmLFSy9RvUkT6p9zDvW7dyeqWzdC7I+nwuDqCOFrEekK1FRVG10LSNu+nUN//03zq64qfDAnx0wvdO0K55xTrnItWgSDB0PjxmaaqCgjkoAAeP55oxTuuMNYHs2ZA7Gx5SqqfxIWBjffbBZmfvmF/c89R/2PPzZ5n5s3N4rhhhugaVNvS1pqAoKCqBcXR724ODo9+CBp//zDjl9+YefChWydPZvNn3+OBARQs00b6nfvTv1zzqFup04E2hzhfourCiEaGAec6VhEXgj8BixS1YKB6ioFyacKZvfDD8b56dNPy/WJcdky44Bbvz78/LNZEz0Vt9xiInFfcQWcfTZ8843JM2NxARHo1YuNqtTv0sWYq370EYwZA089ZSKv3nijWYeoIE/Q1Rs3ptX119Pq+uvJPnGC/WvWsGvxYnYvWpS39hAYGkrduLi8EUTNM88s2j/H4pO4pBBU9f8ARCQEkxPhXOAW4F1HYDq3ciJUBFISEog84wyqRUcXPjhunJliuPLKcpNn1SoYMABq1zbrBw0bunZd375mVHHRRSaSwyefGAVhcYNq1czCzIgR8M8/5kP86CO4/Xa4+26zmHPjjcYxsYKEPQ8MCckbPTB6NJlHjrBn2TJ2LVrErsWLWfWayZkVGhlJvS5dqOs4N7JVKwIqyGdQEXH3mwkHagARjrID8E/XzlKQcfAge1esoM3ttxc+uGGDGSE8+yyUU7artWuNj0G1akYZuLvO2bo1LF5s7ltXXmmcd//9b7+fDvcOjRubdYZHH4XffzeKYepUs3ATFQXXXmuyvJ17rl86vRVHcNWqeYvOAMf27mXXkiXsXrSI3b//TtJPPwEQVKUKdTp2pG7nztSLi6N2hw4EhYV5U3SLE64uKk8C2gJpwBLMlNFrqnrQg7L5LDtyg9kVNV00frwJjXDHHeUiy8aN5ik/NNQog5KuA9SrZ6aZbr7ZmOJv3mwsLa0FUgkRMfNv3boZD+hvv4UPPzQf6tix5gO/9FK47LL8L7ACEV63Lk0vvpimF18MwNFdu9izYgV7ly9n74oVrHnzTVAlICiImm3bUq9zZ+rGxVG3U6dTZxy0eBRXRwiNgVBgM5ACJAOpHpLJ50lOSCA8KqpwMLsDB8yf/rrroG5dz8uRHM7w4ebe8/PPZpG4NISHG6/m5s3NovP69TBpEhQXs8/iIiEh5sZ/2WXG6/m772DmTJPV7b33zNBu8GBzfPBgqFHxjPaq1K9P7ODBxA4eDMCJQ4fYu2oVe5cvZ8/y5Wz6+GM2TJkCQETz5tQ+6yxqt29PnfbtiWjZ0k4zlROuriEMEhHBjBLOBR4E2onIAczC8lMelNGnyMrIYOdvv9Hs0kuRgnMq770Hx46ZODgeZutWeOCBs1CFxESTFKwsCAiA554zSuDuu6FjR7j/fnjySXPfspSSGjXgmmtMOX7caPKZM2HWLGP/GxJiRgyXXWZGEKezDPBTQiIiTppiysrI4MDatexZvpy9K1eSEh/PlunTAQgMD6dWmzZGQXToQO327anSoEHh/5+l1Lisdh1JFdaKSCpwyFEuBroBlUYh7F6yhOxjx4guOF2UlWU8k/v0gQ4dPCrDP/+Y+GsZGYEsWOCZJ/jhw80i9SOPwMsvmynwcePMOoP9H5YRoaFmRDB4sHF6W7TIKIfccCd33GF8Gy680HwZcXEVat3BmaCwMOp16UK9LiaJl6qSnpTE/tWr2bdmDfvXrOHPTz9l4wcfABBWuza1O3SgTocOHM/J4XjHjnaqqQxwdQ3hHszIoAeQicPkFHifSraonBwfb4LZFbTPnDEDkpKMUvAgKSlGGRw8CC+/vJqzzorzWFt16phBz803w7/+ZdZCL7rILJM0a+axZisngYFw3nmmvPyyCcM9cybMnm2GZ088YYIk9utnlEP//sZmuIIiIlRv3JjqjRsT61iHyD5xgtRNm9i3ejX7165l/+rVebnMv37jDao2bEjN1q2p2bo1tVq3pmabNoTXrWtHEm7g6gghFvgKuF9VK222Dc3JISUxkYY9exYOZjdunLlLXnSRx9rftcvMJuzebTyQMzLSPNaWMz16mPwxb7xhTOzbtoXHHzeWSBVsLdQ3EDGjzA4djDLYu9dMLf34oylffmnOa9XKKIcBA4zNcPXq3pXbwwSGhFC7fXtqt88PuHzi0CHiP/uMxsHBHNi4kYPr15P88895x8Nq1z5ZSbRuTbVGjaySKAZX1xAe8LQg/sD+NWuKDma3bJnJezB2rMeG9Hv3mofDpCRj1XrOOWbtoLwIDoYHHoBhw8yawhNPmBA/b71llJTFg9Stm7/uoGpMm3OVw3vvGU0dFGRMWQcMgN69TQymSqCtQyIiCGvdmja9e+fty0xP5+CmTRzcsIED69dzcMMGdi1ahDoiuAZXr54X3C/yjDOIcLwPqYCL+e5il+7dIDeYXXTBYHbjxpmns5tv9ki7Bw6YGYK//zYGKued55FmXCImxjygzp0Ld91llNS118KrrxpfPIuHETERC9u0MXGyjh83DyO5CuI//zHnhYUZk9eePU3M8+7dK/wIIpfgatXyneYcZB8/TurmzRxcv54DGzdy6M8/2fbtt2R+/nneOVXq1yfyjDOIbNmSCMdrjaZNiw5tX0GxCsENkuPjqdely8nBvHbuNOaDd97pEXPB1FTz0LdhgwktUVyk7fJm0CAzzf3SSyZI3rffmqgNI0dC1arelq4SERpqFpUuuMAk99m3D379FX75xUQ5fPFFYzYWEACdOhkFkVvKwTTaVwgMDaV2u3bUbtcub5+qcnTnTlI3byb1zz9J3byZQ5s3s2vhQnKyTMg2CQqiRmwsEc2bU6NZs7zXGrGxFTJmk1UILnJ4+3YOb9lCy6uvPvnA228bC6O77y77Ng8bA5PVq82a9YABZd5EqQgPN0rguuvMaOGBB8y9Z/Ros12J7je+Q506xhRs6FCznZZm3NBzFcTEiWZqE4ytcs+e1K9VyzjKtWpVYa2YikJEqNqwIVUbNswzfwWzeJ22fXuekkj9808OrF/PPz/+aKbsAAkIoGp0tFESzZpRo3lz89qsGSF+PBKzCsFFUooKZpeRYf5gF19svLnKkPR0sz69bJmZovHgWnWpadnSTCH99psxkHnmGZM+4OabjZIorcOcpRRUr27mG/v3N9vHjxsLgQULjJL44gvOPHTIDPWqVzdrD7ke1t26QXR0pbMzDgwJyVtjcCYrI4O0bds4tGULhx3l0JYtZkSRmZl3XnjdutRo2pTqsbFUb9KE6k2aUCM2lqrR0T4//WQVgoskx8cT2aoVVZ2jxk2dalZ7yzjnwdGjxidp4UJj/5/7sOfLiORbTW7YYNYUJk82+vKKK4xFko2k6gOEhprF53PPNTFKcnJY+skndANYsgSWLjWhNnJvcA0anKwgunSBSmrvHxQWRs0zz6RmAS/QnKws0lNSOPz333lK4vC2bST9+CPHU1PzzpPAQKo2bEj12FhqOBRF9dhYasTGojk55dyborEKwQUyDh5k38qVtB05Mn+nqhl6t2tXphP7GRnGSTUx0VjxFJVuwddp3doYvzz7rPFZePtt+Oors7b50ENmGsxGRPYRAgI42rixsUy68UazLyMD/vjDKIfcMmtW/jUtWpj1iI4d80uDBpVuJJFLQFAQNZo0oUaTJmYtx4njqamkbd9O2vbtHN62Le/93mXLyDp2LP/EoCDmNGpEtUaNqBYTY16dtssrAKBVCC6wIzGxcDC7+fPN5P5775XZH+HECRNt9Mcf4f33zdy8P9OggVlwfuwx8zG9/rqZXWvTxowYoqMr5w3E5wkLMwkyzj47f9/Bg2b+cskSWLnSTDvl+kOAWTByVhAdO8IZZ1SYcN8lJTQyktDISOqcddZJ+1WVY3v3GgWxbRvrfv2VCBHSk5LYs3w5WUeOnHR+eN26JymIao0bE9OnD8FlbMFRub8tF0lOSKBK/frUbOOU9mHsWJN8YPjwMmkjMxOuvtpY60yc6DELVq9QvbrxXbjrLmOQ9fLLpn+1ap3DzTcbxdexY6V9wPQPatY8eS0C4NAh81C0alV+GTfOPNmAUSzt2xsHu7Zt80vDhpX+yxYRqtSrR5V69Yjq2pXkunXp6fClUFWOp6aSnpRE+j//kJ6cTHpSEmlJSexavJhju3cDcMVvv5W5XFYhnIasjAx2LlxIsyFD8r0bt2wxIQUee8yY2pS2jSy4/noTqWD8+HKLnF3uBAebfl53nRkFPfNMGuPHh/Lqq2bUcN11Rr/aVJ5+QkREvglrLpmZJia7s5KYPdssKOUSGWm+cGcl0batSfVXyRUFGGURVrMmYTVrUqeIuGhZGRkcSUnxSOwmn1IIIrINk3MhG8hS1S7elQh2L15M9rFjJ08XvfGGMc+7885S15+dbdLzfvEFvPKKR6xXfQ4RGDgQQkPX0r59b7780mQbffxxU3r0MMph2DAzCLP4EcHBZlTQvr3JLZ3L3r2wbt3JZfp0ePfd/HNq1sxXDq1amSmnVq3ME0Iln3pyJigsjIgytmrMq9sjtZaOPqq6z9tC5JIbzK5e165mx+HD5mln2DDX81QWQ06OybL46afGfv/BB8tAYD+jdm0YNcqUbduM4dYnnxhde889ZgH6uuvgkkugShVvS2spMXXrmoVrpxATqMKePYUVxZdfGvf8XIKDjVl3q1YnK4ozzrDOLmWMLyoEn6HIYHYffGCcfUqZ80DVOHBNmWLilz32WOnl9XdiY03myUceMUYun35qEvZ8843JxXD55cYct1+/CpO3vnIjYvI9REUVss5h/37YtMmUP//Mf//99/lrFACRkXRu0MBYPTVvfnKxU1BuI+rwvPMFRGQrcBBQ4B1VnVTEOSOBkQBRUVFx06ZNK1Fb6enpVDtNxpcTW7aw79VXibzpJqp07Qo5OZx9ww2cqFmTlaUIc60KEya0YPr0GIYP385tt20t0e/WlT74Mq7In50Nf/wRyU8/RfHLL3U5ciSIwMAc2rU7zNln76dbtwM0a3bEa/97f/8OwM/6kJ1N2O7dVElKIjwpiSpJSYRs30613bsJ27MHcbLnzw4L41iDBmQ0bMgxp5IRHU1GvXqoD+WHLe/voE+fPsuLnJJXVZ8pQEPHaz3gD+D8U50fFxenJSUhIeG056x89VX9rEMHPX7okNkxe7YqqE6bVuJ2c3JU/+//TDX332+2S4orffBl3JX/xAnV+fNVH3lE9ayzzGcIqtHRqrfeqvrVV6qpqR4RtVj8/TtQ9f8+5Ml//Ljqn3+qfv+96oQJ5g926aWqbduqhofn/2BAVUQ1Jkb1vPNUr7tO9fHHVd99V3XePNXNm1UzMrzTh3ICWKZF3FN9aspIVXc4XveIyAxMNrZfvCVPckKCCWaXG7Ru3DgT7vPyy0tc5xNPmMXj0aONN68d0bpOcLBxbjv/fOPfsGOHCZnx/ffG8W3y5Pwo0BdeaEr79tYJrtIQEmLiqBQIOQEYNbBzpwkZ/PffsH27WbTats3EXJk2zQxHcxExjjSxsaY0agSNG5vX3FKrVoX7A/uMQhCRqkCAqqY53g8AnvGWPHnB7K65xuxYu9YkKXnhBXNnKgHPPmsWj2+7zZiXVrDfUrnTsCHccospmZkmhtv335vy6KOm1KplrJbOO8+8VpI0AZaCiJgfTMOGJ5vJ5pKVZdIROiuK3LJokVnodopXBBiTc2cF4VxiYkwcqMhIv/qj+4xCAKKAGQ5b/yDgM1Wd6y1hUuLjAYjJDUsxbpz5Adx+e4nq+9//zOLxjTfCO+/Yp9ayJjg43yT++efNw+APP5gYbr/+ahamwSiDrl3z4y6de66xdrRUcoKCTErSJk3MELQgOTkmVWFSUtFl3jzzoysYkyg8PF8RRUebUvB9Ka0VyxKfUQiqugU467QnlhMnBbPbt8/YQt54Y4kM48eONXHErr3WhKSwysDzNGhg/Dtuusls79ljggX++qspr7xiUgWAMXvPVRBnn20MVOx3ZDmJgADzo8oN9lcUmZlGKSQlQXKymdNMSTFlxw74/XfjfZqRUejSHjVqmFFFgwbGOqq4Vw+POHxGIfgSGQcOsG/VKtrmugxPmmS+xHvucbuut94yYRuuuAI++qhShZv3KerVOzlNwNGjJmbbb78ZBTF1qhm5gQm10amTKZ07m3LmmdY3ynIagoPNOkPjxsWfo2qyXuUqCcfrnt9/Jzow0CiUhQvNaxGKg9BQoxjq1zce4PXqlWkX7E+8CFLmz88PZpeZCW++aYzf27Z1q5733jOLx5dcYuzp7Q3Fd6hS5WQ/qexss0y0bBmsWGHit02aBLkBKcPC4KyzTlYSTsm3LBbXEDFzlDVrnvQD2pyYSHRBp73Dh41i2LUr/9X5vQcS8dhbVBGk5Aaza93aWB/s2GHuDm7w0UcmneSgQWY9ysfzYlR6AgPNDf+ss+DWW82+7GzjC7ViRb6S+OwzE3wQjIJv3LgLXbuaZ4XcVMctWpTY7sBiMYgY78uICDM8LSesQihAXjC7oUNNMLtx44wZ24UXulzH55+baJ4XXGDCtVirFv8kMDD/Jn/99WZfTg5s3ZqvJBITj7N0aTWccrUTFGSiKuRemxvHrWVL+1uw+DZWIRRg16JFJphdnz7GjnHJEhPMzsVVxunTTeyd884zOUXKIBiqxYcICMiPjHDVVZCYuIbevXtz5IgZTaxfb8q6dSbQ5/Tp+YYngYHQrJlRDC1anFxiY+2owuJ9rEIoQEpCAsHVqplgdiNGmCFbrqnKafjmG5PToFs3mDMHyjh3hcWHqVo1f23BmWPHTCieXEWxYYPxi5o/H5xzoAQGGqVQUFE0b2722wcLS3lgFYITOdnZpCQm0qBnTwL37DGT//feayKrnYa5c022s06djGOUB9Z7LH5IeHj+2oQzqsas/a+/CpdFi8x6ojP16uWbycfGFn6f60xvsZQGqxCc2L96NRn795vporfeMv/au+467XU//2zMGdu0Mc5QNhKn5XSI5FsPnnfeycdUTbDPXAWxfXu+A+3q1Wb0WdAiMTIyX0k0amR8nnKdZXPf2xGr5XRYheBESkICEhREwy5djBPakCHQtOkpr/nlF2NWesYZxlnRer1aSosI1KljyjnnFD6em0Zg27Z8ZZGrMLZsMd7ZBw8Wvi4iorCSiI6GfftqU7WqUU5RUdYirjJjFYITyQkJRHXtSsjs2SZBx333nfL8hQth8GAzbP/pJ/MHtlg8jXMagbPPLvqco0fznWSTkwu/rl1rTNnNgnd7/vOf/Gtr1cofveQ6yDqXqCiTl6ZOHetbU9GwX6eDw9u2mWB2V19tYhp06lR0ECwHv/9uLFEbNjRTRmXsMGixlIoqVYoP/JlLVpZRCt9+u5zo6LhC/k+7dpn1jJ078x30ClKrlvnt161b+NX5fe3aptjRh29jFYKD5NxgdiEhxmbwgw+KjRmyciUMGGB+4PHx5inKYvE3goLMtFGrVmknZbYsiCqkpxvFsHOnSY+8Z0/h1w0bjPXU/v3mmqKoVs38b+rUyVcSRW3nOvPWqmWmumzIl/LBKgQHKfHx1DzzTKp+9JF5rMkNe12ANWugf39j1REfb/5QFktFRsRYzVWvbtbKTkd2tlEKzgpj//6Ty7595vXvv81rauqp64yIyFcQucoid/vgwUZs2mQW1iMizKvz+7Awv4pA7VWsQsAEs9u7ahXthg2Dp54ypQiX0g0boG9f8wOLjzdrBxaL5WQCA80zVb16rof/ysoyy3a5CuPgwfxy4EDh9ykp+e8zM5ufsu6QkJMVRG5EiBo18l+dS1H7qlevHKMUqxAwwexQJWbLFvPrGTWq0DmbNxtlEBBg1gyan/o3aLFY3CAoKF+JuIMqzJ27gA4denLokBlppKaS977ga27ZscP4ehw+DGlpxU9xORMenj9ScqVUq3bqEhLieyMXqxAw00VVoqKo+dVXZqqofv2Tjm/ZYuISZWVBYiK0auUdOS0Wy8mIQHh4dp4pbUnIyTFe44cO5SuJ3HLokClpaUWXXbvMw2LutrP3+ekIDMxXDgEB3ahTx/iKnK5UqWJeL7mk7D3YK71CyDp2zASza9oUOXLEeCY7sX27UQZHj0JCgnE+s1gsFYeAgPyn+tKSnW2UQnp68aWo41u2pFO1ahWOHDFTYUlJ5rwjR8y95+jRwm3t3WsVQpmza9EisjMyiFm50piZOgWjSUkxyiA11awZdOjgPTktFovvExiYv+7gDomJ6+ndu/j5spwcY/qbqySOHPGME2ylVwjJCQkEh4ZSb/16ePnlvP27dhllsHev8UAuGLTMYrFYyouAgPwpI09SqRVCbjC7hkBgkyYmVAVGCfTta0YIP/xQvDeoxWKxVCQqdSrx/atXc/zAAaL//hvuvhsCA9m/32TL3LrVBBHr0cPbUlosFkv5UKlHCMnx8YgIDXNy4NZbSU01HsibNpncBqfy3rRYLJaKRqVWCCnz5hF15AghI0ZwOCCSgf2NJ/LMmcYb2WKxWCoTlXbKKGv3bg4nJRF96BBHbrmbwYNNjtwvvzQRTC0Wi6WyUWlHCBkrVgDQoHNXLn6wFYsXw7RpeevKFovFUumotAohe+FCah47xgsHH2T+SvjkE5MC02KxWCorPjVlJCKDRGSTiPwlIo94qp2Mffs4sn8/1Y+H8eKK/rz/Pgwf7qnWLBaLxT/wGYUgIoHAm8CFQBvgWhHxSKCIfya9ByLMOHAt77wj3HSTJ1qxWCwW/8JnFALQDfhLVbeo6glgGlDmM/rZ2bD+gxmEncjijOceYOTIsm7BYrFY/BNfWkOIBpKctpOBQj7CIjISGAkQFRVFYmKiW41ojpIZXJuc4Pq06vw7bl7uU6Snp7vdf1/C3+UH2wdfwN/lB9/pgy8phKIigxeKUq6qk4BJAF26dNHeJfEeW9eHxMRESnStD+HvffB3+cH2wRfwd/nBd/rgS1NGyUAjp+0YYIeXZLFYLJZKhy8phN+BliLSVERCgGuA2V6WyWKxWCoNPjNlpKpZInIX8AMQCLyvquu8LJbFYrFUGnxGIQCo6nfAd96Ww2KxWCojvjRlZLFYLBYvYhWCxWKxWACrECwWi8XiwCoEi8VisQAgqoV8v/wGEdkLbC/h5XWAfWUojjfw9z74u/xg++AL+Lv8UP59aKKqdQvu9GuFUBpEZJmqdvG2HKXB3/vg7/KD7YMv4O/yg+/0wU4ZWSwWiwWwCsFisVgsDiqzQpjkbQHKAH/vg7/LD7YPvoC/yw8+0odKu4ZgsVgslpOpzCMEi8VisThhFYLFYrFYgEqoEERkkIhsEpG/ROQRb8vjCiLSSEQSRGSDiKwTkXsd+2uJyDwR2ex4reltWU+FiASKyEoRmePY9jf5I0XkKxHZ6PguuvthH+53/IbWishUEQnz9T6IyPsiskdE1jrtK1ZmEXnU8f/eJCIDvSP1yRTTh5cdv6XVIjJDRCKdjnmlD5VKIYhIIPAmcCHQBrhWRNp4VyqXyAIeVNXWwDnAaIfcjwA/q2pL4GfHti9zL7DBadvf5B8HzFXVM4GzMH3xmz6ISDRwD9BFVdthwsxfg+/34QNgUIF9Rcrs+F9cA7R1XPOW43/vbT6gcB/mAe1UtQPwJ/AoeLcPlUohAN2Av1R1i6qeAKYBQ7ws02lR1Z2qusLxPg1zI4rGyP6h47QPgaFeEdAFRCQGuAh4z2m3P8lfAzgfmAygqidUNRU/6oODICBcRIKAKpishD7dB1X9BThQYHdxMg8BpqnqcVXdCvyF+d97laL6oKo/qmqWY3MxJkskeLEPlU0hRANJTtvJjn1+g4jEAp2AJUCUqu4EozSAel4U7XSMBR4Ccpz2+ZP8zYC9wBTHtNd7IlIVP+qDqqYArwD/ADuBQ6r6I37UByeKk9lf/+O3AN873nutD5VNIUgR+/zG7lZEqgFfA/ep6mFvy+MqInIxsEdVl3tbllIQBHQG3lbVTsARfG9q5ZQ45tmHAE2BhkBVEbneu1KVOX73HxeRxzHTwp/m7iritHLpQ2VTCMlAI6ftGMyQ2ecRkWCMMvhUVac7du8WkQaO4w2APd6S7zT0AC4VkW2YaboLROQT/Ed+ML+dZFVd4tj+CqMg/KkP/YCtqrpXVTOB6cC5+FcfcilOZr/6j4vICOBi4DrNdwrzWh8qm0L4HWgpIk1FJASzcDPbyzKdFhERzNz1BlV9zenQbGCE4/0IYFZ5y+YKqvqoqsaoaizmM49X1evxE/kBVHUXkCQirRy7+gLr8aM+YKaKzhGRKo7fVF/MepQ/9SGX4mSeDVwjIqEi0hRoCSz1gnynRUQGAQ8Dl6rqUadD3uuDqlaqAgzGrOj/DTzubXlclPk8zJBxNbDKUQYDtTEWFpsdr7W8LasLfekNzHG89yv5gY7AMsf3MBOo6Yd9eBrYCKwFPgZCfb0PwFTMmkcm5un51lPJDDzu+H9vAi70tvyn6MNfmLWC3P/0RG/3wYausFgsFgtQ+aaMLBaLxVIMViFYLBaLBbAKwWKxWCwOrEKwWCwWC2AVgsVisVgcWIVgqVSISLaIrHJE+/zGOcKki9cnikgXx/vv3L2+mDo7ich7pz/zpGumiUjL0rZtsThjFYKlsnFMVTuqifZ5ABhd0opUdbCaAHel5THgDTeveRsTG8piKTOsQrBUZhbhCBomIt1EZKEjcN3CXI9kEQl3PI2vFpHPgfDci0Vkm4jUEZHYAnHu/09Exjje3yMi6x3XTysogIhUBzqo6h+O7TEi8qGI/Oio/3IR+Z+IrBGRuY4QJgALgH6OqKUWS5lgFYKlUuKIL9+X/NAlG4Hz1QSuexJ43rH/X8BRNTHrnwPi3GzqEaCT4/pRRRzvgvEadqY5JlT4EOATIEFV2wPHHPtR1RyMp+tZbspjsRSLVQiWyka4iKwC9gO1MElKACKALx1P+q9jkpOAyYHwCYCqrsaErXCH1cCnjqiiWUUcb4AJq+3M92qCz63BJLGZ69i/Boh1Om8PJmqpxVImWIVgqWwcU9WOQBMghPw1hGcxT+LtgEuAMKdrThffJYuT/0vO116EydIXBywvYornWIHzAY5D3iggU/Pjy+RgwnA7t3PsNLJZLC5jFYKlUqKqhzDpJP/PMS8fAaQ4Dt/kdOovwHUAItIO6FBEdbuBeiJSW0RCMeGMEZEAoJGqJmAWgCOBagWu3QC0KGE3zgDWlfBai6UQViFYKi2quhL4AxOS+3/ACyLyG2aaJpe3gWoishpzUy8UhtgxvfMMJovdHMx6BI56PhGRNcBK4PWCVkmquhGIcCwuu4yIRGFGOzvduc5iORU22qnF4mVE5H4gTVVd9kVwXHNYVSd7TjJLZcOOECwW7/M2jnUDN0glP8m8xVIm2BGCxWKxWAA7QrBYLBaLA6sQLBaLxQJYhWCxWCwWB1YhWCwWiwWwCsFisVgsDv4fqrbtLvQiScwAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.clf() # Clear figure\n",
"# Plot the maximum tangential winds\n",
"for mtmax, col, dpmax in zip(M_tan_max, colors,Delta_p_max):\n",
" # Plot inner core wind speed (Eqn. 15.42)\n",
" RCV_Mtan_inner = mtmax * R_inner / R0\n",
" plt.plot(R_inner, RCV_Mtan_inner, color = col, label = r\"$\\Delta P_{max} =$ %0.1f kPa\"%dpmax)\n",
" # Plot outer region wind speed (Eqn. 15.44)\n",
" RCV_Mtan_outer = mtmax * R0 / R_outer\n",
" plt.plot(R_outer, RCV_Mtan_outer, color = col)\n",
"plt.title('Tangential velocity vs. radius')\n",
"plt.grid()\n",
"plt.legend()\n",
"plt.xlabel('Radius (m)')\n",
"plt.ylabel(r'Wind speed (m s$^{-1}$)')"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.6"
}
},
"nbformat": 4,
"nbformat_minor": 4
}