From d4b9077d339bf8574099f07e4b45f3470d78bc96 Mon Sep 17 00:00:00 2001
From: Ben <106089368+benflexcompute@users.noreply.github.com>
Date: Tue, 13 Jan 2026 15:00:03 -0500
Subject: [PATCH] =?UTF-8?q?Revert=20"refactor():=20Moved=20implementation?=
=?UTF-8?q?=20of=20the=20report=20functionality=20out=20of=20F=E2=80=A6"?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
This reverts commit d103a560d9731d1357d2ee80097dfaa17c4f86e2.
---
.../report/automotive_report.py | 424 +
.../report/geometry_to_report.py | 194 +
.../report/sweep_launch_report.py | 336 +
.../report/sweep_launch_template.py | 289 +
.../special_features/bet_disk.py | 35 +-
.../special_features/heat_transfer.py | 69 +-
.../notebooks/notebook_automotive.ipynb | 333 +-
flow360/__init__.py | 2 +
flow360/cloud/flow360_requests.py | 14 +
flow360/component/project.py | 7 +
flow360/component/simulation/services.py | 14 +
flow360/plugins/__init__.py | 22 -
flow360/plugins/report/README.md | 311 +
flow360/plugins/report/__init__.py | 69 +
.../plugins/report/fonts/TWKEverett-Bold.otf | Bin 0 -> 178056 bytes
.../report/fonts/TWKEverett-BoldItalic.otf | Bin 0 -> 182960 bytes
.../report/fonts/TWKEverett-Regular.otf | Bin 0 -> 150540 bytes
.../report/fonts/TWKEverett-RegularItalic.otf | Bin 0 -> 155060 bytes
flow360/plugins/report/img/FC_logo.pdf | Bin 0 -> 2127 bytes
flow360/plugins/report/img/Flow360-logo.svg | 157 +
flow360/plugins/report/img/Flow360_logo.pdf | Bin 0 -> 16430 bytes
flow360/plugins/report/img/background.pdf | Bin 0 -> 94114 bytes
.../img/colorbar_rainbow_banded_25.jpeg | Bin 0 -> 1303 bytes
.../report/img/colorbar_rainbow_banded_30.png | Bin 0 -> 100 bytes
.../report/img/colorbar_rainbow_smooth.png | Bin 0 -> 99 bytes
flow360/plugins/report/img/cover_logo.pdf | Bin 0 -> 17283 bytes
flow360/plugins/report/img/flow360.png | Bin 0 -> 7245 bytes
.../plugins/report/img/flow360_logo_grey.pdf | Bin 0 -> 3522 bytes
flow360/plugins/report/report.py | 365 +
flow360/plugins/report/report_context.py | 31 +
flow360/plugins/report/report_doc.py | 247 +
flow360/plugins/report/report_items.py | 2184 ++++
flow360/plugins/report/utils.py | 1216 ++
flow360/plugins/report/uvf_shutter.py | 964 ++
poetry.lock | 58 +-
pyproject.toml | 1 +
.../results/nonlinear_residual_v2.csv | 202 +
.../results/surface_forces_v2.csv | 202 +
.../results/total_forces_v2.csv | 202 +
.../simulation.json | 763 ++
.../results/nonlinear_residual_v2.csv | 10001 ++++++++++++++++
.../results/surface_forces_v2.csv | 10001 ++++++++++++++++
.../results/surface_heat_transfer_v2.csv | 4001 +++++++
.../results/total_acoustics_v3.csv | 2001 ++++
.../results/total_forces_v2.csv | 10001 ++++++++++++++++
.../simulation.json | 474 +
.../results/nonlinear_residual_v2.csv | 10001 ++++++++++++++++
.../results/surface_forces_v2.csv | 10001 ++++++++++++++++
.../results/total_forces_v2.csv | 10001 ++++++++++++++++
.../simulation.json | 475 +
.../results/monitor_massFluxExhaust_v2.csv | 202 +
.../results/monitor_massFluxIntake_v2.csv | 202 +
.../udd_massInflowController_Exhaust_v2.csv | 202 +
.../udd_massOutflowController_Intake_v2.csv | 202 +
.../simulation.json | 1391 +++
.../meshBoundaryBoundingBox.json | 98 +
.../meshStats.json | 10 +
.../simulation.json | 371 +
.../meshBoundaryBoundingBox.json | 34 +
.../meshStats.json | 10 +
.../simulation.json | 342 +
.../simulation.json | 1391 +++
tests/report/report_testing_fixtures.py | 397 +
tests/report/test_report.py | 129 +
tests/report/test_report_items.py | 1869 +++
tests/report/test_report_utils.py | 112 +
tests/report/test_shutter.py | 520 +
.../service/ref/default_report_config.json | 112 +
tests/simulation/service/test_services_v2.py | 7 +
tests/utils.py | 19 +
70 files changed, 83255 insertions(+), 33 deletions(-)
create mode 100644 examples/post_processing/report/automotive_report.py
create mode 100644 examples/post_processing/report/geometry_to_report.py
create mode 100644 examples/post_processing/report/sweep_launch_report.py
create mode 100644 examples/post_processing/report/sweep_launch_template.py
delete mode 100644 flow360/plugins/__init__.py
create mode 100644 flow360/plugins/report/README.md
create mode 100644 flow360/plugins/report/__init__.py
create mode 100644 flow360/plugins/report/fonts/TWKEverett-Bold.otf
create mode 100644 flow360/plugins/report/fonts/TWKEverett-BoldItalic.otf
create mode 100644 flow360/plugins/report/fonts/TWKEverett-Regular.otf
create mode 100644 flow360/plugins/report/fonts/TWKEverett-RegularItalic.otf
create mode 100644 flow360/plugins/report/img/FC_logo.pdf
create mode 100644 flow360/plugins/report/img/Flow360-logo.svg
create mode 100644 flow360/plugins/report/img/Flow360_logo.pdf
create mode 100644 flow360/plugins/report/img/background.pdf
create mode 100644 flow360/plugins/report/img/colorbar_rainbow_banded_25.jpeg
create mode 100644 flow360/plugins/report/img/colorbar_rainbow_banded_30.png
create mode 100644 flow360/plugins/report/img/colorbar_rainbow_smooth.png
create mode 100644 flow360/plugins/report/img/cover_logo.pdf
create mode 100644 flow360/plugins/report/img/flow360.png
create mode 100644 flow360/plugins/report/img/flow360_logo_grey.pdf
create mode 100644 flow360/plugins/report/report.py
create mode 100644 flow360/plugins/report/report_context.py
create mode 100644 flow360/plugins/report/report_doc.py
create mode 100644 flow360/plugins/report/report_items.py
create mode 100644 flow360/plugins/report/utils.py
create mode 100644 flow360/plugins/report/uvf_shutter.py
create mode 100644 tests/data/case-333333333-333333-3333333333-33333333/results/nonlinear_residual_v2.csv
create mode 100644 tests/data/case-333333333-333333-3333333333-33333333/results/surface_forces_v2.csv
create mode 100644 tests/data/case-333333333-333333-3333333333-33333333/results/total_forces_v2.csv
create mode 100644 tests/data/case-333333333-333333-3333333333-33333333/simulation.json
create mode 100644 tests/data/case-444444444-444444-4444444444-44444444/results/nonlinear_residual_v2.csv
create mode 100644 tests/data/case-444444444-444444-4444444444-44444444/results/surface_forces_v2.csv
create mode 100644 tests/data/case-444444444-444444-4444444444-44444444/results/surface_heat_transfer_v2.csv
create mode 100644 tests/data/case-444444444-444444-4444444444-44444444/results/total_acoustics_v3.csv
create mode 100644 tests/data/case-444444444-444444-4444444444-44444444/results/total_forces_v2.csv
create mode 100644 tests/data/case-444444444-444444-4444444444-44444444/simulation.json
create mode 100644 tests/data/case-5555-5555555-5555555555-555555555555/results/nonlinear_residual_v2.csv
create mode 100644 tests/data/case-5555-5555555-5555555555-555555555555/results/surface_forces_v2.csv
create mode 100644 tests/data/case-5555-5555555-5555555555-555555555555/results/total_forces_v2.csv
create mode 100644 tests/data/case-5555-5555555-5555555555-555555555555/simulation.json
create mode 100644 tests/data/case-666666666-66666666-666-6666666666666/results/monitor_massFluxExhaust_v2.csv
create mode 100644 tests/data/case-666666666-66666666-666-6666666666666/results/monitor_massFluxIntake_v2.csv
create mode 100644 tests/data/case-666666666-66666666-666-6666666666666/results/udd_massInflowController_Exhaust_v2.csv
create mode 100644 tests/data/case-666666666-66666666-666-6666666666666/results/udd_massOutflowController_Intake_v2.csv
create mode 100644 tests/data/case-666666666-66666666-666-6666666666666/simulation.json
create mode 100644 tests/data/vm-11111111-1111-1111-1111-111111111111/meshBoundaryBoundingBox.json
create mode 100644 tests/data/vm-11111111-1111-1111-1111-111111111111/meshStats.json
create mode 100644 tests/data/vm-11111111-1111-1111-1111-111111111111/simulation.json
create mode 100644 tests/data/vm-22222222-22222222-2222-2222-22222222/meshBoundaryBoundingBox.json
create mode 100644 tests/data/vm-22222222-22222222-2222-2222-22222222/meshStats.json
create mode 100644 tests/data/vm-22222222-22222222-2222-2222-22222222/simulation.json
create mode 100644 tests/data/vm-33333333-33333-3333333-333333333333/simulation.json
create mode 100644 tests/report/report_testing_fixtures.py
create mode 100644 tests/report/test_report.py
create mode 100644 tests/report/test_report_items.py
create mode 100644 tests/report/test_report_utils.py
create mode 100644 tests/report/test_shutter.py
create mode 100644 tests/simulation/service/ref/default_report_config.json
diff --git a/examples/post_processing/report/automotive_report.py b/examples/post_processing/report/automotive_report.py
new file mode 100644
index 000000000..f95b66149
--- /dev/null
+++ b/examples/post_processing/report/automotive_report.py
@@ -0,0 +1,424 @@
+import flow360 as fl
+from flow360 import u
+from flow360.examples import DrivAer
+from flow360.log import log
+from flow360.plugins.report.report import ReportTemplate
+from flow360.plugins.report.report_items import (
+ BottomCamera,
+ Chart2D,
+ Chart3D,
+ FrontCamera,
+ FrontLeftBottomCamera,
+ FrontLeftTopCamera,
+ Inputs,
+ LeftCamera,
+ RearCamera,
+ RearLeftTopCamera,
+ RearRightBottomCamera,
+ Settings,
+ Summary,
+ Table,
+ TopCamera,
+)
+from flow360.plugins.report.utils import Average, DataItem, Delta, Expression, Variable
+
+DrivAer.get_files()
+
+project = fl.Project.from_volume_mesh(
+ DrivAer.mesh_filename,
+ name="Automotive DrivAer",
+)
+
+vm = project.volume_mesh
+
+log.info("Volume mesh contains the following boundaries:")
+for boundary in vm.boundary_names:
+ log.info("Boundary: " + boundary)
+
+freestream_surfaces = ["blk-1/WT_side1", "blk-1/WT_side2", "blk-1/WT_inlet", "blk-1/WT_outlet"]
+slip_wall_surfaces = ["blk-1/WT_ceiling", "blk-1/WT_ground_front", "blk-1/WT_ground"]
+wall_surfaces = list(set(vm.boundary_names) - set(freestream_surfaces) - set(slip_wall_surfaces))
+
+cases = []
+
+for beta in [0, 5, 10]:
+ with fl.SI_unit_system:
+ params = fl.SimulationParams(
+ meshing=None,
+ reference_geometry=fl.ReferenceGeometry(area=2.17, moment_length=2.7862),
+ operating_condition=fl.AerospaceCondition(velocity_magnitude=40, beta=beta * u.deg),
+ models=[
+ fl.Wall(surfaces=[vm[i] for i in wall_surfaces], use_wall_function=True),
+ fl.Freestream(
+ surfaces=[vm[i] for i in freestream_surfaces],
+ ),
+ fl.SlipWall(
+ surfaces=[vm[i] for i in slip_wall_surfaces],
+ ),
+ ],
+ user_defined_fields=[
+ fl.UserDefinedField(
+ name="Cpx",
+ expression="double prel = primitiveVars[4] - pressureFreestream;"
+ + "double PressureForce_X = prel * nodeNormals[0]; "
+ + "Cpx = PressureForce_X / (0.5 * MachRef * MachRef) / magnitude(nodeNormals);",
+ ),
+ ],
+ outputs=[
+ fl.SurfaceOutput(
+ surfaces=vm["*"],
+ output_fields=[
+ "Cp",
+ "Cf",
+ "yPlus",
+ "CfVec",
+ "primitiveVars",
+ "wall_shear_stress_magnitude",
+ "Cpx",
+ ],
+ ),
+ fl.SliceOutput(
+ entities=[
+ *[
+ fl.Slice(
+ name=f"slice_y_{name}",
+ normal=(0, 1, 0),
+ origin=(0, y, 0),
+ )
+ for name, y in zip(
+ ["0", "0_2", "0_4", "0_6", "0_8"], [0, 0.2, 0.4, 0.6, 0.8]
+ )
+ ],
+ *[
+ fl.Slice(
+ name=f"slice_z_{name}",
+ normal=(0, 0, 1),
+ origin=(0, 0, z),
+ )
+ for name, z in zip(
+ ["neg0_2", "0", "0_2", "0_4", "0_6", "0_8"],
+ [-0.2, 0, 0.2, 0.4, 0.6, 0.8],
+ )
+ ],
+ ],
+ output_fields=["velocity", "velocity_x", "velocity_y", "velocity_z"],
+ ),
+ fl.IsosurfaceOutput(
+ output_fields=["Cp", "Mach"],
+ isosurfaces=[
+ fl.Isosurface(
+ name="isosurface-cpt",
+ iso_value=-1,
+ field="Cpt",
+ ),
+ ],
+ ),
+ fl.ProbeOutput(
+ entities=[fl.Point(name="point1", location=(10, 0, 1))],
+ output_fields=["velocity"],
+ ),
+ ],
+ )
+
+ case_new = project.run_case(params=params, name=f"DrivAer 5.7M - beta={beta}")
+
+ cases.append(case_new)
+
+# wait until all cases finish running
+for case in cases:
+ case.wait()
+
+exclude = ["blk-1/WT_ground_close", "blk-1/WT_ground_patch"]
+size = "5.7M"
+
+exclude += freestream_surfaces + slip_wall_surfaces
+
+top_camera = TopCamera(pan_target=(1.5, 0, 0), dimension=5, dimension_dir="width")
+top_camera_slice = TopCamera(pan_target=(2.5, 0, 0), dimension=8, dimension_dir="width")
+side_camera = LeftCamera(pan_target=(1.5, 0, 0), dimension=5, dimension_dir="width")
+side_camera_slice = LeftCamera(pan_target=(2.5, 0, 1.5), dimension=8, dimension_dir="width")
+rear_camera = RearCamera(dimension=2.5, dimension_dir="width")
+front_camera = FrontCamera(dimension=2.5, dimension_dir="width")
+bottom_camera = BottomCamera(pan_target=(1.5, 0, 0), dimension=5, dimension_dir="width")
+front_left_bottom_camera = FrontLeftBottomCamera(
+ pan_target=(1.5, 0, 0), dimension=5, dimension_dir="width"
+)
+rear_right_bottom_camera = RearRightBottomCamera(
+ pan_target=(1.5, 0, 0), dimension=6, dimension_dir="width"
+)
+front_left_top_camera = FrontLeftTopCamera(
+ pan_target=(1.5, 0, 0), dimension=6, dimension_dir="width"
+)
+rear_left_top_camera = RearLeftTopCamera(pan_target=(1.5, 0, 0), dimension=6, dimension_dir="width")
+
+cameras_geo = [
+ top_camera,
+ side_camera,
+ rear_camera,
+ bottom_camera,
+ front_left_bottom_camera,
+ rear_right_bottom_camera,
+]
+
+limits_cp = [(-1, 1), (-1, 1), (-1, 1), (-0.3, 0), (-0.3, 0), (-1, 1), (-1, 1), (-1, 1)]
+cameras_cp = [
+ front_camera,
+ front_left_top_camera,
+ side_camera,
+ rear_left_top_camera,
+ rear_camera,
+ bottom_camera,
+ front_left_bottom_camera,
+ rear_right_bottom_camera,
+]
+
+avg = Average(fraction=0.1)
+CD = DataItem(data="surface_forces/totalCD", exclude=exclude, title="CD", operations=avg)
+
+CL = DataItem(data="surface_forces/totalCL", exclude=exclude, title="CL", operations=avg)
+
+CDA = DataItem(
+ data="surface_forces",
+ exclude=exclude,
+ title="CD*area",
+ variables=[Variable(name="area", data="params.reference_geometry.area")],
+ operations=[Expression(expr="totalCD * area"), avg],
+)
+
+CLf = DataItem(
+ data="surface_forces",
+ exclude=exclude,
+ title="CLf",
+ operations=[Expression(expr="1/2*totalCL + totalCMy"), avg],
+)
+
+CLr = DataItem(
+ data="surface_forces",
+ exclude=exclude,
+ title="CLr",
+ operations=[Expression(expr="1/2*totalCL - totalCMy"), avg],
+)
+
+CFy = DataItem(data="surface_forces/totalCFy", exclude=exclude, title="CS", operations=avg)
+
+statistical_data = [
+ "params/reference_geometry/area",
+ CD,
+ CDA,
+ Delta(data=CD),
+ CL,
+ CLf,
+ CLr,
+ CFy,
+ "volume_mesh/stats/n_nodes",
+ "params/time_stepping/max_steps",
+]
+statistical_table = Table(
+ data=statistical_data,
+ section_title="Statistical data",
+ formatter=[
+ (
+ None
+ if d
+ in [
+ "params/reference_geometry/area",
+ "volume_mesh/stats/n_nodes",
+ "params/time_stepping/max_steps",
+ ]
+ else ".4f"
+ )
+ for d in statistical_data
+ ],
+)
+
+geometry_screenshots = [
+ Chart3D(
+ section_title="Geometry",
+ items_in_row=2,
+ force_new_page=True,
+ show="boundaries",
+ camera=camera,
+ exclude=exclude,
+ fig_name=f"geo_{i}",
+ )
+ for i, camera in enumerate(cameras_geo)
+]
+cpt_screenshots = [
+ Chart3D(
+ section_title="Isosurface, Cpt=-1",
+ items_in_row=2,
+ force_new_page=True,
+ show="isosurface",
+ iso_field="Cpt",
+ exclude=exclude,
+ camera=camera,
+ )
+ for camera in cameras_cp
+]
+cfvec_screenshots = [
+ Chart3D(
+ section_title="CfVec",
+ items_in_row=2,
+ force_new_page=True,
+ show="boundaries",
+ field="CfVec",
+ mode="lic",
+ limits=(1e-4, 10),
+ is_log_scale=True,
+ exclude=exclude,
+ camera=camera,
+ )
+ for camera in cameras_cp
+]
+y_slices_lic_screenshots = [
+ Chart3D(
+ section_title=f"Slice velocity LIC y={y}",
+ items_in_row=2,
+ force_new_page=True,
+ show="slices",
+ include=[f"slice_y_{name}"],
+ field="velocityVec",
+ mode="lic",
+ limits=(0 * u.m / u.s, 50 * u.m / u.s),
+ camera=side_camera_slice,
+ fig_name=f"slice_y_vec_{name}",
+ )
+ for name, y in zip(["0", "0_2", "0_4", "0_6", "0_8"], [0, 0.2, 0.4, 0.6, 0.8])
+]
+y_slices_screenshots = [
+ Chart3D(
+ section_title=f"Slice velocity y={y}",
+ items_in_row=2,
+ force_new_page=True,
+ show="slices",
+ include=[f"slice_y_{name}"],
+ field="velocity",
+ limits=(0 * u.m / u.s, 50 * u.m / u.s),
+ camera=side_camera_slice,
+ fig_name=f"slice_y_{name}",
+ )
+ for name, y in zip(["0", "0_2", "0_4", "0_6", "0_8"], [0, 0.2, 0.4, 0.6, 0.8])
+]
+y_slices_lic_screenshots = [
+ Chart3D(
+ section_title=f"Slice velocity LIC y={y}",
+ items_in_row=2,
+ force_new_page=True,
+ show="slices",
+ include=[f"slice_y_{name}"],
+ field="velocityVec",
+ mode="lic",
+ limits=(0 * u.m / u.s, 50 * u.m / u.s),
+ camera=side_camera_slice,
+ fig_name=f"slice_y_vec_{name}",
+ )
+ for name, y in zip(["0", "0_2", "0_4", "0_6", "0_8"], [0, 0.2, 0.4, 0.6, 0.8])
+]
+z_slices_screenshots = [
+ Chart3D(
+ section_title=f"Slice velocity z={z}",
+ items_in_row=2,
+ force_new_page=True,
+ show="slices",
+ include=[f"slice_z_{name}"],
+ field="velocity",
+ limits=(0 * u.m / u.s, 50 * u.m / u.s),
+ camera=top_camera_slice,
+ fig_name=f"slice_z_{name}",
+ )
+ for name, z in zip(["neg0_2", "0", "0_2", "0_4", "0_6", "0_8"], [-0.2, 0, 0.2, 0.4, 0.6, 0.8])
+]
+y_plus_screenshots = [
+ Chart3D(
+ section_title="y+",
+ items_in_row=2,
+ show="boundaries",
+ field="yPlus",
+ exclude=exclude,
+ limits=(0, 5),
+ camera=camera,
+ fig_name=f"yplus_{i}",
+ )
+ for i, camera in enumerate([top_camera, bottom_camera])
+]
+cp_screenshots = [
+ Chart3D(
+ section_title="Cp",
+ items_in_row=2,
+ show="boundaries",
+ field="Cp",
+ exclude=exclude,
+ limits=limits,
+ camera=camera,
+ fig_name=f"cp_{i}",
+ )
+ for i, (limits, camera) in enumerate(zip(limits_cp, cameras_cp))
+]
+cpx_screenshots = [
+ Chart3D(
+ section_title="Cpx",
+ items_in_row=2,
+ show="boundaries",
+ field="Cpx",
+ exclude=exclude,
+ limits=(-0.3, 0.3),
+ camera=camera,
+ fig_name=f"cpx_{i}",
+ )
+ for i, camera in enumerate(cameras_cp)
+]
+wall_shear_screenshots = [
+ Chart3D(
+ section_title="Wall shear stress magnitude",
+ items_in_row=2,
+ show="boundaries",
+ field="wallShearMag",
+ exclude=exclude,
+ limits=(0 * u.Pa, 5 * u.Pa),
+ camera=camera,
+ fig_name=f"wallShearMag_{i}",
+ )
+ for i, camera in enumerate(cameras_cp)
+]
+
+report = ReportTemplate(
+ title="Aerodynamic analysis of DrivAer",
+ items=[
+ Summary(),
+ Inputs(),
+ statistical_table,
+ Chart2D(
+ x="x_slicing_force_distribution/X",
+ y="x_slicing_force_distribution/totalCumulative_CD_Curve",
+ fig_name="totalCumulative_CD_Curve",
+ background="geometry",
+ exclude=exclude,
+ ),
+ Chart2D(
+ x="surface_forces/pseudo_step",
+ y="surface_forces/totalCD",
+ section_title="Drag Coefficient",
+ fig_name="cd_fig",
+ exclude=exclude,
+ focus_x=(1 / 3, 1),
+ ),
+ *geometry_screenshots,
+ *cp_screenshots,
+ *cpx_screenshots,
+ *cpt_screenshots,
+ *y_slices_screenshots,
+ *y_slices_lic_screenshots,
+ *z_slices_screenshots,
+ *y_plus_screenshots,
+ *wall_shear_screenshots,
+ ],
+ settings=Settings(dpi=150),
+)
+
+report = report.create_in_cloud(
+ f"{size}-{len(cases)}cases-slices-using-groups-Cpt, Cpx, wallShear, dpi=default",
+ cases,
+)
+
+report.wait()
+report.download("report.pdf")
diff --git a/examples/post_processing/report/geometry_to_report.py b/examples/post_processing/report/geometry_to_report.py
new file mode 100644
index 000000000..8176d1d33
--- /dev/null
+++ b/examples/post_processing/report/geometry_to_report.py
@@ -0,0 +1,194 @@
+import flow360 as fl
+from flow360 import u
+from flow360.examples import Airplane
+from flow360.plugins.report.report import ReportTemplate
+from flow360.plugins.report.report_items import (
+ Camera,
+ Chart2D,
+ Chart3D,
+ Inputs,
+ Settings,
+ Summary,
+ Table,
+)
+from flow360.plugins.report.utils import Average, DataItem, Delta
+
+project_id = None # if running for the first time
+
+# then replace it with your project ID to avoid re-creation of projects. You can find project ID on web GUI:
+# project_id = "prj-...."
+
+if project_id is not None:
+ project = fl.Project.from_cloud(project_id)
+else:
+ project = fl.Project.from_geometry(
+ Airplane.geometry, name="Python Project (Geometry, from file) - for Report"
+ )
+
+geo = project.geometry
+geo.group_faces_by_tag("groupName")
+
+
+def simulation_params(angle_of_attack):
+ with fl.SI_unit_system:
+ far_field_zone = fl.AutomatedFarfield()
+ params = fl.SimulationParams(
+ meshing=fl.MeshingParams(
+ defaults=fl.MeshingDefaults(
+ boundary_layer_first_layer_thickness=0.001,
+ surface_max_edge_length=1,
+ ),
+ volume_zones=[far_field_zone],
+ ),
+ reference_geometry=fl.ReferenceGeometry(area=1, moment_length=1),
+ operating_condition=fl.AerospaceCondition(
+ velocity_magnitude=100,
+ alpha=angle_of_attack * fl.u.deg,
+ ),
+ time_stepping=fl.Steady(max_steps=1000),
+ models=[
+ fl.Wall(
+ surfaces=[geo["*"]],
+ ),
+ fl.Freestream(
+ surfaces=[far_field_zone.farfield],
+ ),
+ ],
+ outputs=[
+ fl.SurfaceOutput(
+ surfaces=geo["*"],
+ output_fields=[
+ "Cp",
+ "Cf",
+ "yPlus",
+ "CfVec",
+ "primitiveVars",
+ ],
+ ),
+ ],
+ )
+ return params
+
+
+cases: list[fl.Case] = []
+for alpha in [0, 2, 4]:
+ case = project.run_case(params=simulation_params(alpha), name=f"Case for report, alpha={alpha}")
+ cases.append(case)
+
+[print(case.short_description()) for case in cases]
+
+# waiting explicitly for all the cases to finish (report pipeline will not wait for cases)
+[case.wait() for case in cases]
+
+
+top_camera = Camera(
+ position=(0, 0, 1),
+ look_at=(0, 0, 0),
+ pan_target=(5, 0, 0),
+ up=(0, 1, 0),
+ dimension=15,
+ dimension_dir="width",
+)
+side_camera = Camera(
+ position=(0, -1, 0),
+ look_at=(0, 0, 0),
+ pan_target=(5, 0, 0),
+ up=(0, 0, 1),
+ dimension=12,
+ dimension_dir="width",
+)
+front_left_bottom_camera = Camera(
+ position=(-1, -1, -1),
+ look_at=(0, 0, 0),
+ pan_target=(4, 0, 0),
+ up=(0, 0, 1),
+ dimension=15,
+ dimension_dir="width",
+)
+rear_right_bottom_camera = Camera(
+ position=(1, 1, -1),
+ look_at=(0, 0, 0),
+ pan_target=(4, 0, 0),
+ up=(0, 0, 1),
+ dimension=15,
+ dimension_dir="width",
+)
+
+cameras_geo = [
+ top_camera,
+ side_camera,
+ front_left_bottom_camera,
+ rear_right_bottom_camera,
+]
+
+avg = Average(fraction=0.1)
+
+CL = DataItem(data="surface_forces/totalCL", title="CL", operations=avg)
+
+CD = DataItem(data="surface_forces/totalCD", title="CD", operations=avg)
+
+statistical_data = [
+ "params/operating_condition/alpha",
+ "params/reference_geometry/area",
+ CL,
+ Delta(data=CL),
+ CD,
+ "volume_mesh/stats/n_nodes",
+ "params/time_stepping/max_steps",
+]
+statistical_table = Table(
+ data=statistical_data,
+ section_title="Statistical data",
+ formatter=[
+ (
+ None
+ if d
+ in [
+ "params/reference_geometry/area",
+ "volume_mesh/stats/n_nodes",
+ "params/time_stepping/max_steps",
+ ]
+ else ".4f"
+ )
+ for d in statistical_data
+ ],
+)
+
+geometry_screenshots = [
+ Chart3D(
+ section_title="Geometry",
+ items_in_row=2,
+ force_new_page=True,
+ show="boundaries",
+ camera=camera,
+ fig_name=f"geo_{i}",
+ )
+ for i, camera in enumerate(cameras_geo)
+]
+
+report = ReportTemplate(
+ title="Geometry to report",
+ items=[
+ Summary(),
+ Inputs(),
+ statistical_table,
+ Chart2D(
+ x="total_forces/pseudo_step",
+ y="total_forces/CL",
+ section_title="Lift Coefficient",
+ fig_name="cl_fig",
+ focus_x=(1 / 3, 1),
+ ),
+ *geometry_screenshots,
+ ],
+ settings=Settings(dpi=150),
+)
+
+
+report = report.create_in_cloud(
+ f"Geometry to report - Report, dpi=150",
+ cases,
+)
+
+report.wait()
+report.download("report.pdf")
diff --git a/examples/post_processing/report/sweep_launch_report.py b/examples/post_processing/report/sweep_launch_report.py
new file mode 100644
index 000000000..cc32a1919
--- /dev/null
+++ b/examples/post_processing/report/sweep_launch_report.py
@@ -0,0 +1,336 @@
+"""This script is used by sweep_launch_template.py to create a report."""
+
+import pandas as pd
+
+import flow360 as fl
+from flow360.plugins.report.report import ReportTemplate
+from flow360.plugins.report.report_items import (
+ BottomCamera,
+ Camera,
+ Chart2D,
+ Chart3D,
+ FrontCamera,
+ FrontLeftTopCamera,
+ Inputs,
+ LeftCamera,
+ NonlinearResiduals,
+ PatternCaption,
+ RearCamera,
+ RearRightBottomCamera,
+ Settings,
+ SubsetLimit,
+ Summary,
+ Table,
+ TopCamera,
+)
+from flow360.plugins.report.utils import Average, DataItem
+
+
+def main():
+
+ csv_path = "PATH/TO/CSV/FILE.csv" # Replace with the actual path to your CSV file generated by sweep_launch_template.py
+
+ generate_report(
+ *csv_reader(csv_path),
+ include_geometry=True,
+ include_general_tables=True,
+ include_residuals=True,
+ include_cfl=True,
+ include_forces_moments_table=True,
+ include_forces_moments_charts=True,
+ include_forces_moments_alpha_charts=True,
+ include_forces_moments_beta_charts=True,
+ include_cf_vec=True,
+ include_cp=True,
+ include_yplus=True,
+ include_qcriterion=True,
+ )
+
+
+def csv_reader(
+ file_path,
+):
+ csv_case_specific = pd.read_csv(file_path, sep=",", skiprows=2)
+
+ case_list = []
+
+ for case_id in csv_case_specific["case_id"]:
+ case = fl.Case(case_id)
+ case_list.append(case)
+
+ params = case_list[0].params
+
+ return case_list, params
+
+
+def generate_report(
+ cases,
+ params,
+ include_geometry: bool = False,
+ include_general_tables: bool = False,
+ include_residuals: bool = False,
+ include_cfl: bool = False,
+ include_forces_moments_table: bool = False,
+ include_forces_moments_charts: bool = False,
+ include_forces_moments_alpha_charts: bool = False,
+ include_forces_moments_beta_charts: bool = False,
+ include_cf_vec: bool = False,
+ include_cp: bool = False,
+ include_yplus: bool = False,
+ include_qcriterion: bool = False,
+):
+ items = []
+
+ freestream_surfaces = ["fluid/farfield"]
+
+ exclude = freestream_surfaces
+
+ top_camera = TopCamera(pan_target=(3.5, 0, -0.5), dimension=15, dimension_dir="height")
+ bottom_camera = BottomCamera(pan_target=(3.5, 0, -0.5), dimension=15, dimension_dir="height")
+ front_camera = FrontCamera(pan_target=(3.5, 0, -0.5), dimension=15, dimension_dir="width")
+ rear_camera = RearCamera(pan_target=(3.5, 0, -0.5), dimension=15, dimension_dir="width")
+ left_camera = LeftCamera(pan_target=(3.5, 0, -0.5), dimension=10, dimension_dir="width")
+ right_camera = Camera(
+ pan_target=(3.5, 0, -0.5),
+ position=(0.0, -1.0, 0.0),
+ look_at=(0.0, 0.0, 0.0),
+ up=(0.0, 0.0, 1.0),
+ dimension=10,
+ dimension_dir="width",
+ )
+ front_left_top_camera = FrontLeftTopCamera(
+ pan_target=(3.5, 0, -0.5), dimension=15, dimension_dir="width"
+ )
+ rear_right_bottom_camera = RearRightBottomCamera(
+ pan_target=(3.5, 0, -0.5), dimension=15, dimension_dir="width"
+ )
+
+ if params.time_stepping.type_name == "Unsteady":
+ step_type = "physical_step"
+ else:
+ step_type = "pseudo_step"
+
+ geo_cameras = [
+ top_camera,
+ bottom_camera,
+ front_camera,
+ rear_camera,
+ left_camera,
+ right_camera,
+ front_left_top_camera,
+ rear_right_bottom_camera,
+ ]
+
+ geo_camera_names = [
+ "top_camera",
+ "bottom_camera",
+ "front_camera",
+ "rear_camera",
+ "left_camera",
+ "right_camera",
+ "front_left_top_camera",
+ "rear_right_bottom_camera",
+ ]
+
+ avg = Average(fraction=0.1)
+
+ force_list = [
+ "CD",
+ "CL",
+ "CFx",
+ "CFy",
+ "CFz",
+ "CMx",
+ "CMy",
+ "CMz",
+ ]
+
+ for model in params.models:
+ if model.type == "Fluid":
+ turbulence_solver = model.turbulence_model_solver.type_name
+
+ CD = DataItem(data="surface_forces/totalCD", exclude=exclude, title="CD", operations=avg)
+ CL = DataItem(data="surface_forces/totalCL", exclude=exclude, title="CL", operations=avg)
+ CFX = DataItem(data="surface_forces/totalCFx", exclude=exclude, title="CFx", operations=avg)
+ CFY = DataItem(data="surface_forces/totalCFy", exclude=exclude, title="CFy", operations=avg)
+ CFZ = DataItem(data="surface_forces/totalCFz", exclude=exclude, title="CFz", operations=avg)
+ CMX = DataItem(data="surface_forces/totalCMx", exclude=exclude, title="CMx", operations=avg)
+ CMY = DataItem(data="surface_forces/totalCMy", exclude=exclude, title="CMy", operations=avg)
+ CMZ = DataItem(data="surface_forces/totalCMz", exclude=exclude, title="CMz", operations=avg)
+
+ table_data = [
+ CD,
+ CL,
+ CFX,
+ CFY,
+ CFZ,
+ CMX,
+ CMY,
+ CMZ,
+ ]
+
+ if include_geometry:
+ geometry_screenshots = [
+ Chart3D(
+ section_title="Geometry",
+ items_in_row=2,
+ force_new_page=True,
+ show="boundaries",
+ camera=front_left_top_camera,
+ exclude=exclude,
+ fig_name="Geometry_view",
+ )
+ ]
+ items.extend(geometry_screenshots)
+
+ if include_general_tables:
+ items.append(Summary())
+ items.append(Inputs())
+
+ if include_forces_moments_table:
+ table = Table(
+ data=table_data,
+ section_title="Quantities of interest",
+ )
+ items.append(table)
+
+ if include_residuals:
+ residual_chart = NonlinearResiduals(
+ force_new_page=True, section_title="Nonlinear residuals", fig_name=f"nonlin-res_fig"
+ )
+ items.append(residual_chart)
+
+ if include_cfl and params.time_stepping.CFL.type == "adaptive":
+ cfl_chart = Chart2D(
+ x=f"cfl/{step_type}",
+ y=["cfl/0_NavierStokes_cfl", f"cfl/1_{turbulence_solver}_cfl"],
+ force_new_page=True,
+ section_title="CFL",
+ fig_name="cfl_fig",
+ y_log=True,
+ )
+ items.append(cfl_chart)
+
+ if include_forces_moments_charts:
+ force_charts = [
+ Chart2D(
+ x=f"surface_forces/{step_type}",
+ y=f"surface_forces/total{force}",
+ force_new_page=True,
+ section_title="Forces/Moments",
+ fig_name=f"{force}_fig",
+ exclude=exclude,
+ ylim=SubsetLimit(subset=(0.5, 1), offset=0.25),
+ )
+ for force in force_list
+ ]
+ items.extend(force_charts)
+
+ if include_forces_moments_alpha_charts:
+ force_alpha_charts = [
+ Chart2D(
+ x=f"params/operating_condition/alpha",
+ y=f"total_forces/averages/{force}",
+ force_new_page=True,
+ section_title="Averaged Forces/Moments against alpha",
+ fig_name=f"{force}_alpha_fig",
+ )
+ for force in force_list
+ ]
+ items.extend(force_alpha_charts)
+
+ if include_forces_moments_beta_charts:
+ force_beta_charts = [
+ Chart2D(
+ x=f"params/operating_condition/beta",
+ y=f"total_forces/averages/{force}",
+ force_new_page=True,
+ section_title="Averaged Forces/Moments against beta",
+ fig_name=f"{force}_beta_fig",
+ )
+ for force in force_list
+ ]
+ items.extend(force_beta_charts)
+
+ if include_yplus:
+ y_plus_screenshots = [
+ Chart3D(
+ caption=PatternCaption(pattern=f"y+_{camera_name}_[case.name]"),
+ show="boundaries",
+ field="yPlus",
+ exclude=exclude,
+ limits=(0, 5),
+ camera=camera,
+ fig_name=f"yplus_{camera_name}_fig",
+ fig_size=1,
+ )
+ for camera_name, camera in zip(geo_camera_names, geo_cameras)
+ ]
+ items.extend(y_plus_screenshots)
+
+ if include_cp:
+ cp_screenshots = [
+ Chart3D(
+ caption=PatternCaption(pattern=f"Cp_{camera_name}_[case.name]"),
+ show="boundaries",
+ field="Cp",
+ exclude=exclude,
+ limits=(-1, 1),
+ camera=camera,
+ fig_name=f"cp_{camera_name}_fig",
+ fig_size=1,
+ )
+ for camera_name, camera in zip(geo_camera_names, geo_cameras)
+ ]
+ items.extend(cp_screenshots)
+
+ if include_cf_vec:
+ cfvec_screenshots = [
+ Chart3D(
+ caption=PatternCaption(pattern=f"Cf_vec_{camera_name}_[case.name]"),
+ show="boundaries",
+ field="CfVec",
+ mode="lic",
+ exclude=exclude,
+ limits=(0, 0.025),
+ camera=camera,
+ fig_name=f"cfvec_{camera_name}_fig",
+ fig_size=1,
+ )
+ for camera_name, camera in zip(geo_camera_names, geo_cameras)
+ ]
+ items.extend(cfvec_screenshots)
+
+ if include_qcriterion:
+ qcriterion_screenshots = [
+ Chart3D(
+ caption=PatternCaption(pattern=f"Isosurface_q_criterion_{camera_name}_[case.name]"),
+ show="isosurface",
+ iso_field="qcriterion",
+ exclude=exclude,
+ limits=(0, 0.8),
+ camera=camera,
+ fig_name=f"qcriterion_{camera_name}_fig",
+ fig_size=1,
+ )
+ for camera_name, camera in zip(geo_camera_names, geo_cameras)
+ ]
+ items.extend(qcriterion_screenshots)
+
+ report = ReportTemplate(
+ title="Sweep Template Report",
+ items=items,
+ settings=Settings(dpi=150),
+ )
+
+ report = report.create_in_cloud(
+ "sweep-script-report",
+ cases,
+ )
+
+ report.wait()
+ report.download("report.pdf")
+
+
+if __name__ == "__main__":
+ main()
diff --git a/examples/post_processing/report/sweep_launch_template.py b/examples/post_processing/report/sweep_launch_template.py
new file mode 100644
index 000000000..4c70de5c5
--- /dev/null
+++ b/examples/post_processing/report/sweep_launch_template.py
@@ -0,0 +1,289 @@
+"""
+Sample Flow 360 API scripts.
+Requires a volume mesh or CAD file that you are ready to upload and run cases on.
+
+This script will:
+1:
+ Either:
+ 1.a- create a project and upload a volume mesh
+ 1.b-create a project and upload a CAD. Mesh parameters will need to be given and mesh will creation will need to be launched
+ 1.c- connect to an existing project and associate with a volume mesh from that project.
+
+2:
+ Create the directory structure on the local machine to organize the data that will be generated
+
+3:
+ -Launch all cases required for the sweep
+ -Save all the caseID and other relevant data in the directory created in step 2: above
+"""
+
+import os
+
+import pandas as pd
+from sweep_launch_report import csv_reader, generate_report
+
+import flow360 as fl
+from flow360 import u
+from flow360.examples import EVTOL
+
+# Variables we want to export in our volume and surface solution files. Many more are available.
+vol_fields = ["Mach", "Cp", "mut", "mutRatio", "primitiveVars", "qcriterion"]
+surf_fields = ["Cp", "yPlus", "Cf", "CfVec", "primitiveVars", "wallDistance"]
+
+velocity_magnitude = 100 * u.m / u.s
+
+
+######################################################################################################################
+def make_run_params(mesh_object, models):
+ """
+ Create the params object that contains all the run parameters.
+ Needs the mesh_object to get the list of surfaces.
+ """
+ with fl.SI_unit_system:
+ params = fl.SimulationParams(
+ # Dimensions can be either in inches, m, mm or many other units
+ reference_geometry=fl.ReferenceGeometry(
+ moment_center=(0, 0, 0) * u.m, moment_length=1 * u.m, area=1 * u.m * u.m
+ ),
+ operating_condition=fl.AerospaceCondition(
+ velocity_magnitude=velocity_magnitude, alpha=0 * u.deg
+ ),
+ time_stepping=fl.Steady(max_steps=5000, CFL=fl.AdaptiveCFL()),
+ models=[
+ *models,
+ fl.Fluid(
+ navier_stokes_solver=fl.NavierStokesSolver(),
+ turbulence_model_solver=fl.SpalartAllmaras(),
+ ),
+ ],
+ outputs=[
+ fl.VolumeOutput(output_format="tecplot", output_fields=vol_fields),
+ # mesh_object['*'] will select all the boundaries in the mesh and export the surface results.
+ # Regular expressions can be used to filter for certain boundaries.
+ fl.SurfaceOutput(
+ surfaces=[mesh_object["*"]], output_fields=surf_fields, output_format="tecplot"
+ ),
+ ],
+ )
+
+ # Add meshing params in case the project starts from geometry.
+ params.meshing = mesh_object.params.meshing
+
+ return params
+
+
+######################################################################################################################
+def launch_sweep(params, project, mesh_object, dir_path):
+ """
+ Launch a sweep of cases.
+ """
+
+ case_list = []
+
+ # Create the csv file where we will store all relevant sweep data.
+ os.makedirs(dir_path, exist_ok=True)
+ csv_path = os.path.join(dir_path, "sweep_saved_data.csv")
+ general_info = {
+ "Sweep name": os.path.basename(dir_path),
+ "Project name": project.metadata.name,
+ "Project ID": project.id,
+ "Volume mesh ID": mesh_object.id,
+ "Velocity magnitude": velocity_magnitude,
+ }
+ df = pd.DataFrame.from_dict([general_info])
+ df.to_csv(csv_path, index=False)
+
+ # For example let's vary alpha:
+ alphas = [-10, -5, 0, 5, 10, 12, 14] * u.deg
+
+ cases_params = []
+ for i, alpha_angle in enumerate(alphas):
+ # modify the alpha
+ params.operating_condition.alpha = alpha_angle
+
+ case = project.run_case(params=params, name=f"alpha_{alpha_angle.value}_case")
+ data = {
+ "case_id": case.id,
+ "alpha(deg)": params.operating_condition.alpha.value,
+ "beta(deg)": params.operating_condition.beta.value,
+ }
+
+ print(f"The case ID is: {case.id} with {alpha_angle=} ")
+ case_list.append(case)
+ cases_params.append(data)
+
+ # fraction used for averaging over the last % of iterations, 0.1 -> last 10%
+ fraction = 0.1
+
+ cases_forces = []
+ for i, case in enumerate(case_list):
+ case.wait()
+ results = case.results
+ avg_total_forces = results.total_forces.get_averages(fraction).to_dict()
+ forces = {
+ "CL(avg)": avg_total_forces["CL"],
+ "CD(avg)": avg_total_forces["CD"],
+ "CFx(avg)": avg_total_forces["CFx"],
+ "CFy(avg)": avg_total_forces["CFy"],
+ "CFz(avg)": avg_total_forces["CFz"],
+ "CMx(avg)": avg_total_forces["CMx"],
+ "CMy(avg)": avg_total_forces["CMy"],
+ "CMz(avg)": avg_total_forces["CMz"],
+ }
+ cases_forces.append(forces)
+
+ df_data = pd.DataFrame.from_dict(cases_params)
+ df_forces = pd.DataFrame.from_dict(cases_forces)
+ df = df_data.join(df_forces)
+ df.to_csv(csv_path, index=False, mode="a")
+
+ return csv_path
+
+
+######################################################################################################################
+def create_mesh_params(project):
+ """
+
+ Create the parameters object containing all the information required by the solver to generate the mesh
+
+ Parameters
+ ----------
+ project
+
+ Returns
+ -------
+ params object
+ """
+
+ geometry = project.geometry
+ geometry.group_faces_by_tag("faceName")
+ geometry.group_edges_by_tag("edgeName")
+ with fl.SI_unit_system:
+ params = fl.SimulationParams(
+ meshing=fl.MeshingParams(
+ defaults=fl.MeshingDefaults(
+ boundary_layer_first_layer_thickness=1e-5, surface_max_edge_length=1
+ ),
+ volume_zones=[fl.AutomatedFarfield()],
+ refinements=[
+ fl.SurfaceEdgeRefinement(
+ name="leading_edges",
+ edges=[geometry["leadingEdge"]],
+ method=fl.AngleBasedRefinement(value=2 * u.deg),
+ )
+ ],
+ )
+ )
+ return params
+
+
+def assign_boundary_conditions(project):
+
+ if project.project_tree.root.asset_type == "Geometry":
+ geo = project.geometry
+ geo.group_faces_by_tag("faceName")
+ models = [
+ fl.Wall(
+ surfaces=[
+ # "*" will select all geometry boundaries
+ geo["*"],
+ ],
+ ),
+ fl.Freestream(surfaces=fl.AutomatedFarfield().farfield),
+ ]
+ else:
+ vm = project.volume_mesh
+ models = [
+ fl.Wall(
+ surfaces=[
+ vm["fluid/fuselage"],
+ # *_pylon will select all boundaries ending with _pylon
+ vm["fluid/*_pylon"],
+ vm["fluid/left_wing"],
+ vm["fluid/right_wing"],
+ vm["fluid/h_tail"],
+ vm["fluid/v_tail"],
+ ],
+ ),
+ fl.Freestream(surfaces=vm["fluid/farfield"]),
+ ]
+ return models
+
+
+def project_from_volume_mesh():
+ EVTOL.get_files()
+ project_name = "sweep_evtol_from_mesh"
+ project = fl.Project.from_volume_mesh(
+ EVTOL.mesh_filename, # mesh could also be in ugrid format
+ name=project_name,
+ length_unit="m", # length_unit should be 'm', 'mm', 'cm', 'inch' or 'ft'
+ )
+
+ return project
+
+
+def project_from_geometry():
+ EVTOL.get_files()
+ project_name = "sweep_evtol_from_geometry"
+ project = fl.Project.from_geometry(EVTOL.geometry, name=project_name)
+ mesh_params = create_mesh_params(project)
+ project.generate_volume_mesh(params=mesh_params, name="mesh_name")
+
+ return project
+
+
+######################################################################################################################
+def main():
+ """
+ Main function that drives the functions to achieve all the steps mentioned at the top.
+ """
+
+ # Step1: Connect to an existing project and volume mesh.
+ # Chose one of two options below
+
+ # Option 1a: If you want to upload a new mesh and create a new project.
+ project = project_from_volume_mesh()
+
+ # Option 1b: If you want to upload a CAD geometry and create a new project.
+ # project = project_from_geometry()
+
+ # Option 1c: if you want to run from an existing project.
+ # project = fl.Project.from_cloud(
+ # 'prj-XXXXXXXXXXX') # where prj-XXXXXXXXXX is an ID that can be saved from a previously created project or read off the WEBUI
+
+ vm = project.volume_mesh
+ # If the project has more then one mesh then use hte line below to choose a specific mesh instead.
+ # vm = project.get_volume_mesh(asset_id='vm-XXXXXXXXXXXXXXX')
+
+ print(f"The project id is {project.id}")
+ print(f"The volume mesh contains the following boundaries:{vm.boundary_names}")
+ print(f"The volume mesh ID is: {vm.id}")
+
+ # Step 2: Create the directories to locally store relevant data.
+ dir_name = "evtol_alpha_sweep"
+
+ # Step3: Launch the cases and save the relevant data.
+ models = assign_boundary_conditions(project)
+ params = make_run_params(vm, models)
+ csv_path = launch_sweep(params, project, vm, dir_name)
+
+ generate_report(
+ *csv_reader(csv_path),
+ include_geometry=True,
+ include_general_tables=True,
+ include_residuals=True,
+ include_cfl=True,
+ include_forces_moments_table=True,
+ include_forces_moments_charts=True,
+ include_forces_moments_alpha_charts=True,
+ include_forces_moments_beta_charts=True,
+ include_cf_vec=True,
+ include_cp=True,
+ include_yplus=True,
+ include_qcriterion=True,
+ )
+
+
+######################################################################################################################
+if __name__ == "__main__":
+ main()
diff --git a/examples/post_processing/special_features/bet_disk.py b/examples/post_processing/special_features/bet_disk.py
index 2501b3e30..20903fe4e 100644
--- a/examples/post_processing/special_features/bet_disk.py
+++ b/examples/post_processing/special_features/bet_disk.py
@@ -4,6 +4,14 @@
import flow360 as fl
from flow360.examples import BETExampleData
+from flow360.plugins.report.report import ReportTemplate
+from flow360.plugins.report.report_items import (
+ Chart3D,
+ FrontCamera,
+ Inputs,
+ Settings,
+ Summary,
+)
BETExampleData.get_files()
@@ -86,4 +94,29 @@
# save converted results to a new CSV file:
results.bet_forces.to_file(os.path.join(case.name, "bet_forces_in_SI.csv"))
-# Report workflows have moved to the standalone `flow360-report` package.
+cases = [case]
+
+front_camera_slice = FrontCamera(dimension=350, dimension_dir="height")
+
+bet_slice_screenshot = Chart3D(
+ section_title="BET effective AoA",
+ items_in_row=2,
+ force_new_page=True,
+ show="slices",
+ include=["slice_x"],
+ field="betMetrics_AlphaDegrees",
+ limits=(-18, 0),
+ camera=front_camera_slice,
+ fig_name="slice_x",
+)
+
+report = ReportTemplate(
+ title="BET results screenshots",
+ items=[Summary(), Inputs(), bet_slice_screenshot],
+ settings=Settings(dpi=150),
+)
+
+report = report.create_in_cloud("BET, dpi=default", cases)
+
+report.wait()
+report.download("report.pdf")
diff --git a/examples/post_processing/special_features/heat_transfer.py b/examples/post_processing/special_features/heat_transfer.py
index d67ba3dfd..1c3e26b7b 100644
--- a/examples/post_processing/special_features/heat_transfer.py
+++ b/examples/post_processing/special_features/heat_transfer.py
@@ -1,6 +1,16 @@
import flow360 as fl
from flow360 import u
from flow360.examples import TutorialCHTSolver
+from flow360.plugins.report.report import ReportTemplate
+from flow360.plugins.report.report_items import (
+ Camera,
+ Chart3D,
+ FrontCamera,
+ Inputs,
+ LeftCamera,
+ Settings,
+ Summary,
+)
TutorialCHTSolver.get_files()
@@ -105,4 +115,61 @@
surface_heat_transfer = results.surface_heat_transfer.as_dataframe()
print(surface_heat_transfer)
-# Report generation utilities now live in the external `flow360-report` package.
+cases = [case]
+
+exclude = ["fluid/farfield", "solid/interface_fluid", "solid/adiabatic"]
+
+front_camera_slice = FrontCamera(dimension=1, dimension_dir="width")
+side_camera_slice = LeftCamera(pan_target=(0.35, 0, 0), dimension=2, dimension_dir="width")
+front_right_top_camera = Camera(
+ position=(-1, -1, 1), look_at=(0.35, 0, 0), dimension=1, dimension_dir="width"
+)
+
+x_slice_screenshot = Chart3D(
+ section_title="Slice temperature at x=0.35",
+ items_in_row=2,
+ force_new_page=True,
+ show="slices",
+ include=["slice_x"],
+ field="T",
+ limits=(285 * u.K, 395 * u.K),
+ camera=front_camera_slice,
+ fig_name="slice_x",
+)
+
+y_slice_screenshot = Chart3D(
+ section_title="Slice temperature at y=0",
+ items_in_row=2,
+ force_new_page=True,
+ show="slices",
+ include=["slice_y"],
+ field="T",
+ limits=(285 * u.K, 395 * u.K),
+ camera=side_camera_slice,
+ fig_name="slice_y",
+)
+
+surface_screenshot = Chart3D(
+ section_title="Surface temperature",
+ items_in_row=2,
+ force_new_page=True,
+ show="boundaries",
+ field="T",
+ limits=(285 * u.K, 395 * u.K),
+ exclude=exclude,
+ camera=front_right_top_camera,
+)
+
+report = ReportTemplate(
+ title="CHT results screenshots",
+ items=[Summary(), Inputs(), x_slice_screenshot, y_slice_screenshot, surface_screenshot],
+ settings=Settings(dpi=150),
+)
+
+report = report.create_in_cloud(
+ "CHT, dpi=default",
+ cases,
+)
+
+report.wait()
+report.download("report.pdf")
diff --git a/examples/tutorials/notebooks/notebook_automotive.ipynb b/examples/tutorials/notebooks/notebook_automotive.ipynb
index 175cada4a..2a1790b55 100644
--- a/examples/tutorials/notebooks/notebook_automotive.ipynb
+++ b/examples/tutorials/notebooks/notebook_automotive.ipynb
@@ -785,6 +785,41 @@
">Note: Importing units (u) here is only for the sake of convenience."
]
},
+ {
+ "cell_type": "code",
+ "execution_count": 10,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "from flow360 import u\n",
+ "from flow360.plugins.report.report import ReportTemplate\n",
+ "from flow360.plugins.report.report_items import (\n",
+ " BottomCamera,\n",
+ " Chart2D,\n",
+ " Chart3D,\n",
+ " FrontCamera,\n",
+ " FrontLeftBottomCamera,\n",
+ " FrontLeftTopCamera,\n",
+ " Inputs,\n",
+ " LeftCamera,\n",
+ " RearCamera,\n",
+ " RearLeftTopCamera,\n",
+ " RearRightBottomCamera,\n",
+ " Settings,\n",
+ " Summary,\n",
+ " Table,\n",
+ " TopCamera,\n",
+ ")\n",
+ "from flow360.plugins.report.utils import (\n",
+ " Average,\n",
+ " DataItem,\n",
+ " Delta,\n",
+ " Expression,\n",
+ " Variable,\n",
+ " GetAttribute,\n",
+ ")"
+ ]
+ },
{
"cell_type": "markdown",
"metadata": {},
@@ -994,12 +1029,308 @@
")"
]
},
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### Visualizations\n",
+ "\n",
+ "Here we will define separate lists of `Chart3D`, which are essentially 3D visualizations (you can think of them as screenshots of your model). These come with a variety of options that can be used to customize the way data is shown.\n",
+ "\n",
+ "One of the common things is that pretty much all of those charts contain either an exclude or include field. They provide control over what parts of the simulation domain will be excluded/included in the image.\n",
+ "\n",
+ ">Note: If you took a peek at the charts you may have noticed that list comprehension can used to iterate over camera positions."
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 16,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "forces = [\n",
+ " Chart2D(\n",
+ " x=\"x_slicing_force_distribution/X\",\n",
+ " y=\"x_slicing_force_distribution/totalCumulative_CD_Curve\",\n",
+ " fig_name=\"totalCumulative_CD_Curve\",\n",
+ " background=\"geometry\",\n",
+ " exclude=exclude,\n",
+ " ),\n",
+ " Chart2D(\n",
+ " x=\"surface_forces/pseudo_step\",\n",
+ " y=\"surface_forces/totalCD\",\n",
+ " section_title=\"Drag Coefficient\",\n",
+ " fig_name=\"cd_fig\",\n",
+ " exclude=exclude,\n",
+ " focus_x=(1 / 3, 1),\n",
+ " ),\n",
+ "]\n",
+ "geometry_screenshots = [\n",
+ " Chart3D(\n",
+ " section_title=\"Geometry\",\n",
+ " items_in_row=2,\n",
+ " force_new_page=True,\n",
+ " show=\"boundaries\",\n",
+ " camera=camera,\n",
+ " exclude=exclude,\n",
+ " fig_name=f\"geo_{i}\",\n",
+ " )\n",
+ " for i, camera in enumerate(cameras_geo)\n",
+ "]\n",
+ "cpt_screenshots = [\n",
+ " Chart3D(\n",
+ " section_title=\"Isosurface, Cpt=-1\",\n",
+ " items_in_row=2,\n",
+ " force_new_page=True,\n",
+ " show=\"isosurface\",\n",
+ " iso_field=\"Cpt\",\n",
+ " exclude=exclude,\n",
+ " camera=camera,\n",
+ " )\n",
+ " for camera in cameras_cp\n",
+ "]\n",
+ "cfvec_screenshots = [\n",
+ " Chart3D(\n",
+ " section_title=\"CfVec\",\n",
+ " items_in_row=2,\n",
+ " force_new_page=True,\n",
+ " show=\"boundaries\",\n",
+ " field=\"CfVec\",\n",
+ " mode=\"lic\",\n",
+ " limits=(1e-4, 10),\n",
+ " is_log_scale=True,\n",
+ " exclude=exclude,\n",
+ " camera=camera,\n",
+ " )\n",
+ " for camera in cameras_cp\n",
+ "]\n",
+ "y_slices_screenshots = [\n",
+ " Chart3D(\n",
+ " section_title=f\"Slice velocity y={y}\",\n",
+ " items_in_row=2,\n",
+ " force_new_page=True,\n",
+ " show=\"slices\",\n",
+ " include=[f\"slice_y_{name}\"],\n",
+ " field=\"velocity\",\n",
+ " limits=(0 * u.m / u.s, 50 * u.m / u.s),\n",
+ " camera=side_camera_slice,\n",
+ " fig_name=f\"slice_y_{name}\",\n",
+ " )\n",
+ " for name, y in zip([\"0\", \"0_2\", \"0_4\", \"0_6\", \"0_8\"], [0, 0.2, 0.4, 0.6, 0.8])\n",
+ "]\n",
+ "y_slices_lic_screenshots = [\n",
+ " Chart3D(\n",
+ " section_title=f\"Slice velocity LIC y={y}\",\n",
+ " items_in_row=2,\n",
+ " force_new_page=True,\n",
+ " show=\"slices\",\n",
+ " include=[f\"slice_y_{name}\"],\n",
+ " field=\"velocityVec\",\n",
+ " mode=\"lic\",\n",
+ " limits=(0 * u.m / u.s, 50 * u.m / u.s),\n",
+ " camera=side_camera_slice,\n",
+ " fig_name=f\"slice_y_vec_{name}\",\n",
+ " )\n",
+ " for name, y in zip([\"0\", \"0_2\", \"0_4\", \"0_6\", \"0_8\"], [0, 0.2, 0.4, 0.6, 0.8])\n",
+ "]\n",
+ "z_slices_screenshots = [\n",
+ " Chart3D(\n",
+ " section_title=f\"Slice velocity z={z}\",\n",
+ " items_in_row=2,\n",
+ " force_new_page=True,\n",
+ " show=\"slices\",\n",
+ " include=[f\"slice_z_{name}\"],\n",
+ " field=\"velocity\",\n",
+ " limits=(0 * u.m / u.s, 50 * u.m / u.s),\n",
+ " camera=top_camera_slice,\n",
+ " fig_name=f\"slice_z_{name}\",\n",
+ " )\n",
+ " for name, z in zip([\"neg0_2\", \"0\", \"0_2\", \"0_4\", \"0_6\", \"0_8\"], [-0.2, 0, 0.2, 0.4, 0.6, 0.8])\n",
+ "]\n",
+ "y_plus_screenshots = [\n",
+ " Chart3D(\n",
+ " section_title=\"y+\",\n",
+ " items_in_row=2,\n",
+ " show=\"boundaries\",\n",
+ " field=\"yPlus\",\n",
+ " exclude=exclude,\n",
+ " limits=(0, 5),\n",
+ " camera=camera,\n",
+ " fig_name=f\"yplus_{i}\",\n",
+ " )\n",
+ " for i, camera in enumerate([top_camera, bottom_camera])\n",
+ "]\n",
+ "cp_screenshots = [\n",
+ " Chart3D(\n",
+ " section_title=\"Cp\",\n",
+ " items_in_row=2,\n",
+ " show=\"boundaries\",\n",
+ " field=\"Cp\",\n",
+ " exclude=exclude,\n",
+ " limits=limits,\n",
+ " camera=camera,\n",
+ " fig_name=f\"cp_{i}\",\n",
+ " )\n",
+ " for i, (limits, camera) in enumerate(zip(limits_cp, cameras_cp))\n",
+ "]\n",
+ "cpx_screenshots = [\n",
+ " Chart3D(\n",
+ " section_title=\"Cpx\",\n",
+ " items_in_row=2,\n",
+ " show=\"boundaries\",\n",
+ " field=\"Cpx\",\n",
+ " exclude=exclude,\n",
+ " limits=(-0.3, 0.3),\n",
+ " camera=camera,\n",
+ " fig_name=f\"cpx_{i}\",\n",
+ " )\n",
+ " for i, camera in enumerate(cameras_cp)\n",
+ "]\n",
+ "wall_shear_screenshots = [\n",
+ " Chart3D(\n",
+ " section_title=\"Wall shear stress magnitude\",\n",
+ " items_in_row=2,\n",
+ " show=\"boundaries\",\n",
+ " field=\"wallShearMag\",\n",
+ " exclude=exclude,\n",
+ " limits=(0 * u.Pa, 5 * u.Pa),\n",
+ " camera=camera,\n",
+ " fig_name=f\"wallShearMag_{i}\",\n",
+ " )\n",
+ " for i, camera in enumerate(cameras_cp)\n",
+ "]"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### Report assembly\n",
+ "\n",
+ "Once we have all of the components of our report, we can create a `ReportTemplate` that will allow us to then generate the report itself. In our case, we will also add two `Chart2D`s related to drag coefficient."
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 17,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "report_template = ReportTemplate(\n",
+ " title=\"Aerodynamic analysis of DrivAer\",\n",
+ " items=[\n",
+ " Summary(),\n",
+ " Inputs(),\n",
+ " statistical_table_general,\n",
+ " statistical_table_forces,\n",
+ " *forces,\n",
+ " *geometry_screenshots,\n",
+ " *cp_screenshots,\n",
+ " *cpx_screenshots,\n",
+ " *cpt_screenshots,\n",
+ " *y_slices_screenshots,\n",
+ " *y_slices_lic_screenshots,\n",
+ " *z_slices_screenshots,\n",
+ " *y_plus_screenshots,\n",
+ " *wall_shear_screenshots,\n",
+ " ],\n",
+ " settings=Settings(dpi=150),\n",
+ ")"
+ ]
+ },
{
"cell_type": "markdown",
"metadata": {},
"source": [
"All that is left now is to create the report by calling the `create_in_cloud()` function, then downloading it. To ensure that the download will only start after the report has actually been generated, we use `wait()`."
]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 18,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "application/vnd.jupyter.widget-view+json": {
+ "model_id": "87aaf851429a4b2ba59e36507a9dd1b8",
+ "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": "69c8a4f4849f48898412f0dc4755f344",
+ "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": [
+ "[14:15:44] INFO: Saved to report.pdf \n",
+ "
\n"
+ ],
+ "text/plain": [
+ "\u001b[2;36m[14:15:44]\u001b[0m\u001b[2;36m \u001b[0m\u001b[36mINFO\u001b[0m: Saved to report.pdf \n"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ },
+ {
+ "data": {
+ "text/plain": [
+ "'report.pdf'"
+ ]
+ },
+ "execution_count": 18,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "report = report_template.create_in_cloud(\n",
+ " f\"{size}-{len(cases)}cases-slices-using-groups-Cpt, Cpx, wallShear, dpi=default\",\n",
+ " cases,\n",
+ ")\n",
+ "\n",
+ "report.wait()\n",
+ "report.download(\"report.pdf\")"
+ ]
}
],
"metadata": {
@@ -1023,4 +1354,4 @@
},
"nbformat": 4,
"nbformat_minor": 2
-}
\ No newline at end of file
+}
diff --git a/flow360/__init__.py b/flow360/__init__.py
index 30fe688d7..c8ad893b9 100644
--- a/flow360/__init__.py
+++ b/flow360/__init__.py
@@ -203,6 +203,7 @@
from flow360.component.surface_mesh_v2 import SurfaceMeshV2 as SurfaceMesh
from flow360.component.volume_mesh import VolumeMeshV2 as VolumeMesh
from flow360.environment import Env
+from flow360.plugins import report
from flow360.version_utils import warn_if_prerelease_version as _warn_prerelease
__all__ = [
@@ -346,6 +347,7 @@
"UserVariable",
"math",
"solution",
+ "report",
"snappy",
"ModularMeshingWorkflow",
"SeedpointVolume",
diff --git a/flow360/cloud/flow360_requests.py b/flow360/cloud/flow360_requests.py
index 4e294ebce..79410ed38 100644
--- a/flow360/cloud/flow360_requests.py
+++ b/flow360/cloud/flow360_requests.py
@@ -183,6 +183,20 @@ class NewSurfaceMeshDependencyRequest(Flow360RequestsV2):
icon: str = pd_v2.Field(default="", description="project description")
+class _Resource(Flow360RequestsV2):
+ type: Literal["Case", "Project"]
+ id: str
+
+
+class NewReportRequest(Flow360RequestsV2):
+ "New report request"
+
+ name: str
+ resources: List[_Resource]
+ config_json: str
+ solver_version: str
+
+
class DraftCreateRequest(Flow360RequestsV2):
"""Data model for draft create request"""
diff --git a/flow360/component/project.py b/flow360/component/project.py
index cab7a6c55..ce4f632e4 100644
--- a/flow360/component/project.py
+++ b/flow360/component/project.py
@@ -80,6 +80,7 @@
Flow360WebError,
)
from flow360.log import log
+from flow360.plugins.report.report import get_default_report_summary_template
from flow360.version import __solver_version__
AssetOrResource = Union[type[AssetBase], type[Flow360Resource]]
@@ -2214,4 +2215,10 @@ def run_case(
draft = case_or_draft
return draft
case = case_or_draft
+ report_template = get_default_report_summary_template()
+ report_template.create_in_cloud(
+ name=f"{name}-summary",
+ cases=[case],
+ solver_version=solver_version if solver_version else self.solver_version,
+ )
return case
diff --git a/flow360/component/simulation/services.py b/flow360/component/simulation/services.py
index c341a7540..b4e7d8c44 100644
--- a/flow360/component/simulation/services.py
+++ b/flow360/component/simulation/services.py
@@ -94,6 +94,7 @@
Flow360TranslationError,
Flow360ValueError,
)
+from flow360.plugins.report.report import get_default_report_summary_template
from flow360.version import __version__
# Required for correct global scope initialization
@@ -1181,6 +1182,19 @@ def translate_xfoil_c81_bet_disk(
return bet_dict_list, errors
+def get_default_report_config() -> dict:
+ """
+ Get the default report config
+ Returns
+ -------
+ dict
+ default report config
+ """
+ return get_default_report_summary_template().model_dump(
+ exclude_none=True,
+ )
+
+
def _parse_root_item_type_from_simulation_json(*, param_as_dict: dict):
"""[External] Deduct the root item entity type from simulation.json"""
try:
diff --git a/flow360/plugins/__init__.py b/flow360/plugins/__init__.py
deleted file mode 100644
index 205e11b5b..000000000
--- a/flow360/plugins/__init__.py
+++ /dev/null
@@ -1,22 +0,0 @@
-"""
-Flow360 plugin namespace.
-
-Report-related functionality has moved to the standalone `flow360-report` package.
-Install it via `pip install flow360-report` if you need reporting utilities.
-"""
-
-from typing import Any
-
-_REPORT_HINT = (
- "flow360 reporting utilities now live in the `flow360-report` package. "
- "Install it via `pip install flow360-report` and import from `flow360_report.plugins.report`."
-)
-
-
-def __getattr__(name: str) -> Any:
- if name == "report":
- raise ImportError(_REPORT_HINT)
- raise AttributeError(f"module 'flow360.plugins' has no attribute '{name}'")
-
-
-__all__: list[str] = []
diff --git a/flow360/plugins/report/README.md b/flow360/plugins/report/README.md
new file mode 100644
index 000000000..23ea80864
--- /dev/null
+++ b/flow360/plugins/report/README.md
@@ -0,0 +1,311 @@
+# Report Plugin
+
+The **Report Plugin** provides automated PDF report generation for simulation results. It integrates with Flow360 (or similar frameworks) to pull data from simulation cases, produce tables and charts, and optionally render 3D images or screenshots. Users can easily configure how the PDF is structured, which items to include (e.g., summaries, inputs, tables, 2D and 3D charts), and how to style each item.
+
+---
+
+## Table of Contents
+- [Report Plugin](#report-plugin)
+ - [Table of Contents](#table-of-contents)
+ - [Overview](#overview)
+ - [Features](#features)
+ - [Folder and File Structure](#folder-and-file-structure)
+ - [Installation](#installation)
+ - [Usage](#usage)
+ - [1. Create or Load Simulation Cases](#1-create-or-load-simulation-cases)
+ - [2. Prepare a Report Configuration](#2-prepare-a-report-configuration)
+ - [3. Invoke the Plugin](#3-invoke-the-plugin)
+ - [Configuration and Customization](#configuration-and-customization)
+ - [Generating Screenshots](#generating-screenshots)
+ - [Camera controls:](#camera-controls)
+ - [Key Camera Attributes](#key-camera-attributes)
+ - [Camera Examples](#camera-examples)
+ - [Examples](#examples)
+ - [Minimal Example](#minimal-example)
+
+---
+
+## Overview
+
+This plugin:
+- Loads simulation data from one or more **Case** objects.
+- Extracts relevant metrics, such as velocities, forces, boundaries, and other fields.
+- Creates custom sections (e.g. *Summary*, *Inputs*, *Tables*, *Charts*) and assembles them into a single LaTeX-based PDF document.
+- Can optionally include 3D images or screenshots of the geometry, surfaces, or isosurfaces for an at-a-glance view of simulation geometry or results.
+
+---
+
+## Features
+
+1. **Section-Based Reporting**
+ - Easily add sections like `Summary`, `Inputs`, and custom sections with specialized content.
+
+2. **Tables**
+ - Collect data from one or multiple simulation cases and display in tabular form using the `Table` class.
+ - Built-in support for custom formatting (e.g., `.5g` or specialized formatting functions).
+
+3. **2D Charts**
+ - Plot data such as forces, residuals, or general x-y plots with `Chart2D`.
+ - Supports logarithmic scaling if needed (e.g., for residual plots).
+ - Can combine multiple cases into one plot or generate separate plots per case.
+
+4. **3D Visualization**
+ - Automatically generate (or embed) images of your 3D geometry (e.g., surfaces, isosurfaces, different camera angles).
+ - Users can opt to generate screenshots of simulation geometry in various positions (requires an external service or local tool to produce the images; see [Generating Screenshots](#generating-screenshots) below).
+
+5. **Flexible Layout**
+ - Choose how items are arranged (one item per page, multiple items per row, etc.).
+ - Control figure sizes, captions, new pages, etc.
+
+---
+
+## Folder and File Structure
+
+A high-level view of the relevant files:
+
+flow360/plugins/report/
+├── report.py # Main file that defines the ReportTemplate class
+├── report_items.py # Contains classes for different report "items" like Summary, Table, Chart2D, Chart3D, etc.
+├── report_context.py # (Imported by others) Manages the context in which the report is generated
+├── utils.py # Utility functions for data extraction, data formatting, etc.
+└── uvf_shutter.py # (Optional) Helper for generating 3D geometry screenshots
+
+Example files in the repo:
+
+examples/
+└── automotive_example_in_cloud.py # Shows how to configure and instantiate a report
+
+## Installation
+
+1. Make sure you have Python 3.8+ installed.
+2. Install flow360 package:
+ ```bash
+ pip install flow360
+ ```
+
+
+
+
+---
+
+## Usage
+
+### 1. Create or Load Simulation Cases
+
+You need one or more **Case** objects that contain all relevant simulation data. The plugin references fields in these cases to build tables and charts.
+
+### 2. Prepare a Report Configuration
+
+Typically, you have a JSON or Python-based config specifying which `ReportItem` objects you want in your report, e.g.:
+
+```python
+report = ReportTemplate(
+ title="Aerodynamic analysis of DrivAer",
+ items=[
+ Summary(text="Short summary..."),
+ Inputs(),
+ Table(
+ data=[
+ "params/reference_geometry/area",
+ DataItem(data="surface_forces/totalCD", operations=Average(fraction=0.1)),
+ "params/time_stepping/max_steps",
+ ],
+ section_title="Statistical data",
+ ),
+ Chart2D(
+ x="surface_forces/pseudo_step",
+ y="surface_forces/totalCD",
+ section_title="Drag Coefficient",
+ fig_name="cd_fig",
+ )
+ ]
+)
+report.to_file('path/to/reportTemplate.json"')
+```
+(This might be stored in a file like reportTemplate.json.)
+
+```JSON
+{
+ "include_case_by_case": false,
+ "items": [
+ {
+ "text": "Short summary...",
+ "type_name": "Summary"
+ },
+ {
+ "type_name": "Inputs"
+ },
+ {
+ "data": [
+ "params/reference_geometry/area",
+ {
+ "data": "surface_forces/totalCD",
+ "operations": [
+ {
+ "fraction": 0.1,
+ "type_name": "Average"
+ }
+ ],
+ "type_name": "DataItem"
+ },
+ "params/time_stepping/max_steps"
+ ],
+ "section_title": "Statistical data",
+ "type_name": "Table"
+ },
+ {
+ "fig_name": "cd_fig",
+ "section_title": "Drag Coefficient",
+ "type_name": "Chart2D",
+ "x": "surface_forces/pseudo_step",
+ "y": "surface_forces/totalCD"
+ }
+ ],
+ "title": "Aerodynamic analysis of DrivAer"
+}
+```
+
+
+
+### 3. Invoke the Plugin
+Use the main ReportTemplate class (from report.py) to build and export your PDF. For example (similar to automotive_example_in_cloud.py):
+
+```python
+from flow360.plugins.report.report import ReportTemplate
+from flow360.component.case import Case
+
+# Suppose you have a list of Case objects
+cases = [...]
+
+# Point to your config:
+report = ReportTemplate(filename="path/to/reportTemplate.json")
+
+# Generate the PDF:
+report.create_pdf(
+ output_folder="my_report_output",
+ cases=cases
+)
+```
+
+
+## Configuration and Customization
+
+- **Report Items**: The core building blocks, each corresponding to a section or figure in the PDF.
+ - *Summary*: Adds a textual overview.
+ - *Inputs*: Collects key input parameters from your cases (e.g., velocity, time-stepping).
+ - *Table*: Aggregates data from multiple cases into one table.
+ - *Chart2D*: Produces an XY plot from data paths in your case(s).
+ - *Chart3D*: Optionally display geometry or color-field surfaces (like `Cp`, `yPlus`, etc.).
+
+- **Custom Data Paths**: Many items accept `data` arrays where each entry references a JSON path inside each case, e.g. `"params/operating_condition/velocity_magnitude"`.
+ - If the plugin finds that path in each case’s data, it extracts the values and populates the table/plot.
+
+- **Formatting**:
+ - Use built-in format specifiers or pass a custom function to format numeric data.
+ - For tables, set `formatter` as a string (e.g. `".5g"`) or a callable for specialized formatting.
+
+- **Multiple Cases**: If you provide multiple cases in a single run, the plugin can:
+ - Combine them side by side in tables.
+ - Overlay them on the same chart.
+ - Generate separate pages or subplots (use `separate_plots=True` or `items_in_row=...` to configure layout).
+
+
+
+## Generating Screenshots
+
+The plugin can produce screenshots of your simulation geometry or results. You can activate this in `Chart3D` items by specifying a field or object type to display. The actual screenshot generation leverages an external tool, but the plugin manages how the images are embedded in the final PDF.
+
+
+### Camera controls:
+
+The `Camera` class specifies a 3D viewpoint and zoom for rendering geometry or visualizing simulation results. By adjusting properties such as `position`, `look_at`, and `dimension`, users can control the angle, orientation, and scale of the view. Below are the attributes you can configure (note that all lengths are in the same units used in your geometry or volume mesh):
+
+
+### Key Camera Attributes
+
+- **Position**
+ The camera’s eye position in 3D space. Typically, you can think of it as a point on a sphere, looking inward toward `look_at`. Adjusting `position` can rotate the view around the object of interest.
+
+- **Up**
+ A vector that determines which way is “up” in the final image. Changing this can rotate the view around the axis from the camera’s position to the `look_at` point.
+
+- **Look At**
+ The point at which the camera is aimed. By default, if unspecified, it’s often the bounding-box center of your model.
+
+- **Pan Target**
+ A specific point to which you want to pan your camera’s center. If left unset, `look_at` remains the default center of view.
+
+- **Dimension & Dimension Direction**
+ Together, these control the zoom level. For example, if `dimension_dir` is `"width"` and `dimension` is `2.0`, the rendered view’s width represents 2 model units. Similarly, use `"height"` or `"diagonal"` to scale in those directions.
+
+
+---
+
+### Camera Examples
+
+
+
+See how changing `position` and `look_at` affects the orientation of the camera.
+
+
+
+See how `pan_target` changes the view.
+
+
+
+See the effect of `dimension` and `dimension_dir` on zoom level.
+
+
+
+
+## Examples
+
+### Minimal Example
+
+```python
+from flow360.plugins.report.report import ReportTemplate
+from flow360.plugins.report.report_items import *
+
+cases = [...] # Load or define your cases
+
+report = ReportTemplate(
+ title="Aerodynamic analysis of DrivAer",
+ items=[
+ Summary(text="Short summary..."),
+ Inputs(),
+ Table(
+ data=[
+ "params/reference_geometry/area",
+ DataItem(data="surface_forces/totalCD", operations=Average(fraction=0.1)),
+ "params/time_stepping/max_steps",
+ ],
+ section_title="Statistical data",
+ ),
+ Chart2D(
+ x="surface_forces/pseudo_step",
+ y="surface_forces/totalCD",
+ section_title="Drag Coefficient",
+ fig_name="cd_fig",
+ )
+ ]
+)
+
+# Save the JSON config or pass it directly:
+report.create_pdf(
+ "my_report",
+ cases
+)
+report.wait()
+report.download("my_report.pdf")
+```
+
+
+After running, check the my_report.pdf for your compiled PDF.
+
+More Advanced
+See `automotive_example_in_cloud.py` for more advanced usage (multiple 3D fields, custom cameras, integration with a pipeline, etc.).
+
+
+---
+
+**Enjoy creating automated simulation reports with the Report Plugin!**
diff --git a/flow360/plugins/report/__init__.py b/flow360/plugins/report/__init__.py
new file mode 100644
index 000000000..e3dd44434
--- /dev/null
+++ b/flow360/plugins/report/__init__.py
@@ -0,0 +1,69 @@
+"""Utilities for creating reports"""
+
+from flow360.plugins.report.report import ReportTemplate
+from flow360.plugins.report.report_items import (
+ Chart2D,
+ Chart3D,
+ FixedRangeLimit,
+ Inputs,
+ ManualLimit,
+ NonlinearResiduals,
+ PatternCaption,
+ Settings,
+ SubsetLimit,
+ Summary,
+ Table,
+)
+from flow360.plugins.report.utils import (
+ Average,
+ DataItem,
+ Delta,
+ Expression,
+ GetAttribute,
+ Grouper,
+ Variable,
+)
+from flow360.plugins.report.uvf_shutter import (
+ BottomCamera,
+ Camera,
+ FrontCamera,
+ FrontLeftBottomCamera,
+ FrontLeftTopCamera,
+ LeftCamera,
+ RearCamera,
+ RearLeftTopCamera,
+ RearRightBottomCamera,
+ TopCamera,
+)
+
+__all__ = [
+ "Average",
+ "BottomCamera",
+ "Camera",
+ "Chart2D",
+ "Chart3D",
+ "DataItem",
+ "Delta",
+ "Expression",
+ "FixedRangeLimit",
+ "FrontCamera",
+ "FrontLeftBottomCamera",
+ "FrontLeftTopCamera",
+ "GetAttribute",
+ "Grouper",
+ "Inputs",
+ "LeftCamera",
+ "ManualLimit",
+ "NonlinearResiduals",
+ "PatternCaption",
+ "RearCamera",
+ "RearLeftTopCamera",
+ "RearRightBottomCamera",
+ "ReportTemplate",
+ "Settings",
+ "SubsetLimit",
+ "Summary",
+ "Table",
+ "TopCamera",
+ "Variable",
+]
diff --git a/flow360/plugins/report/fonts/TWKEverett-Bold.otf b/flow360/plugins/report/fonts/TWKEverett-Bold.otf
new file mode 100644
index 0000000000000000000000000000000000000000..d66f27fa1b3fc049c907368d94e529b180aaf6f7
GIT binary patch
literal 178056
zcmc$`2V4|a*D!u(W(&-~0y``sxCIp}AiY`XDhLSJ>jDd`EJawF6?;#tG4>XF?7hWq
z>?RtG8hbRc#}@D6fP80WmxU*9p1jZV`~Sb+hn~IX_H)iX_uM=8&Y7gNv?OGYDiMW}
zqhez2x+z>^rSL6wObrQh(P4NN(gP@qGMtsydoGxNPHo=Yg}?t>ZJYCjEIny
zBVy9*xYS;eV|KpSg2-7m2(>9sN^R$_T{2-dLb8cae`rRYw($GbDO!Zwnjj=8%F$^v
zJxU(khrGuC!8wp%e^ho1$`6Hj^PD_md1t?{fzTI}w~x-#mKRbzlOYS;0g&Wt^K@^Y
zv>E{Ap*|B*SWs-NO+KlF(-2qhz`Ctw5sg;3*pn?DS3=n>0N=g{}^iTQC)dJtwd$0C3z0|xC(70p;qK0z#GB^
zq$2+YD0?X13AG@2@dav2tVJ%wDhP)FBm*3Q_&2C6`VQcCBqP2=
zUPMRKmh1v~=K&5t{ul_qgK!RnW1y{cD7O@~Bo-klnT8yQ)i50!H4J(%&BzLi0Fr!qb)qt
zej<7zPaY?#0U({|Bb0B>*FkR~-X6)(3&cRW2N92G;x4oiheSjovLOPIo6vS=n4|xB
ze?s5?)L&hj|JX0zAJ9rU!Y>T{!c3`U-LN!toO61fK$k%8gKzo2M0g?ft0RjM$0J;Da
z0(1lj2WU@}BNvc`Eun+4)4{q|3{U`g1%MCe`zRQ*3XsSDqjKm%Z3wJiev}7NQCR?T
zfJvx1<%-l)0pLLZ8^K1-LQ3iz)P!mV>16;c(okMPStT_N;t{9?)dHzSYS6Wjz*|FT
zi6V%WCDeZpzvLG^__P-|FQ-H|g7L;#9_avDMd
za&3&dBD&Tatw&AKF3^S7M00?~eA^)J(WtWalO?Yu?4{HI&akevAclagoPt~dnu7e}
z0o)-^1Hcm5s1<;Sn2S^ZEs0Xt!;~Tw`iPVODs&%t12hJ_NrSj6tbgG=knK=7)z5;*
zK>d@b0~v@qP=iqi^cr~pU03;#mjNjSIlMuwi5Do8a76wDjeN<2s5L-KSc6*wv?Pw9
zcF2U7yw?U69F;+QUTNe0bq-o3J?k*
z0BM9d*%kSb-C!&xL7j1^vjE93w#ks526H71g^)9VzI`Z=Tnps~01aJ`hiEwPzrMDb
z+6?`Ei`u|GOib=Z-N;5aw^Ki6T)p
zKpa3zOUTD30fbvZI|%u-NK0q~VG@840>8cj909(!18)G-VOW11z`orDzDB^@T7vwo
zPzgM51iK@JasLTy-$Kxv)}XI`U|-TfhT*^$gyKf15493(;0e@>2myWVgi`r-LAU%7
z`kaAKZ4U^yn4#u5#9IUIWr0x;egv2UkZ&%_m#^&x&;h^^p_*w3RiB6a(}4FQcEB3+
zJ_Vo!ZMC5`)b5k;
zG-jPAklPc9$RzSG?0H4xbat>SnW~GOAQiM`Zc~c-
zidtH58`MNJ)`HtYnYy|R@)sSjq}!oD(F+T1PtF&sEw~H~mhQ6PEb49JWX6FZHyaB*
z1kh|V&XP_dnax5AP9c@eS_>|M@+U307zNtgwcrxuWn0%K$l12zMvcF(U!b?YufMO4
zCaR#YqDY^eW7MQqXwv#6$N2P#Nl1&2NlDz`qlqvWG(u*vrbt(;D=O7x`tTYPQ>rV{
z8I6$zhRgsTUtiz$d?Np2$?)cL8Why<@!qw4yK8Xf?(nZKf_yTa>FQ$odPx|4M#hfkB(EN!Au>i*yEqKF4S*EDmqi
z&R9{X^C{CAjJdj^e4mVhyt<5hzPxu~U3s63oCIxIR*6AVQmo_qGUn(sSq1sVVogDj
z5XsaUwVDh_(dzScnVK@aF-MY4gtW05|+#K@P{=4V14
zaG0MBBbHyGDakkdg>VoYuw9T{q%F+RXK0EGvW#UQ6(AV~I9;bHEGj70Wt13Ucr~%W
zoi;O5Z-lYZ8Z`R+tb(FEEg#bq8h~(~5S^w@X#Yc+?HW{Q2Q29Fi-EBJvUnPZdU%OG
zBUh7NR8W==VlUSWE6FPa(Szm*A~tA8RA@2_ve8DQfrY^r`Jq7A%rT0MzqiXoR$46b>87c7XHT>VU^G
z_%%R0SCCIWq+~!#dHrb{VxKX*DJnvo+iFkzD-`wjMldKk>vU2>zT)pr;MMS5s<_8
zQUW30!#_#fta0YiE`a*H*6~`d2l?=A^EzjQ7I+=xspO?>p26utyb#h01R62~>ddoO
zgJK0a@SJ7}Jn(e$($GR}jlcsxGDRRgeoX7;Tp{F|C0-{HKK$=>v)#Yxp&hi(bEX4I
z_?ga2?Eh}-Y0y5;!*Gxe-!ngt(*>E8LGNa5FNbg#qdjud2>4jAj!HMv)=IX)Iz-!sKsjouP<>Z)xvq5-kEh)XVJ_w=$Riiel|56
zvkIukk5--i)&U)6-Q;OC&%ys`wfX>cd`mpF8GseRFW)auNxr}{Uz=~W7%VYQtyxC@
zTUpq5@;4!Ha{|u{uxk`Tkz_E0XHg}Rz!Secn6WT4hVUc05@U!jiMhlwVlVM4@jLN?
zs3qOWc4Pn?k?ajA})NQuC>e)M4tFtg~v8YPo8a>VWEq>XPb-
z>btc1`=i8rGe%=8}Q}-vU{ajmJ`{r|-&n}-upNS7}#rE(@U8XsFH2#tH5j^4Xe~(-q
zef8+tqi&DLhrc{r`*1CkhI_P!?;gTA$iv$Ye|!kfT@PPBe2CD4s}HU{IR0P}U+ZD~
z!>9)*9_)bB4G&g6nE#;YK^oxFyXALY-MM_XKSFm~-R?{Gl#Hgn6m;zj-xV>8oC#+S
z!XM>Axl(Ra6G}xjr`o_Nm*tOYLHU9{!SnH7|0o|&-(biKr#euPR5TT9E&+HfoN)f-
z56)*070oBY4V3U_E^kf~(xdBZ)s+#_y9mXFA1bsS2aF;+k`u`oavph_tRW|nePK34
zk=w{HGM@}3W65x`J-L_MMeZRxlhyDHDn~3lDJoHG)CTl10G>}HU@O=Sbw@qXP%vv5
zU=nyEREfr+$!H3kmJB50$U)>0vKY>?e}ePqU(qx43VlKxK@cLsk#Hs&5lTWuc!FjJ
z6Cp%aqFY+Y$nT((UNf+V>Iga!sJ;)Yh7epg5lETxk4a_A6I9>IG6FLt#4Gf2^
zZ7^)(++fq%7tCZD8jO08JJ3KdwIl#1?FoB07nGt7gbdCF#c&tGqDVqPxS(jl6~z*6D28Z^;s|%tm1qK|L25Yt@PN}&
z4N8Pl)c$aioJzDr{RnR~1Wrpc;bbn02uC_14CTOWYcA0l4TBRRJ<$;r5iw{u5seCo
zC^VKxL8FKyR7!LM+u9S2BofhRA{kA98-Z!WKr|hFNu;AWL?)U|WT3f(4mqKRu$`$v
z1BiBTUhx18f^#VY(FGL~v8XrE2F)afz#M5q`GK|ar@Ub%HG-L>qSP>(6qGZigqh`z
zCc&xAc%m1YNc2IraEj3!eSlLU6Py4us1MN=RS*ei2As;wB8HM<$g$)Iauhk590``t
zj*KMFl3U4d$-(3hu=QKW&17$~54o89n*55gr%2Kp>}LSkjBHN6B7Y~JkuS*C{lWuX5V0MCj6OIavU030i9p8siSS+8X9udzf)?
zFx&c|K``6&uwfVhb8Ra65-otUrgdl&+JO$CGw33^2AcX9H0vX(CZvQt%rrLwPLm05
z(9IB{6U?)&L=w@P7(@&sh7+TR3B)v5Sr)(=xq;XQ>*NvQ6mg!oOk5{^ftBVD;!ole
zNsuc2I|?)6^B}7WFIjimDb#M2;diQFD=xC{z?F
zN)YuH4HXR&8AW46Q$=${OGWEM+eHUOr$m=UH${&`FGU|k)nZ007dIAb#4W`>;y`gH
zakMx|+)F$_tQBXA^Tp-jvEs?%S>pNP#p0FX&ElQn!{YD77sNk`Z;Ky_pNjtw|0%AP
zNF{PfV~MB4OX4pHm2{9qO1etoB`K1Ak|B~TNuFecWSnHGWS(S^WQAm%LhiOdP>_${iGeFk0s0473?yGt8#Irov`|
z%`BURHmhtl+w8G9YIDx!n$2CC-)-L7aJFJw2iwNB&1}7GLu?~#<86D{4zbnS7Tb=t
zooTzocCGC;+XJ>IZ7d2%o=7pv!6M~++gl9PnjyF)=pw)Z>O+xxAU-TW#?-bY8P!6ZNp>^s=GiT=TWhz~ZlB#zyYKCOu)A&d*zT1bwyU+5+B@1ewr^(N);`d_lYNYR
zf_=h`o_UunPIeyjZ+`@{Ap?a$l)Xn)iGzWodPDw$Zu
z%G_kFW#O_OvK(2lY@BSiY^iLcY`^TX?6&MT*(=!x2hzdTp^<~i!P}v|LwAS%4hDzO
z4s#s7aX8>`(cyu^J4dOb%u(T}act`t>=@&i;;3~rI!<<+0oIb}HIIh8w&cbe+-rPEhV%bdP(
z`qpWm(=n%uPCq&Q;`FQ2OQ-ivpXFk?om?SrB5xu0led?5mG_co$c^&R@@ev~}JDQ!y&S2-Vi`dodMs_EAfIY^3
z&t7J)vv=5E+2`zg_LDPmmO9Iv70&L?9?q?t+c^h2cXW<+PH;|j);i}nmpNBDPjQ~@
zywG{M^IGTa&ikB?Ie+hb+4;8fL+7W?e>i`1u2F~-j6$wZD!dgP6|ssG#ZZMoF+wp}
zF;}r%u~l(UaYpfz;(_9=qQ=G6Md9M%;_VXV677=UlIk+lCEI0-%aO!h2?D-Y+J?HzH4)QB;s0QIK7bugir*U%wdth!`z@f)o%BZIM8S
zzki&yfv)wXy48~(717OFqMNzVZuN}#21fV>`1*A-!f~vwdwl~5bu{|=1q3GuED8Mu
zL?#M3i9%9LezuSr7$gK?(FxF1QeEG^{{AtL?iUjkm1OQXDW@bq8z|2+XiJPjm6(X6
zVuQ9gr@n`v;1nT}BJ`6g;HepK1Zps7Q|k#78V=*z_w@^njC!LX;$LF*
z3Lq5i^)HzbK5#y(E7BJf^%(;Hx((^>V<^ZlBql^hT61EoShN+3vSP7TteX{!vtnJX
zSa&Pd!-^$Xv3M(%XvKP3v1BWjWW`dfST8GazZL6a#rj&YepXB?XrM(Lb>s_r9pR@nPYA7fdT4cmA`8kH6ObW9
zGOSvjAf
zwm^e_K(@7wocg>RnEZO7f?iPgw3K+Ev8YJBxoN$*Be?be!zXke5W)w(!Gb{?R^OD+
zO+ZL)c9Bj8H!b;@`V4bYltG}&AjrS~ROAcsd?Bu@GmXJQ<01Y9b$o+aGKZ=7$KAyVI3TSg2>=(
z9o#PEXQms9eKNELqn5W>fm)%g85JlrqxCJD<8_n>F=!9gWGvF|lDnA)XcB&|JaVbx_ujEcidf!2jss1`Q<0
z7YWQ7GAwI|(KsZ_IK-$QGOW%3mkoh`TC2%Pw_+JqEX#`NtXQTM%eG=UR&1CR(_67z
zD`v1_1y(HIisf0cLMvvpVntT0*oqCeVr5pW)QXi`u@Wm*VZ}yRG0*{v!8DIa!we9g
zHJWZMmubZ^tQg4J(n^*U%duj5D+WqnDL2fDx8c=74VkVi9~07$XEBA<^bUe#7M-80>2f_&y=#
zZw>;?fuA`LdI$^_dI$^_x(y5tf(q#K|x6o#Edy$^ZBr>0Dg=0eA%1yT8i3R3Y66=dWaDu}~3l(#DWzM+B`eM3Wp
zHh8AY$-+4JhKBN#@u_BNcr1%2Pe6Ebg*>2;4|y8-P@qys5a{DEGsQe6^urH;xy8Dn
z5K^EWJ{062AbeLmVrEpppbtJ2`rr|B%RB~z`}#(j2@l~phZrvh2+W*^m{|@nvuxHz
zU$Ywenw8QwM5yi;Bq*$3kWkVuNT}`?YRMB?@e2~f=ocjDm|u`3Pf$F+AVF{af&}sV
z1qtH!3#uCvVVJDqt{=1z-})S$pZc6YfBwx5ZC)n4p8zXYv2J*Yu9$yI!kQp#E;IQ@
zK6paN6w3T1iQjogMWhK?=_Lk(&M2VpPD)<=I$*A2HG1Yuc-aG%FH6?HMC#P{@4FH<
z5kBS$83hJ-vM}Rj1rsK-Z)mi+kGhtGypU-AHJW^bF3U`#IWa~^6y6XrGXYx@owZCl
zAN$wLT%GY>vxQjQ6y+;GiNCMXu#lzIhDp|X4HGTp>Nv5qFZ9_k%}jE`6l-gBiGt2X
z2u9D>&o2jFU4mK854$cu6CPL1Xr2|d6sV)c5;e=(5~*uB8{T{}=<*60ux(AU%F5im
z*>w1Z6&HeFEh7gr&tQP;*ch`#_l*=Deq#b+`BwOMZw&SFv?f^_jwvy(2X&J)(#&F{
zFk5^hg}LKvUPPctuxlA5MW6>66@Usl5Gf>rl0!mf0X&oQcG8>^0a%uSf2GNq6b%^P
zr8Oa}4vXNiVwfr*r0lxH$WXv41bZqpXB?hCu&`=>n>Sz=q
zR_cUULlUjo4RZus3%h&s(}_7+&#MsUX{)ElP>{`Y0-Xw(7V?Z)<`Tv%ffI<4lskc1yH0euy*DjDllh_`-WKK!V_nGhPj}X6``2WsHKXP
z=>|Db4RHSkcw_@yV7x&*umK(+;8t-o$O#tkh7zgI^f&Wh;kZG%xsq8tLG|eYX4cH<
z=DL}BaIS>`1y#=}D9E))JT^2owk|4E_w^OthaK1|V_@t5?7iH9tx^WI{$JnC6<)a|
zM)B8kk??*Z{}$6+v;{7IhZEz8pW(XfKEc5idInrUj)m*PS8xgTkrKltW*`*-*Zf1N
zOzJCY9d(YnE=m?vh&GE(h%Sj+!L`*u@e1(<@ow=!@k#MDiCW?f*H6iC`BWqs3D-_D
zB`YNxC0iu>B}d@G={L!Hsa)zJRZ2aiUeW+*u5_GqnsmN&nRK;ut#p%gr*yycxb&{{
zR~rRfEVZ-=w~4Vyut|aIqzoIQ&0M%jT55B_<|B-=dO
z3foHC8MgCm7uznkJ#Kr>_7B?+bQ3ySxF|Bv#c)|vNl&I%(!1fB=nQ>Nw@_1(@t)uw)3*G
zSK52nx3>4SkGD^^A89`eu5pgqUly)!KFSE0ovaaD-+0Jc$h>8Nvi7n_S%2AZS%qvg
zT;xoX&66#c?UJ2^tDMKOx3bR;l!Lv4i-XdknL}#_Ux#pqM2Eo+!yHN-#yCuOnCr0E
zVYR~shpi3=9WFWCba?3S+~KXmCr8SWc64%dbyPbxcWmVt>e$_JpkuycrQ=-3)s7n+
zw>s{I%blZ+KR7;g{N2&yLWT*a4L!ESRu~XzU!fA}tB&V59^PCns
zEq7WAmpeP*dgqMOb*INp@8zW2Nv@H%lLyKp<;n8ia;>~TxYSu7Uo2k%7dt!T2j$1)
zXXF>t76--p=<{>mhHi&u>;t2HXE*Z#<0`jf@e9qiQUT{WiPPT
z*k9mk=LKBtnBZzh;%w*a1eZG+=N8U>&K;cNoKu_!IcGWNIu|;ZIFE21>pazYj`L#Y
zHO}8U?{WT4xZ1hl{EPE_=f}>^oc~ZbE7S@vMTjCo5wGZ_7_7)v6e=ndqZQ*6lN8ew
z3ryKnyA<3gE{RK8Kx6(_iIXt=a`4NIQrAos;Y=lN!zaRTU6?`hsTs7I@t|d#s%1Bh
z8a+(8W6tKS-@9MjJTy%`O*&D^H5)C~jN@nn=jhh-e$)pQqwr9gvNIKURDB-{==#3#
zC^d7{n0T?RJLl%d*>fRW$P>;EyWPM1-PH|hMuEF{U!KO
zZRk5!#qN#6BE{Af6{RKRBeJu~m#$j5V%at|dw`oL$|_j8e#3@E+qP~TonKHeB1*q=+}Y>ru#tW*
zVQ-%v+y6Un_qI`qTobOzkBP5dUp#f=m$btj)$EO(^p4fAk1j5T9=r@)dxP{Hu#CR<
z0n7h99Urzq&C8`z2B&(%b#>=*xInHgH;LiGIEu63+H>tca<;f5c6o6cOSdr{l?gy2
zC#@{T;UjQ2?}Rh5gWnp^?*8li+NnPPBmQQm0vyk}zE^$xdQ
zQqB#3J6^nU-pa4nxo=%vnxii+9;Rl?hm|f-vim#C{bBGU_usA_c>E}Jf1Fwts>}z4
zZ4C;$SHjKaI|-NIT$ivjh{SdL;8o;26osts^@&ZCRj>|DrUduRA+p_(LxrOvo&)246MY}>hNcwT`ZTs#3g(VP+f
z`NE)0U?Qp)dz%)EIgR8^^)4}9Cb?!>e64yhkG(PN5_4*a>{D+A?nUFal5!=(y)_k`
z>2tIeb*hB-jKV9?(%<7R$XC6+|3~XrY){;|YM!D?n!e{mHJkEQ!fstM^Xr8x-B&EF
z91gQ~sx*3%IB9ZPd7@j|fvjUHcFw^KYu0|PW_u2Og*$;1PG=o$sbWjVaU&EZqrO_U
zeEEV+U#PcyA{@he{=U_x-OhI38K`0xp^zSUk`fea0a`t4>9P&(>z0oy
zD=lRd`*4!dN5!5ENaz?)t5#IJA&yT7WV
z>6HQpRKJ$8FAw1#-DM=
zewerqCxRoe(ZJWPjgNhxvtU`>%5_Tz9c}xQb6wEE74+h>BO<4ri^r#eL)|`M1YZIOfMB!boJ_Etb*ZfzFH||xCt1P0~;CB
zUA*d^vYk>Dtkfz)X)q|uPtDu2_Nd#5K?j1>>;T~Fi>e|qcLs*Zuc}CbzpN<|
zGuViIVHhSb<+PM}v$7Oxv8%~_GlQMhD4D~n7wuD>00BuxD)FG|DC`a0Mo73|?x+|$
zO0XJl62IlvNoD>t$1MLN1-rSIfGFQ6a1iztbKN8yrKRzoE)Ap1*lZ;;Yw^s*U%N5Y
z)h%h5SDYsf!JbkEdvRXS1^32O^aJ*~LaccE`YnaKn2HonPoKSizvuTp-aV5-LQW?>
zQ~w~+k6eG~;F1knwk*ycIABzP9&GX?+7l=EzWwdHU!*W_2?{R3ODcPh-QHq{_sk0=
zR!B?*IDlqqQQu00CGQzE_K`4&wCrYq@(zHHGR#6A3bq#jhfv$
zan{6HlT>jlPnW%MXP$p|{NguiW|tDr$AR<}rG)u-8)lE%^b_Y@eOB7vIC_YxkP}H{
zV}QK@*cQ9v1$aIfjwajyCVD97Uxoz7XlX`I-vl}TfJsa-ru-4GN?|7)C*~X^ksJ|;
z2{8w{<;2BtN^TIZXAiLaJ|`DrCy9(3qIe}`S7HZR%?@3+Va=HpDmIyGD`972JuUm4
z#^E&QYdXk&V5ZP`ASja`7gdU9Z!!%mV{lviR5DA-uA*V%HHzWdO4%uuoSl?&Dh&+d#BQAIVJfyUs9p^hKPlc_*vqY3zl?xz
zrs|mD*17YyZuL0dv139)$ByR{ZmFM$3XE&ktzW%*{rXkLg2IyG0yWe3nfr@tXK@qU
zB<5^uHG6m_*F(X~zz$ekl1T{Ftt4p_olG
zE}g<8><`D3gJXY`T-dYasG5CnP``P=9F@%7vRzQuuh9&4QxaFf`2Gl!044%%&X_h~
zCx&x+_;bvpNrf*IIfyfdlGsu1_O#lPL^4!aZ7oQi89
z#z_*U%0_Yba7Nc|{WHSUTq~@XFu!D6h*Gh`GPpv8+2~BC$KLN-iXT1b_y8}yKekN1
z0ZUlR+s?9SSFP{fYHUV
z^lBQ<%%IujjHyQPB@ODeveXoHjAka@>|cugcV%3^b(7gBS+;omh>Fn@O4STrpKw0m
zN1RGCxNOEZMe9vtzF|%&Ie+H)4J8v{S}5igLA#!~H>|okw}kV#A5_QQNZ3z$NpAU~
zb(@xKJbysVzGJAJrilvd&lPa~FKN)35u7hB5QFi^;p|{DoXtsbHfO+g*no3@DH*o`
z=xU8soC=<*3^co&nfLPF>}p{zk88wSe6`~-<9B~&zak7+4<$8)agm2#C21#X~Vilao6o^GqrtB$hE2Inf+4OXE0bn9~!G-
zs=atqx{EiZB?@d8bCGk}s%o`P%m!2ju3>P25_lzUV?R?f>_;){8QdBULh+6jU$C}%m&GlB@;%D7(J;B(hBL}BMoZipgQupIL816)YUW*X5IUCzhjOW
zRA3yeVooT}Vqf{nPl(;mz#3lMxpY;jX#g(Z9LrV#y|A<4M6k1&td!lpe(U~i`Kpf$
zXOQ3?uqOjm@W!zF8H*XyUd#|ol|=EeTCR~Y*Y<2&^DQvw#kIgK7(6BO?wcnJ>@6GN
zMhr}J684)?8D8)A#P21{5g6TM?9L4J@=D)({PclcPeAoyBDZG>$9{9>^uo34)_yg-
zXEN*`)iB>WVF#~L?6)i9_P3bsbd$ls4?o~W2ltFF`c}pL%smq`XU1O3EXC#5aNsrE
z;TqT~=72hi=RX(sQhvQ`sG2z`tj5*Vifo#B_x!cjbLJX;bd6wO!zy7e#6Ma4{oxDp
z6@UEx2YVP_WG3oivss-v%9i~G8B
zt5#iJxA^(Z>l;!TSRD8`;NgY4*M8@Aw%bnTCmNrVFl(_5n3W(r*bNVcAq`M6cZfxJ
z9kqyAc~i-}Y4InPy=SmDaS?|xPu|>>AG!-fX`Y(7hlkyhug4c*G#=omOm$NQ`?6K9a}REQfBNZDX35AEc)}9-VpAHs9lvDi
zP7~gLfYStc?wq
zj-M|l_CCW2&nSEy3eTUhV1e6$84G64S1~{1R#?v29B#KjEn7WV?J4Q;)v^5R?y`I(
zY#s05JQ^g{x&1FV;TvYtn#~8@F7`ac@Q;hr)bAxsD=9l{5M$cNdBMaDxdzhWG~zcB
zrf{Hpww{5bBiJej;^VkIm|N;a>Ro}I+XBXqnTsY&to`2#lVdkm1O38z-8BjusISLEZW$edP
zod$cwn5uDef>ESu#0)V<$y5g^%U_YdIKv7=FALG
zT)ZEL!!FG2Hyi^~O!^yhWBUza3tmd$@0d)Q;Nsy48El(UDl;DFt^$($fI@@`rhn{Y
zkZEfPd$aPxM#hv4=1EgVFgFM+|Cxfxf*lFd0_;dJKXU`s+)VBTv+cP0y!6{e^P$WZ
zJdMJenTfemQ~}b>QOe*iQY>LTSWFLWBLgv0E_T2
z+**p6d#6FIT|-ZEjCwG4M1hALl%Cp_9UYyWodVDLu-2%!rko1uXK-m^OiOm*G%;hE
zqrgf?hiyPgOEDL~**ujn-(%Mm>NVI8x>4d}#jO?{;nBRmnrm83z}~Sl4#(cTq>TJW
z9L@y@&!|`0#}iKM+3>*NCg&^&5^(!g)dWY>kRTQVqUZm(XldLwbNOWi;88Y>qC2-(oYkT3S!GBwK(h
z{S~q`*$VxDek9w%H-UCHj?<(Q0alJVeS
zcbiNgd!S#y!z&5C$~6zZrgjI-2j8>R;3t_%rjV&*8cHL3p}XKq*Bic#v;aIO`;q;@
z=W_t+jqZ_yz*F)8IT$@8hl1ao7B0ln!EdcEnMLZrUo0Cw6w)6(B6HBMq@El`R-gg!
z4X?+ffy@K%(*m-P91fn@#o$*~LY9(cWI1@Mj^I7|z^|_oF4o3?cjp9hJb0~625-E4
zatb+>oCY3fGr%Wz7CD>zlKcWZM(2|A!Mk+<_@ylZZ@-1$^S7K_N-iUpkgLg+2+JUm^&SJn;u4bPJ2$P46o@-leET_i7&
zKay9;AINLu&*XLT2Kf{D3way77SoA?&46hDyf$UjLF3EOAzPOK(5@)P*l5#aq=3m)Jk__K>C3OvDWC|mHr
zW55s69{k%KC>eO;If8dMEBIFOzLu_3W0XO;p(WsH*#tZ-m!ma;_oaqvN_l`kd~>1;
z5dr>!QQ*VZf@%r=@onHsdk)}9>kfXa72tkd4IfO)4Fg4jzc$uS{i7@No)CM`cqvR2HSD^1v@~
z7-gXHsX}TvRX`O{M(|TCrHa9CcOBY*R=~Fl;X8)#Rl^Zf1zJj#qakP^>P(HKMpL7x
zG2qcSj^e*PFoK#$ji)A1lc=fGWNHfdInJVHQZuOO)NJZYY7X@Uc;(Nd7Ep7kMby{S
zLbwiIMlGe5P>ZRR)GBHPwVYZ@t)W&^-+*uAI`|U725K|)Ewzc-LT#hAQ(M7bau2ne
z+C}Z8_JS9)6h57h1)iU4seNb>T16e8_EQJJKXfH3gs*E2r@o_(z_$#JQpc%d)Jf_D
zcv_yJPEp@e=iv(aEOmjpOkJcdQ9pqHSyX2^%HfSx`8I3@n{r$J8&`@10JU1
z&_pyAJi8~MDb!8s7wR^6Yu=&mQFp2P)FbfQd`LZ@ey5&Lzfq5=r_>AT8TFicNxcSt
z&OfL()H~`e)sw$5=UvCWVJYYUW+ofDSPf?FF)|S;uq5pP1OEu_&I@459wP693iSm+
zE&<^HLqP0p>Ya!Zxr;i9`iP1}6Ge+f&&5vSP;l=&Eg>awk_D36lGoBO=|E|Yv|PGE
zx?Ort`q`$bjlWHt%_y72Hap?I`?GCp+d;Ouwv%kv*j})GNjIf+bSbz#&ZM`}uNi-4
z5R=P{X4Wx3GgWr1oj*7$_O{EjTW5E}-rnBFevJKY`!}+FaQ{3;HbZtz_SAuZ+vXDv
z_Z^WV={OATjQ2QRgL~q)j@3?5rvRsPxCh?pbj0bb(=&N9xcAk|XUmt#zm=bp|0;ja
zdc&=5SGdnz%&uh*vJcreY^}2c+}RqOM?0@~KH_}K`4im4_EqdvJamz|>~lHpa=sDO
z$f1#{QKv>d8>KfYX*8kH_C}{0J!tg2QLU?=Yq)E)Ym)0Q*HNxBT<5!fF{-ty>#9H0ta^_6
zg2rB>)U?v{)nsZ`YR+k%YHB?g4-b!79w(bNZ`!eGM$;`#Z#8pjmeOoevu~T-Ywq4$
z*}QG@$mYi8tDFDQ{H`bA>FL?Qb1Z!N;;82jo=;nNw;0l5b&D%4Znt>XQq(fM<$#v6
zTE1=N*ebYH|5h_vooV&Hb??@ht&3Z)YJI&m*T%I?P@9A{1KJF4bH2@ow$0l{w_Vis
znU{lCxR=3egx3_W>)vATmfjt_w|L+4aq^kw^OesbpASACeN64xc5T}YYd5#unRXw2
z?R`7?=K5~&z2oQR7wOl_FWYa5-$uW~ejoiwe;0p0|9Jnw{)PVY{8#yJ@jvK)HNYu=
z4e$zx378nLBjC3{Hn45rkid>5APE`CVYMPk@ij6$F(2Z{+srf+CS?M+###O
zmJWwIoaykYW8;n;I_7sA)3LJS-i|Ljwd&Nj)5K1fJKJ{l?VQC&aV
z%kC~$BAP_>j>w8w7;!t|RfH+hJF;Ko=*Y2=<03ajc}B%Xjfq+lbu(HPt&5%?eKf`~
zrghBMG3R1l#0JFni`^dQ8J8V*p=JnpGm%w
z{9E$to|2x*o~?RD_MF*se$OpE_x8NhvnItkrAD3v@
z8RIe*WUS3NmGL3dGgFs&AoE3Ltu8`0MfbICm+sH3_E`y8*;)Fm+^h*%OR|<_?asQ8
z^(d<<+b+9Vc4&4~_K@r;+1s-(XaAbxmeV0eo0FF_GiP1Sc|D-?1bnfc!r1QwJR+%2>z^e)UO
zEG@h{+-rF7@W|l_!}|`;A3k>Yp5eEKzbx`DN-HWa+FW$2m@4jCY%Kn=c#ko{SYcdi
zyj5ab(zaxFNmc2h(r-$)mmV*@TFRBNWgcaH%XDSK%SM+?E!$M~x|}ZWTAo`zx_o{4
z(ej5C4i(CZwiTfjF%{VrJ1X9eupjZ|$QC1$MwX6TI`Z_$o1>^vZAPVx8aZm$sE4D`
z=vJd+M(2&5HhTXUYE1Z;9%D+zx{Td3_P5Hf$N7vKH}3R!W_h{KfH4
zC-j{#bwVv%%Xm$SoHSt4oyp0Q2TdL}xpeY`$+IUfnY?-O$;l5UznNk?rO}k2DUnn9
zPRW^4F=f`2)l&{l`Dx0dDQ~CLPPL!vKGkb#%+!pj#;LQWE}Ob->WQhpO?^L2GA(3U
z_i4SSrB5rEHfq|!X$Pj=osOnEO>a9re0t>czSGA{-#Pv4^heV_&G49!Jfn1`%ZVA-
z$ut+vE$+h8?rQ8yRe8Bo&8)s%GgFFvs~i+~b@g0swL~VX9>V|cuOy~dZKbNVxOme2
zg#$W1YUOQ`m!g2zzDYIiX=
z6V4c_w>>d;B#?-?MKr!#bDHnqIu@DJl73X;l+UryDqg9?4t4FyO+R3N{@Bh`)i*}5
z?wi%$uTXKb;cbfgcOWiet(g7Lk1iNqm^;`l_JZNPaKLTCUn;qaRXeR@fP{`ddHaOG
zbDz9G|_GK^SQOPZ(
zDW~Rrm=ZUl<j>R!{G&Pq;PBZr*JLJqKR&cSbfRCcT(
z%^7?7i~6*S-Ql;3`E7vJID}Ke8v%&DZ7MZ2RCH-GlezbJ7#maduEn>J0XpyQM;{#3#Zr;u)X6i8RIDT?|`{;nHp
zQFt7_=|#)(_$DZm2!{T3jT{pAeDhqKqhvqeXLIOg)nn@mUJ9hSq1EGr1eln78B-U_
zylhr|SxAS}wlFI}A3kD_1{0J0#GMfO$Wua*Rf5J;K1Oj+UK}7{-Tz^z@7ivhJewd+eg*LDcTt7F#h3qF2#;%HL
zxc0)I;6ncsyjMWL3ll49W}1|eex-)Y!ERmoiFp)mo~m;t`U*JbVn4xjhAU=49>U}T
zR$t;Z|8s8XtG7Y7$#9x`5vD`}058x-v2^u3{;U#G|A0>afSf-iG6M~YVCsM!m7q{A
z_>=-p)+fPR1XIFeI2zu}3v5PnGjK;xumm;C7m{X2ETf^;Ufw#^SFJP!Rh?9@dn=cN
z%2nQ=0ejA4KhuWlc?$Lr`iDj=kDZpNt^a1SDt`wib_Pl6tBX}Rgg%{7gO0#_K{{t`g
zgqFSoUTW+_DYSvVkKh||g+{;xD%tNU@vj+5PE^xE!5;mGhM*j?DGC^G?D1We^
zX-eF-<}|QYMBfk=%G-FJ5^j${32uT8JOvfF2^#S9XZm){$?^0LU~%MpkyCr<($Bs5
z1bI3R8M=2U_UC_iTm2Nvs!l85sGNO^vy|*3{_LE6yqE6C9=0Z>LIF~Urlnmqs!%d--eFOIcaBFCxynG^!11*$4Gf!4HjB%s!u9_Qs
zL&ACrP2K(}lILg~Tz10|za)Ejq24kp`TO8Euw1fF%?qpyWblw*?O;OQT!%e@>tw)s
zDe>Z(OXF#7C#+9ZHf6;B;p{u$qPV_)S(lx~QCGuc(Xg3aRP3PG3l{9X8^wlQQLzDn
zU6N>G!D#GViVZ6&AS#wvv3Etq-mrJitanM?_s%Rx%6tF!-sks=%+B0;?m6e4d+sT?
zi>jh0PmdVUpA{2KG)Q#5x(73<5CIp*r3}5YkD0%I2h30jqwq-_eg
zH-Ub5K>$@xVoM03@G
zZ(zPKfz35gKH(tEHJ~1tsR{}}Y95ySE&XVZP7;r@2NVftz822Min!P<$M^FHToKGs
zvED+pw~1s?3Y%d_RU21UCO;=YUvPNNCzyCuxo;rJ$f{H~&<<9`VKmkXGB9`Jom2x`
zp-cp~yKOWnUbENQyKEY>4!tX1sIrv$Wg2FGT)_HtgnK*Wzkbj?`+YxnMGxU^2|)nq&XcMh=vn_OQ0>GOlS5cUo~S
zyb8PbDwRq6seF48vbe`_S$aQ2>)R=zSODkmxs)>VqIg^T0$$!zmD6DU<@{J{s=2tB
z$<9B{GmXX5#G|=QsT0c=Bo57)*h$FT{3WbO!6gP}Kj8BPmw4P$mB&ILaL>Td-rJLu
zB4$6?Dt%nNI(p73+mKNa!%oQW&zyqwnUANV2wSrU)0qlmG6R)XhEgf~{RWrd92hRr
zZ}vvbM5*gGSz4bUON*%}9@hU?c~G5NWr@q9&Xb&rO0BbRl$4b~h$_0B0LmB8ki9=;n2CW*7$lK=uq!kRY#eI>qxvaq3gIt}
zv)QeXP+(uU!Gn2vJ)){q_46b0v{8P?G{F=LVPgfTcCspLP*RQnh|OmOO@x6A(`V62
z8Q)78v-bCI5YW9pOb<7qmN#llqVK4HO6$BLn1b0%1;DoO5p^pu{WD&7l2$iBt4n?w
zBug?=gxx~lIKkM<9!e$6$*1BovpnlzEmpRX?_G=Cb;2suJ=7iMX)wrj3?hS>6|Me}
zYG@`3RP=rcR-@Onev=2=hHo;kE&|p@HiQowI5o&(%46CBuBO_x*avwu5;J8*#9^@*
z=vn=7xtd7p!Mz)ad)0%U(;NnMj#JH}u!d0*=%hJo$HD8bJkh@JjJ38x`s4!HcCrbu
zhuIE;Ab~yFF2nHZ5|~6hC%SMj*mqOSBfx{>4P|6gJf}UpRLXRbUqTlD_lHZ>#7gXg
zxr{tQG{xiC;0?PGdnrH&*T4^*hlfj}MO2_-+OfuI?J(T>`m!C?0i^)Bf1S}=}eX@G;nQ_(-4yjgY
z>SPTI90AUWXrzoO+@+eh%huCpWw82*iP^;ftfn&AfDtgsE=obr0_?C)!{N-n)C@E(
zR-IMqZNoYHa>BXYPb_bqbz-GrZH{}wtf{l7MtQPd$;O(~6dF#oV4d{~uIX!1Sr%di
z2V$;uACNG?ruYlYd#q2Vau(C*7}e@rN56w~R$>P!mV^uMIl$$M@IC1SV{wY%?7loW
zIX5W$k#%A%Ss_6FB6fp|)QSqRRw(x)Z05A4xEYb!nsO9}qpF?ET9d=XH`nA`@nx8(
z?vJ^@kBFLGXpe#
zSu-sq#lSiUtTCI-8q;k4OiP%Zr8vqMs3WvW*arHsWd3YghMKAgl(NMh?JrU(St;B=
z?#P1tMOa61DqIh#Q)!bku-42|{o)RxtY$x|(nsBjKay{rAx>gJQz7|MqNBTkod8!m
zX{1D@o56(E>E<%=aVIn}l&`L4eMR0EFp(y`<>+{B7as2HIPp_^5_1%1}5B>jpM&TgxE`B2FMh~3#{WpWlPFm3%hYG
zXo*toCMiQHPG`+1n4kFxu*j|o8yd$))sI0+WtCYi16v7`abI|$u_ifQ4LyveOjWO=
z)Zqd*j?BJd=GWseUgG?u9)@2-0Qw98=rshO&+0+gSES;qJfhakIo_&owS%ixn;%>G
z=?kY;eDmbA^PSoyER}15$7tPpi}SM!a#>Pq64R>*Ee#}M=K<4Kd<&ig7YgiI-orP;
zf~wvHy3Y>3gy3><#wzTSg>|Uhls8`+`rU^zRE?-tZr;XK4Z}kv)P5
z_7Mh_EHDAjCHS-Ye)RvXtD3Zw%Gk$H>FhBEN`+e!+}K7^zU-0euRPy&3sU-jxwx@4
z5j0F@!V;0~VZtr>W;>b1u|hW~>ZUICvz|s!y!_-UMTpSUV#1w^w@==ozT(b@BBL#L&fU?ae9+S=T5V=)
zG>h75cUUzjqT$Lg3mcT1Os(G5yo-}Z-_5;67QZD66-auo2W7fHk{vUdpcFxS6eB7{
zV3TV!Ir^d(-GdMM8%($_vQ%C`w;L4}K;2S|yHOP6$69{_`FlHxyrGR!pITf69xiWGg+ax
zV*eJzk{xe4m>uTwUrN5lG=WtT|_>iS=JQePRHc#=f
zFt<7l;5z%woyK5mh0b69WJ@zbcR2u{T%(KtW+_EDW+_E($;w5JT1x(0)Cy7Zo1a_*
zYK`Ei#XZ!vrXJiE??kPLfHYqVY{-FHE*!P6^1VcRw2E3vksQ=g5;$%t@wvzamhf-f
zQu61;EhUa4mlB`FO#pKc%VX6cXz%a(z`+zV~Nku1N4WOp0G8M4rY93&ek
znCgqqMCgmtVeDQ4`s0;Q4nuo4397NQ(jva5k>6R79>WhRiuL0~_M^yT6q$k|2h<{|
zPyDD_oY*Wwtp0~6B}Y+uE3_0vuCkOI@C$!)QJWuopk|zx9{P#weGea-6u2s{D>dZoA~bAd~l)#(Uh)
zd%TSvr8_QUVB400-=P=hI`A=+nq?#i-}_x!oZ|Ngui&JndZ^i=1zs`>HsKohZ{5v|6MxFQn;8o5
z?~kGrZAIxoKr|b~yFcKC*-YU8pTYrtKm+(0)nft>hxm^Lo4y*ri>l0*sby
z>YO*PzK=t2B8T8a4#9mKf;3*#CE|uQQq+HGOR=S8ya1Z;cY0+-9g!j(Q`MK3KCAC&
z|2~|n9Ro|tL#J}@2m`JbKFahImf%tZPR~B1GfF9^aMw^vwn-@?MXSl5x`?{_1
z0RO#1P<~X+nEV2cTuKxA!7Occ@q@SDDi`MSs5r@Et)L%BJp9KhV1mX9*6yo}{am1*bH}QFHB}vruxnq@KHdJGBLKVJ<-};i`=h
z#!wsQ2v;g?m72I(DkV$3atBJiSsOFygk)-CU{Lk%WslWyqvq~c<*2%|u8WkSJhHbn
zOWRKI*`>gdR@zqXIcs4=Ya@AU%8O%o4`0l!?5^V-WQEbEuGmQ1%BEzM?uQ=~+K~ttU1Cz7~E*o}BP<7LQJ1Z($x4Ix^#sVe;ohJdjp=(}$ZwdV1
zYnOMFz2Hu)B}%{Mg>V0>Y^yV0OW$%)ssahG^BoU;tVOLCF6v;T{$B4CbjBF#L8`KY
zj}u*ZL^(iTpRT|R>`
z&EZPC@33^z#}y)F;@XiR&IhyyCxx1_Nh_w(_$ZP5pBbg+ui*53^ORA##(V_nNt;f1
zDtFZJg4ZzS)4)7TuY>o<|27J-i}5g^>qdpNZk99{#WvyLHWQzZc=e9>jTatGO(yO
zv8Gq9J6`72Ktos+G%sbPHWT3&%Xzw#rTN9ef6>o?jViZv7C_mVstnpOh9A@Ls-@IUHQ1D=YQ`1wipv8jJLsBb#?Z(?q-Rk~WFA
zq)j3u`;b}6`U#^~zeZwUn{ep=`r_*<^~+8&|Ly(%0Cd(UJJy^N2q0bf8Xv%Ds-&(o
z{m8ph3ejOR&e_w|7!iSNltF2MFa#V~HyEKz+5q=k!BmOs^no;}!&iDav^jVT32V#l16(^J>8D8D
zJ5ECjxa;#uc+NSD&HEXchw*%}fyI%plFxW`j{!cECR$FFWduEFi&pWUw9P(%`@?o?
z*)q&_VO+OUz7{E^Qjg9p#vslR=@IFHd}L|ci-eR@r#5~K!*kO3&o^pd9#gkR%w#oD
z2N8T=i^gS=Rshp@OJ{gXTfi~!8Jra*Qpk?IurXwMX)j9MEO}SG68*psi=Tt;f8rt2O1a*)#b&w
z04>92n;Zbz+jxj?w|E7aJSBj>^)_SB*p4u$s{wpe$pLwXl@t1imvxS(3m2tnjHste
zQ_F~5qD3x1>XpKfr~AwNV~#M@DMRe6iWJxA5}4(4&MdEjSzZIPyx~w>ONhKQRs=+g
z_OLrtz|NHkM#;QfwTJ;t#9lWH?D)Uoa#-;az
zvfVH&trTP;9LM(Y@fBCg+*Ch8Uiv`n!`H8virGuN$@VUnb0_4T!!Cc3Pheusc~U|j
zxw|$`oJKepl;1hi?40e6Ub#HQTZXjq91>M~jPAYiisH>VlgODCM8gI1XGbsK^3<0s
zc@~)J`TaXHUwb+#dXPRRxYS6oRv&~Rv6{Nzob#R@zpwOxUk(FSxEzK8b`5TacjHoc
zE59A~5^+r|<}C?6sTN*BZ*jJm+l`GdOWja@1IjPqSCtiSV&>DSz3U^$8V?i$*$P2V+5k{owmMCvahdD``N~ha_E83Hn6Ohv
z;>|H*(~ai>?^?+7F`_2aV(;PQ__fIz{VbMM#A
z5}YUBS?F_XfU=FMn;-s=#)P$2t?T3fRNfj;MO<0(MkmtYs`qRFtWVXG@HXuvyd~!(
z9FThwE;|23`UlTE^U^p%&g4=Q7jT7b;0m)qR~X?nbdsWuP$lDBP~fU)4F=`{D=$TRZN;k6(>+1nfNl
zzNv<6F*QXJxCl`Tzm%#}peYq(DJe>p!o3gxZLl~VY7%UHHb)%Z@5TEGct2?{s%SFq
zaq_K}A=ZH2XR`~P1lLY>;cR3WxCb8RB&`)Vn<3S5ZL?
z_(h>m3~M34iVC&)w$sn-{PaFmZL6er|{Ylymka?*Pgvz
zSsJNZq}JQeg%sT`43u-G+Lg!+>;inEu?8}
zR_TZ}C#&O%)@Vsc#xs=%dgVJxHp4z*K98n|yP&$vjp`a$O@Ugo`H+-xd0k7uy;3*a
zB6TzHd#m4M-EYc0Jso;bRs0PX;(prC>p=Sn+0~uzFkBRZS-3lf`x3*Q@A$|+0X=mu
zFi<_NK{aIz>zG?fYiLDBIl!@u7k>X)ee({w**4jp{0-KGH|qjnClbr?8LRQ!sQ~?a
z#c{r>NIo#JQpW3*IZBDnumKk7XdQ7n`%`^Zet}lWll8|B9}zIqVOCK1bmK!obQTn1
z--361deQX76BtA!E8}oCew{{|smF`5!u`7jDuvrOS4u*_Bwxh0Yg^%gP^zqfl~Tq*
zsd2_=YARpO368X~Kb9#}l$JYpyfm2T146FsS$SyC@L`qn!(~mLixb;rh#d(Z7e&zp
zFoxx-N@EdMX)NL@jl}`5Vrz--g9A`|v6r(t!ZEcvj;v2RRg@?QF
zF+}C|1<5W>Y={MD0v#s0D8T2A_*tqReq~e>b$G*~P-!J1(W!g&Q1AzYxl=dwb5tb`
z70?o(Ow&?rsv%K_~kMc|_RA1#nMeJkDpq14wiu-5q
zd6!n~9BPc0
zztu&fc|N*zPa$pq#;VSJHYk4%Z(|i;4ANiCVfu;>H2@v2_R0klBsQ)vk9Mxd9N^jt
zPWp>-)a@&xs0;fUHHULpt;m5Dg^6lt8#+KW&2?mgR{b(oeVpSN|D}e$0c9Zg4K&dv
zEEU&UEOcztK62W4%4y?iE^V~eDx#Ss9py%jwpu~Z1&^Qe3UtxRaZY;&%HYjH3qYY5
zty2xauLiC=en#*bbYJiQ!UJMIwC`|g;}FzQ)yBP+cUw(LX`DD4lPl-~=K-BN^m@L$
zI3oE~L9kTmw!J2xo@(qRE4z%+YgXwBOn|W`pg0Vbl+s$`{79t#;=$(omcK_6D!IBU%`Bg
zYpoC-5OOg2042US6+cP~kA|YE{rQxunelO4Hueae5M*V8$_c0yI9mMPxLl3W7QQWX
zt$PgQ2s9D+ef}T#2DLONQtoQO*Ri}IkD0qnkOCVzmiY2GZmaNi$uaqBvHu9DD{XwP
z%zNgT$AjxiOTPx%l{P;6&-ZWnUg-tF{2eQFEbzd|6DzX52HSO!#!v$Ub_wVb%){-5
zI>YTQ{YSXn>0J;oy7QhP2)CQy47WSf8E#h^qXpa@`6HW%fV)e+2HYK_2Hd3par<+h
z*Du~lgK4cSb*0*69YqlBcK!Yr;;BrVG&oUYj1@ZHoOFoZ*X)A(ns+czsSg_xDey}A
z&WmId_pRKJU_W%Lqw1iDq
zp7_g7*)izhmBjP9V{h~{lzN!YcQdx{8`9Y}zHg%x7W9UK4qxfA*ch{9-fGCrkJ+em
zEuO48PsMp`syIFz7ik-fdB0iyeH9NF)nei0uyo7moAFONAU&f?X#FvwXVO2JeEa573u-;xhvK+DYME#hU)F-jYc
z6Eup;VS}+@{SEi+Z>7G3P}29_GH3)hlx7>CWLwUHFVXEAy6tcE_%_4gZ9cI(;i!#v
z^3Rb`(;}@=O@I%dc^kS!%yY^P^GhmhK+MG2#SJv1jH3v=t&DeBF#{Q4_JYn}Bg0#R
z@)c$QrMamC0xk)E$mic&c!Md|fEZNQ+-vmkp;K*Je+=E+%hEG!MDS#r)b$wQ2F}O0
zYss3}*d=OkCkCVNhu{8N}!3W=OBpIXa~UQD=59#?cn
zuyTVQny+*{U8Z$VBN#G)S$9qual*R)@>pTq)cLCs9S)HeIe0E2fFx~j-}3efkFXs`
z0~VMEO^FB#pR#(V6*yIkd`h#%s4gvRv3^kc{&gF+F0o2ox)`NENCmZMLm4#VV8R;R
zB08g|I=&A?HVo6PCzqQbI$TsdnoF1U6p>(l3?r8b;$v2{JVO{w7vbRDgEbAGfo}F{
zv12LPeF0rUxWU@EZ0y^aA`GPkDdrrDHrx|J*f)mv0?_$4*&%6DB!BUO`6-*zh>98!
z#`8Ot&rpe$US=gRoX1%7cIjGBIOhtj{MQPFbbn17au&#q1uOY=(R5BvRVSPg17!9X
zfs)&vM5yHERL4D+ON~MJ@!lgt5a@jBa5b%Lust#lVEa8W?lo~?2^p={a?ChAz&Xn=AVRq!c
z49?JSnj_i1(VS7bzEPU@mL2)Z*p#A|=x}bax4Fe;AWZa$-J6cupa|Ph8WFW4zH&L<
zG;}B?G)=#YG}-PBW@qBeu1-k<4JwUr$QTwd>@F*6Ac6BcWqc05Q^uQf`JFO8FTcaw
zdW{az_}wPX??$Tp4r~g-gK?+E?>J`zzw-sZYsUi@E?v#}-H@EO8VG6}cXB-5fs^Cq
z0bQ3Pobxs0Y%d3B)i}@war`Yx)tVvJH($Cuyetm}uRorLr*RDKn)r`;U)fs@Hu~@4
zsF;_(4b8U)r0PFLymtrd3emR7BUzVVGAq!~A_s4#QMJ!7wQ}GKXP~qXEiC^FWnjuGs;O
znaeMeocwYd_+@x!sIXpvS%TR)l<{k50W=sYzf|UMh6%vcla-mAW2!jmK_2elnj9%3
z?4ntUc?_rB2?$aAn0X*Zad{r7J}>`N7J-4zo*k*O&_t@OY&LhMv7Cc;ILBtI9Ml~g
zROO%Al)*pkew=|e#&!+{>dEw+fnH%HjHUn}&`a1SE6+CRX7Df)+aZ=aLoV3232XrO
z_AVk&C*>LJ;VbZxOJGf%iZBJ`>4yKU3`JrSHm;DV8y~}agzdV=y
zU9*cZ*XKAOlwr0w1QdOL88z{x$jFPo7Yd49WgtP@2hztrZiYmhK-;?WP!!!MfB>4sTp
z&QgqQ4#AGV1u@9kN_b;9%Cg^#rYPzmDv^BHC6$%t@~t|y53{Y=D*AV$D3J=+pqduA
zB1}Q}R76BYs9xHJb8cBG;JycnYCl-A{6W>AwenTwn7!+^sT*peO4IkilTor<*)=F|
zxeb%2LSx>E@Fl7*_NBV=;ZC?W`BTL4DJUU%ph(Yq@v=*evG
zXYC<-AFQ;(3d4>9^t<_S>VS?v44`l2fNi%gpG|)B@NB=vwF83dIww{Z69bC4@1kxp
zf)6yrvb?`}>gHM-tMq$sW?_=0(SRKXCIRB1FdNE?Nvv&j?(C@^K;zS
zd(%9brYJ5njUlkU^wAq$&ve>MeW|=*?dqi~Vm;T7T{C3pm~p`~t;{4mn0Ni?dyAxe
zzS7yh)9@i(2O&&oDMVhOsXh)jt=ZO8`oOakICY@m2o{0J=nTGRyzx)o@~_m|EZJG1
z_QxLnO!LIx-(n+@EnBy(+@3UGX`n6S*(Tw{>e#jG!&VL*JZ|imnYj02g@oI`Upw;3
z@;b9;<0e0hX$PU})3NCnMh19RYtV}sac5pZbl*pPD2w*7$uwEY!W=?im(!=B3QCHJ)F)ESLUptrP?O`r+V7jHQl
zV=h_mk1Rgm;lIU)?y~GiSi5IWNc`ZC&`G1CY%KqxSXyO30FA3YELcFcD8ts}JGULS
zNMAC}{!oIs)aqIS>%#kZOrDQ6K>P)0LWnz%mIF^JM3EBnmS>_HCf2$FN01u#uaev6
zlbsqavjOH(_S&B#R}8nb9^0&bv~Ab2qI&#Tw>z;IG`y?4XbJw|n&412$T
zqd%S;(ZsWM{jM-(kn91p%XQD1xcIeW5Eb13|66r9kpN6N0LTx33X!tC$d?t013=(zdNj?LSTdfsYu#LFtZ@ovzm+Akaocel-1
zw|c8(du-(Bu~R07j(&O~XEM4B1zWDIk%bu6(>?;d%;2GrtI@g|Vqe6modgHn|t78)_TjM8;4jnhP
zy-hlu9)hc@(@QrA2X;;AKO}gpQdrxhgHZ#*)^66udN-=#f*U`4uzb98h+ICl{G+~xC@R&eFylKdQunB2
zTk5v-rmvU0S*c=f-*51M$+nD{u@_ETBsjUOnHai&ugROfS0}eU>IIh;U~F|(^G+jR
zv)hfcPhPY-;bW6t-wS48I00=++_`c6p45$F`VR;h-PegV!}zo~SQ>(e8*;GhQHGge
zmxbI3-pEo$j|^0hWg-I)$7N74UsE6o{96ie0>)nxS#L`YjPVw|F-zGoN0y~=x^1ur
zd{EZ_belHn24&An*)}^)yofcV9%b&hC?~VoeMuzw99SL%d
zjoWQQC6?;iFBoLJW%cgeTR_Iak^NL~Ro1AO52q9-Vq5_fz>qw#G8oDn<)WFn)obZj
z;Y|DICz2DfFV&wjwpnY7bg=L6$Q5gUR*~jz--gRhcKU1~xBJ&{*IBdH6$b0PyD`x$
z+#ED=UDB>i>r>L!g$y1za&&K;oxi8S2;st|ueXXBM-MJ#{W8IziC_aX5v(~)v{
zN3q#w;X*V>bQEL4S(#SzC#47yQQQ-?L==fM5^ZdxPNH;So;(gV5_KN4nw9HNTxNY}
z&x$GI&&u{sgGAOE;>h3~I&8O<^-2K5<2^!2C^gtPBTb_RhU~#SzJ{Yb1guXvi4-o-
zOhiAD+JczWD*2y=^?G#)OtAc9)1lN*I!S$EaT)}mP
zCb*7<;My62Yaa-%M>N5;4HsPbs&Tm$$2!;=x^l@izrdbDa#g#^kzCU`h$_)qm^x~b
z>(8p>`Y86cCb^2&)f5mSmt1c{axLYQTqQ`Zr4ChE&L!7^Tyot5$+bTu*M^W>&m=#3
zbOw@Z!(e~s%yJ~xV~|`CqYIL2l@CtIbzV@e!UnT9=#AFb&yU
zmfG6PiRi3mG@OKthKu(vl?RSRM#K2{#IG3*+4*T}#c6z!AmTY2)Tu?(Jm
zB|lG*lC*2Ut<8NzpV|jII=M~=e
ziO%Z+k*w`G2ChQnJ`G^PN1qSUtV3$%vst*2kfqrHfJv339k5?GPGBhMe;Z1A%{E9B
z`Vk{Z&k{L%n{|QrKkF#D$Wp6A!8v`;I@)mlp&qX-uTS#6s~k=5vX=KvgVFgkED3xCVOKMftF!LDWScn-UgfL+Um
z^S{Eba%Cs%-h#c+3}Bpa0Lic3OIs-fIJcR@Eg4|j3aGVl+~;g>L@7AUOI!Wk8>JDP
z?QI3ytEi@vgee1cKs(bm$WDZzH}ALw_`+y<4^_L?m9%zY?lh30HU8V0&J3+_3*+)+
zXyvRZZtW<4L7E@EAT5$;p4k_xKWmrB=@E;csAjEH`h1k8`ys%gw(Kea{^G}%eY)S;IjdGDSQ2=i
z*~wEm0am*J{h^jvjpbX}InLA*3!G(fyCo@4p^c#b8rIE5{YcY=0NNe$6!
zE#RP4DV=O8aQ>MsZs8hd4rDYMufZ%ZULEw(_`35i2HjUPA={){3!Hn@Jd(5Z{{9@M+}>4J2f-r
z$bO5otAjHw_2v%8qPidWXI^R>s`1>^R%x@}PH|^!PI77r%f8qPCCugjvLvNSY{jd0Av;;X|LGGf>ibr?vFy(4QKvSWkN6S9s?{rfQ9
zAU(DXRWXdZu-?cR4hPnO5$h5*#PfvCIRkpZ<%>*}zKt9Tt7rpk_eOxiuW%4#aR}TD
zmV0rQ(@21h1_^Ax$w^=fn*%CzP~)uQuZq5)k>PHjkxGu^W=#AI;Msyoi+Jfw6~4l6
zA&kCSKpGx=P@bHNm(&`S`P?~fdU51jRLUzab}l^Tm4luO*Li99xiG;?Lr^abL%s5#
zs8=qIdTB^HaDmrXEUn0gw=wPzGUxHE=4<1<0;#?%p^S5wOXChV0RSXC;s6->$l3C;3^cs#5!x*QxTR*g?|@CX
zk1E=(@uD+GI*PQwb*Pb1gES2FMtR}#^eekpTi$kPUrsi+yy@Aud;L8>yA_quZ+AqB
zUVW(9ePaLZgx87u2Mo9k2)Zcjy40_5A-3
zzERLAfA{P#22u=GzC*{w`vaXxS-~SasEJvv(u%%3GwZ))bHzoamKm9ov>_NVgog*j
zB)HywloltdXu1d&Vletly#;D8^NHr4
zKBlV@O4md(w>w|LDy<*Qc273L9_gM8Sy;gDIszrJQ|7*^fU-Cdf=fy=1Q%Wu!b?9~
zSe)TUmuZ(oT31nIXCS%bSt-h&DOP+LNi{Et_66!s$b|T2cR(bY1OR$mL2y9vjK?8i
zsR$i{*KFKhwWBZ@pnZB@{`&h2JvH+Hq6pJ{$Pr~;3VUr3TUV+i<1AL2|EEV|Q)H#D
z_RkgBVMbfZhe(1I;{_zeTWpS3W*JyhA>J{|Al;)fY%z(l$x#aTYRVn*h38$Jzy9ir
zgwA~t#ai*iN7QTJxML1i#==uh2zQJ%m=NF8SD8imalk>;zZz6=H;&5)c|&$|aA?KG
z{VCRidkv|(l9I1@l55#(tPU&B3=PB=04$s}4c*b%*acI9G;bB1oY~6Yry7mg~
z*1!{YX*Jkl6f6v
ztu0N1wkqvRK?4U4ZsG|KLfxwXs$(@d
zngARXbw{bVo_cVrSPE2t#SJqij;Km*gBp~%-X-H*MyczE1KyG`5tk0NS47x<
zkJM4A7i8oA*V6~3v^nG5^->v!8@wUoj*$k<-q46hgb5s@w=dx)q9tju2uR_dQf0u7hqDSWrp6O3krB2J|I~cL
zN(!5)0yjK&3}Y3{c&49aUJy=9I8m%oX6SL-SwhVHrL6Q|j}Q+IR_l)t9?sHn&*FJ+=S(*>tn!G(OK6RJv;2y)F
zy=xC%_cS@m5WM|#=}6~Ck$40WL#%A4$hyl40%yb!?x6Me??X9}ulfV}My_5B8RPTa
zuTJ@FrtEqM8tSWKLolS?^aTj#h+)EW-PYH7#PN&~pO=M@=UB*JCB4=s;msCdXd}aM
zA?c*yQ&~{lWLXLF0F91CG?2;Q*fdIfBY+!$R6pmEuF~l~6efusq(
zy6NPFAj))(yQg|_3U~;3eU-91YW{VKq*l_V#^kKD?P~E2YdZg;D
zohmwizrycAthPBheFq%#-uJ&g@R+qUl@k$1t8v7a4nPzr&EORqC;C%?dQ?%qLB6-B
zfR$Ae>M*>hEZ$(%Gg$%fv78#f8dRx>Dw$N^2CCE}-#S{ATMREMi}hG_f3?aX7Pt|{
zx<$CgD?zzr-X9yk+q#TgdkW0lxp`Szabu~F3C*6!^sRLtxlR`PhEo}NrS-+C^sQ6?
z?o`hDu$k4<8ccbP@e}NUUUfRhzP0vbuFHhovD8yMIqdlatIznhG*gTS3$jjTuKNTs
zA36WpMum;!g}r*4EDF){(Z9S=%X|Pur~v^LzQM6*y`_>h<{Z`X=YAJ=U$K+uuE}Z
z-($j;4#=5+@SX&``3I=|{sn5SJmRjd=h3C#g$$(ly4lZ3uw0q9lM=POTYfOgg|NIrU6tDMdObY~9e9?W&QkP-(#rw@BE
z!CGN_I{7)vIr4Jyl7AM)O<52FtBo^2ko`kf
zb)IKEHgEHqG_VA%g^OsRK5P#f$P0v3%chUD%`|p|*+<~EiHCo&wwpdOqNisIWBBB0
zh$LbK&9bFB}JA#*1#*n(C-R2WApTIoq=sjTRHar2FSEloH7!5;k3m?|ov
z*jvyeRME~c+YgxImSW5cG3G|YwJ}e`m{(xTKf#{p>zIT2n2-Oyd9^xbwFPy|_W2m|
zLSfa?uVdcM$J|yOb4z2`6pZ;NjQPrDjCm^`^YEN8tF5SG{^G5US$%Qitv=qK=^S$<
zjJc9vnj}s_z*h<*cil0iJz`B$n3>!i?GbF;1m>&0w~yD`cd75Q$Mg57k&lPqfiuu7
zGAp6VRq4Vjb==9bh&&OlOF@M;wPLzpA7I^fzh(cpf*tGPFWk{&jM4Y+L9)u=Y(
zsHZk_0^d)d2~Dz_zgPJ|m!$lOUo|aPd&9%0?i*qI^*#X5gXd`&+g98;GKt@pG_;MJm!$f=g?buVb%sqWkifyoemzI
zWp<;?ZWKWl<58?Ka?Iz97L!n4{c94dw6^&y{?@h$ce%XR-Sp_cTU0Y5oiZIanM^B|
z$@j?*g3Vz-3U^e}etT(h7w`S?$uq0vZB}Nr6#@1tFn87gF(o=6Uv-6>Vq+`giA+S{
z^1i|CZR{r+FF!sspmm$B{VLl`%hU#2({~5yyNPJihBk2|&4{v>Iw$ryUv7&P9^v0H
z_n(b@kR2Ki2m|F(3HcV_E&3Gge@8`=gvIK^FY*x@B=b(%_Y7{lM-C!E{)7C|sH{{*
z>BfU|Eod`g$e=csy1mbmK>2oG$UL>ZcPBzH)iUEhHoaUT&P}R{J6~Ac!d=VAij1sN
zik0+a-nS4OdxTZ`6?BJh&B@*PHun1Y7;7w2IP7Yl0
zhv$Jbsw$xmieO=HT`IZPW-3Fy#RO$7GK=smS>IGK0@}+(i1Xi&Nr+$;
zL;Igbj_hD*!@hZhTmn~VIhTjOVOC-BWy6Vu8#f-enAY;E&bpB_5i+_AGn9$gka6L8
z0fDT4=)fdvcZ5)(?k5b$J;J)TH=Hu6KMd**ss$<-Vno<6906{|8``Rm2L1>y^$`e(
zLiwLkxXpBe8u43lKl+VdOv;2*A-xE$lfKwFMVcj7VXlwb@qaeS^T|UzL4mTpJarRq
zrVVP-rpKVV)IfCTDc!JKkam)1ir5{4#tkIW^^h9`@~We#JWQZ!Y!`&?I|OWBVinixi|ujk6=0*dzJdI5yEx0y~47u%M2i6dnw)K96PL
z2KrgJXrFW;dlG+&CA`RogP^b=Q?;qSKAUFq{e&lK&w?rCarm^4Q%>P<&Z7tNu4n5=5VZ`X6;pOPGrs1bj2L(EKZck+^$
z?Y0XV@y8@CJq>dK5PzY_m`L}BROu)u{u=M(?wt7df%qqUZV2L!SBbySV!3l#W|Wmj
zS}GH@Y%L}R+W(fS>XnnA{TAVU`gF214?OZ>tc?OV#gF8){}X>&D!hnIJ96G~B(-0Q
zSR47lF?xt^7At3A-pmNnA7pC+l79-)Z^}s@cj&cuP=0Un2I;?gO!5va5GbE?AbtBc
zD(P1>LbL(tdx7)|X{7%%Cw-)Y=cIp)lfGXr=`Y&ANnf>caFRZ^K%2lxf3r&ZU2{o4
zP$T`o4h94T=RZLD_yO~0PWt{T>EGg{&mRp&JZhvr&`J8oIq4TyNgt5;P9y#N8tLZ;
z=~rg1kK6HoHUWShNR0&2uM5(@u^*)0eNa74`mQ;oZ%GmPpmRwdc0AQasEm8%t9WL|@?EFjw(tZ;EF};jI)cZ>8VxWG{qA{PUaG2HKJk
zxRC~fwa=%52un|QPs(2Woi0t8z}Gs2O3KC8D06hXlsS53uo7!lO)nZ!L;t$e2+xgJ
zqg!=hn}z?ehBt4rh6p2xXM?HtY(m;1<2AyWY#vleima5UH8MRXi%8EIqNN~;CF#iR
z4Y}z~DNFRqHM;Ro2In4&&;a}`{)}|T7YJFw$~C2}8H5yg>>W4Sl!7~!LflpSMxhWR
zErTm4A*4j$d$DY+GE1iny{bQ|lrk%R>|2K)Gb0a;QJFxMkR`IV3MLj&hxfLDrx|5YL!
zNnGrYfY_f?=hNAh54xjsq&piVIG&z#M;2n`DJ%BqIx50lFtiWn2
zpjXUXT~z?aC&5_erMSU>DJ=!RmFHd|QyVfRQ(^q;eo%RU`F=xt&D#>TtlwkVzIDoY
z*rII=Ke}%E4$GPJluhClc0m0|Zg8fh_G0oT(KH$O93i@=z(5GT6mHb*HaAkixAUZv
zhmY=U84|p7RbQLb_o-~E3O58dol=tzY8z5KK6XrJZVGgO&;uIn?M~{pDDx4F35YsU
zV5J>4#4LbZM)zS>jC@U{1(VWgISkPgAYSc<`_WvI_4c#MMl;Lr-@JAQ+wJ!L=Xnmg
z;;f(-giaqS=mjR=mU-sR30tl7%vNNyAa;Gh&oBV5zzTk5C2Wqlj^gG<9S=UZf9&9u
za|b(j-eFrkdi)qer-5C&w6rKQJz91>d~~2~{JPca4M%nyJaW!*M~vs!txdNcAlJqMD_$><6}LH#7@Lm{ZD$RH-;;@=kr)U$
zKM<6j(h$PoOQ~k!p9R54t=V&t4{a``IsaM((=^azF*{5Sb;-&nz>Kb&*=0IrP#Pe&
z&Rf<%pluk&Hnsz&A>|X^H8GlEQ70dit;Lk)k-B;Io_cy6YA#9d4eUpB;wHRhkP7As
z^y?SNz|GF0RUZ#>sQ@RBE1Hvs`(rT(P9DiqAunUZjS*C6_XkRhbdLo~{R!zRVRHz>
z^zZFu5D0JVE9O3*m3KQartMkgf~fvX?+-7O+GMT*nXzLQtTIQ!1l|V9KD@<-=z}$~
z{mqSg9((#|-@(h5_6An1)T@0H7&pf4P=jLjJcJO%&&(v8zeZlvM`m6Z0!V;;EPI9k
zRfm+OP>)&#M^mUrS1R?*y~pl+0-M-F^s8vaXTR`)LAdWbq-=Av;%!_70-%lJ7aF3C
zyRP)DXi`?~kYR!&@MA?!7sSXtWX|nI>($D7BC94@1_LljEXA8&V}gzS#@?7Q7W&O2
zD16tGCvm${Q{x8s`GpS%0$%;Km1ZF22f-rcFR;`ZvON&b9pJ2F)bXQ*UMZ+F<1CfG
zaI)0PV5yaqJ(0S-oTY}?-OTXAfF}lofD$s{;npquQ2%Ez^`%KXm_i98X{Du3=J|pk
zS&%56;7s$f67W#{bxK*RSFS2a=Bk3mBUKK`*`omV2^5Xq>Y>-*zBHWC!5l9
zHl=s6sX5S4MTp6q%N$bfgK>jP!OU?xg~06Qs(dybtmZqgY3e0k27~(gKII&LAB^{=
zGD<@L3$cOFPaF;4tmhYIR0^|D?)&ox6>d$1-S;T>IS&)*pm`f_Jcc)(a(<7k@fE|+
zkm2-$yhP<>@tsD`8*FLWzJJqRTqeRY&?*t?NhpX5g2x6Ux@<7^qJr4yhLDk~SJb5|
z#&^Fz@mx|q&X)?xj^AKs%GV0poVq+%D?H>~tK)pFj`FqY#MkPDy3t{+j$^Gz#PZk+
zDgkD4b1Y_VMrAQWT43!Ny`n@RUZOA-u!_2XC55jG*s%-(xiYt>0+#VTL{N3v4#AJO
z#17{#GU1nLuY6oxxsI
zUM4mezIYXt;@BfL=Ijw6fa}s#Fm?!hDr19je!|-6f`W*25D`I&Gy#z=Ac)cgq)7?AGaF`M{@*)E0G0OM^Y?!5|Np<)
z&zYH>x%bX3=bn4+Ip>~xVef^VfalvaS}T8R@6pw3FI{f(efhY?FBdHIU85sDPv13N
z1FgOO&mcfndtj2-zNa9APp_dfH~5IPaomiVvp!CGF@}IIrS|F>BT_oYD1LPH_m?j<
zSxZWd3KjhBBZEfo{+zs(Hr=B0wD6xT5C)mj(kTR8h;;m
zHxD%>y%m{Y;+L5Z&y#udU+R0shs8&Grq0~3@#7h@eEuSM9_Z9a>78v_JUb$3U$HP~p3S`c#YsuO$Zm_U)>z{Bi}{P$b&;lf|C01Lk-erOd$|K0H@nQZZ@T=K(yJ4fSlUJ;
zH@a#I5X(WjVP-lPwTo!>$jkK~t3y-wEn6I{c+f!~=?&1QaoBz^5Vm8V~(?xO>Qpgk_z_2zn)
z%u1x0%t{=S;@`(x@rm?C-T-eaJ>FmU=
zAAs2PPwVIp(*eA~2pfhf(4MFCwV30tvEdM?jcb>Gk;Na#?q3w(hmWAQA$#j!OJSfe
zKt}3!=baYpf6}Y1)z|53sfW+>cIP$gp^E1WslQkm@O-Lt{>b@r5nr8$b@pntksRbc
z@A+;`SFw#@DgCAVUAwjFkUVRtST>sr4a~I~+u`Fav5kukjnT5%7k)5%=8C8l^LjNw
z?J~rn2Mx95O^toEXr$-xn$mk396S8&%H`8yNUzbf@7w+sb@K9Rtt0Nd*Q8*2*QV3d
z#ow;@WXp%Xogdoz;F-5sy%_f?;RG?g^wo$jwbD4^Wr0a?&s@ZcHN2xg96e_A80x%3
zPFgfh&ks748vU;A%IdAF4n&@)Z~!Slo-$9CYVc*7HT7aVi;FaQvbOGs)bnY!r?pB)
z_RU$d6)zfF;>DhYfu~p*y=ZQ7%hv6CHtOMf>s?#d`F&=78M%D^^hHGN64>Rf(s23y
z&GQy*S>V%Z>8)*@I<`ovALafgJ-fz&AjI)B*%i-(sIyR>N#
z`TSSS(^|!7#hXlDQE^CPE@n=%a%S^np2>v~P>-6=A8Xi9QFj1>pHdg8=cZQ9wB$K_Qg
zFEA{}p0N1aFs}DyR`&w?+@|PU^}hV-vW>1-2g+6J%NjYoCe@vvwZ}EZ9$n^@H8VV|
zNpt1bYx+i?@Lv0Uqb5!9Av|l6+Ocn9)JtuaojAF8>H6iflM>h{D8jDEeM@aXAnYw(
zhM82pvt>0_ZHQcHV@_()$@3yx;Dwb>1ZnR2F0DBWDULt8j`gDY3xOuQh}&R4q{p$r
zuDu9bI<24jCDdj0RcOK0`#)Tw8`RtTCT`hrMum6||wZpM`R-g9md-iEUw%Yh}EYLh?wN=@FU_$2n=*
z)HK{f`_e0U^(7azCG5bY@2coXUvP+dv%Rs_4O_Z-+F8d2{VdDY?Vd>*>nJH8XWKdD+mow|F8i^%t^ed0pDx?EAH)Ha*tfL0@Wf_ZY(Xa~K3(VP0WH
zztfsVBht!-khDzwNe#bD@YQ+V6^PMOZqAYSXoW2r7Tn&ObM(Tu=b!=0Za&{Ho~QT=
zYi}d;zvblOu8femp4P?X5BG2JYR#hq_`Bdiye3=)g$%vw-*r`H`C19~e8T)@sa6&}
z^bN;42D-5Qzn;>ZkyjcVU<)>1jpDj1Mjz+UmD=k{em~bX0!6NFXSbM$@0WixXI0E)
zM^dNugA=0?lBO*4ZS{ZX{c2M0M)3(RJkyZPG95pC@zqx`A}>hb!u-J+H|bk7dSLx8MKvY%Tq$*5zp8?1*F7bU-eK{t=+0R)T?hqz+5)H2i_NM^hve&(ze}(}(mPNv$?P1*{a6aZOuq|`BeXJT4CzL{t9Twn7u*2
zEN~DE)TURZHnMBd*sEH%1;1-j#5_$-ucOJHLpN)&ku~1nTC*AB`}%0T4ZrvG_xaaa
z8-4os*(-mEdFh0$U4xo^pN-NR;Rc`-f=XO#Bg<>VwG`;ps80LV&-s$=y5hImglHvy
zpr;KpY|rhvn$xD$!V%tFog1%Rv%Yicu$Ch-&|Jf*&h1S)``d+}cn`~evnBgRV+1iM
zzxQTKy~OR7;Gz<6#tyS4sSC8ITP_zHC
ze(C0!zUODGOgtEM`pQNzknL#QH}Jc+X8U!g(-v&_>8AxRl`Pr5K~1L6js3x^HFf!d
z<;z#et5*HA?z;_}S6L9efJ#!CvLK^$q73ftj-TLE_P8$l^z+tA&u#nd(2u~jIu$F`
zsZ-H+Onbrmb=O%!%TA*+%mW+&z
z*1Kv?Sn{r5OC#K3Mg%(f8+oSag_np(;*UjtGrf!#PD<<$A-F34IdkiYKu$KM>8ZD_
z6Ycq0_~Y#7ivVmJ$(!Zg`
z76=z2-yq0j3jtZ*p?z!V>o4tPmjL}+oqd1AZUJpvx1LNdYYEKE?w=?e>=m;ICgPmd
zjdf*8^kmTpTFqn8fmHKfS^kJ-(2r1$lZB(KQ#8cge-&s6cp2W&6xtIsh2AQo1tyCI
z&`vL*4Cd=b7VS?hpwG{2gUO-^v_CaLplC)TOcrgR(FW|aMk`FtXoR2cX#;(#Xobl%
zg1?Yxg+PytYW&))Ma@1rIuM^(k(14O)P}z;YBsZ4r=;N6?!3dZCM`N|imh&^w(Xi^
zY`)fmT{k@CvM~pyj7N#qJ?`ek}r*a?ug`NUbcZ2}K&0NTaj8Yrjus&H40;
zEjSdIqK9<1z0|2=+oq9qw{*P*-Ru*MB?jBg>SHsvRTq!+=h-{8?@myfdp~vqeMj8A
z*V{A#y@~g;f3+p>OiM&?A87B(_^pVu=-hk16I)vpv=7kVx4F|dXa(`?73k$p;NwMZ
zzT&MJ?P)FcY?&hVY{^SUBPHHi^mlGLy|MdjdOkaDGjZ3Gz*R$nPnh&|$6~d@5rKSu
zH>tE2I~up1L>c$s;D^msr|bpo*POU~Le|%^v#&)jKQ7vxzFlm94&SbLTI4jtc2n<#M6Y;VRQ>3_8C5?;W7D(Nvt~@&G&zP&ie_ST$6h@<
zvysYH_GEJ-2M>}1HT`AjsYivF2k8x+p}Tfg4pi`qW<><$aRl-Usx8cm%17TWa?@wI
z)h&B^JA@YX?YG(qqtcC8Y{fpJ6$&y)i~
z#x3m^8%lUuufsUUGo^*_a@s1Urx0s1*=`~Pl<`^-Omsxp96xTp;oGbiY3*p*H>FX`
zlREAUwR~e!w3|D$oPDSF)kgM>@O!r73;Q6;{I#}6lQ!!e%ZW|ao3CIC*Vv<4yLvO?
z^fopvsWx6r^%O^Wt^EtVkLvC9a$0PQ1hjz_^iIF`s;9L@f6}{_#t7H8{PcQShht*G
zxGX^=6{~mL-0PUu5mowjt+(Tdp165>IebeJclmNXmYHVz*cqGkw&)l7?7(}_IRBz?
zTAOipz1=uR%#;d7u=DL;RZNIWa1IXEN(6_b*YYm*$O9EgmU2l2q{7nkQj#=Mnk9W9
zt&=uMyXA1Xm|Q~cA`g?lkk`xSlwc*l606iyo>!VH?UY{1>&jAPg|bQ6q5P^Gve+#i
zON^z8rMjiQ5UcAMLl%U0G_$M%M8
zw(UFHCfioq&$eG}hiunu0lU@iw1?PB+dJBm?IY}??Bng9+NayUvTwF;vmdY@v7fS^
zv)^(&;b`jU*!iRmZI$TaYU#H0a5oT0uR71_ZqxG%DzWpz%R-
zgT4<+3%c#JIEy;#IU70?oNb+5oV}a_ogbL4@b%VN7-K`!`PpTKxYpOq33U&ks
z2Nw>m65KerNAMfL?*z{bUK#v-@TTCOgZBj=3O*5hF8FHj%@8>xS4f$ViXqiQo(*{+
zq+v*>kij7zhD-`s9I`BAb;zcW?IC+Zj)tU#TyV?oZ0;iNr`+}2&D54+Rcm)$o*rBFwxJJcJRJ2Y=-xzJ}r8;7#(qT`B)eUPHmJpU0)+H<{tbf>$u#sW!hK&iE
z5H>YzcG&!|Z^D*`tqR)~b|CCbSbBI+c#iNw;pM|Chu02&CA>v=m+*e!Bg03BPYjfhi?i0DSS`(f$&q|*CONycSKY~k%;mUH6rRqyb_TR(JrEMM4yOZ
z5${HP5-~GkUc`!s?;|!v?1*0^r&Y3RaV_=m-}M)wbSxureLr0FVy&Xv-+c0EJ(T?d%jXvcjwR6B
zC9?{>0HUkjzTX0Q9csqr)B3=4Cc$BWD4YHcgv4F>B%8MN3H0zbPW*Ds<}Dxp@S*S6
zSOghyiO*GuQN38z=O}+okbtAta0`^Sf1?eW3(nZznGLI(Ga*y|
zVzf3dp~3hrm>?$?IiTl{V)chZXNP(9wTWL0nG>~X@wcbm^)<66zuvEZm#DreW5)FN
zsn>cJSZ&sM-Z}*;uWNp#kMC@XZTEnw1G_|a>e#1t^Mr|OMOjs^p5NcXhV88u8yF}i
z2=V%8mM~~tggsJuyPM&e;WqZ^KMunxo8xFvc5^2){bhSym#QV7VytQA^EamuogJlI
zY8-SzHTONt+6b_3w8{=`jPO-mICh`u((aNIUtV8c^#fnomc!iqLd#5kAF3)eq!@dZ~25vm_|S~EvKq~-p;jQqs!cx
z1Mfd&J(q?>BYxQ~HW^79j7xE^Sv#U&?6fwMDmE8ukBRo}v$pBG=)i%=cNb{;+}`Uu
zmK@mBegV_R`r~aaKN|4il*ma_$9);&**AMeYVv5`b&RyO40~lo!KiXITRii^H(j<5
z@@*PoTl(6kL;FW2rwnSJ?rq1gEV5A6Dn4may&v)k#$*;&r5YI7aJ91#h|SZyrR=`bB<
z9qwB}Ue(;o)Tvcgd%0p%?TCFFtcjYh9?RNv@8(+9NFZxpSUSwz-C~hab#?hs`Rni5
z_3Qq$yEM#;vp0LIw>o!g8dd*c8!3B_*E)fpzJ2pnelx*W;joshw{yTKi=8iEcN^KU
z;wfnJ2jRBmp8u-mdU{_3Et}~a0I|l5RIR_NRn4ohQ--S3v&1H=#VGObef
zy)4dXz-n?=zq@{Uptet#o2adslCLH
zosaNO_>!s>AF73yN~f-HLz+aS8_^(lElP5D8yWNcLAQ=c39t4j?CxB2lZ
zy3qnEroEc)AiStvzH+SYqeFJr@Dz{SuxEFDOlNJ3t%7}_*zqscK7$vYuczzdY+8By
zLhWG|rAn6c*&vLDnxxA)b1U4G3rg|*+G<5VoA6U(wY`&6VL!(AhyeY$F$J9}6^
z(ICqyQUZ#}g`Os!g{%Iw%-OPT-8(4|r8?61E1
zw0n!@1ABC)CKI)E8yl)8>glg=SJb{iXkK2hFVb!=)^Cfu^{8KpfqRY2{hnt_D*O=s;W0H4mF9Xa7Ylof_Kc7AO
zlX>%J5ANQh$AAueS#z%1bkV2v&99GeB=+sux^>@~i-hcHM1H+518kAiFXGwGEf(v}
zpWHvCPrm^jI`*45b;{(==7D*G^g(JAOYfXjb+z`q)(hMDiZT5GsYtv`%R09R{bo4jK^S1Sb@nIR$xqEum>}(z=X{h
zv#(8>!{`5v0e6j~iI{+?r)^pp6tBBp
z>93hRgHvjr_{ogUG3+99+V*jm58AhhbhA+3wwk`?*|p2u>}>JfJ;|13Zg%$8gi#Ye
zi_|{M-rdzbY3!%{Vzee&8{2yW-x<(9(k+a_gjRo?sqmlPpx5a3QFZDk)@)Au<@Wk9
zhz=aIqu5;F4X`Q$>&Y(S3FWD)1;Ca-*lEqMZ`BPnQg+Y2t=dhSe)BX_-N82PmKGj(
z*WNx*Q!75$wr<+W`8%R^ZEsej*b9xTs1|Kl`a6rEBjIcws##a=S^w^lAjWGN?vmcw
zK4Q#>F~i3MsR!C_i;wYiX{-10b{_cox8Htx+=WXje_7qjTpN485Ml3E@~*v~xtza`
zxt#xN)$}Oz1j-m(EYUy9!Uqn7=>KLH&CAxj`-Bi9EklT5G9T)V)feaaJPVNfl8-qE
zT1UQ;uG$67zCiW1+tBXd<;cqm_AJ_jsJZQn
z?Vc_aS!iM9MK8vngezvHJ(xrl?pvN|L&CVgJsogisF;)hsUc?Ts!q(=Cq(`io{&cd(;}YDcNv+*Del
zyvAcIJYu{9&2@7dHLzP-qa{YOQ|NB}XT
zJAYpHa`oyBUw-DZ>*p$7W)O!^P1Y(k3#%&2h=xAE^}kFx=)^l_nd1MhMHbte;}0+
z=oOcGt;7}0az*X<)~me}qSQ*e>Ocxq3|tCnih85?8Goa`sutz&tkRyC_}gyH5$Rd2
z*VGPqu`*ZhZobDRbrD%x)UCuj-%u@=OWoi)@TC}&RquvqHMRSkgsAwos_+ICp?VK}
zv0#f&`!c`&iKFYlHl2wRS%YGkygTB}wBDNP}uhf!V)p6IO)z_oRmMZ-Ce2y;~6#-4B?$K&O
zx9RiW5)|Fh<&!zR7f1cFNEHi(o`abS!~42AZR#>)QPb?6(&rqFy0>omn$LaL)_^7#
z&5j)gYxJJ^Su43eV(Ja;`5T^z{_JYeV|t}1y*_@gz51v6RJPy>RZVlTIhc!eGR?(U
zA18KL1~tS6wdxcsu30lTMc-^wh3gR+hh7$KE#4_z>yb*U`aIJzqq)`v)}{H|pQqQK
zt7@fxxsV2t)WXxFXqpoG6ROx^vyIqw_s-Uls&|h*v6_SRkJgMk_pMNyR^P5|_M-(X
zADHHF3wzh0*t@Ff#n7N(OQB+hXqpkLEKKYIjx)@b(1t@pBD8K=xLV$v;Hq8iol-uU
zu1=Xc;x)z~_(kMI9lW03?x5xBa}(6PT02dtx;NtUD_X=AT(8wr6Z&BGH+9kSJ7#|{
zVbLW^J=d;Ylf}JMRA`&E@R^!*E5`HK^+u@L&9}9jyR@7`GNu__zl1+xr_<8reqj_I
zO%H5R`W|629IIQIGdJ*Rj=&<^tgBkteOg%+kv27Gb+|tkdlrmOdVR>i-cjB9&!7{D
zEk$R~cx~XM7+5rTS}o4s*?T}JF3X#WOEtTeu_#?P7o{UZj==R!ixLCxJ>r9&>E_G^m;~@7O%p__|$`_yF=(OT7l@xpyDDej`>*$y=A1{Qq^+f
zn6j6OSQFNrxusQC)$x-?;6q7m+oK1YPpf;UYci3=7R45ITo-s2Scr(3I3NC4dlxOp
zPudsi!o9mzb$Wh+(DQqqNL5D*_1^fOhBf@AN>ur}D&`pM+f+g;m-oD0&ZqkD4_d|!
zHJfeI`r6eyn_&D@Q+qN(a8-J3-NOoVz1;Pp=y|=r?^v99-w6^xN%2J)7WmBf6q!-l&sD!s@g+Y1$qQ0^!Ydf(sQt@9=U-cZ+6t?F=@~f!&
zLhwI5NNndW6859J5MJTcY1*mJDQl|6Ro6rH
zu%k6pufM!rI8ff9F1a|vm&cwslqaD~=0cfE?9tz<=5TLY^3@uUoHcwZy`!~?2^O3`
zt1ngiS7+$P)%4uewA?n`Y1i`TLu`7Xt6Cuk=G8-NZu>F~m4nOP63qj$9(E&Y)~g3Q
zk3d89JfIHlghZ1Dm)CQp#;Dcn>pA=S2A3eHwQ!1mCC!*=Mk?BWhp^AArN;Eqt()vO
z4ygV`URAdvHF55v{9KdvsYu=>DHsrkBE?EpDV%2>sR-v{QfCePqJOa7AcSMog0
z3*;r7m&waGuaL>D@`3UJ=gG=sNw$PoLO5r$6p<{JI7=Moik7;Pg`7%~OIojlNy;-d
z>o$_yueM3(zZGkCf7&{tJ1e)J7uAs74K*&m)b^OhXoJWu~S`A(?j;$(?y;l|nP`Y*JX}ot^f~
zh~bbTGVg+D-;8HYDO={9OY%zDyS7j2E{!zLqs{XK^E|^m&oj@<%=23Fyv;oCG0#WL
z^O&#~tDDf9dc^<)8Lv&Q}%cM9JVt@QU^
zZM0NwsiYJqJuNw*9#%6jF7rIRJel4^O+GSM>2GS{-!vd?nC8f?vHt!S-j
zZDCEa4zj*xebf4$b*y!obvF98W!BZ!_0}!cgVrmyAX|*BqOF1LRoiQ}v9=|)mA19E
z&9?2fG=z#)d!qeS`v7~Y{X=wXR~)$=B^@msiH_Nh6^?U`+d;8GPX{##Y8lins7uhG
zpx1)l44Q^EZDG)|ptV6i1>JN8I}16BJD+yeb+&XSJ4ZU-an5rtcJ6Ura=Bb#l0|7_
z+)zF+t|$|Y0GOZ3d8A~Q!z8QRRIZ5Z&wI?nb2qW?o~oxA@ngqpCj~H@hvGzN(G-sirORBD6i4)0C!?lbT#lP3|YxTdB>f)aEf|zVQn+xlvg_Z7wu^QWo*-2fn<8
zdfUdcAB}4kEBzrjvsBujYWP}a4PP*8c!JqB7tLB0F>X+jqm=&|u`W=Wqmr8!Nk
zUn$M+#5ziTZcw80^wU&hA0>(gCqm?gMjCkGp?!X#eSV>2$0*r#N_Lpkeo|_JQMHUC
zq$*lkVAc)gL)vPrv6<4HrF7@PmA%yLPTJ=XU)V)y1>PK@q$g;rU5rqtsA0vp1pcMb
z{(kU5r*19sX`TmzFQMkYGsMeI335>)Ky$@N@%VL@AJ@@$au+%ZeFFOi{}sW)P)P3{45*d$M99uB`Ydxj5f4&*!M^
zU_wMt!#;YCc$0BHP25m0&Sz}9=jkzM*>&Rws8}j=eI#(*5aqInTmj%MPtxRuv_VtJ
zZa%eUKHYiGQ_|p{g)lNLa6`#W3)sObjo4cl!GTX!$y*JbkOAdsvhl)p7!YeU~p1nb+;JKK|+TgCxueIq5QY1E-^*BOWS;sJV79WzJYLsQn#
zHtUTmeE&MRItsR~=ZhPR%f#DDZf=nqdNx?U#W+W9E|QxJDp!eq_q607R{g$P-fiI_r!(j5~B%uY6ca&Wt=oi7?8p#;n~bGr4e@zanDnxUy1uW<@sHy
z1S=9xt3OB0wdUK2(3&35mH|8;2+n7fgZkxd#u#%+3MO5R{9GZ;qojGAG=C?}3#56T
zv`&-OMbbJ!TECFS8M8cxNaF};oB%u0=#}fB=o`R?jmB@Jal@?ZE#z-2^go^!e2y=+
zhB_rurXJ8YQ7;O0dW(P8dEbV7wk7>`(98C`e*?^Y%eYM0d)z(Kt~bn<{+YX-wBI%A
z8s39{gyf8XpPa0umIP+5=j#IJP7vy66MrsHt^gr73QEBFn)y_e_LhlNIIl^8f-3x>|O}ZD%{QOG18<}6)XXar;W*&YdtVmJx5aubQdc};f
zn;4gQeuU?{&2OjiMBrpbisy)b&aBZNp!%Y}o#aUxPc9Q~Ke@PSrgJ(okMuLjDC*z_
z7!XfyZAw3C4*kjCDFR9&?ba5W*dBT#`g#y!?>bt4Jzw8Io;Ok+(L1(sf0TMT2{jOO
zK=2WlN#zo4lr;h#Hs7T2<~nc0h?>P=2%56dq{CuNQKTq&9VtrkQPOkD8k3n&1tTgc
zmlC7|Ns-uP`J^Z;vtp%EN(H5YR9dN`RF}#q@nB7yQd6lVRZyN&>PVH8CQyebmDc3H
zy3#}GA=OewDkG)Z%A3lY(sRma