From bafbd98c3df94acfcd1f571752871c6eadf346e6 Mon Sep 17 00:00:00 2001
From: Sander Roosendaal <sanderroosendaal@proton.me>
Date: Thu, 13 Jun 2024 15:27:11 +0200
Subject: [PATCH] fixing filtering

---
 rowers/dataroutines.py                |  18 ++++++++++++++----
 rowers/interactiveplots.py            |  15 ++++++++++-----
 rowers/metrics.py                     |   1 +
 rowers/tests/testdata/testdata.tcx.gz | Bin 4000 -> 4001 bytes
 4 files changed, 25 insertions(+), 9 deletions(-)

diff --git a/rowers/dataroutines.py b/rowers/dataroutines.py
index 92c83cb7..154d75f8 100644
--- a/rowers/dataroutines.py
+++ b/rowers/dataroutines.py
@@ -1512,7 +1512,7 @@ def getrowdata_pl(id=0, doclean=False, convertnewtons=True,
 
     
 def read_data(columns, ids=[], doclean=True, workstrokesonly=True, debug=False, for_chart=False, compute=True,
-              startenddict={}):
+              startenddict={}, driveenergy=True):
     if ids:
         csvfilenames = [
             'media/strokedata_{id}.parquet.gz'.format(id=id) for id in ids]
@@ -1520,7 +1520,9 @@ def read_data(columns, ids=[], doclean=True, workstrokesonly=True, debug=False,
         return pl.DataFrame()
 
     data = []
-    columns = [c for c in columns if c != 'None'] + ['distance', 'spm', 'workoutid','workoutstate', 'driveenergy']
+    columns = [c for c in columns if c != 'None'] + ['distance', 'spm', 'workoutid','workoutstate']
+    if driveenergy:
+        columns = columns + ['driveenergy']
     columns = list(set(columns))
 
     for id, f in zip(ids, csvfilenames):
@@ -1583,9 +1585,10 @@ def read_data(columns, ids=[], doclean=True, workstrokesonly=True, debug=False,
     if len(data)==0:
         return pl.DataFrame()
 
+    
     try:
         datadf = pl.concat(data).select(columns)
-    except (SchemaError, ShapeError):
+    except (ShapeError, SchemaError):
         data = [
             df.select(columns)
             for df in data]
@@ -1593,6 +1596,7 @@ def read_data(columns, ids=[], doclean=True, workstrokesonly=True, debug=False,
         # float columns
         floatcolumns = []
         intcolumns = []
+        stringcolumns = []
         for c in columns:
             try:
                 if metricsdicts[c]['numtype'] == 'float':
@@ -1600,7 +1604,11 @@ def read_data(columns, ids=[], doclean=True, workstrokesonly=True, debug=False,
                 if metricsdicts[c]['numtype'] == 'integer':
                     intcolumns.append(c)
             except KeyError:
-                pass
+                if c[0] == 'f':
+                    stringcolumns.append(c)
+                else:
+                    intcolumns.append(c)
+                
         data = [
             df.with_columns(
                 cs.float().cast(pl.Float64)
@@ -1610,6 +1618,8 @@ def read_data(columns, ids=[], doclean=True, workstrokesonly=True, debug=False,
                 cs.by_name(intcolumns).cast(pl.Int64)
             ).with_columns(
                 cs.by_name(floatcolumns).cast(pl.Float64)
+            ).with_columns(
+                cs.by_name(stringcolumns).cast(pl.String)
             )
             for df in data
         ]
diff --git a/rowers/interactiveplots.py b/rowers/interactiveplots.py
index af49c155..2a3769ce 100644
--- a/rowers/interactiveplots.py
+++ b/rowers/interactiveplots.py
@@ -1802,7 +1802,7 @@ def interactive_flexchart_stacked(id, r, xparam='time',
         'metrics': metrics_list,
     }
 
-    script, div = get_chart("/stacked", chart_data, debug=False)
+    script, div = get_chart("/stacked", chart_data, debug=True)
 
     return script, div
 
@@ -2129,16 +2129,20 @@ def interactive_multiple_compare_chart(ids, xparam, yparam, plottype='line',
         compute = True
         doclean = True
 
+    driveenergy = False
+    if 'driveenergy' in columns:
+        driveenergy = True
+
 
     datadf = pl.DataFrame()
     if promember:
         datadf = dataprep.read_data(columns, ids=ids, doclean=doclean,
-                                    compute=compute,
+                                    compute=compute, driveenergy=driveenergy,
                                     workstrokesonly=workstrokesonly, for_chart=True,
                                     startenddict=startenddict)
     else:
         datadf = dataprep.read_data(columns_basic, ids=ids, doclean=doclean,
-                                    compute=compute,
+                                    compute=compute, driveenergy = driveenergy,
                                     workstrokesonly=workstrokesonly, for_chart=True,
                                     startenddict=startenddict)
 
@@ -2165,10 +2169,10 @@ def interactive_multiple_compare_chart(ids, xparam, yparam, plottype='line',
 
     if (xparam == 'time'):
         datadf = datadf.with_columns((pl.col(xparam)-datadf[0,xparam]).alias(xparam))
-    
+
     data_dict = datadf.to_dicts()
 
-    metrics_list = [{'name': name, 'rowingmetrics':d } for name, d in metrics.rowingmetrics]
+    metrics_list = [{'name': name, 'rowingmetrics':{k: v for k, v in d.items() if k != 'numtype_pl'} } for name, d in metrics.rowingmetrics]
 
     try:
         workoutsdict = [{'id': id, 'label': labeldict[id]} for id in ids]
@@ -2185,6 +2189,7 @@ def interactive_multiple_compare_chart(ids, xparam, yparam, plottype='line',
         'workouts': workoutsdict,
     }
 
+
     script, div = get_chart("/compare", chart_data, debug=False)
     return script, div
 
diff --git a/rowers/metrics.py b/rowers/metrics.py
index 3766addb..83299c9c 100644
--- a/rowers/metrics.py
+++ b/rowers/metrics.py
@@ -5,6 +5,7 @@ import numpy as np
 import pandas as pd
 from scipy import optimize
 from django.utils import timezone
+from polars import String, Int32, Float64
 
 from math import log10, log2
 from rowers.mytypes import otwtypes, otetypes
diff --git a/rowers/tests/testdata/testdata.tcx.gz b/rowers/tests/testdata/testdata.tcx.gz
index 8058cd3557d950a0012ee8831c8091f80bed9d17..cca9807a60c3a074c80a095ef2d65c35e1781136 100644
GIT binary patch
delta 3914
zcmV-Q54G^1AE6%xABzYG;OvnGAQLJz=TfIGj&VdFt3?7f3MiYmRgpO(f1lo+p6tC@
zUTjup>w9<6z}>x{4(}X1zvx!$)%y6sa@{X49(H~AqT4L@($4G6y}Pga{(OJ2c>DIP
zIqoh_S8H>0cDhJkZjN3pPrJ?H`okA*=<e<m``x>5ILrN`e)VS6zyIZ3zg*`VymEqf
zo7HD$c+2MO;&?#-xcKAA72{{yQ~a=6Z~AV1v`j~BR{Oj2K0cA~7Xi?d2mvvFJ3qVV
z@7?`p|7v-0clTxb`q9fn0(3us{6O#=dEeqbITIamguemo0~{>UP9OJw+@0?|>${8o
z`Ra6e@9tk6dwHRB=<~C_JIRMVTOOUQUv3VC{lVg^t)C8fxZ1rzPnPM&Zqj|`d--DH
zr@cSoJ46fzi_1@+cDSEJtvg<S{$rW${<p4Q{(OGAeKq^{=joSRZho}>dD-=wUzZn8
zR_oXO^56ZHACXS;dw24Bd5AV$WcTN%D=e<L#Fej~bnjNDuTTF^Hz4ks8@T2YSH6C<
z-oC2+`@`*9u}owyw_JUYFPG^79}?!rxcmDG#f#PPak{l-|LW}J;gfEE{kl6z6hGT~
zb?@%zzn|UxcFXzMYTaMmoWDr#z&E^$+uc4t@{3nOj2(VWwg=pjz*P$5zoM`saXo_1
z9xz<<e78GXj$WIuF8ln~KfZY;xONQx@{H-6|LXcaiA_9M?Ed=lnu~PNFMr<kx0qjq
z>+ou+GngdpVinYp4Ne1pdw$RraVJ+1_XH8a1qF9XdQ(B};pW__E_t_75qD>ZQ8mh^
zA%Y8XkF(+KY(?B-l6Na!!#$RicZS(;_Z4xsv>oT2g1Zfd3UbF=aL1~+2dId<B<@M5
zG}^p3ZX(?4lJ~$BarY*8a2kEyxsqla(|ox5dbkT1AJA&#-HVZbiEyuMoUh1^6E+cS
z)Nl_{(te+lyoZXoJL7;rqt`_#s?X0y-VrKN>?9gyotqY+EP1{Kd8&&%RWrrYJ|t64
zbfz4ll=R)OosZ13R7BpIn3Q^z%#)Z<ki5)@e5{CkGUo(2sKPU@s2%U-L*BZ2$RikM
zkyJ9DjQ^6JJ3pU)d0wC*@@y=lu!ekL-(8#$d8&xKGw5_q9CW_qJa3_m_M?ywCX!f9
zcn*Ci>A$-b`8L`Qk|#8dz^g|4P?r3BB6KJfS#@NM0HR96i?O6VKPU1wR75`M=}f3f
zzJgdXN@r<CqditcKAM2wtt#!Vk)Hv1CgUPT75SW;Taf&JjQsXBjrJ^t9eB+=Pqr+1
zJ0tQHMQk`20PvzAFHkZ%Z*Swb=VsLnkx%k#9jel>NMcDJ-mJ)1<ioRU7}08?wB(Cs
z*f;Vck#`2Hl6-*4$b4<{d__@O8;nn{qRO6+V0_c?rasRLnk4cnd!7M`Mr>v@+SerW
zVoiFsKpOIYh?CLnTa<Q~yfa=fXri=CSk#`M6Zwkl`D~%FRn7CvXyQc3W1Hu%C(k+e
zmH|oQ$w@Ct$*gPZn%{mD@?wA?sAOKKXvC(GABVg*AxS)H$UAN_e?#(I4|xP50fMSY
zhUcO(@0)Hv3VAkcQB?DM;G!CyZ!$j$d1o-9sAN8W+M?duCiCNv-{#(`TA<mo6_vLl
zkIjcSL_X)<idE_MwroYE4Nc}pA#Y6@>Z3-tx8#b-d}vnPDCE7Nh)(0Zwd{)K+lMyV
zk3pVc+Y(Eo+gt03>h__<=SLw=#xX=y6FezPeohsYt|IbctT3zQc~`WSc2?x2BJ#l`
z?y3xbOk1|1QlQEF800OQq;Wb@wLC4_>sDq|MJ=`>@@%}*MQLq}C3AA-M82LH8_v16
z&Le6@dqJ3t%r|>}n0z!8AO;or*mV2pM{I<0w-i;BRv4OYKYf`OYeLfPRWfhOW?jct
zOFIgAZ-jLb8<Dc)n-6aU@;1RX=&dU2T3jQ4e^c_dBJyN{*Hu@E$H~O!o6HZB7h?%k
z)m0M8*44EhnB$NSh7eVq6C{?cvGGl}AA>xCiAb#R-a5dRdpmuZM>Z0OG*McI(;?5z
zlQT%(<*aK}CqzeBRJWfK`I<)iXaF&2;`5HFXce_@@%b^x=iJ+XLDTovSt^=)J16pg
z71{G-4LM}hJTEkv+*@w-^MmBQL5rjs?Hv`h=VwITQAOl)=LBL+e|zV-=sZ6s@)ha!
zNn7KHs#VKe;1;Ew0eLY1qUv()BDW~*^hbO4Mnu;M(L*DDQ}Pw*_7)6#QF(Ga+L9XH
zt@7}ihBr)}bN5q~_tqy%Ft!ZL$&(j<<E=wYtj;IhpkS;{a7OlgMe%tXj73IGPL5A+
zi=7Dh+I)DS9`b^@mR7Z<g_olB^K+8<nlwCN<AqS;!}HFS^xMyBp0CKBPx8x=R5do<
zyH-(oL-G|xX(bw90@Y|A++=k7=Di&tFZ+PTD}8uA7F~sWRx&RYk+;TzsAN8WK+(0R
zXC?ElBJ$oK5NNErfMGH+-^dS>kH)3JUNhPU3PtVtIgR$A9`Z=WT8x^k>%g(e{0+&+
zipYzBpd9UO(Ul5k<;kgNv`25eh^kvka79-qos;+0Q$6Gfa&IiwcydD7H7bP~`BBKT
zLDKcW458?%Wpg55(>zaNf?(8tWMGDfMKd#KM)SO{h<q@TWL{Hg5rLY_Pd!#A`;5ji
zsWLMnaw9+Wh)p7I4W)l7@>Fykn_Pd|GBZXX?{glWsvlnDqU+d*ol*DGn#e~J97dJQ
z7tY4At?%s^<SiK|y6UP}I2&tD>(W*frA_iHtkmtDD{9rviF`$#92*ROMFLH}ee}NQ
zXz$uu)MJp}<=*P%`4Ea)b#o$Lktauzg~kLb@-Y<6$#HF-ABDU#!4hktv<Of(Ogks^
z6}fUG=7f6)D)a;m{SEt{Rx{qa-1*i;-FT0{STyDSHt3-;^imOe$qiFLH{e^E40>$+
zt_RTvBS=Up^hq`3i5T^NP3i~H=XR!qtV8c?QQJ=1Tt5nWYpkHkhl}8)N&OAcQ$^@~
zZjt3tg+92+NPT7K+0}zS*=aTqSQYx{i>i5Z8t}0q^ql5;0@YlP5GI4ZwgJB;sdt6}
zNj29amMyoL6MC))easqu)Jc8WRSO&XanR>ZsE$>cS(s+5rs2<j0KMb^Ulp=N6e=36
zpOczbk#WbtkjJ3X@X<AI-c60WP!IYvc7x0P@igdt(LT7h%E*Lf+zq0)27p*6^}$XC
zJvQTR5WP?6bxkyKHdj&f+|ZALKALEiE*%0_)VQPO%^L%KQr08t%5IRd<G)G$<k5?9
zUW=Z}wo03m)I&vo=z|d<RRe6KWD)v?e%dY_pF3k(C!#X$C@Bask@h**mZcv=Zw-50
z`<%Qm#$-AtHM@Qoy*Csj=mvZbMGNb(sd;0d$0WeYf;*&ulj(?yEiH2pJsWhou31R2
z=)$G58t|zi^v-}KT}CbgPiBpqHrEfM4<?GT`i9XKb?LN!xqcM%geDS!DnW~BI%}FW
zAO0|UYXU2)Z<yU=nx?fZoniFecq^>Z@Y$AK=Z0F%-5BU|hNYmczMiFQI~8hb-YDqV
z2w>Dny$h2W@UbOm4Wf6(qHZliWFLz*PNQbmkAglLifDBjekdEh;pWX71HGl(KTGG$
z;~1NoH~j&B---cI7G7C|$;9hxOUAA!Xls+kh7_X;JyDzMr;lF1MCHQ8Xqk#e?wSvO
z9Q16WL(-wQJelRqnj?1@{kCPISE+ecY;);M-;uj*nTTFD*E^X^f>uk*96%qoEfYbv
z<PBOc&A6MoT_1AGL@R<SN3IR6`sSwSLq+IA?x3ZA>#AX6+3*dw0e=kih=y65O6mp5
zE?3VpQuCl9^mf}?L8azN;WkQhs;{>dq4$QYGH5FdMZI~p4ftcACoo8itU_;LGHV+)
zZ{9HaZILyp;`PZQXhqgDkb0N9YD67+M@5ah)>UI1^dU>V*A4idCnNQ3z#l}<Xo9dV
zUN6Ca7S-@u&~_a3);L5|b=8n)n>TOzj$HP}Gg(#8)`6Q$szy`u2GJ)w%{tU&={R(4
zz@K_<UP8}CC_8C1=xI777_E)cFnVV!<U~sidJb*CpS~m4#)O_oRe0rWoJ^%{%hDM{
zFJxHj!*>qaTt5T!Vu-YAo`;r^JN3e=kXt5yV%nxdA8;}aRhpVNfIjAyiGUb&c73Fl
zmN^6TYyv2~c^=9xXAG^E&M4?zMjusny+<xua|g}18wGtd9+i3cp30VBKnvQAfgUM0
zR8dMjmo3477PK7&z2t@}q?39pMI(36)Vxv92g8CvC-vfrrfM`bZxr;Dj5|`cR`AY$
zHAn8$T{<zhOa$$kw%+^6%=OfQwu9(%ku`%VeB*;}cKr;{M<dC)(@A~wld<cmsd+=_
znT%Mg3*Q6?&90vTddY=XphJ&E(dVUt+JHX-`e1-j*AGAB<o$`Dr{>KYLZ4T%1`Im%
z94CXmtzSQko=p(aMefpzQc}auohTZA`Vr7O<Jm@)OD8xi>C%}MdZ`F~%xl_$PU<}t
z9q@fC!59O*kZ~4+DoZDXqLsF{Lf>lchS5u|z6q*~+z^Xa+TLowZ@qMe(TBXItuAsG
zq3B{~v(oTu$<h(nu4E19RgT<<G8vZ+w7z+R=(j6b>(;r8R5oY}ZLS{$z0cKuH%j#E
zC$qX`o9hSB=asA(N#)Ip))%e5q1HEV6!g4Z#9U?EMTu>7%ZYpQGJ2N_uaxLrEGqS_
zQR+D8V=lbX+4V6%Yo0c7Z(c&bUCCM*xr-57XY3iEw<gW>ss&I8p!LYTDSAH^`VU(_
zZa7#x?#?eiJU#1|o5Sx;`X7IPIa!`A*Zto2+fCo^b^RW^^M86eI&y!$c){@9asT7N
zVt2ru-I+fA#{Bv6(`xtq!Gnv$@$BNmzE8XK_FSC3-E`}h-O0)JX{C>=_3_KD@6vg@
z-P@I@d^+yn{_FnL*~OLD<o@~j$?B-<(@(#AF)q*j)ydL#Kfhj`yuAE>xxc6PF&%p4
z;e6Qgr29|0WOf!@{rvTAt99D!<<oxK4VRDhi}N2mPd9z|+wx@DZI+2c{`BF|UAeD=
z$KCq)(Zj=Eo@O@NKX==IO>gnj*QYO*>COOt0NDM1u(<NY?#XW-)8U;h>eV_uhX)De
Y^5S878ZXkN5AXaB9RSb;oy33v0KH!FGXMYp

delta 3913
zcmV-P54P~3AD|xwABzYGaJG>KAQOt3bE#7o$2cO8)gl2K1(eO(s*yP(e;?nSp6tC@
zUTjup>-%@nz}>wc5APhjxad~v)%y71a@{X49(8^9N4Ht-rJdKC`*&aU{rUc4@%HUo
zbKG5=uGZ%0>~xX7+#J1Ho_3qX^@lIs(B0iB_PclAaFz#0{p!uCfB*Ble!0#!c;y7|
zHmlFh@RrTl#qojwaPj-o6(`TPr?|ITZ~AV1v`j~BR{Oj2J~@%_7XjFl2mvvFKR>(Z
z@8A7r|7v-0clTxb`ti#{0(AF4?h(8|-nY0<&O}EX;co!@00)b-)5ravbmx1|`|hHD
zu{vGezx!9mUS232`r@qbPV!;Tmq%yomzzUjf3Wy!>!$-At#)tF(`EXxn{=P~UcT7)
zY44Bt4iUq_;_}m{9Udf6>yDRy|5&EG|F!FvKb@a$U(NpgdHN-nn;)-#T6X>Bm*vIN
z)%tb6{C9ulN2JsI)}6dw9->Vb+5P$H3X5wlapmi$-MiK4>(l?!4T!tu2Clipm9HPK
zx36md{&4$NEECzwEmt4p%Vm1NhlKet?*6_)@nUs+oNjH|zdCz)__SMpzwS;F#m}~0
z-Mc&b?`L<v-Ew}mTK88s=P%Mb@D1<ccDK)u{Nj}mV~1ap?E$wWaFqi2uPE$DT#w+h
z2MpId-|Y^Uqu1uE%f9&K_ivsFt{ua_JYzcNzq-CpViOM*yT87?<|1A6%b$1sE#?>D
zI=ouy3?@mtSOs-tgVVr&o*#5Y+{snMJwZfpLBXAp-c*o#xH)&KOWv(i#N8QURE_d!
zh~R?U<7~J)TM_q|<lTzbaE~SBonbcIeMQ_YZO3`1;BJGVg52>I+_5U|0V?7yiF*<%
zjW+L%n+W&1<UMdj+`S1NoJOB_uA~{qG#~E19_|9h2ecY__hMv!BHSw*=PR<~giQn+
zHQa-gwBP3>@1Y{@&Nv{@=yg$w>htrFcZ7-*JBfx_=cYv{OP+5*p6Vh`)lBiU56P4h
zohgSXC4Dz+=Ognh6_K|lCZ%2_^CTt|Brh`}A1flC%sD|0s_={}YR9|zkhiWL@(9LR
zB$do3<G-Zm&d+Cmo)@TyJR8d>tRY|6cNb?wo+={m3_6_?2c0iD&s%7t{V3#vi6mAN
zo<ko>`tNQ<zK!;S<Oz)<@T$>1lqEl(2pvjARvlR*fT+^&Vk~LT&xyPZ6_HPRIuojr
zuOOC;(pj3(Xpa?<k0u~^t4h0T<Yz#h$+(D7MLs9z79>A^Bfou3qdkjZ2VOJJlPyc$
z&WL<P5gQH$0K90(3zUq`+uQi<xmk5X<dghbhpIFzl33D*H!JcL`S2_oMzoqJE%~Au
z_Ko~V<edSlBp+ZhGGE&~Us06S2IJGKsIuoH7~eF!sn7F*CW*Yto@aof5t|u}_BF}8
zSd*SDkcK>e;$(FD7Ns2~?~E4=nkX$37PaT+M7|<>K3iyPRr5SEnm7^i*yj1`$#c%V
zWkAw+a?*=ZGV9v9=C>b(ycl2zDw!858nJ2Q$06@cND_}4@{XI#-;g}lLmt6MfS_uU
z;kl^H`=;BELY@s<6xBQ*xTuEbo6L_w-WiN2Dw&Uewy5{E$^1Cvx4E~f7HGC?Mdhu?
zWAot+k<YoeVpY1mEn87(LzDSY$Xk<!`l!+EExDpHADUG+3VClRqSJV9ExV%m_MwgT
zV~}Uqw#3rt_SU+hx_xNz`BBJ|aSTz_1W(G6pHoGptBAZ9E6l2S-W9E-ofUbhh<q@K
zyD9^J)0VBM6lgL(26>AnX`GH!El-Q~x|JDKQH!mJJR9$HQCb^g$()=yk*}x5hI8(%
z^N5<!UJxcD^Ua<gCLaw2h(Sd@Hr;;u5gVc0EkzZj6^5qUPhaN6nvisRmCW0+S=X`E
z(vCvj8)03<Mx-qH=EEC-yiKqTdaKI17T3ss-;}(qh&-9#b=6hkaWe7wCiBDO#aKdB
zb(Ms&b#<)=<~ZboAw*T@1c_y9Y<$!0#~_bjA`)x7w+^u7-cDcUk&Of*O_bK*bjWk_
z<P4H`IqO>03DFT2)$Qj*zNXPW8bA!1_`G8(T1D+!e0~h_IrlbT(Dc1^mWt-y&WU_~
zMfN;dLk?Lr&kIc^_m*4z{2+O6&?2cudq+j>`5BRSR1x{yIe}Qy-`;sHI?vCEd_}r_
z($+YlYSl6qxJ7AaKwb=hsJfiH$Sq1c{n4Ji5z%!*^w7xPlzc_Hy#>QwRGu7<wxou4
zt314>;SH1L-2GJLz4gfwj4cCm^5n&Tc<WFTtMf@WC>W~~oRK|WQGDJ8W06slljGCd
zVkbhrHXmN7hrD2}rB$tI;iYK({G4RICJj&6cp=pI@Vs*+{r0n(=PR=3ll*cdRgI1J
zu2odtkbFf^T8RdjKsDM2HyPc&d2a{E%RZp-N*|t&MOPu8mCQ>;<gKwFDwz*|P;@Qo
zS;@Ssh`cul1RAR@V3>@|H}b>eqj71l*NpaoLQ#8uPNRLOhdh$87NaKXI&f?<e?#)I
zBJyG&C`Wr+bfv;sd2%Wm?a>=AqUx3sT+!7@=j6ThR1bNA+#8EEo}5s2jY^?LeiZU-
zkaRsTLnyjx*__DNG|!WmAQ&}&8JHnr(aa2*(LC=fA|H$-nb%ZWM4%?~Q;*fjKBKWr
zs?3at+{jNoVw1>QL+PK2JQZEXCfA>~%#0Dp`<#cT>W3G(=sGrHXVm?)Ci2k)hfyW-
zg|o42>w7x}c}vELuDU7~&c>S4y0jHVX_Nd4D|LJ4iduDZB43dw#|A@xkwBAgAH6R+
z+Pk(E^%&%LxwpD`K7^uH-JHl*<jIj_p)rAqd<;c%a$K9|M<MS_u*8}uEdrDc)6NNf
zMXnr)IpH3H3OzwXf5ZN#)r|KpcfNH|H{K&K7EQUo4SJ{yy;Ou=a>EqR4fvKOgC1MI
z>p}Fv2ojPCeNqj1B1V0GllnpQxt%E?>(Dz})V7m0*N=kU8Y`&s;UaixQh!7AR1tcg
zTV#1up$~2{QePQ*cJ-i7cA5<YR)s$LqH5lp27IgtJ*T;zKsDDRgvp?<ZNRTd>YZUg
zQqA><Wy@{mgq|xxAG3xZby8n;)xw5;9Q3&ps$*4V7N!}iY4|gLKreZ~SA}d5g^EV&
z=cMLUWZZEu<T0o;e00s5cT?jo)Pp{a-QaS6JPmqZv=8pBGBTkVcZ2Ay0U*{%eXx^3
zkIlFnMDG)NT@#I*%~cdVH}s>Rk0x5BONYP}HSVZ+^Tt4*l=X<ZvKyrA_-|4_dGunO
z*P^Ggt<vTs^-vLi`d~yz)c_kQS%ki!pSDZK=gyeciKvV_N(zEZq<s#yW$6deTf<)0
zJ|{1XF`3Ru&8{Ct?+pbBx&hxq(ZYIcYTg*=F$u7;;0`I^WIEzvOUoQY&jy{YYZg*0
zx^U^N27Iary)$4*myyfBlUd`Y&Gp0RgNdT7zG1XQT{>-lt{(+Gp@~GGO3-4O&YEV;
zhd+$on!w8H8)i3|rfDrpXBfRV-U_QUe70rRxuI5bHwOBgVJWDquV*RSPKBD9Hwt<-
z0vL5t@4{pTd~69?gXo>Hs9Vbr*~g-d)2P|?qo9w5B3hk>AIgSrxOwx&KyNAc&(eAG
zIL4;tO@F|Dw_-q)g;y3~GV%J_lCdia+S;VCA;qXdPt@l6>7y4gQMqt2TBf3ryXM0m
z2R)nUkaXxRPiA?u=ExmJzipZ5Rcf9U+gv)+cjRtcCZgBP^-d;}pw-ec2hfLY%S6yE
zd4twVGw!Bt*N5CP(TbqTk!wS%zPTy-P!al&J80>Dx@y>1HhjZvz#juWqG1-Nl6rx%
z%hmIY)I6vNz1_A}P^o!RxQ)`B>g#Pq=)EDU4B84qQE#4Y1O6E32@Db=tI%7R%-V*{
zn>UPpTVzeDczv=6T9Ne(q~7JO8c~PdQBmWrb=4RLeaKSpbpyWV$w+-0@CVT|njox;
z*GsT}MK$~uv>gY%H4YI~T{R@y=FOYFBbUALOjZ@Nb>JqGs?pTELG;N^vkrAxIu2bM
z@TZ=em(a5j%1#;$dYVoNMr)%qjNTawInh#so<ke(r|-zMF`*|?6<#?TCsS$LvUCQ~
z3mMk>@STG;*UtdG7$U8j=b>ffPQCCd<d%tln6~NA2b@epm8RwmppUs_A|OVcT_35X
zWzGOSn*d60o`<r_8AI!(GYWc_(MOeC?~#kv+(9$$MnNBqM`a$qr?Mp&(1NyOphwCL
zRg_ZCWlJ!i1#L$`FS(%#>7?FD(a0S%HE$I3!LVS^NxitDsTxhq8wEWj<BpWA6})qQ
z&5=8Gmrl$r6G6MCt@nO1b3L`7?I8MGWX+%o-}vC0T|Wc#(MYoHbW$JvWbArsYTgif
zCL`AB!Z!gzv+HMoUUK0T=+I+P^m(bEHsFtdJ{Vxs^}`Q2d4D44sd@8;(C3w`0fP=b
z$H|~?>(>vXXA^{Uk-PMwl+^HZCyIuDegyQ+c(zgH(g_Yrx^!lRUMfN#^P0AxlX{Ou
z2Yla3FvdVHWSqsI%F+p;Xr=9~(6^eqVf2!#Z-Oc#H^icqwznGaTQ8kq^dYZltBc%4
zD7x6$tTg;uvUCKtD_H}2l_NKzOva@Ht#95S`t3^Ax^?a%l?~cLo9jnG?{oEkjS@Zk
z$*gYK=K4YOc_nK`QhD>D^+l_1sP)Yo1wC&UF;^LPQDR%&a^l{+jNawKD<ygti%NZK
zlsXRjm<z9Tc6|)cnx{?No0rgUSF%<{?qbB&8G8oktw}SzY5^1iXgzXoir$Zf{=?Rf
z8x9svy7S8q&(8Yg=J30d{)eA`PnM_4b-(xhcGLHJUB3tK{GWL>^LVth0^c3?KO8J}
z2i)16>EmzApD#bHcHbX7yht3+E<Wu0tV?gt#o60Uw|?23oNS*~`nXyjzwG)howwV)
zU5U!4;|?CY?q8i<TzO3%oS&bpj=Dbm^vf6H^4wpYEPeOW>($B2%b)vyTY4YUp;sQx
zhb>RL|D;Q1XTjCaU+=bBr_EkI?YG@<`Dnj5|G|rN(}%w<PnO+gnK<N6A0FM6`#N~i
zt&bl+I{f)rX0!csxBZv&7C(D^`o}We8NfY&-Tw!RD_`uM{Pr;&-r1sFt<!UOm{2Y+
X9;K)8N4oUko&N!PA4q)!#DD<+*L(8R