From 1e7f0601ed9b79ef795ec41ebb2ec19fd8aac5bb Mon Sep 17 00:00:00 2001
From: Jody Klymak <jklymak@gmail.com>
Date: Mon, 24 Sep 2018 10:56:04 -0700
Subject: [PATCH 1/3] ENH: add secondary x/y axis

---
 .flake8                                       |   1 +
 doc/api/axes_api.rst                          |   2 +
 doc/users/next_whats_new/2018-09-08-JMK.rst   |  16 +
 .../secondary_axis.py                         | 179 +++++++
 lib/matplotlib/_constrained_layout.py         |  11 +-
 lib/matplotlib/axes/_axes.py                  |  75 +++
 lib/matplotlib/axes/_base.py                  |  18 +-
 lib/matplotlib/axes/_secondary_axes.py        | 463 ++++++++++++++++++
 lib/matplotlib/axis.py                        |   6 +-
 lib/matplotlib/dates.py                       |   1 -
 lib/matplotlib/scale.py                       |  42 ++
 .../test_axes/secondary_xy.png                | Bin 0 -> 42332 bytes
 lib/matplotlib/tests/test_axes.py             |  33 ++
 lib/mpl_toolkits/mplot3d/axes3d.py            |  36 +-
 14 files changed, 861 insertions(+), 22 deletions(-)
 create mode 100644 doc/users/next_whats_new/2018-09-08-JMK.rst
 create mode 100644 examples/subplots_axes_and_figures/secondary_axis.py
 create mode 100644 lib/matplotlib/axes/_secondary_axes.py
 create mode 100644 lib/matplotlib/tests/baseline_images/test_axes/secondary_xy.png

diff --git a/.flake8 b/.flake8
index c7e9114821ff..007f636b27df 100644
--- a/.flake8
+++ b/.flake8
@@ -228,6 +228,7 @@ per-file-ignores =
     examples/subplots_axes_and_figures/axes_zoom_effect.py: E402
     examples/subplots_axes_and_figures/demo_constrained_layout.py: E402
     examples/subplots_axes_and_figures/demo_tight_layout.py: E402
+    examples/subplots_axes_and_figures/secondary_axis.py: E402
     examples/subplots_axes_and_figures/two_scales.py: E402
     examples/subplots_axes_and_figures/zoom_inset_axes.py: E402
     examples/tests/backend_driver_sgskip.py: E402, E501
diff --git a/doc/api/axes_api.rst b/doc/api/axes_api.rst
index 3e2f7cdd8c34..a3660fac3ba7 100644
--- a/doc/api/axes_api.rst
+++ b/doc/api/axes_api.rst
@@ -188,6 +188,8 @@ Text and Annotations
    Axes.inset_axes
    Axes.indicate_inset
    Axes.indicate_inset_zoom
+   Axes.secondary_xaxis
+   Axes.secondary_yaxis
 
 
 Fields
diff --git a/doc/users/next_whats_new/2018-09-08-JMK.rst b/doc/users/next_whats_new/2018-09-08-JMK.rst
new file mode 100644
index 000000000000..1c3bc954ae7f
--- /dev/null
+++ b/doc/users/next_whats_new/2018-09-08-JMK.rst
@@ -0,0 +1,16 @@
+:orphan:
+
+Secondary x/y Axis support
+--------------------------
+
+A new method provides the ability to add a second axis to an existing
+axes via `.Axes.secondary_xaxis` and `.Axes.secondary_yaxis`.  See
+:doc:`/gallery/subplots_axes_and_figures/secondary_axis` for examples.
+
+.. plot::
+
+    import matplotlib.pyplot as plt
+
+    fig, ax = plt.subplots(figsize=(5, 3))
+    ax.plot(range(360))
+    ax.secondary_xaxis('top', functions=(np.deg2rad, np.rad2deg))
diff --git a/examples/subplots_axes_and_figures/secondary_axis.py b/examples/subplots_axes_and_figures/secondary_axis.py
new file mode 100644
index 000000000000..466b1a2f8e32
--- /dev/null
+++ b/examples/subplots_axes_and_figures/secondary_axis.py
@@ -0,0 +1,179 @@
+"""
+==============
+Secondary Axis
+==============
+
+Sometimes we want as secondary axis on a plot, for instance to convert
+radians to degrees on the same plot.  We can do this by making a child
+axes with only one axis visible via `.Axes.axes.secondary_xaxis` and
+`.Axes.axes.secondary_yaxis`.
+
+If we want to label the top of the axes:
+"""
+
+import matplotlib.pyplot as plt
+import numpy as np
+import datetime
+import matplotlib.dates as mdates
+from matplotlib.transforms import Transform
+from matplotlib.ticker import (
+    AutoLocator, AutoMinorLocator)
+
+fig, ax = plt.subplots(constrained_layout=True)
+x = np.arange(0, 360, 1)
+y = np.sin(2 * x * np.pi / 180)
+ax.plot(x, y)
+ax.set_xlabel('angle [degrees]')
+ax.set_ylabel('signal')
+ax.set_title('Sine wave')
+secax = ax.secondary_xaxis('top')
+plt.show()
+
+###########################################################################
+# However, its often useful to label the secondary axis with something
+# other than the labels in the main axis. In that case we need to provide
+# both a forward and an inverse conversion function in a tuple
+# to the ``functions`` kwarg:
+
+fig, ax = plt.subplots(constrained_layout=True)
+x = np.arange(0, 360, 1)
+y = np.sin(2 * x * np.pi / 180)
+ax.plot(x, y)
+ax.set_xlabel('angle [degrees]')
+ax.set_ylabel('signal')
+ax.set_title('Sine wave')
+
+
+def deg2rad(x):
+    return x * np.pi / 180
+
+
+def rad2deg(x):
+    return x * 180 / np.pi
+
+secax = ax.secondary_xaxis('top', functions=(deg2rad, rad2deg))
+secax.set_xlabel('angle [rad]')
+plt.show()
+
+###########################################################################
+# Here is the case of converting from wavenumber to wavelength in a
+# log-log scale.
+#
+# .. note ::
+#
+#   In this case, the xscale of the parent is logarithmic, so the child is
+#   made logarithmic as well.
+
+fig, ax = plt.subplots(constrained_layout=True)
+x = np.arange(0.02, 1, 0.02)
+np.random.seed(19680801)
+y = np.random.randn(len(x)) ** 2
+ax.loglog(x, y)
+ax.set_xlabel('f [Hz]')
+ax.set_ylabel('PSD')
+ax.set_title('Random spectrum')
+
+
+def forward(x):
+    return 1 / x
+
+
+def inverse(x):
+    return 1 / x
+
+secax = ax.secondary_xaxis('top', functions=(forward, inverse))
+secax.set_xlabel('period [s]')
+plt.show()
+
+###########################################################################
+# Sometime we want to relate the axes in a transform that is ad-hoc from
+# the data, and is derived empirically.  In that case we can set the
+# forward and inverse transforms functions to be linear interpolations from the
+# one data set to the other.
+
+fig, ax = plt.subplots(constrained_layout=True)
+xdata = np.arange(1, 11, 0.4)
+ydata = np.random.randn(len(xdata))
+ax.plot(xdata, ydata, label='Plotted data')
+
+xold = np.arange(0, 11, 0.2)
+# fake data set relating x co-ordinate to another data-derived co-ordinate.
+xnew = np.sort(10 * np.exp(-xold / 4) + np.random.randn(len(xold)) / 3)
+
+ax.plot(xold[3:], xnew[3:], label='Transform data')
+ax.set_xlabel('X [m]')
+ax.legend()
+
+
+def forward(x):
+    return np.interp(x, xold, xnew)
+
+
+def inverse(x):
+    return np.interp(x, xnew, xold)
+
+secax = ax.secondary_xaxis('top', functions=(forward, inverse))
+secax.xaxis.set_minor_locator(AutoMinorLocator())
+secax.set_xlabel('$X_{other}$')
+
+plt.show()
+
+###########################################################################
+# A final example translates np.datetime64 to yearday on the x axis and
+# from Celsius to Farenheit on the y axis:
+
+
+dates = [datetime.datetime(2018, 1, 1) + datetime.timedelta(hours=k * 6)
+            for k in range(240)]
+temperature = np.random.randn(len(dates))
+fig, ax = plt.subplots(constrained_layout=True)
+
+ax.plot(dates, temperature)
+ax.set_ylabel(r'$T\ [^oC]$')
+plt.xticks(rotation=70)
+
+
+def date2yday(x):
+    """
+    x is in matplotlib datenums, so they are floats.
+    """
+    y = x - mdates.date2num(datetime.datetime(2018, 1, 1))
+    return y
+
+
+def yday2date(x):
+    """
+    return a matplotlib datenum (x is days since start of year)
+    """
+    y = x + mdates.date2num(datetime.datetime(2018, 1, 1))
+    return y
+
+secaxx = ax.secondary_xaxis('top', functions=(date2yday, yday2date))
+secaxx.set_xlabel('yday [2018]')
+
+
+def CtoF(x):
+    return x * 1.8 + 32
+
+
+def FtoC(x):
+    return (x - 32) / 1.8
+
+secaxy = ax.secondary_yaxis('right', functions=(CtoF, FtoC))
+secaxy.set_ylabel(r'$T\ [^oF]$')
+
+plt.show()
+
+#############################################################################
+#
+# ------------
+#
+# References
+# """"""""""
+#
+# The use of the following functions and methods is shown in this example:
+
+import matplotlib
+
+matplotlib.axes.Axes.secondary_xaxis
+matplotlib.axes.Axes.secondary_yaxis
diff --git a/lib/matplotlib/_constrained_layout.py b/lib/matplotlib/_constrained_layout.py
index a72b79ca7ad1..39795b544a3c 100644
--- a/lib/matplotlib/_constrained_layout.py
+++ b/lib/matplotlib/_constrained_layout.py
@@ -180,7 +180,8 @@ def do_constrained_layout(fig, renderer, h_pad, w_pad,
             sup = fig._suptitle
             bbox = invTransFig(sup.get_window_extent(renderer=renderer))
             height = bbox.y1 - bbox.y0
-            sup._layoutbox.edit_height(height+h_pad)
+            if np.isfinite(height):
+                sup._layoutbox.edit_height(height+h_pad)
 
         # OK, the above lines up ax._poslayoutbox with ax._layoutbox
         # now we need to
@@ -266,10 +267,14 @@ def _make_layout_margins(ax, renderer, h_pad, w_pad):
     """
     fig = ax.figure
     invTransFig = fig.transFigure.inverted().transform_bbox
-
     pos = ax.get_position(original=True)
     tightbbox = ax.get_tightbbox(renderer=renderer)
     bbox = invTransFig(tightbbox)
+    # this can go wrong:
+    if not (np.isfinite(bbox.width) and np.isfinite(bbox.height)):
+        # just abort, this is likely a bad set of co-ordinates that
+        # is transitory...
+        return
     # use stored h_pad if it exists
     h_padt = ax._poslayoutbox.h_pad
     if h_padt is None:
@@ -287,6 +292,8 @@ def _make_layout_margins(ax, renderer, h_pad, w_pad):
     _log.debug('left %f', (-bbox.x0 + pos.x0 + w_pad))
     _log.debug('right %f', (bbox.x1 - pos.x1 + w_pad))
     _log.debug('bottom %f', (-bbox.y0 + pos.y0 + h_padt))
+    _log.debug('bbox.y0 %f', bbox.y0)
+    _log.debug('pos.y0 %f', pos.y0)
     # Sometimes its possible for the solver to collapse
     # rather than expand axes, so they all have zero height
     # or width.  This stops that...  It *should* have been
diff --git a/lib/matplotlib/axes/_axes.py b/lib/matplotlib/axes/_axes.py
index 57f4cc6ca991..13e8cda4e622 100644
--- a/lib/matplotlib/axes/_axes.py
+++ b/lib/matplotlib/axes/_axes.py
@@ -34,6 +34,7 @@
 import matplotlib.tri as mtri
 from matplotlib.container import BarContainer, ErrorbarContainer, StemContainer
 from matplotlib.axes._base import _AxesBase, _process_plot_format
+from matplotlib.axes._secondary_axes import SecondaryAxis
 
 _log = logging.getLogger(__name__)
 
@@ -599,6 +600,80 @@ def indicate_inset_zoom(self, inset_ax, **kwargs):
 
         return rectpatch, connects
 
+    @docstring.dedent_interpd
+    def secondary_xaxis(self, location, *, functions=None, **kwargs):
+        """
+        Add a second x-axis to this axes.
+
+        For example if we want to have a second scale for the data plotted on
+        the xaxis.
+
+        %(_secax_docstring)s
+
+        Examples
+        --------
+
+        Add a secondary axes that shows both wavelength for the main
+        axes that shows wavenumber.
+
+        .. plot::
+
+            fig, ax = plt.subplots()
+            ax.loglog(range(1, 360, 5), range(1, 360, 5))
+            ax.set_xlabel('frequency [Hz]')
+
+
+            def invert(x):
+                return 1 / x
+
+            secax = ax.secondary_xaxis('top', functions=(invert, invert))
+            secax.set_xlabel('Period [s]')
+            plt.show()
+
+
+        """
+        if (location in ['top', 'bottom'] or isinstance(location, Number)):
+            secondary_ax = SecondaryAxis(self, 'x', location, functions,
+                                         **kwargs)
+            self.add_child_axes(secondary_ax)
+            return secondary_ax
+        else:
+            raise ValueError('secondary_xaxis location must be either '
+                             'a float or "top"/"bottom"')
+
+    def secondary_yaxis(self, location, *, functions=None, **kwargs):
+        """
+        Add a second y-axis to this axes.
+
+        For example if we want to have a second scale for the data plotted on
+        the yaxis.
+
+        %(_secax_docstring)s
+
+        Examples
+        --------
+
+        Add a secondary axes that converts from radians to degrees
+
+        .. plot::
+
+            fig, ax = plt.subplots()
+            ax.plot(range(1, 360, 5), range(1, 360, 5))
+            ax.set_ylabel('degrees')
+            secax = ax.secondary_yaxis('right', functions=(np.deg2rad,
+                                                           np.rad2deg))
+            secax.set_ylabel('radians')
+
+        """
+        if location in ['left', 'right'] or isinstance(location, Number):
+            secondary_ax = SecondaryAxis(self, 'y', location,
+                                         functions, **kwargs)
+            self.add_child_axes(secondary_ax)
+            return secondary_ax
+        else:
+            raise ValueError('secondary_yaxis location must be either '
+                             'a float or "left"/"right"')
+
     def text(self, x, y, s, fontdict=None, withdash=False, **kwargs):
         """
         Add text to the axes.
diff --git a/lib/matplotlib/axes/_base.py b/lib/matplotlib/axes/_base.py
index 2e8359b4288b..868a459538e5 100644
--- a/lib/matplotlib/axes/_base.py
+++ b/lib/matplotlib/axes/_base.py
@@ -2500,8 +2500,17 @@ def _update_title_position(self, renderer):
             title.set_position((x, 1.0))
             # need to check all our twins too...
             axs = self._twinned_axes.get_siblings(self)
-
-            top = 0  # the top of all the axes twinned with this axes...
+            # and all the children
+            for ax in self.child_axes:
+                if ax is not None:
+                    locator = ax.get_axes_locator()
+                    if locator:
+                        pos = locator(self, renderer)
+                        ax.apply_aspect(pos)
+                    else:
+                        ax.apply_aspect()
+                    axs = axs + [ax]
+            top = 0
             for ax in axs:
                 try:
                     if (ax.xaxis.get_label_position() == 'top'
@@ -2544,6 +2553,8 @@ def draw(self, renderer=None, inframe=False):
 
         # prevent triggering call backs during the draw process
         self._stale = True
+
+        # loop over self and child axes...
         locator = self.get_axes_locator()
         if locator:
             pos = locator(self, renderer)
@@ -4315,6 +4326,9 @@ def get_tightbbox(self, renderer, call_axes_locator=True,
         if bb_yaxis:
             bb.append(bb_yaxis)
 
+        self._update_title_position(renderer)
+        bb.append(self.get_window_extent(renderer))
+
         self._update_title_position(renderer)
         if self.title.get_visible():
             bb.append(self.title.get_window_extent(renderer))
diff --git a/lib/matplotlib/axes/_secondary_axes.py b/lib/matplotlib/axes/_secondary_axes.py
new file mode 100644
index 000000000000..83960601c260
--- /dev/null
+++ b/lib/matplotlib/axes/_secondary_axes.py
@@ -0,0 +1,463 @@
+import collections
+import numpy as np
+import numbers
+
+import warnings
+
+import matplotlib.docstring as docstring
+import matplotlib.ticker as mticker
+import matplotlib.transforms as mtransforms
+import matplotlib.scale as mscale
+import matplotlib.cbook as cbook
+
+from matplotlib.axes._base import _AxesBase
+
+from matplotlib.ticker import (
+    AutoLocator,
+    AutoMinorLocator,
+    FixedLocator,
+    FuncFormatter,
+    LogFormatterSciNotation,
+    LogLocator,
+    NullLocator,
+    NullFormatter,
+    ScalarFormatter
+)
+
+from matplotlib.scale import Log10Transform
+
+
+def _make_secondary_locator(rect, parent):
+    """
+    Helper function to locate the secondary axes.
+
+    A locator gets used in `Axes.set_aspect` to override the default
+    locations...  It is a function that takes an axes object and
+    a renderer and tells `set_aspect` where it is to be placed.
+
+    This locator make the transform be in axes-relative co-coordinates
+    because that is how we specify the "location" of the secondary axes.
+
+    Here *rect* is a rectangle [l, b, w, h] that specifies the
+    location for the axes in the transform given by *trans* on the
+    *parent*.
+    """
+    _rect = mtransforms.Bbox.from_bounds(*rect)
+    bb = mtransforms.TransformedBbox(_rect, parent.transAxes)
+    tr = parent.figure.transFigure.inverted()
+    bb = mtransforms.TransformedBbox(bb, tr)
+
+    def secondary_locator(ax, renderer):
+        return bb
+
+    return secondary_locator
+
+
+class SecondaryAxis(_AxesBase):
+    """
+    General class to hold a Secondary_X/Yaxis.
+    """
+
+    def __init__(self, parent, orientation,
+                  location, functions, **kwargs):
+        """
+        See `.secondary_xaxis` and `.secondary_yaxis` for the doc string.
+        While there is no need for this to be private, it should really be
+        called by those higher level functions.
+        """
+
+        self._functions = functions
+        self._parent = parent
+        self._orientation = orientation
+        self._ticks_set = False
+
+        if self._orientation == 'x':
+            super().__init__(self._parent.figure, [0, 1., 1, 0.0001], **kwargs)
+            self._axis = self.xaxis
+            self._locstrings = ['top', 'bottom']
+            self._otherstrings = ['left', 'right']
+        elif self._orientation == 'y':
+            super().__init__(self._parent.figure, [0, 1., 0.0001, 1], **kwargs)
+            self._axis = self.yaxis
+            self._locstrings = ['right', 'left']
+            self._otherstrings = ['top', 'bottom']
+        # this gets positioned w/o constrained_layout so exclude:
+        self._layoutbox = None
+        self._poslayoutbox = None
+
+        self.set_location(location)
+        self.set_functions(functions)
+
+        # styling:
+        if self._orientation == 'x':
+            otheraxis = self.yaxis
+        else:
+            otheraxis = self.xaxis
+
+        otheraxis.set_major_locator(mticker.NullLocator())
+        otheraxis.set_ticks_position('none')
+
+        for st in self._otherstrings:
+            self.spines[st].set_visible(False)
+        for st in self._locstrings:
+            self.spines[st].set_visible(True)
+
+        if self._pos < 0.5:
+            # flip the location strings...
+            self._locstrings = self._locstrings[::-1]
+        self.set_alignment(self._locstrings[0])
+
+    def set_alignment(self, align):
+        """
+        Set if axes spine and labels are drawn at top or bottom (or left/right)
+        of the axes.
+
+        Parameters
+        ----------
+        align :: string
+            either 'top' or 'bottom' for orientation='x' or
+            'left' or 'right' for orientation='y' axis
+
+        """
+        if align in self._locstrings:
+            if align == self._locstrings[1]:
+                # need to change the orientation.
+                self._locstrings = self._locstrings[::-1]
+            elif align != self._locstrings[0]:
+                raise ValueError('"{}" is not a valid axis orientation, '
+                                 'not changing the orientation;'
+                                 'choose "{}" or "{}""'.format(align,
+                                 self._locstrings[0], self._locstrings[1]))
+            self.spines[self._locstrings[0]].set_visible(True)
+            self.spines[self._locstrings[1]].set_visible(False)
+            self._axis.set_ticks_position(align)
+            self._axis.set_label_position(align)
+
+    def set_location(self, location):
+        """
+        Set the vertical or horizontal location of the axes in
+        parent-normalized co-ordinates.
+
+        Parameters
+        ----------
+        location : string or scalar
+            The position to put the secondary axis.  Strings can be 'top' or
+            'bottom' for orientation='x' and 'right' or 'left' for
+            orientation='y', scalar can be a float indicating the relative
+            position on the parent axes to put the new axes, 0.0 being the
+            bottom (or left) and 1.0 being the top (or right).
+        """
+
+        # This puts the rectangle into figure-relative coordinates.
+        if isinstance(location, str):
+            if location in ['top', 'right']:
+                self._pos = 1.
+            elif location in ['bottom', 'left']:
+                self._pos = 0.
+            else:
+                raise ValueError("location must be '{}', '{}', or a "
+                                 "float, not '{}'".format(location,
+                                 self._locstrings[0], self._locstrings[1]))
+        else:
+            self._pos = location
+        self._loc = location
+
+        if self._orientation == 'x':
+            bounds = [0, self._pos, 1., 1e-10]
+        else:
+            bounds = [self._pos, 0, 1e-10, 1]
+
+        secondary_locator = _make_secondary_locator(bounds, self._parent)
+
+        # this locator lets the axes move in the parent axes coordinates.
+        # so it never needs to know where the parent is explicitly in
+        # figure co-ordinates.
+        # it gets called in `ax.apply_aspect() (of all places)
+        self.set_axes_locator(secondary_locator)
+
+    def apply_aspect(self, position=None):
+        self._set_lims()
+        super().apply_aspect(position)
+
+    def set_ticks(self, ticks, minor=False):
+        """
+        Set the x ticks with list of *ticks*
+
+        Parameters
+        ----------
+        ticks : list
+            List of x-axis tick locations.
+
+        minor : bool, optional
+            If ``False`` sets major ticks, if ``True`` sets minor ticks.
+            Default is ``False``.
+        """
+        ret = self._axis.set_ticks(ticks, minor=minor)
+        self.stale = True
+        self._ticks_set = True
+        return ret
+
+    def set_functions(self, functions):
+        """
+        Set how the secondary axis converts limits from the parent axes.
+
+        Parameters
+        ----------
+        functions : 2-tuple of func, or `Transform` with an inverse.
+            Transform between the parent axis values and the secondary axis
+            values.
+
+            If supplied as a 2-tuple of functions, the first function is
+            the forward transform function and the second is the inverse
+            transform.
+
+            If a transform is supplied, then the transform must have an
+            inverse.
+
+        """
+
+        if self._orientation == 'x':
+            set_scale = self.set_xscale
+            parent_scale = self._parent.get_xscale()
+        else:
+            set_scale = self.set_yscale
+            parent_scale = self._parent.get_yscale()
+        # we need to use a modified scale so the scale can receive the
+        # transform.  Only types supported are linear and log10 for now.
+        # Probably possible to add other transforms as a todo...
+        if parent_scale == 'log':
+            defscale = 'functionlog'
+        else:
+            defscale = 'function'
+
+        if (isinstance(functions, tuple) and len(functions) == 2 and
+            callable(functions[0]) and callable(functions[1])):
+            # make an arbitrary convert from a two-tuple of functions
+            # forward and inverse.
+            self._functions = functions
+        elif functions is None:
+            self._functions = (lambda x: x, lambda x: x)
+        else:
+            raise ValueError('functions argument of secondary axes '
+                             'must be a two-tuple of callable functions '
+                             'with the first function being the transform '
+                             'and the second being the inverse')
+        set_scale(defscale, functions=self._functions)
+
+    def draw(self, renderer=None, inframe=False):
+        """
+        Draw the secondary axes.
+
+        Consults the parent axes for its limits and converts them
+        using the converter specified by
+        `~.axes._secondary_axes.set_functions` (or *functions*
+        parameter when axes initialized.)
+
+        """
+
+        self._set_lims()
+        # this sets the scale in case the parent has set its scale.
+        self._set_scale()
+        super().draw(renderer=renderer, inframe=inframe)
+
+    def _set_scale(self):
+        """
+        Check if parent has set its scale
+        """
+
+        if self._orientation == 'x':
+            pscale = self._parent.xaxis.get_scale()
+            set_scale = self.set_xscale
+        if self._orientation == 'y':
+            pscale = self._parent.yaxis.get_scale()
+            set_scale = self.set_yscale
+        if pscale == 'log':
+            defscale = 'functionlog'
+        else:
+            defscale = 'function'
+
+        if self._ticks_set:
+            ticks = self._axis.get_ticklocs()
+
+        set_scale(defscale, functions=self._functions)
+
+        # OK, set_scale sets the locators, but if we've called
+        # axsecond.set_ticks, we want to keep those.
+        if self._ticks_set:
+            self._axis.set_major_locator(FixedLocator(ticks))
+
+    def _set_lims(self):
+        """
+        Set the limits based on parent limits and the convert method
+        between the parent and this secondary axes
+        """
+        if self._orientation == 'x':
+            lims = self._parent.get_xlim()
+            set_lim = self.set_xlim
+            trans = self.xaxis.get_transform()
+        if self._orientation == 'y':
+            lims = self._parent.get_ylim()
+            set_lim = self.set_ylim
+            trans = self.yaxis.get_transform()
+        order = lims[0] < lims[1]
+        lims = self._functions[0](np.array(lims))
+        neworder = lims[0] < lims[1]
+        if neworder != order:
+            # flip because the transform will take care of the flipping..
+            lims = lims[::-1]
+        set_lim(lims)
+
+    def get_tightbbox(self, renderer, call_axes_locator=True):
+        """
+        Return the tight bounding box of the axes.
+        The dimension of the Bbox in canvas coordinate.
+
+        If *call_axes_locator* is *False*, it does not call the
+        _axes_locator attribute, which is necessary to get the correct
+        bounding box. ``call_axes_locator==False`` can be used if the
+        caller is only intereted in the relative size of the tightbbox
+        compared to the axes bbox.
+        """
+
+        bb = []
+
+        if not self.get_visible():
+            return None
+
+        self._set_lims()
+        locator = self.get_axes_locator()
+        if locator and call_axes_locator:
+            pos = locator(self, renderer)
+            self.apply_aspect(pos)
+        else:
+            self.apply_aspect()
+
+        if self._orientation == 'x':
+            bb_axis = self.xaxis.get_tightbbox(renderer)
+        else:
+            bb_axis = self.yaxis.get_tightbbox(renderer)
+        if bb_axis:
+            bb.append(bb_axis)
+
+        bb.append(self.get_window_extent(renderer))
+        _bbox = mtransforms.Bbox.union(
+            [b for b in bb if b.width != 0 or b.height != 0])
+
+        return _bbox
+
+    def set_aspect(self, *args, **kwargs):
+        """
+        Secondary axes cannot set the aspect ratio, so calling this just
+        sets a warning.
+        """
+        cbook._warn_external("Secondary axes can't set the aspect ratio")
+
+    def set_xlabel(self, xlabel, fontdict=None, labelpad=None, **kwargs):
+        """
+        Set the label for the x-axis.
+
+        Parameters
+        ----------
+        xlabel : str
+            The label text.
+
+        labelpad : scalar, optional, default: None
+            Spacing in points between the label and the x-axis.
+
+        Other Parameters
+        ----------------
+        **kwargs : `.Text` properties
+            `.Text` properties control the appearance of the label.
+
+        See also
+        --------
+        text : for information on how override and the optional args work
+        """
+        if labelpad is not None:
+            self.xaxis.labelpad = labelpad
+        return self.xaxis.set_label_text(xlabel, fontdict, **kwargs)
+
+    def set_ylabel(self, ylabel, fontdict=None, labelpad=None, **kwargs):
+        """
+        Set the label for the x-axis.
+
+        Parameters
+        ----------
+        ylabel : str
+            The label text.
+
+        labelpad : scalar, optional, default: None
+            Spacing in points between the label and the x-axis.
+
+        Other Parameters
+        ----------------
+        **kwargs : `.Text` properties
+            `.Text` properties control the appearance of the label.
+
+        See also
+        --------
+        text : for information on how override and the optional args work
+        """
+        if labelpad is not None:
+            self.yaxis.labelpad = labelpad
+        return self.yaxis.set_label_text(ylabel, fontdict, **kwargs)
+
+    def set_color(self, color):
+        """
+        Change the color of the secondary axes and all decorators
+        Parameters
+        ----------
+        color : Matplotlib color
+        """
+
+        if self._orientation == 'x':
+            self.tick_params(axis='x', colors=color)
+            self.spines['bottom'].set_color(color)
+            self.spines['top'].set_color(color)
+            self.xaxis.label.set_color(color)
+        else:
+            self.tick_params(axis='y', colors=color)
+            self.spines['left'].set_color(color)
+            self.spines['right'].set_color(color)
+            self.yaxis.label.set_color(color)
+
+
+_secax_docstring = '''
+Warnings
+--------
+
+This method is experimental as of 3.1, and the API may change.
+
+Parameters
+----------
+location : string or scalar
+    The position to put the secondary axis.  Strings can be 'top' or
+    'bottom', for x-oriented axises or 'left' or 'right' for y-oriented axises
+    or a scalar can be a float indicating the relative position
+    on the axes to put the new axes (0 being the bottom (left), and 1.0 being
+    the top (right).)
+
+functions : 2-tuple of func, or Transform with an inverse
+
+    If a 2-tuple of functions, the user specifies the transform
+    function and its inverse.  i.e.
+    `functions=(lamda x: 2 / x, lambda x: 2 / x)` would be an
+    reciprocal transform with a factor of 2.
+
+    The user can also directly supply a subclass of
+    `.transforms.Transform` so long as it has an inverse.
+
+    See :doc:`/gallery/subplots_axes_and_figures/secondary_axis`
+    for examples of making these conversions.
+
+
+Other Parameters
+----------------
+**kwargs : `~matplotlib.axes.Axes` properties.
+    Other miscellaneous axes parameters.
+
+Returns
+-------
+ax : axes._secondary_axes.SecondaryAxis
+'''
+docstring.interpd.update(_secax_docstring=_secax_docstring)
diff --git a/lib/matplotlib/axis.py b/lib/matplotlib/axis.py
index 127c9c5b0fb3..e1e1ff61cd02 100644
--- a/lib/matplotlib/axis.py
+++ b/lib/matplotlib/axis.py
@@ -1159,11 +1159,11 @@ def get_tightbbox(self, renderer):
             if (np.isfinite(bbox.width) and np.isfinite(bbox.height) and
                     a.get_visible()):
                 bb.append(bbox)
-
         bb.extend(ticklabelBoxes)
         bb.extend(ticklabelBoxes2)
-
-        bb = [b for b in bb if b.width != 0 or b.height != 0]
+        bb = [b for b in bb if ((b.width != 0 or b.height != 0) and
+                                np.isfinite(b.width) and
+                                np.isfinite(b.height))]
         if bb:
             _bbox = mtransforms.Bbox.union(bb)
             return _bbox
diff --git a/lib/matplotlib/dates.py b/lib/matplotlib/dates.py
index dbfcae354baa..915791db71c0 100644
--- a/lib/matplotlib/dates.py
+++ b/lib/matplotlib/dates.py
@@ -309,7 +309,6 @@ def _from_ordinalf(x, tz=None):
 
     # add hours, minutes, seconds, microseconds
     dt += datetime.timedelta(microseconds=remainder_musec)
-
     return dt.astimezone(tz)
 
 
diff --git a/lib/matplotlib/scale.py b/lib/matplotlib/scale.py
index ddd1bcbcd3d2..fb55d33156a0 100644
--- a/lib/matplotlib/scale.py
+++ b/lib/matplotlib/scale.py
@@ -409,6 +409,47 @@ def limit_range_for_scale(self, vmin, vmax, minpos):
                 minpos if vmax <= 0 else vmax)
 
 
+class FuncScaleLog(LogScale):
+    """
+    Provide an arbitrary scale with user-supplied function for the axis and
+    then put on a logarithmic axes
+    """
+
+    name = 'functionlog'
+
+    def __init__(self, axis, functions, base=10):
+        """
+        Parameters
+        ----------
+
+        axis: the axis for the scale
+
+        functions : (callable, callable)
+            two-tuple of the forward and inverse functions for the scale.
+            The forward function must have an inverse and, for best behavior,
+            be monotonic.
+
+            Both functions must have the signature::
+
+                def forward(values: array-like) -> array-like
+
+        base : float
+            logarithmic base of the scale (default = 10)
+
+        """
+        forward, inverse = functions
+        self.base = base
+        self.subs = None
+        transform = FuncTransform(forward, inverse) + LogTransform(base)
+        self._transform = transform
+
+    def get_transform(self):
+        """
+        The transform for arbitrary scaling
+        """
+        return self._transform
+
+
 class SymmetricalLogTransform(Transform):
     input_dims = 1
     output_dims = 1
@@ -660,6 +701,7 @@ def limit_range_for_scale(self, vmin, vmax, minpos):
     'symlog': SymmetricalLogScale,
     'logit':  LogitScale,
     'function': FuncScale,
+    'functionlog': FuncScaleLog,
     }
 
 
diff --git a/lib/matplotlib/tests/baseline_images/test_axes/secondary_xy.png b/lib/matplotlib/tests/baseline_images/test_axes/secondary_xy.png
new file mode 100644
index 0000000000000000000000000000000000000000..2e91f63326eedb185fb2f5637c3801a4b0739dfb
GIT binary patch
literal 42332
zcma%j2Rzs9`u|r(%1#I&Ba*#05eXs6-g_lvZ?dy8GRi0;WRu9=BT}-n$%w2FW%R#3
zJWtO#zyCSE|KIC5&p9V|_q^`=x~}*8eSN}iE6L)Wr8tX1q44BxN~@w!=u0RR+BP;8
z{7KQ3(`@ix7>-hMcd_Be3)?gl{{I>Kn_7-26oC=)AGBPF94q+4TTU{XPIv6gom`C`
zn4xTqob2z}Io-1|zT{%|z|qR?J`X!5J2$(%oud;MFDD<11s6A$sgW^1k09$MOD89L
zVGfRe{tmm{0}BqXyhcwH>Jmy$TJo-2%C`wm-Mh|55*x#~iIQKg3S1{VU+jh^)1+gg
z<;l#l|Lcw}Nxw9A1&Khzv$r)*w2LO{?vY-8N@|(>Y)bpxCr$ZJK_9Q7QE|A<yz~q8
z)5S{qd2nQyiYJ3<p5)lO;UP}Fb=t8v***P(S1_5~LXXf@+CUW1$&Z{on=_j9$xka*
zYUuILP12F)#3#S+LQy)=o&1m+1}L30WI(_7G3=zNX`}i(CV>mMSJ>Hu%gcH0UVRdu
zO-Mx4^^ow?;NalK&Q6ea;rQ=No;kL~_LLb-SMjKddHVvk5EEM=+b%ab%I(=%<GuZT
z1+7y=L@}%+=g&8O_@KPZEgV6vEG<1yEn-%ig;xz7`z@T()78yl4GImtN&4#5tLgQ1
zyPlib#_c-2kwPt{YhD*^x@R<X@H;*~_#C?#%fBRm|H9tEVaD_l9o;X?Y#Q^<h)U6a
z#)j5!@3gb*cH927cJ)9!;gc1|;Fd@-ovb((jeqH{-*Y$q*&M9@|IP9He{^c-SxIvC
z?AhC@st<h*e>J^-U)WX_9E`=l#FUwl0Xss3efI3~Ybn_4&f28S)a%&(0|R%zEH0WX
zeSCq9hgY}4uTiYu_H%1-W^N8$LP7%Pu9de^il9F+lTu@EFTRnHQDp@EB5#a4*LC3z
zTM0?!TekwEqX|5A*CP8eW#{(S8|+GJU%x*8;>8QMUte>G$jCybrwv=%+MHeLVBw~8
z=3a|dmX?j{KD#2fQiVxqeD-MDBPl~hMn=Yp(&OR+o<G0j=H_NNkgGCMVibUxUqPLD
zBeHp8re$-u;fSO=kwf<OZNf>nIlQ?Cy;gKvM!`-+WBf~KXlT$Jy{xQETkmMx@q=V$
z-?jX~+3BW$(CBE{d-qtsZ*7SbYQiCKQ}NVq-TVCdjKuE)o(B&eJWff82slk(ytVLt
zZ>?rsTW@M~l>CFlZ&7X@9xkV)+c+9FVhbPAqd$F8wXD!i5%VN1uK#say+9M$64Z1G
zBR7-z&Bbcl0CA&cjGh#s^1&7&3QEdSix1SST4fj$)z0WBX&ITOi3zG2*Y&kpYh>$9
zd-i(E9T$~)(!?mmJbz$(Jm^|$4I_?u{P^SYFIY=xSeT@x<&}FdajjQsTD5Lm9@|Ss
z2kw0tGK5l>L`0rBKJDyDqD{@e&J@<&-F;nAQ89e|Zr<Iifq{Y3Bc&$hEKy+9y@7;#
zOH*XG$38MRGn=<WtP2J-`Hhc{>(o0^Q&7AsEia)=|At4w8w$NWEcklJsMbEU!{AWJ
z%)xOQRaH|{Sv55=F;P{(C!C5eoZ7qRsHUbSo0bMWKaj$%{XQsn?{HJJzPCQG71P8p
zQ(aCkp(|;VN<&XCGB`NcC0ftZp<un>%^Ol2GWOPZ<~xMv&ON*yf{$kTI|WN{XrN+T
zw1JNTkC#HRhB`BooK8U1+$%&}QX~JF0I}5TDg$gI-Jz$>D??3i>)p0?ch@`n;RTVC
z$M)A~=rdl4UPMFPk@D1%am44GQ|#+2wc&``S>@y9HC!IZD=95)ilP!hCH<(jFUVur
z+}*X~ni?+D8D8rs8CY-XWpHkHwT$=j@<O35U%t$?#=*?Y#cuQd-4hlSEiJ2Bcr={8
ze1(Hzpbu=+*;k%Br0GRN>8EjH@U(a;=sx+$IVaNF)paH+Dr%Lw=EtmNWl>(<nV6WE
z&JQ0tDyMvKG&D3&RO$Vwj;GZL4<9}x-SVD)_hi)v_s*o<gG(V1>$d}_FI>P@R#wi+
z&L&^)|LWA)(<ANXwr5+}Z^;XfjFkgD9i3`yWS7N;(mO3wqe3*6N$<!tD;#bfbFU6K
z=LibES8XF`*R-`iQ*{*82@X|$TG>hBvmt%p=opidg3H0dF;^DVYvpI6sMr18v92zp
z68&Y})R$}c<i2BAbuV9HX%~{_D5qng&XbTx<_~5R6yRD{Q1tcnjZ`_%!)2ijZ0zP}
zz^%}Z?CiW5P6AJ-sob+Wy~ZuB7ATajyBlZP#VOg@N5aA6Q|`vb#v5Hs>6hr~(XVZ3
zP4MvWlo~djf?0Fa1<rRTQVx?z!DW<!e{%X~mIw1cT31xLZDug+undmBbu9J$?NhtI
zR+C#$z#C;HBZCRAOR|<ta3|wjo9(x;N@Y~l{2Q;I4`>4kot>SH;dGmOdzGr*zCB}S
zXXh~8=szZ=n}3ys1rr5lH(Im`UWAoxZNHq(P$m{zGlqdbC*|ww#Zqa0?m3a+@AKV3
zsi~<O8ykdhc*flcZ0%-uUfv>tyUcr}tEV@yYy)p(-@@X`I}`7K$VkQfxhkJS586QY
zh1Y}^Z4_oTb(B0iKA$5uwtDc+?FpfPcR)kxqwsJQ6&3%Z!~K@Y$$Gspbu~3GT9PME
zo`{;lY)o0icP{GxTxLK`j~Wm(emecj^|NCaGi5T|6@C3Ea97)3`%fz=DW(3>)$=YG
zYl#=a1BVs{4lR{$=w<yBtIU1!sTbV^hq1-Q9A4huiryFCA=CQou=kLwDp6^9IIoFK
zpl2^?Yh%MNBqU_O?4zP!P8h4~nRNjn^IGm+XE;3u|JLnrx3As0sK$)_hZ`LR1dpo_
zm&ASgt+^=TkTcZ1p2kEse0z_EY8@Xt{o&B%(%q{oRYO8VL@FALixC!xt9FVLlf~IQ
zSR6(9F?F1t{XjsFplA=bQl+CyY;F|s@Q2N0S~7e6X(oiHl6V+V9OCC#^QRf1`zRHu
z_Hx7x{)*0;YulU4(SjD_Yte(zcrkj%ee5$|KIsE(k8$ywu|(HHNmPH!LCKFzqucrO
zWaHpI2E^sxxgj2OV4g})<)MX(xs)AYl5hr@{-JM+8Rc;k7y0g!CfV7R3Z;UZSDxaY
ziVohJjr>Ym660_@(%2QQ<B?Pnoj07wNf|H7Wa>}P7>SjNs~$Be!Z@xu11AgC`}Cq*
zgK(>lh?kdq<C!7Ny@NySk<Xl@@HWINb&(fGT!sJDb7fUmD8tqC+OTOOdO+=FZ$V7S
zotwx=eXGo3b5HItLnK0Os$2Ew=NQlJDa*!IG<`P#Yld&O<RLJayagV4aRK5NZf!XF
zU@qA1pQ*s}xV5wxw!0W8X<|d~taa=;ZMsg@b@?JoRq$IKUn#`^YLk5r!*1ftX~Ue_
z<6TZDU_M@1%wRY3tJajluen2fw4x=a-}F)j-T(0p_dXo?FAkfN9WotMI<G-%Z6o|@
z`57VUbd;@fRB&cp;eZ)jftav@lRb(V6n1_~p`z38Q22(3_UQh|(HBcWt9UA8uTAU7
zKOjL&!Kcwr);&?mErU<7a98FZMK6T3j5aAPmx1+nN)M%MBEenfatcdLRI<(|hj|b5
zhwBba9ep$QtL9Y5h7UZ^@;tGQ{0F7ASmkfBG7GdrB*ju=N-mwcAN4~uQ3i%%Tt55P
zcsH5`wR;XO;U(rT_f3tRIfFy@F|)=QQvk8m0X6P4BuLT<%CK=xdV3kK?a9v4${<x3
z7tw~(r<;1Du)5@27PDikSFP=SvXfMC883*$8GGBlwJ3>e!^P!sH7jy8R5VWZDa$fj
zZ3pRztlMU0DwEQsJ;Faad<W}zxIPk|<KdXqwJPv=Dzjz|o=|K*SP9nSGMQc(Qv7QV
z^!QGWh3<*Lu#t-w#rJpC_R;z1E^QJ1={1>u?RfQrlRaq5%e89!w&#?>{uZBx3=bMf
z+C>);Av-=)kZ*A^@%VS2Z%ErJhWvQISDV5&Le7KbpEi-r9u8Z6tM$Zq&#J0Tg)+U#
zrmw6zOFgk#>L`ivur?*@?>2Vr$l0BLQC&~uSa}#|67<C0^>s^c64%oD9P3{;XMW&Y
zid@>!{64;!E&`s%$-a0+6W<=gv7mdpyhm)VjIn1g0SRDZ3I_RkjgEfcUVN*GW^r1F
zL0Z-@v^f8-#U*W48Y53gnc7rt$3rUF$390_Q~2KO2@0QTIjD9)Zi^PWx!1~ZPppFi
z%S3SNF7~p>x_VOHOQ(st?JJPijl5XITyu6G)R#}eLhrw{`DCKftdf=X_ThaO{=JPZ
z`#{Vp1>}gkbR{Lx<mKf*78^)Zy{)*nM?pb>gnEx2orP_@bmhuY!fO7jR~r+qQ!h$O
ziFT)aC2H&H*wj1{n<5%`d>$e*{q)haF}C$qR?5|YfK-aHg~%Ij>D0%j$eoasgiS+3
zv%NM!nJF9H4u5^h*xkyCS?c=r&Fy|=dUkgFrlzLQ@Nk>b-cpl}j}5-QwR^MSZd>mK
zh|Zr6L_dwMY0ptW!*jHDyxx1f%XpgAEW}f(tp4rlk$X%xOzu5~haB!<4<-hNVM{Px
z=eu{-&OaSLh@FnSfnKN{4lyg8xHxrbX(<F90!74BRFMXU%Nm`%y*F;$z#`?)k8s~w
zkj&}VvKA<y;HkGg*;n4Lc$z<~Nb9@ma6Kvc#NCX~Gj8w3Bj&E^cZQhJaG@u;kowK!
z$Ve!JrG{OxmzRbLDN&cMUOl_Ewl-4dDS(;=XC>&deI<dz0RM?RCQ7s3O9&ibOl++2
zTdGpkp1n6F$T7$EQ(`3ey5Jg6HDHKPlFTkdw8E=*ym{wf->n^1I&((d-Q9EZ@*vX9
z-4uT!AN2S!;YAS_T$KI6o&)OI&6}YPRN&(Ge*K~tei3EVUP*xXso40YOB>^9b9;38
z7v5Nt@PD8PlMC=VHfqzxsVSPGqN4V(YUh>DhZpAO=TRuv-2^5&M^{$?b`}nfkbr=I
zP2wLzLoa@qv{$BgSvVBu=eIQcKD5)V_aa9U4L>yW9u+boB4;$LfW~(p$~U9aaf~w6
zXV!USrv<Ki$m@|Fp5{r5UsAHztan+Zz^4(9c=hTkF12V?TTV8EGA1S_gl}9&1$~w<
z0c+<ii&&{p(p?<di0%RF9ky;e{SfVHof{B8wr61-e(8#H+wPO?8yLXw_4PG1HJxrF
zGiXi!ePH1@IBd?Ol+qYR%w%X`5mHjZx%cy@RxOsJ)xCSt>guGu=@JB}3&KuU9tH$l
zE_x@ryuAGQ*)tpn`!FytZt5x3kdTpOR#jECc~eAUf%U}Y(rxAC<~D+%M*sTdfqPMq
zkp4ylr_X5xX2^kV8yNU@A5v3OFZJg{V~`|=DCJzYw!SK_pg=@O=;tuk&b(Ifb|%q?
znv^s+T8VeMpC0pez8@Zs)aJ($&n))Ms9SYj(XpM(b77z9xVVTui0`w)-$g@t?yj}X
zMN8b*O@o_f)&H_-Y>aYu+>zi@fB$ADee`s{vV<8rqX|1ZI~&Dz)W)XOQR|QP^>s?P
zw-F8b-d<kM`@i1S(mJD|F-*ELK9Y8{a-SRmQyhaYSFX;PT~9J>4m>j}(C{<T($bPp
z3L64nR)7mSJ3B4=e?atg$`6StvvG+?Nf{xPYHjpKm$`lWaY<z*S#iTrwOllFwvgij
zHWL#Qm)D*>muV;A3-$cA0<H2u{o=%G-Z(bB>bodHh>2AbKhbj3ZO7h*KpEE5u48Xs
zGGLZ`Z?QMExvNWV``hRpNW8FRWo2XH;~AW_R{6d5Hb{#Vch@E`QACuKEY4O`F2A;|
z`f6P5`tmh6xU4=g%&e`k--x6z`@XKooFeETg@K6~=I4jv@@dAH`8n2h9KCvp%Xu9+
z`VY)cs7OgY-1zZ>6or8*bzPq%?N*SNZ-wkuW__}*t!{Uc<=g6No9BAn@QX!m4gRO{
z?BvR>9i{gZo<6<mvNq0)<i+0U_~hKuJxN^GFW_Q9FrA#7ES0S~Q8t+h!DddTiJqRG
zclS#1@JdilP7d-=MeYs`595%LkwFsA6iAqnn|l*QIOVk!>5$4}(Tk0m?V$1Hb+G$V
zNN?hu6Tx`_ced?!xB9Ji9aWw7>J~f?TiZg|u-Wnt;x-Pxuw$ClPI;qe?}fjACJ_Gp
z`-8a;b`6JbZs*;l6A&2A>Ui=l2!|{pPA_z!Q6=y0?A8Z~#FDH1RKiX-F|e>6!3B^J
zSb``$z8bgu*U!cDpu)n!ky`ii&HATLDI9i|)z0K|M{jRCp<O$G_ZJS#h{wDKv2jA-
zR3+@bacmDax6nibAFbJUpVx0BtvCsnE~LLjc6sS=JaKmC1HQc8XT6D4g5Pf)7ax86
zcw0qBhb543HZxlM+>akWW@p>NSAE!pcu^4cSBpeBD-`lsN@33|FJGN>Ureo;bY%hn
z$ZB2FputBkMbusOqL>HR&U)=ul%V_ftSfA6BvRO=#%*V3*4H`M^=q%-Vg>AV8ywL=
zI<9?C$K(Fpl*{FtCS)+oPQs~4{I*Xo@ZU>P^SbH6v`V!iV%AENz|LpxHN7Qg?0Un;
zBECS=mV5=D=0pF-kEnnql&9xD>F(Ir*xX{e@7hv}OD9aXFSEm!!@j|?j}DvmacZh<
zpK#CROP5eCE-vlqOt+FyuCBYJ!&Y}GQk6oMmMm+xdPVFXJczBn_$x?Y(hYsCJ+cKl
ziQE+u7M^@SM{@b|8(Yd?T=Iv2@+j!Zhmeqxf-i2~s~Lhj9e{-=rJ@qk`lLcxPmc<D
z6ESg&H`E_czd0H?AfbN0`9USAT6T7Mt#_b%_ex|6cUq?|9X)+aQc`GH88<fV#!Lu}
zRPnIM_b*+li8st2KYjYBDd03cD=RK)qlf2%;l^~6+iJxy3YD*3*u2!2*#}O2d#Mjy
z=gqw)i*)ay+4a=DJvT7xy}iB0mKOBU(NV-#HJoyY8#4<`trJsyrx!uqu;v?V^y`W*
zv;5uH_#gNoN8er<CV>0QH|e@O!2D%lVe-eV5aHCd1cl2^!371ZsK)1#0mAc}Y}t_A
zlaP=gmjd4Ezz%B9Tg$P+_jZ;CFNk|psRd4tmRoF&m_!@3N04ha_;`n~@6s|cGc#AM
z*Y9sE=NHey`r7X(iv>kQ;O2F6Yr!Uo|NKJSlP+Pfta|GlIeA!C7DI}VBNobjV_MSG
z%*^Q5_SbV7L*fF>2iL-|)_>dMqc*<JUGF=5kT6SG|3<f7*6K^jrZm+~ugH31cQ?-D
z<fOg5eIw29Jz1}x^NHBFxV5YNDZ}qJYo-qIGkz{jlU52$c@d+gr>FgV59Vm9>}O8-
z9Bj&-KYxC@=`_u~%?se}rM0ywz!}(o=_Eodc=J%&^Yu=wOQBhydcKUm4fd<Uc9doB
z-c-o$o0}6u?OAcHwj(8I{Vx@-L#PE44)8tNN?#hQv<p}E^|b=vOS8_Sf?3&JqoSfh
zv&x?K<*nq5B_2m}SYB%=5ej8tVexIUPN$UV?K?m}yu50Sqb>I?g9{Tq*l1arYLH+z
zsE>dvOD^KN);iM?Vkfdb-($%6JSz)hy>6Gfxuu26dHL=|VJB?FSpCsY+fRA#%IO)0
zcpv>-%q=WzgB{IcZ(~qMWS8;r5y!&4fP&DePCvS;O7PZGUb*E@pW7f!_omi+co_8l
zH<|v9C1>-*O3broW3z!`@&gHAt1%aO(mKy009M|eoQ)5@N5rsEtvTsFkHDSHzypXp
z!a5VdTc3YHoYFizyc@N0<&YtrYJPp{d;eOZzQ-A`S!?HsXFRhjUL+|UyuNNF%YBAD
z<e9Z`lHEaac)uVBh)GS21Pjoualte2J>*j_)OH?gxu#yE8-`dng1Zf}5jUO3!^2Aw
zcD@R(ti+@Pzg(>7JaQo)K0I~z?%hk^ckkpV$;iu}g?pLBUa0r<=~GVMBd>#1JKyq0
zJcjIAWm#<>xVX5m&zzAoHl|y6|BOjdh4GwJV}C!ffzQv&1J-<`teSo-8byZhpYlDr
z|FP0`EVwI<sinJH0j%`<kdB>EM-<iG-dr@%#fuRD8@IN%!-~lW#I_#~%vX#q#Fuko
z2n?^NY@p-Y{+1B+O$=6M{`}x4_2pQdcH-pOgP&U?_1+>s)~BKke*YRq`Ke^yn1Ss>
z5U1JMS?!6at*6zoKXgB>C6YUA$P0|P^BPEYS2EyjwX^&V)nZQeqS^*b+i0}R4E0n(
z=BAt+7q$EYp2!>hFKo!($lSR@jJzr7>Cu2Gzi68$iicPj-^a&iV{5A=zqtNUO9;Un
zg!^i`0Y^K-2G{ck=|n`v1eJT!#KMt9L$_{+JHvKAj~>P1Qt)D555`sbW>f&E7~HN4
zA};sBC91*h^a~6c^qb$^a-3<#%*x7|4(F(UTzx(co3_My@G8U(SO5i~{9@>3XZJSe
zG;3T6jVZQweLGtyr_xTLSTC%e$w{lPa(Ta5gqN=Rdj7kZ{N=<5CT=xbTQ=1A19xcw
z0fD))A$$Eh8lNT=KJUN-HE~1qQTnbG8j3H()p+fjfwQ{1Ua%aVhl>mUMJg(y3l~Cb
zYeiZ0YcJG!>@)+O_Fm24+0u?j&r?2}w{PDHx@|n4mEf|8tgAbuLYFQ&I6UNXT^lEz
zGUB>``*CF_x45`ro|q}SZ)K;->nC@$>pHtLGZiT{bri%5q`OZN5}Nz_?_@l>Rq4I&
z0?+G}m;QBBb@~LUO2casSLPW2+8R1!s!IIrfjZi_=8Jm&+p<=tCFi`qzyHSeHu9vj
zN7E2I|LK_z@kWx6<IC1T=M-J;wA?~#YwO#$ZyO#QYTmtzQK*hXKtO<dffH3<Uk|&+
zcJAD{=`UXbAw(dBht%!&;qK+4l$-8l<>i~s`eKs`0QYTQ5ylAll$j*rn&0|_IbDuM
z=U{*8cRNJzNXQKDK?XeA`@8E%(EoPXCJr3C@$aKUx7`U>@Uv%NNJyg~9}FKdVt*!J
z*ADT=bFf}5tuV!B{4-#piFC^61dLm;fxW1x!h^ZZbwtbR>ZZ&du^b-luhUCN&|n3i
zC2|<FHwWRUXlo;Hf}Eaq2%Gk6?EN3B*B?CKMqOfJI*r)ltdf*V0O2Wk9XsGzGPJSz
zE-4jTR#DLgR#Ix!&!}7NbmrRDO5Eeee!!s^fUmsYF@xhapC5<v%O=?F@Sfc1?CuWG
zFYfr-pTlL+j@!oT#mC17P^zH!FJ8ajM>7WM>R15Rf-k=}UPbit=TAbZ>G?#1VC}*i
zrlwh~Ph3y7ea(r~PDDh+>&K`CL?#U56B9EF3&FB9-sgSyT4>Jilkm|oA09rNYhO6{
z{36w!S+FO9lk86MybHuw4v?6s=<1HI*}ivXlcS`hM4L_+`?a&8=%JezK_j7`d-p0t
z(3q&1uV0xbefDj^*FA#_=IRgOCE$Zj`9ZLE$`5f)^TPS^+n;TU&qUs!2h=r4KKmRF
zDJ!SGtBzIu{+bBdHN*ietM<@wa41Jm3Ml0cT0`7H@Ip1O71GZ(xE9&$t5N4Zv4{GR
z$RR<To9&dG5$A<g*gK8^>x#`o9*5ahd_5-Jvv3eLd=`+cvO@@4RHOu?liJAW=oz?Z
zGpnm?_!k8~zF3UycTSjpyy@h6apeuq1>B61lG{j3J}@-&#?$kJ<IAZtIqZ)3Jr#}r
zC?g|7spzA}YqJ71shh454r*#7u!v*yd0oT_6C#BxSFb+I&u8)2S<c>Dx73b<C<jR}
z)6&wa7SkYO>bUnRKi^?<Rt93*@7g=uch%K{fnWmKvioRvK8O6L=4SLlKELQre!b?;
zkC$LyQSTeZ&f?${t{pmFhXnKL)vJIh{Kii=Xwfd~)BbV&jNnj0OX-`7Fxl4@1?%y!
z7F%)~Gn(I4ZE(bPZ%qlh0umDw)zTIlf!_ie=1DbfHC}d4%?Y6UzUlDGUZc*#vCqDh
z74RwnMk4Q28NswQ-BAWOAE|*U$&9jRlqB|aSSgy@W{e>Me#MX&w(H;{NjbOHr%mI7
z&3;yFv>x2`NVRBuqhmk^&Ais^7(a^higP?unr_;3$*)tzM{|GgcJAOY7@2L-u0xG<
zv0Q@AB)`$qdus{r)E(LF=IXo_as=^G_ao=PJBEB|nLn$tDv3T~3}`f=t#>o=z`l0@
zY4NBdBJp?@)|CM{-6fneS15K)#f3}+RMH>!Fa$lY<y(Jr7N!EvikYg!1Z-H-VT^#y
zXBwN&WPbboqxh+H50j#!p1<hwo5J603Z4pS9(uXzKv#Ll^`v5z0N?<<m*%3$N|yr(
zRmI$O*az6>=j0;y>FtJIOANi{K<(!6WT&!p)Drkgd1s;#;Qwf^DQ_{f8>^pA=s!wA
zFE4N&Dx;BaJUJgow-(zqnP?h9+iDyYA)P3H(_<4QNJLa2*B+Fg)F6hw6Nxo0@Db;k
zMYQS3HYyeyvi?+d?2Skyl?BXhZF%wFY|q#aVr$Cns>%kY*^K7rr4=zJsAXSi##N!K
zmqB9imli<TAgKBVI*&fqy=<@Vl(NGbYHG>zoA*Vd-cl=DkJGr9zdS5`>5Wz;v5;^_
zH!{b?Ec2ygbRSxg+twzUWp2LYb-fO1<>f;z0mt<ZksUD`l$HS4${--LO!}8=g$}Jt
zC=D-8i}M%1jie~ZeOiA`MoQW|HD&NBhlwtureycb!`=1a=JuvDJ3mLSOq8F^h+(>3
z9tdFzKpJv-H$pb|GygVS`vMUIC6&wOwsP0iyor~$%3~_5Xgxft74#C+Y`SX`wuzk}
zYDB&F4cv&7m1f|`kuNf7>b-XtK<<GLnCoo5XD8-Jq46L<Ob;6x8o1~Ycz#vdfh6P2
z%gk5xLURqe4Ome<Zd#9DOKRjjEzVh-ooxmf^yUpEy`uJJ<p<B6Ir<8CU8KUdI4zRb
zy{vlQlmQG~D79KQGPm$+mHJpH0FkwnLKh4UUIa|((3i)(M)1*}LOR-!qnlP$+Wmer
zK1Cy3>F<F!iQ%(kn^TLt*E=6*;Q*~N9nFTh92JC%bnUKfh_UZk`UvXQlnTLbRE2@4
zIaUdv8r!~8^n8wuDzm1-+Htxu`w+WI=<_kANt7H)UVIV3(y=i$Y`C@pcy)_j$@<zc
zF|daEdAc57$xmE!Pl5pcRw2`HDn#(;7gc4{-gVXVT#TpElqmi{s82$U8@c3_;qyq`
zylx~z^~=qr<KN(KA|Cx@1T}P4xeV}k#bv6W8a913A@nZo(v2NK2yCC_)Z;dfpt=K?
zPt;a8Z4@Ya4e}&8D({bw6Z?DW*2tZJ?J(<m$gHcAzOM)aO>kXb3(rfBYNP9p$@Jpd
zmA|3p+V99yu`G4fnu1eAoGM45sHH<DAQf=*ser#nTR}t<x#775w$cMbPL>=RE!|oH
z0JRe?jq3dcwfDHKG$D;}a$niC5|4d~(b2_JX7c>@DD0X62h*PoD8xZt9Tp9}snEA$
zKFy`g^4m-Iqk+l(_+d=zpRQV<tIkJ*OHVhxgphakG++tHiXgzYq4M9`T}XV~k=oeD
zK%jT*zZhfcUc47WNab&6K=#jQWXFx}{$JSlb<ekeDI%0aD7oDPvfWzAzix;biwi$&
zai3qhZx^t-)j~A8&iH4Sb`X5kRK`isgH>Me(EyLe^j`WjK8}6j$@M4VskcF7HKvy-
zp5ki8CWXb{F|#kZV=`9TcY$i;`rj+tD+&Muh#|2Ype%9HV5}F>Iu%(+p|cp(o>YH=
zaZ58Ax4W_-=XSnG&aSacxqL%$&1_^Wce9B&le0cO=MS*2cMQy8{yXA3nNA!WhF<%J
zZCHPbEROl#(c8&Xs^I<8?tiR%${8^M)k^`AhVzS~dmSHwN>uVUc9ZEhRt#zWY-SBV
z;t$ROHxEeQ{{$odB0!jQ?$BJz>8Iu6BgxLr1}#=&mBZX~)w~`RN+O~OX>5Nrg9SP+
zF`BC+8!o0-sPQ(RRI`?;_PD8Lz{3<DAO9SH9D7H{8Z&nzBU*eaVGIbZAO|%8zWLl;
zIVYBEGf|!l$E1Ph>5#S)i0@t6_7KCuzZlxna|@6>S|K5FKs6U$#e{}_yf|L$<$O-$
z$eRcTYky5BDxQ%zU>5VlxIOL@10&<?s71OAPzNYO!1CM<Ham0*w%)N68v?ojnag#+
z5FG$IM9={chN#BhGR=A&VR6HWIu75#cP5BzdS=Goaj_Q}5h(q|59u8*6`%E}(GU^6
z5UMgbyb)d~k#-tYu%B*}M7qG_pl13$X)}xWU-dHeHpA=%dNrhSa&m&cN1_NGtgDMl
zL*rAKn3WAw5;B&!BnK-i_*x9W8PQNkLd0uH_~_B29@UFPM0}@r{}D;}cwO<Knbetm
z)%lrgPR@!EkR}vDq33sl{vrf~MRS!6hlQKs%F6o6PqVUtXFcgZIwC?^RrOqNiV!*q
zl9o%5esFsK+EyPP{_2#A5C-&g86smzURdi{DVV4N58exo?HjjmEqnF2l`#U+56NCT
zKtr235{G99Bn~)@O-vkSnuA2{cepF2Dk_otZs|f#m3ozR<dU`)ssZ&*UsVANMMuZN
z^$Q3<r*@l0L*nK*Mn!Aq!+(*T2+V^Jq+k+~+mQ`D$Gt=6?}j~;i1Z;h->^N7kC~kv
z3ndDmb}pdnzYjkAo*u}%gQVC#K1a5dbG0vf87d=o|JI+ljyaD#BN}jF<Xf_<+aY!N
zHsT*GadaVsfTsPKU47Wdh$aXgiU)TzAoYLuPJXWN`SZqjaHD0`>%Sj|2-J%K-?CXF
zscCO#5XO>G6iBW-v^43tF5~4T3L=ch@$ncynAHAUcvh*GzV%c%{31{N^XJdoybRtx
zn719awYzb9%y~~_Nea6_yW*_PP=Vy#yXR4cUpiw1eGUYNiuA%ms6CbvY<Z7mIUV_z
z>6Ync?|d+t5NAUf{wzLi*yk;IPD5Jn3{P5t;)gUbw1ZvGWfUYc;lOdt0+3OO<nY1)
zx;&F@ZP|(zhU+(vM<jMS&KO<0cP>wikKfkOv3X!5Ev+{g)k3!MPUF2VC%gA^2_$_)
z78MMOY6naUGRP||EC3zD-B3O+OKOHk{vl#Xm+JmUTU%Q$p942U#FfadcMA38%a=>E
zv}nL>G>o2Qpj(zd*UvGV&EbfK3SY-*rlGm3Zy)y^j4=#S&gX;H$aO+pQ&Ay$tzB5k
zZ1CEW=*|5J>x0$<ID#ff8{Je{$}RxeZZ5sZ{Z(X1r}BPKW+wf?!2xpBkT+nqHB7VI
zJo4=_J>4?5`1{TdgC7ehRN5~!gM+V?=#Zo^6g{fGYA(vG2qX#?0d)kR*QiFG5268p
zYR***>!YJs4eBp~?xd;2sAU?`FXlTh5zqzDomcT@cgT!vk)?`)sFHI>^4+_4x3#s;
z`uhIn_qLk}2DYiHx;g`{hu#=arq7W0Jx8qL62``ty=%35H%l@Co(U`6AMK1tl!$WP
z*M0i{6L|p7p0(cH?G_qJFi|vN+l>$~D-d}~q^E`;^-N*=2<qnM?w&mNA~du*g`@0I
zMAKpG2x=5lz}Xh=7e6E}2a`ft5fv7qJk1BUgj*p9y%pzFd%^dB7e@*Oswk7e%<BbL
zt03G$R_=^}xVy0bCFs8+i@fx)jsC=!-OMc*x3T>5=sk=k24;@40+z-IraMSDR914{
zmYI}^S|(0%x|GnMmP-tS|F}h)Nqb!0vAg`()NiJlWPd)=!@sL#t!XUg`1dx@gX6Iu
zqz&oSm+X<1{RnKXUC5qc%ldUAn-+r>B*p1~Ml%JpbY*_89p^(QW6<>Cj-&S}i>N@V
zfM=j|+f3RtRZ0<(mQH1zYgX%-O7s791AUEhDd>H<Z1ZG7gAc4YK8jp<p?&na^gx(T
zZ*gO-t;CBrlF^>~o^<S2aVc~t|2VEZ9%N0+3>{~SjMj$IZ}&g%7^$b;Gwi)0IKqYO
z9nw?karnm6`Hz01R+->ls=(QaZ^!8jOgg6*5!QD0{15rY){U0FN%g-yyxa0|(ABdv
zFmDY)fsTsLHv)gv812Vp&SV{DXUtSQX{`Q86p)<4R@nY+k5{+Qs?lCT<69EVN%;F0
z-QJ~&Y;<zd+p-l<K~QNF30&HvCiCKzxL7oCoPi<xsDQ+QrWxV$d`PHhHYkj{{57H#
zlgQrj^`FQjjWrS(?#$VuA>)P9R#m!c(MOX?H#G^-cGv#S2II~ngE`;ss}SxPj}4xE
z=dr!B9MzhH?IvvXCzxQ$JKjh?g(uEQKP4&7vMhh>C>KH06qV%xf8r8`%SfECdRJ8~
zT1)3sjBf>Lou{C&#Vrgk-9OXLL(I@R1cGfyI$D0R@GZ>`;fa{g*+<*_qc@PlN9I(a
z@^fclx|b#2zc?RV+}#@YV4UiY;DYHbVr&CgBz$;DROwMX*~2q)8BHLj*jR5uR1y&2
zi*}uXTaf`mxBo48jip2;?Rr~lI8^Ar2Cs5a$b=?tID6AC+}}Ubfy2_5StE9iEJ9Mk
z&f<@`$rT+>JLAT1Ad^JTHT5Pw$RucKg+aQ({kOF}Ktv!*j1_o_ZVx=MqI3)V!UNdO
zWzZ9TTcP_q<P_~nEgZlce-BWO{cP(QKu>6angL{$EkS^n_Yy=N0C|8M;~J2w(oZH4
zdIMyLGxPI7Kxob#%I6uJg$8>7nW-2UM8gYf4ER4_K`6od53>5dB6JUcA3$6ij#oLf
z598eKwNW2_j)+7y`hVuz4~@*uEn(xX&R=!bdZcmtwqI=QIfS{nfB!yPLWSA=yJ~7P
zHS6`)O-$(58h-n@d3x%+g4d_)97J;%v@gJqqFn=k<xRCHE9))L8Hd6Q9)Zl*$j**K
z-}lgY?mKkI<+fpn(6bHz@t*bc^xOous`cu?MGHtv)TNBf%*b8lJRgOzw~jeLlbr(=
z;5;oYZF66r3d7~gXJ}z1PXC4A^|E7`=gYUTJg2g0b<FDe7^Mq7n0YmDnOR@zf3Uh^
z+yd~h2Ub>AUf?zjy})N3j|>JFX62XrngbJMD$E6D67hbZ3<S23YgdzJ&MtZj^y%C~
zW#X;gM+ZNFXbb`l8{wmXcQS5!L?Y-gi-DScyO<VmBhw##1ulq)h9(+94N!e6=(;N)
z!3VG`VI?Icq`NHHJ__AEJ*~vbqW;7Bo_OG2r9f@~K8W-FxFG7HfRM+3i#btO=J@*g
zzHT}l>v2xWp2&7eP_B91s<mJU_kb+Y($Xqqx9c7cjf!e7_Ss=EwXmQE)o3^wM|-XN
z7O?%S6_k)x8p6q)(-K&MKfe)8Efyic0f?#;Hnh$Jt&Zo9Z%wh66+RZ~2%xS(00a0|
zspszfJo)+{&@DJQ@d6`DB`9>bTZV_2o?Y!G7NRh33n!yocXV<Jtf>*{P<#ns8f0B)
zaMUQj%1Qw&Jc_`K3_6)wb38>XR=ET6sr0{9KkZA6)UH#H*Q`Np00_2uen?@Z+5Hr8
zZwlFH>NZf6Mj9mln+k@Oj!p&%?mMd`8%86Q4UxV9qaPR;xWvmlP^Fu#O3T220bEMl
zt5%%?{W=O@CORN`(ity~|CRkXHufRNVF?hKIknH8k?79A9gsAj8Xe;DL5mCKHq(>L
zgYrWJKrAdc{W;2shc-jG0)}6UZh5(uou|;@d(%r5$>lAU#h$NV**%WLcBo4X46RjB
zuu#yCnsNDpXRC)iK)OsqZXE9)e(oV6!2v=IHN7_ZedD~Vt808>qJqjS8;Ic%tNHfr
zn@+++u1fEYLeu}o!~p<*+Ha0g;J;CCL+gW#b!gq<@0$hQCK}mgwx3pxXx{Tliyazx
z@_K&P(A{0=B5+LUC6#aRK}CtAEO#GSex#(J-~`GK?yw>39tu@a{j<fFRyy1{h_KyE
z6{5Rl;FvZczrYh07q4VZ^gLAkZ~PX*%s8YztNgk=$CzC?>478S@~sj0ya=dP$a?vb
zj-Q_tn22^Bm+n=*|Crb^<ma$hcdtuUXnS{me>Mzf{_8S?e97N_S6f+IOPQK7fQou^
zwQ{O?u;QmGw!qWL*w|RomBje?*PjV>J`WrC;sO2e?b|(8oj16kLvDdP8<UU_Y}TF7
z)ROsggW~NF2T6rC9t89eR3dNMXr^{2Jvo%pM49)8>sZc784kWu+uKV29g&-RCE}vc
zIS_&+ahYI25YN&*++#@ScK?HeSBx62PQ}^vUi_BB0Jl+g1)J{ijm=GexF%+1kPyJ6
ztNqE9<(XEn^niXJD$weD|6T@Yi`v?w{q*9-#%pdK9&~JM@7>bjub7)N0_U22r4kdw
zN<FC}WKb5cvx6{5Ax(yU`o*%|-eRXtoif#Hh=+D6>gq$Y_3*Ak7w?Xo*B?QG>?AzF
zHvGoL(eXNHL%!+r{HHXZ-$NB_7M{NIsMnf{-$u)?aBwVoLeYUXl5RbA4iI-;)BQHp
zKm*)5JIKCY3|hlKK)Rv*ck*K(+XOrgqQ-Ne<TOM9BZY9lmb?^N!yMT%;iN`W_1-s?
zl(3LYV!DMu0wEd9epkE7ym|BHf`pGKzwPHhFl97LpZCB^*#jBXz05}}H66t9VB>4d
zsa&<jVDcAkmLp=Mg)j2+^74>(7HAqYl#GlFFf>RS(?1TCOKot;OFee1!D9g*f(EQ(
zD<lA%ApBr<jRc`M<a6qocZwP(Z@AXnu8E9{MA>)S9{$>q!lng*S1V+eh<F`f31Y~=
zXF<-CFdpQREa;E{?<xOL!RudaUm?37rii=&DZ<$&@|S_oy##s^!iRH=vLLoI2LCDr
zR!1%7fepWXshrMUWYEAK__)Ky9LSUK--&Fe)EUAW7vowd({b2-9q0|#d+%!>iq#|f
zd+D)-_H2bjh?@gJY(NNtbz?w7khD+rTa?T=8|9Pce~5^b9}#>!le{Uz3H=GB>Sy}i
zjk>aCy?S*P<ar><3Jecd2%*WTs=5Gy0Z{!7o5DTv8QZ@^6r7#;Kpqp?*w{!q{J6n?
zj#0wQ%nb2r4Gq5e>sERCJhT{ymj2C=8}$jCatmjxd}OPC0wRGg^Ybr&^MjP0Q=fQ(
zBAJCU+njCx1qL`vsN_=8>3B+RC)^7<oo-FnJ=QfLA)#aX(+8+WIPsiJf-zdS^l(rW
zOO8+DKCP^?_V)II%5KRv(Q*7sM|XFA+grN(jBPfIfsOf7mZdu#TPQz`9DEb|?%O&G
zHz528dPOT_CaU5I6(opG0nV;w{@{v$Zk-1Sh{di0!ZWk_o;Z}5gh!vEC$x8P2Fv@$
zw^3mBY#a;}74ZPHC^c%qio>C&<9<ST%5GI6a@SYzG|Fv$^F<O@b!Dp>Q5Ox!53l?O
zmEGiWHB_L08Mgx?S9))hTOmQtycwtwz^R<~7#!fmYysjdveNf(#Q^oxfRm0+Xygx`
z;&Tr>qaEX>E>R3OCuGWSkHIprK#DTqPzTx?UVi=z@OGrHB7}iS5I&x|@<+QjwqyLb
z<lGpbf#n~*evZg_<QONJAV)PCFcGwmJNSPhJ#ps`L0t+3l7Om>IfGpExR^`gU;`cQ
zqFvr6C}!IH-W<WnQO>&rMIL)!;*>XL!#O4$#8Q<`rz%+-ynXfR-9=?(>4D~bi7zve
zcFJg|G+#?X((y5k*+-%A;eZ)@`I?j74sije{d6OL;Rkm%2`I91+&UPt<o#&L>pK3j
z5CMJ?R-XI*BM(`tngprn^zyQG?a%%FXLA7cMEnhSO#Axspp=!rgtdj<YKma8DFv7;
zMZ^`KMZLgs3&b@r-*w%iqa(J3eduGYT2Y7R;pbjfykD>>obKPhh@{n#qYw?p-#`vB
ze+<-!!``XmVvtKe{EytsRRQy$Qx!KC)d~QT5b*^xUK^4%>|0&#Y&5U{h<864*B9%)
z{a!TonvXtr_16859J*D0@$u(DKJ>e2#hRA_be^xPoa2>pT3)H;YhQH9Ixk<OCx203
zzVwStLXr;=hj%g#O7np)VS-rav^qt7$Gq$-&>U2<>AYzi8X|?VwDvNy?h-H=2vHCe
z63~tJ2zZrn8u+_$(cyI$2&%bZbziW&flmb&(sZ3b%rwg5^R1kIBLFHN9b8xHs!#}<
zdnuaW@%|4aLVEOy<7ig^fkcFNgt`v)s1s@OK>RfxOh6NbP|n~!+9^@Ie0-a(Q9a&i
z&dKB)cfhucKFQ4n5+Y&|&?BSk=#-}triJFd|8bG6bt~Nh9kmd>9<WG5(K^k@jgG)9
zaC4BFM;#p;6vDi!=Ldy^oThgFLL~aVlb&s8s@#GA-n!nE=iQpSFXBpmtt^yR@XXBz
zk@vaZ`h{@_R6iL!I-M;Pp*PR{JJUH6(*%hrj}b(35EU4M5Dgm#2OTB4)x&c^)QwZ$
zbCoZQB{OR8YUO!q0*(L_kA74VQEEq)udC=rBRZj!liNP!n9G+~5`i}=v-9lnx;62F
z_uus=7mrply(^i;VGBgP4EoVpzWXz?y883C57!e~cHPP;_wD5Cp9bWSE1Oz1>`$zJ
z=Cc3`_!Ktym^n6zY>vYq?V)^j`|FQsnmbE`{EOTz02akkN88UuTis#Tdhw-O@5);j
zW{2_kj`3oxN;H!hzIpfD0sbl<FNsNy)pa)<+%q;!Fb>!~I0|_^ETR}0?>_|XO2_V8
zG0X6^LUTtYy^z0)W$*mjDoBzZAZ5<!E<M<av*gAW<sP8qkZ8c-UrvNHxbxnh8d|!H
zqS?l<bE4%u#J_%<9xMz``D?dD<P^z#!a8h!akcgyCQ*gVqc`V&*Hpe0@?~u$Td@di
z4*spi#M0Ix6qd^(*Oy@IFsot$DMBs>KMs=6Yo&B0YSGh`1^!Ql_xeZ*QaONTF4}io
zsW2+y^y4hTo&+<Z`-rrjru3mXt3!G83!n@n1y@l>Ih1Cwv08jbR7a0q&4^BY@mplW
zhl^)Kj_I6PuDP)LZLTPyF7^M5iVdt|^R<PP0f9|&oxR8iV^|MU4E){4e57mxsCR$t
zO6%Ue!25>(uUZx3cH#+y*p5v#!5rpvT~W@p`e<^`fZ&8F{ewSL^FX@Kz#=5`B&nQ_
zYOtCyaQCZ!;ns~8&%}=Lr>{wy(ykqlqnpnCVKL7bp{<SZ+woppwfzTwutsnL2Amvr
z4SJ~PDH`Ob8<@fyOTUQ6g!iW-Tcjff*Eg>z%IEs+D_)|t9O{;u_z0#@eftMmv8F*f
z%KJ}L<sXb{zz2Ato7_+j1EinuyI0N^2KvdvaOkENEIIz51}Ys0NF#S%x3H~;;FJB^
zSJa1fsYvC>+Qi>L7SpmI3})&eGrY{{s6g?EWu<R8U*wz#-SpS2zvyft5a98`Xq_4C
zYC}eV8eJ*Ky-xFCZG?meYqA&rpWRN50&Gd^mYkXxAUU!}gWB6`Goh(a^yBac?dlYD
zLzAK_3;o~Fo-XD&5DKW?xq}FLAZ-Qt3I5e5@2YAcVC(AY0+|0A&>#;VK4i#?TXF$_
zMJ0!&Ve;a*yBXhxvvVVpa^$~M^n8imwCGJ~3Lp^22Hhdx4s<*`#FCPd;H87Z!*M_)
z?*N*`7f_A~gpAwwE~Zk4xf`N~+x@Wv3Ii0>vr5vbuQFbu!D%_X+!prI<wfdzj>L{_
zo;h`#`zVM^-w)op-c*s0FJDZN0%Zwb>Sz46r;a%Z@OvMvC39Nd|E&e6K$w-;xjE3J
zhaVhxby+i3lY@;e^FFDF4E_8Re|6s4O~|H%#T@sPSj}Io+|4L*15HhUcd$ps1PngC
z|D*T)_cCWb8Xp=G5==-L5mD#6vM_A~5KGM{D<X0{DM3n8`Y}R6>|d|1^w*)4nYM!#
zp`z-mQ}tL8<$qSeLTCil&3I5(6S<QAr@lE=RYXM=v99v`BD>acTZtJrqstynm9n!B
zxt11J|E^wAu&8ZlkN|B?WF;^<6D|`Fz8Zr@_W~pVkYB<JO<)Mx(Ij50bMLyl(NIw1
z!uG}isTHA*MNeah3hjCIAF#8s6T+UE8@Lq|uzc<5d=&3(hrL${D-m6$7f^WE`1m&f
z`-gCDW@QBjq(&$|!i+Y4DeT<AgPzzg2=$t*XnUPFZ9L7No{VLB)-XJ92^(X-dffXu
zyi~D#dE5wtA`{>v^540dXz^dzsJ{yR$tziynP;+xZmRFCgB&Nk8cU4yo8ZIrUW1_W
zc2XKySwrB;{{jY?Y3Fv3Tr71n?mH#u&%tko|C8!fflH?aE1j8{NnlT0dIJ_7)3r>w
z<F-6UFPNecwf2L5L|oz`R00Af0#;$F$V_er4AACWhY4tbtDL@`^ZWf_GPOfSZQUF7
zyQc(qpYB-(>umkQu<57fHRc|5?8DxYTs!6o0_AD%KHci(zfw$+n|viWL9hIqI^)~(
z_x&FE(=lE;f9X=dCG?(yGzqs9SyKl?tlpjt_Xho~(-yr#TVfLUi!H9^y+YSO?~0Tk
z!%1FaPec^s@IphW6nregX1FLwm%RfD&l!r+mVR%ttVw5)sZ~%Yx1}`oZLgH<H%wD6
zEqbhD0`A_Ko2=b#H(G`Xfsa&<QtI^V>_e!_1(om_V9C%>wLd=F+*DSU9WK%X)R4f&
z#s(78oZg_h&@2Rre_B|k9gr<TOypaN;#YrgFI5CP;9%@T{qVKxJdZf&?mMDsWI-nl
zpa6VX;m*o|Gsx!<tu~TzgP;pS$d`2%!ZD@Dt$DASJv9AwMEW59DM3*7MV$rQ=4sSl
zaI{#mN4uw+8_3)<n-=TpAfH&Bctk?<HEI!_qBymDKFxX&^51ojPut)KiKkLU80r00
z6}a6R!4yjhx83?&#SO!@wDP?l2k9kck7wev(7DPyTh>*<k;YgVwiC-mgMEIzRY>IK
zv#xc(Bt(|*fc!|z<P!DLSLWl_J&Va)G{ym$sc@eqhW^3w#d`T3{~+C9f*icYmb$tV
zXMzmKO$|@?Qd%IP+fNDRP!Vsu2jm+~1D8Z6i~oY@)4;#ABa8<h6!s79a~KYXv38m^
zdfs`+$^KYF{4yA`bB`O{9HPUZU<hpdjGni5P<acQw9EM?aX$a6cJjEJFz9Ab0oxz%
zbbrat9u3M5k9gncQ=`&34WxSX7kFfH4uw3NOLTJ}9T1VwZDo=zLZK#;K9b!z{Zn3-
z1dygWGe9IcI0*D=T!O;G!(FOxWFlf|P@?04@F*}iI0JrjPW>t#vH?6d1N3!uBOpBj
z0rqo9kC1o8lslPLT~Yd6s9=1QVh9fkRV*=jn&<SsoGyI4Ac$P~r7<LeZfkEBAE&1i
zphha~Xb_70xw*tze+~hFl6`fau27Tz0t)hhd2EBPLS(bn406qNJx(GD&T$Q(hNpXD
z!(LEO5E1@={rVMPjTZRi1D*;EQAL?z2*+GN`X@9LNAzlt(DlP?Hu#7^3y^8Ei;E#3
zQ$W>u?n0h=yKwQ#10jU)g|xQ?asW^{BOokTHJ5qM1pwG?-gRC943n0NUOQMH&jXlw
z^ABl&QD7q<wLsn0*N5Uiv&XZcxz66-6|ne^cIi-tgL7h~hSGkmBt$WHfi|;-C?6kH
z<-TQ}x|ouh3RPGy4@aAuPr=I)`sk4qR4GEyK)<!KJnTwvYlNkc2j3E<!2Q;bFRoWK
zoJW<Emj?jaF&RXpvhqmW19OGoB^v5z=ccPDUs*UZMm&tM78@65WNRDghuX5co4jdi
z*%O+OIjlBL05vEeFE@hHw+lSxxPV7Fz=s2{7B%wj-oE;zf?ArPZnZ(bjdJScDO8$!
z3FVPc^)sZ;0qC<A5gC6{;Fy6jqmaTv+kma)7N7SKOnMV_#R5O9-0tu^bQ=DKrI*=v
z2`mLs4g>qm=`bq|u#ORYMg_h&QtX-xsdmfgsMbCtOIB;$#4-NgL}S9hz)`3K@sF6@
zG{*MC6=2!i7Q@8`9Z+TT2$<zL(52fs>ED?uPTm|i4;W^}k6q|zSj51Q`=JvPu~g^+
zYi^&Vp^=flew_#H<;$24+!c!@0GW?}@?-{Te~QW@<g;JxC`G8|Azx*IIyw+h5A-Nr
z%2@a{dzVyy@4G0<u<sB}bH`0t-#zZ8F$iw~;)pe;=b3zcp+(%}*8R!aGqSQB+8PcG
zg#gb$VjrBCNCujHc=T3o44b{%ugJ3y-CY@S^eiYZFSm)Tkb2=K^&*$Y7p7>vRAwNs
zPwezo7azWD0l<Ltz5<#b5&exze0&v?X6Mghrlh1;i-q2#@cf|L@~YVSq32OaFC&!p
z4?KX4kAnGbKwV&CTN@^%nKc%U_(^N9^EnGuo(9FmLZ_c55@m_hv7aM@FNwg9c0cl(
zuB{bms7Mc#9mS;tEbgWS5jLjK`nb)9=|kf73b}AE;VMBz5a5WXkEH@^Y<NH&A@wp)
zn}~dV;of8|IjGY=&jU98ez{q_$ITk=*p6gokvduc+GoKgQAo~xZ|n_T9FsEtO$TXH
zx6{4zv!fibv>2RX17yWwqamPDbLO1u8i3-j^N1QB5IQu}?}J4Nr0B19XGjM=^U(^D
z_F$+3LP`vv%_ni6!ySz<7J-ZAq7|g((%pS4#M_|nHiY^aC^}$;Vvpx{ufp94R?Al>
ze<e=G{e1n-Q7cbP+zX<{0YjQvu6hoL@0Z))j&DH;(u-0B7QU@A1;zk!<-#!Xz6Vde
zW2!_z#*ajvfL6e=Ke`L)ANc3)9PD{glXj#?f;X7;rrPYyy4E=lduC#C55xlkJ{MPp
zObSeFqsZ-$e{YkoESO8{WOr+>?40>!3m?ybN}pm{bs1^4X?tz)gl`3B);1+4M=Y^m
z6AEE))A6KDhRf!GGhdzYg^9=;=o8PDTv@<<HE_yg-?_0S4&CUO7v*+`*L@3K_s$5a
z1}`cc35m#e)_Bybb3p`0L>LM<N2B1Ck|A@8ai2Bdib<BQQVlIx&3XQ6X!?bI+MfMq
z?-IC4&_ooP{G2u5jv*IzuE5@s5D9v874OE*M&nEJ-KNj6zH8tDUIM*tI=MG|%hdc_
z`6DjZ6}fngS6ms-M8wmC2+u0uoCHW>kOcP7fTQGjz`u{Ra8D?#Or9?oqiUvKtEr^)
zWHcm|90M(4=dG+S@rFYQ#8=eR=(8L}?j=*=>ew9Tl6fZW#1&NF{(2hM-^%xlKREU7
zku~nNW82W%07DDge>|zm7${IXTlEg|&+}BKG)8{Eit)2*y4*cr`nlop(vAm1wkS}q
zm!mCW3FG_|wG3Y-3n+^}B}>Wj&#XT)LN>*?EZ?MfD_q!rMS!KRCV*HRT3P7)GbKHk
z(yiO(M&vJLQ?E(pofR7sF<yLlyJ|%1B;fxVMy@>Y2550^u2Gr!q#stNj#Sm3Kg&Tk
zCJk&JOhc5Irpw^9A@J!8{4!n@WTkSZlP&L>Ux=6sJ8n=-fYh-9;ifVi>?w#IlOrS)
z?i+p1Ia8(Xg_An5N>pRdDN1GgPmGW;hP3h%HN!+^q^=_5H8yHA+<e<j=s4fJ1^IU9
z8ZufZN&G8Lv8!mlx*XSK$#lMHUwDMy`{pE>@Br`Lhv+BCOPinZHKnTD#<H`Z{Dw~x
zp;YAZPZJeyq)FavnS4*Exc~OSvyBI;hQrNG8&(|%OMe_|SYyN4)8AFq4@%Ziy0luX
z>Lz%V5GH*vv3injBY{muJER}n+?Oj~aoVh$(d|@go@k;d#ko#2J#*<#j6^hnh)n5K
zwN-+Ibd+buE6%+xuWawd4Cou|SNUg`K4BuQR9?PnueHs2oZ^kn?xG%g21h3P$FV5l
z4yY^W^pFENHUDHVX?K)4rI$50F7mQ*vsB;MnA|@dS)YaTKUbrc*Zy!d$KiW4o#0j0
zQhL_&-Xg;PoUOEC8&>K~DH1j7z;(e$1D>-+cDDW>7D6t4?2!4FuJDnrRMhe}@Ai;i
zRBXVn`Hi#Y?^E<ITj>Z>ioIr5fu-QeNOxq6cCX;e2)6|x*3ZhgV=3w<VmP>M(BkGk
zQOo0V6Xqw9$X>W^#4Nl0P4dF0wv$g|9Gl4}P{(Gd-z^+6{JCa&uq#Df!1iU&$%Ql6
z)$v3-a@-XzXjN2LRBSNRzX^rgefWKb1&3Tb&6(zg%kIesJ5E|bFsZ_F+Xf?T%9HfB
zTyIdeZm?qiddpJU$=3&@5le}og}GZBQA9R|8>F^qgb9^A?%tSrjpml6?D@}i>ih`C
zQgI0gha=^)vq}a5x>#PEeAv;uEB|uD6&SElY#=E8CwQ&7bBc@oqt3r<CUb%+*)`Ls
zKtfH=OC;&P6V4>#y6vwYZdRhDp4<am!#QOQLdUC#@R2&K4&{$wxVQc#j@|KH3F&GQ
z|0|9KO4DpWv^C`&CTN`KT0jnLdJVqEwbdfy_2A1{l;LV%{hn-SF>nPC&)35@z%~+Z
z35Fa|q2M^QgR=}A$fQuHc=_yjGeYv7hG^5-ITJ0EpBlkwLoYjxwJGTVhRSI2N}kxz
zh?dh5+4R4L@1gj}C2UI4kpM}^+9e?E*|Tt@Zn+;EI&cj+Tf)cH{^H^X)`ZF8Ikb1C
zJ$~OqUF-YMr+F|mtjX%2C~aB|O=iA*W9_n`q?;S#rp1V&7JCFWw%e;?WB`an_9EYY
zYi*SVzB$uwsvbT>pv<CXod%&;Ok!e4kzP$R<ikUbP}MK!v_uFD!vE>(J>a=s`~UF|
z2~9MVJt`_QWoDF&QmB*=QZx{<WklI!q>d384M_+I*@;M^Wbc(IyMfIA^HTRY=l*}c
z-(QdJ$2oULyg%>v`?{{z>$zS`tgP3C4RQybR|j9wP%XV@M=o+Mewa2j>di6HO?uiH
zY(F#==GKum5MFVI$k<peVf`0z-(}#MqpD-yAAR?3<jLz{i<MGvaRa}&Ww3Rgt;>Yr
z)?zx}RBB3?LRauFKIXR1vw7VU8@VilEeG?WLyn3nm1Q5tWo(qRWQSPQ6P@$yM5`4s
zxT16a8to1r;ncu`;1B7>ssLAeph_%%c;I##(%y7*be`M%d}Q&zGyrjK@^F(OfI?AM
zmjzNoQL7GF99Uh62{C4N1zZcat~y?MK@-Q<8|yYW$lOZHtMxMAg08dMYd1n{1_xVB
zd=tB~l7ec_<b(Uil?C#N)-D09?jO}MtA1Vg&FlR_mv`XgP0n}2jwQo8EZbtr+N<{`
z4n&Px90yUoX6;(UcU+qoZtX=uPH=6amw*)XarG>UpE5`skOM)X5Fx>$PPRRFJuI;o
zhrg+lyDC4`wFE&Zc^y_l$@)2KmU~kJ^SYXnB1Gh^OII)(L1sJlGi5>l+rWbpY28Mw
z9+3gDJKwB$pEaP(c!w4>Pq%W9-;XY*&b3WbT2)KC!{nN(9!A;n@_y{ZgBOL_?YUJa
zCxsD;0Qf3-`S>1n{6cowQ;TXUBc%R;KU4q!*0}v=<ZU9)Y<X__^NK|Q$VPa(Uw@20
zm>D}=bys>Vs7LflNsex@91BaVT%(qIV6mmBrqa&g0}}9!*opHO=4T*H@_BOO&A{s|
z^ZpJDjjlqL6_16TCU|g$-eEZBT5GzCRsY4Md({cqcypy(CXzJ<!B8rrp>YP8w`NP0
zE(;3_NZ+KNc_9@h@Vg~XI<CEa)4#Fs7g>al<xyPL4$-#5scu{rQ0&LF=|iFxKapTV
z2%?IRJ<m!wzP5(G)nnLqbn80UsKwVU_uktmu--CwoYe6hQ}bdJqxr)MvI`d&DvE-h
zFxXDFQ`$I<%ZGPdDT^3?6o&*xyrY-w?Nh|t;c@o%_LfQM1wlqoM>KNnU9#+hkMiW^
zb7w6Jr9!>Ki!`VX-yXl=W48xt+w#bF9qr|zK9M7zO=!fgir$3R7fqHBh^>Y$*`yPx
z%#^RPsMS{<r#*6`)mG-@RH_X3tf3UQVQVb!F4n-hEr`2<Yja2LMU%2YmBeGqDQCd2
zXFf{R%}{~qr``6~3B$|E&&?Y(K&U;f`6bm5C?F)1^4(<$4bW0bNKA~Oc?UGHerT-B
z+zWU2o^lH{7FusC%N&{1YV0;?49ZDHDD|Mq;Z+6`=AO}>9U12}34>dcdZ_?0J6pGH
ztC7dnvFhhHrIen*hwQRNlY(B9a}47pNt!<-zGzna``Hx-GzLQW4N1$*@bDAn&W47Z
z`T6;*1`Z%?foqYd6!-Th*IOq?jZ92QPUIzu^qz9|H&(Y$MTA1i+#I4>>P)2#(3mBl
z;T(Oaoiq+nM==uMiZA1F>AEew*%D?|{kbRSg$X^?p%vq=&3mWQr1|elTxCgRyK~3~
zN|Vi-Hya_-vE|FMWy_%W+V5G3EDsl*Y)g*9g?SgW2m>*xCFzzj5L_USaT7$&XlR!l
zwtcb(4LNlF4Mt_Q$%#UyBklX`D(F8NWGL)i0C##bG<0d8if7UM^okqILW(fI^QlLg
zyqS*dYe+ds{GyQa#0x81y?*_g01$vF5ye2y<5jms{DxZmD{PxyfH}HhP36nx*xjQw
z_HJbJ1sc&UFV;4w53I<^e<YonoLmYe8R0|`TeDw3(Gs_t2&<tmc+g;nAC3*&7~TZ}
zIKwBf`AM27YfH>-(EZ(bP%SgPtQjsA?=#IC$8aBwpe08l;2D7vjdp<+?*rK+NiNfA
z{B@OuYFIiKf&g3eC@lP9z<!|Je_2bGc~@z~&W<J_{TGaczJxiSQm0=pw{6=p_(&-f
zBB}yU`=NSayqv2SeSCSVRnaP7)5u%oLH`G1^sPE8sM_A%gyiO6xY>wq49=hD&vzja
z;eQ7bRr6lbrO)2@<cw4F>lu~jX`ipS%WvI^NEER<hbW9_t?a9F#7Z5<qKt6@z~;JQ
z#m)ieldZU7FBp(7bQTt1Cs>DfhX^;0hS!@wTek*sgZswiL1&O=@v6T5I1L?LRH?Kw
zT>AIj#x}j0JauYC(T3H2MOJMhPwE%kt1lZG#thUag+rpSo@AReG&GRl(ly#`FCro~
z!giY|mpxEj)C?tC&AokeWYJDf3=1IH${Wr4)VE5JEy*veaf{nI!7@|rA-i3Q&gt5L
zXL2|v*N4IF4dWgjkg5)B&{*kehldN=&0Vfcf8p=$>ER{zaFrnTZRko%aj1FerG!21
z@ma+#br@%xA7mwdIFjiBkB;L=G>?`1l>TyKb>6Pe0g_AcI2g6HwO>_LsY5{O57pAJ
zE&S9QHmL^H!F2NJjr03C;@^JyeK4~xy|5N<*C0pbdI4>`moMxG!y+JaNqhEe4FGrI
zNY5|$?o!kD+RD3g1zl-I-fgq9u5MBdDZa+-&XF<=K5zW8)WXvn26?SvgzkYW9w9T?
zx|P|bTpa+aAtou>D0M0BfaHL`zyB+Ul;9;eusbC=`4FJD13JgQ6PGaD=zFja^C(p*
zJOf~GX1FDD%?Ssq5Epl5s-!}uD<bTSBfaxm3TK~86ywB`=Ed89@?`~&*VOFy_uf#|
zq22%Mew&Mo^N+_(ndaeOeV*2ew!g66r=+TCo*kc?t0$=c;&4-I>kUAC2ael3U~P@r
z4b6~0v4CTHtbWY)1FM9*LJ$KH0c~b2i0(rZlUQ9`<hUC4X1;yvme(P^oo{b;J~zF+
zRZy^Qeiqu*Wp@toq1@s25*<o^sWZP#yevR+3!tVAC_dK05Fj+KXW^C{wVj_|6)6UO
zqwfvY8k~A!S{>UF>h%7zK7EkAIbuy?J_lqLveam2PAT`u%F2>3HM@Zidt>W20%uax
z(b*_JJ+KK-<bH!7LkY_^zkB<9&mc6c#T2i(y(Pn$<+yYogz4y<JPFcVPOb*~G&ENG
z^>ur%<Q(&0SSt|6EGFAhxulUH-*<P>IFd8HvDD9e`52#;nYomLJl#LX_U*@C%T{s8
zvy-?HES#^h{X;{<S23K0H;-FZmK`!7l^Of?vaJ5Bwl@9z*>9Co2m%V6kEEJxi2;aa
z40S9C+sdqkUi7rFajo7KiHd@({+3MhjXTbGO-!y?w~hu;P!D7!2*inpT|1wTJ&B&@
z(){e@+p^39ubRX<>X18n4<Z)6g#|z9Zxa(`_5zx+1vBh_U=0^7>eI!knQOlU0C`#I
zu-?&lY<e*Df&0Qdu$<(&_<>;r853=r#Kk5Lpgg#6eHZzt#5-bVm&1wffdC`QOVl_d
zsSeQ_2laJJf=(3%jlI`=a5Sa(fm?psX+_OF7tE?|hlE^(UIbP%(}W5WG1$mi=W?Fs
zs)1L`%oK-u=b><XbBy!XZGr~5>*Sn&AflpR9Wyh@;2_Se829<h0iY7!%YI+^)UaT?
zy=j%i*|+P4+H%$+!HiMSs5Q$1*-d=bD>WgDJ8aw6l66UtTa4fF{4Cu9{k3>+-MNE@
zl=$9V3{rt>4C*ryF@TS!{{O&4!*lcb5ZJ@=NF(|>s>PDBvPIt<ie^#g9~Tl2CN)~I
zc+<CDS^fn0&(}VGhJ<#yB<E{+X=#j<udYc*R6Ra(cUvWFUL;o*V5i|sVqBtsVBqUe
zffHk|R*393xZYuTO4JNdyT}eN>ztpE)X#Tff{zG<(5)A_pMRk&hlhs>*d!3{F8J<V
zBeRtxGTEC3ZWj?Lzps2l6viUl&B(cQk2iM*BFB}Kj_zo(Xn1%GBqA}dRY`W{_q(a6
zdoxm#M7$z*1pw(JH)XIf)fg5YqS>zN6R&_E{l;)>HuiRof%ZH}((Rm&xl9Y)9<lhv
z3#W!@ow&i6HMeg2lU*C@OO`~PYojwiD7}5~&I#evWp>t$>%Vw2w)7vW$E!Pz>pn9z
z)=lC`;)N4*u?3_-kVmp{U7Wq8@%qw0piV`>Nw6=~uwV?KKoOs!ZsX{G1<~_3E;$ae
z9n?i4q6i*iN5me0$6LAi&d0U2&%n$C<@ZYz|5Oo^<U{ZJ3m>|yRFL%t6qc(^O-+d@
zkfi-lNI(b_TtUe6?S-&aBGLUrpN%w@HKNmQ>MgVdZrnkQ!@AHbR6LrII<;Ve6^c^N
zJ3G%)Q;@7JPsCG$4JkXH(9Rt?x8t*=oc}j9aTPn;^75T@il<F$x@NOGSAvJ#2bL|f
z_WQSQtaAJ<O$f5<MSPMEqHi5WJ65f@;Zy!H0()ff$B%2Qt*u+WdsM!Cn`YH1zco}X
zfBI|EniEaCaO3F?+Zs+<;6OuJMGrPVetS~${8TRk$-S7u_Zu54VOT;zdY(1nj!35k
zEI91Nt6bzV?7&KZ@U0u7_^Zf?r+k-VX1sG7{(sOK?|LluA@>miY`X2C#MJ;jfCq&2
z2-#3p_O58RWuqu6DuQ4+Y_5m03tP>WCyy^9TIbmI_4@kyX6laC)_X5s&QUsE9}1wy
zHakX1H=cDLb=%mNWTPe>an3e^(K<S+Dk@%ch{PU?IuGAaFZUDk5k$JHa$csv0m!#d
zVT&XaG;&hIJMJOZeUJL>cw#$y`USf*39@G{UrK2yayM3xXnUlcw8L%M4jl`XZ{><6
z%pUli3$U8w8p47nq-|x%$NOyZ`{S+sLA8Xo?$`4zsFF|vioex7ZI_O15b++IzN<N2
zLCWCw2M5RI5l|N|?jIk&#b*_p547pB6HO#Xe#z0JM;X_yeK4nFV75(qOGH5*joR^g
z`Y&y(PkyE?X~^Y!fk%c&vJ{JrFq@(RqXk2b(0A_N^lgOfWido^r#TGkz3H<#bVp42
zt-9G^<JMk;S3#l8PO8@Y%R5oSW>)?}=n~SPpFCr$YnqT9FJ0@@<IitaRpdBocqm3z
zcSfDf*fw>}V!Y=yiGA5&-XQqnYyFB&_XW8th`$2erP}j$u(4^vZnT2~Ym-#opp!{V
z*24T`k(BF{C<*HF7f3b70Z-nWG~~JydoaYEd#V!p(Npp{3UZqXuvR+3-(UC5(ZyWn
zU`HoHv_kmYxa@emoFEJEQoP3qnnjWCwm<se3!_^s)x^Za2+;+GdIOG{wWuM<^|s|-
z*82Lf)cfw5nuW}T<)rrWUyv786L?}K=Z<JFs<>)ncsig5wrtPc!l!X}Yx1I)_6D+c
z)+x+SJy;iY|8`jxKcQRtDI13++6O0XadRaf6f#WXu$w6OA|xp2g<3ETx$JO8_>>-(
z95QURMX(DY9#5nD16AjZ-jPt0&nJ~7`IwoRe_B=!&sE$kr?$<ROwGCIyRA#7anz8d
zlN}rZ+7mha<3XiPrPs}ona<yg5<0M}w{PD*c6MVXsw<Oo=TKL6ej1I}gMA0W69fru
zd_6kTb1FW9In?9+Gc_RIMyfjW>=wz9nNM+kuS0w#Y9hY5%u8acP`0RIG3DBQq4rkf
zv6pscQskDD+<>2i#lq)O?$Nr^T2M6n;;!&v(WUHzxhwCCJEC~u46+~W#y<?oS_BJw
zEUSV0)MQ(oN~{#6XJ81{>$gz*VsmBNIi*G##|4e*CizoET}b-#f%%3?V6$07US8gT
z7)L`cbZUr_Csq2<@H329EC=3y0B}$j`L2i8T1h?R$fdQl_cgkvtJHpsba}n+_kGzM
z_BTjw^0g^WE6V1lGtl-eadF4IV}?19!@E^~XetyQCEi7B*A_%VOTf^5n?u2Bdt#1y
zTgI&R(`8(K!sTz!K7_c-EM8^y*v$6U(tHX##{IH{bN2=xfX>fFo)UH&V7r*e#Saxj
zbGkMwY9DBD_m>ITIb^F20o%5F{##46h>yNo#bzgFe);Y4Lo>Umbd~$YEYABHtk~Dx
z2Q}p%09;HJ-{@jrxVDkrG3uwxVGkbM5ONS&sZ4jiwy<0bs3HiGH?BZXGwkLbE#HX`
zueO&ZKiVQ~x2vd3PRk$N184-L&jG7s?mz>i_AcQ6Lfge0@xzkawx%u{d0wQ}-TinE
zxn#%ee*wQO)Di<G5^9#mPd9|!o|@C#(#Da@lqb55it~^Bi$J<ZP6BCur)a9Y{_{mS
z|L@XUE}XXtj9x+eASvNr8iBen!1lKtwT=hov}*83&(P`*B{QyNzRo$kC`ls(k&*?u
zL#703FB)D~ZNgkqQMP8s>3aji!-)38ttRRR2pDnMChbp-Z5qm`mAafsAGG}75d8|}
z&Y<dlLF6aTZo@2;sv}3@Ri=V!r-s#dhEjyanr9hlreFL{q9$z1^mbuk6<D}Qs1H((
zND3AK^N@V%kD#7=pkd(_82A`}qqKfaJ^Xryf}ZFFD(|lDx-f0uG<WI}!w*}e7vOCo
z5BSIW|3e@^lF>kk6JY_Ydl8aWocJ4Kby}JS3?xFJ3lZBETF)xzf1p$(e7TYFmU7Tk
zqHl_#f4eE6uXvb>h&nQ`0xy8hGm~YMxcq!}9D<ZIn?ALz>z0#*0{A5$MW7?~5USgk
zjg2SXe)NCtv|8p)!1hd9&2qk;!+P@T!fY<y*-7uYr0pLLXdelI^I6UD1ehL~GeT?v
zPoBK$C~!%KPrs?9g_?4%`YwWCPaurzz;Ur7Al522ANzl+UG_<1|3=2D_bN6C0zzWh
zmlLFb@*AGp2fe<-6)bN+xOe4Tw4dDa|B7X%q6y)P@Q7fc+yZ~Hpvja;!=8Yl5uU#v
zefy20iOD}-%-R3U0}u$y!D-Z8QVv5Ku^s!7EMN|qb>rjX=)`C!FLE0SA6-pKN}7=7
zHO49rh}7FOHXD=C<Sy~QifX>RDlezQNzqacT%GOi33@7`r-RuAsK!gcZ(}565>@8S
z%D$G%BZi+RY?aD9Hu!IP5MM)|JpBzaoVWk0495VnmJ%e@&SZe9ca;B|pMR)bB8Px8
z0t1YFd(vF=?c28xndT>kzr!q&U(d?2zR*BP<-h}`sn}GlwElPEZ_d7?)AAg1D*9*P
zX=zDGOV>N^s*78K$Y%J=mz$|qzJKovmmRX^z|IGGdC?5p;zdEHfM5jE^lAUd>#Jl~
zICyBZM~gg~ZtLEx@R)o=b6P(x`|p>n?$-aoy!+@2tD^73t4lyG9vWy^R^;dh0W93F
zFI}_Hd!nu4F|znRzZ`7}xe1!3NopJ8T3Wlsi0j4mq1OVNRGMYG3|v*r5)u-Sg8-@_
z%6Uk2P_iPMxCXZXO?gR63oE>aiMHKRyxXwbJwWvEw5Ms4!|Zm7QAz2Qh+S#7pDPoR
zr6XdzPLD|aQEqSghi8{RZ*1%-QZ%~;a1%Nbm2As)BOo%TGnGip$pdb$mt_lvr3wFk
z5falE*cno)YHFBJwOfo`h6BvW<$K+7F%4O|=l~cX;(I&CStuTju|uABJm)q)1zPF1
z5<V`ahr1;PHs>%f#Bcg`SccNomtG)x{&CUsrQ{m-iuU^ZMf$9gIbde%O4~6vzUJS8
z9ee@L;ulmK^~ZJQ+5_I#C&JR%*{S+1x!C~j#~LP0DIV@}A2~Dwhj^`V<{aZXCo3Mh
z<9t3KL!`4|e__DDAoBX7ngRKBh<WF$w3QRfrg3jYFMQwzz};{-Y4gsk;VCYs-qvG?
zL|||zQBZZfO(s#8(eI3jd=;@%UQUjqJ<qWmm`(cGN@gOzh0zvI3CXY~&!B-idy{q!
zKr3n8c3kS(voAU5U0aaXq5S!(Q4*DBSNzR4fFJhe#}c?a@AUJt0;&1bD1k^40xDQV
zL=uulKuJYK&44xJz!{Z|F&(RN|AsRj=uNV0exf;4A(;Qo9&8%aGE<G4<fbykczHKL
z#r{atJD>>3-%zAxP8Loe;}GCw$nD#ktj{Vv@@1(J`(gP3A)fX>MmxDsGl(n6nUqc3
zMbkP{4dVjCV8HR<$z8cPYws_B*EP7KR;GDJY$@Gp_C3zdRiS;s?9$gj;*esn=dIr%
zhMo7%4%mmhz}5@Qs2A<KOjYX9L<HG9SNQH4n`M0Jqj>V{H+C<|=au6n0-8U>-jZ+4
z`p)TN=(e-7uytUWc%$9kckl3?6zyJgIV2@zFgkWdJadcVXeZd3<=*!@E=NT&J<L1j
ze$DWoKGzF*aT=A$r;0eqj^rC0tZ92H=Fpq@j2s;%M`lGz?LN2hmfLOt62wc+qM5rh
zGr%21b~=J0hnmut0DTr9PF&)Z019kW@ZGlPQC4Yrc|yZ{A#yyO`^_~}ZjP8J(%0Nw
zIyu@IjIEiM;33={d5$7zmc~BpW9UYo*$~rdd(v%Gae8Wsl-=;1IbWn@NgAzQ@J}$l
z9oImAXw^jVwErA;IGYvp5(m^1kG+EO>LbX9DCWco7OhGr9vTRb(#}<}%FYaBw0yx7
zE5OaXDC@M(vy3!;+~H$<`lR(AL6iNy*f>f4O^<zqXY4bIdPaJCJ&_V#CvyS2J#q+%
z)UU9xC_VMS&0oK|U!BsKa4(h{FAP%QLKfXI2wx`8TxUtuPJ-tJD$=9oxnNw-)D9Z>
zkUNfxmnmTRcsGE%T|93u7xaSicp3F-KYsj(jki&zHgQ<`uSnncyfsHsm^N~K3JknW
z9Vk!UX4;)Gtj14RtV~V)p=#q16I&B3H$l8<M3z)qT1uk(kA|z7h($07RdP@6<<Hfc
zeeEh3bU-Oo&*<X9Y9~h<Zoxy`<XJtoeiSnfDCTN(G>44+R~jr1^0oqrEsYd6e<&R2
zNAD2>0oX}o^75Qe<r~nMwue^0Q!SG=_@`+w!z0?_kJm<*{MO~%*;htKL(>JWVC&Yc
zdr+pMb#uD1GlWWX)5`L=RWg|tMQ7KYDeg}=teaB2uRxYb+JcgNNO18wgl#LG;JriZ
zw;vvqBOH%hOtw6*`D(`WDl1k|Ih2;+$A9EId5Gi>A&MHNOioUN16Jq7Zv0EWBk`V{
z-aFp=@6$s`Ddom0nJ-&V-w~LB$afs9_tg!d!~RE7^7v=^-h^`pE4CF63Gp9(@<e5S
zzk^>)$zl~uhjN<NG!6g#m)r@csQBZOAFHIr3S?Dr(Cjuhw1x2bO&zLWzlV_JL=KTY
zCI?L}DpZfu6<JkmG%ibt>Eh4VKyMNGX|dJdk9(ZFW^T8%@yslvREy5PD&@^`5(5u2
z0<J|xFxyi%T4~=8&R*;Aj=_8F;`jMPepbK|HoK1%^@C?4I72<I9bNC7<*0N>XY1YJ
zMKj4_>1ra|6jv*Ec0;JK^D*jge{Hv$w0T**^xmb+?6AKJM87vL@WE6J1cb{K<kk0B
z9$1>|)MVU$t@x7d;eTpUB@<V4E}rv(Io@B6@ZYrUWtPcn-Q>`1<RrD29kHm&c_W2e
z=-jtWKy<j4#$<hHNp>^`8U^o!+kyXj7G-!A^&zM9u0kmELQf}1<7t-PDz1S4lI}3E
zNJ<PGtNa<Y{hgo8Re_~v#~iZuq7fJvc=dN%!rKe`%Jb`Qx9`|FP4n@))%}9otIm|r
zl_W&}Q=~%U1(gUD#quw7buFjj$Lx+l-g9^=r}<w37!w)jI{v$zNH^~`qHzL_O`_OH
zO|E33*A$m+9l!o~M8?m})%}0PMZ2kAa{Rd%$Vl{POd(%Rdv^G=A3JV5lBKo3x{zwa
zZh=xS@B4Y-|LQlC9{I0wW;<@E8~CeLqZ_c--HU=u=r=k#H2!H5nC=Km3|M?fQ4@b9
zzdP{ui5p#y(@KNV4*wV)`0Zp^>^MS9jXgWcjvY&6D(u~LCGRLbo%fb3U+;g(b^x4$
zoO)BLt7bsS<b>4yhQn{Edq&pF^WSN^<xl;#0Eh`F(i4dy(ikA(6M>SFHF^?Z@B7u&
zmEpXXaU2QY+2^x|@3#1`vXEf9>*%P|mG}s9PhR1kF#V4<fpMtw_pG0HYG==GA=&}V
z2jCVJWh|PTU}F~6C-sSQfAIC|uRz8x7FR7TyyD{PNt8C?ujo<ZMjK|RrCcrd3coa~
zwqx-_K(BAHCsfZf!Fa9pRw1tC2c@&#v3M0zSuWl;6Nf+f@!0%GV1a57dolRXF-(-G
z-@Sm%=*Y2SUZAZ>z8a>wk-2QBsO_B#(1r+1#uq;|N4eqmX663rWx0&h6i6)N1xu2T
zLtEq)?Zfr)X~|nf9|pdkDhitBLDEbiA-Qg@h0!4<dJe-a-Yh%M5j>UR0Rope@4?nN
z4cznn`fjwZ@AzOp!s584q&DK_PP(R;WztqU^7np-SM>Y3-E}OOYGcFkY|VB|Ly|hs
z5*Hp8Moqzf`B1RLwc^tyt{L{BnSg7DPjt7$*l+y)*)%Z;g%QVJVC9<coamI;f>|D@
z%U<GCH&bU=zWf!&6t)dKgLWH0Bt#&p=pnPLIy<`ZQEm5SZy!_W)Gu}1u-s@{IM4Y`
znUyX-Z|=?V^z(bT$<%XXV1S;2%rx4=wU41}FR4hpJoy#FGn`XPA`Dfz%mQI>Hhd-7
z|DC6<-f6JHa^v^qo~<_zomSw!A9jAgxOn&{6eMCCV3LTPJIQ?bMn(=|sQ@L{oP}x|
z?tfZzRwQ(OE#ILFmpbgNo|YGB3Uq%M-AivTZd2sCcKwMdTC5^zBykqGd4Jr!aL?Mt
z=Hm<(`@UpcLPFv$R%G-=m*5gk^DLQfS&%5pxB0NuT%4{fmy7YppxBS}xQY|wd=(Xa
z;^7lHGe2l}ye4vdNRG0kRYxgE6oQdU$;c2g$^LWg1)8_c!M>i`H!zEJ4Oh%DPQMGE
zmEcA8QP=UdssKi6kp=#x)(NS5g_kr0kNPLFUHRQoXY_Y<QIMC6KCeD9210>;{qK2~
zU0jw>z=a0E!QTbOKB+dYvwBM8whd1>71yjWrEID5_gpb0rE^{<tR*+W!PDVRm8W5{
zmv2c&)EcjHXP3o&=13Ix?Fu4l!R!#d+O=!eFs@m1|GRT?O3D#KLql^tHx>oq53Stx
zn5;A|8cAz#?=Z;(JiK+?t*7J$ZTq*H>`$As>#bhZ((*7r3g?Q)S;(_KLl}7buj@P_
zlXCgDY&l(gZ7Id0KIsGj5nFy4BO-5NNEaL@ikdDYvHJxCl-|~qPmE5u2W1)W_f2cp
zbd)PEmQyU4lPK$Zf<ck~qEv*BW{^HiF+V4-M9923(enf)h4i~EY*lA>NA?mkE!Ji?
z!n*<BK@+9e8mAc{N+oI8wg>Mg#Rxxo3gG(@1A|pOH2AP&Rv<~<17#?4*7YRkU_;q&
zNRc3mxdGyAc5C^6+ch4QpzD$3`|<i+@s&j_T-90O<_p5DHTuxE5FZV=&F-$Q4+Hrm
z4-h6GAk&a^g9Tr;XEMNQ(6~MF4@HcBsD`J1R_-~}28o-w=gd$38AvmLe^C=@xUFri
zsUaRZ2D@Zb+hMOi(pvju6*dZZJd_cK15W*hhlhuY@Png?V;Gaj4eb`x`k0^k8mwLB
zAW2)h)C}Mg{8M`M58w-+)5Q}T$_)E0;UNT#tr#osW@q?FT7FuZgUJ_ujarv&M%H~B
zZw*`S(o15b-XLPA1Opig$C@15^H9=zWq&u_qc!Py_bxm3wEKWwNOL&u=HaKOYAbf{
z@nd@T*>O4uSRf@z(0+XDT}3p@!nY>Fh@Ew(TCU*%+52hVM|}3CU|BMAq!dj6k;Ume
zzYs#<hL-260pGT52O!Q_1`c$G(Mw9rQ-j_3+5LKhI$0NWJyU<7by!Y$1$j%uCVu{e
zicE<K8As^w)>4!Z(PI}`lp|4oikXn&<bO-fKmC}~peg@LNe4E{7&({gpcUw0dBXIo
z&A<rso*Qb_)y;K5lEjQ|(e7*A0Yo}!5)}ci|2onxw38$KR;DDDU;Ezw9ygmpwTHq;
ztoI)WL;w0~S7hL$SNd7|mOGxC{XA7Uc{dik|9?GO%gm$~swygE91)aAL^BwK&I)V|
zZNG?+&})Q&ikj7Jh1!D*P7ui-qsBcjOCb)nDO?k)ez3RqQg|VB2d!reAE+n(G)!Gm
zch_oSDC0c{hL{VkqCeKT{a7x+erO40Wmcu?=gAtGMj3F2&<~KzQOF8j*z_%jKOr5h
z4s%NLI&pi2jZoH+nHC6rAzWC|?6(z#(_aNHw&!Ms3h?vpGy9tU5?wzAMiANuFInr`
z`df*4J9`)YalJ4dF~^zW28xWNy|iRoF;HMOECa2~cGn$bBkVSO^0Q<6hF>2AT^6XD
zajssk{7uLvrnPHpdosU%lfeYta*Qcp+_>?=2j`I=&EpDVW!Z@okMo>EJ_k>^H$BuV
z$j|Q@7<fFF1Mra;8KK&^QGfACWG>gCSowkmQ#NRj_JMw~9%;+bP7L?cQ_%`u3P_Pi
zAW=)eol+ASf6G5_ckjY)R}HS&+`-cj+&(5$@L)6=c)O{&$xZ-V4{ExQoo`>=4Epqx
zXC=G5>qeqDgj045I;Jx)2!VZcjkX8T-HpW<SUB<Aw3cwR`9D4$#Hns&m&f`959d4d
zLVyXV5SnE6V6aj3kRqSq%%gL6!->~M#nDtdg%YQogj30fOcVe-hr6e8qu%0rWFA*r
zBGht#1b-S%`T8C5<Oy2(=<zmKlB^xeMBM$VE*}IuTZSAIMgXf&sM3G2G*dFML7SD-
zyup73up{!{ppQiU3@Q#xS6O(8&XPO4ab)J;!Gnl}C-bZHYB6k@Ro~sxQWyvm8TbJ&
zDk)7KV_q4~@c=9}vwmB9du24I1%psOyoSIA8Sn07xChdRVqL#2LxND(;K>EIS8Pu*
z(;JU$4|V#Qn2?4H4gjPQFVv@-+{9LhxJOe{bnM)TRf+)`LcT`lD6p82I^nb#`y6F;
z1~FDbN3#>ns;^wi#v~I~5Q-6UE)EC(+OL%IlChjxtM0D={<-?_@q8!qO7QE-PEJyA
z%XT5p4f+J&ucy)86EtSc;yqY34G9H-YKH4S1as{6fp&!pDXOo)zh~=x@0eF^pBwp9
z(2<gF-AjuPKy<OFXnLX@JND}9vm<(QIquT13~LJ$VG+^Wnz~;>#uy{U<cXU%aLdT#
zcARBR%DmU8n`5&DabSUTT&`&~!=_IQ&Ohyb?S?VH-+k@No|L<QUK3AhXoUjF3#P0d
zy!1MSuF_zA03m;5PV8LwWoH&;aU2hVKKG*RO3_;G_vZP4i=gOn^$Mqtn)_HZg*HJ#
zmT~^Q3OwFng1XdkpXGmS4^<;j3>}~t@tGnxREqSnx3kH+O2Jq>C>P#bTug<t3r6K6
z@dBi+7XQqKS~YsNqY9Q|<{-{%=t-djs;dT`LWiKuXtS)0j1uEuxrTwpHF%Qu8n+1w
zl9@8g@dklgH~#glZwk6K0@&mKZW0nwPiSHNW`A@wk+dsPY<;G>L;0OeJ@XPXJ1Stu
zCUDucr7K8)M?U0(&h@62KUi8J9KaA^YH*!y4-%{yL&j_Sb&TFd-WsoS4LmEG832{t
z8IZLtpV=S(S{NI#s$e`)f(LqVG<yMnHO4#yICfm7Bg2;=jU9Zq>0h;qX0%Lwm-9Xl
z-8wq-fhrFSRKtbEPl(w&IyP2^V}b=z#G(fa0Fogsn9a5Do*tVshk$^9m#;4d<Z}LI
z&G#Vb$BPjg@`r{__RhOiuWhX>FcFsX>y5^kAuWLyQF`Y4ewRWWLo!uIY|AH5XqCVw
zI|fMT%!jzWLI!+HnWFnR&Z~zCvuIY?rUoR0Sx=3f#%Sc_Y!Vk=Vudqua&`<XJW#3=
zJ*TvyreEF?5Om8-zKo2_b?h<{QlnU>gqu}Zn9TrV$%GWCma}BC8pg$Jk8v519-aks
zi!*O==nL|qtmmffA%*7tWRpY2dNC^c;W;nQ(HDjd(VOvpwG8Iilwu`&ftyr|KMIU-
zIdoX_bEEE(xduoWV+7zVQ&inP-|s%JIyw{Q?kIeMCVwhr%%$5yC#OI8aFRG}qCGZu
zUTJAvEPOwdQ2vl%Xh2DR9j?!Uui~oe9Icq32GCm=u~eHOZO9BpQE=cZtN<6jrtf)B
za{<!PiJ@&|zz_cGYm|HT&T_^%TozhlhV(ojnUoG3@Sg|~;dOcTA1wgJgRN^l9bi}G
zWWZ7>Z#;w%8aFFU;F{bAlw2F*%l+}8$U4Oh?sg9Ez9Pgi(wNdep@Y_;v87XB`*v3S
zM~N(Z#qF9?f0}*3vsQL+$dYm#x$A_b?1e&O|Ni|KF3qE@yx<<<D38Oo6iwP|tWoo-
z1?=6+pvgh$1_gj%DgF>!&I=;;Lor6^eO!QPe~}Bv!7xE4n;iW&CSt?Ywt3w};gAOo
z%*r^{X763FdK-Jca$)IY#;MKiJ6}W{9mJ?iez4=ujn4pc$6Dsa78q%eDsBF@38lcZ
zj=7mp_k}!~P?$l)`*aYot2dgCe7<7|jbdT$YEhBnkUv_A_I_Cq7w@${EjSa9u9Z|z
z?Lp+1kVdt9(Ucwca_YKD_Z7%qrrnEE*)ji-Ca)5%tqX`_%^hkK<5}r=im3Ux@Sc?6
zKp%87j6d3q{k}BTeSR6>a=cU^!~fu#KZUh$NFr&h7&DnLj*((kJyZQpU3)<t?G}1A
zA86}Q*GAcCwmixYx;J|_%ol86oSKlR{2R7W^<~$*XFE<^HdjmeuQE7u@=n6^fpjkH
z1g0s*SK)ie$gq9TymhI8*nIU0i*2EAUCz0ltZfS@+F;y%J^;3h4KYqV%PuxQQ{Fv{
zHCKj@5lmJ}^=w;ee92ty&o<6MsHgReQx0#9JG|9Sue{(t7SQb@KXV#$-T_<&Lq9e-
z*~si9f--`*I#NvkRM?4ps;a8+4V+B4JlgPA;|{M=!t1*dIs6|VT@f(cNj0jmx=tF;
zF$o~$XjMB+N&&9i+HA0#+Ca>((fv~kubWuRO8(88A66@0|Ni~Exk2ZEMq6s45V}aS
zi71}@N1@%jW_Ya;hVcqZ0CQ=GfCeI!S%@?9rS|M3cEj%2Se1A^&+_T#<>0l9OiveP
zkD+DheH|bvcrJfjC*NtuYED@{6hc%Kh+4=|`yXWf_i?cH2ZA(7U<Z~9(l(&nDJRM=
zlwB;(3MIO;0~Eb16qX*(q-B5<{nK-*MYCdvQpC7VNr0p<QU7B{BLQ#_r0u~dL?~cT
z{@u#m#e2IL29GMs#Um9SRj4M(;1-VE4qp1zFHN8PqaJgL^QDuhdcn{T#;%I!h}Ang
zL?gz_Ore4!mk*T;1a9vKD<ke<FF_#s3d}zDh>eYPa0XTbtilW59KLoltL{VUCB#z_
z?K$%G>vAxy`|3Y_zFQ$!)|SG5#4xeU&vK+m)|Pu+SnjFPAT$*5qP;)yAuRt`LF>*Q
zUZI<5wjRA7ajRks5k<U_m1rHY?z-SU$2eNtT53(XpMRCsNAeNd$4QPg2dWFN=A6{L
zpL{g%o<SR<`0kAo+JC#=ZdGGW=k>XnNfNsc`@nB^36Z4Pe4>ZsZ5+b^xdj9c-dFSg
zT25;~{oi7v8rwx+^fPQ8n>jFE7HUyfA^_De_<vE$!bb`Kay{#`4UJ%OBbH5?nVd7N
zZF|!)BnH|?fMo5b1~^=zMO@^L9pUqqxI~%h$<bhV6&fOA{)aBJ2W6{}_U<FLX?v}d
zH!N&_{SIieP15out_Eoq9o$U<f1|f`MD^4)G?eTH&ZV&97m6iDO!Y<fr_6hdTcq9-
zZTkpAIy2uhv8>=z709`J>N-o!vaQ>blosVH3N%3y+W&xxnRn4<fSNUR*PgX<ndEFg
zcIWY`*#l}?144fp3V9tp`Ba(pq5dbaOGah!yaQjiz8OyPJ*p(drNki;U`xxmKA(HZ
zHEo*p#c%IzU0u?Z@N(T3#p}y0goOGeJ~9hztNXaZV$H4mPOiJ>SXvaL9=5wAythi3
z4T&y@IzMU{Q+Qw^Z*Dd+**H@V)5i3PE)Uqk#cxi%D|;2`6Z#bCtK~JGl~6g;NxA3X
zjaBdwuxH+^L}%RE7D9i{8h+kwoV#NlyDg=o+r8>TR0iaIL=l9%6%5>>TpDSUCTt%5
z9XursWn-WUXr<j*J6Qws>sTvRpAs`Wzu~EV$-7D0OF0bzR7dtU#mKDDb7v+qi$13J
z`hVJsj1#!$&t<Ui-Cl;@4tX3BQLFGTWIo=C75g+af{fWRnggIj3fQ<q@zp0K^>L0L
z4#y_U<<1qxE-?;dnd4nDUXdR{zwE<)4l;70tu%=LRxtF>ymE5wrn*S<K+FzDrIwaf
zI=E)p>466yzS!xYf`t|kKqa5D$$e&Pj}8rC?fTwyb?L3C5<J4Pl-PG{-XC5*>VL<G
zD}9#M>;Knz6oarB36jtr^cXN!ynFYKZ*4i#fxsm{9L5gXFF9A}zEa_o#h!~g6nS+v
z?=kMbKjvTEP?%{<<~a;s$4K6lI1TXPE5eC$uVze3WU~T$oJ&M{nHNix&EaWtIp@hS
zRh5?gN;Hn4G6V(B@#oNd=vq@#19@+l44ghv2kz1xuZtAzonzO3oUEF9({jpgerwyk
zVLp2zKCxN58!4SK)N})!Hu2Q~x}#hrMIxv6J<k##mle~lOJ+5xl?P!cP|z!?T)J83
zs%y<r!QB^?IqqM}>#6<qAa1Z*<lo<F&HKHMwe{0%uMJ!inVI%pd9$C*KjH%=(U|br
zQToH6_24=oIZZGN>`QITo|cAcv}<ti1|$cL_ko)!j&s;a9+T8Bzbn*vZy=o2YpL1I
z5Q*iEK6GlopV}w0@{@whZH}EJ|Fg`O>xv(kY-32+Q=F?Y%7c*E?&#>+12r*l<TX|9
z4Kww|kZ;9~Tc>hZ<)|&!Il9u&^^?!luTh=~rj*QV1_Wdctsd!RNt^|GS<4&)Y_>#>
zEj!ia8Pl{BVPsoZe=s<FKC5S(p-9!R`MBNW<m-vLkbw;>n}6Se+j$(Z*i3Mepk!N4
z!E6UF60w6^gzUCzBT*yh5%(WDWGA42wZcX>U~$YU{mh7OXRp(N_1E$oit)WRuHz8=
z`+Hesw;s{aY5H{67C}@n7&2#L6df1jmt&1Z^O?ur8%snT(bDR*L(4GwlFFHFm(Iog
za|3j)j@>6{e%Tji;oG%o3H<qX{4#}*4>I3aik4fWlJo((y3Iopm`A6RW<*646e%YM
zh1;_*5fPD8<?0EwncMg=lYIl-K#EDOWy!3CLf1f%0MkI@(vw=M+(#M>zUts6e&W-#
z&)0M7?Zn`NS{oj(-t;IrOxA$&gBSj#y085R<K6~=E=M*(vL((V8|HEkrC|L#ZMK*u
z2lnRksV`6UdIrziN{x(g*vzZ~)49Jvwt1Y4BAlG`t+`;QXSs`&iG`)vWDc2(WPAh`
zJvI^iM#Q$NvhvcEv`YYKMY(aJGJ_XqcRM(|-IqMBurVibtjLaiDS2M_K~gF!m#jSd
z0XOGAEC17%AjFF7S5s5Owo&^H#W?ni`xVph_CzSS%R`KBBzA}V&eY@a@~AMm)^6!2
zjqvGhS{o9VylJ+N;#iohT`Q%dPgNvr=6X|%{lUn}KT9yNKt>)yT%33WTp_860Rby9
zsP+s3gDD<}3cmOB>1r~^<I*Kj-1g|K1<Vw_Ms&Nnn86;sgl<46tz}4b=uOE9Z<PzK
zbxsQln@`)hM%1Na%X*_e+&oTtT{zzmiscD9W7n=-1kbkBmlp!0c_P*D=NntBwaFv?
zQRl_J?D;)#5Qs%s{Vy$vyq8)Fw|(V)^^N^Ja=P~J=Gk%&da^c;@N+j>UKGd~QB}2f
z#d{ABljMy71r2te>%29XKA|Gpz~Nh-+dHP?c(;MwywLN?+9-pr=69|`6(Xt|5Lo<<
zwXB%T&)P8XCWv9xD&-4XHE&@4OIG)0>!Q=yEAq#hgQsLC&rdx*?&MN-c=q(+THO=m
zaqBigyZ(@dJ#IlH&!FqIhz9TXDP_d-c0L-kUNCU{XRGGrjqjpxcra2|POR&x9aXz`
zl*&rI)>!S#^FOy{2n6uaTN(;`HYj$*7brJ{X;FtRcj^sMl9hd)@b9mQ)R4p17}bZr
zo*H|zl=AZHr@XN@X@(ztuqE<l<FN7D?=Le<F`&OHStt{z?z^yx;`FP#I6^)@yEKtJ
z1%t&O^+aR$o>MMef*sO!oRqXPX4L0`uZLhJ923)JYO=5X^+ap+hfaAX=0S#-rIgvT
z-}XEV7};dOjUV!j-a!^`lWqcAT#?IwY{QVk5{gszJXM6`h=xCoQl=1ooFsnxzsS^6
z+-xU`##<HyXeqN_I<Dx5%4_ss;VFGuMlRwTZS*;h`R<U|v2O|pSPeV$^~OSE(}{n7
z1#|Qnk3yjY8iDyEG!&g5T-KtqF|Uyb@o5#g^Pkl7uXxOVam3I=5e)?EOtnW()fUW=
zi6|OGq&jk@9Gi0Xm!9e}DELWKCsh3e^+B4=<?RqJ%pdkq7!{^a-f^tGv0;SkWe+5&
zW=kFN2fY~4_ozUfjBzR}3k*AZa$&E6lhN*WAD~v);Vj-qL>ungw}iyHV|Z!AM#JmD
zkA*&UFsXh1Mj^gHDB3`g`PQ2$!(nZB<=0sv1&Mjr|ExV{!+6_y1x0%9zL<@(lb%b?
zYLS!5rYy)Y`t_M&<cD2?I+_rPz|8j_JV-EzNP`9F95F9*9y9<Yp%3=2y>@nXs!Aiz
z4V}7|63q-^prKOFsI&H4w{L{wo}I9JsT%-O03Lf#9FV|1bS1=^sEWzi{NHm-G_Od<
zG~Z$1@wl_e#rRhpP4hZZuC%518f0U962?(d^J~Ts(<K2p2w?)v&;@HYxDo2=>!sVJ
zq^wSLUYQYs1KATG+5ryl{nnq_m^xneQB9-jC*w2A=lEk?*C^5h^>*4Eku<S>_FBR`
zKZ66+A<d87eZI7nPlzhxQxL~!F18smK$gt1wal;G1;F)UXQ4rlN!^=^5}SZ;17O)>
zWE4j}o$6DLZlN|Ay18tq2PqVUFJG|c$Y@@f(LD6x6poa{1UibVfzz@Z3RT`!!Z%$V
zHIvl}SD0nyUnlcbOq<42j2F))KZxgc3pEb@oWhA@Upaj4UBuPtX=hkuUmA(xJ-`3@
zXX}iPOCEZDqc=Bsh)bZusW<TC(@JCY=!z<Z&flYc^M)N#zR;V=1;6Vt#Iy<i7$<5B
zC}xDm{iMDXGSiXA9TcUJ@V&Hl3FYaX6rJhIRLKrkeSg%D#ofGqu|!PkJ(>Ey9dO9x
zDFK0Ujj-d$O9xH^7a=#24NkhoOvpBL2t$=l@}NL9pUjjG&U^+o;55E!3cki5<ZDS>
zZab6S1;Plx5N&m=g2-OQ%e7zZG!#bB^yLmqqzvQ2Lm3rcd|xWayn+~|FrF4xl)Xvy
zSuZSYQx(73$s9Xz0*<?;=#3PQ*=I>j<8Iru<#%l~aLu1zZi5x*{d>ldTtnc9QF*pH
zpbvOl$an-YR6D4TM2?3EoxFrbu=PqHZZ9SU&yF}2nOpZfn~-nlcx#qPoJOF9n1CPA
zQ6fZ7Q={~}oP73iiu(#HNo!r4k_asGfL)Y_kB=4*Cs_|^M#a}AI@ULv%QnD=r%^ub
zFdCN~ziGgyVOBryCzE4RJ(-Pr4>)_L>-Q>)qHQU25nf4a?At>g55lt*i*_;tmup~T
z#0}Bt>0PYRVQhgaxJy}I%`p_U7eir7@CKjMPbdh^%m(Yhg=6~2;lmWqV#E?Sj<7$~
z&!P4#_JA&sgEg?sl#SE84|9QqIbDW!l;_moVqcT-OQ_6@fgBQR1L;U%KP4uECQP=+
zhecwr6NyI)lPZwK*`Gb_osi`~lNnPfvbXi31aTsQWh8cHz+!TDNX@_bRgh!gZ;uWN
zsk>qQpFSBOjv-}U2K3!Un`Utv5}S{N!guF(9^zYz-$M?&FG{fCU)p95Z#yugy*AvS
z=n4l4^jEN+Ee^DJiFT0*@m9c#)Yj{!H!%%Y_BnpgI(l>r)_9^V1t5z9n3jT@BOyI`
z9hdQ!Cl)NM%kVaz#&g{WWE$?0n+WMZz%$~;giiY8=E{)K)Rdcd?)Za<i-RFiG^&ck
zR4m%Bk%4)1#ODO7jNMjJMTEg~wMkySu)^kOVW)<^!~E=JGMe1Epkw?6J}&a`z-z@{
zf2=d)H_`Wv;)a(Xt@+{$K+n-(sStUNr;{j*qny~cNoX8KUl6ZhadGj53m0Cxj284q
z$Zq>xUf3tQ4}9PuwjoBvJ@i*Gl!dg4vLayv?m1!^Q&nrjIizi9C@s6LZpsYcs43Lv
zj60d4t@)W8dv}b=^ZMa!D}|^8gn1Bvm;`eR_wm;oGzKlPJp#T{!B&sA(+4v066g*e
zY;x3t0`!&)b$eOW&VCUxcp+TG%G&M5YofUX4rc^^8pcXYm+2Nr#+Gq=U6NHvDG$1E
ziZc%xR`}J<h}4k}A3kWOuz{Y~>q(o5l!~)vmgmO9p6{~O?cEG8oC)kw&`%6D1{U11
zO*+JU>K6H{D1VRkzUv9h4x=5&j(m}<87q%Y(>C&D4R%zRlhxNst+A(Bjd6?kwZDdy
zU)ZHT+VL>tm;h)<^NTzSC_Qkj5R(%29H@~&9|70eVbQWH@xfZ@Zq>Fwf7w#fx*3@k
zgxp7Cyb0-(P>vws?5_P?6W=AocY??Byj5Ml0ZM1$h(k6Sawpb<a$JFN>U+QV=>`A!
zvajoavoLXR2-&YB4J)cjrB8c*k7B`U28!-um|8R}fleALdlkTC3~>Ro`nc6_@Q$Bf
zJ_8+w87V1oo=Mjcg{_WV;R-v>Yboa+TcLP_3OkoMyP~qP1Tz-`Ae{%dbPozoFuvu`
z49ke31gxO?^Btc3a|kM7=*yL<ZaRYWlObxracB4j8{vh2``0!S<_ch(CQxNbRFx0*
z&t)G`bJRKA5htITq2kS(mxTUAR4oN1`zxT(5Ufs|YwW}>0NQyB#`WEN$=?uh&f(Zr
zbNUVnCwZwokOuI6qvW**y?XVEumA4gPiIWek<3$v7=RcIh-Vr3f)V28%mCc1e#zo~
zpjg(B#NS%Yj!{i0<Mv=L7GyqndvLq2mscD*m}U!9#$=>+{<pURFCf7uRWhWeMCBnQ
z{~mV$!|FIBqc!?#=J3LQ8Dn33atw!)tnx>Bye{bp<Na*7i}xJn_6-nI(NG9Y8Xs=>
zDt7({KS-}1r!oo42twNd$Y^3WH-?5Cy!87aQvJ`iXh;I*y<cFU{YnI@w7sy#a0Jx@
z2bSSIj=b30-(P+(4Ig9M0lI;}RMGEk+i)bavF;xK?07o#QHDc4Mw26b1-6fs<<B4`
zB|Bji=M5Z{K<HMWm9WpRywd#&`V*B0ruF0I`Ux2v=G)K#uyYOul^C~M>(!2L=jnE*
zVDJ2-`V5O8lF+rN3=eRWoF(Py{m<=Bdbe=IwGKPLX1{6cR(j+z52rRdt=7gm$1x{&
z^RQkU*X8ez=EfK28QYZ4<%sjpn3<d71$L|a9@h7Cr{nAmkwQQ4n7CENCSu)~mqB)D
ziV->>LdAj+VxO@cYT6;FlSWy!Y8BbC(1dz|a@hnb^7#A4shz!v#>30ExVh^+h@T6A
zN(Doctn_3W)Y@ai-kf%q{beiUJ2$?0n~dQ{s@#P7y1G3{_4kqPZhnRPg~GQ!^_gDv
z8R|>!pMqxgJw4oF{WLjQ3937Ua4>D$*z{>UYO){ogOu&(m3ZgCU<vO%Gr4tm#&&nd
zbA=8vYg%>Wg~CYhX6wa;#r$VON*!H<J-i$a;g8kg+TK?$?m2j1OtBylMZ#sg_<v!8
z1NmEJD0$2uSMvG4u4F8$_Z?<pTeWamk(t<=9NV}io$93*moNF2uY+HZ*BHMlD8H8^
z{#2}sgP)dU?kgTS!URfU2ZlT`i_AWyA(lO~l_WQU@)F9T9nk;#Yh70PSA`z2F|qf#
zE<QO{L39{m8oHn6jojKA(vKKT0$hpSIYg3uh$$$-JYPZv!|D(5Y2EIF7YgraD0j&D
zh^ljBBozFwfbBNQIX~|ylni8LWxe!c;0V!yf$bY;OpUDzX?Y0K)9;;mpratK))>$Q
z9RsHlwmhU+_+u7~|50Zy2TWk@?YqGta%#xw*G={E6(XV+n?$>Ndp7~z6~QX-!nXk0
zn_HKzFLAz_9gl0P&q@|vjD$fD1~W8|>m;Jozh1CB3DpmUo0yJ+dTyafUzVvS1QdWh
z8wbV#j!Fw?4R}RG4-;AqdII7<pBQNuuuqm?!_cqVNj`3FXkEZ{g~&u!0_NG|JUJqS
zWT9=u?+G(7#-!dNJ-Ccs6sS~k^$>S#!ka4qT&EOn2VVPR?dC$aJggx`3Jrx+e+lI{
znMmXX%;?#OuDT1dq6fTN=++d+)8ztmyA5NS#=pG|3a87l&p)y91JqhuF<ggZveS7`
z2>qMV>C^YT67?OifKR9EuusB4>J5$)Bvy`BVw)p8)HNiL0{2Tfd=h(Vd{FlhU}1~~
zxpqtd!R!Q4<3y5lpDl;_SWqvlrIZ45<F{WIt{y9I0GRC<9wtv|;4NsV)}u)MOG>vj
z5Y%HE%7vgJD?F@qpy6E?<|fzh=M@0i4hae2LBA0xGztaiy~aUhH@5=Nv5=GMfKmh1
z!>g>H`YPxItwX+SdT(DKhX{kspi`<8of}FdG;+SZ^H91i<Dr2VY4f4JT*RkLO;5{=
zK*7KaR10}dp`xKKDBSQE$w(Gt0;qBLUg-^#?0AX4o_?uYv0MzXJkaj#k9fP2WMGjn
zAJooARjr%R%IvWDv>#VOA!+fnD>r)JR_XSbSk!*_;6WrPfOPP=*79qVA`K=Typn%y
zdk1<-?4JSW6h0z^$bP6Lb1}XLwjm0I;54vTIgGS5*BmKWz1wj9_w$m(=X;aVTX4qX
zN(Kro%1aNbC7U&RI8kxqKqaGK)AUlvT+zA^{j4*umy>BN2rvzkx!y{8P>8LF$OhA}
zgx;OYs0(Sk1-l4FwU9H2kei}V-;le0wA|(-)-Fl(0u!vp@mdH<sS?yxq{g7U#Kuoz
z+Az}N7(SzMx{gtHE60M@u8$XXVg8(#Qlc!euDpOE;D>X8@twbm$SsFHN8@vp1ci$@
z<mcl#X<iSSj6Bf8THn~%3!T-=eGGD5SR;wOF9eab3m$GMIpa=4=yyGug0hLUTiQc;
z97XHV^pgS}XZdTKdBl^5V4^O>{@6RKDs4mI9iXzQ!)?+4*c>%9c!WVBViH0wVNX1&
zYS)ZNi6C@D@b;^I-D;@|@(fj04I)pViYq>s!2@Fwp4VSkb^tX|SCM$5psQ~XI|XFb
z4=K?8wK)kOttpZDo7S&c)8iB%>&YRafX(Dgc_kN~*uw#*RMt4)&71l8QIFE0_Xw2<
zI*I+`Fitu=E>|cmfI~xZ#Y|Fm?4L`d<Kdx?F?9QR90k7`N4W)39c*%Aky?b?LcV@6
zPGVUf9}f*m6XH&XN9A3VBm+hW7b7U+;lQ`~E>8vwq%nc~;vgOz;KE|WK#_NjWY8iR
zSLaz_`HgG}g~sfF(%xex%TV!I<C!1iNL7c1r>Hwbm2)pnSe#Jf!S1t5APY(0T^OGe
zLY?ee<He{^>CG5a{>BX-^z`Bf?Z6}DS1xy{a$l{Tw4Nlk0JtGaH{7?xO*yX{8#j=2
zTm#2;Z(txqaWk6#J*Y0d1sH=HG0Gt)9Y`#GlUEdY>c*9{I)coJn{v<M^@ft*@KCu<
zy%{O!F1?qVl$r}+K^UJ#w9J;%Vv{820NU&83j<L=_aJtQU6?&jY^&IBBq0O1hUe8A
zdcrV7#xP1fMwSJ6Cr)F<bR?Muff@`5|I!-T?HkBS$EWh3#x$qq!jJ2FZ?NNF@aGV*
z8gja9TrS|bj$?HFXxzWLLf5IZo>S`PK}ZY7K2*1Z%*6r0^!f|Ca9;cBs4Cjq6?LEL
zx5lTVO+r9E$!7`bd4^!*E|lpSl_;ATnVI>H@J_mcpafHPtS#xiWJdEvG->KjUpK><
zuRJeD@xVKVTS1YOmk7f$g}7R)=hz574+9?kP2^W$^Ps%gVX6(08PF3RYneXd^2!bW
z*(*-a9*`+C`0O)kgp~U=dMVnj$P?+oP_dw>qar+HlN*P8epyoeb`lALSA-vW@TsZZ
z;P{}=Oz4$Q3~kXtd#DWc2->|K#A7g_8-Zh06^OJNN7Y%9bgYbU9zuC|F)9LPT8xAO
zV6*@3J8_3R1K5=O&d0iGf8ENt<Mb;UvTY*^gR!S<67{trq0!gyo9A)4&mMt3QbQ+2
zn<R?uKxRn|bZ|8vK5*MVZw3CpHbEnL4f6CiP)Y$FmSPXXQR<D|pM=5;*vXLHf5;)v
z39X<q*68}yeKt4=n6bm)G5+Ty*cdl6H6<l1xX{`P{Te>vE!;ga6%E=K3O8<O=}D`C
zE%t)JM{F*m+@~T@P*i(J$XrF{n$Q6<XAOXpVp2Vsvj+H}{{>&dM|P=8IJ$9pw`8uX
zVibpA0H>TD79>I6aEg{9=X*PvjE6xD=P+k920D_U&REbo0-t5NPXUw6kb_z7Jrb`_
zb~0j@dIbkHqA3heYWZQOI(WBbcUo_#_uspUXPb=S7Oz(d?w{f@eP`_Zdxd?dHP63~
zxqD&13Xz&ZOG0L~I<{H+)>+Fw0;mNYH)I$Lsm22P!bMIYwPdNC3<rwm_@Lfr*n;%)
zooWWJxU+^LhZ7wq5gFU$u-sQ!+{G}Pcl_;#BOLFk7Q11<lPkxzoWwFDzgW5`m<mmP
u?Wg_SFyRkxQT^#?F(|_0PfLj}om%MC@DO-kB&I>ZuLH^&N>3GyuKj<CcuQdb

literal 0
HcmV?d00001

diff --git a/lib/matplotlib/tests/test_axes.py b/lib/matplotlib/tests/test_axes.py
index 40c23f6c5967..1a848692907e 100644
--- a/lib/matplotlib/tests/test_axes.py
+++ b/lib/matplotlib/tests/test_axes.py
@@ -6073,3 +6073,36 @@ def test_deprecated_uppercase_colors():
     with pytest.warns(MatplotlibDeprecationWarning):
         ax.plot([1, 2], color="B")
         fig.canvas.draw()
+
+
+@image_comparison(baseline_images=['secondary_xy'], style='mpl20',
+        extensions=['png'])
+def test_secondary_xy():
+    fig, axs = plt.subplots(1, 2, figsize=(10, 5), constrained_layout=True)
+
+    def invert(x):
+        with np.errstate(divide='ignore'):
+            return 1 / x
+
+    for nn, ax in enumerate(axs):
+        ax.plot(np.arange(2, 11), np.arange(2, 11))
+        if nn == 0:
+            secax = ax.secondary_xaxis
+        else:
+            secax = ax.secondary_yaxis
+
+        axsec = secax(0.2, functions=(invert, invert))
+        axsec = secax(0.4, functions=(lambda x: 2 * x, lambda x: x / 2))
+        axsec = secax(0.6, functions=(lambda x: x**2, lambda x: x**(1/2)))
+        axsec = secax(0.8)
+
+
+def test_secondary_fail():
+    fig, ax = plt.subplots()
+    ax.plot(np.arange(2, 11), np.arange(2, 11))
+    with pytest.raises(ValueError):
+        axsec = ax.secondary_xaxis(0.2, functions=(lambda x: 1 / x))
+    with pytest.raises(ValueError):
+        axsec = ax.secondary_xaxis('right')
+    with pytest.raises(ValueError):
+        axsec = ax.secondary_yaxis('bottom')
diff --git a/lib/mpl_toolkits/mplot3d/axes3d.py b/lib/mpl_toolkits/mplot3d/axes3d.py
index 922972cd423c..fc4f83ecf697 100644
--- a/lib/mpl_toolkits/mplot3d/axes3d.py
+++ b/lib/mpl_toolkits/mplot3d/axes3d.py
@@ -801,8 +801,6 @@ def get_zscale(self):
         """
         Return the zaxis scale string %s
 
-        .. versionadded :: 1.1.0
-            This function was added, but not tested. Please report any bugs.
         """ % (", ".join(mscale.get_scale_names()))
         return self.zaxis.get_scale()
 
@@ -829,24 +827,34 @@ def set_yscale(self, value, **kwargs):
             This function was added, but not tested. Please report any bugs.
         """
 
-    @docstring.dedent_interpd
     def set_zscale(self, value, **kwargs):
         """
-        Set the scaling of the z-axis: %(scale)s
+        Set the x-axis scale.
 
-        ACCEPTS: [%(scale)s]
-
-        Different kwargs are accepted, depending on the scale:
+        Parameters
+        ----------
+        value : {"linear", "log", "symlog", "logit", ...}
+            The axis scale type to apply.
 
-        %(scale_docs)s
+        **kwargs
+            Different keyword arguments are accepted, depending on the scale.
+            See the respective class keyword arguments:
 
-        .. note ::
-            Currently, Axes3D objects only supports linear scales.
-            Other scales may or may not work, and support for these
-            is improving with each release.
+            - `matplotlib.scale.LinearScale`
+            - `matplotlib.scale.LogScale`
+            - `matplotlib.scale.SymmetricalLogScale`
+            - `matplotlib.scale.LogitScale`
 
-        .. versionadded :: 1.1.0
-            This function was added, but not tested. Please report any bugs.
+        Notes
+        -----
+        Currently, Axes3D objects only supports linear scales.
+        Other scales may or may not work, and support for these
+        is improving with each release.
+
+        By default, Matplotlib supports the above mentioned scales.
+        Additionally, custom scales may be registered using
+        `matplotlib.scale.register_scale`. These scales may then also
+        be used here as support is added.
         """
         self.zaxis._set_scale(value, **kwargs)
         self.autoscale_view(scalex=False, scaley=False)

From 67ef64ee6b6f7855f8373f063b28c3fdabffec4b Mon Sep 17 00:00:00 2001
From: Jody Klymak <jklymak@gmail.com>
Date: Tue, 5 Feb 2019 07:39:40 -0800
Subject: [PATCH 2/3] FIX: more fixes

---
 .../secondary_axis.py                         | 23 ++++---------------
 lib/matplotlib/scale.py                       |  6 ++---
 lib/mpl_toolkits/mplot3d/axes3d.py            |  2 +-
 3 files changed, 7 insertions(+), 24 deletions(-)

diff --git a/examples/subplots_axes_and_figures/secondary_axis.py b/examples/subplots_axes_and_figures/secondary_axis.py
index 466b1a2f8e32..caa27a9c895f 100644
--- a/examples/subplots_axes_and_figures/secondary_axis.py
+++ b/examples/subplots_axes_and_figures/secondary_axis.py
@@ -6,9 +6,9 @@
 Sometimes we want as secondary axis on a plot, for instance to convert
 radians to degrees on the same plot.  We can do this by making a child
 axes with only one axis visible via `.Axes.axes.secondary_xaxis` and
-`.Axes.axes.secondary_yaxis`.
-
-If we want to label the top of the axes:
+`.Axes.axes.secondary_yaxis`.  This secondary axis can have a different scale
+than the main axis by providing both a forward and an inverse conversion
+function in a tuple to the ``functions`` kwarg:
 """
 
 import matplotlib.pyplot as plt
@@ -26,22 +26,6 @@
 ax.set_xlabel('angle [degrees]')
 ax.set_ylabel('signal')
 ax.set_title('Sine wave')
-secax = ax.secondary_xaxis('top')
-plt.show()
-
-###########################################################################
-# However, its often useful to label the secondary axis with something
-# other than the labels in the main axis. In that case we need to provide
-# both a forward and an inverse conversion function in a tuple
-# to the ``functions`` kwarg:
-
-fig, ax = plt.subplots(constrained_layout=True)
-x = np.arange(0, 360, 1)
-y = np.sin(2 * x * np.pi / 180)
-ax.plot(x, y)
-ax.set_xlabel('angle [degrees]')
-ax.set_ylabel('signal')
-ax.set_title('Sine wave')
 
 
 def deg2rad(x):
@@ -98,6 +82,7 @@ def inverse(x):
 
 xold = np.arange(0, 11, 0.2)
 # fake data set relating x co-ordinate to another data-derived co-ordinate.
+# xnew must be monotonic, so we sort...
 xnew = np.sort(10 * np.exp(-xold / 4) + np.random.randn(len(xold)) / 3)
 
 ax.plot(xold[3:], xnew[3:], label='Transform data')
diff --git a/lib/matplotlib/scale.py b/lib/matplotlib/scale.py
index fb55d33156a0..d70049ec3b8b 100644
--- a/lib/matplotlib/scale.py
+++ b/lib/matplotlib/scale.py
@@ -163,8 +163,7 @@ def __init__(self, axis, functions):
 
         functions : (callable, callable)
             two-tuple of the forward and inverse functions for the scale.
-            The forward function must have an inverse and, for best behavior,
-            be monotonic.
+            The forward function must be monotonic
 
             Both functions must have the signature::
 
@@ -426,8 +425,7 @@ def __init__(self, axis, functions, base=10):
 
         functions : (callable, callable)
             two-tuple of the forward and inverse functions for the scale.
-            The forward function must have an inverse and, for best behavior,
-            be monotonic.
+            The forward function must be monotonic.
 
             Both functions must have the signature::
 
diff --git a/lib/mpl_toolkits/mplot3d/axes3d.py b/lib/mpl_toolkits/mplot3d/axes3d.py
index fc4f83ecf697..fe644d9cb967 100644
--- a/lib/mpl_toolkits/mplot3d/axes3d.py
+++ b/lib/mpl_toolkits/mplot3d/axes3d.py
@@ -829,7 +829,7 @@ def set_yscale(self, value, **kwargs):
 
     def set_zscale(self, value, **kwargs):
         """
-        Set the x-axis scale.
+        Set the z-axis scale.
 
         Parameters
         ----------

From 85c1b99efb7eb22221067aebd6fe13e7e8b88837 Mon Sep 17 00:00:00 2001
From: Jody Klymak <jklymak@gmail.com>
Date: Fri, 8 Feb 2019 07:07:15 -0800
Subject: [PATCH 3/3] FIX

---
 lib/matplotlib/axes/_axes.py           | 3 +--
 lib/matplotlib/axes/_secondary_axes.py | 2 +-
 lib/matplotlib/scale.py                | 4 ++--
 3 files changed, 4 insertions(+), 5 deletions(-)

diff --git a/lib/matplotlib/axes/_axes.py b/lib/matplotlib/axes/_axes.py
index 13e8cda4e622..a6238e0c6ae0 100644
--- a/lib/matplotlib/axes/_axes.py
+++ b/lib/matplotlib/axes/_axes.py
@@ -613,8 +613,7 @@ def secondary_xaxis(self, location, *, functions=None, **kwargs):
         Examples
         --------
 
-        Add a secondary axes that shows both wavelength for the main
-        axes that shows wavenumber.
+        The main axis shows frequency, and the secondary axis shows period.
 
         .. plot::
 
diff --git a/lib/matplotlib/axes/_secondary_axes.py b/lib/matplotlib/axes/_secondary_axes.py
index 83960601c260..3cbe8e0f1ff1 100644
--- a/lib/matplotlib/axes/_secondary_axes.py
+++ b/lib/matplotlib/axes/_secondary_axes.py
@@ -441,7 +441,7 @@ def set_color(self, color):
 
     If a 2-tuple of functions, the user specifies the transform
     function and its inverse.  i.e.
-    `functions=(lamda x: 2 / x, lambda x: 2 / x)` would be an
+    `functions=(lambda x: 2 / x, lambda x: 2 / x)` would be an
     reciprocal transform with a factor of 2.
 
     The user can also directly supply a subclass of
diff --git a/lib/matplotlib/scale.py b/lib/matplotlib/scale.py
index d70049ec3b8b..8ce9b68a7165 100644
--- a/lib/matplotlib/scale.py
+++ b/lib/matplotlib/scale.py
@@ -163,7 +163,7 @@ def __init__(self, axis, functions):
 
         functions : (callable, callable)
             two-tuple of the forward and inverse functions for the scale.
-            The forward function must be monotonic
+            The forward function must be monotonic.
 
             Both functions must have the signature::
 
@@ -411,7 +411,7 @@ def limit_range_for_scale(self, vmin, vmax, minpos):
 class FuncScaleLog(LogScale):
     """
     Provide an arbitrary scale with user-supplied function for the axis and
-    then put on a logarithmic axes
+    then put on a logarithmic axes.
     """
 
     name = 'functionlog'