{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Inflow analysis" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Load profile data" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import pathlib\n", "from cfdmod.analysis.inflow.profile import InflowData\n", "\n", "data_folder = pathlib.Path(\"./fixtures/tests/inflow/pitot_inlet\")\n", "profile = InflowData.from_files(\n", " hist_series_path=data_folder / \"hist_series.csv\", points_path=data_folder / \"points.csv\"\n", ")\n", "\n", "# Or you can manually instantiate a profile object\n", "# data = pd.read_csv(data_folder / \"hist_series.csv\")\n", "# points = pd.read_csv(data_folder / \"points.csv\")\n", "# profile = InflowData(data=data, points=points)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Mean velocity data" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjUAAAGdCAYAAADqsoKGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAABO2ElEQVR4nO3de3xU9Z3H/1cmySQQciEJ5EYgKJFbAlEuIVQFJTUWULNrK1AVylJt+ysIpK0Sf1xE7UbborQLW0q3q/3tyuKytWhTio3Ba0m5BSoBgsidkCuRTEjIbeb8/pgwMGaCDIRMMnk/H4/z0Jz5njOfyTDMm+/5nu/XxzAMAxEREZFuzuTpAkREREQ6gkKNiIiIeAWFGhEREfEKCjUiIiLiFRRqRERExCso1IiIiIhXUKgRERERr6BQIyIiIl7Bz9MFdBabzcbZs2cJDg7Gx8fH0+WIiIjINTAMg9raWmJjYzGZrt4X02NCzdmzZ4mPj/d0GSIiInIdTp8+zYABA67apseEmuDgYMD+SwkJCfFwNSIiInItLBYL8fHxju/xq+kxoebSJaeQkBCFGhERkW7mWoaOaKCwiIiIeAWFGhEREfEKCjUiIiLiFXrMmJprYRgGLS0tWK1WT5fiEb6+vvj5+emWdxER6ZYUalo1NTVRWlpKfX29p0vxqN69exMTE4PZbPZ0KSIiIm5RqME+Md/x48fx9fUlNjYWs9nc43orDMOgqamJyspKjh8/TmJi4ldOciQiItKVKNRg76Wx2WzEx8fTu3dvT5fjMb169cLf35+TJ0/S1NREYGCgp0sSERG5Zvqn+BXUM6HfgYiIdF/6BhMRERGvoFAjIiIiXkGhRkRERLyCQo2IiIh4Bd39JCIiItes2Wrj3IUmqi40UnmhkaraRqpaf06IDOLxCYM8VptCTTsMw+Bis2dmFu7l73vN8+QkJCSwaNEiFi1a5NiXkpJCZmYmkydP5r777iM/P5+77roLgJ/97Gf84he/YP/+/URFRd2M8kVEpJtpbLE6gkrVhUaqapuovNBIZW3j5X2tj5+vb273PHclRirUdEUXm62MWP6uR5774PMZ9Dbf+FszefJkFi1axOOPP84//vEPjh07xrJly9i0aZMCjYiIl2totl4RSloDyxU/VzoCTCOWhha3zu1r8iEiyExknwAigwOI7GOmX58AhkYH36RXc20Uarzciy++SF5eHk8++SRFRUXMmTOHBx980NNliYjIdahvanH0olzZq1J1xc+VrZeDLjS6F1T8fX2ICAogMrg1rDg2M/2CnX/u29uMydT1Zt5XqGlHL39fDj6f4bHn7ihms5k33niDUaNGMWjQIF599dUOO7eIiNy4usYWKmuvHJ/SSKWLnpWqC43UN7k3LMLsayKyj7m1N8VVQAmgX2uICe3l3+2XCFKoaYePj0+HXAK62UwmE4ZhOO1rbna+3rl9+3YAqqurqa6uJigoqNPqExGRyywNzRSV1LD/TA37S+zbyXPuLaQc6G9yGUqu7Em5FGJCAv26fVBxR9f/1par6tevH6WlpY6fLRYLx48fd/x89OhRFi9ezG9/+1vefPNN5syZw3vvvaflEEREbrLahmaKSiwUldTwaUkNRSU1HK+qc9m2t9n3ciBxjFMJoN+Xfw4OIMh87TeT9DTXFWrWrl3Lz3/+c8rKyhg9ejT/9m//xvjx49ttv2nTJpYtW8aJEydITEzk5ZdfZurUqS7bfv/73+c3v/kNr776qtMdPdXV1SxYsIA//elPmEwmHn74YX75y1/Sp0+f63kJXuPee+/l9ddf54EHHiAsLIzly5fj62u/fGW1WnnsscfIyMhg7ty53H///SQnJ7Nq1Sp+8pOfeLhyERHvcaGxhQMll3tf9pfUcKzSdYCJC+vFqAGhJA8IJTkulKTYUPoGmTu5Yu/kdqh58803ycrKYt26daSmprJ69WoyMjI4fPgw/fv3b9N++/btzJo1i5ycHKZPn86GDRvIzMyksLCQpKQkp7Z//OMf+fvf/05sbGyb8zz66KOUlpaSl5dHc3Mzc+fO5cknn2TDhg3uvgSvkp2dzfHjx5k+fTqhoaG88MILjp6an/70p5w8eZLc3FwAYmJiWL9+PbNmzeK+++5j9OjRnixdRKRbqmts4WCpxXEJ6dMz5zlWVceXRgIA9gCTFBfCqAFhJMXZQ0y4AsxN42N8eUDGV0hNTWXcuHGsWbMGAJvNRnx8PAsWLGDJkiVt2s+YMYO6ujrHFyvAhAkTSElJYd26dY59JSUlpKam8u677zJt2jSnuVcOHTrEiBEj2LVrF2PHjgVg69atTJ06lTNnzrgMQV9msVgIDQ2lpqaGkJAQp8caGho4fvw4gwcPJjAw0J1fh9fR70JE5LKLTVYOltbw6aUxMGdqOFp5AZuLb86Y0ECSW4PLpV6YiD4BnV+0l7na9/eXudVT09TUxJ49e8jOznbsM5lMpKenU1BQ4PKYgoICsrKynPZlZGSwefNmx882m43HH3+cn/zkJ4wcOdLlOcLCwhyBBiA9PR2TycSOHTv4p3/6J3dehoiISBsNzVanHpj9Z2o4UlHrMsBEhwSSFBdqv4wUF0pSXCj9ghVgPM2tUFNVVYXVam0zcVtUVBTFxcUujykrK3PZvqyszPHzyy+/jJ+fH0899VS75/jypS0/Pz/Cw8OdznOlxsZGGhsbHT9bLJb2X5iIiPQoDc1Wistq2X/mvKMX5kjFBawuEky/4ABGtfa+jBpgDzD9g9WT3RV5/O6nPXv28Mtf/pLCwsIOHc2dk5PDypUrO+x8IiLSPTW2WCkurXX0vuwvqeGz8lpaXASYyD4BjuByKchEhSjAdBduhZrIyEh8fX0pLy932l9eXk50dLTLY6Kjo6/a/uOPP6aiooKBAwc6HrdarfzoRz9i9erVnDhxgujoaCoqKpzO0dLSQnV1dbvPm52d7XTZy2KxEB8ff+0vVkREup2mFhuHy2r5tOS8/VbqM/YA02xtG2Aigsz23pfWy0ejBoQRFRKg26W7MbdCjdlsZsyYMeTn55OZmQnYx8Pk5+czf/58l8ekpaWRn5/vdHt2Xl4eaWlpADz++OOkp6c7HZORkcHjjz/O3LlzHec4f/48e/bsYcyYMQBs27YNm81Gamqqy+cNCAggIEDXN0VEvFVTi43Pymsv30Z9pobiMovLABMeZG4ziDcmNFABxsu4ffkpKyuLOXPmMHbsWMaPH8/q1aupq6tzBJDZs2cTFxdHTk4OAAsXLmTSpEmsWrWKadOmsXHjRnbv3s369esBiIiIICIiwuk5/P39iY6OZujQoQAMHz6c+++/nyeeeIJ169bR3NzM/PnzmTlz5jXd+SQiIt1bs9UeYC71vhSV1HCotJYmq61N27De/o4Ac+lSUlxYLwWYHsDtUDNjxgwqKytZvnw5ZWVlpKSksHXrVsdg4FOnTjnNVjtx4kQ2bNjA0qVLefbZZ0lMTGTz5s1t5qj5Km+88Qbz589nypQpjsn3fvWrX7lbvoiIdHEtVhtHKi44el8+LanhUKmFppa2ASa0l79T70tyXCgD+irA9FRuz1PTXWmemmuj34WIdKYWq42jlXV8eua8YzmBQ6UWGprbBpjgQL82AWZgeG8FGC930+apERERuV5Wm8HRygtOizkeOFvjOsAE+Nln4L0iwAyKUICRq1OoERGRDme1GRyvutC6jIB9DMyBsxbqm6xt2gaZfR1LCNjnggljUHhvTCYFGHGPQo2IiNwQm83g+Lk6p5l4D5ytoc5FgOlt9iUp9ooemAGhDI4IUoCRDqFQ080lJCQ4rZMFkJKSQmZmJqdOnaKiosJp3a3m5mbH3Wnz5s3zQMUi0p3ZbAYnztU5TWR34KyFC40tbdr28vdlZGyII8CMGhDK4Mg++CrAyE2iUNMew4Dmes88t39v6IDrxt/97ne5++67KS0tJSYmBoDc3Fzq6+uZMWPGDZ9fRLybYRicPFfvNA9MUUkNtS4CTKC/iREx9tWoL/XA3NpPAUY6l0JNe5rr4V89NAfOs2fBHHTDp5k4cSJDhw7lv/7rv3j66acBeO211/jWt75Fnz59bvj8IuI9DMPgdPVFPi057xRgLA1tA0yAn4kRsSFOk9kN6dcHP1+TizOLdB6FGi/33e9+l/Xr1/P0009TXl7OX/7yF7Zt2+bpskTEgwzD4MwXF50G8e4vqaHmYnObtmY/E8NjQuzrIF0KMP374K8AI12QQk17/Hvbe0w89dzXyGQy8eWphpqbL//FNHv2bJYsWUJBQQHbt29n8ODB3HXXXR1Wqoh0bYZhUHL+omMm3v0l9hDzRb2LAONrYlhMsFMPzG1RwQow0m0o1LTHx6dDLgHdbP369aO0tNTxs8Vi4fjx446fIyIiyMzM5LXXXqOgoMCxnIWIeB/DMCitaXCaibeopIbquqY2bf19fRgaHUxyXBijWgfy3hYVjNlPAUa6L4Wabu7ee+/l9ddf54EHHiAsLIzly5fj6+vr1Oa73/0u06dPx2q1MmfOHA9VKiIdyTAMyi2NTjPxFpXUUHWhbYDxM10KMJdvpR4aHUyAn6+LM4t0Xwo13Vx2djbHjx9n+vTphIaG8sILLzj11ACkp6cTExPDyJEjtQCoSDdVYWlwXD66NBam6kJjm3a+Jh9uiwomOS6E5AFhjGoNMIH+CjDi/RRqurmQkBA2btzotO/LvTF1dXV88cUXmpdGpJuoqG1wWo360zM1VNS6DjCJ/fs49cAMjwlRgJEeS6HGi9lsNqqqqli1ahVhYWE8+OCDni5JRL6k6kLj5TEwrSGmzNLQpp3JB4b07+MYA5MUF8qImBB6mRVgRC5RqPFip06dYvDgwQwYMIDXX38dPz+93SKedK41wFzZC3O2pm2A8fGBIf2ce2BGxIbQ26zPsMjV6BPixRISEtrc7i0ineOLuianmXj3l9RQcv5im3Y+PnBLZFBrgLHPxjsyNoSgAP31LOIufWpERG5QTX2zffBuyXlHL8yZL9oGGGgNMJcWc4wLZWRcKH0UYEQ6hD5JIiJuqLnYzIHWW6gv9cKcqna9TlxCRO/W3pcQkuPCGBkXQkigfydXLNJzKNRcQZdq9DsQuZKloZmiL42BOXHOdYAZFNGbpLhQx3ICI+NCCe2lACPSmRRqAH9/+1889fX19OrVy8PVeFZ9vf0v7Eu/E5GeorahmQNnLU4B5lhVncu28eG9Wi8f2cfAJMWFENbb3MkVi8iXKdQAvr6+hIWFUVFRAUDv3r3x8fHxcFWdyzAM6uvrqaioICwsrM2sxCLepK6xhQNnLU6z8R6vqsNVR2VcWC+nu5CS40LpG6QAI9IVKdS0io6OBnAEm54qLCzM8bsQ8Qb1TfYAs/+K2XiPVl5wGWBiQwMvh5fWO5HCFWBEug2FmlY+Pj7ExMTQv39/p1WuexJ/f3/10Ei3drHJysFS58UcP6+4gM1FgIkOsQeYUXGhJLUGmcg+AZ1ftIh0GIWaL/H19dUXu0g30NBs5WDpFT0wZ2o4UlHrMsBEhQRcHgMzIISkuFD6Bwd2ftEiclMp1IhIl9fQbKW4rJb9Z847FnU8UnEBq4sE0y84wN77EhfKqNYemP4hCjAiPYFCjYh0KY0tVopLa51m4v2svJYWFwEmso/ZMXj30hiYqJCAHjfQX0TsFGpExGOaWmwcLqt1mon3s/Jamq1tA0x4kD3AXFrMcdSAUKJDAhVgRMRBoUZEOkVTi43Pymud1kMqLrO4DDB9e/tfcfkojOQBocSGKsCIyNUp1IjITdPUYuO1vx1ny/5SDpXW0mS1tWkT2sv/cu9L61iYAX17KcCIiNsUakTkpth+tIplm4s4Wnl5Vt6QQL/WeWDCHJeSFGBEpKMo1IhIh6qobeBf/3yIzfvOAvbBvIvSb+OuxEgGhve82bpFpPMo1IhIh7DaDP777yf5xbuHqW1swccHHksdxI/vG0pob60lJiI3n0KNiNywfafPs3TzfopKLAAkx4Xy039KYtSAMM8WJiI9ikKNiFy3mvpmfvZuMRt2nsIwIDjQj6czhvLt1EH4mnSZSUQ6l+l6Dlq7di0JCQkEBgaSmprKzp07r9p+06ZNDBs2jMDAQJKTk9myZYvT48899xzDhg0jKCiIvn37kp6ezo4dO5zaJCQk4OPj47S99NJL11O+iNwgwzD4w54z3LvqA97YYQ80/3x7HNt+NJnH0xIUaETEI9wONW+++SZZWVmsWLGCwsJCRo8eTUZGRrurW2/fvp1Zs2Yxb9489u7dS2ZmJpmZmRQVFTna3HbbbaxZs4b9+/fzySefkJCQwH333UdlZaXTuZ5//nlKS0sd24IFC9wtX0Ru0Gfltcz4zd/50aZ/cK6uiSH9+/A/T0zglRkp9AvWgpAi4jk+hmG4WP6tfampqYwbN441a9YAYLPZiI+PZ8GCBSxZsqRN+xkzZlBXV0dubq5j34QJE0hJSWHdunUun8NisRAaGsp7773HlClTAHtPzaJFi1i0aJE75bY5Z01NDSEhIdd1DpGerK6xhV/lH+F3nxynxWbQy9+Xp6YkMu/OwZj9rqvTV0TkK7nz/e3W30RNTU3s2bOH9PT0yycwmUhPT6egoMDlMQUFBU7tATIyMtpt39TUxPr16wkNDWX06NFOj7300ktERERw++238/Of/5yWlpZ2a21sbMRisThtIuI+wzDYWlTG11/5kN98dIwWm8F9I6LIy7qbH0y+VYFGRLoMtwYKV1VVYbVaiYqKctofFRVFcXGxy2PKyspcti8rK3Pal5uby8yZM6mvrycmJoa8vDwiIyMdjz/11FPccccdhIeHs337drKzsyktLeWVV15x+bw5OTmsXLnSnZcnIl9y6lw9K94p4v3D9kvBA/r2YuWDI5kyPOorjhQR6Xxd5u6ne+65h3379lFVVcVvf/tbHnnkEXbs2EH//v0ByMrKcrQdNWoUZrOZ733ve+Tk5BAQ0PY6fnZ2ttMxFouF+Pj4m/9CRLxAY4uV9R8eY837n9PYYsPf14fv3X0rP7xnCL3Mvp4uT0TEJbdCTWRkJL6+vpSXlzvtLy8vJzo62uUx0dHR19Q+KCiIIUOGMGTIECZMmEBiYiK/+93vyM7Odnne1NRUWlpaOHHiBEOHDm3zeEBAgMuwIyJX98mRKpa9XcTxKvvyBhNvjeCFzCRu7dfHw5WJiFydWxfDzWYzY8aMIT8/37HPZrORn59PWlqay2PS0tKc2gPk5eW12/7K8zY2Nrb7+L59+zCZTI6eHBG5MeWWBuZvKOSx3+3geFUd/YID+OXMFN74bqoCjYh0C25ffsrKymLOnDmMHTuW8ePHs3r1aurq6pg7dy4As2fPJi4ujpycHAAWLlzIpEmTWLVqFdOmTWPjxo3s3r2b9evXA1BXV8dPf/pTHnzwQWJiYqiqqmLt2rWUlJTwrW99C7APNt6xYwf33HMPwcHBFBQUsHjxYh577DH69u3bUb8LkR6pxWrj/ys4ySt5n3GhsQWTD8xOSyDrvtsICdTyBiLSfbgdambMmEFlZSXLly+nrKyMlJQUtm7d6hgMfOrUKUymyx1AEydOZMOGDSxdupRnn32WxMRENm/eTFJSEgC+vr4UFxfz+9//nqqqKiIiIhg3bhwff/wxI0eOBOyXkjZu3Mhzzz1HY2MjgwcPZvHixU5jZkTEfXtOfsHSzUUcKrXfHZgSH8aLmUkkxYV6uDIREfe5PU9Nd6V5akQu+6KuiZe3FrNx12kAQnv588z9w5g5Lh6TZgMWkS7Ene/vLnP3k4jcfDabwf/tOUPOXw7xRX0zAN8aM4Al3xhGRB8NrBeR7k2hRqSHOFRqYenmIvac/AKAoVHBvPhPSYxLCPdwZSIiHUOhRsTLXWhsYXXeZ7y2/QRWm0Fvsy+L02/jO19LwN9XswGLiPdQqBHxUoZhsGV/Gc/nHqDcYp8eYWpyNMumjyAmtJeHqxMR6XgKNSJe6ERVHcvfOcBHn9mXNxgY3puVD43knqGa10lEvJdCjYgXaWi28usPjvLrD4/S1GLD7GviB5Nv5QeTbyXQX8sbiIh3U6gR8RIfHK5gxTsHOHmuHoC7EiN5/qEkBkcGebgyEZHOoVAj0s2V1lzkhdyDbNlfBkBUSADLp49kanI0Pj6ac0ZEeg6FGpFuqtlq4/W/neDV9z6jvsmKr8mHuRMTWPT12+gToI+2iPQ8+ptPpBvadaKaZZuLKC6rBWDMoL688FASI2I1W7aI9FwKNSLdyLkLjbz0l2I27TkDQN/e/mR/YzjfHDNAyxuISI+nUCPSDdhsBht3neblrcXUXLQvbzBrfDxPZwyjb5DZw9WJiHQNCjUiXVxRSQ1LNxex7/R5AIbHhPBiZhJjBvX1bGEiIl2MQo1IF1Xb0Myqv37G/1dwApsBfQL8yPr6bcxOG4SfljcQEWlDoUakizEMgz99WsqLuQepqLUvbzB9VAzLpo8gKiTQw9WJiHRdCjUiXcjRygssf7uIv31+DoDBkUE8/9BI7krs5+HKRES6PoUakS7gYpOVte9/zm8+Okqz1SDAz8T8e4bw5KRbCPDT8gYiItdCoUbEw7YVl7P87QOc+eIiAPcM7cfKB5MYGNHbw5WJiHQvCjUiHlJy/iIr3znAXw+WAxATGsiKB0aSMTJKyxuIiFwHhRqRTtbUYuN3nxznV/lHuNhsxc/kw7y7BvPUvYkEaXkDEZHrpr9BRTrR34+dY9nmIo5UXABgfEI4L2QmMTQ62MOViYh0fwo1Ip2gsraRnC2HeGtvCQARQWaenTqcf74jTpeaREQ6iEKNyE1ktRls2HmKn28txtLQgo8PfHv8QH6SMZSw3lreQESkIynUiNwkn545z9LNRXx6pgaApLgQXsxMJiU+zLOFiYh4KYUakQ5Wc7GZX7x7mP/ecRLDgOAAP35y/1AeTR2Er1bSFhG5aRRqRDqIYRhs3lfCT/98iKoLTQBkpsTy7LTh9A/W8gYiIjebQo1IBzhSXsvSzUXsOF4NwK39gnghM4mJt0Z6uDIRkZ5DoUbkBtQ3tfCr/M/5j4+P0WIzCPQ3seDeRJ646xbMflpJW0SkMynUiFwHwzDIO1jOyj8dpOS8fXmD9OFRrHhgBPHhWt5ARMQTFGpE3HS6up7n3jlAfnEFAHFhvXjuwZF8fUSUhysTEenZFGpErlFji5X/+Pg4/7btCA3NNvx9fXjirluYf+8Qepv1URIR8TT9TSxyDf72eRXL3i7iWGUdAGm3RPBC5kiG9NfyBiIiXYVCjchVVFgaePHPh3jnH2cBiOwTwNJpw3koJVbLG4iIdDHXdXvG2rVrSUhIIDAwkNTUVHbu3HnV9ps2bWLYsGEEBgaSnJzMli1bnB5/7rnnGDZsGEFBQfTt25f09HR27Njh1Ka6uppHH32UkJAQwsLCmDdvHhcuXLie8kW+UovVxut/O86UVR/yzj/OYvKBOWmDyP/RJDJv13pNIiJdkduh5s033yQrK4sVK1ZQWFjI6NGjycjIoKKiwmX77du3M2vWLObNm8fevXvJzMwkMzOToqIiR5vbbruNNWvWsH//fj755BMSEhK47777qKysdLR59NFHOXDgAHl5eeTm5vLRRx/x5JNPXsdLFrm6vae+4KG1f+O5Px2ktrGF0QNCefuHd7LyoSRCe/l7ujwREWmHj2EYhjsHpKamMm7cONasWQOAzWYjPj6eBQsWsGTJkjbtZ8yYQV1dHbm5uY59EyZMICUlhXXr1rl8DovFQmhoKO+99x5Tpkzh0KFDjBgxgl27djF27FgAtm7dytSpUzlz5gyxsbFfWfelc9bU1BASEuLOS5Ye4nx9Ey9vPczGXacwDAgJ9OPp+4cxa/xALW8gIuIh7nx/u9VT09TUxJ49e0hPT798ApOJ9PR0CgoKXB5TUFDg1B4gIyOj3fZNTU2sX7+e0NBQRo8e7ThHWFiYI9AApKenYzKZ2lymuqSxsRGLxeK0ibhisxls2n2ae1d9yP/stAeah+8YwLYfT+axCVqvSUSku3BroHBVVRVWq5WoKOf5OKKioiguLnZ5TFlZmcv2ZWVlTvtyc3OZOXMm9fX1xMTEkJeXR2RkpOMc/fv3dy7cz4/w8PA257kkJyeHlStXuvPypAcqLrOwbHMRu058AcBtUX144aEkUm+J8HBlIiLiri4zj/s999zDvn372L59O/fffz+PPPJIu+N0rkV2djY1NTWO7fTp0x1YrXR3dY0t/OuWQ0z71SfsOvEFvfx9yf7GMP781F0KNCIi3ZRbPTWRkZH4+vpSXl7utL+8vJzo6GiXx0RHR19T+6CgIIYMGcKQIUOYMGECiYmJ/O53vyM7O5vo6Og2AaelpYXq6up2nzcgIICAgAB3Xp70AIZhsLWojJV/OkiZpQGA+0dGs+yBEcSF9fJwdSIiciPc6qkxm82MGTOG/Px8xz6bzUZ+fj5paWkuj0lLS3NqD5CXl9du+yvP29jY6DjH+fPn2bNnj+Pxbdu2YbPZSE1NdeclSA92oqqO77y2ix+8UUiZpYH48F689p1xrHt8jAKNiIgXcHvyvaysLObMmcPYsWMZP348q1evpq6ujrlz5wIwe/Zs4uLiyMnJAWDhwoVMmjSJVatWMW3aNDZu3Mju3btZv349AHV1dfz0pz/lwQcfJCYmhqqqKtauXUtJSQnf+ta3ABg+fDj3338/TzzxBOvWraO5uZn58+czc+bMa7rzSXq2hmYr6z48yr9/cJSmFhtmXxPfn3QL/889Qwj09/V0eSIi0kHcDjUzZsygsrKS5cuXU1ZWRkpKClu3bnUMBj516hQm0+UOoIkTJ7JhwwaWLl3Ks88+S2JiIps3byYpKQkAX19fiouL+f3vf09VVRURERGMGzeOjz/+mJEjRzrO88YbbzB//nymTJmCyWTi4Ycf5le/+tWNvn7xch99Vsnyt4s4ca4egDuHRLLyoZHc2q+PhysTEZGO5vY8Nd2V5qnpWcpqGngh9yB/3l8KQP/gAJZNH8H0UTGaDVhEpBtx5/tbaz+JV2mx2nh9+wlezfuMuiYrJh/4zsTBLP56IsGBmg1YRMSbKdSI19h9opqlm4soLqsF4PaBYbyYmcTI2FAPVyYiIp1BoUa6veq6Jl76yyH+d/cZAMJ6+/PM/cOYMTYek2YDFhHpMRRqpNuy2Qz+d/dpXtpazPn6ZgAeGTuAJd8YTniQ2cPViYhIZ1OokW7pwNkalm4uYu+p8wAMiw7mxcwkxiaEe7YwERHxGIUa6VZqG5p5Je8zfr/9BDYDgsy+LP76bXxnYgJ+vl1m1Q8REfEAhRrpFgzDIPfTUl7IPUhFrX2m6WnJMSydPpyYUM0GLCIiCjXSDRyrvMDytw/wyedVACRE9GblQ0lMuq2fhysTEZGuRKFGuqyGZiv//v7nrPvwGE1WG2Y/Ez+cPITvTbpFyxuIiEgbCjXSJb1fXMHyd4o4XX0RgLtv68fzD44kITLIw5WJiEhXpVAjXcrZ8xdZ+acDvHugHIDokEBWPDCC+5OitbyBiIhclUKNdAnNVhv/+clxfpl/hPomK74mH/7lawksTL+NPgH6YyoiIl9N3xbicTuOnWPZ20V8Vn4BgLGD+vLiPyUxLFoLj4qIyLVTqBGPqbrQSM6WYv5QaF/eIDzITPY3hvHwHQO0vIGIiLhNoUY6ndVm8D87T/GzrcVYGlrw8YGZ4wbyzP1DCeut5Q1EROT6KNRIp9p/poalbxfxj9PnARgRE8KL/5TEHQP7erYwERHp9hRqpFPUXGzmlb8e5r/+fhKbAcEBfvzovtt4bMIgLW8gIiIdQqFGbirDMHh731le/PMhqi7Ylzd4cHQsS6cNp39IoIerExERb6JQIzfN5xW1LNt8gIJj5wC4pV8QLzyUxNeGRHq4MhER8UYKNdLhLjZZ+bdtR/jtx8dothoE+JlYcO8Qnrj7FgL8tLyBiIjcHAo10qHyDpbz3DsHKDlvX97g3mH9WfngSOLDe3u4MhER8XYKNdIhTlfXs/JPB3nvkH15g7iwXqx4YARfHxGl5Q1ERKRTKNTIDWlqsfHbj4/xb9uO0NBsw8/kwxN338KCe4fQ26w/XiIi0nn0rSPXbfvRKpZtLuJoZR0AqYPDeTEzicSoYA9XJiIiPZFCjbitoraBf/3zITbvOwtAZB8z/++04WSmxOlSk4iIeIxCjVwzq83gv/9+kl+8e5jaRvvyBo+lDuLHGUMJ7eXv6fJERKSHU6iRa7Lv9HmWbt5PUYkFgFEDQnkxM4lRA8I8W5iIiEgrhRq5qpr6Zn72bjEbdp7CMCA40I+n7x/Gt8cPxFcraYuISBeiUCMuGYbBW4Ul/OuWQ5yrawLgn2+PI3vqcPoFB3i4OhERkbYUaqSNz8prWfrHInaeqAZgSP8+vJiZxIRbIjxcmYiISPsUasShrrGFX+Uf4XefHKfFZtDL35eF6Yn8y9cGY/bTStoiItK1KdQIhmHw7oFynv/TAc7WNABw34golj8wggF9tbyBiIh0Dwo1Pdypc/WseKeI9w9XAjCgby9WPjiSKcOjPFyZiIiIe67rmsLatWtJSEggMDCQ1NRUdu7cedX2mzZtYtiwYQQGBpKcnMyWLVscjzU3N/PMM8+QnJxMUFAQsbGxzJ49m7NnzzqdIyEhAR8fH6ftpZdeup7yBWhssfJv+Uf4+qsf8v7hSvx9fZh/zxDyFk9SoBERkW7J7VDz5ptvkpWVxYoVKygsLGT06NFkZGRQUVHhsv327duZNWsW8+bNY+/evWRmZpKZmUlRUREA9fX1FBYWsmzZMgoLC3nrrbc4fPgwDz74YJtzPf/885SWljq2BQsWuFu+AJ8cqeL+1R+zKu8zGltsfG1IBFsX3c2PM4bSy+zr6fJERESui49hGIY7B6SmpjJu3DjWrFkDgM1mIz4+ngULFrBkyZI27WfMmEFdXR25ubmOfRMmTCAlJYV169a5fI5du3Yxfvx4Tp48ycCBAwF7T82iRYtYtGiRO+U6WCwWQkNDqampISQk5LrO0d2VWxp4IfcguZ+WAtAvOIBl00fwwKgYLW8gIiJdkjvf32711DQ1NbFnzx7S09Mvn8BkIj09nYKCApfHFBQUOLUHyMjIaLc9QE1NDT4+PoSFhTntf+mll4iIiOD222/n5z//OS0tLe2eo7GxEYvF4rT1VC1WG//5yXGmrPqQ3E9LMfnAdyYmkP+jSTw4OlaBRkREvIJbA4WrqqqwWq1ERTmPuYiKiqK4uNjlMWVlZS7bl5WVuWzf0NDAM888w6xZs5wS2VNPPcUdd9xBeHg427dvJzs7m9LSUl555RWX58nJyWHlypXuvDyvVHjqC5b+sYiDpfZQlxIfxouZSSTFhXq4MhERkY7Vpe5+am5u5pFHHsEwDH796187PZaVleX4/1GjRmE2m/ne975HTk4OAQFtZ7jNzs52OsZisRAfH3/ziu9ivqhr4mfvFvM/O08DENrLn2fuH8bMcfGYtLyBiIh4IbdCTWRkJL6+vpSXlzvtLy8vJzo62uUx0dHR19T+UqA5efIk27Zt+8rrZqmpqbS0tHDixAmGDh3a5vGAgACXYcfb2WwG/7fnDDl/OcQX9c0AfGvMAJZ8YxgRfXre70NERHoOt8bUmM1mxowZQ35+vmOfzWYjPz+ftLQ0l8ekpaU5tQfIy8tzan8p0Bw5coT33nuPiIivno5/3759mEwm+vfv785L8GqHSi088psCnv7Dp3xR38zQqGA2fT+Nn39rtAKNiIh4PbcvP2VlZTFnzhzGjh3L+PHjWb16NXV1dcydOxeA2bNnExcXR05ODgALFy5k0qRJrFq1imnTprFx40Z2797N+vXrAXug+eY3v0lhYSG5ublYrVbHeJvw8HDMZjMFBQXs2LGDe+65h+DgYAoKCli8eDGPPfYYffv27ajfRbd1obGF1Xmf8dr2E1htBr3NvixOv43vfC0Bf18tbyAiIj2D26FmxowZVFZWsnz5csrKykhJSWHr1q2OwcCnTp3CZLr8RTpx4kQ2bNjA0qVLefbZZ0lMTGTz5s0kJSUBUFJSwjvvvANASkqK03O9//77TJ48mYCAADZu3Mhzzz1HY2MjgwcPZvHixU5jZnoiwzDYsr+M53MPUG5pBGBqcjTLpo8gJrSXh6sTERHpXG7PU9Ndeds8NSeq6lj+zgE++sy+vMGgiN6sfHAkk4fqcpyIiHgPd76/u9TdT/LVGpqt/PqDo/z6w6M0tdgw+5r4weRb+cHkWwn012zAIiLScynUdCMfHK5gxTsHOHmuHoC7EiN5/qEkBkcGebgyERERz1Oo6QZKay7yQu5Btuy3D6COCglg+fSRTE2O1mzAIiIirRRqurBmq43fbz/Bq3mfUddkxdfkw9yJCSz6+m30CdBbJyIiciV9M3ZRu09Us3RzEcVltQCMGdSXFx5KYkRs9x/kLCIicjMo1HQx5y408tJfitm05wwAfXv7k/2N4XxzzAAtbyAiInIVCjVdhM1msHHXaV7eWkzNRfvyBrPGx/N0xjD6Bpk9XJ2IiEjXp1DTBRSV1LB0cxH7Tp8HYHhMCC9mJjFmkGZLFhERuVYKNR5U29DMK3mf8fvtJ7AZEGT2Jeu+ocxJG4SfljcQERFxi0KNBxiGwZ8+LeXF3INU1NqXN5g+Koal00YQHRro4epERES6J4WaTna08gLL3y7ib5+fA2BwZBDPPzSSuxL7ebgyERGR7k2hppM0NFtZ+/7n/ObDYzRZbQT4mZh/zxCenHQLAX5a3kBERORGKdR0gm3F5ax45wCnqy8CMHloP55/MImBEb09XJmIiIj3UKi5iUrOX2TlOwf468FyAGJCA1nxwEgyRkZpeQMREZEOplBzEzRbbfzuk+P88r0jXGy24mfyYd6dg3lqSiJBWt5ARETkptA3bAfbcewcSzcXcaTiAgDjE8J5ITOJodHBHq5MRETEuynUdJCqC43865ZDvFVYAkBEkJnsqcN5+I44XWoSERHpBAo1N8hqM9iw8xQ/31qMpaEFHx/49viB/CRjKGG9tbyBiIhIZ1GouUGv/e04L/75EABJcSG8mJlMSnyYZ4sSERHpgRRqbtDM8QP5392neTR1EI9NGISvVtIWERHxCIWaG9QnwI+tC+/GpDAjIiLiUVo1sQMo0IiIiHieQo2IiIh4BYUaERER8QoKNSIiIuIVFGpERETEKyjUiIiIiFdQqBERERGvoFAjIiIiXkGhRkRERLyCQo2IiIh4BYUaERER8QoKNSIiIuIVrivUrF27loSEBAIDA0lNTWXnzp1Xbb9p0yaGDRtGYGAgycnJbNmyxfFYc3MzzzzzDMnJyQQFBREbG8vs2bM5e/as0zmqq6t59NFHCQkJISwsjHnz5nHhwoXrKV9ERES8kNuh5s033yQrK4sVK1ZQWFjI6NGjycjIoKKiwmX77du3M2vWLObNm8fevXvJzMwkMzOToqIiAOrr6yksLGTZsmUUFhby1ltvcfjwYR588EGn8zz66KMcOHCAvLw8cnNz+eijj3jyySev4yWLiIiIN/IxDMNw54DU1FTGjRvHmjVrALDZbMTHx7NgwQKWLFnSpv2MGTOoq6sjNzfXsW/ChAmkpKSwbt06l8+xa9cuxo8fz8mTJxk4cCCHDh1ixIgR7Nq1i7FjxwKwdetWpk6dypkzZ4iNjf3Kui0WC6GhodTU1BASEuLOSxYREREPcef7262emqamJvbs2UN6evrlE5hMpKenU1BQ4PKYgoICp/YAGRkZ7bYHqKmpwcfHh7CwMMc5wsLCHIEGID09HZPJxI4dO1yeo7GxEYvF4rSJiIiI93Ir1FRVVWG1WomKinLaHxUVRVlZmctjysrK3Grf0NDAM888w6xZsxyJrKysjP79+zu18/PzIzw8vN3z5OTkEBoa6tji4+Ov6TWKiIhI99Sl7n5qbm7mkUcewTAMfv3rX9/QubKzs6mpqXFsp0+f7qAqRUREpCvyc6dxZGQkvr6+lJeXO+0vLy8nOjra5THR0dHX1P5SoDl58iTbtm1zum4WHR3dZiByS0sL1dXV7T5vQEAAAQEB1/zaREREpHtzq6fGbDYzZswY8vPzHftsNhv5+fmkpaW5PCYtLc2pPUBeXp5T+0uB5siRI7z33ntERES0Ocf58+fZs2ePY9+2bduw2Wykpqa68xJERETES7nVUwOQlZXFnDlzGDt2LOPHj2f16tXU1dUxd+5cAGbPnk1cXBw5OTkALFy4kEmTJrFq1SqmTZvGxo0b2b17N+vXrwfsgeab3/wmhYWF5ObmYrVaHeNkwsPDMZvNDB8+nPvvv58nnniCdevW0dzczPz585k5c+Y13fkkIiIi3s/tUDNjxgwqKytZvnw5ZWVlpKSksHXrVsdg4FOnTmEyXe4AmjhxIhs2bGDp0qU8++yzJCYmsnnzZpKSkgAoKSnhnXfeASAlJcXpud5//30mT54MwBtvvMH8+fOZMmUKJpOJhx9+mF/96lfX85pFRETEC7k9T013pXlqREREup+bNk+NiIiISFelUCMiIiJeQaFGREREvIJCjYiIiHgFhRoRERHxCgo1IiIi4hUUakRERMQrKNSIiIiIV1CoEREREa+gUCMiIiJeQaFGREREvIJCjYiIiHgFhRoRERHxCgo1IiIi4hUUakRERMQrKNSIiIiIV1CoEREREa+gUCMiIiJeQaFGREREvIJCjYiIiHgFhRoRERHxCgo1IiIi4hUUakRERMQrKNSIiIiIV1CoEREREa+gUCMiIiJeQaFGREREvIJCjYiIiHgFhRoRERHxCgo1IiIi4hUUakRERMQrKNSIiIiIV1CoEREREa9wXaFm7dq1JCQkEBgYSGpqKjt37rxq+02bNjFs2DACAwNJTk5my5YtTo+/9dZb3HfffURERODj48O+ffvanGPy5Mn4+Pg4bd///vevp3wRERHxQm6HmjfffJOsrCxWrFhBYWEho0ePJiMjg4qKCpftt2/fzqxZs5g3bx579+4lMzOTzMxMioqKHG3q6uq48847efnll6/63E888QSlpaWO7Wc/+5m75YuIiIiX8jEMw3DngNTUVMaNG8eaNWsAsNlsxMfHs2DBApYsWdKm/YwZM6irqyM3N9exb8KECaSkpLBu3TqntidOnGDw4MHs3buXlJQUp8cmT55MSkoKq1evdqdcB4vFQmhoKDU1NYSEhFzXOURERKRzufP97VZPTVNTE3v27CE9Pf3yCUwm0tPTKSgocHlMQUGBU3uAjIyMdttfzRtvvEFkZCRJSUlkZ2dTX1/fbtvGxkYsFovTJiIiIt7Lz53GVVVVWK1WoqKinPZHRUVRXFzs8piysjKX7cvKytwq9Nvf/jaDBg0iNjaWTz/9lGeeeYbDhw/z1ltvuWyfk5PDypUr3XoOERER6b7cCjWe9OSTTzr+Pzk5mZiYGKZMmcLRo0e59dZb27TPzs4mKyvL8bPFYiE+Pr5TahUREZHO51aoiYyMxNfXl/Lycqf95eXlREdHuzwmOjrarfbXKjU1FYDPP//cZagJCAggICDghp5DREREug+3xtSYzWbGjBlDfn6+Y5/NZiM/P5+0tDSXx6SlpTm1B8jLy2u3/bW6dNt3TEzMDZ1HREREvIPbl5+ysrKYM2cOY8eOZfz48axevZq6ujrmzp0LwOzZs4mLiyMnJweAhQsXMmnSJFatWsW0adPYuHEju3fvZv369Y5zVldXc+rUKc6ePQvA4cOHAXsvT3R0NEePHmXDhg1MnTqViIgIPv30UxYvXszdd9/NqFGjbviXICIiIt2f26FmxowZVFZWsnz5csrKykhJSWHr1q2OwcCnTp3CZLrcATRx4kQ2bNjA0qVLefbZZ0lMTGTz5s0kJSU52rzzzjuOUAQwc+ZMAFasWMFzzz2H2WzmvffecwSo+Ph4Hn74YZYuXXrdL1xERES8i9vz1HRXmqdGRESk+7lp89SIiIiIdFUKNSIiIuIVFGpERETEKyjUiIiIiFdQqBERERGvoFAjIiIiXkGhRkRERLyCQo2IiIh4BYUaERER8QoKNSIiIuIVFGpERETEKyjUiIiIiFdQqBERERGvoFAjIiIiXkGhRkRERLyCQo2IiIh4BYUaERER8QoKNSIiIuIVFGpERETEKyjUiIiIiFdQqBERERGvoFAjIiIiXkGhRkRERLyCQo2IiIh4BYUaERER8QoKNSIiIuIVFGpERETEKyjUiIiIiFdQqBERERGvoFAjIiIiXkGhRkRERLyCQo2IiIh4BYUaERER8QrXFWrWrl1LQkICgYGBpKamsnPnzqu237RpE8OGDSMwMJDk5GS2bNni9Phbb73FfffdR0REBD4+Puzbt6/NORoaGvjhD39IREQEffr04eGHH6a8vPx6yhcREREv5HaoefPNN8nKymLFihUUFhYyevRoMjIyqKiocNl++/btzJo1i3nz5rF3714yMzPJzMykqKjI0aauro4777yTl19+ud3nXbx4MX/605/YtGkTH374IWfPnuWf//mf3S1fREREvJSPYRiGOwekpqYybtw41qxZA4DNZiM+Pp4FCxawZMmSNu1nzJhBXV0dubm5jn0TJkwgJSWFdevWObU9ceIEgwcPZu/evaSkpDj219TU0K9fPzZs2MA3v/lNAIqLixk+fDgFBQVMmDDhK+u2WCyEhoZSU1NDSEiIOy9ZREREPMSd72+3emqamprYs2cP6enpl09gMpGenk5BQYHLYwoKCpzaA2RkZLTb3pU9e/bQ3NzsdJ5hw4YxcODAds/T2NiIxWJx2kRERMR7uRVqqqqqsFqtREVFOe2PioqirKzM5TFlZWVutW/vHGazmbCwsGs+T05ODqGhoY4tPj7+mp9PREREuh+vvfspOzubmpoax3b69GlPlyTS81hboPo4VH4GNpunqxERL+fnTuPIyEh8fX3b3HVUXl5OdHS0y2Oio6Pdat/eOZqamjh//rxTb83VzhMQEEBAQMA1P4eIXKeWJqg5DeeOQvUx5+38SbC12NsFhEDcHTBgPMSPh7gx0Dvcs7WLiFdxK9SYzWbGjBlDfn4+mZmZgH2gcH5+PvPnz3d5TFpaGvn5+SxatMixLy8vj7S0tGt+3jFjxuDv709+fj4PP/wwAIcPH+bUqVNunUdErlNzgz2gfDm0VB+D86fBsLZ/rF8g+Jig0QLHPrBvl0QkwoBxED/O/t/+I8Dke7NfjYh4KbdCDUBWVhZz5sxh7NixjB8/ntWrV1NXV8fcuXMBmD17NnFxceTk5ACwcOFCJk2axKpVq5g2bRobN25k9+7drF+/3nHO6upqTp06xdmzZwF7YAF7D010dDShoaHMmzePrKwswsPDCQkJYcGCBaSlpV3TnU8icg2a6uCLEy6Cy3GoOQNc5UZJ/yAIvwXCB7f+94otOAYMG1QegtM74cxuOLMTzn0O547Yt39ssJ/H3Adib7f35AxoDTpBkZ3x6kXEC7gdambMmEFlZSXLly+nrKyMlJQUtm7d6hgMfOrUKUymy0N1Jk6cyIYNG1i6dCnPPvssiYmJbN68maSkJEebd955xxGKAGbOnAnAihUreO655wB49dVXMZlMPPzwwzQ2NpKRkcG///u/X9eLFumxGizwxfG2oaX6GNSWXv3YgJC2geXS1qc/+Phc5WATRCfbt3Hz7Lvqqy8HnDO74MweaKqFEx/bt0vCb7kccAaMg6gk8HX7ry4R6QHcnqemu9I8NdJjXPzCOaxcudVVXv3YXn1bg8qtbYNL7/CvCC43yGaFysOXQ87pXVB1uG07/9723pxLISd+vD1UiYhXcuf7W6FGpLsxDKg/53p8S/Uxe6i5mqB+LnpbBkPfwV1v4O7F81Cy296jc+nSVWNN23ZhAy8PQB4wFqKSwc/c6eWKSMdTqHFBoUa6FcOAC+XtBJfj9kG3VxMc43qMS9/BENiN//zbbPYxOKcvXbLaBRWHaDPexy8QYlIuD0AeMB5CYjxRsYjcIIUaFxRqpMux2aD2rOvQUn0MmuuvcrAPhA5wPTC3bwKYgzrrVXhegwVK9jiPz3HVWxUywDnkxIwCP037INLVKdS4oFAjHmFtAcuZtoHl0v9bG9s/1scEofEQ4WJ8S9gg8A/svNfRnRiGfc6cM7suh5zyA/Y7sK7ka4aY0faAM2Cs/dJVSNzNHTckIm5TqHFBoUZuGmsznD/l+lLRFyfB1tz+sSY/e8+KqzuKQuM1LqSjNF6As3vtIed062Wr+qq27YJjnAcgx4wG/16dX6+IOCjUuKBQIzfkRiaf8w340mWiK/4/ZIBuT/YEw7Df2u4YgLwLyva3fR9N/vbb0C+FnAFj7b1k6s0R6TQKNS4o1MhXaqp3MYfLtU4+1/sqk8/FwhVzN0kX1VQPpfucByFfKG/bLqj/5YAzYDzEpvSsMUwinUyhxgWFGgFubPI5czBEtDf5XJT+9e5tDMO+ptWVsyCXftr2cqKPL0QnXR6APGCs/c+E/jyIdAiFGhcUanqQDpl8zsXWO0JfVD1dcwOU/sN5gsDas23b9Y5wngU5bgwE9On8ekW8gEKNCwo1XqQnTT4nXV9NSWvIaR2fU7oPrE3ObXxM9sU6rxyEHDFEIVnkGijUuKBQ083c6ORzfaKumOp/8OXborv75HPS9bU02gcdXzk2p+Z023aBYc4DkOPGQGBop5cr0tUp1LigUNMF3dDkc9jvHHI1MDd8sAZuStdiKb0ccM7sst9e3tLwpUY+0G+Y8wSBkbdpkLn0eAo1LijUeIjNav9X6o1MPudqfEvfQZo/RLova7O9N+fSAOTTO+1TBnxZQCgMGHPFIOQx9nFfIj2IQo0LCjU30Y1OPhc2yHVwCRuoyeek57hQcbkn5/QuOFvourcy8jbnQcj9h4PJt/PrFekkCjUuKNTcoBuafM5sH8viavK50HhNPifiirUFKg5cDjlndkH10bbtzH0g7o7Lq5THjYWgiM6vV+QmUahxQaHmGtzI5HN+vdqffC4kVv+SFOkIdeeuGJuzE0oKoelC23bht7YOQr7UmzNS/3iQbkuhxgWFmlaNta7nb7nWyedcDsy9BYKjdXuqSGezWaHikPMg5KrP2rbz7w2xdzgPQu7Tr/PrFbkOCjUu9KhQ4zT53JcCTF3F1Y8NDL3iVugvbUGRCi4iXV19tb0H59IA5JI9rqdACBvUejt5a9CJTgZf/86vV+QrKNS44FWhxjDsf3E5wspR9yaf6x3Zzqy5mnxOxOvYbFB1uHVsTuskgZXFtLmc7BcIsbc7TxAYHO2RkkWupFDjQrcLNYZhvxui3cnnaq5+fJ9o1wNzwwdrgi+Rnq6hxt6Dc/qKy1YN59u2C413DjnRyeAX0OnlSs+mUONClww1Npt9HEu7k8/VXf349iaf65ugdWZE5NrZbPYe3ytnQa44CIbNuZ1vAMSMdl6lPDTOMzVLj6FQ44LHQo3Nar9zyFVo+eK4i1lFr6DJ50TEUxprL4/NubSu1cXqtu2CY50HIMeMBv/Azq9XvJZCjQs3NdQ4Jp9zcVfRFyc0+ZyIdH+GYf87zTFB4E4oP9B2jiqTP8SMcr5sFRqvmwzkuinUuHDTQs3Bd2DTdzT5nIj0PE119nWsHBME7oS6yrbt+kQ5h5yYFDD37vRypXty5/tb36Y3KjjGHmg0+ZyI9DTmIEi4076BvTfn/MnLl6vO7IKyT+FCORTn2jew91BHJTmvUt53sHpz5Iapp+ZGtTTab6/W5HMiIm01X4TSf7SGnJ32Hp0LZW3bBfVr7c1pHYAce7tueBBAl59c6pJ3P4mI9DSGYb954spZkEv/AdYm53Y+JogaecUK5eMg4lb947EHUqhxQaFGRKSLammE0k9b77RqHZ9jOdO2Xa++l0NO/Dj70g+B+vvc2ynUuKBQIyLSjVjOXnGn1S77gGRr45ca+UD/EfZLVpeWfIhIBJPJIyXLzaFQ44JCjYhIN9bSBOX7nQchnz/Ztl1gKMSNvTwAOW4s9Arr9HKl4yjUuKBQIyLiZWrLoeSKkFNSCC0X27aLHOo8QWC/obojtRtRqHFBoUZExMtZm+0TAl45CLn6WNt2ASEQd8cVg5DHajHfLkyhxgWFGhGRHqiuyn7J6sxOe49OSaHrdfUihlwOOPHjod9wTYzaRbjz/X1do6nWrl1LQkICgYGBpKamsnPnzqu237RpE8OGDSMwMJDk5GS2bNni9LhhGCxfvpyYmBh69epFeno6R44ccWqTkJCAj4+P0/bSSy9dT/kiItJTBEXC0PthynL4Ti5kn4bvfwLTX4XR37YPLAY49zn8YwP8OQvW3QkvDYTXp8N7K+HwX+zhSLo8t3tq3nzzTWbPns26detITU1l9erVbNq0icOHD9O/f/827bdv387dd99NTk4O06dPZ8OGDbz88ssUFhaSlJQEwMsvv0xOTg6///3vGTx4MMuWLWP//v0cPHiQwED7wmgJCQnMmzePJ554wnHu4OBggoKCrqlu9dSIiIhL9dVQsueKsTl7oNHStl3fwc6zIEclga9/59fbw9zUy0+pqamMGzeONWvWAGCz2YiPj2fBggUsWbKkTfsZM2ZQV1dHbm6uY9+ECRNISUlh3bp1GIZBbGwsP/rRj/jxj38MQE1NDVFRUbz++uvMnDkTsIeaRYsWsWjRInfKdVCoERGRa2KzQtVnl0POmV1QWdy2nV8v+8zHVw5CDo7q/Hq93E1b+6mpqYk9e/aQnZ3t2GcymUhPT6egoMDlMQUFBWRlZTnty8jIYPPmzQAcP36csrIy0tPTHY+HhoaSmppKQUGBI9QAvPTSS7zwwgsMHDiQb3/72yxevBg/P9cvobGxkcbGy3MaWCwuUreIiMiXmXyh/3D7NmaOfd/F8/YenCsHITfUwKnt9u2S0IHOISc6GfzMHnkZPZFboaaqqgqr1UpUlHMSjYqKorjYRYoFysrKXLYvKytzPH5pX3ttAJ566inuuOMOwsPD2b59O9nZ2ZSWlvLKK6+4fN6cnBxWrlzpzssTERFxrVcYDJli3wBsNvs4nCtnQa44CDWn7FvRH+ztfAMgNsV5lfKQWE+9Cq/XbYZ2X9nbM2rUKMxmM9/73vfIyckhICCgTfvs7GynYywWC/Hx8Z1Sq4iIeDmTCfrdZt9uf8y+r7H2cm/O6dbenIvVcHqHfbskJM455ESPAv9Az7wOL+NWqImMjMTX15fy8nKn/eXl5URHR7s8Jjo6+qrtL/23vLycmJgYpzYpKSnt1pKamkpLSwsnTpxg6NChbR4PCAhwGXZERERuioBguGWyfQP74p3Vx1pDTusq5eUHwFICB0vg4GZ7O18zDPoaDJsGQ78BoQM89AK6P7du6TabzYwZM4b8/HzHPpvNRn5+PmlpaS6PSUtLc2oPkJeX52g/ePBgoqOjndpYLBZ27NjR7jkB9u3bh8lkcnnHlYiIiMf5+NhXFh89E6a/Yr+VfMlpmJMLU1bA0KnQO9K+Qvmx92HLj+HVkbDuLvjgJfsinz1jKrkO4/blp6ysLObMmcPYsWMZP348q1evpq6ujrlz5wIwe/Zs4uLiyMnJAWDhwoVMmjSJVatWMW3aNDZu3Mju3btZv349AD4+PixatIgXX3yRxMRExy3dsbGxZGZmAvbBxjt27OCee+4hODiYgoICFi9ezGOPPUbfvn076FchIiJykwX0gcF32Tewh5Zzn8PhLfb5cE79Hco+tW8f5EBovL33ZuhUe2+OBh1flduhZsaMGVRWVrJ8+XLKyspISUlh69atjoG+p06dwnTFCqkTJ05kw4YNLF26lGeffZbExEQ2b97smKMG4Omnn6auro4nn3yS8+fPc+edd7J161bHHDUBAQFs3LiR5557jsbGRgYPHszixYvb3FUlIiLSrfj4QGQiRC6Ery20T/L32VZ7wPk8H2pOw8719i0gFBLT7QFnSLoW6nRByySIiIh0Rc0X4dgHrb04W6Gu4vJjJj9IuBOGto7DCfPeG2G09pMLCjUiItJt2Wz2FckPb4HiLVB12Pnx6OTLASdmtL0HyEso1LigUCMiIl7j3NHLAef038GwXX4sZEDrOJxvQMJd3X4cjkKNCwo1IiLileqq4MhfofjPcHQbNNdffiwgxD7+ZuhUSPx6txyHo1DjgkKNiIh4veaLcPwje8A5/Je243AGfc0ecIZNhbCBnqvTDQo1LijUiIhIj2KzwdnC1oCzpe2inFHJ9ktUw6ZCTEqXHYejUOOCQo2IiPRo547ae28Ob4FTBc7jcIJjLwechLvAr+vMyK9Q44JCjYiISKu6c/ZxOIf/DJ9vg+a6y4+Zg+0Ldw6b1joOx7OT3CrUuKBQIyIi4kJzg30czuHWcTgXrliv0ccXBk1sXZdqKvQd1OnlKdS4oFAjIiLyFWw2OLvXHnCKt0DlIefH+4+0X6IaOhVib++UcTgKNS4o1IiIiLip+pi996Z4C5za/qVxODGt8+FMs69ldZPG4SjUuKBQIyIicgPqqy/Ph/N5/pfG4fSxj8MZ/iAkf7NDn9ad72+3F7QUERGRHqh3OIyead+aG+DEx5fnw7lQBgffhguVHR5q3KFQIyIiIu7xD7TfGZX4dZj2CpTutV+iirzNo2Up1IiIiMj1M5kgbox983Qpni5AREREpCMo1IiIiIhXUKgRERERr6BQIyIiIl5BoUZERES8gkKNiIiIeAWFGhEREfEKCjUiIiLiFRRqRERExCso1IiIiIhXUKgRERERr6BQIyIiIl5BoUZERES8Qo9ZpdswDAAsFouHKxEREZFrdel7+9L3+NX0mFBTW1sLQHx8vIcrEREREXfV1tYSGhp61TY+xrVEHy9gs9k4e/YswcHB+Pj4dOi5LRYL8fHxnD59mpCQkA49t7hP70fXoveja9H70fXoPbk6wzCora0lNjYWk+nqo2Z6TE+NyWRiwIABN/U5QkJC9AeyC9H70bXo/eha9H50PXpP2vdVPTSXaKCwiIiIeAWFGhEREfEKCjUdICAggBUrVhAQEODpUgS9H12N3o+uRe9H16P3pOP0mIHCIiIi4t3UUyMiIiJeQaFGREREvIJCjYiIiHgFhRoRERHxCgo1X+Gjjz7igQceIDY2Fh8fHzZv3vyVx3zwwQfccccdBAQEMGTIEF5//fWbXmdP4u578sEHH+Dj49NmKysr65yCvVhOTg7jxo0jODiY/v37k5mZyeHDh7/yuE2bNjFs2DACAwNJTk5my5YtnVCt97ue9+P1119v89kIDAzspIq9369//WtGjRrlmFgvLS2Nv/zlL1c9Rp+P66dQ8xXq6uoYPXo0a9euvab2x48fZ9q0adxzzz3s27ePRYsW8d3vfpd33333Jlfac7j7nlxy+PBhSktLHVv//v1vUoU9x4cffsgPf/hD/v73v5OXl0dzczP33XcfdXV17R6zfft2Zs2axbx589i7dy+ZmZlkZmZSVFTUiZV7p+t5P8A+k+2Vn42TJ092UsXeb8CAAbz00kvs2bOH3bt3c++99/LQQw9x4MABl+31+bhBhlwzwPjjH/941TZPP/20MXLkSKd9M2bMMDIyMm5iZT3Xtbwn77//vgEYX3zxRafU1JNVVFQYgPHhhx+22+aRRx4xpk2b5rQvNTXV+N73vnezy+txruX9eO2114zQ0NDOK0qMvn37Gv/xH//h8jF9Pm6Memo6WEFBAenp6U77MjIyKCgo8FBFcklKSgoxMTF8/etf529/+5uny/FKNTU1AISHh7fbRp+RznMt7wfAhQsXGDRoEPHx8VftRZAbY7Va2bhxI3V1daSlpblso8/HjVGo6WBlZWVERUU57YuKisJisXDx4kUPVdWzxcTEsG7dOv7whz/whz/8gfj4eCZPnkxhYaGnS/MqNpuNRYsW8bWvfY2kpKR227X3GdEYp451re/H0KFD+c///E/efvtt/vu//xubzcbEiRM5c+ZMJ1br3fbv30+fPn0ICAjg+9//Pn/84x8ZMWKEy7b6fNyYHrNKt/RcQ4cOZejQoY6fJ06cyNGjR3n11Vf5r//6Lw9W5l1++MMfUlRUxCeffOLpUoRrfz/S0tKceg0mTpzI8OHD+c1vfsMLL7xws8vsEYYOHcq+ffuoqanh//7v/5gzZw4ffvhhu8FGrp96ajpYdHQ05eXlTvvKy8sJCQmhV69eHqpKvmz8+PF8/vnnni7Da8yfP5/c3Fzef/99BgwYcNW27X1GoqOjb2aJPYo778eX+fv7c/vtt+vz0YHMZjNDhgxhzJgx5OTkMHr0aH75y1+6bKvPx41RqOlgaWlp5OfnO+3Ly8tr9/qpeMa+ffuIiYnxdBndnmEYzJ8/nz/+8Y9s27aNwYMHf+Ux+ozcPNfzfnyZ1Wpl//79+nzcRDabjcbGRpeP6fNxgzw9Urmrq62tNfbu3Wvs3bvXAIxXXnnF2Lt3r3Hy5EnDMAxjyZIlxuOPP+5of+zYMaN3797GT37yE+PQoUPG2rVrDV9fX2Pr1q2eeglex9335NVXXzU2b95sHDlyxNi/f7+xcOFCw2QyGe+9956nXoLX+MEPfmCEhoYaH3zwgVFaWurY6uvrHW0ef/xxY8mSJY6f//a3vxl+fn7GL37xC+PQoUPGihUrDH9/f2P//v2eeAle5Xrej5UrVxrvvvuucfToUWPPnj3GzJkzjcDAQOPAgQOeeAleZ8mSJcaHH35oHD9+3Pj000+NJUuWGD4+PsZf//pXwzD0+ehoCjVf4dLtwF/e5syZYxiGYcyZM8eYNGlSm2NSUlIMs9ls3HLLLcZrr73W6XV7M3ffk5dfftm49dZbjcDAQCM8PNyYPHmysW3bNs8U72VcvQ+A05/5SZMmOd6bS/73f//XuO222wyz2WyMHDnS+POf/9y5hXup63k/Fi1aZAwcONAwm81GVFSUMXXqVKOwsLDzi/dS//Iv/2IMGjTIMJvNRr9+/YwpU6Y4Ao1h6PPR0XwMwzA6r19IRERE5ObQmBoRERHxCgo1IiIi4hUUakRERMQrKNSIiIiIV1CoEREREa+gUCMiIiJeQaFGREREvIJCjYiIiHgFhRoRERHxCgo1IiIi4hUUakRERMQrKNSIiIiIV/j/AZ7YOvJ+8uhYAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "from cfdmod.analysis.inflow.functions import calculate_mean_velocity\n", "\n", "mean_velocity_components = [\"ux\", \"uy\"]\n", "mean_velocity_data = calculate_mean_velocity(\n", " inflow_data=profile, for_components=mean_velocity_components\n", ")\n", "\n", "for component in mean_velocity_components:\n", " x = profile.points[\"z\"]\n", " y = mean_velocity_data[f\"{component}_mean\"]\n", " plt.plot(x, y, label=f\"{component}\")\n", "\n", "plt.legend()\n", "plt.show()\n", "plt.close()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Turbulence intensity" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjUAAAGdCAYAAADqsoKGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAABM/klEQVR4nO3deVhU9eLH8ffMsAwo4IKyKCru5QKuiGsZRWULt02t1MzSLE1Dr1fL1Fv2s2tZaXpFbdFs0WyxrillmjtiAprmmuKSCooKKMo65/eHN+7lRiYqHBg+r+c5D4/D98x8hrnTfO45Z75fi2EYBiIiIiIVnNXsACIiIiLXg0qNiIiIOAWVGhEREXEKKjUiIiLiFFRqRERExCmo1IiIiIhTUKkRERERp6BSIyIiIk7BxewAZcXhcHD8+HG8vLywWCxmxxEREZErYBgG586dIzAwEKv18sdiKk2pOX78OEFBQWbHEBERkatw9OhR6tate9kxlabUeHl5AZf+KN7e3ianERERkSuRmZlJUFBQ4ef45VSaUvPbKSdvb2+VGhERkQrmSi4d0YXCIiIi4hRUakRERMQpqNSIiIiIU6g019SIiIiUR4ZhkJ+fT0FBgdlRTGGz2XBxcbku062o1IiIiJgkNzeXEydOcOHCBbOjmMrT05OAgADc3Nyu6X5UakREREzgcDhITk7GZrMRGBiIm5tbpZsc1jAMcnNzOXXqFMnJyTRp0uRPJ9i7HJUaERERE+Tm5uJwOAgKCsLT09PsOKbx8PDA1dWVw4cPk5ubi91uv+r70oXCIiIiJrqWIxPO4nr9DfSXFBEREaegUiMiIiJOQaVGREREnIJKjYiIiDgFlZprdDG3gOcWb2P5jhNmRxEREanUVGqu0YebD/Nl0jH+umQ7+1PPmR1HREQqMMMwuJCbb8pmGMYV52zQoAFvvfVWkdtCQ0OZNGkSa9aswc3NjfXr1xf+burUqdSuXZvU1NTr9acqluapuUYDuzRg9Z6TxB08zZCFCSwd1gVvu6vZsUREpAK6mFfAjRO+NeWxd70UiafbtdeCm266iZEjR9KvXz+2b9/OwYMHefHFF1myZAl+fn7XIekf05Gaa+RiszLz4TYE+tg5mJbFqE+343BcedsVERFxNpMnT6Z69eoMHjyYRx99lAEDBnDPPfeU+uPqSM11ULOqO7MfbceDMXGs3JXKP9f8wrCeTcyOJSIiFYyHq41dL0Wa9tjXi5ubGx999BGtW7emfv36vPnmm9ftvi/nqo7UzJo1iwYNGmC32wkLC2PLli2XHb9kyRKaN2+O3W6nVatWLF++vMjvDcNgwoQJBAQE4OHhQUREBPv37y8yJjExkVtvvZVq1apRs2ZNBg8ezPnz568mfqkICarGy1EtAJi2ch9r9p40OZGIiFQ0FosFTzcXU7aSrDtltVp/dw1OXl5ekX9v2rQJgDNnznDmzJlr/+NcSa6S7rB48WKio6OZOHEiiYmJhISEEBkZycmTxX+Ib9q0ib59+zJo0CCSkpKIiooiKiqKnTt3Fo6ZOnUqM2bMICYmhvj4eKpUqUJkZCTZ2dkAHD9+nIiICBo3bkx8fDyxsbH8/PPPPPbYY1f3rEtJ7w716NuxHoYBIxZt48jpyr3qqoiIOKdatWpx4sR/vvWbmZlJcnJy4b8PHDjAc889x7x58wgLC2PAgAE4HI7SD2aUUMeOHY1nnnmm8N8FBQVGYGCgMWXKlGLHP/TQQ0avXr2K3BYWFmYMGTLEMAzDcDgchr+/v/Haa68V/j49Pd1wd3c3PvnkE8MwDGPOnDlG7dq1jYKCgsIxP/30kwEY+/fvv6LcGRkZBmBkZGRc2RO9Stl5+ca9MzcY9f+2zLj9rXXGhZz8Un08ERGpmC5evGjs2rXLuHjxotlRSmzs2LGGv7+/sW7dOuOnn34yoqKijKpVqxoTJ0408vPzjU6dOhn333+/YRiGcfz4caNmzZrG1KlT//D+Lve3KMnnd4mO1OTm5pKQkEBEREThbVarlYiICOLi4ordJy4ursh4gMjIyMLxycnJpKSkFBnj4+NDWFhY4ZicnBzc3NyKLHjl4eEBwIYNG4p93JycHDIzM4tsZcHdxcbsR9viW9WN3ScyGffFTyX6mpyIiEh5N27cOHr06MFdd91Fr169iIqKolGjRgC88sorHD58mDlz5gAQEBDA3LlzGT9+PNu3by/VXCUqNWlpaRQUFPzuK1l+fn6kpKQUu09KSsplx//283JjevbsSUpKCq+99hq5ubmcPXuWsWPHAhQ5/PXfpkyZgo+PT+EWFBRUkqd6TQJ8PJj5cFtsVgtLtx1nwaZDZfbYIiIipc3b25tFixaRkZHBkSNHGDBgANu2bWPSpElMmDCB48ePU7NmzcLx9913Hzk5OYSEhJRqrgrxle4WLVqwYMECpk2bhqenJ/7+/gQHB+Pn5/eHy5WPGzeOjIyMwu3o0aNlmrlTw5o8f+cNAEz+ZjdbksvmIikREZHKqkSlxtfXF5vN9rsZAVNTU/H39y92H39//8uO/+3nn93nww8/TEpKCseOHeP06dNMmjSJU6dO0bBhw2If193dHW9v7yJbWXu8SwPuCQkk32Hw9EeJpGZml3kGERGRyqJEpcbNzY127dqxatWqwtscDgerVq0iPDy82H3Cw8OLjAdYuXJl4fjg4GD8/f2LjMnMzCQ+Pr7Y+/Tz86Nq1aosXrwYu93OrbfeWpKnUKYsFguv3t+K5v5epJ3PYeiHCeTml8HV3yIiIpVQiU8/RUdHM2/ePBYsWMDu3bsZOnQoWVlZDBw4EID+/fszbty4wvEjRowgNjaWadOmsWfPHiZNmsTWrVsZNmwYcOmDf+TIkUyePJmvv/6aHTt20L9/fwIDA4mKiiq8n5kzZ5KYmMi+ffuYNWsWw4YNY8qUKVSrVu3a/gKlzNPNhTn92uFtdyHxSDovL9tldiQRERGnVOIZhXv37s2pU6eYMGECKSkphIaGEhsbW3ih75EjR4pc59K5c2c+/vhjxo8fz/PPP0+TJk1YunQpLVu2LBwzZswYsrKyGDx4MOnp6XTt2pXY2FjsdnvhmC1btjBx4kTOnz9P8+bNmTNnDv369buW515m6teswvQ+bXh8wY8s3HyY1nV9eLB92V24LCIi5Ze+IXv9/gYWo5L8NTMzM/Hx8SEjI8OU62sApn+/nze/34ebi5XPn+pMq7o+puQQERHzFRQUsG/fPmrXrl3km0KV0enTpzl58iRNmzbFZiu6XENJPr+19lMZGt6zMT/9ms6qPSd56sME/jW8KzWquJkdS0RETGCz2ahWrVrhjPyenp4lWqrAGRiGwYULFzh58iTVqlX7XaEpKR2pKWMZF/O4d+YGDp2+QJfGNVkwsCMutgrxzXoREbnODMMgJSWF9PR0s6OYqlq1avj7+xdb6kry+a1SY4K9KeeImrWRi3kFPNWjEWPvaG5qHhERMVdBQcHvFoSsLFxdXS97hEann8q5Zv5evPZga4Z9nETM2gOE1PXhjlYBZscSERGT2Gy2az71IhVkRmFndFfrQAZ3vzRx4Ogl29mfes7kRCIiIhWbSo2JxkQ2I7xhTbJyCxiyMIHM7Mp56FFEROR6UKkxkYvNysyH2xDoY+dgWhajPt2Ow1EpLnESERG57lRqTFazqjuzH22Hm83Kyl2pzF57wOxIIiIiFZJKTTkQElSNl6NaAPD6d3tZu++UyYlEREQqHpWacqJ3h3r07VgPw4BnP0ni6JkLZkcSERGpUFRqypFJ99xISFA1Mi7mMXhhAhdzC8yOJCIiUmGo1JQj7i42Yh5tS80qbuw+kcm4L37SQmciIiJXSKWmnAnw8WDmw22xWS0s3XacBZsOmR1JRESkQlCpKYfCG9Vk3L+XTpj8zW62JJ8xOZGIiEj5p1JTTg3qGszdIYHkOwye/iiR1MxssyOJiIiUayo15ZTFYuEf97eiub8XaedzGPphArn5DrNjiYiIlFsqNeWYp5sLc/q1w9vuQuKRdF5etsvsSCIiIuWWSk05V79mFab3aYPFAgs3H2bJ1qNmRxIRESmXVGoqgJub12bkLU0BeGHpTnYeyzA5kYiISPmjUlNBDO/ZmFua1yY338GQhQmcyco1O5KIiEi5olJTQVitFt7oHUqDmp4cS7/Is58kUaAVvUVERAqp1FQgPh6uzOnXHg9XGxt+SeP17/aaHUlERKTcUKmpYJr5ezH1gdYAzF5zgBU7TpicSEREpHxQqamA7g4J5MluwQCMXrKd/annTE4kIiJiPpWaCupvtzenU8MaZOUWMGRhApnZeWZHEhERMZVKTQXlYrMy8+G2BPjYOZiWxahPt+PQhcMiIlKJqdRUYL5V3Zn9aDvcbFZW7kpl9toDZkcSERExjUpNBRcaVI2X7m0BwOvf7WXtvlMmJxIRETGHSo0T6NOxHn071sMw4NlPkjh65oLZkURERMqcSo2TmHTPjYQEVSPjYh5DFiZwMbfA7EgiIiJlSqXGSbi72Ih5tC01q7ix60QmL3y5A8PQhcMiIlJ5qNQ4kQAfD2Y+3Bab1cIXScf4IO6w2ZFERETKjEqNkwlvVJNxdzQH4OVlu/jx0BmTE4mIiJQNlRonNKhrMHeHBJLvMHj6o0RSM7PNjiQiIlLqVGqckMVi4R/3t6KZnxenzuXw9EeJ5OY7zI4lIiJSqlRqnJSnmwtz+rXDy+5CwuGzvLxsl9mRRERESpVKjRNr4FuF6X1CAVi4+TBLth41N5CIiEgpUqlxcj2b+zEyogkALyzdyc5jGSYnEhERKR0qNZXAsz2bcEvz2uTmOxiyMIEzWblmRxIREbnuVGoqAavVwhu9Q2lQ05Nj6Rd59pMkCrSit4iIOBmVmkrCx8OVOf3a4+FqY8Mvabz+3V6zI4mIiFxXKjWVSDN/L6Y+0BqA2WsOsGLHCZMTiYiIXD8qNZXM3SGBPNktGIDRS7bzy8lzJicSERG5PlRqKqG/3d6cTg1rkJVbwOCFCZzLzjM7koiIyDVTqamEXGxWZj7clgAfOwdPZTHq0+04dOGwiIhUcCo1lZRvVXdmP9oON5uV73alMnvtAbMjiYiIXBOVmkosNKgaL93bAoDXv9vLun2nTE4kIiJy9VRqKrk+HevRt2MQhgHDP0ni6JkLZkcSERG5Kio1wqR7WhASVI2Mi3kMWvAjp87lmB1JRESkxFRqBHcXG7MfaUttL3f2pZ6n99w4UjKyzY4lIiJSIio1AkBgNQ8+HRJOnWoeHDyVxUNz4nQqSkREKhSVGinUwLcKi4d0ol4NT46cuUDvOXEkp2WZHUtEROSKqNRIEXWre/LpkHAa1arC8YxsHpoTx/5UzTosIiLl31WVmlmzZtGgQQPsdjthYWFs2bLlsuOXLFlC8+bNsdvttGrViuXLlxf5vWEYTJgwgYCAADw8PIiIiGD//v1Fxuzbt497770XX19fvL296dq1Kz/88MPVxJc/4e9jZ/GQcJr7e3HqXA69527m5+MZZscSERG5rBKXmsWLFxMdHc3EiRNJTEwkJCSEyMhITp48Wez4TZs20bdvXwYNGkRSUhJRUVFERUWxc+fOwjFTp05lxowZxMTEEB8fT5UqVYiMjCQ7+z8Xq951113k5+ezevVqEhISCAkJ4a677iIlJeUqnrb8Gd+q7iwa3InWdX04k5VL37mb2XY03exYIiIif8hiGEaJ5scPCwujQ4cOzJw5EwCHw0FQUBDDhw9n7Nixvxvfu3dvsrKyWLZsWeFtnTp1IjQ0lJiYGAzDIDAwkFGjRjF69GgAMjIy8PPzY/78+fTp04e0tDRq1arFunXr6NatGwDnzp3D29ublStXEhER8ae5MzMz8fHxISMjA29v75I85UotMzuPge//SMLhs1R1d+G9xzrQMbiG2bFERKSSKMnnd4mO1OTm5pKQkFCkRFitViIiIoiLiyt2n7i4uN+VjsjIyMLxycnJpKSkFBnj4+NDWFhY4ZiaNWvSrFkzPvjgA7KyssjPz2fOnDnUrl2bdu3aFfu4OTk5ZGZmFtmk5LztrnzweEfCG9bkfE4+A97bwob9aWbHEhER+Z0SlZq0tDQKCgrw8/Mrcrufn98fngZKSUm57Pjffl5ujMVi4fvvvycpKQkvLy/sdjtvvPEGsbGxVK9evdjHnTJlCj4+PoVbUFBQSZ6q/Jcq7i68P7ADPZrW4mJeAY8v+JHVe1LNjiUiIlJEhfj2k2EYPPPMM9SuXZv169ezZcsWoqKiuPvuuzlx4kSx+4wbN46MjIzC7ejRo2Wc2rnYXW3M7d+O2270IzffwZCFCazYUfzfXkRExAwlKjW+vr7YbDZSU4v+v/TU1FT8/f2L3cff3/+y43/7ebkxq1evZtmyZSxatIguXbrQtm1b/vnPf+Lh4cGCBQuKfVx3d3e8vb2LbHJt3F1szHqkLXeHBJJXYDDskySWJh0zO5aIiAhQwlLj5uZGu3btWLVqVeFtDoeDVatWER4eXuw+4eHhRcYDrFy5snB8cHAw/v7+RcZkZmYSHx9fOObChUsz21qtReNarVYcDkdJnoJcI1eblbd6h/JAu7oUOAye+3Qbi388YnYsERERXEq6Q3R0NAMGDKB9+/Z07NiRt956i6ysLAYOHAhA//79qVOnDlOmTAFgxIgR9OjRg2nTptGrVy8WLVrE1q1bmTt3LnDpepmRI0cyefJkmjRpQnBwMC+++CKBgYFERUUBl4pR9erVGTBgABMmTMDDw4N58+aRnJxMr169rtOfQq6UzWph6v2tsbta+XDzEf72+Q6y8xwM6NzA7GgiIlKJlbjU9O7dm1OnTjFhwgRSUlIIDQ0lNja28ELfI0eOFDmi0rlzZz7++GPGjx/P888/T5MmTVi6dCktW7YsHDNmzBiysrIYPHgw6enpdO3aldjYWOx2O3DptFdsbCwvvPACPXv2JC8vjxYtWvDVV18REhJyrX8DuQpWq4WX722J3cXGOxuSmfj1z2TnFTCkRyOzo4mISCVV4nlqKirNU1M6DMPgjZX7eHv1LwCMjGjCiFuaYLFYTE4mIiLOoNTmqRH5XxaLhVG3NeOvkc0AeOv7/fwjdi+VpCuLiEg5olIj18UzNzfmxbtuBCBm7QH+/q9dOBwqNiIiUnZUauS6GdQ1mFf+culaqfmbDvH8lzsoULEREZEyolIj19UjYfV5/cEQrBZY9ONRRi/ZTn6BvnYvIiKlT6VGrrsH2tVlRt82uFgtfJl0jOGfJJGbr2IjIiKlS6VGSsVdrQP55yNtcbNZWbEzhaEfJpCdV2B2LBERcWIqNVJqbmvhz7wB7XF3sbJqz0meWLCVC7n5ZscSEREnpVIjpapH01rMH9gRTzcbG35J47H3fuR8joqNiIhcfyo1UurCG9Vk4aAwvNxd2HLoDI++E0/GhTyzY4mIiJNRqZEy0a5+dT5+shPVPF3ZdjSdvvM2c/p8jtmxRETEiajUSJlpVdeHRYM74VvVjV0nMukzdzMnM7PNjiUiIk5CpUbKVHN/bxYPCcff287+k+fpPXczx9Mvmh1LREScgEqNlLlGtary6ZBw6lb3IDktiwdj4jhy+oLZsUREpIJTqRFT1KvpyadDwmlQ05Nj6Rd5cM4mfjl53uxYIiJSganUiGkCq3nw6ZBwmtSuSmpmDn3mxrEnJdPsWCIiUkGp1IipanvbWTS4EzcGeJN2Ppc+czez49cMs2OJiEgFpFIjpqtZ1Z1PnuxESFA10i/k8fC8zSQcPmN2LBERqWBUaqRc8PF05cNBHenYoAbncvLp9+4W4g6cNjuWiIhUICo1Um542V2Z/3gHujb25UJuAY+9v4W1+06ZHUtERCoIlRopVzzdXHhnQHt6Nq9NTr6DJxds5bufU8yOJSIiFYBKjZQ7dlcbMY+2446W/uQWOHj6o0T+tf242bFERKScU6mRcsnNxcrbfdsQFRpIvsNgxKIkPkv41exYIiJSjqnUSLnlYrMy7aFQ+nQIwmHA6CXb+XDzYbNjiYhIOaVSI+WazWrh//7Sisc6NwBg/NKdvLsh2dxQIiJSLqnUSLlntVqYePeNDOnREICXl+1i1g+/mJxKRETKG5UaqRAsFgtjb2/OyIgmALz27V5e/3YvhmGYnExERMoLlRqpMCwWCyMjmjL2juYAzPzhF175ZreKjYiIACo1UgE91aMRf7+nBQDvbEjmxa924nCo2IiIVHYqNVIhDejcgFfva4XFAh9uPsKYz3+iQMVGRKRSU6mRCqtPx3q8+VAoNquFzxJ+ZcSiJPIKHGbHEhERk6jUSIUW1aYOM/u2wcVqYdlPJ3jmo0Ry8gvMjiUiIiZQqZEK745WAczt3w43Fyvf7Upl8AcJZOep2IiIVDYqNeIUejb3470BHbC7Wlm77xT9393CqXM5ZscSEZEypFIjTqNrE18+eDyMqu4ubDl0hjtnrCfuwGmzY4mISBlRqRGn0jG4Bkuf6UxTv6qcOpfDI+9sZtYPv+gr3yIilYBKjTidxrW9WPpMFx5oVxeHcWn24YHzf+RMVq7Z0UREpBSp1IhT8nRz4fUHQ5j6QOvC62zunL6erYfOmB1NRERKiUqNOLWH2gex9JkuNKxVhZTMbHrP3cyctQd0OkpExAmp1IjTa+7vzdfDunJvaCAFDoMpK/YweOFW0i/odJSIiDNRqZFKoaq7C2/1DuWVv7TEzcXK97tP0mvGBpKOnDU7moiIXCcqNVJpWCwWHgmrzxdDO1O/pifH0i/y0Jw43tuQrJW+RUScgEqNVDot6/jwr+FdubOVP3kFBi8t28XQDxPJuJhndjQREbkGKjVSKXnbXZn1cFsm3X0jrjYLsT+ncPfbG9h5LMPsaCIicpVUaqTSslgsPNYlmM+e6kzd6h4cOXOB+/65iYWbD+t0lIhIBaRSI5VeSFA1vhnejYgb/MgtcPDi0p0M/ySJ8zn5ZkcTEZESUKkRAXw8XZnXvx0v3HkDLlYLy346wT1vb2D3iUyzo4mIyBVSqRH5N4vFwpPdG7J4SDgBPnYOpmURNWsji388otNRIiIVgEqNyP9oV7863zzbjZua1SIn38HfPt/BqE+3cyFXp6NERMozlRqRYtSo4sZ7Azrw18hmWC3wRdIx7p25kf2p58yOJiIif0ClRuQPWK0Wnrm5MZ882YnaXu7sP3mee2Zu5IvEX82OJiIixVCpEfkTYQ1rsnxEN7o29uViXgHRn27nb5/9RHZegdnRRETkv6jUiFwB36ruLHi8I89FNMVigcVbjxI1ayMHT503O5qIiPybSo3IFbJZLYyIaMKHg8LwrerGnpRz3P32Br7eftzsaCIiwlWWmlmzZtGgQQPsdjthYWFs2bLlsuOXLFlC8+bNsdvttGrViuXLlxf5vWEYTJgwgYCAADw8PIiIiGD//v2Fv1+zZg0Wi6XY7ccff7yapyBy1bo09mX5s90IC65BVm4Bz36SxPilO3Q6SkTEZCUuNYsXLyY6OpqJEyeSmJhISEgIkZGRnDx5stjxmzZtom/fvgwaNIikpCSioqKIiopi586dhWOmTp3KjBkziImJIT4+nipVqhAZGUl2djYAnTt35sSJE0W2J554guDgYNq3b3+VT13k6tX2tvPRE2EMu7kxAB9uPsIDMZs4fDrL5GQiIpWXxSjhrGJhYWF06NCBmTNnAuBwOAgKCmL48OGMHTv2d+N79+5NVlYWy5YtK7ytU6dOhIaGEhMTg2EYBAYGMmrUKEaPHg1ARkYGfn5+zJ8/nz59+vzuPvPy8qhTpw7Dhw/nxRdfvKLcmZmZ+Pj4kJGRgbe3d0messhlrdl7kucWb+PshTy83F147cHW3N4ywOxYIiJOoSSf3yU6UpObm0tCQgIRERH/uQOrlYiICOLi4ordJy4ursh4gMjIyMLxycnJpKSkFBnj4+NDWFjYH97n119/zenTpxk4cOAfZs3JySEzM7PIJlIabmpWm2+e7Ua7+tU5l5PPUx8m8vd//UxuvsPsaCIilUqJSk1aWhoFBQX4+fkVud3Pz4+UlJRi90lJSbns+N9+luQ+3333XSIjI6lbt+4fZp0yZQo+Pj6FW1BQ0OWfnMg1CKzmwaLBnRjSvSEA7288xINz4vj17AWTk4mIVB4V7ttPv/76K99++y2DBg267Lhx48aRkZFRuB09erSMEkpl5WqzMu7OG3inf3t8PFzZfjSdXjM28P2uVLOjiYhUCiUqNb6+vthsNlJTi/5HOjU1FX9//2L38ff3v+z4335e6X2+//771KxZk3vuueeyWd3d3fH29i6yiZSFiBv9WDa8KyFB1ci4mMcTH2xlyvLd5BXodJSISGkqUalxc3OjXbt2rFq1qvA2h8PBqlWrCA8PL3af8PDwIuMBVq5cWTg+ODgYf3//ImMyMzOJj4//3X0ahsH7779P//79cXV1LUl0kTIVVMOTJUPCGdilAQBz1h2kz9zNnMi4aG4wEREnVuLTT9HR0cybN48FCxawe/duhg4dSlZWVuFFu/3792fcuHGF40eMGEFsbCzTpk1jz549TJo0ia1btzJs2DAALBYLI0eOZPLkyXz99dfs2LGD/v37ExgYSFRUVJHHXr16NcnJyTzxxBPX8JRFyoabi5WJd7dg9iNt8XJ3IeHwWe6cvp41e4uf/kBERK6NS0l36N27N6dOnWLChAmkpKQQGhpKbGxs4YW+R44cwWr9T1fq3LkzH3/8MePHj+f555+nSZMmLF26lJYtWxaOGTNmDFlZWQwePJj09HS6du1KbGwsdru9yGO/++67dO7cmebNm1/t8xUpc3e0CuDGQG+e/iiRn49n8tj7PzLs5saMjGiCi63CXdYmIlJulXiemopK89SI2bLzCpj8zS4+3HwEgLDgGrzdtw21ve1/sqeISOVVavPUiMjVs7vamBzVihl921DFzUZ88hnunLGejb+kmR1NRMQpqNSIlLF7QgL5enhXmvt7kXY+l0ffjWf69/spcFSKg6YiIqVGpUbEBI1qVeXLp7vQu30QhgFvfr+PAe9tIe18jtnRREQqLJUaEZN4uNn4xwOtmfZgCB6uNjb8ksad09cTf/C02dFERCoklRoRk93fri5fDetC49pVOXkuh77zNjPrh19w6HSUiEiJqNSIlANN/bz46pku3NemDg4DXvt2L48v+JGzWblmRxMRqTBUakTKiSruLkx7KIRX72uFu4uVNXtPcfv0dSzfcYJKMvOCiMg1UakRKUcsFgt9OtZj6TNdaOhbhdTMHJ7+KJEB7//IobQss+OJiJRrKjUi5dANAd4sH9GNZ29pgpvNyrp9p7jtrXW8uXIf2XkFZscTESmXVGpEyim7q43oW5vy7XPd6dbEl9x8B9NX7SfyrXX8oPWjRER+R6VGpJwL9q3CB493ZNbDbfHzdufw6QsMfP9Hhn6YwPF0rfotIvIblRqRCsBisdCrdQCrRt3EE12DsVktrNiZQsQba5m77gB5BQ6zI4qImE4LWopUQLtPZDJ+6U4SDp8FoKlfVV6+tyVhDWuanExE5PrSgpYiTu6GAG+WDAln6gOtqVHFjX2p5+k9dzPRn27TUgsiUmmp1IhUUFarhYfaB7F6VA/6dqyHxQJfJB6j5+trWLj5sBbIFJFKR6efRJxE4pGzvLh0Jz8fzwSgdV0fJke1pHXdauYGExG5BiX5/FapEXEi+QUOPtx8mGnf7eNcTj4WCzwaVp/RtzXDx9PV7HgiIiWma2pEKikXm5XHugSzanQPokIDMQxYuPkwt7yxhi8Sf9VyCyLi1HSkRsSJbTqQxotLd3Lg1KUlFjoG12ByVEua+nmZnExE5MroSI2IANC5kS8rRnRnzO3NsLta2ZJ8hjunr2fK8t1k5eSbHU9E5LpSqRFxcm4uVp6+qTHfR/fg1hv9yHcYzFl3kIg31hK7UyuAi4jzUKkRqSTqVvdkXv/2vNO/PXWre3AiI5unPkxk4PwfOXxaK4CLSMWnUiNSyUTc6MfK53ow7ObGuNosrNl7ilvfXMf07/drBXARqdBUakQqIQ83G6MjmxE7sjtdGtckN9/Bm9/v4/a31rF23ymz44mIXBWVGpFKrFGtqnw4KIy3+7ahtpc7h05fYMB7W3j6owROZGgFcBGpWFRqRCo5i8XC3SGBrBrVg8e7BGO1wPIdKURMW8s76w9qBXARqTA0T42IFPHz8QxeXLqTxCPpADTz82LyX1rSoUENc4OJSKWkeWpE5Kq1CPThs6c684/7W1HN05W9qed4MCaO0Uu2c1orgItIOaZSIyK/Y7Va6N2hHqtH3USfDkEAfJbwKz2nreWj+MM4tAK4iJRDOv0kIn8q4fBZxi/dye4Tl1YADwmqxuR7W9Kqro/JyUTE2WmV7mKo1Ihcm/wCBwv/vQL4+Zx8rBbo16k+0bc1w8dDK4CLSOnQNTUict252KwM7BLMqlE9uCckEIcBC+IOc8u0tSxNOqblFkTEdCo1IlIift52ZvRtw0dPhNGwVhXSzucwcvE2+s7bzP7Uc2bHE5FKTKVGRK5Kl8a+rBjRjb9GNsPdxcrmg2e4Y/p6Xl2xhwu5WgFcRMqeSo2IXDV3FxvP3HxpBfCIG2qT7zCIWXuAW99Yx7c/p+iUlIiUKZUaEblmQTU8eWdAB+b1b0+dah4cS7/IkIUJDFqwlSOnL5gdT0QqCZUaEblubr3Rj5XR3Xn6pka42iys3nOSW99cy9ur9mu5BREpdSo1InJdebq5MOb25qwY0Z3OjWqSk+9g2sp99J4Tx9EzOmojIqVHpUZESkXj2lX56Ikw3uwdgpe7C4lH0rlzxnr+tf242dFExEmp1IhIqbFYLPylTV2Wj+hGm3rVOJedz/BPkhjz2XZ9Q0pErjuVGhEpdUE1PPl0SDjDbm6MxQKfbv2Vu97ewM/HM8yOJiJORKVGRMqEq83K6MhmfPREGH7e7hw8lcVfZm3ivQ3J+uq3iFwXKjUiUqY6N/JlxYjuRNzgR26Bg5eW7WLQgq2cPp9jdjQRqeBUakSkzNWo4sa8/u146d4WuLlYWb3nJLdPX8+G/WlmRxORCkylRkRMYbFY6B/egK+e6ULj2lU5dS6Hfu/F8+qKPZrTRkSuikqNiJjqhgBv/jWsK3071sMwIGbtAR6IiePw6Syzo4lIBaNSIyKm83CzMeW+Vsx+pC3edhe2H02n14wNLE06ZnY0EalAVGpEpNy4o1UAK0Z2p0OD6pzPyWfk4m1Ef7qN8zma00ZE/pxKjYiUK3WqefDJk50YcUsTrBb4IvEYd81Yz45fNaeNiFyeSo2IlDsuNivP3dqURYPDCfSxc+j0Be6bvZF56w7icGhOGxEpnkqNiJRbHYNrsHxEN25v4U9egcEry3fz2PwfOXVOc9qIyO+p1IhIuVbN043Zj7bllb+0xN3Fyrp9p7hj+jrW7jtldjQRKWdUakSk3LNYLDwSVp9/De9KMz8v0s7nMuC9LbzyzS5y8zWnjYhcclWlZtasWTRo0AC73U5YWBhbtmy57PglS5bQvHlz7HY7rVq1Yvny5UV+bxgGEyZMICAgAA8PDyIiIti/f//v7uebb74hLCwMDw8PqlevTlRU1NXEF5EKqqmfF18N60L/8PoAzFufzP2zN5GcpjltROQqSs3ixYuJjo5m4sSJJCYmEhISQmRkJCdPnix2/KZNm+jbty+DBg0iKSmJqKgooqKi2LlzZ+GYqVOnMmPGDGJiYoiPj6dKlSpERkaSnZ1dOObzzz+nX79+DBw4kO3bt7Nx40Yefvjhq3jKIlKR2V1tvHRvS+b2a0c1T1d2HMug14z1fJ7wqxbGFKnkLEYJ/ysQFhZGhw4dmDlzJgAOh4OgoCCGDx/O2LFjfze+d+/eZGVlsWzZssLbOnXqRGhoKDExMRiGQWBgIKNGjWL06NEAZGRk4Ofnx/z58+nTpw/5+fk0aNCAv//97wwaNOiqnmhmZiY+Pj5kZGTg7e19VfchIuXLiYyLPLd4G5sPngHg3tBAJke1xMvuanIyEbleSvL5XaIjNbm5uSQkJBAREfGfO7BaiYiIIC4urth94uLiiowHiIyMLByfnJxMSkpKkTE+Pj6EhYUVjklMTOTYsWNYrVbatGlDQEAAd9xxR5GjPf8rJyeHzMzMIpuIOJcAHw8+eqITo29ris1q4attx+k1YwNJR86aHU1ETFCiUpOWlkZBQQF+fn5Fbvfz8yMlJaXYfVJSUi47/reflxtz8OBBACZNmsT48eNZtmwZ1atX56abbuLMmTPFPu6UKVPw8fEp3IKCgkryVEWkgrBZLQzr2YRPh3SiTjUPjpy5wIMxccxec0Bz2ohUMhXi208Ox6VvN7zwwgvcf//9tGvXjvfffx+LxcKSJUuK3WfcuHFkZGQUbkePHi3LyCJSxtrVvzSnTa/WAeQ7DP4Ru4f+723hZGb2n+8sIk6hRKXG19cXm81GampqkdtTU1Px9/cvdh9/f//Ljv/t5+XGBAQEAHDjjTcW/t7d3Z2GDRty5MiRYh/X3d0db2/vIpuIODcfD1dm9m3D1Ptb4+FqY8Mvadw+fT2r96T++c4iUuGVqNS4ubnRrl07Vq1aVXibw+Fg1apVhIeHF7tPeHh4kfEAK1euLBwfHByMv79/kTGZmZnEx8cXjmnXrh3u7u7s3bu3cExeXh6HDh2ifv36JXkKIuLkLBYLD3UI4l/Du3JDgDdnsnJ5fP5W/v6vn8nJLzA7noiUohKffoqOjmbevHksWLCA3bt3M3ToULKyshg4cCAA/fv3Z9y4cYXjR4wYQWxsLNOmTWPPnj1MmjSJrVu3MmzYMODSf4BGjhzJ5MmT+frrr9mxYwf9+/cnMDCwcB4ab29vnnrqKSZOnMh3333H3r17GTp0KAAPPvjgtf4NRMQJNa5dlS+f7szALg0AeH/jIf4yaxO/nDxvbjARKTUuJd2hd+/enDp1igkTJpCSkkJoaCixsbGFF/oeOXIEq/U/Xalz5858/PHHjB8/nueff54mTZqwdOlSWrZsWThmzJgxZGVlMXjwYNLT0+natSuxsbHY7fbCMa+99houLi7069ePixcvEhYWxurVq6levfq1PH8RcWJ2VxsT725Btya+jF7yE7tOZHL32xuYdM+NPNQ+CIvFYnZEEbmOSjxPTUWleWpEKreTmdk89+k2Nv5yGoBerQP4v7+0wsdDc9qIlGelNk+NiEhFVdvbzsLHw/jb7c1xsVr45qcT3Dl9PQmHi58WQkQqHpUaEak0rFYLQ29qxJKnwgmq4cGx9Is8NGczM1fvp0Bz2ohUeCo1IlLptKlXneXPduPe0EAKHAavf7ePR97ZzImMi2ZHE5FroFIjIpWSl92Vt3qHMu3BEDzdbGw+eIY7pq/nu5+Lnx1dRMo/lRoRqbQsFgv3t6vLsuFdaVnHm/QLeQxemMCEr3aSnac5bUQqGpUaEan0GtaqyhdDu/Bkt2AAPog7TNSsjexPPWdyMhEpCZUaERHAzcXKC71uZMHjHfGt6saelHPc9fYGPoo/TCWZ+UKkwlOpERH5Lz2a1mLFiO50b1qLnHwHL3y5k6EfJpJ+IdfsaCLyJ1RqRET+Ry0vd+Y/1oEX7rwBV5uF2J9TuHP6erYka04bkfJMpUZEpBhWq4Unuzfki6FdaFDTk+MZ2fSZG8c/YvdoYUyRckqlRkTkMlrV9WHZs924v21dHAbMXnOAe2duZNfxTLOjicj/UKkREfkTVd1dmPZQCDGPtqVmlUsXEd87awMzV+8nv8BhdjwR+TeVGhGRK3R7ywC+fa47kS38yCu4NBPx/TFxHDh13uxoIoJKjYhIifhWdSfm0Xa88VAIXnYXth9N587p63lvQzIOrR8lYiqVGhGRErJYLNzXti7fPdedbk18ycl38NKyXTzyTjy/nr1gdjyRSkulRkTkKgX4ePDB4x15OaolHq424g6e5va31rP4xyOasE/EBCo1IiLXwGKx0K9TfVaM6Eb7+tU5n5PP3z7fwRMLtnIyM9vseCKVikqNiMh10MC3CouHhDPujua42ays2nOS295ax7+2Hzc7mkiloVIjInKd2KwWhvRoxLJn/7Pq9/BPkhj2cSJns7TMgkhpU6kREbnOmvp58eXTXXj2libYrBaW/XSC295ax+o9qWZHE3FqKjUiIqXA1WYl+tamfDG0M41qVeHUuRwen7+Vv332E+ey88yOJ+KUVGpEREpRSFA1vnm2G090DcZigcVbj3L7W+vZdCDN7GgiTkelRkSklNldbYy/60YWPdmJoBoeHEu/yMPz4vn7v34mO0+LY4pcLyo1IiJlJKxhTVaM6E7fjvUAeH/jIe6csZ6kI2dNTibiHFRqRETKUFV3F6bc14r3B3agtpc7B09lcf/sTbz+7V5y87U4psi1UKkRETHBzc1q891z3YkKDcRhwMwffiFq1kb2pGSaHU2kwlKpERExSTVPN97q04Z/PtKW6p6u7DqRyd1vb+Cfa36hQItjipSYSo2IiMnubBXAd8/1IOIGP/IKDKbG7uWBmE0kp2WZHU2kQlGpEREpB2p5uTOvfztefzAEL3cXko6kc8f0dSzYdAiHjtqIXBGVGhGRcsJisfBAu7rEPtedLo1rkp3nYOLXP9PvvXiOpV80O55IuadSIyJSztSp5sHCx8N46d4W2F2tbPzlNLe/uY4lW49iGDpqI/JHVGpERMohq9VC//AGrBjRnbb1qnEuJ5+/fvYTT36QwMlz2WbHEymXVGpERMqxYN8qLHmqM3+7vTluNivf704l8s11LN9xwuxoIuWOSo2ISDlns1oYelMjvh7ehRsCvDl7IY+nP0pkxKIk0i/kmh1PpNxQqRERqSCa+3vz1TNdGHZzY6wW+GrbcSLfWscPe0+aHU2kXFCpERGpQNxcrIyObMbnQzvTsFYVUjNzGPj+j4z7Ygfnc/LNjidiKpUaEZEKqE296ix/thuPdwkG4JMtR7hj+jriD542OZmIeVRqREQqKLurjQl338gnT3aiTjUPjp65SJ95m3l52S6y8wrMjidS5lRqREQquPBGNYkd2Y3e7YMwDHh3QzJ3vb2Bn35NNzuaSJlSqRERcQJedlf+8UBr3nusPbW83Pnl5Hn+8s9NvLFyH3kFDrPjiZQJlRoRESfSs7kf343szt0hgRQ4DGas2k/UrI3sTTlndjSRUqdSIyLiZKpXcePtvm2Y+XAbqnm68vPxTO5+ewNz1h6gQItjihNTqRERcVJ3tQ7ku5Hd6dm8NrkFDqas2EPvOXEcSssyO5pIqVCpERFxYrW97bw7oD1T729NVXcXth4+yx3T1/PehmQdtRGno1IjIuLkLBYLD3UIYsWIboQ3rMnFvAJeWraL+2dv0rU24lRUakREKomgGp589EQYr/ylJV7uLmw7ms5db6/njZX7yMnXvDZS8anUiIhUIlarhUfC6rMyuge33uhHXsGlb0j1mrGBhMNnzI4nck1UakREKiF/Hztz+7Xjn4+0xbfqpXltHoiJY8JXO7WGlFRYKjUiIpWUxWLhzlYBfB/dnYfa18Uw4IO4w9z6xlpW70k1O55IianUiIhUctU83Zj6QAgfPRFGvRqenMjI5vH5Wxn+SRJp53PMjidyxVRqREQEgC6Nffl2ZHeGdG+I1QL/2n6ciDfW8nnCrxiGvv4t5Z9KjYiIFPJwszHuzhv46pmu3BDgTfqFPEYt2U7/97Zw9MwFs+OJXJZKjYiI/E6ruj58PawLY25vhpuLlfX707jtzXW8q0n7pBy7qlIza9YsGjRogN1uJywsjC1btlx2/JIlS2jevDl2u51WrVqxfPnyIr83DIMJEyYQEBCAh4cHERER7N+/v8iYBg0aYLFYimyvvvrq1cQXEZEr4Gqz8vRNjfl2ZHfCgmtwMa+Al5ft4r7Zm9iTkml2PJHfKXGpWbx4MdHR0UycOJHExERCQkKIjIzk5MmTxY7ftGkTffv2ZdCgQSQlJREVFUVUVBQ7d+4sHDN16lRmzJhBTEwM8fHxVKlShcjISLKzs4vc10svvcSJEycKt+HDh5c0voiIlFCwbxU+ebITU+5rhZe7C9uPpnPXjA1M+24v2XmatE/KD4tRwqu/wsLC6NChAzNnzgTA4XAQFBTE8OHDGTt27O/G9+7dm6ysLJYtW1Z4W6dOnQgNDSUmJgbDMAgMDGTUqFGMHj0agIyMDPz8/Jg/fz59+vQBLh2pGTlyJCNHjryqJ5qZmYmPjw8ZGRl4e3tf1X2IiFR2qZnZvLh0J9/tuvSV70a1qvDq/a3p0KCGycnEWZXk87tER2pyc3NJSEggIiLiP3dgtRIREUFcXFyx+8TFxRUZDxAZGVk4Pjk5mZSUlCJjfHx8CAsL+919vvrqq9SsWZM2bdrw2muvkZ//xxNE5eTkkJmZWWQTEZFr4+dtZ27/9sx+pC21vNw5cCqLB2PieHHpTs5l55kdTyq5EpWatLQ0CgoK8PPzK3K7n58fKSkpxe6TkpJy2fG//fyz+3z22WdZtGgRP/zwA0OGDOH//u//GDNmzB9mnTJlCj4+PoVbUFDQlT9RERG5rDtaBfD9cz3o3f7Sf1sXbj7MrW+s4/tdmrRPzFNhvv0UHR3NTTfdROvWrXnqqaeYNm0ab7/9Njk5xU8MNW7cODIyMgq3o0ePlnFiERHn5uPpyj8eaM3HT4RRv6YnKZnZPPHBVoZ9nMipc5q0T8peiUqNr68vNpuN1NSiTTw1NRV/f/9i9/H397/s+N9+luQ+4dK1Pfn5+Rw6dKjY37u7u+Pt7V1kExGR669zY19iR3RnSI+G2KwWlv10gog31vKZJu2TMlaiUuPm5ka7du1YtWpV4W0Oh4NVq1YRHh5e7D7h4eFFxgOsXLmycHxwcDD+/v5FxmRmZhIfH/+H9wmwbds2rFYrtWvXLslTEBGRUuDhZmPcHTfw1TNduDHAm4yLeYzWpH1SxlxKukN0dDQDBgygffv2dOzYkbfeeousrCwGDhwIQP/+/alTpw5TpkwBYMSIEfTo0YNp06bRq1cvFi1axNatW5k7dy5waUG1kSNHMnnyZJo0aUJwcDAvvvgigYGBREVFAZcuNo6Pj+fmm2/Gy8uLuLg4nnvuOR599FGqV69+nf4UIiJyrVrW8eGrYV14d0Myb67cVzhp36jbmjKwSzA2q8XsiOLESlxqevfuzalTp5gwYQIpKSmEhoYSGxtbeKHvkSNHsFr/cwCoc+fOfPzxx4wfP57nn3+eJk2asHTpUlq2bFk4ZsyYMWRlZTF48GDS09Pp2rUrsbGx2O124NKppEWLFjFp0iRycnIIDg7mueeeIzo6+lqfv4iIXGeuNitP9WhEZAt/xn3xE5sPnmHyN7v51/bjvHp/a24I0OUAUjpKPE9NRaV5akREyp5hGCz+8SivLN/Nuex8XKwWhvRoyPCeTbC72syOJxVAqc1TIyIiUhIWi4U+HeuxKroHt7fwJ99hMOuHA9w5fT3xB0+bHU+cjEqNiIiUutredmL6tSPm0bbU9nLnYFoWvedu5oUvd5CpSfvkOlGpERGRMnN7ywBWRvegb8dLk/Z9FH+E295Yx0pN2ifXgUqNiIiUKR8PV6bc15qPnwyjwb8n7Xvyg60885Em7ZNro1IjIiKm6NzIl9iR3XmqRyNsVgvf7Lg0ad+nW49q0j65Kio1IiJiGrurjbF3NOerZ7rQss6lSfvGfPYTj74bz5HTmrRPSkalRkRETNeyjg9Ln+7CuDua4+5iZeMvp7ntrbXMW3eQ/AKH2fGkglCpERGRcsHFZmVIj0Z8O7I74Q1rkp3n4JXlu7lv9iZ2Hc80O55UACo1IiJSrjTwrcLHT4Yx9f7WeNtd+OnXDO6euYGpsXvIziswO56UYyo1IiJS7lgsFh7qEMT30T24o6U/BQ6Df665NGnfZk3aJ39ApUZERMqt2t52Zj/ajphH2xVO2tdn7mbGfaFJ++T3VGpERKTcu72l/78n7asHwCdbjnDrG2tZvuOEvv4thVRqRESkQrg0aV8rFg3uRLBvFVIzc3j6o0T6v7eF5LQss+NJOaBSIyIiFUqnhjVZMaIbz97SBDeblfX704h8cx3TvturC4krOZUaERGpcOyuNqJvbcp3z3WnR9Na5BY4eHv1L9z65lpW7dY6UpWVSo2IiFRYDXyrMH9gB2IebUuAj52jZy4yaMFWnliwlaNnNCNxZaNSIyIiFZrFYuH2lgF8H92DIT0a4mK18P3uVG59cy0zV+8nJ1+npCoLlRoREXEKVdxdGHfHDawY0Y1ODWuQnefg9e/2ccdb69mwP83seFIGVGpERMSpNPHz4pMnOzG9Tyi+VS/NbfPou/E883EiKRnZZseTUqRSIyIiTsdisXBvaB1Wj+7BY50bYLXANz+d4JZpa5i37iB5WiTTKVmMSjJrUWZmJj4+PmRkZODt7W12HBERKUM/H8/gxaU7STySDkAzPy9ejmpJx+Aa5gaTP1WSz28dqREREafXItCHz57qzNT7W1Pd05W9qed4aE4c0Z9u49S5HLPjyXWiUiMiIpWC1XppkczVo26ib8d6WCzwReIxek5bwwdxhyhwVIoTF05Np59ERKRS2nY0nfFLd7DzWCYALet48/K9LWlTr7rJyeS/leTzW6VGREQqrQKHwcfxh5n67V7OZedjsUCfDvUYE9mM6lXczI4n6JoaERGRK2KzWugX3oDVo27ivrZ1MIxLK4D3nLaGxT8ewaFTUhWKjtSIiIj825bkM7y4dCd7U88B0LZeNV6OakmLQB+Tk1VeOv1UDJUaERG5EnkFDhZsOsSbK/eRlVuA1QL9wxsQfVtTvO2uZserdHT6SURE5Cq52qw80a0hq0bdxF2tA3AYMH/TIW6ZtpalSceoJMcCKiQdqREREbmMDfvTmPDVTg6mZQHQqWENXr63JU38vExOVjnoSI2IiMh10rWJLytGduOvkc2wu1rZfPAMd0xfz5QVu8nKyTc7nvwXlRoREZE/4e5i45mbG7PyuR7ceqMf+Q6DOWsPEvHGWlbsOKFTUuWESo2IiMgVCqrhybz+7Xl3QHvqVvfgREY2Qz9KZMD7P5L879NTYh6VGhERkRK65QY/vo/uwbM9G+Nms7Ju3yki31zHGyv3kZ1XYHa8SkulRkRE5CrYXW1E39aMb5/rTrcmvuQWOJixaj+3vrmW1XtSzY5XKanUiIiIXINg3yp88HhH/vlIW/y97Rw9c5HH52/lyQ+28uvZC2bHq1RUakRERK6RxWLhzlYBrBrVgyHdG+JitbByVyoRb6xl1g+/kJOvU1JlQfPUiIiIXGf7Us/x4tKdxCefAaBhrSq8fG9LujT2NTlZxaN5akREREzU1M+LRYM78VbvUHyrunPwVBaPvBPPsI8TSc3MNjue01KpERERKQUWi4WoNnVYNaoHj3VugNUCy346Qc/X1/DO+oPkFTjMjuh0dPpJRESkDOw8lsGLX+0k6Ug6AM39vXg5qiUdGtQwN1g5p1W6i6FSIyIiZnM4DJYkHOXVFXs4eyEPgPvb1uWvkc3w97GbnK58UqkphkqNiIiUF2ezcpn67R4+2XIUABerhXtCAxncvSHN/fUZ9d9UaoqhUiMiIuVN0pGzvLpiT+G3pAB6NK3F4O4N6dyoJhaLxcR05YNKTTFUakREpLzafjSduesPsmLHCRz//lRuEejN4O4NubNVAK62yvu9HpWaYqjUiIhIeXfk9AXe25jM4h+PcvHfa0jVqebB412D6d0hiKruLiYnLHsqNcVQqRERkYribFYuH24+zIK4Q6SdzwXA2+7CI53qM7BzA2p7V56LilVqiqFSIyIiFU12XgFfJh1j3rqDHEzLAsDNZuXef19U3MTPy+SEpU+lphgqNSIiUlE5HAbf705l7rqDbD18tvD2ns1rM7h7Q8KCazjtRcUqNcVQqREREWeQcPgs89Yd5NtdKfz2CR5S14cnuzfk9hb+uDjZRcUqNcVQqREREWeSnJbFuxsOsmTrr+TkX1pyIaiGB4O6BPNQhyA83ZzjomKVmmKo1IiIiDM6fT6HD+IO80HcocJZiqt5utKvU336hzeglpe7yQmvTamv0j1r1iwaNGiA3W4nLCyMLVu2XHb8kiVLaN68OXa7nVatWrF8+fIivzcMgwkTJhAQEICHhwcRERHs37+/2PvKyckhNDQUi8XCtm3bria+iIiI06hZ1Z3nbm3KprG38HJUS+rX9CT9Qh5vr/6FLv9YzbgvfuLAqfNmxywTJS41ixcvJjo6mokTJ5KYmEhISAiRkZGcPHmy2PGbNm2ib9++DBo0iKSkJKKiooiKimLnzp2FY6ZOncqMGTOIiYkhPj6eKlWqEBkZSXb275dnHzNmDIGBgSWNLSIi4tQ83Gz061Sf1aNuYvYjbQkNqkZuvoNPthzllmlreWLBVn48dAZnPkFT4tNPYWFhdOjQgZkzZwLgcDgICgpi+PDhjB079nfje/fuTVZWFsuWLSu8rVOnToSGhhITE4NhGAQGBjJq1ChGjx4NQEZGBn5+fsyfP58+ffoU7rdixQqio6P5/PPPadGiBUlJSYSGhl5Rbp1+EhGRysQwDLYePsuctQf5fndq4e1t6lVjcLeG3NbCH5u1/H9jqtROP+Xm5pKQkEBERMR/7sBqJSIigri4uGL3iYuLKzIeIDIysnB8cnIyKSkpRcb4+PgQFhZW5D5TU1N58sknWbhwIZ6eniWJLSIiUulYLBY6NKjBOwPa8310D/p2DMLNxUrSkXSGfpTILdPWsHDzYS7mFpgd9bopUalJS0ujoKAAPz+/Irf7+fmRkpJS7D4pKSmXHf/bz8uNMQyDxx57jKeeeor27dtfUdacnBwyMzOLbCIiIpVR49pVmXJfazb+rSfDbm6Mj4crh05f4MWlO+nyj9W8uXIfp8/nmB3zmlWIL7O//fbbnDt3jnHjxl3xPlOmTMHHx6dwCwoKKsWEIiIi5V8tL3dGRzZj09ieTLr7RupW9+BMVi7TV+2n86ureeHLHST/e+biiqhEpcbX1xebzUZqamqR21NTU/H39y92H39//8uO/+3n5casXr2auLg43N3dcXFxoXHjxgC0b9+eAQMGFPu448aNIyMjo3A7evRoSZ6qiIiI06ri7sJjXYJZM/om3u7bhlZ1fMjJd/BR/BF6TlvDUwsTSPivmYsrihKVGjc3N9q1a8eqVasKb3M4HKxatYrw8PBi9wkPDy8yHmDlypWF44ODg/H39y8yJjMzk/j4+MIxM2bMYPv27Wzbto1t27YVfiV88eLFvPLKK8U+rru7O97e3kU2ERER+Q8Xm5W7QwL5elgXPnmyEzc3q4VhQOzPKdw/exMPzN7Edz+n4HBUjG9MlXi6wejoaAYMGED79u3p2LEjb731FllZWQwcOBCA/v37U6dOHaZMmQLAiBEj6NGjB9OmTaNXr14sWrSIrVu3MnfuXODShUwjR45k8uTJNGnShODgYF588UUCAwOJiooCoF69ekUyVK1aFYBGjRpRt27dq37yIiIicumzOLxRTcIb1WRf6jnmrTvI0m3H2Hr4LFsXJtDQtwpPdGvIfW3rYHe1mR33D5W41PTu3ZtTp04xYcIEUlJSCA0NJTY2tvBC3yNHjmC1/ucAUOfOnfn4448ZP348zz//PE2aNGHp0qW0bNmycMyYMWPIyspi8ODBpKen07VrV2JjY7HbK8/S6iIiIuVBUz8vXnswhNGRzXh/4yE+ij/MwbQsnv9yB2+s3Ev/8Ab061Sf6lXczI76O1omQURERP7Q+Zx8Fm05wnsbkjmecWlSXA9XGw+1r8ugrg2pV7N0p1nR2k/FUKkRERG5enkFDpbvOMGctQfZdeLSNClWC9zRMoDB3RsSElStVB5XpaYYKjUiIiLXzjAMNv5ymjnrDrB+f1rh7R2DazCke0N6Nq+NxXL9Ziouyee3c6xLLiIiImXCYrHQtYkvXZv4svtEJvPWHeTr7cfZknyG89n59Gxe27xsOlIjIiIi1+J4+kXmbzpEm6Bq3NEq4Lret47UiIiISJkJrObB83feYHaMirFMgoiIiMifUakRERERp6BSIyIiIk5BpUZEREScgkqNiIiIOAWVGhEREXEKKjUiIiLiFFRqRERExCmo1IiIiIhTUKkRERERp6BSIyIiIk5BpUZEREScgkqNiIiIOIVKs0q3YRjApSXMRUREpGL47XP7t8/xy6k0pebcuXMABAUFmZxERERESurcuXP4+PhcdozFuJLq4wQcDgfHjx/Hy8sLi8VyXe87MzOToKAgjh49ire393W9byk5vR7li16P8kWvR/mj1+TyDMPg3LlzBAYGYrVe/qqZSnOkxmq1Urdu3VJ9DG9vb/0PshzR61G+6PUoX/R6lD96Tf7Ynx2h+Y0uFBYRERGnoFIjIiIiTkGl5jpwd3dn4sSJuLu7mx1F0OtR3uj1KF/0epQ/ek2un0pzobCIiIg4Nx2pEREREaegUiMiIiJOQaVGREREnIJKjYiIiDgFlZo/sW7dOu6++24CAwOxWCwsXbr0T/dZs2YNbdu2xd3dncaNGzN//vxSz1mZlPQ1WbNmDRaL5XdbSkpK2QR2YlOmTKFDhw54eXlRu3ZtoqKi2Lt375/ut2TJEpo3b47dbqdVq1YsX768DNI6v6t5PebPn/+794bdbi+jxM5v9uzZtG7dunBivfDwcFasWHHZffT+uHoqNX8iKyuLkJAQZs2adUXjk5OT6dWrFzfffDPbtm1j5MiRPPHEE3z77belnLTyKOlr8pu9e/dy4sSJwq127dqllLDyWLt2Lc888wybN29m5cqV5OXlcdttt5GVlfWH+2zatIm+ffsyaNAgkpKSiIqKIioqip07d5Zhcud0Na8HXJrJ9r/fG4cPHy6jxM6vbt26vPrqqyQkJLB161Z69uzJvffey88//1zseL0/rpEhVwwwvvzyy8uOGTNmjNGiRYsit/Xu3duIjIwsxWSV15W8Jj/88IMBGGfPni2TTJXZyZMnDcBYu3btH4556KGHjF69ehW5LSwszBgyZEhpx6t0ruT1eP/99w0fH5+yCyVG9erVjXfeeafY3+n9cW10pOY6i4uLIyIioshtkZGRxMXFmZRIfhMaGkpAQAC33norGzduNDuOU8rIyACgRo0afzhG75GycyWvB8D58+epX78+QUFBlz2KINemoKCARYsWkZWVRXh4eLFj9P64Nio111lKSgp+fn5FbvPz8yMzM5OLFy+alKpyCwgIICYmhs8//5zPP/+coKAgbrrpJhITE82O5lQcDgcjR46kS5cutGzZ8g/H/dF7RNc4XV9X+no0a9aM9957j6+++ooPP/wQh8NB586d+fXXX8swrXPbsWMHVatWxd3dnaeeeoovv/ySG2+8sdixen9cm0qzSrdUXs2aNaNZs2aF/+7cuTMHDhzgzTffZOHChSYmcy7PPPMMO3fuZMOGDWZHEa789QgPDy9y1KBz587ccMMNzJkzh5dffrm0Y1YKzZo1Y9u2bWRkZPDZZ58xYMAA1q5d+4fFRq6ejtRcZ/7+/qSmpha5LTU1FW9vbzw8PExKJf+rY8eO/PLLL2bHcBrDhg1j2bJl/PDDD9StW/eyY//oPeLv71+aESuVkrwe/8vV1ZU2bdro/XEdubm50bhxY9q1a8eUKVMICQlh+vTpxY7V++PaqNRcZ+Hh4axatarIbStXrvzD86dijm3bthEQEGB2jArPMAyGDRvGl19+yerVqwkODv7TffQeKT1X83r8r4KCAnbs2KH3RylyOBzk5OQU+zu9P66R2Vcql3fnzp0zkpKSjKSkJAMw3njjDSMpKck4fPiwYRiGMXbsWKNfv36F4w8ePGh4enoaf/3rX43du3cbs2bNMmw2mxEbG2vWU3A6JX1N3nzzTWPp0qXG/v37jR07dhgjRowwrFar8f3335v1FJzG0KFDDR8fH2PNmjXGiRMnCrcLFy4UjunXr58xduzYwn9v3LjRcHFxMV5//XVj9+7dxsSJEw1XV1djx44dZjwFp3I1r8ff//5349tvvzUOHDhgJCQkGH369DHsdrvx888/m/EUnM7YsWONtWvXGsnJycZPP/1kjB071rBYLMZ3331nGIbeH9ebSs2f+O3rwP+7DRgwwDAMwxgwYIDRo0eP3+0TGhpquLm5GQ0bNjTef//9Ms/tzEr6mvzjH/8wGjVqZNjtdqNGjRrGTTfdZKxevdqc8E6muNcBKPK/+R49ehS+Nr/59NNPjaZNmxpubm5GixYtjG+++aZsgzupq3k9Ro4cadSrV89wc3Mz/Pz8jDvvvNNITEws+/BO6vHHHzfq169vuLm5GbVq1TJuueWWwkJjGHp/XG8WwzCMsjsuJCIiIlI6dE2NiIiIOAWVGhEREXEKKjUiIiLiFFRqRERExCmo1IiIiIhTUKkRERERp6BSIyIiIk5BpUZEREScgkqNiIiIOAWVGhEREXEKKjUiIiLiFFRqRERExCn8Pzmh+bHWjXY3AAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from cfdmod.analysis.inflow.functions import calculate_turbulence_intensity\n", "\n", "turbulence_intersity_components = [\"ux\"]\n", "turbulence_intersity_data = calculate_turbulence_intensity(\n", " inflow_data=profile, for_components=turbulence_intersity_components\n", ")\n", "\n", "for component in turbulence_intersity_components:\n", " x = profile.points[\"z\"]\n", " y = turbulence_intersity_data[f\"I_{component}\"]\n", " plt.plot(x, y, label=f\"{component}\")\n", "\n", "plt.legend()\n", "plt.show()\n", "plt.close()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Spectral density" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Set normalization parameters" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "NormalizationParameters(reference_velocity=1.0, characteristic_length=1.0)" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from cfdmod.analysis.inflow.profile import NormalizationParameters\n", "\n", "params = NormalizationParameters(reference_velocity=1.0, characteristic_length=1.0)\n", "params" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAGdCAYAAADAAnMpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAABfa0lEQVR4nO3dd3hU1dbH8W96gRQghRACoXcSagQLqJGigFgoNrjYuxivSpQioGLFqKAoF7teEUWkI1VBIlwCkd4h1CSEkpCEtJnz/nEw7+UKmDIl5fd5njzP7GHmrLU9mFmcOXttF8MwDEREREScxNXZCYiIiEj1pmJEREREnErFiIiIiDiVihERERFxKhUjIiIi4lQqRkRERMSpVIyIiIiIU6kYEREREadyd3YCJWG1Wjl27Bh+fn64uLg4Ox0REREpAcMwOHv2LPXq1cPV9dLXPypFMXLs2DEiIiKcnYaIiIiUweHDh6lfv/4l/7xSFCN+fn6AORl/f38nZyMiIiIlkZWVRURERPHn+KVUimLkz69m/P39VYyIiIhUMn93i4VuYBURERGnUjEiIiIiTqViRERERJyqUtwzUhIWi4XCwkJnp+E0bm5uuLu7a+mziIhUOlWiGMnOzubIkSMYhuHsVJzK19eXsLAwPD09nZ2KiIhIiVX6YsRisXDkyBF8fX0JDg6ullcGDMOgoKCAEydOcODAAZo1a3bZ5jIiIiIVSaUvRgoLCzEMg+DgYHx8fJydjtP4+Pjg4eFBSkoKBQUFeHt7OzslERGREqky/3yujldE/peuhoiISGWkTy8RERFxqlIXI7/++iv9+/enXr16uLi4MGfOnL99z6pVq+jYsSNeXl40bdqUzz77rAypioiISFVU6mIkJyeHqKgopk6dWqLXHzhwgJtuuolrr72W5ORkRo4cyf3338+SJUtKnayIiIhUPaW+gbVv37707du3xK+fNm0ajRo14u233wagVatWrFmzhnfeeYfevXuXNryIiIhUMXa/ZyQxMZHY2NgLnuvduzeJiYmXfE9+fj5ZWVkX/IiIiFQJmUfhs36QscfZmQAwa8NhRv2wmYIiq9NysHsxkpqaSmho6AXPhYaGkpWVxblz5y76nkmTJhEQEFD8ExERUeJ4hmGQW1DklJ/SNF2LjIwkISHhgueio6N56aWXWLVqFZ6enqxevbr4z9544w1CQkJIS0srcQwREalgivLhu2FwcDXMG+nsbNh46DQv/riVb/9zmNkbjzgtjwrZZyQ+Pp64uLjicVZWVokLknOFFlqPdc79KNsn9MbXs/z/SXv27MnIkSO55557+OOPP9i/fz9jxoxh1qxZfynsRESkElkcD0c3gHcg3DzFqamkZeXx8JdJFFis9G4TyuDOJf+Hv63ZvRipW7fuX/41n5aWhr+//yWblHl5eeHl5WXv1Cq0l19+maVLl/Lggw+ydetWhg8fzoABA5ydloiIlFXyN7BhBuACt/0LajdyWip5hRYe/DKJ9LP5tAj14+3B0bi6Oq9fl92LkW7durFw4cILnlu6dCndunWzSzwfDze2T3DOjbE+Hm42O5anpydff/017du3p2HDhrzzzjs2O7aIiDjY8c0w/2nzcc9R0OwGp6ViGAYvzN7CH4fPEOjrwfRhnanp5dwvSkodPTs7m7179xaPDxw4QHJyMrVr16ZBgwbEx8dz9OhRvvjiCwAefvhhpkyZwnPPPce9997LihUr+O6771iwYIHtZvFfXFxcbPJVib25urr+5R6T/911eO3atQCcOnWKU6dOUaNGDYflJyIiNnLuNMy8G4ryoFkvuOY5p6YzY80BZm86ipurC1Pv7EiDOr5OzQfKcAPrhg0b6NChAx06dAAgLi6ODh06MHbsWACOHz/OoUOHil/fqFEjFixYwNKlS4mKiuLtt9/mX//6V7Vf1hscHMzx48eLx1lZWRw4cKB4vG/fPp5++mmmT59OTEwMw4cPx2p13p3OIiJSBlYrzH4QzqRArUi49WNw4tYdv+4+wasLdwAw+qZWXNk0yGm5/LdSX0Lo2bPnZVeNXKy7as+ePdm0aVNpQ1Vp1113HZ999hn9+/cnMDCQsWPH4uZmfs1jsVi4++676d27NyNGjKBPnz60a9eOt99+m2effdbJmYuISIn9+gbs+RncvWHwl+BTy2mpHMzI4fFvNmI1YHDn+vyje6TTcvlfFf/7jCoqPj6eAwcO0K9fPwICApg4cWLxlZFXXnmFlJQU5s+fD0BYWBgff/wxd9xxB7169SIqKsqZqYuISEns/hlWvWY+7pcAYe2dlsrZvELu/2IDWXlFdGgQyMSBbSvUBrMuRmmaYzhJVlYWAQEBZGZm4u/vf8Gf5eXlceDAARo1aoS3t7eTMqwY9N9CRKSCOHUAPu4BeZnQ+T7oN9lpqVitBg9+uYFlO9Kp6+/N3MevJMTfMZ8Rl/v8/m/atVdERMSWCs/Bd/eYhUh4Z+gzyanpTF66m2U70vF0d+Wjezo5rBApDRUjIiIitmIYMD8OUreAbxAM/gLcndc3a/7mY0xZaa6Aff22dkRFBDotl8tRMSIiImIrSZ/CH9+AiysM+hQCwp2WyrZjmTw7azMAD17TmFs61HdaLn9HxYiIiIgtHNkAC8/3EIl9CRpd47RUTmbn8+AXSZwrtHB1syCe79PSabmUhIoRERGR8so+YW6AZy2EVv2h+5NOS6XQYuWRrzdy9Mw5GgXVYModHXFzYqv3klAxIiIiUh6WIvh+BGQdhTrN4OYPwInLZsfP28b6A6eo6eXO9GGdCPD1cFouJaViREREpDxWTISDq8GjBgz5CrwvvYTV3r5el8JXvx/CxQXeHRpN0xA/p+VSGipGREREymr7XPgtwXx88xQIcd69GesPnGLcT9sA+GevFlzfKtRpuZSWihEREZGyyNgDcx41H3d7HNre6rRUjp45xyNfJVFkNejXPoxHezZxWi5loWJERESktPKzzZ14C85Cw6sgdrzTUjlXYOHBLzZwMqeA1mH+vHF7+wrV6r0kVIyIiIiUhmHA3CfgxE7wCzP7ibg5Z6s3wzB49vs/2HYsizo1PJk+vDO+npVv2zkVI04SGRlJQkLCBc9FR0fz0ksvce+999KvX78L/qywsJCQkBBmzJjhwCxFROQvfv8Ats0GV3cY9DnUDHFaKh/+so/5m4/j7urCh3d3IjzQx2m5lEflK5/+jmFAYa5zYnv42mQ51/33388111zD8ePHCQsLA2D+/Pnk5uYyZMiQch9fRETK6OBv8PMY83HvSdAgxmmprNiZxptLdgEw/uY2dG1U22m5lFfVK0YKc+HVes6J/cIx8KxR7sN0796dFi1a8OWXX/Lcc2Y3v08//ZRBgwZRs2bNch9fRETKIOs4zPoHGBZoNxi6PuC0VPamn+WpfydjGHBXTAPuimnotFxsQV/TVFD3338/n376KQBpaWksWrSIe++918lZiYhUU0UFMGs45KRDSBvon+C0xmaZ5wp54IskzuYX0TWyNuP6t3FKHrZU9a6MePiaVyicFbuEXF1dMQzjgucKCwuLHw8bNoxRo0aRmJjI2rVradSoEVdffbXNUhURkVJYOgYOrwOvABjypU2ugpeFxWrw5L83cSAjh/BAHz64uyOe7pX/ukLVK0ZcXJz2l6Q0goODOX78ePE4KyuLAwcOFI/r1KnDwIED+fTTT0lMTGTEiBHOSFNERDbPgnXTzMe3fgR1nNfD443FO/ll9wm8PVz5eFgngmp6OS0XW6p6xUglcd111/HZZ5/Rv39/AgMDGTt2LG5ubhe85v7776dfv35YLBaGDx/upExFRKqxtG0w7/ymd9c8Cy36Oi2VOZuO8tGv+wF4a1AUbeoFOC0XW1Mx4iTx8fEcOHCAfv36ERAQwMSJEy+4MgIQGxtLWFgYbdq0oV49J92UKyJSXZ07YzY2K8yFJtdBz3inpbL5yBme/2EzAI9d24R+7avWZ4KKESfx9/fn22+/veC5/736kZOTw+nTp7nvvvscmZqIiFit8OPDcGo/BDSA22aAq9vfv88O0s/m8eAXSeQXWbm+ZQjP3NDCKXnYk4qRCshqtZKRkcHbb79NYGAgAwYMcHZKIiLVy5rJsHsRuHnB4M/B1zk9PPKLLDz8ZRKpWXk0DalJwtBoXF0rV6v3klAxUgEdOnSIRo0aUb9+fT777DPc3XWaREQcZu9yWPGy+fimtyC8o1PSMAyDsXO2sfHQGfy93Zk+rDN+3h5OycXe9ClXAUVGRv5l2a+IiDjAmUPww/2AAR2HmT9O8kViCjM3HMbVBd6/syONgir+StGyqvyLk0VERGyhMA++GwbnTkG9DtD3TaelsnZvBhPmbwcgvm8rejQPdloujqBiREREBGDRs3BsE/jUhsFfgIe3U9I4fCqXR7/ZiMVqcGuHcO6/upFT8nCkKlOM6GsN/TcQESmzpM9h4xeAC9w+AwIbOCWNnPwiHvhiA2dyC4mqH8Crt7bDxUlt5x2p0hcjfzYKKygocHImzpeba+5W7OFRNW9wEhGxi6MbYeGz5uPrRps9RZzAajWI+y6ZnalnCfbz4qN7OuPt4ZzlxI5W6W9gdXd3x9fXlxMnTuDh4YGra6Wvr0rNMAxyc3NJT08nMDDwL51cRUTkEnJOmveJWPKhxY1wVZzTUnl/xV6WbEvD082VaXd3om6Ac74mcoZKX4y4uLgQFhbGgQMHSElJcXY6ThUYGEjdunWdnYaISOVgtcDs+yHzMNRuDAM/BCf9g3bx1lTeWbYbgJdvaUunhrWckoezVPpiBMDT05NmzZpV669qPDw8dEVERKQ0Vk2CfSvA3QeGfAU+gU5JY2dqFnHfJQMw4spIBneOcEoezlQlihEAV1dXvL2rzyUtEREph12L4NfzS3cHvAehbZySxumcAh74YgO5BRa6N6nDize2ckoezlb9brAQEZHq7eQ+mP2Q+bjrQ9B+sFPSKLJYeeybjRw+dY4GtX2ZemdH3N2q58dy9Zy1iIhUTwU5MPMeyM+EiCug18tOS+WVhTtYu+8kvp5uTB/WmVo1PJ2Wi7OpGBERkerBMGDeU5C+DWqEwKDPwN05BcB3Gw7z6W8HAZg8OJoWdf2ckkdFoWJERESqh/XTYcsscHEzCxH/MKekkZRymtE/bgXg6djm9GmrVZAqRkREpOo7tA6WxJuPe02EyCudkkZqZh4Pf5VEgcVKnzZ1eeK6pk7Jo6JRMSIiIlXb2TSYNRysRdDmFrjiUaekkVdo4cEvN3DibD4t6/rx9uAoXF2rfqv3klAxIiIiVZelEL4fAWePQ1ALGDAFnLDXi2EYxM/ewuYjmQT6ejB9WGdqeFWZ7hrlpmJERESqrmUvQcpv4OkHQ78Gr5pOSeNfqw/w46ajuLm68MGdHYmo7euUPCoqFSMiIlI1bZ0NiVPMx7d8CEHNnJLGL7tPMGnRDgDG9mtN96ZBTsmjIlMxIiIiVU/6TvjpcfPxlSOhVX+npHEgI4cnvtmI1YAhnSMY1q2hU/Ko6FSMiIhI1ZKXBTPvhsIcaHQNXDfGKWmczSvk/s//Q1ZeEZ0a1mLCwDa4OOF+lcpAxYiIiFQdhgE/PQon94B/ONz2Cbg5/kZRq9Vg5LfJ7DuRQ11/bz68uyNe7trM9FLKVIxMnTqVyMhIvL29iYmJYf369Zd8bWFhIRMmTKBJkyZ4e3sTFRXF4sWLy5ywiIjIJa19D3bMA1cPGPwF1Ax2ShpvL93F8p3peLm78vGwToT4aSPXyyl1MTJz5kzi4uIYN24cGzduJCoqit69e5Oenn7R148ePZqPPvqI999/n+3bt/Pwww9zyy23sGnTpnInLyIiUmz/L+bqGYC+r0P9zk5JY94fx5i6ch8Ar9/Wnvb1A52SR2XiYhiGUZo3xMTE0KVLF6ZMMe9QtlqtRERE8MQTTzBq1Ki/vL5evXq8+OKLPPbYY8XP3Xbbbfj4+PDVV1+VKGZWVhYBAQFkZmbi7+9fmnRFRKQ6yDwCH/WA3AyIuhMGfuCUfiJbj2Zy+7S15BVaeeiaxsTf2MrhOVQkJf38LtWVkYKCApKSkoiNjf3/A7i6EhsbS2Ji4kXfk5+fj7f3hZenfHx8WLNmzSXj5Ofnk5WVdcGPiIjIRRXlw3fDzUKkbjvoN9kphUhGdj4PfZlEXqGVHs2Dea5PS4fnUFmVqhjJyMjAYrEQGhp6wfOhoaGkpqZe9D29e/dm8uTJ7NmzB6vVytKlS5k9ezbHjx+/ZJxJkyYREBBQ/BMREVGaNEVEpDpZHA9HN4B3IAz+Ejx8HJ5CQZGVR7/ayNEz52gUVIP37uiAm1q9l5jdV9O8++67NGvWjJYtW+Lp6cnjjz/OiBEjcHW9dOj4+HgyMzOLfw4fPmzvNEVEpDJK/gY2zABc4LZ/Qe1GTklj/LxtrD94Cj8vd6YP60yAj4dT8qisSlWMBAUF4ebmRlpa2gXPp6WlUbfuxbdADg4OZs6cOeTk5JCSksLOnTupWbMmjRs3vmQcLy8v/P39L/gRERG5wPHNMP9p83HPUdDsBqek8dXvKXy97hAuLvDuHdE0DXFOy/nKrFTFiKenJ506dWL58uXFz1mtVpYvX063bt0u+15vb2/Cw8MpKirihx9+4Oabby5bxiIiIudOm43NivKgWS+45jmnpLFu/0lemrsNgGd7t+C6lqF/8w65mFJ3gomLi2P48OF07tyZrl27kpCQQE5ODiNGjABg2LBhhIeHM2nSJADWrVvH0aNHiY6O5ujRo7z00ktYrVaee845f3FERKSSs1ph9oNwJgUCG8ItH8Flvvq3lyOnc3n0640UWQ36R9XjkR5NHJ5DVVHqYmTIkCGcOHGCsWPHkpqaSnR0NIsXLy6+qfXQoUMX3A+Sl5fH6NGj2b9/PzVr1uTGG2/kyy+/JDAw0GaTEBGRauTXN2DPz+DuDUO+BN/aDk8ht6CIB79I4mROAW3q+fPGbe3V6r0cSt1nxBnUZ0RERADY/TN8MxgwYOCHEH2nw1MwDIPH/72JBZuPE1TTk58ev4rwQMev4KkM7NJnRERExGlOHYDZ9wMGdL7PKYUIwAer9rFg83E83Fz48O5OKkRsQMWIiIhUfIXn4Lt7IC8TwjtDn0lOSWPZ9jTe+nkXAOMHtKVLpOO/IqqKVIyIiEjFZhgwPw5St4BvkLkBnruXw9PYm36WkTOTMQy454qG3BnTwOE5VFUqRkREpGJL+hT++AZcXOH2TyAg3OEpZOYWcv/nG8jOL6Jro9qM7d/a4TlUZSpGRESk4jqyARaebwVx/Tho3MPhKRRZrDz+740cPJlLeKAPH97VEQ83fXzakv5riohIxZR9Ar4bBtZCaNkPrnzKKWm8sWQXq/dk4OPhxsfDOlGnpuO/IqrqVIyIiEjFYymC70dA1lGo09RcxuuEPh6zNx7h41/3A/DWoCja1AtweA7VgYoRERGpeFZMhIOrwaMGDPkavB3fY+qPw2cYNXsLAE9c15Sb2oc5PIfqQsWIiIhULNvnwm8J5uObp0BIS4enkJ6Vx4NfbqCgyEpsqxCejm3u8ByqExUjIiJScWTsgTmPmo+7PQ5tb3V4CvlFFh76Kom0rHyahdTknSHRuLqq1bs9qRgREZGKIT/b3Im34Cw0vBJiX3J4CoZhMPrHrWw6dAZ/b3emD+uMn7eHw/OoblSMiIiI8xkGzH0cTuyEmnXh9k/BzfFFwGdrDzIr6QiuLjDlzo5EBtVweA7VkYoRERFxvt8/gG0/gqs7DP4c/EIdnsJvezN4ecEOAF64sRXXNA92eA7VlYoRERFxroO/wc9jzMe9X4UGVzg8hUMnc3nsm41YrAa3dgznvqsaOTyH6kzFiIiIOE/WcZj1DzAs0G4QdH3Q4Slk5xdx/xf/4UxuIVERgbx6SztcnNDTpDpTMSIiIs5RVACzhkNOOoS0gf7vOryxmdVq8Mx3yexOyybEz4uP7+mEt4ebQ3MQFSMiIuIsS8fA4XXgFQBDvgRPx98s+u7yPSzZloanmyvT7ulEqL+3w3MQFSMiIuIMm2fBumnm41umQZ0mDk9h8dbjvLt8DwCv3tqOjg1qOTwHMakYERERx0rbBvOeNB9f/U9oeaPDU9iZmkXcd38AcO+Vjbi9U32H5yD/T8WIiIg4zrkzZmOzwlxofC1c+4LDUziVU8D9n28gt8DCVU2DeOFGx7eblwupGBEREcewWuHHh+HUfgiIgNtmgKtjbxYttFh57OuNHDl9jga1fZlyZwfc3fRR6Gw6AyIi4hhrJsPuReDmBYO/gBp1HJ7CKwt2kLj/JDU83fjX8M4E+no6PAf5KxUjIiJif3uXw4qXzcc3vQXhHR2ewsz/HOKztQcBeGdINM1D/Ryeg1ycihEREbGvM4fgh/sBAzoOM38cLCnlFKPnbAUg7obm9GpT1+E5yKWpGBEREfspzIPvhsG5UxAWDX3fdHgKxzPP8dCXGym0GPRtW5fHr23q8Bzk8lSMiIiI/Sx6Fo5tAp9aZmMzD8c2FcsrtPDQl0lkZOfTsq4fbw2KwtVVrd4rGhUjIiJiH0mfw8YvABdz5UxgA4eGNwyDUT9sZvORTGr5ejB9WGdqeLk7NAcpGRUjIiJie0c3wsJnzcfXvQhNr3d4CtNX72dO8jHcXF344K5ORNT2dXgOUjIqRkRExLZyTpr3iVjyoXlfuOoZh6ewalc6ry3aCcC4/q3p1sTxy4il5FSMiIiI7VgtMPt+yDwMtRub+864OvajZv+JbJ749yasBgztEsE9VzR0aHwpPRUjIiJiO6smwb4V4O4DQ74Cn0CHhs/KK+SBLzZwNq+Izg1rMeHmtri46IbVik7FiIiI2MbOhfDr+aW7A96D0DYODW+xGoz8Npl9J3IIC/Dmw7s74emuj7nKQGdJRETK7+Q++PEh83HXh6D9YIen8PbPu1ixMx0vd1c+vqczwX5eDs9BykbFiIiIlE9BDsy8B/KzICIGer3s8BTm/nGMD1btA+CN29vTrn6Aw3OQslMxIiIiZWcYMO8pSN8GNUJg0Ofg7tjN57YezeS57/8A4OEeTbg5Otyh8aX8VIyIiEjZrZ8OW2aBixsM+hT8wxwaPiM7nwe/2EBeoZWeLYJ5tncLh8YX21AxIiIiZXNoHSyJNx/fMAEir3Jo+IIiK498lcSxzDwaB9fg3aEdcFOr90pJxYiIiJTe2TSYNRysRdDmFuj2mMNTeGneNv5z8DR+Xu5MH9aZAB8Ph+cgtqFiRERESsdSCN+PgLPHIagFDJgCDu7l8eXvKXyz7hAuLvDeHR1oElzTofHFtlSMiIhI6Sx7CVJ+A08/s7GZl2MLgd/3n2T83G0APN+nJde2DHFofLE9FSMiIlJyW2dD4hTz8cAPILi5Q8MfPpXLo19vpMhqMCCqHg9d09ih8cU+VIyIiEjJpO+Enx43H1/5FLQe4NDwuQVFPPhlEqdyCmgb7s/rt7VXq/cqQsWIiIj8vbwsmHk3FOZA5NVw3ViHhjcMg2dnbWbH8SyCanry8T2d8fF0c2gOYj8qRkRE5PIMA356FE7uAb96cPun4Obu0BSmrtzLgi3H8XBzYdrdnagX6OPQ+GJfKkZEROTy1r4HO+aBqwcM/gJqBjs0/NLtabz1824AJtzcls6RtR0aX+yvTMXI1KlTiYyMxNvbm5iYGNavX3/Z1yckJNCiRQt8fHyIiIjg6aefJi8vr0wJi4iIA+3/xVw9A9D3dYjo4tDwe9LO8vTMZACGdWvIHV0bODS+OEapi5GZM2cSFxfHuHHj2LhxI1FRUfTu3Zv09PSLvv6bb75h1KhRjBs3jh07djBjxgxmzpzJCy+8UO7kRUTEjjKPwPf3gmGFqDuh870ODX8mt4D7v9hAdn4RVzSuzZh+rR0aXxyn1MXI5MmTeeCBBxgxYgStW7dm2rRp+Pr68sknn1z09WvXruXKK6/kzjvvJDIykl69enHHHXf87dUUERFxoqJ8+G445GZA3XbQb7JDG5sVWaw88e9NpJzMJTzQhw/u6oSHm+4sqKpKdWYLCgpISkoiNjb2/w/g6kpsbCyJiYkXfU/37t1JSkoqLj7279/PwoULufHGGy8ZJz8/n6ysrAt+RETEQaxWWPgsHN0A3gEw+EvwcNwNo1arwcsLdrB6TwY+Hm5MH9aZ2jUcuxOwOFapbofOyMjAYrEQGhp6wfOhoaHs3Lnzou+58847ycjI4KqrrsIwDIqKinj44Ycv+zXNpEmTGD9+fGlSExERW8hOh9kPwv6V5vjWf0HtRg4Ln5qZxzOzkvlt70kAJg+OonU9f4fFF+ew+zWvVatW8eqrr/LBBx+wceNGZs+ezYIFC5g4ceIl3xMfH09mZmbxz+HDh+2dpoiI7FsBH15pFiLuPjBwGjTv5bDwi7Ycp8+7v/Lb3pN4e7jyxm3t6dsuzGHxxXlKdWUkKCgINzc30tLSLng+LS2NunXrXvQ9Y8aM4Z577uH+++8HoF27duTk5PDggw/y4osv4ur613rIy8sLLy+v0qQmIiJlZSmCla/AmncAA0Jam71EQlo6JHxOfhHj523juw1HAGgXHkDC0GhtfleNlOrKiKenJ506dWL58uXFz1mtVpYvX063bt0u+p7c3Ny/FBxubmbXPMMwSpuviIjY0pnD8NmNsGYyYECnf8ADKxxWiGw6dJob31vNdxuO4OICj/Zswg+PdFchUs2UuoVeXFwcw4cPp3PnznTt2pWEhARycnIYMWIEAMOGDSM8PJxJkyYB0L9/fyZPnkyHDh2IiYlh7969jBkzhv79+xcXJSIi4gQ75sNPj0HeGfDyh/7vQttbHRK6yGJl6sp9vLdiDxarQXigD5MHRxHTuI5D4kvFUupiZMiQIZw4cYKxY8eSmppKdHQ0ixcvLr6p9dChQxdcCRk9ejQuLi6MHj2ao0ePEhwcTP/+/XnllVdsNwsRESm5onz4eQys/8gc1+sIt3/isBtVD5/KZeTMZJJSTgMwIKoeEwe2JcDHwyHxpeJxMSrBdyVZWVkEBASQmZmJv7/uqhYRKbOMvfD9CEjdbI67PQ7XjwN3+y+dNQyD2RuPMm7uNrLzi/DzcmfiwLYM7BBu99jiHCX9/HbsTkciIuI8f8yEBXFQkA0+teGWadC8t0NCZ+YW8uKcLczffByALpG1mDw4mojavg6JLxWbihERkaquIMdsYpb8tTlueBXcNh386zkkfOK+kzzzXTLHMvNwc3Xh6dhmPNKzKW6ujuvoKhWbihERkaosdav5tUzGbnBxhWuegx7Pgav9FxAUFFmZvHQ3H/26D8OAyDq+JAztQHREoN1jS+WiYkREpCoyDNgwAxa/AJZ88AuDW6dDo6sdEn5vejYjZ25i61FzO4+hXSIY0681Nbz0sSN/pb8VIiJVzbkzMPcJ2DHXHDfrBQM/hBpBdg9tGAZfrzvEywu2k1doJdDXg9dubU+fthdvjCkCKkZERKqWIxvMr2XOHAJXD4h9Ca54FC7S7drWMrLzGfXDZpbtSAfg6mZBvDUoilB/b7vHlspNxYiISFVgtULi+7B8AliLILAhDPoUwjs5JPzKXek8O2szGdn5eLq58nzflozoHomrblKVElAxIiJS2WWfgDkPw95l5rjNLWY3Ve8Au4fOK7QwaeEOPk9MAaBFqB8JQ6NpFaaeUFJyKkZERCqz/b/A7AchOxXcvaHPa+b+Mi72vyKx/VgWT327iT3p2QCMuDKS5/u0xNtDW31I6agYERGpjCxF8Mvr8OubgAFBLcyvZULb2D201WowY80B3lyyiwKLlWA/L94aFEWP5sF2jy1Vk4oREZHKJvMo/HA/HFprjjvcA31fB88adg+dmpnHM7OS+W3vSQBiW4Xy+m3tqFPTy+6xpepSMSIiUpnsWgxzHoFzp8CzpnlvSLvbHRJ60ZbjxP+4hTO5hfh4uDGmX2vu6BqBiwO+EpKqTcWIiEhlUFQAy8bB7x+Y47AouP1TqNPE7qFz8osYP28b3204AkC78AAShkbTJLim3WNL9aBiRESkoju1H2aNgOPJ5viKR83+Ie72/2pk06HTjJyZTMrJXFxc4JEeTRgZ2xxPd/v3LZHqQ8WIiEhFtuV7mDcSCs6CTy2zk2qLvnYPW2SxMnXlPt5bsQeL1SA80IfJg6OIaVzH7rGl+lExIiJSERXkwuLnYeMX5rhBN7jtXxBQ3+6hD5/KZeTMZJJSTgMwIKoeEwe2JcDHw+6xpXpSMSIiUtGkbTdbup/YCbjANf+EHqPAzb6/sg3DYPbGo4ybu43s/CL8vNyZOLAtAzuE2zWuiIoREZGKwjBg4+ew6HkoyoOaoeZOu4172D10Zm4hL87ZwvzNxwHoElmLyYOjiajta/fYIipGREQqgrxM896QbbPNcZPr4ZaPoKb9G4kl7jvJM98lcywzDzdXF56ObcYjPZvipn1lxEFUjIiIONvRJPj+Xjh9EFzd4box0P1Ju++0W1BkZfLS3Xz06z4MAyLr+JIwtAPREYF2jSvyv1SMiIg4i2FA4lRY9hJYCyGgAdz+CUR0sXvovenZjJy5ia1HswAY2iWCMf1aU8NLHwviePpbJyLiDDknzU6qe5aY41b9YcD75vJdOzIMg6/XHeLlBdvJK7QS6OvBa7e2p0/bunaNK3I5KkZERBzt4G/m3jJnj4GbF/R5FTrfZ/eddjOy8xn1w2aW7UgH4OpmQbw1KIpQf2+7xhX5OypGREQcxWqBX9+CX14Dwwp1mpk77dZtZ/fQK3el8+yszWRk5+Pp5srzfVsyonskrrpJVSoAFSMiIo6QdRxmPwAHV5vj6Lug7xvgZd/9XfIKLUxauIPPE1MAaBHqR8LQaFqF+ds1rkhpqBgREbG3PUvhx4cg9yR41IB+70DUELuH3X4si6e+3cSe9GwARlwZyfN9WuLt4Wb32CKloWJERMReigpgxQRY+745rtsObv8MgpraNazVajBjzQHeXLKLAouVYD8v3hoURY/m9u9ZIlIWKkZEROzh9EGzd8jRJHPc9UG4YSJ42Pdm0dTMPJ6Zlcxve08CENsqlNdva0edmvbf4VekrFSMiIjY2rYfYe6TkJ8F3gFw81Rz6a6dLdpynPgft3AmtxAfDzfG9GvNHV0jcLHzKh2R8lIxIiJiK4XnYHE8JH1qjiNizJ12AxvYNWxOfhHj523juw1HAGgXHkDC0GiaBNv35lgRW1ExIiJiCyd2wawRkL4NcIGrnoZrXwA3D7uG3XToNCNnJpNyMhcXF3ikRxNGxjbH092+reRFbEnFiIhIeRgGbPoKFj0HhblQIxhu/RiaXGfXsEUWKx+s2se7y/dgsRrUC/Bm8pBormhcx65xRexBxYiISFnln4X5T8OWWea4cU+45WPwC7Vr2MOnchk5M5mklNMADIiqx8SBbQnwse9VGBF7UTEiIlIWx5Lh+xFwaj+4uJlfyVwVZ9eddg3DYPbGo4ybu43s/CL8vNyZOLAtAzuE2y2miCOoGBERKQ3DgHUfwdIxYCkA//pw+wxocIVdw2bmFvLinC3M33wcgC6RtZg8OJqI2r52jSviCCpGRERKKvcU/PQ47Fpgjlv2M3fa9a1t17CJ+07yzHfJHMvMw83Vhadjm/FIz6a4aV8ZqSJUjIiIlERKornTbtYRcPOEXq9A1wfsutNuQZGVyUt389Gv+zAMiKzjS8LQDkRHBNotpogzqBgREbkcqwXWTIaVk8CwQO0m5k67YVF2Dbs3PZuRMzex9WgWAEO7RDCmX2tqeOnXtlQ9+lstInIpZ9PMnXYP/GKO2w+Bm94GLz+7hTQMg6/XHeLlBdvJK7QS6OvBa7e2p0/bunaLKeJsKkZERC5m73Jzp92cE+DhCze+BdF32vVrmYzsfEb9sJllO9IBuLpZEG8NiiLU37772Yg4m4oREZH/ZimEla/AmnfMcUgb82uZ4BZ2DbtyVzrPztpMRnY+nm6uPN+3JSO6R+Kqm1SlGlAxIiLypzOH4Pv74Mh6c9z5Puj9Cnj42C1kXqGFSQt38HliCgAtQv1IGBpNqzB/u8UUqWhUjIiIAGyfC3Mfh7xM8AqAAe9Bm4H2DXksi6e+3cSe9GwARlwZyfN9WuLt4WbXuCIVjYoREaneCvPg59Hwn+nmOLyz2cSsVqTdQlqtBjPWHODNJbsosFgJ9vPirUFR9GgebLeYIhWZihERqb4y9pg77aZtMcfdn4Trx9p1p93UzDyemZXMb3tPAhDbKpTXb2tHnZpedospUtGVaROFqVOnEhkZibe3NzExMaxfv/6Sr+3ZsycuLi5/+bnpppvKnLSISLkl/xs+6mEWIr514K7voddEuxYii7Ycp8+7v/Lb3pP4eLjx6i3tmD6skwoRqfZKfWVk5syZxMXFMW3aNGJiYkhISKB3797s2rWLkJCQv7x+9uzZFBQUFI9PnjxJVFQUgwYNKl/mIiJlkZ8NC/8Jf/zbHEdeDbdOB/8wu4XMyS9i/LxtfLfhCADtwgNIGBpNk+CadospUpm4GIZhlOYNMTExdOnShSlTpgBgtVqJiIjgiSeeYNSoUX/7/oSEBMaOHcvx48epUaNGiWJmZWUREBBAZmYm/v66w1xEyih1i/m1zMk94OIKPePh6mfA1X43jG46dJqRM5NJOZmLiws80qMJI2Ob4+luv919RSqKkn5+l+rKSEFBAUlJScTHxxc/5+rqSmxsLImJiSU6xowZMxg6dOhlC5H8/Hzy8/OLx1lZWaVJU0TkQoYB//kXLHkRLPngV8+8SbVhd7uFLLJY+WDVPt5dvgeL1aBegDeTh0RzReM6dospUlmVqhjJyMjAYrEQGhp6wfOhoaHs3Lnzb9+/fv16tm7dyowZMy77ukmTJjF+/PjSpCYicnHnTsPcJ2DHPHPcvA8M/NCuO+0ePpXLyJnJJKWcBmBAVD0mDmxLgI/97kcRqcwcuppmxowZtGvXjq5du172dfHx8cTFxRWPs7KyiIiIsHd6IlLVHF5vNjHLPASuHnDDBLjiEbu1dDcMg9kbjzJu7jay84vw83Jn4sC2DOwQbpd4IlVFqYqRoKAg3NzcSEtLu+D5tLQ06ta9/CZOOTk5fPvtt0yYMOFv43h5eeHlpbvLRaSMrFZY+y4sn2jutFurEdz+CYR3tFvIzNxCXpyzhfmbjwPQJbIWkwdHE1Hb124xRaqKUt1B5enpSadOnVi+fHnxc1arleXLl9OtW7fLvnfWrFnk5+dz9913ly1TEZGSyD4BX98Oy14yC5G2t8FDv9q1EEncd5K+7/7K/M3HcXN14Z+9mvPtg91UiIiUUKm/pomLi2P48OF07tyZrl27kpCQQE5ODiNGjABg2LBhhIeHM2nSpAveN2PGDAYOHEidOrp5S0TsZP8qmP0gZKeBuw/0fR06DrPb1zIFRVYmL93NR7/uwzAgso4vCUM7EB0RaJd4IlVVqYuRIUOGcOLECcaOHUtqairR0dEsXry4+KbWQ4cO4ep64QWXXbt2sWbNGn7++WfbZC0i8t8sRbBqEqx+GzAguJW5025IK7uF3JuezciZm9h61FztN7RLBGP6taaGlxpbi5RWqfuMOIP6jIjIJWUeMW9SPfy7Oe44HPq8Bp72+YrEMAy+XneIlxdsJ6/QSqCvB6/d2p4+bS9/35xIdWSXPiMiIhXKzoUw5xHIOwOefjDgXfMeETvJyM5n1A+bWbYjHYCrmwXx1qAoQv297RZTpDpQMSIilU9RPiwdC+ummeN6HczVMrUb2y3kyl3pPDtrMxnZ+Xi6ufJ835aM6B6Jq6t97kcRqU5UjIhI5XJyH3w/Ao7/YY67PQ7XjwN3T7uEyyu0MGnhDj5PTAGgRagfCUOjaRWmr4xFbEXFiIhUHptnwfyRUJANPrXhlmnQvLfdwm0/lsVT325iT3o2ACOujOT5Pi3x9rDfXjYi1ZGKERGp+ApyYNFzsOkrc9zwSnOn3QD7dDa1Wg1mrDnAm0t2UWCxEuznxVuDoujRPNgu8USqOxUjIlKxpW0zd9rN2AW4QI/nocdzdttpNzUzj2dmJfPb3pMAxLYK5fXb2lGnprpCi9iLihERqZgMA5I+hcXxUJQHNevCbf+CRlfbLeSiLceJ/3ELZ3IL8fFwY0y/1tzRNQIXOzVNExGTihERqXjyMmHuk7B9jjlueoN5f0iNILuEy8kvYvy8bXy34QgA7cIDSBgaTZPgmnaJJyIXUjEiIhXLkSRztcyZFHB1N1fKdHscXEu1lVaJbTp0mpEzk0k5mYuLCzzSowkjY5vj6W6feCLyVypGRKRisFrh96nmBnfWIghsALd/CvU72yVckcXKB6v28e7yPVisBvUCvJk8JJorGmv/LBFHUzEiIs6Xk2F2Ut1zfv+q1gOh/7vgE2iXcIdP5TJyZjJJKacBGBBVj4kD2xLg42GXeCJyeSpGRMS5DqyG2Q/A2ePg7g19JkGnEXbZadcwDGZvPMq4udvIzi/Cz8udiQPbMrCDfZYIi0jJqBgREeewWuCX1+GXNwADglqYO+2GtrFLuMzcQl6cs4X5m48D0CWyFpMHRxNR2z4b6olIyakYERHHyzoGP9wPKb+Z4w53Q983wLOGXcIl7jvJM98lcywzDzdXF56ObcYjPZvipn1lRCoEFSMi4li7l8CPD8O5U+BZE/olQPtBdglVUGRl8tLdfPTrPgwDIuv4kjC0A9ERgXaJJyJlo2JERByjqACWj4fEKeY4LMpcLVOniV3C7U3PZuTMTWw9mgXAkM4RjO3fmhpe+rUnUtHo/0oRsb9TB+D7e+HYRnMc8zDcMAHcbd9i3TAMvl53iJcXbCev0Eqgrwev3dqePm3r2jyWiNiGihERsa+tP8C8kZCfBd6BMPBDaHmjXUJlZOcz6ofNLNuRDsDVzYJ4a1AUof7edoknIrahYkRE7KMgFxaPgo2fm+OIK+D2GRBQ3y7hVu5K59lZm8nIzsfTzZXn+7ZkRPdIXHWTqkiFp2JERGwvfSfM+gec2AG4wNXPQM94cLP9r5y8QguTFu7g88QUAFqE+pEwNJpWYf42jyUi9qFiRERsxzBg05ew8DkoOgc1QuDWj6HJtXYJt/1YFk99u4k96dkAjLgykuf7tMTbw80u8UTEPlSMiIht5GXB/Kdh6/fmuMl1cMtHUDPE5qGsVoMZaw7w5pJdFFisBPt58dagKHo0D7Z5LBGxPxUjIlJ+xzbBrBFw+gC4uMH1Y6D7U3bZaTc1M49nZiXz296TAMS2CuX129pRp6btV+aIiGOoGBGRsjMM+P1DWDoWrIUQ0MC8STWiq13CLdpynPgft3AmtxAfDzfG9GvNHV0jcLHDPjYi4jgqRkSkbHJPwZxHYfcic9yqPwx4H3xq2TxUTn4R4+dt47sNRwBoFx5AwtBomgTXtHksEXE8FSMiUnopa+H7++DsMXDzgt6vQJf77bLT7qZDpxk5M5mUk7m4uMAjPZowMrY5nu62/wpIRJxDxYiIlJzVAqvfhlWTwLBCnaZmS/ew9jYPVWSx8sGqfby7fA8Wq0G9AG8mD4nmisZ1bB5LRJxLxYiIlMzZVJj9ABz41RxH3QE3vgVetv+q5PCpXEbOTCYp5TQAA6LqMXFgWwJ8PGweS0ScT8WIiPy9Pcvgx4cgNwM8asBNb0P0HTYPYxgGszceZdzcbWTnF+Hn5c7EgW0Z2CHc5rFEpOJQMSIil2YphBUT4bd3zXFoOxj0KQQ1s3mozNxCXpyzhfmbjwPQJbIWkwdHE1Hb1+axRKRiUTEiIhd3OsXcaffoBnPc5QHo9TJ42H7TucR9J3nmu2SOZebh5urC07HNeKRnU9y0r4xItaBiRET+avtP8NMTkJ8J3gEwYAq0HmDzMAVFViYv3c1Hv+7DMCCyji8JQzsQHRFo81giUnGpGBGR/1eYB0tegA0zzHH9LnDbDKjV0Oah9qZnM3LmJrYezQJgSOcIxvZvTQ0v/VoSqW70f72ImE7shu9HQNpWc3zlSLhuNLjZdgWLYRh8ve4QLy/YTl6hlUBfD167tT192ta1aRwRqTxUjIgIJH8DC56BwlyoEWxucNf0epuHycjOZ9QPm1m2Ix2Aq5sF8dagKEL9bX8fiohUHipGRKqz/LOw4J+w+Vtz3KgH3Dod/EJtHmrlrnSenbWZjOx8PN1ceb5vS0Z0j8RVN6mKVHsqRkSqq+ObYdY/4NQ+c6fda1+Aq54GVzebhskrtDBp4Q4+T0wBoEWoHwlDo2kV5m/TOCJSeakYEaluDAPWT4efXwRLAfjXh9v+BQ272TzU9mNZPPXtJvakZwMw4spInu/TEm8P2xY8IlK5qRgRqU5yT8HcJ2DnfHPc4ka4eSr41rZpGKvVYMaaA7y5ZBcFFivBfl68NSiKHs2DbRpHRKoGFSMi1cWhdfDDfZB5GNw84YaJEPOQzXfaTc3M45lZyfy29yQAsa1Cef22dtSp6WXTOCJSdagYEanqrFb47R1Y8QoYFqjd2Nxpt160zUMt2nKc+B+3cCa3EB8PN8b0a80dXSNwsXHBIyJVi4oRkaosOx1mPwj7V5rjdoOg3zvg5WfTMDn5RYyft43vNhwxw4QHkDA0mibBtt/RV0SqHhUjIlXVvhUw+yHISQcPX7jxTYi+y+Zfy2w6dJqRM5NJOZmLiws80qMJI2Ob4+nuatM4IlJ1qRgRqWosRbDyFVjzDmBASGsY9BkEt7BpmCKLlQ9W7ePd5XuwWA3qBXgzeUg0VzSuY9M4IlL1qRgRqUrOHDZvUj28zhx3GgF9JoGHj03DHD6Vy8iZySSlnAZgQFQ9Jg5sS4CPbVvHi0j1UKbrqFOnTiUyMhJvb29iYmJYv379ZV9/5swZHnvsMcLCwvDy8qJ58+YsXLiwTAmLyCXsmA/TrjILES9/82pI/wSbFiKGYTB74xH6vruapJTT+Hm5kzAkmvfu6KBCRETKrNRXRmbOnElcXBzTpk0jJiaGhIQEevfuza5duwgJCfnL6wsKCrjhhhsICQnh+++/Jzw8nJSUFAIDA22Rv4gU5cPPY2D9R+Y4vBPc/gnUirRpmMzcQl6cs4X5m48D0CWyFpMHRxNR29emcUSk+nExDMMozRtiYmLo0qULU6ZMAcBqtRIREcETTzzBqFGj/vL6adOm8eabb7Jz5048PMr2L6esrCwCAgLIzMzE318tpEWKZew1d9pN3WyOuz8B140Fd0+bhkncd5JnvkvmWGYebq4uPB3bjEd6NsVN+8qIyGWU9PO7VF/TFBQUkJSURGxs7P8fwNWV2NhYEhMTL/qeuXPn0q1bNx577DFCQ0Np27Ytr776KhaL5ZJx8vPzycrKuuBHRP6LYZg77X7cwyxEfOvAXd9Dr5dtWojkFVqYtGgHd/7rd45l5hFZx5cfHunO49c1UyEiIjZTqq9pMjIysFgshIZeuKNnaGgoO3fuvOh79u/fz4oVK7jrrrtYuHAhe/fu5dFHH6WwsJBx48Zd9D2TJk1i/PjxpUlNpPo4fRDmP20u3QWIvNrcadc/zKZhVu5MZ+zcrRw+dQ6AIZ0jGNu/NTW8dN+7iNiW3X+rWK1WQkJC+Pjjj3Fzc6NTp04cPXqUN99885LFSHx8PHFxccXjrKwsIiIi7J2qSMVmKYLfP4CVr0LROXDzgp6j4MqnbLrT7rEz55gwbzuLt6UCEBbgzfgBbejVpq7NYoiI/LdSFSNBQUG4ubmRlpZ2wfNpaWnUrXvxX1RhYWF4eHjg5vb/vyxbtWpFamoqBQUFeHr+9ZKyl5cXXl7ax0Kk2LFkmPckHP/DHEdeDf3fhTpNbBai0GLl098OkLBsD7kFFtxcXbjvqkY8dX0zXQ0REbsq1T0jnp6edOrUieXLlxc/Z7VaWb58Od26XXz78SuvvJK9e/ditVqLn9u9ezdhYWEXLURE5L8U5MLPo2H6dWYh4h0IA6bA8Hk2LUT+c/AU/d5bw6sLd5JbYKFLZC0WPHkVL9zYSoWIiNhdqX/LxMXFMXz4cDp37kzXrl1JSEggJyeHESNGADBs2DDCw8OZNGkSAI888ghTpkzhqaee4oknnmDPnj28+uqrPPnkk7adiUhVs28FzBsJZ1LMcZtboe/rUPOvS+jL6mR2Pq8t2smsJHNPmVq+HsTf2IrbO9bHVTeoioiDlLoYGTJkCCdOnGDs2LGkpqYSHR3N4sWLi29qPXToEK6u/3/BJSIigiVLlvD000/Tvn17wsPDeeqpp3j++edtNwuRqiTnJPz8Ivzxb3PsXx/6TYbmvW0Wwmo1mLnhMK8t2knmuUIA7ugawXO9W1Krhq5YiohjlbrPiDOoz4hUC4YBW2bB4lGQexJwgZiH4LrRNt1ld9uxTEbP2cqmQ2cAaBXmzyu3tKVjg1o2iyEiAiX//NaXwSIVwemU88t1z9+PFdIGBrwH9TvbLMTZvEImL93N52sPYjWghqcbcb1aMLxbQ9zdtMOuiDiPihERZ7IUwbpp5i67hbnmct0ez5nLdd1ss9eLYRgs2HKcifO3k5aVD8BN7cMYc1Nr6gZ42ySGiEh5qBgRcZbjm2HuE3A82RxHXg39EiCoqc1CHMjIYexPW1m9J8MMUceXCTe35ZrmwTaLISJSXipGRBytIBd+eQ3WTgHDAt4BZhv3DveAi21WsOQVWvhw1T4+/GUfBUVWPN1debRnEx7u0QRvD9s1SBMRsQUVIyKOtG8lzB9ptnQHaHML9Hkd/EIv965S+WX3Ccb+tJWUk7kAXNM8mAkD2hAZVMNmMUREbEnFiIgj5J6CJS/CH9+YY/9wuOltaNHXZiFSM/OYMH8bC7eYbdxD/b0Y178NfdvWxcVGV1xEROxBxYiIPRkGbPn+/HLdDMAFuj4I14+x2XLdIouVz9Ye5J2lu8k538Z9RPdIRt7QnJrqnioilYB+U4nYy+kUWBAHe5eZ4+BWMOB9iOhisxBJKad48cet7Ew9C0DHBoG8PLAdreupH4+IVB4qRkRszWoxl+uuePn8cl1PuOb8cl1323Q3PZ1TwGuLdjJzw2EAAn09iO/bkkGdItTGXUQqHRUjIraUusVcrntskzlueKW5u25QM5sc3mo1+D7pCJMW7eB0rtnGfXDn+ozq24raauMuIpWUihERWyg8B6teg7Xvm8t1vQKg1wToMAxcbdPddMfxLEbP2UpSymkAWtb14+WBbekcWdsmxxcRcRYVIyLltX+Vubvu6QPmuPXN0PcN8Ktrk8Nn5xeRsHQ3n649iMVq4OvpRtwNzRnePRIPtXEXkSpAxYhIWeWegp9HQ/LX5tivnrlct+WNNjm8YRgs2prKhHnbSc3KA+DGdnUZ0681YQE+NokhIlIRqBgRKS3DgK0/wKLn/3+5bpf74fqx4G2bVSwpJ3MY+9M2ftl9AoAGtX0Zf3Mbrm0RYpPji4hUJCpGRErjzCGYHwd7l5rj4Fbm7roRXW1y+PwiC9NW7Wfqqr1mG3c3Vx7u2YRHe6qNu4hUXSpGRErCaoF1H51frptzfrnus3DlSJst11295wRjf9rGgYwcAK5qGsSEm9vQOLimTY4vIlJRqRgR+TupW88v191ojht0N5frBje3yeHTsvKYOH878zcfByDEz4sx/VrTr32Y2riLSLWgYkTkUgrPwS9vwNr3wFpkLte9YTx0HG6T5bpFFitfJKYweelusvOLcHWB4d0jibuhOX7eHjaYgIhI5aBiRORiDvwK856CU/vNcasBcOObNluuu/HQaUb/uJXtx7MAiI4I5OWBbWkbHmCT44uIVCYqRkT+W+4pWDoGNn1ljv3C4Ma3oFU/mxz+TG4Bry/exbf/OYRhQICPB8/3acnQLmrjLiLVl4oRETCX626bbS7XzTGX09L5PogdB97lv1phGH+2cd/JqZwCAG7vVJ9RfVsSVNOr3McXEanMVIyInDkMC56BPUvMcVALc7lugytscvhdqWcZM2cr6w+eAqB5aE1eHtiOro3Uxl1EBFSMSHVmtcD66bBiIhRkm8t1r34Grnoa3Mt/tSInv4j3lu9hxpoDFFkNfDzcGBnbjHuvaqQ27iIi/0XFiFRPadtg7pNwdIM5jrjCvBoS3KLchzYMgyXb0pgwbxvHMs027r3bhDK2fxvCA9XGXUTkf6kYkeqlMA9+fQN+e/f8cl1/iH0JOo2wyXLdw6dyGTd3Gyt2pgNQv5YP4we04fpWoeU+tohIVaViRKqPA6vPL9fdZ45b9jOX6/rXK/eh84ssTP91P++v2Et+kRUPNxceuqYJj13bFB9PtXEXEbkcFSNS9Z07DT+PgU1fmuOadeGmt6BVf5scfu3eDEb/tJX9J8w27t2b1GHCzW1pGqI27iIiJaFiRKouw4BtP55frmt+bULne82vZWywXDf9bB6vLNjBT8nHAAiq6cWYfq0YEFVPbdxFREpBxYhUTZlHzOW6uxeb46Dm0P89aNit3Ie2WA2++j2Ft5bs4uz5Nu73XNGQuF4tCPBRG3cRkdJSMSJVi9UC//kXLJ9gLtd19TCX614dZ5Plun8cPsOLc7aw9ajZxr19/QBeGdiOdvXVxl1EpKxUjEjVkbYd5j0JR/5jjiNizKshIS3LfejM3ELe/HknX68z27j7ebvzXJ+W3Nm1AW5q4y4iUi4qRqTyK8yDX9+E3xLM5bqefnDDS9Dp3nIv1zUMgx83HeXVhTvIyDbbuN/aIZz4G1sR7Kc27iIitqBiRCq3g2vM5bon95pjGy7X3ZN2ltFztrLugNnGvWlITSbe3JZuTeqU+9giIvL/VIxI5XTuDCwdCxs/N8c165pFSOsB5T50bkER76/Yy/Rf91NkNfD2cOXJ65tx/1WN8XRXG3cREVtTMSKVi2HA9p9g0XOQnWY+12mEuVzXJ7Dch1+6PY2X5m7j6JlzAMS2CmVc/9ZE1PYt97FFROTiVIxI5ZF5FBb+E3YtNMdBzaH/u9Cwe7kPffhULuPnbWfZDrPACQ/04aUBbbihtdq4i4jYm4oRqfisVtgwA5aNh4Kz55frxplLdsu5XLegyMq/1uznveV7yCu04u7qwgPXNOaJ65ri66n/PUREHEG/baViS99h7q57ZL05rt/V3F03pFW5D5247yRjftrK3vRsAGIa1eblgW1pFupX7mOLiEjJqRiRiqkwD1a/DWveAWuhuVw3dhx0vq/cy3VPnM3n1YU7+HHTUQCCanrywo2tuKVDuNq4i4g4gYoRqXhS1ppXQ07uMcctbjJXygSEl+uwFqvBN+sP8ebinWTlFeHiAnfFNODZXi0J8FUbdxERZ1ExIhXHuTOwbBwkfWaOa4aaRUirAVDOKxZbjmQyes4W/jiSCUC78ABeHtiWqIjAch1XRETKT8WIOJ9hwI65sPA5yE41n+v0D4gdX+7lupnnCpn88y6+/D0FqwF+Xu78s3cL7r6iodq4i4hUECpGxLmyjsGCf8KuBea4TlNzP5nIK8t1WMMwmPvHMSbO30FGdj4AN0fX48WbWhHi513erEVExIZUjIhzWK2Q9Im5XDc/C1zd4aqn4ep/gkf5ioW96dmM/Wkra/edBKBxcA1evrkt3ZsG2SJzERGxMRUj4njpO839ZA7/bo7rdzGvhoS2LtdhzxVYmLJyDx//up9Ci4GXuytPXNeUB65pjJe7mw0SFxEReyjTGsmpU6cSGRmJt7c3MTExrF+//pKv/eyzz3Bxcbngx9tbl8mrpaJ8WDkJpl1lFiKeNaHvm3DvknIXIit2pnHDO78wdeU+Ci0G17UMYVlcDx6/rpkKERGRCq7UV0ZmzpxJXFwc06ZNIyYmhoSEBHr37s2uXbsICQm56Hv8/f3ZtWtX8Vi9HKqhlETzakjG+b8HzfvCTW9BQP1yHfbomXOMn7uNn7ebbdzrBXgzbkAberUO1d8zEZFKotTFyOTJk3nggQcYMWIEANOmTWPBggV88sknjBo16qLvcXFxoW7duuXLVCqnvExY9hJs+MQc1wiBG9+A1gPLtVy30GLlkzUHSFi2h3OFFtxdXbjvqkY8eX0zanjp20cRkcqkVL+1CwoKSEpKIj4+vvg5V1dXYmNjSUxMvOT7srOzadiwIVarlY4dO/Lqq6/Spk2bS74+Pz+f/Pz84nFWVlZp0pSKYsc8WPgsnD1ujjsOgxsmgE+tch12/YFTjJ6zhd1pZhv3rpG1mTiwLS3qqo27iEhlVKpiJCMjA4vFQmjohTuZhoaGsnPnzou+p0WLFnzyySe0b9+ezMxM3nrrLbp37862bduoX//il+gnTZrE+PHjS5OaVCRZx83ddXfON8e1m5i76za6ulyHPZmdz6sLd/LDxiPmYWuYbdxv66g27iIilZndr2d369aNbt26FY+7d+9Oq1at+Oijj5g4ceJF3xMfH09cXFzxOCsri4iICHunKuVltULSp+bXMn8u171yJFzzbLmW61qtBv/+zyHeWLyLzHOFuLjAHV0b8FzvFgT6etosfRERcY5SFSNBQUG4ubmRlpZ2wfNpaWklvifEw8ODDh06sHfv3ku+xsvLCy+v8m0NLw52Ypd5g+qh81/XhXc2d9cNvfTXcSWx9Wgmo+dsJfnwGQBah/nz8i1t6digfF/1iIhIxVGqpb2enp506tSJ5cuXFz9ntVpZvnz5BVc/LsdisbBlyxbCwsJKl6lUTEX5sOo1c7nuoUTwqAF9Xof7fi5XIXI2r5CX5m5jwJQ1JB8+Q00vd8b1b83cx69UISIiUsWU+muauLg4hg8fTufOnenatSsJCQnk5OQUr64ZNmwY4eHhTJo0CYAJEyZwxRVX0LRpU86cOcObb75JSkoK999/v21nIo536Hdzd90/l+s26w03vQ2BZf9KzTAM5m8+zsT520k/a97E3K99GGP6tSbUX/1pRESqolIXI0OGDOHEiROMHTuW1NRUoqOjWbx4cfFNrYcOHcLV9f8vuJw+fZoHHniA1NRUatWqRadOnVi7di2tW5evyZU4UV6m2cZ9wwxzXCMY+r4ObW4t13Ld/SeyGfvTNtbszQCgUVANJtzchqubBdsiaxERqaBcDMMwnJ3E38nKyiIgIIDMzEz8/f2dnU71tmO+uVLmz+W6He6GGyaCb+0yHzKv0MIHK/cy7Zf9FFiseLq78vi1TXnwmsZ4e6h7qohIZVXSz291h5KSyToOi541e4cA1G58frnuNeU67Mpd6Yz7aRuHTuUC0KN5MBNubkPDOjXKm7GIiFQSKkbk8qxW2PgZLH0J8jPN5brdn4Qez4GHT5kPezzzHBPmbWfR1lQA6vp7M65/a/q0raueISIi1YyKEbm0E7vPL9dda47rdYQB70PdtmU+ZKHFyme/HeSdZbvJLbDg5urCvVdG8lRsc2qqjbuISLWk3/7yV0UFsOYdWP0WWArM5brXj4GuD4Jr2e/h2HDwFKPnbGVn6lkAOjWsxcsD29IqTPcBiYhUZypG5EKH1sG8J+HE+fb+zXqdX67boMyHPJVTwGuLdvDdBrONey1fD+L7tuL2TvVxddVXMiIi1Z2KETHlZcHy8fCfGYBhLtft8xq0va3My3WtVoPvNhzmtcU7OZNbCMDQLhE836cltWqojbuIiJhUjAjsXAAL/glnj5nj6LuhV/mW624/lsXoOVvYeOgMAC3r+vHKLW3p1LDsxxQRkapJxUh1djYVFj4LO+aa41qNzOW6jXuU+ZDZ+UW8s3Q3n609iMVqUMPTjadvaM4/ukfi7laq3QdERKSaUDFSHVmtsPFzWDrOXK7r4gZXPgk9ni/zcl3DMFi4JZUJ87eRlmW2cb+pXRij+7UiLKDsS4BFRKTqUzFS3WTsMZfrpvxmjut1NHfXrduuzIc8mJHD2Lnb+HX3CQAa1vFl/IA29GwRYouMRUSkilMxUl0UFcBvCfDrm/+/XPe60RDzUJmX6+YVWpj2yz4+WLWPgiIrnm6uPNKzCY/0bKI27iIiUmIqRqqDw+vN3XVP7DDHTW+AfpPLtVz3190nGPvTVg6eNNu4X90siAk3t6VRkNq4i4hI6agYqcryz8LyCbB+OmCAb5C5u245luumZuYxccF2Fmw2N8oL8fNibP/W3NQuTG3cRUSkTFSMVFW7FsGCZyDrqDmOvgt6vVzm5bpFFiufJ6Yw+edd5BRYcHWBf3RvxNM3NMPP28OGiYuISHWjYqSqOZsGi56D7XPMca1I6JcATa4t8yGTUk4zes5WdhzPAqBDg0BeHtiWNvUCyp2uiIiIipGqwjBg4xewdAzknV+u2/1x6DEKPH3LdMjTOQW8vngn3/7nMAABPh6M6tuSIZ0j1MZdRERsRsVIVZCx9/xy3TXmOCzaXK4bFlWmw1mtBt9vPMJri3ZyKqcAgEGd6jOqb0vq1PSyUdIiIiImFSOVWVEBrH0XfnkTLPng4QvXvggxD4Nb2U7tztQsRv+4lQ0ppwFoEerHy7e0pUuk2riLiIh9qBiprI5sgLlPQPp2c9zkenO5bq3IMh0uJ7+Id5fvYcaaA1isBr6eboyMbcaIKxvhoTbuIiJiRypGKpv8s7B8Iqz/GHO5bh1zd912g8q0XNcwDJZsS2X8vO0cz8wDoE+buozt35p6gWrjLiIi9qdipDLZtfj8ct0j5jjqDuj1CtSoU6bDHTqZy7i5W1m5y2zjHlHbhwkD2nJtS7VxFxERx1ExUhlkp8Oi52HbbHMc2BD6J0CT68p0uPwiCx/9sp+pK/eSX2TFw82Fh3s04bFrm6qNu4iIOJyKkYrMMGDTV/DzaMg7Yy7X7fYY9Iwv83LdNXsyGPvTVvZn5ABwZdM6TLi5LU2Ca9owcRERkZJTMVJRndxnLtc9uNoch0XBgPfLvFw3PSuPiQt2MO+PYwAE+3kx+qZWDIiqpzbuIiLiVCpGKhpLIax9D1a9bi7XdfeB616EmEfKtFy3yGLlq99TePvn3ZzNL8LVBYZ1iySuV3P81cZdREQqABUjFcmRJJj3JKRtNcdNroN+75R5uW7y4TO8+OMWth0z27hHRQTyysC2tA1XG3cREak4VIxUBPnZsOJlWDcNMMCntrlct/3gMi3Xzcwt5PUlO/n3+kMYBvh7u/Ncn5bc0bUBbmrjLiIiFYyKEWfb/TMsiINMc/8X2g+F3q+WabmuYRj8sPEokxbu4OT5Nu63dgznhRtbEaQ27iIiUkGpGHGW7HRYPAq2/mCOAxuaX8k0vb5Mh9uddpbRc7ay/sApAJqF1GTiwLZc0bhsPUhEREQcRcWIoxkGJH8NS148v1zX9b+W69Yo9eFyC863cV99gCKrgY+HG09e34z7rmqEp7vauIuISMWnYsSRTu6D+SPhwK/muG57c3fdeh1KfSjDMPh5exoT5m3n6JlzANzQOpRx/VtTv1bZepCIiIg4g4oRR7AUwtr34ZfXoSjPXK57bTxc8ViZlusePpXLS3O3sXxnOgDhgT6MH9CG2Nahts5cRETE7lSM2NvRJJj7FKRtMceNe0K/BKjdqNSHKiiyMn31ft5fsYe8QrON+wNXN+aJ65rh46k27iIiUjmpGLGX/GxY+Yq5XNewmst1e78KUUPLtFx37b4MxszZyr4TZhv3bo3rMHFgG5qG+Nk6cxEREYdSMWIPe5bC/DjIPGSO2w2GPpOgRlCpD5V+No9XF+xgTrLZxj2opiejb2rNzdFq4y4iIlWDihFbyj5xfrnu9+Y4sMH55bqxpT6UxWrw9boU3lyyi7N5Rbi4wD1XNOSZXi0I8FEbdxERqTpUjNiCYUDyN/Dzi3DutLlc94pH4doXyrRcd/ORM7z441a2HM0EoF14AK/c0pb29QNtnLiIiIjzqRgpr1P7Yd5IOPCLOa7bDvq/B+EdS32ozHOFvLVkF1+tS8EwwM/bned6t+DOmIZq4y4iIlWWipGyshRC4hRY9dr55breZuOybo+BW+m+RjEMgznJR3llwQ4yss027rd0CCf+xpaE+HnbI3sREZEKQ8VIWRzbBHOfgNTzy3Ub9YD+CVC7cakPtTfdbOP++36zjXuT4BpMHNiW7k1Kf7OriIhIZaRipDQKcmDlq/D7B+eX69aCXq9A9J2lXq57rsDC+yv2MH31fgotBt4erjxxXTMeuLqx2riLiEi1omKkpPYug/lPw5nzy3Xb3g59XoOawaU+1LLtaYybu624jfv1LUN4aUAbImqrjbuIiFQ/Kkb+Tk4GLI6HLd+Z44AIc7lusxtKfagjp3MZP287S7enAWYb93H9W9OrTV1bZiwiIlKpqBi5FMOAP76FJS/AuVPmct2Yh+HaF8GrZqkOVVBkZcaaA7y3fA/nCi24u7pw/9WNefL6pvh66hSIiEj1pk/Cizl1wPxKZv9Kcxza1txdN7xTqQ/1+/6TjJmzlT3p2QB0bVSblwe2pXmo2riLiIgAlOlOyalTpxIZGYm3tzcxMTGsX7++RO/79ttvcXFxYeDAgWUJa3+WIvjtXfigm1mIuHvD9ePgwVWlLkQysvOJ+y6ZoR//zp70bOrU8OTtQVHMfPAKFSIiIiL/pdRXRmbOnElcXBzTpk0jJiaGhIQEevfuza5duwgJCbnk+w4ePMg///lPrr766nIlbDfHks8v191sjhtdY+6uW6dJqQ5jsRr8e/0h3li8k6zzbdzv7NqA53q3JMBXbdxFRET+l4thGEZp3hATE0OXLl2YMmUKAFarlYiICJ544glGjRp10fdYLBauueYa7r33XlavXs2ZM2eYM2dOiWNmZWUREBBAZmYm/v7+pUn37/3vcl3vQOj9CkTfVerluluPZvLinK38cfgMAG3q+fPywLZ0aFDLtjmLiIhUAiX9/C7VlZGCggKSkpKIj48vfs7V1ZXY2FgSExMv+b4JEyYQEhLCfffdx+rVq/82Tn5+Pvn5+cXjrKys0qRZcnuXn1+um2KO2952frnupa/wXExWXiGTf97NF4kHsRrg5+XOM72ac0+3SLVxFxER+RulKkYyMjKwWCyEhoZe8HxoaCg7d+686HvWrFnDjBkzSE5OLnGcSZMmMX78+NKkVnqF52DOo5Cdai7XvWkyNO9VpkPN++MYn609CMCAqHqMvqkVIf5q4y4iIlISdl1Nc/bsWe655x6mT59OUFDJ25vHx8cTFxdXPM7KyiIiIsK2yXn4wE1vwcHf4LrRpV6u+9+GdmlA4r6TDO3SgKuaqY27iIhIaZSqGAkKCsLNzY20tLQLnk9LS6Nu3b827tq3bx8HDx6kf//+xc9ZrVYzsLs7u3btokmTv94g6uXlhZeXV2lSK5tW/c2fcnJzdWHKnaXfpVdERERKubTX09OTTp06sXz58uLnrFYry5cvp1u3bn95fcuWLdmyZQvJycnFPwMGDODaa68lOTnZ9lc7REREpNIp9dc0cXFxDB8+nM6dO9O1a1cSEhLIyclhxIgRAAwbNozw8HAmTZqEt7c3bdu2veD9gYGBAH95XkRERKqnUhcjQ4YM4cSJE4wdO5bU1FSio6NZvHhx8U2thw4dwtVVu86KiIhIyZS6z4gz2LXPiIiIiNhFST+/dQlDREREnErFiIiIiDiVihERERFxKhUjIiIi4lQqRkRERMSpVIyIiIiIU6kYEREREadSMSIiIiJOpWJEREREnKrU7eCd4c8msVlZWU7ORERERErqz8/tv2v2XimKkbNnzwJol18REZFK6OzZswQEBFzyzyvF3jRWq5Vjx47h5+eHi4uLQ2NnZWURERHB4cOHq92+ONV17tV13lB9515d5w3Vd+7Vdd7g2LkbhsHZs2epV6/eZTfRrRRXRlxdXalfv75Tc/D39692f2H/VF3nXl3nDdV37tV13lB9515d5w2Om/vlroj8STewioiIiFOpGBERERGnUjHyN7y8vBg3bhxeXl7OTsXhquvcq+u8ofrOvbrOG6rv3KvrvKFizr1S3MAqIiIiVZeujIiIiIhTqRgRERERp1IxIiIiIk6lYkREREScqsoXI1OnTiUyMhJvb29iYmJYv379JV87e/ZsOnfuTGBgIDVq1CA6Opovv/yy+M8LCwt5/vnnadeuHTVq1KBevXoMGzaMY8eOXXCcU6dOcdddd+Hv709gYCD33Xcf2dnZdpvjpThj7pGRkbi4uFzw89prr9ltjhdjy3kDvPTSS7Rs2ZIaNWpQq1YtYmNjWbdu3QWvqYrnHEo296p4zv/bww8/jIuLCwkJCRc8X1XP+X+71Nyr4jn/xz/+8Zc59enT54LXVNVzXpK52/2cG1XYt99+a3h6ehqffPKJsW3bNuOBBx4wAgMDjbS0tIu+fuXKlcbs2bON7du3G3v37jUSEhIMNzc3Y/HixYZhGMaZM2eM2NhYY+bMmcbOnTuNxMREo2vXrkanTp0uOE6fPn2MqKgo4/fffzdWr15tNG3a1LjjjjvsPt//5qy5N2zY0JgwYYJx/Pjx4p/s7Gy7z/dPtp63YRjG119/bSxdutTYt2+fsXXrVuO+++4z/P39jfT09OLXVMVzbhglm3tVPOd/mj17thEVFWXUq1fPeOeddy74s6p6zv90ublXxXM+fPhwo0+fPhfM6dSpUxccp6qe85LM3d7nvEoXI127djUee+yx4rHFYjHq1atnTJo0qcTH6NChgzF69OhL/vn69esNwEhJSTEMwzC2b99uAMZ//vOf4tcsWrTIcHFxMY4ePVqGWZSNM+ZuGOZf2P/9xeVIjph3ZmamARjLli0zDKN6nfP/nbthVN1zfuTIESM8PNzYunXrX+ZY1c/55eZuGFXznA8fPty4+eabL/n6qnzO/27uhmH/c15lv6YpKCggKSmJ2NjY4udcXV2JjY0lMTHxb99vGAbLly9n165dXHPNNZd8XWZmJi4uLgQGBgKQmJhIYGAgnTt3Ln5NbGwsrq6uf7m8bS/OmvufXnvtNerUqUOHDh148803KSoqKvNcSsMR8y4oKODjjz8mICCAqKgooPqc84vN/U9V7ZxbrVbuuecenn32Wdq0afOX91Xlc/53c/9TVTvnAKtWrSIkJIQWLVrwyCOPcPLkyeI/q8rnHC4/9z/Z85xXio3yyiIjIwOLxUJoaOgFz4eGhrJz585Lvi8zM5Pw8HDy8/Nxc3Pjgw8+4IYbbrjoa/Py8nj++ee54447ijcbSk1NJSQk5ILXubu7U7t2bVJTU8s5q5Jx1twBnnzySTp27Ejt2rVZu3Yt8fHxHD9+nMmTJ9tmcpdhz3nPnz+foUOHkpubS1hYGEuXLiUoKAio+uf8cnOHqnnOX3/9ddzd3XnyyScv+v6qfM7/bu5QNc95nz59uPXWW2nUqBH79u3jhRdeoG/fviQmJuLm5lalz/nfzR3sf86rbDFSVn5+fiQnJ5Odnc3y5cuJi4ujcePG9OzZ84LXFRYWMnjwYAzD4MMPP3ROsjZmi7nHxcUVP27fvj2enp489NBDTJo0qUK1Hv5vJZn3tddeS3JyMhkZGUyfPp3Bgwezbt26v/xyqmxsMfeqds6TkpJ499132bhxIy4uLs5O1eZsMfeqds4Bhg4dWvzadu3a0b59e5o0acKqVau4/vrrnZS1bdhi7vY+51W2GAkKCsLNzY20tLQLnk9LS6Nu3bqXfJ+rqytNmzYFIDo6mh07djBp0qQLfjn/+WGckpLCihUrLrgyULduXdLT0y84ZlFREadOnbpsXFty1twvJiYmhqKiIg4ePEiLFi3KPqkSsOe8a9SoQdOmTWnatClXXHEFzZo1Y8aMGcTHx1f5c365uV9MZT/nq1evJj09nQYNGhS/3mKx8Mwzz5CQkMDBgwer7DkvydwvprKf84tp3LgxQUFB7N27l+uvv77KnvOL+d+5X4ytz3mVvWfE09OTTp06sXz58uLnrFYry5cvp1u3biU+jtVqJT8/v3j854fxnj17WLZsGXXq1Lng9d26dePMmTMkJSUVP7dixQqsVisxMTHlmFHJOWvuF5OcnIyrq6tDriDYa95/95qqfM7L8prKfs7vueceNm/eTHJycvFPvXr1ePbZZ1myZAlQdc95SeZ+MZX9nF/MkSNHOHnyJGFhYUDVPecX879zvxibn3O73RpbAXz77beGl5eX8dlnnxnbt283HnzwQSMwMNBITU01DMMw7rnnHmPUqFHFr3/11VeNn3/+2di3b5+xfft246233jLc3d2N6dOnG4ZhGAUFBcaAAQOM+vXrG8nJyRcsccrPzy8+Tp8+fYwOHToY69atM9asWWM0a9bMKcu/HD33tWvXGu+8846RnJxs7Nu3z/jqq6+M4OBgY9iwYZV23tnZ2UZ8fLyRmJhoHDx40NiwYYMxYsQIw8vLy9i6dWvxcariOS/J3KviOb+Yi60kqIrn/GL+d+5V8ZyfPXvW+Oc//2kkJiYaBw4cMJYtW2Z07NjRaNasmZGXl1d8nKp4zksyd0ec8ypdjBiGYbz//vtGgwYNDE9PT6Nr167G77//XvxnPXr0MIYPH148fvHFF42mTZsa3t7eRq1atYxu3boZ3377bfGfHzhwwAAu+rNy5cri1508edK44447jJo1axr+/v7GiBEjjLNnzzpiuhdw9NyTkpKMmJgYIyAgwPD29jZatWplvPrqqxf8z+wItpz3uXPnjFtuucWoV6+e4enpaYSFhRkDBgww1q9ff0HMqnjOSzL3qnjOL+ZixUhVPOcX879zr4rnPDc31+jVq5cRHBxseHh4GA0bNjQeeOCB4g/4P1XFc16SuTvinLsYhmHY5hqLiIiISOlV2XtGREREpHJQMSIiIiJOpWJEREREnErFiIiIiDiVihERERFxKhUjIiIi4lQqRkRERMSpVIyIiIiIU6kYEREREadSMSIiIiJOpWJEREREnErFiIiIiDjV/wFVAZhvZ5XgCAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from cfdmod.analysis.inflow.functions import calculate_spectral_density\n", "\n", "spectral_components = [\"ux\", \"uy\"]\n", "spectral_data = calculate_spectral_density(\n", " inflow_data=profile,\n", " target_index=1,\n", " for_components=spectral_components,\n", " normalization_params=params,\n", ")\n", "\n", "for component in spectral_components:\n", " x = spectral_data[f\"f ({component})\"]\n", " y = spectral_data[f\"S ({component})\"]\n", " plt.plot(x, y, label=f\"{component}\")\n", "\n", "plt.legend()\n", "plt.show()\n", "plt.close()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Autocorrelation" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAGdCAYAAADAAnMpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA+3ElEQVR4nO3deXhU5d3/8fdkmyQkmRDMBgmLisiOUhdwQ0URQcWNzQWLtNUHWtG2PtKn2lprsbVW/T31wa1CXZBNgRZXKoIb2oJAwqpIBAJJWJPJQraZ8/vjhIRAEjIhyT3L53VdczHnzJnkG4ZhPjnnvr+3w7IsCxERERFDwkwXICIiIqFNYURERESMUhgRERERoxRGRERExCiFERERETFKYURERESMUhgRERERoxRGRERExKgI0wU0h9frZe/evcTHx+NwOEyXIyIiIs1gWRbFxcV07tyZsLDGz38ERBjZu3cvmZmZpssQERGRFti9ezcZGRmNPh4QYSQ+Ph6wf5iEhATD1YiIiEhzuN1uMjMzaz/HGxMQYeTopZmEhASFERERkQBzsiEWGsAqIiIiRimMiIiIiFEKIyIiImJUQIwZaQ6Px0NVVZXpMowJDw8nIiJCU59FRCTgBEUYKSkpITc3F8uyTJdiVGxsLOnp6URFRZkuRUREpNkCPox4PB5yc3OJjY0lOTk5JM8MWJZFZWUl+/fvJycnh549ezbZXEZERMSfBHwYqaqqwrIskpOTiYmJMV2OMTExMURGRrJz504qKyuJjo42XZKIiEizBM2vz6F4RuR4OhsiIiKBSJ9eIiIiYpTCiIiIiBilMCIiIiJGKYyIiIiIUQojIiIioWzDPFg6FaorjZUQdGHEsizKKquN3Hxputa9e3eeeeaZevsGDRrEb3/7W1auXElUVBSffvpp7WN/+tOfSElJoaCgoLX+qkREJNTtXA3/+Cmsex02zDVWhk99RmbNmsWsWbP4/vvvAejbty+PPPIII0eObPQ5Cxcu5OGHH+b777+nZ8+e/PGPf+Taa689paKbcqTKQ59HPmizr9+Uzb8bQWzUqbduGTZsGNOnT+eOO+5gw4YN7Nixg4cffpiFCxeSmpraCpWKiEjIO5QD828DTyX0vh7OudNYKT6dGcnIyOCJJ55g7dq1rFmzhiuuuIIbbriBTZs2NXj8F198wYQJE7j77rtZt24dY8aMYcyYMWzcuLFVig9mv//97+nYsSM//vGPuf3225k0aRLXX3+96bJERCQYHCmEueOg7CCkD4IbXwCDvap8+jX+uuuuq7f9+OOPM2vWLL788kv69u17wvHPPvss11xzDb/85S8BeOyxx1i+fDl//etfef7550+h7MbFRIaz+Xcj2uRrN+d7t5aoqCjeeOMNBgwYQLdu3Xj66adb7WuLiEgI81TDwrvgwDZI6AIT5kFUrNGSWnxNwePxsHDhQkpLSxkyZEiDx6xevZoHHnig3r4RI0awZMmSJr92RUUFFRUVtdtut7vZdTkcjla5VNLWwsLCThhjcvyqw1988QUAhw4d4tChQ3To0KHd6hMRkSBkWfDeg7DjY4jsYAeRhHTTVfk+gDU7O5u4uDicTif33HMPixcvpk+fPg0em5+ff8IYh9TUVPLz85v8HjNnzsTlctXeMjMzfS3T7yUnJ5OXl1e77Xa7ycnJqd3+7rvvuP/++3nppZe44IILmDRpEl6v10SpIiISLL56Adb8DXDAzS9B+gDTFQEtCCO9evVi/fr1fPXVV9x7771MmjSJzZs3t2pRM2bMoKioqPa2e/fuVv36/uCKK67gtdde49NPPyU7O5tJkyYRHm5f5vF4PNx+++2MGDGCH/7wh8yePZusrCyeeuopw1WLiEjA+uZD+GCGff+q38HZo8zWcwyfr2dERUVx5plnAjB48GD+85//8Oyzz/LCCy+ccGxaWtoJU1ELCgpIS0tr8ns4nU6cTqevpQWUGTNmkJOTw+jRo3G5XDz22GO1Z0Yef/xxdu7cybJlywBIT0/nxRdfZMKECVx99dUMHDjQZOkiIhJoCjbBoslgeeGcO2DoT01XVM8pD67wer31xncca8iQIXz00UdMnz69dt/y5csbHWMSShISEpg3b169fZMmTaq9/8gjj9R77Kabbmr071lERKRRJftg7nioLIbul8Cov4CfrXTvUxiZMWMGI0eOpGvXrhQXFzN37lxWrlzJBx/YfT3uvPNOunTpwsyZMwG47777uOyyy3jqqacYNWoU8+bNY82aNbz44out/5OIiIhIfVVHYN5EKNoFSWfA2FchIsp0VSfwKYzs27ePO++8k7y8PFwuFwMGDOCDDz7gqquuAmDXrl2EHTNPeejQocydO5df//rX/OpXv6Jnz54sWbKEfv36te5PISIiIvVZlt3mPfc/EJ0IExdAbJLpqhrksHzpYW6I2+3G5XJRVFREQkJCvcfKy8vJycmhR48eREdHG6rQP+jvQkREaq18AlbOhLAIuGMx9Li03Uto6vP7WEG3No2IiEjIy15kBxGA0U8bCSK+UBgREREJJrv/DUv+y74/9Gdwrrk1Z5pLYURERCRYHN5pD1j1VECvUTD8t6YrahaFERERkWBQ7oY3x0PpfkjrDze9CGGtt2ZaW1IYERERCXSearup2b7NEJcGE+aDM850Vc2mMCIiIhLoPvwf2L4cImJgwpvg6mK6Ip8ojIiIiASyf78EXz1v37/pBehyrtl6WkBhREREJFBt/wje+2/7/pWPQJ8bzNbTQgojhnTv3p1nnnmm3r5Bgwbx29/+lsmTJzN69Oh6j1VVVZGSksLf/va3dqxSRET81r6tsPAusDwwcCJc/IDpilrslBfK8zuWBVVlZr53ZGyrLD40ZcoULr30UvLy8khPTwdg2bJllJWVMW7cuFP++iIiEuBKD8DcsVDhhq5D4bpn/G7xO18EXxipKoM/dDbzvX+1F6I6nPKXGTp0KL169eK1117jwQcfBGD27NnceuutxMUFzuhoERFpA9UVMO82KNwJHbvDuNchwmm6qlOiyzR+asqUKcyePRuAgoIC3nvvPSZPnmy4KhERMcqy4B8/hd1fgtMFExdCh06mqzplwXdmJDLWPkNh6ns3U1hYGMevUVhVVVV7/8477+Shhx5i9erVfPHFF/To0YNLLrmk1UoVEZEA9OmfIWs+OMJh7N8h+SzTFbWK4AsjDkerXCppa8nJyeTl5dVuu91ucnJyarc7derEmDFjmD17NqtXr+aHP/yhiTJFRMRfbFoMK35v37/2STjjcrP1tKLgCyMB4oorrmDOnDlcd911JCYm8sgjjxAeXr9t75QpUxg9ejQej4dJkyYZqlRERIzbsxYW32Pfv/C/4Ly7zdbTyhRGDJkxYwY5OTmMHj0al8vFY489Vu/MCMDw4cNJT0+nb9++dO5saFCuiIiYVbgb3pwA1eXQcwRc/XvTFbU6hRFDEhISmDdvXr19x5/9KC0t5fDhw9x9d3AlYBERaaaKYnvxu5ICSOkLt/wtYBa/84XCiB/yer0cOHCAp556isTERK6//nrTJYmISHvzeuCtH0HBRuiQAhPngTPedFVtQmHED+3atYsePXqQkZHBnDlziIjQyyQiEnKWPwLfvAfhTnvxu8SupitqM/qU80Pdu3c/YdqviIiEkLVzYPVf7fs3zoKMHxgtp62p6ZmIiIg/2bES3vm5ff/y/4F+Nxstpz0ojIiIiPiLA9/CgjvBWw39b4VLf2m6onYRNGFElzX0dyAiEtDKDtmL35UXQcb5cP1fA3rxO18EfBg52iissrLScCXmlZXZqxVHRkYarkRERHxSXQnz74BDO+yBquPnQmS06araTcAPYI2IiCA2Npb9+/cTGRlJWFjA5yufWZZFWVkZ+/btIzEx8YROriIi4scsC5bdDzs/g6h4mLgA4pJNV9WuAj6MOBwO0tPTycnJYefOnabLMSoxMZG0tDTTZYiIiC8+fxbWvw6OMLh1DqT0Nl1Ruwv4MAIQFRVFz549Q/pSTWRkpM6IiIgEmi3/hH/91r5/zR+h53Cj5ZgSFGEEICwsjOjo0Lm+JiIiAW7venj7x4AF5/8YLvix6YqMCb0BFiIiIqa599przlSVwRlXwoiZpisySmFERESkPVWW2kGkOA+Sz4ZbZ0N40FyoaBGFERERkfbi9dqXZvI2QOxpMHE+RLtMV2WcwoiIiEh7+ehR2LoMwqPsXiIdu5uuyC8ojIiIiLSHda/D58/Y9294DrpeYLQcf6IwIiIi0ta+/wz+Od2+f+mDMGCs0XL8jcKIiIhIWzr4Hcy/HbxV0PdGGDbDdEV+R2FERESkrRw5bC9+d+QwdBkMY2ZBCC5bcjL6GxEREWkLnipYcCcc3A4JGTD+TYiMMV2VX1IYERERaW2WBe/+AnI+gag4ewpvfKrpqvyWwoiIiEhrW/0crJ1jL353898grZ/pivyawoiIiEhr2vYefPhr+/7Vj0Ova8zWEwAURkRERFpLfjYsuhuwYPAP4cJ7TVcUEBRGREREWkNxPswdD1Wl0OMyuPZJcDhMVxUQFEZEREROVWUZvDkB3LnQqSeM/TuER5quKmAojIiIiJwKrxeW3AN7v4aYJLhtAcR0NF1VQFEYERERORUr/wCbl0JYJIx7HZJON11RwFEYERERaakN8+GTJ+371z0L3S8yW0+AUhgRERFpiZ2r4R/T7PsX3w/n3Ga2ngCmMCIiIuKrnE/gzXHgqYTe18EVj5iuKKBFmC5AREQkoKyfC//4KXirIfNCuPEFLX53ihRGREREmsOy4OM/wCd/srf73mSvwhsZbbauIKAwIiIicjLVFbB0GmQvsLcvfgCueFhnRFqJwoiIiEhTyg7B/Nth5+fgCIfRT8PgSaarCioKIyIiIo05+B3MHQsHt4Mzwe6sesYVpqsKOgojIiIiDdn1FcybAGUHwZUJExdAah/TVQUlhREREZHjbXwbFt8DngpIHwQT50N8mumqgpbCiIiIyFGWBZ89DR89am/3uhZufhmiOpitK8gpjIiIiAB4quCdB+DrV+3tC+6FEY9DWLjZukKAwoiIiEh5ESyYBDs+BkcYXPMEXPAT01WFDIUREREJbYW77Rkz+zZDZCzc8gr0Gmm6qpCiMCIiIqFr7zqYOw5KCiAu1R6o2vkc01WFHIUREREJTVvfhbfuhqoySOlrB5HETNNVhSSFERERCT1fPg/vPwRYdhOzW/8O0QmmqwpZCiMiIhI6vB54fwb8+wV7+9xJMOopCI80W1eIUxgREZHQUFkKi+6Gb96zt4f/Fi6aDg6HyaoE8Gm5wZkzZ3LeeecRHx9PSkoKY8aMYdu2bU0+Z86cOTgcjnq36GgttywiIu2oOB9mj7SDSLgTbp0DF9+vIOInfAojq1atYurUqXz55ZcsX76cqqoqrr76akpLS5t8XkJCAnl5ebW3nTt3nlLRIiIizVawCV66EvI2QGwnuGsZ9L3RdFVyDJ8u07z//vv1tufMmUNKSgpr167l0ksvbfR5DoeDtDT19BcRkXa2/SO7mVllMXTqCbctgKTTTVclx/HpzMjxioqKAEhKSmryuJKSErp160ZmZiY33HADmzZtavL4iooK3G53vZuIiIhP1s6BN261g0i3i+DuDxVE/FSLw4jX62X69OlcdNFF9OvXr9HjevXqxSuvvMLSpUt5/fXX8Xq9DB06lNzc3EafM3PmTFwuV+0tM1PzvkVEpJm8Xlj+G/jnfWB5YMA4uGMxxDb9i7OY47Asy2rJE++9917ee+89PvvsMzIyMpr9vKqqKnr37s2ECRN47LHHGjymoqKCioqK2m23201mZiZFRUUkJGgeuIiINKLqCCy5FzYttrcvewiGPaSBqoa43W5cLtdJP79bNLV32rRpLFu2jE8++cSnIAIQGRnJOeecw/bt2xs9xul04nQ6W1KaiIiEqtID8OYEyP03hEXC9f8LgyaYrkqawafLNJZlMW3aNBYvXsyKFSvo0aOHz9/Q4/GQnZ1Nenq6z88VERFp0IFv4eUr7SAS7bIvyyiIBAyfzoxMnTqVuXPnsnTpUuLj48nPzwfA5XIRExMDwJ133kmXLl2YOXMmAL/73e+48MILOfPMMyksLOTJJ59k586dTJkypZV/FBERCUnffw7zJkJ5ISR2g9sWQnIv01WJD3wKI7NmzQJg2LBh9fbPnj2bu+66C4Bdu3YRFlZ3wuXw4cP86Ec/Ij8/n44dOzJ48GC++OIL+vTpc2qVi4iIZC2ApVPBUwldfgAT5kFcsumqxEctHsDanpo7AEZEREKEZcGqP8HKP9jbva+Hm16EyBizdUk9bTqAVURExJjqSnva7oa59vbQn8HwRyHslFpniUEKIyIiEjiOHIb5d8D3n4IjHK59Es6723RVcooURkREJDAc/t7uqHrgG4iKg1v/Dj2Hm65KWoHCiIiI+L/cNTB3HJQdgPjO9hozaf1NVyWtRGFERET82+al8PaPoboc0gbAxPmQ0Nl0VdKKFEZERMQ/WRZ88b+w/BHAgp4j4JZXwBlnujJpZQojIiLifzzV8N4vYc0r9vZ5P4JrnoBwfWwFI72qIiLiXyqKYeEPYftywAEjHocL/0uL3QUxhREREfEfRXtg7lgo2AgRMXDzy9B7tOmqpI0pjIiIiH/I22DPmCnOgw7J9kDVLoNNVyXtQGFERETM++ZDWHgXVJVC8tkwcQF07Ga6KmknCiMiImLWv1+C9x4Eyws9LoOxr0JMoumqpB0pjIiIiBlejz1td/Vf7e1Bt8PopyEiymxd0u4URkREpP1VlsHbP4Kty+ztK34Nl/xCM2ZClMKIiIi0r+ICeHM87P0awqNgzCzof4vpqsQghREREWk/+7bai90V7YKYjjD+Teg2xHRVYpjCiIiItI8dK2H+nVBRBEmnw22LoNMZpqsSP6AwIiIibW/d6/DP+8BbDZkXwvi50KGT6arETyiMiIhI27EsWPF7+PTP9na/m+GG/4PIaLN1iV9RGBERkbZRVQ5Lp8LGRfb2JT+Hy38NYWFm6xK/ozAiIiKtr+wQzJsIu1ZDWASMfgbOvcN0VeKnFEZERKR1HfzOnjFz6DtwJtgdVc+43HRV4scURkREpPXsXG2fETlyCFxd4bYFkNLbdFXi5xRGRESkdWQvgiX3gqcSOp8DE+ZDfKrpqiQAKIyIiMipsSz49ClY8Zi93WsU3PwSRHUwW5cEDIURERFpOU8VLJtu9xEBuHAqXP0YhIUbLUsCi8KIiIi0zJFCWHAn5KwCRxiM/BOc/yPTVUkAUhgRERHfFe6yZ8zs3wqRHeDW2XDWCNNVSYBSGBEREd/sWQtzx0PpPohLs2fMpA80XZUEMIURERFpvi3L4K0pUH0EUvraQcSVYboqCXAKIyIicnKWBV/Ogg9+BVhw5nC4ZTZEJ5iuTIKAwoiIiDTNUw0fzIB/v2hvD/4hXPtnCNdHiLQO/UsSEZHGVZTAosnw7Qf29lWPwdCfgsNhti4JKgojIiLSMHcezB0L+VkQEQ03vgB9x5iuSoKQwoiIiJwof6MdRNx7IPY0mDAPMs8zXZUEKYURERGpb/u/YMFdUFkMp50FExdAUg/TVUkQUxgREZE6a16Bd34Blge6XwLjXoOYjqarkiCnMCIiIuD1wr9+A1/8P3t74AS47v9BRJTZuiQkKIyIiIS6qiPw9o9hyz/s7WG/gsse1IwZaTcKIyIioaxkP8ybALn/gbBIuOE5GDjOdFUSYhRGRERC1f5v4I1boHAnRCfC+Deg+8Wmq5IQpDAiIhKKcj6F+bdBeRF07A63LYLTepquSkKUwoiISKhZ/yb846fgrYKM8+weIh1OM12VhDCFERGRUGFZsPIJWPWEvd1nDNz4PETGGC1LRGFERCQUVFfAP34GWfPs7Yumw5W/gbAwo2WJgMKIiEjwKzsE8++AnZ+BIxxG/wUG32W6KpFaCiMiIsHs0A54Yywc/Bai4mHs3+HMK01XJVKPwoiISLDa/W94czyUHYSEDLhtAaT2NV2VyAkURkREgtGmxfD2T8BTAWkD7MXuEtJNVyXSIIUREZFgYlnw+bP2OjMAZ42Em18GZ5zZukSaoDAiIhIsPFXw7i9g7Rx7+/yfwDUzISzcaFkiJ6MwIiISDMrdsHASfLcCcMA1T8CF95iuSqRZFEZERAJdUa49Y2bfJoiMhZv/Bmdfa7oqkWZTGBERCWR718Hc8VCSD3Gpdmv3LuearkrEJwojIiKBatt7sGgyVJVBcm976m5iV9NVifhMYUREJBB99QK8/xBYXjj9cruZWbTLdFUiLaIwIiISSLwe+OB/4KtZ9va5d8Kov0B4pNm6RE6BwoiISKCoLIW3psC2d+3tK38DF98PDofZukROkcKIiEggKM63W7vvXQfhTrhxFvS72XRVIq1CYURExN8VbIa5Y6FoN8QkwYQ3oeuFpqsSaTUKIyIi/uy7FbBgElS4IekMuG0hdDrDdFUirUphRETEX339Kiy7H7zV0HUojH8DYpNMVyXS6hRGRET8jdcLKx6Dz/5ib/e/FW54DiKcZusSaSMKIyIi/qSqHJbcC5vetrcvfRAu/5VmzEhQC/Pl4JkzZ3LeeecRHx9PSkoKY8aMYdu2bSd93sKFCzn77LOJjo6mf//+vPvuuy0uWEQkaJUehFevt4NIWATc8H9wxf8oiEjQ8ymMrFq1iqlTp/Lll1+yfPlyqqqquPrqqyktLW30OV988QUTJkzg7rvvZt26dYwZM4YxY8awcePGUy5eRCRoHNgOL18Ju78CpwtufxvOuc10VSLtwmFZltXSJ+/fv5+UlBRWrVrFpZde2uAx48aNo7S0lGXLltXuu/DCCxk0aBDPP/98s76P2+3G5XJRVFREQkJCS8sVEfFPO7+AeRPhyGF7bZnbFkFyL9NViZyy5n5++3Rm5HhFRUUAJCU1Prp79erVDB8+vN6+ESNGsHr16kafU1FRgdvtrncTEQlKWQvh1RvsINJlMEz5SEFEQk6Lw4jX62X69OlcdNFF9OvXr9Hj8vPzSU1NrbcvNTWV/Pz8Rp8zc+ZMXC5X7S0zM7OlZYqI+CfLglVPwttTwFMJva+DScsgLsV0ZSLtrsVhZOrUqWzcuJF58+a1Zj0AzJgxg6Kiotrb7t27W/17iIgYU10JS6fCx7+3t4f+FG59FaJizdYlYkiLpvZOmzaNZcuW8cknn5CRkdHksWlpaRQUFNTbV1BQQFpaWqPPcTqdOJ2aTy8iQehIISy4A3I+AUcYXPsknDfFdFUiRvl0ZsSyLKZNm8bixYtZsWIFPXr0OOlzhgwZwkcffVRv3/LlyxkyZIhvlYqIBLrDO+FvV9tBJCoOJsxXEBHBxzMjU6dOZe7cuSxdupT4+PjacR8ul4uYmBgA7rzzTrp06cLMmTMBuO+++7jssst46qmnGDVqFPPmzWPNmjW8+OKLrfyjiIj4sdy18OY4KN0P8Z1h4nxIH2C6KhG/4NOZkVmzZlFUVMSwYcNIT0+vvc2fP7/2mF27dpGXl1e7PXToUObOncuLL77IwIEDWbRoEUuWLGly0KuISFDZ8k+YM8oOIqn9Ycq/FEREjnFKfUbai/qMiEhAsixY/Vf48GHAgp5Xwy2vgDPedGUi7aK5n99am0ZEpC14quG9B2HN3+ztH9wNI/8E4fpvV+R4eleIiLS2imJYNBm+/RBwwNW/hyFTtcaMSCMURkREWpN7L8wdC/nZEBEDN70Ifa43XZWIX1MYERFpLfnZ8MZYKN4LHZLtqbsZg01XJeL3FEZERFrDt8th4V1QWQKn9YLbFkDH7qarEgkICiMiIqfqP3+Dd38Jlge6XwLjXoOYjqarEgkYCiMiIi3l9cLyh+3puwCDboPRz0BElNGyRAKNwoiISEtUlsHiH9sNzQAu/zVc+gvNmBFpAYURERFfleyDNyfAnjUQHgU3/B8MuNV0VSIBS2FERMQX+7fBG7dA4S57XMi4N6D7RaarEgloCiMiIs2V8wnMvx3Ki6BjD7htEZx2pumqRAKewoiISHOsnwv/+Cl4qyHzAhj/JnToZLoqkaCgMCIi0hTLgpUzYdUf7e2+N8GYWRAZbbYukSCiMCIi0pjqClg6DbIX2NsXPwBXPAxhYWbrEgkyCiMiIg0pO2SPD9n5OTjCYfTTMHiS6apEgpLCiIjI8Q7tgDduhYPbwZkAY/8OZ1xhuiqRoKUwIiJyrF1fwbwJUHYQXJkwcQGk9jFdlUhQUxgRETlq49uw+B7wVED6IJg4H+LTTFclEvQURkRELAs+exo+etTe7nUt3PwyRHUwW5dIiFAYEZHQ5qmCdx6Ar1+1ty+4F0Y8DmHhZusSCSEKIyISusqLYOFd8N0KcITBNU/ABT8xXZVIyFEYEZHQVLgb5o6FfZshMhZueQV6jTRdlUhIUhgRkdCzdx3MHQclBRCXag9U7XyO6apEQpbCiIiElm3vwaLJUFUGKX3tIJKYaboqkZCmMCIioeOrF+D9h8Dy2k3Mbv07RCeYrkok5CmMiEjw83rgg1/BV8/b2+dOglFPQXik2bpEBFAYEZFgV1kKb02Bbe/a28MfhYvuA4fDbF0iUkthRESCV3G+PVA1bz2EO+GmF6DvjaarEpHjKIyISHAq2GxP3S3aDbGdYMI8yDzfdFUi0gCFEREJPts/spuZVbihU0+4bQEknW66KhFphMKIiASXtX+HZfeD5YFuF8G41yE2yXRVItIEhRERCQ5eL6z4nb3gHcCAcXD9/0KE02xdInJSCiMiEviqymHJPbBpsb192UMw7CHNmBEJEAojIhLYSg/AvImw+ysIi7TPhgyaYLoqEfGBwoiIBK4D38Ibt8LhHIh2wbg3oMclpqsSER8pjIhIYPr+c/uMSHkhJHaD2xZCci/TVYlICyiMiEjgyVoAS6eCpxK6/MDuIRKXbLoqEWkhhRERCRyWBZ88CR8/bm/3vh5uehEiY8zWJSKnRGFERAJDdSX88z7YMNfeHvoze52ZsDCzdYnIKVMYERH/d+QwzL8Dvv8UHOEw6s/wg8mmqxKRVqIwIiL+7fD39oyZA99AVBzc+nfoOdx0VSLSihRGRMR/5a6xV90tOwDxne01ZtL6m65KRFqZwoiI+KfNS+HtH0N1OaQNgInzIaGz6apEpA0ojIiIf7EsWP1X+PBhwIKeI+CWV8AZZ7oyEWkjCiMi4j881fDeL2HNK/b2eT+Ca56AcP1XJRLM9A4XEf9QUQwLfwjblwMOGPEHuPBeLXYnEgIURkTEvKI9MHcsFGyEiBi4+WXoPdp0VSLSThRGRMSsvCw7iBTnQYcUmDgPugw2XZWItCOFEREx55sPYeFdUFUKyWfDxAXQsZvpqkSknSmMiIgZ/34J3nsQLC+cPsxuZhaTaLoqETFAYURE2pfXA8sfsafvApxzO4x+BsIjjZYlIuYojIhI+6ksg7d/BFuX2dtXPAyX/FwzZkRCnMKIiLSPkn12a/e9X0N4FIyZBf1vMV2ViPgBhRERaXv7ttqL3RXtgpgkGD8Xug0xXZWI+AmFERFpWztWwvw7oaIIks6A2xZCpzNMVyUifkRhRETazrrX4Z/3gbcaug6xz4jEJpmuSkT8jMKIiLQ+y4IVv4dP/2xv97sFbngOIqPN1iUifklhRERaV3UFLPkv2LjI3r7kF3D5/0BYmNm6RMRvKYyISOspOwTzJsKu1RAWAdc9a/cRERFpgsKIiLSOg9/ZM2YOfQdOF4x71e6sKiJyEgojInLqdq62z4gcOQSurvaMmZSzTVclIgFCYURETk32IlhyL3gqofO5MGEexKearkpEAojCiIi0jGXBp0/Bisfs7bNHw00vQVSs2bpEJOAojIiI7zxVsOx+WPeavT1kGlz1OwgLN1uXiAQkn+faffLJJ1x33XV07twZh8PBkiVLmjx+5cqVOByOE275+fktrVlETDpSCK/fbAcRRxhc+2cY8biCiIi0mM9nRkpLSxk4cCCTJ0/mpptuavbztm3bRkJCQu12SkqKr99aREwr3GXPmNm/FSI7wK2z4awRpqsSkQDncxgZOXIkI0eO9PkbpaSkkJiY6PPzRMRP7FkLc8dD6T6IT4eJ8yF9oOmqRCQItFtLxEGDBpGens5VV13F559/3uSxFRUVuN3uejcRMWjrOzB7lB1EUvvBlI8URESk1bR5GElPT+f555/nrbfe4q233iIzM5Nhw4bx9ddfN/qcmTNn4nK5am+ZmZltXaaINMSyYPX/wbzboPoInDkcJr8Pri6mKxORIOKwLMtq8ZMdDhYvXsyYMWN8et5ll11G165dee211xp8vKKigoqKitptt9tNZmYmRUVF9cadiEgbOlJoL3b3n5fs7R9MhpFPQrgm4YlI87jdblwu10k/v438r3L++efz2WefNfq40+nE6XS2Y0UiAtiL3H27HLLmwzfv243MAK56DIb+FBwOs/WJSFAyEkbWr19Penq6iW8tIsezLNj9lR1ANr4N5YV1jyX3hisfhrNHGStPRIKfz2GkpKSE7du3127n5OSwfv16kpKS6Nq1KzNmzGDPnj28+uqrADzzzDP06NGDvn37Ul5ezssvv8yKFSv48MMPW++nEBHfHfjWDiBZC6BwZ93++HTofwsMGGcPVtXZEBFpYz6HkTVr1nD55ZfXbj/wwAMATJo0iTlz5pCXl8euXbtqH6+srOTnP/85e/bsITY2lgEDBvCvf/2r3tcQkXZSss8++5E1D/auq9sfFQd9boABY6H7JWpgJiLt6pQGsLaX5g6AEZEGVJbC1nftsyDfrQDLY+93hNuzYwaMhV7Xak0ZEWl1fj2AVUTamNcDOatgw3zY8k+oKq17rMsP7EswfW+EuGRzNYqI1FAYEQkWlgX5WfYYkOxFUHLM+k8de9gBZMBY6HSGuRpFRBqgMCIS6Ap3Q/ZC+zLM/q11+2OSoN9NdgjJOE8DUUXEbymMiASiI4Wweal9FmTnMT17wp3QayQMHA9nXAkRUcZKFBFpLoURkUBRXQnbl8OGefDNB+A52qXYAd0vts+A9Lkeol1GyxQR8ZXCiIg/O7Yh2abFcORw3WMpfewA0v8WcGWYq1FE5BQpjIj4IzUkE5EQojAi4i9K9sPGt+wQsveYVa3VkExEgpzCiIhJlWWw9Z3GG5INHAdnjVRDMhEJagojIu3taEOyrAV2Q7LKkrrHjjYk63cTdDjNXI0iIu1IYUSkPVgW5GfbZ0DUkExEpB6FEZG2VNuQbAHs31K3Xw3JRERqKYyItLamGpKdfa0dQNSQTESklsKISGs42pAsaz5se18NyUREfKAwItJSakgmItIqFEZEfHXShmTjIa2fufpERAKMwohIc5Tsg41vqyGZiEgbUBgRaUxlKWx9Vw3JRETamMKIyLGONiTbMN9uSFZVWveYGpKJiLQJhRERy4L8LHsMiBqSiYi0O4URCV2Fu45pSLa1br8akomItCuFEQktakgmIuJ3FEYk+FVXwLc1Dcm+eR88lTUPqCGZiIg/UBiR4GRZsOvLuoZk5YV1j6khmYiIX1EYkeCy/xvIXmCHkMJddfvVkExExG8pjEjgK9kHG9+qaUi2rm6/GpKJiAQEhREJTJWlsPWdmoZkH6shmYhIAFMYkcDhqbYbkmUtUEMyEZEgojAi/s2yIG+DHUA2LoKSgrrH1JBMRCQoKIyIf1JDMhGRkKEwIv7jyOFjGpJ9XrdfDclERIKawoiYVV0B335Y05DsAzUkExEJQQoj0v68Xtj9lRqSiYgIoDAi7Wn/N3YAyV6ghmQiIlJLYUTa1tGGZBvmQd76uv1qSCYiIjUURqT1qSGZiIj4QGFEWkdtQ7L5sGWZGpKJiEizKYxIy6khmYiItAKFEfFd4S47gGQtgAPb6vbXNiQbDxk/UEMyERFpFoURaR41JBMRkTaiMCKNU0MyERFpBwojUl9tQ7J5NQ3JiuoeU0MyERFpAwojYmuyIdmtdghRQzIREWkDCiOhrNGGZPH25Rc1JBMRkXagMBJqGmtIFhZhNyQbMFYNyUREpF0pjIQCTzXkrLRnwqghmYiI+BmFkWClhmQiIhIgFEaCzeGdkL1QDcnklFmWRXFFNfHOCBz69yIibUhhJBgcbUi2YT7s+qJuf0Q09BqphmRSz5FKD/uLK9hfUm7/efRWUsE+t/3n/uIKDpRUUOWxSIiOYEBGIv0zXAzMcDEgI5F0V7QCioi0GoWRQNVUQ7Iel9gBpPd1akgWIqo9Xg6VVrKvJlTsPyZUHBs29hdXUFJR7dPXdpdX89n2A3y2/UDtvtPiohiQkciADFfNLZHT4pyt/WOJSIhQGAkkXi/s/tIOIGpIFvQsy8J9pJr9JeV2yDguVBy9HSip4GBpJZbV/K/tjAgjJcFJcpyT5PiaW1z0CftcMZHkHCglK7eIrNxCsnKL2FZQzIGSSlZs3ceKrftqv2ZnV7QdUDJdDOhin0lxxUS2wd+MiAQbh2X58l+YGW63G5fLRVFREQkJCabLaX9qSBZUyqs8DYaKhrYrq73N/rphDugU5ySlNlwcEzSO2447hXEg5VUeNu11k5VbSHZuERtyC9lxoLTBMNTjtA4MyHDRv4uLgZmJ9O2cQGyUfgcSCRXN/fxWGPFXxQV2Q7Ks+WpIFgA8XotDpZWNhop97vLa/cXlvl0mSYiOOCZURDcaMpI6RBEeZmYcR3F5FRv32AEla499FmX3oSMnHBfmgJ4p8falncxEBnRxcXZ6PM4I/TsWCUYKI4GoXkOyFWDV/FashmRGHJ1NUi9YNBI2DpZU4PXhnRQVEdbkmYujZzdOi3MSHRmYH9SHSivJ3lNE1u5CNuQWkb2nkAJ3xQnHRYY76J2eYJ89qbnMc2ZyHBHhYQaqFpHWpDASKJpqSJZxnn0Jpu+NakjWiiqqPRwoqWwgZJSfEDjKq5p/mcThgE4dojitXqiIbjBsJESH5nTZAnc5G3YXkr2niA0141AKy6pOOC4mMpy+nRPqDZLt3qkDYYbO/IhIyyiM+LPahmTzIXsRlNYNAlRDspbxei0Ol1U2OdDz6HbRkRM//JoS77Qvk5zWxFiMlJrLJPpt3jeWZZF7+AgbagbHZuUWsnGPu8EZP/HREfTvYs/cGZjhon+Giy6JMSEZ6kQChcKIP2qyIdnNdghRQ7JalmVRerQnRu2tvH7IqO2JUYnHh+skkeGOE0PF0bMYcXWXSk6LcxITFZiXSQKV12ux40BJTTixA8qmvW4qGhjM26lDlD1AtiagDMhIJDleU4xF/IXCiL84chg2LbEDiBqSAVBZ7eVgaf2zFo2d0ThS5fHpa3fqENXgZZHjt10xkfqNOoBUebx8U1BcM3vHDijb8oupbiCApruia3ufDMiwpxm7YjXFWMQEhRGTQrAhmWVZFJZV1YWKBrp7Hr1/uIExAk3pEBXexEDPujEZSR2iiNRlkpBRXuVhS5673hmU7ftLGpxi3L1TbL2zJ/26aIqxSHtQGGlvTTYk6wsDx0G/W8DVxVyNLVBWWd1oqKh/mcRuHd5cEWGOBsPFacf2yajZ7uDUh4Y0T0lFNZv22OFkQ649UHbnwbITjgtzwJkpcccMkE2kt6YYi7Q6hZH2sn+bfQkmgBqSVXu8HCytbDJk7Cu2z2yUVvp2maRjbGQTl0nqzmIkxkRqZoS0i8KySrJyi+wZPLvtgbL57vITjosMd9ArLb5ugGyXRM5K1RRjkVOhMNKWmmxIdkNNQ7KL27UhWaOtwxsIG4fKfGsdHh0ZVnc5pKGZJAn2/U4dnERF6D9u8X/73OV1Le5rzqQcKq084bjoyDD6dnbVW4Onh6YYizSbwkhrqyipa0i24+N2a0h2tHX4viYukxyouV/paX5PjPAwR73BnvVbiEfXCxwdosI12FOC2tEpxnb/k0KydhexcU8RxQ1NMXZG0K+LiwGZdpO2/l1cZHTUFGORhiiMtIY2akjm8VonzCZpcCyGu6LB/wybkhAdQUpC4y3Dj946xpprHS4SCLxei5yDpbULBGblFrFpb1GDjfCSaqYYD6jpgzIg00VKfLSBqkX8S5uFkU8++YQnn3yStWvXkpeXx+LFixkzZkyTz1m5ciUPPPAAmzZtIjMzk1//+tfcddddzf6e7RpGLMu+9JK1oOGGZAPH22NBjmtIdrR1+D530wM99xdXcKjU99bhJ1v87GjXz0BtHS4SCKo9Xr7dV0JWbk2L+9witua7Gxy8nZZgTzEemJlY06zNRWJs6EzhF4Hmf377PE2htLSUgQMHMnnyZG666aaTHp+Tk8OoUaO45557eOONN/joo4+YMmUK6enpjBgxwtdv33YO77QHoWYtgAPf1O62YpIo7XkDe7pex87oPuwvrWT/ugr2F2efEDIaasrUGLt1+IlrkZzYiMtJ/CmssCoirSciPIze6Qn0Tk9g3Hn2vvIqD1vzi8muCShZuYVs31dCvruc/M3lfLi5oPb53TrF1q3Bk+GiXxeXZouJcIqXaRwOx0nPjPz3f/8377zzDhs3bqzdN378eAoLC3n//feb9X3a6szIjl27qd64hI7bF5N8aG3t/kqi+DzifN6qvoj3y/tS7UNmiz+6wupJLpMkxap1uEiwKq2oZtNe9zFnUAr5voEpxg4HnJlsTzEemOmifxcXvdMTdIZTgkabnRnx1erVqxk+fHi9fSNGjGD69OmNPqeiooKKirrVPd1ud+sX5vXimn0pnaxD9qblYLW3D0u8F/G+53yKqRuIGhUeVrc2yXGh4vgzGvpPREQ6OCM4v0cS5/dIqt1XVFZF1p66NXiycovIKyrn230lfLuvhLe+zgXsHjxnp8fTv0tdk7aeqXFq6CdBrc3DSH5+PqmpqfX2paam4na7OXLkCDExMSc8Z+bMmTz66KNtW1hYGFtcF9PZncVX8cP5JmUEUUmZ9IpzcvGxQSMumoQYXSYRkVPjio3kkp7JXNIzuXbfvuLy2hb32TUB5WBpJRv3uNm4x82b/7aPc0aEHbeKcSKnn6YpxhI8/PJi5YwZM3jggQdqt91uN5mZma3+fS6e9jJEODm91b+yiMjJpcRHc2XvaK7sbf/CZlkWewqP1FuDJzvXnmL89a5Cvt5VWPvcOGcE/bok1Iw/sUOKphhLoGrzMJKWlkZBQUG9fQUFBSQkJDR4VgTA6XTidLbDypsRWt1TRPyHw+Ego2MsGR1jGdk/HbCnGH9/sLSuxX1uERv3FlFSUc2XOw7x5Y5Dtc9P6hBVO3PnaCfZlARNMRb/1+ZhZMiQIbz77rv19i1fvpwhQ4a09bcWEQl4YWEOTk+O4/TkOMacY69tVe3xsn1/CVm769bg2ZLn5lBpJau+2c+qb/bXPj81wVnX4j4jkQFdXHTsoCnG4l98nk1TUlLC9u3bATjnnHP4y1/+wuWXX05SUhJdu3ZlxowZ7Nmzh1dffRWwp/b269ePqVOnMnnyZFasWMHPfvYz3nnnnWZP7fWLDqwiIn6sotrD1rxiu719zRo83+4rbrCnUdekWPpnuI5ZxdhFnKYYSxtos6ZnK1eu5PLLLz9h/6RJk5gzZw533XUX33//PStXrqz3nPvvv5/NmzeTkZHBww8/7L9Nz0REgkRZpT3FeMNu++xJVm4ROQdKTzjO4YAzkuPqushmJtJHU4ylFagdvIiInKCorIqNe+vW4MneU8SewiMnHBcR5uCs1HgGZtpnT/p3cdErLV5TjMUnCiMiItIs+4sryN5TtwZPVm4hB0pOXMXYGRFGn84JtWvwDMx0cfppcZpiLI1SGBERkRaxLIu8ovLaDrJHm7QVl5+4cGeHqHD6dbHX4LEv8ySSmaQpxmJTGBERkVbj9VrsPFR2zCrGhWzc4+ZIleeEYzvGRtbO3Dm6WGCqphiHJIURERFpU9UeL9/tL7XHn9T0QNmSV0yl58RFQ1Pij51ibC8WqCnGwU9hRERE2l1FtYdv8ktqA0pWbhHfFDQ8xTgzKYYBXepa3PfrkkB8dGT7Fy1tRmFERET8wpFKD5v21l+DZ0cjU4xPP60DAzMS6V8TUPp21hTjQKYwIiIifqvoSBWb9hTVGyDb0BTj8KNTjGvCyYAMTTEOJAojIiISUA6UVJB9zPTiDblFHCipOOG4qIgw+qQn1F7eGZDh4ozkOMI1xdjvKIyIiEhAsyyLfHc5G3bXrGC8p4gNuwtxNzLFuG+Xug6yAzNcdE2K1RRjwxRGREQk6FiWxc6DZbUrGGfVrGJcVnniFGNXTGTN2RMX/bvYTdrSEqIVUNqRwoiIiIQEj9fiu/0lbKhZIDBrTxFb9robnGKcHO+s7SA7INM+k9Ipzmmg6tCgMCIiIiGrstrLtvxisvbYa/BsyC3k230leBqYY9wlMYaBmTVnTzJc9MtwkaApxq1CYUREROQYRyo9bM4rqh2DkrWniB37T5xiDPYU42MHyPbt7CImSlOMfaUwIiIichLu8io21lzaycotZMPuxqcY90yJqw0oAzMS6ZUWT1SEphg3RWFERESkBQ6WVNjhZHcR2XvsKcb7ixuYYhwextnp8fROS6BP5wR6pydwdnq8LvEcQ2FERESkFRw7xTh7z9GFAosoOlLV4PGZSTH0SbfDydE/MzqG5krGCiMiIiJt5OgU4815bjbvdbMlz83mPDd5ReUNHh8fHVEbTo4GlJ6pcUHf6l5hREREpJ0dLq2sDSab89xsyStm+75iqjwnftSGhzk4I7lD3VmUmks9pwXRVGOFERERET9QWe1l+76S2pBy9M/CsoYv86TEO+uFkz7pCfQ4rUNAtrtXGBEREfFTR8ehHHuJZ0teMd8fLKWhT+XoyDB6pcbXCyhnpycQ54xo/+J9oDAiIiISYEorqtmaX1zvLMrWvGKOVJ3Y7h6gW6fYerN5+nROoLPLf1reK4yIiIgEAY/XYufBUrbkFbM5r8j+c6+bfHfDg2UTjg6WPeYsSs/UOJwR7T9YVmFEREQkiB2qGSy7pWZGz+Y8N9v3lVDdQMv7iDAHZ6bE1Ztu3Ds9vs3X5VEYERERCTEV1Z6awbLF9cajNNYTJTXBWRtObv1BJj1O69Cq9TT389u/R76IiIhIszkjwunb2V5Lh8H2PsuyyCuqGyy7Jd8+k/L9wTIK3BUUuPfz8bb9XHF2SquHkeZSGBEREQliDoeDzokxdE6MYXif1Nr9JRXVbMt3s7nmLEqvtHhjNSqMiIiIhKA4ZwSDuyUxuFuS6VLQcoMiIiJilMKIiIiIGKUwIiIiIkYpjIiIiIhRCiMiIiJilMKIiIiIGKUwIiIiIkYpjIiIiIhRCiMiIiJilMKIiIiIGKUwIiIiIkYpjIiIiIhRCiMiIiJiVECs2mtZFgBut9twJSIiItJcRz+3j36ONyYgwkhxcTEAmZmZhisRERERXxUXF+NyuRp93GGdLK74Aa/Xy969e4mPj8fhcLTa13W73WRmZrJ7924SEhJa7etKy+j18D96TfyLXg//otfj5CzLori4mM6dOxMW1vjIkIA4MxIWFkZGRkabff2EhAT9Q/Ijej38j14T/6LXw7/o9WhaU2dEjtIAVhERETFKYURERESMCukw4nQ6+c1vfoPT6TRdiqDXwx/pNfEvej38i16P1hMQA1hFREQkeIX0mRERERExT2FEREREjFIYEREREaMURkRERMSooA4jn3zyCddddx2dO3fG4XCwZMmSkz5n5cqVnHvuuTidTs4880zmzJnT5nWGCl9fj5UrV+JwOE645efnt0/BQW7mzJmcd955xMfHk5KSwpgxY9i2bdtJn7dw4ULOPvtsoqOj6d+/P++++247VBv8WvJ6zJkz54T3R3R0dDtVHNxmzZrFgAEDahuaDRkyhPfee6/J5+i90XJBHUZKS0sZOHAgzz33XLOOz8nJYdSoUVx++eWsX7+e6dOnM2XKFD744IM2rjQ0+Pp6HLVt2zby8vJqbykpKW1UYWhZtWoVU6dO5csvv2T58uVUVVVx9dVXU1pa2uhzvvjiCyZMmMDdd9/NunXrGDNmDGPGjGHjxo3tWHlwasnrAXb3z2PfHzt37mynioNbRkYGTzzxBGvXrmXNmjVcccUV3HDDDWzatKnB4/XeOEVWiACsxYsXN3nMgw8+aPXt27fevnHjxlkjRoxow8pCU3Nej48//tgCrMOHD7dLTaFu3759FmCtWrWq0WPGjh1rjRo1qt6+Cy64wPrJT37S1uWFnOa8HrNnz7ZcLlf7FRXiOnbsaL388ssNPqb3xqkJ6jMjvlq9ejXDhw+vt2/EiBGsXr3aUEUCMGjQINLT07nqqqv4/PPPTZcTtIqKigBISkpq9Bi9R9pPc14PgJKSErp160ZmZmaTv7lLy3k8HubNm0dpaSlDhgxp8Bi9N06Nwsgx8vPzSU1NrbcvNTUVt9vNkSNHDFUVutLT03n++ed56623eOutt8jMzGTYsGF8/fXXpksLOl6vl+nTp3PRRRfRr1+/Ro9r7D2icTytq7mvR69evXjllVdYunQpr7/+Ol6vl6FDh5Kbm9uO1Qav7Oxs4uLicDqd3HPPPSxevJg+ffo0eKzeG6cmIFbtldDUq1cvevXqVbs9dOhQvvvuO55++mlee+01g5UFn6lTp7Jx40Y+++wz06UIzX89hgwZUu839aFDh9K7d29eeOEFHnvssbYuM+j16tWL9evXU1RUxKJFi5g0aRKrVq1qNJBIy+nMyDHS0tIoKCiot6+goICEhARiYmIMVSXHOv/889m+fbvpMoLKtGnTWLZsGR9//DEZGRlNHtvYeyQtLa0tSwwpvrwex4uMjOScc87Re6SVREVFceaZZzJ48GBmzpzJwIEDefbZZxs8Vu+NU6MwcowhQ4bw0Ucf1du3fPnyRq8RSvtbv3496enppssICpZlMW3aNBYvXsyKFSvo0aPHSZ+j90jbacnrcTyPx0N2drbeI23E6/VSUVHR4GN6b5wi0yNo21JxcbG1bt06a926dRZg/eUvf7HWrVtn7dy507Isy3rooYesO+64o/b4HTt2WLGxsdYvf/lLa8uWLdZzzz1nhYeHW++//76pHyGo+Pp6PP3009aSJUusb7/91srOzrbuu+8+KywszPrXv/5l6kcIKvfee6/lcrmslStXWnl5ebW3srKy2mPuuOMO66GHHqrd/vzzz62IiAjrz3/+s7VlyxbrN7/5jRUZGWllZ2eb+BGCSktej0cffdT64IMPrO+++85au3atNX78eCs6OtratGmTiR8hqDz00EPWqlWrrJycHCsrK8t66KGHLIfDYX344YeWZem90dqCOowcnRp6/G3SpEmWZVnWpEmTrMsuu+yE5wwaNMiKioqyTj/9dGv27NntXnew8vX1+OMf/2idccYZVnR0tJWUlGQNGzbMWrFihZnig1BDrwVQ79/8ZZddVvv6HLVgwQLrrLPOsqKioqy+ffta77zzTvsWHqRa8npMnz7d6tq1qxUVFWWlpqZa1157rfX111+3f/FBaPLkyVa3bt2sqKgoKzk52bryyitrg4hl6b3R2hyWZVntdx5GREREpD6NGRERERGjFEZERETEKIURERERMUphRERERIxSGBERERGjFEZERETEKIURERERMUphRERERIxSGBERERGjFEZERETEKIURERERMUphRERERIz6/1+W8KCn7qCmAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from cfdmod.analysis.inflow.functions import calculate_autocorrelation\n", "\n", "autocorrelation_components = [\"ux\", \"uy\"]\n", "autocorrelation_data = calculate_autocorrelation(\n", " inflow_data=profile, anchor_point_idx=0, for_components=autocorrelation_components\n", ")\n", "\n", "for component in autocorrelation_components:\n", " x = profile.points[\"z\"]\n", " y = autocorrelation_data[f\"coef_{component}\"]\n", " plt.plot(x, y, label=f\"{component}\")\n", "\n", "plt.legend()\n", "plt.show()\n", "plt.close()" ] } ], "metadata": { "kernelspec": { "display_name": "cfdmod-XMkUSlb0-py3.10", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.12" } }, "nbformat": 4, "nbformat_minor": 2 }