diff --git a/Druckrohrleitung/Druckrohrleitung_class_file.py b/Druckrohrleitung/Druckrohrleitung_class_file.py index 7cf4c12..44e759d 100644 --- a/Druckrohrleitung/Druckrohrleitung_class_file.py +++ b/Druckrohrleitung/Druckrohrleitung_class_file.py @@ -112,12 +112,16 @@ class Druckrohrleitung_class: self.p[0] = p_boundary_res self.p[-1] = p_boundary_tur - def set_steady_state(self,ss_flux,ss_level_reservoir,pl_vec,h_vec): + def set_steady_state(self,ss_flux,ss_level_reservoir,area_reservoir,pl_vec,h_vec): # set the pressure and velocity distributions, that allow a constant flow of water from the (steady-state) reservoir to the (steady-state) turbine # the flow velocity is given by the constant flow through the pipe ss_v0 = np.full(self.n_seg+1,ss_flux/self.A) - # the static pressure is given by the hydrostatic pressure, corrected for friction losses and dynamic pressure - ss_pressure = self.density*self.g*(ss_level_reservoir+h_vec)-ss_v0**2*self.density/2-(self.f_D*pl_vec/self.dia*self.density/2*ss_v0**2) + + # the static pressure is given by static state pressure of the reservoir, corrected for the hydraulic head of the pipe and friction losses + ss_v_in_res = ss_flux/area_reservoir + ss_v_out_res = ss_flux/self.A + ss_pressure_res = self.density*self.g*(ss_level_reservoir)+self.density*ss_v_out_res*(ss_v_in_res-ss_v_out_res) + ss_pressure = ss_pressure_res+(self.density*self.g*h_vec)-(self.f_D*pl_vec/self.dia*self.density/2*ss_v0**2) self.set_initial_flow_velocity(ss_v0) self.set_initial_pressure(ss_pressure) diff --git a/Druckrohrleitung/Druckrohrleitung_test_steady_state.ipynb b/Druckrohrleitung/Druckrohrleitung_test_steady_state.ipynb index ff4719f..299c49b 100644 --- a/Druckrohrleitung/Druckrohrleitung_test_steady_state.ipynb +++ b/Druckrohrleitung/Druckrohrleitung_test_steady_state.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": null, + "execution_count": 9, "metadata": {}, "outputs": [], "source": [ @@ -22,7 +22,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 10, "metadata": {}, "outputs": [], "source": [ @@ -35,8 +35,8 @@ "L = 1000. # length of pipeline [m]\n", "D = 0.9 # pipe diameter [m]\n", "h_res = 10. # water level in upstream reservoir [m]\n", - "n = 50000 # number of pipe segments in discretization\n", - "nt = 12 # number of time steps after initial conditions\n", + "n = 50 # number of pipe segments in discretization\n", + "nt = 1000 # number of time steps after initial conditions\n", "f_D = 0.01 # Darcy friction factor\n", "c = 400. # propagation velocity of the pressure wave [m/s]\n", "h_pipe = 105. # hydraulic head without reservoir [m] \n", @@ -69,24 +69,39 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 11, "metadata": {}, "outputs": [], "source": [ - "V = Ausgleichsbecken_class(area_base, area_pipe, critical_level_low, critical_level_high,simulation_timestep)\n", + "V = Ausgleichsbecken_class(area_base,area_pipe,critical_level_low,critical_level_high,simulation_timestep)\n", "V.set_steady_state(initial_flux,initial_level,conversion_pressure_unit)\n", "\n", "pipe = Druckrohrleitung_class(L,D,n,alpha,f_D)\n", "pipe.set_pressure_propagation_velocity(c)\n", "pipe.set_number_of_timesteps(nt)\n", - "pipe.set_steady_state(initial_flux,initial_level,pl_vec,h_vec)" + "pipe.set_steady_state(initial_flux,initial_level,area_base,pl_vec,h_vec)" ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "0.8\n", + "0.8\n", + "10.0\n", + "97309.32104824147\n", + "96532.05564888839\n", + "96532.05564888839\n", + "[0.8 0.8 0.8 ... 0.8 0.8 0.8]\n", + "[1.25752054 1.25752054 1.25752054 ... 1.25752054 1.25752054 1.25752054]\n" + ] + } + ], "source": [ "print(V.get_current_influx())\n", "print(V.get_current_outflux())\n", @@ -100,7 +115,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 13, "metadata": {}, "outputs": [], "source": [ @@ -131,7 +146,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 14, "metadata": {}, "outputs": [], "source": [ @@ -154,22 +169,2011 @@ }, { "cell_type": "code", - "execution_count": 22, + "execution_count": 15, "metadata": {}, "outputs": [ { - "ename": "KeyboardInterrupt", - "evalue": "", - "output_type": "error", - "traceback": [ - "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m", - "\u001b[1;31mKeyboardInterrupt\u001b[0m Traceback (most recent call last)", - "\u001b[1;32my:\\KELAG\\KS\\KS-PW\\04 Digitalisierung\\KSPWDEV Server\\Digital Trainee Projekt\\DT_Slot_3_Project_Repo\\Druckrohrleitung\\Druckrohrleitung_test_steady_state.ipynb Cell 7\u001b[0m in \u001b[0;36m\u001b[1;34m()\u001b[0m\n\u001b[0;32m 40\u001b[0m fig1\u001b[39m.\u001b[39mcanvas\u001b[39m.\u001b[39mdraw()\n\u001b[0;32m 41\u001b[0m fig1\u001b[39m.\u001b[39mtight_layout()\n\u001b[1;32m---> 42\u001b[0m plt\u001b[39m.\u001b[39;49mpause(\u001b[39m0.000001\u001b[39;49m)\n", - "File \u001b[1;32mc:\\ProgramData\\Anaconda3\\envs\\Georg_DT_Slot3\\lib\\site-packages\\matplotlib\\pyplot.py:548\u001b[0m, in \u001b[0;36mpause\u001b[1;34m(interval)\u001b[0m\n\u001b[0;32m 546\u001b[0m canvas\u001b[39m.\u001b[39mdraw_idle()\n\u001b[0;32m 547\u001b[0m show(block\u001b[39m=\u001b[39m\u001b[39mFalse\u001b[39;00m)\n\u001b[1;32m--> 548\u001b[0m canvas\u001b[39m.\u001b[39;49mstart_event_loop(interval)\n\u001b[0;32m 549\u001b[0m \u001b[39melse\u001b[39;00m:\n\u001b[0;32m 550\u001b[0m time\u001b[39m.\u001b[39msleep(interval)\n", - "File \u001b[1;32mc:\\ProgramData\\Anaconda3\\envs\\Georg_DT_Slot3\\lib\\site-packages\\matplotlib\\backends\\backend_qt.py:409\u001b[0m, in \u001b[0;36mFigureCanvasQT.start_event_loop\u001b[1;34m(self, timeout)\u001b[0m\n\u001b[0;32m 405\u001b[0m timer \u001b[39m=\u001b[39m QtCore\u001b[39m.\u001b[39mQTimer\u001b[39m.\u001b[39msingleShot(\u001b[39mint\u001b[39m(timeout \u001b[39m*\u001b[39m \u001b[39m1000\u001b[39m),\n\u001b[0;32m 406\u001b[0m event_loop\u001b[39m.\u001b[39mquit)\n\u001b[0;32m 408\u001b[0m \u001b[39mwith\u001b[39;00m _maybe_allow_interrupt(event_loop):\n\u001b[1;32m--> 409\u001b[0m qt_compat\u001b[39m.\u001b[39m_exec(event_loop)\n", - "File \u001b[1;32mc:\\ProgramData\\Anaconda3\\envs\\Georg_DT_Slot3\\lib\\contextlib.py:120\u001b[0m, in \u001b[0;36m_GeneratorContextManager.__exit__\u001b[1;34m(self, type, value, traceback)\u001b[0m\n\u001b[0;32m 118\u001b[0m \u001b[39mif\u001b[39;00m \u001b[39mtype\u001b[39m \u001b[39mis\u001b[39;00m \u001b[39mNone\u001b[39;00m:\n\u001b[0;32m 119\u001b[0m \u001b[39mtry\u001b[39;00m:\n\u001b[1;32m--> 120\u001b[0m \u001b[39mnext\u001b[39;49m(\u001b[39mself\u001b[39;49m\u001b[39m.\u001b[39;49mgen)\n\u001b[0;32m 121\u001b[0m \u001b[39mexcept\u001b[39;00m \u001b[39mStopIteration\u001b[39;00m:\n\u001b[0;32m 122\u001b[0m \u001b[39mreturn\u001b[39;00m \u001b[39mFalse\u001b[39;00m\n", - "File \u001b[1;32mc:\\ProgramData\\Anaconda3\\envs\\Georg_DT_Slot3\\lib\\site-packages\\matplotlib\\backends\\qt_compat.py:262\u001b[0m, in \u001b[0;36m_maybe_allow_interrupt\u001b[1;34m(qapp)\u001b[0m\n\u001b[0;32m 260\u001b[0m signal\u001b[39m.\u001b[39msignal(signal\u001b[39m.\u001b[39mSIGINT, old_sigint_handler)\n\u001b[0;32m 261\u001b[0m \u001b[39mif\u001b[39;00m handler_args \u001b[39mis\u001b[39;00m \u001b[39mnot\u001b[39;00m \u001b[39mNone\u001b[39;00m:\n\u001b[1;32m--> 262\u001b[0m old_sigint_handler(\u001b[39m*\u001b[39;49mhandler_args)\n", - "\u001b[1;31mKeyboardInterrupt\u001b[0m: " + "name": "stdout", + "output_type": "stream", + "text": [ + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n", + "96532.05564888839\n" ] } ], @@ -180,7 +2184,6 @@ " # set initial conditions for the reservoir time evolution calculted with e-RK4\n", " V.set_pressure = p_old[0]\n", " # V.set_outflux = v_old[0]*area_pipe\n", - " print(V.get_current_pressure())\n", " # calculate the time evolution of the reservoir level within each pipeline timestep to avoid runaway numerical error\n", " for it_res in range(nt_eRK4):\n", " V.timestep_reservoir_evolution() \n", @@ -189,7 +2192,6 @@ " \n", " # set boundary conditions for the next timestep of the characteristic method\n", " p_boundary_res[it_pipe] = V.get_current_pressure()\n", - " print(V.get_current_pressure())\n", " v_boundary_tur[it_pipe] = initial_flux/area_pipe\n", "\n", " # the the boundary conditions in the pipe.object and thereby calculate boundary pressure at turbine\n", @@ -222,7 +2224,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 16, "metadata": {}, "outputs": [], "source": [ @@ -258,7 +2260,7 @@ ], "metadata": { "kernelspec": { - "display_name": "Python 3.8.13 ('Georg_DT_Slot3')", + "display_name": "Python 3.8.13 ('DT_Slot_3')", "language": "python", "name": "python3" }, @@ -277,7 +2279,7 @@ "orig_nbformat": 4, "vscode": { "interpreter": { - "hash": "84fb123bdc47ab647d3782661abcbe80fbb79236dd2f8adf4cef30e8755eb2cd" + "hash": "4a28055eb8a3160fa4c7e4fca69770c4e0a1add985300856aa3fcf4ce32a2c48" } } },