From 9f07c1bb184731880e93ae99028c38218cd3ddce Mon Sep 17 00:00:00 2001
From: Bryan Herman <bherman@mit.edu>
Date: Wed, 10 Sep 2014 21:55:22 -0400
Subject: [PATCH] added more content up to toy problem

---
 docs/source/_images/loss.png               | Bin 0 -> 29146 bytes
 docs/source/_images/prod.png               | Bin 0 -> 18894 bytes
 docs/source/methods/cmfd.rst               | 353 +++++++++++-
 docs/source/methods/cmfd_tikz/meshfig.tikz | 628 +++++++++++++++++++++
 4 files changed, 980 insertions(+), 1 deletion(-)
 create mode 100644 docs/source/_images/loss.png
 create mode 100644 docs/source/_images/prod.png
 create mode 100644 docs/source/methods/cmfd_tikz/meshfig.tikz

diff --git a/docs/source/_images/loss.png b/docs/source/_images/loss.png
new file mode 100644
index 0000000000000000000000000000000000000000..37d4169a18bd6cb5719dffac3320fa78365c6fbe
GIT binary patch
literal 29146
zcmeFZcUV(t_bwbeD$NceO|hW}BS;Y;5gZFrl&XLfr3gqDq98>@5Jv<>L^`N|Ap)Wx
zB|uO_suXEbLa2h2L~4W(@~yoCGw*wy-<<FKeXcWe%{7{~pZz?m-0NQV+QH|vPw}kV
zz7B)I@SHh)LKlNsF@*fbwH#jg^ilmi{BQN;v!_mA*vP+$MXBK!j2Pz3iR1b=Vn->?
zE~c$FMHk*%a7C}&wBp>Rs4e^R-gVfV__ANEa#L#B0fPvlWJuHZl#*=*F`BwMUy?cr
zazCZR7tVKz+s>7Tw2M3Dp6Tn+EazGKK{#Z1U-<W@FYS({@7;Vp?Q%q(eB5q@lnO!q
zj-z6K;z}0GO7?heBb_y)<&zkrzN)3$3cO3amyH=O#ILFIaUvKJS79)T-%rl5cwo&6
zd-GOcaLZTXF__c--MkpgTK>*u7_nt{lrb0yoXZ-F|4J?z1`~$)D2%~u6#IW4{=aiK
z)VFs^n-|YF>FH=`YdcYkSq-N|efmNa=2*hZFu1+R#_26|dQ$7N$VfVwy-*z~Eqd&R
z3v*^Hh*#Ka@~eWm_gIcHz3@Gix*dbbo>5^fte)$}E|jxx{HTeFH20lN3KCf-KL0(G
z_1(K`GTv`inmv%JI!*BPbHHHg+hq(pu@jF&Z%q@;)Ac14%x;cml^4(CWcsj1%;|l2
zuvcO3E1}G`_qKR?*ClpUyoz^ioHC=u#CdweWI4vaHt>%;^L#h!%a>m)L<9#18>fGG
z`?mgPX2`<;-Ji;S+u2?#vGb)=rOD9w)BUn5Ff#msS>-JGX8U9zwrP5xE}^uwt2jTu
zu$CLcT_isGSVW>!pmXn@J!#GMtOM)UuTR#syLD?m#V@y~)J*bTbF(6~(T8!*$ovwv
zGhfy;kAJX+vOgk_c8Ci*AtUffREc{)XY7>7YTy1hehM@N^-!KG#2OlfR2Q#e!mb+0
zB>UE>Fs@Mwj&uau=bV$gx9il-TRN%Q_Z}B8Ug$Q!`fx>-+_X^a!i{9zn*$L>ExZRJ
zsw%_xCfw)^j+xJCXlTjAvaZLD7Qg!zMj8C^DdW-x?e&-voIruO%sqZMBL<D4jK%u>
z>i(FOm8BlUd#Uiocukb7vghc~{I;{nx*dKOnKSp=+PJ8lDOx;Z3R&9Tr1^Dj0>2*z
zedb?i^6j_mK1FMzCUXT&u}If{ev_>`Hu>$Ww3L*T%CC|Yby)U%5u?kyk!&CtVV|sf
zmEn9ab+oJ1mmxlHo+o-!fV(GjN8m*2%r3qaShu=`Lw|O&s=Jj%8N3RfpVzf%wfT4*
zja3!y$!_E6VSi#(lv;54^~*fP=D=5fKlYx%@4sU?>Pl~^(d<ap*j;-c&w>CXx3TlR
zvgaE;&rC!2O4^Sa42D~Q{6}SC&o)gxy`K5m59UWzMUG;Frnhb(P=zt4qYa~)2zk^D
zH>FM=d`GkW5pDPCF9pnSos)`E7z`n6yehFwb^4L=r-%cNY&C60S{i8urj3x7s$Ev3
zJUMCE@!<nei9$?+-?I4q{f2eDXTnkWtr=#H%+@LR<f_Ep_@i4G0$aB!H&ebwEyE-_
z$eRgvS|y~>crlxwKky#b@9fLBqflZ6F&L-F5Tn((**asUcm><ln4wKu7_heTmVBBM
zSR5v|?8||>+)RDAP29)V$O=?MksIWX=#Gzxv0>~JC+#QIeJ#shwH4!kXG7;?lM?b*
zL9Q$C?Gpa_;in%RJp(;xc<e339g6i|@zcU{lx6;prqG}AiyNWoV1*y;2zJ?mY~+fs
zaJ7sXg77dcSyw%)P-qw~D<nxH=(*A(%F0|Al=MjZoWRRB6`Cl^Fsq-!(l^Wy&^>NB
z8zX-&mr5J%pH9`qufc5RY*%j8fQ)RTm!0+Qc8OE&HjoD&!Iv6l2uwy$TU-5v{st<q
z2Og972;CLuwo|Gb)$Z#RsCpiJlL5Y0_O0S{_Hdpr9^=34a-T=_$&<yC(g(;kQ$in{
z#^mOTV$QxMttfe`>#7s?fGI^9KW=_;Ep<}GwKo6VH?Pmk90z2t&gr<2t;54d*oE&l
zo;=kT>*<>1_}*K)&BsaGa}DP7qi)~oAA>@?k&!Q(De7F_j1M=m48g4r)fvRa$ikV3
zNg2l--WQ>fQ(5UUP&IdU1CyeB%-DlkETb&cvk|t_;{H!q07LpN{Z?m;$cxQ?kU}43
z3$Ia95!Oz4<QIG)broivjtlyU=}@LHzTW`vthCEb>rGfXScLi+cuR!QaTqH~EB0o|
z^u|Vc)*0}K*L=Mf9LyCdgRJ2A9p;#1q<fUC`{;$#QM2^GOQfBqWne%4O<c4~SJO!4
z(a2jP;<65s5vIea4{VMxOyxGl2M#pReW^Y@%|;y=k3xg)-QzN9pu<(wLlHY9J<}sI
z33@!Kj=Je5iKJF4FA=Od@%6RB!$a9RUkS7>%RKLMaGr<HA$#mCAv8iePcyh=%H8oH
z=hD}oX3ZLZl9e&I25<wHNgsPXiqGVax7nL}K7navt?SqFz}ML<_y70;b&1z*(Xx#k
z%P!%6J6_do$RKE~n#i7`oI0t-s(THWUVRU>1oyXj&#PmjVrEFE2&Bx6Xb61gaNtE&
z9jVBlHwSY*g;+Ue9J#hKQy<LB|1xKpwR1;o_h?+1`=TTZ=YyI140hc?O`Wl>)qJx&
zUr83Cynd%%YH!yjpD0;bg>8rjL^L=5@Wq2M`p3VA2+-kUZ~miZhw2XA9rs<fz;BDW
zf+I{WSo8PvsUF!<n_!L|{aBO+e$q`lspm<}m9D9U%VEN>A5}cM$lSfRI^EMq3XW@s
z2O^Ou;D4QqMriEoau*sTo%)HFjug^1$uxAhs!I=qZ()3P<}uXsckJ2Q%D2mg5f#;P
zHpe)%<AUtPi))j=&GmGe#-37LhN*wt-L+S#nG!1O5JSnl^p~vFp69R~^=rUF?^6Yl
z!(3olddgp$H_UZh|JPp|YRO-~N=`ozL81zy8KT_K)qqA3<*kNBLfeg7OpCT_Bwvk_
z(qhpfG;T~(DuS_DtPslNo@h_~p^xVizt;B6tFyo2^o5JT?rULhP6wbHzg{fWm%;yO
zdZec@!O7LcUN78BLB54D4biW^1orwCgD?xubhSbL{TFwGWfoV_&{ZTeK_ms$dv5J-
z)F~G{@)C?cOcpNYm8J7Flcsw3+vVDt=nH%}1N<8Q<$9>See5>^U$M3M3L9`+2Jz<p
zFYx<aV|tTBQf@U(-0Zd~MHsX0Oqq&@B<-FYR&<7T3H$H@-kD+A%&(TbanGKerfP_V
z@**ZZt$Kt!m?nbl^7OcM?BKx@c@GMdF`FO4Jr(*1PW^qb4<jfs+8q+9qaq_)FJQVq
zz{$bE*C&K}YkM$?tq>y_5`|;->Rd0Ksw6w$*3m-~53oNn+P<Yn5?SO4-OlcXbZ4KD
zO$on0uTlMVJKZGrnMN{p?j_cK?hZgC>~xb}czAe6uCw1%!zt>GzK{#VdwPe0Dl4Z;
zFSPl%W<fI3cdMlB)&8s`g((T!@8p8T(>lld%ISH1kN9qNuMjedjEuxiCGo`_yE&nk
zFxNDw%E<p~o3!fh?z@qPE)DzUB!aiOH;y8C1b4jYfxWf0^_Xb9-@+{aCk68C*SY4G
zr$Phy&3tDEa!S{4JM{3;BUR?;drGEfJ=^b6cz6zf3;SC8=E3Za4@Y7@{H*(ZBbNB^
zp<JxzV67rVG$$=~I5%a#Q<;y|t4795?_bhg9<7)7RHo|NdjNom`u^_OW_BZq*{oJg
zpcc7Pij3h{F+-%@)r{KOcomX2gEnV0Guq}(>&(yX@*-9@zZf}7=7eQ9uM#PZBPecr
z@rQ=Cre>3U&b@p04jFB#&nzGz{=y2kccDk>3KcD=BEBAqzRGWKzIR<#^7j7)pA}X&
z*S1p*@b{dmN|ScCfh>*d*|IN#-yrRy4*(+Dy7+Xzf$3Ir@?gC82PaDVfW3xx+^`PW
z`)OD@1SPS*^zSnJLa}B9eKP=EIi(uOZWH~Ld*g4sU^aiAN~Zca9f};y0DG=CBfKRm
z5Ss1Tvsm^Z|0m~#AI>5wUQ@COvo#DqHiP~>t515Nj4Ysb!OX0Sgs0EX&wsqTru@qn
zUg4wH$9v0UZEHurXT7X$=G91ku`p85_Xx>3Jt<l}W4XIxJ+B=!Eo%4l;D97aNV(0;
z)1`6BrT5k&Fi%d>=}Aw^W(LcfIhL#7+|@bN8_ND2>esX0`~q|EY~&%8$r@SdqN#82
z&kDX{zvd_FsNVD<h2nQH&7SL6t_jiW2J13KXZDcF=W-jKQW#_F-@kVS`?+P56waMn
z!i+rLfTj5mNTQ2j5gei7VnaB$LLaMI>~<(~5&5G2W0I+F&(%q}T1Cq;cE$Mft8F|J
z_z2wW)wpXN<dvoedeUTc&&FRpwz?|aOgNk!(HI-$mfw>-ZJ(o9O1X=iv*TLY-K%}>
z2B(n3ZId#s`(};oJ;O<gH58>EeJqYNQz9e#njj#qmz7N#kPR0SJ=qZv(#ERX`!a3^
zMobF=<LOM<gJ`&_I6ZPP%p})y`|<E#fE4#7izL){ZDHJou<6ZZ>5()2(g0#=o2;IB
zdn=?!ZUu-<x|Gx(5P}qX5wx~1UO$|hn%t|p5(1g@p$|^a3C(JR%}DlGi+~$7E=gL{
zt0Eo8G?H`Ij%Uw$V{cYJYe}8Bwg<ECNp?85Bu!lc;)Qs?NUKnaPS;yQe5kG%0{aBK
zAh&UUQq5tw+o%;T&xh3weG-$1>PhKUwQ`G-ZLEd(cU~3=E$W&sqN+j(#O()2x!xJ}
z_Ze54A#ue=0V0Uqib?#nl9801MN>Vxy(inD&cHF%IAqg1h1w`2=W{Ph*B0b9w!NV3
zZj+MoX?D`l=*p3l7)T(3JjPE1RbhDt(QsC37o#uPFo%NIz+eDuck^`|z2@xRD9LIK
z%<t*zv$PcdWDtF_4M{#@JT^}R_rQ0ZUH^zwl3m2d8|Vf{kQ?JgZfrvJNa`#N`Dzwm
z6F_WI0{s%QJM=r(iucRg+PIGI@xUYd=M;imi>&k2u2ypp_2~LdcN45u$30N?;6R6X
zYOEk2)e;PVbc+i!1f{Q0Kefn6M5zDrRwC+Sa&PTwcrD+K?2RmzMFRv|f&-{wkvBtt
z?JKbm%#e`1uit4`ilC>R$j=Bc2w>9mj%_x!yG;LtY%<J!<%!@yN!oST$gaD-j7@uc
z<x-sAkVYVYJ=vY-4gvpptSNR>%QGQn{i#PM@~XZFK0Xgdj`7lV8QeOROk0?0%tzA%
zYJdVWwQsu^q(}`Otzk+%owuZ3XhtGjbqn6<qOZHnunvL)_C3QRD-v^s{ra977-mB8
zv6qa+QoZarcxvNuim-pup2}C+WpNMPh5<j6)@Y6-4-8zzKOIKaW!)EX+h@2Sc9^yZ
zt<Bv&`%S{_bMHR02XbTQGkB}~5d@5Lk~tS7BVHJkH$x?`vh9YoMMC?*)B<pU$saAL
z?FBiFyNqPd<JVYJoIVW5KXp9-4t1{T>hx2E<|G@41~RdbOW*1^N|E79jvY2Nkgp#~
zAT?EE&g|@4JbljcGrYrNTbu+y`pID0a*R6EF)Nr&+Tb_NZQnc}F{iF8`W*eV8|nyw
z<xeyLi3UrMeHkOICJHqH{6Bv5ZmvL9_X&_#iR$zZ_!=a8*aD1GutfOLE%+BznY1oi
zewuhde!gO;A!2I#WWZfV^_&fS>kL$U3&-pM=|Og1q4DGlVDlBC3~m}px^QlI&hpJD
z3gDbcP!?|8y3;gyKdIclY?*)Y6D!?8R&`0aN(}6kPbFeJuXNz+B6zoOh!OdnF711o
z?VmsA07CyeU_G!p&whQpsL*<&%nz3k<idmQMEHC7iNF3rYsKa?u+I{I&ZlPxkbKBz
z(4F1;W3aTVXe3C5h}zcbdtm-E1Pato{tJ(-bW?k^`KF1~OK{t#6VSI@q1?e?Ima@%
z`VCqm)^9>h=F2NcKSjGy15vGXuC3{}LQC5#1rSUU2hkG?$&~hB{QT)McRhdzzZA9^
zwTvkP^Erenc)V*v#tZn7sW9|fH|-Xq*v{Aw>XagVdoa<L1eypjmJgG~SwALvm{wh#
z-%4=+s$xJFpMcg&6?f6^&sGtrtYYRTV&(@kIIRq-uzxS-rY!rRBpKE(wHs-+=P8?T
zdh;eg*N$}#vOa%IxAx<w)jXv;u2?hNlxk|S&C;5WQocl24%@re)nr^kl95=Ro+6YX
zFJErZ_8z^^G3q;J959(VZfbxLQ-E+^FW90CA)t+nt*H$RalYJ{+YFAi?#z=jms}6f
z%oipi{pM2r=IR#Ouq4Nh4~)V=6)bxa>o+?wFejaNYj#3UN9Bl;Qoav!26(Qi5QVP$
zIU4K%IoluKv)&f8*FP*y(O;LWn_asVDL8g)#`jHk-{o7NvYCF^`5r#14}$`KWYc$@
zdenJyFwV9%cD9So>f--2F)_hnP?780UBid)X!S|z&c(Ao95RbM2jlRb`R%uT>@$z>
zTn+_56Y$ph?8rvn&q}hpcJbW11(zC9IJ4(w&WpLb*xz^OO6&X<CdC8H+0*9c%E!E>
zu5$GZM4Fo^7rOL(8X6jUU#Hr^vS0J(LhF&Yv|`fMPjd-cF)!j~8xt}@9?8ev8q*<c
z9^06#TN;&``Z7)s^Df*<*QnDl3NQ=%7OUqjpQ?b`?>_88Ury;Z^IPK?B%0SyQcw#o
zK+xG=rF`UHpVm_{bq~TeuNA*F-`6uYy*2aHE%|V;bKJ2q9-3D1L_}r)zwc0x@71wi
zqiyd2>143)sCduBD&uMCxn09<5uYmU0PcIf7sHMRD-I%A@V-Pyx9>nzg!%7#=Fj8f
z7lIRh2Mc9t0{Zj&`uX#TbVHL17lzy-cG{<rq@_<)gL_HQ40x#Nq*H5wIbwZY3#BH;
zNqzZ8R(p<v=)$j{ZQHil+S)2Owl`6oQ<_ftk%>uR05SND>!L1hI$#Bv562gL;BMN`
zN7mf*2zE9cJB#eVd>nN)j#11gQZP$+{`|EeaoRlp?b~yb5nsy7Nv30OC`HRD^A9`I
z-13v&@4ig-4mxw#-92fXA)9Mwoowmt(i;FS=e9G0J56m3BRM%a)OYHo`D~Su??|q%
zgO^u<uigKMIVJCX$ot+&7vUs!t|a4-a7v5@GvYzKAf~%zi$WfmHPgXl(6Kis`nMJT
zx`$K&TEzdE3ZxnBFNlTmRk`IQS^&71JWvi|+xpGe`mtL5SXV@u-*<%`bL%TZZqu(v
zdj1}2pIt@XgM)*7BF7Zv<voY$lc1i6S8}&0e<XY<!?QbU)vCjs%D1yHe|kkk!7<11
z3w-C>vk6H`wj4Sn8%k7(@PgiAKku1QJNv6^#pmbi*z=F+yx)uel`QN#K)UEyUS4#8
zmbc3^|DEsi>z)bmys0{;fR4-FTluP%8F@R@yGQQcs8-CU4(C&fW8W|2dF6}^|Gshm
zu#&##n#5d)xi{H|f4NG(NC}HG^b{3bJDu?Q!oV%e)jtuxch<gV@&r!4tn34y0Z%@K
zXb*9uJ|2xzWz9|Wctc7N!S*Di*js3DwH1?7G4in9{qWIA#W5|QJOY*ERg7|UksxES
z3(YfsWk4*2ErvMCAA{zk)fHa~y!qxR2O<<%Gzd@#(RLJ#%_+wrwq@@`WAp9C$N~hL
zpw!wO<-4y<d4Zb>QmLk$SM$^ZMDdsV`h%D-#gB2znj;ICf7iW2S3Rn5wimqSzR`=F
z&7a4Gpp@MSPQbq{&nbLIwJ5C5!7lx$1BgQa3Wk9hu}FwcD{i9jiSn-D<<z^6;TEBs
z?e2?zOCY!`%p(c0>w3G_5|x4OeL1*t>4~M(nLVJ7r=RsS!mklihMeKmDYzYce?N!A
z_)EUB0U|7lVZW#H-r+!yzH$@kKi@nel2+8OLemU)6l#!LZEVyCO_vF}!7d9&5Jbk0
z=6DYm$fA7W4<80R6Md3yq8yS~c><23gR^N#PnRJ2eTGBp(GcG{coP%&1U<mJdjN+t
zMmYmhwwN|+`Ky6Xz)}k=Q~pL)0*BF~;7yqyrkQ*QNKW3GiL8we1s{T5Xecu!cNJt7
zA>?1|m+JsCW#rJbDQs1gAiGKJIck%60RSp;aj-KFlmqzP;!$jlKmvpw5m<gl`Igat
z&^RB?HrpQprtuI$4P;|*mRz)Rds7RT?d}F01jD|7YY^Bi677(q{r<LsK;?kGzu5_h
zBkf@b5a1Q+m$OYVb=o(mPzi<#=ES`YHN&8@KuKg}xlG%47{CJDb|FS4+>ZQ|Yu2<2
z;|OYJoLk=oA6R`4wbh}kB`8aXd&jv$vn8nWK)Le{oYTQW^2a6!NR^a`5VkXfQf)Ym
zr)b&_GLg4HfpysCkHYWvMG8fm7#e;;c%?2m5*?vHfW}9-RsR4kTG)A@1aIbV^I&M~
zM$9=D>?&*>)FDIMS1xI@Xj-X`hf$ssHD`6m(y*_9RRNI<0)2kvNbW!A;1$j#s3Rwf
zTC_$**4@3Q;xU)e#Y>K~sNioW`G6O4);#=c{T}1zoFdi<DZl`_ls>G`B_uT@3u&py
zqKysUxDszH#f#LOlrK}S-h021_rg3l1s69I95VhJlwZC66CvlqxX`Ws`HFxOaj4bQ
zJD^>3A(_Shsp9l&_BBW%GWSu)r!E0p<^@LDvuaByRKHEQ#UiXwvjH4&VmG?Y13%oN
z>BWA2&*ONcu)Z6x;ZO^4%c0)E-~45E3DUTqpAXdCZ|qSB#EC9axRqHV;9fFU5-lF3
zOpFuk5&G(LegLcv1-K8C`^A<j%lNtL>N(MUXrXitb4R)VrK!^mjQtyL19h&+Nyiw<
z*<{5!gG(3fu2v)QY#+35-NGR_9vM4X$mtyww0F`GHp3^mdH;$y%g#q&5dkwSVB4+~
zAXh}<0$dR(7Li)kq%LahxQIUzQ#S{q*3wx#|Dz@w=(zwN=BLqX0cW95l;P*>l1{&Z
zZ3@L576atg&1?i4pTnRA`YjzXWB}-&IUL=Qmivsyj|1)GA(HBXVWW9Kd4s0by~)8=
zK0t^mjq<?}l<6AGD@|At0`l>_6;L~xzateW<Br1he!2$5SHmm$ot}B#NJ^dL-;JDy
z4tgHDQy^ApIZkO^-hz}H8zoTZnX|e=QK>?b#al;{j%Ki1CE0l5I?OB@s3LyP>EeMN
zka^0v9rNjAnprw@h;VM4-;@Sd>PM=E2g;jqWR8~|eIV3H21qm?LLOJ}1Q7KIUw8$w
zYABB<D{4s#fL-`aXD`rO{XPy<M<#d;6Dlh!$2xMcb3-qk3&&5sF&Mh*rlhHzAv!Z^
z)_m+8GkB9La6V1@Y0jFe?Aba5l+pdzOuz0mzTICHsPoO{tY(em38d8uXo^?-ptLmA
zkIjmdHdnpzt9xdwgHQPA)Q=eF;+b{_eAi9eNTTrXRh#lj^x|K`#F3ygjG7x6lDLnj
z4_6X*_l7PEg{pA@tYx!uO8<>G467#}+9<DdalIw0E=9{#o8)K$RKmi;(DMAMiqoZb
zc-ljEW6WiUXb9NcL-%SW6%TqB+cu_XIRr@;|L(nnr%czY^Q8*|DLB~^B<kC9Z(D)+
zJ=GB-pyE*e7w7aqq%-ppX=cpXkIBdW#mA;mS)`5+j?Acp*{GaW=%b?4VrJ)i+q!)B
z!7i_n^g~8TY)G4`3Xx|azczwHo!_jIhn=X{MP>L<83>>@pYPq~$ISI(uJQY|rqzt>
zHc)k<Pi*+_e(`|hW;OF+n!EdkLmO&<%r}J}@SAy`-~%n#g;KxYd~{z1Nz{8ZE9K+i
zA|v~kWa^f?1ONz+PqVckFD3eYjyP^G)xDsOol+N_?_R@hq0X5|ljkSCfB!Dc`X-%T
zr25An=hF>G$ZnRFN&1o@p`rNv_KPI%RR!~ZBrR3deDA!|w`wFW<@cdYuGZ(asN_-K
zdW^IHkH=uoC$Qgxp)+W3^DTZK|7j=ZR2%Gjrk6UAMrzINm7D#%cc<x2KnVWN!7z;s
z4Tt%pFBM!LZXzIX&>`sGFlCYp$~<(~FpQ-#;{_p|<XtL@D`h9~`3%uW&}L$GF>3L}
zlQGo!9|#={j3cX-K3{9TQ0=D_@V2n95PI~pv$IhDJPHW`!gHv%>@O&g+On>k0_(eR
zZ5iu&xo}Ta1d!l$Htl5W4s4f~-AIdesAJWupKY4T_`vnsdSeCsO@Ii=$fHwjvCKB4
z5`shX`?c3^VlPw|^9`b`2J`tcbHoBdr1Sm_b_!YN2#ez9t|2rF+4uY1QIYW#)J4v-
zL7TTk3d|H`1L78{3ABZ%x_}jHj6YB1iQe$Q|3ekTspT#DOV&GJP^?~g=!%N+WyQPI
zKc%Fmmi}yninZkH65`T6kD{DwZ5y9GGcGifJhGDJ1ORdhNR5JC?TBUVRS>JVpqx5K
zy!SOG##{l&=Q%BxAi04l)8-@WIrCg!*3VR=$&HX`zl6fM=kEj~g_I_j<+V4{2P@{i
zY`{D#q9DoB%~z|^HlO!p(Z--k8h4Uxrt{XT1Ki23!u2i^4;#E~j_S7z{s4Xp$Sv5`
zJoP#4xGEx*1$`ZooWi!H_DSdVkHV|ffHzcK7>9<gXO8jL-6X(IxC3a#eX;?%jF8wA
zQQoCo_9uP~x`r)y0WnUyO3^VKfR~P{AfO<beVj<i18>Eg)R70SY2Pt43(B}4G~ZN3
z>K+LBVzy{jc=7~-^q5FZHo)`55;W+R`UAJ%<NC+7juuY?0D^A2@=Cky?VX@mhOz>U
zqaxP)Oh~Bjlrg2SeYacm@#bc<KPg4Jl_PDMfI4M9XCJjrlzpfV>jvzLTyEV-G-h+(
zMzM!?6WvTQT@<J_iwh`4rqK=?E8v8}Ma!$#-jGtgGt+?DlDV}KFR(|;-R4=O_hZVH
zd!s3%zP$lZvIRcq=2cIDZo}n!2+9+<jPBHz?yL5ibS?5C0cLO#8lUWq&^YVd&e?&d
zTemTPnKtL=qh9_7X?>RPze9+oh~tfriGuioR;)przbNO1Rz0>DEly5UaOORGb{q)B
zCF*?%@d%uNg#-EiJ5TmGkTz#*85bO<Fn3SGEsNj<5cIo05~wJoy${U0IuY%g3Y@#6
z%&!Klm!2Jf2c5Weva;!sk`d3_iWbaZfmZ=<@7&ONn_tbz-R2QO!1H{Z3)I;CX}0xE
zhCJ2?PR9TCOB6wH7b6IQr?l1iJJ5?<(VDQAL&(J+D4$D(qeLpLFG`cc_9GkUtr;HM
zS@CB8@F(}?r(w7^qTow#D>M{4LZbh$T`lH^9>Q%N=s3QdA_ehR4Ycsc<{0mb=sjkH
za6PCnVTQmRi9-Es;Vq&g`Dt2G({RWVQB#~98<O+*Fa>=MdnXXp)FVql;X~kg`e02!
zCxRbjI#K+v=4h2it)_C8$O|jLH!DuKtZ_==lx)@24KxEsq|ZoiNaOGi2CKx*-%)m1
z4_5qs#<bZ}GD7B$4Ly=j41&#tiE&`>cf{tL@|e-%<8R;c_ejoBZYM#N3$Q)U8a$$U
zAh1!MU<i0b>~)^)<%D(O{f7U+FAhbc0AB1Rc*TG6i`+`+K3QIZGzE?KVOq&KpfjOz
z=zO-w$vUH0MCLQ<eD`c%H@=CZmmN|g&r#rd0;;IPsW||G+0Q9Eu78F7tl*H3^}ztu
z%1{bg`2J_)^>7sIxwr#;uA>vwn2OTp0#`Z5^BPUL+fs3WByqRb(s0)nY4eMoST_l1
zAG||~TN8_Ai4MM3p4>K<IldZIl1!J1Ow~Nq;0ZZIxjG!xw%Aw4v6o^fZ`)KcRRvl{
zRN>PJuovlvaZVF4vU=2$=14ysAUjko>!3@#f~lV;jD#3t*I`FWyWL8yJ+Lbnl!N3&
z!7D5%%C6}8m@y!AgkN8dQnXu~%%P@YF-u3UU7f|&+-Exgt-JaTy6uU3-O@;+Nc&tJ
zDPr{>$XB{vp+L6@;5No73f+(TbU2zM*6S$rd<QoU({R0L*jl0T6k(q}^<FO3OK^H^
z2&w=J87>;mIqqmgH3En^=zMn_n32aF<$8CeX;P-$MJS7pv92F#10LAltYP@+?OK`p
z1aKO(F)&zT@*!KLx!4D~{-1V7&vB9tGy}u-LOik4N5ykbHIftN2nlI^;}ZkbIx5pm
z{r&xb|A5_WQ0Tg`t3|SYy_N?9dI=r~0}p6kJXUc-Dy}%o7EuqlJ;S;9jhAinW8nSp
z*RN+Ulno3FK*m?!s}n!p1X8{HIv*M&az1E1C!dWhFkj`BYFbOA78UR4vO*f85D(!<
zCJYXI0L<1w@*JDXF=yt0{ER$_U#Jaf4HZ$@u6AQ~AhMNL+B8q*=+U{+_jR7+_T%Gh
zkK<4N<}Z7Z24GLNmfx51KGS4<d&RM->#>c?`=J)U)Np{t9uj4}yhQS3)zr;vHjex|
zoLuE41r>PK6{yn2>~ppwDbCshi5p%o>p!xZ2zu-p>Ox1I@67vTp0PCO4O5D|@vH+t
zXej~zY2sYeJRn)WW)}FMg>0?M)=1}SC>Yw>5fKrgKBHGe*}u!FfTJ8ui%=BF?)#7a
z_3ef4piL-zKQN=yp+q{hKg_w%$iQ=r*a!IjV1nNQoKM1h|6V35L~dVB1D(DeO}a;I
z(drr6TsZc9;42e9vS7bKDA>6dJ6j3Ywa+n5H&o|3dlu-1>9JfYuQ;%#Z;Eb|e)^<F
z(BCd9N_We7@5xwqowaF&ctGbIN8S)G=x*b$y1{6*!STD^`KHBQ)5CQOqjlQ1T6=qY
zo#&yT*3`W9^Q`$g@zJg#Ca}2E|J@{)thWUVt*Z-i>}bz8SwovXA_b($XTADvUjU=N
z`3=fl<B!NQtqbju<oN*xnavUpC>QlB)Fl*7R2^_IHz(+8#~+*hyh-Pcu;=@AQ;6G@
zm6=%~J6dt-uud_h*oS!vJJ*o2VN^T-jIPL=QCnAs$V;zX8wK`!o3z<6_a8Nngp~_S
z$4c(rNlHcAnj=9m3PH)qT&cpGt_%egV)j70{K&DD_De+NC8P|gf56mo>?)l73km>W
zy?@up*3D1pu%}2}(-~A2g<ABukaTg?gBJu^ap4YIq$}y2`~-waoAyrCVCr|HAvcK+
zGAgQmd*sqy>%kZw=eJgUQ=oiKVkzgU`vG^E{f>zc2K19}40rgCLO@7vmt4X9xU9_d
zN6plz8A21>|A<Ddp;ToWmr=$W(G(Dh4Mzbynmdl>%CP$(lN4{4Ot9L)Wv{guX@-(e
znBojMVVe!jaBtaNhg-Cg!J9d4p2p>ixB{n(_G*fT@L1QKJCrs8N#nuDSAz5EJj%I?
zeRGFxfBfok>Ba>Rp<G8WDF;KI6DvvpFmyz&;JSJpH_DiyNJ2de;iT-&p5yfK(VABT
zFLv=CW^jFs{GWLz5MHikC0@*$6aOE2WZLKNu(2_gPS5yqQ%3K57)2-6B~ub3456?W
z$>~;%E9LDr4ErO1Bw!^aQkyz&Pr=QBFiOYo=u-`;XI-<kiYs=ia~(Lq6Zt%5J;}2S
zQq$csm34x7$t%`L0WV=$iYXWRQeOH=Efygk(8Zw{ejmb%2k?KmX2-g32=;#+@=7Bd
zL>T&VhPm^l>zU2lqnk*UknV-o(n>7QKE<<K;*UsJDAK;=xuc8%)lF)Oo^LFTT-1@G
zOfi?}^oE4A6u6Fudha<}Dz<xenv_d2I`WxB&?hU0tw8Zzi9gc0>~U(kaQ8c;-0vWx
zUu~p4Nl^ze=Mi^zSCRCB@itA!&4Dr9yxiGcSfe8W=hDBdZ36&t?i5jPHT<vvcpaT8
z`bgB7w4M0JC9*RA-)9x}mPuQmSzO>NaDb$sl96-f%_*W+h9v*prNa#S0Lck$ZQ?w^
zJg4$jsP6)CgqsqAr2*9spoSP2ie@qyF35JFA@_)&9R!ABl^6;a1P((|-L}nH7g4PG
z(@>6i-QAbcNDm(^i~w|SC;tOFA|62WcZ%Ws;#JVzpuF+}Wtq=&%?Qne^F;0UlLhoW
z1CfPp!CMSP!A7p*AeG}SKpb-V#Rg53+a{b<hy~QcF~g2+Zv^zry^Mo?egPY^M2~ym
z6L+C|j`08};@f*d5rG}hab6wGS<69Sa!_3?fG|kXLQC#i${)qy5WiFbKbdZXmT7X<
zj0@?CXeLJyWKho&JxY-@l9j`(zNMB6Cx(CpJw>@c-Sn-4Nohd(mtaujz9sNH2T_do
zrAJ=o!Slh@5OV;T;H^%})12A~AqCmQjpNuOT5T%58JN={OJ^oihy)lUq$%2p2-*oZ
zmUHa7f(Ews<V=DH-ruBNj_Pw$NqbSqavz)r5Ow|u5CnFba)f<}L7eThj6}Hf9!_WW
z2FO~(L^!*YdoXn{Es$S3&c0*RO!7|<mk_djuhw$o?9nf-Qmc4RE4Dt{(Q*?2N&(<N
z4j=bR80F5db#>0@XE{lFWJ+=$fFCrr#=JWaS=!$}*HfjDEp#R_5(J0CWd(%meMM<r
zHW*zCl+&yi1M;3UD$*Mgd@XHJgb}-SpUmQt9Ds^zDN0LNxZaU$e?(LBMBWKKq{Z>p
z!6=|{d6>{g;l%ZD`Ptv)AcyE^zTXD-<A42k$eGo?A6I5S1Gtr97`sPlx&~sv-Lq&A
z#kC?arjuxVNFnyanVjvZ`N;&^6r{D|f0=VyL$4s*-*oRt8jX6$G5#--VKE!ls&iO>
zoCFG%pvkj{zaumbG19(lS6PcYhFbD+g;ssB&}Mp9K6J~eK2{P<Y$I*5GWS-_O`VrI
z^y9E?iaXp{Y!pflud7DmezzfXnNt)HKeE7~RQsBW<4|)u86%!%;eyto>#RWZil!F*
z&Sh@07O0wSDgw3Wk}-(2`&M%<qjc5>xN0)OgeGv&FteP)u&jFrqPC(e+dWlMZ!uO9
zc`K$)fX?-^7X~3>S0$q58ZOrhONEZV0B<9YyN#GDIK~4LSN6BaVtxCm)P<&mH1@xV
z7%?~7F5v4(wkb##1$8*2)Fbvf6&6n=ORE6qMxHo8em6q;`qbpZwOORao8xC@W>(Bz
zV8WEiEtpFH&QlbqA)qjCjpek`3MXq187aTjS=D>G01-S$z%iPozb^h=CLPGX&<do8
zl+nK;$j5AtWx@GOnlaaipS)w$K=!mBgzCuG&Gyh%w~dc~<Y^<OgWO(+q<d*I(a1r_
z^M6A`t{S?B`Z!OzB1qJ<bF*wNgL_~RG7)AUu3WhSdQSGQAnbgI^R2!wBTLZHh}QOg
zuruG?xoFxjl88}3I*}Z4LxHE1@h5rU2L};?n9eMwm)u@Ko+Fay?n&S35uLy5=K?xI
zQD&AX>Be|(`*6eoqdnUbqz(gDNNnD^b!&m?s-d`)0w(ktIEZ!+k{u-@xCYG|pF*|*
z<#vUE>s=<3%3df9;=|7NheJ~zyU^)g$3k`vnd@<0_)7js#~YxeW!PJ|!5v6_V`9tD
z-CvG~u+Cw6s}d!hT%0F9^7;M(o^OFgMOa`UM*(-TS1k>}EQ&1a2XxkO=KzVK&P61q
zcQ{MJXn{-#@ShQyE_~vJ(8E@c{|od0s@^l#Si}ECQTnS^fWSm8&$r0qB-t}v+y$==
zoX=j^8>&;R{d^30HikOalIuHl%J2HvukmuxOnu2rz>@8j$+}9=72l+Dtbrwo^8?Ya
zxhc2Q0mppRm@^b%K-aZBh|ZQhm9WqV(>~bwA}X^evv_Wp%x16`YL)4?rkjn(Y$h4|
zf5cLTk(L^#wzBtTLTw-XzLugDvt2s}ZOnqlATC(u!(M16&l!=&g~;QGLFg>>>v)cS
zA4fKx-4Dn;179@qWt(Zi^`B9;@i%`aM?QP@OyoxQZKM$B@}`HpYkw4~T)$pP1uy`}
zK`t2E7i~UZrC^P(Gq|@C%2<@06droKl=9KJ1E-j^;+mXd`~~Qmj!M-ZfyJPi=5|*l
zU=K?PQ&btX`3XWgcEzrfX@c^(5fxKu?3Pnapvgn>4MSAmqk6Z^fcsn|_K<e_goK?3
zexCp*!mN5mbvGE%>3C%vB3r{LgLlPlUF6HwO`U?Am(mcfHAN~d;S}BHZh$(F-UaL~
zymLzR2=unj$#6m@>nUVY%H61zA$-s>wCO7>wgGUC^r1kcx-cXHC$Wmt2CMgjGJcFh
zZiGGcMQ3b=5UnOG_XSuyB4z;y5|#+23>^2OF6{PYz|rqGtVY;D0NbdrafX1~0D`Ps
z4)=%aoio6P9_ELt*Wbhj%}?h{|BKx4QbTh{39e#MVYSH9!FBRhaqi^LR)8x#r9CCL
zxs1-^gJL0&quq9135&m}W)Aba2r1AUZYb#gQTLF79MbM=Xy8h{WO^-c)d#cmminP@
zAt~`u`ARCX7zK#qD^fu=1O5(}Enm+w)GfGP5<#8(YA~C*x)Brzpl(AaZju{+S7-vs
z61>T>H3OQSr=z-gWk3^7_Z_|PJZf#S?qA6ctvh4=zE)|J?uH=k;%>jk=_lqJ(h5X+
zo8NyYeAeB6Yw9I13$^fp+pnQ!L7~r@CkT)6R{R201`*J?A#F&s$<MJeeSqpH5$Rf6
z2qGrEs10k&khXNFatS}uo#XJ^IahC*d;bwSm<9wV_koYX7K%_!*l3X?Sg3Xqpkz|N
z!X82BNQD*7&mj{x)12butQ>`xXuJCGl2}hzV5LgD?E5{26%a@Z5O7H3E(vQ>r21T!
ztBq3b3;~P`;*J!glgiD<F_Vvv;Rb-=8SeuL?3RkJQMWll17VAOkb-T-^3%rjB`q)S
z{A2(dvnVqtNwh0?DQ8Ibh8plqx$KdgrV)Hl%-?P(fCC(2choUF_8lo+5G6O!dc#5n
zP3j0R$+PAlZtfxuQbgv+n(bqQ_JO}~^Ke#n9KkHSd^YQT)ZlbD7(gr_tOMyrf;-Ls
zsF|+_J)tOtd(45|D6x7y1F--V4w>hF^sm3ZvS=!=Nk&B1VH%MLe?e$34z(q)6;}CE
z3z<E;s*aj622mFLNI?a-?P!N{i@+hU`9V2E4+DHDA#X+f7qF`W6lwh2)Vb_x=Avc=
ziW`lJF2!AkH4gBiAme%yB5pZ8id=9*=d!#gP@epwEaYCt*k{%KjMw7zf6~<(v^#pt
znso_1s~er<D8z!r7<xKr3I3M69bv^r&A)?G#4QSu2asi+HPlwE35)*0!-B$*gHgG}
za|^toL9gP0zt5?eF@m67Syam!OVzqC5^1~=nl5X^NNA>Ok%xWisocBEbSFT>S5+Jl
zf&Xf)uToiFHZyd-_wOzH{G4~m*Q|865jeS2BCL>usRoq*nZw_90s*~PXAa?Q5y>o6
zlele*4HozfHAsUz6h#GnV#A>%;4B}YPz{qlV1y^5ts^NjKzEBl3)c$NG1QeFX&$AT
zOT9j1$SoK&%q6N@5|J(nDqm?iyVwa(qr{1GN+SQ$pp!Kxcg@TaV5-knvmc2DKf^JL
zw&LQB$XI%L-1YIciGN2Wm$T8deId_56!76})|V6V)=SKR+Ooxaxt6LuFaa!gyDD0X
z1|m~wzf27hU&COyTNqRrSyy~25Q4j`o8z!Jc?jG(NJpY0<42S1OKRNIq9rT*!*BP=
z<-&$<)@-=TJPRWU04<Q2Vfw8>W$drpd<a3L@a)k4G<H19zN|pRoX|To_>FFM2<4s(
zCG)3hVijR`hhD0&zw}dCBtM2bmf;?HY@*^ysLBV>dA6E)k4dwV$#Sx=zc769{|7ih
zO5)#v#JNRrjk`@EDBB8+&-*}o^Hw=fTa1M+szbil)KCNer)p_(q34aBPxl%c8f+7M
zJ}{ekAIxX2PaN8_*FpRZyDfWvtzE98gQ8<&!uYtNk)WWR;#h4Q^zRYHZbqI^rr(4P
zvFfVck``d05)~qTOFrAz*GD80cZJl|Fe^@}PBdrw=I9bSZwzFWi!j?;i)W@Y{fb~(
zWK?fi!ME#rgaYNMXGmV3$zvm`Q)tK*BMnj@UE8zm8c9?Z&3rb(Jb&Cy8Td??J05|7
z<9L<9+Bol#^jv_GV$x{R^&G}S?@?egXgLwYwE%tgUy**Zk?e7h6OdUQWR}rFmHEw?
zHd}TZmTBuiG|DiIx;oD{{*O*`i~kNJ@<}6+FEJPCGaDI^S!2%Z_DoiLx#~<7LLnzG
zI|oHTTk1CsQ!!Ol8lq&Pu8`0CGt5TW(@03#CxNuM74oRA<-+0tfTq~PAx0uM`yT`H
zJQS(UwQ19)X8W9lk!@Rng8BCQ_l-zmTn@>od&W=OUwmmA>^Ws5ccs=bqlhGN&}0JB
znj=3o$Jf@5&7>q3EH%qn|Aj;`UOi0&(_z}=By?wv8^AOdXhXxr1E_Q9)PnIU#9~)q
zf33h%eehg893ZL?U#|b(geER@ibz#~q8hjpr2{nmMC0@sP|qf)`po4I#>YfQAN8I8
zJvusybg>hd`3dV0i3XiQ!u$Fzbbq`H^#qL27Xa{&lC=fk*f3zfCUD=f&-ld_`PZ1z
zZ@K}Ep79&<7X^W`JI_3T!&f7tIt{Y<T(bF`4*Mro>fpiod@N(ow`ZGgk28CW>~m!r
zObaHrYv4|@v$Lm4U>t}&KM9lD&QF^?Z$V|*5IIZBLKKOiIlKLL{M~#Bi%dZv2GTQK
zj6a~Br*{e-GQIC&d|wL@TY7gc6W9q>gfBLl{+U{e$jaZn+c>;t47$P-90-Ri_A<vy
zM{IvEKskp`{l)2yoD-OBri`G|5S+r(KXxq*v$^uiVA><-4mt)wS2xf5>JjQC&o#Y}
zUG_asSw4$YKHeNUv1~g{#Mqj#H{l`F&qJuDP5?&-tvRQ6uU5@|erjI2&=-$r!8z?d
z|M)7qVouv$6WaFc*XNsldLOk9vxJzJA*}*<AEtOhuTr3Oc_rrXd3>q`EcmDEn#2<M
z+NfZ+OnU^;wCw82yj1O4`_LDkxSoRqr%s{cuq|a19eV81S99Ckqv;v8fgEW+FRllH
z8d|KWN+D#@ANUriz;WD*Uv2z;L7j0E)THW#%WpOorOi=(uHm3=iOVo&gAg^qRZ?U&
z^{Sk(^T_Gh2t=@zEect2>;X_%UGn4`|8gQ)R^5Ct_U@kre6rZBgkXJXF`?{3dX-C^
zdt^5DAL19Y>Pe%!0Xhj8FW43rza)qdPP~A$>0=OOu?o_;3{}QdBFcsQPA_O{a(J35
zQ0yavIk;1tJKB66n2hMES*@BV@#Bz2fxS3w3NLNlR{Fuw(Fl5BUbQF-;@>-bmZAdO
z-Y~6y5u&WH1R6*KMPA7pQ>i{zC_s9?z0FCaf2LA}+c}WG1oIQV+7b)Olthk7ciA^c
z_vv80dai3Mx~GT;fe^MAg<8vWp?H_$lo2KV(FuIs43Jh`tJQppXg^C<bC&%o^*=`@
zHbtQHBxbeoe~nBualn<B0pO;8M&te|(Znp5I-6DuT>_M~FTrP!d85AAC`%}II7Aw*
z9oW{zf#GIKr8XF;2p$u75;m!7AW(kJE5}$sIK(A5ktz&WUqBL`3W1D^-f#2a$PJNk
zVY@2OD7lKvu9Q%7fdL)_+I`QGk=>+`y!C8r#-EDp?H>3=CBkYMMNso_B)4qa4GwwU
z_4HcD#r-%}g%si!E;;(7-H(%k{xv8!25WG)ITV!4RdkLA6|umEvN)nse{lr6B2&TX
zLN<s=-^`G75dx49Xacv;(nMc3tb+70As+}_<+Oa4aleU;CvHEgYJ8++#km+So1I;&
zsxAJ;)D~gWzO|b7x4oMtxvvMGvx`^;bA+|-HtkU{f5{?cn!s|-F#!%n2cP=N|DJLO
zKV77ML2;`j%Nd5_*|cdx2An&$9LEL%?-eMcDhL#-`->yg45$h#_XDC<M+-B~8~>2K
zXZHAbbuL>Xdl3i^P!!;IFb6wZ2q04FUlu8eH0?v2xQ06gnERqBcX&8z=JFiJv$#PZ
zg{yL08UzEwz*vmxAZ(q)X%|r9=>116SY+o(%{jp7n8%LtRzmFSeuhpydY{L`m@okz
z`2Zz^9j7;pQL9L&7X8eIy1Nl6{UvFbu=Aj_ot#vGv9X92zG;&BQ8<+frz>P{XPQCC
z#&<?WG2S*WvKv+?0lM%00OA%#w^{Rrke?7#)_y}f&Y0*UM^ezCMnqG43(Pp_Op%Zl
zi`LlKGxbx^HYIHFX194^|BL*){fbkm;2c05+jY!nv5JFH;($3eUe?ye<!>m8Gg)+#
z3vICe>76q*a@Kei*}o^a@R*~O5U!oyk!EQq7Ypuj9~r3R^zrqHf5!9!U%fdVo>eq1
z_Z1{RkHR^l@*v4qNDJiQ%!u4kgK27ZjNhd5Gx-~ACe^QR{l5TBF94b@Oe8=)ql}1p
zd6fcNg!m9(!IB3Z!&6i?(~pge$p{5ee}7T%9pC0R4w;=fV7@T$e>Cq?D}!+|0O$A=
zD0F5~8sh{-g3mdQ=u#A<g)piI<J7sHH~tN0Zk%gfg=15E@KOQOSZKH@PkwNO(YtY5
zD~a~q8p(mOIy$G6Wf31e6cW-pJUpB*E9>`z0)oW@nT79WssL0z&h(uoMraItX+1Xc
z{faq)D;K6SL4XLNV@w+a5g128RLhR%y%cN%XW*o_<dzq8Azj_c!8&%64y!4|XlbZB
z6CTC$kMeSvSTxjG0sN1I9fa#0asu~EI|1xp6@JnM3&<JNkI>?EDuhNk@@$^D%qW<*
zX|@MO@lwWeoF?KojT8E_NRK@+b0Zgxx-gT=uJ)Ui7KQ1aTxV(*oxc0pj_dDJsMKwk
zRey&cSueCBC+-cZjJBd*2=|C{!MgT-28MTpX{ELV8zW$*T$doHs@gSRZ#zcYl1CBc
z_3t@{6csIeMusn}p{L&>yHmX|jB-5~1W!6XJ4JyD)ey@QRq?uRZ$Au=E=u!r{hz0K
zF#VTZv+}7bhMsg6h2c)qtE0R{tg~{{Ae&RK1Jz~%XcZa0bVfWVeN(9PV++6gkRu3-
z`j6pq15CwP=+=B5<;*(5lWIJv-(}&!2}K4|1RNp&z_L==QzPc=Id6D&(p)P%6X1XI
zfFHM<73xnh_ZNAnyqoS9j32nPFYzEe?$@H1prkQEbX`iKl8w_n;WDEj$0f6J+2kT7
z-;;tuI%eTug)ly`>vqI(mD3U0Jc&Fj6tCP7aM=AieS3)7l_+^*;av+?n~sJXJCqJR
zv^o4u?m^0<08cf{nWJ`rPqH>hT^W@yP|8XCcB10rs(`4;Lx;<Iggd809UocKgAOuS
zU9`!fR)_h!l)-tLY8Qb%=$^p3Wlkd%I8%$gY5BH&Uqa&I;wFE5H7UL|JN2`1(?#L?
zqIk@c$KP%-UhF#l_$d-1#>1>G4>FCGV`&(PlzmowI-Z%t^}#i&hpK4rZSF~#Y|J~j
zKm4DM>`yLW=DGHCe_UpyR|S;GAMkafpQ75U{A-52&6`iC65G3cYlbR!?)NbSwaDp&
zyGCnc^_HEEFHkQ1^GPsB*)UHH@`q``-0q_$cydqR%0ug?%0CiU))FwUk`E%)%t;_l
zJs}>STV`zhX|(N>D6bBsDE}f^`NZeN&z&VJ+dZ+m4D=;t3(NysSn!(dcF|sCq4=5<
z&b#<!%PIrB&MO=w-dA5ZU*k@A*F?!As~%&RF3)^HGkYu%^&My%+@F#lRhS-WOoQ>#
zAaRK8iMmBI#?NJXhDv)XX(1s!!FxB}h<GvD_5FE`RD#7y<P&!ylBuj(Wj$p2RYuM@
zP9ss4+kfIT=#7S(pzFpyo8S3|*06dKwdnHc!&v~U{2xT>!t~V7=ha(>e=4P1P~8mF
zmD-L%jqccWI3*rGiV6WyR~7AC{P^itP!$sPA2e=a>zbHcPS^*!RP`U-2FnNaF;&uV
z2)$)J`Gaa?cFfsdY@cD8_3p~gwT#S_zMacnbN%^M>Wno!dF1tZR{7MFmyONWbtWGq
z+^tm~No-Q4k~jVaRId99nIzIT3aYL3jn5<qjy@yqK985)1ltFhTTE>~I<+{}@CV*!
zxVZt^{7)kxVYu4cYYY^_K=HF@%VH?*`TDf4*}g^0v`^(5;Y-HwrOJ_9nn<V7ma3CP
zy|vc^mz@g^BvVx+7)iP;TIp$y1)Mn4?8`9CU?+e1VugtCTsjvwVn72XlFoVZRfYt^
zwD8s^nxg=K#WsB-lK8Hh!;?x(OkQW|bR#rtUpHqjBrka}3@+t5l^A$1PA~2q_Zc96
z+jGiIKb-dRfM;KwE8!TxXK~oKCIS$>%Z!>S2jXEfwm^oDeF{pSXf?+%f1uC94#Q*m
z-o2*#7!cDo=3Vn<mRm7SIQbhtrG!K4@-dvf*vkjQv`d#Nn<#7QXkWh)nkYg-fXmah
zch?K9#@*5d_)u_d1B<rh4CiL8j5Sy^DXCJB%HCk^Wzim9<5;Vh;r?bydC2-Te7t#S
zG&4(%HI{JI!~W524JbwtkE<5fo&K7pYNI(cs@eVSZYw2RGHhEj#TGqNm$t_zIG1@S
z00vdFxQ|;mZ>i`T#)}@2+`K@(K=$6h^y7_P@i4OyI)>M=64S)zD#<>)u(%rKs>8N>
z2w$$d+i2hCY>dETcqt|8kK2)4KWb#ulx{?nC~fh7Be<X>3)WpWktAcA0)koOh%8nL
z!K*4N<(sR6UeJ-58mnlI<5#3S{ZXSCx7~AW*l+(wq)}F`wYcCLU@dm`M}NG#=jHXp
zj=(2=9lwnMaX%cyp%{KOJw4O~SH5Sy_>(2=;lVTWN6V9`u_(K%r+-wPI|JAZDBt$j
zs2qpV7k`#>o70;c9+yl{)<jvZtRL229%qTMGS`3yMI_ivnXT2n`-iQ~InHAJm!r0Q
z9nUcTr1SPYJRs`uSI$O1+>E}XZv^is*zWOJ{Jg3z!Z!XMdbuEE<K?AxXj}Y^owwnU
zQ-Ke`==~SKQ};M0{_z;!6{YD3ng1a2>N4cjfchr6?t>M|bWIzvhf1<mj|@{c4soxr
zkdTa)4YwdOaj-nDR&&vK??bmE_eIOfzGG5$qd^ICci054<5F9LYsu6}`J>1gspT4w
z9}W#SgD*1R(N}Wvfpc>!&EUB`7scVJlI39C-V4QPE9@N3BC7LEdpRaqayn9%%#2m)
znQmQUnoeAwp1!uQ>OLJG$qAr4a-R;<GK=K9{VNdvI{a1YM-9Kyn#_%gsmpO(gl0;r
zHV?^@Cofw6_N|qzOUv7BZ@2cnt*Fm)eTq*QGQjo21}xi!V@v)obwLm#w2fC}hl>H3
z$<J|1jP_}8gt4)^NGr;6rTZ_*+wPK>&3gFIAdbCRa6!n1{d=7H;x6;Y*BImxEREXA
zV^8(8qW}d&A8}Y`#Gq5jR+UfIZrbxu6zdC2BE4`@TbGEzt<^2a9PVx_%2$dQEk+`e
zZk8m>0>jYG`mG0l1E0hm7iGIPrx}Fr)miA>M*dwy_mL*IH#g^4e7vI-Bk$8EM7@yj
zHUP7|$g~6vdelSkvz2i4lUZ%B`MWf9k_+g?Ot<__0Amw;e4?tI_4RLIl6=0@fM|G+
z-|W_S*?cMRQbq;_6qqwnHy1v!o87nEe<DUEX>Y(wzYvkQ8fyoK(EMu+<{loEg;lT5
zmz&PKdH*E;+V01YMf0~XU>g3Rh{_y1MDi=PhMmvMO%R)H2=%@GBg)pyt#bElkZqht
zOLC;yT&d`S3#_kX?P1&1-}Li8vYweZ5z^C(`W`8J(z}ZM*bC(%-lG&ewTMY0O%h@9
zC+3eUCwp4!VwK!sq6!$Z?Y)v1m)@#x^sbg?niY`kLqg^#gnbc0ogkLfaRMF`;DcU@
zW|W-Ww|CanEpINmmu)fXH^Rb`ljg33woXq^tEi}?r>Co%uEGSqI_SZGV7~g(b6NFO
zL?8lTVhVX?Qe&zPyNA#3cOU3I1BIOO#=i(vchx@;!ri29!RSAOgxs3Rr*rIa3Gp-5
z_M~ZLn#6*XM!6Art~ok7nny@p-^eJ4f9oOB@1326Frwm+k-sP4oVhu*rN!{~G+aF2
zG$SL!=|1u9h-=4kjG~B?lx1z~uelL6TPA@L<m;R=GxL%5_4;6(cz{c}u)0qK>~Ow<
zUNP0D`dN~hXCC)vBP+&t1E(?;=t1LQh+8#h5UNDbB&f{Rs9gE-5IFHIzO0sXlR9|p
z6Rb%*V>zbs##qkbCJMb!Hll#!i<iU=tY*m8D)zy|zT|w#5Tc0<*HHx(fex>@sYkQH
z*RNk=6r5|Nk0>ZOw`czxCT7-RcT`DY!k$mfQPxYF^k<MZ-Uzm+<=(Q_D62@+p*4fm
zh-EeYb*NzILx&SFBg2gB%PQ~~LAu^B?WGt{-`sozh-qkq!(VT7=D8vaK3FdG*5p@s
z$|awx2Y<flzyjD_)e+ara9}U_PK9bGLO_!7u>YPzWG9=_uD-11-!cO=NA1xAqg`<|
z5kYYhn@A597FbMZ@k2XcW~Rd#@9C)>BR>XzBL7T@ImwzCCH)Jw^k#P@z2W@~99L}O
zi?BF-PxC^>HzuuY-8*YZ`$*<~J~wFjRvE_%ArZ|kM3f+(%Y(pk*NU-{of#KgK!S@`
zjk&@j^pcZz*TnlwSpr^gXbs$(z-?D|wt=p_%7qg!xD7NsIM62YyMz7e9ikW0Dn@9F
zk=pX<H{e#(@Obu>oytca?SIAx%!T!Ch>oR1G^Ibuck}iJrI&9;OJlg9r}MR^6hV1N
zH!`DWQV(vXIFnUwO`XTqF{T1c52bs{Qe4!RN8Ss4*-Th@v=k$D=I1noHU$fFW^VBo
zn8d#wdH*f-5@a;G`Pz8^lb>{>=_0Ea-q>?lw!8s3!i4nWCTqQ2e}glzn?IW=Yoazp
z=~p~TqaCXJjL?GT41a!&x;as~EONV$x{uFtX@Fa$6`Y(D_UuIg^Z6f1?`3aS^SYTg
zhne0kd`&YWsK_qf^SUs&XEX_}$dFfn$4Wv7y{?ERYs|NZ|0VQmlrczWQAmoBJP83A
z9;;Df{a1F*ZbQ8Eyy6>Q{F6nLDZmT(WQ`zA<m3qoa)W@}&som3z!N<YT;%;aUDzcj
zT}n<);+Le)AV>)IiZ4e0o^V*u#}Qgoye!-|3~XRtpN70DScuSyTWUDBCDZv}UB4^r
z<@PP>!qxNt&U;(l6sbfg$%=N+<4;qHNFA|)8~Zs(Zy%q#%@eW3AB3mk`akeK1C>Q$
zY)zCPOa&no15Tm}9t!gK%TzJah1J`+n8juz;E9I1t0vynMnRlKs)iHEplx^}8HVQ|
z4);z3O~`9$S$+z5$rmI=!j7L<p7|{oek-)(@Ny8!t)0KK79&}Kf=Bp7&-*lhK%z3j
z4_?%znJpPN$niXj{dWKA<HbbmPM)RxdFqbZPiw~1U=3tI=yzgH5b3uVC%kN;qaWp8
zONx)bvi(xvr&vLa#XEd`qSPD@Q2g3TtuukZk#o|bm>XEkHmPAh!pmJ*I!~izOI;Up
zr+x-_&|lkv2IJSWR{O8iJAR(*ebfWLxEFh-(n&@JyItCH9Mtjnih;R{RRbG|q;1nA
z1%)6XCje1a_c<6X@mn&d7GPkV4uAAzJpJ3BxX&Q{+g}Z*SfSBBhfX^ny!|7@o!Xv!
z^D&q-5_#RG66tPk@NEyg78PzrQ3*I_(!rhtD5L?;K^V-<J{p0u6~KJ9`pu-MyD8wb
zl2B&>TM9c*m?1Hg5s(CR6;9y)Ywt?qn!eI*I%>7js)$w*lwm5Or2z{fB0>~b6tN>>
zWR1280tEy_kU%1CRIOV<Hi5WNWf24lVUL1KSriunvIG#>mju}YBqVcg#LoL=e!utI
zd>KBtBoO}h-e)=Y+~+yh*SWLgqiM|H40Q{q91>q?SRb|NuSc<TpF|w{1;IBYAogRf
zD_rGfQ{3Ms4R2q9KZEbr^Z6%~v)x=V8Iwz*e0?pX2mKIOY709vHt#MIY7Z&z`}dN0
z??@CnMf5raDy$u(+5X4dz}8#8Ie}KM86F;`4Vy}1qm?Ts-EtbdV-mgFedg=@5nxlk
zu>|C-tDl|>va{sMN61XcxQ*<XG4jdnw%oynV|rRXH!CsUxA-Gsr`Q_&*8M4pDW^Cn
zPW9qGF3;buZ`tz616Q`dmazz*Oy<vwl)o_9={aRJ@ES-x22Y&~mIoQO-ILq+fNL8*
z%_U`Kc!ChS*@ruDTCZ;fiSjjjC($-Gytxh#NuKJIuslDF!Kvu9`1V9oRQN-^D9Zk+
z(FWQDvhi1;7Q+?uvRdYqwOw*Os(Sube6($5aRe3dC`i7&`PpL1mI#bcfg6W;;ZDbB
zJt+>pnM+j+E&iB_So`am&6>=G8O)u}_iOw_ElsaAU2qa_ni6T<o_dj!&9#a{llC61
zW^l#~9raPls)kyvbMNI$c{t<RZb9&emi`9cn@vZr?r$CS5ST77-<3gk4zReT1<HSa
z>vl!af!iMWfjnM5C;rB(+Z+A4%}Y$~>r)6y9s^tz3T5?PG1tR-qoYkBP|;pt(S9=g
z)Lfr0$Y$owjx?c@<Knd@lH+DL>iHOFDVOX6{lu5e7CopnO=KD<PeNph=1^k(Wz{$A
z^q@|<@R}l2I&P*nK(V{+f?iQhc6N4A5et%qJaC0d{}I1ZPOV7F_xD2(_Hs3L9zU2C
zBoM58cJ;o`Y$qDJ*7fo_Z`o!#PoP=dhRKuwk0m81tEs6yyR4yM-r3z91AL6cP2?mY
zeo1<^iO4vba132v-_cB+Rct(;vc_!`O?z_`TtMtw8?ETh;XO%blonXUuK7?hn^r*!
zf?7Gy?d{8EG3*>Se$t#`3?Pb48l{txl7MLaJCje{Yh^o7)~i}Zzg8)J#oKb(s<-S+
zBz4@qIz^@M&QcX_nP~vOriKyNmCz)YiJ6?}wSL`qEW%^}6(*;T2>fa$q#ZObPfwv^
zVas_P8Bu#A1nC_5^mnl{Df8fiMFHOXu%GAW<-M%4>(;`4$j0+B6ADFsetz!^pa&qx
z^K+L~5=bkV_w@vYqU=uY2q5sb8xr0Ay>Gpt$wP2j_zh8PG%GYC(yg!V!QodYl07_b
z-qZ}CYkmuqOK3qlBREUK)wt0Jk`(l@$y)MdK0_m@f5SvNcj>*fG*MTVtB1#8VCD6$
zh}vVDT<P>5L6&a&;ZCmYF%4JjG0+n}g#E-5bmeEEahn_tm;R<THizc6$GFMYlX56%
z2W9<O0-<xde#&67xqIi+Gm?BxNp>0u{`k2z1J$bW-cME0-LGcgTyc%=WXoW1M4ZI5
zG%t5|9qD*=b+v;ossfB2^B;XRWE*l*w1Dyp&H(P+HyvLIOZ+NcP+vFh?C-6WjHd~m
zvm~hh5s)EK3(d_vUUK#;Y-Kn`K;y1Ye^te-=}23kvl<^fQfp{vNWlA=4Sjuy4fL!|
z42sulwan%yDp0Gcs$z-ULHzOdB%9=HuO?ZW9@|(hwV^xmfyIvxwo6Oe+uN~oatt2k
zZ)D(Oyf;q>K|75ctBos*n>qAU%`2$Kz9|M;o%jBEn<#>k?tgBBPdOI_4pKKdy#~!v
z_I7qlV71=wpOUI)m|8Hz@X<22?lFJE!dp`*>XC`oC8>^EH7=>DtTj=|>767WD+y{f
z{j-T$`Z(2x$A<kEjJ02fOJEUMa+np6bADfT{Y)A?{Xam~NjkJSVE%Lev=hC+Vv_{g
ziMjLmg0MWJH_f%d%}fIe{q(^Joe(gmH&$)d%ulHW8gfivubEHz`;WjiA8h%`r5}^R
z!*k2Q{#bts*6Yj83zGv|x7|zH-khOf@9%YZTeHo*g$RM`$IrLJ-M`veY3b{m0}|Bf
zOc~08v=)Z&_Gh$ReHgLz7v&sx`z2qr)DK+gs2Jk4Jz{=j@xB}t&z%j{>cqR<5t5|U
ztIul3<3pEom$LR!C~K^!)Mr(K!P?+Xj*W+xMv9Qt5x23>Z2m(p%k3-w#(s3L{F@N!
z<vHIRsmWe9cT+)|N}R#q_x4Mc=6A?nY?FH{o4*7~`QdY*(G>+hYhE>Fk!LklP?4?=
zdzTdNzu&)S9SE`;mH>ibWx@+z5A&)B?Xp6Jyc#tABHH5Mz<6EOgp1fgH{^o#vRAA-
zx7B@E)t{Y+8lL)GsBADRW9lqp%6hxNj*9P`)T>i02Fm#So0)@hdCE*x1B{IF%Lyc3
z{0X{|qO7<}1A0MF3s&5vF$2Y&lsTW_w^9C~xJzRypxeK{dB+26eYJK(SVpZpARYYx
z!QLQ<#+6N6+MbO)UA{1{r&viZB{8wI%?EqTwDU##`%PyZ%5rc*z|sru+QAbud31Yw
z<_qJF@1XS{mra~uTq&gTP-w-LK4S+}+qzs_Y;p|c=5$9tesWro-RQL1dY_v;srkB*
z<r7aB(#`L}YU>A-^a`$+RJ-1d+pBm#?V_5|3l5EqN4krXj57<Ia^OLY3B7&eG483Q
z!`_dkEf`{L7|B*Qe+%Ej4nMFCKkQSs!nfmvW2W=YcW%qL(t&b4$jm*-%F4p#8=_;#
zC~U`^m{O|m$-p&cp!$ix!%9?mqqgc8q9m!1?z+p7$t0pkN$|CHXy2KkZj{u5@6gR%
zdinuAyqNz;ye6ElpEufafi|4Di;8{&uq066V|#e?iDW_%yfbu+Z|ONk%A$y(S{i$z
zCwcGvK$MD;O)+u<eyEjhklW@YAMui}28hOMG^hoNg7w4fWu?1VbIk*!AD-37tBJ^4
zBR<Snga^a!_*m^B=se3uxV%Da@3Y)%+9UOId*tJ3ayv#dZqWEjbKTI(jf!8-r!(#+
zp$^vNMK?d0-G53$TFQ}H%s4RI%~!ln>Vu-b|4BEe`J^aV#6;PS<;IOADCe4`zUB(!
zd2596fs>LCc3EyVGu}N}z^I@#pkdX<hL4LF>>K@UlE4hK${#@4u;@{+!YF=igh*j=
zXfi!*mgYKWK97Y)H~3tYifoFYQPZm798f-emipnthXcUQJiVw!U)X>t2s$F?$j3M=
zG9B^p@o^^p=(`+2JAlyNFgahBln}N0<idpT0Y2qofj>hx5&@bad-tv+RBh2Ba2-+X
zs8T9ge2=s|K36zDf%G0nmYYJ~hCG12w&f9;6DiX}HN5R;1=f@WV-4?_`0}5;$p`mo
zWy^DrS)?UbAGG1}inyb%Ey!<@9XQJk!avyKZ#Y7sP<N&f1M)iJ<>e*9|LibD%)mgB
zKCpbqI`Jp<H>ITr^J=bsvlTp$K3gVjbkWakOABh#(Z*-q+!z|ji8l6aKRHclfc>I1
zqNlqX8i;wtF66t>1T5tk<Q}q$S|Wa*1sX-wz<HjT0UA*;c=feM?>Ts_@IEE-UY>mI
zJ(3=psgZ8RLRA}mwr+hEGG#ha4y&BXt*WkW1gkODk(t*<Mh&xJbtm5S$#r<OdoNer
zJPHENImjNDHP~o_DMI!wtaR}!r%iP2BDZYdDAQ=wD#ED|jd6%mH@<lj#(3VH$^^N%
z1NDC2N~K$tabxAXsu%k=q$B@KNW?T8?}F(Qf<B8QhGSx4g5p9Yjw-JI024Gw#;P!;
zaY$Nr_w*={Wc?8-P;2bBH#9TV%iPX@c9|`MAqRpw8~eTVN*9#@Yrh-oQ$`NI!_dGy
z*1@5P8+->VU3|*7UQ*8;#&s-EFC46-@kH?Vpt2Qr8&J_FN2~202TIOLL_9W70!KD<
z7pIrk$&=wC28e@$k3_3U6Jh4am^EbYG|(7PIF>ECB5)ML;`x}NRGw^Hv3$E&17X)!
zd(`2z@D<j}cE%78V(4D&;QeI2!jaJSf|O3r&(9y6Ek>IWqQoV|*OgaP-(st8#cC%G
zyngk{8}psb<q|<9C>`mF>n#M(1m)M&<$=J`<^&C#*GUbKjdh@E(q3qEEdp{)AaCXF
z42fNC)1Uph%_o%tw_a$ad9p~?Jjz9Tccw>I2?W8-4?|Ex$<A$$aW|hq*jrR{4kOe}
zAdXc^c~$VbQK&Q4Am(}^ptiB}y2uSWL?>uhaJRM2-jsWw%($Kch>`}pw$N_lbw&Sy
z=mHI=F{07g95Q{@l69j5&uFx4y|~*kA&PVEz0D-6XzA|mt{aK<*7<8GYoP$Jl_QsS
z5CwWDh0CK>d)SEvRT~GOTxtZjph+@mOpM~$MK!Szd(a5h$|Y^sW{6#+PT-Qd@$2g9
z2zDr9w1<aBC}ak08+;h33gGA;+AeCSOz}JpXV*MUr0=_f+cU~_@f7ZEZZm5IevJ^_
zgx2J~RUNAkOuOkkJ8SDqm;;heAPDfLu^0aBl^h-hC?!VoW~UK};V2bwiS>Sp@CVfA
zORm7Fsp*iwec_PJtODnXsC4*bu?uwE15uw^Q&kn-L)G$B1ipf`CMLZ8GkT=Fm?ryh
z+2whRO-?wUx*#brF&br=daMr@q$Rx`xrvEgsK6wPWh=Mu>+8cV0G_~x9F_Uk+=evt
zQi(^z@WwGuLF4ZujkFG>Fk1`=fbw^%xEM}ry-t`t3nJ_gH?n;@Xo|m}ecC$Z@G>}d
zK>wrZ0Ea}G3OEx~@Xt*}!HosS#}ZOJqg7UYsd3`@SyBhBd0Gzb%nDjX<(@p~(w9k6
zt@iP(4}ISB{)fHg_mY!)208MeojZ4mm^$I5R1LGBoc?$4BpR|5sybK6S2FEo&%IcK
zuhXAiS}Y+V4amN3uC86BEhHS~<>k%P@O=9~^E7(bW-EfWDNoQ2HZ~Aq^)6<<<6U<<
zv?~t!)qo90KfjZ1qy%-p^|v~Z0DuDe(JS8zt=ZuAOwL4asWvB8u4liD>u?YYq<I>!
zpk3qAX{xVus`x#>w)S`<EDD<e9eQd-)01X@R);_Q+gd2Li9x`YI<6bwc{!&4_yi_F
zf_OW9x?-5iX*=5xJE#TcEKm;1G7PgX>dh$kglaAZcAeyE=+Z7)v<PuMP+k&ON?ovG
z(e&t2DqHBwk)U63h)9HyTn-Z42`Lfh8h}y6oEqM|c>|X$)ryQo0R|MWM$z1M3Pe8U
z1E&;i^vt{-q&c1>pQofL9f;(87ibn{|Ds8I_|^qoyTZU>qGz!A!q@YcY{z^&efspe
zXP)bm^YeGnCWd~~+JWsRNb|aUf?<lMisu%fzmbTHA(ObdyD#`|-vbz_5AXJG*!vAn
z+%HF>sFh4pWR14ntG>|j>Uy0}NB~;ezVsFYEXd<?DsbIO<BnK0uY&Z>)usSvS8p{p
zhCQLOal#X2$~@-iD6r6ut9tqJ2x_@e+s*A)%HXyYK#a=>N?dWPTz6vzu#rR-m**Za
z49myN%N}c+eecy+ph)YmMTeuWHmhh`tpmpo=NQHs&%4xD*%(1KpO~vEKHr8;2l^UA
zFj}5xym2JE;*0X(d8?Fg*WI<!P;1-|_y5PV|NG|ui}3$98~?r82$4?-uo$Vl{JHOa
RTzp8e*le{a)qMZSe*t^T5OM$j

literal 0
HcmV?d00001

diff --git a/docs/source/_images/prod.png b/docs/source/_images/prod.png
new file mode 100644
index 0000000000000000000000000000000000000000..a4c2887d8fa4aceaa71351dab6ac73a8fb02b705
GIT binary patch
literal 18894
zcmeHvc{r8Z|Mr7~1{5kqwgzKSh6=?#D2mLYga(mRhRV*o8xd)j5He)QP!T2blrm)=
zGtbJjcjlS*v-Z~cp5O01=e*bT{`Fq(_3K=x9z1(JYprLk`@TQ-=d*ZfX{fAUvvm!I
zVe5~pD(hev!%y@-Gd+CgbH36C_-FM6H5FxS8T}WY9pj5(IPAFcA9}Z64!60Hj4!z;
z{Q5Dr)6AFAZs*3vRdqy#pvOWS1s=!vPreepe7EpSi^dr@w)OlCQL?9USxbj2!gVC~
z?wX0*C--hcT=Q{VkI47=>uOGn8eC`~D&8&sY)%kp(&Mbz)9bpNwky5K?4Z-sl5C|+
zV&g9rvdi@sH|#3f=W{B@9psIK{7%qem};tSkPGxOuY`q+VeEAIyco6vm&SnMSFIvo
zm<O?w6~i{Lx6)xa`nz%%ri`~?!idbw|8w*Im$QM_xHsR~*_nLl+iT6R_U7j1>guG6
zMLy<y8!*f}T1&()$YtH*1BVYEc3WD=&(E)@sAy<tNJvOH;xu)&<<qr7Z|-t=DGWQ5
zR12%><gdrp&KHLo-UT1DHc7KN&)jv=+Od{#^Vwft#a&1HgRQSmy`sbNJo{`5T?T5R
z*^W$TnVNPywMfVqj5l2zG<b4k%8=s;Y_395(bnnq_V%f%sZxvnijd43KPniQSfz$Y
z7<TTe!#eW()QIRNFTYAoudh0BQt+q2)XI|Of#M5^=o}^{3=@0?4Y0`CgtIH&o@&w(
zlA71TFo7NJ(lw`_6?K~I*)&b*Fr6){s?vI;_2R&_W<$4kT8A7zY+ln<^=gjyyd&qD
z^cKdh=i@l{7#=%^1`Qdumv%L|MlJVV<&PPP)I!%<T8w7Prbb#UDnr|9Y9ty=GbbA&
z-NgJ9U)hvNVwjdD^wMjyIC*HOCHa!t?0HjDQ%%jvXS>xj#v1*_+4y;M1wYtb3h`mb
zuw%t-^t#OFZ{NP%p5@>)O&%E@ekptFih+TFo7-~3JENB7=E5Z}@jV(Lb|(&vJ@KCm
z;#MO$bUwea!`))!`FNi2!C}g?)iT+x3=9l!PQ7dXXjz+ok7@HB!}sRpK{7Y`Urt!_
zX1Cq&ZM!kULH@5}Dp(S;?3~+eQ#KO%WH_`S%qHd~lg!~h(zg%xrSODpQ-tI94mYGV
z9=oI?-dKrY_%@GVeZ$PDpJ9`MbXa)k+ZVfR?lQi<HT3xrmqRRk>q3m4fPI1~ePla%
zHkQDIMF^x7X~e1T<&nMlF*nl@5M~hn<?V(hUX1AYOx#K(j2pglEbXnZi#dkHvN>^X
z+|uaJ3SSWQl4Yb(ODxw3=*#CDS}(@vy_O%=XT-u^uPV)W+D3=P@>(W`PN1F~Q!bIM
zwOEZ2&%9URzr%8UpbEZFr3==*_+xIjR{cmq<vOfN@2=DE^HQlt+3^^r_lNUU9O{PW
zGu~}h%3<72`gGV9zWnSik@>Rw><PIh_Doo3)UdW{xVe!mddAj@{<`n!t}Jyk_Knc^
zWnBId`xBRficIuj!Jag{DuP-Yeaabp5)G>G4;xK7G^j(|FsRrb+Mr}V&<1rYw+&6;
z;4zygAv%FlNd{aPQQPC8@XkyI%w~t(`MIBx3|OA+s%gD5Fxgp7Xf#-`uJ7due^%%y
zu{jCaF>lpX&jyHgFS$D)03T#mWN+Pxg9}_*(I3zUVk<}M&ixUUWu<=Qa-0|zl-J+A
zDkyie5&Z02L{OMntPl5etaognHiZ=<ZF^>VRm95Ni~yf|f8iPa(WV!5Wf;a1EmX3P
zzu&4S3T6;r!A!2Mfxh50;mwl@83m*zvAmh6{UxOP+|$Ja&}kBzoyh{4u3|b-0hTDo
zJIBk{Vf$X|333<a4_`fn#znHc>-3O%VSGw`qsK++EjBxBcX1!HC_Hk5`r#dMVqVk>
z9%x1RQ(lz$YAlxh^s}x*{TDu{p~0=&&?<>$n)nv=_U&5rE3is#`>t~jw5_C0e79gh
z-#vNQi%&nm4HMKGa77WC66qHG(weAo5C6_k$vGJ0my$=v>6=cU(}CX2`iAWCH85|z
zaU3=NUG=I%G#a@c8gM>!z@%;5>Jl;+L!~>Yjoj}Hb)q(c%g$eniwQlGZr2wbVA5Pu
zqa(x_db`_mgZLI6p0l5>HcbrF8W|bM$;#Rff7gCf?R+}tFpq3hJ1>UotJzj<Pi}Q+
zG$v&KJSv_sQKR;Rk9Yg_oUu|-Cr8KYIkWw8<h01#+~u)v;(S}S#JWnwpogob>zp~V
zyfEIU@b>fvf&1^R8WT#Z$z(G7=l=5GZ)^lg`w`2UsQNbs1`4tdrIv5h=QXE9E573V
zympyuO}RV{3-^}@HM7HIbG(W1mKn&j?I{eDv=BR+_I|;p;q7U!O?#W)pI4fUoQ;P0
z$9=tU(mXpmdt_w9ud*yej_{)*go%|W)=r8-?Mt-bz4;uM`A1(b^htibkfbNFe*JoV
znMddq2y<H5So?n#3UWCsjhedP@0IpUEk=4C6Nok3L+;4U9SeN*>XqWFeb<^lil`Kw
zDy^h!latS!hpozfue$ru@b*kY3!Rph7AGes(Xb|X%y82OqcdmXdHP~YgSf{*o}A>e
z;xx*-Z3~p*zt^agZGLXf(b3U=&uN`SA(c2ZS&tR3SSg%{t|8=%hp>M>GTmC+)|P#}
zEn{=;%-3tzu35W|2P-ib6cmj2m$zgRtSu~l)Wzv$W@OCG&!=8*Yi?*b)%SDDng<h0
zH|iUOWA2x7<o5Z$zi@QiXEpZ1V?iMBttw^GiJbmJLt$}ZgQkhgWIdt7g)m$B54g<Z
za7j8cvSf0^B(QBQ<0Khl4>>CssBW{kv3Pz&ERz9WHWq5u9b*<~yKOWxunso16lwVA
zMzq$XJh^;V?b(-&hs~v|`ph$gTlZP*&W%`jm$I0EG9cIGJq});$p9<ickb_ve{Vlh
ztq^eHl@@bqdvjJKY)%OPrT5W2s~X;|6+QL+%^GY#u;fd{q*~sHya=EcZgN-PcVWL*
zn-JW_O4504MQ}TKJR`y$z3SL))sKP@4us1CK<4@T*jB?AI>(_SLjlB=BWVBOF43Bx
zB6h5>?IS8;Us3UopdxnchCd!r5j%VntqB&?o(_JciEI%VCSLcrM3SzEQY7l27Gsy^
z>SxdkffIHSNvIbu*poaTBGAXlP)8GK4&Nm{^fkKxnBQrixF32Eeuvfz2Q-i^2emz4
zqb3SfG$NW%6I+=0JvkA$cAC+8A&Yu3sQ$YbWYl8DQ!}8ohg6h4$9Z=H6{TN#=qXcC
z+J@(;haUP=r{R>)u5)TCN(cs}cdtxm`zkb$$EwtUWHiA*5KMa#hgkupIhde}0MkT|
z`>;#!IZcry;$Gg9uK{rM{`ftaT&NE$v_5=CGtv93O6fk}G|p|;_HG2rB*|`8Dl8Xs
zP$y#NE%XuF-xHyU$n_i#bt3p(r{d``8C*AsC0Pz5-5Ohe8W_SGwi*5tGuce_Qdl7S
z=gG!zWDIZqo!XlXXnKzEuMB7VdT6Wj&@*))Qqob>R^ZMmk~V5<Z_T<U1eGNIyK;e_
zV8uC1hr^fb`)PeT1ve9SaVSWVdx@X}aIPEI%_6dT%3)f(9;nkJa9A-Ev4XOijiAM`
ziNMq8$qw%H7>QW^mF2rv2rtYdNdX>QBiG})yf~e^yhw5D6E&5&RQi0aSXyVsP{KAI
zp102Ve<NqlJbvb~^NKmzLU94@;wAy6uCvQal;xZuV~Q-mr`zl|({{(9car91-!pCv
z1zLo1D{k4oJygMMsk8H(g@uKcm6d@(W9Sj5PuE%wTrBhki0X8@^7^gO=2LI?QCRV=
zsuI<r(e*o@9M+k&S)7~Xzu#h=o4xQpGgI#Vdx!D<v4XZtf@YAE<>X*}{SzPl^5wJY
zAU)S~Jy>4zz?^SAR`+xlN9M#S=B&B;#02NtvwZ>X?4izMC1VATw;il_DRZOi^MgM3
z59Z}(K3=YVxPDta-)>X1fOkA?@09<1k4=7y;_VSDQzO^!fRXL(CbM_bZ~P#)rjO^{
zVbJ$^CWk;{zsJrrN=quGhi`f2>t@phA3|17F+%Dw>Hom>3GyCOGC`6K8z7wdz}DFt
z>N3?t$R09A(@n{w<mhT1Nl9nYjeQyPPAj>7r+twH0R)Lh?Pc2s-=kg3N71yl&sGqY
z7Jeu!FHPsVbv@X;-63ZM<o%SClpZLa2Z^fh<8$p}HhuleotE0CB`>z~5-)p-p7tBO
zZ1<}(pO8D>#wI`0iT?Ex&q=e31=>kXO&tpyX-TaNJ(3w69bKwDrUvLox-d|_yArV6
zsbju9G&B^tIA5bC&bvS3_GoubPR?GVPtW<LCMPE+Cia|8j+)aR(*Rz_`PH<zFDvoQ
zRgI)}EIobV-1|s0;m}y`gUu#0bB($3_xC}_W8I#W?KEv~Z?6#|oATz(o7`MCh;4u<
zpzryKdyCiR2D1?NJTtox&^~mKKyKZ@xR?L_P(miTp`bSGUoqb{^T#*gm|`+!=DCYQ
zu2GyW)rqV}E)8g-^LUBYK|`9-I6{%AhvFtIGoeAcplym*9T%`SXYVpYuU!Xk<8B-K
zHjst??b)U*B|6Nm6co&+7h@VQCHrCHTNLsZUQ;foZfu6*+euV>@ip9kJdd73TuYJx
ze;?o=Q2c%`YHh2AjXMiM@XLEZZ5f`sunT0}7bO57ZbIx`giMje+*-T{s~&>bA~Fgb
z0}?0Jq377MUqM<kykZq0b_2S@F`0UY2I`3=LP2l2y63137ga7TE3^)8o%_0rZms&R
zaA@N>;F^UP`=AnP<2AIg4YkpnL~TPr5!x8Q(b@<^ZR{;(*K$K`=!2$%e#A#W8?WIN
zSoE`PJ&M{`b5~A^1h?iXpQa8r<rUmvPNL4pR_d#l@uUkvUB>ZK2OCicGPep3z*08;
z$2i0(Wm%c%^FN(OPvlMO4Svp&!`37C^Olop@o*pjPv?1hFq7zA4{~5O5_jdI?sb%}
z!n}{i>C%f7#LB_X>mI{fl_Zk{2HWUB)-#eM03U63NYF)9dSu7LH7-r*#5-fBiApq-
z12F;@)R@+%TsJmM`V7Ium7^=jA0hzIOWaCD{?#?*@FU!P_~a^pet7+d5LTxD=;7P(
z3Q)E;_{bWT?=^q;!m`10!MlY(C9Ul^Bh2Otc33zkF5wOU0OQ4Gi`|(^$bh8s6n+hO
zJicP)z58Xf*c>iAQKbKLlOzk@-FF@4<81D#R&o@2ci#74K+zKH3b562OwM^q0<d;T
zx8&RSMG@1ek5{vDM?f#k0%*Nh1>Zf#`MU{a^m2W6mVD8OTW9tZsD>mXU%?*<!zkVq
z)<-f%y?8i-b5}=Pp8R?HgH=}ta3|oq?206_ZYq}xUVtxqL($&wFFicddC%HyzQv+E
zNZMhn_tU3OU%fUfEc`0)`+K3SU5rEhwi{UNPH|nFTfgnVS-X<mKSG_z0O&aYxBw4}
zKUEYmX1cv-`tL_29}my)@bJ}Z*G2{hn?G1o0`8laz0d~z=d8KgvA>PZ4m5&^@J7wU
zqsq)(7y0H5XU`J|1m^Rg9sArf36xMN%bLOZH*ffOm6Vsy?u%#lSuy|{;xcxb2VO_@
zo`Brm-QqGMA8R#^OxFD~%0KW%lhXZss&i@hL#6!UMC}v#xu5O|3thsxubqG0<J+Uq
zv4Z%ee?@%7!o8;38xoG3p=?()1R_;oz`PBb@6B@=oV9H@ZIRtCotQq|?vy*7O+nYf
z3c1!DITw-)XG`HP5x<8IAI9@M7&gjjg;nRF(X8;!I&<U`VCllQ+=Xx7Oqa&qBxGDs
zz8F|2>K&1;!pOm&$zb+8?U_E3^cvkMi7wY>+0+8WT<YJLOvY)t47NK_*7Ji~vp7F(
zkYjiGa;%n+dhmgFrPWD(!(Ro|Sb1feEndGwL9s{2W!t0nPmHB*T;02O?~VSl*~vsV
zW7qF1p3B|8GC*-aY1DOypLF(Q^VL);EIaoV+_`fH=&fk8$?fyAUEN>3APRuUOULzl
zV`F1&?FX}xXM-yI1hi_mM5VGV`x%%xEY1Efy;*(KBFtqf!8ntg*5d<>P~4X11B~rY
z|GO^}&%=Aq)Ch0>4hex@nSogH$!DI%j)@)DQrkYkT2tBX(2Gb6qO<my!h^13zF{xr
z2;`@m*R5Ok_Dsr^+SezU*tVa2e=a{r+N!@SP}jG*I<2kbP_6+NuJ-pXE|G8I#bv@*
z(W!;Jm+6n}pu@d$7GWsbT5VAXuSb<BpZAc)e}r4r%@ZLTwYZP@06>sres98FAJ}-@
zJE)s*pDfxB($4D`E`l1y@j9<&5R8QF8wKXCHwu}nu2Q8eQ#vS7eDCjm0B90#Qx!SG
zuJef^5QMN9<M=_k1?3BhkzS(D+1!=e$8XiRJ%qr9c6A#~(C&K%&mN^=aO_<Wf0sGK
zlJU0y04_80CE@GgB32<R5x*Bf+AV=J25E^uU5jE|uouVs&D%5pyNMqW(iif^#%++6
z1Z)jI=KPvI<<x0OuubrfY2=1jg&KMUvM^9tC000}m&z*P_7K9KyTgK51uqMNT(XC?
z4zY^XTKPZj%%Oc2Z%$<<EiT}{lZd;rHdhIufpD!Ngx^7n#*@ZwvFqWEavC#L*`Pb@
z<i+opBdNS`<&OPu$NF8+#-7bphvAOAG!`;m7XfY0McPxufz<~!sVs|QM2-874lqZd
z`QIl%rKY<Xx{?f92N5+{(pe&a^;WcD(m4ijdT0afjzVp?!wTJz#t?^M0iJIb{G?N;
zM~t*P1c2?ecA`7Z(G8%;;y(K-0*)T;AnF7v@7ec<B*ezCK2*lS?}BdNI-JMoMiKS!
z&`4I55r)D6OQ*t_iyiTE@$P(H4@9!PPpDz&9~of%sVdS&at2x^3CjN<zCcuTo%U9u
zo(FLgZP5y_U*pm7N~x$Ua|nLAX|~&sNQom|0YEK8bLA~Z;@a8YR-eg0+AIH^xDG@K
zad<3lwV4+@&Fr~O%qqQye@KCj=TU*Dr3>M_;F=_3>-=J`amGh2p)86E<&CqQjg3uH
z@}=4sO=|$~pFh7RT~v8@?9Lui8(pzKvwX>c>#Z~Q_!iF9L>>3jWMyMBb?Eaua^%Rt
z4(D-RLThzLha2c`wsq_D6xNuz!_qZVZB{6BozCodXPDGI9=5zN5EEuHyKpdP^zl2R
z)CAMqB?i8kJg4??pp<lr<)uZ_h3?I7_!uQFjWkEPO#v);oCfS&o_AZOWQ|wMcQPp)
z>Yn{4VGkdh9lJOOCgiDvtL@cT-j`_EG=MaPrLllq=ay@!Tp|iLm6-V_rlvlnrW(jg
zg7gG4X8R$V@KSA@mA|+Fg%!U#YPxo4_D9%qv5I`M^X!<HSjNwpYM1#bF0Zep7GTlD
z>22M*HQtW1SfGiQSRB$Y%&?^%FlaRslfra*rSv@B$9AFBEifUL{f?5m&lEx2g0R?j
zXV2EX3GA8AS5Gf3Ep5A=*2T9x$`|IgGz+jE5O(<Xug;~#d9b)N<wibUNV?rSoBXG#
z!>%YDd6A=v|9Vkn<~K6F-+Pn6s=4pcH<M-#=FXN&RvsRI#pV)w`#Zv2IZD<C!#(X_
zp77s?6aXN9X>Ky~z3GER>c@|IGHc8b*nAO=&Ub;(XqwMzthdB%p~%e{!a`G^cK-Vy
z2&N||uQeq{FU&y@0CB-?HUH;)rzBTC8pp0SqGDw=l-aQeZW(#3psTA3!nBdL%+wDb
zZVcAN85<ifE-o6SUcY8#^~Ubzx(A9<%V+aSc;g6c{G4m9m=nuD{07u;Py5OOCsVSg
z+rNqDI@ibRqyC$YJquVsFZo&48*auiMcA(q-ChRtIWdZ>0|`uaRAnuI>k#qHa3Ekh
zWh$H`<RG#$&~|LunN+d2hM-%E4}eqv3;a5v_|PO7|C)xqHZNfCQ7@kCCawV~Y1t>q
z(rJkH+zELS=7V5-;38GX;Sm9tp^c3K47pw(c}e_3st>_HB>>g+e2N$%09B?iJVqq|
z;m#hWb3A|n;Z0OKgtb`O#hh3S&%Lp@(I%4t9y1CM9HX}y#{v{`c_W_;#sX?m>R#nG
zA_^tKyPUsc>K@QRDSxUIq~ihchZFhDLU}CECuM{bI^4uTFf;H<Geo-BpkY@0;_OYj
zSfp2DCHeMZ(g^J)&8P#x1VQ1FF8tX>8nhAhjH*bn6muY=Tb(_t!dZbuu9RD_;9cOy
z_>;3ZIx6bN@>1;w232+FhxBhMyMx>&RXVCL#s({aMdSMm7{d9Xjrbm4hH!20AHb9+
z?pj44u8tp;<WNjoa)gayo2<f5XCzD93h;P&<L|c<j*O?R8njWFzlYD=ch&=Xgi&pj
z7sP^4mZR{b4bx?0t!?3>btrZ%OzlH}M5@4L8dMGA4zGaz3H!kdfyv$C%LXA9K-QKJ
zO-7;=LPBk3GA8`tPiI9OQXQx>h0xPjY3Pi3E1k}9$OsgT9mhzgU~zCc1F63lX3aY<
z&w$k+w_K6=t{i^VpY962_WRkF#O&{QNdRj-gsI#FqAD#VK$2{qrt5~_9N4wjn+o)F
zalo$Bc*5KYy8vqqF9uoU4KSBDHyXLBZkCqBZ=~^&GWhhw!=qws?2<5Uy8>#2N_-9`
zJ2jrTA&HFw>Ef7Rj@4m4aYHLj$Vp4){X}wyP`H>gPZs4FdU>NF?jMXD`Y$KoDZ{ew
znwV@hn$Fb(Q~-L6g;gHv@!na<M#J`d_wE6<-=LW|eF3I9*IGVtOPT|au&1V26+Ge+
zFUt)0km*1mUKNV2hrPT_aRxG$wx9paN~hzA)Dn8|>$BSf#Sc$bt1`m;Z+IIaKVBvo
z-7}a36mtl)>w~qM{rM^YSF0nM%e5na*8&{fcHnBPD<^I(tefJOADN~Z*KSe=lOZK#
z(wO4d*z-iewRaX^59L1%D*ocMbK-WQkN8ilp8ilFYxVUJ7s{K0<w4eT?&bgvgM4c+
z*wx=6S2CoD;+wppbu2+}mueHe4N=0OcRHTl?3m-h7!!cXfU|bTXDba%T_!(2SSy}6
zoa{E^?v^(2IymUI-h<MCH*1PEk7H4vN1NJA^ouj9yVoRnnNRh-vv#49%dY~q0KQEt
zs{Y(xRsG{n2H~Y27-A!-<~)g^`^2fQ6H|+=H!68&S2xkLyiCSlF`Xn0a!7(do6-hu
z&e^^KWTZnj#j$5whoNiQKbSeht7z>(WLKWb%gZb7GN5KUozWQL&>yHUTcWu%5Vm9q
zDQa|W0+j5|{?Oh1M|p56<V7Hp9j6;(f9uxd^t5$bM#{T)5RB}fIhv1eZFna)ml>$Y
z`~akCF&DhofoxZ|<!Lv!A0Y(!rI{Yk^AyOGw*&pUkZ*)Jm8hB6^)mj0;ZtXJQ|dYa
zw^l>usJ**0fG2IcOGG61+#hC52V+4jTN*3_Bycv2eHEqa(7mG$qD}_l<%&IP@g}e_
z*}3w0aYc}vQ)Cu3I6%?D4zW;~f<76NYrQBomn3892P4cO93bkAk=GGXV>G7**Hc$U
zWe&+0`n>Z1;88C+S@8&GYiO1cF$)CZX`1-LiU4MCE?JuyY<8eqS9X1Rr1@$IWdshK
z*fGM-s{P|+_^H$rsyVf$_qhXsW=^%z9SYnqP6lNjlpidP^??10|KRY0@tiX(+RGx*
z8Z5Ev7_=y?h(j!LRuOlBc$4xxcoy3TEEJTVxfd$|gS^V7{R-U+Qi|1+r6!6weXIfv
zFFUskH9X)F9oBzf2qZ@1qn9vHb{2R;CoW;JD8O0YHNr)V2S%Fwk<CRs1E>;Sk(ETR
zwG;YV=SyQG8_1yGAN|Q`-Z1l|6M$IHOEE1g@JQ<J9Mpv@wLMd;G$jhFbBDQ<?z_fG
z%z(KGx2N$U@ykTuow`54`vuO$9LGol@9-r^HN*y`YG7Lxp$(tulNfTQh$6RZS5ae#
z7ecfFg|CJ-1i;;)t%z7|7xVD0l{UDc4b|dpyRnlyG8uMt8E+)^!_&{TW*w!k(u1*h
zKcS69Ycq^R>CgG7dnzxfWkDMYi514;UDC1xi-X|YyTWM)hC+xup)7uSk_5Ez(u3AU
z6%N|q{j&`zFnqLhSK2_Pz1}YYX8a?VQ+EfZRk${ImCXD=M9O%elVDwrq-ch4OBhM~
zpED{2L6sK5LwGA5pr7$R4&NAXBFS<jgLVq$@fdW7n8ODKoFa5$mZnv(D9{Wf(r>MT
z9|R8X5rDlh5-kyY;F{x*Q3vce$EhX@Q^$d5qfb|5B%}@<O1s*WJWv~x`stJX&6~sD
zlOo5&LE+)0VDjvC)|+PKjQJ&=bs9+l&*Q7Ufw44X3IO@D`}SJT)+#L4vVZ>d>(|2U
zI1>|-oSYoUR!H?tL0a)QJIJo?(52HU&5}cOm92ATq&WpCO+N67G!Hemwec!`&vLji
z+eL266g(O9(Q$HcP^E~}?<PjfTHEzv`FeVE6Vs|V*vDeoYL5u5wneX$+~wz+Dj%}5
z@*X54|5L6NQ;tinRA=s*8v-{(LI9G=$nH71+4v(E3gbX;b-J1^kg!;a7I%redD|rh
zjB&8X$0ngjrVTuxI}A+mJdC&#-vZoW4>e@Euf!2`CXV;!!M>mT8TrI@>O<IKnT7m#
zHUvy;%fHwR5=~rcPagDBnI35|*t1F0)O61#wdfisknxMezFe8LbJD-FP7X+RU%Zk#
z>k-hr%L|pu&gIgLP%Du&Rw@XI?zN^fp9SSdJyD9_w2FT0%WrXzB4EJn26rdJAmwG-
zi2iPiV_u>#AG#aF9jm<*ncXH-+-~*;iYBC8fUpSy0c%C(uNQ3L+AooB>sArYHQ<nd
z{=OPH!u4XtfnMuIe_P2x+ud(!=1oqCDJ5o=l)B0HJ=ZObgc<Nl?D$Jwm=U*D#G*+p
zT1yCICOk5KF@)l5!Q6!;|5FR7`MC)M6AH1*?`VN{)<xaL#RW)I&3})u!dr7rqW6|0
z1AWkqb^{{?rXqfu_MVxj4O^Ocxb6Dwrf(+pkgZ=F%3T@)v&ViUrSakanO;&Du~Eq7
zcJ71S5trkdyEtO4fc&xfmX7UPw>l!XZ1joSg81^7_!Do|fAfgwB6ejeF+-9t3}sw-
zaUl0pB}Ub(AIo8N-mt?KvTo3kkO}jYqm`~0vUbE5!naP~FKxkZiX{jp0?q{5GT|%$
zfo;dGGvb$u$@u$O`gA0<WI6V26kV*>4;s;1MZoN?Du)I7i3SsdFv0UTZ0V;iFqrwu
zOa7*v-gb{5ZU|099IWYy7g7BP&X1<!5c7~MkrWidsLUe_&JJQ8;zycAf|y5X6Tk)-
zF;wPJ<)bKrPY#_HrdcEyFwetvC`GUYWl}K^1#E>XJa4GX!11m;B)6`NT>)kjv@bPs
z2`Zmw1$czc>(UYc*cEe+005Q`j2Oy-;C%tRNMK&zC27_K2DJhF2T(WD81n7|?nwC)
zxUm;P(9hMheqtJsW`Byzi|ZH`5309{9eSY$ayop8)(c`%Y#@IdZBB_`%dL(TqzGeL
zDC>4>WhCBSLEO{?c_7p|p{Pn0YAnN90Zp%zQ-v|6daQh%C30E3;T0l~9J$RzM|1>%
z6H7CD=+E7S8KsA=I#Za4$gqkHf@;n65FKgp(gZ>HnJea(x6`JBN;U$e^NbLS$xG(l
zUm316Oer--ZX}OM75uhhUV78W9RE0w5QB2&usVqS-1{jB?Yvt;4l-kxMIhS+^-<vV
zl>;jm{y|N5V6i(V$#sD=1wvvSY`e4HK^H?|%bjWm6^tmw4I9XwiFMo_YCurrq5!rJ
zk2x)@tHP}?dh{oFAT?n@cNykAR+z>exOV{LBDdFKL{8tfI|U014EPLCL&(2^I#qc|
zVtJhB_Zm8U?N3Ri`v0^b0WI8x9z{!F>}*c@{6#hM3~m5aoIW=|@A5eOU4pnY*Q8ts
ziwYQKb?|HF@?D-wTy`8Q5qxskd115&$;Wn^g7smG!xsPN{J@@lbfirx$xd351z?s|
ziyrL;$hi4R>wEeQi+{}#U}mjQoWt%jOYt`fa#4>8fo#X2oZGyzH-B;^-u{YAca+Qw
zJDqWJK#A>0U1KAmSK7e%R(+!D<eQ@HPH&NZ4PAqO1k<4ii=<n~J9Bwye))&n(hnDk
z6VzEcIyphEK-BcVI7XO9!C2p58^1yK3=4%!4wbz%7`9Ln#^YcPX!_rh!)vhpPuutp
zYFJIty=<-3R9LiLp0Wm;A9AACrjqU}_L+yT*%XRh4*G2wNlVTiICGZYe#5G12pNpP
z<x#&BqR#y4?`nPo!Wl@%sP2#qd}hPXi<|TBo*@vXznM%WLmYKVekE-9OWCDPpAe>O
ziA*+Bc*%Y}I|ULQbtd4oLD1wjQ@|$%QsE2X3(XGVwUns+OBR6cH$@Ima_lZa5l?LP
zM*?F;7odCAXdz7K@ZGgyvCgjlR+owVe}NO2a;IvxQ`OFA7ePg-O-_pXoq&{4sI*lB
z1KVMng`atOdGY#h5amE_8$vzE2AG3!CC$zIxP_N^sF?iH<G1p`s!44!;Mrt;Y=oL@
z4(-6YNrJojuZ!z2Zm6v=^P+{}?19vZ%PZzWziT|0a3==aSonbgbVV4j(}V?zKI4vF
zmCcwP-rjv?Hl2~^3}Sf45BE)&9XK)`fi&%hpO^?N*$>-|>M~&xieZ6gX$B8A28E3t
zBtAWGdz*MK^u%ytV-RV3JmS;FilMRwx1E+Ez__C-)vy1l9!aO?%6g&91MNOyc04<s
zf<K05{=?q=<G3UX#smIy)C)dUEF7}Rq*E&`q=#3kANx~1l0d6dT`OoI&l9*?xW}Cx
z+wx$T)v{X<F4{mTPQ@#Z1B~ar__J2i{WZ-D!HS{SMdEku3@`C)t%ZOcYz7eb;JYV^
zNcn+fb*|NnaCy6eHHJr_03rmu;_svd{vNID1+#gc9iN1#r;ry?sFUO>L)Su7E!y_#
z&H|AR(<imP(q6bpht&hzZv|T(i$Cx?X;Hj~mb1VLBP!K1RtT=Hm`apdx~>tB;EsD6
z__XcSYZ`rEg-6TPMc}^iBU`BC5Y7q^hsfd18JY!;$U!{_W*Jc>A_pxwS69fj136sz
zLZc5v4pC_4I~?h#<PZ;ZOC<-rBHA5@9D-;b57b3L<DO%+^i*=##}9>=@AFn+Kn|}_
z3^MDp7EC=P$$?;`JmKO#ZRO!d;sPx>X#_)j4UaD@s$32EqPeoXLn{PR_cSO_7Ol%|
zdodC#{NM!)1}6>At6pJbsk=debBIr2YE>cZ=~Q0Wg<%80tUSNF`l_z1?l~HaDIM`6
zJVpe)q(^}M?qG&ngSlTFar6h+LZ>Z4q@gq&ImUK3|JL@0&r@v=sIycIvTAi(=;1>;
zGv{Z7m;6-Hlarx`(=#*Y6_$Dw8WK&ZV99?U+~C>qA8f3wQ2YVuWpD{bz<Y*zcPJ2w
z1vglv01LMi`2l$fP1v~*8Z3@<EM0r0wPni|Aq8j9mB@th@$nfe_?KEtI(*?fV<QpH
z6Z~1fOkY(EGdy|%`+h>~04nGN<72SNBFxy;aB2KifobBw|IJTY@sKdpozjxMCXSRD
zHUK=`*XpN}E_{WE1cDlc<=IMwvErR5VEKjwGAi6m41@e5pd%a4AuT{plF3{d>Egh8
zvfD2vw<<`GTGqZEuT*dY*W_UP?Vb=qyF<Oc!QY1k`1k=yC6nlbY2cqsUbA*;fik57
zSHJ&iqqt5Vp*U2UE_^kW$_|g-_o#n67uM?$6zIFbAbIekIbUAv5}C4;5FOn|&ZNlq
zcuOiWgM9<y8x(b$bmLhLMWIbz^V5_>=f&S8`e*ziFN!vKt;Z~FoT2`@|2?fJbi!3>
zu}|Dh267;FSHztPwzdC)@W@Y+0Cl3)*4BMPsn#9xPE*7F{{A8=bi_k?4P`F~f=o3|
zUbz#|;%%R<8TCD^MiXQa3Z@aN`!f7b38~V54F@ofdqLcG0@*GzXaOxXg)KH2wU7B*
zC1y`~OIFM@q~_d${=3Z8hLMMd*YU4oW2=ON5R~@pkdP4F7!9Z-7^F}-GHzKxs0pL_
zWwjQ=J=T!ft>=pGuv0Tkif~ZzH36*Baaly6Nu}G;c<w@<INSEax1oqi%rs{KrJ01?
zbJZa?^>^9h8@-fMlCoP<<*+EdOZL-F`YON6_E@oj-I-zD2$#=ccmeO3Hub@_5}s)(
zjKl>Cz;b2#Z_tvITmzds$CDmwics0b-FIp=CZS?4=+#<;;WstCInkfyNk{ZXMd@KE
zw7PF(eXYG3pF(p*Fpsi~Tj#&7ank?xryN6t9QWCG1<ru4jU`u@=u6~3UW@e-hcy;j
zQ!u<dh?U;fWDI{DF3Ns4W8yPk_9vy_kM$;3QoAmJ;g`O@IdQ-BHHH@}A%7`60>zLJ
zrzn424;PRY`BFeF9)aTUQc&lA(pvFb$Z62ZXs|pJ{T-dJS7FTj#pI|Exs6!oJ{NP<
z{i!=4VG0+bG;ap8<BxPP|IZ*)BpNM_QT=_4$dP2=_uO_ZJtSl81y8oRVOV6X#lt$A
zaXM;bQj@>ZhnJ;rs^C9;NGYa{3DNUD`KkCCz3YcCsL^b*+)mg=C8DZb*QGJy{fWR=
z+$;0Y$+%$(^_yl8VLlrm7KH}ER>1^HP4u@Gk1#@6Y*WZFTWdN@53&qfKR<@=Ky>lj
zG{d&&5B_B;F~O(Z@{%ibSO<n#-Q!>usJk*F3>Ji}0{d&OOpEFIPK+p9W%2OL2Wol&
zBd*eEpbnlIXkw~+pnaGb^gd$+o%G{2q_Sf8kP>p;-FYYoUZQ4B{CWZNi4S9cXHG(t
zSB4ZWa&klxS(~WT7^n>IJBQvhnw2WUsdx$g>zwd8-?rlw#`IX#I#`>CfvR@F#eJG&
z!Ribd{yL;v)ey|B3`r_DJKp<V3q6KoG$Uvf9-*TW&^`f(T5!@(4I@k0J`vcpWkq*X
zF;6zV5dnLqUXsjq_*dSg0OJE-5nV`{X`jxSE)fiS%ftJ3h5OnV(Oid8;TbjU)d)g-
z!ap^uuiHy}DL8jFxC3e)5pZ%_yc`9K<!wreG_{uSACFI=TEb~PEZp|W5sFxea^mw-
z$jAnB#1%XNe?C-3DbBg(VB1aIe-0jkmeXO$|7H@czGbLQGIq?p4Y}T<UYn4)AU{@y
zWHcl_7N*^np-ky7$FHiPlxNT&cdq{8lOw3W-C2_n+xpHf@d|U7zY(14;aF5SE@7)?
zT~Y_9icnsq)IvZrGz+}@u`obc6j7+aNr06zT98!1NIO8oQ414amDFr&H60IXJPhSo
z?yLbf)~{PP+?+B5)CP5>5Ka9}ngYVT)2QyKw8p6`F_WBOx;)n?YLGBYp}0XX7O>ao
z6zgb=)TH@d?+8!#RhU_aM@PAoInPFeWDziV5JWzPH&Db&$O44V{`Gk+#3>VR5_2*P
z_xKM%>bo0q)S{*eQc`+CoUcy?!;rrpYZ292oP`}7oosGiHY}E4rtf&qiMdNMY__aD
z+bQ*0O9Mn((`^KWxu211H|-H(OG!<go|-z9a~lQ#TK~GGWt<%+z<FI+uSs(vK~3$O
zmFavloTacQP(S3h*a$^TuA|*wfm0wT5mp!;8X6iHIGt>Ed=m0P1JV)ZuUE?RU%jen
z;ra>Xv$K0sTgM(A%<6istKsk~pDA|<N}(XCSt>Qi(b>)m1<D7l+o6BpUdCvIG}qVH
zcXZ_58fpN^0tL$9o*S_$XKIKo0HsN=WuvF(z{4DZ(udAUZ~)OeK*3{5qY0R#ji<<O
zK))EKT&agyVcmW>?Ogs{2w}<sC1V64BW<B}+Vbhbbo`0m<~t^=(<Efh(_{>bt=S3=
zj2=p|oet#RMys2|aJLGLsXP6kB!+51e7AvVg75uL79(~i&9lI2hRmeV4lYZR^e<fs
zmie#e?#A{Xt1J8dCP*eh2*bnqAzwsR#a;Bix)lsDFF9lmTfvz%joB~9;sQ&$6buSy
zwy%a-bci9I18|UJ;NmGrSG%;lUUXFMhcatrv+$2Qy7Hc2R;`KFX<tg0C<Mm^f9x@M
zICiLy_HY#_5$)ivm|~Zjr}voWfEOc{9L1@gtuhSa1+VU00telsZ(X_F&n-#TMv@%c
zphf-r|K{AR6FuZ7zxnSCmUnT46s3wvkq@t2kc+Y>6b9<T!q0-kKycD!&4bSznuqsy
zGhPh4C>>dD_b(;ze@to2<}>j*qm6IX@V3NPFf}D*uSwQSmpfmAGyF8YrY2f#?Ss$r
z0Cc5-9|)zev`xDVK~QV-txYDBs$6aUH~>jL=ni$O4nzlBZ%Xbf)|#rA3$jWO=K1kV
zO-yF&?r?}5^LgkxoD-5oYFu{cE$HSBhaVw+k%5Q%n6}^UO#UVoEOF^XlA$wL5c9ph
znx}ls$DRj<Z|u1N!Ej?X`4)HwsU5EKp!{xCAKe@Rc|Xu=mQ%L7$MYE2Uc#3v*-0hl
zO#2ZCxaWD^a9@k~CYCu2yI~Rb`uDcBMNduV3#HG6!6y-=0;h_kLn#6dFnbAZU29E$
zlF2!c)ij6VaX1hK9&SuuwR-5-22`#t?=m-W+*e?+-gL44sUUn13L~wVyVM|?oXXmg
zX4~@=ZnbL)05Qh?844(q;8<m%bDt;@tEhNoOiQ40PGNamAy>cab@25y{H<UXQJXJd
zQ#cHdgB`qns4;Q50uCNQ(f!W*@4?9*yDXqNtD-Ph4aG#RU`#nfqI0IXF~`2$;fFB@
zK0vLT4pP1`W%pUgJ@tQ=ZRp>wuY-ynXW8H-%gCoaaq+45K2ZJN4<YRHjXA97VBAM5
zQhlg73lEfsI;5Pbl*E&ebe2f77<>LLe-WL6VnOo;`QpX|V_!kdB`At@{9OGu^j5T{
z+gvRa?*CY}rjQjsyya8yF`ceVkqK+wx?9NT+4i@sQQE5q&Yq6bWjirKm&{tTwRcyX
zx{y)7kI3h-aN(o1^c#|7w?2$+y0e|@a@P*wy?dRv{1mvzHqw=|%(OW-cs^J0>YIk?
zxoi`a%mfN+0$cs<R5&p)sOKD5Tmn<48LhKf?5u}4)9TeH)YM>E_Iy*;NO<>dudr~5
ztI7Zrc*W@r0*-l!bI1gq=O7zc&5ref<}T0v0GVUR5j|-b_pai5IPZ7Rr+1=b=hEy&
z_~SgAyDvR8^`N-8HJsTPy>wSF5h{MbnHjk0sG25dXLq0#QbIKY6Qo*YjS4dfIL|ga
zI}1@%;7e;Iv)#9DHNfeR)2C0LJb7|TV~3!^{O|_{2L}kWn&h~+xYlkK^H$L4?Tv_t
z;1WI6+b%3B8hpgbe!+H!;Mrq=$KYtC75cpZl2=W7TvN^%8Cg`k40MW#i0FnxXV6CE
zq^euGs5MW?KyJ~~ohlJR@o(RHi*p=KkBN!-@L~LmCkLc*>R@%AN;Ki^&>ZYd%w2K>
zHi9TsEm}v(=b5a`qX-B&Mth1hU+g;%^C74iTGkFX7PzxjLP-L=V>m_uhgkx_nSv8@
z5`pT@L3Q<~-636bqhFm&vEF58x3Jp^C?j`aOwGusWe6UzINbr0Bqo2dw-*Yl?%uuI
zDhKDqkGL+}Mgx<)`g!q+amLNKw6sIh8n(8!0=|~k)_m_&=Dry_^r&hq?-n!8nE1R_
zTu)#Bt9X$2#v=~-cNu*JPKt<&L)2S|j%AyeZ12$Q?S<3k85tRffP$olrei!`qMejt
z)d;l?nNS!E6%>b{-V8!c2xd2j3E0_j9tF?sh;8>5;QHpOs`qZoOQ?nbir}7@PVw$0
z+-iV3i=OZ;S1GYwqU^F8D~F;Gz#tGZetv$)2S@?}Bey~J$5WF-COs=_0a`i{Q#;xR
z_M_|EKn#?T!5^VlzTBU(s$H=e_X#*uVw#in=FJcsvw@l99x0k_f%M_K?CeG0O;qd$
zC%&gB9SQ(x)vmF*-;AuY>l+$AL9JP}E6J?_$_<~oSUsI-ET_QvB`78|7d`N_X5c0q
zJ#^=DodEz;*8XVI^%>RJA*6%DnBB{f5MC}IBz<eR|NQ6s)vo@zjW86n(xLg3xVX6R
z@N+N~z#V1ZzIiRb64Y^VS%hNrFbHv>>JUj1NZ{FiS^lArlUb^@6q0u4>({TKJJSix
zEiJlndfM96mygt#l1sO#*vSsM&X1&C`sU9oXOBwIgZE!NHsyL((BRamhPt}LQ{OL|
zI*)#x8SO!t*%{X{0WcV@L)H#L-|bLD|MRCu*!#r9#K_2t=g#3yboX|=iHn1+@GXa8
zp(7aZ{JFCB%?_98OB6V+9v&Il17~El?F@>h$ZaicZAW@tG&p2t;jq%4J$tTNTBd`a
z4V#|;$7gkQbxlqA;JuH5wYqiNw&39481+C2cnMYs>E<ABVML=t78_UK0BWnD!XusF
z!>8a1c}}iPn|vB}F7_p4+3t#oe3O_+0lqwJ*INu&d`4g2&*^zxN&1?#Ypv7!{NX+D
z=cu>$LOCf^_$$oz`hHyosg(jJ^}D;f;nxx<AAh;vzIyfQ<U|VWu9^Cbrp+l93O>2}
z1MBLOpm+zufybkAv*x?)#-OCpKwJCid=LfBhb=;`*_KFMxqvL2#Z4BbTCCj`CmX@6
z;1JZ%GcYJpIUZ8BZ;>Gv3tk7a?-B76stmowd#3dqVPzkij(a%5)}h}|eNB&`><5Ws
zL&Fy?pU%Q|*|(46<gm0bi%uj#`Tb*Y--7@G-`tmS;Ws>po%#5)>Xi~av}~8Hoq+EE
zg;JRfGwd00xu@{)=%pQk7E5wBU=cORL0UN*O!02Qi{pcX>H;?)1+FJD3Fna6(!+{N
zz#u}WWZ{r<GJy<x)00E6(6nizjDtsif4|Tpz@yEkOXBMNgDlf;>KhvkPM<~x_TXnT
ziCfOXU}4Ylj{KOr0gHaY;Ba?X>H%qKD8roo)m1P$fYACT=U2*W9Rq`ac@=e^?MJqv
zSALgqEudaDNL`rT>&+kIK7N!^P*9kEAotYQeao0V$@68_<?0BOl!o(N=(J$t+ta>b
zVWwx!z$xBM2rUGBkH^%G|ESEJYPL$f+B7y=#D|WMXmi5$-79|*jdZ+JG*MZ50;;iC
z(MU%-j)#EWh9d%S0wWSC`*%&doq*N(vpN!f>%x^QSDGRF3^#|Khz8hA2E;0R&oZsh
z1r-KxX5%vKDFEwDVrRcMH_J9_PQ*ar(ehG^7^F|iT|ZMt`f;6u;9Elo5T-bu9ZZjl
zvxEbDyLazqWMqW3X9~y4NeZ>TStUyASf<rC3g*uC2Sei7NLyRm#KgqX($d%%rT^^|
z_#J#N28f#`BqoL(cCd4Gbv=6YE2u$|Q&4d|FdH1cgXQaXJE*t>a1_p~hzbh_9<sfl
zuV3#($qtzRv{(}C3?;Cj<TeA1adH9x!=AM)=Ce?3H+#6(uf!n<hhe`cy#c!LTVk-!
t@Ec_4ap)Ju{Lkh8%)|eii{O`A7C|@Rwt#K`{OTj@_)!hz1jTcA{}=l28v6hM

literal 0
HcmV?d00001

diff --git a/docs/source/methods/cmfd.rst b/docs/source/methods/cmfd.rst
index 09f1328bf..67967e933 100644
--- a/docs/source/methods/cmfd.rst
+++ b/docs/source/methods/cmfd.rst
@@ -195,10 +195,361 @@ however, it is very simple and usually fits into MC tally frameworks
 easily. Different methods to calculate more accurate diffusion coefficients can
 found in [Herman]_.
 
+CMFD Equations
+--------------
+
+The first part of this section is devoted to discussing second-order finite
+volume discretization of multigroup diffusion equations. This will be followed
+up by the formulation of CMFD equations that are used in this NDA
+scheme. When performing second-order finite volume discretization of the
+diffusion equation, we need information that relates current to flux. In this
+numerical scheme, each cell is coupled only to its direct neighbors. Therefore,
+only two types of coupling exist: (1) cell-to-cell coupling and (2)
+cell-to-boundary coupling. The derivation of this procedure is referred to as
+finite difference diffusion equations and can be found in literature such
+as [Hebert]_. These current/flux relationships are as follows:
+
+* cell-to-cell coupling
+
+.. math::
+   :label: eq_cell_cell
+
+   \overline{J}^{u,g}_{l\pm1/2,m,n} = -\frac{2\overline{\overline
+   D}_{l\pm1,m,n}^g\overline{\overline
+   D}_{l,m,n}^g}{\overline{\overline D}_{l\pm1,m,n}^g\Delta_l^u +
+   \overline{\overline
+   D}_{l,m,n}^g\Delta_{l\pm1}^u}
+   \left(\pm\overline{\overline{\phi}}_{l\pm1,m,n}^g\mp
+   \overline{\overline{\phi}}_{l,m,n}^g\right),
+
+* cell-to-boundary coupling
+
+.. math::
+   :label: eq_cell_bound
+
+   \overline{J}^{u,g}_{l\pm1/2,m,n} = \pm\frac{2\overline{\overline
+   D}_{l,m,n}^g\left(1 -
+   \beta_{l\pm1/2,m,n}^{u,g}\right)}{4\overline{\overline
+   D}_{l,m,n}^g\left(1 + \beta_{l\pm1/2,m,n}^{u,g}\right) + \left(1 -
+   \beta_{l\pm1/2,m,n}^{u,g}\right)\Delta_l^u}\overline{\overline{\phi}}_{l,m,n}^{g}.
+
+In Eqs. :eq:`eq_cell_cell` and :eq:`eq_cell_bound`, the :math:`\pm` refers to
+left (-:math:`x`) or right (+:math:`x`) surface in the :math:`x` direction,
+back (-:math:`y`) or front (+:math:`y`) surface in the :math:`y` direction and
+bottom (-:math:`z`) or top (+:math:`z`) surface in the :math:`z` direction. For
+cell-to-boundary coupling, a general albedo, :math:`\beta_{l\pm1/2,m,n}^{u,g}`,
+is used. The albedo is defined as the ratio of incoming (:math:`-` superscript)
+to outgoing (:math:`+` superscript) partial current on any surface represented
+as
+
+.. math::
+   :label: eq_albedo
+
+   \beta_{l\pm1/2,m,n}^{u,g} =
+   \frac{\overline{J}^{u,g-}_{l\pm1/2,m,n}}{\overline{J}^{u,g+}_{l\pm1/2,m,n}}.
+
+Common boundary conditions are: vacuum (:math:`\beta=0`), reflective
+(:math:`\beta=1`) and zero flux (:math:`\beta=-1`). Both eq. :eq:`eq_cell_cell`
+and eq. :eq:`eq_cell_bound` can be written in this generic form,
+
+.. math::
+   :label: eq_dtilde
+
+   \overline{J}^{u,g}_{l\pm1/2,m,n} = \widetilde{D}_{l,m,n}^{u,g} \left(\dots\right).
+
+The parameter :math:`\widetilde{D}_{l,m,n}^{u,g}` represents the linear
+coupling term between current and flux. These current relationships can be
+sustituted into eq. :eq:`eq_neut_bal` to produce a linear system of multigroup
+diffusion equations for each spatial cell and energy group. However, a solution
+to these equations is not consistent with a higher order transport solution
+unless equivalence factors are present. This is because both the diffusion
+approximation, governed by Fick's Law, and spatial trunction error will produce
+differences. Therefore, a nonlinear parameter,
+:math:`\widehat{D}_{l,m,n}^{u,g}`, is added to eqs. :eq:`eq_cell_cell` and
+:eq:`eq_cell_bound`. These equations are, respectively,
+
+.. math::
+   :label: eq_dhat_cell
+
+   \overline{J}^{u,g}_{l\pm1/2,m,n} = -\widetilde{D}_{l,m,n}^{u,g}
+   \left(\pm\overline{\overline{\phi}}_{l\pm1,m,n}^g\mp
+   \overline{\overline{\phi}}_{l,m,n}^g\right) + \widehat{D}_{l,m,n}^{u,g}
+   \left(\overline{\overline{\phi}}_{l\pm1,m,n}^g +
+   \overline{\overline{\phi}}_{l,m,n}^g\right)
+
+and
+
+.. math::
+   :label: eq_dhat_bound
+
+   \overline{J}^{u,g}_{l\pm1/2,m,n} = \pm\widetilde{D}_{l,m,n}^{u,g}
+   \overline{\overline{\phi}}_{l,m,n}^{g} + \widehat{D}_{l,m,n}^{u,g}
+   \overline{\overline{\phi}}_{l,m,n}^{g}.
+
+The only unknown in each of these equations is the equivalence parameter. The
+current, linear coupling term and flux can either be obtained or derived from
+MC tallies. Thus, it is called nonlinear because it is dependent on the flux
+which is updated on the next iteration.
+
+Equations :eq:`eq_dhat_cell` and :eq:`eq_dhat_bound` can be substituted into
+eq. :eq:`eq_neut_bal` to create a linear system of equations that is consistent
+with transport physics.  One example of this equation is written for an
+interior cell,
+
+.. math::
+   :label: eq_cmfd_sys
+
+   \sum_{u\in
+   x,y,x}\frac{1}{\Delta_l^u}\left[\left(-\tilde{D}_{l-1/2,m,n}^{u,g} -
+   \hat{D}_{l-1/2,m,n}^{u,g}\right)\overline{\overline{\phi}}_{l-1,m,n}^g\right.
+   + \left(\tilde{D}_{l-1/2,m,n}^{u,g} +
+   \tilde{D}_{l+1/2,m,n}^{u,g} - \hat{D}_{l-1/2,m,n}^{u,g} +
+   \hat{D}_{l+1/2,m,n}^{u,g}\right)\overline{\overline{\phi}}_{l,m,n}^g
+   \\ + 
+   \left. \left(-\tilde{D}_{l+1/2,m,n}^{u,g} +
+   \hat{D}_{l+1/2,m,n}^{u,g}\right)\overline{\overline{\phi}}_{l+1,m,n}^g
+   \right] +
+   \overline{\overline\Sigma}_{t_{l,m,n}}^g\overline{\overline{\phi}}_{l,m,n}^g
+   - \sum\limits_{h=1}^G\overline{\overline{\nu_s\Sigma}}^{h\rightarrow
+   g}_{s_{l,m,n}}\overline{\overline{\phi}}_{l,m,n}^h =
+   \frac{1}{k}\sum\limits_{h=1}^G\overline{\overline{\nu_f\Sigma}}^{h\rightarrow
+   g}_{f_{l,m,n}}\overline{\overline{\phi}}_{l,m,n}^h.
+
+It should be noted that before substitution, eq. :eq:`eq_neut_bal` was divided
+by the volume of the cell, :math:`\Delta_l^u\Delta_m^v\Delta_n^w`. Equation
+:eq:`eq_cmfd_sys` can be represented in operator form as 
+
+.. math::
+   :label: eq_CMFDopers
+
+   \mathbb{M}\mathbf{\Phi} = \frac{1}{k}\mathbb{F}\mathbf{\Phi},
+
+where :math:`\mathbb{M}` is the neutron loss matrix operator,
+:math:`\mathbb{F}` is the neutron production matrix operator,
+:math:`\mathbf{\Phi}` is the multigroup flux vector and :math:`k` is the
+eigenvalue. This generalized eigenvalue problem is solved to obtain fundamental
+mode multigroup fluxes and eigenvalue. In order to produce consistent results
+with transport theory from these equations, the neutron balance equation must
+have been satisfied by MC tallies. The desire is that CMFD equations will
+produce a more accurate source than MC after each fission source generation.
+
+CMFD Feedback
+-------------
+
+Now that a more accurate representation of the expected source distribution is
+estimated from CMFD, it needs to be communicated back to MC. The first step
+in this process is to generate a probability mass function that provides
+information about how probable it is for a neutron to be born in a given cell
+and energy group. This is represented as
+
+.. math::
+   :label: eq_cmfd_psrc
+
+   p_{l,m,n}^g =
+   \frac{\sum_{h=1}^{G}\overline{\overline{\nu_f\Sigma}}^{h\rightarrow
+   g}_{f_{l,m,n}}\overline{\overline{\phi}}_{l,m,n}^h\Delta_l^u\Delta_m^v
+   \Delta_n^w}{\sum_n\sum_m\sum_l\sum_{h=1}^{G}\overline{ 
+   \overline{\nu_f\Sigma}}^{h\rightarrow
+   g}_{f_{l,m,n}}\overline{\overline{\phi}}_{l,m,n}^h\Delta_l^u\Delta_m^v
+   \Delta_n^w}.
+
+This equation can be multiplied by the number of source neutrons to obtain an
+estimate of the expected number of neutrons to be born in a given cell and
+energy group. This distribution can be compared to the MC source distribution
+to generate weight adjusted factors defined as
+
+.. math::
+   :label: eq_waf
+
+   f_{l,m,n}^g = \frac{Np_{l,m,n}^g}{\sum\limits_s w_s};\quad s\in
+   \left(g,l,m,n\right).
+
+The MC source distribution is represented on the same coarse mesh as
+CMFD by summing all neutrons' weights, :math:`w_s`, in a given cell and
+energy group. MC source weights can then be modified by this weight
+adjustment factor so that it matches the CMFD solution on the coarse
+mesh,
+
+.. math::
+   :label: src_mod
+
+   w^\prime_s = w_s\times f_{l,m,n}^g;\quad s\in \left(g,l,m,n\right).
+
+It should be noted that heterogeneous information about local coordinates and
+energy remain constant throughout this modification process.
+
+------------------------
+Implementation in OpenMC
+------------------------
+
+The section describes how CMFD was implemented in OpenMC. Before the simulation
+begins, a user sets up a CMFD input file that contains the following basic
+information:
+
+* CMFD mesh (space and energy),
+* boundary conditions at edge of mesh (albedos),
+* acceleration region (subset of mesh, optional),
+* fission source generation (FSG)/batch that CMFD should begin, and
+* whether CMFD feedback should be applied.
+
+It should be noted that for more difficult simulations (e.g., light water
+reactors), there are other options available to users such as tally resetting
+parameters, effective down-scatter usage, tally estimator, etc. For more
+information please see [Herman_Thesis]_.
+
+Of the options described above, the optional acceleration subset region is an
+uncommon feature. Because OpenMC only has a structured Cartesian mesh, mesh
+cells may overlay regions that don't contain fissionable material and may be so
+far from the core that the neutron flux is very low. If these regions were
+included in the CMFD solution, bad estimates of diffusion parameters may result
+and affect CMFD feedback. To deal with this, a user can carve out an active
+acceleration region from their structured Cartesian mesh. This is illustrated
+in diagram below. When placing a CMFD mesh over a geometry, the boundary
+conditions must be known at the global edges of the mesh. If the geometry is
+complex like the one below, one may have to cover the whole geometry including
+the reactor pressure vessel because we know that there is a zero incoming
+current boundary condition at the outer edge of the pressure vessel. This is
+not viable in practice because neutrons in simulations may not reach mesh cells
+that are near the pressure vessel. To circumvent this, one can shrink the mesh
+to cover just the core region as shown in the diagram. However, one must still
+estimate the boundary conditions at the global boundaries, but at these
+locations, they are not readily known. In OpenMC, one can carve out the active
+core region from the entire structured Cartesian mesh. This is shown in the
+diagram below by the darkened region over the core. The albedo boundary
+conditions at the active core/reflector boundary can be tallied indirectly
+during the MC simulation with incoming and outgoing partial currents. This
+allows the user to not have to worry about neutrons producing adequate tallies
+in mesh cells far away from the core.
+
+.. tikz:: Diagram of CMFD acceleration mesh
+   :libs: shapes, snakes, shadows, arrows, calc, decorations.markings, patterns, fit, matrix, spy
+   :include: cmfd_tikz/meshfig.tikz
+
+During an MC simulation, CMFD tallies are accumulated. The basic tallies
+needed are listed in Table :ref:`tab_tally`. Each tally is performed on a
+spatial and energy mesh basis. The surface area-integrated net current is
+tallied on every surface of the mesh. OpenMC tally objects are created by
+the CMFD code internally, and cross sections are calculated at each
+CMFD feedback iteration. The first CMFD iteration, controlled by the
+user, occurs just after tallies are communicated to the master processor. Once
+tallies are collapsed, cross sections, diffusion coefficients and equivalence
+parameters are calculated. This is performed only on the acceleration region if
+that option has been activated by the user. Once all diffusion parameters
+are calculated, CMFD matrices are formed where energy groups are the inner
+most iteration index. In OpenMC, compressed row storage sparse matrices are
+used due to the sparsity of CMFD operators. An example of this sparsity is
+shown for the 3-D BEAVRS model in figures :ref:`fig_loss` and :ref:`fig_prod`. These matrices
+represent an assembly radial mesh, 24 cell mesh in the axial direction and two
+energy groups. The loss matrix is 99.92% sparse and the production matrix is
+99.99% sparse. Although the loss matrix looks like it is tridiagonal, it is
+really a seven banded matrix with a block diagonal matrix for scattering. The
+production matrix is a :math:`2\times 2` block diagonal; however, zeros are present
+because no fission neutrons appear with energies in the thermal group.
+
+.. _tab_tally:
+
+.. table:: OpenMC CMFD tally list
+
+   +--------------------------------------------------------------------------------------------+----------------+---------------------------+
+   +--------------------------------------------------------------------------------------------+----------------+---------------------------+
+   | tally                                                                                      |  score         |  filter                   |
+   +--------------------------------------------------------------------------------------------+----------------+---------------------------+
+   | \ :math:`\left\langle\overline{\overline\phi}_{l,m,n}^g                                    | flux           | mesh, energy              |
+   | \Delta_l^u\Delta_m^v\Delta_n^w\right\rangle`                                               |                |                           |
+   +--------------------------------------------------------------------------------------------+----------------+---------------------------+
+   | \ :math:`\left\langle\overline{\overline\Sigma}_{t_{l,m,n}}^g                              | total          | mesh, energy              |
+   | \overline{\overline\phi}_{l,m,n}^g\Delta_l^u\Delta_m^v\Delta_n^w\right\rangle`             |                |                           |
+   +--------------------------------------------------------------------------------------------+----------------+---------------------------+
+   | \ :math:`\left\langle\overline{\overline{\nu_s\Sigma}}_{s1_{l,m,n}}^g                      | nu-scatter-1   | mesh, energy              |
+   | \overline{\overline\phi}_{l,m,n}^g\Delta_l^u\Delta_m^v\Delta_n^w\right\rangle`             |                |                           |
+   +--------------------------------------------------------------------------------------------+----------------+---------------------------+
+   | \ :math:`\left\langle\overline{\overline{\nu_s\Sigma}}_{s_{l,m,n}}^{h\rightarrow g}        | nu-scatter     | mesh, energy, energyout   |
+   | \overline{\overline\phi}_{l,m,n}^h\Delta_l^u\Delta_m^v\Delta_n^w\right\rangle`             |                |                           |
+   +--------------------------------------------------------------------------------------------+----------------+---------------------------+
+   | \ :math:`\left\langle\overline{\overline{\nu_f\Sigma}}_{f_{l,m,n}}^{h\rightarrow g}        | nu-fission     | mesh, energy, energyout   |
+   | \overline{\overline\phi}_{l,m,n}^h\Delta_l^u\Delta_m^v\Delta_n^w\right\rangle`             |                |                           |
+   +--------------------------------------------------------------------------------------------+----------------+---------------------------+
+   | \ :math:`\left\langle\overline{J}^{u,g}_{l\pm 1/2,m,n}\Delta_m^v\Delta_n^w\right\rangle`   | current        | mesh, energy              |
+   +--------------------------------------------------------------------------------------------+----------------+---------------------------+
+
+.. _fig_loss:
+
+.. figure:: ../_images/loss.png
+   :scale: 50
+
+   Sparsity of Neutron Loss Operator 
+
+.. _fig_prod:
+
+.. figure:: ../_images/prod.png
+   :scale: 50
+
+   Sparsity of Neutron Production Operator 
+
+To solve the eigenvalue problem with these matrices, different source iteration
+and linear solvers can be used. The most common source iteration solver used is
+standard power iteration as described in [Gill]_. To accelerate these source
+iterations, a Wielandt shift scheme can be used as discussed in [Park]_.  PETSc
+solvers were first implemented to perform the linear solution in parallel that
+occurs once per source iteration. When using PETSc, different types of parallel
+linear solvers and preconditioners can be used. By default, OpenMC uses an
+incomplete LU preconditioner and a GMRES Krylov solver. After some initial
+studies of parallelization with PETSc, it was observed that because CMFD
+matrices are very sparse, solution times do not scale well. An additional
+Gauss-Seidel linear solver with Chebyshev acceleration was added that is
+similar to the one used for CMFD in CASMO [Rhodes]_ and [Smith]_. This solver
+was implemented with a custom section for two energy groups. Because energy
+group is the inner most index, a block diagonal is formed when using more than
+one group. For two groups, it is easy to invert this diagonal analytically
+inside the Gauss-Seidel iterative solver. For more than two groups, this
+analytic inversion can still be performed, but with more computational effort.
+A standard Gauss-Seidel solver is used for more than two groups.
+
+Besides a power iteration, a Jacobian-free Newton-Krylov method was also
+implemented to obtain eigenvalue and multigroup fluxes as described in [Gill]_
+and [Knoll]_. This method is not the primary one used, but has gotten recent
+attention due to its coupling advantages to other physics such as thermal
+hydraulics. Once multigroup fluxes are obtained, a normalized fission source is
+calculated in the code using eq. :eq:`eq_cmfd_psrc` directly. 
+
+The next step in the process is to compute weight adjustment factors. These are
+calculated by taking the ratio of the expected number of neutrons from the CMFD
+source distribution to the current number of neutrons in each mesh. It is
+straightforward to compute the CMFD number of neutrons because it is the
+product between the total starting initial weight of neutrons and the CMFD
+normalized fission source distribution. To compute the number of neutrons from
+the current MC source, a subroutine was implemented to sum the statistical
+weights of neutrons from the source bank on a given spatial and energy mesh.
+Once weight adjustment factors were calculated, each neutron's statistical
+weight in the source bank was modified according to its location and energy.
+
 ----------
 References
 ----------
 
+.. [Gill] Daniel F. Gill. *Newton-Krylov methods for the solution of the k-eigenvalue problem in
+          multigroup neutronics calculations*. Ph.D. thesis, Pennsylvania State University, 2010.
+
+.. [Hebert] Alain Hebert. *Applied reactor physics*. Presses Internationales Polytechnique, 
+            Montreal, 2009.
+
 .. [Herman] Bryan R. Herman, Benoit Forget, Kord Smith, and Brian N. Aviles. Improved
             diffusion coefficients generated from Monte Carlo codes. In *Proceedings of M&C
-            2013*, Sun Valley, ID, USA, May 5 - 9 2013.
+            2013*, Sun Valley, ID, USA, May 5 - 9, 2013.
+
+.. [Herman_Thesis] Bryan R. Herman. *Monte Carlo and Thermal Hydraulic Coupling using
+                   Low-Order Nonlinear Diffusion Acceleration*. Sc.D. thesis,
+                   Massachusetts Institute of Technology, 2014.
+
+.. [Knoll] D.A. Knoll, H. Park, and C. Newman. *Acceleration of k-eigenvalue/criticality
+           calculations using the Jacobian-free Newton-Krylov method*. Nuclear Science and
+           Engineering, 167:133–140, 2011.
+
+.. [Park] H. Park, D.A. Knoll, and C.K. Newman. *Nonlinear acceleration of transport
+          criticality problems*. Nuclear Science and Engineering, 172:52–65, 2012. 
+
+.. [Rhodes] Joel Rhodes and Malte Edenius. *CASMO-4 --- A Fuel Assembly Burnup Program.
+            User’s Manual*. Studsvik of America, ssp-09/443-u rev 0, proprietary edition, 2001.
+
+.. [Smith] Kord S Smith and Joel D Rhodes III. *Full-core, 2-D, LWR core calculations with
+           CASMO-4E*. In Proceedings of PHYSOR 2002, Seoul, Korea, October 7 - 10, 2002.
diff --git a/docs/source/methods/cmfd_tikz/meshfig.tikz b/docs/source/methods/cmfd_tikz/meshfig.tikz
new file mode 100644
index 000000000..1e5bbfa47
--- /dev/null
+++ b/docs/source/methods/cmfd_tikz/meshfig.tikz
@@ -0,0 +1,628 @@
+
+    % these dimensions are determined in arrow_dimms.ods
+
+    \def\scale{1.0}
+
+    \def\latWidth{0.2808363589*\scale}
+    
+    \def\RPVOR{3*\scale}
+    \def\rectW{0.75*\scale}
+    \def\RPVIR{2.8694005485*\scale}
+    \def\BarrelIR{2.4547472901*\scale}
+    \def\BarrelOR{2.5293848766*\scale}
+    \def\ShieldOR{2.6040224631*\scale}
+
+    \def\bafCIRx{0.9829272561*\scale}
+    \def\bafCIRy{2.1062726917*\scale}
+    \def\bafCORx{1.0119529842*\scale}
+    \def\bafCORy{2.1352984197*\scale}
+    \def\bafMIRx{1.8254363328*\scale}
+    \def\bafMIRy{1.5445999739*\scale}
+    \def\bafMORx{1.8544620609*\scale}
+    \def\bafMORy{1.573625702*\scale}
+    
+    \tikzset{Assembly/.style={
+        inner sep=0pt,
+        text width=\latWidth in,
+        minimum size=\latWidth in,
+        draw=black,
+        align=center
+        }
+    }
+    
+    \def\tkzRPV{(0,0) circle (\RPVIR) (0,0) circle (\RPVOR)}
+    \def\tkzBarrel{(0,0) circle (\BarrelIR) (0,0) circle (\BarrelOR)}
+    \def\tkzShields{(0,0) circle (\BarrelOR) (0,0) circle (\ShieldOR)}
+    
+    \def\tkzBaffCOR{(-\bafCORx, -\bafCORy) rectangle (\bafCORx, \bafCORy)}
+    \def\tkzBaffCIR{(-\bafCIRx, -\bafCIRy) rectangle (\bafCIRx, \bafCIRy)} 
+    \def\tkzBaffMOR{(-\bafMORx, -\bafMORy) rectangle (\bafMORx, \bafMORy)}
+    \def\tkzBaffMIR{(-\bafMIRx, -\bafMIRy) rectangle (\bafMIRx, \bafMIRy) }
+    \def\tkzBaffleC{ \tkzBaffCIR \tkzBaffCOR }
+    \def\tkzBaffleM{ \tkzBaffMIR \tkzBaffMOR }
+
+    \def\tkzBaffCClip{\tkzBaffCIR (-\RPVOR, -\RPVOR) rectangle (\RPVOR, \RPVOR)}
+    \def\tkzBaffMClip{\tkzBaffMIR (-\RPVOR, -\RPVOR) rectangle (\RPVOR, \RPVOR)}
+
+    \def\highenr{blue!50}
+    \def\midenr{yellow!50}
+    \def\lowenr{red!50}
+    \def\lightgray{black!25}
+    \def\darkgray{black!80}
+
+      \begin{tikzpicture}[x=1in,y=1in, xshift=3in]
+\scalebox{0.6}{
+        % draw RPV, barrel, and shield panels
+        
+        \path[fill=black,even odd rule] \tkzRPV;
+        \path[fill=black,even odd rule] \tkzBarrel;
+        \begin{scope}
+          \clip[rotate around={45:(0,0)}] (-\RPVOR, -\rectW) rectangle (\RPVOR, \rectW) (-\rectW, \RPVOR) rectangle (\rectW, -\RPVOR);
+          \path[fill=black,even odd rule] \tkzShields;
+        \end{scope}
+
+        
+        % draw assembly row/column headers
+        
+        \draw[red, thick] ($(-7*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {R} -- ($(-7*\latWidth,4*\latWidth)$);
+        \draw[red, thick] ($(-6*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {P} -- ($(-6*\latWidth,6*\latWidth)$);
+        \draw[red, thick] ($(-5*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {N} -- ($(-5*\latWidth,7*\latWidth)$);
+        \draw[red, thick] ($(-4*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {M} -- ($(-4*\latWidth,7*\latWidth)$);
+        \draw[red, thick] ($(-3*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {L} -- ($(-3*\latWidth,8*\latWidth)$);
+        \draw[red, thick] ($(-2*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {K} -- ($(-2*\latWidth,8*\latWidth)$);
+        \draw[red, thick] ($(-1*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {J} -- ($(-1*\latWidth,8*\latWidth)$);
+        \draw[red, thick] ($(-0*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {H} -- ($(-0*\latWidth,8*\latWidth)$);
+        \draw[red, thick] ($(1*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {G} -- ($(1*\latWidth,8*\latWidth)$);
+        \draw[red, thick] ($(2*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {F} -- ($(2*\latWidth,8*\latWidth)$);
+        \draw[red, thick] ($(3*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {E} -- ($(3*\latWidth,8*\latWidth)$);
+        \draw[red, thick] ($(4*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {D} -- ($(4*\latWidth,7*\latWidth)$);
+        \draw[red, thick] ($(5*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {C} -- ($(5*\latWidth,7*\latWidth)$);
+        \draw[red, thick] ($(6*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {B} -- ($(6*\latWidth,6*\latWidth)$);
+        \draw[red, thick] ($(7*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[above, anchor=south] {A} -- ($(7*\latWidth,4*\latWidth)$);
+        
+        \begin{scope}[rotate=90]
+          \draw[red, thick] ($(-7*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {15} -- ($(-7*\latWidth,4*\latWidth)$);
+          \draw[red, thick] ($(-6*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {14} -- ($(-6*\latWidth,6*\latWidth)$);
+          \draw[red, thick] ($(-5*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {13} -- ($(-5*\latWidth,7*\latWidth)$);
+          \draw[red, thick] ($(-4*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {12} -- ($(-4*\latWidth,7*\latWidth)$);
+          \draw[red, thick] ($(-3*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {11} -- ($(-3*\latWidth,8*\latWidth)$);
+          \draw[red, thick] ($(-2*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {10} -- ($(-2*\latWidth,8*\latWidth)$);
+          \draw[red, thick] ($(-1*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {9} -- ($(-1*\latWidth,8*\latWidth)$);
+          \draw[red, thick] ($(-0*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {8} -- ($(-0*\latWidth,8*\latWidth)$);
+          \draw[red, thick] ($(1*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {7} -- ($(1*\latWidth,8*\latWidth)$);
+          \draw[red, thick] ($(2*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {6} -- ($(2*\latWidth,8*\latWidth)$);
+          \draw[red, thick] ($(3*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {5} -- ($(3*\latWidth,8*\latWidth)$);
+          \draw[red, thick] ($(4*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {4} -- ($(4*\latWidth,7*\latWidth)$);
+          \draw[red, thick] ($(5*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {3} -- ($(5*\latWidth,7*\latWidth)$);
+          \draw[red, thick] ($(6*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {2} -- ($(6*\latWidth,6*\latWidth)$);
+          \draw[red, thick] ($(7*\latWidth,\RPVOR/\latWidth*\latWidth)$) node[left, anchor=east] {1} -- ($(7*\latWidth,4*\latWidth)$);
+        \end{scope}
+        
+        % draw fuel assembly nodes
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-7*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-6*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-5*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-4*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-3*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-2*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-1*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-0*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 1*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 2*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 3*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 4*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 5*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 6*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 7*\latWidth,8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,8*\latWidth)$) {};
+        
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-7*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-6*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-5*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-4*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-3*\latWidth,7*\latWidth)$) {}; % L1
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-2*\latWidth,7*\latWidth)$) {6}; % K1
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-1*\latWidth,7*\latWidth)$) {}; % J1
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-0*\latWidth,7*\latWidth)$) {6}; % H1
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 1*\latWidth,7*\latWidth)$) {}; % G1
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 2*\latWidth,7*\latWidth)$) {6}; % F1
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 3*\latWidth,7*\latWidth)$) {}; % E1
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 4*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 5*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 6*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 7*\latWidth,7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,7*\latWidth)$) {};
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-7*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-6*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-5*\latWidth,6*\latWidth)$) {}; % N2
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-5*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-4*\latWidth,6*\latWidth)$) {}; % M2
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-4*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-3*\latWidth,6*\latWidth)$) {16}; % L2
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-2*\latWidth,6*\latWidth)$) {}; % K2
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-1*\latWidth,6*\latWidth)$) {20}; % J2
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-0*\latWidth,6*\latWidth)$) {}; % H2
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 1*\latWidth,6*\latWidth)$) {20}; % G2
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 2*\latWidth,6*\latWidth)$) {}; % F2
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 3*\latWidth,6*\latWidth)$) {16}; % E2
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 4*\latWidth,6*\latWidth)$) {}; % D2
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 4*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 5*\latWidth,6*\latWidth)$) {}; % C2
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 5*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 6*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 7*\latWidth,6*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,6*\latWidth)$) {};
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-7*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-6*\latWidth,5*\latWidth)$) {}; % P3
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-6*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-5*\latWidth,5*\latWidth)$) {15}; % N3
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-5*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-4*\latWidth,5*\latWidth)$) {16}; % M3
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-4*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-3*\latWidth,5*\latWidth)$) {}; % L3
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-2*\latWidth,5*\latWidth)$) {16}; % K3
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-1*\latWidth,5*\latWidth)$) {}; % J3
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-0*\latWidth,5*\latWidth)$) {16}; % H3
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 1*\latWidth,5*\latWidth)$) {}; % G3
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 2*\latWidth,5*\latWidth)$) {16}; % F3
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 3*\latWidth,5*\latWidth)$) {}; % E3
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 4*\latWidth,5*\latWidth)$) {16}; % D3
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 4*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 5*\latWidth,5*\latWidth)$) {15}; % C3
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 5*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 6*\latWidth,5*\latWidth)$) {}; % B3
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 6*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 7*\latWidth,5*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,5*\latWidth)$) {};
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-7*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-6*\latWidth,4*\latWidth)$) {}; % P4
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-6*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-5*\latWidth,4*\latWidth)$) {16}; % N4
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-5*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-4*\latWidth,4*\latWidth)$) {}; % M4
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-4*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-3*\latWidth,4*\latWidth)$) {16}; % L4
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-2*\latWidth,4*\latWidth)$) {}; % K4
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-1*\latWidth,4*\latWidth)$) {12}; % J4
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-0*\latWidth,4*\latWidth)$) {}; % H4
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 1*\latWidth,4*\latWidth)$) {12}; % G4
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 2*\latWidth,4*\latWidth)$) {}; % F4
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 3*\latWidth,4*\latWidth)$) {16}; % E4
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 4*\latWidth,4*\latWidth)$) {}; % D4
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 4*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 5*\latWidth,4*\latWidth)$) {16}; % C4
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 5*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 6*\latWidth,4*\latWidth)$) {}; % B4
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 6*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 7*\latWidth,4*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,4*\latWidth)$) {};
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-7*\latWidth,3*\latWidth)$) {}; % R5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-7*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-6*\latWidth,3*\latWidth)$) {16}; % P5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-6*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-5*\latWidth,3*\latWidth)$) {}; % N5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-5*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-4*\latWidth,3*\latWidth)$) {16}; % M5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-4*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-3*\latWidth,3*\latWidth)$) {}; % L5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-2*\latWidth,3*\latWidth)$) {12}; % K5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-1*\latWidth,3*\latWidth)$) {}; % J5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-0*\latWidth,3*\latWidth)$) {12}; % H5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 1*\latWidth,3*\latWidth)$) {}; % G5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 2*\latWidth,3*\latWidth)$) {12}; % F5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 3*\latWidth,3*\latWidth)$) {}; % E5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 4*\latWidth,3*\latWidth)$) {16}; % D5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 4*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 5*\latWidth,3*\latWidth)$) {}; % C5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 5*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 6*\latWidth,3*\latWidth)$) {16}; % B5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 6*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 7*\latWidth,3*\latWidth)$) {}; % A5
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 7*\latWidth,3*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,3*\latWidth)$) {};
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-7*\latWidth,2*\latWidth)$) {6}; % R6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-7*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-6*\latWidth,2*\latWidth)$) {}; % P6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-6*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-5*\latWidth,2*\latWidth)$) {16}; % N6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-5*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-4*\latWidth,2*\latWidth)$) {}; % M6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-4*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-3*\latWidth,2*\latWidth)$) {12}; % L6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-2*\latWidth,2*\latWidth)$) {}; % K6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-1*\latWidth,2*\latWidth)$) {12}; % J6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-0*\latWidth,2*\latWidth)$) {}; % H6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 1*\latWidth,2*\latWidth)$) {12}; % G6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 2*\latWidth,2*\latWidth)$) {}; % F6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 3*\latWidth,2*\latWidth)$) {12}; % E6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 4*\latWidth,2*\latWidth)$) {}; % D6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 4*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 5*\latWidth,2*\latWidth)$) {16}; % C6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 5*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 6*\latWidth,2*\latWidth)$) {}; % B6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 6*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 7*\latWidth,2*\latWidth)$) {6}; % A6
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 7*\latWidth,2*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,2*\latWidth)$) {};
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-7*\latWidth,1*\latWidth)$) {}; % R7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-7*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-6*\latWidth,1*\latWidth)$) {20}; % P7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-6*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-5*\latWidth,1*\latWidth)$) {}; % N7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-5*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-4*\latWidth,1*\latWidth)$) {12}; % M7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-4*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-3*\latWidth,1*\latWidth)$) {}; % L7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-2*\latWidth,1*\latWidth)$) {12}; % K7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-1*\latWidth,1*\latWidth)$) {}; % J7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-0*\latWidth,1*\latWidth)$) {16}; % H7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 1*\latWidth,1*\latWidth)$) {}; % G7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 2*\latWidth,1*\latWidth)$) {12}; % F7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 3*\latWidth,1*\latWidth)$) {}; % E7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 4*\latWidth,1*\latWidth)$) {12}; % D7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 4*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 5*\latWidth,1*\latWidth)$) {}; % C7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 5*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 6*\latWidth,1*\latWidth)$) {20}; % B7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 6*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 7*\latWidth,1*\latWidth)$) {}; % A7
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 7*\latWidth,1*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,1*\latWidth)$) {};
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-7*\latWidth,0*\latWidth)$) {6}; % R8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-7*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-6*\latWidth,0*\latWidth)$) {}; % P8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-6*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-5*\latWidth,0*\latWidth)$) {16}; % N8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-5*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-4*\latWidth,0*\latWidth)$) {}; % M8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-4*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-3*\latWidth,0*\latWidth)$) {12}; % L8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-2*\latWidth,0*\latWidth)$) {}; % K8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-1*\latWidth,0*\latWidth)$) {16}; % J8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-0*\latWidth,0*\latWidth)$) {}; % H8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 1*\latWidth,0*\latWidth)$) {16}; % G8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 2*\latWidth,0*\latWidth)$) {}; % F8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 3*\latWidth,0*\latWidth)$) {12}; % E8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 4*\latWidth,0*\latWidth)$) {}; % D8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 4*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 5*\latWidth,0*\latWidth)$) {16}; % C8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 5*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 6*\latWidth,0*\latWidth)$) {}; % B8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 6*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 7*\latWidth,0*\latWidth)$) {6}; % A8
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 7*\latWidth,0*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,0*\latWidth)$) {};
+
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-7*\latWidth,-1*\latWidth)$) {}; % R9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-7*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-6*\latWidth,-1*\latWidth)$) {20}; % P9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-6*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-5*\latWidth,-1*\latWidth)$) {}; % N9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-5*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-4*\latWidth,-1*\latWidth)$) {12}; % M9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-4*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-3*\latWidth,-1*\latWidth)$) {}; % L9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-2*\latWidth,-1*\latWidth)$) {12}; % K9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-1*\latWidth,-1*\latWidth)$) {}; % J9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-0*\latWidth,-1*\latWidth)$) {16}; % H9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 1*\latWidth,-1*\latWidth)$) {}; % G9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 2*\latWidth,-1*\latWidth)$) {12}; % F9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 3*\latWidth,-1*\latWidth)$) {}; % E9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 4*\latWidth,-1*\latWidth)$) {12}; % D9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 4*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 5*\latWidth,-1*\latWidth)$) {}; % C9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 5*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 6*\latWidth,-1*\latWidth)$) {20}; % B9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 6*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 7*\latWidth,-1*\latWidth)$) {}; % A9
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 7*\latWidth,-1*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,-1*\latWidth)$) {};
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-7*\latWidth,-2*\latWidth)$) {6}; % R10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-7*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-6*\latWidth,-2*\latWidth)$) {}; % P10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-6*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-5*\latWidth,-2*\latWidth)$) {16}; % N10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-5*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-4*\latWidth,-2*\latWidth)$) {}; % M10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-4*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-3*\latWidth,-2*\latWidth)$) {12}; % L10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-2*\latWidth,-2*\latWidth)$) {}; % K10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-1*\latWidth,-2*\latWidth)$) {12}; % J10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-0*\latWidth,-2*\latWidth)$) {}; % H10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 1*\latWidth,-2*\latWidth)$) {12}; % G10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 2*\latWidth,-2*\latWidth)$) {}; % F10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 3*\latWidth,-2*\latWidth)$) {12}; % E10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 4*\latWidth,-2*\latWidth)$) {}; % D10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 4*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 5*\latWidth,-2*\latWidth)$) {16}; % C10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 5*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 6*\latWidth,-2*\latWidth)$) {}; % B10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 6*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 7*\latWidth,-2*\latWidth)$) {6}; % A10
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 7*\latWidth,-2*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,-2*\latWidth)$) {};
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-7*\latWidth,-3*\latWidth)$) {}; % R11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-7*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-6*\latWidth,-3*\latWidth)$) {16}; % P11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-6*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-5*\latWidth,-3*\latWidth)$) {}; % N11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-5*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-4*\latWidth,-3*\latWidth)$) {16}; % M11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-4*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-3*\latWidth,-3*\latWidth)$) {}; % L11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-2*\latWidth,-3*\latWidth)$) {12}; % K11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-1*\latWidth,-3*\latWidth)$) {}; % J11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-0*\latWidth,-3*\latWidth)$) {12}; % H11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 1*\latWidth,-3*\latWidth)$) {}; % G11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 2*\latWidth,-3*\latWidth)$) {12}; % F11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 3*\latWidth,-3*\latWidth)$) {}; % E11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 4*\latWidth,-3*\latWidth)$) {16}; % D11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 4*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 5*\latWidth,-3*\latWidth)$) {}; % C11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 5*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 6*\latWidth,-3*\latWidth)$) {16}; % B11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 6*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 7*\latWidth,-3*\latWidth)$) {}; % A11
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 7*\latWidth,-3*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,-3*\latWidth)$) {};
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-7*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-6*\latWidth,-4*\latWidth)$) {}; % P12
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-6*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-5*\latWidth,-4*\latWidth)$) {16}; % N12
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-5*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-4*\latWidth,-4*\latWidth)$) {}; % M12
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-4*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-3*\latWidth,-4*\latWidth)$) {16}; % L12
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-2*\latWidth,-4*\latWidth)$) {}; % K12
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-1*\latWidth,-4*\latWidth)$) {12}; % J12
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-0*\latWidth,-4*\latWidth)$) {}; % H12
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 1*\latWidth,-4*\latWidth)$) {12}; % G12
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 2*\latWidth,-4*\latWidth)$) {}; % F12
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 3*\latWidth,-4*\latWidth)$) {16}; % E12
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 4*\latWidth,-4*\latWidth)$) {}; % D12
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 4*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 5*\latWidth,-4*\latWidth)$) {16}; % C12
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 5*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 6*\latWidth,-4*\latWidth)$) {}; % B12
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 6*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 7*\latWidth,-4*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,-4*\latWidth)$) {};
+
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-7*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-6*\latWidth,-5*\latWidth)$) {}; % P13
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-6*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-5*\latWidth,-5*\latWidth)$) {15}; % N13
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-5*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-4*\latWidth,-5*\latWidth)$) {16}; % M13
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-4*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-3*\latWidth,-5*\latWidth)$) {}; % L13
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-2*\latWidth,-5*\latWidth)$) {16}; % K13
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-1*\latWidth,-5*\latWidth)$) {}; % J13
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($(-0*\latWidth,-5*\latWidth)$) {16}; % H13
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 1*\latWidth,-5*\latWidth)$) {}; % G13
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 2*\latWidth,-5*\latWidth)$) {16}; % F13
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 3*\latWidth,-5*\latWidth)$) {}; % E13
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\midenr] at ($( 4*\latWidth,-5*\latWidth)$) {16}; % D13
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 4*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 5*\latWidth,-5*\latWidth)$) {15}; % C13
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 5*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 6*\latWidth,-5*\latWidth)$) {}; % B13
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 6*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 7*\latWidth,-5*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,-5*\latWidth)$) {};
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-7*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-6*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-5*\latWidth,-6*\latWidth)$) {}; % N14
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-5*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-4*\latWidth,-6*\latWidth)$) {}; % M14
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-4*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-3*\latWidth,-6*\latWidth)$) {16}; % L14
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-2*\latWidth,-6*\latWidth)$) {}; % K14
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-1*\latWidth,-6*\latWidth)$) {20}; % J14
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($(-0*\latWidth,-6*\latWidth)$) {}; % H14
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 1*\latWidth,-6*\latWidth)$) {20}; % G14
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\lowenr] at ($( 2*\latWidth,-6*\latWidth)$) {}; % F14
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 3*\latWidth,-6*\latWidth)$) {16}; % E14
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 4*\latWidth,-6*\latWidth)$) {}; % D14
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 4*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 5*\latWidth,-6*\latWidth)$) {}; % C14
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 5*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 6*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 7*\latWidth,-6*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,-6*\latWidth)$) {};
+
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-7*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-6*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-5*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-4*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-3*\latWidth,-7*\latWidth)$) {}; % L15
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-3*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-2*\latWidth,-7*\latWidth)$) {6}; % K15
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-2*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-1*\latWidth,-7*\latWidth)$) {}; % J15
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-1*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($(-0*\latWidth,-7*\latWidth)$) {6}; % H15
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($(-0*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 1*\latWidth,-7*\latWidth)$) {}; % G15
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 1*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 2*\latWidth,-7*\latWidth)$) {6}; % F15
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 2*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\highenr] at ($( 3*\latWidth,-7*\latWidth)$) {}; % E15
+        \node [Assembly, fill=\darkgray, opacity=0.7] at ($( 3*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 4*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 5*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 6*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 7*\latWidth,-7*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,-7*\latWidth)$) {};
+        
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-8*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-7*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-6*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-5*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-4*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-3*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-2*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-1*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($(-0*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 1*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 2*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 3*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 4*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 5*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 6*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 7*\latWidth,-8*\latWidth)$) {};
+        \node [Assembly, fill=\lightgray, opacity=0.3] at ($( 8*\latWidth,-8*\latWidth)$) {};
+
+        % draw baffle north/south
+        
+        \begin{scope}[even odd rule]
+          \clip[rotate=90] \tkzBaffMClip;
+          \path[fill=black] \tkzBaffleC;
+        \end{scope}
+        \begin{scope}[even odd rule]
+          \clip \tkzBaffCClip;
+          \clip \tkzBaffMClip;
+          \path[fill=black, rotate=90] \tkzBaffleM;
+        \end{scope}
+        
+        % draw baffle east/west
+        
+        \begin{scope}[rotate=90]
+          \begin{scope}[even odd rule]
+            \clip[rotate=90] \tkzBaffMClip;
+            \path[fill=black] \tkzBaffleC;
+          \end{scope}
+          \begin{scope}[even odd rule]
+            \clip \tkzBaffCClip;
+            \clip \tkzBaffMClip;
+            \path[fill=black, rotate=90] \tkzBaffleM;
+          \end{scope}
+        \end{scope}}
+      \end{tikzpicture}