{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Hinge Torques Monitoring\n",
"\n",
"**Note:** The settings in this example are by no means a validation setup; they are crafted to showcase the capabilities of Flow360 and we have intentionally reduced node count and example FC cost. For rigorous validation, modify the settings as needed."
],
"id": "212cac3f"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
""
],
"id": "db0eec14"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"This tutorial will show how to set up a case where moments on specific surfaces are monitored throughout the simulation. The content will cover:\n",
"- creating a project from a geometry,\n",
"- meshing using the beta mesher and automated farfield,\n",
"- setup of a steady case,\n",
"- setup of `SurfaceIntegralOutput` using `UserVariable` to capture hinge moments on control surfaces,\n",
"- results postprocessing."
],
"id": "ae363631"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 1. Imports and Project Creation\n",
"\n",
"To begin working with the **Flow360** API it is necessary to firstly import the `flow360` python package that will be the interface between you and the software. The files necessary for this example are being downloaded through the `TutorialUDDForcesMoments` example. \n",
"\n",
"All simulations are managed through the `Project` object, which stores all the runs associated with a certain root asset. Root asset can be a geometry, surface mesh or volume mesh. In case of this tutorial the `Project` will be started from the geometry."
],
"id": "876f7dc7"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"import flow360 as fl\n",
"from flow360.examples import TutorialUDDForcesMoments"
],
"execution_count": 1,
"outputs": [],
"id": "b8542815"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"TutorialUDDForcesMoments.get_files()\n",
"\n",
"project = fl.Project.from_geometry(\n",
" TutorialUDDForcesMoments.geometry, name=\"Monitoring hinge torques\"\n",
")"
],
"execution_count": 2,
"outputs": [
{
"output_type": "display_data",
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "ab48cdf3467f44208ef23db1f780ff91",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"Output()"
]
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/html": [
"\n"
],
"text/plain": []
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/html": [
"
[10:16:16] INFO: Geometry successfully submitted: \n",
" type = Geometry \n",
" name = Monitoring hinge torques \n",
" id = geo-ca5ad8b9-6320-4438-ba13-9b22704743f3 \n",
" status = uploaded \n",
" project id = prj-c1642540-6862-4460-89f3-220e400aa909 \n",
" \n",
"
\n"
],
"text/plain": [
"\u001b[2;36m[10:16:16]\u001b[0m\u001b[2;36m \u001b[0m\u001b[36mINFO\u001b[0m: Geometry successfully submitted: \n",
"\u001b[2;36m \u001b[0m type = Geometry \n",
"\u001b[2;36m \u001b[0m name = Monitoring hinge torques \n",
"\u001b[2;36m \u001b[0m id = geo-\u001b[93mca5ad8b9-6320-4438-ba13-9b22704743f3\u001b[0m \n",
"\u001b[2;36m \u001b[0m status = uploaded \n",
"\u001b[2;36m \u001b[0m project id = prj-\u001b[93mc1642540-6862-4460-89f3-220e400aa909\u001b[0m \n",
"\u001b[2;36m \u001b[0m \n"
]
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/html": [
"
INFO: Waiting for geometry to be processed. \n",
"
\n"
],
"text/plain": [
"\u001b[2;36m \u001b[0m\u001b[2;36m \u001b[0m\u001b[36mINFO\u001b[0m: Waiting for geometry to be processed. \n"
]
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "cdcad31d5909419090ba954685f233dc",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"Output()"
]
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/html": [
"\n"
],
"text/plain": []
},
"metadata": {}
}
],
"id": "c24464c9"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"If the `Project` has already been created before it is possible to retrieve it using its unique ID."
],
"id": "61daeb53"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"# project=fl.Project.from_cloud(\"prj-e8aa4e7a-a6b8-47f5-97c4-7e64bbf76aeb\")"
],
"execution_count": 3,
"outputs": [],
"id": "6909cc0e"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"`Project` contains all the information about its root asset, in this case - geometry. To be able to correctly reference it it is necessary to group the faces and edges by the attribute that was assigned to distinguish the faces in a CAD file. "
],
"id": "b630c2fc"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"geometry = project.geometry\n",
"\n",
"geometry.group_edges_by_tag(\"edgeName\")\n",
"geometry.group_faces_by_tag(\"groupName\")\n",
"\n",
"print(\"Geometry boundaries:\")\n",
"for boundary in geometry.entity_info.get_boundaries():\n",
" print(f\"- {boundary.name}\")"
],
"execution_count": 4,
"outputs": [
{
"output_type": "display_data",
"data": {
"text/html": [
"
[10:17:08] INFO: Regrouping face entities under `groupName` tag (previous `faceName`). \n",
"
\n"
],
"text/plain": [
"\u001b[2;36m[10:17:08]\u001b[0m\u001b[2;36m \u001b[0m\u001b[36mINFO\u001b[0m: Regrouping face entities under `groupName` tag \u001b[1m(\u001b[0mprevious `faceName`\u001b[1m)\u001b[0m. \n"
]
},
"metadata": {}
},
{
"output_type": "stream",
"text": [
"Geometry boundaries:\n",
"- vTailLeft\n",
"- fuselage\n",
"- wingLeft\n",
"- wingRight\n",
"- vTailRight\n",
"- rudderLeft\n",
"- rudderRight\n",
"- aileronLeft\n",
"- aileronRight\n"
],
"name": "stdout"
}
],
"id": "ff6f95dd"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 2. Meshing\n",
"When the projects starts from the geometry it is necessary to set up the automated meshing workflow.\n",
"\n",
"Firstly the AutomatedFarfield object is created to indicate that the geometry only contains the simulated object and the farfield has to be created around it.\n",
"\n",
"Then the defaults for meshing are defined. Those are the \"global\" settings for the mesher."
],
"id": "b1ba1819"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"farfield = fl.AutomatedFarfield()\n",
"\n",
"with fl.SI_unit_system:\n",
" meshing = fl.MeshingParams(\n",
" defaults=fl.MeshingDefaults(\n",
" surface_max_edge_length=0.5,\n",
" curvature_resolution_angle=10 * fl.u.deg,\n",
" boundary_layer_first_layer_thickness=2e-6,\n",
" boundary_layer_growth_rate=1.2,\n",
" ),\n",
" volume_zones=[farfield],\n",
" )"
],
"execution_count": 5,
"outputs": [
{
"output_type": "display_data",
"data": {
"text/html": [
"
[10:17:15] INFO: using: SI unit system for unit inference. \n",
"
\n"
],
"text/plain": [
"\u001b[2;36m[10:17:15]\u001b[0m\u001b[2;36m \u001b[0m\u001b[36mINFO\u001b[0m: using: SI unit system for unit inference. \n"
]
},
"metadata": {}
}
],
"id": "13df1c5b"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"To refine specific parts of the mesh the `refinements` field is used, where user can define surface (such as `SurfaceRefinement` and `SurfaceEdgeRefinement`) and volumetric refinements (`UniformRefinement`).\n",
"\n",
"It is worth to notice how the specific parts of the geometry are referenced using the `geometry` object. They can be then referenced using face names directly or through a wildcard entry (`*`). "
],
"id": "95e7d8b7"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"with fl.SI_unit_system:\n",
" box1 = fl.Box(name=\"box1\", size=[2, 6, 3], center=[6.5, 9, 0], axis_of_rotation=[0, 1, 0])\n",
" box2 = fl.Box(name=\"box2\", size=[2, 6, 3], center=[6.5, -9, 0], axis_of_rotation=[0, 1, 0])\n",
" box3 = fl.Box(name=\"box3\", size=[4, 8, 3], center=[12, 0, 2], axis_of_rotation=[0, 1, 0])\n",
"\n",
" meshing.refinements.extend(\n",
" [\n",
" fl.SurfaceEdgeRefinement(\n",
" method=fl.AngleBasedRefinement(value=1 * fl.u.deg),\n",
" edges=[geometry[\"leadingEdge\"]],\n",
" ),\n",
" fl.SurfaceEdgeRefinement(\n",
" method=fl.HeightBasedRefinement(value=5e-3),\n",
" edges=[geometry[\"trailingEdge\"]],\n",
" ),\n",
" fl.SurfaceRefinement(max_edge_length=0.5, faces=[geometry[\"wing*\"]]),\n",
" fl.UniformRefinement(name=\"box_refinement1\", entities=[box1, box2, box3], spacing=0.2),\n",
" ]\n",
" )"
],
"execution_count": 6,
"outputs": [
{
"output_type": "display_data",
"data": {
"text/html": [
"
[10:17:18] INFO: using: SI unit system for unit inference. \n",
"
\n"
],
"text/plain": [
"\u001b[2;36m[10:17:18]\u001b[0m\u001b[2;36m \u001b[0m\u001b[36mINFO\u001b[0m: using: SI unit system for unit inference. \n"
]
},
"metadata": {}
}
],
"id": "34a45fee"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 3. Boundary Conditions, Flow physics and numerics settings\n",
"\n",
"To describe the physical behaviour of the flow the Boundary Conditions or surface models are used. In this case the ones needed would be the `Freestream` condition that is assigned to the farfield and the `Wall` condition for the object. The `Fluid` model is a volume model that enables access to the controls of the fluid solver.\n",
"\n",
"The `operating_condition` is responsible for defining the fluid behaviour at the freestream and is used to set the velocity, angle of attack and the thermal state of the fluid in this case.\n",
"\n",
"`reference_geometry` defines values that are used as a reference when calculating various coefficients. \n",
"\n",
"Finally `time_stepping` describes whether the simulation is steady or transient and how many time (or pseudo time) steps should be performed."
],
"id": "397fe4d3"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"with fl.SI_unit_system:\n",
" reference_geometry = fl.ReferenceGeometry(\n",
" area=60, moment_center=[5.7542, 0, 0], moment_length=[1, 1, 1]\n",
" )\n",
" operating_condition = fl.AerospaceCondition(\n",
" velocity_magnitude=50,\n",
" alpha=10 * fl.u.deg,\n",
" atmosphere=fl.ThermalState(temperature=288.15),\n",
" )\n",
" models = [\n",
" fl.Fluid(\n",
" navier_stokes_solver=fl.NavierStokesSolver(\n",
" absolute_tolerance=1e-9,\n",
" linear_solver=fl.LinearSolver(max_iterations=35),\n",
" ),\n",
" turbulence_model_solver=fl.SpalartAllmaras(\n",
" linear_solver=fl.LinearSolver(max_iterations=25)\n",
" ),\n",
" ),\n",
" fl.Wall(\n",
" surfaces=[geometry[\"*Left\"], geometry[\"*Right\"], geometry[\"fuselage\"]],\n",
" ),\n",
" fl.Freestream(surfaces=[farfield.farfield]),\n",
" ]\n",
" time_stepping = fl.Steady(max_steps=5000)"
],
"execution_count": 7,
"outputs": [
{
"output_type": "display_data",
"data": {
"text/html": [
"
[10:17:20] INFO: using: SI unit system for unit inference. \n",
"
\n"
],
"text/plain": [
"\u001b[2;36m[10:17:20]\u001b[0m\u001b[2;36m \u001b[0m\u001b[36mINFO\u001b[0m: using: SI unit system for unit inference. \n"
]
},
"metadata": {}
}
],
"id": "60efc886"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 4. Definition of `UserVariable`s\n",
"\n",
"In the next cell, we define several custom `UserVariable`s to compute hinge torques for the control surfaces (right/left ailerons and right/left rudders). \n",
"Each variable calculates the local contribution to the hinge torque by taking the moment (cross product) of the force per unit area about a specified torque center, and then projecting that moment onto the relevant hinge axis. The resulting user variables represent the hinge torque per unit area at each surface node, \n",
"with units of pressure times length (Pa·m), and can be integrated over the relevant surfaces to obtain the total hinge torque values.\n",
"\n",
">Note: The variables are defined **per unit area** because `SurfaceIntegralOutput` integrates over the surface area, not over individual nodes."
],
"id": "2b9a8efd"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"right_aileron_torque_center = [5.7542, 7, 0] * fl.u.m\n",
"right_aileron_torque_axis = [0, 1, 0]\n",
"\n",
"right_aileron_hinge_torque_partial_dim = fl.UserVariable(\n",
" name=\"right_aileron_hinge_torque_partial\",\n",
" value=fl.math.dot(\n",
" fl.math.cross(\n",
" fl.math.subtract(fl.solution.coordinate, right_aileron_torque_center),\n",
" fl.solution.node_forces_per_unit_area,\n",
" ),\n",
" right_aileron_torque_axis,\n",
" ),\n",
").in_units(new_unit=fl.u.Pa * fl.u.m)\n",
"\n",
"left_aileron_torque_center = [5.7542, -7, 0] * fl.u.m\n",
"left_aileron_torque_axis = [0, -1, 0]\n",
"\n",
"left_aileron_hinge_torque_partial_dim = fl.UserVariable(\n",
" name=\"left_aileron_hinge_torque_partial\",\n",
" value=fl.math.dot(\n",
" fl.math.cross(\n",
" fl.math.subtract(fl.solution.coordinate, left_aileron_torque_center),\n",
" fl.solution.node_forces_per_unit_area,\n",
" ),\n",
" left_aileron_torque_axis,\n",
" ),\n",
").in_units(new_unit=fl.u.Pa * fl.u.m)\n",
"\n",
"right_rudder_torque_center = [12.01, 0.861, 0.861] * fl.u.m\n",
"right_rudder_torque_axis = [0, fl.math.sqrt(2) / 2, fl.math.sqrt(2) / 2]\n",
"\n",
"right_rudder_hinge_torque_partial_dim = fl.UserVariable(\n",
" name=\"right_rudder_hinge_torque_partial\",\n",
" value=fl.math.dot(\n",
" fl.math.cross(\n",
" fl.math.subtract(fl.solution.coordinate, right_rudder_torque_center),\n",
" fl.solution.node_forces_per_unit_area,\n",
" ),\n",
" right_rudder_torque_axis,\n",
" ),\n",
").in_units(new_unit=fl.u.Pa * fl.u.m)\n",
"\n",
"left_rudder_torque_center = [12.01, -0.861, 0.861] * fl.u.m\n",
"left_rudder_torque_axis = [0, -fl.math.sqrt(2) / 2, fl.math.sqrt(2) / 2]\n",
"\n",
"left_rudder_hinge_torque_partial_dim = fl.UserVariable(\n",
" name=\"left_rudder_hinge_torque_partial\",\n",
" value=fl.math.dot(\n",
" fl.math.cross(\n",
" fl.math.subtract(fl.solution.coordinate, left_rudder_torque_center),\n",
" fl.solution.node_forces_per_unit_area,\n",
" ),\n",
" left_rudder_torque_axis,\n",
" ),\n",
").in_units(new_unit=fl.u.Pa * fl.u.m)"
],
"execution_count": 8,
"outputs": [],
"id": "92e728f8"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Below, each hinge torque user variable (representing the torque per unit area contribution) is integrated over the corresponding control surface using `SurfaceIntegralOutput`. These integrals compute the total hinge torque for each control surface (right/left aileron and rudder) by summing their local torque contributions defined above.\n",
"\n",
"The `output_fields` for each surface integral specify which user variable is integrated, and the `surfaces` parameter specifies which part of the geometry the integration is applied to."
],
"id": "e68a0299"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"outputs = [\n",
" fl.SurfaceIntegralOutput(\n",
" name=\"right_aileron_hinge_torque\",\n",
" output_fields=[right_aileron_hinge_torque_partial_dim],\n",
" surfaces=geometry[\"aileronRight\"],\n",
" ),\n",
" fl.SurfaceIntegralOutput(\n",
" name=\"left_aileron_hinge_torque\",\n",
" output_fields=[left_aileron_hinge_torque_partial_dim],\n",
" surfaces=geometry[\"aileronLeft\"],\n",
" ),\n",
" fl.SurfaceIntegralOutput(\n",
" name=\"right_rudder_hinge_torque\",\n",
" output_fields=[right_rudder_hinge_torque_partial_dim],\n",
" surfaces=geometry[\"rudderRight\"],\n",
" ),\n",
" fl.SurfaceIntegralOutput(\n",
" name=\"left_rudder_hinge_torque\",\n",
" output_fields=[left_rudder_hinge_torque_partial_dim],\n",
" surfaces=geometry[\"rudderLeft\"],\n",
" ),\n",
"]"
],
"execution_count": 9,
"outputs": [],
"id": "b50144da"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 5. `SimulationParams`\n",
"\n",
"Finally all of the previously defined fields are input to the `SimulationParams` object that defines a whole process from meshing to performing the simulation.\n",
"\n",
"After preparing those parameters the simulation can be kicked off using the `run_case` command."
],
"id": "e2f05891"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"with fl.SI_unit_system:\n",
" params = fl.SimulationParams(\n",
" meshing=meshing,\n",
" operating_condition=operating_condition,\n",
" reference_geometry=reference_geometry,\n",
" models=models,\n",
" outputs=outputs,\n",
" )\n",
"\n",
"case = project.run_case(params, name=\"Hinge torques monitoring\")"
],
"execution_count": 10,
"outputs": [
{
"output_type": "display_data",
"data": {
"text/html": [
"
[10:17:26] INFO: using: SI unit system for unit inference. \n",
"
\n"
],
"text/plain": [
"\u001b[2;36m[10:17:26]\u001b[0m\u001b[2;36m \u001b[0m\u001b[36mINFO\u001b[0m: using: SI unit system for unit inference. \n"
]
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/html": [
"
[10:17:27] INFO: using: SI unit system for unit inference. \n",
"
\n"
],
"text/plain": [
"\u001b[2;36m[10:17:27]\u001b[0m\u001b[2;36m \u001b[0m\u001b[36mINFO\u001b[0m: using: SI unit system for unit inference. \n"
]
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/html": [
"
[10:35:38] INFO: Saved to C:\\Users\\user\\AppData\\Local\\Temp\\tmprp_0gk5h\\5874a8b1-0f4e-402b-99bd-327224fd4adb.csv \n",
"
\n"
],
"text/plain": [
"\u001b[2;36m[10:35:38]\u001b[0m\u001b[2;36m \u001b[0m\u001b[36mINFO\u001b[0m: Saved to C:\\Users\\user\\AppData\\Local\\Temp\\tmprp_0gk5h\\\u001b[93m5874a8b1-0f4e-402b-99bd-327224fd4adb\u001b[0m.csv \n"
]
},
"metadata": {}
},
{
"output_type": "execute_result",
"data": {
"text/plain": [
""
]
},
"metadata": {},
"execution_count": null
},
{
"output_type": "display_data",
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjwAAAGxCAYAAABmyWwBAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjcsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvTLEjVAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAZT5JREFUeJzt3QeYU1XaB/B3emWGoQ5IFwVEKYIiFhaVIriuqOvaPkWXxRVhLSgqFpoFBUVZRV0LomtDd1fWFUQQRVEQFCmCgIAgojRpw/SW7/mf5Ny5ySQzyUzKzc3/5xOHSW6Se3IzuW/e855z4hwOh0OIiIiIbCw+0jtAREREFGoMeIiIiMj2GPAQERGR7THgISIiIttjwENERES2x4CHiIiIbI8BDxEREdkeAx4iIiKyvcRI74AVVFZWyq+//ioNGjSQuLi4SO8OERER+QFzJx87dkxatmwp8fE153AY8IioYKd169b+vLZERERkMT///LO0atWqxm0Y8IiozI5+wbKysoJ6EMrKymTRokUyaNAgSUpKEruxe/tioY12b18stJHti348hnWTl5enEhb6PF4TBjwiRjcWgp1QBDzp6enqce36QWvn9sVCG+3evlhoI9sX/XgM68efchQWLRMREZHtMeAhIiIi2wtpwDN16lQ57bTTVN9as2bNZNiwYbJlyxa3bYqLi2X06NHSuHFjyczMlMsuu0z27dvnts2uXbvkwgsvVClpPM64ceOkvLzcbZulS5fKqaeeKikpKdKxY0eZM2dOKJtGREREUSSkNTyfffaZCmYQ9CBAuffee1XR4Pfffy8ZGRlqm9tvv13mz58v7777rmRnZ8uYMWPk0ksvlS+//FLdXlFRoYKd3NxcWb58uezZs0euu+461Q//yCOPqG127NihtrnpppvkjTfekCVLlshf/vIXadGihQwePDiUTSQKCN7P6Ku3EuxPYmKi+vKB/bMju7eR7Yt+PIa+4XyfkJAg9RXnwCD2MDlw4IDK0CAQ6tevnxw9elSaNm0qb775pvzxj39U22zevFm6dOkiK1askDPOOEM+/PBD+f3vf6+Gjjdv3lxt8/zzz8vdd9+tHi85OVn9G0HThg0bjOe68sor5ciRI7Jw4UK/qrwRbGF/QlG0vGDBAhk6dKhtiyXt3L5gtRF/Znv37lXvSavBvhUVFUlaWppt56GyexvZvujHY1izhg0bqsSH599vIOfvsI7Swg5Bo0aN1M/Vq1erk8mAAQOMbTp37ixt2rQxAh78POWUU4xgB5C1GTVqlGzcuFF69uyptjE/ht7mtttuC1vbiGqigx0E/OiatdJJFxNv5ufnqy7l2ibuilZ2byPbF/14DH0HgoWFhbJ//371O3pu6ioxnAcTAchZZ50lJ598snESQIYGkZsZghvcprcxBzv6dn1bTdsg8tPf6sxKSkrURcN2gOAr2N0N+vGs1o0RLHZvXzDaiC6Uw4cPq2xmTk6OWA0+UEpLS1X9m5UCsWCyexvZvujHY+gb/m4RQ6BXB5+h5u6tQD6XwxbwoJYHXU5ffPGFRBqKqSdPnlztekxMhm/fobB48WKxM7u3rz5tRO0IUrH4g9XBtRVhena7s3sb2b7ox2PoHT4/kcBAja550BKyP5YKeFCI/MEHH8jnn3/uNvUzTgL41oVUvznLg1FauE1vs2rVKrfH06O4zNt4juzC7+jP88zuwPjx42Xs2LHVZmpEQXUoanhwohw4cKAta1zs3r5gtBGFspjFG6MVU1NTxapr0dh5LTm7t5Hti348hrV/juJ8jvpf8+doIF8iE0N9AP/2t7/Je++9p4aNt2/f3u32Xr16qRMIIjYMRwcMW8cw9L59+6rf8fPhhx9W/XeofwCcfBCYnHTSScY2KCo1wzb6Mbylx3DxhH0J1Uk7lI9tBXZvX33aiC4tnGRRO2LF+hF8cwK9j3Zk9zayfdGPx7Bm+LvF36/n53Agn8nxoe7Gev3119UoLHyzQq0NLkhLASqrR4wYobItn376qSpivuGGG1SggoJlQNYFgc21114r69atk48++kjuv/9+9dg6aMFw9B9//FHuuusuNcrr2WeflXfeeUcNeSeiuuvfv39Axf/z5s1T82Chj90Kgwauv/56Nf9XMNsYqv1o166dPPXUUyHdDwq/SZMmSY8ePQK6D07s+Fuqzc6dO9W2a9eulWi2dOlS1Y5Qj2INacDz3HPPqZFZ+EBBZbW+zJ0719jmySefVMPOkeFBqgrdU//5z3+M2/HBie4w/EQg9H//939qHp4pU6YY2yBzhGHpyOp0795dnnjiCXnppZc4Bw9RmP31r39VU0ygC+/BBx/060RPTl9//bXceOONlno5GIRJvQOVO++8U/VihAJKMTA3nR4IFKoAzC5C3qVVG/TFzZo1S118adu2bbUuK08IqtasWSNWVVJWIbuPlkr7Js4JF4nsBsO+0fWMKSFatmwZ1BoqK3eXmrss6wOj+OwKtZoYkRuL7cNUCLiEAhIBupbVikotdtzt15ltUfe8t1HOfXypfLfbORcRUbTBVA74tnrcccepmdL79OmjUtGAn+i2hvPOO08FAPgS8uqrr8p///tf9TsuevvaUvTIAv/ud79TX4gwe7q3b6Xo/kEGwhx4oHscAyCwVA26uD2/dBUUFKgMMU5AyDYjGxxIOwHL1uA53n//fdXdjq511B364/HHH1fPi/1Dt7x5SK1nNgWvAzLVl1xyiRo9esIJJ6jnNMPvuB63X3TRRer19uwawMjYc845RxV8IiNwyy23qNehNjh+P/30kyoN0MdP+/e//y1du3ZVbcd+e76OuA4ZPrzWqLfUmSu8dphnDfuLduF+5gEr3jKC6G7Ee8pc64KRtsjso03I6v/rX/+SQLpO0CPQrVs39f5C+YR50tqDBw/KVVddpY4/9hPzwL311lvVXhsMxsG+NWnSRAX5+r2IduE59O+e711k8jAAAvdDWQfe599++63UhWeXlm4fMkq9e/dW+3/mmWcaSzrNmTNHjVBGeYg+pnoZJrxnsEIBAm8cM7zm2M7soYceUrW0+FvHtvfcc49b2/TxQ90tvvR06tRJXf/Pf/5TTj/9dPX+w/VXX321Ma9OODHgCZNdh51D5zbtse6wZAofNZlWaXnYL/WZWB0f8Jjk8+2335b169fL5ZdfLhdccIFs3brV7UMVJ0Ok2XEy/tOf/qS2we+4YDt/4IP01ltvlU2bNvndNY2TJz68Z8+erU7yhw4dUgMmzBAEYaZ3BGGYhgInCM+TTU3tNA+Ffeyxx1RAgglQ9YCKmqBOcfv27eonAhPsa21r/uHkhNcQ+4GZvq+55hrVLr2kDroPcYJBdhsnmwceeMDt/ng+7DtKBvAYCCTx2qCNtUFpAUbVonxAHz9ArSX2CbPZf/fdd+qEjuf1bAuCOwQj2DfcvnLlSlWziefGCfrcc89VJ9BAIdh57bXX1Iz7eO0RkKHUAcfVX1iPEe8XBB84wSNY1MEnRgNhQI2evR/BGmpIPUcL4xgie4FlkLAveCx45ZVX1Gulf/eE0YLDhw9Xx+Grr75SASuObTCHo993332qfd98842aFuPPf/6zuv6KK66QO+64QwWr+pjiOsD7HEEIVjfAMcbalOeff77xfsMXDwQyeN/jdgSuKFvxhGALnwUoMUE5CuC1xXt52bJl6n2FQA3v13AL60zLsay8wnmiOZBfNeEhxa6isgo5acJHYX/e76cMlvTkwP/skcHABzl+6u4qZEGwdAuux7p2+qSPmdR1mh3fwJExCTTtjm/OWFMvEMiOYMoJfT+chDDIwdzlhmAIAynwQa5PWuapMvxpp/4Ax+AInND9hQnTnnnmGdUNgRnlsf4fTg4jR470eR+cFJBtADz33//+d3XiRRDzj3/8Q32Dnj59usp6IHOEwRt6H3VwgCBJF2Xj5IrHQFYBJ6uapknAccS+4tu8+fjNmDFDvX46uDrxxBPV+ojYD/NJDBkCnFw1bI/9RtCp74f1Ef1Z/kfDewnt+/jjj41RuB06dFDBA14PtMsfEydOVFkW83sAwTECOWR2cMw1jDTG+wgDYZCl0PBaTps2zecSCL6Ys1XwwgsvqPvoJZeCAYGJfi3w5QHvtWLXsG5kN/XcYBpeP7yvEPDowUAIWFGPhOwZgr6nn35aBawYWAQTJkxQXxrwd2WGrCi+CJi7shBw6XnIkD3CexBrbOrZz8OFAU+YA56D+aXhekqioME3eXQZ4STleQJC90ywIR0fCAyOwLdVdD9p+FDH4+isFjIiqCkwb4OTuk67B9JOfJijSyQQ+FZtniEWAQqerybm58CJBCcL3RWAb9E4aZh5/o4uCWR28O1cw+uBkw9eD6xbGChk3S6++GK36zCDPgJOvHa6jZ7HEPdDd48ZgpZAAp5t27ap7JoOVjQcVywz5C/zlCX6PYD9A7QBQRUCnF9++UU9No6/56S0yALVBeaIw0hjZBdxLPF8aBMK/YPF/L7RSzHs379fZWW8wfsEwYfn3zJGVCNLqN9vN998s9vtCAA/+eQTt+vQBehZt4OMEIJMZPbwt6qH4OOLhZ5eJhwY8IRJuesA/8YMDyHzkZSgsi2ReN66wIchTmT44PJctTgU39BwcjdDQbBnd1woljPxt511WYTUs/Aa99cf/MG8j2d7MHIOdTuefJ38QnUM/VHbcdbZBHQ3IRNj5m1utbpApmrmzJkqgMPJG+1AhgyBT33bB+jOQp0QngMDcrDfCMA8H78+zO8b/T6trOF9g9cVgZG3GjvPpZ9q4/m6oF4M3dKYYgbZLNQ27d69W10XzDb7gwFPmJS5MjwMeEh/CNWlaylS8O0Z30TxLREFsP7CNz3cr75QZ4E5vHAy1B/g5rlHUPyJD2zUiehuAUw/r2sRAEWuOBFgG32yxxpnP/zwg5H+r2s7IwFZCc/Rq6jZMEPb0d2EuZHqwtvxQ1YIdStm+B1ZMc8g0fN+eO3NUMPieZzNBcT6OOsTuLlI3N/uK2/wvJ7vAZ3tQluQwUJdkA4UcLs/mQjsZ23vdzw+ukNRtwPI7Pz2228SLslejineJ/j7QlbUPBDA8/2GuiQUomu+6pTMMDceAjx0r+LvFFnKuhZp1xeLlsOkopJdWhS9cDJDLQg+7FB0iO4Q9PnjQwzftn3Bhye6VJAOx4d6XbMyGBWDhQNRM4EUO6axQHGlGYqcH330UVV3gA9ZpN/No5WQoUEtAQpWkYbHiRU1J+bh5HVtZyQgc4N23n333eqEjBoU1KOADgpxG+pkdKEwCq9RsO1P0bI+flgSCF07+qSMuhzUHmEUFp4Xz4naJHPdizfIMqH7CrUh2A/cx7M7C/UtCNpQlIxt0A1iDoBQT4TnQaEynhfvBZw8UV+i2+4PFGKjDfo9gBFTenQYanNQcIvXDd1ceJ09ly6q6fXC4yJ4QCDlDR4fo5bw2AgA8X7ztgRSqLRr1069r/F+wDFFd92AAQNUlgmvAepyUFSM9qP4WQfRqGV6+eWX1euMY4OCc/xt15bpRGCJIAvHG4+LwQx470QCA54wKXcFPMzwULRC0S4CAZzw8G0PH474hldT1wgKcrEt6jnw7d0zM+AvfPvGt2IEOigURhDieYLFfmE0DboM8OGNk6NnzQgCJmRuMCoHH/Jnn312tVqMurQzEpCxQkEpAjMMDUZBNoq2zd07qOVAMSwCE7QbGSwUm/o7TxICA5ykjj/+eGOeIGQDUN+CUWyY8A6Ph+1qG3WD4d8vvvii6srBMcSJFbUsZujmQHEzCptRj4SRS+aMAuBkiW0QhOJ9gUJoBKOeSxfVBIExAmQcewQn//vf/4y6E+wT2oh9QaCN4l5/J8/EyCgESxh+7aumCEEDgiE8B96vCAT9GeUXLJdddpl6zTBKDscUQ+4RtCBbiOwoipIR+GMUHqYlaN68ubofAjO8v/B3h31H0IRjXtv6gHgOjODDexXvAfwNIuiNhDhHfcap2gQqx5FqQzFVKBYPxRtp0ro0OVxYJgiGtz08VBLi7bGAoW4f0rNWnhwukm3E6Ah8OOAD2YqLh5pHT9hxnalYaKNuH75FY7RSMAtgQwknQtTH1LakQLCOH2pUcKJHwBFobUqoReN7dODAgSogRMYq1O3z9TkayPk7eooIbJLhQXh5qKBUmjYIToEdEcUuZL2QCcGQdwzVxjdnf7uriAKBkWSY6gGZL9RqITOE9xwyWtEiOsJIGyivqKqQZ7cWxSoM99VT7XtehgwZItHKV5twwWRroYJaChTYomsJwQ5mmsZEgP7AftW039EIC0n7ag9uizZW+nuJM3V7oSsQ3YCYZBRdw9GCGZ4wZ3iAc/FQrMJJB5O7eRPOws1gq2m1as/h08GExZdxqUt3AeqqIr3KNmpAgjnjLmqJfBVP47VBrUw0VXFY6e8lLS1NZXSiGQOeCAQ8zPBQrMIkb7jYTV2HfUcSTmDRuN81QUATzgLgULPr30uksEsrDBDrmL9UMOCJTdH0zZKIyG6fnwx4wsA156DhNy4vEVP0yC4U/RERUeD052d9RgOzSysMTL1ZCjM8sQUjGjAEVq+BhDV5Al2WIJRQ/4Ep3jHsM1qGwwbK7m1k+6Ifj6HvzA6CHXx+4nO0ptm8a8OAJwIZnoNcTyvm6JWJddBjJfhAwSKBdVkfKlrYvY1sX/TjMaxZbavQ+4MBTxiwS4twksVaTyioDMWil/WB/cHyARhuaufJI+3cRrYv+vEY+oa/2fpkdjQGPGFgmoJHYZdW7MIfbTD+cIMJ+4OFNjF7qR2DgVhoI9sX/XgMQ89+ndlR0aVVyhE7REREYcSAJ4xFy8kJzpe7tKJS8orLw/HURERExIAnvBme9JQEyUxx9iKyW4uIiCh8mOEJY8CTGB8vTTKT1b+5vAQREVH4MOAJY5dWUkKcNMl0rpLODA8REVH4cJRWGDM8CfFVAc+BYyWy47cC+WLrAfnsh9/kcGGpNM5IlkYZyWo7XOLjnD8T8W9cF+f8mRQfJ2nJCZKenCjpyQnq33jc1jlp0rRBii3nGSEiIqoPBjxhDHiSEuKlSQNnl9bE9zeG5Lly0pPk3M7NZMjJLeT8zs1UgERERBTrGPCEQaXDGXQgU9M5N8u4HqO2erRuKL/r1FTaN8lQMzAfKSyTCodDKisdaoV1878rXb+XlTukqKxCCksrpKisXApKKlTGaM/RIjlcWCb/+fYXdRnQpbk88afukp1mv3lHiIiIAsGAJ5xFywnxcuVpraVLiwaSnZYs7Rqnq+uCpbS8Ur7ddVg+2rhX3li5Sz7etE/+8MwX8vLw06Rjs8ygPQ8REVG0YdFyWLu04lSA06ttIxWABDPYgeTEeDmjQ2OZeFFX+ddNfeW4hmny08FCuePddZzokIiIYhoDnjCO0kIBcrh0a9VQ3rv5TFXUvO7nI/LRxn1he24iIiKrYcATzgxPfHhf7mZZqfLns9qrfz++aItU6MiLiIgoxjDgCYOqGp7wj5ga2a+DKlretj9f/vPt7rA/PxERkRUw4Alz0XK4IdgZ1f949e/nPtse9ucnIiKygpCegT///HO56KKLpGXLlmoyvHnz5rndfv3116vrzZcLLrjAbZtDhw7JNddcI1lZWdKwYUMZMWKE5Ofnu22zfv16OeeccyQ1NVVat24t06ZNEyupqHT+xISBkXB1nzaCuQh/PFAg+/OKI7IPREREtg14CgoKpHv37jJr1iyf2yDA2bNnj3F566233G5HsLNx40ZZvHixfPDBByqIuvHGG43b8/LyZNCgQdK2bVtZvXq1TJ8+XSZNmiQvvPCCWIUr3glr0bJZVmqSdGreQP37m58OR2QfiIiIbDsPz5AhQ9SlJikpKZKbm+v1tk2bNsnChQvl66+/lt69e6vrnn76aRk6dKg8/vjjKnP0xhtvSGlpqcyePVuSk5Ola9eusnbtWpkxY4ZbYGSJDE8EurS009o1ks17j8k3Ow/L0FNaRGw/iIiIYrKGZ+nSpdKsWTPp1KmTjBo1Sg4ePGjctmLFCtWNpYMdGDBggMTHx8vKlSuNbfr166eCHW3w4MGyZcsWOXzYGtmMSBYta73b5aif3/x0KGL7QEREFJMzLaM769JLL5X27dvL9u3b5d5771UZIQQxCQkJsnfvXhUMmSUmJkqjRo3UbYCfuL9Z8+bNjdtycpwnerOSkhJ1MXeLQVlZmboEEx5PjwZHj1awH99f3Y9zdmlt/DVPjhYUqYVHg0G3J1LtCge7t9Hu7YuFNrJ90Y/HsG4C+ZuOaMBz5ZVXGv8+5ZRTpFu3bnL88cerrM/5558fsuedOnWqTJ48udr1ixYtkvT09KA/X4VrLa09v+yWBQt2SSQ4HCLZyQlytFTkxX8vlhOygzsnD2qs7M7ubbR7+2KhjWxf9OMxDExhYWF0rqXVoUMHadKkiWzbtk0FPKjt2b9/v9s25eXlauSWrvvBz3373GcR1r/7qg0aP368jB071i3Dg9FdKH7GaLBgR58fzv5Y/bt92zYydOhJEimL89fL/A17JSH3RBl6rnOoejDahz/QgQMHSlKSPRcptXsb7d6+WGgj2xf9eAzrRvfQRF3As3v3blXD06KFs6i2b9++cuTIETX6qlevXuq6Tz75RCorK6VPnz7GNvfdd596s+gPMnywoSbIW3eWLpTGxRPuH4oPQ71aekpSYkQ/bE9r30gFPGt25wV9P0L12lmJ3dto9/bFQhvZvujHYxiYQP6eQ1q0jPlyMGIKF9ixY4f6965du9Rt48aNk6+++kp27twpS5YskYsvvlg6duyoio6hS5cuqs5n5MiRsmrVKvnyyy9lzJgxqisMI7Tg6quvVgXLmJ8Hw9fnzp0rM2fOdMvgWKZoOULD0rXe7Rqpn2t+OsxlJoiIKKaENOD55ptvpGfPnuoCCELw7wkTJqiiZEwY+Ic//EFOPPFEFbAgi7Ns2TK37AuGnXfu3Fl1cWE4+tlnn+02x052draqvUEwhfvfcccd6vGtMiQ90jMtm3XObSBpSQlyrKRcdh4siOi+EBERhVNIu7T69+8vDlTL+vDRRx/V+hgYkfXmm2/WuA2KnREoWVWlRTI8CLha5aTJ1v358uuRIjm+aWZE94eIiChm5uGJBVaYh0dr0TBN/dxzhEtMEBFR7GDAE8aAJ5IzLWstslLVzz1HGfAQEVHsiPwZOAZYpWgZWjTUAU9RpHeFiIgobBjwhHEtrUgXLUOLbGZ4iIgo9kT+DBwD9GrplsjwZLtqeJjhISKiGMKAJ6wZHisEPMzwEBFR7GHAE86i5fh4y4zSOlZcLvkl5ZHeHSIiorCI/Bk4Bhjz8Fggw5OZkigNUp3TL+1ltxYREcUIBjxhzPAkWKCGx9yt9Svn4iEiohjBgCcM9OKhVpiHB1i4TEREscYaZ2Cbs9I8PMDCZSIiijUMeGJspmW3DA+7tIiIKEZY4wxsc1ZaS8stw5PH5SWIiCg2MOAJ4ygtyxQt6+UljnB5CSIiig0MeGK4S2svFxAlIqIYYY0zsM1ZtWj5WEm5HCsui/TuEBERhRwDnhjM8GSkJEqWa/LBPczyEBFRDLDGGdjmrFa07D4XDwuXiYjI/hjwhHNpCYt0aZkLl7m8BBERxQIGPGGt4bHOy90oPVn9PFLIGh4iIrI/65yBbayy0npdWllpSepnHouWiYgoBjDgicGiZdBFy3lF5ZHeFSIiopCzzhnYphwOh1RKnOVqeHSG52gRu7SIiMj+GPCEWLmuWLZYDQ+7tIiIKJZY5wxsU+W6P8tqNTyprhoeZniIiCgGMOAJsXJdsWy1gCfNVcNTzBoeIiKyPwY8IVZmyvAkWalLixkeIiKKIdY5A9tUhauGJy5OJN5CRcvZHJZOREQxhAFPmIqWrTRCy5zhKS6rlJLyikjvDhERUUgx4AmxsopKy83BAw1SE1XWCTgXDxER2Z21zsI2pEdpWS3Dg+61zBRduMy5eIiIyN4Y8IRplJaVRmhpLFwmIqJYwYAnbDU81nupqyYf5NB0IiKyt5CehT///HO56KKLpGXLlhIXFyfz5s2rtuzChAkTpEWLFpKWliYDBgyQrVu3um1z6NAhueaaayQrK0saNmwoI0aMkPz8fLdt1q9fL+ecc46kpqZK69atZdq0aWIVVu3SMq+nxeUliIjI7kIa8BQUFEj37t1l1qxZXm9HYPL3v/9dnn/+eVm5cqVkZGTI4MGDpbi42NgGwc7GjRtl8eLF8sEHH6gg6sYbbzRuz8vLk0GDBknbtm1l9erVMn36dJk0aZK88MILYqkMjwW7tIyh6ZxtmYiIbM75FT9EhgwZoi7eILvz1FNPyf333y8XX3yxuu61116T5s2bq0zQlVdeKZs2bZKFCxfK119/Lb1791bbPP300zJ06FB5/PHHVebojTfekNLSUpk9e7YkJydL165dZe3atTJjxgy3wCjSo7Ss3aXFomUiIrK3kAY8NdmxY4fs3btXdWNp2dnZ0qdPH1mxYoUKePAT3Vg62AFsHx8frzJCl1xyidqmX79+KtjRkCV67LHH5PDhw5KTk1PtuUtKStTFnCWCsrIydQmm4lLn4yXGOx/fSjKTnUHYkYKSOu+bvp/V2hZMdm+j3dsXC21k+6Ifj2HdBPI3HbGAB8EOIKNjht/1bfjZrFkzt9sTExOlUaNGbtu0b9++2mPo27wFPFOnTpXJkydXu37RokWSnp4uwbTpMLqyEqSwIF8WLFggVrL3Z+e+bfjhR1lQvq1ej4UuR+1oqUh+mUhumojFph8KWhvtyO7ti4U2sn3Rj8cwMIWFhdYPeCJp/PjxMnbsWLcMD4qdUQuE4uhgStywR2Tzd5KTnSVDh/YVK9m/4idZuHuL5DRtKUOHdqtzdI0/0IEDB0pSUpLkl5TLOdM/Vz+TE+PljPY58uxVPSQlKUGilWcb7cbu7YuFNrJ90Y/HsG50D42lA57c3Fz1c9++fWqUlobfe/ToYWyzf/9+t/uVl5erkVv6/viJ+5jp3/U2nlJSUtTFEz4Ig/1h6IhzpjiSExMs90Gbk5GqfuaVlNd73/Rrd/hIiQp2oLS8Uj7felA27S+U09o1kmgXiveHldi9fbHQRrYv+vEYBiaQv+eIdTigGwoByZIlS9wiNdTm9O3rzITg55EjR9ToK+2TTz6RyspKVeujt8HILXM/Hr7JderUyWt3VrhFxSitIM7Dg7W5oElmspx8nDNbdrTQnnUTREQUPUIa8GC+HIyYwkUXKuPfu3btUvPy3HbbbfLQQw/J+++/L999951cd911auTVsGHD1PZdunSRCy64QEaOHCmrVq2SL7/8UsaMGaMKmrEdXH311apgGfPzYPj63LlzZebMmW5dVlYIeBIsPA/PsSAOSy92LUSampQgOenOQnKOAiMiokgLaZfWN998I+eee67xuw5Chg8fLnPmzJG77rpLzdWD4ePI5Jx99tlqGDomENQw7BxBzvnnn69GZ1122WVq7h7zyC4UG48ePVp69eolTZo0UZMZWmFIOpTrxUNjZFh6cVlVwGM8Puf5ISIiOwc8/fv3V/Pt+IIsz5QpU9TFF4zIevPNN2t8nm7dusmyZcvEiqzcpVUVkJSr44TjUV8l5c4ALyUxvmqtLi5dQUREEWa9tIPN6AyPlZeWKK2oNGpv6qvELcPjWo2dGR4iIoowBjwhVmZkeKz3UmemJIqOw4LVraUDp9SkqgwP1+oiIqJIs95Z2GYqjNXSrZfhQRdWsOtsSlxFyymJphoeLl1BREQRxoAnXKulW7CGB6rqbEKR4dFdWsEb9k5ERFQXDHhiePFQqKqzKQ/uKK3EBNM8P5yHh4iIIsuaZ2Eb0aO0kmIkw2OM0jIPS2fAQ0REEcaAJ1xdWhas4QGdhQlWYbHO8LgNS2eXFhERRRgDnhArr6y07EzLUBWUBLuGxzQsvbhMKl2ZLiIiokhgwBO2Li1rvtRVQUl5UEdpmYelY+7J/FIWLhMRUeRY8yxsI1bv0gpVhgfD0pHlQddWMB+fiIioLhjwhKlLy4oTD0IDPXQ8WMPSTRkez+UriIiIIsWaZ2EbMdbSsmiGJyPFGfAUlDgDlfoqMdXwgDEXD0dqERFRBDHgifGJB3XAUxikGpuqmZbdMzxcXoKIiCKJAU+M1/AEO8NjTDxoZHiCWyNERERUFwx4QqzM4jU8GcnOwKQgSBke89ISUDX5IGt4iIgocqx5FrZhhifJohme9OQg1/CYFg+FbGPpCmZ4iIgochjwhGm1dKtOPJiRkhDUGp5qGZ4gL11BRERUFwx4Yr1LyyhargjKbMhVS0u4anhYtExERBZgzbOwjVi9SyvD1aUFRa5gJRiLh1YvWmYNDxERRQ4DnnDNw2PRYenoeopz7VpBSXlQFw8F83paREREkcKAJ8TKK6zdpRUXF2dkeQpK65fhcTgcNWR4GPAQEVHkWPMsbCNWn2nZXLhc3wyPDnbMRcvZrhqeYxyWTkREEcSAJ8YnHoSM5KrC5WAsKwEsWiYiIithwBO2xUOtG/CkBynDoxcORWyX5GqvXksrv6Tc6N4jIiIKNwY8IVZmjNKy7kttTD5Yz7l4zAuHojYIGrhqeHTQQ0REFAnWPQvbhNVHaUGmnounnrMt6wyPHqEFyYnxkuYqYObQdCIiihQGPDE+0zKkB2k9Lc+FQzVduMyh6UREFCkMeEKszFW3YuUurYxgFS17DEnX9Fw8Rzk0nYiIIsS6Z2GbiIYuLV20XN8aG89JBzXOxUNERJHGgCdsEw9GQw1PfQMeZ1tTqmV42KVFRESRxYAnbBMPRsMorfp2ablqeKpleNilRUREkWXds7BNRNNMy4X1Llr2XsOTblqRnYiIKBIY8IR4hJbDGe9YuktLZ3jy6zss3UcNT2pigltAREREFHMBz6RJk9QkdeZL586djduLi4tl9OjR0rhxY8nMzJTLLrtM9u3b5/YYu3btkgsvvFDS09OlWbNmMm7cOCkvL7fMCC2rd2ll6gxPkIqWPTM8el0tfTsREVG4Ob/aR1jXrl3l448/Nn5PTKzardtvv13mz58v7777rmRnZ8uYMWPk0ksvlS+//FLdXlFRoYKd3NxcWb58uezZs0euu+46SUpKkkceeUQiCd1Y0y49WdasW2ec9O1dw6O7tDwyPK4ASNf4EBERxWTAgwAHAYuno0ePyssvvyxvvvmmnHfeeeq6V155Rbp06SJfffWVnHHGGbJo0SL5/vvvVcDUvHlz6dGjhzz44INy9913q+xRcnKyREpiQrxc0rOlpOxZK0kJ8bav4SkxurR8ZXjYpUVERDEc8GzdulVatmwpqamp0rdvX5k6daq0adNGVq9eLWVlZTJgwABjW3R34bYVK1aogAc/TznlFBXsaIMHD5ZRo0bJxo0bpWfPntWer6SkRF20vLw89RPPhUsw6ccL9uMGU7IrFsPioYHup7l9BSXOf2PiZvPj6IRPYUnwX99wiIZjWB92b18stJHti348hnUTyN90xAOePn36yJw5c6RTp06qO2ry5MlyzjnnyIYNG2Tv3r0qQ9OwYUO3+yC4wW2An+ZgR9+ub/MGARWexxOyRagDCoXFixeLVR0sxv8TJa+wRBYsWFDn9v2wA5FNvPy8c4csWLDduO2H/SjYTpCff91b58e3Aisfw2Cwe/tioY1sX/TjMQxMYWFh9AQ8Q4YMMf7drVs3FQC1bdtW3nnnHUlLSwvJc44fP17Gjh3rluFp3bq1DBo0SLKysoIefeINPHDgQFVXZEUHC0plypqlUloZJ4MvGBLQul/m9n0x/weRvb9I184nytD+HYxtKtbvkbe2fyfZjZrI0KG9JdpEwzGsD7u3LxbayPZFPx7DutE9NFER8HhCNufEE0+Ubdu2qQ+n0tJSOXLkiFuWB6O0dM0Pfq5atcrtMfQoLm91QZCSkqIunvBBGKoPw1A+dn01zKiqLyqXeElNCvxtgbaVVjjH4Gekurc1IzXZKGq26msQ7ccwGOzevlhoI9sX/XgMAxPI37PlKmnz8/Nl+/bt0qJFC+nVq5dqzJIlS4zbt2zZooaho9YH8PO7776T/fv3G9vgmxwyNSeddFJE2hBtMG+OTurUZ2h6iY+lJfS8PCxaJiKiSIl4hufOO++Uiy66SHVj/frrrzJx4kRJSEiQq666Sg1DHzFihOp+atSokQpi/va3v6kgBwXLgG4oBDbXXnutTJs2TdXt3H///WruHm9ZHKoOcx9lpCTKseJytYBoszq+SMU+lpbQw9L17URERDEX8OzevVsFNwcPHpSmTZvK2WefrYac49/w5JNPSnx8vJpwECOrMALr2WefNe6P4OiDDz5Qo7IQCGVkZMjw4cNlypQpEWxV9MlIdgY89Vn+wZhpudrEg655eDgsnYiIYjXgefvtt2u8HUPVZ82apS6+IDsUzaN/rCDdNRcPhqbXe+LBahkezrRMRESRZbkaHoqMDNdsy/XL8HhfPLRqLS12aRERUWQw4CG32ZYL6jHbctVMy75qeDjTMhERRQYDHlIy9HpawejS8rF4KFaPNy+oSkREFC4MeEhJT9EBT/2LlqsHPFW/s1uLiIgigQEPKRnJ9V9AtNhHl5b5d87FQ0REkcCAhxTMwwMF9Sla9tGlhXl+qiYfZOEyERGFHwMecsvw1LWGp7yiUtXomGt2vM7Fw8kHiYgoAhjwUFBqeMwjsFJcw9C9z8XjvWh5f16xvLTsRzlaVMYjQkREQceAh4JSw6OHpHur4XEbmu6jS+vFZT/KQ/M3yZ3vruMRISKioGPAQ0Gp4dFD0pOxEKleidTL5IN6O0+/5Zeqn4u/3yeLNu7lUSEioqBiwENKup5puY41PMYsy16yO/4sL2HOLE16f2O95gMiIiLyxICH3GZaxmrpdaFXQvdcOFTT1/uq4TEvafHr0WKZ9ek2HhkiIgoaBjzknuGpZ5eWtxFaUNuw9CLX8w7o0lz9/OrHgzwyREQUNAx4SMk0RmnVtWhZd2l5z/BUraflq0vLef3xTTNc+8H5eoiIKHgY8JCSbozSqqhnl5avGp6au7SKXJmfpg1S6tW1RkRE5A0DHnILeBB4VLomEAxqhifRv6JlHfB4W7Xd4XDI9gP5ddo/IiKKbQx4yK2Gx5xtCcayEtVmWvYZ8LhneNC1hgDH7OUvdsj5T3wmb6zaxaNGREQBYcBDRrFxnGv6nLp0a5XqLq3ahqV7mYcHgY0uWm7mCnjKKhzV5uzZui9f/Vz2wwEeNSIiCggDHjIW+ExPqvtsyzo4qb2Gp3owVVpRKeWubqommc6Ax1sBdV6xc9mJ7345yqNGREQBYcBDQVlPywh4ahul5SXg0dkdPeNzmmtbz/04VuwMgPYcLZYDx0p45IiIyG8MeMhL4XLgGZ5SvbREQm3z8FTv0tJdaEkJcZKUEC+ZqYleR2odc2V4YAOzPEREFAAGPFStcLleGZ46dGnpgEc/v54TqHrAU/X7+t3s1iIiIv8x4CEvK6bXpWhZd2nVEvB4KVrWXVo6w6SXuahew1P1+3e/HOGRIyIivzHgIUNacv2LlrFaeqCLh+rn08+fkew9w6OLloEZHiIiCgQDHjJk1GM9rVqLlhN9z8NTWOae4WngquExZ3hKyiuMLBLsP1Yi+/KKefSIiMgvDHjIy/ISdSharvCvS8tzbh23Lq2kRGOklmeGx1y/c0KzTPWTWR4iIvIXAx4ypBu1M3Wo4SmrT5dWhXuXVg0BDwqae7RuqP793e7qdTyrfzokF/59mXy76zCPLBERGRjwkEGPkqrL0hLocvJvHh5vGZ5ytwyTt5Xb9ZD0rNRE6dYqW/17vZeh6U9/sk02/pon737zM48sEREZGPCQQQccnqOjgtGlZczD4wqMzAo8MzxG0XLVtnlFzn1qkJokxzd1dmn9crioWlHzl9t+U//+6WAhjywRERkY8FD1iQfrUbTsu0vLn3l4XBmeVN8ZHhQ06wVGD+S7z7b86eb9ag0uYMBDRERmDHio+sSD9VlLy1eGx6jhqay2CnpVl5aeeDDBZw0PAh693taRwjK3kVsfbdxr/PvXo0VGNxsREREDHvIySqvuEw/WluHxNlLLM8PjrWhZz8GDLq3stCS1DAUcLCgxMkdLt1Stoo6Y6udDzi6vVTsOyfLtzq4uIiKKTbYKeGbNmiXt2rWT1NRU6dOnj6xatSrSuxRV0sMwD4/a1qNwufpMy966tJz/zkpLlPj4OCPLoxcRXbb1N7XfLbNTpXNuA3XdrkMFKmi69uWVcs1LK+WLrQx6iIhilW0Cnrlz58rYsWNl4sSJ8u2330r37t1l8ODBsn///kjvWtTQSzrUK+DxsZYWMjLxzqRMtcLlqmHpzkCngZeAx5zhAaOOxxXw6O6sQV1zpV3jDKOO57vdR9W+IeNz+ztrq62yfqSwlF1fREQxwHlmsYEZM2bIyJEj5YYbblC/P//88zJ//nyZPXu23HPPPZHePdtPPFhWy2rpcXFxqlsLwY1n4bIx03JS7fPw6FmYm3pkePS8O/07NZUVPx40Ah5zjQ+2vePddfLidb1UJur1r36Syf/bqPZrcNdcOatjY2mYlixHi8rk40371IrsZ3RoLFf2Pk6W74uTmTO/kJJyh1xzRhsZdFJz1YX21Y8HpU2jDDmtXY7KPP18qFAFWE0ykyUrNUkN8dcBXXycM+iLUxf8O04QAzp/d71O+C/O/XXTv6rtXL85/111P3WL8Rju96tNeUW5bDgUJymb90tiQv0/Esz7H26+nru8vEI2Ho6TtC0HJNFHFrJez1vbq+3Ha1LbJjimvlSUl8umI3HSYNtvNR5Df45NbW3x7zHqu4H7flRUlMvWo3GycschSXC1Ly7Er2nVNvV7Dn8eA4+CY7jjmMiaXUckITEx6PsR509ba32Mur9/ysvLZXeBRJQtAp7S0lJZvXq1jB8/3rguPj5eBgwYICtWrIjovsXaaul6gkFvqgKeSr/m4XEPeHxneFAErYeod2iSKb8ecS458dPBAiMgurTncTL/uz3y+Q8HpP/0pWrywg83OLNCZRXl8q/Vu9XF086DhfL215jTB/vmHOo+beEWdTGb/eUOiW4J8uKWtWJvCfLC5jViXwny/KZvxb4S5JnvvxF7S5SnNti3FCMxLkFuvDyCzy828Ntvv0lFRYU0b97c7Xr8vnnz5mrbl5SUqIuWl5enfpaVlalLMOnHC/bjhkJSvMMIQPzdX72dHhEVL5U+76tHcOUXlUhZWapxve66QryD+ya79gOBUVFxiSQmxEtekfMx0xPj1DY56c637r68ItlzuEAFXMieNMlIkOOyk9VtO38rkDLX/EDDeuTKoJOaysT/bZI9R4tlz1FnsDN2QEfp3TZHBT/bDxTI0eIy9Q3lrOMbyynHZcn/1u+RxZv2S2aiQ24+9wRplJkqLyzbIdsOFEjvtg3l3E5N5ZcjRbJm11FJTIiT1jlpKrA7WFCqslJpSQmSlhSvvl0hMEPLKh0OqcQ/HM5/6zFr6HZTv7muUJu4RrQ5/139+qr76X87H89jIFyNcJ+jR49Kdna2X98Ca3wsY08iwFFbG/MkOzur3m0M4Gldz+3PYzjq9RhoX17eMWmQ1cDnN2y/jkwtT+TPY9S6r348iudj4D4F+QWSkZmh2heM/fCnNbW3JXivB7YrKiyUtPT0apmUeu+HH29CR/0fosbHwHu0vKQ4ZOfYmAl4AjV16lSZPHlytesXLVok6enpIXnOxYsXi9XlleL/iapLa/78BQF1TTgXBY2TLz5bKq54o5qKEmRJ4mTpF1/K7qyq63877Lx+3epVkr9VxJkscr41581fKIhtdu9zbrNlw1pZ8Msa2b8HO5cg3239Sd4tQHYlUbKTHLL4o4VySMWyiSrDU+n68P91w0pJSxS5q4vIiv1xsv5QnPTLdUjbgs1y4HuR3vEivc3xcvkhqfhJZGi2yIDTRNBTl3Bss8gxkTHHi1R0EEmM/03k2G/SKkGkT3uPxlbFc1HkkNif3Zccqb7cir04v5za2zGxs8VBPhcWFhbGVsDTpEkTSUhIkH379rldj99zc3OrbY+uLxQ4mzM8rVu3lkGDBklWlulMHKToEwd44MCBkpTk7I6xKmRaHlj9iTgkTs4bONiY+bi29i1ctNgILIYMGigN072387kdK2T/3mPSs9fpcs4JTYzrH9nwmUhxiZzX72zp2tL5+t/zzWI1ieCZ/c6Vlg3T5Ikty0QKiuS8s8+QXm1zJG7DXvn3zvWS1KCRtOnSWmTDd9KxZSMZOvQ0qah0yCPrPhZdKtShSbpc9oezjef7g42PYV3YvX2x0Ea2L/rxGNaN7qGJmYAnOTlZevXqJUuWLJFhw4ap6yorK9XvY8aMqbZ9SkqKunjCB2GoPgxD+djBkmUqdix1xEmWn/trnlYnMy1Fkkxz7pjpAKrMEef2Wuii5QbpKcb1qOM5XFgmJZXObfUyEzmZaer3FjnOkVi/5ZfKr87UlLRulO58nfHvnHT58TdnhVz31jlBee2j4RjWh93bFwttZPuiH49hYAL5e7bNsHRkbF588UV59dVXZdOmTTJq1CgpKCgwRm1R7TDKCDUngS4vYQ54fE08aJ6Lx3OUluc8PJ4jtdD3a15awnOU1m5XwXKrnKruyLaNq/6tFxslIqLYZYsMD1xxxRVy4MABmTBhguzdu1d69OghCxcurFbITLXPxYOh1IEsL6EHXSXGx0mCnmzHCz2CyzzTMoaNl6sK3qpRYp4rpmN7vUaWDniauEZpYeHRH/Y5+7xb5aQZ92+r5uJxzrzcrVVDHnYiohhnm4AH0H3lrQuL/Ke7nQKZfLDcVZrvax0tz+UlnAXOUi2TZM7wmAMePekgiqgzXEFRRjJGPzmDM8yXI65uLM8MD4IwXRdERESxyzZdWhQcGXp5iQDm4tEJm5q6s9xXTK/K8BSWlRszMSeZJi3UXVoY2p1XVG4EQeh2Awwt1nPx6IyROcPTpUWW0Z1lXseLiIhik60yPBS8DE9durR8raOl6QyQuYbHWFbCIygxZ3h0/Q5mLjZDwLPrkHNIIrrSWmRXjQXHDMnPXXOqnMTsDhERMeAhTxmuDE9RXbq0aphl2S3DY1pLq6pgOdHrul6o0fFcVkLThcuQm5WqJig0G3JKC7/bQERE9sYuLXKTXocMT3llXI3raGk6IHLr0vIyQstzlJaxUrpHhqdJg6oZDs3dWURERJ4Y8JAbHXgEkuExurRqy/B4GZauFyr1nORQr5ieX1zVpVU9w1PVhWUekk5EROSJAQ+5STdqZ+oySish8KLlWjI85lFa1QIeV9EytG7EDA8REfnGgIfcpLuCEj16KqBRWrV0ael5eMw1PEbRcrUanupdWnqldG8BDzM8RERUEwY85DXDE8iwdL+7tLzOw1NuzKvjdZRWqamGB6t/+gx4mOEhIiLfGPBQ/YuW/Zx40Fi2wtuwdB8BD2p4qrq0asrwMOAhIiLfOA8PucmoR9Fyci01PFV1ORV1GqXlWcPTrEGKGo6OOXjwk4iIyBcGPORG19Jg/ptAa3hqy/DouXUQxGg62+M5D0/VxIMVklfkPcODmZkXj+0n8XFx1ebgISIiMmPAQz4yPHUoWq4l4GmQkmR0U1Ublu4x07IOjg7kl8jhwlKvGR5vQRAREZE3DHjI+9ISgRQtO+L8yvBkplYNNa+tS0tvW1HpUBcsAHpq6xweLSIiqhMGPOS1dsZcWOx/l1bNNTxGIXJpuVRWOtRCoFVLSyRUWzbij71aqUkHrz2jnZzVsbFaMJSIiKguGPCQ91FapixMsLq0dMDjcGCenwr1u695eBDcPH55dx4dIiIKClZ6khtdPFynpSUSa594ECOqzHU8uobHM8NDREQUTAx4yGvRMubhcSAVE8R5eJC1Mbq1XBkkX/PwEBERBRMDHnKjA49Kh0iJ7qsK0rB08Ax4dCYpw6NLi4iIKJgY8JAb83w4Ovvif5dW7Vka8wzKcNQ1x052GoeXExFR6DDgITeosdGZGn8Ll40urVrW0jIPN88vKVNdZkcY8BARURgw4CGfQ9P9zfCUV8b5tVq6e5dWherWwhw70DCdGR4iIgodBjxUjZ71uNjPuXj8XS3dLcNTXGZ0Z2E4u15JnYiIKBQY8JDX4eOBTD5YNUrLjxoeV40QsjtHCp0BT0PW7xARUYgx4KFqUgPM8Pg78aA5w3OspNxYFJQFy0REFGoMeCh4XVoBDEtHQbQuWGb9DhERhRoDHvI5F08ourT0iucYls4h6UREFC4MeMhnl1ZRaWXQu7T0CDBzDU92WjKPAhERhRQDHqp3DU9durQQ8DDDQ0RE4cKAh6pJC2CUFiYPrHDE+R/wpJoDnlL1b9bwEBFRqDHgoXoVLZea1tvyp0urgWlpCWZ4iIgoXBjwUDWpumjZj5mWzQuM+jUPT6qXeXg4yzIREYUYAx6qJtUVuBSX+5HhqXAGPHFxIkkJzq6tmmSYJh7UGZ4sTjxIREQhxoCHfA9L92OUls7wYB2tOEQ9fg5LLy6rlIP5rhoeBjxERGTngKddu3bqJGm+PProo27brF+/Xs455xxJTU2V1q1by7Rp06o9zrvvviudO3dW25xyyimyYMGCMLYitmt4SlxDtPwpWDYPS4d9x4rVT860TEREts/wTJkyRfbs2WNc/va3vxm35eXlyaBBg6Rt27ayevVqmT59ukyaNEleeOEFY5vly5fLVVddJSNGjJA1a9bIsGHD1GXDhg0RalGMFS27urT8KViGpAQsFOrc1uGasJABDxERhVrV1+0IadCggeTm5nq97Y033pDS0lKZPXu2JCcnS9euXWXt2rUyY8YMufHGG9U2M2fOlAsuuEDGjRunfn/wwQdl8eLF8swzz8jzzz8f1rbYRUoAw9J1l5a/GR49F09xmbM7CxjwEBGR7TM86MJq3Lix9OzZU2VwysvLjdtWrFgh/fr1U8GONnjwYNmyZYscPnzY2GbAgAFuj4ltcD3VL8PjX8BTUaeAx/zvxISIvw2JiMjmIprhueWWW+TUU0+VRo0aqa6p8ePHq24tZHBg79690r59e7f7NG/e3LgtJydH/dTXmbfB9b6UlJSoi7nrDMrKytQlmPTjBftxQ8mV4JGikvJa97uo2Hk7Rmj528aMlKrh69lpiZZ/baLxGAbC7u2LhTayfdGPx7BuAvmbDnrAc88998hjjz1W4zabNm1SRcZjx441ruvWrZvK5Pz1r3+VqVOnSkpKioQKHn/y5MnVrl+0aJGkp6eH5DnRzRYtflTxX6IcPHqs1gLwDYcwMitBigvy/S4WL8lHROWKqkqLoqbIPJqOYV3YvX2x0Ea2L/rxGAamsLAwcgHPHXfcIddff32N23To0MHr9X369FFdWjt37pROnTqp2p59+/a5baN/13U/vrbxVRcEyCSZgy1keDACDAXSWVlZEuzoE2/ggQMHSlJSkkSDjb/mycyNX0l8cqoMHfq7GretWPeLyJaN0rRRQxk6tI9fj//fQ2tkW94B9e82uY1l6NDeYmXReAwDYff2xUIb2b7ox2NYN7qHJiIBT9OmTdWlLlCQHB8fL82aNVO/9+3bV+677z71RtAfUvjQQjCE7iy9zZIlS+S2224zHgfb4HpfkD3ylkHCc4TqwzCUjx1smWkpxlw5te1zpV5HKynR7/aZJxpsmJ4cNa9LNB3DurB7+2KhjWxf9OMxDEwgf88RqxZFUfFTTz0l69atkx9//FGNyLr99tvl//7v/4xg5uqrr1bdXBhyvnHjRpk7d64alWXOztx6662ycOFCeeKJJ2Tz5s1q2Po333wjY8aMiVTT7DPxYKhGabkmHwQuK0FERLYuWkaG5e2331YBCgqIUZyMgMcczGRnZ6u6mtGjR0uvXr2kSZMmMmHCBGNIOpx55pny5ptvyv333y/33nuvnHDCCTJv3jw5+eSTI9Qy+4zSwsKglZUOiY/3PYNy3YalV0XkXFaCiIhsHfBgdNZXX31V63YoZl62bFmN21x++eXqQsGhJwbU62mlu9a/CsbEg5BpGqXVMK1qygEiIqJQ4QQo5HPxUH9WTA90aQnPeXg46SAREYUDAx6q/qaIjzMCmGJXl1VtGZ7AanjMRcv2LSAlIiLrYMBDXqXq2ZZry/Do1dKZ4SEiIgtjwEP1WkCUAQ8REUUDBjxUr6HpGMkFKaa6n0CGpbOGh4iIwoEBD9XYpVVbhqfUtXhocqLvoes1FS2zhoeIiGy/eChZf2h6bTU8peUO9TM5gBXPG2ckC6b2QbeZOfghIiIKFZ5tqMYanlq7tFyjtJICCHhyMpLl71f1lKzUJImL8z8zREREVFcMeKheRctldZh4EH7frSVfeSIiChvW8JBXqa6iZSwgGuwMDxERUbjxLEU1zrZcW5dWWYWu4WHXFBERWRcDHvIqLdm/omXdpcUMDxERWRkDHqq5hsc17DzYNTxEREThxLMU1Rzw1DosXWd42KVFRETWxYCHvEpJCqyGh11aRERkZQx4qJZ5eCr969LiKC0iIrIwBjxU41patS8twaJlIiKyPgY8VM+JB3WXFmt4iIjIuhjwUL3W0uIoLSIiigYMeKjG1dJDsZYWERFRuPEsRXXu0nI4HJxpmYiIogIDHqqlaLmy1vod4MSDRERkZQx4qM5dWro7C9ilRUREVsaAh+rcpVXmGpIODHiIiMjKGPBQrRke1OrUNEIrThySEM9h6UREZF0MeKjGGh7EOiWmTI6Zvp7rhhIRkdUx4CGvUk1RjK9uLZ3hSWRyh4iILI4BD3mVmBBvzJ7sa6SWHqXFSZaJiMjqGPBQnUdqGRkevouIiMjieKqi2ldM97G8hK7hYYaHiIisjgEP+cQMDxER2QUDHqo1w1PComUiIopyDHjIp9Tkmmt4StmlRUREsR7wPPzww3LmmWdKenq6NGzY0Os2u3btkgsvvFBt06xZMxk3bpyUl5e7bbN06VI59dRTJSUlRTp27Chz5syp9jizZs2Sdu3aSWpqqvTp00dWrVoVqmbFlLSkeL+KlrlQOhERxWzAU1paKpdffrmMGjXK6+0VFRUq2MF2y5cvl1dffVUFMxMmTDC22bFjh9rm3HPPlbVr18ptt90mf/nLX+Sjjz4ytpk7d66MHTtWJk6cKN9++610795dBg8eLPv37w9V02KvhsdH0XKpa1h6Ypz3mZiJiIhsH/BMnjxZbr/9djnllFO83r5o0SL5/vvv5fXXX5cePXrIkCFD5MEHH1TZGgRB8Pzzz0v79u3liSeekC5dusiYMWPkj3/8ozz55JPG48yYMUNGjhwpN9xwg5x00knqPsgYzZ49O1RNixm1raelu7Q4LJ2IiKwuMVJPvGLFChUMNW/e3LgOmRlkhDZu3Cg9e/ZU2wwYMMDtftgGmR5AYLR69WoZP368cXt8fLy6D+7rS0lJibpoeXl56mdZWZm6BJN+vGA/bjikuqZQPlZc6nX/i0vLjGHp0di+WDiG/rB7+2KhjWxf9OMxrJtA/qYjFvDs3bvXLdgB/Ttuq2kbBChFRUVy+PBh1TXmbZvNmzf7fO6pU6eqDJS3rBOyQ6GwePFiiTa/7UECMF7WbdgiC/I2Vbt97V4ERAlqaYlobF+g7N5Gu7cvFtrI9kU/HsPAFBYWhibgueeee+Sxxx6rcZtNmzZJ586dxcqQEULdj4YAqnXr1jJo0CDJysoKevSJN/DAgQMlKSlJosn3i7bKsn07pEWbdjJ0aPVjuvfLnSI7flBFy9HYvlg4hv6we/tioY1sX/TjMawb3UMT9IDnjjvukOuvv77GbTp06ODXY+Xm5lYbTbVv3z7jNv1TX2feBkFJWlqaJCQkqIu3bfRjeIMRX7h4wgdhqD4MQ/nYoZKVnqx+FpU5vO57hTi7vJDhicb2BcrubbR7+2KhjWxf9OMxDEwgf88BBTxNmzZVl2Do27evGrqO0VQYkg74BoZgBsXHepsFCxa43Q/b4HpITk6WXr16yZIlS2TYsGHqusrKSvU7CpypftJd8/Dkl7pPFaCVlbsWD+VsTkREZHEhO1Vhjh0MJcdP1Nng37jk5+er29F9hMDm2muvlXXr1qmh5vfff7+MHj3ayL7cdNNN8uOPP8pdd92lanKeffZZeeedd9ToLw1dUy+++KIa1o7uNBQ9FxQUqFFbVD8ZKc54uLCkvObFQ52JHiIiIssKWdEy5tNBEKJh1BV8+umn0r9/f9UV9cEHH6gABRmbjIwMGT58uEyZMsW4D4akz58/XwU4M2fOlFatWslLL72kRmppV1xxhRw4cEA9H4qcMcR94cKF1QqZKXCZroCnoMTXPDwMeIiIKMYDHkwi6G1WZLO2bdtW67LyhOBozZo1NW6D7it2YYUuw5PvI8NjLC3BLi0iIrI4nqrIp8wUZw1Pga8aHiPDw5mWiYjI2hjwkE/pybpLixkeIiKKbgx4qM41PCxaJiKiaMGAh2qt4cFq6RWV1butylyLh2JpCSIiIitjwEM+ZbhqeHzV8ZRw8VAiIooSDHjIp+SEeEmMj/NZx8MuLSIiihYMeMinuLg4o1urpoCHw9KJiMjqGPBQnQuX9Tw8nGmZiIisjgEP+VXHwwwPERFFMwY8VOfZlktdo7SY4SEiIqtjwEM1ytCTD3oZpVVa7uzm4kzLRERkdQx4yK8urXwvNTzGPDx8FxERkcXxVEV+dWkVclg6ERFFMQY85OcorRqGpXOmZSIisjgGPOTXAqLeurQ40zIREUULBjxUo0x/hqUzw0NERBbHgIf8G5buZZSWLlpO5LuIiIgsjqcqqpGvpSWwerpeQZ3z8BARkdUx4CG/ipYLPWp4dHcWcFg6ERFZHQMeqlF6sp6Hxz3DU2oKeJjhISIiq2PAQ/4NS/eo4dELhwKLlomIyOoY8FCdanh0l1ZSQpzEcZQWERFZHAMe8ivD49mlVVbuLFhOYgEPERFFAQY85FeGp7is0hiVBaUVziLmZAY8REQUBRjwkF9Fy551PKVGhof9WUREZH0MeKhGKYnxkhgfV62Op6qGh28hIiKyPp6tqEZxcXFeC5cZ8BARUTRhwEMBFC5XVBuWnsxJeIiIKAow4CG/63gKTRkePfEgu7SIiCgaMOAh/xcQdevS4rB0IiKKHgx4qE6zLVd1afEtRERE1sezFdUqIyWhWg2PeaZlIiIiq2PAQ7XyNkqLNTxERBRNQhbwPPzww3LmmWdKenq6NGzY0OeQZ8/L22+/7bbN0qVL5dRTT5WUlBTp2LGjzJkzp9rjzJo1S9q1ayepqanSp08fWbVqVaiaFZMykr0EPLpLi/PwEBFRLAc8paWlcvnll8uoUaNq3O6VV16RPXv2GJdhw4YZt+3YsUMuvPBCOffcc2Xt2rVy2223yV/+8hf56KOPjG3mzp0rY8eOlYkTJ8q3334r3bt3l8GDB8v+/ftD1bQYzvBU79JiwENERNHAeSYLgcmTJ6uf3jIyZsj+5Obmer3t+eefl/bt28sTTzyhfu/SpYt88cUX8uSTT6qgBmbMmCEjR46UG264wbjP/PnzZfbs2XLPPfcEuVWxKdNVw+N94kHW8BARUQwHPP4aPXq0ytp06NBBbrrpJhW4oGsLVqxYIQMGDHDbHoEOMj06i7R69WoZP368cXt8fLy6D+7rS0lJibpoeXl56mdZWZm6BJN+vGA/bjiluiYXPFZcarSjyBX86B6taG5fLBzDWG5fLLSR7Yt+PIZ1E8jfdEQDnilTpsh5552n6nwWLVokN998s+Tn58stt9yibt+7d680b97c7T74HQFKUVGRHD58WCoqKrxus3nzZp/PO3XqVCMDZYZ9wL6EwuLFiyVa/bgfAU+C7Ni9RxYs+EVd9/3PiHTiZd+vv4h0iO72+cvubbR7+2KhjWxf9OMxDExhYWFoAh50ET322GM1brNp0ybp3LmzX4/3wAMPGP/u2bOnFBQUyPTp042AJ1SQEULdj4YAqnXr1jJo0CDJysoKevSJN/DAgQMlKSlJolHS9/vlze1rJaVBjgwd2kdd9/2irSK7d0iHdm1EZGdUty8WjmEsty8W2sj2RT8ew7rRPTRBD3juuOMOuf7662vcBl1TdYURVg8++KDqbsKoLNT27Nu3z20b/I6gJC0tTRISEtTF2za+6oIAj42LJ3wQhurDMJSPHWoNM1OMeXh0G8qdEy1LSlKiSEV0t89fdm+j3dsXC21k+6Ifj2FgAvl7Dijgadq0qbqECkZi5eTkGMFI3759ZcGCBW7b4Fsarofk5GTp1auXLFmyxBjdVVlZqX4fM2ZMyPYz1mSlOt9QeUVl1UdpYablqsFbRERElhSyGp5du3bJoUOH1E/U2SCYAcylk5mZKf/73/9UJuaMM85Q8+cgkHnkkUfkzjvvNB4DRczPPPOM3HXXXfLnP/9ZPvnkE3nnnXfUKCwNXVPDhw+X3r17y+mnny5PPfWU6hrTo7ao/rLTXAFPcfWAh4uHEhFRTAc8EyZMkFdffdWtRgc+/fRT6d+/v0pDYcLA22+/XRwOhwqE9BBzDUPSEdxgm5kzZ0qrVq3kpZdeMoakwxVXXCEHDhxQz4ci5x49esjChQurFTJT/TM8xWWVasJBZHVKXX1aHJZOREQxHfBg/p2a5uC54IIL1KU2CI7WrFlT4zbovmIXVuhkpla9TY4Vl0njzBRjaQkuHkpERNGAa2lRrRLi44wV0/OKnfPvlLmWlmCXFhERRQMGPOSXLFeWRxcuVy0twZmWiYjI+hjwkF8auOp4jrkyPEaXFhcPJSKiKMCAh/ySlaa7tMrcVktnlxYREUUDBjxUp7l4OCydiIiiCQMe8ktWmvcurSTXwqJERERWxoCH/NIg1b1Lq8w1Dw9reIiIKBow4CG/sEuLiIiiGQMeCrBo2dmlVWIULbNLi4iIrI8BDwWU4cFMy5Bf4gx80pMT+AoSEZHlMeChgObhySsql6LSCjnqGq2Vm5XKV5CIiCyPAQ8FPA/PnqNF6t8ZyQlGMTMREZGVMeChALu0ymXv0WL179zsVImLYw0PERFZHwMeCmxYehEyPM6Ap0V2Gl89IiKKCgx4KLCJB0vK5ZcjRUaGh4iIKBow4CG/mGt1tu3PVz9bMOAhIqIowYCH/JKSmCCpSc63yw/7jqmfzPAQEVG0YMBDAQ9N336AGR4iIoouDHjIb1mubq2yCuc6WrlZLFomIqLowICHAi5c1ljDQ0RE0YIBDwXcpQUpifHSMN09ACIiIrIqBjwUcJcWtGyYxkkHiYgoajDgoTp1aXENLSIiiiYMeCjg5SWA9TtERBRNGPBQnSYf5Bw8REQUTRjwUJ26tJjhISKiaMKAh+pUtJzLhUOJiCiKMOAhvzHDQ0RE0YoBD9Uxw8OV0omIKHow4CG/ZbtqeJIT4qVRejJfOSIiihpVX9mJatGhSaZc0vM4Ob5phsTHx/H1IiKiqMGAh/yGIOfJK3rwFSMioqgTsi6tnTt3yogRI6R9+/aSlpYmxx9/vEycOFFKS0vdtlu/fr2cc845kpqaKq1bt5Zp06ZVe6x3331XOnfurLY55ZRTZMGCBW63OxwOmTBhgrRo0UI914ABA2Tr1q2hahoRERFFmZAFPJs3b5bKykr5xz/+IRs3bpQnn3xSnn/+ebn33nuNbfLy8mTQoEHStm1bWb16tUyfPl0mTZokL7zwgrHN8uXL5aqrrlLB05o1a2TYsGHqsmHDBmMbBEl///vf1eOvXLlSMjIyZPDgwVJcXByq5hEREVEUCVmX1gUXXKAuWocOHWTLli3y3HPPyeOPP66ue+ONN1TGZ/bs2ZKcnCxdu3aVtWvXyowZM+TGG29U28ycOVM9zrhx49TvDz74oCxevFieeeYZFeAgu/PUU0/J/fffLxdffLHa5rXXXpPmzZvLvHnz5MorrwxVE4mIiChKhHWU1tGjR6VRo0bG7ytWrJB+/fqpYEdDZgaB0eHDh41t0EVlhm1wPezYsUP27t3rtk12drb06dPH2IaIiIhiW9iKlrdt2yZPP/20kd0BBCqo8TFDZkbflpOTo37q68zb4Hq9nfl+3rbxVFJSoi7mrjUoKytTl2DSjxfsx7UKu7cvFtpo9/bFQhvZvujHY1g3gfxNBxzw3HPPPfLYY4/VuM2mTZtUkbH2yy+/qG6pyy+/XEaOHCmRNnXqVJk8eXK16xctWiTp6ekheU50w9mZ3dsXC220e/tioY1sX/TjMQxMYWFh6AKeO+64Q66//voat0G9jvbrr7/KueeeK2eeeaZbMTLk5ubKvn373K7Tv+O2mrYx366vwygt8zY9engfQj1+/HgZO3asW4YHI8RQQJ2VlSXBjj7xBh44cKAkJVUtvmkXdm9fLLTR7u2LhTayfdGPx7BudA9NSAKepk2bqos/kNlBsNOrVy955ZVXJD7evWSob9++ct9996kDrT+E8KHUqVMn1Z2lt1myZIncdtttxv2wDa4HdIkh6ME2OsDBC4DRWqNGjfK6XykpKeriCfsQqg/DUD62Fdi9fbHQRru3LxbayPZFPx7DwATy9xyyomUEO/3795c2bdqoup0DBw6omhpzXc3VV1+tCpYx5BxD1+fOnatGZZmzL7feeqssXLhQnnjiCTXUHcPWv/nmGxkzZoy6PS4uTgVDDz30kLz//vvy3XffyXXXXSctW7ZUw9eJiIiIQla0jCwMCpVxadWqldttGEquR1Ohbmb06NEqC9SkSRM1gaAekg7oCnvzzTfVsHPM4XPCCSeo4eYnn3yysc1dd90lBQUF6n5HjhyRs88+WwVJmKiQiIiIKGQBD+p8aqv1gW7dusmyZctq3AbFzrj4gizPlClT1IWIiIjIE1dLJyIiIttjwENERES2x4CHiIiIbC9sMy1bmS6iDmQ8v78w5B4TI+Gx7Tgc1u7ti4U22r19sdBGti/68RjWjT5v6/N4TRjwiMixY8fUi4HJB4mIiCj6zuMY+V2TOIc/YZHNVVZWqhmhGzRooEZ8BZOexfnnn38O+izOVmD39sVCG+3evlhoI9sX/XgM6wYhDIIdzL3nObmxJ2Z4UMgUH19trqBgw4esHT9oY6V9sdBGu7cvFtrI9kU/HsPA1ZbZ0Vi0TERERLbHgIeIiIhsjwFPiGGR0okTJ3pdrNQO7N6+WGij3dsXC21k+6Ifj2HosWiZiIiIbI8ZHiIiIrI9BjxERERkewx4iIiIyPYY8ITYrFmzpF27dpKamip9+vSRVatWidVNnTpVTjvtNDURY7NmzWTYsGGyZcsWt2369++vJmk0X2666Sa3bXbt2iUXXnihpKenq8cZN26clJeXixVMmjSp2v537tzZuL24uFhGjx4tjRs3lszMTLnssstk3759UdM+vOc824cL2hStx+/zzz+Xiy66SE0whv2dN29etQnIJkyYIC1atJC0tDQZMGCAbN261W2bQ4cOyTXXXKPmOmnYsKGMGDFC8vPz3bZZv369nHPOOepvFpMVTps2LeLtw7IDd999t5xyyimSkZGhtrnuuuvUhKm1HfdHH33U8u2D66+/vtq+X3DBBVFz/Pxpo7e/SVymT58eFcdwqh/nhmB9di5dulROPfVUVczdsWNHmTNnTv0bgJmWKTTefvttR3JysmP27NmOjRs3OkaOHOlo2LChY9++fZZ+yQcPHux45ZVXHBs2bHCsXbvWMXToUEebNm0c+fn5xja/+93vVHv27NljXI4ePWrcXl5e7jj55JMdAwYMcKxZs8axYMECR5MmTRzjx493WMHEiRMdXbt2ddv/AwcOGLffdNNNjtatWzuWLFni+OabbxxnnHGG48wzz4ya9u3fv9+tbYsXL8aM6o5PP/00ao8f9uG+++5z/Oc//1Ftee+999xuf/TRRx3Z2dmOefPmOdatW+f4wx/+4Gjfvr2jqKjI2OaCCy5wdO/e3fHVV185li1b5ujYsaPjqquuMm7Ha9C8eXPHNddco97/b731liMtLc3xj3/8I6LtO3LkiDoWc+fOdWzevNmxYsUKx+mnn+7o1auX22O0bdvWMWXKFLfjav67tWr7YPjw4er4mPf90KFDbttY+fj500Zz23DBuSEuLs6xffv2qDiGg/04NwTjs/PHH390pKenO8aOHev4/vvvHU8//bQjISHBsXDhwnrtPwOeEMIH0ujRo43fKyoqHC1btnRMnTrVEU1w8sQf72effWZchxPmrbfe6vM+eBPHx8c79u7da1z33HPPObKyshwlJSUOKwQ8+OD0BieXpKQkx7vvvmtct2nTJvUa4EQTDe3zhGN1/PHHOyorK21x/DxPJmhXbm6uY/r06W7HMSUlRZ0QAB+cuN/XX39tbPPhhx+qE84vv/yifn/22WcdOTk5bm28++67HZ06dXKEk7eTpadVq1ap7X766Se3k+WTTz7p8z5Wbh8CnosvvtjnfaLp+Pl7DNHe8847z+26aDmG3s4NwfrsvOuuu9QXUrMrrrhCBVz1wS6tECktLZXVq1ertLp5CQv8vmLFCokmR48eVT8bNWrkdv0bb7whTZo0kZNPPlnGjx+vVqPW0Eak35s3b25cN3jwYLVezMaNG8UK0N2B1HOHDh1UmhxpVsBxQxeC+dihu6tNmzbGsYuG9pnfi6+//rr8+c9/dlsrLtqPn9mOHTtk7969bscM082jG9l8zNAN0rt3b2MbbI+/y5UrVxrb9OvXT5KTk93ajbT94cOHxWp/lzieaJMZuj/QndCzZ0/VVWLuKrB6+9CNgS6OTp06yahRo+TgwYPGbXY7fujmmT9/vuqW8xQtx/Cox7khWJ+d2Mb8GHqb+p47uZZWiPz2229SUVHhdlABv2/evFmiaWHV2267Tc466yx1YtSuvvpqadu2rQoY0J+M+gL8wf3nP/9Rt+Pk463t+rZIw4kQfcL4YN2zZ49MnjxZ9Ylv2LBB7R8+TDxPJNh/ve9Wb58Z6giOHDmiaiTscvw86X3yts/mY4aTqVliYqL6sDZv0759+2qPoW/LyckRK0CdBI7ZVVdd5bY22C233KLqHtCm5cuXq0AW7+8ZM2ZYvn2o17n00kvV/m3fvl3uvfdeGTJkiDrJJSQk2Or4wauvvqpqYdBms2g5hpVezg3B+uz0tQ2CoqKiIlWjVxcMeKhGKD5DEPDFF1+4XX/jjTca/0a0jkLR888/X31QHX/88ZZ/VfFBqnXr1k0FQAgA3nnnnTr/MVnVyy+/rNqL4MYuxy+W4Rv0n/70J1Wk/dxzz7ndNnbsWLf3NU4+f/3rX1WxqdVnmb7yyivd3pPYf7wXkfXBe9NuZs+erTLLKDyOxmM42se5wcrYpRUi6CrAtxLP6nT8npubK9FgzJgx8sEHH8inn35a62ryCBhg27Zt6ifa6K3t+jarwTeSE088Ue0/9g/dQMiK+Dp20dK+n376ST7++GP5y1/+Yuvjp/eppr83/Ny/f7/b7egqwMifaDmuOtjBcV28eHGtK7/juKKNO3fujIr2maGrGZ+j5vdktB8/bdmyZSqjWtvfpVWP4Rgf54ZgfXb62gbv9/p8IWXAEyKIynv16iVLlixxSwHi9759+4qV4Zsj3tDvvfeefPLJJ9XSp96sXbtW/USmANDG7777zu0DSn9An3TSSWI1GNqK7Ab2H8ctKSnJ7djhwwk1PvrYRUv7XnnlFdUNgCGgdj5+eI/iQ9J8zJD+Rm2H+Zjhgxh1Bhre3/i71AEftsHQYgQW5naj6zPS3SE62EHtGYJY1HjUBscVNS66K8jK7fO0e/duVcNjfk9G8/HzzLric6Z79+5RdQwdtZwbgvXZiW3Mj6G3qfe5s14lz1TrsHSMEpkzZ44aYXDjjTeqYenm6nQrGjVqlBreu3TpUrehkYWFher2bdu2qWGTGHK4Y8cOx3//+19Hhw4dHP369as29HDQoEFq+CKGEzZt2tQyw7bvuOMO1T7s/5dffqmGSGJoJEYd6KGVGG75ySefqHb27dtXXaKlfXpUINqAERxm0Xr8jh07poax4oKPrhkzZqh/61FKGJaOvy+0Z/369WoEjLdh6T179nSsXLnS8cUXXzhOOOEEt2HNGGWCIb/XXnutGnqLv2EMjw3HkN+a2ldaWqqG2bdq1UodD/PfpR7Zsnz5cjW6B7djmPPrr7+ujtl1111n+fbhtjvvvFON5MF78uOPP3aceuqp6vgUFxdHxfGrrY3mYeXYJ4xM8mT1YziqlnNDsD479bD0cePGqVFes2bN4rD0aID5A3DwMR8Phqlj/girwx+qtwvmX4Bdu3apk2OjRo1UQIe5MPDGNM/jAjt37nQMGTJEzRGBYAJBRllZmcMKMMSxRYsW6rgcd9xx6ncEAhpOkjfffLMa/ok/vEsuuUT9YUdL++Cjjz5Sx23Lli1u10fr8cMcQt7elxjOrIemP/DAA+pkgHadf/751dp+8OBBdYLMzMxUw2BvuOEGdZIywxw+Z599tnoMvDcQSEW6fQgCfP1d6rmVVq9e7ejTp486IaWmpjq6dOnieOSRR9wCBqu2DydMnABx4sOwZgzNxjxRnl8OrXz8amujhsAEf1MIXDxZ/RhKLeeGYH524rXs0aOH+ozGFzLzc9QVV0snIiIi22MNDxEREdkeAx4iIiKyPQY8REREZHsMeIiIiMj2GPAQERGR7THgISIiIttjwENERES2x4CHiIiIbI8BDxHZyqRJk6RHjx6R3g0ishgGPEREdRQXFyfz5s3j60cUBRjwEBERke0x4CGiOuvfv7+MGTNGXbKzs6VJkybywAMPiHOdQZFnn31WTjjhBElNTZXmzZvLH//4R+O+lZWVMnXqVGnfvr2kpaVJ9+7d5V//+pdx+5w5c6Rhw4Zuz4dsCrIqZo8++qh67AYNGsiIESOkuLjY7XY8z5QpU6RVq1aSkpKiursWLlzoV/tKS0tV21q0aKHa0LZtW7XP0K5dO/XzkksuUfukf4f//ve/cuqpp6r7dOjQQSZPnizl5eXG7dj+ueeekyFDhqi2Yxtz24koBOq9/CgRxazf/e53auXqW2+91bF582bH66+/rlZIfuGFFxxff/21IyEhwfHmm2+q1ZG//fZbx8yZM437PvTQQ47OnTs7Fi5c6Ni+fbtaDRmrPy9dulTdjt+xarTZe++9p1Zn1ubOnavu89JLL6nnv++++xwNGjRwdO/e3dhmxowZamXtt956S21z1113qRW5f/jhh1rbN336dEfr1q0dn3/+uWrDsmXLVHtg//79xkrRWA0avwO2xfPNmTNHtWvRokWOdu3aOSZNmmQ8Lu7XuHFjx4svvqhWdL///vvVa/X999/X63gQkW8MeIioXgFPly5dHJWVlcZ1d999t7ru3//+tzrx5+XlVbtfcXGxCoyWL1/udv2IESMcV111ld8BT9++fR0333yz2zZ9+vRxC3hatmzpePjhh922Oe2006rdz5u//e1vjvPOO8+tfWbYF+yT2fnnn+945JFH3K775z//6WjRooXb/W666aZq+z1q1Kha94mI6oZdWkRUL2eccYZbN1Pfvn1l69atcv7556suIHTXXHvttfLGG29IYWGh2mbbtm3q3wMHDpTMzEzj8tprr8n27dv9fu5NmzZJnz593K7D82t5eXny66+/yllnneW2DX7HfWtz/fXXy9q1a6VTp05yyy23yKJFi2q9z7p161QXmrldI0eOlD179hjt99xP/bs/+0REdZNYx/sREdUIJ/pvv/1Wli5dqgKFCRMmqCHjX3/9teTn56tt5s+fL8cdd5zb/VBnA/Hx8UYtkFZWVhbWVx11ODt27JAPP/xQPv74Y/nTn/4kAwYMqLHeBm1Dzc6ll15a7TbU9BBRZDDDQ0T1snLlSrffv/rqK1WonJCQIImJiSpAmDZtmqxfv1527twpn3zyiZx00kkqsNm1a5d07NjR7dK6dWv1OE2bNpVjx45JQUGB8djItph16dLF6/NrWVlZ0rJlS/nyyy/dtsHv2Ad/4DGuuOIKefHFF2Xu3Lny73//Ww4dOqRuS0pKkoqKimpB0pYtW6q1CxcEcd72U/+O9hBRaDDDQ0T1gqBl7Nix8te//lVldJ5++ml54okn5IMPPpAff/xR+vXrJzk5ObJgwQI1YgrdQxhRdeedd8rtt9+urjv77LPl6NGjKhBBgDF8+HDVVZWeni733nuv6k5CYIORW2a33nqr6nbq3bu36qZCt9nGjRtVN5o2btw4mThxohx//PFqhNYrr7yiAidsW5sZM2aoEVo9e/ZUwcq7774rubm5xugxjMxasmSJem4EcGgnMlm///3vpU2bNmpUGu6Hbq4NGzbIQw89ZDw2Hgv7jbZjX1atWiUvv/wy341EoVLH2h8iIlW0jOJfFOCiQDknJ8dx7733qiJfjGjC7bguLS3N0a1bNzWqSsM2Tz31lKNTp05q1FTTpk0dgwcPdnz22WfGNigI7tixo7r/73//ezX6y/NjCwXJTZo0UaPFhg8frkZhmYuWKyoq1Aip4447Tj0Pbvvwww/9Onp4vh49ejgyMjJU+1CQjNFm2vvvv6/2LzEx0dG2bVvjeow8O/PMM9V+436nn366eiwNbZg1a5Zj4MCBapQZRnGZXxsiCr44/C9k0RQR2X4eHmRNnnrqqUjvSlRBkfd7770nw4YNi/SuEMUM1vAQERGR7THgIaKY9cgjj7gNHzdfMAsyEdkHu7SIKGZhtJUeceUJSz54DpknoujFgIeIiIhsj11aREREZHsMeIiIiMj2GPAQERGR7THgISIiIttjwENERES2x4CHiIiIbI8BDxEREdkeAx4iIiISu/t/bgXF1kJmuEoAAAAASUVORK5CYII=",
"text/plain": [
""
]
},
"metadata": {}
}
],
"id": "10095f31"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 7. Report generation\n",
"This section illustrates the process for generating summary reports of simulation results with Flow360's reporting tool. Automated reports simplify compiling, visualizing, and sharing key simulation outputs, including tables, charts, and calculated quantities, without the need for manual post-processing. The example below builds and customizes a report that includes convergence history, force coefficients, and hinge torque summaries.\n",
"\n",
"The object necessary to construct a `ReportTemplate` are imported from `flow360.plugins.report` namespace."
],
"id": "7287eb64"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"from flow360.plugins.report import (\n",
" ReportTemplate,\n",
" Chart2D,\n",
" Table,\n",
" NonlinearResiduals,\n",
" DataItem,\n",
" Average,\n",
")"
],
"execution_count": 14,
"outputs": [],
"id": "b3b7f1a1"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"As the first entries to the report summary tables are created containing force coefficients of the whole aircraft and monitored hinge torques. \n",
"An object representing data is the `DataItem` that is defined using the `data` argument which is a \"path\" to a certain variable as would be accessed through the `Case` object. \n",
"For example to access angle of attack:\n",
"- through `Case`: `case.params.operating_condition.alpha`,\n",
"- path for `DataItem`: `\"params/operating_condition/alpha\"`.\n",
"\n",
"`DataItem` also allows to perform operations on the base data provided. In this case an `Average` operation is performed to get the variables averaged over the last 10% of pseudo steps. \n",
"\n",
"In the cell below necessary `DataItem`s are defined and `Table`s are constructed. "
],
"id": "8d5dedd7"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"left_aileron_hinge_torque_final = DataItem(\n",
" data=\"results/monitors/left_aileron_hinge_torque/left_aileron_hinge_torque_partial_integral\",\n",
" operations=[Average(fraction=0.1)],\n",
")\n",
"right_aileron_hinge_torque_final = DataItem(\n",
" data=\"results/monitors/right_aileron_hinge_torque/right_aileron_hinge_torque_partial_integral\",\n",
" operations=[Average(fraction=0.1)],\n",
")\n",
"left_rudder_hinge_torque_final = DataItem(\n",
" data=\"results/monitors/left_rudder_hinge_torque/left_rudder_hinge_torque_partial_integral\",\n",
" operations=[Average(fraction=0.1)],\n",
")\n",
"right_rudder_hinge_torque_final = DataItem(\n",
" data=\"results/monitors/right_rudder_hinge_torque/right_rudder_hinge_torque_partial_integral\",\n",
" operations=[Average(fraction=0.1)],\n",
")\n",
"\n",
"CL = DataItem(data=\"results/total_forces/CL\", operations=[Average(fraction=0.1)])\n",
"CD = DataItem(data=\"results/total_forces/CD\", operations=[Average(fraction=0.1)])\n",
"\n",
"forces_table = Table(data=[CL, CD], section_title=\"Forces\")\n",
"aileron_torques_table = Table(\n",
" data=[\n",
" left_aileron_hinge_torque_final,\n",
" right_aileron_hinge_torque_final,\n",
" ],\n",
" section_title=\"Aileron hinge torques\",\n",
")\n",
"rudder_torques_table = Table(\n",
" data=[left_rudder_hinge_torque_final, right_rudder_hinge_torque_final],\n",
" section_title=\"Rudder hinge torques\",\n",
")"
],
"execution_count": 15,
"outputs": [],
"id": "1dfe8404"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Another useful object to be used is `Chart2D` that allows to visualize time histories or pseudo step progressions of variables through the simulations. \n",
"A different set of `DataItem`s has to be defined not to include the `Average` operation."
],
"id": "449a2632"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"left_aileron_hinge_torque = DataItem(\n",
" data=\"results/monitors/left_aileron_hinge_torque/left_aileron_hinge_torque_partial_integral\",\n",
" title=\"left_aileron_hinge_torque_monitor\",\n",
")\n",
"right_aileron_hinge_torque = DataItem(\n",
" data=\"results/monitors/right_aileron_hinge_torque/right_aileron_hinge_torque_partial_integral\",\n",
" title=\"right_aileron_hinge_torque_monitor\",\n",
")\n",
"left_rudder_hinge_torque = DataItem(\n",
" data=\"results/monitors/left_rudder_hinge_torque/left_rudder_hinge_torque_partial_integral\",\n",
" title=\"left_rudder_hinge_torque_monitor\",\n",
")\n",
"right_rudder_hinge_torque = DataItem(\n",
" data=\"results/monitors/right_rudder_hinge_torque/right_rudder_hinge_torque_partial_integral\",\n",
" title=\"right_rudder_hinge_torque_monitor\",\n",
")\n",
"\n",
"aileron_hinge_torques_convergence = Chart2D(\n",
" x=\"results/monitors/left_aileron_hinge_torque/pseudo_step\",\n",
" y=[\n",
" left_aileron_hinge_torque,\n",
" right_aileron_hinge_torque,\n",
" ],\n",
" section_title=\"Aileron hinge torques convergence\",\n",
" show_grid=True,\n",
" ylim=[-100, 100],\n",
")\n",
"\n",
"rudder_hinge_torques_convergence = Chart2D(\n",
" x=\"results/monitors/left_rudder_hinge_torque/pseudo_step\",\n",
" y=[left_rudder_hinge_torque, right_rudder_hinge_torque],\n",
" section_title=\"Rudder hinge torques convergence\",\n",
" show_grid=True,\n",
" ylim=[-30, 30],\n",
")"
],
"execution_count": 16,
"outputs": [],
"id": "bb623054"
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Finally all of the objects are compiled into one `ReportTemplate` and submitted to the cloud. An additional `NonlinearResiduals` object is added that automatically adds the plots of convergence to the report.\n",
"\n",
"When submitting the report using `create_in_cloud` the list of `cases` has to be given. In this example only one is specified but when performing multiple simulations more cases can be specified and compared. The tables would then be populated with the data from all of the specified cases.\n",
"\n",
"For an example of utilizing the report API for more `Case` instances see [Alpha sweep tutorial](./alpha_sweep.ipynb) or [Wall roughness tutorial](./wall_roughness.ipynb)"
],
"id": "16f18174"
},
{
"cell_type": "code",
"metadata": {},
"source": [
"template = ReportTemplate(\n",
" title=\"Hinge forces report\",\n",
" items=[\n",
" forces_table,\n",
" aileron_torques_table,\n",
" rudder_torques_table,\n",
" NonlinearResiduals(),\n",
" aileron_hinge_torques_convergence,\n",
" rudder_hinge_torques_convergence,\n",
" ],\n",
")\n",
"\n",
"report = template.create_in_cloud(name=\"hinge_forces\", cases=[case])\n",
"\n",
"report.wait()\n",
"report.download(\"report.pdf\")"
],
"execution_count": 17,
"outputs": [
{
"output_type": "display_data",
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "d5802357730e421f99d09e9c56ca9258",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"Output()"
]
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/html": [
"\n"
],
"text/plain": []
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "fecf693752cf429189ca0a9e37341bd2",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"Output()"
]
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/html": [
"\n"
],
"text/plain": []
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/html": [
"