From eccb20b82fb9d4415fd6507d077363c328c3a683 Mon Sep 17 00:00:00 2001
From: Sander Roosendaal <sanderroosendaal@proton.me>
Date: Wed, 10 Apr 2024 22:17:07 +0200
Subject: [PATCH] fixed dataplep

---
 rowers/dataroutines.py                |  33 ++++++++++++++++----------
 rowers/tests/testdata/testdata.tcx.gz | Bin 4002 -> 4002 bytes
 2 files changed, 21 insertions(+), 12 deletions(-)

diff --git a/rowers/dataroutines.py b/rowers/dataroutines.py
index 38f2ccce..aa562ba2 100644
--- a/rowers/dataroutines.py
+++ b/rowers/dataroutines.py
@@ -569,6 +569,7 @@ def clean_df_stats(datadf, workstrokesonly=True, ignorehr=True,
             except KeyError:
                 pass
 
+
     if not ignoreadvanced:
         try:
             mask = datadf['rhythm'] < 0
@@ -1418,7 +1419,7 @@ def getrowdata_pl(id=0, doclean=False, convertnewtons=True,
                   checkefficiency=True, for_chart=False):
     data = read_df_sql(id,polars=True)
     try:
-        data = data.with_columns(pl.col(data['time'].diff()).alias("deltat")) # data['time'].diff()
+        data = data.with_columns((pl.col('time').diff()).alias("deltat")) # data['time'].diff()
     except KeyError:  # pragma: no cover
         data = pl.DataFrame()
 
@@ -1781,7 +1782,7 @@ def read_df_sql(id, polars=False):
     if polars:
         try:
             f = 'media/strokedata_{id}.parquet.gz'.format(id=id)
-            df = pd.read_parquet(f)
+            df = pl.read_parquet(f)
         except (IsADirectoryError, FileNotFoundError, OSError, ArrowInvalid, IndexError):  # pragma: no cover
             rowdata, row = getrowdata(id=id)
             try:
@@ -2018,11 +2019,15 @@ def dataplep(rowdatadf, id=0, inboard=0.88, forceunit='lbs', bands=True, barchar
     if 'driveenergy' not in df.columns:
         if forceunit == 'lbs':
             df = df.with_columns((pl.col(" DriveLength (meters)") * pl.col(" AverageDriveForce (lbs)") * lbstoN).alias("driveenergy"))
-            df = df.with_columns((pl.col(" AverageDriveForce (lbs)") * lbstoN).alias(" AverageDriveForce (lbs)"))
-            df = df.with_columns((pl.col(" PeakDriveForce (lbs)") * lbstoN).alias(" PeakDriveForce (lbs)"))
         else:
             df = df.with_columns((pl.col(" DriveLength (meters)") * pl.col(" AverageDriveForce (lbs)")).alias("driveenergy"))
 
+
+    if forceunit == 'lbs':
+        df = df.with_columns((pl.col(" AverageDriveForce (lbs)") * lbstoN).alias(" AverageDriveForce (lbs)"))
+        df = df.with_columns((pl.col(" PeakDriveForce (lbs)") * lbstoN).alias(" PeakDriveForce (lbs)"))
+
+            
     if df["driveenergy"].mean() == 0 and df["driveenergy"].std() == 0:
         df = df.with_columns((0.0*pl.col("driveenergy")+100).alias("driveenergy"))
 
@@ -2067,14 +2072,14 @@ def dataplep(rowdatadf, id=0, inboard=0.88, forceunit='lbs', bands=True, barchar
 
     if 'wash' not in df.columns:
         data = data.with_columns(
-            wash = pl.lit("0"),
-            catch = pl.lit("0"),
-            peakforceangle = pl.lit("0"),
-            finish = pl.lit("0"),
-            slip = pl.lit("0"),
-            totalangle = pl.lit("0"),
-            effectiveangle = pl.lit("0"),
-            efficiency = pl.lit("0"),
+            wash = pl.lit(0.0),
+            catch = pl.lit(0.0),
+            peakforceangle = pl.lit(0.0),
+            finish = pl.lit(0.0),
+            slip = pl.lit(0.0),
+            totalangle = pl.lit(0.0),
+            effectiveangle = pl.lit(0.0),
+            efficiency = pl.lit(0.0),
         )
     else:
         wash = df['wash']
@@ -2145,6 +2150,10 @@ def dataplep(rowdatadf, id=0, inboard=0.88, forceunit='lbs', bands=True, barchar
         data = data.with_columns(
             workoutid = pl.lit(id)
         )
+        # cast data
+        for k, v in dtypes.items():
+            if v == 'int':
+                data = data.cast({k: pl.Int64})
         filename = 'media/strokedata_{id}.parquet.gz'.format(id=id)
         try:
             data.write_parquet(filename, compression='gzip')
diff --git a/rowers/tests/testdata/testdata.tcx.gz b/rowers/tests/testdata/testdata.tcx.gz
index 62b2f5931143a4b2872cff35338cfaf73f70a23d..4623d9967757276973f7be061272a75a8eb56885 100644
GIT binary patch
delta 3803
zcmV<14kYoSAEF-yABzYGb@7n~D-1S3(k<Gc+qcEkRvnQ&HWC+`Zna*m&mS+>{c`i9
z>$_jOk%d(jo<I9(cZx4o>uukyPnYSa?doWM-e-~UB!64ipYHcx-dqa*+NA~o-S^P_
zZ#+DCbZ~I|?g6Bqkw3oQ@cY3}*H@eVtHXcxFPEFc{g>(Mr)MVw=)M5?g5U-6k;Nl&
zCOYB>e*-uII9{Zk{@DLncYW|v-);IAtIOqAhkthL&4tpTFRuFTA|LkC^7LwbwmlK{
z$BR$4z6d|y$!h-wJztZX0V99x`sKIRm%CT9|NT1slAFy>*WWI?e*436^L({_)i3|u
z-}({hG(UG2ua+lh(?#}wzQ4laj!WG7`g!+ub@}S@|8xW5?zn+FE^+JYr|aFT+W&sC
zdn=ZS?9G<j5Atl89`FfaevJFSZ&BQ=&d<}WE&G>OXD83Q^{eh8QT%^k>+QX})Bk>O
z_q#3ESF3e@dvpFGy#s&YUEJ;V;gKJ`5@PJ|YqC4wfdp<-ApZ%4J&C&!eDHwbj_14E
z;b!#Oesb9tKm78SXM#J&@Q=@!&iSvd?~~ZX<Hi23kFU8%7ybC>-F=JsMYs#EmO6t;
z(k@m(9ogVCu;&L|5qE!b6>(1x5nNDkr=&L(<R0$No$8WzD;05fh8R_&d>SISAon;M
z?#@=kJtld#;x*i3NqJ|O4R>D=cT2l*-YK}-V5lH>d;oW>ihF>HxJ%-mgi52$d*dd;
zy)JnVToHF~f(NJ3=bbBQ#xc!@yRV14fbjvXM&7*`nF#mF#`%AW>^NZ)!A1@DASLbh
zImvsdh`Tcm2sC<Kl%o3leB>RWBE?RkVb-~65z3P12au<_$Wt{_JRLzY<wR%7AxcT#
z4cqz1JWEC7t%*sgSIIny2?fc^jL64|$R~47kb^2b<BHnxZa(C#tA{*-aTZA>^U3%x
z>ACasndb#6BF}%uG74+R7xvx78Ih-o$UB2h=fpwhOV0Ba+Gsxt`CuZ6)r9BJhm!ug
z2a#{1{UCWl;|RQJv=3#;&nH5MQjt|h)(9Y~G`tu~+VgWFZ$m}olb+6ms^lw(C8Knf
zW;EJkMdYIi2;QpF?i%?SkY_S3VpNgO$+-o|&&Y3I(`bLsV%UM#%=2W+lD9J=Us1${
zg8=|98u9`qqx1G5etT|K-4OXCzt*8D4T~g}^x@5ld__Jy%Z3rHCQ3`bXoh_wKN5Lo
zz$(cHn2gNVHqTcSrM1EM^eU?C`3S~04R7l6yr4-Uud?SEplHNqMx%XAGB4JoXA7hu
zk2o3KzD0j&hsitR1%oC^%Y;Sk`8koV$ezy@8e7#o&x|Hcggmx+{%-P|b8i`tG@hLF
zqLj?Kwyyc@M<Fi;7=lXXg^EUO8u@X^dlQnxqlUcWCiC|s&-IW;FcKiBnq+t`D)YYS
z_M?zz!xlv~&j&85;rS->qmXw7BZ^Aqqb=&aZ8Cp94*6Z~t*Ql@En89fAoAFJcthlK
z?yXprZg0z0RNByFeiZW7q@g})bbCv#sLY3E)r~^l8;a;O-doGAXuf@Dqx~4<8FnqP
zG`hXDuBdJwT6}&K@?;!CR5ihqvgGGfQRyloFUAV9YMysRYiVaiUMeCVOyaJ}z_ev6
zDg}R<%#T6dqDdO3BUQ`OqP=cqMpe{eD<aRvJ6)94##k~ZXHMknsj=ald+R)+X0#WC
z$;f=O=ZDEhLjht?k&jKcpMJzfD0fRyMQMeh>Gso?d9fxW-Cia0wrtjQY_+tbkoQJd
z7qJm3OTPK=Mj&q!Y=hpavaZE7^7kcgD<XeSCU{+Sm3W*?e7?#2FnKYSP*q(ep=@1U
z>w!5A`Ctf9)j2_8*%}+)bo(*LBbbQ98t<(GY`M47mw99(fk+djbvPaJ+&nph<Xz6X
zR&_#jghh4xIgzhvw2uZ5gC;)jn2J_W`xc)cgM7}t4Hz_iZ=I#0xwms7Uy(ge){uWg
zR?YK5lgYj1RzE*T-W#+?s?pw2QG0$y<Q-K+K6g$a*7Ub`o{P@&b0S}nZlAO@j;LC-
z%mr>y+8K}+10bp{=Pq)K(oTQ0XKzGwoe(`V^7kcQk#29nuosml$D=K&;XNo1uW5M0
z<T-agRe5iHvIJwxz??jJG2S}V#Oi;1(hUm6>I7$G&sP+mx4~Fs)a2y&^tRZEkgv^$
z7wRD|m}_ZOYg%|IT0cK0nXgI16E<E5H9kD=TuHzEtmgTO?D-_W97$DU<GpJYmG>lH
zQIuAq0VYt5_Q6d?w{PCt0rGMLXuQ&g=VQ@T$Y&+<QW1G;EQm_x0~B40dRBij?<yki
z4FZA2stXt<BlC^?F!^X)8tgTreV|a(o}bfbAL=2GWUR%g$+`|4o6O&ne5{DP7zoPI
z-WFY{a8{n2ibi|%#*3)Br36=Wb<#O`Z#~sRo*?(eVvQ#!lwG4zsF5FqJR2lk56lpX
zu39!H@-@x#Bqj((O$KI&STuh#gJv|(`-;d1BT42pl@<}G$^6t~b#laLER!lTBO*8Q
zQ;*mr^43uLry@^9*Rjd<r!6yM1oA%T;i>xJMJ~FIjo2A=Kdp&;G{Iq1$$a5#EZh3t
zjzQj%aiXiPiiNYW=Cm$tMN!%$zrsr0-npVy-JHl*<jJwYP$bag+ed%zi;niLtwlWs
z`F-xKZk`XJs8u&7@)db<Bw1)opdueb(VQID=J`>`I}<FiCQ6F{Wy7>{LSK<9M`BL6
zhoC}F(9qwr|7kVjy~~|%UDS>D2#iHj?jM34Dnl<7p_kk+1#|<xrOBYj*6(@{eK3NA
zq(YxmL!O9H-=uyJeQtkeO2|6&&K9-pq|NoCptr^fs(iQzUYgY36FpUg-scus9#!ar
zn~c;~hMrwL=#!mh1A$eckG`myH>Uv~D?-m{t|w5<^$1}y=xZDBYm$0r7?4zRJ!08%
zn>nH9iqOZb;YXd+mtD26p&ti*?u6=Cm6?TU#%db=4A4s+@Kt{yTSTFv(fT>5c@-IV
z91M93Dh(f9^XA>xxC`~5Ph&T@+#gSa-WTnI`=E?WXvW<jdTRiPby6SfWYA+X?gr8O
zgkIM~BWH6JMb8cWDCnb!R_W3qa7B$fYTmpt&?jX*qOR-)DLej~)K4D080WR<scfsX
zIY~WKggzJ%QZ;|TMoJc;Z|JA((($=7rgb7J<BpPoAQNeygKb&*LG;$J*R{{d3u8>C
zb5gVGhtYdOL4t0;_fWL39-Eps26{{atSq=g3OJdLxY*J%2hp=Zr|X)96pJogI;#Pn
zDnjoJSkh(WGVo;9xM_3!F#2GkD64N6ZBds_o9jnGPiTK45vUThn5MI)S@Yo!qqioo
zvigSEO{QsD%hDM}?~S*@Dh;1)*>!HH)!dDNK4(}8>gwxR%C=LXrsj=;o{a!Toz%N9
znE@YLg4Q5<XDsU0GDP;VXyY_$cKs;mqoIgar{RaP;Tvw=yfM&Q%Kfu+-aL-6sd>{M
z@U0jSW#NC7MVL&yzP4oSih{N_X>3R_s?Za)xqkZS1x!>fT#S~fXymT>@W(;VCORY?
zddrhp-mE!thtcm^CVG{cXT>&`&h#C*yOxRQb#uLw$s}mCw9EnYVb?MdbW7f#_0o*H
zsoV7-w@kDmsB+}m(5i3li#}9@KI9Hsx~>{FmJNU3a2xQ)K#ypc#i^uTpzLz>JR>y^
zDnf5}trb*io)m7QG^hG{TM>G1$SQ-j!cf$kXWM{326_U6#K<c27ACW{Ve{q<qu&)-
zlPX@HEP_^KJp-wCxvNIhp?6f&xNBWC#z7yl)O+24?|Cv(-v<0a^o%A5>*DniY*7ur
z1#N%FL2r#iL{(P}iMDz3rtipQZ#<J#1#KO;$)svDHE$4oveT?XU6zhR*9QEl=jJ8!
zY=p9tMuVQFQ-aajC=H``#zIcC)S&0k2K?zea&1iLiByGG&c?}9+O{m6LG(g~wLX03
zpw0C&Kre<!tLAxV8M#w0yb8HxBBpIR^Z|b-(@>?Uc?0NUZkY&(QD@gjYH67>K+h(C
z(wpa@>~hA?dg+XU-evSrW!HP;qBVEWjJr|LN8?eMhwrIu2?n&F?HK5hazhoR)N|Pq
z3}`{yQP4|ns6slaw^B562Tjcz1${6q7<5uEu4t-8Q}aeaPsz9=WorfRTyx}3-KBpM
zbIU}~u4(JNpUhlOEoeK4J{MUtsKPfs_-5D70DUx)tUI05M?V?6o|>9Bgr3QWwYu<4
zfY9vv8K9S3cm+E2SQLF;DyR+kBcKlk7<K*dLr&hG2zqMXydm^?C2PQ-L(g$C=-c}B
z!|2%rAzkDyy(lF${M?D6p&tRgGoF8KRJnA5!;&tYS)rGT(8s){E$F1)W6=TMw-Sso
z&<h!7F{rY1LMU2k`ylkK=583h<m#KC%E%3|Xr=9g2K?4bXBd6RYuf4}cM*y%b~Y;w
zzm_Z=f!#{hfL`UujVP0G=|Jn7H;8_>lC^G~yGUh&w$SGKQPBHbeWOIrelmZnTei7=
z5Pe?BnvqoAyl8#V>Kkf(^F~3>yG6`Z#$A-yR=1qEH!q`ix$sJf-o>I)-x{TkgFfcM
zE1g{*1GMI86ZhsN^t+X;m65v`v316t0eWlFOs`r1g#cQQ-20;UW1)Yy^^Xn5i)Y>S
z&4=%=`sMcIvy1-A?=F^?%XNQ$@cC}j=LcPX0B`+o(fL!?ot^y(pPlz#9xwI>Jldb>
zkKdR--+Wr_zdwGwNgS^>@Amz^OK;ES>dm%WpLG`(yQh^tuGZ&gUEihicKf$0QF(vd
z@i(vfmsguxugN#p*B7hPu1`Px=Eb-<_a`Sy-+lXPb#Zp{bAL|nV>*BI*28(X<$3p?
zbjj>2xc&Ll-B#<g*~|O=b{lRU?MLT7evxka<j3X3vfC~bhy3Z?qq}uq$IrU;`O_yS
z-+iCi?Ec(u{~^7_-@m&2b(!uA;0u8L|Hq44U+ka!?lGM_+M!;p({p&7P%bx5($n}g
RUHatF{{Y}vM=ivF0RR_n!D#>h

delta 3803
zcmV<14kYoSAEF-yABzYGoad1SD-5<k(k<Gc+qcEkRveK%HWKHXZna*m&mS+>{c`i9
z>$~5&k%d(jpFRC~cZx4o>uukyPnYSa?doWM-lvi9B!8moPxt#TZ!U#@?NWn)?tAF|
zHy)lmIyg9f_W;t*$RFQt`2FDL>#I%w)!{$;m&?uJ{>$|B_h%;r=)M5?g5Wvwk;Nl&
zCOYB>e*-uII9{Zk{@DL%cYW}4-);KmtIOqAhkthL&4tpT&#(IKA|Lki^7LwbwmlK{
z$BR$4z6d|y$!h-wJzJBT0V99v`sKIRm%CT9|NT1slAF!nufJV({r1P@=Gkifs$c%Q
zzx5;1X@2Q0UM)}1ri<+Ve1C<-9hbQE^|S8n>hjg)|LF$A-EjkVT;kT(->-MCYXAGm
z?yXoRvNv09KghFXdcY@y`7!SQzD04fIzLahw(MVCot-@E)~~vYMDc%vt+)5?PXGJC
z-S4(sU#-^t?ald%^bY)mcX7Adhev+&N{F$;ugUI!2NJkVf&3>F_9X5`@WBIyJD%@u
zhnvxB`^jaW|M=@)o(b+8!#_S_I_JN-zE5Hkj~DyDKECE6UG(FhclRyk7vV0vTIviY
zNxN7Db!3Cnz@8s;McjYMRm43(L~udDos!;EkbAg4cdAR?tyIL_8Ddn8@@a_Rg52Y5
zxI0@B_n73}iq~+DCFPxAHr#zh+%4_Kd8goRgQ0@l@d4biD((R);x37M5-N>0?~R)X
z_qyaga7Em`2_BqApLedL8OJmq?!F%G0>%fl8hQ6(WFp)v8|Qy3vg3qJ1RFKngOs%2
z=Opi;BJR#OAkgS_QHtvG^O1LiiWEDEhFRyPMJP+2A3&b!B2U#!@pJ^qloOpPhbSd|
zH*Du4^DGsSw<ac~UM2G+CKMzuGa?@=BA?7TK@O_$j4Nu#yZMl}t{(CT##tnl%qQc&
zr034hXPy_Rh&+EA%P6cNU)XmSXGESVBJT`3of8M0FFDU!Xruio<b#PMRui5>A4>Y~
z9z?#4_JiaJjU(`?(LR(VKc5I4N<~&3StEd`((qy|Y0uAzybTqRPkK5Ns*<lDmW<L_
zn$c*F6_Jl7Ab6`vyKCfUK%U9Ch*3p8C+8L<KO?_=O{0H3i(v;|GtZMPOWw|id_@r(
z4h8_cXvhncjLzGK`0cq_bwlKn{91>qG%S)>(uX%I@)h~;EE`6&nkX&#q8av${7B@T
z0jne*U@|ga+dN-Ul-360)2pbm=OY;3G`y+L^MWRcyvm+ufT9tb8IATe$-G#To-L4u
zJmO?@`xbwt9VYLL7Yv#xEfW^C=jTMeB6~htXlzyUJTsa&5%SpP`Mb$;&b?(o(s**x
zi&8S{+PdbqABDUaU<fLi7b+UDY2?Qt?@dS&j~eoho6O&nJl8`W!AO9hYLel(sLcDO
z+mAw?4O<k|JRi8IhUc5ik3!xVj3_FZkG81yw#k3|IOKP^x2hIswroY^gUDm^;SG_`
zxwm3fy1gx1QE5Yy`BBJQlZN`J(d{j{qB0+vRW}NGZz!VEcyBGcqWSirjrL=ZXV|sG
z(&+Zqx}v&$Xz}?`$dhplQPl)b%95W`MWw5VycjFYs(IcOt)-n6d8vqeFp0Y=1JjnR
zs1$!_GCu}+izaEDj#MpAi}t#e8C6k>t%y7u?{raG8)M0woH>!Nr^beJ?yd8Pn$cbm
zCL{CBo*yP34F!lnMLss&e)<s`q1-J+6{Qu1rrS?n=Ea(jbbFP|+p<~LvDMO!Lf#u;
zUBpJDEcxcc8-ct{unl^v%DNWU$lsT|t%!d-nc#KRRpN0n@%bk6!{o(SLREE@gtB#Y
ztq0~f<bxqZRp$hWWovAF)9uF~k6<DaYrMA(u;t!PU*?gG1R_n8*5P!>bMxd3l6N`l
zTGa{B5f;_$=S04y(LNeL44U}7V=7uj?OS|)4DvbmHek^7y>*s~=HAYUd`0#=Swnvg
zSvAiKO(yr2TmAeXd2i4nsYZK8MeX?+k#|%P`P?~ySkvF$c`iE7&xw3Rx_#2tIHGFR
zG8ec-X=gxQ41lP*oV&;^N<010p1l##bwc#e$lsTIMY_EO!(LRL9FMl7hWDU6yr$s|
zljq$1ROP+($r6k$19S4^#dzyb6RUsoNjE4Ms}r1&Jzr6L-UeflQInJ7)7xSvLcTU1
zUZ{t>V6LTAt!d$<X#M=0WWFX1PuO@N)cEkcb0z)uvzq5CvgechawJuajrXoqRNj+(
zMNwLb2ADuK+6Olo-M)Em2gu72pz%r{o{vRWA)l4ZOGV_Zu^=j$4^VV1>REruysL=3
zHwXk8t1e)ejLbLk!{noJX|UIf_JKlCdwx!%eW-^#lCc(}ChIzIY%+gO^06ZFVjw6-
zds}p+!dZE8DjMz48!w{jmJ(di)k){%z4cTNd4k*<i#48{P<D+<p+<fb@@$ZFJupKk
zx@y^+$k#N_lb9eFH5r&8V$pxh44TnA?<*o7j3k-YR9Zx!Ci7E|)yWZ~u}rGWjELOG
zPd#Fj$Xi3{pNc#cUB@QZpSH}55y<<Tho|a?7rE#<HezSg{j?_X(FBK4CG&-|v25#m
zI|g}6#)+=FDi+Shn$x<p6-8;2{0b{|d*_N;b#o$LktfFnLy<s}Zy$fXFFM-0wifjm
z<oCI^x_Lf?qE_9U$XDdakz}DUfr@+#MRRgoo99O%?@X}7nkX#-lnv9)34KMb9EmyM
z9)b!zK|_Df{-@Q9_bzw7bx}9oBQO?Cxqk?Hs0_VSgkEyP6wnR$mL`K9Tfgf;^uY)c
zk_vrN4S6C)eUth@^tpeXDIx37J6qJYlQ!3ng5DY{sPf?=cxh69PxMp~dY@Zlc~qef
zZZc9|8G3g0pig$14Fpz&KKi0+-kb(}tOz})xt>5Z*CT|<ps#JfuSx2iVL(#N^@wH5
zZRUiYD?%T$h97lOUv|~PhJGCMxf7~mRc02Z8LMgdGe9qSz*m2TY!QWuM(gLK=2c|e
zaWLdDs5E?Z&6{^$<1W;LK8@Yra(_GxdSA2;?t?Nip&56B=&b=D)=7P^lR=NoxEn<8
z6M9_}jhxL@6g@Zeqo9u_TBS>ez!f#_sCo0oK%bQLh`O>Hr0n=_Qa^e0Vw~5ar?Rcm
z<|Oq{5&B?6NY#G;8!1_YzM-GCOULKVnAVA?j5|sSf=r}+4z^|K2hm%@Ue`V+FN`sn
z&PmO#A4cyD1qr$V-$T*DdTeUm80awxu(IF|Dd1!};$ln797N9sovv#ZQY^Y~>8u8P
zstCO^U`dye%fORa<EG8^!{~#FqO87Qv_)MyZLS{$J)wVzM4(F0Vw%pHX3d8`jNY2S
z%IX_tH<_ksElXz@y*J(pt2BJJW!Jf(R&zH7`kY}YsH?AMDceqknwmEXdNu+WbyDxb
zWCna}30i~bow2A}%MjVeqK(t2+4ZBKkA@;zorWLEhHtoe^Tt4LDfiFPdGk2Nrshq5
zz_(&Rl!bp+7GW~+`r4AQD+=1$q_H8zs6tQF=KATQ7cfz|a4}k@qLI7i!ygAdo9K{q
z=q*oXd9&uo9Y()vndnt&o)z0%I@5RL?ph|I*Uj}#CX=Am(lQ6ohh57=&@Fj`)=M+)
zrf%1V+%nOMpvsYJL#w{IFZxgs`j9(l>AGszST=ur!)?GH13jW)7N?SWfwIfh^NiFy
zs0h8?wN_B6c~ZEI(wyq+ZAIw4A*&493PVwEo^1pE80ZNM5+ke7TbRt+hRvHdjDA;S
zO{#c(vIts{^$eum<*pi0hu%?9<F0kp7zcgGQtx#GzURqEeH-ux(KDJLtc%x6uthce
z7PNmI2fZ~85mj9^B--Z9o4zBLz41&|6|{BWCX=es)Vx9T$xgElby+$NT^sPHo|~7@
zvk}Tp8V!1yP6<Y9qcn`(84EelQiGmD8}O&^$h9${CsGw&IU6TaY1^`N2GI){*81?B
zgErUC0KFI@t(xbdW#mr1@G9h%iI}$O&<B5<Ohc8X<_(~axn&|CMx9+BsikGk06m)k
zN^hQrvdbAm>!mXadY92hm0j<Vi`LvhGwwz~AB{(49=@ltB^c0xwqu}2$_-VNQqN^e
zFrWo(M?o*Sp$h4w-b&HP9W*s>6!gKcV9-guxT2{VP0bqxJtgCgl&uxKbIp-Eb(eom
z%q<f^yQZ!8ell}CwV>@F`dnnqpbFpk;G11P1N6~IvhH+JAN^$PdTMIk5PBvf*6PAH
z0YbCuXMkRE;T7o6V^Q>Zsh~FCkAOZHVAS=)4>@^%BIv1k^M=spm8=1S4n4=opl|Ef
z52I%jgmjU+^rDp1@N*}MhJFO}&Uk;eQRUJJ4okXpW`$lVLLc**wxE-Gk3|Q3-%2pX
zKrdvR#h}X43884E?Ss&_n!922lB;imDkC?<qLsD}8t_{$oniDLuW74`+(jt5*x9T!
z{93Yf1a>P~1A3JsH=<0&r30;R-XQwjO4hn{?jn^9+CrP_M?vp%^^Foe`^kT-ZrSGg
zLG*bgYerIe^P=@dt8b|F%^L+h?-nsv8Fx`)TitTv-n@+7<-#i^dKZgIeQT6D4*Hl2
zuXJ{O4A7dVP28K8(C=2VRz~h(#MT*m2I#FxGrej76ar{Ha_@`YkA?o-);~5JFP?VS
zHy?hu>X+M-&o25ezq?poF4uqk!RNb8pC5Gn0lf7uURaP9r)SF-@Y#9)<?&*Fz@z<{
z{`igg^UbH#{`=#{o5b;I^KRcCy7cyJuHJ0B^;vguv3pwS<7$0=*7aRFZ?}KD5|#JI
z9e?wxe|fdJ^_qNheSNVy?fUf7Z(fX>bANKO^xe0wRu^YCKlhjPKBj*|Z#|rMTb_0Q
zNtevdg4>@z-EFl_o4vf>Z@1y*(SCIP<LBw7PkvfnEW7P8amb(EJ-S=>b^NqjpMU@4
z<hvg-o86!L?LVfs_=i`Qzb(_90ek_l|NnS#>x=!9-#w<2M?2K3b$Sku6Uyc0NqQQ;
RrAwbY`XAY=VbH{Y0RXz*z1{!-