{ "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 = \"June 6, 2025\"\n", "LIFE123_VERSION = \"1.0.0rc6\" # 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, os\n", "#os.getcwd()\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, check_version" ] }, { "cell_type": "code", "execution_count": 4, "id": "a14ca0a0-a287-4bba-a20f-ccd09b5e82de", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "OK\n" ] } ], "source": [ "check_version(LIFE123_VERSION)" ] }, { "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": 5, "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": 5, "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(chem_label=\"A\", conc=10.) # Same across all bins\n", "bio.set_uniform_concentration(chem_label=\"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": 6, "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": 7, "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", " \n", " \n", "
SYSTEM TIMEABcaption
00.010.050.0
\n", "
" ], "text/plain": [ " SYSTEM TIME A B caption\n", "0 0.0 10.0 50.0 " ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "bio.get_bin_history(bin_address=0)" ] }, { "cell_type": "code", "execution_count": 8, "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", " \n", " \n", "
SYSTEM TIMEABcaption
00.010.050.0
\n", "
" ], "text/plain": [ " SYSTEM TIME A B caption\n", "0 0.0 10.0 50.0 " ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "bio.get_bin_history(bin_address=1)" ] }, { "cell_type": "code", "execution_count": 9, "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", " \n", " \n", "
SYSTEM TIMEABcaption
00.010.050.0
\n", "
" ], "text/plain": [ " SYSTEM TIME A B caption\n", "0 0.0 10.0 50.0 " ] }, "execution_count": 9, "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": 10, "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": 11, "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": 11, "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": 12, "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", " \n", " \n", " \n", "
SYSTEM TIMEABcaption
00.010.050.0
10.117.043.0
\n", "
" ], "text/plain": [ " SYSTEM TIME A B caption\n", "0 0.0 10.0 50.0 \n", "1 0.1 17.0 43.0 " ] }, "execution_count": 12, "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": 13, "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": 13, "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": 14, "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", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SYSTEM TIMEABcaption
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 caption\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=0)" ] }, { "cell_type": "code", "execution_count": 15, "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", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SYSTEM TIMEABcaption
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 caption\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": 15, "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" ] }, { "cell_type": "code", "execution_count": 16, "id": "490fdc0f-fa2a-48d8-ae04-7c80d298842a", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "{'A': 23.9931640625, 'B': 36.0068359375}" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ "bio.bin_snapshot(bin_address=0)" ] }, { "cell_type": "code", "execution_count": 17, "id": "f22258b2-5181-4ff9-b379-f6a12ad5c8fb", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0: A <-> B\n", "Current 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": 17, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Verify that the reaction has reached equilibrium\n", "bio.get_reaction_handler().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": 18, "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": "iVBORw0KGgoAAAANSUhEUgAAA4EAAAFoCAYAAADpfv6vAAAgAElEQVR4Xu29C7AcxX3o3eelc/Q4Qg9AL67KhGAFg3gIGcd1lcIprh1TuPhk7v1K6DoiOMYR2BCgKEwqpkqlr3DKmHIB1zggX7AJGAPXN4Q4TvCHL19QLgmGiIcB2Vd+YWP0QkLofXR0pKOv/7Pbq9k5s7szuz0z3bO/pQ7nnN2e7n///r1H89vu6ek5ph+KBwQgAAEIQAACEIAABCAAAQh0BYEeJLAr8kwnIQABCEAAAhCAAAQgAAEIBASQQAYCBCAAAQhAAAIQgAAEIACBLiKABHZRsukqBCAAAQhAAAIQgAAEIAABJJAxAAEIQAACEIAABCAAAQhAoIsIIIFdlGy6CgEIQAACEIAABCAAAQhAAAlkDEAAAhCAAAQgAAEIQAACEOgiAkhgFyWbrkIAAhCAAAQgAAEIQAACEEACGQMQgAAEIAABCEAAAhCAAAS6iAAS2EXJpqsQgAAEIAABCEAAAhCAAASQQMYABCAAAQhAAAIQgAAEIACBLiKABHZRsukqBCAAAQhAAAIQgAAEIAABJJAxAAEIQAACEIAABCAAAQhAoIsIIIFdlGy6CgEIQAACEIAABCAAAQhAAAlkDEAAAhCAAAQgAAEIQAACEOgiAkhgFyWbrkIAAhCAAAQgAAEIQAACEEACGQMQgAAEIAABCEAAAhCAAAS6iAAS2EXJpqsQgAAEIAABCEAAAhCAAASQQMYABCAAAQhAAAIQgAAEIACBLiKABHZRsukqBCAAAQhAAAIQgAAEIAABJJAxAAEIQAACEIAABCAAAQhAoIsIIIFdlGy6CgEIQAACEIAABCAAAQhAAAlkDEAAAhCAAAQgAAEIQAACEOgiAkhgFyWbrkIAAhCAAAQgAAEIQAACEEACGQMQgAAEIAABCEAAAhCAAAS6iAAS2EXJpqsQgAAEIAABCEAAAhCAAASQQMYABCAAAQhAAAIQgAAEIACBLiKABHZRsukqBCAAAQhAAAIQgAAEIAABJJAxAAEIQAACEIAABCAAAQhAoIsIIIFdlGy6CgEIQAACEIAABCAAAQhAAAlkDEAAAhCAAAQgAAEIQAACEOgiAkhgFyWbrkIAAhCAAAQgAAEIQAACEEACGQMQgAAEIAABCEAAAhCAAAS6iAAS2EXJpqsQgAAEIAABCEAAAhCAAASQQMYABCAAAQhAAAIQgAAEIACBLiKABHZRsukqBCAAAQhAAAIQgAAEIAABJJAxAAEIQAACEIAABCAAAQhAoIsIIIFdlGy6CgEIdB+BMz9ypfrYhUvVnWuv7b7O02MIQAACEIAABGIJeCuBN665Rz29fsOETp046wS1/om7nUh3o5OvR598Rt1218PqyhUfVzdfc7kTsUaDWLF6rXpj05vKJZ42QZkcxNWZ1wmzjI/oI6+2bbJMWleWMvLci6+r1V/4qrr1hlVq5fKL6kIyfyvWfeUmteyCxUnDdbKcvC+37diV6m9cUu5ZcjL5CUN96pHb1cIFc5zkTFAQgAAEIACBshPwXgKjJxLmxNqFEwxfJfCtzdvVxZ+6pTb2bZ88G8FMW685zkZujQRGpcGcrGYpY+ZkO9pGmPvGZx8s3d+epDLSTseRwMbUknLPSgLj3mt33PuYevDxHygb7+V2xgvHQAACEIAABLqdQOkk0JwMujDLlvTky7VBaE7QRNJkdsWWEIVnb9s9+TMiKMzarUOObSSB8lo7sy1Jc9isXVOHtP/4ujVJq/SmXJbvh2YS6A2gBIG2Mzaz5J4gZHXhZderuSfNmjCm5fkli09nmWoSiJSBAAQgAAEIWCbQVRIYXUJ61qJTY09Mdu7aU4e50ZLIuCWFRpjilvqZehqdsCZZMhU+9sVXflq3JDZuKVw74yV80makq5OZKRvyF+1HpzLYTMak//LIYllxO3UnGRfhmZVV1/2VCo/hOFmOzvZKf6PjPG58R8eBkZKHv/aXdbPH0bqSvh/e3rojmCGSh7yXLjjvjGDpdPQRHutxfMzxch1co1mnLLhKfsPsky6njh4Xl4/wmA/zaPVhiJHAKMvoB2VxnJLmt9HfGTPO4j6UM38XOvnb0s7fN46BAAQgAAEIQECp0klgo5N7c5IVPuGIOyGX5770F1fVXTsUV86cwIRPwOSER07AjTw0+gQ+TgLDs2/muqW458InruETYVvLq6KxJZm5avRGykL+mslgmpPJRv3qpL+t/qCYE+I0M6tJx4UpJzGEx2TczFGz94jInFynFbc0MG7MN7p2NG7st3o/SOxRWZBYn3zquboPa+LibzYTGPfeyIKr9C/8wVKafCf9u9PuTGBYiuXnuBUTjSQw7trgpLOLafPS6j3E6xCAAAQgAAEI2CFQKglsdNLVSJBM+VYzaOak05xcJ1161uqkN9xuo7LRk75GbTf7xD3NUIn7dD56ctuqvvAMUqtZilZ1JX097UY2zTaGiYpU0hhalWtnqXLScdFojJt+hq+/bHUC32x8R5fwNZKSRuMoToCTvp/CfKNxpJUN21xbtd/OBlDRvzvS/3YlMI57NEfNZgKjs+JJZ/Gafahi64OrVu87XocABCAAAQhAYCIB7yUw2qU46Wh24iQng42WRUXrNifSSU9ekkpgqxPI8AYKWUtgXMxxM0DN3kztSGB4JsvUHbdct1G77Upg3AcAJpZWHw6k/YOSVgLTjIukEphEuJqNb+EsD3PNYl4SGLecNDw+0rBKUzYp1/Dy2rQbHplxFPcekNfC9dmUwOgHBEhg2nc05SEAAQhAAAL+EvBeAqOzc3HX4MRdbxNOmfmUPHwiF5bJ6MlS0k/Bk0pg3GxN9MTQnAhmKYGtZsfSbrbj43JQw72da/da/RlIu/NomnGRVFaa1Wnib3T7FfN6+D2WtQQ22kk1utlIGrHLgquwibvOUp5vtUw56d8dqcs3CUyTl1bvH16HAAQgAAEIQMAegdJIoCAxJ3fRGaTo7EUjfElPpMs6E9jsBLMTKcpCBrPcGMacbMu1UK1O4NO+FdNwTHMCnXTsppkJTNL3rCWw0QcpnUhgFlzjxkGjv0fRsklzZ1sCo+1mMRPIxjBp/0JQHgIQgAAEIJAPgVJJoCAzwhGetUoqbY1m+KIzB61OIs31P422QI87Puk1SlnNBLbaxCJuI420Q9SGDHYqfybmVhvApJG1NBwMx2azqsJJdrWUR9JxkUYkmt2/8j9+8Cw9o/VOwxuvR/uaRgLTvB+knWYCEZXAZtf3xrGxzVXG03+Yf/KEm9EnmblL+nfH/H17+fWfp9q5tpO+pslvo/cBt4hI8xeCshCAAAQgAIF8CJROAs2n5TKLE76mK253UHNSJVunr1x+Uey948JLJKPX5kgb4WWjcjIXPkFrdHLXbHfQuB0/w+1mJYGtJK+VJKYZrkbk0l47ZY6zsdlMq/sESm6TxGfqSbPjZ5qbxcddn9hsZ8som7ilj3F1RjchaXTjcCkntyYxkppGEtK8H8x4arSDr9yGITrj30g2mu0O2ur9llSuDb+w3Cd9z8SNxUZ/d5Is542+F+M2dko6htLkt9HfgLj+Jf1gLs3fFcpCAAIQgAAEIJCcQCklULpvpC98ct7oWqfwsrfodXHhe5VFpSBuI4foErrw9Yi27xMo4moene4OmmTmy8Y9A5MPzWxLtrr+MclSSImwHQmU4xpdPxY3Qxi+LYihEpW9pLJijo+rMzq+GzFKslFJI+FL834Ii2D0vodyK5boDcijTM17vxEb21zj/h4kvY42zd+d6N+xVh+KiATK357o/U+TjCEbEih5TMI623c8tUMAAhCAAAQgECbgrQSSRghAAAIQgAAEIAABCEAAAhBITwAJTM+MIyAAAQhAAAIQgAAEIAABCHhLAAn0NnUEDgEIQAACEIAABCAAAQhAID0BJDA9M46AAAQgAAEIQAACEIAABCDgLQEk0NvUETgEIAABCEAAAhCAAAQgAIH0BJDA9Mw4AgIQgAAEIAABCEAAAhCAgLcEkEBvU0fgEIAABCAAAQhAAAIQgAAE0hNAAtMz4wgIQAACEIAABCAAAQhAAALeEkACvU0dgUMAAhCAAAQgAAEIQAACEEhPAAlMz4wjIAABCEAAAhCAAAQgAAEIeEsACfQ2dQQOAQhAAAIQgAAEIAABCEAgPQEkMD0zjoAABCAAAQhAAAIQgAAEIOAtASTQ29QROAQgAAEIQAACEIAABCAAgfQEkMD0zDgCAhCAAAQgAAEIQAACEICAtwSQQG9TR+AQgAAEIAABCEAAAhCAAATSE0AC0zPjCAhAAAIQgAAEIAABCEAAAt4SQAK9TR2BQwACEIAABCAAAQhAAAIQSE8ACUzPjCMgAAEIQAACEIAABCAAAQh4SwAJ9DZ1BA4BCEAAAhCAAAQgAAEIQCA9ASQwPTOOgAAEIAABCEAAAhCAAAQg4C0BJNDb1BE4BCAAAQhAAAIQgAAEIACB9ASQwPTMOAICEIAABCAAAQhAAAIQgIC3BJBAb1NH4BCAAAQgAAEIQAACEIAABNITQALTM+MICEAAAhCAAAQgAAEIQAAC3hJAAr1NHYFDAAIQgAAEIAABCEAAAhBITwAJTM+MIyAAAQhAAAIQgAAEIAABCHhLAAn0NnUEDgEIQAACEIAABCAAAQhAID0BJDA9M46AAAQgAAEIQAACEIAABCDgLQEk0NvUETgEIAABCEAAAhCAAAQgAIH0BJDA9Mw4AgIQgAAEIAABCEAAAhCAgLcEkEBvU0fgEIAABCAAAQhAAAIQgAAE0hNAAtMz4wgIQAACEIAABCAAAQhAAALeEkACvU0dgUMAAhCAAAQgAAEIQAACEEhPAAlMz4wjIAABCEAAAhCAAAQgAAEIeEsACfQ2dQQOAQhAAAIQgAAEIAABCEAgPQEkMD0zjoAABCAAAQhAAAIQgAAEIOAtASTQ29QROAQgAAEIQAACEIAABCAAgfQEkMD0zDgCAhCAAAQgAAEIQAACEICAtwSQQG9TR+AQgAAEIAABCEAAAhCAAATSE0AC0zPjCAhAAAIQgAAEIAABCEAAAt4SQAK9TR2BQwACEIAABCAAAQhAAAIQSE8ACUzPjCMgAAEIQAACEIAABCAAAQh4SwAJ9DZ1BA4BCEAAAhCAAAQgAAEIQCA9ASQwPTOOgAAEIAABCEAAAhCAAAQg4C0BJNDb1BE4BCAAAQhAAAIQgAAEIACB9ASQwPTMOAICEIAABCAAAQhAAAIQgIC3BJBAb1NH4BCAAAQgAAEIQAACEIAABNITQALTM+MICEAAAhCAAAQgAAEIQAAC3hJAAr1NHYFDAAIQgAAEIAABCEAAAhBITwAJTM+MIyAAAQhAAAIQgAAEIAABCHhLAAm0kLot745YqIUqfCYwdahf9ff1qD0HxnzuBrFbIDA0qU9NGexTu/YdtlAbVfhMYFJ/r5o+dUDt3DPqczeI3QKBvt4edeIJg2r7e4cs1EYVPhPo6VFq7szJausuzh3nz57scyq9jx0JtJBCJNACRM+rQAI9T6DF8JFAizA9rwoJ9DyBFsNHAi3C9LwqJPB4ApHAYgczEmiBPxJoAaLnVSCBnifQYvhIoEWYnleFBHqeQIvhI4EWYXpeFRKIBLoyhJFAC5lAAi1A9LwKJNDzBFoMHwm0CNPzqpBAzxNoMXwk0CJMz6tCApFAV4YwEqgzsWL1WvXGpjfrcnLlio+rm6+5PHjuxjX3qKfXbwh+PmvRqerxdWvqyiKBrgzn4uJAAotj71rLSKBrGSkuHiSwOPautYwEupaR4uJBApHA4kZffctIYFUCl567qCZ9YUSPPvmMuu+h76n1T9wdPC3CGC2LBLoynIuLAwksjr1rLSOBrmWkuHiQwOLYu9YyEuhaRoqLBwlEAosbfUjgBPZxYmcKRV+LSqGUQwJdGc7FxYEEFsfetZaRQNcyUlw8SGBx7F1rGQl0LSPFxYMEIoHFjT4kMFYCw8tBw0tBL7zsenX1FZeqlcsvCo577sXX1eovfFVtfPbBWj1IoCvDubg4kMDi2LvWMhLoWkaKiwcJLI69ay0jga5lpLh4kMBiJPCOex9T3//h87WVfcWNgM5aNh7y1CO3q4UL5nRUGctBI/gM3HVfuUktu2CxOvMjV6pbb1g1QQIN/D/69h+ps09aolYv+ZyaM3VuR8ngYH8JDOj7gfXq+0CNHj7qbyeI3AqB/r5eNdDfo0ZGGQtWgHZQyfj4seB9WdRDTvwH9X0jDx46UlQItOsIgR595j9lqE8dGGEsOJKS4sLQf5KmDQ2o/SPcV3h4yoC1PLy1ebu6+FO3TKjPnK8XKYEyobRk8enqzrXXdtxfJLBjhM0rCC8BbTUT2LO2coIx0Dugli/6z+r6C25SZ885N+MIqd41AvKJv5xrHhobdy004smZQH9fjxrQIjjCBwI5k5/YnBMSONCrDvKBQOFjoegA5N+HKYP9aj8fCBSdisLbl7PGaZP71T4+EFDDmoONhwjeg4//oG7SRuoVYfraA08EGzoWKYE2+mjqQAJt0oypKyyBra4J3LBlg/rSs3eof/zl36mx8cqnOkvmfFBddc616pLf/aTq77UzwDPuMtV3SIDloB0CLNHhLActUTI77ArLQTsEWKLDWQ5aomR22BWWgx4HaONm8WYGMLxqLy5FRgLlEq/b7nq4ViR8eZc8aYTSFIhbDSjPheuQ1YO/3fJO3XPh5ZriEvPnzq6bCYzemcC0I3uPhOuWOMKXqSGBHb4Bw4fL4LnzG9+tJcbAN8lLujvoOwe2qW+9vk49svGb6t2RHUET86bOV1cs/jO16qyr1MyhWRajpirXCCCBrmWkuHiQwOLYu9YyEuhaRoqLBwksjr1rLSOBdiXQSFtU5qJ5N+U+duHS2jm/iJg8zK3forOFRjDNJWJGwMK3izO3kQvXK89t2fZurd6oBMoqw7knzaq9HnYR8Y63t+6o3bEgepkaEmj5HS3X/YUfJtnmuTT3CRw9Oqqe/Nnj6v4f/7X6yc7XgiqG+ofUZYtWqqvPu16dNuP9lqOnOhcIIIEuZMGNGJBAN/LgQhRIoAtZcCMGJNCNPLgQBRJoVwKjwtUox3HLQaMTPeID0Q1XpH55yPV8cQJmJo/CEho3gWRmAuM2mGw1LkUaP/HRDwdiiAS2opXz6412B31+8/9WD/z46+rpN/9RHT1W2SRi2Sl/qK469/PqP73vYtWj/+NRDgJIYDnyaKMXSKANiuWoAwksRx5t9AIJtEGxHHUggW5JoCy9FIEzchU3yszMXxoJNPVKfeGZwCTXJkaXpEodZkkoEujY34FWt4jYsv9tdf+rX1eP/fRv1J7R3UH07zvhd9Snz75G/dcPXKmmDEx1rEeEk5YAEpiWWHnLI4HlzW3aniGBaYmVtzwSWN7cpu0ZEmhXApNIlbTYaCYwKoHNbr2QhwSaawXDM4vMBKZ9l+VYvpUEmlBGjoyo7/6fb+vZwb9Wv3hvU/D08KTpasUZq9Rnz71OnTK8MMeoacomASTQJk2/60IC/c6fzeiRQJs0/a4LCfQ7fzajRwLtSqARs0Ybw4j8yTLKVhIoUUVvCxfNuw0JjFs+Gm4nLgYk0OY70HJdSSUw3Oyzb/0wkMF//s3T6pj+r7enN1gi+plzPq+XjH7EcoRUlzUBJDBrwv7UjwT6k6usI0UCsybsT/1IoD+5yjpSJNCuBEptZu+OqAiKUJmlnEkk0NQTng0U8fvbf1yf+prARstBJd7ofQPDG8OEhU/KGmlkOWjW78w2629HAk1Tv9nzpvrGq/9NzxA+og6M7Q+ePmP2WVoGPxdsJjPYN9hmVByWJwEkME/abreFBLqdnzyjQwLzpO12W0ig2/nJMzok0L4ESo1x1/SFd+xMIoFST9z1eNHdQcOS2GhjmGYSaERw5649NRhGYKM3vZc+vPz6z9kYJs83aZq2OpFA087+sX3qOxsfVN967V711t5fB0/Pmnyi+uMz/1R9evHV6uSpc9OERNmcCSCBOQN3uDkk0OHk5BwaEpgzcIebQwIdTk7OoSGB2UhgzmksRXM9x/SjFD0psBM2JNCEL0tDZTfRB/RGMv+6eX3w9EDvgLrktE+qzy25UZ150jkF9pSmGxFAAhkbhgASyFgwBJBAxoIhgAQyFgwBJBAJdOXdgARayIRNCQyH87NdP9VLRb+m/u5nj6lDRw4FLy2Z80F11TnXqkt+95Oqv7ffQvRUYYMAEmiDYjnqQALLkUcbvUACbVAsRx1IYDnyaKMXSCASaGMc2agDCbRAMSsJNKG9d2iX+vbGB9TfvLZObT2wJXh63tT56orFf6ZWnXWVmjk0y0IvqKITAkhgJ/TKdSwSWK58dtIbJLATeuU6FgksVz476Q0SiAR2Mn5sHosEWqCZtQSaEI+MH1H/9Msn1f36BvQvbXsheHqof7LeQOZydfV516vTZrzfQm+ooh0CSGA71Mp5DBJYzry20ysksB1q5TwGCSxnXtvpFRKIBLYzbrI4Bgm0QDUvCQyHunHHj9V9r9yt/uEXf6vGxseCl5ad8ofqqnM/H9xqokf/xyM/AkhgfqxdbwkJdD1D+cWHBObH2vWWkEDXM5RffEggEpjfaGveEhJoIRNFSKAJe+fBd9SDb3xDPfzG/Up+lsf7Tvgd9adnf06t/MCfqCkDUy30kCpaEUACWxHqnteRwO7JdaueIoGtCHXP60hg9+S6VU+RQCSw1RjJ63Uk0ALpIiXQhH/46GH19z//H3qp6F+rN3a8Gjw9PGm6WnHGKvXZc69TpwwvtNBTqmhEAAlkbBgCSCBjwRBAAhkLhgASyFgwBJBAJNCVdwMSaCETLkhguBsvbPlX9YC+bvAHv/oHdfTYUdXb0xssEf3MOZ/XS0Y/YqHHVBElgAQyJpBAxkCUABLImEACGQNRAkggEujKuwIJtJAJ1yTQdGnbga2BDH5n47fU7tH3gqfPmH2WlsHP6c1kVqrBvkELvacKIYAEMg6QQMYAEsgYaESAmUDGhiGABCKBrrwbkEALmXBVAk3X5B6D3930iPqmXioq9x6Ux6zJJ6o/PvNP1acXX61OnjrXAoXurgIJ7O78h3vPclDGgiHATCBjwRBAAhkLSODEMTB/9uTSD4xHn3xG3XbXw+rKFR9XN19zuVP9RQItpMN1CQx38V9++0xw3eD/9+sfqGP6v4HeAXXJaZ9Un1tyozrzpHMs0OjOKpDA7sx7XK+RQMYCEsgYiBJAAhkTSGB3SuCK1WvV0nMXqe//8Hm1/om7nXojIIEW0uGTBJru/mbPm+qB176uHv/Jw2r/2L7g6SVzL1BXnf15dcnvflL19/ZbINM9VSCB3ZPrVj1FAlsR6p7XmQnsnly36ikS2IpQ97zOctDjue6GmcAzP3Kl2vjsg0q+r/vKTWrZBYudGexIoIVU+CiBptsigI/95CH1rdfuVb/e86vg6XlT56s/OXu1Xi76GTVzaJYFQuWvAgksf46T9hAJTEqq/OWQwPLnOGkPkcCkpMpfDgnMRgJ/PnpYbR6r3Dc7z8fpg5PUgoGB2CZlKeiLr/xU3bn2WnXjmnuCMvKzKw8k0EImfJZA031ZGvq/fv2Uuv/Vr6vn3v7n4Omh/sl6A5nL1dXnXa9Om/F+C6TKWwUSWN7cpu0ZEpiWWHnLI4HlzW3aniGBaYmVtzwSmI0EXvf2ZnXPjndzHzhfO2W+uvakE2PblaWg133msmD277kXX1df/PL9Ti0JRQItDJcySGAYwy93/0x945X/pv7npkfVoSMjwUvLTvlDddW5nw9uNdGj/+NRTwAJZEQYAkggY8EQQAIZC4YAEshYMASQwGwk8O4dO9WTu/fkPtD+XAvgJ2ecMKHdtzZvV6uu+6s66bvwsuvVl/7iKmeWhCKBFoZL2STQINkzult9e+MD6sHX1qkt+98Onj51xml6R9Fr1MoP/ImaMjDVAr1yVIEEliOPNnqBBNqgWI46kMBy5NFGL5BAGxTLUQcSmI0EujY67rj3MfXg4z+YENbHLlzqzJJQJNDCqCmrBBo0csP5f/rl3wf3HPz3rc8HTw9Pmq5WnLFKffbc69QpwwstUPS7CiTQ7/zZjB4JtEnT77qQQL/zZzN6JNAmTb/rQgK7QwJlKejyi5eplcsvqnVYloSu/sJXg41iXHgggRayUHYJDCN6/Z1X9C0mvq6Xin6n9vSHF/yB+r3ZZ6qzTzpPnTPnfLVo1gcsUPWrCiTQr3xlGS0SmCVdv+pGAv3KV5bRIoFZ0vWrbiSw/BJoZO+pR25XCxfMqRugsiT06isurZPDokYwEmiBfDdJoMG14+B29e03HlDfen2dendkRx3FaQPD6vx5H1Lnz/2QOm/OUvXBeR8OZg7L/EACy5zddH1DAtPxKnNpJLDM2U3XNyQwHa8yl0YCyy+BvoxfJNBCprpRAsPYZFfRl7a9oF7a+oJ69Z2X1IGx/ROovn/WGcFModyLUMTw7JOXWCDvThVIoDu5KDoSJLDoDLjTPhLoTi6KjgQJLDoD7rSPBCKBroxGJNBCJrpdAsMI5VYTm979iXp524vqlXc2BN9/tuunavzYeB1puf3E2SdrKZxzQVUMP6jmT1tgIRvFVIEEFsPdxVaRQBezUkxMSGAx3F1sFQl0MSvFxIQEIoHFjLyJrSKBFjKBBDaHeHDsgHp5+7+rV/RXIIfbNyhZThp9yE3qg5nCuR8M5PCck8/X9yocspCh7KtAArNn7EsLSKAvmco+TiQwe8a+tIAE+pKp7ONEApHA7EdZshaQwGScmpZCAtND/O3e31SkcPuLWgz/Xb2x41U1enS0rqL+3n694cxZWgi1FGo5XKLl0NWb1iOB6cdAWY9AAsua2fT9QgLTMyvrEUhgWTObvl9IIBKYftRkc4RzEii75uzcFX+zR1e2VI2mAgnsfHCOjY+pjTtfC2YKK7OF/65+vedXEyqeMTgzmCk8T8RQfy2d9/tObDqDBOh8JUAAACAASURBVHY+BspSAxJYlkx23g8ksHOGZakBCSxLJjvvBxKIBHY+iuzU4JQEyj015s+d7cxNFJMiRgKTkkpXbvfoe2rD1h9VlpLq2cJX9TLSvYfrPyDoUT36Bva/G8wUni9LSU9eqj5w0tmqr6cvXWMdlkYCOwRYosORwBIls8OuIIEdAizR4UhgiZLZYVeQQCSwwyFk7XCnJPDMj1yp1n3lJrXsgsXWOphHRUhgHpSV3nLmmPrFe5sqM4VaCl/Scrjp3Y1KbmYffsimM+fo3UeD2cJgGekFSq43zPKBBGZJ16+6kUC/8pVltEhglnT9qhsJ9CtfWUaLBCKBWY6vNHUjgRFashx17kmz1OPr1tReuXHNPerp9RuC389adGrda/IcEphmyNktO3JkJJghlGsLRQxFELcf3DahEbPpTOUWFR+0vukMEmg3rz7XhgT6nD27sSOBdnn6XBsS6HP27MaOBCKBdkdU+7U5JYGyHHT5xcvUyuUXtd+jDo4UAZRHWAIfffIZdd9D31Prn7g7eE1iXHruInXzNZfXWkICO4CewaFb9r+t71tYlUIth6/rTWcOaVkMP2TTmTNOXBwsHw2WkurdSH9n5ultR4MEto2udAcigaVLadsdQgLbRle6A5HA0qW07Q4hgUhg24PH8oFOSeBzL76uvvjl+2vCZbmvTaszAvr21h1qw6ubarN9UemLSqFUigTmman0bcly0Y07fhzcmiLYeEYvI/3V7p9PqMhsOiMbzshsYZpNZ5DA9Hkp6xFIYFkzm75fSGB6ZmU9Agksa2bT9wsJ7A4JFF+47a6H6wZI3GrC9CPI3hFOSaBcE9jskdXuoGHRu+Pex+okUGYHr77i0trspIjq6i98VYVjQQLtDci8apINZl7a+kIghCKGsqRUNqIJP2TTmd+ZcXqwG6lsOiNyKLOHcZvOIIF5Zc79dpBA93OUV4RIYF6k3W8HCXQ/R3lFiAR2lwSGfUE859YbVhW24jE6xp2SwLzegOF25Ho/edy59trge1QCowkzEvjUI7erhQvmFBEybWZEYNO7m9QLb7+gfvT2j9QLm19Qr21/TR0ZP1LX2pSBKer8eeer3z/l99WHTvmQ+oOFf6BOnnpyRhFRLQQg0AmB0bGjanAg352CO4mXYyEAAQhAoBwEzExgWALjLikrsrfOSaCRrDCULHcMbXRfwhNnnRAsS2UmsMjhWWzbh44cUj9+56VglvDFrf+mfrz9JbX1wJYJQcmmM0sXfEhd/oGValLPdNWr/zt5yhw1b/gUNVnvVMqjuwgwE9hd+W7WW2YCGQuGADOBjAVDgJnA42Nh/mx750g/3/VztXnv5twH2umzT1cLhhdMaDcqgW9t3q4u/tQtyqVJJKck0AALAzLQ8po+jc4Eck1g7u8npxvcefAdfWuKF4OlpBu2/UhL4ssTNp0Jd2D6pBPUvGkLgq85U+ep+VoM5fu8KfpLPzd32nw1e/JJTveZ4NIRQALT8SpzaSSwzNlN1zckMB2vMpdGArORwOueuk7d82JldV+ej69d/DV17QWV1YThR9w1gfI6EtggO9FZN1MsbjOWrBIclUB2B82KdHnqfUPvPvrGuy+pX+7+mXpr99tqy/7Natv+Lfr724k7uXD6+wIhnKtnFU8ZXqhOkplEkUf9uzwvz/HwgwAS6Eee8ogSCcyDsh9tIIF+5CmPKJHAbCTw7hfuVk/+nyfzSGFdG3/+oT9Xn/y9TzaUwPByUNcuKXNqJrDRzeLjNmPJKstRCZR2uE9gVrTLU2/cxjByc/sdB7arbXoJ6bYDW9VWkUP981YtiJXvFVncP7avJYjenl514pST1dzQDKIIY2WWsSKP86edoqYMTG1ZFwWyJYAEZsvXp9qRQJ+ylW2sSGC2fH2qHQnMRgJdGwNx1wRKjI1cp4j4nZJAF2YC20kCu4O2Q61cx3SyO6hIoBFCkUMzkxiIogikFkVZhipS2eoxPGl6IIRmqWmdKFYFcvaUk/S+pz2tquL1NgkggW2CK+FhSGAJk9pml5DANsGV8DAksHslsJEYFjXMnZJAF64JbCcRSGA71Mp1TCcSmISE7FK6vTqbaMQwEMeDFUmUn+X10aOjLaub1DdJb1wzt04U50+tXJ8YFsiB3oGWdVFgIgEkkFFhCCCBjAVDAAlkLBgCSGB3SWB05HNNYJO/BXnvDmrjzxISaIOi33VkLYFJ6ewa2RkIoRHFbQery1BDS1D3jO5uWZ3MFM6afGJFFGUjm5jZRZFG2fiGRz0BJJARgQQyBqIEkEDGBBI4cQzY3B2UEZaegFMzgenDd+MIJNCNPBQZhSsSmITByJERtXWf3sDmQOUaRZlJNNcovnNgW7ChjVzDmOQxWd83cZ6eRTQb2HzgxMXq7JOX1B06qXeSGh6crmSpqvlKUrevZZBAXzNnP25mAu0z9bVGJNDXzNmPm5nA40yRQPvjK02NSGAaWg3KIoEWIHpehU8SmBR1cG1iaAObdw5uU5v3/ba60U1lg5tDWijbeYSFcMbQTDV1YJqaNmlYzyyKKJ4wQRpFIqcPTNdl9Hf9s8u31UAC2xkR5TwGCSxnXtvpFRLYDrVyHoMEIoGujGwk0EImkEALED2voowSmCQl7x3aVZtNlGWoOw5uV++Nvqf2Hd6r9o3u0d/36a/qd/373rG9amTsYJKqE5WRHVOHtTwOa0EcHtTyKD9riRRRnD44Q03TchkWzum6jMimKTNtYDhRO2kKIYFpaJW7LBJY7vym6R0SmIZWucsigUigKyPcCQmU7VLlZvC33fVwUy7he224AlDiQAJdykYxsXSrBLZDW3Y5lesS941qURRZNMKod0mVn/dqWdyv5XGvlsfgu8iklsdAKEUkdZn9+ks2y+n0IbfeCEtibdlqIJXHl68aeZRrICvlj4vkCYMzldRjHkhgp1kpz/FIYHly2WlPkMBOCZbneCQQCXRlNDshga7AaDcOJLBdcuU5DgnMP5cHxw7UhFDEUH7ffUhmISszjwf0157DWjYj8rhXy6dIpM1ZyaH+oWB2UWYZZwydoJb9hz9QAz2t79koO7AO9g2qSfor+Llf/9yrv4Lvk5Rcc9nX0xe8Lru6mrLB77q8lMtiNjP/bJazRSSwnHltp1dIYDvUynkMEogEujKynZLARjdQlFtH3PfQ99T6J+52hVtdHEigk2nJNSgkMFfcVhuTJa1mNnL/2P5gBjI8Gykb6bw38u7xWUszexnMYO5Vu/XyV1ceMks5ILKoBbIikhWhlN8Htaj2izhqgRzUZY7LZuU58xW8FkipHFORUVOX/C6yWqmj/jWR1MpXpa7J/ZNdwVJYHEhgYeidaxgJdC4lhQWEBCKBhQ2+SMNeSKC5bQTLQV0ZNsQRJYAEdveYMLOSMgt5+NgBveD1sNq+V/88fliNHT0c3L/xcPB1OPguv4+Nj9V+Dr8mxxw+ostWv4+O65/171K+Vk/oNZvXWGaZxSEthTKrKV+9vZXvffJd9db93qufCZ6vlpOltuEy/b39ckTkGF0mdIzUKr/3yLEJ26nEpuuVduW4ILb6WKXdunZMf6ptRGMdHBhQw5Mnqb0Hjuibruj/9Nlf9HuFecxrlWdrx+gfGtZRK3esRx/Uulxde5F2JMZqREHbsbFJG7ot0x/TpsyI84gngAQyMgwBJBAJdOXd4IUE3nHvY+r7P3yemUBXRg1xTCCABDIoDIGirgmUGcuKTBrZ1DJZFUgRSiOQgUwGklmRy+C1QDKNoFa+175EVqvlR6sSeySo45AarQlrpS75XaRXZlZ5QCAtgYpEhoS4or0T5DlOdMPi2nNMXLhybGIpbtBWWNwnCHGonbj2wjH16vonDfSpsSPHTEspYqv/QCD+w4QIp5Ckh2Or49Tkw4fE3CIfCLSKzVw/Xf14IXmeqtSSxhX0M4jt+FhoFVvt9chYSNKmtJOknGSyV8c1e3hIvbtvNO1bpHTl/8vZf1S6PvnUocIlMO7m8HEA133lJrXsgsVOsmU5qJNpyTUoJDBX3E43VpQEugxFBHV8/KgK/tPfjx0bV0eP6Z/193H5Lq/p7+Eyx3/XZeuO0eUmHFMpE62jVk6NV9pM0M7xY6px1h3Toh3pQ9CnSv+O6a9evWfQ6JEjenZY/3cs+H/tu/5h4vP67DhcTv8W/C6PuDpqZfWrpmyzcnVtRtpqGVslgloclTF3TN8q5pDLw4/YIAABRwkcW1P528ajGAKFS2C4242uCSwGTfJWkcDkrMpaEgksa2bT9wsJTM+srEdwTWDrzBqBFLWsF9SJ0lyR4fjnK2p6XJQTC7Gewpkg51JPK2GPkfJmbfZq4Z4+tb86+9OkD9V2AxYhSQ90O+5Dgyi3iKgbsW/Jo9pWy3LV9ir6X/+hQKIPOlJyC/pczVHS2CojQUo3/+Ak/kOY+twkabP2EU2D/Ji81drTkcnfhsNHxlu/QUpe4t+u+peS99Dt7jklgW6jahwdEuhr5uzFjQTaY+l7TUig7xm0Fz8SaI+l7zVxTaDvGbQXP9cEHmc5fzYbiNkbWelrQgLTM5twBBJoAaLnVSCBnifQYvhIoEWYnleFBHqeQIvhI4EWYXpeFRKIBLoyhJ2SwLc2b1cXf+qWhmzYHdSVYUMcUQJIIGPCEEACGQuGABLIWDAEkEDGgiGABCKBrrwbnJLACy+7Xn3iox9WHz7/TPXFL99f2w10xeq1avnFy9TK5Re5wq0uDmYCnUxLrkEhgbnidroxJNDp9OQaHBKYK26nG0MCnU5PrsEhgUhgrgOuSWNOSaDZGGbhgpPVquv+qiaBsoNoWApdgWfiQAJdy0j+8SCB+TN3tUUk0NXM5B8XEpg/c1dbRAJdzUz+cSGBSGD+oy6+RSclUG4FIUJoln9ys3hXhgtxNCKABDI2DAEkkLFgCCCBjAVDAAlkLBgCSCAS6Mq7wSkJlGWfS89dpG6+5nIV/pmbxbsyXIgDCWQMtCKABLYi1D2vI4Hdk+tWPUUCWxHqnteRQCTQldHulARGochsoHk89cjtauGCOa5wq4uD5aBOpiXXoJgJzBW3040hgU6nJ9fgkMBccTvdGBLodHpyDQ4JRAJzHXBNGnNaAl2B1CoOJLAVofK/jgSWP8dJe4gEJiVV/nJIYPlznLSHSGBSUuUvhwQiga6Mcqck0GwMI9cE+vRAAn3KVjaxIoHZcPWxViTQx6xlEzMSmA1XH2tFAn3MWjYxI4FIYDYjK32tSGB6ZhOOQAItQPS8CiTQ8wRaDB8JtAjT86qQQM8TaDF8JNAiTM+rQgKRQFeGsFMS6Pr9ABslDQl0ZTgXFwcSWBx711pGAl3LSHHxIIHFsXetZSTQtYwUFw8SiAQWN/rqW3ZKAt/avL3u/oCuQGoVBxLYilD5X0cCy5/jpD1EApOSKn85JLD8OU7aQyQwKanyl0MCkUBXRrlTEhjeDTQOkLlvoCvwTBxIoGsZyT8eJDB/5q62iAS6mpn840IC82fuaotIoKuZyT8uJBAJzH/UxbfonASu+8pNKroxzKNPPqPue+h7av0Td7vCrS4OJNDJtOQaFBKYK26nG0MCnU5PrsEhgbnidroxJNDp9OQaHBKIBOY64Jo05oUEPvfi62r1F76qmAl0ZdgQR5QAEsiYMASQQMaCIYAEMhYMASSQsWAIIIFIoCvvBi8k8I57H1Pf/+HzzAS6MmqIYwIBJJBBgQQyBqIEkEDGBBLIGIgSQAKRQFfeFYVLoJnlawUkbploq2Pyep3loHmRdrcdJNDd3OQdGTOBeRN3tz0k0N3c5B0ZM4F5E3e3PSQQCXRldBYugWEQRd0sXm5N8camN2uhRIXzxjX3qKfXbwheP2vRqerxdWvq8ocEujKci4sDCSyOvWstI4GuZaS4eJDA4ti71jIS6FpGiosHCUQCixt99S07JYFFQREJNGJnZiafeuR2tXDBHBXdlEbKLj13kbr5mstr4SKBRWXOnXaRQHdyUXQkSGDRGXCnfSTQnVwUHQkSWHQG3GkfCUQCXRmNSGBMJsIzklHpi9upFAl0ZTgXFwcSWBx711pGAl3LSHHxIIHFsXetZSTQtYwUFw8SiAQWN/ocnwm88LLr1c5de2L55LE7qEjebXc9rMxMoMRz9RWXqpXLLwpiitupFAl0ZTgXFwcSWBx711pGAl3LSHHxIIHFsXetZSTQtYwUFw8SiAQWN/oclkCZdZs/d7a6c+21ufMJb1ATviZQZgVvvWHVBAk0kiiB7js4lnu8NOgWgYH+XtXb26NGDx91KzCiyZ1Af1+vGujvUSOjjIXc4UcaHB8/Frwvi3rIif/gpD518NCRokKgXUcI9Ogz/ylDferACGPBkZQUF4b+kzRtaEDtH+HccXjKQHF5oGXl1HLQojaGCY+DtzZvVxd/6hZlRDDJTOA+/qh3/VtJPvGXc81DY+Ndz6LbAfT39agBLYIjfCBQ+FBwQgIHetVBPhAofCwUHYD8+zBlsF/t5wOBolNRePvysdS0yf2Kc0elhjUHHsURQAJj2IfFj2sCixucPrXMclCfspVtrCwHzZavT7WzHNSnbGUbK8tBs+XrU+0sBz2erfmzJ/uUutLF6pQEinAtv3hZbellHrRlGejzL22s7fZploWamUB2B80jC/63gQT6n0NbPUACbZH0vx4k0P8c2uoBEmiLpP/1IIFIoCuj2CkJFAH74pfvV+ufuDtXPrIMNfzgPoG54i9FY0hgKdJopRNIoBWMpagECSxFGq10Agm0grEUlSCBSKArA9kpCYzKWBRSHruDtpMYdgdth1q5jkECy5XPTnqDBHZCr1zHIoHlymcnvUECO6FXrmORQCTQlRHtlAS6AiVtHEhgWmLlK48Eli+n7fYICWyXXPmOQwLLl9N2e4QEtkuufMchgUigK6MaCbSQCSTQAkTPq0ACPU+gxfCRQIswPa8KCfQ8gRbDRwItwvS8KiQQCXRlCCOBFjKBBFqA6HkVSKDnCbQYPhJoEabnVSGBnifQYvhIoEWYnleFBCKBrgxh5yRQdgh9Y9ObAR+zQYtcK/ixC5cWchP5JIlCApNQKncZJLDc+U3TOyQwDa1yl0UCy53fNL1DAtPQKndZJBAJdGWEOyWBIoDz584OZE/u1felv7hKLbtgsYrepsEVeCYOJNC1jOQfDxKYP3NXW0QCXc1M/nEhgfkzd7VFJNDVzOQfFxKIBOY/6uJbdEoCZcbvqUduVwsXzKmTQHPvPnYHdWXYEEeUABLImDAEkEDGgiGABDIWDAEkkLFgCCCBSKAr7wanJFBm/x7+2l9OkEBmAl0ZLsTRiAASyNhAAhkDUQJIIGMCCWQMRAkggUigK+8KpyTwjnsfU9//4fPBzeLNctCFC05WF3/qFnXlio+rm6+53BVudXGwHNTJtOQaFBKYK26nG2Mm0On05BocEpgrbqcbYybQ6fTkGhwSiATmOuCaNOaUBEqcZulnOOZbb1ilVi6/yBVmE+JAAp1NTW6BIYG5oXa+ISTQ+RTlFiASmBtq5xtCAp1PUW4BIoFIYG6DrUVDzkmgK2DSxIEEpqFVzrJIYDnz2k6vkMB2qJXzGCSwnHltp1dIYDvUynkMEogEujKynZLAG9fco55ev0FFN4DhFhGuDBfiaEQACWRsGAJIIGPBEEACGQuGABLIWDAEkEAk0JV3g1MSKNcBXn3FpROWfrIxjCvDhTiQQMZAKwJIYCtC3fM6Etg9uW7VUySwFaHueR0JRAJdGe1OSaDM+JkbxIcBcYsIV4YLcSCBjIFWBJDAVoS653UksHty3aqnSGArQt3zOhKIBLoy2p2SQGYCXRkWxJGWAMtB0xIrb3kksLy5TdszJDAtsfKWRwLLm9u0PUMCkcC0Yyar8k5JoCz7vO2uh2s3jJdOv7V5e3CLCJd3CGVjmKyGpz/1IoH+5CrrSJHArAn7Uz8S6E+uso4UCcyasD/1I4FIoCuj1SkJFChxt4iIWyLqCkCJAwl0KRvFxIIEFsPdxVaRQBezUkxMSGAx3F1sFQl0MSvFxIQEIoHFjLyJrTonga6ASRMHEpiGVjnLIoHlzGs7vUIC26FWzmOQwHLmtZ1eIYHtUCvnMUggEujKyEYCLWQCCbQA0fMqkEDPE2gxfCTQIkzPq0ICPU+gxfCRQIswPa8KCUQCXRnCzkmgbA6zc9eeWD7R+we6AhEJdCUTxcWBBBbH3rWWkUDXMlJcPEhgMvaHjh1Tx6Ro9bv8fKxHf1We0t8rrwdlYp6rldNn16bs8XL1x1fqM3Xp0tJO6LlqJLUy5rWgPhOHlA/aqn+uUiZSrvpcb2+Pmj61X72773ClD+Gyde1XYxMWtTaq/Q76fpxFPZPQ8zV2EW5N2YX4NitXi77ahzp+oRiqfarjFvQ5wq0Ko45HrX3ThhkfYW5VzhF2wbN1+akvF2ZfGSMRnkF/Yrg1YB+UFasL5zQyZs34rERSGdvyt+HwkfFKYrv48W9nnN7FvS++605J4IrVa9X8ubPVnWuvLZ5MigiQwBSwSloUCSxpYtvoFhI4EdqIPlEa119H9UtH9WmQ/CynP0f1CeG4/v2ovKbPjoKfq2WO/94TvC6vjesyR/UJVKWu6jHysz4JC+rQxwav6TLBMVK+etx49blamUg74zqWI6aNoI7KyWQQT7V9iTe2nWqZoJyOxbQrp3y9fb1qVJ/s1clB9SS3cuIYOXEOiUxFjCIn59Hypo7KeegEoTkuUiHxMeVCkhHEUSdI8SfZFWkyqlRpUySOBwQgAIG0BI6dd3baQyhvkYBTEtjoPoEW+5tJVUhgJli9qhQJ9CpdmQZblASKaB0+Nq6/5LuqfNfqZH4f1Sf5h8fH1Zju/WjwWuV3mZcY1dYyJr/r40f174f17+HXK3Udqxynv0SW5JigzupzwbG63TH9+z79Mw8IpCUwVJ1REZ/t0T9rD9faHfrSv+hng2prz1efqzpw5Xk9BqUqU7K+bKWO8HNBfUFb1eerx1bKHC8rBx2PR+I73o4pG9RVa9uU71G9+vlJA71qTH8gYPoV1F5XZ7W90HOV2KSdUMym/6Hn6pnEsTOE4pnW+lXtkyEdx9k8V+OhPxloxjjMrbeWFZkTPM6nFn8ce3muOhZqccVwi2UvxzVhPIFbbSxUYwuPybpxF8mHfBASyVF9n+rLnzh9SM8Ky1/b7n78l/mzuhtAwb1HAi0kAAm0ANHzKpBAzxPYYfgHtUzt1ScB+8ePqtF+PVOiZ3+2HxgNBEmkqF6WKjIlMiavj4qoBdIl8iVSVhWr6u8VeauUHxXJC16vHqt/F/nz4SEn+H36tK9PB9srP+sTs8rverZMvlef69UnlOY1vYKu8lqtjFL9ulyv7nJQR3CclJG69Je8FvxcaUdOnM3zde1IDNVjpay0E47l+O+mnWrdut0gNtOO1F+LpVpW58P0b7CvTw1P7ld79481PEk2uTt+olo9WawTjJD4NBCU4PhgKMSdkNcLzvET3ziRqUQUyE/VG+rEx7QfEjRz8mskzofxmHeMXBOYN3F32+OawOO5mT97sruJ6oLInJJAWQ66/OJlauXyi7xCjwR6la5MgkUCM8GaS6XvaXHbp+VLvvYHMqe/9HP79Um4fDevVb4fVQe0kOw5qn8+Jq8dU7v1c648hnt61YA+wxjUJ+qT9M+T9M/yNahP6Ae12fTr78Hv8nz150k6+Em9uqz8rssMyu9ST/WY8Oty3IB+fbBWt5SVukxb1d/165ONQbgCp4A4uCawAOiONokEOpqYAsJCApHAAoZdbJNOSaDcI/CLX75frX/iblf4JIoDCUyEqdSFkMD80xuefRNxk993V0VN5OyAlrQ98rt+TQRubyB5R/X3ynJFETpbs2gyAzKtV8/6aImaoWd/puuZwB59YdqkQMYqElYRLSNJImIiWiJUFTE7Ll2V3+Vrsj5GZppEzkTSTNlA8AJ50+3qn3m4SQAJdDMvRUSFBBZB3c02kUAk0JWR6ZQEyjWBzR7sDurKsCGOKAEkMPmYkBVre2Jm2PbpGbaKrOlZOC1qtVk4EbaqyIVl7kjyJhuWFH0a1qIlIjU9kLjeypf+PfyzvDZNdveTMsFrPfqr8rNIX1jDirom0AIOqrBMAAm0DNTj6pBAj5NnOXQkEAm0PKTars4pCWy7FwUfyExgwQlwoPlulUBZSrnt6BG17cgRtVV/33H0qHovWCZZmX0zSyhl9i2YhbN8DduJff0VIdMaFkhZVc6ma5Ez4hZIW1Xqas/pslImi1k0JNCBN6QjISCBjiTCgTCQQAeS4EgISCAS6MhQlE2n3NpVQJaErv7CV+v4rPvKTWrZBYtdYTYhDiTQ2dTkFlgZJXCLyJ3+2qoFT76/owVv85ExtW38+PPtbg0fnmmboYVsaiBkeqYtkLXQjFwwM1eRO3ltWo/IW4+areXP1QcS6Gpm8o8LCcyfuastIoGuZib/uJBAJDD/URffolMS+OiTz6jb7npYPfXI7WrhgjlBxG9t3q4u/tQt6tYbVjm7YQwS6MpwLi4OnyRQroPbqmXOSF5U9LYc1aKnhS/JQzb/mKeFTL7m6q/5/QPqhOpyymCmrTojV1s+WX0tSd2+lkECfc2c/biRQPtMfa0RCfQ1c/bjRgKRQPujqr0anZLACy+7Xl19xaUTZE/k8L6HvufshjFIYHuDr0xHuSKBu7S8ybJM+QqWaFZn7bYdqSzb3KoFTzZLafWQLeFn6WvdArnTM3Hz+kX0BgLRm9svz1d+FtHjUU8ACWREGAJIIGPBEEACGQuGABKIBLrybnBKAhvdLN4sEWVjGFeGDXFECWQtgbIJikhdReSOL9EMlmZWr8fbrp+X+8+1esiukycbwQukTs/g9VZm8uRnM6sntxrgkZ4AEpieWVmPQALLmtn0/UIC0zMr6xFIIBLoyth2SgKZCXRlWBBHWgKdSKDsfBkreIHwVWbwZ/bVMgAAHklJREFUduqv1nqnd7rU182FRc78XBG7ysyeXE+H3qXNcPLySGByVmUviQSWPcPJ+4cEJmdV9pJIIBLoyhh3SgKLuiZQ5HPnrj21nEQ3orlxzT3q6fUbgtfPWnSqenzdmrr8sRzUleFcXBxxEijStqO2LLN+F02zXFMEb7+WwFYPWXQpO2GapZnBrJ1cixeauZuvf5/C8sxWKDN/HQnMHLE3DSCB3qQq80CRwMwRe9MAEogEujJYnZJAgZL37qCy8cyd3/iuunPttUFOjIiapafR6xFXrF6rlp67SN18zeW1HCKBrgznYuJ44/Ah9cb4mPrlkcPqN4cO13bTlI1Xkj4W6g1VKmLXp07RP5+kl2eGBe8ULXs8/CCABPqRpzyiRALzoOxHG0igH3nKI0okEAnMY5wlacM5CUwSdJZlzG6kZofSqPTFbVKDBGaZEbfqlmWZLx0eVS+NjqgNo4fUj7UANrtNgtzWwGyqMkcLnuyeKd/nBZutiPj1OX27A7fo+xENEuhHnvKIEgnMg7IfbSCBfuQpjyiRQCQwj3GWpA0kMEIpuglN9DrFuE1qkMAkQ82/MiJ3P9ai96oWvRe19InwyTV60YeI3AenTFXnTx5SM4/qWybIEs1gJm9AyS0UeHQXASSwu/LdrLdIIGPBEEACGQuGABKIBLrybnBCAs0SzLh7ATZ7LQuIIn2f+OiHa8s9ZcfScFxGAsP3MswiDurMn8Cm0VH1woGD6kf6S76/NjKiogs6ReqWTpmiPjR1ivr96teCgYH8g6VFCECgJYHRsaNqcKCvZTkKQAACEIAABLqNgBMSKEsu58+dXbsuL5oE2Zhly7Z3J2zIYjtZIoBLFp9eFwczgbYpu1HfXn2vPFnS+bJe2vlydbZv93j9LJ/M4Z2qZ/OWDE5W5w8OqfMmDaoPTBpScaeUnewO6gYRorBFgJlAWyT9r4eZQP9zaKsHzATaIul/PcwEHs/h/NmT/U+oxz1wQgIb3R/QcM3jPoFxAijtc02gx6O7Grqo3UYteq+MVYTvZS1/vzoyNqFjM/R1eucFsjeklmjhWzo0ObjlQpIHEpiEUneUQQK7I89JeokEJqHUHWWQwO7Ic5JeIoFIYJJxkkcZJFBTFgm9csXH63b8NPDZHTSPYWi3DdmV8yWRPhG+wyPqdT3bF928Rfba/D0tektE+LT4yddp/ZPaDgQJbBtd6Q5EAkuX0rY7hAS2ja50ByKBpUtp2x1CApHAtgeP5QOdkECZhfvSX1ylll2wOLZ7MhP4xS/fr9Y/cbfl7h+/JUS04o9duLS2LJT7BFrHbq3CEb15y6sifHrTFpnhk5m+7ZFlndLYPL15iyzrFNk7b2BInaO/D1nctAUJtJZS7ytCAr1PobUOIIHWUHpfERLofQqtdQAJRAKtDaYOK3JCAu+49zG14dVNDa/5a3XNYIcMOj6c3UE7RpioArn5+i/GDgei94qe4ZNbNWzSX9H9OkXuztEzfMGyzuos3zx9D74sH0hglnT9qhsJ9CtfWUaLBGZJ16+6kUC/8pVltEggEpjl+EpTtxMSKAHLbKA8orN98vzOXXuUuXl7ms7lVRYJzIa0bNQi9+J7Wc/yydJOmfHbe2y8rrE0m7dkE2WlViQwS7p+1Y0E+pWvLKNFArOk61fdSKBf+coyWiQQCcxyfKWp2xkJlKBlRvDBx39QF394WWaajuVZFgnsnPaYXta5MZjlqyzplOWdv7a8eUvnUTauAQnMkq5fdSOBfuUry2iRwCzp+lU3EuhXvrKMFglEArMcX2nqdkoC0wTuUlkkMH02fqsFr7JxS+XrDb2sc1SLYPhhe/OW9FEmPwIJTM6q7CWRwLJnOHn/kMDkrMpeEgkse4aT9w8JRAKTj5ZsSyKBFvgigc0hHtT35AuWdAabt8j3UbVD7+AZfZjNWyq3aLC/eYuFVDesAgnMkq5fdSOBfuUry2iRwCzp+lU3EuhXvrKMFglEArMcX2nqRgLT0GpQFgk8Dkbm8jbp+/FVZvkq9+X7mf69/ko+FezMGd68Re7PNz/jzVsspBoJzBJiSepGAkuSSAvdQAItQCxJFUhgSRJpoRtIIBJoYRhZqQIJtICx2yXwf40cCDZwkZm+l/Q1fXLbhujj/QOT1NnV3TpF/s7V9+gr04OZwDJls7O+IIGd8SvT0UhgmbLZWV+QwM74leloJBAJdGU8I4EWMtGNErhTL+d8aP8e9fD+veqdyNLOaT296oNDk7XoDamlWvbO1z8P6+fK/EACy5zddH1DAtPxKnNpJLDM2U3XNyQwHa8yl0YCkUBXxjcSaCET3SSBb+jZvv++b4/6nwf21sgt1Us5ZZbvbC188rVooFyzfEmGCBKYhFJ3lEECuyPPSXqJBCah1B1lkMDuyHOSXiKBSGCScZJHGSTQAuWyS6DcjP2pg/vVA/t2qxf1ck95yMzeimnT1WeHZ6pT+rO9EbuFFGVeBRKYOWJvGkACvUlV5oEigZkj9qYBJNCbVGUeKBKIBGY+yBI2gAQmBNWsWFklcI++Wfu39ZLPB/XM35bqks9TBwbUp6fNUCunTldTesu9xDPN0EAC09Aqd1kksNz5TdM7JDANrXKXRQLLnd80vUMCkcA04yXLskigBbplk8BfHjmsvrHnPfW3B/fVNnlZNjRFXTU8Q/2nyVNVjwVmZasCCSxbRtvvDxLYPruyHYkEli2j7fcHCWyfXdmORAKRQFfGNBJoIRNlkEDZz/MZveTzfj3z978PHQyoyG0cLtMzfldPn6FO659kgVR5q0ACy5vbtD1DAtMSK295JLC8uU3bMyQwLbHylkcCkUBXRjcSaCETPkvg/mPj6jG9w+eD+3erN8fGAhpy0/Y/0bN+fzztBDWzt88CofJXgQSWP8dJe4gEJiVV/nJIYPlznLSHSGBSUuUvhwQiga6MciTQQiZ8lMDfHBkLNnp5XAugiKA8lugdPmXJ5yVThxVbvaQbGEhgOl5lLo0Eljm76fqGBKbjVebSSGCZs5uub0ggEphuxGRXGgm0wNYnCfwXvdTzgb3vqWf0d1kCOqC/LpkyrD43faY6s2Q3cLeQ2sRVIIGJUZW+IBJY+hQn7iASmBhV6QsigaVPceIOIoFIYOLBknFBJNACYNcl8NCxY+q7+r5+39Qzfz8bOxz0eJZe5inLPT89fII6uY95v06HARLYKcHyHI8ElieXnfYECeyUYHmORwLLk8tOe4IEIoGdjiFbxyOBFki6KoHbjh7VSz7fU9/RSz5369s9yOMMfSP3z2jxkw1fBuUvEQ8rBJBAKxhLUQkSWIo0WukEEmgFYykqQQJLkUYrnUACkUArA8lCJUigBYiuSeALh0bUA3qjlx/o3T5F/eRufnJrh8/o6/3kVg887BNAAu0z9bVGJNDXzNmPGwm0z9TXGpFAXzNnP24kEAm0P6raqxEJbI9b3VEuSOBhveTz7/V9/e7XSz7fODwaxDfc06tWTJuuPjs8U53Sz5JPC6luWAUSmCVdv+pGAv3KV5bRIoFZ0vWrbiTQr3xlGS0SiARmOb7S1I0EpqHVoGyRErjz6BF9e4e96mF9fz/5WR7v6x9Qf6pn/VbqJZ9TemUekEfWBJDArAn7Uz8S6E+uso4UCcyasD/1I4H+5CrrSJFAJDDrMZa0fiQwKakm5YqQwI16tu8+vcvnP+jZv8rd/VSw1FNu8SBLP7naz0JiU1SBBKaAVfKiSGDJE5yie0hgClglL4oEljzBKbqHBCKBKYZLpkWRQAt485JAmef7pwN6yae+3u+l0UNB5EP6r4ls8nL19BnqtP5JFnpDFe0QQALboVbOY5DAcua1nV4hge1QK+cxSGA589pOr5BAJLCdcZPFMUigBapZS+B7emfPb+vlnn+jr/fbqnf8lMe8vj51xbQZapXe6XOmvt0Dj2IJIIHF8nepdSTQpWwUGwsSWCx/l1pHAl3KRrGxIIFIYLEj8HjrSKCFTGQlgXJPv29o8fs7fY8/udefPJZMGgqWfF4ydVix1YuF5FmqAgm0BLIE1SCBJUiipS4ggZZAlqAaJLAESbTUBSQQCbQ0lDquBgnsGKFSNiVwXMfzw5ED6gF9vd+/jo4E0Q3or0umDKvPTZ+pzpw0aCFiqrBNAAm0TdTf+pBAf3NnO3Ik0DZRf+tDAv3Nne3IkUAk0PaYarc+JLBdcqHjbEjg/mPj6jv79qhv6WWfbx2pbPUySy/z/ONpJ6hP6yWfJ/cx72chVZlVgQRmhta7ipFA71KWWcBIYGZovasYCfQuZZkFjAQigZkNrpQVI4EpgcUV70QCf6OF7xt61u+7esOXA1oE5XHGwKC+sfsJwYYvg/LXgofzBJBA51OUW4BIYG6onW8ICXQ+RbkFiATmhtr5hpBAJNCVQYoEWshEOxL4rCz51LN+/6y/y9V+cjc/ubXDZ/T1fnKrBx5+EUAC/cpXltEigVnS9atuJNCvfGUZLRKYJV2/6kYCkUBXRiwSaCETSSVwRG/u8j/0Ji/f3Ltb/eLI4aDl4Z5etWLadPXZ4ZnqlH6WfFpIRyFVIIGFYHeyUSTQybQUEhQSWAh2JxtFAp1MSyFBIYFIYCEDL6ZRJNBCJlpJ4JajR9T9WvweO7BH7RmvLPl8X/+AvtZvhvqvesnnlF6ZB+ThMwEk0Ofs2Y0dCbTL0+fakECfs2c3diTQLk+fa0MCkUBXxi8SGMrEitVr1dJzF6mbr7m8Lj83rrlHPb1+Q/DcWYtOVY+vW1P3eiMJ/De9u+c39fV+/69e8llRP6X+4+AU9Vl9Y3dZ+snVfq68DTqPAwnsnGFZakACy5LJzvuBBHbOsCw1IIFlyWTn/UACkcDOR5GdGpBAzTEseVeu+HidBD765DPqvoe+p9Y/cXdAPE4UwxI4qpd8yn395Hq/nxweDY4Z0u942eTlai1/p/VPspM5anGKABLoVDoKDQYJLBS/U40jgU6lo9BgkMBC8TvVOBKIBLoyIJHAUCYuvOx69YmPfrhOAqPSF5VCOVwk8B295FNu7/BtfZuHXeNHg1rn9fWpK6bNUKv0Tp8z9e0eeJSXABJY3tym7RkSmJZYecsjgeXNbdqeIYFpiZW3PBKIBLoyupHAFhIoYnj1FZeqlcsvCko+9+LravUXvqo2Pvtg8PuGgyPqtre3qb/Xt3gwjyWThtRV+nq//2vqsCt5Jo6MCSCBGQP2qHok0KNkZRwqEpgxYI+qRwI9SlbGoSKBSGDGQyxx9UhgCwk88yNXqltvWDVBAp965Ha1cMEc1fPKa0ENA/pr+fTp6vpZJ6qzBwcTJ4CC5SAw0N+rent71Ojhyiwwj+4l0N/Xqwb6e9TIKGOh6FEwPn4seF8W9ZAT/8FJfergoSNFhUC7jhDo0Wf+U4b61IERxoIjKSkuDP0nadrQgNo/MlZcDI60PDxFzp55FEUACWwhga1mAj/2i1+pc/TM3+qZs9QcbvFQ1DguvF35xF/ONQ+NmS2ACg+JAAoi0N/Xowa0CI7wgUBBGTjerBMSONCrDvKBQOFjoegA5N+HKYP9aj8fCBSdisLbl4+lpk3uV/v4QEANaw48iiOABLaQwKTXBBaXQlp2gQDLQV3IghsxsBzUjTy4EAXLQV3IghsxsBzUjTy4EAXLQY9nYf7syS6kpGtjQAJbSGDa3UG7diR1eceRwC4fAKHuI4GMBUMACWQsGAJIIGPBEEACkUBX3g1IoM5E+BYRJjFm4xf5vd37BLqSZOLIngASmD1jX1pAAn3JVPZxIoHZM/alBSTQl0xlHycSiARmP8qStYAEJuPUtFSjm8VbqJoqPCGABHqSqBzCRAJzgOxJE0igJ4nKIUwkMAfInjSBBCKBrgxVJNBCJpBACxA9rwIJ9DyBFsNHAi3C9LwqJNDzBFoMHwm0CNPzqpBAJNCVIYwEWsgEEmgBoudVIIGeJ9Bi+EigRZieV4UEep5Ai+EjgRZhel4VEogEujKEkUALmUACLUD0vAok0PMEWgwfCbQI0/OqkEDPE2gxfCTQIkzPq0ICkUBXhjASaCETSKAFiJ5XgQR6nkCL4SOBFmF6XhUS6HkCLYaPBFqE6XlVSCAS6MoQRgJdyQRxQAACEIAABCAAAQhAAAIQyIEAEpgDZJqAAAQgAAEIQAACEIAABCDgCgEk0JVMEAcEIAABCEAAAhCAAAQgAIEcCCCBOUCmCQhAAAIQgAAEIAABCEAAAq4QQALbzMSNa+5RT6/fEBx91qJT1ePr1rRZE4f5QiBNzh998hl1210P17rGGPEly8niTDMWwjWacbHuKzepZRcsTtYYpZwmEH2vb3z2wZbxhsfPxy5cqu5ce23LYyjgPoG0Y+GOex9TDz7+g1rH+Lvgfo5tRbhi9Vq19NxF6uZrLrdVJfVAIDUBJDA1MqXkD/19D31PrX/i7uBo3sxtQPTskLQ5l5O8G//s/1YLF8wJenrhZderT3z0w/zB9yzvceGmHQumDnPczl17FCd7JRgIugtvbd6uLv7ULeqpR24P3utyUr/h1U1NPxSUfy/kwQeH5RgDphdpx4IRRvOhQfT3ctGhN4ZA+AOgK1d8nHMChkahBJDANvBHpS96UthGlRziOIFOcy5/+OXBJ/6OJzpBeO2MhedefF198cv3Bx8cnfmRK5HABJx9KBKVvqgIRPsQHgc+9I8YkxNIOxbSlk8eCSV9IMAHwz5kqfwxIoFt5FjevFdfcalaufyi4Gj5h331F76qkiwDaqM5DnGAQKc5Z7bYgSRaCiHtWIj+fUACLSXCgWriPtxpll858f/+D59XMhtsHrfesKr2b4kDXSKENgmkHQvSjIyVE2edEHw4JGNjy7adfFDYJn/fDkMCfctYOeNFAtvIq/zhDv/DbU7yzJKgNqrkEMcJdJJzlvk4ntyU4aUZC3EzQ0hgSuAOF5cPd+bPnV134h4dH+HwRRS2bHu3thTU/NvB8mCHk5wwtLRjQaqVY7bt2FX7UIBxkBB2CYohgSVIYgm6gAS2kcS0MwFtNMEhjhFoN+dGAPmAwLGEdhBOmrEQ3Sgi3CwzQB0kwZFD087+xJWPjidHukYYKQmkHQvRmT8+EEgJ3PPiSKDnCSxJ+EhgG4ls55qgNprhEIcItJNzBNChBFoMpZ2xEG6emUCLySi4qrTXdcVtHIMEFpxES82nHQtxlwg0m0W2FCbVOEIACXQkEV0eBhLYxgBod3fANpriEEcItMq52erbXBca/d2RbhCGBQJpx0K0SSTQQhIcqaLVjpDRvwOmvFn2x1JxRxJpIYxWYyE60xcdG3xoaCEJHlWBBHqUrBKHigS2mdx27xPWZnMc5gCBZjmP/oMuf+DDmz+Y8FkW6kAiLYSQZiwggRaAO1xFs3vDxX0YZGSAvwkOJ7XN0JqNhbjlnuG/I9Ik/z60Cd6jw6I5l9DZVNCjBJYsVCSwZAmlOxCAAAQgAAEIQAACEIAABJoRQAIZHxCAAAQgAAEIQAACEIAABLqIABLYRcmmqxCAAAQgAAEIQAACEIAABJBAxgAEIAABCEAAAhCAAAQgAIEuIoAEdlGy6SoEIAABCEAAAhCAAAQgAAEkkDEAAQhAAAIQgAAEIAABCECgiwgggV2UbLoKAQhAAAIQgAAEIAABCEAACWQMQAACEIAABCAAAQhAAAIQ6CICSGAXJZuuQgACEIAABCAAAQhAAAIQQAIZAxCAAAQgAAEIQAACEIAABLqIABLYRcmmqxCAAAQgAAEIQAACEIAABJBAxgAEIAABCEAAAhCAAAQgAIEuIoAEdlGy6SoEIAABCEAAAhCAAAQgAAEkkDEAAQhAAAIQgAAEIAABCECgiwgggV2UbLoKAQhAAAIQgAAEIAABCEAACWQMQAACEIAABCAAAQhAAAIQ6CICSGAXJZuuQgACEIAABCAAAQhAAAIQQAIZAxCAAAQgAAEIQAACEIAABLqIABLYRcmmqxCAAATaIXDHvY+pBx//wYRDNz77oHr0yWfUbXc9rG69YZVaufyiujI3rrlHPb1+g3rqkdvVwgVzVLN6nnvxdbX6C19tGp60IQ9pL+5hYjAxnbXoVPX4ujV1RU0Mca+Zgmd+5MqmcXzswqXB69I387hyxcfVzddcrt7avF1d/KlbgqeFT/gR7qN5zcTarD/t5IxjIAABCEAAAs0IIIGMDwhAAAIQaEjAiFxUaOT5C847IxA/I1bhMkZ41n3lJrXsgsUqST3hIKT8lm3vTpA4I03ReMLHhsXKCKh5/cLLrlc7d+1RzSQwrq649kwfw22EJTAqxoZBWBCT9IfhCQEIQAACELBNAAm0TZT6IAABCJSIgMyKmVmuZt0SuVqy+HR159prg2LR35PWY9roVALve+h7QTzyMDGJcMnzc0+aFTwfnSWM618zSWsmgcLs+z98Xq1/4u6gWiOHMosoM4jRmcBmUlui4URXIAABCEDAEQJIoCOJIAwIQAACLhKIylyjGMMzf8+/tLFOgOKksFVfbUjgw1/7y2BpppmpW7F6rVp+8TL15FPPZS6B0uaq6/5KXX3FpbXZ0i3bdqr5c08MltYiga1GAK9DAAIQgECWBJDALOlSNwQgAAHPCcRds9Zo1iq83NEsAzXdT1OPHNNKAuOwhq8JlBk/mYUT8Zs/d7b6z5dcqL745ftrz8nxWc4EigQ+/r1/rsmwzISa5+IksFl/PB9ChA8BCEAAAg4SQAIdTAohQQACEHCRgLmezsQWFT15XmRHljyaJZhx/UhSTysJbHVNoJFAM0Mp1wDKLKBcwyhimIcEymY4wkPaFhEVJtHrJ7km0MWRTkwQgAAEyk8ACSx/jukhBCAAAesEROTk2rrobFrS5aMmoEb12JJAaUek741Nb9aWYOYpgUb6jDAjgdaHIhVCAAIQgEAbBJDANqBxCAQgAIFuICCbmciSRrn1QfTRSKTiJLCdemxKoMy2vb11R60feUqg9P3m/+e+miwjgd3wzqGPEIAABNwngAS6nyMihAAEIFAIAbOj5YmzTqjtcimBmCWMcctBG0mgbNCSph6bEphUYOMgt7s7aPTWFKZuJLCQoUyjEIAABCAQIYAEMiQgAAEIQKApgbibpzeSnGbLQdPU00oC4wKO2xgmrpyNmcDwJjjSRvRm8WklsFl/GJ4QgAAEIAAB2wSQQNtEqQ8CEIAABCAAAQhAAAIQgIDDBJBAh5NDaBCAAAQgAAEIQAACEIAABGwTQAJtE6U+CEAAAhCAAAQgAAEIQAACDhNAAh1ODqFBAAIQgAAEIAABCEAAAhCwTQAJtE2U+iAAAQhAAAIQgAAEIAABCDhMAAl0ODmEBgEIQAACEIAABCAAAQhAwDYBJNA2UeqDAAQgAAEIQAACEIAABCDgMAEk0OHkEBoEIAABCEAAAhCAAAQgAAHbBJBA20SpDwIQgAAEIAABCEAAAhCAgMMEkECHk0NoEIAABCAAAQhAAAIQgAAEbBNAAm0TpT4IQAACEIAABCAAAQhAAAIOE0ACHU4OoUEAAhCAAAQgAAEIQAACELBNAAm0TZT6IAABCEAAAhCAAAQgAAEIOEwACXQ4OYQGAQhAAAIQgAAEIAABCEDANgEk0DZR6oMABCAAAQhAAAIQgAAEIOAwASTQ4eQQGgQgAAEIQAACEIAABCAAAdsEkEDbRKkPAhCAAAQgAAEIQAACEICAwwSQQIeTQ2gQgAAEIAABCEAAAhCAAARsE0ACbROlPghAAAIQgAAEIAABCEAAAg4TQAIdTg6hQQACEIAABCAAAQhAAAIQsE0ACbRNlPogAAEIQAACEIAABCAAAQg4TAAJdDg5hAYBCEAAAhCAAAQgAAEIQMA2ASTQNlHqgwAEIAABCEAAAhCAAAQg4DABJNDh5BAaBCAAAQhAAAIQgAAEIAAB2wSQQNtEqQ8CEIAABCAAAQhAAAIQgIDDBJBAh5NDaBCAAAQgAAEIQAACEIAABGwTQAJtE6U+CEAAAhCAAAQgAAEIQAACDhNAAh1ODqFBAAIQgAAEIAABCEAAAhCwTQAJtE2U+iAAAQhAAAIQgAAEIAABCDhMAAl0ODmEBgEIQAACEIAABCAAAQhAwDYBJNA2UeqDAAQgAAEIQAACEIAABCDgMAEk0OHkEBoEIAABCEAAAhCAAAQgAAHbBJBA20SpDwIQgAAEIAABCEAAAhCAgMMEkECHk0NoEIAABCAAAQhAAAIQgAAEbBNAAm0TpT4IQAACEIAABCAAAQhAAAIOE0ACHU4OoUEAAhCAAAQgAAEIQAACELBNAAm0TZT6IAABCEAAAhCAAAQgAAEIOEwACXQ4OYQGAQhAAAIQgAAEIAABCEDANgEk0DZR6oMABCAAAQhAAAIQgAAEIOAwgf8fwOpIAnh4uiYAAAAASUVORK5CYII=", "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": 19, "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": "iVBORw0KGgoAAAANSUhEUgAAA4EAAAFoCAYAAADpfv6vAAAgAElEQVR4Xu29C5Qd1Xmg+/dL/dATSSAhYRlsE/EQNgYZm1i2mBAcsB1GZuIAgyFMTCwYg4Hhgr1i1mLpLpxrYLGACwkoIRMCg4H4DmEYJzAQEjQhQ0zABsTDgHkJBAi90aO71VLr7r+6d6u6+pw+p07vqtq7zndYh/Oq2vvf379Pq76zd+1q2Wtuwg0CEIAABCAAAQhAAAIQgAAEmoJACxLYFHmmkRCAAAQgAAEIQAACEIAABCICSCAdAQIQgAAEIAABCEAAAhCAQBMRQAKbKNk0FQIQgAAEIAABCEAAAhCAABJIH4AABCAAAQhAAAIQgAAEINBEBJDAJko2TYUABCAAAQhAAAIQgAAEIIAE0gcgAAEIQAACEIAABCAAAQg0EQEksImSTVMhAAEIQAACEIAABCAAAQgggfQBCEAAAhCAAAQgAAEIQAACTUQACWyiZNNUCEAAAhCAAAQgAAEIQAACSCB9AAIQgAAEIAABCEAAAhCAQBMRQAKbKNk0FQIQgAAEIAABCEAAAhCAABJIH4AABCAAAQhAAAIQgAAEINBEBJDAJko2TYUABCAAAQhAAAIQgAAEIIAE0gcgAAEIQAACEIAABCAAAQg0EQEksImSTVMhAAEIQAACEIAABCAAAQgggfQBCEAAAhCAAAQgAAEIQAACTUQACWyiZNNUCEAAAhCAAAQgAAEIQAACSCB9AAIQgAAEIAABCEAAAhCAQBMRQAKbKNk0FQIQgAAEIAABCEAAAhCAABJIH4AABCAAAQhAAAIQgAAEINBEBJDAJko2TYUABCAAAQhAAAIQgAAEIIAE0gcgAAEIQAACEIAABCAAAQg0EQEksImSTVMhAAEIQAACEIAABCAAAQgggfQBCEAAAhCAAAQgAAEIQAACTUQACWyiZNNUCEAAAhCAAAQgAAEIQAACSCB9AAIQgAAEIAABCEAAAhCAQBMRQAKbKNk0FQIQgAAEIAABCEAAAhCAABJIH4AABCAAAQhAAAIQgAAEINBEBJDAJko2TYUABCAAAQhAAAIQgAAEIIAE0gcgAAEIQAACEIAABCAAAQg0EQEksImSTVMhAAEIQAACEIAABCAAAQgggfQBCEAAAhCAAAQgAAEIQAACTUQACWyiZNNUCEAAAhCAAAQgAAEIQAACSCB9AAIQgAAEIAABCEAAAhCAQBMRQAKbKNk0FQIQgAAEIAABCEAAAhCAABJIH4AABCAAAQhAAAIQgAAEINBEBJDAJko2TYUABJqPwJEnnCtfWbpYblhxYfM1nhZDAAIQgAAEIFCRQLASeOlVt8gjq54e06jZM6fLqvtv8iLd1Q6+7nngMbn6xrvk3NNPlssvOMOLWJNBnL58hbzwypviE0+XoGwOKpWZ1wGz9o/kLa+6XbKst6wsZeSJp1bL8iuulysvOVvOXHbiqJDs34qV114mS447qt5wvdxOv5cfrN+U6m9cvdyz5GTzE4f60N3XyIL5c7zkTFAQgAAEIACBshMIXgKTBxL2wNqHA4xQJXDN2nVyylnfH+n7rg+erWCmLdfu5yK3VgKT0mAPVrOUMXuwnawjzv3Fx+8o3d+eemWkkYYjgdWp1cs9Kwms9F277tZ75Y77HhYX3+VG+gv7QAACEIAABJqdQOkk0B4M+jDKVu/Bl2+d0B6gqaTp6IorIYqP3jZ68GdFUJk1WobuW00C9bNGRlvqzeF49doytP77Vl5Vb5HBbJfl92E8CQwGUB2BNtI3s+ReR8iy9LSLZe7+M8f0aX3/mKMOZZpqPRDZBgIQgAAEIOCYQFNJYHIK6aKFh1Q8MNmwaesozNWmRFaaUmiFqdJUP1tOtQPWeqZMxfd96pcvj5oSW2kqXCP9JX7QZqVrIiNTLuQv2Y6JyuB4Mqbt11sW04obKbuefhEfWTn7oj+ReB+uJMvJ0V5tb7KfV+rfyX5gpeSum/941Ohxsqx6vw/vvr8+GiHSm36Xjvvs4dHU6eQt3tcr8bH763lw1UadsuCq+Y2zr3c6dXK/SvmI9/k4j1o/hlgJTLJM/lBWiVO9+a32d8b2s0o/ytm/CxP529LI3zf2gQAEIAABCEBApHQSWO3g3h5kxQ84Kh2Q63s/+sF5o84dqrSdPYCJH4DpAY8egFt5qPYLfCUJjI++2fOWKr0XP3CNHwi7ml6VjK2ekatqX6Qs5G88GUxzMFmtXRNpb60/KPaAOM3Iar39wm6nMcT7ZKWRo/G+Iypzep5WpamBlfp8tXNHK/X9Wt8HjT0pCxrrAw89MerHmkrxjzcSWOm7kQVXbV/8h6U0+a73706jI4FxKdbnlWZMVJPASucG1zu6mDYvtb5DfA4BCEAAAhCAgBsCpZLAagdd1QTJbl9rBM0edNqD63qnntU66I3XW23b5EFftbrH+8U9TVep9Ot88uC2VnnxEaRaoxS1yqr387QL2Yy3MExSpOqNodZ2jUxVrrdfVOvjtp3x8y9rHcCP17+TU/iqSUm1flRJgOv9PsX5JuNIKxuuudaqv5EFoJJ/d7T9jUpgJe7JHI03EpgcFa93FG+8H1Vc/XBV63vH5xCAAAQgAAEIjCUQvAQmm1RJOsY7cNKDwWrTopJl2wPpeg9e6pXAWgeQ8QUUspbASjFXGgEa78vUiATGR7Js2ZWm61art1EJrPQDgI2l1o8Daf+gpJXANP2iXgmsR7jG69/KWW/2nMW8JLDSdNJ4/0jDKs229XKNT69Nu+CR7UeVvgP6Wbw8lxKY/IEACUz7jWZ7CEAAAhCAQLgEgpfA5OhcpXNwKp1vE0+Z/ZU8fiAXl8nkwVK9v4LXK4GVRmuSB4b2QDBLCaw1OpZ2sZ0Qp4Na7o2cu1frz0DalUfT9It6ZWW8Mm381S6/Yj+Pf8eylsBqK6kmFxtJI3ZZcFU2lc6z1PdrTVOu9++OlhWaBKbJS63vD59DAAIQgAAEIOCOQGkkUJHYg7vkCFJy9KIavnoPpMs6EjjeAeZEpCgLGcxyYRh7sK3nQtU6gE/7VUzDMc0BdL19N81IYD1tz1oCq/2QMhEJzIJrpX5Q7e9Rctt6c+daApP1ZjESyMIwaf9CsD0EIAABCEAgHwKlkkBFZoUjPmpVr7RVG+FLjhzUOoi05/9UWwK90v71nqOU1UhgrUUsKi2kkbaLupDBicqfjbnWAjBpZC0NB8txvFFV5aSrWuqt3n6RRiTGu37lFz+3yIxofVj1wuvJtqaRwDTfB61nPIFISuB45/dWYuOaq/anj807YMzF6OsZuav37479+/aL1a+lWrl2Im1Nk99q3wMuEZHmLwTbQgACEIAABPIhUDoJtL+W6yhO/JyuSquD2oMqXTr9zGUnVrx2XHyKZPLcHK0jPm1UD+biB2jVDu7GWx200oqf8XqzksBakldLEtN0Vytyac+dsvu5WGym1nUCNbf1xGfLSbPiZ5qLxVc6P3G8lS2TbCpNfaxUZnIRkmoXDtft9NIkVlLTSEKa74PtT9VW8NXLMCRH/KvJxnirg9b6vtUr15ZfXO7r/c5U6ovV/u7UM503+V2stLBTvX0oTX6r/Q2o1L56f5hL83eFbSEAAQhAAAIQqJ9AKSVQm2+lL35wXu1cp/i0t+R5cfFrlSWloNJCDskpdPHzEV1fJ1DF1d4mujpoPSNfLq4ZWH/XzHbLWuc/1jMVUiNsRAJ1v2rnj1UaIYxfFsRSScpevbJi969UZrJ/V2NUz0Il1YQvzfchLoLJ6x7qpViSFyBPMrXf/WpsXHOt9Peg3vNo0/zdSf4dq/WjiEqg/u1JXv+0nj7kQgI1j/WwzvYbT+kQgAAEIAABCMQJBCuBpBECEIAABCAAAQhAAAIQgAAE0hNAAtMzYw8IQAACEIAABCAAAQhAAALBEkACg00dgUMAAhCAAAQgAAEIQAACEEhPAAlMz4w9IAABCEAAAhCAAAQgAAEIBEsACQw2dQQOAQhAAAIQgAAEIAABCEAgPQEkMD0z9oAABCAAAQhAAAIQgAAEIBAsASQw2NQROAQgAAEIQAACEIAABCAAgfQEkMD0zNgDAhCAAAQgAAEIQAACEIBAsASQwGBTR+AQgAAEIAABCEAAAhCAAATSE0AC0zNjDwhAAAIQgAAEIAABCEAAAsESQAKDTR2BQwACEIAABCAAAQhAAAIQSE8ACUzPjD0gAAEIQAACEIAABCAAAQgESwAJDDZ1BA4BCEAAAhCAAAQgAAEIQCA9ASQwPTP2gAAEIAABCEAAAhCAAAQgECwBJDDY1BE4BCAAAQhAAAIQgAAEIACB9ASQwPTM2AMCEIAABCAAAQhAAAIQgECwBJDAYFNH4BCAAAQgAAEIQAACEIAABNITQALTM2MPCEAAAhCAAAQgAAEIQAACwRJAAoNNHYFDAAIQgAAEIAABCEAAAhBITwAJTM+MPSAAAQhAAAIQgAAEIAABCARLAAkMNnUEDgEIQAACEIAABCAAAQhAID0BJDA9M/aAAAQgAAEIQAACEIAABCAQLAEkMNjUETgEIAABCEAAAhCAAAQgAIH0BJDA9MzYAwIQgAAEIAABCEAAAhCAQLAEkMBgU0fgEIAABCAAAQhAAAIQgAAE0hNAAtMzYw8IQAACEIAABCAAAQhAAALBEkACg00dgUMAAhCAAAQgAAEIQAACEEhPAAlMz4w9IAABCEAAAhCAAAQgAAEIBEsACQw2dQQOAQhAAAIQgAAEIAABCEAgPQEkMD0z9oAABCAAAQhAAAIQgAAEIBAsASQw2NQROAQgAAEIQAACEIAABCAAgfQEkMD0zNgDAhCAAAQgAAEIQAACEIBAsASQwGBTR+AQgAAEIAABCEAAAhCAAATSE0AC0zNjDwhAAAIQgAAEIAABCEAAAsESQAKDTR2BQwACEIAABCAAAQhAAAIQSE8ACUzPjD0gAAEIQAACEIAABCAAAQgESwAJDDZ1BA4BCEAAAhCAAAQgAAEIQCA9ASQwPTP2gAAEIAABCEAAAhCAAAQgECwBJDDY1BE4BCAAAQhAAAIQgAAEIACB9ASQwPTM2AMCEIAABCAAAQhAAAIQgECwBJDAYFNH4BCAAAQgAAEIQAACEIAABNITQALTM2MPCEAAAhCAAAQgAAEIQAACwRJAAoNNHYFDAAIQgAAEIAABCEAAAhBITwAJTM+MPSAAAQhAAAIQgAAEIAABCARLAAkMNnUEDgEIQAACEIAABCAAAQhAID0BJDA9M/aAAAQgAAEIQAACEIAABCAQLAEkMNjUETgEIAABCEAAAhCAAAQgAIH0BJDA9MzYAwIQgAAEIAABCEAAAhCAQLAEkMBgU0fgEIAABCAAAQhAAAIQgAAE0hNAAtMzYw8IQAACEIAABCAAAQhAAALBEkACHaTuvY29DkqhiJAJTO5ql/a2Ftm6YyDkZhC7AwJdk9qkp7NNNm3b5aA0igiZwKT2Vpk2uUM2bO0PuRnE7oBAW2uLzJ7eKes29zkojSJCJtDSIjJ3v255fxPHjvNmdYecyuBjRwIdpBAJdAAx8CKQwMAT6DB8JNAhzMCLQgIDT6DD8JFAhzADLwoJ3JdAJLDYzowEOuCPBDqAGHgRSGDgCXQYPhLoEGbgRSGBgSfQYfhIoEOYgReFBCKBvnRhJNBBJpBABxADLwIJDDyBDsNHAh3CDLwoJDDwBDoMHwl0CDPwopBAJNCXLowEmkycvnyFvPDKm6Nycu7pJ8vlF5wRvXfpVbfII6uejp4vWniI3LfyqlHbIoG+dOfi4kACi2PvW81IoG8ZKS4eJLA49r7VjAT6lpHi4kECkcDiet/ompHAYQlcfPTCEemLI7rngcfktjsflFX33xS9rcKY3BYJ9KU7FxcHElgce99qRgJ9y0hx8SCBxbH3rWYk0LeMFBcPEogEFtf7kMAx7CuJnd0o+VlSCnU7JNCX7lxcHEhgcex9qxkJ9C0jxcWDBBbH3reakUDfMlJcPEggElhc70MCK0pgfDpofCro0tMulvPPOVXOXHZitN8TT62W5VdcLy8+fsdIOUigL925uDiQwOLY+1YzEuhbRoqLBwksjr1vNSOBvmWkuHiQwGIk8Lpb75WfPfrkyMy+4nrAxGq2HvLQ3dfIgvlzJlQY00ET+CzclddeJkuOO0qOPOFcufKSs8dIoIX/tZ98TY6Zc5wsP+a7Mr1z+oSSwc7hEugw1wNrNdeB6t+1J9xGELkTAu1trdLR3iK9/fQFJ0AnUMjg4N7oe1nUTQ/8O811I3f27S4qBOr1hECLOfLv6WqTHb30BU9SUlwY5k/SlK4O2d7LdYWn9nQ4y8OatevklLO+P6Y8e7xepATqgNIxRx0qN6y4cMLtRQInjHD8AuJTQGuNBLasGDrAmNY5Tb5zzH+WCxdfKjO7Z2YcIcX7RkB/8ddjzb6BQd9CI56cCbS3tUiHEcFefhDImfzY6ryQwI5W2ckPAoX3haID0H8fejrbZTs/CBSdisLr16PGKd3tso0fBGSq4eDipoJ3x30Pjxq00XJVmG7+y/ujBR2LlEAXbbRlIIEuaVYoKy6Btc4JfGLNE3LlYytk1Zp/iErq7uiRbx35bfnuMf9F9u+Z2DBtxs2keIcEmA7qEGbgRTEdNPAEOgyf6aAOYQZeFNNBA0+gw/CZDroPpouLxdsRwPisvUrpshKop3hdfeNdI5vET+/SN61Q2g0qzQbU9+Jl6OzBd977cNR78ema6hLz5s4aNRKYvDKBrUfXHomXrXHET1NDAh1+GbXz3PDnPx1JjIVvk1fv6qAvrH9Wrn/qR/Lom38ve81/nW2d8vuHny0XHXu5zJ/6MYcRU5SPBJBAH7NSTExIYDHcfawVCfQxK8XEhAQWw93HWpFAtxJopS0pc8nc2+2+snTxyDG/ipje7KXfkqOFVjDtKWJWwOKXi7OXkYuXq++998HGkXKTEqizDOfuP3Pk87iLqHe8+/76kSsWJE9TQwIdf6v1vL/4zSbbvpfmOoGvbnpZbvi3/0d+9uv7ZXDvoLS3tstpC8+U7x17hRwy45OOI6c4Xwgggb5kovg4kMDic+BLBEigL5koPg4ksPgc+BIBEuhWApPCVS3PlaaDJgd61AeSC65o+XrT8/kqCZgdPIpLaKUBJDsSWGmByVp9U6Xx6ycdH4khEliLVs6fV1od9M0tr8tNT18jf/vqvbJ7cLc5X6xVvvrJZWZk8P+SRfsfnXOEVJc1ASQwa8LhlI8EhpOrrCNFArMmHE75SGA4uco6UiTQLwnUqZcqcFauKuXfjvylkUBbrpYXHwms59zE5JRULcNOCUUCs/6Gpix/vEtErN32jtz8zHXyNy/fJf17+qOST1hwklxoZPD4+V9KWROb+0oACfQ1M/nHhQTmz9zXGpFAXzOTf1xIYP7Mfa0RCXQrgfVIldZYbSQwKYHjXXohDwm05wrGRxYZCfT122ziquc6get2vC9/8ewtctcLt8v2gW1Ra46d+3m5aPHlctLBX/W4dYRWDwEksB5KzbENEtgcea6nlUhgPZSaYxsksDnyXE8rkUC3EmjFrNrCMCp/Oo2ylgRqVMnLwiXz6UICK00fjddTKQYksJ5vVkHb1COBNrSt/VvkjtUr5fbn/lQ29W6I3j5s1pHy3WMvk39/6DelraWtoFZQ7UQIIIEToVeufZHAcuVzIq1BAidCr1z7IoHlyudEWoMEupVALc2u3ZEUQRUqO5WzHgm05cRHA1X8/vvfrUp9TmC16aAab/K6gfGFYeLCp9taaWQ66ES+dRnum0YCbRi9u3vl3pf+Wm795Q2iU0b1tmDawXLBMZfK7x92tnS1d2UYMUW7JoAEuiYabnlIYLi5cx05EuiaaLjlIYHh5s515EigewnUEiud0xdfsbMeCdRyKp2Pl1wdNC6J1RaGGU8CrQhu2LR1BIYV2ORF77UNv1j9GgvDuP4iuiqvEQm0deuiMQ+89jfyp89cL7qyqN5m9xwgf3jU+fIHn14uMzr3cxUm5WRIAAnMEG5gRSOBgSUsw3CRwAzhBlY0EhhYwjIMFwnMRgIzTFlpi27Za26lbV1ODZuIBNoQ9dqCeo3Bm5++Vn6x7t+it3s6JssZh58jyz97sRw0dUFOraGaRggggY1QK+c+SGA589pIq5DARqiVcx8ksJx5baRVSCAS2Ei/yWIfJNABVRcSGA/jX9c+EU0Tfeyth6MLz+t5gl/71Dfku8f8Fy4v4SBfWRSBBGZBNcwykcAw85ZF1EhgFlTDLBMJDDNvWUSNBCKBWfSrRspEAhuhltjHtQTa4l/f8qrc+osb5f5X7hm5vMQX5y+V84+5RH7r47/jIHKKcEUACXRFMvxykMDwc+iqBUigK5Lhl4MEhp9DVy1AApFAV31pouUggRMlaPbPSgJtaBt710erid65+i9kS//m6O3fmHm4nPeZ78rvHXaWdLZ1OmgFRUyEABI4EXrl2hcJLFc+J9IaJHAi9Mq1LxJYrnxOpDVIIBI4kf7jcl8k0AHNrCXQhtg7sFN+8tId8hfP3SzvfPR29PbM7tly9qLz5NxF35EDJs910BqKaIQAEtgItXLugwSWM6+NtAoJbIRaOfdBAsuZ10ZahQQigY30myz2QQIdUM1LAm2oe/bukb9//X+Yi8/fLM988PPo7Ultk+RUc53B880iMofPWuSgVRSRhgASmIZWubdFAsud3zStQwLT0Cr3tkhgufObpnVIIBKYpr9kuS0S6IBu3hIYD/m5D38ht/3yRiOFD4hebkJvx8//kvzR0RfJVw75mrSY/7hlTwAJzJ5xKDUggaFkKvs4kcDsGYdSAxIYSqayjxMJRAKz72X11YAE1sdp3K2KlEAb2Hvb18p/fe7PzHTRv5Kt/Vuitz8+/RD5AzNN9Iwj/kCmd85w0FKKqEYACaRvWAJIIH3BEkAC6QuWABJIX7AEkEAk0JdvAxLoIBM+SKBthp43eO/Ld8pfPv+n8uaW16O3u9q75bSFZ8gffvoCpoo6yHelIpDAjMAGWCwSGGDSMgoZCcwIbIDFIoEBJi2jkJFAJDCjrpW6WCQwNbKxO/gkgTY6vb7gP739iNyxemX0OLh3MProcwceH8ngKZ/899LR2uGg9RShBJBA+oElgATSFywBJJC+YAkggfQFSwAJbC4JvOeBx+TqG++Sc08/WS6/4AyvvghIoIN0+CiB8Wa9u22N/NXzt0UjhFv6NkUf7d8zR844/Bz51qJvy0FTFzig0NxFIIHNnf9465FA+gISSB9IEkAC6RNI4Ng+MG9Wd+k7xunLV8jioxfKzx59Ulbdf5NX7UUCHaTDdwm0Tezb3Sf/47W/kTueXynPr/9l9LYuHPPlBSfKt478drSQTHtruwMizVcEEth8Oa/WYiSQvoAE0geQQPpANQKMBO4j0wwSeOQJ58qLj98h+rjy2stkyXFHefPlQAIdpCIUCYw39Xmzquh/e/G/yt++ep/sHNgRfaTXGdTRQb3u4LwpBzkg0zxFIIHNk+taLUUCaxFqns+ZDto8ua7VUkYCaxFqns+RwGwk8LX+XbJ2YCD3jnRo5ySZ31H59CqdCvrUL1+WG1ZcKJdedUsUmz735YYEOshEiBJom71j13a5/7X75L+98Jfywvpno7d1dPA3539Zvnn4t+TrnzpNus3CMtzGJ4AE0kMsASSQvmAJIIH0BUsACaQvWAJIYDYSeNG7a+WW9Rtz72g3HzRPLtx/dsV6dSroRd8+LRr9e+Kp1fLDH9/u1ZRQJNBBdwlZAuPNrzQ6OLljinz1k8vk940Q6vUHue5g5Q6DBDr4IpWkCCSwJIl00Awk0AHEkhSBBJYkkQ6agQRmI4E3rd8gD2zZ6iBD6Yr4nhHAb8yYPmanNWvXydkX/cko6Vt62sXyox+c582UUCQwXa4rbl0WCbSN0+mhP/v1/fLTX90tT679Z7PO6N7oI11ARkcHv7nwrOgahNz2EUAC6Q2WABJIX7AEkED6giWABNIXLAEkMBsJ9K2HXXfrvXLHfQ+PCesrSxd7MyUUCXTQa8omgXEk721/V/7m5bsiIXxr6xvRRzoa+Pn5X5TTDztbvvbJb8jkSVMcUAy7CCQw7Py5jB4JdEkz7LKQwLDz5zJ6JNAlzbDLQgKbQwJ1KuiyU5bImctOHGmwTgldfsX10UIxPtyQQAdZKLMExvE888HPzeqiP5X7X7lPNvcNzbvWC9Gfeujvycmf+F358kG/Jd0dPQ6IhlcEEhhezrKKGAnMimx45SKB4eUsq4iRwKzIhlcuElh+CbSy99Dd18iC+XNGdVKdEnr+OaeOksOiejES6IB8s0igRbVn7x7553f+Uf7WyOBDbzwoOwa2Rx91tnXKbx60VE46+Kty0iFfbaoVRpFAB1+kkhSBBJYkkQ6agQQ6gFiSIpDAkiTSQTOQwPJLoINukksRSKADzM0mgXFkeu3Bf3z7f8k/rXlEHn7jf8qm3g0jHx8260hZ+rHfli8aMdTVRss8SogEOvgilaQIJLAkiXTQDCTQAcSSFIEEliSRDpqBBCKBDrqRkyKQQAcYm1kC4/h0AZnnzPUHH3vrYfmHN/9eVptLTthFZSa1TZJj535evvSx34runzngGGlraXNA348ikEA/8uBDFEigD1nwIwYk0I88+BAFEuhDFvyIAQlEAv3oiWaNj73m5kswocaBBFbO3Mbe9Wba6D/JP79r7mb66Npt74xsOHXSNPncgcfL5+d9Ub4wb4l8Zs6x0tFa+WKbIfQLJDCELOUTIxKYD+cQakECQ8hSPjEigflwDqEWJBAJ9KWfIoEOMoEE1gfx15tfkX95d5U8YaRw1ZrHRs4ltHt/cf5Ss+rokuh6hMfMOc4sOtNVX8EebIUEepAET0JAAj1JhAdhIIEeJMGTEJBATxLhQRhIIBLoQTeMQvBOAnXVnA2bKl/s0ZclVZPJQwLTd2ddXObZdU/L42v+wdwflefWPSP6nr3pIjM6fXTofMKl8kRybLQAACAASURBVNm5n/N6pBAJTN8HyroHEljWzKZvFxKYnllZ90ACy5rZ9O1CApHA9L0mmz28kkC9psa8ubO8uYhivciRwHpJVd/uo11b5Yl3Ho8WmfnX956QN7f8etTGeimK4w78TbP66Jfl+HlfksUHfmHilTosAQl0CDPwopDAwBPoMHwk0CHMwItCAgNPoMPwkUAk0GF3mlBRXkngkSecKyuvvUyWHHfUhBqV985IoHvieh3CZ97/uTyz7qno8ZcfPi07B3aMVNTe2i4LzeqjRx9wbLTIjN4Pm7VI9P0ibkhgEdT9rBMJ9DMvRUSFBBZB3c86kUA/81JEVEggElhEv6tUJxKYoKLTUefuP1PuW3nVyCeXXnWLPLLq6ej1ooWHjPpM30MC8+nOL5jVRn+x7t/kFx/8mzxvViF9ZdNLYyr+7JzFcszc4+SIWUfJov0/Y+5H5xIcEpgL5iAqQQKDSFMuQSKBuWAOohIkMIg05RIkEogE5tLR6qjEKwnU6aDLTlkiZy47sY7Q3W+iAqi3uATe88BjctudD8qq+2+KPtMYFx+9UC6/4IyRAJBA97mop8TegZ3RZSj0shTPffiMPGvub215feSyFFqGXpriN2YeEQnhUUYIF83+jBwx+yjp6ZhcTxV1b4ME1o2q9BsigaVPcd0NRALrRlX6DZHA0qe47gYigUhg3Z0l4w29ksAnnlotP/zx7SPClXHbRxVvBfTd99fL08++MjLal5S+pBRqIUhgnpkav65tuz6KFpx5YcNz8uL65+WF9c/Ja5t/NWanQ2Z8Sg4zcqhieNjsRXK4mUq6YNrBDTcECWwYXel2RAJLl9KGG4QENoyudDsigaVLacMNQgKbQwLVF66+8a5R/aTSbMKGO5KDHb2SQD0ncLxbVquDxkXvulvvHSWBOjp4/jmnjoxOqqguv+J6iceCBDroiRkW0bu7V3618cUhMfzw2ejxZfO6z7wfv+m1Cw8z5xnqXUcLVQx1WunkSVNqRocE1kTUNBsggU2T6poNRQJrImqaDZDApkl1zYYigc0lgXFfUM+58pKzC5vxmOycXklgzW9OBhvo+X56u2HFhdFjUgKTCbMS+NDd18iC+XMyiIgi8yAwuHdQXtv0mjy/7vnovnrd6ujxrS1vjZpOqrEcOOVAOXTWoXLozEPlUzM/NfJc3+s2q5ZygwAE/CTQP7BHOjva/AyOqCAAAQhAoLQE7EhgXAIrnVJWJADvJNBKVhxKliuGVrsu4eyZ06NpqYwEFtk9869bVyB9acNqeWnjC/JrM41Up5PqKOLW/i0Vg5k7+UA5eMYnI0FcOHuhHDzlMLNq6REyb8pB+QdPjV4QYCTQizR4EQQjgV6kwYsgGAn0Ig1eBMFI4L40zJvl7od0/WF/7Udrc8+xDgjMnzp/TL1JCVyzdp2cctb3xadBJK8k0AKLA7LQ8ho+TY4Eck5g7t8nLyv8cMcH8urml+WVjS+bEUTzaJ6/ap5v6d9cMV6dWqoyqOcdfmLGoZEUzpv6MXM/SA6cPM/LNhKUGwJIoBuOZSgFCSxDFt20AQl0w7EMpSCB2UjgRQ9dJLc8NTS7L8/bzafcLBceNzSbMH6rdE6gfo4EVslOctTNblZpMZasEpyUQFYHzYp0Ocpdt+P9aOGZN7YaKdz0K1n94QtGFF8SXaBmvNvHpn1cDpq6wCxGc4gcPP0TI/dP7HeoTOmYWg44TdoKJLBJE1+h2UggfcESQALpC5YAEpiNBN7085vkgV89kHtH+97nvyffOOwbVSUwPh3Ut1PKvBoJrHax+EqLsWSV5aQEaj1cJzAr2uUpN7kwzHvb341k8FebXpQ1W9+Steb1++a+dts7VaeWWhr7dc2SQ8wUU51qun/PHJnTMzd6PGDyXDkgejww+oybnwSQQD/zUkRUSGAR1P2sEwn0My9FRIUEZiOBReRyvDornROo21dznSLi90oCfRgJbCQJrA7aCLVy7ZNmdVBdrfSdj94SFcU3zXUN12x7S97Y/Jq8tfUNcx7iK3WD0WmlB5kRxQOnzI8ubzHXvJ5vppzOH556ul/XzLrLYkN3BJBAdyxDLwkJDD2D7uJHAt2xDL0kJLB5JbCaGBbVp72SQB/OCWwkEUhgI9TKtU8aCazVch0tfPujN0Wnmq7fsU7W7fxA1u9cJx/q3ZybqI+b+zbWKib6XAVxVvf+5j5bZnWZe89s0ZHGoef7y/7dB0QjjBO5RmJdgTTRRkhgEyW7RlORQPqCJYAE0hcsASSwuSQw2fM5J3CcvwV5rw7q4s8SEuiCYthluJTAeknoKOL7O9bKux+9HU0z1SmnOrr47kdr5L0d70rvwM56i4q2m2FGDnWa6Zzhu31+gJmOqtNQZ3TtJzM695PZPQekKrfZNkYCmy3j1duLBNIXkED6QJIAEtgcEhhCz/dqJDAEYJViRAJDzZy7uIuQwFrR6+I0m3o3ysa+DbJp5/rocWPvhuicRB1N1OebetebcxXXygdm1DHNTRev2a97ZiSOKoZWEGeaEceZZqRxthl93H/ynH3PjUQ2yw0JbJZM124nElibUbNswUhgs2S6djuRQCSwdi/JZwsk0AFnJNABxMCL8FEC0yLV6afR3Uw/VSlcb6ae6qPeN+78MLocxpY+c69yWYxa9akwqiDO7DZ386jTUvc3o4o97T3S1d49fO+SrjbzvMPc27pkcscUmdY53dxnyHRz7zbb+X5DAn3PUH7xIYH5sfa9JiTQ9wzlFx8SiATm19vGrwkJdJAJJNABxMCLKIMEpkmBnpO4eVgIPzIji5vM60gQ+zbJBjO6uDG6bzDnMOqI4/qGxbFSTDoddfokI4VdM0SvxxhJ4qTpQ8/Nfap5HX9/ihHJKeb9KZOGHrMWSSQwTU8q97ZIYLnzm6Z1SGAaWuXeFglEAn3p4V5IoC6XqheDv/rGu8blEr/Whi8ANQ4k0KdsFBNLs0lgI5R1cRuVxc29m6KpqZvNVFV9vWNgu/Tv6ZO+gV7p1UezeurQvc98ts1MX90qH+3aKhvMaKSrm45KTp00dUQWO1onyaQ2czePHfrY1jnqdad5rfepRjZ16quuvKp3Ow1WF9+xNyTQVZbCLwcJDD+HrlqABLoiGX45SCAS6Esv9kICfYHRaBxIYKPkyrMfEphPLnWxm61GCLeaKakfqRyau577+JG5bzPv7xjYEY1GbuvX9/T19kgit0efb6t7VdVGWqOjjyqFs8x01xnd06VDuqVn0mQz3dXcO3rM1NapZoprpxmJNNNfzWsdkdTn+tjTMVk627si0eyOTYfVabL6GbcwCSCBYeYti6iRwCyohlkmEogE+tJzvZLAahdQ1EtH3Hbng7Lq/pt84TYqDiTQy7TkGhQSmCvuCVWmI4wqhduNIEaPu4ZGIncN7pJde/TeLwPmsd+81kd9re/3m0eVT3tepC6wE02LNSObjZ4nmaYhKplWFDvN+ZJdw9Ko77W1tEt7a7t5bIvu0fNW81wf9T/zOPR569B7up2+Z/az27Wb96IRUB0Vbe+UjtaO4dFQ8zwaHbWjpcOvzefR+8Ojp0PPh97TkVWNr9lvSGCz94B97UcC6QuWABKIBPrybQhCAu1lI5gO6ku3IY4kASSQPrHZjECqIPbu2Sq7pVfe27pJdu7eITt3mfvuneaSHTuk10xz7dXn+mhGNfV5nxHQ6NFMh9Xn8emwfXuGpsWGflMhVFG18qniqdLZGhfW8T4fltaR/Y3MqtRG+8dENhJgLVfLH/7cSrGtq1VazH7m/5EQawytov/pa32u29vPdduRbYb3ibYx/41sY/czZbUMlx2Vae6d7R0ypbtDtu3YM1RPrIyheodiMZWPfD6ynUY1vM+ofRN12Ljj7dA4kHC/vjVIoF/5KDIaJBAJLLL/xesOQgKvu/Ve+dmjTzIS6EuvIY4xBJBAOoUlkMU5gSqCOlrZHz32G1nsHXmur/cM7pbd5q6Pe2Qwej44uEd27x16f3Cvea6f7R0c2TZ6z3y+x2yn+w3K3khAB8wIaP9uMxo6ODA0CqojpMOvtS7dT7fbZT7XkVLdXkdKh7YfGj3Vabjc/CagU42tOI4nqHHJbBmWVaOslQU1Jrbx/VSsVUxHCasKsX2vmtjG6hkTR4Xy9sn76Pgi4TcxFHlrbRHp6WyX7X27iwyDuocJmHSM/VEm8YNQ/EcY80vNyI82lb8vFX7kGf4+JPu97j97Wrds3Nbf9Pn4vU//TtMzKBJA4RJY6eLwlYCsvPYyWXLcUUWyqlo300G9TEuuQSGBueL2urIsJNDrBtcRnI56RsJpBDIS1UhQja4Oy6kV1pHPjayOltfh/XT/qIwhqU2KbPTZOJ+r6OrneyMZ1voHh+6qwPpeFNPQe/o62n54u+i1bqOabT+37+k2Zlu7TbSdKUvMuzoC1Ld7YHS5w2Unyxl6vXcknpH44vWZekbFGGuDvq+Czg0CEIBACAT2XrU3hDBLG2PhEhgnW+2cQN/pI4G+Zyj7+JDA7BmHUgMSGEqmso/T13MCdTrykOQOS3BCaofEeLTUjojn8D6RrMaFOSa2Y+R2ZJ999dlt9pohmbjsjtq3yn5DsY2NL9mm0WUVe7DJSGD237c0NVT8ISX2g9CYH2CGf3wxva5Cv4/9cGR/lLHfjdj28e9QR1uLmWExmCbkUm77f87736VsVyiN8koCQ4GWjBMJDDVz7uJGAt2xDL0kJDD0DLqL31cJdNdCSqqXAOcE1kuq/NtxTuC+HM+b1V3+hHvcQiTQQXKQQAcQAy8CCQw8gQ7DRwIdwgy8KCQw8AQ6DB8JdAgz8KKQQCTQly7slQSuWbtOTjnr+1XZsDqoL92GOJIEkED6hCWABNIXLAEkkL5gCSCB9AVLAAlEAn35NnglgUtPu1i+ftLxcvyxR8oPf3z7yGqgpy9fIctOWSJnLjvRF26j4mAk0Mu05BoUEpgrbq8rQwK9Tk+uwSGBueL2ujIk0Ov05BocEogE5trhxqnMKwm0C8MsmH+AnH3Rn4xIoK4gGpdCX+DZOJBA3zKSfzxIYP7Mfa0RCfQ1M/nHhQTmz9zXGpFAXzOTf1xIIBKYf6+rXKOXEqiXglAhtNM/uVi8L92FOKoRQALpG5YAEkhfsASQQPqCJYAE0hcsASQQCfTl2+CVBOq0z8VHL5TLLzhD4s+5WLwv3YU4kED6QC0CSGAtQs3zORLYPLmu1VIksBah5vkcCUQCfentXklgEoqOBtrbQ3dfIwvmz/GF26g4mA7qZVpyDYqRwFxxe10ZEuh1enINDgnMFbfXlSGBXqcn1+CQQCQw1w43TmVeS6AvkGrFgQTWIlT+z5HA8ue43hYigfWSKv92SGD5c1xvC5HAekmVfzskEAn0pZd7JYF2YRg9JzCkGxIYUrayiRUJzIZriKUigSFmLZuYkcBsuIZYKhIYYtayiRkJRAKz6VnpS0UC0zMbswcS6ABi4EUggYEn0GH4SKBDmIEXhQQGnkCH4SOBDmEGXhQSiAT60oW9kkDfrwdYLWlIoC/dubg4kMDi2PtWMxLoW0aKiwcJLI69bzUjgb5lpLh4kEAksLjeN7pmryRwzdp1o64P6AukWnEggbUIlf9zJLD8Oa63hUhgvaTKvx0SWP4c19tCJLBeUuXfDglEAn3p5V5JYHw10EqA7HUDfYFn40ACfctI/vEggfkz97VGJNDXzOQfFxKYP3Nfa0QCfc1M/nEhgUhg/r2uco3eSeDKay+T5MIw9zzwmNx254Oy6v6bfOE2Kg4k0Mu05BoUEpgrbq8rQwK9Tk+uwSGBueL2ujIk0Ov05BocEogE5trhxqksCAl84qnVsvyK64WRQF+6DXEkCSCB9AlLAAmkL1gCSCB9wRJAAukLlgASiAT68m0IQgKvu/Ve+dmjTzIS6EuvIY4xBJBAOgUSSB9IEkAC6RNIIH0gSQAJRAJ9+VYULoF2lK8WkErTRGvtk9fnTAfNi7S/9SCB/uYm78gYCcybuL/1IYH+5ibvyBgJzJu4v/UhgUigL72zcAmMgyjqYvF6aYoXXnlzJJSkcF561S3yyKqno88XLTxE7lt51aj8IYG+dOfi4kACi2PvW81IoG8ZKS4eJLA49r7VjAT6lpHi4kECkcDiet/omr2SwKKgqARasbMjkw/dfY0smD9HkovS6LaLj14ol19wxki4SGBRmfOnXiTQn1wUHQkSWHQG/KkfCfQnF0VHggQWnQF/6kcCkUBfeiMSWCET8RHJpPRVWqkUCfSlOxcXBxJYHHvfakYCfctIcfEggcWx961mJNC3jBQXDxKIBBbX+zwfCVx62sWyYdPWinzyWB1UJe/qG+8SOxKo8Zx/zqly5rITo5gqrVSKBPrSnYuLAwksjr1vNSOBvmWkuHiQwOLY+1YzEuhbRoqLBwlEAovrfR5LoI66zZs7S25YcWHufOIL1MTPCdRRwSsvOXuMBFpJ1EC37RzIPV4q9ItAR3urtLa2SP+uPX4FRjS5E2hva5WO9hbp7acv5A4/UeHg4N7oe1nUTQ/8Oye1yc6+3UWFQL2eEGgxR/49XW2yo5e+4ElKigvD/Ema0tUh23s5dpza01FcHqhZvJoOWtTCMPF+sGbtOjnlrO+LFcF6RgK38Ue96b9K+ou/Hmv2DQw2PYtmB9De1iIdRgR7+UGg8K7ghQR2tMpOfhAovC8UHYD++9DT2S7b+UGg6FQUXr/+LDWlu104dhSZajhwK44AEliBfVz8OCewuM4ZUs1MBw0pW9nGynTQbPmGVDrTQUPKVraxMh00W74hlc500H3ZmjerO6TUlS5WryRQhWvZKUtGpl7mQVungT75zIsjq33aaaF2JJDVQfPIQvh1IIHh59BVC5BAVyTDLwcJDD+HrlqABLoiGX45SCAS6Esv9koCVcB++OPbZdX9N+XKR6ehxm9cJzBX/KWoDAksRRqdNAIJdIKxFIUggaVIo5NGIIFOMJaiECQQCfSlI3slgUkZS0LKY3XQRhLD6qCNUCvXPkhgufI5kdYggROhV659kcBy5XMirUECJ0KvXPsigUigLz3aKwn0BUraOJDAtMTKtz0SWL6cNtoiJLBRcuXbDwksX04bbRES2Ci58u2HBCKBvvRqJNBBJpBABxADLwIJDDyBDsNHAh3CDLwoJDDwBDoMHwl0CDPwopBAJNCXLowEOsgEEugAYuBFIIGBJ9Bh+EigQ5iBF4UEBp5Ah+EjgQ5hBl4UEogE+tKFvZNAXSH0hVfejPjYBVr0XMGvLF1cyEXk60kUElgPpXJvgwSWO79pWocEpqFV7m2RwHLnN03rkMA0tMq9LRKIBPrSw72SQBXAeXNnRbKn1+r70Q/OkyXHHSXJyzT4As/GgQT6lpH840EC82fua41IoK+ZyT8uJDB/5r7WiAT6mpn840ICkcD8e13lGr2SQB3xe+jua2TB/DmjJNBeu4/VQX3pNsSRJIAE0icsASSQvmAJIIH0BUsACaQvWAJIIBLoy7fBKwnU0b+7bv7jMRLISKAv3YU4qhFAAukbSCB9IEkACaRPIIH0gSQBJBAJ9OVb4ZUEXnfrvfKzR5+MLhZvp4MumH+AnHLW9+Xc00+Wyy84wxduo+JgOqiXack1KCQwV9xeV8ZIoNfpyTU4JDBX3F5Xxkig1+nJNTgkEAnMtcONU5lXEqhx2qmf8ZivvORsOXPZib4wGxMHEuhtanILDAnMDbX3FSGB3qcotwCRwNxQe18REuh9inILEAlEAnPrbDUq8k4CfQGTJg4kMA2tcm6LBJYzr420CglshFo590ECy5nXRlqFBDZCrZz7IIFIoC892ysJvPSqW+SRVU9LcgEYLhHhS3chjmoEkED6hiWABNIXLAEkkL5gCSCB9AVLAAlEAn35NnglgXoe4PnnnDpm6icLw/jSXYgDCaQP1CKABNYi1DyfI4HNk+taLUUCaxFqns+RQCTQl97ulQTqiJ+9QHwcEJeI8KW7EAcSSB+oRQAJrEWoeT5HApsn17VaigTWItQ8nyOBSKAvvd0rCWQk0JduQRxpCTAdNC2x8m6PBJY3t2lbhgSmJVbe7ZHA8uY2bcuQQCQwbZ/JanuvJFCnfV59410jF4zXRq9Zuy66RITPK4SyMExW3TOccpHAcHKVdaRIYNaEwykfCQwnV1lHigRmTTic8pFAJNCX3uqVBCqUSpeIqDRF1BeAGgcS6FM2iokFCSyGu4+1IoE+ZqWYmJDAYrj7WCsS6GNWiokJCUQCi+l5Y2v1TgJ9AZMmDiQwDa1ybosEljOvjbQKCWyEWjn3QQLLmddGWoUENkKtnPsggUigLz0bCXSQCSTQAcTAi0ACA0+gw/CRQIcwAy8KCQw8gQ7DRwIdwgy8KCQQCfSlC3sngbo4zIZNWyvySV4/0BeISKAvmSguDiSwOPa+1YwE+paR4uLxVQJ79+6VveY+2CIyuNfcxbw2mAb1vejR3M1nQ+8NfT70aN5rGXoefTb8/shzaRn73nBZWnZ8n6gsfc8cEQ/VOVz38PvRe1F8Nrbh+qJ6h+pPxqxbxts0Kk4TW5G3VlN9T2e7bO/bXWQY1D1MIOq7UV8b3b8HTf8eem9fn4z6lek/lfrvmH466ntUoX/rd8fU0dHeKrt2697Nffs/hx/a3AAKbr1XEnj68hUyb+4suWHFhQVjSVc9EpiOVxm3RgLLmNXG2oQEjuXWOzgou81B8B5zdBXdzSZ6KKwH6bujg/+90es95nHo85Z975ttovejz/T50D3az9yT5Uaf6/vDB3e7R/YbOpDTcvQgTx+t2NgDv/h7Q5Jkth+WEysce+JiMnzAp3VaMbHSpLGKOdhra2uVPnOwt0+YrOwMCdPIAWdcqOIHoiPvD8cbSVNMxKw8DcfVpwVygwAEIBAAgb2f/XQAUZY3RK8ksNp1An3HjwT6nqHs40MCs2ccSg1FSKAe+PfvHZRd0aOYx0HR93YZaP1GwHYZaejX1/aefK3SZLbtM/sNRGXslQHzetfwvrqfvu43hqPC1Be9L9G2A6asXeb9aHt9be47TDnc/CbQZWSyVe+mv7SaUPW5jpXpiFWredZq8qjP9V37nv3cdIGhbey26rsjz81++jr2eZuWY/bZ996++qIy43cbR/SeqSMeW/y5ia9F44/HbOPYOza2NkYC/e6QOUen43qj++m+/jaqnw73Y+2nUb+P+nns+zL8PPoexfuu7ZtV+unsaV2ycVt/zq32r7rfmzfTv6CaKCIk0EGykUAHEAMvAgkMPIETCF+n1u0c3GPEZ0h+BvSI19zf27lLdprXO/eYu/ms1zzXbUcejUjp8z5Tt5UvHXVSmdJRJB3Biu7Do1p29Ct6barYafYP5abCoQfhbSbgNn1u7u3RQZh5jF6b982B+xC6oW2jz6PtYvvpc93GlKP76cGY/TzaVvc1fGwdttxoW7Ofbts6XK8e1OnnVjTsgZ+WPXIwZ/Iwsk0kREOftVU44NP94kKjMUwyo4BTuttl2/aBMSI1VO8QA72NHFzag04Vmbjg2PorHogOS4+K1vCBqDLn5g8Bzgn0JxdFR8I5gfsyMG9Wd9HpaOr6vZJAnQ667JQlcuayE4NKChIYVLoyCRYJzASrk0KjEbHoPjR6paNZu8yEuqH3hkavdhiJ224etxuxUpHT+zZ9NNtuN4/6vgrddvNat9Xn+plu48ttakurTDIm0WkkojN6NK9NcPp8knk+9CgySeXEbNsZPYp0mMdoW/M8+szsr593jHrdEm1nPx96PrSvbq+PHVq/Kb/D1NNh3kNCDE9z3s+0yR2yYSu/+PvyPSkqDiSwKPL+1YsEIoG+9EqvJFCvEfjDH98uq+6/yRc+dcWBBNaFqdQbIYHZpFcFbuuePfKRka2tRsS2GgGLnpvRNX2+zby3bfi9bZHMDcqW6H2znQpbTpI2o7VNJpt/2Se3tspU83yqOfjvMH7YY173GDnqGf6sy8hRt3neZd7Xx27zWfRoXrcbhNEIWDT6pSNdQ6NdOkKmkmVHz/Qxej1cVjbkKdUFASTQBcVylIEEliOPLlqBBCKBLvqRizK8kkA9J3C8G6uDukg5ZWRBAAmsTvWNgQH5YHC3mRa5R3rNZnY6pJ5Xps+3Dkudyps+32y222Lu68xrV7fJOkqlI17R6JaOWJkRrGi0TIaFrDWSuClG4KYMP1qpm2IETfefOix09n0VO5W/5K2IcwJdcaIctwSQQLc8Qy4NCQw5e25jRwKRQLc9qvHSvJLAxptR7J6MBBbL34fam1UC1+7eLR8awfvQSNuHe4Ye39k9IO/sGZA1Rv7eM+9N9DarrU2mGQGbYURsmpGu6Lm5T29tN2LWEr2ngqbTIfUzfT6tZeg9fZ33DQnMm7i/9SGB/uYm78iQwLyJ+1sfEogE+tI7vZNAnRK6/IrrR/FZee1lsuS4o3xhNiYOJNDb1OQWWFkkcKORtk3DUyo3msfN0es9sslMv1w3LHnrh6VPR+7quc1ra5eD2juiqZLR1MfhaYw6BVKnS+rImo6oRXezkMYMI3BDr81nBUhcPW0abxskcKIEy7M/ElieXE60JUjgRAmWZ38kEAn0pTd7JYH3PPCYXH3jXfLQ3dfIgvlzIkZr1q6TU876vlx5ydneLhiDBPrSnYuLw3cJfNOMyr1v5O0DI3Lvm5G696PH3bLBiNxGc99sJK9eqYtTPsAI3v5mpG6OPhpx08f55v7xjg45qK1DDjGPzXZDApst49XbiwTSFywBJJC+YAkggUigL98GryRw6WkXy/nnnDpG9lQOb7vzQW8XjEECfenOxcVRhAR+YKZe6sjdBiN366PnRubM6J2Knsqdjtyp6G2uc8RO6U03I2/7GZmbae77GbnT5/uZ0bnZ5nF/M/1yjnlv/2HxUwHkNpYAEkivsASQQPoCEkgfSBJAApFAX74VXklgtYvF2ymiLAzjS7chjiSBKxlwRwAAHv5JREFULCRQz61bY86te3dgtznHztz1tZ5vt3uXGcmrbyqmjfNAI2x6n9c+9DhXH43UzTXPVfZmGvmbjdQ56dhIoBOMpSgECSxFGp00gpFAJxhLUQgSiAT60pG9kkBGAn3pFsSRlkBaCVShe9vcdeEUHa3TUTsdwdMRvfeM+OkoX62bnjenUzEP0JE78zhbR+mM2M3S943k6WdW/mqVxefuCCCB7liGXhISGHoG3cWPBLpjGXpJSCAS6Esf9koCizonUOVzw6atIzlJLkRz6VW3yCOrno4+X7TwELlv5VWj8sd0UF+6c3FxVJLADTp6Z+5vDOyS143wvW4e3zKPL+yq78LROjK3IFpUxTy2T5KP6evhc+0+0YTn2hWX3XQ1I4HpeJV5aySwzNlN1zYkMB2vMm+NBCKBvvRvryRQoeS9OqguPHPDn/9UblhxYZQTK6J26mnyfMTTl6+QxUcvlMsvOGMkh0igL905/zjWmVG818z0zDf27pY3zePrfbuGLpFg7r3m4uXVbnqu3aEdk8x0zDaZZxZQOcA86tRMXVhFHw9G8vJPpqMakUBHIEtQDBJYgiQ6agIS6AhkCYpBApFAX7qxdxJYNBi7GqldoTQpfZUWqUECi85a9vVb2XvFjOa9akbyXjWS9yvz+JG5uHm1m17+QEfvDjaXR/ikEb5PDj8u7Ogs5Pp12VOiBiWABNIPLAEkkL5gCSCB9AVLAAlEAn35NiCBiUwkF6FJnqdYaZEaJNCX7jzxOHaY1TVfGuiXl43svWbuOnVzPNnT694daqZqHtndJYd1dcoBe4z4mVE8ncY5i4VWJp6QAEtAAgNMWkYhI4EZgQ2wWCQwwKRlFDISiARm1LVSF+uFBNopmJWuBTjeZ6lbW8cOKn1fP+n4kemeumJpPC4rgfFrGdZRLJt4RkDH717p65fVfX3yfG+vrO4denx714BUmsTZY/5qH97VFcnekcOPRxjpO3jSJHPpc24QgICPBPoH9khnR5uPoRETBCAAAQhAoFACXkigTrmcN3fWyHl5SSK6MMt7H2wcsyCLa3IqgMccdeioOBgJdE05//L03LyXzIjei2aET0f2XtrVJ78cZ3GWT5mRvcON3B1upm0eaaZxHjqpUz5upnKOd0u7Omj+FKgxLwKMBOZF2v96GAn0P0d5RchIYF6k/a+HkcB9OZo3q9v/hJU4Qi8ksNr1AS33PK4TWEkAtX7OCQyr928z5+g9298XyV4kfEb8XjXTOivddPXNI4zkLTTCd4QRPj1X7zNG+Bq5IYGNUCvnPkhgOfPaSKuQwEaolXMfJLCceW2kVUggEthIv8liHyTQUFUJPff0k0et+GlhszpoFt3OTZk7zfl7zxvRe86M7D1nZO958/jmwMCYwnW6pq62eWQ0stcpi4zofdrcXV4cHQl0k9MylIIEliGLbtqABLrhWIZSkMAyZNFNG5BAJNBNT5p4KV5IoI7C/egH58mS446q2CIdCfzhj2+XVfffNPEWJ0qw5xwmC/7K0sUj00K5TqBz7A0V+LwZ4XtahW9Y/HThlko3HdU70kieyt4ilb7OTpnS0tpQnfXuhATWS6r82yGB5c9xvS1EAuslVf7tkMDy57jeFiKBSGC9fSXr7byQwOtuvVeefvaVquf81TpnMGtItcpnddBahNJ/vtFcZP0ZI3vPGPF7ZldvNMUzed09HeE7xJyrd/SkLjOy1yVHG9lbZB51xc68b0hg3sT9rQ8J9Dc3eUeGBOZN3N/6kEB/c5N3ZEggEph3n6tWnxcSqMHpaKDekqN9+v6GTVvFXrzdF3DxOJDAiWflDTOq96/9vfIvfb1m0ZY+edtchy950+mbx5jRvWM7jfBF0teV+QhfvS1DAuslVf7tkMDy57jeFiKB9ZIq/3ZIYPlzXG8LkUAksN6+kvV23kigNlRHBO+47+FRbY5Py8waRqPlI4HpyOklGF42o3w/N9I3dO+TD83IX/zWbl4cbkTvWJW+rm7z2FVzhc50UbjdGgl0yzPk0pDAkLPnNnYk0C3PkEtDAkPOntvYkUAk0G2Parw0rySw8WYUuycSWJt/r1nE5fG+nfJo7w55xNw3D+4ZtZNO4Ty2s1u+YEb3vtDZE43yFTGts3ZLKm+BBDZKrnz7IYHly2mjLUICGyVXvv2QwPLltNEWIYFIYKN9x/V+SKADokhgZYjrzejew0b4VPz+xQhgn7len71NMwu1fM6M8H3BiN/nzV2lL+RLOiOBDr5IJSkCCSxJIh00Awl0ALEkRSCBJUmkg2YggUigg27kpAgk0AFGJHAIoires+Z8vn8w0vdY705znb6+6D29TTJ/9RabKZ1LunrkS91mpM88z3a9TgeJTVEEEpgCVsk3RQJLnuAUzUMCU8Aq+aZIYMkTnKJ5SCASmKK7ZLopEugAbzNLoI7uPWak73/t3C7/aEb74tM8DzeXZ/hSd7csNdM7P2/kL6TpnWm7BRKYllh5t0cCy5vbtC1DAtMSK+/2SGB5c5u2ZUggEpi2z2S1PRLogGyzSeAuI34qfA/u3BZN9dSLtuuty/xl+00je7/dPVlOMvd5ZjXPZrkhgc2S6drtRAJrM2qWLZDAZsl07XYigbUZNcsWSCAS6EtfRwIdZKJZJFBH/P7OjPj9z97tI+KnF2H/as8UOdncv2zO7etuLdMkz/o7BxJYP6uyb4kElj3D9bcPCayfVdm3RALLnuH624cEIoH195Zst0QCHfAtswS+Zq7fd/eOrfLAju2iC73oTcXvK2ak79TJU+UEM/LXUcDF2R2kzWkRSKBTnEEXhgQGnT6nwSOBTnEGXRgSGHT6nAaPBCKBTjvUBApDAicAz+5aNgnU8/p+uv0j+amZ7vmSuZ6fvf2uGe373Z6p8jXzyG00ASSQHmEJIIH0BUsACaQvWAJIIH3BEkACkUBfvg1IoINMlEECdRXPVWa65092fCSPmCmfA8NcPmMu1v7NydPkG2bUb0ZryBdxcJDocYpAArPlG1LpSGBI2co2ViQwW74hlY4EhpStbGNFApHAbHtY/aUjgfWzqrplyBK4wUzx/IkZ9bvHyN+a3UPqN7etTb7RM03+45Rp8omOSQ4Ilb8IJLD8Oa63hUhgvaTKvx0SWP4c19tCJLBeUuXfDglEAn3p5Uigg0yEJoE66ve4GfW7e/vWaHVPPdNPx/j+nTm/71tTZ8hvmfP9GPNL1zGQwHS8yrw1Eljm7KZrGxKYjleZt0YCy5zddG1DApHAdD0mu62RQAdsQ5HAreZcPx31u9PInx31O9BcxuGMKdPlLDPqp8+5NUYACWyMWxn3QgLLmNXG2oQENsatjHshgWXMamNtQgKRwMZ6jvu9kEAHTH2XwJcH+uX2bVvMCp/bRC/urhdx+HdmtO9bRv5OZNTPQQ8QQQKdYCxFIUhgKdLopBFIoBOMpSgECSxFGp00AglEAp10JAeFIIEOIPoogQNG9vR6fncY+Xumvy9qpS7scqZZ5OUPp81oqgu5O0hxzSKQwJqImmYDJLBpUl2zoUhgTURNswES2DSprtlQJBAJrNlJctoACXQA2icJ1Ms7/LURv7820z4/HL6u3xEdnXLu1OnyH4wAdnFNPwcZH1sEEpgJ1iALRQKDTFsmQSOBmWANslAkMMi0ZRI0EogEZtKxGigUCWwAWnIXHyTwVXNR9z838ve3ZpVPnfKpZ/f9jrme33+aMkOO7+p20EqKGI8AEkj/sASQQPqCJYAE0hcsASSQvmAJIIFIoC/fBiTQQSaKlMDHzOqeer7f/+7bGbVkemur/MfJ05ny6SCvaYpAAtPQKve2SGC585umdUhgGlrl3hYJLHd+07QOCUQC0/SXLLdFAh3QzVsCe81I30/NiN9ffrRFfr17V9SCQzo65Lwp+8np5qLu3UYEueVLAAnMl7fPtSGBPmcn39iQwHx5+1wbEuhzdvKNDQlEAvPtcdVrQwIdZCIvCdTz/XTU76+3bRV9rrcvdvbIH5mFXn7brPLZ4qAtFNEYASSwMW5l3AsJLGNWG2sTEtgYtzLuhQSWMauNtQkJRAIb6znu90ICHTDNWgLf3j0gf7Z1k/x/O4cu8dBp/oIs65kq50/bT36jY5KDFlDERAkggRMlWJ79kcDy5HKiLUECJ0qwPPsjgeXJ5URbggQigRPtQ672RwIdkMxKAp/b1S//79aN8og572/QxKmXeDjHXNvvPLPS5ywu7O4gc+6KQALdsQy9JCQw9Ay6ix8JdMcy9JKQwNAz6C5+JBAJdNebJlYSEjgxftHeLiVwrynv4Z3b5bZtm+Xp4ev7zTfCd97U/eTsKdM4389BvrIoAgnMgmqYZSKBYeYti6iRwCyohlkmEhhm3rKIGglEArPoV42UiQQ2Qi2xjwsJ1Iu7/3ez2MufmcVeXh9e7OWISZ1y/tQZssxc36/NQZwUkR0BJDA7tqGVjASGlrHs4kUCs2MbWslIYGgZyy5eJBAJzK53pSsZCUzHq+LWE5HA3sFBuWv71ugaf+8PX9x9SVePXGDk7wSz2Au3MAgggWHkKY8okcA8KIdRBxIYRp7yiBIJzINyGHUggUigLz0VCXSQiUYkcEu00udW+Ssjf/pcL+pwirm4+/fMYi+LJnU5iIoi8iSABOZJ2++6kEC/85NndEhgnrT9rgsJ9Ds/eUaHBCKBefa38epCAh1kIo0EfrBnj9z60Sb5iZn6udOMAnaY+k8z0z0vmjYzutYftzAJIIFh5i2LqJHALKiGWSYSGGbesogaCcyCaphlIoFIoC89Fwl0kIl6JPDNgQG5ycjfA0b+BkydPeaC7mcZ+fuOGfmbx0qfDrJQbBFIYLH8faodCfQpG8XGggQWy9+n2pFAn7JRbCxIIBJYbA/cVzsS6CAT40ngqwO75HpzmYe/Myt+6sqfepmH/2TO9/u2uczDfuY5t3IQQALLkUcXrUACXVAsRxlIYDny6KIVSKALiuUoAwlEAn3pyUhgLBOnL18hi49eKJdfcMao/Fx61S3yyKqno/cWLTxE7lt51ajPK0ngywP9ct2WoWv8qfzNbWuT70SXeZgejQJyKxcBJLBc+ZxIa5DAidAr175IYLnyOZHWIIEToVeufZFAJNCXHo0EmkzEJe/c008eJYH3PPCY3Hbng7Lq/puinFUSxbgE6gXedeTvMSN/elP5u9Cc73eWkb9J+s3nVkoCSGAp09pQo5DAhrCVcicksJRpbahRSGBD2Eq5ExKIBPrSsZHAWCaWnnaxfP2k40dJYFL6klKou6sEPmMu7K7yt6pvZ1TiHDPV87vTZ8q3jPx1In++9PfM4kACM0MbXMFIYHApyyxgJDAztMEVjAQGl7LMAkYCkcDMOlfKgpHAGhKoYnj+OafKmctOjLZ84qnVsvyK6+XFx+8Yer1jh/zg7ffkX/p7o9f7m0Ve/rNO+zTn/HUjfym7Y7ibI4Hh5s515Eiga6LhlocEhps715Ejga6JhlseEogE+tJ7kcAaEnjkCefKlZecPUYCH7r7Glkwf460/PL5qITZ7W1yyX6z5I/2m4n8+dK7c4yjo71VWltbpH/XnhxrpSofCbS3tUpHe4v09tMXis7P4ODe6HtZ1E0P/DsntcnOvt1FhUC9nhBoMUf+PV1tsqOXvuBJSooLw/xJmtLVIdt7da345r5N7eHSaEX2ACSwhgTWGgn8+utvynGd3bLcyB8LvhTZlYutW3/x12PNvoHBYgOh9sIJtLe1SIcRwV5+ECg8F15IYEer7OQHgcL7QtEB6L8PPZ3tsp0fBIpOReH1689SU7rbZRs/CMhUw4FbcQSQwBoSWO85gcWlkJp9IMB0UB+y4EcMTAf1Iw8+RMF0UB+y4EcMTAf1Iw8+RMF00H1ZmDer24eUNG0MSGANCUy7OmjT9qQmbzgS2OQdINZ8JJC+YAkggfQFSwAJpC9YAkggEujLtwEJNJmIXyLCJsYu/KKvG7lOoC8JJo58CCCB+XAOoRYkMIQs5RMjEpgP5xBqQQJDyFI+MSKBSGA+Pa12LUhgbUY1t6h0sfiaO7FBqQgggaVK54QagwROCF+pdkYCS5XOCTUGCZwQvlLtjAQigb50aCTQQSaQQAcQAy8CCQw8gQ7DRwIdwgy8KCQw8AQ6DB8JdAgz8KKQQCTQly6MBDrIBBLoAGLgRSCBgSfQYfhIoEOYgReFBAaeQIfhI4EOYQZeFBKIBPrShZFAB5lAAh1ADLwIJDDwBDoMHwl0CDPwopDAwBPoMHwk0CHMwItCApFAX7owEuggE0igA4iBF4EEBp5Ah+EjgQ5hBl4UEhh4Ah2GjwQ6hBl4UUggEuhLF0YCfckEcUAAAhCAAAQgAAEIQAACEMiBABKYA2SqgAAEIAABCEAAAhCAAAQg4AsBJNCXTBAHBCAAAQhAAAIQgAAEIACBHAgggTlApgoIQAACEIAABCAAAQhAAAK+EEACG8zEpVfdIo+sejrae9HCQ+S+lVc1WBK7hUIgTc7veeAxufrGu0aaRh8JJcv1xZmmL8RLtP1i5bWXyZLjjqqvMrbymkDyu/7i43fUjDfef76ydLHcsOLCmvuwgf8E0vaF6269V+647+GRhvF3wf8cu4rw9OUrZPHRC+XyC85wVSTlQCA1ASQwNTIR/UN/250Pyqr7b4r25svcAMTAdkmbcz3Iu/Q735QF8+dELV162sXy9ZOO5w9+YHmvFG7avmDLsPtt2LRVONgrQUcwTVizdp2cctb35aG7r4m+63pQ//Szr4z7o6D+e6E3fjgsRx+wrUjbF6ww2h8Nkq/LRYfWWALxH4DOPf1kjgnoGoUSQAIbwJ+UvuRBYQNFsovnBCaac/3Drzd+8fc80XWE10hfeOKp1fLDH98e/XB05AnnIoF1cA5hk6T0JUUg2YZ4PwihfcRYP4G0fSHt9vVHwpYhEOCH4RCyVP4YkcAGcqxf3vPPOVXOXHZitLf+w778iuulnmlADVTHLh4QmGjOGS32IImOQkjbF5J/H5BAR4nwoJhKP+6Ml1898P/Zo0+Kjgbb25WXnD3yb4kHTSKEBgmk7QtajfaV2TOnRz8Oad9474MN/FDYIP/QdkMCQ8tYOeNFAhvIq/7hjv/DbQ/y7JSgBopkF88JTCTnTPPxPLkpw0vTFyqNDCGBKYF7vLn+uDNv7qxRB+7J/hEPX0XhvQ82jkwFtf92MD3Y4yTXGVravqDF6j4frN808qMA/aBO2CXYDAksQRJL0AQksIEkph0JaKAKdvGMQKM5twLIDwSeJXQC4aTpC8mFIuLVMgI0gSR4smva0Z9K2yf7kydNI4yUBNL2heTIHz8IpAQe+OZIYOAJLEn4SGADiWzknKAGqmEXjwg0knME0KMEOgylkb4Qr56RQIfJKLiotOd1VVo4BgksOImOqk/bFyqdIjDeKLKjMCnGEwJIoCeJaPIwkMAGOkCjqwM2UBW7eEKgVs7tUt/2vNDka0+aQRgOCKTtC8kqkUAHSfCkiForQib/Dtjt7bQ/pop7kkgHYdTqC8mRvmTf4EdDB0kIqAgkMKBklThUJLDB5DZ6nbAGq2M3DwiMl/PkP+j6Bz6++IMNn2mhHiTSQQhp+gIS6AC4x0WMd224Sj8GWRngb4LHSW0wtPH6QqXpnvG/I1ol/z40CD6g3ZI519BZVDCgBJYsVCSwZAmlORCAAAQgAAEIQAACEIAABMYjgATSPyAAAQhAAAIQgAAEIAABCDQRASSwiZJNUyEAAQhAAAIQgAAEIAABCCCB9AEIQAACEIAABCAAAQhAAAJNRAAJbKJk01QIQAACEIAABCAAAQhAAAJIIH0AAhCAAAQgAAEIQAACEIBAExFAApso2TQVAhCAAAQgAAEIQAACEIAAEkgfgAAEIAABCEAAAhCAAAQg0EQEkMAmSjZNhQAEIAABCEAAAhCAAAQggATSByAAAQhAAAIQgAAEIAABCDQRASSwiZJNUyEAAQhAAAIQgAAEIAABCCCB9AEIQAACEIAABCAAAQhAAAJNRAAJbKJk01QIQAACEIAABCAAAQhAAAJIIH0AAhCAAAQgAAEIQAACEIBAExFAApso2TQVAhCAAAQgAAEIQAACEIAAEkgfgAAEIAABCEAAAhCAAAQg0EQEkMAmSjZNhQAEIAABCEAAAhCAAAQggATSByAAAQhAAAIQgAAEIAABCDQRASSwiZJNUyEAAQg0QuC6W++VO+57eMyuLz5+h9zzwGNy9Y13yZWXnC1nLjtx1DaXXnWLPLLqaXno7mtkwfw5Ml45Tzy1WpZfcf244WkdetP6Kt1sDDamRQsPkftWXjVqUxtDpc/shkeecO64cXxl6eLoc22bvZ17+sly+QVnyJq16+SUs74fva184rd4G+1nNtbx2tNIztgHAhCAAAQgMB4BJJD+AQEIQAACVQlYkUsKjb5/3GcPj8TPilV8Gys8K6+9TJYcd5TUU048CN3+vQ82jpE4K03JeOL7xsXKCqj9fOlpF8uGTVtlPAmsVFal+mwb43XEJTApxpZBXBDraQ/dEwIQgAAEIOCaABLomijlQQACECgRAR0Vs6Nc4zVL5eqYow6VG1ZcGG2WfF1vObaOiUrgbXc+GMWjNxuTCpe+P3f/mdH7yVHCSu0bT9LGk0Bl9rNHn5RV998UFWvlUEcRdQQxORI4ntSWqDvRFAhAAAIQ8IQAEuhJIggDAhCAgI8EkjJXLcb4yN+Tz7w4SoAqSWGttrqQwLtu/uNoaqYdqTt9+QpZdsoSeeChJzKXQK3z7Iv+RM4/59SR0dL3Ptgg8+bOjqbWIoG1egCfQwACEIBAlgSQwCzpUjYEIACBwAlUOmet2qhVfLqjnQZqm5+mHN2nlgRWwho/J1BH/HQUTsVv3txZ8h++tlR++OPbR97T/bMcCVQJvO/BfxqRYR0Jte9VksDx2hN4FyJ8CEAAAhDwkAAS6GFSCAkCEICAjwTs+XQ2tqTo6fsqOzrl0U7BrNSOesqpJYG1zgm0EmhHKPUcQB0F1HMYVQzzkEBdDEd5aN0qosokef4k5wT62NOJCQIQgED5CSCB5c8xLYQABCDgnICKnJ5blxxNq3f6qA2oWjmuJFDrUel74ZU3R6Zg5imBVvqsMCOBzrsiBUIAAhCAQAMEkMAGoLELBCAAgWYgoIuZ6JRGvfRB8lZNpCpJYCPluJRAHW179/31I+3IUwK17Zf/37eNyDIS2AzfHNoIAQhAwH8CSKD/OSJCCEAAAoUQsCtazp45fWSVSw3ETmGsNB20mgTqAi1pynEpgfUKbCXIja4Omrw0hS0bCSykK1MpBCAAAQgkCCCBdAkIQAACEBiXQKWLp1eTnPGmg6Ypp5YEVgq40sIwlbZzMRIYXwRH60heLD6tBI7XHronBCAAAQhAwDUBJNA1UcqDAAQgAAEIQAACEIAABCDgMQEk0OPkEBoEIAABCEAAAhCAAAQgAAHXBJBA10QpDwIQgAAEIAABCEAAAhCAgMcEkECPk0NoEIAABCAAAQhAAAIQgAAEXBNAAl0TpTwIQAACEIAABCAAAQhAAAIeE0ACPU4OoUEAAhCAAAQgAAEIQAACEHBNAAl0TZTyIAABCEAAAhCAAAQgAAEIeEwACfQ4OYQGAQhAAAIQgAAEIAABCEDANQEk0DVRyoMABCAAAQhAAAIQgAAEIOAxASTQ4+QQGgQgAAEIQAACEIAABCAAAdcEkEDXRCkPAhCAAAQgAAEIQAACEICAxwSQQI+TQ2gQgAAEIAABCEAAAhCAAARcE0ACXROlPAhAAAIQgAAEIAABCEAAAh4TQAI9Tg6hQQACEIAABCAAAQhAAAIQcE0ACXRNlPIgAAEIQAACEIAABCAAAQh4TAAJ9Dg5hAYBCEAAAhCAAAQgAAEIQMA1ASTQNVHKgwAEIAABCEAAAhCAAAQg4DEBJNDj5BAaBCAAAQhAAAIQgAAEIAAB1wSQQNdEKQ8CEIAABCAAAQhAAAIQgIDHBJBAj5NDaBCAAAQgAAEIQAACEIAABFwTQAJdE6U8CEAAAhCAAAQgAAEIQAACHhNAAj1ODqFBAAIQgAAEIAABCEAAAhBwTQAJdE2U8iAAAQhAAAIQgAAEIAABCHhMAAn0ODmEBgEIQAACEIAABCAAAQhAwDUBJNA1UcqDAAQgAAEIQAACEIAABCDgMQEk0OPkEBoEIAABCEAAAhCAAAQgAAHXBJBA10QpDwIQgAAEIAABCEAAAhCAgMcEkECPk0NoEIAABCAAAQhAAAIQgAAEXBNAAl0TpTwIQAACEIAABCAAAQhAAAIeE0ACPU4OoUEAAhCAAAQgAAEIQAACEHBNAAl0TZTyIAABCEAAAhCAAAQgAAEIeEwACfQ4OYQGAQhAAAIQgAAEIAABCEDANQEk0DVRyoMABCAAAQhAAAIQgAAEIOAxASTQ4+QQGgQgAAEIQAACEIAABCAAAdcEkEDXRCkPAhCAAAQgAAEIQAACEICAxwSQQI+TQ2gQgAAEIAABCEAAAhCAAARcE0ACXROlPAhAAAIQgAAEIAABCEAAAh4TQAI9Tg6hQQACEIAABCAAAQhAAAIQcE0ACXRNlPIgAAEIQAACEIAABCAAAQh4TOD/B/Cc+OTNBhqYAAAAAElFTkSuQmCC", "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 }