examples/simulation.ipynb

{
 "cells": [
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Imports"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from paraflow import SymmetricPassage, get_flasher"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Passage"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.plotly.v1+json": {
       "config": {
        "plotlyServerURL": "https://plot.ly"
       },
       "data": [
        {
         "name": "Contour Top",
         "type": "scatter",
         "x": [
          0,
          0.024242449999575005,
          0.04800040799326811,
          0.07128917372863348,
          0.09412404695322528,
          0.11652032741459765,
          0.13849331486030483,
          0.1600583090379009,
          0.18123060969494,
          0.20202551657897644,
          0.22245832943756422,
          0.24254434801825767,
          0.26229887206861086,
          0.2817372013361779,
          0.3008746355685131,
          0.3197264745131705,
          0.3383080179177044,
          0.35663456552966877,
          0.3747214170966179,
          0.3925838723661059,
          0.41023723108568705,
          0.4276967930029154,
          0.44497785786534516,
          0.46209572542053057,
          0.4790656954160256,
          0.49590391758535984,
          0.5126460913394928,
          0.5293474657668148,
          0.5460641399416909,
          0.5628522129384864,
          0.5797677838315667,
          0.5968669516952969,
          0.6142058156040425,
          0.6318404746321685,
          0.6498270278540405,
          0.6682215743440233,
          0.6870802131764825,
          0.7064590434257835,
          0.7264141641662912,
          0.747001674472371,
          0.7682776734183885,
          0.7902982600787087,
          0.8131195335276967,
          0.8367975928397181,
          0.861388537089138,
          0.8869484653503216,
          0.9135334766976344,
          0.9411996702054416,
          0.9700031449481084,
          1
         ],
         "y": [
          0.2,
          0.1919526120239622,
          0.18558013691077682,
          0.18079643761021705,
          0.17751537707205603,
          0.17565081824606682,
          0.17511662408202275,
          0.1758266575296968,
          0.1776947815388622,
          0.18063485905929205,
          0.1845607530407595,
          0.18938632643303774,
          0.19502544218589993,
          0.20139196324911912,
          0.20839975257246857,
          0.21596267310572137,
          0.22399458779865064,
          0.2324093596010296,
          0.24112085146263126,
          0.25004292633322894,
          0.2590894471625957,
          0.26817427690050466,
          0.27721127849672905,
          0.286114314901042,
          0.2947972490632165,
          0.30317574261807634,
          0.3112068269566033,
          0.318888903225938,
          0.3262221712582713,
          0.3332068308857941,
          0.3398430819406972,
          0.3461311242551716,
          0.3520711576614082,
          0.3576633819915977,
          0.3629079970779312,
          0.3678052027525993,
          0.37235519884779317,
          0.3765581851957036,
          0.38041436162852144,
          0.38392392797843755,
          0.387087084077643,
          0.3899040297583284,
          0.39237496485268475,
          0.3945000891929029,
          0.39627960261117384,
          0.39771370493968833,
          0.39880259601063733,
          0.3995464756562117,
          0.3999455437086023,
          0.4
         ]
        },
        {
         "name": "Contour Bottom",
         "type": "scatter",
         "x": [
          0,
          0.024242449999575005,
          0.04800040799326811,
          0.07128917372863348,
          0.09412404695322528,
          0.11652032741459765,
          0.13849331486030483,
          0.1600583090379009,
          0.18123060969494,
          0.20202551657897644,
          0.22245832943756422,
          0.24254434801825767,
          0.26229887206861086,
          0.2817372013361779,
          0.3008746355685131,
          0.3197264745131705,
          0.3383080179177044,
          0.35663456552966877,
          0.3747214170966179,
          0.3925838723661059,
          0.41023723108568705,
          0.4276967930029154,
          0.44497785786534516,
          0.46209572542053057,
          0.4790656954160256,
          0.49590391758535984,
          0.5126460913394928,
          0.5293474657668148,
          0.5460641399416909,
          0.5628522129384864,
          0.5797677838315667,
          0.5968669516952969,
          0.6142058156040425,
          0.6318404746321685,
          0.6498270278540405,
          0.6682215743440233,
          0.6870802131764825,
          0.7064590434257835,
          0.7264141641662912,
          0.747001674472371,
          0.7682776734183885,
          0.7902982600787087,
          0.8131195335276967,
          0.8367975928397181,
          0.861388537089138,
          0.8869484653503216,
          0.9135334766976344,
          0.9411996702054416,
          0.9700031449481084,
          1
         ],
         "y": [
          -0.2,
          -0.1919526120239622,
          -0.18558013691077682,
          -0.18079643761021705,
          -0.17751537707205603,
          -0.17565081824606682,
          -0.17511662408202275,
          -0.1758266575296968,
          -0.1776947815388622,
          -0.18063485905929205,
          -0.1845607530407595,
          -0.18938632643303774,
          -0.19502544218589993,
          -0.20139196324911912,
          -0.20839975257246857,
          -0.21596267310572137,
          -0.22399458779865064,
          -0.2324093596010296,
          -0.24112085146263126,
          -0.25004292633322894,
          -0.2590894471625957,
          -0.26817427690050466,
          -0.27721127849672905,
          -0.286114314901042,
          -0.2947972490632165,
          -0.30317574261807634,
          -0.3112068269566033,
          -0.318888903225938,
          -0.3262221712582713,
          -0.3332068308857941,
          -0.3398430819406972,
          -0.3461311242551716,
          -0.3520711576614082,
          -0.3576633819915977,
          -0.3629079970779312,
          -0.3678052027525993,
          -0.37235519884779317,
          -0.3765581851957036,
          -0.38041436162852144,
          -0.38392392797843755,
          -0.387087084077643,
          -0.3899040297583284,
          -0.39237496485268475,
          -0.3945000891929029,
          -0.39627960261117384,
          -0.39771370493968833,
          -0.39880259601063733,
          -0.3995464756562117,
          -0.3999455437086023,
          -0.4
         ]
        }
       ],
       "layout": {
        "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": "Flow Passage"
        },
        "yaxis": {
         "scaleanchor": "x"
        }
       }
      }
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "passage = SymmetricPassage(\n",
    "    inlet_radius=0.2,\n",
    "    area_ratio=2,\n",
    "    axial_length=1,\n",
    "    contour_props=[0.2, 0.5, 0.75],\n",
    "    contour_angles=np.radians([-20, 15.0, 15.0]).tolist()\n",
    ")\n",
    "passage.visualize()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Inlet and Outlet Condtions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "flasher = get_flasher(\"Octamethyltrisiloxane\", \"gas\")\n",
    "inlet_total_state = flasher.flash(P=904388, T=542.13, mach_number=1E-9)\n",
    "target_outlet_static_state = flasher.flash(P=200000, T=293.15)"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Run Simulation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Info    : Meshing 1D...\n",
      "Info    : [  0%] Meshing curve 1 (Nurb)\n",
      "Info    : [ 30%] Meshing curve 2 (Line)\n",
      "Info    : [ 50%] Meshing curve 3 (Line)\n",
      "Info    : [ 80%] Meshing curve 4 (Line)\n",
      "Info    : Done meshing 1D (Wall 0.00183667s, CPU 0.002025s)\n",
      "Info    : Meshing 2D...\n",
      "Info    : Meshing surface 1 (Transfinite)\n",
      "Info    : Done meshing 2D (Wall 0.00102575s, CPU 0.000764s)\n",
      "Info    : Meshing 3D...\n",
      "Info    : Done meshing 3D (Wall 3.2042e-05s, CPU 1.1e-05s)\n",
      "Info    : 10204 nodes 10407 elements\n",
      "Setting up SU2 Simulation for /Users/afshawnlotfi/Documents/paraflow/simulation_out/config1.cfg\n",
      "Running SU2 Simulation for /Users/afshawnlotfi/Documents/paraflow/simulation_out/config1.cfg\n",
      "Cleaning up ..\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuMAAAGGCAYAAADRrI/QAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAACXwElEQVR4nO2de5wU1Zm/v1XVt7kDAjNAUFDjXSGBMMHoqptJRuMa2SQGTVYQEXdV8tOMJitJBDSueItBIwnrhWASEeOuMdloiDoJuq5EIkiiRkkwKKgMghcGBpmZ7jq/P7qruqq6qrqquvo634dPfaa76pyq090z3U+/vOc9khBCgBBCCCGEEFJy5HIPgBBCCCGEkKEKZZwQQgghhJAyQRknhBBCCCGkTFDGCSGEEEIIKROUcUIIIYQQQsoEZZwQQgghhJAyQRknhBBCCCGkTFDGCSGEEEIIKROUcUIIIYQQQsoEZZwQQkhorFy5EpIk4fnnny/3UAghpCqgjBNCKhpJkvJuixcvLvcwHXn88ccxd+5cHHfccVAUBRMmTCjatSZMmABJktDR0WF7/O6779afM8oyIYRUBpFyD4AQQtz46U9/6nhs8eLFeO2119De3l7CEflj1apVePDBB/Hxj38cY8eOLfr1EokEfv/736OnpwdtbW2mY/fffz8SiQQOHDhQ9HEQQgjxBmWcEFLR/Mu//Ivt/nvuuQevvfYavva1r+GMM84o8ai8c8MNN+Duu+9GNBrFP/3TP+Gll14q6vU+9alP4Y9//CMefPBBXH755fr+N998E//7v/+Lf/7nf8Z///d/F3UMhBBCvMM0FUJI1fHyyy/j//2//4ePfexjuOWWW0zH+vr6cOWVV2L8+PGIx+M48sgjceutt0IIYWonSRLmz5+PRx55BMcddxzi8TiOPfZYrFmzJud6b731Fi688EK0trbq7VasWOFprGPHjkU0Gg3+YH2SSCTwhS98AatWrTLtf+CBBzB8+HB0dnbm9Pnzn/+MCy64AIceeigSiQTa2tpw4YUX4t13381p+9Zbb2Hu3LkYO3Ys4vE4Jk6ciEsuuQQDAwOmdv39/ejq6sKoUaPQ0NCAf/7nf8auXbvCfbCEEFIDMDJOCKkq9u/fjy9/+ctQFAWrV69GPB7Xjwkh8PnPfx6///3vMXfuXEyePBm//e1v8Y1vfANvvfUWvv/975vO9cwzz+Dhhx/GpZdeiqamJtxxxx344he/iG3btuGggw4CAOzcuROf/OQndXkfNWoUfvOb32Du3Lno7e3FFVdcUcqH74mvfOUr+OxnP4vXXnsNhx12GIB0usyXvvQl2y8GTzzxBP7+979jzpw5aGtrw8svv4y77roLL7/8Mv7whz9AkiQAwNtvv41p06bhgw8+wMUXX4yjjjoKb731Fv7rv/4L+/fvRywW08/5ta99DcOHD8eiRYvw+uuvY+nSpZg/fz4efPDB0jwJhBBSLQhCCKkiLrzwQgFA3HfffTnHHnnkEQFAXH/99ab9X/rSl4QkSWLLli36PgAiFouZ9v3pT38SAMQPfvADfd/cuXPFmDFjxO7du03nPPfcc0VLS4vYv3+/57GfeeaZ4pBDDvHc3i+HHHKIOPPMM0UymRRtbW3iu9/9rhBCiL/85S8CgHjqqafEj3/8YwFA/PGPf9T72T2GBx54QAAQTz/9tL5v1qxZQpZlU18NVVWFEEI/f0dHh75PCCG+/vWvC0VRxAcffBDa4yWEkFqAaSqEkKph1apVWLFiBc4//3zMmjUr5/hjjz0GRVHw//7f/zPtv/LKKyGEwG9+8xvT/o6ODj1yDAAnnHACmpub8fe//x1AOtL+3//93zjrrLMghMDu3bv1rbOzE3v27MHGjRuL8EgLQ1EUfPnLX8YDDzwAID1xc/z48Tj55JNt29fV1em3Dxw4gN27d+OTn/wkAOiPT1VVPPLIIzjrrLMwderUnHNo0XONiy++2LTv5JNPRiqVwhtvvFHYgyOEkBqDaSqEkKrgb3/7G/7t3/4NRxxxBH74wx/atnnjjTcwduxYNDU1mfYfffTR+nEjBx98cM45hg8fjvfffx8AsGvXLnzwwQe46667cNddd9le85133vH9WLyya9cupFIp/X5jYyMaGxs99f3KV76CO+64A3/605+watUqnHvuuTnCrPHee+/h2muvxerVq3Mez549e/Sx9Pb24rjjjvN0fetzO3z4cADQn1tCCCFpKOOEkIqnv78fM2fOxMDAAFavXu1ZSPOhKIrtfpGZ7KmqKoB0RZfZs2fbtj3hhBNCGYsdn/jEJ0xfIBYtWuS5pnp7ezsOO+wwXHHFFdi6dSu+8pWvOLb98pe/jGeffRbf+MY3MHnyZDQ2NkJVVZx++un6c+CXfM8tIYSQNJRxQkjFc9VVV+GFF17A7bffjo997GOO7Q455BA8+eST2Lt3ryk6/uqrr+rH/TBq1Cg0NTUhlUo5LqRTTO6//358+OGH+v1DDz3UV//zzjsP119/PY4++mhMnjzZts3777+P7u5uXHvttVi4cKG+/29/+5up3ahRo9Dc3Fz00oyEEDLUoIwTQiqaX/ziF7jzzjvx+c9/PicX3MrnPvc53HXXXbjzzjuxYMECff/3v/99SJLkux65oij44he/iFWrVuGll17KSdHYtWsXRo0a5eucfvjUpz5VUP+LLroIiqK4LoqkRbCtEeulS5ea7suyjBkzZuBnP/sZnn/++Zy8cSGEYxoMIYQQZyjjhJCKZceOHZg7dy4URcGnP/1p/OxnP7Ntd9hhh2H69Ok466yzcNppp+Hb3/42Xn/9dUyaNAmPP/44fvnLX+KKK64wTdb0yo033ojf//73aG9vx7x583DMMcfgvffew8aNG/Hkk0/ivffec+3/5z//Gb/61a8AAFu2bMGePXtw/fXXAwAmTZqEs846y/eYvHLIIYfkTWtpbm7GP/zDP+Dmm2/G4OAgxo0bh8cffxxbt27NaXvDDTfg8ccfxymnnIKLL74YRx99NHbs2IGHHnoIzzzzDIYNG1acB0IIITUMZZwQUrFs3rxZn/BnXE3SyuzZszF9+nTIsoxf/epXWLhwIR588EH8+Mc/xoQJE3DLLbfgyiuvDDSG1tZWrF+/Htdddx0efvhh/PCHP8RBBx2EY489FjfddFPe/hs3bsQ111xj2qfdnz17dlFl3CurVq3C1772NSxbtgxCCHz2s5/Fb37zG4wdO9bUbty4cXjuuedwzTXX4P7770dvby/GjRuHM844A/X19WUaPSGEVDeS4GwaQgghhBBCygLrjBNCCCGEEFImKOOEEEIIIYSUCco4IYQQQgghZYIyTgghhBBCSJmgjBNCCCGEEFImKOOEEEIIIYSUCdYZz4Oqqnj77bfR1NTE1eUIIYQQUtUIIbB3716MHTsWslz+mOyBAwcwMDAQqG8sFkMikQh5RGVAEFe2b98uAHDjxo0bN27cuNXMtn379nIrlvjwww9F22gl8GNoa2sTH374oadr3XDDDWLq1KmisbFRjBo1Spx99tni1Vdfzdvv5z//uTjyyCNFPB4Xxx13nHj00UdNx1VVFddcc41oa2sTiURCfPrTnxZ//etffT0PVRcZX7ZsGW655Rb09PRg0qRJ+MEPfoBp06bZtl25ciXmzJlj2hePx3HgwAHP12tqagIAbN++Hc3NzcEHTgghhBBSZnp7ezF+/Hjdb8rJwMAAet5J4Y0NE9Dc5C9K37tXxSFTXsfAwICn6PhTTz2Fyy67DJ/4xCeQTCbxrW99C5/97Gfxl7/8BQ0NDbZ9nn32WZx33nlYsmQJ/umf/gmrVq3CjBkzsHHjRhx33HEAgJtvvhl33HEH7rvvPkycOBHXXHMNOjs78Ze//MVz1L6qVuB88MEHMWvWLCxfvhzt7e1YunQpHnroIWzevBmjR4/Oab9y5Upcfvnl2Lx5s75PkiS0trZ6vmZvby9aWlqwZ88eyjghhBBCqppK8hptLO/+dWIgGT/oiK2BH8euXbswevRoPPXUU/iHf/gH2zYzZ85EX18ffv3rX+v7PvnJT2Ly5MlYvnw5hBAYO3YsrrzySlx11VUAgD179qC1tRUrV67Eueee62ks5U8W8sFtt92GefPmYc6cOTjmmGOwfPly1NfXY8WKFY59JElCW1ubvvkRcUIIIYQQUlxSQg20AWmhN279/f2errlnzx4AwIgRIxzbrFu3Dh0dHaZ9nZ2dWLduHQBg69at6OnpMbVpaWlBe3u73sYLVSPjAwMD2LBhg+kBy7KMjo4O1we8b98+HHLIIRg/fjzOPvtsvPzyy6UYLiGEEEII8YAKEWgDgPHjx6OlpUXflixZkv96qoorrrgCn/rUp/R0Ezt6enpygritra3o6enRj2v7nNp4oWpyxnfv3o1UKmX7gF999VXbPkceeSRWrFiBE044AXv27MGtt96KE088ES+//DI+8pGP2Pbp7+83favq7e0N70EQQgghhBATKlSoAfoAuXP64vF43r6XXXYZXnrpJTzzzDM+r1ocqiYyHoTp06dj1qxZmDx5Mk455RQ8/PDDGDVqFP7zP//Tsc+SJUtM37DGjx9fwhETQgghhAwtUkIE2gCgubnZtOWT8fnz5+PXv/41fv/73zsGZjXa2tqwc+dO076dO3eira1NP67tc2rjhaqR8ZEjR0JRlIIecDQaxcc+9jFs2bLFsc2CBQuwZ88efdu+fXtB4yaEEEIIIc4UkqbiFSEE5s+fj1/84hf43e9+h4kTJ+btM336dHR3d5v2PfHEE5g+fToAYOLEiWhrazO16e3txXPPPae38ULVyHgsFsOUKVNMD1hVVXR3d3t+wKlUCi+++CLGjBnj2CYej+d8yyKEEEIIIdXLZZddhp/97GdYtWoVmpqa0NPTg56eHnz44Yd6m1mzZmHBggX6/csvvxxr1qzB9773Pbz66qtYvHgxnn/+ecyfPx9AukjIFVdcgeuvvx6/+tWv8OKLL2LWrFkYO3YsZsyY4XlsVZMzDgBdXV2YPXs2pk6dimnTpmHp0qXo6+vTa4nPmjUL48aN05P3r7vuOnzyk5/E4Ycfjg8++AC33HIL3njjDVx00UXlfBiEEEIIISSDCoGUz0i338j4j370IwDAqaeeatr/4x//GBdccAEAYNu2baZVSU888USsWrUK3/nOd/Ctb30LH/3oR/HII4+YJn1+85vfRF9fHy6++GJ88MEHOOmkk7BmzRpfK4NWlYzPnDkTu3btwsKFC9HT04PJkydjzZo1+qRO65P4/vvvY968eejp6cHw4cMxZcoUPPvsszjmmGPK9RAIIYQQQoiBIGknQdJU8rF27dqcfeeccw7OOeccxz6SJOG6667Ddddd52s8pnNU06I/5aCSiuMTQgghhBRCJXmNNpa/vtKKJp+L/uzdq+KIo3dWxOMolKqKjBNCCCGEkNpCzWx++9QKlHFCCCGEEFI2UgFyxv22r2SqppoKIYQQQgghtQYj44QQQgghpGykRHrz26dWoIwTQgghhJCywZxxQgghhBBCyoQKCSlIvvvUCpRxQgghhBBSNlSR3vz2qRUo44QQQgghpGykAkTG/bavZFhNhRBCCCGEkDLByDghhBBCCCkbQz0yThknhBBCCCFlQxUSVOFzAqfP9pUMZZwQQgghhJQNRsYJIYQQQggpEynISPmcxpgq0ljKAWWcEEIIIYSUDREgTUUwTYUQQgghhJDCGeppKixtSAghhBBCSJlgZJwQQgghhJSNlJCREj5zxrkCJyGEEEIIIYWjQoLqM1lDRe3YOGWcEEIIIYSUjaGeM04ZJ4QQQgghZSNYmgoj44QQQgghhBRMOk3F5wqcNRQZZzUVQgghhBBCygQj44QQQgghpGyoAVbg5AROQgghhBBCQoA544QQQgghhJQJFfKQLm3InHFCCCGEEFI2UkIKtPnh6aefxllnnYWxY8dCkiQ88sgjru0vuOACSJKUsx177LF6m8WLF+ccP+qoo3w/fso4IYQQQggpG6lMzrjfzQ99fX2YNGkSli1b5qn97bffjh07dujb9u3bMWLECJxzzjmmdscee6yp3TPPPONrXADTVAghhBBCSI1zxhln4IwzzvDcvqWlBS0tLfr9Rx55BO+//z7mzJljaheJRNDW1lbQ2BgZJ4QQQgghZUMVcqCtlNx7773o6OjAIYccYtr/t7/9DWPHjsWhhx6Kr371q9i2bZvvczMyTgghhBBCykaQtJNUZgJnb2+vaX88Hkc8Hg9tbADw9ttv4ze/+Q1WrVpl2t/e3o6VK1fiyCOPxI4dO3Dttdfi5JNPxksvvYSmpibP52dknBBCCCGElA0V/idxqpm+48eP11NKWlpasGTJktDHd99992HYsGGYMWOGaf8ZZ5yBc845ByeccAI6Ozvx2GOP4YMPPsDPf/5zX+dnZJwQQgghhJSNYKUN0+23b9+O5uZmfX/YUXEhBFasWIHzzz8fsVjMte2wYcNwxBFHYMuWLb6uwcg4IYQQQggpG9qiP343AGhubjZtYcv4U089hS1btmDu3Ll52+7btw+vvfYaxowZ4+salHFCCCGEEFLT7Nu3D5s2bcKmTZsAAFu3bsWmTZv0CZcLFizArFmzcvrde++9aG9vx3HHHZdz7KqrrsJTTz2F119/Hc8++yz++Z//GYqi4LzzzvM1NqapEEIIIYSQsqFCggp/i/j4bf/888/jtNNO0+93dXUBAGbPno2VK1dix44dOZVQ9uzZg//+7//G7bffbnvON998E+eddx7effddjBo1CieddBL+8Ic/YNSoUb7GRhknhBBCCCFlw5h24qePH0499VQIIRyPr1y5MmdfS0sL9u/f79hn9erVvsbgBGWcEEIIIYSUjWClDWsn05oyTgghhBBCyoYqJKjCZ5qKz/aVTNV9rVi2bBkmTJiARCKB9vZ2rF+/3lO/1atXQ5KknBqRhBBCCCGkfKiZyLifzW8pxEqmqh7Jgw8+iK6uLixatAgbN27EpEmT0NnZiXfeece13+uvv46rrroKJ598colGSgghhBBCvJBv2XunrVaoqkdy2223Yd68eZgzZw6OOeYYLF++HPX19VixYoVjn1Qqha9+9au49tprceihh5ZwtIQQQgghhLhTNTI+MDCADRs2oKOjQ98nyzI6Ojqwbt06x37XXXcdRo8e7alYOyGEEEIIKS0pSIG2WqFqJnDu3r0bqVQKra2tpv2tra149dVXbfs888wzuPfee/UC717o7+9Hf3+/fr+3tzfQeEvBZ+Rzyj0EQgghhDjwhPpQuYdQFQRJO6mlNJWqkXG/7N27F+effz7uvvtujBw50nO/JUuW4Nprry3iyAghhBDiCUmGJEuAJGfuZm7LmaioJEGSZf02AECW07dlCVKmn7G9/tPaz3TMus98Tae2Qh+DoS/JSwrwHelOFWcoZaFqZHzkyJFQFAU7d+407d+5cyfa2tpy2r/22mt4/fXXcdZZZ+n7VFUFAEQiEWzevBmHHXZYTr8FCxboqzIB6cj4+PHjw3oYhBBCSO1hkGaTMGeEVNIF2SDAGYmVnETX2NYqzLJkFmVtn7GNsZ/dcSArz4A5cde43+G2p77EE4yMVwmxWAxTpkxBd3e3Xp5QVVV0d3dj/vz5Oe2POuoovPjii6Z93/nOd7B3717cfvvtjoIdj8cRj8dDHz8hhBBSbUiK4k2wDZJsikY7RaILlWs5V5BNUWlrP0tbU3sNq9tZj1vu++5PHCnFCpyVTNXIOAB0dXVh9uzZmDp1KqZNm4alS5eir68Pc+bMAQDMmjUL48aNw5IlS5BIJHDccceZ+g8bNgwAcvYTQgghNY8kQ1KUtFhn7msSnCPWYUu1XSqIdsxwP69Q28l0noh2EGnO6QPAMYvCQbpryBVJkakqGZ85cyZ27dqFhQsXoqenB5MnT8aaNWv0SZ3btm2DLPO3nxBCyNDEKZJtyqs2CHd6lzniDUXJts2XRw2UTrbzibaDlOeItZ07BxFqpz4MiPtGQILqM2dcsJpK+Zg/f75tWgoArF271rXvypUrwx8QIYQQUmRcJTvfxEVrHrZT2ofTZEW7tsZjKKJkexFsyaGNpwh4zi7X9rltnA/ZRteJLUxTIYQQQkjFIEWiniPbeaU7X2Tbh3S7CrLk0NbY3mtk2066Awh3jmg7yXEeZ/Ys1bXjhiVHFRJUn/+l4Ld9JUMZJ4QQQkqFMW/ba4TbWmnEmrNtJ892kyDLna/tJttAVoqdZNtJtF3SUGxFOo80B4po144XloUUZKR8fpvx276SoYwTQgghxSCTVqKlmHiKcvsR7xLkbufNyXYSa6/ibbyWX/G2CLDbRE1HwfYo0YGCsExT8Qwj44QQQggpCElR7KXbq3CHnMftJS/bU/Q6n2y7tQkj0i3ZnNdyfSfJt3U1F0H27Xa144KkzFDGCSGEEI845nNLUrYKiVG6/eRx54tw+41u27W1a2+5LSQpvImTdsLtKM+5Y/Qk2o5SDzOydYcHAadwlwQVMlSfaSd+21cylHFCCCHEAVf5ti52oyje5DvsaLefSLfltlvud1Dx9hXl9irdXoRbtmljkemiR8tJIFJCQsrnk+23fSVDGSeEEDLkyVs60C7VxCTZMuyWeneNeIcl3V5zufNFu+3E1zb32+56MBw3X89OtK2SbSvYsvMxuzF4bWf7OKzXdsCpmp6kuvcj7jBnnBBCCBlq5Ktqkq98oN3ESv14yOkmPvLAfYu3sa+bePuRbpfIdk7E2irc+WTdr4zLufty2hrOYSWI70nCf5+hjhAyVJ91wwXrjBNCCCHVQaCoN2AWb2tE2y7yrVgEHTAJecHS7SENxfMEyhzJNYu0m2gXJNg2fWzPZ5XpQJLunkPuml+OYCUOJWEwcUq5Z1KQkPKZoO+3fSVDGSeEEFKzyLGYv7SToPnedpIOy2TIYsq3VVLdcrldxFtIkjfxNQq3k2zbnCdHnmWP7WzP53Tf5ouD4bj1tt90FSG5Rb8lSngAVOE/7UStoeeZMk4IIaQm8CXeQG7k2494G9tlbnsW73zS7VSj20m4XWTba4TbJLuKSxvjfot4G6PYtset/W1v213bMDYXEc+J0Dv0ASy+bHXAPCJuJUfMheUnIXmgjBNCCKk+/OZ82024DBL5BnIFHEjLrvV88BH1lm32uch3PvH2It2QLRIumdtY7+tjcJNt07lsRNty3L199rZthByG47C0sTuHtT1s5NpOxF3k3AQlPDBqgJxxv+0rGco4IYSQiidwtRPtGGCWamvOt9fIN2AWa8nhuAcpt5Nu22i3VbiNUW4HMdZFW++nXdPhvuzWzu6nBKuAO7bXHqNV2A37jWPPJ+F2km0XFXcqb2hqD4c2XgU858SWn8QTKiSoPp90v+0rGco4IYSQiiNHvhXZVrwB2Mu3l6i38bgx6p0vpcQpqu3Sz5pm4hjlNgh3TnTbGKWWXcTaq2jL9pKt53G7yLhpH5zb20XC7SLottJtve8k4fo+kd0Hcz/TfuMxx/sBbVq7CGXcF6wzTgghhFQInhbZcYp+W4U5DAG3i3DbCLq+36t8m1I2DPJdgHgLxbLPTrZzrmM5nnMeB9k2Cr7dcQf5Tj9eyzHYtJOEpwh5TlsjpvML037b2277vMJqKoEoRZrK008/jVtuuQUbNmzAjh078Itf/AIzZsxwbL927VqcdtppOft37NiBtrY2/f6yZctwyy23oKenB5MmTcIPfvADTJs2zdfYKOOEEELKSqCKJ0Hl2y3f29gO9vncTv2E4dxC1tpp57PJ55YAIecRbtlyWxNg2Tm6LWSDaMuwb6edB5b72nGrjDuIdU4/kzwL+35Osg2tjTCcw9xe258j2IaXzm5/Vq4tdmyRbsktGu4m6JZuooYitqVCRYBFf3x+a+rr68OkSZNw4YUX4gtf+ILnfps3b0Zzc7N+f/To0frtBx98EF1dXVi+fDna29uxdOlSdHZ2YvPmzaZ2+aCME0IIKRtyPOGYgpK35KBT7rdL9NuTgNtExq377SLfQh8rTH2y8myRb6Ns61JsL95G8c3elgztsj+NIm4SYhsht57T2i5HtmUXybZGqPXjwl6urWJteLn0PjD0QaYNDOKsv4R5xNzYxnjc5lgOPkQcoIxXKmeccQbOOOMM3/1Gjx6NYcOG2R677bbbMG/ePMyZMwcAsHz5cjz66KNYsWIFrr76as/XoIwTQggpCaYUFKOAK4q3qide5dtjtRPXyLedeLukmwhJgohYpFsxRrDthdt0P0eqM9czHMuRb6tAW2XcoX2unGeFOR0ZN0eubaXaKNRGmTa2Q1aqJb1NVpz19pnjkqWtfhv63ex9/WfmZczZby/fskW83URczpNrYhedpYz7R8D/BE6Rad/b22vaH4/HEY/HQxvb5MmT0d/fj+OOOw6LFy/Gpz71KQDAwMAANmzYgAULFuhtZVlGR0cH1q1b5+salHFCCCHFwTABU8oIt10aivFYjoADWcE2Cri230v0203ArZFv432jgNulnVij3hGLfOtCLOmSKwzHdCG2irdRpI0RbjvBtmlrTVPJtjcItx4JF+Y+kgAyMm7cJ8nZfUbJ1gTaTrDtjwmTVBuFWjZJubtoa5JsJ9l6mzzSbT2u7w+Y7G2USb8pF0MdVQRIU8m0Hz9+vGn/okWLsHjx4oLHNGbMGCxfvhxTp05Ff38/7rnnHpx66ql47rnn8PGPfxy7d+9GKpVCa2urqV9rayteffVVX9eijBNCCAkdPQ/cSwqKopjF2yrfgFnW/cq3Je0kb+TbIN5ClhzzvIViaC8DqiKlhdsm+q2lseiSbIqMmwXcKOUmYTdGt00SLnKj4nKmsogm0Ub5Noi1JJsFW5KFvs8o1bJskG/Yi7RkEPD0y5NtI0PkCHdOmxzRtux3+Gns43jcRrydZNzunHY4RXIp4v4pZALn9u3bTTndYUXFjzzySBx55JH6/RNPPBGvvfYavv/97+OnP/1pKNfQoIwTQggpGFMpwoxce56IaSlVmJOCYhMhDyzgDvINGCLYkpRehVKyRrYli2hnq5gIGVAj1oh3rpBbI9+mfXZpKDn7Ra5064KeEWttn2wV7vRPWTILtySr2Qi2jVg73bcTbP2YScKzcizb3YY36bYXbjVnn1NbK8a+YVFLC9GUkkIi483NzSYZLybTpk3DM888AwAYOXIkFEXBzp07TW127txpqrbiBco4IYQQ/0gyINIyI8cT5hQULdLtVgXFKNiKDNf8byf5tksxsWlrlW/jREuhSLlRb0Uyy7R2THGS7Mz9iDUS7iLcBik3HxPZqHgm4m2KbsvZyLaUua+JtiRnhDlzXM7ctxNsRVYd5TqqpGyFWrZuyJVrYx+rMHuVasW2nb10A4BiiWBbJdtJxhWH/YVQS7WvS0m1LPqzadMmjBkzBgAQi8UwZcoUdHd36yUSVVVFd3c35s+f7+u8lHFCCCG+sNYCl3T5zhMBN0q0Uby1/sZjQE5EW9uXs99GyrVSg6a0EznbLzuhUtInWppFW5P1TPqJnXwrZrlWI7nCrSrwL96yUbwN+doWyZZkQ1Rb2wdAzuwzCrciq6ZItnZM1iU8K9ORTFs7+QbMAu0m31apthNqxeWY8Rzp87hLtpLnePac4UfEraQYIa849u3bhy1btuj3t27dik2bNmHEiBE4+OCDsWDBArz11lv4yU9+AgBYunQpJk6ciGOPPRYHDhzAPffcg9/97nd4/PHH9XN0dXVh9uzZmDp1KqZNm4alS5eir69Pr67iFco4IYQQdzITMEVyEHJdnb8UFDsBt7azRsDdUk88pJ1okW8939sq3nZR70h2n2qJgMOwDzKg2om2kpVxk3ArwrAvkz6i2Mi2LPTotpZaIslG8c6mgCiyMAm3jPQ+Taiz+1TbSLb9vqxYRyTVJNTpNubotmK3T5dva1vz/XTbzD7k7ku3tRdmaxTcitfUk3znCUKqDJHaWqGQNBWvPP/886ZFfLq6ugAAs2fPxsqVK7Fjxw5s27ZNPz4wMIArr7wSb731Furr63HCCSfgySefNJ1j5syZ2LVrFxYuXIienh5MnjwZa9asyZnUmQ9JCBH+b2QN0dvbi5aWFuzZs6dkOUle+Yx8TrmHQAipYfQIuDUH3K4UoRf5BuwFXJH8ybcu2/YpJ/p+q3hrsm0Qb1PqiZb7nZFso3QLi3ibIuVKWrRFNPPTGt1WMpKtp5cIs2xnRNsk2LIwibYWyY5ImZ8G0Y5ITtKt5tw3yrVVrLX7USllEumsbJsF2hphtopwvvQRuzZubc398kt3vgmahZwbAFL6xARnLj3y94HGUEwqyWu0sZyxZh6iDTFffQf7BvCb0++uiMdRKIyME0II0dEnYlpLEXrNA/cYAQcMaSWZ/n4F3BT51uXaRsCV7H496m1MPcnIdk7EW7GRb31fRrqVTHqJIoCoJdJtFW6LeGsRbkUWUKSsiCuy6km6tQh2PuE2ynb6mJoj206i7SbYweTbXnTdxNmLHHtNPQlattCOKFKhnWuoU4rIeCVDGSeEEJKOghsFXFHyp6AYq6DoUW0fKSjaPsV825R2oi2cI0u2ke+0XBvEOxP1Vk0iLWVyv61C7XBbvy9Mt40pJkIRkBSRTS9RBORIOqItKxlB1qLbmZ9KJoVEk20tjUTL0c4Vb9V0W5PsiJxyFGwvkWw9BcVGrO1yuPPlYtvJsp0chynCpcj7JqWDMk4IIWRIYhsFd5Jwo3wXkgNuFwFXJAhZzi/gipQb+VY00TammUhmwc5EvW3lW3GWbz3dRNFSTQAoaflWFFWfQCkraYGORZOexFuLcGclPKVHtyNyCookEMmkimjSHZVVz8JdiGxb87yN19DbFCHf2gple2gh4L86Si3lWFPGCSFkiOEpCq5YJmgWKuBKroxrAi4isp5qoqWdqIo551uPgEfMsq1LuJLN7c4Rcdv9IjfirUl3JvotySqkTL53Wr41wU7f16Q7qqSgyCqicipHtiOSqsu2FuGOWn9aRBuAnjpiF9W2CnXuffsJlKYJkgWoTK2KsrVaCyGlgjJOCCFDgUwE3DQh0yrhxhKD+STckqKSLwIuZDk3/1vJ5HpHpBwBV40R8EhWulVj/nfEHN1WFbOEm2U8Ld/Qo9/28i0rmTSTTI63LKuIKOlNkdIiHlFSttHuhJLM3E8fj8opU1qJFvW2RriN0W0n6bYTbifZDiratSrZXlFtJmRS0EsD01QIIYTUNHIs5l/Atci41wi4LNumn9jKt5Z+kkkzUaPZY2nBlnSxNkW+IzbibdwXyY16p6XbkOOtCFvpVhQVESWVlu1MeklUSUGR1LRs61HubLRbi4BH5RQiUgp1yqAn2bYT7XySHUSwq1Wuw6zRXehzYCfoJHwo44QQQmoOSVHSom0sS2isiGIn4bKcGwG3SLgxwp13EqbsIOB6NDwt3qmYZCvgmlxbI96miHhGtNUIslHviDXqnUkzUdI/jfKdjnhnU02M8h0zpJ0Yo91aPndUSplSTKJSylG8cwXcLN2FpJBUqnRXwsI3hY6hUp/bWoMyTgghpGaQY7G0SCtyfgG3VkMxCrhBzB0j4NrtfPKtTbw0VD9Ro9kIeFrGM5IdsUa5rbczUe9IZoJlxCDdunBnot6ZvG4tzSSiRbsz0h3LpJtk5Tupi3dU/5mynTypRavthNso23b52oA34a4EESynUIe13Hkh+fGV8IViKEAZJ4QQUr1oueCKkhVwScouzONFwPWUFEMKinEhHlPKiZydfOkm3xHJMNEynYqiR74jhvzvCJCKGUTcsN9WvCOqHu2WM7eVSFa6FSUTvVayaSZG6U5Hu9P340oSESkr3tYUEy3nO59wA7nRbSC/BJZatosplmGJczFQIZWkAgwJjhAShE+59tu+kqGME0JIlWKsiqJLt3F1TGNFFKuEGyXdLgXFKQJuEHBNxlXjREw93cSQ/61NwIwYRNsQ8U4lLJHvSCbfW0s3iaQj35KSzvVWIumod8SQahKV07e1iiZ6iokh4m0n3xEpFSja7beOdinEO2zZrmTB9gNFvPJRIfn+fauV30+AMk4IIVWHHgmPRLITMhXZXcBlxRANz07QNOZ750TADfKtT8B0k+9M9Nso4UJBOiVFzka/9f2RdNQ7lchMtIyYxVuJpNNNFC3iHUnLdiySLSWoRb31n3ISETmFuJw0pZrYpZhE5aRv4Qac5S5s4Q5Lrv1Ii/GxhS2xpZYnSjipFijjhBBSRZhywqNRfxJuzAM3CrgWEZdlg3jL5hSUSPanLt9atDuSiX5r6Sc5wm2+r0Yt0e+Emq5yoqjplJNM2kk0ko56x5RUTtQ7K+FJxDLiHc+IuJ18p/flRrujmci4X/EGwpHvQoU7n+DKEBUjpZUyDlJ5MGecEEJI5WKXE65VSIlGsvKdR8BzUlAcIuDpn3J20qVFvvWJlyYJN6SeRG2i39G0dKsK0tHvqAopki4xqEQEYvFBk3jHIkk919sk3nJSTzWJy8l0aopNmoldxDuodAcV7qCSHTR6TNEtL5Uw2baaGeo541U3TXjZsmWYMGECEokE2tvbsX79ese2Dz/8MKZOnYphw4ahoaEBkydPxk9/+tMSjpYQQoIjKQrkaASStmmpKZH0fUQyW0bU9fuRTEnDiAIRUSAiMhBJS7yIWraYAlXb4gpSCQWpuIxUQkayTkYyISNZL2OwXsZgg5T92SBhsAFI1gPJBqRvNxpuNwgkm1QkG1WkGlNINaaApiSkxkFEGgYRaxhAXWM/Ghs/xLCGDzGifj8Oqu/D8Lr9GB7/ECMT+zEy0YcR8f04KNaH0fG9aE3sxYhoHw6K9qE5cgCNygHUywOoVwbQpKTvN2X2NSoH0CQfQJPyIRrkfsSkZE75QOOmoUiqafNCSsg5mxtafqzd5pdKinzXEtbfg3wbKQwtMu53qxWqKjL+4IMPoqurC8uXL0d7ezuWLl2Kzs5ObN68GaNHj85pP2LECHz729/GUUcdhVgshl//+teYM2cORo8ejc7OzjI8AkII8Y4Uj2drhVsnaUpSWrq1KLheRUU21/3WI98OEXCt6omW/x2xRr0lPdKtb1r0O2pNPclsUQEpkk49UaIqlEgqm/MdSUe640oSUSWddtIYGchJN9Ei4FrUW6vhrUW7s2kn6eMxKWn7HIYR8fYa5S5lTnQlCjillARlqEfGq0rGb7vtNsybNw9z5swBACxfvhyPPvooVqxYgauvvjqn/amnnmq6f/nll+O+++7DM888QxknhFQcenUURTZPzLRKuJb7HYnYS7ip+olBwiOyIe/bkIYSsZHwqJSTfqJt+qTMiNAlHBEBRNMCHskj4FrqSUxOISYn0RDpTwu4nMzmezuknESlZKCKJl5E0U9qSS1VcrCDYk1KiQgQ6aaMl4GBgQFs2LABCxYs0PfJsoyOjg6sW7cub38hBH73u99h8+bNuOmmm4o5VEII8Y1cV2deMVNLPdHLFcrmfHBZTqehaBMxtTxwWZt4KWcnXRoj4NFM7ncknf+tTbpMRSVTvrcWAc9GwYVe/UREM3nf0fSky0g0BSWSLjWoiXciOoi4kjSJdyQj33HDFs3kf/uNegedWOlFuCtBtMOOfFeSXMuonLEQUglUjYzv3r0bqVQKra2tpv2tra149dVXHfvt2bMH48aNQ39/PxRFwQ9/+EN85jOfcWzf39+P/v5+/X5vb2/hgyeEEBtyIuHaRE1NtvWIuGyWcK0sYVTJVkaR5WzkW7EKuKxPwkzFLFHwKMyTMGMw74tmo98iKiBFVUiZiieRaMpWwBNKMlv1JJNyEjOIty7fmZrecXkw70TLIPI91MS72MJNiSbFQgAQPv8MKi9RKzhVI+NBaWpqwqZNm7Bv3z50d3ejq6sLhx56aE4Ki8aSJUtw7bXXlnaQhJAhhRyLZcXbmpaibdaouCbgmnjLclrGlWz1EzXiEgGPZtNPUsYIeNS6iZzotxxJpVNPFBWxaDIr35FBJCIZ+TaUGdSj3rKKuDSYE/m21vVOyIPp58Xjx6tROr1OmCw3hQp3MUSbck0qBRUSJC76U/mMHDkSiqJg586dpv07d+5EW1ubYz9ZlnH44YcDACZPnoxXXnkFS5YscZTxBQsWoKurS7/f29uL8ePHF/4ACCFDGscouCJnhVuLgrsIOBQt91uGGldM+d9qTDblfKeihui3ddOj4Ab5jqWj3nJURTSaQsSY7x1N6lHv+shgJuUkpdf21gQ8IQ861vQ2pp0UIqdW+a6UD+VCHlNYsl1qwVakWopPknLBCZxVQiwWw5QpU9Dd3Y0ZM2YAAFRVRXd3N+bPn+/5PKqqmtJQrMTjccTj8UKHSwghOnIsll0t0yrhkiH1xCDhpoV4InI2F1zJpqOkEkom4m2OgGv537p0a7djmoSLbPqJIfUkGssKeDySRCySQiwj34lIOuodRMD1Sigh/cdypcg3EEzAwxDvYko3BZuUGlVIkLjoT3XQ1dWF2bNnY+rUqZg2bRqWLl2Kvr4+vbrKrFmzMG7cOCxZsgRAOuVk6tSpOOyww9Df34/HHnsMP/3pT/GjH/2onA+DEDIE0CLhav8BSHV1lgopirOAK5bJmBnxVqNyNg88kl56PlknZyPgBgFPRQERBVKxdDpKKiYgYgJqJP1TiqqQoylEoyoikRRi0STikfRiO3FNvpUkEsqgSbqNaSd6aonLUvKFyHc+4U4JuaSTEv0+lkLHFrZsV4NgKz4fc6r6lkohDggRIGe88n+lPVNVMj5z5kzs2rULCxcuRE9PDyZPnow1a9bokzq3bdsGWc7+cfb19eHSSy/Fm2++ibq6Ohx11FH42c9+hpkzZ5brIRBChgDGSLhSl8hNR9EmZCqZiHdEMUm4ULK54OlqJ2npFlFZj4KnYsBgvZyTfpIV8GwEXMRUICogRVNQoumJl/Fo0hT9TkQ0CR9EnTKIiJRK1/xWkohLg4jKKcQzEXAAoQu4F/m221csIS9lxDtM8S61dPsV6Fq5NiFhIglRS98twqe3txctLS3Ys2cPmpubyz0cE5+Rzyn3EAghFpTGxnQqSiTiKuBQrOkn6ZUyhS7fmUmZUUlPRUnFpEzkW4IaAwbrM6knWgWUWCb1JCaAmApEVCixlEm+jbnfiUgS9cpAWr7llF7tRPuZMMi3Jt7p28Hl2026gywhX3gEujTSHYZwF1u0KbfhM/PwP5Z7CDlUktdoYzlm9Teh1PtLEU7t78dfzr25Ih5HoVRVZJwQQioRU064tmy9lo6irZJpJ+GZiZhWCVcj6WooqYyIaxKevp2Ngg82ZgQ8BqhRNZ2CElPT1U8y+d+JWFrCjZVPEko2Am4UcGP0OxsBL0zAg0S8/RAkOl6KlJNC5btY4k3hJpUIJ3ASQggJhBxPQIpGDJMz5eyqmJbbtgKuR78z6SeZaiipGDICnpHvWDoPXIuCaykoakMKiKmQY+n872g0iVg0hbrooK18NyoDegRcl3BNxDP1vYG0eANwlW+/kygLlW4jfuTYj3iXOs0kbOEuh2jLFbSYEKlehvoETs5+IIQQP0hyOhIOQErEgVgUUiwKxKJALAZEI6bbIqZAxCIQsQjURARqPAK1LoJUXQTJOgXJhgiS9QoGGxUMNsgYaJIw0KT9BAabkP7ZDAw2Cww2CySbU0g1J6G0DCDe3I/GpgMY1rgfI5v60Nq4F6Pr92F0XR9GJ/ahNdGLtkQvRsX2oTnyIVpM2340yQeQkAb1ZeijUgoxKYmYlFma3iCzKiR9c0Or/W3cCkWRVH1zQ4YwbWGd13wN1bR5QZFEzlYICtScLWzkTG14t42QMNAmcPrd/PD000/jrLPOwtixYyFJEh555BHX9g8//DA+85nPYNSoUWhubsb06dPx29/+1tRm8eLFkCTJtB111FE+Hz0j44QQ4hljJDzS0pxZkj5ijoJrS9TrOeCyXo5QjaQnYqox2ZB+YoiE61Fvw894ZhJmXNWj4LFYOge8ue4A6qMDevWTmJxEnTKIukweuFZSUIt8p6Phg/rjMQqcNQKeT7rDjHTb4VWOgywU5JUgEe+wot3FjHJTokmlkZZrv2kq/q7R19eHSZMm4cILL8QXvvCFvO2ffvppfOYzn8ENN9yAYcOG4cc//jHOOussPPfcc/jYxz6mtzv22GPx5JNP6vcjEf9qTRknhBAXpEg68q3u3w+poS49MVMX70g6J1xLQ7EIuBqVMxVQZKixdA54Kp7NA0/Fs+knqXhWvtWMgCOuQoqmEImnEI0mTfnfjdEBDIt9iDolXWqwThlAvTKgLzGfFfCknu9tRJNTGUKXbqt8exVua4m5ICLpRZaLmXJSLvEOW7rLJdpDIRedpRSLRylyxs844wycccYZntsvXbrUdP+GG27AL3/5S/zP//yPScYjkYjr4pNeoIwTQogDpkh4Q30mCm7IDY8orhKeiqYlXI3KevRbk3A9+h3Xfop0KcK4CsTTdcAjsRQSsUEkokkkooOoiwyiPjKI+sgAYnISw6Ifol4e0Jeb1xbeseaAa9gJqlHAg8q36Ro+pcxL2klY5zKft/R53tUq3kNBtL3A52Foo6oq9u7dixEjRpj2/+1vf8PYsWORSCQwffp0LFmyBAcffLCvc1PGCSHEBqWlJRMFV8ypKIqSrguuyXdUyVZBiSmZiZgSUjEZqbg2ATMTBTdFwDUBVyHi6fQTLQJeF89MvszId2O0H4lMBDyeSUVJyIOolwd0+dZWusx5HA7Clk+8/UQB/UhKWPJdyeIdprQVU7gpl6RSEJnNbx8gXR7RSLFWUr/11luxb98+fPnLX9b3tbe3Y+XKlTjyyCOxY8cOXHvttTj55JPx0ksvoampyfO5KeOEEJJBj4THotl0FF3G80h4Jhc8HfnO5IHrP7MCrqehaBHwmDkCXh8dQF00EwFXBtAQ6TctvlOvDOSNgGsYhdV20ZwC/9s9LAmvFAEvd9S7WOJdbumuhtU/ASBVQ9U5qo1C0lTGjx9v2r9o0SIsXrw4rKEBAFatWoVrr70Wv/zlLzF69Gh9vzHt5YQTTkB7ezsOOeQQ/PznP8fcuXM9n58yTggZ8sh1dZBisWwkPBo1iLhikG8lnQNuEPBUPD0ZM5kw5IHHgFQik4KiC3g6Ai7FVETiScRi6RzwuuggGmP9SChJ1EcG0BAZQJ08iLiSRERKISEP6rW/7Va9dMMo4GHmu3qRuzAi4JUq34XIbdjCXSrRrhahLoSh8BgrlgJC49u3bzct+hN2VHz16tW46KKL8NBDD6Gjo8O17bBhw3DEEUdgy5Ytvq5BGSeEDFnkeAJq/wFIdYm0iEej2Uh4NJLJAVfSIm4j4enIdyYSrot3WsRTMUBNpCUccRVSLIVYPIl4LJP/HRtEU7Q/Uwt8EA16FDy7CI+2AqYmcJqEx6RkqGkmXqk0CS+VgFeKfJdCvCmkpCwEiIwj0765ubloK3A+8MADuPDCC7F69WqceeaZedvv27cPr732Gs4//3xf16GME0KGHHokPBaFHBmWrg1uFPBIxFnAMxHwZFzSpTuV0CTckIKSSC9Fn4glkYgPplNPogNojA6gPjKgL8JTl0k7sa5+aV350oiXiLd1QQzZIFnGY3Ie+congJUs36WOfIch3sUS7kqQ7EJXJS0VKqumlJwgdcP9tt+3b58pYr1161Zs2rQJI0aMwMEHH4wFCxbgrbfewk9+8hMA6dSU2bNn4/bbb0d7ezt6enoAAHV1dWhpaQEAXHXVVTjrrLNwyCGH4O2338aiRYugKArOO+88X2OjjBNChhxSXQJSNJpNR8kszoNMFDwt3wYJT6Sj4Ml4NgqeNKahZCRcTaiQ4ik9DaUhPoB4NIkGTcIz9b+1OuAJeRAyBKJyEgkpiaiU1MfolIbiVb7zHaeEW/qUScCLId+lFO9qEWw/+H1MlPfq4Pnnn8dpp52m3+/q6gIAzJ49GytXrsSOHTuwbds2/fhdd92FZDKJyy67DJdddpm+X2sPAG+++SbOO+88vPvuuxg1ahROOukk/OEPf8CoUaN8jY0yTggZEhij4aivA6JRiGgEiGYEPKqkF+OJKVDj2YmYyYT205CCEgdSdZkUlIQKOZ5ENJ6ehFkfG0SdHgEfRFP0AGJyEvXKIOpkLQqehAKh1wEHzHKsCZpbqknQpaDdJLxSBLzY8l0O8Q5Tuost27Uo2MWEz1fhlKLO+KmnngrhEk7XBFtj7dq1ec+5evVqX2NwgjJOCKl5lKYmSLFYepn6aASoS2TSUBSosUhWwjNR8KyAS4YUlEw0XMsDT6Qgx9N54In4IBpiA/pqmI3RfjRFDiAqqahT0nXAjfnfWipKzjihIgXZJOFBpdtKUAkvZSlCP1JT6QIelnwXU7wrVSKDrJYahGKvIkt8ICQ9B9xXnxqBMk4IqUnkujpI8Xg6Eq6lpMSjaQFviKflOyan64EbBbxOQlKLgNelU1BSCQGRUIF4CpFEEnXxJBKxtIA3Rfv1hXgaIv16LfB6eSCdhiKpehUUOwE3MigU1+N+V7osloADpU9DKYV8lyvyXQzhLrVol0qgwyTMMVPsC6MUOeOVDGWcEFJTKI2NkOLxdBQ8I+IiHklPyoyn88EHm6K6gKej35k0lLp0FDxZJ9IR8Lp09DuaSEe/tRSU5li/qQyhVgXFWAEFgKkMIZCNcsuSyIl4+61+4iR/Q1XASynfQcU7TOkutmxXo1yXEz5fBVLIqj81AGWcEFITaJHw1AcfIDKsBYjHIGJRiFhWwtWoglRCRv8wJUfC0z8zUfCMhMcTSdTFB9AYH9AroTRG+vVJmI1KP6JyShdxOwlXhZQj3tr9oOUH/Yp4peSCp89VXgkPIuDllu+hLN7FLOdYjPKfJBilyBmvZCjjhJCqRo+Ex2NALIZI68i0gMciUOMRpOJKpiShnE5BSUjoHy4hWZeR8Lp0FRRkBLwuMYi6+CDqYwNoivWjMTKAxugBfRn6RqVfn4Rpzf22ioObdKtC9i2GdmJSzlzwcgh4Lct3MaS7HKJd7hU/vRL2OCn3JCiUcUJIVSM1NWZSUmLpdJR4NFMRRUEyoZgkPFmXnpB54KBMNZSMhCuJJOLxtIQ3xfqzpQgj/WhUBhBXkqgz5IDbTcD0VgFFdr3vJI5hRsJLmY7iVS69SKwfcaoWAQ9Tvksl3dUi2uWAz02B1FDaiV8o44SQqkJpaoKUiAOZaLhIRDNR8CjUuIJUQkGyzizgyTogWZ/JCa8TUIcPIJJIoiGRXgmzIdZvygNPV0JJr4RZr/QblqEXjqKXgpwj116pFgmvJQEPIk5BxTss6S6mcFeLSHqt3uMVFbWT6lDNME2FEEKqhMhBB6XTURJxiHgMImGWcC0KPlgnIVmfzgfXJbxehVqXLknYPGw/GjJ54LqEK9lqKAl5EHF50JOEuwl4ChKUPPJQCSIeRiS8GiS8VAJeyfJdbukOW6YLpZDxUORDhBM4CSGkclGamiDV1wGJeEbCo1ATEaQSkbSA1ylmAc9EwZMGAZcTSUTrspMx2xr2ZiZiDqAxMoC4NGhajh6AScBVYR/1TuX5MHYT8WqQ8LAi4WEKeK3Jd9jCXSrZrjSpLgd8DsJEymx++9QGlHFCSMUSGTUynY5Sl4BIRNNR8LpINhWlXjZLeENGwutUqPUpyIkUEnXpXPDGeL8+IXNUYm/OiphGCdewCng++faKnVCWulRhtYl4MSW8WgW82OJN2SQlg5FxQgipHJRhwyDV1yH59g6I0SMgEjGk6jKR8DoZg/UyBuvTAj5YD6TqgWS9SOeC16eg1CURjyfRUNevL8rTGB1AU/QAGpR0Gkq6JGHSUgkl/c7uJNxO6SjOEe7cT4qwouFufYot4bUg4KWS7zCku1jCXQmiXcklFb3AhX5IWFDGCSEVQWT0KKAuAdQloNZFgQmtGKhTkKxXMFifyQVvSAt4siEr4MgIeF1iAA3xATTG0/LdHDuAOmVQF3AtF9xpJcwUJN8TMO0j3IVLeLpPuCJeLRLuVcBrTb7Dlu5SyXa1C3UhDOXHHjqMjBNCSHmQ4wnIjQ1IvvsuMLwFaiIGtT6KVH0Eybp0Kspgg42E16u6hCcSg2iID6ApfiC9KE+03zQZU1sV0zgRU6vx7STfxui4U943RdzYvzYkvJQCHqZ8F1u8KZ2k6AgpvfntUyNQxgkhZSEydgxQXwdRH4d0+Dj010eRykTBdQFvyBXwSN0gmqwCHulHQ2ZlzHp5ALIkdAEHcvPAU5CQEkreMfqT6+oT8UqQcMCbiHuVVz8CXir5DkO8iyXclSTa+SoPVRphzSEhgBDpzW+fWoEyTggpOZFDJ0DUx6HWx5CqjyLZELGX8IashEfrBrOpKFpdcCWdCx6Xk6iXBxCVU5AhEJWT+rUUiJwPTS0iHmRxGDv8iLgbtSTitSTh5RDwsOW7XNJdbYLth1p+bCWHaSqEEFJcIgcdBDTWQzTWQa2LY6Ahms4Fb1Aw2GgW8MFGAbVeBeqTiNUNojGRLkfYGOtHc1RbmKcf9cqgXpLQKuAaqnDXLjcpt37QhhERT7cvzWTNoCJeCRIetoBXsnyHJd2lkG3KJykaTFMhhJDiETn4I0B9HdSGOJINmSi4UcIbkRbxxnQ1FNSlEKsfRH1iwFSOsDGanpBplXBNQhzzvy37rdLiNffbDr9yEqjudYA+lSritSbh5RTw4q7GWV3SHdYiS0ZUsFJKKZFEevPbp1agjBNCQicydgzQ1AC1IY6BxhiS9REMNqYFfKAxLeBJg4BLmSh4gybghnKEdfIg4koSkcxETAC6hKeEXFB5sUInYYZFmAv6FEPEK0HCiyHgpZDvQsU7/EWBymMwxRDmYlJt4yXVDWWcEBIqkY8eBrUxjlRjHIMNNhLeCCQbVaTqVSj1g6irS0fBm2KZKHi0Hw2RAdTJg6hT0pMxo1IKiqSaxEaTcLclqctRSznIpE2/1JKIe8GLiFeShBfyexemfJdKvCmupGCYM04IIcGJjB4FNDVCbaxDqjmO/RkBH2iSMNhkTENRoTakEKlPImGoCd6sCbgykK4HnomCawJuxCje+dJP7DC2KSQ9xW9bN/xGxStFxIsdDQ87Eu5Hcksl3+EsClQcI6Fgk5LCnHFCCAlG5LCJaQlvimGwMYrBpoyEN0oYaAIGG4FUgwrRkIJSP4h6h5rgRglXoJoE1WsaSkrIFVWmrZQMRRGvRgmvVPmuJPGu5r9hrshZAIyME0KIN5QRwyENa4HaVJ+OgjdGMNCsYKBRxmATMNCUzgNPNqSAhkxN8LpsOUKtJnidkk5BicpqjoCrQoKaiXioLrXArRHlnImZcD9eSvxWUPETFQ9bxAtJS6kUCff6WleyfIcp3aWU7WqW6UIZyo+9YCjjhBDijcHjJiDZFMVAk4z+prSADzalyxGmGjPlCOsHMSwzEbMhOqCXI6zLRMA1AQeyMpSCjFSeN9awlwu3IwWp6ipJOBG2iBcSDS9Uwssl4KWS7zB+54ol3BRMUhIo44QQYo8ybBikEcOgttQj2ZxA37g4BpoyKShNwGCTgNqQglSXRLw+Ww2lITqAxshAjoRbhSEFWY+CW3Gru50zTg/CYE1jUYVc9CopQQgrKu5EMaStGkS80iS8UAEvxutYieJdLV+OuRonKQTKOCHEFvnjx2KgJYHB5gj6m2UMNEnYPxYYbEpPxFTqkkjUZSdiajngCWUwnQMuJyFLAookTOKgQobq8PkaVDCqMV/cz5cNvzhP8vSfI16utJT8Nc/DjYT7EXD/1ViCv9ZhSncp/0aqRaLDYqg93tAZ4hM4Pc82ePvtt4s5Ds8sW7YMEyZMQCKRQHt7O9avX+/Y9u6778bJJ5+M4cOHY/jw4ejo6HBtT8hQJnLQQYgcfiikqccj9ekp6DukEXsPjqH3YAX7xkvYNx4YGDsAefQBNI/sw+gRezGmpRdjm3rRVr8XIxP70BL9EI3KAOqUQV3EgYyAQ8agUJASkr4VE7uSh0NhglWgKG0Rvhi410gXjteUkfs/KLn9vUXCPa0GCuFZxBVJ9fElQOibH2T9ryX/8+B6/cxYjVuYGB+f3UaIH7RFf/xutYLnyPixxx6LZcuW4Stf+Uoxx+PKgw8+iK6uLixfvhzt7e1YunQpOjs7sXnzZowePTqn/dq1a3HeeefhxBNPRCKRwE033YTPfvazePnllzFu3LgyPAJCKotI62gkd74DTJ+ED5tjGGhW0N8iY6AJGGhJT8ZMNaT0XPBxzfvQGOtHfWQwk4IymClDqOpRUlkShgmYuUJsrpQiBRZBVUg5VVfsSiG6iZZdqopd3rjXdsUiDJHymyceNCIeNBpeynSUYkXBg/w+FBr5Lka0uxpkupj/s+QFp/Q6EpAhnjPuOUz0H//xH/jXf/1XnHPOOXjvvfeKOSZHbrvtNsybNw9z5szBMcccg+XLl6O+vh4rVqywbX///ffj0ksvxeTJk3HUUUfhnnvugaqq6O7uLvHICaksIuPGQjnuSAx+dCz6/2ka9k6sx56JUeyZIGPvwUDfISoOjBuEaDuAxKj9OGjkXrQN60Vbw16MTuzDQfE+NEUOICYnTSKuChlJVYEqZNNmxOuHmLVdysPblZfIdy1Hx4Okp5SSYoq4n0i4F4JEwb1SSOQ77Gh3OaPasiQCb+WmWsZJsjz99NM466yzMHbsWEiShEceeSRvn7Vr1+LjH/844vE4Dj/8cKxcuTKnjZ+MDSc8fypdeuml+POf/4x3330XxxxzDP7nf/7H98UKYWBgABs2bEBHR4e+T5ZldHR0YN26dZ7OsX//fgwODmLEiBHFGiYhFYnS0oLIoRMgTzkOANDbPh7vTh2Bd6bU452PRbB7EvD+8Sn0HTsAcWQfGib0YsxH3sfEMbvx0VG7cPiwdzGx6T2Miu9FUzQt4doHTQoSBoWSTkGBZNo0rELuhAp3cbcKuRexd1uh08/YSkExJxnaXi/EqLibDLulpeQjn3R6kXAtFSXfc+ZVcv3Ka9DUkzBSTUqRTlKtQk0qBwkB0lR8XqOvrw+TJk3CsmXLPLXfunUrzjzzTJx22mnYtGkTrrjiClx00UX47W9/q7fRMjYWLVqEjRs3YtKkSejs7MQ777zja2y+JnBOnDgRv/vd73DnnXfiC1/4Ao4++mhEIuZTbNy40dcAvLJ7926kUim0traa9re2tuLVV1/1dI5///d/x9ixY01Cb6W/vx/9/f36/d7e3mADJqRCiBz1USSH12OgJYaBYQr6p5yIgWZLNZT6JOJ1g6iLD6IhNoD66AAaowNIKIOIySlE5ZQuMvkEuNgfsinIJvnykq7i/xreUlWqgaD1xG3PVYR8dPdJpYVNEPUTBffUzmf0OwiF/u4WbUVOyjOpcs444wycccYZntsvX74cEydOxPe+9z0AwNFHH41nnnkG3//+99HZ2QnAnLGh9Xn00UexYsUKXH311Z6v5buayhtvvIGHH34Yw4cPx9lnn50j45XKjTfeiNWrV2Pt2rVIJBKO7ZYsWYJrr722hCMLjvyxYyCiCsT6F8s9FFJBRCYcAvWgJgwMT2RzwFu0coQCqcYUUJ9CtG4QDfFB1MUGUR8dQEN0AIlIEjE5iZicgiypukjJEFAhmRbhMYqOX0kzCrQ2kdM42TN9flVvm76GQbh9CnmQ3HEvVFJd8rByh50XKHKOiDuPKZiIl0LCK0HAK2VRIKAyZLsUawmUEi+pdSRDBVZTWbduXU7wtrOzE1dccQWAbMbGggUL9ON+MzY0fJn03XffjSuvvBIdHR14+eWXMWrUKF8XK4SRI0dCURTs3LnTtH/nzp1oa2tz7XvrrbfixhtvxJNPPokTTjjBte2CBQvQ1dWl3+/t7cX48eODD7yICEWBkGXgxMlQIxJEVEYqKiO25o/lHhopIZGPjIN6UAvUP/0FH85ox0CzjP4WSY9+JxtUoGEQSiKFWHwQ9YkB1EUH0/KtJJFQkohkIt/p/z7OfiBqomxMOclKs73gmkTbIKtW4c0VaPNkThWyuSSitX0BQm472dMyvqCibT/ZU64o0fCbnhJ2RLyWJNxv2kkQwkwnKTaV9Htebvhc+KCACZzWDIZ4PI54PF7wkHp6emyzMXp7e/Hhhx/i/fffLzhjQ8OzjJ9++ulYv3497rzzTsyaNcvXRcIgFothypQp6O7uxowZMwBAn4w5f/58x34333wz/uM//gO//e1vMXXq1LzXCetFLAUiIgOyBKFI6Z+SBEkAA2d8ArHfUMhrHeXoI5AaXo/9LVH0D1PQP/1EfTGeZFO6CoqcSCJWl0QiNoj62CDqogOojwwipqQQk1O6hGukl6I3RL8tH95OUuYWWS5EyItJkHSWcqSqFBI9LfbEzSCyETh9o0QiXikSHsqqnCUuWUlIYAqQcWvAdNGiRVi8eHEYoyoZnmU8lUrhz3/+Mz7ykY8UczyudHV1Yfbs2Zg6dSqmTZuGpUuXoq+vT8/VmTVrFsaNG4clS5YAAG666SYsXLgQq1atwoQJE9DT0wMAaGxsRGNjY9keR1iIiAQhS1kh1/6XR5Zw4Ox2CBlQIxKklED9w8+Ve7ikACJj2iBGNKdzv4fF0N+iYKA5E/1uBJJNKlL1SciJJKKJJBpiST39JB5JR7/TqSfZyLdWgtA4wdH44S1D2KeIGITZGG02iqq1n1XI08eybZ3Pn5uyUsp0FS/R8eARdClHlrxG0MOYvGmlFOkpQSLilSThXgW81PIdpnRXkmxXSgqYH7gSZzCC1A3X2m/fvh3Nzc36/rACqm1tbbbZGM3Nzairq4OiKIEzNqx4lvEnnnjC14mLwcyZM7Fr1y4sXLgQPT09mDx5MtasWaP/F8G2bdsgy9kcrR/96EcYGBjAl770JdN5qvFbkx1CkdMCrqSj4kKCLuR6+pUA1KiEved9Ek0P/AEDZ3wC9RvfSNeWJhVNZPQoJN/ZBQDY//GD06tgNsvo19JPGgXU+hSQSCGSSCIRS6IuPohEJIlEZDCTfqIiIqt6CgqQlWHjKph2Yuwmrl6E3Ho+q7Ra5d1rhLzY6Sr58BIdr+TJnsX+n4cwz1+oiJc6El6qlTnDku9Si3c1yrVfhsJjLAoFRMabm5tNMh4W06dPx2OPPWba98QTT2D69OkAgmds2FEdsy8NzJ8/3/FBrl271nT/9ddfL/6AyoiIpCUcEiBkCUIGIEuZGkEAMsc0er86HUIG9p9+GIR8GCABw1f4m2RAioPS1ITU3r0AgGTHVAw0K5nc74+mq540Aqk6FWr9IOR4So9+x6JJXb5jSgoRKS3fgPkDOz350hABN02+FK5ibNxnbOcm5Onz5p4vqJCHnT9uJEfWQ8odD4sgixmFhd+oeJA64k6vS7FFvByR8HLJd7Gkm+JJQqMAGffKvn37sGXLFv3+1q1bsWnTJowYMQIHH3wwFixYgLfeegs/+clPAAD/9m//hjvvvBPf/OY3ceGFF+J3v/sdfv7zn+PRRx/Vz5EvY8MrVSfjJIuQtWi4pKenCNlmgrHN/5pJIt3uvYtOhFAAKQXU7U6hccsHSL20uQSjJ0Cm7OCwehxoMa9+mZ54KTICnoIcSwt4NJpCIpYW70RkEBFZ1dNPAOfVL03pJ5KwlTm7aK5Zsr0LufV85RDyYhIsnaV44yu8HF55ovjFEPFSpqQUU8ILEfCwX09KN6kFnn/+eZx22mn6fa1Yx+zZs7Fy5Urs2LED27Zt049PnDgRjz76KL7+9a/j9ttvx0c+8hHcc889ellDIH/GhlckIQT/ylzo7e1FS0sL9uzZU5T/BimEUz53c1q0ZUuKSiaHXEiAUJCOmEtZUReyts8g9NZjAKL7Bep3pVD/ei9SL/qbGUxyUU44GsmWBAabo+mUkyYZg40Z8a4XSNUJqAkVUjwFJZZCNJpCPDaIeCSFqJJCLJN24rRghlVQjIJu99PYxyzrquG297aKTVu381k/4J3a5Z7bLBo5beF2HsMxOB+zjiffeO2O27dxH7tdG7uxAebx2x/3Hr22u2a5ouKuCw4VORoeZiTcr8AGle+wxLscwl0J5RSLzfHj3yz3EHKoJK/RxjLxuv+A7FJ22g71wAFsXfjtingchcLIeDVjI+Km1BQtXQWGaLmHuSVa1HywUcIHzRG8f9gIKCefiMQHAg1v9yO2oxfJzVvyn2iIEzn4IxAtjUi2JNA3LoGBxox8NwLJeiBVL6DGVSCuQoqlEImlEI0mEYumEFPSAh6VU5AsK9alK57kTqzUXtt8wuKWhpK+bR+p1iLfTm29TuwsdoTc7bEaUz4KTfUII53Fbey1Tlh10b2eLwwRrxQJD0PAiy3fQ/X3mgSkAuuMlxLKeBUjIulfRFOqipQb4YZkaAcYBF0y3be/SLp7Kgb0tUnYNyYBIScAaTTkJBD5EIjvEUi8m0R894dQN7wU8qOsfJRjj4DamECyMYrBxkhGuiUM1gPJBiCVAAaGpSDig5CiKpSYCiWaQl0khWgkLdyKLKDIala8DR+UQkhICQmqJdXEiJaeYkxByaaVZH/atQVyK6cY8779yHsQIQeyYhCGkAfNHy9F7riXVJUg6SxeJ6IWOrkyrKi43/SUQtJSSiXhXq4FlF6+a3FxIFKDlCBnvJKhjFcxJgkHzNFx4/3Mbb2P4b5nDO0lkf4bUCNAfwvQP1xC74QopGQUysknItIHxHtVxN9PIvZBf9WvECrHE5CHNSO58x3Ik46B2hBDsiGCZL2SFu96CckGIJkAUnVAKpGOeKtRARFTIUVVxOoHoSgqopEUIooKRVKhyEKXbyPGcoN6agkEhJAg+fwgtAq54z6b3PD0bdkk5PnaAf6FHHCuRe5HyI3kE3J9vyU67lfI/R6vRCqllF0ling5JbwSBLySxLtS89ZZyjAcCiltWAtQxqsYoVVOgVnC9Vxxg6xrUXO9LXJvw+k2kLZvp4mhmWNqFFBjwEATsH+sDCHFACkGnHkipCQgDwJKP6AMAMoBIHJAILZPRaQvhcj+JOQDg5D6ByF9OAB8eACivx+p994P/Pw4ERk1EojFgHgMIh6DSESgxqNQ4wpScRmphIxknYxUXEIyISGVAJJxQI0fjlQs/RjVqICICoioChFVIUUEJCUd8ZYVgaisIqKokGUVkgQomfKCRtTMC5NC9kPPKOeSTbRbE3K7qLcx+m1MWXGNjBv3WaLkbkLu1g4Ipxa5XyH3W2FF359HyItNJaaq+PkCEUZUPGwRr6SUFD+vbRABr9TFgYxUqkiHQS0/tpLCyDipVrSJmYA54q1XVDGkruhtHIQ7UOqVnaAbjmXKnKfHFAGS0XTahnHyKKAAkgyIGCQV6S2V2VRAHjDcF8i2yWz6OIwY0nS0a6ly5guKZVO12xFh2C8gFAEoKkREQFIEIAvIigpJEZBlgYiclmZJ38xCbXoqBJBSZV26teOy0Nqbh6/JWY54uwi5Ebs0FO05ype2Aphzw7PHzULu1i59W0Uqs1/JEXhvpQ+DCrk2lvT57FM+/CwGFHZ0vJhpKMXEz5iDLndvPkd5RDyUkoge5bbU8h22dFNECQkHyng1Y0lRAbLyac0ht0q4sBpg0CEI/yKvpblIqtZXAiQBNQK9qku2VrrhzV6YxdLxc8DQ1/wlRQDa8yOnj0lK5v/G5HR7SRKQZAE5c1+WDYItZUVb/0+BzAWE4XmQpNzJlUZBt0s3cRJwr0JuPIdbGkohaStuQm5sl74dzsTOfPXP9efPY8qK6ToF1B4vNqUs01hLlFvEiyXh5V8cqPKluxhjZApKCQmQplIFv5aeoYxXMdoEzexEzcx9g4Db/jRiFXUEiJJrfxBhvG8ZI+oCmUWMMlKdEWLjbf26Dsf17xxaeyl7Ickg19BFOnMOAUiQoKra/vQmkJZgYeybEWM7yRb6sXR/o2SrkpYLnnvMq5AD0D/W3dJWnCZ22vWztgfshdytXfp2cYXcGHlPn9dbykoY6SphR8eDpKrkqwbj9gUlLPxMCPUzaTNIVNxNxiolGu5Hwktdm7yQ61biNcKgWsZZEwzxNBU5fxNSqRjTMOxyx12xE/MCkYL8Mbmdyw0h5V4rsy89jvRtoabFWm+vjzFzXEhpqTWey+NjEIYn0Hw7/dOaPmLdL0R2oqbdMbf7KnKvbV3sx3h9az/bYy7tze1y3zbs2lnbpoz7Yd/G2DfloU3uOLy1y47J/FiM48o5t6Uto2b2FPMLQDFSdcKIhueTYAWqZxFXIHxJoHb9oCKuXc/vdf2etxjXIDWECLjVCJTxakbLjc6IdXpCZ256Su5Pe4nISWcJOizDH4pR0O32Z/dJpuOOOP1Xltc/yhyBt9y1SrWNZBvlXQgpmzWTEfv0lt6nGvYZz2EVchXpPkJkqqkY+jnd1/oZx20n5Dn74HLMob1VUjUxNbU3tCulkBvPmT6v07mc+ju/DebIukvbfLKee9x6bsn1uNs4qxU/UfEgeeIyVNf89XzndI225xFgTcC9SLhfUQ0q4GGJsZtoU7hJECQRbKsVau/dfQihR8GlrITnTUtxPI+XCxY03GAYRT3nmGQv78bouLbLcl+PjuvHXaLjNjJtbWMUctM1kRvddoyCe/wm5BQldxJyu332Ee88+yA5SrZdFLqcQu6EVXCz/Q3ntbwOrrJeQdHxSpb1Yk4+DSJ+xc4NL0YUPAwBD0K5RVuWqmMjpBAq992b5EVbyt4o4bYibtonmSd1VsqbiE103PEzx+6YIcKuYUxXybaz3Ef2vibkulALg3xbhNzcRsobITfu084B2At5vgi5qb2LkKuWPm4Rb9fIuLBKuLFNrhxXQoTcLV1FE3L3dBeXYyFGx0uN1y8tpSCMqHjQiLj7uAqPhufDq9QGEfBCxLmY4l3rklvt4yflhTJe7VikWlh+2lZTKTd2qSp52+eJjufdZxPVdoyE+xdy62msEXLVdD+/kLse9yjkRlzzyW3E2drP2tZ8rtJEyJ0IIuT2fb2nq5iOFSDc1ZyqUuhqnt6v43/CZuBrFTnq61fC/Z430P8ShCzelFISCBFwqxFYTaWaMUa/DfucjttFxQMFyoSh4omU3afdFYafxn2Q4KkUoiQMf2PayYzHJMMYLNfXL2gap6SfNH06yXwRKXNCkW4rSXa3tWZS5rbIyHJue/3xG6P0WsUSmyfTtpKKXjGlsiqt5KuyYtqXeSHCWBxIloQuvQry1CoX3hYFssNYQcWtukq+6immc7ocqzXsotFeU1TcKqgUOgYvY8kXEXcjXzTcTyqKV4KKd6FUkliX8u+qklLSapEgOeC1lDNOGa9i7IJkORFx/bbkHB03CLptukumDWAj01YpDvRAjPIu5f6F6Xctx6yy7dBWEuY2WtlCU58QhVy7qCbkTqUNtXKI2qCz8orAQq499PR53FfsBGAr2U4iruEm5Onj3hYHKnbZQzchdyt16CbkRtxqj/s5FkaZQydKUd4wCJ5zqgOkp/g9l9v5gMJEPPwVOv39XhS2QFDgrp6oti+o1TbeqmQIP8WV83+eJDBCym4A7EXc0DbwF/ww35wz/8Xk+jnkcDxwukq++9YJnba3kbPfS8qKtcKKuX/uExs0ZcWIl4mdelvL9eyOubXNm55ie15/KSth41RdxbVPwMmchUTV3FJV/Ew2rWTCioo7Ue0i7jeNJGjaSdhpJay2QjwjAm41QnW+cxMABgG3bObodjZaavINYzuf1ww+4Myl8wi4bX64yB6z7je1tfsjtWljzAM3t88KurWEYfZ2dhxeJ3Vqbd2E3HHipgchd5roWYzSh3ZtzcdkXaKdctPt8s3zCXmxJnTa9/UmvH4mc7odc6ryUgrKee0g+I2KF0PE3SZqepFNb7XJvYtrUNEtNK+b5Q1JWLC0IaleLDJtio7DLOKmPhYKipZ7wcsfjVMbv9FxS1/b6iraYTuJN953uW0n5Kbblgi5t5KHZuk2kk/I7fAq5Hb78pU+zNfWrp21bfa4vZC7PSagMCHPnsPbW6DXspNhRce9RuprmWKWQwQKE/Eg59TO622lzuJEzNNjCE++S4FcRRspAEbGSbViTE8xpqVYI+JaWxja2H3Why3lUr4/lsxxJ7EuODpuOJfTfVP9cb2fOcJtd9va11rGMKeP3sZbyorfCLnd4kDGNm5C7ifq7dgnYIQ82ydXyItZ8tB8jsIWAypFdLyQVBXrl5BKKm/ol7Ci4qUW8TAj4f4XCAom4MWKdNey4Fb7+En54O9KLSDBFCW3FXFjW8Nt28/lMD6r87xvO30u2Qq8nbDbyLo18u0lXcV4PvN+b0Ju7Zsj5JnbAnb70z/DiJC7kU/Ibds6SLbXfPN8EXK7fXbi6lfIw8Bruoq5j+HLBKPjrnidvFkOiiXi+fATCfdKodHvQqkFuSalg2kqpHqxCLi26ceQGzE37TOex0JgD8j3B+JHwu3ae2mXEwmHc7qKMEivsZ9RyOEs5NYa5DkRcsN4jUKer9Z40Ai5dj9IhFzDTcjz9nGIkHtpZ2wbdFGgMKLjVpyE3OtkzqCR89xzVt/bddBUEz/9woqKl0PE/eaE5yOMCLhfKkG6Fam4GykBTFMh1UqOgGvYiTic95n225wnLBw/lxz3mwfgOTrudG67dBXtplPE3DqhE3a3nc+XI9fGy3gUcjuKnUPu1CdnfxGEXCNohRUvQm57PEB1FafruuEeOfch50WsnFLulULt8FVNpEpEPB9+01H8EJZ8F4NKk+VKG09NQhknVY9k2WAj3ZKlvdbGZr+n69nhJts2x5yi4Tn789zXP/OcxNumf7aPUbZhiphn+1iE3PZ29nxWybYTcmP+uLFtviorxmN+yh76jZA75YW7RdDtzm9u503IvVRYsfYB/OWPO1VX8ZI/XsrouNcvBYXmjVdj1B1wryseBuUS8WJEwoNEwIsR8R4KMltLj6VUME2FVC9Wmc6gS7ZFxD19rlvalPJ33U7O80e+cwfsmBduae89f9zbbaca5OY2uS+Ck5DbEVTI3fATIXe6bdsnT+UUr+3sI9qlyR/3gtfJnKY+AaPj1SrNlUCQqHgQwhBxb9fxNp5AFVdQPPEmxBZGxkk1I6TczYRDpNwYFc/x2RDfMJ3KFfre5zU6brxvlWqbaHjecofZM+r7/FRYMZ3TIOT6Loe2xpzvfNFzvwsDmSLdHiLkGqrItnfNMbeJZtvlhRerwkq+42FHx90oZ3S82BSrKotdikkYKSq+xxEgKu6e7uIuxd7qk/uPhHulkOh3KSPcCqSK2ggJA8p4FWMr3tbN2M5GvnPSWezw8H7jVbo9tbOkkQSJjpuOuaSr2LVzTVfRmvsWcmsbb0Kuty2CkPtJWTFS6ULuJV0le83w01W8Snexo+PVshpn0NU3/aSo+JXnQlfW9NvHq4R7uY7nFBcEE/BiSHc1y2+1jbdiYWScVC024m3FUcTt2ksOt8PCg5x7WRzI6b7ruYqQruI0odO2v7GNMZ3FcsytrGHYQu5GECF36mt3LMg+u+s74TUiXQq8ps4EGacf4S603jhTZILliYch4l7wmwfuh7DEm8JKnGDOOKkNLGJuSlnxKOp5z+9nOEH+SFxk2ppPnjvJM0903Lg/X7qKMKaWWCLrhnQVAJ4ndLotCOQm5G4lD4FwU1ZM44O9kLtFsr2WRwx7Qqef/HE9em46T2mj40a8i7rHdgWIjVW2y/0lJgi+SiL6jIqHLeL5ItheouFBIuFeKDTyXQnS7VRy0c9GSggj46RqyRMV19tk8CLfTsHggnGJaOe97ybalvueo+MeouF588edxmQVcrvz6fvshdxIOYXcOgbrbes+P+UR08eKO6HT7XixV6J0SknxLupDJ1WlmBR72fagIu5+zvzXLUYkPCz5LgblEGkKfOkoZWR82bJlmDBhAhKJBNrb27F+/XrHtqeeeiokScrZzjzzTL3NBRdckHP89NNP9zUm/u7UEFo03E6686arGN8/PbyXFuIx+VJT3NNLLO2CRsdt9udNVzFdzxIRzxx3yh83CrRd/rjxkn4XBbI7HpaQ56SnuES8ndqYjoWUP25t53Tc8wTKEpQ6dML7JE/7P7pCUlWKTaV+GQgrKl4OEfc00RPe5TBIBDxs+a4lya3WcZedEkXGH3zwQXR1dWHRokXYuHEjJk2ahM7OTrzzzju27R9++GHs2LFD31566SUoioJzzjnH1O700083tXvggQd8jYu/K1WMUb5zPqftRDwIIQQ5fH17dWnrS9K99PMTHXeU+dxxOuVk2+aYW4Xepr2fhX3CwOuiQFYhL2X+eND64+7ncY+uF0rY0fEwKPb/DgxFii3i+fA7GdMPYch3NYs2qX5uu+02zJs3D3PmzMExxxyD5cuXo76+HitWrLBtP2LECLS1tenbE088gfr6+hwZj8fjpnbDhw/3NS7+HdQaEpwj4nmi4vk+lz1/brt9XgQVc7f8GQ/RcWtOuCkP3LDfcTEg22u7RcQt/W0EP8iETtV0nuJEx8197SLU+SPepc4fz7bzlq5iV+rQiVLnjhcq6mHmjfulVJF3u0oqBZdEDCEqHlTE8+WH+4mG58NPJLzQCHg5xVuRpNA2UgIKiIz39vaatv7+fttLDAwMYMOGDejo6ND3ybKMjo4OrFu3ztMw7733Xpx77rloaGgw7V+7di1Gjx6NI488Epdccgneffddzw8doIzXBhJsJdzOMbwKdXGzLB0u4nTRAFHvQmZZ+6uu4hztzrltaBM0fzz3PMUTcjvyVVnxfCyk/HG3dkY8T5T0kK5SCIVWeQmSqlJIqoh1jNbnoRoj634W+SnVQkGF4icn3FO7AiLgxRbvcgg0xb34SAE3ABg/fjxaWlr0bcmSJbbX2L17N1KpFFpbW037W1tb0dPTk3eM69evx0svvYSLLrrItP/000/HT37yE3R3d+Omm27CU089hTPOOAOpVMrrw0fEc0tSeTj83TulrJj2O0XF7c5ZyPuLMPQ33rZeQpjHp9037dd25Du3QztJACJ7I3OhzHHjfpvzCAFIUmYfRLaf5fpCSJD0PpnbWrvM+bT9InNdYxv9mNYl8+C19pKUFjLZdB7J9bgKCTJcjjvcz/Yz79fQ21naG7G2MR3LtDceU4UMWVLt29vsSwkJitZXP1/6HNY+KUhQDMdNxwznyYe5n6xHZI2P37TfMB6vGPukILvWvbbtb/NaZI/Jvmp0l5tiT8BMX6O4z0e+qHiQfnp/L9f38R7uV8LDlO5akNlaeAxlwfg/1n76ANi+fTuam5v13fF4PLRhGbn33ntx/PHHY9q0aab95557rn77+OOPxwknnIDDDjsMa9euxac//WlP52ZkvAbQouA50XDDV0dbEc9HEd9TXL3H7phLFN1TdNxLyothv+d0FU3kc86l9bFEwi2RaNNtax/j6fJEwIPmlxdS8tCIW/64WzTdU755EdJVXB9LkUsdOkXHwy5z6JVqjG6XE79R8WKJuJfoc5B0FC8UGv1mVJlYKaSaSnNzs2lzkvGRI0dCURTs3LnTtH/nzp1oa2tzHV9fXx9Wr16NuXPn5n0shx56KEaOHIktW7Z4e/CgjFc1TqkoALyJtENU3PUjoNjvlx4+57xNyLQfqO/JnJ7P7SzX2f02p7UIuB1eFgTKh1O6ih+EkFzF2notvxKcPVaedBWnyZyVgJc0mbBSVbi4Txq/EXm/7QsV8bzj8fgn7icSHpaAlwoZUmgbKTIi4OaDWCyGKVOmoLu7W9+nqiq6u7sxffp0174PPfQQ+vv78S//8i95r/Pmm2/i3XffxZgxYzyPje+6tYgx8GsXLbe0qVQ0cZacRLiY0XFhiY5rNx2j4+bbOZFvm2tY5TpnASHjJWwi3aUud2i9bpD88ULKHVrPacTPYkD5jrldz2903OmaXqLjXnLK3cTZ7YuJny8bxV78p5AJp27YCW2pc72dc87999H75rumz2h4PoJEwYsV9S6nQFPga4Ouri7cfffduO+++/DKK6/gkksuQV9fH+bMmQMAmDVrFhYsWJDT795778WMGTNw0EEHmfbv27cP3/jGN/CHP/wBr7/+Orq7u3H22Wfj8MMPR2dnp+dxVZ2M+ynW/vLLL+OLX/wiJkyYAEmSsHTp0tINtJRIli1DjmPYvDd4DZAG+Z/snBU09ZPZ37Zt46WtE4VEx/O2tYt+2+xzio67pZF4ipanf3oVcq94FfJ8/V1F0IO82/bzGB3Pd7yU0fFi1dkuRunFoYqffHE/dcWLkevuRcS94EfC/RCWeNeK3FbjmMtKEaPiGjNnzsStt96KhQsXYvLkydi0aRPWrFmjT+rctm0bduzYYeqzefNmPPPMM7YpKoqi4M9//jM+//nP44gjjsDcuXMxZcoU/O///q+v3PWqmsCpFWtfvnw52tvbsXTpUnR2dmLz5s0YPXp0Tvv9+/fj0EMPxTnnnIOvf/3rZRhxkXH5O3cVcZcsC2Ft4/ZeUinvMwKuk0TNk0CNbaWsbdvtd2hrmsyZ01+C9izmTuBEZpJmbn9tX04fw6m9Tti0fYqEt8mcbqiQIAuYJ3lafua0R5422jGbCYd2Ey3dJn3aTea0u5b1/HbHch674bhxQqXTZM5C0CaaWsfoZSKn2xjcj3mfwGodh5++dv3DIIzn3XS+Ik4aDRoVD0PE/eSEe6VQ+aawEsCcA+6nTxDmz5+P+fPn2x5bu3Ztzr4jjzwSwiE/tK6uDr/97W+DDcRAVUXG/RZr/8QnPoFbbrkF5557btFm11YiXkW83ME1xz+kvBFq44NwauPj2kGj4z72OeWRG48VWu7Q6XzFLneo4Rbx9tLGy2TOfPvyTeb0QhjR8UImcvql0nK8rc9Zpa7CGRZhRcWLKeJ+U1LyUUj6Samj3dZa6YVupEj4jYoXEB2vRKrmXTKMYu1e6O/vzykgX+loeeG2+eFeRNzZE8MnyEX8pLXk3HZOA8k5ni8NxS533La/sZ3TbefxhbE6p9fqKkGEPN/KnKZjHgTb/Vi4kzmLnTtebPxeN2jeeKUJfiUQRr65U1S82CKeD78S7oew5bvc8kxxLw6FVFOpBarmHbfQYu1eWbJkial4/Pjx40M7d9jkyLeRQt4H8vQtpnv4iW776Z/3jzZAdNx5ZU5rO8l2v2mXRZ5tBTynbfAXIqxJc14maDq19TuuMKPj+coIBo2OF6PMof/Ifvhv68WexFlLBJm06USpRDwffiU8LPmudsGttvGWFUbGiZEFCxZgz549+rZ9+/ZyD8k/dn//eaLibufR2wd9X/EjwtYhOEbCHSLadvsc5M+uWorn6HjesVja2Uh2kMmctg81YHQ8SLqKXd1xKzmC7kGwqz06Xm68jsW9/KH3P/BKeuyFYJe77ifFJIzcckd5z9fP5eXyIq7e6pR7l/BCBbzaxZsUxlCPjFfNBM5CirX7IR6PV29+eQFR8tB+p4W36+Vt5/E8jpMzg7Q17rM5buqvNROZyZjGPl7Gnmlj6m85p5cJmU6TOUtBzoqdHq7tpa3r6pE2/fKtzJltZ14J0zhR0o4wVuUsBn4nQLo9n0OJQv4uvPb1GxUPKvP5RDwfXiTcD4VEwCtBuq3jV/n3QkpM1YQ3CinWXrNIls3peAbbnHK7PpWIx5QQwKlEYomi46ax5Y+OW/vbTuZ0iI67pasUKzpu7GOHYw1ylwh1OaLjbv38rMqZ97wVnqoSVt74UFzFs6gVV9yuW+BTHeaHfpBIeLGi32Eu6MMFgsoA01SqB7/F2gcGBrBp0yZs2rQJAwMDeOutt7Bp0yZfS5RWLdZoq/X9IU+WRzHeT+yCS1732R5zSevIuW3b3+YCeUTebmy2FY8s6Sq21VR0yc69TvbckqWt8wtTinQVq1z7qbjiZ2VON9EuVu64bU3yPLnjXhYBKiZB0maCjs3PF4JKW8G0nIQZFS80RzystBS/0hmmfFeDAFfquCqaIS7jVZOmAqSLte/atQsLFy5ET08PJk+enFOsXZazbzdvv/02Pvaxj+n3b731Vtx666045ZRTbGtJ1gwh/v1XQsDLLj3Ez3Ff6StufUxpLFKumWv7PKSp2NUTtx6zS2GxDsWarlIJeKlB7tTWrY3tMZsUDGs6il07P3XFgxx3unYhKSx+a44HwW/NcO/nDf6486UT1RqBl5svUMS9EiQSXsrrkeqklHXGK5GqknHAX7H2CRMmOBZqrynyiZ/XqLjbecJ8P8wnuk778uR0G/fpgp1XoH0ch+XchuO24mzp6y7gyF1IyKGfLuB2T4GA7WI++RYKyrcYkHmRHfNCPXZtrFj7+GlrPpa7EJC1n5FS5I47LQKUD7NcOy30E1xEjQJcjLzxQr4QFOvLRC1TSJ54GDnifiPhQaF8D1GCRLprSO+qTsaJAb8S7qGP5zaVgkeBNmIXCc8/udMlEp5vn27rlqYFRsdz23rvU2xCn9RZoui4XT/zOWTIHiWy2BM5ba8ZQHLdHlM5pLkcz5uG3RefwiZ+er8OUJy80WoR8ZIsAFSEN8bUUAj4lQBJCEg+n0u/7SsZJvbVKF5E3HNU3O68Ybynefk7smkTKPfc69+s3UROm/N4rjtudXXbiZ6W/j5yx90eVjFyx3MW/KnS3PGcNrbXce/nJ1/b70RO52sa2zs/D/mo9VUxi02p02aCRsULFXE/ec9+c8JDXwzIsCqo01YMSn09UpvwHXmoEOB9wdNnfJHeb6xybRuc8rrPdNxmwJ7PbdfX4z6Hc7rLuVmO3QhjIaCw8bMgkLWP72Mu0h4GfqqFhLEip19J93ROj+cJWhllKFZUMVKooAf5bS0kHcSLiHsdg99xFEO+K5FKHVdFIgJuNQJlvIYQUnYzkc8/JYfbPq/tfAF3fP0vsGdJtrlt2Oe5iovPfeWIjttRrMoqprYhRsddyxq6/FJ6jY5br+veP7wyh04rcpaKsCPgtbQSZ7ly1v1Ke9BShu7lEQsX8UIi4UGoBvkmwRnqi/7UzjvrEMQo336CUr5/f/2879XCe6TH6LnrG4GXb0RaU1epzh8dFzltK+9FCCs6nm3jt+62VcLtyhQW9rz5q8ttn6oSlEJXBi1GKcKhXt7QKV88LApJTynmta2EJeCkxmFknNQsEvLLcSW9xzn9YfnNEbeLgJsi1lLuOVxywe365t1nPYflvv3CP9oxm/NoTYoYHQ9Cvui4H/n2EtH2G0H3sghQvui4Rr6a4/nO6YViLADkeC2P53GT+yDiP5TxO3EzSFQ8jDxx1/4ePzQCLQhUpOi3XMJ/JBiMjJPaxHvasm0fP6nQrqcM8Y/FcxqJ2zGv4ykgOu4s024C7vyElyo67iVVxS5tJR9+VuV06mPfxvnty8tETq/RcT8TOe0od6qKkXJev1YmjVZCzfOgueJugptPnr2mpfheEChE+a4EQaacB4SRcVJzeBXxUn4uB/ijcRR5r3njfvGcH+7xm4pTxNzDFwg/0XE33HK/wyZI7njOOQqd5FnkX+p8kxSDpqoUkyDXGeqTMf1SiKD7jYoHTU8pRHb9RMN9nbdAAaf4klqBv721hEtaSt6PCqf3xJA/k23TRtxwaReoxKF+XptUFdf2/vblyLRtBN2DVHuQaGuZQy/CHXQip23bAqS70DKH1ja2/TzIqNdUFdPxEFNVvKSk5MtpLzRv3Ctuz4sfia+GFJegNcaLnS9ue02XY4WkpoQdDS8kEk75rl2YpkKqGwmuEg44uJ5b+wqPmAc+r99rlWsip2tU3OVaDucLIze8UPzkkFv7uLcpzUTOQp+7YlU1CVNoh/qEy2LjJ3oeJCoelDBE3CtBI+GU7yEA01RI1eLhfa0ov6tFcLq833D95oa7HQujvZ+0FJfreJmAWWiZQ6e+QaLj1mh4GGUONUo9kTModtFfX4+zgmZNhzGJ0w2vgl9Jz0k1EjQqXukiXkwJty46FMZGCmOoRsUBynhN43keYSW9h+T5A/MXhc7Tz2+006fYe0lV8XRMP5+H9BNL20osc6jhJf3E2ja0a3uMfgeR+ErMGzcStvhWQ7pJpTJUnrmgIh4mpRBninkBCBFsqxGGynvBkMOziFv7OeWcV9p7i89ccrd+7mUSi5iq4iGP3O69Rgj7tm6UciKnlUJW4vSUS15gVRWv5FsAyGsf57bhljI0SnK1TeKslcor1UqYUXG/Ih5WNLzcUWtKuT+YM06GNm7vFVXwPuJXvEPFT/1H1wi6l3wjD1HxAsoZ5ktVscNrqoofwoqA+4n+hrUAkJ8Sh0Gxk3SjdJdzXoAXgj4H1tfD+iWI8m5P0BSVfBRTxAuFElylMGec1BKuv59hvz+V6v3ORw6252M+I9+2CwB5aO+cquJ83UpckdN4vTCj6kFSVfxWVckeM7/deRX2fAsA2ZEvdaNSc6Q953iHVFGFeMcpWl2pH+J+RLzQaDjTQ0i1U6l/x8Qneb8k1sJ7lCcp95k37kKx2+clpDzybNv0z6FWVSXsa4WNaVGgMkR4GVUOl0LLGgZZddMvxY6K+xXxQsZCAa8NJDXYVivwXbiKKfh/aizvYdUS0PIivWHkjWfbe4yeu7V36lNgVRU/58k9r2T5md7vVlXFKTUlzKoqGqVcAKiQEodh5o07XT9sSQ8SmS/GRM2hMvmzHKt2Bi0j6L6oUGWIOCW8BmGaCql57N6z8r2PVcP7XDHSV/yey20SpdP5PYypElNVio2XBYCsbW2P+UhVCZtS541n99nnjRc6iZPkpxyiXakUW8SLKeGKJIeykWCUcgLnsmXLMGHCBCQSCbS3t2P9+vWObVeuXAlJkkxbIpEwtRFCYOHChRgzZgzq6urQ0dGBv/3tb77GxN+cWieE9y2TX1Sg0xVlRnWIEu8szV4mbno5f3WlqmgUK2XF1xgKKHFYq3njRrzmf7s9Vi4mZI+fZyXMfPGgtcXDlGC/Ih62hBdToinlASlRacMHH3wQXV1dWLRoETZu3IhJkyahs7MT77zzjmOf5uZm7NixQ9/eeOMN0/Gbb74Zd9xxB5YvX47nnnsODQ0N6OzsxIEDBzyPi78xtYyP966yulnm7ylbTrCw8/g/FjDP3C3dxXrOAlNMPL3nhJiqEiY5KSxhLwQUoMRh2JMMCxVORq2r48tJtRA0RcWJMNNT/Ih4GBJerug1pdwfpYqM33bbbZg3bx7mzJmDY445BsuXL0d9fT1WrFjhPDZJQltbm761trbqx4QQWLp0Kb7zne/g7LPPxgknnICf/OQnePvtt/HII494Hhd/U2oVt/evKvrMCxT1DijXQRcU8vVNplDZ9SPaBVzGrcSh37zxQq5fKIXmjdu3K07+u99xFAtGsmuXQlfcdOxbJBEvFMowMTIwMIANGzago6ND3yfLMjo6OrBu3TrHfvv27cMhhxyC8ePH4+yzz8bLL7+sH9u6dSt6enpM52xpaUF7e7vrOa3wN5XkUipZ92OLAcyyKDXIw8gzz3eekPLGc9tW7rewIAJfirzxIPXGw6RU1w8SmS/3l4ZapBSVVPyST4orTcSZv12lFDCBs7e317T19/fbXmL37t1IpVKmyDYAtLa2oqenx7bPkUceiRUrVuCXv/wlfvazn0FVVZx44ol48803AUDv5+ecdvC3tRbxGRWvYEfLJew0ijDPF+hc4Tz51ZY3HiRlxdq3UinWJM5KqKjiFbc0IK+Pn6k77jg9O2GnqIRBKUScAl7dFJKmMn78eLS0tOjbkiVLQhvX9OnTMWvWLEyePBmnnHIKHn74YYwaNQr/+Z//Gdo1ACAS6tlI+am89+GiIAkfXyIEgj0vbv1cjjmOTUjO4XoP1xICyPs5G+CxCiFBkoThZ/o6qpAghzQ7Vjt3Ief00lcVMmSH4rNO/VVIkC3fpOza2p07JSQoBTxHdtf2QgoSFAjTmFKQoSD3sRv3p4QMpYDivE7XIMWjkPQRI0FSVMKIinsV8UIknNQAQSZkZtpv374dzc3N+u54PG7bfOTIkVAUBTt37jTt37lzJ9ra2jxdMhqN4mMf+xi2bNkCAHq/nTt3YsyYMaZzTp482fND4W9xrSChKOUKixGEDLX6SVipLiFO4sx7Tpc+niLcHtroqSn5z+b5XCok51rkRcgb16j0SLgb5aqjXc3PmRNhp+wMlRrn5aSYIl6qSLjs4x8JTiGR8ebmZtPmJOOxWAxTpkxBd3e3vk9VVXR3d2P69OmexplKpfDiiy/q4j1x4kS0tbWZztnb24vnnnvO8zkBRsarmxA+m8L4zK6Ez31fkfKg53GNbAeIegeMvOtNMtFmL/hpWyn4iaK7tXWKPttFtd0i63nHABlyAVHjQqPWxcJr9D9o1JzR9vJSzKh4sSiWhBcq1Fp/lb/P/jHkgPvq45Ouri7Mnj0bU6dOxbRp07B06VL09fVhzpw5AIBZs2Zh3LhxeqrLddddh09+8pM4/PDD8cEHH+CWW27BG2+8gYsuughAutLKFVdcgeuvvx4f/ehHMXHiRFxzzTUYO3YsZsyY4XlclPGhgtf3Ta/55hUg4LZ4kNiCxT2ARLunmEhwelfxItGe0lcs5/PTp1RoQu0ndSNomocftJQQ03UDpNv4+3IR/EtBGBT6xaLY5yvWOUl4FDs9JSyKEdFO/2byd7MSmTlzJnbt2oWFCxeip6cHkydPxpo1a/QJmNu2bYMsZ38n3n//fcybNw89PT0YPnw4pkyZgmeffRbHHHOM3uab3/wm+vr6cPHFF+ODDz7ASSedhDVr1uQsDuSGJEQxKgvXDr29vWhpacGePXtMOUmVwPFXft97Y6eMCet+yeW4w+18bfTjdn2k8NrmthPe+hraOz0eYWlj7pd7zK69JFn6mx6n0G9LljaS9RyG80r6T0NbSz+9i0MfWRKux4z7ZQjnYw73NVnO2W8QU+s+rU/Ofrs+OW1Vz20Uy3HjMeO1FDif2xg1NsqhftxmX7qfmnNtxXT+7G2769vvc7iWYb/pGg7Xtj4WxeFc1i8MikOf3PPZ97N+scpt53xOa3u7lTGtX4asEfmc89t8ebK2ye2T08V+LLnNbKup+F3wxyli7SS9xYqMF0vGw4yKFzu1RIWKSNuWol4jCJXkNdpYTuy8DpGod3kFgOTgATz724UV8TgKhUlOQwGvIl6DFJyZETBf3PG6Hhb/KZgA5yll+cNiV1Lxs/iPvr+IK3GGSTHLCdZaecNyl6Uk4ROWiJcqx5t55D5QRbCtRuBvCiFV9Pfs5/+xChHrMMsfaqJbyHiqeUJivomCQSe5hiGblbbqZaWNh9gTVr54OVJUKMgVigi41Qj8rSQVTyjzDgs8RzVWVPFwKe/nMFzPb0WVQvC3cqVzWyfJK2VUu1gYo9NG8Tc+H8b9parfHUatceKPSqkxXgzhDSMqThGvXLQsTl9buQcdIvzNrHUq6be1mCkjYVBpYubh8fpb7KfCHp8LYa3EGejaBYhioUvJV4OkDqWSgNX8PzJGyj1JshKgiFc4Wp1xv1uNwN9OkqWS3q89RKIDRcwL/kLglvMd4AkM6b2kVOkrxaZS5ccuJaTYq4dWcj52qaiOLyeV+Ttb6fj5glBIVJw1wEk1wN/QWqZAb6xQLyo6lVaOu1jyHGZeeDnwNpnTYcJmFUheJVEL6TzEO4ykk1JTyKI/tQA/kUhFUUi020/fclRZcY9ee/l25KFJQZM2S/8BXEhllbAJP93F/Xx+or7liL4GSbkpdgoLo9DED4yIVxGcwElqEn5mDW0KeJMqZXnDYvUtdeS7VsobFpNC8+lLAWU/P2HUFy9VigqpHiQhAm21An/LCakkPEW/iz8M43VqrbxhoRVaiiHDlVLSr9pqjVeD4FcyTgv+1AKMilcZasCtRuBvK6k9/MhqEdNVAi38EzJhRrmFkAxSnd5XivKGfghyvVJEtYudvuFU3tAJo3RXyheBMPD7ZWIoVYapRoJGxSni1Qcj41XGsmXLMGHCBCQSCbS3t2P9+vWu7R966CEcddRRSCQSOP744/HYY4+VaKSkZijD33sl1hovF6FFxEsknZVU3tCp1nglE+brZE07qYS5CUOdYk8OpYhXKcwZrx4efPBBdHV1YdGiRdi4cSMmTZqEzs5OvPPOO7btn332WZx33nmYO3cuXnjhBcyYMQMzZszASy+9VOKRk7wU44+q0j94K1XyK4RKqzVejvKGlYpXsR+KC//U0ArdvglThJkrToYSVfXbftttt2HevHmYM2cOjjnmGCxfvhz19fVYsWKFbfvbb78dp59+Or7xjW/g6KOPxne/+118/OMfx5133lnikVcWNeAC4VPMD9AiP+HVJNiVRDVMhKyGMQ41OMmzcmFUvIrhoj/VwcDAADZs2ICOjg59nyzL6OjowLp162z7rFu3ztQeADo7Ox3bA0B/fz96e3tNG6luwqpF6nqeYr0nFOm8QphzvyuBIGUOC6k1XgrCvHa5JbBctcZrNYJOnAkSFaeIVzesM14l7N69G6lUCq2trab9ra2t6Onpse3T09Pjqz0ALFmyBC0tLfo2fvz4wgdfahi4qS7KsAonsafUkxkrXTQrsdY4KR/5yhrmg4sJEUeGeGQ8Uu4BVBoLFixAV1eXfr+3t7dihfzF73293EMghBBCCCkISU1vfvvUClUj4yNHjoSiKNi5c6dp/86dO9HW1mbbp62tzVd7AIjH44jH44UPmBBCCCGE5CdIpLuGIuNV8/+JsVgMU6ZMQXd3t75PVVV0d3dj+vTptn2mT59uag8ATzzxhGN7QgghhBBCSknVRMYBoKurC7Nnz8bUqVMxbdo0LF26FH19fZgzZw4AYNasWRg3bhyWLFkCALj88stxyimn4Hvf+x7OPPNMrF69Gs8//zzuuuuucj4MQgghhBCiEaRueO0ExqtLxmfOnIldu3Zh4cKF6OnpweTJk7FmzRp9kua2bdsgy9lg/4knnohVq1bhO9/5Dr71rW/hox/9KB555BEcd9xx5XoIhBBCCCHEQJAVNWtpBU5JiBp6NEWgt7cXLS0t2LNnD5qbm8s9HEIIIYSQwFSS12hjOW3KAkQiCV99k8kD+P2GJRXxOAqlqiLjhBBCCCGkxhAA/FZHqaFQMmWcEEIIIYSUjaGeplI11VQIIYQQQgipNRgZJ4QQQggh5UMgQJ3xooykLFDGCSGEEEJI+Rjii/5QxgkhhBBCSPlQAUgB+tQIzBknhBBCCCFlQ5vA6XcLwrJlyzBhwgQkEgm0t7dj/fr1jm3vvvtunHzyyRg+fDiGDx+Ojo6OnPYXXHABJEkybaeffrqvMVHGCSGEEEJI+dDSVPxuPnnwwQfR1dWFRYsWYePGjZg0aRI6Ozvxzjvv2LZfu3YtzjvvPPz+97/HunXrMH78eHz2s5/FW2+9ZWp3+umnY8eOHfr2wAMP+BoXZZwQQgghhJSPEsn4bbfdhnnz5mHOnDk45phjsHz5ctTX12PFihW27e+//35ceumlmDx5Mo466ijcc889UFUV3d3dpnbxeBxtbW36Nnz4cF/joowTQgghhJCqpLe317T19/fbthsYGMCGDRvQ0dGh75NlGR0dHVi3bp2na+3fvx+Dg4MYMWKEaf/atWsxevRoHHnkkbjkkkvw7rvv+noMlHFCCCGEEFI+CoiMjx8/Hi0tLfq2ZMkS20vs3r0bqVQKra2tpv2tra3o6enxNMx///d/x9ixY01Cf/rpp+MnP/kJuru7cdNNN+Gpp57CGWecgVQq5fnhs5oKIYQQQggpHwVUU9m+fTuam5v13fF4PLRhGbnxxhuxevVqrF27FolEQt9/7rnn6rePP/54nHDCCTjssMOwdu1afPrTn/Z0bkbGCSGEEEJI2Sikmkpzc7Npc5LxkSNHQlEU7Ny507R/586daGtrcx3frbfeihtvvBGPP/44TjjhBNe2hx56KEaOHIktW7Z4fvyUcUIIIYQQUj5KMIEzFothypQppsmX2mTM6dOnO/a7+eab8d3vfhdr1qzB1KlT817nzTffxLvvvosxY8Z4HhtlnBBCCCGElA9VBNt80tXVhbvvvhv33XcfXnnlFVxyySXo6+vDnDlzAACzZs3CggUL9PY33XQTrrnmGqxYsQITJkxAT08Penp6sG/fPgDAvn378I1vfAN/+MMf8Prrr6O7uxtnn302Dj/8cHR2dnoeF3PGCSGEEEJIzTNz5kzs2rULCxcuRE9PDyZPnow1a9bokzq3bdsGWc7GqX/0ox9hYGAAX/rSl0znWbRoERYvXgxFUfDnP/8Z9913Hz744AOMHTsWn/3sZ/Hd737XV+66JETAJYyGCL29vWhpacGePXtMEwQIIYQQQqqNSvIabSwdh16OiOJv4mUy1Y8n/357RTyOQmFknBBCCCGElJEgi/jUTiyZMk4IIYQQQspHkBU1ayixgzJOCCGEEELKhyrgO9IdYAJnpUIZJ4QQQggh5UOo6c1vnxqBpQ0JIYQQQggpE4yME0IIIYSQ8sGccUIIIYQQQsoEc8YJIYQQQggpE4yME0IIIYQQUiYEAsh4UUZSFijjhBBCCCGkfAzxyDirqRBCCCGEEFImGBknhBBCCCHlQ1UB+KwbrtZOnXHKOCGEEEIIKR9DPE2FMk4IIYQQQsoHZZwQQgghhJAywTrjhBBCCCGElAchVAjhLwfcb/tKhjJOCCGEEELKhxD+I901lKbC0oaEEEIIIYSUCUbGCSGEEEJI+RABcsZrKDJOGSeEEEIIIeVDVQHJZw44c8YJIYQQQggJgSEeGa+anPH33nsPX/3qV9Hc3Ixhw4Zh7ty52Ldvn2ufu+66C6eeeiqam5shSRI++OCD0gyWEEIIIYR4QqhqoK1WqBoZ/+pXv4qXX34ZTzzxBH7961/j6aefxsUXX+zaZ//+/Tj99NPxrW99q0SjJIQQQgghvtAW/fG71QhVkabyyiuvYM2aNfjjH/+IqVOnAgB+8IMf4HOf+xxuvfVWjB071rbfFVdcAQBYu3ZtiUZKCCGEEEKId6oiMr5u3ToMGzZMF3EA6OjogCzLeO6550K9Vn9/P3p7e00bIYQQQggpEqoIttUIVSHjPT09GD16tGlfJBLBiBEj0NPTE+q1lixZgpaWFn0bP358qOcnhBBCCCEGhEhXR/G1UcZD4eqrr4YkSa7bq6++WtIxLViwAHv27NG37du3l/T6hBBCCCFDCaGKQFsQli1bhgkTJiCRSKC9vR3r1693bf/QQw/hqKOOQiKRwPHHH4/HHnvMPHYhsHDhQowZMwZ1dXXo6OjA3/72N19jKquMX3nllXjllVdct0MPPRRtbW145513TH2TySTee+89tLW1hTqmeDyO5uZm00YIIYQQQoqE76i4GqjO+IMPPoiuri4sWrQIGzduxKRJk9DZ2ZnjmBrPPvsszjvvPMydOxcvvPACZsyYgRkzZuCll17S29x888244447sHz5cjz33HNoaGhAZ2cnDhw44HlckhCVH+d/5ZVXcMwxx+D555/HlClTAACPP/44Tj/9dLz55puOEzg11q5di9NOOw3vv/8+hg0b5uvavb29aGlpwZ49eyjmhBBCCKlqKslrtLGcKv0zIlLUV9+kGMRa8Qtfj6O9vR2f+MQncOeddwIAVFXF+PHj8bWvfQ1XX311TvuZM2eir68Pv/71r/V9n/zkJzF58mQsX74cQgiMHTsWV155Ja666ioAwJ49e9Da2oqVK1fi3HPP9TSuqsgZP/roo3H66adj3rx5WL9+Pf7v//4P8+fPx7nnnquL+FtvvYWjjjrK9N8NPT092LRpE7Zs2QIAePHFF7Fp0ya89957ZXkchBBCCCGk9AwMDGDDhg3o6OjQ98myjI6ODqxbt862z7p160ztAaCzs1Nvv3XrVvT09JjatLS0oL293fGcdlRFaUMAuP/++zF//nx8+tOfhizL+OIXv4g77rhDPz44OIjNmzdj//79+r7ly5fj2muv1e//wz/8AwDgxz/+MS644AJP19X+44BVVQghhBBS7Wg+U0mJEUnR7zvtJIlBALl+Fo/HEY/Hc9rv3r0bqVQKra2tpv2tra2O8xN7enps22vFQ7Sfbm28UDUyPmLECKxatcrx+IQJE3J+sRYvXozFixcXdN29e/cCAKuqEEIIIaRm2Lt3L1paWso6hlgshra2NjzT81j+xjY0Njbm+NmiRYsKdr9SUzUyXi7Gjh2L7du3o6mpCZIklXs4Or29vRg/fjy2b99e9pyvoQxfh8qAr0NlwNehMuDrUH4q+TUQQmDv3r1559uVgkQiga1bt2JgYCBQfyFEjpvZRcUBYOTIkVAUBTt37jTt37lzp2MxkLa2Ntf22s+dO3dizJgxpjaTJ0/2/Dgo43mQZRkf+chHyj0MR1jxpTLg61AZ8HWoDPg6VAZ8HcpPpb4G5Y6IG0kkEkgkEkW/TiwWw5QpU9Dd3Y0ZM2YASE/g7O7uxvz58237TJ8+Hd3d3fqK7gDwxBNPYPr06QCAiRMnoq2tDd3d3bp89/b24rnnnsMll1zieWyUcUIIIYQQUvN0dXVh9uzZmDp1KqZNm4alS5eir68Pc+bMAQDMmjUL48aNw5IlSwAAl19+OU455RR873vfw5lnnonVq1fj+eefx1133QUAkCQJV1xxBa6//np89KMfxcSJE3HNNddg7NixuvB7gTJOCCGEEEJqnpkzZ2LXrl1YuHAhenp6MHnyZKxZs0afgLlt2zbIcrbQ4IknnohVq1bhO9/5Dr71rW/hox/9KB555BEcd9xxeptvfvOb6Ovrw8UXX4wPPvgAJ510EtasWeMr2l8VdcZJLv39/ViyZAkWLFjgmB9Fig9fh8qAr0NlwNehMuDrUH74GhA/UMYJIYQQQggpE1Wx6A8hhBBCCCG1CGWcEEIIIYSQMkEZJ4QQQgghpExQxiuYZcuWYcKECUgkEmhvb8f69etd2z/00EM46qijkEgkcPzxx+Oxx4KtaEXM+Hkd7r77bpx88skYPnw4hg8fjo6OjryvG/GG378HjdWrV0OSJF9lpogzfl+HDz74AJdddhnGjBmDeDyOI444gu9NBeL3NVi6dCmOPPJI1NXVYfz48fj617+OAwcOlGi0tcnTTz+Ns846C2PHjoUkSXjkkUfy9lm7di0+/vGPIx6P4/DDD8fKlSuLPk5SJQhSkaxevVrEYjGxYsUK8fLLL4t58+aJYcOGiZ07d9q2/7//+z+hKIq4+eabxV/+8hfxne98R0SjUfHiiy+WeOS1hd/X4Stf+YpYtmyZeOGFF8Qrr7wiLrjgAtHS0iLefPPNEo+8tvD7Omhs3bpVjBs3Tpx88sni7LPPLs1gaxi/r0N/f7+YOnWq+NznPieeeeYZsXXrVrF27VqxadOmEo+8dvD7Gtx///0iHo+L+++/X2zdulX89re/FWPGjBFf//rXSzzy2uKxxx4T3/72t8XDDz8sAIhf/OIXru3//ve/i/r6etHV1SX+8pe/iB/84AdCURSxZs2a0gyYVDSU8Qpl2rRp4rLLLtPvp1IpMXbsWLFkyRLb9l/+8pfFmWeeadrX3t4u/vVf/7Wo46x1/L4OVpLJpGhqahL33XdfsYY4JAjyOiSTSXHiiSeKe+65R8yePZsyHgJ+X4cf/ehH4tBDDxUDAwOlGmLN4/c1uOyyy8Q//uM/mvZ1dXWJT33qU0Ud51DCi4x/85vfFMcee6xp38yZM0VnZ2cRR0aqBaapVCADAwPYsGEDOjo69H2yLKOjowPr1q2z7bNu3TpTewDo7Ox0bE/yE+R1sLJ//34MDg5ixIgRxRpmzRP0dbjuuuswevRozJ07txTDrHmCvA6/+tWvMH36dFx22WVobW3FcccdhxtuuAGpVKpUw64pgrwGJ554IjZs2KCnsvz973/HY489hs997nMlGTNJw89o4gZX4KxAdu/ejVQqpa8IpdHa2opXX33Vtk9PT49t+56enqKNs9YJ8jpY+fd//3eMHTs2502YeCfI6/DMM8/g3nvvxaZNm0owwqFBkNfh73//O373u9/hq1/9Kh577DFs2bIFl156KQYHB7Fo0aJSDLumCPIafOUrX8Hu3btx0kknQQiBZDKJf/u3f8O3vvWtUgyZZHD6jO7t7cWHH36Iurq6Mo2MVAKMjBNSJG688UasXr0av/jFL3wti0sKY+/evTj//PNx9913Y+TIkeUezpBGVVWMHj0ad911F6ZMmYKZM2fi29/+NpYvX17uoQ0Z1q5dixtuuAE//OEPsXHjRjz88MN49NFH8d3vfrfcQyOEZGBkvAIZOXIkFEXBzp07Tft37tyJtrY22z5tbW2+2pP8BHkdNG699VbceOONePLJJ3HCCScUc5g1j9/X4bXXXsPrr7+Os846S9+nqioAIBKJYPPmzTjssMOKO+gaJMjfw5gxYxCNRqEoir7v6KOPRk9PDwYGBhCLxYo65lojyGtwzTXX4Pzzz8dFF10EADj++OPR19eHiy++GN/+9rchy4zJlQKnz+jm5mZGxQkj45VILBbDlClT0N3dre9TVRXd3d2YPn26bZ/p06eb2gPAE0884die5CfI6wAAN998M7773e9izZo1mDp1aimGWtP4fR2OOuoovPjii9i0aZO+ff7zn8dpp52GTZs2Yfz48aUcfs0Q5O/hU5/6FLZs2aJ/GQKAv/71rxgzZgxFPABBXoP9+/fnCLf25UgIUbzBEhP8jCaulHsGKbFn9erVIh6Pi5UrV4q//OUv4uKLLxbDhg0TPT09Qgghzj//fHH11Vfr7f/v//5PRCIRceutt4pXXnlFLFq0iKUNQ8Dv63DjjTeKWCwm/uu//kvs2LFD3/bu3Vuuh1AT+H0drLCaSjj4fR22bdsmmpqaxPz588XmzZvFr3/9azF69Ghx/fXXl+shVD1+X4NFixaJpqYm8cADD4i///3v4vHHHxeHHXaY+PKXv1yuh1AT7N27V7zwwgvihRdeEADEbbfdJl544QXxxhtvCCGEuPrqq8X555+vt9dKG37jG98Qr7zyili2bBlLGxIdyngF84Mf/EAcfPDBIhaLiWnTpok//OEP+rFTTjlFzJ4929T+5z//uTjiiCNELBYTxx57rHj00UdLPOLaxM/rcMghhwgAOduiRYtKP/Aaw+/fgxHKeHj4fR2effZZ0d7eLuLxuDj00EPFf/zHf4hkMlniUdcWfl6DwcFBsXjxYnHYYYeJRCIhxo8fLy699FLx/vvvl37gNcTvf/972/d67bmfPXu2OOWUU3L6TJ48WcRiMXHooYeKH//4xyUfN6lMJCH4/1SEEEIIIYSUA+aME0IIIYQQUiYo44QQQgghhJQJyjghhBBCCCFlgjJOCCGEEEJImaCME0IIIYQQUiYo44QQQgghhJQJyjghhBBCCCFlgjJOCCGEEEJImaCME0IIIYQQUiYo44QQUgZSqRROPPFEfOELXzDt37NnD8aPH49vf/vbZRoZIYSQUiIJIUS5B0EIIUORv/71r5g8eTLuvvtufPWrXwUAzJo1C3/605/wxz/+EbFYrMwjJIQQUmwo44QQUkbuuOMOLF68GC+//DLWr1+Pc845B3/84x8xadKkcg+NEEJICaCME0JIGRFC4B//8R+hKApefPFFfO1rX8N3vvOdcg+LEEJIiaCME0JImXn11Vdx9NFH4/jjj8fGjRsRiUTKPSRCCCElghM4CSGkzKxYsQL19fXYunUr3nzzzXIPhxBCSAlhZJwQQsrIs88+i1NOOQWPP/44rr/+egDAk08+CUmSyjwyQgghpYCRcUIIKRP79+/HBRdcgEsuuQSnnXYa7r33Xqxfvx7Lly8v99AIIYSUCEbGCSGkTFx++eV47LHH8Kc//Qn19fUAgP/8z//EVVddhRdffBETJkwo7wAJIYQUHco4IYSUgaeeegqf/vSnsXbtWpx00kmmY52dnUgmk0xXIYSQIQBlnBBCCCGEkDLBnHFCCCGEEELKBGWcEEIIIYSQMkEZJ4QQQgghpExQxgkhhBBCCCkTlHFCCCGEEELKBGWcEEIIIYSQMkEZJ4QQQgghpExQxgkhhBBCCCkTlHFCCCGEEELKBGWcEEIIIYSQMkEZJ4QQQgghpExQxgkhhBBCCCkT/x+kdXcPY6YF1gAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 800x400 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from paraflow import run_simulation, SimulationParams, SimulationResult, get_frames, display_frame\n",
    "\n",
    "use_cache = False\n",
    "\n",
    "if not use_cache:\n",
    "    sim_results = run_simulation(\n",
    "        passage,\n",
    "        params=SimulationParams(\n",
    "            inlet_total_state=inlet_total_state, \n",
    "            target_outlet_static_state=target_outlet_static_state, \n",
    "        ),\n",
    "        working_directory=\"./simulation_out\", \n",
    "        id=\"1\",\n",
    "        auto_delete=True,\n",
    "        verbose=False,\n",
    "    )\n",
    "    sim_results.to_file(\"sim_results.pkl\")\n",
    "else:\n",
    "   sim_results = SimulationResult.from_file(\"sim_results.pkl\")\n",
    "\n",
    "frames = get_frames(\n",
    "    sim_results, \n",
    "    property_names=[\"Mach\"], \n",
    "    num_pnts=500,\n",
    ")\n",
    "display_frame(frames, \"Mach\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.11.7"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}