From 9493b7ffbaecb9e77ecb4c5fb2905a106b28ad9b Mon Sep 17 00:00:00 2001
From: Lawson Woods <ljwoods2@asu.edu>
Date: Sat, 15 Jun 2024 21:15:49 -0700
Subject: [PATCH 01/13] initial commit

---
 kerchunk/hdf.py              |  34 ++++++++++++++++++++++++++++++++--
 kerchunk/tests/air_linked.nc | Bin 0 -> 75948 bytes
 kerchunk/tests/test_hdf.py   |  33 ++++++++++++++++++++++++++++++---
 3 files changed, 62 insertions(+), 5 deletions(-)
 create mode 100644 kerchunk/tests/air_linked.nc

diff --git a/kerchunk/hdf.py b/kerchunk/hdf.py
index 46c325c8..0339be4a 100644
--- a/kerchunk/hdf.py
+++ b/kerchunk/hdf.py
@@ -113,7 +113,7 @@ def __init__(
         self.error = error
         lggr.debug(f"HDF5 file URI: {self._uri}")
 
-    def translate(self):
+    def translate(self, preserve_linked_dsets=False):
         """Translate content of one HDF5 file into Zarr storage format.
 
         This method is the main entry point to execute the workflow, and
@@ -121,6 +121,13 @@ def translate(self):
 
         No data is copied out of the HDF5 file.
 
+        Parameters
+        ----------
+        preserve_linked_dsets : bool (optional, default False)
+            If True, translate HDF5 soft and hard links for each `h5py.Dataset`
+            into the reference structure. Requires h5py version 3.11.0 or later.
+            Will not translate external links or links to `h5py.Group` objects.
+
         Returns
         -------
         dict
@@ -128,7 +135,17 @@ def translate(self):
         """
         lggr.debug("Translation begins")
         self._transfer_attrs(self._h5f, self._zroot)
+
         self._h5f.visititems(self._translator)
+
+        if preserve_linked_dsets:
+            if h5py.__version__ < "3.11.0":
+                raise RuntimeError(
+                    "preserve_links requires h5py 3.11.0 or later, "
+                    f"found {h5py.__version__}"
+                )
+            self._h5f.visititems_links(self._translator)
+
         if self.spec < 1:
             return self.store
         elif isinstance(self.store, LazyReferenceMapper):
@@ -247,10 +264,23 @@ def _decode_filters(self, h5obj: Union[h5py.Dataset, h5py.Group]):
                 )
         return filters
 
-    def _translator(self, name: str, h5obj: Union[h5py.Dataset, h5py.Group]):
+    def _translator(
+        self,
+        name: str,
+        h5obj: Union[
+            h5py.Dataset, h5py.Group, h5py.SoftLink, h5py.HardLink, h5py.ExternalLink
+        ],
+    ):
         """Produce Zarr metadata for all groups and datasets in the HDF5 file."""
         try:  # method must not raise exception
             kwargs = {}
+
+            if isinstance(h5obj, h5py.SoftLink) or isinstance(h5obj, h5py.HardLink):
+                h5obj = self._h5f[name]
+                if isinstance(h5obj, h5py.Group):
+                    # continues iteration of visititems_links
+                    return None
+
             if isinstance(h5obj, h5py.Dataset):
                 lggr.debug(f"HDF5 dataset: {h5obj.name}")
                 lggr.debug(f"HDF5 compression: {h5obj.compression}")
diff --git a/kerchunk/tests/air_linked.nc b/kerchunk/tests/air_linked.nc
new file mode 100644
index 0000000000000000000000000000000000000000..c8d4136e1ecfe1a3265d9cff71c1c66b689adf03
GIT binary patch
literal 75948
zcmeFa4OEqPp67XkXi<b%#Mne_23eY<XtXIz(`YPXC`6m6MWJCZOQS)!_j2Fw6$!B^
zW6^3e=-Oz}YJ1RT(Q=7#DdSqk<*4IenpK;@6obYl#wF+>n&x1Z#T1JfTle$*KVV|I
zx@LQ3ws!lB-1G2q`QPVx?*I4x{k_~R*;$L8ikuc177<|^`+6|?)tvN%c`9}EAKh&~
zz51i<tmR=(bv+T5V4fToIr_h`*L?G5WSGsIBsy9Pr;W3gFB*MS)OYpEI9rr0!WM3`
zEeHlG9D#CoFtFD44F8zG5yr9LM?YAe^=j0MHqJGUXU*brce&4Dmc@-OOIsIjU+wn&
zX?)HP7sQu2RtHx5{&cN-ZG6(K7t`kae$;xM?BTWSj<qFa?vO5KPVy(CC&^k}zB-;4
zIm?26b0T-ZUG84(b5}a-Upr}JQrf!j*;oJR@2AAa=eZnf9r3Hn9P#BYM|^6^i!c6B
zYW!Mvseg66&r#v<#aFEMz2R7!67_p}uj6drw%$O<x4PUJEc2Voj~~7K4}*4xFRH}G
zc8s$<X)P>sI2~n<K#9ZJl?9_qUFGGWx!?QVn{U3E;#gbeO9=#5uTCiqR(vnC*3PY!
zul`<XncHr6SbwDO`cTSOyOQE6_xsEa+t!-jx2FH)cP#4Ko_Jz1*RX{>@x&8P>KDGo
zevX?gdOjULUG(BB-u7&S?b)C3+xSTRXI!%NeCtnsC)O;1L*BNYvc0G`$=B<nKO-jT
zuW?U&-(1aRo2ftfiWzTSDazZ{=-Y4oYx|RzfApgT*-OUW!vxzm`FiYci?9Fg+tw}q
z;E_AMci`@e){g$}*LD;(G+$#s|8sWKZ+6aR`^hL-{-@EOV>^0kbVmoFh`zqQquZ>f
z{{Qp+&%A&2jN5LDMgfiSlt^#g`YXx*^MB-zMa#2u5e7-HQ6$#gj}hUE(YN{Mql`6r
zp0V>;i}^mLtR6n`PyE7djI!+5Z6t9!`B7ub&%My}RAydw7I(y*KF!niq`w=zcI3mW
z{KB>Nv|06St$NB@^)2fhU#*(>;_XX|Ugcc2v3D{0W=1z)?3%izT-!QZ5>HqwzCU_D
zzdO2Q?B`F;KNFdmw{#JAm#&rQSSiw_hvK7c|6Xe|the{Y*6E41Ki3LgslEL-|2A{B
zZ3VyBW?0u~yS00s?MbcR3BA7y>%O<#woxn4xx6Xc>;J=_ZM8M(306qCJEG?QrER5N
z!i}0Y>u&g8Y>j%tXG8PBi(lM6WsBh#+Yj{3W-Gba`S;)67`Mav>+$#h8arU<|LsS%
z`4P%KLc~W1_y~(ULi|U#=}{Vgor@kJ_Vy1~L-xmCzwI56@v=-);*-|3Ge&P=f)TwZ
z9_qJocp(!d%KDO++H9UbtbOwzdjh{>{%5m|{ru0dTM_@Ab&`AX?)a)7{d>%Am5$o2
zJnP>!8%EN)r?Dq~+du!o{?lK6wV55Gn;9Iv4Qq>wWdZp5f%U^ux1+z;{&K;#MX$cP
z(AGNoR-d@@gzYJ_-@o@a5gWc~MaIQjZ}%(v@YJ{e-T7}?k#PzCXw{MIn8UiysD;wh
z<7^Yr=+OA}oCUdM!H}b@+&$_&CU9fAe1)TIExzwdFU$-$%1i9dm(o&ROnG6Z%kKQ6
zmu98RN=x}e>&=ZD4FcX7dvoUJ*`Cf_=539c+$XSJQRZzmd>Kkvl%19N(CsrX`0L*-
z{_Z&2WITOH{CbYLWZHsYu*~ibtS;x?9x?wop4K~jl>85U{{>5b@an6^Z<aTp%Vxdu
z!-YAozVf4-*I!-mgQW}Ef$`Q3<ow`=3(ZP#yjH^RW`Clr{ju&O%kBR)E(ohx`Xslk
zn_lGh`SRt-|62D~&hkIL=hmM*B*ohD)dBnJGW+X+)qX2(vHTfxNAB`B?2d7`4u)#J
zV4(CLFI&2rkjE?=jey=DGO<DeITw0wc1LNM!?E`DK(MUb^>7PTuX9_=bbng*-;c9>
z$Cj!$wS3{Ch07P_ELiyZ(pO&1v!(N2?Wm%MwT^O|?OFcNZQ9*_M?m~yi_rfq?GbMx
zgm1p-d@bN>?0M_2$KU^3+JT4GBTer*Y}^x1j72-63iP2B`KO*A7bOOgQV_$6E5}FZ
znQ_Vb`-!SCi}C$==HIpppIAS~BDU|?;zn)E(b1n*Mt|Cz*2>?;_e;(rnbOK*MpquP
zzQ(O~myOBE1j(4~iE#PX`U)>c)`~_~7XRtr{GYsj?A?w9JY#WP6wZ|ph5z@Z*}6CU
z{^jm`y`yj5QQqi+-*`vAIf}bAsxiJlcl2CiKFZk7(Wn2*qPTI1V~T0zJFi3-8#Z>q
zC%A*DhTR^_FaLMC+Hbz$lcP8E{n*i8#%^@%=jhYl8U6X6`oYT{zS6o$g`M9wGNjj@
zI{6RwZ#oQrF?#LqzysD-lE1zAn-0TY2qNR_t5v`IcH5tP^IVe>tW{q*A06TQZdlL#
z+V8tFZML6sY~%Ph92+?*IYJyxjzW&*9N8R|3v5;Efv$eCz}8yD&y5Rg8Jl=+Gk^a%
zfBzf)uHo-3{JoXGxAFI1@OLeLZ|Coy@^>A7|0RF#<nP@JGTJ|VCB)Il(aB-^_n*F!
z%2C9zgQJyWkR!h1(^qmisyL2v^m0UB{PdL!ju1y9M<<8v176Qj#Ib{;m1B@2{+GO-
zql)7wM=wWoC$HxSaWryta@hU@ujeS@*ul}tF~||$#p^k$IF54kaztO^^&BCNMvhJn
z+htzQQN*!>qm^TjBmP5P&r!v3l%tm;`U<b-2yrxWbaL4Kj@NS(aqQq|<rw6M@8<O!
zRUAh-dO4y$;`JOMjz*474%<~;&r!s&gQJyWkR$$MUe8g*ag?K%Bl;Sz=Lm5$a&&Uo
z{&!x_QN*!>qm^TjBff{%b5wC0<>=*z{)E?agg6>GIyr2;yq=?oV+Tho#~?@ibzaZ$
zYgzE|laF`c@eVxRfyX=WTi*c-wH|YV#=NSrpJTW4P1Jf)qV=CKw0i7k@q-s$97U~P
z*9z3VPGwN*W%}(8EL<h(TI|2F{azSPkVI58vHj@>aW?np6N1E7y)Zb_c8_0de`q;W
z$$uOFecO3G@ker(ZMK#5?|$F5K_{JIfz<!{hjFjkepmnEiHZC$F!uMfBEdqw8*ctz
z@7NCUw{7}pNFrSOAGG3mE3xqV!>^ySy`v`p+!ZsQ%J`$pwoX04tqb|?I9+<f_O9*}
zw_PE7$zT5T?=6gW)VK)>_x{Ou`uX2~i(g+yt;fv$*O}o_0(_mKkI>upb;8-c&ODD0
z{Sjh4La0aR`UtZ<!f21s@)7p=3(tBn)8nsyF+1>^BRw!%^EC$IK8({I#%tF9{|iIY
z*wM?_Vkg)C<eM?v_vegVW9%x={m1%WSxR{9{>ENE_S1TWwBVzuL;ENUGp2mUe!lYe
zw>}<K!n1Xbq~DY}tkw!%DP8v)Qin7DS*gR<by5`QxA<X`irud#(3eu6zczJPf2Rrk
z`uOX=)(%*Z$SADv1m=B_QDs%X|37|1rTy`OA1xJe3}Va=Su}l#AM{l;x?c8FUt`aY
zWt+wTi?QdYjPbR#+;(Ve!RX(Qe>~oS$2;(N2OjUh;~jXs1CMv$@eVxRfyX=Wcn2Qu
zz~dcwyaT`e9WX_HV^uch^8eevuFv^r75R-J3S&PFoWQAsk^d9baLJ?pPFt5{dVPN2
zE{iXB_$hH*UH%3wIIGJaX}zHr=V8_0SM_p&?Q{+5<(B9=rYbR1<|uKS8n>^P9sRnx
zKD7}=STFkH<*#JEK3Y1cGYSp)buW4HtM_AG@~{X{Z(T3U`(f_G&Y+6bWv{D0$<!?h
z1<`_UyiwhODHa^7)l+$spcbuF1Jv?|4NwmoiA*=rS7#Q+rq!>$QRZA-@>TgEm5Sx=
zwdHQ>`eWyOC17_~xb1JO_Qn6Gc<sYJst|oof%4b&nnxO$enc-+**HKUQ#`r4q?{MN
zUbZ?=`cIcn>R$9U_H#Jt#RV6m_AhvNmbE3TeU8_;*YaSQ?dkZJM<JyQTiTPhv<K|u
zuiN<_Pd@U;*DE)_tSYH!)jWPxDr#yTr*cJ=ERK~r>KtSIIQCO~`Tfzvm0+%~|E>M`
z`ro3Lia!4z)R&I2JFoZ``1SQ{%@0QJ|Lff&|EwN^QIy#Bb=cYV2pIfzxcL!O__YY<
zW2Er^1OBm64`Wq(W4t`}bL@t`nR<9I+PXHDav1yhuRKqEI`*2!KmR4{z%~0d`$Wf&
z9bY(n&Q9l(rOBn)rOwjU(xtA+?m_oiPmZ_Pd*3tBbJ3OIc-E2O__-s;8CJTd^j_(B
zS9|I5(mv<g&X1gxrEb?ccbR9BH_Lm+ljEH0%yDjY_BpdlPp~f2wbYg3`mD6JG_CZ>
z(s*87<SKD*_iXXT^ZHI+-^p2uovF?R&eNQ0sp|$;zu_u(jdu-|mb-SjBfPopCRd{C
zQfV`1Sza1e`jIoHw7s;+waC5BJ<r|ka=WIvX1j)686Kzas&|q%&ok)0>z?8X_ndL>
zavye|cb{>WxF_={Fn`~4S9>A@!+wwdfNwgFLLOnhDsQa!tS7;<$i3H9=ZfPsGd$Iv
z`<{;jX95L*Zhw>i1OGFDxImkKiGRqq$JfW|<L)HRy2#z=&hSKd>w`OjF9&A^R|VH|
z3<d%LtuNq6^xyW*Wp|EvT%LMQmiM*L>CjFdS3}1_Uj!3^+X6|f+8Ef)*+af#-b3C7
zZw{|$@FuT0vF7toR;VK64OO%5nZR5AR^KMyJzt8y$9KTj>TCBk^SJD5w0GGv9376=
z*vAp*G~)amng&a!xgNOgy03UvK+P+jF1OpY+kV(S({bPNxwEx2-&Nx(Vy6p7sTqU<
zx$bn&S<e;StD_$7sB*@a?l-i|hnk7dzOeLjsCEoGTq+%KC3v3m_Bk7zK4|!bquMdi
zan1g$;}b_%X`<^QblK~QcGbdJL#4CacRZ2oZ=$QY^rg~H=pO5Qp8IVt?RIU1HjCUB
zxxx_Fy69TTStohxp~zh*u^w9F@w=#cm_0x0YIDuw6^&4Lz2`mlqlaC|_Bq*$R$rwr
z%6G{7oOiKjvU?ibR|C)8<Sg}`8h^JxF0e3=8<@vlcKd7mnf_>hyKlt1(o^RughOlL
z&;rgm56T4s9|v9zrU#F+w#J|0p9<~g_zQjaJ?GteA2smm3{QA)F*IBi><zXAFF?UL
zP$=7%#{KN!`yuGq;aTZf?CJ2#40Z(91YN;Nq0~@#u!%Dc!GBq-==3IY&m(-_#QhKX
za_#f%0Vw*4<9$c0^EGD```hPS&d%I$b-52BZFks@B=?EZKeMl6-43V}ES=3>??o1x
zk>x(9zrVB&3374edhYn=oa=e0mJQWryEee(HSG3hP^ORbzkwv~FKu%Tx;s20$mKbn
zyYE=#IK%mWVXt<abKZkW+WF7erB2S$<XZ1J$ys#0n9{A#@)N@?LF9H8UH9{#=4{sq
z?xKjZPV{6zziaMm?gVJr;@Qsq&p=WG?lxCGcN1r*dDy+!lj_~ceq8nTu^W}}xYHNo
zyXxKIJ%ViBLkmcQ=zNPk*}nIDOZ+GOZAj-n9yR_V9{K(~zFtqaD~QHfUiumDsmq<u
ze{E3I4kh1Z#RUIlUn0+yL(MzfMUs2As}ik}<c>rRw*@W)rXax!1Le>%#68?Y!eV?|
zxQhz!LC-eUA9u$EZo(f=1>X;TEA(<OGLY|2^xuPiSG_BsXfe`zm3M$x_<mu(fJ9Uy
zbE4!{_WBLz8HTQHay7zhchEEIJ-ggRt~0EB#ql)%>nYs;k2gV~p3+Kav>%#9vU}^H
z<`H&Mx;L9UUg|pP8sI;}u5M_N33c|cTI*AI{tTMAnLRtfv&+#5lG3echFZ>clK;1H
z{`S(%+{sZ_xM!H>+^&1PHppW$KR+wY<Q^U%Yq{=OM%QK{o8hDKI0fsl*mKa+>*+T<
zmVw5d<sN`bN}%R0bZD8U%6kk-rXdTr%r5Rh+a~)u(Y8}OI%6apa2)B(^R%<iJ-#}B
zZs2&}9b_*F`P_$mzQ^hoq<#S2J_%<Hvh&xG<J`cP0T=R=ht3TT-h^fw(dvo*>Ap$c
z8SLgZY{1RHAZPnBa3uI*=*7^7W_2}E_%kGSW55P4Wcq7;r=VDM@CcHg7fcO}@Tdy)
z2kVjNFM=N;*&m~s8_oIKIR7e#3tRX!l<Y$@rXg#|XtO=+cO-VF15LXEO?%#*@2ZF9
z9gbpWGB&e`|3)G`l88z}oyqWUFIF}Sj@mrB+q03$He|dTs+MCpdpK_;I_WdGVIJ@2
zif0<ybUY8yIl1&0oV1g>=t0L_<TaVDRw%az34}k9KR3EotM9Sm1XNCB&17h#JIUni
z?MTxu^hvcR9Vxhp=C6m69mwZF^#3(hws9}>*h^W`#hzH+n<zTsoey7!e9L_Ep{A_v
zO{{hecfJm)pXZ$|gU%tOCLj46K=y3V<~qM8`bWIgo_X#K?1UT6I`1Cf*(F%xVgE)X
zN;Fu^kr(U=eBf_}_V>~GxoDFv_j;)I0_WL{R<4GY=}`Vi@cGaVp4A!NMfyd}x!$YT
z%Y<MfXZsMnTn&Gm2_!(zHNjkHZbJ_QjAbtJ4`MZ2@N7Dv(sJk+gKZ1L!bNj$6Oqd*
zB>t>tqr0b6vYCN4jz<ei=1#JM`_YA;JJYaj53qv~tX+Z4+6ixTItzKWoK+V&>X0+(
zu0rEuT{0GFG3$d|Q;&Ek8S6Y~bi_x_9=Jv4PvQK_q5cN+XFhbl1P%A|kTnfLyK*CG
zi`aQLR&x)}&vqS#PK|g*>!5;*EBK(7-otZPxLkJ~cQV;6Te}Nd4{}taeX6|Ek<ZIW
zUk_5+>^tX8N3%zxmm^ts9=<!vnY8M%uN*qrk<T+cb_WL7?-o47O8l})-qjhm>|V2P
zg8wbNr8CIq5oB#W6up4Vw)-l)*YH4ExzAj8roRrFy#;r?32k;jr96%!(6}3!kHRDC
z#YRlwp{q$&<B-8TR;~(8;fzN(k^<5GO029C-LeFoUVzp71>Vg;<g*w}w-b$b4$AhF
zUUXf<r_D0^FziY$y@nU`sr|=}1<>;nS}v0vlqHm2t!4kuyIb%;?x1NKIFHUD&o>iq
z$qg;X8xQ6qtZAq7CFtC3ta3AJbYz#FhMM*G)}J_zIh&!=MdJkp;jZPJ<2c;E6tA@h
ziPZWae6g9c*WoQmipIP4z=_%LnWX73w3M$a>6cHQ<xS<jGO+%$;i+<Hl*@l*Tjlkp
zV_8$ZVTJ~YXtzwf+_~OH<ahu{pU7@5#&5lbEcbeRth?p21-b$=@m$l<p#y%W@2n?a
ze49AV*@34~;f+F)TcP~6z#4X4e(oA1Y9TV$=^ey7KM8+sfjgY6O2G<01D#%mO7-Z^
z4k%xTuUd(gi12PlvlPM|x45ncIeZg|*R|KM?tLV_2zxvqE_(~jKZM=<)SidcY+&~*
z*{51|X%Dn4M<XS0|0~#wdhWKTv>uAbLb((C?*@{45jrXAkWJs<3Lpuo@NW*9r;umI
zLybuGR5tAa`fV?qE<dpcsduwa8##Nurxm`4MoTI72pW1cqvKEV$|PQO(nz#u)(A%@
z@@m=G{ai1aYhN@vvk|>n!Zj4nBywy-me%3D$}5&%7Rx=I!b>>Coz)|a$I)@up!YWD
zl!13N1@5aJrDimg^x&@~vO~pKj56rF9%{=5hqGGo+Es6*FVi21|MCv9xeYJ(BpS2<
z4cfu~rP1r5=T&brRC>#w#Lj;!^jatuyZAn|x{1CW!lP0|5rGdj-8YprlD%;Le=JnE
zrgBZx8gHlrt@{qPwvlVZ1r8bA(d6%f;(2J>3SyemSkRa?ifv{>efdZ?S>J;H`30J`
z0XluoL$TlI__UVvv~d49*m-H%Fyp;CNBz_dP-h6rCnG)5%ksHnxxWflZ@@blV6WGq
z(@LO=q$dUaQi$&t2F=>hNd-v7LA?8u@Wm{sF#vroA?t-uYZ_FkgXYtCO|9XT!}uh7
zc_yESWIZ3vAbOS{iHai!cttdFw*k6O!OIxt?;?2aJXSQ7dphJjhvzrRTL%3yjNhMO
ztZ5!`%yDR&00m;$<3!>FMH4;P$tZlw>Bz?lBxEHtbs_0<;f+qOlXWN2d$~}v5v$dP
zhv{T(l<&5QD#}=wjVwiTrVsq@AfLU#eiMnL2MdT<mY`!#`h94<2)Jx2XOee&Ca{<o
z<WT6ZLz6--D3ru`yOCO3pn`X>0_&ddSL~CHr2h=sKM$9DE3_EfSwzfo9X(JG*ouv7
zFcci&`%bj&0(_AX=cUqN*A(cZ7`7KrxwW)`=Z|3{Tali8<moe{DXjDy^e@DF4S)s|
z<5Bp~4e?N{6_4Zw(Lkj8B>O6xw+BkSgx-`#dl4C$f&be&Drat1I60T>Z7tDkH{7t-
zaKL6XYdQ8v@#<!%EsMDs8ecN{AOlIrLk|vf4f#fHbn+r-JDW%~h?kU$pS{9++gs?X
zg+AxlzbnwB7b)t1Po+gI|7sQ%v;r#iAS?3zn~irH1v;Q8{;KylJQI$_Eklm`Jv*UQ
z4|Yme!gU^XSghNlbWZTxgrAOJp-v%Dm8|plUoakSHC8d5y`SKp&-01?BXIL}?&TPA
zT+UfUQ_(XAISvm#1KO|*JVPFL3;v5D`95z7R(2sA^pw%G@`xItWj3623thX!AA`Ke
z|B+NAgRpER=6nPHwiU}Z5sqyz@mG>NjlF)G=L=0_P{e+<W2ZV<FC9INecO&VF~YlV
z@hnICN=|N|?JgoEvUlxhNj;<`Ymmm9tiNM)@({crJz9%DCBE7WAIP%BVP7PjpOv13
zuY@9nnYgZ*XNEX`K2j~LW}>GN%~}IB<y~xa%c8d%zDa?{)?>?eV(sQ5X;tXMJIIX>
z>oJM%!VzY<g*a_Mic`Jwk<Msh{dQJgg@*Il$vt>A@A(GM0@CRP+}9cR5h$>Um2cs9
zo`j|zw5Mcq6FNAFsOCH|%~|Ml1>b24oK)^V9@v3~Y{A|Y1m+N>-t!InCv){i^nC=o
zR*B`C1EmH7lR$@_$Ah&6?tvvF6B#Y>$FV9LF1g}KLw|Ss8#!B6Xltk#jeH!c>t2U^
zA8;><NnV1Mwa88-zkLqIavD6K+28@|;fe^ai|A`JJ2`@6%C6Pm)9%Hl^%z^X911ld
z6K6r{It}e*3+oI;a@prL)|c}Tif|MP&oi2>o_(F<T@HPhqHWut=VwOt^tbj-F<CB>
zeS`BHcGr94$y$#S&{D6N$G#>YnHfk}m$8*M!1;u9Sl)0gJeGv-f`?CRJs&=AKwl{)
zk#*1EUc0zT5<KIE=dMA`H0UP3R;bW?sFRFj2=%M<N#cY@oab;M?bA)9KY{&|e`^mk
zA)QVrKh>}3d?k{fN@S4(A4xAnV;ffmbI_u%h1L^G#$gR#4o>&wV4Hi{wJ_**0BODM
z|1yvo`X%dPiBrOnyW{B9mS7=hPb@loF41u<GWe!BTk;ymnyYB#YP{e<d?O(Ym)IL2
zoGI|OaI=rlwmY$$dyMZPZ!Z<UCdvI0x~&J@rI_Oe@4gOc5)vdkn~bz5LXq7p_OAEH
z>UFb|i{N}Abn-@rcy&AH7>`aIaBV}gT|sIRxo$Z$K8o&cC))iS)Nj9$x)P`vV6O*|
zz+K4kQLZN|y2r@3?mD0IY#jAzO1Mr9R2k+j@~{hg;k0Jn*8nmndpC*Q5JI{g?IvB@
zg1$V*Ugm)*Jb=%R64~7G+(v_@LG|10?{t2PAwJv!?W53j)zBz`)$*QV`F|nQJ3##T
z0@TdKSItJ6)8GWrr~?{32Q{DLA<r=68}P3R7UQ`Thw_LfyWxTG;B|aw+0F)dV>)qV
z6<1#vs0vlCd16g$Xbn<#9X~Y=iVoxTZUVV&fH%_MqanQNj$q}Qzg=@0T7C?Y^e!u(
zfvcvmyLEUi@**ms&1=pGr?7-<_`k%x5($!T-|gxHjcY~bti9i0XesS0{OEK1+>7XC
zS=kCC*~yOQv$MHIx<n~aQ#eC5k6L&_7_j{1F#hgAvu<F86TQ6?D_luLA}K3?mXg^T
z&Qp)3Z2-xS<=&p>=i6wB&%lTZkk4yS^f;?;z&DeT%}iHz>1)KUJMiuHm&!XkYw(~l
zw3R$sMI?$JJG}$&nLICPoXP017-%#dT8^-5hu}+LA(F9EUU^~CAmzN5I3lqaDA^7z
zlki>M!A}|DRg!vPJ&o}2E-cA4?5C`W54-qg;Qiozyq7`nz6r!Mt-gJD!Vy^R#U`fT
z!gY?obBX?r;0W@z8~RMevlf3vqk%6&_YUsa2Y2~+#Q0A7e-_*t>J9D=?8DC2azAsR
zdA{*KrP0cHXUU+TsboOjW>4$Tp7qB5o&kf7fM!DD6|WuTHNxCC6ZQAuF;zmljaZYN
zXaXOo*iLk{6Pj80oG{u%c;z#s&WVK<&NLaHLYhMv6=6d|P`VTUE0tI7FWrk4l{YMe
zS~fKfu8J>x19bMs;JG{S!QVh1j7RbgbDnVY>_)7tVxTsVy&jN=Jck|B_n@P)bRP8W
zfSPxZ$nET%eA9X)v>)3y%sXyyUUBq0o-B=Y-}g*6XizQsvWUGrfc}d@;?AL`XTZ<8
z$87HYqU#8_Z8UqAODx|DHn0&-#EBeV#g9t!<#`^!SIWUW;BF?t858g~JMeF+!Em~<
z%F@$C_)a-qAs@TQy{zZ1W+46Nype%=gYgt`ePPBiTq!qD>6_w-gJ<gvUmoNgb$Xlp
zGl^FU&<nNLg<69nPh$`EvI{4X4@u_;lvI|Z9Sr#lS)*RG(p@Z(FrOH9-iihzks8@e
zMW^y!=fUyM!NZ+!Xa-(TA<t$SYIYmBS5zUp*~FEE#Y|*JBj9@JLm{=z@QFA>xi&X)
zDO7GTn(ZcBE^S(m#O>kUwu0r#ua<QkhTHPco@cqkBWS_|@X|)?WH<M9+Sw1yN1XD8
z<%Q-TYn>b`kO4g=5gSw>aTn3vpX1vNl^%A_MWbf3E0*7A!&@md`fn$)Duhd`m+~IY
zyNi*wBK%o9*6)j8EK$M-_@Bqvd11MQWH+|MW4n;f4$oC@yYF416c==l1!0{M9Donz
z!3P(h!-$FLp7ZuYbwyHpeD=W4@Rx-5D2{9-UKJjDCXk8mG#@lxevI^YnD4rO2eIrC
zD7qN$Y6^Dt<=}p7nXI8~m%Ka4+X!e+GcjQRuF+0S^xPopqimV7BIU@$EE7d){Uszz
z2-r1tJQn)&<CEOuIpJdRdiHYGrC3K<NXc0p^c0#s$of>TBF7?p&>pycr_rN3$<UNT
z&-1KMw9&>n7J-mA@%j_SOY3C!MbB~*v#v+tWYNyERuNGG9&R06`8HZRSbD>i&OU8H
zYLqvUPupT7wI7|ch0Iq89<b2NX|9{@L+s^cyv=fKRRI`}rRmP0Z{{PjTZlQ6z%c6I
zk7BPfDsO=yJca%gHLs&L)8K?me>%TAv7OtIfxCEy=dhYX=*0_A={`|ncrXqGp#blu
zI5db%H)DT0*|%6@WzP{)3WpbRw3~=id6t)hjx|-GkMXo_p(|FP|F^(h6>!N$yrES@
zJzwHg<$@A@9FVU+9;uXn?=;+>gIzji_H-NmeFR+W;(bWVD306@?e}8Wgrh3PYsc;_
zg+9t}6|-Nl<Oq_H4P{FU;fc>Ud!mW;YmiJKH1grkAbV%!qhY6G$!@;SAm^E6!>iGu
zCx~>fA$5z<ph*TXn+FC^SX%6CAw!@XK`QiAhEzDe_MrvKZe>Vvk;G}PV_<~Oa>mon
zy)Kuh0o${QM;bouBw}ajgdBEb#M_O$H{d7CB`RukpF$&+L!*7f03VRcdJmtr&*+;p
z-(s}qJ@!YwV6HpM+e#GqE>?38y?BJ3mq*jd+WAPQ4{N9#-BEZxgSfL3uhxTf?m#+^
z;Ir(*qge)pzeMXUgYHL+Pd*WzEJGezeY*qek<G<;LGNNOo8>1Xr%lLa7*b!2R@Qq+
zg+KG*iYfTK?<1K7SnqCsC0bFqcQ!Ut@!kUbj97e|m(Y}!7j&QNB)Nk`DCa<5RXHH6
zrUt1}eyS3R%|?4p#4_#V&MVMKK5WrJ#}&MgKI~^Z^qhudWg?k-4Lz0rzlNWh<$c=u
zrG1h8S0y&Gtc~^rax}-FW*dl30v4_dYq)_|&BO25>73}8Z0{~9u&;7#1!0i<PDF3q
z=h4ACoQNlzi$ve?RF+md=Go)GD>s*RxhvS=5HUswYjp@ZPU0x$E>B_OPU74Bm@N5y
zN22Sj=Q0snHyO+U{M1~ch&HUJ>I*ip^9RY6pJ0E&c<)m@$wUd)4Z7z7k?|5A=7E1t
zh5ka@Pk9rF=#*{P1{bbCzb7NJlA{jjxeA{)np|eS-xZt}H~<bWyiwSMY|9nST?gVb
zh(v}5hoPn_3?_wshD0}W|1D^6<p8pH2j}3@d%if3ru)e9w!jNKtjOU6_U|R^Ruz2u
zJo54p$oENbajTN#3Y?#UPHjeh?;(jBpk)yfcL{3lfm@4^w~1IMQF9AdJOO72?`wz7
zLiywkb|cTj$gtuTSwiXX2=C|Ew;$uz+(#xqM<3k7PYpmxQAu(nDh-<4)j4MwwF%+S
z`*Ww||K!1=IcS3lgR3hZTaR~B>>a|RU4>-s#A@boPsvcS8Hx;{2^DiH-rT}oIDJBx
zPV;{IoHtxq-g~|#kbrlwX@g)gPXz-Yu+d0Fh$!?nd;Wm;E4`5lA`@qF8a6PVNmPRL
zLC-IVBRtUmnZRYV;Y9M?l7%X4X$)GShZrTCUEfN?sa)k6Y~gw+kdJiwkj}*hZ-i&j
zpKW+OXP{;^lwXBqR6guOWI72f{(wP1vhb{Axs^R12EFYKzJ?$CArb98?ob(~mykUn
zJMt^PqNcQGG4xpnvX_Z0N^i>kDHl--#e@%y;DyTnZ2=oQg-=`UQ4NR{ohtTlvwA$1
zQP12%Ps@{dfb^||Vkfxf8=OPfZ6fP3k@{}vc^3*UhK>bjtc!Sky{uG)LNT&lgS|X%
z>}m$SkMJ?Yg5hAvxk!Af_ad=Q1DgI)={h7;+2}%Kvl6|z1*;)!T(5b~yUceNJubhm
zmRgl8_EFU#J@}Asu}jKEPo+Mm4-VC<Q^6XVxW72gxra5%aNPD@MmAnXf2JclbCBc|
z{}enTC#%!R(73otv9X5@CKhe;pFuOegAC;o*HnaNB7GD5F+?@8{8M-rD?kS(7!FE7
zk8d=Y%j1FNYqGGhqGm0$??)e3;RUsVHwd%PJ*2}M37omv_}yW6^$XCCKX$x=Wl~nE
z9iKgeyV;4wErTx=lUcF5sGo*qoJ+*e!rEedlBHy#s*(RHG}%N?%&5O4Kd1&8^%#vD
z<Vva;=z%ZVpx8vb@^r_i_A+Wyl$jKgp2;g_L6epEypqcjyr8{2ljm8*%4)|)&ZW?@
z3!nEccXtr&9*f;nwM#u-hw>-=9_7}2L`}o4GI%`(D#qYFj3A9d^i+4T#n{bkUkV=3
z=UDeTYQy@VruIg@{sgS1>O@kACKRPq`g-AKRX(`Uo|BMv)trpL4?S>z@{YN9O><bS
ztknpz)&(DmC+pFjcks8Cp)teYjQ7yy&tN%ISUrd~Pr(DOLhh$<_sX`;1u?vbEDJ?a
zOse{Zss0vff_8%?Wct&fys+I0e6vF6Y1#H#UpEqZ6WUiHiK+M@JMbH}*(cj?lthxp
zeBY6brk#P;sk*VnATg&|6$=H6*(p_J{h0Ni+P}2Br~pbv*DB|)z3N3P|D5PRapO`b
zQ3sW#A)C>BSEYe|=A(OOAzy>;4@%bAuQ;6Sr|L~k;yDQUQC*d?OX(&XauPir?pcRb
zZE()UmWLZYPJkDL0mxfd#Bl`a+y*tTU|~Ne-#tJE@hqC<6y8p<L2lBJrUrINw3H{5
z>{A3O-*}txJte^*?Bp$|sY<v5=;uvPtq)qPB<@K;%LqBl!b4VMD%nWH0@gy)a-xn2
zNV5E$D&+D!Xrz#eyHKbXtU;N!J)r)#p;0C)Yq{1_!F{~q6!#*3Pnhi#WV4-fT_)N|
z<<4ul7v;}e$YWlIzwW_nIn)mHgGh*)xA2fQnX^WN#eYoJ^=X5y%!NYF<K2CN?-F5<
zpiE>SkKIy5<wQ>#vbhyKEKfDr)SHBvIPFQSqv}3#xs$U<aw^xl>Fz`yrl2urVIz~k
z03P5a%|Z{Vf?_gMQ+9L_?>F2t7h8GKU;u;st%`x;MrK>s1KGiwyha(bn{MGc51`OW
zB9bG9<Av4<JJVHT$&I+6=R~Y=7&=|KFl8-I;rS)=P*t8H``bJt^ijzBbf0X$?tLY^
zz6{zOfO@q)A@6gLgE{Ek_t?KEXjFlxu8IX|s}W+3W>yWcms63lay;N+Bvh1b_2ptW
zRjpBN?o&Iz1$ist|Ee@ltwD{+hh6u_hHes}Bx7ZdfE<K_GR}p5C$Zb|z~17#C-KLW
zFKJ$L4nJ!lxwHd#LMzby&2UXOSG9pH#UYd7SlN2)@+*!6M>tY;#=gy7V9%uz;hMeP
z@wW4`(m41^I=Tf;bGth6cU5nvZF!9tM)WVn;>Ez-8`zN<?Eeb(XcFjwxLC2~MMJkF
zXt>Bc=Ha8NdP~@Q0eoA=%1G#^$jm~ug=|}ulZta^kR@;9Op|%8)gZx&7WJPi#!HGI
zew<|J*$*wdz))2!F%69nfV;|Aci(Vv8j_Tal%2y8<QUpN$8Xub^Kh1UPUu82R9a?e
zdCySk9`rhZ1QqhyJ^Y-Dr@GCc6VD-=+0b(;J1#%uUFajNtlVfm8c%j}67r^+h3(v_
z54&53ejPS?a{|^ilUFSaJWrfzgZAgJgsQZxhMJX7b0JY?4QIQI9vwnbp9%ae==%Dg
z9iJ%-o!^3Q<>cx$TwRYQetQ=Wdm+_9_o)ooN_=??TBTw+RXv%C1dd0)&EQ@lz@`rq
zV>j?zJd&sg@*YTBtFgB)@%=OOWC3gY*^vtL#uT#EsvDYRB&`eWcibT8szMm%Kk|r`
z%Tm?xNi?VMOx3&I#7}QTLT{oIL^T)J72bH%XmM4mPK1Bu{khnsdh-x=`5ZJA_986s
z8g^3^e~Y-{an5^%Xw$NsvYW-+T`yNng{I0;3W2-G`#g-l;zM)wV6*h-LF0D9nTi60
zzdJ!x<&nu7zl*omz-|r^V_0@g@*9Q(D-SXOna+mt>x~sxHQ^w3He@22IwR}KSw_JN
zsxp}eeWWSR@=lclba}!!>viZOZQKn1b;4nZ{(2<uV`MPXf0e55^Tc)|#Ky{apEMa(
zReCOgqH%#Y(4LB;MoeAnT%&hY*L)8zLiONh;qh(6CtJyQ4jI35HkNM;Lw^bSwi+x>
zn2R#}E8t319moy{m0QOi7x0jlm1GJXYBRB-vPK2aT1e|PR;x-uI$1PX48^7x$#k>&
z3^b}CzFCB4Au8Quk0i^|PLk8JtRCb)8=+q+XyFigeghoS1)|Yl=n;dC-3cx0k$@|#
zkv$FhRBx^Pp7^%-AvJFyizCLiW|6H8Lr-t=Rf90i;vF=?Y0|awYTrWFb)ye^ptjCG
z0)Mn4#Z~wz6WQ~Pte1Dt>!~E_NJ3^+A79J4PvPNCMeAOMt`+cN2JxwzNG`!MAH2?v
zN0rH0E6qJyo>PrM)<e9upWRvsN}Gwys#b3v=T&{}5_s#Rza809?UZUl6rmKOO=b6H
z<BMB9*+HbH4P;L{J(>N!0q+Y*7mhAmZ?=h(RbwtJbPJr82M6G<vcty>MYplj%2$fE
z0ggtf;X;<P@PEtjcIvsmyT)#cTB?!L8CJ3<jqrgglsk}O)d?r?tc(AZVN><6Y^Z7p
z%2<69E=ctr<^Gb8^9<y>7{6maQj^R+DAqV;bgio4ee8Xg(f>kwl~GZ4GLik3XKMY{
zjEpMJxy+{u+ije4B`X^^qbd*N`Fwz7Qx%9h6T103lQSgyo<pBWcjyiy&<*+6Le;!>
zaj250fcQlD$O5EozHcX<$ztw(29|0n5<3%Q=K_>^W)#=!1?x$`YBnLS$#_9`*tZyB
zp%?IbW`YuZNqj8yKbELh{UQ5_Hq{T%Psc$S-faYOs@j8VUz5$9PI)j^m9wfHv%zgU
zk!$su+%q1e@~-F5PW9}!DhJO)IZ5M1q)Vv2q;dnAQ<c^8ppf#_H<2%CSyg7mnXFSc
zl-R|yGoWvPXEyN4UF2typzC2?^_5mWjJ-=>zYcOtM9;<{+t=I?UUi&EZk73x-8{x!
z>7go@Nq9>g@R1@fr%%>)1|ITaBs>RyNfppTXj+Sw?;%IBn`v-E6kP`yTqT2ZPU2j9
z$Ug}Iln1gA`p71FjJ`XG?v!7%-tgQJY~gftrSR59C~1f4*G)Buu%UP8Cvh78)$*ER
z(W@U2lZuvc&{7DajkBxn%N2YRG^3T(^58o0tlp&tK~>jpf;A{xEncze8u#F*ik|ne
znjOR(9}#=sK+|qF_G=@Ri^rVuE+v(!A{h??t~wI=Mdv^{Yd{Y=i0-bkOGlu?8K@)8
zDrt^1`n8-`Qjv<#hdQVkjy%PBtvdLzdOj-$Jv@r*o#%Q&pj3NGCkT8{ght5Z>Pz|E
zl0Ru$MQa7%0K%QNLCpfLcmoYG3n?#yV}#^YnJ6<Ak4{o86iU>bfxHQu>VjJ=G`AA1
zS&R2^0C^l?&pM%8Cz7xpEuocXkp4U5MVg^dKC&nrP&k7!Xt#)0C-_rA0<OVz3C4b?
zB2!nIW01gMR;!aqp3?_JDV=y^%87KqC0)j2=%CtC_0H;aXhL(U*2hErbQ~5|-6_K7
z)lYC19W)<nJj{7iU#)uj-JlY}f}F_NP9hYc+8K^KD7qC2?PXW<urnD%&P!ckL=vjr
z5gyczG>euwNb?nFKjNH*Hv>M4MW}+BqF?vuLq4fh8Ldia)n*Ri1J6Q>E5|6zKwgnm
zm7<!KR)h6uH_pQe!u9jbx$BVZCeD@5J*Z|=-C^NKiXxe@I(Bu=WO;XCqi?WM9+tYc
zwxDBU45k;$y;-_h&xU-8e^k32;9ivLP!*734_QBXz=}>Apyo2)MsQa*@8~={d(K4o
zDQH#25ZxTou&Q{K4|J2eR!zi2bi+z)v@Eo|9uNDj%AGd8cSEnFKsNe2k9(JH5qhAy
z&~S8wvTB*w)LZa?2To~21GF1?Z-I*P=(j_YJQJI$D@OL)j^w`yHI=K1#;cK!D$A)l
z^|$eYdZ3=-EA?0%$8MfL^VaZ=tURSKchx(q3a=8hPqog9)f<q5vvgQRBA2SGl%_oi
z9cmsv)UPp&f7)d{S9z?%*kje32rZ0;w(7IfPRua%bV9F%Gacm`4~S=UBw<Sr!w(j+
zC+SorrQ-EGy!SIm#!ci>diV(PuefOod)kNQ*n$Tvuc@DPAz!-Dv#R?`g=Vs^(yHz7
z`~dqNjTTgGkP}&Taen1oTiC-WqDsXgs@3c^axutW2G~tI+EE>K7RU6%Nfmfe>J?H&
z-9mJsY}}VvyF9SK9yrYE3R*$DEE$dEPSc1@h5I$~=tAZaOzu*Vj4E4hv%})AQ}`&;
zvB$f?^jsi|lc>cMwl{>N3A0l@h5GT{K;No1L|Nl`yaQDbABLLE;NeARK6!!CcG1vx
zIj=Z|rhN&^soD_fPGwvI?6&+^VF!vytjcCxrIN>fR?Ck(&+Z6+ALO?Jcz+sIdy?I0
z(60o_Y=g^Hhg$@FGtgejs@%n<E`moQc_;$9hVP(fwN}}-^X})c+oG>k0oxBQpn94r
z_9@3;4BL%8Q<ZHpJzI+2RFx$hOa7Aj=jOBf_pqM_*uSg1T68YP_c+KpRS75;rAScr
zwH&%9fxtJRYgJ94So0M6vjZtu$yy7w5pI-X{%ghi1=RXnz!z;bUTX!~vzN20W+?)W
z6B=}#v#QUm!1zqPR5u7ScB21PwUo?v<r-v>2guG$A%dyKA376=HnmgQWj$2sv79;p
zb>l{|(`D!}VP?bZw7P=Yp@i_dqsW@7ywZqlV(^Fak+zl4vI6W#8OEh(LPb#X(0NM@
z-L1+%)zP2e9MZKnIeQ|KDs-V6S|*dhi~yN<2|iN2&pIfq%Cqsv?s|GXgg~^hu9@qp
z8edetiEpgvWD$KuvZZIxv~j#LmzBcr&!GR615-ZJ2Q`nd$1B-s<&=_r>f$_XIPWfX
zB&yPykJlq#sC6{zkUWoI_g10{uNnXRra=LfC7XkeQ_e-YZV8fnkG(CYX0MF3UHI+i
zv9Yqxx6#x+NL(8+MF9_azsjcxq3vL`?xzE5+{;dP7-S%t^T=DyMIT4PE2^g7g3Rmx
zci<-FfS%*E%DBsSt3yWB_wW?82ias#KEgW{)>F)GszTX{P{Mh=>T*?kP)Wu_l|ia?
zpT?f5`chqbs>@O@?jnO6DR)ejG-uOu%7LkRE*`7wL*hHp;HtGfLEPGn-&E~UR(vV4
zcoe!RJC=t}y#T~ql|l2MkTR&B5&62vp}1L<4`-xd<<*^{ehJl_D7I5S?p|}wSg-mj
z>fq6rxXul-k!RTf^(|WYHAMpIN|kRSjU73fK|KiWE}%=@VDf-et9KPTs5iF9cbw=-
zy=c-Xy`CtZk<}Ddr%LK}Bu}W_5+3*Qkmb=>)D+?##(pcFQKiv)$g?n<51_3o_R5TP
zwXk_vfLN}XZ8WfI`9fIRBqTHiB)<h7$uRWD;@)>6|MG)ZVq@yj{^V+@S*Q-Gjv~xF
zp5!-G7pUScg(zYj?_1vA9W=%+YKYZ!os3*4!n8Vttp2-1Y~5tu<#xQsLum90Zxs|)
zH@xuJHrA@LL7qzy@*<x{*!ng&H<sE1QBxjMCj6qFqcEhg9H|p3T7aHbCDBIk7kR?n
z{GM+zGAH1uy{y`1G_|bcVb;rsJw(lOF1s%3w&4v6S(3jcG`F93lFs{A92$eBoPoTl
z&QrZZicwV6uG-d0=qVd39H|m&^$_JL2atf(QJj3r#HMo$np<S3syw4)NWN`0TCkrz
zmiMfj?=p~vlTbz7eQo@9-9(>fh*!pRtq<*Shr5#YU&5L>M4?-tSTeMK0i5<@;*fXo
zvJRM-NgQWsku*?S<&P>rK9;aDA*h^)kiBzAjwTTFp*fvUtOmXh=bpsZ>+whGz}?f3
zWAzPe#=a%P0e#Rnp1-2mw?-sQ700sEwMf)>tYs5By41w8(yv3TlfNaT<|w@4LffhG
zR`t)KW-aI2$$h8-P>9eh@`;J8mrtwSiQ_y9(4fk&sWMj{(*r22p0EI`li+|dP!jF>
zHD0f(U-e0<-#}<a1{(RWi7_IGJ8I$kSP(d^RklNshw`(1yaQ!0<Ua~wxMk$dX)qgg
zj;rpi!Q?NLK@nmj>fd5FRfE#z|Jcx|4;^LIe5hZn*?3LD>S~exWT+7h?JKC^iwnpX
zQvHoAxN3ZFVKsXUb}Jox#bABk3e{67-G=ll>!vgG@eX3lZ^}WcBSFaiGXdqpVzF2E
z$zt`Bt^5L%_HDX~QsDfF@Jkjp=MFp>=T^V0q-Q6dQVs~<Cq$;7IOf6us{Ob|7WXRe
zQd!~6&|Mw9nqQF#1-1L?4a-Eb%AxQBwCr~FF4;GYcv5(w;!0J_N=K_Mdu%pE6Y{u`
z_~AAZCplGrgDUr|C|7-K>L(BaEQ~M@X-`Jt8sUz6@LeSmEJRr7cQu+yHPnmns3e)9
zz(I6$9Qs@w)o65!<Wc*R4)+}9x}8wHm6fNcjS}804KN6A-DdX`$!OX{EVLI^XZ5+y
z@vGzZ3KF17YSjhZ21{2>r@9;DHA|jWphZ70$X*y{Xu*~(Hq@*lV<wa~*-*HeeQ)7f
zLTJ@1T*%t{yhkCZ$KVqA+{*9y_&Ws0yz5u}NCVbV`LO$-K&i;8>}Dboqu%ObG|pV~
zdoGgzmQxA$wDafWHmk`~tGE6%vZ?t1$|@DJlY8L4={`*>QKx~tszScEva9k@WJlwW
z8fE2$kflLQ^;XLUcX9_E@RMqOl%-Ofq1=kF*!ytkEd1MA>JN5M$MYp!!Jm*pJptV(
zL(6UGQ!84LuIj~8-c$MoHuw{=wDQb_8Xd<DE3O@;if;uHn2Zg(N}em9NV=0VCGidv
zQwW!zLByO1*M%F7Q&r|=pK1>VjGairLlBBN!hUXHXVj-HPgjx9h<6`R{xQyy09yTk
zTwgsJOStX<<VihYb2&6y!{>bq`i7yyg|;-ZC+ZAd!O9fSB;|=^N0pC1<xL`UcGY-6
z^42@J`dqHoiBGEzHeFfhl<Lkz&7EYol-ty9rb1hF%eHYR>#>BVjOK}jTUBlQIkZ-d
z-F+&!P9wq6$g0|1Y49#pHK;pL9WAOesO7N-TB?sIjys-$zSG^@=N%q1v{%oN<x>eK
z4&%@yibX^tnuem-UH0QP5+dyHB6J_&4vNvo>M}fOsGfl2sR~@H`;eI&(3m*)XT}2%
zPbzwqrVS@@N`mrb@SQpUg!trm!&%!6#qWX##hTiL9QLgXDINd|%VAw5JZcZzLq4w<
zZ9N4WTa1NOZ^KpemXI$Obl1vSpSsso@28m|iiT7tkj%bJp!TEzDpYZ8MJ#jtAL4_~
z!((kjvmV3Jwm`W~Z#3FkvZ<(2vo(gmhH^qJU_IN2X}jS7>2{&cPISfrezU5*JXqwo
zKp(nOHOI<X7DHE6R%k9k9GXiUcaFPGMvvd2uScB<qGmrbT#SEviEbj*$yqaXcJh8&
z@hzwtG*na%w4zvfPAO<SX~Jsr5XY%Yd6TaR3u<-t$_o;9@wut3v|?A~x|GS>!ky0m
z)7^-4sFzRGgX&hBhNK-fS62qD1KFxTs|okNi<c`DQd2HeMRXUdq1P+BlfXLlHOoTv
z;pb05Cu+8VI8FZVBsAX{^shSA@(eBZ_}&jj25O1*SE6sN3N?Af(~0LNfjz0RRCylN
zTFzh752BjrKZndI7M9hHK#G++ISWVKHG1hPvh-I}$*!TQAdXma1i4jRa6felt?ajE
z<j5;hkIH+#KPU^WQ9qcvd26UWv7+E*{71T1SXV7b!Ux1l>Wh~*)5k-V1gd@(8lyQi
z$_T|6eXXk8b5MN&+PK4@c^^?Zr1zkD4%Hg#%E{bGB3|7HQO;bbF0I{%=5&LdNqW^m
zA!Jf9(=j---FK8&N0`5AKZ|)BGto&G^&sk?-HA<5#q}|0r2J(Y9@Ud{Fg(wF$iAu<
zNb$UQ#OfDXf!{I-O(wh}g>!F!BJx<(C9QsUWwkUnM?IbI!S#EH5`+t?e)kHpscHsw
z1689x2awKP_;#utc%S;ZB~;~TmXI*07S_t=+yspdk`a)9ps4t^_dbYB0cb+V*MTlj
z-Gp@9A^0&08z4I^B&-a7$BV>WAfB8`&E`t%w7T~*@PR9lYt?;aaL4LbjPY%yXY6AT
z`UB8e9-87U)#=K+)yxugS_(US*YCmxzHerdJc-{~ZZI0vPxYf^<)N$BP<ikepX!Fb
zOXW}jb&yYES&xF|DO;$HEA@+R@hVfKSx2$-fn1}y=o%GV#q4@F6(E|4p!%JD)+?tY
zzhMux9GkJlpJ3x=5^1W}<sQ{#$_PIPC9OIEQB$&6;2uww^%dp{sou*6l3a&}q`3uC
zkULAh&Js0tVwnz;Cs2)W50UFEbi8scWza#oHkZHe;`20k-(gBm7`v{1>mx|ebH<|)
z+NS7EXH+J22LE}Unz`$|vqF3r)vYMsEo@D_c$yir4jKvHpuQP|O}%;_VEwYNs><9d
z>l$XV(we5AxUS8$*we>sE=imEjBX>h;?_m(v-pyy&=~4~+Q|J5;$6-+weEGqD1~TX
z#aZ$a)X6x)U1(0nZDc(WpErX(o%h-Emxv-~!NuabxrRTL)0J$hw=If0yG8}Gjarju
z@qRARt>Z%4a@YZ(QU{^tPN><AeYouFqlRIk;{p@k`td{dVt<4_S8={nc3Oz!T=c^T
zdT9}}MzTwD!24Cf9?v~2;=R|yOWHB5Rt77TwW>teOKs&F_+bm6dJo8ea+u1#Wk69?
z)`UYbqFuDqE_O!IZJUX8y13UXSTOmumTz?zU9rNukUHf~?q8BCTP`}P?y!L!5^AD{
zGL`g*Pz8H|3dm(dDq7LWPAVExZ~PXhE9}3;_yo^GqY6+{)o`V;YrR}o6_a{}?4}|S
zAv#&!UxuDobI8<IyU^pCun<<RD~jD#E=uUKI!`NoC)PxlMWMl6pu6gzuQ8PgwP<h6
zAyX~Gdj?yc7f2_ksVO7Vc;AC?tn5aFk)VEdv4Q*56rEFWYr3P0SzhUQKrzIC>(EUd
z>`^`XK=N^f-BBOhWv-E2I@eiEEoVJ7&-+1grtpq5M@)5+d90OoRNYA*aYqr8u#RE1
zeNd>`pgY$_X{7mE%66ZJ7R6||ZJ;)?t*YQv#inega+AtwD}E6we;6;kkL#SrYO3z{
zB6eMMW4hjbD5I>YvXsh(M|cBhuMjpw^M#aG7RObwkJ?vh9QjtF9+e?vkUs`#%S4uy
zTTsu56+52ezruQ5@QSb;P4p>d!o)GEeriAxlv!10cN)KCv9IDNVa~nC=152r^rFiO
z*K|NnRZ}{V2i5jThl!sROXhetuSqPMP<Fy}-+n8U4r+Uy=vMU!**^73T{Y432gJM2
zvgbbZrm#^-g?RWZKa~f#&n|XC&k^1ewG`BW4N@=pB(&TJZ)@cd_C*yvQ@Bf2KdD32
z>3bSYn~t<Kkn<iuAFDc6*oeBIClaqKD=Gh~(030%R&~5rkjy?4aTLH)sw0uzkmu56
z^u;mnI%t_z`U%m=HKyD&ppj>BRmG9%c%Yg~lFM<ESkQ+hD?uwKgKSBH#SQYMgwLqv
zdON>Kt8O8h7~;KMF<he>BGE}`)**IK6@!vW)jCyq3$SnUHA84R$z~W5qY6FM4++at
z52<FXO!B@!wbcbO(d*Hn-B>zlLuEr-`Azd;G-E?Kg_TU?+ekj;XVj@zkc%$z%TG|W
znL;$8`i*qnzc{bnbDUL^Ligzw>xHWn1F7>toqp%IQ`Ih>;+5rS+KJQ}EWrP1hn}j;
zQdRFpEVpVM<s({%y%5C$EanB~?CnI}<nwBsYD3i7rn+@;th#rm`(8ncClc*Fi3ADR
zoMOB?Nrb9E<h6)vvY_W_=ljTT2KVqhQEnwsuJT{1d|A(G)jKO+qFUXPbO6OV&)5e_
z21<Tm?{_RGHd<$VRaqukm|Xl=VE~$8d6BwjRTU=V320u1uBr2>3Q6&dZ1_$h?q%3P
zVF<gRg{DAhW{2YO4)#(rWsdNy?1HNDRLd?O(!y`mjoX8Ksz%GINs`|0_dEl>HvkH*
zc`r>SQkl-)3&mCJlmVYLfuX3bBNqu%@58s4h4?aZES;^2f|Kw-7ArU6L#eiLhDXyH
zkHaDF<JEo)#`hk(d=I2po!k|0qb68I8UN-Bs8dJ;af5e#llP>`3h{)hkJT@$xaSbm
zCznn=wB$GlE!8VhZET9_e-e1#>XIymKB{5bK~?7}Jf(PpD9NK&1&`Hhe$!}A)ze8g
zzfJ6_{-f8xXVsH=gEI>W5EdsLs(n^nNk37(Y-R^GURjpeXjM&skl&{}7W%4sKGo3X
z!TlScy_1MEiR)MyvTELsdU02BU1{A|_D+%4T&S7H{fQfI^AO5%2#P7IrOH4>$19-b
z5+rjD_DWR=!t#U=4M9)&nW~>%4;RTBP;JLXvY6FOf%y{3svaf9WLgrDsS&QH2?4^^
z7h?s#q$j8<6h&u(9Zf4*UZ5^SQDruKCWPPsC}19aU@g@B9Ao-Sz44gan0+t^g*9VS
z^Ak654b4YU&$%Y0DZ1CZM|JV2YI2Z0J<I#(M5}Zn&D1tB`>6rSt0PRfrt*Rl@lCAC
z#EGny28y9-NvK~ek<T7u2fN@d?X>bl+gM$|cUkBx-qYvEPXrOiR<wgAMp?+Kx=>VM
zIuXB4vY)|y%*I=KgSD~PP0^_dN~s@RzRyIX!8Kz-wg2k3Qx?veNi_qS9O0p8P8CRt
zxlh&L%EQlo*wG^|NIqX5??CT&h0(OCqmgW`pyOFNuS_KLo{8Qokp$Ix)S(wu5fTee
zNGpK+5}{0ilBsw?Q?Q2R<V@noK5M@100{jaXs;L`6+c<i4=3RZIhf+Kn6+)l^v~!d
z(}b`h5ZOLx{~WgaLnt|kX^HdKgx35uKJz=A>m53u6s-<(Pd(UL%{m(L1;Kr4h{^(J
z*W*w$fWA_V`YyOyv!9Ofn2X$BbO{kxMVBfsl!sE^vu0>rgu3bxN`*SAjeE{phDZB~
z@s{K#cX1}g5UOQY7F!eW7PEWumphR63;0$ukvVn0_Yif=<Lt7uiV#&@p$Q4}f`Ab}
zPoG!<Xa0hIlpf+9AqVnkC7r6pxJ$%z$UBP+K`U{orfF(cN;T9|tWys+$zrW%&(+Nu
z1ucb;s+v%J`|`4dh->Qbb6{Vxgqj(m+QMk`WH!D`jPEM%QnO@Q@!s>%pXwS>6{T|7
znky5DJSt<Nj-Ubae<AY2)DsLtOX2h}ROc$E`2y9G!zM?f`G@mBo7Bx)0No=&AA|xg
z!{ZX=uOo@7ZAd`ZsxSXt&Sd4)_MtC&@MtwlB^u2k4o-p^il7f;dj{ZoVJO3B-UQ><
zX=X~mJ%sF75+P|>2OTdWNt&sqxupfHRsUumni&oS2Z(`7)Q_rpUZSF`n(`dVO<Og_
zvdJ0HJp)~PAM{>RI5ojy42!cWFg2~~EC@*s9kzqSM>b;KPwh{m@7oNzBu%fH8r3Qs
z#@{(cUZx$ar~@?ld9*_kF`zPn#ZXi8h8Lk16ai%+g$Ln;UEq2`d(#b`s%*8WITx>Q
zCGTE-rTit;ZYXoD9xHV!sfI`se#(*8%Sh{VQ%9*dO`Ri4jMmef;d@YB5xY8fBZ)yz
z;&%>>BDU%mQoX!tgt}R)s@wxq-WEXPOkSyO81-GR@aCKA6k(4S!Xb0e-l|4Z-qj#)
zylRLDMg6Ag*s1~XJP8FgZB4Ug^vnPfXmx#*p?f3U>eo}fow6V`ynoGU&_k0dRLhnJ
zHHFZtM_ly?l~A(?IuD@3rCZM+qtf6)#^mM5gQTMdovdk%J3)ayLg%+5o61Q_k5^+m
z<%=Z1@yIcKFhcyR@s*z^ni;}dZsbgwizF{u`J=^1$QEW(l%r?k$=04gkDMQ^K2YD0
z&ZdqfMPWHmM{___w;*m74&H}!DkBz#hEbP;y3!SUZ$UQqU>{X2ugt&DHq|F+j&Ke7
zww%2W!2^m$;-F0~l9X$_;C=jGeFf?nu^yU(sx0>$<gF3QH;CoboP$PWx(=GF{xTPx
z`xbw7Q=wo@2fu(<CCSWUuHKhqO5?b~ZlrY~meYoet3Fp#5;Oy-2>Pg-OI2K&+Fit7
z${Z?A)9k=&$e}VJnmD0OY1N6`fS;`jODij)F7KVZn+kTQ2g*mYstzrz&UlNO*PxO-
zd38?7UlxL)nIxJ%tiD4Fb#Di4(DX`Ge8?)x0tgw)MtAms=P4g&RYzoU<};j2QOp_R
zTOEcAe9Vne=1g8tt3iu2yHfQ^W%yL;S`seD?0IOGL(JepIyK!<{c!RTMvMe1Hojx<
zDB)&Tz1rPGq^F8kFc$}%s0j(tJg0|df(r}R4BL8Qe9ahshkb5?(yBm@W8YPE`vT8h
zM}}0jx6glw8N>2tbEs-ouX`?Z36Ptf31tV+;!`<PvHOAl9CN!W@PPgrtl>TtBkw>(
z&G^&Y<48Q`I;eh}T4{S=Gd_@J21&af<<NIo!T>ujUf0~Cn^?!q#Q3N2Jgu2)n$)eD
zJmH?XXimi_^6)#+v{yW;ven!lWu+e=n~Hm6|Ku%OxnX4@)k#-uFxQ<#Cw*AuxB00V
z7L91f7JSh2oL$~o3w;-wAYV*mTFv7VIzF@U#)KP5)1KiTbQR6KJmlTLyJ=^24wP@E
zH%fVE#o_h54^_@9Yny{!2!R%=GO&cYicD<65m0kY3sM!Sx+L?k+`<oJ-#4Mz6njPB
zV@*V#3Psmce08m;4@G#QRkaWaV)QOLu^jJ6Ipr{{-(`cLs(#w)_f^+!G_}pc$gL(V
z=%G2(Atb&D8LmUtE0J~8ZmOQ_CbaVqv9|kau(QvQA5(p1IW%0zITEQ<P;HCqsI5vL
zO*#{Ys$XOl)+LVLcd`G<dMIzuMtq>2FlC~~L-9D)X$NmWFHNyjUC1Q%LE5I@(~lL@
zEPqWr6LzgGdDYnEn)53kDX-%gK4k~K?@lP$2^w^T_@|$&rMlB}W%XpK;!pIr>Q#n1
zhcz=DhwZ!R<(!L;qbe?WSgQI{C64g)5#+rOANfAkSUJ#XNP0C?QPxskfu!(0akTK0
zbx386sfSP&ZYBP(A}r+$m1`Ed+lxIo1T}l0jHFu`k+E!qCL-c@KrO{Js`!({CeT|g
zEMAXeMpvIAikt(5>Y&|Q$Xg;dQuSdk(p_K9>s5t&3ty^%^&$MAeXJ~JO*6V(S*aBN
zFUg+_^S3JO)jd`T{|IH5uiS;6odf5nZuGieod?b6#eAb}Rr#vg9eF{ro8nLPJ56Sy
zX$QV!F{qj~#X>gIs;||=Ur|%}=iBfoNEbBfWSwUDw?R?myk;R4s`XJVhtM8nNp72!
z>OLDm?ybD#W_Vy8)KS!-XkXc^vq+7yb_>YGRh#<lPw;A0%cpb6e^aC{KT#71qp9k9
zlKQ`=;gA^Yh-MiHCljimnvs>zNF9Iked@H^1_mgbP{2K@&p}yI)msW}D8t68<M)vF
zI=ibGnwl(H%o<6ex<WPEUzpG)_SeGYG~r2gcDd|wIr}bDPZfJjXjsk5R4uh;D>Mc+
zulXhDuM31Hi=W$&=_UMcnEw@^lhrGvIm9nuNo{Lu$~3L5A{5Wmwn_ATs#9WkbZX7d
z(9+hFy}^JF{if=zAl6GAMw$mB??pCA)@?3!{~$<ciYpbZDCvA@l=`B2IgdJ|ZzFr5
zm{>BO<5f*U5Iw7II{C7Tj8rN*k#w$OM;9Br+2^%1@;%mClfhKYsfqv6%ZhyEH6^gi
zbnW5QXl|owS=1NY2e;IjSjp=Bw0h?wkataUyJWI#@%TZS1lEP+6)v8Ko)-!%ghxGL
z4ZJti4k$A%tWt==DI`P{-imnDO{NZ{E%-L?ve&neL;Y9zbJ21cx>tJn6o1W!Zi-#9
zm?p8Fh(moqMQBg;jp`E{ZgI#9m5$Z-Zoj9jrtI>XKO}<G%#o*vH<gLLj7O}V-FJD&
zUyjDgZd{vQzG!VLQE3HH=Y8yECXzapI-qcaDog@X_!8UvLf{K@mi0L_nlh+*4@-_z
zMX%hY@XazJl<_Xj5mGgtKCxyRJEp2y`N9|3{ar|LD(`U__nzu~pUMPHITbb4#b;^S
zc}6-9v!C)ygft|RWou!^!e+dtX8a+2>ViI9qX}78FRPUoaluj9K2=8ioX!eO70Tmx
z^;Xmv{UHg~1Pf)l)+5DrWI@kSaiWR!KH^g0K|;`bIivat5A(PBQ?5c6Vq@N$Vlu5#
zKXC)JktZikvMjD<Y+2rpJXhgj5%{{A4O9y~w}3OKI#+S|7LJAZo0_E9g*Q5cG~eOc
zs$h^$d<czP$Zx81zU7llrb3@%WG~N9W!WhkMU<D!N=*`<gZwE1D)%?9nLxg}k9j6}
zX6lS&Ms+Netcf!dsH*U#e#R6(Rrd%L(C1yW@t)Lk*T_R%lGc1deI~|a_kMJjYO9np
zEJwE%u!HIxQeICJNM~R-&vAug?1Zv1nm6=-Xk!stT9uQ^aLBtAexr(}R3itw$$-8M
z{l*jThn41l^QjY9lc=RbH5FPluNlN0)zB~j{C78xYs?I~@8~o+7xf{kel`=G;^sc?
zgLB-2FJ>}n?-%wjn6%eqXs)QVfQkb3fhvEkDl_H4&f*oQyJ8XAK)#oJR$&Km_zJmP
z*~(s3VwZ&Dsp{#fLH;Vxw%JIXkL->-Yjtu93(}_tscQBlapXR5zzV1x%f6^5s2F<5
zBNOU0-*<uNWFvH6!(8dbOm{j3XUTTU3Ts}@Bvw@-#ZCV8tf;1gVHFae1j-=3kWCOh
z<zeO03(-QHQy<i)p_Ji$yaBc*i?{*me6{bYs96c0PA1;0GPNeBp_3*A$*a{TE@_HE
zBy?2{LcPfHX|0TiVn;>mi=eqY-*rgEHTF@Mt}u6fa>bYSpF>SucN%`E=G<6SCyHN%
z2B~vU7{)mwQ`HoHK}LJ7;VRvQDhs5$<*BaYA?!Q{pSA~`oI~Y<`WBBu6U`(pGkR9q
z-ug?}kEYul<lQODpm<rarB#iT1qLIlsZMNV%0l?ALVo0vYVy!!sG|s86*`)Fn*|3b
zFT4c{Ya>!fU~j)*_ciNOvrY_7gLPH!pt{a9^VY@gw*+&rg@p#$6FRGviC{WI*kYZ%
z*?2`G%zK}|CL3Ju56N4vr&48~(brR<^g@F=1+cE_Z_zZg2=+&_-1V7Xn(d?x9>pk{
z{wZqSgnLy(t15NPmsKZ_A_7%KXabBrl|&QNeOS&@=$3PMAF^1(cp!?+tr%q^5-$|Z
z%DKq5i1A%9`5S$LgQn=|bB~fi9R~1Kh1Q+n8se3k?kw*DthOc@srFZM-&7Z*I$mK*
z>Rr;z;3Dp0J9X6iv9Ic<`J8H(vC5?ksI~|#F@TRNTrvYn$};aVuU^D6(z9E*E6qGQ
zh=mn4Va+4r9ijg;5lo2t9=uSYJeQGIO-mDIC{#~*ync_msJF4R>0s(pOl~<JdN$)N
zEyFKX$H7^+v>%Lr6|{UGJhzW5=P@i_5}iY;G}44)bsK8R>0EYnlW+5y%(B0QmVX7^
zz2rBaK`(34e~R&i)CD7JEKe)H<YdV_bXPU@Q<djEcp>V5eSp?dwTyD_Gw9^LRH|uf
zKS#q=8Ct1AQ608*tX9OL%%ZA8qI`3`lA^gthN`UdjejgOPcseFJFkjV)j`P0b`fh_
zK|0f@?o5UTG3bX4?6)G_YLJ@(=wwmTMOH?2pPz6Cso3u@Y{OA}7kNuvaGnq#&9z-g
zHCi(G{QG3tp60yj=@fpd3RlIis>f7aXF5_R?`!}bNaS4{HutaE5GNE8rIm${d}><f
zMD9k@?IfvH+_$2Nd>#}0QNAN^o!*_YZAZwJ&quB_XG`6E19-9dM2hOsZomd9no<Xw
zVxMR5dzA^PM~aVN^SX&_^eG3Y&~WM-R{y_hdMdmxo6P2iAoO3rAp>wyp`p5>Q+4ew
zp(a9;QWE{O#kIxF#Y4rfm;9>au=%u;FuYsMEfvC6#x5vdzKc4WGW?3&Onq%cH+Psk
zntB-Wk%=bak8m`Ve7g}i>JzB0NK(jaHi(R-J;@>~AFB7Q+R=kplM3&jF$<*{{zyaa
z3r!D&D)BT?)9Otd=Il$6ClCU1<*GA^MK+&ABWs?9CjaV_N;H!|_>Ve6Z@M-8X%$%h
zr({8$qci*TnOIa$Kp}OdDesrgGnyUN%2oAwE)QJW@j!$CXy+@iVcPQ}?7Hgs)Z42r
z2}$)4Bvsa`4$o7Q^krogNy&O@N|-WF`c#G`NKF3ds#YZPoQV#^q4Fp4JVNdFnQ54+
zY*5$ybgZz@3L(PkgFOuI38fb9BaQkL>+Iyfnu&0w*(RXX-{F7zK&j+UQyWw=jmMml
z4J9|2#FT(kDdRNFP}7o4)k$=L?P;1;BlJ`?Ne6aT{=cFb^@j+{R+VWd_p!noj<@nG
zlDq&9RUYtcql;zT^iZ#Xss`lwH<JIZH=Lmf{+dp$4jAQVX0YRf>|h-4T@`<;sAiB&
zXyCf~#D!`2GxDXa-g$ivhBa9z0hIU+P=+!%>TS?m#lq@Ek$0}Vt|m<i0am4}at8ZN
z+@tw<Q=qSEXH=ya#=dGEqw0iI#i!Xd%FM3gj`W#Ks;W}Q-w?J?lS7W-lc)zzGqy&M
zPHTpgC@LIRNbXd;qsvI`GFB^R90Qi4&+0hI1poC|&~&DCDT)-5blu=Vo4^r^jsK*4
z%m*Mp`fQaih)^~ATgdNjGT0yDA1w^1{;OHmx#V*3f#S_bXlwBwlyoDX_0(xtGbY20
z1Z5dCGnWZn365Xb?e-UVtYeB-3sFKgR!$zS&;vey5&t_Ie@ipl)YYKO<!i{2l}FRm
z7uiiqgK2_EnD2S;fLLsyHB%~$s8^rwcbKyYwU~hp6aCf6ei}QN1D!s?zPB3tt_<f`
zPF7P*tZLcCpv0;kQsqGyyrF)1VeQIP3zbt8Y1KW=WumE~8TD0NGPUliBbC3YN#L5I
zr#VAWM3kzKR6V*rb6_4+RbQU4T4g+=(W8nv(y)rk2oCX@=}5QoXNAVQRbS^m<WN~i
zp$v+0RZE?Ol&RZbBR-Qpk1-s%`v7?_!rLn6KZRIiVT)0|<3u|`k9MP#W1;6N5Q<zV
zdV#;yqmoN3Bx<4yN?tGdDb*RdR7|QPU<Xuqjop<usLX;j_s@lm(JZz>);vqhm(C-D
zoRc~Y^vO{3(2S|(-Q0$kukhbN>NL}_b!AADK1=Q;cJ=`lQJ?d527RTP);{k+GkGe3
ztM+5%MEem_@uRN1^~k0=JBRU!=xT*S%DASBl{@j3nu)07d#QR-)<|6+%A6kbsG}v8
z`KQw1Va9GMz7%il;(Y7S^1bK<<+;>}pIw^DnNK@6Ae-`W^`ARfEXlm6CR+AUV_ASM
zSFh`M?nlz8PQDP2+t6Q=71TkXPcl-ay|_oU5M`XPf^~``4!{@c>{axkiZCDAcqj6&
z*--h=RI?N0FDDvm%6r~JM5C-n1)U2XDw3EBhzF|~#KMy_mH#O$sp<o7;s;GY`ZN81
z!)(|BA~WR~SMkVWzYdzH?;d!mW;)H~nltFfecwz)$Y5PMQW{RZ)w8VcC;C^sf1G{H
z;yo%O7E8rmhKchBv5V5C8PF=m_@EYAE2Q=yXVzzjm6=B#yrMJyoX_9Xr>m$A@HifA
zCU!FpRNzVGlvYC_O&bfxb5eD54N<S^T~+%!kNB_&uG&tma*C_U*=3L9UHO>Ir|8L=
zs9A#6j)R+2CDTHLa**1ab|^2)U5;(mXFINEt++__>v~VB>s4RDE^0*!v9M|Q>e3ZL
z5L6}EVmL|N6Vl~d;D+N|MV}O*>2sRu63gp_xN4SC77tZ?s^{_;ws|AkTl5kFrCLGF
zS)Gb3D?4@%eXp459Ch5w_)W8qF5_P*ldYaneR9zb=HuQ+u7%)h_M1LOY$kNt0S)z;
zjG7j!I-Pj-KF;u`X4p+(YTy_AUweKH`UK#B0F=rg0z5~povO|y$Ld7f?%9e?)Euc6
z>IAy*UJkQPvlWF7S$M4?kQS&Z^hX`Q&ztN>9@kz#<e|?Ev*u^07f5r0TIj+K($m|)
zoyvNi1_{bX(+XEorMqTDDSOt16o)ghGguk{jZ(1~H`se&TB<M9lu`B8yP%m+e@$Cc
zC!;>gXBM_u6{U(Wd(p?TlP(kM=dhohP;D28u4YkfW2Nfi70*;j6Tt<_)t=)Z*8{St
z+CY6auQH0tvuUQPlcNw#d%);rC-f*nm#%>NdhJ!@KMIQL(~_Fdz2T<9T6m%|YVSa+
zTAw=bG}X`gRKHE2>>pDh_SaxZhsc<{NaegM_$8l^x`R57D|r8^flYv3`V2=^Wo{q>
z$Ya;5k&=9u_FNNaT*%NiCO*GlWLgy?Ru_XZ7^(xE!skh;6D1pbEy9`Z*v-1#hMw}4
z$CD4Q!GAe}yxpNwshySj6qeVZR0b0N3Un7e_4zjHB~dn@fQ+*yltjCPrfZU1ADVl=
zsqo8Xo%(hkAenX0w~YNeh+e$JIr~OG3t1Ue*}VsxOY$pAn!(lOv86K+Quwsya0r`K
z-Q{9rMrfpJ8Pvx)$?KugS<~4b;Li#7wI=9k9&!R}rD>!+h2iA}4xmBRVSWv*E<K?7
zH(?8!Co&g$hoGsZA3qiJ;Ki!S+Q+WxIvuQ3&*%g+sSTYP4tjfm9OjEu+D~PUZ4}y8
zlOKfy$@1=lo_SzIdBm-XUe<sG>2t_082bFG<Y-Av$v$G0Mk*~&W5JRPe=Vjmf0_%%
zVJGh+y<N=N{%J`jk9{Tc(2lX}t+JA`*{V)e^q5Uu&I0Oc3+zoL9%9;U#&?y6tC=yH
zH=#;wP3+qbrWb3FxL=h>JAUFQM03mcjbvV#u;X?39kPr!=qP%P`hyk|8>@R;9W%lJ
zwi!yPgGQgomWrLL;Zq{?iRMiP$5So5Fdt2bQjM9q&U=~7GXg!-%ebHW*NlXI9;%5y
zP8@m>Us~~=(7qz5c?aotgEJMu<LjYt6?CXVTUrr^;*-U^+QvKseQLWt$6S*aN_bu~
zgBrM!q+b)rHIMHtW=R+M6@zNhu;wJ}LI(*e)+{+EQuqR1&LDoQi?}2*piV|*0fhPK
zA)G3Pe5X*LE;7}F;1xT9htY%E@KnFc%=DU)r6qsB(NwaCDC9LPQyaD;fyhfZ?c3l1
z=g3B;^2ssbj&)GfhJ5ZuPwSJ(qmc_$J8CLl8W^2w35w|&iL+;x3>7ae`7Y9_jstb#
z>ywey>wD0nPw1Qm3Uu1(VrqM0Nl&q-<P6qMpD?IcQ?radgQt#w6z@Wl9wnmC?C$$W
za0A+7Hv1Jx^rTN`5ue;(_P0KXuO1pjqi=;7syazi|K?Hs6-`W{&U@ugnE;Q6IhWaA
znhCuZ4Pc?I_0VD#w9s63p(o0?yhTJ8!R#7+f~>M<n(CxnimK{wm}puSP_(!}zh5WU
zSEfz#1e#cVkQ{n3v)$x_i=Okr+~4*0g3zgQR#|z?=TODR8Mkm^(e1XkoLE$!F#mBN
z3sim+h=%IkgkA~5)Fk+`WOL_0_2W=fQOvf$Zv2X0A!AERCX{>zitq<U-omkm@-7E>
zN6Jt%(W9eJb3BHX)1+Z}jEAwCgUrJ3Kr72*O=34R*QLa5eJYuHci%^jzl%R&M}if}
zs&2%ZCnA3|1K&oUi5d-V_y%&UnW~@K?^F4#nI)RZRm1LWU_Z8dl-bgh%Nt;T`s7sg
zZqG&=>XW6<@a`0C4)B-y2XgS+5{%5K%0*f1iO6O+*n#qKvYg6>-tj25Xut>EXwchq
zw5&dt>OPcOZ18|7Jfq&TNn}Qrfco@vR@FCIJz|<&qp2dQELA>E)t&b}JE*D4hxR+s
zQmQ9Zt(ZPrXvFxEnnA0&1)=i4qyjS+jBXDR%@!!A7)pJ|F>sjXZ01v)@Dw<EUN9a!
z?@fA-g!mS+R&ldxofY3EV`men+N>rg_GRD$B(WTwITZ=&DZW)awdAdm$?WOhL(T8t
zmmPNh2X;Wdf_gmli3}b1pwB|3U)cXG>(rsJ!1+7a!*BDRgz6r4|1o>5*)wlA|5wNV
zVgC_*WsATZ{v5ydC-^U-;18gdkQRLw`z1JGn(O~=@7&|6F7thllA$S)8IgmG89GL6
zLx;#5w<+T-^FoA|9COzCt^1WFBB?nJkz+(^jFA~bN~Vk&LuQQIlt>vP#?aKr%+$z?
z$jlrgBQx^6KfmSdy?gC*&g-1l-sk+Y$DfPade-mqe4p>_^SwL^$>P$!7Zb2?x+8zp
zXG&Zv?KsX`icD6o4B>*RchBWL<-aIG@CDpw^~88iurA*4q%kXv%wQVQQY%<?&mrzi
z`J$KMtCu~)p}%^C)%mo~$y&1W4lzRbLY=RaduoJ&pL&ObFEy|_OF7zn2~bmAxJUUS
zIe(q6!?c1_eN<}Vdas*T5Kc6V{rSGZ9=(xHU2eEs<F(7?T1E=DYQ*kUh4ZYyI?&mv
z_?TOulB(q+SwVAo_n@h*e&b6D6oMlti>UwHKjjRY(_&(HE6|(j<LBl+rHeP9{ojY4
z%0!x<eaW`M*4I{J_>IYiVQe*0ZRc(K?Z+G)&LsCVPX=qQpRd|mf&G!}Jmi?~7~;5W
zKZ9mXaM+<|w9}1$n&%qtp6HqBP4iXx7Ge8McC|SV^Y_!vMrW~eBs5BP=J5MuBI9S_
z<6Y1y0}k2l%X0Q(HEeZNLf=71(N6AqS6}G6(sRLc2AZ94uW+}zyLp%ThPmQhMXn|8
z9iCdQ6URN-z4N`*%y@yfn|Co)bd~PCygS<$^ai2C1<yd<Tk1Q_ahNBs;xhxixt_)5
z-y?k^`TJ3SY+$4RurG_b>G0+GOZ@r%JmzY<Z<)7{kr`Y`)NJ<c39cko-R|#Tq~<Z+
zKiR*=f0XfV-&pSvPl9(DD<YYZ=d$hD8apCTf*hXhYw;dqEw1q&;5xP5TJL(FpHFV|
zxBKHl7qe$(p9$_qhpq^8aNj%7vdz%5)Hl;x?#=TL4esE*X~B_T1a+avP(k+PU{heU
zKMlG!n6sVb+syaWKtlGR>{?z<+l|oDW?XN4Wc|oGz}RFwY+DQ>^n-Ie+_1nK$(jjz
z)9E(+K3sd3{SqAfq3t>_hquwLFXE>zg6B@5>C~s&jcf^ePh%y|buPvJ?dqD0racHI
z(Aj+oS}L;fG@SgI`v@9&l<)UYQ`yd^@DV0Ky>P~$IH&@C`wY2!8#Zql``B#9mY&O%
z?sxA8i&f^l@Pg$~PCJN&YkEh($A{PzL-@`qe3&nwtt#HM%R({Fp)QQ7bkY8S<?!9B
z=sDf*a3~@9tX(<PQRiN~>n7IJ5l<v`fKUKMo0njJw(;9+cvD@}g!<L7c5^wa;%+Y@
z&*_8%UN;R--UZhC9-mZRxEoL0tmHu-B*W<g_~a(?=hbmW-S-xo+&K!{U^w>G7*>cR
z&SO3o@i6Ta@H)D4es%^KCf9)GW(EvXI|SkM23GPn9=Gp>KnAwZ^Vyq>Va7hI$9ju3
z+L~;A40?8gRvGN5_N`00P&KLJm?2O++qcFC##y-c1zUgn9>yp~-Q(PjMWQ|s>hY|c
ztDtuj^7T7puuur)g?$F+h4)_#?_Z_7S@|8GdbLN<QmFYfUU?H7`!Ta3+fLre5omh?
z?Wmm0Lta@vs?U$buDO*@Dr2dX*;1~$;$}}Vr!m+OpCU(oOT79?_X%vGTi`hZ?PA~;
z%2)r`bHto?$-Giu2?)Sr=Bg07vk4iu3H}VA!~TZt+YS0mHF@(Rtf0Hd>bL`bOMw!~
zF$mWmrGVDRt`MgvUfqW}6ye}wL1%X1{Xb&TTvqBq_*yxO@A>NSE&B4_@6Fx9riLaE
zKe>izn5v}1_0rml?G`YFY2@6hPWxwbM$c&`e)cuTouM<_i8=-;?>^vfL{_x<{uC%6
z>ivh%0He;DV;u=4M_PYlZL>aY%(2abH%~i9xN|*0cvHP!r?cmR&FEtcFdj0#HC`a%
zGSEbLeh)2Aqc_D{vVVoS(7y*9w*Xp?!i!x?lr0@92|a%bd8qy1pLNH;A=+C(-Hp|;
zC?CG<3<c8BgF@2P8R9-H*citD3)GAx&ZU*`nu*iKU`48<jB1}gh0lhVlvl^(q10X_
zGr#-ct!Ft~5WWcGc;dMGptLMp;l+p0U@rt_vCid9Xpb#b^2v`g&`Q4|W~8+&d*@a7
zyC6FOr1W9DedXe+ioeJ_y81um+S>Ia5{v02){bpX6&~O-U~x&%QuTPs=6wjN%8OOT
zjegU5SeMhC*mKmJ384(~!B(Kbg-9G^J^coYxE#x`2Gma#3_Isq=B_4_M7!L4hiv{H
z52my5aPBk>nWKFm{uDSvr{XGWn)Q%nz9oZWnI+D8+S<uB(O&Bq?b>E`LR{p_VRi~V
zM_m)F7cG(2+13tgy-{f^vllo<q8*31)?l?%Lir4yD2J7D$XjU(*|yoz(XI`S?aq4F
z38bSP{WiqA1?=Dw7SlA(NzXDgax`2VjaE(Ln9aI4>q_^`_fGIFMXzm!{u4cW(T0Vd
zbKXjKrhBWa13g@7O7`>af$+n2IP{`-ruU4ekd?9v?S0I97!7IhsQx~W8t}15)B((p
z-=E{_=1t_v(x=B*FKxa-%waQ9t&7LvDTlgw*e6GsX}52SXO?>*l2%uc-kh0TLG*F5
zKg)NHs~+}E=8X0Jtw`Pi??g``YcS5Y-n5Q#&~A%^d*QyxepwMq(bzrF#smBp(dC!9
zzwO*lKG!)N>PPhBa_|VUIfVZ1V5GEnJU-J(&q~i`o@9+bDNqT{Ffx=CTEy4+*@eLq
ztig8f)Q_A$<{juwFe!i3pB|i<J;Ug*mcp40mQ|K=sMTmaWu%+btaP8@+F8B{%xkXa
zfGgTs&4^=`Ea-XGh_R1#xScax?XFqSvI|$2_DRDNDEE%AwIh$Wu_`W5Yq=d7^~T~2
zdNWwZ)1c;g_aJ2A5Nwz%Y$P`pQZmsQJN!Hw+FN*z30ze=VUTBudlsB*N26STwzJ$b
zTub>(v9pK_%w??KIz|?F_jr2aUG|0Y1>W(<tTf+5)?bM$6MO2I_q1;#Gr!$8+}j&&
zKf-!Z#db56>J~UC67MISD;)N1^>+$%LYijzay_x`@lbQNY2hBmwwalo9^B}U#y)9b
z9TY*$HM|b<yE-TjK?ca)O7XA9w%m`clII`E6KK`Q+M5BbYLN_kJk{KD6f0&#U>IDf
zRninV1qGM-qB%6Ok^4OFzJ%?U4Bh+U%@iQL3y|Iyf|lTFY}z%f20t?2?aMH&wFKs?
z-M=4esm6E&&H0hF7_Q5)?z0vdPulj_`?0Rop-Z@lGQ=l&-gLcVxz94h@+-?bmSRTj
zwp~N6RUTH!3G9~&yoXirrsVx{?`_ysHAWBH!)W05?Z>GBy5GcY>}W$_V^1;uWBiVC
zuX?DC#w(2@Crf+7Ny0D0dQlChtgi|a6;Su+uBL^bi>EZ0e7w2vad)h=*O|qsNP}Bh
zZ%s(*F346F&-)(Qv;w@|YD(6>!FSr%MLGM*7gdJZ>*&n{R{AXbj$6=^%J7v&tmo(o
za;{uqJMvf9&;b9d#GszSvN4dG;aHyXE!RRLp?%7H5e~22V;&)Q?scT$)u48&>{UP6
zaGjm@iW?4v)Rja1HgC+?2x7j7s=YVxIh8rG7E4$Cp-a&T+UIbFe+BuN?+4e+d6Ue`
zDM;nG&>Q5CsCPjU+<z9$P--H|9qfknHn>odt-I~4alqKfaR82LwH>jaWBy0rpNxe*
zv$5Z6;rz2MHxjwPVn-fdwA7g#T<wU$pJ;V2#-Et)JqI5zLP8dJbL}zqCfhpXZKiD@
z9_$cDsk4(i5w7S3hwXun3q3oq>gQvV@4%AX>n?Hia~=VQ^4ga<QjpBO;M;QMskdi5
z(d1o7RX<Ogf-N`DyUXo&6~oaz@rw%ai<WuEBEef3x7D@KJ=5DMaKK!*$5>x>B$0fT
zTJK2Y?^tL)lxrWs`x%S(GBdj@=nM>lFLL0D!@R2C8^8ZF9Jd!qyON%}>&?FnqbtRc
z;JQFQoYRbtQ3J2n<24%oLrCv!?h5>-WM4m?e|*k>IkQ9MCMS3JM#H~(jM~U68JT{;
z92W@}9`+5IyD;bM9D_3_lCL+}-xKL^*q6`PQh2Tx(xVtpXboI^dF~eSg3iu49a@z=
zl=XMi^rO~r4;@JUndq<_-zonkxOjvu$=1nMZxrEgPBu!6zP3xY1Y~Fpl5B`K4y{p)
z&$JRflmj*GmIO;T%Wz9I{8?`dvX{e4-I1I-;E6a>M`VEoEryzV@nchMM~yYCkrk{E
zH(p2=c=8DJE9BUZo~`v-_^Ax(yU-o)TII0txjpvT*k6+IE$G?ap1$t0u3c!qW~7*&
zSkfhiwPm=rqS-H?QD-xYLGLtb;}5u2BWbtd+Z8dgiu-Z1x}u?D7Ed9KK9*UMrnw9?
z7y0%F%Yx^Tt)u->*dH}~uHJtN8jbcJ1Cc$C^|g%kEJ?WvE2x68ir%(CqkKlE2IR>{
z`?`3NSa16o8;7iohLWrBC(4lFr@49@{E)&q^Zc8U$*qi&oNqz;H-{QSb)hOU2?ho8
zSf}Usd?VCM@z=8I%FyH!xaTU?PBa>J0-CcH45OZV?Z6VN<F3|o_Tgx+;b^#g`1n%b
zM%z6|#3dXLbG!g8pR;F^AG8$@Qe8^qLH>%Jg4k8j`Vz8uDz@*vcrsB?X|L@I`?bz^
zyoT@5F|twJNZR4v=RlWy;7hHX`xRpf$V4hW?6uB(B2-<_T-v`i7mKmMyAZ!*I?`_v
zvBT5OnP9rtp=+z0Pm<XY1o=K;;)y5NtxrCtx)`a4vuf6M<B5*P_tDO|c1CsqWq*@M
zYtXCi1qO32YgFBgUcz?>*J%$&eyLNYVlMZA2Klf*v>T)PzI33mHbKo#IA*{{ALAbh
zsa0LuC^}AP&l&AxTaTS;z;ViGQ|9dT{xXj|>4n&VlFWZazhz`sgPJPe`8^ZuSMJx1
z>}UNIah?0IyZ3_>6PXINfUn*OigN~>sSaP@SFD{BuKs8532c!g_+>YtN4ol-H>=AE
zSR;1@g)OM6{T8fQ?R2Bv*VGBpAU5-EAPSDmhl4i4mnC39v9`VF%}A#KN1P$j5yy%?
z>B)3wIoq(1>a5$X>#a+zi>y`Fc;l3D-Zl_EE&<D_V9iBik@SN%hr(kq;0(KLvrIXj
z${HEy$Z@XWnOm`GdSi76_lWcrx$0d>?j_)N3&DZ%z<UP3bE}bY8O(VG60|pTs^GbT
zrj50pSx>=cknMBKyO}ZN=1~B@twf40W=;8_$yRt_GCXF$MJ1+GTmx4vMHjd5ME&3c
zd5C9--;3kMV}~5Y0!reig@M%pCw5dLbP!)j3Kzq1E!pKz>0~e$WI$C2>NT)8aDvYq
z_JuV@0<vb2Z#I0C9CE`+%d#za`dg8eeFII5tznJyLbnJxDraQRIUQv6bc9ZW@)QPB
z1AgwgFXLA8Q!=xjh1_l9dGVF!<jv^`8ng)REW?i-ALz~+Y&Y$1J>voYcr4%ZfkDB~
zt$pAg8^;WeM?eR*+1|Ft!AIIlxSTm$fJJb#+XHVNzH;p2dlVek15`*>&eurW{&<w)
zpEPXRd1%1vz<u|U4fq0I@4^;RgmWKgg7W2s6uiv3PxkdRbDg*1MSjWgt(lFt4th({
z<if>Ig5zqXXn&n=kSOAnG%)84XaIFBQ5K0Zg1&*fD@_k`6}#aol5gPYbzr|3NL}T~
zu7#&mU7{X07kMhxD5`?z8~Aet7DW%daYYc5h#@9nzo=fr4IRV%*0rzcolw+4c9r7A
zf5H;UMAANt&GHg{<WBB(C(<O@_W;(4joI%*-s?Bad=1xKOVlI@E-puczQwv#XUIpr
zH<Rn!2`}<#syCNsCk8)2MwAfU>>Fsr16IYes$Z4A@W`AebN135;Z1P**UcVB!a20R
zi~4%2`%oh&eP2F#8pNQ@+6>=S!F6@kVPHYM@g+O3EqAce&tZ>Tz>^#1D#nkAG>?mx
zbNDi8e8O*pY!z^_AKG`r8qYyPZT31%j%?>ldL<#nr4NISi|}zB9_mpv-VW0bS>@R5
zIEUYw13tJE#5SJyn%sM^ZR}iiDE`R^w>Y)MTLVHd0?LZgJFq4)pvY3FC0w-;UONL1
zS7K?|@d{Lzsbi^$sLo)Hj<S;Lp=<|#sp6e;jGFB`5ts(AB?c?;;GFP8CC9eFF8K2R
z{NKpfX7tl~zIVezDd2TESOH?YI5-skJRe984njMkRnfzHxTjS-e{FUo`0F_^)A`v)
zf|mlZcsytQyEuD??*jIo(;bTx+T|YZtqL{cLG}w36J^fGZZbX2c+*~L@LmFAX?6F;
zPLF|;i{W868a9s%gJtCBDdw`>lpQV1lOm{9zF~Z_-T#C2dNA4_@J))r_dMvvDMqI4
zkM@5eT75e!L(vBzdr!O8qcNkcf3q}Lp0_-TcY4g?u?{xwCL;M)$6z#;{0cWS^=I$5
z#?x@_6tJUx)<3XD`askEAc3QZgyzG8E@C7V_z8dWenEWNfnV|>U$1eTB=#-eNIUBU
zO>H<7Bwm=#qh8e;4kqg$oRcq_UPf|+I{~@jIPr(}7F^EVsso|&a@5OQnFFdvn1OYn
zxS^03`5DS~z6Z_s7QCpHu3hKDys4U>vnGn8b?^|M`vtootHSPmG^cc>It9JUD-kZe
z8!AN-H&PTu5%2_b;W1v%Q1PeSD;J)=YJ{|V)dG04KN?@Ti90=UAaN~>yOZyKMI+X;
z);hqMW{~61pL$tUhd&Kc_N-euea|yeyV5@oHQRZDxm+`dXD0rYpC%kC;EzXd`jN5m
z#&)aFo{JW@b+WO}HXQC;hs}5h?Y0V?6ziUCZL!o^HiIC{$Er=R48&KNV4Z<v9^$BT
zW<sB8tT06lI`CbSIpX+4HJ^w?GIs(E+|0NVIBYN8k5+O&UlSO^d0V1A(S8CB4%trF
zPTD7snKBfwVl#SZJ~rqDI5pX~6b!M=Ilwi{HN!-ocHy05U>l|(50@el7op{<ks}lE
z1Lotk$i^Rte#*kml!s7*#2v<b#WNG~yxY*f<I!iW$jeU1!&9Ix<B5EXfHT+e>JG0I
zp>?$iCi{;ByJZids(A<gy9Mqnggejk*(P{;BOXl-*6ReMf4sjoyCwT#c0af-Idl$P
zn}DX97ECu~?s{ZGIjb~-`;>-thvLwM#+-V1E)EZ<3ph^~M%IH$2!+n!nJ)4~IiXU{
z*Bna93FS<t8|Hp^K7>4(#ohMhu9tAftGQ<~-VG<2w5_&xnaxhNB4Z31RGjo8834jy
zlF@|AiJ|1XHkq2Wh$zb%`0gW6-Ko6#;M2^*UVV){h*u%K)JN_BctjLO8=QWp^>yfZ
z74ex)=+ANRXFQ0kvX=iwOn#f0Y5yU2aJTId<8x&G0mBD9za(!jyibQ08z_jaD9k+`
zE&IR~Tvs*q)A6B%pQ~~}yFsePM6>W1+@hL<8+cd!^JIgnqmw#3zQ-C?q(qr?*I?DC
zBmNY$<N)X$&VyB7{uTJ=JxpCHl%U=J9h{*Eh7aCU&ynXi3dxbv@9I>eJzWiCIhHHj
zGas9_g+uY<&G6{Y`Rqv$0`>ZLGoR{VL?$HG>@&opqtKjwFoyzi)|FTJHhNrl^bw<O
z;Ynuuo<<9+lI+_#Bg~$*jw|=^A@_L)8d-bttUyD&5U9a>coIuxn{5u%eA={xZUd7$
z?pO<b%0Y>R7Wel|b`3TTGuvaJ<(p8{hID=vJ*$|@T_AyB+N%rKr`|`P<44vX7`@(l
zlXVqd>Sm;}pLO#~@Sa-g!cMuSV=p)Q?k2W+561z{y~?ul*E@gxj-|#pY2OZij>iX)
zZKS-rV&7>p2D8X{c*n8KF_V137CQcYix>422)xiG?MQVy^mIe3Q6RE0XwfbpPRF6I
zV$1SuWHYJ9Un4JN1pE>0_W=}E?rS5tv+72o*nbjV!<AoS;LH@(!P{nJ>IEd`EMjA~
zL&?Qv4zqeRDNd%inQ|>ZV1=lMf%-~>Gj`^K3+)J=!M|A&7(kz~Pq8Y)`#2V`I-Wv)
zf0z9i`0h?}e<UZ~g({mESB@N2F3nT$$lq8QzhaF{MbiES>HI!cScw_Clg3EnP`>0R
zM8H2_jz5GR7W~8O$a~*ydYmshx|;F9i>&Rb<am8gMq?j%NAb_THmB`DqM~o29rfE+
zwj0>JTGG*ubtRkj2XyeOjC}~2e~x$hsU;e{{DJWy`gkB{?6v66`^|XOjkY!57De!S
zvGq~QQ@^(U`d5}Q#&h=H6MK1$45X#x@tt=5&RtHf(J#pVTYz?b-@b)>?mg65{GQwt
z?QZlSIM8e8R7JgmUh%-AQ1eMMPpBO2tjb#L|F)C4`Gys$D*QCo{`*8DHW9H<b&2+r
zl9p6nWxZ)HO4~jH$0?&WkxEQin`e1XJ+C$A+65;O-s}tB*9OnGafVO)PeHZ5yk|9O
zsQz4u$WC<;QV(HOPAHdqF44GK$l23w0dEAw$;;8QswTOc6;lMxtg3)<$d_^G#0SV{
zcn*1;Km=2Juxh{SDQJxpI9R@ns!^vPS)L6H0|ju}ny@-!@b<UaDv`P?!FXiT9s!Y=
z0G2WYG_B3q2&XlpCyyFQpgghQfk*L9SAz?bJKE5jg`NU*tRfouM9zkpUaTK3bzA#c
z!$|#kFroQawo8ejjRy~0=G|gXx1Z;lNwyhAmbJh#(~@jGYMimx!M&;EA@qcvRUq*d
zSUn*S_fldlan6g5LylvPEUe>HP~uBix=CPzJ&CX!aIN%g3#{??gFA&bMIzZ|^8|;`
zqys^yW5E0;f-LSqMvnvw8U&h}g7w-4;&K$;Tm$v;_$m~!8Glc&I$|J~nad_qH_d{U
z=L19WI!~D22VoJPLKBZhyDA2g1Cm}A+(TA+0a~;6isoEQLc|FZbpiu2(fJyGY);QP
z4UCG$s;#4!TXbk3-tQh_D(RqiHLR_C?nd{uA}2E^Wlm{mJHF{@zD9>Ga6jr=u?PRE
zm6dfVkjv-J1Usy2jPY0>>2T*X`%pZ^TzfZ1B=fBpN()v*KhT|B?!(R$ETuf;>OxR~
zUEp@JIM&&Ql2bAP+vXBBcN*6GY|w)hwh<h&Kxxk+xwQ&9z;!R#x<Q33d=&NEyo7~S
z?G<%e96?8gJ!I=?T*PKtWF29f306A-TmF<Qh0LPEP%|Bl>V-wP*PTl4Sd=S=oQ2uW
z6p#l|b2T=<K?JOhT$OF^B&1Lk*upYw4Pkm|*dCj)Jx+o|%p%U>#0oo)<$uH@TVf-8
z70&&UXKFx|(O65{`AiTS>=<Vdwx<|ZDz>EhDyCx9cEJug4`*7S>IkBq%2g3MxY54|
zEVY6;wtyO*fu_T;S<1*3*$<DV!gVRwOv)=N@~7r(MH^OUFN5L@SWeYgOc#O~d}3&@
zJ0lG+y!G5$ls}ZSEhiR?-XOj?8;tKF7MY~CVk*0^!%{g`gVLr4*IRpn0d;5Ymp~^c
z$7=k|%l2e^i4#OPWm)!un#81y6s+F^aHi;(3B59*W+GWa7aha!P!!`Ug_;@O^Zb1Y
zW7@HoOX199yxFtF*@DphFf)FX(S_JH7rpr|KQWZ;&T>Z=`(&uO%93Nvw3R!yy3T`u
z?S)EyBH2#R&5Q6^12}6LKG-SfSq$xjLFIb35S`rSZs$B5#K~ISMQ}+mkpkscD;^bx
zm0b&MyAgp7!No!o&!bZ-JRvmX3=pC;tj8!gD<4XxK+Drem>u}9!r&D52zr~)XtCfr
z${iB2lL$3i12NdO?fykz^xYV>3d+~{w&0O2W|oISy~5xSd{lJ?Ud(T?fi@F4z6eS#
zo1r>eR_`Gwxfx`4Gx<b&u2{b%P*J#Y5VTm4`nFIBh`W&5*=S8Ww)GaySIOPyfk#In
zTPCn-C!;t2fcG#M&Gn}733(!KLak2VwG-_#$t6*Z)i*rP%b?=hod?j8n~Ya@J%Z)i
zhhq%CWkUObrrw-{SE!m0I((67dy*`oZD_?}>pOU;9wfje+t1B9lb&WKY%ZhId{2V3
zeoxLoz*XyLwb_hUK`-{&jyu#fXee}xrGvpiR+V;2m<pE!k)J;y8#E72Y-YW=ph~!(
zNq2m+e5xZVJnHP9j6C_m%+3pU;CmE`d`xuhX7U$8T(b*j+38V~Qkgz!SV+QQ9wDM7
ztLHs5=}xHm6ZGs66Ol>6#yyK&ycRoW7PVorf#hw<qWuj%k#-Tify|RM-!!vsED`RU
zhAotUEu?NS>T~#T;B#z>pF(-<dtOX6uDWp@!|tsHGgKeGKk^x6tgPk!7l7K{z|4Ql
zI)0rBg(rxK4a4$n!5UJ%-_>SrgR(q@*D4NML;gV(QqXVIL&I2X+EeJlcJyc>QLG%|
zC?}Ymlcvt<aE{0R)R+U%el(g<oO;Swh>RYB{4BxeKZ)j?&wM8neZGkNEJS|hp*f>K
zfmf5QAT%foJva<Wn}}BH=NkmB6Fz@WSG8jcvR?Uu!-#+jonHqvh2u^mQWDAfi^euc
zfNLzoB378$d|o)Sg5NKJmSrMChayXwi7uQ&OLl`Fli|AA_<~uy%JG_46Z_t4YR)sr
z`6#4WJW?5pnHYp3D2hREg}1YO3FHA*5M%F!gzg3Q+wOP5p9|5Ntz<^-$LADYD}S>H
zPKx7-$Kbeg-b<jRh4>ls(TD};&5OK-AZ-gxTxbB8xlk&}&w<czChymQ^o$HC=T144
zV?mH4y=%A=MbnD7uiESh_@)EseR^O{V<-(T^$?%EfYw}$=8xgbm(b)DXwDh%S|_}{
zOe~uoHu2|;$XbitL8ZZ)P_vBH`CGD=e{c@6_s8np%g9&Y%Q}3?HtR#iZP>sUz;$gz
zx7GLgQ6e;xL1P1CQGe&yX<x$VA`|_6oIKrZGs|x~5t=3>^C4!u9Dngq@^x=>Z*XSQ
z8G0PjSQ(wa1Jm0`R(n@+{F>a4nrMgmH!43`JN=Eoe|Z+!kqd(VFHrqyGdu8E_kBd6
zwG+t#WM8;XjgVYr%m%?em0eUxo|qNREC;($1-&}gmf~Fq&sC;Zk3c=LOVQO^`HbSm
zH?YfqvQa*PrhV~7r@?D~gg)A<aTlDl0f{85S+ZTbg{hZMA-+dJb{up%236J9Qv2h)
zgnTLId5_~OEQcZyW;dFv@J~D-!M9-FMB{h<GT81{r@!~$IAs;7dLf!D>A`beo3jz$
zbTGZWCU9QWdrm~!s~^}?*hJ3-#f$C8*An<GimZ|p_$Xxag4D|AYr^u}LIy`BR?vR?
z5qlYy?|7SGtj9`SbR|P*C1aQ3D=3?FEOS`l@e?CcUdI6Br^UVkjHeYeH-r4!W1zDA
zoa>?264uEo@JMOqRm4ihyX(j%5%)&%>1}Y3@_^!zsFNA1xa551_A>F%K|~(6L6h_7
z-UM(J;ki!63MI;PC%_|-<innX>&L>I;V889P&E7|>3s<fo(SKq;*62VDcLzWP^%tV
ziCU>pvx=B|FRb%w?rkhoibFDNzmi3pY0|0;Uf6+T96=1_7@9E1=nU4TV&zA{m7Lg0
z7OWsyIT@zavJf=A0LitI86A#B4>6(%o+{+10Bf2LZyv%|8J=Bfj=W%I-^dPa0<Vzl
zJ`=QO2U!~pV8PL$2K473qV4tAy<6#Z5rftJ3Y7eW;}yKaVOMma6W(0r*k{%*D7Q78
zxo&nokKg}{Bf;Ubuf|q=8WiDMvQB;p(yGk(&RB_xNvbE${m98lXsOO<+h^dnx8ceH
zYL~u)mfw-JSp_v0n$+A4Ua6i!zvT+5C%cM@(0DYZb~ri+x2bc6atxI3*A1O>1YW$I
z_~=_uvxy9rou=e?iz_|gI&Id&yv#i)C-NKKxesb7`#ufaKw5VQxrblDnUgqvYhL%y
z@T^{L^rot*55j-ir>zdOIf$RC$f)Y$!+G<n>iq~$@nJM)BJwa1S=%}AD*jJ7c$E5Z
z9;X*ZDOy;Zu6?jxpdMa1v)@1o?P4<*&g{kvb>rLt<h<&4)$2rFW^YC+D|Rz6>4%8Q
zYj-Beh6ZmYo}_TRZP?L5h!c?vKZmNfvO}POme=mfF3v(dHyF?oq~}?v83iKqxzz?b
zG!bfQcZDw4RNC|MK_u-}U@;!h7zg!X@0tk0HrsS)sm}d**dx=CNXj(MWuLX(Am!=E
zwYB8%w6F@UB1W<g9N;}Nrp0aFQj7O8Xv1w>K^<;$@%5j85BHgsqfNxuzQ^`d9bi0C
z`**}=ny_BAtK1RpNVv~X=&OBzUgwHlthB92<)48WR`7T2{yvU-IS!RnNvEjw=ghJC
z9x2}^mqU2rH{3xvvvm{i8PGBbo8}XIMrElgS4bWCRM`~?8Wx00RcG@O*+WV6JP^0M
zix$=1)2fvC9Irz?hLm%uJ^YHyJj)xI|36|2sdHW-c9444=*-$}=UtOFS7Y~_hZcRI
z*45}+Rr!4bR#A_v)s7!R0mmX!dqF!%>>XHH>QY>2dbANBPxFuv@3I<HNtlcLyv)mD
z--VU45sz{!cB;IIDst^EnO4wZW~&WaHo`aR_#5rCW3z6-Vo3mJQ$2JZcFQzocQZ0)
zCD|Xlpyo6(2iox^CSVc8;T@KOn^$q1CVn=8_s23jsc@JD*^@{f$N*OsXLg%WF<CT8
zSU?8yb1T|=l~)xE6Oprvq2pfUcqVqwP%{pti1sXWY&=ISmRLGt!kI@x-DN*TV*`mp
zmD5?qySdm+Wl&QYkr!AW%Hdsvg};?U*?mP={cikhH`29=`8<Zi&0y8cL^i7$A|JfI
zl1OzI=-|Y=ZH6C9@wJy?Z583)RH5(G6{w1wle1)aSb_twZVH178QaL(SjIf1`N~Y*
z?gD?0M>;FdC<zN{Cf>>xklIt&*2(xj8EB4Xe3J%bc6Ii^>|CrKMa!E4UaX!SW<)P!
z#>Cr+w|k+<N&61_M5vjHy;_47U193YDv*`Errxv=XB^;)M&HFF{RX%iq3lIR1IIZ>
z8rhIhW+sINEmpwQy`B_T4H~QxtFx1PGStt3k{x7g=EG+OucOG>X=Vhzk#&%XCX2xn
z%yX@WmddGB-q>Pn+7R4XLR_;Iye^6>*CTV+u~I7edV(=Cp;R1mvjuI~7c{5@BtF6G
zgwjd;zJzNi`m!FRXMpKN_QMxG0`<=DQofb4R=un`WjeS)gESw?La1T%B4)iG`ZR~H
z7tx%T&Cd*kN29@pVxeX$mShuAjM0pZW=u5e#SPjf>$w(tJ<GJN20%~M_iX`xSOYE9
z4|hM3Hqw`l&J_;T6HXqD%~MB=a!KF_l<!Th=V3mr)z?dQ0JOZo%2C{8Jc!Z(?sGL7
zS8>*LM0=ExdDI`m=6)Z%vjya|1ujdVa&|g+?P9R;caXGS;`Nmw-<|@&e!z7VQZfJ~
zz63#4ztDYVZT$wM<Imt9?O3diQA40+GZt=VtVdM|s_$%gOdi&0KEFS}_dB4UDg(|D
zQ|XGH9D!B(1Xl2Nbe}Ts!v4vZ#7w^Cn!DZ0y)8)etB7b1hPv8|`FrB`s&mjhPe)E4
zh1w7B&Ty~n>7F1GKcFdcsc~~7XZx_ym7^0x*4}{TQeT2Y)THj_9#2t+*dOlvh$_po
zSV8K=uWmhxb==EK_V7;Aa#ogpf8@A&S<gXEYFCXHOf9QAef4)!-nsgysIuu+=0bUQ
z+No|L)V!VSpD|>--^vVCa9?wfrrNW!4;fzS#5{_%5TNJyNzU5CtT0=}o&J<_sB^_r
zQ1c+NSNOUjw4bm>to-c`?o|=o&REjV1_mKzGaRcp?2cZJJ>c0Z;Z>nFqv4u#Qx26u
z&B?AI;Bu-UGN{$cGkLTFPCCtB)T6HnnXFu%nXLF0B<oT*b_Va{f%fFVyJyMO6vuSp
zNQJvBW(0FS)=@O5jsZ6xfcH+ped4+^9EC{KcIaP&Wlf$X_qU4kHX|`A;ky0Ug+lJq
zm<i?OoX2NTHfYdmfC@DuM>?FbSiOoQwjo6oscdAlyu!_Vq5xU91z*=oUd0Nswe5H|
z$|!1w@3h*5s%^&)t%FCDA1+Hr6%t2~(2J0#NpOnza~BwqkU+`WVvZ=}u4L{acr%Sx
z72JB1J`{`5bgSXTew?eExU3uWpaTCUiDSH(q0kc<T|}0`4lLY-fu8uh%Gq1QT{MSA
zgVJtB&MGtB0zaN59=i@%TZXSW0zArzMLQ8~x`hnAYsofNZ;&{Agqc*&{xe?Tuc4-D
zu3qN;y~w<-t_1#?MFxNk+Pua~H6)|W-gg@rsr^^gEi{PlI2UVa19TqDh<uQrH?exO
zzqx$GSn!w!p;#HCKlLigrMe__238D2JwGbZhZWdHUZV7!v7^3p^(T7P)3wO#Y*S?-
zK>g8sU%;8K!wsr{-)Z7T-y%!O*&RYQ(GT$Ex30&Dh$+)gyAx?&%``@TYEni#tA2|w
z<E`ESuU{sXramFshrBcTO%XKpV^Kx3I?QX;EP$tX!kah2>*|2{C#bms-fU#n)x$>7
zL+t?2f$g+`pB1Tj3>`WKj(-<j_X^rr`GVmcVSM<9iFl{s?u*KHP#t%TUlCM|JkGOe
zXArt>(0QcU^hf^&w)6&iv!45U1daSXcI+6gCp%j^?M#Ogv+*9T#aca%v>k`M{W)Dn
zvamcK#7>nD{W82c*fo#3T~SoEU#HQagUv%7DUV+15TZPVa*$M=r+{-*Bd^sZ<1eg$
zZ;_w7@hslrt9BMwWa?AoXDt%+uZ&#AaU8$(2{XrD{U7c#cO2@>8t{T{J<t2uth1l<
z{nt>m(yU^V%_=MLK75Z5REdG^`AFJ*XpU#V8*e7RC5;(Yv`c$^7c<9cj8~2BWAr2W
zEp~T42*p9@J{38q3{35LTVqB&<&U&6!^I}f_6?Xr6L;@LN`6E~lO!aY@=&fPLwX#P
zR{ojxO4I&H@{2Aasnj7-)wrGE;fIN0X_r^^Xq8WN7QaH>i27pDN?wZ6+KqP`lz)Xg
zl@?YHtYY8W)QV|$1?6_01W#;dt<lGaPwrwh9YpSa;D3u=wO0q<N6yyb(f=AB!b2yZ
z@1Q`quUa0u(290bl<h0DnyRq(;5kaqs!NnIb(9UD4#5l2nnSpkQ9MI8s!NVTAN60)
zhS!&|Qq&o7JY2V%m#V&ESmWC1L|0V*q2sKg_gPoVz${|ToPtNW?o;s1&1kk<kmse0
zx(eB`fa4rGJ=~{DS5%F|0@jdvB&1X6c?x+@MW4N8P&1aO-c4|asx4KCtf;a&^NX4r
zs5beD=S#+B;TKJU>Z)8=3q|D}sK<;vmygW#@+cN<nddFm%XsjoN2ybo!Kfs1H(p_+
zb~0B-UG?3)iMvsE@^U^Qt#=<kO~wDY5q^4%Ts(DX76PV@AHRUB#Q&;o7{hs%vxYyX
zrn(0{tSXAj;f`xr8`_EI6Qn^V)Y}IAGNHXZ;Q{=u1baYn^1iGrdb=}+|IED0xBm`0
z)k1}8?2Yiwf3j?ZUw83rb=}~24nM|%jkFUhDFfa+2j6XgqUzq`Wt{f<Ye2VTfiw(s
zyiJbuU7&%BIe!+mjOsZ1Gk-#}m3e$0SeOS|{u6a#4lvq5*k1#|`X2-#lmu<yc~ptn
zg|XD&GwwUMSX?g;_f6!l9d3}X{u{V$xrr=%;Ys$X&&DFC`6ARj=(55|r98L6e93OC
zgSMrt`($h*?arXx2h=t2dt|gaz-Z4)?RFrmN%<uMhzP%jCA6Fz(9g{H%NL$e{9U_r
zwu8WntCZEC{*u@8skxl{5IngYEJ2+!bv@OFX2OBm=T+XT{MHV}t3#Y5rm9bF1n*Uj
z#3x|q+VwhM^54CzzW4B9Ut#UsO>OSa$!M96&nfv{!#TGh?d7Yg+xRYQt?5^4q|ZUY
zbLMO<fVNkm%eG=ms%BVT_E*TGfy}tiq~;&VkG%nUJ_9{HaOl=6)Rd%sk*AJxW?)ZE
zr!U4{av84U3ff)hHninXyu}#kc{_UVr=0l@$OLg_f@6@~30mB0b5Lve6D$q&jTuMh
zw10t{e~VPkLr2GBsrF-~WFtNP0^j9giH1A!-GmNQJ?<OS-#&|W&NF)`c%af}^mS5?
z1Jz6`>-{&-axT=@j<>3D`W_wI%!*uScG}Y}rzOmtfiBcOZz)i-k@=M798NxsvRh(I
z-uwlY<p*$_s%)i=wRd+TG?LYH2wAFFiSobQ=HDd?526{>MM&6}cD((Z@xS5S@J@>A
z)A%ZXSqqKsfmR)?5Xq+o)>wF_`!S|%909i5n=E=+wn^}*3$9frn)*zi4_pU@_Mi_J
zm|E@~cvHP8)YIt{Ga(C3I|qJ-hY|s%cO4n*{rSY(cInn@&ARA(EF8sJWYtJ-R(SqO
zz1Kt&pI5!oZO}6t>`yk;;~d%#$;!NoH_x#mGLW&YRF&0{tyyHuv%lx;0+p2cyd3HF
z8P>8Y_TQtTzQyF>>*zycp^`9<+!c*eR+mtA)yqE)SC1#YHrKff>F^A4`!+I%!d3Av
zF?Vy}WOb>~TG>FCDi^C%{R8(Q5iYY5ROzKXX~#grQlzZ%iDaox!bABTV^7nm=t;M_
zYN*z|Gxlhhn(F*@3zT~m9IzZI5yL7}mZ>tp=R#r0<+X68a;)}aSuWrVs;E<bS|qFG
zA>K_PzOsSrBxS5<Z%5&sxnz<^Hz?yv`v)aKvEC+Ura(<~|ElK>g)V82XIG#Q52cr5
zEd42zsa(y`V8*RFOw4@`)KnDe5*|vq;|!b`L*_?s@)J(l?NHr)CHlO=^f(&0v(5PO
zS<WGjY4&x*&(_cxqlDas#YD6X(DOm2w=SG$2YPnCx7b<BtQXoV8L52TI%_JK2V2P;
zQI50DrHTx7$cS~PfGW*o>?NW$i_x%c<htyEQvE=BB0)`*Cz6XtGL~AUdRH>?VKZY^
zfDWlbOjwYjR;uI5!P{>^b6r3}%pwk-hTkJ(UmXojBSn=(qDWdJKP$U3#~0!*GTm$N
zV<)>-FmgRU&l)6T5;FNPXH%VK9QRR-E>>LN1aW~A$c~v{3U0DdO8m;`ET+a`56?ag
z|0@p7C>)^zOl%~!R!^el=lIDF&UA{A8ARlIv+mCj4OAw$YUW0ARb^b(!<jv~ekXpC
zCv%ABs-kkYQ(&6C0?J&5q6y6FaJ;&5Jizs!Z}K*KnRZz>Z!6RknwM)&gqq5LE#wom
z#6Sn&CoV*SULY2djE8iN*h7wUK2rITZM<!(NzF_mC~3&!R{WBarZ0IJ-?SR9>ok$j
zRgNRn2X2HnmlCVqW~+4^HL0l_mIBtfcEO@v8!y%GN^@SYcOmz7F}b<%SZOPf2+??V
z>$&nWWJm%?p5iS{?w+m&o}`#cLkrgu6_wj2|4$WG@;VMvJ=lvU^J51occL3}*6#`p
zQ_PtDFeGk@nI%#N78rsz2jCOOu)=4P-ID?jAK=J0>tc$b`6AAKlKbuASpgEI%DGiU
zK-!pRW$6mPi)ZdPqhV^XX_c#VhImQ>*Hx{{1Xe?bPn9Uro5q!>scK<VH<m`Wa3rfO
z4=z{UQg_C6!Sj@_stlLz=J!haJ(ZCSk^m($;85AO>+w!v@x=F#XWYUumJB$9iacq~
zq0Hqvj`QS#*|BbhxjL~9z2-fsezPBuio?!h^r0C_th$Cc+b|;Fy^)NQ@hS?TWj8d+
zK-OhBUVJ_}Yqqn2b+QY3c97?nO0={mGgiP#9LW>(Vl{VTWFisig`jsi&LE>NkXg1K
z9N`2rrfMcpei>sG1G6(r`-zdon^c+whm2+ADARfftNa`?PPpnK{Crg;?SMP;@N`>2
zk<?9R0`pqM`C9m`6r5G@?0CP<J&4M_eDIB-@c()~?S~gK;iD-31ne5c*3{3Yfi<Nn
z(L`jkq(Tl9s>d@_Es(s}Xe5FKyY84b1%IU%a<AOv-=$E^&GnW*qa}f(T(=bNRjxp_
z_YAm$s^V1{5`|u0LlpZg^_uEGvkR=F2xQRUuO)tEj@(TY<wIgyTd+Rgq1(1%Qcs#l
zoq8s#SAKu$H=n^u?QMV3b{pB3cT+Dk5lg5A&7=;sw_~BJ54P$&F48~bpOABp+jkTH
zR@O<8aSz;d+@U_0>eelkTODX$#qUstYyoH9YuiM&=10~VV>;Q<+KFC$VC;Bj<4o@?
zpL^<z#j+1;HGzIos=!x&DrH1BptTO6fz;1rDmqU4Kzzp7czhM@nxf8ZLf*8mfchdU
z1}v{d{p)Wh;-g*apEdWpQkO67zA6hupOVJ99xeJkUe(Q5zw*oOfQ!}Ftq4ohMF!`+
zcqorULG|S*rc2CFG;A!@hp%D9KSgJcQbx|f_LP^H##+0b`taY89X^Hr2I+V`>QAGc
zVO3e!2i<oU9%Uw!e-&F`I@;riXJYUZ|0*oP{dm)}P0g?F^1=%f$PBxkozYd9p`IDB
z$io$Iniu|@kBlvamK(`z+{>CC$eI?O97GeHa&30(M{BNN{CKE;9{Nbux{<Z*<YZmO
zCu~Iew815f^oi=_$VDbQiL`AcPbCxD2dS_wA=YvU`8XUNS%iF?$vO9M<^5!0I?<he
z^3A+R<#w|4&Omk5VXMy7?c0nNj053s<POz`ss*H>h}uBaqXb#sst6bh510B<k#a)(
zLww(hrj|TV#CafGt_XHN<dt&lhWWdo;kHrj>qma7p1Xmk6>m}-WXixHtWY6q@&0y@
z$r9IL?4>lWTIf-2^kjU5qu}%EDz%0CR17-Jmku7Oe3h-}*wN^MlZ=f5hnqzdTxfm{
zvg;B$^@2y-c7f<ZI{GmNU8)M$^}Yjs)$TPRnQK7v)%T;9cP05E1N{5RpDBY=Cm>f2
zV^<wT!?qHinP}#Fbpq`}T98@So&3AiWYm_?nc^amuo>X4yG&}T%2Ab)XNaQM9c9ob
zWM<NKA|kq%dGCeIv<SI89$CAJnH+#rn(d1tVs-}2+67sA6wX{j4z()qtI@M5Naj|a
zQ59Z|zCzYVBKB`Onllx-tUB!d_K{#cA*i0tiVy3R^}cS%*95ru5ZR&X4K^FgWg~Xg
zP$-nZdQ)y;1C(m^<%3Nfg+A@(yIWWdikJ+<o7bzF)vY|)W~7=ZJro_6N=#-Yvi%&?
zl&l?zoLvUDL^5ydIXdWXqr9VfY#_ziFW}WSvl{%!H`Q?Gn4eK6l{TMju~t^mLU4xF
z;K{?WjzrB&;vYg({pdsK?>f#M&RI}*k0HohWmhahDvO%u(Tb^N<P?<FHxx9s6B&Q%
z_#zuYnM~WUmolMxBGPsx*zF*yShfX}hj*49L05yehcg655rdozJ-<g993ozzo|_Ah
zwR`X%SDUfz7w{tAg+ABeQ@%|e$y_qizrt>JBPmj_$3{7S?pTae{*ao8+laEgf;9gV
zygU=1W;-~73+bo64}JVEB9B`++cxq9J0X?*WF<XrzmxO1kVcBaRGMSke6{#7LM|dq
zFIPy1hr3iQ%1$JeFgoS2H{i$R!xLYFK8*(RO-B-IXLjwU^aWJau1eY=pcE-T1Kv!>
zVo{CgTObzNX-XMzY2Z2X|CPJbhNLZJ=9L#%3C4I3d0$D*`BW%81|Q)e@X}J|xR{lx
zPQ_i>%e6li?Po~l`;iW}n^mnZ;div-aVlOS<VIBv+ILQRM7!RVqfNqltsVDv&^Pgp
zfc9MY1=zv{Miu+!kp&;_HznDpXm~L`Sutl73MsEsJEJUC^c>!lgfqxZ-w0hZ$&@a`
z#yIRug|A!4`y7B(*FqHMBvxuBS$!Mfwmgm^Vxh6vL%Dbh>Lh|Ki}dj`x(X|32pP}4
z$?-bN$QgV_nVQOUk49D=MUK{*eOl^FnQUQ|Y$S_%0#<E_F%3z(28((Y+NV3xvI0)X
z_ebOPEJikLBQJjt@!0uD+8S(^9Z>5ec77qR;Ygt(zmU3Wu5V{VU+k{2aIWgIwEl&b
zO<?Ag^?$&h$$YD#+aPMN1y6vwEqF5;f2RT7Y+!YshOe6O9x~yy@x*9O!SPMt5$ek~
z3|%!Ev@jW6G}5H`VMZ==D;ku_z2|e@R3uyhqlaOk41+7x5obGGxrj5*#IoCk#eD|p
zB1}cSjKuZBIIFzfX0Cl2s?We)jVJcf51FF6@@CE=%UyksO3gbgBDOh^=t{5by<{;b
zaVR%Egg#u&Y^9U!o{Wx5M0##VB1LhoGEbsi`OfvwP+nwDvQ8SQmaj2((e+|G)Eth~
zFCumy<v(Xe*Obq&+Ezj(;cV)5HyS5wm&vZoAyRe#Jv*Fd?8y#9he06Xpj0clnr`Sh
z+DL?&>PNAH2%72_2I2ek^p{}$_Ctd2;H>FpKE86YL-?YqJ>So922Ugz`4dZ?#4;$K
z3&m5Y_upY!(Xyl^k<}?ecwGvX{v~ARQr1&HB635qGF4+2iOy7RR{?g_Qf5!t?|y$O
zx%2~|j$-F0h`QA=E1KC3W?j|yLh42{a)v+3H4>jdy~(1$8rof|-B<sqB2x#cucRzR
zNpJO5iDz}}LL&Et`U~-C_Tz_Y{Y0~#R4Ffwpgt%CW(N#qjn<-5<Y&lEP^U)ae^=qp
zE<uAPWPinKe-dluZX})RgdRci3cbky*<OhCrV3NpsgHYBQv1Gyao<ACG1#AVSfFE&
zp_gpqz}Z~bLF&F5Kt9}w-C01i<8i22M8>m+nxYz`)fNE~xZR9s39FNq%|{arCZ65G
z*iN=K&iS;_pAI+DPRbO^1v9cE0mrcByR+KggYr7Fx*Q9)Q+Idu|B<{>t<YoO+v#Ah
zUxWYM09vKoCUs2Jd2S#Q^=DH?sz2fwsJRm<+tqa(TkbE;kl7bN*^%L0d{v*_$n%VW
z@>Be~$?)t-lubS8)RRbkgoNJ;S5oEKD5xq_*YMwgXRa*m2cgexXsaR2p!Pf;gN#+L
z;EVjOO5lM14W##O@UiLG--?4z!m^bmsU1tyKQNhdDu$ywZ3GunKB#&`4mY#U6n~q|
z)nDQ+W`JTS%joK0D;n)dvmeRB=tNb2InAz!HxPIFk~q<);Dn8wyTY>!Z}ThO71gC}
zYsmCV<<mYi=hx0kx82;kNxhSfvgSsC+f29LL$y#HG@4;Nh~9h#*{8iUh9VP3u{ve>
zKMR$9=4xT=o6!9+Mz-4CaXf)MD?<-zS2@*S&hV=$*+n<%Eb7EAm^g7Y2>okjw%7vZ
zNO7)Yo^2ELVP8O{)1XY%W)+zWjrSIwQMxIMujOgPnFcf9#t(_{<beZjCnn?NYRVT;
zu9~`Q9wK8;agH&dZb1_<z5`43H6nI5Fzc#oRVR1FAzx&E)&E`y`ekUOh+hnui|TnL
zY){#`zd*l!%1C)i>Rqj!Tt5ck)H4p{epQ{QS`MMI>X@s%GG$gsdw&@CoP2s^`c-gV
zFXI$B9mPxWnror!=j>Fuf(Ts~EKgx&t=OomO{v=#52T5{BypfO)6j3FaI*AfB!2r5
zq}(>F&w8}1x?ZS)<%B6|SG(tvJFcpmo_;?vCJP?&+DF*p4b@p!S+`QTFcWM>U4FK}
znewXR{mM2_4oIBSYu{+(FtWqC&34Rj3JIl*s#;{8Z1O1oYUGeA8QM%udl@OJ9w(Bc
z%gp(zVqVi&tNHZhSp{a>-CYH;dm3y$5xY2yKy|}v9LZHxGo6n%l3YCldN6?1mkB0a
z&9ziZB<on&2D8yhqI|x8I5=A+nPlsU8`ijXpxqX6z4g4mp1D!xNHvt0&3%@GGAKt$
z9RZ`z1mVs*(tg{)LWUsi&wv(PFeU0`@|63s*0#Z$E9i2tg=a2C|Kxz{?cujVaE6oU
z^Cs+B<t6rpl1+@92sg|p%R|}A2^=%Ib_c!LW(B@rXB+iIi@<|_ocPsVICV7^?iW~=
zZ&EWW#9EQJ7|#S^EX$~$oM)l}s^!~9w!1pHePQ1~jpco4lM=N3d%n)}joyk~SxsbQ
z6O|~*)N^hE4|<-e(3@F>JCTj54Bg<HZ003A0mrR`Hy@+cvybtv?K?*-dA5HB6Ht_J
z2>3;*FBXe37K>#mI#Kz85yX2`kNXp>Mx9GF-G|`#8FUx86|X#&de0w_1lou41l;nK
zDY?~~LAzwCzWM|DS1uslc*3+kB`*#^MMXcvp)VuB_3SUf%fEm(wX4yCU;^P<7gbEE
zYmcgHZ-ryt0HG*{V;aH66_<Gs$@~Z6#WrNe8{iw-k5_x^y$mJ2*wO0i*CQ|rRQE|x
zi3f-c^(3CZ8;eie8HptP9x1J?rDxHl6N4pqo3~*_yhe_OYIL+CpZYqex<=hgl%IP&
zQJS~u<NmmLEXLpH51*>H`%Y$lABfCvkrc{i_!U|>8)WYqa-E;UQ@Pvr5Vi6zQ9peE
zo?8Hw)**k=;P?*abpuh)4P>U&q78f4>gauOFSRbc(T7i9-&>h0b>mKBK2%BjEfMd2
z&OP=_+Y~Z-4;YVAS>1%CaWAVsl64f;c9NgFd2a2Jnm}y(ugKYh=&+yi9NN7>7QZ^2
ze$EW90&jT@TB>jE@1Wvz=1@M0I*BV2Ky^ooITjPcC^NgnS290#X!$xcc@&Bcfkq3E
zVX{Hh^}i3AOdPL1smk+G99x~GN|_t&<e)ym%NhR?`tWmPgq`@uZhWj0_*dGKRhcZ0
z5iPuxr`X9cj#+(|Sk`0sejh{M802#o=(61G2&jB(^%~L6YwA!jGL(>A!`d6pb=1E~
zyP2scuJ(dapTWmin>AoVD+1&08S~4CuKgnZ?O(yV`3oZ>Ziwi?E5b6xInVBJ-7z{8
zn(Odp8}ppG5!ZE%h`7oee{)3VPW-#_*}s49%<qvAS4VW!sL_+sv$7(3@u%x}>8&vn
zCXbzxH6de)C2Ms0q_GU`5*|7weezgyq)x7pflxN@Udx^Em%rc17#5`C$M<k$cw|I(
zZdo@y&Kn3!Hv)Oih@Ln1{nvYq>lzW4ud)A=dya_Mov%xo_nbS=$aWZW9hPh(=rr%B
z+rQnB+dD7M;f##9COnpVbpOY3llbeEaaZOcFWWmW*F0B`E9Y`J-E*DJT#M7lo%gSE
zV9>Y39Hq(B96S&i5fd?3Pj&y;abxcvJ7x4(%cKceQzP!=T}`ROm}lfV=S4)s@e7oR
zaCn2x>|8H%*H!QT?<h501Nf&4{o~(%%Nh9P)Jfx5HmO=~i4kyZM8u>`H(wV~s6YHH
zyrOn)`0nP2eVU$8;q~<O6Mg$f%+Mc_BO)~Ljn-%WIiiPt8W;gPM@0N}LCXCRkMcuA
z-|+IXEGu{-Vx9iLM<V&=>~VKV#B_aRNJK=G-tey78d0S`1R^5Da5sK^wkG06y_2pl
z5fS$N2Twn{{qETCZ~yrEUoitF{k#2Wntw#u9})3K1pLv8_!04cw3>cQ<A1-benjkw
zU!M5>zv6lRagP7H88G#K@cMuEA?s~3|NNiQ|6<|)l>V0hkxH6}Ci_1%0Kymje`Wyu
z-6YX1{V#O?--k>!hx`Y%FF~e<@`Uo{=H8s^&9*yl9&*><!Gi`5;eR9&&$K=4Mu2);
zW3FigNS*65e>X_}f7k+C?jQgD`^-Sb=zAt!Z7SJLSA|VQm=~k=rfb4;aP?30y;F})
zvJCXEYrvRi$ql*YT@{hSk5|X*kDabEHTd7(_Fz0;tA)$ItNZ9HC(vTnzi#3o^9)@g
bQuOnc5m!m)>tA?)&S90Iz`vLB{D=MroykEN

literal 0
HcmV?d00001

diff --git a/kerchunk/tests/test_hdf.py b/kerchunk/tests/test_hdf.py
index a0cadce5..fa91c9bf 100644
--- a/kerchunk/tests/test_hdf.py
+++ b/kerchunk/tests/test_hdf.py
@@ -6,6 +6,7 @@
 import pytest
 import xarray as xr
 import zarr
+import h5py
 
 from kerchunk.hdf import SingleHdf5ToZarr
 from kerchunk.combine import MultiZarrToZarr, drop
@@ -92,9 +93,11 @@ def test_multizarr(generate_mzz):
         assert set(ds) == set(expected)
         for name in ds:
             exp = {
-                k: (v.tolist() if v.size > 1 else v[0])
-                if isinstance(v, np.ndarray)
-                else v
+                k: (
+                    (v.tolist() if v.size > 1 else v[0])
+                    if isinstance(v, np.ndarray)
+                    else v
+                )
                 for k, v in expected[name].attrs.items()
             }
             assert dict(ds[name].attrs) == dict(exp)
@@ -331,3 +334,27 @@ def test_inline_threshold():
         fn, inline_threshold=1e9
     ).translate()
     assert inline_0 != inline_1_million
+
+
+@pytest.mark.skipif(
+    not h5py.__version__ >= "3.11.0",
+    reason="'h5py.Group.visititems_links' requires h5py 3.11.0 or later",
+)
+def test_translate_links():
+    fn = osp.join(here, "air_linked.nc")
+    # choose a threshold that will give both inline and non-inline
+    # datasets for maximum test coverage
+    out = kerchunk.hdf.SingleHdf5ToZarr(fn, inline_threshold=50).translate(
+        preserve_linked_dsets=True
+    )
+    fs = fsspec.filesystem("reference", fo=out)
+    z = zarr.open(fs.get_mapper())
+
+    # 1. Test the hard linked datasets were translated correctly
+    # 2. Test the soft linked datasets were translated correctly
+    for link in ("hard", "soft"):
+        for dset in ("lat", "time"):
+            assert z[f"{dset}_{link}"].shape == z[dset].shape
+            for key in z[f"{dset}_{link}"].attrs.keys():
+                if key not in kerchunk.hdf._HIDDEN_ATTRS and key != "_ARRAY_DIMENSIONS":
+                    assert z[f"{dset}_{link}"].attrs[key] == z[dset].attrs[key]

From 5a48720f3cd017c53bca626f50a49fecca6c463d Mon Sep 17 00:00:00 2001
From: Lawson Woods <ljwoods2@asu.edu>
Date: Sat, 15 Jun 2024 21:24:24 -0700
Subject: [PATCH 02/13] require instead of create group/dset

---
 kerchunk/hdf.py | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/kerchunk/hdf.py b/kerchunk/hdf.py
index 0339be4a..0783ed5e 100644
--- a/kerchunk/hdf.py
+++ b/kerchunk/hdf.py
@@ -462,7 +462,7 @@ def _translator(
                         )
 
                 # Create a Zarr array equivalent to this HDF5 dataset...
-                za = self._zroot.create_dataset(
+                za = self._zroot.require_dataset(
                     h5obj.name,
                     shape=h5obj.shape,
                     dtype=dt or h5obj.dtype,
@@ -510,7 +510,7 @@ def _translator(
 
             elif isinstance(h5obj, h5py.Group):
                 lggr.debug(f"HDF5 group: {h5obj.name}")
-                zgrp = self._zroot.create_group(h5obj.name)
+                zgrp = self._zroot.require_group(h5obj.name)
                 self._transfer_attrs(h5obj, zgrp)
         except Exception as e:
             import traceback

From 2d73e5768860ac50bb70d743578593b2a8676245 Mon Sep 17 00:00:00 2001
From: Lawson Woods <ljwoods2@asu.edu>
Date: Sat, 15 Jun 2024 21:28:24 -0700
Subject: [PATCH 03/13] spelling

---
 kerchunk/hdf.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/kerchunk/hdf.py b/kerchunk/hdf.py
index 0783ed5e..e578bae6 100644
--- a/kerchunk/hdf.py
+++ b/kerchunk/hdf.py
@@ -141,7 +141,7 @@ def translate(self, preserve_linked_dsets=False):
         if preserve_linked_dsets:
             if h5py.__version__ < "3.11.0":
                 raise RuntimeError(
-                    "preserve_links requires h5py 3.11.0 or later, "
+                    "preserve_linked_dsets requires h5py 3.11.0 or later, "
                     f"found {h5py.__version__}"
                 )
             self._h5f.visititems_links(self._translator)

From a631bcb23c42e6f2d4c0853663a53c5979f0f0fd Mon Sep 17 00:00:00 2001
From: ljwoods2 <ljwoods2@asu.edu>
Date: Tue, 25 Jun 2024 09:38:08 -0700
Subject: [PATCH 04/13] test dataset content rather than shape

---
 kerchunk/tests/test_hdf.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/kerchunk/tests/test_hdf.py b/kerchunk/tests/test_hdf.py
index fa91c9bf..1e0ec06b 100644
--- a/kerchunk/tests/test_hdf.py
+++ b/kerchunk/tests/test_hdf.py
@@ -354,7 +354,7 @@ def test_translate_links():
     # 2. Test the soft linked datasets were translated correctly
     for link in ("hard", "soft"):
         for dset in ("lat", "time"):
-            assert z[f"{dset}_{link}"].shape == z[dset].shape
+            np.testing.assert_allclose(z[dset], z[f"{dset}_{link}"])
             for key in z[f"{dset}_{link}"].attrs.keys():
                 if key not in kerchunk.hdf._HIDDEN_ATTRS and key != "_ARRAY_DIMENSIONS":
                     assert z[f"{dset}_{link}"].attrs[key] == z[dset].attrs[key]

From 89567e85d80ee6599f9252cf30c393eb1bcb7bc3 Mon Sep 17 00:00:00 2001
From: Lawson Woods <ljwoods2@asu.edu>
Date: Sun, 30 Jun 2024 12:02:56 -0700
Subject: [PATCH 05/13] improved version compatibility check

---
 kerchunk/hdf.py            | 12 ++++++----
 kerchunk/tests/test_hdf.py | 45 ++++++++++++++++++++++++++++----------
 2 files changed, 42 insertions(+), 15 deletions(-)

diff --git a/kerchunk/hdf.py b/kerchunk/hdf.py
index e578bae6..3885060a 100644
--- a/kerchunk/hdf.py
+++ b/kerchunk/hdf.py
@@ -2,6 +2,7 @@
 import io
 import logging
 from typing import Union, BinaryIO
+from importlib.metadata import version
 
 import fsspec.core
 from fsspec.implementations.reference import LazyReferenceMapper
@@ -139,10 +140,12 @@ def translate(self, preserve_linked_dsets=False):
         self._h5f.visititems(self._translator)
 
         if preserve_linked_dsets:
-            if h5py.__version__ < "3.11.0":
+            try:
+                h5py.Group.visititems_links
+            except AttributeError:
                 raise RuntimeError(
-                    "preserve_linked_dsets requires h5py 3.11.0 or later, "
-                    f"found {h5py.__version__}"
+                    "'preserve_linked_dsets' kwarg requires h5py 3.11.0 or later "
+                    f"is installed, found {version("h5py")}"
                 )
             self._h5f.visititems_links(self._translator)
 
@@ -275,10 +278,11 @@ def _translator(
         try:  # method must not raise exception
             kwargs = {}
 
-            if isinstance(h5obj, h5py.SoftLink) or isinstance(h5obj, h5py.HardLink):
+            if isinstance(h5obj, (h5py.SoftLink, h5py.HardLink)):
                 h5obj = self._h5f[name]
                 if isinstance(h5obj, h5py.Group):
                     # continues iteration of visititems_links
+                    lggr.debug(f"Skipping translation of HDF5 linked group: '{h5obj.name}'")
                     return None
 
             if isinstance(h5obj, h5py.Dataset):
diff --git a/kerchunk/tests/test_hdf.py b/kerchunk/tests/test_hdf.py
index 1e0ec06b..51f1d63e 100644
--- a/kerchunk/tests/test_hdf.py
+++ b/kerchunk/tests/test_hdf.py
@@ -25,7 +25,9 @@ def test_single():
     m = fsspec.get_mapper(
         "reference://", fo=test_dict, remote_protocol="s3", remote_options=so
     )
-    ds = xr.open_dataset(m, engine="zarr", backend_kwargs=dict(consolidated=False))
+    ds = xr.open_dataset(
+        m, engine="zarr", backend_kwargs=dict(consolidated=False)
+    )
 
     with fsspec.open(url, **so) as f:
         expected = xr.open_dataset(f, engine="h5netcdf")
@@ -84,7 +86,9 @@ def test_multizarr(generate_mzz):
     m = fsspec.get_mapper(
         "reference://", fo=test_dict, remote_protocol="s3", remote_options=so
     )
-    ds = xr.open_dataset(m, engine="zarr", backend_kwargs=dict(consolidated=False))
+    ds = xr.open_dataset(
+        m, engine="zarr", backend_kwargs=dict(consolidated=False)
+    )
 
     with fsspec.open_files(urls, **so) as fs:
         expts = [xr.open_dataset(f, engine="h5netcdf") for f in fs]
@@ -132,9 +136,9 @@ def times_data(tmpdir):
     lon = xr.DataArray(np.linspace(-90, 90, 10), dims=["lon"], name="lon")
     time_attrs = {"axis": "T", "long_name": "time", "standard_name": "time"}
     time1 = xr.DataArray(
-        np.arange(-631108800000000000, -630158390000000000, 86400000000000).view(
-            "datetime64[ns]"
-        ),
+        np.arange(
+            -631108800000000000, -630158390000000000, 86400000000000
+        ).view("datetime64[ns]"),
         dims=["time"],
         name="time",
         attrs=time_attrs,
@@ -162,7 +166,9 @@ def test_times(times_data):
         "reference://",
         fo=test_dict,
     )
-    result = xr.open_dataset(m, engine="zarr", backend_kwargs=dict(consolidated=False))
+    result = xr.open_dataset(
+        m, engine="zarr", backend_kwargs=dict(consolidated=False)
+    )
     expected = x1.to_dataset()
     xr.testing.assert_equal(result, expected)
 
@@ -178,7 +184,9 @@ def test_times_str(times_data):
         "reference://",
         fo=test_dict,
     )
-    result = xr.open_dataset(m, engine="zarr", backend_kwargs=dict(consolidated=False))
+    result = xr.open_dataset(
+        m, engine="zarr", backend_kwargs=dict(consolidated=False)
+    )
     expected = x1.to_dataset()
     xr.testing.assert_equal(result, expected)
 
@@ -201,7 +209,9 @@ def test_string_embed():
 
 def test_string_null():
     fn = osp.join(here, "vlen.h5")
-    h = kerchunk.hdf.SingleHdf5ToZarr(fn, fn, vlen_encode="null", inline_threshold=0)
+    h = kerchunk.hdf.SingleHdf5ToZarr(
+        fn, fn, vlen_encode="null", inline_threshold=0
+    )
     out = h.translate()
     fs = fsspec.filesystem("reference", fo=out)
     z = zarr.open(fs.get_mapper())
@@ -240,7 +250,9 @@ def test_string_decode():
 def test_compound_string_null():
     fn = osp.join(here, "vlen2.h5")
     with open(fn, "rb") as f:
-        h = kerchunk.hdf.SingleHdf5ToZarr(f, fn, vlen_encode="null", inline_threshold=0)
+        h = kerchunk.hdf.SingleHdf5ToZarr(
+            f, fn, vlen_encode="null", inline_threshold=0
+        )
         out = h.translate()
     fs = fsspec.filesystem("reference", fo=out)
     z = zarr.open(fs.get_mapper())
@@ -336,8 +348,16 @@ def test_inline_threshold():
     assert inline_0 != inline_1_million
 
 
+def has_visititems_links():
+    try:
+        h5py.Group.visititems_links
+    except AttributeError:
+        return False
+    return True
+
+
 @pytest.mark.skipif(
-    not h5py.__version__ >= "3.11.0",
+    not has_visititems_links(),
     reason="'h5py.Group.visititems_links' requires h5py 3.11.0 or later",
 )
 def test_translate_links():
@@ -356,5 +376,8 @@ def test_translate_links():
         for dset in ("lat", "time"):
             np.testing.assert_allclose(z[dset], z[f"{dset}_{link}"])
             for key in z[f"{dset}_{link}"].attrs.keys():
-                if key not in kerchunk.hdf._HIDDEN_ATTRS and key != "_ARRAY_DIMENSIONS":
+                if (
+                    key not in kerchunk.hdf._HIDDEN_ATTRS
+                    and key != "_ARRAY_DIMENSIONS"
+                ):
                     assert z[f"{dset}_{link}"].attrs[key] == z[dset].attrs[key]

From 561c829a98d7adfe8a657dffba4d676b4f754804 Mon Sep 17 00:00:00 2001
From: Lawson Woods <ljwoods2@asu.edu>
Date: Sun, 30 Jun 2024 12:08:36 -0700
Subject: [PATCH 06/13] undo black autoformat

---
 kerchunk/tests/test_hdf.py | 30 +++++++++---------------------
 1 file changed, 9 insertions(+), 21 deletions(-)

diff --git a/kerchunk/tests/test_hdf.py b/kerchunk/tests/test_hdf.py
index 51f1d63e..ae80c5bf 100644
--- a/kerchunk/tests/test_hdf.py
+++ b/kerchunk/tests/test_hdf.py
@@ -25,9 +25,7 @@ def test_single():
     m = fsspec.get_mapper(
         "reference://", fo=test_dict, remote_protocol="s3", remote_options=so
     )
-    ds = xr.open_dataset(
-        m, engine="zarr", backend_kwargs=dict(consolidated=False)
-    )
+    ds = xr.open_dataset(m, engine="zarr", backend_kwargs=dict(consolidated=False))
 
     with fsspec.open(url, **so) as f:
         expected = xr.open_dataset(f, engine="h5netcdf")
@@ -86,9 +84,7 @@ def test_multizarr(generate_mzz):
     m = fsspec.get_mapper(
         "reference://", fo=test_dict, remote_protocol="s3", remote_options=so
     )
-    ds = xr.open_dataset(
-        m, engine="zarr", backend_kwargs=dict(consolidated=False)
-    )
+    ds = xr.open_dataset(m, engine="zarr", backend_kwargs=dict(consolidated=False))
 
     with fsspec.open_files(urls, **so) as fs:
         expts = [xr.open_dataset(f, engine="h5netcdf") for f in fs]
@@ -136,9 +132,9 @@ def times_data(tmpdir):
     lon = xr.DataArray(np.linspace(-90, 90, 10), dims=["lon"], name="lon")
     time_attrs = {"axis": "T", "long_name": "time", "standard_name": "time"}
     time1 = xr.DataArray(
-        np.arange(
-            -631108800000000000, -630158390000000000, 86400000000000
-        ).view("datetime64[ns]"),
+        np.arange(-631108800000000000, -630158390000000000, 86400000000000).view(
+            "datetime64[ns]"
+        ),
         dims=["time"],
         name="time",
         attrs=time_attrs,
@@ -166,9 +162,7 @@ def test_times(times_data):
         "reference://",
         fo=test_dict,
     )
-    result = xr.open_dataset(
-        m, engine="zarr", backend_kwargs=dict(consolidated=False)
-    )
+    result = xr.open_dataset(m, engine="zarr", backend_kwargs=dict(consolidated=False))
     expected = x1.to_dataset()
     xr.testing.assert_equal(result, expected)
 
@@ -184,9 +178,7 @@ def test_times_str(times_data):
         "reference://",
         fo=test_dict,
     )
-    result = xr.open_dataset(
-        m, engine="zarr", backend_kwargs=dict(consolidated=False)
-    )
+    result = xr.open_dataset(m, engine="zarr", backend_kwargs=dict(consolidated=False))
     expected = x1.to_dataset()
     xr.testing.assert_equal(result, expected)
 
@@ -209,9 +201,7 @@ def test_string_embed():
 
 def test_string_null():
     fn = osp.join(here, "vlen.h5")
-    h = kerchunk.hdf.SingleHdf5ToZarr(
-        fn, fn, vlen_encode="null", inline_threshold=0
-    )
+    h = kerchunk.hdf.SingleHdf5ToZarr(fn, fn, vlen_encode="null", inline_threshold=0)
     out = h.translate()
     fs = fsspec.filesystem("reference", fo=out)
     z = zarr.open(fs.get_mapper())
@@ -250,9 +240,7 @@ def test_string_decode():
 def test_compound_string_null():
     fn = osp.join(here, "vlen2.h5")
     with open(fn, "rb") as f:
-        h = kerchunk.hdf.SingleHdf5ToZarr(
-            f, fn, vlen_encode="null", inline_threshold=0
-        )
+        h = kerchunk.hdf.SingleHdf5ToZarr(f, fn, vlen_encode="null", inline_threshold=0)
         out = h.translate()
     fs = fsspec.filesystem("reference", fo=out)
     z = zarr.open(fs.get_mapper())

From de48c1cd3a577dd674851c42c097419035a2d2cf Mon Sep 17 00:00:00 2001
From: Lawson Woods <ljwoods2@asu.edu>
Date: Sun, 30 Jun 2024 18:03:16 -0700
Subject: [PATCH 07/13] has_visititems_links suggestions

---
 kerchunk/hdf.py            |  7 ++++---
 kerchunk/tests/test_hdf.py | 10 +---------
 2 files changed, 5 insertions(+), 12 deletions(-)

diff --git a/kerchunk/hdf.py b/kerchunk/hdf.py
index 3885060a..2ee8081d 100644
--- a/kerchunk/hdf.py
+++ b/kerchunk/hdf.py
@@ -140,9 +140,7 @@ def translate(self, preserve_linked_dsets=False):
         self._h5f.visititems(self._translator)
 
         if preserve_linked_dsets:
-            try:
-                h5py.Group.visititems_links
-            except AttributeError:
+            if not has_visititems_links():
                 raise RuntimeError(
                     "'preserve_linked_dsets' kwarg requires h5py 3.11.0 or later "
                     f"is installed, found {version("h5py")}"
@@ -673,3 +671,6 @@ def _is_netcdf_datetime(dataset: h5py.Dataset):
 
 def _is_netcdf_variable(dataset: h5py.Dataset):
     return any("_Netcdf4" in _ for _ in dataset.attrs)
+
+def has_visititems_links():
+    return hasattr(h5py.Group, "visititems_links")
diff --git a/kerchunk/tests/test_hdf.py b/kerchunk/tests/test_hdf.py
index ae80c5bf..b352dc99 100644
--- a/kerchunk/tests/test_hdf.py
+++ b/kerchunk/tests/test_hdf.py
@@ -8,7 +8,7 @@
 import zarr
 import h5py
 
-from kerchunk.hdf import SingleHdf5ToZarr
+from kerchunk.hdf import SingleHdf5ToZarr, has_visititems_links
 from kerchunk.combine import MultiZarrToZarr, drop
 
 here = osp.dirname(__file__)
@@ -336,14 +336,6 @@ def test_inline_threshold():
     assert inline_0 != inline_1_million
 
 
-def has_visititems_links():
-    try:
-        h5py.Group.visititems_links
-    except AttributeError:
-        return False
-    return True
-
-
 @pytest.mark.skipif(
     not has_visititems_links(),
     reason="'h5py.Group.visititems_links' requires h5py 3.11.0 or later",

From 5ebb0f4aa6efae73018627eb8f557db59abe5ae5 Mon Sep 17 00:00:00 2001
From: Martin Durant <martindurant@users.noreply.github.com>
Date: Sun, 30 Jun 2024 22:28:39 -0400
Subject: [PATCH 08/13] Update kerchunk/hdf.py

---
 kerchunk/hdf.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/kerchunk/hdf.py b/kerchunk/hdf.py
index 2ee8081d..9dd52744 100644
--- a/kerchunk/hdf.py
+++ b/kerchunk/hdf.py
@@ -143,7 +143,7 @@ def translate(self, preserve_linked_dsets=False):
             if not has_visititems_links():
                 raise RuntimeError(
                     "'preserve_linked_dsets' kwarg requires h5py 3.11.0 or later "
-                    f"is installed, found {version("h5py")}"
+                    f"is installed, found {h5py.__version__)}"
                 )
             self._h5f.visititems_links(self._translator)
 

From 3bccb5ac47f9411421bc8ac885cfa5d853422ea5 Mon Sep 17 00:00:00 2001
From: Martin Durant <martindurant@users.noreply.github.com>
Date: Sun, 30 Jun 2024 22:30:35 -0400
Subject: [PATCH 09/13] Update kerchunk/hdf.py

---
 kerchunk/hdf.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/kerchunk/hdf.py b/kerchunk/hdf.py
index 9dd52744..b1d06fd7 100644
--- a/kerchunk/hdf.py
+++ b/kerchunk/hdf.py
@@ -143,7 +143,7 @@ def translate(self, preserve_linked_dsets=False):
             if not has_visititems_links():
                 raise RuntimeError(
                     "'preserve_linked_dsets' kwarg requires h5py 3.11.0 or later "
-                    f"is installed, found {h5py.__version__)}"
+                    f"is installed, found {h5py.__version__}"
                 )
             self._h5f.visititems_links(self._translator)
 

From 6964b405456d52ff5c3c63dc03762f3e38e209a3 Mon Sep 17 00:00:00 2001
From: Martin Durant <martindurant@users.noreply.github.com>
Date: Sun, 30 Jun 2024 22:32:36 -0400
Subject: [PATCH 10/13] Update kerchunk/tests/test_hdf.py

lint
---
 kerchunk/tests/test_hdf.py | 5 +----
 1 file changed, 1 insertion(+), 4 deletions(-)

diff --git a/kerchunk/tests/test_hdf.py b/kerchunk/tests/test_hdf.py
index b352dc99..a20be2ad 100644
--- a/kerchunk/tests/test_hdf.py
+++ b/kerchunk/tests/test_hdf.py
@@ -356,8 +356,5 @@ def test_translate_links():
         for dset in ("lat", "time"):
             np.testing.assert_allclose(z[dset], z[f"{dset}_{link}"])
             for key in z[f"{dset}_{link}"].attrs.keys():
-                if (
-                    key not in kerchunk.hdf._HIDDEN_ATTRS
-                    and key != "_ARRAY_DIMENSIONS"
-                ):
+                if key not in kerchunk.hdf._HIDDEN_ATTRS and key != "_ARRAY_DIMENSIONS":
                     assert z[f"{dset}_{link}"].attrs[key] == z[dset].attrs[key]

From 6e48c41b06bbc4e97a7627a421bc41800873abbe Mon Sep 17 00:00:00 2001
From: Martin Durant <martindurant@users.noreply.github.com>
Date: Sun, 30 Jun 2024 22:41:58 -0400
Subject: [PATCH 11/13] Update kerchunk/hdf.py

---
 kerchunk/hdf.py | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/kerchunk/hdf.py b/kerchunk/hdf.py
index b1d06fd7..bdab55a9 100644
--- a/kerchunk/hdf.py
+++ b/kerchunk/hdf.py
@@ -280,7 +280,9 @@ def _translator(
                 h5obj = self._h5f[name]
                 if isinstance(h5obj, h5py.Group):
                     # continues iteration of visititems_links
-                    lggr.debug(f"Skipping translation of HDF5 linked group: '{h5obj.name}'")
+                    lggr.debug(
+                        f"Skipping translation of HDF5 linked group: '{h5obj.name}'"
+                    )
                     return None
 
             if isinstance(h5obj, h5py.Dataset):

From 9a5499aaee71fc5bd1745649c0e75d1c7057c227 Mon Sep 17 00:00:00 2001
From: Martin Durant <martindurant@users.noreply.github.com>
Date: Sun, 30 Jun 2024 22:45:26 -0400
Subject: [PATCH 12/13] Update kerchunk/hdf.py

---
 kerchunk/hdf.py | 1 -
 1 file changed, 1 deletion(-)

diff --git a/kerchunk/hdf.py b/kerchunk/hdf.py
index bdab55a9..8fad60ef 100644
--- a/kerchunk/hdf.py
+++ b/kerchunk/hdf.py
@@ -2,7 +2,6 @@
 import io
 import logging
 from typing import Union, BinaryIO
-from importlib.metadata import version
 
 import fsspec.core
 from fsspec.implementations.reference import LazyReferenceMapper

From 3f7cacb6a7cfaec95f75fb23ac7fae141ea21963 Mon Sep 17 00:00:00 2001
From: Martin Durant <martindurant@users.noreply.github.com>
Date: Sun, 30 Jun 2024 22:45:58 -0400
Subject: [PATCH 13/13] Update kerchunk/hdf.py

---
 kerchunk/hdf.py | 1 +
 1 file changed, 1 insertion(+)

diff --git a/kerchunk/hdf.py b/kerchunk/hdf.py
index 8fad60ef..b93a8937 100644
--- a/kerchunk/hdf.py
+++ b/kerchunk/hdf.py
@@ -673,5 +673,6 @@ def _is_netcdf_datetime(dataset: h5py.Dataset):
 def _is_netcdf_variable(dataset: h5py.Dataset):
     return any("_Netcdf4" in _ for _ in dataset.attrs)
 
+
 def has_visititems_links():
     return hasattr(h5py.Group, "visititems_links")