Newer
Older
{
"cells": [
{
"cell_type": "markdown",
"source": [
"Well use a GenericCascadeCase and we will load the example dataset, generated with OpenFOAM.\n",
"We don't necessarily need a complete dataset for postprocessing procedures."
],
"metadata": {
"collapsed": false,
"pycharm": {
"name": "#%% md\n"
}
}
},
{
"cell_type": "code",
"from ntrfc.cascade_case.solution import GenericCascadeCase\n",
"inlet_path = importlib.resources.files(\"ntrfc\") / \"data/openfoam_cascade_case/boundary/inlet.vtp\"\n",
"outlet_path = importlib.resources.files(\"ntrfc\") / \"data/openfoam_cascade_case/boundary/outlet.vtp\"\n",
"blade_path = importlib.resources.files(\"ntrfc\") / \"data/openfoam_cascade_case/boundary/blade_wall.vtp\"\n",
"fluid_path = importlib.resources.files(\"ntrfc\") / \"data/openfoam_cascade_case/internal.vtu\"\n",
"case.read_meshes(inlet_path,\"inlet\")\n",
"case.read_meshes(outlet_path,\"outlet\")\n",
"case.read_meshes(blade_path,\"blade\")\n",
"case.read_meshes(fluid_path,\"fluid\")\n"
],
"metadata": {
"collapsed": false,
"pycharm": {
"name": "#%%\n"
}
}
},
{
"cell_type": "code",
{
"name": "stdout",
"output_type": "stream",
"text": [
"ratio of cleaned sites: 0.8360508152445734\n"
]
"source": [
"\n",
"case.set_active_blade_slice()\n",
"case.compute_domainparams_from(0.01)"
],
"metadata": {
"collapsed": false,
"pycharm": {
"name": "#%%\n"
}
}
},
{
"cell_type": "code",
"text/plain": "<Figure size 640x480 with 1 Axes>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAHJCAYAAACmFmJFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABam0lEQVR4nO3deVyU1f4H8M+wDyCgLK4oqERiqLlmXnMlEVu8lhZimKn5K0nTNvNmaqaZuaV5rbxdTVOMW1lWuKBmairirmGICu6KIILKKjy/P6YZGZiBWZ5n5pmZz/v18pXzcOaZMydkvpzzPd+jEARBABEREZGMOVm7A0RERER1YcBCREREsseAhYiIiGSPAQsRERHJHgMWIiIikj0GLERERCR7DFiIiIhI9hiwEBERkewxYCEiIiLZY8BCZCdmzJgBhUKB3NzcOtuGhITgxRdfFO21d+7cCYVCgZ07d4pyv+zsbCgUCqxatUqU+xlj1apVUCgUyM7O1lzr3bs3evfubfG+ENF9DFiIiIhI9lys3QEiIrnbunWrtbtA5PAYsBAR1cHNzc3aXSByeFwSIrIzubm5GDZsGHx8fODv74+JEyeipKSk1ufcvHkTb775JiIjI+Ht7Q0fHx8MHDgQx44dq9H20qVLGDx4MLy8vBAUFIRJkyahtLRU531TU1MRHR0NX19feHp6olevXvjjjz9Mfm87duxAz5494eXlBT8/Pzz99NM4deqUVpvz58/j1VdfRXh4OJRKJfz9/TF06FCtnBS1P//8E3379oVSqUSzZs3w4YcforKyska76jks6pydpKQkzJ49G82aNYOHhwf69euHM2fO1Hj+smXL0LJlSyiVSnTt2hW7d+9mXgyRkTjDQmRnhg0bhpCQEHz00UfYv38/lixZgvz8fKxevVrvc86dO4cff/wRQ4cORWhoKK5fv44vvvgCvXr1Qnp6Opo0aQIAKC4uRr9+/XDhwgVMmDABTZo0wZo1a7Bjx44a99yxYwcGDhyITp06Yfr06XBycsLKlSvRt29f7N69G127djXqfW3btg0DBw5Ey5YtMWPGDBQXF2Pp0qXo0aMHDh8+jJCQEABAWloa9u7di+effx7NmjVDdnY2li9fjt69eyM9PR2enp4AgGvXrqFPnz64d+8epkyZAi8vL3z55ZdQKpUG92nu3LlwcnLCm2++iYKCAsybNw9xcXFITU3VtFm+fDkSEhLQs2dPTJo0CdnZ2Rg8eDDq16+PZs2aGTUGRA5NICK7MH36dAGA8NRTT2ldf/XVVwUAwrFjxzTXWrRoIYwcOVLzuKSkRKioqNB6XlZWluDu7i588MEHmmuLFy8WAAhJSUmaa3fv3hVat24tABB+++03QRAEobKyUggLCxMGDBggVFZWatoWFRUJoaGhQlRUVK3vJSsrSwAgrFy5UnOtQ4cOQlBQkJCXl6e5duzYMcHJyUmIj4/Xeo3q9u3bJwAQVq9erbn2+uuvCwCE1NRUzbWcnBzB19dXACBkZWVprvfq1Uvo1auX5vFvv/0mABDatGkjlJaWaq5/+umnAgDhxIkTgiAIQmlpqeDv7y906dJFKC8v17RbtWqVAEDrnkRUOy4JEdmZ8ePHaz1+7bXXAADJycl6n+Pu7g4nJ9WPg4qKCuTl5cHb2xvh4eE4fPiwpl1ycjIaN26MZ599VnPN09MTL7/8stb9jh49iszMTAwfPhx5eXnIzc1Fbm4u7t69i379+mHXrl06l170uXr1Ko4ePYoXX3wRDRo00Fxv164doqKitN5b1RmS8vJy5OXloXXr1vDz86vxXh555BGtmZ7AwEDExcUZ3K9Ro0Zp5bf07NkTgGrGCgAOHjyIvLw8jB07Fi4u9ye04+LiUL9+fYNfh4i4JERkd8LCwrQet2rVCk5OTjpzONQqKyvx6aef4t///jeysrJQUVGh+Zq/v7/m7+fPn0fr1q2hUCi0nh8eHq71ODMzEwAwcuRIva9ZUFBg8If2+fPndb4OALRp0wZbtmzB3bt34eXlheLiYnz00UdYuXIlLl++DEEQtF6z6j27detW4366XkOf5s2baz1Wv5/8/Hytfrdu3VqrnYuLi2YJi4gMw4CFyM5VDy50mTNnDqZNm4aXXnoJs2bNQoMGDeDk5ITXX3/dqJkQNfVzPvnkE3To0EFnG29vb6Pva4jXXnsNK1euxOuvv47u3bvD19cXCoUCzz//vEnvpTbOzs46r1cNkohIHAxYiOxMZmYmQkNDNY/PnDmDysrKWn+j/+6779CnTx989dVXWtdv3bqFgIAAzeMWLVrg5MmTEARBKxDKyMjQel6rVq0AAD4+Pujfv785b0fzurpeBwD++usvBAQEwMvLS/NeRo4ciQULFmjalJSU4NatWzXuqZ4JqkrXa5jb7zNnzqBPnz6a6/fu3UN2djbatWsn2msR2TvmsBDZmWXLlmk9Xrp0KQBg4MCBep/j7OxcY1bgf//7Hy5fvqx1LSYmBleuXMF3332nuVZUVIQvv/xSq12nTp3QqlUrzJ8/H3fu3Knxejdu3DDszfytcePG6NChA77++mutwOPkyZPYunUrYmJian0vS5cu1VrmUr+X/fv348CBA1r9Wrt2rVF9q03nzp3h7++PFStW4N69e5rra9eu1SwbEZFhOMNCZGeysrLw1FNPITo6Gvv27cM333yD4cOHo3379nqf88QTT+CDDz7AqFGj8Oijj+LEiRNYu3YtWrZsqdVu7Nix+OyzzxAfH49Dhw6hcePGWLNmjWarsJqTkxP+85//YODAgWjbti1GjRqFpk2b4vLly/jtt9/g4+ODn3/+2aj39cknn2DgwIHo3r07Ro8erdnW7OvrixkzZmi9lzVr1sDX1xcRERHYt28ftm3bppWLAwBvv/021qxZg+joaEycOFGzrblFixY4fvy4UX3Tx83NDTNmzMBrr72Gvn37YtiwYcjOzsaqVavQqlUrg5briOhvVt2jRESiUW9rTk9PF5599lmhXr16Qv369YWEhAShuLhYq62ubc1vvPGG0LhxY0GpVAo9evQQ9u3bV2M7ryAIwvnz54WnnnpK8PT0FAICAoSJEycKmzdv1trWrHbkyBFhyJAhgr+/v+Du7i60aNFCGDZsmLB9+/Za34uubc2CIAjbtm0TevToISiVSsHHx0d48sknhfT0dK02+fn5wqhRo4SAgADB29tbGDBggPDXX3/VeM+CIAjHjx8XevXqJXh4eAhNmzYVZs2aJXz11VcGb2v+3//+Z1C/lyxZIrRo0UJwd3cXunbtKvzxxx9Cp06dhOjo6FrHgYjuUwgCs8OIiCypsrISgYGBGDJkCFasWGHt7hDZBOawEBFJqKSkpEZOzerVq3Hz5k2W5icyAmdYiIgktHPnTkyaNAlDhw6Fv78/Dh8+jK+++gpt2rTBoUOHeLAikYGYdEtEJKGQkBAEBwdjyZIluHnzJho0aID4+HjMnTuXwQqRETjDQkRERLLHHBYiIiKSPQYsREREJHt2kcNSWVmJK1euoF69eizEREREZCMEQcDt27fRpEkTzYnx+thFwHLlyhUEBwdbuxtERERkgosXL6JZs2a1trGLgKVevXoAVG/Yx8dH9PuXl5dj69atePzxx+Hq6ir6/UmF42wZHGfpcYwtg+NsGVKOc2FhIYKDgzWf47Wxi4BFvQzk4+MjWcDi6ekJHx8f/qOQEMfZMjjO0uMYWwbH2TIsMc6GpHMw6ZaIiIhkjwELERERyZ4sApaKigpMmzYNoaGhUCqVaNWqFWbNmlXj/A0iIiJyTLLIYfn444+xfPlyfP3112jbti0OHjyIUaNGwdfXFxMmTLB294iIqApBEHDv3j1UVFRYtR/l5eVwcXFBSUmJ1ftiz8wdZ2dnZ7i4uJhddkQWAcvevXvx9NNPY9CgQQBUZ28kJibiwIEDVu4ZERFVVVZWhqtXr6KoqMjaXYEgCGjUqBEuXrzIGlwSEmOcPT090bhxY7POz5JFwPLoo4/iyy+/xOnTp/HAAw/g2LFj2LNnDxYuXKizfWlpKUpLSzWPCwsLAaiiwPLyctH7p76nFPem+zjOlsFxlp69jnFlZSWysrLg7OyMxo0bw9XV1aqBgiAIuHv3Lry8vBiwSMiccRYEAeXl5bhx4wbOnTuH0NBQrQJxxvwbkcXhh5WVlZg6dSrmzZsHZ2dnVFRUYPbs2Xj33Xd1tp8xYwZmzpxZ4/q6devg6ekpdXeJiBySi4sLGjVqhGbNmsHd3d3a3SEbUlpaikuXLuHq1atay0pFRUUYPnw4CgoK6ixLIouAZf369XjrrbfwySefoG3btjh69Chef/11LFy4ECNHjqzRXtcMS3BwMHJzcyWrw5KSkoKoqCju9ZcQx9kyOM7Ss9cxLikpwcWLFxESEgIPDw9rd0dT1p3HskhLjHEuKSlBdnY2goODtb53CgsLERAQYFDAIoslobfeegtTpkzB888/DwCIjIzE+fPn8dFHH+kMWNzd3XVG966urpL+cJD6/qTCcbYMjrP07G2MKyoqoFAo4OTkVOe5L5ZQWVkJAJo+kTTEGGcnJycoFIoa/yaM+fchi4ClqKioxiA4OztrBomILKisGEh5D7h5DvBtBsAJyM8C7uQAXkFAg1AAlUDBJcAnGFAo9H/dt7nq7/nn7389MAyI+hBwU1r3fRKRTZFFwPLkk09i9uzZaN68Odq2bYsjR45g4cKFeOmll6zdNSLHoA5S/twAFOXpb3fjFJD9u/lfT/sPoPQHvIMYxJBd6t27Nzp06IDFixdb5fWzs7MRGhqKI0eOoEOHDjrb7Ny5E3369EF+fj78/Pws2j9TyCJgWbp0KaZNm4ZXX30VOTk5aNKkCcaNG4f333/f2l0jsi9lxcDW6QB6AMnvAIpy4K9faw9SpFKcp/pTNYjxCQZa91HNzjRoySCGZE/fh/4PP/xg1eXA4OBgXL16FQEBAVbrg9hkEbDUq1cPixcvtlokSmTXqs+eOHkA7XsAx9YClSXW7p22wovA4dWqv5/doQpifJsD4w8wcLEjFZUCDmTdRM7tEgTV80DX0AZwdrKvpNkGDRpY9fWdnZ3RqFEjq/ZBbMxSIrIXlRVA1m7gxHfA6W3AL5OAeS2BOY1UH/zWmEURQ8EF1XtYFKkKvsimbT55Ff/4eAdiV+zHxPVHEbtiP/7x8Q5sPnlVstf87rvvEBkZCaVSCX9/f/Tv3x93794FoFq6ef3117XaDx48GC+++KLmcWlpKd555x0EBwfD3d0drVu3xldffYXs7Gz06dMHAFC/fn0oFArN86rfNz8/H/Hx8ahfvz48PT0xcOBAZGZmar6+atUq+Pn5YcuWLWjTpg28vb0RHR2Nq1f1j0t+fj7i4uIQGBgIpVKJsLAwrFy5EoBqSUihUODo0aOa9snJyXjggQegVCrRp08fZGdn17jnnj170LNnTyiVSgQHB2PChAmasbI2WcywEJGZ0jcCm94Bbl+xdk+kow5cmPtiszafvIpXvjmM6rU0rhWU4JVvDmP5iI6IfqixqK959epVxMbGYt68efjnP/+J27dvY/fu3UadVRcfH499+/ZhyZIlaN++PbKyspCbm4vg4GB8//33eOaZZ5CRkQEfHx8olbq/F1988UVkZmZi48aN8PHxwTvvvIOYmBikp6drlo6Kioowf/58rFmzBk5OThgxYgTefPNNrF27Vuc9p02bhvT0dGzatAkBAQE4c+YMiot1B/UXL17EkCFDMH78eLz88ss4ePAg3njjDa02Z8+eRXR0ND788EP897//xY0bN5CQkIDXXntNFisgDFiIbFlZMfDtCODsNmv3xHJ05b6ExwCxidbuGdWiolLAzJ/TawQrACAAUACY+XM6oiIaibo8dPXqVdy7dw9DhgxBixYtAKhKZxjq9OnTSEpKQkpKCvr37w8AaNmypebr6qWfoKAgvYmr6kDljz/+wKOPPgoAWLt2LYKDg/Hjjz9i6NChAFT1ez7//HO0atUKAJCQkIAPPvhAb98uXLiAhx9+GJ07dwagOtZGn+XLl6NVq1ZYsGABACA8PBwnTpzAxx9/rGnz0UcfIS4uTjMzFBYWhiVLlqBXr16YO3euJHXOjMElISJbU1YM/PrG/eUeRwpW9MlIBtYOs3YvqBYHsm7iaoH+nCkBwNWCEhzIuinq67Zv3x79+vVDZGQkhg4dihUrViA/P9/g5x89ehTOzs7o1auXyX04deoUXFxc0K1bN801f39/hIeH49SpU5prnp6emmAFABo3boycnBy9933llVewfv16dOjQAW+//Tb27t1bax+qvj4AdO/eXevxsWPHsGrVKnh7e2v+DBgwAJWVlTh//rzB71cqDFiIbIEmSGll+zkpUsncAsyPYJ6LTOXcNizB29B2hnJ2dkZKSgo2bdqEiIgILF26FOHh4cjKygKgKmhWfXmo6vk2+pZ4pFB9V5FCoah16WrgwIE4f/48Jk2ahCtXrqBfv3548803TX79O3fuYNy4cTh69Kjmz7Fjx5CRkYHQ0FCT7ysWBixEcqQOUNb8U5VsqglScq3dM3m7c1k1VvNaqcaPwYtsBNUzrJS/oe2MoVAo0KNHD8ycORNHjhyBm5sbNmzYAAAIDAzUSmytqKjAyZMnNY8jIyNRWVmJ33/XXV9Iffpw1fNxqmvTpg3u3buH1NRUzbW8vDxkZGQgIiLCrPcWGBiIkSNH4ptvvsHixYvx5Zdf6u3DgQMHtK7t379f63HHjh2Rnp6O1q1b1/hjzinLYmEOC5FcaLYf/yiPwMQzAHgwBqJWui28BpQaPh1vsqJcVYDHbdGy0TW0ARr7euBaQYnOPBYFgEa+qi3OYkpNTcX27dvx+OOPIygoCKmpqbhx4wbatGkDAOjbty8mT56MX3/9Fa1atcLChQtx69YtzfNDQkIwcuRIvPTSS5qk2/PnzyMnJwfDhg1DixYtoFAo8MsvvyAmJgZKpRLe3t5afQgLC8PTTz+NsWPH4osvvkC9evUwZcoUNG3aFE8//bTJ7+39999Hp06d0LZtW5SWluKXX37RvK/q/u///g8LFizAW2+9hTFjxuDQoUNYtWqVVpt33nkHjzzyCBISEjBmzBh4eXkhPT0dW7duxezZs03up1gYsBDJQWKsKg/D2jwDgLaDpd15ow7McjOBSweBcom3TKp3FzFwsSpnJwWmPxmBV745DAWgFbSoU2ynPxkhej0WHx8f7Nq1C4sXL0ZhYSFatGiBBQsWYODAgQCAl156CceOHUN8fDxcXFwwadIkzVZlteXLl2Pq1Kl49dVXkZeXh+bNm2Pq1KkAgKZNm2LmzJmYMmUKRo0ahfj4+BqBAACsXLkSEydOxBNPPIGysjI89thjSE5ONqu4nJubG959911kZ2dDqVSiZ8+eWL9+vc62zZs3x/fff49JkyZh6dKl6Nq1K+bMmaNVUb5du3b4/fff8a9//Qs9e/aEIAho1aoVhg2TR36YLE5rNldhYSF8fX0NOu3RFOXl5UhOTkZMTIxdHWQmNw45zmXFwLKuqg9VCyl38kBy+y8Rg9/hqii3flVZrbOLmgOnU1RLO1KxQOBir9/LJSUlyMrKQmhoqMmnNW8+eRUzf07XSsBt7OuB6U9GGL2lubKyEoWFhfDx8eHhhxISY5z1fe8Y8/nNGRYiS6r64Zx7xqKBimb2pPcMYNsOIOZjQA4fpm5KYNAC7WtSBnLqGRduhbaK6IcaIyqikd1XuiXxMWAhkpo1c1N8g4EHBmjPnlTZASFbbkpg0glpx069FTouSdz7Up2cnRTo3srf2t0gG8OAhUhK1shN8WsBdBsHdBkLuFg/s98s6tmXQQvuBy8nfwCKRarVod4KPeEQc1uIZI4BC5HYqh82aAmWSJa1tqrBi5iBoHorNJNyiWSNAQuRmCw5o+LpD7T9p30HKfrEJoqf58LcFiJZY8BCJAZL7PbxbQ4EtLbujh45qZHnItKMVkYy8M0wYARzW4jkhAELkSkstdvHkWdRDFU9z0WMwPHMFmDhQ0BCGsedSCYYsBAZyxLLPsynME3VWRdzA5fCi1wiIpIRVtoxQklZBab9eAIvfJWKaT+eQHGZ/rMjyA6VFavO9ZE6WAmPUX3oMlgxnTpwmXpNdXyAOTKSVUEqEVkVAxYjdJ6zDWv2X8DuzFys2X8Bbd7fjLGr06zdLbKExFjVb9tSLP34Ngda9QW6jFF9wPK3efG4KYHJJ4HWA8y7T0YyUHJHnD4RyciqVavg5+dXa5uZM2eiZ8+elulQLRiwGGBC4mG9X0tJz2HQYq/UJybPayn+rIqn//0AZdIJ4IUNqhwMzqpIY0SSaubKHB83B9I3itMfIpl47rnncPr0aWt3wyDMYalDcVkFdmTcQHRX/W1S0nNQXFYBpZuz5TpG4pM6kdYRaqXImblboYUKIOkFYMhXQLtnxe+fI6msAM7vBe5cB7wbAi0eBZzk9fOzrKwMbm7yKrwoRZ+USiWUStv4ecQZljrMSU4XtR3JlHrJJ+0/wNkd4gcr4THA22c5i2JtVXNbvJuado8fRgPrnhe3X44kfSOw+CHg6yeA70er/rv4IUlnr3r37o2EhAQkJCTA19cXAQEBmDZtGqqe/RsSEoJZs2YhPj4ePj4+ePnllwEAe/bsQc+ePaFUKhEcHIwJEybg7t37J4z/+9//RlhYGDw8PNCwYUM8++z9YPa7775DZGQklEol/P390b9/f81ze/fujddff12rn4MHD8aLL75odp+qO3bsGPr06YN69erBx8cHnTp1wsGDBwHoXhKaO3cuGjZsiHr16mH06NEoKSmpcc///Oc/aNOmDTw8PPDggw/i3//+dy3/B8TBgKUO2XlForYjGZFyyUfNtznzUuTITQm8mQ6EmZjbcnqTqlYLGSd9I5AUDxRe0b5eeFV1XcKg5euvv4aLiwsOHDiATz/9FAsXLsR//vMfrTbz589H+/btceTIEUybNg1nz55FdHQ0nnnmGRw/fhzffvst9uzZg4SEBADAwYMHMWHCBHzwwQfIyMjA5s2b8dhjjwEArl69itjYWLz00ks4deoUdu7ciSFDhmgFSYYwtk+6xMXFoVmzZkhLS8OhQ4cwZcoUvaeIJyUlYcaMGZgzZw4OHjyIxo0bY/ny5Vpt1q5di/fffx+zZ8/GqVOnMGfOHEybNg1ff/21Ue/NWFwSqkOIvycOnDOsHdkQS2xN5nZY+YtLMv17gbVajFNZAWx+B4CuD2wBgALYPAV4cJAky0PBwcFYtGgRFAoFwsPDceLECSxatAhjx47VtOnbty/eeOMNzeMxY8YgLi5OMxMSFhaGJUuWoFevXli+fDkuXLgALy8vPPHEE6hXrx5atGiBhx9+GIAqYLl37x6GDBmCFi1aAAAiIyON7rexffLw8KhxjwsXLuCtt97Cgw8+qHmOPosXL8bo0aMxevRoAMCHH36Ibdu2ac3gTJ8+HQsWLMCQIUMAAKGhoUhPT8cXX3yBkSNHGv0eDcUZljpczi82qN3UmAiJe0KikHJrMnf72KbYRODNLNOeq67Vwm3PdTu/t+bMihYBKLysaieBRx55BAqFQvO4e/fuyMzMREXF/fIUnTt31nrOsWPHsGrVKnh7e2v+DBgwAJWVlcjKykJUVBRatGiBli1b4oUXXsDatWtRVKSabW/fvj369euHyMhIDB06FCtWrEB+fr7R/Ta2T7pMnjwZY8aMQf/+/TF37lycPXtW7+udOnUK3bp107r2yCOPaP5+9+5dnD17FqNHj9bqw4cffljrfcXAGZZaqBNu3esI9vuGBzDh1hZINavCIm+2z7uBakbM1O8Pda0WBqn63bkubjsJeHl5aT2+c+cOxo0bhwkTJtRo27x5c7i5ueHw4cPYuXMntm7divfffx8zZsxAWloa/Pz8kJKSgr1792Lr1q1YunQp/vWvfyE1NRWhoaFwcnKqsTxUXl5udp90mTFjBoYPH45ff/0VmzZtwvTp07F+/Xr885//rHNMqrtzR7W9f8WKFTUCG2dnaT8HOcNSC0MTaZvW53KQbGnyVFpJE6ywyJv9iE0EwqJNf35Gsur7jXTzbihuOyOlpqZqPd6/fz/CwsJq/ZDt2LEj0tPT0bp16xp/1Lt1XFxc0L9/f8ybNw/Hjx9HdnY2duzYAQBQKBTo0aMHZs6ciSNHjsDNzQ0bNmwAAAQGBuLq1aua16qoqMDJkyfrfB+G9EmXBx54AJMmTcLWrVsxZMgQrFy5Ume7Nm3a1Birqo8bNmyIJk2a4Ny5czVePzQ0tM7+m4MzLLVgwq2Nk2pGhduT7Vfct8Cmd4FUE3c8fNZFVaiOamrxKODTRJVgqzOPRaH6eotHJXn5CxcuYPLkyRg3bhwOHz6MpUuXYsGCBbU+55133sEjjzyChIQEjBkzBl5eXkhPT0dKSgo+++wz/PLLLzh37hwee+wx1K9fH8nJyaisrER4eDhSU1Oxfft2PP744wgKCkJqaipu3LiBNm3aAFDlpkyePBm//vorWrVqhYULF+LWrVt1vo+6+lRdcXEx3nrrLTz77LMIDQ3FpUuXkJaWhmeeeUbn/SdOnIgXX3wRnTt3Ro8ePbB27Vr8+eefmjwcQFVIbsKECfD19UV0dDRKS0tx8OBB5OfnY/LkyXW+B1MxYKlFiL8ndmca1o5kQnNy749AUa7492cirf0b+BHg7ALsXWL8cwsvAp91BUb/Jn6/bJ2TMxD9sWo3EBTQDlr+zi2JnitZPZb4+HgUFxeja9eucHZ2xsSJEzXbhPVp164dfv/9d/zrX/9Cz549IQgCWrVqheeeew4A4Ofnhx9++AEzZsxASUkJwsLCkJiYiLZt2+LUqVPYtWsXFi9ejMLCQrRo0QILFizAwIEDAQAvvfQSjh07hvj4eLi4uGDSpEno06dPne+jrj5V5+zsjLy8PMTHx+P69esICAjAkCFDMHPmTJ3tn3vuOZw9exZvv/02SkpK8Mwzz+D//u//sGnTJk2bMWPGwNPTE5988gneeusteHl5ITIyssY2bbEpBGP3WMlQYWEhfH19UVBQAB8fH9HuW1xWgTbvb4a7s4B5XSvw9gFnlFYoarQ79UE0c1hEUF5ejuTkZMTExOjdclcrKWZUfJsDAa2BBi3tZkbF7HF2FPfKgH/3AG4aXwW03MkDye2/tLsxLikpQVZWFkJDQ3XuRjFI+kbVbqGqCbg+TVXBSsRTRt2qsrIShYWF8PHxgZOT/gyH3r17o0OHDli8eLFpfXZwho5zbfR97xjz+c0Zlloo3ZwRFRGEXRn6k8CiIoIYrFiTFDMqXPIhAHBxAyakAWuHAZlbTLvHjg+BAbp/k3VYEU+pti7LvNItyQ8DljqsiO+CV1anArhW42tREUFYEd/F8p0iFSlmVLjkQ9WZU6sl9XOg379UwQ/d5+QMhFr/MD2yLQxYDLAktiOSk5PxfJdgnMsrQYi/J6bGRHBmxdKkPOuHW5OpNrGJwJ2bwHwTdkF83kNVXI6sZufOndbuAomAAYsR3hsUYVfr0TZFysq0nFUhQ5haqyX3tKqM/4gkafpF5CBYh4XkS8qzfjz9WZGWjBebaNr5Q2e2qHJh7IQd7NUgCxPje4YzLCRPnFEhuYpLAtYNA04bmYibuUX1vOG2O9OinmEuKiqCUsnlUzKc+sgCc1YpGLCQfJT9fYT5p+2AO7WdOWICT3+g7T+584fEMTwJWPe86tRmY5zeYtPLQ87OzvDz80NOTg4AwNPTU+t8HkurrKxEWVkZSkpKTN5uS3UzZ5wFQUBRURFycnLg5+dnVvl+BiwkD4mxQOYOoP2XQNFNce/NGRWSwvD1wJF1wE+vGPc89fJQnG0GLY0aNQIATdBiTYIgoLi4GEql0qqBk70TY5z9/Pw03zumYsBC1lO9hoqTiYWodOGMClnCw8OBtBXAlcPGPS/TdoMWhUKBxo0bIygoSOdhfZZUXl6OXbt24bHHHuOGCAmZO86urq6iHIzIgIWsgzkqZC9e/k21zHPGhJwWG18ekvp0XkP6cO/ePXh4eDBgkZBcxpkBC1mG1DVU7Kx8PtmYEUmmBS1ntqiCdwbYRHViwELS42wKOYIRSaaV8c9IVgX0DLSJasW0apKGpoZKKwlqqASwhgrJU1wS8IAJdVqWdBK/L0R2hjMsJD7OqJAjG54EfDPcuOfcuWyzSbhElsIZFhKXFMEKZ1TI1jz3tfHPyfw7n4WIdOIMC5lPnVB74zSQvUvce3NGhWxZq35A5q+Gt2c+C5FeDFjIPJLMqDQA2sRwxw/ZvmFfA5+2Vy35GGpZV2DSCen6RGSjuCRExpPyUEIAmHgcGLSAwQrZhwmHjGtfcMGuDkokEgsDFjJOYiwwpxGQ9h+gKM/8+/k2B1r1VeWovHHG/PsRyY2bUrW0aQx1UTki0uCSENVOyoJv1fNTrFzmm0gysYnGF5ZjUTkiLQxYSD9J8lN4xg85KFMKyzEJl0iDS0KkTcqCb+ExwNvnmJ9CjisuSbUMaozPukjTFyIbwxkWuk+KGZXQXkBAGGdUiNTGH1DlgRmq8KJNH5JIJBbOsJCKFMFKeAwwciNnVIiqMiUJ9wyLyhExYHFk6uWfr5+SJlhhsiCRbrGJQJiRZw6p81mIHBSXhByVVCX02w7m8g+RIeKSgPkRLCpHZCAGLI5Cyu3JAGdUiEwx4ZBx+SwFF5jPQg6LAYsjkGI2xbc5ENAaaNCSMypEplLnsxjz7/PMFp7sTA6JAYs94mwKke2ITTS+PksmgxZyPAxYDFBSVgEAeHnNQTRr4I2pMRFQujlbuVd6SDGbArDgG5GU4pKMr4SrLt/P5SFyENwlVJuyYuxYMAKH5qm2IB48ew1r9l9Am/c3Y+zqNCt3TgepghUWfCOS3ggTispxuzM5EAYs+iTGQpjTCH1v/4wezn8CAA65/R++cFkAAEhJz5FH0CLV1uSQXqoDCade4/IPkaWMP2D8c7jdmRwEAxZdEmMhZCQDQs0vPe58SCtoKf57ucgqqp6cnPW7ePcNjwFeZME3IoszpagcoMpZI7JzDFiqKyvWzFQoFNpfUj9+3PkQ3FAGAJiTnG7Zvv36BrDmn8CiSOmWfzijQmQ9sYlAayOLyh3/Tpq+EMkIk26r+/s3FYWeL6uDln+5fIPp915Cdl6RZfolVX4KtycTyY+xJzuX3uKuIbJ7splhuXz5MkaMGAF/f38olUpERkbi4MGDlu/IzXMGNQtRXFf9199Tyt6oSJlMO+kE8MIGLv8QyU1cknHl+9VbnYnslCxmWPLz89GjRw/06dMHmzZtQmBgIDIzM1G/fn3Ld6ZBS+DsjjqbZQsNAQBTYyKk6Ye6lsqN00D2LnHuydkUItsSlwR8FAKU5hvWPnMLsG4YMJwzLWR/ZBGwfPzxxwgODsbKlSs110JDQ63TmagPVUmsegh/J+LOvjcCURFB0tRjkerkZOamENmeds/U+jOphtN/b3Xmv3eyM7IIWDZu3IgBAwZg6NCh+P3339G0aVO8+uqrGDt2rM72paWlKC0t1TwuLCwEAJSXl6O8vNy8zihcgHrBwN0bmkvlTh6a/6oDlsSAlWgX+735r1dWAvw2C8jPBuqHALcuAed2AH+/plma/wMIaAX0mQa4eQDm9lVi6rE0e0ypVhxn6Yk6xr1nAIe+Me45mTuAu7dV/+7tGL+XLUPKcTbmngpBEHRs3rUsDw/VP6rJkydj6NChSEtLw8SJE/H5559j5MiRNdrPmDEDM2fOrHF93bp18PS0QE4JERERma2oqAjDhw9HQUEBfHx8am0ri4DFzc0NnTt3xt69ezXXJkyYgLS0NOzbt69Ge10zLMHBwcjNza3zDddpy7+Aw19rXSp38kBK5BJEnZgA18qS+1/oOBIYMNvwe1edTck7BxReNK+v+oQ9Djz7X2nuLaHy8nKkpKQgKioKrq6u1u6O3eI4S0+SMf52JHBuu3HPadUPGPZ13e1sFL+XLUPKcS4sLERAQIBBAYssloQaN26MiAjt5NU2bdrg+++/19ne3d0d7u7uNa67urqaP5j5mUDVoKTq/StLtAKWiht/wdnQ15Nqpw9gd8m0ovx/pDpxnKUn6hiPWGfCIYm/At/G2f15Q/xetgwpxtmY+8liW3OPHj2QkZGhde306dNo0aKF5TvToKXBTXOunTesoZTBCrcmEzkOY7c6AzxviOyGLAKWSZMmYf/+/ZgzZw7OnDmDdevW4csvv8T48eMt35le/9JVkV+nGxW++r8o1Rk/ABDKc36IHFZcEvCAkUELzxsiOyCLJaEuXbpgw4YNePfdd/HBBx8gNDQUixcvRlxcnOU78/tsvVVuq7vrq2c2RuoZFQYpRI5teJLqeI6CC4Y/Z1lX1WwskY2SRcACAE888QSeeOIJa3cDuJFpUDNBALoW7gA2TgRQCRRcUi0n5V9UTcFKgcEKEamNP6A6/NRQBRdUsyxcNiYbJZuARTbuXDWomUIBOJcXAIdX3b9oQIVcg4UNAPyCVUcF2EkyLRGJSH2yszGzuZxlIRvGgKWayuJC6yf2cCaFiAwRm2jczqGCCzwkkWyW1T+b5ab87i3Lv6hvc6BVXybSEpHx4pIA76aGt1efN0RkYzjDUk0hlAiE7joskuBsChGZa8Ih4/JZTm8Bvhlm9/VZyL5whqWas37/kP5FQrgtmYhEpM5nMcaZLarlISIbwYClmvZjPocg3D+VWXThMcCLG1nkjYjEFZuoWl42RiaLypHtYMBSjdLbG0e9HgUgQdDC5R8iktL4A8Y/h0XlyEYwYNHh4bc3aYIWk4UNUC37MJmWiCzFlKUhQLXdmUjmmHSrx8Nvb0LxnTs4ujIByiJVbRaDJ1w4k0JE1mLsVmdAtd05MZY/t0jWGLDUQuntje6vrUJ5eTnOJydD8cYZYOcMVFw/jZzrF5Ar+KHYKxgPN/eD653LLPBGRPIQl6TaunzaiKBFvTTEn18kUwxYjOHmgbF5zyPldM79a4UArgJREUFYMaiL1bpGRKRleJJq67IxR4V81gWYfFK6PhGZgTksRpiQeBgp6Tk6v5aSnoOxq9Ms3CMiolqMSFLl0xmq8CKLypFsMWAxwo6MG7V+PSU9B8VlFRbqDRGRAeKSjNvufJr1WUieGLCIbE5yurW7QESkzdjtzpl/V8IlkhEGLCLLziuydheIiLS5KY0vKsdKuCQzDFhEFuLvae0uEBHVZEpROc60kIwwYBHZ1JgIa3eBiKgmU4vKcaaFZIIBixH6hgfW+vWoiCAo3Zwt1BsiIiPFJhq3a0gtcwt3D5HVMWAxwpLYjoiKCNL5taiIIKyIZx0WIpK5uCTgAROCltM8KJGsi4XjjLQivguKyyowJzkd2XlFCPH3xNSYCM6sEJHtGJ5kfPl+gNVwyaoYsJhA6eaMWYMjrd0NIiLTxZlQCRdgNVyyGi4JERE5KmMr4QKshktWw4CFiMiRxSUBrY0MWk5zuzNZHgMWIiJHZ8pMC7c7k4UxYCEiItVMi0+wcc9hYTmyIAYsRESkkmDCifOcaSELYcBCREQqplbDzWSNFpIeAxYiIrrP1Gq46hotRBJhwEJERNpMrYb7Gat9k3QYsBARUU3DTdjuXHiRSbgkGQYsRESkG7c7k4wwYCEiIv1MKSzH051JAgxYiIiodiNMqNHCargkMgYsRERUN9ZoIStjwEJERHUzp0YLgxYSAQMWIiIyjKk1WpjTQiJgwEJERIYzJQkXUOW0sBoumYEBCxERGceU7c4Aq+GSWRiwEBGR8VgNlyyMAQsREZmG1XDJghiwEBGR6VgNlyyEAQsREZknzoSgJXMLsPAh5rSQwRiwEBGR+eJMqIZbeBGY04i7h8ggDFiIiEgcplTDBVS7hxi0UB0YsBARkThMrYYLcMsz1YkBCxERicfUargAtzxTrRiwEBGRuEythlt4EZgfwZkW0okBCxERic/Uarh3LjMRl3RiwEJERNIwtRouwERcqoEBCxERSWd4knmJuHduitsfslkMWIiISFqxicCbWaY9d34oZ1oIAAMWIiKyBO8G5s208Pwhh8eAhYiILMOcLc9nWMrf0TFgISIiyzF1yzNwv5Q/D050SAxYiIjIskaYkYgLqA5O/CRMvP6QTWDAQkRElmdOIi4A3M0BEl8Qrz8kewxYiIjIOsxJxAWA7N9U/y0rEac/JGsMWMxQXFaBaT+ewAtfpWLajydQXFZh7S4REdmW2ETzghYAWNCaeS0OwMXaHbBVY1enISU9R/N4dyawZv8FREUEYUU8D/AiIjJYbKJq98+STqrS/KZQ57W8lSlu30g2OMNigurBSlUp6TkYuzrNwj0iIrJxbkrgzXTTtz0DqryW1c+I1yeSFQYsRiouq9AbrKilpOdweYiIyBRxZu4gOreN9VrsFAMWI81JThe1HRERVRObCEy9Brh4mPZ81muxSwxYjJSdVyRqOyIi0sFNCbx3HWjZ3/R7ZG4BPmzI2RY7wYDFSCH+nqK2IyKiWsR/DzxgRl7LvRLVbAsPULR5DFiMNDUmQtR2RERUh+Fm5rUAqgMUGbTYNFkGLHPnzoVCocDrr79u7a7UoHRzRlREUK1toiKCoHRztlCPiIgcgLl5LYAqaLlzU7w+kUXJLmBJS0vDF198gXbt2lm7K3qtiO+iN2hhHRYiIomIkdcyPxRYFMm8Fhskq4Dlzp07iIuLw4oVK1C/fn1rd6dWK+K74NQH0XjhkeboGRaAFx5pjlMfRDNYISKSmrl5LQUXuIvIBsmq0u348eMxaNAg9O/fHx9++KHedqWlpSgtLdU8LiwsBACUl5ejvLxc9H6p71n93i4K4P1BD1a5Uony8krRX99R6BtnEhfHWXocYwsYuhbl370MACh3MnGZ6OzvwOwWwOsnATczlprsnJTfz8bcUyEIgiB6D0ywfv16zJ49G2lpafDw8EDv3r3RoUMHLF68uEbbGTNmYObMmTWur1u3Dp6e3J1DRERkC4qKijB8+HAUFBTAx8en1rayCFguXryIzp07IyUlRZO7UlvAomuGJTg4GLm5uXW+YVOUl5cjJSUFUVFRcHV1Ff3+pMJxtgyOs/Q4xpahNc4/jAHObTfvhl6BwIQj4nTOjkj5/VxYWIiAgACDAhZZLAkdOnQIOTk56Nixo+ZaRUUFdu3ahc8++wylpaVwdr6/68bd3R3u7u417uPq6irpDwep708qHGfL4DhLj2NsGa6urnAdsU61A2h+qOk3un0RSHxelSNDNUjx/WzM/WSRdNuvXz+cOHECR48e1fzp3Lkz4uLicPToUa1ghYiISCfvBubXazm3DVjYlruIZEgWMyz16tXDQw89pHXNy8sL/v7+Na4TERHpFZuoKhCXkWz6PQovqXYR+TYHxh9Qbacmq5PFDAsREZFo1EXmfJubdx/19udvuP1ZDmQxw6LLzp07rd0FIiKyVW5KYNIJ1dLOxyFARYnp9zqzBfgkDHgrU7TukfE4wyKC4rIKTPvxBF74KhXTfjyB4rIKa3eJiIgAVeAy7bp5heYA4G4OsPoZcfpEJpHtDIutGLs6DSnpOZrHuzOBNfsvsEQ/EZGcDE9SzbbMC1Gd4GwKdUJuwkHmtVgBZ1jMUD1YqSolPQdjV6dZuEdERKSXGGcRqRNyeR6RxTFgMVFxWYXeYEUtJT2Hy0NERHIT/z0QZuYSkTohl4GLxTBgMdGc5HRR2xERkQXFJZlfswW4H7gkxpp/L6oVAxYTZecVidqOiIgsTKztz4Cq7gtPf5YUAxYThfgbdsiioe2IiMgK1Nufp14DXMw8sTlzC7DwIS4RSYQBi4mmxkSI2o6IiKxIjIRcACi8yNwWiTBgMZHSzRlREUG1tomKCILSjecgERHZDDEScgHmtkiAAYsZVsR30Ru0sA4LEZGNEishF2Bui4hYOM5MK+K7oLisAnOS05GdV4QQf09MjYngzAoRkS2LTVQt6SzrqpotMUfmFmB+BDDhEAvOmYEBiwiUbs6YNTjS2t0gIiIxVT2PyNzA5c5l1RJReIwqGCKjmbwkdPjwYQwfPhwxMTF45513kJWVJWa/iIiI5KHqTiKfYPPuxSUik5kcsDz33HN44oknMHv2bDzwwAMYMmQItm7dKmbfiIiI5MNNCUw+af5BiuolIu4iMorJAUv9+vUxfPhwPPzwwxg9ejS2bduGt956S8y+ERERyc9wEZJy1UtE3EVkMJMDllatWmHBggUQBAEA4OfnJ1afiIiI5E2sKrkZyazZYiCTA5bS0lIsX74czZs3R3R0NB566CH0798fly9fFrN/RERE8iRWbgsPUjSIwbuE9u7dCx8fHzz00EMAgB9++AEAcPfuXZw4cQLHjx/H8ePHERsbiytXruDMmTPS9FjGuL2ZiMgBqXNbvhkGnNli+n3UgQt3EulkcMAyfvx4JCQkaAIWtWvXrqFt27Z45JFHRO+cLRm7Og0p6Tmax7szgTX7L7CAHBGRoxiRpMpJyUg27z4ZyargZ0SSOP2yEwYvCWVkZKB37941rm/btg2xsY6dNFQ9WKkqJT0HY1enWbhHRERkFercFnO3P5/Zwu3P1RgcsPj4+CA/P7/G9Z49e2L//v2idsqWFJdV6A1W1FLSc1BcVmGhHhERkVWpl4hai7D9mUGLhsEBS3R0NObPn1/zBk5OKCsrE7VTtmROcrqo7YiIyE6MEGH7c+YW1fIQGR6wzJo1C7///jueeeYZnDhxAgBQUlKCjz/+GO3atZOsg3KXnVckajsiIrIjYmx/5vIQACMCluDgYOzfvx/FxcVo3749lEol6tWrh59//hmffPKJlH2UtRB/T1HbERGRnam6/dnUwCVzC7DOsYMWo+qwtGjRAsnJycjOzsa3336LjRs34vTp0+jevbtU/ZO9qTERorYjIiI7pQ5cTF0mOu3Yy0MmndbcvHlzNG9uZnU/O6F0c0ZURFCtibdREUGsx0JERCqxiaoCcUs6qUr0G0O9PBTneFueTa50S/etiO+CqIggnV9jHRYiIqrBTQm8mQ6EmbCTyEETcU2aYaGaVsR3YaVbIiIyTlySasYk08gKuQ4408KARURKN2fMGhxp7W4QEZEtiUsyrax/5hbV+UPjD6hmbOwcl4QkUlxWgWk/nsALX6Vi2o8nWDiOiIj0G5Fk2vKQ+vyhRPuvOM8ZFgnwXCEiIjJaXJJq6/JpEw5QzEi2+yUizrCIjOcKERGRyYYnmV7S386TcRmwiIjnChERkdlMXR4CVHkwdro8xIBFRDxXiIiIRBFnxkxLRrKqzoudYcAiIp4rREREojFnpmVZV3H7IgMMWETEc4WIiEhUcSae+Fxwwe7yWRiwiIjnChERkehMPfHZzk55ZsAiIvW5QrXhuUJERGQ09cGJxi4RZdpPEi4DFpHxXCEiIpKMKcm4dpKEy8JxEuC5QkREJJkRSaqS/AUXDH/Osq6qGRobxoBFIjxXiIiIJDP+gKokv6EKLqiWhmITpeuTxLgkREREZGvclMbvHrLxpSEGLBbEAxGJiEg0sYnG57PYcH0WLglZCA9EJCIi0Y1IUm1dzjTwwEQbXhriDIsF8EBEIiKSTFyScTVabHRpiAGLxHggIhERSW78AePa2+DSEAMWifFARCIikpyxSbjqpSEbwoBFYjwQkYiILCI20a6XhhiwSIwHIhIRkcXY8dIQAxaJ8UBEIiKyGDteGmLAIjEeiEhERBZlp0tDDFgsgAciEhGRRRm7NJTynjT9EBELx1mIvgMRAWDajyd4SCIREYlHvTSUkWxY+5vnpO2PCBiwWFD1AxFZ/ZaIiCQTm2j4qc65Z6Tvj5m4JGQlrH5LRESSM3RpqOCC7PNYGLBYAavfEhGRRbgpDU/AlfkWZwYsVsDqt0REZDEBrQ1rJ/MtzgxYrIDVb4mIyGIatDS8rYy3ODNgsQJWvyUiIouJ+tC49jLd4syAxQpY/ZaIiCzG2Oq3f26Qri9mYMBiBax+S0REFmVM9duiPFnmsjBgsRJWvyUiIosypvqtDHNZWDjOilj9loiILMbY6rcp7wGDFkjbJyMwYLEyVr8lIiKLiU0E5rUCinLrbiuzcv1cEpIRVr8lIiLJtR1sWDuZleuXRcDy0UcfoUuXLqhXrx6CgoIwePBgZGRkWLtbFsXqt0REZBGGbnOWWbl+WQQsv//+O8aPH4/9+/cjJSUF5eXlePzxx3H37l1rd81iWP2WiIgswphy/TKqySKLHJbNmzdrPV61ahWCgoJw6NAhPPbYY1bqlWWx+i0REVlMQGvDTnH+80fg8bmSd8cQsghYqisoKAAANGjQQOfXS0tLUVpaqnlcWFgIACgvL0d5ebno/VHfU4p7q7X098CBc4JB7aTshzVZYpyJ42wJHGPL4DiboX4Y4LS37nYld1D+3cuA52BJP18NoRAEoe5PSQuqrKzEU089hVu3bmHPnj0628yYMQMzZ86scX3dunXw9GQ5eyIiIltQVFSE4cOHo6CgAD4+PrW2lV3A8sorr2DTpk3Ys2cPmjVrprONrhmW4OBg5Obm1vmGTVFeXo6UlBRERUXB1dVV9PurTUg8jB0ZNwxu3zc8EEtiO0rWH0uz1Dg7Oo6z9DjGlsFxNtN3LwGZW+tsVu7kgZTIJZKMc2FhIQICAgwKWGS1JJSQkIBffvkFu3bt0husAIC7uzvc3d1rXHd1dZX0m1bq+y+P71br1ubqNqXn4tXEo3ZXn0XqcSYVjrP0OMaWwXE2UewaYF5LVSl+A0gxzsbcTxa7hARBQEJCAjZs2IAdO3YgNDTU2l2ymhXxXXDqg2i88Ehz9GilO4enKm51JiIik7X9p7V7YDBZBCzjx4/HN998g3Xr1qFevXq4du0arl27huJi+ez/tiR19duWgd4GtedWZyIiMomhNVkAoKxEun4YQBYBy/Lly1FQUIDevXujcePGmj/ffvuttbtmVdzqTEREkjKmJstvs6TtSx1kkcMis7xf2Qjx98TuTMPaERERmcTQmiz52ZJ3pTaymGEh3dQnN4vVjoiIqIYGLQ1rl2fdwxAZsMiY0s0ZURFBtbZp4OmKOcnpTLwlIiLTGJrHUnjRqmcLMWCRuRXxXWoNWm4WlWPN/gto8/5mnuZMRETGs5GzhRiw2ICqW50beOnfs56SnsOghYiIjBfQ2rB2N623LMSAxUYo3ZwxNSYCN+/Wfu4C67IQEZHRDM1jMbSdBBiw2BBD662wLgsRERnF0DyWgovS9qMWDFhsCOuyEBGRJNyUQNiAutud3mK1xFsGLDbE0HorZ67fxrQfT3BpiIiIDOcXbFg7KyXeMmCxIYbWW7laWMqdQ0REZBxDE2qtlHjLgMWGGFKXpTruHCIiIoPIPPGWAYuNqasuiy7cOURERHUyNPE23zqJtwxYbFDVuiyNfNwNeg53DhERUa0MTbw9Y53EWwYsNkrp5oxZgyMR1rCeQe1/PnaVibhERFQ7GSfeMmCxcYbuHLpVzBL+RERUBxkn3jJgsXGmnNTMRFwiItLJwITanOw/Je5ITQxYbJwpO4cAJuISEZEOUR9CACAI+psIAhB47zqK79yxWLcABix2wZSdQwATcYmIqBo3JfJdG0Gh0N9EoVD9OboywXL9AgMWu1F155CfUv+JzlUxEZeIiKq7WNnQoHaetyybx8KAxY6odw492b6xQe2ZiEtERNXVq7wpajuxMGCxQ0zEJSIiU91x8Re1nVgYsNghJuISEZGp7vqEGtTOv/yyxD3RxoDFTjERl4iITNFh1GcQhLp3CjURblh0pxADFjvGRFwiIjKW0tsblxWBstspxIDFzjERl4iIjHXT3bAS/V53zkvck/sYsDgIJuISEZGh7nq3ELWdGBiwOAgm4hIRkaHC4z6tNY9F/bXwuE8t1icGLA7E1ETcPvN/Y14LEZEDWbTrArZWdAJQM2hRP95a0QmLdl2wWJ8YsDgYUxJxrxWWMq+FiMiBZOcVYdy9NzRBS3VbKzph3L03kJ1XZLE+MWBxQMYm4lbFvBYiIvsX4u8JABh37w10KvscAPBHRVt8fa8/wktXYdy9N7TaWQIDFgdmSiIuwLwWIiJ7N6l/uObvZVDNxo+7NxnT772EMrjpbCc1BiwOzNREXIB5LURE9mzRtgxR24mBAYuDMzURl3ktRET2y9DcFOawkEVVTcRt5ONu9POZ10JEZF8MzU1hDgtZnDoR97c3+5j0/JT0HAz/ci+XiYiI7IChOY6m5kKaggELaTEnr2XvuXwuExER2QFDPguiIoKgdHO2UI8YsJAOpua1VMVlIiIiEhMDFtKpal5LYxPyWgBufyYislXFZRVISc+ptY2lf8YzYCG91HktO0zMawGAvtz+TERkc+Ykp4vaTgwMWKhO5uS1XOX2ZyIim8NtzWSzmNdCROQ4uK2ZbFrVvJburRqYdA9ufyYikj85bmt2sdgrkV1Q57UAwNjVaXUmZemy91y+Zgt0VEQQVsR3EbubRERUi+KyCsxJTkd2XhFC/D0xNSZCa4uy0s0Zgd5uuHGnTO89LL2tmQELmWxFfBeTgxa1lPQc/GPuDoQGeqGlvwc6KkTsIBER1VD95/buTNT4BXLs6rRag5VAbzeL/7LJJSEyi7ll/QHg0q1i7M7Mxfq0iwCACYmHxewiERH97aWV+/X+kqnOMzRkS/ONO2UWX9ZnwEJmM7esf3U7Mm4wOZeISGQvrTqAHRl5tbZJSc/BBz//adD9LLmlGWDAQiIyZ/tzdSnpOYhbsY/JuUREIhi7Og07/rphUNsfjlwyqJ0ltzQDDFhIZGJsf1b74+xN1nAhIjKTIUs8VZXeEwxqZ8ktzQCTbkkCK+K7aDLQz924gz/O3jT7nlWTc3VltBMRkW5SLd1YckszwICFJCLG9ufqLt0q/jtBt2ZGOxER6fbr8atGP6eBpytuFpXr/bqltzQDXBIiCxBzmaiqlPQc/OPj7Xjhq1TmuhAR6VBcVlFr4KHPoHaN9f7cttYvi5xhIYuoukyUnVeErNw7uJRfYvZ9L+WX4FJ+CWddiIh0iFr4u0nPUy+7F5dVYG7ySQDZeL5LMKbEPGS15XgGLGQxVZeJAPGWiqpS5bpsR2igN3NdiMihjV2dhku3io1+njWWewzBgIWspuqsy9mc2zh03vzkXAC4dKsEl25x1oWIHJexO4PU+ob7a1W7TUnPgbuzgI5dgfVpF/H1/ktcEiLHVHXW5ZXVqQCuif4anHUhIkcTtXCnSc/bkZGnKSNRV0VcSwctDFhINpbEdkRycrIk964+69I3PBBN6yv1HvxFRGSLissqELVoJy7dMj1H0JCZmZT0HBSXVfDwQ3JsB6f2x7yUTFGTc6vbkXG/4iOXjojIHkiRF1ibOcnpWnmJUmPAQrLjYYHkXF3U26RDA7h0RES2QZ0H+OvxqyZtXzaHpUvzM2Ah2auxJfrGHbOmO2uja5v0kuc7al6bgQwRWUvVn4Mh/p64nF+EHRm5VusPS/MT6WCJLdG6pKTnoM37mzWPNTkwDwaiqR9zYIhIWpoZlBNXcfPu/RmU3ZlW7NTfWJqfyABSFaIzVNVTT5nIS0Ri0hekyIk1arUwYCGbZa1ZF310JfJyJoaIdKm+vDOpfzgWbcuwSJCiLrlf28/L2tqwDguRmaw966KLrpkY5sUQ2T99AYmun03qnw2WUDXY0PdLXtU2LM1PJBG5zbroojcvhktKRDZDrgGJPs38lEiZ3EvrZ0r1X/J0/dxRujnjvUERSE7OxnuDIuDqar2fSQxYyK7p+gdp7cx6fXQuKYX7o2l9LwYxRCLR9QENoM5rcg9IalPbEk71X/LkjAEL2T1d/yAttU3aXDsy8gDkAah7SQmo+UOXwQ3ZK9VSRTo6KoAPf03HlJiHANQeeFT/t64r0DD0mtw18HLFoMjGdvVzgAELOSRbWDrSR9+SUlW1BTdvR4VZustEBjNkBkQ9S1r9UL6qbDHIEIu9Vu2WVcCybNkyfPLJJ7h27Rrat2+PpUuXomvXrtbuFjkAfWu5E9YftplARhddwU1S2nnM66r6rfRcXonWbIyh0+X28hsbiU9fToRBgcitYp2J6lXJof6InIiR+1ZXHotcyCZg+fbbbzF58mR8/vnn6NatGxYvXowBAwYgIyMDQUFB1u4eOQBdS0e6c2Du/r1UY9vWp11EaYVC86EQ6O2GG3fKNF83duaGwY08GBN0it22en6Ysd9bZJgGnq4Y1E6c5Z7qs8tyPltNIQiCYO1OAEC3bt3QpUsXfPbZZwCAyspKBAcH47XXXsOUKVO02paWlqK0tFTzuLCwEMHBwcjNzYWPj4/ofSsvL0dKSgqioqLg6uoq+v1JxZbGuaSsAvNTMnDhZhGaN/DE1fwi/H7GNoIYdycBszpXYtpBJ5RWKiR/vb7hgZj3THut8XozKhwebs41xvHNqHAAMLith0yDobq+l/W9F3PH4+3vj2klb9cmwMsNuXfL6m5oZFtLsvT3sjXVV7piwEONRP2+n5B4uNbvl77hgVgS21HSn82FhYUICAhAQUFBnZ/fsghYysrK4Onpie+++w6DBw/WXB85ciRu3bqFn376Sav9jBkzMHPmzBr3WbduHTw9LXu2AREREZmmqKgIw4cPNyhgkcWSUG5uLioqKtCwYUOt6w0bNsRff/1Vo/27776LyZMnax6rZ1gef/xxzrDYMHscZ3N/C5aCvf9Wqm82wNjr5tzb3sdYLmx9nJv6KtEiwBPNG3hifO/WWLbzjMVmED/8NR3r0y4a1LbqOPcIC8KS2I6i9aOwsNDgtrIIWIzl7u4Od3f3GtddXV0l/aCT+v6kYk/j7OrqipmD22tdWx7fTRa1YUorFSitsL0f8nW5XFgOoOb7Mva6GPe21zGWG7mOczM/D4QGeteo46Ivv6v6zwopncsrMXrMSisV2JSei1cTj4qW32LMz3pZBCwBAQFwdnbG9evXta5fv34djRo1slKviKRTV20YXTsmiEhemvkpERroZXBAIqcCbSH+niYnOqek56C4rMIxDz90c3NDp06dsH37dk0OS2VlJbZv346EhATrdo7IQgwJYuxhqzWRtTWr74HQAG+TK93aQkBSl6kxEWbVqZmTnG7x9yuLgAUAJk+ejJEjR6Jz587o2rUrFi9ejLt372LUqFHW7hqR1Ri61VrfD10GN2Sv+oYHIri+O6oeygfUHXjo2wqs68O3+jVbCkjqonRzRlREkMk/H7LzikTuUd1kE7A899xzuHHjBt5//31cu3YNHTp0wObNm2sk4hKR/vM/DAluruffBlBzqal6rQwisRjzvaXr/CxAd9BRXl5e41A+Q/5dkMqK+C4mV/kO8bf8jlzZBCwAkJCQwCUgIpFVD25UP+ST8XyXYJMr3XLmxvYZE0QYF3DorrxqbhVlBh3SqP5LTTM/JRIN2D2k/v9nSbIKWIjIclS/lWpn6Jszc2NscCPVB6atEmM8jK1ALFVVXF0M/d4iy6v+/yb3bmmtv5BERQRZpYI1AxYiMom5wY25H5j6jkjQ92Fu7HVdxLiHvqBCjABC/SFiaGBgTBDBgMNx1LZUZM2S/bKodGuuwsJC+Pr6GlQpzxTqKfSYmBi7qQ8iRxxny7CncTYmEDL2OmBcoFD1ekt/D3RUZKNv/wGYl5LJs5UkYk/fy3Kk/p6+dPMOBvvnoG//Aajn5SHqaxjz+c0ZFiKyWfp+6xfruqn3UCeDeuhpS2QL1N/T6sDQ2md3MWAhIiJycLUtMcoFAxYiIiIHVj1fZXcmsGb/Bavmq+jiZO0OEBERkXXUVoclJT0HY1enWbhH+jFgISIickDFZRV11lNKSc9BSVmFhXpUOwYsREREDmhOcrpB7eanZEjcE8MwYCEiInJAhp4HdOGm5c8N0oUBCxERkQMy9Dyg5g0sf26QLgxYiIiIHJCh5wG9GRUucU8Mw4CFiIjIASndnBEVEVRrm6iIIKsXjFNjwEJEROSgVsR30Ru0NKvvgZLySnz4q2HJuVJjwEJEROTAVsR3wakPovHCI83RMywAzfxU5wVdyi/B7sxcrE+7CACYkHjYmt1kwEJEROTo1OcGebg64dKtEp1tdmTcsGohOQYsREREZHAhuWIrFZJjwEJEREQGF5IztJ3YGLAQERGRwYXkDG0nNgYsREREZHAhOUPbiY0BCxERERlcSM7QdmJjwEJEREQAoNnSrE9URBCUViok52KVVyUiIiLZGLs6rc4dQn3DA7E8vouFelQTZ1iIiIgcWF3BStO/Z12WxHa0VJd0YsBCRETkoAypvXJZTyE5S2PAQkRE5KCsVVPFFAxYiIiIHJS1aqqYggELERGRg7JWTRVTMGAhIiJyUNaqqWIKBixEREQOSunmjKiIIIPaTkg8LHFvaseAhYiIyIGtiO+CvuH+dbbbkXHDaic1AwxYiIiIHF7T+l4GtbPmriIGLERERA5O7ic1AwxYiIiIHJ7cT2oGGLAQERE5PLmf1AwwYCEiInJ4huwW6hseaLWTmgEGLERERATVbqHaghZrH37oYtVXJyIiItlYEd8FxWUVmJOcjuy8IoT4e+LtqDDs2LbF2l1jwEJERET3Kd2cMWtwpOZxeXm5FXtzH5eEiIiISPYYsBAREZHsMWAhIiIi2WPAQkRERLLHgIWIiIhkjwELERERyR4DFiIiIpI9BixEREQke3ZROE4QBABAYWGhJPcvLy9HUVERCgsL4erqKslrEMfZUjjO0uMYWwbH2TKkHGf157b6c7w2dhGw3L59GwAQHBxs5Z4QERGRsW7fvg1fX99a2ygEQ8IamausrMSVK1dQr149KBQK0e9fWFiI4OBgXLx4ET4+PqLfn1Q4zpbBcZYex9gyOM6WIeU4C4KA27dvo0mTJnByqj1LxS5mWJycnNCsWTPJX8fHx4f/KCyA42wZHGfpcYwtg+NsGVKNc10zK2pMuiUiIiLZY8BCREREsseAxQDu7u6YPn063N3drd0Vu8ZxtgyOs/Q4xpbBcbYMuYyzXSTdEhERkX3jDAsRERHJHgMWIiIikj0GLERERCR7DFiIiIhI9hiwEBERkewxYNFh9uzZePTRR+Hp6Qk/Pz+DniMIAt5//300btwYSqUS/fv3R2ZmprQdtXE3b95EXFwcfHx84Ofnh9GjR+POnTu1PufatWt44YUX0KhRI3h5eaFjx474/vvvLdRj22TKOAPAvn370LdvX3h5ecHHxwePPfYYiouLLdBj22TqOAOqnx8DBw6EQqHAjz/+KG1HbZyx43zz5k289tprCA8Ph1KpRPPmzTFhwgQUFBRYsNfyt2zZMoSEhMDDwwPdunXDgQMHam3/v//9Dw8++CA8PDwQGRmJ5ORkyfvIgEWHsrIyDB06FK+88orBz5k3bx6WLFmCzz//HKmpqfDy8sKAAQNQUlIiYU9tW1xcHP7880+kpKTgl19+wa5du/Dyyy/X+pz4+HhkZGRg48aNOHHiBIYMGYJhw4bhyJEjFuq17TFlnPft24fo6Gg8/vjjOHDgANLS0pCQkFDnWR+OzJRxVlu8eLEk56DZI2PH+cqVK7hy5Qrmz5+PkydPYtWqVdi8eTNGjx5twV7L27fffovJkydj+vTpOHz4MNq3b48BAwYgJydHZ/u9e/ciNjYWo0ePxpEjRzB48GAMHjwYJ0+elLajAum1cuVKwdfXt852lZWVQqNGjYRPPvlEc+3WrVuCu7u7kJiYKGEPbVd6eroAQEhLS9Nc27Rpk6BQKITLly/rfZ6Xl5ewevVqrWsNGjQQVqxYIVlfbZmp49ytWzfhvffes0QX7YKp4ywIgnDkyBGhadOmwtWrVwUAwoYNGyTure0yZ5yrSkpKEtzc3ITy8nIpumlzunbtKowfP17zuKKiQmjSpInw0Ucf6Ww/bNgwYdCgQVrXunXrJowbN07SfvLXJRFkZWXh2rVr6N+/v+aar68vunXrhn379lmxZ/K1b98++Pn5oXPnzppr/fv3h5OTE1JTU/U+79FHH8W3336LmzdvorKyEuvXr0dJSQl69+5tgV7bHlPGOScnB6mpqQgKCsKjjz6Khg0bolevXtizZ4+lum1zTP1+LioqwvDhw7Fs2TI0atTIEl21aaaOc3UFBQXw8fGBi4tdnP9rlrKyMhw6dEjr88vJyQn9+/fX+/m1b98+rfYAMGDAAMk/7xiwiODatWsAgIYNG2pdb9iwoeZrpO3atWsICgrSuubi4oIGDRrUOmZJSUkoLy+Hv78/3N3dMW7cOGzYsAGtW7eWuss2yZRxPnfuHABgxowZGDt2LDZv3oyOHTuiX79+zMvSw9Tv50mTJuHRRx/F008/LXUX7YKp41xVbm4uZs2aZfBynb3Lzc1FRUWFUZ9f165ds8rnncMELFOmTIFCoaj1z19//WXtbto8qcd52rRpuHXrFrZt24aDBw9i8uTJGDZsGE6cOCHiu5A/Kce5srISADBu3DiMGjUKDz/8MBYtWoTw8HD897//FfNtyJ6U47xx40bs2LEDixcvFrfTNshSP58LCwsxaNAgREREYMaMGeZ3nCzKYebD3njjDbz44ou1tmnZsqVJ91ZP5V6/fh2NGzfWXL9+/To6dOhg0j1tlaHj3KhRoxoJXffu3cPNmzf1To2fPXsWn332GU6ePIm2bdsCANq3b4/du3dj2bJl+Pzzz0V5D7ZAynFWfw9HRERoXW/Tpg0uXLhgeqdtkJTjvGPHDpw9e7bGTsRnnnkGPXv2xM6dO83ouW2RcpzVbt++jejoaNSrVw8bNmyAq6urud22CwEBAXB2dsb169e1rl+/fl3vmDZq1Mio9mJxmIAlMDAQgYGBktw7NDQUjRo1wvbt2zUBSmFhIVJTU43aaWQPDB3n7t2749atWzh06BA6deoEQPUDvLKyEt26ddP5nKKiIgCosVPF2dlZMyvgKKQc55CQEDRp0gQZGRla10+fPo2BAwea33kbIuU4T5kyBWPGjNG6FhkZiUWLFuHJJ580v/M2RMpxBlQ/jwcMGAB3d3ds3LgRHh4eovXd1rm5uaFTp07Yvn07Bg8eDEA1y7p9+3YkJCTofE737t2xfft2vP7665prKSkp6N69u7SdlTSl10adP39eOHLkiDBz5kzB29tbOHLkiHDkyBHh9u3bmjbh4eHCDz/8oHk8d+5cwc/PT/jpp5+E48ePC08//bQQGhoqFBcXW+Mt2ITo6Gjh4YcfFlJTU4U9e/YIYWFhQmxsrObrly5dEsLDw4XU1FRBEAShrKxMaN26tdCzZ08hNTVVOHPmjDB//nxBoVAIv/76q7XehuwZO86CIAiLFi0SfHx8hP/9739CZmam8N577wkeHh7CmTNnrPEWbIIp41wduEuoTsaOc0FBgdCtWzchMjJSOHPmjHD16lXNn3v37lnrbcjK+vXrBXd3d2HVqlVCenq68PLLLwt+fn7CtWvXBEEQhBdeeEGYMmWKpv0ff/whuLi4CPPnzxdOnTolTJ8+XXB1dRVOnDghaT8ZsOgwcuRIAUCNP7/99pumDQBh5cqVmseVlZXCtGnThIYNGwru7u5Cv379hIyMDMt33obk5eUJsbGxgre3t+Dj4yOMGjVKKyjMysqqMe6nT58WhgwZIgQFBQmenp5Cu3btamxzJm2mjLMgCMJHH30kNGvWTPD09BS6d+8u7N6928I9ty2mjnNVDFjqZuw4//bbbzp/ngMQsrKyrPMmZGjp0qVC8+bNBTc3N6Fr167C/v37NV/r1auXMHLkSK32SUlJwgMPPCC4ubkJbdu2tcgvjQpBEARp53CIiIiIzOMwu4SIiIjIdjFgISIiItljwEJERESyx4CFiIiIZI8BCxEREckeAxYiIiKSPQYsREREJHsMWIiIiEj2GLAQERGR7DFgISJZGTVqFN577z1rd4OIZMZhTmsmIvmrqKjAL7/8gl9//dXaXSEimeEMCxFZRGJiIpRKJa5evaq5NmrUKLRr1w4FBQUAgL1798LV1RVdunTRtLlw4QJGjhyJhg0bQqlUon379tizZ4/F+09E1sWAhYgs4vnnn8cDDzyAOXPmAACmT5+Obdu2YdOmTfD19QUAbNy4EU8++SQUCgUA4Pz58+jatSuKi4uxceNGHD9+HAkJCfDx8bHa+yAi6+BpzURkMb/88gueffZZTJs2DQsWLMDu3bvRtm1bzdcfeOABLFq0CIMGDQIAxMTEwM3NDT/++KOVekxEcsGAhYgsqmPHjvjzzz+xdetW9OrVS3P91KlT6NKlC3Jzc+Hh4YHz588jJCQER44cQYcOHazXYSKSBS4JEZHFbN68GX/99RcqKirQsGFDra9t3LgRUVFR8PDwAAAcPXoUbm5ueoOVb775Bl27dkVkZCQGDRqE0tJSAECPHj2QmpoKABg9ejQWLVok3RsiIothwEJEFnH48GEMGzYMX331Ffr164dp06Zpff2nn37C008/rXns6uqKe/fuoaioSOf9Bg4ciAMHDuDEiRNo0qQJdu7cCQCYNm0a5s6di4ULF8LJyQmTJk2S7D0RkeVwSYiIJJednY3u3btj4sSJmDJlClJTU9G9e3ccPHgQHTt2RE5ODpo1a4YrV64gICAAAJCXl4ewsDA8+eSTmDJlCgRBwK5du9CvXz+0bt0aH3/8Mb7//nuUlZXh4sWL2LBhg2aJqVOnTvDz88PmzZvh6upqzbdORCLhDAsRSermzZuIjo7G008/jSlTpgAAunXrhoEDB2Lq1KkAgJ9//hldu3bVBCsA4O/vj59//hmZmZno0qUL/vGPf2Djxo0ICgrCqlWr8Ndff2HXrl04duwY6tevj4iICABAWloabt68CV9fXwYrRHaEheOISFINGjTAX3/9VeN61eJwP/30E5566qkabXr06IG9e/fWuP7nn3+iR48eUCqVWLZsGYqKihAYGIjLly9jzJgx2LFjB5555hmcPHkSDz30kLhviIisgjMsRGR1//jHPxAbG2tw+xdeeAHz5s3DI488gqysLERGRqK4uBhDhw7F0qVLERoainfffRezZs2SsNdEZEnMYSEiIiLZ4wwLERERyR4DFiIiIpI9BixEREQkewxYiIiISPYYsBAREZHsMWAhIiIi2WPAQkRERLLHgIWIiIhkjwELERERyR4DFiIiIpI9BixEREQke/8PRmB9kFKJ0JkAAAAASUVORK5CYII="
"case.blade_loading(pressurevar=\"p\", densityvar=\"rho\", velvar=\"U\")"
"metadata": {
"collapsed": false,
"pycharm": {
"name": "#%%\n"
}
}
},
"name": "stdout",
"output_type": "stream",
"text": [
"2.6096421199420345\n"
"avdr =case.compute_avdr_inout_massave(densityvar=\"rho\", velvar=\"U\")\n",
],
"metadata": {
"collapsed": false,
"pycharm": {
"name": "#%%\n"
}
}
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.6"
}
},
"nbformat": 4,
"nbformat_minor": 0
}