From ddfffa4914dc861c72c08afd1b5e864c2fe324de Mon Sep 17 00:00:00 2001 From: "Kalyn R. Dorheim" Date: Tue, 25 Jun 2024 12:34:03 -0400 Subject: [PATCH] Addressing a handful of open github issues (#741) * Github Issue 718. Add papers to Hector manuscript lists, including V3.2.0 manuscript * github issue 727. rename FCH4 to RF_CH4 to match all other RF output * Add a new release checklist * github issue 734, add output streams * quite build message * github issue 740, add matilda refernces * github issue 735, add ref temperature period to plot --- .../ISSUE_TEMPLATE/new_release_checklist.md | 22 ++++ NEWS.md | 2 + README.Rmd | 3 +- inst/include/component_data.hpp | 2 +- inst/include/forcing_component.hpp | 2 +- inst/input/hector_ssp119.ini | 20 ++-- inst/input/hector_ssp126.ini | 14 +-- inst/input/hector_ssp245.ini | 20 ++-- inst/input/hector_ssp370.ini | 20 ++-- inst/input/hector_ssp434.ini | 20 ++-- inst/input/hector_ssp460.ini | 20 ++-- inst/input/hector_ssp534-over.ini | 20 ++-- inst/input/hector_ssp585.ini | 20 ++-- man/figures/README-tas_plot-1.png | Bin 56572 -> 57577 bytes output/README.md | 11 +- src/simpleNbox-runtime.cpp | 10 +- vignettes/articles/Hector-Wider-World.Rmd | 108 ++++++++---------- vignettes/articles/ex_GMAT-GMST.Rmd | 5 - vignettes/hector.Rmd | 8 +- 19 files changed, 175 insertions(+), 152 deletions(-) create mode 100644 .github/ISSUE_TEMPLATE/new_release_checklist.md diff --git a/.github/ISSUE_TEMPLATE/new_release_checklist.md b/.github/ISSUE_TEMPLATE/new_release_checklist.md new file mode 100644 index 000000000..dcb6e62dc --- /dev/null +++ b/.github/ISSUE_TEMPLATE/new_release_checklist.md @@ -0,0 +1,22 @@ +--- +name: Release +about: Checklist! +--- + +Since the dev branch only includes materials that have been peer-reviewed at least once this PR review process should be relatively smooth should. That being said, before merging this branch into main and the new release is created several things must be completed first. + +* Open a PR into main branch, it should be titled with the new release number +* Let Leeyabot generate a report detailing how the Hector output has changed! +* Add text to the PR documenting the major model changes. The PR is going to serve as an important record +* Pass automated checks +* Change the version number in the DESCRIPTION file +* Update the NEWS.md, are links additional materials necessary? Does the release need a new naming file? +* Check to make sure the internal package data is up to date (fxntable, inputstable, unitstable) + +After merging PR + +* Create the release via zenodo +* Add the zenodo links to README / NEWS & push directly to main +* Create new dev branch & update the DESCRIPTION version number to support the pkgdown documentation for the developmental branch! +* Add archived outputstreams to release assets! + diff --git a/NEWS.md b/NEWS.md index 101ca4ef9..718d4a85d 100644 --- a/NEWS.md +++ b/NEWS.md @@ -1,5 +1,7 @@ # hector 3.2.0 [![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.10698028.svg)](https://doi.org/10.5281/zenodo.10698028) + +* Fully documented in Dorheim, K., Gering, S., Gieseke, R., Hartin, C., Pressburger, L., Shiklomanov, A. N., Smith, S. J., Tebaldi, C., Woodard, D. L., and Bond-Lamberty, B.: Hector V3.2.0: functionality and performance of a reduced-complexity climate model, Geosci. Model Dev., 17, 4855–4869, https://doi.org/10.5194/gmd-17-4855-2024, 2024. * Correct aerosol forcing coefficients based on Zelinka et al. (2023) * Enable permafrost module and recalibrate model's default parameterization diff --git a/README.Rmd b/README.Rmd index 484777714..a3667f452 100644 --- a/README.Rmd +++ b/README.Rmd @@ -58,7 +58,8 @@ hector_tas_results$scenario <- ifelse(hector_tas_results$year <= 2016, "historic ggplot(hector_tas_results) + geom_line(aes(year, value, color = scenario), linewidth = 1) + theme_bw(base_size = 15) + - labs(color = NULL, x = NULL, y = expression("Temperature Anomaly ("~degree~"C)")) + + theme(axis.title = element_text(size = 12))+ + labs(color = NULL, x = NULL, y = expression("Temperature Anomaly Relative to 1750 ("~degree~"C)")) + scale_color_manual(values = SSP_COLORS) ``` diff --git a/inst/include/component_data.hpp b/inst/include/component_data.hpp index 84d850f8a..afa24bc02 100644 --- a/inst/include/component_data.hpp +++ b/inst/include/component_data.hpp @@ -256,7 +256,7 @@ // methane component #define D_CH4_CONC "CH4" CONCENTRATION_EXTENSION #define D_PREINDUSTRIAL_CH4 "M0" -#define D_RF_CH4 "FCH4" +#define D_RF_CH4 D_RF_PREFIX "CH4" #define D_EMISSIONS_CH4 "CH4_emissions" #define D_CONSTRAINT_CH4 "CH4_constrain" #define D_NATURAL_CH4 "CH4N" diff --git a/inst/include/forcing_component.hpp b/inst/include/forcing_component.hpp index 414ba4342..be920c26d 100644 --- a/inst/include/forcing_component.hpp +++ b/inst/include/forcing_component.hpp @@ -120,7 +120,7 @@ class ForcingComponent : public IModelComponent { // Aerosol parameters for aerosol-cloud interactions (RFaci) see equation // Equation 7.SM.1.2 of IPCC AR6 - double const aci_beta = 2.279759; // see Dorheim et al. in prep + double const aci_beta = 2.279759; // Dorheim et al. 2024 double const s_BCOC = 111.05064063; // (Tg C yr-1) IPCC AR6 7.SM.1.3.1 double const s_SO2 = (260.34644166 * 1000) * diff --git a/inst/input/hector_ssp119.ini b/inst/input/hector_ssp119.ini index 276fd795a..5f28e33a1 100644 --- a/inst/input/hector_ssp119.ini +++ b/inst/input/hector_ssp119.ini @@ -37,7 +37,7 @@ preind_interdeep_c=37100 ; intermediate and deep ocean C, IPCC AR6 Figure 5.12 atmos_co2=590.33 ; initial atmospheric carbon pool Pg C of CO2 size used in spin up, consistent with C0 C0=277.15 ; preindustrial CO2 from table 7.SM.1 IPCC AR6; CO2 radiative forcing is calibrated to this value npp_flux0=56.2 ; preindustrial net primary production, Pg C/yr from Ito et al. 2011 -;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 +;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 ;nonpf.npp_flux0=50.6 veg_c=550 ; initial vegetation pool, Pg C @@ -73,16 +73,16 @@ f_litterd=0.98 ; Fraction of litter flux to detritus (balance to soil) ;nonpf.f_litterd=0.98 ; Anthropogenic contributions: emissions, direct air carbon capture and storage, and land use change, Pg C/yr -ffi_emissions=csv:tables/ssp119_emiss-constraints_rf.csv +ffi_emissions=csv:tables/ssp119_emiss-constraints_rf.csv daccs_uptake=csv:tables/ssp119_emiss-constraints_rf.csv luc_emissions=csv:tables/ssp119_emiss-constraints_rf.csv luc_uptake=csv:tables/ssp119_emiss-constraints_rf.csv ; CO2 and temperature effects on the carbon cycle -beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.beta=0.53 ;nonpf.beta=0.53 -q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.q10_rh=1.76 ;nonpf.q10_rh=1.76 @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6) delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6) delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6) -rho_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep -rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep +rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024 +rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024 ; Miscellaneous radiative forcings are by default zero, but can read in from a input table to ; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails) @@ -178,7 +178,7 @@ RF_misc[1750]=0 ;------------------------------------------------------------------------ [temperature] S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6) -diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations see Dorheim et al. in prep for details +diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details alpha=1.0 ; uncertainty scaling factor for aerosol forcing volscl=1.0 ; uncertainty scaling factor for volcanic forcing qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6) @@ -232,7 +232,7 @@ delta_HFC32=0.0 ; tropospheric adjustments scalar unitless (7.3.2.4 IPCC AR6 molarMass=52.0 ; grams HFC32_emissions=csv:tables/ssp119_emiss-constraints_rf.csv -[HFC4310_halocarbon] +[HFC4310_halocarbon] ; This is equivalent to HFC-43-10mee tau=17.0 ; lifetime in years (Table 7.SM.7 IPPC AR6) rho_HFC4310=0.000357 ; radiative efficiencies W/m2/ppt (Table 7.SM.7 IPPC AR6) diff --git a/inst/input/hector_ssp126.ini b/inst/input/hector_ssp126.ini index c8dc92e2c..5482e89e9 100644 --- a/inst/input/hector_ssp126.ini +++ b/inst/input/hector_ssp126.ini @@ -79,10 +79,10 @@ luc_emissions=csv:tables/ssp126_emiss-constraints_rf.csv luc_uptake=csv:tables/ssp126_emiss-constraints_rf.csv ; CO2 and temperature effects on the carbon cycle -beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.beta=0.53 ;nonpf.beta=0.53 -q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.q10_rh=1.76 ;nonpf.q10_rh=1.76 @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6) delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6) delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6) -rho_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep -rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep +rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024 +rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024 ; Miscellaneous radiative forcings are by default zero, but can read in from a input table to ; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails) @@ -178,7 +178,7 @@ RF_misc[1750]=0 ;------------------------------------------------------------------------ [temperature] S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6) -diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations see Dorheim et al. in prep for details +diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details alpha=1.0 ; uncertainty scaling factor for aerosol forcing volscl=1.0 ; uncertainty scaling factor for volcanic forcing qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6) diff --git a/inst/input/hector_ssp245.ini b/inst/input/hector_ssp245.ini index b511b93aa..d1eaf2b37 100644 --- a/inst/input/hector_ssp245.ini +++ b/inst/input/hector_ssp245.ini @@ -37,7 +37,7 @@ preind_interdeep_c=37100 ; intermediate and deep ocean C, IPCC AR6 Figure 5.12 atmos_co2=590.33 ; initial atmospheric carbon pool Pg C of CO2 size used in spin up, consistent with C0 C0=277.15 ; preindustrial CO2 from table 7.SM.1 IPCC AR6; CO2 radiative forcing is calibrated to this value npp_flux0=56.2 ; preindustrial net primary production, Pg C/yr from Ito et al. 2011 -;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 +;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 ;nonpf.npp_flux0=50.6 veg_c=550 ; initial vegetation pool, Pg C @@ -73,16 +73,16 @@ f_litterd=0.98 ; Fraction of litter flux to detritus (balance to soil) ;nonpf.f_litterd=0.98 ; Anthropogenic contributions: emissions, direct air carbon capture and storage, and land use change, Pg C/yr -ffi_emissions=csv:tables/ssp245_emiss-constraints_rf.csv +ffi_emissions=csv:tables/ssp245_emiss-constraints_rf.csv daccs_uptake=csv:tables/ssp245_emiss-constraints_rf.csv luc_emissions=csv:tables/ssp245_emiss-constraints_rf.csv luc_uptake=csv:tables/ssp245_emiss-constraints_rf.csv ; CO2 and temperature effects on the carbon cycle -beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.beta=0.53 ;nonpf.beta=0.53 -q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.q10_rh=1.76 ;nonpf.q10_rh=1.76 @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6) delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6) delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6) -rho_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep -rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep +rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024 +rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024 ; Miscellaneous radiative forcings are by default zero, but can read in from a input table to ; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails) @@ -178,7 +178,7 @@ RF_misc[1750]=0 ;------------------------------------------------------------------------ [temperature] S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6) -diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations see Dorheim et al. in prep for details +diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details alpha=1.0 ; uncertainty scaling factor for aerosol forcing volscl=1.0 ; uncertainty scaling factor for volcanic forcing qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6) @@ -232,7 +232,7 @@ delta_HFC32=0.0 ; tropospheric adjustments scalar unitless (7.3.2.4 IPCC AR6 molarMass=52.0 ; grams HFC32_emissions=csv:tables/ssp245_emiss-constraints_rf.csv -[HFC4310_halocarbon] +[HFC4310_halocarbon] ; This is equivalent to HFC-43-10mee tau=17.0 ; lifetime in years (Table 7.SM.7 IPPC AR6) rho_HFC4310=0.000357 ; radiative efficiencies W/m2/ppt (Table 7.SM.7 IPPC AR6) diff --git a/inst/input/hector_ssp370.ini b/inst/input/hector_ssp370.ini index 2cbb90dc9..58a13243b 100644 --- a/inst/input/hector_ssp370.ini +++ b/inst/input/hector_ssp370.ini @@ -37,7 +37,7 @@ preind_interdeep_c=37100 ; intermediate and deep ocean C, IPCC AR6 Figure 5.12 atmos_co2=590.33 ; initial atmospheric carbon pool Pg C of CO2 size used in spin up, consistent with C0 C0=277.15 ; preindustrial CO2 from table 7.SM.1 IPCC AR6; CO2 radiative forcing is calibrated to this value npp_flux0=56.2 ; preindustrial net primary production, Pg C/yr from Ito et al. 2011 -;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 +;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 ;nonpf.npp_flux0=50.6 veg_c=550 ; initial vegetation pool, Pg C @@ -73,16 +73,16 @@ f_litterd=0.98 ; Fraction of litter flux to detritus (balance to soil) ;nonpf.f_litterd=0.98 ; Anthropogenic contributions: emissions, direct air carbon capture and storage, and land use change, Pg C/yr -ffi_emissions=csv:tables/ssp370_emiss-constraints_rf.csv +ffi_emissions=csv:tables/ssp370_emiss-constraints_rf.csv daccs_uptake=csv:tables/ssp370_emiss-constraints_rf.csv luc_emissions=csv:tables/ssp370_emiss-constraints_rf.csv luc_uptake=csv:tables/ssp370_emiss-constraints_rf.csv ; CO2 and temperature effects on the carbon cycle -beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.beta=0.53 ;nonpf.beta=0.53 -q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.q10_rh=1.76 ;nonpf.q10_rh=1.76 @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6) delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6) delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6) -rho_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep -rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep +rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024 +rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024 ; Miscellaneous radiative forcings are by default zero, but can read in from a input table to ; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails) @@ -178,7 +178,7 @@ RF_misc[1750]=0 ;------------------------------------------------------------------------ [temperature] S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6) -diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations see Dorheim et al. in prep for details +diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details alpha=1.0 ; uncertainty scaling factor for aerosol forcing volscl=1.0 ; uncertainty scaling factor for volcanic forcing qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6) @@ -232,7 +232,7 @@ delta_HFC32=0.0 ; tropospheric adjustments scalar unitless (7.3.2.4 IPCC AR6 molarMass=52.0 ; grams HFC32_emissions=csv:tables/ssp370_emiss-constraints_rf.csv -[HFC4310_halocarbon] +[HFC4310_halocarbon] ; This is equivalent to HFC-43-10mee tau=17.0 ; lifetime in years (Table 7.SM.7 IPPC AR6) rho_HFC4310=0.000357 ; radiative efficiencies W/m2/ppt (Table 7.SM.7 IPPC AR6) diff --git a/inst/input/hector_ssp434.ini b/inst/input/hector_ssp434.ini index d2865e09f..217ab12ad 100644 --- a/inst/input/hector_ssp434.ini +++ b/inst/input/hector_ssp434.ini @@ -37,7 +37,7 @@ preind_interdeep_c=37100 ; intermediate and deep ocean C, IPCC AR6 Figure 5.12 atmos_co2=590.33 ; initial atmospheric carbon pool Pg C of CO2 size used in spin up, consistent with C0 C0=277.15 ; preindustrial CO2 from table 7.SM.1 IPCC AR6; CO2 radiative forcing is calibrated to this value npp_flux0=56.2 ; preindustrial net primary production, Pg C/yr from Ito et al. 2011 -;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 +;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 ;nonpf.npp_flux0=50.6 veg_c=550 ; initial vegetation pool, Pg C @@ -73,16 +73,16 @@ f_litterd=0.98 ; Fraction of litter flux to detritus (balance to soil) ;nonpf.f_litterd=0.98 ; Anthropogenic contributions: emissions, direct air carbon capture and storage, and land use change, Pg C/yr -ffi_emissions=csv:tables/ssp434_emiss-constraints_rf.csv +ffi_emissions=csv:tables/ssp434_emiss-constraints_rf.csv daccs_uptake=csv:tables/ssp434_emiss-constraints_rf.csv luc_emissions=csv:tables/ssp434_emiss-constraints_rf.csv luc_uptake=csv:tables/ssp434_emiss-constraints_rf.csv ; CO2 and temperature effects on the carbon cycle -beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.beta=0.53 ;nonpf.beta=0.53 -q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.q10_rh=1.76 ;nonpf.q10_rh=1.76 @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6) delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6) delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6) -rho_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep -rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep +rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024 +rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024 ; Miscellaneous radiative forcings are by default zero, but can read in from a input table to ; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails) @@ -178,7 +178,7 @@ RF_misc[1750]=0 ;------------------------------------------------------------------------ [temperature] S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6) -diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations see Dorheim et al. in prep for details +diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details alpha=1.0 ; uncertainty scaling factor for aerosol forcing volscl=1.0 ; uncertainty scaling factor for volcanic forcing qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6) @@ -232,7 +232,7 @@ delta_HFC32=0.0 ; tropospheric adjustments scalar unitless (7.3.2.4 IPCC AR6 molarMass=52.0 ; grams HFC32_emissions=csv:tables/ssp434_emiss-constraints_rf.csv -[HFC4310_halocarbon] +[HFC4310_halocarbon] ; This is equivalent to HFC-43-10mee tau=17.0 ; lifetime in years (Table 7.SM.7 IPPC AR6) rho_HFC4310=0.000357 ; radiative efficiencies W/m2/ppt (Table 7.SM.7 IPPC AR6) diff --git a/inst/input/hector_ssp460.ini b/inst/input/hector_ssp460.ini index c1ba357db..50f4e2c7d 100644 --- a/inst/input/hector_ssp460.ini +++ b/inst/input/hector_ssp460.ini @@ -37,7 +37,7 @@ preind_interdeep_c=37100 ; intermediate and deep ocean C, IPCC AR6 Figure 5.12 atmos_co2=590.33 ; initial atmospheric carbon pool Pg C of CO2 size used in spin up, consistent with C0 C0=277.15 ; preindustrial CO2 from table 7.SM.1 IPCC AR6; CO2 radiative forcing is calibrated to this value npp_flux0=56.2 ; preindustrial net primary production, Pg C/yr from Ito et al. 2011 -;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 +;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 ;nonpf.npp_flux0=50.6 veg_c=550 ; initial vegetation pool, Pg C @@ -73,16 +73,16 @@ f_litterd=0.98 ; Fraction of litter flux to detritus (balance to soil) ;nonpf.f_litterd=0.98 ; Anthropogenic contributions: emissions, direct air carbon capture and storage, and land use change, Pg C/yr -ffi_emissions=csv:tables/ssp460_emiss-constraints_rf.csv +ffi_emissions=csv:tables/ssp460_emiss-constraints_rf.csv daccs_uptake=csv:tables/ssp460_emiss-constraints_rf.csv luc_emissions=csv:tables/ssp460_emiss-constraints_rf.csv luc_uptake=csv:tables/ssp460_emiss-constraints_rf.csv ; CO2 and temperature effects on the carbon cycle -beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.beta=0.53 ;nonpf.beta=0.53 -q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.q10_rh=1.76 ;nonpf.q10_rh=1.76 @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6) delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6) delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6) -rho_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep -rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep +rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024 +rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024 ; Miscellaneous radiative forcings are by default zero, but can read in from a input table to ; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails) @@ -178,7 +178,7 @@ RF_misc[1750]=0 ;------------------------------------------------------------------------ [temperature] S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6) -diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations see Dorheim et al. in prep for details +diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details alpha=1.0 ; uncertainty scaling factor for aerosol forcing volscl=1.0 ; uncertainty scaling factor for volcanic forcing qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6) @@ -232,7 +232,7 @@ delta_HFC32=0.0 ; tropospheric adjustments scalar unitless (7.3.2.4 IPCC AR6 molarMass=52.0 ; grams HFC32_emissions=csv:tables/ssp460_emiss-constraints_rf.csv -[HFC4310_halocarbon] +[HFC4310_halocarbon] ; This is equivalent to HFC-43-10mee tau=17.0 ; lifetime in years (Table 7.SM.7 IPPC AR6) rho_HFC4310=0.000357 ; radiative efficiencies W/m2/ppt (Table 7.SM.7 IPPC AR6) diff --git a/inst/input/hector_ssp534-over.ini b/inst/input/hector_ssp534-over.ini index d5a76b1d5..8c97b0d23 100644 --- a/inst/input/hector_ssp534-over.ini +++ b/inst/input/hector_ssp534-over.ini @@ -37,7 +37,7 @@ preind_interdeep_c=37100 ; intermediate and deep ocean C, IPCC AR6 Figure 5.12 atmos_co2=590.33 ; initial atmospheric carbon pool Pg C of CO2 size used in spin up, consistent with C0 C0=277.15 ; preindustrial CO2 from table 7.SM.1 IPCC AR6; CO2 radiative forcing is calibrated to this value npp_flux0=56.2 ; preindustrial net primary production, Pg C/yr from Ito et al. 2011 -;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 +;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 ;nonpf.npp_flux0=50.6 veg_c=550 ; initial vegetation pool, Pg C @@ -73,16 +73,16 @@ f_litterd=0.98 ; Fraction of litter flux to detritus (balance to soil) ;nonpf.f_litterd=0.98 ; Anthropogenic contributions: emissions, direct air carbon capture and storage, and land use change, Pg C/yr -ffi_emissions=csv:tables/ssp534-over_emiss-constraints_rf.csv +ffi_emissions=csv:tables/ssp534-over_emiss-constraints_rf.csv daccs_uptake=csv:tables/ssp534-over_emiss-constraints_rf.csv luc_emissions=csv:tables/ssp534-over_emiss-constraints_rf.csv luc_uptake=csv:tables/ssp534-over_emiss-constraints_rf.csv ; CO2 and temperature effects on the carbon cycle -beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.beta=0.53 ;nonpf.beta=0.53 -q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.q10_rh=1.76 ;nonpf.q10_rh=1.76 @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6) delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6) delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6) -rho_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep -rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep +rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024 +rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024 ; Miscellaneous radiative forcings are by default zero, but can read in from a input table to ; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails) @@ -178,7 +178,7 @@ RF_misc[1750]=0 ;------------------------------------------------------------------------ [temperature] S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6) -diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations see Dorheim et al. in prep for details +diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details alpha=1.0 ; uncertainty scaling factor for aerosol forcing volscl=1.0 ; uncertainty scaling factor for volcanic forcing qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6) @@ -232,7 +232,7 @@ delta_HFC32=0.0 ; tropospheric adjustments scalar unitless (7.3.2.4 IPCC AR6 molarMass=52.0 ; grams HFC32_emissions=csv:tables/ssp534-over_emiss-constraints_rf.csv -[HFC4310_halocarbon] +[HFC4310_halocarbon] ; This is equivalent to HFC-43-10mee tau=17.0 ; lifetime in years (Table 7.SM.7 IPPC AR6) rho_HFC4310=0.000357 ; radiative efficiencies W/m2/ppt (Table 7.SM.7 IPPC AR6) diff --git a/inst/input/hector_ssp585.ini b/inst/input/hector_ssp585.ini index c2a011186..d41fbe7d6 100644 --- a/inst/input/hector_ssp585.ini +++ b/inst/input/hector_ssp585.ini @@ -37,7 +37,7 @@ preind_interdeep_c=37100 ; intermediate and deep ocean C, IPCC AR6 Figure 5.12 atmos_co2=590.33 ; initial atmospheric carbon pool Pg C of CO2 size used in spin up, consistent with C0 C0=277.15 ; preindustrial CO2 from table 7.SM.1 IPCC AR6; CO2 radiative forcing is calibrated to this value npp_flux0=56.2 ; preindustrial net primary production, Pg C/yr from Ito et al. 2011 -;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 +;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978 ;nonpf.npp_flux0=50.6 veg_c=550 ; initial vegetation pool, Pg C @@ -73,16 +73,16 @@ f_litterd=0.98 ; Fraction of litter flux to detritus (balance to soil) ;nonpf.f_litterd=0.98 ; Anthropogenic contributions: emissions, direct air carbon capture and storage, and land use change, Pg C/yr -ffi_emissions=csv:tables/ssp585_emiss-constraints_rf.csv +ffi_emissions=csv:tables/ssp585_emiss-constraints_rf.csv daccs_uptake=csv:tables/ssp585_emiss-constraints_rf.csv luc_emissions=csv:tables/ssp585_emiss-constraints_rf.csv luc_uptake=csv:tables/ssp585_emiss-constraints_rf.csv ; CO2 and temperature effects on the carbon cycle -beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.beta=0.53 ;nonpf.beta=0.53 -q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details +q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details ;pf.q10_rh=1.76 ;nonpf.q10_rh=1.76 @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6) delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6) delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6) -rho_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep -rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep -rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep +rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024 +rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024 +rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024 ; Miscellaneous radiative forcings are by default zero, but can read in from a input table to ; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails) @@ -178,7 +178,7 @@ RF_misc[1750]=0 ;------------------------------------------------------------------------ [temperature] S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6) -diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations see Dorheim et al. in prep for details +diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details alpha=1.0 ; uncertainty scaling factor for aerosol forcing volscl=1.0 ; uncertainty scaling factor for volcanic forcing qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6) @@ -232,7 +232,7 @@ delta_HFC32=0.0 ; tropospheric adjustments scalar unitless (7.3.2.4 IPCC AR6 molarMass=52.0 ; grams HFC32_emissions=csv:tables/ssp585_emiss-constraints_rf.csv -[HFC4310_halocarbon] +[HFC4310_halocarbon] ; This is equivalent to HFC-43-10mee tau=17.0 ; lifetime in years (Table 7.SM.7 IPPC AR6) rho_HFC4310=0.000357 ; radiative efficiencies W/m2/ppt (Table 7.SM.7 IPPC AR6) diff --git a/man/figures/README-tas_plot-1.png b/man/figures/README-tas_plot-1.png index 8da4207b1b3c62f425830f8b022f478e504fff49..f0428ed4c5ffd29c785c0fd9f02c8cf94fa986e3 100644 GIT binary patch literal 57577 zcmY&g1z40@*Bu5>KsqG^gh5&mln?}@yBjH`JEUVs0qF)Q=?>{`>5%S_mhOi6-?>-4 z-|zGKFwZda#yR`!z4lt`5F{@vhJi+c1_FUFB*cXkK_EonC4v{q1K|LnXNtC{d>hfbyO3h~jD7Vs%aooE|1EY~-g5LZOskFdsafG(Fmpzw zdE=S}_FO z6MQe%_Ck7l(}Ozo{t~s5$0C_>1-GLLr1m%!>#%${F_JA%SA>bKAmqMtF7V+d^|>a} zL;J(IFy2>r#q?DDBE1}iUu3DrntX4d)glY#opZi&;T5bN_)Cbj=kW}7aX$mF(^S^x zKL(aQH?8=N4Ku7WMLnVhX*V?b+o2utb$(EBYb4wHm|ws-(#n`So@k7{Fv#@pAm zFMM0|LN^)B!Ba6AeFKXew7kN=T;_T~O-%4M^i-az2oa7NC@o}X!CQ~C1BOV z6ON$7GIiI|%gjY1*)mDr$~H;AoLJsnYNx-t2h}K>-S_Nljj#&lz^D&@6to|hs8Bbe zd%BnuI{ZL2vY{S-A;BUM<$Fm+6}`qv{4|p?indqDAjDwgwIt@N?l>xYgOj?EaTf$? z#Wc>~#^5cdw_W$f1_{fR!`~&-II*XPr-#k^VyCQhUVB)n1viMQrZr63|N00qCNX3F z*fTqxVmAzH*WH{d)7|9MM8ALNi3r`9Q>}Tm^Kw8H&7z*~cDl9{ z_Gw~f<1Pm;&G$<)zoBiI#w*OH32II=^U@jkziIDiH;|IsBKYKos?Xx;$XSRNN~GQJ zpu5c4pwOepizwfI*Q%&u@+ZHtXxsx@_ng%Oc0^gLRX=J8eqt(r%LZ4!y|WxgPpLPv zrQTikM;5%u@7;CXU-IGex7+v>iEV{VMOrTnRBvsjT|4{kfxeZS_4`DH{VJ z`g-eEx}Ey12C|VH%&K^kbx@fb8&|uv#Qh!7oJ+Q$ecA)ATUwqA&1!C}aPW_zwbszk>g~g9!SH_~$i( z7yOrp*sJeBAOVns@LOePgxwV61Z5S{?xj7$U+VQd*|CIKq^Q5b@s;HxVz803l;j13 zac~~V%hxLjWMv7}lM+_1Ho(~D^&MvtW>RKS`tO%E6H|CNWp6k3Q%+J(K1umK6MTw* zOZh*CI3d!zXc3kfiU0R5C4z_-8ps;y|GgQWOrtMSOl~ zTU#40Ir%-CLy*Vc=WrHn{r+8~v9Yno7fGDb8|j@+j&v)tW@SWwJWFs=k{0-za5GtN z)RCRTt3@Bw+0chMk;0pyj-Dr<#BC{z+FkBjz(WB zNKqp|{~cW1z|6_sBY$&&2!Q`P^yI+Gpk0Y*!^r>nY=l=d-dFJQnl5&d*#CT10DKM+ zKyFQasX_SPd4X#Hr_K!m<6~n}zHvKO{>j0d82emX^AqM1v%289|1R~4q!Zh#?LSc8gEAKCUjdj~d$Yltw8~?w+u~7xf zt|Hp(&vVnOmwhF<-I>8kAt0F(`G5B;N`qZW;5Aoc-G6^~vp)UpgWc_=ZJpaii9v6a z$Y(-khEP(@xJ!fW4m7W%q$KGC*0M1hZ= z#kNG%pYsMe5%N_0?SA0+#<@kIp~lajKYgm_k^&5$giVV+_|P+!D~GJ16bg~o3GLEq zX=&43;y|;XU&W%}3ppqLS>%E8=#i34$4Ot{1zNc%=^i)AvZx2)ic6r|ce)I@W#S_Z z6Z|d2`0O{NQ86*&c4sQ2O5W9V_V-H-C390xmFlXAsN~7N zWfdq@E&fs^BIWx4V;qr$4xckmMc}PS(PP|#G%;WXbA*n z$ha#VO0Ha5oMC@lubQq0&DqXWJaEQtZp!X~flnBuNQln;6TZzJHw-j-gI``r_Lddx z{zw-x6J4}hilJ9h_5i}5W+6;-E@ zhp&3k)`AhH2b;#jsj)4vdo6blj}JvSSW}+~7z?Hje|8w|&(|H#7^n44=|Z2ZzT)zF z>yRlLs>AiLuPc&TY^^Vrv5&v?*st8MU(-{0PTir^7Z;J(4b8LYDcxL0M@Om%PbCwb z)phGoI>Ahx4U^M{2&BpL`{#S}7!_gCv@r79})c0Dq+!6lT@nUDi_etiq7o8o(`n$M?M3&Q7Hv6jH zxUEMLPFPEoI!#D=*;hpr19|%g>U0r^>0Pb-{z5~pZs{Tsiw;wh+l6XaYb2)nRl0yT zn8D}AXttE$QnLpGTUUlia7?el&H0{g8LXyUJyrD2YK}lbA~(;e*9r(_olpeAo=W zI$n={>-#{W!e|i405m>Z{MrqURR!p&!hbLH5^1t|PE9a1;dcG_r5T+I=_~OspUPB$ z+*636zK-Y>#rl}vcEyiv%w;u{lZt@PZcRL2A@frs?EYr7&ffT~B8yfE6#E)awlR?H zXlsO;;QhEsm`CK-L&V>C6M_ithTe+|73&nbZf9S?RLJ?C5ubk1h9G>vD>Fl0Uilte z?s8&8%xRU}c(suk8b0xz9I|75tW8M3;O3_NIJM+ReDg_V8DFchilsK{DXulEG#TF)TZ~Z18yY#ginav z@@Vy*-P+uNve^PQ>E_(XPNNdW;bXA{`z^i+kh z903iTVPA}>ar=4npW-SkL`g*@IN3NR*6u4bK2qM~Dmlauo5>yB0IbC1M6U@n&}#(R zB^+weJT()Oykvf;&hc7rhji1>Ek=dq@8zaP3l<~vq*Zg=Rm{#F_#(2Hz_ln#Q}GyB z2%QkUXt3T4YY_^5PDeNmuGQ}=w@OL1!uGXi*bEZtpbo~;1P32~4g530dSqZ#9t6@{ zDGU0))6s-#LGn@=yuYCKeFG>il9lG>2Mp@us1X8GaHuE#RRiPqXCKG5X7BN=ZDo8f zxZ0<0%x+7b5#2+ZlL14%0-tQ1`x7GKl%in$r%%hz7n(O#Rm9MvPp?jr`0W7**XoR9 zp$L*qd7a(v-j0wU_-7Fw9^*~(7=D@OaglU*`C)`g(dDo;zGZCpDQe>yeCVIar7~3* z5678Ll}IeqJNQesN_iTWx&1pty{WB2%u~Z9!>MA|rd$_~Y)?+5&)=bVA%)1c7SIU< z-UJKll5$wc--&Sjm?%_}0t#H{BD6Z~4}PV50Nzt-ZtqSdO{;NIQ4-mx&esz z5aHGW0uT`q(b3V-Kftf~>4|l^pnu4?T=if0+xh^Ut5y_aUpSe^*f^RU`&x;1#^#kA zfTIC)rdkZN0pg%ZAW^tnpXN6nb&)dMl(IX8SW^GN*?|<^*ih)q4_I9>wD>4ZqHiE^ zRqb&p74ZGnU{yvC0i$XR|NW(q9;P1w0@|NxYwec+Il{tOv{h63RlhaH?A)FA&nTRK zMv9_GD3tbMnz1>wmhk-p$PVprY70ef=Xf0I^OZ`TYYxCzIz>+I;8~%O7wMAX?$|wk z7xnkKAqlP5&36#qPW*jRcGwJKjGRp(6EL@}btv)(HYEqX@H`L4`rZ&~o0X|-JKGag zm6mYpXZhS#7Ex9H1V9*ag?UZzfx^svJfuWD!ET#$UQ38b>m6?`Gsf1le%i*n9t>#+E9boW~dYO(PU7@Ac~9yewzoz3o7gRRneQ4@dMLJ?In)B zCLcG*i-YK8S~-K&0N+7pohdTt8^zq}prm^;;isajmit_s;*`YL&iVKtUOEbXPbMsK z-m)FFT$tX?)fKTM!VhTC-e5#dtAdPze#Cx1K*|EIW;#&W6Hd#2%4B!v|rSDyu6S3^E(Tkpo94H$j`r zgGrpRidF|p@5R{lt{@S~DC$*Ksc+mIf&3RXUZkNE_guT+f$_Juv-J*z@I*_(X~m>b zZV{{g+ZmBdJmtZqCu%K{qJ7KvGsv3-$y#xs8b2`Z%xz!UU@wWNhg2$(j; z0I`Y!FD0j8l!%A?u-nWB=RA(P5g)q4dI8L~4opbjZn}6R6`RF0tLa>Iu_F5E&ZX?1 zWVGy#6mjIV+WK_k)5I4p*+vYbFGg0poIuVldIlt1FGU>WTmTqvXHS81IjCorxA!)9 z{wy71a6^|2b_=5#fzNKKO3dg&{WD{-T^>{jDOmLi)G6ZY3#&yr;`;q@E+X%n&7t0A zQ+DZ4!LB71V^W5UIQtO;bMaPnq}>eAmHGYsR#eM#$uQ}s@R&%4t+WdCacBa$yN6oIK7uQdlffWg$5M@wQTl~QDsDJB|mn@~U z)br}2jCaOm;!LxBto7q{DbRdZWJ%h4;);P-l_~=y;W8qG zLIM!_HLYxN9vecpGIu#b%k_34yf36z$Pkv!wdMaa!U9xywe?PEWH)V^Y*kJV_*?3L@84kx5`-6i z4)-6rkalJ5*mPeSnUZ%wikytVjM?zPQ5F#J(Ty%AeE>@Su-BF;7QP1bN`0J^lXxx> zf9ui+%I(bpRoyq3R{@Jm>b#;$&k3DC&?Jf~0A810A{B8UY>#9HRGLr8ULm|R9ThG& z=%pfNGYDKmWCf701V9ZhbNT;0+b>k4tE=A0pA<=?ipJXz+l%vmA)~-CHGrpC0a#0~ z(^z|Yu%x>&lyU>67e?Hkttvz(<&aj)mAz}X-TzDPh`WOux@97AIj!?Q+qdXwf3C8M zy$sb01)g74q}2==RBSaV$vL?OU z28E%n+%wubM7pNQ$qjK$I{}r&7>*M@nwgGfMFPOSa0;MR)&QIw#bK$*0D4Uc4`?`| z^}@rRerb~D_b?zYF2=&e9C~6pS(Ion_(f&? zL95$Xj*N<<-igHge~*LHAMZ@h6sT#k zOG~fV;lGaJ^ma}JxN%3B@7HM@mS0FN5%mjIiV+t224OQhVCgc8oqPeB0m?G z`u*L7Zg0X1>fYb07Dyln-r8xQL@*vFmBFj+$+&iz#3WkO<_D@0U|E^gXUHIK5Sw8i z1mFdTuq-*!)k<}6ao>M_4uX2sJM2gTWY3y(@`{15TKGwPLOwK0O-z5-{I*erMjd)&%R9TJ?P8-8ULK!zBkH12tKMT z*kzgs;y^(xKL`kId7g0RWN)uT0h9{+-XzY_Ea^l=ezH<0VcgkDQ@p{z6_IeFe;F%? z8aTe%GJdS3RCz&cdYI$`pG*nev^5a1Tqx?6;EyYvs7DtU;mTwKULU=61q2_$=%DR> zsG~L3x^K&007$aQ`A_VEP&)-Esd(RUU_Q3s)3jW(`E>VkY#Hd`1eC$5rJxu_HCelj z0eb%1Ltn{Tf+|Wk_f(+dt_`JVk#Jbh@!y>)6^RU`@#|uTjVne#TZf0`$W3vucj?p9 z(+g?=8Yn>~byuS~vkz^`6T7`q{m1#n&%E>P_4Y(Eu)qA$WH|SUo3OCovf<0cdJ3$2 zsLk+Wi61&hw5+2B-uJw1O_r5z$d=vq^14^A+M6<##sJF7r|zDgg!BN9|CX)>^NDJ1 z)`;rrBd1@B*FIzC|>-(4(;*`I#b=DOG#$&ADlaPh*k8;|`gPD`_l*vFB<^*JC@ zL(j_qf-JD9ue|u!K%M=^kDFgi2Phf+iskgb@ew)K-{PjvydmXSYfbh!3ZE$sS6MJp66L9)#tY3`610ocFS#>sQ zGAvem1rSTU0LIIFGn8^ft2c&19@bo|5OGF?c;>X|BIbI!l}?5Dy#v*Z1c2{$8w`l^Bw>WC|EfD}yI>F6dpc={cWVYM2nHqED%BjK} zE2!2*ihzREBj>W(l;-&`*q}M(I?88a<5Zj(k@mfV!_;01=VL z?Zie{i@AT?`2NKur{dfZB8ojebrtb^I^7DNx1D)qwJ+yA%bAHm+(l%CRmVG2G%p;N zRgJ*5yhwa|=|usA);I4fZ%E;!dgtVRw^~V$w?>6bO~VF<8g%kheb{Ln_B};H38L5W zOzCLWKx4afBTx?P>%&!0%;Hw}J(wsSpv5M)t+tKqRCTfYgLK{X(cY5w)AO!Dkcu3^ z)fOxdu3GNSlQiSJamzez9D({3wALFfHku>DsGdtY#SD5hAbqhK z%efprA)K#&2^!&NWrow@s1wpwl-c)F8s(<DAsN;C~i+n;S&H{*GAA3pwSOtWT zaVuH^saczlMW?9RdZ9kMq|M9+kR)CLwGzlCrC~pZOSBsRgv|i9r@joJ5`=V_(YcLe zh-{qbErvGA0$@&|oAsx8zD<$(H^IcXOl>%n5e+3l?efF##w&9Hs*yzPi{(ck8W8(M zAiiIGaiUju2yv#-EptCPheZNDqpCEs#M7(Cf04;Dx3@X;Abyqq3)sjx>oihxwVRJW zL0l{SV!t=m@yRQwFRRLj>qYBesgs-_Ww%jf-x&w=2s{Qrt5HQxYI&>(m2oEuMGd!nOgJPYeLAakSpST!m$z9tKoT2lDGIE#J-njM8sm0zUYlCN!qv z_6^2nN8OnUIJTAYQ%5XMBaJ zuPKkiL&@K0(qm7fo{RvE-WN6O@gPDdD>r)XGLXpW@fr!4HG3=v)sWYP5CDNC;PW+# zKa!6(zcjk-6L;Rw4LRcx2O$g1Vh|-ak;TzJXzU-!K|*S|vs7w*R@Y61T5+@IdF%rW z@zVzjDZDDFMH;!=4l*f{cK6qrwPuK1pI1@GOSI?SNhGovDo5ZB)QZt2)TQxx)PMFt zRw5U+uY5Tc$KLx@OZgrt|1iwBe+;>9MF z-&oS1R@?(Dze>ZhwJ@?k?rMS)&=kzG2GLJ`LX;Wc`1bq{31|jR*Ng~Z779p8QZR2f zUP%{`k&ZPlz|zulEuBs<=mRygyX(s*A|p z@`dhpFu?%H0ofCDi0F?RNiP13dtNN=UDW$IpiQPH>3;o}fu+K`blYHZ`FIpHn^-K! zO+O=1(IzL1+T^BIqB(olGAA`!?ZJqn6;Y}3bW1`pcwOLiNp&mw!}D`_&B-DSY5?R^ zv|ixe;EcQ|XbnGk`_kxLEuWc=cS00c@SDKaXBAvoX=BFz2QIMCd(FWg{XgB> zQPs0Sln}fqHjSHWV9Dbxvc>(?y0C*<&%^m+m%Ba9bVnwcY1FES>{o>ATZbt2!;&cF zrmU5->vJTS>*DT(8=ySw)yB3x(co@PuDku23bBZ3{65ev^2af&fpH=1#kNlxMut8aIYV@QW%I6DbvdumaZfXr_;LCIP7n%aDKEmh zLS3PxxiN^Bgp@C+$igGlK$t^4N?rs=pI1?8P5NJ}y6v_Glj%ZmdAL?7DL5wqWi3Yg5Nx*W_k<{s z5z&&;!XmiJ)ANpUG3|;rL>@bWWM01NKH(szK5#rK<18SgF1!ET|Q@V!V%lv8L}1 zds0ypHK&RzFg*mA1`4x^N8YfflPx=ed~LN4TM*mqIhw&^R7G5^y`1&Ngh(S29g9-8HvE}b`{`T$iT*G#ejg~;m$5mg%fnaLza#~R_eHAL4NC|#>AqL zBP*L6RzcW+Lisl54u9>2nxsE~U2TXo4r5v2-f53LjGsUb%K=ltVL7y8m*U?} zT_1IK#AqoUNdS$*itcJ>QX-~V*u+RoZPAQH|6;%J+Py$6mT9%?ndCx|AHd%o0lXGh z>I?Hf9eNTfFPFV)*nsZW68kfo_}62&!$;>EU!UTZU0AEsXc_7DUY+3Dehb9>U~lY< z01Sl^AW1AM>pmp?r{4owzGvVDQBSJCRz!N3m3434`kDl~OXbJJ`Lwg_dlY0jF8wKdQ&_IwjndfIY~8O`jU6oUNEGsaR)c7)zD| zz$^xw(SGkwO~0ShyOS2XWboRWjCJ~*Er#V`8kL9ZiLsca-_1>btg0o={6SCwl}g8T z6g|!f23LVDQ#_LRvIzJ{5fq_%-OthzgZboj55SnL19TEV2jRj+hw%axDIXsnouVMO zzvRRU#OAQDbbbE_Z@cl-{}KI38W&4>2kMkt3+>zPu)-+buUhuLpbySd=w#YI*JCB$K$$kccEdb z&b4>)541r*LE>|#{`>@w- z-xp*GLY>^_?|!<eGV_;^L5HY|G5$`x^Y9!c+h_BE*N7;sCaeP2}wI7?F6lsaJ-j zt3Zbhw`{~w1DsPF7T=|S)qE|EcI13~Ap911d3j^@V1B502H!su_=pRXVeq(~=ACYh z2mujCbGFUx*+qpz{EgxF{n&7M%OD1bJB#NPe(-}u~$39?YNszX)&WX&%~zoK?LM&061@~o(GUh2vO?j z3m6Y1q(8W^ZobMfs|pSdMnnF>sP_D7JzlqSYKqg;wL|kyBJx7QK0lc+cvY~z_-jl8 zgZb8I*D4=$sOP`Xe61M+)7lO1hz zvFCQbEHiN4xH$Tin3bl{WqakNQ~i?EXgxQfL&+Pzw*HJ_VTEJ5Pe=hu6Ve- z+u@xxNlYLDqO>Q4I%-ceZMgK8ml91jvhRGM=0o&qCCRtc&q9P1$=O> zn`v3p*eu9jO-3@zmczIBoi!xyuFsJ51>G)WNk|2qGgl4eMmp5EGc4Sy+08UeSZ1xP z$rX0KT^Og-jLeqQBsK1Nh;{0ayL1nAr;44yB&c6Z%@x!vr1I{nO!Dgv#Wv8&7Ob$D z8S$^bJfWy0wR9A2@EuNzFMQGv$QwM6pnnLGm}pGY_GJ8_qy6#x+PIlE88qOOBb!oS zbO>lSbQq)@_yfy3xKp-wyzdHwYierXVpHPH*3Irw*f%(HK%?92-e-n{@E+)2b5x13 z2?C=k;o5PIbRrcnK`O=Hqhv4k7ifSsCsGp%yZ4AiLH}Z)h_dn{0|V&OznPz13+=Y~ zKK%4Ks0ijfUuWN4!b_%|Fht6Ztv^q!1Mc7!L&Ch8SQpHYVGv`^Uz{Zi)F(9$!)q;He@+36>!}2p}|% zE3q(Ldc!OUCQD&A8lc;%Q9uggD7!$!qeh(#X`>^@R*020-hb$MsJCQta5FsO(X|Zz zOc-aPlqP;9I(Xik=b=1X+5kpA^$QctKM7H2T(zeqHIrna#+1r65fYpoOrQYHN9t z)R*y|(ud5EoJ}Enzh5>9nLv;n#SNUQBz`#*uTmj5Tk2c&ukHaM0+Zfg#RRfO4{LtMx|ReS zG*2cSVs{Id0yd5-h7z9f#WDgZqvJ6WEAbxRx;4Vpxq=j-Pv{EV82~v; zOW|=E31sFj^6wQPNVqm5ah|CSyR>CxS55gLrQ>U9XapUaN6brB}ArUd9Ak4Tmz``n<2)~79w`RmG=PS{b$~K8--U@M?Nl! zIbu%zMunDcmmp!2winMQ%<1H6ER+mgn-nFF4`zOqIm6-rPs2aeS}9)uP9QE}!8@5N4oSytI_^Sc zNAa__nx_PQx8pQpf%;MPrg(SMOVRWQkd14=oP}q(!Jqh^L%3$O3hd?0KPulL+VGVl z>F_&|ZqH(qTQ*#V+87`f`j{O0Tc-sqOE$Lss@!!S>k3BM9T;3bj5-G>Q{=3l*WMbD zqF%7mU`uJVlo*iT1qp4xPw%1v>6QZ6g}5(mMBOiRY}OgL{c6VdL+ASoTK-=JogxLm zK{N%*`5EmJPw9#NGKvEI+l7WyIun&e*B_yedpj+6)+pXj z7}a9+q(KQ~B;%2`<&?sjhSAT{gY&3U?S|-KgaEK>uE;-7E#dy z${?iOJ0d1^N4fxxrAYs>*->`7kUp`dvaYM7t=kbsPO9tPM9)6{) zo8M+eZLn8MiO3i&jB?mG-ICkTfrHK;aJcR@3ozob#6M`$+^Jltegd{&DZo;t{MuCb zVaD%PLBQ@Ac8SP<)>2KHxL=T9tiBLU@dk^rbzS-dPA$xfVUk(OsJCjD$oA~C(&-fw zOCqkonph+6hq~PMj*eI?6N^Ar9fQ2fX%w3M7fCW+KX?;B%zz~oz&-p|tKs}HLU7_n z#dUwn9TB%iswnjLnv6T6A0DyRXw#qzm=x0FIQ+Q8Uv*f8YTVu^qWi}DzesxSE)33q!r)Q`4&UM=FsgBGX zC@#FX`e=J`%m%IQ{se1oS%0SbI8|aYKbpocbo*MYIahb-G`y2?Cl(xG_!Q()9^SSekEBId7bA$SLh!v@&-;xzo->+ zfs69Fz@~ z4lMo{K9B1jppnc1bQOu&7pi}(NkKB8u$G<0u_*UnVJ#FGkslp7Vw<>~^DeEI^6NLg z@TIH)Q(kMvF&)$S^zi1BhuuNjDw3Kimo8@hk(q5$U)v?xY7x-$pc>5sbN>8K&to1^ z8^-T?D$d%l1cWT{>yLSgl-SeRfL#`@x%VD2na<~?d6Mi^~bpd%cf1!4AfW(2ub{GnAMCv0sQxn3*#Y<*k9XJ-0|)#sk01??BXUi5$`VDf_WRAp0~~_{15hle+)pxb^;C1OLznHATJzcF-k00 z_Lb=<F8=;c_gaS|1?WPDr)_4YfgMAdcPkIDOkTh{hDP zQw!np=AZ>xmaOIWICyj$;1p=>H-}U-v*J^pq(l+onA84V2(_q`s|NtP{30zdq-B3` zabp;D0fpC2xTNJStS(_7lRt7CgEWDkcH-00-j1DFKN@$#hbkogGrVDCQyMY^ZD`8wcn>QUXu?{X|WxJ%8p z#=%U?)a3;j7;(wSn36awqk#P#4AC?(S4Z6x{l;k?IZUD%C6v>Ow8Ju7-p=$OKqORr z%KwOvi4FiN(XhJ%Xrx6W?jqA)mb{Cq%;kqXs;dDT2K{_vzV;KSP*omd{E`JV&OHK} zq*?K;H1#JM`mdz4SiM1KQtN#`7ryk4u_w_Mi?qp+_ohTYq;bQHUNw&57)6N1+FDe$ zT9T?ol@DsKv#`PJCT(6F$ggKw()zE$z}-$@x$2|12PN9p8468k5ym#8B5I;hi3cIBQVwWMY}@!I1X zaMII6#7hF`3Q^4(E9RbVq@W$HAz-hU68qWoC(0t)0WXWl?r@4nMf}gT>ml&=^ByT~ zidvY?r&ApbA1fh<`rBh5$N^T@@X(uFlhTJ{jbsHg#7B7rWqA z4mK-uzp?|yDp-|RTz|hZwYF(smWcvA?9LtVU$(y#zRpeeJV+OWR@ACl4Aj$ zvo>-M_3N>sFd9gB3KX_jV24fqR&Zcz;LJJDyP&zV;o#6!ueu6h58|LUQzm>|s4|cL z(mNI2eHLc^<)vxHwht9}4wk1_VI*S+_B(oVi%+>TS**0@hsX|p-eQnm9V+aHJV+1v z2?^_k#JBg~0e}>t@OZ`_hk}q3x3NoK$4_3_rE1p}{|`+i3Wg3R(>R5v$;Rv%kM}d6ftFOv}G$gzwZYvI%4Qgz`nXw54ChZL}OO zx*hxqSalaura-@qq(k4}b^hg(2t~jDCh7fxtWTs<+%LSF=K+irBLl(=--E+6)9R)b zPs(Yvl>~HB_~|0Sg!B(r0JSuaEmHq63BAQ=9%vew{n5YdPZ_&6+V^zX;%q@0f*9^SW4nCQM6(xkNCdb z#CH+)O}JVyeFNB>)kVX}i*X&=k{{?##WQJSkl;+&Hf76`zjoNp%HVoGN*G-tPhLc@ z=H)Yx`bMjXsnk~vgzwev|EM3>zm)|9sU(iTmYjOj?e$LS1_46?|L&V*#|N_6EW=aE>D;H4bnTxNf78pY#HbcqoC1QWUvTvZ ztuNx%3hzE?e<9AHQL3N&3B^tKSM$p<8yCPn>mC7&qY`9wy)96dPp-t2vWK zhc5LX_Sxe%`?k`Ak9!x@p{dJ!oU=_n)zphl7SFy8ox?u6A*}=&|MWj2T{}HsRJkJm zV(FPZlDT`F(*EnhD(M{8@uv^NoOVBLgyVkvO77hQZbd7y1=?=a=lCsJYSb#LTD6t< z=HqYK2Np)wPg^4j;08c^u1$kDyJ=ifa#rj_npk&W2RK7clE|`<^%a zmJWQr_JjnRZtel*t0l+v0FxQHsBm zhTLbHRCTbW7xvJy$miBL^vb0v20hk&s^;cP7H{|!*WtSpm*wG)Qbtq}a*a4fCijR( zY@k(td8t5!F^RV~^Bl8ZMBD&D&NeIelcIiNnIL})kf-eA`Z-MK%eNLo`l-{kIzL?XE37HYg%7LX5zl5?^F8+irg*H1cH{JWRe99U%#vfzB_9 zY9e?=`MZ%_aO~Ssqo(?clNeNPQqG2QtP)$X%)GJ8;T_q9J#Eg`ja1{`Vl2^of*PC4 zcpuTKq5CCRg5S70t6075b|#gy?Otx}GHt0JstGO)L>ZGr?URdUu8A=l(JqfKf1~)w z=>71*sl2 z%2OnV{!GA40eekXC?IPDc~s8ey78?kJTz)#P2iaI!bikNAf8SVQ{F$nM`uT@eB!G> zjS+E5!Iya?+K`T`)N3n7#x-`YCv*U#YY)cTgiT<%O1HRd?ud%Ci?YT@0OwGvz-#jsIw=neR5Anta3`9+3*WeXZtm zOP&~fE}J6TSM8~VHer%=C-+@0aU{)GzHjgQ+`NDlBmwkHI$|JJ5I(hWc@fT3&)Nit zrH*?*#>LCG0A+~k{^;;S$s<7#{Rspkksd{{-yFGj*0VBT9}zD3`(q47m@G^0k=b`7 zPI~F!N~dfQb5!ViHhr0ivgRyeQi8 zH}3;&l`6Xo;6-8k(@TPv$63xitn7^%LG)VQT8ydHqlPBHsEyi8H1t zU*^eX@4uMf)6-?>MG2rvN7gm{vjwcubTDnT!Pf0k+@;*gw|@SWie#_9bjz>NAYNCi zRJ^5f+l8~$hm{&I`W2^}`nm-CUu=DaPp&1Kr_7QrY~kDQudv~w9}|C+J+n@uIsnBe zF-Ts@b6ltJIK=@BGfH`Xc4#-E?Lfi{hG%;K=a{$E!AIx|BY)!xNwfBS^Y0%Pu+ao3 zJ-~?da?7!(z1p6L#NC^HR$CnId8Gnps~9#la&!-Hb8M1OnaO0EZ=lzTLuRtgMkA^O zpC>$(X0FHPyQTd}ArZT9CSktWLrZTT2%aa%kJvBD#sa`z=pXT(Ny~p005FH84=Q^a z&^s19?!&4=qI*F@Qjss8D*amybfmPg-QnP_H~GfUJ2)prVm9J}%iyiKGGM;2#6h&m zV^^EyQsIAPFUV0sKEnuQPA;N3sW*TA%%kL1F?$aCh*=M9$CR%Xz8}><3{J@4xMCx$ z_{4pM06cu3t<%OpqLGcwWbl=jO640Msr>R+b?RE`7pw_X+8;znP|@l8)beV9(y+B4Tu~K|H_fNu_)? z2W&T~|EDOxLP+RRi)c!vE=BD%RdYL%zo-9jPLb%T57W(8q=~H~LTv;T0xFV{h9ED^h^Q5m>)tFPd5 zqA}0wpspISg=>14Mh|^q%5hFgzpx7>yQSUIlF}hc`ysf{UVs=~9yk=G8`Wy8qH{3t>g-$d^te+^lQ=~lxLb39S;H%h`cji7o0%Z(9^-8hmwW^C0(W0_A zi-6rb|1hMx;-D8W*q4{zl{1Vo!xLaunnNLeJ%b3oozNXe6+iRYmr-u(cD7Oi0Gt9=G@ii#v^R5iG34B*ZtMIPiaW1f}Nycey5RBB=I2JGGL3pTH<@^O+4vJepN|X)R6s{xi{BvPWS(qI_s#a*6-`X z0pw89AxJk!iKIw(NQZPvNq4t&H%NC#NOveDprlB*gmm|LAMU;1-+0Gx4E}NG;p}Jc zz1CcFepbppxFx}RGR5ciqRbrzS8Km2Dq>3EecfQqncqF^=u|`gc-zH&LRsOlQK)vK z72p*&VwkzL8G2Xb7#Jro)*H5z^JYNpV~RNXs$d5yVT@aPtgU`Lsgu-jN2K|PndzXO zQu*b?NRlJ%i2RT)1y2G#JdIC>?|rN6gLZ8`;{ei=qodI~ppNJ(!2y#Wj?Kgm=U*VdPbtVT*d!TaosM4nAy_HGqQ-Ihh;o$V5m)fbS*)UTB zU5uiV=H2B%B8fH*hu+Ir66P`zt|_6pJ73!w{@FH**5D_|8CBC%cS#?9gQxl~qix2b z(B>u$yZMjj*LSBsk@(V4@tJ5q7*n++>1?Ur^*QGi_;S>PP9De7n@)l6E`bFh843Bj zKP~e2F(7_~Zd7TZ2*(k~1?*fZ3fE}TFLYS%w;D4jmVcHm;PK_r zF8h>R@CAt5D53fbTFlHpj!lNZkj-o?SLB8@8HQ3U^jM8mApAl2$(;-L7%4E6tOu41 zkqJ$CEWhJE!UPVJ{?d|DC#uKL^k6b`D&M*tRa5Nre{2lleB9qhzJ4g4XIj%b(vf4T zl`a+TO_aKkRL+~Dzi>R)iygzfybl5YE_GLd(`DFX!PykWL29u6L`T3o6 zyGl{GW_tS|H$)_pUiRwqdZQlR=q#Dp+*~D`VR!61*zvkTg)&hywH*$mUU-eail+cL z)L$BQ!=2XxZ%P8aB3b$?%LLDzt?k2_eZ*P)hSSIGG$m&$jW>%WTf59t% zZ0M)!bw$30JGW`IpHpRpJ*jm#Q+t0i$VIxBVWlq@ygAF{yQ;KJ@pX{sguDn85^Gi{ zTo^??H7u5;B4(ny9F46>j9&{;QR(*H+lCJ#u;95rxgr17&(t}96^+c7$R4y@Ya;as z3&`l18IJr=r!l)KL1y9WSjdCUjhfuY8ulzx12qB4Mmy*)tghqpK~;2qHS zB%$wdV2Bjaa{wB8X`oDO0*hfPsND2Ej1rnDzv9^!rqBm^)WaBF{23Y#a!POpw++){ z;=vmjj#YT!pDMU+<3|4P+3SuCo$y}zb$;FLZKE*yRw%j=mCT^!wNJ6Y%>(k}CfPxP z^(@<38PY_JM6YO|+45St|v(pQU4IKGGN zjXuh~pE)X)xQO%pXK42 zQ9ZyE`l!fgAV}%Re*>dIJ7*xfo-3TkJb$b8rUow!ZMxdPVDvY2hGL83=9H12pC8v> z7*I4B>&VhBbf-pwD4^u1ocX|3lbRRJgb{@&sPlR9m7eiqZtU?7mJkHs9;zE7oyJ)@eFh3|Vn1Y%EKK1@E7Tj5bbKkua=N&GVSyFv;% zVyAhDDr(ZO1-_@*cBV?ts?h+pP$+||j+uo1`#}d)1JiKEju8Z9G&cd2NxEENnioOy;c}o-TKP!u4F0`P+DCwR=IfYg7+tl5;-rZj=c`1@gR{Eje8C77K>$N+Ju=;0kS<4=; z^_}bO>RrM{K%{JAI8||9ESREFl0p`%k;iW7g-L%@yr#Z~gaU;WlRr65mj5WKFt`<` z6Q1HoY}dg&E9|UbXZAP&Qva9Xj1pa`4o`eE?UlDYLnTZA@Qi_q9qhRmd(%&|Xs@O$ zT{@()7rq}XOJq{iNs2=5&9?529k08ep6;ZhB`0L~41_7RM*Z6?J+VL8Ux;+NEM{pv95#)P$$qUL7z;wZc6_8(6X4G&{$_JCHP@uPw0?CUsz z#kDtaJWBXe34;$ zm?4Lf#tIGkB;^&H(x)Kd_G|`eVtCG6MZxQDcOlf{2z%eK*0}(V>Bso%#=!q7TVqOp z5j1(fM<~cUGtY}Sf$EvC&4@J2H@bI03z;t$odb)dH|~7%+*2EWm~PFDcbVOOQl$M2 z3vXQ-{i^!CexDkt{8#MEiV&56eENCt^LI>nB|MTslZWZD3uqmA9afQEJ=$v-FISn- zIVM-Jy>P>bhv==Y9Y{lQ?WZ19qrKzV$;rfV%6<5WA79%o zw;)rjFri4cjQ@>nBq`uZmMhvWON)!Yr>I*g`ZnSinylYn-m}cCiXF*&3tekX)O)MK zi2mL_3g`g_e=KpW=t7Nh(9<`gG5f=%K33N16E61Hk2>&Nt;FIuUd>pDskatNvmHo~ zgh?IKt))~)G&djfyylNtkbZb1D<3FQO+mqK5u1Y{sOJRPn? zQoUWj;}-V;6QtK0b6m+gJX8!?_{qgCa>t}xDs0^SZr?6p)~k5ec9))~eD)`wr`;;u zOJwN7zS>5Y&OUHv(p=?8;9M+U5V9RhQ2@jTlHpe-`*r@bWc?q&lKMY!}2n;R4aKI` zJ$@#m9;uBp@5hGR5sdic`Z!zFU-n<9xGuJqMEyEX%t*w0B{gy&wsKVKyrTkg@4*o} zsBCuMvO*NzRUw)xsD+$j{_j)Ab%z`HCfE5@X>hrNBy`$+lLVwb00v9}o4Ye;->R+5 zP;AE4L^;h|UFQ zSJ$r*Z(nQ_lrXi_NW~{MW6N&^hF!))oilnlB>4l#AA3IG#~}J+!kX(emIjag@OmW( zQ!Pr1_nxYOv)7O626IB0`ss`RSusAv3;U;>JKrOdXJl%4eHWMc+{Y8EASx!NUg`Wz zt_$P`3?+1Bb~LFBibwAT=(L`7Ggeu8_o;gAI0EdVfnJyX^*E4wG>&?MwLwQU60aH; z&tow|h~o)~klP4jw&nXCzGc1G*U`*9uh`~G%KRwK7*l|c2xo3UWA9IuD+f~Fg`nwi z1Tiywo<<;`|BG69bTRm=;{4_{thrk5ESO8VJ>Lx+6yiV;)iY>ymI>+TUwOkSLV-@x z=R7mwZg_!^NS+4Y^L%?Js14^i9u2ELquNB}_IFW7ItC*5owfZ{x}7^V>5R-an)^h~ zwF1lsE*)RD(e#-LM>ZCXE|>9-XD_e8h&ZnOah_z@k1(7i{rWR_z_3$8`rs_H(Fj>9 z$z>5XH+m-Zu7BIl zy*BUq6&#(%ji*U@Y(A25z|Kvq#BPJEu}tm(EG-d?xBj+KcX;VFGw7H6gf4Mvl)Y$q zf}zBDL3AdS90%t&%V8I_z-Ju69lYMc^ebV)u7kcO&xom#Thtn$Sjuj4Uuel1N~bNH z4YKe!xhS>05hIh4JY*kIDQ!0$DlprezA>6r9Ja>wn*TzYG9SM9^cUb0V!Z~V7Jmwx zeh6gCK<`Wo2B01;89hSp&DnJj=n1_d9_CvwNpF+cyg#Iu}C;0OwFyM{v<-}K!CMO<;>%(^bHRZOH!h_j2 z{pQKA`+3KsHkp%+U-+EK5EpH^9)25|R|2A+aSF%h{0R|lpzR*A5XBF1B>^HQgY;bj zGI_LCjviH7h+c)sBy+~!qE!jS>ovb8F)G+Z(lI8rOo-2od8{v5T)I(F8gtRC;5|-e zjIAxF*mG;M{%)JtSDj41WCfP1=_+-jl~}WPSAvF?GYo;`0_geC1PS;y7nhfsLKcI3 zMtI|XM|(anzE;@H2sNI$9PRS>`CNfFI{T3j`C?_oo9}Mo&xo|B(jW&mK2&|!=e+sy z0%#7gM2yR58s$cBgUl~$%O*F6@F#@@A_+d19BeC&ylt}I2K|uVXaK9T?^V}ej!@Tk zMNiQZi-_)X?PZui9Z_rKB$zS#YMqMUGsoI0PLCG6T^rnXoA znotZsBYRl+f@qudf-&}i%+gQm`-hYIUau2g6OWeLm2oti;I)_;E|%wIO({LRz~j0b zV(w?mZ2d1BbK1Jz+=yoTE*H9TTR-r|wA%5i;diWXAU+)R{m~pIi5!3`fewrvCtyDG z>F~o(^9s7Q<`&s?B!jXZA&nU9&8L3^Qv#n|0}ewu3at;Ea|7_4{%$hJ`jf&~u<>wu zNrVi8?JrTZb!E}4<`tuV72M$*BDtrNF7}pqQK38!-FjBP?u7E-l_%Z2E$5*vdRQ84wN*GcdX`hSzI$Or>Yw|)*2 z8hGdBzghFG0PdX|^ety-r4Rudclv|`3tG6*GFeP`z zhw_-buD6-j7-h}PYUNY|c83wURL|RiXcGDcD+!^mc{yJJUE`lGK^h&hI>*#~0iZHe zBc?jA-Ozg)ooh#?D0=VTFBF#9gPxM}pl1MI?wtEp-SNrxK#syF2;S%p`{66vm{F7) z7M~v3Eq{*}B1p$H8_^NB^5Yfpo{)GO_Y9kisppnZW~_0CDgbkID>9{dxX#fIg()xO zL_YZ;gy=zQbg)W-r56pNDl^i$UKG09?^~%0+}rOua?uUVD9nF+dXz?_eZIEf;nV<^ z75_+A*Tapq*;5|xEBi8#&flAFQq0l-7OPWUrLLzH7CC+oEovatyx%Je5>1~U12$h2 zBJn&n2C+SWvvrDt37)^9(5nfM@esx*;I-^*op0@<>d^WBA>)9@H4Rrlfm3I~EikFS zFC5NM)F(eJ%QzMS#uA#XSqKtp={rC1qgo(frVQ^Z=VZ8<<=lo8RvHm+q`6s<}ChCjul+w5JZR>+<4ScuN^<$^Nqc6b7z+M{{0OFLy#S0+W8BpMMwxf z(4+v@{kb6!`z=z+6P>6r5jfn0MkeD^2QZt9ID5UYx!eq?aT zf%jYfr#~#!Tdr$3VRuQ9c?Tofj)?&6G+>|f2s(Rw`*!o8EC_kMefrZqd!wS0+g&Ih z1FL6mTBY77K<>%D|4IM4`J=_; zk4aI6bAsdN<(OCV;6_g1K7y%<{`|n66N+&u+R#_A#NaTj#y;VzBV=fgxh&Amz8(6~ zzf=1oen{b`zmXyAhnNLr8RMeScbT1)Mg zJDJRJPx!I`ni7P!sR4>18DKceUm6}D;8MtuD^{|K7@q_!08|Sqa%{RUIl1aDO9mV; zuwd7TS-Am^$)Qr^0>tkUV4hMFnIQi7Fij}6sZ9B=M8^>~yl?ZN-FrvVhLgqrU0o5= z!<6IdXV&YcBqnv)zWrEjmXD*NaHZQfcueN&-3_jY4rC2)?Q z$XADAkkD0kTx5J$?$o=vJgnAmd}L`kkV?g2Y}i^YG}jWJnP~X_e^Hh%lc77Af>)s) zvz{y>i-6EtD3eH+02))m3{e%d7a zPk*IdTQzfuxmiz^Y`4O|`5B>C|B@7@8dceQ^BwL(9fy~GLm{0($Z%_p?R z-Tv^q)K7Mc=d_ zdJezLywIS0si?&WOyO^^?CROe-PWvx~!`V`Lga>YNv)7^7T1A3?{RF<( z<5u&3?c*|&^wG@k{3?#{0CDfJB#MB;gNSBv8Zve{LWCkCmYY2RAa)-&XOqM8N#Zo6$GhrkV{L}v@-sQ*`3ObZzsJV;7W)8{KWMy6^Iv*0x<=Xx^FSnSjvWzr$sF@ zbXY$aHS7_$?5PbFdN{rO?_#$C}IVG~>J+TR+lqQS;**D94DH%+XSSEgo z%f!=lIJRT()o4L@EbHn;ij(0#ZD$k(yYrB-9kj^DB>#y_v`6s8bWeWqP9@!p!Ya-T z!j=|AAq06|Y}5&FMnTudh|(5cs$TMUNnMo7sXSOV0yQ3a;#X3rMHKKVs&Sq}8^T3g z0lPOzK`39K4@fu&s98miOz3pncL!g&pWM5vf7Vt{|9F4&6*C?Wjv?+blMqG3L6uOA zRu#gSZl=bxCss(oTl_ctuqHO))uZjK-fmgnBWS&NJxX}O+u^^62!-ckKsY6_>xDY9 z9@k{*hYqiD8{YQb1$fSEjv<8658F!{mVdtT5wK$A$X2Hdg@%R>b@p*?o~0!J@o5>> zOeH>yQi=ENJ;+Y%7DAT$z1cN7(Dzk}^YYn??CSHe5uZ*?vtjHAta)X@`3YEvgz-x{P68g6uOvX{_-XSk>MJ z9wLwQ0zwfq>2=t{8y)@~V_N8{z^eEM#>Y%p#%L(DqgsHdahd#D^TWGKjxbLQA$qca z!TuX`{%vs!SA;MT)1Ru9fky)jXgr2 zURQUWdfRW^?hGU}>ghPW%M%Is;__=yItXTbt}XYu)1*%!Gdv{maMI(P3gclBZ6;C8 zr!ib}Q~&MjD;~F_*(dq%Jp#~iqg-~QKg%!Mc|Nq&GQn*reRCaPf~R$Wf{C-(qf$p_ zfE_{VaEU#Q9MiFyv7@a!aIW!~-(ObaJOOBL99~2*0#;oHFsoL3Z+Ihf6%&sD_pX)& z>AjdI9^>%+Z%ILh1!{{_)+54hCS5g7s(apL+*sjuJpV`h<}8}%u=qB8a7=G>5J_wd zY53)4E!{u0RKOr&BRVOauVz(w79%QtQk3=v8Tw89!mwjM=ZShoQ4GP;W5;3c40!KS z{-j2+*F$$h1Aiv5{x~gu^L5;evrNH9W@S8-(`~4yxJZ z_GwJv*!u?#hQM_@c^uEfY+c(0-ClA5wUWY86`3zO;^rRSe={n{ZnGEO*vZoU**&1- zx{XR%JJ*pLJK#GxtGL1Oz4uTfik+~Ts_X*#5r?{E^ZklujiJbCIowD{F2yzk7iUD& zZlPJ35i%$uTq^{hrVku)zXihj<@vpNL3}fx2~2{PsJ^WqPL(B6FkybUN@dY0eo(Z3 zuKdv#cJIaqWF2i3(fSvaMBNJ|#(fc<*t?{t!a&re-k%gQ6oAqqFr%`Q1g*c!Cy-F5 zIy=42Zbn?&L@=Bc{~d3H0_8#7L*mf!e+KuvU{$h50lZ-8hUJ<8EBBvahbt{sIc<~n zVnI`j{IUiMO+zE(=+TeA*Zc=TnY*_kU>6^RYpy;qQddiLwtwFt>_mubBYbPg@eo=n zZg3UFp3`uvh_yVj)8}U-htl?iybbnq$pubQ+4~-OCpr!h$I~M3aw+m`De@JgY~3+q z?uvo3a73;#f?AC~9wUT6Sh-O4up?ENLYQIa&GH@Htcg>L+ySfGkd3=ys!d%TXV?5* z=L5~=DU+=2;#R!L^{?-F0xkV+WN4@nt8JLfaE&6e9q><2>K7HO3;zc9A&qp+F$J}~ zcOkP;;{5>TM=D~Nl6pCh(FQ$#EKImz$bvU|d~0rY-YM2ArQoMOp*LgIUu?FC1yeYo zTBh*+&!EJ7>=Jn1_@oE3dS=oeHrGgE*}L)-3?`jMESE||K{6A@MKuP z7iOjUD%0)n`*NBd;>TGKM>wE#R`uqKjPD^tZTB*IgCR|4>~kWr#oxUdmyxc}Y9SY2 zH_B`AoPxnVn@JAXTdv#0t*d>H?CaF-!h3`nCN0hG7ouZEOy}!&K zrkjR$1;A(LfXR~A98s|UCOWOIAGV72#X1`v@zapr-1CTh0e3}7dERqKjze=YbaP$w z(>!Y1?vo}}r<1F;H;<%0wjm->fH_h?+*hgX?mI=*66Un%wd$Ax{$5D*Ax?P`ae>#Q z`9$JoSrzWawlM9?y?mAB0eb4?hoTO-EtQDD8-x?acQJFFH)NahT$#q|9`~mncvLgC z8@OOh4|)FlA1LfUu%N(0OnxJ0hZqD zB}2q|fB%xw_j&-C&VrsiU}1tNj{}Fbjb_tX7R=ymie_kxNCf}SmG}_SomDVk0b3|> zc)v;CC37yZMsL`-p zppc;{(7L?C%1ct5{!WNfCA1$*AXK+|7LG|aHu__ck#R4tJxeFZQ_c7qGl?h`VsNN3 zxKlFiW!>;sep6Ia)zP%n*5d8CKE=#as;9zQaB2Vi%VRGgbSQ9Pa8g!ms?ZRTq>0ty z`YrDNyIWeIcjp(c@q9VB{>CPD3Ax-k3-)yM5EX^JO6EaE_slq`bEZhXv7x4=p&e_C zrQ315+Zs;IA+9f#>Z@kOV->Mv6=Kz=>e%xBF6VI{z4~q78NB4mtZ+9Cyf8WHMs^m5 z$rC-iWL*3nxE7f56#}0LA&a&+2p3FtMO}iw`tQO2Ky7cb?91P};-n{p@n@v4TYllK za?;7raJBM+-|Z&XHJYdD${z#9x8CrunH5rzqD2@kKcTJHYqe0EV7vnS~ns$6rW-sla$p3=9aseLh+rgBF|fc;~@O zc0a|I#>%+E2r{5hc2J3_vZyATS_?5sk9qa}CKT`T@8K5(L34Vh1|!nG8~2I#%8Qyt z4_Ug4r5>|dCw}!(GpUKtlXIeTW=Y2p5Rl^#EY4eDmlf?~m~!LO7JhFXF-6B~ktD13-&+;i0$ zZ)s@gCJP3&JaI@B74z2^n79L|OmG?bSp(K^5*P$2({U0L7#aCOjtm5kbArSK9bHEh ztKPzY>lffVZIJ?dARGZdtr%0LRPW^72cT#u#vf16XXVqaZ$Hb z0CnF1?32O<>lu|m=y$uBs;mVb4Wyff8T|qv4W*-f+r=7L#Hd4pQ@e3;buG3QML?h< zy-oq)#@tLKN7#BEd#d2WmB|fgdn|x`M-4a%8PzDx8D2bP>o|L5P}g$r==JYuhTo^g zy5B&tA*1a?vYhMxBG%WGuc=so8~I>Iz6db-GfOtxh7tPkPwx^lg%kfJU!=jjuW>U@O+J{|~06psyzpRLHCFA@%S)N?hGhn*U& z)&h!O&5ivjT2F-%ye8_ALrH!oncBk7AX$Qwd|^K}CR`dfX{L;ZmO3_U7YU% zxctTd2l6|_|DA=S$v37MfHnp`#Ks=aUFXm*6 z=e4fL&c~Vmp1Z`7D@Wg_ot^e?rk+7*)tbzf<16Bma!}k3C!*gjI=4LN{hcnlWxFUo zt;j(|e;{+1Dm}7s9AP7}G{b!bk(#o3V*Sx7E zD=Ndb=;Ga>8$2AtjDKcNX@_Zk;_O7FxH_>;r|Y${w!XI^{I+*n>qt{2e6{}BGOT4n4H&`&d`1e&kUYTdmw`Q?mmtLub*dcBxSlbdOc)Q0bDk;f6 zH}(x^u76uI&0Kw95BF!^nXAVDdjf@(Z}1N^^`L3VGqSs@%o1#F0oX;Rowb5yALgyq z0|dKR5O8J)-LATtN7$KJ^T&e!(*^LVXVaPnJ9DR!2UMUYP?EJUFkDkUx&Nw;oF0G$ z1$0c3I&(D;O&7|5kSkD*q9n8f{vmXjI5b5AJFFwnv@pzHm{ zHUb~s_}#@&q~YQ3+4v}0pT%yGNB%ncm`Mt?{v1u}{#D!tQ(H$hCr@12fgpHG*V=vk z5c^pBJEzI+r@Bd2U3Pp}`$cOwHUEegVe>Uxwe0M(G-t@vVLb@!c>WwZpcF2)GEb==6bC!}+1K)H2 zHW39_lq8vVbnNckLtTb*~P(jT-Nti}PVR{V|US?>BkY$70i zwwnR$@6_DSOrBV<=crf+{rZP`8tmxE{^saM`kzt@>2QtOGb2Ety0+;Kk~V`%u#zq*(QO&5qyrjO7dqtR_54B?)vAk#PFA#e&+Mg3h;_zRpeqPZkR_=9oWO5w^0W5#|>fpEEY@ zKR1Fuq8rbvif-Vh_@v)PC^Jy~CLfoFYYq+|e`F+r6D+~=-U0oq zS;9-F#1Wxj{}DzJZ{*L1x7B`%QTY}DcgDy0twoc!_g^F6Ce}__khA-QR8#j0>O}pA zSTiz{FP!MXvK;CueZ+A1bazy^8k==GjPAGOnn=9C2#lX%U2N3Gag!KWM zk{-%BCcsw#_?vv2FN3TZd_gGqMF_Hc7(?(f`Aoqw4y`OhsD}+bROX`oI7m+^UA;8|t!ukY9O!pvz(mV&&xt|EF}_v6KWupP{A-qNa#4b=3X^=vyN zH9NWOo-HhH9|Q;Bc{81BKSw|D6`#JJNNddanYiu*qprRq@Pnux7F-@Sj6vu8D3dM> z+NDBS7fS||%YJcnyl>Sroi9^sr%-2QOu-QpCD~KhjPyXm@Z4=aVnzIQ!;0VfbkAPQ z1Jl?-+lz3H^3$`Ah<>+^lmHCDFEp;ZQ!-XZ$sgZkY0j-`ha7$02M%*mpjwd4HQHoO z&E*KB)5$YqL|5k3boktIQ>M!SRF0WE{Z%)BE2_3a0ortKFxi`z2fJo_KXk=F+>2VT zgi^cieRN;4Kej^WhuN2wo`E{?9J{leCKt$s&pqP}ocWIl#y(O&R1cPMzVwxmRbt$h@U1(TrsJHt5Ao#)g_Jl((vAokr zt#-vE-@A;Pxn334 z^b6HACg`8#@g8??lBE#$w|LfVwVnGr%~Wv03^Q}5HcT;ZNRo>_;v>l0*w&L~BRN?D z<=grG#-&Uh`6q$|M5X3y33*8%C|^gcIU}r)Mq;t+ zz;%=JINm>IYo5ei5smdT2%7$hkv4k3PMZKGD<`(Ci~ruD9&iCPt7HwfD~7T$GA{5h z7TOsbvQAQZ(nG>U4KV+VByJ^02wNXDlQ|#$3aM{oRI3<6VA7%Zr1EKDsAJjK5!K41 zhFld>5YP9{kE}yPCw8H*!(5C(A8~}3ioWOf46CNr%Loi6#4sDbbj4*Q9l>H2#D3SR z>sT#0ig=CXTYj;!o{}nGyo_XbFOm%(GQT_hc?g1N>%=VNM#J>- zI<^joM#oeO>dCpPwLfKOFv^|>-VII9?kQHbh^=4zB!cAl-jt0zAh?ly(>D%Ke|nx4 zCm=EO)O4>MLiqB9E7*8&5U9LdXO8KzJtlIqn*&B~mUJ_N7eS3cAEoB9%O=P3by)0NKR#+>&#-!&1GB5JN7(-wok zHlJ$WH2)~Q>t)TekO34Pa*^v(T2-<5osv?iSDUBRG=V-sbW-G}fmq7+*0Ezp)$?BO zEN<_IdUdVhdC!{xE(_Cx@T~!Sd1i&bgI?ttNXFI&y(j%O}I2rLX zzb<7Re{~c}dJq;Y#5o(v@bM3zrX5dgqBFjpF8rqTxB+uuufT%&3vv%m0GBQa7-2K& zd5<;aGJJJg7lOY>}GcQQBn2chmR$d-GLy) z1pNhoPQ3zbI`mWu(JUl?^qQl88|f1~2zI{oSPC{9@E1>iF5Zh1Q9Ha2c=m_L3)J^{ z)x7a7vh1s=o;bj~s#SD!N?=NY!+qPaV zVl%Ry4-HI)@c+Un{BbU)nkXwVB!hVJk^{mQ0(W=6%75B}OfG#%%zh8dy%eUBh-NPO z&*+1ZX;ejY6@djFEPHWLaaa0b-tZ9_I#EmA-6R@<`o#n6+1mm2NxZG-JH4MH%O}Umvoa{b zJGa%1(MDD9dmqr&%>9z-Bbk)As(WFn94@|`GKT{mF7p?gK57oJ{L+m~8q-fxt(0R` z{9*3!XFGd);#I6aHV&8~Ts@%HqK;JfC&$J4@FBrCBO-oUt}U^b6sV#``i<^{RB-2D z&44T4RdzsVXL08uNn7hzH*whHWeYOF)FQ-d*JYVTo-_vOZ5w%Z9EqpX8b@`MVn+f7odV*A2 ztwGbN?|%!r+{3e6pQ7LU9`Vx_;(t~lzK$qwz0W{W{>$`mCT$_5nEght6HNnz$WB<}F(fsm=?$#bAG=_N5Gv%F0nPAz2Gl!LBPPz+e{5RpF5 zc@%`&L?537f9MK5<^bB^-Gi|!mo*vse=@fo1r~FXEXfn}V!YM&D5w^b8FM>DckQEes=E<&JjXEId&vP0ohx3G5Be=Xkm5R$1l{&4y% zeLF%zpe)!6vWVIF)_~3&vvt?uTfT{jWe}nnev0;B+!vF^+kESE8?VBprw^Z-D zMOsgJp@6aT$v#nh1flYYJ$3}XT)q9Xh^UyN(AZu^(PrE?;Zb?WnAC+Xl$_Oz1FOo; zN8e~$N(T<8bR(HknMJ0VkyN510@UgpEz4;MmIL*7-?!hRQ~2xNGaM!BjlwwL)U$K) z!9~Op31Bsa18PnvI@x$NlK;e1n7H}>$+a)V3`fUS@ti>HYfYEN|4G$bjQ46)@L#>W zrYs^B+LNhJy(WsdlNlMKuP>b=Qb)SH;O!jA>!Y4W%xJJqes;LOm~`6-^{bu;iL+mH z_z?C5e&hLVSm4i`;yBnuj@%5p4e|lcSY;~O3t2O{CrXC3qH6gTluUmPnf06&7AtjR zNRE|uMGd7--C;(&Aq!p%da<^PD7aP9mm84b?oAT6qV8NFS17^2Prby*8(l_@2Y&V8 zh@$7ypY>yQBSS(UT%sO!B(?g}60UYUl3r)LE`APvYNx7Heyu2pO1$X3u4vlnO`dnY_|CG^>U^`N z@GK|SuV*w4IoNEyF{o7)>%;Ce16weBZ5S6KqBRwxW5p=!GTR_F<&mH%yqpK$rKw^Wj9N8%_U~rFzg}0X8X1`V z;Br=t4v#Lt6=mRc+(ji7WJyU(Ejkw_HHU+TCqR55RJqJr9bzA>1%V=*2hU?wBB}?; zT4`8J!1t0wb#n$k&VW<`I!}MR4o;W~w`Ba4OKSmga4*n%XtQ(j=)+3B|1(ufhu#Tn zAVqX~D(G*%*yR=IIxtV?(cht4yA)#B=;OfCn2W9AW{I8XJW(@hvaot?pSoamPB}wH zuG)E5mQ%ByDrzyo1Dgo(a-iW?Y0dlg+oI>npfK&WDTG;(;{3UQKYvp>^Md*J!nT}) zK=Th@N2f!Gr^_h_&33V=R;~u0K@l{!$A;vbx2Mc@=n6j#5I3O5Km*{1c`8(YcH5MH zgueX9u2dCF?qa8!TCa{R@kR_3hbMxMXRTED5=cwfvLbGOe?>o%%~;?;1f>{6eJ5*8 z2B6nw>Hm#jh~T*(^Fj|vwl+5`=W$23hg5eD9qN6j`Iiq;J@+{!w?#fXo+hhX+5KaL z4zc+kK>;02pUSODi_qorQ^wP8R=<_>?~7pgRkEaC+>C-TRJMBg8GX8ChjRO+9{XRY z3_82UMGbe4-~64|&$gW^H{xnx(CYTtj>CdJMKcnvC7}2vtRNlhNtA+E>(AduKesu*_PLG2HxMqx2W`d7oxXm1FxOP97uN7FU<{0rmoiPed z4}0YPq3m!suJrzbrhFnV%U0-Oim{avXjd3##hK@Fe~E7X(XL=*EiM^v)x3y&nl8vQ zE>8s`>{;ZWzO4$pR@|rJV~k%gCL6QIej!k4ZotI&@~m6d2f*VN-?)GEMqCk`YEdU` zVbEOeB%xdUXo$q?c~f5hghC7(I-79E+gA#i=B#~U7nX0>&rosNU$ws8_*mmo(9i@c zp?|nKT74$@SmjwQwetZ5+8iyZCkMi(I5o(rr%4S2VZDbmrAKN!loI$L);ydM`V8no(RV z5(KZ*XRz(G@Zr^{`X>r(5cIm5f|cympomSPyvkuulLTsNLC$YyF<|Nt6trWb&|^xl zBDiF8I2AYHtMmR`a(@izUjF;wLwR`Sz@&dA>}fnYAx+{M35# z(*V+oMnsKu{Efw*yH4;1j!5XfB!74SVAdp+ zyv$4zy=MCWKwf6@I|Bmu)}8!S%zwX+Yk0RERp!t8thUYHCoR0Hc0G4qH`hb=g#)og ze_#=cI3XGy-|tc{edRY%=~Ys@DkjUmu`^*!#httsRCu$qwIlB&`KsWl3&OtR4*>>A z3N{LMj}uc!Y>DF772CTn!C)!Ti^#8Q-T(`08eePRdo4kh`Pz*Yq;&nkno%ZesUGz^ zK)XBxp+el7O)X^~H78f0QKyOaOsZSVmk(#V(~+nEZ|0{}2p{5n3FJ&wqxk87@PCV6 ze{9}|;V1mp#}HC1g{>dOU!qBd?oIP3+n0)7>)FTRjlfhfe2W4xhlW{Vh;7Q`G`%yG zd%jFNk95_U_3@|h)0Q!I` z&;U&h@4&Fc?1aFD{P$9)&ueh1n8aYbP47pE%Z_|jLkMe_Q389x$#a4902-(S>1)>{ zf&7VLuge1!2HwX+QP@mNu>ULpTokUaF4$iDsNPNj0@rr!*a`(M)*WndanZ?=$y?{Y zy1`3yWvNK+Yzw-IlcDi}sm9O1>e{pXuZ5$xz~c5|>W@ENibk6`e^#q^NhICYPsg_t zeTFznpYv5L_{>c+$sA(ZO4Fiw%}&1Zcmb3EBOX>|V_o@1^L>r(P*gldOAR6WE?K!v9SSa89Dr~~jgIUJ?DG64IE&hKoTuPX2;OW{KW?Vn*y&ylprWEG{QDE^ zx-11T!DszY1z+$)+fKM5$pPE!+6f|7fR4H|Cr0ESmsbo148wY^xT*ySujj$D{ksyS zhq2NBm4Z{yZz%f@N0_7W zmAgoc{Z2$gBw<|Ai~S&M*cpq_MdV#pYViYeMXo8u^x{{@!l%LI#dOF*R%w{->quhLepQmIuR`s=Z;3A1wsMb($>SlMHgypB6!&CpEpaLJ&LCZOmMq;*iL36Gh}{rQ4=kgw z+I`ITrtg81F>JU1!XzphTgfCT49Bll$Oa$~NZ~nLx>xS`KZ5FhxS_!Urt`(#yRz}n zznF7Dm#QVfp=e4#TipS)*e^cV$}5-}5Uw7x=d$aMNd(MOS5^wcMmV^x8j!-$$?!bN z<>2_&S;G!Ng={^{Nro0HXu6hNBn_z*cxkqPm^4h#-iTg;KtUZR;96*3JHR)_LNcje zY6})KVJv}o2|%kK9N?w7JV2kKS#Rp zH)B-V@C!=iV@EzbVxuD*E!6LL_C>jrs5r5$wF%Xb$q(wnOQo5x< zgdwF=8YyX%kVXb1zC9k#`@X+Fzw6=;uIqS)XP#$2d*An3>t1UH9`sW=z>>+>)BZaN zNgYf54kq_w-1``9B6_|O65x08`xe!h|29_9fy{$G6-{A9y&7OR8yS=Db(yHGhWB1Q zy+uH=mQ*Kyv*N1qtH#F$)t!}FF851<8lPvl3SQ!eLlJn(onM1G7$5Tr!tw2DnG3QO zoT(|%v|V6&a1i)S-1;JY_{KEO#|{LEag50r*>3b~G(AcMQ-)e~E;CC)*hXouUz#9p z4v`4OZWpf>@eYO6b|1xB3!`X$2mYq0BoeoeGbSaEItQymy?_F8A0#s0x)AYjnGSZ{ zzAsTQctZ)0O9{VEG&b3B%c$x$5X;zuqiY|4eM{fEk3Ro<*axAryVt6UR34}y)JG@Z zzvfr}Qf|cck#6_MuKIV28;wr4&L&UpQa;)uV?yGc}d!8Yjn3)uQ= z+Kbuv1rJ}qbWISXql!q8lV7Du;ey*Qn$PdVGkhsA6H!4DyF}IprBWJdTOnP-4dFJk zN&ZM~cx1&_#84e8Mik1PDKF>FheY!h%)M;fASq_lN_#nQ4U2YysiVb9#`Ktnxw#(W_UatCiAIv#~AVlBEkP%$O!`!GWt`FulIKK>d3wuhuZD5S3 z?k9rvTFpOG)RknF=ly5dq z+(W*Y_@kDDFJ_zh_fI1=Z&OkaysH+j5nNReC-$TtS?MErgk=+vtq?OWd}}DIYo$sH zgO$L?L(u+^Wk(ME`L&mZ?>f4D>X}!mO(V;~uNu})tNPJGyOE@qun6@qN139fbvf~f z90K0_2R&Y2IAGg{8DJJj7XuQ3bVa8@!J7%7ZPjaFUp6P0x9E;z7o`UTZ*Miw23H+j z?LeU+=sDp6r-{l+woU~&pX@krU>$zmW*N|rcGU~qZhiF|fuOtPp_-zCI17)C7Jo+B zuL#zn*t2i?Jm%TYKwb_)yGZ;d*$WOLlPW)xgWM<+bt?ke-uic?x$pQF$A*$>qUxS) z7`K^v-l_iOZ5Mj@Zp$y{t*C4Gz|l>BLSN$ zdv9r44ZZM>hBIw*3w4NuFA33`nP|E!vM|0>9e+#p<5>=pii7$`k)c_nY9{jGH1qDm z&l<}g-V}O6?TFiq1>`KJ`XUl3A9q=ml$P$efWY^Wi%jp&tDZ2I$(R?-Yg>hrrk#Uo z(WWcIalGv2xsUFe-lILHwf0w5;XG)4)G5Yt*;tAiaZ$*qGKbVyKW0zTowA}EmI;6O zgM!5&4fWIJ&|+w7e8gv|H6)e?PUA=wXk~tBk2-hVf_0df>*B?3sV4fmM1(LEma0wB zz5w3v;WNJQ^`)z|zxM3bkCc^^di+H|u&dF;;Mif6?Ky_YZ!f~HQ2ckTo|$Gkikd)s zi|_h7XE$NT!$^1bbZWcr`7+0J?i)E)MVYfj}G!S%RvZ_ybXXnN$hZ&n2Z)$2nL|$&I%B&}6aBTnhm=zichn!lp5&i_>u6aOa zm9EJN;mQW~-{rFcD!+CzE#BczR6a0?raPS^rZ9;zlD=0uUpkj1-jwjL_?-b;)uD5# zG*9ndR?S2HF9FY{DS{ph=A8xcD)BO*&q|H84>eW>-Cxdfy-q#X%yvzMN^?{6oXGYY z_!HC=Dh_8nK!4^i(YO;AD6yIQ?elXlyQdTsJSs(6Uv-4mo+31x*%wHK!cz!&wq&yT zDR&+<@CL&bCL&!(4mkz~nMTR}$m)}sh6E8iV_Hu8<>oUw3 zjj@EHZt6Zd0qqwE&{2ZFq(V@x0@sTXr9?+4gq8Czs(XaEh*0l}_&mcaXFo3I!f|T$ zMZDNyZQz|Q*O3@k1jwI>4)P}3*VOv(7dDWN zkk+~Bc6&+T#HvTn<7x3aXbOn%2g#?S<02nT-_B7Dam%blN>-rK*(C(=d+R;Yst{t% zle)Cqm7y!Lb}FCUtFGcYSoh887Y=EFbTmokrxjuZsTFy^vZ?C%$H^A*!;gPpp)`() zM|0-HvS}$o19J8F;wJCn^J1BA9}ov?KQcc1t$B&p=WKlxk*#4y)%h~mk@%$z0@cZ< zooTshdR%D2Q^eYC?8+?=aVK9ezeCcre2$@nuv)mz{yB&et1E$o1yX_}8uOLeE?FA)FMlra z18V#CEB7B`d9&JKoV#9xpMSp^>-Pv-qlw9SRqZ`n*AuJ1sWm$i0u0>r*mFA7S$!SKN{3cBd$3s;L=>n#B$h=yEcIB7$BI9 z8Y1rWrwe|)@PN-NiG=-cdJ7Bj3+Pu_aBpczTP2&aj2e4v&8r@^g` z_P^=Zsd`zj+1nGRV_b6D`u2g}5lX`sZ``3Yh1~w}SCD&44uzd#pS--Ox3reV(z}(e-_si5wdV~E2+h`HqTTPqA;6RJKq?F=90qy0sz6ti zz193X2DmCqKLeIYHWT*P^DBQV*Cl88rdflU9~p~WkAt{5bFp4NEZt&hC7E#}av~#p zNVj@MHRXCiuIi_hDh9jVtH-s;ET )J??Zr8l_E|f*f4~FVIcyogO{ov{?NXTJ8Vlz-meL=Gh1?)mvh))lVckRBms#So41BK)3DT;iuZ z92B zXfTz_#F50%x%O)GrC#89N~QU>V=BtuhaIf_OxCARt#3Xcy6q8&BZ#7b`yQ|IGzjSx zy(+;Qz|H(}HGTSLjb+Ngsgx)?Y|yOD=2Dto<4o6Lo&pWfx5p1E}iTr%rw zRHAp2n2R=5fJ}rHVV`v?y`!3@COwz4VEO7NVVc88_PVueO*iMAWqwZdr5oQFk=|Bs z5t&#fe&6Uet2&-tJ9ld^W||A|uq4e#$AH|;0+_g_s;cVK_si>l7uF+&fTbmC<8SZT zq$NEJV=iV?{Cf7pN%5XTzbpsza~I*@zA&xGUPfPw-kopX%<4~`2rjBzVS1hA{Neu4 z=zO%j&-bEz+cc^X#d_H2;bSg4!l%%_p4pDcL*o`-X;8xPwFfQ|XCXG29yG~mo^Wu} z+(T)Az{)#zhAja4N+~CgZZtqXEi-FLx9}a4@jn#4@0|6lNwAxP#6K$_v;P_B%9kvE zNU>SrWCi||1GF$l?Ls9{*T^v zTzLZ6wi5aPKcLvb#)aBC=m z?e)6y%1j}b=d--V(~kUb4Ha#Ml*SXs-tXNRlq{-Et%xN`I56Pj8*TxG{%-<1oItmy zLj-gN+0PEvdafb~z(}0xYoWW@nVF~=Hz48Fjw{g4C7x-?NV%u)9xo9E^yB`5d{aW` zT4v_Mzso6Z53ylXofPvWz4-k*_bKkB)z|k!W7HmILGmN^hjt$%AC0qjFh!w1vsVBM z%p${g9PCP70Nnaf^+iVAeXsvBXll(wH$YF92O?#`gnG!uVPC*TKk<2U?O`BJL>t>$ z^i~6m(sI^A`1)l(^Cz--1Q$4f!DRSKi{JU#T6?lfSZFAf0BNUE8`L^=1K>zaY(f!q z4{w5T^?aaA%v;B6!|&PyReRr}mX?-w&^<8bKB4j@U7EFcikPIbr^Si@WWO&^XA=MfyAQc&)zKEsX4i?6`LE$Yi8VOc2@yR zprlj60oL_nXYu32y*W!Op+DF#NY>T?E)4cSEwMlS_{$(JobiSJe!{3H_6OTPDouIN z+a=sI5up?LN}*IMjHLGz88&|e;u6$dTwLPc>ED;jB0HIDH4FmTHV)v%0qaSrwaFTO zoalu0Baq7O#WH^f$eG*6%1Q7>#vRT8yzvtliEJbXzoo8?m8#>;`eYv9R-*>e&~7)D zAv0V7;2=-Mh!CW|(IxJ&ZN93`#$A~NjQ%okcD?06Y_i%@cq$@$s-77ae0#|NEGcMp z?|XkTCw{!|>G2@Q%~)D3>55f0swo8a`A_oYr+{O_kKoH2CXrvCItpHFVEeP4#bPj@ zI8_|><}6;tC-K@XXjKt=CaRy_3hwmOe&2h$VCYp3V+I#R24y!ZY#=>o^^*|0hyGw? z>f|`UPIA3Il^9hmahXrEgfOKC1qDri0WwQA1SMh)zm*C}Pp=B{I1)&v{m|P}Y!WWa z;-JbAodS!I!;m=cMCrEkWHr0_!umQtgYz;ooGT|>rG-oFl_t4e1A!HVGAoGyOD6a$ z>H&)&%Nd2B28m-Ykb0+~IK>&GK3EZjimNJ#$-aRapvb%ark_rJNW!y07% z5Xd2WNstPq8{J`9MJPBOMH%ld>fp0XxjeG)(FxvV@8!N0V&rUuaD9N^d1 zt8n3xZk(|@jbBJ6B#-GkX3xST1B^b(gM!f=hq$RJ(>s3yE25{G)47RVjV6M?et(oC z258wd@2_20VH*!D3?qrP4EX^@BuBm~v{M9*!rEd*g`;MGd#30~g@-xxEk%GMs5EvN zPl7g-x}WC66=LOnXd(7~?eCA+5Qy5G`;U1o1@@BU5P@P7OtoTAoq&$jS1%$CXcVST5y+r5UYW@XmLFg&tA5(;`h z@vH}iFp$YvCL>VSx^LC$DKb((f^WcyAOOI~unoR^qep!+`U*uP-XvyJ>Gr@I&RVpe zw|$TYWFIR325L#1c|b_f2pVQsI?@{Lx+QlWdSAQn?kf-n$~VP-o)lpQ@+@Ad=fuPD zAF?<3-(qOtbABRGF&GMH*jqec);`YnW^cPc8}$lnQati&X3|?HF_+E6zbiluQNk+u zhLuhIwBozFJ7>#PJ;ly)xLVUsWviCsDAS?-ZT&_T_jZImBNSe}!p&}f;lH?=xgmCs zKjml1XRD35aZ}UuyseOPS?o2W$2Eb1^&D=R0KG~qDv4GdU_DQ6XFgOx3($KbNVkdQ zoFRY3iA~sr$f2Q3844&22oN%wkD0p)?kEsTpq6^V=!gMu1y_TJ20m`ZW#|}-{#f%J zQR*YzDaxD6-&B~9sJ)H*uKT}hel<@Bp|Ntn5>f-G0^geQ$J4MfFZq5FhP4AS$?oUF zM}B|trK{+jPK2Ftzqz6_c73`kMTA))xpX(nrGBV^WGBKG<+MR(cT4gqhF zycC?|4_f-wN%SJQ3=&NKo`dMAmoBZ>TCoe*?iZZ7U95=&2QWbyxIjH~;i7KhZ+wit+2%Amth!sCklHto>@n_kXn=zvbq7(?3*i{d9kVx9dEW6hHH_aCh79H%)e2N5nI+YoG)blHgRn@7V}iF7+y);&hsEg5 z-0o==9MS5U9}ovgrw+ug-mRz54k{wLC2NVw4A$Sv?S6x20V9>YPmA%P+#x~Ij{F?) zv6#|5*)P}GN)`3hnM9RA_QAB3rmJhkF|Gc@RTGH6CxP||K;W*BLgeVf84dbX85sX^ zyg=uKMwupI#}8T1d^<#sq|(KUsqn99@6>0H+!Z+|=lBtWqvtkAK-dze4q1{>{WMCP zn|=mBRowx$ZugSODwi=K1ZAS#78GpHL|RbTWh=Hee3*6C*Uqcj&n~&h@$T1Hv=by> zKs?O2MqUg{(#x6bClB}%tcj-%_5S(dE@FrL2UGIHk>Kl@#Fv9TSX4N<&Sp z07io1S}{pHk4`uT#_5wIGTyB}oPM>FIz(`E9v~eo+~IAQx=U20)6i=*yWj5)8ayVQs6z5MlqTq3Om98$EjT-iNAKy4gg`n zmpNrc)7cTJ#_T72tq-XDz8Z5F?$1#WJL8E8+CTS~s8$>Bb_m&i7wE^I{&U(x1D0?7 zqI-VtIGT6uD-luZ;}?)!CzI!DHxjdA)ajT$?Zg5ch!Wo95R!Q_cvke=jIlSJ<+Q8C z0P(yu!UmOe*`$f2Sg)A-w6P{IOz`7#e?yV;E^S3d$@A|Se!Xa;`b~_+Ibdr<>&FJJ z?U-a!NCkEIs#nViK)$seLu*2h#SW9G9{V3QF0qP2b=Y){yp0z8AXC|n z<=+d8vnH|$$>8%z5j`uJGz`kWDyAPfVk`0%G;2laS*Dj^ft3CM_UrxL#iV4Ify-|5 z1Pb4e=>MWFHi(jp%aNOrN8h`m$64M^l)I7gve3<~`!} zPm47#(Cyt}YIO=&jbk1>s$P{rF!}m=l50q)3=5k8JwN<|~;o)*_t{PyNz< zNg^`Jv^DH+-?l4?yI2UYqh%0nW8+J#_(44#*na^Mr+gofP@?bQqNFxT;nw5A>aisf z-K^9XE+mXD+oowyF4f=4e15Im`2%w0z- z!f)7AgM7cL(FcyDz6sx%CLA6SZvh9sKGTULJXmj@Pc!S_DZ8u1Afv;D2MANcYGqZ6A|=@4*VIvB5vQof zsq)=UFVgIEi?`$uu!-clK(jmLJ5lH#LmUp24+Izbt!&`}a$xca6VGJyI|cY6_4K;0 z`D*3B8fv1KJM`+r?U-)yR*RICl~qQhT_0buHHl&%O-p>ViYFY!Fsap%YItL-wA7D6 zL}kqP%t_i-guViX`JkbGG!sCwtsL`2&A6FvlPm$%Rs#NXZdVtRI;`}2th_nBgv^-nu{l4R zOzC_y#2W!Y)>rEzNTXP~N^(mq23F)D5KiYlJNyve?Q~t!rsMrp1>jKkC9xvdoj6HS zkSN0A_1twRr?YyK^Ul4F_sCi1(hD%1>D}tkIzeBB&6kuF25@+2U3TCCeBX#!6{wm~7VzHeN= z_iU9b&^=TJr){{~u`ky%?*l+=D=VkrJTl-79bD@&qXwMWu8c*+{-y1|LIQB)9LSw^ zn{eq(;8(%TUIY4fa?WL%3+zTqdT^i@l1;X7Uh00anrB$+WYjE(PfSl|A0>4ssW%$+yaXPCQJxWM zxm|pC5r441R~*#|89VHik!&UO{3SomGX#CX+*(hQy8rsN_}xltqWJ*pe>YcY*mjSS z+tU!tNjsU2BSYG&oxK#w3U_$=biLej{en+(3qNnCw>eKT9|P0M)#cY_+}jqe#S;yW zd2T0(1T-m`g|`A(3;k28z4>T`$@rZmen@DjnmPuc zXnf@ocHgGQjuD0zA#*B%`Xx!N(;ab)3i9$xp6e8^y{!L%svt04Cvh-mB(J3yf_xL^ zB+I;AjB0aK#vCI`eonlWkIm?Q`cYiD2|d&7b$sZk7G{nO3JgTtqeF>o7Eou6taO9E znF2vsKmvp5u>I|uBeNSoEJB#~M;tjIx&zwg_(T|XKM3;? zDQyQFl}X>)24-R2f?^ZodW1c(N=)04y$(!9gwP3>^x}Q*hAKxw;=uXyX@F0r>o5L| z57$(I_}1NDVZ(13E4$m*-*RDlS%@C$r3B&g5eC7C;#k+ypbTj-EY3SY;qzLPsIZ$L z7tKXQpL|^ZIt=<(l>O>_AB7yWCe4*xjnF9wxn%^#^LsYRjjBR%!8UgzQAq!J!dCeT zFuT`wM0d|a&KPQb+5i0@7#_A}I^#J+icG*n0BC}XtX@{#*$L>#FSmJy-Odu*LC61R z&B(yE_f)u^htcC5VVQ)?_O~A%JojgL4nvxs>@2ns zWBzu@=6B(9sgbcCs4Bt^MNECS2XSiofH!iPGdz2bSCo}EVY5Eyo94K5!4o#OS@M+aCQ0z%a7qYE#boL+~0tVtNSOCeTXqV^*X^ZTnR)`_HZ6 zDLpLI)(GLM4tO7*SX9^7v()N$ix#rYi8ED_hh$vA`pOMi{r8O%`Q84VEx;26A8*Z| z&GW)B;~rO$*!^J*q6yl0wLO4-_2A$71#mdt>(`Bp@vht`&jb!?GT{A4omNU@6WNd0 zrQMVbzYM_xR-m2}+)`4ewstt;XIl7P)3uo<&zEzQ*$Aq7*-Fw@aRmI{;P9|c9!UaW zhinD&#fukzO!9v8Oxr^k(%?a!57X~auwCZ5aQCxD*j6`gHKM-(hXhY#0bw1R{fv3c zEWpUFoO^ENuP|JOC38S*1cLp+_}g>Lt9C9zqZZ9WJmUi}0GY=57+G|hO<_$F*oB+; z+7AyNVY zx4K}rGGy)=)6VI&y?T#Ga)N47P!BCfFvV3O+x4}z^<{!+!THmxEi)4J>(RKBlaN#-hkNG z4l-+aWaDsKciQunk0EoGDCM_H8933t%5ORu~GpNwo5I znpZluazrE&pn*b>3`y>c@%+I|H2nkGNB5H^;`R(2JN-pge*CbPSG# zmdjo%A}M+Nx}8e$jtfAgcem*Bj7yhe2)_Xea%|+|Mg1iTl@)TSDUYsbN`X=!C(8w( z(OL?o`V7cpmL{)+;yNFbi%13k z*}EAS*;aXfjFq5t-|=3z<8RydkomCSU#nm~@3n0{5TN$$>umEVL}5U>n6*}T69LVa zFdT*x@CFUx2C>kY1DaRP2fA7{B&sML-ieCc!nFv$0wDDFJ75fNOFp-|&(iP3%Cz=_ zkWSG)cXFEMT%{G1Lci>BKpUh`Dfx(y@s|Rq$-7U9-Pv@*@brMSEuANM6@?{UTuQ`W zIV(mA4g{HW@-J9)_jC(D#j=~$2q5{}Mjd}nJKTuTAC?QW=Q)H9q$wsK@ZR%4l~r z+m!QuuemBue4I_j(u0s+e!zS0$@ZMw{gHwfVmg3gF}BUc#U>X5Ri`IDTRP#NChV#b z#7GLf=={T2J_+DO^1-`2hE>{_;4|Mx&p*Ti`Y94!LMJ6rP4guw?P#~AixNxrq{sJrzXP5f8O&;;DFH|=(ltxGKjx!e}Q{`)R zX;wJwue1l;#6;Zj43KE$yfaX^OyV zl`&3qs3p>{zAr#%7q8o>|XfUJz5TTZu44F)6nF|0CEv$za_PYhywIxU{DYjdG{7Y9h*s7s1~9Zx!hP|coxKl#{&|HX zJpSd&)eyRSnFR&B*FS>}ZO|WzxAClB>e}C?pefp4iPj&8r+bSM3J9XEK{UGWLm&tL$@suW zAmRHBc6Pa!Tk?jZJI&S%H4kF-MgDegpu6`Q0X`jChvZg0zElbTYpeqL^Y!fgH$b-CCWnN&+bx#*ay|Y2pX0V(zSQl6~))hm+7QT}xvh-E<*Ri`N$YlmI7s2F&x~D$@|0jfBFqZz{ zmLfw{V){dylHRXwvf%_*MQlr>O*|=J_wc>Q(4Q2n;px}m{rai6*q-~GS6!x)pd~(Ge>ZR;Z zU!_y`9ut{A0+r_x4K zQ?zse4h&LNaka1Yh97pG@E9nXNV_C)10m-4_@17gl%@eVnEH>i8`MhyvwyB zj`8W~pi20^c5xd}pcL~5!vzB%8B!k1)Z~HUqPIE&aF-nREx=|%ZD(iamxu%_-k*bI zLLTz&HZ{+kk%O!yE0wM*O?xB3IS62u5+A0iGY`zGG&5G5|6Z*!#s%kXBM!lVzInKjqf|Gt6#x{G)$US9(6`Qh+9v?{(FKM!24+- zao^b(MRIC~dFmlq#6W%Y=|E#g=FfddL7RicJd&@$mUguTMD$gDCc0D^nC`WCy+K)k zhL|l7i6yB+X@U^YoJ^EE_8)Q<_i)H5aku0#z2uVuv1q@8QG*c>mB)uP{qQUSCk?|9 z{PTDb_#+_m>rN2Qi#b2s|MZRW>jy6hx7DE#OH0co&H3JXAe}MV&}6s)_XjJh_ijB@ z&Jb7mCk%_4#JPn11GM6{_Zx};s4_o0Ia*p>y$J-SuSfkx_adUAZW-*R@0k}5&f^{y zfQ)`^&ZvWJBoDWvD_#`)TK|xUpR8+1LO5r$6bOkaC|-v#UUgE%`A0hyR}Y%GgSj%}S|`3zyn;}2d8%MsgWAHk z^%pV{fjtcahZCfq100yYwq`YAS*3eL3IRE)4TOYOn7`xw`3g(G{DG1Z$NbB)3+wU| z%C*OCNCWWcWr8;z9-a;J%lJ?E0Hx91qmX02BG0(w{>@W`ZaugbBv4+yZ_!;MEb~Xv zhx<-TVz@1UHwI6qZg+5mO>EchaD80kT;!iCisMA!LL9j-X};%gY;06}&DvBqM;<(P z*Urwa(Q_}4Sw~Y73HWR?Zi6kQ01;k!fEwMx%{@_aDlc_UW4ZraF0S8t`e$VoZUGNM zLD_T>A6&=P-X`ZdaccC(;G_Z$wLgAQ#irK_i00KBFNWD1A^^uvWbBK)bodQKBDnni z{d*WkfH~^)eEb*&0(BK)*<|y8+9kE-JP-EqUvJGvEI8W(*CV-KZ1ydgmvO~RK+|(b zk?-3dY`asOg51?JuF*StdoT4|SeTiYI%gU^cBthp)5 zY!aBT;H(eFgtR{Ir2Io4gLa!1=dT0d1*+QDqrP>{Q_A2?c|9RMVxC3Cy=#FhO#Kg# zK)Ri9wpP&=aGh{?^T!Wfx}3UT*&DnvGWYIp3Y#?x+||*XRx8wYEFDA}`mX~Kui+du zg)%d=3Jjo0jE!t8x=Je!%GIx<2mEJc#IUX!ingqVuV@y(D7t+#_;@m zN(4wn0jM4Wuon$5z!2XDvPxb+S>ezsNF0A9?<)8Amu=QBSyJvpFza6IKJeeh8kn0}UZU5mOGNNvmbcEdQPj z;O<4>93@QQnnWHzueASb>VXdMdd~QjSeuI_u)Ail02Sdpw)LVGUbc<$=SD(;NIyx5 z4w-^Hr^XSS0oBKFlgCusJ-~dDB`Z)6{DQQA%Xxy)dVpgZ;UFC(jw}8sIA8JSisBkxN1T->u>q&aV;G zO+kw1?f61rDkZo#HC)i2O&GQd$XA|4vJy{<%}AVw4KnJ5nJ=*BgRjsr*k!&8t7T4h z*scql-t#D%_c{Y-dT*wYPFHsQ(JKL?__-GHpI629%}Xu_-Y9_^qOP^GAa2Rhd+Sz* zWZV!tNY;89^7OQs-5!Ol#T|e&aDQASVy@I35obhNC_XyBQQDrnO?%Ms|*@He&AW;X%Bm^2IPWORWDz|*QxV-sNP$cDugu~Fq zW3cwQhOM4yOvh7Jkin3*3`z{;wo3^YOY8eEcL^td_Xpa3$&eLxw6ne{iml?zBk`F$v|&vo0#3$1D7;~r)fZZHsE)D$`cf;b+-;+2+CnWQ3||eHNdUzd3jbS&JItw zkK?g8e1Nkg#fJK|1U5tO~?3#k3a7#7hua`}Li?H5-KuKAm2 z0VLFhQFi~Ji|GL}I1oJ4UgFfXGta532y`TwEp3W5&DCDMz`PS6F^GTH{Oq8JbRl?! zF2x+)ykHcVNY|QUQf*vdpQyO?AHcgMo?2I zQ1(Bsj>AdNj=BWr%A;o|XyGk~0N}lS9M&#Ppl5Q0GpD>G@MMq8dt1)- z_HbHzQ(Pgj%XORp5!ZvApB=V5ao%@FonQ{&oA$sh^p%U$u`}sUE2nNryw`Lh$H2Dr z_0DDUGdT=1v{LqDJ4k9A6d!4{Q<{FB%Z9*oGWR_yy&O!9-L!z{Ly^P}a1o=8r>)u9 zNKwT6-ECk=xpXs2HJRt1k_`vIOv{7%0zJzk@jfUlj}ETWoPyD1zebs?u|I>FQz^{| zMc#sobtN}g7W|fQGwEpxY#xn@C>#d?WAx8N>IO2~w-_h%UpM7gxT zL5ZjQ9WpMIqY0M-&K-oM#o`Dn=rF~syU5Z&S&(QSX%MNHYn4{I1U97QXeBjS-?Z~V zFil*yZVeDLKAEtb%fkyG95p#0g7KDU*6W`)Xiv6L4c0v{gYA)95G8ve`y2#gh+0At|89Fv^R?2p!{6+K~V~O;Uh9to3SZE8H z$J&Ogl{u_nOLC*Kt^vS_IK8afTz?U|V*mEsSJTbGJ+-!DYV2j_tYEs=^42(*Fi}45 zHQZtSvgTRYby_k(c*b3mNQ1Yir|RgzzqlgKxh;UgId?;X=+@n|F3@k17Mu_{=W|VX zVLAelBl+Apv>ob9w-TCD2LYN3v|&1HEAh#21>2|q_3^KPEO7Mqs{x)6YKlBx&=(Kg zRZ7T}yNp$L-J+$T*{E*sOf(1!F^asclyd)0(;uDdPcH5$cEJ^~XP` zhkCLv@t<_EyJ9F<&t}YvvMv+KMUnb&8@)J9F~I1HXCS2MQM!VPR` z3ErUzzvs4C$C5#97H~MiW%d5jBR#GNz ze4-L~w7`gf8F${`*`D_jfi(b(0S393gOemCT-~{xbnPJrn6r!i=5u&=wU0#kGqbj6 z+z-c{S3crF#rq_;v*IljuM-GLn7%Ce$wQup2t*v~cjZi9;Tfx!SP(=~1PI(jhDTAt z_}BxErN~cy>P!oYTySq*0Clsndp_Qvo?Vm?;-z{TTb5Q|HyYBM!e0@n=*j|eG38Jep8|_W1QWnIjH*JWKuHU5F_~P zpV+Inwzoq1La(zSM66Jz<=i&21M2YWT?a(Bn38T<-+9N)$L3_8Q`L?N7z}Q?QQ3oL zB44QoFpY8;aVdp&PHP1RT_Zc}Wh6wK#!H@(PXId?;7VS`XZ~dVC3%;bnZ4);54}yN z4ZhVUVhuzv^42TZ?fP**&E&+^iPBEzI}GOyvA|o)BU0nDek)+0E7+OVX=xz96-6>L zTGt<-BPpbSPeBU%1DWa`^J$eB=9{32xx-U2U)NqmG8Fp&L--VUf$m-_F-kO$*Spz-`AhijO(XjSStPnhyU6|U zjJ@=?-Jc(t+xoR6E?5??i`bgSJcu9L?s|{h5>%K)o6IN7_+EfQkeIHG$wcS`p_G6CMUJ>SBbwpZ9R7G_7I4Hb{U13D8Q)sek?ojqiFmnzI z20pw=DI`Gf-V~&!MwlxO*L*x!Xiw#!pQ34n=L4XOyS+xv~2RLoNP07wA64b1r*JG5}|58OF4)5$!bZPc~d#Tz}|q2();zi zpN!Wprx^!q_1kPd@1fD)z&it)oCxtNrNpMPPQDLGtZsf9@sRyPcvXue4L56A$3R|*{bD{>nf&Z(5`J9%eGGOS#ZUhgp?uT941 zSk;}*I$KVYSDS!?lWc+bQEKp2h2vn(w=pRgEXSuDc6nkof%}&LiB0+>II=&ET_LBK zoPOy1cgy+9=9O>@|ML%$KWGaz!AalF{KfAc7~ci8%n4#TSolKx)Y<%t7~@chdpOPj zf`{Mb0NOIRBc!_B{Ili<>a(T+nxLI4ujfeEn3)BBhiF{4Bu1J>p+={_t9p-%wI!LWs93ant0)>WYn4aA53MjodN&LGxz-20Cl`um_3H{;`1N z47#3?pV8?zKi9p7tfc1C zgCsRznpwFZS=66XkoO+SITR-Dg9fkO(6mNI;V8l%5&q8)I*~}s5DFpSue2lx=CrY% zUcyGP+!*WJZ2nn074dL@2}aIyYas-effGXB9#HcR!Qv_}CpsDXn9qVw>q)pkP!EO0 zmFKOo*NLj=co9jnnW&0*hwBX2ALRf|M5`N?Xm!;foJ;TuJP9H$mZIJ8#pvr}UdTO7 zNRuYF{MBR4tbPRL_^kVxkEF8_8m{>74 zr_jh7^I1B3@syY7iFKg3FLs!}@%U6su9qDDe&;Q>m;)oN4k!eFJ@DH=t*7@zKpIP|4 z1~LkG#U7iiZf0rLvglHJ_I+O=8Wb^Nh z|MLX8sc_GK)am!()&p^IaidG`dpZ95nc${2a3ZZ?ic4(zX9<6o^1o|3Clx(Fif)W7 z$fx_y!~FBhHarL7K{ffm-yGN0zVN>v5WtL#gT)1@ z=~MpiLRG|2$HC@q+W+r;MZvD4Q1~dqgc6Sb^T+&oQ23aa1N#5-=+N*-Xp}vBh1f%S zz+3>LGA}#!1mHpw057%@G^woy9jUI^)_q%jUYaPo?+wEjL8Vv&z^o}yu)6?)M}7G1 z6i^pC03nXr=#_THm5NZ$6l;*AGzJu8x<0}V>eze0`g=s$E5&1Px`{{1%Q*OH++lF- zjc8}~FKJ!i_qUAB?9K^_aK*V6lhqDay?$*($o7NKLuVl2y9&fs zfCRtWxmT_qsw?Q{aQrT`6OTg5yuA6i=-}tnIEaok7L?NKdHwnBq=yD5F}XF;L6F^H zEokiYLey(SLnK&o#Q0XgoN4pzH`fLMTxkm|YYS;pik4E9zXAPPH5do_3K}#gpcE@0 zetq9&%`bg;-*>OyF5NS)#${$Q=K4xf#M=NRfOlR>FBhc=XrfddcQup$sass-F=6w1 z0*K)pan;k^Ffpjn+a~W>pTq0lzzP}zY5?V57J?C}gl091U+`m*(d5u9LPHI3|_{W9%|~Y{3PmxovzMOf+$6h%c5b5es0r zV*3?j`&s~4t`2ak`jTGD2J&;Dg|DTMJsC(e^YvT@+NvHs7X1_qfXOW2;f7!^)fwPe zWw?(e4KTRp*{MWh;&*z`XFAG9i45pMqlwX01^CiI|MMdHBLp~LqEKVS%!AAS=Vch6 zM=WUrZ4?u&Y>_&`N}p)w@*Dma1}S^03oLDg+73hgm-i<;>n!fQ4z4I!|AswMTc-YU zCHgn}PbEZgj@Zrd;cgd;^akxF1sRk6PlACAZ4pYy^_R#*c89g!q<5*Iln9=-RC)mT7EqT(O#5>nc*mj!y5SwvrCwqE0<30txt$A@tWdv}*nGKKLTn5-+9hkxT zDp7`5ak*ycV?YL{Yn)YBhWA1 z7u3WBs!Rar`+zmm^$fVWPiWeSp2hY z*0*!CSiw+cr-Qoe6Q4fj+?`x{Uvuhyap|;{s!Z;o-Op>m?oItAGmpRk^84fG859-_ z@!aY8zrfCPNdaYW1@c6Bm3e-^s2NE{sG|l^^5AlA4W?eFNxX%MU&u z7qd6r#KC)w7jzc4Bz}KN{ldhLD5O45HkujoW%k)0`~@D&#KVs)-#KB+z(6@&wsK1a zz*K2#zgkANfq7Ju;^{%Eoe=wp-GzpR$KPoIF=-P1%73B3{%IyMTAMCgYiyxmD9v86 z$aSqSe@yZNz{lvaa|a#M`i|wOw(RwDDu8h&-P{J(+B9A5^rJp`8fbS~5N+u3QslVX5F(^>xH9kQxf4dRrFlj$uTO@tjr)Do>^WbxL5AqvL!V_RT z@nlP^rFhDz>b4+^T|i-Iq*s>^>;xA$<5C&?!Tj{#(^(&P33=l?ewGuvKUA7L_&7FB z;0qSAL^Ic>*feK-^pYS|s@S%_L5kX(Gl}M(nU*7hnGVwVWbsp6- z5Uj;#2a6`1?D*YBR>wjATi|8@w&qH|iij79MUQmT;=za6tfM#sef;#BO?)@<9Fj0f z#0dwh!yiwDx{`w{NN-Yq(fJK-L;`0V6B}eh4zUJcVlIayWI@GNUjS=)M%JG?%j+Kgw(jn1L?Z0a|_X6bHWXlhR!LcIZVJoSuzrQ|MvhFdJ1fwzH-J& zEqGg8oq3{Tx4H0h?B!a=xK3gnFe$oxwy-;9@?Zi?0$2^I7ToyMNZcn-6{Q1PShXkF zzATU1Zm!R=ZJaL3n;e7r&&kI1vh~RkWMVZ#bvu?>=?b)LFJgk4}BuH=Nj;Wd`4u(@e9_vs_pBH8&T)kTQo$zp{f%X+kR(K{75xeQtX+i zSG&2%x(p|R^IxA$q>H-;)rV5M!+Nz!;s`0K0;?Km3U@iH4^a_Ij7U47J7x0hc7Xpq zpAcAl)-&b@T>Q%Rk-hmZU!{3crHq-7*##Ry@Y1*6>#I$1V?PJ3Un$LXPVr}^(st+y z<%M&@@g?FFtJ#-PhnDsEiuPo;%Z=>-$%?LWcn*7;_Klo!gifbRKHmUh6Ml6{ma`b7 z8*Fd4OG=DO7JX|)XbyTX}nsbcmpUKFHveh zb~KK9`SpWFhV$$^o}Gorixz`4K^PUJIQ*pA+4XPe{KE2uT%RgCP>ll>xq>ZOXwc@w z9no~wIhb&rd{lzt4aMs7ejFUDS?-zoit8(fqo!QswKaTM!MIb~)K{1|Y5b3>pNi*0 zcXc@D^QHP=Iu*?D$I+e}M>H7L{G;9IzjwC|)&mbW!$dK5U5vi4MSjiVd~Ckx9dD4K z#qmDYWDME17IOC9T3eW{z`2@zHS~@^ul>c83)^-!VZgPd>WT{t;Y5-lm`_Qay-?p8 zvKw%AP)3GBl)0fV6CC)!p<)s95@|4D__5^K$=S0tP!C1&@4Q+gE|sd*!-GwYt+otG zZ)lbHhX4g0?ZE~6+pqQ<^S_B#i-ipiAXk3+SQ>bWVzJ(PI5O)rfRHfMz=*Z zM{#dXO8YVFP4R_@s>6Ef`1+EVik9v(lvlJ4GA>2b-@#-n9)$6Pg10}*nGyN;bw)#1 zMCw4>(oXD2?c1gNbD1cVg+3OuJ4!XH-rUPu8AUO@U_^RtT?uHsjKjnt z&<84Kd}8jZ#V+5mYtdZR?m8hg_Cx)1_tPFpg36>-9^oDfQ;W7$wcBFNX5U&6m0c!6 zFD`uOp5PCcVl_MuEoV+|`U*4~(B=fhH5&(P{G;b>-8veju`JDp_lIyR9jy{^jw+WF zqFl8S4|g2)P*kMlILYgeP4I8jjDcvHEFVpNh>+nEQw8-zI92qncj7upLQ^@99_TTH zSfA2{4A+uHr-)@UJn%QO)&a|<$P2T?n9PP;8@X3y=?P}}8XMs10wqGL?MkmfM9^QF z42yu|Pv3QVJ_&;E5RGf#rq-Ab?MeZ0gz?h|cJsxr~{W#-ry z?>rGJEK{}*qpSH{oENzqOXo`sm8cBbqKE8H02qRdLajIcVc?9XO9kkB4M#7Q0$ZP+K;LcuyUR0IwIu9!N7E43qk( z7P_>iTT?Nm`Z+O3pKigq8!6;q7^p{aBAYSexpGsiS+Zf(v^iZkn%9*#^ zFHR^6i%W_x9#E#kYWl6v<>rYoJ8{%Dp=N{8Tfm9!p=^u*>R~quNrzk;S{AOI+F+H^ zp}GVwQ}oP%2(LxQ>GhGB;knP0Jon^^AvYEK`ntL+1LHHLxOkJ}GlTEn*5wc0zm;kLiobH0>0su!T5kc&?sgfKFC>o z5dUCfX=ATot7l*&Ze(p_Z)PAbE(`*3MMf%VnGz^r3Zynx(^HHE%;e?8p**vHt+AA! z3|!iMHmbH-VTiY^c^@-)_%q#%itX29jv1{0=6kfxg7x@yutp~}sd@2HIOC5RxgXF0 zG`AD?;#pJJ@CS`-)y?q@7Olf8cfZq4ymwX|j#m+58W=PZS1BPG@X8psr#GwnfA{q8 zICgqXAbU+P(nlrGB~U@&|LSUAIyxcgTb;kIzX6*Zl*x%hg#HX#DB&6++cW|_zv->g z@`Zl?FaSTVe+%EkYneu}iP2pPQo+eUJ*^r~j%5wi5oMq$ig@l>4u#zzzgCBabvRv) z;$tr;r6n5{9porZ`9wb58t~v#C;HQ(XE{JNx|+oca}BEDI*HCcaU%pRQ)z3(F|?e? zthyHsY}9B5e?|_{YHkU#M>rGcF;MboA=q^+EaIH(pwFL4HbGkpaCPeB8|ph2d8~b_ zlZN2rEuV(;8zL8RjAwp=JJNw+4C@W$bFi#j&wlW3L>a-W$ws3`{DS|Jp+GL9lBHgC_7637<$Flt zVUSm&Ycrh5D~-pkn?b=ii@}>XVk)0st6KsaV-A`2nRaWw16u1$_B!@DcEk>u;^s4M zC+{qd%d;yHew!tIUda?zUtr{;F%07G>h3a>mfd|590;$+?C#82{5ld_tNB%bg^gjc zSBWonp~IndN%h=@prTmfD{9ZYUngu(pR9icl1mEWk_A~cy9XX!Uh2rz!dDV4ZvLsZ~J_BHesEG&DVGduDU#%yg z5`jF|Yk61d9{n{)uGxU@P8C_a@nYhGD-Bm`%2m*UXR9{&;mO_W0%~Rr)bJe)1*UXG z>@}Dh1`-69k*cJzj0}hl_!$la5n={{e#ss$sRQ^20zrO<{QC(M=sVQkpCNo-etC-a zOCJOh0!fOzS9FCqOow$<6vYmfb@!X#g^pQmzFZCWhY@|B6RH?Xri2ColOq{VMJ6Lx z@kBS(OhWSG^2a+;G-=i#Z=E7Hx7@R#ALRU)z)M-hu;KArZR!T zkQ74ybHp&gHROdWSH=FHPlX`AfRI6}aR2ASm?*djkc7AnN2=8STtEo&9mo|m1_b}# zLx_h=2#SQzYtfSKe?A4lPk`3>|M%WPkb9vZWM9Wiqw2!{b0;Av-<$uvAP7E(9P-$# zwKOX$%d9oblK9VMa}Xh!n3y)Vw-cC_hQ#80`QdZSpw2IJ4V z0^pd*!^-Mo9izT{k>77>!n5c?aL1Q1d0OT8vDzRRap^ImhJVBj*&5l0cfWaK`fjdi zRqT;}Q=l?r_2y(PX*ypP3ONW`f&>P7ss;kYlt8bN2>i+wY9}@ zJN(XJi#H3iPHgSAy-ssQy&Ows%h!5_Un*X# zxBiK-Lfs(F$IpMH6gfKF>T;6d>grmsmNSZ2ui)3&+37#_u`dF*cYJ*OQ|&^V>PS&= zV4y^u<)WodL%LK%M8wyF`SS4Q=H`2ur9!-BOne+18n=tzwbF~%beDaR_!#KuS*ihz zjJcYP_H_(VMzKOol)5FPvyKn!l^g%(wd>*ob7|QfoBQ9X1ZKgP6qt6Fo_P^UNndD* z{~R40(59e;U~*|v2n>?{9v{%YL+FNpquf1SjLT$R zMP~4ONOcDzUqyRT|Mv~qU-yWS?9S>3-Jp4xUIwo7J*j^QQuURjavM(q%5i{10@kE3Vrm9C+R=_99Kfnx_Ita5fj z$5KxyVu;|lN+Ccu!~<(SU`FP)+aLRvVBlR&wfp#}RT+lyU5v0*wq-utoD4TP9q1p? zE0=t`y}w70pI`oP=@WQOvz%eH6duCwB(^9vvn)Xagi!fa!#R$%z`+>X+)c*P7UkGLICuSC{vJI)f135fX515f z5Ekwk9hG^a%H-M2;hiY#KABuz)?m>v!}N=V`5Adh&}o?_fET%z6APUUR;xBxIA4J;(=2nl>ilgD4d9g%$!?gu>G67QxkxsWtM z^YioF3sGO)20ScvcX!L(!0|}1fOJi9k5+@xbt%$2@${p?uu->eIAKoTu94gvAEFve z#&VGC$&o3zHk%CzDb^XF)8+lBGiGz(<9wsh72;* zXm++LR4Ps+NkPTJFdEO`Ulp^Z{=2S)_R$(jOXnY~9{cl#4rEdpjs@zF+w|zt@PP;5 zAYGq%1yT>)Wo~rXsgpENgZum6e68S1PJYUgeK3eSO*olSJD)E%DcEwZ%k7+!=)i$a zb4K=kd-_@WPdbxmhmK;B3F_mUa!AH#I&{h)Mfk`5eB90 z(EvI(%LAKJRCsf$?nZfdxb%HIsOEytvjTj&l9B&qh`*f&_LOWO$kH8oF+Otx}o1j8d*`jNkTRrrxmIY^q&ITnilRo zW`=9$o5z-BXStl~#_(^x#Jc{fH?3JCzYN0pZ1c2QH*eEFx4l1Z2^2>M+btWrhl00x zLG3|!oN=4XHP!cO-yPe|}q4w+jqF0_w)Ch!}Z z_(=oVYG;GLcT*rJ9&!CReF86D+2(ZEKrn^<+|(=`8QrIT#eeW>xYgqZ>55>|?IIm| z3d--D4F~{p5@TuHR`*()#2h@%0s;a6PBczUv>NZ47!>Mf^0+>{F6qlN%C5FMT4L{{ zX)0|OPo=|?)tlfIk0Ow|aYf0K7630pb)7t8vgxN^OphwA^0DNI9JHNs{k6b3wn)T}S*~quaf8*RN z`em`-q*saljh?rDKZ#{xVEt5P@1C#rB^sm3H9Q$*DIhRds33~dFV>;jO~#41y&oI) zobJwzwu`@U+GrQ5Rp-@MEh`r(7IujW^|f#ou>X_fh483t4=&B_41QbAN?~)Vn7uVI z#biOrse%4Q1 zTtHQ-UK@w4&86W+xBn|*xA;Y7^!vN%!k z^otPEeTbogoY`r|kb%iAhEuJ9!;$3Yaoqhbb_GOavhI!`DYAOKH6yRf@z*u)2btb6 z*BSYEV^y1k7N@5ENf-hU;l`j%eJ+RJ5hNGJnvDRg>=1=R z_ow+2DngKB;<^t5lNeZ+LbS^Fv)A$!-uycXad^P^77~Q2kU%*;Az+`uI7*{~xl&o$ zI-u2f0qX$Qi$N%Y6B={kdz+5EaN#$|~4G*q}ao(R_auzW}I5J2C!hE%* zu|&OA3YZ;AZEbDYNjinRh|68{U*e&)kU%%+vrq8fOWE@I?e3Z`^6AD%IDY(D8;l$~ z0)xn-jXoy-4p~}MC=fYV^)Qv!#Tu8*I0D+dE%<-l z5T&i>?7Y3(8nG0-AQ%!7$`@$bV+7z;EkO*b6+%Xr>4HAqaj~97%B4kdKAlFUi&aK@ z@eCR=e<0c?VG_*;pBAgzc(fH6Dp^ca9Ubb)D084FbEWw9zUmMJKa%e*?2M+)lM(8; zZNEuiK^pA%?@KNeU%wetpxixwlPmPR#GZS2P@bE?T9OQe;4HkVAt=xc{)2x&K(0s_ z1{VF{(Gg=Zt5F<2r?sjl^2@UQ_Z3$#`aJ+`R9j*6{8nS=yG-`(x;GW?fM0f_15tc> zd8vr#Q*X6Qr7B_=M!@U*F$hXjgaqL)QUd_qLQ&(h2p^JCLXuB)^;WSmyA!5fDG&&s zK6H>lftKdVb#%(KTj_kjPxYchFDaF!YDqzX)RhT5)0A zFF;17fM5UHt)22}PXIneOvcM(yfi>)%ixzbu~eP?&9(^uQ|))g2*W}{W8$f0hg{aY zsBPMx>RRr8r()pY(Y@g65+HM#v{~~08EFmzMBr0P71+0kSK+-+;aNuFKa(v4i3f4B zB_|ZJs)C5mVKK*m%VGI#Po9&LGb|$Fo!3|ffAhyw4w)wBBlD&Hmj(Rq(T@h_b1PRD zNnlk}az_2S=R;~j{y!iQj$#Iu&rSc2_2DGuVq9OT1iC_&ETtmF*tt>- zX@FCZ01EuSxrpQq++n7`2Y)7*f|QH-lCg@+!c)QusTWW?$O^g)RxHz0JN@y`)nCTK zh3)sfL+I`)p7yPjnVw+km@5I~Oyf)r2k^EFB|0|QL4^g`MR=1T`?JN=At;1WR!XXZ ze_jF?qhETq;Mdc2v|OF!c^(PB6!tV_6*@5>=y{-#^gw(wp16W{Tk$ibJ?ej_g#^0H zJiss{iAiCp$;|`a6HNOu?Py89q)5JZ7dwb^#BzTI0zhf1ifw*IA#ufzU1N9wuOxo^ z!@L$q7Jpgqj1c*pAQU1#8jU*3SO}k#jyIF%`%FUdhkkP&h@~C zq*m82;6{62eU-v}347^)DoZ7kLgDproj|~CPoZx0y>8d66R0a=lil&F%`VJy&tJWT zR-q(RRlhAX9e&+UT)2XpMfe$+!-laabP5mzSq_SdUsu|FbTl<A-!92fW8T$qK(Pwt<%pVHxxAXczkko{>|Q_aYg2fac%daKpg-zZX>`aH@`F5I zaTp&ALd3>+qKb(*OtonXbNW=N->bL^iGy*C_C%b;KU8Oue7&hBo2~}}P*{o1+aG-H zmjj{Z^w}$>UoN;e8*D}j#?hrSyr|nW>0D}~I4PtDG30l*u?e3?BIb~yz8~{b>Vt_l z+4KuUZjP)9Akg}MhLwat%Vz=!+qYDserP0_4JJK>bJexGEUXQMfQuxR(aNj$Nrz)A9FW*Fe!D16 z7(gaUOC)4@ST1cRX6`gi(t;kga=rzqe z*-WNfsf0MuaIEB$-yDCl8!x2GcxIB|&WO^Hvzt61mtV5l{p2G9#A`S&1@)g4|Z&6Bc-#-YM>*x~SU?hy;ETfe#AeVbQw&Zth3 zAsU_{;6UyKr`u5jvLx z&4!A{>t&n5&+G?(Yj-}Gt^O5WW6KF@ZhWQ>MB!W4?;ozhC+7CG+)LvnL_X%fwK%D# zyeKZ`QiN5;67&N&z=0$kZ(MMv>$U!?4(0Pjlgz~;nnVK=NUU+tf*LB8>|gFk7s3yY z5sBUt>}1BbXJ9u^7R4>~b0(tm66$Fm%P_&UH~#Tt*M6rxqKTg-jc z@xJR-ynB9ncwzADfi|EP6}^Im9E81ZxmZ;Q7tRHHeladsToAPoTs->a4Sfaee3_QQ zye5O>WtOt!acy|&L9o9RSzQK8w;j)E-Qik=xj>%xXi9`$j@Ljv zC@gy~PSmLe2p#!IaECo95YPedwF_>cKjT-@xE?Ll(!us4GPiseZprq<0W~`x#VMDn zCl7t-C+0S_O{auE-7JoOV4w!b7kg#C_&;k+KB1!wwO?}Ki6r;KV3Kx>|J6Rv$4m)E zS`(p=D6$RoqA>>bg-HgK3@=+%gY*`%i>q(nB!!HR1IeKDww|YD1>aD(2hyWzT95p6 z+~jye{xB}^P*Qji{ma|i+aKleq#aCj;Mnj?LBSHoqTmPZ+0@}G97r|*REo==_6Ku} z?E7HgoFcD@H+Ok7*W|#%d0VQn{z<9}5gnuP)Y~Vb;Yx008AMZtL-EO=RHZnH_Yf&- zKDWKn`$*;I+Y~25wZyBcdmd z2E<}_fDGZ)$>R$M0^7oh6}s=8nzGM7B_AZN)|r8F>TS}F{^gTcDEt^d1WRB3D&U=6 z-J@6gWp`7d(rx`IlN8naPU%hg<8%9yx~7-sp;Ed+(lS6>D9(?!kf(oTwe9w6=f~Ou z=#NimV_uK9XG8oZ>mB*o;+;aF;~kw>c9B?|ob|57{RM)?loiSynlTn63=AJ%REd>1 zU_9h3A{9b78z0*5_vzg(Mma~BH|o;i_>z+(b0EE7xU0_8={Ghu%!kYdVyNxHOBuuG z|Fp_N?qs&s7g*fQ)&b^K09z^Va`}K_;T@W2%aFJR_HW$0THsvKld2#qhWp-)?DLT6 zeJ>eV(LGL#9xxFrXuIx9r@%>dxyuJf{spFa_hHW;-K~UD;!^}zK38h6&FkTkY=+mE z|A?w)UxjLZf2ssn=)Aevhf}DtNiiovQuGv?Yf83u zaEVPjr5caf0f;7;WW7|>q^nDv;T|Bp%K2PRg>a!%nSEZ3V20}#wVLEojjmOf|L}6bz6!i$C zOfVpgrp&>R6rllGf`Jvp$H%wTPwexRQaZ6<4_DQxVy-WWkiJ-@T$%IDW0!3_$>mIu zQt{0GrHYQ1()xH1^F!lKh61N-5|h^9GFjrMzvBz-$m-P@ip=g;46jO{VN!ayw3j$g zYyopN^M~E&WVLKb=6i!tAguGH&>CSk1keXz!`)&_v7cB2(Su0n(6BIEj+T_DWHU<6 zVRRjdvbZYeKpXrUJkU@O=m87~FN}@TF3=QHs|rpu*cr_3oh;Xv)9H5w!qaWhM0Lu* zBL;HXtTNcrPhDWDs$-8!3qW8W!4jpX4eRE1qJ4X75S+&E!7(X)MprB*6n-LKZ;_&z`jd{y0)5#a4d)oVA&vasMzLzEp64dW9!>T6-XMq~d+AEmh;SN-EQEHi+3MU|&--JjN8 zEfeW2Q|X&W;{FJz4BRu(>-V2v_?j!3d^fou2cL~LrG*(51-$NQpg>qZo`vVQS144k zB1arUPt0?CFrvE0QvSdh5;MpU!Bms|M{M(%4!v~nx>oz()#O{d-&sa&qKFLpkJLZO zGGf)CS!lvlpW!7v(=Wa{K3=)Az+t(ZRd^$v49I=7fpDkxxVdU)<t=1|{#z|VB4nhn?M-~Yo2K$p(pUR8#ZknR!SPLCkR2tzTb%iA%$64peC<{9a~zjylh^Nz zgwrmKDj|42CAgX{>bwvFW~sijuR4G%X7qa`dGq9?T%KPAfJeD;eFBD>IsbW3)la(p z1#dB)NmRAun!4%bzdJJyGuLg7hy$E$0#WNnmEU-R%B&pLhW8`L7$68kG_l~@9*RPT zAKAp^rm`hyL=H3%#mK&Z!QursZLN0(&@tXr=yh-X9!+I6ogiVSpD5Gfy#wgsASPjy zCjm9wf2P?N0$%g!!G^oVF?nibbXhv(V}|LxvIq?9aBdQ>#Ed}zxk0LO+y@OCAyNnj z$j|E9p*xvA+CX1v*2X*MtiSVkIV~$+rpbLaB|JD+q$CZD@%4eaH$ZJFn?Q1?Gy>R| z)$QCc`JGwZ3&uCM>KIC(-v+Lc2yhh=J7ej&KuMvqsJL^ATIw4$lb#$75G}S@uJZIq zfbohfs2=N6orHzgX5_9mLb zc?@_(3B|6}#-oOS>~4R1s)sA+-SVQD06p|C+Q?=C$quSlzGbI?nZ^hJ^;Hf#pSE8? zuZfOeb}*XGm8hL>gb)t_+4pjq@H`4P6yPBjYM7oN?oKy8bdK6!`mMSyriAb>msj@* zqLc3-Q%&Gj>tA1Q)KdC9WnQK8x*+Ad6kzdhcLu=L_2|u4U3&1U;)nHC_t5BX?ZB;I zrk!!yKhIy;o6XzoW1S8rD@a}L5=qNy96Z{dVZc8d%ZxDY?9t*1WK!3YlZh11P}OHW z6Eu1~)MaImbC`MJ(L%2h;VVaPJ-(py9jJf3BzSq^h zQuvrlUi^2qOU}6U=Z@A$k6sD?wI6e_9C1UFDV0A|*V=v75xH7^i=6uFrUSo)FY4HF0LVUy6QFA&+nBTZFaTcWU)b}fmJ>PX|h|2 z+_BhVIr^*Pw-H`x2S9v_CDD8)TbiHeY|`)jGPy)zQ9?!Gtf|qCt5s`0y9))2B4WYB zJqb*6^fAzmTE?VpNoOceCRN0v;Y&1DegmMsOqFVInyob0F%*Mw8RORh zheCqI1VDS{TFjTd>-bxke!jjn=woPfdG_Co$(J%JoY|AML^4LJ|jnj9oq8^ zV5Mg3rS@+D3TS3&l*4Li;B=!~=`@j6A+G@4O~hd&kui13&8t!|p z;u@TfQs~%7NJu8D+$Nr(OE^HmD1>x7? z`jQt4H^yG@~ zr51-9F{gDJ&<2aI!<^4)*bD+(!uI+0uz4cHXFr-kX`q4$)e|VFc+LHVn|)2^zUjz1+{y>ES9#_Xaf%4XM1Wzjd$t7v?+_c zi5zhuJ&OoDX{@Dt6ECfVMn(sD6E)ll|*ByzF*tu5Ab%pE9Yn~U8 z?<^POluJ}eduu#UnWVqYPyz}CNRY|P$j4qW7!}7_?!V&5$8SGnPEaD?Q?nonxeY4} z=QItENcDEcB!za$@|x_C@X3;wYHP)n_VPwJf~xb)H8z@~->dQ;-f#AX=K@r96(FeQ zs>t`-$;-*bd>Jk0gho$Y&$F^91R8801iKXv0eut(eOW4r$!e)Ta)wm5vYM2E+p5>C zO6N7zM~U>r%D+&-Rv5Bhxh!tjhDOr6oO#aT0QF3(JYtD0wSXY zpH+h}*qP0*_GgD+;Jl=y_Ax44O^L!q0HrnF(kdK_G3C(072+?|C<>{kUiQVRdySJF zxlE8~X+4FXqGFsE!+p1cZi(MV;0$b>+rX#IU><$1n)d|%l|fA&7xZG}1o)RPSa}K@ zmJ8pj4x7vo2jY3kEJ6p~113&74nZtTavs^WHm|F04mW@cYV85~bbhmp@^=nSD(LYEJqR7qGPVFQ-MeA2M6*$zouTlT`+*)`%UK^j-AC#b zyS9g8UL(HUeyg*x?GSl)6nu`4+8#&jwK%Rlq3GvB^yLdbToYr3A(R1X{BPL_T6=lf{i-qUg zl0OirN_=xUR%mn;cr1p_*cp%ari*{G&SKWb54t2;4B_5Q`u z7qk$u9CI6Gv(oTc^zl88%Zkcw_V#$D;BhVo*mG}+^6ciO-mQ^g!u^GQ>{01`A@kk~ zvpinRvOX6o=*q5#M1MMeaf1Trht>IT;g)Gf;1BYZ;DqKt#lx7sn^)HqmF`GT+*!k8 zK~F*V+1pl;Ne7RVyV#lnAEmoRIzexWorh^i$@imwnMP0d+3g+kS+SwM64bB)F{o9E z09MotQ^c5(AMzO4_+eNbf9<9*uizdZ4tTgP4%}xkPAt`OJNtiGtxv~#h%jn^(x3n|JC!g(mPmfHDGF|{`s+D`;KT?fOt{TRd2l1|W=I2RhzA-K@Q6z32A_<-E2H-Jf&Ad!-+IKE9w> zoDxQzwib0ZqqfZDy=Y&hR~=A{5xAbumFJpCFWksR_uGxW3RKb}U8EOydr$k<8tmsN zy53eN=P({ft>&#Ul1y91fzo>?p<#o4YFtk}*cNOt*rvU8=+V?PxuA+60eDA{Kq75} zvE(6XJWfN3gh(Jm*AH8-wGglvXsAD>O$eU~wdsT6(5o zB#U>eDo2XMXI)H^fv;x4kVu6tVFmd@hI^pHn1Y|bSq6WT_tzzc-k^%j!#jZfE1<`$ z!add-I8 zUe5yv{)Msv@?AWc=$8SYlbGq!F$9Ux$mJB}{pnQzUJ5NYsC)Hv45Sc$;Kd^KSaW4x z_f^B;_fTtsPBEDlAAb+Hd3H17iyMmNB1fVq&{OHW^ixx!{`TS7o8Xi$Lu?*- zx5_?q|GIf7VsIdeFF@brB92p|PgM%#MPrVT!~Loa2Y(ntV-%QvgCj0&}-6q@wlDf;rTK=A^yYDX`Hsw?+p<=$!nj8IWfEb*BX1Fxr|qGFud zNi5A=Z!e|t;$w_EbV&w;goF_%V-o}tL<7`ceT1YHeSya^C-5h#{3z4i6ry~N8(S)A z=$vLaq_&q98Dx18YiqG$V#|eKwb~smEgWHsFZsiEW@NE~JW{z68RYH)V2Zeys@|E- zvOwR;bvl?cFU>D6FIOy9mf^5jC6In7PN}?~l4h24W;_uER)Ep&L>8TkVoCfhsS;d{ z=_?ST7b+klL#TgaQw=$T*9@XlF2PwNJ^o7wyc1oYI_9^x?3J3Oa?v3+<1k2(S$1(6 zV1d80il&>E8KyY|KbhNIUDufr4l*8owv9=z}%fN&l516s#A6GrQg>I^{8%Xbv(=?A_S3`CC$;8$Nm!|TQ3 zIPXp64uMoF^@B)kFnPc*2$*t4w|94vKs%^}mJFz9f=Wy-Tv0$@FI2AGs<_kcadTW0 zkSW(7I|9m4+~rJ^_iy>v^VzWq{rqT|q3J^q$7XS+FYhp*o2rECz7?wJ)?8oMqWP;| zst;>aCX=wy?rFr1yUsd6^U3+#YxLW)Qc}?)$2y<4(P6PfQE(D^VT+@Ly7qpK(WU>-dzlr zw?ESXUD}LG>T2xao2%V)9|EBrBJLl+){){djW*9lJeoOP@B(}i=~l##9D$eYrRtb3 zo2jfI3|O)2H0X7;yn()^%Im>7Y~qo~E*lw8V`7F<{%+?P69ysrxJ2Pv$ipJ~{A_Lc zn5j^OUg%W7)OJ}U62hn0dt2+FZEsqQ7?nuhk20{4nKWEG9b%w#w7E|eRLuExYII4u zOy5WeQ|NL0Q}Ovx6Tkh%rxCp9uQip5laNyd#*Kz=DxHti_vXtpBM7+VSvFB@BK)9W z5in|i4nQf6$NOZJf4k{to(w)nH>vUQ z`QgNVqYDO?)i87^GV0l)yfx1l#6oX-@mt1hxvsPb`#_WyfaG*gB!B^apX&XstjADK z*lS8YK21OsrviW(q+m7Buot^{%neul$y_Ad&@pH8HK#Yd&hj-z;ai#+>y493WFx+x zY@}`%9rg*vF?uP%=aeXsc^I~?nFRN17>#bnh7A1tUTjif!Q;hIKEamPN2onTTnU1n zRj1ohUP_0)OGfaYkDatXnlJC z-~Ia5WrOPS^t8f_{uvPG-h6VL{2?k|q)4kI5CNf@3&*~Ha$iqs@ zDI=N9VpZj9`+ zzoG?|R|)J)(BhP6Wk0=r3F#RgpJz|A$sa-IgdPW_^*P_`CM4^HMr-|*!VCuAG(41O zx2iMJOqTZ+Ilz@@{q1;iIQ%56bj?4kJyWp*PUJOXM4td}HTohYim7ReEpgT=;Y*xz^N?q<%q9T-|D6kbha+ExH;SbTIP)UXv2G* ziIDKr*9gI|$MBC>FDu*}3KT9gamTDtS29zqGG1d#VePO#Ury|zQYDCr!y=@Xq^kG4 z>T%NR@dR*D(8t2RF7Lq{;WJKUPUeNw0hK>#HJX?&e7Wo6lIyBy_N^ zH4=`b3;8pV)f0i6?cJBhw`O1VEQ-1H(gM3GJS}2Nur;C3PTop=8G7(a;6!2u(5b65Yvm%`+SRG7)vA|r z4jY0rc^aevwe%VFWc3}P{mpq|+^BXehK%+U87?Jw721LY_MnCfN*6K#0e|;Fmjv_q z;rP@1^r?VcNgCwE_?q%oFqW_rXqnY-(LNK>y!}C-X2a;pmb?f!qVbx@=G>J{Hymy`^Bl$5W z7wdW&%V_P^@Jrz1WjNh@@3J&2lQvB<0;LaOIvq)@;0+}6P}_E<)H*zypwH+P*7H+> z&M|4*lCX{i3|pqFvw`nN;=9^925syYasV(*(gB_#g|MbmFiZv(fSVMlshl*Kk^8#w zZCbL7Wq@sqrNT~qM=EF0X3-6r<<(z|JRENi}n8nv?M`wmCvwxEkGJ(y+V=}08hcyht=nC$kluQY5f znQ5ZKQjud&l>KaPC)Gbn*>>;2?QZaFcn5TjHw2r=_fL{EDEpibQUja|hBqaZE%q0B zoTqw`^}IlK$cZ=b2Ct>?EF{(V`s-e5=6A>+pw74eRGd8D47cC7Va&`Rv8@uj`1jY1 zM#G7NfPuwI9Vr4m4e;O5Sj^YqsN(%8o2p{E6l9GV2$Od&X zBY*Z0csOL+16*|HN&)YE3pFXe{uRx9v*A}};&6`A@43bRSY49>PR7EeHd2$xewUOS z&Rml;u>#1kk|d z00se6pj!C@pv|ZUQ#Le=BS|K)vGPbn**E)y0f&1oR*b4f9=rP~-X9iZ0=}dl@+S-9 zR5X>6`*qACd(lBiiFf;>vg}I-+qC)*n|6(FSNYk0RNCj{bN0(ZXWdJp1hdIEFQW}A zt<)Wh<{>m})Ijgg6q2hF?-j)Qr^$4@kc*c_>cV99&P@0VzK_FFzCAdVy@&vtS@sOPp4$C|6ivJdVb zTC{f_7gP9Lc1$E+d6X*OCPV#ZZrafoe=is`Cr|brA?SrQADPti>!-C+#|)#WK=}?> zjo^R7u(XfGSvfv3CYx2{Fl2oA<1&~-gPXJ;wZvszaWn&bOg+5PSkopJ97L9jbn{{@ zHcjtNRu(rkp45O3Ihlx`K3p?AiwGa{!5<}8o8lHu8w45-sw^z&yAi_22;q^%09JKG z`;AOu&etqmXL^Dd3OKZ{-0|+U7yr2Z@i7*790U1UaDRy@yzKItr|n~(CGNZ30{!`0 z084Z;Kh{vxabw@AeW||Qqg7kmx3;p|ZwjptKi}cJ*qWImmf3D_b1kaKFw$tDV^SZG z_y;!mz6hUyiB@qFKF8l00B=%Kv7-3MozDQ+89ZF9R$4JK+^H4+Cg<(-M_XNog9IO~ zP>LrqC@3|WRo+#ea#EfHH*cz?JSoAYrIojSx+ylosgZm}e;XdJrTH^QznMvX z8D|JXV@V6<(Yq(YMzS@vG3a zP9&s*F})O8zCllMKokeRI{w+5kM)>Hufp`6X8pd(W&k?ma%o(AAB2YT0;kBZgY@?#+F>xKV$L9ftiwQ(?;wCv~Y& zPS`IY_2p*7oDw9t&_90s*d0oj5B=oYs>B8Y2urLMz<~QzH@>UJY#P{KP=wS8pV0&L zclvUm=dO;aW23viIw#b?G{GFGos8Y@jwI)^J{cd)%@!pRYk(t5GH=U(!5wdtUwQbD zT)DT6kF_g_xP8`yggbxnISziMjR$s&OTEY{7r#gCWZAreRp914`j8PW00~(d(3mpZ zVPIiX%eKY#(JIVefr|7un6Yj*U%4GO?&zCE*Fpfky?71vg&MI)yk4F@)W45+2#vz5 zX8z^wB$Yn@T?Xr*c>;NW|CW4n66O9ezf!PGfP0)t-=$#^Hr6>nn!lA~B6=Kkv{DBD zm^w{Ai{p-NH`6h<771FxN=nTaWyO}1VOMoQfFsHA0Hk9I;BNz*1Fg3-;EGh##lC#m zKuxVwtQ@Jt!*c-uuB6@FUDb!P7fY^uouyhK;JyBTOr3R9mFf5PDUlHAl9cXL5Cjh0 z-AGFcNF&`H(kb0t(v5UTiqZ|z-TmHY<~zUjt~Gx-GcFG2e(q=QYhTxAml@gO7D6j5 zG#_v+ma<-g{~e#{v|Oj@yYa|3H8EJf1Of~Ua$tP)9(pIa`1kvl|0LoxLC_0Se|+T5 zTTrH z=FUiQh43?ZWlS5Z@yEa48+e(P zUm`qSbvzA%Sp~uo&`IKw;7D8u5Ww5sa|nh25Z9cqsd$e^w}yy|NW)5c}JXIV(a};N!cR7-!AR^|2M0M5HvW;dX43F z7pYt)8E&e4c`;q4e%IF+MYZh>Dc#fcXzhLSUJh-+Bd6rCzr2z=y^AfzD(<}rZu~^@ z_(3;VvYMAURz>#v*w;uGsu+zZ8Pyfc&=-6(G1KyrQ%;{9_4hVtsabw#rYiCl|8iMC zfdl6&4(R#CzJI`6s?`5q6wMI0QelKmGNuQBkfO3d)TQ)^CiH85adL_xAYFeMeqQTBtKi zy)?Y%?V|S$3KXpS2?H|EYvrGmCc02lcI7R)sf}3dM$+`VTE9!xre1{S$Cc35MFnRr zQB9y5IR}(a48FbNTstJ=?6Tgu0PUk1%O4N~P_=%0+(!1!HtR#uY@MZs)b2aKMq)k} zMKJWt<$Zt}h`|h55*S;jnn!`*gC1V+sCLio;r}$24Jn*5@{u|)6NBa&Yv6C2Zr{J~ z0j6L&?i_HWN(3ktLZIcyAjlMXy;@cY6%2a;qq_wX6D7|mZq3^Nmt6LlKne&BK*D?6 zRrcN|7(D&SdE0N?7Hq+}Vfzr%Z#^*^{4$YgFj)$z&0C}1`_fnEz2}T-rT1pDa^*DL z$@!qcT+$>6rSwy+Win&gnYM@`nk(JC!O!gJ>7^}g{L4UJ!4qa#`O#0l%<_-$+P7YP zodqMrN8va4CEG~DXu&Skj$e>lcdQN}Vn^Sx=~k*c8E zREmdtEc4I(gG5Tsy}1{*@=bn2ZI52M%Dylx+2g>|rOfr9KK0s^cxbyfwNQc%#|B0N zoLUdFKDFJ5w%)>NJg!E^0&>O}{3WfAIp41WITgxlKxu$_n~O0$+WzK6?rp7p&~0sN zQ!_}q?y83XljmIhXGS2}sJW&)9=$BP(5bUX7#xvc9lmb^^EKrTFVCHYrkc60?FJyf zng`rmr$E0camx9-M# zPgX6zqgTuhqE4c9kr(Y0LWbPtiW2KHPW<3r;~%3Ik+)<0pud1Jjd#^mH#OkZN?UN2 zp7FLeg!>1#itLZK&ZpCoSw~s5Q?P-^Uy&1f}4vZ3vf2kkS)W{BFAFd?0cX#t^`<79GIG)yIJb8?>tQ zU=24g>4IXx^nnXT;Ccl<&ITViZ<^Nh526~Y{#VLj*{TDNR-xzOeB0ZG1j+n0C^ewn zqFRLHL&+VOok$-ie%^JQmi1CTEyjrtSIOmhDvn^^m*ZfAk%234kM`ni=MG zlXv^+HXcE3UZ=o+J_v49H8VHhAu74fi(IF=KHAH@xsDjmw1uDC1lEVi|H(4OAWq;{ z2iKg{?F|N}H>bm_zh#!Au7HKJsJGku{v?=#7=xt9;D=v<3I-M8q4Yk{x9GgNAT>l7 zkLJgREB&3{rAY>1!T*WS;?%YVmu6VkX0m^7Wj`d3B9eWrDWP{xWQm)v5eRSk^IHtX zqfWL}nwEcv=#owI%Lo}OZ1B~$U1rT^kB_;Vf79IyJQJ-%C4yuE>LuUquV`SozPy1W zWff-A-KVe)v zD)Dm=&v=u5JIvZ4>EhD(iR%YwCSKD1BPyWHu2wat{(qlVfoOH%=0`Big*ISh@YTu7 z*i(bl1M(Yf!8938iu(+2mcfsQ`YlY_W)acuzE3Hm<=4U*#pl8e{g01{gfL21@M?vS z`;{RVT@RhQ3C}1L)^LfY`u+1%Ng|tEs0L+q-L} zT>bI7ByW;LFXm&pHg(I?@o>`SWcUH~ty09Y7fpE@F=en!03P!5-0yDp4)>L!J$Sls z3dBOtr0%ix1MVG8;CuxCisqs}KI1x~l~;S!rK z*Ah!r%Z%;~WZ{mlBqHlzEI#_AKqt) zsZpFNfGac4+7w~s2vi=e`z ztUoUsK_hh_DPQbcQ>J>IA7Vc-SYlEXS$XtKz;-p%ojBYTL}TfJ_({jc#;Ee6>+OjbrlNe{jHQ;S5Rta8Deb78lr70E)L!4m*D8=( zQ)M|}li`^Sses@#xAhuwjZ34OC{nHruFnb6_ng8ai1U8igMwAZxq za~oX<&1iOZ)E&@414{{Y0sPE_K&1HD)yg29v3#?Uw3<7c@mb<-t}fFJ*nG30FSl=B z$(fO4$DDIb&xW{(cmp%5?nTn*((JNRyZA5dt&+V;?=OFy9dxTbz4o|xnymtrXXS{V zeXcV&RJnRDzJ$4Z3?%6iPSOTXmbsUeafLm5%baUYtB}F@8??dln@1y6HGPavUL$Rt z;oqJ9N+LJoVbsj`sxh=pYTsLOschzUqvOjp1{%)mWSu?OIyYtJ7wAmt4g5)^zYp4& z>ay8lQPJiw9?RCkVHXd&aVl;@Ya5v@tQRQLkKf}0CKL8kg{4tBAo`P6o4yvMi&bX& ztsZNH>uS*qV;F|_s}7J_e1qy>*z6-HnwRST8jfA=@$Ol|?rSBL4+k9*@I+K{8tgS& zsb!23S;9sXaC=2<3lfSs{@LFDQMg2xv-i1Y3+aB6)iV8}j^8cSR)dsJu&b`D0uiiCVOD6l`zEA!cJlSSLt>rv6e*58+t&`5wv%aA zc0o53D{-h>lF@%h;Rn^STi*~t{=}c)4@O6l@|XDWB~OP}O@l1%S3CFLaDr*2dfHuN z&X1gLEGYtO2@&u|-i0X0+dOGK>IjGs)kbOIq_Laned0C$PClsNf%SN#?pkY-bFB-( z-Luuq7ROR;FI})3q|r3C6u1t!6mo7IU+tstYTRD)*0(@dt=ghk(Nu)SaXT&G=(%9n zq`#X_P~5DAv&tR{@YruCuI5KGykflkC!i0_C(I#Olk(Ht6^w-_iO5X-cAJ3vebW!~ z2ejd8vQ7>p>^XX1?5R?4XFd}ABz*T0%j?mPccezt?nM90sWsi(`0~Nw!h@&Xu&#JO zr*~Lp9TvN1e$494)5!$1zi?c432m+FcNE~Fa7e`Lz$23^$hi=6IFB(SCF`}qHzi;Wv`ta9+=BHT zd)7%zr^YzZ5AsBA{cTuIC)5Rpp#?7@Dct*ji-h4rVc^tZMX#{)w6M=NsVRS=Q}cOR ztViw_fwW3Zx;6N?l{AzarinwbZ*BBy7(a4h^5U`nlx_4HT^rv~u2O(syEs7d_=-D8yXNmex>&f; zm!_t$r(LzRi+6UhH?FUpZHj-kH?X8}9aP-gqZiLy%Kb9?r;M9<>^(Iv~#D| zE%Nfj;+-g8Nh=FW#axFu^S-$tQ}U#GpPA)l z{U+qz&z8y;C-1;Dc>;xQkYAhw6NuWMca7KQhl|?Hjyn_}A*_m1wUYFIm$~@LQg8U_ zbo?ipx^H7C)DrUPvBMt_^>x|4>E%n~;=1cE7j{5q8l@`QxmBO&Tm{l;2vfN^J^yeY zBBUiTXTFSE(PsChUEt%(Kb)uBOnD`Dt}^A`)1-^lY+Wn57)B_m@_wE<&Cq9I}e$P-}F-RHPCm-v8z)d0a zgBTdHBJ31OL=eQRoWIt%hoGgn|9iZYqdE&7;i=c`^73O@JLohx?*?^iMj^Z%+f4oV z{97^gTr0}Kk>5DIBBDJ1G%UyNAGM#U8i9H%{&S6l%yoA%Mur|Cr%aqL_A{@iX3@oI z@eS_Cgb-}Y3*FttHb=`7YXNZw#?z=eNIu#5vmdq(LU-re;%ta3TnInz#6|?WP>~1N z7-2|tSfiD1KdrkFgx7CR&z7Xs3jOB!uxe5K3dD7lxnVi(KGv>`a4)Ow!JG$5BqnyB z{r7YzA;^1OwWU2cwUvxgf=qntX)OH^IxiOxE3l~QBP!3xvYUqQukORKoQ~R?ml_N` z2j4qis+45T)+lkw5q;dfHyhFwoQu#h*uAHO8_)F&-}teBSk3zcG=lJ>)Y9Io-(p-x z_x);EpT)!tYwN(s=XAejjNyvxY?OpTC4JIV989{XpflGR#CtNbB;(s_kHZBXLZAw8 zAY#gu05!V0Yr*S$LV}@%mE*q8v(9>P4q*q(ssBMLVblQ{1E!9I-Aky0^k*`=2#BUo zma3P7uj5`Bn6$I?2$(y+`m%k+&~*K;`;+_(m;FKK5g4;Miet7$xT(74wLQfP^ct%u zq~$&QLLf;PCtLeX$!Fnsm`_LQ;o)D>n}{D>o9-VT@C$SaPHZ@AV_i&s=3R;*k|qR0 z&|~PBhoYqTJl;Z{u6k-rqX;10zhL{VI=9H+U0r@KBe4i4!?KMT`5=k77|M~EBB09_F4D6letvZiqKN7pS7 zOG}q$LDYK#WD+K!-Ne2>x`5)*eViNA(ws#hoC6Y{3T;4q&U@htG$W`PM~)dv7r6DS z?hN_wu`3ocxC^=36k1VJ^9-8Ir&5U%x0s)v)&Rlt1oPla>E%bP??y?y%WoZM=+B3v z)+(ANx3TiS?jC7h7T?#b<9jCIR4hsEUJz3sKSaIsrUPevJ`hGJJv;UnO9u$ zI-L9@vl5X0m?Pj2j?rp0fC0ly#0}dK7(enHjPI}41D&1$nu&V&*9#$t#O_ctu?P&O zAT;$2@+C;Tw-^Lrm2H>b$R{jM)_U`SLcD8Ip%g{Rp9hUfnT-6iF{ z#3a*ywv~6Wo%z5vsgGJW+Cak+68o#6(B3F8(;X1Hun+E+?;Q&k9@d5loGwBj`7-?; ze;z4*Zq=t8_3z#%eV!G1bxu;s`-HjWJnZlKmD=h^SY>OO+5z@tr%}h@>So~+?^6LH z|B~REd>O7!I0mnmEg=Vf*=xZ8RYwd={y{#+k{j2BdREu~ZPsG3tgg0*b^`Im5$fu7 z(DAD*Gq?B8)mbv;z9eyP+}-}iI*2U-BHfU#VPQ50yi8!}&_zSZwI-Os0Bqu<@tbLA zQelYB$4*E}7jEzyY>0wd1gYcpxKNwHEVFL2Y8x;H)4D)iURBmmS`~Py2*O2@V-$b= z@O1mdru>XVMLRNSF0|oqMOY*f@6%tB*v4Cg@&~mFG!pB%bhAFk!I_+<9S8GqUFOL# zTc|VZr@%-X<+3^d5p61E2WP2b7mEdt z5;Zw>cBs^-AiiES@LHz6H2Q-i9@Mmw;kS(# zc*Ab}K2b0X4aUJ~l{Ul)?)bp(ZIQxVa%T7E|BV-_XnpyyAu-yt3dIb3Zy{G{m91qS zv#S}=np+S6s{03L_`sKFh<*r9Q+jnCL1`@At1X_3;7cHd$_6CA~`S=`V57(vw< zndN>ctP%Q86}K-@f7EbBr~)th&D{l_O*wXWUP5xUyd*_!@-=>@yHJkb$5?w4^Eu7a z(Tq>LB%y6w3S>0dS4Tyj8SRAmgINkx)11AVtIbMSUIxWGIulqfuD{91GeUNb@H&Qj z;^Mwmn+Px-12A!u?W*55-0^9YPLl1$K2oIsj4$zfFb8LCUCd@cQ@21Mn{ZkQkHYlO zjpeP#47POgReGV3-tl-Vvuxt_1ghF8iW8{b)EkDIpsnz`q}q%*kyXFK6w!bwIlRd* z=SboSgVmEuYI)4LAG<0XcIDc%O8lRX>3c3u&dVU}bub~5L&|C}GMAHLT%sr+*gEyT zMAxszcufrdm;n<?ITWpcu%#(NBeWz1l=wX;gYW-*1m~tPmu~4Vh_Q6;D4DVmn)ONFL zUr0x&kgSr|X}74VTHo^`A(;rq`EVfyU|Ga~aGuY5`FEPqgDx@^^$=S${AWBpjdbEa zM;(B2`%g_j|9c&Lp|CX=)q!$l)S6rVevd93CbIZJNijpcVed0KuvuX<*!{pBjQn6fO92@kIm4Jz82QI57oBa(GeerPMC{jw8BsrX5^UtQi+eYLri( zOr{{kmB`k#+Z*eF5#2zq;5N~SZ1H(l{5`^vPdihuxP-;yYPJ1r-lk&^Or`S{$}z1S zz=JAX2?Wx-tgI|e&X3oRj{azS{Ks8xgt}FIFoq7*Au~EqBEs8mQB~P?)F1%FFkgG) z)fgdk{uc3Xf8&8B0K5_pr`G-Jlb-SI)Ga@2dU%*9h&{VFe2q&=+tO9!C~UvpZ3mq* z?^cBu246oNm%HFGmk4CeT$kThGtlI{{<6WMrSz)ptssT>A!F2|`y{6Qh@Ps4UK+FY z4P)`%*FMAMmr>akH%(C z_c<=!a0f>?;o0GPeIwiYGK55aPa6T#xMmQ;OmSs3Qzf)nCrJAq#GejAh|fU+rIZ6h z35oK|Ew%y_DkZb9KRG&)-jy&%O2m^jPOS1e{P>2>R1u}%c6C@T?DbI`Zcm)e-N#%Q z2|tSP8Mt=<47mAT>tE`gAPO?^{XJ?^qh~Nf?+9yd{3~8f2Gycy^Wq;J2Z@X}7jM$z zj5gQ?=jy__KHPX;j>5GU?}e`yvee?xE7vx*v~YU7Amf-ZFk-WMCrjXM?5Rs~?cPo3UKbC2qoT+;8zQF&sT zJYf9L@m7s>O;U+*i8d_c%E|e7K$cYZZ;V7K>)W^`6=#vw94#P*lCggbkPq#XQQg?M;I=(+wr|3OW>(r*O@m8%}_*O0Qen{XGT(7vq z#3)&NR;t9+b|lMmD?cvUqVU91+t28=_zpzZ+NVh3HizBFbzGY04+yHw@{PBGH{4SV zvIi0nfH;x`w@E&6wno#<#oe;%J=*I6_-dGBEuYjVc4QoHZz<##Xn9`)wF%|~$2r11 zF(fBmz>r1&uGL^LdM-1RPUifq6gvCZe4}uALg>Ky&u7?4>>MI(qoa8vyOimiAU$QZ z5klqg&>v#ZpUNCaH0B?CN*q;^K!4SBgF|u>!R5_=pMT1#jH~|XPb|t}vPiYB^o=w5 zw42Ca%@&9H_hnXCk^cUk@a8TiiqtOQ>*8!SyHLYXi*6}eA68d<(&Y0vm;)n-$)MRy}c1w{0GvzE0* zcVbTt<#z8s$^it3+K}oF`5WgiPPua?o^)oN9J<)@HzlEh%W`G+dm+bR>d54oPe9X& zBfgp5ZL{l$lK^-#1HfesP@cSX{rEdShR6k^u{nn3+04p#gY;AWcjPps4O|C}5t zg*+HEauPE=V~thD3gin15=z&j?7rl*GA1e;I5^79o(KB=5~66?#r z5b0W8VzWSo1BVqNuN`i@X|fFA=+@K=WC`Vem8_19f6#3hc!pG5sJMu>Q-FoN0V8ff z*UexXN@8Qt;zg>Z(tqj&ehuIY0rdATqiJG7EmuQCj_2ojTxWl2ipK-)>!T^!U3W&~ z_6m>azuS_|z1g7*^SITlUtdk|wJ>QCCNx_6iei5ei|eT?{YIgK{ALvVWiFqg+*I;( zY-v6+UdIq{Nqxp4jX-FGd5vC{?}u~_8c#iWjozLfMd05kwpnhk%**~>>nRZR@3&Gb zM!r31GIzD}PPdTJlPmEwfxg_`rV*?86X)W=yne&1^FCcziEe1UEFlxTcXjvhu)l?g zfDB>Udg(3iEgmMlfHyU*(u&O9>4p|gV1U@9Q8gc06rK!DphTadi)uuOr;EH&^%%W6 z)0o&kY0uz;7D&M-=Yqke+#Fpd!ur}3+wpV;jiCK9ceY#{kysll6~r-A8I8L z&>C47TeEldK^zB55eXSNW?qq}0_>+jzD&78MC<|RA6c**`I12AGs z+t~MSnqFD|=kUu4f-N@&AH^C-&vo+Bf{cr`(!WKxrNB)c%TvgwW%;~Jmojg(PttCw z$N63S>>!y-w&1=;oQ&0s{nR^hAsw6xmL(n4i*&J<6_=?)uTbt6B@92hieHfht3+!a z-`YPA^1EOOp=XFI?w`5zE=9{pFWtE`D@V#^Y$Uyf-^7#A)3+H)T1-XR7?6DCU@9y1 z%Yy|xJjJ5WGYgvkT9aYCSHOeu;#!Q};`$MQYt<1EAvnc<3<%y>Obz5{3)!-BypT)b z5HlqdrljNYy#$*t#t?nPRBx!D=gH>_veh@XP*M_J3Y$);YYyjYPuSPJS?r;BLS5XV z&0|{%BAeD*yE)==new&He>B6~4Ka)dF(lsB=uv1CZmkhKx3^^T)7g6=T`OlaUHzLY zVoty5EIBf*MJ$TK`q$SiTF0|T6%?eG`4ZVN_045k2LgCCD4Hxwk*r~+F3O{_hv%Py zIi@JIxw>`lGr*~O4>H}bKD%~ckd|=ssI-Yjc^O@2IX(2o3UHt=qZK1G%1pa zr^|l7sB-O=1nU~{+aS?6$)*j1*QB&|r&zoBPr@*rLM9@b# zo8gfP>abtAulaZOCv8IY8k>#JBJy6%!VS(UMMQTm>*HM7(_73&)IHY4wUAOS*ko;I zwx9TnBJ$;(!A5O&l}p&GLYhABkR7e=%C>ju($+MzI4%#y|485ni*aW~$3Jy+dX3Ia zyYqLQV!<^na6XD7FC&}~zVVruct46pwmHgI3~Q0GoP44Q5%ye+dn&KKRXff2$4#uV zWwcOon^oQ(=m~zns=8*yA@l)8(naTGDaZq$hl~SP(ctoOCNRb(ti2b%`NVW<7$V{K^k*^_WqC)Nk*IDFv3*DE1W6 zuH}P`Z8S3(&i4&#(FZQPD?dLEr@~E&lCMv9?;Dl>y?ie%6jy=UY$!pjaoPAXZ$E{P zE>-c=GcY!DM56WMLoXaWpb@~YM%oHBKO7S-qP&db#BOY`yvmrAW0zUY90^w>LsB!= zn&gJP7qHxMX7TiX^>;bl|9U0%fE=lp7g)&y7#yHA_mc(vJ1i6_Kq{oCpb-9pXY1r7 zH~z{3v-RSSHgY{Ki%v3|rOzx9-oKaZIE}Av*TvRpl!NgWAc+@zdCwJ>gn5N^pB9Kl zAgfe5-M3tH{fJHq@u^6k%WpYp)l(?8uP4v03K{p&NubZEPFAzMnmB0qiz0qCw&*6? zUjfO%{>bCX8ITb&;i$4k!GujW{y592=ymdKY&!0=Ygr)t`AtZShdcfPH(coF*CC92 zAVncE&u^7*SJ(#KX2S(M>_q#wKv&SGfNg)XMVdgJn=wd;DJr(t#*^wR` zk+B);rNkiwrdJd`rw4$6U4A!JbB(1Bs*q7LYnQ44h_73j)PQeJC!JNTI) zu=9Vwyyws9xK5o4GuR;Wj=C0baB%n~-LSb0xy|^nOsR2om^|FG{UxT zV`6JvpW$iE6^ZmWqwlFc+qw}$(pyc2C)EN<9_~lKRubl0!Z%-B+2L_%7hFgu9#MrH zS3Y8pv|!0%e$3;e_Z8M>urXzzPh4~5P^w-K5I(s|oOQ8X0^sxVN1{Yk@3Lxi&dlVm z-B)B|Fq;16b9KIBr)^hDaW%(5Z7x@f&YR9NK%>u>fAeKu(_*fc9@OvNLBZ1Nuzf;Y zltFil&(DHXUmHwtli_>=1uO8Sr56Dz`Wu5zm@e~{YIHn^75$B!Ki0F}N5E4E=v?oy zKsOLy1EQoP4BGEQ#gLzKVU(tn@vsI5_1FtO3Ck-IL@HMwb~(L7Yb{__!c{CGK8lit zFr95<$=RPbB_=IKvd0u3Rb=Sv#QRS?=->rqyyLD9yO&KY&Uu7Ao@{QiGFieE9$y|5 znO0_qpFpl3CEO@>;RV^sN#2wZEw@K&}d+#W{eW1 zH18i2oO7e_`I6a`ZtTONY?*j6CXI_?l8f@JMyUy|J)S$&IW{R6Hu_ar+`cG>^BOnS z7}6tYdOEW*XRNE3z=AqLbRcHO^Pukl{tINeB-!WsP4*^DM9XHn&5qf?XNI}b-Vo5x@rFHXUwGWSEBsmEy-Zu0IH@9(C}izIJEts}83__!haCO%n*B|N=H@7_-OYF4x!f1zr)lKIAEl;w7>!?Kd@N$-$ls!E6S)vs!f z^2b{0^qXp;^4mOW&p-JLbtNA<^G+*UVD${)+z_CBG!H=-YyrTo{3_Er3996!Dt;iq z3GDI$)Ta`_MVGC6ng0O1s=JPhdmzqa4nY2uwNdaZhw%91`S(UdYl5$XIIGzxvhIz> zlY708G_Ahjv+AXBpz6!2#Lk4Y6`!r#`Uca`(Mci_Lh1HX^sSlyR#213z9Qj?;}?*m z-sP`IRG)9=D>q(;DZU!2@x0qnpC+YZ(3)_F8To4(CGH4e(!>0r1Eb@Ga>JMX?n)MC z@l@NMeqy`|xn?OiWGdC0tW~@e@CLh+{30^eDejjT@dxFdq_=W-sU4I)vCkd#2vvz| zeR5s@5q>3Vnio$iH9)<_E4?ui%Y;MVD2A2DZq%m;{MAU0Hz3{`<6kM+W^Do4S(ye$yuQLdn$I?Xr)1BfeMv?WaGCep(HgtzyGhItAtZ@N^$qxp9maKj{dN zy=BM@wm{8?_!C^1Ni8D?Jgvi|w;?-!kC*D73;g7k$vbF)fC5>t1awBGG5@uf+>E$2 zCZNeBh(L1D6lBAW1`7{vcS|cb{5~I_`d9)ocN=gkI?TiuWi4Tw;(RwA6ox!r=N`Tz zzlTTs1z} zm`OKw?&tMwLWU^n^f=MT2c>=gV4S^u$)x2nSp^-V4ZD{gH#(#*-Q)WiBhu_Os2WDk zM46QC+cjQsd|#bg?5{PR7W?oT=gB}xX_+?CdTglY-hPVIzv?d%Vn9ru35XK{&OW%K z1OYI~_@2fT`NcZt$fvldYJS+3w=gT+eV~kiHcplJzTS{0Khk{&xy_ss0@)YM=H#>& za8gFS18Z_h+l;rjXkm6vIWo+{D&|4g(OeS2TdMs1{6jl85ozIE)KM=XNYmsRIg0Dn z@58Jee>&V=4LE(MeSkZFRO*CG>0wdcjS8=_3LaSCP|G4=|>a>B_mqEkm0>ba$U}Z0WT5Xk_T%--oJAj%2NFb1y%uKSB zxylN&cN9^wVafc0H18tjy7kF+c&k_3F!Hqf={gMgXyUz?z_XQ;#^W$a=6x?@`_|8} z%63&)S3L_bu?HQqLqs|w;BSXxzd%P@VDx3F&Z2>kWB?-N|9cK_d>JV=x2ow`B{(fT z-=^G)9%r@l;m`hrTiO-G%;B)Dr3-Vab|rNCOLFvOohtt68_$s<7kW%VlQO!!4SOgt zu^6v*3^fCL6aHz${|izw(wV8QtdfY^fpM?OeBCV=fO z2UF=DKZ4lXD<}vg$Vqg?%p?*XA zOv(W8BUox!g`0w3FXRS?jIl>UMPGbQk`oyD6*I!@XtKh-qls6QK=W( ze=PrRm1h9AN)c^aj`a5{c5&G2rVnp&o8ckSFo-C<2YPMvR(oHB!DtQ?SFej_M<*-U(#xP%W|DY38T#K- zElW-hH|*YCD08#(Hz|(y@o}Cpc$r8bmXuXXIwqq+_=Z&~X6f*h06+VF`m#kNWzlG$ z?ZVN%Jo>uqtm2pNc{ilce7Wg;Wv%|gSh&5r=&!cDq))fCU*B9#{L+YIulvzIANSq# zw4f}dZBU}|V8t`&#o@-zp*ZPb&UUKpi5+h4`M}G(kIw@hsUxB_!|?jhbppC5bLeCd zxPpU%!U+R0-0(0((1kOGW&_z)<>QlAxUVfb`J1d}64=U(4vDpjtWNG#3pBkL^m71f zYgIw){-IG(<0IRF<9iOq94;)Z?5vmA@n?X5j3NHp5a!!lr8R1bzakOR4j-msL=0#n z?uKzPBxj3RZB_iQ@9y?Eo*V&fB67`$296!62%W=jO>}3?zg{U7IlMhf9mg*LZGIXO zG_4Zl(~BthF+bzPw;R7qwOIKSHn1)vI6m;-av<4<9j-+Tv(o3SckdAvitru;FFNV? z?VqHT9nPH}5n4Dqo8nukiCJ`&E}P3E>_HY=%_0IBkre$;9`^L77WL${HGBK~-x?km z(pxSxJGAVvn!n2Qh@iZ44zoT|Aq>DGj$M({TeVlwiKxII%2a7y?0i71c&Aw9fOJKy zk-eg=qWD~5S%a64BK8m=aQq-E*l|yuZ4qajXfRSoO-XR;_yOin_ z@=$8Ml{3t5@A}YFI{|uL%Y3}J=@%3;F<#OShJY&G6D6L|jS-U#r}h$?GXM+n{@@ad z56jF}%lN|2;b{_KZfFdrZ~q~cL$rf^&9~iBjFiuXTq5u|!n_!ZtstXL6Pe(qMrrf+ z{QDh?s>}0xgFJL!KP1H3AcPEFIn>hMIa_P)xXRX3wiXvA9$_NCy#{i6ab|My#X>Nj z&g7&O9fMpN$earaO4glIIc{BIB&Gv5(h75U#jct~@=ns~cj?lb9|arG7^{_HaK3E4 z3S4ZVUx>`X-tQ&?d~2Qx#gZ%!W@hnG8C2(2l`ksQ#XL8S8Lipgz{tKOA;e`i!Wq>? z&O1x8k!r~2#J8CYv}_o$*X@cB!!<-wWzkM93CXvZrPxSXbQ~dj@XDTJEE!$D8=Fh; z!Sq34g~y8}#V_g)@H6N`ASDdCi>V|d^j_&g>*jd(E*Y!_`n>O0Q!f-K$^_%J)uQ)i6^wZ3Ix3?DTt{?rng}FZ$bUNQD zx`?CmoRY1-=D0;Vv5ia|>S(Z$KizK+N%KIZeRp)6?|!~{N*RzNcZ-Gn`uiz@UbdM~ zgJ^nHl64uuY^81!!tXN8h>`?WX_LUNqkSg0khe5!IxA;Ufq<3jiolxLVjnJ_(IM17 zrklHNRyo05O|Bi&zVeXhNxPX-Z}bJePY~Deg+T!fCnB**)^^@vsUZzncSiDD6!}U# z-BTl}X?XNRu$Vjac4I>^wmOAThm<`*lc3kEK{RyO<@|gLM*qGVs&IWxj@lFq@mUA| z>h?qkC*gY2B(tMJJObe6oO8_dy0-Kq?)zKmz9*~(o*pH%gs47}wCN+_=rxlFXFk>k zU+TAJX3ZElFVD7jvh0E)>%Rz!S9P>-VWJq}t_M=!c`ey-2%buLPDqCN_;b;ydhIppv-5Sag z*83rZWWm#a>}Mg}s(4voNNNuovf)kRfe6*u$rrI_u7)cebrH^`F%RKK{*L7TbgMv$ zNn)A!idViG7=U_9XZ6kruxn~!Z-Fq`6|O7u(ZmM#0402f?uY$%05@n z%ju@A94`*%HXecN#i<+!3fLTTa^EnmVp3uJu^@wY2Hhe-w$wq3X6MxKX>3dQ&aSV;0^=g$4htogkpnYz}ot>vw z!0ym6(Tky^I4nhcC4A$}a(^pgqg(J%i=de`YV{Wq;}6g`fiLRBrqg0S9l{;&cnj?5GT;9^u2WJ z(v`G?!y|$L>pYrN53y0{#ExmgR2SMkQ-_jRGyhx-PMiH(<-I4zpoa@dYS@z;>t&>J z_&qdb+V@6x?4y|A)mPK;IN1Hf}-xV8k9$4kGZUHpRGtv!4 zo&&RsDs)dLay_BTq+3ZH_7_{!aB6)1(O4-ReQ-e{NK@$_4 z(|LE}hK3l#Jf0~P&srmH75#iAA+qTOjA_J4+(clNO4a^I`vMSf)*ToIk6dosg@ z6IdcF61teQSHgFpVjw5ucSv|i3ZYfWDLH_if1d|Twt1LFs(W^&HaX-;DN@|um`Z+g zifp&u_t$n!vmc7rUR3@#`uDI1Y9oC4nI^`qQ673!Ym5K6+JMvq4|n~w-N*Px^{omV zmjfoQr#&n>+OS6o2d6iB_}EcuHso0(QB z+F|1&EAVVMx?z5p?>?g<@?2`GWU5@qE{1>GWTmn~2`C!>^b5HNKe=EqW?HNld3nJS zax0|GOr_1~WWOju1DwGU&1wn&9~s*Ybftt|6F$hTqzJ-da>;5WEs8qE@vW;?7rV-( z^l^^d-p)OA&8-gnx}_=s+kq6B-INhS#_ocj>l3}L%l~5mB<0HxmMZnbM5A{0EH@nrPxk4s}93B!^dKJSz> z(0N_PZr$XdG@ZaY=^f^OXT6yDNG0l=Uyt2MUDDK)urkv|vPS2#T1rphr_{x7E zjB$cfZYL_APodvp!9Ra-=jo0ofv&ev>Q(A0*|IekR`7 z{`;dzvLgKXN+5UTe300%zVcW4An{2`5A4ReLf^Uz>`7uM1At5nSneoIJy2 zPz?d0$6}zF)q+$b!jXgENpYB zhOU{6R4$6bh)$F#d!qVWmMNy%4u={ie#`D;B=;@;s@}$-MmTBOD;3orA;}rsI1NLNM@U z)~xb@t_A;u06eVza}72M+|~=^Q$#K_$`4<*$9jt`jRwn`vMhPG7153w53Q%S>ca<^LMqZgFI$7 zShfc|5qaqKxK2W*pet;vYa z;n#`AUkrg!@B7*j2*lB|u_=CVS1FLq21}F1foQhrh;H!=?3Cq7VLn^^R}088fdd=a z@-DP-?ly{B(TB^TrIM$u=JuJqH*r=ihDR(t7 z`O`bwE9CMMG4#_;({-G^1sG)lXm8@8@4ko(YF`6 z*-)a`m&e%}ARffnf}it&Pal(mIcIzqoNFvWDT^n-=PLqpERIOPNcgs~x5udKyi#IB z2NLtSfF9pEDG&*=T`5y{is;hFFP^!0(fJ#IEb$Z|N&1~uLKl@z(!B!<80LB^KZ$nc z5W-wv0%2Or6lQlKH+bI-=q%;{&GZJ0JQzOI(%}LwqJS~TPqU~MD^q{9j$$4DSDPx9 zGx!81`sZC@u3B5NDcc8ir2}FSmx1AP>@5;kso}U` z4{>y>TvUk?nC$1|6yN*;8I>^bR@EIfn4G(Ew(B*kUW#-k8=Q?4#@Np~-JLK7cVf^p z=_H|;xmAZ;P<$leM3H(X*FK41Q$+;#V|rA@(Zz#$FWzy8Qbm0zp3(dJH$kqxD1qA$ zFc20QZ&}HN7!P2ao}6UswYd#ieU&;09=?pBlqdprbTAaA`y^-}cm?jl_={B41DYu} zzr8WQ?VtnuCo0b|DQ^Iqg=Q+`7HG4ZP-am+uuAMj2D8IZis;*KMRMN*iLrEZzaXx1@*xpbuKXY#04&^#DuHXLy}nJu7onA`{OG zux1f|tO5TOCPm=+Wt6*;;_5Ywf#QsuiMfN02}$%eu>j5iD~_WQ!cSc$H3t_Jw2lZ> zt=7bgi&*LQ!xfy0PQ?MH2t=l595jnF4y#WQn(#jndKm989pR;);*gpj-w5|ngha#V z!okCZXWOmNo4z){!NEmDd*PR($>FXnWVq1STI6p~E3}mS7vXXwRT|w!Pu>oKGW|)u z?k++~{k`J#svepc5nL}xQ2DQLi@Ac!((oUWg1pX1@#lb7IfB%acjKc$JxuN+*z*>nM-Ddjy*5!Sa^MVwg_Q(w^M6k+ z$nPi86nq`cgjDWvy?{aR>tuZ))=gdno3)F${qnz8r(_ZWi`BgRlP8C~vkh)vMBew$ zACxu{hHFg(HBc_2o=*B9~GZba-N#V$#PeN#x8$nt8jd zCgH3^R7Q?8qcBC(sl2(l5qKz_OoD{*V{kf9$W{%3PCu@4)t&^=uX@(sjD2j!*fv?m z5iJy#7E@kPutxxAWM8N^Kq%TjrwC|QX<}9Qq_+>77gN+QSLy$X$CNGdKqN+{yp{$X ze*dZtwe5s&>8KLB$a4@=s&U*Da6QfA-AK^*No4?KBkaT}2*lh4qY;TiN!m4tY>MY- z*uP;++p53-8@$KAi9=U3jcZ(|zeZKOI$Yd*gHr}Us=qjxPh#_|b-%I(H)h;IhWeA3 zSm=&6qr%%_a~IJ@%J~>|(Tn*xJ}cL!k}*l{f_yK3-_!8l}P?`aY>?J&|UHnR$Lrvc{qP-YN9`t(73|^YcO&CC(7|?pk)bGzyuayd} zj2V^Z&V>*q91j;#D)c*8Bmk>WpY{K-_Lgx`b>aH(kPe*!(jn4FcPb4^D%~X@DJ?N{ zNJ}aR5=wWYFm$O%BOxK(ARWVhjryE(-tVtJKN-xPz1LprzOVa=lx@JV>@zb0KInXP zrh;iK>a=z2ro1k%exR)jQXEI@hR@&UY@|g;Yo-{E*dmDQbw;0!ONL{BjBYLW{RP=3 zHn4WxQe}2}9h8(JuN>uE^;2vAxc|3R0?CFJG8`{Az8KzjAm+Vqb^%7I))J89U*UMB z7oH!w%bJ^;H*M~gdic+hmztt~G82;HK(&O1>rar2Hr9Q2W<{j4%Xw@7$0vvPgUwR8 z^XY64=j#{==%nsSuy`vNfz;U%$Me%jL@sL6Q*;$+1T9ft1ogB5QOvBYvW$r{Gw)J{ z{*PQd`FVoWw-TY^vcja7#PhGWdRx`an|=juFEZD}*?(trD-&fUxFKBzDbP+2cqQ#>#gv>_RePxuvklD~8q?C%RE}nH~lexiE1Wa2xJ^5}=st zN|try)z)BF#bi$zf30HsRmQq&iOc56h)TjaMh}@@h>!21D4zKtCGFrdWPCO`fpo|OzNS4*PFHWZ_Orq8_PDLuAkz9?} zJeP7cn~l>afpZg+g5#|fg~G6caMO;Ud>+kSNO*gg;?Vkf*`+(jK+v~W9$hpP*qbz& zgynL2M;Lb!1}-SLiMvpo$q=V4m=y!bW>WKegmT)km+&=f+RO)BU;{`Qw z7#qKY2S;;+R^&1(92KvBiRXX4;g9^0c=*7)S8G!Isqg@HCqw}yCa8vyYg0eb{%AYI zmXyT($0K)YzyUZrPVo}v?S>^&vssd4!Zm>*W-ZhX_dBnO# z5fdhinfzS)a-4?o1)Aj3Y{_2Q&PjJ?R9Ux@#5!5Q^hN>5$_p127lgA@65zDW? zXj`i%zP7Dwy}>QMrtRmQJ)+>iC%BL$+girR{$UGTNO0=sPoD&CqY)AjF<3$@XeIK` z$L+Pg&b!f4>0UJTKUG#dsMOx>6E|KMp9}=oN(EhZB<{1Cl1i@C41&;b>}V@ zvHFIO25-L zM_C+C9i1$k?Fqt_YaMtZ@WJ8=jUsgB&jbNY05CZ_jXwL0>a_$x9%XUxduUb3XLhgbsZ5pl&IL7zA+FiC&#%VoEm`&OcE_X=Q9Z810ifqsd>Y+=H+Z)KLR5FJ`Ef;o4&V9AI zZs(X~5O+~MZ@JA#{`>_=z(bNj+dNuEOjvg-lzvs|wo_fjJ6C6|w_wl(2eQm_U51iK%%eOGSNR!-w! zP@CKvoMs_QDG#Kr_`|0F5Ftrpqa}-Dz#bZLO1PCRzYd%FofkO51dE#}O}`~@p7`dG z3-0`U+PZcFsDxwqg|9hY8@YVseX!VD3J);H-EO%~e;iH!i&$#9pE0#nnyXd(o>`_$ z^z4z-&rfVHD}YH$G+$TvhknA2j?2IO3qAs%L!|cNG8rA`wBdWbz4W z&AD5yjd1Qh)7omdRee|yMlh=2?@q=o?Lu`?_}s~0tn{v*kwBn*3A#; zdlVCF+D8?~NDPHO4IayReeY2&(3cYq$hEY5e6(vjP8VSP7m8gb@xFT)c^@bREac*x zBOQK9BsaMUy}#-mLZ6M^FZdavkYrl6+3%b=^auLm2@+s3JHUnK$c8}6rGYf(17OJ_5IHDuhE+DMdY%Y zemV?Xvm>KC{DJ-y@mfBG;LpFGx&?rI7?J~|NarsH*t`KNk25nfp&Qbg%%{MlykdI| z@9$d-p5?pVo{_z>X|45hoa)zcCU^Q8l-aH3b?Op0Fa(;{)P>(=}1?Aij)xYTl0{Z-oNJcrLN^g>PK2@;qhh$R6FC$joO%F~0Nmh|ev{R2SsW>2z}7|!U?dxDH;7T1 zyr&>ONOs=;l*ankG$D3W|0L+-H368m0rKQ|LE`T%IgTtJt%!kd_wI98xa4!jYFZ?kF1 zd6_(W#@e31s$YF~4$uc1vUHQ~!R|l|qvg3}*1rP_EhH!*0Ma%eDBJ11O2@aKspmo> zsT1YK9I+4ynSX_f**yBLNu$ix{DwR_hH0c%irexcxgtT{=fWyejN&RYsdk!3Fkla= zW!?t{KuYA#^Udy?wn}_`Bgi2eku)9OsK&b8_Z`T_9~~Yar>IJAxdtJCj$7V|00-wa(>g|I3*$5_Ao}Q7#cb9ONaB#7%m8pX^y_ZT zKP!51AEevfYUH;bkOLglkKlP4?j%(xADhyokMV*LLYmsgcuUs`@aHXJU+pUZEZVpM zR=yTMyjI2toR)!VUh6<{&t~Wymo2n(Lfa|wFSp?N1A2qAbi}#M%V_^Yo%%f!9%&!) z)BH=KpJVUs`PoHxRNoqU9dr|z+;^P9bB@bOkxBd9&-aIIF}a}8E!k^93m6qje~nt(f%(Gmhf(y z4)+Z}8G?_3BX7E&<@~iRDI3j9zayrSbVZa7VT}v#?Pm4ewaGg`^;A zrZAwB5i5}d$CrC;9|7rJX^VhP!Uppn+%7|fPIkJkY);NS@J%dg+cCpowcv8`-H}tv zZt~6!_jDuHE;ES|Gpa_;`=bv!gc5I!x}A{jzn^V@djIfHV>T-*3+o3t!-tc>kgmYK zrAUYg)CTo%ZAcxdT{T)jpasnbv;|xj0=q=5AE1+$GXDUD{yQ?43pjv{YNug!H0-xla zK?$Sn!AEyAyJrT-djf!;(u5b@4Uj9VN`9dj3SLI#*2B98bT)#(e=ev!s9DdOdjZAr z@y%D!4^DH0COHuYo9KJHP#o^0vW%&B7U=a+a`@NNvL6AgRbY7mH8r2TOKR|0Ns-}WF})0O%A z%o{<6d!NNB%aH7kF&ZfZ4{45j=P#zGjkhf3GzbYUfHp_TD<8nd%JLMwxFBJUeTU9( zKWLHd&|gh0hI+WSu+^Y)2V~s3<(shGh)}3V8hMX#qAzLiTEr6tKyA3pt_m+7%qK?E zVoQ_0;M6!WciEY1(R~+klkgdu5tp1ks6PLh#!1h_#3Y{&En_tBMVzL8^|(W{4P-wE zY>=u3QTzQ`!TOLh9al_x`Ov@0uoK0FYNcZzvm-s@-fDqz+4dpH^*nRSz7_?%IZxy? zZh;wvHCBULo1KDJy&NDk`s{Bgo|pm|I`wsLpU@}X1u{!OD1EmWDKYZmK~T6%M95Ob z{byACD82D?lk=RD?~g3o11isf>7|vJ1_dUV3XEA}G&8fm>O96^`b|+6Ax>;$e1wsQ zh~(ngri+u6fK0>azd?90D30l*{6xz1s)B)o1Sj5|?FhcO6J9ZA)<$V&qy=Tk;kZ%2 z#XD}7h+sqF81hKbYvDDDAL)2EHB%6nAQXqR^=3cofx9p#y`Kma4zC?P|Gj1h6uo@5 z_9}ZQsVp{>xFM=q)7;rZ$copS21|lyey&8!{X=8nGoMS76_sJj?3knBKG$XArjf14 z(7g<-ba}Xjo#5784%$$Zr;U5S^VhJ)n2~}Eab<5cGZx^Eu`BB2=%?Md(h+K;`jzF;9bu>P9m~oU_me$^ZJ-+m{(9uwL^PjKY z=~;ID(*g}5rQ>ogi_MB|vyNkfw(r8oM-xO;iLK43adi*`Z@yf9t*$jfn*C>n;Jt&D zGqN|U@Q6)pz13>*4gW29(yyx5J&)=44p&TGO5p;fDhCQ~u|>{pF$z>z6}{3alD8L< z|Ggz?x?oXBIq-1?j>`r;qBk+27-IlP$UayFG<7&Pb@=>u}9f=wBaqc4>JZwvD+IV>MV%_95+|qJSZ*qAZ^~PMJ_c#!MJJE>y{KxyxiTBWorP$BOLKXw(m`d z5Yoy<@k-Mjdjz-%Z;S*78Bw4pVPShYv*v>1ZGtA`4mvd&(*oe1Yc++xBAD)h3MU&~ zV<65bTziGK)zxw`JF_tB&=&0RJ5fX6jrdw2n>c<5-nXzHI27nA38J?!d%)g+C5V%N z0zJUtAv=oK?gbhTzj0UP&u@?~dyp=k6(C}h5)vkiF>_7ZCZfOyFlF2JdsWo-87wbW zKu;Y-W8o%FDevyTSIarkbTVnmOG{R8;KbKW<_;^L%%%yV~XF>14i{~15~SI z;N7TXLpi%-(#A{HZ=5wFB#Y)&XaQ$H$u|Z=e4}SMh71vse(O)Za~Wms)AnGXS`jTr zd*mPKILWde)%hOVXWWL$c233VW3ZOkcNrUVG*!5u1#5x?jXN1nh^vSry<}$kV)E`c z-tx|aKhtDdor-?}Ji^o#7Z=@Z(0Xp|l87RR*59V558ki(T{V}9c>4nP@Fqku=6;u- zo*{KtvQj2i*rU z462~KpSf_SUG3ddbSAm|>#q6yr|6<^-ivTeTB}ASyr3a&tw~usZ_}8}ubZQ9#qelk zG3?u@Mft}GuBN0&Qea$IW&-2&Wha&sVn`W<+LE)DYA}(;mk$og3oLO#4=NR+RUA}6 zL|cJrgW$}Vy`jr-&|$V-w36i|m@TuSS@Vj#T4NUx5gBebDvWrBHeaTp7mF3weu^Ei zhKpt&OOrk8$i7O?9M}-aB9pSzwaZW?b>~;wR8^Ft@W-1EaY3z7^H>fyT>1ct%^{&c z*=KktC@2SeXbvYpyXwlRwS04{)l7pudEpi%!hRL1aP###QPu^@so>vEz)7m48tu5= zVA&YQ#tk(oo2y@}D?O)`4%Vw3v$f7%dIENgwQ#9ld$NH#9dEKfa!q&qMhQy#`0j4ZpE3Q?$E2cnFjfnJhk-Rk78}96o~!&5;M(#DtE-o8@?CVCv|u2LSwu3L^9YG z&EnNO+KFDODOxeEi?=!YsFjY z>{lY{%Ie+`cbgugC?qjkL}!aSf|o+OcWMBw#{z9(r(r%duRlGtJy_aJkCkG0guF8m z2R#)&|Am@Y_gvAi55*xpo2;FKM+EoM>Ra7cQB0MH)a1J!>h(+SryIICNX5BMsQ9@J zNXyu$&0^Z*DYy+tC?EqAl~3m1S29uEN)ChuXz&!LCAAeLXF8Nz3XF4;YIB_!zVh#} zMbS0?(+0l=h=L4wLAk0glH)V#Tyxl~%&V}j))|cLEJiI;Pf(c5(jRxY#{^Jvsm`dD z{f47Ohd1%_6pdHpm*+`GmouN@&yplRJJMI>hfq|7B!zbH^o-;&AL|#aNFjx66Ym)7 zRwUca+UJ=_T@70d-`lU?YmSoy+Y+}{_GjUzvAnw3mke%oC|YNFXD2apW;nKV8d(?E zIc95DqI)bEJWu%!rSi7*;I#Zxg+@gd!9!>LLBG@qpaz1_%fFUX2Fz^hanWn3MFz#A zW)-G_w)@OmXKo=KL3EA|4OvM68DiuEm;unbG20{Err~Er^JkpCv(2C*pc`Ma1z&+nkYv<~&?HDE*}1=%g)qM5FoHu`4i4jX$!3OEQ_^Sky% z^9!*-=w-)u+Tzc4>%->x(kt~IL$a|(%&K`%cH~$NZM*-XRsD=Gy)uBVt$z!`u>oUv zk3FA%lw=ciq}LdP{OWBywJ2j=7jK7!T0#dDrQ<)8*cTn?(*8!oHFq}{nK;ND(cr{X z@pos0J17I8GUWCT4bMA8<4H-NAE-h}Af^A=i~QJ1K9qdbo@Z^uq=OmaUgf@n|K$Gr zZ%BYfMkYm&^+sttiU`r=N6L!^|A(s;J%somj#H4d&E$J(=?P-wYbF0W)q6-BfWD5r z_z2wU*9(?!D073Sw0;@*;URZvfM?RfsL=SqNjm94!(UlA`X5xMCX)?j*$}=^dd0S) zC5#L5C6B*HS|uJzmt~#v{5z!QVhdyzrguKMxPW9dJy(KcwAUsdPKUD!jMmEQ!}b{( z?8l=6Zmzu^k#&rG`PfPG2aQRmhcQNi24$S>M#J164h^^K1dsIEX_Ju14&01US;Ml{vn zoL$-;BBJikrdOHc_lhVmpETfTod`LR>Iq`}O>xx0Z~JY1=>ga{mB7#>7j$YH!b~{y z_om>^T3O&Ngh;vTbV<4_KtN~3F&tMzV4wTU`xbJ50=GatOYDy4Qdp#S>AwP*lZC*@ zoUEarg%5dB%0R+&s08hk3#k@CYr6qdf?DT$Q&oAP;V95A3J+L9)1XIPFc{rg9rhe3>j!W9F!6F_9Yras|<%#C$8F)2julC6f(Qol-sCE#pu_ z3veOapxJN+h8ESB+KG(v>i{Y#0B#ps7jMP_aHn%c=#vj+LokOyp7;vqQlVyE;e#;j zW}(NKBC-IdO5KY+KeZi5R({%~c-xbif;7CaG=YbazS!7urw{F{O@P4+6bFaLfWuw@ z^d^J>k^dS*dJpIo(3Bqb(&afM0JgizO13eVc7(KP7?uPvx7YCEA{*x7J5he?FInRH zJEevVaZ(rSOqBdKG43U@k3KLw3{QxUCongg7z<$CW?h&DDe=K0@j5#|iCO@{)n3QZ zNps(NL8R+`gYeICtpsuD>;Sw0gR!^v6QyjIZv`EmZh}YhC|kHN0%Tn4p~fc(Z{i;4{F?l|DIQw10{e0W&%;!+s&V0iVFZ4P^ZPL-liz>}mGp|TX}-$x7Vzzh2g`~%!EN3L78|zlL|3$uzI`+Nb*8-V(Ns_jGAe8=4iY>x2c5MW zUF*wAd|8GmFojaF4g}`x$cl-$94`cdEsf5P4qCIj+lq6nUr^0@APQXoJzV)SZ>$6< zymRto^k;}p0uNx}rUO)xc^)W)q``FT{}!?mEy~$ghc0=F-$6)8K3G?LXWzQJ2? z914O=_0S(^k!@{lI#(-i$#32}enE#nf}-knfemm_zMGa|m1HMJW*eYvh2TxN|1ID= z*q1p{S4JFUDR?)9O#Xe-Y5C4_IyO6#(p|`My88A#k?)<%5;~7cc#+3($SWpZty2)b z*>?CV)0^W{1VEY`NS+BCt2zj0l%61G;aNbEJ_;)5f(|@u%Bl=P?MqkO=v&QAP3BJq zfHzJ772!{~E}bg4m%qGP&*<7eN8-WgxgqCQDTncn6#Mcvd?|MhK_}z}(UX?K>dkqI zWz?@fl8P(^n%q`c0UC7grD?7l*6ZGVl)o7Aei9tenVxM=Rq>5iPe;GunD_e~`|akb z;~Mf_n$V$w%`N#y;DX91K(m_QgQ)X(F8xn`1QE|vDDBL*Nz10vu^wssWN^8WjX@zj z`_qNFJe;kz#sSFnZxkU}7*{a?f>o(R{OqC8xsw{m|zk{_jW!n+iYjz)O4Yx#kw2A3ixGBPru^+sbK z*ky3V%l0rdp3|p76tI@i%tCXxv4TSw(%ZY9qff&xfDPhC;LT-rroPO){zy8kZ}5!; zGEY4-@+6KVHTyY-fPg@^qd^`@Vz1C+iN1+v?6Opi^PtV?4GRsuJLRC)hFNsUL``&CfoDszj z0&C99Gv)7GNVrDg+dC&UE86In`wKym9{pCPa-3mSTuZEKv<}X<7VZqcf$?8(DfpKC z^w50Sx9MnK)()s^JyD51=aLM`2$94rkev6kc}B@zx)u1)GvUHS7PTnW78VF~f!J1r zx4A-`j)O7lS!Ut=+)Hj^*LSWRbu8_?j6`Zi$iN}qXNvhJF(inaX3PCc7TjcvPm`R z9r>LQj2*-mDbetO{`j4SNPZiY+n&9=*-q!4y6U+X?-DE#E9N851)=jUJ z+fBDq*0xTdz(OP!i|Yxj7LtIvSP|39zG!HlL=hsc7t$n%alW2tpA?qQp<$z`3YjPG zSq2Z=sxc5)8S?!a0*iL}1v92UbXTGW-Irz>?&YAxJ@2ish#jhv+d=4Ez_v~cHBwlp zN^Q+EZAx^nZ(nLt^IO>e9ozwt4!ojgZGXrQ64OLuW3`5%MwWl+-=Ad{ICMZz+0!9L-cp|SHDIj~iOG+iyy zW{vepYx_z<^_Ef5+n$U;z7O5lN8(6XGU_SP+BU|||Ctz?*rI~!CkNz>Bw)7#dP#Ps zpnq1z{%7!)W-n1AMcK7K#O&Pf(ye@nN(M4ZDbP_LAS2M$@j8@XIE2{p{~-esDCmkH z|7nH~LSZO}N$^JUrO>Z`4E=iv_J0^|Ex<}*6bwokuXW|cTF_|}!ZW>o`Y770%Eama zI~^ifAPsb+i$Nd&ZD3E+^D@e|^XBhhQwC(6fl?U-nppD&6Q7#ZzaN<}hNA^eyi_?8 zJ^go-VT_Xff6}M1s%S;yQ%w)f8^Z7{cG2{+=tP|gQjmOW_+flJ>0^s#_Zy?EYWn_o(X>zL|VUXx{6 zQ|KuL#^sKt4g)*a$;{7735zah#ynkiOJT8c?%_QjZ|}Zd8Py@7!1(-KQ*Ok)x+6bW z99Ns-!V`EHT3iy9i2w&egNlOY$0t_8%-H- z#TYg%(3ml70|#%3HsU~n)2LCS8Ii*$T*(Lv?Se#wi_+zdVbo z+pO_bsGat@zFhFQmq3f9M~#+1YrZe}XYO9A45|456woO(?9Up2*c>_1bTA2?0afxJ zceRA)sLvsHsqSR|d4xRikm`J3@1aZ^^M)0GOuTz?Qqh0UuXGq5R{9QVs1&R8!17V= z=;&(@_JTavOjuC+j&G0#mWG}$ULYyl`(BG%N09x`f6D&tbeU>g82Dv%%OU{G1qGaD zMZ^{0#(&8B>FQ?So*vBy-e_h^`CBDeTG%q8Wf1$_0_LSApI4A)_Me4ir$Vof2`ee# zNMY4v)Tj9dPKD~YwU9qU)DsoRkC#ASASp@}Eye<9$Q7J^&!*;qMa3512L_Vl0K&?} zsz3Hp$e$)-Pok8H0=q(-ng89+J~^T1><8|*-)g3$-`}aw%)yMgSi!FPtsxD}( z(Tzw60G`JKOk9ctGwHVfdGjaGryxpO>aYp`ox_YtKo~Xu@%r#X-MLy2gcZvEO9dI4 z{B&h2UK&v-egZm*NuBmmbm0!?KeqzN`iRZY=;s1z{Vc&H3paUqFIT==@+nD$dc}RyzBHxwDkB?uN{`ubZj4p^PALMkI0ZxvLpr5WZNAr34 zlBMx3$jDhmYyV?8vkr_}4ke93dHx47P=nKARic>1XJem6f$MooBqRSKG8XDS)P|Q+i9ZQi#cqCDwZUrp^ z7i(}30|qu}4(P)MJrMt#^uTuoj}%TTeS4P#?4nPZ_OH%g7_YuN_|>ZoX5{h;3L}7Q z4rMxSyFPh~d*ldar^cUspO^qS$tR=-+~M!vs?~HLd-DU(Srb>^u*$yu_N~9Ts_^mu zzClF~V0bW{_AlgM8Pan>>0he@<#iDliMRtE+7FDgoh^=&>&w%Rpl=xkuF@+3jVCha zLjYgUH>MSHXN{$meB|+1M9~Y>>X3TKX&?qml-twptjBj^V;v==yi8A zt3?3B-Tm^!(H#RZP!mx_vtMYEmGv;l|E%&;MeWJDJ_ zN@ZU>K{wsn22Oom0JM4qjF&Y7(;iLSH*S{u?0cLGM#$^aibM}tZ^G>)_-(!@A+2wj z5PiQ*L1#WEiX=@SPfay?NfPizO9nRkId6ahz6m^5wzxzE4}j%wT~fr z2T(bAbGH1nu1heD98C0%Sz^i&tTH-Lo)bU(!SbQBGg2nb4~lBNNRZ9&5cBZxFbQYY zzuI<)Tk+D|NW$+-qZtFa0iED5ZR@#6m<+c8yZi4*C-b*@Hy{}ye=By@Y%!3Rmp3%J z1-aQi*%W+^T|=g{7zHu=SNtIF4g<r`16c~{(zA(Xq+pBUmNlzbY})@F5!hC>WE|dm`0~FMS;2(TA_Q4jtu$*# z6%CFY9s)|2?*KZ=NBu1ppyxL);*+vgi{|lPJd)E*D>OPPK6CgU7auDhgvRphTw_-u$w<4 zk9JL8)vPv;bwc{)iYdg+^#;@|?+K*N;LFS|1TJ?BK#Yn&3}KBK#Vd$)8yM_X8{Cf} z%Rk~YZm#XGKZpF!pnHuhDpx_#Sh|P+L^J8tv)N zD2b~;#QDeP$3O_jYX^ZV$?Wd~C@tnBAEfm0S*imj&$4z_%S4(zRZLjf4)%#k|RNA$a8o+KemUeBu_`1VZQPKTo}y;Ej$o`v&l^{-YQ` zZf3Iblue~%oul}wppXGh%XAQ+n8@z0MJ`LYgTy^!NTU_l0)6VFi*zYUmC~YOL>mbt zMgcaG_)(^)S4pV8-LU|uq>>Fjz22-5JQaLCm5>lklD!6^y{s<20z)vCMs#n+Thgdr zL?Wg4c`A8c1(9ZP8P?~6YQK*HG{mGT%gUaH-sPprQZp|=CV$dl$= z_#|z)gSS&6F$wA~lyDHG;gGKtlNb7NGTrFC$+K-dXJ&dlMi$yZkUiu zQ0n9|)g^SU9t~5P_Tyz2L9Vq)_$B_e*IDIi#P09+brd zkqW%d)joKL%S@Yh=JFQSQ9_C%fSZ;9lK+UjJEZ4k5f-FW(d?wMtm|itiVQ3icx(qX z;Y-pW_#@+vLR&UAYPl|mM=hf69HR1%*t$jPbd>)>eE9?^GYrljFYk~vJJdtFOMnR@yDaOL z`t`d9mk-G}29_;Oc#xPk7QCdaL_MU6F`qmE=GwlfK>@t&-!g^?%7)zkX;w=VUCZ8V zwM}8@46vqbwIXP?oh3ix=EYRnSU^l|$J!^bviWSqk{$EIbhh6^I|8XT*vzfkzyAr! z{cY1`AqQnTi2McP&AktG>wYF?WL{Zb9&71Bn=IzH+ou|@CKQPRONuu$Xj~U~_FV>1 z$b^Qg^^B5)^>t$}h=x#VprU?7_oOx?aigO3BKm`RHus#{(NHg|V27usi?(sFn zUCLz-gBdoiB&)KEZgqDbO|+Bw>aD;uP&tZm0vIwEX%<;!(d|ByR1@5$C?SmW#9VlZ zt7;Gsuh)S+AE-w-lzY3KKi1aJKmJ+4{5ngkq88tJ64`0=-g7^%Uw4LYq|{%v?hAr$ zhUYQdDUfS%^92B$b$#)|NV}tOF_^kQFQZ<0M{0(g+o;laTEnkNc2Oo;o^>!4+pQds z8I3M}KFn~OD9ylidQ>WTwT=J`JtWDlD1k~GzliYt_dGs!Y=-tVPyvD62BZkR2M6K+ zMsJ2zqe~!?=6k&TL%ZG_U%j!F{?#C;A(#6klf`v%x z$$5L6Y6~6=W|NqZfU)Hj``P>6LTIvE{sTEe)>hV(=U}*i zoSwX>D{d^kg3U!BJJCVnap%`X09UCgpqV~ z+g<0wBNZ2QmqjJ}P4aR%a40bc;-G;H?;Z9zB|-vNj4Z{Ir>MX70e3DfS~}X~ePar5 zkVQPFwu_)_PEqkw91P+Oo3H)a&ZWjbZQej)v?Y4=7DmF>$a~oSarTkTlNfRbhK)oe znp%C(#AdzfUF^lj7^OvT9dvP!_0W!4D52MnxtF&`5fq_CkcYv2o5El5=JjDk_+d0$ zhgoBR<BXbU|KK;*NJieSeK(T42n08 zUBPt%Bj1JBaSJ*RgAQYc?frw^V$5nU+{BSao$k%2VRizEZ9$t@Of-PDwrpg10HY@! zh>`~vo?O>IZ;I2$JX>X&qF@tA?4{+_y$CLD(CB;Ohq;8Q5Kmi%!x*lazbU$&ZMz|s z>_NeVhKZX1+-Ydojibl5Bbs)KOBL)f;j_QbME)z2z|r?C1$o_&=*&e5qTWKjSzR6y zNJHu}yHEX@%R+UI-`22Huq5xMGwXg{_|mK7*Dm_9m3+yM8>`HYXcZvnr(w8kgl8oW zGZdWM1q~sS^g5DZ(g6l%I>C~-5Gbj$KN=*c7MlEIK$>q0{<-I=)N(b+l61^0(**24 z8?Zk;0o)JO(-E@Vml&c^ie)El4x!H(F+pcK1@oJapjlNuC2GuuGea?oD`n2*>Dd%o z7pqZEY4StWI3ak5>-0s~TN1uq_>LkL-)}QIzwuF(F(uLRco7)D8wbnj`3aW@bedJs z$(R-M~!UdmDLs0#Brz@`Z7hlHf#1$_&+MmQ%{G%HD ztu4%GQ65o1R_9bU>B97}9Iy|^n4a9Tkp$`H^sy-U#kIF&w^6Ej+jp?|A>LuG!Ss$$ z8;69DdD>XPmkYiApeOhj43G}{Ob8lh6R9ZjYnkT9==By=EV8Ajs5QJ>L03D4@RI7= z?cEyy#7M-b-!2-OCI3IqnFqkwpA76z$1S~aND|O>syMs&>CuSM%O3%PT(x~T-4C9d zKSYK9KF*8;+3lh`Ts$XDaZL=+Hce=P+Nu46o&Q}>tt<==;^_C-b6EB4lxbcX`i>{eC&Gcu2X>*---AyFefM;IJ&W~!E-~K{d}97Dg2|^1 zm~A(}N8v}Yk5k&eUhVHH2}hs8=nHpgvim>(4eT3c&Y+>KY)}uq?EKeb{=2~EO5g&i zI#xdX)gS!teIx?+@uw8gE&;M3{qK)Ki_ZYyj_1&{w|})CNWAaQ8aYt_Kj-D8T%+{= z`N!aY!odYWDjAHN-u&Oo1i%s%_&M7RB=Yh9%NKB|Z@|}Ykp|wo z`C!j5tot!Fr57C)#lXVC;^sMl-vNxkNiKGmmMoso)<5Uu=8ntF&FyP%4>B?~cJyH@ zX^0I75d@Q(sJOU$5-jTtd3l~axAsq=-;l@o{;ATYzm@@lHYz6OuCQ>!Scr8P$6^b? zmxCSGgb2!_DVf5i%OOI&+0_m!-!fs)LbL@ZA??m8h!zLd3)6;?esNtnN}u^*tTp|q4=G(P>?$9}_gFhc2LxAI zGP+!zdpcftYi{I^Zr3`_jvFZB=H zv-e@dpFz1Z52mGMR?^Qa2?u~f?w(i1-X$nV8`4aC`)%I&*vz<(K5`tAJjE?ex$s+m zA0aWSnWtcV7#qKfCp>G$q!7ciiZHb+NO``N+F&6oE2EGG1tJDq`m{2 zr*Ix95m@^A(~ zU?f=pL>WB+Y#{E~!j)9RZ~O$&;OP}g9`-SvfC<8n4xIbj@k z3-EBEpmu?yWQ0Xqyt~PDANZ2Rr~FxE?f?+F&L~Q`8fHg4=R78@C@#JcU> zZJf|Bj2{*wMQyFmF99;qIOq5KYirzsanF)r1hre_no0Jg&9Rs-uZLXB<7+Q_Q*^n7 zd7{E2C?xh9&izl{O?*~TxJj!nGF>{bP7ZqME{DKsn!F`Q<1wa*;5XU_rqdUYRy?68 z_(-7kqhs~?O>h&cwI7Rim(&n5bCe>VfDi~mOpvpdvTHPTAC@o*y86ZpA z+ry$M0HT7QKo>I@N`UutL}WGHGA3>k@Qi+mTiM|l$2?9lixPfzHCn}-(_pgx{hjwx zlu$4Ai?b)Gc1IetU>iWl9RvU4W0Ycd>X8-S^OZDp_n)8-Z&?V&dH(r$q6QFtrf=!G zIeXg_(^pC^MKWyzMp(UNyyQepy8p=!|MM##$dU|BUu$E{s^$q!pwk;ZUa$FG8urJ{oh zZnIBx5K~~y?!`A@!sY|Pe2OINB6h31rM!MvfognKRZ4Wl#Gvl0&`jy7$%FnXKIbd2 z)at+&pR1V$i_^Y5F}uH<%I;vz$dd1cU$*~A5smNqdqn!%tPk}Q1`gxzz*m*C1x{B8 zG2v?Bve!mNOR+C^n?9~m3JB_HhtCa4yq*En%d^8%fatLF8h@;<|0cFI)=1t52yzZi zIytKfeZ*N8c<-|Lzlkut@NSs%ENS}YwyLM^RYvA>XnIgJe4k2m^4T?7ykd=Y)RDhDG<2jaOb%xEGe$XhU zR-N23-Q1j6F1H<0IXC3`<+o|_V=N*T{f;NU*Yo1*%Nl-)iS}S>7cP}| zo=>f1SI-rHi+jZK`=LsGO*fe#Mc$lUkYLBrJ0+9PwAzYY>u2{zrHKfGIY_AWG9ergq6r0bE}sStvCoYrglV$9>29&# zqaNd&qNPz#T4}3kI$ff!ylJcDnY?Or{q-fb^{kV3C{R;uV_vh~(Fpb6R2|3qi+;i9 zx7gNPtLGoT75EWfE-9dX3JPk(vuJrBN%Q*gg;%nxc-Mpc-zQU+@#kVeUi3-26}@9b zykgcBiX$;~T&S;;WSZu`b1Bu^UquY=Z|(d9#S1w7r0`+n;&{2=%?c7jZ^~cDSAl(a zZa=Muj<{`s#{gCrDlpQqj{2ap9{WLn-IXv&aIyIS&b`PN>bz_AD_-W2DHZRSSVYN0 zde%da2Cp)7>wR3sb=jE;`JgLgSQC;1C!jXj`-P4#woadCjx@odBTb*kGfySky@kgg zTw8K)-@RN~2tb7yFF1~VfZgTQND-y;n%IhX=`C9CoFU&Yztx)*^o$aQIo1M}MqSFo z16&oqb?1keq%4C=Jr3uNfwaR?bSom7^PkWVsIP9m#Q6)JGmHc^9r1jbZ^-jE| z?O3P){i|VVakaCQP%jOPTlUHjdO?)>!De^B7uF$kPUGaikk1R~rSTFVf3;-0e-608 z{5xE|NanuG(sg`o7K{Eo-uM-1wY^*ToR4Bf-0)L_EMxE|b&Q|!@YeJO z<8ZGVvx3u`yG^|mqx4R!dGYkZCYkr}`)&{DX%x?ZhG~7LRO|i|q5=9?6>^EqShRNe z;45MF*T=|({AAG|RJYS9x+Zq+R(E{xMk1NkMm*31wgpZe{>uKXSL=GSh+yE|@b1RG zV)DeXgd98tbEI!w)hpK0$kU_$k~WJRg(zJFl@jGyM{G0!4r`z+1>X96TOi+_G)<`8 z$E=Cf%opEsmK1R-=`lE@@4sHiqfPDzfEHu-mRw=o5C3(J^WX~OzP$;tOIKO`&j#vL z^D-DkR)XiabjCLlRS!Z1gJ-k9b9>A1q}RcFF04O@`uTg|nn^n)mi3Rlt3W#8obPNn zy@MA~vM?F5S+x5t+tqg)1(tLF; zPsUBEaIbo@^@_v98LYm1pI-0v-m7jU355~|BEK_2O0oW6Fiyc?_T z$MCSAJ0b9T$vlR7;f5|4O}{i+U#d|Ex`nTi@qsIB(df2x$Y6h8A`|})7I^%fWI~Z8 z%PbT3SOxL}a&xB(imWf>I{$Aa?MJTw>#Q|C?-G3v&#l~_`Db(7QQ*Zpvo9+8F@#&x z+^J-HA#?rwkHd94PL%V~3kS&e^1d#I%_uVacFI)usE<$8xW6F8t`fH~Rh$ z_5FtD_s>?k-Y*kytocHxv~5+%vUBHFXCK;{`uE6=qdE>cCr-@L7QXeP-P;G06T^<3 zn$Gd(I?ENtj3V|A&!4mQ`)gba6q|b}scU`7j4rDGC*;B~@})~WAyz1FM$f5(b>s=}E?Aw42S$y~v|jD#b8D4!7N z1Qkk}#wD*}*ZiO2{YabX`G-zH#btd8dEFM`uZpC0?#N2Id(1)DO+Fc@Ec?^346!zy zw5`>zRxUT%6PaS*2F|B5-cRZSB--pKXvyM&g_ zoj{IxOCI)rxRh|zqp7))@!_xL?hC-G>dkW|{(uje>U9lMoi}cJl3?3j+2&g>0>omd~qN#UIrwD`p?V$Z>*CXt8j~>3^N6 zv)ru2Y|>VEaci%z1GYZ)6`!_}`0?lc#^Z9GTrxBNcBL~M{N&HfIX^q@7DwtlzMHk* z?>?92-od}n`T0&Rr`U$sF3Jn#5}rn>fYSc*i}{*|4cGiP@tS?ILr{6o?MPtaZ=YH% zyzJyr(aF2!bILXySpu~8n`e6V?A&dcvvk8G_&3L}E=~Jd1xh8t``qUP`*=qM1K1Us zR7I6~gPK&GAFJ0&KHuO4?0c?$8KC)h%>vhI|9m3wg-d?Bb^4G|HJ(mT{oF-;3OB39_ zOm;=;0|M22;#|Qz;k<_wayL~eKnmDE8ZQY1urDy# z9waG 1)) { + if ((has_biome(SNBOX_DEFAULT_BIOME)) && (biome_list.size() > 1)) { H_THROW("Cannot have both global and biome-specific data! " "Did you forget to rename the default ('global') biome?") } @@ -157,7 +157,7 @@ void SimpleNbox::prepareToRun() { // Zero the cumulative tracker of CH4 release from permafrost cumulative_pf_ch4 = 0.0; - + // A flag that lets run() know the very first time it's called has_been_run_before = false; @@ -477,7 +477,7 @@ void SimpleNbox::stashCValues(double t, const double c[]) { // simpleNbox's point of view). In order not to trigger a mass balance // issue, we track it and adjust in the mass balance check below cumulative_pf_ch4 += rh_fpa_ch4_flux.value(U_PGC); - + // Permafrost thaw and refreeze if (!in_spinup) { // We pass in the annual fluxes here, because want annual thaw and @@ -547,7 +547,7 @@ void SimpleNbox::stashCValues(double t, const double c[]) { } // Add in C that has exited the system via thawed permafrost CH4 release sum += cumulative_pf_ch4; - + const double diff = fabs(sum - masstot); H_LOG(logger, Logger::DEBUG) << "masstot = " << masstot << ", sum = " << sum << ", diff = " << diff << std::endl; @@ -744,7 +744,7 @@ SimpleNbox::compute_pf_thaw_refreeze(string biome, fluxpool rh_co2, fluxpool rh_ch4) const { H_ASSERT(!in_spinup, "We should not be here!"); - + double biome_c_thaw = permafrost_c.at(biome).value(U_PGC) * f_new_thaw.at(biome); double pf_refreeze_tp = 0.0; diff --git a/vignettes/articles/Hector-Wider-World.Rmd b/vignettes/articles/Hector-Wider-World.Rmd index e693c2a24..ff96a20d3 100644 --- a/vignettes/articles/Hector-Wider-World.Rmd +++ b/vignettes/articles/Hector-Wider-World.Rmd @@ -12,96 +12,84 @@ knitr::opts_chunk$set( ``` -This page is meant to serve as a reference for places where Hector has been featured. This list will be updated with time. - -## Hector Publications {#hector_documentation} - -There are three documentation manuscripts for Hector. - -### Hector v1 -> Hartin, C. A., Patel, P., Schwarber, A., Link, R. P., and Bond-Lamberty, B. P. (2015). A simple object-oriented and open-source model for scientific and policy analyses of the global climate system – Hector v1.0, Geosci. Model Dev., 8, 939–955. https://doi.org/10.5194/gmd-8-939-2015. - -The original manuscript for Hector v1.0 was published in 2015. Hector is an open-source, object-oriented simple global carbon cycle model. The model takes in CO$_2$ and non-CO$_2$ emissions (e.g., CH$_4$, N$_2$O), converting emissions to concentrations where needed, and calculates the global radiative forcing and the global mean temperature change - -### Hector v1.1 -> Hartin, C. A., Bond-Lamberty, B. P., Patel, P., and Mundra, A. (2016). Ocean acidification over the next three centuries using a simple global climate carbon-cycle model: projections and sensitivities, Biogeosciences, 13, 4329 – 4342. https://doi.org/10.5194/bg-13-4329-2016 - -Hector v1.1 contains an updated ocean temperature algorithm to better match the CMIP5 mean. - -### Hector v2.0 -> Vega-Westhoff, B., Sriver, R. L., Hartin, C. A., Wong, T. E., & Keller, K. (2019). Impacts of observational constraints related to sea level on estimates of climate sensitivity. Earth's Future, 7, 677– 690. https://doi.org/10.1029/2018EF001082 - -Hector v2.0 incorporates a one-dimensional diffusive heat and energy balance model called [DOECLIM](https://pure.mpg.de/rest/items/item_994422_1/component/file_994421/content) and a global sea level module called [BRICK](https://doi.org/10.5194/gmd-10-2741-2017). - -## Hector in Literature {#other_manuscripts} - -### Journal articles -Below are journal articles that use or reference Hector. +This page lists the key model documentation manuscripts, articles featuring Hector, and other Hector-related resources. +## Documentation Manuscripts {#hector_documentation} --- +### Hector V 3.2.0 -[Climate module disparities explain inconsistent estimates of the social cost of carbon in integrated assessment models]() +Dorheim, K., Gering, S., Gieseke, R., Hartin, C., Pressburger, L., Shiklomanov, A. N., Smith, S. J., Tebaldi, C., Woodard, D. L., and Bond-Lamberty, B.: Hector V3.2.0: functionality and performance of a reduced-complexity climate model, Geosci. Model Dev., 17, 4855–4869, https://doi.org/10.5194/gmd-17-4855-2024, 2024. +* Hector V3.2.0 uses the AR6 radiative forcing calculations, includes the permafrost feedback and carbon tracking feature. Hector's NPP accounts for LUC emissions. -> Wang, Tianpeng, Fei Teng, Xu Deng, and Jun Xie. 2022. “Climate Module Disparities Explain Inconsistent Estimates of the Social Cost of Carbon in Integrated Assessment Models.” One Earth 5 (7): 767–78. +### Hector V 2.0 +Vega-Westhoff, B., Sriver, R. L., Hartin, C. A., Wong, T. E., and Keller, K.: Impacts of observational constraints related to sea level on estimates of climate sensitivity, Earths Future, 7, 677–690, 2019. --- +* Hector V2.0 incorporates a one-dimensional diffusive heat and energy balance model called [DOECLIM](https://pure.mpg.de/rest/items/item_994422_1/component/file_994421/content) and a global sea level module called [BRICK](https://doi.org/10.5194/gmd-10-2741-2017). -[A permafrost implementation in the simple carbon-climate model Hector v.2.3pf](https://gmd.copernicus.org/articles/14/4751/2021/) +### Hector V 1.1 -> Woodard, D. L., Shiklomanov, A. N., Kravitz, B., Hartin, C., Bond-Lamberty, B. (2021). A permafrost implementation in the simple carbon-climate model Hector v.2.3pf. Geosci. Model Dev., 14, 4751–4767. ttps://doi.org/10.5194/gmd-14-4751-2021 +Hartin, C. A., Bond-Lamberty, B. P., Patel, P., and Mundra, A. (2016). Ocean acidification over the next three centuries using a simple global climate carbon-cycle model: projections and sensitivities, Biogeosciences, 13, 4329 – 4342. https://doi.org/10.5194/bg-13-4329-2016 -This study adds permafrost into Hector as a separate land carbon pool which decomposes into CH$_4$ and CO$_2$ when thawed. The thaw rate is controlled by region-specific air temperature increases. - --- - -[The role of climate sensitivity in upper-tail sea level rise projections](https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019GL085792) - -> Vega‐Westhoff, B., Sriver, R. L., Hartin, C., Wong, T. E., & Keller, K. (2020). The role of climate sensitivity in upper‐tail sea level rise projections. Geophysical Research Letters, 47, e2019GL085792. https://doi.org/10.1029/2019GL085792 - -This study uses Hector-BRICK to investigate the effects of Earth's equilibrium climate sensitivity, or long-term temperature response to a doubling of atmospheric CO$_2$, on sea-level rise with a focus on the high-impact upper tail. +* Hector v1.1 contains an updated ocean temperature algorithm to better match the CMIP5 mean. + +### Hector V 1.0 --- +Hartin, C. A., Patel, P., Schwarber, A., Link, R. P., and Bond-Lamberty, B. P.: A simple object-oriented and open-source model for scientific and policy analyses of the global climate system – Hector v1.0, Geosci. Model Dev., 8, 939–955, https://doi.org/10.5194/gmd-8-939-2015, 2015. -[Calibrating simple climate models to individual Earth System Models: Lessons learned from calibrating Hector](https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019EA000980) +* The original manuscript for Hector v1.0 was published in 2015. Hector is an open-source, object-oriented simple global carbon cycle model. The model takes in CO$_2$ and non-CO$_2$ emissions (e.g., CH$_4$, N$_2$O), converting emissions to concentrations where needed, and calculates the global radiative forcing and the global mean temperature change + + +## Hector in Literature {#other_manuscripts} -> Dorheim, K., Link, R., Hartin, C., Kravitz, B., & Snyder, A. (2020). Calibrating simple climate models to individual Earth system models: Lessons learned from calibrating Hector. Earth and Space Science, 7, e2019EA000980. > https://doi.org/10.1029/2019EA000980 +Other journal articles that use or reference Hector. -In this article, the authors use Hector v2.5.0 to emulate the multiforcing historical and RCP scenario output for 31 concentration and seven emissions-driven Earth System Models. +1. Wang, T., Teng, F. Damage function uncertainty increases the social cost of methane and nitrous oxide. Nat. Clim. Chang. 13, 1258–1265 (2023). https://doi.org/10.1038/s41558-023-01803-4 --- +1. Protection Agency, E.: Report on the social cost of greenhouse gases: Estimates incorporating recent scientific advances, 2022. + + Specifically the [Supplemental Material](https://www.epa.gov/system/files/documents/2022-11/epa_scghg_report_draft_0.pdf) + +1. Wang, Tianpeng, Fei Teng, Xu Deng, and Jun Xie. 2022. “Climate Module Disparities Explain Inconsistent Estimates of the Social Cost of Carbon in Integrated Assessment Models.” One Earth 5 (7): 767–78. -[Evaluating climate emulation: fundamental impulse testing of simple climate models](https://esd.copernicus.org/articles/10/729/2019/) +1. Woodard, D. L., Shiklomanov, A. N., Kravitz, B., Hartin, C., Bond-Lamberty, B. (2021). A permafrost implementation in the simple carbon-climate model Hector v.2.3pf. Geosci. Model Dev., 14, 4751–4767. ttps://doi.org/10.5194/gmd-14-4751-2021 + + This study adds permafrost into Hector as a separate land carbon pool which decomposes into CH$_4$ and CO$_2$ when thawed. The thaw rate is controlled by region-specific air temperature increases. + +1. Nicholls, Z., Meinshausen, M., Lewis, J., Corradi, M. R., Dorheim, K., Gasser, T., Gieseke, R., Hope, A. P., Leach, N. J., McBride, L. A., Quilcaille, Y., Rogelj, J., Salawitch, R. J., Samset, B. H., Sandstad, M., Shiklomanov, A., Skeie, R. B., Smith, C. J., Smith, S. J., Su, X., Tsutsui, J., Vega-Westhoff, B., and Woodard, D. L.: Reduced Complexity Model Intercomparison Project Phase 2: Synthesizing Earth System Knowledge for Probabilistic Climate Projections, Earths Future, 9, e2020EF001900, 2021. -> Schwarber, A. K., Smith, S. J., Hartin, C. A., Vega-Westhoff, B. A., and Sriver, R. (2019). Evaluating climate emulation: fundamental impulse testing of simple climate models, Earth Syst. Dynam., 10, 729–739. https://doi.org/10.5194/esd-10-729-2019 +1. Nicholls, Z. R. J., Meinshausen, M., Lewis, J., Gieseke, R., Dommenget, D., Dorheim, K., Fan, C.-S., Fuglestvedt, J. S., Gasser, T., Golüke, U., Goodwin, P., Hartin, C., Hope, A. P., Kriegler, E., Leach, N. J., Marchegiani, D., McBride, L. A., Quilcaille, Y., Rogelj, J., Salawitch, R. J., Samset, B. H., Sandstad, M., Shiklomanov, A. N., Skeie, R. B., Smith, C. J., Smith, S., Tanaka, K., Tsutsui, J., and Xie, Z.: Reduced Complexity Model Intercomparison Project Phase 1: introduction and evaluation of global-mean temperature response, Geosci. Model Dev., 13, 5175–5190, https://doi.org/10.5194/gmd-13-5175-2020, 2020. -The authors use fundamental impulse tests of CO$_2$, CH$_4$, and black carbon to understand the fundamental gas cycle and climate system responses of comprehensive and idealized simple climate models. +1. Evanoff et al., (2020). hectorui: A web-based interactive scenario builder and visualization application for the Hector climate model. Journal of Open Source Software, 5(56), 2782, https://doi.org/10.21105/joss.02782 -### RCMIP Project -The Reduced Complexity Model Intercomparison Project ([RCMIP](https://www.rcmip.org/)) is the first systematic intercomparison project for reduced-complexity climate models (RCM) and provides a standard protocol for one-line models and simple and reduced-complexity models to perform experiments from a common set of assumptions. In short, this creates a standardized test of RCMs to reproduce Earth System Model projections. Hector is one such model being used in RCMIP. +1. Vega‐Westhoff, B., Sriver, R. L., Hartin, C., Wong, T. E., & Keller, K. (2020). The role of climate sensitivity in upper‐tail sea level rise projections. Geophysical Research Letters, 47, e2019GL085792. https://doi.org/10.1029/2019GL085792 + + This study uses Hector-BRICK to investigate the effects of Earth's equilibrium climate sensitivity, or long-term temperature response to a doubling of atmospheric CO$_2$, on sea-level rise with a focus on the high-impact upper tail. + +1. Dorheim, K., Link, R., Hartin, C., Kravitz, B., & Snyder, A. (2020). Calibrating simple climate models to individual Earth system models: Lessons learned from calibrating Hector. Earth and Space Science, 7, e2019EA000980. > https://doi.org/10.1029/2019EA000980 + + In this article, the authors use Hector v2.5.0 to emulate the multiforcing historical and RCP scenario output for 31 concentration and seven emissions-driven Earth System Models. + +1. Schwarber, A. K., Smith, S. J., Hartin, C. A., Vega-Westhoff, B. A., and Sriver, R. (2019). Evaluating climate emulation: fundamental impulse testing of simple climate models, Earth Syst. Dynam., 10, 729–739. https://doi.org/10.5194/esd-10-729-2019 + + The authors use fundamental impulse tests of CO$_2$, CH$_4$, and black carbon to understand the fundamental gas cycle and climate system responses of comprehensive and idealized simple climate models. + +1. Willner et al, (2017), pyhector: A Python interface for the simple climate model Hector, Journal of Open Source Software, 2(12), 248, doi:10.21105/joss.00248. -[RCMIP Phase 1](https://gmd.copernicus.org/articles/13/5175/2020/) focuses on RCMs' global mean temperature responses while [RCMIP Phase 2](https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2020EF001900) explores the extent to which different RCMs can be calibrated to reproduce knowledge from specialized research communities. -## Other Places +## Other Hector Tools and Resources ### Mimi Framework MIMI is a Julia package for integrated assessment models. MimiHECTOR is a partial implementation of Hector in Julia. Github found [here](https://github.com/anthofflab/MimiHector.jl). ### pyhector -pyhector is a Python interface for Hector. -[Here](https://joss.theoj.org/papers/10.21105/joss.00248) is a paper describing pyhector by Sven N Wilner, Corrine Hartin, and Robert Gieseke. -Github found [here](https://github.com/openclimatedata/pyhector). +pyhector is a Python interface for Hector maintained by external collaborators. Available at the [pyhector](https://github.com/openclimatedata/pyhector) Github and described in the [pyhector JOSS manuscript](https://joss.theoj.org/papers/10.21105/joss.00248). ### HectorUI -HectorUI is an R Shiny web interface designed to be user-friendly and to provide an alternative to the command line for running Hector. This allows users who may not be fluent in C++ or R to explore model scenarios and outputs. +HectorUI is an R Shiny web interface designed to be user-friendly and to provide an alternative to the command line for running Hector. This allows users who may not be fluent in C++ or R to explore model scenarios and outputs. Access the HectorUI interface [here](https://jgcri.shinyapps.io/HectorUI/) and read the [HectorUI documentation manuscript](https://doi.org/10.21105/joss.0278) for more details. A YouTube overview of HectorUI by Stephanie Pennington can be found [here](https://www.youtube.com/watch?v=fBHXS7pjZcI). + + +### Matilda +[Matilda](https://github.com/JGCRI/matilda) is an open source R package that provides users with a reproducible-transparent framework for conducting probabilistic analyses using Hector. + -The HectorUI interface can be found [here](https://jgcri.shinyapps.io/HectorUI/). -The HectorUI paper can be found [here](https://joss.theoj.org/papers/10.21105/joss.02782) -> Evanoff J.E., Vernon, C.R., Waldhoff, S.T., Snyder, A. C., and Hartin, C.A. (2020). "hectorui: A web-based interactive scenario builder and visualization application for the Hector climate model." The Journal of Open Source Software 5, no. 56:Article No. 2782. PNNL-SA-156861. https://doi.org/10.21105/joss.02782 -A YouTube overview of HectorUI by Stephanie Pennington can be found [here](https://www.youtube.com/watch?v=fBHXS7pjZcI). diff --git a/vignettes/articles/ex_GMAT-GMST.Rmd b/vignettes/articles/ex_GMAT-GMST.Rmd index e1eaaa0bc..fb900e9d9 100644 --- a/vignettes/articles/ex_GMAT-GMST.Rmd +++ b/vignettes/articles/ex_GMAT-GMST.Rmd @@ -130,11 +130,6 @@ ggplot(GMAT_d) + ``` -## Conclusion - -GMST allows for measuring global temperature fluctuations and is used to report -global warming trends. While Hector does not directly return GMST as an output, it -is straightforward to calculate the metric by using land and ocean surface temperatures. Remember to shut down the core. diff --git a/vignettes/hector.Rmd b/vignettes/hector.Rmd index 74db7b4ed..dab5890e2 100644 --- a/vignettes/hector.Rmd +++ b/vignettes/hector.Rmd @@ -204,7 +204,13 @@ However, the effects only become pronounced in the latter half of the 21st centu Hector runs quickly, making it easy to run many simulations under slightly different configurations. One application of this is to explore the sensitivity of Hector -to variability in its parameters. +to variability in its parameters. This is an example of a basic sensitivity analysis, +where we change Hector parameter values and look at the impacts on Hector output. For +users that are interested in more advanced sensitivity analyses or may be probabilistic +Hector runs we recommend using [Matilda](https://github.com/JGCRI/matilda)! Matilda +is a package that provides a probabilistic framework to the Hector simple climate model. +Now to continue with this example. + The basic procedure for this is the same as in the previous section. However, to save typing (and, in general, to be good programmers!), let's create some functions.