From 1660bae22a252eadbf535b1dc7c8667a161f9f43 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Sun, 8 May 2022 21:15:32 -0400
Subject: [PATCH] add examples import one image to manual

---
 doc/src/create_atoms.rst                   |   5 +
 doc/src/img/marble_race.jpg                | Bin 0 -> 40294 bytes
 examples/README                            |   1 +
 examples/mesh/in.marble_race               |  43 ++++++
 examples/mesh/in.mesh_box                  |  48 ++++++
 examples/mesh/log.8May22.marble_race.g++.1 | 156 ++++++++++++++++++++
 examples/mesh/log.8May22.marble_race.g++.4 | 163 +++++++++++++++++++++
 examples/mesh/log.8May22.mesh_box.g++.1    | 146 ++++++++++++++++++
 examples/mesh/log.8May22.mesh_box.g++.4    | 146 ++++++++++++++++++
 examples/mesh/open_box.stl                 |  72 +++++++++
 examples/mesh/race_track.stl               |   1 +
 11 files changed, 781 insertions(+)
 create mode 100644 doc/src/img/marble_race.jpg
 create mode 100644 examples/mesh/in.marble_race
 create mode 100644 examples/mesh/in.mesh_box
 create mode 100644 examples/mesh/log.8May22.marble_race.g++.1
 create mode 100644 examples/mesh/log.8May22.marble_race.g++.4
 create mode 100644 examples/mesh/log.8May22.mesh_box.g++.1
 create mode 100644 examples/mesh/log.8May22.mesh_box.g++.4
 create mode 100644 examples/mesh/open_box.stl
 create mode 120000 examples/mesh/race_track.stl

diff --git a/doc/src/create_atoms.rst b/doc/src/create_atoms.rst
index 3edbe46444..c32a1510ee 100644
--- a/doc/src/create_atoms.rst
+++ b/doc/src/create_atoms.rst
@@ -126,6 +126,11 @@ the specified coordinates.  This can be useful for debugging purposes
 or to create a tiny system with a handful of particles at specified
 positions.
 
+.. figure:: img/marble_race.jpg
+            :figwidth: 33%
+            :align: right
+            :target: _images/marble_race.jpg
+
 For the *mesh* style, a file with a triangle mesh in `ASCII STL format
 <https://en.wikipedia.org/wiki/STL_(file_format)#ASCII_STL>`_ is read
 and one or more particles are placed into the area of each triangle.
diff --git a/doc/src/img/marble_race.jpg b/doc/src/img/marble_race.jpg
new file mode 100644
index 0000000000000000000000000000000000000000..f3d238e0f25140ae295f7f983c4f0f29cb3d0c5a
GIT binary patch
literal 40294
zcmeFYXH-*NxA+@+5$R|G0fGtAm2N^4P`Y#oy-4qn&_R$QC>ToUVCcQ~E=uoBs#NJn
zRjE?MKkqqr+<V6P@;>kVbU)mcG1tf*YmJ>T*WQ_H&flEt@7KR80BU7<C3yfI9sq!M
z`vLsD0a!iLl9SO<*Aj#1BOzYaXm@*87l_a!0sfbK4xY{yuJZOiJaR5}5LX+BAU{7`
zKu|yw!fEf~VU4zTX5-?Kf8qSn-4O!O^nB@L@8JuPggg=9g=jmwx_a2zyV%}d^OOHw
z06YVb5)+dU6OodTkdTp)l2g#%rMPp40z^wkMbFI0%EHVD218(vI3N#rp<poQ6E0qU
z0byZbRt|(DTu|bXkg(vtZh}WfMn-Xm0(kc>P>>DGCiow>zwH1TQv5*tUj%p$0QfX`
z1T=Vmy8*1Xt3ibKUk&hoH9UL*LLy=kQZn*8w-f5A0r+?X1o(smL_~yyx3dFo#{q;i
zM6_%IGQ@P6FG(Jt=mmqZ*`!d}>JA32FF#>I7U&Q%@_YAzAV&6w9GqO-!Xlzi;0Up2
z&*kJ56qS^<b##$>`UZxUR@OGQcJ>bL9-dy_KE8gTZ^FXgzKe*Ai%&>QN=`{l%gN2l
zFDNW3E~%-lt8Zv*YHsQ5>h9_7>mL}Ln4FrPnf*F9zq+=*vAMPVeP{RA(ecUY+4=8_
z%YW&G2O#)wwQl?WR_uT2MRThcJ|Q6iA<4h=!o&Bz6#^PUA~peHS{Y4}mngaig2ANp
zve@kE4l<~a)=vft^cV7bFyU49U;k3=zcl-2iiP~Y((M0M?Eleg9za2Wcl+=NXaG+E
z7dMfdS@{3I{C|86xTo$V-&eg#F_opH?CEEXQ7-D$_KIw!vv<+WX@S3bjNW`v)B8SJ
zdO#$f^{o><#RX<)mTz{?G@Up3zP28dU&;`js%LYL?&s=7;Bd~e+2e|Xx4_GFhWYaP
zV5Yz!0MG2XWlSN%gC(MVNjwOv;bHhQhHK?aw4hH#{|29j(Owazx*I=3evMd8)p_-*
z+1Ba<H6A8u`{@f7$GrP+Ur^vmFKR)Ab^hL(C;3-%vi|bigqP4EuZ8OYpGZt^^Usa(
z1nX25iR(TW@~f-qo}VMxK{mdxh64@-G7H4~1M6S@XnYdq16zqyd0*}DRr`phud&_{
z%z{s;wb^p&N4u-*rN^yrha7;G&n*x7;JpL>p|qc!XGlAUyJt81RM1w!-!|4^aJ1m-
z08Y`D<dY;flm8DfbECt`r!fvvm7y9s&G%ZRKKr;&eHIx)?fBUTu{Au)P(qAt<ED93
zxhIy7)CFEmFC|cok`y-o1?*`Vo2C>jfL5FF^ZgU<rHqZtB7W#SE>%`p)QgFu-Bt-0
z5aD~&ySBI=0ZpJv$Rfye6PrpBO%gC7q=vF=PSHAql^<RZ2n<GtPEDAYzRBeDu^jjN
za!|Evx+G_2zM0v9HaHm)3XCJ{HtPYdDxKc`<m>oABNZ|>Ytb(N1?!IMJ`wc;PG@qu
zV;1)xy5*I3RKg{=Hf)P?gt{i3IUoD?Gk!4IS1+_rLb)7U(l)He@KdHY?My_etJEg+
ze&LL$9qj@31&Y}&=eD+zJx(bM#r$#5KcrQ$p+FCdUFxuhyLljxp&x&Ypu7F#Q`i6t
z;qc7}#M5YpH;l}=ru{0lvQ%U@91}@Q6BGLuJ`$6MkGh3?BYm6#i6eMn>mWFfLjV0b
zLt|AHnYyD;!%V&-Vo(Ih1y=QK=TBio(&}~Bl04v=>S3XfZ&H43azAuO#Wi3nHfVin
z^70e8t}pIh=?3^z$Jw1vLeiamGG~CLXNT?Pac6GZ&p*=@<QCQ(zlE^BUmMOw9uHv0
zGoI01PSuzsHp%vd2*vnq`+AY<ELEgm?VB|erH&5dI8gU()mP6W#MKn9@kX@u*{yl8
zLYz13+9tThz9ExBK!eatU?TrV#g)PJ=*|KIQP9w0(S*tEszdZs%fD<->G?FwFtI4{
zt+HQOXn?0he;#H&+D0Zp6v&w%fp0Hnd5^dU{NvPr3hNi0DzM6VX!}^W?VEe}v=rHN
zSswLRjo7<8XNul`!Q5=Ue!OGq`HI1hwfnuwPthpFAsg+SNS=hMoJdZ<KbOz;1t!bK
z6bOQ-?v=oLkDbxHX7h^cH}lH-*>Wy&tF1e*2dY>q<g4wYn|jw98olGE>zfj7pZ$jr
zBvPq=+)b(mtK;=_zKy95kIWXd;CdKlGVeUw40>ezv{*C5edjM=Hp|D5@4~T==Q}mN
zB;xk<Lw@_3@@?~4^Ql4M;zm(Li5Ypt+82|=o820ZYnh@CopGc`jmi^$0iLnJheBy@
z<Ja}RUh7>%S#X^rN_d=n+%>-H!<=u%QiKeS3>Lfd23#K$!50Tn7qjdy*!`4`L(Lwn
zYYEOJa^^&C^WJuV|6J-%cwV{}7ri_x{6V!kdHSF}TWO}45w+Ax$AvDK%g_7mYog^p
zs)H#rdHVF{U%;uJbm`q1LOC7HaZUhLXFf#<{(AiQTFoifnPL$Y;Sm*(;aunE2|s;6
zPo>^;hE*YV%~YDvs18?`EvsRPD?w@sEeJ67*_21XCeHT#5YkHycECKap5{}Vnq_gs
zO>AMH0NXgJ09kVDct?J|7AUO6cyl;TIL!@Q3|aebgUQ1xMGfd-z8Lr48KD<C;iKO+
z<$2BAuKI$1Yo}(>!th}f`UF8j(c4dz2d+$OR8$xaPoC2(sxq|VF|~q|AaqKO%oDu^
zP*|gO6M41dT(mkHztTzoS7WjxnHf*|w@yh1MvSNXZx5Erf=!?DLq>0}&IZMe6PBr+
z{ViYZSHVF8YzpdrNk(}pkY=XY%--Dkwxj9vvEe0>e04sC!zWwfLumQhVZIOXI5$Ai
zYPEhyUd13PgUYa?QG6R_J46!m+96<bCxN@@qLo@)M=XJF%vo;@T=O-tj3!?&SLi{s
zjKodCRzrhbKhb*@63HE`u@b$}bC-ZXZQZOm`IJw+@b3bDxSF-zBcoUCrShA~)*8jx
zNea+pTiAWe#Pe>!_p=-lQ}^SCSJ4<sH$egIXRu%0fT$L$@pO(yY9{w%{TCKvdz>uV
zNl@RPu>LAgu%)NE>bY7xmC$`|uR6d+DYjP7l4cvT71i`3(sbg3@u*lO+uQ-#!1=+4
zcOAg(V&PH>7BiyYAIFU0%Sy<nO?{HT(!w@h*I}UUQx@^Xs_arG*CAt$@$cmf8fyi^
z!R}-9wexSaXqH-!LzRU0znGWeBS>8Roa>kA_;B&xgrCS?Cw+O_)-3F`A=mw<H8PUP
zH_3axosx76RCF^z>(+chHfEYJEz;B9+5a|{U!f>}M5*-bj-qw9b_7imU*EGS?d@X>
z)WeFVUxLx^MUBzvMSp+Jov;tw*Onnv<z~E0HLs~FZr0eF`5dcjDB&>gLOIl@rPv0h
zy+x109)6Aic>l_H4h<B_okUU$o;02pc?<~*$sdxeeeFZj7s!qBPr)`?ds~$7wYa`Q
zXi-V~6kO2yI^hwaeK<&jrwyxRH{XjwKM!3mast+8D)PRJf%Nbag9kI8&58_PU_dJP
z%_ey(E0^Hu1;oOLu0*EG>B2Mc!e+8ILaU+3M4>AJa$F`Y4^A0a`9oE?b5f!E1shq_
z_1N&_-C|5KPaIF;FNr(MC+LzOE3*$Qxu18Ra<7fn%xe^$Mk+7OJ~I}ve=-0~zqN;3
zl;i)z7#4EZicasTmNNs{rp=E@@3ekCuR>kUBxE0O%iph`myNG+`Zg&c2Gz}2cO1aH
zfkFZTR9$*~x;zghxk1}@FWFHUJR%iKiw(!${m8G1c%_KB>1VI$hF+bT%vF@-NBefP
zk2yHS=Sgp1Xfjl-6vX2e^fS~@BQhv1^t0x6LKf8XTEc9PdLpI8rN%^SX`harH0O=w
z7e>k0R^YYc!{yCH0M%q<^wCC3n9Ue8GSxCJcZ)pvqigT32CaJ!O%*DBK=N20yR;@j
zrU)$TK`%%dPzo3-u-<nj+D|PO8AcJbf><l-RcI(vUC`os2#M?k_vRIh`wlYqqm}gY
zME+sA4=v*^(cxzAgzKskjAaGr--G~GnJP$rd-b@VE-N|yDYk3~gtJq*Ee5^5obDc=
ziAn`BcSL4(ezJVhS)!|}&m(UXSm?r1DC}$cJv@@L67rp#GbiXjlJ63aLyMaD#uD5l
zeHq4P&cmTp#h;WeEZ}`pdR_)M=ju~~jX&q6j0b96CI15QZA^2P1(Ted_oDv-wDl^f
zU$v&JOV!I;MB?`fs>$Q+I_{yO&$YfD@0_!F^FQCTT&M!zM^n`!(>7egyz#<@j%S%m
z<>3ZnKg>$exz`RJjg8AA^pFM;h*ietD3mg8qs~kI;{Y(npyE-Q&$?G10RiUn3-~^j
zkhrt_nEOOSL7m;&BQ3?03<56Dm|8uzkIUs+f6B)DkpXjHY#`3ctBixR&p`Q)vT7SP
zj5)q;*eW(dnn?aNl_B<I+TcQv_C7pNP-wnEc!8-+%943lKOVsW!T9?lKjjOmT3Pt0
zNw2o={4l0sARQddNWc8xFFJFHbC<LGEoYVziwC<rEe2&ZurtD3R$?jTY`eB4M!z?k
zE3Ta6C-FpZm<gbE;&C^+(~Rda!j`2+La3RM<e%7{ZBD$1Ar2>vcm<zt@CzL)nmt29
zdj7XdZl$o)6&C4hDp2_nl26KCx`19O;Z+5I(3@G@w%&bAP$QE&I$cG_NU!;=ZSL-U
z0Xk|7<w8n`LIvPyx%TI}0D&MZxiZ34{+zg8_>aA1O#E9gfa=vcL$PSOYlYE%D>gA=
z5v0hZCfVcOCs_R2*Z$((_lmcycRP(~tv3}2GEV7CE*@EG_sRS~x}I2179M>XHxx=Q
z@AWx|G0bbR?#>&I9f2wPVs)?=k%yPvi}cB+1;!ELbp2oQUXDj~#s7E_{<>Zf4@wYI
zT*MwM3s$`Nt^m5x`D)3JOJbkJdF5-BKmRj5POqn3T{XdL2v#<ZLs)cMet+6q?db?F
zv-vZbM0lNdda2b{WFq5a+{N9iz}qi<{VOu-_){T@ULSV<Amh`w?KWSXU!G8&Mh_1$
zl%up!7;Ya1;F8!39;9%Vh}`G<DCsbha_^U|D~!_5lcOSFPwZRgnS1}*mVII)b4lLa
z&ss7go@G6T2M)IT3f~(VUgRjPMjHTMQxj600e=1dvUZ1AxNFNu7Atn-RYbu2n?YSC
zst%z8>sjN#s}ZCWV?Ca)KCDOB%Im^H5F+g*qq>O7;!1Ai4qBN_&OS*vOpC9E%;=HM
z3N>DdX5G8?AjJ`N-ilQN-9}0=7RTu*efUzJCCUAO0SZbEiW+=3UnpxBZJ%I&v@{>i
zoTXkuK-~)(O)q7GCR?%mq_ktix4IDO0D44|PXaR!j&Yb!OPkbma36#6mO^<byJ=#2
zGinKhD5A>9SbEEq?yx;D4u2LRpvvFbqn~vKFV{u5q&|7B!jba687&kdkYWF;F#bXR
zN3Me%WlIQiriYjC`-d|FgnpqqiD_Glrr`<q{cW@qVW9?6uqK0zV0v^`Twv1=L(z%C
zc|3E{j$O%8@$O09N`P8v=8oa|)-t(wvhGEWJ&&#-_<&y1jMF=eUQ?Rr--+qu%9AHw
zx+w{sEKB-qAD4)@jC$9%M(b<_UwvBQQGUdo=X=)jIV;~^oH<}hkd)sVhA)&hViyTc
z;@|+^FrXcmewvqZal|5vOD|GMFf2l5no)^}ryi?$X^Uy3!R{x!++5Y$)c#-_Ec*j_
zT?v{2`n3|p%UQ-q$7ixbkeArLiFdriHRhHLN(|^3cJCqg&>lnGqeDdwL<j2&9i8O1
zm-*RR&U`aZ1>8_mc#(r?0`ZMd{Fc-LlR^RlhT+lnp@(90SY48lCJ~Ej3~4_!10no7
zU7P#0m;uu(Sa1(7X&$YD-#krS+XTye(r8}8nN^W3Hx+>HFJQMqUfvdcFm*jwaWp%0
zFNi_AqaRu+(m8s6G*%%!IQm<GPnmm>s{Dt9iq?RhmuamvklSsvU99szMK_HSqN6R@
z^Le45yU|_IuO*PouD?yZ9cYVPs#%`Wm^3?ldX(OpI&;G$`q)LKNe*(nr$J;_Wd?HX
zW5||3b;<y2B^|=Xey{_@<q;yM@OOiP==iCX(C_%KGlZssln8K_$mQ7V`+>yU2#&q^
zYHS=mySxA<`ZV8_r@V=YQW4mvX&T7(N!QEy?#NcN5>=cIwLq&xj6+py3gjF88|=71
z8Iuy}qCR!|MrenKjfuJg?x%op(fCY=ic#3ZH4zTr5J2(TZ<?%b&ZhO~AYYSaBgB`W
ziKT8TM8=6d+F0mSujr}NscDy{Yn?Er+uYng&0O$kRD~i{8k*(g77){>C)a2<;1in!
zyvogYZpf<<P;%Q{388{Sj)YH;6?onaBMgpntjU00w1AYEF~5^Ym#XiBRzehx3Ws2V
z>Dfl^sPnPkGsRUQ1G%<y4+PdZ;C=OG{$DMlIpc(BVlwSC0Em^CzK2Z?+67AozZnbm
zEepdAh@wxTt>S;YweTI2yh7Wbtvdjr#|C3e(gS~pJfPaoLn-hG{bmBvr?^ni2i%92
zV;2VoSBwFXTQ)??IIFf`m*?NJ4!bIzST;hjCVa_>^w7IU!9h(2$j22qEe&?eGW{#j
z&HD;%?mvig2_MQg#Abm8V|`BB5qaOW=R1b`Om>h$yF2JPjusEX?GOb&&#lZtzUfxJ
z7TEhZMQM^oW~z0ASYs8DiGv}570v@gZ$F!y?;N^%5Om2p{G@Py-I08MiuO5Qs&P1l
zz?g&{=h=+c`9&yCBZNDMMB-PEq6>$GwmqSnw;Pe34<?>87xS@u8-io)nkf9Ly)R|E
zxES&aKTOZ9FS@N%i>lV?B3Fxavn4LDKCy78#vm(hm-#S25kl2<k=E=3bwt0s08^~W
z<4L6y6ix1f$C^VcWeDIHE-UxgJ1nG<Jm8((T0;HCcnlyfey8tZ$B>G7^NNgCD7b$B
zkl9&x(!{^~YJ*I%y;qQArY;2I^}Mx~(GkJ7rss%!NBO4UDY=776OjfjYQAM>dp<*n
zboZKh@z)Q7@tr_fq345ou!CNKjJXkyvS%vK71-X>iPU9@(OU@gm%ni=3>D~JLLdm9
z2A2rEu)!t|Ev|+HN6X(e5;B&-^VTyJ+I_0Ljl-qS*&E45e*Mr|7-uWM$)p?{r}#Ag
zz8EM<ZxF{Y_6uYAi-Do-3|x&EfNd3ITnlEP4rGgAbuR}a?|W&_TYn=GI8V+JQD+;I
z3{0%_^nle`TVfl*GmHXSc}GGQ$f}afT=w1tx*m!c!7~5Lvl73Phj&tDz1&$*(WkZ6
z{01EyMIz~I_Ru6B(SV<YUXtRsC<UaP)}gF7^F0k6_~pNy@aTW!gqPptusy>uU0U7v
zlwH=_(=D1PfNd%ILd3ntRt?9Ub>AyrmXjrSj*kA!Y8WNa9I;Sye}J783EWy<<?FQ`
z3?UIo>+k={;iYbP`O0aNqs$4VrNlR2Zf?~lOQMWTI3y_)0rNyC9Xi;=Zv)Px5(g%=
zpUEN&ABA#E6j(PUQ#28JtUt_BGf{V%`KVRWgQ3U>uPJL=?X#p^Qz_dI#5-YjD6PO7
zTZSMr^10(eG216)a~H4}<%u|CQ+_jjCa*0s9G)k;k>18tep5<PiJuqN8C3`>w=9cr
zlpU}8*i_>jSLVz57@X;JRw0>?>+Enab0hA^H`<QZXEwQKY*k^MxOnQtcagR=YN^WY
z^5}3r?d8JEhOnlbkFBDn96jJa5#+7+%van@7*ZWadQD`0>r59MojwHi`}D(aSQIX<
z9b5StXLi0dC%UfIX15{5y`O;gc^<8ByOAOG_A;iuHJTP>gOv7qG&(5mP^DxheM$z2
zaMeX~+kbU40N#V&iB+#;y(yUpUwN=`i89D>5^S^}6LA=XLMV58$krjjx9F!*u*Rss
zmFCzqz8?;469bJW+HSsu(DiK^Y3|G<Qy-*?QIf);JCVK3lRfD}FVFj~5qmMu3{E&_
z>syY$5pGIvGgM_-^E1Q&&#=7!hS8m!B=-B897JKReb)3}oCOt0!Qk<r&bVHiDQ3Q}
zPAjE|c0lx>%MD~$bRyin#piN^CXVW;qySOwYQo%QoA`1OSzeoxOLR%oL67%LKiuC&
z8$u+r`c@sa|GKlx{tVS0!9HLkDH|F@(#GC%aST|kB7<bgX!k&9cz7r>NV7;i{fWZA
zjybZ4%Nvl5v1H8%Yl}SxDY=>%?m1_YdUNt%X!i6Omp*DJ)-DvuS%z!2k6Dm7NQ-vv
zWLS)0qX8mqf3iyNjjwFD8yBQSxx1kks6C;@XAI+_ADqs%X<M4W;O(bjq!fioz{hiI
z8m%>IPmPgn`P*qY0l?3IJxUMVyRMIQ5fp7Qr`Pr~v&0Jyhl?bTbr3yX^=n%)b%YTg
zE)foy1|FUbl1Kx3WW0-ApN)gfkrs<|AVe2I-{eG)n1Kk9{#iJoMR5m{@jX5+{jWH@
zS%x<*uz2o|V7I$wwQc0`{*-d}@%PUj+vMjY`|(3I&_dsqYs(HzHc0u!Y5jE<-jDxc
z0I>%HC>=pSqe05+BoGI1dm?s(sleYiq4u5@+;v+kT@}fX*9Sh<sc@yWlZ}Nl#DDgG
zv#6#HKB33UN*+2J0KdM6I4u`?RlxTqgmmkJrK26SOBnDrtg*U9YM;sTiwya3clb$N
zNI8J#+w~C6gvKRKKpZeu<~w(i{-BHt9-O#RwQzABs<>)`iqqI;&ER3rx7ZAqS5*x+
zip}KNdO4X@g2h#l<P)XFz_3&9DqW%qhD|)QBRp*PGwxu1`cFTmC4IHvP8;S#z2F}$
z@W+1qBHT*HDGoI$qHPqh)`CGSl5fY%v<P?_bqq?hA@j6je$X|?$kFSKrU%Q>YvBDm
z(G$M78S#q4#y+iNP^;Xjg3Z6t>?eCRa9wLSzI9f;xxCvwxLC2RgbDyBH6cl0-et`&
z+9%9-<^HX9at*SW5OI&P2u;mdFvSVO$=7MEn?KRquNuaa2(Bi?0Rl+!Xvq;7%B5}Q
zw<2tJn`qic(B(AQ0XH&zQd3`0j}{+KrUeo&f1Ch_((M=ylB;WJR~zJdziQ@KoQiTf
zR71lL<n@~%w0PSDB53h?qGSLqvAgj!D^Hz@Vd{NFORe(R>VhuDrh%=(DVUxUWgkx0
zZm}4yC1<a+v@Dy9tz-zxYw+<&I0V8@S=PsA8F+0FAoujVn8B8$w4eHj#zqc14#fHP
z-bsa0O0A1aeKbnUSz9wIy|LjhU<=7C>C*c}j$w_N?cc59{~kBG0TNYi0kA1CrfqWO
z$n&r+W~H+XEXyr2SW5Z}Xq!A0EWZ|YSu~Pk;-MzvC@x!w8iA9eFl+T9cYlNN#>y@F
zyOj}pW(v_LOtThz$D}mvp6Hj4O%Q$p4v|moOb@|s#3b1aVFfOfHlSWVmbE(DX7F{M
zzs1MFs3+Mia`B!uc=f@AR#&XbeXT8VO2mtClxVIDONk+J+<jF1EW>!24BxA!AYFob
z*(obS4ZZrPw}5Sq@$Ixp-(jUX%jnZk@WlM0utT7pM9gc$I$DK+IFnLL@C%}(U5+L_
zrwtVDCrO(+v?fd+N9|`@Z!`QWj7?2x*WQ^S+tWQxyfQ26?!$U1BXFEz7pvJ#y>9gS
z{!x*elIa_71_zaHK$vSXevM0ddMM{gU>(00aod}|W*DC023_R9l*CmPf8KLJZ)#r=
zhPJStE%GqA-8(Fw!%Q6SB)CPd)-KXTN?XtOJ5N>;Zd@^byv%jl3<nfH@)$T+6%4|p
zh*!DyRX(VdtD7i1<Vz1;chiy{do=-^X*p)+^`G7|5nNj-3j@Ru2C0!<i%lu0hR?7k
zZra8}HvyKCd!{o7vQ`bGU=3o=k+7Cr1(?l~F$%J$x6g;k{M>c~p#$%;{Bct5taf0c
zAgfgN!7;3%WjN@<U@wx;lVI=)p%A4^6T~lYq4(r9J@Q<RI+k?tQ<o%<tcmQCHI7Ab
zyb&7*<avfc`uK&<cS=Q;JWV$yiu`;)bsw7QX&7HOpo;CBhX#jXtz?j2Hpgw)E@!aY
zZ|cV$sVVJC?og<__kQw#^FufYBq`~xxAOC{hQug-3ok>*@?vtiySr2nfhJ~>%E&2G
z@VW;fF*EeSj~n63trG;(%#wsi%26Bp#@F8KCt=7e9kg!J)9B{xfS_h)ihEG7$4zJ|
zORe_Zhv~A;+KK07mifLN8{sk-92c0X?=PTSlW8{gP4D9Q!@BHLB=!9e5{{YoVHgmU
za}&tG?ItHeGTL0qm=LsxPb^QY?ov?tcvyOj!M>o|ap&;8o0x}4=2}3MFz;u7LU4f3
zO})dQu>Gt*HNPKY{+gm?q@AfCdQfIt^64MXISkr|%&y$?BA46itf>9J;ONRffup6!
zUfpI)Ct1KZ&eM|m?0m7FBH2R}ZsjV2$B-v2PupvnjxV(p#Xc5GaO=)Cm*CzF!7plq
z_u=LK1BvB6eVY2ma2UY<*h~cz!~2N`r?WKEb<p=i>9m)-3%;0jIwIrTyv0qdh+h$a
zNg{f;H{ZPOZ)bGOM=$8H;y3xi+VTfxI*1%lLUx9q8}If)kXq?bz_8WbH^Q{16{s7u
zVs+%MTW=Q|?g<dYgm6_Y`ofX$6A{D9y%BaD9JQchRIlFD#BL8U=uHD!sj!`R5<5a`
zE*CsyDct`B7YB@9)YH-9O~fo0ErzT@>xkX~C{^9!EQc+KiPd;27}e*?6n($CfPmkr
zoUZg(`gP~@kLtrcxy!Y+k;uuPJQTS<46wJBY$g5W|F%^B?%HS>82d19<u8EVimF1U
zjB)8_|I%wtfbF1UU`}hz&>xrE09o&b>4^wvmd%f)K32_!5Oc-&rLEE0yeEP~B=zC5
zuu5-cE<__VDGoV4Tbu0MxHP(k8WngZJSx2&AFVEf<TYotRKr`3ld5MP&|iQC-Ss{C
zrfsFdl0wX3U^J{nI?As;Q7LlG4zQqkOGvj|q~xK+KxiIp(37!-Bn=zrlaVD*7%@1M
z(N6yu<M(;9RiQ3}v;#Qm$T;&PT+2$UOpIW;MdpDo_jegWOO}LjV7Y@9DnyY3z~bgG
zp}x7P*hNVyS!yu~L9`ESq*?!|FYkUg!a7VTrH%Q6k_N`nLB#KK%uXb1;+GZ=d{<ZR
zthZ{E(|GRCv;hz>nDquEhWTg)A@bU%>I1M*dEx8814HBncM=M_xYgWhw0@_g&N_ib
zQ37iiw|npE4lG-*f@~v@9<X=SMWa8pY9`^QbJj?X4aElc3K?%Ss``86K+_8eprG?|
zACG2O%!Bp2m~~gy(IB4PNiH`4;b$PTyK2KEhX;wThQhkRIhjb?wFO#GP@am_a6bHx
z7Wj~@RfAZ**l{Dw8k>^TLDJwq?TydbAPS@zYf~vvvqhuT5d@748O-_CT7JLMo{0)i
zM=+=n14<j8a9di#7yuYn<<4*eOJjv&u*WFZvg4AA{Gj_}<+faVxvTS?TKqfUIFkCO
zLy^lG-&@uF#zU3AW;ki5uXq!XTys}PQ6g7p=Q;q`n86_4{FNhRFty{O-DjU^J_)+n
z$7LA|_HNt0+|k?5QXFrskOATc@jNep50xG*(&^3-#xa+Wab(<DuDDg!k(Z={r|ptd
z`Rv#`4`&IfeF@E{s(5>BBD035JAPm?M@vmO=14iab6N2A9kmR!ZSXFne;B9=yNl(T
z2md}+!i0UXdJjJs5Xs2JMsHPq3#F~-N^lX{(V>3TuF2dz#6P(AAV^-F*`c70-6(;U
z;-%nI+Si>8!W$T4{6Q)b5Hoj7#22|KL9VNQd8xO50hO0}NbdKS359`0k?230(Sf(=
zgN%QMoGJ{YT`IMhdt~7769%<%gBLUBTn{j~lX5QfmE2hUDSf3gW!hG|f(t3O30gd2
zQ)dn2j`L<Y&6WHf2gui_Rj|}acpGrNszwyR@_5c{5x#S5r1htyq`4ggA1F>uffHcU
zDjnnY{dIF>7Gews#ky_k6U<jK3ga9cPvwJ(0qadA8jBZ-@`B<t^Es`bzTKM|RqFM#
zs7`p8XXwyFoD4(ZJc}MmE5X6Rx6#JMNrCa*GU*ero}r$ejZ?*517DsuYm_MBp46QH
zt|I}T7SA0eS{^Ixa@a5b>`ARjhxg7FibK?tfSXT={k65t9K{u9OOUI3w<(uk3V|tB
z!&O+jmj*fr$GA0`t;y}8)zYx^zDuMWdzW6*lx&J#Qx5Q-2*c;qkfEY}lIvrHSvQi%
zgN4D#Y>-Nw^DjWa<ihM~^QkuXUqF}5a^Bp+5pKK9LRl$3L&(T5A)<XzD6QuesCmcq
z>}9#?sG&IHOpOYyg;1=#q9j*b-f@3+RS<R}wa^5H_SO2C=eC@f^B3SFh6_WmoNnna
z#u`T9qP~J)lg^2{sX9ah+XM{!3C2Y<;if@xT=8>EJ}-*7c-2HRPP%2hN>V!pn0t`4
zgT>AIm`lc5>pLccD$(>Fk|U0id5YC>#7Y3aKpi|7;~Qs}*1rI6PVmNxws}jqN_;Bu
zI4Ss4gsnbGBZ*ZMURD^PI{zT}7|fc_FZ}qi9-%G#q)dn_Z5@u_2sV*R&(o~%B3Rt2
z_0+e#(n#a%jP(V00`u@?ta6x-5I_)oN2WSr#m1Vwx=#{umVYIE#U&Gp(=B>lK#L!u
z_+1rXV}G`HFS}8K8kK20D?33N@50|zT;cG2V7p>mb83TCnMtD+0yWRr1`6ZV)A62L
zH-y$%!rV0$I|!4Nb@>&R34#LawU}k$S)?)sB*@dE*L837l(*T{%?ir_t|9lpPE3mJ
zWlZ1CU-a;k0ZvNe9J~&jm58ZiG7@pebrn4_>`esoK|+N2WMhOVy|VIo3+9xSGiKDD
z*5bliw=nkAI)@uGV#Ov|E#0elWgVjhLTj9T!RObZRNRZ~^Gz8&s^pd?iJbh@ta)b}
z2*Dt^Ula%}#6QEIx>xMB`Zd#>+i}6G{;h|T@;Y|shiX;2E`11`-fS_3tpgHAiXU{Z
z`)(=>H8AEEqdhO32%yZP6c5!g|D-=o(%<(=tA4!Mjo6DIKI4=xn7M1Q&R@q;0?cC1
zn{DL=BT|ywT@b?oD9J$Z6`hGhyDHLtFe5f>N^=M=D!~m7Et0AqEgr(la5S*8pDMZ~
zznPW!Rzw2BgRkGjz9y%edeP=k`*CQqLh%PqSJh86a!bT}2zbF7b8%g&+=&7Y21IH7
zqTk~kf3sWS4fc|Jime#dd3nkuzUN`9`q{x;Q;y-d@V5W=$ORdmM8B))3Ia*68jh8E
z#SZXI*@L<ccl|Buv*t_hC7wPXTD3LvDst`xtV9m$Q0%rvAc_y!nxbi%yB-B_m&su}
zByhjX<f1>s|8_!6AfNE@cFA|+1!JbXsd>>41<@$=x>l1ivOYC#rikl-`=PUuD+OtO
zA5<&`oZp9Pz~x_re#<KUbh$C5)+DV^?ni~RK=OnD;6+<2RX^Id@@?_L<+-_ToHD-V
zm6wM<dvs4C#&|6y|D3i>d@}P7l4`l1_}x1VnO=9<>9?|7MrK2W&6tw(&yP7cc8mkO
zEpEj9`9>ExM^I|L@|L?jrSXU;k^Mz$Z)a8hkv)%Rx3O@z$3M!s{Lh?>ix29V1Fh9P
zYs|qdthvn#ygYLdv~shA+*sa*aIjw<KNkO$=gkcOhe!|0><ZjcaI2sQ67bAn&Y<%8
zOw|@5KG<|wid@LcR)oXDC`9mU`*!4XEhhHMHXPB{T-Hi!auy$zg}Xw&p?$}J46B*z
zAjz&Euryqvr_BoL%CK*Y3kO)n3IywLBq0;1$-$T^Ngg$E&t4;m=)yL9p}MPJM{$tO
zf#o<FwWsG&nH6jCWI!T_>OqlNQB<(@`1AtLbAgLc!uSw54XXR2KUdRL@S@RoweT(7
z9s=U>Sa7;?=5nd>GR2_vFf#2pdD5&qcZy|T(k@?D&;_Nw9oxr3VW1AlJzhCCB`Ymr
zAC8F%gAl83T&x@3fg*WuP#esO6Vd0eU=&8pu#TN=lPfkV=y@5`ZjLqJjXz*A0x^>a
z2jN<ng4CfTB1A(aRnQ_4zKl$f4$CciPzTKFMwf3nhP@^N^I^)fzn2sEwLHHy6mqw%
z(fd-D=K~_nT3$^IWx^h-iQgb)Y;KW@<DG<X;7-p9UxKfU4rIMnqE8*wRQud%)g0B%
zIKx^Y1ccy>tr{oWxo{vaG3mp6%v@@sg+K@?G=nFj{fQHsH~<d+u0Bz=?C`WR(wikz
zM487Qn@U4==vzRX=SD=7+877tsJDKvbG)Hz9JNR=PTAC>_&B9lH4nZL`gKY`?Tf(S
zM*QgrK-A8Vfj6(1fCEf5SlU~o7I@~qI1}5Dv7-p%glvn@)3oKi5!;{?6l0Xual9A9
zyX&T)MP(QX3`$Of*kVq~{olXgZ&1S{)UoP03&6Pvs=u1BBSVKkQV9cighhAK9(<L9
zb2PIgR1d+Ts^zKBe%d1u5#B;p8A?ACE%P{I8%7lgF`iP>?<x#s)wPsF3i2-A+#@8<
z<~OS!Xav|kEfh{#F{`qh@m$RM>4%wqO8pV7O4)M}g6(@@dkZ|(8>2Q-b??~DmpW<)
z-?z&@zngW^Vd_{&r)8!c*N|+O+VR*f^<Mxf$nzilX^bbg@0~|LGZRo3(_)_sy5r5d
zYp7z6L%%UJuzr~DOZLyPb<uOrYyKFm;znSbmO~n!*g@_$y|-YO^i{Bk<U!GcuwF4b
zGX>1{*7>LyvLUS1I2Fhu-6!JohR||K*Wcm^&I>$Mme&M{71Hl}u4(Zird5-+2jFU$
z@>2e;z}%|}Nk=0`v>$a1CU{TRg-lDqG*Re~g7mr=IadHXLK8jT<xGgT1_6Z)`ku|a
zTcZ^Et&^OzMLz~1h`TFC|AvwttEq-cs3%Sn(SA{%Hsth<%skL?Du5>fOWtgzyt~oV
zYJ!*t0uM;G;6MU|)b9$}s!?{;D-qJUg@(=5$(jAz>b+l-vHv7W{hxCx-gb81*{Pi!
zyvb?aUHZ8!^^J`q5p&RIn$vpJ=VbC^%i8F-)6$fi0LHrvp(uzleDdjK>^)-&mcM{M
z4kNF~WF#n#)6S}ow?2$-<&{P82`=~xVfz7idp9V(`ERturln)+8go*8M2e#Vr=>S`
z%ddJmBn9@2Nx#~^8K}dCEVXgZ?7dsr4AyL4L|pG+N2(c2$mZR#JE|W)1@(F<3hZRq
z(de2d9e%&xF6~FTRi|9X(*O<e3!W%7GU<)Ybz`cOKt`3auM1@nkPzDDYtT0|nKOIS
zE#s1!Qwgj>9QJ!fQpuR4;lewECHVP@gw!&7(zL#$)wrM1D@c!bQZp&=F=NK`Q<_Kj
z_)TH}9(k&w1~EA61qiD*n94y~|D&{X6)6OR({)sF`^5vc!dOj<M@e@xa?!s9Q;TWM
zo7M3jyiAn%DJ0E!h7H54lh{M1^lFO3I##_hxYGvEKenDL6D)wnX0n^_(<<g5nFR3w
z<{E?=W@Cg*40otxQAobW247o)TY-J~1}b*lpgxggchHrheieLsK9=sA3GLekd4a_l
z@em*dO)vNEd@a8@sQ)Z4l|dh$RI&ru^H8up&$mfDAwHTXfU;kZI><)?Wia!BBCoiA
ztu8MS|E1b}a3po~YfbLPphBjgQAzuv56Z3Xp^?$Yf=i;4L4{D^FdU?RY+bI`S&p9u
zFF4aM$hbItHl(Seq(P%D_tCJ-_vW?o#kt<9c<T8T+u*55e9i;5)Cf1_?PxWU30#B5
zb-(*$Xj-MSPwBlUK~i-&rz#Dkeaf+H&+j#|12^cf8Ea)qf-Q;bB)5ZO?2kGUqk$A$
zEQv^3%Y-rvC}G9+Vgkh1las^BoBn%c8*x^qb9EJ?<!+=kI9jGdS4Y?P%Bu{1;o?74
z-u`vcT7L`F-d|!0#lo82R4&gyW8J)c6C-<Vm@dtK;dEJPNaco2HRipJV=7r^T7RLF
zrwk`%%~GR;+lnJ@?a4@b_P>J-{~AY@U60i(Y-}|6Fd!>;;z?bxnM{6GGj&>?4?Z;H
z4cWLX*7pawY7qN{>b^wFP<7cxPJT=`oI{27mnEAU)9ei}ocJ{5^Rc@}e`pbQmZ7bd
zuw~#)=By!-bnhn+Q3EZoD;QH`j%H}MXsIblIsy-}O+1B}F`>c1n%wJY?&(46nGR9L
zRlWR<dM?i199dz#tLl;86*$b@#3m)m4@vfrOsYx**3r@{Tg|(5OR0IvX;Ga6G=LDJ
ztdfT7N2J@rjT#rQVi<PuQU?VU=(pf5sxYuG$}eta=tj+jbaMv}aA5TAvfYIWJ{s>{
z-!bu-2<wTU60%NM@gn&))6_s@WM4G8safK1Tf+Y>=;wf=jp%KtNx1tTRjK`}s_nZ3
z@1O`7nnd9AWFt^_Y*ubnqu|o%FJSEHo@uz=&4t?1aOq#b-4Be>Um6KLMZyG#O?kXE
zrAZ-_^=}_ZjEd_ur(8R2=2lAx^#i%F9sL~jb{`7tp^sKpLKI6{obIrByLjQnIBj2Z
z4=p5tpmtKzTOZQOd9`!OUOb^0rE++u&t*;?BEH@B%<g9gDzp<Dq`0EH<I*`=Aj}sU
zBCf$UVNle_K>+OKb`&?qLfLp?1?3bZe(-np>VPTLm6K>B_!*Fc(cQKT(gu-`_SRQl
z;SY;VAJm4K;z8U3cFShZ`ISXih<8c*OSbGQObotql0Y7V$BC!gN_rqGRzJD-<9zLF
z2@rkEd*mjvX#AGLm*d6w_50%~iC>f8pog&_u|!P4OYs6Iw9;&}<h?h6Ux>_U`N9Q3
zuZAs2y&w(i+)o9F<&HFEZJcO55=cQzDk)Rq7aEes#T3zgIi>zUMjSxaKYITbkm_#Z
z1TJ(!aU0`AMYf4+qy{>8m??33H62ldD#5^d3K|u+?ow*96>^o<Cyi`F_u7eLnXtW8
z1?pgwfmpM=oP}G*+dVyJ_M|O_<@VrS)F@|5_+$3iAT+{dD;}gwYUy+C(uLQ-01g{^
z!-qQ~h?mR9mF1r2JuZh;jX0yh0<o4w%v?c39~j7Cvff~ul$l4)14s}799sV3Nby%V
zx7UVM2Dfec3n0YpsTiQxp+-KJTRocr7M}`)aBzfT$Ni?Nh6(lC5kGYEQ!ZLV;s6MU
zmbrqfC0LsQK<!IZ($gu+zW{)iYdMbF^p661rPp$rKMFqTkAPk21hcj1pWp2Q;K*~p
zpbP~N1^TW*$bfRP`&jm_dbvP(yJajY<1X9#t&mndJb&yYLwrZ-*-xRqMv;P<qZ)*%
zG}!$cGmjv@w>vI=Id1U|KN&@TNf!q+Uw7t#PsOwqG=b23`?x%cKgZZ=9mbfuCL6xR
zfxKwI#$?}C>)#~KM&=}ujb1K+I?X;-2jAM7v3~|S^*pXn?mpORadloftx_cV;}gzb
z6-HI%oJ}QI^!!@l{=P>eQnSycEbp`z?>2K!^OK6W(%3|Np~!9y?H}Wi-BfG!E3{qY
z1YW2*w5R}n?5u*^CBV7&NWQ-aWftmrc1v`AudL&j300gj9~+dS6Q@4$yNSo^<<?*d
z^ig{{8e7s|clyuEiJI_v33_Y|T{F5O$Fl(8KVI&9%{X=5UcY@k&e1%h0`sMd;vrtu
zeFc7+{Tt08J9eU%a&a!==5!ZRH2p<hXM4{=EGPo}v~x7%)!mJ>RmC5eyjL$K<yDfi
znHK|~GTm=$G~Ww&@j8BfJ9tHn)X=FQg@&fBurhtHRqhkR%+{h4r!$8vszw~G1aGYk
z>3?o*zWp5RcirG#f+EHwQ8x_1CRE|<`ai7oxbG-*NU%^rFK-yWB^m8`T>KX=uFiiA
z&JKR0ZYt1YrA!)6I5U%tVxp!$5}K&a6VDqJKdH23$tl<$b<nG>x22lPL5?khoYKZ!
z`8v64VwZwTkQFX`aSwskK&kcW!nV%|pB&Q53L&%PqP&YnkA~e_tBTjc;W-ga0P{c*
z!})7<!Sa4`f%AxVj!WGsj0&3unnXtWy%tfVpbunaE9Kk`SW_mi0EvvsMkZPURtS21
z#wId_NwJo!0&7JS-FN&0=>w(-|ESu3V5gxAk}Tom#k_z?KY0wp&QyRHR|9qC)k!3U
z6e$#<l2#rB$$)L^&7#i~E%lca5$({)<qyJA&vWDvk~fKues879^a9rD@d7L~vKMWk
zp*kG$_5vVo|5$DXU-~x@HCgUhTmZKtYP0y5^xG_UCWyyy4rASHp{8^>>3RUw*QYcB
z(;#6V^~osk^du~+iwHxuJk=t7L6WKfAvJu`_Jo&jpGDt&lr1P&4YOqM(g=CV%gGBb
zV_Hyp)8)gaE`!62<1c*s@G@*b3m|z>sN%g{+zXIu+pSp;8ZFXaLHm?z7PLdelZ$SP
zsEUBG#wDHj>x5d+x3P^(I@E#nrt|HzyGBgvm71!z;7ShT1S^J6t1H_F&28%s6+;15
zx<2i8YTwgax2sr{^}4e9I3U$McGCUTws2Qvb1OinO#QoO3(Iv_El$~D5-4s;NcbfV
z)juK!h+B2zC#kG|Ncf@%<^;e)5L`_1u|-vYi@+K6sNCoW{8^zL9Py_m;|tY^!!7_=
zENL-_Z<0IP8sH3q7fp9bGgSGEZR;MsS<D79_iyFm%bAxkxmHj;h{v}XCM81<+CDYx
zwK;xd9Q`JaNA_~%4mvdCD*$OofQwXd6W+JfOMfQNOQl>`*#Xw0sdEan7H=L?PWm=_
z93%(Bos=*ZyZcsmES8QPS71yO|7i9lS8Qmj&zo*J-%?DPe?&1Evh!m_2p>=40vx5v
z=cH%vD*@r#h;lp8lAa!uCA!C9i@TKucUx9sK0d{kMLeCMpgi^4novBLJ<q%KZ+8qf
zK|p->_^H`17MmQ8ef2jA9qsfRU-c`XB2!~swNP;mJf2;+j5Y_8Z*USg2?ojIto)d+
zaR$5L&O*W=>4sPQ#LBCnMF8sXle+}H<W6RmtvARX%#l$yhl-=BS02`FCeF^}KuLHM
zvKUEec!LD+h!NZHqDfy^g}_On*cybYQuz5Qn#y>Z&$_7PUS4%;Q<Zs5*kBwz<)^uW
zOvClk24^)MBgx^V`$1yy2gZ)g-h44_#;JpOI0W-j{mrN27vf2c(A%_Lx0vjITC7=r
zEkRTNZhGLV(6zbmw@nTAbXmEpJ9uC~{_bDE%uV7a#+-aZ29|+>7qh~K3YJCfkput|
zkqCHqQ1EhtOTm+SN@CX~yFh|Np=hxg?QBSV{Ap!jUp|+qxTG4Ro>3Z-Tca%Rr-Ito
za(u1Lc4=t>Y~<+Vv1F4h$zlIm1N7}S5`$UKuH<7TGB`+(8-2nuEmHT|Rwsd$_^Fst
z`o8j;i<m-Xi?SQVuowc|`K88#j^orM>EFyU_#Z?zj?7EKAsZReD^oAGBk_pwG{?{J
z>ns|jgBvzzrnK>hPN=B6oSe{9#MzBgE2_Vm*lJ@@9fiI$*M|y&lyC-m(GP``gCxdO
z$^h%`Uat-<>Qk6k(PjkpZ5+DYQ@fW(&C!E~CD2vbv#o|#*$p2%r9pIuCl`!KhJz#+
z8-3lv1odg60eiZYD@d*TOs$%QbEOK+E%BilC~5i!?-rE}wR|M-{m5U9Qwy;iA}E@H
z<y_@D#t=ekb(d@1YDNVbW-?%3{XKjpM6Qo+I2l2#LqRSl%}*((5Wwk(iZeOc?MWGG
zHA=BB5}I|TL2gtW6`Ta3;<3rfhi|x6rspbO&dc@#!0Qq>Iz#=I^SbitPP#^L0;Pzz
z!&HcqBH+{~yp)=%f%VY-{5O9ZG*Cy|P{|Ramd3X)=LMo=7&4v3j2bBUed+PC#b8Em
zSv0eYI%M}#w(Z$SG`Q0bb3P~+=x4xpE~el6=O6U~DWV-Jh#arXgA#1mqAAt)o<CGf
zumW2a6)c&M!DuRn(PVw~MnuEX`~CO@?J}UN=O1IRz7U+%<c9T|*eoK;$%!n?3zh!r
z^t=;7+oO^Rkc$neV7A|7yrOk}ynppdkihQ=Q1f0tgy29!t<nKUJv0f8cBo9*8vt#e
z7CVvb#zEy4V8427dV*&PKEd`SwlnSu!aE3>OA0DZ=)SgD!=ttQEowykgWUs%KX;05
z&ocNPEHf!RqL0SeR2jT)SXW^+bV#rp@D1avu!a50{fGUtgDQ>IZqA2s>Z@Q~<<I$V
zp8R^-2ih}wIq7MUvws+o*?H))Xj;~=i?<Cq;tmejr}QDN|3NF@L`ndOUY?GB_@q*;
z{sdkx^4NEA{noUoxMm8}_OCP8ZZ3W{3?XKU7WXP9r+d%{&6vh`7~UJYxac@Yn|N3B
zZo}$w#*sQVgNyLU`2t9@qg)2g_n@(?=)4%h*h)&4<o6A#A$~;q^Ojwu9BU4eNPQ&G
zQ-b#taEaG&sr?S_1Fg?|pm%7GBQf;IwX$RudLi$oHVrJjGJIbU-eNF6V6M4tEOvw%
zo82=OubAbgl&+}I{S+WMtQmYn7rTf(_z}6mlZ?Qb)hdfl`WtSgSe;JYGMS&F|1p!f
z<ue2|<By0ABEl|wq%=<ueZK?g_nMx&bb(pq&}ZK925JxAtekAPDP&g+vOc=$c_c{Y
z3Hl4LI=_wzFg1G@{pfAYYp~vOH2L_M`_}j2EAiy&rf-dC){md3kZt8tPvWsm^slT4
zI(n?sn#8G$*K)~!RbiDV#Hiim(7Ji+e0QOIlof=K2-wi$XX-=}6Fum+pYtv4=W*fc
z4}FPezAJtY*waw7TmnH#!IWZ1Il9i+(P4(*K1qLV4zETn>Zsi61BJP6^C}pvASt}$
zOAp`hFf6u;pa5`VrVMCfgukI}@9hRphvQK@LpH2v=1ateTwUvp8E$qCw7PdhN@}_8
zej{89qHbMu0TSTu9lWq#o<>5C3*4wsnt$@k<`U$U1vb134f~Tvzj0}`r3LZ6OL%m?
zf$c#?+C<w%^XKbZZC4jADt>uLzbgIRx1^!M(pAbF|1(MEZ8~`I8Jljb6c8=2r}@m6
zZz}wr2AHyWCC1QWRwMDk3VZ?rmCZO8qMwYaQ;@v)CCxZQYQSISO*1BGZ)-R6;Nz_U
z*wc9<*#9I|(RxEyY4d}Rd>6i(8^41QsL+*>P%jWI6z5RCUqW6Jt4|FWb5$;R)c*kZ
zQxVLyLN0Tp^2v!>0nfc7-;F(yMR9bCoOIt;L>5Z_!kb!V{a0yTT`t)qv%2zdTcbdh
zr(`ec+r!wN_@Ovm=wawUpJ7XE*QQmpx|0s_yVNKdV8hB^!^WH4E2`b|JC`x3?NM%v
zYA6Yb)G6vvAqft!sI`3;SJuaebMv!Fw6qzua?UMOtE*^ls!e<g?|@~_6kn9Z)o$Z8
zU9ROOAqk0{qJI|+S62gYMBllM@Q#}koS6j$9X{3&VK%k_Td9&0!Ec~n$;+PCQp97t
z=;P2|EkF5w@Pz|c{4|r0cLtb89yFUJH{~4aL@pL__`k-@-0YigevLA!8d^9S|KK?)
zR&Vu1&eyY4xvV+Q6<AdIuj3LALjK7HTyCfUOXjQ=xSR1Q5m~V_`R*~zHzq&uSvhA>
z3s&pDLL`?(IW)v%pOm<`%r(Q!chAuKiQ?G`KGlQbQj!j4yg!cLQ34^!Er+Piv#Egk
z)(CydKUcaDD3)nJe1>Aq^z9y2fGCxtPb4IW6PuDVu@fbd=tz^u2AfAW?`btJbW(L@
zTl7e~Aj<2SB}qSyg9$W@Uv*>qlxqf`;OQ1R=S8<6@_yK&VxI~Gxi{e$JzxNtJtcv=
zKxnkzq-BgZLegGA$&Sj>+)RNe*4rGv+^xb}oMxg`egA{~zE`tF4;mQXu`h>uk<Vmb
z;dC(l=bO}w*Y(qt`{A22U3Rx-Bs%MVV%vZEz5I@<yg$dVlk#@L7|HeJOaOz&t}-qm
zYkxdV)gV1O-QBrFsA4IuXr~o+wUxALxOni2W-bhk9fHx3LD3a$AXU49RyX;F&y&Zi
zI#`n#4km_oMl<=W?Y+PMQM1>uUlwQl^?zu4&!{HbZd*9i0MZ2o1*C))5NUygB2{|t
zgdz}nlhBJGU5KHBKn%V2-b56X-X)<Ik=~Ie54PvzJ>NLv?6c3izq8-{oi9I<JNLLo
ze%!gny4IR&&bi{8NMlgO-9D5;09wV7^`BH4-d-ZdpN7CgA}}>5=k*Vi!w}D7OG=Xb
zJdbjd5`hZ#XS4LBxOtU;04LYGlVmKkpE!TZg9i&FMe4v%k+7g&-R97l`WX_2W%4(S
z85iX8!6UW2nroj0wN+}mvU2r^b_%FJICoHqai^sthaqDs{@=5O_*VM^SP$fB+{S8`
zG`g0#8D$U;hUdH>U22`TxNp#fRoCSz_tM9!)##vdb55@-+&OuU?^A-Z{Ll}E+8!;d
zo!fUNCLWLwgT2Y#><9BdjGBa23`v@o0Z5U<&_fg3&2js8)VI;Q=-*)7&wUR@l}09&
z%8;*lH884**c98H#LuKwgm%!o*$~#nJ5K_p43zdww?6nTTN5Kkj4>i)EtUN|i=14U
zcBz!W^J*iFmaj1k8IiI$9>=V=$#i~S)^a7S?)$)ef_XY-$~PY-IH*o>z2qUT(e+)9
z4cFX!oLLt*$O$3|6uYWKIMh3V=+#8SW;}0AK~<?tqJf5!?r&iEhJ6dAxAXsGG2t@K
zi#YQm(@f#VMS2>bARyBl;;osS&+X4k8riEpkB$+q*El+2?RHx38mNIZD%caYP>$D}
z_1x6Z!B8Xot!7xzHh~VT>OCiF(8)<FG1u~1T`1}J2IC=OanHI|1TO79t}V!U$!TIV
z&PxhDw!JxU%c&esNJ{Dr*Okf)u&)8R;sv~k4iuR#pGOicfonUoZZlvzw5pu~L8n7A
z#lbmBKi?)P^^8eN17ND`UC+K%qX*vbHoMCRp~}>y1Krb=o-mmq{?lM%cZG6`ZSe;%
zjAyRO^o{Puwp^<ATqfL<d+9v^to6S(Yb>Q2Es<>+Ffyf3WS5u(yKNPW7@V!U%r6xt
zCT_okZM=_R+}J*;VQW40(!FI=lq^GJB+924LPD@DdBU~u1Ai)T!GtK&ia05>7&hg^
z>-7$jOH9Dx&nA?6!RXM-I*}8uVi@~-9kx(>kmx(d6C;VJ{Mj-`rbvr#A(HN7aOo)?
zwg}US<Gf$5w$<2b+Q>smG}+z+Qz2E-r>pog1|HK(ox>*l)lT5_-|_iY-I=)6;&_!f
z1yMG-(5(=@?yJru;g>n)u6lmB3Rpc^`^inHxl{n~uI`2H;mMy5%HJlhr~gp`vi_$^
z!2gKo!t8!4*K12XX;JI?NP(*633Oy&O=SFDR~ujZjGZpA^v}AKvl^+6c0ZMPg!lTl
zTaFRIe*z{AZ5qdg1Z{1v8?OBL)gD^9`X}szVw}c3p$DHn=*i%5i3~nx-(=^7%;Sp!
zd2W4=cbZb``}2JMiQLCK?#LI9{sPH7Dv^vblc!)BVKz@oL$y_`_Jz@9^9^*e-XJ4a
z|EerAWn>yOZ+ft4BOYr5fp(WGdufe7hsVD3^d@2@K^{oJscP<yAN4;NHdBg~aMHA0
zzlUGojbzFs26CREwzEDLUgnEh-ogd|*}}o^vWa=dsz|R-|C&uWrcB|QK1++4t6b~<
zT;<o150%5b<ReQIy(^$Nd_ZI1wv{x-bNes@W5qJQs6)M3rYmW$CKah}>mKsbKFwMz
zLxQF^SyaY&%ct<h>4OOPRd*G*1PYr_yY^+~@ut(cDjS!*C5g61Z~Z_i*8VlZf8b8$
zykga(f$qQNS*oWAxDD!6IZ((Fw2$^fQ)d|H`*v$p?(Un6b+PbP^CzI$UzZp>^W9zI
zF&jMki37b|J+Fw^sWP-aWss_U%FX{%H4fQHbZpRa2PT<mCZc;X+DM>cME<Kz9^gOO
z+(V&p!>l(U)y#n(=~``y#d--T`%O=6NKltVhfvqd>Tr^G=U@=smV5eosyf2h#xkMM
zi1q#V$7h$-vbhwdo{qZ(@E_DBWmBg-=@_$VT)Tr!k2-4&nWB@to_pIXa%w$^Y5cfV
zaNi~*J>sM+K)CuDe*>FcP$hp-2`E25OIX6|xeMFbld(jsh|<VT*S4S^Q+g9=5FeWt
zD|<wi?ZN&IFS0Fos>g%3I;Ulby+-^y7hAXbs={RTLoIk=D0eAVP3np~)*O`D?}xHv
z61HPV4xJg`7Ff~^cp~DZvo<i|U*3gFkf!OZH}U7_{b|lRJ&zEnZnR`C%FXhD{)ve0
zB6pbRdOOeIvv;Agc<S{+ACIG(a6z^bR7FGova0Ampu`Zl&-MqUwve=0c>g-A?P)WO
zY$V}yL2Y={q#u~$HY7Z$M`?W{o;^2j%Q8V#7%|$qK!Ur=QAZl7h_fn$&E(yARIG0(
zL5bi#tN$i#6888~_h5k=Z-zD|Nu_HDO+M)xgpwlW8mvaZVtk?|XC38W<#JyoYoet_
zfltY+C0v6<tUwILA~yp{nqL9&Jzq7G0wm{>w-l#hpLss=461Z>g4u^}p#jhVlU$AG
zE)DRxPK@{h8vwe=oGqo(brfo{Y?uZh+U1|C?Nxu>%cu9E%$$cNW~Lra4_XS9zNp)4
z2<;z3fY-&xm!nv(DutdJcHB@$BAC;xp(wpP^*bJvHr9k4+WTM1X#SD;`~xGuD=<n4
z8+_bm)+(7EeG~M9E!+c6`L;=O*~sZzi^Im=J1VNin{HG|H<I)$CX&YE3$qlV4D|V<
z)?7*wJSa2&;gylm?8<h9$DLZN8d@EMjWUueTMc5vUS3vu&X1>}%|0w_h;3)OG!b`&
zSEHIAz5bE4sL5ue7N@tNIi!O}FGm<J#o3uBx6Ncao~TvkJu_9KRT?UDUe%cy*MqXl
zR|n~&9U6@E*#&QzRGT4oxKgtdZXfoRB`UGXYnEkt87NK33nuO`B2lw$4>k9^z`g#*
z$#?DbbNnaOA-(2L<+F~K+?K&9*<6ZSDLUUbx~R4%D+FcS2A9IOzou&*5JGsci@|@e
z%l;+(+dqE8+=V^+jT=)^w2t=Wo)EgeL=C?s3VD~<<zO@VDF(ex9RWFXViv5fXO^jx
z*L$w)Ev?M*J)O>g29C?S^8y=rAC0L#XLSY%8fX%XRGtS@o3}&j4hH<(n#30N&6t#x
zIFV_m1|wogBP5(?dZ2ql1?+h#TGKW%e{D9<7ZuBOKXhL932QcE?!q8SX$bVax^Z#+
zB4Z!Z?K<q@J467(4?-$*`LhL=NQjo_A=D)@CyS;pl__e&k_&13V^C{v{5PXoX}fC2
z-XL2M_3jp|Q<Gx)Nxgk1VD*H`zP_o@`e0_wWs^ISzTTTd-hTt&Q*RqxbHoLCaBcde
zgHqHV66)xrVB8=CbT1Tj=wG$Lb*d0NmD{)9G)~YKgwB{9wKKHJF4*P!neyP^PNQOz
zGMI~s=$gBZ^`6xOxvi%2xC(U&<r8@^3<&4F>2vQ6Y`?uJIJD-dxodN!?1|nt-Me3s
zvI|b|%o)0};)JrNl;}Uxi2fgI%#PMQUS%q3oOYUObs-cMPhwxnMjX4TXD_K2c%aN%
zn!5bFgbk|1@7!-1{ti*C5oDFu^Yu{PIfvO9Q1HVJk4;$L4U+nk_Z(?eOq@uh-)%3e
zk=Taue=vp@c-Z8Kh+Wl&`D#n=n@I_#MYK&Fb+tXYxGN6IY~0swB}#R5>itM%_pk$4
z7Vt!1Pi5YaItPix`z37^l!G0J?7)<hGwQZaY~QY^>L~<PhwfsyI0=7mjB16ydnonU
zp(dNwCy?QTSBhbGxAB{8&D7(Rx8hj0u>SP|Sj>rTTq!2R?8!~WDWO3(*9qUAoVCL3
zy3gUheb}?TpKV2ODfmw!+y7Bo^Kb7ncYHl>sA3NM`8l=g#INA2Rx69WB*4+3l&ywV
z$P@AuVFtO^lPi%T375S;tg2TFWbLai`Rzkev<V4TbZJ3pvb2bkWSmIMR7>kce1lGe
zOG)MZzO2n1(7cpFv>22$4j8@FtpH@ClS|B%zcTJ8hM8;!k6u^r7$W_#^j)L`TMXU%
zMLcSjrUfT4!^Cwu{n6Zp=<p4h*(z=#ZF&kGl<I>03dJyi1L-VThJ$VSbDEk7Ti>xJ
znv$Pzdl$B|jHlK7UT&?L!^ws=v2eCjf&r-Ux0gc{H=|`u_H<~XBVNvxbZ7fw0|se0
zB`#XGf{yhr#6wt!^GgJWZ1FQ)P_TqqeJW!QzP_7ih0ZY9Lu)2@PL5t_z|2aFXr(mc
zSJffgI8<4^&fmCS3)_{;813HQtXIBGi;MIPKGyGg$M_CeL8U)Vi2D|fDz7<z%Nh)K
zk=E1Aby##Q)^Gik^i*xSvFPPYc2{VG%gMqi&gM^|EW1b1g3w=8UtH77iH7S3lFE49
zRL2%IJ%_h&5Ba(T?_B;lukn$mgOzrE`%HC@vgOaGlmms2ukZS2c0S9&VrRd;5a<<9
zAhG^d7E@W!n)K%QZvdbJoVJMly^|H-*uC|{acWIvz<J*7ad88Q$d{aG(Gc>ImZrww
zVOMU4MSh+N9PV3vXqTx>R*269ZoJK=a?N#fufBTs)8`2u`FGvqP&HEY@6gSQ`p}&v
zZOd_^daARl>jvrDQ4R5P#kdbDuBPbkLIk?&u<O6wPXFwhxu{$%DY)nc;KE#B<~h5x
zt@<Lgt?<}PWFDXUT=?J-4_dL5==$Y3B%5#bFy;Ad)y$T1_xH5?l6huR9%n>%eK6#7
zqb%gP6wI!51w%{;JzgrK(1|R7Ylk_LCBWo0n4RcvUVTUorM<$y-Nk!9)?0j4KQvzY
z$V-ykQ!QwchPDd_Te-a}WWplKn6b5Wjj931)nH&sykW*4=4=5{jvMY)?~{4&iQbH*
z77!?29#;G+Fo@Dy+_32`4pv^$IJ}%z00)uwcgbZKJwemiUbT8}m&f)biacIfLb?_&
z75QN=lDFur&3SA#ZPjVwcX?d(UVdN~imBH*U_ck*>yO;xy($U`PwT%RTl^m(@1J}E
zwF?h;OQejT*3^SmDV0BdL_pkCg<z{O`D!vG)^S0oZ>x)-rcDo#@V?6*s2_4~TboSs
z1mzJO6~9Q5Tf?zk*+%)I;-S|+@Th)zr7cWc=@{Bqsjo_HMq=TzcK5z4jrdk}_xEB!
zAj=@Mve?|&kBcjX2gnv`V$gZpZ+IWm1wt|F0fDf3OCbqWbv8dWwyF_%rJd@@(|koh
zz6Pw0M=rubU_F{ZV(RW9O=buIt>eYA%AUdo$pP!8m36?^nj`>|CO8t8HtFVw<^DQI
zvB^_ZW#_cg-_+Slg1wD;J=1TkK}$}wE4TATP}&$Q2%X|*Zu58*EMApLE*|a{Wv2a%
zTalK=Dv1Y7$s<3YOKN{UPINUDY8I~`No{v~+E{Bt!Z{K173~ff)SoMs>c@(*Ue%c=
zLwOH0yz{#%9(MVCNR~BSEn;c#;y&wR`sDd3y<{Jr&NkM<GQ_{HeGgO!5H0!Or*}jx
zl|qe7RWIXOt7zDTcLY+UpCo;wmB@FeJYwgeWsz<dGbxRn;tyJ@sG$xRuuPsFEs%)R
z-#1;Gd6HUa@Ye~n+EjSJ!S%!m?X!C^wQn%R6dD7>iaM>!Z1Nvs5<h>1rZVjBxjW?$
zMUwf0-_OXk%m(sgXRWei-$`HAZrzTLF#52Is<7cv6`=rqRj+_Y@9HYPndfH@e-bk5
z^T^V8VF?l{L51qzEg{JS{dR8=;_-4{&$JGL=xq`(v$X$^Ir<dmOm9Q(=_Ti90{)&y
ztRK*JV>l_~sUY2!B67_axE&;Ya>Bf>!!WN<R4uNZ*R8L&(O?Z%VbfB45cZ&`cArVO
zke293=*xiLMc5Fptd;U==G0CWKY{Ku*9tDTPQ%i+WjGX$YV5)I+8&mjm%&)gHutU1
zbw3Wbd0_OZSOQZi<#w7aH3HG!p&Xrcce|5Jj`$zv)P9%EX*$H>3QJ2F7ZN^pzl2c#
z4G`q~&Rs6N9Cc_gF0QAi39Zla@=V~u#w)6`U6l{9%I&gVSBw%D$4nFaRx)WtedMeC
zcdlc~<8CB#yd{6#c_CVFasmZbZgf^3iaaXUg<u4*Uy&c!Xcw>R^*S#SjeK9Wm<7&X
z{?_&RH5WtxT;3rH|JCyG|7|_^ciz@(-|Yyesqr36rJ*V1D5KG<7G*bxbHGOh8H$Q1
zMMjYgh4S#wnd)z9dSQ;UHpG7zZo{?urH+TKD&#*jY?@ajjb)>sm-K|oA*e2yExfO7
zU6@xu8^a?blVXQj;~HP5rg4_>E0^+TnzKGI$S3Hfb5+Vun0>mQ4BA+gT#A+n%Z8pE
zW_L<L>>c7mQCZ)MF9_y(pAk9-hzEkk$sSV2W50?t*hiSlw+%yvcsNnzpF>V@JQZ`L
z<V)bePdw?cx;welN{W4_`Bz%Ex0{T_mHMXK+a#M<zGJC_24o|LPhX3gwxkU{(ItKV
zylB`Cs1Ki5=ph1v&S|XLvFd}0+<Cp-{*|@$wZRAX^5?;v0HOlJwW^UF8$r`vAR1HO
zq)EN(=$!|UJQ2G4ctWx{5x_FUBR886SR2Ag5A7?ypI0@%b5$*USBcji>w-1zBpXM&
zP(`Nt=*tUsh#W*VE`6}Q9CoB$UhevZ)O`NOSJ>5;X!v`X4?8*FfSf$!R5gP~tFsrh
ziF$Q3=`ft?NT=8O*}GhPk*gO9x*iZh`;Us8m%X?CrdS<x$BP8;*iMg#jqSJKnRPXq
zVW>>qFX){eN5v1EabvAMO>MELNfoCrSjyL8a8wgUN|nC0=^Jb0w^qwSAH(OfI_XiI
z71;W^zcXh>y(yuOW3VvzwNgBDi+pJ|1r@*|6pOU1V@mkRu2m3<N=-UyQ_($f*7n@T
z3Ae%97d3`?FjbI7&*t9;8#VFmO_J!JxtGSdTJ3Y?^IU`j<$u?in08C0iKGtRlCG$%
z&MGOF8G7E`Il6Av%6u?jyOi;edm^Mfx5VIce%GPNx1;qe7(d#q4p$Y2OzHUz&f#q2
zX=-@Ee5b5uQbV5dh-a}gUvq2ri`er4V2o)UOA|5~YnmD}ZWIE-_sP0Szb~D-H!Kgn
zlj0r?W35^prKWYi6Fv-J+3oSX33)eyPiS0&RbK<C2h>%JTP^Jm4JOXFIAE-0^&PCX
zMsNs!mAv3*x@_pF>XsK-XpD~HZt6;?k=_h#JWiK!e@lye>Jc?vIjo+z`&N>ON0MYp
zt5o30t!WOKdz**ZvRRs17_4VqQLtlrp#4AeBmsoX`2XzT^{=c6xEI(&*=7}TcE9jr
z#^;yMit24q2Ol>M>?<Z{OHInlYbT49gr*Xvv*g1fOit)H*RR=Woc2Q#nu9D=Bt2z)
zm1<g627GQhoSfxD{sbHhmVxoVizBWPK-MSAMdDg#odI$&#^J3~d=Iv#jIpS4>V!^t
zeYqcDMp8(u<8Kd)j}xD5ZnOZ=2*Af}dR47)>p^ZPYA7t$K+(wReJB;$1}euEinOpO
z(V4nPf9wz%rJ~70%ar10!t+M_d9UD+dT6J&hGJRr<JNjvnjiq#2_3|Jj#qB;;aoFs
z=Ex(`P5C6N_zQ=~Az&FX1d;PekSM|S6Kp``;j(0JBM^B!SJ}1@c;5=x1@WbcvV23<
zqbu2)S5J%v{-C|=qxYG#5>|ik)o8$Lr-{I5xL7dp8`_4Nz<|lYJ9)i`ug1RE2Q!!o
zsQ2!9$*8kK=@}lS$sN2g{8hhb<5d(vL2z^{T_(+Rb(>V}6W{{#T$~h;uf*{VFQIoU
zW*=itlmVxSn^R`FBE70)oYYkO<aL}ep>vTNa!J7KH5el7)0Wah(ea|+Gnk&c)Q3N@
zMj?cvgGr(ca&$N2CjonIb)<8<&kPr0%o#0fXGk8rE{rFsJT+2pbLdCWsw;!@-lQw8
z-^S<9Mp*JNWI;i+optF9uY#y*Ct8vC7jY&6nviLwM@!^`2Il1yhnnFP%3P1`j;_0q
z6xEPqc6F>?JhJ_Bm1lMFXf$}I`DBCQ>r1si0iP1fo`v62wxP;2e1gRMZyVCKu|rx~
ztBXsWz{U%yhWHM;MDMdvNBeld*$C(TPwQOU$-&F-PCP0*VyveI#UYN5JToTrk-?xO
zMQYfTilN$aFJegk@JoM~^dq0`xZ^7C&x+3M?#~LwIr6H7M19wChdE|Fn@^v*8`$ni
z$sySAzGMsXNw7X|EFC;Cvyyeox$S*vYD@b!BqsjF8yIY)1{A#s-8XfCB9awhB|W_;
z(~sZ+I8U15a{d!gX{)<_193vC=v0r+vx&mBrPtmN(-?hOt4AFq$uHTQyQoglK2kXd
zC!2DG90YrhM^wG_0pxc1wopy@z*o%NFySctHu+x$w*C)oJ(7}{Q$)T5O`kC*iJL;d
zLcfuo)JXenY{6TPd}@=efEtFC5w?$t7JmcYh0d3;y?{XUikBOV?`lo%eZ~lfZ!d-Y
zxL4ITTm{0v>+iD-2RCcgei0M8s;V)|mhSv0met13!$!T6`g%{3xZs`yx{6vLQ#>UX
zsXKV^HC&BY5em$|yqz|em#+Bqsk@oU5R5$8fOVMANOs3iZdphE;^@a)k*~?%h<1T>
z+GZpr@jt1=eFzpYO!KIT-jN@3j=n@JD=)P2rIO<PV)k6OgIL{JttM&v=|fR1JGK3B
zd$KNN&L2k6eNnyq$?l~=eyZCsF7isJx<x#2SLIHMRW*HPTMI2wXRo!jOkA*|sZXh~
zGq#=lQ-O~Ebrs&Bi7J;(N-!#jmL)1V8k;uqXl-|`4$_dt0{#>zu5{$Ssb$p!<JLAa
zefkDEAk>7>A;LOhzn5EmJQoyR=-W@8ZSlCLRtTmhnyRYXB=$CITS$y#LS&rkqi=96
zaYRakiZveg>|@!WLiwn}s|rn*?|fudg89FzkGNIt<){b}U{p5L#D!B<yQjR5=CKvA
z34<3D#;iV;#4=pE2}=9z78O_7i|Dj7sBRqIL|+d$t*E~hr<e$kK_~LQIUOR2Igv4k
zaIInzpFo@?9a2r!XH2e1Znr7f%9;hr5dV9`2oBdOd78BL()-0!xrbE!;s9vgFky0r
z{P@f6d0pdib(7muBbRS2W$#q!8%SW)ku+nGvO9S)dX?P^3<X-roV|b$&@SYXCD%60
zVk>&O$NlEfki3`~BjY*xYU=T<?V``PY`z~nG?-DsSDX%l-ojH^rMG1ybO-_$=E6P<
zlf!^%`&Q}0#G8c!cH>{mK>naKy(~>>XGV>v)O1ljch+Fg1M7RarCjmulz$u|-kH7i
zaM0MkF+t>q8q<G(lCV}lZ<siG_Rr|aDG8V~cLQe_-MhB;89L6iMao+ZzG71fT3d&0
zop?kCE>T#b1rT~RfHWPzs$u>+j}P)m+Qo-JAsGC}qxt_~SOnPO^TJn^2iLAE6qPlC
zE^@it;=m6~H6F5>jUse6r{+f$i?4E-U2DYQW31c0@|98<q%_-SKc@F}Wi4Nd^CY_@
z{HoIn%LN@Kca<lszUA8djWH%IbS#PR$6?IO8(l4uuW4_giOlkEs6(hJFDINs^s%1L
ze(qq~XF61W4GOtB2tBcM*N)Y_?MIXX68r^zrzV{ne?&WbZ((#Oug2(M6kD$7i^V$Y
zL*4StCMnc`rR&Ff(?MCh;pj?UVG2Uu%Al8c61f*T8`~;$jktAH*Tu^pX7Q%fbj{k{
z^koCnEjQLED9QA>hmzNuzY4w)m)_Y#H`OtVY+caz-BZ6W2>R(zX7U9sg2W|*`+bgJ
zs-Pd6eC5N4zP|b#P&;qd%$i~lT7YGvZh(05ZvYIGOpt)q9ce3bs5)q_`vU4-Y49Y=
zH4-WI-7*%Eq8QL%6dkUFT`~I#1>`D{at+dco~zY&>RQq;=8PE(!%JP?qUWq<p{7@(
zfJm{AJaphfy#x`fYdKfl=dJM^yIl7T#UpKRjHEez-O>*<I~}*8DM_a90ZJJ1tp!-F
z5i@137`~<BtjBo=cwt_LEu#EVBWrr5%Aw~{U^S<Mz+|r=wHttALnU+c3e{eHy^h&-
zsR1n%@J+wWPc;LJZx>IRUKV-r4qRv1?KuoM0=I=QM@V_M6H1~&LbmlkJ*>?AZw7j%
z`A_mGt4&MFJ~pp9=-5%O&n-l_n8|go<Su8eEO)+5wlJoq1j>h*h&B4u(vhH=0=(|2
z#R6ahiq+Ljv@c=Imn?CQUj%Q()fKtpGd_Se70RsC7r8SX=-Ny>?4%L<s8m65Vog?g
zAxK+?K)&0Kj}n2N!|qOs($w0Gq}w%1y+fwXQbS)pOB(4D!koGEs#UFRrR-rwF;K4I
z(DLF}X^IMFwJS?X%I7Y-5Lp6P=>%!6=dD@JSF0|z$|5f>kj!hJu`)9%#d%4d%NFZM
z2Ay)BwR*G0UmlNr8&Hcu8lhOqRS1CWA4tEsW$<-{u0tBLVu1ulcX*s{zz|l8*2gfZ
z_1FXR71z@Jf46KD9F1EhE{*5m+6$cRM|1Wv*U5;4Uh<9(__sQU>(RAow7*KpNp~L<
z6;y=Gv&k{gdd3_z?MD5;kfmd(=VAkEO~iot`1j_omm~L=I+E|uS<&x|DlBexkr3Jm
zq`FLD^dj;a2ksx^ID_~vs?sNYb1f#tu#V1exs>zx{k7kuSwACTHU8_x)d3_iw)-L@
zS`2L((<w$uuTFX|v<PJ9m@i~juL;*cM{X-t+^7~#m03%*2&{b0>-nqVlM`r{3h+0;
z^|GS879vp`WxBwIfbAQz7@$cmYcAY;m8Hdo|9oS<yXXUjDQF?L8}mjFy<=TKIQXnD
zC27o3{eZWlrvM_k*R=9-{<V(S646%E)OmF_ZyTAP^pVCV5Ms~uj$2HmxNMqL`u>D>
zVJbVFkGnP$!T5Z*Ct)MLYB>kKHb~Kt!McE4*o@Kp3Hoa)b07@sBh;G74Nq%ji$bN(
z((@hBl(g9gc%h)PxJRtg2|17I>9ly7$0{G;f0(=NP;ASZ&?|uuWrpd4_j3theF#9W
zuIy??^dlgXA`(n6URo=`dmb;jZdu@^e0*e5LhG}2JuMwcz1Gua(zQ#Vt5~uvm#(tL
ziArWws{~S7G0ULSt(8LKxaV7e!_}L{cBS;6IA+y$+?+ZswSLv~vhPdXk%3Buge8Cw
z=|UWTmcyqVNrgz#>0s;!t_|iRy~s&u<$L@GH?lXb1{zbLy2%5wRsy@Q%odmAQp9-`
z?~FszN1jGwJRO|z%&s%Mgfm?&m~G~FYilnSnGvO}gnljE_~3;v$NS<Mxxu|GBjx1M
zpzjqnI^}$%n=t05$CT7OegHtFr{lfYNXuqv>&VAlwS14W9Co(rc^x#D^MQrcew`Xq
z4N8`J)=i#nHkmuHOuUJqtyNHvG-3Cr{x69yyyYTdCSOCgJ6v`aj+=4PI2HXx&N6aU
z*GM3x8VBNY_wRV<f_8Z6IK-xRN>Y4GSFaNxfyem%PVQ<UZz^@q`Dq7)+o@kV?V-bL
zb=te@N%v<G*Ci0LDq*V*=!AYN@Tv8j_Sflf?eCu52@XUD5^H(S5162K6-&j*vw6On
zL$<P}9-X;T(YvIZ`A7>{E%@DLAx)q2cGkZx_f5B1fX^D98k@;Ar7E!Dji1;09omS+
zWpH(;YTV+!xG+hLBm^!T9}l3R0j{y>yd}n3lbk-SkAQo0Zs2mU=#X%_2pMxqo2<bM
zwjXYM>Tbj*Pv4+7%N4XXJu+mn%mw^vgk0U?4O<9Yg#TOkuaRm#B5s$sq^5l&P9zlb
z)@L^>b~R0$ZyLaKrbX@pwCyz=Ze7HBdkS10c;E~*Su(4Sv)XS+I-$U*<nr2q3QnK?
zbfDaIO)ALE=Wc4h51m#wMJ}gogeg0$awmRugDVqwMM7Sl{`Ky)tN6_j`D-EhiJqL8
zAI5M`?pT8mky?Yzrkz%rPap+wBY$!p8P%S35eytk)oNH~N)-dxUBTfg5S${maxwF3
z=>*fud#Y$e6-gUhu~`d62ZeFOI@Id&XRtgW6<fy<(}{|QxDpNTCTA*n@m34TY*TJU
zp_<&>-=#~e6oh(hMKgmb!NCArCxeYOxf)y6AfWp7_vVR(2>AmmSwa5r3r2vPz#r$<
zq~;|~A5CljDXq$hx;O_c!wv7W0O8!3N6j?1I*>3@%B8#Kb=|gp594PJI%Wjl7^+z(
zOJE&deK8|xK-k36%Cgb64z|Y(Q$IE*i581Fa<{;$B?*6dPWSf9Odhp|^aSVnMbxyR
zPZM^-C8b<!c%<SJFYeI@<>MRT@(;nssidS;oGv(s_mZi_Cp&6t^bE1UHiO#o)G5Ul
z13Jv)^*14*J0Qn&jjL*zx)ClmGN2ETq_Scl_~QzlsQ(B_+g;*Y-FA;CoJTu{e<|cC
zadvbRumVY<<}v!x$Jq#M+`5?j6g<Nqj+7`%?P{0j!wpltcyKzU`Hqp&DVl2Py`Xn3
z6F~O%gzhtcIV?E~L@lJGtJ7qKiXeQQ6qF)0n=9YW?7sblBg^^wDkDTJoocf|E~q?l
zv;Auj@^jA<<#}Bs=M(!@wsaXHt@S68mNo;;=J}n~xVKhAJ3<^QE)BcJ?!KOp!<{Rl
zkGr6?tS6v1_oqk-*~zpKO%5Krne;?}v-FVar2M?MI2F6U0ZEWIYh}YkE8zSW_sS$r
zJ!zG3<!|z;-P?T3YsQa~5cdgj!g@RI-$OFCxVpHBX)`FvPBkR(<JY>A*~rOYNE_`_
z503D9orcjP?-h6mvPz`ie^B?yy*N|Ax+A{h?kkWxz%vD5pt_uSJ2j*iax(w@`Jb<2
zx0a2rDoFxmOsXGDkoq5i??ZYK$u|QQ7vxv$L++aV?CN-Y@(u)T->9+0F<Vfn1L|d2
z-@>HA<$G_`Bye*zW1iEg4?r=NYM7m7#53jhyVSH`SiE0Gx*nN+sx?(KB6?vL2*EWw
zGH6IuS6!7Rn;N$6CSwf+g;Dx`GkJ_v2Ku2q5ne>cRlIA(#`^a6x_ARm2xC21`MP9I
zon5VT{C%WYzP$ji*JQ}V5*&oIeuNL^IM*H#dg2ZA`-EuvuZ_a^XX|mX4!kVKmRVIq
zU#QZ<%=MVV#uJa++zf@Y)uA9k&Is;ylnYfU+6QTI$0_}UPpxl-=$9Uie_Bb+9iev|
z`tpa<O{r;C(F=VW=~k-gK17o+1$Y+eOqbJW+7L&<7OvmRi+i-UQI6pYN1GWv&O0qm
zSkcIbR>tJ~qoJy6@5za*Nf3Xin5ly-^}X;|5hhBBPG5=reI2qpT&5O8eOsGzbvN7s
z!!VkqbQhyC6MDLDm}s@m^AY<<8ArB_NCAzVwXI>aFYK4wX&5MuwPa6;Uq0L|by-vL
z4}x4Me-kz1L2NB(D^^ekuu0aSCMwS0Y>=DtJZiB96H@>|yEAV@usSFdxyIHHD$C2=
z_9Pp&jLGmMfcrwJgG97+RiBYNyIO$8GN&UNAlS~JcJ|T}i%yr>sW|{gOY!i&)XC56
zjb^ddWHskXAWd89sN0|_LG}*8cORZo8`m^%tCI!|c)l+#=5zp?pd<CUUs;ghTLcW=
zwJSczv(RMbG^n6Ki~qPc8Gi*LU@t=t0`wmw8NPo<b8XS2uzCRes<m$WZq=7cV3?FK
zT+rNrOIahwiDj6F)bwe)RZetC)S&1P@VxSQyPXpYNg=Sa%!e<BwTUOMkS!9nJ6DYP
z!Nt0g8^E#$E4HePX3Oh5(&uhk+klBg`Nhh4UzQj^ggc~pL-;|CXWEsToMr-Bk5n*@
z(hNf^OMe4I57bpmAh)JJi&v56(@936LJS>ly=%#*0AgORZYFa&*eS7cUXeGAr3bT=
z$tCN#Es{;#>OPBp6|GOYTc%ESUhG@?z`eMaHglrazj(B-N`8CO#y*6K7J1;FYgAGx
zS=Jb`k#kN4B>`Q_3khlo`DXPCgKJy1zAD@{k_qXVJgENM%%S!KX@<?uqv1{+6)5g{
z^-niF?(;0o{%T71poe(9#V1gQ01vN)-`kE)()T_SN_*#qJO0^e#r?6gSB4rZkwnvf
z17!RUeHcg@72?$IWeHG@fUxu3iCMdC$reJ~znfOONm7x&s{NIAMM6e)ikn~g(aYxR
zTBR@cQFjp<n*$TPqoYlA02Av}r6Gbhr-o?GA;$GxngLx;RI<KSnIF%n5<h<<JDH9>
zXmE)~EV?&lq|%8eHFy`35n={NUZ`}w=?L9L?dw`DM5X_#?%d7K>=)5z`Bi1@#N>j#
z&__p#rBW0L4eL8-j|mOo*!!!Q_>Loe6<OteE68jew!XzP-Xb3v5sZ|;r=d0uUqO=~
zMsmpnWACpknE&3F^83HimlD$Y9P=bHg!zrx!|n63?+Rg%I@WVtamV^3lVQ$_v}x3V
z*U7b&3YX20D=4auFgpjn{W+docJ6(4<-@nF#>5C1=T+a|0IDm)#2#-YV=%SdIUUnj
z;uxE-r<FXQ+QYY5?+&Gs^A>HtS^5aQa=OR}utDXoS!-oo@f&*w(k8auHuc1vIHeop
zN&tbw$CJiF?MQ&lMQ`jiK;*zEMV6%d1Dqj2GBLv-OeYGIYtjV{Axa#$dk(I-!?O@N
zXTHp5PG<4QPAX<-Q<v|q3|X-KN<9Erv_}S#zR&xxKMkM++qa-W@|^D?Ywnmk^SLfK
zfi_`#q0Ti~sGX~t0iTo4jgx-S6o34e=gB2fMah)sHpDs|*?g2Lw?$=$El!5GQbN<>
z^V4e`Pd>tihbPIDf+#U|apxbh#GcIoK5vP;paaUSY(MdiTu?sYv@3uMfBec@kzSP5
zh%Bfx|EU~xF2xdH{Q}$*A_Hb?>}B(_@BrtI-2prZtvopDAA9d9QUH}5=_6P<)rt+~
zie0h!S@LO2@y63={w`GZTA?5zk+0HXj}jBnuw;%TirPeFHAJBwsic&-Xd`!|kW+Q$
zSu4qfDO%PwmDZTA{gaO!6el&RYl$z=zun#}Akzh<X?G=UI9HTI_$6^sbC_|WjZDYN
zA--~B4Mt+ceH~$>ZlKVy>Gd7hgx};lLwDZD%@1F7gN^|FOC{^X<a%K<9ity?ucpA~
zVRuJUtBZ!e5JcGaibA6WGjYaIja&5MNVCmYxiD+CEj#~TZbBWEtfR3|_ON@;n*?u=
zBQ?HD30QhxI=`ioG>_l`NY%pw#29N*qpy0MIvl44QcAO|En|s4v<`5CD)4BGZuuKN
zas1zSd}i@k(CBY>_da=t3q5Xatz@*S)etA%>G*gN=6ZhW`wugtz!OJT9BIH5ys~bG
znWuv1s^;#UCN-2I5H#ykWJxeH7zkUak9ZdR4}GrtzRoT2CcrM|k2Cm;RvC<kxUtSs
zZ$n~;u)oV3)3DEXNiRVgU3q+0PcCSE;YKFww{?A#E9EdR7^s_e3OMHNg5~DYpQquF
z4=qj8I(dD-rCmnIVe$?>lf~ZN0S6eX({cBMI9*tmjaaL=N|(td0K-VHb_&;wFtt9H
z&`8g^ubf^|=3Vz9ILWK1#IFM0p36rMq|orbA%P5-n}$Q?zayCA`onz2!|7~=vg(C5
zOty{K(%SQCJznba2NGB|!Mlq8to!$0xJMMYv|hBZn-Vcw*s#JDv*$roo*a}C>Dz|W
zrr)UJ=DqbMQHb!eiy`ow%@x|bHSRB3&GR&DTWAI+4Q_mZXAwB(bOQl{t+wM610~OX
znr(@SHLCoXko{?U^xpDQ7(>do#Yg+U0rG-t!KWhXGi}B#Z*QwDyuUCPj)+*_`eT9<
zx}7g}S)1n1<g-p!rqR($Q*I#tO(?)rAu(;#Eo6#7z1N5RR^PU42rb;4d^(eMVnW^e
zF5ESM`|V31DV#;O#Yb;iAPG;Jsfb8<-nFd%9-z1Ov#;2vjK6UA_Cg87yidqQ42MV5
zb>ig?#!2Bz?O+%1=U@(j<9Si}{j<upb2cWXisY<3S518aGa(_iS8`RxIIKoAG;JcY
z_q@DeDPD%q{FR$l^NS8Y6q@Qg=)dNiDP5P+fd1Gi+K_uW?0Qq0FCG#722>g_%->Ww
z`#F-UWq&EN{^Z^6T-j{4Qe|OM2bAVxMvR$btIY=<mz8`PGYX>NUsVJ<4Zh;@*GHk;
zl=Jo^>^BfbErie5O(_247Ik1H$fi@iM{q$E*MT9TQhNb*SIT4yDj$M90jnwPa%8!~
z%=8@C8<lFE?la}uObc!nmfv`Oj+vBqx9srmHN?8~SR_)~lt8c91q{OChCOH149~@e
zr&P>T7z67y_l<rhXmxS!cJ%1FcrxlXEL<|z+VPEcFpM^FI&!+ejxQ@MRj3*Mn)_8-
z8)rqPIK<PYxb)merSmyjon`;VT-8W5M=wWN=DJ$i&?$(vxOGcS-O1BQle|W&a4Uek
zyzX)E!tEBZ$|`oMF6+z6!_j{L1OIC(+vYe^H_8%EXwZmUnYH>?Sey)=Zm_V{H6C0!
zG<azie4KD|%1u)lLWkjzO`MFdnhwu=IN*t~4kLveHQ_LSWxTd4%n`Tzgi~5>7XnSa
ziA0pCxc1jDhF#SF9_4fn&@WnWZ04#^<DdO(WJMOJ%#`}*>t&JV0=%o!_0?kIwF5u}
zBE&KK6Q*g%{FOg#2mmJ9Iyx*ebPkC^IZ68piX&`yCI3RUA}7R+L`#tq6zE=22MO~#
zH{FRx9yZo$7?XTwv28^DF9h?)dHb-kmpIFJJ0dlmB-5K%ht3~61k0dFwu2>CTUNEn
zK0=Q=;a~SW{P)GtKlx6(vR3})|Dgg9Ik1jfdkxnR6tIX+hAZZRCmvL|_=Iexv14JK
zO=>`Sw(_zH-DH)a?tBKV$y4a9PySWRv{3K1Tw}^oHsD4*WcEtg=Lmj2!oL)kA{_F4
z-_&e-ex@S+)Pud(bg=fm`Lv#Jvs87|)LglN*!hRArthcGo6*4Ms4$Jl=M|8*hCyx_
zRlX89AJR~KoM{dKnLmq7hglP=$;{>~BR{<lf5~j9xyPq-1`gmP86K(bCn?A-Fg9m|
z9F={@(@MRwuAWjk;913|FNS$I<{Mn|y&d)W#F|lNvrf<JODBNNSU--Y<*sc-G{taL
zm8!d{5xhEG`A+~suXWAR(T;M9Q9&zHZ)P|j0T7pniZVm(auJ5~VAG7M7DVVcKF6{4
z4v|1?Kw^yd<hH&kZkg4+d&V+3(Hl4*<(lr69Vs#KRKR*dW8xEn%=|;0Hjj8-A$jis
zc)ZdJU#KAAU!B<V)YbYyMENTVHTJOw^5><Dp$w?XPd#1-O$}K+nL<^)y@$sZ-s(z<
z1hRj7;819E?dpsAU;hBlnad~U2wT8}`rp!nz<&z7`QK8FE3*9!h%=~r$jj$ZI|wDA
z!F-`%27lHiq_xL)+*g)vlzhhhh9FEHt>B9l{Pi=w+o6~I6*EOiuK>uuA|S^bOeVi^
z7y^g+rx}r+OsL6MG471D8cVvYmU+Hf?5P1|1&za-+wcqtGM2*S*6JQ2N&Iq`cIXv^
z&`(KH91J5dG%L+$W-}7si7Ltfe&epu9#o?ksgr&o^<<df&{e)N#A-&u7DU`^YH3|A
z5y@<F7Nd{Sp|(3a))a)2)YAgL>cu9Nl=&O~)ynnS`=tEO_zMw!&iWZ7j)!^JJSfQb
zRk@Z0nq~O6P;9;D!9zWctIJ2!m!X8}-&SzupHW8tYxm3;i9@`Bg_3ydyi{q29tDIn
zS8$CI@uh*F@<zMK1EUaB_x<V5O+7SLx2D*2^|+Qq6(4><%|N}qVs{bJh!_6&H07Dc
zObB+rnbv4W1k3RN=PMvditVfiWj$ejroBV=ynilJ^J4LN>De;Gv+uG@>d(5_-vG0?
zx>0@E6GPO*eK%oU;n!)En{~EUDoKkE;$U?ew&ij0Np+)aCqHvHCg_x3P}w(+Y%1ir
zmy*)@4Y3Q3(aA+eWIZ*B*?Z(V{-O!HDOrqTCz_gTzc(*jg5Yr9k^+Z-ALkZBDcWgk
zL>Wh>xPTjAL+7!+DID|r_F<M^g>-fD$D1}qiYxl|R>Dk1U5{sTgq(PY5>n>NG~R}}
zq(6?q?3pd*_>_hRm!W#I9mavv&J70%w$q<}QrH$^(u98|3~gEsj`Ff?5p+j^pW|2*
zedHy~uPiER?!w|z-=0lPn^w*cPH^Fo2YMd*qbpxcLx5>!n=R?$tn=IQzXgzJDf_^2
zKhoq*0ymVzfPy`kA6j3%0hMbIE>2fBoNGMcbd6e}d(IY6S|>l$L{>|#3F{@vIPeIR
zwOABi-+dr>U0b8Voo^8P=evyQ>#m>|f)>ZMTeX9)qFpy^O257&rJBGTJzS&$uAVMq
z@~t2OEE2#eNNxbJ-~~@w3<t%8#+Xw7Otjzmxqc5XKFVhcjW!kFcTJ^<i?K`Il{~-o
zcz2yPa>kSDM@b^rl32PLWvyBqo>*~cuEFNP<!e@6b*247(`|m&UInc(vavv+0|U&*
z5WOp^wO+g+p{gUF*8C&V)^uw_l!Eg&n3VR}tMx4A*v3{X{Vw-XR71%McRU%q;P)$J
zEB)E#dx-(h5$Jg%(@y59fR3l<{8PMWc4+QTcsu}CzSjG;Ug4FzyJ4~3NT;~V&kw&t
z+jWw}^D4acnF~9(Zj;e^fT6;Y4HL-G+POh)*KFE1*PM^DZoByj9rKE<{0&I6I$I@r
zu%Zy>vZ8U!Lmw;1lW|pHP8NqYPC8+zlDYzmq)5~URp6NTQiSfTi~>6#4b0OP-y9U;
z7*;leM({soYiVyQg~8+H;5zBr+ox<{_wXI<e<nQ{&ZWdPPfs;T%e-w+sgzT7+Cjd2
zDGrY|-1iKMv_wiHK!T&bfnn6|IxQ`gQ+*_c>dcZReMHsk%omKK7dPeI2b$50y-@{|
zO?l=<=1v0VHQZM<WGZR+Gc&~VjDnuaN^U0u3G{4Eq;*%F{q$?n)K5I97HXyFra{eh
zj#+xDUCEFQL1OGCD|wd1d&7Z)r&-L#EJQDn7PnU!XE=c{PRs}7i+lZA>_cp3qP!i_
z4aQfsJokIFJU!^rG~{;5V`*Z_-TC^UwsbGt?-@u$M-;!c!`m#<8N2@pp4WYYd~h(l
zyKl|z<&GOQgKBu8rColm&AzgRz2{#x)A~XM`Sj`a<08Tto$38wi$4F}FWt<mWMv>l
zNl(l-P;vF;0EL0J?^>JGm(Yk6gWHh3j~G<5nD?fb)@=|`z8;xH^lEu{VI*J#5+K}8
zJQl5Z;(o-P(7qh4XAxm4+NZ-&To+=WQSrDC&j(z5N31l(+{(*5)89J!n=J~HMd3iQ
z&W2B2dY24fnSxK-JldX|%$tHgpL1v%Rc2$DRDX_(eQdvp4Xn|8#shoxZ1u$L9IkyY
zo?U4Q(jZ6X{Dh`iL?`*KJdr}3Rq`x*gf+tfjV?#N7BG;62PU}?8E>F;;};DKHLPxM
zZUr*lc2IFl2iCirw5uMbhh@NDQTy`wQpiIq!jqJ-et((6QF~9ur}!jNvTyL+Oak}C
zb!{}5^25T6G0ecOYK<q*x|<P21tPukzGp=>!KUlP`fX}8uX}C7X|yqLYvDmLszNeQ
zKD3O(0i6kb+~&K+IYZy;h`g)`A_k?{_mMYe1x*xbfd|DL)OGwO+~9-V;u<V$W?{@L
zCTGdHRfs$WQR`5I9LCf*Fk5LzIan0eO%H?wZkU+#Du_r#dr)bAr7n93I}$AL2%BF<
za~{UbHMO*jc14VQ&V40KihbaDW^CIG=_colghEW~vPafEK|Ur9Ez`r-<bo`z6(sx(
zMOB>>+95?2zLuJDUWX=It`ncZKFCygcg<oY76X9&FXbjD;Yo{j3Mj{{;j7d)O;rH_
zTx{6<9NQ~fjiGerOcz!Kc>Mj!*}%@RK<_r1*#*v8$MoP~mI3WoJyDWS>;5+u>9r+M
zXjF9A`GK+s<%zB_pupB`n}XfXbd13!jeSu#1awu-EY1Nq9FjbzB*^m_Xu5&Sx3g6x
z;}Z&41{Lm^MD>A0A826(+=FdndMIhqSf&63P=KoG#s@f4ccZ%A62VlYe6@c()7Hex
zK%M6cl-w~N8mHTwsdowbc70h!fAVu2xhG(VQ@bmW_bC!cXpU7V^aHNzZZkG?mgkkG
z63H1M!!4}Ko0=PRQs6pv;nNZ;rrRCD66wIQM~m&74+=t`r{5}pD);XO(p6RVdreD)
zPq;A|EU@8;$*+Gnlw<~hA&tX`-l3=PES`^Oy9e859Q*>JrK!OiwSDEIU&Q)AE-d@%
z-Rz&tj8QAyK+oe#_Cvk36)mZACJy`_r(6kug~Z%5&D?iOkS!&2YRSz}mZ2S+c*g*n
z+<#S@7{9N{x{3PowqOk1LJvK(kUHGenpjNayjO{DZo%vEYRB4?RD$r-eN0i81D>nP
zg9ql`6fC4>!kyO0Vg8G1G`@aXQVOy{n7c4LX62M_S4;kRr1~=#f6xnfk0&qHB!Kgq
zPzxRD(p&L2;nWQb{fT~eCi88@E){-<eyr)}5MNp89=U!wHsomvzBzKTm~9Xou#^fk
zeuZYVCxjX?!rT93Qv6r1;bwcJTIxv4wmQb+MD;fRucljzw?l8mXrqF^*9m!uSF_^p
z9uOVd#-Zcg$XM5;I3?`R*MolM{KGfoGFME+^QnD`PFoDdkHFfF?Z^}@VKR$Nn}vt1
z`K`<V=1mg*$8Gj_qIgBZL7yO;)LIGh37?Q98#%6<<3s~lwr2Nji0RQ(m^knEBNQmt
zMAlwAi5r~V`veIgr7U?$!@dgpTMj7a)W=Ls!~C69C-LT_)I@qx9|K8KDYg)3liQRI
zc9x5of&3^pf1*8(4N^6hB~j#;nX>dT2U;2Q*^{`_Y{3h2fpFP-ZOEYef=i3?ku_OT
z`LX4u&16b6snYZk=DyXGlR!Vrgs?9`+gFvF8_6OYdRd<Z(Yq6t=FyB44r`=~fDULd
zx$xy6R46#NHq#9IL=*ytIAUiiOh*=z2K$@9M<SbS=ivvL7g;$W2vQCjTaxknDu6s>
zetAiEN20x4<hmongrBw=d?8t(b##07<L+{sMMu24gQF)ltza;HG<c$VnA~qrUhvLu
zo`(;zCw=M4bWdwodGQqH_Z97bnI`v2tE8Db=bNun3m$0Hiu3<8E=*l#b3~rPjjN0E
zM%=f{O#xRMWsR>{%9n*AYsKk`E~{~H=2)>H5&ejq%BMms&)H}{-I&6v%2VqQ-Wz9t
zUQpCbvty^@TWX)8ml5~^uUOSCWeS3KYivfYb&CiPx+`a<`MAt_CV#K_t98@<fZ=nF
zSm)~d*;8yJGua{hayVUeyOp$+Xf^lST%aF6aj(vyJJ|r5g{@h?b|IP+5g&bMRf(A!
zak5|EB{PMxo;fPl(KUmlYCToVV$uesky}EwV2k{MpGnmyuYEZa=WBgdpy$Tb?5_ut
zB)e(RaidBG>WBJv+3Dd3BcS(hSN;un_$3W<M^1Wbex|`HXY>O7z4|)GsFMd)iC#S!
zoQt-zL!oOpUZONqWH*5%_3nQI7CX4S7I>QK@-&V{S_G^8;z2)&^uR<vOX<oi+Y%)5
z7S`AC*Y2cG6oKpq;dS{5&&wm`9GG=2toR)#Ph#=~CBIub3HZ<{!IPP90|4sAst2-r
zltE%^yg-(|(&a^U$>ovuRd#G70OCC%UH_*`4pU*i`gX;-R7`cvp{BBLlCH<CI#!2X
z#Z<PfZ-pi4U5a~s01TC?j74^1`a1Xn&H7t3(pHH52u%3W>ldN|8YI|caJSavz?%;s
zoKt`k004Yp3|Gd4CN>Kt?i8F#h{yEuRp54#AC5_ns;ev&`U3c8YH|T1y}#Vc$Qm;q
zt=SGnHGWNj0(K#2N@wArdJU!*3<z|7%nU&;JtD99t6C6;_ev0v@*&qRD(jjvyYEK|
z$VEo4MqOP)oMjcnbY4cvFs#grY%+>Yx(=Av#iI;yn6G*aZ{`aL{grnTmH)N|{r_?;
zWV^D%{J6Lo{QK1v`fSBq?Vedd)ou7aA2P-wPdicmea!ySqB8lQ5-SUYFdb%nF0SXe
zU$DtSJUJ>JdD);hC@tMjx%JX)D_Onl&TBXE7J}pp^D+6)wxcsOg}h0!P{a(3_@+se
zY7;gDl<f7IEq%Vej+iNIZK5iB110D;ntb`fM>TG(RBp{ICC)HXk}364!wM<t?aYZR
z-=>nQhuSysZ?hvu22>)P#j74pJ0j2z;A>5WkZKJLL{vdTXdsDz@Z53wa<<&r1ZXBC
z#&xQzPPVa7jFZSxj4Y6Ut+j$_QJk6nk@2b4ctyFNR$ygiOv4Are{1hd!=d2XI6jt4
z*_S3m_EF44Av;MFMhqiknN%42z6+71Qkj`-Lrsh&`;u+!6f$HlTe8Jt9YU5NOT6=Z
zd#>yK_P+1O_uKh&KAh`1_qp%;-2dNy7~97N)*O-Xq4xk$-xDy{Y|`nam)$iVl%Eh<
zfzcdV_fU=&TXSM30J1wVM^vd#VuyQ@*+Ls!W3XFbr?XMFvsj!{;xMhPxtASQ+sYP2
zZ1BBpc2_BrB*t2CKsN+xh)>qIA&KQO?60glqq6aOO_Oigjn-adnsCwF8K{Iqxq|q!
zUMCHBQ9vy~O&$u|IifHiVyL^8q2kL3+Y6>{0eDB~!SKx)I%V3BAN{m%BJQ(JI{v(H
z6ieZ_@p;2zJBHWu{+;C|V{LDOodL8(5?zlWpxB%e<&XXgLr&B?Wg5(Jj?EV3e<f3O
zPpx2#bQePUo|E~P@;W;}Q|NH8fq_$hG|ThT5p-2$y_=+uXolAM@|1Va3`)KR!{ED0
zLRBObW^VKi*v}84+5(Cr>Q+f@C#S&@JBO+f>VVKFQ>ZE(6pT!vJ5W5;O@3oxKn5Vk
zqhe4!^56lgVQ#EIUL%96Z~D7tcx1Esh^>{-Vw&~a`cmEtDM-x@rfcpN)5E6*z&lFM
z5}O?4=(J>tD-Hnx%TN$%V5xocz6Z!GC+$jH6QfUf5xJ+v)GDc8t-x`ZsoOx7ETX!K
zY$}tVgtMxfx2njp>?#HPO4NXgiM&$^wFmk5MqVP~UMr)XL-IIeYArdQx~4_puUe{J
zYtPWNAGA;Go?wHBoL+&}ref?^^QEu}Saan)n8fVn@13GOukA+#1_E7L)d@N`drRxj
zFb60(ZoD^~eco9hG{UjKF911PC=WURe6KIikFwmm-gRbjD&(5mB){KXCNg95FfsE-
zm{`Vp*Z1?b4n)Dq;8XQkGL<&mr&%?&iO)MQcdDdT2XiW17BzwazG-}Y==At2c2VQ?
zL+)n1vEhB&82mKhJZPnzVa(jGTe2>^uEXsjMA$FX(#SC|NLM%SZP3|(4SY45d@H~o
z?|K`yr3R^i7UfO=v*y1B?J&GFSJ<|cc$#C&x2gEm?Lu4jLUTB*YDK|y5bgzGpZooa
zN-^`zt}i;07)9&*8Q(@e<olU$`MkFb*NwFeo3-0AWr{e-SpQn|%O)F(i>RgEUFvc`
z8!Icl_4PkWcm7V3%=E?a+MhJd1MVRDmpMK=az>?;yu;!jmYkN7`=GBoS}opyxABjm
zPB+6;pDH=mW3w>Ya#D#5r8pH!<#IAF$0m)LCt<D6_7c-HN$BxIOE~bx<tl*$!<&ZD
zXC}1=7Kr4A(zKUy6zf+y@-Z~3&HDvfzT?CFs3R8>Q+2jyKpBGpJENLL(pYWs40OJ5
zp3gdH@TZ5yqx{-tbT>(ttE15zbt|C6F-&3pTawAA%C#LKbMa`QM23Um3-Q6C={o&v
zLUN2568P4zlyMJiDs!k4QJBj@t>|g2(;sLaim6bMA)vc?K-AEgROF{V0+HSs7)|OE
z=)aD(WaB4pUy20m6|IR+#b2crE!yCt6Y?fTB9B~^<9~MbRIxjYw3r>e%=hPYdsBCR
zEhATD4L~}jTPQLIPjRrd#Bz0H$-oQHGx0tIKBM^{QNwOj_Ot9B$VaIoSGMFOG0b#P
zaU)OC{!nG8eN|vth0VM1pJ{kRI16-BYcrukJ;oIDIR_9k1WL|V3luaU25A`d!wc+W
zG7p@y@<D2=i7uw0KCrV7kV$-me*FbwGx*A>t5TV{U7c=?Vv<gsfcZZ#gV;A64~6w$
z<y}hlqYor}dK*)-7@w7u3-a6*d$ov+X`vNWm8Bf&M`#xhgm|2-tP@y^A@&BlI$qWT
z-pDW%KQFIrkPk_ru9-xhxxhw}ClPFxH0hQ(xR{vFY7~#%$_rO*?{+}?^}=CKW+6W6
zG2?lj#sYqyTPwZ*ZkZ;GTgr?9d|2F7DU+Uw-bQ~n(V&uBAJcWvngN;h%S{%`i6~jp
zyR(XMc%TxU*gYV!!&|ph%FbrdD|p`*?a;YU&_c)aN1q3OZxyC4hCORb{Y;}^q-%Dd
z3cZ8*X%ivG)_G?bz|0OiAl5Yjj&I*er5yzT!K-*e6NvC4EHV9)&<s5(Xb#)kFcKq7
zf(2bYL?SR$<KokyJr%X$;<3KBWqY>@Em9#L0et)Bl8KX7q^iZRQ7Lqy*4B1_@J~g8
z5NEhS&Q#1+OrJr)8I;KU@G}MZs6=%D?45DNJ?;|E$)#_KHSLP6COxEl0g7uUCYf^y
zfQ`J~Q65Ma&k>~`qBIj$V4FM#WDq=jiUREaiH~-Ft08T&kpwAV<~12A^|MT6S)ndL
z%8U*IK-O<wue*7MSWUKB<-Bvh@MW3+;_&|aD>7t}e!2hqGUlIqAk0v;vdS+3m9=Jt
zv)eR~T_9uThDmo-K`92E6XZ~r=Zd=-tcwuL-qkw+e0Y37BiwSS2((0osS-~B<B^ZK
zwoU+o$Wp>`zDO+HUdL(&CFAK>#b5QJCa^{SAT4aY@7~b*xg$BwV(DX1<jOVo)LiC5
zI%8s#gt54~NGh&rBY>+Nv%Xc656hJ^(4}%jd0Z2Tg~gqh?eJ}|tWa;g@IG)mHwoE5
zbWB+E?Ah~?KPrB2jpG;QGg<f(!QzHZnjGKCV^!cX#1?NMe328YhqmT*v0<Dl#eNal
zCAuY6O%Iwlaa8Pab_2$e4SC7D_WcvpzJN>0Gc$|`lhW%w*Qg0A=WrqQb3XNLwL+ou
zL4Vv)Gm)PbWnPr~H~LVF>0soyPm;>vb?ev!Q7y6p&llQXW!AZ+UGXwSy<Afq69>cm
zIB!pdOIi_jjozvPOIielEb!1St~eA2-?n<9=R15H@I!jN!#@h0H$B<9BokL9`j{S?
z%*e$4+i#G$OozF62HMimA!>jN+%2DkUMonRFphs+4#m6@<FhDx+Dz70koIa3qKgHh
z2N%CvOWqrs=Zoe-=#LUVb&elu6vt{T)RVr>VELkbO`sjB`PV2pme2EyyN-<V#X{+B
z-B0fkWY1lZGv-Eg`Ycfg6tncME-ttm(|{LCoN^p022}ImB0!-CC9Mo4?0|%3l$EmH
zHue?WECUww?gF8byAtT%G2)T2P0K0N9XmugdYdjj8m+Ll|Gb)Z_;p@Xj{Ug0n0?%%
z90<36;ADIHXeD|@Fmy!&yeM1uBU)SS_&zOX*(^HFAHDy=5OHwYv*kW{Y_@xz_@yxe
zf#+>=dHdap_b~&Kx`z&^nY}lMN32RQA4|pjMl~oBP@I8#aYUQi3e$Ebc9G-iE6G_Q
zlnWHv6H8x!Zsv8meZm|*J!Gdi>7Hx0@T#H|U1XGB%AJ2}J8acN&I|2=G>96Aa1{Yf
z%3F9`380E=VRF3dP5>oyA8km`#2?hJ)7c9Gz-hbgS0-h%@U~Y~W=j?kb|Lpi9h|sX
zT<~K2ve=R0;weRKYljyvmW$N1#)Pj6lBVZO>S97Q4MVxdnE57Hs&TA{M@4SC>rznk
z?p4v#zQclkS$2u)#{}UiT4V-T0Y(bq<SL_&^{{>=e;D$TI00A-^0e~WnXG>MRZ<8;
z!BqWASyL4(D~xWc2<mCa4x@g*Cl4vUbdC|@%m@78(*2F%&^x9e+9h;P07k#`m*`V0
z&^xi$)epPkl7h-5V%IuvAY5YPwEh&`7QSPgI&Ukk`M9LxI@8cL1syh}cG*;^-DVFc
z7H%*Tp-3Idl7>Pvz2M^%<rg5hu)Y#+`VzpFJ{NQXFfGR^=fsJ7MQf^T&7h30Ta@eY
z8jNd2>>+>ei4tz<>nNdA6ra^kdfVR-oGaC5nfj@oHgedk)214CCgD}KAN+u8)oK*|
z<9klhBk@mZn~|0i(|P1yJ8m5c#C`XDPMBX@kmAj-T+<@>H5TaBL^t<>Z)a$4-6s10
ztcgb^3=!ewzc%qQdIR0iU?RH<%UEYr;c6fC@0^v$b()}Nj==Eks|9ykFA8bCKBj_;
zOUuSeRPh2#n5anL5^=$$Oj#f|G#FVs>9Krc-KCGNm45uK0rHqpz6F$i3-n1nsbE9J
aA9eQM<u(6|ivM5#Z(o`J;5>mR-~R=w^az3g

literal 0
HcmV?d00001

diff --git a/examples/README b/examples/README
index 8d0d7cf1c9..ac3fa4ecc2 100644
--- a/examples/README
+++ b/examples/README
@@ -86,6 +86,7 @@ mc:       MC package models: GCMC, Widom, fix mol/swap
 mdi:      use of the MDI package and MolSSI MDI code coupling library
 meam:     MEAM test for SiC and shear (same as shear examples)
 melt:     rapid melt of 3d LJ system
+mesh:     create_atoms mesh command
 micelle:  self-assembly of small lipid-like molecules into 2d bilayers
 min:      energy minimization of 2d LJ melt
 mliap:    examples for using several bundled MLIAP potentials
diff --git a/examples/mesh/in.marble_race b/examples/mesh/in.marble_race
new file mode 100644
index 0000000000..42eecc856b
--- /dev/null
+++ b/examples/mesh/in.marble_race
@@ -0,0 +1,43 @@
+
+units real
+lattice         sc 5.0
+region          box block -110 60 -30 220 -90 130  units box
+create_box      2 box
+region          particles cylinder y 0 -30 47 130 190 units box
+create_atoms    1 region particles
+region          lid cylinder y 0 -30 47 190 200 units box
+group           mobile type 1
+
+create_atoms    2 mesh race_track.stl units box
+group           mesh type 2
+
+mass * 39.95
+pair_style lj/cut 8.76
+pair_coeff      1 1 0.2339 3.504
+pair_coeff      1 2 0.2339 7.008 $(7.008*2^(1.0/6.0))
+pair_coeff      2 2 0.0    1.0
+
+balance 1.1 shift xyz 10 1.01
+
+neigh_modify exclude type 2 2
+timestep 1.0
+
+fix track mesh setforce 0.0 0.0 0.0
+fix pull mobile addforce 0.0 -0.05 0.0 region particles
+fix dir  mobile oneway 10 lid -y
+fix move mobile nve
+fix load all balance 1000 1.1 shift xyz 10 1.01 weight neigh 0.5 weight group 2 mesh 0.1 mobile 1.0
+
+minimize 0.0 0.0 1000 1000
+
+reset_timestep 0 time 0.0
+
+velocity mobile create 150.0 54634234
+compute ptemp mobile temp
+thermo_modify temp ptemp
+
+thermo 1000
+
+# dump 1 all atom 1000 race.lammpstrj
+run 10000
+
diff --git a/examples/mesh/in.mesh_box b/examples/mesh/in.mesh_box
new file mode 100644
index 0000000000..28ce2031f9
--- /dev/null
+++ b/examples/mesh/in.mesh_box
@@ -0,0 +1,48 @@
+
+units real
+atom_style hybrid sphere bond
+
+lattice         sc 5.0
+region          box block 50 250 50 250 50 250 units box
+create_box      2 box
+region particles block 110 190 110 190 110 190 units box
+create_atoms    1 region particles
+region          lid block 100 110 50 250 50 250 units box
+group           mobile type 1
+set type 1 diameter 7.0
+
+create_atoms    2 mesh open_box.stl units box radthresh 4.0
+group           mesh type 2
+
+pair_style lj/cut 8.76
+pair_coeff      1 1 0.2339 3.504
+pair_coeff      1 2 0.2339 7.008 $(7.008*2^(1.0/6.0))
+pair_coeff      2 2 0.0    1.0
+mass *  39.95
+
+neigh_modify exclude type 2 2
+
+timestep 1.0
+
+run 0 post no
+
+fix dir  mobile oneway 10 lid x
+fix move mobile nve
+fix load all balance 1000 1.1 shift xyz 10 1.01 weight neigh 0.5 weight group 2 mesh 0.1 mobile 1.0
+fix rot  mesh move rotate 150.0 150.0 150.0 1.0 0.0 0.0 500000.0 units box
+
+reset_timestep 0 time 0.0
+velocity mobile create 150.0 54634234
+
+compute ptemp mobile temp
+thermo_modify temp ptemp
+
+thermo 200
+compute ke all ke/atom
+#dump 2 all movie 200 mesh.mkv c_ke radius size 960 1440 zoom 1.5 box no 0.0 view 120 180
+#dump_modify 2 bitrate 4000 framerate 12 color orange 1.0 0.5 0.0 amap min max cf 0.0 6 min blue 0.1 fuchsia 0.2 red 0.4 orange 0.6 yellow max white
+#dump 1 all custom 500 open_box.lammpstrj id type mol x y z vx vy vz
+
+
+run 5000
+
diff --git a/examples/mesh/log.8May22.marble_race.g++.1 b/examples/mesh/log.8May22.marble_race.g++.1
new file mode 100644
index 0000000000..b5f5d1c2d0
--- /dev/null
+++ b/examples/mesh/log.8May22.marble_race.g++.1
@@ -0,0 +1,156 @@
+LAMMPS (4 May 2022)
+OMP_NUM_THREADS environment is not set. Defaulting to 1 thread. (src/comm.cpp:98)
+  using 1 OpenMP thread(s) per MPI task
+
+units real
+lattice         sc 5.0
+Lattice spacing in x,y,z = 5 5 5
+region          box block -110 60 -30 220 -90 130  units box
+create_box      2 box
+Created orthogonal box = (-110 -30 -90) to (60 220 130)
+  1 by 1 by 1 MPI processor grid
+region          particles cylinder y 0 -30 47 130 190 units box
+create_atoms    1 region particles
+Created 3601 atoms
+  using lattice units in orthogonal box = (-110 -30 -90) to (60 220 130)
+  create_atoms CPU = 0.001 seconds
+region          lid cylinder y 0 -30 47 190 200 units box
+group           mobile type 1
+3601 atoms in group mobile
+
+create_atoms    2 mesh race_track.stl units box
+Reading STL object Georgs Murmelbahn from file race_track.stl
+  read 9472 triangles with 1.00 atoms per triangle
+Created 9472 atoms
+  using box units in orthogonal box = (-110 -30 -90) to (60 220 130)
+  create_atoms CPU = 0.040 seconds
+group           mesh type 2
+9472 atoms in group mesh
+
+mass * 39.95
+pair_style lj/cut 8.76
+pair_coeff      1 1 0.2339 3.504
+pair_coeff      1 2 0.2339 7.008 $(7.008*2^(1.0/6.0))
+pair_coeff      1 2 0.2339 7.008 7.8662140345520858986
+pair_coeff      2 2 0.0    1.0
+
+balance 1.1 shift xyz 10 1.01
+Balancing ...
+Generated 0 of 1 mixed pair_coeff terms from geometric mixing rule
+Neighbor list info ...
+  update every 1 steps, delay 10 steps, check yes
+  max neighbors/atom: 2000, page size: 100000
+  master list distance cutoff = 10.76
+  ghost atom cutoff = 10.76
+  binsize = 5.38, bins = 32 47 41
+  1 neighbor lists, perpetual/occasional/extra = 1 0 0
+  (1) pair lj/cut, perpetual
+      attributes: half, newton on
+      pair build: half/bin/atomonly/newton
+      stencil: half/bin/3d
+      bin: standard
+
+neigh_modify exclude type 2 2
+timestep 1.0
+
+fix track mesh setforce 0.0 0.0 0.0
+fix pull mobile addforce 0.0 -0.05 0.0 region particles
+fix dir  mobile oneway 10 lid -y
+fix move mobile nve
+fix load all balance 1000 1.1 shift xyz 10 1.01 weight neigh 0.5 weight group 2 mesh 0.1 mobile 1.0
+
+minimize 0.0 0.0 1000 1000
+WARNING: Using 'neigh_modify every 1 delay 0 check yes' setting during minimization (src/min.cpp:187)
+Generated 0 of 1 mixed pair_coeff terms from geometric mixing rule
+Per MPI rank memory allocation (min/avg/max) = 4.803 | 4.803 | 4.803 Mbytes
+   Step          Temp          E_pair         E_mol          TotEng         Press          Volume    
+         0   0              737062.81      0              737062.81      21986.781      9350000      
+        67   0             -2063.91        0             -2063.91       -5.0227698      9350000      
+Loop time of 0.518872 on 1 procs for 67 steps with 13073 atoms
+
+99.8% CPU use with 1 MPI tasks x 1 OpenMP threads
+
+Minimization stats:
+  Stopping criterion = linesearch alpha is zero
+  Energy initial, next-to-last, final = 
+      737062.806250078  -2063.90998808136  -2063.90998808136
+  Force two-norm initial, final = 689296.27 22.226599
+  Force max component initial, final = 336546.89 0.90593277
+  Final line search alpha, max atom move = 1.2850327e-11 1.1641532e-11
+  Iterations, force evaluations = 67 393
+
+MPI task timing breakdown:
+Section |  min time  |  avg time  |  max time  |%varavg| %total
+---------------------------------------------------------------
+Pair    | 0.38371    | 0.38371    | 0.38371    |   0.0 | 73.95
+Neigh   | 0.038019   | 0.038019   | 0.038019   |   0.0 |  7.33
+Comm    | 0.001586   | 0.001586   | 0.001586   |   0.0 |  0.31
+Output  | 0          | 0          | 0          |   0.0 |  0.00
+Modify  | 0.033492   | 0.033492   | 0.033492   |   0.0 |  6.45
+Other   |            | 0.06206    |            |       | 11.96
+
+Nlocal:          13073 ave       13073 max       13073 min
+Histogram: 1 0 0 0 0 0 0 0 0 0
+Nghost:            590 ave         590 max         590 min
+Histogram: 1 0 0 0 0 0 0 0 0 0
+Neighs:          66778 ave       66778 max       66778 min
+Histogram: 1 0 0 0 0 0 0 0 0 0
+
+Total # of neighbors = 66778
+Ave neighs/atom = 5.1080854
+Neighbor list builds = 5
+Dangerous builds = 0
+
+reset_timestep 0 time 0.0
+
+velocity mobile create 150.0 54634234
+compute ptemp mobile temp
+thermo_modify temp ptemp
+WARNING: Temperature for thermo pressure is not for group all (src/thermo.cpp:520)
+
+thermo 1000
+
+# dump 1 all atom 1000 race.lammpstrj
+run 10000
+Generated 0 of 1 mixed pair_coeff terms from geometric mixing rule
+Per MPI rank memory allocation (min/avg/max) = 3.678 | 3.678 | 3.678 Mbytes
+   Step          Temp          E_pair         E_mol          TotEng         Press          Volume    
+         0   150           -2063.91        0             -454.27256      2.8467694      9350000      
+      1000   197.01825     -1960.8576      0              153.32871      13.788868      9350000      
+      2000   216.32291     -2037.8959      0              283.44712      13.915645      9350000      
+      3000   239.06547     -2059.8437      0              505.54782      14.550975      9350000      
+      4000   266.60476     -2076.3375      0              784.57583      17.457495      9350000      
+      5000   299.6816      -2109.6562      0              1106.202       19.375766      9350000      
+      6000   335.17037     -2129.5487      0              1467.1364      24.636144      9350000      
+      7000   367.9265      -2101.7855      0              1846.4029      33.591291      9350000      
+      8000   404.7304      -2144.8541      0              2198.2739      43.134333      9350000      
+      9000   435.75368     -2180.3183      0              2495.7179      53.466409      9350000      
+     10000   457.96804     -2194.6681      0              2719.7486      64.521177      9350000      
+Loop time of 12.1468 on 1 procs for 10000 steps with 13073 atoms
+
+Performance: 71.130 ns/day, 0.337 hours/ns, 823.265 timesteps/s
+99.8% CPU use with 1 MPI tasks x 1 OpenMP threads
+
+MPI task timing breakdown:
+Section |  min time  |  avg time  |  max time  |%varavg| %total
+---------------------------------------------------------------
+Pair    | 9.4112     | 9.4112     | 9.4112     |   0.0 | 77.48
+Neigh   | 1.1197     | 1.1197     | 1.1197     |   0.0 |  9.22
+Comm    | 0.042153   | 0.042153   | 0.042153   |   0.0 |  0.35
+Output  | 0.00041578 | 0.00041578 | 0.00041578 |   0.0 |  0.00
+Modify  | 1.2618     | 1.2618     | 1.2618     |   0.0 | 10.39
+Other   |            | 0.3115     |            |       |  2.56
+
+Nlocal:          13073 ave       13073 max       13073 min
+Histogram: 1 0 0 0 0 0 0 0 0 0
+Nghost:            590 ave         590 max         590 min
+Histogram: 1 0 0 0 0 0 0 0 0 0
+Neighs:         112709 ave      112709 max      112709 min
+Histogram: 1 0 0 0 0 0 0 0 0 0
+
+Total # of neighbors = 112709
+Ave neighs/atom = 8.62151
+Neighbor list builds = 129
+Dangerous builds = 0
+
+Total wall time: 0:00:12
diff --git a/examples/mesh/log.8May22.marble_race.g++.4 b/examples/mesh/log.8May22.marble_race.g++.4
new file mode 100644
index 0000000000..a728e292cf
--- /dev/null
+++ b/examples/mesh/log.8May22.marble_race.g++.4
@@ -0,0 +1,163 @@
+LAMMPS (4 May 2022)
+OMP_NUM_THREADS environment is not set. Defaulting to 1 thread. (src/comm.cpp:98)
+  using 1 OpenMP thread(s) per MPI task
+
+units real
+lattice         sc 5.0
+Lattice spacing in x,y,z = 5 5 5
+region          box block -110 60 -30 220 -90 130  units box
+create_box      2 box
+Created orthogonal box = (-110 -30 -90) to (60 220 130)
+  1 by 2 by 2 MPI processor grid
+region          particles cylinder y 0 -30 47 130 190 units box
+create_atoms    1 region particles
+Created 3601 atoms
+  using lattice units in orthogonal box = (-110 -30 -90) to (60 220 130)
+  create_atoms CPU = 0.001 seconds
+region          lid cylinder y 0 -30 47 190 200 units box
+group           mobile type 1
+3601 atoms in group mobile
+
+create_atoms    2 mesh race_track.stl units box
+Reading STL object Georgs Murmelbahn from file race_track.stl
+  read 9472 triangles with 1.00 atoms per triangle
+Created 9472 atoms
+  using box units in orthogonal box = (-110 -30 -90) to (60 220 130)
+  create_atoms CPU = 0.069 seconds
+group           mesh type 2
+9472 atoms in group mesh
+
+mass * 39.95
+pair_style lj/cut 8.76
+pair_coeff      1 1 0.2339 3.504
+pair_coeff      1 2 0.2339 7.008 $(7.008*2^(1.0/6.0))
+pair_coeff      1 2 0.2339 7.008 7.8662140345520858986
+pair_coeff      2 2 0.0    1.0
+
+balance 1.1 shift xyz 10 1.01
+Balancing ...
+Generated 0 of 1 mixed pair_coeff terms from geometric mixing rule
+Neighbor list info ...
+  update every 1 steps, delay 10 steps, check yes
+  max neighbors/atom: 2000, page size: 100000
+  master list distance cutoff = 10.76
+  ghost atom cutoff = 10.76
+  binsize = 5.38, bins = 32 47 41
+  1 neighbor lists, perpetual/occasional/extra = 1 0 0
+  (1) pair lj/cut, perpetual
+      attributes: half, newton on
+      pair build: half/bin/atomonly/newton
+      stencil: half/bin/3d
+      bin: standard
+ rebalancing time: 0.001 seconds
+  iteration count = 17
+  initial/final maximal load/proc = 6175 4316
+  initial/final imbalance factor  = 1.8893903 1.3205844
+  x cuts: 0 1
+  y cuts: 0 0.5859375 1
+  z cuts: 0 0.36376953 1
+
+neigh_modify exclude type 2 2
+timestep 1.0
+
+fix track mesh setforce 0.0 0.0 0.0
+fix pull mobile addforce 0.0 -0.05 0.0 region particles
+fix dir  mobile oneway 10 lid -y
+fix move mobile nve
+fix load all balance 1000 1.1 shift xyz 10 1.01 weight neigh 0.5 weight group 2 mesh 0.1 mobile 1.0
+
+minimize 0.0 0.0 1000 1000
+WARNING: Using 'neigh_modify every 1 delay 0 check yes' setting during minimization (src/min.cpp:187)
+Generated 0 of 1 mixed pair_coeff terms from geometric mixing rule
+Per MPI rank memory allocation (min/avg/max) = 4.506 | 4.54 | 4.57 Mbytes
+   Step          Temp          E_pair         E_mol          TotEng         Press          Volume    
+         0   0              737062.81      0              737062.81      21986.781      9350000      
+        67   0             -2063.91        0             -2063.91       -5.0227698      9350000      
+Loop time of 0.375696 on 4 procs for 67 steps with 13073 atoms
+
+99.0% CPU use with 4 MPI tasks x 1 OpenMP threads
+
+Minimization stats:
+  Stopping criterion = linesearch alpha is zero
+  Energy initial, next-to-last, final = 
+      737062.806250145  -2063.90998808134  -2063.90998808134
+  Force two-norm initial, final = 689296.27 22.226599
+  Force max component initial, final = 336546.89 0.90593277
+  Final line search alpha, max atom move = 1.2850327e-11 1.1641532e-11
+  Iterations, force evaluations = 67 393
+
+MPI task timing breakdown:
+Section |  min time  |  avg time  |  max time  |%varavg| %total
+---------------------------------------------------------------
+Pair    | 0.0026793  | 0.10176    | 0.30773    |  39.0 | 27.09
+Neigh   | 0.0083945  | 0.011206   | 0.01705    |   3.3 |  2.98
+Comm    | 0.0076602  | 0.21941    | 0.32088    |  27.2 | 58.40
+Output  | 0          | 0          | 0          |   0.0 |  0.00
+Modify  | 0.0035482  | 0.0090038  | 0.018889   |   6.2 |  2.40
+Other   |            | 0.03432    |            |       |  9.13
+
+Nlocal:        3268.25 ave        4314 max        2222 min
+Histogram: 2 0 0 0 0 0 0 0 0 2
+Nghost:        1010.75 ave        1101 max         792 min
+Histogram: 1 0 0 0 0 0 0 0 1 2
+Neighs:        16694.5 ave       52643 max           0 min
+Histogram: 2 0 1 0 0 0 0 0 0 1
+
+Total # of neighbors = 66778
+Ave neighs/atom = 5.1080854
+Neighbor list builds = 5
+Dangerous builds = 0
+
+reset_timestep 0 time 0.0
+
+velocity mobile create 150.0 54634234
+compute ptemp mobile temp
+thermo_modify temp ptemp
+WARNING: Temperature for thermo pressure is not for group all (src/thermo.cpp:520)
+
+thermo 1000
+
+# dump 1 all atom 1000 race.lammpstrj
+run 10000
+Generated 0 of 1 mixed pair_coeff terms from geometric mixing rule
+Per MPI rank memory allocation (min/avg/max) = 3.481 | 3.598 | 3.711 Mbytes
+   Step          Temp          E_pair         E_mol          TotEng         Press          Volume    
+         0   150           -2063.91        0             -454.27256      2.8467694      9350000      
+      1000   197.01825     -1960.8576      0              153.32871      13.788868      9350000      
+      2000   216.32291     -2037.8959      0              283.44712      13.915645      9350000      
+      3000   239.06547     -2059.8437      0              505.54782      14.550975      9350000      
+      4000   266.60476     -2076.3375      0              784.57583      17.457495      9350000      
+      5000   299.6816      -2109.6562      0              1106.202       19.375766      9350000      
+      6000   335.17037     -2129.5487      0              1467.1364      24.636144      9350000      
+      7000   367.9265      -2101.7855      0              1846.4029      33.591291      9350000      
+      8000   404.7304      -2144.8541      0              2198.2739      43.134332      9350000      
+      9000   435.7537      -2180.3187      0              2495.7178      53.466393      9350000      
+     10000   457.96586     -2194.6411      0              2719.7522      64.522003      9350000      
+Loop time of 4.62401 on 4 procs for 10000 steps with 13073 atoms
+
+Performance: 186.851 ns/day, 0.128 hours/ns, 2162.627 timesteps/s
+99.3% CPU use with 4 MPI tasks x 1 OpenMP threads
+
+MPI task timing breakdown:
+Section |  min time  |  avg time  |  max time  |%varavg| %total
+---------------------------------------------------------------
+Pair    | 2.0661     | 2.4819     | 2.9751     |  25.5 | 53.67
+Neigh   | 0.18188    | 0.2993     | 0.52875    |  25.5 |  6.47
+Comm    | 0.40938    | 1.17       | 1.736      |  51.7 | 25.30
+Output  | 0.00028188 | 0.00049757 | 0.0006548  |   0.0 |  0.01
+Modify  | 0.22514    | 0.32517    | 0.53165    |  21.2 |  7.03
+Other   |            | 0.3472     |            |       |  7.51
+
+Nlocal:        3268.25 ave        6890 max        1643 min
+Histogram: 1 2 0 0 0 0 0 0 0 1
+Nghost:           1701 ave        2074 max        1456 min
+Histogram: 2 0 0 0 0 1 0 0 0 1
+Neighs:        28177.5 ave       34170 max       21435 min
+Histogram: 1 1 0 0 0 0 0 0 0 2
+
+Total # of neighbors = 112710
+Ave neighs/atom = 8.6215865
+Neighbor list builds = 129
+Dangerous builds = 0
+
+Total wall time: 0:00:05
diff --git a/examples/mesh/log.8May22.mesh_box.g++.1 b/examples/mesh/log.8May22.mesh_box.g++.1
new file mode 100644
index 0000000000..ba8ab1935d
--- /dev/null
+++ b/examples/mesh/log.8May22.mesh_box.g++.1
@@ -0,0 +1,146 @@
+LAMMPS (4 May 2022)
+OMP_NUM_THREADS environment is not set. Defaulting to 1 thread. (src/comm.cpp:98)
+  using 1 OpenMP thread(s) per MPI task
+
+units real
+atom_style hybrid sphere bond
+WARNING: Atom style hybrid defines both, per-type and per-atom masses; both must be set, but only per-atom masses will be used (src/atom_vec_hybrid.cpp:130)
+
+lattice         sc 5.0
+Lattice spacing in x,y,z = 5 5 5
+region          box block 50 250 50 250 50 250 units box
+create_box      2 box
+Created orthogonal box = (50 50 50) to (250 250 250)
+  1 by 1 by 1 MPI processor grid
+region particles block 110 190 110 190 110 190 units box
+create_atoms    1 region particles
+Created 4913 atoms
+  using lattice units in orthogonal box = (50 50 50) to (250 250 250)
+  create_atoms CPU = 0.001 seconds
+region          lid block 100 110 50 250 50 250 units box
+group           mobile type 1
+4913 atoms in group mobile
+set type 1 diameter 7.0
+Setting atom values ...
+  4913 settings made for diameter
+
+create_atoms    2 mesh open_box.stl units box radthresh 4.0
+Reading STL object Open Box from file open_box.stl
+  read 10 triangles with 512.00 atoms per triangle
+Created 5120 atoms
+  using box units in orthogonal box = (50 50 50) to (250 250 250)
+  create_atoms CPU = 0.001 seconds
+group           mesh type 2
+5120 atoms in group mesh
+
+pair_style lj/cut 8.76
+pair_coeff      1 1 0.2339 3.504
+pair_coeff      1 2 0.2339 7.008 $(7.008*2^(1.0/6.0))
+pair_coeff      1 2 0.2339 7.008 7.8662140345520858986
+pair_coeff      2 2 0.0    1.0
+mass *  39.95
+
+neigh_modify exclude type 2 2
+
+timestep 1.0
+
+run 0 post no
+WARNING: No fixes with time integration, atoms won't move (src/verlet.cpp:60)
+Generated 0 of 1 mixed pair_coeff terms from geometric mixing rule
+Neighbor list info ...
+  update every 1 steps, delay 10 steps, check yes
+  max neighbors/atom: 2000, page size: 100000
+  master list distance cutoff = 10.76
+  ghost atom cutoff = 10.76
+  binsize = 5.38, bins = 38 38 38
+  1 neighbor lists, perpetual/occasional/extra = 1 0 0
+  (1) pair lj/cut, perpetual
+      attributes: half, newton on
+      pair build: half/bin/newton
+      stencil: half/bin/3d
+      bin: standard
+Per MPI rank memory allocation (min/avg/max) = 5.641 | 5.641 | 5.641 Mbytes
+   Step          Temp          E_pair         E_mol          TotEng         Press     
+         0   0             -1778.6527      0             -1778.6527     -27.271044    
+Loop time of 1.544e-06 on 1 procs for 0 steps with 10033 atoms
+
+
+fix dir  mobile oneway 10 lid x
+fix move mobile nve
+fix load all balance 1000 1.1 shift xyz 10 1.01 weight neigh 0.5 weight group 2 mesh 0.1 mobile 1.0
+fix rot  mesh move rotate 150.0 150.0 150.0 1.0 0.0 0.0 500000.0 units box
+
+reset_timestep 0 time 0.0
+velocity mobile create 150.0 54634234
+
+compute ptemp mobile temp
+thermo_modify temp ptemp
+WARNING: Temperature for thermo pressure is not for group all (src/thermo.cpp:520)
+
+thermo 200
+compute ke all ke/atom
+#dump 2 all movie 200 mesh.mkv c_ke radius size 960 1440 zoom 1.5 box no 0.0 view 120 180
+#dump_modify 2 bitrate 4000 framerate 12 color orange 1.0 0.5 0.0 amap min max cf 0.0 6 min blue 0.1 fuchsia 0.2 red 0.4 orange 0.6 yellow max white
+#dump 1 all custom 500 open_box.lammpstrj id type mol x y z vx vy vz
+
+
+run 5000
+Generated 0 of 1 mixed pair_coeff terms from geometric mixing rule
+Per MPI rank memory allocation (min/avg/max) = 6.157 | 6.157 | 6.157 Mbytes
+   Step          Temp          E_pair         E_mol          TotEng         Press          Volume    
+         0   150           -1778.6527      0              417.60814     -14.721534      8000000      
+       200   179.71663     -2353.7693      0              277.59468      14.591893      8000000      
+       400   181.98778     -2346.3965      0              318.22101      14.694157      8000000      
+       600   183.48983     -2373.9255      0              312.68467      15.032091      8000000      
+       800   183.99836     -2366.1103      0              327.94574      15.154418      8000000      
+      1000   186.19063     -2405.1258      0              321.02886      14.881141      8000000      
+      1200   185.74206     -2396.0437      0              323.54307      15.729931      8000000      
+      1400   188.03786     -2424.6865      0              328.5147       15.077756      8000000      
+      1600   187.4202      -2437.7998      0              306.35792      16.323387      8000000      
+      1800   186.98057     -2417.0238      0              320.69695      14.445045      8000000      
+      2000   185.81788     -2404.4322      0              316.26462      15.382214      8000000      
+      2200   185.60774     -2386.4973      0              331.12284      15.033058      8000000      
+      2400   187.14304     -2409.5005      0              330.5991       15.226519      8000000      
+      2600   186.63304     -2410.3607      0              322.27162      15.010459      8000000      
+      2800   185.96347     -2421.8039      0              301.02465      15.577503      8000000      
+      3000   185.22976     -2387.0491      0              325.03669      16.111731      8000000      
+      3200   186.32095     -2414.4926      0              313.5702       14.41138       8000000      
+      3400   186.75847     -2414.3018      0              320.16692      15.224392      8000000      
+      3600   184.82282     -2391.1485      0              314.97889      15.96507       8000000      
+      3800   185.71262     -2381.6131      0              337.54256      14.870266      8000000      
+      4000   186.58931     -2399.6637      0              332.32827      14.434356      8000000      
+      4200   185.27225     -2400.3577      0              312.35024      15.429373      8000000      
+      4400   183.90694     -2367.1261      0              325.59133      16.558529      8000000      
+      4600   185.84295     -2403.8288      0              317.23517      14.371372      8000000      
+      4800   184.58178     -2359.4293      0              343.16897      15.516075      8000000      
+      5000   186.15467     -2412.7678      0              312.86031      15.366761      8000000      
+Loop time of 9.53365 on 1 procs for 5000 steps with 10033 atoms
+
+Performance: 45.313 ns/day, 0.530 hours/ns, 524.458 timesteps/s
+99.8% CPU use with 1 MPI tasks x 1 OpenMP threads
+
+MPI task timing breakdown:
+Section |  min time  |  avg time  |  max time  |%varavg| %total
+---------------------------------------------------------------
+Pair    | 5.3496     | 5.3496     | 5.3496     |   0.0 | 56.11
+Bond    | 0.00036644 | 0.00036644 | 0.00036644 |   0.0 |  0.00
+Neigh   | 3.1301     | 3.1301     | 3.1301     |   0.0 | 32.83
+Comm    | 0.03837    | 0.03837    | 0.03837    |   0.0 |  0.40
+Output  | 0.0010388  | 0.0010388  | 0.0010388  |   0.0 |  0.01
+Modify  | 0.88529    | 0.88529    | 0.88529    |   0.0 |  9.29
+Other   |            | 0.1289     |            |       |  1.35
+
+Nlocal:          10033 ave       10033 max       10033 min
+Histogram: 1 0 0 0 0 0 0 0 0 0
+Nghost:              0 ave           0 max           0 min
+Histogram: 1 0 0 0 0 0 0 0 0 0
+Neighs:          98536 ave       98536 max       98536 min
+Histogram: 1 0 0 0 0 0 0 0 0 0
+
+Total # of neighbors = 98536
+Ave neighs/atom = 9.8211901
+Ave special neighs/atom = 0
+Neighbor list builds = 374
+Dangerous builds = 0
+
+Total wall time: 0:00:09
diff --git a/examples/mesh/log.8May22.mesh_box.g++.4 b/examples/mesh/log.8May22.mesh_box.g++.4
new file mode 100644
index 0000000000..a255fdd55c
--- /dev/null
+++ b/examples/mesh/log.8May22.mesh_box.g++.4
@@ -0,0 +1,146 @@
+LAMMPS (4 May 2022)
+OMP_NUM_THREADS environment is not set. Defaulting to 1 thread. (src/comm.cpp:98)
+  using 1 OpenMP thread(s) per MPI task
+
+units real
+atom_style hybrid sphere bond
+WARNING: Atom style hybrid defines both, per-type and per-atom masses; both must be set, but only per-atom masses will be used (src/atom_vec_hybrid.cpp:130)
+
+lattice         sc 5.0
+Lattice spacing in x,y,z = 5 5 5
+region          box block 50 250 50 250 50 250 units box
+create_box      2 box
+Created orthogonal box = (50 50 50) to (250 250 250)
+  1 by 2 by 2 MPI processor grid
+region particles block 110 190 110 190 110 190 units box
+create_atoms    1 region particles
+Created 4913 atoms
+  using lattice units in orthogonal box = (50 50 50) to (250 250 250)
+  create_atoms CPU = 0.001 seconds
+region          lid block 100 110 50 250 50 250 units box
+group           mobile type 1
+4913 atoms in group mobile
+set type 1 diameter 7.0
+Setting atom values ...
+  4913 settings made for diameter
+
+create_atoms    2 mesh open_box.stl units box radthresh 4.0
+Reading STL object Open Box from file open_box.stl
+  read 10 triangles with 512.00 atoms per triangle
+Created 5120 atoms
+  using box units in orthogonal box = (50 50 50) to (250 250 250)
+  create_atoms CPU = 0.000 seconds
+group           mesh type 2
+5120 atoms in group mesh
+
+pair_style lj/cut 8.76
+pair_coeff      1 1 0.2339 3.504
+pair_coeff      1 2 0.2339 7.008 $(7.008*2^(1.0/6.0))
+pair_coeff      1 2 0.2339 7.008 7.8662140345520858986
+pair_coeff      2 2 0.0    1.0
+mass *  39.95
+
+neigh_modify exclude type 2 2
+
+timestep 1.0
+
+run 0 post no
+WARNING: No fixes with time integration, atoms won't move (src/verlet.cpp:60)
+Generated 0 of 1 mixed pair_coeff terms from geometric mixing rule
+Neighbor list info ...
+  update every 1 steps, delay 10 steps, check yes
+  max neighbors/atom: 2000, page size: 100000
+  master list distance cutoff = 10.76
+  ghost atom cutoff = 10.76
+  binsize = 5.38, bins = 38 38 38
+  1 neighbor lists, perpetual/occasional/extra = 1 0 0
+  (1) pair lj/cut, perpetual
+      attributes: half, newton on
+      pair build: half/bin/newton
+      stencil: half/bin/3d
+      bin: standard
+Per MPI rank memory allocation (min/avg/max) = 5.457 | 5.476 | 5.493 Mbytes
+   Step          Temp          E_pair         E_mol          TotEng         Press     
+         0   0             -1778.6527      0             -1778.6527     -27.271044    
+Loop time of 3.34625e-06 on 4 procs for 0 steps with 10033 atoms
+
+
+fix dir  mobile oneway 10 lid x
+fix move mobile nve
+fix load all balance 1000 1.1 shift xyz 10 1.01 weight neigh 0.5 weight group 2 mesh 0.1 mobile 1.0
+fix rot  mesh move rotate 150.0 150.0 150.0 1.0 0.0 0.0 500000.0 units box
+
+reset_timestep 0 time 0.0
+velocity mobile create 150.0 54634234
+
+compute ptemp mobile temp
+thermo_modify temp ptemp
+WARNING: Temperature for thermo pressure is not for group all (src/thermo.cpp:520)
+
+thermo 200
+compute ke all ke/atom
+#dump 2 all movie 200 mesh.mkv c_ke radius size 960 1440 zoom 1.5 box no 0.0 view 120 180
+#dump_modify 2 bitrate 4000 framerate 12 color orange 1.0 0.5 0.0 amap min max cf 0.0 6 min blue 0.1 fuchsia 0.2 red 0.4 orange 0.6 yellow max white
+#dump 1 all custom 500 open_box.lammpstrj id type mol x y z vx vy vz
+
+
+run 5000
+Generated 0 of 1 mixed pair_coeff terms from geometric mixing rule
+Per MPI rank memory allocation (min/avg/max) = 5.998 | 6.028 | 6.059 Mbytes
+   Step          Temp          E_pair         E_mol          TotEng         Press          Volume    
+         0   150           -1778.6527      0              417.60814     -14.721534      8000000      
+       200   178.85669     -2312.6876      0              306.08541      15.832907      8000000      
+       400   182.01107     -2334.4183      0              330.54024      14.457841      8000000      
+       600   183.04428     -2348.8534      0              331.23317      15.166984      8000000      
+       800   185.68203     -2396.0112      0              322.69659      14.155191      8000000      
+      1000   185.58372     -2376.7148      0              340.55362      14.125122      8000000      
+      1200   184.96212     -2401.9634      0              306.20372      14.791738      8000000      
+      1400   186.67578     -2392.327       0              340.93097      14.859731      8000000      
+      1600   186.98761     -2418.1463      0              319.67739      15.173567      8000000      
+      1800   188.23692     -2431.082       0              325.03382      15.02044       8000000      
+      2000   184.29434     -2378.8865      0              319.50307      16.504785      8000000      
+      2200   187.1493      -2395.48        0              344.71119      14.464979      8000000      
+      2400   186.60933     -2409.0279      0              323.25719      15.165523      8000000      
+      2600   185.29773     -2391.0001      0              322.08084      15.967212      8000000      
+      2800   185.05884     -2376.7398      0              332.84337      14.263097      8000000      
+      3000   185.03409     -2385.6071      0              323.61377      15.053931      8000000      
+      3200   185.37823     -2400.0415      0              314.21823      16.076039      8000000      
+      3400   187.03675     -2424.3537      0              314.18955      14.596162      8000000      
+      3600   185.55628     -2397.9284      0              318.93826      14.036493      8000000      
+      3800   187.10605     -2408.018       0              331.54002      14.527143      8000000      
+      4000   186.2403      -2410.4382      0              316.44375      15.342402      8000000      
+      4200   187.32974     -2389.1104      0              353.72276      13.931896      8000000      
+      4400   189.27742     -2435.0111      0              336.33944      14.316321      8000000      
+      4600   189.26108     -2453.7231      0              317.38821      15.615617      8000000      
+      4800   187.9187      -2429.9397      0              321.51688      15.226135      8000000      
+      5000   186.91414     -2428.2233      0              308.52478      15.780318      8000000      
+Loop time of 3.39374 on 4 procs for 5000 steps with 10033 atoms
+
+Performance: 127.293 ns/day, 0.189 hours/ns, 1473.303 timesteps/s
+99.3% CPU use with 4 MPI tasks x 1 OpenMP threads
+
+MPI task timing breakdown:
+Section |  min time  |  avg time  |  max time  |%varavg| %total
+---------------------------------------------------------------
+Pair    | 1.0734     | 1.4274     | 1.7321     |  24.7 | 42.06
+Bond    | 0.00031015 | 0.00032026 | 0.00033017 |   0.0 |  0.01
+Neigh   | 0.65409    | 0.83138    | 1.0016     |  17.4 | 24.50
+Comm    | 0.33062    | 0.82642    | 1.379      |  52.2 | 24.35
+Output  | 0.00070506 | 0.00086089 | 0.0010669  |   0.0 |  0.03
+Modify  | 0.21589    | 0.24797    | 0.28078    |   5.5 |  7.31
+Other   |            | 0.05935    |            |       |  1.75
+
+Nlocal:        2508.25 ave        2930 max        2111 min
+Histogram: 2 0 0 0 0 0 0 0 1 1
+Nghost:         1038.5 ave        1110 max         985 min
+Histogram: 2 0 0 0 0 0 1 0 0 1
+Neighs:        24668.2 ave       32893 max       17327 min
+Histogram: 1 1 0 0 0 0 0 1 0 1
+
+Total # of neighbors = 98673
+Ave neighs/atom = 9.834845
+Ave special neighs/atom = 0
+Neighbor list builds = 371
+Dangerous builds = 0
+
+Total wall time: 0:00:03
diff --git a/examples/mesh/open_box.stl b/examples/mesh/open_box.stl
new file mode 100644
index 0000000000..f9ef56d64a
--- /dev/null
+++ b/examples/mesh/open_box.stl
@@ -0,0 +1,72 @@
+solid Open Box
+    facet normal  0.0 0.0 1.0
+        outer loop
+            vertex 100.0 100.0 100.0
+            vertex 200.0 100.0 100.0
+            vertex 200.0 200.0 100.0
+        endloop
+    endfacet
+    facet normal  0.0 0.0 1.0
+        outer loop
+            vertex 100.0 100.0 100.0
+            vertex 100.0 200.0 100.0
+            vertex 200.0 200.0 100.0
+        endloop
+    endfacet
+    facet normal  0.0 0.0 -1.0
+        outer loop
+            vertex 100.0 100.0 200.0
+            vertex 200.0 100.0 200.0
+            vertex 200.0 200.0 200.0
+        endloop
+    endfacet
+    facet normal  0.0 0.0 -1.0
+        outer loop
+            vertex 100.0 100.0 200.0
+            vertex 100.0 200.0 200.0
+            vertex 200.0 200.0 200.0
+        endloop
+    endfacet
+    facet normal  0.0 1.0 0.0
+        outer loop
+            vertex 100.0 100.0 100.0
+            vertex 200.0 100.0 100.0
+            vertex 200.0 100.0 200.0
+        endloop
+    endfacet
+    facet normal  0.0 1.0 0.0
+        outer loop
+            vertex 100.0 100.0 100.0
+            vertex 100.0 100.0 200.0
+            vertex 200.0 100.0 200.0
+        endloop
+    endfacet
+    facet normal  0.0 -1.0 0.0
+        outer loop
+            vertex 100.0 200.0 100.0
+            vertex 200.0 200.0 100.0
+            vertex 200.0 200.0 200.0
+        endloop
+    endfacet
+    facet normal  0.0 -1.0 0.0
+        outer loop
+            vertex 100.0 200.0 100.0
+            vertex 100.0 200.0 200.0
+            vertex 200.0 200.0 200.0
+        endloop
+    endfacet
+    facet normal -1.0 0.0 0.0
+        outer loop
+            vertex 200.0 100.0 100.0
+            vertex 200.0 200.0 100.0
+            vertex 200.0 200.0 200.0
+        endloop
+    endfacet
+    facet normal  -1.0 0.0 0.0
+        outer loop
+            vertex 200.0 100.0 100.0
+            vertex 200.0 100.0 200.0
+            vertex 200.0 200.0 200.0
+        endloop
+    endfacet
+endsolid
diff --git a/examples/mesh/race_track.stl b/examples/mesh/race_track.stl
new file mode 120000
index 0000000000..e05acc5019
--- /dev/null
+++ b/examples/mesh/race_track.stl
@@ -0,0 +1 @@
+../PACKAGES/machdyn/funnel_flow/boundary.stl
\ No newline at end of file