From 6902d886f237739b2ef82c46b3ea1a4fe13a8d40 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Fri, 13 Dec 2024 15:06:12 -0500 Subject: [PATCH] update the "Reading and plotting manually" portion of docs this adds how to exclude ghost cells. --- docs/source/manual_plot.png | Bin 20382 -> 21175 bytes docs/source/output.rst | 35 ++++++++++++++++++++++++++++------- 2 files changed, 28 insertions(+), 7 deletions(-) diff --git a/docs/source/manual_plot.png b/docs/source/manual_plot.png index 61619fc2be8cb0102da576283cd6cc3accec7307..358368709f5fc5ad32f70be36a50d178aa99ec4c 100644 GIT binary patch literal 21175 zcma&ObySpV7(F^PNDd9s0}4uaBLYKrt27J}($Xa$If#^`pnx<2BGR3bA|(w1BHdm0 z9nbpZ{o~%XjtdX-&D-Cz_p_h<8l|P7NO*_#4g>-rR92GLfk4n<5D40PTx{^0R&;j` z@I%Z)!N^0`8R6k&;bsF-xA1UrboOwxvt;tLadWqGb`szgEs300#m|i~9abkVb4{BlHZ`}guZ!ypA6DzfJ1oD@KYD;Q`Savy=jxr42UskE zI+(?r+Sa-g71qq0X^sN!h(A|1sbcpN1Kzp2A%+^W;&6z2t?!+$bgT^c&1avkT|-%w z6CTOS%jc0L&SNmKu;h$&n}4;5Mag)neSi*E#WF9 zgiTLBIoh6~x~F@29+!s;AIg@AySnmkijW0e@`<|tA-Fi*UD#R6y#qHWd-7Vp#9(Qx zNEZ^!XVty9-^19`=yN#Tc$+Td;4&pRS)j<_dIC%z%Z3wLM?{`7;p)4v~CdptMuNwlJNfH}8bAGFTXF@|S3v zuP|`5w6q{KGqFe{L8(~_L#|?+(_u3Vj%>|vZjRbt?#a#1#~zlx}k~d zd*6F;-#aPJgqw;)MJg(Xdi$+DEXY}^$6PQ|_l1k;c#qgP_L_b`lwd*p8`-z4%F%!D z{hCaZl9NAzwsxCN8BBw&U-}<>zrT&m%7r6jKgK@m+Cy*D?C;ZX>IZ&Z={WVppr+iU zfpUF)eS{zdbx?l(_z^6z+ejuPBt$>q52r_d6?B;3B%l{3FK;|hkv}!7bz$bD#Gm)w zEc4tQH9)Y+$S{tNkJnRGP*PD@oor3v>5?aHi55&RHWtUZ7ZVPZT#(vnfGPN%?#$Wo zl6qywKVtg9>GC&nVYdgY3UAQO716G3&5T9O%+ooi>L)L=f63EE#Q&c=)1C>xxCY1S z9k)re&*2}7)?l>4va+t}Y6oo{oshu{5ex_=*ORco>x++tg}9KtFK|%3%Bo`-NEiL4|gcn@T>@nEaBIaWAl~|Gazk7#lY;1$lrj?dh z2Bk*OnO?At>}-2urYtNhKPf_?|L2blw`wQot4GhaUqto>H#S2|Z+q}E>&8BuY8A%t zQGeuq`89_UOJQd3hr}me;@6R+0}IU$dV&HysBdlO0-GWpn;_e8UluK-(8;~;iMX~g zJYu*5WlL2_xeKXi*}BtQR66b_I`*0xsg)w;UWCRNXiOCL&OS1-eY`%PuVrD*&tS4r zYtoYUvB{rU`uj`7WM<38>A_lP2QJ>0)0I> zW+6R2_=Hd)*Vi^`C7X@xV!pDI`&yG0prPNx4Vt;yM^+wLWFaJtTvsL}(0G6Kzf*(Q z5O3SzGz|t>c$;~c;2Pz%QPm+iWbwX01WN~{1dGP7mCGN`kNCz zdSevS$hUuPY!{556`$q=#R(2DX2hr|@=_Na@(MhbS6>8nXCaYF0n-4ytil1sO1rpCemly&#@=Y7A-B1}O?4u2>p9@ziR;fH|$+6&gOz(_xFE8l|G(gQW1?Qx__T6m=O3+_>)I<>?Hpoe_IiGD!>cqBV} z<5HU1y!vSGoM1R=u_ywel@Z*WX{9vGIh5(1qk=Z@6%k$rY}mg~*?2~9(anj=ZV}S;Q-p%7-Bu?4D3nd_i=cp+=WP=2xcP7RiLGqUudqNeUI`Fkj`5 z-RGlbR{u@~t@GCJ<}Q^dfAiE&_b%8Fz1&~o0wEP02lg~b5@_%f&AaqQQ%9UO;54aM zJ9^)951+rTz<&t+_}+%LZen2Yu~&x`@8fiPAMRLAIph2#g>d%qbM)Z9#7>9K?rGs^ z(EmObudF%Oa>aV)3eKxc?Nt1<~qk_Z8Uvh0q zkL)+(9(R-f2ga5+SLY~%|6ITEK=S%%w70jnb!V<&rtvn0r#TV;C=Id96)!xN_VsFi ziq*l&cklDnlyO=aHqlSP)i#JuhGH5EL>T*Oog7x_kbw|Pok!czAK~~8cW3k-fN`&d z0IK8TV|I4-=ZDe_p4*QtEG(*K{tUg3A-pT>aaJwGKrG1+cTm~CiVlg)!PO?WqIa8J z*)3b6PEZdO5j8ZBwHJ1_`)Hf>thG-}pDMd#V>#1aj$DCFAi+sUD8o!PfeAO56|I_( zAuny_R(y;#KC2}D4fBTiO6w~O&};C_HXIJeo+vkq8_89y@mYwDzO%BjBG|xn|2_`D z9{RwzoE#R;r{7A)x1WA0pJNzEj`1zw<*$?-$1zo<{9a_#TjvCfK1_NV`W-(xx zz5}WC?86+${#7%zmVNbmcTYHyJ=Ikze9ii=VBboD(y+ddO|F184oo!{7$yoO3Isy2Q}sJ3Yx0yNb2_h47PVT5NNA7Y_U_vYwjFZvNgv4f6*8IR5kVQ z?*Qkq+McdH8!r!hZFX~zf_jOB9e)K^TJ;p)CTcy_nQp5|ronq@@!|-tOj=RCNNxxz zN?)dXAz3FPBu{~C-0!|iAK09impA?bpNhM){WX3rBMqCBazZ;ltjrKopiWIL=|DaR zV##zCo;NsvJ6xcQCz%m0xO|Mmd5?R5OLJQx=G~t4Gk`7v zPR^sW&QJRV2aqK3{zi!UL4ugNV^tkpL>NFUkhciJct%py(Mw-EvtBU=zmqRg-D%I9 zf*W8$wn}M~@lb8hQq8pP|PdHu~oy^WAU2W}GaZ$`R)(5w)f65E_4m~z4a zJg>@zHf;L`Y#Fe}^!2SVQ^9^D=vzf-sZF5)tMpxd3J3Q%HHurQ*%q9HY_`dBa9!;+ z?0tY`b~%QO(@Cl&zTW8uDXk#P?70zOJ1^oFDpIlVDBW@T>nL{$|1kv`F{~L2L?^(* z96dH$`I#aU{a>phT%}wBoM0i+p7WzfmWtFt##@yzn{iOfK^SOz0(}3iut1y|6w^995oOqu9AAeiRX`3<89UE zIZ5Nx76{f@t|h`aFn^YW8#HqOKc9rKVDw%F)m2!*ODIATw)@7eJP6W0ly5Nz&TqT~ zkC#DfU&jFw^`*ywI95XV^rytAZ`yFA#p61_f1$V5g7bNgWxYAN!k5(4I#Ne-1w|JB zO*l$v=o1?1_Q?oXJW6Fh9TA@LyoJqtJooaRjPD=I1z(n9R!lY{+PVWD1uJp+;3AX+%7~-_LMyY}tgo)O3lfqXL<(g_{$RtzVV7b02b= zUd7T1-SL`r#d*hLHvcCe%l|<6KL&z=Gr&MNJ3eurtU0Uw(tXkX0^_>2$Qc|Os%HLW z>$2e6-3K!Fjx5WeIePoM4TF2i=vk>-a+-dZrzQ57S=`mF7Vq}VMpNmP?-3Fb)*cQ? z0Lp<;Fgmq!ZP|z9H2#9g7c`8GW^mf2;(oLA9V1Ey{nnx6rss)adTB?BKdPA@qe%oELZ~T(wfYz7>cDO zB#;s?NQ!Q}I=>F26LoFGfrK$3X-?j833Ng;GKrCvLfH&Ft9)RFwqaIazEny5fCyF1`?j_O)dL0OJl zZ4Q%kFTPo%`@te4=BV`&qIUK)@X`{Ieg-q#xQaZ~XYC7EP1M4%A6>lqb@LKqHfntULK|HM%D-NIH8-n-&7#J~vRZ6}Dg zO9na>@tmWW0#2wXP9Ljn4TrQ^Sg3l#P^nU6K1Kd>-KI5?@B^jY16?WQY13jt_v-r^@RF#5cKMM4z z(QG9DMIHk+*u3W{b}Ly52O0GIlpqL;$I&6>y;GB~x+6hCSOnX!p+&$6P~-f}AwUY! zQd9<+D<~DS*ozxhjR!XR-#AB(sAY{t1|!AQ@&0V4Rl~v{`TD0B_yG z)$WV@h6ro0L;jofU7Nv3W89rav`v6|@KPa|Wc6e!fS`@Wo5e1cp97Ihf+lZM!EzD0 zC3Yyyw25HB29W<2)8uQvpt+YqPSWxpe!|xP=!;`K9eCPl8TTtBbH)}=6l&D0YWS)7 z_)*gVm|7^Hj#V1JQ@`v)i`v5CL>|f!Od!0Rd%Iu*hyp~7oenIS$+I*l4=Q*+3S5Tz zNO)uaV~TqB%{aFz$B;lDvi8oAfs8OpG9{`HM`;mn!gzY%{`Tc!-m!a;V09zG zN^$3>$7}Jn=iAD6-)~tb31b%*5)iO_j@r)=Pjf5*aC*AviwRY$(wT7|%g0%ul~e-I z9`~^{`VX!nJ?|5?-!YncHU4iwEb-whb0|>+`I)cpz2NhAG)qG5k*KBm&0!8wc$VeF z1&`kA*DP)Mw%vf->J^Yzz^ke!yVF6CTzX@b7A^$q1}9tZeu4J3whi1{m{7Yo@)JUS z=qomBYtL5X5@1N7^yXOSuIpoAr?xowLS8U?N7{vmTRmA;H8fYZoRZvfh?=6Qe z^*28&<^#%0UVXkaEWTxPX~HyT45~U92^bW|SVEXWX?F+cjcP&G2Q_H!cXNUe^BT z9%vu;*V(`i)PpPLtm)FQr8KB=g5ce|&~=`Z^!OkAgaU=(|N0%CFm;hi?cG}eYE~&K z2cXu_Yg5wJhJJ3ou(w7aSO5y+Q*yR~gjgbAwqdt9CLX5PzT{#)TkHC(*40YVDj#$5 zYqdzUIV-0NU>EU~-(PQl1u)iY>%E`tAJP>=Wd#&xc5d!)_UKdrTcVB4&6$9Kq4yF% z0{2>+_uK6PU|MP3hR$hJ0Wk|WPNrQhz|7+ul^^-hYbngph96@Th2;8!ZXgf}&HW

iG}JxcIdpn|L@Vx}pU17)Pr+f@?I1F)$=f76n&>-^KVeHZS1Y+C=rxWf3QQVF18*Z}O?+u3r6Ixkkxr&uA1O--83N4kNljR!3}N(g z18qu7B(0&u(vdU*2y+1appKD~gDxf}Ht5v<{kyD!ffzo45b z&nE%kXV~aIqw4d)zG{nFi)%-uBN)lNG&75eGrv8`Vo4k1`u-`s&I;U*S zje4(2;IZ9TA+Jn?=dJ$-?DJv4RNry7MKA_NZqL-@2HoD!&x?qNP$efP4?2~UmdXQe z*brAavt2Rsv|zo=pk#+ELupw_9*sFQusTh>V#*OJGh|=QGRWXB9*y?uu!Qu#CCP5g z%t&)G{KKE%^U^4E5|^RSO7LPjnsH}f6l=y90TG3d5 zwT%rk7Z)L<;`s>6(dHyE*j7j|4udyDWYz^E%Wp?))~%n*sK$vN9KTRQvjF{=>jONgsA=&uD>W52(%TZpY1&GcxO0}83-OR{dZ~J6)1MX2RG7y7{ zYv0+e4B&!NP)Bq@aCOl-5E61xmq)Ng-#cvM5+6$@O@z7q{bIfAz4u!Maif##HYN4_n2 zIBIcj25|3YLuvGxn7Sry!_Rekv|ec(z#~7Mjfgs( z5X|#oC=`a;`h&tTqUQO{h`0y~C{1)GP1H*f!gc{zqcrs$)U+D=1(&Fbo}AD!^OFnR zF-m^72LKQCNpcV40)Ond0!{ ztgWP)yuaVx&8`HEv#*FyAre237xjbyc8KDq_vXWG`N3u*3yQqQ(Pc?2Z^Zw{xX%}h zy*CtE!4k^pEO(lQHBIfGd240Gqk8Y03DBSTnh_3|EBxcvkhok^OrT((L>b!gcXohG z^3w@DX@mKvg*Vv!H_g~^Tgclk)PA@Y{kX%hC17a;PKD5cMDkgR-LF^4z0Z@-Cq$pL zg<_S412IOAY?s-G-v(8WVd>DiW3)VotXNNDR=6&%47l{11(oZ+Gc`PzvaR8ScL7uY zX%FyJ|2`Io z6pGiPDh&27wh3E|>Sm>nwh%FXj`DfsPpa=2ae=AH@4qh8P8D@!kz^q5enTGv$PB3f zV&M5{^Z(5OHd5+ST57Val&oX378{M!A9uQs&1ZU6rD^XTFYCO3@xQ3_a!}L%hkCTx z>DWAw4;j_x$IR0{La?*^9@xcs(8!zrG=Ra-co4=*O`B+S3p8H$EVpfzx4cOcu>GBV zdv*VwX4dQC1!ZDt1cb_-UjvIlZ6w3!eCB7#W-`O#4PK`{p;W%FDdD-o_x^+9 z;R6gdNMowNp*V72P|qpO)3kxHCDJf%OW{ck!9eyxV7OE=+fR3W)2Y8|DO{sr8Nx4U z0YMqrS=;uM*z)u%nfag4FBr&CoB6xKI43LnoqKvR!~HKh3u5d%9_l<+bic96S#Uoa_U6u zX)cskm`(BW$a6w)|^*{;XVrIYBuk5TAdPvlOI=E4rfW# z0^vV7Zgu(%VBU8fCQ4;}e8fPZC21HKOs+0lZY5BC9UXZf8D$t4gRWew%4NGh7~mK{ zT*OzpZ5{D6leUr=V4nDD`DdOz-_xMP0

m%RM$f7Fx`#&*^wON~djmK7O%FykZ- zlpJ$&^Ky=8_%fDno}j{<;%J@V{v?m!F%3{&SwTHPFV_iG6#xNaBC6^6$gApHk?y2t zs5muNMuW#ST8qlAx9s)pv4|lj6$2%IZFw!pQ88j(%MSrU6#KJSY>4B-#kGiLfoG#0)xczzvrrgg?eZjpPLk)9)!pxtrS1|Xx=$hU**ZUMJX&TK#K}& zLJ}1xN;U2FM9r6VT5Ib^H}Xlf+++1<$^5YA-+5pEz$*4o^u`=o_Ku;-J ziWqO^&d8!cajw_XWs5dbfKnxWP5P}*v8h19fCSkxFF2hKfJ_n?)S$m=&YCpc)C%(E zyivWzV$&wvcCS&C5D@|<*Pg45U{&Ncw*f-`I=d8HuapYSYm}N4G~Lyu#rRe&LSel8 zp;;#&T}l35K`{cY1hrB6blhHO5crC?{&l^ZZ%zFWeo*0@B$N$j`<-1#SH=HWm59e} zy?>iGuwhLCjwG0^jBj0vvP_Jk5~2+H1?RTf%=7Z&CoywT$KtI3EccKfv_@Hl*HrIX z6wXfbJq>I=&nS|>#+>-Q7oU|q2?3nliE2}_(8GrbcY|k?$+77i9VE&LeOM{N{em4| z1l~Fcy%i--Tcaz;d4>y^RTyGTW$dM-ABEb` zRw=jE1{y~Z`4k2rqmA5z^B@raq(8lt^AiN)CfNBHwzUxzww~%Vp?#O+ZHy_zGDYr;udH^Gf zPTUtYa40Y>Uc$^WZFm#5fScXcin7E8?k6gW2uf!NAXiL$6X15xGb^Pf|t$*oHd%^C}rT_;EtOqplHJP)`X$?mDO z_9Z`r)!zaHa_#4;S?s>+EN^GW!Qj6qI}1z$zZ!kqIhQ|Sx;1wbrp(~PD7RK@5vaFJ z24}V5KbiW3t=>Wm?JEB=L0B#WiBbMhT45s0oJ4;fqEM>LweijS@XufeaEh`t-)@YT z@Okfb(-wEU#{*n@Ep0=;i%&c5md*O`O4L)Q-c~8~+~S=MDs1 zUilt3Qt64CBJe-fomaOkGM}q@PW$grR%R7{9z7O&%gOTIFq&>n(5;|2bWQ(wcj>Q3 z^Y8bQd_HZZTSLp^XOo7W?X0-KD3%(wcw$z?^e8tr5KkGhH@7eLC`1xSG@*amloS-O z+}zyK=iVhIVgTOL-3`zC^r_lc4X9=BxSmdIoNzt;#w2CZaNMKh@2cE-K}J#khI!mT+yIs4`t z>Tx2}ZX;e=>**4gs#}@Eo>N;(U1d8QK5qZwG|-AiXz?_{YccwH_n^q|tM5sVonKm# ztT*bch!(`poXn#wI*R3MQolN%B;MI;wjbE+48f+8+ygi7%ZUE6NZtzRRkt9k$Z7%_GwR*%f2aW$%* zlX_gnG2_R>4{GgtP9x=+s<}3CPgc=x|r7!-%P`)3OwzE^Wwg>p-{q1k_`5@}Xf9^Jg$w%#OtuF4` z=Ff%ftUnvtP+2NhwM<}B%=t`{CgzTah)Av;H+lJIAPxAmgN|}lAGw{VZZ^}7(fY-b zut}c}J`NoA7$h3b`}Vqp`y81?qx^CJ`8NVki_ZHzxqnU_{X)_0dhCf#3}s0X@Uv)I zwh8btU@N~R$ix%FyGl$nNibKdquL_dqWbs=3qF^8r>ZG_LQLCqqp_}HETd3P_k@_0NcPH13*MZE2)GWW0PF&a`AStt9#1WSm5q9d zKUO-GMj2X4GXm@hIh`+y=3_2?klMHGaGH`<0_cGvG+C3PA& z3l-~~4C$xV4()x0W9(5|0yz(X6rl3V5u;p&do&OC1T+U^h$z*-#2?vGaE}i333BG} z`e{CaoVY#nyj+b|`LSc8#$uU>Od91c{jkt@cBO3H#ipAsxSey_N@ z1T3$gdT@z*M`PcoL6;B1^)V|&RT+h2@O+Hh$ayt)*;F`fA72MyA6fF|P~kYhk;7yR zK;VFjl-^}lk<5_^J~W(#q=3jLfI-Qd21#u=3Q7qu-;I&M@MdEoNNcX3XWBy6B`ogSVP@wyxycdJ0R9&i!pReq=O42VIX)6^sHi(s?xMM5~c8 z-yQbO4w$bUCTLP3>lwU$vJ&m&E2FvL7yt@1cf1=PuGWABoM%178%GR8+A7N$>75^s zXqIGsx>t#|^DVNW*LWep+HTN?DfW>b%aqrm+9t%j0tqDGw4rSYBp7wGOV|S!)%q$} zTi|>4^{ir+NPMw)+${~`c(g9iZqCZw%dsT43V3K<${rgnL*hr7Z=srKjyTyy@CX|? zE5~UhH)Ud4a`~$iKH_8SOC0ee3EKM!>w9pU z68SG%)_o)b`Y_&U=$X@__dPh^Tg(CV>r!;~+VDk)er+vf>_*kKb;JZFxmxa$(H2Q~(Wx_cLS4 zbJD0Q+AbT~(1X3P`7i$R&Q7TPEP^xcG;;5JMzCa?Vsc?ePv6a=?I}}r(Xbk zd*eZaCT0~2%#D8};=Li63b3j53x=~6G|~avnAVE9^=VX$a&>tYds=`mC-oH-!78bG zjvdN)+W8hptOH2aB-iUNZ`UhYS5KBo+#Lxu$0$7Mi4njdjgce@nBfGHI&r4Qb1$mb zku(Xtp1!A6PcEQrl&QtS*qKZ}wL#co@2Sxk5M+-92{FtHN=@@@=SlJsj=QM!Gkv4|TV zi36N;7-`Wq(7RT^eP4fsQk(xmko;)ayh-VC!m#Ec>qs>(Fsg3XEWj-+n6tcp<895< z5TiKYY@NHdzW(?Q>WWQYe}A$JO0DslH5p=j%XAUi6%Xv?teGDJ_}?yn?byU;k4=qK z3P@n(NP{}9tE;(jyqf7vi(ktQVmMt0k^vyXfMk%O1vN9cudn#L|KXq~BMcYa&D}j& z22`W=-!%05blm6y6ZLN>CYr!ws%}kVk;i#yiQ~tGE_D_H8N4sGh?^qhUMO{=8>k9; zo%lBWYVhI#a#8UtfXNx0a971*vGo8 zyw46eO@d*;RdP$MY2*JuB&QnpyhYU6}~_g=ejD3k=2We-l#D z4^?C{M?Y+v;Z!XqVq3Q%A)TUTd%#_qxWFUgedITor108)3lZTvmGH zg&74bJ8>J&|71|h1WMZ5+wW~v4O?Cut&4i^F`+I}{l7K05dvc&LrIAFE9JL)r)5uH z(uBrSP*L?<3(bDK&Ba)=-G);D<6ASJiLC9dk=E&=s@2JSh$M==Jc`+ zh4eG=Zx{mtu2_YgBtLH&(bCeg)Mx)m*#z!yP*`GT)XbKq&xOafKP+zy0d{Io;{2qV zl2cKkg{nUOFZ=J9Z|@87fJ(GqH5dkDynZK+MrEZ;hN0}ypyX=sLgUu|k1kz<7P0${ z{6O}YVKU!6SoK=6XaPiPcYkR_iVcEZyZdgvnC@+lkN!gZ@AAs+D~^?u{+_Gz?efPS zSpS^ASg9nxpi$JH=+5TNC`b%c#$mKZxi^;?S)Sa!g|(!8PBUC>#93NWXu1E$Xre zIXF1TIiILc#DDe6lC^S9m@k#^{K!bC*nSQ~!(J_Eb}1al*-mF*F*ZvcH?UojQdGdN ze?BoXikvj{BMiD)W^_HPruuvnZvkRDS{^5Dx| zTK|=z+jzvEHGNkXZ8L#Q71S>igk=nMbbf9vQ#>><2Lf(oEcYvKfM&Ec$l(rMMuzqj zTL+h`$Migz@(<2F97(p8D(Sd@O4_o0A9649LVi-B24nKTxM<12f&Y8yC^Us}El37wDxc7~v9mQ)H3>Yw;^#@~Iu z89j)fY4eMH)|y^yB#`h75jP}TKZyLIO86`2Pc*7+`yddlbzIrgpx~}bM|%89f>V|G z&|h=yn4cnL?tryp8 zE;}l@s0LM9c<{cF-OjHB_Nv(Gez6AH4~A8C{E7fjQpgF`vCJN|mt_C-lg7~I7x!+) zo_V9!#tM;XX09gAR=an?x+Tz_GJ(s{hMPfjFPTjQc;8DeuL!uue8rn&phCtJebKZ0 z%{y9xwF^e%EOoHfMeOEMHk;G~%Sa^J-`B4Bb1R2M@%=Btea5g&m4GAtuNm z;Y8zj>)*K5p~OnTJb>3|B@IN1iQQOnD3M@t*72p%pSLt<NW;M3l zazYkuKbmJHL~wCZE5EA;D+{$mb2uQyM2nH58k@XzMNj6=`lC7)zu+W(XM|@dy-kV_ ztPU8kHh9K>Ikn3c(P=)*;s8&Wx^neNO+(WNk@@zXl7rg%L7Ifmf?GCrNPC-E%-zs}6 zo`SBsWWS?yL8C^MY^c+=I)bQ#rsaHI6w-Qoy=HfH{XBVIj^(FG=ww8@mkM~XQbmNR zs^r$hwGRLl0VF~zm4@%lx}4uK&Lk1IOmb_x9Y|MFJPP)yKAsitc7KOTR7X2zo<>;m ze2&kTI7;uwS2dI=oIL<8Y8AeB=GDEpAuFY|xxb_sG~I!Cp7h0Z+LPcev0ke^WyJ-K z?%QRLzf!C^uh^3Qt^3N~pWQS9D8i>NlNtiVA(xzb2G#j#oNO7fK)a~1FY6_DB+agb zAm-;UtoUAPX@YfpP{{`zaZXw4R-#P;Ti!26$n%T`&w510y^pRj_U?ZdBa+{r4*1P- zXMqC?GHb$1ax4BK7|fc&HLW|?|0ei=rOR`#n?K)~(qL8^ibxDTJ;~|{bxwI~9Uexo z!w4{!_;!zmwlVOTNTW#L@%s_fG!ZC0w6^0$AdQAEE2HPn?^TI=ZI-B@nAww1eVLH2 zO``qv8+}n{Q_w6%4Rz}PG+=K!r>G*`$3lSigof|Ko*)&l-rLRBN;cL%zXK?wFk@+F zlf%iPTa)3g1+-|RsmXSKS!Vb3BvtFyh?w{4XO2GCLu<|nEy{)Wtj?upKWWj#25qWR z&42U9TFeRZ<2cU;2>mQYY>67g?lykKFw1JnIjxBN{b&DtzyBC#k8w8;4V-t0<-@(p zpbUbB*RLYAg@mL)<&}43ATrbr+HpCtH}8FGJIYsRC_ZFj*LmAuPc?s0N|4amZlI#E z_F3@vWr>%?woE_BWaO$ow+G#BvT_enXU={|OoL*?!&D_gwVi zx4$&p^7^i#^`?~I28#-sLxMHjn^Xt0RxL4I?H-y;odr80wy*8yX-g+yoXbVT#1?aF zHs0PlN&~dXgtChgYASTe+#;s47liXSR|HpRcyo>jNw-6#tHtS8yaIOW?i6>zk1Vww*VhAV+?3Gx zO`M&j_swmH;eDvh0)1%n+~^Sl!#dN{1-O6^6|9J=_21J^Q$xSmPJeYaRNHLSbS-+b zx30V;7C__Ge4>fm&3=7ZemBI&%H6_3o9X4pWvY$6?l*~N*-f^|RPQ-3{1&N1;yS;! zIc{h5C6r;HB4w>{s|vI6y1~`&Kk;WvN3z~NGklFzYXr@bPeZd?su{2A3zp2vhkDwJ z&U*)cTUg~0jxe@PbR-*{p3@}JHGPFIlgmnOP$&EQ6$>+MT9-WM*K*vxra{|Y8_Zr0yzXzhy(!|92K>hGX-OqhYa=V^BG z(xGn3{{GSffB)J$G@kf>kV9>#$jnzN^3&0uNZT6Rzz1T#qqNbh+v@{KfI5cKC>VdO zd%1k6d&pPoC1TQ}9Hk6Wzkx^;*aCZu1f79Evq+r`>dzbf@!xacWm}?J74`ZAMG@`O z)9la=`*!5VCfCtIM}I$Sq7=76ikiPYw)c4~X5t@^GZuH&J&f&J1LomosRws!)VCr~ z)3u{{83|qm9ea4HQvohP|4`>opWMBeI#^}-7XQs@%ptq9UeuwGlcSW*>oUdSJ?iH6 z)MgFPzHqZpBk{cz+e#s?J9jWy={~zyfB=RHVrN@pY~CH7GQYU)f%rr0M3WsemFpgF zvBs9c(vc4Ox&YxlEXSN5FIBwoem>Dg)cqC-Y*1~e+*YFa#_+lGM!?WzuMr=p9(xNj zZ`v$`Ec(C^1sr9PX#Ms6bwfwDTILcEU;m$+V_~VBcqp97sQa>m*OsSd-JJM z;y9CpyB(jTqVYh9=$QFp-&aj+d~FV=cRzG{m+>DZ=Tt9KvPna7U-8R zyn28j<6$)|ZGnkSs-ZWz++St5w>m}Z?<1UgZ9nDF2Nnxpwc6b=veAg+?iDntR4mx| zRC|^VgBjD3*KwBz+tWHXw!qVdg^wtf3(`|no~AIn(fn$^Th#T=mYd^p*Nrg_0!C@N z+heoa*XM_8Jpw!2N(^6;#8dt<3pV_tkCaQO9_*PaH)7 z4A<&=Fn<~NlS2r_E(KFBd~{*ObmbI5EE#h6C9%e z(G}(=i9fSrEZ>0_Cx)t@Og)L?WlW(wOA;$fb=%dA{eQwUzTyL*Wmx09)4T72> zW55v zt@ZQ|sljl;3>p9+c{7}CK{JC&RZq;TaCtM-04X#YK;*zP7e;%tEH(l>a>g~*Xt@?Nr_NJbgjoE zVT%q1p_AF<681TF;Ekp{sN2&Z(|%Bsrgm>2T?ijsM7BTV`@gRvf4vL<#vsY8d@lxW z@NX3`MYG6jJD$3ZTyzYL4xa4Sthx9QLRwjr`b3rJ51zScXYJ;2)F_k0UPN7XB-KBYz+Fk!Lu4oiWjh zjo8~nyU>CSa;K&xd*ncWNL8#?x=t%G-&?u{DpIU+w2j2Ip9tZPXn~R~tGHi&N(>{f zD$GA@^f$)ew7|hXvtel1v3m-;e;=$Z&~ilGR+KX0iV6##AKrk=%`fH|A^*3xXoR3n z&Ve}9wV31Z5``_A#cNvzc}^^=W)T}(8A+8 z&yO$|WEKf~gYf%$ZXhI&_MBycgRg2Z&)Arl)*CP_LNH}<5l`B0gmEf!KYNSw|?J*Ni^+*g*K#31i?MV{raZK`&ZyoPRj|3$q<|qDJh(VIm#~+TCXM zxW6}24w~j#`Sn)7v9zVG`S{&$9@PD4;s5)7^w^PD=)uEy{XB%k_E}1l*xBVoS*nzu zaOp(n|I^Bu1~rwXVHjvz1`vb>QCtYhq6Rby7#3M<8kJxvpe%+hiV8*~0U=>ngLW&j z$*{{VX)+OMKq0b42sYS&hy(-?NI+1;giV@F!;<+fPIZ<2b87zNR;q4p&OP^>b5FhB z`@EjC6RK-vue*T{%1d)xaLG5?w4XPh&-aaN#*)dAw4(L_$%gH&va+Q&rNvXvCu74) z*K%+|PKlRZ7VR;i7CpmwEvdTmQ4#FoY_%k)ZXtO%>|pZUc5g>#XDqv}y!?lnni@X~ z9k47vSo$uCPm1P~$P{+Z6{q>u+RLJU-JilRv%A%$|6;|F7@U*p-EzGf$l}M9y4E7v zRwqW=;gWBR#F9j6d}0P=6v477N9O zaIMrJAZkdQ6{>r+zht82sxcakJ`gZseYKy3+k1MSuI|sR-KqYWfW^Ho^@t;^2wgsp zVrVhLjnNkXg~GkMfTSL=HvL_V5LSW~;x+vH%}3$75X+)(OdW$fpIsi{z4@OG&$_Qa zF2)s9G-;a0p4EltSx_LkKO>HL;=wOpvv0lz;b4+JI7af*3C#m3Do<_v{Psy`SXvUG zBB7+@kb1?}Hz_TxXKL*Gq6&X)$#qhaIhHi_%*;&1bw5@GRIk-o-+|H1-!vC=NSnV$ z>uBSw;DVxs$$<$)etPmQ126qdE|k^kq6?5;h^4p5e~mYShYVRq1o{i&aWy41xwA31 zZ`-&)!Tjriey(a@rb*eZfmAhauG0ASo6w~{Oisp1IU^mmAIwzW{O;Rtu4mybyk36m zy6DF1zNIzm^CK4U(TlW2e(*t<^~T#3;-(?H5jT3og??7BH}Qtf1X!C>2P{J$!)PP6I%~0M{FKS`1lGw%!$Z5yWHY z?l29q5dl%JjW2N)K!SsZ#e(&w$(+HEvD_{&Mlc}c;zaaDRC*-PbhdarzN@dohaEl< zqt_pjmTMMrTG(5bAH=8h;;Q&mH)yAt^+pbU#xELK}{ zUJE_9m$N*IN#>I>fa2%a?6%8Gg**X0aNFbHhV3&F%^m!B)efN(kD3cXP$x)fln^-e zjZ68)0amtlb{VpAa*Us%6>ylkx#Wx3YzDwAl)k%%8_#LIU=d#x0x_VAT3U%LA`jJ{ zc5raWUKFn`axi?r0@_LJ4_Us2kgV05m3RmbZQb`J++weR!5Pszl=%Meo)kO1S>f!L zLV#>flyO#TgWzc!Y$5LKH}QCUCC%_+i(>2eJQX>QHCy}787~3EgIuAyoj+e-2*_)D zb5{H=nPci9@0Exk#KXU>EDSJe0>`eH)wTp=>NJhddrc44!N^Fzmfd>s_uudp4}pxt zqAK9idr>+HPOS_RJ1V#3RCaT@5)LYz6+U+I$FtnMa(kJ1p|X8_ebta=!M6MVJ{{<%rO$T}>EcN=z!!F|;U zI%=N9$s_lkl-hfE$U|sIPEL+@)!6&@Z(EWy{Ve`Xr_gS zdh0hIm-^f*vvRSklIzl88Hz(-VH0d*UpUP>d^c(m@7B7&2Xwn2d`&RzSgpcA; zB|Ebj>&*#fkDlq51rLVohRVTA2{}S8S?lT8>@ZU6kNna`pmyF*VAl1b|8j%~uvj$% z&D4i^As^I`D2egjvgpQPZkREy7*Z&L`IOC`VUvgEruc5t!&{DNC>{aw+B0I_35nIH z8j3uu9yW0hHfgTY1=^2c>@9yV;;b17kb=)NM{I#q{q+d-?$ZTHx9+OR2>5KH7v2VF zUEOTHCmtWe6%_p5MT z%>ML{mH9GS0(+`f4h9-Z0s+GAmz6aqM$T!PExfsKAoH`jsq18&vq|pDUs__LGLv*5 zKLnArVTz{>NL`lzjs}Ln4J_)6y*+ngRl?rgoj0}CW?2c&$|U(w_?K(&o{@o8PSQJd z5-*QYpj?A!{t%9mW;wJn6&&y`!HCXcsQ~5ev8=3>cQ&sWJlYMwD5o|kNYt;2eBY0b z9#oWxx#YQR>sHtAw_j^${13UxCPv%v9I%i!niS9~}g+(M~O|s9=bfJH*{HOdGS# zm*dKIW;-u=d)t!h=F@{lnq-J=>HE641wS#sM6zVMsuIJftaR?@&^RobI8r8RZh5xj z#M0t)29zZB;gYq?VJqL|5zJ@s2!s7gpBgu=t$k2@Y%QsanPOykCa>SHMFW;rVn2)6 z1^)N64!(b*d>21>h(ED9{bcq|P>+)f4sK+tfp`2|=W&$afUrrmJx5`?vAq;m*a7V; zc=b`)h1JpHzu{QF>SWFoujWiGEJ|8z{(gBied|z>fT^9|!EA&XY&eUKrO&-jv9pEc zE^WI)bZUM+24xqZ*ZNCAz0e5WhsA}4Z=R9OWugb6&{4LNa|H=HnLa=Hqt)- zeCk(nw#LEfqp%Y(^=khI1jY%#jr#vP2*uKY&pdJJ7VWkWA}N78ddj-g>fDY00Mq43 A00000 literal 20382 zcmeGEbySq!_dg6@Go*x|fFe=~f`EjSNGm9ygmjl64T3ZX3zuGV>PplUOcV%$P^%~_ zY9R<3gCM9$a#Hw4toQp*_#k$_uA)N@e;$%sKZE~2=BjMyjvzFaguke4xh#7GVMkOH zuj%+CEdKC+IPfd+=wQp6>;3yD%*Z)za_(#Nv^VIfanBEWwJLO$_ssC4dQSRFqb=h` zn>c>wf`ZDkRRuTL-ugM1){d1L7(X?2tWds@XJH)WgkpHXK1+N)>It(=^!t^A)}K>Rj)AXO6%b+sS!G1CBM7FI7C|A1bO=I@Ad{E>-&g)`ON@(qjL5ar;8LbS=<%H5 z;yWRh*8UO*lES9iucNK~nwOXN=bCTXshbG%Bf=AadV(hu%%zziqpA(i%;^n({_kG=`&%!!3 z|39^k!Rn5V#~0d$Hm&^`5hpw6&T9>SmwTOPh?L(aIME=*kCl~Fy=4#dZ+y5M@Zjcq z-jLaQ{>q5cw6azFvQO?!- zpnBa*o(jbek{*wS$f}7@v{B`=cJ9;Np5arMB8V|R2E-~a9M-QG+&-m*PCrigIZO5o z66xVvLWXIkZX0^CC194*wR*@_pgalf^5k=5`zuGJ-S3X~ zAVe)q%~!d(x$AWk2ZcjpMJ6>`A(kF>Pp+p@MKF4m?aj4)PR3^==(!J-sg6iGn%vv9 zMPRbTs&(OTZ)GTkN836sGw^UN86W3Z*o!ML77N&2XI;_O(MeB~^5$w3r}FajHQeH{ zAdJ8v2l@2WgGF!hK5`rWbsmGmU~_a^f&4t zD36UP9*cfa9wzJ`lglLEbw!y57Y6eNzkkn%A)DRz{`*rUoK4Q`&6_v1zvCRsty&}c z16N!n^1pof;&NJIVsGO|mIti-rA<c(jXiT!gE1qaA06g zBWU+^*2V?^IKPGjVe&XfPgmg?KzuZVTL*nu&8Tm5P$dYByT?W#ty+I+Vp_S4j1LI1Ox zX`_|S@$KKf^-J>(zp5be>C3j{<$qa4beDx4;n^{|qsQF$q>JXM=;M&H)gRE+P9wsm zzQ5iuU2#c>i;F9F802*xF1{w%xc!0C{?Bv^lZegB+`POCDqG|JMFKTHPpTzJ8*-~B z^w0#Z5-t6x)4wwD&|N%;+B(=iE=&wVr^Jw*e$i~|mxKfotusa&#D9^G!hZkE9}4!? zj~}T)&mx(xerS1d8oxT#%qt*(>v%7+G5tbtVTSFfn$v!wG2rh`bW7JvYb5@yv5Tg& zbHUA*jHx4KwvIDj!W&$eg8c(}m^HPHL#Y&)f^Ih6dNIu)AMu!oU8N$tAbMIGiz|(y z@ji|ud}Tnb+RE)2CI>ZHHix*F3z|2Ye0wLTG4pJa9=E;{JLOC9Ygace8##%fIERq8 zr@4m^@k-A*#oLO_YeXZ=OhTN6rbOiHp&-EN>}_uUJoHz3Aw(s2=^_6`{^LA?QTtKK ze2fw-DB7gTx^JYkD!e7)!9Ap;m|7?F{JyK58c{z5sk##KB+HQBBA+6q;@Fcl4C;y0 z2T2_Bs~9Q_T~`VN)%Re8T%!z!s57(3k$jFr%br=j<6h0-ID?J~z>$5FkI-mdbH8wb z4^#dqxc>w)W!M}0-k{Y=4RxB>UG&oAJ5H9VdJSy9L!IgoQ^+}$Fg+n1WXgkB-UpXT zSMc29GyUqG&Shyqw}1y6`Xa^ag{`~OMiaN(=K=?rM!C&HRry?6n8;KY z!q_o?D0-1+F#&BuM?PsEMq%N$g~YNwWBV{T%G$Xdi_tlUd@B2tC^())EDx3{CLehl zbP*wSpD8+nyM&%pgn8;ofB`zNzcXrxlDzbw-ux7Funm8 zcY;PWh&_56-7i~kVib#*4Q$aDgQ8NfC%k~Qq}b{sH;TRv7GDDo>d!@6jUuRmF^63a zo?b|B6?57UAGC(PnzhdoFF#32kaP7Kovbz%k53ymSr*u4P16-SC#aUskVP^bnCl=~^61cDwdIG%R>fL2C`rFjfdive`X zPmuq#c0OPo)nQA)xGIXTnb1GW2-6fMOtYvAL+5M(J~1gK8`&&B#`z8*pM~yUtM)Ih zU9wO|qw4uiiM)11{NvD@WygrSU~6($D0O~L_LrutF55cIAdmI0)!XzdM7%ap;8F== ze_b2+U+;5&B8VQ0#1~%vdBek)nnE`erL_W7~ z+Y)VCL}&yiyR}bvCyLEU#v>BP;4%97KJ*tM^?=~+jk95=55;jPbW$I#byqHDFCT?R zrgT7nTIRFJ?PyEJ_N+6alPWlO*{0BWbHeXZ!TVt>iyjrzflqypgLsY9`3a`+&i;$#pQc8zNpFU$`6m%l`cPW zZ@f5V-}~XlXk&nXh*E^p)gSc51HVLRc`=x*@*93&*P@+BJMlybZ44;asZUR97jpAN z@3qq>dg2Us$F~}@9z8-g{P}WLX8bTUl(mmrZb#*Z?`)#vY=XPTM8MYkBT9ju$e{fh zsz~OMlP)U<%oewwQLqz))0lNi=#=9ohQ~tU4IkPJDPrpdd-MtLr5Wo7)OtpB? zpu>HaA%25m`FL^Xg*E@sW=k<~aj!b;`}4j1MW}Mw17rSkh>yZ;IbIB^izuYIh6`~S zt5hQ)Apt*JyMZ&_SS;&y7%jJ#8LzCY)XT`o_)*Bm&#&w3?3^r@R$FV-($W%k``*2K z?O!50N_jB{hD6F|8L`HryJe(NYV58OL&Vq|?X29)%$q^G{mJXU1cMB(jQgaG`{4^> z4qx8;_BJOq^_<&CQ&7;YCOKj5@4b;clZdkgA_(j2o2YbslhY9Oa|6{iu-D`gdW;wy zC+%CxX}`HRvir^a@Ixr8pK_$k<}gncKOdjAsi|odeqmujM^6to+fDv2mdH6PTFTXz ztj3;?B9_kaL>S~9#Zj(ILpq4urBg3>n2&F5d+X}z8Z`w5RN*019Ky$4tJe5HW8t!* zfx^>c#Q&`$#HtNVvYg_A25aQGu!`ZPXKmn)3+j)dz}%wx*gxbxln2-;3BTw462{Xq7v3Q|W6A*!#J5 zVPbB!D=`mxpO2r|&ul#0C~|C%;V~++3AgEpeN|wNI&&)Xgeb;ry@*=XmJle1KDdGG zk8}s(W+(0s7f;3D3mZ0#O#PQC?3qq13YNAc2rDCyLu{bM~XWb+X==${V*R@ehQB;}(fD08)K&#))NA=kf*+5r za4+f>7g|J0ZCZ)28GA`GPc5TvbMEprCZxI*Hji;NJXeI9R^tl;Tsd3L5sWE>u7(K3g@+aUB? z#ugT_nc6ZSS?tIs_mgTVdvu8V4CO~dB6uWd2>HEDWk^Rpb+L%Xm6b4+CK-^6J+cm} zlqz;IC_uuLq(|uG&zoRcUD*i_$ju_g`Fga(^4PORV>8;=EMuaJ7vvO6Nomy;kteh3 zTg(<=FuOZ|3 z;5|wE$ZOw#tv=>VM2G~(bo7Eg#{V$JvEg0m{>ZDd-B0+S!?_HsQ3b1bY+j9pl8F!kuY!F%bFj*3dyDSjhD zQoQM#{Dx1zUFfCDyE_+SQcrNiM~A0q9CgdZ&b;n^c2pDYiNWMu5XY=r9VCCI7fR@a z{NjziS#R`(+k{lv?nkYs!~Oe1+_Ha;iCwv}<=gb(!+93jhZ1fhb>8EPcKgoyi7KIDN-GI`RZQz+aS~o;LGI9UJ1)FC+}9(rM&FWZ-&Q!+5t46=cyF z{yxR;zaqSt=$O-AHQhPeeX*Et)a7WqL@q_xxj&y%OoobyiYiCnj9?y_T~bey&iOOb znpaw?^sNJ>iyj`pqa&@|6{fBN+Emb-CM-i61RC3^4ypQehL&YrC^osxx?UboBM{QWHDsw_0Cx+-4>jiUh6PrY}}Rh`)E`BEP{>uukSuOgdNn$ec~`R~KyB zFZ~k1B9*fD`~RCc7gZTnI=yFNVmkTmU~h}Y^Xu0~s3O!N)pLHZsf9SEuqN>*!^CJB z*Bc}~qG|!Z_{@9Dx@GgI&x+eke-0x$WQEuP(1(6qebYg_GXM~eK2F{Vp*UMCTd-vd zd*_$89%4^KRBil1jO2Sa0`WaN_7-V&kdc!)wvy_a|g`9=1L2hK2V0#$S~U zu?fe`mpuS3RJhIDRErnq*KWS~R`CN-Ji0RTu^kLOeuH^sg}zu9_0m<9lNs{ko^th- zWpLVOYXce{?mXZ(ZVgE4&;tK zfvqb5on$yR3*w+*^IXA@W+W9S`A?JSWh6n!OP8I{RFIv$; zWMP>dgz7O;Vh(@5z)t+h5{490D21U;3qlnO#qjYMvLKV{IQRPOWK^*uA5X;N$Ni9f z@r1Ep4-yXFL=R#zD>utu9+TgEp9vlm7uwiAKdQoh9eL9H?=g4GUJXj|7Af`%!O<84 zJox;}bk?m6Zwl~;F*7f^Iq#sJ{A-cO8O3oS5A`a7k=h*ft`v6!U}Mu5R@qei5ws)( zhOS+{8}&B*kwqy}k0JxD>jlR!CCiIh6Zpp$)QN%#sTvQF44J`med2uD;o){kd2Bcn z;E==7b0o?nw8s@dqEzb;8$+fKnN0Xe2CT)AIL09QpSgSrzVyu8RqWdjU7i?l7m1{! zL_H@7D*XR_Rq8C~RRP$NTkxM@!^s`GZL>NJBW$>A2NB}Z*2PSvxF1Tn4c1$|F zWKBoTOo*S(bbxb*CznWnkQJ{D5WE@_tnq~sQddMtMy2iP6u5}*JXcAiv9cHlHWUal zc?Nb^A}FMRHZ~RIf}S5G{UX{K6Gnr5!o1N2*elvMkhTdi}{l)??3SKwd|xYuJ3sAIciF(#pZJ`Q$hxJ%h=zR#1&%; z5X}Q?#y@iXR7|d?^T;f&2`PnaD6V8#BR?OXOZc}SP9#0EB099zEbNmvA1|RU0$O6k zCV(PD2il~@sF?H*v<(RA3X>C!e9kuk3y=@gI;8Iv1Gs=!x5&hxu9ToHdf9m~xDet0 z3a`cQ)$!wyH}f%K6Qi}UkFya72L5M}Q?KbF_N3I+y>~KUtA5HMY=+M&>5QH>4)^v zN@jO1&7Ms}xZCK+mQf1}p#-Md`fTh+eT1Vo_QPfEzzMASNmOYIKW zd+Ag^i}&}`+ky=`A~qcmNmo=<@bdF7d70YV=Qubz`Zs7uFhicD{_w}}u*5ow!wOX6(vjN6Ws-gB_CQ#RRLnMjUy;}D64A?c&iAy<;kv_Ec7_betK z-A75QKK&fAi_3JWTM{SW_q)w(hbdun+)10V|2xUdA`}W%=%`3_b#w-n(6hz3S*i;3 z@Br>y{{lA=w)Ydu{S50cWX7mV z?JxAnB=2U9j#{--QZ4Tx#Waq)Lqh3~u#@8_VRwX?XrdlMA!MMZCk4E3PJO+!o9{wS z*2$`os{6^u*<=&8{g%Ji5XQqAwC|SaJ?^`^Q4~}l=7?LbaU<>$AE1~l#+d7cl4S;? zyo-~*kpwgEkdAz6ku-c})8l$lUkBL}VW|ReTD!1rbJWpusbVP8Qn5I#@|+%_>;m>8 z+?mC}34iApAA+==Z6~F@{WTZqqhdNf{hWfmhfSdXx#dKR`fhWP-{O)UNpUy@J2QY} z-Nn3^kJp%qY*PV25@=(=Dik3Pn3)s)`%~&QalrdJQ_d31c2xAgEknp9I3T;lhKe39 zAxnw*PVoNf;N)z}FW`{lsp!~NCV8xymLbaew9vnXz~xRhYn2}3>?Ml=m<}C{PfBdTB6{A%p3jV4-0Q%YF3fEnEcK zPe^&1jZ#={-5@C@Bn_mv(3(e+lYE$K@FQV}kHLS?+-?!D$V1zS#gipqV6w^(&p8k~ ztDQx;R%X(#qsnoJQRieJ0Fb|7A0YIX zC~zkr1WHVA>@BcuOp@Mf0LY*0#(z|HK~E{1OL$CJ2N+5+j;vXTD*S7t1_O&@M7p&J zQkvGtb4i8Cl##h_tHezE{5tB%V+a7wl8nF?Sqyvl3ld z5BiIXzajIlsR9($_Uj@(OC0LcIYLK@l(FYPH7%*xCm;64=FmtIS+vpeVPun2F@scb zpumJi$g?IHC0;Nf@t;^6!dLCtTt8TKo@h~R6H0&8j8-@dOdv&P!CKIPK3<~w&3P(7 zn}j%FURk?7a`sH~qxL^E2c92|$QXZ1TSySYT1MHS6CY8eSTeTRPre2WBH^O$S@}z0 zW*9FZAd1~ywOR6B4vUEv=OPvU&7798zcTih(@XwN_^dsr+v+cZ{b9 zDWzj{Enn9~ppHTcI+e=*GyusE$O;xBp_0(5#UUY(E&OXCU*ng84B2qB1oacZYF=rv zA_Z#%@B2?4I!VOmFT)$AxkCskd}&D#l8!z?r;SxB4@M=-*>Nz)`w$?Hs6sI)$Z$Mw z>4^wwF`3A@%QV9mWoowk_2~bS5 zaZQckY@&CPUcSNQl)a_zHQJ#xZ0aVrZr$n=y>6`yBmTV`5Voc;LXMXB_DbNaJXN!D z+g7=(i5BXKcKYwMNNN)(E7tyP=OcUA6diq%noU-k8uAOyF0{HPhaFF&&%S}Xq7qZH zWZk8E6Q|_kGW?*7vhS(_wr47iV$y<8edQGt93tR-4Trn|Qza78?_IxPfK*TxcwL2Ve4V6i_N-)r1A#{F09 zcmDp`?LOM?W?d_<5D=Ut3J?Bn3BH9M1pQ5zb_kAKHP5?Ba$!%Aikj1uy ztM2Qgjzt^0TOTd%e6(0F6NNPgb1R0^%;;y+iUI$$DzW{L;jnKnrdp{knHEU&3v&3$1LfynsPWJ zcc>7(KKYo&0W#V+*??MD@%8@XL!F@QX14Xci6dF#CUoD$5mK%V1}ujY8pYf~fw^!a zH6}m<-+T|k3c{kh0yJAY{haTlxv}GhOZ}gl1wcK3b;kkRdNPPDOmUqc144g)*Ve?` zf2%gE_uc-{NcNrMGD$+aDK>$ON+QBOuFVw7lZdtfZ_xs6g7k$u$0L6l)jYUYRkTuz z7x>{hWNOdKun5YZ&_`iW@7FVcORe0rYga#mijt)+SaOK`GQoF~fmmxZV4jFP3k-(z zITi-&hdUAIM8RDh6w0_G^#n)0k$-x| zLm#2(r=%^e@Yq1g4cLUiRvshAxXR#LBw{itQp6HHpxA;b$iOp7+>ZKBAT)H_o5WRh zT0#P3mB^y1ipr`&e^p_3MH7a+#DMMEP*jThqvV1KSV4Z}iP~5oVL%@2<#CZi8w(ze`JX7i6E~|-@$X0!jfQR!@ zPuvMwqemhX`$Y@`a4HBwtfC@oX%=Cu2EY4M>aH+VPPs&EWwKSnmtsmrb%I8~s-=n+ zqL@#0Xu6pda)Jr=vf}S#)TLF$tIFS{W`|rDfqu=PY^3$T?W$BT%D?2TjZV57(KDz) z{6}l(!0Gy8_bs7FUKt(6ia80D5<|i?Ug{zR;p`B+7+;N%AmMC45sD!Y-&7|JzUKP@ zaf+MkA~C7cg%LDYCrtOXJUOOV4un@wQef(QMuI~9R^>Yhp9F+_ScS_Exn-(B2u3%d z+?#hWl)Tv{n$dQFgvm6jag!!m)Jt){^A{ANi{Ava8gX7oPMa5&mP=tyPIXmMHQr>^ zz;~lV0IL|H@K)UTiTCGJIoj@35vJlOkh?#fZ5cgXhk1`iFQ!7z;GCnQyL&x&$bZir z)mY?n&Zlk!-6p)>E{0l#%GcnhrYqm~TZT+GY2vB1vC8kwG+gaUI4@<_oY>i?c+K%(`9MGu}64b>$_l zt5=dm@#>rTW)#-{Z4aC{>{VhL{2WHNfhD3*j_$fT%iI8Yj8Wsmnp2kyq6#g~E8v>k z6<1iC^(=*yhPuzZD@m5J9}pT`luP}w0U*?IFYiFc_?T4WVoM&T$#+fu1BHbD{)2^ta?(P*ZB7UBYI_s3b6yz?Y@b7k~s#qUf+6mF!oKdk3RFp+i)LoL{C?D9;gI> zeAIg4b#!$H@d5vE^L24L)=PMmH*Bl6)~$Sg;+JC6N8_e&S`GF6@p{8Do6f#`!#_XtZod;4%r`7sFSs9YJKAG=bxM&jCLI8a zg{r250sqJtFvopE$@J>s&kO zSL=P(c1>)1`+3+PmgN8|m#Fh>dCCM3x)Jj7+PS_wse(C_^B}ZHxT98rVUazm@3@a141!nQ?L2`q^x7XRoC>g5fLbS8Rwz&h@YmkMt?s+7F3kIZ{kV5@3AnKIK`4Q3( ziWHU}O4_die#j`Lv&CX4%tjS8MyZXJ_i9LvMgOL!l~Ptmm67%ckYfZa19Cr5zm!p3 zc-RjNts6_4W~&%gFub0V68sj6oWAtTAmM*AHN{qH8PV$mVv8bTa~jynfw7FqsyCnb z5X-|r>Mokg9RnVBb;~~YEe%zFZ|~!_`*Onb^Jhq* zZERd31q?zQ2JuRoBH}2+%WXOKu1a4e5wSe`>q_syxj0SDgkP_)pXiGZ{;#Q-s;*FB z#sI@Ob59rPW-*R|2cqqWfayhzPPyGf`Xw-*Vlv`!`12MCI>xV0AmfTq2#*HQserq$ zgn@#i=+tb;fcbP57=VtFZ*MPe=(zNO4irhYN%vpTq77&bu+J*?q5Y3{F4j>yuz zkP5nOi^r<3!`bB}A&4Gd_>MeSHsHm0-R)GbP*U@Bm5k;Gp`_qH81Z5@g)`FcpqexQ zyBcl1Ny}ca$?B-7U9IquIhm^r^;xgB@k)g!EwY3`=Fq062T9#Qfw+U2le+{WlgUdD zEUextL8+?ZaAooD`Cr?<{N=9;SJW1XE6I}130TOye{@wgLZ&qJ4AfDD&)1-yYy^4D ztev6wY0v99RvdT5AJQ4J9-CAbOYG*)3|jgYm*Eh82USjlPNC$?#Ow?>8lTpXR){YD z)3AFflGv3SU*p-9{t904vtrXW<#b-17?;7SWYO(I29?-qC**&A8&l_HEM@(f765y8J#UgFa6eVd56u8IqRSPAE=*-V^y+o4qtm}mCHjUlO_kK! zg@%#4-UM3!10%ZggFVwGf z^+PGZ{yVA}xcwv6~vQ-y0d2ZtPNXgHU}Ja~>5XlqcG|yN6zD$d2FxsQ%INr?rfX;y#|p zI$noKeCj16?;v!52P9aigGt`l^#<}%@%(7%y%)c`las5vX4mWw_jiEDrd_lp<@-1N zOC(!Sd;2X@_fE0V5-TOHkoX7AqK~{{B2~l4u=-xCnoj}^uqowAR(`*Jq(!IKydFWO ziE%e?fGz#V5SH55Xx3}D9z$XF--<bO$lGp%8K8Rg(?NG$kk1F`S?U6rpFDFGyjrfA||)a z)cV=faLxxEZruTkBJmn!%2sx zMkN=QLWs{ZOzXTh{Nn!?Kkq2&OCoPR$|LFNtgdo2u0dBxTgL)3I-0h75V4JK*f4o? zP(nFu1{E!$V2n+7Qo#m(h`aDb%Uz(>H9);%cSGhwZ+4~gSmxiqf6eayPx0d}8Sar0 z4v}olZ=Yfkk=E#Zm2EaYvZwBeWQ`WDWvhos>1ZTs`*rGJu6yT`z)$eU^)0(-ZxXOYDBHk*I$TSSk^gBv>ft8oR2^q~i={H{> z$sZ@iHKu@eTFI!w58zn)pNRYwvwA7>b}A((-!K>WLJ@YVR|Ch(O2d>UrfLN;GthIF z{@73YkUai@l|Pc+8fuQ7j+di)Ja5iDzsov!lV%{F1CcZtYAL_gh4xB#nz&?Pb}*(1 zh(5vWHi(s%m@lM{ST ziAC|`9oB2SvQwRnIX$1>_X#X7(a*IrVYJFfhb%g!`GPSP_OZ^DAnRyO`(Ixg>Mx=j zo_Ty;UnyF7#rb3XRT!K{B7LDJvCpt{Jn|L0d5n#|#b%F*Y)GzXR9? zFF(aXzv_|Xsm>HaAn9tFg_%;V2~0km>hrlV6DGH;neVQD1tZZfy78IK`*w7r&%7du z!OEfmfaB$l+9Z>19B;ksy1csO&a2DV(ObyDT`^^!#s z34|qr$6{;y)AwV9G!{J4!3H)A8Z3UgM0OUOAkELMGrQ;G<(LIyn(f)Wt~t0$L zKA-Xp8@D;y{3Io;hC_cEQ|%vUYFLfnoYxllCT8Pf_Ezah>fO`D0kjB9*Ft8OfQ9$= z!hRY3*FW{uQPeaL$oT^jaP zRr{a<(rujb>Fb+%--DmmwIk?MWhA|dyD6_uY_lSB?bo;N0XSVE zIna}t>@@slX$mncnf_%e>GMI?Dr9#4+Zvaxzl8D$$4XwHKNu>(8jfPCHO#O@b#h{A z(S)6`5g6e9R*Bdex3|$5me}MS-K@V^TQ(9 zZye}8{GfNF!_MuMThWbmb|GENrcu96oEx+wvSI@zD3iVMRN`#^8(DYDj6;`AgBI2K zql2P?e^Edbrd$G~QB^d~+aFVytPLk6@;f%8qDX-_OQZF;`Kck%oI{8u^k_4;SbBXh z`Q==t?)?YPBu%mkWO6SwJ=AVB9u>1veFfu{8B@3(2jLS9Y@Ev^i-EqVCf1#OCkA5o z^>b*lDTZ3x%{V92aH-<{05O*(Rszt_a~$T7!(P;x(3!XcKm z&ax6h%!YNZL>Kc0T^7yj>AzD5@ z)GlDrUG4GrqKsFw6Un~pfs21Z<8I#31>#R7=UFZv^b3v;gd{34v~SF}37Vw8;Ck7y zSsZ2au?rc#g>*>5oHbGFZ+*3IuYa8tyi+*QbF{y5E?{1+D!ZgWW~>=anto_)kbE?s zTz;6yLGgQg{b%^;`DhheYkz~U7{#-UIK+kp)iY&4G$(U3o1j&&y+}fMD%3k-?#Sp* ztrwfj`z1@%LFtnGW7&gM$KAEPRvp@o-fowAWq0n#b{NY!h!6=fG^DhDkKTN%p(@y4 zbU!ig?{w6-)M&n=UVY)l{h@pQ2J%I))B5AL@GGo8vK6SE<3t@|ahGaz*xp~?UONc^ zCmE;=aF3XIF=p<|!*>fBmb5xq<*utav((Pd#(U28ALWJUWj?VK3{@{*VADIi;eEJU zMB)9}qh+Ai&#NsX&g6QUZ4vNVfG>h0SO+o8xSkecU#_tLPR>M9Gvu61!!L|%$v#8F z`l&&AYjaVXh12&J2f5!i<-YS~)y1y#CsJlEWL+^!pRWJzL_=d31&nTLRopLYe}u7q z(@vG;vbNR(N11x=*qNcKQ|kS?vv21M&5v~J@bR8&K~6X&uE@TjLTO_#+P41BRvG(6 z39f4Q@Wsm9kJ|0s0m?7bLt95KnGj6d_at^Ld65j*E{L67>NNX- zS-1h_zR|94`|h{aGa{36K*ypiS8u)?2m{w4cjT|EvF}zjp5q)%AEdS$8mT%^@GiN{ zL-1s6qkct&k}fT0h&5$o{bNI*jYULW@$k^TDr`=3cyQ^*!x2uI^27Un^GlKq>zgMC zf7OaE_Dz`UJSq z55|yGI1)1z!B}8X`VyJ5Xgej8f27cD^O6Tg_d;;MUrA#jE9iULM7^)P)N!$e7Ohj zWMvC<*0SC#aA@Dmx%m04#5Cx*RE64Zvy{?$-Bs6eOU!IsXE2vnqsZB%m_UJ zd`%Bcs!1}Fkv(amGqF5bS(3Cb3eCQ2`?DBB08|(b+$!d2ki3XzLLAAyLQgnX(-u?9Uo?yXgLG!Z@+!2)(3ihFQG*U!Ywa-_X;=Y7+Bu5Ij9j?J@D;o1A-RqdX3|9W6RnfRh~pCL*zlO~QMaf+Rl{6JR0Hi9+N{@0B-JEPtvZoYr8@F`NnElbg(~A{$n|7}_tth!yKz%}VmS zMZh4_B4u+4DJhNYb_xOR21FyCsh_m46ZBHhMajrn&&q6QYPNRF*UL3ztKJTg!Z#M2* zKDPDs+lJ$(Y|C8cSAY!q@(Lc$dFYeGD|uvW#Q0MbgA!8H1_(#;drxjN_ozp9s&1(I z=c9-J(qYG|_b96JK3S-Lcp%mALAk$`eja*unw0+{MB+pfVjpKyN%o53wE4plf!ewA z2ivW3Z#i4;UWPjS0-yQL^Ty4JqQnP#gk%D06@G9c?c|W@gSHEjNBe6Vb7_&Acw2FR z->a()sHv#|-@er;D2yEKd@=jcU{fc1wLvvN#H@VeD;OB}+9^qei zbN|||HJ{MmC@HG7HZKqhH(Nd6+j~dnt^>k7tjpK z$2=)dT_TV<{0^>8xf)ElN)U(#Gqg(MWv46IaFk12wdG}kO^Rl#Xv0!q7W~N5hMUAau=~f_t{sRDX{vC%CjRvfS3c%JnPXz7%R5ANvfn4GGG}Zwx1Z;p5~~yM&`SSj&1O-zfnwGp z*;69A6U^$wSUUa_iFVHV3b*S=;iPe4<)+WmOMxnS-$v6)d>Hs`a+`CrJ$EEei43y7 z$f>#-A?s)tu%9vV%{(sirwG-~!{~`;`m6`C96Pqnxyrd?o{;1GCHg~QYv!`w1R>k? znm~rxUk+fI0ZtMvEt8alW#aftR37xm%5cSS?6|ny*vF9m@j%!Q#8^>)PC8;B@%PMQ zA}AtIA%Drug0@m>L@rinl}@vGM`G#tUW2&(*}gSgHS$%c>2DA{fL9DXE1zGq3u zU2$|giu(Dy#VS}{zRb28x4VI-8-l%4VBTZ~{7IQ{0&#L|>4H7fTpGq%5d_N*G{u_aO?j*IdW|XF7{BC0wU4>G1H$b`A??*Xcl-e z6>#(+U@d-PV5sTPlVD*TarthlH1cpA%f|<5Y*?q1NVu(~t9h^LHgJgHPLZSkKkoZ5 zai=#CHwev!`$TMfIiZoJ{4Hlewh+RFB$=69Dq28KeX9hmlji2+sK9B=5!Rb}uBTRM zKi|E~BI%(F=Q5yP>o{ZQHhfh3AOGj0Ss*gcKF-0xK?|r; z4`uHchNPh-h-H&C;O>gfO`(&|pFj6Wy+y=2G%#~njpn}}QqqNH1^b5eXZ#5_&%mK& z!p$>>_k;GvCb~A+g**USf(nqf`th>eoX<0He$tuqAOGjWfd&XTgv)Dy#r=P}wx;JC zP_I&t(Q#d*rJSI-dzsIom=y+$|yC3o1Z_R-?Xa>1m$~n<(pdzc|{Izgf}-kJIay~PFZ%^ zuasFflbX7IBpjEGgQKZ>g(mM0gU@$8St>AV7^Nac{|3a;etEd0x&GH%;|@5=84q|H zob}xpGQV;V*t7%i)Wk!4K~)EEwc>36e5-LDtE7tjK2Ew@9k^MRET*4tVCnLscJLwa zdOHOlPEp#pxs~`%VuRtKM$+9_BRC6d6>&!7OZ{T$w+_d`8ZL!2&xl~z!YS9``qfXA zo#67@9|F6!L&kqI_N=&*kfVcxu=ImHra3FZ8QMN7(#O9AgLZDaL$Ulf_+Y!Y@7~1| zN{S-zkRQN~o%bKkoc#^Bqtz|pY%5i6AMr!%v&UR=cVt)puTtAGOVfi(I2*>%LOwY zaWG&uTyL-CsZ*!opBXQFSrtY*vr9Ok1SYBv4o2x18Rec4v6(YJ z+AwDbzbPW-WPG;6N)qmEaULnXa9OP$WGDkZpKu;rSlH&ubOUdlU9+?M#>U3;n`!^} z@Bx4>kHcATxVpyJcx+|~O7#8W>RE)2ERX&Q9({m)scnX>cMyu4cPwhiOg>n&1uM|@ znU=*~C=$Q?qi!QqZvR^H-uKyc04sq#XS$mew98#lRfSLTUvB$c?CcD$67C5Z92&xN zuv-zkIn_QCSHTJ;B=?YoxpwEc^GMqXVM`zlx-Q zO~a?xp?uQ%3T{vFXr{h07sj2`ldYLnD#NTF)^eAZho^TYQmzoLFu~X23;VA${-Gv- z`6H}xR)mXIj&_TVa+CLeauW!Ogo|e82$kpMx|K#J8fr#*UX1t^hyTy6X#D(J2VRfA z53(K~SP`sQljE@c`=Wz3>!!*7`_R7Y|G(c~FL~={0?!Y6DYI-b&lkhm_j|v80iFkx z4a~8fuYu>osA#IIyIYsPyD~F<-pqTB)&~;)tp%oaZFBSRy8ZTrcP*dKS$ykOl&g{; zD^uf^zMJpZKF_#U`~B{&hi%d&K-YOsI%oSm2DrWPwRillN&7zb=AU@Jc%D(f-o4?X zqN0I+Pc!^GdMidj0K5yW?(eT7rNE7audlBU2UcDE-&d|$^{TAeGirNY?0Sg9Cw&2C z{`l{A%jd^cJYq(VPG8%W~NI72sybnKJuwd7s^X z1x(+ofMfKkOj}j<0{3wM z>*I3SnjMRuo|^g$c*;*|-AnIyBjD1;yldC4&03Z@%luA3^Utm8Yjq2M1KZno&zwAI z2wFKgD}VprZI$P3zh7Wm#?D+D^M3#TI%5Hj8T)~I&Ci18c6@z(Z~hX^Nng zj&z`+?{9CXgEBH`{2Dk5;`bR?P#b%Dd7XJ*_kA}^vF-Oe$sk3*{P92a-MziH{{Vv{ ze_rLYna|SgKF`eg%eC$Dl`AQrvhNOX>x3{14l$#PBU| zLI3Mxz0&{xm%Y7l%u#{k@CrRfbH`044$m$uoe<30`Feqf*CK`q$1b{-pH#K;{$Ksm zYcXfo{IKZnzxV#yxBgm`?fmn**1bL(T)%Dqn=dx!qpoi~R=DQ&+pxn1S2pRC${e@% z`8V@h?f` method. The following sequence (done in a python session) reads in stored data (from the compressible Sedov problem) and plots data falling on a line in the x -direction through the y-center of the domain (note: this will include -the ghost cells). +direction through the y-center of the domain. The return value of +``read`` is a ``Simulation`` object. + .. code-block:: python import matplotlib.pyplot as plt import pyro.util.io_pyro as io - sim = io.read("sedov_unsplit_0000.h5") + + sim = io.read("sedov_unsplit_0290.h5") dens = sim.cc_data.get_var("density") - plt.plot(dens.g.x, dens[:,dens.g.ny//2]) - plt.show() + + fig, ax = plt.subplots() + ax.plot(dens.g.x, dens[:,dens.g.qy//2]) + ax.grid() .. image:: manual_plot.png :align: center -Note: this includes the ghost cells, by default, seen as the small -regions of zeros on the left and right. +.. note:: + + This includes the ghost cells, by default, seen as the small + regions of zeros on the left and right. The total number of cells, + including ghost cells in the y-direction is ``qy``, which is why + we use that in our slice. + +If we wanted to exclude the ghost cells, then we could use the ``.v()`` method +on the density array to exclude the ghost cells, and then manually index ``g.x`` +to just include the valid part of the domain: + +.. code:: python + + ax.plot(dens.g.x[g.ilo:g.ihi+1], dens.v()[:, dens.g.ny//2]) + +.. note:: + + In this case, we are using ``ny`` since that is the width of the domain + excluding ghost cells.