Skip to content

Instantly share code, notes, and snippets.

@bungernut
Created February 2, 2022 21:06
Show Gist options
  • Select an option

  • Save bungernut/f8ea76ef78a8d4c53108a7d532ddeaf2 to your computer and use it in GitHub Desktop.

Select an option

Save bungernut/f8ea76ef78a8d4c53108a7d532ddeaf2 to your computer and use it in GitHub Desktop.
ThermosiphonCalcs.ipynb
Display the source blob
Display the rendered blob
Raw
{
"nbformat": 4,
"nbformat_minor": 0,
"metadata": {
"colab": {
"name": "ThermosiphonCalcs.ipynb",
"provenance": [],
"authorship_tag": "ABX9TyNqIm7GitZj3LQjbCSqbgBr",
"include_colab_link": true
},
"kernelspec": {
"name": "python3",
"display_name": "Python 3"
},
"language_info": {
"name": "python"
}
},
"cells": [
{
"cell_type": "markdown",
"metadata": {
"id": "view-in-github",
"colab_type": "text"
},
"source": [
"<a href=\"https://colab.research.google.com/gist/bungernut/f8ea76ef78a8d4c53108a7d532ddeaf2/thermosiphoncalcs.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"id": "h-sQck14OOk9"
},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"source": [
"NA = 6.06e23 #atoms/mol\n",
"xe_Hvap = 12.64e3 #J/mol\n",
"xe_Hlat = 21.01 #J/mol/K\n",
"xe_GasDen = 24.22 #L/mol\n",
"print(xe_Hvap/xe_Hlat)"
],
"metadata": {
"id": "NGkPUm3iOSFt",
"colab": {
"base_uri": "https://localhost:8080/"
},
"outputId": "3e1a95b7-3eae-4100-955b-e84af71417e4"
},
"execution_count": 2,
"outputs": [
{
"output_type": "stream",
"name": "stdout",
"text": [
"601.6182770109472\n"
]
}
]
},
{
"cell_type": "code",
"source": [
"x_SLPM = np.linspace(300,500,100)\n",
"#x_SLPM = 1\n",
"dT_gas = 20\n",
"y = lambda HX_eff: x_SLPM/60./xe_GasDen*(xe_Hlat*dT_gas + xe_Hvap)*(1 - HX_eff) #Watts (SLPM*(Min/60s)/(L/mol)*(J/mol))\n",
"for i in [0.9, 0.95, 0.97]:\n",
" plt.plot(x_SLPM,y(i),label=i)\n",
"plt.legend()\n",
"plt.xlabel('Flow (SLPM)')\n",
"plt.ylabel('Condenser Power (Watts)')\n",
"plt.title('LZ Condener w/ different HX effiencies');"
],
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/",
"height": 295
},
"id": "KKJUpa3KPDGu",
"outputId": "885607ef-9bfc-4413-ca5d-1cdd95bd2a81"
},
"execution_count": 3,
"outputs": [
{
"output_type": "display_data",
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEWCAYAAACJ0YulAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3dd3xc5ZX4/89Rl1zULdmW5SZXGWNAGIgppllACCRZQkkooSwkgW9IIcHON7+EZJMl2S9Lwi7ZZEmjpDgsgYQlYJtmerMpwXLDGGPLQbKa5aIund8fzzOjkawysjUaaXTer9e8PHPbnBmN77n3uc89j6gqxhhjDEBctAMwxhgzfFhSMMYYE2RJwRhjTJAlBWOMMUGWFIwxxgRZUjDGGBNkScEYY0yQJQVzxERkrYhcF+04IkVEpomIikiCf/2EiFwVMv8HIlItIhX+9adEZJeIHBCRY6IV93AhIqki8r8iUi8i/+OnBb8zESn031V8BGM4RUS2RGr7scSSwjAiIjtE5Kwepj/h/9OEPpr8jqqwl22NF5GfishOv/z7/nVO5D/JyCAiK0TkXwe6nqqeq6r3+W0UAl8H5qtqvl/kDuAmVR2rqm8NXsT96y9Bd09wIdPvFZEf+OfHiMg+ESkKmX+ciOwVkWmHEdZFQB6Qraqf6f6dqepO/121H8a2w6KqL6jqnEhtP5ZYUhgB/E5obOABpAOvAver6s7uy4tIEvA0UAycA4wHTgJqgMVDF3nkdN+pHaaPA48f4TYKgRpV3RMybSpQdjgbG6TPdUR8Irsb+KU4icBvgO+o6o7D2ORUYKuqtvnXPX1nZrhQVXsMkwewAzgrjOV+DPwdSOtl/nVAJTC2j23MA9YCe3E7sAtC5t0L/Az4G7AfeA2YGTL/bGAzUI/beTwHXBcy/xpgE1AHrAamhsxT4AvAe/69fwbIANa90a/7QQ+f6T7g6/755MDy/vVMoBaI868zgT1AfA/biccd7VcD2/17KpDg56/13/FZQCPQARwA/uj/VeAg8L5ffhLwZ6AK+AD4csh73QY8BPwO2Oe3mw78GvgI2A38IBAn8HngRR9fnd/euX7eD4F2oMnHcXcPn21a6Gfp9jf/QcjrZP83vgH4LvBS4Lvr5fc0F3jSf8dbgIv99O8BLUCrj+mGbt/Zvd1jOtzP7+dnAb8F/uHn/8VPXwqUhyzX199kMbDO/z0qgTujvW8YykfUA7BHyB8jjKQAXIjbmc7qY5mVwH19zE8EtgHfApKAM3A7/zl+/r10nlUkAL8HVvp5OX7Zi/x2vgq04ZOCj28bLukkAN8GXg55bwUeAzJwR4xVwDkDWPdJ/x8/tYfPdQ3wv/75Z4H3gT+FzPtryLKXAn/s5fv5Am6HOMW/17P0kBT88y47m5A4i/zzOGA98B3/Xc/AJZpSP/823A7zk37ZVOAR4L+BMcAE4HXgBr/85/3y/4xLXl/E7QCle2y9fLZphJEU/LQluN/aPmBuH9scA+wCrvZ/t2NwCXV+yGf8XcjyXb6z7jEd4ef/G/AnXNJPBE7r/p5h/E1eAa7wz8cCJ0Z73zCUj6gHYI+QP0Y/SQF3tFsH/FM/23kS+FEf808BKgg58sMd5d7mn98L/Cpk3nnAZv/8SuDVkHkClNO5k3wCuDZkfhzQgD/i9//5Tw6Z/yCwfADrnhHG9xMH/AJ3VBrYEdwHfC1k2QcC//F72M4zwBdCXi/j8JPCCcDObvNXAL/1z28Dng+Zlwc0E5L0gMuAZ/3zzwPbQual+ffL7x5bL59tml9+b7dHC4cmhXTczv2lfn5vlwAvdJv238B3Qz5jWEnhSD4/MBF3BpLZQ4zB9wzjb/I87gwnZ6D/h2PhYdcURggRScE1M/xGVf/cz+I1uP8gvZkE7FLVjpBpH+KaXAIqQp434I6YgusGZqj7X7QrZNmpwF3+ouReXHOChLntcNYNfa8uVPV9XLPNIlziewz4h4jMAU7DNXMhInG4JrBVvWyqy2fEfTeHayowKfCZ/Of6Fm7nF9D9+0sEPgpZ/r9xR8wBwe9PVRv807EMTI6qZgQewB96WObfcd9ZgYhc2se2pgIndPuMn8PtqAfqSD7/FKBWVevCeI++/ibXArOBzSLyhoicfxifY8SK+kUtE7af4dpgbw1j2aeAH4jIGFU92MP8fwBTRCQuJDEUAlvD2PZHuP98AIiIhL7G7eB+qKq/D2Nb3YWzrvazjedwTVtJqrpbRJ4DrsI1J7ztlzke+FBVq3rZRpfPiPtuDtcu3PWPWX0sE/qZduGOlHO088LsQPT3/YTF94K7AJiPa0b8rYisUdXaHhbfBTynqmcPwlsfyeffBWSJSIaq7u1nuV7/Jqr6HnCZP3j4NPCQiGT38n8p5tiZwvCTKCIpIY8EEbkGOB+4JMz/KA/gfvh/FpG5IhInItki8i0ROQ934bgB+KaIJIrIUuATuGsR/fkbUCwin/Y9Zb5M1yPCXwArRKQYQETSReQzYX72I1k34DngJlwTALjmlJuAF7Wzy+N5/nP05kHgyyJSICKZwPIBxhDqdWC/iNzq++vHi8gCETm+p4VV9SNgDfDvvltxnIjMFJHTwny/Slwb+WETkTHAPcBXVbVaVR/HNUn+pJdVHgNmi8gV/veUKCLHi8i8gb73kXx+v+4TwH+JSKaP49QeFu3zbyIil4tIrj9gCiSXjh62E5MsKQw/j+N6ZwQet+EuuGYBW3u4X+GU7htQ1WZcz5jNuP/M+3D/EXKA11S1BZcEzsW1Gf8XcKWqbu4vOFWtBj4D/AjXTDUL1zMlMP8RXO+olSKyD9jg36dfR7JuiOeAcXQmhRdx7c7PhyzTX1fUX+J6Pr0DvAk8PMAYgnwiOh/XpPUB7vv+Fa69vjdX4i6AbsRdI3mIvpsDQ90FXCQidSLyH4cZ9r/iriGFnrF9BThXRA45G1DV/bjrLpfizkIrcH/H5MN8/yP5/FfgLkRvxvUu+0oP8fb3NzkHKBORA7jv81JVbTzMzzLiBK7YGzMqiEge8BYwWe3Hb8wh7EzBjDbpuHsZLCEY0wM7UzDGGBNkZwrGGGOCRnSX1JycHJ02bVq0wzDGmBFl/fr11aqa29O8EZ0Upk2bxrp166IdhjHGjCgi0usNmdZ8ZIwxJsiSgjHGmCBLCsYYY4JG9DWFnrS2tlJeXk5TU1O0Q4mYlJQUCgoKSExMjHYoxpgYE3NJoby8nHHjxjFt2jRcrbbYoqrU1NRQXl7O9OnTox2OMSbGRLz5yBebektEHvOv7xWRD0Tkbf9Y5KeLiPyHiGwTkb+LyLGH835NTU1kZ2fHZEIAEBGys7Nj+kzIGBM9Q3GmcDNueMXxIdO+oaoPdVvuXFxxtVm4QTB+7v8dsFhNCAGx/vmMMdET0TMFESnAVaT8VRiLX4gbiF5V9VUgQ0TCrYxojDGjwo7qg/z3c+/z8vvVEdl+pJuPfgp8k0Nrkf/QNxH9REQC5XUn03UEqnK6jrgFgIhcLyLrRGRdVVVvY6RE36pVq5gzZw5FRUX86Ec/OmT+hx9+yJlnnsnChQtZunQp5eXlUYjSGDPcqSobdtdz55otlP7keZbesZbbn9jMC+9FJilErPnID2G3R1XX+0FcAlbg6q0n4QbyuBX4frjbVdV7/HqUlJQMy2p+7e3t3HjjjTz55JMUFBRw/PHHc8EFFzB//vzgMrfccgtXXnklV111Fc888wwrVqzggQceiGLUxpjhor1DWbejltVllazZWEF5XSNxAiXTsvjO+fNZVpxHQWZaRN47ktcUlgAX+JG+UoDxIvI7Vb3cz28Wkd8Ct/jXu+k6BGKBnzbivP766xQVFTFjhhsA69JLL+Wvf/1rl6SwceNG7rzzTgBOP/10PvnJT0YlVmPM8NDc1s7L22pYtaGCpzZVUnOwhaT4OE6elcP/OaOIs+blkT32cMctCl/EkoKqrsCdFeDPFG5R1ctFZKKqfuTH9v0kbnQtgEeBm0RkJe4Cc70fXu+wfe9/y9j4j31HsolDzJ80nu9+orjPZXbv3s2UKZ35raCggNdee63LMkcffTQPP/wwN998M4888gj79++npqaG7OzsQY3XGDN8HWhu49nNe1hdVsHaLVUcaG5jbHICp8+dQGlxHkvnTGBs8tDeORCN+xR+LyK5gOAGUv+Cn/44buzcbbjxg6+OQmxD5o477uCmm27i3nvv5dRTT2Xy5MnEx8dHOyxjTIRVH2jmqY2VrC6r4KVtNbS0d5A9JonzF06kdEE+H5uZTXJC9PYFQ5IUVHUtbgB1VPWMXpZR4MbBfN/+jugjZfLkyeza1XnNvLy8nMmTu14znzRpEg8/7Ib+PXDgAH/+85/JyMgY0jiNMUOjvK6B1WWVrN5QwboPa+lQKMhM5YqTplJanM9xUzOJjxseXc1j7o7m4eD444/nvffe44MPPmDy5MmsXLmSP/zhD12Wqa6uJisri7i4OG6//XauueaaKEVrjBlsqsrWygOsLqtgdVkFZb4Ze27+OG46YxalxXnMnzh+WN5zZEkhAhISErj77rspLS2lvb2da665huLiYr7zne9QUlLCBRdcwNq1a1mxYgUiwqmnnsrPfvazaIdtjDkCHR3K2+V7Wb3BJYIdNQ0AHFuYwbfOm8uy+flMyxkT5Sj7N6LHaC4pKdHug+xs2rSJefPmRSmioTNaPqcxw1lrewevbq9hdVkFT26spHJfMwlxwkkzsyktzmfZ/DwmjE+JdpiHEJH1qlrS0zw7UzDGmAFoaGnj+a1VrC6r5OlNlexraiM1MZ6lc3IpLc7n9LkTSE8duRWMLSkYY0w/9ja08NQm13X0hfeqaGrtICMtkbPn53POgnxOmZVDSmJs9B60pGCMMT2oqG9izcYKVm2o4LUPamnvUCamp3BJyRRKi/NZPD2LhPjYG6fMkoIxxnjvVwV6DFXyzq69AMzIHcMNp86gtDifhQXpw7LH0GCypGCMGbVUlXd31wcTwbY9BwBYWJDON0rnUFqcR9GEcVGOcmhZUjDGjCpt7R28saOO1WUVrCmr4B/1TcTHCYunZXH5CYUsK85nUkZqtMOMGksKEbJq1Spuvvlm2tvbue6661i+fHmX+R9++CHXXHMNVVVVZGVl8bvf/Y6CggIA4uPjOeqoowAoLCzk0UcfHfL4jYklTa3tvPheNavLXLG5uoZWkhLiOHVWDl89ezZnzcsjc0xStMMcFiwpRMCRls5OTU3l7bffjlb4xsSEfU2tXYrNNbS0My45gTPmTaC0OJ/TZucyZoiLzY0E9o1EgJXONiY6qvY386QvNvfy+9W0tiu545L55DGTKS3O56QZ2SQlxF6PocEU20nhieVQ8e7gbjP/KDj30JHUQh1p6eympiZKSkpISEhg+fLlljCM6cPOmoZgjaH1O+tQhanZaVy9ZDqlxXkcMyWTuGFSbG4kiO2kMIz1VTr7ww8/ZPLkyWzfvp0zzjiDo446ipkzZ0Y5YmOGB1Vlc8V+VvkaQ5sr9gMwb+J4bj5zFqXF+czNHxfzXUcjJbaTQj9H9JFypKWzA8vOmDGDpUuX8tZbb1lSMKNaR4fy5s66YNfRnbUNiEDJ1Ey+/fF5LJufT2F2ZIanHG1iOylEyZGUzq6rqyMtLY3k5GSqq6t56aWX+OY3vxmNj2FMVLW0dfDy+9WsLqvkyY2VVB9oJjFeWFKUwxeXzuSseXnkjov88JSjjSWFCDiS0tmbNm3ihhtuIC4ujo6ODpYvX97lArUxsexgcxvPba1i1YYKnt28h/3NbYxJimfp3Aksm5/HGXMnMC5l5BabGwmsdPYINVo+p4l9tQdbeGpTJWvKKnj+vWpa2jrIGpPEWb7r6JKi2Ck2N1xEtXS2iMQD64Ddqnq+iEwHVgLZwHrgClVtEZFk4H7gOKAGuERVd0Q6PmPM0Nu9t5E1vsfQ6x+44SknZ6TyuRMKWTY/n+OnZcZksbmRYCiaj24GNgHj/esfAz9R1ZUi8gvgWuDn/t86VS0SkUv9cpcMQXzGmAhTVV9srpJVGyp4d3c9ALMmjOVLS4soLc5nweThOTzlaBPRpCAiBcDHgR8CXxP3Fz8D+Kxf5D7gNlxSuNA/B3gIuFtEREdy+5Yxo1hHh/JO+V5Wl7mmoe3VBwFYNCWDW8+ZS2lxHjNyx0Y5StNdpM8Ufgp8EwiUGcwG9qpqm39dDgT6ak4GdgGoapuI1Pvlq0M3KCLXA9eDqwtkjBk+Wts7eP2DWl9srpKKfa7Y3Ekzsrl6yTTOnp9PfvrwG57SdIpYUhCR84E9qrpeRJYO1nZV9R7gHnAXmgdru8aYw9PU2h4cnvKpTZXUN7aSkhjHqbNy+UbxHM6cN4GMNCs2N1JE8kxhCXCBiJwHpOCuKdwFZIhIgj9bKAB2++V3A1OAchFJANJxF5yNMcNMfUMrz2ypZPWGSp7bWkVjazvjUxI4a14ey3yxudQk6zE0EkUsKajqCmAFgD9TuEVVPyci/wNchOuBdBXwV7/Ko/71K37+MyP5esLhls5+9tln+epXvxpcbvPmzaxcudLqH5moq9zXxJqN7vrAK+/X0Nah5I1P5qLjCigtzueEGVkkWo+hES8aN6/dCqwUkR8AbwG/9tN/DTwgItuAWuDSKMQ2KI6kdPbpp58eLJtdW1tLUVERy5Yti9ZHMaPcB9UHg8Xm3trph6fMGcN1p8ygtDiPowsyrNhcjBmSpKCqa4G1/vl2YHEPyzQBnxmKeCJtsEpnP/TQQ5x77rmkpVlNFzM0VJWyf+xjTVkFq8oq2FrphqdcMHk8Xz97NqUL8pk1Yax1HY1hMV3m4sev/5jNtZsHdZtzs+Zy6+Jb+1zmSEtnB6xcuZKvfe1rgxq/Md21dyjrdtSyusyNQ7B7byNxAsdPy+I7589nWXEeBZl2YDJaxHRSGM76Kp0N8NFHH/Huu+9SWloaxShNrGpqbXfF5ja4HkM1B1tISojjlKIcbj5zFmfOm0D2WCs2NxrFdFLo74g+Uo60dDbAgw8+yKc+9SkSE634lxkc+5taWbulilVlFazdvIeDLe2MTU7g9LkTOKc4n9Pm5DLWhqcc9ewXEAFHUjo74I9//CO33377UIZtYlD1gWae8sNTvrSthpb2DnLGJnHBokksK87nYzOzSU6wrqOmkyWFCDiS0tkAO3bsYNeuXZx22mlR/BRmpNpV28AanwjW7XDF5qZkpXLlSVMpXZDPsYWZxFuPIdOLsEtni0gmMAloBHaoakckAwuHlc6O/c9p+qeqbK08EOw6WvaPfQDMzR/HsuJ8zinOZ95EG57SdDrs0tkikg7cCFwGJAFVuLuT80TkVeC/VPXZQY7XGNOPjg7lrV17fY2hCnbUuOEpjy3M5FvnzWXZ/Hym5YyJdphmBOqv+egh3BgHp6jq3tAZInIccIWIzFDVX/e4tjFm0LS0dfDq9hpWl1Xw5MZK9ux3w1OeNDOHfz51BmfPz2PCOCs2Z45Mn0lBVc/uY9563CA5w46qxvSp8giu/mEGqKGljee2VLG6rIKnN+9hf1MbaUnxLJ2TS2lxPkvnTCA91XqomcET1oVmEVkCvK2qB0XkcuBY4C5V/TCi0R2GlJSU4E1gsZgYVJWamhpSUuyIMFbVHWzh6c17WF1WwfNbq2hu6yAjLZHS4nxKi/M5ZZYNT2kiJ9zeRz8HjhaRo4GvA7/CNSsNu+4xBQUFlJeXU1VVFe1QIiYlJYWCgoJoh2EG0Uf1jazxdxS/9kEt7R3KxPQULj1+CqUL8lk8LcuGpzRDItyk0KaqKiIXAner6q9F5NpIBna4EhMTmT59erTDMKZfbnjKClZvqOCdcjc85czcMdxw6gxKi/NZWJAek2e7ZngLNynsF5EVwOXAqSISB1hDpjEDoKq8u7vedx2tZNseV2xuYUE63yidQ2lxPkUTbHhKE13hJoVLcOMqX6uqFSJSCPy/yIVlTGxoa+/g9R21rPHjFP+j3g1PuXhaFpefUMiy4nwmZaRGO0xjgsJNCl9V1WAhIVXdKSLFEYrJmBGtqbWdF9+rZnVZBU9tqqSuoZXkhDhOmZXLV8+ezVnz8sgcY8NTmuEp3KRwNm5wnFDn9jDNmFFpX1Mrz/oeQ2u3VNHQ0s64lATOCCk2l5ZkVWXM8NffHc1fBL4EzBSRv4fMGge8HMnAjBnu9uxv4smNlawuq+SV96tpbVdyxyXzyWMmU1qcz0kzsklKsB5DZmTp79DlD8ATwO1A6CDD+1W1tq8VRSQFeB5I9u/zkKp+V0TuxXVlrfeLfl5V3xbXzeIu4DygwU9/c4Cfx5iI2lnTEKwxtH5nHaowNTuNq5dMp7Q4j2OmZNrwlGZE6++O5nqgXkTaut+oJiIPqOoVfazeDJyhqgdEJBF4UUSe8PO+oaoPdVv+XGCWf5yAuzfihAF8FmMGnaqy6aP9wUSwuWI/APMnjucrZ86mdEEec/Ks2JyJHeE2cna5qCwiCcBxfa2grhbDAf8y0T/6qs9wIXC/X+9VEckQkYmq+lGYMRozKNo7lDd31rF6QwVrNlays9YVmyuZmsm3Pz6P0uJ8pmTZ8JQmNvV3TWEF8C0gVUT2BSYDLcA9/W1cROJx9ZGKgJ+p6mv+OsUPReQ7wNPAclVtBiYDu0JWL/fTLCmYiGtua+fl92tY44vNVR9oISk+jo8VZfPFpTM5a14eueNseEoT+/prProduF1EblfVFQPduKq2A4tEJAN4REQWACuAClwp7ntwPZi+H+42ReR64HqAwsLCgYZkTNCB5jbWbtnD6rJK1m7ew/7mNsYkxbN07gRKi/M5fU4u41LsHk0zuoTVfKSqK/wgO7Nw4ykEpj8f5vp7ReRZ4BxVvcNPbhaR3wK3+Ne7gSkhqxX4ad23dQ/+LKWkpMTKhZoBqT3YEhye8oVt1bS0dZA1JonzjppI6YI8PjbTis2Z0S3cKqnXATfjdtRvAycCrwBn9LFOLtDqE0Iq7l6HHweuE/jeRp8ENvhVHgVuEpGVuAvM9XY9wQyG3XsbWVNWwaoNFbzhh6ecnJHK504opLQ4n+OnZdnwlMZ44V5ovhk4HnhVVU8XkbnAv/azzkTgPn9dIQ54UFUfE5FnfMIQXIL5gl/+cVx31G24LqlXD+yjGOOoKtv2HAjWGHp3t+v9PDtvLDeeXkRpcT7Fk8ZbjyFjehBuUmhS1SYRQUSSVXWziMzpawVV/TtwTA/Tezy78L2ObgwzHmO66OhQ3infy2pfY2h79UEAjinMYPm5c1k2P48ZuVZszpj+hJsUyv3F4r8AT4pIHTDsBtgxo0trewevf1DLqg0VrNlYQeW+ZhLihBNnZHP1kmmcPT+f/HQbjMiYgeivS+oi4B1V/ZSfdJu/YJwOrIp0cMZ019jSzvPv+eEpN+2hvrGVlMQ4Tpvthqc8c24e6WnWY8iYw9XfmcKvgBkish5X6+gl4BVV3R/xyIzx6htaeXqz6zH03NYqmlo7SE9N5Mx5ruvoqbNySU2yHkPGDIb+7lMoEZE0YDHwMeDLwAMiUgG8pKpfGoIYzShUua+JNRsrWb2hgle319DWoeSPT+HikimUFuezeHoWiTY8pTGDrt9rCqraAKwVkTeA14AlwJXAORGOzYwyH1QfDNYYemvnXgBm5IzhulNmcM6CfBZOTrdic8ZEWH/XFD6LO0NYhCtwF0gMJ6tqReTDM7FMVSn7x75gItha6UplHTU5nVuWzQ4OT2ldR40ZOv2dKfw3sAX4BfC8qm6NfEgmlrV3KG/sqGV1WQVryirZvbeROIHF07P47ifmc/b8PAoyrdicMdHSX1LIAI7GnS3c5u9N+Ah3N/MrqvpMhOMzMaCptZ2XtlUHewzVHGwhKSGOU4pyuPmsWZw1L48sG57SmGGhvwvN7cCb/nG3iOQBnwG+gitiZ10+TI/2N7Xy7BbXdXTt5j0cbGlnbHICp4cMTzk22YanNGa46e+awkLcWULgkYTrmvqfuO6pxgRVH2j2w1NW8PK2GlraO8gZm8QFiya54SlnZpOcYMcRxgxn/R2q3Yvb+T8BfFtVd0Y8IjOi7KptCF4feOPDWlRhSlYqV540lXMW5HNMYaYVmzNmBOmv+ejYoQrEjAyqypbK/aze4M4INn7kxl6amz+OL58xi9LifOZNtOEpjRmp+ms+ehQ3dsFqVW3tNm8G8Hlgh6r+JmIRmqjr6FDe2lXH6jKXCD6sccNTHleYyf89bx7LivOYmj0m2mEaYwZBf81H1wNfA+4SkVqgCjfIzjTgfeBuVf1rRCM0UdHS1sGr22tc09DGSqr2N5MYL5w0M4cbTp3JWfMnMGGcFZszJtb013xUAXwT+KaITMONkdAIbPV3OpsY0tDSxnO+x9DTm/ewv6mNtKR4ls5xxeZOnzuB8TY8pTExLew+gaq6A9gRsUhMVNQdbOHpzXtYtaGCF96rormtg8y0RM4pzqe0OJ+TZ9nwlMaMJtZRfBT6qL6RNf76wGsf1NLeoUxMT+GyxYHhKTNJsGJzxoxKlhRGicDwlGvKKnin3A1PWTRhLF84bQalxfkcNTndegwZY/pPCn6M5ftV9XMD2bCIpADPA8n+fR5S1e+KyHRgJZANrAeuUNUWEUkG7geOA2qAS3yTlTkMqsrfy+uDxeber3LDUx5dkM43SucEi80ZY0yocEpnt4vIVBFJUtWWAWy7GThDVQ+ISCLwoog8gevN9BNVXSkivwCuBX7u/61T1SIRuRT4MXDJgD/RKNbW3sHrO2qDTUMf1TcRHyecMD2LK0+axrLiPCamp0Y7TGPMMBZu89F24CV/38LBwERVvbO3FVRVgQP+ZaJ/KHAG8Fk//T7gNlxSuNA/B3gIV2tJ/HZML5pa23nhvUCxuUrqGlpJTojj1Nm5fH3ZHM6cO4FMKzZnjAlTuEnhff+IA8aFu3Hf9LQeKAJ+5rexV1Xb/CLlwGT/fDKwC0BV20SkHtfEVN1tm9fj7p+gsLAw3FBiSn1jK89u3hMcnrKhpZ3xKQmcOS+P0uI8Tp2dS1qSXS4yxgxcWHsOVf0egIikDeT+BF9ldZGIZACPAHMPK8qu27wHd5c1JSUlo+YsYk9geMoyNzxla7syYVwynz52MqXF+Zw4I9uGp5E2B0AAACAASURBVDTGHLGwkoKInAT8GhgLFIrI0cAN4Y7RrKp7ReRZ4CQgQ0QS/NlCAbDbL7YbmAKUi0gCkI674DxqfVgTGJ6ykjd31qEKU7PTuGbJdJYV53PMlAwbntIYM6jCbWP4KVAKPAqgqu+IyKl9rSAiuUCrTwipwNm4i8fPAhfheiBdBQTKZDzqX7/i5z8z2q4nqCobP9rH6rJK1pRVsLliPwDzJ47nK2fO5pwF+czOs+EpjTGRM5A7mnd12xm197PKROA+f10hDnhQVR8TkY3AShH5AfAW7gwE/+8DIrINqAUuDTe2kay9Q3lzZx2rN1SwemMFu2obEYHjp2bx/50/n2Xz85iSZcNTGmOGRrhJYZeIfAxQ3730ZmBTXyuo6t+BY3qYvh1Y3MP0JtyobjGvua2dl9+vYU1ZBU9u3EP1gWaS4uNYUpTNjUuLOGt+Hjljk6MdpjFmFAo3KXwBuAvXQ2g3sAa4MVJBxaIDzW2s3bKH1WWVrN28h/3NbYxJimepH55y6ZxcxlmxOWNMlIWbFA4M9I5mA7UHW3hqYyWryip4cVs1LW0dZI9J4ryjJlK6II+PzbRic8aY4SXcpLBBRCqBF/zjRVWtj1xYI1d5XUPwjuI3dtTSoTA5I5XLT5hKaXEeJdOybHhKY8ywFe59CkUiUgicAnwc+JmI7FXVRRGNbgRQVd7bcyB4oXjDbjc85ey8sdx4ehGlxfkUTxpvPYaMMSNCuPcpFABLcEnhaKAMeDGCcQ1rHR3KO+V7g8NTflDtKn8cU5jB8nPnUlqcz/QcG57SGDPyhNt8tBN4A/hXVf1CBOMZtlrbO3hte60fnrKCyn3NJMQJJ83M5pqTp7Nsfh554214SmPMyBZuUjgGOBn4rIgsB94DnlPVX/e92sjW2NLOc1urWOOHp6xvbCUlMY7TZrvhKc+cm0d6mvUYMsbEjnCvKbwjIoGieKcAlwOn0XnjWcyob2jl6c2VrNpQwfPvVdHU2kF6aiJnzptAaXE+p87KJTXJegwZY2JTuNcU1uEGy3kZ1/voVFX9MJKBDaXKfU2s8TWGXt1eQ1uHkj8+hYtLplBanM/i6VlWbM4YMyqE23x0rqpWRTSSIba96kDwQvHbu/YCMCNnDNedMoNzFuSzcHK6FZszxow64SaFFhG5EwgUwXsO+P5IvVfhofXl3PI/7wBw1OR0blk2Ozg8pXUdNcaMZuEmhd8AG4CL/esrgN8Cn45EUJF2yqwcvvuJ+Swrzmdyhg1PaYwxAeEmhZmq+k8hr78nIm9HIqChkDc+hauXTI92GMYYM+yEe/W0UURODrwQkSVAY2RCMsYYEy0DqZJ6v4ik+9d1uAFxjDHGxJB+k4KILAKKcIPe7AZQ1X0RjssYY0wU9Nl8JCLfAR4E/gn4G3CJJQRjjIld/Z0pXAIsUtUGEckGVgG/jHxYxhhjoqG/C83NqtoAoKo1YSwfJCJTRORZEdkoImUicrOffpuI7BaRt/3jvJB1VojINhHZIiKlh/OBjDHGHL7+zhRmiMij/rkAM0Neo6oX9LFuG/B1VX1TRMYB60XkST/vJ6p6R+jCIjIfd92iGJgEPCUis1W1fQCfxxhjzBHoLylc2O31HT0u1QNV/Qj4yD/fLyKbcGM89/VeK1W1GfhARLYBi4FXwn1PY4wxR6bPpKCqzw3Gm4jINFz57ddwg/XcJCJXAutwZxN1uITxashq5fSQRETkeuB6gMLCwsEIzxhjjBfx0p8iMhb4M/AV33Pp58BMYBHuTOLfB7I9Vb1HVUtUtSQ3N3fQ4zXGmNEsoklBRBJxCeH3qvowgKpWqmq7qnbgejIt9ovvBqaErF7gpxljjBki/SYFEYkXkbCvJYSsJ7hBeDap6p0h0yeGLPYpXKE9gEeBS0UkWUSmA7OA1wf6vsYYYw5fv3c0q2p7aN2jAViCq6b6bkjxvG8Bl/m7pBXYAdzg36dMRB4ENuJ6Lt1oPY+MMWZohVv76C3fFfV/gIOBiYEmoZ6o6ou4bqzdPd7HOj8EfhhmTMYYYwZZuEkhBagBzgiZpkCvScEYY8zIE1ZSUNWrIx2IMcaY6Aur95GIzBaRp0Vkg3+9UES+HdnQjDHGDLVwu6T+ElgBtAKo6t9xJSmMMcbEkHCTQpqqdu8e2jbYwRhjjImucJNCtYjMxF1cRkQuwtc1MsYYEzvC7X10I3APMFdEdgMfAJdHLCpjjDG9U4X2VkhIGvRNh9v7aDtwloiMAeJUdf+gR2KMMaar9jao+wCqt7pHlf+3+j046UuwdPmgv2VYScEPkPNbYD/wSxE5FliuqmsGPSJjjBltWg527uyrtkD1Fve85n3oaO1cbtxEyJkFCy+GguMjEkq4zUfXqOpdfjS0bFz5igcASwrGGBMOVThY5Xf6W7se/e8r71xO4iFrOuTMgdnnQO4c9zynCFLSIx5muEkhUK7iPOB+X6eopxIWxhgzunW0Q90Od6RfvSWk2WcLNNV3Lpc4xh31T1vi/s2Z7Xb+WTMicq0gXOEmhfUisgaYDqzww2t2RC4sY4wZ5loaoGab3+mHHP3XvA/tzZ3LjZngjvYX/JPf8c92r8dNgriID2kzYOEmhWtxg+JsV9UGEckGrPSFMSb2HawO2fGHHP3v3YXvpQ8SB5nT3A6/6KzOHX/OLEjNjGb0AxZu76MOEakE5otIuInEGGNGho4OqN8Z0rtnS+dF38bazuUSUt2OvuB4WHQ55Pomn+yZkJAcvfgHUbi9j34MXIIb6yAwxoECz0coLmOMGXytTZ1NPqEXemveg7amzuXSst3Ofv4F/iLvbJcM0qcMyyafwRTuUf8ngTmq2tzvksYYE20Ntd3a+n2zT92HBJt8EMgodM08M04LafKZDWlZ0Yw+qsJNCtuBRMCSgjFmeOjocF05u9zU5RNBQ3XncgkpkF0Ek46FhZe6I/7cOW5aYmr04h+mwk0KDcDbIvI0IYlBVb8ckaiMMSagrdn16KnutuOv2QatDZ3LpWa6pp4554b07Z/lzgbi4qMX/wgTblJ41D/CJiJTgPuBPNz52j3+Brgs4E/ANNwYzRerap2/7+Eu3L0QDcDnVfXNgbynMWYEa9zb2cwT2tOnbgdoSA/49ELfv/9k379/jksCadlgt08dsXB7H90nIqlAoapuCXPbbcDXVfVNf1/DehF5Evg88LSq/khElgPLgVuBc4FZ/nEC8HP/rzEmVqjCvt2HNvlUb4UDlZ3LxSe55p38o2DBRb69fzZkz4KktOjFPwqE2/voE8AdQBIwXUQWAd9X1Qt6W0dVP8KX11bV/SKyCZgMXAgs9YvdB6zFJYULcXdLK/CqiGSIyES/HWPMSNLW4gq5hd7UFTj6bz3YuVxyutvZF53tu3f6R8ZUiLfe79EQ7rd+G7AYtwNHVd8WkRnhvomITAOOAV4D8kJ29BW45iVwCWNXyGrlflqXpCAi1wPXAxQWFoYbgjEmEpr2+WYe37c/cPRfux20vXO58ZPdzv7YK7qWdBg7wZp8hplwk0KrqtZ3K3cUVpkLERkL/Bn4iqruC92GqqqIaK8r90BV78GN7UBJScmA1jXGHAZV2F/R9YauwNH//pBjtrgEyJrp2vcD/fsDTT7JY6MXvxmQcJNCmYh8FogXkVnAl4GX+1tJRBJxCeH3qvqwn1wZaBYSkYnAHj99NzAlZPUCP80YMxRCa/d3KenwHjTv61wuaZzb2c9Y2vVCb+Y0iE+MUvBmsISbFP4P8H9x3VH/CKwG/qWvFXxvol8Dm1T1zpBZjwJXAT/y//41ZPpNIrISd4G53q4nGBMBzQfcHbyByp2Bi76123uo3T8bFl7SWccnZw6My7cmnxgm7rpuBDYscjLwAvAunU1N38JdV3gQKAQ+xHVJrfVJ5G7gHFyX1KtVdV1f71FSUqLr1vW5iDGjkyoc2HNoHZ/q93qv3R+4qWsIa/eb6BCR9apa0tO8Ps8UROR/6bwn/BD99D56kc5xGLo7s4flFTcWtDEmXO1tsPfDkP79IQXd+qzd78s5RLl2vxl++ms+usP/+2kgH/idf30ZUNnjGsaYwdfSENLkE3L0X7MN2ls6lwvW7r+os4hb7hzX+8eafEwY+kwKqvocgIj8e7dTjf8VEWu3MWaw9VS7v2or1Hev3T89Jmr3m+En3AvNY0RkhqpuBxCR6cCYyIVlTAzraIe9O3sYrnFrz7X7pywO6d8fW7X7zfATblL4KrBWRLbjrhNMBW6IWFTGxIJg7f4tXUs61GzrVrs/J6Rv/+zO/v3jC2K+dr8ZfsKtfbTK358w10/abGMrGOMNuHb/0s66/aO8dr8ZfgZSXOQ4XGXTBOBoEUFV749IVMYMN4fU7g/p5hlauz8+2TXzTDoWjr6sa5OP1e43I0C4BfEeAGYCb9N1OE5LCia2dK/dHzj6t9r9ZpQI90yhBJivkbrTzZih1rj30HF6e6rdn1HomnimndzZ3JM7B8bkRC10YyIp3KSwAXefgpWdMCNHoHZ/9zo+VVvg4J7O5brX7g9077Ta/WYUCjcp5AAbReR1ug7H2esdzcYMmbYWV7enp5IOobX7U9Ldkf6sZZ3VO3PnWO1+Y0IMZDwFY6IrWLu/W9/+ug+go61zuUDt/mMu9wO3+CqeY3Ltrl5j+hFul9TnRCQPON5Pel1V9/S1jjGHZaC1+yfMDenf78s6JI+LXvzGjHDh9j66GPh/uJHXBPhPEfmGqj4UwdhMLDukdv/Wzj7+vdbuD7nQa7X7jYmIcJuP/i9wfODsQERygacASwqmb80Hut7QZbX7jRnWwk0Kcd2ai2oAu//eOAOq3T/D7fznntdZvjlnFqSMj178xpigcJPCKhFZjRt1DeAS4InIhGSGrWDt/u7DNW7tWrs/aWxI7f7QJp/pVrvfmGEu3AvN3xCRTwMn+0n3qOojkQvLRFXLQXcH70Br9+f6BGC1+40Zsfobea0IyFPVl1T1YeBhP/1kEZmpqu8PRZAmAlQ7a/d3b/Kp39m5nMS5i7o5c2DW2Z1VPHNmQWpG1MI3xkRGf2cKPwVW9DC93s/7RG8rishvgPOBPaq6wE+7DfhnoMov9i1VfdzPWwFci6ut9GVVXR3+xzC96l67P7TZp7Guc7nENHdXb+EJkHOlG6PXavcbM+r0lxTyVPXd7hNV9V0RmdbPuvcCd3No0byfqOodoRNEZD5wKVAMTAKeEpHZqtqOCU9ro6/d362KZ6+1+y+02v3GmEP0lxT6ah/osw6wqj4fRuIIuBBY6cdo+EBEtgGLgVfCXH/0OKR2v3++dye91u4PDtdotfuNMX3rLymsE5F/VtVfhk4UkeuA9Yf5njeJyJXAOuDrqloHTAZeDVmm3E87hIhcD1wPUFhYeJghDHOB2v3BI/6Qo/+Gms7lElJc/Z7Jx8Giz3b29LHa/caYw9RfUvgK8IiIfI7OJFACJAGfOoz3+znwL7hD2n8B/h24ZiAbUNV7gHsASkpKRnYp72Dt/h6Ga+xSuz/LHenP/XjXC71Wu98YM8j6TAqqWgl8TEROBxb4yX9T1WcO58389gAQkV8Cj/mXu4EpIYsW+GmxIbR2f+iF3u61+9MLXft+oHZ/oMnHavcbY4ZIuPcpPAs8e6RvJiITVTVQ1exTuHEaAB4F/iAid+IuNM8CXj/S9xtSgdr9PQ3X2GPt/oUhtftnu2lWu98YE2URKyIvIn8ElgI5IlIOfBdYKiKLcM1HO4AbAFS1TEQeBDYCbcCNw7bnUffa/cFmn35q9wfa+612vzFmGJORPMJmSUmJrlu3LjIbb6r3zTyhwzX6Qm6h+SpQuz9wxB8crtFq9xtjhicRWa+qJT3NG92HrKG1+0NLOlRthQMVncvFJboePRPmdfbvD4zclTw2evEbY8wgG51J4f1n4envuzOBlv2d05PHu149M8/oOlxj5nRr8jHGDBlVpb65nurGaqoaq6hurKamsabL82XTlnHxnIsH/b1H554uMdWNzrXoss7SzVa73xgTYc3tzVQ3VrtHg/+3qZqqhqquO/2mGtpCh5j1UuJTyEnNITctlziJTAWC0ZkUCk+Eqx6NdhTGmBjQoR3UN9cHd+jdH6E7+/2hLROeIGSmZJKbmktOag4zM2YGn3d/jEkcg0T4wHV0JgVjjOlHU1tTjzv20EdVYxW1jbW06aFH9akJqcGdeVFGESdOPNEd5afmkp2aHXyemZJJQtzw2RUPn0iMMSbCOrSDvc17g801gaabno7wD7QeOGT9OIkjMzmT3DS3Yy/KKCI3zR3VZ6dmdznCH5M4Jgqf8MhZUjDGjHiNbY2HHsU3VFHTVNNlWm9H9WkJacGd+ezM2Zw06aTgDj6ws89NyyUzOZP4GC8tY0nBGDMstXe0U9dc547oeziSr2rsvDh7MPTGUS9O4shKyQo218zJnNPjEX1Oag5piVZNIMCSgjFmSDW0NvTbdFPdWE1tUy3tPRQ2GJM4Jrgzn5M1hyWpS7rs4ANJYDQc1UeCJQVjzBELHNX31HQT2PEHXvd0VB8v8WSnZJOd6h5zs+Ye0nSTk+Je21F9ZFlSMMb0qqG1oUtzTU83UVU1VlHbVEtHaMVfb0zimGBTzdysuV2abUJ74WSmZEas370ZGEsKxowygaP6nppuAu30geeNbY2HrJ8gCWSlZpGTmsOEtAnMz57fpZ0+sKPPSc0hNcEGexppLCkYEwNUlYa2hs7mmqaQO2b9XbM1jTVUNVRR11zX41H9uMRx5KS5nfn87PmHtNPnpOWQnZJtR/UxzpKCMcNYW0cbtU21/d5EVd1Y3e9RfV5aHsXZxZ3t9L75Jjctl+yUbFISUqLwCc1wY0nBmCGmqhxsPdhj0033HX9dUx3KoeXtxyWNC+7UF+Qs6NzBd7tbNj053Y7qzYBYUjBmkLR1tAW7Wgaaarq30wd64fR4VB+XQHaK27FPGjuJo3KOCva6yUnNCTbt5KTmkByfHIVPaEYDSwrG9EFVOdB6oNe+9KE7/d6O6scnjQ/uzBfmLuzaVu+bbgJH9ZEudmZMfywpmFGptaOV2sbaQ5pueno0tzcfsn5CXEKwt03B2AKOzj36kLIIged2VG9GkkiO0fwb4Hxgj6ou8NOygD8B03BjNF+sqnXiDo/uAs4DGoDPq+qbkYrNxCZVZX/r/i616rs33QR65dQ11/W4jYzkDLJTsslJy2HRhEXkpOQEi58Fe+Gk5jA+abwd1ZuYFMkzhXuBu4H7Q6YtB55W1R+JyHL/+lbgXGCWf5wA/Nz/awyt7a2HFDYL1r3x3S8DO/6ejuoT4xK7HNUfk3tMjzVwslOzSYpPisInNGb4iFhSUNXnRWRat8kXAkv98/uAtbikcCFwv6oq8KqIZIjIRFX9KFLxmehSVfa17OuzTn1gR7+3eW+P28hIzgju0AsnFB5yp2zgXzuqNyZ8Q31NIS9kR18B5Pnnk4FdIcuV+2mHJAURuR64HqCwsDBykZrD0tre2m9ZhMD8lo6WQ9ZPiksK1qcvHFfIcXnHHdJ0k5PqbqJKjE+Mwic0JrZF7UKzqqqIHNpVo//17gHuASgpKRnw+mbgAkf1wTtlQ3buwaP8BtdeX99c3+M2Akf1uam5FI4v7FoDJ6TNflziODuqNyaKhjopVAaahURkIrDHT98NTAlZrsBPMxHU0t7Sb/NNYF5Pg4gnxycHd+zT06dTkl9yyBF9TmoOWalZJMbZUb0xI8FQJ4VHgauAH/l//xoy/SYRWYm7wFxv1xMOj6r2OYh46GNfy75D1g8MIh66sw9tsgk07eSk5jA2cawd1RsTYyLZJfWPuIvKOSJSDnwXlwweFJFrgQ+Bi/3ij+O6o27DdUm9OlJxjVTN7c1dd+q+uSa0+2XgdU9H9SnxKcGLrzPSZ7A4f3Gw6Sa0J05mSqYd1RszikWy99Flvcw6s4dlFbgxUrEMVx3aQX1zfY8li7sf1e9v2X/I+oGj+kBzzcyMmZ3NNmk5wfIIuWm5pCWk2VG9MaZfdkdzBDS1NfXZxTK0Bk5PR/WpCanBtvmijCJOnHjiIbXqA0f1CXH2JzTGDB7bo4SpQzuoa6rr7HnTFDLsYMh4szWNNexv7fmoPislK3gkX5RRFKx7k5PWeXE2NzXXhhs0xkTNqE8KTW1N/TbdVDe4o/qeBhFPTUgN7tBnZ84mZ1JnkbPQNvuM5Aw7qjfGDHujci/1QvkL/Nsb/0Z1YzUHWg8cMj9O4oKDiOek5jAnc06XIQZDu13aUb0xJpaMyqSQnpzOnKw5LEld0qWrZaAsQmZyJvFx8dEO0xhjhtyoTAoLcxdyx2l3RDsMY4wZdmycPmOMMUGWFIwxxgRZUjDGGBNkScEYY0yQJQVjjDFBlhSMMcYEWVIwxhgTZEnBGGNMkLiq1SOTiFThxmU4HDlA9SCGM1iGa1wwfGOzuAbG4hqYWIxrqqrm9jRjRCeFIyEi61S1JNpxdDdc44LhG5vFNTAW18CMtris+cgYY0yQJQVjjDFBozkp3BPtAHoxXOOC4RubxTUwFtfAjKq4Ru01BWOMMYcazWcKxhhjurGkYIwxJihmk4KIpIjI6yLyjoiUicj3/PTpIvKaiGwTkT+JSJKfnuxfb/Pzpw1xXL8XkS0iskFEfiMiiX76UhGpF5G3/eM7QxzXvSLyQcj7L/LTRUT+w39ffxeRY4c4rhdCYvqHiPzFTx+S7yskvngReUtEHvOvo/r76iOuqP6++ogrqr+vPuKK+u9LRHaIyLv+fdb5aVki8qSIvOf/zfTTB+/7UtWYfAACjPXPE4HXgBOBB4FL/fRfAF/0z78E/MI/vxT40xDHdZ6fJ8AfQ+JaCjwWxe/rXuCiHpY/D3jCr3ci8NpQxtVtmT8DVw7l9xXy3l8D/hB4z2j/vvqIK6q/rz7iiurvq7e4hsPvC9gB5HSb9m/Acv98OfDjwf6+YvZMQZ0D/mWifyhwBvCQn34f8En//EL/Gj//TBGRoYpLVR/38xR4HSgY7Pc+nLj6WOVC4H6/3qtAhohMHOq4RGQ87m/6l8F+7/6ISAHwceBX/rUQ5d9XT3EBRPv31VtcfRiS31d/cUXz99WL0N9R99/XoHxfMZsUIHhK+DawB3gSeB/Yq6ptfpFyYLJ/PhnYBeDn1wPZQxGXqr4WMi8RuAJYFbLKSb755AkRKY5ETP3E9UN/SvoTEUn204Lflxf6XQ5VXOD+UzytqvtCpg3J9wX8FPgm0OFfZzMMfl89xBUUzd9XH3FF9ffVR1wQ3d+XAmtEZL2IXO+n5anqR/55BZDnnw/a9xXTSUFV21V1Ee6oaDEwN8ohAYfGJSILQmb/F/C8qr7gX7+Jq1NyNPCfRPCIpZe4VuC+t+OBLODWSL3/AOMKuAzXHBIwJN+XiJwP7FHV9ZHY/uEKI66o/L76iCuqv68wvq+o/L68k1X1WOBc4EYROTV0pj/rG/R7CmI6KQSo6l7gWeAk3GlVgp9VAOz2z3cDUwD8/HSgZojiOse/73eBXFz7ZmCZfYHmE1V9HEgUkZyhiktVP/KnpM3Ab3HJFUK+Ly/0u4x4XAD+e1gM/C1kmaH6vpYAF4jIDmAlronhLqL/+zokLhH5nX/faP6+eoxrGPy++vq+ovn7QlV3+3/3AI/4WCoDzUL+3z1+8UH7vmI2KYhIrohk+OepwNnAJtxO5SK/2FXAX/3zR/1r/PxnfCYeirg2i8h1QClwmap2hCyfH2h7FpHFuL/ZoO9M+ogr8AMU3Kn0Br/Ko8CVvtfDiUB9yGltxOPysy/CXfRrCll+SL4vVV2hqgWqOg134fgZVf0cUf599RLX5dH+ffURV1R/X73F5WdH7fclImNEZFzgObAM992E/o66/74G5ftK6H+REWsicJ+IxOP+cA+q6mMishFYKSI/AN4Cfu2X/zXwgIhsA2pxP5ChjKsNVwb8Ff+be1hVv4/7YX7Rz2/E9WyJxG3ovcX1jIjk4no1vA18wS//OK7HwzagAbg6AjH1Gpefdynwo27LD9X31Ztbie7vqze/ILq/r978Psq/r75E8/eVBzzi/1YJwB9UdZWIvAE8KCLX4v6eF/vlB+37sjIXxhhjgmK2+cgYY8zAWVIwxhgTZEnBGGNMkCUFY4wxQZYUjDHGBFlSMKOCiLRLZ2XLt0VkmriKl4/1v/aA3meidFbaTBNXnfRdcdVJXxSRsX7egR7WvU1Edvv4NojIBSHTVUSKQpb9ip9W4l8/Jb5ipjFHwpKCGS0aVXVRyGNHhN7na8Av/fObgUpVPUpVFwDXAq39rP8TX9LjM8BvRCTwf/Rdut7b8BmgLOT1A7hKrMYcEUsKxhCsU/8XcYXZXhWRhX76uyKS4e8UrRGRK/30+0Xk7B429U90FpubSEipAVXd4ss59EtVNwFtQKCEwl9wlTARkZm4gnrVIas8iqvTY8wRsaRgRovUkKajR3qY/z3gLVVdCHwLuN9PfwlXH6cY2A6c4qefBLwcugERmQ7Uhez4fwPcKiKviMgPRGRWuMGKyAm4qp1VftI+YJe4YoCXAn8KXV5V64BkEYlU5VUzSsRymQtjQjX6ZpnenIw7ykdVnxGRbHG19F8ATsWVFPg5cL2ITMbt/A9228ZEOnfiqOrbIjIDV7fmLOANETnJnwX05qsicjmwH7hEVVU6h11YiUsIpcCZHFrKYA8wiQgXcjSxzc4UjOnb87izg1OAtbid/kW4ZNFdI5ASOkFVD6jqw6r6JeB3uPo0ffmJv+ZxSkh564DHcGMh7OxW3z8gxcdgzGGzpGCM8wLwOXDj8ALVvkzyLly7/ixV3Q68CNyCSxbdbQWmBV6IyBLpHEM3CZiPO+M4LKragCu498Pu83zlznzcEI7GHDZrPjLGuQ3X2+fvuCqTV4XMew2I989fAG7HJYcuVPWgiLwvIkWqug2YCfzc77DjcHX5/+wXTxOR8pDV7wwnSFVd2cus44BXtXPUaYWwXwAAAF1JREFUN2MOi1VJNWYQicingONU9dtD/L53AY+q6tND+b4m9tiZgjGDSFUfiVIPoA2WEMxgsDMFY4wxQXah2RhjTJAlBWOMMUGWFIwxxgRZUjDGGBNkScEYY0zQ/w8kuWG374pifgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"cell_type": "code",
"source": [
"N2_mu_liq = 1.64e-4 #kg/m/s LN2 at 78K\n",
"N2_rho_liq = 806.4 #kg/m3 \n",
"N2_rho_gas = 4.4 #kg/m3 \n",
"g = 9.8 #m/s^2\n",
"C_pl = 1040 # J/kg/K\n",
"Twall = 77.34 #Ln2\n",
"Tsat = 81.0 #at 1.5 bar\n",
"h_lg = 199000 + 0.68*C_pl*(Tsat-Twall)\n",
"k = 0.139 #W/m/K (1W/1000mW)*(100cm/1m)"
],
"metadata": {
"id": "8FpNEpLiXZWX"
},
"execution_count": 10,
"outputs": []
},
{
"cell_type": "code",
"source": [
"h = lambda L: 0.943*(N2_rho_liq*((N2_rho_liq-N2_rho_gas)*g*h_lg*k**3)/(L*N2_mu_liq*(Tsat-Twall)))**(0.25) # W/m2K\n",
"L = np.linspace(4*2.54/100,15*2.54/100,100)\n",
"plt.plot(L*100/2.54,h(L))\n",
"plt.xlabel('Tube Length')\n",
"plt.ylabel('h (Watts/m2/K)');"
],
"metadata": {
"id": "opzEo9Wk8RDc",
"colab": {
"base_uri": "https://localhost:8080/",
"height": 279
},
"outputId": "ada50b68-4645-4e3c-f17e-77eadd2830d5"
},
"execution_count": 11,
"outputs": [
{
"output_type": "display_data",
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEGCAYAAACUzrmNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3dd3hVZdb38e9KgYD0FiAEgjSllwgqDCoqIqJgQbGPOo+94IzvjE51io+OOtgfFcuoI4qOWHAsiL2gQOiEGimSANI7BJKs94+zGSOSnBByslN+n+s6Fzv3bmuLZJ17383cHRERkeLEhR2AiIhUfEoWIiISlZKFiIhEpWQhIiJRKVmIiEhUCWEHEAtNmjTxtLS0sMMQEalUZsyYscHdmx5sX5VMFmlpaWRkZIQdhohIpWJmK4vaF7PXUGaWamafmNkCM8s0s1sK7bvJzBYF5fcWKr/DzLLMbLGZnVaofEhQlmVmt8cqZhERObhY1izygF+5+0wzqwvMMLPJQDIwHOjh7rlm1gzAzDoDo4AuQEvgQzPrGFzrMeBUIBuYbmYT3X1BDGMXEZFCYpYs3H0NsCbY3m5mC4EU4H+Ae9w9N9i3LjhlODA+KF9uZllA32BflrsvAzCz8cGxShYiIuWkXHpDmVka0AuYCnQEfmZmU83sMzM7JjgsBVhV6LTsoKyochERKScxb+A2szrABGC0u28zswSgEXAscAzwqpkdWQb3uRq4GqB169aHezkRESkkpjULM0skkijGufvrQXE28LpHTAMKgCZADpBa6PRWQVlR5T/i7mPdPd3d05s2PWjPLxERKaVY9oYy4BlgobuPKbTrTeCk4JiOQA1gAzARGGVmNc2sLdABmAZMBzqYWVszq0GkEXxirOIWEZGfimXNoj9wKTDIzGYHn6HAs8CRZjYfGA9cHtQyMoFXiTRcvw/c4O757p4H3AhMAhYCrwbHlrnde/O5571FrNq0KxaXFxGptKwqrmeRnp7upRmUt3rLbk4d8xm92zTkhSv7EqkciYhUD2Y2w93TD7ZPc0MV0rJBLX5z+lF8sXQDb8z6SbOIiEi1pWRxgEv6taF36wb85T8L2LAjN+xwREQqBCWLA8TFGX8/tzs7c/P463807k9EBJQsDqpDcl1uOKk9b81ezceLvg87HBGR0ClZFOG6E9vRMbkOv319Ptv27As7HBGRUClZFKFmQjz3ndeDddv3cNd/FoYdjohIqJQsitEjtQHXnNCOVzJW8dmS9WGHIyISGiWLKG45uQPtm9Xh9glz2a7XUSJSTSlZRJGUGM9953Xn+217uOsdvY4SkepJyaIEerVuyDUntGP89FXqHSUi1ZKSRQmNPqUDRzWvy28mzGPzzr1hhyMiUq6ULEqoZkI8Y87vyZZde/n9m/OpinNqiYgURcniEHRuWY/Rp3TknXlrmDhnddjhiIiUGyWLQ3TNwCPp3boBv39zPjlbdocdjohIuVCyOEQJ8XE8eEEvCgqcX74ym/wCvY4SkapPyaIUWjeuzZ1ndWHq8k2M/XxZ2OGIiMSckkUpndenFUO7NWfM5MXMz9kadjgiIjGlZFFKZsb/nt2NxkfU5OaXZ7EzNy/skEREYkbJ4jA0qF2DMRf0YPnGnfxpYkyWBRcRqRBilizMLNXMPjGzBWaWaWa3BOV3mlmOmc0OPkMLnXOHmWWZ2WIzO61Q+ZCgLMvMbo9VzKVxfLsm3HhSe16bkc1bs7UUq4hUTQkxvHYe8Ct3n2lmdYEZZjY52PeAu99f+GAz6wyMAroALYEPzaxjsPsx4FQgG5huZhPdvcIsY3fLyR34+tuN/O6N+fRMbUCbxkeEHZKISJmKWc3C3de4+8xgezuwEEgp5pThwHh3z3X35UAW0Df4ZLn7MnffC4wPjq0wEuLjeHBUT+IMbnxpFrl5+WGHJCJSpsqlzcLM0oBewNSg6EYzm2tmz5pZw6AsBVhV6LTsoKyo8gPvcbWZZZhZxvr15b/2RKuGtblvZA/m5Wzl7ncXlfv9RURiKebJwszqABOA0e6+DXgcaAf0BNYA/yiL+7j7WHdPd/f0pk2blsUlD9lpXZpzZf+2PDdlBe/OWxNKDCIisRDTZGFmiUQSxTh3fx3A3b9393x3LwCeIvKaCSAHSC10equgrKjyCun204+iR2oDfvPaXFZu3Bl2OCIiZSKWvaEMeAZY6O5jCpW3KHTY2cD8YHsiMMrMappZW6ADMA2YDnQws7ZmVoNII/jEWMV9uGokxPHohb0wg+vHzWTPPrVfiEjlF8uaRX/gUmDQAd1k7zWzeWY2FzgJuBXA3TOBV4EFwPvADUENJA+4EZhEpJH81eDYCiu1UW3GnN+TzNXbuFPjL0SkCrCquC5Denq6Z2RkhB0G901axGOffMu953bn/GNSo58gIhIiM5vh7ukH26cR3DH0y1M70b99Y/7w1nzNHyUilZqSRQzFxxkPjepFw9o1uG7cDLbs0nKsIlI5KVnEWJM6NXn8kt58vzWXm16epfUvRKRSUrIoB71aN+Qvw7vwxdIN3DdpcdjhiIgcsljODSWFjOrbmrk5W3nis2/pmlKPYd1bhh2SiEiJqWZRju48swt92jTk//17Lpmr1eAtIpWHkkU5qpEQx+MX96Z+rUSufmEGG3bkhh2SiEiJKFmUs2b1khh7WR827MjluhdnsDevIOyQRESiUrIIQfdWDbj3vO5MX7GZP741n6o4MFJEqhY1cIdkeM8Ulny/ncc++ZYOyXW5akDbsEMSESmSkkWIfnVqJ7LW7eBv7ywgrXFtTj46OeyQREQOSq+hQhQXZzxwQU+6tKzHzS/PYuGabWGHJCJyUEoWIatdI4FnLj+GukmJ/OL5DNZt2xN2SCIiP6FkUQEk10vi6cvT2bxrL1c9n8GuvXlhhyQi8iNKFhVE15T6PHJhLzJXb+VmzSElIhWMkkUFcvLRyfz5rC58uHAdf347U11qRaTCUG+oCubS49L4btMunvpiOSkNanHNCe3CDklERMmiIrrj9KNZvXUPd7+3iOb1kxjeMyXskESkmovZaygzSzWzT8xsgZllmtktB+z/lZm5mTUJfjYze9jMssxsrpn1LnTs5Wa2NPhcHquYK4q4OOMfI3vQr20jbvv3HL7K2hB2SCJSzcWyzSIP+JW7dwaOBW4ws84QSSTAYOC7QsefDnQIPlcDjwfHNgL+BPQD+gJ/MrOGMYy7QkhKjGfsZem0bXIE1/xrhmapFZFQxSxZuPsad58ZbG8HFgL736c8APwaKNyCOxx4wSO+ARqYWQvgNGCyu29y983AZGBIrOKuSOrXSuT5K/tSLymBy5+dzsqNO8MOSUSqqXLpDWVmaUAvYKqZDQdy3H3OAYelAKsK/ZwdlBVVfuA9rjazDDPLWL9+fRlGH64W9WvxwlX9yC8o4NJnprFuuwbtiUj5i3myMLM6wARgNJFXU78F/ljW93H3se6e7u7pTZs2LevLh6p9szr884q+bNiRy+XPTmfr7n1hhyQi1UxMk4WZJRJJFOPc/XWgHdAWmGNmK4BWwEwzaw7kAKmFTm8VlBVVXq30TG3AE5f0IWvddn7x/HR2780POyQRqUZi2RvKgGeAhe4+BsDd57l7M3dPc/c0Iq+Uerv7WmAicFnQK+pYYKu7rwEmAYPNrGHQsD04KKt2BnZsyoMX9GLGys1cq4WTRKQcxbJm0R+4FBhkZrODz9Bijn8XWAZkAU8B1wO4+ybgr8D04POXoKxaOqN7C+4+pxufLVnP6Fc0LYiIlI+YDcpz9y8Bi3JMWqFtB24o4rhngWfLMr7K7IJjWrN9Tx5/e2chtWvM5d5zuxMXV+x/ahGRw6IR3JXUL352JNv35PHQR0tJSozjr8O7EnnzJyJS9pQsKrHRp3RgT14+T362jKSEeH53xtFKGCISE0oWlZiZcfuQo8jdV8DTXy6nZmIctw3upIQhImVOyaKSMzP+OKwzuXn5PPbJtyTExXHrqR3DDktEqhgliyogLs64a0Q38guchz5aSnyccfPJHcIOS0SqECWLKiIuzrjnnO7kF8CYyUuIM7hxkBKGiJQNJYsqJC7OuPe87rg793+whAJHNQwRKROHlCzM7Ahgj7trrokKKj7OuG9kD7BIDaPAndGnqA1DRA5PscnCzOKAUcDFwDFALlDTzDYA7wBPuntWzKOUQxIfZ9x3Xg/izHjww6UUFDi3ntpRvaREpNSi1Sw+AT4E7gDmu3sB/HdBopOAv5vZG+7+YmzDlEMVH2fce2534s14+OMscvMKuP30o5QwRKRUoiWLU9z9J/Nhu/smM3vb3ScEM8tKBRQXZ9x9TjdqJMTx5OfLyM0r4I/DOmtqEBE5ZNEmErzjYIVmVh/4AOBgyUQqjrg44y/Du/CLAW15bsoKfvvGPE0+KCKHLFqyGGBmdxUuMLNk4DMir6ikEjAzfnfG0dx4UnvGT1/F6Fdmsy9f05uLSMlFSxZnAT3MbAyAmXUAvgKecPc/xzo4KTtmxm2ndeL204/i7Tmrue7FGezZp05tIlIyxSYLd98DnA2kmdnLRBq7/5+7P1EewUnZu/aEdvx1RFc+WrSOK/45nR25eWGHJCKVQLHJwsx+CdwETCWyQt0soK2Z/TLYJ5XQpce24YHzezJtxSYueuobNu7IDTskEangor2Gqht8koCHiSSLuoU+UkmN6JXC2Ev7sHjtdkY++TWrt+wOOyQRqcAsskBdETvNLgImufvG8gvp8KWnp3tGRkbYYVQKU5dt5BfPZ1AnKYEXruxLh2R9BxCprsxshrunH2xftJpFKvBvM/vCzO40s35WwlFdZpZqZp+Y2QIzyzSzW4Lyv5rZ3GBN7g/MrGVQbmb2sJllBft7F7rW5Wa2NPhcXrLHlpLod2Rjxl9zLHkFznlPfM2MldV2eXMRKUaxNYv/HmRWFzgFGAL0BRYC7xOpdXxfxDktgBbuPjM4fwYwAsh2923BMTcDnd39WjMbSqR9ZCjQD3jI3fsFo8UzgHTAg+v0cffNRcWrmsWhW7VpF5c9O401W3fz6IW9OaVzctghiUg5O5yaBQDuvt3d33D3a9y9F/A3oCnwQjHnrHH3mfvPJ5JgUvYnisARRBIAwHDgBY/4BmgQJJzTgMnuvilIEJOJJC0pQ6mNavPatcfRMbkuV/8rg5emfhd2SCJSgZQoWRw4pYe7LwCed/fTSnh+GtCLSK8qzOwuM1tFZILCPwaHpQCrCp2WHZQVVX7gPa42swwzy1i/fn1JwpIDNK5Tk5f/51gGdmzKb9+Yxz8+WExJap4iUvVF6zp7kpllA2uC9oW0Qrs/KMkNzKwOMAEYvb9W4e6/c/dUYBxwY2kCP5C7j3X3dHdPb9q0aVlcslo6omYCT1+WzqhjUnnk4yxu+/dc9uZptLdIdRetZnEvcJq7NwHGApPN7NhgX9SG7qBGMgEY5+6vH+SQccC5wXYOkQb1/VoFZUWVS4wkxMdx9znduPWUjkyYmc3P/zmNrbs1BZhIdRYtWdRw90wAd3+NSAP182Y2gh/aGg4q6DX1DLDQ3ccUKi+8dNtwYFGwPRG4LOgVdSyw1d3XAJOAwWbW0MwaEhkcOKnETyilYmbcckoH/jGyB9NXbOK8x6ewatOusMMSkZBEm6J8n5k1d/e1AO6eaWYnA/8B2kU5tz9wKTDPzGYHZb8FrjKzTkABsBK4Ntj3LpGeUFnALuCK4J6bzOyvwPTguL+4u/p3lpNz+7SiRYMkrvnXDM7+vyk8dVkferVuGHZYIlLOog3KOwVY7+5zDihvANzg7ncd/Mxwqets2ctat50rnpvOum25/OP8Hgzr3jLskESkjJW666y7f3hgogjKt1TURCGx0b5ZXd68vj/dUupz40uzeOSjpeopJVKNlLTr7DAzm2Vmm8xsm5ltN7Nt0c+UqqRxnZq8+It+nN0rhX9MXsKtr8zWNOci1US0Nov9HgTOAea5vk5Wa0mJ8Yw5vwftm9XhvkmLWb5xF09d2odm9ZLCDk1EYqhENQsig+LmK1EIRHpK3XBSe564pA9L1m7nrEe/Yl721rDDEpEYKmmy+DXwrpndsX8tC61nIUO6Nue1644jPs4474kpvDlLw19EqqqSJou7iHRnTULrWUghXVrWZ+KN/emR2oDRr8zm7ncXkl+gCqhIVVPSNouW7t41ppFIpdW4Tk3G/aIff3l7AU9+vowFa7bxyIW9aFC7RtihiUgZKWnN4l0zGxzTSKRSS4yP468juvL3c7sxddkmznz0SxasVoc5kaqipMniOuB9M9utrrNSnAuOac0r1xzLvjznnMe/4q3ZascQqQqizTqbCODudd09zt1ruXu94Od65ROiVDa9Wjdk4k396Z7SgFvGz+bOiZmauVakkotWs/jazN40s2sPmJ5cpFjN6iYx7n/6cdWAtjw3ZQUXPvUNa7fuCTssESmlaNN9pAOjgx8fNLPpZvaAmQ02s5qxD08qs8T4OP4wrDOPXtSLhWu2MeyRL/gqa0PYYYlIKURts3D3Fe7+hLuPAI4H3iayHvcXZvZOrAOUym9Y95ZMvLE/DWrX4NJnpvLIR0spUPdakUqlpA3cALj7PmAW8KK79wWujklUUuW0b1aXt27oz5k9WvKPyUu44rnpbNyRG3ZYIlJCJZ1I8FMzq2dmjYCZwFNmNsbd1dVFSuyImgk8eEFP/jaiK18v28gZD3/J9BVamkSkMihpzaJ+sH72OcAL7t6PyKsokUNiZlxybBtev+54khLjGDX2Gx77JEuvpUQquJImiwQzawGcT2SVPJHD0jWlPm/fNIDTuzbnvkmLufyf01i/Xa+lRCqqkiaLPxNZ9zrL3aeb2ZHA0tiFJdVB3aREHrmwF3ef041pyzdx+kNf8MXS9WGHJSIHUdJkscbdu7v79QDuvgwYU9wJZpZqZp+Y2QIzyzSzW4Ly+8xskZnNNbM3giVa959zh5llmdliMzutUPmQoCzLzG4/9MeUisrMuLBvaybeOICGtRO59Jlp3P3uQg3iE6lgSposHilhWWF5wK/cvTNwLHCDmXUGJgNd3b07sAS4AyDYNwroAgwB/s/M4s0sHngMOB3oDFwYHCtVSKfmdZl44wAu7teaJz9fxnlPTGH5hp1hhyUigWJnnTWz44iMrWh6wPoV9YD44s519zXAmmB7u5ktBFLc/YNCh30DnBdsDwfGu3susNzMsoC+wb6soDaDmY0Pjl1QgueTSqRWjXjuOrsbP+vQlN9MmMsZD3/BnWd2YWR6K8ws7PBEqrVoNYsaQB0iSaXwOhbb+OGXfFTBVCG9gKkH7LoSeC/YTiGyIt9+2UFZUeUH3uNqM8sws4z16/XeuzIb0rU574/+GT1aNeDXE+Zy/biZbNm1N+ywRKq1YmsW7v4Z8JmZ7Xb3ewvvM7ORlKCR28zqABOA0UH32/3lvyPyqmpcaQI/SKxjgbEA6enp6odZybWoX4txv+jHU18s4/4PFjPzwc3cP7IHP+vQNOzQRKqlkrZZjDpI2R3RTgpmrZ0AjHP31wuV/xwYBlxcaF3vHCC10OmtgrKiyqWKi4szrjmhHW9c35+6SZHG7zsnZrJnX37YoYlUO9HaLE4HhgIpZvZwoV31iNQKijvXgGeAhe4+plD5ECJrep/g7rsKnTIReMnMxgAtgQ7ANMCADmbWlkiSGAVcVLLHk6qga0p9/nPTAO55bxHPTVnBl1kbGHN+D7q3ahD9ZBEpE9FqFquBDGAPMKPQZyJwWjHnAfQHLgUGmdns4DMUeJRIu8fkoOwJAHfPBF4l0nD9PnCDu+e7ex5wI5FxHguBV4NjpRpJSoznzrO68MKVfdmxJ4+z/28KD0xewr58dbEVKQ/2w1ugYg4ySwwmEawU0tPTPSMjI+wwJEa27trHnW9n8sasHLql1Ocf5/egY3LdsMMSqfTMbEawNMVPlLTNIs3MXgsG2C3b/ynDGEVKrH7tRB64oCePX9ybnC27Gfbwlzz+6bfka34pkZgpabL4J/A4kXaKk4AXgBdjFZRISZzerQUf3DqQQUc14+/vL+Lcx6eQtW572GGJVEklTRa13P0jIq+tVrr7ncAZsQtLpGSa1KnJ45f05qFRPVmxcSdDg1pGntoyRMpUSZNFrpnFAUvN7EYzO5vIYD2R0JkZw3umRGoZnX6oZSxauy36ySJSIiVNFrcAtYGbgT5EejldHqugREqjWd0kHr+kN49e1ItVm3dz5iNf8sDkJZqUUKQMRBtnMQf4KvhscPflwBXlEZhIaZgZw7q35Ph2Tfjz25k89NFS3p+/lnvO7Uav1g3DDk+k0opWs7gYmA2cCkwys5ygV9StZtYv9uGJlE6jI2rw0KhePHN5Otv27OOcx6dw58RMduYWO5ZURIpQonEW/z3YrAmREdSjgbbuXuzMs2HROAspbPuefdw/aTEvfLOSlvVr8dcRXRh0VHLYYYlUOKUeZxGsJ5FuZjeb2StERlafCjwNDCr7UEXKXt2kRP48vCuvXXs8R9SM58rnMrhh3EzWbdsTdmgilUaxNQsz20Vk+o3HgE+DNosKTzULKcrevALGfv4tD3+cRc34OH49pBMX9WtDfJzWyxAprmYRLVlcCBxHpAdUPjAd+Br42t0r7MyvShYSzfINO/n9m/P4KmsjPVIbcNeIrnRNqR92WCKhKnWyOOAitYmsXHc8kR5RNdy9TZlFWYaULKQk3J23Zq/mb+8sYNPOvVx+fBq/PLUjdZMSww5NJBTFJYtiu84GJx8B9COSJPoDxxBZue6rsgxSpLyZGSN6pXBSp2b8fVJk+vN35q7h98M6c2b3FlrKVaSQaK+hZhFZeCgDmBJ8vnH3HeUTXumoZiGlMWfVFn7/5nzm5Wzl+HaN+fNZXeig2WylGjmcNovuwDw/lP61FYCShZRWfoHz0tSV3DdpMbv25nPlgLbcfHIH6tSMWgkXqfQOZ4ry7lEu3M7MBpQ6MpEKJj7OuPS4ND657UTO7d2KsZ8vY9D9n/LmrBwq2XcmkTIV7etSY2C2me1fIW89kAS0B04ANgC3xzRCkRA0rlOTv5/XnVF9U/nTxExGvzKbF79ZyZ1ndVGvKamWovaGMrN4IgPw+gMtgN1Eljd9z92/i3mEpaDXUFKWCgqcVzNWce+kxWzetZdRx7TmtsEdaVynZtihiZSpMuk6W4qbphJZJCkZcGCsuz9kZiOBO4Gjgb7unlHonDuAq4iM6bjZ3ScF5UOAh4B44Gl3v6e4eytZSCxs3bWPhz5aygtfr6BWjXhuObkDlx2XRo2Ekk7eLFKxhZUsWgAt3H2mmdUl8hprBJHEUQA8Cdy2P1mYWWfgZSJjOVoCHwIdg8stITLNSDaRgYEXuvuCou6tZCGxlLVuO39+ewFfLN3AkU2O4PfDjuakTs3U1VYqvbJYg/uQufsad58ZbG8n8uoqxd0Xuvvig5wyHBjv7rnBtCJZRBJHXyDL3Ze5+15gfHCsSCjaN6vLC1f25dmfR/5NXflcBpc9O43Fa7Wkq1Rd5VJ/NrM0oBcwtZjDUogM9tsvOygrqlwkNGbGoKOSeX/0QP4wrDNzVm3h9Ic+547X57F+e27Y4YmUuRJ1HjezmsC5QFrhc9z9LyU4tw4wARjt7jFb59LMrgauBmjdunWsbiPyIzUS4rhqQFvO6ZXCwx8v5V9fr2Ti7ByuO7EdVw04klo1KuQs/iKHrKQ1i7eIvPrJA3YW+hTLzBKJJIpx7v56lMNziIwW369VUFZU+Y+4+1h3T3f39KZNm0YLTaRMNTyiBn86swsf3DqQ/u2bcP8HSzjp/k95bUY2+QUanyGVX4kauM1svrt3PaQLR1r7ngc2ufvog+z/lB83cHcBXuKHBu6PgA6AEWngPplIkpgOXOTumUXdWw3cErZpyzdx1zsLmJO9laOa1+WOoUczsEMTNYJLhVYWDdxTzKzbId63P3ApMMjMZgefoWZ2tpllE5n6/B0zmwQQ/PJ/lcj6Ge8DN7h7vrvnATcCk4g0kr9aXKIQqQj6tm3EG9f355ELe7Fzbx6XPzuNS56ZyrzsrWGHJlIq0eaGmkekq2sCkW/5y4BcIt/23d2LnQ4kLKpZSEWyN6+AcVNX8sjHWWzauZcze7TktsEdadP4iLBDE/mRw5lIsNj1Ktx95WHGFhNKFlIRbd+zj7GfL+PpL5azL7+AC/u25qZB7WlWLyns0ESAkAblhUnJQiqyddv28NBHS3ll+ioS4+O4on8a1wxsR/3aWnRJwqVkIVIBrdiwkzGTlzBxzmrqJSVwzQntuKJ/GrVraDp0CYeShUgFtmD1NsZMXsyHC9fRpE4Nrj+xPRf1a01SosZoSPlSshCpBGas3Mz9kxbz9bKNNK+XxE0nt2dkn1RNVCjlRslCpBKZkrWB+z9YzMzvttCqYS1uHtSBs3unkBivpCGxpWQhUsm4O58uWc8Dk5cwN3srbRrX5qZBHRjRsyUJShoSI0oWIpWUu/PRwnU88OESMldvIy1IGsOVNCQGlCxEKjl3Z/KC73nww6UsWBNJGjec1J4RvfR6SsqOkoVIFbE/aTz00VIyV2+jdaPaXH9iO87p3UoN4XLYlCxEqpj9r6ce/ngpc7O3ktKgFteecCQj01PV5VZKTclCpIpydz5bsp5HPs5ixsrNNKtbk6sHHslF/VprcJ8cMiULkSrO3fn62408+kkWU77dSMPaiVzZvy2XHZemaUSkxJQsRKqRGSs389gnWXy8aB11aiZw8bGtuWpAW5rV1YSFUjwlC5FqKHP1Vh7/9FvenbeGhPg4RvZpxdUDj9TU6FIkJQuRamz5hp2M/XwZE2Zkk1dQwNBuLbj2hHZ0TakfdmhSwShZiAjrtu3hmS+XM27qd+zIzeNnHZpwzcB29G/fWMu9CqBkISKFbN29j5emfsezXy1n/fZcOreoxzUnHMnQbi00wK+aU7IQkZ/Ysy+ft2bnMPbzZXy7ficpDWpxRf80LjgmlbpJ6kFVHRWXLGL2NcLMUs3sEzNbYGaZZnZLUN7IzCab2dLgz4ZBuZnZw2aWZWZzzax3oWtdHhy/1Mwuj1XMItVJUmI8FxzTmsm3nsDTl6XTqmEt/vbOQo6/+2P+992F5GzZHR61IJQAAA9pSURBVHaIUoHErGZhZi2AFu4+08zqAjOAEcDPgU3ufo+Z3Q40dPffmNlQ4CZgKNAPeMjd+5lZIyADSAc8uE4fd99c1L1VsxApnTmrtvDUF8t4b/5aAIZ2a8EvBrSlR2qDkCOT8lBczSJmQzzdfQ2wJtjebmYLgRRgOHBicNjzwKfAb4LyFzySvb4xswZBwjkRmOzum4KHmQwMAV6OVewi1VWP1AY8elFvsjfv4vkpKxg/bRVvz1lNepuGXDWgLad2TtZst9VUufytm1ka0AuYCiQHiQRgLZAcbKcAqwqdlh2UFVV+4D2uNrMMM8tYv359mcYvUt20alib353RmSl3DOKPwzrz/fY9XDduJifc9ylPfb6Mrbv3hR2ilLOYJwszqwNMAEa7+7bC+4JaRJm8B3P3se6e7u7pTZs2LYtLilR7dZMSuXJAWz697SSeuKQPKQ1rcde7Cznu7o/4w5vzyVq3I+wQpZzEdKYxM0skkijGufvrQfH3ZtbC3dcEr5nWBeU5QGqh01sFZTn88Npqf/mnsYxbRH4sPs4Y0rU5Q7o2Z37OVp6bsoJXpq/iX9+sZGDHpvz8+Dac2LEZcXEar1FVxbKB24i0SWxy99GFyu8DNhZq4G7k7r82szOAG/mhgfthd+8bNHDPAPb3jppJpIF7U1H3VgO3SOxt2JHLy1O/48WpK/l+Wy5tGtfm0mPbMDI9lfq11PW2MgplnIWZDQC+AOYBBUHxb4m0W7wKtAZWAue7+6YguTxKpPF6F3CFu2cE17oyOBfgLnf/Z3H3VrIQKT/78gt4f/5anp+ygoyVm6mVGM+IXilcdlwbjm5RL+zw5BBoUJ6IlIv5OVt54esVvDV7Nbl5BRyT1pBLj0tjSJfmWsmvElCyEJFytWXXXv6dkc2/vlnJd5t20aROTS7sm8qFfVvTskGtsMOTIihZiEgoCgqcz5eu58VvVvLRonUYMOioZlx8bBsGdmhKvBrEK5RQBuWJiMTFGSd2asaJnZqxatMuxk//jlemr+LDheto1bAWF/ZtzfnpqTStWzPsUCUK1SxEpFztzStgUuZaXpr6HV8v20hCnDG4SzIX9W3D8e0aq/ttiFSzEJEKo0ZCHGf2aMmZPVqStW4HL0/7jgkzs3l33lpaN6rNqL6pnNenlZaBrWBUsxCR0O3Zl8+kzLWMm/od05ZvIiHOOPnoZlzYtzU/U9tGuVEDt4hUGt+u38Er01fx2oxsNu3cS8v6SYxMT2VkeitaNawddnhVmpKFiFQ6uXn5fLhgHeOnf8eXWRsAGNC+CeenpzK4SzI1E+JDjrDqUbIQkUote/Mu/p2RzWszssnZspsGtRMZ0TOFkemt6NKyftjhVRlKFiJSJeQXOFO+3cAr01fxQeb37M0voHOLeoxMb8Xwnik0OqJG2CFWakoWIlLlbNm1l4lzVvNqxirm52wjMd445ehkzuvTioEdm5KoRZoOmZKFiFRpC9ds47UZ2bw5K4eNO/fSpE5NRvRsybl9Wmkyw0OgZCEi1cK+/AI+Xbye12as4uNF69iX73RuUY9zeqcwvGeKRopHoWQhItXOpp17eXvOaibMzGZu9lbi44yBHZpwdu9WDO6cTFKielMdSMlCRKq1pd9v5/VZObw5K4c1W/dQp2YCp3dtztm9Uzi2raYY2U/JQkSEyCy43yzbyBuzcnhv/lp25ObRon4SZ/VsyYieKdW+fUPJQkTkALv35jN54fe8NSuHz5asJ6/A6ZRcl+G9WjK8Zwop1XDdDSULEZFibNyRyzvz1vDW7NXMWLkZgGPSGnJWzxSGdm1O4zrVo2E8rDW4nwWGAevcvWtQ1gN4AqgDrAAudvdtwb47gKuAfOBmd58UlA8BHgLigafd/Z5o91ayEJHSWrVpF2/NzuGt2atZum4H8XHGgPZNOKtHSwZ3SaZuUmLYIcZMWMliILADeKFQspgO3Obun5nZlUBbd/+DmXUGXgb6Ai2BD4GOwaWWAKcC2cB04EJ3X1DcvZUsRORwuTuL1m7nrdmreXvOanK27KZmQhyDjmrGmT1aMuioZlWuR1Uo61m4++dmlnZAcUfg82B7MjAJ+AMwHBjv7rnAcjPLIpI4ALLcfRmAmY0Pji02WYiIHC4z4+gW9Ti6RT1+fVonZq3azMTZq3ln3lrem7+WI2rEc0rnZIZ1b8nAjk2q/MSG5b34USaRX/ZvAiOB1KA8Bfim0HHZQRnAqgPK+x3swmZ2NXA1QOvWrcsuYhGp9uLijD5tGtGnTSP+eGYXpi7byNtzV/Pe/LW8NXs1dZMSGNy5OcN6tKB/uybUSKh6U42Ud7K4EnjYzP4ATAT2ltWF3X0sMBYir6HK6roiIoXFxxnHt2/C8e2b8JfhXfkyawP/mbOGDxasZcLMbOrXSuS0LskM7daC/u2bVJk5qso1Wbj7ImAwgJl1BM4IduXwQy0DoFVQRjHlIiKhSoyP46ROzTipUzNy87ry5dINvDN3De/OW8urGdk0qJ3I4M7JnN6t8tc4yjVZmFkzd19nZnHA74n0jIJILeMlMxtDpIG7AzANMKCDmbUlkiRGAReVZ8wiIiVRMyGek49O5uSjk8nNy+eLJRt4Z94a3gsSR72kBE7t3Jyh3ZozoEPla+OIWbIws5eBE4EmZpYN/AmoY2Y3BIe8DvwTwN0zzexVIg3XecAN7p4fXOdGIg3h8cCz7p4Zq5hFRMpCzYRI4/cpnSOJ48ulG3h33tr/vqqqUzOBk49uxuldm3NCx2bUqlHxE4cG5YmIlJO9eQVM+XYD7wWJY/OufdRKjOfETk0Z0rU5g45qFuo4Do3gFhGpYPLyC5i2fBPvzV/LpMy1rNueS434OPq3b8xpXZpzSudkmpTzyHElCxGRCqygwJm1aguTMtfy3vw1rNq0mziD9DaNGNwlmdO6NCe1Ue2Yx6FkISJSSbg7C9dsZ1JmpMaxaO12AI5uUY/BnZM5tXMyXVrWw6zsp1VXshARqaRWbtzJ5AXfMylzLRkrN+MOKQ1qcWrnZAZ3TuaYto3KbCyHkoWISBWwYUcuHy9cxwcL1vLF0g3k5hVQLymBQUc145TOyZzQselhNZArWYiIVDG79ubxxdINTF7wPR8vWsemnXtJjDdO69KcRy/qXaprhjKRoIiIxE7tGgmc1qU5p3VpTn6BM+u7zUxe+D0JMVoiVslCRKSSi48z0tMakZ7WKGb3qLwTlYiISLlRshARkaiULEREJColCxERiUrJQkREolKyEBGRqJQsREQkKiULERGJqkpO92Fm64GVh3GJJsCGMgqnotGzVV5V+fn0bBVDG3dverAdVTJZHC4zyyhqfpTKTs9WeVXl59OzVXx6DSUiIlEpWYiISFRKFgc3NuwAYkjPVnlV5efTs1VwarMQEZGoVLMQEZGolCxERCQqJYsDmFm8mc0ys/+EHUtZM7MGZvaamS0ys4VmdlzYMZUVM7vVzDLNbL6ZvWxmSWHHVFpm9qyZrTOz+YXKGpnZZDNbGvzZMMwYD0cRz3df8P/lXDN7w8wahBljaR3s2Qrt+5WZuZk1CSO2w6Vk8VO3AAvDDiJGHgLed/ejgB5Ukec0sxTgZiDd3bsC8cCocKM6LM8BQw4oux34yN07AB8FP1dWz/HT55sMdHX37sAS4I7yDqqMPMdPnw0zSwUGA9+Vd0BlRcmiEDNrBZwBPB12LGXNzOoDA4FnANx9r7tvCTeqMpUA1DKzBKA2sDrkeErN3T8HNh1QPBx4Pth+HhhRrkGVoYM9n7t/4O55wY/fAK3KPbAyUMTfHcADwK+BStujSMnixx4k8hdaEHYgMdAWWA/8M3jN9rSZHRF2UGXB3XOA+4l8a1sDbHX3D8KNqswlu/uaYHstkBxmMDF2JfBe2EGUFTMbDuS4+5ywYzkcShYBMxsGrHP3GWHHEiMJQG/gcXfvBeykcr/K+K/g/f1wIgmxJXCEmV0SblSx45H+7pX2G2pxzOx3QB4wLuxYyoKZ1QZ+C/wx7FgOl5LFD/oDZ5nZCmA8MMjMXgw3pDKVDWS7+9Tg59eIJI+q4BRgubuvd/d9wOvA8SHHVNa+N7MWAMGf60KOp8yZ2c+BYcDFXnUGgLUj8iVmTvC7pRUw08yahxpVKShZBNz9Dndv5e5pRBpHP3b3KvPt1N3XAqvMrFNQdDKwIMSQytJ3wLFmVtvMjMizVYnG+0ImApcH25cDb4UYS5kzsyFEXgGf5e67wo6nrLj7PHdv5u5pwe+WbKB38O+xUlGyqF5uAsaZ2VygJ/C/IcdTJoLa0mvATGAekf+vK+0UC2b2MvA10MnMss3sKuAe4FQzW0qkJnVPmDEejiKe71GgLjDZzGab2ROhBllKRTxblaDpPkREJCrVLEREJColCxERiUrJQkREolKyEBGRqJQsREQkKiULqZbMrHHQRXO2ma01s5xCP9c44NgTD2cWYjN7zszOO/yoi7z+iWZ2fKGfY3o/qZ4Swg5AJAzuvpHIWBPM7E5gh7vfH2pQpXcisAOYEnIcUoWpZiESOPAbuZntKLS7npm9Y2aLzewJM4sLjhlsZl+b2Uwz+7eZ1SnhveKDNRymB2s4XBOUn2hmnxZad2RcMCodMxsalM0ws4fN7D9mlgZcC9wa1Ip+FtxioJlNMbNlqmVIWVCyECmZvkRGwHcmMt/POcEiNr8HTnH33kAG8MsSXu8qIrPjHgMcA/yPmbUN9vUCRgf3OhLoHyzm9CRwurv3AZoCuPsK4AngAXfv6e5fBNdoAQwgMtdSpR3tLRWHXkOJlMw0d18G/53SYQCwh8gv9K+CL/81iEz1UBKDge6FvvXXBzoAe4N7ZQf3mg2kEXnNtMzdlwfHvwxcXcz133T3AmCBmVXl6cylnChZiPwgj6C2HbxmKtzQfeC8OA4YMNndLyzFvQy4yd0n/ajQ7EQgt1BRPqX7d1r4GlaK80V+RK+hRH6wAugTbJ8FJBba19fM2gZJ5ALgSyIruvU3s/YAZnaEmXUs4b0mAdeZWWJwbscoi1EtBo4M2igIYthvO5FJ+ERiRslC5AdPASeY2RzgOCILRO03ncjMqAuB5cAb7r4e+DnwcjCT79fAUUVc+8lgFtJsM/uayNK9C4isbTCfSHtEkTUId98NXA+8b2YziCSIrcHut4GzD2jgFilTmnVWpJIwszruviPoHfUYsNTdHwg7LqkeVLMQqTz+J2jwziTSIP5kyPFINaKahYiIRKWahYiIRKVkISIiUSlZiIhIVEoWIiISlZKFiIhE9f8B7SM7lzciLx8AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"cell_type": "code",
"source": [
"R = 0.5*2.54/100 #m\n",
"L = 9.5*2.54/100 #m\n",
"A = L*2*np.pi*R\n",
"print(\"%.2f W LNTF TS at 1.5 Bar\"%(h(L)*A*(Tsat-Twall)))"
],
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/"
},
"id": "2RFON3jE9msO",
"outputId": "74404941-1391-4927-b15c-1828d8ea4815"
},
"execution_count": 12,
"outputs": [
{
"output_type": "stream",
"name": "stdout",
"text": [
"146.61 W LNTF TS at 1.5 Bar\n"
]
}
]
},
{
"cell_type": "code",
"source": [
"from requests import get\n",
"import pandas as pd\n",
"from bs4 import BeautifulSoup\n",
"url=\"https://webbook.nist.gov/cgi/fluid.cgi?Action=Load&ID=C7440633&Type=SatP&Digits=5&THigh=180&TLow=161.4&TInc=.1&RefState=DEF&TUnit=K&PUnit=bar&DUnit=mol%2Fl&HUnit=kJ%2Fmol&WUnit=m%2Fs&VisUnit=uPa*s&STUnit=N%2Fm\"\n",
"raw_html = get(url).content\n",
"html = BeautifulSoup(raw_html, 'html.parser')\n",
"table = html.findAll(\"table\")\n",
"xe_data = pd.read_html(table[1].prettify())[0]\n",
"#print(xe_data.columns)\n",
"raw_html = get('https://webbook.nist.gov/cgi/cbook.cgi?ID=C7727379&Units=SI&Mask=4#Thermo-Phase').content\n",
"n2_html = BeautifulSoup(raw_html, 'html.parser')\n",
"table = n2_html.find(\"table\", attrs={'aria-label': 'Antoine Equation Parameters'})\n",
"n2_data = pd.read_html(table.prettify())[0]"
],
"metadata": {
"id": "hDcozeb-Bqbn"
},
"execution_count": 13,
"outputs": []
},
{
"cell_type": "code",
"source": [
"Twall = 77.34 #Ln2\n",
"R = 0.5*2.54/100 #m\n",
"L = 9.5*2.54/100 #m\n",
"A = L*2*np.pi*R\n",
"T = np.linspace(78,100,100)\n",
"n2_P = np.power(10,[n2_data[\"A\"][0] - (n2_data[\"B\"][0]/(i + n2_data[\"C\"][0])) for i in T])\n",
"plt.plot(T,n2_P)\n",
"plt.ylabel('Pressure (Bar)')\n",
"plt.xlabel('Temperature (K)')\n",
"plt.title(\"N2 Vapor Pressure NIST via Antoine Equation\");"
],
"metadata": {
"id": "PciGT4ZQBsdl",
"colab": {
"base_uri": "https://localhost:8080/",
"height": 295
},
"outputId": "a44d39e1-c7eb-4e60-c330-6fcc4407e675"
},
"execution_count": 14,
"outputs": [
{
"output_type": "display_data",
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAEWCAYAAABsY4yMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3dd5wU9f3H8deH3vtRpHdpUjwRVBTFFrvGHkUsoMFYYhLTo0mMMbEk+ZlERZEigihBY0WNFRCBo0jv/ahHOerBlc/vj5lL1uPKAre3d3vv5+Oxj9ud9v3M3Oxnv/Odme+YuyMiIomnQrwDEBGR2FCCFxFJUErwIiIJSgleRCRBKcGLiCQoJXgRkQSlBC9SgsxsgJktj3cc0TCz583s1/GOI5bM7Htm9lG844gVJfg8zGydmW03s5oRw+4ys8/D943NbIKZbTazdDObbmanF7CsG8PlWZ7hlcIyLovpyhQhjO2Qme03s21mNtrMasUzpuJiZo+amZvZ9RHDKoXD2oSfR5vZYxHj7zSzZWa2L9we75tZbTP7INxG+80s08yORHx+/ljicvep7t75BNdtSLgeNxzjfN9a36K4+z3u/vtjj7DIOIaYWXbENsx9nVTcZeUpt0243SrlDnP3V939wliWG09K8PmrCDxQwLhawGzgVKABMAZ4r4DE+BZQDzgnz/CLAQemFEu0RbBAQf/ry929FtAHSAZ+lc/8lY6aq4QdZwy7gN+aWcUoln8O8Dhwk7vXBroAEwHc/TvuXivcTq8Cf8797O73HEdcJ+o2gnUbHIeyi8uMiG2Y+9oc76ASjRJ8/p4Efmxm9fKOcPc17v6Mu29x92x3HwFUAY6qlbl7BvA6R38RBwPj3T3LzN4ws63h0cCXZtYtd6KwxvW8mX0c1iq/MLPWEePPMLPZ4byzzeyMiHGfm9kfzGw6cBBoV9gKu3sq8AHQPZzfzexeM1sJrAyHXWZm881sj5l9ZWanRJT3UzNLDeNcbmaDwuF9zSzFzPaGteJnwuEDzWxTZAzhEcX54ftHzWySmY0zs73AEDOra2YjzWxLWNZjRSTvKcAR4JbC1j10GkHSmRduj13uPsbd90Uxb+Q6VA23T/eIYUnhkVLjvOttZj8zs9XhdltiZlcXsfzWBBWGYcBFZtY0YtxAM9tkZj8KjxC3mNnt4bhhwPeAh8Pa8jvh8C7hvrLHzBab2RURy/tvjb+wZUes91NmtiH8Pz9vZtWPZdtFLKu3mc0Nt8lEM3stIo4hZjYtz/RuZh3C95ea2bxwf9toZo9GTPpl+HdPuA36511eFN+p31tw1L7PzD4ys0bHs44lRQk+fynA58CPi5rQzHoRJPhVBUwyBrg2d2c3s7rA5eFwCJJqR6AxMJeghhjpe8DvgUbA/NzxZtYAeA/4P6Ah8AzBkUTDiHlvJUgEtYH1RaxHS+ASYF7E4KuA04GuZtYbeBm4OyzvBeDt8IvdGfgBcFpY+70IWBcu42/A39y9DtCe4AcvWlcCkwiOgl4FRgNZQAegN3AhcFch8zvwa+ARM6tcRFkzCRLmb83sTDOregxx/q9A98PAZOCmiMHXA1+4+/Z8ZlkNDADqAr8FxplZs0KKGAykuPu/gKUE+0ekpuGymgN3Av8ws/phRSTy6OPycJu8A3xEsP/dB7wa/j/zk++yw3FPAJ2AXgT/n+bAbwpZj3yZWRWCI99XCI6Q3wC+ewyLOECwjeoBlwLfN7OrwnFnh3/rhdtgRp6yo/lO3QzcTrC9qhBFjognJfiC/Qa4z8ySCprAzOoQ7Ii/dff0/KZx9+nANiC3ZnY9sMLd54fjX3b3fWFieBToGf4I5HrP3b8Mx/8S6B8m40uBle7+irtnufsEYBnBj0eu0e6+OByfWcBqvGVme4BpwBcEzRS5/hjWZA8R/FC84O4zwyOXMcBhoB+QDVQl+CGo7O7r3H11uIxMoIOZNXL3/e7+dUHbMx8z3P0td88B6hD8AD3o7gfCZPkX4MbCFuDubwM7KPyHAHefClxD0FT1HrDTzJ4p4gihIOPzxHVzOCy/ct9w983unuPuEwmOlvoWsuzBEcsaz9FHh5nA79w9093fB/aTz9FlqB9Bk+MT7n7E3T8F3uXbP05FLtvMjGD/+GG4v+wj2I8K+9/0C48acl+5+0s/oDLw17CcSQRNolFx98/dfWG4PRcAEzi6ibQg0XynRrn7ivA78TrBD1qppQRfAHdfRLCz/yy/8WGN/B3ga3f/YxGLG8v/voi3hp8xs4pm9kR4iL6X/9V6Iw/7NkbEtJ+g7fWk8JW3Vr6eoOZ01LyFuMrd67l7a3cfHu64+c3fGvhR5JcSaAmc5O6rgAcJfqC2h4fUuSfM7iSo2S0LD3mP5cRy3vIrA1siyn+BoCZVlF8R/DhWK2wid//A3S8nqDleCQyhiB+GAnwG1DCz0y04odsLeDO/Cc1ssP2v2WsPQRNZvof9ZnYm0BZ4LRw0HugRHkXm2unuWRGfDxIk8fycBGwMf0Bz5d2HIhW07CSgBjAnYj2mhMML8nW43+W+2kfElOrf7gWx0KPPSOE2/8zMdphZOnAPBWzPfETzndoa8b6wbVsqKMEX7hFgKHl2+PDw/S1gE0GTRVFeAQaZWX+CGkpuM8zNBInkfIJD3za5RUTM2zKi3FoEyWdz+GrNt7UCUiM+n2hXoZHzbwT+kOdLWSOs5eDu4939rDAmB/4UDl/p7jcRJOI/AZMsuELpAEFSyF23ihydEPKWfxhoFFF+HXfvRhHc/WOCJrThUa10UPv7BPiU8JzEsXD3bILa3U3h69382vItaE9/kaB5q6G71wMW8e3/f6TbwnHzzWwrQbNS7vCoQsvzeTPQ0r59Aj7vPhSNNOAQ0C3if1M3PCl9rLYAzcOjgsiYcuXdb5rybeOBt4GW7l4XeJ7/bc+ivg/RfKfKFCX4QoQ104nA/bnDwnbLSQQ79G15aj8FLWcdQRPIBOBjd8+tBdQmSFo7CXbax/OZ/RIzOytsm/w9Qc1nI/A+0MnMbrbg8r8bgK4ERx2x8CJwT1hDMjOrGZ7Qqm1mnc3svPCHL4Ng2+QAmNktZpYUbqc94bJygBVAtXAZlQlq2QW2e7v7FoK24qfNrI6ZVTCz9hZc/RKNXwIPFzTSzK604LLW+uH69SU4tD+WJqVI44EbCNrI822eAWoSJJ0dYQy3U8APiplVI2jeG0ZwRJD7ug+42aK7ymgb3z7ZPpOgFvqwmVU2s4EEzRGv5TNvgcL/7YvAX8yscRhvczO76FiWE5pBcJ7l/jCma/h2k9U3QDcz6xVuk0fzzF8b2OXuGeH/8OaIcTsI9r2CLjgo6e9UzCnBF+13BF/EXGcAlxGc4Ms9G7/fzAYUsZwxBLWDsRHDxhIcAqYCS8g/mYwnOJLYRXBp5i0A7r4zjONHBD8QDwOXuXvaMa1dlNw9heBo5u/AboIa8ZBwdFWCk2xpBIewjYGfh+MuBhab2X6CE643uvuh8JzFcOAlgvU/QHBEVJjBBCe2loQxTAIKOyEZGf90YFYhk+wO128lsBcYBzzp7nlPekfF3WcSrNNJBCfS85tmCfA0QVLbBvQAphewyKsIfjjHuvvW3BfBie9KBNu5KCMJzpPsMbO33P0IQUL/DsH/7p/AYHdfFuVqRvopwT7xddjc+B8KbvuH4FxS3uvgTwtjuoZg39pF8CM5OXcmd19B8J38D8H/alqe5Q4Hfmdm+wjOo70eMe9B4A/A9HAb9IucsaS/UyXBXA/8KLXMbDSwyd2PujZdpLzQ9+D4qQYvIpKglOBFRBKUmmhERBKUavAiIgkq7p1IRWrUqJG3adMm3mGIiJQZc+bMSXP3fG8qK1UJvk2bNqSkpMQ7DBGRMsPMCrzTV000IiIJSgleRCRBxTTBm9kPLehjepEFT0EqtLMnEREpPjFL8GbWnKAPl2R3707wlKRCu3YVEZHiE+smmkpA9bAjpBoEvbWJiEgJiFmC9+ARcE8BGwi6AE1396OeXm5mwyx4pFvKjh07YhWOiEi5E8smmvoEfZ23JehRr6aZHfVsTHcf4e7J7p6clFTY8wFERORYxLKJ5nxgrbvvCB8XN5mgq10REQmlrNvFS1PXEItuY2KZ4DcQPHexRvh0lkEEDwkWERFgbdoBho5NYdzX6zl4JLvYlx/LNviZBA9kmAssDMsaEavyRETKkp37DzNk1CzMjNG396Vm1eLvWCCmXRW4+yMETyMSEZFQRmY2Q8emsDU9g/FD+9GmUc2iZzoOpaovGhGRRJeT4zz0+nzmbdzDP2/uw6mt68esLHVVICJSgv7w/lLeX7iVX17She/0iOqRwsdNCV5EpISMnLaWkdPWMuSMNtx5VtuYl6cELyJSAj5YuIXH3lvCRd2a8OvLuhJcXBhbSvAiIjE2a+0uHpg4n94t6/G3G3tTsULskzsowYuIxNTKbfu4a8xsWtSvzsjbTqNa5YolVrYSvIhIjGxNz+C2l2dRtXJFxtzel/o1q5Ro+UrwIiIxsDcjkyGjZpF+KJNRQ06jZYMaJR6DroMXESlmh7OyGTY2hVXb9zPq9tPo3rxuXOJQghcRKUY5Oc5DE7/h6zW7+OsNvRjQMX695KqJRkSkmLg7v3t3Ce8t3MIvL+nCVb2bxzUeJXgRkWLyz89XM/qrddx1VluGnt0u3uEowYuIFIfXZm3gyQ+Xc3Xv5vziki7xDgdQghcROWEfLt7KL95cyDmdkvjztadQoYRuZCqKEryIyAn4es1O7pswjx4t6vHcLX2oXLH0pNXSE4mISBmzKDWdoWNSaNWgBqOGnEaNKqXrwkQleBGR47A27QBDRs2idrVKjL2jLw1K+C7VaCjBi4gco63pGdw6cibZOc7YO0/npHrV4x1SvpTgRUSOwe4DR7h15Ex2HzjC6Nv70qFxrXiHVKDS1WAkIlKKHTicxZDRs1m/6yCjbz+Nni3rxTukQsWsBm9mnc1sfsRrr5k9GKvyRERi6XBWNne/ModFqen8/abenNG+UbxDKlLMavDuvhzoBWBmFYFU4M1YlSciEitZ2TncN34e01al8fR1PbmwW9N4hxSVkmqDHwSsdvf1JVSeiEixyMlxHp60gI+WbOO3V3Tju6e2iHdIUSupBH8jMCG/EWY2zMxSzCxlx44dJRSOiEjR3J1H31nM5Hmp/PjCTtx2Rpt4h3RMYp7gzawKcAXwRn7j3X2Euye7e3JSUvy61RQRyevPHy5n7Iz1DDu7Hfee2yHe4RyzkqjBfweY6+7bSqAsEZFi8Y/PVvHc56u5+fRW/Pw7J2NWOvqXORYlkeBvooDmGRGR0mjU9LX/7RnysSu7l8nkDjFO8GZWE7gAmBzLckREisvE2Rv47TtLuKhbE54sRT1DHo+Y3ujk7geAhrEsQ0SkuLw1L5WfTQ66/f2/m3pTqRT1DHk8ynb0IiLF5IOFW/jRG9/Qr21DXrj1VKpWqhjvkE6YEryIlHufLtvG/a/No2eLurx0WzLVKpf95A5K8CJSzn2xYgf3vDKXk5vWYfQdfalZNXG66FKCF5Fy66tVaQwbm0L7xrV45c6+1KlWOd4hFSsleBEpl2at3cWdY1Jo3bAG4+7sS70ape+BHSdKCV5Eyp0563dx+6hZNKtXjXF3nU7DWlXjHVJMKMGLSLkyb8Nubnt5No3rVGPC0H40rl0t3iHFjBK8iJQb32zcw+CRs2hYqwoThvajSZ3ETe6gBC8i5cSCTXu4ZeRM6tWszPih/WhaN7GTOyjBi0g5sHBTOre8NJO61SszYWg/mpfSh2QXNyV4EUloi1LTuWXkTOpUr8xrw/rRon6NeIdUYpTgRSRhLdyUzs0vfk2tqpWYMLR8JXeIcWdjIiLx8s3GPdxaTmvuuZTgRSThzA+Te70alXltWP9y0+ael5poRCShzFm/m1tfmkn9GlWYWI6TOyjBi0gCmbV2F4NHzqRR7apMvLsfJ5Xj5A5K8CKSIGas3sltL8+iad1qvDasH83qlu/kDkrwIpIAvlyxgyGjZtGifnUmDEv8O1SjpZOsIlKmfbJ0G98fN5f2jWsx7s6+Cdtx2PFQgheRMmvKoi3cN2EeJzetwysJ2uXviYhpE42Z1TOzSWa2zMyWmln/WJYnIuXHv+encu/4eXRvXpdXh56u5J6PWNfg/wZMcfdrzawKUP7uNBCRYvfarA38/M2FnN62AS/ddhq1Eugxe8UpZlvFzOoCZwNDANz9CHAkVuWJSPkwevpaHn1nCed0SuKFW09NmAdkx0Ism2jaAjuAUWY2z8xeMrOaeScys2FmlmJmKTt27IhhOCJS1v3js1U8+s4SLuzahBGDldyLEssEXwnoAzzn7r2BA8DP8k7k7iPcPdndk5OSkmIYjoiUVe7OEx8s48kPl3N17+b883t9qFpJyb0osUzwm4BN7j4z/DyJIOGLiEQtJ8f5zb8X8/wXq/ne6a14+rqeVKqoW3iiEbOt5O5bgY1m1jkcNAhYEqvyRCTxZGbn8NDr83nl6/XcfU47HruqOxUqWLzDKjNifer5PuDV8AqaNcDtMS5PRBJERmY2Pxg/l/8s3c7DF3dm+MAO8Q6pzIlpgnf3+UByLMsQkcSzLyOTYWPn8PXanfz+qu7c2q91vEMqk3TxqIiUKjv3H2bIqNks3bKXv1zfi6t6N493SGWWEryIlBqpew5x60sz2Zx+iBcHJ3PuyY3jHVKZpgQvIqXCym37GPzyLPYfzuKVO0/ntDYN4h1SmacELyJxN3fDbu4YPZvKFSswcVh/up5UJ94hJQQleBGJq8+Wb2f4uLk0rlOVV+44nVYN1WVVcSkywYc9QN4CDACaAYeARcB7wDh3T49phCKSsN6ct4mfvLGATk1qM+aOviTVVl/uxanQG53M7APgLuBD4GKCBN8V+BVQDfi3mV0R6yBFJPGM+HI1P5z4Dae1acBrd/dTco+Bomrwt7p7Wp5h+4G54etpM2sUk8hEJCHl5DiPv7+Ul6at5dJTmvHM9T3Vr0yMFFqDd/c0M6toZp8VNk3xhyUiiehwVjYPTJzPS9PWMuSMNjx7Y28l9xgqsg3e3bPNLMfM6qq9XUSO196MTO4eO4cZa3by04tP5p5z2mGmfmViKdqraPYDC83sY4JufwFw9/tjEpWIJJSt6RkMGTWLVdv388z1PbmmT4t4h1QuRJvgJ4cvEZFjsnzrPoaMmsXeQ5mMuv00BnTUcx9KSlQJ3t3HxDoQEUk8X61O4+5X5lC9ckVev6c/3U6qG++QypWoEryZdQT+SHCJZLXc4e7eLkZxiUgZ99a8VH4y6RvaNKzJ6Dv60rxe9XiHVO5E+8CPUcBzQBZwLjAWGBeroESk7HJ3nv1kJQ9OnE+fVvWZdM8ZSu5xEm2Cr+7unwDm7uvd/VHg0tiFJSJlUWZ2Dj/710Ke/ngFV/duztg7+1K3RuV4h1VuRXuS9bCZVQBWmtkPgFSgVuzCEpGyZm9GJsPHzWXaqjTuH9SRH57fUZdBxlm0Cf4BoAZwP/B74DzgtlgFJSJly8ZdB7lj9GzW7TzAk9eewnXJLeMdkhD9VTSzw7f70XNVRSTCvA27GTo2hSNZOYy5oy9ntFfvJaVFUZ2NNTKzR8zsfjOrZWbPmdkiM/u3mekJuCLl3LsLNnPjiK+pXqUik4efqeReyhR1knU8UBXoCMwC1gDXAu8CLxW1cDNbZ2YLzWy+maWcaLAiUjrkXinzg/Hz6NG8Lm8NP5MOjXVarrQpqommibv/woIzJevd/clw+DIzuzfKMs5Vh2QiiSMjM5ufT17Im/NSubp3c574bg91GFZKFZXgswHc3c0sb5LOiU1IIlJa7dh3mLtfSWHuhj08dEEn7juvg66UKcWKSvDtzOxtwCLeE35uG8XyHfjIzBx4wd1H5J3AzIYBwwBatWoVdeAiUrKWbtnLXWNS2HngMP/8Xh8u6dEs3iFJEczdCx5pdk5hM7v7F4Uu3Ky5u6eaWWPgY+A+d/+yoOmTk5M9JUVN9SKlzYeLt/LDifOpXa0SLw0+jR4t1KdMaWFmc9w9Ob9xhdbgi0rgRXH31PDvdjN7E+gLFJjgRaR0cXf+8dkqnvpoBT1b1GXE4GSa1KlW9IxSKhR1meQ7Zna5mR11r7GZtTOz35nZHQXMW9PMaue+By4keFi3iJQBh45kc/9r83nqoxVc2eskJt7dX8m9jCmqDX4o8BDwVzPbBewg6E2yDbAa+Lu7/7uAeZsAb4YnYCoB4919SnEELSKxtXnPIYa9ksLizXv5yUWdGT6wvU6mlkFFNdFsBR4GHjazNkAz4BCwwt0PFjHvGqBn8YQpIiUlZd0u7hk3h4zMHF4anMygLk3iHZIcp2j7osHd1wHrYhaJiMTdqzPX8+jbi2lerzqvDUumQ+Pa8Q5JTkDUCV5EEteRrBweeXsxE2ZtYGDnJP52Q29185sAlOBFyrntezP4/qtzmbN+N8MHtudHF3amYgW1tyeCqBO8mVUHWrn78hjGIyIlaM76XXx/3Fz2ZWTx95t7c9kpJ8U7JClGUT3RycwuB+YDU8LPvSLuahWRMsbdGff1+v/2BPnmvWcouSegaGvwjxLcpPQ5gLvPN7NouioQkVImIzObX721iElzNqm9PcFFm+Az3T09z3WwBfdxICKl0sZdB7ln3BwWb97L/YM68uCgjlRQe3vCijbBLzazm4GKZtaR4NF9X8UuLBEpbp8t386Dr83H3Rl5m65vLw+iaoMH7gO6AYcJHgKSDjwYq6BEpPhk5zjPfLyCO0bPplndarz9g7OU3MuJImvwZlYReM/dzwV+GfuQRKS47DpwhAdem8fUlWl8t08LHruqO9Wr6OEc5UWRCd7ds80sx8zqunt6SQQlIiduzvrd/GD8XHbuP8LjV/fgpr4t1Z9MORNtG/x+YKGZfQwcyB3o7vfHJCoROW7uzqjp63j8/aU0q1eNycPPoHtz9d9eHkWb4CeHLxEpxdIPZfLTSQuYsngr53dpwtPX9dQlkOVYVAne3cfEOhAROTELN6Vz7/i5bN5ziF9e0oW7BrRVk0w5F1WCN7O15HPdu7u3K/aIROSYuDtjZ6znD+8tpVGtKky8uz+ntq4f77CkFIi2iSbyeX/VgOuABsUfjogci8gmmfNObszT1/Wkfs0q8Q5LSolom2h25hn0VzObA/ym+EMSkWjM27Cb+ybMY2t6Br+6tAt3nqUmGfm2aJto+kR8rEBQo1dXwyJxkJPjjJi6hqc+XE6TOtV4457+9G6lJhk5WrRJ+umI91kET3a6vtijEZFC7dh3mIden8/UlWlc2qMZj1/Tg7rVdZWM5C/aJppzj7eA8E7YFCDV3S873uWIlHefL9/Oj9/4hn0ZWbpxSaISbX/wD5hZHQu8ZGZzzezCKMt4AFh6/CGKlG+Hs7L53TtLGDJqNo1qVeWd+87i5tNbKblLkaLtbOwOd98LXAg0BG4FnihqJjNrAVwKvHTcEYqUYyu37eOqf3zFy9PXclv/1rx175l0aqIHYUt0om2Dz60qXAKMdffFFl314a/Aw4D2SJFj4O6Mm7mBx95dQs2qlXhxcDIXdFUPkHJsok3wc8zsI6At8HMzqw3kFDaDmV0GbHf3OWY2sJDphgHDAFq1ahVlOCKJK23/YX72rwX8Z+l2BnRsxNPX9aRxnWrxDkvKIHMv+sFMZlYB6AWscfc9ZtYAaOHuCwqZ548ETTlZBDdH1QEmu/stBc2TnJzsKSkpx7gKIonj02XbeHjSAvZmZPHwRZ2548y2euKSFMrM5rh7cn7joq3B9wfmu/sBM7sF6AP8rbAZ3P3nwM/DAAYCPy4suYuUZwePZPH4+0sZ9/UGTm5am3F3nc7JTevEOywp46JN8M8BPc2sJ/AjgpOmY4FzYhWYSHkxb8NuHnr9G9btPMDQAW358UWdqVpJD+WQExdtgs9ydzezK4G/u/tIM7sz2kLc/XPg8+OITyRhZWbn8OwnK/nH56tpWqca4+/qR//2DeMdliSQaBP8PjP7OUGb+oCwTV63z4kcpxXb9vHQ6/NZlLqXa/o059ErulGnmr5SUryiTfA3ADcTXA+/1cxaAU/GLiyRxJSd44yctoanPlpB7aqVeP6WU7m4e9N4hyUJKtquCraa2b+AjuGgNODNmEUlkoDWph3gx298w5z1u7mgaxP+eE0PGtWqGu+wJIFF25vkUIJr1RsA7YHmwPPAoNiFJpIYcnKcsTPW8cSUZVSpWIG/3NCTq3o1V1cDEnPRNtHcC/QFZgK4+0ozaxyzqEQSxLq0Azz8rwXMWruLgZ2TeOKaU2haVzctScmINsEfdvcjuTUOM6tEPo/wE5FATo4zZsY6/jRlGZUrVuDP157Cdae2UK1dSlS0Cf4LM/sFUN3MLgCGA+/ELiyRsmvNjv08PGkBKet3M7BzEn+8pgfN6laPd1hSDkWb4H8K3AUsBO4G3kc9RIp8S1Z2DiOnreWZj1dQtVIFnr6uJ9f0UVu7xE+RCT58YMdidz8ZeDH2IYmUPUs27+Wn/1rAwtR0LuzahMeu6q4OwiTuikzw7p5tZsvNrJW7byiJoETKiozMbJ79dCUvfLGGejUq84+b+3BJj6aqtUupEG0TTX1gsZnNAg7kDnT3K2ISlUgZ8PWanfxi8kLWpB3gu31a8KtLu1C/ZpV4hyXyX9Em+F/HNAqRMiT9YCZ//GApr83eSMsG1Xnlzr4M6JgU77BEjlJogjezasA9QAeCE6wj3T2rJAITKW3cnXcWbOF37yxh14HDDDu7HQ+e35EaVaKtJ4mUrKL2zDFAJjAV+A7QleAh2iLlyoadB/n1vxfxxYodnNKiLqNvP43uzevGOyyRQhWV4Lu6ew8AMxsJzIp9SCKlx5GsHF6cuob/+2QllSoYv7msK7ed0YaKesqSlAFFJfjM3DfunqUrA6Q8mbV2F798cyErt+/n4m5NeeSKrrphScqUohJ8TzPbG743gjtZ94bv3d31TDFJOGn7D/PH95fxr7mbaF6vOiNvS2ZQlybxDkvkmBWa4N1dzw2TciM7x5kwawNPfricg0eyGD6wPfed15HqVfQ1kLJJp/9FgPkb9/DrtxaxMDWd/u0a8vurutGhce14hy8fPNYAABFTSURBVCVyQpTgpVzbuf8wf56ynNfnbKRx7ao8e1NvLjulme5ElYQQswQfXkP/JVA1LGeSuz8Sq/JEjkVWdg6vztzA0x8t5+CRbIYOaMf9gzpSq6rqPJI4Yrk3HwbOc/f9ZlYZmGZmH7j71zEsU6RIX61O43fvLGHZ1n0M6NiIRy7vRofGteIdlkixi1mCd3cH9ocfK4cvPSRE4mbT7oM8/v5S3l+4lRb1q/P8LX24qJs6BpPEFdPj0bCr4TkEXR38w91n5jPNMILnvdKqVatYhiPl1MEjWTz3+WpGfLkGM/jRBZ0YenY7qlXW1TGS2GKa4N09G+hlZvWAN82su7svyjPNCGAEQHJysmr4Umxycpy35qfy5ynL2bo3gyt7ncRPLz6Zk+rpZiUpH0rkjJK77zGzz4CLgUVFTS9yomav28Vj7y7hm03pnNKiLv/4Xm9Obd0g3mGJlKhYXkWTBGSGyb06cAHwp1iVJwJBp2B/mrKM9xZuoWmdajxzfU+u6tWcCuo7RsqhWNbgmwFjwnb4CsDr7v5uDMuTciz9YCbPfrqSMTPWUalCBR4Y1JG7z2mnrnylXIvlVTQLgN6xWr4IwOGsbF6ZsZ5nP13F3oxMrju1BT+6sDNN9DxUEd3JKmVTTo7z7sItPPnhMjbuOsTZnZL42cUn0/Uk9X8nkksJXsqc6avSeOKDZSxMTadLszq8cmcPPTJPJB9K8FJmLEpN509TljF1ZRrN61XXCVSRIijBS6m3Lu0AT320nHcXbKFejcr86tIu3NKvtW5UEimCEryUWlvTM/jbJyt5PWUjVSpW4L7zOjD07HbUqVY53qGJlAlK8FLq7Nx/mOe/WM2YGetxd27t15rh57ancW1dGSNyLJTgpdRIP5jJi1PX8PL0tWRkZnNNnxY8MKgjLRvUiHdoImWSErzE3b6MTEZNX8eLU9ewLyOLy05pxoPnd1IXviInSAle4mb/4SzGfLWOEV+uIf1QJhd0bcJDF3SiSzNdyy5SHJTgpcTtP5zF2BnrePHLNew+mMmgkxvz4Pmd6NGibrxDE0koSvBSYvZlZDJ2xnpemhok9nM7J/HA+Z3o1bJevEMTSUhK8BJz6YcyGT19HS9PX0v6ISV2kZKiBC8xs3P/YV6evpaxX61n3+EsLujahPvP66imGJESogQvxW5regYvTl3D+JkbyMjK5pLuzRh+bnu6naTELlKSlOCl2KxLO8DzX6zmX3M3keNwZc+TGH5uezo0rh3v0ETKJSV4OWGLUtN57ovVfLBwC5UqVuDG01ox7Ox2ukFJJM6U4OW4uDvTV+3khS9XM3VlGrWrVuLuc9pz+5lt1KWASCmhBC/HJCs7h/cXbeWFL1azePNekmpX5eGLO3NLv9bqBEyklFGCl6jsP5zFxNkbeXnaWlL3HKJdUk3+9N0eXNW7OVUrqdtekdJICV4KtXnPIcZ8tY7xszawLyOLvm0b8NsrunHeyY31oA2RUi5mCd7MWgJjgSaAAyPc/W+xKk+K1/yNexg5bS3vL9wCwHe6N+WuAe10c5JIGRLLGnwW8CN3n2tmtYE5Zvaxuy+JYZlyAjKzc5iyaCujpq9l7oY91K5aidvPaMOQM9vQor6uiBEpa2KW4N19C7AlfL/PzJYCzQEl+FJm14EjTJi1gVdmrGfr3gxaN6zBo5d35drkltSqqlY8kbKqRL69ZtYG6A3MzGfcMGAYQKtWrUoiHAktSk1n9FfrePubzRzJyuGsDo147KrunHtyYyqqfV2kzIt5gjezWsC/gAfdfW/e8e4+AhgBkJyc7LGOp7w7nJXNBwu3MnbGOuZu2EONKhW5PrkFg/u3oVMT3XEqkkhimuDNrDJBcn/V3SfHsiwp3MZdB5kwawMTZ29k54EjtG1Uk99c1pXvntqCutV1/bpIIorlVTQGjASWuvszsSpHCpad43yxYjvjvt7AZ8u3Y8CgLk0Y3L81Z7ZvpMscRRJcLGvwZwK3AgvNbH447Bfu/n4MyxRg294MJs7eyMTZG0ndc4ik2lW579wO3NC3Fc3rVY93eCJSQmJ5Fc00QFXEEpJbW58wayOfLttOdo4zoGMjfnVpF87v2oTKFSvEO0QRKWG6Bq6M27jrIG+kbOSNOZvYkp5Bo1pVGTqgHTf1bUnrhjXjHZ6IxJESfBmUkZnNh4u38kbKJqatSsMMBnRM4pHLuzKoi2rrIhJQgi8j3J0Fm9KZNGcT/56fyt6MLJrXq84Pz+/Etckt1LYuIkdRgi/ltu/N4M15qUyas4mV2/dTtVIFLu7elOuTW9K/XUNdCSMiBVKCL4Vym2Amz01l6sod5Dj0aVWPx6/uwaWnNNN16yISFSX4UiI7x5m5ZieT56UyZdFW9h8OmmCGD+zA1X2a0z6pVrxDFJEyRgk+jtydxZv38vY3m3l7/ma27s2gdtVKXNKjKVf1bk6/tmqCEZHjpwQfB2vTDvD2/M28/U0qq3ccoFIFY2DnJH55aRcu6NqEapX1hCQROXFK8CVk855DvLdgC+8s2MyCTekA9G3bgNvPbMulPZpRv2aVOEcoIolGCT6GtqZn8MGiLby3YAsp63cDcEqLuvzyki5c1rMZzerq0kYRiR0l+GK2Jf0QHyzcygeLtjB7XZDUT25am59c1JlLezSjTSPdXSoiJUMJvhhs3HWQKYuCpD53wx4AOjepzUMXdOKSHs3o0FhXwIhIyVOCPw7uzsrt+/lw0VamLN7K4s3Bc0y6NqvDTy7qzMXdm+qyRhGJOyX4KOXkOPM27uajxdv4aMk21qYdAKB3q3r84pKTubhbM1o11IOpRaT0UIIvxKEj2UxflcZ/lm7jP0u3k7b/MJUqGP3bN+SOs9pyYdcmNKlTLd5hiojkSwk+j217M/hk6XY+XbaNaavSyMjMoXbVSpzdOYkLuzZhYOfG6ipARMqEcp/gs3Ocbzbt4bNl2/ls+XYWpQbt6S3qV+fG01pxfpcm9G3bgCqV1AWviJQt5TLB7z5whC9X7uDz5Tv4YsUOdh04QgWDU1vX5+GLOzPo5CZ0alKL4LGyIiJlU7lI8FnZOXyzKZ0vVwQJ/ZtNe3CHBjWrcE6nJAZ2TuKcTknUq6G7SUUkcSRsgt+46yDTVqUxdeUOpq1MY29GFmbQq2U9HhjUkXM6JXFKi3pUVGdeIpKgYpbgzexl4DJgu7t3j1U5udIPZjJjzU6mr0pj2qq0/17G2KROVS7q1pSzOyVxVodG6vNFRMqNWNbgRwN/B8bGsAwyMrO5YcTXLNy0hxyHGlUqcnrbBtzarzUDOjaiQ2O1pYtI+RSzBO/uX5pZm1gtP1e1yhVp16gmAzslcWaHRvRqWU9XvIiIUAra4M1sGDAMoFWrVse1jL/c0Ks4QxIRSQhxr+q6+wh3T3b35KSkpHiHIyKSMOKe4EVEJDaU4EVEElTMEryZTQBmAJ3NbJOZ3RmrskRE5GixvIrmplgtW0REiqYmGhGRBKUELyKSoJTgRUQSlLl7vGP4LzPbAayPdxxRaASkxTuIUkjb5WjaJkfTNjnaiWyT1u6e701EpSrBlxVmluLuyfGOo7TRdjmatsnRtE2OFqttoiYaEZEEpQQvIpKglOCPz4h4B1BKabscTdvkaNomR4vJNlEbvIhIglINXkQkQSnBi4gkKCX4KJjZD81ssZktMrMJZlbNzNqa2UwzW2VmE82sXD3stYBtMtrM1prZ/PBVrp7EYmYPhNtjsZk9GA5rYGYfm9nK8G/9eMdZkgrYJo+aWWrEfnJJvOOMNTN72cy2m9miiGH57hsW+L8wtywwsz7HW64SfBHMrDlwP5AcPjy8InAj8CfgL+7eAdgNlJveMgvZJgA/cfde4Wt+3IIsYWbWHRgK9AV6ApeZWQfgZ8An7t4R+CT8XC4Usk0g+O7k7ifvxy3IkjMauDjPsIL2je8AHcPXMOC54y1UCT46lYDqZlYJqAFsAc4DJoXjxwBXxSm2eMm7TTbHOZ546wLMdPeD7p4FfAFcA1xJsH9A+dtPCtom5Y67fwnsyjO4oH3jSmCsB74G6plZs+MpVwm+CO6eCjwFbCBI7OnAHGBPuNMCbAKaxyfCkpffNnH3j8LRfwgPK/9iZlXjFmTJWwQMMLOGZlYDuARoCTRx9y3hNFuBJvEKMA4K2iYAPwj3k5fLW7NVhIL2jebAxojpjju/KMEXIdz5rgTaAicBNTn6UKtcyW+bmNktwM+Bk4HTgAbAT+MWZAlz96UEzXYfAVOA+UB2nmkcKDfXJReyTZ4D2gO9CCoIT8crxtIiVvuGEnzRzgfWuvsOd88EJgNnEhw25T4wpQWQGq8A4yC/bXKGu28JDysPA6MI2l7LDXcf6e6nuvvZBOdlVgDbcg+vw7/b4xljSctvm7j7NnfPdvcc4EXK2X4SoaB9I5X/HenACeQXJfiibQD6mVkNMzNgELAE+Ay4NpzmNuDfcYovHvLbJksjdlYjaE9cVMgyEo6ZNQ7/tiJoax4PvE2wf0D520/y3SZ52pOvppztJxEK2jfeBgaHV9P0I2gC3ZLfAoqiO1mjYGa/BW4AsoB5wF0EbWKvETRFzANuCWuu5UIB2+QDIAkwgsPxe9x9f9yCLGFmNhVoCGQCD7n7J2bWEHgdaEXQFfb17p73ZFvCKmCbvELQPOPAOuDu401gZUX4jOqBBN0CbwMeAd4in30jrCD9naAp+CBwu7unHFe5SvAiIolJTTQiIglKCV5EJEEpwYuIJCgleBGRBKUELyKSoJTgpVQIb2fP7V1wa57eBktVT51mNtDMzojh8qub2RdmVtHM2uTpgXComc0xs/pm9pSZnRerOKTsq1T0JCKx5+47Ca6NxsweBfa7+1PxisfMKkX0NZTXQGA/8FUxLS+vO4DJ7p4dXBL932XcCtwHnOfuu83sWYI7QT+NNg4pX1SDl1LLzE4Na7JzzOzDiDtlPw87M0sxs6VmdpqZTQ771X4snKaNmS0zs1fDaSaFHV4Vtdy/mlkK8ICZXW5Bn//zzOw/ZtbEzNoA9wA/DI8uBljQD/61EXHvD/8ONLOpZvY2sCSskT9pZrPDjrbuLmDVv0eeO17N7HqC7mQvdPc0AHdfDzQ0s6bFtMklwSjBS2llwLPAte5+KvAy8IeI8UfcPRl4niAZ3gt0B4aEd48CdAb+6e5dgL3AcDOrXMRyq7h7srs/DUwD+rl7b4K7lh9293Vhmbn9mU8tYj36AA+4eyeCZwaku/tpBB2yDTWztt9a6aA5ql1YTq7WBHc2XujuW/Msfy5B30giR1ETjZRWVQkS9sdhM0VFgp4Hc70d/l0ILM691d3M1hB01LQH2Oju08PpxhE8pGRKEcudGPG+BTAxrOFXAdYex3rMcvfc+S4ETomo7dcleKhD5HIbhbFH2kHQl/j1wF/yjNtO0KOnyFGU4KW0MoLE3b+A8bn9/uREvM/9nLtf5+2Hw6NY7oGI988Cz7j722Y2EHi0gHmyCI+GzawCwY9Bfssz4D53/7CA5QAcAqrlGXaQoC/1qWa23d1fjRhXLZxH5ChqopHS6jCQZGb9Acysspl1O8ZltMqdH7iZoMll+TEsty7/66b1tojh+4DaEZ/XAaeG768AKhewvA+B74fNRJhZJzOrGTmBu+8GKppZtTzDtxN0PvW4mV0UMaoT5bc3RimCEryUVjkE3TH/ycy+Ieid8lgvTVwO3GtmS4H6wHPufuQYlvso8IaZzQHSIoa/A1yde5KV4EqWc8Ll9efbtfZILxF0NT03vPTxBfI/iv4IOCvvwLCp5wrgZTPrG/5QdACOq6dBSXzqTVISUni1y7vhQ8HLFDPrA/zQ3W8tYrqrgT7u/uuSiUzKGtXgRUoZd58LfGZmFYuYtBJ63J0UQjV4EZEEpRq8iEiCUoIXEUlQSvAiIglKCV5EJEEpwYuIJKj/B78aC7YO5tzbAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"cell_type": "code",
"source": [
"h_lg = lambda Tsat: 199000 + 0.68*C_pl*(Tsat-Twall)\n",
"h = lambda Tsat: 0.943*(N2_rho_liq*((N2_rho_liq-N2_rho_gas)*g*h_lg(Tsat)*k**3)/(L*N2_mu_liq*(Tsat-Twall)))**(0.25) # W/m2K\n",
"plt.plot(n2_P,h(T)*A*(T-Twall))\n",
"plt.title('LNTF TS Power vs Pressure')\n",
"plt.xlabel('Pressure (Bar)')\n",
"plt.ylabel('Max Power (Watts)')"
],
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/",
"height": 312
},
"id": "yaFtNrpioxNx",
"outputId": "254126ee-7969-4725-b6fe-6dbb8031614d"
},
"execution_count": 15,
"outputs": [
{
"output_type": "execute_result",
"data": {
"text/plain": [
"Text(0, 0.5, 'Max Power (Watts)')"
]
},
"metadata": {},
"execution_count": 15
},
{
"output_type": "display_data",
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEWCAYAAACJ0YulAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3dd3yV9fn/8dcVEvaGEEaYsmUTBRcqiKsqKq6quEv1q63a4fp1aFtrba1VO2ytqODCrbhFVFzsvQWZCSSMAGGHJNfvj/tOGhHIScjJSXLez8fjPHLu+9z3fa5zxPs6n23ujoiICEBCrAMQEZHKQ0lBRESKKCmIiEgRJQURESmipCAiIkWUFEREpIiSgoiIFFFSkEMys9VmdtpB9p9iZm5m/zpg/5dmdo2Z3WNmO8PHXjPLL7a9KDzWzWxXsf3bDvI+kVxnhJnNNbMcM9tsZp+YWcdDfJ5nzCw3PD/bzCaaWffy+bYqBzPrEH63hd/TajO7K9ZxSdWhpCBltQsYZWYdDnzB3f/o7vXdvT5wIzClcNvdjy52aN9i+xuX9jpm1hkYB/wcaAR0BP4J5B8m7j+H10sFNgLPlOGzHzEzS4zyWzQOP+cPgd+Y2ZkxiOGwLKB7UCWj/yBSVtsIbqi/jWEM/YBV7j7JAzvc/TV3X1vSie6+G3gB6AVgZj3M7DMz22Zmi8zsvHB/x3BfQrj9XzPbWHgdM3vWzG4LnzcyszFmtsHMMszsD2ZWI3ztGjP7ysz+ZmZbgHuLx2Nmrc1sj5k1Lbavf1j6STKzzmY22cy2h/teiuQLcvcpwCKgV1jCSzezO80sE3jazBLM7C4z+9bMtpjZy4UxmFltM3su3L/NzGaYWUqxz7PSzHaY2SozuyLcf6+ZPVfsMxSWXBLD7c/M7H4z+wrYDXQys+5hqS3bzJaZ2SWRfDaJDiUFORL3AyPNrFuM3n820D280Z5qZvUjPTE89gpgjpklAW8DHwEtgJ8Az5tZN3dfBeQA/cNThwA7zaxHuH0yMDl8/gyQB3QOjz8duKHY2w4CVgIpBN9dEXdfD0wBRhbbfTnwqrvvB34fxteEoJTz9wg+o5nZCcDRwJxwd0ugKdAeGB1+1vPDz9Ea2EpQ2gK4mqAE1hZoRlBa22Nm9YDHgLPcvQFwPDC3pHiKGRW+dwNgEzCRIEG3AC4D/mVmPUtxPSlHSgpSZu6eCfwb+F0ZLzE7/AW6zcweK8P7rwROAdoALwObw3aDwyWHX4TtFyuA+sA1wODw+Z/cPdfdPwHeIah6geCmf7KZtQy3Xw23OwINgXnhL+izgdvcfZe7bwT+RnCTK7Te3f/u7nnuvucgsb1Q+J5mZuG5L4Sv7Se4kbd2973u/mUJX89mIBt4ErjL3SeF+wuA37r7vjCGG4H/5+7p7r6PoARzUfjLfj9BMujs7vnuPsvdc4pdp5eZ1XH3De6+qIR4invG3Re5ex5wJrDa3Z8Ov5c5wGvAxaW4npSjmNYpSrXwIPCtmfUtw7kD3H3Fkby5u08FLgEws2OAl4D/B9x9iFMecvdfFd9hZmnAOncvKLZ7DUGygSApnAekA58DnxH82t0LfOHuBWbWHkgCNgT3cyD40bWu2DWLPz+Y14C/m1kroCvBjfeL8LU7CEoL081sK/BXd3/qMNdqHt50D7TJ3fcW224PvGFmxT97PkFp5lmCUsJ4M2sMPEeQQHaZ2aXAL4AxYVXQz919aQmfr1Dx76E9MMi+29EgMXxviQElBTki7r7FzB4huGHFOpYZZvY6YTtBKawH2ppZQrHE0A74Jnw+GfgLQVKYDHxJUELay/+qjtYB+zj0zRjgsFMSu/tWM/sIuBToAYz3cBrjsFT2IwAzOxH42Mw+L0NSPTCGdcB17v7VIY6/D7gv7FDwHrAMGOPuHwIfmlkd4A/Af4GTCDog1C12fku+r3gM64DJ7j68lJ9DokTVR1KSpLDBsfBxsB8SDxPUK/c4yGtRY2YnmtmPzKxFuN2d4Bf91FJeahpBo+cdYaPuKcC5wHgAd18O7AGuJLiB5QBZBPX/k8NjNhDU+f/VzBqGDbhHmdnJpYzlBeAq4CL+V3WEmV1sZqnh5laCG2vB908vtX8D94clHcws2cxGhM9PNbPeYWN5DkF1UoGZpVjQFbgeQSLcWSyWucAQM2tnZo04dImt0DtAVzMbFX73SWZ2TLE2G6lgSgpSkvcIboiFj3sPPCC8Sf6ZoAGzIm0jSAILzGwn8AHwRhhLxNw9lyAJnEVQF/8v4KoDqkMmA1vcfV2xbSNo7C50FVATWExw434VaFXKzzQB6AJkuvu8YvuPAaaFn3MCcGvYpnKkHg2v95GZ7SBIqIPC11oSfIYcYAnBZ36W4L7xM4ISVjZBI/VNAO4+kaAKbz4wi+Cmf0juvoOgQf6y8HqZBFWStcrhs0kZmBbZERGRQiopiIhIESUFEREpoqQgIiJFlBRERKRIlR6n0Lx5c+/QoUOswxARqVJmzZq12d2TD/ZalU4KHTp0YObMmbEOQ0SkSjGzNYd6LarVR2bW2MxeNbOlZrbEzI4zs6bhjIjLw79NwmPNzB4zsxVmNt/MBkQzNhER+b5otyk8Cnzg7t2BvgQDYO4CJrl7F2BSuA3BwKEu4WM08HiUYxMRkQNELSmEQ9yHAGMgGDXq7tuAEcDY8LCxBNP2Eu4fF86LPxVoHE4MJiIiFSSaJYWOBHOlP21mc8zsyXCulJRwnhgIhrSnhM/b8N3ZE9P53yyVRcxstJnNNLOZmzZtimL4IiLxJ5pJIREYADzu7v0JZk/8zlqx4QyQpZpnw92fcPc0d09LTj5o47mIiJRRNJNCOpDu7tPC7VcJkkRWYbVQ+LdwacMMgrnbC6WG+0REpIJELSmE87+vK7ZU4zCC2SMnECzzR/j3rfD5BOCqsBfSYGB7sWomERGpANEep1C41m1NgrVpryVIRC+b2fUEq1sVLtL9HsFyhisI5ra/NsqxiYhUKQUFzqL1OUxamsUZR7ekR6uG5f4eUU0K7j4XSDvIS8MOcqwDN0czHhGRqmbnvjy+XL6ZT5du5JNlG9m0Yx9m0Kx+raqXFEREpPTWZe9m0pIsJi3dyLSV2eTmF9CgdiJDuiYzrHsLTu6aTLP60VmHSElBRCTG8gucOWu38vGSjXyyNItvsnYC0Cm5Hlcf356h3VNI69CEpBrRn8NUSUFEJAZ27svj82828fGSLD5btonsXbkkJhjHdmzKpce0Y2j3FnRsXq/C41JSEBGpIOu37eHjJVlMXJxVVC3UqE4Sp3ZL5rSeKQzpmkzD2kkxjVFJQUQkStyD3kKFiWDR+hwAOjavxzUndGBY9xYMbN+ExAqoFoqUkoKISDnan1/A9FXZTFwcJIKMbXswg4HtmnDXWd0Z3jOFo5LrxzrMQ1JSEBE5Qrv25TH5m018tCiTT5ZuJGdvHrWTEjixczK3DuvC0B4taB6l3kLlTUlBRKQMtuzcx8dLsvhoURZfrNhMbl4BTeomcfrRLTm9ZwondUmmTs0asQ6z1JQUREQilL51Nx8uyuLDRZnMXJ1NgUObxnW4clB7Tj86hbRK1j5QFkoKIiKHsWLjDj5YmMkHizJZmBE0FHdv2YBbhnbhjKNT6NmqIWYW4yjLj5KCiEgxhT2GPliYyfsLN/Dtpl0A9G/XmLvP6s4ZR7ekQwzGD1QUJQURiXsFBc7c9G28v2ADHyzKZF32HhIMBndqxtXHd+D0ni1p2ah2rMOsEEoKIhKXCgqcWWu38t6CDXywMJMN2/eSVMM4oXNzbjm1M8N7tqRpvZqxDrPCKSmISNzIL3BmrM7m/QUbeH9hJht37KNmYgInd03mjjO7MbR7Co3qxHZEcawpKYhItVaYCN6dH1QNbdqxj1qJCZzarQVn9W7JsB4p1K+lW2EhfRMiUu0UVg29M2897y0MEkHtpASGdm/B2b1bcWq3FtRTIjgofSsiUi24O3PWbePteet5b8EGsnKCEsHQ7i34QZ9WDO3egro1dcsrib4hEamy3J2FGTm8M38978zfQMa2PdSskcDJ3ZI5p08rVQ2Vgb4tEalyVmzcwYS563l7/gZWbd5FYoJxYpfm/Gx4V4YfnRLz6aerMiUFEakS1mXv5u3565kwdz1LM3dgBsd1asboIZ048+iWNInD7qPRoKQgIpXW5p37eG/BBt6au55Za7YCwcji357bkx/0bkWLhvExoKwiKSmISKWya18eHy3O5M056/lyxWbyC5yuKfX55RndOK9va9o2rRvrEKs1JQURibn9+QV8uXwzb8zJYOLiLPbsz6dN4zqMHtKJEf1a071lw1iHGDeUFEQkJtydeenbeXNOBm/PW8+WXbk0qpPEBQPacH6/NqS1b0JCQvWZfbSqUFIQkQq1Lns3b8zJ4M05GazcvIuaiQkM75HCiH6tOaVbC2omVu31CKo6JQURibqcvft5f8EGXpudwfRV2QAM6tiUH5/ciTN7tYr7+YYqEyUFEYmKvPwCvlixmddnZ/DRokz25RXQKbkevzyjGyP6tSa1iRqMK6OoJgUzWw3sAPKBPHdPM7OmwEtAB2A1cIm7b7Vg6aJHgbOB3cA17j47mvGJSPlblrmD12an8+acDDbu2EfjuklcktaWkQNT6ZvaqFqtUlYdVURJ4VR331xs+y5gkrv/yczuCrfvBM4CuoSPQcDj4V8RqeS27c5lwrz1vDIznQUZ20lMME7p1oKLBrbh1O4tqJVY9Rawj1exqD4aAZwSPh8LfEaQFEYA49zdgalm1tjMWrn7hhjEKCIlyC9wvli+iVdmpjNxcRa5+QV0b9mAX5/TkxH9WtO8fq1YhyhlEO2k4MBHZubAf9z9CSCl2I0+E0gJn7cB1hU7Nz3c952kYGajgdEA7dq1i2LoInIwqzfv4tVZ6bw6K53MnL00qZvE5YPacdHAVHq1aRTr8OQIRTspnOjuGWbWAphoZkuLv+juHiaMiIWJ5QmAtLS0Up0rImWzJzef9xdu4KUZ65i2KpsEgyFdk/nNuT0Z1kPVQ9VJVJOCu2eEfzea2RvAsUBWYbWQmbUCNoaHZwBti52eGu4TkRhZmLGd8TPW8tbc9ezYm0e7pnX55RnduHBAG1o1qhPr8CQKopYUzKwekODuO8LnpwO/AyYAVwN/Cv++FZ4yAbjFzMYTNDBvV3uCSMXL2buft+auZ/z0tSxan0OtxATO6tWSS49px6COTTXKuJqLZkkhBXgj7H6WCLzg7h+Y2QzgZTO7HlgDXBIe/x5Bd9QVBF1Sr41ibCJSjLsze+02Xpy+lnfnb2DP/nx6tGrI70YczYh+bTS4LI5ELSm4+0qg70H2bwGGHWS/AzdHKx4R+b7te/bz5pwMXpy+lqWZO6hXswbn92/NZce0o4/GFMQljWgWiTPuztx123h+2lremb+evfsL6JPaiAcu7M25fVtr+co4p//6InFi57483pqbwfNT17J4Qw71atbgwgGpXH5sO3UllSJKCiLV3LLMHTw3dQ2vz05nV27QVnD/Bb0Y0a+NSgXyPfoXIVIN5eYV8MGiTJ6bsobpq7OpmZjAuX1ac8XgdvRv21htBXJISgoi1Ujm9r28MG0NL0xfx+ad+2jXtC73nN2diwe21cL2EhElBZEqzt2ZtiqbsV+v5qPFWRS4c2q3Fow6rj0nd0nWuAIpFSUFkSpqT24+b87NYOzXq1mauYNGdZK44cSOXDGoPe2aaa0CKRslBZEqJn3rbp6duobx09exfc9+erRqyIMje3Ne3zbUqak5iOTIKCmIVAHuzozVW3n6q1V8uCgTM+OMo1O45viOHNOhiRqOpdwoKYhUYrl5Bby7YD1jvlzFwowcGtVJYvSQoxh1XHvaNNaEdFL+lBREKqHsXbk8P3UN46auYdOOfXRuUZ/7L+jFhf1TVUUkUaWkIFKJrNi4k6e+WsVrs9LZl1fAkK7JPHRxR4Z0aa4qIqkQSgoiMVbYpfS/n69k0tKN1ExM4ML+bbjuxI50TWkQ6/AkzigpiMRIXn4B7y/M5L9frGR++naa1qvJrcO6MOq49lrfWGJGSUGkgu3OzeOVmen894uVpG/dQ8fm9bj/gl6MHJBK7SS1F0hsKSmIVJDsXbmM/Xo146asZuvu/Qxs34Rfn9OT4T1SNOpYKg0lBZEoS9+6mye/WMX4GWvZu7+A4T1T+PGQTqR1aBrr0ES+J+KkYGZNgNbAHmC1uxdELSqRamB51g4en/wtE+auB+D8/m248eROdG6hxmOpvA6bFMysEcESmT8EagKbgNpAiplNBf7l7p9GPUqRKmTeum3867MVfLgoizpJNbjquA7ccFJHWmuwmVQBJZUUXgXGASe5+7biL5jZQGCUmXVy9zHRClCkqpi2cgv/+HQFXyzfTMPaifx0WBeuOb4DTTVltVQhh00K7j78MK/NAmaVe0QiVYi78/nyzfzjk+XMWL2V5vVrcddZ3blycHutaiZVUkT/as3sBGCuu+8ysyuBAcCj7r4mqtGJVFLuzidLN/LYpOXMS99Oq0a1ue+8o7n0mLbqVipVWqQ/ZR4H+ppZX+DnwJME1UonRyswkcrI3flocRaPTVrOovU5tG1ahwcu7M3IAanUTEyIdXgiRyzSpJDn7m5mI4B/uPsYM7s+moGJVCaFyeDRj5ezeEMOHZrV5S8X9eH8/m1IqqFkINVHpElhh5ndDVwJDDGzBCApemGJVA7uzsTFWTwSJoOOzevx8CV9Oa9vaxKVDKQaijQpXApcDlzv7plm1g74S/TCEoktd+ezZZt4eOI3LMjYTodmdZUMJC5EmhRud/c7Czfcfa2ZHR2lmERixt35+tstPPTRMuas3UbbpnX4y0V9uKB/GyUDiQuRJoXhwJ0H7DvrIPu+x8xqADOBDHc/x8w6AuOBZgRdWke5e66Z1SJovB4IbAEudffVEcYncsRmrcnmLx8uY+rKbFo3qs0DF/bmooGpajOQuFLSiOabgP8DjjKz+cVeagB8HeF73AosARqG2w8Cf3P38Wb2b+B6gt5N1wNb3b2zmV0WHndpxJ9EpIwWrd/OQx8u49Nlm2hevxb3ntuTHw5qR61EdS2V+FNSSeEF4H3gAeCuYvt3uHt2SRc3s1TgB8D9wM8sWDpqKEH7BMBY4F6CpDAifA7BSOp/mJm5u0f0SURKadXmXTw88RvenreeRnWSuPPM7lx9fHvq1tSgM4lfJY1o3g5sN7O8Aweqmdmz7j6qhOs/AtxBULKAoMpom7vnhdvpQJvweRtgXfi+eWa2PTx+8wHvOxoYDdCuXbsS3l7k+zbm7OXRScsZP2MdNWskcMupnfnRkE40qqMOdSKR/iT6TqOymSUS1P0fkpmdA2x091lmdkrZwvs+d38CeAIgLS1NpQiJWM7e/TwxeSVjvlzF/vwCrhjUjluGdqZFg9qxDk2k0iipTeFu4B6gjpnlFO4GcglvzIdxAnCemZ1NMLNqQ+BRoLGZJYalhVQgIzw+A2gLpIdJpxFBg7PIEcnNK+D5aWt4bNJytu7ez7l9W/Pz4V3p0LxerEMTqXRKqj56AHjAzB5w97tLc+Hw+LsBwpLCL9z9CjN7BbiIoAfS1cBb4SkTwu0p4eufqD1BjoS7896CTB78YClrs3dz/FHNuPusHvRObRTr0EQqrYiqj9z97nCRnS4Ev/oL939ehve8ExhvZn8A5gCF026PAZ41sxVANnBZGa4tAgTdS+9/dwmz126je8sGPHPtMZzcNZmgr4OIHEqks6TeQNC1NBWYCwwm+EU/NJLz3f0z4LPw+Urg2IMcsxe4OJLriRzK2i27efCDpby7YAMtGtTizyP7MHJgKjW0BrJIRCJtaL4VOAaY6u6nmll34I/RC0ukdHL27uefn6zg6a9WUyPBuO20Lowe0kndS0VKKdL/Y/a6+14zw8xquftSM+sW1chEIpBf4Lw8cx0PfbiMLbtyGTkglV+e0Y2WjdSjSKQsIk0K6WbWGHgTmGhmWwEtsCMxNXXlFu57ezFLNuSQ1r4Jz1x7rBqRRY5QSV1S+wHz3P2CcNe9ZvYpQXfRD6IdnMjBZGzbwx/fW8K78zfQpnEd/v7D/pzTp5UakUXKQUklhSeBTmY2i2Cuo6+AKe6+I+qRiRxg7/58nvh8Jf/6bAXucNtpXfjxkKOoU1NzFImUl5LGKaSZWV2C3kLHAz8l6DaaCXzl7v9XATGKMGlJFve9vZi12bs5u3dL7jm7B6lN6sY6LJFqp8Q2BXffDXxmZjOAaQQjla8CzoxybCKs3bKb+95exKSlGzkquR7PXT+IE7s0j3VYItVWSW0KlxOUEPoB+4DCxHCiu2dGPzyJV/vy8vnP5JX889MVJCYY95zdnWuO70jNRK1tIBJNJZUU/gMsA/4NfO7u30Q/JIl3Xy7fzK/fWsiqzbv4QZ9W/PoHPdXFVKSClJQUGgN9CUoL94ZjEzYQjGae4u6fRDk+iSObduzj/ncX8+bc9XRoVpdx1x3LkK7JsQ5LJK6U1NCcD8wOH/8wsxSCqShuA34HqNuHHLGCAuelmet44L0l7N1fwE+HdeH/TjmK2kn65yVS0UpqU+hDUEoofNQk6Jr6d4LuqSJHZMXGHdz9+gJmrN7K4E5Nuf+C3hyVXD/WYYnErZKqj54huPm/D/zK3ddGPSKJC7l5BTz+2bf889MV1KlZgz9f1IeLB6ZqAJpIjJVUfTSgogKR+DFn7VbufG0+32Tt5Ny+rfntuT1pXr9WrMMSEUquPppAsMLah+6+/4DXOgHXAKvd/amoRSjVxp7cfP760TLGfLWKlg1rM+bqNIb1SIl1WCJSTEnVR6OBnwGPmlk2sIlgkZ0OwLfAP9z9rUOfLhKYunILd742nzVbdnPl4HbceWZ3GtROinVYInKAkqqPMoE7gDvMrAPQCtgDfBOOdBY5rN25eTz4/lLGTllDu6Z1efFHgznuqGaxDktEDiHiFUjcfTWwOmqRSLUzfVU2v3hlHmuzd3PN8R2448xuWvRGpJLT/6FS7vbuz+ehD4O2g7ZN6vLS6MEM6qTSgUhVoKQg5WpB+nZuf3kuKzbuZNTg9tx1Vnfq1dI/M5GqosT/W82sBjDO3a+ogHikisrLD8YdPDppOc3r19IUFSJVVCRTZ+ebWXszq+nuuRURlFQta7bs4vaX5jJ77TbO69ua34/oRaO66lkkUhVFWq5fCXwVjlvYVbjT3R+OSlRSJbg7r85K594Ji0hIMB69rB8j+rWJdVgicgQiTQrfho8EoEH0wpGqYvue/dzzxgLenb+BQR2b8vCl/WjTuE6swxKRIxRRUnD3+wDMrK7GJ8jM1dncOn4uWTl7+eUZ3bjx5KOokaA5i0Sqg4iWsTKz48xsMbA03O5rZv+KamRS6eQXOH+ftJxL/jOFGgnGqzcdz82ndlZCEKlGIq0+egQ4A5gA4O7zzGxI1KKSSmdjzl5ue2kuX3+7hfP6tub+C3ppmgqRaqg0I5rXHTCtcX75hyOV0VcrNnPr+Dns3JfHn0f24eI0TXEtUl1Fugr6OjM7HnAzSzKzXwBLDneCmdU2s+lmNs/MFplZYbtERzObZmYrzOwlM6sZ7q8Vbq8IX+9wBJ9LykF+gfO3id9w5ZhpNKlbkwm3nMglx7RVQhCpxiJNCjcCNwNtgAygX7h9OPuAoe7eNzz+TDMbDDwI/M3dOwNbgevD468Htob7/xYeJzGyZec+rn5qOo9OWs4F/dvw1i0n0DVFHc9EqrtIq492lnZEs7s7sDPcTAofDgwFLg/3jwXuBR4HRoTPAV4lWBPawutIBZq1Zis3Pz+brbtzeXBkby5JU+lAJF5EmhQWmlkW8EX4+NLdt5d0UjhFxiygM/BPgrEO29w9LzwknaD0Qfh3HYC755nZdqAZsPmAa44mWOeBdu3aRRi+RMLdGTdlDb9/ZzGtG9fhtZuOp1ebRrEOS0QqUETVR2GVzg+BBcAPgHlmNjeC8/LdvR+QChwLdD+CWAuv+YS7p7l7WnKy5tYpL3ty8/n5y/P47YRFnNw1mbd/cqISgkgciqikYGapwAnASUBfYBHwZaRv4u7bzOxT4DigsZklhqWFVII2CsK/bYF0M0sEGgFbIn0PKbv0rbv58bOzWLwhh9tP68pPhnYmQWMPROJSpNVHa4EZwB/d/cZITjCzZGB/mBDqAMMJGo8/BS4CxgNXA4XLeU4It6eEr3+i9oTom/LtFm5+YTb78wsYc3UaQ7trzWSReBZpUugPnAhcbmZ3AcuBye4+5jDntALGhu0KCcDL7v5OODJ6vJn9AZgDFF5jDPCsma0AsoHLSv9xJFLuznNT13Dv24vp2LweT4waSKfk+rEOS0RizCL9MW5m9QkSw0nAlQDu3j56oZUsLS3NZ86cGcsQqqTcvALufXsRL0xby7DuLXjksn4anSwSR8xslrunHey1SNsUZgK1gK8Jeh8Ncfc15ReiVJStu3K56flZTF2ZzU2nHMUvTu+muYtEpEik1UdnufumqEYiUbdi405uGDuD9dv28rdL+3JB/9RYhyQilUykSSHXzB4GCifBmwz8LpKxClI5fP3tZm58dhZJNRJ4cfQgBrZvGuuQRKQSinSai6eAHcAl4SMHeDpaQUn5emXmOq4aM52UhrV58+YTlBBE5JAiLSkc5e4ji23fF8ngNYktd+eRj5fz6KTlnNC5Gf+6YiCN6qhBWUQOLdKksMfMTnT3LwHM7ARgT/TCkiO1P7+Ae15fwCuz0rloYCoPXNibpBqRFgxFJF5FmhRuBMaZWeG8B1sJBppJJbRrXx43PT+bz7/ZxK3DunDbaV00oZ2IRKTEpGBm/QgmtLuMcEoKd8+JclxSRpt37uO6Z2awaH0OD47szaXHaNJAEYncYesTzOw3wMvASOBd4FIlhMprXfZuLv73FL7J2sETowYqIYhIqZVUUrgU6Ofuu82sGfAB8N/ohyWltSxzB1c9NY09ufk8f4O6nIpI2ZSUFPa5+24Ad99iZmqprITmrN3KNU/PoFZiAq/ceDzdWmqFNBEpm5KSQiczmxA+N+CoYtu4+3lRi0wi8vW3m/nR2Jk0q1+L528YRNumdWMdkohUYSUlhREHbD8UrUCk9D5duo7Y1GEAABGSSURBVJEfPzeL9k3r8twNg0hpWDvWIYlIFXfYpODukysqECmdDxZm8pMXZ9OtZQPGXTeIpvVqxjokEakGIh2nIJXIO/PXc+v4ufRJbcQz1x6rUcoiUm6UFKqYt+ZmcPtLcxnYvglPX3ss9WvpP6GIlJ+IehOZ2fcqq82sefmHI4dTmBDSOjTlGSUEEYmCSLuYzjCzwYUbZjaSYMEdqSAT5q3n9pfmcmzHpjxz7THUU0IQkSiI9M5yOfCUmX0GtAaaAUOjFZR813sLNhSVEJ665hjq1lRCEJHoiOju4u4LzOx+4FmCdRWGuHt6VCMTACYtyeKnL86hf9vGPK2EICJRFukazWOAo4A+QFfgHTP7u7v/M5rBxbsvl2/mpudmc3TrhjytKiMRqQCRtiksAE5191Xu/iEwCBgQvbBk1ppsfjRuJp2S6zH2umNpUFvdTkUk+iKtPnrkgO3twPVRiUhYvD6Ha56eQctGtXn2+kE0rquBaSJSMSKtPuoCPAD0BIq6p7p7pyjFFbfWbNnFVU9Np36tRJ67YRDJDWrFOiQRiSORVh89DTwO5AGnAuOA56IVVLzauGMvo8ZMJ7+ggGevH0SbxnViHZKIxJlIk0Idd58EmLuvcfd7gR9EL6z4s2Pvfq55agabduzjqWuOoXOL+rEOSUTiUKTdWfaFayksN7NbCJbl1F2rnOTmFXDTc7NZlrWDJ69Oo3+7JrEOSUTiVKQlhVuBusBPgYHAKODqw51gZm3N7FMzW2xmi8zs1nB/UzObaGbLw79Nwv1mZo+Z2Qozm29mcdG7yd256/X5fLliM3+6sDendmsR65BEJI5FlBTcfYa773T3dHe/1t0vdPepJZyWB/zc3XsCg4GbzawncBcwyd27AJPCbYCzgC7hYzRBG0a19+ik5bw+O4PbTuvCxWltYx2OiMS5w1YfFV9l7WAOt/Kau28ANoTPd5jZEqANwcI9p4SHjQU+A+4M949zdwemmlljM2sVXqdaen12Oo98vJyRA1K5dViXWIcjIlJim8JxwDrgRWAawZKcpWZmHYD+4TVSit3oM4GU8Hmb8L0KpYf7vpMUzGw0QUmCdu3alSWcSmHG6mzuem0Bgzs15YELe2NWpq9WRKRclVR91BK4B+gFPAoMBza7++RIV2Uzs/rAa8Bt7p5T/LWwVOClCdjdn3D3NHdPS05OLs2plca67N38+NlZtGlSh39fOZCaiZE27YiIRNdh70bunu/uH7j71QTtAiuAz8IeSCUysySChPC8u78e7s4ys1bh662AjeH+DKB4pXpquK9a2bkvjxvGziQvv4AxV6dptLKIVCol/kQ1s1pmdiHBYLWbgceANyI4z4AxwBJ3f7jYSxP4X8+lq4G3iu2/KuyFNBjYXt3aEwoKnNtfmsuKTTv55xUD6JSsXr0iUrmU1NA8jqDq6D3gPndfWIprn0DQdXWBmc0N990D/Al42cyuB9YAl4SvvQecTVAa2Q1cW4r3qhIe+2Q5Exdn8ZtzenJSl6pZ9SUi1ZsF1fqHeNGsANgVbhY/0AiaBBpGMbYSpaWl+cyZM2MZQsQ+XpzFDeNmMnJAKg9d3EcNyyISM2Y2y93TDvbaYUsK7q4W0HKwavMubn9pLr3bNOL+C3opIYhIpaWbfpTtzs3jpudmkVjD+PeogdROqhHrkEREDklLeUWRu/OrNxeyLGsHY689VrOeikilp5JCFL0yM53XZ2dw67AuDOmqhmURqfyUFKJkaWYOv35rISd2bs5PhmoKCxGpGpQUomB3bh63vDCHBrWT+Nul/aiRoIZlEaka1KYQBb97ezHfbtrJs9dpOU0RqVpUUihn7y/YwPgZ67jx5KM4sUvzWIcjIlIqSgrlKHP7Xu56fQF9Uhvxs+FdYx2OiEipKSmUk4IC5xevzCM3r4BHLu1HUg19tSJS9ejOVU6enbqGL1ds5lfn9NBEdyJSZSkplIOVm3bywPtLOKVbMpcfW3UX/hERUVI4QvkFzi9fnU+txBo8OFIT3YlI1aYuqUfo6a9WMWvNVv52aV9SGtaOdTgiIkdEJYUjsHrzLv7y4TJO69GC8/u1iXU4IiJHTEmhjNydu19fQM0aCfzh/N6qNhKRakFJoYxemZnOlJVbuPvsHrRspGojEakelBTKYPPOfdz/3hKO7dCUy45pG+twRETKjZJCGfzx3SXszs3jjxf2IkGT3YlINaKkUEpTV27h9TkZ/HjIUXRu0SDW4YiIlCslhVLYn1/Ab95aSGqTOtwytHOswxERKXdKCqUw9uvVfJO1k9+c01NrLYtItaSkEKHNO/fx6MfLOblrMsN7psQ6HBGRqFBSiNBfP/qGPfvz+fU5PTUmQUSqLSWFCCzZkMP4GWsZdVx7OrfQDKgiUn0pKUTgj+8toWHtJG4bpoVzRKR6U1IoweRvNvHF8s38dFgXGtVNinU4IiJRFbWkYGZPmdlGM1tYbF9TM5toZsvDv03C/WZmj5nZCjObb2YDohVXaRQUOH96fyntmtZl1OD2sQ5HRCTqollSeAY484B9dwGT3L0LMCncBjgL6BI+RgOPRzGuiL09fz1LNuTw89O7UjNRhSoRqf6idqdz98+B7AN2jwDGhs/HAucX2z/OA1OBxmbWKlqxRWJ/fgF//egberRqyLl9WscyFBGRClPRP39T3H1D+DwTKOzw3wZYV+y49HDf95jZaDObaWYzN23aFLVAX52Vztrs3dxxRjfNbyQicSNmdSLu7oCX4bwn3D3N3dOSk5OjEBnsy8vn75OW079dY07pFp33EBGpjCo6KWQVVguFfzeG+zOA4nNQp4b7YuLlGetYv30vPxveVQPVRCSuVHRSmABcHT6/Gnir2P6rwl5Ig4HtxaqZKlRuXgH/nrySge2bcGLn5rEIQUQkZqLZJfVFYArQzczSzex64E/AcDNbDpwWbgO8B6wEVgD/Bf4vWnGV5M05GWRs28MtQzurlCAicScxWhd29x8e4qVhBznWgZujFUuk8gucf322gl5tGnJKV7UliEj8Uef7Yj5clMnqLbu56WSVEkQkPikphNyd/0z+lvbN6nJmr5axDkdEJCaUFELTVmUzL307PzqpEzU0LkFE4pSSQuipL1fRpG4SFw1MjXUoIiIxo6QArN2ym4lLsrh8UDstsykicU1JARg3ZTU1zBg1uEOsQxERiam4Twp79+fzyqx0zujVkpaNasc6HBGRmIr7pPD2vPVs37OfKwdpvQQRkbhPCs9NW0vnFvUZ3KlprEMREYm5uE4KSzNzmLduGz88tp0Gq4mIEOdJ4ZWZ6STVMC7of9ClG0RE4k7cJoXcvALemJPBaT1SaFqvZqzDERGpFOI2KUz+ZhPZu3K5OE2D1URECsVtUnhrbgZN69XkpC6aDVVEpFBcJoWd+/L4eEkWZ/duSVKNuPwKREQOKi7viBMXZ7J3fwEj+qmBWUSkuLhMCvVrJTG8ZwoD2zWJdSgiIpVK1FZeq8yG90xheM+UWIchIlLpxGVJQUREDk5JQUREiigpiIhIESUFEREpoqQgIiJFlBRERKSIkoKIiBRRUhARkSLm7rGOoczMbBOwppSnNQc2RyGcaKpqMVe1eKHqxVzV4oWqF3NVixcij7m9ux90NtAqnRTKwsxmuntarOMojaoWc1WLF6pezFUtXqh6MVe1eKF8Ylb1kYiIFFFSEBGRIvGYFJ6IdQBlUNVirmrxQtWLuarFC1Uv5qoWL5RDzHHXpiAiIocWjyUFERE5BCUFEREpEjdJwcyeMrONZrYw1rFEwszamtmnZrbYzBaZ2a2xjqkkZlbbzKab2bww5vtiHVMkzKyGmc0xs3diHUskzGy1mS0ws7lmNjPW8ZTEzBqb2atmttTMlpjZcbGO6XDMrFv43RY+cszstljHdThmdnv4/9xCM3vRzGqX+Vrx0qZgZkOAncA4d+8V63hKYmatgFbuPtvMGgCzgPPdfXGMQzskMzOgnrvvNLMk4EvgVnefGuPQDsvMfgakAQ3d/ZxYx1MSM1sNpLl7lRhYZWZjgS/c/UkzqwnUdfdtsY4rEmZWA8gABrl7aQfKVggza0Pw/1pPd99jZi8D77n7M2W5XtyUFNz9cyA71nFEyt03uPvs8PkOYAnQJrZRHZ4HdoabSeGjUv/qMLNU4AfAk7GOpToys0bAEGAMgLvnVpWEEBoGfFtZE0IxiUAdM0sE6gLry3qhuEkKVZmZdQD6A9NiG0nJwqqYucBGYKK7V/aYHwHuAApiHUgpOPCRmc0ys9GxDqYEHYFNwNNhFd2TZlYv1kGVwmXAi7EO4nDcPQN4CFgLbAC2u/tHZb2ekkIlZ2b1gdeA29w9J9bxlMTd8929H5AKHGtmlbaqzszOATa6+6xYx1JKJ7r7AOAs4OawarSySgQGAI+7e39gF3BXbEOKTFjVdR7wSqxjORwzawKMIEjArYF6ZnZlWa+npFCJhfXyrwHPu/vrsY6nNMIqgk+BM2Mdy2GcAJwX1tGPB4aa2XOxDalk4S9D3H0j8AZwbGwjOqx0IL1YifFVgiRRFZwFzHb3rFgHUoLTgFXuvsnd9wOvA8eX9WJKCpVU2Gg7Blji7g/HOp5ImFmymTUOn9cBhgNLYxvVobn73e6e6u4dCKoJPnH3Mv/CqghmVi/seEBYDXM6UGl71Ll7JrDOzLqFu4YBlbazxAF+SCWvOgqtBQabWd3wvjGMoA2yTOImKZjZi8AUoJuZpZvZ9bGOqQQnAKMIfr0Wdo07O9ZBlaAV8KmZzQdmELQpVIlunlVICvClmc0DpgPvuvsHMY6pJD8Bng//XfQD/hjjeEoUJtzhBL+6K7WwFPYqMBtYQHBfL/N0F3HTJVVEREoWNyUFEREpmZKCiIgUUVIQEZEiSgoiIlJESUFERIooKUi1Ymb5YffdhWb2ipnVjXVMkTCzVoWztJrZKWa2Pfwc883sYzNrUcrrJZtZZe+qKpWQkoJUN3vcvV84E24ucGPxF8MJwypEKd/rZ8B/i21/EX6OPgRjPm4uzfu6+yZgg5mdUIoYRJQUpFr7Augc/vL+wswmAIvDSfv+YmYzwl/iP4aiX+ufFytpnBQe+0y4vcDMbg+P/czM0sLnzcOpMjCza8xsgpl9AkwKRyA/ZcE6E3PMbMQhYh0JfO+XfThCtQGwNdw+1symhNf6unCk8IHvG57+JnBFuXyTEjcq7FeTSEUKf6Wfxf9utAOAXu6+KpxZdLu7H2NmtYCvzOwj4ELgQ3e/P5xHvy7BCNw2hWtwFE7jUYIBQB93zzazPxJMn3FdeO50M/vY3XcVi7UjsNXd9xW7xknhbLPNCCaRuyfcvxQ4yd3zzOw0gtHBIw9833B7JvCHyL4xkYCSglQ3dcKbKQQlhTEEk4NNd/dV4f7TgT5mdlG43QjoQlBN81Q4EeGb7j7XzFYCnczs78C7QCRTEk8sdmM+nWDSvV+E27WBdnx3bppWBNNLF/dF4YI/ZnYn8GeCqrBGwFgz60IwhXbSId4XgunLW0cQr0gRJQWpbvaEU3cXCWpg2FV8F/ATd//wwJPDaah/ADxjZg+7+zgz6wucQXBTvgS4Dsjjf9WvBy59eOB7jXT3ZYeL+SDXKG4CwWy5AL8HPnX3C8J1Nj47xPsWxrXnMNcV+R61KUg8+hC4KSwRYGZdw7r/9kCWu/+XYCW2AWbWHEhw99eAX/G/aZ9XAwPD5xdxaB8CPwnbBjCz/gc55hugw2GucSLwbfi8EcHykADXHOYcgK5U4hlUpXJSSUHi0ZMEN+HZ4c16E3A+cArwSzPbT7Ce91UES6A+bWaFP6DuDv8+BLwctk+8e5j3+j3B6m7zw2usAr6zDrS77zKzb82ss7uvCHcXtikYsB24Idz/Z4Lqo1+V8L4Ap0ZwjMh3aJZUkUrAzC4ABrr7r8rxmp8DI9x9a3ldU6o/lRREKgF3f8PMmpXX9cwsGXhYCUFKSyUFEREpooZmEREpoqQgIiJFlBRERKSIkoKIiBRRUhARkSL/H07Rw2Q64mwIAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"cell_type": "code",
"source": [
""
],
"metadata": {
"id": "2R9nv2Qbsmy1"
},
"execution_count": 9,
"outputs": []
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment