From 72ee866258e6374ba56f7d29a185a08813d23a95 Mon Sep 17 00:00:00 2001
From: DylanAdlard <dylan.adlard@lmh.ox.ac.uk>
Date: Fri, 19 Dec 2025 12:39:21 +0200
Subject: [PATCH] verbose

---
 README.md                  |  16 +++------
 demo_files/output.pdf      | Bin 44306 -> 44088 bytes
 demo_files/output.txt      |  15 ++++----
 gui.py                     |  18 +++++++---
 src/ecoff_fitter/cli.py    |   4 +--
 src/ecoff_fitter/core.py   |  26 ++++++++++++++
 src/ecoff_fitter/report.py |  30 ++++++++++++++--
 tests/test_ecoff_fitter.py |  36 ++++++++++++++++++++
 tests/test_report.py       |  68 +++++++++++++++++++++++++++++++++----
 9 files changed, 180 insertions(+), 33 deletions(-)

diff --git a/README.md b/README.md
index 548bf7a..02fcc36 100644
--- a/README.md
+++ b/README.md
@@ -23,15 +23,9 @@ Demo input files are provided in `demo_files/` to illustrate basic use.
 
 ## 🛠 Installation
 
-### Install from PyPI
-
-```bash
-pip install ecoffitter
-```
-
 ### Install from source
 
-Assuming in project directory:
+Assuming in project directory (after git cloning):
 
 ```bash
 pip install -e .
@@ -127,8 +121,6 @@ If `--outfile` ends in `.txt`, the tool writes:
 - ECOFF estimate
 - fitted mean & variance (per component)
 - mixture weights (if multiple distributions)
-- percentiles
-- likelihood values
 
 ### 2. PDF Report
 
@@ -137,8 +129,8 @@ If `--outfile` ends in `.pdf`, the tool writes:
 - histogram of observed MICs
 - fitted distribution curve(s)
 - ECOFF location marker
-- table of model parameters
-- censoring diagnostics
+- fitted mean & variance (per component)
+- mixture weights (if multiple distributions)
 
 ## 🚀 Command-Line Usage
 
@@ -147,7 +139,7 @@ Once installed, you can call the CLI.
 Example using demo files:
 
 ```bash
-ecoff-fitter   --input demo_files/input.txt   --params demo_files/params.txt   --outfile demo_files/output.txt
+ecoff-fitter   --input demo_files/censored.txt   --params demo_files/params.txt   --outfile demo_files/output.txt
 ```
 
 Instead of using a parameter file, you can also specify parameters directly.
diff --git a/demo_files/output.pdf b/demo_files/output.pdf
index 4ba591a18c56d01c8ff9a5a3e73fb630a54e2778..acee096c9370f4a78921eefc352842e18b0adb54 100644
GIT binary patch
delta 14307
zcmZvCbzD{5)-~NogM?BNg23s+K|&;@yF=-a4(US*h%_kO-3@}2(nv^4hje$t$8+x+
z_l0l&wSVi(xyG7f&NbKBWACU*w7LnjvOrX%v<;nlRKV#mr@c8(O#N=4kEVvb_S#3n
zpH}Zoey0S8XFy@k?4~QGkjN&oT8(&3G0I*gxDVC~A12N?R~ymkATl;4<?!aZi_UPK
z(wmE|?xUlInWg37-qV*t*KOmQ-LgXa$5VfEO+6ZC5rT~~qiMS~Tkd9|sn$)myVj@3
zpF8`jIu|LNDcco^4Fk#t4m{Iw2JzBF-@Y`i^=XkSB`d$KY2MuXqjv2&`GQhA#j^O%
zjeTjagT?S%Y0hoJtZcXtQBLdhYuP`?E2JG~>$S;=(tmmZ=I*W5<5yew^cTX0yGwMQ
z+uigKy<m?c2TOsB(Zta{T40aRAJ3NpNdF7Ao(!bHf|3ql$W2IBgJ57(W%A5qN%a*Y
zA@3L^x<zt`d~18Vx7gprL<dJ;74y%rgFwi<##C*CH;*MQaK#UmQ!!{d!Xe1YazgyZ
z8w+E%+NRqH$k6`nn+EPjKh%(}K+yQC8}!`%RQx6ai_S6x6*?erxpQgI!nBx+9CO&w
zu2dTeRJqX@u5DPG+)8$EMNO+_?=Nwm3wy5>La%U5rmHCQ4Ux4n<tt>V_)JI2icNe*
zl3wZDP*$<&BjSYg+E;*65p}oWFx$%N&OzzgP&dDZyAQEGpkzg*DDAL%vr?;8+)k8t
z6I^l{i5Z(m8En@KwPY50Y@3PvvVIWl95}+e-nt_3%=i`v`-E#*<~wp5kDNK`vO|7$
zz)L?z@fb%KsHNyN%?P+QxnE^ltR7uqGfNHVZ{NT?7fn_TZXGXn*U$Bxk#0?w&ZZ5@
z)`tXoA1ml4q#RFop5YFT<k5>!>E-zF)UqpFCh7JumVVk6b`<xC3f@C*y&jvxr(j=D
z@tI0<Bbh0F74(jK^f+PrM+iC2x#R8Byah=z0+~^)oF!d_&oePMS8;CMR3d5$<52*Y
zUT^)thsozAre=+?PVLO{i@`e|`#CTl^OSjzxonCY8Ot>;Z3#2CT%{NY)Q*iKz4$zD
zyQjv(n=fvLl@r6X3Ehy`%?c&jbegLKs?_N_b}Q#N^$tI!Z0#%*0)FO%tc``LGE~{L
zqCX#@l$PHAHG96(=bo5<R(4{1DCTInlcjmNu$r<$gNC<Pp7c`q!|4fU)Q1VA$W{OI
zdGr_Un7iM5_m03cNp{`?uFlmIc}E$oObrr#Tc)<(SM&^zO(8hQRgXclM0<RJ5{QP;
zsTo5{z`^AgU-?xC?-IJ<O!)~NC81iyAq|H5>X+|@v4(l3kRVf!S5=3Nl)LIQFDl$2
zD6ae#Z913Vu4|v3x#p@qH5Zs$Nstg#CEfL45^B6Tp5fUY=kz3QB?;Mc=*Vf7rDll#
zk)$y!Tgtk<oZBpjObip}NXg-Giaf}kr56p<?A&P}Qf8My0<(P;x#&Gui_LN8`HtF_
z<d;R99FG}S2IW<C280&3n2W1?6J~z@4&X2+EGlH)Ryn5O;Eg6N$X$%xOOx`RKxDD6
zE8&6)#;GIJ<u$WpA8EZ3%|UB<1P*9tx14DlW-T`l`n~c6soT9|G0K)|5o=>$X-B}7
z=8jUFOTL0qk~10A2TH?9pPzo7sc2(Jc&3`uWBUdxOF=NS6+0sM>ARi`p>`*2Cfg=d
zavtS(<vNBqtR0hjB1AuI!YdEz5)L`@m!1vIb;`YqiSv!musBxvR+>Fo7@Qu_e=c|V
z20=9Q&3f^*FfyCL==Yu^y}-(oo~%K>HtO6{t^KsnmqS}IB0-wZYC$qz$C<b!XT$9Z
zYNPktKPv6PK3OFTCmF1=^I?;WCuot-3m3{fB1Gw;L&We2v4VO(lGAAuC6F%XE$;^e
z2S29VAF+3^{>UCjm&$ggFIHkp^n3Zbh#BjkNq*>aCKuB6*DF4P&qxw{QDE;E8<a<O
z+U-Kc_N{2LPe;^`=t94${TOk!`qUlAOo-YU6AFy~KELxUTyhq7#LlBwF+^#>(IdYy
zBo8qryew|!+|Yyw>#Ewqz{U$Vep}r52*^Q(^*j#NhGGpy#2{4GS-1D!Y%IAs+E%X>
zksCoo>~xTGXM^91dFtsu#gc9iD`K5~jYgO}xp>Tu`HCNMEJogp28Nd6KKBJ_RdfIP
zp)o9xCO^bJNWq<@1My_<cxEsheDTcFgG;2OB`uHBStE6tv*uB}$b{<?gNW}O_TKZF
zGNi-`i!b{YWORbm-O~9MMrE)mWF<0a;gKe?9^RtJ%8=^DO<pr5g|v(J#s(R<$r-6Z
z@HjJDk?Epa^o6v}l9x{6LHWD!{LYX~UM4b74m2W7)X04q1AAB|@&!)GV%i@LDcS^E
zlspL<;xUG&=9Wq|HtDZoD)$P)rMvY>i+r9v{i$p(GHq4%3G<l}PvbKOGDjf@uNmoJ
zHaATfl71<IUo_<{@b6Y27lnGXAGzbJ81zDzu1P6%Xr%>aQ6Th12%23=CiwmG8jYmc
z7eD$m40Z<g0g_)eT(|?|nXgo$nW-EGU{U-2wq!c?G~c{?y4&jF%+OhgmFqsaJ}ZdB
zCH&NumT{iY@`I&2@LDC?wDQ`+2Fb8ZuaD0rkAyz^gY!}<1CJrWDWaB`*A#zVR+%yn
z7O)m(LMba1feJ)<yeMO0wI)z&wj%D6(-VpoSNx(QGa-~iNaGzL(<}^X&=7}1URb2a
zYhVyYP&Dr+#o5k>Nv6In&%KuSoMR><#G9G!vuS?4Qh3Kg<nDYD>pBfb=Y)*aYkGdm
zJ*1VvBWGR%x(%pbzkrq<OPu&6u__U>Dt}=e+HE@?d<UgEfz7tEhaAO3-S2%FAPp4q
zk;ODZOaMKhMY$)>N4+dT`s~wFqYNEeIoU@!PT4sb6xPLP9;oqceJzr_Z(M_LKUWk6
zV~$Qk_Ae=eB@YYA=ynw^&pq3?IZGlZu#xUD^s>RTC+Z)}L{vRwrX`G^oF(M8)5Au9
zx1(Kh!+xTuFj0z_H_r)qjWv&T$ajM4Q6$z4*(EhULM13gUkGVNN@h{;PWU`(L%r4o
z0b%6SrPZ|3f!9WbgS2nalq%eXIyh;B6(QMASXipd?Gw@8_Qo49xkR5>dh=xgSCONR
z^vkc<e)5>wS>qX!eV2e!qHf}fNo2Fo2IazPj{HelvGdTo>_5wXjZWaTW}M${dCT`@
zE=z(O|2dzYgrB;Vv?|u<GIxQhUF)8tj2S*HNDoW(vG@1a7>l37K%*YDXQ`q|jNIt-
z-k&Nyv%{OeQNigvxhB$=vTKXH7z~hfox&|F%fhg~4B_up-h7W`m*dPU7Im7fSZBzH
z_}sx{wU=X-N76tc&R9o<oBEj|tb(w?Hcan(=;R=_su-SucDn*(#LN@KSd{)$5*^U;
zxq-wN-LA8k{F_1Ugi{T>g6ahygJE7kRY*UZs{GeVV|8`+Ue~i=H(XQ8>aBQWGnIf1
z$2j$^+LB8F+}1Mgs~JriamT1|b^+Sk6qH<rSb<|92p$tlnl+tD<lk2aI`%JEFGgTA
zK?U~sy_tv=x%W?WbO-@)IN3M~UR|GP8!9M9pVl$_-qYLnt1Y-#|INODu_|#2GTR>t
zT61G+q@h00`zaBGwwNaTLTbAmDJ!6yELt*u<VCD<;aq@X+>sIMJKIs&mz_>RfhADk
z>Nz?j=~tpcYtLJE_8MO0$4+gaM2sUp|AIp|IBq8}DlI5Za&+PM$uaqff%CH+f;FFW
zK`n*|@>d_k9rk4G%z1V=Vs#m_oVn?a72CA;p`p$gHZE~Yac^ut3c#eakanNwt1x%W
z`tjUA+n?<e^zC8;irRFu2SE9&9?bn2-m(L&qG5F&@AKlAJ{3h&DpF6r<C$ngZiqYq
zfR?QIa0<9_w{MrLjDL~)_K}Y?$z=QvO#gw|=f-C>W(oNy5i7}**1^s|n(4(MV4w9=
za&q7B>8rm{;FS1wDGBYG*GMKioqkT6$RP2LEeCcXV$X1XBgS-Ff(Q<n@hD%A?shTe
z5|X#_fn+MJO>iK*Riv{CM-=4HOf)NeF^KAX@Tr8xh{~C`9+!fs-y>~We=AI}MQoLh
zzu{?J6(^4a{STe?Dxd_BaljpHOwPwl6JS7O4+(*}8L_RR(2w*yH6)GRlW<$J!$+-n
zf#eC<pYwarMBLoaR;Pd7-DD9a_3tWL`Y_sDJJ5jxDG2hz-l_6^Oh++%zJXB#cm}sK
z1@4spV4M%8bFxC|?|UXb!fdHmA;z}$uFtQhq~%Q}_Ok1gZ!~e3V~8=noS3@;?3yxZ
zKA&WR(&nuhkzB&l79{@?p$Va*((r&{C`p^I9BGKsOwj3CmM}6$`cT_(C@jhE<!iR1
zsh9EsSo`u9Enl9Pf&f=$IU|e3%q&7B0txKbK=}*9l%r1m^|^N`RaB`mG4DI0gQNO<
z<*=m;`yKfeX!HcA*{P#a31TY52D@LGk90@-C;_kdpR#9P$06UnR45S4Nj#0m_&dgm
z*}fIW7CGqq-qOofqly)E{wI;}e0CwY4!$DeHXDE*dStw4|4Y}FK7I7Pln9YdX9QL%
z%dtjI)morc9+8Th0xsFUR+SRPD2r<KggwrD;<U?B?q+_VWrrsU<G3|-12LEL{_hDE
zYI(kt`WWI4fC=)YQ*f*OuFJ-2AkE_T2n!7lUwr*DW<oaQuh%XcMjHl_ULdFEJ7f9M
zw9E-h`Mr-BljU?f04qitpf)yTD&<SZEl>`S&G2h^SNBAF4lWVGx%^#cXPEB$hyxpt
zvnHYPm#>92guUCPn#J@Uv$r+8R0ykbWO(nCcqx1!7J>Y{!`wZxwahRGv3gn*xP3Tp
zs_)5X?)KKk*W;isdZ0O>>FV(C`0$<kua~+Paq2&(k~%K86r*Db;A9yboslZ>v(}Lb
zp(p)xiQF@Dk;y&orS7MMEU>CFSI-=4!t3E-o2Dy!tFQf1WCui^THk9_WGLL+3^tkO
zv>LwH^dBK>rQc^`{5AM6XTCkWS$q58dS_#XC%U=MMCdv!7tK=Sjl@N@?fX^X?@$Mp
zn`^ZV<CZ4Rh7CF9pHHwag@=DtHxj=rTN{}dYpQnPy1iz#T%?!@fOT3$VfB7p`qrKN
zMG>A>Az?Cu1BY?_&;}nS3G4A_#)8K@Jg+x4`hHs?p_^2u+qZMmhm!{eQ3>R~gWst(
zh&4e!Gf@&99P@F0YkckQVskRfth4Qu`+M$c;Igp!c4=euXsJQ}Y@t2z^uu@NIEMCo
zOj_c66^zD|-|EtV(c|tZo=&v>shZJY%vn`lrfZJBCB9yLvlU{;_`E!AJ#3_K9j_1!
zX?(Big!HGVa2-PoAFW&*AKfg6z|~V_!`UBuUk3DAx}oeM>eib!M;jYl^x`SRk<RHf
z8>I$~13l)GJhfMRd;&^+gNo_B3q!XN(Qk5%FH`%r)n>K>)6#La`lqaT=H!ab*Fp+Y
z^{!S+=Nes3Z#FkSjE8%9o@R!~r^xB+$j_yWBj<9G64gzF<ZgF`&KEBF3Qy;dpIC&s
zoUVDLXK&9gwnOl@(){pd^j+!~hK4CPaS%6T^6Jq>r*Cr$->rIT9Hi$Ko{XIbB6J0k
z>WDJcoX&nMjWYXRl5BifDG8^i-lt0C?z2|c-Y-Syv>d!ss&x|})5tulZMenup#l8H
zdunWL<lty;{LTvR?#{-*0uKmc05JS>0|6o+(Tk*Q(VtKsA@#hCa1g(o-a2|xb>Qnl
zx`?Az41D$nP7hbs;|n%+*3D5rd>6(`U(dQo`V}2-(-iM?bhbV!)O>xhF_7vB61aQ%
zb~!&fZb(jTS{f+qmYv&Q&Fz*?-CrH;-(T-5!MeVlym?f-KioHmoZD{tBXAQZF}L4l
z+OV~<a~N|qjl|cNDgTunsMD)=<i>@U-~XI(3b}bwS$uI_rTL<Ue_X5bdZmGTAYb;j
zF5$-1?WMpCZyc5bKGb+TQTp*xZmP32w=vJptZnoL-NQ;AkF{DFcly;eOZgkMC9;g4
zyYW4#TTO9KU#~NJ4)T@TJ}>zwwkf94LzmFI(-K#G#)uph8lgD0{FfFXb}jZ-I#IVU
zy{1@`_FRO<Hg7vo*VeOyr7g~eP$T&DUwlO1M9oL5zeHxQ_Vg~l=un|6B_)Q$W;H?I
zlAJ57Olp_1>Nhn;3T2@NP8z)Ga~Np4WMNTf^0E+EZGu{|cuFo4K|L!&O8RCIU#a{x
zLl(E{k*C`&In`&Vqe4Rz$I*!*7EN4hm6T#XvNfWtL+c|2vkn3iGCT+!wpcDO-PJPC
z20vgAz5lluA)=dK_Kz4buJrC-V#H;iP4T~E4QAb4W}K%4R*{cSwJTu<+UV9JAl;_n
z9b#%Wfr%^cxGkX$ao{^bwTO$kL!>@jz}`cs0}(Ljfzq@|2jrC3Z~o(X;Bk7U9F8`8
zX!Ds)0Qv25^T;PLIpq51y3O#R(3aP&$7y@NbME8{^>+eeyZ=-pi2hL{*vRud`%6|L
z)Z0(g{;fvbi4nX1h!Gk-$_kA^J8BGVfUg^FeJArsF+6{d@3+tEv_ws><g1)MIR>B1
zV_sg8CDN~MN;#Sw4h55Jn-b}jx$@0ku~)<XIG?UOuSaw|PNL06`Yb25XutgZ6>;4E
zj)(?z<&=E$TeTqa#OPwf4_^B8cE6A)+=}l#&l=U$>u+*&M#*3mWB>Y9z<|VHGOd~~
zeyK18(i*NS$_7{S8ElH!ib=~WuOc7VO5UCQUK{_3$Xp+$FUsb?);AgMuw!1rn1({%
ziy*3RsMU}BK6MRm+f70fY*EtL>u;c&L=g7V-4i3s9X$FPJ2nQ)i2IX;Gc9|#j_QWq
zM$TmVpiuSfYU2qe56}NXQa|1dt6nWLD~n1VXQ<<={i|DCm!9&QAq9b+k=$&#@0I<J
zvVk7anJNCOI&})^jyAVj)x%u8I;C?Mlqa>AEdQ5nT&U@KVa6>xRm`@H=_;vwy;K=X
zG6%pm<~ZE@#{GPKAXRn>h*vzVqr<rH)uJ}YdX=_;2buEv-Gv|WImKs-^LenjtTf7a
zI@{A`4*G)L;($02A+`Ot&ZzCgd8p)}VbCAH(D31&v3jHD2v$-ysU#dDdYND5yy>6X
za56Pn)K=CAgcGtD(Z5I}^7{W#l=qFikrSZzER%PzuW?L}_LU`i-ZJFz$t;)ogHc@Y
zrA)_a+c(O1w)vkoNNNFd=A4mVtkM|5az|k+Zc-hN8HcP`0KTMJc?Mt05+`MvXO^*5
zzl4wVtiJi55qJlz<C_5T@lE(+fnpTP&&U7N1BTHHNIi~lt8yJ#={pvL8QpVdzvhd|
z;9h~suVf;vl=CZXO6tmMq*N_JO(<G=(jxOEeA|O&^kz83kR!D7sH|WK%pAGUXHJew
zClKtS=oelVR*&Aw$y~<2VMdcqb&Ge9XABi{!J_}|-<?UpK+R0KSmDC=gvO)Gf~&-s
zw%h2J*`q_r*Ca6-rvru(N`o>_ceCOgsvXi#P~{Uwj(7`|k8242cI*e+x0qYz3VM8=
zd4|e{RVD_kqC;-eanj@M%)9<z*q^!fm((S|Cc^*=jF?%Mv`OS;zv&+iEvut$a7?v$
zA<mT?Zrlc%mx=R}Vr2c1(O{8iUzA=g{x(`-@XKU0m&nO8@4OFwRXIHdy;W&rfonDS
zk1PVN4@pcZy=vZt8@7I@swKcK7I`|@%`ay9g|{V6CIx|10aLv9el1)kOGoqFqe#!B
zh>cVxF~)-Wk&UnEZsBwJi7}E$3^CR`u{GTO$CrYFT5MJJWg-gJDzCK;9=%e1+RhGQ
za$fACn`*_O6Z0;kv>O{nL7AR5+--~J7yE|)Q#qL(ZLk6X`@rPL9~R&e-RN7X@VvdE
z^icPL?1KLX^6n9|9ShEF)}7N7G2EXwn%J$0750Yja@t!duo0Halxb+D*apTisQ894
zlJxSAs_X2gt66IVydtZdcwNB+#}?*~QfGE~Zj-1mQLR=0FE^+amGd`jhtHr1WBdA?
z*)pl1MJ{&^C+WXADkXcD$T>aPfndP<lW&2QiqT(*#BE(pGkq_L52$Ahn-q6A^+ilj
zWCp9UvsWl%ZL_zpDB(xHKHScgiSu1(9YK-P$v8)A4XR4$=y{XOw{FhP9tS6HGh?dZ
zobt!1;|N35Vv>JZEp@gLz4TdDCZKd4+C@0kr2!bT9;>Dj*@@{UR&KM<uUvU1l4(Jc
zi{0jre2S;=gS+aHt4w1V<Gea8tLj-+fOgCLt3i`Oc#D<v$_KgJnjI*UGEujUJ5819
zHYy6S1hXMbfq7G!q=<?2HuIeu1Bs<?Yq+AA!_rlTl$-?~pm5U2b}H<JoD9%<Dq3If
zm`#30K9B6>yCS9C3KV+T$;Ry2N%;*+5zCjfNR2mDqSziSG(TqW+!zLmTf|XRxRrKG
zFfuHjb7y!*ap0d#eQOwePl8H`!K0Z>C7$pS0u@<tBCe}K#W;Zs7)9X{kBJIhl2-&j
zSGBc>#9a33^ZGo?u4*Kz@}A<LErW&Q1vAPJdGKGIauiI#bwo>eLFLEe%-Uews>yR{
zMP!%|VK0t=Bc?%AtlL_POI?uu0vYEqCIP4=v`uYndZ2J*Sd-Nc2c#W@=|xFv1@W8G
zXsb6QMSXQGB@-V!^6jriDL!6T>Q;wK*j2iOy|G2j8yZ2gPg0Aea%tQEb*@En_(4=S
zTDB6yX-<&05Y8f-9bAHmML9cPuOX}UIdVjYS4OHPpSY=^WJq?+2lY1*f=vanYM_DB
zUiyVng#bt$hi=bH<?|kwCSXtTs{5~0L^;76A&ucgSJ;6h)qDFU=^t)3_6d0_9iXx&
zs|MYyJg=$qQ<IecKAQ`i#^;fdLsy9wM%YWaLH<f{P>AK4Z#c@$zE)utoe#g3d5Rp3
zcot3Olzy4v)Z#)NC0mj(<MZ3xb{0U)XAXW9g)?NgcxK%x5{HleF5da7?Mv;REQdR8
zs+wyEI?u|1D9HC3lTu=tV3aml%BScfvmc2|G(SG&mWp7xQ%^aLfm&6fM*HiS^g<1V
zTW#7sbOah2uKa7*n}nJ?6-XJg#3gcOQ9oZc|6RCZi>fBpkHlLCT<7Ao5FaXC)!B+i
z=+1HKa~1eLUrq9{w(wXyL3<fb58A%y<h@zYR-=m{8`HpPNQ+`51iqUotPrsc&@(PV
zB~myX@>i0}>fOkuQi3R!Wu?<5bO;Few6bGYghrhhf2zT0ZZ!>1f~aM>Yepg|@!(gb
z_50jmC;=5X!&sqon7aI$f3td0Rl0cff)pDQEVx3Y)HeV2%(H4m9=r)}{Wj+HQ`Xu%
zC?R;d$V$6BVyFAuL<Pynn81tb@eV-wL#ixd<U<668K+$~VW<a36_eL}!k%nCJ(yxt
z0P)T>w;C!D>~%{sVEIaZj<FmG@~AUectJY+xw;%jYl>ga($+&@CtfB{jO6nqi8b-X
zm+mxueA|syqXH5N;gUER>7NM!^IgssYkLH8m?lrkyAVEeGlahAPR<()2<1A-3U$fW
zuE(Cvsb&ObE>e=(tRg<aWz`EjeJhM9Y1(5iS89mV4pL`gftmvLFMh5gM^ZjrC-W+c
z;5`9}%_j~$df)Aq+OX@<N>kv?XBDdwVGbCgc!SowmcH57D=C*g<YjAQL59sZPHRoM
z=qL5^Pjtty#jNdz+OSmC(JAxjwP&IA5CksOnZRwv74fqU;9Kl5pFp1=r_>QRk669|
zIVE&GgNH)&v4OLCWH34pvfyy3k7H^O?GzXzF~E<5>rJNYv`!azNnt({$BmiDoxhnF
zV_mzHJ=Xzi4z=<+pXk`)e>yb{tOi;;L#oX`a!V-Fx?IQopp&su*!?@g`gycgBr$18
z!{}v|k4?R@2rXSBJ+9S9*xNHNw+iY>l{%erGTz6N-ZOX4nIefPNagSG{Y4*P<uQ~w
zOX+>vp(%dpvK{SgIj}X2Pzrh(E1ct)P8Sn_Y9ClkBP=mrR?1Zp4yL_$3&)b1iYFxe
z<qhw9&taaME)f5be#*-v8X}cGnhZf<Y%`!fc+H1bx(yY<-n`UErDIM>R0tK!#_0;Z
znkZ)2<j>NQ&%_-<LfP$S_|T7y2?uJpt`*OgE`ikA`}E6{MvA-mwaSTkzETERT<D~%
zdtH0ypU_Cojc~Kx+A>s#Dhr=1cg92~ED*wnRM^XnXzXLu@+sP&{f5LZGT?buQeP~O
z`Vq4-$a$?X6dd;M<W`76Xd2n!FeZOBrr)%GfuMA#&tE+xF(M^!E4r#J#I|nxk|Z?5
zdm$<_xKfG#eejpRP2Rjty;VV+Bc-0?zXB}YJR06V_O^7@;DnSAI#E_EZ)?%c&(yw`
zz9uhX`N%b+^pYY5-7vlgS2aZI<Hm1u=|^#~QtAtjPmIDKxyHL419WFb-jwVZ5h8p%
z)om)rjpdxA$k<VQnsYlhc`+8B9UQa{FHSvZfi%0eUi<i4prLiyis}ZT1I>5i_Cg%2
zq9+n{3#w&HP4wVmS5@@v&d>H@MI)q_6XC%W!Xy$etb-_i@sE1R6!DHJx8V`c;#L}~
z8I9uvb5)Y=3W*iFJ%QLs7jt}z2+xmgGBrf8k5TC?MJ}1Ur2+L8zbEI2j@2p>Y`s&p
z{5I}v)BtorU(HWU-cp_tR|~a^Y#ftNtt>Z}BPFV%o9h+$D!dO!CIx9ok+4qKiY|xm
zNP8xdtn%+jFObaZWdb7gK1I=^GpAcIVg7wn|LCP^UlK}J!HtaJoX$qyuV4i;0Ea9r
z{qn9`A<xy>*2Gl2&###*#(ZgB_p#>5<cR+k5JIV1AB*qX>id%6UgnIM*0t`x6o>H9
z1t5tHrWyDBrj@c^eoyu;z_W+v48^z>OQH~ypH0UCMK7%W9!+(bU&x0gsrh9s)6CT{
z@_DA)Tz|U0L58%5$Z2oc5V%dPjjlXR%AmdKLoQ)O&&p@@QG}H>y6$Z?rwCs5R_5Mb
z@6+eXrwZKKw&7r`_=V>nt+<2!H2SU&xnsUAN$L<v2ufg4*9{;ce<FolM<?B<(va1i
zBOUGfwgl@>LfhO?v7&JOw^G&Glbf6FRKaUUtFJmaR1`v|Za*fKe#momMiz&7SZ-@I
z*&(;(zfMR&d)k{kbF1Y`=C0<BrwVI@WNu>myxhXSf-Obz!^zui%~p$DELWtiZqznT
z%;)r08-}}V0}YX5uRd7uJjO=lgSkZ7aOw5V%^Bswe(9PyOtYGoc_8Qprp@C&_8V~1
zR$BYhCQ>eR|Kc}|(5lhu*El>s?kyb0K1J4z!K=6FU1uERy@2_DLMG}4emK8wf)?c~
z%yf-!7X5WLy*$Y{Y2^NwOoY6c`=EY%wmRCoK2<BU_cLUCcc(5Ufk9=&WVO<3s)cbg
zWq*5UVNLK4<tS15)vw7h#oS}#>ZNvi-Gj{^*P=-lFKeeD)J91EUUREM<IPD-<RkjY
zKOzY^<Kv|{AM~4DJdb}LN$Td6QgMg7FBJ6xXHLJHNDm%W_7s&}Wa5vq)1Z!A7A56e
z+xOm%aXkv6xsqXbxLn(1HH2|Wp6ij?sLp!f)Z3FM>6$uj@s+P+=Bm#%p7n<7H@RHy
z9B$U$sYD)^#UZKu@hN)p`bpzL$SKZsd)VZx#BHhI*_~9?`$69`h(+&n^9@LP*7I_s
z?@lHPTHL;Qz57wRk<-~bv^nELp#i@@%m=7JGm*5H&j)jrryUJ$RQB}Z#T)wpJK3xs
z?=+%od-yZ;r(F%QA4ack)i*W@Cj5A7@mPCk+M6e?X@!;w4*yYw?)0JPC>APISabo)
z4`dX1C&kCpYw~aYm}0X$qArZDh}r!qa9;5EQ#}@cY{30R97mw3Ia+hrX?<e@tJf;f
zkMUzXpK`2g$?dtk+MK7$$yCF*>6}ng=~Mo>=_e(6`Zk~422QoJN8Ou*I5M9eUOQ)g
ztDktrDKxMtp9pTAxw`5!IUHWuQJ;3}dlHU>3x8<GSe8|1B66c@(0DIdPL&)kHkRA(
zZgj0YwOPM0U{Mj*n=bu4;3C}AY0X2eIVr!nK-%e)QyaQgAuG5RRo(eGQ_$Lx%J5DP
zS|?1Dg(U~;#1waQ)b+=!eE(y?PpbK~{;{u<!R6kb=;wn|_`t~XkvT=mYXSVuyyz#$
zBYM$}eaq*F$$VY1^XG(dci)rX>eJX%(-J%uEq9z`G3^~GDL;$Z7;tKw(9o?XFLroW
zc;jo<+C^4i6&mRn%1CJ8K5ncdsc18*N8CrxUp(m$pzt!F1VdSd+EIsiTgzUmdKFLi
zeNun&Z>>O7BT7;HFG&j9PXZLCHv2E$BM-c}q%(A-<?>H-7)Ru499a0k<Wo^*<59{5
z0_s7t`JPnYJqQ~e)T`ch<BrZMt||2Dmm3wHgo0H4PAk~_MSFMNBbXmwXZJdyjlIQG
zyfD2P-w0cAUFEZ}f>F9f@&y8g8=EGkk>%emgPL+jLsB`qUH456GwZWo7VCa1M&9J;
zc~qWXw6FZh*^w(GcrtIgM4S?7cDpko{8bw3I&<1335%7RJwl4g5ziUm`AkcZ;bVVH
znpFf+f`fRw1WkJWWzY%d+vIJ7A(jBQB}Kf6vW8nyI>L<DvuLvbL0a9${h;f<A1e)`
z2-4#Bk9f+(PYU5r%fHRdc;vdVLXkB{hh4|jq)M$?U!fldZx%!S%rXW=l63uWY(BR)
zRMm_Pvk{H-YE#+ldR-ozl*k?eQD<$ioU7)$Mlf9@q@Dc0<6yJhk3eS%l$o-12CeDG
zVfsdMX2JtW$~7^=J+scYD`I0H;qVtvMGF&64b=Hf+|_?+oovnN%yWpuN+1)7GDiDa
zc>!-4whX)Mj8E$-UiOV@=k|5&Gyl1`iQ0nEDDf;qHxGMmS{Nv^cwxv0F0`+^6(NZ3
zP%M+TpmY791oXV{Yns1kPA7xkBC-6oys?gUWJufxzMG6gG~NhVdYDGnn3MZOzSl*!
z#8I=Cv*pi$5~|rpoY(atHjs9tVMf@Hx1yt(pIA^a3fXgV0j!(Tkrc8!N_jUJ1ly>@
zuC~Vu2jiRjQmJD3{;hG`Z%4y9Pz-vUQ=rvunMORBO~@h>31FT%)X?>*U1+hHyLg4`
zjVqq!v8-y3n!RZeZe%%qt{HXCl=>pmdIVU`G-p=P&dd5W*_SMi*GS4t3OQ_}Jxe)u
zypH7-vlKClbS82Cl3~5ddK`0BB@I=YWU^1Ksp+y>c$pCy5pv0*QsbQ6FmsbK=Pk@=
z3}F#n6u;2Z$elnorae+onslf9uKLCAf1}KY@??oufvF8BLd50;9^Kj$uRL5bWo7Zs
zz$YIi@qo3x!t&k8WXLB`Xbki1zZ+P&qSE!m%vJxyas(8Lp&7>zIH_Ig$iGbswC8!L
z{zOE#)y!k#2}SG9D{?y=Rhb*Ja9OF9y$C6s<Hw?hW(!i9)<j=eW8Hu+-MCRSN5<3Z
zyUgr`42w%yR0(e9x=<(PO>%XxzATe*tlZqh=`=JGFEMbX7waP(og4W-8u;%93C_lC
zaCs4X(s~Z3=*X$MbFW!sp)BL3w&~V)drJIDR?a{PUMwzMWEdFwtU_rLVp(9rz>X0T
zP%y4s0d4CT<c;Qe_1(W~!D%W)`d3QVh)adkFyI5}XThVZNx54oQwDm$YW{Z=!y%fQ
zY;#csWl2F?S-slGy&~hcMoDY!?P4OpQAcRD*yJLv_OlO`*2AF%UT2|Y^WkbZ?NC`_
zwrLx!K8yu4CMR8=ukC~5Ukg>S>zv}&7y^|(51AH*Qi^Zne{Hh=^K*=ugW4o-t&Fc~
z`_l~yxe^xL52OnofOIYW(WjjAXgGGd#2|1m6<R%Q8<?$wl#OK|m|55G2}HE0D2wTU
zsI{&sCMuh!oER-VbUagBN?6tQvZ6yU)N`JA6e{z`iTW-wrDVrQF-el~*2u+t=~NE3
zo{`3nYYiBP`oOQ&PAp##U5}@<`zNagOF<-qrlM@b{N-q!9jZzzn(FJuU$d%>IR3*F
zts~?tqsoPUknCrK#HO)%dXvt&S}roUY0(CbMT4vEcKsmcp1#RWtOTr%%}VqM;i9$f
z+5#u?RfUr`Kosh2krLRC3?Eiz3NeqCwH4O3-z<ok?2lIL*>AK@V?TX*ec(|^XD1%V
zk=2`LkjSO+D#1e8(*L{l;B$#4m7?O`T(IGGH@gvG2Bd5Te%Sbd^jW5cJ|}h?s#0Oh
z%P8ZL{QMZE6k<E9yjTAF-+J>cg!$|%JTH`S?Rc%?ua@&F9;1E*oUQbc2>MWFxr<bi
z2$DE^7b9H#<$2>@NLpj3zKYN}AWwZ2tt8J+BZxppVL|O#YpYbAz>8fs9rh-j$o)yh
zIk*+nMTm?uzre4JK*ADlX1Ur}$^#ya!#+bmEhwP3?Cw))4Aml7=P_piYI<5_krMN~
zy$vK0oti_|pP~t#{{8&r<Db7r?6UYV=k0@S1?&9MohH9lopG%4GCy?@PS9?0us@3W
zsz(#*T}tX+9PQA_kqaM01Bl=Srb4sV>nE@XFw~H8{Z$8!Q5wI_f&!%9y@sY#`KOt>
zjz=u7J%<py{^j*(@dZn#WbV`uNZq|aUrO!~zE*Su<CLk5h4v4QxyNXWe=Rqou-ag<
zS}sv93fsk&1HUl4uuZtFJ^Rx+BrC)XGF!ihCNFeys%?|PBKeS!U|BV-rnM8GaF{i7
zaKPh%JjI{8E<K!KsmH&_#jH=v)XZp5v7OCs;LRVMj@vPrS=x=+HW)*jA_P(h?6P8H
zJ1?1oah0fIwu<5h=X_Y{P5Ll6Vr>^g>ugwyWSDW_x=T6nT?;k<uiE%|^QpDYtoXOc
z8k5|_SL>frfhOi`@fg5EXVQuM7gc0ZjNDGh=)uD@S9t4W-H5Ecxsl_n+^0}@sNs;D
z!5sd|^|KnB-4n9f2?guPU$=@P=O8kVc-)&0CcY%Y^m?Xph>nAhzJ0s7g#k6Y=r3Fs
z@Z22hBj*D_KxN8VMq@WBsnPlG=XuE~zJlk6KiP!%UC+LjHmQcEe6K6cn(7Wo)jl0T
z?nc-8AAeSU;@+|D-*$Gw2zb>Jb#NdT${#O15V?YB4y|s=Fd6o#?faE2iTbo5+xM+b
zP4uanxy9}{QIq?T+54~lGLN1Mp5{5e&xa5WIdqN*3syV7NcBHi0@e)mzb=yBtv+cT
zZ5Jx#9^gi*nsl3F3NE}gdC%_8?CSafjI4IxOi4hHBI{#wLK7sY%}+5-Dp_~R8x`K9
z9+$NGd06ag+^n-&y4qxiSN)L<85whnm*;oq0u!Lsxmb*a+8PWq=oI~e`Z1s+E7UfQ
z@i)3TI19QLLT-|zt_L@K<S_|ZjTr;lk5%KS{wn_j1=Yrek5hMVnyS~(Yakx~K>pTq
z553vh9XD;hay}7nY`BWm*Gtylbw%u&H-A>I{vh31c2FJx(wf})K9!JhLD!=`_T^Ne
z%CVv5EZMLrY&tx?1dv(0m7@}{<zC=r%lf3Ygbx0!>!Kh!hOgPJutdL^PpBbxr)B(Y
zs06pX%M-(~G1<VKz$WUC0jkF08c5GwN$3V~+#}6=mFDR^Q#JxpYXy#;_cqtcP8GDm
zJp-yr4IbA4>mr%2JgW=-sFtrJ2b%8cU+*oP7xr%Es()4gAgFWqmhA@pcuz_|dz}8q
zPXc=N%5CK%r0ZsG;URpS$Hp)?GW<~hav*B+gAI?~lXuwYp-()m4F&2?Cf20zUbb`y
z!Rludl*QEh0gv8wuS%#(*vOOVBJ}wE6*Ij&j9iw<_W-Xw4IIq}**;mD+}N)Osn+u)
z*hkH?GP8X*BuCbqQ)N)tA4OAxXfeP50HBf=)M$&10D{0U5G494L)H|805aVD001yp
zGvNK}1_A^Cc^{6!9@+t6fCqL!-iLM|04TaF6M+tcz@oP@S&&c!;+`J>0*8R_jR8S0
z80_8{00INu#{?iC0O)-?1Qc=~<8D6Qf6a&dk6`bdLH9xK3_<rN1VIqs`@kRw9CUxt
zAP5i!xwjY)6v6wSc6ayy01yNRJOt*2Jy<mu^w16r1>Oe%g5l8n%mcv)zymu7013T&
z@&C{DKPd-7!0?A-5X5~JfOquYFAX3Bj=0YZ@Q$qqfj~hZ_`h(G{{ajI!ybS@p@;{o
zhQh%QfMLJ~b$UlW_`!TIDDVGS;6MBWVZ6K#K;Zu$s^KvBLpucGLF#ycum=qCLg5eX
z;Gnzs@6pf8`;aaK5dOgKuF@Y&h~Nd>&pQYJK-^!<-5!Men|=WDfAS6j+;M$xF(3e(
z_klAI@BkPD1cDy;0`L0D{rQ0KdqsW6Kj`lJ(FcYg5bymQgFs;T!!hWC)PX<<#D7Ki
zzrlgQpa;x=p)g+L{S|`Y!23l50`oq|9T)-m=aZR#uKxi6BR~(f7?>CKAe|5Z1p42p
zL;iQC1;L;Xsucu--%bC2LGD@w_`jY1#}^KU+%HTJ9QvU9g5YrQJ!X&q5S*9yLBzaw
z8NCn83x3ezLA($90*Ds|c{tz09soiB;Bbb2^}qk65Oinv02uNAVFrPK+-Cs{06l1C
zV8C6X?gat{!2frbAph&20|W2s@V?<)>v}i_gWOLE7<fm+eP0j|{-1sRpCp1o;0LvU
z1m7j<UQX@|d3hg<fdLP;A{hLz2Y|sa$o*Y%XZGMc07Jl#JEs4||NmAn1oB|(fFZoV
z2a!Ml;0LrqLC^=lQ1Jgb?f%6A6#gHqBLDqA0}MqxsOr1+@sI@=>_43U%NGWJ2m(X!
xK1c~1@?bZCU%}xI&U`R0;2{gVcPIC~97FCp95lM6ND}WZbUYRoF?n&k{|8>OC#L`a

delta 14483
zcmZv>by$_{(mt#RNH-!S-4ZK^B_$;x4bmOb-FZ{eoeBuj4I&`j-Q6kO9nv5Bd4GGq
z?}l$3>yLHJn(LaGb7szaaZ)t}7d->_GvIklI}SGB<n1Rn^;|#)uK{)CODLjdV3d7n
z)6MO5ol8wlevOf<T8rcyLUeoIs$y!43;3j;vn-qpa6CI8YpA{2DWe-5wQ`I|NNsSd
zJQ<C+abvIE61}^9cezD)ueC?QU%xQ^qh#`&y}ZLIpya-6)}ztwh9RP1CCI3D*C-)p
zbwdmIs}XEiT)eJ!p1xB-pUL6Y-4S5vS<Vh5PJC>S--zD%+e@eUt@X5|V5>2Qb96<<
zcaN`9R>85p+K<BX%3)(&{v6yl&G&Pboj)Ljs*yW)$0J!G6t-YoE9Q=7XU%{Rln`Or
zMt^5;advuL;8{qYpTo1U{qZPdKcuUc&u9Y(h-XTgSPoFQQ~6Y068B@XB0%I?WEEz)
zVOMb0m0xl5XN-U1?2N@aoS8xu&WHx3v>#NR(mGKXcGEnDxiH8Mi`N{Jv9XM6w2%X^
zFJ+A;KXz#yuZ<F4w&nqaOxu%!uVeSG)9)yM9mp@tlp~QzAL%n|AEduHv~_gpo~{7s
zZR@5tieOAzzImSMboq2I;XW+J`)P)x7ZDr>iE#dtnL=#I^Jk|wnz7o#2<Arjb?hcI
z9$e!QXPG~0pc^9!5+4OLg*V?@D{fT2Dr(2fz6r{{D1Ck~zboOCy}tdBO4!xQ3dsbr
z3hvByaZHZI=*S~nouqQA)P!Sx;7FU3*Vg@A2W(@8#|Al5B(!JLsl-Vp<=_oaJObmr
zUoajzNvLg<AxjGlxn0>J8$D^<y*qWkI2$<6&R>W5<MiayvZmKMT}X?(odCb7lYF6V
zDfi>Nx++pgD{h>4m=HmNw(C2u;}LeOb$*k@P+j1@68ZbIz+^4EF<$}(0PA?`uCjTI
z^4$ucCzX`2dAfvfeXY*H-#;F=*|(mhA8OQ?{*rZy=ewn6v`Ugjt=kD6DOx5%<;Dtu
zN&|9T1+Is8JdKR)Ym;Z&Zoy}BYT9>~Tc>-XCEC1q%d-(O7Vk`{8DXxzM3p4^hWyGO
zez=7nyWh-JXs+?-3GK9o1DrnUoM$)k-1z1)DV=%}cYO_)17U7S1tr!;N#6;&i(na%
z<KYS$wdB1Y%{I<EAWjwydbR=coU8kD=X}CL8NS*q{??ZFwPAEZ8-96AJId*9=SD%4
zhb*~RgS%t{Q;Y8znqE?Rmfcg<7~LZh*`8Mn5vR`ITVDI-1cdVZ1a=#7zAy57Su=b_
zDP_Xzt-XkFnkD-IYsc*>xFs0oz4nFl$86an>}d?N7v1vjl@CxAZO+cQ7E*26@~76w
zEZ!ZjZ#eD5@p<(Hg8sr<U*GW=LY$;W`KQ*_GQZ!PP8y94JGYGv*mOC3(jXY_U)8(l
z56H#MZu~0IEEb2E0yyjaK`Rh9yw5ona;UDrvfIo-Db69ZIq$wqk^MEymARB5ZWcN~
zGBcfJ#JCVUl5DHEPPrD8Y||thLW+Xgl_@RvK25kbaJ@(`Vn#*unJ(5*o>!j03Ve`w
z5?*w`>7Q|l7v}XmA2<Vyc4AAB)o|*9uKSKEpoSvc&v`m<(}5Qh879)PblO2XW>gB$
zfwWZyg&BJA{FO<1j$)Tkl5GHZ`zsAblORf#%$IRND-$tX_gp%himunI*%iyhIQt^b
zGV0+CI0~e8%M4c{Q^r<suY>-C9GXPFam7FR83`4$u0liJP<0mO$P%Zb(;Lp&Ka!AX
zj$fMc*csc%a{@?)6f9{i?IJR5=4C#&s3a|9RI+IP39>^>@{tk@)=TjDSRF^klvS5h
zcZ5YJa84~<Rr_kM3RLXSIqa1vCjbl2_gSLJX0f5UIZO-+niV4trg*0>$G89TbrACu
z$^_W$?)1o=E_YK-_xOt#Ty%!Xi}&Y@j8>AP8@{g_n+}13_)Y>XU#u8A_sH}GsfDs5
z+ip}Fu{s-$FOh}l<B|*XiUuFRF8LH;RmKKJMd1@tw#hSfq;Q^1wgZzBta0>}Wd}sT
zpD?Dz3cpK2c)6>Wn82&SDJuP?r1P>`qYQ;K`<Fkv8D=XGgLGfEsfBW&+irt=cs_n*
zRPTv=$+HgxcYYCOU(Ae`l-_usL_sTJiA7-(s!hU^Gx2rpI2(nK`3D>?j?IdFVj@<X
zVh2eZofKVOxSup)c>cw2KGo4OGX$3|`84&ekYMgkL0!rNVv2--s--G58<b8pHJSRP
zs*numvshbk?+Bq}xIt=LoBXBp6B9I0wAT<*5|=69J@OhF)z~IYM}tgmTHMajrjPjg
zaL12RFLZtk5zF>-eYWS^&!NJ>czm`!2H$sDFR4Mgy*9t`zP$GHkjQ>>C$UNwABEs2
z7Az!*h#U7F(G7L|E&5-VaAt=p{A$160!)r5<je-TLInG1CnXT52d9Oa<L7I}=Y}DQ
z-nt)Y03>13T;m{=ax6RM6Evfhe3XWukRF|lF+W_d)*<+@Y^h%uRKZs9*o^S>6bHY8
zbdNJWi_S>AhkA!p`3i>^2Etaws?Tuj{<bq*Wtb9vi0lp&&n_f|oe?JN6aTe*?a&a=
zNU_;ifgznD(UpfQA1j2<9}8DobIPlJgVC=n3#1_*7rkY)C4Yqg2R4Yq$dsiNY)7M^
z4t^&aG;B-lqNx0wO_5uT-8P;lHuGCx7C1D0Cb<1mY#9#WZc+XRS^d7qCdG8EJ2~pf
zsBc3QJPQy)F|kao=yvmRb4<;wZ~Z2$=Q26phTs|SEa7!Y%uq|XRW%ir%lM2@@^O@A
zW`Oh-iBHD**yB-(^tnnZ+{2bKFU&|sy}08uS`u_=a;K5T@Mw>b!j~@eX_kY1Y2}0M
z({-T?$cTwv-gtzcBz!CBBQHRCW!w1?F1xw1#^EEpN6&In_HO)ERcB0f;>geZ;iG=$
zL~^E=*ibqbe%;bM(4-eKT(+AUxy2Iwodck>=b1^%kKUP(*GCp<Eq@qYlEi+`g34IL
zgA4XW!ahlywf>2t@H3HyGZ_V%0Z$Y3{HLYtGZi?IPZ_?&au_d3E^yx(6Z1{P3uBa;
z>(zbF<^Lp0q)YDcc4QsPUI`q4NvJJO<neqvvzRR{XyNPY5{?1JK=WsI4u#LB6tLa_
zOHzC)OG>!Pud%Wa=oByg@5Loo)=`2=oe0-wC=JB)25B$kdC?{16_o_2eu*jdz8IHR
z+)=eNZBb_23e3SL!TKzLJGuYrUi3BOD`(U5Wcpq_+mj?WfmYZ~BOVg7-_lE*=o^#)
z7BmaSXdkKkthg0?{q+l3oFq)AQYCT@Kz!!N{!A!%Sb&^RykTS)EvQ#AY+Q$bT4z!A
zM}bCEx$YPOc@-u@n@>4Mn^@W`>J&adonN4J?~BYKzBocJ?*$%JDrA_qqTt1+OZxp+
zjqu#W1XvNvTc2nHWY^R6<p9jCOWqyK{2I-6B1*vo9len>S-6-My!K$?ey^W+z{HDj
z%DVCPDgOdmflNPnG@)rHWH>h>G-v#^ifuHSfJvAQ-7q1R?U+qa?NKYTRkN-;5i8O(
zs^5j_4xA`$48}a-887ysfIy1fEC+=E1k>4!a1T)&H5^(LPAW0|E5VE<T|BR&wrcnr
z$D2$Yaj}Rl;x5wuxB|-$k}4Bn07QL;kXc-6;M&t)DbJ|Ny6kfmNg@f|O{F&;d~PU*
zPd>V2_jkuSl;RyBZAN(7VVobIoAy3oLF7`Hr{9&g?yhcXEdNNt{#F>6I&q6z>JJ+4
zO2phAJuDsvKVK3+Q4O!yy_6u1f5C?0Y-Jez!BLUoOvfveDIThMdpG2Z0O($HgV9zQ
zea}|}HYJGMztFUDS&Sd)S;W{n-3ywJe+vseXSf-NHI&5hk;36Jrz8@qBc#z7%G^*-
zJT<S^FPpBP`J0(bM~(KO4!=)NK637W+X?Q=P`vsv`a5Q>Hi~!;28}<m;~$3%uz$Bu
z^zc<o7sD4?BSrprKl!_vA_~B00bAH04ajw&ht%L>7z&71%v80c(ef*rHNPT<vm<BC
zIMAlR7XNr2POcq^)U#(hpiY9rEsKg1ZXH-i<$BS7i4uUs^4ye1lIff)vjIK)t>5C3
z>`zBgkFHk%xiR4_bph1~0Z5#()Cp9tXEe}MtdV9Ih*l$X0yjU{En)(P_Ffe7!c>0V
z)l9F`_hYav`%skF+NP1xH2Hie(9dj1(Kq3Q-evhU=qvqv7PW-@mM^gxYw0wf;*FPa
z6Yj2TD`kN)lk>=u_rBZG9-Pd$>l)X9yEoljNeAznw(FyCcbp5^hNc)bH&8S)4@DT>
zg*&_q&dbRUvTqMF^Z42T8@xD#lyWMS@5`@Jf{6e=rQe6}n(Mo@stw=hqGoNzN>^Ou
zNyiWE)1k$UHz2ty{q)mL;jQ}%;qG`oQu>Om3Gb1rQhOd*A_5?MLp<K8^CWZKfAHa>
z?^$YJO#X2C!}&uhukdy8<kX7eSr_>VK99Nm?(`9k&&c)@{y_RFJ_ZGkMT)`R?=a-R
z(0nsQ4$nVkB$MX!NzlK%)1z2JOa4{Hq4J3q$df^LY;ozzMJAi1b5-t{A4$@JRk44B
zS(7fOm-Q_yn9!Sd;3V@H@@_V~5ff-Me%5|DIwyaYE=}yZ+^hRn@b!gW@gId|A8`Jb
zgSGyOz?<1!Lpb0Y{5x+lamz_)R;Ue1cT$Ez62jh+?<;ho2H~afc#U3m&YlW?Rsp-L
z2CJ%TE~WYu@{~UWek9|{qX<FGra8n>`B8Y%W`|US@K~}#bkAav9hAOM4!<Y=gXvi^
zsjSZdzlm=ha)?7r6WHOOX#9-9x*8O3yz=7*mV8{_DAo%gyp*gVQnHQGMTaV)?!uMe
z)}G6dpS)Xt6Y;0_@bGmjS5{a~t(Ni1TI*rEtU_yogdw9Le?GEqpH!Q4a`ZGJx46?P
z^#t+TsO{U=T}gYZ*wf3g75c`pJ5}bw?`+;HUV-$8hH|$;kb*v5Lu}e7&J$JJ{pmW9
zF@Cq|6>fq7BdgAS<XLE2!S3=Nd3rvRyF;id7H7t4EKg5$WrM#pU1$uMRT=b(olC+r
zMwcE>QTH9Q0htWYlxfH=pvod;F&jHuEK{Uyv*6a5WRyIys6BA}ySJbp)gU}sO}M_Y
zSHMSCNxS|WeRi^8u;H(BaQ>MoD*x4XcaFAX>is2fdBuLXbv?J$cddbLcG6Jf`{lyY
zi4Xo+=Sb9Wj(piLOH1jNZi!K+BV*%j{+G}ihF{hA3YGqMTU*MhK*R5k;qt&`8ms?g
zyx_;@!;X48l#6P$@mB6vJ2bx$>-8#R7ESg>D@*P#_L2%X`DN}nl5S7mus)ygchz)}
z&o$oz>aeZ)HG&owC9IU~T17tibCX;Pob*=IC%)l!GCh2_tuDT?v0M7noTjne9PL(j
zZX0+*y0;Dy>8PvF*IIv+vW^`xNJ#x5WronK!TtJp`)c#N?~0@FHudv<^_RM@W9Ig+
zE>x#1qsRz9kj!;Wik3ys)=wq%{=BHMu9*A|d;rb{os{pdF2;Me0yz0k`a^b{oHTAf
zT4UEcdTE$o5Tsis6eeKMRXV45*bw`rsw&Kj81EDa<b?3mtPXy;&SPRUa!R80-Ik$U
zA^7vlVSLr(95}bm-n!Z=*qV1@EiU2v_Qr5i+l+vvBx%9L;qs_qqhMy+Eml?cl#U=D
zc%IeWyjXb%-AhpSEE@9W@34tGSs%V1p@@<0J?Eb;Rd=Yg@~JPh{am9z8N^QJ;dr(Z
z(s+B^AHsR?Q11F(hW5hE(sE4BsF77p$>0v=*bw*p#(nSE{dva0qE@5(#a^dM`g75e
z#{Ot0{yNv&WjO9J2F=}h>pZ7+lsR50<U7xWadq<xorO<Uu4M7<xlZ*qmjOq`Ol6I%
zACw*b3SNw~R$OCm^V&^zF1|h@a8qWU-o)|5+;F{6N3R~L=0I2T{X1EE|3U<f+%85p
z4iA-_joeDl6a@+eam2vmv;iW~{s(7F$4D=;makuw@p?hKpd7oeeI3u)`UFCVn7?Lf
zePG)iT--RzE$%DaE@o*SauWa#&o9Vf@Bi5^Y25$D<LP=fc+gON)mt{0O{cy0hLzvr
z#B(F(;PS%NcJ$)r+U0l@a#z0q^K`8(5YS$nnj9ipoOi!HnBSkBUjbxPuTt;_9j6Cb
zWV#3I53BKdH8lGxc15fv(j#YfX^6?*-59&wu@zMYGY`^C9wo@sfCIM3lt%&SgY3>!
zL~akuwCr71J7l4o)02VztX;#eWtu<DJTZ7hJsmQNE^MSh>&~+0incclqvvJw38PlU
z2SH%^<m2Q5)^=7RRv=w|JrIM!3QhGRne>^y$Z);V@{h?kujUD&SS<$90%R*WELJgd
z?n;ZToSM(7X{bLM7YbAShS5&iO(j-D@qfK)zOJV2{2g+|Fi>>t6toi874WWBOYl6J
zHFcqjO!`(|Q?4<7`f|iWdrVvST(d_#>D8ZB2G7E^SxtK#jRAm#m9X}eNHOPi9lH^Y
z2iv&JsqopVSHtFP^iMOn@>?{e->*p9FmrB7v+1>8{N1>$rEd&G`7v<mVm2aiWt~s9
zBaZH1GLV%M8&6adIOM&QJ7C4E7D6U@i7RzL^=zSbI&NRGe>{it*!POl^c&0Kj=ghm
z%o~!1HPRKjCXOvYn2&W0SZcI%UHxkxYBx~fvWz1Z@;WZZYGI$xt<>!XCek6r+Q7bk
zvOK<TH7ZKKqId6nb$z-VQOjR9<xyQ$#*pt4p62%}zbj;Of3#kRm&}Y1m36c!K^-sG
zgN0=^QrYl~{9VTTZusglLjr99<Q$Ya#Qu;?O9Vl*vDi7lEL_bw3XUQ*1t6AWPKVsl
zP1rjJNu#oIEiT&rH7PC@Ho{;f)#us3XRcd1w0He0u4iGpx8JH@&ldS=F2zS%=aTkt
zLi!XzEuN=aH}iS!0mFoyo5#^j_vx9al{Wv~=p<R{jQtq5qV1^{QKg3d8qdFEmg0ZO
zEDjoLe87Xwl7x%tqmKC?v(O1-*d4zOVYq~BtaD@TTNF2ZX4_w16|prNpf_+>9=l}e
zIv$y+E<5Ve3@G_J*S|$<>$LAoM!7lM%t|1{<#adS5tC-uG8T8nG|anECW7rupcrU0
zfN8XA^#gC|1vY<jW8gll(c6E{gSoubXs~hi4p?eML;75q{GhY&{~w)Y@PBodlLwt;
z<G*#5-SR=Iz5lDT6a%goE7xbzjsMYEIR4RD7GZ1SnK3t~2}Yxil<8VUP{ww+1EyhJ
zmD20&ga@64KN0vJp~bl+jZi73=T2rVVB6M=V&pi5(HR;-P9QHLo+gwokj)0~xS#UM
zhZb)js_CN*MDhZPH=)L@h)$73ZnNg=5_pY<AYC8aEUNACQ#u~&EbryYS$$Q~G+Z0A
zUnmiBEwQMWSmfzkHEwb3`#dRBnJp$>TcQd^tbJMlo!VZIZ%e+G4n>N{X!_4|<shuX
zlqO8&n=k5>H@9d;{9c{xItFQxO^@z=HU6X|T6X-<GyJ75A9I-{B=P$g%7?NgpG~Or
zC%jRE%dF7I_Ms9xHm=0RB||UD%9C9c#NU-jHcD<qkjc_h+8Ga5)Fnbx&KHrU38T91
zW519GvJ;I^lo>vA7o$Gs`qINqvazFuCbUQMs?C6#L3z{@XXIlENa{PP@UX7*hzf<P
zaTVPZeOI$A-m~oV*i7M*Z8|rigt(05ujaJRtJOtXBlSgty^RmK0@Yi;M=R_xGWu~{
z4v5A$H+CcA=1bC%=_pzcPD#3<$uQGKsT&J?0ZhVgnc{re_L{_>372w*os~1aCn}~5
z@0+ChW$c50`}f@vD0Rfz)M9E`kVN>`SV};mwf5kdGqcRG+vYmICC6`-WX|1mQ>fA8
z2e@IXFWF~)&*iKjQj|b(Ipz3lQ%6V+I;=WUv+Y&)+H%T;?N^y@Lp0~)elKb|kd9k@
z1mLk5v#laS|A9c^r=6QP%`AI)x&~F1lw~#dsVGb=e$rJhT=QxuCkU=xt6!_Q4FPp!
zB5i4$y=d|<d&R7qn>;JM+jVXq{?BijoMVH8_nw-K(xg|<XfVpn)|q4b=If@0pD$UP
z3vDYH+I;W+vc@7BM-LZ-we3$-udlz$SPKX<ueZosXnjckpj@RsPd*<QEo6%mDBO|5
z=QOwGqnRWa1pjJLKK)!iLld$jzD?)j3U(j<T9#Dz9KN_T@0Sn#@xsIl&fp^zS||oO
z=X5etV3VL}uzYcvK52nCN!S5a$jupc&ky*mWcX_QF|$|;T-@v?1{In&gtp=Ffhd4N
zF@xN!I1OBJ_Gi853}pOlJ&J_4q^kxVn0Wg97-oj%aH<=Urk3kOq?T_)an2|V)2VQ+
zs6{KK=r@{pFEQR6)6mu0$u^Y7B0vA>&C{iC(Nf<%l-MsJ9slw<+LE#eWS}SP^i9W)
z$e&b`U9PX5u{5VgAtkz9$j!4!Sz`+TP>{H9ufC#76X}TdG9_*Qem)Hwp`mXZ*H6w8
zJsk1qcc0znq<fijM}N65zJ&j&uCYCBszzOvKOHD!+)U4kQ#&#HRX;kcKL(<YBJpba
zOeCEm!)XeWIT;cP<FDefr)m8CS0?3H#`9hKC#7~L*q@o=V+ZZ_%5Vh-ju)(sz)P$a
zy+~1aT)}F&Mywaq<Ld!q>OPcCwcnzz0&nOKyS_6KC#0l*ipv#5T!=xmmFv!6AM)$#
zF^d<%O(gpJ;Zq?=zP)dkLe1OqbVeom1}JJgb>ml+zd;OAX`j5D)eFBCrg1Y}hQF^6
zAjMt{Y?;9?^Qu%$xAZn!wo}a$Is({h^K=v`SY*!@w&2Y5dtf4x;|!=8Gl9RwG;ERt
zxnJQC&6vM@MURJ_S)NJP1pDn!N@DE(3(6~G!yy%(@8FzCjD@wA7OJY@Er5?5m|YN_
z(b7RHy^eGdBq(8y=*x(J{{GYWoDB*F(MXTT%$V9~*bmnBaQH=Wt9%AqOJ2Z+o0g!#
zH(#3M9n)kWUrY?&7~)j;OZ921v0<c!4a9vE@|V4@YhE1UyHfseIBtD2mp+u_^D4L-
zQcBRlhgT`ANwP!yMQ4txY0%-6Zz85bk3R7fTi}_<Y{{UGHCr`-a-d-&P|}@@V!r$(
zUpx$3Ax!EDO6!lY|0>{4)rlWC@1E$6GdLIY7bLo-A4})e)<Y(tQqhz#PjnZys+HLk
zt`v;ofyb3-4RE+$RbWdfizhQB*EVp02r}7VJ8Z419dZ%h(O{E)g_D~toH1Zhi<B7x
z3)r^fq>c2qB;H^lRv6<V=1Z`%<qT|IH5uZg1WJJT8U(2kxWujM)-iyi#!Q7jBZz!b
z5?@<SolcTITAQ_NqE&W9)})d6Z&Fk+mtvXQbTzqrkt;rCo)Q}F@90n}7t};42t5h@
zfdeSQ%oH@t)}h}*XR_4BCGbKx&Jk0E6`nD7DV1jP7{3ks=;N&vpc5v_?uRUuruSxy
zF{r;Mr5JA+^X+((&rb+o@JihqRW69uTB^dNF|o}+^jTDom~e&CZxa6=6^7dQEhdx~
zN_h|>2ir6@HoMRxUle#<lwB-&$`(uBxWE@(2p%tQ6G{18L-84&%lP@aCEAeFZ~!<T
z*+8Z7CeMA};p7{YzNGCidl)g*<d^qo3s%(<1ka_qv(c9sf4uYr;+?H1qXh58eMOEL
zg46zZ_UibPyyEH#ULW9S7<v==CCyQu@CWtO#NXoesw|F7wjr+H#T@krB!(Z;iTRef
zf7P91=Pz>$L!;8nOcOi4;dzxvVHb8tIW+%%#hw9gePl|xDYlGV+8m+Lfi?0+kA;Rb
zmFbti;~3e)ztsk6P5_dk{ShZ06&*c-Te~rT0nfD5^S-&6RZ3T$Dr}Rut=Z^y(@$r7
z>5@{n8MBbQVRTu8>XH_b(A8X``e?uLL^n1U9&?&1oHUWI46SU_DI+gjAb06GY(xT4
z;9VS2sv0{msP1iIEnX0wpZQPODz+rr_sFMp<QQ<TP3@J)n*qFHt}2q3fyQ>@q%Uk{
z&Te(6L;FH*me`k#H|XaZL)XHB6n+NbHe|huVh@NKJ!nl{41hJxvvuzBldix1nSx>W
zA&U3gTJEcPaPJ@M+RhbK=^0F#a3ffTbaJR|)%omaSP{I_?Aw4lLw@e*0b4azq3;>0
z81I)m(LQtcrmX;^RQY>v1KtPef3Su2!s^<}zOPxIJ5%XHda>=JF@;hezEn%@yqi^(
zF1RLN7_hA}nE%`THpijZ)H=QYZNL~&fAnWq7@P=_4?GcV#LbSNK@H~b;luj(k%)^v
zce2oqm6y4BWI5Y%il1OeIt(t8o%TZT-#-V$UTD&A?}`F+s33VL?bnNxF)wL6C%6|R
zMr7%Fg_BcA$EY?NwzBn;Y!-AdIp>X+e_p}#WIR+b(S}$s^j^Rj9uGOV(IX7cPBOzq
z6Qx}SPFBW7G(&ul&;$EjExl{j%@$8+kC}6Iu`o04(9he^k5$nb&UDClPz(zvGZ=MG
zsq3hdmcIhza@5a}&akz<5?ySw-t;cC`bJ||rc)$mN9pt4T5{OS`9OP^ey~JZasfZT
zd@&O*j?W^CV(!+Y`@0fvAy9m(_MRjr@{8q`6l%YDM>SIqO~i5<P6Pu%bL;ITRZkXQ
z2rXVp#fC;^Kt<AoxN^epQwfAZ1k-`9eAtZ|@5||c9~7e&;mFv8QU_V0Spwgk2>Db&
z3{erp`IMnAyV=^l%2=-S;WtytojuzL)Whe+q%)hrZkJD;Lu|2Rctxu9+tc)8=?`|d
zySg-$#Qs=?U=inwKOyMtK3iw5sO_dUk}Q>IAG4T}!EXcQVKgdu8Ks@^x(Hz@(uDEt
za}2HNz%&@qx$S3cBL@<>%NT{<@YttBp(OoC7C!ajyaLpT8{>w|VP)k(mKr#=2yAh@
zBcygjUxUx=bs=z?4lJM1zs+|98nkcyeT_ua5gObpogPuc=Mbvzm1Ze@HXPCG*i4eA
zi-$x<$#(UQzs;E)&y75s3U{Hxc!(=lmYyUN0Vu5oIm?S@pQ@{iz9jBgVo(rzKdTqp
zo<bVlj9dAg;Wn^>9K%Nq>XC6gTWh4FI)|A8?R_rqo9KrUMg*Zs2bcQlRWYT8l=9a5
zEK_z-eHMIO4HBYCsS@Wcu`=-s+&a4E+P<(EyFcJEW&f<UusM7Fp|XCq((}&q;NGr@
z2*|m;u)=WVUOQV4<7(%S{n}Y8$q3$UaHXR#Ulez~JCI4eKf4^g|9dJkx-qEjZtnf#
zf{^dZ;A*kIKriQjMAB#hhGn_t@sV&1+l^KHz5QF#AzMavg*`(ljC$D4bZK>>_4TJ~
zqBCdnlXa!}nd<;L;O^$6Li0GmbyoZ87-*T1uW%Ty#8oG_oS3O_Uztf|Eq#D#-DzNk
zZqn+0<7$4qmW&Qa>i4|zGU6~3s@LirB{oy%##}D?#&i^*ZbG_z-rm3+hUUa|PSEc<
zd2SlGJ(t%U;%;`NcK3%CRj90XXuV2mBnA{_IPg8*gmQi*N|m!EOe*AXqL0T2h)gNV
zaqLKBE!Nhohq<<z-FE!GZPd7^V`%&;+>g1S)#F+V()mXoB3a!jQ0?BE-)lvkeUOJ7
zCcb~TzieF4zFkRGF3CPvxb5O(U94+dIa9nc;+?ZQoi`n~*#0^FF1gO+x5$}3J%3qc
z_Gwy?wm$d#%`gqX!ZSFnxDCvij2cAVc4uB->ou5bOTEzfw6{wmHJY?DIeKlS^&o2)
zxWWT(8?>tGjSaQxbKMrN^69WI=*wERWmpu2l4<^4bEIV7I`)>}3m%SgX>s$tz2%I1
zR&{H8gBk5yaC|Q1McCrmZ{L}R>&EBmdUo5}cy~J9>vYgq<>_?0nFBa2R{x!j+1ffg
zU30k|ccgYoy=R@g--&3b_dMCYog~u_nuO1|d-rhxxVYNuMcob5Wc{PP`Hw(^F~=Kp
zr@K%Td~_<i@Cnm{EdB=NU?htt)JgN(tp}(kzmuRQ;Z5m$ge+Nntwi?AW6M#<pBwco
zOY<XTpP?@U71)zR$Zd^tw`BD1GbLjLJaJ<hs2~0eis_+l0c`O^U-@60;Z^$~h%B%Q
z_rdv5Sbw$Ac|F+T-0U^<eiN*svwpO}t(sMwU%b7I7{JYg{h->C8E#bHZwwYTTpylZ
zIF-l%t~c+$^H}X;lxS&leyHELmNYp_IT-B8?Q~CNPgEw5`gMF#cUKNOne870IQwVU
z_dVjZ%9d?#0m%a2<80=r3m1afOa-x0g6|O`6B`qOg>IwkhVv)BL{@9I#`<FO1P5cK
z72A~xXEF|X`wNwy+>$Jg1riE*!(X3O4_k!Rqd8@xCekKN#4r7w4Xf^`n&2jw>zEPH
z)3p3PNu$!%s!pee1`P2i01aES{(3c-7Gks?AX9vSKq%@jp>o~jQDJwUl}f4ht62Kr
z=wQ-9GU3wl4%LFuvAoj1&K*|juEQ#tGrLh}hGpbx)<?!4)b(dp*=I36xJ|CNDZdHT
z!m?;Q`uf2t$#<CH_e`5!kIo+-<NeN~&;^c-bQ>ZSRDpmr2!SxFCrSV2H`wbvcBBo-
zH~9|YP~c+4`JZFmaOVCc%`FJSa}g_q^%fXZ^4A<~%k^X*RyE3(I)_73zLi-Om;5aS
zU3GjlZFU<*-Wp0LjZ?rGuGbrIo%X&6=_dQDUBI2_#69UB`Aau=daQ<iW%o?vWJ$#C
z0y{&dfWV3=%;hNiQe=LlXrg#dsWG2mk#iBKH%?*y#Hj!(ANW8l3UR>No3Mng6Lw-1
zp<<RbwIRHaH_7<jjFA!{%!JtXo_D8acL__Fn^Onb`sJ|Yi^4LrAX$0qUx(x`%)jnv
z`=kk!M*}^Qv7$MKf1Y3U`y9I%g$pzF=u&XJJ{3`Q_Oj(y=G*x$9<Zmh|J~pHbAcNl
z;Ni6q(?U!)EbJV|2@Rme4i;OfDgQk=N=LW7Ruoe9pgz^{UYP8VM{W+IqrCM4rO9*_
z522@Nl%#5I>)}yTkr2oE8psuyXpZAsB<Yh+Qmous_B-Zmw<rz!!V{1C8hmp2Y$MrW
z_*ze+AnAKPdi5ywkEz@@)CKurrA1%t0OwS_Q>K~~%X!IByt3ay@uKMtd$G>qUT-_G
z!eweJHs)uV`BHPIDYmJGDOHe_z9t)_ndnol1R+;UfL$xE($!nPO<~S(M)UUAnQyV+
zdRRLaE%ogOmvy`tK%1sGr5djBq&R;^BlLIchW)m-cc2E&K$z5jDMHxo)b3a^50E|X
z-W=bQ`FJ2ibv+kFDrISh?FU}cqTdO$u|3;5UO98lDbbM@RFotpXKc)sowumJ!e_RM
zH3o0zbeircAyBOanl1+_nhPhPMyy{>e$^kp7^hYzBym9l39Y~x^mN{7FN=(L5k0Ro
zv|Ku;3Cy-(i)Q?cvf8HG21hL(0>mv5>T^89L;B`2fs*?w0CVE12HFERD}EVP$d;k@
z0xy<Qo^>Uki&b8D=2}{fCkdPCFC|a<Z%c|67GW5scwQ#fT6DTFNk2`?e2g=p#C~SC
znzPq2!e-tCxgiM-ZmQz5`LHGpI_7r?8j67y6LQzvB9^|t2UHOlSEJ2yxq(<Ef~M($
zug@js-Ow0%MUpaC`_@r5++mAmuD%-H)p`a;)jrP}P-ksnayhSje|9RuF56|TwnR`r
zSB;d=4eQ8f*A)6G(aN(9CNj#LhMTyglIRFala!xPhL0f9Txn@Yc;U%d<)Zm9uIibV
z^jvKUR(hCB4<&eB{{ZDC-~!jX+-A|e!{bcUfiHcjE+#hg6y_0CDE7Y9K)xApH&-&!
zl})7QL&&kB88K|>^u(Fh_I(_0_y4L1B5O{e`ld@9;fb`1hi5uowO5Z5Zj+h{s;=p?
zh6m~OuwTal{kb3!N9w@}`OFzth9cZA*B@2o$j{RIM#Q~SkTlVi$^mSU$9CSb1}aDN
zGO^96`2qnazHHWMa0ocul#cC6f_0#ZQbc<;yUns9q18LSpmauG%Sk;WWgYLcG)$mM
zv$CdqgtB3ykXE@gtDM9ER&&p_TX!CyD@N$mWoKzGV)}Y_xOtR8@5*x*+giACWl7vB
zC_^BLB`|X+K`bD63IJ)nyZ1HAiu1s-@hX{UAh78pIgIhs$WdXm|GK39o_0+5@Qu2w
zZa57c6T`%6=h!}4WZ<T?;3a(@^rFDT*-Nar6@0q}U;fVKl$>OO05508o1$V`#8eX@
zKjekl{)EpPTaPWp^fxHqwYVv4X-DkMv1E^gvjoEegR*Ia8Gw&@S;e&`&y@Xs*}R9J
zaMor;nEd(vRjZ(0aYK`Z5aP^KA{a-gfI{x-8+P+VERjf|nIuIY9^O$qAr(!WtBEbB
z+avpjxFe~G2zk7MawrZh*O#*TauxEN^Os}4V=aw%maNhS{A%p>w5vN1cudTq#l+ko
z#xi^KBZ+E6>1M!KCQWR4VxT}$vYFPQ7X?i_Wn3FPn?nqmMAQq5@OdODLvN40QAC?6
z*>5*%KT=+khs6KbFq)#?DhPEAmO7OTCRxat5)RD^SvH!A{bn$1t{X!x&PeJHd1fGq
ztM)#+MS5C_BE(s>R~D(BYkN5j*ZyUv{GtPb*yjrT(Cj<l*_uWV-jY~W*&9o>_GGMI
zQg|r#Uen(+Dk~4OmkpD#5MSuMUhM7mm^f|{lJ017PF~LZ^JA+%v{g%shl_OQ+m{r`
zC#D1Rh0jc#I>8ERHLg8626)ZiqCO%Jd|UX=69qNLD)B|ItvU01mfr#sn+fP1TZ7w1
z{enzQWx(J(34D2<(vIbLnA4It@LfhRCxH&yAQ)aT?c_#hi~#xN_~y^YdFKQ>duk^U
zwLzC3KWF`5M1hhPyn&Kip9IsIKdBRH+M~@R0?8d^IA#ZOR(T=m3GWa;)(YCQj1?j1
zRrnG)AoR-BT1qtp(gu`flxubNSs1RrS#gou^-lz->VZYE18tvEN6;*7eRBulmeh<>
zt6RwF=XnHM3D+;MIqN^zTW9te^R)GZ?gN}7-luRIX{bHMKgF_*^l7@<z1y6@N<4ht
zDx4=Pp9}wDG*&XO5??<Pog-pudUt+0#Q%m#$5eEYPJzbRUm{ToKQPJhBrQ4yV#BAn
zy6$O`Jp@GBm}EOYb2G_?pH2G3uJF3g#U|IFo<~tn`Kl}}y0-cqF(YTIK3L)ASQxn%
zxBZ@9GS*7e_u<{xH>@Lk86cXV5MO$VW#_(9c@ZODIQ$E^4;)iuHi{7iRI`2KLt?I0
z!-BtSt<P*RZ^kt8yC0v#1l(8M`?{X4wp-Bn0@C|X%Xx+S?>tz^i?;c3>e9*#lA;62
zc>L;GP-$zokFMcRR<VMDGgYXBl|Ao*V?y&xA~trc>ONgW?#S5UfwZUlsqt>Y>PMC3
zNv(u%-nizaZNu}iNddky$G54O=O+{m;PPj)b?bZG<Ly8(iwrA2W7%q%oVGP5#Br<(
zI1zP$a-P{yrPRQ1fa-*7Po$1FtC{mrQzh@5B5uDFg6!5YFpu*(^)C}@(P*Rv$TrK8
zIJS116u6CYj6Q@vx9lw<{?nEtP&M_<45Pbya0A~eMXMghyXcw?y?ZZl3%{KbafilY
z2s4@pnGOkh;p^TD+SZ{W3;0Z3LMA^AxPuL?P(*7S-`%b>7RvzFx2IvQPM1;ITHmz~
zoG~yfyrA+~0u4=EyNTu&MuWeza~sgtnB0id0GFF3^M=dsLU*LOGRlnA*;08Z(-(@j
zp_1>+wTaZ|Wwg+v<u{a#wy(|#PDb6frI+aOPrPSiw2arVvORVvtV~itRD}(I(+|Ry
z@xtlKQu$lQ*3>M&^0*@Eze!m@fq)~G*rhCuL(io4DpJSvubSfF{Hc?H*3pX+&%4u`
z@wKzjg5Cn|RF*}PizQ*+I{Cqgy@XkTl8d!_>uc==S2u^d^{X`t9`6WSom>4BjYjTT
zX9q8`WW%23>572Q%2VD8qnp~o$(mOErN5Xi@UDNZ?vxc|RCab3hZjvNXJojz>$swd
zcHbhMeVP8!LCR38iedO?&{PxFOiEL2iwb>yWYmBEEM(VsH-;ZX3=f9Epk-R1dh3^5
zATS#nJ6jB8x)dWj8;qRoUoUb?a+H5QKtW&#$5R_;%y0oNBPZ9BZy>OzH$cH`F%#*q
z7hD|B7@Q0mfC~(IbdC)Q<6wVe1A)1?ppR^9P%bF=u?@xnejJkx#>okJ<Oc%7*#5Z>
z_&?>bJ;Y^u;>ZPkd>G6Qf<7)5423}-7XpStxE?16<^Z$(U%!X=Y-|tNfFD-}g0OLM
zK28P#{r~afggg!cg1|T)qX2qn-V;A47~p(TF_azl)W!*VjKTv3(8t37VrOG}j0}h!
z#QvlZb}00JS{x&qg${5%xtslA5FZ!M4r6~3j|1{#AVC})>`(6F;DkNNixb55KQVy+
zop2E6{|~~PT(Bp}JdEv=j&VV_o&d+i347`XgZx)A|C#|1jO{5(4}SmU_>UjQlM#4W
zS<vHlgV`XU$H{=%pidSH%*M$HJZ=vd^bqM0{NRW2d0Yt?1c5$1%l-rfFbMW!mciiv
zpauMgV8CE@w#Pkwu)&_%xc*=7{<n2t2pII3Hoy>e$P;8BoLv8X-vj<&2n_OgXuwdm
zC%t=EEXZSIATa2Y&3b6ulXgPcI5?jm$ic}D0Ukrm0RuhmJD3y9{^T$xJICX}1aope
zAzd&hjQt5RTp+f87YzKH-@y-r`UI|rLY^eU^^^`C2%6*H7yjoy*wfX2pgP#2WB@iW
zjN{3}1`K<E<`FUwHW0^S6d-KSCmaT0<9Na!5H{G8ErEa@noRx=2mkLNK|l}8@HjAt
zi~TVQ5b#4oANxHF(*LFRhdpA0fZ3nS0_0!S{rfEN{~oW9hn4-0@jo`mlU01+yeF#&
zfpBsB_f$UwhJv3w1|U#&_Qx#-9=7<&>Oj~Z806m>{yQ0VF#F?ygRnz6p9E&-_^*fD
zKPWtyo}kJB`qxnZga3ax0P^q<euBcoMm#yp`9OP*6Xbm0vnMu=C))&hVCDavBgp^!
kx`A*(o}$17<$Qc06vWO2iCHfYL*ax%P-tjGWImw$e+|vo>;M1&

diff --git a/demo_files/output.txt b/demo_files/output.txt
index 837b7ad..d0b1ca3 100644
--- a/demo_files/output.txt
+++ b/demo_files/output.txt
@@ -1,6 +1,9 @@
-ECOFF: 0.29
-99th percentile: 0.29
-97.5th percentile: 0.22
-95th percentile: 0.17
-μ₁: 0.049303392476164935, σ₁: 2.1439340245815863
-μ₂: 14.283983662174911, σ₂: 1.354070842559145
+  ECOFF: 0.3963
+  log scale: -1.3354
+
+  Component 1:
+    μ = 0.0835
+    σ (folds) = 1.9534
+  Component 2:
+    μ = 3.8329
+    σ (folds) = 1.4984
diff --git a/gui.py b/gui.py
index 07b2547..92791a5 100644
--- a/gui.py
+++ b/gui.py
@@ -55,6 +55,16 @@ def __init__(self, root):
         self.percentile_entry.insert(0, "99")
         self.percentile_entry.pack(side="left", padx=10)
 
+        self.verbose_var = tk.BooleanVar(value=False)
+        verbose_frame = tk.Frame(root)
+        verbose_frame.pack(pady=5, anchor="w")
+        tk.Checkbutton(
+            verbose_frame,
+            text="model diagnostics",
+            variable=self.verbose_var,
+        ).pack(side="left")
+
+
         # OUTPUT FILE
         output_frame = tk.Frame(root)
         output_frame.pack(pady=5, anchor="w")
@@ -125,6 +135,7 @@ def run_ecoff(self):
         tails = int(self.tails_entry.get()) if self.tails_entry.get().strip() else None
         percentile = float(self.percentile_entry.get())
         outfile = self.output_entry.get()
+        verbose = self.verbose_var.get()
 
         if not input_file:
             messagebox.showerror("Error", "You must select an input MIC data file.")
@@ -163,16 +174,15 @@ def run_ecoff(self):
             text = "ECOFF RESULTS\n=====================================\n\n"
 
             global_report = GenerateReport.from_fitter(global_fitter, global_result)
+            
             if len(individual_results) > 1:
-                
-                text += global_report.to_text("GLOBAL FIT")
+                text += global_report.to_text("GLOBAL FIT", verbose=verbose)
                 text += "\nINDIVIDUAL FITS:\n-------------------------------------\n"
 
-
             # Individual fits
             for name, (fitter, result) in individual_results.items():
                 rep = GenerateReport.from_fitter(fitter, result)
-                text += rep.to_text(label=name)
+                text += rep.to_text(label=name, verbose=verbose)
 
 
             if outfile:
diff --git a/src/ecoff_fitter/cli.py b/src/ecoff_fitter/cli.py
index 11dd165..ef7cd7a 100644
--- a/src/ecoff_fitter/cli.py
+++ b/src/ecoff_fitter/cli.py
@@ -117,13 +117,13 @@ def main(argv: Optional[List[str]] = None) -> None:
         global_report = GenerateReport.from_fitter(global_fitter, global_result)
 
         if len(individual_results) > 1:
-            text += global_report.to_text("GLOBAL FIT")
+            text += global_report.to_text("GLOBAL FIT", verbose=args.verbose)
             text += "\nINDIVIDUAL FITS:\n-------------------------------------\n"
 
         # Individual fits
         for name, (fitter, result) in individual_results.items():
             rep = GenerateReport.from_fitter(fitter, result)
-            text += rep.to_text(label=name)       
+            text += rep.to_text(label=name, verbose=args.verbose)       
         
         if args.outfile:
             validate_output_path(args.outfile)
diff --git a/src/ecoff_fitter/core.py b/src/ecoff_fitter/core.py
index 0b30c3f..33406db 100644
--- a/src/ecoff_fitter/core.py
+++ b/src/ecoff_fitter/core.py
@@ -326,3 +326,29 @@ def compute_ecoff(self, percentile: float) -> Tuple[Any, ...]:
             mus_sigmas.extend([mus[k], sigmas[k]])
 
         return (ecoff, z_percentile, *mus_sigmas)
+
+    def model_summary(self) -> dict[str, Any]:
+        K = self.distributions
+        n = self.obj_df.observations.sum()
+        wt_idx = int(np.argmin(self.mus_))
+
+        summary = {
+            "model_family": "interval-censored log-normal",
+            "model_type": "mixture" if K > 1 else "single",
+            "components": K,
+            "wild_type_component": wt_idx + 1,
+            "dilution_factor": self.dilution_factor,
+            "boundary_support": self.boundary_support,
+            "n_observations": n,
+            "log_likelihood": self.loglike_,
+            "converged": self.converged_,
+            "n_iter": self.n_iter_,
+            "pis": self.pis_,
+        }
+
+        # Information criteria (if identifiable)
+        k_params = 2 * K + (K - 1)
+        summary["aic"] = 2 * k_params - 2 * self.loglike_
+        summary["bic"] = np.log(n) * k_params - 2 * self.loglike_
+
+        return summary
diff --git a/src/ecoff_fitter/report.py b/src/ecoff_fitter/report.py
index 08515b4..e24bc75 100644
--- a/src/ecoff_fitter/report.py
+++ b/src/ecoff_fitter/report.py
@@ -108,10 +108,10 @@ def to_text(self, label: str | None = None, verbose: bool = False) -> str:
             lines.append(sigma_line)
 
         # Verbose model details
-        if verbose and self.model is not None:
+        if verbose:
             lines.append("")
             lines.append("--- Model details ---")
-            lines.append(str(self.model))
+            lines.extend(self._format_model_summary())
 
         return "\n".join(lines) + "\n"
 
@@ -175,6 +175,28 @@ def _make_pdf(self, title: Optional[str] = None) -> Figure:
 
         fig.tight_layout(rect=(0, 0, 1, 0.95))
         return fig
+    
+    def _format_model_summary(self) -> list[str]:
+        """
+        Format a structured model summary from the fitter, if available.
+        """
+        if not hasattr(self.fitter, "model_summary"):
+            return ["  (model summary unavailable)"]
+
+        summary = self.fitter.model_summary()
+        lines: list[str] = []
+
+        for key, value in summary.items():
+            if value is None:
+                continue
+
+            # Format arrays nicely
+            if isinstance(value, np.ndarray):
+                value = ", ".join(f"{v:.4f}" for v in value)
+
+            lines.append(f"  {key.replace('_', ' ')}: {value}")
+
+        return lines
 
 
 class CombinedReport:
@@ -212,13 +234,13 @@ def write_out(self) -> None:
             lines.append(f"\n===== INDIVIDUAL FIT: {name} =====")
             lines.append(report.to_text(label=name))
 
-        # Join and write file
         text = "\n".join(lines)
 
         with open(self.outfile, "w") as f:
             f.write(text)
 
         print(f"\nCombined text report saved to: {self.outfile}")
+        
 
     def save_pdf(self) -> None:
         from matplotlib.backends.backend_pdf import PdfPages
@@ -239,3 +261,5 @@ def save_pdf(self) -> None:
         print(f"Combined PDF saved to {self.outfile}")
 
         print(f"Combined PDF saved to {self.outfile}")
+
+
diff --git a/tests/test_ecoff_fitter.py b/tests/test_ecoff_fitter.py
index 192dc44..dee18e7 100644
--- a/tests/test_ecoff_fitter.py
+++ b/tests/test_ecoff_fitter.py
@@ -236,3 +236,39 @@ def fit(self, **kwargs):
     assert fitter.converged_
     assert len(fitter.mus_) == 2
     assert np.isfinite(fitter.mus_[0])
+
+def test_model_summary_basic(simple_data):
+    """
+    model_summary() should return a structured, self-consistent summary
+    after fitting.
+    """
+    fitter = ECOFFitter(simple_data, distributions=2)
+    fitter.fit()
+
+    summary = fitter.model_summary()
+
+    # basic structure
+    assert isinstance(summary, dict)
+
+    # required keys
+    required_keys = {
+        "model_family",
+        "model_type",
+        "components",
+        "wild_type_component",
+        "n_observations",
+        "log_likelihood",
+        "converged",
+        "pis",
+        "aic",
+        "bic",
+    }
+    assert required_keys.issubset(summary.keys())
+
+    # internal consistency
+    assert summary["components"] == fitter.distributions
+    assert summary["n_observations"] == fitter.obj_df.observations.sum()
+    assert np.isclose(np.sum(summary["pis"]), 1.0)
+    assert summary["converged"] is True
+    assert summary["bic"] >= summary["aic"]
+
diff --git a/tests/test_report.py b/tests/test_report.py
index 0d8714c..007037f 100644
--- a/tests/test_report.py
+++ b/tests/test_report.py
@@ -2,7 +2,7 @@
 import sys
 import pytest
 from unittest.mock import MagicMock, patch, call, ANY
-
+import numpy as np
 import ecoff_fitter.cli as cli
 import ecoff_fitter.report as report
 
@@ -205,15 +205,29 @@ def test_generate_report_to_text_two_dist_verbose():
     fitter.dilution_factor = 2
     fitter.mus_ = [1, 2]
     fitter.sigmas_ = [0.2, 0.5]
-    fitter.model_ = "FAKE_MODEL_DETAILS"
 
-    r = cli.GenerateReport(fitter=fitter, ecoff=4.0, z=(10, 8, 6))
+    fitter.model_summary.return_value = {
+        "model_family": "interval-censored log-normal",
+        "model_type": "mixture",
+        "components": 2,
+        "wild_type_component": 1,
+        "n_observations": 100,
+        "log_likelihood": -123.45,
+        "converged": True,
+        "pis": [0.7, 0.3],
+        "aic": 260.9,
+        "bic": 270.1,
+    }
+
+    r = report.GenerateReport(fitter=fitter, ecoff=4.0, z=(10, 8, 6))
     text = r.to_text(verbose=True)
 
     assert "Component 1" in text
     assert "Component 2" in text
     assert "--- Model details ---" in text
-    assert "FAKE_MODEL_DETAILS" in text
+    assert "model family: interval-censored log-normal" in text
+    assert "components: 2" in text
+
 
 def test_combined_report_write_out(tmp_path):
     # CombinedReport now writes to its outfile
@@ -235,8 +249,7 @@ def test_combined_report_write_out(tmp_path):
         individual_reports={"A": report_A, "B": report_B},
     )
 
-    # Act
-    combined.write_out()   # NO ARGUMENT NOW
+    combined.write_out()  
 
     text = outfile.read_text()
 
@@ -253,3 +266,46 @@ def test_combined_report_write_out(tmp_path):
     report_A.to_text.assert_called_with(label="A")
     report_B.to_text.assert_called_with(label="B")
 
+
+def test_generate_report_format_model_summary():
+    """
+    _format_model_summary() should format key/value pairs from
+    fitter.model_summary() into readable text lines.
+    """
+    fitter = MagicMock()
+    fitter.model_summary.return_value = {
+        "model_family": "interval-censored log-normal",
+        "model_type": "mixture",
+        "components": 2,
+        "wild_type_component": 1,
+        "n_observations": 100,
+        "log_likelihood": -123.45,
+        "converged": True,
+        "pis": np.array([0.7, 0.3]),
+        "aic": 260.9,
+        "bic": 270.1,
+    }
+
+    r = report.GenerateReport(
+        fitter=fitter,
+        ecoff=4.0,
+        z=(10, 8, 6),
+    )
+
+    lines = r._format_model_summary()
+
+    # basic structure
+    assert isinstance(lines, list)
+    assert all(isinstance(line, str) for line in lines)
+
+    # content checks (formatted, not raw)
+    assert any("model family: interval-censored log-normal" in line for line in lines)
+    assert any("model type: mixture" in line for line in lines)
+    assert any("components: 2" in line for line in lines)
+    assert any("wild type component: 1" in line for line in lines)
+    assert any("log likelihood: -123.45" in line for line in lines)
+    assert any("converged: True" in line for line in lines)
+    assert any("pis: 0.7000, 0.3000" in line for line in lines)
+    assert any("aic: 260.9" in line for line in lines)
+    assert any("bic: 270.1" in line for line in lines)
+