pmagpy 4.3.8__py3-none-any.whl → 4.3.10__py3-none-any.whl

pmagpy/rockmag.py

@@ -22,10 +22,9 @@ except ImportError:
 try:
     from bokeh.plotting import figure, show
     from bokeh.layouts import gridplot
-    from bokeh.models import HoverTool, Label
+    from bokeh.models import HoverTool, ColumnDataSource, PointDrawTool, CustomJS, Div
     from bokeh.embed import components
     from bokeh.palettes import Category10
-    from bokeh.models import ColumnDataSource
     from bokeh.models.widgets import DataTable, TableColumn
     from bokeh.layouts import column
     _HAS_BOKEH = True
@@ -3384,18 +3383,22 @@ def process_hyst_loop(field, magnetization, specimen_name, show_results_table=Tr
     sigma = np.log(E_hyst / 2 / Bc / Ms)
 
     # plot original loop
-    p = plot_hysteresis_loop(grid_fields, grid_magnetizations, specimen_name, line_color='orange', label='raw loop')
+    p = plot_hysteresis_loop(grid_fields, grid_magnetizations, specimen_name, line_color='orange', label='raw loop', 
+                             return_figure=True, show_plot=False)
     # plot centered loop
-    p_centered = plot_hysteresis_loop(centered_H, centered_M, specimen_name, p=p, line_color='red', label=specimen_name+' offset corrected')
+    p_centered = plot_hysteresis_loop(centered_H, centered_M, specimen_name, p=p, line_color='red', label=specimen_name+' offset corrected', 
+                                      return_figure=True, show_plot=False)
     # plot drift corrected loop
-    p_drift_corr = plot_hysteresis_loop(centered_H, drift_corr_M, specimen_name, p=p_centered, line_color='pink', label=specimen_name+' drift corrected')
+    p_drift_corr = plot_hysteresis_loop(centered_H, drift_corr_M, specimen_name, p=p_centered, line_color='pink', label=specimen_name+' drift corrected', 
+                                        return_figure=True, show_plot=False)
     # plot slope corrected loop
-    p_slope_corr = plot_hysteresis_loop(centered_H, slope_corr_M, specimen_name, p=p_drift_corr, line_color='blue', label=specimen_name+' slope corrected')
+    p_slope_corr = plot_hysteresis_loop(centered_H, slope_corr_M, specimen_name, p=p_drift_corr, line_color='blue', label=specimen_name+' slope corrected', 
+                                        return_figure=True, show_plot=False)
     # plot Mrh
     p_slope_corr.line(H, Mrh, line_color='green', legend_label='Mrh', line_width=1)
     p_slope_corr.line(H, Mih, line_color='purple', legend_label='Mih', line_width=1)
     p_slope_corr.line(H, Me, line_color='brown', legend_label='Me', line_width=1)
-    show(p)
+    show(p_slope_corr)
     results = {'gridded_H': grid_fields, 
                'gridded_M': grid_magnetizations, 
                'linearity_test_results': loop_linearity_test_results,
@@ -3627,7 +3630,7 @@ def plot_X_T(
         temperature_column (str): Name of temperature column.
         magnetic_column (str): Name of susceptibility column.
         temp_unit (str): "C" for Celsius.
-        smooth_window (int): Window for smoothing.
+        smooth_window (int): Window for smoothing, if 0, no smoothing is applied.
         remove_holder (bool): Subtract holder signal.
         plot_derivative (bool): Plot derivative.
         plot_inverse (bool): Plot inverse.
@@ -3861,6 +3864,174 @@ def plot_X_T(
         return tuple(figs)
     return None
 
+def estimate_curie_temperature(
+    experiment,
+    temperature_column="meas_temp",
+    magnetic_column="susc_chi_mass",
+    temp_unit="C",
+    smooth_window=0,
+    remove_holder=True,
+    figsize=(6, 6),
+    inverse_method=False,
+    print_estimates=True
+):
+    """
+    Estimate the Curie temperature from high temperature susceptibility curves in multiple ways. Automatically calculates first derivative minimum, second derivative maximum, and
+    second derivative zero-crossing. Also has option to estimate Curie temperature from inverse susceptibility (E. Petrovsky and A. Kapicka, 2006). Uses Bokeh for the interactive plot.
+
+    Parameters:
+        experiment (pandas.DataFrame): MagIC-formatted experiment DataFrame.
+        temperature_column (str): Name of temperature column.
+        magnetic_column (str): Name of susceptibility column.
+        temp_unit (str): "C" for Celsius.
+        smooth_window (int): Window for smoothing, if 0, no smoothing is applied.
+        remove_holder (bool): Subtracts holder signal to better isolate specimen signal.
+        figsize (tuple): (width, height) in inches.
+        inverse_method (bool): Use inverse susceptibility method.
+        print_estimates (bool): Print estimated Curie temperatures.
+    
+    Returns
+    -------
+
+    temp_of_first_derivative_min_heating : float
+        temperature of first derivative minimum for heating cycle
+    temp_of_first_derivative_min_cooling : float
+        temperature of first derivative minimum for cooling cycle
+    temp_of_second_derivative_max_heating : float
+        temperature of second derivative maximum for heating cycle
+    temp_of_second_derivative_max_cooling : float
+        temperature of second derivative maximum for cooling cycle
+    heating_zero : list of float or None
+        temperatures of second derivative zero-crossings for heating cycle, or None if none found
+    cooling_zero : list of float or None
+        temperatures of second derivative zero-crossings for cooling cycle, or None if none found
+    """
+
+    warm_T, warm_X, cool_T, cool_X = split_warm_cool(
+        experiment,
+        temperature_column=temperature_column,
+        magnetic_column=magnetic_column
+    )
+
+    if temp_unit == "C":
+        warm_T = [T - 273.15 for T in warm_T]
+        cool_T = [T - 273.15 for T in cool_T]
+    else:
+        raise ValueError('temp_unit must be "C"')
+
+    if remove_holder:
+        holder_w = min(warm_X)
+        holder_c = min(cool_X)
+        warm_X = [X - holder_w for X in warm_X]
+        cool_X = [X - holder_c for X in cool_X]
+
+    swT, swX = smooth_moving_avg(warm_T, warm_X, smooth_window)
+    scT, scX = smooth_moving_avg(cool_T, cool_X, smooth_window)
+
+    dx_w = np.gradient(swX, swT)
+    dx_c = np.gradient(scX, scT)
+
+    temp_of_first_derivative_min_heating = swT[np.argmin(dx_w)]
+    temp_of_first_derivative_min_cooling = scT[np.argmin(dx_c)]
+
+    dx_2_w = np.gradient(dx_w, swT)
+    dx_2_c = np.gradient(dx_c, scT)
+
+    temp_of_second_derivative_max_heating = swT[np.argmax(dx_2_w)]
+    temp_of_second_derivative_max_cooling = scT[np.argmax(dx_2_c)]
+
+    # Physically consistent zero-crossing logic for heating
+    min_idx_w = np.argmin(dx_2_w)
+    max_idx_w = np.argmax(dx_2_w)
+    lower_w, upper_w = sorted([min_idx_w, max_idx_w])
+    dx2_slice_w = dx_2_w[lower_w:upper_w]
+    temp_slice_w = swT[lower_w:upper_w]
+    crossings_w = np.where(np.diff(np.sign(dx2_slice_w)))[0]
+    temp_of_zero_crossing_heating = [temp_slice_w[i] for i in crossings_w]
+
+    # Physically consistent zero-crossing logic for cooling
+    min_idx_c = np.argmin(dx_2_c)
+    max_idx_c = np.argmax(dx_2_c)
+    lower_c, upper_c = sorted([min_idx_c, max_idx_c])
+    dx2_slice_c = dx_2_c[lower_c:upper_c]
+    temp_slice_c = scT[lower_c:upper_c]
+    crossings_c = np.where(np.diff(np.sign(dx2_slice_c)))[0]
+    temp_of_zero_crossing_cooling = [temp_slice_c[i] for i in crossings_c]
+
+    if inverse_method:
+        if not _HAS_BOKEH:
+            raise ImportError("Bokeh is required for inverse_method=True")
+        bokeh_height = int(figsize[1] * 96)
+        title = experiment["specimen"].unique()[0]
+        swX_arr = np.array(swX)
+        inv_w = np.divide(1.0, swX_arr, out=np.full_like(swX_arr, np.nan), where=swX_arr != 0.0)
+        mask_w = np.isfinite(inv_w)
+        T = np.array(swT)[mask_w]
+        inv_chi = inv_w[mask_w]
+
+        # Creates initial fit endpoints (choose two points in the linear region)
+        # 0.7 is Magic Number which places initial fit near expected linear region
+        fit_x = [T[int(len(T)*0.7)], T[-1]]
+        fit_y = [inv_chi[int(len(inv_chi)*0.7)], inv_chi[-1]]
+        fit_source = ColumnDataSource(data=dict(x=fit_x, y=fit_y))
+        # scatter and line for data
+        data_source = ColumnDataSource(data=dict(x=T, y=inv_chi))
+        p_inv = figure(title=f"{title} – 1/χ",
+                       height=bokeh_height,
+                       x_axis_label=f"Temperature (°{temp_unit})",
+                       y_axis_label="1/χ",
+                       tools="pan,wheel_zoom,box_zoom,reset,save",
+                      )
+        p_inv.scatter('x', 'y', source=data_source, size=8, color="red", legend_label="Heating – 1/χ")
+        renderer = p_inv.scatter('x', 'y', source=fit_source, size=12, color="blue", legend_label="Fit Endpoints")
+        p_inv.line('x', 'y', source=fit_source, line_width=2, color="blue", legend_label="Fit Line")
+        # PointDrawTool for dragging endpoints
+        draw_tool = PointDrawTool(renderers=[renderer], add=False)
+        p_inv.add_tools(draw_tool)
+        p_inv.toolbar.active_tap = draw_tool
+        p_inv.legend.location = "top_left"   
+        # Div to display Curie temperature
+        curie_estimate = Div(text="Curie temperature: --", styles={'font-size': '16px', 'color': 'darkred'})  
+        # JS callback to update fit line and Curie temperature estimate
+        callback = CustomJS(args=dict(source=fit_source, div=curie_estimate), code="""
+            var x = source.data.x;
+            var y = source.data.y;
+            if (x.length == 2) {
+                var slope = (y[1] - y[0]) / (x[1] - x[0]);
+                var intercept = y[0] - slope * x[0];
+                // Estimate Curie temperature: x where y=0
+                var Tc = -intercept / slope;
+                div.text = "Curie temperature estimate: " + Tc.toFixed(2) + " °C";
+            }
+        """)
+        fit_source.js_on_change('data', callback)       
+        show(column(p_inv, curie_estimate))
+
+    if print_estimates:
+        print(f'First derivative minimum is at T={int(temp_of_first_derivative_min_heating)} for heating')
+        print(f'First derivative minimum is at T={int(temp_of_first_derivative_min_cooling)} for cooling')
+        print(f'Second derivative maximum is at T={int(temp_of_second_derivative_max_heating)} for heating')
+        print(f'Second derivative maximum is at T={int(temp_of_second_derivative_max_cooling)} for cooling')
+        if temp_of_zero_crossing_heating:  
+            print(f'The second derivative of the heating curve crosses zero at T = {int(temp_of_zero_crossing_heating[0])}')  
+        else:  
+            print('No zero crossing found for the second derivative of the heating curve.')  
+        if temp_of_zero_crossing_cooling:  
+            print(f'The second derivative of the cooling curve crosses zero at T = {int(temp_of_zero_crossing_cooling[0])}')  
+        else:
+            print('No zero crossing found for the second derivative of the cooling curve.')
+
+    heating_zero = temp_of_zero_crossing_heating[0] if temp_of_zero_crossing_heating else None
+    cooling_zero = temp_of_zero_crossing_cooling[0] if temp_of_zero_crossing_cooling else None
+
+    return (
+        temp_of_first_derivative_min_heating,
+        temp_of_first_derivative_min_cooling,
+        temp_of_second_derivative_max_heating,
+        temp_of_second_derivative_max_cooling,
+        heating_zero,
+        cooling_zero
+    )
 
 def smooth_moving_avg(
     x,

pmagpy/spline.py

@@ -89,7 +89,7 @@ class Spline(func.FuncOps):
         """
         Simulate a ufunc; handle being called on an array.
         """
-        if type(arg) == func.ArrayType:
+        if isinstance(arg, (numpy.ndarray, list, tuple)):
             return func.array_map(self.call, arg)
         else:
             return self.call(arg)
@@ -127,7 +127,7 @@ class LinInt(func.FuncOps):
         """
         Simulate a ufunc; handle being called on an array.
         """
-        if type(arg) == func.ArrayType:
+        if isinstance(arg, (numpy.ndarray, list, tuple)):
             return func.array_map(self.call, arg)
         else:
             return self.call(arg)
@@ -165,8 +165,8 @@ def logspline_interpolate(x1, y1, x2):
     Given a function at a set of points (x1, y1), interpolate to
     evaluate it at points x2.
     """
-    sp = Spline(log(x1), log(y1))
-    return exp(sp(log(x2)))
+    sp = Spline(numpy.log(x1), numpy.log(y1))
+    return numpy.exp(sp(numpy.log(x2)))
 
 
 def linear_interpolate(x1, y1, x2):

pmagpy/version.py

@@ -3,5 +3,5 @@ Module contains current pmagpy version number.
 Version number is displayed by GUIs
 and used by setuptools to assign number to pmagpy/pmagpy-cli.
 """
-"pmagpy-4.3.8"
-version = 'pmagpy-4.3.8'
+"pmagpy-4.3.10"
+version = 'pmagpy-4.3.10'

{pmagpy-4.3.8.data → pmagpy-4.3.10.data}/data/data_files/PmagPy_plots_analysis.ipynb

@@ -1162,7 +1162,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 18,
+   "execution_count": 3,
    "metadata": {},
    "outputs": [
     {
@@ -1293,7 +1293,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAASUAAAElCAYAAACiZ/R3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztnXd4k1X7x78nadrstGwUEAWEskGcgALK8EVAXMhQXxUZIiiviAP9iYIighNFBERAhhNkKRsEAUFkFsrevIwCbWaz798ft30LspI0yZO053NduZ6O53nOnTb55j7n3EMQESQSiSRRUCltgEQikZyPFCWJRJJQSFGSSCQJhRQliUSSUEhRkkgkCYUUJYlEklBIUZJIJAmFFCWJRJJQSFGSSCQJRYrSBkgSAyGECkAVADcAsAAwATADMGk0mnStVls6JSWllEqlsgCwEJEpEAgY/H6/wefz6bxeb6pGo/FrNJr8lJQUl1qtdqpUKhsAGxHl+f3+PI/Hc9br9eYSkR2ADUDB8QiA/UTkU+TJSxIKIdNMShZCiDQANQBkqtXqTIvF0iQYDNZxOp2VjUaj74YbbvCVKlVKpKenq9LT01MyMjI0FotFbTKZYDabcamjyWSCwWCA2+2G3W6HzWa74Hj+11arlfLy8ry5ubn+vLy8QF5eHh06dEh95swZrcFgOJmSkrLbbrdv9Hq9OwBkA9j9t4hJSghSlIopQoh0AJkAaul0uvp6vb6xz+er6XK5ylSoUMFVu3ZtatSokaF27drqzMxM1KxZE2azWTF73W439u7di+zsbGRnZ9PmzZudWVlZwSNHjuhTU1PtWq12n8fj2exwOLYA2AUWrFMkX8DFDilKxQQhhAHAXQaD4V8pKSn3ud3ua66//vr8unXrqhs1amTIzMxErVq1UK1aNaSmpiptbsgEg0EcPnwYu3btQnZ2NrZt25a/detW7969e7VCCLtarV5stVrnAlhORDlK2yspOlKUkhQhRAqAmzQaTVuTydTZ4XDUqV+/fn6nTp2MrVu3Vt10001ISSm+S4ZEhF27dmHJkiWYM2eOfe3atWlarfaox+OZl5+f/wuA34koX2k7JeEjRSlJEEIIANWFEK0zMjI6O53Optdcc42/Q4cO2nbt2mmaN28Oo9GotJmK4fP5sH79eixevDg4Z84cx65du7RGo3GrzWab7ff7FwPYTERBpe2UXB0pSgmMEEIHoL3ZbO4cDAbbajQaXbt27dC+fXv93XffjQoVKihtYsJis9mwcuVK/Prrr94FCxZ4cnJyVDqdblVeXt4sIppFROeUtlFyaaQoJRh/e0SNjUZjX7/f37Vx48aBLl26mFq3bo1atWqBfy0Jl6NHj2LZsmWYNWuWc/HixSlarXaJ1WodA2AZEQWUtk9SiBSlBEEIUVqtVvcwGo0D0tLSyvft21f75JNPqq+77jqlTSt2nDt3DjNnzsSYMWPsx44d8wQCgfFut3siER1U2jaJFCVFEUKoAdyTnp7e3+1239O+fftA37599S1btoRKJYPt48GWLVswfvx4z9SpUyk1NXVbXl7eJ0Q0Wy6SK4cUJQUQQtyg1Wp7qtXqXpUrV07t37+/qWvXrsjIyFDatBKL2+3G3LlzMWbMGPvGjRvVKSkp3zocjrEANslYqPgiRSlOCCFShRCPWCyW5/1+f90nnnhC9cwzz6Q2aNBAadMk/+DIkSP4+uuvA2PHjnV7PJ7TDofj00AgMJmI8pS2rSQgRSnGCCE0arX6Ca1W+279+vV1AwcONHbs2BFpaWlKmya5CsFgECtXrsTnn3/u+vXXX4mIRrvd7g+JyKa0bcUZKUoxQgiRIoToZjAYRtarV884atQoY9OmTZU2SxIhe/fuxeuvv54/b968QDAYHOHxeD4lIofSdhVHpChFGSGEWgjxiNFofL9mzZrpo0aNMrZo0UJpsyRRIjs7G6+99ppr0aJF/kAgMMzr9Y4lIpfSdhUnpChFCSGESqVSPWAwGEZff/31pT/44APj3XffLeOKiinbt2/HK6+84ly5cqXP5/O95fP5xhGRW2m7igNSlIqIYDXqaDQaP6hUqVL50aNHG9u1ayfFqISwefNmvPzyy861a9d6PB7PG36/fyIReZW2K5mRohQhf0de32s2mz8sX758pdGjRxs6dOggxaiE8ueff+Kll15ybty40eV2u4f8vVsni9ZFgBSlCBBC1DWbzdNKly5dfdSoUYbOnTvLYEcJAGDt2rUYPHiwY+vWrXaHw/EkES1S2qZkQ4pSGAghUnU63f+pVKr/jB49WturVy8hxUhyKRYtWoTHH3/c5Xa759tstr4yATh0pCiFiBDiZpPJ9N1tt91WftKkSfpKlSopbZIkwXE4HHj55Zc9kydPzne73T0DgcBPStuUDEhRugpCCL3RaBypVquf/uKLL3SPPvqoXDeShMWaNWvQvXt3Z25u7mqbzfYUEZ1Q2qZERs49roAQoqXRaNzXrl27p/fu3avr2rWrFCRJ2DRt2hS7du0yPPvss3frdLo9Go3mKSFfSJdFekqXQAhhMZlMn6ampj48efJk3X333ae0SSUDnw9YsgS45hqgYcPQrnE4gMOHAY0GUKsBjwcIBACLhe+TYCWBt2zZgq5duzr/+9//brXZbN2J6JDSNiUaUpT+gVqt7qDT6SZ36dJF/+GHH2otFovSJpUcevcGxo8HVCrg2WeBVq345xkZQP36fMzOBr77DrDZgPnzgX37Ln8/gwFo0gTwegGdDqheHahWDahRg4XL4QBuugmoVy8+z+9v/H4/Ro0aFRg+fLjH7/cP8Xq9Y2ShuUKkKP2NEKKs2WyeYDKZWk+bNk0vU0PizMaNwO23A927s8fzzTfsOZ1PSgoQDAJEfE6DBkCnToUiEwyyx6RSAVYrsHkzsHUrkJYGuFzA/v3AmTMXj92uHT+efz4+z/Vvdu/ejR49ejj37Nmz32azPUpE2XE1IEGRogRACNFUr9fP69Wrl+Gdd95J1ev1SptUcti1C1iwABgxgr2ZbdvYI3I42CtKSQFycoDt24GzZ/n7/v2BMmWASJZl8vKAvXtZ1IxGYPhwFkCAxS7OBINBjBs3jgYPHpzv9Xqf9Xq9U+JuRKJBRCX6odFoehuNRteCBQtIEkdWriSqV4+IpYCocWOiffuUsaV1a7Zh+HBlxieirKwsuvbaa51Go/FzACmUAO8NpR4l1lMSQqSaTKZxpUqV6rJ48WL9jTfeqLRJJYfVq4EWLYAbbgAGDADuvx+oXFk5e5xOXhS32YCdO4HMTEXMOHfuHDp37uzasmXLVpvNdh+V1IBLpVVRiQeAcmazeWObNm2ceXl5JIkjx48TVatGlJZGdOaM0tYwb7zBnlLDhkQ2m6Km+Hw+euGFFzwGg+EEgLqUAO+XeD9KXJySEKKxwWDI6tevX/1ff/1VL3fX4sz77/OC8/z5QOnSSlvDBP/uUXn6NK9tOZ2KmZKSkoKPPvooddy4ceX1ev16tVr9gGLGKIXSqhjPh1qt7mYwGJw//PADSRRiyBAitZpozx6lLSnE7yeaMYPo9tvZY7JYiAYNIvr5Z6KjRxUza+PGjVSmTBmXwWB4F4CKEuA9FI+H4gbE5UkCaqPR+FGFChWcW7duJYmCDBvGL7tZs5S25GKCQaLFi4keeohIpaL/LcI3aEA0ZgyR1Rp3k06ePElNmjRxms3mxQBMlADvp1g/FDcg5k8QyDCbzb/dcccdzpycnIv+6ZI4EAwS/fYb0b//TSQEUatWRG630lZdGZuNaN06olGjiG65hf7nQY0cyZ5VHPF4PPT000+7jUbjIQDVKQHeV7F8KG5ATJ8ckGk0Go/169fP4/P5LvHvlsQMn49o/nyiCROI/vUvfqkZjUT9+hE5nUpbFz7r1xO1b8/P47XXiAKBy5/rdBKdPRt1E7744ouAXq+3qVSqtpQA769YPRQ3IGZPDGii1+ttkyZNCl7yPyyJHcEg0eOP0/+mP0Yj0YcfJqcYnU8wSPTYY/yc6tUjmjyZ6PTpwt9brUTt2hU+7+bNiaZMiapntXr1arJYLC6NRtONEuB9FouH4gbE5En9LUhz5sy5zL9WEjPsdqJnnuGX1osvEh0+TJSfr7RV0aNgUbxGDX6OQhB170504gTRvHn0v9CCoUOJbryRv3/00agumGdlZVF6enqxFSbFDYj6E5KCpBx79xLdcAO/UQcPZs+iuBIIEG3cSPTSS0SpqUQmE9HHHxOVLk1UpgyRx8PPf8QIfptpNER9+hBt2RKV4YuzMCluQFSfjBQk5XA42DMoXZpo1SqlrYkve/YQ3XMP/W/aVrr0hd7h7t1EvXpxwCjA4rR3b5GHLa7CpLgBUXsiUpCU5c03+eW0bJnSlihDIEC0YAHRe+9dXnDOniV69lmilBT+W7VuTfTjj1deNL8KxVGYFDcgKk9CCpKyLFvGU5iHH1bakuTg+HGit94iqlyZ/hcH9csvEU93i5swJX1CrhCiiV6vXz5z5kxTx44dlTan5LFgAdClC3D99cCqVVx2RBIagQDw/ffA668DBw5wIbuUFH507gw880zIqTg7duxAs2bN8p1OZ0+v1zsjxpbHlKQWJSlICnPqFFCrFgvSggVAxYpKW5SceL3AV18Bn33GeXhuN3DoEP9u0yagUaOQblNshElpVy3SB+SUTVmCQaKnn+b1kexspa0pfjz5JE/tKlQg8npDvqw4TOUUNyAio6UgKc/IkfzyGTRIaUuKL1On8t+4d++wLkt2YVLcgLANBqrp9XqrFCQFOX6cSKcj6ty5SDtHkhB4/nl+m65fH9ZlWVlZZLFYXMmYkpJU9ZSEECaj0bj0/fffN8o1JAVZtQrIzwfuvZeL9Etix9tvA9deC/z731y3PETq1KmDefPm6bRa7Q9CiOqxMzD6JM0rSgihMpvNPz700EMVn3322aSxu1hSvz4X7n/pJWDLFqWtKd6YzcCUKcDu3UC3bhd3eLkCzZs3x+jRow1Go3GpEMIUQyuji9KuWqgPvV4/vHHjxk6PxxOOFyuJFQcO8LSia1elLSkZfPYZ/707dw5r4ZuI6KmnnnL/XY8pKQrFJYXHoVar79fpdP9ZsGCBPjU1VWlzJABw8CAflSz4X5Lo1w/4+GNg9mwOHQiDL774Iq1GjRp3GAyG4TGyLqokvCgJIepotdppv/76q65ChQpKmyMpQKPh4+HDytpRkhgwAGjbFnj5ZWDlypAvS01NxYIFCww6ne6FZKj5ndDBk0KIUgaDIWvs2LEVHn/88Qg6D0piylNPAT/+yF1npQcbH/LygNtu442GHTu4oWaIbNy4ES1atHA6nc7biCgrhlYWiYT1lIQQKWazee7TTz9dWgpSglK9OmC3h7X4Kiki6enApEnAkSPAc8+F1dW3SZMmGDdunN5gMCwWQpSKoZVFImFFyWAwfNigQYNGH3zwgfwITlQOHgTKlgUMBqUtKVnccQfw5pu8KzdhQliX9ujRQzzzzDOlzWbzfCFESowsLBIJOX1LTU3tUbZs2S+3b9+uL1UqYQVdcuON7C398ovSlpQ8gkGgXTtg7Vpg2zbuNhwifr8fd999t2vTpk1f2+3252JoZUQknKckhGiSmpr65aJFi6QgJTL79wN793IApST+qFScxJuSAnTtCng8IV+akpKCn3/+WW+xWJ5MTU19PIZWRkRCiZIQwmwwGH6ZOnWqvm7dukqbI7kSP/zAxxYtFDWjRFO5MgvThg1cPsbvD/nSjIwMLFq0SJ+amvqFECIzhlaGj9KBUuc/jEbjxMcee6wYVZkvpvz5J5d2bdVK5r4lAh9+yIGVkyaFfelnn30WMJvNWwGoKQE0gCiBirwJIZpnZGQsOnDggC49PV1pcySXw+UCMjN59ycnh9NNJMri8wHlygENGwLLlwMi9M3qYDCIW265xbl9+/YhHo/nkxhaGTIJMX0TQuiMRuPMSZMmSUFKdLzewmlCXp6ytkgYjYZ341auBAYNCitMQKVSYfr06Qa1Wv2uEKJqrEwMh4QQJYPB8E6rVq0y7r//fqVNkVyN9HRg9WpApwMeeSSsN4Akhjz/PNC/P/Dhh8C774Z1ac2aNfH666+nmc3m6UKE4WbFCMWnb0KIm0wm0+q9e/fqypcvr6gtkhAhKixZ8tdfQOPGytojYYJB4IkngGnT+IOjWbOQL/X7/ahfv75z7969/X0+39cxtPKqKOopCSE0JpPp288//1wKUjJx7lzh1zJsI3FQqYBx47hW+quvskiFSEpKCmbMmGHQaDSfCiEUTTJVVJR0Ot2rN9100zU9evRQ0gxJuGRkABYLULs2L3yfPcvb0mG8CSQxwmAAhg0Dfv8d+PLLsC5t2LAhBgwYkGY2myfFyLqQUGz6JoTINBgMf+3cuVNXpUoVRWyQFIHp0zn36vzF7lGjeKFVoixEwC23cAunTZvCutTtdqNWrVrOo0ePPhEIBH6KkYVXRBFPSQihNpvN344cOTJNClKS0r07sGcP0Lt34c9at1bOHkkhQnCVys2bOeo+DLRaLaZPn27QarUTlUraVUSUUlNTn6tRo0a1vn37JsTunyRCypYFvviCd+QAjpMp6FcmUZaHH+bjjPDbvzVt2hRPPPGEzmw2j42yVSER9+mbEKKqTqfbsXnzZn3NmjXjOrYkRpw6BQwfzhUR77kH+O47uQCeCLRuzd5sdjag14d1qcPhQLVq1VynT5/uTESLY2ThJYm7p2KxWL587bXX0qQgFSPKlwfGjOEyGitXAnXrAn/+qbRVkhdfBI4eBdq3D6sTCgAYjUZMmTJFbzQaJwshNDGy8JLE1VMSQtyUkZGx+vjx4zqdThe3cSVxZPNm4MEH2XtasoRr/0iUY8YM4PHHgc6dge+/DysFBQCaNm3qXL9+/Qt+v39ijCy8iLh6ShaLZfRbb72llYJUjGnUiGv8uFwcMyNRlm7dOML7xx95xzRMRo0aZdBqte8KIeJWbDFunpIQ4tYyZcosP3bsmD4tLS0uY0oUpEcPLm/icBQ2GZAoQzAI3HwzkJvL/ePC/H+0aNHCuWbNmpd8Pt8XMbLwAuLmKaWnp38wbNgwnRSkEsKdd3Ly7uTJMj9OaVQq7rR78CDw7bdhX/7+++8b0tLShgshtDGw7iLi4ikJIZqWK1du8dGjR2XftpKC281VKVeuBHr25J5lspa3chBxZ2OVirsah7m21Lp1a+dvv/02xOv1xry8SVw8JYvF8uGIESOkIJUktFpg6VLOwZo4kWt5f/qpTEVRCiGAgQO5nvdvv4V9+fvvv2/QaDRDhRAxXxCOuackhGhRsWLF+YcPHzZo5NpCyWTNGs5e37+fv1+/ntMgJPElPx+oVAm4/nr+H6jVYV3+r3/9y7Vs2bI3PR7P6BhZCCDGnpIQQlgslg9HjhwpBakk07QpB/E1asTf3347MGKEXGuKNzode0t//cUNB8Jk5MiRerVa/YYQIqbz8Jh6SkKI1pUrV5594MABQ0oEfwRJMcRu53y5mTOBF14APvigsDZTtCFij2D1ag5RqFiRAz1TUgCTiTvNlrQlhfz8wujuqVOBxx4L6/L777/ftXDhwnfcbnd4leTCIVbFvwEIi8WyfcaMGWEXM5cUcwIBogEDuNh9q1ZEBw9e/lyXi2jHDqLffydav55o926i3NxLnxsMEh04QDRxIlGfPkQtW/IYl3ukpxN9801MnmJCc+5c4d9mzpywLt2xYwfpdDo7ADPFSjtidmPg3uuvv97u9/vDetKSEkIwSDRhApHJRGQ0Er33HlFODv9861aiN98kuukmopSUSwtK1apEXbsSffkl0YoVRD16EFWsWPh7i4W///BDorNniXw+oqNHiTZtItqwgd+MzZvzuW+/rfRfI/7k5xPVr09UuTKRwxHWpQ899JBTq9W+ScnUzUQIIcxm886vvvqq1kMPPRT1+5/P6dO80yx3m5OUw4eBPn2AhQt5KlW1Kq8/CcFrUXfeybl0pUtz3FNeHnD8OK+LrFrF6SwAYDQCHTtyWkuLFlyA7mrb3n4/L8DPnMnR5716xfrZJha//w40b845cqNDX7ves2cPGjZs6MzPz69ERNHvHhELpQPQrHLlyvZAjHqC+XxEM2cWfig+/3xMhpHEk6wsoueeI2rRgmjsWKKTJ69+TTBItG8f0U8/ER07Ftm4LhdR7dpEd90V2fXJTu/eREKwtxkGDzzwgFOtVj9PyeIpWSyWma+//vojL730UtRXMM+evbjV2DPPcGxemNUZJBKme3dg7lxg3TquETVnDr/IVq3iLrQvvqi0hbHD6QQaNOCduc2bQ96VW758OTp37nzQZrNVo2iLSLRVDoApLS0t/2Qon3TncfIk0YMPEvXrR6TXEzVtSjRrFpHXy8sAffoQffcdUbVqhR7Sc88RqVT8dfv2/MEpkYTNoUNEZcqwx3Cp9avevYt3J+CffuLnOWxYyJcEAgEqV66cA8BNFG0NifYNhRBPt2nTxh7u3+XBBy/9eqhWrVB4AKLy5YkmTy4UII+HqH//wt/37Hn5zRmJ5LKcOEH0xhtEI0bwQvtvvxF9/TV/8gFE993HPy+uPPggkcHAmwIh8uabb/qNRuNESnRRSk9P3zYnzG3GAqEueDRvzru/Wm3hB9WePUTz5hHZbBdf//nnF15ftSrR/v1hmSCRXJpgkBctAaJbb1Xamtjx44/8HD/7LORLDh8+TFqt1glAR4kqSgBqWSwWl8/nC/mJbd1KVKpUoaBUqUJ06hT/Lhi8tAj9E4+H6MgR/nr1ar5P9+4Xe9y5uUSnT4dsmkTCBINEnTrxC6u4ekt9+vDze+ONsC5r1qyZTQjRlaKoI1FdiNbr9b179uyZEmr09tKlwK23XhhUO3UqUK4cfy0EB95ejdRUXo8EOEhXreZ6Vk2bAgcO8G7zM88A11zDj44dgZYtgSFDONhX5ohKrogQha2w161T1pZY0bs3UKcON4LIC32X/7nnnjNZLJYXompLtNQNgEan01mzs7NDUtiNG4nKliWqU4fo8OFCT+nw4bCE+pKsX3/xlE6v5/Wm/v2JKlUiatSISK3m33XqxLFkEsllmTePXyxr1ihtSezYtIkX+198MeRL8vPzyWAw5AOoSok2fQPQsUGDBiFMtojsdqKMDJ6qZWcTLV1aKB5ud8h/j6vSty9R69b8erqU6OTlEb3/fuHYx49Hb2xJMWP4cH6RhLEQnJQ89RSRRkP00ENEL7zA8WNO5xUv6dOnj1ur1Q6nRBOl9PT0pV999VVIz/vLLy/80Nmyhb+vUCGky6NOQRrQO+8oM74kwbHZeNv3jjuUtiT2nDhBdMstF04z6tfnhfAJEy45ldm8eTMZDIYcACqKgpZEJXhSCFFBp9MdOnXqVJophEWgu+8GTp4EsrIKMwGWL+f29DfdVGRzIqJFC85YyM5WZnxJAnP0KFClCtC2LS9QlgSsVu5+8u67wH//yyk+AHDttcCxYxedXrNmTfuePXseJKIlRR06KgvdGo3miQceeCAYiiA5ncAff7AwnZ+a1KqVcoIEAA89BOzaxWlXEskF+Hx89HiUtSOeWCy8O3TwIC98r17Nu0nHj3MDgn/w3HPPGS0WS/9oDF1kURJCCJ1O91yfPn1CKpP5/vtc2qZr16KOHF06duTjnDnK2iFJQIxGPq5cCezbp6gpiqDTAc2acRMIgBN5/0H37t2Fx+NpI4QocmvkaHhKDQ0GQ0bTpk0ve8L06dx+6scfORm5SxcuPphIVKnC4QmffBLWjqikJPDRR3ysXj2iio3FhoISxn/8cdGvSpUqhdatW/uFEA8UdZgii5JarW7bqVOnFHGZMhG//84twH76CXj4YfaS3n67qKPGhldeYe906VKlLZEkDLt3syi1agXs3culVUoqBdPYAs/xH3Tu3NmQnp7euajDFFmULBZL53bt2l22mduIEZx4vX174c9uvLGoo8aGgpI+LVoobYkkYdBqeS3p+uuVtkR5TCZeCHa7L/nre+65B/n5+XcJIYqkK0W6WAihczgcDVtc5l1MxGs1w4ezEP3wA1eISEQ2bAAmTOC1vdRU4OuvC5cPdu3iVlky8rsEct11vAvyww/AuXNKW6MsRPzmcDgu+evKlSujTJkyBKBhEccpUsDkPfXr17dGPVYizgQCHIqh1XJ11vPz8G69tfB7k4mrFrzwAodz/Pmn0pZL4sK2bRzpPHCg0pYoy6ZN/EaYOfOyp/Tp08etVqtfpiLoSpE8JZ1O96+OHTteeoL5D4YN41ikROSvv7hHX5Uq3Gxj8GCujHrkCDfDGDKEc/KeeAKoWZMLylWsyO3ZL+PJSooT9epxl99PP+XgupLK0aN8rFHjsqe0a9cuzWKxFGmxu0jBkxkZGfvmzZtXrVmzZlc9V60GXnuNxSnRGD4ceOMN/vruu4H58zkk48wZFp/9+3lXFGCfaejQwsX6Pn3Ywz90CEhL44KFZcoA6en8P9y5k9dKGzUC3nqL19ckSciZM0BmJrdoWrcutEzx4sYnn3BbrBMngAoVLnmK1WpFuXLlvF6vN52I8iMaJ1IXC0AZrVbr8Xq9V/X6vF72+oYPD8tZjBu//caVK199ldN8hg0rnLL9/POlrzl58sJI/HLliMzmC38GcEONZs24KUd6OtGnn3KNcUkSsnQpVxzs0kVpS5Th8cc53eYqJV7r1atnBXAPxTv3DUCXVq1ahZSAe+YMj/TppyE/fUUpEJTZs6983quv8lrU+TXXPR5eb9qx48IKmDt3cnIwQNSwYfEty1Pseecd/id++63SlsSfdu2Ibr75qqe9/vrrAZ1O9yHFW5TMZvO0jz/+OKTnsn8/jzR5ckinK07VqmxvtMsyB4NcZbN8eU7EfvllLkoXzcoIkhjj8/EbU6/nf15Jwe3mOuYPP3zVU1etWkXp6el7KZ4L3UIIEQwG27Zu3Tqk861WPprNkYwWf9as4XWgaHeTFgJ44AFeK33wQWDkSG67lZ7ODTPs9uiOJ4kBKSkc11KpEv8zjx9X2qL4sHYtr6v16HHVU2+77Ta43e4qQogyVz35EkT6tqum0Wj0mZmZIZ1ckL+XLG3br7kmtgGeZcpw/8MqopBMAAAgAElEQVScHGD2bE67+fBDXkf95ZfYjSuJEhUqAD//zC/szz5T2pr48PnnHEh6111XPVWj0eCOO+7wAGgVyVCRekqt27Zti8ullvwTtZqPu3dHMlrxpUwZ4P77Oc9x3TpOzO7WTQZpJgWZmZykunix0pbEnj/+4Dyx117jF2kIdOzY0WQ2mztGMlxEopSRkdG5ffv2Ibd+3LiRj/feG8loJYPbbgNeeomnurJ8SpJQvz6H+18mwrlYkJ/PcS/lywMDB4Z8WevWrREMBtuKUD2X84hIlDwez82hxCYBQCAAjB3LVQFCnO2VWJo04eOffyprhyREOnTgDPPi6i1Nn86Bo9u2ARMnXjYR91JkZmaCiEwAyoc7bNiiJISw+Hw+Q5UqVUI6f+5c7ihSnDsfR4vMTMBgkKKUNNSrx7sXmzcrbUn0mTaNF7UtFi4ydt99YV0uhED16tXdAGqFO3QknlKtqlWr5qtC3JqaMIEraHbqFMFIJQy1miPIDx9W2hJJSJQvz2/Wzz67ZDXGpGXRIs6patGCd906dIjoNg0aNEgFEPb8KCJRqlu3bkjXEQG//cY7pyW5NlaoHD/OlQlCnBlLEoH//IerAq5dq7Ql0WHrVqBdO6B2bc63SrtsVaKrUr9+fZ3BYKgf7nVhi1JaWlq9hg0bGkI59+hRnnLLtaTQKAgHaNdOWTskYXDDDXw8cEBZO6LFJ5/wccIEXksoApmZmdBqtY3DvS5sUTIajY0zMzNDWlGfPp2Pt90W7iglD58PGDWKBbxuXaWtkYRM5cpArVpcbynZ+e03LiT2/PNRedPWqlULHo+nerjXhS1KPp+vVihBk3Y7V53s2JEz5CVX5ssvudrq++9f2OVFkuAIwbEuf/7J6xXJyvLlQPfuvE42YkRUbnndddfB4/GYhRChb9shTFESQqS6XK5y1atfXfyysliYevYMZ4SSycmTwJtvAi1bAu3bK22NJGz27OFo52QVpWXLgNatebq2eHFhnZ4iolarUblyZReAmuFcF66nVL1cuXKutBAWvwqaOtYKe0OwZEEE9O7Na29jx0ovKSlZv553c6KdLBkvZszgxNS//uKA0ChSt25dgTB34ML9K9aqXbt2SB8HmzZxrFXBOqDkYoi4quXcucA770gBT1pKlwbOnlXaisjw+4Fff+XqhmEER4ZKo0aNjKmpqXXCuSYsUVKr1bUbNmwYUnrJX39xx9uCvDfJxQwdytP3nj25oJ8kSWnalBeJ/X6lLQmf2bO5kuRjj8Xk9pmZmcJkMoXV+zosUTKbzU3q1KkTUsTRwYNXLOVb4rHZuDLAww8D48cnr+cvAfCvf3GsUrI1DDx3Dujbl3vZhRmxHSqZmZnw+/0xnb7VqRXCHMPjAU6dAjIyknftL9ZMmcJ5nC+9JNeRkp6CumLJlgO3di1PO0eNitmUpkaNGnA4HBWEEJpQrwlZlIQQKofDUSWUcICCom6jRnEKjeRCAgGOUbvtNu6IIklyjh3jY7LFvsycyQ0QYli+Q6fToUyZMm4AIa8uh+MpZWg0GrKEUE/l/B3FRG2rpBQeD8em7d8vk5SLDX/9xcdka+k9bx7w6KNFjty+GlWrVg0AuC7U88MRJbNerw9pJc9k4gV9QKaY/JMJE7iIX9++QOcid12XKE5ODvDqq7yVfscdSlsTOsEgpxF4PDEfKj09XQAIuSdVOKJkMhqNgVBPLqg33bRpGCOUAH77jT9Qx46VO5NJz08/8T/z9GlOz0imf6hKBTz5JE/h8iNrzxYqFotFjRiJktlkMoW8bD17Npd7vfXWMEYoAWzaJP8mxYLTp7neUP36HDzZOOy8U+Vp0YK9pR07YjpMRkaGGkDIbUPC8pTMYbQj+eUXznuTJUsKIQL++1/O4ZQkOT/+yD3bx49PvgXuAvbv52OMI5wzMjJShRCxmb79PTe8KoEA78BdF/LSVsnA7eZHjNcVJbHG5+PCbnXrcvXJZOWPP3j6WapUTIexWCwqjUaTHur54fgx5vT09JDOLwhslTFKF6LT8et41SqlLZEUiaVLObkz2cuVnDwJ1AwrVzYiTCYTtFptyD3gwvKUMjIyQhKltDTO49qwIYy7lxDuuYfbKcVh00MSK5YvBzQajuROZtxufh4xxmQyISUlJWRPKZzgSXNGRkbI7STr1OHuM8mYDhRLGjXi14Ksw52keDzA1KlAmzaAPuQuY4mJyQQ4nTEfxmw2QwgRfVHSarWlTKaQ16rQsiVXCG3bVgrT+RT8LYpQ+liiJIsW8c5bv35KW1J04tSl4m/diP7uW2pqaqlwdt/69QM+/ZQ93VmzQr6s2JOTw8fSpZW1QxIh337L/efvvltpS4rODTcAR47EvOyK2WxGMBiMviilpKRkhOMpAcCzzwLVqwNvv81TFgmvLRqNMSldI4kHM2dyydjUkFcyEpfHHmPXfdy4mA5jMpkQCARC3nMOZ6E7PRxPCeAA108+4dis114L69Jiy8mTQIUKSlshiYgzZ/h49KiydkSLGjW4REWMPQaz2Qy/3x/yAlxYu2/hekoAb1D8+98sxgXVA0oyp08D5copbYUkIo4c4ePzzytrR7RYuZLjdgr6xccIk8kEn88XcuHvkEUpGAzqtFptREb17cvpNcke1hENHA4uhyxJQipV4mO1asraES3mzeMdxBg3GtRqtfD7/dGvp6RSqRzOCLcPb76ZPcVvvono8mKFy5X8O8klFpOJ43oWLFDakuiwZg23Y47xVrDT6YRGowk5Mi+cOCW7vSD1P0yEAJ55hiOZs7IiukWxID+f23JXqaK0JZKI0Om4mPrixckfaOb18mJvHPL27HY7NBpNyKUIwpm+5UYqSgBXSdBogEmTIr5F0vP777ym2KaN0pZIIqZfP16HmTJFaUuKxs6dnMPXsGHMh7LZbEhJSXGFen7IouT3+/NsNltkVoHLmHTuzEXOisvmRbgsXcrCfOedSlsiiZjrruPI4G++Se7kzoJqmTeF1WgkIux2O9RqtSPU80MWJY/Hc64onhLA7YQCAY5fSub/Z6QsXQrcfrusEpD0PPggz8MPHlTaksjZtIlfiHFYtLfb7RBChCweIYuS1+s9Z7VaiyQlN9wADB8OzJ9f8hoKrF7NrwNZArcYUODqrl6trB1Fwe/nVuNx6O319wwr5ICgkC0iInteXp43EqPO5/nngebNeWp++nRR75Y8vPsuB0326qW0JZIiU6cOkJ7Ou1fJyu7dcas2aLfbEQwG80I9PxyZtOXm5hY5tVat5tZLdjvXqy4pHDvGLZVkOEAxQKXif+batUpbEhkbN/JWeKdOcRnOZrPB7/fnhnp+OKJkz8vLC7lxwJVo1Ig9x5JS7Mzp5PikOFSJkMSL22/nHazckN9ricHKlZxMXKkS8NxzcRnSbrfD7XaHnPUblqdU1DWlAlJTgQ4duNXQ119H446Ji9PJIQCHDgE9eyptjSRqdOzIuzUffaS0JaFhtwMTJ/KLsUwZjk8pE3IxyCJhtVoDPp8vJtM3uzWKyWuTJ3O346eeAkaOLH5VBHw+YNs2bnixdi3w5ZfAI48obZUkajRsyIujS5cqbcnlCQSAjz/m9jkWC0cw167NraHiGMGbm5vrAxD93TcANofDEbWu93o9MHcu8MADwCuv8CLw7bcDH3yQ3Avg773HoR9GI9CgAfDzz0CrVixOkmJGqVJAEWL3YorTyZ/6AweyR/fGG7xesnlzXAImzycvL88PIOQ/VFieksPhiGq3vbQ0TtJdsoS9iEAAGDSoMHQgxj3yok5WFjdL3bSJdxmnT+dQlmXLeA1NUszYuZOTOhON1av5U3HFCn4jbdgAvPUWe3Yian5FyOTm5gYRhqcEIgrpAcCSlpbmDQaDFEuysogeeIAIILrhBqIVK2I6XNTYvJmoVCmismWJlixR2hpJzLFa+UU6bJjSlhRy6hTRO+8QpaSwbUOGEMX4/RoKjRs3zgPQikLUmrCmb0II7+kYz63q1OEp77JlHD7QqVNiT+f8fqB9e95R9Pv5Q+qee5S2ShJztm/nYyJ1xm3TBhgyhKdt586xl6SAZ/RPDh48qAGwP9TzwwmeJL1efzA7Ozsiw8KlVStec7LbY16ts0iMHMndgIcO5YajcWijJUkECjZ9YtzIMWROngS2bgXeeYdfkBkZSlsEALBarXA4HCkAQs54DSvGPBAIbN+1a1fYhkVKrVo8DZ45MzFz5Xbt4qn6I48Ab74Ztx1WSSJQvjwfjx9X1o4CLBY+BqISShg1du/eDaPReJSIgqFeE5Yo2Wy2v7Zv317kVJNw6NaN3/wF3nKiQMSL2Xo9MGaM0tZI4s6pU3wsW1ZZOwpQ/70HlQDTtfP5e2a1I5xrwhIlIsreunVrXPfEOnbk47x58Rz16syfz7W+hg6VNbdLJAULnQUlcpWmoIJHgrXJ2blzZ8Bms20M55pwU4R37dq1K6TW3dGiYkUupzt7djxHvToTJnBpneLQk1ASAQXxKrqQ6+HHloJllQSrH75582ZnIBAIayE6XFE6bLVaUx2OkOs1RYWuXbkm1Y6wnMDYcvAgUL9+XFqxSxKRgszqREloPHmSjwnWlWLnzp0CQOxEiYgCRqPx+O7du8MyrKh07w6kpCRWBdKKFbkteSIuwEvigOvv6q6JUrFv0SKO0K1bV2lL/ofX68WpU6f0APaFc13YFZ5UKtXOeIUFFFCuHIcIJNK6UqdO7LmNH6+0JRJFKPBI8kLOM40d585x1cQePRKqH/y+ffug1+tziCjkTiZABKJks9n+ys7ODnl7L1q0bcvT5oJ+gErTty9/KE2dqrQlEkVo0ICPiVDobdo0XuOKUymSUNm1axc0Gk3Y06qwRcnv9+/YtGlT3CfSBR1Ali+P98iXRqXiNmAyp62EUqcO73RMnapsbFBeHjB2LFCvXqFQJgg7d+4kp9P5V7jXRVKgN3uHAivO11zDx0TwlgsoX74wXEVSwhCCy1usXg2MHh2/cb1ejtpu0QKoWpUjyvft4/IaCcaWLVtcbrc77AjDSERp74kTJ/Q+ny+CSyPHbObp8h9/xHXYy0LEguQvcoFgSdLSpw/P4X//PT7jEQEPPwy8/jpP15o140C5Vas43y3ByMrKCgAIOwUk7JgjIso3m81n9u7dW7527drhXh4xKSmczjF5MntL6elxG/qSrFwJrFvH1TMlJZiqVbkAezywWjkh9NZbE+fT+TL4/X4cOnRIhwhEKaL+KiqVatXKlSsjubRI9OrFHxATJsR96Iv4+msWxiefVNoSiaKUKxe/OXxB1Patt8ZnvCKwceNGaLXa/xJR2AsuEYmS1WqdM2/evKJ1poyAhg15Kv3VV/Ee+WK2bwfuuCNxAnolChEMxqV3GgDg7bc5x+2FF+IzXhFYsmRJ0OfzRRTEE+lfc9mqVatS/QosqLRoAezZU/ihoRTHjydO2pNEQcqWBXJyYr8D98cfXPh/4EDg+utjO1YUmDNnjsPlcv0SybURiRIRndRoNCc3bgwrzy4qtGzJ631K5sIRcWedBIpTkyhFrVq8I3boUOzGyM/nXvfXXss1chIch8OB7du3awFE1EQtYr/T5/PNW7JkSdyDKJs35w+K77+P98iFHDvGu27XXqucDZIEoaCqX6yyHLKzOX1gyxbeVUmwKgCXYtWqVTAajVlEFFE8Y8Si5HK5fpkzZ058M3PB4SFt23J33ThHJfyP9ev52KSJMuNLEog6dfgYi9i9WbM4IHLVKt7diVNH26KycOFCr91u/znS64uyQrdq+/btunhXDAA4D87h4G4x8SY/n7vVXHNN3DvVSBKRrCw+FkT3RoulS4EuXfiT78gR4Omno3v/GDJv3jyPz+dbFOn1EYsSETmNRmPWKgV6bzdrxsfffov70Jg3j3PwvvySW0RJSjB79vB6AsDdI6J53y5deGq4cGFSVRE8ceIETpw4oQYQdnpJAUXay7Tb7bMXLlwY1/K4AJcNqV2b+8XFm99/51I6bdvGf2xJglHQ1vnpp6PXQGDNGp4KqNXcyTTB6iNdjaVLl8JgMKwhooi3I4skSj6fb/G8efPCKksQLZo3B/78M/71jA4fBqpXl8XdJCispXTbbUW/V34+F6Rv1oxjn5Ys4RdakrFgwQJXbm5ukfbGixr1tfHEiRMpJ06cKOJtwqdBA043OXgwfmMSAUePysoAkr/ZtImPb79dtPts3w7ceSfw7be85b93b8Jl/IcCEWHRokWCiIo0hymSKBFRQK/Xr1m2bFlRbhMRd93FHu4bb8RvzOXLeXH9qafiN6YkgSn4RDx6FFi7lrdld+xgd3rmTK7hfCUcDuDll7mh5aFD3IV16NDEqWYZJtnZ2fD7/U6E0XjyUhS5CUBeXt6s+fPn39GjRw99Ue8VDrVrAy+9BLz3HjBiBFClSuzHLAgF6NEj9mNJkoAXXuD6Nf/+N9C06aXPuflmIDOTs8mPHgX++1/etj14EPjoI04NePJJ4P33k75x4JIlS6BSqRYRFW1RRRTxegghbjCbzVlnzpzRaeK80LJ/P0+7P/0U6N8/9uM9+ig3H7XZYj+WJIlYswZYsYLzjlJSONz/5pu5lMSiRewxXSov6vbbuQ7S7bfH3eRYcMcdd9jXrVv3FBH9WJT7FFmUACAjI2Pr5MmT63dSILircmX+kPr22+jcb9Ei4LHHgE8+4WTslBRed/zhB2DwYH6MHBmdsSQlhLNnuaV2WhrvkKSk8G5d1apKWxY1Dh06hMzMTKfb7S5DRO4i3YyIivwQQjzZunVrOynAM88QaTREWVnRud/XXxPxkvbFj+bNiTye6IwjkRQn/u///s9nNBrHUzT0hKLgKQkhjFqt9vSBAwd0FStWLPL9wuH0aZ6y16rFlUmLWkXC6+WpYEGXkjZtOLq/eXMuMphgXZElEsUJBoOoUKGCMycnpzkRFTnPIiqFYIjIkZqa+vPUqVPjnqBbrhwwahRvfkyfXvT7paZytPbq1bwpsmwZtw6vV08KkkRyKZYvXw6v13syGoIERGlNCQCEEE0rVaq08MiRI0YR53dvMMjxaydP8qaGWh2d+x48yAvpL77ImyMSieRiHnzwQdecOXNe9fv9n0bjftEsmbfWarVa165dG8VbhoZKxQvQR48CixdH777XX887bqNGceiBRCK5kNzcXPzyyy/qQCAwLVr3jJooERG53e7Pxo0blx+te4ZDx44c5jFpUnTvO2kScO+93EBCdi6RSC5kxowZSEtLW0JE56J1z6gWF/b5fFN++uknlRLlTFJTgc6deUs/mnWW0tLYY7JYeCdXIpEUMmbMGLvVah0TzXtGVZSI6ERaWtraH374IZq3DZl27ThGLdolTQ4d4ngoiURSyJYtW3Ds2DEvgKjmmUW9DUNeXt4nn376qSJl/e+9lz2ab76J3j2JuBqBLOgmkVzI+PHjvYFA4EsqQpmSSxGL3jC/7N69m3bv3h2DW18ZnQ647z6ewkVjU9Hp5HzJnBxO4pZIJIzb7cbUqVODbrd7YrTvHXVRIiKfEOLriRMnKlJBu3lz7g24Z0/k9yDissi1avHO2zPPAE88ET0bJZJkZ+7cudBoNNuJKOrFg2LSRc/lcn05YcIEv9tdtBSYSLj3Xj7+HGHZ8p49edH8rrvY81q9mqO7oxX7JJEUBz7++GOH1Wr9JBb3jokoEVE2gHVffvllnOtCcgmTOnVYTCJBo+Gt///8B9i4sbAeuEQiYX7//Xds27bNRUQx2dGKWb9hq9U6aOjQoW6XyxWrIS5LQdWIcxFETrzyCh/LlEm68sgSSVx46aWXHG63+1Uiikl9/piJEhFtJqJVY8eOjXs+XK9evEi9YkX41153HZe3mTgR8ChSfVwiSVxWrFiBrKwsRyAQmBqrMaKW+3bJmwtRz2KxrD927JjOGMfOnvn57OnceScXZQs3FW/RIo55+ugjLi4okUi4zNHNN9/s2LRp07PBYDCKgTcXEjNPCQCIaDuApWPGjIlqHMPV0Om4ENvCheHlrO3bB0ydyh2S1Wpg27bY2SiRJBtLly7Fnj178ohoRizHiamnBABCiEyTyfTXsWPHdOY4LtIEg8CDD/Iu3JQpwOOPX/l8IsBk4mkfAFSrxlVOy5ePva0SSaJDRGjQoIEjKyvrmWAwGKU6r5cmpp4SwDtxQohfP/7447ims6pUXML2rruAZ5/lsslXQgiu/w6wt/Thh/HvKSeRJCq//vorDh06dJaIvo/1WDH3lABACFHDaDRuPXr0qC49PT3m453PwoUcuzR//tU7K3u9vMhd0M4L4ORemYgrKckQEerUqePYvXv3vwOBwE+xHi/mnhIAENFelUo194MPPoh78Y/mzbk++8CBvAB+JVJTgR9/5GMBs4vU61MiSX7mzp2L48ePnwoGg3F5N8TFUwK4FZPBYNhx5MgRbalo9V0PkSVLuNb26NEsTkJceUfu22+Brl3569xcIM7OnUSSMASDQdSqVcuxd+/e7kQ0Nx5jxsVTAgAiOqBSqX4YOXJk3HPiWrcG2rYFBg3iXbXmzYErZcA8+iiLUrlyvPgtkZRUZs2ahVOnTh0DMC9eY8bNUwIAIUQVvV6/+9ChQ9qyZcvGbVyAE3QbNy7cXRs5kkvoXo4ffwQefhiYMaPQa5JIShKBQADVq1d3Hjp06GEi+jVe48bNUwIAIjqiVqu/GTx4cNxjpW+8kUuQnD3L2/3TpgGBK0RPPfAAtwb/JCYphxJJ4jNhwgQ6d+7cPgAL4zluXEUJAOx2+0s//PCDY0UkOSBFRKfjxqTvvANs3w58f4XNTZWKO+9u2CBrc0tKHseOHcOgQYPcNputO8VzOoU4T98KUKvV7StUqPD93r179Xq9Pu7jB4NAo0Zc5nb3bqBChYvP2bMHqF8faNWKU1UkkpICEaF169audevWfeB0Ov8v3uPH3VMCgEAgsMDhcCwaMmRITLKMr4ZKBXz3HWCzAY88cunAyoMHOSG3f//42yeRKMnMmTOxYcOGky6Xa7gS4yviKQGAEKKsXq/ft2LFCvMtt9yiiA3vvgsMHQrcf//FUzmbDahYkUMJZs4EtFpFTJRI4kpOTg6qV6+eb7PZ7iKiP5WwQRFPCQCIKMfj8fTt1q2b0+tVxGHCa68B99zD6Sj/bMBiNvPu3M8/AzVrcjnczz5TxEyJJG707ds3PxAITFBKkAAFRQkAAoHAzJycnD/eeecdxZaSW7bk41NPXfy7N98Eli/n9koLF/JUbmFc9yEkkvgxb948LFq0KM/pdL6qpB2KTd/+Z4AQlXQ63a4NGzYY6tatG/fxn30W+OIL/vr4ceCaay59nscDVK8O1K0L/Bq3iA2JJD5YrVZUq1bNdfbs2fZEtFJJWxT1lACAiI75/f5B3bp1cwauFDgUI85f5N637/LnpaWxxxTN7rsSSaIwcOBAt9fr/UFpQQISQJQAwOfzjT98+PDOjz/+OO6lc3v0KPx6wIArlyvJz+dYJ4mkOLFixQp89913Lrvd/rzStgAJMH0rQAhRXa/Xb9u2bZuuWrVqcRvX5+OOurt2cY83gKO5f/zx4qTd++4DsrP5cX4lAYkkWXG5XKhevbrr1KlTXQKBwHyl7QESxFMCACLaFwwGhz722GPOeAqlRsOL3O+9V5hSMmsWN7T8J927AwcOAOPGxc08iSSmDBkyxOt0OhcliiABCSRKAOB2uz/csWPHfiV241Qq3l27/nr+/q23uOjb+RR0i7rS2pNEkiwsWLAA48ePd9pstt5K23I+CSVKROS32Wz3jhgxwr5gwYK4jy8Ep53ceCN7Q716Xfj7zp35KHPhJMnO7t278eijj+a7XK77iChHaXvOJ2HWlM5HCHG70Whc9ueff+pq1aoV9/HPnQNKl+Yo7r/+4moBBXTuzI0uv/mG15gkkmTDarWifv36rpMnTw70eDzjlbbnnySUp1QAEa3zer0D2rRp48rLy4v7+KVKsccUCHAL8I8+KvzdO+9w+kmHDkDv3pF14ZVIlCIQCOChhx5y5ebmzkxEQQISVJQAwOPxTMzNzZ3+0EMPuZSIX7rxRuDPP7li5X/+U1hCt0wZbr2kVgPjx3OxOIkkWRgyZIhvw4YNO+12e1+lbbksRJSwDwAas9m8YfDgwV5SCJ+PaNAgIo5gIhKC6Jtv+OsGDYjy8pSyTCIJj++//54MBsNpAGUpAd7fl3sk5JrS+QghyhoMhqyvvvqqXJcuXRSz49VXL91tNxDgnTuJJJHZunUrmjZt6nI6nc2IaLPS9lyJhBclABBCNNTr9WvWrFmjb9iwoWJ2EAGTJwMbNwJjx/LPRo8GXnxRMZMkkqty5swZ1KtXz5WTk/O03++PaXfbaJAUn/FEtMXj8Tzdtm3b/Jwc5XYvhQCefBIYMQKoV49/JtsvSRIZn8+Hjh07uux2+5fJIEhAknhKBZhMptF16tTpu3r1ar1Go1HaHPj9snuuJLHp16+fZ9q0aX/YbLa7iSj+O0YRkBSeUgEOh+Pl7OzsDQMGDIh7N5RLIQVJkshMmjSJpk6dmmOz2e5PFkECkkyUiChgs9nunzZt2sm3335bxlVLJJdhzpw56N+/v8PhcLQhovgH+xWBpBIlACAiq8PhuH3UqFEnpDBJJBczZ84cdOvWze5yuVoRUbbS9oRL0okSABDRCYfDcasUJonkQv4hSBuVticSklKUAClMEsk/KQ6CBCSxKAFSmCSSAoqLIAFJLkqAFCaJpDgJElAMRAmQwiQpuRQ3QQKKiSgBUpgkJY/iKEhAMRIlQAqTpORQXAUJKGaiBEhhkhR/irMgAcVQlIALhOl4nz59PD7ZQVJSDCAifPbZZ8Fu3brZiqsgAUmWkBsuQgiL2WyeXbt27Vvnzp2rL1u2rNImSSQR4fF40KtXL/esWbNOOByO1kS0X2mbYkWxFiUAEEKoDQbDewaDod+iRYt0StZjkqrQn5EAAAiKSURBVEgi4cSJE2jfvr1r//79q2w228NE5FDaplhSLKdv50NEAYfD8dLZs2efbNq0qeu7775T2iSJJGQ2bNiAevXq5e/evfsDm83WvrgLElACRKkAv9//ncvlavr000+fHjx4sE+JZgQSSThMmTKFWrZs6czNze3qdDr/j4iCStsUD4r99O2fCCHKms3mBU2aNKnz008/6dNl6UhJguH3+zFw4EDP119/fdbpdLYhoh1K2xRPSpwoAYAQQmM0Gj/LyMjosXjxYr0SDS8lkktx9uxZdOrUybV9+/ZNNputIxHlKm1TvCkx07fzISKf3W7vferUqedvvvnm/Pnz5yttkkSC7du3o169eq6tW7dOsNlsLUqiIAElVJQK8Hg8Ex0Ox91dunTJHTZsmL8keo2SxGDWrFm4/fbbXadPn+5tt9tfSKbytdGmRE7f/okQ4lqTyfRLnTp1qk2bNs1QrVo1pU2SlBByc3MxYMAA96xZs+wul+teIvpLaZuUpkR7SgUQ0XG73d54y5Ytb9evXz//gw8+CMrdOUmsmT17NqpVq+aaM2fOdJfLVU0KEiM9pX8ghKhuNptnXHfddbVnzJhhqFu3rtImSYoZJ0+eRK9evVwrVqw453A4uhHRaqVtSiSkp/QPiGifzWa7bdeuXS/ecsstzjfffNPv9XqVNktSDCAiTJkyBTfeeGP+8uXLv3A4HDWkIF2M9JSugBDiWrPZPKVs2bK3zZgxw3DLLbcobZIkSTl8+DCeeOIJ16ZNm47b7fZHiWiT0jYlKtJTugJEdNxms7U+dOjQMy1btrQNHDjQ63K5lDZLkkQEg0GMGTMmWLt27fz169e/Z7fb60hBujLSUwoRIUQZs9k83mg0tp02bZq+ZcuWSpskSXB2796Nbt26Offt27fPZrM9SkS7lLYpGZCeUogQ0Rmr1frAyZMnu3To0OHsk08+6c7JyVHaLEkCkp+fj2HDhvkbNWrkysrKes1mszWWghQ60lOKACGE2WQyjQoEAo8/99xz6sGDB2tKly6ttFkShXG73Rg/fjwNHTo0PxgM/m61WnsR0WGl7Uo2pCgVASFEFZPJNCwQCDwycODAlEGDBqXIBN+Sh9frxVdffUVvvPFGvt/v/8Nqtb4k140iR4pSFBBCXG82m98NBoOdBg0apBk4cGCK2WxW2ixJjPH5fJg6dSpee+01l9fr3ZSXl/ciEW1Q2q5kR4pSFBFC1DCbze8R0b2vvPJK6oABA9RGo1FpsyRRxu/3Y/r06XjllVec+fn5WVar9UUiWqO0XcUFKUoxQAiRabFYRhLRPW+88Ubas88+q9Lr9UqbJSkigUAA3333HV5++WWn3W7fY7VaBxLRb0rbVdyQohRDhBD1LBbL+0KIu4YOHart3bu30Gq1SpslCZNgMIiffvoJgwcPduTm5h60Wq0DASwn+eaJCVKU4oAQoqHFYhkN4I6ePXum9OzZUyMLyyU+J0+exJQpU4Kff/65y2q1HrPZbAMBLJJiFFukKMURIUSmXq/vRURP1axZU/Tv39/08MMPw2QyKW2a5G98Ph9++eUXfP75547Vq1erU1NTf7bZbJ8DWCvFKD5IUVIAIYQGwL/S09MHuN3upg8++GCwb9++ujvuuANCCKXNK5FkZ2dj4sSJvq+++sovhNhrtVo/IaIfiMiutG0lDSlKCiOEqKDRaJ7QarX9LRaLpV+/fvonnnhCVbFiRaVNK/bY7XZ89913GDNmjH3fvn1BIpqUn5//JRHtVtq2kowUpQRBsIt0u9ls7uf1ejs3a9Ys0K9fP2P79u2h0WiUNq/YQET4/fffMW7cuPxZs2aptFrt6ry8vE8BLCQi2d89AZCilIAIIYxCiIctFsvzwWDwxnbt2tF9992nv/vuu3HNNdcobV7SYbVasWLFCixcuNAzZ84cv9PpzM3Pzx/j9/unENEppe2TXIgUpQRHCHGDEKJNenp6Z5fL1axixYqBDh06pLVt2zb1rrvuggzOvBiv14v169dj0aJFwTlz5jj27NmjNRqNm61W6+xAILAYwBa5aJ24SFFKIoQQagA3aTSaNiaTqbPD4ahbt25dd6dOnYxt2rRRNWnSBCkpKUqbGXeICDt37sSSJUswd+5c+7p169J0Ot1ht9s9Lz8//xcAa4jIrbSdktCQopTECCEMAJrr9fr2qamp93k8nop33nmnr0OHDsaWLVuiRo0axXI9KhgM4siRI1i9ejXmz5/vWrx4sfD7/Q6VSrXIZrPNBbCCiM4obackMqQoFSOEEOUB3G2xWDoFg8G78vPzy1SsWNFVu3ZtNG7c2JCZmanKzMxEzZo1kyI2yuPxYO/evcjOzkZ2djZt2bLFuX379uCRI0f0Go3GkZaWtiE3N3c2ES0BcEBOyYoHUpSKMUIILYAbAdRSq9V1LBbLTcFgsLbD4ahkNpt9NWvW9DVs2FBXr1691Fq1aiEzMxPly5ePe6xUXl4esrOzsWvXLuzYscO/efNmV3Z2tionJ0er1+tPazSaXXa7/S+v15sFIBvAbiKyxdVISdyQolQCEUKoAFwHoJYQItNsNjdWq9X1XS7XDUSUqtfrfUaj0W80GsliscBisQiLxaLOyMhQZ2RkpFosFpXJZILZbIbJZILJZILRaER+fj7sdjtsNtv/jjabjfLy8ny5ubm+vLy8QF5eHv39e+F0OtVOpzPF7/fDaDQeAZBltVr/CgaDOwHsArCfiGQrmRKGFCXJBQghLADMAEyXOgohTBqNJkOr1ZZOSUlJV6lUGURkIiKjECIfgJ2I8gKBQK7H4znn8XjO/e3V2AEUHP/59Tk59ZIUIEVJIpEkFLJxgEQiSSikKEkkkoRCipJEIkkopChJJJKEQoqSRCJJKKQoSSSShEKKkkQiSSikKEkkkoRCipJEIkko/h9hiUgm0b/laAAAAABJRU5ErkJggg==",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAASUAAACZCAYAAABpGCzNAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAApE0lEQVR4nO2dd3iUVfr+7zM9vUwykwSQ0IvUZelSFFEEabIKqyArNpCVtigLiAqWBQsq4FcQsPyM1EUpikA2VCkSiUAIRSCBAMnMJDMppE2mPL8/DrMBl4RJmJl3kpzPdb3XO8m85R5I7pz3Oc95HkZEEAgEAn9BJrUAgUAguBlhSgKBwK8QpiQQCPwKYUoCgcCvUEgtQFB7YYwpAKhv2uwArDe2chKzKIIawMTPTf2AMRYIQAsg8sZeq1AoosLCwhpqNJo4uVyuI6IwAGqn06kiIrXD4VABUMhkMsYYk93YM8aYjDEmUygUUKvVpFarSaVSweFwwGq1wmq1MpvNxpxOp5OIyOl0ktPpJADOG5BcLrfKZLJyxlg5Y6wMQCmAHJvNln39+vVrxcXFRgDmP2wFROSU4t9P4DuEKdVyGGMMgB5AYwANAgMDG4eHh7dUKBRNHQ5HA6fTGa5UKpVBQUFMq9WSTqeTxcTEKPR6vVqn02mioqKYVquFVqtFcHAwNBoN1Gr1f/cKhecH006n02Ve/91KSkpgsVhgsVhgNpthMplsBoOhzGQylRuNRkdOTg4KCgpYeXm53eFwlCqVSgOAy8XFxefz8vIuENE1AFcAXCaico+LFvgMYUq1AMaYHEATAK0CAwNbRkZGdmKMtbLZbHEqlUodGxuL5s2by5s1axYQHx8f2KBBA7i2sLAwcN+qO1itVmRnZ+PatWu4du0aMjMzbRkZGcXnz5+3ZWRkoKSkxAYgXy6XXygpKTllsVhOEdF5AGeJqEhq/YKqEabkR9yI0bRijLWLjo7uqlaru9hstqYajSagadOm6Ny5s6Zt27YhzZo1Q9OmTREXFwe5XC61bL8kPz8f6enpSE9Px7lz56wnT54sSk1Ndebn5ztkMpkRwMnc3NwjVqs1FcBJIiqQWrOAI0xJIhhjwQD+FBQU1DUyMrKvzWZrr1arg9u2bYuePXuGtG/fXtOmTRs0bdoUSqVSarl1BiKCyWTCmTNncObMGefRo0fzk5OTHWaz2a5QKC5ZrdZDOTk5BwEcA3BFBOt9jzAlH8AYUwJoHxAQ0CsqKuohm83WMSwsLKBHjx6K3r17h3fq1Im1a9cOAQEBUkuttzidTly6dAknTpzA0aNHi/fv31+SkZHhlMvll61W656cnJy9AI4SkUVqrXUdYUpegDEWCqCHTqd7WKFQDFQoFDFdunRhAwYMiOjevbu8Q4cOUKlUUssU3AEiQmZmJpKTk7F///7CvXv3WnNyckoVCkWy0WjcarPZfgaQIUZTnkWYkgdgjAUBuC8mJmYEgIFhYWGhDzzwgHLgwIHhvXv3hk6nk1ihwFNYrVakpKRg79691h07dhT8/vvvpFAofsvJydlktVoTieiy1BprO8KUagBjTAags1arfUytVo8MDAyMGjhwoGrIkCFhffr0QWhoqNQSBT7C6XTixIkTSExMtG7evLkwPT3dplAoDmZnZ69xOp27iahQao21DWFKbsIYC1cqlY/o9fqn7Xb7n7p37y57/PHHtQ8//DCLioqSWp7AT7DZbDhy5Ai2bNlStHXr1tKioiKT1Wr9t8ViWQ+ekiB+4e6AMKUqYIw1Dg0NfTw4OHhcUFBQ7OjRo4NGjBgR2LlzZ8hkYtmg4M4YjUZs377dmZCQYE5LSytnjO0wGAzfAPiZiBxS6/NHhCn9AcZYk4iIiLEajWZcw4YNw8ePHx8xfPhwRcOGDaWWJqjllJWVYffu3fj2228te/bssTPG9mRlZX0OYD8R2aXW5y8IUwLAGIsJDw8fHxAQ8Fzjxo3DJkyYEPnYY4/JtVqt1NIEdRS73Y79+/fj66+/zktMTLQxxpKysrKWAjhS3x/x6q0pMcYCAwMDnwgPD385MjKy0cSJE8NHjx6tFPEhga+x2+3Ys2cPPvvsM8uhQ4fKiGi9yWT6lIguSq1NCuqVKd1YvNotNjb2H3K5vP/48eODJkyYENi0aVOppQkEAICSkhJs2bKFPvnkk9zMzExjfn7+4tLS0g1EVCy1Nl9RL0yJMRYaGho6ITAwcGrXrl1Dpk2bpu3fv78IVgv8mszMTKxatarkyy+/LLbb7bsMBsO/iChNal3epk6bEmOsbWxs7DylUjlg0qRJIc8995xGPJ4JahsOhwM7duzAokWLci9cuGC0WCxvWa3WTXU1OF7nTOlG8bGHo6Oj34qPj2/8+uuvRz300ENiVCSoE1y6dAmLFy++vmHDhiKbzbbcYrEsJaI8qXV5kjpjSowxZVBQ0Ljg4OB5Dz30UPjs2bPD27RpI7UsgcArFBUV4YsvvihfvHhxgdVq3WIwGOYT0VWpdXmCWm9KjLGAiIiIv6tUqunjxo0LeeWVV4LFWjNBfcHhcGDTpk3O119/3VJYWPhzdnb2K0R0QWpdd0OtNSWXGanV6hkvvfRS2NSpUwPEmjNBfYWIkJiYiJkzZ5pzc3MPZWdn/+NGtc1aR60zJcaYMiwsbGJAQMCcyZMnh02bNi0gODhYalkCgV9AREhKSsKMGTPMOTk5+wwGwzQiuiK1rupQa0yJMcaCgoLGBAcHL3rmmWciZ8+eHSRGRgLB7SEi/PjjjzR9+nTz9evXNxuNxn8SkVlqXe5QK0xJoVB0j4qKWj148OCGCxcuDBMxI4HAPRwOBxISEuxz587NKykpWZyXl7fY37u9+LUpMcZiY2JiVjZv3rznihUrItu2bSu1JIGgVlJSUoJFixaVrFixIjc/P39iWVnZT1Jrqgy/NCXGmCIyMvKV4ODgGcuWLYscOnSoSDISCDzA1atXMXHixLxjx46dMBgM44koU2pNf8TvTEmhUHTRarXrnn322dh58+YFiWL6AoHnSUpKwvPPP597/fr1xbm5ue/5U20nvzElxliAXq//OC4u7i/ffvttpEh8FAi8S1lZGebNm1f0zTffXDUajY8T0SmpNQF+YkoKhaKrVqvdOHv2bP2UKVM0YkmIQOA7Tp48iSeeeMJsNpuX5ebmviX1qElSU2KMKXQ63b9iY2MnbNy4MbJFixaSaREI6jPl5eWYN29e0VdffZVhMpmGStmVRTJTYow10ul0Pzz33HPN5s+fH6RQKCTRIRAIKjh8+DBGjx6dk5+fP62wsHCNFBokMaXAwMDBWq32y7Vr1+ruu+8+n99fIBBUTn5+Pp566qn8Y8eO/Wg0Gp8lIqsv7+9TU2KMyXQ63XvNmjWbsGXLlojo6Gif3VsgELgPEWHx4sWl7733XobJZBrky6UqPjMlxliITqfb9tRTT3V5//33g+VyuU/uKxAIas7BgwfxxBNPGHNzcx+zWq2HfHFPn5jSjfjRng8//PCesWPHKr1+Q0HtobQUyMwEjEbg2jXgwQcBMYL2KzIzMzFo0CCL0Wh81Ww2r/b2/bxuSkqlsmN0dPRPmzZtiu3Zs6dX7yXwQ2w2oLAQcLWrOn8eWLEC2LULuHSJm5L9D1Vd/SBNRXArxcXFGDp0aP6pU6f+Lycn5zVvtoHyqikFBATcr9frNyQmJkaJ6f56ysiRwPbtwFNPAbm5wLZtgFwO9O8PtG8PKJVAu3ZARAQwbBg/58oVQDT/9DvsdjueffbZwp07d24xGo1/IyKnV25ERF7ZgoKChrZo0cKclZVFgnrKmTNEfNxDFBZGFBBANGYM0cWLtz9+6lR+7NGjvlQpqAZOp5NmzJhxXafTbQGgIC94h1eSg0JDQ0fEx8ev3r9/f2RkZKQ3biHwd3bvBsaMAYKDgYsXAXfKzfz2G98vWsRjS0OGAI0a8e8RAT/9xONPw4YBcXHe0y6oFMYYPvzww+Dw8PAHly1b9gNj7FHydFcVT7tcUFDQo23btjVbLBZvmbXA33E4iOLiiFq0IDp+3P3zrlwhGjaMSKHgIyalkmj1av7exx9XjLqUSqLnnyfauZOoqIjI6SQ6fZro/Hn+WuAT5s+fX6zT6X4AICcPeohHDUmlUvVt2bJlbm5urjf/LQT+zubN/EcrIaFm55eXE507R9SpE79O69ZEsbH8dVIS0fDhRCoV/1ouJ4qMrDCsMWM8+lEEVTNnzpwinU63Fjfi057YPBboViqV7ePi4pIOHz4cHSeG1vWbb78Fxo4FjhwBunev+XXy84Hly4FvvgFOn+bB74wMQKHgM3qHDgE//wwkJ/Pjd+3i+5ISQJS88QlEhBdeeKFw69atq41G4wyPXfRuNwAxer3+2tmzZ73qyoJaQFoaUfv2RMHBRAaDZ65ZVkb03XdEV69WfdzOnXy09PnnnrmvwC0cDgc9/PDDeRERES+QPzy+AdBER0enJSUliYf5+k5CAv13pm37dt/f3+Eg6tePx6RSUnx//3pMcXEx3XvvvWaVStWHpH580+v16+fMmTNs6tSpGg8M3AS1lfR0oGtXnm+0YwfQvLk0OgwGoE0boGNHYO9eaTTUU65evYpu3boZsrOz/0RE2TW9zl1VUwsPD3++d+/eD0+ZMkUYUn3n3XeBsjJg507pDAkAYmKA+fOBffuALVuk01EPadiwIRISEvQ6ne4nxliNF7fWeKTEGGseHx9/ODU1NUo0g6znZGYCrVoBTz/Nl5BIjc0G/PnPQE4OD7bfc4/UiuoVs2fPLv7iiy+WGo3G2TU5v0amxBiT63S6Ez/88MO9Xbt2rcl9BXWJ8eOBf/+bz5A1biy1Gs7Jk/wR7qmngIQEqdXUK+x2O7p06WJJS0t72G63/1rd82v0+KbVamc9/fTTjYUhCXDlCrBhA08B8BdDAoAOHYDhw4G1a4F166RWU69QKBTYsGFDpFarXccYU1X3/GqbEmPsntDQ0Olvv/22eGYT8LhNWRnw8stSK/lfEhKALl24tvx8qdXUK1q1aoXJkyfHRkVFvVbdc6ttSjExMauXL1+uVavV1T1VUBdxdS12JS76E8HBPACfmytm4iRg1qxZgcHBwZMYY42qc161TEmhUPRq06bNnx566CFWPXmCOovLlE6elFZHZfTqxUul7Nwp6jT5GLVajaVLl2pjY2OXVuc8t02JMcaio6OXL1u2TCz7F1Sg1/PZrdJSqZXcnsBAYNQovlzltWo/SQjukiFDhrDY2Nj7GGP3unuO26Ykk8n6de/evWFb119GgcBFfr57pUmkYu1aXkblo4/4o5zAZzDG8PHHH2tjY2M/dPcct01Jr9d/8M4770TUTJqgzpKZyRfH+nNlUZmMj5JKS4ERI6RWU+/o06cPdDpdF8ZYM3eOd8uUGGNt4uPjG997r9sjMEF9YflyHrMZOlRqJVXj+tk9eBB44gk+YyjwGfPmzdPGxMT8051j3TKl2NjYf86ZMyfq7mQJ6iROJ+BwABaL1EruzLVrPH9p40bgySf/t2GBwGsMHz6cyeXy4YyxO9aUuaMpMcaUAAYNGjTII+IE7nHuHPD++7UgBNKqFd//9JO0OtwhLg44cYLHlr7/Hnjzzf89xlUuTuBRFAoFxowZE6BUKofd8eA7lRGQyWRDJk6cmO+bAgiCrCyif/2ropAiQBQdTeS3xTyPH68QumqV1Grcx6W5uJiX1N28mWjJEqImTXg9KFFW1+OcPn2aGjRocJDutp5SgwYN/n3gwAHff4J6yC+/EKnVtxqSa9NoiFq1uvV7hw9LrfgGzzxTIaq20Lcv1xsScvt/8OxsqRXWSeLj400AAqmqgVBVoyjGGHM4HPeJJpLe4/33+YL2hARg4UKAMeDYMZ7+A/Dwx9q1QM+egMl067lTpvCQjuSsWsVrGAHAmjXSanGXffv4Y1x8PDB6NLByJa8JNWsWf1/hlUY/9Z6hQ4eqZTLZA1UeVJVjAWj5wAMPmKRw1PrApUu8FZpcXvEHeuZM/t477/CvMzMrjnc6iWw2vp84kf7b2MNslkb/LdjtvNB/69ZSK7k7Bg0iuuceqVXUWX766SeKjY1dSTUdKalUqn6DBw8O95xPCm6mcWNe895s5l83asRbngHAnDm8/n2jm1YNMcb/gDPGY7QDB/LSQY884gexWbmclzA5exZYskRiMTUkI4NXzXz2WamV1Fl69+4NAP2rOqZKU9LpdA90795d6UFNgj9w773Ad9/x19u38zw/F1U15NDr+RrYF14Ajh4F+vTh9dVc/RwlYfJkXi5k6lRg2rSqnfLkSe6sI0cCffv6x4JeV4mTp5+WVkcdJiQkBCqVKowxVunzcZUPzg6Ho1P79u09r0xwCyEhfF+T+NBHH/Gy1Fu38rxAAHjgAWDZsoowj89QKoFNm4Dnnwc++YQP6T74AEhK4q2QiouBli351xs38nNateKZ1sOG8en67GweQPN1FQqnk5tShw48ziTwGq1bt8bly5ebAzh72wOqeraLi4sTUxA+YP16Hh+qTjPZm7l+nTfvOHOGaO7cit6Mp055VqfbOJ1EU6bcOpvFGJ9CBIiCgohmzya6fJkff/jwrccGBfFuub7k22/5vb/+2rf3rYfMnDmzEMAQqm5MiTGmDAgIEFMQPuDYMb5PT6/Z+cHBQOfOQOvWwNtvAz/8wL//f//nGX3VhjHg4495E8kRI4Cvv+ajoZIS/iGvXOF1jly1s3v0ABYs4NUr//xnPqJauNB3es+cAV58kev46199d996SqtWrUJUKlWTyt6vynSidDqdP0w412mMRv6o1bkzD1h7gp49ebWOTZv4U5Qks9uMcZMZO/bW7zep5Gdx3ryK1w0a8KL/vuD0aR7XCgjgj5RKEUL1NtHR0YiIiKi0dnJVge7wyEhROsmbWK08/FJWBqxfD2g82Khq2DBueEeOeO6aPuHqVSAryzcdSC5cAHr35uv2Nm3ibcEFXic8PBxKpTK6sverMiW1RqO5q75wgqr58ktg2zbgjTc8X/lj2za+37yZB8IlTxlwF9cIyZUn4S2MRuChh/h05y+/8OlLgU9Qq9WoamFuVabDZDLhSd4kIQFo1+7WJxdPkJvLOx4BwIcfArGx/PHw00+B8nLP3svjuIZ2vXp59z6zZvG8pLVrgaZNvXsvwS0wxgCg0pLaVbmOtbS0tLb8fa11lJbyKfxTp3jlD08SHAz87W+3fu/aNeDvf+dxq2vXPHs/j+LKi2jd2nv3WLmSB99nzOCjJYFPKSsrAxFVWj+5KlMqyMvLE6bkJW4OPicne/baGg1/NHzsMf51Xl5FUuXu3Tx00rQp7z7kDzmLt/DMM3xo9/e/e6fu99mzPOP0kUf4VKXA5xQWFsJms1X6fF6VKeWaTCbx/OYlHA4gNJS/9laS49q1vHx2eDg3oqNH+fc++ADo1g0oKgKGDOHf8xsCA/lI5uRJbh6eDoZ9/TWfGVy5suqUeYHXyMnJQWFhYWalB1SWwES8bIlInvQSv/7Kc/VWr5ZOQ35+RQWPIUOITP609Pr117mwX3/13DXz8nhtmGHDPHdNQbWZNWvWdQDDqSYLcp1O5/Xi4mJPG6UAFek6e/ZINzMWFsbbob39Nl/5MWSINDpuy4sv8r0nK1pu387zMKZM8dw1BdUmLS2tFMDFyt6v0pTkcnnqqVOnPC5KAERG8mB0QgLfpEKjAebO5VUJkpOB69el03ILcXE8u3vDBs9dM/pGaszp0567pqDapKWlAcDvlb1fpSnl5ubuTk5OFtXVvcTq1TzY/I9/8HxBKWnQgO+l1nELjz4KpKV5Lru7Y0e+9/u8iLpLSUkJSktLi4io0v+EKk2prKxs3/bt2/M8L00A8Ly9b77hCcWffCKtls2beUC8mVuduXxEw4Y8RcBg8Mz1XOkGvq5AIPgvR44cgUwmO1jVMXeaXTv922+/kdMvaq7WTdq04UtCvviChzuk4IcfeAb4P//pZ1Vgg4P5vrDQM9fT6bjRff11LUpxr1vs2LGjyGAwfF/VMVWaEhE5FQrFsWOuZewCrzBxIs/C/uIL3987M5OXqL73XmD6dN/fv0pciY379nnmejIZv2ZKinR/Aeo533//fanT6fxPVcfcMQ/JYDB8sW7duiLPyRL8kYEDeWG2WbN4VQ9f8uab/PdzyxZApfLtve/I1q1878ns7sJCnjnqydXPArfIyMhASUnJFSKqcuh7R1Oy2+0/rl+/vlQ8wnkPVy6f3e7bJOPTp3nm9xNP+FksCeBp6JMm8eD0sDv3L3SbS5eA5s09dz2B23z11VeleXl5n97puDuaEhGVEtH+vXv3ekSY4PY0bcrXoG7bxqt3+IKUFL6fM8c396sWW7bwZSaff+7ZQFdBQUUqvcBnOJ1OfPXVV8WlpaUb73SsW8tIsrKy/rVw4UIv15IQLFrEl3706wckJnr/fq56Zjab9+9VbRITgZgYoGtXz10zLY1XBvBFrSbBLSQmJsJms+0mojtmwrllSkR0LDU11ZiRkXH36gSV0qUL7/BDxOOxw4fzsj/eQqvle79JmHThcPCmlg0b8mdbT5CUxIN30dG824rApyxYsCA3Ozv7LXeOdXvBrcVimf3mm28W1FyWwB169eKxnoUL+WChS5eKGt6exlXCxO8KLrqWNp04cffXslp53e0HH+TB7V27/PAD121SUlKQkZHxOxG5tTzEbVMqLy/ftnPnzpxLly7VWJzAPTQaPhN36BAfNXlrqZZrpi8uzjvXrzGBgXzfsuXdXcdq5cXK160D5s/nbt+u3d3rE1SL6dOnW0wmk9sJJ26bEhGRxWL5+4wZM0SGt4/o1ImXF/rlF8/lD97MwYN8IsrvZscVCp7JefZszZeYOBzAgAHAjz8Cy5cDr7/uhx+07nPgwAGcP3/+hN1uP+ruOdWql2Sz2XYdPnz4wpFaV42+9vLgg/z3y1P5gy4yMoD//Ic38vBLxozhH9zVPri65Odz1x03rqLigMCnOBwOvPDCC+bs7OyJ1TmvWqZERGQwGMY988wzZrtdrNP1BT178j/wns7I+OQTnhc1YYJnr+sxOnTg3XNXrqzZkhCtls+yieUkkvHZZ59ZLRbLOiKqtCLA7ah2ZUkiOmexWP7fBx984IVapYI/olYDjRrxkc3dYjTywDlj3JS6d7/7sI3XYAwYP55H+WvyCOdw8GoAoiKAJFy5cgVvv/12jslkmlXdc2tU7tZkMs3++OOPDWfP3r4VuMCzKBR3X6760CGeoJmSwtN/3nmHh1v8umFNSAjf1yQtYMcOXl3gL3/xrCbBHSEiPPnkk3kWi2U8EVW7SmSNfiSJyGo2m0eNHDnSYhULG71Onz7cVGrKlSs85ykmBjh8GMjO5lncrjwlv6c6pmSx8Bm3oUP5B370Ue/pEtyWjz76qOz8+fPflZeX767J+TX+O2mz2X7Lycn5cNq0af6WelfniI/ns281HS1t2sSrEHz/PdCjh0eleRe5nO+r8wj22ms8OD5jBh8WiuYAPiUlJQWLFi26ajQaJ9f0Gnc1eDebzf/6/vvvk9esWeOPCxXqDBERfJ9Xw2SMqCi+z872jB6f4Uqgqk7pBFfPqHff5a2aBD7DbDZj5MiROSaT6REiqvEj1F2ZEhGR0WgcOX369KsprtWdAo/jWj9a01wlV2nqWjcR1a0bD6i9+657x+/cyWcE5szxwzosdRubzYZHH300Ly8v7wUiunA317rrMCcRFZpMpoHDhg0zXfXV8vZ6husppqb/vAsWAHo9n22rVcTG8gLmW7cC5jusBz91igfOWrUCXnnFN/oEAHhge8KECYUXLlz4uLCwcPPdXs8jcy9EdDEnJ2fUgAEDLHk1fcYQVIorDlSTtICrV3mQfNq0isfAWsWAAXxfVRvhnBy+vi0sjCd0hYf7QpngBm+88UZxYmLi9tzcXLcW3N4Jj00IW63Wn7Oysp4fMGBAnugV51liY3lLppq0QJs7l++HD/esJp9x3308WL1t2+3fv3aNT0+eP8+7MOh0vtVXz1m2bFnZihUrjhiNxrFEHgoQVNalsqZbRETEcz169LAUFxd7sqlmveeVV4hkMqKCAvfPSU/nTWanTvWaLO9z5Qr/EH/saltaSvSf/xDFxRGFhBDt2yeNvnrM559/btXpdL8AUJMHPcTjpkREiIyMnNizZ888YUyeY/Vq/r91+XLVxz32GNHIkUTJyURDhxKpVETnzvlGo1fIyeEfHCCaNImoc2eili25QwNEjRsTnTwptcp6x8qVK63R0dHJAALIw/7ByEtTMpGRkc81adLkvd27d0eEhYV55R71iQ0beNeRtDSgbdvKj/tjnuErrwDvveddbV7F4agoh6tU8kqUej1vv9KuHe817mrFJPAJH330UdnChQt/M5lMDxJRiaev77UuXxaLZVVoaGhhr169PktKSoqMiYnx1q3qBa6qG0V36CtTWAg0bsxzmr77DhgxwuvSvItczlcOnzvHF+r5VWO6+gUR4bXXXitetWrVYZPJNISq6HJ7N3h15VNhYeGG9PT0J7p3754r1sndHa65gzsNbENCKpZ7abWeqyYrKXI5Hx4KQ5IMm82GJ598smDVqlXfmUymQd4yJMDLpgQApaWlSVlZWff3798/KzExUfRpqiENGvD9wUoaHmdmAn/7GzB2LLBxIx9U9O7tM3mCOozFYkG/fv3ykpKSPjKZTOOJyOHVG3o6SFXZBkCv0+lOLl68uNTpdHo06FYfcDqJ/vxnonbtbv9+YmJFPHjwYKKjR32rT1A3OX36NDVp0iQnJCRkJPnKK3x1I+LGpNbr9etGjBiRX1RU5NF/vPrAokX8f2zLFm5SLubPrzCkd96RTp+gbrF27VqbXq+/rFAo2pEPfcKn1XSIyGowGMbs27fv1c6dO5tPnz7ty9vXelxVOIYPv7UvXFlZxetevXyrSVD3sFqtePHFF69PmzbtgNFo7GCz2dzqQuIpJCnxZbFYPs/IyLj//vvvz1i+fLmVvJSWUNdo2xbYs4e/Tkqq+P7NS73EOlTB3XDu3Dl06tTJvGnTpneNRuMAIvJ5WzWv5Sm5dXPGAvV6/eoOHToMSkhICNeJJQJu0bcvcOAAr5HkKtSm1/P6ZuXldWTGTeBTiAiffvqp9a233jJYLJYRNpvtuFRaJDUlF4GBgUPDw8M/X7p0afSoUaPkUuvxd1ym07MnX2y7cycwaBB/dKtsdk4gqIwrV65g9OjRlosXL242mUwvkxcSIquDX1RoLikp2Zadnd1u8uTJOwYNGpSXXeuqkfmOm4swNmvG964cpvJy3m9RIHAHp9OJJUuWWLt27Xr1119//YvRaHxWakMC4NvZN3e2gICAITExMVc//fRTq91uv4u5g7qJ00k0Z07FbFtWFpHDQfT++/zrUaOkViioDZw8eZI6duxo1uv1nwIIJD/43XdtfjFSupmSkpIfDQZD6wULFqzu0KGDObmqOjr1EMZ4JxIXCQm8I0nPnvzrVq2k0SWoHRQUFGDSpEmFAwcOTEtLSxtgMBgmkz+Mjm5GalesagPQRq/XHxszZkxeVlZWDf8m1E3y8/nISKUiys0lGjeOf71rl9TKBP6I3W6n5cuXW+Pi4gxhYWEvAJCRH/yO326TXMAdBQJMo9GM1Ov1ma+//nqxSLqsYOtW+u9jnGt74w2pVQn8jcTERGrRokXujUe1UPKD3+uqNr+YfXMHxpgqMjJyulqtnjFnzpzwF198UaVUKqWWJTlz595aV79RI+DyZZEWIAB+/fVXTJ482ZyZmZlsMBheIiIP9Fn2PrXGlFwwxkJ1Ot28gICA8W+++WbE2LFjFYp6vnr81CmgfXv+OiOD94kT1F9SU1Mxc+ZMS2pq6u/Z2dmTiOi41JqqQ60zJReMsSidTvemRqMZ/cYbb4SPGzdOUZ9HTnl5vOFH8+ZSKxFIxYkTJ/Dqq69aTp48mZ6TkzPNbrfXyqy1WmtKLhhjUXq9/g2lUjl65syZYc8//7wqMDBQalkCgU8gIvz888+YPXu2OT09Pd1kMk2z2+130eRdemq9KblgjIVptdppSqVy0oQJE0KmTJkSqNfrpZYlEHgFh8OBzZs304IFC8y5ubnHsrKyZhHRCal1eYI6Y0ouGGOa4ODgcUFBQbP79u0bPnfu3IiOHTtKLUsg8AgFBQVYsWKFdenSpYV2u32bwWB4u7YEsN2lzpmSC8YYk8lk98fExLyl1WpbvvrqqxGPP/64XK1WSy1NIKg2KSkp+OCDD/J3795dVF5evjQvL285EdWwkbt/U2dN6WYYY/dER0dPl8vlfx01alTASy+9FNq2qpYgAoEfUFhYiLVr19qXLFmSn5+ff9pgMLzldDp3E1GdLitdL0zJBWNMqVAoHtXpdK+Eh4c3f/nllyPGjBmjCBdtngV+gtPpxP79+7FkyRLLoUOHiu12+xdms/lzIsqSWpuvqFemdDOMsdjIyMgXlErl39q1axc8efLkqMGDB0M83gmkIDU1FatXr76+cePGMgA/Z2VlfQjgENXDX9B6a0o3wxjroNfrJwIY2bt3b+WECRO0AwcOhEqUcRR4kTNnzmDNmjUla9asKS4rKztrMpk+sdvtPxJR2Z3PrrsIU7oJxpgMQPfY2NjnnE7n4F69eimffvpp7cCBAxEUFCS1PEEth4hw/PhxbNy4sWTdunXF5eXlv+fm5q6wWq1b6mrQuiYIU6qEGwbVVa/Xj2WMDWvWrFnAuHHjIocOHSqPi4uTWp6gllBWVob9+/dj7dq1eYmJiTbG2HGj0filzWbbLozo9ghTchPGWPPw8PDHAwMDxwQEBMQOGTJEM3LkyJBevXqJxzzBLVy4cAHbt2+3rV+/Pj89Pd0qk8l2Z2VlfQ3gABHZpNbn7whTqgGMsWCZTNY/Njb2rw6Ho198fLxqxIgRoQ8++KC6U6dOkMtFmfH6RHZ2Nnbv3k1bt261HDp0yElEFwsKCtYWFRX9QETpUuurbQhT8gCMsXs0Gs2D0dHRo2w2259atmwpf+SRR0L79eun7tKlixhJ1SGICJcuXcKBAwdo165dloMHDzrLy8uzrVbrVrPZ/COAY2I0dHcIU/IwjDEGoJlKpeqj0+ketdvtXaOjozUDBgzQ9O3bN6Rbt25o0KCB1DIFblJSUoKUlBQcOXKkPDExsSA1NZXkcnlGcXHx9ry8vCQAvxKRVWqddQlhSj6AMaYF0DM6Orq/SqXq73Q6G7Vo0QL3339/SNeuXQM6deqEuLg4MFGZTVKKi4uRmpqK48ePO/fu3ZuXnJzsLCkpua5QKJJNJlNieXn5YQBn63pGtdQIU5KAG6OpJoyxrnq9vo9CoejucDgaxsTEyLt166bs0aNHeJs2bdCmTRuEhoZKLbfOYbfbcfHiRZw9exa//fZb6aFDh4rOnTtHNputUC6XnzCbzftLS0t/AXCivucMSYEwJT+CMRYFoGNISEiX8PDw7g6Hoy1jLEKn07GOHTsqOnfuHNaiRQt506ZNER8fL7LPq4CIYDAYkJ6ejvT0dKSmpl5PSUkpvXDhAisvLy9VKBTpVqs1JScn55cblRkvEpFDat0CYUq1ghtm1VqhULSJiorqqFQq29pstsZyuTwoIiJC1rJlS9a6deuAJk2aBDVo0ACuLTIyss4+EpaWliIrKwvXrl3DtWvXkJmZWX7+/Pnis2fP2q5evcqsVqtVoVAYAfyen59/sqioKA3AGQAZwnz8G2FKtRzGWDiAJgAaazSaeyIiIloqlcqmDoejocPhiFAqlarQ0FAWExND0dHR8piYGGVMTIxGp9OptFotbt5CQkKgUqkgk/muHSARwW63o6SkBBaLBWazGWazGRaLBbm5uU6DwVBiMBjKjUajIycnh0wmEysvL7cTUbFCocgmostFRUXnCwoKLhLRVQAXAWSJuE/tRZhSPYAxFgIgGkAkAC0AbUBAgD4kJKSBSqWKlclkeofDoXU6nUFEpAIg4+WoZOymPZPJZDKFQgGNRuNUq9XQaDSujQUEBDCbzUalpaUoKyujsrIyWK1WlJWVsfLycuZ0OomInDf2rhfOGwtO7TKZrEwul1sYYzl2u91QXFycVVhYmEVEZgA3b7nCcOo2wpQEbnMjQK8EoAagubG/ebMDsN7Yym56bRVGInAXYUoCgcCv8F3wQCAQCNxAmJJAIPArhCkJBAK/QpiSQCDwK4QpCQQCv+L/A8yth1nCo9ATAAAAAElFTkSuQmCC",
       "text/plain": [
        "<Figure size 360x360 with 1 Axes>"
       ]
@@ -1307,20 +1307,26 @@
    "source": [
     "# retrieve continental outline\n",
     "# This is the version that uses cartopy and requires installation of cartopy\n",
-    "af=continents.get_continent('af').transpose()\n",
-    "sam=continents.get_continent('sam').transpose()\n",
+    "af = continents.get_continent('af').transpose()\n",
+    "sam = continents.get_continent('sam').transpose()\n",
     "\n",
+    "# define options for pmagplotlib.plot_map\n",
+    "plt.figure(1, (5, 5))\n",
+    "Opts = {\n",
+    "    \"latmin\": -90, \"latmax\": 90, \"lonmin\": 0., \"lonmax\": 360.,\n",
+    "    \"lat_0\": -20, \"lon_0\": 345, \"proj\": \"ortho\",\n",
+    "    \"sym\": \"r-\", \"symsize\": 3, \"pltgrid\": 0, \"res\": \"c\", \"boundinglat\": 0.,\n",
+    "}\n",
     "\n",
-    "#define options for pmagplotlib.plot_map\n",
-    "plt.figure(1,(5,5))\n",
-    "Opts = {'latmin': -90, 'latmax': 90, 'lonmin': 0., 'lonmax': 360., 'lat_0': -20, \\\n",
-    "            'lon_0': 345,'proj': 'ortho', 'sym': 'r-', 'symsize': 3,\\\n",
-    "            'pltgrid': 0, 'res': 'c', 'boundinglat': 0.}\n",
     "if has_cartopy:\n",
-    "    pmagplotlib.plot_map(1,af[0],af[1],Opts)\n",
-    "    Opts['sym']='b-'\n",
-    "    pmagplotlib.plot_map(1,sam[0],sam[1],Opts)\n",
-    "    "
+    "    # draw Africa in red and return the axes handle\n",
+    "    ax = pmagplotlib.plot_map(1, af[0], af[1], Opts)\n",
+    "    # overlay South America in blue using the same axes\n",
+    "    Opts[\"sym\"] = \"b-\"\n",
+    "    Opts[\"ax\"] = ax\n",
+    "    pmagplotlib.plot_map(1, sam[0], sam[1], Opts)\n",
+    "\n",
+    "plt.show()"
    ]
   },
   {
@@ -11862,7 +11868,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.9.13"
+   "version": "3.9.12"
   }
  },
  "nbformat": 4,

{pmagpy-4.3.8.dist-info → pmagpy-4.3.10.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: pmagpy
-Version: 4.3.8
+Version: 4.3.10
 Summary: Analysis tools for paleo/rock magnetic data
 Home-page: https://github.com/PmagPy/PmagPy
 Author: PmagPy team