From 3e393e0f0d139ccc6b5b3f890e38ca666985599f Mon Sep 17 00:00:00 2001 From: mhidas Date: Tue, 31 May 2022 16:18:12 +1000 Subject: [PATCH] fix ValueError in velocity_hourly_timeseries and update unittest to include this case fixes #124 --- .../velocity_hourly_timeseries.py | 2 +- .../IMOS_ANMN-NRS_ADCP_LAT_LON_DIMS.nc | Bin 0 -> 86972 bytes .../IMOS_ANMN-NRS_ADCP_SINGLE_TIMESTAMP.nc | Bin 0 -> 106434 bytes .../test_velocity_hourly_timeseries.py | 10 ++++++++++ 4 files changed, 11 insertions(+), 1 deletion(-) create mode 100644 test_aodntools/timeseries_products/IMOS_ANMN-NRS_ADCP_LAT_LON_DIMS.nc create mode 100644 test_aodntools/timeseries_products/IMOS_ANMN-NRS_ADCP_SINGLE_TIMESTAMP.nc diff --git a/aodntools/timeseries_products/velocity_hourly_timeseries.py b/aodntools/timeseries_products/velocity_hourly_timeseries.py index 35cc60f..c13272e 100644 --- a/aodntools/timeseries_products/velocity_hourly_timeseries.py +++ b/aodntools/timeseries_products/velocity_hourly_timeseries.py @@ -54,7 +54,7 @@ def append_resampled_values(nc_cell, ds, slice_start, binning_functions): :param binning_functions: list of numpy function names for binning :return: end index of the slice """ - df_cell = nc_cell.squeeze().to_dataframe() + df_cell = nc_cell.to_dataframe().reset_index().set_index('TIME') # shift the index forward 30min to centre the bins on the hour df_cell.index = df_cell.index + pd.Timedelta(minutes=30) diff --git a/test_aodntools/timeseries_products/IMOS_ANMN-NRS_ADCP_LAT_LON_DIMS.nc b/test_aodntools/timeseries_products/IMOS_ANMN-NRS_ADCP_LAT_LON_DIMS.nc new file mode 100644 index 0000000000000000000000000000000000000000..057873738eb2e1526956d169a1a81440677d150f GIT binary patch literal 86972 zcmeG_31E{&))OemsUQj{C?DccBuSHAKrx>*J!qvDZ2<)|q{-JdnkFSl3#cf97bxqE z_pvI<@2-sg($dW)#*eeu6HGcWnRuT3vXSKSOW^L_Kpym|9x=FOWo zZ)Vn%6c!I3vi}f#zyc*`q?Sny)_mdMDO%{F+C*lD@zKi`whP(t=~aBLjS#F zyxMy4@UHZSHHh~i{TR<^-|>J7WBv2k=mFx(M(IB#eO8STQZ7Y-6a`Wg=sgPTxZ<)A z#nmNMY!JZghwM;rKE{go7P&Mm7XBB+iB^htmM;A&;ZK*Zs}B$l1V;wPJ>X+;Rwm-4 zL!RsM#d8|_clOSw-L)l!fKfljM#49Lda>Zx@SG>%+;-aMD|oEI7;CtQwI}6Fu@23< zX;X0xz+#qt#kavSJr)-oVVTCN*fL%u5%Tfh#P>$wG4emzG-Io~wyYSiD;)~oj%I8h z@b>tL4PaH6|6;}t1TS*lr~#~GG(17()4*H*)Tq7LXdG@9TbKggy7CR9nYBMW&0vt2 zv8Ji7OkkLTdGLp^VAdNpmWiLr7()u2wfPJ?8$Ar3Y8j82{+nUB?37Uu0J&so3tjAZOmYNcW59do?y5reBxdViY%F-O(fb!OhLwF3V~2tFUfQNv z>@Xz5s}_b48*1OIW8WbjP>|u^X=6M-mL>}~6$0M)_!=MEe-C&9dW-~b$r`_({T2xX zm?G&9`tv}IU51Gv#(3~==6tfo3yyFnFUz z)||v1!!b3Bu|vQM-TBf=HXn-*c-z2x?ZZ1xV^3gyfN2hR3x9lKH9G}SpUD`av*eZc zSF`&u2>=B!Y|g5OPG_GYVv8BWrhC)2Oe;duQIHdGaB z!A}5^7vTT!@RX7mSj821}jnJ&{oFRToVR0r z+c1=p-G@bn8pnz^R2|F246bF2_mzL(ODSVqP8p3EyvBI|GgKyv$eP7Qa{G(NM!@@D zcJA8J3aE|C_W>f5@P5SBfAYt(ih0#RBP!-phlcakJ&zRqFB9rlT6W3*L&#^2JlB zb{At4za7cMMeyGOakxvquzd>{R3eDxnR^_FTc0xLjyzmsOr^EnVzY=wd1^s5;ckA{eQ55%*C8J1IdeVga>!% z`(yU|NIa%~e`OSqRy8_XdJWIHY((%az*F@**YbZE6l`_8}a4sInzC zA0S2mG5~dEWjHf4oYNG?v^<9c{#%uBOU#v-mEp+Aa86eoS$SC)JXxs^dK(%vzv6T_ zGnIJ-bE_rXsIoK=Yvi95x7Y6h)V#f#*W5nLYknW*RZvq^*jF?xC@k%>`YP%}8s@ue zN(-EQMMa)lW}oJk)u*{-_f>8M<@vMvtlBWMKC3m%tj{V9GwZWDt17DQyE?dLY{B40^S zeqm|FtiCDMyvnlLzR9VusH(Ov>WT}%C<^r10EU_M88gJp`ivQ3W}P#`S;26fH(2cr zHQ;=x38pPJCrnhDU;^aS%}jEpJGwS8=?1(mO+MOzmm|>2Nz9uNyDm*H+CocG)-FvT zl0wVMHIc4MbA-0wY?*dh*lYmU?j$rZ42g+hJ3yQgh7;zuNx+!Kn53XOhot8e0!+2m z5D15iZEQM~D{ER8%cYyTYt6AnXsahygrT{5gXP{x1J1SuZuy~xpyrzDP+Y2ltvPeL zJ=2jR@y^Y(XM;;L9A`I1*M#n%n+S;JnwzQ5R!^qnv;1mMCRElJG|Ju?Ym8{o#&FQ@ z$~in5^A~EZ36R;Hbx%*M(Hrx4Bbui<9MeLvfHxRKtw|)J`C?vBiOAaOdX|~lN+(^o z_w&q|n$uZN-wl~6ov?%6%d$&-JZZDYjW>(-;I!Mxdk9P4enrEn3HY_LaGf?!3x<85 zkgL&R&`Cq?K&U2gqULhgb^JQZkvZMjS%d8<*iPDHPr-I+eCa+BI%!!w#kNZWPHi9< zGgeo9AQFvXaBZB`AYYAYpk4JF!!v=Owb#mY>PvzAZplSg$-6NZU7hd7Ty#xfcjl5% zQ;J+l!jZsFf; zM(@ndMj%1oqDUkhF_Kn`L_7kVQ^yT6ld{692ngqOFibau!+Pz1mqBMcbGyc7y9^Q5 zEp+v#7hQS<*(9)yiyz%VR@a4l8QD$+k#05vpscP@_cF4%$zb=8ve0Ie3^BdZ43fY$ zB3l5mYq__<*u{KJ-Ue5$PQ=gd;>5Mf(7Mz_y9O<5de>M<*Pua3lQf_n(yg;v;0yJZ zAX`2t(*?5N~H*P+l;?l;eOVDN;+!E>+(!5R>L&FDYcZYU_U2)Je zv2$`K`Xp{XI8$!-vaqzY#n{zMWbzS)ThZ8p@@=k^TEHjh%aC}$sPGy4taH_G(H`MFtSBy$~Gz^*mevk6NcgY z&betc+k)kT0L=S~wG$2&_Uh6e`G5}&u(BW~0ULOQ?0BBB=TQQquTPYi?Us4Q{C#mi#bGrgH+X}2o-w6Bdmgp1wy{8eS~31(Pl{?Q!k@24%eLxjnB##4x~ z+@Q?o8Jn;dOL3mhGqNxMHe}%*Q9g_K13ciA0LT8b{_378CLRt9W$eNC)*r#VG6_!J z_&lQ+KY`4TnR(7RZ7t)~009`0PduK+7(2y$hfh4t!%vXQVElj+558x1kcDqQ9{t7k zd)N;c12DxYB%gQ;#sJuRhFlQm8I#^0_c}bLTzeGI&ogA*F?}@|s)*_9oayU|=^GH! z*8tNu97B^}hPu@^I76KpYW(AC`HsN$N|Fle2?c~8QQ6NZNMHv=M?MJ>M@mJ^OqM;A zW$%QpKVSq#!uf<$(LkzeW~K3YhdzV^iW{XsLElKA*oU17yH3%Kl{{LNKKmnKA-<=6 ze*qL&(L8-}$3yW7SqeW@Gb!xIm*1sku+uWPukxEeoEp!;W2&z=51)6+yod6^1GD%@ zTVHJ+KJSFV1zjT6fTp`<3Y8y~f^s2Wq-d-1tKA~L;7(|CG!&}{1AyX~m>7yD>xSYm!q5-JAxt8*gnZJ$Rlfm`3`LJ1DILm;$J}?i)1=ZqPVwfG* z^xYuyU|djZjx}k<=6+%@rc%Unk-3p)UVjU35npH!PWT5blY{X?7{jR>jIT!E`oZ`s zjG-TlZ^VF!9Sz1a5e5J?4#tP;2IFRgp&yKwU=mn6<0nPoG}vVxjJHd~XfVEXf8UV^ zCFOcVfuw`6tR`I9j8_rU*E!R-8K!R@OkeTLP^J22#q?EasOXp1FId~7(xeh~mjYri zE|=3<>|5w)9^9FYP;!O(pX8uR?2+cbgBaODnEhcBF=q=%`*_TVcUd5c?|7^TVnvuO z7zQRrxL`koM+=5S4W5vAXG!tn9-?8uSn!0}!8|D(O%jL&S}7$@ zir7+8e|TwVSQhd(3!Ec~I7PI2Euw{d8XgSN(!^8Xw2e|kHLntFiPU?2+Ax?*79H2`mu$zd17T}hcS4$`mi$WBPits6qreQEk8wjsM;5?SGtQ8Jp z@q1wRAVfjeo>>=$)5-~fA)l?BQH&7xFhdN4Fd4XRLuA)&+5Fqve=i!=zW;S&wYy|P zl;df9WoIb+>q1e&%0T9>vm?hQiaAqL7M6AXYr`#+gRVoGLLcU0WEYfHxwswp7 zl8{iv`|S%y-zz260)#PJC#n930i1=wvMpMQeT``#02gc$mSQCs3a443rPyyI4CKLF zOEduEf*d*oJib2j3dT5Ch^|~IbPK#2?7&KqNkA^B;N{CPmtX)$ODLFfJ1|Zch=HJf zypw%;%GR(X6$GSQcPXHkRApu2ah1Fhki<-1MTSPe^lipa`OHwy`j*7>^~>~i%k+&U z^WmnOG4GVy?G%tJIR6Dpv0=RZh=ZT8m?X;%VG&vRX&XiH50^gMfD9$&QWQv0;Ac*O zG0zT|Em#Xok_=X6al=qq0-6-JT^vrn3_MY?4&n>MI!3dqCIR+s=?{_}0rnyllI`{*(ogwHKn@X<77+Dm@!%mZmQo{|sZOaXamet+@I z*p?m6G%j!sz<4O^`qNLBdG968d55q?tFIc)V-3Vum^Z`NYMyvSh^15t%7uK9 zg2Yk=E}EWm6g!h!n2W=;IDC~J=&*mg)Zg-xD(&X4>Cfpk{-P~%Sa zJ09WBr5jIy&nefN6yWub8gAi`X^fT41|C5T1V@a=o49BIaoIq$lK6~Hv>%+5a3@+o z9JFW~*fY5k6+MVvMQgyh2X~@*;dF*O(SqW0MbTb{(O1y`htsQQEPK+cXup@g_(6<9 zQofD?V(j>f#ZrQSmNe*T;Arfi+7qF~PdfijV@+oY*tVNAX@uF6J<`HaVaK!z|5WLd ziPEETZ2E?f`%J|9AN$Sv?PE8~B6qBk8FiG1W+&C?2$_sJq4eV&sS`p%zmIo)Q3beo zI8ly$7B**(g7s#V%u1Pe!d!Owx(nYot7KO4MH~I*y^RC!EfHhiwUYY~Pa>a?HLGO) zwncn7L`ES0)3CR+s4*AvsXCZ?}rrf(9=P?h=y z$MjVwjA~^C*Y3Ryms3)%w+UX)qz{57T44iG35sF$pX^*w|xBl z{{2B2@q9ZLJJO^Szd?#`uox|zN$iLV>~Xm5KeR7)eDsruZDEgI?iQ12d0`<=roWbX zC(8Dxm!3~IiydFfRp>;lxut)3Lx>#?+1seokxw|9#SZ&c@rA?<@?)jPJt@VGlaZ1~ z3$f$1@gs(cE!+iSxaV@nIVfMj3Jg+^%OSMm8^aSzLxu1LjxP=IZ8w>^?LQ?L7l6j) zkQ-2%fVCnN^Y~iP0(fBCbcYB$kB9-nd^UzgOjlrOYs88aKSbaqi5N*4Y!U)pDNZKf zi${5wWU!@BoQbF9bbm!E2eb?;r~Sg&93hc6iuLf^mAdRge(*!{-e4&EXLNxK{sX>cT;N@y0lYqp9nByN-c z*oc9rOR;+ej=OEd?n0f|eJ--q-VnP1Xk6^R5hdgf#qJk66}w-Sh{eV3gF6$uZ+ZOI zrvyb(zr8AU%l5+bt*ITcTebjZs9SyeV){B|`ub%07PINBbEc$faMvgx1iv}GE%J4F~xTFg1HR)SU7Vqp*Ia`pgDKABV6a{`J6tMhfY>i+p`}F_0 z%_=ubR%4!34q-2E7UVljvY2vWHDzWY>t{md6dJ!!3cPjsgcC&nnEf*td~ZD}+2DKc z-Gaas_Xs~GeRe@=@co50y%gHNUJCqu#CtWf3M$KZcb0k!(Y3|28%tiIiSpBsMr)1~ zP92-Vr_2)ZaiWcm!YMw^x6x5J#V0#C3cs^L#3OS@3MV-9PCkW`0VIVJd{WI*IK`Km z{Wkav6WPbZ2_LfUGQ%PLjyas{peda2F{_CPpOL0~DgNH3aKeY|SqMJC*~b*0!qLtN z9mOa3qeXB$oZ^$Ev^oA*Q+&c_KT|lt**{76I8*rBst6cw3MV)Rn8GRkfu?XO*922I z!BI@%1ZSctoZuW}3a9c;GKCZTgOh|GVhX2x4>g4ooWo4v6hGY*PIyi>g%ccFbEJAC zIHVMzaDrns!J+cnOyLCIo+O-%90)$;>qrvcX$q(KuolZ5!67`eOyPtNC<+bn2~LhF zoZ#e|!l}GdP2mK8T9WYTrf|wP&lFB@W|+b${t>2d!qa66Cpfa(#QQJ7p`%gEbem~{ zL*>mkg%f;tlJEjkIOPj^JMJi7f>UG)r}*4z5WZ1-Dz6+);^73RBng~SQ#ipn$`nrJ zooxyy_@KS!j_{%Ia#J|rFNdo(cvdEfUu6m>Jm;9gDc@>SIOPjQCEO971gF*%PUV%s z43ElBxfBIb6i87ZMS&CrQWW^#q(H~k4K({E`{lUV8cpvd(P(zv84V#>Jf+YypyAu* z4=J5>$g)q(#%hrp#MgvL+dqQe+-SCeh-@&&dZD~El$CTD85ZDjjlV%Q&s+L9Y*)IZJUVv=Sd zHgmIKe!dwpDx~-4^-N>l^RT^UFoM z#xL^_E$C};e#uvr^!#!lVyfqte2bT!U-I1@dVYBivRECze16uIkBiczetUxgJ>ZwJ znJ|6TFn!gEoqMucnZ8*F12S2^V%1T$k2MwFo!1-XNx|qI1-iyBvu~A9cMHEPeOo}c zu*34szOiRK?&>5w?yZCG1KygD2|@3MY^tX+VJ0R9onf5EorVE=9`_S6EIp6QK&ID* z$F0=Kgm;*DT;1_6Erc{F*Gm-W0gsc-gz2k>>8qCMtCi`Sh2U{3rK()tK31Mug)@be z>thNCj&cX~RB$9MkuO5qN=^4C;V4(#HE6V+qs-kkj*^9FRp~g&YBNXKY33+=t8RxJ z#u58Z`3H zdFZ$t$>YAo82LILcO3$U)cX$WIwLX4_7@$MZEl5)L70X>f^lrZ(R!u0jd z^v#6n>zf&>R^ND-zDk+CI+?zGwEpxwh|?+8*A(cO$E}l^v@7nBcve#0Ewm3o;F ztU(ZRPiL79ObnJR@mVHUi|P604~U0;mT5z7sFw-(zP*(R_9(7Q(`nNF@9GUZAx+Bl z5(V`9GGD^%WaU9N6GPQ7ebq93wK9FP5c@J$%2}qqeO&t4{MUM^NGUMgpa4He>I*Kw z0Kr<&WpB;VeA9Uy#7D7+xQncP1EutFc=Qca(zE`953kU|9JLM9%+fR&*Wti zk?!ozj^YCECzGGv#3RvV_m;cBZe2uryB|A{L$Y6d(!Y2lIz=mS6qjXkKlT!jG-kpR zt9hjHtB@K%B@jZt6OhvTv6p$IhIg*_^GLMFW0}Fk(j2p#*q@!rakAWb%pE)uX*Mmu z8O6Ga<@o+=J&$z4r7zyeBhk4s7O-CYv&b$UsrA*B6+9A2`xX$wKSplM9+>qEMz+9@ zXaRSw{mZoCl0rB8o)`q?V}JD7V=Qneb_^=Mx~RCQx~QU{$WvBY14^jT00>7`xX;QP z)nf3}v4h&C%uL~N=sGGC^3joh`U&D)G*6JfBgXfZ zla248GmY<$uQk4Zw9fc0-cr9Q3jA6r!23VlEx`Mna!=;CpKJdo6$SY_OjEH3$ZlR% z*HuRw>uvjD<2&&W3MKwQ@5Ddoo%jd66aS!h;ve*$@>3M}c~gM*f7*!4Q7tNbI%0vV zGls$V?7EP1_$7M0riV>W~Wg36U_t8z3^+9igr%CgM0-=T|z(7fJJ^WF` zim+!+fu~uE)Q2NYnqQd}4*T)8Os`TE3D{ z4x&$z;uAX}ALQ-+X#)oCF$gYpPJZv`;W|1jX@QG5W6KrHQ2|;gEzL{u0#LXj)U(tZ z(;}W$Eg1Ge6boD#f){?u6j1C?ED{c)--@&dhXk%jY5XHfmzG!7D0HzX&nIeu`#{I2NTH&gxkUhOMWn?T}y43D# z@D{9#4WJ51&c?=_}(kmk|x{!d3!@R0p4cxsjfE{nyM;SKqK0eB+QtQcgSI#T))Bt@dtC20?Ew*O1)ZII0@fM7ZdM47Jik7q@J8B7C1km>C zUYs&g&P__x=IphLd)aj}O|W+9%HG~=y0c`Vsp*!O-f@$jzG)xD-}mJ0%TuAneE$0E z!a!5N58-WQju4h|$tl45Kb?Zb<^Rvm{sXLp^d4ivEvWk1!>y0b~Ei92dIhe43K~`#}D;0kM(sp=%@mhr*FqV`>WUOPd0uLtRpxc2tGNS-^i$ zlqoJLb;nP^M#Z#0$?8h&X)!(qBA z?VtvKHZlWh4MA#%kN?>JOaF<#;J+$-k`~C{Z}SJAaIO)#Sw!GcLv9%x(_N{78$19E z^7Xm{G+NQWqXJd`j;1Urm!d$50{;yP82dkBF!z9YZ~=BoYS`^kLpBgQsnn2IWNeOT z(P#@StomCtPb{2ReoU=8{%ltrGjvtkg_`^?N)cRYh>!o+|4)4vf5CrS_#`cmKfd*| z&p6j`UoRr?-_Msn3|7f0uCd#>28mNh^)_=2STc__ruf3omM@HiN+Xuf>&n>kySKnc z!Z!JOp&G4T^zTDW@5A_ek9C~9|NBn^An4Ep;)io5)?2_p0<8<;)+A^()Y8OFP6#24 zAB(hT3fIJTv|c|RViK)7QrntV(7nqr3vk^vKwwPCH^JJ4)p%uVcjr(Y8Ioi{dGysu zq62AKb4=fJ`FoG67gDEOiUP?gVC?@GZ*J{R&_cUsp){{u3nkT#Ik%d`6*c}Rcwg}h`DzUIaNx;+)*6)-?Z&ay6v$LkK zKTwYrYSFNQ){9CY23%>WNDMY;5Q>b7rxNCD%Mqei0o`!2$gzFWX(>^^b zE6cAnL^LgGS1Ky=^AC$Ev$QbaqO?x6J0Q<`)Cy|K6s2@lb+I!O(b0l{U8^S0<3^Ey zgvbFo70(N&vxeKQiLK)SQbXCBfTBV1buzcM1OBbKG9Aw88IG(BN4AoV0SZ|1#TqpQjC=fG*9qn2 z=6b$xs}}J#Xl7%*csQEs*)41wHy6WFd81LX5QS>=NP|Jq>P9L@LMuQwlVuNO**oPZ zM`m}_9?0MJ3A&&D_6xRf>-vDVFDD_jYcQ0I_XQ+a&q>z>+c(TW>Kgs0i5LF6`E&8I z`4i#OY=QiDChu`N=NgaDS(Z57eG7Aq-jy200aR)e9q4UM1O2a3Xux|jkQ5{&D5m#` zruT#RdrD{kGNjz@pn$ReV-82Z^A((_+dZjaZc1uMNeyCONlI#9AWn60!27QY4Pu2d zKK^6>|L8mX1^+AIle9qovER(z&bdZ>HzdZ}=NOL{575SWLP1%64Q`_t1%J3bgU0^E z_Kli?{IVhlAHd^R`(GXK_pMU9}0yuLJ*~V3!BJgPnj0i zBg8bo6RyWSWt1Bv+=BeF=|zJ$*Ql)*5%Bll?q!ErNXNzo6c$z0mJq%cSnw2X zUmkAWHbq;QEouE^yix7!grfVLENP^X9YkBG=t?I`8euRvzNw0?b+V)p8HZO?mX}uK zmw7lc1@lCnmb8H~^FdYBMKzx4qB7b_ZAsfhei}=k%tfGG)|NCes|KwtvOKz*%94hx z2wfV=SF{J)l7@O)bo;a?cG1>sOB(8z(WRkIhjwUN(ooHeZlCs4J+yJ#k~Un#VrjT% zjrMO_()N_W`?R-}t&?rty#D{XP(ncOpHpU8x*}_$Jqk64WYokCiJRV}$^w{QJwJM* zyE^B)4m=3-#L4&M9iqb}*`7Q@W5Ff3zSxa7NrzY$7C1X%=iLU?mfPz`1b>|uwn`iL zZ_j7Q}PjOuXBq#v!6J(eY?(dnn|hK9@E%V{y!)V_e_c zkuFUy63w4gGc8+b(83K7Z!^e+AV|`NU?sXd8q=CO>alvMRhMLr4bCtJfFeabm7m`XyPN8>_59pnQ0D7lkp5zr8&*_~8 za(br`o8D=trgs{X>0Q+Pf903=|1(QL}Tzg+ZBHahpDx6Yf{yS;zjmS0-$Cn8htr+4c8^iI8> z-l_K!H=y26@6`M0oq9jLQ}3sD>izVt>wA6w$oqf%3_Z+$+dS)2z0v>mXKb`H3xqb` zJtY!WM`x^HL3=u(H+e!V8e{YDh@7ZSAHvE z^GAwrqwWX#gClJb`J)fNWB}(H7ojy3xnY6WfUal__JstLce^^I_ImI)xj{W8PJX4@GI!qiBrJ95(N&R;;h8i0$OJ<81crM z!qMhNEdmpWIA}bi$?F6280y3D#TUSHMupi5loL*M0#3mTj0PdUK=p@yV&nD2fL5Y7 z9Wfh&k|Mk`)DgQOb4_f8pk7fxtzPf-f%ZK)&q3uOH9$_Y<>lGcUEhr!|GJ<%r{SdMR|p-q8c8GYUeHRfQ8jGXQnf&)&c7> z4#xrye5$UjwdEHURC$W$Ih^_c8^#4T0O!huw4f~LUPCV}VFIuSMiFMyy-ha)@7QOiRd3+y4E886HkW;Y(taCl%Xtp=hDscmdg@|z<{ zmP2u7<~gS3Wo2VI9av5{ybb!a&1}N_yGnI1lvYl}%VQa93!ZvT6PqDQss8hG1YTSw%ZYv)*tfX^Jv>T{6YbL(7 zBW_x)Ip(VClO{FYk&~HS>u_d4jfy7aobJfYgeH}pndO+;UX$t*ys#g2($t5bL8Hx) z3*7*)&CATnbLK)T6de42a0G6L9@CK?WM7yvwIv11FU!=nZM|kXv0lA((E!wGpiT?A zz0rWLs>ojkZPW2pB!2Y2DZID6ojBg@=Ge5&hAslt@ z$wzg&$6b_PUenn$6q2;VER)PmVPbIlNka8&sq$rT0K#Bq>3Zy9T3!{K>{*Rje2e-c6U}5{URrPvW!#;0a3m(Kx zg+!(9Z1^L+033r{P*|)K1OuQ2*SIRKI1Jh_SOHhk3yLQzhF}QVik3lsc!oc!?=5U}yKuPGD8$F~ z`UQ^ogdKGga>8z52PB7~!v|WM zwNh0uzqBGZ+rsuv6i>f-4a)#7DT#DC4bm95Q(fJFk;d_;D1gi4Wxk}o!uM?ZRC^{q zq+G92z}Ww>Bdx#Xx4$Q??g`rRrC=r#rR=aJOoPTr?K7c&i)-g{L3Javf9=gjtZgP9 zB63Ypb!kzJHGTp=Jkb&SR@h9cGY+*X+_JJff0mVNlRW&1Ym=0MnyNypQczf0EkcV> z$HeD(1y1Q_NhqP{QWEAzH z%%nRsnqdd9FVGAy18@&VG!|&^f|^wHV7(Prp*K>u+^UqK>AhfW?Wgo;P2NBdbHSPy z26iG(u-do7B4FL}nRG)2>5Q8r;p5>36rN%-M5eY<8j8X0U>FwspgIOTHKno+ZqjJQ z)zcbofu*FQh7&Hg7|Pi|d_F9jAEKZ?3q*HS2Qo=}IMM(Yb&v^3G!WBZudQElR*sH- zl~Qn03xs1u@Yr%YFkZtEYZS%ok^TD0Sx_N+@G?ZhEfF7}fViW@G7+O3tk0m}2Cf8I#@-+0ZtSSe~tYgS9$Il(f=F^kOavB z_Hxl?SSEu8MT_6uxWBtrzb`NRy-F>uJ(&8XYQ&!wF%UyPUpW$Y#NaYMZ1cuF@E(2h zz%&Z{5qaBh8W~zsAiytcp(d6v3xLGn5YGPaM{ueOOO^I{93DA?6iPtwz&>`q8$19F zJpvmc5Ma-Gv~tx<6kAO2QFr8lD};D0_o!jAY=FC;8E!kxzK3|TVt_ApyX)u@g?&3g z4LFjHEW1D#A(`sog9;BCXs897Hzzv$U`^^}2WNJc6Rt548!u^CGfT9Xs1E+9`uS}| z&E5#erMUha0w0V*YUDLRQMe}v4jpK0sk6i0{tR!}AId;RWRvn&MjhN>l+gsTs20gU zQ2c8Zfs5K>%fK!aWn8d#4S)ekRO4t*(P1DPvM-`RFKmVFOHEB+K!U6;Lneu4)Gg1* zwmULtsk9k%w`hh!%nu8tUtQISknqLW|IyOT^TO#jQC~&Ll*CB${tAf7(Y*m|N{C={ z1w?q{! z5`pzXyHX401spXzkpdL6&5g@JF%0Rsy^2(ulBJ6RiehL2xau!HFNJD1LAV^I$({@9h zJu!XLVDBG5fp*4Ib4tQ9y*Gko6VM{4!F^$z%&=DkBEl&im=})?x_GehMnteOK>G~; z@go(Buknef0t}KE;w(id(Z+B~5cJfr-_PrZd1fm}DSX}tTJZeJ!U|4Xct(J{6Jy80 zs#^}$7h=TM4MTv4r-@9myUSSL+&1*0|9f>2FuwNTl4*9&;dU?-1Jd0ciYu!#-Pm=d zGxmQRcMn|bgDoI(2WJ4T34^Xgv5}Xogg@%)}ru&5B!VQ=Yta zn-rvjJcaoVV(JlQgOH=(g$L{M&ToNpX28zK8&WELaQC1vCIk%$%&1OK(y^I>s4HgP zI#>yjTf)>UHxTEUIMBp}-Qk=9NEr_wBs4$aqL%^NtpCD#HE9cUx)m-iV;$@B)QuP1 z5QAhVu(R|9y`XK8V=ks!3V+l#3V}XWjn48W*W~w*h{HW1HY$z+SPzLRSIN zi9?X<*?jox>a-gOAwWZ41yL8-K4k3!DGrQ_(N+@7DM2^`#6dU%Mic?THUeE6BsMEZ zY|YJZ))`wPL-&11NIE(1epCIZciy({ExS zV?6lCt%l?Fe;9k6OF&;YUHYJk=ixAY;=WD#FnXaQYvQi@>-Ot7-8r?Ex+QN|-5-Bv zb#JxV+(m0`?j5(=+}FQha}WR8=AQDC&3(p;Hutm#ZSKQX+uRp=Z0_Hu+uZYBvbq;s zX?6d7nbp1II;(r&Sym8qZ0;LAu)5#<*6N0Bs%@QX$gG(Y&bR^o{x$Qw`$y)VcfS+< z=HwGoCFzbahx?+jc6V{E-Mu_ucW->y?!IP+-F@Bu8Sd%X8E*fK40p?@4EI1whI_@| z?CzVcw7VDh?e6g2c6S;en{|iHJ?S-@d(BoG#I?KMIN0t!GRNLNQjK%x&VAw1muDh% zUk3lOTW-1KRq!$V4H>>ohQBSt-;v?(%kU3n_;wlosSN*IhJPu;zmehJ$?zX#IJ->b z+YfxqcL4Z^&mb8-1pF_Wx4UXKY;_H?UE>N5-Q+ql_(#`(&{o%oyFPF|=04B$g$D{yF!;->w2bnr7QQb8Lq>>SmgRgtyMkt+kMpUZ+*jc z;__!)Gtb)O+H~LFU9aBlRyX9$Q2W0;TD@@ae_ZRef4QdGesmpk#6q=DX;S}uX_>m? zn5k;cq!H>9$Bs}h7`jS*CGQM%^YJIEulvi?4PzYYwa*`*mQ>ZM&QCpR>%L3Xi`UFk zxAo6amyVdE&V0e89#~kb`W|dhUw+7^UUgrAnsMoY>bZYdqBdS|nmXVQYt%OmT%)F4 z(4^jKFH@IIYE(DAuw4D?>1U~5N6uDXJ9@Qx+2`}r_YMfFZ+?54Y9D!;`pu5BRNLCK z)kT*rQ~f8cREIyeMqM&?o%-RF3)Q0Hi`A{~ovGTNuTuB?)}!7w{1mleTdVqjJ)*Wg zvrv8Y=sD`TKg?H0P4TJahkMkcE%oYO=9Z}MtMk>p#sS`!`qd{6X;Sl}9(8b9xq9)Q zHEQswW7KmGh^o`}kE(?;W9n@mMpVb`hpO8S$y7%zEmo)9SD|iLG+RA#&=G3opd9s+ zVa4iWcbBR|Cse4Lj;vK#V}&~Vw*KnP{YI&L67^c$<*8$3Yuo-j)NxbA7!?UUYg zHKo7i+WO_^u6s^;#?=`7ldEkv>dTOM8^)=h-sDsN^XA3sJrgchH{X4=y5fuL)NhvT za7}%Dnz|v^uim!fboH7)u2qYkI!oQ=`por^GeezsPOG{xx=wxN#4FSp`(CbE-+9-y z{)frx=A)Zb<Rc%?M9%?^Bz2VWT)$4sXt7Es`s@_}ih0FRz zx_aTr+HZ#Af`wdbpgzSyWLvCZnBQ{Qo|xpJU- zzUEZt{%(~z;l8WXcdoxteQWz$uB=sitM9*>p+5WJYPBKvCe@yKi~9OU-@3jzVt@6J zYhCKi<5sAfzTT*=_-TupG5Z}rYW;*;)G6vP_2Ga5Fs7<& zo+wqdOIE619(T2R=vGY87?A*Jr3V4fCrH|Lr_=@V~daUi$M`_1I@K)z>a6 zR@Y2uR!=|hEOlp1uKMCDZuQ-FXRDXraJ+iK##QRXBhON2U)N79zha#F(u5<`StA#! zo2`Cz@gG;JE6zJoeV}2sy7zPQ)yx%5>e8uasukPLSI@t64|Pk)Q1#r8$E*7vFhw2l z%uIFRkfYVT7rpB`b?&#WEv5UZ&wY8YYW>Nj?mxIj?aMaSHm{qmeDI!oW25fHzoVyT z-#g@`N&GLo_v*XDd-uNw{UI;i0Dk{_gufb}c^r&?)0KtyLLB^q`23HjfpDJ2z4#p- zAJb4i_zZFH!Ecy9z9XzJeyl}!JRRmmIGuNY3&zibZ!^H3fxdjsgR~Haf2+BVyd2YF zIS?-9fp8!`hmZNn@<7;w90ySjfd|H?KTJn3Fb@6^EX)UB-HpE9>F_~%NZdGnL?508 z;@yMk5g$x1pE<1o*Z4R5j`)y|x@a7=vGb{E&URyD5_jp7K)^NKxRIL;*g9 zi}!!TeAHtC|Ek&TPqkid`Q1X!HR2mTk!$>+@Yx+Ht|2Iw`b|;b*G&P=HRAmr%eoE@ RCPsH=W&2#?%$jc}|33qq3o4ng(Mq$z&n(xDn?!t|{dm$I;)jG| zCnqED!6Yjbd2&(b<8#FC9Cq~Q_a0r*y0{808pqglybI)>DmXrVm)dZCQ@`)^JlD3A zYdeu^nygvodhgnK7uK`@7AxFIyxZ0;Pu-;;TbaX}*c#p?ky5r^ybp{P;`Oodes08!TAI3sOFXpk230ZZFQNzD~&nfxr^lk8~l`-m2?%7jI z*vuXA3$^4UI`)dzQubSF$9aqyhu7zcjs1caG>kEow08FS*}tg)sI~;ri*H}yXCKfY z0Uy-IRi_66Y@XdxYTxYj6D^k5PUM({=(>tyyV&mu`~g@v5&gw@A-2674x+8Uo^e`& z^-RaFV#ekm+F^RjaqJdKC}M0MM3H-+Igve0h4*93jp&6p?^?$mriP)VD8-KaV#7(S zoSMFnF>2qc=Qf_i3MdC&7a}_Bd%rsQ{uj3Y3_k}O`1o(@7%J9Hh1c5jq$1DJN5^=x zB5HX71LYvhq;bB9F+R)M=uHA-MQ?sdZ$LtvWL0u-bd0SdFc{=1R7M>w_A$M=l4c3Q z1yn&TV|=d6B$zLZV_eQr$B)$#7z`$1n_X=`5RpXZgs-Nw|Kd>5THAop*g*JP$=L3Q zZg`~nlvDZ8F*fT|nj;8%-1}_Pr{bsn`!&LVFN7~2I=^@=a=W+rh1)B@Jae)Uel z{SFRRaK!nKv6&0;>VQdX9-?DUnDcA4j#f7;F<{Q@jgA$+W492DLl~Qj=-pLyPqH2D zIkg)S&Yu3&pB#Ls==SoL*bT(1(2m)Nu37T8jcg;qz}mMzqGw;f{Y?kEd#maB_Z=+n zmG`dtz`-hi{Qiv3*pUSEP{#H~^u@U^f97D#cfYsxtKQ|it-?DI`2642vHqXRnKRIR zDy@D9`)>N+I9la(`1)bMUjJ3`#vzqCgWowMSqC?DFg*uk)40naB{K*`8q_I1( z91vrG44~ek0&ih~_dvsQV425*e{Lh%omfyv! zg~qbVrOh_nxUo8z=;W`B3NzqC)65aAYv~x)wR{Zgs%&Yh8Y>)DR@IK#cvX)f4$CWA zYAd~Cg+*Rl;h5G|G^Vu`k5z4z_2rAkY}lx>F&j0iY|I9YDjTygYpQM@yD_86#%#=} zvN0Pos%$LAthm@{t*&nxgBgvQnwqOye9hH$V^dqp*sPq?uiDni#beW{=Ek})Lf_)* z@~YZ~MPt*gWsP;MV^dRAbyMqDj1?DvNfa2f1&k^iGiFGYjTtkf$_8eJi$c*hGt_KG zI!GT1Lv6|PLQx4r0rD!E$-Gj};1rV~z#EkEkpbR3(Y$#ny@^^kDD|Q*v@Crclmd|z zT2VIS7F^aYn!m~VHs7yvUWp{Zp^D25q;cx@P7=(kzGq^&VoK@AK^@SWvU*Y=iwByXzC@>)@R>2o*A-1zkwnl8g~)0Wi&_4J36qFDT7$2$u-F)23Xgo9 z^XAVRXr_-srN#j4VC1SC)Ev)B79~lt$OI?jBp)U$m3{@`v;+fIU9`4 z1SV;uA{c229&0V|FtuI4I!T^*t+|~|i!D$VK63*dj z560pNN^VWE8se*Q3)a~4=6V@tBZ9wuf%vsspy5g~&tu$_&qgzQ>x zb~vY)KWug^C{aZG;z25|gN8OJ5*-}0qSC>!lEFbkPqPrvF!46fDDaJX*RA9rf{^Z* zx_0d==ZZJ_@7EIp;yB0p@h;k@XD*l6YL^@CG9JGp?ZHns_Oktw4Yu~3!lSxGNvmz^H4u#KLus0a^} z$OC|lA2%+rk4U1w1X`XZ*^Yfku2TpnlanpNEq2Uaw(Ath`me^o*AN50`D@`B&a(-! zlk!~GJ>-~-@CP|BZ0|eFI1%iadp_Gq{Y994Vs0aon3Hri;VkLKT|Bk_anBK)AZ8|elr zvbXKK9=zfX7N=_KMPGmCJK?vilp26`kaPYa$KUpQb_ylH!IOs8fBv~|u!RJ}E!^b4 zxA*i<^2YDhiWD6(VZioHBu86j52j;O90TAW0;bZ81L+P*a8LnrjIT;#j&XI2se@|v z-R)r8Yd%wFblQ+HAVi+RoB>20F4$}j?kKq2BsO7N`bnY=rL4=DK#%-q|3VIBp)jD3 zbAelArN8uF3OSS6oPunpI+I&&G)_*kyGzLZw1GilC-HC#GY?r1mNe(UJ z0Trfk*%YQpTooVhdaPJ51Fn~$OOhGPV7(Qhhf^v=?6qUoO;Ugk8UtgFfggUhY<2%? z$h0M*QLE$`p;1voF8M^%;$lrF{Aqf+8gjEeQBn%Vu9w6s)76lhZJ98YEAZ>{*Qdj3 zyH3vL`|1ZjT$&E6?S}MR;m*rW6KS>@@;FU%9B-b0h^ryD(=4{7eXfQSQ36hH zkxs#tjGy2KX$-XXfBzR1^MzhwpRjoB^@|Q>rd@<|JFbS*P!XshRGIJ0b!V_+?F7om z^^zwj;VYpIalPamdWBlHr5B`^@Z-0WT>QN2%zwV~bM`sqKubvt;d;rolz>BK#05zW znX_@$U+`0hOc+qr5PR&HI!ulsV(OSPb*z|5KujG0rV@^WB$#7tRf2Pjsbj?7zLKB7 z%+wQ|HXIBHL85M)iyL5QVp5QT!?$ev%tBT?g%y8=sXk$Qn!V#7IPd~I_Z5qBxSFD- zNCDz6wg4gDNdaOTCQU=>@l?!qVi6Ni_atVn3C8J}Ot92C*}`v61+|j>VlvMbibgto z5i@L2UIfZcRm$wJ=LQqq0c$E20J;ijdtNGLwOcVO;|oIp`b(oP>2wPF=rBSY zW2ZEW_{Vn4jv|3r2aSPI!$9c57ccK$^M2i~`;U?+C7+eC{odXu9I&M*mMkeI1Z@Cb z@}4aCDEl1Of2A%bqpiI;E#mDIarF`&+Wx$1NxHT?XK8w_FLQkNhz&s53&{;+ENe~Q z09+ChZ|EK^3NpT9*@ic4&AXCd94dOlHSaGdfpa)aZrT9k)&wMwqMh0RT#P3Kz5#f) z4TD%TSG4~@GaiaI?bUI;`Z>z6kI;#^8hs0%9ou`gR0Oq<)XdfB3n>BeA;}qhdk^>H z@Eio{i+kDoUmWtJK}k9c3j?a=Z4V}0L6k%>Q-{bw1WctFrVc)HjI&BfOdVgQjxAG( zWw)L8xPDl@(S?SJ0U>Whw!96%e^${97ZM!(k`l&b)^G)X!CtoGiMXqdhwF(q>7X&7 zF|dVWpmExHg@W@y#Tt~iWnaruc5C1JS*$WM+4X*bnE3$_A7#gx?TKvVSfFwIV^7`0 zrB>S0;8*I(TBS{LXFK8hN7+B5CM9D#nQo;UOn}ImEGd|*K~n2Gvg`*Xr6+KCbOPlg zSvZjUoL%=f#O0qz+o>@BU7l+q<)XTg%O_x@a{cj@gMKTdR@n1Ng zBF7Eu#T&ZYkAjRZ`N{sJTxyvxH}BW8*hP4?JG=hyS4-G&_&JN+xq5sVdjv4{WA}ER zQpre={wBM9eM=21K%?k;SnqxX7|oYer~q};!K!ovJ#boAfS;Nwp_Ya;A~XcwDo@DoNPT+N&$aX|n{*+5VU=70!-en@BGK~O*vw4e=| znLLPrCgp?&K@EIUke&s3k+#EwprE8GDd=UYd@BfWJNYPxWtw~x^!vl7-$=+59V-|R zYsa}RS1u*G#$l#W%>UZzWrpWbwvBNHH-6msed<%qAwnWjlhW%n))5k|?6H@kf7{R5 zO?vDGQDSTOJL81$J9Aa44{E)#h!`#O1A`<`kSFjBdS zOHVvWtbJG7+(+^x8H*Bf?IU6qZF`vn(~!E?Uc@aC4NA^I)wwqOMBmbCOZbt zWaOdRxsm4sO6U?!9Nc+gZ+cac?mV&B#n)BtJaL!ZbmEm{^~rs|>{NoHLk0|}78{TK zzJ_9Zm>eVG7<2XzWwR-QspG^{YGUdbGnFKmV<=UEW9kqJUkQ7FTW7pXTgp0&6b8iF z7aTvdwQu~Gxp2#+ZZcbYt?2XeN)ex6$CKnrxnenS6Tw|3+A%wen1K!&16v9PDog(t z%XmvQleG z%2h9Kb4R*tvMMUxrbMvd=d<<`vPqLY-{iEDu_)13HYvDIypb}AjJuCJ(zIoh;|V#3 z3)$p_-M612PWNE;Pe3{sR?eg@)G5M>JPbtXsqLjIJR$JyCGK)rsJMPErn~@}6jn|m z@DN2v7UZJH5%|GgR?e6GnHq)$bCr$!I4`B%VndG7`T4SE?1o8Egv}wKE5)HTz5!K6 zMUX(z#wV}*Y?|lu0DkHK45+e+J$6hT1g0X-f&Ao&L&Cp1MrkBNlxY%r zPwv#``i57I6);k{zIoechl|y*!CnsKswiVoqVMY1dY^a$v6TlImoB}24_n?of!3D8 z)8*|dMczJ(nD0o)+W?xBw@)U~xLY$ZJS=0N=+H$ZJPp^V{BFWV(OSObv&6$F|m7(Ytn8ck$&$j-KyJ&M{wxcuD*n#f>CP4nmx{M3 z*Sh;&ZoMwuxw+v->A4QMdf_VJY$1CgZ5uL{wWeQ{b6iZknPq2?@mYslivwSEHBg>? z1mh6lGR%FMx&B-(3UJru<@EcP!a;8Ghi_J}ov6qX;SjmY#%&H|yOZstRQNLUb5{K7 z3b;~`4THNbZ=@n#;lXjLb;nI?S2Dwlr`sACyY8U&C)iSgfgN;G=ec9@ z9{u5Ir5pYaDLRZ022^JYd$^f8)Q(|xkO@aIL4Q$^syxF+hc^fq`?dU z1O4;Mn{9KVTi!;KZsxPiZkxw1mt97aoLGzdqm6UY`Q=)>?i9bg|A|E_)A{Avt3{_% zxlULT+fndK?$m?V(&j5;QKGMNT5lHTw3GoL}<8AS%C{ zNG(e);=1e)hENr;axt1H<5#HkmMWXqY;*OdVFHl7-l4w8vMt zt=aUkKWnpB_08{-hA?;x431xJf1M3+RQRQFvE-NLQTulekH=k+g~#_uCYHJ_a*g3<#7e+JZ>GqSgFW_cct;T+aEdFhb%gb5C(?9<7_fv>d-KCXqh^! zOeG7Zj<18q9sAbjmyA%CG?*b^AYHb4{uWU5n3d4kLtt zVelQBOqe<}OdVRL4l7g1f~n){;5)ug;%ARgmo%6mV4#25%D(@Q-6uU)rR7b=c1@e7 zi*9;MHn`@cRcq6A(aAp%Z&RG`{_g6(rSr?lx2NZNu(9$+!7pVmh?P)~v3TmsFB4~o z9MWbdJO}|hJ@mzk*^??LqWr~3oez_64kPah-0hM1?*+A7% zVd}Uul}wmA-pnzyD)BINNSQjEOr?)r_rvo?s!tlu;4#oYzr4-1-EEXk%cG);-p%hw zCoXXMhChr?=a*%7iQc66F$d zlWl&f%9qUwzZ}M1hlgat`u^n0&xsqVdTn01e94a_u1v8T@{(| z(Tepv)87~F_5#l&Zy<01L!XLFi^sA1c&5kJtnu+oawo{O7h|uawdwEkuK18=lIG7$ znXbL@jzXSEUT@>VzTc~9kUo36C3->YvRF&S z3NF`WxxUKSFF!Bb9x^+8VatpJBOP(TAkpt4{FfOUqQ1i_eO*?pJsJyJ0b@}#8lW#5 zn?_SC+Ge(edVL+qqzY5K^CddXgwKpwzOHD(iX?(&2>JZcSj_S#P{uN3PJ_4I=kGKl z9e~~tP58Rw8GqCV<8idnhtCvuTP}2-y3bqBQ>g6_eX?|4vLhMeWcK)+2@|)Ov~7PM zd){Bc!>Vq49y!r#P>lv~VRUz`rmr&NbEm%5_;z#5*JFjEeq?ds119v9X}biv9ZAHZ zA&Pql7tx4lD`A@cOH9|+H?|n^YtOu%xC>t=$poq+1k{I8#U+Re7n{>hfj&*(^*^e# zp*+?RHIiSjHpEyJnnrHivW(6|qN}W+VD;+N`TnpO%lF5Fv1q>8U0|YAK_I-!Ux2ZT z`2taY0ceIW70>TXghO-PhS}8>3SwNN5rY#?Fk*;L&mtij3H2JeZB}BnWkn3pH^WRM zV!^iVgcUaukWGE~I&3J~3n4Y(1CC+|17F4N^Tvt<=1Rhrm{h-;y|2>Mtm8*xkm!!{ z_R})EJ?N~ew!XTdrM9ub7PEH+AQH2s_$;RxSvn}lm>dISdTbEYhM`Obx@oRB)?V`o4zVeF3W!1hG;HR;f*b;q9Je@5~H~N0t zmii*`CYvea_va*P`8CCgqU5mFb&pb zeSUZ(Nq2_bX-Ue0ZB66!<7==I{y=+iRWKYR-M)`zh%_CBj{!da<&|0-%%|u7;2ORD zkHPw%2XxE1-e>)fL}Q6gEdp$9B0#P(ru_}!Xh{KNWF?xCtW@qzAt`g<^Iv{@gJ=2u zqHEXd`9B8p-wkwUo&S~@Pw45twN3y1lp*bB1V?N3(}6H`_9yjGKCv_Y;qza9B$?;9 zVg5CL(DQ!`=0CiEvd{mN;@=gs;_+@6I|AL7FA+_d1hge!i?sw4C__7lk=z`9AXgAI zmaD(mLd!VsWg?ⅆU;9B3`rjl)d%*AA|W%78*HfedK=_q!OKahHtSme7Z7te5^)m z4wJ2qOig4=8wFzspa1eREPVg3>B`S?T`)hA?JjAnC7XPAm>YCSsE}lvh;>^A*K_)_ zhV8DD=uat1vrlb?t;;A+gS8p9F5}^mv~^9jaUogscB4DIfUbkaKz0l`=Rf89``lk$ zI8^)pvW)hB;vi)I57WL1=N%7$aY)_~Haj9#BIvKO{GnjPq@Am#rZTyI6BN68Vz(%U zojD^BHDnQN-beeIajW%_{(Lk`B_J870gZ2iEU; z4YaKbn?QjSRg?3G{t)Xb7 zL%j2P=j9iBOG;#hY8VLI248v>6c>k$++ZY*UJ|O&_b@OT5^v}7y)5u=YeAvMTUy{L zD)1B=xs+gFpDfX788EE}uxpC`^35=Rw8x5>9aj2QSTdb9!tz~g7T=7aUYYT@+<(O| z`t0pEL3PuOZ-ELZEM&z~Sn*ff_Y8F45|kk6LB{XBOX150x8nVy!7q3_aTH3r#x`Vf zE)c?okGetZEKvpc9O>fPisWb5c*>@~@Q7`CR>WVWg$31jCpg!TAE)P1!``Q^dun8* z20D==X`&xnO4>sIwGSgNH4 zEj0+QEiE-Lh*N`H@c!dMgOD1M>p#tZ_dn<_;#WkRCNz9Ct?H`lIM+z}8dJWH1CQYR zh7eXN>&jbNkVf7gv~4NZ{*-f6OJ#XoHPR<=y`rJKz8cwhe${P;{E3FOI!eIinjrU7 z5^gA|hP1kH;~v*&%n&4C(>_}uS|GPxs725fI-+owinw@orwY`B46pjZ=3qc zt+N40wmV?<#tpb>`Ypp-I?t2m@#3H1@s#m@#?scx?9$-Sln~zfe}9bDPv%&HwiXNk zLGl30u&}{6>K0#fb)7t@;L71g;*f6_ z`6|5&@(hD3M`+cU+QjnY2Rd9i#EK~7kbEUiIJk1;X7G-g#4dU8!Ii@g>``g5>B#dB zt{k#kQrIz*sYf1*aOIFWlR^%8gv!$quAFHimF<|xQS`c#?YEgDKK>saro6mq{oxY_ zBWtEk65`Gp8`0_5B-<{zl-ec>Fu!(wux+R=rT6VQ8+PLC9~!j#b$R_R2{jKl-oek&Fi=NFK{4|_0d(~u%96GEfk3WZ?C zuA&n*N%J|Zy|Ba6AsVAC7|9ER7s4doP3K`?5Dvkd9--qo+(e%=u~Y3FnlN-j0FCno zCX7xi*wG2m)NF$t8M0|4j`J*JX@=C(w0W0CERITY9{nTst-uXnVO^K7EO|b*BuEa;>aNzEo%>_v6(G{Wu17? z|I$jgC7$V)h?8+rVxbX2;8LSw{$KKuc(Vyd^C_fDjaNT-+NGRpd^k^J;D0+#K6c?? zq=uc4Tec;m44G*L|4S*v_he`*tBQO$|~VnKfjG#!C5{n#Fpa~PHSgLH~tY+Iq9I5P!K z!8=SFLR+EuBlqKpL`;7IxDwszkJ&hcl;~=x`PscCNSv0`1o6s%THS8?p?zoPIix(u zA{@0_C(B`1Jv3Cu3E68kmQs>Tqk$Sk;PZc()Wa5{#5*@%{J5OD*nmw$W4)YtPe~r? z#LYvX?RvAnQ#ErawL9)Bs?J?*$rLVl9z9Fxn89%C9 zmmT55_UeJ<^#42Syo>=LjEOH8-~Al~+|Z`D&JVylO%o`h{c$RRHh%`@BlMqZCWO0}DGV>hX| zZoY81rP7%q91NbjkRCJ0n{M~K^x+#5vMTk=D=cpHc#ANif~>ryo{~b4RdHdFXMQHK z8k2U>TpEDfCzD4*o~Hy;0j(`7EGqMsfC>c*|1VsEGvK3?@F5)xw6D$9>1)$>TnoK4 zu4Y?2h%pVeS)mFu9`rXg<%i3R@`{#POo12##efhZDWy`l9K}f$Bjp8xX*Wc_aCayX z6sPCgdJRv0@x1)LwH1ll0{*UwV5A=UAyh+A7u<+Ug4kx6VL*(=0_4jD1r{419pr=r z(KG2qiX-(=qY5`k`O3*?>Z_;Oh zVRj>AwTq+vNnIFQEYMm5`GXKKbQX;PjRB1TjRB1Tje#u;15W*)Ed49@yYe(?|I6U! z%-c`KTy<97)omhmU?1|+ENpYsf`hPad1QrkzqZ(a$(7UMbmGE+hDRJ#JG@)BSM%j zeZdHP>3!`%4#0*%a+I{(RF!pyj$IhsgBtaxIDVjfvR){qs=tyPs|I?2a>X;ia~ zqqzSn#xD>}M{8j|{NkZf8%>qVYa2@N9e^1rl&Bu-XbFY%o+XD^v12&M#5u%6fRV%7 z(NNi#ETirw6z9*+PrJR4=hJ^i1Ov|bPxJ7JYoB;cT6K@HdRIeLiOZ^pC_v9iQcaoQ zk_NgGSZ;*%FOxdr?xS~z%q`W;wbd=|&g<1mTh9O5>(XGMMaA~!cC`S>qx zOfo84nyTDJWmRpnNUbIl6Ib&py>?t^$3=ErjQoIr*A(F`MR-XO-ZfP>(>sdrjv~CH z2=7{}>znXkJ50EZ7G8T(ZENLXx6#~KSLX%>%Bx^i6s=j-Sl5dFlN*i;T!vUe#DuF< zhtKRzz>tYc>r_0e)^jraNmF(G& z!-N4<%yIRYVwD+d>vbEo+5=GQ?;cxp0IFahq{j&<)(0)CADOm zDya$T5n<(iJYu8y(O3s=$D`dw!gxHGuq5WF!BkphHv)UavdO80=eAD?(cO26PAW4@ zC%+RHEp)M6MJv=TR9{cjZ?<*Q^>G5aXw{*3kVwMiUX@ifMrA07-DQh!an(ez%YEp{bFg9-P^R(|1QvAWYowK?ZM?eeVuvJSVYiQ(_Cu-%1i zGEk1({Kj39zzR$zRGP4-(rJM2@`&v%+H%K{EIO}7S$SUa#FBlIn*x$ZmHa++jHwPo zDVKSC3#@Q3$geq+bduWgM_~klvm5@)`kSk2jm-Ruj+4{qry%%nMvb5J3W8c=aB?Jo z!*uoh(L$1F(%_R>$cM`}xlAxTr^#Z#-8CI<2HG8FyIBq3i<4{$UgGe{4Ow6-3R*S; zFJgA6@L`<;@kP2xn&L>u9~w#8v63U{EiUqcWW}Khxn!l!Dk{{b6e6BtbeS>C4cc1< zC#G|YE%pVQaon(p1BIR4ZTUF2S71g1kpki-dA68U(1yD|3&NOxR;++P@wZJ9Ts%Lq zCLuXrC?^UCoFq}h(UlrIz;stHJ1;zQE0uu6wVsN!7u)!tj zW1|skc$j9!LbSpT&iPOCGyKopAXJL^PeQ7cDAPRx2>>DWTKJU6kvPo0o)um$blGYM zOr3M3)EUo2XDr&?(P_ww9B?AC8}2x8*TDkew*v6phU;=STySt`73(6-kJ8eBLbU%U z{B?vois9JAj~P8d+QEm7HQt4FhVQ;(W1e3W3$LDNTTd|B9XG643|ofzMl0M4c+==a zD(tddoxQLNqdfOlv8|@k)5ULjx;xwSVt^AlPx+hu#dv!EK9tEnU*ghExLou&jHcsFKed!L{4og_rkJcT4pTx{sRi+5MG*N62Ikh zCp?>gi;#uecsZKkD~gOb#Y=VZ+%Uz%`<@y?g8}0+{?SX>ExyGkvI;PSV&qv(P~x4@ z?hwq>@b@5DHXnr}do44T7j4fY3dkd2Arv(ON4bq7 zTS=d$AN-B!_v0?7yrs=`ND>d=0rbwI(jw)k{Mi^xA#0?ho;X@3+9|L2qvy&R98fQD zQzP&Ly00ikmHC9{ALZiR+Xu1iUU@f;io;dNj2Ml6T+S{00U?Y*DLTf;B}#)>D)eXt zR!Tg!VXEd0@)2A}>q$~D-HwZd)Hxqt3#@=_OPYzrxr2)xXQ}%usy9~(XqePG6_tRy zQ%DgBu9E(c2^)vKv{Jd%_%Hik1^l>)F~eDs92GGGxoD2j5&LC)REcLHPfsv_L7EM= z9gNTBEDAguWx%e?7+5@R?Z!Jy%Jv|cH3ai8X$gi?B}uIu;(VOpu(F(=I0WmCdm|y+`Ao4)qm*KwzFX+~ghA^Wx4tKnO z+gM)iHj3wmae)N1=12#WYRgadDNq8ulJKiEANS&0_A(5r2KJLm9xh5AA`#; zAft5f<6^vXccN-pdpSs#Zc@h{ZM(7>5xK{#Ro%fD0MnWr0#n1Htbh{>S+qL_u|_VV znBFAHI`^CT`~lyjb&7O)GXYjfl#4k7F^EI3x2E}0HdtL5h!CJLX(8H@gn_F4kjUUD zOn#4Wjf9|uhC|Q-M~Xp^jbV;M`f@}1>gvM1A?O{pfvH-B4KAW9FCrvDw8Uq>w8Uq>wUr!8VlNaO>@AO{h%##v=cq5-=oSt7l zYYb=%Xbfl!Xbfl!Z2cHmFnQ6!y-v9n|6g5r&Lh*y&w0d)|A&>!%^YdyWz18ta_Q2g zPhI@%!Zoe~?)&kpIV|U6*8^)@ewY@=HPpnU#e*XI3z5BcOKAnB}!PaYk znb5P6v3urDUE|1jn@IXge|_IWdY_trD#Oqz08?g8)pb^JJX=CUadDH*r^kw8Uq>w|7i?llNYk>FYKB{ zUeNms|LKv`^=k}h3}_5!3}_6D9tN_>3wsPu+2lq##MEE-)&O#&wmbEQJ^fB&Kx05-Kx05-Kx05-VC%y`Ho0-n0lKne z0J%}y8wWt3eyuT}F`zM^F`zM^F`zN9m0=*8+&GZEv2y^qQQI50GD6VJ)fmti&=}Ad z&=}Ad&=?p11KH%pf$WX(0pvz)ZyW&e`nAS@#(>6v#(>6v#(>7aR)&FWa^pbu#&ZXd z8@0W0DfsOikH&z;z*dF<{*|lb{HMlmdPsEP?lXk$D8DczPn zkFF9KIY&7Bwp*UnT*F3)jx+|o))?SiBRT)6tDAl!I{fC@nYqRr&qV{AYplK2#tg|3 z8*a>9tht7b5FKd@e62CSxkhsSQ&%@_5FIwo&&)M`dEn#+IM-;n!p017jkjNX<}J-N zY=r1YW8iCz0nRm&^Pjr9=?|jA%P!8$HD>J+{FHNz59isKA+C{sL83u(4I3dk(ir$! cV}Nswh($ literal 0 HcmV?d00001 diff --git a/test_aodntools/timeseries_products/test_velocity_hourly_timeseries.py b/test_aodntools/timeseries_products/test_velocity_hourly_timeseries.py index 5f734b7..9c3ee4c 100644 --- a/test_aodntools/timeseries_products/test_velocity_hourly_timeseries.py +++ b/test_aodntools/timeseries_products/test_velocity_hourly_timeseries.py @@ -75,6 +75,16 @@ def test_all_rejected(self): self.assertRaises(NoInputFilesError, velocity_hourly_aggregated, [BAD_FILE], 'NRSROT', input_dir=TEST_ROOT, output_dir='/tmp') + def test_size1_dimensions(self): + input_files = [ + 'IMOS_ANMN-NRS_ADCP_LAT_LON_DIMS.nc', + 'IMOS_ANMN-NRS_ADCP_SINGLE_TIMESTAMP.nc' + ] + output_file, bad_files = velocity_hourly_aggregated(input_files, 'NRSROT', + input_dir=TEST_ROOT, output_dir='/tmp') + + self.assertEqual(0, len(bad_files)) + if __name__ == '__main__': unittest.main()