From 1afe2c140411006c4fa0e8ac97ed7c1656ca5ad3 Mon Sep 17 00:00:00 2001
From: Alex Guinman <alex@guinman.id.au>
Date: Wed, 13 May 2020 19:51:27 +1000
Subject: [PATCH 1/4] Move csv export to function

---
 nemreader/outputs.py | 33 +++++++++++++++------------------
 1 file changed, 15 insertions(+), 18 deletions(-)

diff --git a/nemreader/outputs.py b/nemreader/outputs.py
index 9e4e8a8..e0fbaf3 100644
--- a/nemreader/outputs.py
+++ b/nemreader/outputs.py
@@ -94,19 +94,24 @@ def output_as_csv(file_name, output_dir="."):
         last_date = rows[-1][1]
         output_file = "{}_{}_transposed.csv".format(nmi, last_date.strftime("%Y%m%d"))
         output_path = output_dir / output_file
-        with open(output_path, "w", newline="") as csvfile:
-            cwriter = csv.writer(
-                csvfile, delimiter=",", quotechar='"', quoting=csv.QUOTE_MINIMAL
-            )
-            cwriter.writerow(headings)
-            for row in rows:
-                cwriter.writerow(row)
-
-        log.debug("Created %s", output_path)
+        save_to_csv(headings, rows, output_path)
         output_paths.append(output_path)
     return output_paths
 
 
+def save_to_csv(headings: List[str], rows: List[list], output_path):
+    """ save data to csv file """
+    with open(output_path, "w", newline="") as csvfile:
+        cwriter = csv.writer(
+            csvfile, delimiter=",", quotechar='"', quoting=csv.QUOTE_MINIMAL
+        )
+        cwriter.writerow(headings)
+        for row in rows:
+            cwriter.writerow(row)
+    log.debug("Created %s", output_path)
+    return output_path
+
+
 def flatten_and_group_rows(
     nmi: str,
     nmi_transactions: Dict[str, list],
@@ -185,14 +190,6 @@ def output_as_daily_csv(file_name, output_dir="."):
         for row in rows:
             all_rows.append(row)
 
-    with open(output_path, "w", newline="") as csvfile:
-        cwriter = csv.writer(
-            csvfile, delimiter=",", quotechar='"', quoting=csv.QUOTE_MINIMAL
-        )
-        cwriter.writerow(headings)
-        for row in all_rows:
-            cwriter.writerow(row)
-
-    log.debug("Created %s", output_path)
+    save_to_csv(headings, all_rows, output_path)
 
     return output_path

From 8049da2cd53267b978754ffc27a20d22123fadbe Mon Sep 17 00:00:00 2001
From: Alex Guinman <alex@guinman.id.au>
Date: Tue, 26 May 2020 19:38:21 +1000
Subject: [PATCH 2/4] Add more examples to readme

---
 README.md             |  33 +++++++++++++++++++++++++++++++++
 docs/plot_cal.png     | Bin 0 -> 10717 bytes
 docs/plot_profile.png | Bin 0 -> 14205 bytes
 3 files changed, 33 insertions(+)
 create mode 100644 docs/plot_cal.png
 create mode 100644 docs/plot_profile.png

diff --git a/README.md b/README.md
index fb8b42b..ffdc62a 100644
--- a/README.md
+++ b/README.md
@@ -51,3 +51,36 @@ period_start,period_end,E1,Q1,quality_method
 2004-02-01 00:30:00,2004-02-01 01:00:00,1.111,2.222,A
 ...
 ```
+
+## Charting
+
+You can easily chart the usage data using pandas:
+
+```python
+import matplotlib.pyplot as plt
+from nemreader import output_as_data_frames
+
+# Setup Pandas DataFrame
+dfs = output_as_data_frames("examples/nem12/NEM12#000000000000002#CNRGYMDP#NEMMCO.zip")
+nmi, df = dfs[0] # Return data for first NMI in file
+df.set_index("period_start", inplace=True)
+
+# Chart time of day profile
+hourly = df.groupby([(df.index.hour)]).sum()
+plot = hourly.plot(title=nmi, kind="bar", y=["E1"])
+plt.show()
+```
+
+!["Time of day plot"](docs/plot_profile.png)
+
+Or even generate a calendar with daily usage totals:
+
+```python
+# Chart daily usage calendar
+import pandas as pd
+import calmap
+plot = calmap.calendarplot(pd.Series(df.E1), daylabels="MTWTFSS")
+plt.show()
+```
+
+!["Calendar Plot"](docs/plot_cal.png)
diff --git a/docs/plot_cal.png b/docs/plot_cal.png
new file mode 100644
index 0000000000000000000000000000000000000000..71072554f1f3402ad55fd6cd9dcf856dfbf481c9
GIT binary patch
literal 10717
zcmcI~cRbbq`#1TZk`W>$%9d<Vwz9Wmb12y>BbyT{D-<%4Ep+Vcb?mH=l@+pOosPZT
z*ZF*ZkKg^c@B4B8asP4G!*RUNd%mvgdOn}m^E#oL>I$S6uU*8$!y{EvlzoJUhhGEl
z5$B2ERY#})GyM9;N=8iv53eMO<oF38{C(40@sS!Ho;MpF-pf~bc!#j*<uV@L3qCx&
z6*D|M(F8m^+UK8Yw8Y^9qNgefvUq2>|7rC(F|dW$O-W6Tc%JYYEk0g_hz1EB9{Zw_
ztjuH2@%3acw-3YBXHu9%!QfY~D4HqGXHjr0e~=-Rl};ZiL>L(sSC;G5r+q-x>lau2
zEN?Z;smZ4<zs)_hTEo3+>tL528`phvM6vnitMs%0X-k3Nc@}mnUnFvWqP=D2ty8a(
z`s```6@$gjqm3EAUazV8y{$<jEgc;>ruT|0(NDr0i@9#XJ4z+-_W|m-Z7tX6mRTr6
zWd-v8A2vM9cgjWGW{4t=Y7hF6!4(#$p7StgFrA~ztkM8?(o|-qi*6{?tea@jesaHc
z-EC_xafVOo?2eUUzx~`qZx99H*CTidQe<IgqE{gV$=wT^rwth#S0lc0oqk()pL>4q
ze0gR<`&`q^b2GlH=ZKU;al1ZMQ~kNzs>n0oAZWrfawR|;cdcm8%Q1IFD<>k_=c;Et
z!ZM8cpK9yugh?;ERR37rS#}d(l*jFc&72NU*BQ?HIe&?WfGB7tOnM3H`MlbAs^Q8=
zr!+h6M9-^?c_JC5BHQO?KH3#FSqsfRzjx37_?~xtib1Jm`p7B)Vgkc4C-rD>b?0h?
zywjV{``DBEMiD`6fig{<EAf+&rI;?ps87*xiY(5JmZjM5f_=r7J=a!NR`7E3^2~ca
z-WA{b<vMC#kv8(p7(sD6^O54E<DC0PK`0b{jG3%i`NX=61JnI)jPE%z+p^($g)(ox
z;wSW2{&K8jOUjKxHmtabs<OpcCZSNH>$`aqHUj+y+TC$F9y4V-UPz3)sKN1%lRI^1
z(|=Z1Yga`12Rv#g^(?S(lBq-Ml0R-WFUK0gGR%aN?oJ%e`Sc$ynC;)Wbfe+YQYq7f
zQMg0yfdyYu=d%Y|)II+9uTFnPdx|rBEiJXNv&*cj`_$b0psTB^qpORFi;HVM?F+2u
z+tISO-=E(7NFpRK^v4aMW_mM;LrVV#sRfk<k)1}G-NP{^PQmEIe@gAjU*8_IsnqsZ
zp^ltFU@*q_s(!g~rXMQ&644ZEBT1x5!uVzGNO1deu27+p!!oSDwq;He60FD1j+?Cy
z*<QcYaA}eFDGj$_5=P?0uWw|!d~|jp{*)uEdEm1OCl{AFcD9XP)GdTnB~DXMk5(YX
z`eO~X6DLo6w~>L&J>qS@!FQX18W)Di4IP||1TTti>z#14*q0{DV*8tN6*9{Fa)<0V
zDp|EwE|cl#B#OScO%wjQz^F1hHnx=({r7co5L?f*GEY27a7)0XV*^uF_2fq0n;s@?
zt7+g(4`oOmN~Nq^l(MbuZP%q<i|#mXW8}$BR#w)ldY;8jjt9mYM`E|w_85Ce+uo|2
z?S5NHDXD0yXVi3_t6OOgsy!3YxJmRiifY|}-8Mr%<0n^~LumWn+c<BQYfOD}i%Cw4
z2TQ3|{TXujf@de8o)NtsCpV+jv5a2`4O`#pq(M63FGb2#SQUF&JR#TJlNQ!m(i|bI
zajFtHGN;<v#JFmWSUWs9BIaYy(fqniW9*m8{kztLNB+tc*2iZdM?{<Dr=z}pEQZr0
zmQ)smwkjM+PJZR%L<78jN4nwMv1i*k<V0VQOxPmL3Nt^7Ubb<%?|WtA5P-my6<f%O
zTofv!u8v+uM@dO3y}CN_J+p#<`L|1#-u#|X)ndN1z|?ra<4$BWDs=Tyf+`#4EZjyt
z08zAOiKq02&%e*aV)^@y3QKf)%S=n~Scy6D<mBY>x)Tx?)A{-N=Pjfn#h%#N;-Ksh
z1-DbN<i=q-rCHOn2YyAQENmF^M+m;Zn}MELw&ZO>j-0%T^bN6_fiv8^2@*cSrB;2}
z6%}zGKVD<-{EZ(L9v&L&-$~*dypzs^dfx2GnyTp4)R^tu&a7k9Mbez}A=oY_me$wZ
zq1av7tLaA^tu^anAXfTSQu|q*AoT)aD(!zo358^}t$*WUMERKk`*iuUgfb36DG5Kq
z{V8##$F8;fal?9M&5+4?H7&iCRl>$?2DTRI<xfO&pyL(KbIO$MM2P+f4Q%Zf?^h(<
z_O9P!kJ~RI538P*whKm!Ra>i(ki{!=b?uO@rQ4&CnQd%`jc58<W|V=N!e=#ym4lUM
zW#h+~iD@g!4X;r5L(J?TqfznpUqxTucXf?>iT|vm2E3!a^vd*bZ(X{BWOk<X)Pg+A
z^zcwGOL+R|&89>_g=aU(KM<|&AMbK8z2CpG7ycN5mwKc4n))Z&uhMd%UW=b`fcr+I
zJ&JMZSG8-te#tpGt<iEjT2fNd<B5io>!Qh#en-3bxX9!cMtI7yKl*ayv-wi7Qb%X!
zVu1<r{9|^kB~{nYQ~-f<Y71|^Hd#428QN!<W^xKrPCUn?>i-g9IH6ib2VqrbhcVc{
z-7`Z{C-qJl`}X-8T*YJd#~U6E=8*IHcW9=3FcOjU_nY+^4t4I(kT^A*2)izRzeRDG
zSTw%-iq>o2K{dhvQG_aT*W;S@vgNtZpFuWKuhJ&h-HC^K&y#AVWm-NMQ)Qz9+KKOe
z@xh%fR!xig4o$da(vNTZQjGP8s4@GW_FJCf_(y%#LQnaw3XQ#$&-75KQgpzh31`jK
zW&fnM{NZ>9A(~7U-j>m6DV&5abc87LA7f0tGLT~J=gxOBfb5nwHA?HN^F%bE$zjQx
z8y*}W_)JpA?B$C_6xMyhGTz0;QV-Frh=#YFqj=uPbgYusGl#X_tbgwRM@38PGEq=K
zSC_Jji%YTo9o1CHWJq2Cn*kO>L&F3~Khb~?jvs{MDykuO{FaJ?Oh{6b`FKQ?e+g<?
z5I?kV*0sCNR8G1hOJ{l`mc~eYoTlhic7XB4WsNgg7l$e<O7aH!Wm&<3k=DVpGJkcV
z*{zhgIC+16$qyesZ13#6kB)ABODA~GdG;Tu<dM}xaW9ufJ7zLBH0foh{P(`azO=e+
z8wBVk$26%l$aE$;Z-|QJ_k<<|Mn73QJ1*^P)pP(mhlz^*$;qB9rI@SyPi5f@Kl=}e
zYBX~~f}5qNxV5GZS>fR5eQH`-SM02u<l!=Nwpx-r0>MBQ95h+$AqYWD8<8s?`HT_D
zisMu`dq1N|-zi7=DU<YfsL(J|b+l~1fYJ?<*%1K93!l{@6&NplkmSC7yLsmMj~_o8
z%WMYgod3*s#vBhyooN491r35uRCwu;uKjkcnefi7J9|>Gccbl?za;frU#iM6<4D>c
zu@)k7aJmWXLonuh#G1gZ>2*@llrz}bw3|L?OLe`8XuRXf4||d!OuBsr52WujDbd5m
zo^R#juS~)*HrVdxlGl5AnXzNiC*zDs%MAg4Ri*`qyR77@v42%=wG2EZ**>8P?b~gq
z;ccMC5!Cub<s8;C_sf^2IqctGu6_K-Q(xUe<4#z94i1i+S1+$F(~~JQs-`H$Xq}5#
zRAtjvjn;aw-bL9(vu0ROq1g2{m>aR9ROZfRC`NRMez>04ajI0;r&rhYcF7N)2T|AY
zXE4Xc(Q$F4;q<P@w4cbq!GY?hd+%pwXNw84LvnaXiCT2T9D=VSf|+#1aGYgoYWjm3
zV<ybWyWCKBQ*-)FDXEB8fbKLrGd7ymdS}e@$b#d7%ssCF<n)`b7tX@U`!1ecQ}mL#
z_eX|hGf)$qGfd;;)%5BDRRZ6WbLHjbYrm>v2uWyUU#M!dc5*zmz80?C)J7*C`du@z
z^yfGcC+}5buO{ijl1pymzFB{;Cthr)yRtiD`rRZ?XT%=rASjehf-Vmu;-7BlEze|(
zv9@`f>N}pX;uKl{(~4Sh+btBDJ?sc=5rXs^*=WsJdw#!38Yl4{68E%k?_1f)eNDSZ
z4}9N88!Eucyg(RgTBKhg*F|z}uBa_CHX)f$QXk(Tx059B16?l-k25>M=2yT9^Ge1r
z8|=%@`(wo?bT%HHTB(0%&oX`G5rP8gh41_EpuxPIOg&YyC$9*ttgcQrz5h*wNm)5`
zC4m7Ms6>43g`+sZlZ*}$BCehLDSL9HWn_8M=RTX0xBD~RMvFW?wBj4TAG`7~*Y9`t
zLU+7Dtvh#HTbr@R<O@)pi{BF;LK)2b5XzvMb(-F(*8H$zk~vu8r_YyFchkSNbq@P-
z@}izw>P2}HR1?br;Va-KwP6PJ&2t5mf$1!~=h<#uwe|}5+pVbn>#FVMIT>Z;i#R1R
z_SHi6;bRcKI>tQlA{}zx5sHTRbq4XF-nIcw^}(X=zhi^d3iRsFsyCgO>dwpr6X1-i
zyoM8Wwu5i<`*ODxQUrd@6PBQJbiHKiSRRHl@oG}PE-d5`cAEOWyFNj8*Oq#3(u1+O
zy81~`#LDksX(y+Crob%E@MzFPp)B+$6i!105|ZdwJ--DA_xA1Ekiv+2164Qy)pGKs
z7?i)d)qnVtEVEw`a8S_ld&o*nE*WCyM#`bx+yk|wkS{7F*~F;@TcqXEdj~PJx%j72
z{MoF;fwE&wvRx!hd3qT??~flUd&!vjl4By@KtqA)d2KB|ev_5P%`eCVIYKUrDv)YX
zLr`=PVU_0OwHmEwP>tClX%rKW*PE2*e%c<1cal){@ILG$X)f2!Q|lrjGDVn`=vSTW
zP9)xQ&Jb~zt+}8iyLs;H?ME{@{`bFqE(Ris+MfxKifL8<+nR)<P_c1wn>kV^p^=ef
z1A~JH<0I&A#|K;H9^Nu{xbj|0uLeC`5i(5jT)nLQ^)y5)lqM;I1Ah;n)FP-c=b`LV
zs=b|5!(E3Q-K4Eeixu~^*13sb26BTDL`z7b>7{W@hGbKUQd-5UuWhr5t(nM|?@q6D
z(X5~Yjs*~s*9)hq`gzN@qkdw-(&yu66JLD|J-pN^_($XHpRVaDw?CA;=zXR<-GbSg
z)xu0PI;E^kh)HJf!ki(clB~kMs7h|R#F<^Px|L>qG))qsh#XZ<=^Z02%zuv5k}8r$
zRe`pH$nyrLq0ys9lz8enf@Xgy3~z7?Pkd&J6hA({F@-4VnR~w377-=*bdi=|H=rXM
zjx)U$;O5q3pJt-KzODR#W-dH~gN2pV=4&XkM&LpGj}5}o1NE6B$t5#?x-Ozro9MFf
zH_@gjl(&x2E7J2TK24I9Gxnj2@j$?4E@@8E7cHqysQycLqVFZEd;9p%-FIQRfB*hl
zMk&eJ#I_KPtF(+njXvs$PQ_c(_-&?@mlb_PhA6dUkDf{lUt94vbl^u_4MRlD&Z&m6
z?|#gewY)hYn%gHT(?#F46ZXQ0vTcv@b$_OUYx@n4VvBAX|D%a%0q5D)kV?F8j)8v+
ze^WPg|8aw+*7+X&J$kY=cQN<c$)hu}VCMOWL|hWE%<x@{`1NAMq;9i|>ldJ}o15F;
z*>GMW`pbt@*mqrKe64Hr_utzHQcp4c+jH-auDhf1^6DpEPI}Ho2#JU&WMLw2)5-nA
zUSLRcRn+JuY4A<@n0U%$vNo$QgMPo7^exR>8AEfohEimo6zBY!E&qen?2jR33Eezc
zUR#@wZ#ZWCYeWB~sA$RbbiuT4b>Y@T+QRm^>ANDPmpyYMS1{SV0q64cc?WGtUaSOK
zTa{Iq?UCwbtLF^bll+uwE*xI07_*jgzBue5W2~b@ClVK1Ldd_qcM6$hMOow@6>LiO
zrte!G@|8&*_oyH7VIqrd1AADX{Ge3Rb~Y|5Q{JoGQ=W6AhDJu5w{Oq(D=;oV|H4N9
zah?#(<S?x!apIwVVY^88!HYM#*C_-n^z(X5%@}(Zj~dHI19vny6vA4wMOHdFJL)!b
z8FXI1WqBNk9JBtj_q;Ewq0jE-kpu&Drjiq>r<bqGw6^dPBqsnf7Fk7yoG|p#<U0}P
zAv2UV>>L~{MhZ>xO&e0iO0Bewj8s5G$tFc^8!6D{Nz+_>uTHP}AgzrpGA=|HaWhgp
zw7l*nC#O=WnpSL8mh+%pn+clam*j%C__F19h2*#!0y<(p`dOHC#76lgn(Tj^2|4Br
z(g{Tvg{!H?w+`^OX$dyX6T$}#LK~h3tgRY|-xyPLr$du3Z>*^;*5E3A19dwyVKl|7
zQv6$f8hL7Yrt+oOkikD$wY5pmD@7+KtI5d7*e?*G_c2C)kW<DS$0NJ1h~~D`XS+!D
z-6YhpUWH*7Iat2s`yDIALSM_Pc<kL)!yP*76WYAopJ`C-!dg_nd+p1YFVt7Alss*_
zn5~*fxjUZw%KkwWHw(e+D8ntT^DOUn<j&KzMg4q5gL6}S(YR!F?UTSIYd)g><Kp7t
zXg7)MJJ)Js8S>w*qRC(}yfmu8;kazZmD<S-lla=x!&L(Say2!zk`J`9-VDj@bJ*T7
z@&%jDbiP<diozA**L^*8pIJCU1s1YDPXTw{2;CGH6B9dV<Mp4Jv@C{htjKqPt0aat
z<Qc{E`;YqZL4A=ESk@m}lBbyxnf+HcL&#0feYmxfHcZ1IrZi*{OIM)gs~>35L)aa1
zQ8!i59=pP5tUeT(ts%%Q<w_JgbV;%Ed<~slu3B>Ww#J<lCF4VO^%I5y#?}mESu7{z
z!)A*K@Dd1EZT=}es1S_h9^D(&-0GIajwv49d}Pc+yVRGu1SSIHu6)F`wpdQ>y3U~V
zko^P(cw&bvBn8!A+-a@R(U}SIT=Ks~ZWv-4q5d1WoSMOrpqiamu|qkS;3E330ttGZ
zJx35oD)L&sV`Kw~s`__-f>5?5&)hT%zN>ZGiGiNPV%$q}Z6vp<DI7@I!o;8d#&fSt
z{;U$eYgeN8tm#l4KYY5)mTJQ+#Nqes{gZqYYVY81gKfs8x^OhzijQLOZ_s<UD}`(N
z2YUBp1IMK$i1}6Kv7TUc5!BqTE4nvjE%e5B1AC~xdG<YY{p#G<8Pf!Vq*0`D>^~nA
zL-#rNC5xr4x(}G--O;aV<++5DB9>B-DFno1mt>KzTOBJ)o#*6r3k=)*_9olw_oh?~
z3=A3@B)qpd82$I3Gz2yb{+SPDPx;*|hRG@{Bzw8B;mF{<d<)lJ^O&OPM=Bhr>Jl;j
zd$rEoIqB&EwHR~~<Jo>YqW(KnA}&vP>$eBeda3KSdd2KPoiuj8sGY=R;PF8plI&1~
z(Qo_m=g*&=*BSGMyZZch1%>TLDW2=S*EcU68?A7-PA;)cVK-8^P}0sYH}pl@nU^-L
zq~v`(!aKCwZp7v6WN!;nTFc1j`d_CTWW>C6f6u|$_>0<CiyC^nK5*UlI~2mR>isnH
zUV$Oh=*+{zW2=KX>Ti#ef6Ni4AvKuMKVbyn6A>2|R~1Xk%OvSJj1+yOql4*XJ9Xlo
zRWqvdYNv{7fB53X-N9_NOcW~Y<Y;ejWTgFjqF|z!2M;vG^*(#hKvsETnLzD;&=Yi-
z@2DDrn|M9p)R0`{zmo@g@6n@2_uQ5#6e8(wadSr`IiV8+1JB{X0)DS>m?$<ZxPz!m
z=!{`k&`FtC&n_s4m}`%6T^r@Uc=2Lxe*XOLrXVQuC*_9sy*63cHBwkOIYW~iYeGY*
z*j)}bEq2$&9?HuT4G$0BF{z1lp8F=J<L=<V4Xqxq&F$N_aY~F`{Zks0A{s6EF)7KQ
z%=)r+j{0ZV&=O1@1Dzr4ZhCe&qNJe_l}lecsSYh3lu1ra%}2N^*X19#PO(u)z85Qx
z$0sHZeg~1~78J0Gh{PK?R=?ifcH7;URC)5`mYkei)uf|4W=hT3ng8|c*Ze+H+}zw{
z<>gwsx^ECgR&byC1-iSK={{1CU((>(!A;?%6&6NbcdQCX6twy_?@)zCGlaJ_qy3IL
z`ucLJsuI3^lLv>km%H-UxmGFe<F&C;AkD$y;Wo&wgIU^!jsODUj^19&!|gw1Wo4R%
zhFuPolW#IO)}|X$qohs~z{C9IC<4kDZdYV>6INsY1cU55+*zfQ@)y5*_ik2pHn_f?
zkVpsnb5V2>-a;zzJgls&ugN9%@q$}q%*~ncOb-X*Gj41A7#wUV?UM|jShs9>L%O`<
z?SHarx;2=38k>4{yfkVrj0*&0;$V&jJ1ou(G56=YYX-NuxDE#5(X9a6<^xL(4i5BU
z?r)`FJ$+|LOw%eW<58&F;tb?EIy!5UwY2fwZ6F~zIXTnw^DmW-eJ>wDkVxCuaCUZf
z8oTz~z5HRxrP|ubX~}5%kVa(Ejf$4X*bLmq$2vMCP{*R8qNAl&<hpqu?z=8N@3}jY
zHnRGqu+VZE?dRs<AqW;7II7e!yLPOoAU(g!+UJ{EXJ)@w2P6zo1XQ5WQ(|RPUoA+W
z@NeOrZAqQAsPNXm>RxyXhucrq+%0@jGO-R31RWGi&bTNAyC|uv<NUUsB$2_%$p~2X
ziOQ(BI2v19ThEQ^C82s&gF7@dcl4zV4Cws){BWOoZ7se}Nm19(AdiTMc=qhskI~T%
z65cv#D=YTo<Ye)kA#M9<bn^WCGc|3eV6qz-X=yLTHmZWoPNvV~bacA9yAL`w{2xF)
z#`BoW{>e$r2Ky1hi+g%y6%{2FmGBi3eH|V7NczrvgVNCkq!hS|pTxbQuRGOUj_H*3
z@ezaTyuHdFxVh<qK%TsKtg`^2sXN-J#e-GEf_RZH3wj_FLR$k4MX3h{@msr~4|y_z
z>>nDEwXk4~l04*q(!wddwbH(Up&{y+uX47#Vwk!X>|8rO4hV2JtEQ)?_bcdA`&>FW
z>CmE=`X0DIegsXdF9OPfy-vUHCWL;mid#)g!<i(G&QYkZPScVM$YbYg*REyU*0%N{
zB_e>x1K0#ZB2W44nYgcy%R#LJkjzAgO#O{KMR?8-+6!N|u&}V7sxycohiA~z(rQ{6
z&cnId;L|*TDa0o+k4QlLzc8P84zL1fWNmBvKu0G!CZ-eYeGw7bm`>_af01hucTCZr
zkjHb3pG04DiAjtP4YfkXtH5&KxN$@L?0D<3uC5DMIXKk^R4O4Mal9zf02c856#+S;
z#opv@0&?-j67zNrF0Q!`dWiJTp9vw00R_3BjMe@8Bs3&;sF!+^TkM6kO--ei?_ayD
zaEsrYB+K-|8w`G4ll$Y)D*&p|(UkB|5GX(}2y3E<69Wuo+QMl|y?1Ou>onUH)yu;Q
zw`PVZYEcs=^?YMI{YxBqL+o3}d*(0s`7B~$*TdT$z)FdEY|PKLg)h5_;_NF}`ta_B
z*`w(*$%Tal--B+`NR%qp(im0@=%DX439Wvy8Gfw?hD|~u86asB+s-IxJH!FY-o~El
z`@qqWhgtsZ@1*HtHb}aNn&kTWdguaOmXbYi)?J0uv^u0nScc}KN6iod{-mCKliIj}
zfky&@f;eok`<egbdjkK*l=bG<B)Bog+35)$c(UP(i<X;PTey(Dc<};fg8y#B*DXTy
z6jV-iE%qc;Lq|tU&X84H{0@?=ecS8kXWyQMOnw9k)dlWb6_9IZN5^lm=>sw#yIZ`x
zyp|S6y}AE*--@7JUtfRx_;Io0r083DW?Wkcmd%vfJ2}7By%qp{y|_CD%_n=)8gg=k
z3Q-IZHIQvqUm;tC9mZv${flAOums+&+ne$MAJ1Gk2`*#?h%lu!AI++=52me(j?R_U
zp9PuO*>5s3GIF~Lii+5Q<+L*8y^tbLTh9v!2-GoGt0f7)PflL%QQ)|J`zhGGTF8`^
z8TVRF9v(KJGjI~IABubl>k)lTLB!6rlc41pAakOS$L3`8_Zo5Q(3>pw1?<YwQeb2%
ztToTov9U22%blL?HMq{UUH}pEy6tItUS2p*a&&xrFbXBU*}<HZljAcoNE`zm!8?k5
z9@Vhv;94FU;)*vdGtd8ScJsRa8DNEs3WT)ST<8V@a}!PWsNoY79IVwb$iM95f-}{)
zbn_b*T%7wP_I|sG);oD~=@qg_NijyXw<PeJ*VQIOUFCnUWvJ)vEea;?BEZsOx9O>X
zmzNjhU2m#XpOgRA^>1(K9$H#vS8BMpx&lryOG%k(dtJG5B{n`D6h&~q^*rP*y|^dE
z%ep%8N7~x`f83r~Sq+YjDOp;wB_<|*%#M(4prN60Io!7K$A+8!1tLq$&Su;8mhjoz
zSXi%|wj3|B+1}k97#UIU@DKsFDeLoR?CBB4b$eM0qD<LOcG2+f*`<(p86BO@o*rO1
zLO?`_F=Vkbmr;m2@vO?i-<cMiZwcXo=-dY#2)Bs9@)>igZmg_q679czd+F=fuNnSW
z=m{2lrjHnKE(35DFvJI(ECa6yL;bn)8BA<;OhHhzoAG&e(}Qiq7Naz1b%&WbXle?P
zuba2$0gQk^b_4!EvF@)eCz>slwuRwD6X|tfO^^&tW@!-2{QNsG0$_^?{q%`mR7{Lc
z%>52noVZ)u-F1fv&r!WS^RDQ|+FI+nyAQ!Kpkhaz*Vg`WW>iOKCpN{5KjDITV6N@W
zz__bpaxP@=D$~x2>RoFm`88Cvs9&d~FyICaH7`~Z?%L5<^`+d1Rm*I1zG_)Y2iXUL
zQ%mOqXxi)GRQdSS?mhHB-ggEC1>;tns{%J`qMZnMg;Qe<#|yl749hx>h+do?dwF|%
z?@#Ta6<l2f#oX8a?uKI2FZ6%4>SKuGHSLA=qknZLHa7M-NGmqgghnVioJ9<vsIA=s
zX&7*7!G=O9L*Pkh?xcYXJJ?#<+Ni~3pl$>8S|0CX#rJ;);f7|JWo6VLP{IBiF-1x=
z98F>rm6e|Y(>QxgWjNxH*0$UQ>_)Eav^hFAYPnkgY%m@J*2z&%&dAH-a$L9oNDh#<
z+j|}~4eoBB0s*;(hKFZh9uATdupb7mqb24z{n;JQiyIIP{roAvjs%1*tEhOSso4lZ
zr5ne4IPw8L1qCy+YI-KBo)ZI04A~6Q;nO|mb9uwcfHOStDuadwVEz61^9Sc1-?jUh
z#&HpOb~Ih0f9IF?YJvSkMFc2Rev2+@SV+F49vCCm)Yi@ooPoa>AK#9H10Zx6n7___
zi@>3%kdV;3PoKWSI@Qxv;+6taw?&3RxlH@*$xUL=GPUL9AE1)sp|)^<he4kCIM%KU
ze~V(&P*5NSl67-;-)tgg)Pw*)O@k(P?hFlill!nBfuVo1jdd{w1_9rv#Rmrl5(KTt
zj4JHY0gR#ipcr1}<#@~{x?_-mhn46ZZUD|!a^(I)NLxTCT%fLAB04-g+=q@wU0a(9
z9&7<@xxvQ_jr-U+uhoJ|ctGF%nP3159Ji%N_>5LLGsCgaqhhA*{JOz1ngduqZbu#m
zi+JMfeVD@x9a)^^nVy6l9{w9Nh{=85{o+Mw_LTVcxhDYGxMcKUOH;%ft9ii%Vt8NY
zq~ZC@1BCy91tcLu=`%beh@b2v$BeSFSjLlO1-L#af_h#`^bU5~I=x{5zI_nKZS)O>
zHrDp`S*4{hFwk<6YPGwz6{8`U4g~-`n0F4otE;QmW+(f>67tjvD6kjQwKIItV(ilu
zE?u1F_#4D{+;}8#b=|5M)V`jn-<6ooS16P*+BXpw`N2P)TRkG2yu5onPkh(yOqqcX
zfL*bzR2rHlCOshNvu(!jv9kvQE);`ef)4SedJg3ESKu@t$0ZZ<K4<|3vXpRE98N3&
zS4BaRqfgeI6q47jT$;vwg&92-dRW}p2;kwK%lvr_aqK_PvH`w8LY!-a^*dPLMRtZ<
zW=^>8YWaA9Gwgwr6K_$&sc5bHdb4FIK8`|_m6dY~3vmd~#>`A0r&To<DV89GlmzCq
zfJ}U5Y5B<5`1^Rd9r24pfeh|O%0!UE{;&YgfjjtocI_n(=BZ#Z33v{p7#PR_-2+3#
zBqoOTSY|?J3+-E*ui`q4jPmlBfkiJf%X=cH$H$e`)JPekf`9!ohNjTk!J*Bv6dFP>
zCRINQajx8b2p9nB3Jf8=SgUeF8wUqgDXCPrHYmxCo}RV^m-jI-bAGD~V(#`sBn#jI
zK-$2W>a+uxwswF!a9!tco}PKRA=uwcaC^`&433X?!I&q*9c)kahYtyW&!)amV1Fg6
zvx5t3T=hHvz#d0ypytinBirnS2Vr@G$R#_>qa+r4M8g+K{UI+vF#Ih%f%fcOQd0Qf
z9~^N)V+B<Q$H?UjI)k3XjH@G;*Vh+jf*BX%j2xe|g<S!q`Lx^+3bSd!MN3mt8iYH*
z#@4nQ2LJ&KyexlJIV0%acQ0fBJc6k<Fg2wP`@moVtlEntBv?=JDen|1WQx@B47vEx
zS}Bwb|Evp`1fYLFz2l{-D@4)>63J=h78T``lw1bs0fn!qqA~|wLncf&p)JLSWd%fQ
z$Sz7sO8Qbzu;X78f(Ajg0P<~ndmA38yUKZPd3ALT7!dk=HjNZjXns@=(Qv|K^(7PD
z{4Zb9!S%Bl%!Y7*d}&^H7ktwE4<;OqLg5}!=vrp4D6B8&_Wqwgp>t(nXNQ5vRZP8u
z@O9|!8g0v?<Kn&nP(ZKs4X}KsnGmp~qr2ODZM1~nra||OTUsqN9Y&RopCN|qDV~9!
zK7DE`99a4YFo?r?|J~A5Qg@wuFIK2|8;%7s<s)X{3Bzp`jUc}pphlpVCSM&X8WEx#
z0kKWugg#))(MVnFWKVr-6sVRDUk92S)-#3n3xz^eZPct*L4BDY8vB;O(0*A#>giOS
z*DvSx8$px4oi_!u2u#2|fa9}^i}@_I2`JZxhK9^a$2fR+7}ijlMn>I0c0jgqi!2Fd
z5UNN92`H+IAVvsk1eL)R0KGB4u%Pw$vC#6E{d=av{LIWC+!svz4&4Bq=5)a#S)7ZK
zViOb;<jFMb{P9$gSoGg-i#7A_Cq4YnWJ%#!S`7G@oyGD*|9*DK?*Ge%Q^qIC+`C?n
z3OnFSfMyo@N|s7WxAD%y9|U-m_yl<P@DBIGkADu20)7j?k4^AD|89E#?{WXMt)1LF
zEbQ%X^6?1r@!Z2Zhj%XR!oPcXoP+cH*ZVjB*%@z%CvFtpvbxFXyFGj2W+7_sVgbMK
z1o#B_xcTmK3kW~v7Zl~=7v<yQ;^%``g<q2<|H}amPS0#Cz5ee92nSfqzyXxFD_D5i
zSc=NQi-x78tA!h$u&{`S)f4#i>VJQ#Y2#^OuP<xkVBzR$<K&1ZEXdEOcenok+{fU5
t?!zxAwC5+>1pD0h@AGIoS-N>Vak2Q{*W-N-@4qXhB&RN0_`odie*nMKLKXl3

literal 0
HcmV?d00001

diff --git a/docs/plot_profile.png b/docs/plot_profile.png
new file mode 100644
index 0000000000000000000000000000000000000000..1ec2fe346206fde32b7dd2e66ab362227f16808c
GIT binary patch
literal 14205
zcmb_@1yoh*x9=7tM37QJP(o6ryCkKg#Z4$BjdV9iDbjHeX#@lXr5gm1mJ$%@knWVO
zH`h7u-uu7rjrZ=jV_eTTvW3OkYkg~ebN*_-R#lcIz@xxJp-=?!at|J%P#5Uo^(-zn
z{LN{qy$WA0nn@^0pio7TmrtHy!M|yY<Q^%ZP#%mZ)QdnA>IlB|ViARM=0KsA4Nxed
zw<r{uU2?UW2>bxsP(k(q>Kyqm{Yy?Pd<Vx#UP&5f7K?)H0zunL?s620_oe&;iN|gu
zYvZn7j|WrE&j(X-Z^T9O-2IRqf|7TlCDQnMO$S$d(14Pk>a4nO&xWtvl=|3v)yOxg
z%a!_R@yC^`hD=w=XZ)>`3;!4vT5uWoNd!lz#;dYtCiHfWdHVVqgh$I_UozsE>ft7&
zipa`yT>rwo^Gk4H!ZRXcr8<I0GVns9`Nj((N&F~r63IY<OZ2o4hW>fk58uZn^b>zE
zkf&)>F<Q0LExfd{VsyN>)|o6A%WL`T$Kc?r4CbutY=`5uiV~kwkDFI;TqoR_P){pt
z-hce~8i(2|Dlw6w+-CgY<HuL`Cw;|;iHUy>3~c>=MLq3M^;X24%W=N5!GJ?mUA^%Q
z>r?;TmZm0E4GkF+69&Gd@KKwZ@Si^)KSs|r5{k_Z6&Wce@r6^`RQi>d^OakV%E-#%
zxVyWLxy^*5tS0MHa&zw*8W|;s`HCv%IF_mukTY7)UYhcvIX&$0ou52EowP7F|GMT<
z@_gX#^2&;Kz0Z9M3yZ5)uSUhi1?A>`(y59e)Ax#AU$@a=BV(^W;VP&<6VfY)hc|=6
z-8a7MZ@P(~ehv=K>>SF;$xZFopWlb2Xq?&F{C1iB$K+(sQh$#7T%0Bx_E=LhHY!S7
zLE*C4@o>-L#zxMkPhqe!uPQ4A^t^Y`ckkXcx3FlKbBK<K`Sv~T-f*$W<<ryCpFe+c
zn0As!P>XeM&vjh7h>_{-(aa%#Sss<XB6R7}CDf;<6%n56RXCZw!5PdCWMpDeQo{R}
z=()H^GJE9=4R7b>I?T2avs}O2+uQrL%k(k&eO6XrX(`Xyoc>vmU(G|)o>V1y`AZ=g
zcB@02xw*M^0yo`6s8ZdgFO7_j9`;811P-p4w8t<ohe%pmv!4G^IS=fZ%YbXIa$MNn
zY$WU-AMaWpuhA~Ix`Bm-1xHCqVdf~OOy}3EPq#CrH5qWgVm8CIZ+U!qAT1rfcjtN|
z4J$pdF$ej<uZp^6KTMntA3sW{UA-#(?*029{t34Lv6J6~R-@$v^^=jqR@omu1fL%7
z@AS%2b6St^uw3tu?$Edz6w9RYNL^jZ-kx)zJEifr9T73H1Qb+t{>-k2+Epy2{-II}
z(wWvsCN8cRcz!W4F(^Bs6enlr3+c>w;o;$rA3s*r(Yg7lV`IO>enuu+G2!+zCxYnf
zF}Qs^;mOH2yHkW)tDRR~?6Bti`QvnUwsD@?)bs!!8`GwOQmOx~kn`Pp8lS|`=v41L
z?e85O>`#1*jOgmmb_=vC?P#F{OG`_Ma4-JgAwZ*hoo{9K;wIa^za{GR*7u~<m;Ux`
z40v2U3z%Ow_Kmnwd{<mtT(+TNI4pm$Jbd`jzjbRq$+mBJxShL4gnJ|h9rmI$s)+_W
zz)w2#y8YDG3(L#Ps#mde2Lp83mHPelnds?NA3gf^NxMYu;X}!Oo++>3s&ShWA2F&(
z9Nj?+TiXRDvBOL7bjIvkZ!K0U%gA7_t*t@XI!-rTx=ly-`t|E)P`0)E6P_ZcfBNZ}
zm@uQFqa|R+d3$>&i@3+Eu3FXVLH*CM+Y8(j+5C#xaOwfYKlH_uzen_#1uCnp5qI#)
z459BqBXz`{a&T`>?p-)zW?^B}Yf>icR^z(EF56;@e#X1^?hRKtGRa0#nS2YpbW_AF
zYI@qxrglpvTO}=`b>{b1Kg=<gQR~MV8rv)S=ZTrJk&m7}?FhU?ioDlG?amdKnw0Sx
zS5ji)hS~PmrS<jswTj8?j0}GqDz_I%Rj*B+g{P)c!Q+#4cNbL8QD&f{3$Pw7FM@~s
zN#kCm$Tsfr{>H-lc)i;~LWy5I_cL<uZf$Ly?1-In%Vb7#n|0f5&&o3@B?b;#QLU-7
zyp)JG8@`xqJK*HvLUQd|P_yB?f`UkRhh$}@&!0cLebIB$Z8a>4iiwd_Q6c2y<ZK}H
zpdm`HsYyCGa639Wf;|njje7Lx5whhoGwFsx9E4JuT3SeOFPhufv{>=sJk>@@cx5-k
zg2J>U0d7K(k@Wuk`%X?yQP41;sM+E}%gf6RT-n*!uDoiLP**4Qsi#)z4-O8dVPIew
zx_0%q5{5i3F`qF$@;7}{Aki~A7PjMxva+pi;f;$Mp2U(mVqT)uNCkMTe2yREW5M{C
znHk`?(7p8)lQ3*AAD1?tMq2GdZtnEv@5Vu8nhW#{3~E|h0wXoFftBAnO}Tk_c^gOe
zGW@bugy7PkK|ovQOCr!6Tv}d+ikqI9flU|LDg|p$R8*w(Sp)h^W+rCdz(Ykvd=&i3
zftM&8#7T}X^@}0!U>m~^m75#!<TND1fGfd2Z715RoSi}rZWP)7uwP*bRPOIyHS0+=
z8OTMW#@rWD6cf0DCQsL0VGnb<tZGsDp4ju*PIT6JA9ExqiflEJe#rf9S@1=MhRC2h
zS#ZpGP&d=1IIuFcN{7wE(=#CQR~WWZKL<H}Q)_GF!C67*dfB+m!n;>g-=3Fkr=1_B
zB?~&fL<+^A0V5+bQ@g~J;LDdUD6jM5$ptg9;~O|vZcb%JicTk(r9SIPy}vqCgjqk?
z14Wo1>Mg|n<qygG_wNO6Mi|`Ziehbw`@Q8B{Q7l+O@-r`k+WbLYsxzcUb5q4y&=3l
z=M)$pCqCzK1oY6Mt*lr)_t%@Df)W!Sy1McwBqYG<!QLexAb^cXN=X^ozm%1qe+7<}
zl*EAXbW^|%1s#aTf&|)s16yk2@AT~KmusUH>gwtN<2L7a^_ghvu3weDq)a)p^|A7F
zIrz3b&Z;RLpMrT5lxuZOjT)L)U@oj<?)L55{%vh-62}}r3p?#}igh*h^txf@y3Iz*
z$|)<in3mYDj|swJa3(0yGcr2tt!hIX(ACvdOcPC!dUd6DbX4K#Q%WfMix?N4Jb5Cg
zq9WVR1G@}mZ*Sk3dY>|<zdeE5EPwqTQ<L-h=TqksAI)owjhl3L(w9!sY-bBgdG5`c
z0WdJ>&sLh9pO1khaB<;7hU`pB_{`x#TK(-~uU>Izchh9`2UzXx?OvM=I7k;8^V;e9
zr1ynFTCH-ch4RkbyGY@l?E9YcL-`?4Bs-hY=Xi7yc6X%c5p!<tb_=C#q3^j*NJvO<
zXecHPBjd}rZ^@y|%F4<#x89Q&I=P?RrJ9`bGfQd??(&*Tty*z7{_uRQ(dwoT%^^C4
zt-`h~4Z7;~&Q8@v-7!Z}7tC5{(el#L&DTa3La`sgMm;|&I1icJkX5^CXlVHMz9-*}
z8#l}oWu8C33$0}#%{LXEd`xU?(dgQkW6%B0MBdQ0HaQYf(hF|uRSTOQU$S5T{-o0!
zdLYp6ltQI~TJztB$rxOez@amLu<#vOK(lKtDm8q$*;~Hz^P5cs+L~7N*RaAuxH7zX
z?&^zZ!2U$~Zhd_k%&_aMPaZfq{ke7vv(B6*L6OnH!9lX_{7@B|cmVcC*6Wvthlc~o
zZ*ZG-UK_Qb&Chm;kF57OMJ;&LlEVZYl#h5l#jju*MQ;7&6PNI@S8`SQ?zE<}(6LuM
z`WIR6g<H3DseKF<)6Q&}Lozsg&%EBSXp_{&xct`9)|4C?Bfrx$!1u$N;-GRFf7}|+
zT!{CUH;;=uCv*DJ6v2m!r5)!qK6SeeMQ_^Yij%^dWS?-{cqR7rrYi5M{ts)Oq^{MG
zve3FZq-)$^W;W<b<V7|p0VQ9dL4!a0{jSi9F_rm$c5$>ED~UDF%_K6Hgt)n9-lh`s
zn~%mFuN*XXmynBDzyFZi9qU*kRjq`#P-BiAq=NaNHJzDOvw%l8XHdA(^3@e$O*bp^
z%Ijoih3R&k!3D+5!>>{)Ps>MbV0aOs(XXM&vav<HSBQu943!@5ZL%gwA~#*#H$$-5
z-^FCHx_$BK%?s!m>h-teByuU<qRJK3`C1hp$LD6M_(fMn8V-LA#7?%=Wf(~{UbEyU
zwXApGHUXSUdh6Ex%7~$%d(}P2=`v6<0sg&q4F>r(<PHY=4_%Wm5|y2g&2)tip0~WL
zYumN6m&1?`zDHkbxjzs_;Qj0r^K-G99=W5zwS*2q8-VEfbTAvCbsilb%UV-EMmHO9
z==}PIM<OSwBfIJ%%ouC&_^CYf+p;I@mQhbJ<aJ1SbX(^L?0(0&1YoEU$NM@??|wHd
z(s8=>+m6wXKL64h8+A{;`-?ZJXH#9T<vBHcc7M!GCNqDFgq3H!$dgj?IZdyD>qT}-
z{la3$*uVXN`9(vz0f(0dR!6!K2Xj2{Xe6$FF@6luZ7VsNOSIo!UPNBDk;jmuPcMY~
zoKL;}{>hY1BlhI;iBPpUjYVS1JE>O!J|_bFlV03cUIm0x^7jFy$hJTxS7}9s<K~ov
z6(1!uTi_VQu%WwC?<)Wb6ciK$K%)*s;o(C(n2MRzxBYS1@9ht6wT_zw^nb$UNqPT%
zScH$svD8&5%s=$(MjsDB!C(gf1wr?1`s<-VpFXj)gxx?mhh5D$QhX(5JpdaeVSvTO
z#Uah+#kG4HKgPy7^9t%hdKQd!=DSuW>&3FDw4G;H^2P^xw&o`!&PQ}J_lKk^?T;T%
z2><wMNj<Qc6Zzb2Et3NN`C^pcV<q1^j!n>L7OHuu-d8L!`xy@fGD4_*PdE|0%Ea`F
z)2KyNU%w|`w}zbC6u&j%7J{P?fJ4fptg5EwH@Y@I;kjYb@rD&}5|>%`O%!y<A%Il_
z0|W89W1*FM1S0mfC#6AIE895BlH*ij=T*PS33OAZ2;Nx~@y*W)6+O*AyQ!dG*0C(S
zHt{7XBja|M4BcfiHsjylf~-bLU&A(wdiyrG*BszGGN5<t3B5jBJSj|GTr|hS!z<Kc
zfyIv8t0KC+ia+Md%aIb2;ix<5>)amQxVPs5qre<^)y&q9p&>b4-CNDg&44<A7wqot
zCJVb>N$QeOR=&cdl9~lHVq!v%QSt4I3}&xCy|RGH@e?hD`S{*g^k;KAFF&@R)s=O+
zrE_}0<4eH#E>0=Cer098u7he#ZS74o`eS~6c6K&)$#CcAiiw9X9;@@!b#=R-f`TI>
z7l47>hBLt=`2PKS<?(t=dso+7-SN87Xt}kqiHX{iClrr*b~%dvI#Be~7n^inA${Ue
z+pr#=yCkBP*qq70rNvq|B#qf_zyU1><2s8L$Kr8SRn^zSx--^_v2XA{*W!Zfr1yS$
z0T7<=KR>F=9&0tDp+zdEth_Jk==f@<D&B9g64&tD-rnI<SI3^VCxtd(z|d`Wo%`Pz
z^}mqoe|cH8%KoxzMt#_0^`=c)NJet~U$zzGRjgt+Mx76I;wq2mqnZkx-A_*Hon-?Q
zq%DQAT4Xr5@9Dv=CMJ<H9FmL=`PCz~DL5j`xOwn+LTti+lnl3yD)CBe-_5JpEn$1>
z1WB+F2mRP0)2Kgx#xK|VeL^R4Zg(6QB`!2FBqc4iZ|2&YX;ksH(%G=@LxmGVRAA+m
zz;eBY5BR5*-PzwweQK93n$8@irEGYNGJg+?1))Gk7c+L+yZmBRgaj>1^#qBWf2Pv?
z$9zsW=ZVYJ%Fhx5{Qt5vTb0nErjNRq6v4AInWSuRRC<Y8)!E1H_~Ipc`vKkU7P+(@
zG98nPNb!u1hj+~3DSw)Gso8Q9X_(wPKCRTrpW0Zgj_8>Cwy+bvw-rLm{Qc`equIG}
zE)gdM+lVCIMT}4S^{LP`wpUAfV&dal5yAyMFMwMjRlwdKC=MZDhV!%I_Rh}Pxp@8T
zva-03>NyTQ_xB5dpFpJm0Twwsm;t~Ov1fmo_S+BXr*c5Ol5DD9cXlcQd;$`VVDr?}
z?mtWYfM`GI*1SP<3HP1(qGyfRE^8y89f<67a(kU@wcZqT2-K}{DRNqR1VC!I$Vj%(
zE3HS?%#0}|B_)9IoJ`W|zODn_=+KZBFtH1T01cp|oLpUZzTYb-G;Svuv*CxL18lpu
zHu~}NXA=OPSvfg?3l)HlGseCKc3NCqEa16kRd_No^5j>(F7g9_N}1kl`NFM<9s(mj
zNox}`mV~hGaL#rNMOA(~$(0cw5UY8d6VU#_fG=Xzewdixc~Y0DP1PLW#@}#PJUl!G
zMOuq%YYarV@FF-o{C!15{MnfgTnzxM54j<OEAMl1aB9Y>SlIdb<LggH>*o%(W<YcC
zF6~)zHk%>C*HEB&eQ#}i+e7bOUc<qXrxMzkDkZjEs)S4Cakn(4{GjvR>Tshhwf6%c
zG-_&Uuu>NzTLU`g43GXSEmn6^+SdAITXd!f$NT$ZykXbxX7@epnBtO!W_&WwU;hp+
zsFli-MXPRmY3d`(=hJ~oyFv;9--6fp8leT^>J{;rWho8UqxA<WKO4+`NmXE!030@X
z*zL+e5e`7>S;Z(+m9eogkgTsKhdb@xze}oJ-QM1YAG16y4+C-(9UCj<<ivw;Mzii@
zT7UsXMy(NhRiu*m0j&kaGsY{2!?#c=J0-%Mlxe|TA_oyR0T_${zx}TT6OfI$33`zc
z&`Ybeb$3&hR?fWTG{Pn`;EGO56I=BB`0*o*FYNk=kdPP?*jg=4Ri`JudG`=Bk&`4J
zUbYp#@UTd*E>Tv|*q@KKDe7vs>oHFC8Rq$U-s&HP{}`nIGE@KCFIyNe?-1dFN=2Q#
zGqSNiv*qTmpO+glx|Ur~K!`>MOzwDmX+X;Zzg^h{wkmV%7DV)HH6||tk(^kcmT&&)
zSBXtY=|V^s)H6^M7VY9d5C%(JH}qN}Zgue2A1VRdb>AJ(w0QnJ1B4vVFSgE(H%crA
z*`}wb-BycQp;1$>&D5E?pWZb=ueo{O#f<&=Kr5`X-Jawb)woGi&_OX<)v}z0#U1zq
z+s>e5_;1-rUFCjrZ>H?Wf`S~_!NI|}rrpWypa&i<W=1Z8BDvJ}{;#+=HEr#7*f;J6
zzh4z;v3S83sb6$2b2`=*&g(QAxP9-M4X1r~pI4bXncTm-=x{9SIT(1*hrA&S6Oc6)
z_|KenqZ(fkg4X#uG7_plFj>es1aAJvz<~5Wbj{ez4B3rJN(8VQK`8;0{Q>yeo5`O)
zo6b&;VMeuV#1sBSr0?2jlZdZ)#rNdE{?{k%#ux?#Xf=u{LUG7wN7@Uh9foC3Ku^6K
z!Q0%_l%D(C4ofGg>pDxn72k}VKyXCFl{W>;LOy89KPm%%dPP+vgY%}~d904t#Csp_
z^^S~$0HyUkTkwUNn}_WPGK0gx<}(Tk3IXpU9>7eq3ky5bq3qjv`llg4J+fhJ{ce1v
zUFRhT-ALqQv$0{Dl!OGKRPstnCIE0ywVD%Ola4Jpy@j;h#QWVpxUmp+{?&XTuyR7z
zVSpG29z2MkUUPa*&gg|y&nqC=wHqV6vbI(PzP#nN)UX?f&PnwMfic_?Ngdhvd#b3F
zY6dP&AVnD=B+KNL6gaAbOQt<be2QsoIV7sQb)0O^N_%37rtenH=F)+fHaj<mRUaQ0
z$7`Mh0ECMRn;knoAwk~M)O2$NSez-4gHN}(nDH0`Xon=Hx#HD9@bH}7e)FF}=YM;H
z{DW2h{<Ri{1t~l;sV;g8p0>vfk~(<+;f&GfYk+7FT>>;*L@<PcY;S9$frbdJQuJ_^
z^s4l;kx~mSmg~xCqDtQ0_d$UHWd#bWA(~EBPEpYqL@-U`N9Pfu)M6)-4@Zu@Fhgyp
z!b&$AZLg6~TijUN5i+d4k|-)*KoH`uh{aH+%tj%XoBIfjZeA&<?~YXwl~+?EX0JQA
zt(+5!Y)Ke;2q;!g7Vyi><>BY&cHcJE^Eou>>+2gTF(U-XI1{I-*I>@m2RL@RIW$4Y
zc{zS!>+FS$Hd>fKmkK|`KYJysd1R5$Jkwj5y(5i~Z_~zIw=gK@xzIjY_}s(|i_9;U
zk~+L(*12yL8RvV_dctmS|8PzLP2Wj(Vc0CqWUSJDb-3geoovMBQ~+VI+m<1+5-?cX
zdU_T>J!x3jv0wZ_J2Eojw%yKjhmCC>M8Ae@oa=u%)r<#Z7u~uNv@_1#k98h~9b{Q1
zp*B9zT0X;SCD!|WkMy&Z@PwqZ=cM_@8zRkt!q%v((hW0Pp#Aza7<@G?!P3kN4-OXZ
z=K*;Mq|2t}W~_xVmnt3da2e3Bz&?Pcj@9uAWG>(Hecz9Hd2|d6CZ#>iAGny?etwy@
z+=vMf1kE(QLt0al3^WYrH+IWE?;sRIPmc=N$oKCLuiO#}1brQzF_;U=sUq@lD3A{r
z(Cx8Ig3ilVp|>{_J32TZoL;Ni=}vTXG>qavqRiC{R(jB-%dIH^r#E2HEOaHc!c$b%
z(6|m%gOPzD2qr(7wA#U5MZFTR7&xl-a6W0wsb2~9tt0p;I{^EmVO-U^@A%JAhRXnG
z=SyO(WJ@BNURjB`H#p&QyvCwch!UWtQO{O<c6M?IcMh60Y@w4zu@kAy1(b@&*2Uah
zL#_nXt3d!*IH9FYZQ#;2sj>;aH0BDjG2~0a@xD*=023cuN^aBZ%r!qd(S@IP&AUP4
zZT^a}5M8^|h^%#5Akn-0B}%$iE9_}J2o2`&tv<-oA2>MkzvQ*>tCn4R`~84R<wqOJ
zn=iF5MH@FMkmGH+>xYW>Ip0#kL<{;_`T=bu-g<tp(}3+$en`(alcNxjCp7t#<19Du
z9noN7fiz7xqwb&@?|rw{vfL$2R{WsKUvf0M>E(?|N;RH8W+*gb{XQ-gX(<M@HaXpL
z5>>d4<=!b3%&vg3gyAbV^i_p@QO}}U`T3`CbkD?Ba@hxVHQl1r=eFaVLPtj+R_YGE
zjx62oxxDkb20Ja2^tKM1T>MyP|M}EP&ki<^OX_n1_lR!Fw;Us{XqtGHFqXJ<4sMPg
ziccr{L?{Q_(YoHj`!_cG|7k2O!5dP5NxvT6k-<(i&_rL8R`c(VGEcldhK)^SGa#9L
zhgdTBayWAoNX2eDUHlf-*6(4S1zjd9ww)Ap(lvwKglI*N(de!C`V*|s=5MfTLBK)0
zsg@QQL^_{`)-%$|5{Ro#w_H4I#dn=uuTzTf{sS5!*b-{MhoIRefEb0!)A*cGSqVU-
zr4WocmyHP}7~7KornD<;Zo!|1Z62k@js5<~NwI0PY#`BK8ZlN(g8>Z)<$x@~Zn?w4
zG7CCd7F;1vVQ}aIM%PTYXWOrmk@Z7A0j=*X$1{vrMx}37eEkCh*H|{-d$r(&)t+p3
z{Ml`f5tH;94|h|2+Rr0;v`7cmQ#i~6aQ^0(a?)DNc|hwA&0sz_^?}s@QFq9S5A;Sq
zaT<Z1y%#W9Q+<vssJ#EEX*)l_<=mE7tF9ITk17Fx71$YoivytRf>8%7IjgQNwPL~p
z(T@NZ$b;cwTX)E4ZEY=Z^9tx~VtHf6x9{E!`S5}HpL2<WENBRj0e%b!S2>3?vb7Ob
zT^0}{Fa~ynn``970Zj(MFqDiv8;ndKsLx7<WscthD0<%HH|r#1FO-5BTHBLY`P6Xu
zwR$N#3ESHi7Dq8v#Jw~UjY~W~+amX0(7u0RfdAFY@D4oaCjR$AQM5bHR;op3_XCNh
z4P|(!V+L2|5MonPBP?*^skZiwxw$#$A0_6!w{!FIw7+=p!mS){)ThDZ)cE{F8paa;
zxDzR&nRs_f26OOcT#xBlm|5)2G&VJLXe8kO(bxA4KonwZ!KrQvIR&HhG}Kg8#Q__y
z{Q7JJYB&!cpF;<0MWM&;b0J}2sqb-1mO(@u-^7anbj?ssrk9k&4(j`+6c=-W5%QMH
z1P_!)1n)zqh4}@B?>ksgFlzuEBl85FeSpR-Fs=}f6p+pDokN*p9b8(O;oxd3@QXzb
zrUFnW%bNAm2*_z|#awu*0u18o;ehphDnjgl3!%;aJquE3^kej`QjgETRf>wD=jTmt
z7?DVZdCiRwik%QC^?xWRSO9-y_H=)8p=?kuv%LIG&3g6a_>Sh9jXDKeTXyK7wBQ~g
zCM_5~85tSi@#4OTfAglrdw;wYSmO{hV)L}iv`x&EYBUL5rXaJ8kB^J&4H?0L#{xlw
zDenfd27rPofW2@x_;~~0LHPms7^<9tlCo+hLKrw7?cBo^Iou0$`1GO8B*htWHa2WP
z`Wgz13=Jz!_eOE7)gM28pr+PqPMt?|yR$t!S~9qh){j0AmlEW?|B08J8HZ6kN}2E8
z`GJB6&j|5Kot#uPHF3mEO&Njmm6nwuW<XzmKQte(5ulV_hlJc_WmRI_C6mOrXi^`t
z8LK1)YHMkEuSe`$B#F;Dqqa5~_6#`iK$%aD2J6A?GU7`@tYY7@oiqd)uk3FBkRs%?
z=*#LAwV{1LPN5dh3LCSP+BX?MBObVy;E@?~y;~8=%_X^UqrD*j2e#SC6pk1u1c;uG
z7~{nz9f;WOxwm@h4L^tTT2QDm9&KHn%umTjSFu1G%2rAOGlmvM_Ux>&Y3ExB*c7()
zr+o7A^046n8Ng)!^-JULX&Met4qR~mIW!twWHT=G+|p8BQSs)QBW)nL3EgwrrD`ge
z<J!d+le%7x+c0r+$A*RBS65fVf&tv>|Mg1-?q6G58)Sb+N5nR3Fo{LMs<|(|dsPVt
zxp3I(A2eX~fE8>VHjR88dEgpgDX@iJXPs_(P5M`SJSpx)E!b_)$Y$H3=}xxe>*?9q
zBOy!B51w$ywR_<FAt?It8?&+}aCKQFC5}HltR;gBa-d;DM@7sMKXH^m{Wm!W2QKh5
zo1T{$R*b4$4a!OLp{N*D0<QwB#>1Tj^+%6LbY$dlX<u?>h;xzSgM%e8$qnYlg9rFW
zsRPfOpO@+VK|0IM|7i32CsIP*>m?W8G3QX^vi1bphIXOB1*m$>LZ>;$^DjeujjWcH
z_50{w_Y+&a8J4MRsvhkx4tn!WE|9sRpU6Eew?di*F!F6M%cplea8A7pcBB2y3zjJo
z^D|e{mq*gma4Mu^`<-}{qy{DWdwf50H<dTC?l_KGB=Q+2XlAbWimbCg-47a-AOEet
zKj4+7ACjTpZ5#eW#rH^kzB57U*)!VnUH$Xs$<y)4=)eAY0Dzy8&uVI{%JJu~UtC5<
zPFH?9ITYr{zL<<`ouR!e^zz1yBE_CCss8Sb+ze*KoY3qItPv9jSCp8K(q4^kJuz7_
zYv8KL4eWq#4T;KEemU!wcD^eQMXW;JqD_CU0*&!Z_c8tK^@h=p2dAPRg*#Wy?oKXn
zbIn*-Tha*%CV-11DuS`JAa*VY&(#>_%gI5j*nG|8DI3`MV9Wxp_FpN!6XGv%A-1;-
zYG$~|kx3!V4Flc(<K>V-tOVH6XO)InJM-NPu44ix7kG_pe!mg?=~A6fN}jQur?ziS
zwIBa^)L0d~(ntHekM=bWMF*WzU|IlWuly^I>m`r<^aD`2q(Nr_klb)sX3a+l!pN;{
z*CP?Vn3&e6m^VJ{+c@P4@=7&}E{&|}%WtaRuNRcBvl~LP12O^~M_t6^!Y(Y3Z~3#o
zk4z@$5$uP;(Aq&3xngPm0<spMy+fPS%u@XL6n?*cq`JNhy;OXOw#lx|8)<s>=$0s%
z4Cq*p2LP}k4NfZJ3cP<$u;L0+KqK#3z%#q<<?4uBU@8lQ{Gep<P^kUYu&)gARc<ah
zKDLL42XJbt;@8VtZlL)yhkS!Dik!SWEfdp1Mu#7vlXhlChUA`Yx=lcXI3$Z~Dw49Y
zvLq}=3Bao(4{uT6OTj5$+t>jA|Hs6HGQ>#$>(h=73=Wo9kMccv`V`R^p*UNXoP4ga
z7aZU@ySO+&d_zb`2xJd4#BHFEJIm+^$hLq1YRE$YRvf-)8ex0~@2qF$=g;BY-O3;u
ze*E-_UQ8?vl0d}Njp!Rv6#3AjK@d}hT@0cNq7^paP;*$1ltwi(JL|SC?7+Oly@&xz
zyR%-ifp{2?AOB2`tp?DCv>RVvF=$X$!^K#y8aWpmZ{S{Z18$*R<HE|MoQxQ4(CcBm
z6OqsWV4QZXKfcu(Nxd-cI{6-Q7J-2mL2i_cSO9zgBHf11zA&HlNMP&Cvy+3Z0KF5q
zWB`^R+)l%T#)POv4vgM!-@ZjLv}oq*wCAGpAh`k|252yCl$W}3rl6g{80CI-<<@in
zp=f$;Zpa<=?27^J2Gr+J*oGYBA3l7*J#z&bE%j0ldMk+cMxZjn1|*@T7Bg=1J}Uzb
zdlgDYi{*`km^wNR$W$LV)8N}cRBLZ<Pom)W3xLW9i-TJic3p>V_%ciuG--W(YB&U1
zKO}fytiTDPc4oKSpztANP-NCa1r-7;9pC-SsgD={WK=B($A87e5#P8W@8xxpz6M|s
zhmZ>Ct|h}Sy>ziFb=X{1hx{POG7~9=#8UqJ`9ph-TUJ{7H{xC&lZX5l>{Hq_Ael%g
zYH4v1K)O9R>wvIY!9tMgJ`CbfGJ>ZdkkTDBXb-pgNvp_zd)wivv_wBox$R^(a@6u!
zm#t6Z#23%uGD|YR99}@BJKu87WH=Kte0l;wipyv;0@x4|`%kuL%AraUyB*r>VY`YC
zMza15O*9%Ao?ZNtBIY%3-dw*zCDN859g6r!kk<fzCpImOx>y%<E%0xw#;fDNa6<AV
zKSo9%qx3!98cdIe&=LD{R9a0-I6;7hmPzU$0!q^AXhkFl+E6H3wQiigNBwEccka*%
z3cCNW*VNE>1l|_BYbWHD+<pg<nOj>ESln`Pbw#v9nD4_#a0~4Jy-{1jswcznZqPm}
zCVhl02lMA}#11YXK|Fx+SvkHgSXmVnz?hN10>o<;B81ls)<(*peuMrR9v$saVE=+}
z^*eW7!}HhkSYkMvw>^8gxVSjwBAONnavQeo3p6?oB+H$Xot{^<OC8_8YeDb>W@BI9
zqvDd1;Tl)=NU^hn^fjmix5K%32(kUy-HlF2Xb0U)YI%KibyimF2#Sxma@+>OO~?!b
zk%f+t(Nr~olHUdqc)*aLf}qUeVmc&;K!}OjZ)pjqgt<Kn(;n2r5{J2msSx3l@Ee+#
z=w|oXx(s;X2EWZs{xLg*x=X$61h}puEsX^a8v0|#iqJ;=xhRMQu5E(8Ag+Q3aTB6t
zP{*K8Hi96>%lsEW5@=E<Fuiy|Zu+ygw_TEU8kV1*FY$X~V4xXUnQV)O9%R9RK*_X;
zQh`wcxHoBROh?c6)EUuiR$V;XT%hcS)?C)V?CbGa4SfZub+BgZs~{`e%5bjbt*3<r
zV37oQTu8(Why_qE2v0qL@3xggN-3$U1)-wG#`KUGMJ(Off~0H<I3eEU%QwV)L<(xR
z@nA>Jf-Of`nv%Z)k`2&Pe27mV^8wB$PDAv_)yGE^3=PDX1kD9f4Lq|+-R3+<t`YTy
zP4^+tJ7~NIUUd?q4D|HiS1_})=K%HuJ2i9VD4p=hCKh)><?ok3=AU7JVz&yGCpe~2
zn}ECl$|C&=GLE2j!lMmxx<^V%`Z^+F9(k_tt`NfEwx4dIeaZ~051^gOYf}QkPw<bE
zlLawC7h4=Cx2>JD^W<;eBq3%A_X%rHtG@&2%??6xTs|jmdY)^%pd0bbUP}bXhop;O
zTR?Tp0_~lzs6VwXOtSlJuad}f33w+I1@s~o-D;<{8EEeym%YoqJ9#{@FW|OGbM4wS
zT0dj11gPP!a}I!;>4CWlxU7=fRCub=e9Zk0(pwcYDWq`kY3758_0lZOn{T$v3dmGx
z>0kFf_c*+doZvcPUQ=!-TPC&u&utB5uWR$ix~v5=077j6dOLo324bg6DuDH&6G9?Y
z($H`^eOrWzu&b*J(IemqV)%(~M)oZ*FfxKU0CA!D*G#EXtzu_eQH@L>wPj_|gJR$g
zc^!bA(8Gx<B4|3Rtsn>kb`A*758<K0l4n*{y0zs?J$P^tpdT1ppmT#N9!4!jmD4Zy
zA<3vn3(RKFOrhnr!^$ANQ%uZPH5futBqSuru;%6_1*opY0+ud;c72zr5o|fdM4s1S
zVPUitMNd8h#0P`<D=<Yk=xx`06jzriH@q1D0FcH&H?2R)tH*QaCb;Peq8o(TUNsVb
zI^VHn2Adl4vrt)`N&Jy7Ut-U9C6R64eAjMl*q0@bsJ@VZ1u;Piylj`6`&pe}LLdYi
z))2;?8N_fP(ggc}YY|rYCHxkIXc0FnCWZvcyY2gTyu-`F95S2$rWloz+d#yDOU{Dj
z_~8OI1W6&FEiG`%aMj7+T*3^M@C{O8Bn4I~G3YFMFog*@h|AYm$Q{P=VQStIc7X`l
zcWA2S&!68W!VN4<U2_oz$bz_!p!6BUm&3VLs~}?sOc*pH1$Y2J`T;Nz!yx^Tguh^m
zan_-K24LYTsDi&iTt*Z`+8H)zhhTw}fVT;H_C&2aOgZq$uuttF7Y6DE^p6%u3$(#o
zzZG`1X``lr*dHtk=s>sc-%o+KHnclg<eoj<K&1iQA0lFqRRL-UjsP`DKCIrsnTx<p
z<`9yT@ps+F&21*WoVQ^S;6h<wAw>v&YkWLe9}hKomL{Uzef;o&Mp#&Qwo8NqTY%^>
zgt&Nkd0_^vN{lX5D}(X{oRJwXHy2USpw&oK8)AYBvU%tQZJLXlTk46tPkkEnDj<K*
zWR8wJ0Ifj^!$VHHIv)X(!2Y3z!LtC10z2+w?p>dHv7tD~?g7q)iXeVek-q}V0Q(34
z1RUn1OL$4@1;|Z_iHWcm7=gnE;1I^d$MY}hJWp(a+!C^TVY!ju8bGfgfO>&{I6%G<
z^xR{C{z(T(1-Oi4WivT9ULKy!Ur%kH6>0h34}mwdU-?{dWyB_Z9Z&#_-pbuyPlv0W
z??G;hd^>WvIS?JM6j%W_11JV8O^9Fv^8{gL#uK_mx;Vv9_nHglg4hRyy1=*e<Y;(7
zt~vT+F?FIf1@T;e-L%VxUUgAi{5E540Qbfv{G6VRV;a$89v>X4LXu+&RJ_G%(GTlA
zIDHF|!bG?LW{#q`F#S*viX`8z6w`NAQc^0k$?R=|)4jAx+1=lF1c@KSJ0QRpG5lzI
zOaS*F2};Oh7(&_=nu?mHW-C-TPTYb+mD^S`S%Bs*D)(7pB+`kvv7MbaOV{AT0^%MX
zLI7fFkCyY0EeFD2lsLr)Psph>*pGk!2wR)RkG@U_6u{a)eG2e#00Izya&w_6>+ZRr
zkb4F{8w{iic3?M4V4jiPD=8~;{<&P&*Q$#~@4jEABEOdcVM<a)#e;7}eSlx#p-!Hk
z?F-sXp$fk2Q2?IlOcN8!QbeP{BhG*(3K9XpZ&BmRHkHnVuP(-rmxGXqOh}uWH9pvs
zkU|0Yygle!IeGy=Spa-Kf!6I#6=;&n%g;Z02AT+^C_hwC=;aBA-$QO1NwV9728V^w
z3kgwxgz|SDztHP&o3<4`FrW}X=8WC)=fuai&86|PrvSu!!`c>XGbrEYPCHZxiw#DI
znya34yg77+AA(!E=3-Pb3>X$NH-MZ#7o**;4abIjg)~HAU(i0q65qr$s2F|aTv2u4
z0qN$ZYXg0KP%}GZXF6!~+2$)TuixJC;+qju0G>mC&Fvo})dmAE(|c$dFPf9yULHX2
zEHqS@mz<9Cih%Qm;0SPB5U;{Ku?27`5*dZAO--!=F6qi(0bCU~L<FH94prLI!H`K7
z^`=1I5@{cVE;<CCb+9{AKA7u>=Z2A1ce-W^N=gfGdvd@C!0J}Vs$!uTfI@>j5%>n8
z0Mcr{ZYq;~A~zj_U}unw6d^o{=sNZ1ryKu>OT~uGc<{jrgkT`5aCLQJNC+Mhk@&~o
zsfC?{B)u?ku0*AzD8W+8-5Ra1&4G5hHPcE;MkWgxp%T-s+k%2Opr?6xi6Ey(LdnqG
z0U_hq;;{m-rKNpVsuEj?0ihfQZI<hxFg1kpPXy=Yf-y^0F&dtmJ3BjzLgp@#EAKMB
zEXtq<C$23yKt8~dED&uAA}I`umx?GcSWe)=NKRzaVi&L&+%hnp%>4WasLaHqBoN&U
zp#{SUiQw{FT=s&~5=ef_L-42XD}gV93t_?zmzc$ZTmk#WPUD4No+u`vb#_ipiBann
z<TElrXo7-*ZVJ1E5sDtX0Qmk&hvCuHAZW)xXCRr+0QM~yU7&nH-sC-`AWMc7X^6zt
zuEMo2#$JkTZvvkaTx!a-e){|ug&uE|0(0Pa00&@zV_m*(2ki>m46!%qSy{t@FhY3)
zL<JzFuA`#}y@@LFm96*6Q}~<>5<c+afq0wUk4!u$3PfCn06f@=DzHE~Dru%j+$B33
zQ1iDe`B-LVW+aRWb%D(kAqJK<;Inw=$xvK&hhcr4J|LH{hh~8VR2rI4!To_UkZC~y
z*}f#(I&JLtK@eOl^&947=cAjKb9Qoqcp7=FAjkdtr<bW?uIW?%qooFnfVi-Aj``C3
zo0&al9ekt<WniKuZz?Z;7ljM|z(kQ;z(ieuSIGaI7cfvn@RbYxw|I5&@9#;%Yvhl*
zxvi6{iIo)%2OBpBn*a&}g^_;g?<3?e;5`5Pbt`fj_;b7CCvqYNCuuDw<7ZALLPicI
z@CC)i!NtMK!OP0U_n4Dgh=WszllvYghY$xxb!w~LfBFS$TVo4T_y6}VbTs--!!M8^
z?_lC)VJajIf2o+7I+{44`1tr;&7Q$e$^ZSQsupf0R$31%tW9hjEo^O2yaGHCk=g10
w_c7Z4aSTl1A1g*eaLlcLKaaYtsgvt72b2GEJ)(2i|6VD1Y2^p`k_P_&3#V_;mjD0&

literal 0
HcmV?d00001


From 0b109c674f9cac509a2752c89b14e18e616df6b3 Mon Sep 17 00:00:00 2001
From: Alex Guinman <alex@guinman.id.au>
Date: Tue, 23 Jun 2020 17:28:09 +1000
Subject: [PATCH 3/4] Linting changes

---
 print_examples.py                             | 81 -------------------
 setup.py                                      |  2 +-
 tests/__init__.py                             |  0
 tests/test_actual_interval.py                 |  7 --
 tests/test_cli.py                             |  4 -
 tests/test_consumption_data.py                |  4 -
 tests/test_file_outputs.py                    |  6 --
 tests/test_header_datestamp.py                |  6 --
 ...t_powercor.py => test_incomplete_files.py} | 11 +--
 tests/test_multiple_quality.py                |  7 --
 tests/test_open_examples.py                   |  7 --
 tests/test_optional_columns.py                |  7 --
 12 files changed, 3 insertions(+), 139 deletions(-)
 delete mode 100644 print_examples.py
 create mode 100644 tests/__init__.py
 rename tests/{test_powercor.py => test_incomplete_files.py} (87%)

diff --git a/print_examples.py b/print_examples.py
deleted file mode 100644
index 28ff00b..0000000
--- a/print_examples.py
+++ /dev/null
@@ -1,81 +0,0 @@
-import nemreader as nr
-
-
-def output_file():
-
-    from nemreader import output_as_csv
-
-    file_name = "examples/unzipped/Example_NEM12_multiple_meters.csv"
-    output_file = output_as_csv(file_name)
-
-
-def pandas_df():
-
-    from nemreader import output_data_frames
-
-    file_name = "examples/unzipped/Example_NEM12_multiple_meters.csv"
-    dfs = output_data_frames(file_name)
-    for nmi, df in dfs:
-        print(nmi, df.dtypes)
-        print(df.describe())
-        print(df.head())
-        print(df.tail())
-
-
-def print_meter_record(file_path, rows=5):
-    """ Output readings for specified number of rows to console """
-    m = nr.read_nem_file(file_path)
-    print("Header:", m.header)
-    print("Transactions:", m.transactions)
-    for nmi in m.readings:
-        for channel in m.readings[nmi]:
-            print(nmi, "Channel", channel)
-            for reading in m.readings[nmi][channel][-rows:]:
-                print("", reading)
-
-
-def print_examples():
-
-    print("Example NEM12 - Actual Interval:")
-    print("-" * 10)
-    print_meter_record("examples/unzipped/Example_NEM12_actual_interval.csv", 5)
-
-    print("\nExample NEM12 - Substituted Interval:")
-    print("-" * 10)
-    print_meter_record("examples/unzipped/Example_NEM12_substituted_interval.csv", 5)
-
-    print("\nExample NEM12 - Multiple Quality Methods:")
-    print("-" * 10)
-    print_meter_record("examples/unzipped/Example_NEM12_multiple_quality.csv", 5)
-
-    print("\nExample NEM12 - Multiple Meters:")
-    print("-" * 10)
-    print_meter_record("examples/unzipped/Example_NEM12_multiple_meters.csv", 5)
-
-    print("\nExample NEM13 - Actual Read:")
-    print("-" * 10)
-    print_meter_record("examples/unzipped/Example_NEM13_consumption_data.csv", 5)
-
-    print("\nExample NEM13 - Forward Estimates:")
-    print("-" * 10)
-    print_meter_record("examples/unzipped/Example_NEM13_forward_estimate.csv", 5)
-
-    print("\nReal NEM13 Example:")
-    print("-" * 10)
-    print_meter_record("examples/NEM13#DATA_14121801#WBAYM#3044076134.V01", 5)
-
-    print("\nReal NEM12 Example:")
-    print("-" * 10)
-    print_meter_record("examples/NEM12#DATA_16081001#WBAYM#3044076134.V01", 5)
-
-    print("\nZipped NEM12 Example:")
-    print("-" * 10)
-    print_meter_record("examples/NEM12#NEM1201005Scenario1#GLOBALM#NEMMCO.ZIP", 5)
-
-
-if __name__ == "__main__":
-    # output_file()
-    pandas_df()
-
-    # print_examples()
-
diff --git a/setup.py b/setup.py
index b58109a..8d00d69 100644
--- a/setup.py
+++ b/setup.py
@@ -35,7 +35,7 @@
 setup(
     name="nemreader",
     version=about["__version__"],
-    description="Parse NEM12 (interval metering data) and NEM13 (accumulated metering data) data files ",
+    description="Parse NEM12 and NEM13 metering data files",
     packages=find_packages(exclude=["contrib", "docs", "tests*"]),
     include_package_data=True,
     long_description=long_description,
diff --git a/tests/__init__.py b/tests/__init__.py
new file mode 100644
index 0000000..e69de29
diff --git a/tests/test_actual_interval.py b/tests/test_actual_interval.py
index ff4ae1b..b12f89a 100644
--- a/tests/test_actual_interval.py
+++ b/tests/test_actual_interval.py
@@ -1,11 +1,4 @@
-""" Test Suite
-"""
-
 import pytest
-import os
-import sys
-
-sys.path.insert(0, os.path.abspath(os.path.join(os.path.dirname(__file__), "..")))
 import nemreader as nr
 
 
diff --git a/tests/test_cli.py b/tests/test_cli.py
index df1540b..fa8c812 100644
--- a/tests/test_cli.py
+++ b/tests/test_cli.py
@@ -1,9 +1,5 @@
-import os
-import sys
 import pytest
 from click.testing import CliRunner
-
-sys.path.insert(0, os.path.abspath(os.path.join(os.path.dirname(__file__), "..")))
 from nemreader.cli import cli
 
 
diff --git a/tests/test_consumption_data.py b/tests/test_consumption_data.py
index 969693f..836eb8c 100644
--- a/tests/test_consumption_data.py
+++ b/tests/test_consumption_data.py
@@ -2,10 +2,6 @@
 """
 
 import pytest
-import os
-import sys
-
-sys.path.insert(0, os.path.abspath(os.path.join(os.path.dirname(__file__), "..")))
 import nemreader as nr
 
 
diff --git a/tests/test_file_outputs.py b/tests/test_file_outputs.py
index b4073a8..82476b5 100644
--- a/tests/test_file_outputs.py
+++ b/tests/test_file_outputs.py
@@ -1,9 +1,3 @@
-""" Test Suite """
-
-import os
-import sys
-
-sys.path.insert(0, os.path.abspath(os.path.join(os.path.dirname(__file__), "..")))
 from nemreader import output_as_csv
 from nemreader import output_as_daily_csv
 from nemreader import output_as_data_frames
diff --git a/tests/test_header_datestamp.py b/tests/test_header_datestamp.py
index a74496f..0b41119 100644
--- a/tests/test_header_datestamp.py
+++ b/tests/test_header_datestamp.py
@@ -1,11 +1,5 @@
-""" Test Suite
-"""
-
 import os
-import sys
 import datetime
-
-sys.path.insert(0, os.path.abspath(os.path.join(os.path.dirname(__file__), "..")))
 import nemreader as nr
 
 
diff --git a/tests/test_powercor.py b/tests/test_incomplete_files.py
similarity index 87%
rename from tests/test_powercor.py
rename to tests/test_incomplete_files.py
index a341571..a35647e 100644
--- a/tests/test_powercor.py
+++ b/tests/test_incomplete_files.py
@@ -1,11 +1,4 @@
-""" Test Suite
-"""
-
 import pytest
-import os
-import sys
-
-sys.path.insert(0, os.path.abspath(os.path.join(os.path.dirname(__file__), "..")))
 import nemreader as nr
 
 
@@ -44,12 +37,12 @@ def test_correct_records():
 
 
 def test_zipped_load():
-    meter_data = nr.read_nem_file(
+    nr.read_nem_file(
         "examples/nem12/Example_NEM12_powercor.csv.zip", ignore_missing_header=True
     )
 
 
 def test_missing_fields():
-    meter_data = nr.read_nem_file(
+    nr.read_nem_file(
         "examples/unzipped/Example_NEM12_powercor.csv", ignore_missing_header=True
     )
diff --git a/tests/test_multiple_quality.py b/tests/test_multiple_quality.py
index 889b2c3..ee4b4b9 100644
--- a/tests/test_multiple_quality.py
+++ b/tests/test_multiple_quality.py
@@ -1,11 +1,4 @@
-""" Test Suite
-"""
-
 import pytest
-import os
-import sys
-
-sys.path.insert(0, os.path.abspath(os.path.join(os.path.dirname(__file__), "..")))
 import nemreader as nr
 
 
diff --git a/tests/test_open_examples.py b/tests/test_open_examples.py
index 3c32a19..622fa21 100644
--- a/tests/test_open_examples.py
+++ b/tests/test_open_examples.py
@@ -1,11 +1,4 @@
-""" Test Suite
-"""
-
-import pytest
 import os
-import sys
-
-sys.path.insert(0, os.path.abspath(os.path.join(os.path.dirname(__file__), "..")))
 import nemreader as nr
 
 
diff --git a/tests/test_optional_columns.py b/tests/test_optional_columns.py
index bde6c23..da810b0 100644
--- a/tests/test_optional_columns.py
+++ b/tests/test_optional_columns.py
@@ -1,10 +1,3 @@
-""" Test Suite
-"""
-
-import os
-import sys
-
-sys.path.insert(0, os.path.abspath(os.path.join(os.path.dirname(__file__), "..")))
 import nemreader as nr
 
 

From b76747b83eb0242ccabbd46e6a0ea6189891f9c8 Mon Sep 17 00:00:00 2001
From: Alex Guinman <alex@guinman.id.au>
Date: Tue, 23 Jun 2020 17:50:31 +1000
Subject: [PATCH 4/4] Warn instead of fail if interval row empty

---
 ...NEM12_incomplete_interval_daily_totals.csv |   4 ++++
 .../Example_NEM12_incomplete_interval.csv     |   8 ++++++++
 .../Example_NEM12_powercor.csv                |   0
 .../Example_NEM12_powercor.csv.zip            | Bin
 .../Example_NEM12_powercor_missing_fields.csv |   0
 nemreader/cli.py                              |  17 ++++++++++++++---
 nemreader/nem_reader.py                       |  12 ++++++++----
 nemreader/version.py                          |   2 +-
 tests/test_incomplete_files.py                |  18 +++++++++++++-----
 9 files changed, 48 insertions(+), 13 deletions(-)
 create mode 100644 Example_NEM12_incomplete_interval_daily_totals.csv
 create mode 100644 examples/invalid/Example_NEM12_incomplete_interval.csv
 rename examples/{unzipped => invalid}/Example_NEM12_powercor.csv (100%)
 rename examples/{nem12 => invalid}/Example_NEM12_powercor.csv.zip (100%)
 rename examples/{unzipped => invalid}/Example_NEM12_powercor_missing_fields.csv (100%)

diff --git a/Example_NEM12_incomplete_interval_daily_totals.csv b/Example_NEM12_incomplete_interval_daily_totals.csv
new file mode 100644
index 0000000..b67e0cc
--- /dev/null
+++ b/Example_NEM12_incomplete_interval_daily_totals.csv
@@ -0,0 +1,4 @@
+nmi,meter_sn,day,channel,day_total,uom,quality_method
+VABD000163,METSER123,20040202,E1,53.32799999999995,kWh,A
+VABD000163,METSER123,20040201,Q1,106.6559999999999,kVArh,A
+VABD000163,METSER123,20040202,Q1,106.6559999999999,kVArh,A
diff --git a/examples/invalid/Example_NEM12_incomplete_interval.csv b/examples/invalid/Example_NEM12_incomplete_interval.csv
new file mode 100644
index 0000000..3f0f6a2
--- /dev/null
+++ b/examples/invalid/Example_NEM12_incomplete_interval.csv
@@ -0,0 +1,8 @@
+100,NEM12,200405011135,MDA1,Ret1
+200,VABD000163,E1Q1,1,E1,N1,METSER123,kWh,30,
+300,20040201,V,,,20040202120025,
+300,20040202,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,1.111,A,,,20040202120025,20040202142516
+200,VABD000163,E1Q1,2,Q1,,METSER123,kVArh,30,
+300,20040201,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,A,,,20040202120025,
+300,20040202,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,2.222,A,,,20040202120025,
+900
diff --git a/examples/unzipped/Example_NEM12_powercor.csv b/examples/invalid/Example_NEM12_powercor.csv
similarity index 100%
rename from examples/unzipped/Example_NEM12_powercor.csv
rename to examples/invalid/Example_NEM12_powercor.csv
diff --git a/examples/nem12/Example_NEM12_powercor.csv.zip b/examples/invalid/Example_NEM12_powercor.csv.zip
similarity index 100%
rename from examples/nem12/Example_NEM12_powercor.csv.zip
rename to examples/invalid/Example_NEM12_powercor.csv.zip
diff --git a/examples/unzipped/Example_NEM12_powercor_missing_fields.csv b/examples/invalid/Example_NEM12_powercor_missing_fields.csv
similarity index 100%
rename from examples/unzipped/Example_NEM12_powercor_missing_fields.csv
rename to examples/invalid/Example_NEM12_powercor_missing_fields.csv
diff --git a/nemreader/cli.py b/nemreader/cli.py
index 8a2b9f0..1ec905e 100644
--- a/nemreader/cli.py
+++ b/nemreader/cli.py
@@ -23,8 +23,13 @@ def cli():
 @click.option("-v", "--verbose", is_flag=True, help="Will print verbose messages.")
 def list_nmis(nemfile, verbose):
     """ Output list of NMIs in NEM file """
+
     if verbose:
-        logging.basicConfig(level="DEBUG", format=LOG_FORMAT)
+        log_level = "DEBUG"
+    else:
+        log_level = "WARNING"
+    logging.basicConfig(level=log_level, format=LOG_FORMAT)
+
     nmis = list(nmis_in_file(nemfile))
     click.echo("The following NMI[suffix] exist in this file:")
     for nmi, suffixes in nmis:
@@ -48,7 +53,10 @@ def output(nemfile, verbose, outdir):
     NEMFILE is the name of the file to parse.
     """
     if verbose:
-        logging.basicConfig(level="DEBUG", format=LOG_FORMAT)
+        log_level = "DEBUG"
+    else:
+        log_level = "WARNING"
+    logging.basicConfig(level=log_level, format=LOG_FORMAT)
     for fname in output_as_csv(nemfile, output_dir=outdir):
         click.echo(f"Created {fname}")
 
@@ -69,6 +77,9 @@ def output_daily(nemfile, verbose, outdir):
     NEMFILE is the name of the file to parse.
     """
     if verbose:
-        logging.basicConfig(level="DEBUG", format=LOG_FORMAT)
+        log_level = "DEBUG"
+    else:
+        log_level = "WARNING"
+    logging.basicConfig(level=log_level, format=LOG_FORMAT)
     fname = output_as_daily_csv(nemfile, output_dir=outdir)
     click.echo(f"Created {fname}")
diff --git a/nemreader/nem_reader.py b/nemreader/nem_reader.py
index 796d6ec..6261399 100644
--- a/nemreader/nem_reader.py
+++ b/nemreader/nem_reader.py
@@ -29,7 +29,7 @@ def read_nem_file(file_path: str, ignore_missing_header=False) -> NEMFile:
     """ Read in NEM file and return meter readings named tuple
 
     :param file_path: The NEM file to process
-    :param ignore_missing_header: Whether to continue parsing if missing header. 
+    :param ignore_missing_header: Whether to continue parsing if missing header.
                                   Will assume NEM12 format.
     :returns: The file that was created
     """
@@ -156,9 +156,13 @@ def parse_nem12_rows(
 
         elif record_indicator == 300:
             num_intervals = int(24 * 60 / nmi_d.interval_length)
-            assert len(row) > num_intervals, "Incomplete 300 Row in {}".format(
-                file_name
-            )
+            assert len(row) > 1, f"Invalid 300 Row in {file_name}"
+            if len(row) < num_intervals + 2:
+                record_date = row[1]
+                msg = "Skipping 300 record for %s %s %s. "
+                msg += "It does not have the expected %s intervals"
+                log.error(msg, nmi_d.nmi, nmi_d.nmi_suffix, record_date, num_intervals)
+                continue
             interval_record = parse_300_row(
                 row, nmi_d.interval_length, nmi_d.uom, nmi_d.meter_serial_number
             )
diff --git a/nemreader/version.py b/nemreader/version.py
index dd9b22c..7225152 100644
--- a/nemreader/version.py
+++ b/nemreader/version.py
@@ -1 +1 @@
-__version__ = "0.5.1"
+__version__ = "0.5.2"
diff --git a/tests/test_incomplete_files.py b/tests/test_incomplete_files.py
index a35647e..ce92764 100644
--- a/tests/test_incomplete_files.py
+++ b/tests/test_incomplete_files.py
@@ -4,7 +4,15 @@
 
 def test_correct_NMIs():
     meter_data = nr.read_nem_file(
-        "examples/unzipped/Example_NEM12_powercor.csv", ignore_missing_header=True
+        "examples/invalid/Example_NEM12_powercor.csv", ignore_missing_header=True
+    )
+    assert len(meter_data.readings) == 1
+    assert "VABD000163" in meter_data.readings
+
+
+def test_incomplete_interval_row():
+    meter_data = nr.read_nem_file(
+        "examples/invalid/Example_NEM12_incomplete_interval.csv"
     )
     assert len(meter_data.readings) == 1
     assert "VABD000163" in meter_data.readings
@@ -12,7 +20,7 @@ def test_correct_NMIs():
 
 def test_correct_channels():
     meter_data = nr.read_nem_file(
-        "examples/unzipped/Example_NEM12_powercor.csv", ignore_missing_header=True
+        "examples/invalid/Example_NEM12_powercor.csv", ignore_missing_header=True
     )
     readings = meter_data.readings["VABD000163"]
     assert len(readings) == 2
@@ -22,7 +30,7 @@ def test_correct_channels():
 
 def test_correct_records():
     meter_data = nr.read_nem_file(
-        "examples/unzipped/Example_NEM12_powercor.csv", ignore_missing_header=True
+        "examples/invalid/Example_NEM12_powercor.csv", ignore_missing_header=True
     )
     readings = meter_data.readings["VABD000163"]
 
@@ -38,11 +46,11 @@ def test_correct_records():
 
 def test_zipped_load():
     nr.read_nem_file(
-        "examples/nem12/Example_NEM12_powercor.csv.zip", ignore_missing_header=True
+        "examples/invalid/Example_NEM12_powercor.csv.zip", ignore_missing_header=True
     )
 
 
 def test_missing_fields():
     nr.read_nem_file(
-        "examples/unzipped/Example_NEM12_powercor.csv", ignore_missing_header=True
+        "examples/invalid/Example_NEM12_powercor.csv", ignore_missing_header=True
     )