diff --git a/rocket_twin/drivers/mission.py b/rocket_twin/drivers/mission.py index 5544079..e9f54eb 100644 --- a/rocket_twin/drivers/mission.py +++ b/rocket_twin/drivers/mission.py @@ -43,11 +43,10 @@ def __init__( # Init and stop conditions init_fuel = init init_flight = { - "rocket.stage_1.tank.fuel.w_out_max": 3.0, - "g_tank.w_in": 0.0, + "rocket.controller.is_on_1": True, } - stop_fuel = "rocket.stage_1.tank.weight_prop >= rocket.stage_1.tank.weight_max" + stop_fuel = "rocket.flying == 1." stop_flight = stop # Fueling diff --git a/rocket_twin/notebooks/Demonstration.ipynb b/rocket_twin/notebooks/Demonstration.ipynb index 720c08a..319f2db 100644 --- a/rocket_twin/notebooks/Demonstration.ipynb +++ b/rocket_twin/notebooks/Demonstration.ipynb @@ -18,7 +18,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 1, "id": "e207e054-e7c1-462f-ba12-717bfdde929a", "metadata": {}, "outputs": [], @@ -46,7 +46,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 2, "id": "311d6d17-b7a0-41f9-b772-a532de531b4a", "metadata": {}, "outputs": [], @@ -76,10 +76,67 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 3, "id": "0b4dcf43-7c56-48f6-be1d-6c35042f3f5e", "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "\n", + " \n", + " " + ], + "text/markdown": [ + "\n", + "### Child components\n", + "\n", + "- `controller`: StationControllerCoSApp\n", + "- `g_tank`: Tank\n", + "- `pipe`: Pipe\n", + "- `rocket`: Rocket\n", + "\n", + "### Inputs\n", + "\n", + "- `inwards`: ExtensiblePort\n", + "\n", + "
\n", + "\n", + "| | |\n", + "---|---\n", + " **n_stages** 🔒🔒 : 3 | Number of stages\n", + " **fueling** 🔒🔒 : 1 | Whether the rocket is fueling or not\n", + " **time_int** 🔒🔒 : 10 s | Interval between fueling end and launch\n", + " **time_lnc** 🔒🔒 : 1e+05 s | Launch time\n", + "
\n", + "\n", + "### Outputs\n", + "\n", + "- `outwards`: ExtensiblePort\n", + "\n", + "
\n", + "\n", + "| | |\n", + "---|---\n", + " **stage**: 1 | Current stage\n", + "
" + ], + "text/plain": [ + "sys - Station" + ] + }, + "execution_count": 3, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "sys" ] @@ -114,7 +171,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 4, "id": "74278944-3801-40fc-971a-afee35419490", "metadata": {}, "outputs": [], @@ -129,7 +186,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 5, "id": "e7ff7cc9-b0a0-4485-88d3-911d52bd398c", "metadata": {}, "outputs": [], @@ -144,7 +201,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 6, "id": "62f1a722-25d0-4ed2-bc8f-fd332548b5ee", "metadata": {}, "outputs": [], @@ -159,7 +216,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 7, "id": "fe0f658e-dad8-4f28-ab76-eebaa9e4c54f", "metadata": {}, "outputs": [], @@ -174,7 +231,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 8, "id": "2918b401-bad6-4877-9442-f7ff6e6b70c4", "metadata": {}, "outputs": [], @@ -189,7 +246,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 9, "id": "390b25d1-a7fc-4154-9a7d-8dd0435c4e72", "metadata": {}, "outputs": [], @@ -206,7 +263,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 10, "id": "7a53323a-97a5-4a3c-94e3-be556381eb54", "metadata": {}, "outputs": [], @@ -223,7 +280,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 11, "id": "cb0aaa03-cd47-42cc-bcb8-069c501b3130", "metadata": {}, "outputs": [], @@ -256,7 +313,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "id": "07a77b38-c572-4dc3-a501-5fa697ab780a", "metadata": {}, "outputs": [], @@ -270,7 +327,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 13, "id": "4c091eb9-71a8-4e56-a778-d78655f6a419", "metadata": {}, "outputs": [], @@ -284,7 +341,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 14, "id": "ace3cc11-5ebb-4ae3-83e0-8f196a657fdd", "metadata": {}, "outputs": [], @@ -298,7 +355,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 15, "id": "71f7cecf-782c-48e6-a77c-d657d3195e64", "metadata": {}, "outputs": [], @@ -312,7 +369,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 16, "id": "ff07b0d0-7e88-4f98-9738-0ccff2466b02", "metadata": {}, "outputs": [], @@ -326,7 +383,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 17, "id": "adafef6b-cbcf-49d2-bdf5-55ba8fc646d9", "metadata": {}, "outputs": [], @@ -357,7 +414,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 18, "id": "8f0ae3db-00d8-498a-9026-fc60a7eaf5d3", "metadata": {}, "outputs": [], @@ -371,7 +428,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 19, "id": "ec3bed35-64fa-4ba1-a67b-7f1b14efffe0", "metadata": {}, "outputs": [], @@ -385,7 +442,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 20, "id": "8015a193-b30d-42c1-a74e-c2e6c0bfe02c", "metadata": {}, "outputs": [], @@ -399,7 +456,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 21, "id": "6009fc0d-123a-4f0a-a34a-f5520f362eab", "metadata": {}, "outputs": [], @@ -428,7 +485,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 22, "id": "c697aa0c-b08c-4291-9252-233ae6396f6f", "metadata": {}, "outputs": [], @@ -440,7 +497,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 23, "id": "3ff7653a-3869-4f38-bc48-04db10913f26", "metadata": {}, "outputs": [], @@ -452,7 +509,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 24, "id": "5321dac7-7b27-4b70-a40b-9b2b1d459312", "metadata": {}, "outputs": [], @@ -464,7 +521,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 25, "id": "a646f6ab-1d9e-4627-8f86-dcf38ba1a980", "metadata": {}, "outputs": [], @@ -495,7 +552,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 26, "id": "6f19b613-e749-48c6-b53c-d5ae4a909fec", "metadata": {}, "outputs": [], @@ -507,7 +564,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 27, "id": "1acf4530-c11e-4bde-a287-13db8887abb1", "metadata": {}, "outputs": [], @@ -519,7 +576,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 28, "id": "c828bce0-3c5c-42f2-a904-a672c7829053", "metadata": {}, "outputs": [], @@ -531,7 +588,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 29, "id": "f4d195e2-52cf-4e73-9c5e-503c1e1e35d5", "metadata": {}, "outputs": [], @@ -543,7 +600,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 30, "id": "0445c458-d6ca-4446-852c-1a4641913985", "metadata": {}, "outputs": [], @@ -555,7 +612,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 31, "id": "6e84cfdf-d918-4815-9f6a-f591baf024df", "metadata": {}, "outputs": [], @@ -567,7 +624,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 32, "id": "7741d88b-c6f1-4fa2-8be9-c81560efcedf", "metadata": {}, "outputs": [], @@ -579,7 +636,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 33, "id": "aebd952d-df10-4293-9ea5-2db4c44867e3", "metadata": {}, "outputs": [], @@ -601,12 +658,25 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 34, "id": "9fee58e0-11d3-4ff1-af73-a7384dfec86b", "metadata": {}, - "outputs": [], - "source": [ - "use_fmu = False\n", + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "Connector source 'stage_1.outwards.weight_prop' has physical unit kg, but target 'controller.inwards.weight_prop_1' is dimensionless.\n", + "Connector source 'stage_1.outwards.weight_prop' has physical unit kg, but target 'controller.inwards.weight_prop_1' is dimensionless.\n", + "Connector source 'stage_2.outwards.weight_prop' has physical unit kg, but target 'controller.inwards.weight_prop_2' is dimensionless.\n", + "Connector source 'stage_2.outwards.weight_prop' has physical unit kg, but target 'controller.inwards.weight_prop_2' is dimensionless.\n", + "Connector source 'stage_3.outwards.weight_prop' has physical unit kg, but target 'controller.inwards.weight_prop_3' is dimensionless.\n", + "Connector source 'stage_3.outwards.weight_prop' has physical unit kg, but target 'controller.inwards.weight_prop_3' is dimensionless.\n" + ] + } + ], + "source": [ + "use_fmu = True\n", "\n", "model_path = r\"systems\\control\\station_controller.mo\"\n", "model_name = \"station_controller\"\n", @@ -636,7 +706,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 35, "id": "c7370468-2245-498a-88f7-c86f8d1640c2", "metadata": {}, "outputs": [], @@ -678,14 +748,14 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 36, "id": "46800442-0c21-44fa-9483-0bcd7632d814", "metadata": {}, "outputs": [], "source": [ "# Time interval and step\n", "\n", - "man_sim = False\n", + "man_sim = True\n", "\n", "dt = 1.\n", "T = 50.\n", @@ -702,6 +772,7 @@ "driver = sys.add_driver(RungeKutta(order=4, dt=dt))\n", "driver.add_child(NonLinearSolver('solver'))\n", "driver.time_interval = (0, T)\n", + "driver.set_scenario(init=init)\n", "\n", "includes = [\"rocket.a\", \"rocket.geom.weight\", \"rocket.weight_prop_1\", \n", " \"rocket.weight_prop_2\", \"rocket.weight_prop_3\", \"g_tank.weight_prop\"]\n", @@ -731,7 +802,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 37, "id": "3fba5efd-0cd2-4021-ab22-298aa03bb26a", "metadata": {}, "outputs": [], @@ -771,12 +842,12 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 38, "id": "380230af-5678-418c-8ba6-6d0995173fa1", "metadata": {}, "outputs": [], "source": [ - "fly_sim = True\n", + "fly_sim = False\n", "\n", "dt = 1.\n", "\n", @@ -784,11 +855,12 @@ " \"rocket.stage_1.tank.fuel.weight_p\" : stage1_tank_mass,\n", " \"rocket.stage_2.tank.fuel.weight_p\" : stage2_tank_mass,\n", " \"rocket.stage_3.tank.fuel.weight_p\" : stage3_tank_mass,\n", + " \"rocket.controller.is_on_1\" : True,\n", " \"fueling\" : False,\n", " \"rocket.flying\" : True,\n", " }\n", "\n", - "stop = f\"rocket.stage_{n_stages}.tank.weight_prop < 0\"\n", + "stop = f\"rocket.stage_{n_stages}.tank.weight_prop <= 0.15\"\n", "\n", "includes = [\"rocket.a\", \"rocket.geom.weight\", \"rocket.weight_prop_1\", \n", " \"rocket.weight_prop_2\", \"rocket.weight_prop_3\", \"g_tank.weight_prop\"]\n", @@ -806,12 +878,12 @@ "id": "29b6343b-7ca7-4db4-acff-f33ff3bc745a", "metadata": {}, "source": [ - "### Mission (fueling + flying)" + "### Mission (fueling + flying) (bug avec le temps)" ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 39, "id": "0c5e8840-180f-430b-b37a-3592300355a5", "metadata": {}, "outputs": [], @@ -828,7 +900,7 @@ " \"rocket.flying\" : False,\n", " }\n", "\n", - "stop = f\"rocket.stage_{n_stages}.tank.weight_prop < 0\"\n", + "stop = f\"rocket.stage_{n_stages}.tank.weight_prop <= 0\"\n", "\n", "includes = [\"rocket.a\", \"rocket.geom.weight\", \"rocket.weight_prop_1\", \n", " \"rocket.weight_prop_2\", \"rocket.weight_prop_3\", \"g_tank.weight_prop\"]\n", @@ -851,81 +923,29 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 40, "id": "742acffd-710e-478e-bcf6-d6ad9a4a7fe5", "metadata": {}, "outputs": [], "source": [ - "seq_sim = False" - ] - }, - { - "cell_type": "markdown", - "id": "b9a7d713-d977-4bdb-9a4e-f876c423a4ec", - "metadata": {}, - "source": [ - "### Initialisation sequence" - ] - }, - { - "cell_type": "code", - "execution_count": null, - "id": "06f203b4-74b8-4198-92df-918bb0d6c262", - "metadata": {}, - "outputs": [], - "source": [ - "init = {\"g_tank.fuel.weight_p\" : 20.}\n", + "seq_sim = False\n", "\n", - "seq_init = [{\"name\" : \"start\", \"init\" : init, \"type\" : \"static\"}] \n", + "#Stop conditions\n", + "stop_fuel = f\"rocket.stage_{n_stages}.tank.weight_prop == rocket.stage_{n_stages}.tank.weight_max\"\n", + "stop_fly = f\"rocket.stage_{n_stages}.tank.weight_prop == 0.\"\n", "\n", - "if seq_sim:\n", - " run_sequences(sys, seq_init)" - ] - }, - { - "cell_type": "markdown", - "id": "8051e989-27cf-48c3-866a-f9e9ec366691", - "metadata": {}, - "source": [ - "### Fueling sequence" - ] - }, - { - "cell_type": "code", - "execution_count": null, - "id": "7c094628-fb44-4f43-9cb7-25299e40624c", - "metadata": {}, - "outputs": [], - "source": [ - "stop = f\"rocket.stage_{n_stages}.tank.weight_prop = rocket.stage_{n_stages}.tank.weight_max\"\n", + "#Sequences\n", + "seq_init = [{\"name\" : \"start\", \"init\" : {\"g_tank.fuel.weight_p\" : 20.}, \"type\" : \"static\"}] \n", + "seq_fuel = [{\"name\" : \"fuel\", \"type\" : \"transient\", \"dt\" : 1.0, \"stop\" : stop_fuel}]\n", + "seq_fly = [{\"name\" : \"flight\", \"type\" : \"transient\", \"dt\" : 1.0, \"stop\" : stop_fly}]\n", "\n", - "seq_fuel = [{\"name\" : \"fuel\", \"type\" : \"transient\", \"dt\" : 1.0, \"stop\" : stop}]\n", + "seq = [seq_init[0], seq_fuel[0], seq_fly[0]]\n", "\n", - "if seq_sim:\n", - " run_sequences(sys, seq_fuel)" - ] - }, - { - "cell_type": "markdown", - "id": "e71d863c-4b2a-4aae-982f-b46ee56bad86", - "metadata": {}, - "source": [ - "### Flying sequence" - ] - }, - { - "cell_type": "code", - "execution_count": null, - "id": "a7f0f296-a164-42c8-9251-fda27d9038ed", - "metadata": {}, - "outputs": [], - "source": [ - "stop = f\"rocket.stage_{n_stages}.tank.weight_prop == 0.\"\n", - "\n", - "seq_fly = [{\"name\" : \"flight\", \"type\" : \"transient\", \"dt\" : 1.0, \"stop\" : stop}]\n", + "includes = [\"rocket.a\", \"rocket.geom.weight\", \"rocket.weight_prop_1\", \n", + " \"rocket.weight_prop_2\", \"rocket.weight_prop_3\", \"g_tank.weight_prop\"]\n", "\n", "if seq_sim:\n", - " run_sequences(sys, seq_fly)\n", + " run_sequences(sys, seq, includes)\n", " data = sys.drivers['rk'].recorder.export_data()" ] }, @@ -939,7 +959,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 41, "id": "b14eb1c5-c5a8-4a0a-a3a5-25bdf048b163", "metadata": {}, "outputs": [], @@ -964,10 +984,21 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 42, "id": "be8b8507-a9ff-4643-8fe0-c2d488affedb", "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjMAAAHFCAYAAAAHcXhbAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABzIklEQVR4nO3dd3gU9drG8e+k94T0BAIJVTpIE5SuSFMQVBBp4lGRckBEj8hRwaOgvAqK2BWOHT1KUSmCAgGkCCISQJGSUBMCAdLr7rx/BFYjLYEkmw3357r2MlP32ZGwNzPzzM8wTdNERERExEE52bsAERERkauhMCMiIiIOTWFGREREHJrCjIiIiDg0hRkRERFxaAozIiIi4tAUZkRERMShKcyIiIiIQ1OYEREREYemMCNSCfz3v//FMAwMw2DNmjXnLTdNk9q1a2MYBp06dSr3+iqLTz/9lFdeeeWCywzDYMqUKeVaj4gUUpgRqUR8fX15//33z5sfGxvL/v378fX1tUNVlcelwszGjRv5xz/+Ub4FiQigMCNSqQwYMICvvvqKtLS0IvPff/992rZtS/Xq1e1UmePIzs6+ou1uuOEGqlWrVsrViEhxKMyIVCL33HMPAJ999pltXmpqKl999RUjRoy44DZTp06lTZs2BAYG4ufnx/XXX8/777/P38egXbVqFZ06dSIoKAhPT0+qV69O//79ycrKsq3z5ptv0rRpU3x8fPD19eW6667jySefvGzdp06dYtSoUVStWhU3Nzdq1qzJ5MmTyc3Nta3TvHlz2rdvf962FouFqlWr0q9fP9u8vLw8nnvuOa677jrc3d0JCQnhvvvu48SJE0W2jY6Opnfv3ixYsIDmzZvj4eHB1KlTL1hjp06dWLJkCQcPHrRd0jMMw7b875eZzl36W7VqFQ888ABBQUH4+fkxdOhQMjMzSUpK4u677yYgIICIiAgmTpxIfn5+kfcs7ucQuda52LsAESk9fn5+3HnnncydO5eHHnoIKAw2Tk5ODBgw4IKXSBISEnjooYdsZ202bdrE2LFjOXr0KE8//bRtnV69etG+fXvmzp1LQEAAR48eZfny5eTl5eHl5cX8+fMZNWoUY8eO5aWXXsLJyYl9+/axe/fuS9ack5ND586d2b9/P1OnTqVJkyasW7eO6dOns337dpYsWQLAfffdx7hx49i7dy916tSxbb9ixQqOHTvGfffdB4DVaqVPnz6sW7eOxx9/nHbt2nHw4EGeeeYZOnXqxNatW/H09LRtv23bNn777Tf+/e9/ExMTg7e39wXrfOONN3jwwQfZv38/CxcuLOb/EfjHP/5Bv379mD9/Pr/88gtPPvkkBQUF7Nmzh379+vHggw/y/fff8+KLLxIZGcmECROu6HOIXNNMEXF48+bNMwFzy5Yt5urVq03A3Llzp2maptmqVStz+PDhpmmaZsOGDc2OHTtedD8Wi8XMz883n332WTMoKMi0Wq2maZrml19+aQLm9u3bL7rtmDFjzICAgBLX/tZbb5mA+cUXXxSZ/+KLL5qAuWLFCtM0TfPkyZOmm5ub+eSTTxZZ7+677zbDwsLM/Px80zRN87PPPjMB86uvviqy3pYtW0zAfOONN2zzatSoYTo7O5t79uwpVq29evUya9SoccFlgPnMM8/Yps/9Pxk7dmyR9fr27WsC5syZM4vMb9asmXn99dfbpkvyOUSudbrMJFLJdOzYkVq1ajF37lzi4uLYsmXLRS8xQeHlo5tvvhl/f3+cnZ1xdXXl6aefJiUlheTkZACaNWuGm5sbDz74IB988AEHDhw4bz+tW7fmzJkz3HPPPSxevJiTJ08Wq95Vq1bh7e3NnXfeWWT+8OHDAfjhhx8ACAoK4rbbbuODDz7AarUCcPr0aRYvXszQoUNxcSk80fztt98SEBDAbbfdRkFBge3VrFkzwsPDz+v2atKkCXXr1i1WrVeid+/eRabr168PQK9evc6bf/DgQdt0ST+HyLVMYUakkjEMg/vuu4+PP/6Yt956i7p1617wXhOAn376iW7dugHw7rvv8uOPP7JlyxYmT54M/HkzbK1atfj+++8JDQ1l9OjR1KpVi1q1avHqq6/a9jVkyBDmzp3LwYMH6d+/P6GhobRp04aVK1dest6UlBTCw8OL3H8CEBoaiouLCykpKbZ5I0aM4OjRo7Z9fvbZZ+Tm5tqCD8Dx48c5c+YMbm5uuLq6FnklJSWdF7IiIiIuWd/VCgwMLDLt5uZ20fk5OTm26ZJ+DpFrme6ZEamEhg8fztNPP81bb73F888/f9H15s+fj6urK99++y0eHh62+YsWLTpv3fbt29O+fXssFgtbt27ltddeY/z48YSFhTFw4ECg8L6W++67j8zMTNauXcszzzxD7969+eOPP6hRo8YFawgKCmLz5s2Yplkk0CQnJ1NQUEBwcLBt3q233kpkZCTz5s3j1ltvZd68ebRp04YGDRrY1gkODiYoKIjly5df8P3+3p7+9xBVUZT0c4hcyxRmRCqhqlWr8thjj/H7778zbNiwi65nGAYuLi44Ozvb5mVnZ/PRRx9ddBtnZ2fatGnDddddxyeffMK2bdtsYeYcb29vevToQV5eHn379mXXrl0XDTNdu3bliy++YNGiRdxxxx22+R9++KFt+V/fe8iQIbzyyiusW7eOrVu38vbbbxfZX+/evZk/fz4Wi4U2bdpc9HNcCXd39ytu3S6psvwcIpWNwoxIJfXCCy9cdp1evXoxc+ZMBg0axIMPPkhKSgovvfQS7u7uRdZ76623WLVqFb169aJ69erk5OQwd+5cAG6++WYAHnjgATw9PbnxxhuJiIggKSmJ6dOn4+/vT6tWrS5aw9ChQ3n99dcZNmwYCQkJNG7cmPXr1zNt2jR69uxp2/85I0aM4MUXX2TQoEF4enoyYMCAIssHDhzIJ598Qs+ePRk3bhytW7fG1dWVI0eOsHr1avr06VMkNJVE48aNWbBgAW+++SYtWrTAycmJli1bXtG+LqcsP4dIZaMwI3IN69KlC3PnzuXFF1/ktttuo2rVqjzwwAOEhoZy//3329Zr1qwZK1as4JlnniEpKQkfHx8aNWrE119/bbvnpn379vz3v//liy++4PTp0wQHB3PTTTfx4YcfEhISctEaPDw8WL16NZMnT+b//u//OHHiBFWrVmXixIk888wz561ft25d2rVrx4YNG7j33nvx9/cvstzZ2Zmvv/6aV199lY8++ojp06fj4uJCtWrV6NixI40bN77i4zVu3Dh27drFk08+SWpqKqZpnvc8ntJSlp9DpLIxzLL6TRQREREpB+pmEhEREYemMCMiIiIOTWFGREREHJrCjIiIiDg0hRkRERFxaAozIiIi4tAq/XNmrFYrx44dw9fXt8I+tlxERESKMk2T9PR0IiMjcXK69LmXSh9mjh07RlRUlL3LEBERkStw+PBhqlWrdsl1Kn2YOTcY2+HDh/Hz8yvdnaekwOuvg4cHeHmV7r5FiisrC3JyYPRoCAqydzWVVkpWCq//9Doerh54uer3XQQgKz+LnPwcRrceTZBX6f79k5aWRlRUVLEGVa30YebcpSU/P7/SDzN5eeDuDoGB4ONTuvsWKa6MjMJg7esLpf1nXGzyXPJw93Yn0CsQHzf9vosAZORlkJKVgq+fL35eZfP3T3FuEdENwCIiIuLQFGZERETEoSnMiIiIiEOr9PfMiIiIfTnhhCuu9i5DyoAFCz5OPuTn5pPjlFOibV1dXXF2di6VOhRmRESkzIQTTqhTKE6GE+hRX5WO6Wxi9bCSciyFM05nSrx9QEAA4eHhV/0cOIUZEREpE+GEE+kaSVBwEG7ubnpwaSVkMS1YTSuBnoG4Ohf/7JtpmmRlZZGcnAxARETEVdWhMCMiIqXOCSdCnUIJCg7Cx1+t7JWVk+mExWrBw8OjRGEGwNPTE4Dk5GRCQ0Ov6pKTbgAWEZFS54orToYTbu5u9i5FKjCvsw+czc/Pv6r9KMyIiEjZMIr3wDO5dpXWnw+FGREREXFoCjMiIiLXoKysLO4bdB81I2oS6hNK6pnUUtnv8OHD6du3b6nsq7gUZkRERP7iRPIJHh37KM2va061wGo0rNmQu/vczZbNW2zrhPqEsvSbpXassqhZM2bRs2tPaoTUoHbV2sXa5vNPPmfzhs0s+X4Jcfvj8PN33LHd1M10FRL2pBFlGpTOI39ERKQiGHHvCAoKCnjt7deoEVODE8knWLdmHWdOn7F3aReVl5fH7XfcTsvWLfn0w0+LtU1CfAJ16tWhfsP6ZVxd2dOZmSu0+esDLPkogR8z69i7FBERKSWpZ1LZvHEzTz37FDd1vImo6lFc3/J6xk0cxy3dbwGgRYMWAAy/ZzihPqG26fgD8QwdMJQGMQ2IDoumW4duxK6OLbL/40nHGdR/ENWDq9OyYUu++uIrWjRowduvv21bJy01jUfHPEqD6AbUjKhJv5792Bm385J1/+vf/2LkmJE0aNigWJ+zb/e+vDn7TTb+uJFQn1D6du8LXPiMU+2qtZn/8XzbdOKxRB4Y+gB1qtWhQfUGjLhnBAkJCcV637KiMHOFQqr7AhCXG0XcmSA7VyMiUvGZpkl2nsUuL9M0i1Wjt4833j7eLPt2Gbm5uRdc57vY7wCY/dZs4vbH2aYzMzPp2q0rX37zJat+XEXnrp0ZctcQjhw+Ytt29AOjSUpMYuGyhcz9ZC4fzf2IkydOFjlGg+4cRHJyMp8u+JTv131P42aNubPXnZw+dfpKD/155n06jyHDh9CyTUvi9scx79N5xdouKyuLO3regbePN4uXL2bRikV4eXvRu1dv8vLySq2+ktJlpitUs1kIN9wSzqaVSaxLrkaA3zGi/LLtXZaISIWVk2+lw4s/2+W91/6rBZ5ul78pwMXFhdlvzebRsY/ywfsf0LhZY9rd1I6+d/alYaOGAASHBAPg5+9HWFiYbdtGjRvRqHEj2/SkZyax9JulfLfkO+4feT979+xl7eq1rFi7gmbXNwNg5uszuaHpDbZt1seu57ddv7E7fjfu7u4ATJ02lWXfLuObRd8wdMTQqz4WAFUCq+Dp5Ymbq1uRz3A5i75chJPhxKzXZ2EYBhbTwsw3ZtKwekPWrFlDt27dSqW+klKYuQrXdwjhdOwv7MmL4Lv4cO6sd4QAj6t78I+IiNjXbX1v45but7Bpwya2bt7Kqu9XMWfWHGa9PouBgwdedLvMzExemv4SK5etJCkpiYKCAnKyczhypPDMzL69+3BxcaFJsya2bWrWqklAlQDb9K/bfyUzI5N61esV2XdOdg4J8Qml+jmvxK+//Er8gXhiwmOKzM/JyWH//v12qkph5qoYhkEnn99JzfIjKcebb/dHcGe9I3i4WO1dmohIhePh6sTaf7Ww23uXaH0PDzp16USnLp2YOGkij4x+hBnPz7hkmJk6eSqrf1jNlOenEFMzBg9PD+4ffD/5eYX/yL3Ypa6/zrdarYSFh7Fw2cLz1vP39y/RZ7gShmGcV2dBfkGR+po2b8ob779ROG1asVgtBHkFERkeWeb1XYzCzFVyMaz0iIznf0fqkZrrxnfx4fSufQxnPfRSRKQIwzCKdamnIqp7XV2WfbPMNu3q6orFYimyzuYNmxl470B63d4LgIyMDA4fOmxbXqduHQoKCoj7NY6mzZsCcGD/gSLPd2nSrAnJx5NxcXGheo3qZfmRLigoOIjjScdt0wf2HSArK6tIfYsXLCYkJARfP18spgWL1UKod2iJx2YqTboBuBR4uRTQq2YiLk5WjqR7sf5wsL1LEhGRK3Aq5RT9evbjf/P/x66duziYcJCvF3zNnFlz6N67u229qBpRrFuzjuPHj9tatqNrRrPk6yXE7YhjZ9xOHh7xMFbrn2fq69SrQ4fOHXh07KNs27qNuF/jmDh2Ip6enhgU/gu4Y+eOtGzdkmEDh7Hq+1UcOniInzb9xPSp09m+bftF6z5y+AhxO+I4cvgIFouFuB1xxO2IIyMjo0Sfv33H9sx9ey47tu9g+7btPDbuMVxd/wwp/Qf0JzAokKEDhrLpx00cSjjExvUbmfDIBNvlNHtQmCklwV55dItOAkx2ngwg7kTZnw4UEZHS5e3jzfWtruftOW/T59Y+dGzdkRf+8wJDhg9h+svTbetNnTaV2NWxNK/XnK43dgXgPy/+h4CAAHp37c2Qu4bQqWsnmjRtUmT/r7/7OiGhIfS5tQ/D7xnO4OGD8fbxxt2j8GZfwzD4bMFntL2xLeMfHk/bZm15aPhDHDp0iJDQkIvW/eJzL9K1XVdmPD+DzIxMurbrStd2Xfl1268l+vxTp08lslokt996OyNHjGTUuFF4ennalnt5ebH4u8VUjarKfYPuo0PLDkwcPZHs7Gz8/Oz30D3DLG6/moNKS0vD39+f1NTU0j/QJ0/CrFkQFAQ+hUPcb0sKYOOxYAxMbqutDicpBxkZkJICjzwCwTorWFZOZp1k1sZZBHkF4ePmY+9yKjx33KnvUp9qUdVwdbff5YeK7tjRYzSr14wvv/mSDp072LucErvay0w5OTnEx8cTExODh4dHkWUl+f6265mZN998kyZNmuDn54efnx9t27Zl2bI/r0mapsmUKVOIjIzE09OTTp06sWvXLjtWfHnNw85QLzANE4Pl8eGcztEvsYiIFFq3Zh3LlyznYMJBftr0Ew8Of5DqNarT9qa29i7Nodk1zFSrVo0XXniBrVu3snXrVrp06UKfPn1sgWXGjBnMnDmTOXPmsGXLFsLDw7nllltIT0+3Z9mXZBjQuXoyEd7Z5FmcWbI/gpwCXc0TERHIL8hn2pRpdGjVgeH3DCcoKIiFyxYWuS9FSs6u3Uy33XZbkennn3+eN998k02bNtGgQQNeeeUVJk+eTL9+/QD44IMPCAsL49NPP+Whhx6yR8nF4uwE3Wsm8eWeaqTmurE8Ppzb1OEkInLN63JzF7rc3MXeZVQ6FeaUgcViYf78+WRmZtK2bVvi4+NJSkoq8jRBd3d3OnbsyIYNG+xYafF4uVroVSsRVycrR9O9WHc4hMp9d5KIiIh92D3MxMXF4ePjg7u7OyNHjmThwoU0aNCApKQkgPMesxwWFmZbdiG5ubmkpaUVedlLkGcet5ztcNp10l8dTiIiImXA7mGmXr16bN++nU2bNvHwww8zbNgwdu/ebVtuGEWvzZimed68v5o+fTr+/v62V1RUVJnVXhwxAVm0rZoCwPojwRxK87zMFiIiIlISdg8zbm5u1K5dm5YtWzJ9+nSaNm3Kq6++Snh4OMB5Z2GSk5MvOSjWpEmTSE1Ntb0OHz580XXLS/PQM1x3tsPpO3U4iYiIlCq7h5m/M02T3NxcYmJiCA8PZ+XKlbZleXl5xMbG0q5du4tu7+7ubmv1PveyN8OATn/pcPpWHU4iIiKlxq7dTE8++SQ9evQgKiqK9PR05s+fz5o1a1i+fDmGYTB+/HimTZtGnTp1qFOnDtOmTcPLy4tBgwbZs+wr4uwEPWom8r89UaTlurH8wNkOJ2UaERGRq2LXMHP8+HGGDBlCYmIi/v7+NGnShOXLl3PLLbcA8Pjjj5Odnc2oUaM4ffo0bdq0YcWKFfj6+tqz7Cvm6WqlV61EvtpTjaMZXqw9EkKnqBNc4hYgERGRMpGVlcXof4wmdnUsGekZ7D2yF/+Aq29UGT58OGfOnGHRokVXX2Qx2fW8wPvvv09CQgK5ubkkJyfz/fff24IMFN78O2XKFBITE8nJySE2NpZGjRrZseKrF+SZR7eYwg6n3Sf92aEOJxGRCuVE8gkeHfsoza9rTrXAajSs2ZC7+9zNls1bbOuE+oSy9JuldqzyT4cOHmL8qPG0bNiS6sHVadW4FS8+9yJ5eXmX3O7zTz5n84bNLPl+CXH74/Dzt/9tGVfKrmdmrlXR/lm0q5rChqPB/HgkmAD3fGr4Z11+QxERKXMj7h1BQUEBr739GjVianAi+QTr1qyzjY5d0ez7Yx9Wq5X/m/1/xNSM4ffdvzNhzASysrKYOm3qRbdLiE+gTr061G9YvxyrLRu6Y8NOmoWe4bqgwg6nFfFhnMpWh5OIiL2lnkll88bNPPXsU9zU8SaiqkdxfcvrGTdxHLd0L7xy0KJBCwCG3zOcUJ9Q23T8gXiGDhhKg5gGRIdF061DN2JXxxbZ//Gk4wzqP4jqwdVp2bAlX33xFS0atODt19+2rZOWmsajYx6lQXQDakbUpF/PfuyM23nRmrvc0oXZb82mc9fORMdE071Xd0aNG8WSr5dcdJu+3fvy5uw32fjjRkJ9QunbvS9w4TNOtavWZv7H823TiccSeWDoA9SpVocG1Rsw4p4RJCQkXP7gliGFGTsxDOgUdbbDyerMkv2R6nASkcrNNCEvyz6vYj6C3dvHG28fb5Z9u4zc3NwLrvNd7HcAzH5rNnH742zTmZmZdO3WlS+/+ZJVP66ic9fODLlrCEcOH7FtO/qB0SQlJrFw2ULmfjKXj+Z+xMkTJ/9yiEwG3TmI5ORkPl3wKd+v+57GzRpzZ687OX3qdLEPdXpqOlWqVLno8nmfzmPI8CG0bNOSuP1xzPt0XrH2m5WVxR0978Dbx5vFyxezaMUivLy96N2r92Uva5UlXWayoyIdTnmuLDsQwe21j6rDSUQqp/xsQl5tYpe3PjFuB7h5XXY9FxcXZr81m0fHPsoH739A42aNaXdTO/re2ZeGjRoCEBwSDICfv1+R5541atyIRo3/vK9z0jOTWPrNUr5b8h33j7yfvXv2snb1WlasXUGz65sBMPP1mdzQ9AbbNutj1/Pbrt/YHb8bd3d3AKZOm8qyb5fxzaJvGDpi6GU/Q/yBeN57+71LXmKqElgFTy9P3FzdLvnstr9b9OUinAwnZr0+C8MwsJgWZr4xk4bVG7JmzZoiQxCVJ4UZO/trh9OxDE/WHg6hU3V1OImI2MttfW/jlu63sGnDJrZu3sqq71cxZ9YcZr0+i4GDB150u8zMTF6a/hIrl60kKSmJgoICcrJzOHKk8MzMvr37cHFxoUmzPwNdzVo1CagSYJv+dfuvZGZkUq96vSL7zsnOISE+4bK1JyUmMfCOgdx+x+0MHj64ZB+8GH795VfiD8QTEx5TtL6cHPbv31/q71dcCjMVwLkOpyX7I9id4k+gZx5NQ1PtXZaISOly9Sw8Q2Kn9y4JDw8POnXpRKcunZg4aSKPjH6EGc/PuGSYmTp5Kqt/WM2U56cQUzMGD08P7h98P/l5+UDhJaQL+et8q9VKWHgYC5ctPG89f/9Ld78mJSZxR887aNm6JS+/9nJxPuZ5DMM4r86C/IIi9TVt3pQ33n+jcNq0YrFaCPIKIjI88oreszQozFQQ0f5Z3Fg1hR/V4SQilZVhFOtST0VU97q6LPtmmW3a1dUVi8VSZJ3NGzYz8N6B9Lq9FwAZGRkcPvTnkDp16tahoKCAuF/jaNq8KQAH9h8g9cyf/3ht0qwJyceTcXFxoXqN6sWuL/FYInf0vIOmzZoy+63ZODld2f0KQcFBHE86bps+sO8AWVl/fhc1adaExQsWExISgq+fLxbTgsVqIdQ7FFdn+zWy6O6MCqRp6Bnqq8NJRMRuTqWcol/Pfvxv/v/YtXMXBxMO8vWCr5kzaw7de3e3rRdVI4p1a9Zx/PhxW8t2dM1olny9hLgdceyM28nDIx7GarXatqlTrw4dOnfg0bGPsm3rNuJ+jWPi2Il4enpiUHhvQcfOHWnZuiXDBg5j1ferOHTwED9t+onpU6ezfdv2C9aclJhE3x59qVq1KlOmTeHkyZMcP36c48ePX3D9S2nfsT1z357Lju072L5tO4+NewxX1z+/i/oP6E9gUCBDBwxl04+bOJRwiI3rNzLhkQm2y2n2oDBTgRgGdIxKJtLnzw6nbHU4iYiUG28fb65vdT1vz3mbPrf2oWPrjrzwnxcYMnwI01+ebltv6rSpxK6OpXm95nS9sSsA/3nxPwQEBNC7a2+G3DWETl070aRp0RueX3/3dUJCQ+hzax+G3zOcwcMH4+3jjbtH4c2+hmHw2YLPaHtjW8Y/PJ62zdry0PCHOHToECGhIResec0Pa4jfH8+62HU0rduUxrUa214lNXX6VCKrRXL7rbczcsRIRo0bhafXn5fovLy8WPzdYqpGVeW+QffRoWUHJo6eSHZ2tl3HQjTMi13EqyTS0tLw9/cnNTW19A/0yZMwaxYEBYGPT6ntNrvAiS9/L+xwivTJVoeTXFpGBqSkwCOPQHCwvauptE5mnWTWxlkEeQXh41Z6v++VlTvu1HepT7Woari66yzzxRw7eoxm9Zrx5Tdf0qFzB3uXU2JXe5kpJyeH+Ph4YmJi8PDwKLKsJN/f+oqsgDxdrPSslYirk5VjGZ7EHg4p7iMSRESkAlu3Zh3LlyznYMJBftr0Ew8Of5DqNarT9qa29i7NoSnMVFBBnnncGpOEgclvKf78mqwxnEREHF1+QT7TpkyjQ6sODL9nOEFBQSxctrDIfSlScupmqsBq+GdxY7WTrD8SwoajwQR45BOtDicREYfV5eYudLm5i73LqHR0ZqaCaxKSSoOg1LMdTuGkZLvZuyQREZEKRWGmgjMM6BB1gkifbPKtTizdH6EOJxERkb/Qt6IDODeGk59bvm0MJ4v18tuJiIhcCxRmHISHi5VetY7h5mQhMcOT2MOh6nASERFBYcahBHrm0y3m+NkOJz9+TQ6wd0kiIiJ2pzDjYM51OAH8eDSIhFTHHOdERESktCjMOKAmIak0CE4FdTiJiMgVysrK4r5B91EzoiahPqFFBry8GsOHD6dv376lsq/iUphxQOc6nKr6ZJFvdWLJ/giy8/W/UkSkNJxIPsGjYx+l+XXNqRZYjYY1G3J3n7vZsnmLbZ1Qn1CWfrPUjlUWNeTuITS/rjlRQVE0qtWIUf8YRVJi0iW3+fyTz9m8YTNLvl9C3P44/PztN7bS1dI3oINyNqB7zST83PNIV4eTiEipGXHvCHbv3M1rb7/Gxu0b+eiLj7ix/Y220bErohs73Mi7H77Lhl82MPeTuSTEJzBi8IhLbpMQn0CdenWo37A+YWFhGIZRTtWWPoUZB+bhYqV3rUTcnC0kZnqy5pA6nERErkbqmVQ2b9zMU88+xU0dbyKqehTXt7yecRPHcUv3WwBo0aAFAMPvGU6oT6htOv5APEMHDKVBTAOiw6Lp1qEbsatji+z/eNJxBvUfRPXg6rRs2JKvvviKFg1a8Pbrb9vWSUtN49Exj9IgugE1I2rSr2c/dsbtvGTdI8eMpGXrlkRVj6L1Da3554R/8vNPP5Ofn3/B9ft278ubs99k448bCfUJpW/3vsCFzzjVrlqb+R/Pt00nHkvkgaEPUKdaHRpUb8CIe0aQkJBw+YNbhhRmHFwVj3zbGE6/n/LjF3U4iUgFZZom2QXZdnmZxfyXnrePN94+3iz7dhm5ubkXXOe72O8AmP3WbOL2x9mmMzMz6dqtK19+8yWrflxF566dGXLXEI4cPmLbdvQDo0lKTGLhsoXM/WQuH839iJMnThY5RoPuHERycjKfLviU79d9T+Nmjbmz152cPnW6WJ/h9KnTfPX5V7S6odVFx3ya9+k8hgwfQss2LYnbH8e8T+cVa99ZWVnc0fMOvH28Wbx8MYtWLMLL24vevXqTl5dXrH2UBY3NVAlU98vmpmonWXckhI1Hg6jinkdMgMZwEpGKJceSQ69Fvezy3kv6LsHTxfOy67m4uDD7rdk8OvZRPnj/Axo3a0y7m9rR986+NGzUEIDgkGAA/Pz9CAsLs23bqHEjGjVuZJue9Mwkln6zlO+WfMf9I+9n7569rF29lhVrV9Ds+mYAzHx9Jjc0vcG2zfrY9fy26zd2x+/G3d0dgKnTprLs22V8s+gbho4YetHan33qWea+PZesrCxatG7BJ//75KLrVgmsgqeXJ26ubkU+w+Us+nIRToYTs16fhWEYWEwLM9+YScPqDVmzZg3dunUr9r5Kk87MVBKNQ1JpeLbDaWWCOpxERK7UbX1vY8feHXz4xYd07tqZH9f9yM033lzkUsuFZGZmMvXfU7mpxU3Urlqb6LBo9v6xlyNHCs/M7Nu7DxcXF5o0a2LbpmatmgRUCbBN/7r9VzIzMqlXvR7RYdG216GEQyTEJ1zy/UePG80PP/7AF19/gbOzM2MeHFPsM1LF9esvvxJ/IJ6Y8Biiw6KpFV6LRtGNyMnJYf/+/aX6XiWhMzOVhGFA+6gTnMl15Wi6F0v2R3BnvSN4uVrsXZqICAAezh4s6bvEbu9dovU9POjUpROdunRi4qSJPDL6EWY8P4OBgwdedJupk6ey+ofVTHl+CjE1Y/Dw9OD+wfeTn1d438rFgsVf51utVsLCw1i4bOF56/n7+1+y5qDgIIKCg6hVpxZ169WlWb1mbP1pK63atCrORwbAMIzz6izILyhSX9PmTXnj/TcKp00rFquFIK8gIsMji/0+pU1hphJxNqB7TBJf7qlGaq4byw+E06fOUZx1/k1EKgDDMIp1qaciqntdXZZ9s8w27erqisVS9B+LmzdsZuC9A+l1e+GltIyMDA4fOmxbXqduHQoKCoj7NY6mzZsCcGD/gSLPd2nSrAnJx5NxcXGheo3qV1zvuUCSl1uy+1iCgoM4nnTcNn1g3wGysv68baFJsyYsXrCYkJAQfP18sZgWLFYLod6huDpf+P6c8qCvuUqmcAynPzucVqvDSUSk2E6lnKJfz378b/7/2LVzFwcTDvL1gq+ZM2sO3Xt3t60XVSOKdWvWcfz4cVvLdnTNaJZ8vYS4HXHsjNvJwyMexmr985kZderVoUPnDjw69lG2bd1G3K9xTBw7EU9PTwwK26I7du5Iy9YtGTZwGKu+X8Whg4f4adNPTJ86ne3btl+w5m1bt/H+W+8TtyOOw4cOsz52PSPvG0l0zWhatmlZos/fvmN75r49lx3bd7B923YeG/dYkZuI+w/oT2BQIEMHDGXTj5s4lHCIjes3MuGRCbbLafagMFMJVfHIp/vZDqc9p/z45XiAvUsSEXEI3j7eXN/qet6e8zZ9bu1Dx9YdeeE/LzBk+BCmvzzdtt7UaVOJXR1L83rN6XpjVwD+8+J/CAgIoHfX3gy5awidunaiSdMmRfb/+ruvExIaQp9b+zD8nuEMHj4Ybx9v3D0Kb/Y1DIPPFnxG2xvbMv7h8bRt1paHhj/EoUOHCAkNuWDNHh4eLPl6CXf2vpN2zdsxftR4rmtwHYuXL7bdRFxcU6dPJbJaJLffejsjR4xk1LhReHr9eTbNy8uLxd8tpmpUVe4bdB8dWnZg4uiJZGdn4+dnv4fuGWZp3x1UwaSlpeHv709qamrpH+iTJ2HWLAgKAh+f0t13KdiR7M+6IyGASc+aiepwqqwyMiAlBR55BIKD7V1NpXUy6ySzNs4iyCsIH7eK9/te0bjjTn2X+lSLqoaru/0uP1R0x44eo1m9Znz5zZd06NzB3uWU2NVeZsrJySE+Pp6YmBg8PIre11SS72+dmanEGoek0ujcGE4J4ZzMUoeTiIg9rVuzjuVLlnMw4SA/bfqJB4c/SPUa1Wl7U1t7l+bQFGYqMcOAm6JOUM03iwKrE0sORJCV72zvskRErln5BflMmzKNDq06MPye4QQFBbFw2cKLPtxOikfdTJWcswG3/qXDadmBcPqqw0lExC663NyFLjd3sXcZlY6+0q4B5zqc3J0tJKnDSUREKhmFmWvEuTGcAPac8iMxs2QPkBIREamoFGauIVF+2YR65QDo3hkREak0FGauMS5Our4kIiKVi8KMiIiIODSFGREREXFoCjMiIiKVxIznZ9C5bedLrrN542Y6tu5IZEAkQwcOLbX3NgyDRYsWldr+SkLPmREREfmb48ePM/ul2az8biWJRxPx9fOlZu2a3DngTu4edDdeXl72LvGKPTPpGRo2achnCz/D29vb3uWUCoUZERGRv0iIT6D3zb3x9/dn8jOTqd+wPgUFBezft5/PPvqM8IhwuvfqfsFt8/PzK/zTfBMOJDDs/mFEVo20dymlxq6XmaZPn06rVq3w9fUlNDSUvn37smfPniLrDB8+HMMwirxuuOEGO1UsIiKV3b8e+RcuLi6sWLeCPv37UPe6ujRo1IDb+t7Gp199yq09b7WtG+oTyn/f+y9DBwwlOjSamS/OBGDeu/No1bgVVatUpW3ztnzx2Re2bQ4dPESoTyhxO+Js81LPpBLqE8qPa38E4Me1PxLqE8ra1Wu5pf0t1AipQc+uPdn3x74itc5+eTYNYhoQEx7D+FHjyc3NvejnOve+p06dYtzD4wj1CWX+x/OZ//F8aletXWTdpd8sJdQntMi875Z+x8033UxUUBQtG7Xk/6b9HwUFBSU8umXDrmEmNjaW0aNHs2nTJlauXElBQQHdunUjMzOzyHrdu3cnMTHR9lq6dKmdKhYRkStlmiZmdrZ9XsV87PmplFOs+WENIx4ccdFLMIZhFJme8fwMuvfqzprNaxg0dBBLvl7Cvx//Nw+PfZi1P61l6IihjBs5jvWx60t8zKY/O52p06ayYu0KXJxdGDdqnG3Z4q8WM+P5GTz5zJOsXLuSsLAw5r0776L7qlqtKnH74/D18+W5F58jbn8cffr3KVYdq75fxah/jOIfD/+DdVvX8dLsl/j8k8959f9eLfFnKgt2vcy0fPnyItPz5s0jNDSUn3/+mQ4d/hwK3d3dnfDw8PIuT0RESlNODik33Xr59cpA0PrvwNPzsuvFH4jHNE1q1alVZP511a8jJ7fwoaMjHhzB0/952ras/939GTR0kG165IiRDLx3ICMeHAHAw3Ue5ueffuaN2W9wU8ebSlT3pKcn0a59OwD++eg/GdR/EDk5OXh4ePD2G29zz5B7GDx8cOG6z0xi7Zq15OTkXHBfzs7OhIWFYRgGfv5+hIWFFbuOV2a8wj8n/JOB9w4EIDommn899S+e/fezjP/X+BJ9prJQobqZUlNTAQgMDCwyf82aNYSGhlK3bl0eeOABkpOT7VGeiIhcI/5+9mV57HJWbVjFdfWvIy83r8iyptc3LTK9d89eWrVtVWRe6xta88eeP0pcR4NGDWw/h4UXho+TJ07a3qdlm5ZF1m/Zuuh0admxfQcvv/Ay0WHRttejYx7leNJxsrOyy+Q9S6LC3ABsmiYTJkzgpptuolGjRrb5PXr04K677qJGjRrEx8fz1FNP0aVLF37++Wfc3d3P209ubm6Ra4ZpaWnlUr+IiFyGh0fhGRI7vXdxxNSMwTCM8+5NiY6JPrub8/dzoc6mv4chE9M2z8nJ6dxMm/z8/AvWU+Rm4rO7tFqtl/oIJWY4GeddhivIL3ovjNVq5bHJj9Hr9l5F55tW3D3O/y4ubxUmzIwZM4YdO3awfn3Ra4oDBgyw/dyoUSNatmxJjRo1WLJkCf369TtvP9OnT2fq1KllXq+IiJSMYRjFutRjT4FBgXTs0pH3336f+0fef0Wty3Xq1eGnDT8xYNCf319bNm2hbr26AAQFBwFwPOk4jZs2BmBn3M4rep+ff/q5yPts3bK1xPsJDg4mIz2DzMxM2+fduaNoPY2bNWb/3v3UrFWzyHyLacFitZT4PUtbhQgzY8eO5euvv2bt2rVUq1btkutGRERQo0YN9u7de8HlkyZNYsKECbbptLQ0oqKiSrVeERGpvF6c9SK9b+5Nt/bdeOzJx2jQqAFOTk78su0X9v6xlybNm1xy+9HjRvPA0Ado3KwxHTp14Ltl37Hk6yV8+c2XAHh6etKidQtmvzybqOpRnEo5xQvPvlDiOh98+EHGPjSWZtc3o03bNnz5+Zfs+W0PNaJrlGg/17e8Hk8vT6ZNmcb9I+/nl59/Yf4n84us8+gTjzL4zsFEVo3k9jtux8nJid07d7Nr1y4e+/djJa69tNn1nhnTNBkzZgwLFixg1apVxMTEXHablJQUDh8+TERExAWXu7u74+fnV+QlIiJSXDE1Y1i1YRUdOnfguSnP0bltZ27pcAvvv/U+o8aN4omnnrjk9j1v68lzM57j9Vdfp32r9nw490NefetVbuxwo22dV994lfyCfLp16MbkxyfzxNOX3ueF9L2zL48+8Sj/eeo/3Nz+Zo4cPsLwfwwv8X6qBFbhjffe4PsV39OpTScW/m8hjz1ZNKB0ubkLH3/5MbGrY7m146306NKDt+a8RbWoS5+AKC+GWdx+tTIwatQoPv30UxYvXky9evVs8/39/fH09CQjI4MpU6bQv39/IiIiSEhI4Mknn+TQoUP89ttv+Pr6XvY90tLS8Pf3JzU1tfSDzcmTMGsWBAWBj0/p7ruMLPyjKscyPLk1JpHaVTIvv4FUfBkZkJICjzwCwcH2rqbSOpl1klkbZxHkFYSPm2P8vtuTO+7Ud6lPtahquLpX7IfIyZU7d5kp1DsUV+eS/3/OyckhPj6emJiY8+5HKsn3t10vM7355psAdOrUqcj8efPmMXz4cJydnYmLi+PDDz/kzJkzRERE0LlzZz7//PNiBRkRERGp/OwaZi53UsjT05PvvrPTne8iIiLiECrUc2ZERERESkphRkRERByawoyIiIg4NIUZERERcWgKMyIiIuLQFGZERETEoSnMiIiIiENTmBERESlnM56fQee2ne1dRqWhMCMiIvIXYx8aS6hPKKE+oUT4R9D8uuY8Nu4xzpw+Y+/SbOZ/PJ/aVWvbu4wKo0KMmi0iIlKRdLmlC6++9SqWAgt7ft/D+IfHk5aaxtv/fdvepckF6MyMiIjI37i5uxEWFkZk1Ug6d+1Mn/59WLNqjW251Wrlpekv0bRuU6oFVqNz286sWrmqyD6OHT3Gg8MepG5UXaJDo7ml/S38vOXnC77fwYSDtG7SmsfGPYbVaiUvL4+p/55KkzpNiA6Npnun7vy49kcAflz7I/8c+U/SUtNsZ5BmPD/jgvs9dznr0w8/pfl1zYkOi+axcY9hsVh4bdZrNKzZkAbRDZg1Y1aR7d587U06tu5IdGg0zeo14/Hxj5ORkWFbfvjQYQbfNZjroq6jTkQdmjZpytKlSwE4ffo09957LyEhIXh6elKnTh3mzZtX4v8HJaEzMyIiUi5M08SSd+kx+cqKs5uBYRhXtG1CfAKrV67GxeXPr8x33niHN197k5defYnGTRvz6UefMuTuIazbso6atWuSkZFBn+59iIiI4MPPPyQ0LJS47XFYrdbz9v/brt+4u8/dDLh3AP+e+m8A/jnynxw+dJh3/vsOYRFhLP1mKQPvGEjs5lha3dCK5158jheff5ENv2wAwNvb+5L1/7DiB+YvnE9CfAL3D76fQwcPUat2LRYvX8yWzVsY9/A42ndqT8vWLQFwcnLi+f97nqgaURxKOMS/HvkXz/77WWa8UhianpjwBHl5eSxcvhB3D3eSE5Lx8SkcTf6pp55i9+7dLFu2jODgYPbt20d2dvYVHfviUpgREZFyYckzWfDkLru8d79pDXFxL36YWblsJdFh0VgtVnJycgB49oVnbcvfePUNxj4yljvuugOAp//zND+u/ZG3X3+bF2e9yIIvFpByMoUVsSuoElgFgJq1ap73Pls2b2HwXYP556P/ZPS40QDEH4hn4f8W8usfvxIeEQ7A6HGjWbVyFZ99/BmTp0zGz98PwzAICwu77GcxrSavvvkqPr4+1Ktfjxs73Mj+vfv5bMFnODk5UbtubV6b+Rob1m2whZmHRj9k275GdA2eeOoJHh//uC3MHDl8hN59elO/YX0sVgutGrXC1dkVgEOHDtG8eXNatizcV3R0dLGP+5VSmBEREfmbGzvcyIxXZpCdnc0n//2E/fv284+R/wAgPS2dpMQkWt/Qusg2rW9oza64wrC2c8dOGjdpbAsyF3L0yFHuuu0unnj6CUaOGWmbH7c9DtM0uaHZDUXWz8vNIzAwsMSfJap6FD6+PrbpkNAQnJ2dcXJyKjLvxIkTtun1set55aVX+OP3P0hPT8dSYCEnJ4fMzEy8vb154OEHeHz846z5YQ03drqRwQMG06J5CwAefvhh+vfvz7Zt2+jWrRt9+/alXbt2Ja67JBRmRESkXDi7GfSb1tBu710SXt5etjMp016axh097uClaS/xxNNP2Nb5+2Ur0zRt8zw8PS77HkHBQYSHh7Poy0XcO/RefP18gcL7cZydnfl+3fc4ORe9tdXb5+KXky7GxbXoV71hGBecZ1oLLwEePnSYQf0HMez+YTzx1BNUqVKFzRs3M37UeAryCwAYPHwwnW/uzHfLv2PN92to26YtL7/8MmPHjqVHjx4cPHiQJUuW8P3339O1a1dGjx7NSy+9VOLai0s3AIuISLkwDAMXdye7vK70fplzJk6ayBuz3yApMQlfP1/CI8LZvHFzkXW2bN5CnXp1AGjQqAE743Zy+tTpi+7Tw8ODj7/8GHd3d+7uezcZ6YU32DZu2hiLxcLJEyepWatmkde5y0qubq5YLJar+kwXs33bdgoKCpg6fSotW7ekVp1aJCUmnbde1WpVGXb/MN775D3GPzKed99917YsJCSE4cOH8/HHH/PKK6/wzjvvlEmt5yjMiIiIXMaNHW6kXv16vPJ/rwAwevxoXpv1Gou+XMS+P/bxn6f/w84dO3lw1IMA9LurH6FhoQwbOIzNGzeTEJ/AN4u+YcvmLUX26+3tzSdffYKLswsD+w0kIyODWnVq0X9Af8Y8OIZvF3/LwYSD/PLzL8yeOZvvv/seKLx0lJmRydrVa0k5mUJWVlapfdbomtEUFBTw3pvvkRCfwBeffcEH739QZJ1/P/5vVn2/ikMJh4jbHsea1WuoX78+AE8//TSLFy9m37597Nq1i2+//da2rKwozIiIiBTDyLEj+fi/H3P0yFEeePgBHh77MM88+Qwd23Rk1cpVfPTFR9SsXXhpys3NjS8Wf0FwSDCD+g+iU5tOvDbzNZydnc/br4+PD58t/AzTNBnUbxCZmZnMfms2d91zF1OenEK75u0YcvcQtm3ZRmTVSKDw/pxh9w/jwWEPUj+6PnNmzSm1z9m4SWOefeFZXpv1Gh1bd+Srz79i8tTJRdaxWCw8MeEJOrTswOD+g6lbry5vvPGG7bNPmjSJJk2a0KFDB5ydnZk/f36p1Xchhmma9umTKydpaWn4+/uTmpqKn59f6e785EmYNQuCgsDH5/LrVwAL/6jKsQxPbo1JpHaVTHuXI6UhIwNSUuCRRyA42N7VVFons04ya+MsgryC8HFzjN93e3LHnfou9akWVQ1Xd1d7lyNlxGJasFgthHqH2rqZSiInJ4f4+HhiYmLw8Ch6n1FJvr91ZkZEREQcmsKMiIiIODSFGREREXFoCjMiIiLi0BRmRESkbJiFD5ITuZjS+vOhMCMiIqUun3ysppW83Dx7lyIV2Lnn47i6Xl3Hm4YzEBGRUmfFSrI1GZeThV8zbu5uV/0UXql4LKYFq2klxzkHi3Pxn0hsmiZZWVkkJycTEBBwwefvlITCjIiIlIkkkiAfCo4X4GQ4gbJMpWOaJlbTSrpbOs5OJQ8kAQEBhIeHX3UdCjMiIlJmkkgi2ZqMK3pwXmWUlZfFmZwzDG82nCqeFx8h/EJcXV2v+ozMOQozIiJSpqxYySXX3mVIGcgmmwxrBq7uruc9wbc86QZgERERcWgKMyIiIuLQFGZERETEoSnMiIiIiENTmBERERGHpjAjIiIiDk1hRkRERByawoyIiIg4NIUZERERcWgKMyIiIuLQFGZERETEoSnMiIiIiENTmBERERGHpjAjIiIiDk1hRkRERByaXcPM9OnTadWqFb6+voSGhtK3b1/27NlTZB3TNJkyZQqRkZF4enrSqVMndu3aZaeKRUREpKKxa5iJjY1l9OjRbNq0iZUrV1JQUEC3bt3IzMy0rTNjxgxmzpzJnDlz2LJlC+Hh4dxyyy2kp6fbsXIRERGpKFyuZKP8/HySkpLIysoiJCSEwMDAK3rz5cuXF5meN28eoaGh/Pzzz3To0AHTNHnllVeYPHky/fr1A+CDDz4gLCyMTz/9lIceeuiK3ldEREQqj2KfmcnIyODtt9+mU6dO+Pv7Ex0dTYMGDQgJCaFGjRo88MADbNmy5aqKSU1NBbCFo/j4eJKSkujWrZttHXd3dzp27MiGDRsuuI/c3FzS0tKKvERERKTyKlaYmTVrFtHR0bz77rt06dKFBQsWsH37dvbs2cPGjRt55plnKCgo4JZbbqF79+7s3bu3xIWYpsmECRO46aabaNSoEQBJSUkAhIWFFVk3LCzMtuzvpk+fjr+/v+0VFRVV4lpERETEcRTrMtOGDRtYvXo1jRs3vuDy1q1bM2LECN566y3ef/99YmNjqVOnTokKGTNmDDt27GD9+vXnLTMMo8i0aZrnzTtn0qRJTJgwwTadlpamQCMiIlKJFSvM/O9//yvWztzd3Rk1alSJixg7dixff/01a9eupVq1arb54eHhQOEZmoiICNv85OTk887W/LUGd3f3EtcgIiIijsmu3UymaTJmzBgWLFjAqlWriImJKbI8JiaG8PBwVq5caZuXl5dHbGws7dq1K+9yRUREpAIqcTfTHXfcccFLPIZh4OHhQe3atRk0aBD16tW77L5Gjx7Np59+yuLFi/H19bXdB+Pv74+npyeGYTB+/HimTZtGnTp1qFOnDtOmTcPLy4tBgwaVtHQRERGphEp8Zsbf359Vq1axbds2W6j55ZdfWLVqFQUFBXz++ec0bdqUH3/88bL7evPNN0lNTaVTp05ERETYXp9//rltnccff5zx48czatQoWrZsydGjR1mxYgW+vr4lLV3+wmJe+J4jERERR1PiMzPh4eEMGjSIOXPm4ORUmIWsVivjxo3D19eX+fPnM3LkSP71r39d8GbevzJN87LvZxgGU6ZMYcqUKSUtVS7A3z2PYxmebD4WRHXfLDxdrfYuSURE5KqU+MzM+++/z/jx421BBsDJyYmxY8fyzjvvYBgGY8aMYefOnaVaqJSOdlVT8HPPIz3PlWUHIrAoy4iIiIMrcZgpKCjg999/P2/+77//jsViAcDDw+OirdNiXx4uVnrVTMTNyUJipidrDodSjBNkIiIiFVaJLzMNGTKE+++/nyeffJJWrVphGAY//fQT06ZNY+jQoUDhmEsNGzYs9WKldAR65nNrzeN8uy+C31P8CPTIo3nYGXuXJSIickVKHGZmzZpFWFgYM2bM4Pjx40DhE3kfeeQR/vWvfwHQrVs3unfvXrqVSqmq7pfFjdVOsv5ICBuOBlHFI49o/yx7lyUiIlJiJQ4zzs7OTJ48mcmTJ9vGPfLz8yuyTvXq1UunOilTTUJSOZ3jxq6T/qyID6d/vSMEeebZuywREZESKfE9Mx9//LHtZz8/vyJB5rHHHiudqqRcGAa0jzpBVZ8s8q1OLNkfQVa+s73LEhERKZESh5kxY8bw7bffnjf/kUceKRJ0xDE4G9C9ZhL+Zzuclh8IV4eTiIg4lBKHmfnz5zN48GDWrl1rmzd27Fi++OILVq9eXarFSfnwcLHSq1Yibs5nO5wOqcNJREQcR4nDTPfu3Xnrrbfo27cvW7duZdSoUSxYsIDVq1dz3XXXlUWNUg6qeORza0wSBia/n/Ljl+QAe5ckIiJSLCW+ARhg4MCBnD59mptuuomQkBBiY2OpXbt2adcm5ay6XzY3VTvJuiMhbDwaRBX3PGIC1OEkIiIVW7HCzIQJEy44PzQ0lObNm/PGG2/Y5s2cObN0KhO7aBySyqmzHU4rE8LpV/cIwV7qcBIRkYqrWGHml19+ueD8WrVqkZaWZluup/46vnMdTmdyXTma7sWSAxHcVe8IXq4We5cmIiJyQcUKM7qx99ribED3mCS+3FON1Fw3lh0Ip2+doziX+A4rERGRsqevJ7mgcx1O7s4WkjI9Wa0OJxERqaCKFWZGjhzJ4cOHi7XDzz//nE8++eSqipKK4a8dTntO+fHL8QB7lyQiInKeYl1mCgkJoVGjRrRr147bb7+dli1bEhkZiYeHB6dPn2b37t2sX7+e+fPnU7VqVd55552yrlvKSZRfNu2jTrL2cAgbjwUR4JFPzYBMe5clIiJiU6ww85///IexY8fy/vvv89Zbb7Fz584iy319fbn55pt577336NatW5kUKvbTOCSVU9mu7DwZwMqEMPqrw0lERCqQYj9nJjQ0lEmTJjFp0iTOnDnDwYMHyc7OJjg4mFq1aqmTqZK7KeokZ3LdOKIOJxERqWCu6KF5AQEBBAQElHIpUpE5G3BrTBJf7anGmbMdTn3qHMPFSXcFi4iIfambSYrt/A6nEHU4iYiI3SnMSIkE/KXD6Y9TfmxTh5OIiNiZwoyUWJRfNh2iTgCw6VgwB85427kiERG5linMyBVpFJJG45AzAKxMCONElpt9CxIRkWtWicNMdnY2WVl/jqR88OBBXnnlFVasWFGqhUnFd1O1k1TzzaLA6sTS/RFk5jvbuyQREbkGlTjM9OnThw8//BCAM2fO0KZNG15++WX69OnDm2++WeoFSsXldLbDKcA9j4x8V5btj6DAqhZ9EREpXyUOM9u2baN9+/YAfPnll4SFhXHw4EE+/PBDZs+eXeoFSsX21w6n41kerD6oMZxERKR8lTjMZGVl4evrC8CKFSvo168fTk5O3HDDDRw8eLDUC5SKL8Ajn+41z3Y4nfbl5+NV7F2SiIhcQ0ocZmrXrs2iRYs4fPgw3333nW34guTkZPz8/Eq9QHEM1Xz/7HDafCyI/afV4SQiIuWjxGHm6aefZuLEiURHR9OmTRvatm0LFJ6lad68eakXKI7jrx1O3x9Uh5OIiJSPEoeZO++8k0OHDrF161aWL19um9+1a1dmzZpVqsWJ47mp2kmi1OEkIiLl6IqeMxMeHk7z5s1xcnIiLS2NRYsW4evry3XXXVfa9YmDUYeTiIiUtxKHmbvvvps5c+YAhc+cadmyJXfffTdNmjThq6++KvUCxfG4/63DaZU6nEREpAyVOMysXbvW1pq9cOFCTNPkzJkzzJ49m+eee67UCxTHdK7DyQmTvad9+TlJHU4iIlI2ShxmUlNTCQwMBGD58uX0798fLy8vevXqxd69e0u9QHFc1XyzaX+uwylRHU4iIlI2ShxmoqKi2LhxI5mZmSxfvtzWmn369Gk8PDxKvUBxbI1C0mhyrsMpIYwTWe72LUhERCqdEoeZ8ePHc++991KtWjUiIyPp1KkTUHj5qXHjxqVdn1QCN1Y7SXW/TApMJ5aow0lEREpZicPMqFGj2LRpE3PnzmX9+vU4ORXuombNmrpnRi7IyYBuMcep4pFHZr4LS9XhJCIipcjlSjZq0aIFLVq0KDKvV69epVKQVE7uzlZ61TrG/36PIvlsh9Mt0ccxlGlEROQqXVGYOXLkCF9//TWHDh0iLy+vyLKZM2eWSmFS+fi7F9CjZhJf741k72lfqnjk0SritL3LEhERB1fiMPPDDz9w++23ExMTw549e2jUqBEJCQmYpsn1119fFjVKJVLVN5uO1U+w+lAoPyUGUcUjj9pVMu1dloiIOLAS3zMzadIkHn30UXbu3ImHhwdfffUVhw8fpmPHjtx1111lUaNUMg2C02gaegaAHxLCSFaHk4iIXIUSh5nffvuNYcOGAeDi4kJ2djY+Pj48++yzvPjii6VeoFRO7ar+2eG0dH8EmXnqcBIRkStT4jDj7e1Nbm4uAJGRkezfv9+27OTJk6VXmVRq53U4HYggXx1OIiJyBUocZm644QZ+/PFHoLCD6dFHH+X5559nxIgR3HDDDSXa19q1a7ntttuIjIzEMAwWLVpUZPnw4cMxDKPIq6TvIRXXuQ4nD2eLrcNJYziJiEhJlTjMzJw5kzZt2gAwZcoUbrnlFj7//HNq1KjB+++/X6J9ZWZm0rRpU9vAlRfSvXt3EhMTba+lS5eWtGSpwAo7nBJxMkz2nfZlq8ZwEhGREipxN1PNmjVtP3t5efHGG29c8Zv36NGDHj16XHIdd3d3wsPDr/g9pOKL9M2hY1Qyqw+F8VNiEAEe+dSpkmHvskRExEGU+MzMX2VkZJCWllbkVdrWrFlDaGgodevW5YEHHiA5OfmS6+fm5pZ5TVL6GgSn0zS08JkzPySEkpypDicRESmeEoeZ+Ph4evXqhbe3N/7+/lSpUoUqVaoQEBBAlSqle4mgR48efPLJJ6xatYqXX36ZLVu20KVLF9sNyBcyffp0/P39ba+oqKhSrUnKTruqKdTwy8RiOrH0QAQZ6nASEZFiKPFlpnvvvReAuXPnEhYWhlGGz6MfMGCA7edGjRrRsmVLatSowZIlS+jXr98Ft5k0aRITJkywTaelpSnQOIjCDqckvtxTjdM57iw9EMEddY/i6qS7gkVE5OJKHGZ27NjBzz//TL169cqinkuKiIigRo0a7N2796LruLu74+6uSxSOys3ZpFetRL7cE8WJLA9+SAjl1hiN4SQiIhdX4stMrVq14vDhw2VRy2WlpKRw+PBhIiIi7PL+Uj7+2uG0/4wvWxID7V2SiIhUYCU+M/Pee+8xcuRIjh49SqNGjXB1dS2yvEmTJsXeV0ZGBvv27bNNx8fHs337dgIDAwkMDGTKlCn079+fiIgIEhISePLJJwkODuaOO+4oadniYCJ9/uxw2pIUSKBnrsZwEhGRCypxmDlx4gT79+/nvvvus80zDAPTNDEMA4vFUux9bd26lc6dO9umz93rMmzYMN58803i4uL48MMPOXPmDBEREXTu3JnPP/8cX1/fkpYtDqhBcDqnctz5NTmAHckBCjMiInJBJQ4zI0aMoHnz5nz22WdXfQNwp06dMC/xyNfvvvvuivctlUOUbxa/JgdQYOqmGRERubASh5mDBw/y9ddfU7t27bKoR0RERKRESnwDcJcuXfj111/LohYRERGREivxmZnbbruNRx55hLi4OBo3bnzeDcC33357qRUnIiIicjklDjMjR44E4Nlnnz1vWUlvABYRERG5WiUOM1artSzqEBEREbkiVzXQpIiIiIi9KcyIiIiIQ1OYEREREYemMCMiIiIOTWFGREREHFqxupnS0tKKvUM/P78rLkZERESkpIoVZgICAi47BtOVDDQpIiIicrWKFWZWr15d1nWIiIiIXJFihZmOHTuWdR0iIiIiV+SKbgBet24dgwcPpl27dhw9ehSAjz76iPXr15dqcSIiIiKXU+Iw89VXX3Hrrbfi6enJtm3byM3NBSA9PZ1p06aVeoEiIiIil1LiMPPcc8/x1ltv8e677xYZMbtdu3Zs27atVIsTERERuZwSh5k9e/bQoUOH8+b7+flx5syZ0qhJREREpNhKHGYiIiLYt2/fefPXr19PzZo1S6UoERERkeIqcZh56KGHGDduHJs3b8YwDI4dO8Ynn3zCxIkTGTVqVFnUKCIiInJRxWrN/qvHH3+c1NRUOnfuTE5ODh06dMDd3Z2JEycyZsyYsqhRRERE5KJKHGYAnn/+eSZPnszu3buxWq00aNAAHx+f0q5NRERE5LKuKMwAeHl50bJly9KsRURERKTEShxmOnfufMlxmlatWnVVBYmIiIiURInDTLNmzYpM5+fns337dnbu3MmwYcNKqy4RERGRYilxmJk1a9YF50+ZMoWMjIyrLkhERESkJK5obKYLGTx4MHPnzi2t3YmIiIgUS6mFmY0bN+Lh4VFauxMREREplhJfZurXr1+RadM0SUxMZOvWrTz11FOlVpiIiIhIcZQ4zPj7+xeZdnJyol69ejz77LN069at1AoTERERKY5ih5kDBw4QExPDvHnzyrIeERERkRIp9j0zderU4cSJE7bpAQMGcPz48TIpSkRERKS4ih1mTNMsMr106VIyMzNLvSARERGRkii1biYREREReyh2mDEM47xhDC41rIGIiIhIeSj2DcCmaTJ8+HDc3d0ByMnJYeTIkXh7exdZb8GCBaVboYiIiMglFDvM/H3cpcGDB5d6MSIXU2B1IqfACQ8Xq71LERGRCqbYYUYt2WIPTkbhjeenc9x4f0dN/NzyCfHKJcQrh1CvXEK8chVwRESucSV+aJ5IeYrwyaFBUCpH0r1Iy3O1vfaf8bGtcy7gNAk5Q6Rvjh2rFRERe1CYkQrNxcmkc43C5xvlFDhxIsvd9krO8igScBJSvbij7lHCvHPtXLWIiJQnhRlxGB4uVqL8sonyy7bNOxdwfjlehcPpXiw9EMFd9Q7j42axY6UiIlKe9JwZcWjnAk73mokEeuSSle/C0v0R5Fv12AARkWuFwoxUCm7OJr1qJeLhYuFEtgc/JITxt4dWi4hIJWXXMLN27Vpuu+02IiMjMQyDRYsWFVlumiZTpkwhMjIST09POnXqxK5du+xTrFR4fu4F9KiZiJNhsv+MDz8lBtq7JBERKQd2DTOZmZk0bdqUOXPmXHD5jBkzmDlzJnPmzGHLli2Eh4dzyy23kJ6eXs6ViqOI9MmhU/VkALYmBfLHKZ/LbCEiIo7OrjcA9+jRgx49elxwmWmavPLKK0yePJl+/foB8MEHHxAWFsann37KQw89VJ6ligOpH5TOqWw3tidXYdXBUPzd89XhJCJSiVXYe2bi4+NJSkqiW7dutnnu7u507NiRDRs2XHS73Nxc0tLSirzk2tO2agrRfplYTCeW7o8gI8/Z3iWJiEgZqbBhJikpCYCwsLAi88PCwmzLLmT69On4+/vbXlFRUWVap1RMTgbcEpNU2OFU4MKS/RHkW9ThJCJSGVXYMHPO30fmNk3zkqN1T5o0idTUVNvr8OHDZV2iVFDnOpw8XQo4me3BDwfV4SQiUhlV2DATHh4OcN5ZmOTk5PPO1vyVu7s7fn5+RV5y7fJzL6B7zSR1OImIVGIVNszExMQQHh7OypUrbfPy8vKIjY2lXbt2dqxMHE2kTw6d1eEkIlJp2bWbKSMjg3379tmm4+Pj2b59O4GBgVSvXp3x48czbdo06tSpQ506dZg2bRpeXl4MGjTIjlWLI7ouKJ1TOW78cryww8nPPZ9wdTiJiFQKdg0zW7dupXPnzrbpCRMmADBs2DD++9//8vjjj5Odnc2oUaM4ffo0bdq0YcWKFfj6+tqrZHFgN0SmcDrHlYRUH5buj+Cu647g61Zg77JEROQq2TXMdOrUCfMSd2QahsGUKVOYMmVK+RUllZaTAbdEH2fBHldSctxZuj+CfnWP4Oqsu4JFRBxZhb1nRqQsuDmb9LR1OLnzvcZwEhFxeAozcs0pHMOpsMPpQKoPm9XhJCLi0BRm5JoU8ZcOp5+TAtmjDicREYelMCPXrOuC0rk+7DQAqw+GkpTpbueKRETkSijMyDXthsgUYvwzbGM4pefZ9Z54ERG5Agozck0zDLg5+jhBnrlknx3DKU9jOImIOBSFGbnmuTmb9KxZ2OGUog4nERGHozAjQmGHU8+zHU7xqT5sOqYOJxERR6EwI3JWuE8OXc52OG07HsieFD1pWkTEESjMiPxFvaB0rg87BcCqQ6EkZXjYuSIREbkchRmRv7kh8hQx/hlYTYOlB8JJy1WHk4hIRaYwI/I3f+9wWnpAHU4iIhWZwozIBbg5m/Sq9WeH00p1OImIVFgKMyIX4etWQM9aiTgbVhJSfdh0LMjeJYmIyAUozIhcQrh3Lp1rnOtwqsLv6nASEalwFGZELqNeYAYtznY4rT4USqI6nEREKhSFGZFiaBN5ippnO5yWqcNJRKRCUZgRKYZzHU7B5zqcNIaTiEiFoTAjUkyuziY9ayXi5VJASo47K+PDsarDSUTE7hRmRErA162AHuc6nNK82XRUHU4iIvamMCNSQuHeuXQ52+H0S3IVflOHk4iIXSnMiFyBuoEZtAwv7HBacyiUY+pwEhGxG4UZkSvUOuIUtQLOdThFqMNJRMROFGZErpBhQNcaxwn2zCGnwJkl6nASEbEL/VPyCq1IWMHXv30JkQfAzQ2cne1dUqUWiR8NrKE0MEOoaQbiUkFyuOvZMZz+93sUp3LcWREfTs9aiTgp04iIlBuFmSt0KP0QsckbwcfelVxDzuZFD9OFumYwDcwQGlrD7B5wfNws9KyVyMI/qnIwzZuNR4O4sVqKXWoREbkWKcxcofZV2xNU4AYrVoKPN7jrBtCyYsFKvHGaXU7J/G6cINPIZ4eRxA6SwDkOKAw4QyzNGGu5AYPyPy0S5p1L1xrJrEgIZ3tyFQI986gflF7udYiIXIsUZq5QvcB61LMGQdoecA0CV52iKXMWsGJy0DjDbiOZ3UZykYDzrstWgkwv7rU2tUt5dQIzOJVziq1Jgaw5FIq/ez6RPjl2qUVE5FqiMCMOxQmDGLMKMWYVelHPFnDmOW/jFZcNzHBZR3R+ADeaNexSX+uIU5zOcWX/GV+WHYjgrnqH8XMvsEstIiLXiopxF6XIVXDCYITlevpa6mM1TB5zXc4B45RdajEM6BqdTIiXOpxERMqLwoxUCgYGTxV05nprBOlGHmNcv+UM2XapxdXJpGfNRLxcC2wdThrDSUSk7CjMSKXhhjOz8ntR1fTjsJHKBNdl5GOxSy0+bhZ61iwcw+lgmjcbNIaTiEiZUZiRSiUQT17L742X6coWp6NMc4nFxD6nRcK8c+kaXTiG06/JVdh9UmM4iYiUBYUZqXTqmEH8X8GtGCZ86byLT5132K+WKhm0iih85kzs4VCOpauFX0SktCnMSKXUwRrDo5YbAZjhvI71xkG71dIq/DS1A9JtYzilagwnEZFSpb9Vr9RP70LsS+CSDulOkFFBcqFLMLhWBbeowv86B4FRQWorZ0MtzdlvnGah824ec13OJ/l3UdMMLPc6DAO6RCeT9ocryVkeLNkfwZ31juDmrLuCRURKg8LMlcpNh8wkMADz7KsiyEuFvP2QeXba8CgMNa5VwbUauFUD5ypgVP7/9YUdTp04aJxhm9MxRrt+w2d5dxOAZ7nXcq7D6X97ojid48538eH00hhOIiKlovJ/o5WV5oMhuAV8/AkE+IOXl70rAtMCBUmQdwTyj0L+MTBzCsNN3v4/1zM8IXAIuNexX63lxBVnZuX3ZJDbFxwx0njEdRnv5PfBlfIfGNT73BhOe6pyKM2bDUeDuanayXKvQ0SkslGYuVI+oRDiBIQUXsqpKMMZuFUHr9aFP5sWKDgO+Uf+EnAOgZkNpz6CkLHgEmLfestBIJ7Mye/NYNf/sdXpKM+7xPJMQWe7jOEU6pVL1+jjfBcfwa/JAQR65NEgOK3c6xARqUyuzZsprhWGM7hGFoabgH6F4SViGrhWLww0KfPAmmXvKstFbTOIGQXdcTINvnLexcfOv9qvliqZtD7X4XQohKPp5X/ZS0SkMlGYudYYrhA4HJwDwHKi8AyNaZ8Hy5W3DtZoJpztcHrJeT3rnBLsVkvL8NPUqZKOFYNlB8LV4SQichUUZq5Fzr4QeB8YbpC3D1IX27uicjPU0ox+lgaFYzi5LGefkWKXOgwDutRIJtQrh1yLM0v2R5Jr0a+jiMiV0N+e1yrXSKgyCDAgayNk/mjvisqFgcG/CzrRwhpJppHPGNdvOW2nMZxcnEx61krE27WA0zlurIgP0xhOIiJXoEKHmSlTpmAYRpFXeHi4vcuqPDwagm+Pwp9Tv4acP+xbTzk51+FU1fTjqJHGI65L7TaGk7erhV61EnExrLYOJxERKZkKHWYAGjZsSGJiou0VFxdn75IqF59O4NkCsMLpj6Ag2d4VlYsqePJ6/m34mG787HSM/7issdsYTiFeudwcfRyAX5MD2HXSzy51iIg4qgofZlxcXAgPD7e9QkIqTivxhkOpWBz9soBhQMCd4BZd+Eyaa6jDqZYZyIyCW3EyDRY67+ZD5+32q6VKJm3OdjitVYeTiEiJVPgws3fvXiIjI4mJiWHgwIEcOHDgkuvn5uaSlpZW5FUW5qzay6AvfucFs0aZ7L9cGS5QZVjhk4EtJ6+pDqf21mgmWm4CYKbzj6x1irdbLS3+1uF0JsfVbrWIiDiSCh1m2rRpw4cffsh3333Hu+++S1JSEu3atSMl5eIdKNOnT8ff39/2ioqKKpPaagR5A/CuWZUvsivBZQFnn8KWbcP9bIfTIjAd/bRT8Qy2NKX/2Q6nx12+qzAdTkv3R6jDSUSkGCr035Q9evSgf//+NG7cmJtvvpklS5YA8MEHH1x0m0mTJpGammp7HT58uExqu61pJOPaVgVgcno4m7PdyuR9ylWRDqdNkLXB3hWVCwODyQWdaGmtautwOmXnDicf13xO56rDSUSkOCp0mPk7b29vGjduzN69ey+6jru7O35+fkVeZWVcu6r0Mk6Sj8HI44Ecyi//8X5KnUcD8OtZ+HPq15Czx771lJPCDqceVLN1OC0hz44dTj1rJeLiVNjh9OMRdTiJiFyKQ4WZ3NxcfvvtNyIiIuxdCgBOhsFLxj4au+Rw2urM/UlBpFsrwTDI3h3BsyWFHU4fQ/610eEUgCdzznY4bXNK5D8uq+3Y4ZTHLWc7nHacCGC3OpxERC6qQoeZiRMnEhsbS3x8PJs3b+bOO+8kLS2NYcOG2bs0G0/Dyrv+RwlztrA335V/Hg+sJB1O/cEtprDD6dS11eH0f/mFYzgtcv6ND51/sVstNQP+7HCKPRzCsXQPu9UiIlKRVegwc+TIEe655x7q1atHv379cHNzY9OmTdSoUbE6iMKdC3g3PAUPw8rqbA+mn6oE/4o2XKDK0L90OH14zXQ43WTW4LGzHU4vO/9IrJ07nGpXScdqGiyLjyBNYziJiJynQoeZ+fPnc+zYMfLy8jh69ChfffUVDRo0sHdZF9TEPZ+XQs4A8F6qL5+nedm3oNLg7AOBI852OO2H1IXXTIfTvZam3GlpiGnA4y7f8Ydx0i51nOtwCvHKIaegsMMpz1IJLmWKiJSiCh1mHE1vn2weqVL4XJt/nwxgU6XocAr/S4fT5mtqDKcnCzrSylqVLCOff7ouIQX7XGpzdTLpWTMRT5cCUnLc+T4h7FrJlCIixaIwU8r+GZBOb+8s8jF4uFJ1OPUq/Dnta8j53b71lBNXnJmZ34Mo0982hpO9Opx83Ao7nJwMk/hUHzYnBtqlDhGRikhhppQZBrwUcpqm7nm2Dqe0StHh1AE8WwEmnP4E8o/bu6JyUdjh1Btf041fnBJ51o4dTuHeuXSuXthZ9nNSIHtP+dilDhGRikZ3E5YBDyd4JyyFPkdD2ZvvytjjgcwNT8HZkTONYUBAv8KbgfPiCzucQsaCk7e9KytzNc1AXsrvwcOuX7PY+TdqmYHcZ7neLrVcF5ROSrYb25Or8MPBUH4+XgWsVigogDl/gMulh/uQK2exFhCS2RpnJxfyDf07UATA1bQSYi3gj4KTBN9qv2diKcyUkTAXK++Fp3DnsWBisz2Ydsqfp4JS7V3W1TnX4XRyNlhSCjucgh4onF/JtTOr83hBe15wXcss5x+JNgPobK1pl1raVk3hTI4bCWnepGS7/7kgKccu9VxLXPG1dwkiFYrT2VdOWr5d66j830J21Mg9n5khpxmVHMT7qT7Uds3nHj8Hf17LuQ6nk3Mg70DhGE7+/QvP3FRyg6xNOGA5xRfOO3nCZQUf5t9JPbP8/yXiZEDPWokkZnpQYDUgOwfS0qBfP/CvBI8FqKDSctNYsHsBvu5+eLrqmT8iANn5OaTnphHTeqBd61CYKWM9fXJ4ND+Nl0/78dTJAKJdC2jrmWfvsq6OazhUubfwUlPWZnAJA5/29q6qzBkYPFHQgYPGGTY7HWGs67d8lnc3QZR/G75hQKTP2TMxThmQcwrq+EJwULnXcq04mWWSe+oUPl4GTm5qJxMBMPMyyM06hW+I++VXLkO68FsOxgSkc7t3FgUYPHw8iIOVosOp/l86nL65pjqcXs7vQXWrP4lGOuPt2OEkIiKFFGbKgWHAjLMdTmesToyoTB1OXq251jqc/PHgtYLCDqftTolMdVlltw4nERFRmCk3Hk7wblgKEc4F7M93ZczxQAoc/fvPMMD/DnCreXYMp7lgybR3VeXiXIeTs2nwtfPvzHPeZu+SRESuWQoz5SjUxcq74afwNKyszfbg+RR/e5d09WxjOAWC5RSc/hDMAntXVS7amdX5V0EHAF5x3sAqJ7VFi4jYg8JMOWvkns+s0NMAzEvz4dNKMYaTNwTed3YMpwPX1BhO91ibMMDSGNOAJ1xWsMc4Ye+SRESuOQozdtDdO4eJVQqfOfP0yQA2VJoxnAZTOIbTT5C5zt4VlZt/FbSnjbUa2UY+Y12XcNJOYziJiFyrFGbsZHRABn18/uxwiq8UHU7XgV/vwp/TvoWc3+xbTzk51+EUbQ042+G0hFyujUttIiIVgcKMnRgGvBh8mmbueaRanbg/KYhUS2XocGr/tw6nJHtXVC7+7HBy51enJHU4iYiUI4UZOzo3hlOEcwEH8l0Zk1zZOpxyCx+sd410OEWbVXg5vzvOpsE3znt43/lne5ckInJNUJixs792OK3L9uC5Stnh9ME10+HU1qzOE2c7nGY7b+QHp/12rkhEpPJTmKkA/trh9N80Hz6uNB1OI8DwKBxlO3XBNdPhNNDahIFnO5wmuazkd3U4iYiUKYWZCqK7dw6Pne1weqbSdDiFFY7hhAFZWyBzrb0rKjf/KujADdaosx1O33KSa+NSm4iIPSjMVCCjAjLo65OFpdJ1ON1W+HPakmumw8kFJ17K7060NYAkI4NxrkvV4SQiUkYUZioQw4AXgk/TvNJ1ON0EXm24Fjuc5hTchp/pzg6nJJ5Rh5OISJlQmKlgPJzg7bAUIitdh1NfcKt1tsNpLlgy7F1VuahhBjDz7BhOS9ThJCJSJhRmKqC/dzj9p7J0OAUOAecgsJy+psZwamNG8WRBRwBeddnI90777FyRiEjl4mLvAhyVaZqY2dlQUAD5+YWvUlTfKZ9XqxxnbHIg80+5Utdw5R4/R39Mvhv4DoGUtyArAaxfQkC/wjM3ldydXEd8wQn+57yTKeYKqhZ4U88MLp2d5+cX/jnMzoYsR/8zUnGZ2dk45xbg7JyPk7V0f99FHJVzfj7OuQWYdu5WNUx7V1DG0tLS8Pf3JzU1FT8/v1LbrzUriz3Xtyi1/YmIiDiq4B9XEBIUVar7LMn3ty4ziYiIiEPTZaYrZHh6Um/lCpgzBwIDwcenzN4rxwr3JQWxPdeNGNcCPos4ib+zg59QM01IXQzZW8Bwg6CRhc+luQYcMs5wn8sC0o1culnq8B9LVwyu4lJbRgacOgVjxkBQUOkVKkWkZKfw2uY5BHkF4u1adr/vIo4kMz+DlKxTjPXwsGsdCjNXyDAMDE9PcHEBV9fCVxnxAuZUzaDP0RB+t7gz9rQz88JTcHX0W02C+0HKKcjbB+kfQvA/wbnyf0lEE8KLRi8ecl3MN277iC4I5kFLqyvfoatr4Z9DT0/wqgRPj66gDLKwuLtgcXfF6lZ2v+8ijsTi5IrF4oJh53sfdZnJQYS4WHkvPAUvw8r6StPh5Hy2wyn4mutwam1WY9LZDqfXXDaxUh1OIiJXTGHGgTRwL+CV0NMYmHyY5sNHqd72LunqOXlB4H1/juF05qtrZgynu62NuLegKQCTXVay20i2c0UiIo5JYcbBdPPO4fHANACmpPizPsvdzhWVAtdQqDIEcILsrZAZa++Kys1Ey03caK1OtlHAWNdvOaExnERESkxhxgGN9M+g39kxnEYlB7I/rxLc+uRRF/xvL/w5bSnk7LZvPeXEBSdm5HcnxlqFZCOTca5LyNEYTiIiJaIw44AMA6aHnKaFey5pVif+cTyIM5VhDCevduDVlsIxnD6F/GP2rqhc+OHOnPze+JsexDkd52mX7zWGk4hICSjMOCh3A94OP0VVlwLi810YnRxIvqN//xkG+PcBt9pnx3Cad82M4VSdAGbl98DFdGKZ817ecd5q75JERByGwowDC3a28l5YYYfTj9keTK10HU5n4NR/r5kOp1ZmNSaf7XCaow4nEZFiU5hxcPXdC3j1bIfTx2k+vHemknQ4Bd0HhifkH4QzX14zHU53WhsxuKAZAE+qw0lEpFgqwZ2jcot3Dk8EpjH9lD/PnQrAxYDh/g7eFeMSCoGDIeV9yP4ZCpLAuDYeVPYoJglBsN6jgLHOn/NZskGo9RL3RFks4FwAX20El2vjGNmDvzWfB9KP4erkgrPhbO9yRCoEi2kh31qAh2cYtPun3epQmKkkHvTPIM3qxOtnfJmSEgBUgkDjXrfwHprUhZB/1N7VlBsXYEaSweDIcA64uTIuMId5icl4XOrslBOQmFheJV6TXIFoexchUkFlpdm3YUNhppIwDJhYpfD5M+cCjQnc5+iBxrsduNWAglP2rqRc+QJzcrIZ5PIzO93hmcgWvJBV/8JjOOXkQHo63HYblOLI8FJUWm4a3+z5Bl93Xzxc7DsOjUhFkVOQQ3puOj0b3I49B1NRmKlEzgUaA5M5Z/yYevYMjcMHGteqha9rTBQwsyCSB10Xs9QtmTpOtfiHpeX5K1oywEyBWrdCcHC513mtyMs6ye6U3QR5BeHjVvnHEBMpjoy8DFKyUrg1uK5d69ANwJWMYcCjVdIZE1B4lmZqSgDzKsOwB9eoVmY1JhV0AGC280ZWOx2wc0UiIhWPwkwldKFAM1eBxmHdbW3MAEtjTAOecFnBH8ZJe5ckIlKhKMxUUucCzdizgeZZBRqH9q+C9rS2ViPLyOefrks4Tba9SxIRqTAcIsy88cYbxMTE4OHhQYsWLVi3bp29S3IIhgETFGgqBVeceTm/O9VMP44aaUxwXUY+FnuXJSJSIVT4G4A///xzxo8fzxtvvMGNN97I22+/TY8ePdi9ezfVq1e3d3kV3rlAA/DaGT+eTQkg3epEf58sO1cmJefDvy19eMTzc7Y6HeVJp7UMzGsNZnbhc2aOH4LcVHsXWWml5qRyIqOAvIJsvFwqwVhoIqUgqyCb1JwCjqWeJtjLfg0IhmlW7EertmnThuuvv54333zTNq9+/fr07duX6dOnX3b7tLQ0/P39SU1Nxa+021ZPnoRZsyAoCHwqdneDacLM0768dkatu47O2ec3PKt9iGFU6F9dEbmG1PO4jS8HTCvVfZbk+7tCn5nJy8vj559/5oknnigyv1u3bmzYsOGC2+Tm5pKbm2ubTktLK9MaHcW5MzQmMC/Vh4ILPa9EHEPmdViO345zyAow8u1djYiI3Z+KXaHDzMmTJ7FYLISFhRWZHxYWRlJS0gW3mT59OlOnTi2P8hyOYcDEwHQmBqbbuxS5atWh4B+FP2ZkQEoKPPKInjNThk5mnWTWxll6zozIX5x7zswjbcfZtQ6HuAHYMIqeRTBN87x550yaNInU1FTb6/Dhw+VRooiIiNhJhT4zExwcjLOz83lnYZKTk887W3OOu7s77u7u5VGeiIiIVAAV+syMm5sbLVq0YOXKlUXmr1y5knbt2tmpKhEREalIKvSZGYAJEyYwZMgQWrZsSdu2bXnnnXc4dOgQI0eOtHdpIiIiUgFU+DAzYMAAUlJSePbZZ0lMTKRRo0YsXbqUGjVq2Ls0ERERqQAqfJgBGDVqFKNGjbJ3GSIiIlIBVeh7ZkREREQuR2FGREREHJrCjIiIiDg0hRkRERFxaAozIiIi4tAUZkRERMShKcyIiIiIQ1OYEREREYfmEA/NuxqmaQKQlpZW+jtPT4fcXDh1CnJySn//IsWRlVX45zA9Hdzc7F1NpZWelU5uZi6n8k6R46rfdxGArPwscvNzSU9Lx62gdP/+Ofe9fe57/FIMszhrObAjR44QFRVl7zJERETkChw+fJhq1apdcp1KH2asVivHjh3D19cXwzBKdd9paWlERUVx+PBh/Pz8SnXf8icd5/Kh41w+dJzLh45z+SjL42yaJunp6URGRuLkdOm7Yir9ZSYnJ6fLJrqr5efnp1+WcqDjXD50nMuHjnP50HEuH2V1nP39/Yu1nm4AFhEREYemMCMiIiIOTWHmKri7u/PMM8/g7u5u71IqNR3n8qHjXD50nMuHjnP5qCjHudLfACwiIiKVm87MiIiIiENTmBERERGHpjAjIiIiDk1hRkRERByawswVeuONN4iJicHDw4MWLVqwbt06e5fk0NauXcttt91GZGQkhmGwaNGiIstN02TKlClERkbi6elJp06d2LVrl32KdWDTp0+nVatW+Pr6EhoaSt++fdmzZ0+RdXSsr96bb75JkyZNbA8Sa9u2LcuWLbMt1zEuG9OnT8cwDMaPH2+bp2NdOqZMmYJhGEVe4eHhtuX2Ps4KM1fg888/Z/z48UyePJlffvmF9u3b06NHDw4dOmTv0hxWZmYmTZs2Zc6cORdcPmPGDGbOnMmcOXPYsmUL4eHh3HLLLaSnp5dzpY4tNjaW0aNHs2nTJlauXElBQQHdunUjMzPTto6O9dWrVq0aL7zwAlu3bmXr1q106dKFPn362P5y1zEufVu2bOGdd96hSZMmRebrWJeehg0bkpiYaHvFxcXZltn9OJtSYq1btzZHjhxZZN51111nPvHEE3aqqHIBzIULF9qmrVarGR4ebr7wwgu2eTk5Oaa/v7/51ltv2aHCyiM5OdkEzNjYWNM0dazLUpUqVcz33ntPx7gMpKenm3Xq1DFXrlxpduzY0Rw3bpxpmvrzXJqeeeYZs2nTphdcVhGOs87MlFBeXh4///wz3bp1KzK/W7dubNiwwU5VVW7x8fEkJSUVOebu7u507NhRx/wqpaamAhAYGAjoWJcFi8XC/PnzyczMpG3btjrGZWD06NH06tWLm2++uch8HevStXfvXiIjI4mJiWHgwIEcOHAAqBjHudIPNFnaTp48icViISwsrMj8sLAwkpKS7FRV5XbuuF7omB88eNAeJVUKpmkyYcIEbrrpJho1agToWJemuLg42rZtS05ODj4+PixcuJAGDRrY/nLXMS4d8+fPZ9u2bWzZsuW8ZfrzXHratGnDhx9+SN26dTl+/DjPPfcc7dq1Y9euXRXiOCvMXCHDMIpMm6Z53jwpXTrmpWvMmDHs2LGD9evXn7dMx/rq1atXj+3bt3PmzBm++uorhg0bRmxsrG25jvHVO3z4MOPGjWPFihV4eHhcdD0d66vXo0cP28+NGzembdu21KpViw8++IAbbrgBsO9x1mWmEgoODsbZ2fm8szDJycnnpVIpHefumNcxLz1jx47l66+/ZvXq1VSrVs02X8e69Li5uVG7dm1atmzJ9OnTadq0Ka+++qqOcSn6+eefSU5OpkWLFri4uODi4kJsbCyzZ8/GxcXFdjx1rEuft7c3jRs3Zu/evRXiz7TCTAm5ubnRokULVq5cWWT+ypUradeunZ2qqtxiYmIIDw8vcszz8vKIjY3VMS8h0zQZM2YMCxYsYNWqVcTExBRZrmNddkzTJDc3V8e4FHXt2pW4uDi2b99ue7Vs2ZJ7772X7du3U7NmTR3rMpKbm8tvv/1GRERExfgzXS63GVcy8+fPN11dXc3333/f3L17tzl+/HjT29vbTEhIsHdpDis9Pd385ZdfzF9++cUEzJkzZ5q//PKLefDgQdM0TfOFF14w/f39zQULFphxcXHmPffcY0ZERJhpaWl2rtyxPPzww6a/v7+5Zs0aMzEx0fbKysqyraNjffUmTZpkrl271oyPjzd37NhhPvnkk6aTk5O5YsUK0zR1jMvSX7uZTFPHurQ8+uij5po1a8wDBw6YmzZtMnv37m36+vravvfsfZwVZq7Q66+/btaoUcN0c3Mzr7/+eltrq1yZ1atXm8B5r2HDhpmmWdj698wzz5jh4eGmu7u72aFDBzMuLs6+RTugCx1jwJw3b55tHR3rqzdixAjb3w8hISFm165dbUHGNHWMy9Lfw4yOdekYMGCAGRERYbq6upqRkZFmv379zF27dtmW2/s4G6ZpmuVzDkhERESk9OmeGREREXFoCjMiIiLi0BRmRERExKEpzIiIiIhDU5gRERERh6YwIyIiIg5NYUZEREQcmsKMiIiIODSFGREpd1OmTKFZs2Z2e/+nnnqKBx98sFjrTpw4kX/+859lXJGIXA09AVhESpVhGJdcPmzYMObMmUNubi5BQUHlVNWfjh8/Tp06ddixYwfR0dGXXT85OZlatWqxY8eO8wbmFJGKQWFGREpVUlKS7efPP/+cp59+mj179tjmeXp64u/vb4/SAJg2bRqxsbF89913xd6mf//+1K5dmxdffLEMKxORK6XLTCJSqsLDw20vf39/DMM4b97fLzMNHz6cvn37Mm3aNMLCwggICGDq1KkUFBTw2GOPERgYSLVq1Zg7d26R9zp69CgDBgygSpUqBAUF0adPHxISEi5Z3/z587n99tuLzPvyyy9p3Lgxnp6eBAUFcfPNN5OZmWlbfvvtt/PZZ59d9bERkbKhMCMiFcKqVas4duwYa9euZebMmUyZMoXevXtTpUoVNm/ezMiRIxk5ciSHDx8GICsri86dO+Pj48PatWtZv349Pj4+dO/enby8vAu+x+nTp9m5cyctW7a0zUtMTOSee+5hxIgR/Pbbb6xZs4Z+/frx15PWrVu35vDhwxw8eLBsD4KIXBGFGRGpEAIDA5k9ezb16tVjxIgR1KtXj6ysLJ588knq1KnDpEmTcHNz48cffwQKz7A4OTnx3nvv0bhxY+rXr8+8efM4dOgQa9asueB7HDx4ENM0iYyMtM1LTEykoKCAfv36ER0dTePGjRk1ahQ+Pj62dapWrQpw2bM+ImIfLvYuQEQEoGHDhjg5/fnvq7CwMBo1amSbdnZ2JigoiOTkZAB+/vln9u3bh6+vb5H95OTksH///gu+R3Z2NgAeHh62eU2bNqVr1640btyYW2+9lW7dunHnnXdSpUoV2zqenp5A4dkgEal4FGZEpEJwdXUtMm0YxgXnWa1WAKxWKy1atOCTTz45b18hISEXfI/g4GCg8HLTuXWcnZ1ZuXIlGzZsYMWKFbz22mtMnjyZzZs327qXTp06dcn9ioh96TKTiDik66+/nr179xIaGkrt2rWLvC7WLVWrVi38/PzYvXt3kfmGYXDjjTcydepUfvnlF9zc3Fi4cKFt+c6dO3F1daVhw4Zl+plE5MoozIiIQ7r33nsJDg6mT58+rFu3jvj4eGJjYxk3bhxHjhy54DZOTk7cfPPNrF+/3jZv8+bNTJs2ja1bt3Lo0CEWLFjAiRMnqF+/vm2ddevW0b59e9vlJhGpWBRmRMQheXl5sXbtWqpXr06/fv2oX78+I0aMIDs7Gz8/v4tu9+CDDzJ//nzb5So/Pz/Wrl1Lz549qVu3Lv/+9795+eWX6dGjh22bzz77jAceeKDMP5OIXBk9NE9ErimmaXLDDTcwfvx47rnnnsuuv2TJEh577DF27NiBi4tuMxSpiHRmRkSuKYZh8M4771BQUFCs9TMzM5k3b56CjEgFpjMzIiIi4tB0ZkZEREQcmsKMiIiIODSFGREREXFoCjMiIiLi0BRmRERExKEpzIiIiIhDU5gRERERh6YwIyIiIg5NYUZEREQc2v8DL7xmV30L1qIAAAAASUVORK5CYII=", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], "source": [ "plt.plot(time, fuel1, label=\"Stage 1 fuel\")\n", "plt.plot(time, fuel2, label=\"Stage 2 fuel\")\n", @@ -985,10 +1016,21 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 43, "id": "b1c58307-4649-4236-bd2a-e0f0893515bd", "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAj4AAAHFCAYAAADyj/PrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABR/ElEQVR4nO3de1xUdf4/8NeZK3eUO3i/oKXgJS2VSjHT1nuZtV6+ZdllSytNrc3cLexr6NrmWlnutpVp37z8Kmsr1wubipqX1DQVzdU0wBRRQG4DM8PM5/cHzoFhBoRhhrm9no8Hj5Uzh+HNWS+v3p+bJIQQICIiIvIDCncXQERERNRSGHyIiIjIbzD4EBERkd9g8CEiIiK/weBDREREfoPBh4iIiPwGgw8RERH5DQYfIiIi8hsMPkREROQ3GHyIPMDbb78NSZKQlJTk7lJkqampSE1Nddv3X7t2LZYvX273NUmSkJaW1qL1eBOdToe0tDTs3LnT5rWPP/4YkiTh119/bfG6iDwBgw+RB/joo48AAFlZWThw4ICbq/EMDQWfffv24fHHH2/ZgryITqfDwoUL7Qaf0aNHY9++fYiPj2/5wog8AIMPkZsdOnQIP/30E0aPHg0A+PDDD91ckWtUVFQ47b0GDhyItm3bOu39vJHRaERVVVWTvy46OhoDBw6EVqt1QVVEno/Bh8jNLEFnyZIlSElJwfr166HT6Wzu++233/Dkk0+iXbt20Gg0SEhIwMSJE3H58mX5nmvXrmHu3Lno3LkztFotYmJiMGrUKPz888/yPQaDAYsWLcJNN90ErVaL6OhoPProo7hy5coNa23s13bs2BFjxozBxo0b0bdvXwQEBGDhwoUAgHfffReDBw9GTEwMgoODkZycjKVLl8JoNMpfn5qaik2bNiE7OxuSJMkfFvaGuk6cOIHx48ejdevWCAgIQJ8+fbB69Wqre3bu3AlJkrBu3TosWLAACQkJCAsLw913343Tp0/f8OcHgD179mDYsGEIDQ1FUFAQUlJSsGnTJvn1n376CZIk2Q2wmzdvhiRJ+Prrr+VrZ86cwZQpUxATEwOtVoubb74Z7777rt26P/nkE8ydOxdt2rSBVqvF2bNnbb7Hr7/+iujoaADAwoUL5Wf3yCOPALA/1JWamoqkpCTs27cPKSkpCAwMRMeOHbFq1SoAwKZNm3DLLbcgKCgIycnJ2LJli833bczPQeQRBBG5jU6nE+Hh4eLWW28VQgjxwQcfCADi448/trrvwoULIj4+XkRFRYlly5aJ//znP2LDhg1i+vTp4tSpU0IIIUpKSkTPnj1FcHCweO2118TWrVvFF198IWbNmiW2b98uhBDCZDKJ3/3udyI4OFgsXLhQZGRkiA8++EC0adNG9OjRQ+h0Ovl7DhkyRAwZMkT+vClf26FDBxEfHy86d+4sPvroI7Fjxw7xww8/CCGEeP7558XKlSvFli1bxPbt28Xf/vY3ERUVJR599FH567OyssTtt98u4uLixL59++QPCwDi1VdflT//+eefRWhoqOjSpYtYs2aN2LRpk5g8ebIAIP7yl7/I9+3YsUMAEB07dhRTp04VmzZtEuvWrRPt27cXiYmJoqqqqsH/v3bu3CnUarXo16+f2LBhg/jqq6/EiBEjhCRJYv369fJ9ffv2FbfffrvN1z/44IMiJiZGGI1G+ecMDw8XycnJYs2aNWLbtm1i7ty5QqFQiLS0NJu627RpIyZOnCi+/vpr8e2334qCggKb71FZWSm2bNkiAIjHHntMfnZnz54VQgixatUqAUCcP39e/pohQ4aIyMhI0b17d/Hhhx+KrVu3ijFjxggAYuHChSI5OVmsW7dO/Pvf/xYDBw4UWq1W/Pbbb1b/fzXm5yDyBAw+RG60Zs0aAUD8/e9/F0IIUVpaKkJCQsSdd95pdd/06dOFWq0WJ0+erPe9XnvtNQFAZGRk1HvPunXrBADxxRdfWF0/ePCgACDee+89+Vrd4NOUr+3QoYNQKpXi9OnT9f/wojpMGY1GsWbNGqFUKkVhYaH82ujRo0WHDh3sfl3d4DNp0iSh1WpFTk6O1X0jR44UQUFB4tq1a0KImgAxatQoq/v+3//7fwKAVbiyZ+DAgSImJkaUlpbK16qqqkRSUpJo27atMJvNQggh3n77bQHA6ucvLCwUWq1WzJ07V752zz33iLZt24ri4mKr7/PMM8+IgIAA+XlY6h48eHCD9VlcuXLF5hlZ1Bd8AIhDhw7J1woKCoRSqRSBgYFWIefo0aMCgHj77beb/HMQeQIOdRG50YcffojAwEBMmjQJABASEoIHHngAu3fvxpkzZ+T7Nm/ejKFDh+Lmm2+u9702b96Mbt264e677673nm+//RatWrXC2LFjUVVVJX/06dMHcXFxdifDOvq1vXr1Qrdu3Wze58iRIxg3bhwiIyOhVCqhVqvx8MMPw2Qy4b///W+9378h27dvx7Bhw9CuXTur64888gh0Oh327dtndX3cuHE2tQJAdnZ2vd+jvLwcBw4cwMSJExESEiJfVyqVeOihh3DhwgV5uGzq1KnQarX4+OOP5fvWrVsHvV6PRx99FABQWVmJ7777Dvfddx+CgoKsnumoUaNQWVmJ/fv3W9Vw//33N/KJNF18fDz69esnfx4REYGYmBj06dMHCQkJ8nXL70HLs3Lk5yByJwYfIjc5e/Ysdu3ahdGjR0MIgWvXruHatWuYOHEigJqVXgBw5cqVG07mbcw9ly9fxrVr16DRaKBWq60+8vLycPXqVad9rb1VQzk5Objzzjvx22+/4a233sLu3btx8OBBeS6IoxOgCwoK7H4/yz/YBQUFVtcjIyOtPrdM9G3o+xcVFUEI0ajvExERgXHjxmHNmjUwmUwAqufW3HbbbejZs6d8b1VVFd555x2b5zlq1CgAaNQzdZaIiAibaxqNxua6RqMBUB14AMd+DiJ3Urm7ACJ/9dFHH0EIgc8//xyff/65zeurV6/GokWLoFQqER0djQsXLjT4fo25JyoqCpGRkXYnpwJAaGio07629mRki6+++grl5eXYuHEjOnToIF8/evRog3XfSGRkJC5dumRz/eLFiwCqa2+u1q1bQ6FQNPr7PProo/jss8+QkZGB9u3b4+DBg1i5cqXV+1m6RTNnzrT7PTt16mT1ub1n6m6O/BxE7sTgQ+QGJpMJq1evRpcuXfDBBx/YvP7tt9/izTffxObNmzFmzBiMHDkSn3zyCU6fPo3u3bvbfc+RI0filVdewfbt23HXXXfZvWfMmDFYv349TCYTBgwY0KSam/O1FpZ/uGsvpRZC4J///KfNvVqtttEdoGHDhuHLL7/ExYsXrYZl1qxZg6CgIAwcONChemsLDg7GgAEDsHHjRvz1r39FYGAgAMBsNuP//u//0LZtW6uhvREjRqBNmzZYtWoV2rdvj4CAAEyePFl+PSgoCEOHDsWRI0fQq1cvuZPiDI3pYDmLK38OIldg8CFyg82bN+PixYv4y1/+Ynd35KSkJKxYsQIffvghxowZg9deew2bN2/G4MGD8fLLLyM5ORnXrl3Dli1bMGfOHNx0002YPXs2NmzYgPHjx+Oll17CbbfdhoqKCmRmZmLMmDEYOnQoJk2ahE8//RSjRo3CrFmzcNttt0GtVuPChQvYsWMHxo8fj/vuu89uzc35Wovhw4dDo9Fg8uTJePHFF1FZWYmVK1eiqKjI5t7k5GRs3LgRK1euRL9+/aBQKNC/f3+77/vqq6/i22+/xdChQ/HKK68gIiICn376KTZt2oSlS5ciPDz8xv+nNMLixYsxfPhwDB06FPPmzYNGo8F7772HEydOYN26dVYdGaVSiYcffhjLli1DWFgYJkyYYFPHW2+9hTvuuAN33nknnn76aXTs2BGlpaU4e/YsvvnmG2zfvt2hOkNDQ9GhQwf861//wrBhwxAREYGoqCh07NixOT9+vVz1cxC5hHvnVhP5p3vvvVdoNBqRn59f7z2TJk0SKpVK5OXlCSGEyM3NFdOnTxdxcXFCrVaLhIQE8eCDD4rLly/LX1NUVCRmzZol2rdvL9RqtYiJiRGjR48WP//8s3yP0WgUf/3rX0Xv3r1FQECACAkJETfddJP4wx/+IM6cOSPfV3dVV1O+tkOHDmL06NF2f65vvvlG/vo2bdqIF154QWzevFkAEDt27JDvKywsFBMnThStWrUSkiSJ2n9dwc6KpePHj4uxY8eK8PBwodFoRO/evcWqVaus7rGsjvrss8+srp8/f14AsLnfnt27d4u77rpLBAcHi8DAQDFw4EDxzTff2L33v//9rwDQ4Gq78+fPi+nTp4s2bdoItVotoqOjRUpKili0aNEN627If/7zH9G3b1+h1WoFADFt2jQhRP2runr27GnzHvX9/whAzJw5s8k/B5EnkIQQwk2Zi4iIiKhFcVUXERER+Q0GHyIiIvIbDD5ERETkNxh8iIiIyG8w+BAREZHfYPAhIiIiv8ENDOswm824ePEiQkNDPXJ7eCIiIrIlhEBpaSkSEhKgUNTf12HwqePixYs2JzwTERGRd8jNzW3wwGYGnzosBy3m5uYiLCzMeW9cUAC8+y4QEAAEBTnvfcn76XRAZSUwcyZQ59Rwcq4CXQHe/eFdBKgDEKTmn0OilqYz6lBprMTM22YiMsi5f9+VlJSgXbt2DR62DDD42LAMb4WFhTk3+BgMgFYLREQAISHOe1/yfmVl1cE4NBRw5u85smFQGaAN1iIiKAIhGv45JGppZYYyFOgKEBoWirAg1/x9d6NpKpzcTERERH6DwYeIiIj8BoMPERER+Q0GHyIiIvIbDD5ERETkNxh8iIiIyG8w+BAREZHfYPAhIiIiv8HgQ0RERH6DwYeIiIj8BoMPERER+Q0GHyIiIvIbDD5EPqTCYHJ3CUREHo3Bh8hH/Ovob+j56hZ8deQ3d5dCROSxGHyIfMSOn/NhFsCxC8XuLoWIyGMx+BD5iHNXywEAVWazmyshIvJcDD5EPkAIgXNXqoOP0STcXA0Rkedi8CHyAVdK9SjTVwEAqkzs+BAR1YfBh8gH/HK92wMAVWZ2fIiI6sPgQ+QDzl0tk39tZMeHiKheDD5EPuBc7Y4P5/gQEdWLwYfIB5y7UtPx4aouIqL6MfgQ+QDLUnaAq7qIiBrC4EPk5fRVJuQW6uTP2fEhIqofgw+Rl8su0KH2Qi52fIiI6sfgQ+Tlas/vAbiPDxFRQxh8iLycZQ+f8EA1AO7jQ0TUEAYfIi9nWcrePTYUAIe6iIgawuBD5OUsmxd2iwsBwKEuIqKGMPgQebHah5NaOj4c6iIiqp/XBJ+0tDRIkmT1ERcXJ78uhEBaWhoSEhIQGBiI1NRUZGVlubFiItcrLDeguMIIAOgaYxnqYseHiKg+XhN8AKBnz564dOmS/HH8+HH5taVLl2LZsmVYsWIFDh48iLi4OAwfPhylpaVurJjItSwbF7ZpFYjQABUAHllBRNQQrwo+KpUKcXFx8kd0dDSA6m7P8uXLsWDBAkyYMAFJSUlYvXo1dDod1q5d6+aqiVzHspS9c3QwVEoJADcwJCJqiFcFnzNnziAhIQGdOnXCpEmTcO7cOQDA+fPnkZeXhxEjRsj3arVaDBkyBHv37nVXuUQuZ5nf0yU6BCpF9R9nruoiIqqfyt0FNNaAAQOwZs0adOvWDZcvX8aiRYuQkpKCrKws5OXlAQBiY2OtviY2NhbZ2dkNvq9er4der5c/LykpcX7xRC7yS62Oj9rS8eEcHyKienlN8Bk5cqT86+TkZAwaNAhdunTB6tWrMXDgQACAJElWXyOEsLlW1+LFi7Fw4ULnF0zUAiwdn85RIVApr3d8uKqLiKheXjXUVVtwcDCSk5Nx5swZeXWXpfNjkZ+fb9MFqmv+/PkoLi6WP3Jzc11WM5EzGU1m5Fw/nLRzdDDUCnZ8iIhuxGuDj16vx6lTpxAfH49OnTohLi4OGRkZ8usGgwGZmZlISUlp8H20Wi3CwsKsPoi8QU6hDlVmgUC1EnFhAXLHxywAM7s+RER2ec1Q17x58zB27Fi0b98e+fn5WLRoEUpKSjBt2jRIkoTZs2cjPT0diYmJSExMRHp6OoKCgjBlyhR3l07kEpZhrk5RwVAoJHlVF1C9iaFG0fAwLxGRP/Ka4HPhwgVMnjwZV69eRXR0NAYOHIj9+/ejQ4cOAIAXX3wRFRUVmDFjBoqKijBgwABs27YNoaGhbq6cyDVqL2UHALWipoFbZTZD470NXSIil/Ga4LN+/foGX5ckCWlpaUhLS2uZgojcTJ7YHF19Rlftjg+XtBMR2cf/JCTyUpbDSbtc7/ioag1tcYIzEZF9DD5EXqr2UnaguutpCT88qJSIyD4GHyIvVKwzoqDcAKBmjg9QM9zFg0qJiOxj8CHyQr9cH+aKCwtAsLZmqp5lgjMPKiUiso/Bh8gL/ZJvvaLLggeVEhE1jMGHyAudu2pZ0VU3+PCgUiKihjD4EHkheQ+f6xObLWqOrWDwISKyh8GHyAvV7OFTT8eHQ11ERHYx+BB5GZNZILug+nDSLtHWHR95jg87PkREdjH4EHmZC0U6GExmaFQKJLQKtHpNxRPaiYgaxOBD5GXkw0kjg6GscxCpSmEZ6mLHh4jIHgYfIi/zyxX7S9kBQK1kx4eIqCEMPkRexrKUve78HoDL2YmIboTBh8jLnGug41NzVhc7PkRE9jD4EHmZX+Sl7LYdH7WSR1YQETWEwYfIi5RWGnGlVA+gno4PDyklImoQgw+RF7Gs6IoK0SIsQG3zumVVVxVXdRER2cXgQ+RFzl2tf34PwFVdREQ3wuBD5EUsHZ8u9QQfruoiImoYgw+RF5HP6IqyndgM1DqklKu6iIjsYvAh8iINbV4I1J7czI4PEZE9DD5EXsJsFvi1oP6l7EDNUBeXsxMR2cfgQ+QlLhZXoNJohlopoV3rQLv3cKiLiKhhDD5EXsIyv6dDZLDc2amLk5uJiBrG4EPkJeSjKqLsz+8Baub4cDk7EZF9DD5EXqKhoyos1NzAkIioQQw+RF7iRpsXAjyygojoRhh8iLzEjTYvBHhIKRHRjTD4EHkBnaEKl4orAdS/eSEAqK6v6jJyVRcRkV0MPkRewNLtaR2kRutgTb33cR8fIqKGMfgQeYFzV288sRmodUgpOz5ERHYx+BB5gcYsZQcAlYL7+BARNYTBh8gLnPitGADQPS60wfu4jw8RUcMYfIg8nBACR3KuAQD6tm/d4L01Q13s+BAR2cPgQ+ThLhRVoKDcALVSQs+EsAbvrRnqYseHiMgeBh8iD/djThEAoEdCOALUygbvlTs+nONDRGQXgw+Rh5OHudq1uuG9cseHQ11ERHYx+BB5uCO51wAAfdu3uuG9nNxMRNQwBh8iD1ZpNOHkxeoVXX3bNTyxGeCRFUREN8LgQ+TBTl4qgdEkEBmsQbuIwBvezyMriIgaxuBD5MFqlrG3giRJN7yfR1YQETWMwYfIgx25vqLrRvv3WKg5x4eIqEFeG3wWL14MSZIwe/Zs+ZoQAmlpaUhISEBgYCBSU1ORlZXlviKJmsnS8enTiBVdAFd1ERHdiFcGn4MHD+L9999Hr169rK4vXboUy5Ytw4oVK3Dw4EHExcVh+PDhKC0tdVOlRI7LL63Eb9cqIElAr7bhjfoadnyIiBrmdcGnrKwMU6dOxT//+U+0bl3T/hdCYPny5ViwYAEmTJiApKQkrF69GjqdDmvXrnVjxUSOOXq929MtJhShAepGfQ3n+BARNczrgs/MmTMxevRo3H333VbXz58/j7y8PIwYMUK+ptVqMWTIEOzdu7fe99Pr9SgpKbH6IPIETdm/x8KyqotndRER2adydwFNsX79evz44484ePCgzWt5eXkAgNjYWKvrsbGxyM7Orvc9Fy9ejIULFzq3UCInqJnY3KrRXyPv48Pl7EREdnlNxyc3NxezZs3C//3f/yEgIKDe++ou+RVCNLgMeP78+SguLpY/cnNznVYzkaNMZoFjF6o3LuzTiI0LLSw7NxtNAkKw60NEVJfXdHwOHz6M/Px89OvXT75mMpmwa9curFixAqdPnwZQ3fmJj4+X78nPz7fpAtWm1Wqh1WpdVziRA/57uRQ6gwkhWhW6xoQ0+uvUipr/ljGZhRyEiIiomtd0fIYNG4bjx4/j6NGj8kf//v0xdepUHD16FJ07d0ZcXBwyMjLkrzEYDMjMzERKSoobKydqOssy9t7twqFUND681A46nOdDRGTLazo+oaGhSEpKsroWHByMyMhI+frs2bORnp6OxMREJCYmIj09HUFBQZgyZYo7SiZymDy/pwnDXIB18DGazAhQK51aFxGRt/Oa4NMYL774IioqKjBjxgwUFRVhwIAB2LZtG0JDQ91dGlGTWFZ0NXbjQovaQ11c0k5EZMurg8/OnTutPpckCWlpaUhLS3NLPUTOUFxhxNn8MgBAnyas6AIAhUKCQgLMggeVEhHZ4zVzfIj8xbEL1wAA7SOCEBXS9In33MSQiKh+DD5EHqb2ieyOkDcxZPAhIrLB4EPkYWomNrdy6OstwYdDXUREthh8iDyIEAJHLROb2zdtRZeFmkNdRET1YvAh8iDZBToU6YzQqBToER/m0HvU7N7Mjg8RUV0MPkQe5Ehu9TBXUkIYNCrH/niqFJbzutjxISKqi8GHyIPUTGx2bJgLANRKy+RmdnyIiOpi8CHyIJbg09SNC2uzLGc3co4PEZENBh8iD1FpNOPUpRIAji9lB2otZ+eqLiIiGww+RB7iRH45qswC0aFatGkV6PD7cFUXEVH9GHyIPMSRi9XHVPRt1wqS1PgT2eviqi4iovox+BB5iCOXrgefZkxsBmoOKuWqLiIiWww+RB7i6PXg05yJzQA7PkREDWHwIfIAeUKDi6UGKCSgV9vwZr0XDyklIqofgw+RBziKEABA97gwBGtVzXovNVd1ERHVi8GHyAMcEaEAmreM3aJmqIsdHyKiuhh8iDzAEVHd8XH0RPbaaoa62PEhIqqLwYfIzaoEcOz6UJczOj41Q13s+BAR1cXgQ+RmP1dpUQklQrVKdI4Kafb78cgKIqL6MfgQudkRY/UuzX3iQ6BQOL5xoQUPKSUiqh+DD5Gb7b8efPrGN7/bAwCq6xsYGjnURURkg8GHyI2MAthlCAYADOnUvP17LFTs+BAR1YvBh8iNDldqUCqUiIARfeKc0/GRDyllx4eIyAaDD5EbbdcFAABSpSIonTC/BwBUCh5ZQURUHwYfIjeyBJ+hKHLae/LICiKi+jH4ELlJjlGJs0Y1lBAYLF1z2vvyyAoiovox+BC5iaXb019dgXDJ5LT35T4+RET1a/JpiMXFxfjyyy+xe/du/Prrr9DpdIiOjkbfvn1xzz33ICUlxRV1EvkcS/AZpi0DKpz3vtzHh4iofo3u+Fy6dAlPPPEE4uPj8dprr6G8vBx9+vTBsGHD0LZtW+zYsQPDhw9Hjx49sGHDBlfWTOT1ys0S9ldoAQB3acqd+t7y5Gau6iIistHojk/v3r3x8MMP44cffkBSUpLdeyoqKvDVV19h2bJlyM3Nxbx585xWKJEv+b5CCwMktFNVoYvS4NT35iGlRET1a3TwycrKQnR0dIP3BAYGYvLkyZg8eTKuXLnS7OKIfNUOyzBXUCUk56xil9UMdbHjQ0RUV6OHum4Uepp7P5G/EKLWMvagSqe/P4+sICKqX5NWdW3duhWTJ0/GuXPnAACPPfaYS4oi8mVZBjUum5QIlMwYEKB3+vvzyAoiovo1KfjMmzcPY8aMwaOPPooLFy7g5MmTrqqLyGdZhrnuCNQjwAUbSqi5gSERUb2atJw9PDwcU6dOxcCBA/HEE0+gqqrKVXUR+azvrgefu1wwzAXUXtXFjg8RUV1N+u/NkJDqQxS7dOmCmTNn4scff3RJUUS+6qpJgZ/0agCumd8DsONDRNSQJgWfv//97zCZqneYHTNmDA4dOuSSooh81U5dAAQk9NQYEKdyTUfGMseHh5QSEdlq0lBXx44dAVTv1yOEQN++fQEA2dnZ+PLLL9GjRw+MGDHC6UUS+YoduuubFrqo2wPUrOqq4qouIiIbDk2tHD9+PNasWQMAuHbtGgYMGIA333wT48ePx8qVK51aIJGvMApglwuXsVvwyAoiovo5FHx+/PFH3HnnnQCAzz//HLGxscjOzsaaNWvw9ttvO7VAIl9xsFKDUqFApMKE3lqjy74PDyklIqqfQ8FHp9MhNDQUALBt2zZMmDABCoUCAwcORHZ2tlMLJPIVlmXsQ4L0UDp5t+baLKu6qriqi4jIhkPBp2vXrvjqq6+Qm5uLrVu3yvN68vPzERYW5tQCLVauXIlevXohLCwMYWFhGDRoEDZv3iy/LoRAWloaEhISEBgYiNTUVGRlZbmkFiJHyKexBznxKHY7LKu6TJzjQ0Rkw6Hg88orr2DevHno2LEjBgwYgEGDBgGo7v5YJjw7W9u2bbFkyRIcOnQIhw4dwl133YXx48fL4Wbp0qVYtmwZVqxYgYMHDyIuLg7Dhw9HaWmpS+ohaopsoxK/GNVQQeDOIOfv1lxbzaouBh8iorocCj4TJ05ETk4ODh06hC1btsjXhw0bhr/97W9OK662sWPHYtSoUejWrRu6deuG119/HSEhIdi/fz+EEFi+fDkWLFiACRMmICkpCatXr4ZOp8PatWtdUg9RU1i6Pf0DDAhTuDaQqBU8nZ2IqD5NCj4JCQl4+umnsXnzZkRERKBv375QKGre4rbbbsNNN93k9CLrMplMWL9+PcrLyzFo0CCcP38eeXl5VkvptVothgwZgr1797q8HqIb2V7rNHZXkzs+HOoiIrLRpOCzdu1aBAUF4bnnnkNUVBQeeOABfPLJJygsLHRVfVaOHz+OkJAQaLVaPPXUU/LeQXl5eQCA2NhYq/tjY2Pl1+qj1+tRUlJi9UHkTOVmCQcqqvfvceUydgseUkpEVL8mBZ/U1FS8+eabOHPmDPbt24dbbrkF7777LuLj45Gamoq//e1v+OWXX1xVK7p3746jR49i//79ePrppzFt2jSrg1IlyXqpjBDC5lpdixcvRnh4uPzRrl07l9RO/mtPhRYGSGivqkIXtevPt7NsYGgWgJldHyIiKw6fDd2zZ0/Mnz8f+/fvR3Z2NqZOnYrt27cjOTkZSUlJ2LRpkzPrBABoNBp07doV/fv3x+LFi9G7d2+89dZbiIuLAwCb7k5+fr5NF6iu+fPno7i4WP7Izc11et3k37bXOpT0BjncKVS11srzoFIiImsOB5/a4uLi8MQTT+Cbb77B1atX8b//+7/QarXOeOsGCSGg1+vRqVMnxMXFISMjQ37NYDAgMzMTKSkpDb6HVquVl8hbPoicRYia/XtceUxFbepa8+54UCkRkbUmndVVV35+PvLz82Gu81+V9913X7OKsufll1/GyJEj0a5dO5SWlmL9+vXYuXMntmzZAkmSMHv2bKSnpyMxMRGJiYlIT09HUFAQpkyZ4vRaiBory6BGvkmJIMmMAYGuXcZuUbvjw+BDRGTNoeBz+PBhTJs2DadOnYIQ1n+xSpIkn+DuTJcvX8ZDDz2ES5cuITw8HL169cKWLVswfPhwAMCLL76IiooKzJgxA0VFRRgwYAC2bdsm7zBN5A7/ud7tuSNQD20LDHMBNTs3AxzqIiKqy6Hg8+ijj6Jbt2748MMPERsbe8MJxM7w4YcfNvi6JElIS0tDWlqay2shagwhgC9LgwAAI4JbZpgLqP6zoFJIqDILdnyIiOpwKPicP38eGzduRNeuXZ1dD5HP+KFSg+wqFYIlM0YGu/aYirpUyurgY+SSdiIiKw5Nbh42bBh++uknZ9dC5FP+3/Vuz5iQCgS7eLfmuuTdm7mcnYjIikMdnw8++ADTpk3DiRMnkJSUBLVabfX6uHHjnFIckbcqM0v4d3kgAODBUF2Lf39uYkhEZJ9DwWfv3r3Ys2eP1enoFq6a3EzkTTaVBaJCKNBZbcQtWkOLf3/V9RPaeVApEZE1h4a6nnvuOXmFldlstvpg6CGqGeZ6MFTXIpsW1qW+vrKriqu6iIisOBR8CgoK8Pzzz99wV2Qif3TWoMJhvRZKCEwIaflhLoAdHyKi+jgUfCZMmIAdO3Y4uxYin/DZ9W7P0KBKxKjc03HhHB8iIvscmuPTrVs3zJ8/H3v27EFycrLN5ObnnnvOKcUReRujAL4oqw4+D7hhUrMFV3UREdnn8KqukJAQZGZmIjMz0+o1SZIYfMhvZeoCcNWkRJTS1GJnc9lj6fhwHx8iImsOb2BIRLYsk5rvC9FB7YZJzRaWOT7cuZmIyJpTTmcnIuBKlQLbr5/N5c5hLoCruoiI6tPo4LNkyRLodI37y/zAgQPYtGmTw0UReaOvyoJQBQm9tQZ001S5tZaaoS52fIiIamt08Dl58iTat2+Pp59+Gps3b8aVK1fk16qqqnDs2DG89957SElJwaRJkxAWFuaSgok8kRC19+4pd3M1gNoy1MWODxGRlUbP8VmzZg2OHTuGd999F1OnTkVxcTGUSiW0Wq3cCerbty+efPJJTJs2DVqt1mVFE3mao3o1zhjV0EoCY0Na9kBSe1QKdnyIiOxp0uTmXr164R//+Af+/ve/49ixY/j1119RUVGBqKgo9OnTB1FRUa6qk8ijfVYaDAAYFVyBsBY+kNQeTm4mIrLPoVVdkiShd+/e6N27t7PrIfI6FWYJ35RVH0j6gAcMcwGAWsnJzURE9jgUfADAbDbj7NmzyM/Ph7nOX66DBw9udmFE3mJLeQBKhQLtVFUYGNDyB5Lao1LwyAoiInscCj779+/HlClTkJ2dDSGs/2Ll6ezkb/7f9WGuB0J1ULhx757aeGQFEZF9DgWfp556Cv3798emTZsQHx8PyR3HTxN5gByjEvsqtZAgcL+b9+6pjUdWEBHZ51DwOXPmDD7//HN07drV2fUQeZXPry9hvyNQjzYqz+l08sgKIiL7HNq5ecCAATh79qyzayHyKiZRE3we9KBuD1BrHx/O8SEisuJQx+fZZ5/F3LlzkZeXZ/d09l69ejmlOCJP9n2FFhdNKoQrzBge5P69e2qT9/Hhqi4iIisOBZ/7778fADB9+nT5miRJEEJwcjP5jTUl1ZOax4foEOBhp95xHx8iIvt4OjuRA342qPAfXSAkCEwL84y9e2pTc1UXEZFdDgWfDh06OLsOIq/yblEogOqdmru4+UBSe+R9fLiqi4jIisMbGP7yyy9Yvnw5Tp06BUmScPPNN2PWrFno0qWLM+sj8jjnDCp8W169U/PMVqVursY+7uNDRGSfQzMTtm7dih49euCHH35Ar169kJSUhAMHDqBnz57IyMhwdo1EHmXltRAISBgWVIEeWs/r9gC1h7rY8SEiqs2hjs9LL72E559/HkuWLLG5/sc//hHDhw93SnFEnuaCUYkvy6qXsHtqtwfgUBcRUX0c6vicOnUKjz32mM316dOn4+TJk80uishT/aM4BFWQcHtgJW4JMLq7nHpxcjMRkX0OBZ/o6GgcPXrU5vrRo0cRExPT3JqIPFJ+lQIbrp/L5cndHqBmOTsPKSUisubQUNcTTzyBJ598EufOnUNKSgokScKePXvwl7/8BXPnznV2jUQe4Z/FITAICf20egzykFPY62PZwLCKGxgSEVlxKPj8+c9/RmhoKN58803Mnz8fAJCQkIC0tDQ899xzTi2QyBMUmRT49PqGhc+0LoWnn8vLIyuIiOxzKPhIkoTnn38ezz//PEpLq1v+oaGhTi2MyJOsKg6GTijQU2NAaqDe3eXcEA8pJSKyz+F9fCwYeMjXlZglrCoJAQA808rzuz1AzaquKq7qIiKy0ujgc8stt+C7775D69at0bdvX0gN/O3/448/OqU4Ik/wSXEwSs0KdFUbcU9wpbvLaRSu6iIisq/RwWf8+PHQarXyrxsKPkS+QmeW8GFxdbdnRqtSKLzktz1XdRER2dfo4PPqq6/Kv05LS3NFLUQeZ11pEArNSrRTVWFcSIW7y2k0NVd1ERHZ5dA+Pp07d0ZBQYHN9WvXrqFz587NLorIE+gF8P616jlsT7cqhcpLuj1ATceHc3yIiKw5FHx+/fVXmEwmm+t6vR4XLlxodlFEnuCL0iBcNikRpzTh/lCdu8tpEhXP6iIisqtJq7q+/vpr+ddbt25FeHi4/LnJZMJ3332HTp06Oa86IjepEsDK692eJ1uVQutF3R6g1gaGnNxMRGSlScHn3nvvBVC9j8+0adOsXlOr1ejYsSPefPNNpxVH5C7/KgtEbpUKkQoTJntZtwfgIaVERPVpUvAxX58o2alTJxw8eBBRUVEuKYrInXRmCW8WhgEAHmtVhkCF94UHLmcnIrLPoTk+58+fb/HQs3jxYtx6660IDQ1FTEwM7r33Xpw+fdrqHiEE0tLSkJCQgMDAQKSmpiIrK6tF6yTv9961EFw0qdBGVYXpYWXuLschKh5ZQURkl8M7N5eXlyMzMxM5OTkwGKwPbHTFeV2ZmZmYOXMmbr31VlRVVWHBggUYMWIETp48ieDg6jOUli5dimXLluHjjz9Gt27dsGjRIgwfPhynT5/mDtPUKL8alfJKrj9HFiPAof80cD/LHB8jl7MTEVlxKPgcOXIEo0aNgk6nQ3l5OSIiInD16lUEBQUhJibGJcFny5YtVp+vWrUKMTExOHz4MAYPHgwhBJYvX44FCxZgwoQJAIDVq1cjNjYWa9euxR/+8Aen10S+57WCVjBAwp2BlbgnyDt2abaHh5QSEdnn0H/PPv/88xg7diwKCwsRGBiI/fv3Izs7G/369cNf//pXZ9doV3FxMQAgIiICQPXwW15eHkaMGCHfo9VqMWTIEOzdu7fe99Hr9SgpKbH6IP/0XXkAtusCoIZAWmSxV5zJVR95ObtZQAiGHyIiC4eCz9GjRzF37lwolUoolUro9Xq0a9cOS5cuxcsvv+zsGm0IITBnzhzccccdSEpKAgDk5eUBAGJjY63ujY2NlV+zZ/HixQgPD5c/2rVr57rCyWNVmoHXCqq3Z5geXoYumio3V9Q8akXNH21uYkhEVMOh4KNWq+WzumJjY5GTkwMACA8Pl3/tSs888wyOHTuGdevW2bxW9wwxIUSD54rNnz8fxcXF8kdubq7T6yXP90FxCLKrVIhRmvBs61J3l9Nslo4PwOEuIqLaHJrj07dvXxw6dAjdunXD0KFD8corr+Dq1av45JNPkJyc7OwarTz77LP4+uuvsWvXLrRt21a+HhcXB6C68xMfHy9fz8/Pt+kC1abVauXDV8k//ValxIrrE5oXRBYjxAuXr9dVO/gYzWYEQunGaoiIPIdDHZ/09HQ5XPzv//4vIiMj8fTTTyM/Px/vv/++Uwu0EELgmWeewcaNG7F9+3abHaI7deqEuLg4ZGRkyNcMBgMyMzORkpLikprIN6QXhKFSKHBbgB7jgr3nINKGWA11seNDRCRrcsdHCIHo6Gj07NkTABAdHY1///vfTi+srpkzZ2Lt2rX417/+hdDQUHneTnh4OAIDAyFJEmbPno309HQkJiYiMTER6enpCAoKwpQpU1xeH3mn7yu02FQeBAUEFkZe8+oJzbUpFBIUEmAW3MSQiKg2h4JPYmIisrKykJiY6Iqa7Fq5ciUAIDU11er6qlWr8MgjjwAAXnzxRVRUVGDGjBkoKirCgAEDsG3bNu7hQ3YZBfDq1eoJzQ+FleNmrXdPaK5LpVTAUGXmsRVERLU0OfgoFAokJiaioKCgRYNPY5bkSpKEtLQ0pKWlub4g8nqri4Nx1qhGhMKEOa19bxsDtUKCAez4EBHV5tAcn6VLl+KFF17AiRMnnF0PUYvIr1JgeVH1eVx/jChBuNL3uiKWYyuMnONDRCRzaFXX//zP/0Cn06F3797QaDQIDAy0er2wsNApxRG5ypLCcJQJBXprDXjAC09fbwz5oFIeW0FEJHMo+CxfvtzJZRC1nMOVGmwsC4IEgdcir0HhIxOa61IpeGwFEVFdDgWfadOmObsOohaztiQIADAxRIfeAUY3V+M6lr18jJzjQ0Qkc/js6V9++QV/+tOfMHnyZOTn5wOoPkg0KyvLacURuUKhqXozv1sDDW6uxLXkg0q5qouISOZQ8MnMzERycjIOHDiAjRs3oqysDABw7NgxvPrqq04tkMjZSszVnZAwhW93QlQKdnyIiOpyKPi89NJLWLRoETIyMqDRaOTrQ4cOxb59+5xWHJErlJirf9v7fPBRco4PEVFdDgWf48eP47777rO5Hh0djYKCgmYXReRKpXLw8e1AwFVdRES2HAo+rVq1wqVLl2yuHzlyBG3atGl2UUSu5H9DXb4d8IiImsKh4DNlyhT88Y9/RF5eHiRJgtlsxvfff4958+bh4YcfdnaNRE5jFIBOVP+2D/Xxjg+HuoiIbDkUfF5//XW0b98ebdq0QVlZGXr06IHBgwcjJSUFf/rTn5xdI5HTlJlrfsuH+njHh0NdRES2HNrHR61W49NPP8Vrr72GI0eOwGw2o2/fvi16dheRIyzDXMGSGSof3bjQwrKBIYe6iIhqOBR8LLp06YIuXbo4qxYil7Os6PL1YS6gVseHy9mJiGSNDj5z5sxp9JsuW7bMoWKIXK3UTyY2A7U6PtzAkIhI1ujgc+TIkUbdJ0k+Pn5AXs1f9vABao6sYMeHiKhGo4PPjh07XFkHUYsoMfnTUBdXdRER1eXwWV0AcPbsWWzduhUVFRUAACH4Fyx5NnkPH6Xvd0HkfXy4qouISOZQ8CkoKMCwYcPQrVs3jBo1St7M8PHHH8fcuXOdWiCRM/nXUBc7PkREdTkUfJ5//nmo1Wrk5OQgKChIvv773/8eW7ZscVpxRM7GVV1ERP7NoeXs27Ztw9atW9G2bVur64mJicjOznZKYUSuwFVdRET+zaGOT3l5uVWnx+Lq1avQarXNLorIVfxpqIsdHyIiWw4Fn8GDB2PNmjXy55bzut544w0MHTrUacUROZs/DXVZlrNz52YiohoODXW98cYbSE1NxaFDh2AwGPDiiy8iKysLhYWF+P77751dI5HT+ONQF8/qIiKq4VDHp0ePHjh27Bhuu+02DB8+HOXl5ZgwYQKOHDnCIyzIo/nnUBc7PkREFg6f1RUXF4eFCxc6sxYil/OvoS4eUkpEVJdDHZ9Vq1bhs88+s7n+2WefYfXq1c0uisgVzAIo86uhrusdHw51ERHJHAo+S5YsQVRUlM31mJgYpKenN7soIlcoFxLM8J/gwyMriIhsORR8srOz0alTJ5vrHTp0QE5OTrOLInIFyzCXBgIBzTqsxTvUrOry/ZBHRNRYDv31HxMTg2PHjtlc/+mnnxAZGdnsoohcodSPzukCALW8qosdHyIiC4eCz6RJk/Dcc89hx44dMJlMMJlM2L59O2bNmoVJkyY5u0Yip7CczO4Pw1wAOz5ERPY4tKpr0aJFyM7OxrBhw6BSVb+F2WzGww8/zDk+5LH8aUUXACgVXM5ORFSXQ8FHo9Fgw4YNWLRoEY4ePYrAwEAkJyejQ4cOzq6PyGn8afNCoNbkZq7qIiKSObyPD1B9KGliYqKzaiFyKX/avBCoWc7OfXyIiGo4NMdn4sSJWLJkic31N954Aw888ECziyJyBX8b6rJ0fEyc3ExEJHMo+GRmZmL06NE213/3u99h165dzS6KyBX8baiLk5uJiGw5FHzKysqg0WhsrqvVapSUlDS7KCJX8L+hLi5nJyKqy6Hgk5SUhA0bNthcX79+PXr06NHsoohcwdLx8Z+hLsuqLv8IekREjeHQ5OY///nPuP/++/HLL7/grrvuAgB89913WLt2LT7//HOnFkjkLH7X8eEhpURENhwKPuPGjcNXX32F9PR0fP755wgMDETv3r2xfft2hIWFObtGIqeoCT7+EQR4SCkRkS2Hl7OPHj1anuB87do1fPrpp5g9ezZ++uknmEwmpxVI5Cw1Q13+EQR4SCkRka1mHdW4fft2/M///A8SEhKwYsUKjBo1CocOHXJWbUROJXd8/OSsLq7qIiKy1eSOz4ULF/Dxxx/jo48+Qnl5OR588EEYjUZ88cUXnNhMbmUSwDvXQrGlPBArYgrRVVNl9XrNWV3+0QHhIaVERLaa1PEZNWoUevTogZMnT+Kdd97BxYsX8c4777iqNhu7du3C2LFjkZCQAEmS8NVXX1m9LoRAWloaEhISEBgYiNTUVGRlZbVYfeQ+hSYFHsmLxPKiMPxsUOP7Cq3V65VmwAD/GupSKXlWFxFRXU0KPtu2bcPjjz+OhQsXYvTo0VAqla6qy67y8nL07t0bK1assPv60qVLsWzZMqxYsQIHDx5EXFwchg8fjtLS0hatk1rWT5VqjP0tGrsrAuq9xzLMpYBAsOQfQUAe6uLkZiIiWZOCz+7du1FaWor+/ftjwIABWLFiBa5cueKq2myMHDkSixYtwoQJE2xeE0Jg+fLlWLBgASZMmICkpCSsXr0aOp0Oa9eubbEaqeUIAawtCcIDF6PxW5UKHVVV6KU12L239nEV1xc7+TzLUJcQPLaCiMiiScFn0KBB+Oc//4lLly7hD3/4A9avX482bdrAbDYjIyPDrZ2V8+fPIy8vDyNGjJCvabVaDBkyBHv37q336/R6PUpKSqw+yPNVmoEXrrTCy1dbwwAJw4Mq8HXbfLRTVdm9399WdAE1HR+AE5yJiCwcWtUVFBSE6dOnY8+ePTh+/Djmzp2LJUuWICYmBuPGjXN2jY2Sl5cHAIiNjbW6HhsbK79mz+LFixEeHi5/tGvXzqV1UvPlGJWYcDEan5cFQwGBFyOK8Y/YwgYnLfvb5oVAzXJ2gBOciYgsmrWcHQC6d++OpUuX4sKFC1i3bp0zamoWSbIexxBC2Fyrbf78+SguLpY/cnNzXV0iNcMOnRZjfovBSYMGkQoTPokvwIxWZTccvvK3zQuBmg0MAR5bQURk4fAGhnUplUrce++9uPfee531lk0SFxcHoLrzEx8fL1/Pz8+36QLVptVqodVq632dPINeAMuLwrDyWigAoI/WgPdiC5Ggatxmmf441KVU1B7q8p/AR0TUkGZ3fDxFp06dEBcXh4yMDPmawWBAZmYmUlJS3FgZNddJvQrjf4uRQ89DYWXYkHCl0aEH8M+hLkmSag4q5couIiIATuz4tISysjKcPXtW/vz8+fM4evQoIiIi0L59e8yePRvp6elITExEYmIi0tPTERQUhClTprixanJUlQD+cS0Ey4vCYISECIUJ6dHX8Lvgyia/lz8OdQGASqGA0WTiXj5ERNd5VfA5dOgQhg4dKn8+Z84cAMC0adPw8ccf48UXX0RFRQVmzJiBoqIiDBgwANu2bUNoaKi7SiYHnTcqMSc/Akf0GgDA8KAKpEddQ7TKsc6FPw51AddXdhm5qouIyMKrgk9qaiqEqP+/XCVJQlpaGtLS0lquKHIqswA+KQnG4sIwVAoFQiUzXo0qxv0hOjQwR/2G/HGoC6h1UClXdRERAfCy4EO+7WKVEi9eaYU913dgvj2wEkujr6FNE+by1MffzumysKzsYseHiKgagw+5nVkAn5UGYVFhOErNCgRIZrwUUYKHw8qdtsuyZajLX05mt5A7PpzjQ0QEgMGH3CxLr8afr4bjR331lgJ9tAYsiy5CZ439HZgd5a9DXSqu6iIissLgQ25RYpawrDAMa0qCYYaEYMmM51uX4pHwMqhccJZW7bO6/EnNUJd//dxERPVh8KEWJQTwr7JALCoMx1WTEgAwJliHP0UWI87BFVuNIQ91+VnHh0NdRETWGHyoxZwxqPDnq62wv7J6WKuz2ojXIotxR5Depd+3SgDlwk8nN18f6jJyqIuICACDD7WAcrOEt4tC8WFxCKogIUAy49lWpXi8VRm0LhjWqqvUXLNBeYifdXxUCnZ8iIhqY/AhlzEKYH1pMN4qCpWHtYYHVeCVyGK0Uzd/iXpjWYa5giQz1C0QtDyJfGQFl7MTEQFg8CEXEALYXB6AN4rCcd5Y/Vusg6oKr0QWY5gDx000l78eVwHUdHyM3MCQiAgAgw852f4KDRYXhuOn60dNRCpMmNW6FJPCyqFxU7elxE+PqwBqLWdnx4eICACDDznJzwYVlhaGY7uuetflIMmMJ8LL8ESrMoS4udPir3v4AFzVRURUF4MPNcsFoxLLi0LxRVkQBCSoIDA5rBzPtipFjAuXpzeFfw91cVUXEVFtDD7kkF+NSrx3LRQbS4NQhep/XEcH6zAvogSdWnDicmOUmPx3qIsdHyIiaww+1CRnDCq8ey0UX5cFwnw98NwRWIl5rUvQJ8Do5urssyxn97dzuoBa+/hwjg8REQAGH2qkLL0a714LxebyAIjrgeeuoErMbFWKfgEGN1fXMP8e6rre8eGqLiIiAAw+dAMn9GosLwrFf3SB8rV7girwbOtSJGk9s8NTlz+v6uI+PkRE1hh8qF7FJgkTfouGARIkCIwJrsAzrUvR3cknp7taqR+v6qoZ6mLHh4gIYPChBpSYFTBAggYCm9vmo4uXBR4LDnUBVVzVRUQEAFDc+BbydypJeG3oATjUBXBVFxGRBYMP+Tz/Huq6fmQFgw8REQAGH/ID/jzUpb6+gSGHuoiIqjH4kE8TouZ0dnZ8iIiIwYd8WrmQ5I0W/TH4KBVczk5EVBuDD/k0yzCXBgJaN50O707y5GZuYEhEBIDBh3xc7RVdkh8GH8tydh5ZQURUjcGHfFrNOV3+2fHgcnYiImsMPuTTSkz+u5QdqJnczFVdRETVGHzIp/nz5oUAoFLwyAoiotoYfMinlfrxHj4AoGbHh4jICoMP+bQSP97DB6g5pJRzfIiIqjH4kE+zLGf336EuS8eHwYeICGDwIR/HoS5uYEhEVBuDD/k0DnXxyAoiotoYfMin+ftQFw8pJSKyxuBDPs2fT2YHau3jw44PEREABh/ycf58MjtQs6rLyI4PEREABh+qhxDAeaPK3WU0G4e62PEhIqrN+/9lI6cQAjhrVGF/pRYHKjT4oVKLfJMSAKCVvPcfzRI/P6tL7vgw+BARAWDw8VtmAZw2qHCgUosDlVr8UKFBgVlpdY8GAn0CDHgorNxNVTZPpRkwCP8e6pKXs3Ooi4gIAIOP3ygzS/hJr8HhSg0OVWpwRK+R97ixCJDMuEVrwIBAAwYE6NFHa0CAFw+GWro9EgSCvbhr1RwqDnUREVlh8PFBQgC/VSlx+HrQOVypwSmDGmZIVvcFSWb0CzBgYIAeAwIN6KU1QCPV86ZeqFSe3yOg8KGfqylqhrrY8SEiAnw0+Lz33nt44403cOnSJfTs2RPLly/HnXfe6e6yXKbELOFYpQY/6dU4otfgJ70GV0xKm/vaqKrQT2tAv4Dqj5s0Rqh8OBD4++aFQO1DStnxISICfDD4bNiwAbNnz8Z7772H22+/Hf/4xz8wcuRInDx5Eu3bt3d3ec1mEMDPBjV+uj5c9ZNejV+Mapv7VBDoqTXiFq0B/QP0uCXAgHiVfwUAf1/RBQCq660uk1lACAFJ8uGkS0TUCD4XfJYtW4bHHnsMjz/+OABg+fLl2Lp1K1auXInFixe7ubqm0QvgtEGN43o1Tug1OK5X47RBDSNs//Fqp6pCb60BfbQG9AkwoqfGiEA/3bTPwt/P6QJqNjAEqld2aXy5xUdE1Ag+FXwMBgMOHz6Ml156yer6iBEjsHfvXrtfo9frodfr5c9LSkpcWmN9dGYJPxtUyDJocEJfHXb+a1Cjyk7ICVeYa4UcA3ppjYhS+m9Xoz4c6qpZ1QVUr+zScOsuIvJzPhV8rl69CpPJhNjYWKvrsbGxyMvLs/s1ixcvxsKFC1uiPADVE48vVylwyqDGSb0aJw1qnDKocd6ogrATclopzEjWGpCkNSJZU/2/7VQmcMTixkpqTW72V5ZVXQD38iEiAnws+FjUncfQ0NyG+fPnY86cOfLnJSUlaNeunVPrEUJgSWYOTph64NTVEBQK+489WmnCzRojkrXG6rCjMaINQ47Daoa62PEBgCqu7CIi8q3gExUVBaVSadPdyc/Pt+kCWWi1Wmi1WpfWJUkSMn4pwjm0AgSghEBndRV6aI24WWNED031/0b72eRjV+NQV/XvPaVCgsksuLKLiAg+Fnw0Gg369euHjIwM3HffffL1jIwMjB8/3o2VAU/dmgDztm3o0UqFbuFqr94Y0Fv4+8nsFqrrwYd7+RAR+VjwAYA5c+bgoYceQv/+/TFo0CC8//77yMnJwVNPPeXWuh5Mjgb+kw+oIwGF7fJzcj4OdVVTKxXQV5m5ezMREXww+Pz+979HQUEBXnvtNVy6dAlJSUn497//jQ4dOri7NGphJabrQ11+vuJNxfO6iIhkPhd8AGDGjBmYMWOGu8sgN+OqrmqWlV1c1UVEBG7qQb6rlJObAdQ6oZ3Bh4iIwYd8Vwnn+ACodVAph7qIiBh8yDdVASgXHOoCAPX1oS52fIiIGHzIR5WZa35r+/MhpUCtyc1czk5ExOBDvskyzBUkmaH2852v5cnN3MCQiIjBh3yTZWKzvw9zAbUnN7PjQ0TE4EM+qcTEic0WKiWXsxMRWTD4kE/iOV01VApuYEhEZMHgQz6plJsXytRKruoiIrJg8CGfxD18asj7+HCODxERgw/5Jnmoy8/P6QJqVnVVcVUXERGDD/kmDnXV4KouIqIaDD7kk6rAyc0WXNVFRFSDwYd8GoMPoOaqLiIiGYMP+TQOddWe3MxnQUTE4EM+jR0fQMlDSomIZAw+5NMYfGpNbuZQFxERgw/5tjAOddUcUsqODxERgw/5NnZ8uJydiKg2Bh/yaQw+NZObuYEhERGDD/kwNQS0krurcL+aoS6GQCIiBh/yWWEKMyQGn1pDXez4EBEx+JDPClPyH3qg1s7NXNVFRMTgQ74rlPN7AAAqBTs+REQWDD7kszixuZpaaTmdnc+DiIjBh3wWj6uopuIcHyIiGYMP+Sx2fKqpLUdWcDk7ERGDD/kuBp9qNYeU8nkQETH4kM/iUFc1y6ouDnURETH4kA9jx6eaWsFDSomILBh8yGcx+FST9/Fhx4eIiMGHfBeHuqrVnNXFIEhExOBDPosdn2ryqi52fIiIGHzIdzH4VOOqLiKiGgw+5LM41FVNPqSU+/gQETH4kO8KU7LDAQAqDnUREckYfMgnSRAIkfgPPcChLiKi2hh8yCeFKASub1/j92oOKWUQJCJi8CGfxInNNVQKdnyIiCwYfMgnMfjUUPPICiIiGYMP+SSu6KrBDQyJiGow+JBPYsenhmVVl9EkIAQDIRH5N68JPq+//jpSUlIQFBSEVq1a2b0nJycHY8eORXBwMKKiovDcc8/BYDC0bKHkEcLY8ZFZ9vEBABMnOBORn/Oa4GMwGPDAAw/g6aeftvu6yWTC6NGjUV5ejj179mD9+vX44osvMHfu3BaulDwBOz41LIeUAlzZRUSkcncBjbVw4UIAwMcff2z39W3btuHkyZPIzc1FQkICAODNN9/EI488gtdffx1hYWEtVSp5AAafGqpa6/pzCnUI0ijdWI17FVXoUVahhRoqmIz++xyI3KXcqEJZhRZVblxl6jXB50b27duHpKQkOfQAwD333AO9Xo/Dhw9j6NChdr9Or9dDr9fLn5eUlLi8VnI9Tm6uoa7V8Rnxt11urMRT9HN3AUR+7+HeesSFuud7+0zwycvLQ2xsrNW11q1bQ6PRIC8vr96vW7x4sdxNIu83PKgSJ/QaDA6qdHcpHkOpkDC2dwK2ZdX/58B/CFSZqwBIkMAdLolamoAAINz658+twSctLe2GoePgwYPo379/o95PkmwfpBDC7nWL+fPnY86cOfLnJSUlaNeuXaO+H3mee0MrcG9ohbvL8DjvTO7r7hI8wlXdVfxt398QGRSJEE2Iu8sh8jtlhjIU6ArQPjLFbTW4Nfg888wzmDRpUoP3dOzYsVHvFRcXhwMHDlhdKyoqgtFotOkE1abVaqHVahv1PYiIiMi7uTX4REVFISoqyinvNWjQILz++uu4dOkS4uPjAVRPeNZqtejXj2P6RERE5EVzfHJyclBYWIicnByYTCYcPXoUANC1a1eEhIRgxIgR6NGjBx566CG88cYbKCwsxLx58/DEE09wRRcREREB8KLg88orr2D16tXy5337Vs9Z2LFjB1JTU6FUKrFp0ybMmDEDt99+OwIDAzFlyhT89a9/dVfJRERE5GG8Jvh8/PHH9e7hY9G+fXt8++23LVMQEREReR2v2bmZiIiIqLkYfIiIiMhvMPgQERGR32DwISIiIr/B4ENERER+g8GHiIiI/AaDDxEREfkNBh8iIiLyGww+RERE5De8ZufmliKEAACUlJQ4941LSwG9HigsBCornfve5N10uurfG6WlgEbj7mp8WqmuFPpyPQoNhahU888hUUvTGXXQG/UoLSmFpsq5f99Z/t22/DteH0nc6A4/c+HCBbRr187dZRAREZEDcnNz0bZt23pfZ/Cpw2w24+LFiwgNDYUkSU5735KSErRr1w65ubk8Ld6F+JxbDp91y+Bzbhl8zi3Dlc9ZCIHS0lIkJCRAoah/Jg+HuupQKBQNJsXmCgsL4x+qFsDn3HL4rFsGn3PL4HNuGa56zuHh4Te8h5ObiYiIyG8w+BAREZHfYPBpIVqtFq+++iq0Wq27S/FpfM4th8+6ZfA5tww+55bhCc+Zk5uJiIjIb7DjQ0RERH6DwYeIiIj8BoMPERER+Q0GHyIiIvIbDD4t5L333kOnTp0QEBCAfv36Yffu3e4uyavt2rULY8eORUJCAiRJwldffWX1uhACaWlpSEhIQGBgIFJTU5GVleWeYr3Y4sWLceuttyI0NBQxMTG49957cfr0aat7+Kybb+XKlejVq5e8qdugQYOwefNm+XU+Y9dYvHgxJEnC7Nmz5Wt81s2XlpYGSZKsPuLi4uTX3f2MGXxawIYNGzB79mwsWLAAR44cwZ133omRI0ciJyfH3aV5rfLycvTu3RsrVqyw+/rSpUuxbNkyrFixAgcPHkRcXByGDx+O0tLSFq7Uu2VmZmLmzJnYv38/MjIyUFVVhREjRqC8vFy+h8+6+dq2bYslS5bg0KFDOHToEO666y6MHz9e/seAz9j5Dh48iPfffx+9evWyus5n7Rw9e/bEpUuX5I/jx4/Lr7n9GQtyudtuu0089dRTVtduuukm8dJLL7mpIt8CQHz55Zfy52azWcTFxYklS5bI1yorK0V4eLj4+9//7oYKfUd+fr4AIDIzM4UQfNau1Lp1a/HBBx/wGbtAaWmpSExMFBkZGWLIkCFi1qxZQgj+fnaWV199VfTu3dvua57wjNnxcTGDwYDDhw9jxIgRVtdHjBiBvXv3uqkq33b+/Hnk5eVZPXOtVoshQ4bwmTdTcXExACAiIgIAn7UrmEwmrF+/HuXl5Rg0aBCfsQvMnDkTo0ePxt133211nc/aec6cOYOEhAR06tQJkyZNwrlz5wB4xjPmIaUudvXqVZhMJsTGxlpdj42NRV5enpuq8m2W52rvmWdnZ7ujJJ8ghMCcOXNwxx13ICkpCQCftTMdP34cgwYNQmVlJUJCQvDll1+iR48e8j8GfMbOsX79evz44484ePCgzWv8/ewcAwYMwJo1a9CtWzdcvnwZixYtQkpKCrKysjziGTP4tBBJkqw+F0LYXCPn4jN3rmeeeQbHjh3Dnj17bF7js26+7t274+jRo7h27Rq++OILTJs2DZmZmfLrfMbNl5ubi1mzZmHbtm0ICAio9z4+6+YZOXKk/Ovk5GQMGjQIXbp0werVqzFw4EAA7n3GHOpysaioKCiVSpvuTn5+vk3iJeewrB7gM3eeZ599Fl9//TV27NiBtm3bytf5rJ1Ho9Gga9eu6N+/PxYvXozevXvjrbfe4jN2osOHDyM/Px/9+vWDSqWCSqVCZmYm3n77bahUKvl58lk7V3BwMJKTk3HmzBmP+P3M4ONiGo0G/fr1Q0ZGhtX1jIwMpKSkuKkq39apUyfExcVZPXODwYDMzEw+8yYSQuCZZ57Bxo0bsX37dnTq1MnqdT5r1xFCQK/X8xk70bBhw3D8+HEcPXpU/ujfvz+mTp2Ko0ePonPnznzWLqDX63Hq1CnEx8d7xu/nFplC7efWr18v1Gq1+PDDD8XJkyfF7NmzRXBwsPj111/dXZrXKi0tFUeOHBFHjhwRAMSyZcvEkSNHRHZ2thBCiCVLlojw8HCxceNGcfz4cTF58mQRHx8vSkpK3Fy5d3n66adFeHi42Llzp7h06ZL8odPp5Hv4rJtv/vz5YteuXeL8+fPi2LFj4uWXXxYKhUJs27ZNCMFn7Eq1V3UJwWftDHPnzhU7d+4U586dE/v37xdjxowRoaGh8r957n7GDD4t5N133xUdOnQQGo1G3HLLLfJyYHLMjh07BACbj2nTpgkhqpdMvvrqqyIuLk5otVoxePBgcfz4cfcW7YXsPWMAYtWqVfI9fNbNN336dPnvh+joaDFs2DA59AjBZ+xKdYMPn3Xz/f73vxfx8fFCrVaLhIQEMWHCBJGVlSW/7u5nLAkhRMv0loiIiIjci3N8iIiIyG8w+BAREZHfYPAhIiIiv8HgQ0RERH6DwYeIiIj8BoMPERER+Q0GHyIiIvIbDD5ERETkNxh8iMijpaWloU+fPm77/n/+85/x5JNPNureefPm4bnnnnNxRUTUHNy5mYjcRpKkBl+fNm0aVqxYAb1ej8jIyBaqqsbly5eRmJiIY8eOoWPHjje8Pz8/H126dMGxY8dsDnQlIs/A4ENEbpOXlyf/esOGDXjllVdw+vRp+VpgYCDCw8PdURoAID09HZmZmdi6dWujv+b+++9H165d8Ze//MWFlRGRozjURURuExcXJ3+Eh4dDkiSba3WHuh555BHce++9SE9PR2xsLFq1aoWFCxeiqqoKL7zwAiIiItC2bVt89NFHVt/rt99+w+9//3u0bt0akZGRGD9+PH799dcG61u/fj3GjRtnde3zzz9HcnIyAgMDERkZibvvvhvl5eXy6+PGjcO6deua/WyIyDUYfIjI62zfvh0XL17Erl27sGzZMqSlpWHMmDFo3bo1Dhw4gKeeegpPPfUUcnNzAQA6nQ5Dhw5FSEgIdu3ahT179iAkJAS/+93vYDAY7H6PoqIinDhxAv3795evXbp0CZMnT8b06dNx6tQp7Ny5ExMmTEDtxvltt92G3NxcZGdnu/YhEJFDGHyIyOtERETg7bffRvfu3TF9+nR0794dOp0OL7/8MhITEzF//nxoNBp8//33AKo7NwqFAh988AGSk5Nx8803Y9WqVcjJycHOnTvtfo/s7GwIIZCQkCBfu3TpEqqqqjBhwgR07NgRycnJmDFjBkJCQuR72rRpAwA37CYRkXuo3F0AEVFT9ezZEwpFzX+3xcbGIikpSf5cqVQiMjIS+fn5AIDDhw/j7NmzCA0NtXqfyspK/PLLL3a/R0VFBQAgICBAvta7d28MGzYMycnJuOeeezBixAhMnDgRrVu3lu8JDAwEUN1lIiLPw+BDRF5HrVZbfS5Jkt1rZrMZAGA2m9GvXz98+umnNu8VHR1t93tERUUBqB7ystyjVCqRkZGBvXv3Ytu2bXjnnXewYMECHDhwQF7FVVhY2OD7EpF7caiLiHzeLbfcgjNnziAmJgZdu3a1+qhv1ViXLl0QFhaGkydPWl2XJAm33347Fi5ciCNHjkCj0eDLL7+UXz9x4gTUajV69uzp0p+JiBzD4ENEPm/q1KmIiorC+PHjsXv3bpw/fx6ZmZmYNWsWLly4YPdrFAoF7r77buzZs0e+duDAAaSnp+PQoUPIycnBxo0bceXKFdx8883yPbt378add94pD3kRkWdh8CEinxcUFIRdu3ahffv2mDBhAm6++WZMnz4dFRUVCAsLq/frnnzySaxfv14eMgsLC8OuXbswatQodOvWDX/605/w5ptvYuTIkfLXrFu3Dk888YTLfyYicgw3MCQiqocQAgMHDsTs2bMxefLkG96/adMmvPDCCzh27BhUKk6hJPJE7PgQEdVDkiS8//77qKqqatT95eXlWLVqFUMPkQdjx4eIiIj8Bjs+RERE5DcYfIiIiMhvMPgQERGR32DwISIiIr/B4ENERER+g8GHiIiI/AaDDxEREfkNBh8iIiLyGww+RERE5Df+P1/5lAY6kEsRAAAAAElFTkSuQmCC", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], "source": [ "plt.plot(time, acel, label=\"Rocket Acceleration\")\n", "plt.axvspan(0, 15, color=\"RED\", alpha=0.4)\n", @@ -1009,7 +1051,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 44, "id": "742bbdab-a4e7-42cf-998d-22283bcdd99c", "metadata": {}, "outputs": [], @@ -1035,6 +1077,14 @@ " display.DisplayShape(shape)\n", " start_display()" ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "1be4b214-37ab-49aa-86e0-14e0897a324b", + "metadata": {}, + "outputs": [], + "source": [] } ], "metadata": { diff --git a/rocket_twin/systems/control/rocket_controller/rocket_controller.fmu b/rocket_twin/systems/control/rocket_controller/rocket_controller.fmu index 1f4dfaf..2d526fd 100644 Binary files a/rocket_twin/systems/control/rocket_controller/rocket_controller.fmu and b/rocket_twin/systems/control/rocket_controller/rocket_controller.fmu differ diff --git a/rocket_twin/systems/control/stage_controller/stage_controller.fmu b/rocket_twin/systems/control/stage_controller/stage_controller.fmu index 6cdf9f4..bdc15b6 100644 Binary files a/rocket_twin/systems/control/stage_controller/stage_controller.fmu and b/rocket_twin/systems/control/stage_controller/stage_controller.fmu differ diff --git a/rocket_twin/systems/control/station_controller/station_controller.fmu b/rocket_twin/systems/control/station_controller/station_controller.fmu index 6e558f4..437ec02 100644 Binary files a/rocket_twin/systems/control/station_controller/station_controller.fmu and b/rocket_twin/systems/control/station_controller/station_controller.fmu differ diff --git a/rocket_twin/utils/run_sequences.py b/rocket_twin/utils/run_sequences.py index ab9432b..c72b682 100644 --- a/rocket_twin/utils/run_sequences.py +++ b/rocket_twin/utils/run_sequences.py @@ -2,7 +2,7 @@ from cosapp.recorders import DataFrameRecorder -def run_sequences(sys, sequences): +def run_sequences(sys, sequences, includes = []): """Run the command sequences over a system. Inputs @@ -13,7 +13,7 @@ def run_sequences(sys, sequences): the commands to be applied """ rk = sys.add_driver(RungeKutta("rk")) - rk.add_recorder(DataFrameRecorder(includes=["rocket.a"], hold=True)) + rk.add_recorder(DataFrameRecorder(includes=includes, hold=True)) for seq in sequences: print("sequence ", seq["name"])