{ "cells": [ { "cell_type": "markdown", "id": "49bcb5b0-f19d-4b96-a5f1-e0ae30f66d8f", "metadata": {}, "source": [ "### A simple `A <-> B` reaction between 2 species with initial uniform concentrations across 3 bins,\n", "with 1st-order kinetics in both directions, taken to equilibrium\n", "\n", "Diffusion NOT taken into account\n", "\n", "See also the experiment `reactions_single_compartment/react_1`" ] }, { "cell_type": "markdown", "id": "a8445ac5-5fca-43f7-9858-cdd35c84772f", "metadata": {}, "source": [ "### TAGS : \"reactions 1D\", \"quick-start\"" ] }, { "cell_type": "code", "execution_count": 1, "id": "e3fee905-f0fe-4537-8352-2f07345d25cd", "metadata": {}, "outputs": [], "source": [ "LAST_REVISED = \"Jan. 13, 2025\"\n", "LIFE123_VERSION = \"1.0.0rc2\" # Library version this experiment is based on" ] }, { "cell_type": "code", "execution_count": 2, "id": "7245be7a-c9db-45f5-b033-d6c521237a9c", "metadata": {}, "outputs": [], "source": [ "#import set_path # Using MyBinder? Uncomment this before running the next cell!" ] }, { "cell_type": "code", "execution_count": 3, "id": "3cddd49a", "metadata": { "tags": [] }, "outputs": [], "source": [ "#import sys\n", "#sys.path.append(\"C:/some_path/my_env_or_install\") # CHANGE to the folder containing your venv or libraries installation!\n", "# NOTE: If any of the imports below can't find a module, uncomment the lines above, or try: import set_path \n", "\n", "from life123 import ChemData, BioSim1D" ] }, { "cell_type": "code", "execution_count": null, "id": "e049ec5d-7502-48f2-a79d-c647fe732e00", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "id": "8af818bc-9bac-4ec7-9672-b76223876595", "metadata": {}, "source": [ "# Initialize the System" ] }, { "cell_type": "code", "execution_count": 4, "id": "de47e1df-67cf-4200-a580-fbcf50ebd1c0", "metadata": { "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "SYSTEM STATE at Time t = 0:\n", "3 bins and 2 chemical species:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SpeciesDiff rateBin 0Bin 1Bin 2
0ANone10.010.010.0
1BNone50.050.050.0
\n", "
" ], "text/plain": [ " Species Diff rate Bin 0 Bin 1 Bin 2\n", "0 A None 10.0 10.0 10.0\n", "1 B None 50.0 50.0 50.0" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Initialize the system\n", "chem_data = ChemData(names=[\"A\", \"B\"]) # Diffusion NOT taken into account\n", "bio = BioSim1D(n_bins=3, chem_data=chem_data) # We'll specify the reactions later\n", "\n", "bio.set_uniform_concentration(species_name=\"A\", conc=10.) # Same across all bins\n", "bio.set_uniform_concentration(species_name=\"B\", conc=50.) # Same across all bins\n", "\n", "bio.describe_state()" ] }, { "cell_type": "code", "execution_count": null, "id": "7b653cda-902a-40d4-b1d4-33005a8b4274", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "id": "e648c35b-0ec5-47b3-b759-67d14f7865b2", "metadata": {}, "source": [ "## Enable History" ] }, { "cell_type": "code", "execution_count": 5, "id": "db8a132d-febf-45b4-8ca2-4e3cb5eb2dc4", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "History enabled for bins None and chemicals None (None means 'all')\n" ] } ], "source": [ "# Let's enable history - by default for all chemicals and all bins\n", "bio.enable_history(take_snapshot=True)" ] }, { "cell_type": "code", "execution_count": 6, "id": "64b4eaff-5e55-4a15-a148-d271ccdd9d11", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SYSTEM TIMEAB
00.010.050.0
\n", "
" ], "text/plain": [ " SYSTEM TIME A B\n", "0 0.0 10.0 50.0" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "bio.get_bin_history(bin_address=0)" ] }, { "cell_type": "code", "execution_count": 7, "id": "e07f8c5e-270d-4703-b30d-1e019f611b46", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SYSTEM TIMEAB
00.010.050.0
\n", "
" ], "text/plain": [ " SYSTEM TIME A B\n", "0 0.0 10.0 50.0" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "bio.get_bin_history(bin_address=1)" ] }, { "cell_type": "code", "execution_count": 8, "id": "d7c855b6-6a52-4700-b8e9-04a55c4b3b0c", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SYSTEM TIMEAB
00.010.050.0
\n", "
" ], "text/plain": [ " SYSTEM TIME A B\n", "0 0.0 10.0 50.0" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "bio.get_bin_history(bin_address=2)" ] }, { "cell_type": "code", "execution_count": null, "id": "1c2ccfa6-4fc4-4340-8d38-09ce125acdaa", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 9, "id": "25155c63-e53c-41e1-be1e-41577158a4ca", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Number of reactions: 1 (at temp. 25 C)\n", "0: A <-> B (kF = 3 / kR = 2 / delta_G = -1,005.1 / K = 1.5) | 1st order in all reactants & products\n", "Set of chemicals involved in the above reactions: {'B', 'A'}\n" ] } ], "source": [ "# Specify the reaction\n", "reactions = bio.get_reactions()\n", "\n", "# Reaction A <-> B , with 1st-order kinetics in both directions\n", "reactions.add_reaction(reactants=\"A\", products=\"B\", forward_rate=3., reverse_rate=2.)\n", "\n", "reactions.describe_reactions()" ] }, { "cell_type": "code", "execution_count": null, "id": "3210ae02-198a-4d4a-9232-184265d9147c", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "id": "0b46b395-3f68-4dbd-b0c5-d67a0e623726", "metadata": { "tags": [] }, "source": [ "### First Reaction Step" ] }, { "cell_type": "code", "execution_count": 10, "id": "bcf652b8-e0dc-438e-bdbe-02216c1d52a0", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "System Time is now: 0.1\n", "SYSTEM STATE at Time t = 0.1:\n", "3 bins and 2 chemical species:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SpeciesDiff rateBin 0Bin 1Bin 2
0ANone17.017.017.0
1BNone43.043.043.0
\n", "
" ], "text/plain": [ " Species Diff rate Bin 0 Bin 1 Bin 2\n", "0 A None 17.0 17.0 17.0\n", "1 B None 43.0 43.0 43.0" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# First step of reaction\n", "bio.react(time_step=0.1, n_steps=1)\n", "bio.describe_state()" ] }, { "cell_type": "markdown", "id": "7dc56592-179d-4e4c-b75a-8eb81dcafe71", "metadata": {}, "source": [ "NOTE: the concentration of the chemical species `A` is increasing, while that of `B` is decreasing.\n", "All bins have identical concentrations; so, there's no diffusion (and we're not attempting to compute it; didn't specify diffusion rates) " ] }, { "cell_type": "code", "execution_count": 11, "id": "499da263-4575-466f-841a-05490c75bb3f", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SYSTEM TIMEAB
00.010.050.0
10.117.043.0
\n", "
" ], "text/plain": [ " SYSTEM TIME A B\n", "0 0.0 10.0 50.0\n", "1 0.1 17.0 43.0" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "bio.get_bin_history(bin_address=0)" ] }, { "cell_type": "code", "execution_count": null, "id": "9a753a5f-25be-46e5-9fa0-ae1e05777aa3", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "id": "82a62165-425d-4e8d-abac-01d579dfd1ae", "metadata": {}, "source": [ "### Several more reaction steps" ] }, { "cell_type": "code", "execution_count": 12, "id": "cf6a7337-8e2e-4c02-9bb3-85052f37144f", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "System Time is now: 1.1\n", "SYSTEM STATE at Time t = 1.1:\n", "3 bins and 2 chemical species:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SpeciesDiff rateBin 0Bin 1Bin 2
0ANone23.99316423.99316423.993164
1BNone36.00683636.00683636.006836
\n", "
" ], "text/plain": [ " Species Diff rate Bin 0 Bin 1 Bin 2\n", "0 A None 23.993164 23.993164 23.993164\n", "1 B None 36.006836 36.006836 36.006836" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Several more steps\n", "bio.react(time_step=0.1, n_steps=10)\n", "\n", "bio.describe_state()" ] }, { "cell_type": "code", "execution_count": 13, "id": "363053b6-e9cf-45d5-a77c-b0f4560cb51e", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SYSTEM TIMEAB
00.010.00000050.000000
10.117.00000043.000000
20.220.50000039.500000
30.322.25000037.750000
40.423.12500036.875000
50.523.56250036.437500
60.623.78125036.218750
70.723.89062536.109375
80.823.94531236.054688
90.923.97265636.027344
101.023.98632836.013672
111.123.99316436.006836
\n", "
" ], "text/plain": [ " SYSTEM TIME A B\n", "0 0.0 10.000000 50.000000\n", "1 0.1 17.000000 43.000000\n", "2 0.2 20.500000 39.500000\n", "3 0.3 22.250000 37.750000\n", "4 0.4 23.125000 36.875000\n", "5 0.5 23.562500 36.437500\n", "6 0.6 23.781250 36.218750\n", "7 0.7 23.890625 36.109375\n", "8 0.8 23.945312 36.054688\n", "9 0.9 23.972656 36.027344\n", "10 1.0 23.986328 36.013672\n", "11 1.1 23.993164 36.006836" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ "bio.get_bin_history(bin_address=0)" ] }, { "cell_type": "code", "execution_count": 14, "id": "0e059c0b-e956-4b5f-bcbe-c60025df3717", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SYSTEM TIMEAB
00.010.00000050.000000
10.117.00000043.000000
20.220.50000039.500000
30.322.25000037.750000
40.423.12500036.875000
50.523.56250036.437500
60.623.78125036.218750
70.723.89062536.109375
80.823.94531236.054688
90.923.97265636.027344
101.023.98632836.013672
111.123.99316436.006836
\n", "
" ], "text/plain": [ " SYSTEM TIME A B\n", "0 0.0 10.000000 50.000000\n", "1 0.1 17.000000 43.000000\n", "2 0.2 20.500000 39.500000\n", "3 0.3 22.250000 37.750000\n", "4 0.4 23.125000 36.875000\n", "5 0.5 23.562500 36.437500\n", "6 0.6 23.781250 36.218750\n", "7 0.7 23.890625 36.109375\n", "8 0.8 23.945312 36.054688\n", "9 0.9 23.972656 36.027344\n", "10 1.0 23.986328 36.013672\n", "11 1.1 23.993164 36.006836" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ "bio.get_bin_history(bin_address=2)" ] }, { "cell_type": "code", "execution_count": null, "id": "b7c73b92-2d91-4644-8aeb-66f833a2fc25", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "id": "962acf15-3b50-40e4-9daa-3dcca7d3291a", "metadata": {}, "source": [ "### Equilibrium" ] }, { "cell_type": "markdown", "id": "809b4afa-fb2f-4ac3-92c9-083fc487c81b", "metadata": {}, "source": [ "NOTE: Consistent with the 3/2 ratio of forward/reverse rates (and the 1st order reactions),\n", " the systems settles in the following equilibrium:\n", "\n", "[A] = 23.99316406\n", " \n", "[B] = 36.00683594\n" ] }, { "cell_type": "code", "execution_count": 15, "id": "490fdc0f-fa2a-48d8-ae04-7c80d298842a", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "{'A': 23.9931640625, 'B': 36.0068359375}" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "bio.bin_snapshot(bin_address=0)" ] }, { "cell_type": "code", "execution_count": 16, "id": "f22258b2-5181-4ff9-b379-f6a12ad5c8fb", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0: A <-> B\n", "Final concentrations: [A] = 23.99 ; [B] = 36.01\n", "1. Ratio of reactant/product concentrations, adjusted for reaction orders: 1.50071\n", " Formula used: [B] / [A]\n", "2. Ratio of forward/reverse reaction rates: 1.5\n", "Discrepancy between the two values: 0.04749 %\n", "Reaction IS in equilibrium (within 1% tolerance)\n", "\n" ] }, { "data": { "text/plain": [ "True" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Verify that the reaction has reached equilibrium\n", "bio.reaction_dynamics.is_in_equilibrium(conc=bio.bin_snapshot(bin_address=0))" ] }, { "cell_type": "code", "execution_count": null, "id": "77741c9d-a432-4911-979e-9b36856a3ce5", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "id": "cbf6c9c7-8cec-400f-9e70-49ff1a9f485c", "metadata": { "tags": [] }, "source": [ "## Plots of changes of concentration with time" ] }, { "cell_type": "code", "execution_count": 17, "id": "71c3c2e2-0411-4c3e-a607-31c8196d1c7f", "metadata": {}, "outputs": [ { "data": { "text/html": [ " \n", " " ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/vnd.plotly.v1+json": { "config": { "plotlyServerURL": "https://plot.ly" }, "data": [ { "hovertemplate": "Chemical=A
SYSTEM TIME=%{x}
Concentration=%{y}", "legendgroup": "A", "line": { "color": "darkturquoise", "dash": "solid", "shape": "linear" }, "marker": { "symbol": "circle" }, "mode": "lines", "name": "A", "orientation": "v", "showlegend": true, "type": "scatter", "x": [ 0, 0.1, 0.2, 0.30000000000000004, 0.4, 0.5, 0.6, 0.7, 0.7999999999999999, 0.8999999999999999, 0.9999999999999999, 1.0999999999999999 ], "xaxis": "x", "y": [ 10, 17, 20.5, 22.25, 23.125, 23.5625, 23.78125, 23.890625, 23.9453125, 23.97265625, 23.986328125, 23.9931640625 ], "yaxis": "y" }, { "hovertemplate": "Chemical=B
SYSTEM TIME=%{x}
Concentration=%{y}", "legendgroup": "B", "line": { "color": "green", "dash": "solid", "shape": "linear" }, "marker": { "symbol": "circle" }, "mode": "lines", "name": "B", "orientation": "v", "showlegend": true, "type": "scatter", "x": [ 0, 0.1, 0.2, 0.30000000000000004, 0.4, 0.5, 0.6, 0.7, 0.7999999999999999, 0.8999999999999999, 0.9999999999999999, 1.0999999999999999 ], "xaxis": "x", "y": [ 50, 43, 39.5, 37.75, 36.875, 36.4375, 36.21875, 36.109375, 36.0546875, 36.02734375, 36.013671875, 36.0068359375 ], "yaxis": "y" } ], "layout": { "autosize": true, "legend": { "title": { "text": "Chemical" }, "tracegroupgap": 0 }, "template": { "data": { "bar": [ { "error_x": { "color": "#2a3f5f" }, "error_y": { "color": "#2a3f5f" }, "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "bar" } ], "barpolar": [ { "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "barpolar" } ], "carpet": [ { "aaxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "baxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "type": "carpet" } ], "choropleth": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "choropleth" } ], "contour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "contour" } ], "contourcarpet": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "contourcarpet" } ], "heatmap": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmap" } ], "heatmapgl": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmapgl" } ], "histogram": [ { "marker": { "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "histogram" } ], "histogram2d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2d" } ], "histogram2dcontour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2dcontour" } ], "mesh3d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "mesh3d" } ], "parcoords": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "parcoords" } ], "pie": [ { "automargin": true, "type": "pie" } ], "scatter": [ { "fillpattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 }, "type": "scatter" } ], "scatter3d": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatter3d" } ], "scattercarpet": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattercarpet" } ], "scattergeo": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergeo" } ], "scattergl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergl" } ], "scattermapbox": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattermapbox" } ], "scatterpolar": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolar" } ], "scatterpolargl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolargl" } ], "scatterternary": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterternary" } ], "surface": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "surface" } ], "table": [ { "cells": { "fill": { "color": "#EBF0F8" }, "line": { "color": "white" } }, "header": { "fill": { "color": "#C8D4E3" }, "line": { "color": "white" } }, "type": "table" } ] }, "layout": { "annotationdefaults": { "arrowcolor": "#2a3f5f", "arrowhead": 0, "arrowwidth": 1 }, "autotypenumbers": "strict", "coloraxis": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "colorscale": { "diverging": [ [ 0, "#8e0152" ], [ 0.1, "#c51b7d" ], [ 0.2, "#de77ae" ], [ 0.3, "#f1b6da" ], [ 0.4, "#fde0ef" ], [ 0.5, "#f7f7f7" ], [ 0.6, "#e6f5d0" ], [ 0.7, "#b8e186" ], [ 0.8, "#7fbc41" ], [ 0.9, "#4d9221" ], [ 1, "#276419" ] ], "sequential": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "sequentialminus": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ] }, "colorway": [ "#636efa", "#EF553B", "#00cc96", "#ab63fa", "#FFA15A", "#19d3f3", "#FF6692", "#B6E880", "#FF97FF", "#FECB52" ], "font": { "color": "#2a3f5f" }, "geo": { "bgcolor": "white", "lakecolor": "white", "landcolor": "#E5ECF6", "showlakes": true, "showland": true, "subunitcolor": "white" }, "hoverlabel": { "align": "left" }, "hovermode": "closest", "mapbox": { "style": "light" }, "paper_bgcolor": "white", "plot_bgcolor": "#E5ECF6", "polar": { "angularaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "radialaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "scene": { "xaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "yaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "zaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" } }, "shapedefaults": { "line": { "color": "#2a3f5f" } }, "ternary": { "aaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "baxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "caxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "title": { "x": 0.05 }, "xaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 }, "yaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 } } }, "title": { "text": "Reaction A <-> B . Concentrations at bin 0" }, "xaxis": { "anchor": "y", "autorange": true, "domain": [ 0, 1 ], "range": [ 0, 1.0999999999999999 ], "title": { "text": "SYSTEM TIME" }, "type": "linear" }, "yaxis": { "anchor": "x", "autorange": true, "domain": [ 0, 1 ], "range": [ 7.777777777777778, 52.22222222222222 ], "title": { "text": "Concentration" }, "type": "linear" } } }, "image/png": "iVBORw0KGgoAAAANSUhEUgAAA00AAAFoCAYAAACL/RCAAAAgAElEQVR4Xu29C7Qd1Xmgue9bQg9AIIEQDk7HEGKQbYwBd8xyswyrPRZxM6OMwYMJDd1qMD3MmhGMjdeoTSZ20wvwGHqm6YVxmBFNZBKgQ8d2LGcm4NcibvMKjgUBTOIABgTiYYQe6Er36s7+65xdt07dOufUY1fVrl3fwdf33HN2/fvf31/nqr67q3aNzOmH4gEBCEAAAhCAAAQgAAEIQAACiQRGkCb2DAhAAAIQgAAEIAABCEAAAv0JIE3sHRCAAAQgAAEIQAACEIAABAYQQJrYPSAAAQhAAAIQgAAEIAABCCBN7AMQgAAEIAABCEAAAhCAAATyEWCmKR83toIABCAAAQhAAAIQgAAEWkIAaWpJoRkmBCAAAQhAAAIQgAAEIJCPANKUjxtbQQACEIAABCAAAQhAAAItIYA0taTQDBMCEIAABCAAAQhAAAIQyEcAacrHja0gAAEIQAACEIAABCAAgZYQQJpaUmiGCQEIQAACEIAABCAAAQjkI4A05ePGVhCAAAQgAAEIQAACEIBASwggTS0pNMOEAAQgAAEIQAACEIAABPIRQJrycWMrCEAAAhCAAAQgAAEIQKAlBJCmlhSaYUIAAhCAAAQgAAEIQAAC+QggTfm4sRUEIAABCEAAAhCAAAQg0BICSFNLCs0wIQABCEAAAhCAAAQgAIF8BJCmfNzYCgIQgAAEIAABCEAAAhBoCQGkqSWFZpgQgAAEIAABCEAAAhCAQD4CSFM+bmwFAQhAAAIQgAAEIAABCLSEANLUkkIzTAhAAAIQgAAEIAABCEAgHwGkKR83toIABCAAAQhAAAIQgAAEWkIAaWpJoRkmBCAAAQhAAAIQgAAEIJCPANKUjxtbQQACEIAABCAAAQhAAAItIYA0taTQDBMCEIAABCAAAQhAAAIQyEcAacrHja0gAAEIQAACEIAABCAAgZYQQJpaUmiGCQEIQAACEIAABCAAAQjkI4A05ePGVhCAAAQgAAEIQAACEIBASwggTS0pNMOEAAQgAAEIQAACEIAABPIRQJrycWMrCEAAAhCAAAQgAAEIQKAlBJCmlhSaYUIAAhCAAAQgAAEIQAAC+QggTfm4sRUEIAABCEAAAhCAAAQg0BICSFNLCs0wIQABCEAAAhCAAAQgAIF8BJCmfNzYCgIQgAAEIAABCEAAAhBoCQGkqSWFZpgQgAAEIAABCEAAAhCAQD4CSFM+bmwFAQhAAAIQgAAEIAABCLSEANLUkkIzTAhAAAIQgAAEIAABCEAgHwGkKR83toIABCDgLIH/49a71ea7v6u++40b1K+tOcrZPEkMAhCAAAQg0BQCjZKmF156VX3iM9cksv03/8vvqf/hvz27du4XfPYP1Cs73lQ/vO//7MnF5H7yib+u7v7a79eeZ1ICf/XwNnXZ578avOUKT9ug/sn6/1m9/ubOxLBVHGCag9l4AlX0bZtlmnj9Pg9ptk3TpsmftzTjkzbmc5nlM5lWmsr+vST1eeLpfwiH+k//yWnq5j/4H9MOnXYQgAAEIAABZwg0Upri//CaAwQX/kFu8kHcxt//j+r/++Ejwc5pW+5MjS694BPqf73igtQfAHPAaKu2Ik3yiEutObgrS16iwh/vw3DPyiY1xBobIk3F4TdVmpI+ayeddYmy9VkuTpYIEIAABCAAgfQEvJAmGa456H3yB3ekH30JLcs+SCwh5TCkOaCRF0SebAjEH//ZA+rf/vs/CvrI8pfy6DijM2BFD7j6SVOeA9MstejXr4kh/f/Xx/42k1Bm6b+utmV/HsqOXxe3pP0/y+cn7UxTWeMzn/uv33i1+sjpa8NuzOs2freUlTtxIQABCEAAAkkEvJempFP64v+QR2dYopD6/cOedIqXtP3cl7/WcyqKiWXi9Psra5pTWMy2p59yYighEr+oRJgcowc58pqcpldk5sOGLMV3WBvy1E9eTL5ZDkzT/krJG1tqHn3Ec4uK3sOPPx3OEg4S1H77bvS6l3ib+H4QPSD/vf/p3/Wc7hj9zMT3636fh//+3I+Gp4UeueLQYBYwKU/znomTNn78dLBhXCV+ls9bdF83uaXZj5K2i9cuus9H94Vhn/tBNYr/YSn+eynPfhX/LJjfqfG+zO/jIr9b0n7uaAcBCEAAAhCwScAbaUo6GDYHJdF/oJMOYOUf+DVHH9nzV37zj370IND8gx8/YJG2IjNyTdWgv3wnSZO8Fj8Vrt9rcUmyOTsSz3vYzEi/nbAMWRokT1kPvvqNK+9403wYk/alQduZukbHNui1+IF20ixDmn036fqWpO2i12XFJSl+Pd+wz4PknnQqqNTj331hQ88sRVKNhsWPflbTcjXSlObzZlhE/xCTtt7yWfmzv3iw5xrHpN9PeT7ng2ok1xhFZaafNKXZr/rtx1nqkuYzRBsIQAACEIBA3QS8kKZ+Byn9Znak/V9v+/mC61qixTAHi9G/GKc5FSjLwUK/U2iSDpL6jUUOJI9etaLQ4hJJf/1NOhgctLNGZ/TS/JXdxo6fZ+GKQQtBDPvrfd6cs5462m8fiu8v/Q6mk+pZZN+Nn1LVb79NOiUry+dhGN+kU7uyxE/L1UhT0v4Q/7z1iym5fuS0k3OtXCd9fHDtCeGCCUWkKT5bnrRvDJppii6uk2WWaNAfIbiuadiezvsQgAAEIOAigUZKUxxk0sHNoAMNc9AX/Wtrv9Ngon/tT/OPvY2DuKSDtjKlKcvMRL+dOK80xU+Vkvjx0yf79ZlXmiRefCEIwzx+CpiND21WaepX6/g+nUWa0uy7w/o1dalKmvqtNBjdP7J83oaNLyr7aT9v5g82RRZOSfoMROPZlCbZn+PihzTZ+JQTAwIQgAAEfCfQSGmKSpI5II3PbvS7XiBaUPNXWDP7EI0R/6tq2r+yZjmIGzRLFD9FL+1BXJ4ddtDsi8TLsrhGU0/Pk3HmvfZoGPMsK/MN2s/ip5ellaY0++6g5fzN+Mzno2xp6rfSYJGZrCxck/5oYRgkfWaTrolMc9qo2S7+R594H02Upiy/B4d9fngfAhCAAAQg4AKBxkuT+cup3Hsn+hdoc6AxbNYi7YHnoAOpaCGzHCxkaVuWNA06ICsiEWXIU5kLQUgNk655sfEhzcox7YyI7X1X+k1zsF+2NJUVPy3XrNIU30fMHyGG/e5J+5m2LU3xfsuYaWIhCBu/OYgBAQhAAAIuEfBCmqJ/mR50gXMS+H5L42a5LkTa/tUjTwQLQQy6Xip+cOLCNU3DLlpPWpQiyw5sQ55syFJ0pkCeJ52el1VusnAwp2D1m7WTfeieb/0gWIwk7bU3WaSpX8zovpvmuicZcxapyfJ5MDz7HXAnfVazxE/LNYs0Sf/x1fnSSM6gma/4TFPS9ZXD9r0iv1uy7Ff98mDJ8WEV4n0IQAACEGgaAS+kKTpLEL0mJWn1PGkrByGyVLIcOCetDhaVsKQVzKKn08QPaPodLPQ7EJOD6fh1NP1Wz0tzYXrWHXCYFA2TqrT9uX5zW3OgmPbaFDObkPZ+M1lubmtySbvyW/zU1KQD8qSY8X03afW8qGz+0X/434KFDbJIU9bPg/SXJK9R+Y7O4GSJn5ZrFmky+0FUhtN+ZpIWSzDx4vth1gVf+i3kkuZ3iw1pEobc3Dbtb0faQQACEIBAEwh4I03Rgy15bg5m+12rEZWheBuRGDlA/MRnrkk8XSl+4Xb8lKb4NQ4279MU/6t21oOp6E6ZZmalrFPW6vpwDLp+K82paVGJkFNC00qT2a7fPcGSZqDi+9mg+zSlXeVs2L5rDnZlbNFH9CA+izRJjKyfh/hnWX4WgTT3KBt2n7Vhn7dhXLNIk2kb35/TXgcY3x8ld/mDTnxFzPjvqGErPZoaye+yaC2Ttivj9DzDI8096Or6XUC/EIAABCAAgSwEGiVNWQZGWwhAAAIQgAAEIAABCEAAAjYIIE02KBIDAhCAAAQgAAEIQAACEPCWANLkbWkZGAQgAAEIQAACEIAABCBggwDSZIMiMSAAAQhAAAIQgAAEIAABbwkgTd6WloFBAAIQgAAEIAABCEAAAjYIIE02KBIDAhCAAAQgAAEIQAACEPCWANLkbWkZGAQgAAEIQAACEIAABCBggwDSZIMiMSAAAQhAAAIQgAAEIAABbwkgTd6WloFBAAIQgAAEIAABCEAAAjYIIE02KBIDAhCAAAQgAAEIQAACEPCWANLkbWkZGAQgAAEIQAACEIAABCBggwDSZIMiMSAAAQhAAAIQgAAEIAABbwkgTd6WloFBAAIQgAAEIAABCEAAAjYIIE02KBIDAhCAAAQgAAEIQAACEPCWANLkbWkZGAQgAAEIQAACEIAABCBggwDSZIMiMSAAAQhAAAIQgAAEIAABbwkgTd6WloFBAAIQgAAEIAABCEAAAjYIIE02KBIDAhCAAAQgAAEIQAACEPCWANLkbWkZGAQgAAEIQAACEIAABCBggwDSZIMiMSAAAQhAAAIQgAAEIAABbwkgTd6WloFBAAIQgAAEIAABCEAAAjYIIE02KBIDAhCAAAQgAAEIQAACEPCWANLkbWkZGAQgAAEIQAACEIAABCBggwDSZIMiMSAAAQhAAAIQgAAEIAABbwkgTd6WloFBAAIQgAAEIAABCEAAAjYIIE02KBIDAhCAAAQgAAEIQAACEPCWANLkbWkZGAQgAAEIQAACEIAABCBggwDSZIMiMSAAAQhAAAIQgAAEIAABbwkgTd6WloFBAAIQgAAEIAABCEAAAjYIIE02KBIDAhCAAAQgAAEIQAACEPCWANLkbWkZGAQgAAEIQAACEIAABCBggwDSZIMiMSAAAQhAAAIQgAAEIAABbwkgTd6WloFBAAIQgAAEIAABCEAAAjYIIE02KBIDAhCAAAQgAAEIQAACEPCWANLkbWkZGAQgAAEIQAACEIAABCBggwDSZIMiMSAAAQhAAAIQgAAEIAABbwkgTd6WloFBAAIQgAAEIAABCEAAAjYIIE02KBIDAhCAAAQgAAEIQAACEPCWANLkbWkZGAQgAAEIQAACEIAABCBggwDSZIMiMSAAAQhAAAIQgAAEIAABbwkgTd6WloFBAAIQgAAEIAABCEAAAjYIIE02KBIDAhCAAAQgAAEIQAACEPCWANLkbWkZGAQgAAEIQAACEIAABCBggwDSZIMiMSAAAQhAAAIQgAAEIAABbwkgTd6WloFBAAIQgAAEIAABCEAAAjYIIE02KBIDAhCAAAQgAAEIQAACEPCWANLkbWkZGAQgAAEIQAACEIAABCBggwDSZIMiMSAAAQhAAAIQgAAEIAABbwkgTd6WloFBAAIQgAAEIAABCEAAAjYIIE02KBIDAhCAAAQgAAEIQAACEPCWANJkobQvv/GOhSiEcJHA5PioWr5kQr2+c9rF9MjJEoHVKxarV371jpqbsxSQMM4RWLFsUu2dnlX79s86lxsJ2SGwZNG4Gh8bUTv3HLATkCjOEZD6rlg2pXa8tc+53MpM6JgjFpcZntgpCSBNKUENaoY0WYDoaAikydHCWE4LabIM1MFwSJODRbGcEtJkGaiD4ZAmB4vSopSQJgvFRposQHQ0BNLkaGEsp4U0WQbqYDikycGiWE4JabIM1MFwSJODRWlRSkiThWIjTRYgOhoCaXK0MJbTQposA3UwHNLkYFEsp4Q0WQbqYDikycGitCil1kvT7Xd9R9389XsXlPzJH9wRvnb557+qHnx4W/DzmaevVbfdeHVPe6TJ308M0uRvbaMjQ5r8rzPS5H+NkSb/a4w0+V9jl0eINGlpeuSnTy8QIVO0L918p3pp+2vh+yJQa1avVNduvDisK9Lk8i5eLDekqRi/pmyNNDWlUvnzRJrys2vKlkhTUyqVP0+kKT87tixOAGkaIk3rLrpGXfeFDeqUk48PaD/+xLNq0/W3q61bbkCaiu9/zkdAmpwvkZUEkSYrGJ0OgjQ5XR4rySFNVjA6HQRpcro83ieHNMVOzzvu2KNCIdq+4011zvlXqfvvuUmtXrUi2BmSXmOmyd/PCdLkb22jI0Oa/K8z0uR/jZEm/2uMNJVb46Szqcrt0V70pONze9E7kVovTXGgssPIQ65bSiNNR9x4hPo3Z/6B2vCBy23XhngOEBgbHVFTk2Nq774ZB7IhhbIILF08oXbv0/d24T5NZSGuPe7iqTF1YGZOzcwerD0XEiiHwIS+r96o/p09zb24ygHsQFSp72L9b/Kelv2bvOyQCSv05ZKTu7/5vZ5Y0cmCOqRJzuC66Mrr1JZbNoVndeUZLNKUh1rBbUzxZCGINNI08gcjQY/HLD1GXfbBf60uP+VfqyWTSwtmweauEECaXKlEuXkgTeXydSE60uRCFcrNAWkql68L0ZGm/FU46axLVFSQTCQRqWOOOkJtuPBcVYc05R9R75ZIky2SGeJEpUk2G3ZN013b7lI3/fj/Uo+98lDQy9LJZerStZ9Vl6y9XB2tRYpHswlwel6z65c2e07PS0uque04Pa+5tUubOafnpSXV3HacnpevdiJGP3nsyZ7r8ZMiGWmS98yMVD/Ris5YRS9jkePm9es+qu7b+iP1/IuvBt185YtXqNVHrQhmlMzDbJMkO/EZMdl+3dlnqKSZMrPaNdKUb9/ItJUUN7qog/z84VNPClfHS7t63v3PfVfdse1r6vvP/2XY/z9fe5m69H2fVccffmKmnGjsDgGkyZ1alJkJ0lQmXTdiI01u1KHMLJCmMum6ERtpylcHmWXaeNmngtmkQQ9zix0jKdJ22HHx1gceUp/78q3KyIu0F1kyUmTej58GKLGTLoWJC57I0B9+48+D43J571995nfCdQaGXVKTj1b/rVp/TVP0HkyC6YLzPtaznLi8luU+TQ9t/yu1+WdfU99+9k9D6utP+HQgTx88+nTb9SNeyQSQppIBOxIeaXKkECWmgTSVCNeR0EiTI4UoMQ2kKTtcMwOT5pqhpNPz5H6m5tY8/WZzRJSuvHR9MBtkZpqMoCVtIzFlJkomLaLvy+hkAbY0uUpbEbJbNt+3II5ZvC07rcFbtF6abABNWj3vqdef0PJ0q/rG324Ouzjn3Z/Qp+5doc467hwb3RKjAgJIUwWQHegCaXKgCCWngDSVDNiB8EiTA0UoOQWkKTtgm9JkLmFJysLMTvWTpqgI9ZOdV3a8EZzCZ2atkvoxM1nR9/qtQ5CdFtJkm9mCeIOWHH9x1wvBzJN8Tc/uC7b98DFnBjNPv/Oe9aXnRgfFCCBNxfg1ZWukqSmVyp8n0pSfXVO2RJqaUqn8eSJN+dhlOT1vzeqVPWdcRWea4tf995MauaYpPtNkQ5pkHNEzwqKnBnJNU759o/Kt0tynaee+t9Tmbbfq655uU6/t3RHk+N4j1mp5ukJdeNIlledMh+kIIE3pODW9FdLU9AoOzx9pGs6o6S2QpqZXcHj+SNNwRkkthi0EIWLUb/W8pNPzBp0+V2SmSXLvd3pekrAhTfn2h1q3SiNNJsGDcwfDmad/2Pl3wcvvWn5cMPMkq+5Njk3VOhY67yWANLVjj0Ca/K8z0uR/jZEm/2uMNOWvcdKS40ZEzCIRw65pkt7NCnbRU+hErE593wnBfZaKSJNciyQ5/PLlHeEibWYhCFkAIi5UMiZ5cHpe/v2i8i2zSFM0uXuf3hII1N/s+Ovg5UOnDgtmnkSgjly8svJx0OFCAkhTO/YKpMn/OiNN/tcYafK/xkhTsRonLdkdnTVKI01RcYpmE109L+/peWYBh/gibSZHkbObv35v2K1cR2VW7uP0vGL7RmVb55Umk+Bf/OLbgTw9+OL3g5fGRsfCmad3H/oblY2DjpCmtu4DSJP/lUea/K8x0uR/jZEm/2vs8ghZPc9CdYpKk0nhxy/9UP0/f3Or+u4vvhVmdf5v/V4gUO9beYqFTAmRlQAzTVmJNbM90tTMumXJGmnKQquZbZGmZtYtS9ZIUxZatLVNAGmyQNSWNJlUfrbjcb1oxNfUPU/9UZjdJ/7RPwvk6SPHnmUhY0KkJYA0pSXV7HZIU7PrlyZ7pCkNpWa3QZqaXb802SNNaSjRpiwCSJMFsralyaT03M6/DxeNmJ2bDV4WaRJ5EoniUT4BpKl8xi70gDS5UIVyc0CayuXrQnSkyYUqlJsD0lQuX6IPJoA0WdhDypImk9ob+17vytOt6q19vwpefv+qD6pL1l6u5PQ9HuURQJrKY+tSZKTJpWqUkwvSVA5Xl6IiTS5Vo5xckKZyuBI1HQGkKR2nga3KlibT+f7ZaXXHz25Tm5/4mnph53PBy79+2HuCpcpl9ml0ZNTCaAgRJYA0tWN/QJr8rzPS5H+NkSb/a4w0+V9jl0eINFmoTlXSFE31rifv0LNPt6q/fWNb8LIsUR7c60kvWS5Ll/OwQwBpssPR9ShIk+sVKp4f0lScoesRkCbXK1Q8P6SpOEMi5CeANOVnF25ZhzSZzv/87+4LTt37ycsPBi9Njk515emzwU1zeRQjgDQV49eUrZGmplQqf55IU352TdkSaWpKpfLniTTlZ8eWxQkgTcUZqjqlyaT/wxfuD+TpL5/bGo7oMyddqq97+qx675FrLYyynSGQpnbUHWnyv85Ik/81Rpr8rzHS5H+NXR4h0mShOi5IkxnGX7/ysLpj223qT5/543Bknzz+d4PZpzNWf8TCaNsVAmlqR72RJv/rjDT5X2Okyf8aI03+19jlESJNFqrjkjSZ4Tz75tPBvZ7+07avhyM869fOCa55Oufdn7Aw6naEQJraUWekyf86I03+1xhp8r/GSJOfNV530TXBwLZuucHpASJNFsrjojSZYb26Z3t4r6fdB3YFL3/w6NODmaf1J3zawuj9DoE0+V1fMzqkyf86I03+1xhp8r/GSJN/NX78iWfV1+78lvrlyzvUdV/YoE45+XhnB4k0WSiNy9Jkhrd3Zk8oT9t3vxS8fMKK3wrk6eKT/5UFCn6GQJr8rGt8VEiT/3VGmvyvMdLkf42RJv9q/KWb71Qfet9vqkd/9kwwuGs3XuzsIJEmC6VpgjRFh3nntj9Ud+h7PT3zxlPBy0cvOSZccW/JxFILRPwJgTT5U8tBI0Ga/K8z0uR/jZEm/2uMNBWr8bd3vl0sQIGtP3no8sStTzrrEnX/PTepV3a8oTZdf7vTp+ghTQV2ALNp06TJ5P1ffn53MPv02CsPBS+JMHXu9fTZQKR46CXcx0fV8iUT6vWd0+DwmADS5HFxu0NDmvyvMdLkf42RpmI1/vRzL6i7f/VWsSA5tv6Td/+auuDwhfcQNafm3Xbj1UFUEagtt2xy9hQ9pClH8eObNFWazDgeeP4vghvlfv/5vwyHJqfsiTzJKXxtfiBN7ag+0uR/nZEm/2uMNPlfY6SpWI1FmvYePFgsSI6tP6OFKUmazKl5684+I4gqP8vD1VP0kKYcxfdNmsx4Hn75x8GKe9969j+HQ1z/m59Wl+p7PcniEW18IE3tqDrS5H+dkSb/a4w0+V9jpMmvGsvMUtLjyR/c4eRAkSYLZWn6TFMcwVOvPxHMPH3jbzeHb8ky5bJcuSxb3qYH0tSOaiNN/tcZafK/xkiT/zVGmvypsZyal3QNk8un6CFNFvY/36TJIHnx7ReCmSe57ml6dl/wstwg9+O/8TvqQ0d9WJ26ujOd6vMDafK5uvNjQ5r8rzPS5H+NkSb/a4w0+VPjyz//VbVm9coFp+K5fIoe0mRh//NVmgyat6d3BjNPIk+vvbMjJDY1tkidevQZ6kOrP6xO018fOvrDavnUoRaIuhMCaXKnFmVmgjSVSdeN2EiTG3UoMwukqUy6bsRGmtyoQ1uzQJosVN53aTKI5tScuvepLerR7T9Rj77yE/XMm50ly6OPk1e+P5CoU4/SMqUl6tcOfbcFwvWFQJrqY19lz0hTlbTr6Qtpqod7lb0iTVXSrqcvpKke7vTaIYA0WdgT2iJNcVRv7Hu9I1BdiZLvB+d6V2V51/LjAnkyM1JrV37AAvHqQiBN1bGusyekqU761fSNNFXDuc5ekKY66VfTN9JUDWd6SSaANFnYM9oqTXF0MwdnghmoqEj9at+bPc2WTS3vnMqnr4kKZqT06X2LxhdbqEI5IZCmcri6FhVpcq0i9vNBmuwzdS0i0uRaRezngzTZZ0rE9ASQpvSs+rZEmvpDfPqNJwOReqQ7I/Xczr9f0FiWMxeBkhkpEapVhxxtoSp2QiBNdji6HgVpcr1CxfNDmoozdD0C0uR6hYrnhzQVZ0iE/AQql6Z1F12jnn/x1cSMXV2XfRhepGkYofn3t+9+qed0vp/ueGzBxr9x+AmBQAVfeoW+Om+wizSlr22TWyJNTa5eutyRpnScmtwKaWpy9dLljjSl40SrcghUKk39lhcsZ2jVRUWa8rPeO7MncjrfQ8HzPQd29wQ8cvHKzkxUdzZKvo/o/6p4IE1VUK6/D6Sp/hqUnQHSVDbh+uMjTfXXoOwMkKayCRN/EIFKpcnlG1YV2U2QpiL0Fm77Nzv+umc26uXdL/Y0Gh+bUKd1r4kyInXYosPtJtGNhjSVgtW5oEiTcyWxnhDSZB2pcwGRJudKYj0hpMk6UgJmIIA0RWBtfeAh9bkv36q23LJJnXLy8eE7MkP24MPbgp/PPH2tuu3Gq3sQI00Z9rgcTZ9/+xfhNVFyfdRTrz+xIMpvHXly93S+znVRxy3/Rzl6WrgJ0mQFo/NBkCbnS1Q4QaSpMELnAyBNzpeocIJIU2GEBChAoFJpEvk47+NnqnVnn1Eg5XI2FWG6ZfN9wfVWUWmSOxO/tP21UJSSTjFEmsqpSb+ov3rnTfXYq51T+YJFJvTXzOyBnubHLD2253S+96/6YK4kkaZc2Bq3EdLUuJJlThhpyoyscRsgTY0rWeaEkabMyNjAIoFKpenxJ55Vm66/XW3dcoPFIRQPZYRJ8oqfQigLV1z3hQ3hzFPSGJCm4jUoEkHuDSUCFYqUfnoYKX8AACAASURBVP76O6/1hFwysbQrUfqmu93row4ZXzK0W6RpKCIvGiBNXpRx4CCQJv9rjDT5X2Okya8aJy0O5/KicJVKkwjJoEcdoKLCJLlFpWn7jjfVOedfpe6/5ya1etWKIPWk15Am9z7Ez7zxVOeeUfrrsVceUn//q58vSPIDq07tmY1avXTNgjZIk3u1LSMjpKkMqm7FRJrcqkcZ2SBNZVB1KybS5FY9imYj0nTlpevDM9DiZ3cVjW97+0qlyXbyRePFhSmvNBXNg+3LJ/DyrpfVj3/54/DroZceWtDpe1a8R/32u347/Fq7am35idEDBCAAAQhAAAIQaCGBuDQlHZe7hKXV0iRGe/c3v5dYj6988Qp1ytrjmWlyaW+1mMu+mXfmV+iTGantD6ld+9/u6eHwRSvU6av/sfrIcf9YHTaxWq1Z9i61Zqn+0t9HR0YtZkOougkw01R3Bcrvn5mm8hnX3QMzTXVXoPz+mWkqxvjbz3y7WIACW3/yNz+5YOu4NMm6Aad94ES14cJzC/RU3qaVS5NcE3TRldf1jCi+Wl15wx0emWuahjPytcW2134auWfUT9SLu17oO9Rjlh0byNMxSzrfzZcsQCFiVdYS6L6yr3tcSFPdFSi/f6SpfMZ194A01V2B8vtHmoox/vR//rS6+8m7iwXJsfWf/O6fqAtOviBRmmQBtuhj42WfQpoEiFnSO+kaIZnZcWFVvbg0sXpejk+HJ5u8sPM59dPX9PVQO59Sf/f6c+ql3b9UL+36pYrfNyppuEsnlmmR6opVMDt1rDrGyJX+WeRqfHTcE1LNHwbS1PwaDhsB0jSMUPPfR5qaX8NhI0CahhEa/L5I094De4sFybH1Z9Z+pq80Ra9pktByHO6KE8SHWulMU3waziTj0jmMSTfg5T5NOT4hnmyStBCErNYn8iQSJQIVPO9+vbxH//z2Lxec6peEQxaeCGarurNT8dmqwxd3Fh/hUT4BpKl8xnX3gDTVXYHy+0eaymdcdw9IU90VsNt/khe4fIpepdKUJCSC35yyV8fqeTbKz+p5Nii6GSPv6nlvT+8MxcrMThnRMj/Pzc0NHLQsk94Rqc4pf4FcdWerglkrfWrg+NiEm+AalhXS1LCC5UgXacoBrWGbIE0NK1iOdJGmHNAc3iQuTWaFamaadNGaMNOUZ99CmvJQa8Y2eaVp2OhEmEKZ6p72F/9ZxGvY4+ilx8zLVHeRilCu9M8rFh8xLATvawJIk/+7AdLkf42RJv9rjDT5VeOk+zS5KkxCvtKZpiZc05Rnd0Sa8lBrxjZlSVOa0ctqfuHslD7lr+dUwO6pgXKq4KCH3MA3nK0ypwJGF67Qs1UTY5Np0vG6DdLkdXmDwSFN/tcYafK/xkiT/zV2eYSVSpOAcH31vDzFQpryUGvGNnVKUxpCfWerutdY7Zx+a2iYo5asXjBbdeQhK9WyiUPV8qnlatmk+b5cLZ1cNjReExsgTU2sWrackaZsvJrYGmlqYtWy5Yw0ZeNFa7sEKpcmu+m7EQ1pcqMOZWThujQNG3NntupFPUs1v/JfuHBFd7Zq9uDssDDh+3J/qmWTy7VMHaqlqvvd/JzwPRAu/brr4oU0pd4FGtsQaWps6VInjjSlRtXYhkhTY0vnReJIk4UyIk0WIDoaounSlAbry4FU6VX/ItdWvbnvDb0C4E719vTbwUqAcn2VfN9zYHeakAPb5BYvLWqBgOnZL1nS3eYDabJJ081YSJObdbGZFdJkk6absZAmN+vSlqwqkSaz5vrnvnzrQK6snteW3a4542yDNGWphsxKBRKlhSr8rsUq+DkQLCNa+rtuV6p4aXlarmeygpkv811mwAbNePURL6Qpy17QzLZIUzPrliVrpCkLrWa2RZqaWTdfsq5EmnyB1W8czDT5W2GkyX5tQ/GS2asD88JVtXiNqJFg1kqEa8Uhh6nxkSk1Oaq/xqbU1Hj3u/55anyRfm0yeH3R2PzzKf1c2k3p1+W94H1pO9ppK9tNBfEmuzH0axJvotOPzLjxqI4A0lQd67p6QprqIl9dv0hTdazpaSGBSqWp332aXLq5bZ6dBGnKQ60Z2yBNbtZJVg00pw2aWa/oaYTmdMJQxA5EZsC6pxzu3r+r1sHJqoWBVI3r7yJj+vmirlCF0iWyFoqbSNe8gHUkTrbvSFggaRIzELvOa33jTeg2I5NqbGxMjY+Oq1HV+e7zA2nyubqdsSFN/tcYafK/xi6P0Alp4ua2Lu8i7c4NafK3/lHxmly0T7381ttqemZaf+1T+w/uV/v1832z+vnstP7aH7w+fXD++X79fJ+01e9Jm2ndVraXbeW5bC/tg9eCGCZeJ9ac/s+1R1SgRKjGRrRMjYyrMS1UwXP5Hj6X9zuvy2vyXrC9aSevm+djul1XzMz28/G620djBTHn4xqpm48//57ENf0H2+g448F387yT0xHLFqvpA3PqgP4aGRnR6GWuUf+nn3a+m58730dHZSaw+1rC+9H287Ei7ZO2Ge3to6fPuc57klA8H9f2E1fzQZpcrYy9vJAmeyyJlJ2AE9J0+13fUfdt/ZHauuWG7CNwYAtmmhwoQkkpIE0lgXUsbB3XNB0QEYtKlUjYAZEuI27yvStgoXTNC1hH4jpyF0haIGWR9wOJ0++beIHYzQvd/rn9anZ2Vs0cnFEHVec7DzcJLJQoLXX69M4k2TMiJiOJi2Dws8hZRMxE+EQQ+8bS7QMpTJC5+W06fQWx+uXVRwjDGJG84uI4n9+8xi4QVd3/+JjWa93/zGxHjJMkNMwvMqYkUQ2FuptXEoMwr26ssE2McVL9FvDU9KLshjKP5RXuufpvMZ1azD+iPwdM+rwnf8cZtG38/Z628jegAf0G20b6jjXtec/su0ljGtV/eFi6aFzteify+yqedzyPDGNeMP74ttGfY/0O27aHT8Zt/8Vp57v5y6llWZUuTUn3ZUpivOWWTeqUk49vJH6kqZFlS5U00pQKU+Mb1SFNLkKLCpQI1eyclqm5GTWrhSp4Lt/D5/J+53V5Td4Ltjft5HXzfFa364qZ2X4+Xnf7aKwg5nxcI3VB/0E++pVufxI3fC7b6Pdngu/meSenuaCd3k6f2jk3J0d4eq5Pf+8+6z43r3W/S1vzX9A29n735/lYYevEttKZzHB2+o3FGun8HMSKvefivkJOEIBAdQTmft+9MxOqG707PZUuTdGh9rumyR0c+TJBmvJxa8JWSFMTqlQ8R6SpOEPXIzT5mqaFEtURrIP6v0QBEzsL/tdfzkxMI3yDZE53MkQedZSIiPbLK0kIw7F1pTHQ2JhUHtSy22nXjSBPpH1MLqcm9IzZ6Jzas28m5BKPJdv0jDWhv3mmur+ozCbIbjRWKM+RsSTVp0eyQ7adsZmxJraJ9h/LK/r5i5/626lN96EnmqI/97zXVfaeWAO2lXa92/ce2C/oN/gTxfwjGnpgzpG89ESTmpwYU+9M986MLzjdORI8bWyTWQ+feM7Dfh7Qb7AD96lDnGU85+/+3rdc/xXbivwqlSZfiSJNvlZW/3IeH1XLl0yo13dO+ztIRqaQJv93giZLk//VsTNCrmmyw9HlKFzT5HJ1/M8NabJQY6TJAkRHQyBNjhbGclpIk2WgDoZDmhwsiuWUkCbLQB0MhzQ5WJQWpVSpNG3f8aY65/yr+uLl5rYt2vMaMlSkqSGFKpgm0lQQYAM2R5oaUKSCKSJNBQE2YHOkqQFF8jjFSqVp3UXXqPXrPqpOfd8JatP1t4er5V3++a+q8z5+plp39hmNRM1MUyPLlipppCkVpsY3QpoaX8KhA0CahiJqfAOkqfElHDoApGkoIhqUSKBSaTILQRy96gj1L6+6IZQmWWEvKlEljreU0EhTKVidCIo0OVGG0pNAmkpHXHsHSFPtJSg9AaSpdMS1d4A01V6CVidQizTJ0uIiUOZ0PG5u2+p90OnBI01Ol8dackiTNZTOBkKanC2NtcSQJmsonQ2ENDlbmlYkVqk0yWl4p33gRLXhwnNV9Dk3t23FvtbIQSJNjSxb5qSRpszIGrcB0tS4kmVOGGnKjKxxGyBNjSuZVwlXKk1xcjLbZB7333OTWr1qRSPhcnpeI8uWKmmkKRWmxjdCmhpfwqEDQJqGImp8A6Sp8SUcOgCkaSgiGpRIoFZpKnFclYZGmirFXWlnSFOluGvrDGmqDX1lHSNNlaGurSOkqTb0lXWMNFWGmo4SCFQqTWYhCLmmyacH0uRTNXvHgjT5W9voyJAm/+uMNPlfY6TJ/xojTf7X2OURIk0WqoM0WYDoaAikydHCWE4LabIM1MFwSJODRbGcEtJkGaiD4ZAmB4vSopQqlaam34+p336BNPn7iUGa/K0tM03tqK0ZJdLkf72RJv9rjDT5X2OXR1ipNG3f8WbP/ZlcBpMlN6QpC61mtUWamlWvvNky05SXXHO2Q5qaU6u8mSJNeck1ZzukqTm18jHTSqUpulpeEkxz36amgUaamlax9PkiTelZNbkl0tTk6qXLHWlKx6nJrZCmJlcvXe5IUzpOtCqHQKXSVM4Q6o+KNNVfg7IyQJrKIutWXKTJrXqUkQ3SVAZVt2IiTW7Vo4xskKYyqBIzLYFKpanf6nlbH3hI3bL5PrV1yw1p83aqHdLkVDmsJoM0WcXpbDCkydnSWEsMabKG0tlASJOzpbGWGNJkDSWBchBwQpoef+JZddGV1ylOz8tRQTYplQDSVCpeZ4IjTc6UorREkKbS0DoTGGlyphSlJYI0lYaWwCkIOCFNt9/1HXXf1h8x05SiYDSplgDSVC3vunpDmuoiX12/SFN1rOvqCWmqi3x1/SJN1bGmp4UESpcmM4s0DP6WWzapOm56K8J289fvDdM78/S16rYbr+5JV5ZKf/DhbcFrSe9zet6w6jb3faSpubXLkjnSlIVWM9siTc2sW5askaYstJrZFmlqZt18ybp0aYqC6ndNU50wRYiikiQ/r1m9Ul278eIgrS/dfKd6aftrYZv4+9IGaaqzguX2jTSVy9eV6EiTK5UoLw+kqTy2rkRGmlypRHl5IE3lsSXycAKVStPwdOpvITNPj/z06VCS1l10jbruCxvCWTCZOdt0/e09pxIiTfXXrawMkKayyLoVF2lyqx5lZIM0lUHVrZhIk1v1KCMbpKkMqsRMSwBpipESSfrwqScFM01yM95zzr9K3X/PTWr1qhVBy6TXkKa0u1vz2iFNzatZnoyRpjzUmrUN0tSseuXJFmnKQ61Z2yBNzaqXb9lWLk0iJc+/+GoixzpXzzN5Ra9ZSitNu/Ye8G2/YDxdAmOjI2pqckzt3TcDE48JLF08oXbv05/jOY8H2fKhLZ4aUwdm5tTM7MGWk/B3+BPjo2pU/86e3j/r7yBbPjKp72L9b/Kelv2bvOyQiZZX3o3hVypNSdcDuYFhPovoNUxIk2vVqT4fpKl65nX0iDTVQb3aPpGmannX0RvSVAf1avtEmqrlTW+9BCqVJhcXgojvEPF7RnFNU7s/Mpye1476c3qe/3Xm9Dz/a8zpef7XmNPz/K+xyyNsvTSJFG3dckNYI5kNk4dZUY/V81zefcvPDWkqn7ELPSBNLlSh3ByQpnL5uhAdaXKhCuXmgDSVy5fogwlUKk0iJOd9/Ey17uwznKlL9B5MkhT3aXKmNE4kgjQ5UYbSk0CaSkdcewdIU+0lKD0BpKl0xLV3gDTVXoJWJ1CpNCUt1+0DfVbP86GKyWNAmvytbXRkSJP/dUaa/K8x0uR/jZEm/2vs8ggrlSa5pmnQo87V84oUCWkqQs/tbZEmt+tjKzukyRZJd+MgTe7WxlZmSJMtku7GQZrcrU0bMqtUmnwFijT5WlmlkCZ/a8tMUztqa0aJNPlfb6TJ/xojTf7X2OURIk0WqoM0WYDoaAikydHCWE6LmSbLQB0MhzQ5WBTLKSFNloE6GA5pcrAoLUqpcmmKLryw5ZZN6pSTj1dy2t4F531MXbvx4kaiR5oaWbZUSSNNqTA1vhHS1PgSDh0A0jQUUeMbIE2NL+HQASBNQxHRoEQClUpT9Oa20fsfbX3gIXXL5vt6lv4ucczWQyNN1pE6ExBpcqYUpSaCNJWK14ngSJMTZSg1CaSpVLxOBEeanChDa5OoVJpkRun+e25Sq1etUFFpit9QtmnVQJqaVrH0+SJN6Vk1uSXS1OTqpcsdaUrHqcmtkKYmVy9d7khTOk60KodApdIkovR/33TNAmlipqmc4hK1OAGkqTjDJkRAmppQpWI5Ik3F+DVha6SpCVUqliPSVIwfWxcjUKk03X7Xd9R9W38UnIZnZpqOXnWEOuf8q9TGyz6lNlx4brHR1LQ1M001ga+gW6SpAsgOdIE0OVCEklNAmkoG7EB4pMmBIpScAtJUMmDCDyRQqTRJJuZUvGhWX/niFWrd2Wc0tlRIU2NLNzRxpGkoIi8aIE1elHHgIJAm/2uMNPlfY6TJ/xq7PMLKpcllGHlzQ5ryknN/O6TJ/RrZyBBpskHR7RhIk9v1sZEd0mSDotsxkCa36+N7dpVK05duvlPd/c3vqSd/cEcPV5Yc9303a+74kKbm1i5L5khTFlrNbIs0NbNuWbJGmrLQamZbpKmZdfMl60qlSa5juvLS9QtOxWMhCF92J//GgTT5V9OkESFN/tcZafK/xkiT/zVGmvyvscsjrFSaZEbJ3NA2CoUlx13eRdqdG9LUjvojTf7XGWnyv8ZIk/81Rpr8r7HLI6xUmphpcnlXILckAkhTO/YLpMn/OiNN/tcYafK/xkiT/zV2eYSVSpOchve5L98a3uBWwGzf8Waw5HiTV9BjIQiXd/FiuSFNxfg1ZWukqSmVyp8n0pSfXVO2RJqaUqn8eSJN+dmxZXEClUqTpJu05HjSKXvFh1ZdBKSpOtZV94Q0VU28nv6Qpnq4V9kr0lQl7Xr6Qprq4V5lr0hTlbTpK06gcmnysQRIk49V7YwJafK3ttGRIU3+1xlp8r/GSJP/NUaa/K+xyyNEmixUB2myANHREEiTo4WxnBbSZBmog+GQJgeLYjklpMkyUAfDIU0OFqVFKVUuTbIYxPMvvpqIOH7/pqbUAWlqSqWy54k0ZWfWxC2QpiZWLVvOdUjTnE4x+JqbC77LQ74f1P8fvmfahG3NNt3v0naku52OY2LEtz94sBtTt+1t03n9YDeHOTXSyUX/XxBDt59/z/RpYsz1vifbdHOJjmvhWEb0dgeV7kqPfT5mtJ3pc+F4Yn1GuC0YRze+iTE+rtToyIjau392AV9pE3JPyknnu2AcOlbQZ5dtEvvOOOZrFLbRGyZxNftBJxeJ361bJKcg17Be83WK708dtvr/JM/uvmH6j346gre6+0XwvnnB7CtmWx1n/n15pn8eMRl394mkWN3tOtua7Ux7s31nrPOP+M/dbYfkOaL7mhgbVftn9P7VoseP3/ueFo3W3aFWKk2Xf/6ras3qlerajRe7SyRHZkhTDmgN2QRpakihCqbZRmmSw5dZfYQzow9kZoPnSs3on2f1QUvweuz5QX0AJW3ldTlcmdEHknJoKtvO6G3lebDNgrjShxywduIH2+s+grZBe4nb6TOIJTHNcxPPxIznJ6934wbbSFzJK9xO+uu0kYP9WX1gLHn0HLzr9sHBXvQrIinm9bBN9+C256C22960iQpBwV2TzSEAAQiouVPeBwUHCFQqTf3u0+QAh0IpIE2F8Dm9MdLkdHmsJVelNO3Xh+fT+kB9Wn/frw/up/WB/LQ+EO+8frD7mryu3zev6ZHKz9I22EZ+lrb6eydGJJ5+XeKZttJG2gb9RfoWaeFRLgH5o3nwJX+w7/6pf7T7LHzPtAnbRbfpPpfJBP2fzKLII2lbeS/op1vWaJvREenVbCfzG52HtO+XT6efTp894whyMTlG3uvJa05vNxqYaGfsvV/RmPHxBPlE+zS5dmPMb6vH0c0lmBDR70+Nj6nR0RF14IAofqxfiWlyj+U0zzU2Hv0ZkTiGbRL7MB9t5SFXwzY2DpOTxOlw78SXZ/P7iTBL4NrNfcE+ENm3OqgkZmS6pptLdJJH4nenhLr7gbTv/DUgeC98yGuxWJF3zT7UG7uTw4JYwjIWe3ienVCdfjp56EkmtWzxhHprt/xma8/jd9cc0Z7BOjxSpMlCcZAmCxAdDYE0OVqYgmnJzMMuLRJvH5xVu7SIjB8yrl7ZNd2RDyMYRjxENLpiIoIi4rHPCIvIj5GW7ntJItQRpI7IuPKQQ5AxfSAyrg+wxoLnSo3rn8d0nsHrseejegzSVl7Xxy1qXB+cjnW3Hdfbzj/Xr3e3DWLIl95W3pefg+11H/Pv6efyftBOnpttum0kN5NrPL/ue2FOJu9urE6fne2PWDqppvdr0dRfvULROYLuOciOHKDOH+h2DkbNwW2PTJiDdzm401/mPVdq3ZY8uKbJ/0pzTZP/NXZ5hJVKk5yed97Hz1Trzj7DZSaZc0OaMiNrzAZIk3ul6hWeuY74BAJ0MBCgXfrn4Lm8Fv481yNJe/R7dT0m9WH1lD6on9LSMam/y8+L5Ln+Lj9P6b/QT+kj70670e5rnW2C10b1+zp503ZSbys/m7bReNJ2UgtJJ6bpQ6lF+mcRjTY96rimqU18XRgr0uRCFcrNAWkqly/RBxOoVJrkHk2brr9dbd1yg1d1QZq8KmfPYJAmu7WVU8LmZaYrOIHYdKUnEB5po79HhUc/78jRnNprYbZGZkqWaaFYrr+WaYE4cmpCjeqLeoy4iFRMRsSlIxxaTPS5IaH0yGuh5HTfN5ITyM28CM3Hi5/+Ypcv0foTQJr83zuQJv9rjDT5X2OXR1ipNMk1TYMerJ7n8q7SztyQpvm6d4RHz+zMdU5p68zodAXHzOwE37vCE4hOVIgOWhEeOY1rmRaS5aNjapk+USz4bn4WCdI/L9c/B9+1xIgUBXJkJEl/XyLXW0QeVV7T1M5PUv2jRprqr0HZGSBNZROuPz7SVH8N2pxBpdLkK2hmmnytrP83txWpeXl2Rr00e0C9dGCm83zmgHpTi1FpwqOlJxCdHpkR8YmIzYgIT1d8hgiPjb0PabJB0e0YSJPb9bGRHdJkg6LbMZAmt+vje3ZIk4UKI00WIDoaoskzTTIz1BGijggFzw90vpvX5bS4tA+5mL4zWyPC0xUf83N0ZkdekxmeYMZHvs+3jc/wpO277HZIU9mE64+PNNVfg7IzQJrKJlx/fKSp/hq0OYPKpUmua7royut6mG+5ZZM65eTjG1sHpKmxpRuauMvS9IYWnlCGtBQFz/WpcUaMXtFyNOyxVMvNmolxdczYuFozNqG/j6k14xPqiNHx7iluRpL0KW1agHx9IE2+VnZ+XEiT/zVGmvyvMdLkf41dHmGl0rT1gYfU5758q7r/npvU6lUrAi7bd7ypzjn/KvWVL17R2FX1kCaXd/FiudUlTfv0imdGfkIxmtVCFJkxkjaDHrI2mghQIETjXTEKf9avazE6TG56wUMhTf7vBEiT/zVGmvyvMdLkf41dHmGl0rTuomvUlZeuXyBHIlO3bL6vsavqIU0u7+LFcitLmmQWKDxdTmaHzExR91Q6mUUa9jhczwqJ+IRCNNEVpMiM0bAYvN8hgDT5vycgTf7XGGnyv8ZIk/81dnmElUqTrJ6XdCqeOWWP1fNc3lXamVseaZLrhOQ6ovDaoe51RIEYdV+X640GPSb1aXPmVLnOTFF0xqhzCp2r1wc1cU9BmppYtWw5I03ZeDWxNdLUxKplyxlpysaL1nYJVCpNLs40fenmO9Xd3/xeSPXM09eq2268uoey3JT3wYe3Ba8lvc9Mk92d0qVocWmaDRZX0DNDerU5s9JcfMZIVp0b9lgls0RafNZoIQq+60UTwp/1jNEq/TOP6gggTdWxrqsnpKku8tX1izRVx7qunpCmusjTrxCoVJpcvKZJRC56s135ef26j6oNF54b7CEiVS9tfy0UKRGoNatXqms3XhzuQUiTfx8mmQl6ZPod9fiBafXzuQPqH/btD64x2p5icQWZAerMDHVmhJJmjOSmqDzcIYA0uVOLsjJBmsoi605cpMmdWpSVCdJUFlnipiFQqTRJQq6vnnf7Xd9Rj/z06VCSRKKu+8KGcHU/yX/T9bf3iBbSlGZXc7vNG/q6okf2vaMend6nHt3f+d7vBLpQhOTaIb3Awvypc7K4wpiSa414NIsA0tSseuXJFmnKQ61Z2yBNzapXnmyRpjzU2MYWgcqlyVbiZcWRmaTTPnBiMNNkVvZLWu0v+hrSVFY1yov7D/r6okf1TFIgSvv3qWcO7F/Q2drJKXXG4sVqYnxM/dPxxd1V6CYUc0Tl1aWuyEhTXeSr6xdpqo51XT0hTXWRr65fpKk61vS0kADSFGEis0w3f/1eZRakSCtN7FjuE3hk7zvqx3v26K+96q927wlOtYs+Ful7EP32kkO6X0uC74cyY+R+YckQAhCAAAQgAAEIVECgEmky1zIl3Ytp0HsVjD/sYtD1Vsw0VVmJ4n3tnTvYOc1OzyQ9un86mE2S16KPlfo0utMWLVYfmlqkPjS5WJ2qvyc98qyeV3wERKiaADNNVROvvj9mmqpnXnWPzDRVTbz6/phpqp45Pc4TqESakhZPiBYhvthC1QWKzzBF++eapqqrkb0/ueeRuRbpES1LP9Vf8cdvjE8GgnSaSNLUYnX8xGSqjpCmVJga3whpanwJhw4AaRqKqPENkKbGl3DoAJCmoYhoUCKBSqSp3/2ZzLjqvE+TCJ084suMm9xYPa/EvS9n6Gf27w+uQwpmkrQg/WJm4fVIp0yKHIkkLQ6+ZInvPA+kKQ+15m2DNDWvZlkzRpqyEmtee6SpeTXLmjHSlJUY7W0SaLU0mWuWkoBGb8LLfZps7nLZY5nFGowkva5Xuos+luolvj+0SE6z64rSokPUIkvLNSBN2evVxC2QpiZWLVvOSFM2XZQl9QAAF+pJREFUXk1sjTQ1sWrZckaasvGitV0ClUhT/BS3+BCSlvG2O8xyo7F6nj2+cmNYmUUSUZL7JMnzA/qeSdGHLPMtp9gZSXq/lqWyHkhTWWTdios0uVWPMrJBmsqg6lZMpMmtepSRDdJUBlVipiVQiTTF730UT27YNU9pB1NXO6QpP/lfytLfwal2Ikp71ZMJS3+fOKmvR9KLNXSuR1qk3q2vT6rqgTRVRbrefpCmevlX0TvSVAXlevtAmurlX0XvSFMVlOmjH4FKpEk6l9kmeWzdckNPLvL68y++Gi7z3cRSIU3pq/bkgenuynZakqb3ql/OzPRsPKp/kmuQgpmkYGW7RWpFzuuR0mfVvyXSZIOi+zGQJvdrVDRDpKkoQfe3R5rcr1HRDJGmogTZvgiByqRJkjSr1EUTvuC8j6lrN15cZAy1b4s0JZdATqt7dL+cZteZSZLnO2d7l/4+TN8fKS5J4yPu3D4Waar941VJAkhTJZhr7QRpqhV/JZ0jTZVgrrUTpKlW/K3vvFJp8pU20tSprCzQ8Ki+FsmcbicLN/RejaTUcWMTnVXt9GINMov0W/rUO5cfSJPL1bGXG9Jkj6WrkZAmVytjLy+kyR5LVyMhTa5Wph15IU0W6txWafqFvv5o/iay+9QzCdcjrZ2Y6qxs11244djxCQvEqwuBNFXHus6ekKY66VfTN9JUDec6e0Ga6qRfTd9IUzWc6SWZANJkYc9oizTJynZ37347OM1OZEluKht9TMnS393FGszCDcv06XdNfiBNTa5e+tyRpvSsmtoSaWpq5dLnjTSlZ9XUlkhTUyvnR95Ik4U6+i5NL+oV7u7YvVNt3rVT7ZubvyZp1Zgs/d1Z0U4k6VT93bcH0uRbRZPHgzT5X2ekyf8aI03+1xhp8r/GLo8QabJQHV+l6Wl9ut3mXW+pLVqYzOOcxUvUf6O/5HS74yfcvh7JQmkV0mSDovsxkCb3a1Q0Q6SpKEH3t0ea3K9R0QyRpqIE2b4IAaSpCL3utr5J08N6MYfNe3aqb+3ZFdL57w5Zpi5ddpiXs0mDdgGkycIHpAEhkKYGFKlgikhTQYAN2BxpakCRCqaINBUEyOaFCCBNhfB1NvZFmh54Z4/avPst9f139oZULl56qJalQ9UJekGHNj6QpnZUHWnyv85Ik/81Rpr8rzHS5H+NXR4h0mShOk2Xpj/buyu4XkmWCJfHEr2gg8wqiSwdra9bavMDaWpH9ZEm/+uMNPlfY6TJ/xojTf7X2OURIk0WqtNUafqjYHGHt8Klwo8eH1eXLOnMLC0dHbNApvkhkKbm1zDNCJCmNJSa3QZpanb90mSPNKWh1Ow2SFOz69f07JEmCxVskjTt1avfiSjJzNL27pLhsqDDpUsPU/9cyxKPXgJIUzv2CKTJ/zojTf7XGGnyv8ZIk/81dnmESJOF6jRBml7VgiSiJMK0u7ts+Af1EuFyGt56vcgDj2QCSFM79gykyf86I03+1xhp8r/GSJP/NXZ5hEiTheq4LE1/NyPLhu9Ud2hZMo+zFh+iZelwdc6iQyyM3u8QSJPf9TWjQ5r8rzPS5H+NkSb/a4w0+V9jl0eINFmojovS9Pj+fYEs/emet8MR/s4hS4PrlT48hSylLTvSlJZUs9shTc2uX5rskaY0lJrdBmlqdv3SZI80paFEm7IIIE0WyLokTT/cp5cN17L0l3r5cPO4cOny4DS897Z02fAiJUaaitBrzrZIU3NqlTdTpCkvueZshzQ1p1Z5M0Wa8pJjOxsEkCYLFF2Qpj/XN6LdrFfD+0l32fBJNRIuG/6u8QkLo2xnCKSpHXVHmvyvM9Lkf42RJv9rjDT5X2OXR4g0WahOndL0x7JsuP56cv90MJIj9VLhnXssHaYOHR21MLp2h0Ca2lF/pMn/OiNN/tcYafK/xkiT/zV2eYRIk4XqVC1NB+bmOivh7X5LvTBzIBjBu8cmwpmlsZERC6MihBBAmtqxHyBN/tcZafK/xkiT/zVGmvyvscsjRJosVKcqaXrj4Gx32fBfqbcOHgwyf9/kVHCPpfP1dUs87BNAmuwzdTEi0uRiVezmhDTZ5eliNKTJxarYzQlpssuTaNkIIE3ZeCW2LluanusuGy73WJrtZvARvVz4pUsPVZ/QK+LxKI8A0lQeW5ciI00uVaOcXJCmcri6FBVpcqka5eSCNJXDlajpCCBN6TgNbFWWNG07MB3MLN2tr1kyD5EkWTb8IywbbqFyw0MgTcMZ+dACafKhioPHgDT5X2Okyf8aI03+19jlESJNFqpjW5p+PL03kKWte3eH2X1qiSwbfqh6/+QiCxkTIi0BpCktqWa3Q5qaXb802SNNaSg1uw3S1Oz6pckeaUpDiTZlEUCaLJC1JU1/8c7uQJYe3Lc3yErWvuushHeo+vXxSQuZEiIrAaQpK7Fmtkeamlm3LFkjTVloNbMt0tTMumXJGmnKQou2tgkgTRaIFpWme3e/HayE9zfdZcNlqfBLlx0eyJIsIc6jPgJIU33sq+wZaaqSdj19IU31cK+yV6SpStr19IU01cOdXjsEkCYLe0IeaZrT/Zplw39xYH+QhdyE1swsyc1pedRPAGmqvwZVZIA0VUG53j6Qpnr5V9E70lQF5Xr7QJrq5d/23pEmC3tAFml6Wy8VLqvgyczSa7OdtfDeOzGlLtGzSp/Rq+HxcIsA0uRWPcrKBmkqi6w7cZEmd2pRViZIU1lk3YmLNLlTizZmgjRZqHoaaXpR34Q2mFnSwjStZJ5JqTOmFgfLhn9yyTILWRCiDAJIUxlU3YuJNLlXE9sZIU22iboXD2lyrya2M0KabBMlXhYCSFMWWn3aDpKmp/R1SiJK39jzdrj1OYuXBDekPWvxIRZ6J0SZBJCmMum6ExtpcqcWZWWCNJVF1p24SJM7tSgrE6SpLLLETUMAaepS2vrAQ+qWzfeprVtuWMDt8s9/VT348Lbg9TNPX6tuu/HqnjZJ0vTw9L7gFLxv7dkVtl3fXTb8gywbnmbfdKIN0uREGUpPAmkqHXHtHSBNtZeg9ASQptIR194B0lR7CVqdQOul6fEnnlUXXXldsBMcd+xRC6TpSzffqV7a/looSiJQa1avVNduvDjccaLS9IBeLnzzrl+p77/TWTZcHhfrU/BkJbwT9LVLPJpFAGlqVr3yZos05SXXnO2QpubUKm+mSFNecs3ZDmlqTq18zLT10mSK2m+mad1F16jrvrBBnXLy8UFTkaxN19/eI1ciTf9l767gNLzH9AyTPJaMyLLhnXssHT027uO+04oxIU2tKLNCmvyvM9Lkf42RJv9rjDT5X2OXR4g0dauTJE3bd7ypzjn/KnX/PTep1atWBC3jr936+hvq37/ymvp5d9lwESQRJREmEScezSaANDW7fmmzR5rSkmpuO6SpubVLmznSlJZUc9shTc2tnQ+ZI00FpWnk8Z8FEU6cnFKXHX642nDY4T7sF4yhS2BsdERNTY6pvftmYOIxgaWLJ9TufQdUd2FLj0fa3qEtnhpTB2bm1MzswfZC8HzkE+OjalT/zp7e37mdBw//CEh9F+t/k/e07N/kZYdM+FfMBo4IaSooTe964in1vx+5Up2/nHssNXD/H5oy0jQUkRcNkCYvyjhwEEiT/zVGmvyvMdLkf41dHiHSNECa5K201zS5XGRyy0+A0/Pys2vSlpye16Rq5cuV0/PycWvSVpye16Rq5cuV0/PycWMrOwSQpiHSlHX1PDtlIYorBJAmVypRbh5IU7l8XYiONLlQhXJzQJrK5etCdKTJhSq0N4fWS1N0yXGzG2y87FNqw4XnhntFnvs0tXeX8mvkSJNf9ew3GqTJ/zojTf7XGGnyv8ZIk/81dnmErZcmG8VJurmtjbjEqJ8A0lR/DarIAGmqgnK9fSBN9fKvonekqQrK9faBNNXLv+29I00W9gCkyQJER0MgTY4WxnJaSJNloA6GQ5ocLIrllJAmy0AdDIc0OViUFqWENFkoNtJkAaKjIZAmRwtjOS2kyTJQB8MhTQ4WxXJKSJNloA6GQ5ocLEqLUkKaLBQbabIA0dEQSJOjhbGcFtJkGaiD4ZAmB4tiOSWkyTJQB8MhTQ4WpUUpIU0Wio00WYDoaAikydHCWE4LabIM1MFwSJODRbGcEtJkGaiD4ZAmB4vSopSQphYVm6FCAAIQgAAEIAABCEAAAtkJIE3ZmbEFBCAAAQhAAAIQgAAEINAiAkhTi4rNUCEAAQhAAAIQgAAEIACB7ASQpuzM2AICEIAABCAAAQhAAAIQaBEBpClnsS///FfVgw9vC7Y+8/S16rYbr84Zic3qJpClll+6+U519ze/F6ZM7euuXvr+s9Q5GnXrAw+pz335VrXllk3qlJOPT98hLSsnkKfGJ511SZjnV754hVp39hmV502H6QlkrfG6i65Rz7/4atjBkz+4I31ntHSOgPw+vmXzfWrrlhucy42E/CeANOWosRw4v7T9tVCU5Jf4mtUr1bUbL84RjU3qJJC1lvIPcPSXtfy8ft1H1YYLz61zGPQ9hEDWOptw5h9oOehCmtzezbLWePuON9U551+lECW36xrNLmuN4/82x7dvzsjJ9PEnnlUXXXldAOK4Y49CmtglaiGANOXALgfK131hQ/hXZ/kwb7r+dj7EOVjWvUnRWt5+13fUIz99mpnGugs5pP88dY7+RVNmI5Amt4uctcZyAH3MUUfwBw+3y9qTXdYaS/srL10fzh4yS9GgYvdJlRo2v4ZNHgHSlLF65q+T999zk1q9akWwddJrGcPSvAYCNmopf8k87QMncuBVQ/3SdpmnzvF/mJGmtLTraZenxlJT+Yt19NSt6O/1ekZCr/0I5Kmx/FHr5q/fqzZe9qngdzS/r5u/fyFNza9hk0eANGWsXp5f3Bm7oHlFBIrW0vyDzDnyFRUsZzdZ65z0jzLSlBN+RZtlrbFpH5095PNcUbFydpO1xtE/aEblmN/XOQvgyGZIkyOFaGkaSFPGwuf5xZ2xC5pXRKBILc3iAPxluqJiFegma53ji31Eu+b6lwKFKHHTrDXud3YAclxikQqGzlpj6S5eT8S4YBEc2BxpcqAILU4BacpR/KznVefogk0qIpCnlvzDW1FxLHaTp87R7jmYtliMkkJlrXFSTalzScWxFDZLjfNIlqU0CVMiAaSpRLiEHkoAaRqKaGGDrCv45OiCTSoiMKyWcg68PMyS8vGfK0qTbgoSyFrneHccTBcsQAWbZ62xtP/JY0+GC/jIH0Pu2/ojFvSpoFZ5u8haY/ncXnDex8KVbalxXvLubIc0uVOLNmaCNOWsetZ7ReTshs0qIDCollFJMn+5TEqJldUqKFTBLtLWOakbpKkg/Io2z1rj6KmYLGNcUZEKdpO1xtH7cFHjgvBr3Dy65LhJwyzwUWNadN0yAkhTywrOcCEAAQhAAAIQgAAEIACBbASQpmy8aA0BCEAAAhCAAAQgAAEItIwA0tSygjNcCEAAAhCAAAQgAAEIQCAbAaQpGy9aQwACEIAABCAAAQhAAAItI4A0tazgDBcCEIAABCAAAQhAAAIQyEYAacrGi9YQgAAEIAABCEAAAhCAQMsIIE0tKzjDhQAEIAABCEAAAhCAAASyEUCasvGiNQQgAAEIQAACEIAABCDQMgJIU8sKznAhAAEIQAACEIAABCAAgWwEkKZsvGgNAQhAAAIQgAAEIAABCLSMANLUsoIzXAhAAAIQgAAEIAABCEAgGwGkKRsvWkMAAhCAAAQgAAEIQAACLSOANLWs4AwXAhCAAAQgAAEIQAACEMhGAGnKxovWEIAABCAAAQhAAAIQgEDLCCBNLSs4w4UABCAAAQhAAAIQgAAEshFAmrLxojUEIAABCEAAAhCAAAQg0DICSFPLCs5wIQABCEAAAhCAAAQgAIFsBJCmbLxoDQEIQAACEIAABCAAAQi0jADS1LKCM1wIQAACEIAABCAAAQhAIBsBpCkbL1pDAAIQcILA1gceUp/78q0LcvnKF69Q684+Q6276Jrgva1bbljQRt571zGr1G03Xt1pMyTWSWddMnDMxx17VNDPl26+U939ze8ltt1yyyZ1ysnHq8s//1X14MPblPnZNH78iWfVRVdep848fW2YVzxQmjzWr/uouvnr94abbrzsU2rDhedm6jfNOJzYCUgCAhCAAAQqI4A0VYaajiAAAQjYIWAO6u+/5ya1etWKMKjIzzf/3wdD6RDJuOC8j6lrN14ctrn9ru+o+7b+KJSptLHichOXHnlfYr20/bW+0iNtjDTF8zKvD5KmKD0jWUl5JL2Xpd8047BTSaJAAAIQgEBTCCBNTakUeUIAAhDoEhAZMjMog6DE5WH7jjfVOedf1TPLkzaWTWlas3plMCNlpM/kJSI1TLrS5NFPmtL2izTxUYMABCAAgTgBpIl9AgIQgEDDCMRPrxuUvgjATx57MphZktkWEYfozFOWWNLPoBmeNLIhOZz2gRPVy6++oY456ojg1DmZ/ZKHvFamNKXtN804GrbLkC4EIAABCBQkgDQVBMjmEIAABKomYMQl2m/SaWrm/ei1QE/+4I6edLPGGiZNaa5pEnk59X0nBNcwST6Sn8w6/eE3/rx0aUrTL9c0Vb1H0x8EIAAB9wkgTe7XiAwhAAEI9CUQP8BPOm3PiI5ZJKJfsCyxilzTJNJkFmeQXMzsV5YZnjzXNKXtN0se7JoQgAAEINAOAkhTO+rMKCEAgRYQkNPcZOW4+GxS0rVMw3D0izVspmnY6XXm9DyRJrNqnxGwLLJSRJqG9Zslj2EceR8CEIAABPwggDT5UUdGAQEItISACNB37v+vwUxN/GFkIL6qXj9pyhPLpjRJ/nJNlVkWPYusFJGmYf1myaMlux3DhAAEINB6AkhT63cBAEAAAk0iYARIco7PKCUtMS7tBkmTrKaXJZZtaYqyzyIrRaVpUL9Z8mjSvkOuEIAABCCQnwDSlJ8dW0IAAhCojUDSjV77XbM07PS8LLGGSVPahSCSZsqyyEq/PMxphaYw0Zvbmmua4kWL98tCELXt1nQMAQhAwFkCSJOzpSExCEAAAhCAAAQgAAEIQMAFAkiTC1UgBwhAAAIQgAAEIAABCEDAWQJIk7OlITEIQAACEIAABCAAAQhAwAUCSJMLVSAHCEAAAhCAAAQgAAEIQMBZAkiTs6UhMQhAAAIQgAAEIAABCEDABQJIkwtVIAcIQAACEIAABCAAAQhAwFkCSJOzpSExCEAAAhCAAAQgAAEIQMAFAkiTC1UgBwhAAAIQgAAEIAABCEDAWQJIk7OlITEIQAACEIAABCAAAQhAwAUCSJMLVSAHCEAAAhCAAAQgAAEIQMBZAkiTs6UhMQhAAAIQgAAEIAABCEDABQJIkwtVIAcIQAACEIAABCAAAQhAwFkCSJOzpSExCEAAAhCAAAQgAAEIQMAFAkiTC1UgBwhAAAIQgAAEIAABCEDAWQJIk7OlITEIQAACEIAABCAAAQhAwAUCSJMLVSAHCEAAAhCAAAQgAAEIQMBZAkiTs6UhMQhAAAIQgAAEIAABCEDABQJIkwtVIAcIQAACEIAABCAAAQhAwFkCSJOzpSExCEAAAhCAAAQgAAEIQMAFAkiTC1UgBwhAAAIQgAAEIAABCEDAWQJIk7OlITEIQAACEIAABCAAAQhAwAUCSJMLVSAHCEAAAhCAAAQgAAEIQMBZAkiTs6UhMQhAAAIQgAAEIAABCEDABQJIkwtVIAcIQAACEIAABCAAAQhAwFkCSJOzpSExCEAAAhCAAAQgAAEIQMAFAkiTC1UgBwhAAAIQgAAEIAABCEDAWQJIk7OlITEIQAACEIAABCAAAQhAwAUCSJMLVSAHCEAAAhCAAAQgAAEIQMBZAkiTs6UhMQhAAAIQgAAEIAABCEDABQJIkwtVIAcIQAACEIAABCAAAQhAwFkCSJOzpSExCEAAAhCAAAQgAAEIQMAFAkiTC1UgBwhAAAIQgAAEIAABCEDAWQJIk7OlITEIQAACEIAABCAAAQhAwAUCSJMLVSAHCEAAAhCAAAQgAAEIQMBZAkiTs6UhMQhAAAIQgAAEIAABCEDABQJIkwtVIAcIQAACEIAABCAAAQhAwFkCSJOzpSExCEAAAhCAAAQgAAEIQMAFAkiTC1UgBwhAAAIQgAAEIAABCEDAWQL/P07D3xFchzm9AAAAAElFTkSuQmCC", "text/html": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "bio.plot_history_single_bin(bin_address=0, \n", " title=\"Reaction A <-> B . Concentrations at bin 0\")" ] }, { "cell_type": "code", "execution_count": 18, "id": "46c41c61-aeed-42f4-a4b1-6f49a5970399", "metadata": {}, "outputs": [ { "data": { "application/vnd.plotly.v1+json": { "config": { "plotlyServerURL": "https://plot.ly" }, "data": [ { "hovertemplate": "Chemical=A
SYSTEM TIME=%{x}
Concentration=%{y}", "legendgroup": "A", "line": { "color": "darkturquoise", "dash": "solid", "shape": "spline" }, "marker": { "symbol": "circle" }, "mode": "lines", "name": "A", "orientation": "v", "showlegend": true, "type": "scatter", "x": [ 0, 0.1, 0.2, 0.30000000000000004, 0.4, 0.5, 0.6, 0.7, 0.7999999999999999, 0.8999999999999999, 0.9999999999999999, 1.0999999999999999 ], "xaxis": "x", "y": [ 10, 17, 20.5, 22.25, 23.125, 23.5625, 23.78125, 23.890625, 23.9453125, 23.97265625, 23.986328125, 23.9931640625 ], "yaxis": "y" }, { "hovertemplate": "Chemical=B
SYSTEM TIME=%{x}
Concentration=%{y}", "legendgroup": "B", "line": { "color": "green", "dash": "solid", "shape": "spline" }, "marker": { "symbol": "circle" }, "mode": "lines", "name": "B", "orientation": "v", "showlegend": true, "type": "scatter", "x": [ 0, 0.1, 0.2, 0.30000000000000004, 0.4, 0.5, 0.6, 0.7, 0.7999999999999999, 0.8999999999999999, 0.9999999999999999, 1.0999999999999999 ], "xaxis": "x", "y": [ 50, 43, 39.5, 37.75, 36.875, 36.4375, 36.21875, 36.109375, 36.0546875, 36.02734375, 36.013671875, 36.0068359375 ], "yaxis": "y" } ], "layout": { "autosize": true, "legend": { "title": { "text": "Chemical" }, "tracegroupgap": 0 }, "template": { "data": { "bar": [ { "error_x": { "color": "#2a3f5f" }, "error_y": { "color": "#2a3f5f" }, "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "bar" } ], "barpolar": [ { "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "barpolar" } ], "carpet": [ { "aaxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "baxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "type": "carpet" } ], "choropleth": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "choropleth" } ], "contour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "contour" } ], "contourcarpet": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "contourcarpet" } ], "heatmap": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmap" } ], "heatmapgl": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmapgl" } ], "histogram": [ { "marker": { "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "histogram" } ], "histogram2d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2d" } ], "histogram2dcontour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2dcontour" } ], "mesh3d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "mesh3d" } ], "parcoords": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "parcoords" } ], "pie": [ { "automargin": true, "type": "pie" } ], "scatter": [ { "fillpattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 }, "type": "scatter" } ], "scatter3d": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatter3d" } ], "scattercarpet": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattercarpet" } ], "scattergeo": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergeo" } ], "scattergl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergl" } ], "scattermapbox": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattermapbox" } ], "scatterpolar": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolar" } ], "scatterpolargl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolargl" } ], "scatterternary": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterternary" } ], "surface": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "surface" } ], "table": [ { "cells": { "fill": { "color": "#EBF0F8" }, "line": { "color": "white" } }, "header": { "fill": { "color": "#C8D4E3" }, "line": { "color": "white" } }, "type": "table" } ] }, "layout": { "annotationdefaults": { "arrowcolor": "#2a3f5f", "arrowhead": 0, "arrowwidth": 1 }, "autotypenumbers": "strict", "coloraxis": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "colorscale": { "diverging": [ [ 0, "#8e0152" ], [ 0.1, "#c51b7d" ], [ 0.2, "#de77ae" ], [ 0.3, "#f1b6da" ], [ 0.4, "#fde0ef" ], [ 0.5, "#f7f7f7" ], [ 0.6, "#e6f5d0" ], [ 0.7, "#b8e186" ], [ 0.8, "#7fbc41" ], [ 0.9, "#4d9221" ], [ 1, "#276419" ] ], "sequential": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "sequentialminus": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ] }, "colorway": [ "#636efa", "#EF553B", "#00cc96", "#ab63fa", "#FFA15A", "#19d3f3", "#FF6692", "#B6E880", "#FF97FF", "#FECB52" ], "font": { "color": "#2a3f5f" }, "geo": { "bgcolor": "white", "lakecolor": "white", "landcolor": "#E5ECF6", "showlakes": true, "showland": true, "subunitcolor": "white" }, "hoverlabel": { "align": "left" }, "hovermode": "closest", "mapbox": { "style": "light" }, "paper_bgcolor": "white", "plot_bgcolor": "#E5ECF6", "polar": { "angularaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "radialaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "scene": { "xaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "yaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "zaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" } }, "shapedefaults": { "line": { "color": "#2a3f5f" } }, "ternary": { "aaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "baxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "caxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "title": { "x": 0.05 }, "xaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 }, "yaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 } } }, "title": { "text": "Reaction A <-> B . Concentrations at bin 0" }, "xaxis": { "anchor": "y", "autorange": true, "domain": [ 0, 1 ], "range": [ 0, 1.0999999999999999 ], "title": { "text": "SYSTEM TIME" }, "type": "linear" }, "yaxis": { "anchor": "x", "autorange": true, "domain": [ 0, 1 ], "range": [ 7.777777777777778, 52.22222222222222 ], "title": { "text": "Concentration" }, "type": "linear" } } }, "image/png": "iVBORw0KGgoAAAANSUhEUgAAA00AAAFoCAYAAACL/RCAAAAgAElEQVR4Xu29C5gd1X3gefqlVuuJBBIIyeAXBgeJIDuITMznsIZdx2I8zDAGsiATP2QwCZkJMDb+ljX5Yg9ZHgtkNgQMwYbBMgmQkJAEObMLXsLiBzi8X8IY29iAQAIJPVutVrf2/Oveul1dXbdvPU7VPefU7/I1996q8/if379uq359qs7t2a8figcEIAABCEAAAhCAAAQgAAEIJBLoQZo4MiAAAQhAAAIQgAAEIAABCLQngDRxdEAAAhCAAAQgAAEIQAACEJiGANLE4QEBCEAAAhCAAAQgAAEIQABp4hiAAAQgAAEIQAACEIAABCCQjwAzTfm4UQsCEIAABCAAAQhAAAIQqAkBpKkmiWaYEIAABCAAAQhAAAIQgEA+AkhTPm7UggAEIAABCEAAAhCAAARqQgBpqkmiGSYEIAABCEAAAhCAAAQgkI8A0pSPG7UgAAEIQAACEIAABCAAgZoQQJpqkmiGCQEIQAACEIAABCAAAQjkI4A05eNGLQhAAAIQgAAEIAABCECgJgSQppokmmFCAAIQgAAEIAABCEAAAvkIIE35uFELAhCAAAQgAAEIQAACEKgJAaSpJolmmBCAAAQgAAEIQAACEIBAPgJIUz5u1IIABCAAAQhAAAIQgAAEakIAaapJohkmBCAAAQhAAAIQgAAEIJCPANKUjxu1IAABCEAAAhCAAAQgAIGaEECaapJohgkBCEAAAhCAAAQgAAEI5COANOXjRi0IQAACEIAABCAAAQhAoCYEkKaaJJphQgACEIAABCAAAQhAAAL5CCBN+bhRCwIQgAAEIAABCEAAAhCoCQGkqSaJZpgQgAAEIAABCEAAAhCAQD4CSFM+btSCAAQgAAEIQAACEIAABGpCAGmqSaIZJgQgAAEIQAACEIAABCCQjwDSlI8btSAAAQhAAAIQgAAEIACBmhBAmmqSaIYJAQhAAAIQgAAEIAABCOQjgDTl40YtCEAAAhCAAAQgAAEIQKAmBJCmmiSaYUIAAhCAAAQgAAEIQAAC+QggTfm4UQsCEIAABCAAAQhAAAIQqAkBpKkmiWaYEIAABCAAAQhAAAIQgEA+AkhTPm7UggAEIAABCEAAAhCAAARqQgBpqkmiGSYEIAABCEAAAhCAAAQgkI8A0pSPG7UgAAEIQAACEIAABCAAgZoQQJpqkmiGCQEIQAACEIAABCAAAQjkI4A05eNGLQhAAAIQgAAEIAABCECgJgSQppokmmFCAAIQgAAEIAABCEAAAvkIIE35uFELAhCAgLUE/s8b71S33vld9d3vXKkOW3qwtXESGAQgAAEIQMAVAk5J0y9fe1N94uxLEtn+73/0afW//vuTus79zC/+iXpj0xb1L/f8t0mxhLEvP+o96s5v/HHX40wK4PuPPqPO/fI1wS5beJoG9dun/Wf11pZtic1WcYIZnszGA6iib9Ms07TX7vOQpm6aMi5/3tKMT8qEn8ssn8m00lT27yXJz7Mbft4a6v/y28ep6/7kD9IOnXIQgAAEIAABawg4KU3xf3jDEwQb/kF2+STuwj/+C/V//8uPg4PTtNyFOfrsmZ9Q/+X8M1N/AMITRlO5FWmSR1xqw5O7suQlKvzxPkLuWdmkhtjFgkhTcfiuSlPSZ+3oEz+jTH2Wi5OlBQhAAAIQgEB6Al5Ikww3POl97sHb0o++hJJlnySWEHKryfCERjaIPJkQiL/6+wfUf/2zbwd9ZPlLeXSc0Rmwoidc7aQpz4lplly06zdsQ/r/4WPPZxLKLP13q2zZn4ey2+8Wt6TjP8vnJ+1MU1njCz/3N191sfrIqhWtbsLtJn63lBU77UIAAhCAAASSCHgvTUmX9MX/IY/OsEQhtfuHPekSLyn7pa9/Y9KlKGFbYTvt/sqa5hKWsO6qlUe1JETaLyoRYYzRkxzZJpfpFZn5MCFL8QPWhDy1k5cw3iwnpml/peRtW3IefcRji4reo09saM0STieo7Y7d6H0v8TLx4yB6Qv7pP/zTSZc7Rj8z8eO63efhU6d8tHVZ6EEL5wezgElxhvvCdtK2H78crBNXaT/L5y16rIexpTmOkurFcxc95qPHQqfP/XQ5iv9hKf57Kc9xFf8shL9T432Fv4+L/G5J+7mjHAQgAAEIQMAkAW+kKelkODwpif4DnXQCK//ALz3koEl/5Q//0Y+eBIb/4MdPWKSsyIzcUzXdX76TpEm2xS+Fa7ctLkkmZ0ficXeaGWl3EJYhS9PJU9aTr3bjyjveNB/GpGNpunphXqNjm25b/EQ7aZYhzbGbdH9LUr3ofVlxSYrfz9fp8yCxJ10KKvn406+snTRLkZSjTu1HP6tpuYbSlObzFrKI/iEmbb7ls/L3//zwpHsck34/5fmcT5cjuccoKjPtpCnNcdXuOM6SlzSfIcpAAAIQgAAEuk3AC2lqd5LSbmZHyj/+zE+m3NcSTUZ4shj9i3GaS4GynCy0u4Qm6SSp3VjkRPKQxQsLLS6R9NffpJPB6Q7W6Ixemr+ymzjw8yxcMd1CEJ3+ep835qyXjrY7huLHS7uT6aR8Fjl245dUtTtuky7JyvJ56MQ36dKuLO2n5RpKU9LxEP+8tWtTYv3IcctzrVwnfXxoxQdaCyYUkab4bHnSsTHdTFN0cZ0ss0TT/RGC+5o6HenshwAEIAABGwk4KU1xkEknN9OdaIQnfdG/tra7DCb61/40/9ibOIlLOmkrU5qyzEy0O4jzSlP8UilpP375ZLs+80qTtBdfCCJkHr8EzMSHNqs0tct1/JjOIk1pjt1O/YZ5qUqa2q00GD0+snzeOo0vKvtpP2/hH2yKLJyS9BmItmdSmuR4josf0mTiU04bEIAABCDgOwEnpSkqSeEJaXx2o939AtGEhn+FDWcfom3E/6qa9q+sWU7ippslil+il/YkLs8BO93si7SXZXENVy/Pk3HmvfeoE/MsK/NNd5zFLy9LK01pjt3plvMPxxd+PsqWpnYrDRaZycrCNemPFiGDpM9s0j2RaS4bDevF/+gT78NFacrye7DT54f9EIAABCAAARsIOC9N4V9O5bt3on+BDk80Os1apD3xnO5EKprILCcLWcqWJU3TnZAVkYgy5KnMhSAkh0n3vJj4kGblmHZGxPSxK/2mOdkvW5rKaj8t16zSFD9Gwj9CdPrdk/YzbVqa4v2WMdPEQhAmfnPQBgQgAAEI2ETAC2mK/mV6uhuck8C3Wxo3y30hUvb7P342WAhiuvul4icnNtzT1Omm9aRFKbIcwCbkyYQsRWcK5HXS5XlZ5SYLh/ASrHazdnIM3fUPDwaLkaS99yaLNLVrM3rsprnvScacRWqyfB5Cnu1OuJM+q1naT8s1izRJ//HV+dJIznQzX/GZpqT7Kzsde0V+t2Q5rtrFwZLjnTLEfghAAAIQcI2AF9IUnSWI3pOStHqelJWTEFkqWU6ck1YHi0pY0gpm0ctp4ic07U4W2p2Iycl0/D6adqvnpbkxPesB2EmKOklV2v5s/3Lb8EQx7b0p4WxC2u+byfLltmEsaVd+i1+amnRCntRm/NhNWj0vKpvf/vP/LVjYIIs0Zf08SH9J8hqV7+gMTpb203LNIk3hcRCV4bSfmaTFEsL24sdh1gVf2i3kkuZ3iwlpEoZ8uW3a346UgwAEIAABFwh4I03Rky15HZ7MtrtXIypD8TIiMXKC+ImzL0m8XCl+43b8kqb4PQ4mv6cp/lftrCdT0YMyzcxKWZesdevDMd39W2kuTYtKhFwSmlaawnrtvhMsaQYqfpxN9z1NaVc563Tshie7MrboI3oSn0WapI2sn4f4Z1nei0CG31HW6XvWOn3eOnHNIk1h2fjxnPY+wPjxKLHLH3TiK2LGf0d1WukxzJH8LovmMqleGZfnhTzSfAddt34X0C8EIAABCEAgCwGnpCnLwCgLAQhAAAIQgAAEIAABCEDABAGkyQRF2oAABCAAAQhAAAIQgAAEvCWANHmbWgYGAQhAAAIQgAAEIAABCJgggDSZoEgbEIAABCAAAQhAAAIQgIC3BJAmb1PLwCAAAQhAAAIQgAAEIAABEwSQJhMUaQMCEIAABCAAAQhAAAIQ8JYA0uRtahkYBCAAAQhAAAIQgAAEIGCCANJkgiJtQAACEIAABCAAAQhAAALeEkCavE0tA4MABCAAAQhAAAIQgAAETBBAmkxQpA0IQAACEIAABCAAAQhAwFsCSJO3qWVgEIAABCAAAQhAAAIQgIAJAkiTCYq0AQEIQAACEIAABCAAAQh4SwBp8ja1DAwCEIAABCAAAQhAAAIQMEEAaTJBkTYgAAEIQAACEIAABCAAAW8JIE3eppaBQQACEIAABCAAAQhAAAImCCBNJijSBgQgAAEIQAACEIAABCDgLQGkydvUMjAIQAACEIAABCAAAQhAwAQBpMkERdqAAAQgAAEIQAACEIAABLwlgDR5m1oGBgEIQAACEIAABCAAAQiYIIA0maBIGxCAAAQgAAEIQAACEICAtwSQJm9Ty8AgAAEIQAACEIAABCAAARMEkCYTFGkDAhCAAAQgAAEIQAACEPCWANLkbWoZGAQgAAEIQAACEIAABCBgggDSZIIibUAAAhCAAAQgAAEIQAAC3hJAmrxNLQODAAQgAAEIQAACEIAABEwQQJpMUKQNCEAAAhCAAAQgAAEIQMBbAkiTt6llYBCAAAQgAAEIQAACEICACQJIkwmKtAEBCEAAAhCAAAQgAAEIeEsAafI2tQwMAhCAAAQgAAEIQAACEDBBAGkyQZE2IAABCEAAAhCAAAQgAAFvCSBN3qaWgUEAAhCAAAQgAAEIQAACJgggTSYo0gYEIAABCEAAAhCAAAQg4C0BpMnb1DIwCEAAAhCAAAQgAAEIQMAEAaTJBEXagAAEIAABCEAAAhCAAAS8JYA0eZtaBgYBCEAAAhCAAAQgAAEImCCANJmgSBsQgAAEIAABCEAAAhCAgLcEkCZvU8vAIAABCEAAAhCAAAQgAAETBJAmExRpAwIQgAAEIAABCEAAAhDwlgDS5G1qGRgEIAABCEAAAhCAAAQgYIIA0mSCIm1AAAIQgAAEIAABCEAAAt4SQJq8TS0DgwAEIAABCEAAAhCAAARMEECaTFCkDQhAAAIQgAAEIAABCEDAWwJIk7epZWAQgAAEIAABCEAAAhCAgAkCSJMJirQBAQhAAAIQgAAEIAABCHhLAGnyNrUMDAIQgAAEIAABCEAAAhAwQQBpMkGRNiAAAQhAAAIQgAAEIAABbwkgTd6mloFBAAIQgAAEIAABCEAAAiYIIE0mKNIGBCAAAQhAAAIQgAAEIOAtAaTJ29QyMAhAAAIQgAAEIAABCEDABAGkyQRF2oAABCAAAQhAAAIQgAAEvCWANBlI7etvDxtohSZsJDCjv1fNmz2g3to2YmN4xGSIwJKFQ+qNrcNq/35DDdKMdQQWzp2hdo+MqT17x6yLjYDMEJg9s1/19/WobbtGzTRIK9YRkPwunDuoNr2zx7rYygzo0AOHymyetlMSQJpSgpquGNJkAKKlTSBNlibGcFhIk2GgFjaHNFmYFMMhIU2GgVrYHNJkYVJqFBLSZCDZSJMBiJY2gTRZmhjDYSFNhoFa2BzSZGFSDIeENBkGamFzSJOFSalRSEiTgWQjTQYgWtoE0mRpYgyHhTQZBmphc0iThUkxHBLSZBiohc0hTRYmpUYh1V6abrnjPnXdzXdPSflzD97W2nbel69RDz/6TPD+hFUr1E1XXTypPNLk7ycGafI3t9GRIU3+5xlp8j/HSJP/OUaa/M+xzSNEmrQ0/fjJDVNEKEza1667Xb22cXNrvwjU0iWL1GUXntPKK9Jk8yFeLDakqRg/V2ojTa5kKn+cSFN+dq7URJpcyVT+OJGm/OyoWZwA0tRBmlavuURd/pW1auXyIwLaTzz7krr0ilvU+nVXIk3Fjz/rW0CarE+RkQCRJiMYrW4EabI6PUaCQ5qMYLS6EaTJ6vR4HxzSFLs87/BlB7eEaOOmLerkMy5S9991rVqyeGFwMCRtY6bJ388J0uRvbqMjQ5r8zzPS5H+OkSb/c4w0lZvjpKupyu3RXOtJ5+fmWm+0VHtpigOVA0Yect9SGmlafPVi9af/09XqzF87y3RuaM8CAn29PWpwRp/avWefBdEQQlkE5gwNqJ179He78D1NZSHuertDg31qdN9+tW9svOuxEEA5BAb09+r16t/ZI3wXVzmALWhV8juk/03eVbN/k+fOGjBCX245ufPe701qKzpZ0A1pkiu41lxwuVp3/aWtq7ryDBZpykOtYJ0webIQRBpp6vmTnqDH4w49Xv3Rqv+iPnnEqQUjoLpNBJAmm7JRXixIU3lsbWkZabIlE+XFgTSVx9aWlpGm/Jk4+sTPqKgghS2JSB168IFq7VmnqG5IU/4RTa6JNJkimaGdqDRJtU73NN382M3qTx+6Qr2y/edBLx9ZdqL6/Q9dpE487OQMvVLUVgJcnmdrZszGxeV5Znna2BqX59mYFbMxcXmeWZ42tsblefmyImL0o8eem3Q/flJLoTTJvnBGqp1oRWesorexyHnzaas/qu5Z/5B65dU3g26u/ur5asnBC4MZpfAR1kmSnfiMmNRffdLxKmmmLFztGmnKd2xkqiXJjS7qIO9/88NHt1bHS7t63l8++efqhsevVZt2Nw6Qk9/9iUCejj/0I5niobBdBJAmu/JRVjRIU1lk7WkXabInF2VFgjSVRdaedpGmfLmQWaYLzz09mE2a7hF+xU4oKVK203nx+gceUV/6+o0qlBcpL7IUSlG4P34ZoLSddCtMXPBEhv7yO/8UnJfLvi+c/W9b6wx0uqUmH632tWp/T1P0O5gE05mnfmzScuKyLe33NI2O7Q3E6YbHr1M7R3cE1P/dEZ8K5GnFomNN5472KiCANFUA2YIukCYLklByCEhTyYAtaB5psiAJJYeANGUHHM7ApLlnKOnyPPk+0/CredrN5ogoXfDZ04LZoHCmKRS0pDrSpsxEyaRFdL+MThZgSxOrlBUhu/7We6a0Ey7elp3W9DVqL00mgMZXz9s+sk39xWPXqBueuFaN72/cdHzGBz+tvrjyj9SRCz9ookvaqIgA0lQR6C53gzR1OQEVdI80VQC5y10gTV1OQAXdI03ZIZuUpvAWlqQowtmpdtIUFaF2svPGpreDS/jCWaukfsKZrOi+dusQZKeFNJlmNqW9dkuOv7Hz9UCcvvnUDa06Zx/92UCe3ntA43ufeNhNAGmyOz+mokOaTJG0tx2kyd7cmIoMaTJF0t52kKZ8uclyed7SJYsmXXEVnWmK3/ffTmrknqb4TJMJaZJxRK8Ii14ayD1N+Y6Nymt1+p6mX2x7Wd34+J+pdc99sxXbOcu/oL74oT9Sh897T+Xx0mF6AkhTelYul0SaXM5eutiRpnScXC6FNLmcvXSxI03pOMVLdVoIQsSo3ep5SZfnTXf5XJGZJom73eV5ScKGNOU7Hrpaq5M0hcG9tGWD+saT/0399fP/vRXv53/994OZp0PnLOvqGOg8mQDSVI8jA2nyP89Ik/85Rpr8zzHSlD/HSUuOhyISLhLR6Z4m6T1cwS56CZ2I1YeP+UDwPUtFpEnuRZIYfvX6ptYibeFCELIARFyoZEzy4PK8/MdF5TXTSlMY2AtvPau+8cSfqb958Y5WrOcd+5+DmafFsw6uPH46bE8AaarH0YE0+Z9npMn/HCNN/ucYaSqW46Qlu6OzRmmkKSpO0Wiiq+flvTwvXMAhvkhbGKPI2XU3393qVu6jClfu4/K8YsdGZbWzSlMY2DObnwzk6e9/clewqa+3T33h1/9QnXvsH6qDZy+pLH46QprqfgwgTf4fAUiT/zlGmvzPMdLkf45tHiGr5xnITl5pCrt+4s0fqxu1PN33079rRXPuyv+kztUCtWTOUgMR0kReAsw05SXnVj2kya185YkWacpDza06SJNb+coTLdKUhxp1TBFAmgyQLCpNYQiPv/GouvmpP1f/+NLftqJa++t/oGee/pNaOvddBiKliawEkKasxNwsjzS5mbcsUSNNWWi5WRZpcjNvWaJGmrLQoqxpAkiTAaKmpCkM5ck3/1Xd/OSfq3tfmrhu83PHnB/I07vmHW4gYppISwBpSkvK7XJIk9v5SxM90pSGkttlkCa385cmeqQpDSXKlEUAaTJA1rQ0hSE9vfkJdfMT/5f6u5/c2Yry08vXKllx74gFRxmInCY6EUCaOhHyYz/S5EcepxsF0uR/jpEm/3OMNPmfY5tHiDQZyE5Z0hSG9uxbTwXy9Lcv/lUr2k8ddXYgT8csWmlgBDTRjgDSVI9jA2nyP89Ik/85Rpr8zzHS5H+ObR4h0mQgO2VLUxjihrefV9966i/Ud56/tRX1Ke//D4E8Hb/kIwZGQhNxAkhTPY4JpMn/PCNN/ucYafI/x0iT/zm2eYRIk4HsVCVNYaivbP+5uvXpG9U3n7pBje8fDzZ/7PCPa3n6A3XiYScbGBFNhASQpnocC0iT/3lGmvzPMdLkf46RJv9zbPMIkSYD2alamsKQNw9v0jNPNwTytGt0Z7D5uCX/Rv3einPVf/jAmQZGRhNIUz2OAaTJ/zwjTf7nGGnyP8dIk/85tnmESJOB7HRLmsLQd47uCOTp1mduUpt2vRFsfv+CIwN5kp++nj4Do6xnE0hTPfKONPmfZ6TJ/xwjTf7nGGnyP8c2jxBpMpCdbktTdAi3P/uX6r9reZL7n+SxcOhALU7nqd9bfq5aNGuxgdHWqwmkqR75Rpr8zzPS5H+OkSb/c4w0+Znj1WsuCQa2ft2VVg8QaTKQHpukKRzOP/30HnXbMzerH772UGuEa47+nDp7+edZcS9DzpGmDLAcLoo0OZy8lKEjTSlBOVwMaXI4eSlDR5pSgnKo2BPPvqS+cfs/qF+9vkld/pW1auXyI6yNHmkykBobpSkc1vdffVDPPN2s7nv571sjlcUizj7682r1+041MHq/m0Ca/M5vODqkyf88I03+5xhp8j/HSJN/Of7adber3zjmSPWvT78YDO6yC8+xdpBIk4HU2CxN4fBe3PKCuuO5b6l1z35L7RkbDjbLfU8y+3TW0Z9VswfmGCDhXxNIk385TRoR0uR/npEm/3OMNPmfY6SpWI7/cdv2Yg0UqP3J+fMSax994mfU/Xddq97Y9La69IpbrL5ED2kqcACEVV2QpjBWWWXvjuduVeu0QP10a8PqB/tmqjM+uEadcdQa9aFDVhkg4k8TSJM/uZxuJEiT/3lGmvzPMdLkf46RpmI5/t1f/FLdufWdYo3kqP3X7z5MnbnggCk1w0vzbrrq4mCfCNS66y+19hI9pClH8uNVXJKmaOzrX75Xfee5b6oHf3l/a/NvLPlNLVCfVmccebYa6JthgI7bTSBNbucvbfRIU1pS7pZDmtzNXdrIkaa0pNwthzQVy51I0+7xxvd7Vvk4WwtTkjSFl+atPun4IBx5Lw9bL9FDmgwcNa5KUzj0DW8/p+7c8G1194Z1auvwlmDznBlzg9mnU99/uhKRqusDaapH5pEm//OMNPmfY6TJ/xwjTX7lWGaWkh7PPXiblQNFmgykxXVpiiK464Vvq7s2fGfSqnsfPGi5OvWI0wOBOmz+uw0Qc6cJpMmdXBWJFGkqQs+NukiTG3kqEiXSVISeG3WRJjfylCZKuTQv6R4mmy/RQ5rSZLZDGZ+kKRzqM5ufVH/34p3q3pfuVm/ser1FQFbe+533/jt10uG/ow6du8wAPbubQJrszo+p6JAmUyTtbQdpsjc3piJDmkyRtLcdpMne3GSN7LwvX6OWLlk05VI8my/RQ5qyZjmhvI/SFB3m//OL9eren9wdCNT4/olrYeWyPZGnk9/zCfVrB64wQNK+JpAm+3JSRkRIUxlU7WoTabIrH2VEgzSVQdWuNpEmu/JRt2iQJgMZ912aQkQ79m5X8qW5D/zin9X9r/yzGh3b26L3vgM+oE56jxYoLVEfWXaiAap2NIE02ZGHsqNAmsom3P32kabu56DsCJCmsgl3v32kqfs5qHMESJOB7NdFmqKo9o2Pqvu1PD3wynfV/T//Z7Vp9xut3QuHDtQzUJ9QJ7/7d9RJ+meof5YByt1pAmnqDveqe0WaqiZefX9IU/XMq+4RaaqaePX9IU3VM6fHCQJIk4GjoY7SFMf2yMbvqwe0PIlEbXj7+dbuHtUTiNPJ7/5E8HzoHLfug0KaDHxAHGgCaXIgSQVDRJoKAnSgOtLkQJIKhog0FQRI9UIEkKZC+BqVkabJEF98+4VAnuQyvh+9/vCknXLv0/FLP6KOX/IRterQ31IHz15iIAPlNYE0lcfWppaRJpuyUU4sSFM5XG1qFWmyKRvlxII0lcOVVtMRqFyaVq+5RL3y6puJ0dm6LnsnlEhTe0Kbdr+p5akhUPH7oKTWBxZ+UK1a8luBSK065LfUsnmHdcJd6X6kqVLcXesMaeoa+so6RpoqQ921jpCmrqGvrGOkqTLUdJRAoFJpare8oOuZQZrSZ1Au43vkte+rR4PnH6jd+3ZNqvye+e9vzkQ1ROrwee9N33gJJZGmEqBa2CTSZGFSDIeENBkGamFzSJOFSTEcEtJkGCjNZSJQqTTZ/IVVmajFCiNN+ek9/sajKhQped4+sm1SY/JlujIDtUoL1IcXr1JHHXR0/s5y1ESackBzsArS5GDSMoaMNGUE5mBxpMnBpGUMGWnKCIziRgkgTRGc6x94RH3p6zeqdddfqlYuP6K1R2bIHn70meD9CatWqJuuunhSEpAmc8fkU5seV4++/oPmTNT31dt73prU+Kz+2eqYxSvVisXHquUHHRu8lkv8ynogTWWRtatdpMmufJQRDdJUBlW72kSa7MpHGdEgTWVQpc20BCqVJpGPUz9+glp90vFp46usnAjT9bfeE9xvFZUm+Wbi1zZubolS0iWGSFN5aXr+7Weal8QmwoQAACAASURBVPP9QD315uPqle0/m9LZ3Bnz1PJFWqAWNWRKREq+N8rEA2kyQdH+NpAm+3NUNEKkqShB++sjTfbnqGiESFNRgtQvQqBSaXri2ZfUpVfcotavu7JIzMbrhsIkccUvIZSFKy7/ytrWzFPSGJAm4ylp2+CW4bfV05ufUM/Iz6Yngte/2v7KlPLzBudrifqQWiEypSXq6IOOUe9bkF2kkKbqctvNnpCmbtKvpm+kqRrO3ewFaeom/Wr6Rpqq4VxVL0mLw9m8KFyl0iRCMt2jG6CiwiSxRaVp46Yt6uQzLlL333WtWrJ4YRB60jakqaqPV3I/m4c3BQIVStTT+vXrO1+dUli+ZPeoA49WH9T3RR25oPF81EHL1YEzD2o7AKSpu7mtqnekqSrS3esHaeoe+6p6RpqqIt29fpCm7rEvo2eRpgs+e1rrCrT41V1l9FmkzUqlqUigZdSNC1NeaSojNtosRuD1Ha+rxzY+ph7f+Lh67PXH9KV9T6lfbvtlYqPL5i1TKw5eoS/xWx48r1isf/RzX09fsSCoDQEIQAACEIAABCCQSCAuTUnn5Tahq7U0idHeee/3EvNx9VfPVytXHMFMk01Ha8FYtg5vURu2PKs2vP28/nlWvfD2c+rFLc+rnXt3TGm5p6dHHbnw19TRi1aolUuOUYcMvVffJ3WEev+CI1VvT2/BSKhuGwFmmmzLiPl4mGkyz9S2Fplpsi0j5uNhpqkY03988R+LNVCg9ieP/OSU2nFpknUDjjv2KLX2rFMK9FRe1cqlSe4JWnPB5ZNGFF+trrzhdm6Ze5o6M/KtxC+2vaxeeEsE6rlApDbon59ufTFxmCJTsshEKFDvW3BE8F5k6oCZC3xDU5vxIE3+pxpp8j/HSJP/OUaaiuX4d//md9Wdz91ZrJEctf/6P/61OnP5mYnSJAuwRR8Xnns60iRAwiW9k+4RkpkdG1bVi0sTq+fl+HR4UGVkbE8gUD/d+rz6+fYX1bNvvqBf/0SJYLV7LJ51cEOotEiJRMnCEyJXy+YdzqV+lh8TSJPlCTIQHtJkAKLlTSBNlifIQHhIUzGIIk27R3cXayRH7bNXnN1WmqL3NEnTch5uixPEh1rpTFN8Gi4MxqZrGJO+gJfvacrxCfGkSnwhiNGxverld14KZqJe1hIlIvXyO/LzUuJlfiGGQ+cs0/J0mFo2R/9oiWq8flfjee5harBvpifE3BwG0uRm3rJEjTRloeVmWaTJzbxliRppykLL/rJJXmDzJXqVSlOSkEhKw0v2urF6nolDitXzTFC0s40sq+e9tuNXDYHa2pQqLVI/0z+yvdNDZqmWanl613wtVCJW+rUI1eJZh6iDZx+iFun9Pfo/HuUQQJrK4WpTq0iTTdkoJxakqRyuNrWKNNmUjeKxxKUpXKGamSbN1oWZpjyHANKUh5obdbJIU7sRjY2PqVd3/lLL0y/Vq9v1jzxrkQqe9fvX9L594/s6AhGxWqwFSp4XaZma9H623jfU2D97YE7HtigwmQDS5P8RgTT5n2Okyf8cI01+5Tjpe5psFSYhX+lMkwv3NOU5HJGmPNTcqGNCmtKMdOPO15oy1ZSqpkxt2vWm2rz7TbVJ/6R9zBmYOyFXIlMiW03JCmRLtsns1dDitE16Xw5p8j7FCmnyP8dIk/85Rpr8z7HNI6xUmgSE7avn5UkW0pSHmht1qpKmTjTG94+35GnTrje0RL0Re6/FSgRr+E19k+euTs0F+2Xp9GDmSs9QyeV/wevmTJZ84e/8mQfoFQEXqvmDB6gFgwvVvMH5qdp1sRDS5GLWssWMNGXj5WJppMnFrGWLGWnKxovSZglULk1mw7ejNaTJjjyUEYUt0pRlbDtGtqtNwyJWIlL6OZyt0kLVEi697a3hzVmaDe6pkmXVWyKlX8+fIe/1s/5ZMKiftWCF++V5gezT22xf6AJpynQoOFkYaXIybZmCRpoy4XKyMNLkZNq8CRppMpBKpMkAREubcFGa0qKU+6jCS/8CmdJStTmYxXpTvT38lto28o56Z88W9c7I1uBHZCzvY6h/VlOotFw1Z68CuZohktUQrwO0dM3TgjW7f7aaM2OumjWj8Tx7YLaa2TeUt+tU9ZCmVJicLoQ0OZ2+VMEjTakwOV0IaXI6fc4HX4k0hWuuf+nrN04LjNXznD+evBuAz9KUNVmyoEUgUiNbguete7RMaakK5Wrb3sa2bbJdS1a4fat+LUu1F3n09/arWVqeAonqnxOI1OymUAWCpUVLFsCYM0Pvkx8pI8/NbbP0cyhjsk329fX0tUJCmopkx426SJMbeSoSJdJUhJ4bdZEmN/Lka5SVSJOv8MJxMdPkb4aRJjO53b1vl3pnuClSWrreCWexRLICAdPypaVr+8g2tWt0Z/CdV1JHnnft3aXky4ZNP2T2K5AvLVEHzJqn+tWgGuwf0rNag2pm8Dyk38/Ur/WPfj1zQL+Xffr1kH49U3+3llx2KNvltdQJ3gd1JrbNHJipZvQOmg6f9jISQJoyAnOwONLkYNIyhow0ZQRGcaMEKpWmdt/TZNOX2+ahizTloeZGHaTJjjztGxtVu0KJ0lK1a6/+2ad/9KIXgWDp50C2ZJ/8jGjZkm26zE5ddncgYs19Uk+/lsU1qnrIohsNqWoKmUiVliuRrEEtVAN9M9SMvgE10Duj8br5PNA7oLfrbfp98By8jm3T76N1GuUnb2vUb7Qftidl+vr6lMzi9arGs88PpMnn7DbGhjT5n2Okyf8c2zxCK6SJL7e1+RCpd2xIk7/5D0RLi5gI1cyZo+q1rdvUnn3D+mePGh4dViPje9SeUf0zpl/rma49etseeW6WmXgt+0ea25uvdR0pP6zLyr69+seFR1SgRKjkEsb+nn7Vp4UqeC3Prdeyv7Fdtsm+oH5YTraHr/t0uaaYhfUn2mvWj7YVtDnRbih1E+1P7JN2w/6DOrqd/uA5fN2I6cC5Q2pk7361T38lWo+W2N7enmAFyfBHb9XyOPFetks5WQAlKCP/9Tb29/TI1kh9vS9os1le9vOongDSVD3zqntEmqomTn9RAlZI0y133KfuWf+QWr/uSiezw0yTk2lLFTTSlAqT84XKvqdJZrVEvgIJE5GaJGQjwT1fo+OjWq70s/7ZO75Xf+Hx6MQ2/T7c1yoXbpNyzTrBc7Ot6LbRoKy0F+5v9Dc2NhZ8sfK4ajzzMEMgFKpQpEIxm9geEy5RsIjAteRsf2x7UKaxrUfvi0pfTyh0EckLpS+UvMnPIo8N+RNhDP7T76NSGG4LyjTLttqI1JmunIwl6CEYS7OveFsJ74OYZHuzXpRd2F+rTLPsYH9fIMMje8daiWyMb+IRFdrp9jVqtK/b2DuxPy7KndqeNi495nYxd4xrfyOvucbcoa6K7Z/Uj94nx0jbuKXuJF6TP2vT8or0K/mdM7Nf7Rie+vsq3sbEQdCIK9zfijvSbuK2MN4O5RqNN8c+HaOkXy8JzJOKfe64M8z8cqKVQgRKl6ak72VKinjd9ZeqlcuPKDSYblVGmrpFvvx+kabyGdvQQ9nSZMMY08QQFSgRqrH9Wqb271NjWqiC1/Lcei37G9tlm+wL6oflZHv4ekyXa4pZWH+ivWb9aFtBmxPthlIX9B/Eo7c0+5N2W6+ljt6/L3gOXzdi2i/ltLzKvv26/vj+/cElmvKzX/6T1/Jfc1uwvbkvuq2xPagxUVbXC+q32pIzSB4QgAAEzBDY/8f8TjFDslgrpUtTNLx29zQVG0L3ayNN3c9BWREgTWWRtatdpMmufJQRTZX3NIVCFchZU8SikiXS1RIxEa5Q2kJhC6ROb+9pylxE5ET2lAhbXPB0nSSRa7TdkLzJz7qV5nZpL7p/vNlWq7xMMcTbiNVptZXQl8Qc/hf0lRjP5O2tcgGDifgmi+nkff39MlOm9IzuxP2KUjP6kLYmvY/sj++TGNLWlXLR+sHryMTLdP3G6wbv9bjDh8m4JKYpcUYGOS2vWN2pcU/DS+pO4TnRcdo8SX5nDPSp4ZGEmfGE+IIemjmP9hEyCLe1mETaaJWPboscP5NiDrd3ZJTw261d3JGi3/30P5Txa5E2MxKoVJoyxuZMcaTJmVRlDhRpyozMyQpIk5NpyxR0ldKUKTAKGyPAPU3GUFrbEPc0WZuaWgSGNBlIM9JkAKKlTSBNlibGcFhIk2GgFjaHNFmYFMMhIU2GgVrYHNJkYVJqFFKl0rRx0xZ18hkXtcXLl9vW6MhzZKhIkyOJKhgm0lQQoAPVkSYHklQwRKSpIEAHqiNNDiTJ4xArlabVay5Rp63+qPrwMR9Ql15xS2u1vPO+fI069eMnqNUnHe8kamaanExbqqCRplSYnC+ENDmfwo4DQJo6InK+ANLkfAo7DgBp6oiIAiUSqFSawoUgDll8oPr8RVe2pElW2ItKVInjLaVppKkUrFY0ijRZkYbSg0CaSkfc9Q6Qpq6noPQAkKbSEXe9A6Sp6ymodQBdkSZZWlwEKrwcjy+3rfUxaPXgkSar02MsOKTJGEprG0KarE2NscCQJmMorW0IabI2NbUIrFJpksvwjjv2KLX2rFNU9DVfbluLY83JQSJNTqYtc9BIU2ZkzlVAmpxLWeaAkabMyJyrgDQ5lzKvAq5UmuLkZLYpfNx/17VqyeKFTsLl8jwn05YqaKQpFSbnCyFNzqew4wCQpo6InC+ANDmfwo4DQJo6IqJAiQS6Kk0ljqvSppGmSnFX2hnSVCnurnWGNHUNfWUdI02Voe5aR0hT19BX1jHSVBlqOkogUKk0hQtByD1NPj2QJp+yOXksSJO/uY2ODGnyP89Ik/85Rpr8zzHS5H+ObR4h0mQgO0iTAYiWNoE0WZoYw2EhTYaBWtgc0mRhUgyHhDQZBmphc0iThUmpUUiVSpPr38fU7rhAmvz9xCBN/uaWmaZ65DYcJdLkf76RJv9zjDT5n2ObR1ipNG3ctGXS9zPZDCZLbEhTFlpulUWa3MpX3miZacpLzp16SJM7ucobKdKUl5w79ZAmd3LlY6SVSlN0tbwkmOH3NrkGGmlyLWPp40Wa0rNyuSTS5HL20sWONKXj5HIppMnl7KWLHWlKx4lS5RCoVJrKGUL3W0Waup+DsiJAmsoia1e7SJNd+SgjGqSpDKp2tYk02ZWPMqJBmsqgSptpCVQqTe1Wz1v/wCPq+lvvUevXXZk2bqvKIU1WpcNoMEiTUZzWNoY0WZsaY4EhTcZQWtsQ0mRtaowFhjQZQ0lDOQhYIU1PPPuSWnPB5YrL83JkkCqlEkCaSsVrTeNIkzWpKC0QpKk0tNY0jDRZk4rSAkGaSkNLwykIWCFNt9xxn7pn/UPMNKVIGEWqJYA0Vcu7W70hTd0iX12/SFN1rLvVE9LULfLV9Ys0VceanqYSKF2awlmkTvDXXX+p6saX3oqwXXfz3a3wTli1Qt101cWTwpWl0h9+9JlgW9J+Ls/rlF139yNN7uYuS+RIUxZabpZFmtzMW5aokaYstNwsizS5mTdfoi5dmqKg2t3T1E2YIkRRSZL3S5csUpddeE4Q1teuu129tnFzq0x8v5RBmrqZwXL7RprK5WtL60iTLZkoLw6kqTy2trSMNNmSifLiQJrKY0vLnQlUKk2dw+l+CZl5+vGTG1qStHrNJeryr6xtzYLJzNmlV9wy6VJCpKn7eSsrAqSpLLJ2tYs02ZWPMqJBmsqgalebSJNd+SgjGqSpDKq0mZYA0hQjJZL0mx8+Ophpki/jPfmMi9T9d12rlixeGJRM2oY0pT3c3CuHNLmXszwRI015qLlVB2lyK195okWa8lBzqw7S5Fa+fIu2cmkSKXnl1TcTOXZz9bwwrug9S2mlacfuUd+OC8bTJNDX26MGZ/Sp3Xv2wcRjAnOGBtTOPfpzvN/jQdZ8aEODfWp03361b2y85iT8Hf5Af6/q1b+zR/aO+TvImo9M8juk/03eVbN/k+fOGqh55u0YfqXSlHQ/kB0YJqKI3sOENNmWnerjQZqqZ96NHpGmblCvtk+kqVre3egNaeoG9Wr7RJqq5U1vkwlUKk02LgQRPyDi3xnFPU31/shweV498s/lef7nmcvz/M8xl+f5n2Muz/M/xzaPsPbSJFK0ft2VrRzJbJg8whX1WD3P5sO3/NiQpvIZ29AD0mRDFsqNAWkql68NrSNNNmSh3BiQpnL50vr0BCqVJhGSUz9+glp90vHW5CX6HUwSFN/TZE1qrAgEabIiDaUHgTSVjrjrHSBNXU9B6QEgTaUj7noHSFPXU1DrACqVpqTlun2gz+p5PmQxeQxIk7+5jY4MafI/z0iT/zlGmvzPMdLkf45tHmGl0iT3NE336ObqeUWShDQVoWd3XaTJ7vyYig5pMkXS3naQJntzYyoypMkUSXvbQZrszU0dIqtUmnwFijT5mlmlkCZ/c8tMUz1yG44SafI/30iT/zlGmvzPsc0jRJoMZAdpMgDR0iaQJksTYzgsZpoMA7WwOaTJwqQYDglpMgzUwuaQJguTUqOQKpem6MIL666/VK1cfoSSy/bOPPVj6rILz3ESPdLkZNpSBY00pcLkfCGkyfkUdhwA0tQRkfMFkCbnU9hxAEhTR0QUKJFApdIU/XLb6PcfrX/gEXX9rfdMWvq7xDEbbxppMo7UmgaRJmtSUWogSFOpeK1oHGmyIg2lBoE0lYrXisaRJivSUNsgKpUmmVG6/65r1ZLFC1VUmuJfKOtaNpAm1zKWPl6kKT0rl0siTS5nL13sSFM6Ti6XQppczl662JGmdJwoVQ6BSqVJROmb114yRZqYaSonubRanADSVJyhCy0gTS5kqViMSFMxfi7URppcyFKxGJGmYvyoXYxApdJ0yx33qXvWPxRchhfONB2y+EB18hkXqQvPPV2tPeuUYqPpUm1mmroEvoJukaYKIFvQBdJkQRJKDgFpKhmwBc0jTRYkoeQQkKaSAdP8tAQqlSaJJLwULxrV1V89X60+6XhnU4U0OZu6joEjTR0ReVEAafIijdMOAmnyP8dIk/85Rpr8z7HNI6xcmmyGkTc2pCkvOfvrIU3258hEhEiTCYp2t4E02Z0fE9EhTSYo2t0G0mR3fnyPrlJp+tp1t6s77/2eeu7B2yZxZclx3w8zd8eHNLmbuyyRI01ZaLlZFmlyM29ZokaastBysyzS5GbefIm6UmmS+5gu+OxpUy7FYyEIXw4n/8aBNPmX06QRIU3+5xlp8j/HSJP/OUaa/M+xzSOsVJpkRin8QtsoFJYct/kQqXdsSFM98o80+Z9npMn/HCNN/ucYafI/xzaPsFJpYqbJ5kOB2JIIIE31OC6QJv/zjDT5n2Okyf8cI03+59jmEVYqTXIZ3pe+fmPrC24FzMZNW4Ilx11eQY+FIGw+xIvFhjQV4+dKbaTJlUzljxNpys/OlZpIkyuZyh8n0pSfHTWLE6hUmiTcpCXHky7ZKz606lpAmqpjXXVPSFPVxLvTH9LUHe5V9oo0VUm7O30hTd3hXmWvSFOVtOkrTqByafIxBUiTj1ltjAlp8je30ZEhTf7nGWnyP8dIk/85Rpr8z7HNI0SaDGQHaTIA0dImkCZLE2M4LKTJMFALm0OaLEyK4ZCQJsNALWwOabIwKTUKqXJpksUgXnn1zUTE8e9vciUPSJMrmcoeJ9KUnZmLNZAmF7OWLWZT0rRfdys/4/p/+/V/4/JanvV7eVaqd+J9T2O7lNuveiaVk3pTt+tt0o6uFy0f9NNsP2ivN9J/a3uP2t8TxtaMSzcWxivP8gjehzE1y0fLyD4pNR6rK/E0yjVaitcJ903pT/rS/TTamxxPo6+wrei+sK+JfdHYJ8em6zXH0d/bq3o0m72jjUw06sgraa/5LJ326PfBc4NF+NzY1twY3ddqS2+MtjmprVibYZ2kMs0xh2MKnhtNN+MK25p4lkDD+OLPrXqtOJsvWjE0xhg+gvqyIdqn7JSxhzEEr6RMSDKML6GtZr1G3bBeZDzNtsJjp7EnzEk0rs5x9ui+Bvp61d598gmqz+MHv/b++gzW4pFWKk3nffkatXTJInXZhedYjCR7aEhTdmau1ECaXMlUsTjrKE1ybjOmz3D26ZPJseC1Uvv0+zF9nhRsj70e1yc5Ula2y+nKPn32Ptasu0/XlddBnSntSh8iFI32g/q6j6BsUF7abfQpdYN4wtdhe2Gb8fhke7PdoI60K3G16kl/jTJyUtrYPiEgcoInp9f79YlYKD7CJRSUhrg0ywTbJ5/2FTvqqA0BCEAgHYH9K49JV5BSpRKoVJrafU9TqSOsoHGkqQLIXeoCaeoS+Iq7LUua5AR8z/5xNaJfyHPwo8/C9+gT9cbr/Wqvfj06Pt541iflo/oEPnjWP3vDZ91OsE2fwre26TYb25rldJm9uj15P6rbbrTXKDPRTnNbU3AqxuxFd/K3+F6dI/2Pp/7DfI/q1Yzlfa/eoec5AsWS5+C9NjV51qWDv6u3tst+XXLq9sYf+1v1grYb5cL2g/f6GJq6XcfTmDyZ6EfiDPqZmGhovQ/GMLEvWq5HZmwi9RqvRSEbschjyv7mxEVjTJP7C3i1xhup25zpmLIvYBvtY6LPoK3YvnAcgzN6VZ/eOTyiVbzVtryQ2JvPzZkkyZ1Mi8hzfF8429QoEx1LY9al2fS0bTZKNdoP+2k9x/gEbTYndNqVD8oktdUcw8SHKxzrxJaARWTCKGwnLBG8lwLBRFs4OtnbnCaMfHKjexsxxdsO68XayhXnRNthXHqSSc0dGlDv7Nwbicr/l/9x6YH+D9KBESJNBpKENBmAaGkTSJOliSkQlsw27NY/u7RYNH72q8FZ/WrjjhE1rOVlj34/EkhNKDZaenR/w03Jae0PywblJiQoKkTSjq0POdeRE8x+fXLXF7xWql+/79MhB9tjr3v1WKSsbJcT/n591t7XrNuv60681tubdYM25EfXlf3yPqiv+5jYp1/L/qCcvA7rNMtIbGGs8fia+1oxhXE322r02ah/4JwZamTvuBrVP6GAyAmf6IEIUCg0wqWxPxSQZplge6gMtma13nFxT5P/+eeeJv9zbPMIK5UmuTzv1I+foFafdLzNTDLHhjRlRuZMBaSpu6kS5QjERgtK47khOyI9O8fGWtKzO9jf2LdTv94dEaJdY5H6uozIT5WPmfomi5n6DFyeB/XMw5A+Iw+2yXt9cj5Dn7TP0NsG9POAPimfoYOT5wG9TfZNbOsNyoRlZ+i6QTnZJs/SVvBelwvaaWyb3E5zW1NgquTQ7b5M3dPU7XHQf3sCSJP/RwfS5H+ObR5hpdIk39F06RW3qPXrrrSZSebYkKbMyJypgDQVT9WO8TG1TcvOO/vH1DtadLZpqdmqxeWdsX36Z1xt03eZyPNW/bNdlxX5aUiPliD93vRDZgtma2GZ1XyerZ8XDmrV0P0PaploSI5ITY9+r2Wn+XookBzZr390udbrYL/ep7c1hKjRhgiRPPOwgwDSZEceyowCaSqTrh1tI0125KGuUVQqTXJP03QPVs+r62Fo77iRpkZuZG5GhOcdESAtNCI474xr6dECJBIk2+X1Vi1C27TsBO+b5eUm/yIPEZzZWkyCZxEdea3lZJZ+P6e3r7FNJEhvk+c5zbKyPyw/R18IH76X5/ijrHuaioybumYJIE1medrYGtJkY1bMxoQ0meVJa9kIVCpN2UJzpzQzTe7kKmukPkvT21psNmux2aRFR3426/eb9PvN+vXboSCJCOnt2/Vz3oeIzAF9feoALTPz9esFzdcH9Par+Vp0DujrVwc0t88LhGdiJkikqIoH0lQF5e72gTR1l38VvSNNVVDubh9IU3f51713pMnAEYA0GYBoaROuSZNc1jYhQk0p0jNCIkKhEAXPWoRkyee0D7k5fr4WmgO0xMjPAj1zM7+nIUKB9MhzRI6kjJRfoJ/l/hrbH0iT7RkqHh/SVJyh7S0gTbZnqHh8SFNxhrSQn0Dl0iT3Na254PJJEa+7/lK1cvkR+UfR5ZpIU5cTUGL3NkiTqE0oPQ3xGW0IUDAzJEI0MUO0I8MiBwv1jM9i/bNIz/gs1uKzSAvO4n79Xr8+MCI9oQCViLnrTSNNXU9B6QEgTaUj7noHSFPXU1B6AEhT6YjpYBoClUrT+gceUV/6+o3q/ruuVUsWLwzC2rhpizr5jIvU1V8939lV9ZAmfz9jVUjTW1p6XtPy8+q+UfXaPv2spSh4Hc4O6ZmitJNCsljBokCGRIJ6tQQNqMVagBb1R8RI9uttskwzjwYBpMn/IwFp8j/HSJP/OUaa/M+xzSOsVJpWr7lEXfDZ06bIkcjU9bfe4+yqekiTzYd4sdhMSNNGLUKBFAUy1JSj4PVYIEdySV2nh8wCtWRIC1FrhkjPDAVSFGzrV3JPEI/sBJCm7Mxcq4E0uZax7PEiTdmZuVYDaXItY37FW6k0yep5SZfihZfssXqeXweXD6NJK02yZPbL+/aql0dH1M+0IL08Kq/1j5aiTvcOySIJy7TwLNOzQsv0jNBSfbnc0ubrQ2TGSP+gQuUeTUhTuXxtaB1psiEL5caANJXL14bWkSYbslDfGCqVJhtnmr523e3qznu/1zoCTli1Qt101cWTjgj5Ut6HH30m2Ja0n5kmfz9AcWn6eSBGo+pnWoYaYqSf940E9xi1e8gsUShBy/pEjBpy1Hjdr+YyO9T1Awhp6noKSg8AaSodcdc7QJq6noLSA0CaSkdMB9MQqFSabLynSUQu+mW78v601R9Va886JcAmUvXaxs0tkRKBWrpkkbrswnNaWJEm/z5jP9Ny9MLeveon+jK6l8f16+E96qdakNrNGsl3/7xPy9D7BmY0nmfIs/4ZGFBynxEPuwkgTXbnx0R0SJMJina3gTTZnR8T0SFNJijSRl4ClUqTBGn76nm33HGf+vGTG1qSJBJ1+VfWtlb3k/gvveKWSaKFNOU9/LpfAhbEewAAGu9JREFUT76Y9fnRPWqDFqQNeubo+b0jwXO7+4yW6pmh9/UPBjIUSpGI0lJ9CR0PdwkgTe7mLm3kSFNaUu6WQ5rczV3ayJGmtKQoVwaByqWpjEGYbFNmko479qhgpilc2S9ptb/oNqTJZAbKa+sn+n4jkaMXZBZpRL/Wz7/Sl9klPUSOjtIzRStmDqmVc2epxaM9gSjJF7Xy8I8A0uRfTuMjQpr8zzHS5H+OkSb/c2zzCJGmSHZklum6m+9W4YIUaaXJ5gTXNbYnh4fVY7uH1eP6srrHdu9WT+nnPXpWKf4Y6u3RYjRTHTM0pFYMzWy8njWkDtSLM/CAAAQgAAEIQAACEICAEKhEmsJ7mZK+i2m6fVWmaLr7rZhpqjIT2ft6UV9O98zePeppfWnd0/KsZ5FGlHwl7OTH4Xrm6Cg9W3TUjEH1Qf36g4OD6v36ebpH2tXzskdNDZsIMNNkUzbKiYWZpnK42tQqM002ZaOcWJhpKocrraYjUIk0JS2eEA0vvthCutDNlYrPMEVb5p4mc5xNtCQr1okYPaMvs3taX24nsrQzYQbpvfo+oxUDg+oYLUbHaElaPjAz13cYIU0msmZ/G0iT/TkqGiHSVJSg/fWRJvtzVDRCpKkoQeoXIVCJNLX7fqYw8G5+T5MInTziy4yHsbF6XpHDq1hdWanux1qKHtkzrB4ZGVZP6deycEP8cZherW6FFiORo2Nm6EvstCwtMHR5HdJULIeu1EaaXMlU/jiRpvzsXKmJNLmSqfxxIk352VGzOIFaS1N4z1ISxuiX8PI9TcUPtDQtjGhJemRkt3pUX2YnovSoFqX4Et+H6gUaRIxEkESOVujXiwwJUlKMSFOazLlfBmlyP4edRoA0dSLk/n6kyf0cdhoB0tSJEPvLJFCJNMUvcYsPKGkZ7zIHbbptVs/LR3TX/nEtR3sCOQpkaWTPlIZkBul4vYLdqhlD6sODM9UhFS/tjTTly61rtZAm1zKWPV6kKTsz12ogTa5lLHu8SFN2ZtQwR6ASaYp/91E8/E73PJkbbjktIU3puMqldXKZXUOS9Mp2CZJ0rJ45On5wqClKM41dZpcuwqmlkKa85NyqhzS5la880SJNeai5VQdpcitfeaJFmvJQo44pApVIkwQrs03yWL/uykmxy/ZXXn2ztcy3qYFV2Q7S1J72y/q7kB7QS37fv2eX+v6e3VMK/oaePQokSf+s0j9ze+36HiSkqcpPUvf6Qpq6x76qnpGmqkh3rx+kqXvsq+oZaaqKNP0kEahMmqTzcJW6aCBnnvoxddmF5zidHaRpcvr+Vc8gPTC8S90/vFM9r1e7iz7+zcxZgRwdr2Vplf4ZsvzLYpEmpz+aqYNHmlKjcrYg0uRs6lIHjjSlRuVsQaTJ2dR5EXil0uQFsYRB1F2axvUCDg/oWSSRpAeGd6uNY/talObpL4/9n4fmqpOGZqmTZs5WcyybSep0TCJNnQj5sR9p8iOP040CafI/x0iT/zlGmvzPsc0jRJoMZKeO0rRlbFyL0k592d1uffndLjWsF3UIH4f3DaiTZ81WJw/NVh/VM0suP5Aml7OXPnakKT0rV0siTa5mLn3cSFN6Vq6WRJpczZwfcSNNBvJYF2l6fd8+9U+7dwSiFL8/SRZwOFlmk4bmBMuB+/JAmnzJ5PTjQJr8zzPS5H+OkSb/c4w0+Z9jm0eINBnIjs/SJN+TdK++7O7eXTuC+5Sij481JekkPZv0Lv0Fsz4+kCYfszp1TEiT/3lGmvzPMdLkf46RJv9zbPMIkSYD2fFRmmQm6e92N2Rpd+TSu0/OmqM+oX9O0pfezbF8EQcDqVVIkwmK9reBNNmfo6IRIk1FCdpfH2myP0dFI0SaihKkfhECSFMRes26vkjTdv09Snft2q7u2rlNPRdZ9U6WAz919lx1ql7Q4YA+u5YEN5C+aZtAmsombEf7SJMdeSgzCqSpTLp2tI002ZGHMqNAmsqkS9udCCBNnQil2O+6NP1If9Hs3VqW7tazSmP6cjx5HNrfr06fPU+dOmuuOnJgRgoKfhZBmvzMa3xUSJP/eUaa/M8x0uR/jpEm/3Ns8wiRJgPZcVWa/krPKH1Hy9IT+nuVwodcdneGlqV/qy/B46G4PK8mBwHS5H+ikSb/c4w0+Z9jpMn/HNs8QqTJQHZckqbX9XcordOydMeObWrz+Fgw+sV9/VqU5qoz5sxT7+uv76xS0qHATJOBD4gDTSBNDiSpYIhIU0GADlRHmhxIUsEQkaaCAKleiADSVAhfo7IL0vSIvgTvjp3b1d/omaXw8eHBmepsPat05pz5Bij42QTS5Gde46NCmvzPM9Lkf46RJv9zjDT5n2ObR4g0GciOzdJ0n14B7zY9s/QDvRpe+Pj3+j6ls/Ws0m85/sWzBlLXsQmkqSMiLwogTV6kcdpBIE3+5xhp8j/HSJP/ObZ5hEiTgezYKE3f1rNKt+3cqjbs3RuM8IDeXrVGzyidrX8O8/Q7lQykckoTSFMZVO1rE2myLyemI0KaTBO1rz2kyb6cmI4IaTJNlPayEECastBqU9YWado6NqZu17NKMrO0Sd+7JI/36nuUfk+L0mfmzlf9PT0GRluvJpCmeuQbafI/z0iT/zlGmvzPMdLkf45tHiHSZCA73ZamV/aNqm/teCeQpX3NJcOPmzmkztH3K52mf3jkJ4A05WfnUk2kyaVs5YsVacrHzaVaSJNL2coXK9KUjxu1zBBAmgxw7JY0yaV339qxNVg2PHx8XC8Z/ntzD1C/zf1KBjLLkuNGIDrQCNLkQJIKhog0FQToQHWkyYEkFQwRaSoIkOqFCCBNhfA1KlctTU/u3aO+uf0ddc/uHa3oP6VnlD6nL8H79RkzDYyIJkICzDTV41hAmvzPM9Lkf46RJv9zjDT5n2ObR4g0GchOVdL0w5Hd6pv6+5W+q1fECx+f1vcrfU7/fGDGoIGR0EScANJUj2MCafI/z0iT/zlGmvzPMdLkf45tHiHSZCA7ZUvTQ3q58Jv1ZXj/73Bj2XBZ0OFz+hK8z885QC3r7zcwAppoRwBpqsexgTT5n2ekyf8cI03+5xhp8j/HNo8QaTKQnbKk6cE9w+rm7VvVv+zZFUQ5t6dXfW6eliUtTAf29hmInCY6EUCaOhHyYz/S5EcepxsF0uR/jpEm/3OMNPmfY5tHiDQZyI5pafqenlGSmaX/r/mFtPN6e9R5cxeqL2hhmq3FiUd1BJCm6lh3syekqZv0q+kbaaqGczd7QZq6Sb+avpGmajjTSzIBpMnAkWFKmh7QknTTtq3q+/reJXks6OtT5+pZpS/onyFkyUCmsjeBNGVn5mINpMnFrGWLGWnKxsvF0kiTi1nLFjPSlI0Xpc0SQJoM8CwqTQ8M71I36svwfjgyHEQjl97JrNK5cxeoQb6Q1kCG8jeBNOVn51JNpMmlbOWLFWnKx82lWkiTS9nKFyvSlI8btcwQQJoMcMwrTf+iZekb+ktpZaEHeSwKZGlBMLs0gCwZyEzxJpCm4gxdaAFpciFLxWJEmorxc6E20uRClorFiDQV40ftYgSQpmL8gtpZpekHWpJkZul7TVmSy/C+qEXpi3pmSVbG42EPAaTJnlyUGQnSVCZdO9pGmuzIQ5lRIE1l0rWjbaTJjjzUNQqkyUDm00rTo3o1vG/oBR7+h55hksdcPbMUyJK+FG8m9ywZyIT5JpAm80xtbBFpsjErZmNCmszytLE1pMnGrJiNCWkyy5PWshFAmrLxSizdSZqe3Duibti+Rd3X/FJaEaQv6svwRJZkGXEe9hJAmuzNjcnIkCaTNO1sC2myMy8mo0KaTNK0sy2kyc681CUqpMlApttJ04bRveqGbVvU3+7eEfQil96dp7+Q9vz5C9WCXmTJAPrSm0CaSkdsRQdIkxVpKDUIpKlUvFY0jjRZkYZSg0CaSsVL4x0IIE1NQOsfeERdf+s9av26K6cgO+/L16iHH30m2H7CqhXqpqsunlQmLk2/3DeqZ5a2qm/v3NYqd66eVTpf37O0uK+fg9IhAkiTQ8kqECrSVACeI1WRJkcSVSBMpKkAPEeqIk2OJMrTMGsvTU88+5Jac8HlQXoPX3bwFGn62nW3q9c2bm6JkgjU0iWL1GUXntM6JEJpemtsLJClm/R9S+HjnLnz1e9rWXpX/4Cnh5Dfw0Ka/M5vODqkyf88I03+5xhp8j/HSJP/ObZ5hLWXpjA57WaaVq+5RF3+lbVq5fIjgqIiWZdeccskuXr5rV3qL7Qs3bD9HTWyfzwo96nZ89Tvz1uojhxAlmz+AHSKDWnqRMiP/UiTH3mcbhRIk/85Rpr8zzHS5H+ObR4h0tTMTpI0bdy0RZ18xkXq/ruuVUsWLwxKxrddvWmz+j82vqm2jjdkafWsOcHM0srBmTbnndhSEkCaUoJyvBjS5HgCU4SPNKWA5HgRpMnxBKYIH2lKAYkipRFAmgpKU88TTwctfGzWbHXhQQeq3x6aXVqyaLh6An29PWpwRp/avWdf9Z3TY2UE5gwNqJ17RpXaX1mXdFQxgaHBPjW6b7/aN9b4AxcP/wgM9PeqXv07e2TvmH+DY0QBAcnvkP43eVfN/k2eO4urlmz4CCBNBaXpPc9vUFccdLBaPWeODfkkBsMEkCbDQC1tDmmyNDEGw0KaDMK0tCmkydLEGAwLaTIIk6YyE0CappEm2ZXmnqZO39OUOStUsIYAl+dZk4pSA+HyvFLxWtE4l+dZkYZSg+DyvFLxWtE4l+dZkYbaBoE0dZCmLKvn1fYo8njgSJPHyY0MDWnyP89Ik/85Rpr8zzHS5H+ObR5h7aUpuuR4mKgLzz1drT3rlFbesn5Pk80JJ7ZsBJCmbLxcLY00uZq59HEjTelZuVoSaXI1c+njRprSs6KkeQK1lyYTSLk8zwRFO9tAmuzMi+mokCbTRO1rD2myLyemI0KaTBO1rz2kyb6c1CkipMlAtpEmAxAtbQJpsjQxhsNCmgwDtbA5pMnCpBgOCWkyDNTC5pAmC5NSo5CQJgPJRpoMQLS0CaTJ0sQYDgtpMgzUwuaQJguTYjgkpMkwUAubQ5osTEqNQkKaDCQbaTIA0dImkCZLE2M4LKTJMFALm0OaLEyK4ZCQJsNALWwOabIwKTUKCWkykGykyQBES5tAmixNjOGwkCbDQC1sDmmyMCmGQ0KaDAO1sDmkycKk1CgkpKlGyWaoEIAABCAAAQhAAAIQgEB2AkhTdmbUgAAEIAABCEAAAhCAAARqRABpqlGyGSoEIAABCEAAAhCAAAQgkJ0A0pSdGTUgAAEIQAACEIAABCAAgRoRQJpyJvu8L1+jHn70maD2CatWqJuuujhnS1TrNoEsufzadberO+/9Xitkct/t7KXvP0ueo62uf+AR9aWv36jWXX+pWrn8iPQdUrJyAnlyfPSJn2nFefVXz1erTzq+8rjpMD2BrDleveYS9cqrb7Y6eO7B29J3RknrCMjv4+tvvUetX3eldbERkP8EkKYcOZYT59c2bm6JkvwSX7pkkbrswnNytEaVbhLImkv5Bzj6y1ren7b6o2rtWad0cxj03YFA1jyHzYX/QMtJF9Jk92GWNccbN21RJ59xkUKU7M5rNLqsOY7/2xyv787IifSJZ19Say64PABx+LKDkSYOia4QQJpyYJcT5cu/srb1V2f5MF96xS18iHOw7HaVorm85Y771I+f3MBMY7cT2aH/PHmO/kVTZiOQJruTnDXHcgJ96MEH8gcPu9M6KbqsOZbyF3z2tNbsIbMUDiW7Tajk0P0cujwCpClj9sK/Tt5/17VqyeKFQe2kbRmbpXgXCJjIpfwl87hjj+LEqwv5S9tlnjzH/2FGmtLS7k65PDmWnMpfrKOXbkV/r3dnJPTajkCeHMsfta67+W514bmnB7+j+X3t/vGFNLmfQ5dHgDRlzF6eX9wZu6B4RQSK5jL8B5lr5CtKWM5usuY56R9lpCkn/IqqZc1xWD46e8jnuaJk5ewma46jf9CMyjG/r3MmwJJqSJMliahpGEhTxsTn+cWdsQuKV0SgSC7DxQH4y3RFySrQTdY8xxf7iHbN/S8FElFi1aw5bnd1AHJcYpIKNp01x9JdPJ+IccEkWFAdabIgCTUOAWnKkfys11Xn6IIqFRHIk0v+4a0oOQa7yZPnaPecTBtMRklNZc1xUk7Jc0nJMdRslhznkSxDYdJMiQSQphLh0nRHAkhTR0RTC2RdwSdHF1SpiECnXMo18PIIl5SPv68oTLopSCBrnuPdcTJdMAEVVM+aYyn/o8eeay3gI38MuWf9QyzoU0Gu8naRNcfyuT3z1I+1VrYlx3nJ21MPabInF3WMBGnKmfWs3xWRsxuqVUBgulxGJSn8y2VSSKysVkGiCnaRNs9J3SBNBeFXVD1rjqOXYrKMcUVJKthN1hxHv4eLHBeE38Xq0SXHwzDCBT66GBZd14wA0lSzhDNcCEAAAhCAAAQgAAEIQCAbAaQpGy9KQwACEIAABCAAAQhAAAI1I4A01SzhDBcCEIAABCAAAQhAAAIQyEYAacrGi9IQgAAEIAABCEAAAhCAQM0IIE01SzjDhQAEIAABCEAAAhCAAASyEUCasvGiNAQgAAEIQAACEIAABCBQMwJIU80SznAhAAEIQAACEIAABCAAgWwEkKZsvCgNAQhAAAIQgAAEIAABCNSMANJUs4QzXAhAAAIQgAAEIAABCEAgGwGkKRsvSkMAAhCAAAQgAAEIQAACNSOANNUs4QwXAhCAAAQgAAEIQAACEMhGAGnKxovSEIAABCAAAQhAAAIQgEDNCCBNNUs4w4UABCAAAQhAAAIQgAAEshFAmrLxojQEIAABCEAAAhCAAAQgUDMCSFPNEs5wIQABCEAAAhCAAAQgAIFsBJCmbLwoDQEIQAACEIAABCAAAQjUjADSVLOEM1wIQAACEIAABCAAAQhAIBsBpCkbL0pDAAIQgAAEIAABCEAAAjUjgDTVLOEMFwIQgAAEIAABCEAAAhDIRgBpysaL0hCAAASsILD+gUfUl75+45RYrv7q+Wr1Scer1WsuCfatX3fllDKy712HLlY3XXVxo0yHto4+8TPTjvnwZQcH/XztutvVnfd+L7HsuusvVSuXH6HO+/I16uFHn1Hh+7DwE8++pNZccLk6YdWKVlzxhtLEcdrqj6rrbr67VfXCc09Xa886JVO/acZhxUFAEBCAAAQgUBkBpKky1HQEAQhAwAyB8KT+/ruuVUsWL2w1KvJz7/94uCUdIhlnnvoxddmF57TK3HLHfeqe9Q+1ZCptW3G5iUuP7Je2Xtu4ua30SJlQmuJxhdunk6YovVCykuJI2pel3zTjMJNJWoEABCAAAVcIIE2uZIo4IQABCDQJiAyFMyjTQYnLw8ZNW9TJZ1w0aZYnbVsmpWnpkkXBjFQofWFcIlKdpCtNHO2kKW2/SBMfNQhAAAIQiBNAmjgmIAABCDhGIH553XThiwD86LHngpklmW0RcYjOPGVpS/qZboYnjWxIDMcde5R6/c231aEHHxhcOiezX/KQbWVKU9p+04zDsUOGcCEAAQhAoCABpKkgQKpDAAIQqJpAKC7RfpMuUwv3R+8Feu7B2yaFm7WtTtKU5p4mkZcPH/OB4B4miUfik1mnv/zOP5UuTWn65Z6mqo9o+oMABCBgPwGkyf4cESEEIACBtgTiJ/hJl+2FohMuEtGusSxtFbmnSaQpXJxBYglnv7LM8OS5pyltv1ni4NCEAAQgAIF6EECa6pFnRgkBCNSAgFzmJivHxWeTku5l6oSjXVudZpo6XV4XXp4n0hSu2hcKWBZZKSJNnfrNEkcnjuyHAAQgAAE/CCBNfuSRUUAAAjUhIAJ03/0/DGZq4o9QBuKr6rWTpjxtmZQmiV/uqQqXRc8iK0WkqVO/WeKoyWHHMCEAAQjUngDSVPtDAAAQgIBLBEIBkpjjM0pJS4xLuemkSVbTy9KWaWmKss8iK0Wlabp+s8Th0rFDrBCAAAQgkJ8A0pSfHTUhAAEIdI1A0he9trtnqdPleVna6iRNaReCSJopyyIr7eIILysMExP9ctvwnqZ40uL9shBE1w5rOoYABCBgLQGkydrUEBgEIAABCEAAAhCAAAQgYAMBpMmGLBADBCAAAQhAAAIQgAAEIGAtAaTJ2tQQGAQgAAEIQAACEIAABCBgAwGkyYYsEAMEIAABCEAAAhCAAAQgYC0BpMna1BAYBCAAAQhAAAIQgAAEIGADAaTJhiwQAwQgAAEIQAACEIAABCBgLQGkydrUEBgEIAABCEAAAhCAAAQgYAMBpMmGLBADBCAAAQhAAAIQgAAEIGAtAaTJ2tQQGAQgAAEIQAACEIAABCBgAwGkyYYsEAMEIAABCEAAAhCAAAQgYC0BpMna1BAYBCAAAQhAAAIQgAAEIGADAaTJhiwQAwQgAAEIQAACEIAABCBgLQGkydrUEBgEIAABCEAAAhCAAAQgYAMBpMmGLBADBCAAAQhAAAIQgAAEIGAtAaTJ2tQQGAQgAAEIQAACEIAABCBgAwGkyYYsEAMEIAABCEAAAhCAAAQgYC0BpMna1BAYBCAAAQhAAAIQgAAEIGADAaTJhiwQAwQgAAEIQAACEIAABCBgLQGkydrUEBgEIAABCEAAAhCAAAQgYAMBpMmGLBADBCAAAQhAAAIQgAAEIGAtAaTJ2tQQGAQgAAEIQAACEIAABCBgAwGkyYYsEAMEIAABCEAAAhCAAAQgYC0BpMna1BAYBCAAAQhAAAIQgAAEIGADAaTJhiwQAwQgAAEIQAACEIAABCBgLQGkydrUEBgEIAABCEAAAhCAAAQgYAMBpMmGLBADBCAAAQhAAAIQgAAEIGAtAaTJ2tQQGAQgAAEIQAACEIAABCBgAwGkyYYsEAMEIAABCEAAAhCAAAQgYC0BpMna1BAYBCAAAQhAAAIQgAAEIGADAaTJhiwQAwQgAAEIQAACEIAABCBgLQGkydrUEBgEIAABCEAAAhCAAAQgYAMBpMmGLBADBCAAAQhAAAIQgAAEIGAtAaTJ2tQQGAQgAAEIQAACEIAABCBgAwGkyYYsEAMEIAABCEAAAhCAAAQgYC0BpMna1BAYBCAAAQhAAAIQgAAEIGADAaTJhiwQAwQgAAEIQAACEIAABCBgLQGkydrUEBgEIAABCEAAAhCAAAQgYAMBpMmGLBADBCAAAQhAAAIQgAAEIGAtgf8fPlrSETMQh3gAAAAASUVORK5CYII=", "text/html": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Same plot, but with a smoothed line\n", "bio.plot_history_single_bin(bin_address=0, \n", " title=\"Reaction A <-> B . Concentrations at bin 0\",\n", " smoothed=True)" ] }, { "cell_type": "markdown", "id": "37879680-50e8-4564-a872-a1f94cedfd22", "metadata": {}, "source": [ "## For more in-depth analysis of this reaction, see the experiment `reactions_single_compartment/react_1`" ] }, { "cell_type": "code", "execution_count": null, "id": "da56d751", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "jupytext": { "formats": "ipynb,py:percent" }, "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.13" } }, "nbformat": 4, "nbformat_minor": 5 }