From 22da9155612554cfff5f04036d66c954e076f52b Mon Sep 17 00:00:00 2001 From: nileshsawant Date: Fri, 27 Dec 2024 09:15:48 -0700 Subject: [PATCH] update to documentation. Update to sod notebook --- Docs/sphinx/LBM.rst | 5 +++ Docs/sphinx/VandV.rst | 51 ++++++++++++++++++++++++ Docs/sphinx/figs/kinematicViscosity.png | Bin 0 -> 23755 bytes Docs/sphinx/figs/sodShock.png | Bin 0 -> 37811 bytes Docs/sphinx/figs/thermalDiffusivity.png | Bin 0 -> 32449 bytes Tests/test_files/sodTest/sodTest.ipynb | 51 +++++++++++++----------- 6 files changed, 84 insertions(+), 23 deletions(-) create mode 100644 Docs/sphinx/figs/kinematicViscosity.png create mode 100644 Docs/sphinx/figs/sodShock.png create mode 100644 Docs/sphinx/figs/thermalDiffusivity.png diff --git a/Docs/sphinx/LBM.rst b/Docs/sphinx/LBM.rst index a83b15d..b77e165 100644 --- a/Docs/sphinx/LBM.rst +++ b/Docs/sphinx/LBM.rst @@ -51,3 +51,8 @@ where :math:`u` is the velocity vector, :math:`\rho` is the density, :math:`w_i` u(x, t) &= \frac{1}{\rho} \sum_{i} e_i f_i(x, t) MARBLES implements both the D2Q9 and D3Q27 LBM approaches for two dimensional and three dimensional flow, respectively. + +For the thermal version of MARBLES, please refer to additional theory related to the additional energy lattice and cross coupling, +in the PhD thesis +"Sawant, Nilesh. 'Kinetic Modeling of Reactive Flows.'. ETH Zurich, 2023. https://doi.org/10.3929/ETHZ-B-000607045." + diff --git a/Docs/sphinx/VandV.rst b/Docs/sphinx/VandV.rst index fd55618..557482a 100644 --- a/Docs/sphinx/VandV.rst +++ b/Docs/sphinx/VandV.rst @@ -24,3 +24,54 @@ This demonstrates MARBLES ability to simulate flow through complex, porous media .. image:: /figs/pine_box.gif :width: 600pt + + +Viscosity verification +---------------------- + +This test case demonstrates MARBLES ability to maintain the correct viscosity over a wide range of temepratures. +Viscosity is computed from the simulation by tracking the maximum of velocity over time arising from an initial +condition consisting of a sinusoidal velocity perturbation. +The temeprature is not restricted to 1/3, but it can be set to any value less than or equal to 1/3. +This feature has been enabled by incorporating product-form of the equilibrium distribition which +improves stability due to maximisation of the entropy and the incorporation of an extended lattice +Boltzmann approach which adds a non-local correction term to the pressure tensor to correct the insufficient +Galilean invariance of the standard lattice at the third order moment. + +Input and post processing files: `Tests/test_files/viscosityTest/` + +.. image:: /figs/kinematicViscosity.png + :width: 400pt + + +Thermal diffusivity verification +--------------------------------- + +This test case demonstrates MARBLES ability to maintain the correct thermal diffusivity over a wide range of temepratures, +as well as at different Prandtl numbers, adiabatic exponents and mean molecular masses. +Themal diffusivity is computed from the simulation by tracking the maximum of temperature over time arising from an initial +condition consisting of a sinusoidal temperature perturbation. +The capability to solve for temperature has been enabled by an additional lattice which conserves the total energy of the fluid. +The new lattice representing energy has been two way coupled with the existing momentum lattice which had been augmented with +product-form equilibrium and extended pressure term to compensate the third order error. The new energy lattice adds an +energy equation, while the multi-relaxation on the energy lattice enables a variable Prandtl number and an adjustable adiabatic +exponent. The mean molecular mass is allowed to be a variable because the product form of the equilibrium on the +first lattice allows the variance of the equilibrium distribution to be RT instead of T. + +Input and post processing files: `Tests/test_files/thermalDiffusivityTest/` + +.. image:: /figs/thermalDiffusivity.png + :width: 400pt + +Sod test tube +---------------------- + +This test case demonstrates MARBLES ability to solve compressible flows with realistic parameters of practical interest. +The correct position of the wave strutures confirm that the speed of sound is correct in the thermal version of the solver. The model +is the same as that built for the thermal diffusivity test. The sod shock test marks the completion of the MARBLES solver to the +extent of a contunuum limit, single-component, single-phase, transient computational fluid dynamics solver with mass, momentum and energy equations. + +Input and post processing files: `Tests/test_files/sodTest/` + +.. image:: /figs/sodShock.png + :width: 600pt \ No newline at end of file diff --git a/Docs/sphinx/figs/kinematicViscosity.png b/Docs/sphinx/figs/kinematicViscosity.png new file mode 100644 index 0000000000000000000000000000000000000000..d60d7f1139d4d9b119425341dae0eab1af8877e6 GIT binary patch literal 23755 zcmce;cRbdA`#ye210_ThNhqQaA}&HU*<@sf?91M}?44Bhi0l=zSHl(>_LVJr@7?!& z)%(6b_x-u=&*yu8e}DWu-fwZa#_RPQ=XspRah%T}1vv>q0!jiD3PmU-DXN4*9aBc3 za4PUl!XuBE^>4#Jy!K-1_R3a9_D+v%AERU+*;`v!*;|=xEs8^yO6LgV4PID z&6}J^Pb!Q@OhQ5uDRn6?5rzD4*MMD2R8;gTP6!hneBrChz@sEm@AZQN0|Oj7=tYnlao99 zV`wNp@hp7E#gwcpyt=x&3eQ6{KR^HHD$8)0i`i4%N!8VSzb}=~AczT_;4Sx$7w;Sw zk#$iiTO}(m@0}y1sHkjHYb(BOSx7 zIZFrk*&!)Q+f(P^R#8>8_RvJMn+}gw;o-yJr%#Weydxsyqb!pf$rL$XD`nsIdl4NS zjoh%AQvSozV<<8*vJ{W~Em_yL2+}V>aY1psN(BW4s-3=pEoifFFWI}sFOZjhQg##V zR??W?Y5rda_hv?53xiBNHzCp84TUHxr=69tm^)TFQC87`aYX+$A7=Q{(%nR~?}~y6 z9P(o>KivW^YSmAp%J*$CHRY<{H8ieUS8gQp`V1EogLkt?6IZ6zrRpo zf-PgIsi|3+%Pkb{b`3i^I#RE7E$wvO-<;3X=z3@R{jJrAK%Pd)Nsibv@ZZ3|fa^Lw z5s`$bsPELWuIKIv9UYyloE(n!`D=9hjrze^*Gz^U9#*((P-zwkQ@d?x43#Znb0?;| z%nepo=lYf9XUjIX)tcoKzBF-TI5Y?2?^yZA$5WlX%pdr(J(`}Io4nd-xmROpW+u?* z6w&yQ(@0~^gWQ#`5;PZUj=0Qw^;Skd#BiCDudS~~CMOf=>gv{F_?(tPYijsAM(D*u zgCm6Vlby`G46-&}Fz5H&2>+gc3i({Uq2eW`VG%G0hUciUR{;ITb)3bnGiiM;!j zs8(V!GT$_b2>qU%2k`E%?#*wntauB*_Rh_{fp%Qz*TZ%u@=d&df*V)=nB$_BmX`JG zw`-)j)FfOKmtyJ6@Pk8oh^(wI3+d6_MWyBrSaPzg3GFHq1H0)Q0~)V9EH|Gf_40r+ zRXlfxh3WPpR&CmfhK8miOa^DPpnBRnT_xRHeUk>5~1F0#3$A zNv_>==XF+AaYaRvhvk-G6gtj0^8BmE(GGKl-}VBtc1D}-NHFi*Xk!<3=pzAAQ=fkD43>8`btQ~8e& zdb;YbbRpR!^x`S9@!Y*~K7m(V*IF5E)~95B{r&ZqepEUyRZg!gd;QWZabF#8tk_?w zUVV>Cr1zu3#(Cy-z?Xc@JBbdTlLcvPD=&Tc@WBM`u(R*zyaCJb@UXR`V^gKw^joiA zp0YM#qE|`F<*vUiKHT3nXbe6tdH=qMr>8)f)rih?PRnxrF9Fp&^(O}h?i+JGr*qW` z#nshg6tM%`X41#cv79(4xb7xc)2km&(dZ2nKSp`Kdw09!XNkoC zN6pdRypbjigkCN-M=s+Wh}Wj3rWfwIzLRs`nEAT1*{91=+|d|9k@ohj!g@<!#0lH#*?mHYk&)TEeMADUB!zH6v&r+jD#~VXV78-w!Vdt@(#FLhmPElbl z^4NDgLwQ@$O4Hcbn1_$=ep@_p$!yik8D`yeE`&UQ+;meXYizl z4Doy&YFb*M^}miZ*PT>V$vHSUMrUV(ShcIJFf$VmO$|lfbIw@JE+a3< zE(3m+q#;*68OP_u%E8HLv$OKZJX6?E72wD)ZC<10?Oj>^0o;FtPTN(SrZZ*tb{XV_hU;-;U9cXO*GSZS)}H* z4i^E`@F8HnPoV@!=fW*Ll;`eGikrH1mn0%CH_8hh`a&qsgKcj#qGP~+_wLsW*?5K# zOOd8<>JZqQlhe}%NbIn;Z@?wuSIA6+TPw2|cnbN(O#kuY$2)87Y+uvwQ?$C^(LRfe zeDb*c6|JBkO{QGp$G49@L1224mv<8-lSgX>!vU)9vJjHwOzZjnaQ^-`ft6%V96nwI*&d%Wj*!zVj?%T3?q^(;PEBr!3}2@ zmq*=Mir1N#Ozj8qGk@mf=B{pGt>~{@L3MU^&U{NhJ?uO|m}VQd7W_Hw3ta#I-99poYq zd;nppv7N0xG$)6-LsG-cj2ZHZ^>m`s9RD_K5j;Xd2x6yeYHN$kdYD^xyU8IQrQeU% z|Mo_Tu&kP^es!YR0AeoQnKPn#dRLI^dS^LQHN1UU5py%+OoQ&a;{BAu6tiIAxw)R) z9pvU=e^pdg_Gk>3VKi&IZZI+3gTXmVBM=Jf@coAmIFPACVIojhV1gn9k9dqe$2@_@ z^p=NfV{V(C^Cc!^{ca2)-+6zrya?90p^*_CE$zK&Zq=}#`j(B2jXUEZY{bOGC^oi2 zbUWN_INYt-(a~g-{4Il~Q1QNMSmEN!4sowuos5<7uk+Y9(}x^~N|}vXuG>F5(Hs%r zUJ3*&kl{!U%Qt6wNIJGW4A{kTcyAF z@o9*4w;-rQmEU6typl#nPELSIef#zc;Lh0!m-UR6&!5TeSPj3ZDd_K4#}T?tMt~y( z+cb{$(;hXyvjU1(N(vd(=_EH%Lqo&b_R|Ce$5A#LGm5mn#1NH_c|+1H_+IwT($i%PD5h0qFNLgih*;g#iErLqB_Sj09yN~_^s1{k+`@8s{etMPKl=VjDkLxA zbR~RRUzlgI5s&S}jE;`@JAO*pWV-IFehjJ)9$XuIXG$+{v@?OmR@nXw=JD9G_lfxS z?c2%q<5wV;7B~&$@8Z)oBzqk({rK_2SBsV6^5x6p0W_W_)`x~}l<2b$*Vg0dl|@yW zpI^J5qHtY7Q4t3~ND)ak<+xYcy|A#rA9Ik(VmKZXSk@mpR#aD~nfDc-yw|6@a8bg- z!qxyfQst8coF^kCU`w9_G?oALaXH#%{CHb?yAAv`E+K>$BKfh47caK7w7iGl!|u9a z;u{b!0a4#!(ACXt;%I-l8S-zD&4l<5uV05NTU(a)R0wj?^*Zd)i-)9XXl$J3FyCu! zV{@-RyVjL+iZLlo1*O_SDPLyf{4Id6QG;9{@TIx8U|0h84? z;GdbvcY zfCP}wZ!`wsIWI5oWb{zmGor;NTlQ0T?%a{3C84AYh8=d{+O=!&CO*SkoJODDLtuUD zvTo!fw|bHJ?pt^Upz2%Iv$%3ez+~y&J2`f4o6;tTL+La4y#ha z)fG)fZ~|7=@Jul>Ea9s84}08#PKwyE--(qudFs@uMBT+6Zb)4!Dk>>}YGCWLw{t_1 z2$}CIoT%Gh5)l&eo;r8D)C1y(GbAPl{n^=B2ZHwPVK*Ycf)eo5c?;vMg@HK5ioMIY zr_SD<<(g@ir=dj_lq_;1+)h>ifi8PQb(kMfQB&uymsi&QI@-U$VTdOzB0@q%6#~H} zH7ADvikzlXH0}z<-E(DFOP}ayirzw9BB%)xV5@a`^LKmHICtIOd*ir9Q;?C3 z0+4VhSbhQNac<7cVz}latf(e6T@ODZfgMzKcD8HZ#Oy3S%$ARBL4HNW1=u0wQ_GMz zAy%2}71bT+cmVCMd|HU**?{pQE@jxX=iye#o;-O1 zRXn%U-0_V#sp(BwZYzsQUabo2d!G~e)OL1v@%@Bp(N0S@FT1bYi%dwk>6^tD1fX*$ z6$a(wz<|aElvN8u)oRAyY>bRDE^r!!0>n=QR6`E$JYKs!Y;GT0hsg%uX+NtpHa*P< zC`B}yRl6bHs?OA!)n7*Fl)|(kr^Z3=Pt4#u5}%WV=O6>;mL&Nz6yM>-@t9}|OMM;T zFBU=`1PHUQ`T!u`IRLmlL5X}$awuZKgB2gZedP+*$&d)C57}w6pjoB0wUuxyt?s~= zSidRsV%#=l*Tn|7nj4ObgO9fkT8i_ya1t9|MDCai6hBZ3W)eA_o^d^^r7ny{P^SiB zF1pC7{!rD~)%770z4qQ(+SI3qjb!rcB&F6(S`sV0^MF51|zPdx4ZFskk#@Xd34#jc;8 zQiXr7S%{aGmU`s^B-8Tp&S+?8%nX*^fRZ;4=orN)^BEqK)-Ed@$|ne1{V!UGgV@R* z0lIr{V>rG*_>7E=4FbH2RF-sf#k8=);p?j2f6g(*-*Vrv{v3a&YtZyUSmTQ)&R^xs zq*R#CZ!Halu+<-&L+TCLurT}WAtw=rsCx`iAAsVyckiC0q@>An48$E8$lDTt=K!0o zble%auu%gzP1jDa1PAxzm21};-*ERD&Y5?$w|}ZX+GGQS^`WQd{NdiF$is(mMwRvT z^+sO3p0M`%vJM~uBEVpv#JEp`+qm;miy;T_G96hC+#d0tLdLI|# z%T5Q3^|ZcTP(1X)6DH*hy@o)NmF;Z_C8g)3$|a|T@nCx(MORB#S97*;uFch>=t)lw+sTVnWxgfZ%0;)1Q8ZI%8SlQShd-?JuIX2^30>6tAikOLs zDXxFS>lX_u5D=3%9<8XV>MPWhx}!aF5^{0|G6KIHprJk|okl@zxFwGD^Y;&iF`sX6 z-I$>zXL)F-67O%nJgoi0p-(&Y-8&W1+=2pfr(u^&=h^h=GY&r&l3WPhXK*>!h25CE85&BnL;;}ES$snPbtg|%C8o!nHK&5$myBhLkH z(ig&Q1%&%XxX1jc_55ZZNU`F5no%PqT>>}g=;_5Epor~nEh1&SsEtig3^yjuWT>hX z!B&mo)O_|yq`;EwYNjg%l@Gh#x1aBo1w=SDKY!yDtM)@HnmhKF-c3&5k;n+d^WN*$ zF*LmDOYHEfrvZS4agnko#36AVousK{ge0k{xdaiKd?1O58#8!3R;KqiT7l_8;v{K% zN5=`2{md5}6cj%=sJ|iSOM;%4fClS*=e6E(=PO+0Nt!8iqpG6f6F@qNi@5J+zNW3t z7uHk$$Y6)D6?1j1>O4Atp#+Y|4-&(JOpJ2f5dxYaJ2ucq@IT#J#7>0sExq?Aq?*I% z!W+TpUjkHNvQev$Jo`P<0K#fqz2NR7DUbtca@CO2uCTL{K$+GA+Jq?e5+lNsM3+M~N{md!r`$a}2;1@* z1g)&BE}`x5U{)_BxlCQ;Uu&U9@^knl-ly3RN^F!e@Oae1QPD@Msab!f7|JT9#W`P$ZF8!ZvFQ41WNK5BLP2^8?C9@$gwqwNW za5MS>k9EpSUtuy75yE)1`uh65@$vD1)4x@R1IZBhJyZV9PS)Lz$}|^lp`)Us8Bu}V z-71MYR8}>#$8d28Ee!&{d{I~(c3TF34n-r0bh=WWW_f`epIwcoh)y-JFksKqK)Vz< zuj<=Owveoj1<(*sP*4;B&Z}L7N``sTwz8$Y-TdIhNUD+H+_&`mNQyOl{5WV~!D4KD zoDO)CH=G;~w6tzQ+-acDb${*RDJJGO)14iqUgus(EPOZEZa!-;<(^LiXaj^g&G`g> zD_vw4dcV?%^tqXPLa;;15}8Ew1Jd@cE(%n`ufvgFJp-`dB9X*R#ch#QuWa=wKa#1L zp>P=J@SoGu@ui7HmejFrfY?jT`)(Pw5Ff03BElmgDqK$h5}#H;;4&$LjF6NRKGdj4 zeg67EuByXeI3*WBnt0gr7Zb2N*lorePGX@7qUGbG%GWFxfAHYEA;-nb9vg~C&IR=T zrKSw<88PtTu6<*FK#o);<$SLqIjv}73@wK8+=UB)QH=7(P_Ni@&jJME-|SX)807jn zUszDk1a(_vNeHm2Uj+3F>-Q9uv(py1yCh+=P zyZ;K|Z9f5T%0;gPw3P$`JEp8mEX${W9UA-j(+@(Pi9T1F3>R;o<9Gu;NEQ>YFnZzx zSjJ2R=#Hfb;Ki3D76Gyf9Us>xJAeLlRTU*$&C0PP=W%>Qzk#avyPKY%qwA>XmJXu8*ApfQ> zGASuMTE`V1DjXiSt((&K?-!YVXGl~g^*UHRjH^=m1*w##`VYT7Isf7vt78_jh)#@bIWJ_LdmS*OJb6CGrJ8 zVR@0$xbsEL&_zxoJt$N9Yg|f(mEFg^GJC0R=p`m3*uQEw?Qyw8v>&9PE5QE{2_D=+vK~V4L+#NF zY_QEjKYI=24wN@+?lb^5sBSU!kvfl6tW`bj?b~FK*BYXAJ%bQ1Xv`W?)cN1a7RJ5O z^XJba%h~CNjYREy>|Hl*sM09i&s^g9GRbeQe4Pw9H1AG+DineR{Lfs}x26*54wYef z_zn*+6peWJrR;KTh~7|7nTfZXC7>M$V9p!{Eq_#yEMm)>54T6yT&iJ*iHV7wfU5Eq z&@kwB++)X}uCpF}k0YB`Tzn4bRX*3P1=BCwl$y;zI9tPttJ#~&t$-TpE6fd2gm=0^ zrtj2}vYRvt1SFjX+p{0w!a%Pn6Mu(Ef*E`GA z*d4-FhwMALYu8QyOyz-{@=P9ihhRhqeA>#&O4glbxMu?(m3AQAZLk)` zx(n#N_jc`D1$q)Vse98Z&G~NS_`nCLP?8|Er z*gqbLByC7726F;m4#OdPz1RztS_K!GdbX|F;PSip*U#v<&R&x}0t#Nj)HHioIc@KD z*;yhYB7{T)e7mx``eAo{y2#{fDtdl(d%GEU53bC@yLazin>NCOfj~&Z)8scn{(27| zY6D6aM6e%uSoTv;PrP9_{dN_(V1QR~x+yr`O4$vNosmaD`hai-!lq zFpa-Xt!s~|7c6hsy{1Pn;dja`hk^l=j6!(Wfhc*)Wz8^8z35o^kYl_NH4<(rcb2}U zNkF-ryZdWzDhMS68!k~XbkNMRM`fQH2VE1U*QNPm$*K<7*-Ee5W_7eL7Tq1bg)$ILju zv=DD^R7vd=Wmv5v%0cHBU}_+lokwb2NC8ElUHovwVoqA66*(Y~8xdk%x0k{IDb#>GhY72b4vnUyVATT6IYRM-yK-dE-tgud5O+ubYAUF!_s<9( zPt8zKRrQC0yW@OLI)DYx)$V~*1u24iqMi+7(|RKm-WigA%g0zYVC+ zAI3PtBeFwOM&=AskJY+u&uM#Z&OLx7HU@w=i_<2q<6 zAzxJz=RJ^B2!jmb-=~og&>Xg67>)pjj00K}GBLRhnE#2_ul*5_K1380h>=7J)2%n& z5Cr^sg;B%o&gw**vIP*A==v-^C_H>JN3tG)6x;xpJ+FNP`xMBL7zZ7~NS4u-7f;3~ zCP-*#h+azjLFPqJTsSY_=$nr==XxM9aDMy#JwK6vOSzXa7Q6~qlcwl@$NGn8_z9_C zMgGSMmWy+52{gm*bLt@I-$)_lA|>!UcJ}Pq3hS|BP@0|c6INXI*aBsXl=e7814KFI z;o*U6yFI(5Aw$AN3eQDKY>-oJ#y*}f&lC~up48ITM$`!&p4*E-cVOFf1PZ?vg21Bi z;6b`wv=d$GlpJgHr>|eXGSjcYmBn2IhKrGrk%yO811j!R@+aUw_)iM%lfe-Y64utY zP7+cYE~>yQTKrpli9ZfX^+)}8IspE7Mh{B=@c+aZ@;{CiEcXA4r{u4fTEAGEY-^L& zs&OVsPfthGS`c|UUQ)mtL~7r1HOv9NX}&lA6D+5(u`wJkFE3F0(turASzrGIk>48R zv-^|MP}>>Wg)(OV8F&JoW6}Q>2Nu1klG1Mu6ZGkWkza>2I&O8)yN~quyPEln6-!!W8X)`BvUP zPaY9JuiSumJX$VmtgT}wb=M!Lsewjt0$6=ume3~R<#jPRCi=m+vWF;L;9ufOcu1V-v9Vf7)J=S zkcYQ!&=sGDDu52w?bWMS11Bd96(2l!yoCYjpCDa{31xB_N*DMg&O?)haKHCJ-@xDm zYPH`WVs;VReA!I~MGRVIuiS%xYl8J1sxm_rHX;c;Ho^d~LxyYJ;%kzCNC-N0*;#03 zX9o|Z54E?$AefL&!%{$s{JrgEHpo;~-^u`{-gMoZd)#dY#;4Ltva9zJVq&g?pQ)1@ zg%>Ns0M&uu#x=c>hjOwXR7NI>nlG0)20N(#nN`_j0csZ)7r>=rP=6cMy&yiu%E}65 zB-i889at%madE+*zD+=QYl4D!^@|9R=_{+Ds`H5KyuVx@0hmjFsLEkyCN1I_5NA+a zp-{iI{k1jp{BWr`1vod3qd-O>`1tW7VupFt@{IQ6bmlJCd%x4x%l#zJTDwm(Gkdp5 zzyC^N?bs{u`iLoL(nChfP)*VJzLwOe&KS zJn(Lfxa+S;6j}UhaT{zBkG|m&-L4U&a!Cf)%g+8T?Hx9vJwVjp2ZKU0TmJs)1w}8c z=2oqAIzpz|+S$!7gAXV2Wn%u_pUnoo=-!q(HD2Bi^_u0b`{XEhcp0_e*~;A9j!X16 z_srdubKkD4ShL$r$*c~V7Y1$4_l3}S?%V_N4dK#bbH9A&<2X-A$!ye>lutVUH7!8> zok`?}d6pYDNC9t*gO`OX3-jpF>#4}7sI&qL99jZ`cZAWb=Z_KVPst3Ouc*qI(9js{ z&v`J{|AMIsQ(N1!YS8?yXmh@&2%pp*ySJM?5pdSjAla)RJd83O%gv)sEic#GRA!~t z?0@DfhR-?$eRQ}z*k7PMr!E<2iFI32Y_8vvrsA^}qCSoH1;b~VfUdJi<9?M?5biMd z2J5(H=-|Hi5VO$DKRy z!JxY8EA-S{4%Tr>#=(7736qy1R4r&85=!>KTcONUdAZyYwjK~SVm+y1bJG&ya|LNe zwPBa&fpwueaMWFwEZ8k?H}H+aX~s+2SH-(bZ+6=@R8ptt+&C$G954UQ(>Pkytonm)dfOxsC{baCs30&EbW8v`%q#b( zrB|k4#Z%PI$EfSDirN$rKYmtmnpzjC|uBi;_RHH+57|zFU+_QKJ4usb# z%$s#daeb5X^X8Gd#N18`n(f)MSpOqTzy1g)#MrD5$0Qz`W1t;h;B&l*wx98afJ#$h z)bY`->u!M8;k@q5NWCDEOTjkw2zUd;CKjsUkp}D#ZFKv*mQ8I+3JNIQTk~ME+k#L8 zg?$N zOoD-ovT?o;Q(uGe4{ejKlzr4V`bO$``B$YPD;9cI*0Y*r79R{E?w*7{m3)Rf0K9zx z{4n|(`Nj)z1I<-et86Cxf!m<~$@7@EdZleT{L=TiZ?O1rh_$WlIna)a5VOMiR0pZ0 zr6nC7U$xZQ1Q1B`Qj(IlELfF=pK6q4E~Ayfs!Lueps1{M>0P@1NSu7f2|jby77g8z z;lWPw{?7^7844p>et{1KNy8!{D1c!QeytD)$uR^nHJ64bpufD#>?w|(CL(=0E`!ro zHE(2DtfeNhycA>G6B8xOsU$Xl96Riy?7(U(C29oN3Aa+qvb%O`e0)4bDVvy-)S$BA zj3EG3a%FLSb}>V309!HqE*z*CG8MEM!?e}~x+*__3t@5I!_RGwp(8P)Z}jPkxh}T@K5zKx{fyqAEqYiW^ulbJaI^*0_K-g}5lU$d1NYgRsAo`Z>VV1O){lR-@7B>43ue zqu-1lV(0_{cww;gG}wbL0@mt?Wkc|`EcHwx4~9wX7xtA@Oqu@5)^a~>}-iU z06M>xnyX1k!5feaoQ?{~-eXPB(b75tL_`#6@q+^{q&}WC*jO|K)QUnO%n+b^d%1N` zaxd{a6N8%*bdITxZgIrzD=NCxoMFirm6UW==q`r?7-c{}>EQ4aegYvD%%3Cf7tbLd zbSo_cD11ZI`;FP}$R-Z&1k0CVh|yZ%F1O_mlTLn^jKGEGijCkEAcK`PIyIFM%l9A| zC@pDe#J1OR`N9L4^vq0iOMNX)CuthUQ!F0@7imCCjL?^8hn(11yLx(B?g>{&C(BE_r2LnXW;TWUR{BEjAVi7U}j+v z13?dk0@e8=y#46tsCV^J)oJkYa%w(>&lWBFJgLtNKC%(_bwXQP+qb}3gkHW&i2{!K z{dZk>Lbra5mW-H$gu&6lE~uEzpl1XDWrV8`p)_REOl@=@@mIFIlB&A>5lo?A%8r1c z=HbbjFNS6EVQ7f@uIt7v&F5vlzf01E2-%^TiZURE5XU{b2y_Z5e7IVY0aB^PZpaZx zWra}m&sq&<21*!{T<4g;NgH!R5Eb{0UaLHE<4Eb?z5)Rr7@rjsm0% zyD1Y=(2C*7-95NQ=OCg zbAZVkmF4gHftWzY$ao6uV~G39s52o5h=?BvyuVl7%_y&q?M}XBAz)%!8IwhUi9x*w zejg$2hF#}bO0G%5fE%`7)qY2+Fp{ z{Q19j$^S2#Nw*SKxnjs8TlVwR)M3EDBD5E%*Nh-ho)Nk$?a)l^6vJ{4Zj4ToVHT}S z!KNb$Tpl>tMn_W$+Chsx4vVk$&OOVs&WXsB>>HHy#m@B?(L;wq7U^k1!j~@~h(RY0 zdrnMnaImnJ)+I33fi~Kshy9nhl~!D=wEk-|veocgG4hecq|cW-`o3BJwZbN+%UkQG z;srd=V0Xgp)kosM^H?nz?MmAdh+qYrJL1rXS}KN}1xpI^K02y5+#L&{z+MvF;dHp? z;J(?T24)N)$fZx8K3ysFIy|kMA@c-`6&Jt1qehs1Q9~~zt(NgB!*)RoxvQ`#Mj;|2 z4JFNHn6tn8d-e^Y>@v)9VaXvLTW}x&_Mk%WTC*8Kaakn2(P!r6y$Fa;IN&Tb;@m~J z<@#TTY!J>qzeLhUq%7wy!HjhJIkf_<4>9ZpZ(&y;W@ISRr#^=v;cE)+4Px8`I3zwm zQ-o2cJ3Jb=5Zv0>4ew-SH-%Ut4WmDE;3UH1S?K^Q;DP_%ETw15UJCUC=zNS0ESxD;7pZ}f}3I+-U zQAx_m=99|H%8G%XP)5~Kp@r-Ac6HKTOv{3h-C~mr1NSc#Xls0hi~<8 z3Uq3NK~_Y3Z>Z7n@!xs~s2hKa-mSjG!fga#3*`Es^@4$SsBl znI3R2a3S1bxYqCS3jfE1=@0w;Bo>%=B_(A2h;F^l(nk!0fOg`36-LM~M7=tX1pXI5 zd4Y5NrbcnY=yXG4BOJAPS#_Y9`*?AaI<0) z;U)(52(VY6G?VNGq8**G~h0~o^OSYhUVU$GxCkLv~iHb0|S zpt~g{@xT>fWh?xb8wvsgR}TgQzLUlpm-Q9oYQg%bUT#SZekjCm_D{4~*U%r@%om#b z#zB*K50-5osVIg=Ym=?;GM}QLW0*k{SJdK-zi6{8kNo`S=Qh1GOOKgNML$r@!M6Vt zT2$mW_#10x)tkedQloYFLp5K<+-;pfK|9}*F#MOBTb7EDMyw6yF-MB(>#$;uO%A_f z*0vmV7>80`HzZf>io?7Z-hpf@@|b32nALpZQK<(&MZSTG74WV;9J zEgkLOZ+|7@HE5=+U*3EtySKzA<2fGWUq2=vlj4Q2YTe~Xv2M03pS>Klw`EL^QC3o# zJAA8VrF!=i7x2npf|+JdeO*?*;Fh*eZc+xc2$s8eS=nt3->UHR|lvev+GuPI{Rhmj5t4TsAfp+l$qt7Ic)M(yZ5^j=y!UJU=-Z?K$&| zduw-ku0Kn$3ghF`vpQ*rvDkNXuxKG3D#JaMGcn>-7`8Z=9`5HmQ&u~D*`dbq7d1~x z_=)2)Wz=`%s7ERksqV|GHKQHAVfb8%(RGgb+>UcE&{JctvD|dR)Vg(ARP7xPnn4f7 zcE+ZK*KfsfJ1)kdPvR*mz7Fe3E)4HTE)091VuIDF&Ip6QqPdO3m$oiZ_pJWV!dUFu zI^3iaF3%4LKzq);;4X8_M{~OrTVmxByH_{tBmW|V6HcIjv%v?%e~!;lrTnh7(Eh<+ z#F$epE$!S+f%ahk_DBe`(DlTGIg7Ndvfm6=mtbx(&858RVa2rW)h@xlrnA&a%?Ahe zSh<*#@Iu|fFlxSxa2L1vvP2FsYTZ(4s`fUeW-V>Fk%a7UHtn3S^JFI5#+b;+x&C%G zbA!gl!7})Nr=2a=ATqLQKFbu_!VYZ@uZKPJ4G9qX|eWQ3#z;4n0%#zaOUg=OQ&VDE39!4ZyVrI@GYd8Zi8s74@3kHgv7LY8MshD-X~~} z`w}}(#Tp})`EG1oO^w#{i$|cA!r&`9j~)BL=QwF9kThj(J(8=!{N~_@Z|&Ev#NFuG zPCnYg0`7L=qhavT1S+O|iJ7G4yXC%N5+%Rxu{&TK!(kwvhLzemu;V}aBEsKs>VkWm zJVddPE+Q8S))|)I`OouX?5Rhv>GtDzl|cvC+{R3IwzPHjcFu#v#vw1#qtMp&sTVY2 zZZED*Jgcr6{A?-^(&@tOAUO{iJk@HX4p~s(Q13830pClsDRI8s?Eoto7`y(7A=4-u zxB%AMGmmGJ4dd-8KOF~m@ktD07uv^J&J_T+oc? z@1VGNF{nVhdK_Fauu*RyW^tHA128K0S34DE=VG4yv*G5J=P};M6Fmo4_j~c*a*>un zh`MT(mMv%ShYRaY4az$&JX&AS%qbA38?i9ra}V@9+>1^BhPxpz<#5vZY1WGQ&wl>N4=oW~|C#88>Zr2$QpZa+C=bmp*L+k&)q7oE2OhfXNR~ z;UhXd|H^w@#J4gEOcE5N2sX^f>z!I9nCXrItzIAe$YEvBtkJ7;Q-N-H$}jLH3j@PR z5DT+38~zc1U2K83e-4Eg7Sgm07E|8Cr&!7hChjoy8ODK{uoGYiFV4occB_I50rD&9 zi4!M)BkI*K3ZIu5c>bs1UnKpxE54}N640W(fJmSP4*c4)WGY;`w^a*uI)sA*vtn<# z7|vxTdwv9>%LZJLRX1L?#$mpGh-Bh;QAdRfbFB|D{ecuWi6Uu9>vwk?931fXp8@2i z+0B)tE)&Vl6CD)T!rgDk{u8jWLn^YUXmMw#dyq~KXcT~Eo5>PC1p7L6OfqKy_C~xYrxDO^ z_Eo097yfZ92a6bi>KV7-vlU#L6D0PCu(d!$3*QeKygW$D7LUKN~_=S5a12wN$;~L0K_|e z^+$d{-q5qL5usKFjFZQJ;xLj5hgY=pLJ4#YIk1hI63UTF8eob z0zL`M$iZM?BL;ism)&b4&uB#u;-%y&RCJ<-^Y8?@B>o_v_^ahW4X5ra)Qx^lij4UX zLN0*&2nZZ)Mo+MeA#JDsQ0O2Hb>Nb=SsA?tG9a{m;6N|g3BZ66P`Si#O2iZYsWYK# zxCJ@5fh3YJ(%#e;%q=Gn0!TySBI-T#KY}p_X^|Q@_l~RP_Az+YVG)$Glk+k!U%f)y z@jX2~K#+=qd7$HEyzK8v?oT1QcK(Yw9QuHCZ%noLA!D*6$(#WNun&Zz6i^A^mFRW> z@G46M!Tpd7-~wF zz@GZyC*%T#YlicGBL|}TV$BK~5#A}19SjhN=1)O^w0Vh&iS=ORfPxbIwHq-xi4<-Y ziBPyf^9YzuxqivRrFMV2_&}ydptSFcX&pem@0@T~bj9@Ix|CzU|Cw65Kw^b<76V&k z<%(`K%;=l7KEDQ9VcPIx#Qg+Cwt|YvRbF=Fb#8k0?3=LI9)V*OF?`vN1Ji?SU?@u4 zboby2ffm_FCt;9J3WeJ?XnVl~&F?g~z^*_>L1D1aUxe_WAe=*ek5;EsL*6sq;G07R zLa%`}a2(2)kKlq0LENy=#?)VKMFVs3IguW@`3{cV)I{*JA#P!$ClE~7G~fdefsok| z8!tPyXb1x!Nt0&9$O$I4#CoTNejtmVq5V{#Ut?pzg4l{t1lF~e#hfDts7z=taPUqT z-A6o7EyEYIWqO=6j>83|zG|7dPyd^MA^*)={D()V{J%M6I27F7`N8+WZt~^IUmTk# zgD-Kf+{Yxa?tb|62^ab`kOLF`9It|cV+HIC|F;ud{@*?x8Tr2-{__9frSHEE;v@TZ z5%~WnTC7NC=WTX?j4xi0fi>?_bMqs$v(*&6{K>QLq*G6`R)qoDsFPEe5zw47(Q^1KSp`NN)|8UQl27Ur49J zFSqy09#Vp~{1${fc6LRGjESBnu7EQI81e;O(#D2{?pbAz2q00+Gjl#crsGzVftp?( z0Rc#_FAAyb3coW52Z1?1naIZkc$i-V zDCH2&TfS8YUt)FLTuunmJg*Gj|8m=4EIU=eV-NHsVs@=<2l_qOWqoSwrJBw?h$O(d zAO}^zbeRvHxo`^wQ&g~ySmNpqw=|ERJV|%!mV#NT(BfFf?_+CfEUTZN0~-nzm{pqF z+GMOWALfMj32f^=NDm*XF?Gep-|1P+VaCX$LtSl}&F}&y2av@4nG)Q=vC%hd_1)T5 z8nLdxEe)zWl$Lp!-C#Qfr%Y72kb+Ou!9M}VM1P%*{xvhR((;W(;YHWYMR|ik%HrjY zRe`VNg77_iAi4bA18`8U?M_oIVO`c!a$GQ_TeY38;7>u_vK)zD5y*D(qWs<)<1r$s z`BM6^`bh6Wl#y#F3MtKfJr%EsxXu;*--qKM?O~X!tMg zXYh-w2M@G!az0Ip%m(3u9kiZAbuiO6q>6={c{mNn$I$cD;>hv#;R3b80q=^s1l5|w!@gmtF6D*y9ld(%QH|y{1l!3I+@96t`o1Q|MWtq8 zVSDVbT?G25M{ThG3GSS_kT>y1Z>s?B6$k4thmJfh#je!meP#xka@hu#FIVx3isIn} zPG%-^FT|g_(qLX_n(Vb0=UVFxr1Rj;!R}oD!S0~Cm|xjqsrlfJd11DJ-P8NNUAV*`MXz@GvC3%IXy+KQ4Uzgxi{e;tQk<<$${>m*)l>rg=V5lUz+On zq^YZ?FAtUdLyc{SyLaAOvcN_6)T#Lbr|#D=a-T*B}u+TgEkdA@bQXJ%tg5t5;3>*p0(o zsrkikJO4JJtW01QyH+L*b1&0B>~2vJ0r%&#SFilRdoteuKrRgFAy09*2}%wa`jFFa z5RU*bcHrEr`mhDLCN2PW-*m zadFUiSFJ$(cEr?m!I8JaT<@q&cXxJTN~ul#HkSyYO+Qv9!>JOOxI*D4Zy7qx&CM`SFRfLvtYn^46WI`vO-qDL4~qh zieFQ385!f_d_By0!FrOz+cY#-pYX%#GGbUB-d+Sxh0$X7#<*$SuhrjY(=>w!>M#z4 z5JED(EH0L3ewnk1k1b`ai@SdFX5Q+Up02L1pA;4GtM;GJpcZD7SZ;}nLMoy8tr+h= zpArYBT0lAN4-TUFhL_UU;OC}!87aVT0W6wjSD;pdIE-j3wzeEl~_Z$wNAVAS^bN9aa+50V3ZXIz{f)gPGYAPI|JmB(?U%+Gnjb$Kk{?FDBbJomV znYi|=_mn-=SKgZM_e&QuKHXaMl#NHO!gl+URAm9|Valnx=mW%c>8A{os zV5Wd1T)UnDMokc>;QRp|U7x@}u!H>1_K@G^XTVe7ZiVKY$J>o$DYYEDsKlo^M&maa zP0oRBUU6np(mu1b#+8fBqCKV4&Vyxp^}* z6Rh4@8c4;FrqL(vD8(9W00*}r=T!x^unzQy@%%WdFJ6a+H7FDjU9(gf(KV?) z;?Ibih5wezNCvL|u7LmL;J&|Ox5rohv~q8Z zrIUm<{06WCG=Cmy5l4>T`Og7gn^27jBV5Pqrw~WVAc{pWdlU=)c_=dM-wS5{`{{yz zzep&-e`{^HtHbW_;~S#o4&$rC@xL7wE59ESFz z1^Y6oy{`0!9igeY3{5YE@d#y~zsHSQjPeHfQfaHeSp#a+YcSxP{Aqa|jzw&R!*U4x zK2GMSvPM%K5WOiXBV;MGimj`7%?{FW$gc_R46ir*(H`cN)5TUS+gc)nl-~n38wbjp zN*|3#TVY>#%#Z3F&5kE{TRw&3(70ResLZin*7jb_hEP$O6a3!855w~}Fd~@#@^9E; zD_J+o`!kwJQ)%MB82?p$kpIf){`Ln&{@L{p!>{eXqF3&|4$!k^YtS4zJIfrIJI?-r z*`I&~h8tBcwgQ6h_dJ6M3%>rZX9XClhKGkklPnEr`vG7~c*vJ%lnIOuz^52+Xp}0O zf|IA>$vvNwh*t@m`$XIp{d$Wp;=-Wi3ON=LHr)&yJ%u#lnP+zA$t3WQKnoHTRG}7& zDhtCPb4Ng<4r1k^@z{v`+BKa{(vkKQ40&da>>vO`vpKaw8ftTGhMG-zJ??monm9d z4vnCpFrONH)YX8%{XkT|%E@kwz_OMKxVP11V+5D3>e{25IXwG*Fkoa0Y9! z1^}^(N&)TwuHqov^fEZMkPky4XHr6AS_1^H{CE{k4ycIWV8^CG)Bak-$p!4=4D>F- zaXmiZAA1CKv3R~uqsr-GC=veKkzJV`h+7?aWH>D71n$W*Anpe%Vy`nZfAq$|Mb0sR z_}wKFi>=F;!GgoV01i1YYU}sBEjFS9o+jvkhPFv{P@`cKJ6_ac%!BLlSK%Ggpyd=) z?sFw&M6v_rEQ#H9bxho7!Frt z?RsR+SbXYbF}b=1H9yFb4B zZQsxD`|$p}-!E0A|B90%2IUntJUaT_{tBLW!!JmvA!QeGracPrB*MAH$G3{lPLDwu zT!C7 zW@II{Bk~*GQ2}nF-vhxL7E?;3U*y>ooA7Lr1^k20~0!`Z6 z$LGbYhJ~IWWthr{N}?t_VsAq$W!Kk#azA2CkJj)6%r*c#YD-J=aP>`vL5+Yc?~EOe z7WbqL!~FU#Gw`Mzt$M655tJJ(u)jcgns{R z^8g`9$(Oe;9NOD&>45)l1e`>0KI^~y?jD4uDqx=`+2ID=j1DQ))j?Z7_C)ER_Ao)) z$gc2*!mbZi+7Wtr=$v|Wb(GW-M75^Ws~$Bc%ze@Vh4s^l_q-){E>8r%rly$U**EJG zv_E&W!9r{nu8CP2=Qsk7E}&7D_#$)87lvm*V8I9xC8Cyf^Sad8$R87%dmIkH6;1H_ z9=zu`t}kXTHa@u?hU3q|UO?+C-g4$JM23G>ze}0t*!cOF&o_E7Xn)gvnJB3F2v!FI z0NNBMGIN%#mLCkt>I{=Iga3v9F?WpUtS5CU2%lEVmKP9d3Aj`gED!HRF{Em7)S(AP z_k!zMp~n>R>c-O!``Xbq_$<(?fzZ6uEfu(dJFb9?eJs68W%vHi6Ex)-CS6Li=6PG0 z%f<32d+ajUaAF~p<3KlUmStJg8YG5n_8bsxWZf6_YBjT5oSrhcDqo^@H6Q3B7XwsQ z3A!uTIQZ}yFqK$Q0Oa9>3U#T!%k#glfQ7Hxe;gJ9(zD?PRFeY{FT*k5Nf9~9rS1iz z)$4uGf3qH+Dt~!^_mPbOv{mIaxQjYi_bRihB0a-QFl?8O~@)-@7& z$lWjlHSkb027D(YDpjX4G;1Scpw0IRCbt2qs2~r+XRax4oy{l0D1jM8zHVIgVLpyB}xf=Rpw;# z%{f4NrjwF)?>>nQbOMAXOMZG<8qF$^K>&-rJG3`toV@gKD*C0x3Kp_){)@?aqbO(x zpU+$exasH;;-e`vmmp|$g&qMgXCK^aR22lzLcvs_@k#)~V+F+;g**1>%S{->s2WMP zS^^`LeuSnJz8Au(iTVpQvyYvCaV_EofO2;hCyHn2mssGx3_Ac7T=mtG=Hj>YypTQ5 z@pws79_Iyud15Tk0Te`P)aEi^5`Fq8QmXAC=ZCV~NY8@KE`Y_EU9Yo(hMlR|e*=_d BJ;49~ literal 0 HcmV?d00001 diff --git a/Docs/sphinx/figs/sodShock.png b/Docs/sphinx/figs/sodShock.png new file mode 100644 index 0000000000000000000000000000000000000000..1c2850a83bfb8ab71ebe96ff73186516ec84d97a GIT binary patch literal 37811 zcmeFZc{J7S+djNiDk|;0{Ft@n?&wOY#F`?EivYdFv2JdWeMu4x=rS;5TB zOrcO#&<^a^qEHqx;qT$4i}4%7vOP!eL&{0%h?BP62`3lRGZqvzQz!e=c21`)&DNi_ zIOAw(XS-8GN@R!d`cqC$_Kwn`qBj3}gNWUklcF0AFKoa?nCuT2I8rFArsOX}fq{>ZU0^!*Ra0kbFS>Vb>%-y;xk^+uL_&emTKP6&5G<>pRitX>D2`ySDBe@h1V!; z6zI7s>38}y$DW&eZ|a}dRAJX(*vrMfwBqVUfejm|PiPmGG12O$R1O`qh^}|q^2;kc zo-esE+xCH~Vsdro7a@OBE??XVrB$7wkcIri$9n^f{QPqto%|@iL}kDaA12xY^50eT zMdZh3Dn^7tIdp)H{J49AN`7qrzkJZAjJ?Jelv6{PRTj$Ji%FgjH2lpJbo=+Sx61 z@_K)Yr1#9Ie&AY(!_K{Rx~EPh7-rk1xO_QR9_1-v^O@>js`afG7f(}Z9|^8hm@^u0 zc8RzzZu#(a$ks>31#UepMMY=(KX(jq4NSkK1oPrySblXy2h}8MZ?QOWq9^u_@yE&0 z-d*lvr=16y47Td0jySurG-la^hI@_~L>SfFRws>sBn+T5=!K0*Y4Om?5kEq;~53F52cy#C|+TyJt$i6t)xQe zc8mAiu@E8MpUzx2o}TSDwE6USuU*khmTl3@S(;W5ubk7jij{J+6YbnPj$cjI>s2bZ ztVy^j?bOxUZ!c#y(3G2``{;IJf7x0)of6*_qc4_;CY-Dc=N}kMj95b%jHXkHH^d)) zpcTly>k!kvLLZ+OMXFx2zfLxLP4zRgvg%_32c8YeTO2>m-jnUzt9Q($**VEDGf7A{ zDdzA4>F($td7+xYsllRpG8R{_F1>u|QUWGFZ))IKNJxn8>9@H50+FKWGnsZR#^phL z^*bwNqs#fcjp%q{^N7#smere^9V_LmK1A+~68jm?v6`QsUn#rkjQ@)lFW7l`Q=bjw z9o%v_Z2PE`>rmuCPCH9&lJ3C`8#Y*dewz6G{ez>Y-zxFAj}J8Gjn*C2O)=bc%%ODk z;Xu!?FPY?qiWmC%`Dqj;7%j(SgqH3J3o)T^}|yyAKvv zVbI7kxsCQDniP3`&bEsvDL=dUuuHN*hO%ywPL)k|^9ie>@h|SqeV@eoKc$zX{5a+? zSTK&kuu?s6Kw#tx)}ft&Fh|PK?Dri%mTg{2*NIhKz011pkf5L-XLAdm_X&!hGCzZl zBH0n8{+!lt?~3IOzn7`avN_7g$hf-E8$&>5Q8|Ds%5!Q=P+s1cJ~G*x);M+gOTK|b zll?`*tka472X6~dsr!e&ztb-Ant5tg%DiIDj{WqfE`!F--5-^9S=D^NaKvHZeWV+I zn1`m7apK)mr*5ry@68LBu++3Zzqre7tz=o@SHf6mObZI%vPvr~DCdN|oB>!jar%WirYrFg|^q^oKYYlQl$E&a-hwQ|GJ)hiu>*Cc z#`;6<-c2YATx&A&>HSMVny;yVghYLmW8rq%Y|Gqie4JV+hicf-tkag>McVO)IeCWL zuJF};h*H!|)Z)kdX{yVO2Io0-f7Os>4CIl~n;IXiSD5|5N*2u6SjNf9o+=f-TC$+O z6u3Y0c(3%TqT<6q^D|LHr5m#2I~Qs&FQ?cSzq++0P~=!17kxQRSwrJS6=I6{=ckrs z4H+l6ZA_bV%Z zo*L_GsLmPrcrdWNJcx6x$*hBVn55}+ugTnw3l}d2ikX*j<>uzbd5IQHKiTj~(qrOr zqUXd=o3fLWQ+hi`?QXAWSM8~zb_Es(cH9;YRzOCFO zE?y()Iq8(%>DcBF?lsxHZ%|75ECNe_Lwi|>>pGRnnLvpz0wWn)=Qjq=&0Dw3kDolb z7J)%%WLRY2lagTKAcnl(J56*f&wdW;YfXXsxQO_dbo27%w6bNQ`5I=2p1BO}F1yZS zja^eV@zq~6px(s0`4o1lX{KSIy8N$i4?A_jrYpN&6q%Tp1&ARworsB#f9CvNDHl5- zu()`SdFT5FoPqC}Bi8tQR3y74Nk1)h+tI93`2A#GM&(n-&cimJpVEH(_+gl1A5}Oz zXd#K&+{+1^e2Or3|97doap$kDwf~-|m`H4uZcC*2gnTESgf7&khdaC0qm4bGj`mnxwEZ4Df z9q-hA>(r``u?KYv+?Bz(u%Epc`%yU6nS1;)wtAO)va{C6bflo9HVrym(X>{chac_FK8B+{DYw~Cu% z@~Pp&mkAD#o0}fKkA*?D(Ur2FWtGc)T!4V=4Y2F6&mRK6^mxxrMQBfbH<^x&ja}_F z_$l3-9>cxMs_V%ckz2J^wN=q7Jlx~A)a9#gY*4JpwAS|=YdGn5{rcHcMISzV*vPr} zj?;#S@bFZmpN4G5GiS&kk6&Cplk^_pngI39HEnu9>rg_cz5ooBHt#jWdv-`(GJ z@YSnVldTt7#T~xxCkP|LT_xMy;F_rWv^`qbFt-;vwc{kph9PjV3IkckW#<$7r9xQLoC2U-LsV{C20?_j@NtzTHxnUKg6}b%G)rM8&TC`0=CPbe)1X z^8uPRo1a$h84HsCAL*ynjdWM{x7_`FdiOSyLZh2EZ;}jPk$$A6z}=+)pn(jurK+~( z!|3Qti#Hr(w~XzN(u`If8(*d;{&Ux<4|@oT0N6VnLBILn!LhFo{`}(ByRDkte{QO1 zPW$lT>!}4}hGbvLxs7O;nLQ*)v$gflcrST=2j@PLt8S|BR-KISj6*p6_+{XAZs+fq z=WmvNgn8TkIRM1q`}dQ4SFDxxJgqXBCp~%ss}QmJre2Dng~@SInqA?qxSqO{M`Qhs zGUDGad;RL*oat1U)9kHFS@G6HF1s(XGlSl9)tJ7%7LWIL#T1X_IcXw6?!t0Fq|;ECo6eA*XkEC? zC`Uv0k>0LAWLX{?zzKv4y=>d&r=}%~rzC>)k2!}>Es=$|8S1Eg6QK3OFYdMcg0`1mjf+R@<_x$;4OJ2wAqkZPNX9j-x zO`h9rRZ|hw&e37#Jr#Oeh^NX_k8)Iw zkoR;dSs~H&dO<$&q#dsW z!Y}~@0G>&M0>CS0R!oAZ?S?D5&~>QjzGG^)OI zR1c)YDC56lWM+g>EJME1``R3c$F~^it*=bbiuKdXw9;)KYqXOfMSylo$raYMq^8nL zGnU&VAW)HDkP%0KPW!LUh@!5_{G5o`X!8gem-k$9ZcnR_d^ZUTzs>nZKNtY0uvREv1G9mZ@Hpzo67d~kIA=~m-?-~JlWBrcItJo;8>Prb=(@j?%gj+ zN>)6owW>)rsOCau{DBm+TH2Rt-t{4hrF@t7tjE;JHX4AZ=#Rq(e)&nacwC=-OPRP` zx_vfmdMSl6{T}OUG`na<4L6Xfzh7iC5V+&_x0)wTY;JkK(4I-%XfdTY$fhYr6RFx+ za_@~8V2vb{vDBp$BS3qK&zU$xaI%L^yuI7jt3e6Y{Ht)jYfI|r|D)m>pZ|YWmH+=A z{NGR(|E~`Y(}*)=w&FxaH!uN-cVIx$N(jXw`(pd@2XgLugxT1$=f~&*`lXHQWX>Ma z)YL3wS%^0b`2`2FD=8~;Ov%&stxVEO5fU7{!pXt$z5NtkrBqhDeXFCR6FlRVEh(?g zhO8cc$HdFP($doEOMWiGR&46AXB-U@r~-jbS~9KcMR;UgHGsRT0MD66k7inhguC?! z&P=oi<$d+rD(>8)vtpg(hqk*vCAJdsVgpMu>YM2w7c&%|O+$0KSBiiv3= z%V@*ZIT!kE=ILr$v@!UpIM)Tk-03yr0BaGJ1|au zcEa4f5j9GJRc(?b0su>C$h~`sLqql;CnnSLy?Hm>}$^bu+bm;QN}gPeWxgF zEw)XNl*a?>i356iw@Dp>C>r{?JYbEvQ*T}W!zq`+FU?#sXHQ3sB_s%7bk@XWt=V~k zYqYoCMKUWOV7Ve4iMee^KAbM+U3k+voy+(5H$&mFFIhGRn`SGwWpmBtc5uZ1>-!of zHz)?|Jo!F&vSkatdYP2lgWe3)7pRhS5JWwHb=*POd#hmPGNfuu>PK?xc_IC;9q97H$C&ke2dCDL_0YI1ZP_2>ax=cZcD5H_g5VB zOPH3R8SR^__vmAFwOapQ71o>fy3}Kh-D2(S?IAaBntI2^#N>E*WQ$u-Lxwc`H!?1g zEDDk+8!EdyyTf-$*;AUEiq^0SU&u%Sq@uqJIh((1ie&*YOUA9d}Q0wF%+hS?&qN+*}rH}<{ za}OJBL94bzAa`t&$T6p`CMmTr$*TG`4D*4FG>XS-v-DGS$p)a7-Fj29M?iDQN+AAD ztiz|w)6>(rrOw0*++vqG^+%-|J|nL}r%ar;WK~FB)ScG&@MjabTSa_3UxL@9`21VE zE>anWmIQLYS}HJ9d9!mfKE($ve0uNk7-1%+vHcrVw6r{Ae_w4cuI9gT({ZUY`zuXc z9&jJXo2gc}$Y}83SxuoN_@MCiW?Pf!-r}KsQb@lJf!ZYHzFg?hhYKh^TI5dd0Z?>z z$4At8(rvH#x8F`cOs+8iUi9F(RP+iNz|}PeBo8~%La)6BgT>s)tz>;8)57+fwL<*- z>4)*8Id7|ajfC(3%<^Y0ktD;yRKW(q& zwHHp^wZgNjGGd(4lml+?3)1HD2IY6%5J$P`1kA~|gpooKT7VBmC@EVF<)2IXn3GO$ zqZ*&Y_4)5f+l22i0{=>q(-c*1>U-$?x4H7r7HBLU$^WUkjXs;#e)HYtKyD=^+^Q?S zOoL;c?2!q+t=ao=d1I&&!HmCeaU*?Eq_Du(+X?<==D!yOGk%jiF?>5=35W294i0~1 z&(V;qm8ZamUVLI;W`4dF?C%XYow3@6EV;+j`eet0A!<7$S5< zJ$}619%dv)cLby5#kgJ4)pxAmIy zM6upfuSqbY4?g#1<0j`uko66DnTI3fQ+~9UcQ&(xG?K}Z$O)sO?otLTqUYg}kC(U6 zJ{P!WYiVoCajf(4U$%UC4Z<1wAAdxC$Zol{+Ubdu|0hDcZ^>}mqBHHcgXWWgaKPG7Eu>1I3_1j31$PK2xEG_L0 zcOR&FesKv3TPwT4>iD3-iLc9o7?%_~tXjfYz)_K z|I+Q$sW%*d{5d^TceE_XpL6Y0eVT+;__?yjjEg?C$T$zbi|Y?Q8J>G8!r#W-Z28x( zP1ABB-g4RI?d6Tx{i2U*w_^v+tY!9zut%Vc-6Ep0e?KV&HA8pX=%B>uZ+@l^u#{Vu zW|D?7M1Y^abhzP=;MZ$!1Tt<5YYxnezp>(aeRD*V=5p|f>g$l>H*-BsO*h;Y{he$@ zvih%yUG;Yy9?zK`Z9228&o0{`+NnQNsbFE`g&#eu?tp?z>kM~4+*uiMhl-rUzh#T{ z*O%A0#ZTR^Ub)b|eLUy2@M{`8aM z@22V_Hg#7w+P3fYr-Xg~XcvtgVbh$4`a0#Kzi45C3h%l7;2dkQmvzDC#FyQdunRXc zGh^rGPWtreh)u(j*q5`A81%p^?o(8B2GX$x*>-vIEm$f@cy*{d9|B?p@+laDOK=nt ziF|8alQ6`1rPQp*^+v7G3O?H}hYwV^htzbMHHQCVF}3aMPA`v%p{f%E3Z%0zi7&icIxsr{|EFU^jSzF1KK%mbx&OJ3^K)s#zM#4tQ zx#_MBUmliDm-_p2c5eH#Dd$JpTxYO~|~FW|x+RM&%&h6wowQpmXRk zWb|7+#~;W~_3d76i5-n*oyKO7QtIpg1$E|&NNkwTSr-4@z}KFsmHG~gv561 zOE34}VWB9r)n0xB+G)3(Wy9)*h6ci7X!vukW4z}| z=Q|?Do}y;(9Ie%-$6)MeWb)M{n!lY%-gA#*%{>?WRcwnQ_klALv0h9mXUuy-=piQv zY7m(8tE#Fl*KO$gQZV(&YVj~L3$rV(10!8}wMcld-jk`($=gKYZ zlt7j=Bo%}T!VrVWsd;&Am46`r)~yE(3=GEl-@bk8U&c!DQFx=yW^{Kug_6q8v~cnA zp(7N=qVDa>qhH!7`1N?3PC)=7)eNY8_5Q}}dN9$1LaR;GR)O-PbJP5A=x()2`Ke@? zc;p;P@LCamCU6xFQ^iovZq!l&`|B zUAk{>Z@U%jvL(9S*j>%R!69#Mb}AD!+BSnH`?u;PKMrb0Gf{X4@({vGy1E1r37m%C zNxZr5EozDJp_HBxq#emP5okFDWkB`Np&HaJ-8D_K8x(!JAf>wjtp?w|9mlcF@D!c^ zwDxT(XmCL(sfI5Kn>QvA}qOCZ&sFH?3a3-A87uWU`1w${E`2_`yJ4;fg zFX>0z(g*`3OH|&`Z#Mu*8D+OH_{hlKrcOh#cf&$n{c2$-z@${%S_bGO@Ir33aNZ*tSQ&AHqY6**`b_=5MHwXzVHxg?yWm_YSW+EwaAy2mL{SS1|=E6 z&dJ$XMFA)nPJqC06EKOTgt*x zH&`?$LYo|^R!4y=IUtP3!wWFP zN5&|Lnh(v08JGg9csv3u@lnuYkYtWXd~uCdSqL;l(a6L5YSK(h=6+RrGiO$YNmcjz zV9cg|-*w0(qVBd7dXDvD+Dx!RKOpxe_Vt|xZlF+%6z4~ZpD9{7@FT0j^m-b3baAIH z4Z_+F7EYyZK6JO03Zo_F56_kkRN*`zzS3m2jUOXEu55W7KL))&o=iV_way z`}ezw?*XKf6`J2(^zBerh$bB*`|~Pj-x%yzRr(kAK{;fESJ4a@uTH>A%ux0olK8p@ zlz)MkwlGN1%f5&gDXlUhx@^?$xA%4$=DVgNDmc$;QAkS8(-RJV4Flv>)Vk6>w(Jkbs#+7$nt{Rvm|_YF z)Fb34A-B%2FniXJMR*s`AG`wyXbJoJM|*0u-Q4mT_hLxLw>K}v*Q&B1M}t5yNSJ{o z0P1Jl{Hv599t;z1vp2Uhyc*Nb-{E=mCZfr;Ra=NfgLlOr6pCb!>H>;OPB8UQg#0mh zGc=IbMj#O+fH$@J^f(H;Bne8No3t_pg32MX5aZ+`!Mto)bzREQQH;Me@XC{h(^idH zDXB%XZtCS&6-|3E99zQts+>Y(C*5fd(usNgI_6-h|7GZj*9VDf#R+O0)I6RPYP^7p3H<@lg^ARB^an^5peqAhsi5Nb1UN>5X6 z#u-OPi%i1;w)L1=p7}XTqTis?AA+{~oSl~Yotts-V#ITEUe@I(;;wfqBk=3cZ$M`S zRiu6N&p#tERI=Rc_&95IAl83kx3aHk7x20!fOSDX!VfoS`R?5b(A~z2E@5CS$Ub?I zNHv9f;PPjzIgubnL^IioWI`oFsQEw2)_(EBZ(JGXZ^%jwBQLP_NY3(N9>Qb?D;|y zP(OP%i)druOu7k00t?#HQ&k%g2-BV(UWS6#e17pbyKGVq6*KoG4x5{sR8;ANB?Zp2 zoSps!`q40Ib{U7T#?0g>l94b?3+xLK(bcfZFD_m|rxwgk_DXwBJ_TNi18$^g!GF;m zh7IB6i}BN+NT%O-@Bq(rHc>y#IBpDq(-K}Xr2+7vWBX>XIZ~h(^dN`QW8g^G!y|LH zrUZp4TvM5+s`f#oy?`b08}8h@cQ1~ULa|tbX&4wXpxgiaM(g2~+nMKf(#uN(el!L>oHY4lo6&%%x{aSo6#$b%d}u&6;@-2q?-j`7a(4O-n61r zq?NliUWfkSaIwLD8)%>WZI- z9lVm4Kj*m&rXthckvvg;oo5uzr}|ZG3>1wxl9?`74id>pHi(1`kWdJ0^TR?MoB~%! zo+i38rQ?rfcNlWOXMu+Z007wf`}O&dI=|t~NyNA?|N7Nk$c?`*_0nk*Wnev^SmZil zAzb}CHjuaezSKY7s{HTQ|5NPveX0E<3;zDre}@@!U(*91G}lUepAnQQ>VGXiE4|pk z!FHD8PSFmH16YL{V9BZ&auz^fJpTRdJ&1g!Q)9j>|M-(j>q|t*6`AY;)M}4pQA3A? z7L{hYre}LGvcCZj=5MNTMOXmPtF5i=HPq-MNA40irq~j;uevwwW!+b)y^Q_7RR2|t zx!$%*B-4Y6uK&IW|07*hRr~?iWO_WAG;y>>{xi&dXe|#=_4}3K_u+DnhK7dBS~>@ml^!>o z#@4(`^eO+?4;er5MKHPuGpaON+koooTlcwNavfXK{&w$hd z3=E}dvXOL+OFljnrOD?Px54x|QL4H3N5lM^S4Fs!&pQ12nIW2Z;YQG9B=>*1kSC&0 z8L1%y!w*`>!WbJ+S@`v9M%vwprmZ6KQ+|(g)o9|9k z6Rij&6?h{slDKWxOQu{z1}ByiAD^$ij8rBiNenq6hX)UNyw=t5nEC|-2#!N6t z*ub^Jv%;zG8bDBK8agoe+{$4TqC|QeeX`%JmEqNiTMHJ%T)LD$wye>2#Vw%h@@zRV zUx_kdtAI<3EMH#zIG4LnpjsIGUm?qj7A!U)Sy`hsYu2Dpqr6^tDR8d zJ%T()JkN`lb4AW;Yn5OF;PGNEC`8$-<;P8FT?Wkh0hONf#E_{4>Z`7jE2|CG@4fL; zwTrQlSzfD)&`%|6%{ioe-)-iRMy&KTQJ44({|ExV=`k=OuV1~YLSn2zq1E3~WCD&! z#GwV>cK1m%Dn&OLBZ?0zDMUVc(m9qasxZrV_4vroLYtz+R*9n(uGnPwkjfC`O>D+GxPBZAbPBzJaE!Ekt}{EUUQrbhhs3J*!2< zmKIiOFTjY@Bxr5|IiikCCGdodQ4Ks{#9jq%P8+s3@pC_q!)g@$$syy(li-jLi@!Jx zcBG@o<+Otbt9hJzt;{v;V7aC3vQ}NOcGmg1$H=KK;$hRfPyU|uX~zIgmFqQ7$t>S0 z%q23d+Uz=2&KenD@1(Zx9*u2zi3Vu6QB3zKFchU@GO%P4cDb~yObC?E++=;x(?6dt zAa+jx_^oFmh<6IS&FQGKK%fFLFqt_yH|I$|Trm7IMNd#L@r7ro7d+Wk$c}D5N;yYg zu2P?nUPtkIQAK^8R!=qRBF-;~7LRCBO@2y)^VO;;C*AWI5d~?qDp=F})!6<;;iOT4 za#>2Y2Mbv~Eq$oOz@TMoYkc9cz3JtHe^3S2ueVcnm^EBMu~)u#fQ{F4MALKV`4Uov zv3I`^on+yWKK1KoM`%;`fhMqN)ndxTnxexP zi4SnnCjoYz0_X&qcCZY4`EoxPtMZa62Ff&Tw)AnS*{g%gg-)HhQTem;!F7_sFPm&( z0xWq<%HaYmZelM0t;SD7Zp(uShuWdG`Y&mdO2d?D0matjpL4W0dD0s9MvybsSVzV% z^7zwvk_xx-y)r_SmIn?%t=}&Hhc8hgY3G*J(r?}BAH)ZH)603at*nwna=)>`erCky zIdXGytpBK?=|YM}u?~#_z-s%Wm-_ssodTR&R-ePSn=Z(9TS(DyMJR4%De}6@O`*Jh zNWZaVuCo|0B9*|1;AKY3L+9^czO>jts&E1i1b{J56WsC_LLjk*crPihR|_$n`}hAC z3&kV|5c#V}`5TC6h&;A$2p|xRd`5%i?{CVYw%U1mh9|TY4gcO2H1JNV~hK6Wx_#tNBfF*Ay5}(9V2lFox^x?+TzV@`gu*9#e z(aWX0rgPvH!vab71OxM@xaG%}HE`kyM57L1;oGulQ^}3v_>;44D@doAp%PYK>Lyv) zruwLf%fvB354beezhcw+;_Wf{3RhTNY~6o+&DTHoo5=UK{$XY}!z^$+x;~eqL_p@& zYE~ZU7&zm=hI>N$R=Vk#=iD2&+vf8r*bv9E$X0-my&Z!8POe`!c8nSop z+_`%7>R4z-phVBj?x#?A$5BNkU^Am(r5lyfx8JkZXLVWY9;^e;&W|8r&Ke>EKj3b;B zG&C(xcOrR}#D|Rlqyf1rJ6mp|_2TD@lR;T_Edw7Mcpf}rqzv&bVOggEUl1V4t_Sj- zpX=Ku?nJ8$$BB!dzU#90Yv0*(7%XUp{*Qu*Po|#sLXlN21}#-bL1q?mOW2;aW5uLG zIoc>_m!yHXb9>v-&4Pm*38cV+AH`{IX51PA2_ZW4l8;1TE~K1{NOo{Xx_94Y^Ente zN40hM@|!EHd}rMQKxu8=v`G<6cnS(ebMF+GG{0;?#RebsSpVI-T&Umm$P>cZF&gAO ztq9t}L2@sFQ9{ybx@_>jkaY^tunRtV(vl{@(!Ls-GqUd&>PL7FcOHYJhqjntQdbmx z8&-V*4^(ZO+S(adb@f4zgOd{=s>&CnCj4RjVMvD`Wf>vh6B85N+5;tmVFkw`OgK0E z{XXal$RceiiVG=1>ut|7Q^MOvk47$Dyx0mMW3p&&N?~R=Vy*@RZ9+M(`yt|}CcRg; zZ-4IxGXb*hc`xk+YUR-nimZvyS;z%>-TTCd?Q0TF*AQ5p!S}h#ETZgTBrv&&vn%EN z$+vf-!1pSW`NIayJl(Kkw7*fG!CU^97)se?oFY1y*8J5AD3J?(w9B|H+;guzlWGbc zB@dOh6KRYBW>FeLGgaYjVS}$J^Y@ZdY)yLsDiuV(r+dfx2vaVkT}|*Ut5NYob93KH zZsG+F{Bs8+(5Q3)@mjn(IjenIdIhp?`H)R{6+Rom2Db))bJmY`wKsfoy3ZmafX3nac~ys!DlLe>}3a9pLB`w z-O%tprqIH!2+XFqee0ep<v+&6Db-TuG?r_upsTW`;tsVJZ#a} zv3G~t%K~d*!k;>e#S6ymAs&P?lzQm5aCtiMc^EB@X9nN(sP$m@IZCVyAHHuzMu{yi{jVV&pW z4hSkQ!o5w1;^K1?Slmg(v_py0|0_h( zPzK$G-60Zg0Ydh@xg%PSs_dy%tx$F)HbT+gchN??^dxwJ1trx7+6v>wAwu_mkd$#5 z&_T06U*oUp3%;y;sbpXM6{yEEGBQR#nFk3%IJ{$Y<{BxN2sww%JX&lpU)V;}(=WkJ z>Uum3HutQT=i@A!?iqlMIrU zq1G4i=$aZGJ6;a5a@LX+KzubY+HXbi3Lgv6Q=KsemOve3vv@{{q39;y7W+z71KlT)9f6zKSck5l#&G$A%XdjfxZ5prJ8?y zlI-}uohHHK?+&9_+b#j#L|rHr3F(*!S;W!szt-v)ZZj{QTdvEw1^kk1lVej1n`!#t-%y)@$# zxHCXkl)ikKfUv(~`9-bml;ZZkLhZaNff1`ipBBkBk5cngfnw{yz>{$yEneaBz~@MX zKVA7vQM{2VYNX+0M50O%Z``>2+_d8gNHnMstHIg!cO#raETYI3w})IJJsz*NLSlBvszWHTv#nq$ji78he#-1I` zwcn;-P?&!TK!M_mvfLkgrKZo+bRoSp8GHd072~)pAR24%1k&iKxH)y)!a{K8 z&gzi`l;CZ@AVg6Hwz|;v{ZHUob(0S6rdTxO8FjCakyYU*vn=d4c)35BfXp%uT1ja& zAmm-;kEvZu?@6R4zGEBWJBIRwu2AD|@}T{2%79?t5E;NC{Q)5+qNxyfg-vgLn*Ou1 z{m>o0MZ+^0ONjcM%E4QgOuT1w-N*W>5Vxou&DfdYa;k?9-wu}0?E-$w8s`I&gkOjV z<(?p;+_23vkH2D*8~bElxRj5gq56AOk5>nhf#Q{O`;^`O^(&PIb*dgMo&a7a)VZjF z|0UnwHU_ba;&TyGor4Pd;Ply3|0MF=#e^t6MW*wMwncp&4rt^NYP@5?e54|^@$a)3 ze9p{24_fUg8Y1?~bP#C%`QrKcwfon)n!f-uE@00G#y`^NPjY=0LjM0g-eFH}@jT%A z=LP<8^~WFwWv|SWh4T#Ie{q2ZmyraK#%KDYC~zsV8$Q%yYqj3yKNz!wS7S-2nhLeR z{_CZ*O$&T>q#e?kZ&U!zzU=EO1c?b<+Jaew|AC?l;G$>scuF_gc0v;fw|}r?sOtrSfmyOjbF%ZG*s09Bdkdkga;~U`7!dFLhK`x))bPe=H01 zo>3n#Bc<+dlem-wCER5wQf_a>$&eVuD;{ht@c;NmJo|p%IslHD(9dYUGdmO}nO-rw zR>D@Cx~6FQr~jMsw-puJso*Kshy=AgDfRz%l`DGA{fAlS_}|?IjbBTBZ_>@UBD}-L zPZnOVQ5t2?$PW%KnoXJc_4?+7#3D*5aU_~?%}$NngCXK?$HQWlb*Jq3WX>j|_>D$z zG(?s&t=VfjHH zish&6z%8$*^u*WiQ9W?OZT$UI+*8d>f;)Dk2|2Xw_Q~G`-v%u=7w$;2LjjxNH%^OQ zKq(5I@5LDMH}0Gqbt0ZFA;^#@oa4rj*rE7w=BodNjR5H^VMrj%8^K1NAnHna|3g{r z#>*n(`T-IB`?7Zrp{h)8Fh0#^7togGaIU|J`CA=lc77m^;tV^dh}f5 z=^bA@hSMxp?G8OX-Xa?9<(YU=eN9!>4qw(g52Itl<%VaCp1kyqKY4E6ExL}!@#RuM zL0V(nVcmr^zPe;VichF&Y--y>Bzi@9x=rJ5?Wy5P1v_!urcZ4Sq;Yr+%9O2RGnfjp z494_TRsHPZ*21ltwZOu$k#X)=&Dm7a7;?D(Un@;LU)8S!F4h?N@vt9_VN%{RJV{;J z`9fO6ez|Uq;}5fv?aSNIEs3Rt&^q3eWt<^4MR~sw7ohJgS z(0*{G{M_uX&OkU)B%4|~t=W?6^5=dTq@0_(C~PS+=Y9A6-^-lzLX-&H(BWcWxcpu7 zS?H0l`UuDLgIvRAjZ)s5`3iYS*?qF??ahFqOj{eHa~&az5fOrujgo0MV&I|Limizw zCNi?Y!yF?XP2RfKG9&@3+xzlO0WYM%U)5sRZ%kg3i6^B0H&Kfp}lRmB?XzJ1izq3bR>qcr5DLQkN>Rv)#52QSJFY&O}QiNOW?h>Ewg9D7#j_(oR)hV3%wsC8JA=(G8 znU$7%^V?PW7Hu+DfCE@atm#e`$6n)%LACy*u!JnwAGh9{+8d!jg$7r=d3s&hIDNrZ~t%>5U_wG$O(n;Wv8(N5&bMWGEO2t;m<07 zU;O!>jvT^U{%7IpzdEbv|LB9X$h3t+D!W6~P@QQ5TA_iRxNMR3R>FaUE~2EP*{KT_ z_@!~)pt85Y!-I2ZaDor{PW_me7=ax-o_rdovH590=A)fc$cNRMPi0=}Q$BGb4%A{_ z{j<^yG=Y7y?8JsWpPGYZI7786k*e|M|K@p2M+l z*&H$t+4)gXQF$;AXa4#8!ah&JFI?gCIsz94@%WH7J@gW~^|Dju&i=(E<9k?!Ea{KW zbf}WEnn)uc8t?PY_9p>i3P?%iR?MPQ)B}-Ji?a`qw|fBDYwn7flX}oS0(r36zGMYX zejBDS$&;chw{iQ;}qULD1 zytwl^@4JJ`PiGZx;qKdeI6BqZwCqgE@P2L*Kp+uMWJctr4zlCZ6)s-;Ycq~?fI*1 zAS`sJxrp)U*=OGA=C@-Vhs!FL?74g!O^hF$Q!9m-$eKl&IRa-ENRGo}W1PmWd&!Xo zEbAmSPDN}@%hw_jWVpvqnmhbUjWxNBClQ4Ykp+OM4Gs>Dkr~X_Cdb;a@X98TwLuIr zaNQm|E=&&j!pRT~`}%d*G&qQ0MC{R`Q)8B+@cbDMg>QJJ_2TS6L?gptlB>`};dtV$ z5p7ebgg$?$TG5h4r+%P`Wn#Qg_(7;j*<(Wsn(D)cls$Bmdc|4+i3d3ROP|_sDMw}b zj}|>{K9P|K1q%&*aBY)SOoAFwPCijLpm!Ud9^x541m{Aole|wFIcC5O6g)qT^?+eQ zzIRtHz54}ihlx?b)p7s72QOAZ+2Hycp@vlpXTv$gl!=Gq52VKUiT|qqasSOJFzwGO zX3w8LpJZH+M9GX z0qD||1Gen{$t_hq61@yWTDry@NJHYp8dG`=-&q|Po8!;^O#e=dOrWRP@2#hFsCm)1 z!x`2L3u=#Zi)W9lB|T4J?5?msN5?kf+tf$8x7$k88X5dKjy@iIC1VcU{D#n*1Lhq|9l&vQ1sz0cF0Z~7iM?q(j_h|yI-TW#gZK;S4foY#b|tO2w=O{|2A!L!i*SJV~WkBY2vI^Wy@ zAGX!h;XDNRke|I2L-%w7Aq2vgvI=PbjOWH-U$HO2Aa=aIx>9TK%Q=6?wY!g7t)RI2 zpdFMFo&g7X?JaTfjT;9|RgHwB9$pRjHt5dap$U;x{wVyM#D3~r8mA_`vHNzt@f-Kg zKk591EX#a?))byw$L6!+I=a%fcuh<@zQb`jMFY2NNl$HuE;~g>@-KOnyxbN2f>rCa z_f<9@W`9)J`p{WBUdF3j7s&oIy7qqdjaA6{OlXsn`tf|A*=-Nq;kjbWQ~M46QsR!*;6|+rQ7{U%$5=P-w5up;Ds2~*x#Hin%O>v=ub*+ zq!iAp!tkjDFXWWpzXD%#jC>6cQY2AHh!VVKjSD$h3?2tM31LXE#Vy4A43C@74qJ#u zJC-Z1Si4&Xq_`((k=grOXZG14;$5{vJCMRos$lX&7VCP!;`qeFS7&hqPE#CmxZxVi za%-~pFcLFt*e$&o**;h;b$vZNcZgv`@G(>*0W${`}y!*ZwL>zQqOF*U@= z+@AxMB}m_!&f_3PAM1cp|B%OC1{EAXW@N&aJfiu=vE%-Z=tS~2_AqaBA^E=p-WY7XaQUy z+`C8AWGT}eWXK^}zkXP8KrP8YTk2P%-lM1XQS58F-y~oek`F$AVG+4pRb9Q_zdXn~ zH)c?NT0_?`d-Gn31{E(Gk?ZdLJtuCD^pbvN7yXiirmyh@gDY_HC*Q-m*n^yWDhrBD zer(@0FR{OlTo&;ZCgxVl!&{heeAfI4b(n6umkV1rdf|jeKmt{0>yf`b0`a#lkp?F> z=mn&Utk-v~?FVF>XprZ2FBeBH#a7U;wM`{Qx4_i|k5qA`7@1d?$B9OX`8){fha6)k zj)Ox;+Yhnvp9QHL$i-mfGrE+rZho$fjBkBN(%niqBdE<9z<00Z7IoS|ZO5~@)!TjOt=!=H!sP%0kgV(mX_sSq( z=D%pG(P6gt5ejMKfTC*9=*i^i$WdG5>(E+h{V|3{P68B_E+d7HZIh(heR0~gYuAWy z19+g_=XT}4H9zY{6_O?a^B}o6WTHGbKc5yznN_Ia__1CE=1mIrlTu7Lnyn7!)kJ|X zt+2|3P2cpslw%xWI|i#E;lKsm|InOhkeieP$cb9~#4+*_0uS839>7*98grz79yC9p z_Nj^JO(%_Fq^lHUbHfybwx0JggEbIE$S8HedjCFrbL*YLU*FfxK&m6BNn+fv(|mUL zkc}mIjFy8QHs_(%z3(CvOhEP4f}>VrxN64q-rsWHy0Tc4UtC)E-2zg$)Q> znzN_nnTs7P?%HUvg!T-;LV6fr>~8xY1?lBCAvs-jIa_|k(pHu-#0QzMRO&Qx4*HWs z($za&hm7J-%-~O*@JXYWKo`fwCCY8wxDTqi(!j0oaO2~1V0#~+bq%CiuV`#PQ0?c%g_PJ2>k{GNxk7~{g zH{t@1$dQ)*9NVJ@@+P}Zc)6)u44n$Gmv*Q9-Vxu-E{tX~3 z(C>BmGcjb{zuz+?2`%iuZ&#~_jQ$9`JWxR8Y~i~Tx2U9_02~=P?1MRQ?QdVw74sz? zZQo$dO1sW?N0%;k{A%!x^-_|tEPG1H@q04s4zm$o!IXE}(d@(U6^->bH5rl9%g_yS zzGMjt3k&J@+D&(po!-I-Tg`ePORNh>)2=)g85ZA3t(&+&DR( zDsu0_0vsH6m`4^KZ=5rtjSXLuZmvR%>rzq%0Q_KdR|x#o`$i4ql`|ytIFVxhats(g zIera}fekcL0AV%`m6HP)`4cr{G~(9DCiq;Ki>`bZ!oq_i1m?qDQ=2{ghJwu&5WAp| zGqoLO{dXR(_=DKPA!3BVi-s=jWu&8_aJ%%jZSO|jUi!;T01HMvxC7m6+hgzU?JRei z(vCE;syUELPv*=rNYN_sm1kbM$ACnntzpKu>|sg_T08``z8hFT=w|_BFAl3-U;q-p zh9;B(8G?7ki{$QKMcq~tF?T9O&$4YFcnUwI`RP(+gB?103Aj@farbwgJbcX)N1e3a zET&9bVE4Ek`>4KRtyTL&rE{wnvEeveNU_DMS@@Z-Tx%Ys;yqQ0d1Hnxix);7B5&ux+c^)5{TzY7#aXUc10urm%%rImMYSGEWE9RUH2gc^ zMkbTS4t>c7)FvFhKs{P9G`7>N-%7i53WnwNi5ec(;IXA)}IS)2)*=Z740Iq9z+1XRh`#>? zaOI7O+5iF{>G(xbLBu7_DAye#6f8#yrI36f`p8!BNbA6>@)7443&E#|Qw*R-bYU7b zp+pv}+@jG&$cD~}mB%>$h-A5Wvyf%Uo;9V!bZ0*~piDeN@IQ@0&hA0vB`u2c2O*;O zgd2P#x#5_nXqH<^i}-vxsVVv1)8p6D6oFvoJIIi1ah~TY61ou$Gv|pyayBDIudk}r zhq<>)GxO-v19j6jp=#W_RHfZ&?r1UyTqI9-xtez2-3fHMYDA4iMy`iZ7SpEVx`b6I zvUTp0qhl<1`BkhwJFZ}WTrkNv7EPN}nD9mS*&tJslQ{`>9JF|ej?O)pv$lOYl?emt zRr3oZh1?4@b0uniEKIT4lhjT|l$+X6NuIdLhia59NJU0A=Mpxd6wt(uPwV=ph` z=h(i`{ttY5fD8@{?BtLI64+qGF!nXY(Pt;}&;V4hl22m;>5ET0$oy3f0u!FuD~#y=^+w0Rp}6x9^03JYt6jRwOw zPFb}^u%x?Q1&A72afN;f9yV$5yQOA0rK2erC*3N{PAH<4Ngtko`SdANM4xv~+6&Og zS+Lpdl%aU6EEIVs@~}FFnBoU*gGHgb6Za0?P^d^Jz-IBVl%{>&yD}9wl9yX zX;I%_MVDmkgorZiXb$Nhl}xcCr6Q$-2%$NW(nc8$iLHdBh(;ADMH-Yup}7=IQYlg@ z4NCox915cA=HM&HAlN)bz zz;e>lY&FMklNXa(FnIH%!PqSlcklCL{{T@rcV>cRm7x4X`P%!?29?I7kgpu%933Z@v6T{*FV;K<|*cS zOaB*Oh6k(_9aN+htpwrX4KOkMymkc?gY%8#@XW71h@ueUh^+hFJQCQ6iFS!MoXHV9*ER ztqM5u7kV%)T#tcP!M74Kvwd;1?)5r zx*1?oM0lDWs20%0*0+laxSfXH|Gnaao8KoEDX_^Gn9>jVm5KLm%mwxULv$-Qz*&pz zM{wQ>W0L()Kr9lvq2*@oLjtGVpj@snTabLd;Z?jL@xXp^1|~E0j~FDR6@Bz{CgDyCKSTA49#3O;Mv%S6zs2J^SjR%In7m& zxbItYWeQK+7Ybl1x9|ad_TG|p-~n|!K-F!8hwCxEo)DCfrt>d>!Ft_FD4^_-V?^HY9{U_vPlPLFBEreBe36M%0WS&mtZ&Ir+5Gt68 zXGA*jSwDaA;+|r$bwPB1!-cGg%)(=v1Nm)ErnA$?_XIJc!5U_!cZ!#-ku#n5kA)-6 zwz#)qW9Jtr{Q_SVa{DJyAd>0^$RCt*}g{cDk6;8B#UE$dKoz)kUsTbNU zw0;xeS<8Vd={<80r)u#CTBY1_QD=u+atyD%n5xzYtb<`#$ZqNvel*3qNe=I}o6K`mL;zG55xC{vpzJ_?>R6jR#K==x! z1V<0B8T?GiW$U+d;%ec<&WhX!2j&XkkW!h0S4{V})@Dl#9AuHyig-#f5l?_gE4$hJ zK);;Cs1yQHq&1GvFAXC=Yu;Rx*yNq`v%UxFqiU5u+MD}y%o2kIe%htqK1Em3 zs58HyU+_L5xu($Z&7A$-vDoZj9P^Y$IE=^{bQcexCj&<$HK4LiK=j0XOM@+o87INP z&wbofwe*FBQhaTl*2&S`vtU-oIfKJ_9JzkxBtC(@Fkb7Fd>W$--o1*u$NJ3K!V@56 zuj_m8^TN!7cN}MXw<`Pc|Iw`Ij&Yn8if4Ub^;C7EVQ`XlNu4hAI>2xUTNZ5mkm+@ITmC9b~vu0z2^e7-SnG_`)ir0%8xJ1_mDJ^;Dzz15weDT6ZC+USuI zkN$zk0>a}Q0(vaPlr6pUZY&?x(}5521R9iA7y_>4Gsvi?o&zm@K;NAgs_%~mk zIaYka1jUZxXwG(?1tx2{tZepOtFk+GN^Cfpptl7++XRRH zPepui?d`twQGJatTtKEr6d8XM>*Fc(hyF>S|83yJFqnrd^ksYN>F2ABgoRuC(pGbe zWU^bTC9Q_Z8MCr)uKY=DNm&%{K2#v%sYu^TzSK4R{QQMSa;K$qsWmG8R1@K>F7|$15gaZ-1_?loECIczq*q%IP{a16&*-Xmgu>uFlOy}ettU1T~Y7t zh2flzgef)vKsMvC{12nrWVR@#w6bc-`!07*gkRVn30^JZ-TrVeKJ#d_a+P%l9SA*R znEwtUMKoomG(*-P;2z!%4z?_f*~x9SI$8!O2k4v@S+PLR5(Q@1URcgz>T{m82QUhm z&;+mZC*i~=)SD&15aVz+Q-r;dfY~~tFSDiUqodVU>1PHwo#kZcPaq}uY5Hg-D7UX% zwQ3C)u+iF_wWm&>&V#U18KzUO{1+Pni>0dkGoVmlk8DDg;gxhJp|HK%%=l%WbpbQJDC95aKPg?$> z*ZzKg!N}zR_^5ChKUTLlYlw{C%HhL@>)sk8e$f=7d(o7IGE$IepUFwy{@y-1fhZ?t zM>-@ojBUO5S2gky5NvNs>jz`)&C`+_f!==+hx_>CTptvMBeCJn-}h4CFaL!HK{W@G zFS0%U{Qvz$jCoQefOXI)9p#`;#`dpB?>MqmqJL!9ByZ7x9v(?RT1F z;dijTayBR-pqoh6{%C_05;&4#F#ZtI0>Nk) z&}&oUMmIO1_OyQij$S}8b~pdKCA3szfn6n!sd>mvB?eXx6+EV*uRS_xDcOHWXnECh zk(n)2dp_UqXaejYM!r;GI?WYHo)~efs9OBAhp=&y(0QS~dJgi{iRqIbaG06)w2`Hh z|LyGlZk-x{-g~~-e~PSYdTXV$zu$Ok56ItQDQg3)EPEm~0+wOri*JHIdf^NBC#=b9 z`f{xbah1X;ou5mem6qPu-d{hWT{;$Pk?IoaUd}z@r<7TE$uV6Nf9z1mEJhR;s;-*TqWrWZCZ%TX> zF4T9gbl*I2r)k!ch4oRUawJ;2(`t{&MvO|g`N(GhKe}965e7*;v`EHQo|7abq9OL| zk%dR+X{DAxP0!Zkod+~EIBoQ%6Q9X#N|TA5I(Ax6s_4Y>Ofhq(tuK$s{pz>Tu~`A< z-B_ys+qp}ZpS9%gkO?xM(0i|~0rS))i5HJ1d1!#}<)lStZj@A-DoJDVw-RDX1g(z)o4ub9SwBq*SjWbW=m z39+FaFIc`Z9$!1Su$bhn)!TB+1|ZK~0BXD?9W^F)FC;O}^JM{p(MHTDK+{;1;^syg zbmAzJa0I>RVi8|*rq-D_Ym{DvNRaS^uxv}82iCqvekfbZ7z=u>gX@>g3W(?7z!pY} zh6a-(?Hxi2)aH^pL<&PMK-^rGDJli4Hmnx zHlb>f5Z+Ke+DQK(WHmReey#NxY45W@ts$*-+Y{+E=h`2>foB9c$6T3fFFC3ycj@e9 z%kaF=R2pGNhvB3-tve>RT3FY@q(x*``kY%IcNZ)#9kFFAUu{tww^bsp6kH47<3%U) z7*79oZGPQxAEJ#;o%NQK(Ymzq%nt8Zs~enSm1DtSnlu3d$$&vGE0&zwm@cU=4k1_)FrlOGLpNDHb%&$pY($zjy!i~-c8{i zJOe)ovwVVr-`{@BYs|lWhIN{M5t9HLq?NFj%JLK{mVa*PGeN=Ix$}X}dH`5^vbaAk zg$ca9M&=^C%MpBwG|jMyxE}=*Hrn;Krau&G8qKwN=N%>DcvQf9Iawc)R)_3B$cmNt zpIwsps)wwZJBvvkJBe+Fd+%BrggzP%w|3hvRT8ipkYM1iX|1NlL!OA~PtuCJo1Bst zGEO5vI14nP(OE*yZC3SPT?D_BCqhTnm{)~&V!YKu0xq_IRqzp`FhjV3bGclaael(; zsRG0o+7F4R_@$UFKchPXSOTXWllN0+&J>^?4FWU=w}HZn3Y6KrA=@Zov1&BftNh(Q zaMdw>Ce?`aP4*c6{8f`iuc8=Pyl%E2@x^HeVOxC;97yYCVYcwg>~)yg9#<8g#xUmq z3P6JEz@3AAbp8u%myrhoO%8Eh5uk>VjRlJ9Mvyh={AOU90ekaccy=^`l3#+Kg5t99 z!T!y(9IE2f7_g(^#>eE?%CLdK5l59o6wvpCqXMn?jU>YyuZ7BgG!ERv{EgNj@|HGb z(DnuN@^?5K229o{-2lLu!1YC(FcuL3vL1- zrP|G?X5LocL+g73nQ1rHk#TRKiQR)>dv%FhGSKRcp>80~&xa_PO-z+7pA#S$18{m^ z;b;cLr#}jEo;0wwrE>Vh7ugVZsj2KX!2hddhrc!6M*M!L+<`fcWxxb-D*aI-(XWeW z`&(0hkI!w;ZF&x@&pC?B1?~fhc{Y=nNA1&5tp#X3aG^;z`TQt``C%Xhz(ezL86hITm|BklA9oIiBdqv-YEs+!!7^aOsB^$r zLmyS0!zHI@8C#&X-g5;{6J5vqd*=L%7E-vL21k1RbDtgi-{c;g%crrZ3Il43SCb1m zFt4Cpln53DVbY6#xA!IW@nKqJQHZ>_D_!Q&Hc|3Mpk;kto`OuvFBdN8EjtE6m$zTU z6n0K=SdiPCqT$TD_YPd$LV;djyWyD1?ImOhL%MWxc*O2Ot|rpZV6G-hae&p_)?mKD z#y}KyM}uU%g#1p_y(UoPtxfuvacf|%T>@F?}vK<9v1gMAB}f#%b6-sH&;QAe%-Eu*fUx z$P1LzA+n@$RgyRPb)W|Y$CT`_OKu9lI-tAX`aP%$RqCJypmM&-#Vx_m_8O3i3~^_v zP9{?Zxc@qEbMNoc!K`KcA+bVC3NU1=oOzUwsTzrPyK2$8C2QqD1%Z>4RN)ac_mI9i zv0z7s;XERn?z;u^Q~JI~AV`*v9?md^0|CnE(R#?SZ#_0kqs3~p_sUHXGAZJ-VP>XA z!4V%_3!=&8%Tm(KVpspkoS?8d@4(4bMgK82L6>?XS%|nR(%I<m^x-*(yVHt@!}K? zPE~U@fDvjqmbx|kfOmb^BLRMgVq#R+t)GuXIjmnz))UBQlE zI6f zh_ti~dDU^=F?}C6raNGxh*-ibAP<_Xx(k6gzyW&%yiB7jB;I9Bxa{|EFmJvg~e5Bj5Su)I&jY;HIM@Z>JTC}85PCgeNvxFUsrk! zT! z|79<3wX8%$!NB^%zD);!(84@)3`{mF z>>g-^G{vF0HGY@miU_c@G@E^>fP+0PWS`Tc-aW>y2OI&S$uwcl6pDq8VGn8lG_}x6 z07+pKZU(YZvigs9ed?e}^#vI3=|2l@7_&vAOq9mTDl}rgEZnQII8E|01Y@De6@whp z_nZFNhZYF8ymWK*9)W->xo{%tacc%XN7nl2^nI}SiH&o*NReL~r->a2-aN7j3R`5w zgB=D4H=GeK8tPdDx{dDi2UEI$uY1UKT;rWN;7JhAExwsMjH$f>d#`nd!dZgJ)IDNXmE(iJVcDjI4w5&tEL=& zmg(i7?xTI{sFvV{M_f<10_<>3fqRa)5HjQ{wDss7RM-n`-@6~_{kn?MX$O|j%t)UD z<4A8mzfV8S=#k7|uA~ErSsqrdk`q9S*w5=jrU7K@jPTG-#^eC_UN(4S!^gZ3tf{1N z=tC6XQAeGR|uuyo|{&SoeHJqPbjcUxlTxb5p4qbprX3BuSYf_1d*rl z3=-yxHZzrP&{m7SuSSq02uQB-1^(~&h`aTE={xKv$CvD8@2s8Xfd{QVjF&+Uz8%6l z&F3Kf(7JN~z&ZdO^H7Vh*u2&5knB)G<|xsbeL)X+85}9vze+pG=zi@W?+#chckHIg zJEIkVm=1e-vb-Tl*v5#WiVEU9;>iZR3tuf2b+{%M-3to~XGq+#xL-;I1?YbG?bL(F z7UbN5_Ef~ut+Tw}y7{$vkWA;Cb4!%`9}HRTH(YOohn9!&`MYoQ9?V}69Xrl_>-_Bo zD?D}wosi$EG4j-O;rJ7ha=xa5TYWO7g)P(GF?;tt$*02AdP7GT$hkjw{`zs}n0>Xm z7Y>{#E$cqQJ-sNTc}l=VGvn9R74{X~#<(=?Y}9P38Eg5Ou60Si8m^#z%m>keEWe2d z4@c(@b`NPuE=pQ!vL%THFV)S04u353?K&G+ZjNJj>rRT~;AUrKXJ3SYD_N0Z4+^m1 zcv{iD&Fq26*DDJ9rOcB+J zHY!(2pF7>~Kj?^RI@HK|U0EGzT|hwK5eP=Fl5r1{L873bAl>*lXr~f6@ZmP>+|FuPA{k5Kw46ZyTorj(E2_Z$WXse{LgSq6!1;(F#UO`ByM7xZ=f09zmr)>OBE zY>+n{?dG8T0TX>UwK@=d3^0^y1J~?UdYattx2 zd6yglcm23$)?lxoy@Q>wLyx%mSxlBhBqW;ss_@AgGba|l^@KNMAhEn)Axi6=sgei0 zBL%)`e&kyY)?n0L{R+*a3hX55N0GDkklQFSC!C2kDOq5wUxrn`ko#hN=`Gv<@Fote zyUcK&o#x%KfeJ^#9w!9M{iCG(fa|lN??sPh2{7U#8ihkq`)Yk#o3{7aY8-yHWNB0Wsd4$?8~XPt zZgKfE3N2Z)d1F~e#$0|Mcw|D_#TX9JSe_{yYIS$03^QL9H=NcPuM%7KxcH_8+9Q=}VbZd=lqwJ|@>YpJo^-oR|ExGS*@5h@_~wj{AGwi}GdH z+Oe^vo=N>>z97TVoIY=!p6r1p?!C?ZdlKy*y?e9I(JO54+liUObrpPH&2s!K#k35? z`iH z*Ml)2j*LD#sX*%7Qia#@x9m!rp2?X`Zb=#S4HoS<7zi5` zB`n1&cpgq6VBlYzGRA#oZiG79-HwTw_Z+t!C6eUi$WInCt`1!g9v;qOqv{p1!^twK zo%nANF|nTLCG}T^E-+E>rG0|HlCmu#;%-J(ZrHdn96=J%gq5?xeUha_a;U4*%`eqw z%1%aaEvXO<3J;f|hp8h`t~3kBcp%YKaQfXT8`iEJ>9#OUByNmZp}Ok&&X$!;>4tOm zy2Vd&w^T*3fDb-~+=YV{&|Wqi;7>6qe5g97ewG?O=T9GULhrEnfG7JCWL^9bA}cJ0 znSi?zf8^*vW0+03tf4Q6sR6dqBL8)SaKUft;M5e*o$()84Gy#^2a#@su3p#G*>P08 z!$T*OgfB2Ppe`3YDdE)Emhn1n^%1nnXY|SCK#N+eHM%@Q;!>Qy5yl*jwjDf-Pt>l% z5PXRDgyau$Qo?Jb1n^=wVvB&3yvt57M<)XYL%TB*Ls=5dzwUJZ;COv$ zDsA?iHE&+NoKehEut)tMR=A2pb_Il+#HnU_dfUNP;*Q zYqlIR?x7ziN;@gk>OkJLfJX+s9mFzgvSJ*P;@fxK@^>wSZ;*}$Q%LpC91(!Grd%|1DTvq z^&!FV(4j*XsQoQqZllR9g;d)#^JJr=xc0KKw?abv&_iY?BgH9!=Nw-{0-B3l?)K77 zU`yCQ_`OShW^Z~ySz-b7o=WD=Mw+cUE^%mcCyYG(Ab`oib|C**H;%T(b`Y)cw5O*s zE}_s0)vgO6asAldvKIK@&Pphww9znC&P(fWEMU3g!{fyS-k_BtQabum;|`dU5I4xG z@;+J|;~LaX2D-X|F!<_`mvzbZE314_YZn$C&LitAnUOZ;abypQa|akDp^Pi999KQN zqVNbnLW~y+pia4iwp?a#BvG5~J7nA>_4nFMJ>Ncs4nICSWtRmI~lP>Hkq@IO92J`Y`WjEt_l{kp*Foi1Jl&h+Eb#N-<7re-p* zYe00*$0(%U>r~4S>&7FP#ZZR@XG8DwCjkv7vi!eXvvUFH5iuLQ4Avt2|l~ zqw1T>ankRB;F9c;XsLXZ$-|q#`*s5eN7=k(+Kic}mVHz4y5c?!BlS!B&DTb{I(i%rvauBhl2$S@Sin7bYhpY@%`wRV`f;D(;Er8*~5uka6RePWY|lkjskkTT#Rt ztXQ!kQb!n0Gn*sM+$qLDA--nn=#$2=yQ|~xw8&dqn#lw01PGlE`9$RG>rQC_bF0m- z7boH3p)!MMY9Ty7aZzd%LJ?v+ZmjY@aC-8g4gpQRU?tkGf{JpzxF}dx?U;S zdCwg%_x;|#-+mwb{PUY83texE?bus-e_dAqDed4&Fs>)7{i5YR-k9%i^+c z-M)PVn}+DU3^r`oiLzgN+qN4QcHqugyKURHSFQ1*YE`^uO9A^F#WE22u=LY!nrOP7t6+o&PiPz{tY{rcT^p( zx$g9RvD@OS?k8Rqlu2M^tSXk^1Y>wMYx?vjYcJ1vV)DpAIyJ0X=z8&yFYg|bqbKiG zb@d92zR%EPEj7$G8ZCTgV&B_x-Ke&G&*ltR`ASf({bQ2)pQ(p_lvm$ja|kaL{~^2a zGt;X2dfRm3tjdjnwrzU_{G$V0EiK%wqqfguIIPmMb7#Go&M?)pQL+-_5U1`91{eO% z5|-uw>JQuv;?sB&giJp46})S26KV$mPd;WH?HF;QRr=5335>fM&5m6;wz(vY9L84h63&IQw;+v1$!>+dB<0Swn`(FI3Bbc2SG&SmtSTk zmPa9)>_Nmrr?%~G4bH=SF!3w>C4v9%ZxOHlTju}cFMf5?40K=Q8dgCoN*MOabt|Hl IZ9Vq?05seJ;Q#;t literal 0 HcmV?d00001 diff --git a/Docs/sphinx/figs/thermalDiffusivity.png b/Docs/sphinx/figs/thermalDiffusivity.png new file mode 100644 index 0000000000000000000000000000000000000000..d173123857b1494b92a66ca9b8395a2033efc146 GIT binary patch literal 32449 zcmbrm2Rzq(zd!z^B_cDkgh&z@MRp=WB4i}9MHv|xSrszVAdwN0SyCc0N@xh#mCR%o z$&Bdte0N>-Ai(PpF=*20bkoErmj%->0c|kV2u> zr%LzXmPRHFmtX+;#_FKC-pLTLPZD+Ib>@gQtJ10jev0Y-4 zq8m@TxjDPaiHke@^95o~E+@qKbeDa=o2+)$G@PQsWO!_>?jmJv3+X04LuXS zb$Xc?9$Q%TeSG`ZOyNH2hfmU;+@>D8J5zEi{}immPxrx!=BD=)rt(#A18 zy|^dOdz=2AW-ayQtP;+f}77+7Wp$tIE=UF{s3S zR8{D@k%ueZP?JuTM#1aH^2?VmZUab?VA7mrYt4 zD=vK>Y0Y`b9~Kgl87$Fn$fGN|87I)R?5|I%?Cf0M)z#J7$%=1EQYWyar%g?{z5n=; zf7`amOOyC$`;a$tivtP5-UGrhTKxk9hdiW0LS6;DefxG}Y;BHPqG;+++_hV`!sq5Z zg=5}~IMOmPahfd7Pd66WHH$pGED$d2eVm5xN_6z{VsnNQCr(^Eayl|3Bs-Ro+SJrk zQ&)E-MOIeU@U^$8nq5|p%`MuX)RJ4O-ftQjG!7kF$FqLDwS&Wv&N@DR{$&pzKBT5> z+PpdK%;y>-my=PyS=*Z;IHed_SS~(z@Bp7YoOPSBQ${9l;07a;{xRXXRW2457BRU| z>a2>co(+5WQ$J>9&yI_x;;cfZo{>uuqPg~;UutQWext+9WUJ=a(pHd<%5Ls|eO3Uh0)=I4E1UZY?;~&R}hAO+!mt z87^La@`{5HzvNDw%}9&B;=zNowx`-MZKO_UxpK=Ks%BxPwLO16`+fA5rcVQPMiZyn zX>FgL<;?x?)MGq2GLpVNSve{ofRbyHPyI5AH@4K#-ac${u;KG$gYrwP*Y96qE$z#G z(R1&-ZJcLubEgB>wQJYNEyeweEFP>+-1B>BQR4Lb?dIh^oMU5SCWWUt+j31WeyNQ$ z_BDW+niM&%zq|j|-PRut4i46ijtp4PxPk?5&HerTOQ)K&O&(ZLQF?G(HCyvasd;&M z-M)Q%OnrEt+R;*3X7%W+Dc-?}E@>e{;JmbsKv>z)qtWfEHpMdNf<9&7w*`~|!E z>f5W-Q%}2F9{V#GdwY3h+2-fx4|SEv%ejx-nf!qbMtfP!E;Ls5cKpDnVpqYWBD?0j zz9%O4Mso3O*|KlZ+uK{N@AGHDaS0h2@mXF5Z2fCb)0KVZ7^$eJu>6Q6}|Mem=S5%~4%S4*Dq?_ZO1^P_oZ2HtY2X=uD?Z#O=;hALZX zj@2Uzk-f6$sAIdbOP#~nDVjHaUNA(|*92gGuKY#7nf!4i6cpr+4H#xqEn9h8 zJ625a{xWqn!&@ccPvuL?%5t7= zQBzYpy4(c!qCFJFFg5jmuj&8Gzj7wxtW$Iowuw$R6z}Y zn5G785Os!bVr5acU${Yr?(P&lLyznBcgV$)j9Y_W$B!R(c6EIr)9qt0Ji0jDZYg2k zA}lK>htARa>Z*WQnbYJz-H7~)bKc&4Z{wuAzZh=zEm}_*>wJ32*wU8`J2mj};~gb# z!xz4Pcc!S|#jU=kQxdi{Ysxz0eCOTN_U_H*MSW>MXL`)E%zFo}fPRtV5oyOy9BCR+ zdrzOv-LzR*nU~k+JC#|%Nh(rqPJdwDdMG#GU{;H_;`Pl3f^~$NpR6A*on$yF$%$%Z z-dLB#YGnGJN6Cj{tMUEZYwn!X;llju$=7waWL`cyzf9sxpV}H8Sw`|@?~6-0PwjSu zF;ZW>desC+fQrJ%&VH%fccHSqoqb|rf~)&_aO71Vj^fPTK~uLE_2SS~TFmvjzK$0ik@pB)z$-L@_0)~#DloO>_4FLBS|*kNy& ze!obvs-}jPLe{{@`$%q?<&-o1ub7y4dBrIq%T_Uj_$|(_c6D`~==A93EL~4gb9Qzv z_QJ+x$4AFjeSCas^Ey$J+XlL(9}DaRyl1{{L6dEfNI$p*eUERT?p8*tq+80jLOQ;x zUq2_+yUTqaTNEZIC!cK1j?oi|G5f+kHJh^Ml3In|VqkiD`cT`wTOm?tWJl^q?2 zO-xNgu`xqhS`6=X=qK#_wog}g9i?ixIiofIs(?nAOqb`LYcglQnoNx3xZ<<!>}o|x;_sM zX87IL6Rz5v?89&N*k15$U0u+lM~~))(xS8$6KX#;7}t!;XrvZYaWg7+u+pwT_1zdK zJYlx`bbZ3En6oFxQAvloE0lxjxr1wC#a{q?0TR53-*M`J&(9xPrWxgDzkOU0r*r7I zuBWFg_LZE+Sng#(At7o?ytLhgSFbcF*-Qa~+BbTD={O{hh8*Rr9(;S!o1@R!%+XN{ z+eACp|J4<~xy9M;$MstpmKGNPXwG;0E$m<4@zg^ceT2XbO_v<|4mzLTi?iQ<{HUCp z_}uv9j4k>nTSuWo^G&R=aK&>oJ&`h9!n5b}Zr`quaq5adMeB`HUQ)5Q7ZKVYA4oRC z`Sa&x-A7pOn-^dB)=_lIyarW)Y)nsY@5;J5`sTa)Ux-)yy72a`4tAjx7FbAV=$4{w zSp1ZYoIgi8io{9(fAU0@qGD*sTA!pC&Lpt+3@hUQrq+twC=^B+YwOG&XT zEiKW|(Y?T(w6?Lifby5Q*2wYw11dC^TpbDG8+WATaP-@h~#eE^XvtE$w{G54-t82A6Z^@i9H5wC9l z-wT$Fxu#ORy}f+#GEQR08bXDL3Z?`EjMa#DiJb-YWFO63<$>)PB) zUVqV#g$_CHmv@$UjEn7b{3IO6$S);j0BR9>?{LwLS2@`Z;#Stye4?VFRh2<>)#=)? zRqU2Nd?-HEXfE3QW8D>&ufDjP!MonOu`Q4h>z9R|OPX=*1_iyBXal0==H_TF!E!qR z2LVR%&TxjW+s2H=z_R&3607pB?-xBs?-d1NCsHsi3H(uuQdnQIJom8|?eHj!V z)TE6^CF8d5*uf?lr;}^;SZnkfORDpx5rDE~De zAU~h!y3qcEFVTUPQzmeYqD_pwy-R?t7|DkCm13{gDxiPY=tu!mtJ=PO`&QO9tN#rO zv0(PDmy(gW#HcDs%eH=fB-)yG7dQU$?3oJ6)IBx?V5>&oRm^TQdyj=PT@-DKf2p)T#zRoZgw z{f7@m9~C(Fhec)6@M%>QyxyXX18lt>zl@sN>cojlTMX_-zsebLFcyiy#y*(L*Z%mF zv3Z#{+u_2~&rju*3BB$t_mvM`%gdmlp>Z$P@?AHMZNx-XmD^vkl)3q~1UBG&muF=hvv^@22+*w$EO8shOemZ9kbL zMtfEZtZ>BS^l4GD2R?ABYiPv2AI9%#0FV^$d>Sw--6JXvoh59YCGLA(US3at3+37O z-WRQ-+<+Uu78e^IojC9}0rQ-wX73#t3A=s!Hm!q!n%Z)-%}XfNpGw^6KylwcJ{9YI ziHnUb2(2)>I3io_&M_yawZ-O_fn~J&UcY`lGGl)|jb27Z2CbZ5LV}5)3Xm?)Ywdoc z>0Sju<0ns_YTKa5ic8on)hEi!p6OeDS2r=c`A1r6>N_0DK2b?YW<~EQJ>`b_`n~4n z3G)l)pU|_mJNKx_%FDlKZDmW8bFajak^RCWe|8P2fQj-RS`IHr(L=uwkd|iE)z!V$ zP7Z$s-8#ap_5!1Bw`811ochbd_b&v1Dkpy)PM4>}RFyW9W%M`;)$FvO^XV@ay`cKI<=NI^Yn<4s6N z&{}T<DOzAUtC(t#K=gYqoZ@zp*atJoa)d~c+JGzoClTeU7^Ea z!UJWUdnlCFwYqwGjHnEZ*m+p0%P99PD=q`5QS{iXaF#X2u7mH&&Uw5|R^~-B+_p6S zd+C5=#vWtfvP#G$5}x0!3LHDld`o=KL;mDg_ryQ&YIZ=Vu2^Bk%D^ob`a zDe0P*5XJe|+?lDd&cmZT4y+Ms+e{u%F0NWz1E!oOk(XUqrjo4WtCcXr%g(;a*w}dA zX)M{3pPrn_N>h(NIWi8e*@GH!?R{&Gaf+UBPj7F^nO#tYb_B2Cyns{e`|+iAuEBrl zA{O_vjt*`B_ZQJy476A;1dNOy1yGUoncMF7^9y^y$;MPQT}UK4I(L*~U5d>!w0T%} zpyooN;g^wN-*bs!+k`jTB*?QECsS>0?ThMg_QrDGbNAW}x?Bz#7!YRX^Zog%x%P^s zAQYg_plgCRZ1M8Aal6c!m<*8g2M-@sqIuF6>g@U4tH8evmH#x_=jWQpi#kHrb+(HO zQ`+(@^Uk!~P6gMX4pKdjq8q`hyur`U&)MC*`thlD=fV2*PfAMq!1tV8TrOnU^(IfD z)x5)zB*=R(P=I{?yeGGeV_;U+795y_|B_EBAOI;MO-)UD8M^Cl#>YdDq8}I-7#VFd zF*5^PTf1}TPL$|?d-oi&zR@!<)Sz8W3?}>A`;=K7EY8TtsAU_|cXc!m{jh8MiE=pjjqpzkwrb8$PIx_AZ~l8t#U z6oI;EF1)d|vkUz>Irtfs;1Wa*0A{S}4FqSdJaQn2Q~i@SBQEwkiVf1-SzYP z-m_c6@0aKFB(~c`tFdJ$fNv69k`mMr5I?G6@6EgrnRbf zw(DG8OX1k3VuR?%r%WGMoj-{?n3d=B>50A@7EDX~*6f1o%fpf^mVp0hu>)T(P%m2+ zl#`RQ9lHYc{CWa)o`bGe^8Vw;@1p1&^DAdd87Bm1t@5vpa%aB^1Nx{vJ^M{fPwz_q zwC&5wt10|DcOLp&qTO+J|KcLvH4KeT+Xm9!I`OipYiGH$m+enpMag#^)OFVdZumYq z$s#5uR@u~a5S;m-s70wse}BK`$J}Apt$N$@OQoyD&wQU6HXww>+R4dsKyZ1O?p)i6 zT~OZGQRCUyty`DsI#}Nar9x0wE+jINja-03hYpQs^9*$;;hcz8wl8;RUiwI1Pg;}` z^ZonxD*>~gJb4nmE2|=V^Nrf250s{ZP>w))uHLv2{Fg|qxeC-2{LI)=-^>yU1QgBP z2M;py@bJ9Lv#ikKJl|9e^s|OyQtFw%N$C3P2PYf(`v(U<#~52Y2Wa`;cGzhjLq{33 z9u7G-;e!Vce*fOPN7VcYbFl_?Ckh1>qZ16yV#iL-)*v6;f^l9IL(2 zi}L|rh^}epRq_cT3sOWRyDG3Z{(h;$X34;ZGm|kB3w7w!a@T6Ve9`s3Z;ySbuBjO^ z|4|VaATTs^HQ@+gwJJi_t@~eH^#stMqCjh?!or1290|x`4M=rU>IB2ewHtCRES{Ts z7dds`;5ig}`Et!rV_J^ef_*q{DbNCadie%6wiV#U7wfl5to!01ufM<2_ZCO9k5Fik;aL6l7knKw=SoP^E|Eq-jUlirPzHB(kwBz_& z4(XHcD5#q`V|N*unDptVb=K(=zBihZb17skh`+ zEQRT(?S0FcOt)$(Dk?(Hi_3Nr2E|isju`9mf=u!K`}b#~c@>vq zExfa;M*pFf;3uYrn^)nkR0G5?v9K&-XJ@AbWMpi1adFWqmD#k3I1qA=H*emo@h}%# za(xjFAjDoEjXRVJfJbS5NFZXp{GrdI5I(BWkG8vy7#28nZy0UQyM+2tc=|(X{pz)A zHQe0XLNDE4!Ai-sDC2@?H-dkC)YY+9P3@!UV|%udmaH^C!1?TwTWVZHA9o#in>BXa zkUMQ_>KkGrK&*n^qQ=ZCd6dQvluyfa5M>=J_^#ilvU5cT+xAT5dW{tiZoK{)zpRcQ zUkTVm0H@JI$0Hfb$zr{@R!Iew6>vKU#K;B`BZP{P1k}ZxqKqLT?!=VsGHX6wqIL$b zcsJU!ex8M5OOEkc0*Q!N@cMrg#7+l&{LjK!-n$=K23;p?hWr*$GWlX_iTu5+pnQJ) z-})C$>lT-hTKb(kDk*xm_gxpd<_VVJ{CmmwUaSL&5yd2xOeJ%<5)WzJVR$9;c zs>2D>Cp%#7doOCd@U&vG-~5*Hxjwqn9}8Bu-#@ZPMTJ5{d8kC$uJ}mk_jmkwl>K-p z+fj)x7IZMu-qpry_l?ei4)$I@p!6Jw;mkyzx&rj$?-LWxfmvF2LE^$SlyV;Xw2S=k z&!0abM9@e|NliQEBqpvmD{{Q&2cj}umisVzHjbmj(&Ov<^QMoE_ktrF^Bc6Mm6v~d zVnJ=2brmI~r^jk~5}a){rEj*|Kepd@s!;>d!9A^y(Kl`gQ9=p|3PenKIqUyJC$kH> zk!C0oR8&+%*$B-t=wOz#9$oF(rC(l+(15J%?6!RR^r;H1xC(fqI#WL_`{-gYD>wI6 zN{G~n`b=9H$4`tD(viTHODP)}cw{$1=HCn|~$n6c6mNn_3$ z=Y(S9`8iT%Q+~ef%P?%%$1SsSKYxazw%IK%qW?q=eRz}sh2Y@JXsvXVnM*fioRSXJ zzEhuA%n{>~KDn{{*N@tRmIwZ!emyp`*zM0N=cYL|+I|3NCd{$3BozJ`6Bm~VIjNYK z7>?~$0ig1&IdvbWF78YYy}xKiIBw?QqWg1znxrpQ0V-Tha4FW zcoh~Au|~q{$5DNNCn6f0nHp-OEz_Xkd-&*4H4wiVKzd+MP!Qo=L{U6i{OEWcGtlNA zGQU4wnhDL=I&^0K=g+L<-&4&8QV%36RkpPq!mSM@+8yB9&@a>BoZLm9W54{`Gp`DR z-M#wqgE>tXltyrG#+q#~F`@W>pP9MeGJR!&ZuM&YD2w2rAihnT?liS~Vg2CPc~Oa9 zz$WS)8@r4iVO;)UD1?&_g`BP_jbz)en4JoiUXS}Q$?07 z7NI(Dq0(H>MOvq)e8Q}>UptD{qM#08(Pw!P&Fc2_kp6Ce5OL0&L-#5H$j(F8)QOi~ zP4*cu#+i|pE#}28Dy8M6p^2oW-fLYf-v3XbG zernBnIo7qjO3|Y20ER$%8rc1Za!pb;L^f~wpM)_@2H$p(m0MqmP)YyFxsA)w_pPY= zM%Bk60F(-T`cMJ=9kG0SQ)&060C7?__Pcn+MC_d&Y1RHUsQg=Bxx2!@0(u1#>Z|_C zCBWe@9EVodZ1?XaW9d_EI-Z`MAt=G&z;v<2NuXTepazFB?_FBUE9P<36L>Baw|`{N zf9Yoiv#_TnJ`rMiL znY{F3pPvtjBm_zY8$9gy;<$hMfJlrMF<=^THU3Ykmu27L|EKDucyDQ}c0p(ATRU~x zx!HASrR=ISS7T#0fBp8YdgAj-!Cg4;Sf8_~u$476d$H2?t=HDlVy0AK?WN3~la!PU z#Cm%F)FV%B^M6QQk+p29G)_)VghlN;7k*id&3Lg7rKP3sAX_e9v4W3}FQrLnJu{n;8!N5eVr{9F zv^_v}%)@#JnZI!1O&%SOlsgW4mM`!qtB9N&7dU%tsSN%CaTqOQ_wnOGIGu^to0Fft zd875YWFhnjWWxyQl8(Guzm9uxZx#WFR4`M-+O@VDg#%~q2v(B76(^4*B5@W ztJEv2_Q9h^&F7MICiIu{sl^G<5p+ZLbDd20ZO~oPCSeWr3=P37=w|+jefSb+$6ZO_ z3HaN{>JP*(ee~!Lo$fa}+A-ZbFc6IQSMdAA0juRgUY@vavJ$$&1D%`OspLU8CS;7=0uTLX5zDz> zaMv!=@!9_@j=g;Qb`5b{v2Ucs*eL+*n?|>dVNaPnJh~d1(52{T)@d3;h zB5Cpo9_Zdg*oT_2f-<+TkcO^iQsTCSVufpL0)UH>a2_&pPrIdmHTXYtEVHa{eZ30I zs#dSB30|rFa_xz^*vo5z+R_ZQ+^WwR>S}AX*M~<&s$&-k%D)C1W|q$y8anAV+_ZuM zV0#`sE{iS`9ViuObLf4(fRJtijf#8d)eF<~yh=P+VrZ6|q3~UL@L)T+wU7za_4M?D z3b*j{Q=u?Xt9}u&_$xFW#63D_ts>N@U&+lSdS)3%mtM*fW8E zR=wcj)qqIa{awYbs{o*$fonYeB~6PE4d2e4>q0|Au~O*)35@#|?v5842ZG`PsnCmy zixZ@}VZ#Ra(0f8l2~I&?L*wL0VG!H`3qP=R;);c;@7o)832YsG=wLvgO_Kxb9bxHvq6T!Ha^*w%c-+}zcs=Y+YT3! zb>KNQbahiZsN0Xyo;7W$4OE&8SBK}5E-V^({l5LK;L}5fI!LaxSy@@{3mv$LT@OWV zncjNZdyx^2o9v$jP<>hAD4VR8_Z&?0$j($3kvjhN8PM64q$GyJhdV@$$Z%fjOVe8+ za!y(V&f||g^I&YtmoL|(72;d&*i{;e!l`9*=FAzgS0N@rZfTix*r4bgc>DH7;Dg;_ zv7b?jC;_NWY0!fxhAZ9!dOSG$&8&3hdnH6I5kAc-W!ed$Oeg{64o6K@n%C+fL5w$d zoLYC`;ziPoA(5$$9<%$OL^pjK&VPt*C_Msbh`-UZ9@xHQDE+;(0C$tE%X@TNWZLnP zOwEb`u!@|#^2!s*zK3j-y0pd10#UGSK)%@Qlkbyd1Qtnji^ZtS(ozNUQcwNyQKj^| zcl)r+VbyL0E?Y&Bi_~mW+iR@G#F=(MUzml3B~m#-NJxn6UTA-*?Sah?$tD6VN1cK> ztpM||2eLADY%jJqQPZ%idiwj98#Grc4`X+w_%6&4qncv%^T!vGqkx}7gSv9(?p>Ru zUq3w&8)$@HbSFKXWD#JdsGFH_0qGK4NC`+!7sMx6+1f6nqYIb(&Bw&U69L@Z*P#5n zWzPHg9T(h|$UF^o_2;Pdl5g&YmR~9;h+}?W*DSPr)vC|PMo|oZmmnDv1*7`{!+JSz zSZYEKxt&V+Obg)~;HY~?N4E~sIXrggSO;KT4f}>^{rV&JN#z|DU2dyy_ZNzr70{q( z5E7ZB;B^6KL}Y@wpM%M*wQ>JV5F38<`0=*scW2IDSQ7JoFAU8aiK$qg`pB>_RTxn_ zl5{H|S-c+X%re-F67kF`d7sy!*)OTBW6(9Kh#r3IO*FH=Jb3sED4sc`ArJHOp926! zqLOD7?SQ2Z#hz;|->Eb-sgv&?@LjlYftV11fy?hLq}-J(fjm(o3?XUt>ecBCwp(Mg z{?O9kT-=UGwylS(0B46jKHiGz%*f14jnXh}9!ARJ(7W9~%JCQAysTs}A}0Uei!zkk zD0EL;zR)TtC~SB9B>M{&^V%^abhH_cSQfje`<_JG-g>Ow0oj+Wld9n>PEfVNqy9Xq zudja+j4o&}FM98-EmT{!Y(Z4;?)*={8VS5lDGI{L?sCO}`o!>}F7KQw^_2$_WX~gC z!(452{ztvJD91i@%DXlCC1D~?CCl|(IML}QljDh=h-~6!SfP$l0)Sv(3ok#> zkW_ts*(#_S+g@5G-hwI0+`6zoxSl9)0Hv-&)?qzfRQ`j=vE2wRhmMWmH7UO_4%=@x)GGIU1h&}B% zXcVP(MwgZ8>Qy40(Fvz=1^zSEY8T4hvFWcgl${-jKj=%NXTCaq)O!YLF#JcNq~uS< z|Mq1Qb3?tr-#;6AVh}4uU;m$RONe6;ZxD|n^8U1KYqn8q=N%CVt7jCPW=6@a-Evk; z>UYwy@ytuyvRX8t#=Fbp@BKI8&lpUf94Ebz;o(p~x9IqeQZH9``EaG&FsgzcEkEEbTqC5fzj?`eA|KLip9c6olL2Y3*{-_s`);hL^A zoyA7^!vYE&ZsFE2!FK^!xZ6^$4Y2{5`N&JO?-dj%3)%I%yr(jWv@46C5oIHqBz@DN zJ=vPv;^t2_lH>uxCZcEwoAQR|y|}I;_z;F8Bh`zElXJ^nNRR)p_6LdPO!R)`rMWbc z_gTfj!{|QJ!UUAi?&WK6S10S;;Pva*(yUTk6r5O-MvyODa+-wW7Ca zHz4p3^Ufs(l_ng5I`Q%W$c4mkFolktB8Wd4*^TCkA}mQAAnQ61A$Rt34Ly|#&n_1= z*#C5tJ_O2Ak^idv^vT?3w$o#!sRv2^S1b zWirBstZ)!mP-^(4r42pLouflAj;vWw)9hBs@U8dq!pMdHcNwRVmSE^I?6rS}Zkpq~Eehc4ZNV+S-{oog=dt)nGSpO9- zk2o?Y+QQwRXwnw9etuyEB>>04Ns_R{aYmSe@(ds_7y|h|{pk-UlC7~C9uyU6pgEG1 z2CUJPmgxl9Gt}hNo{pF2cpOH28RvB9csjY`WX}>O3$}%_s+y6`}|l<7|o7TZ5y|4rNe#~sEOi@f@+hBgft0) z!|gG0ln-S5K0B*{cwl3`t%kb~&?kxOVl^Vf0GsX(qD$Z z5EKvi8P;D}z5#I>+wvNN~lCh_pgVWD(MN zhIL2_5Kb^jHM^zg;{a;Iv=!(vg??kQ+hy!s(nq9kgy{fI<>cRbgPt@c%L_|hb1F}1#(@+Ai zyoGdb1|T;Oy)ZlQ3S zlaoscj6@b5@*yn(TtKOiu~p(~-wi&vQC6l1%`lWH5?&~C+Lr3h{~+q_-B#PMVgHq= zuPQdT;Lpou7!0~}NdSW{2w#A~1mP?p(nNx!@Z7z7S5KIgZC-&+OLnIP@fJga@B(cxPG`?Vxr8(VENKZ}}fPdoGWW%YqIoU6~ ze$Tg8IDdJ=w(lz-lm^z(c8H3|fZ$h|@&OfB%Wm97qB2+AwLOSN`w6mu;ZHIO*^ zijnMe2$XZEWS6ikd*G94KsqNj^7GElBcJF0D_j3(pye|?)YnmM9)h+0$-Kwj|0m{c zwkJtK`=Eqop2-gn#P1yX~0`(FEsRkpSol_FfgFDFMb3s9)B@Gi~&2KUMG8;4*|Ort5xDK(_)T zBO~V@v8skZ$_OL6#iDtN9u16GqSQv#`uyl8t}iq+Ch7o_r)Tfsi+97)GsvTnJUG0X zCpg0;ZWCIEk0w4(!<*gCChbf{Z2yMy`uq#$eZ72~@VN1c-zzCC_YSk8IiA0C zX%#^AB`DeY4YnrH|Ip?PWp~BI``x{nP?G+lg;&!Q{erSiEaj9~@FwQ;X&}6BZWMDbOZP3YJGN}6bDUtXlh&GK=+4PP62c5T^%0sb6`2_`6U~`fjFD#Bj zcqigWF{<7MfqmjWx{hK6G(eL4ZP)Gh9sBU`7)f)Z*riWQlLikE3NCrgd_HBTuv16q zn)JV{K3Co78*sarI5?0!hC`iV-`Cf7&_Vs67ivJ*l`D*-l2Aw}@#f8&JurPqLJv`P z5Qwa#mj$@J3CjK(+S=Mk_74sEBF18@zV|#I`N9_Af_&;>mmB78Zqj?BHXKRsFy=;9 zG52k)HcI?+MDW0R^C}M+MS@*t-tE(jzAVdZ zib2Bm(9{zIN9NwUM{<-mnm;FgrHdTpaQSh^VXV+yvaWKHRuh zuV3p|4D|O0^6d%kfzc8D$s_!g2{QUbja9g#p=MZ%qcHkFV@Q=4mI))P_j8o;UdZSZ zec|FZfUrb!&dtpwNgW7FpFyrkpMa(vN^$_>4eCclb#!z*t6iW4Y)^DtkVF6!_`B1@VB`ix zu2^x)otgsn7XZOF3JI;m*v&uYOTF3k0o0#i57#l z2n5pt7m*Ht(wGJ>5;?zT@I4@Y9(}cw85NRqeOY*ByTv(r7|dkW1|i4GNX4@NEbhk1 zGchyU{PTsIZyY{+_%nnwvee~fa3PQ{V}j94BtZH$P20*xf96~_q{W+bfe(_x$c)$u zq&)5iPTl3SG7&hkS2u@!WycA7f)_zzrYNo40403TJAs{+B`k*C38||LgmSbGHFsYcI2z5mMyLSCLF+?!bVui7b)=p8$+bDYk z4EQshqdU)lv5(Xg{FMJrWL*cn5On*;@)KYvaISqm3B zZ|wM0s$t)b*!VXsEv%3vG<0;)Ymhh%;t9OWoncXK$HsW zoe}_d3#7o%HT=~(f(|jV%qJnC&s^a#z5ze;QZ;*>L<#>0vE2|-A)aNNPB4kYzxn_Z ziHzLg#6ABW?z=aNcn?xc0I_Y0b#ZYjjsCj_2C{8-?Gl`nl%h(1xcx?KMi?9Kwnjmj zt=eqH!nC$Oe`ddL&X^xm^34t#9AvLyrR{Wm^zcschY#7F7D;}^TE}_^4GeWJi`=-i zwXXFBV^D2tOR{?GVT=1aeVszaL+>t0SIyqJt-N(!aZ6V++hS65+bn)({dVApi~I0` zO)38 zI*B4j0vfV-ERqhuBVan#If8 z1d5cqp0F0>S7mc`iypYloe}2r(U1d1VfAT;^y2dJY)_Q{Io74-1Jo+TQAoM*ZBRT| z#>>brQ{)UGAlt@ua9xdyE1fH^qDX}QT#Cp{y?j$)Tzqz4XlPJ=0LAESyK(%fQ#oNT ztJ!LHIvhlR;rf7_Ty`H5Q&4`Qyiq=rz@6fIhjaQ6O06kh&}AF#wu~R`O$oac$Xes~ zOP_7>$F+f(v3p_ciwkS$jjBI&rWapanbn6}X^oh&mh^0*hvCFpWzmS?Cs8!5Hr>lq0)vBHc=d0G z%`oA41L|@pW#Jfdu|kfJ5`a$o z8Kv$$_}3u_xfGhUYsFpHgGh6L`LDw$GTLZt+*2?%QOk1km-jaEUxTZ?no_K@QvIqv z?XCbpN-zzTsr!!0YlQZ7lx}j6dVSMqNKT)at&dKug@zaeU`k*W1kLkM8t}*ikuKS| zH7GS4ir&(cjqrhhV+g(hj!CPZoSY=)HTda~F-dS6*)#k=5f4s(c<3_e>f&<6;Hb$C zIH8z0(_9!8Z#ngm^~&O-j@mjSdd~QG>UHaa=!7{e3LaAHEba1+kIi@!TZe%K^FJ8aPUE;A^|{xGHWf&JM+$6`rND~BUDZzbOIuKUCz zJm)xcEyr*j1AE7N)^8hI3vb3{{h%^0@nlNs>JmPH;ldL4Z)URQryY(wR1OzqSB*~0 z+276GcP=$cnSMfX%TPm@-+{2%6Y3S4)2Dq`RBXPs)FNmcTP}X>+}D@95~o{)qX!1l zG$aaoc&~iY7)r{~?jBuGY_#D$GcuSg;ao)%y&2;q%V=i|u=D;>+M59kM;3PSP`vLC$d1 z`IVI`N^Z@+CD`Snj_3g`EmciT8W;-@sU&QwXvm}#d3ZsQVkU%rsE!AnqaAeIaG<1+~|_%BWJ5}rf> z`cV4A^PxP@H<|ni*tvYF@V8_}lyEQ9)HM$;H!c73xDY+}y(>6(fJ~xU{s{YWR2nIE z2nxZ-@JGitgXf#DD*Ei%(S{9k2JRT@>OL`7;Hg0DBHOR;^IPA`Np8X+f#L1RGVCLX zZ!#a*ZXq%PjhrAd%#Kc?e16QdZNd&pen*wLdsXoDkzkkUbiFb%jK{H6yY6Ve<6 zq9BW_ijf`+c5lbDx~~jC=c}H6mtF86^*ZZrYu%aEj0UDcdYwzgQ8YB;?%hO|{yZe0bz_VKbGir9 zF+`3KD~c6PMIvn(qCjTnl;t)M6D!**qR%<0qC%O{>XS~gZ4853SZP6JEI$7t_E+jh zwLyLnSl|KnW6kjHtSp?bN>DJcLK4~Vbm-ev@hX0SJFZM@TzG5_ap&_vu}%zJStou* z1&~SIXsA%UfuC>|d|I}zZZtmaZ)te?Q)W}{k*9^K^f(lL$>VZt`hbW~3=Mw!+ZB|p zAkQt)-(UU0OAgr`987$Mq+i8!In@Nof@3VV&ot{fNE=d?3{FntBjh11Oq^7}(1KLU45sw!goGYWGOjPPzqvJp2Q&S(%Q0{4SZ0vl^yTDhh& zv#=2QVuQgsNEUJ!YNXFMNqe#D4=#+-uL~6+%!pmwg5>tkEX+7a+ci}&h6|E2Gh-7; zRo$~k0-`^@yBu-_BCld3V*;bSUGLkU)7O&go=FctaGDGYAxaHNpAC7|3(ySsI*uE9 zcCh2EM6||t>F?gXn?xy)!L=IN$+_bn*LJZ`T%C+DthA&2xo}nIh&hR;`OK1A4Nf?o zm?7}@iz7zl^pSm4KXmBE{7)ES5BwIrf$TW$SpPX;opJs*5c4d3r`ae(jVF4((!y6M znVO(w$z!Z2^e@H_rfIA{W4CV79TYZ*(4*r*yf_DBX0@wo= z0JVunj46mJ9Ko746&xhL?Hv!MRwO;)h%G1K3cP4LR5yx>i;LuDa~3R!dl?=~kKKk3 zlE-fV%R-@4B*8vPwtf~gGOpb7f6gFNHCFbC|8FIXzhWGXiiO2U@oaD93Sv}X%t*>a zAaeDe1!_3S(duAHxsAg<`LQ5tQyZua?G9l7kW8U>PeaV|ZPerP#E+C1cB|vVnK%`3 zzMycgCy|MPWA_4J{2rqJ?&^kCRvWPx;R_J=nLKo3jdw>ZL`$6{1x^wRytD86n#U$- zmH${Voipm!A=O$V%a6EBC?C5v`6;(GW@=Q~ho6D=aC+yHQ`xJITVCz-jXzKYg`s?D z{?r=oow=2N?->v?$lKQyWH45X)zX9d0P+0LlfUm{2LA|5G*TnLhY343bs$$+(p#=B zbOhcCL?sFZu`M-wd;4^0GTQO#{ji}L31lczyJHpO=|@_2f1U=F-#qO6Cnf}6^96Vb zgc^*TB235S>eN|s?8?eN6jlI)9Vsi!Rjfl);TeJ1A}SUZN%Oz81+KpqX+tSQwaaM!;8pb(y-p+7fAsGQ;v!e~=+`N+RH93Po2^ItR0FO+MI&8PCr)w&vVA*_ za+50-ak;C=_a|hQ?HJ52ogc{=p>vLpkGFXfy@kYLpL1!FcZ(45$;Wi84vBKA(qJNv z3)e5${%|@zadndal2`wiFL?Zm6&`p2J1xF#d@vz19=Tf*>Patmgf$2LFLD|qKBr&- z;xQ5z@7&phCuWgHpI}9`cIKWq1gfpI7dwWjK_(nE;+?lO#`&VEl9nXJ4HM^r$k63FX z%u)ibnwwc5&NkYukK$DhhHQ*n=J%zAaYW-SW`mC(r&jmyD=yA9(bauj`c_9JOEluk zmn@rV{~`C+b#D$E1gXw`$`!S9JX01Rf(%8RPhzZG5z`Ia$i#VnwCq4)a}%-7NbT?Lx-)?e;%3b?#v4A4+_0oDk!W z$tC?GO|x7sdzNOEaj~f=WCZ!g;_N34ic!VCQsL?|e;xv&$9kF@#OP_%c4_OQ|3CQe zX45vL!^!+p^!Ha*A*MmiXft9sOr3krXku=z5^2Qj z_f?@y=I72CDae&yaZy6Ljpfw&!LtZ7#0d=i1q6H1Y0TrHS(ck+!Er z((Sd{Y`-U7R+?yrhz#`d)qU%>WP|Hb`ToPBY?INprcX=)h7*bMbhb=P**DYDq&+^# z8ey`w(!@ML`v+c{8o8bh(&8mAlr?XkD0|t!oa|8OnyAGkJrzHqxlRkQ+9RvssI{#`Y%eU zDs0~bkZ3ck=(e+~)SFV0&wIEAnQ0hJ;?vX6KS^o130;M_6upgKljq)z*t2pPUMuw4_wo2O ziDo++FQe`N6QS#Tva;myPb(=HBCq`Tkqfz=XUOQccp+9qTnp$nwMPa31nH4`Apil! z83tt7%b~%zjeRn;^~1v^2^;!a}hR$?7-g_}_+qG+`5l?_oSI=Vj@@o&m`L`~us<(VXG4)z}?W;CiTryQw-p#UYhAG2uiF27o zgo+n|t+%a|t9{jCz6yxg+IR&ys+e<`K97C-4A=orR(pr1674Gm!2#MTKUAx!qf>=Z z%8h{mjMNCNkq82L)`1D@B0u|s{4Kw1sSOhizh5#`$di?NR9W9HQ#2Ht9f0|HcD?RT zss77X^t;qKGz%6D@!nG^77~hg-6ncRKK_Ea!|S8*q=31H6?%ct(SZ%fU|DL(N|pomh23~D*3C0lE3uvSwB}2?$-I%o0j*t`Fyg8O5QL+F@BJf|W;NI7=V6KSIJ^VDP zHtc>#=;{U+SJ$gouH4SYRZTw_3UvWbMKGRxCDSU57EY4LJ1jQaV;`038g5ISt=%X8-hd(@wQ~eIH#`3UYtTpz|Z!|442C{0Wa5oUZ!urYhW_ZB41oF3wwc z>pr|)dU+RUHYO*Htiv-OPK14_N&h&g#s(eY(zkoe@i-&dHa&?n6(SUsP?Iq|4utg*Fu3Q-mmmaf_%OPUt8r7wx1?+I&c~P4)j{ga zoNUm@@Nf@f52;LANDP>#y^Y*tjOIB12%gmRj5rYJtB7we7ZMU8Q*opqU;yis;_9*| z;q2lz$VH3{t~jxLetzBt(%MZ}s$f+4Hr2Ec)KhSoWWA?^@eIJ}PS$%!k^ro+rO)=3 zg2Hda)*@J#*!g%)2ze(GZy@oC+&*K@dxD;_QN$*aK!^a=KbfT=Q;6hRBIAQed<~F$ zJO;t;N)6tXiE}Hs!dj&{O5AVrU!zDDR+4AVV3?J_7vit}158aRC1kutJ;X!4E}D5I zUQcyRh5Pu|hpo(tc;U5jTWZpZ1Qq}}g7FfBDX;45Yw+n;6B5k(_+7-y=WK$w{-Eo( zfO&pPKTnW31}^cF!q;DWZlM6PX1BUMG016wQAR8(T!`IShFTj$Zv*jV7B)W((bClo zuBh-&?|xKog0=Tw9-s=6LBi!(1{rJKvwq%*=k3Ad+9)cTi-(#LD62C5XcFj*Ox%*# z4_vkjuq67Yaa*|WoUaM~lag&!*z10_i#%cq_u%=PH?)ZCTr=m^61v{{@z2OPO(s+p zV~p+~1@o_=?Te==S}J|IO!j5KOj#5dId=16sp zn;D93W$Xmw#hO(qewU!Lm_+I_-+`V(XdH4dC<}x>uxZXAcIjXhE)aSM*U-RnfI}Sl zxFmX;|b#igDX&+BbB?3+6Ri0?~r9UAc~9s z$y<@M1NxWh7W~DkvYXFC3ib#zc+nOUX8_T)7kvXq_yeeH(E5)aT(n(8n77$wu#X9V1XNA5GhcuOS2;E169#)0P2-uP46oy{~0S0SHbn` z1sAb^3Z(+8+Ubv9Y0SFwVPL>!wA3;av<*J$(V+%v+yj~hGClPbQ`l)%v?fji_hKs0lR#Kx#(jvoTpw1%OHU$I_Zo4_|?hzu^(ob(<>4y0)Sy94T(D|Wse=J z`=0$>GvA8yQF%nG`kw$q{}w3ybBOvUiXYbM(tVDRMTI-2WgIVdXLn`p#M` z5TxcURd!%v!F~!qmyP#@%D__(q+XFw02LwD?_yTmsI#^}c?(pfr3X-N)5XQB^}<*c z&;Tm?ecVRXDex=I!J?MEie#=8XBtdFUNu!!+;oJasGOCzY2&im_SRyQx1U3uR~|B5 z8r2;uZiGI!dLzj*&&19i1%gO zNjHfz6^cuKoXjK+H|?ey2Jk*On@#@UJiIg*Rp_l)YVUMAY1dFz?J8pE4wn^w;h zDMtt7JLc6(ZW0#x_(uW5PPE${A<}8X{xl&FxSrVJ5;OrEAccnUjpWclpl^~r*}Lbl zMHfMKhaMGO}gYfbGE8Pz;zuv~idiYxKwfm>ESof2E zu3zu%ZLO(xQCWEA*|7c|h1>n_Ei3j6=P3*Gr-u3LTIMk{U{g^#l&$>h$Fs6cMfqyp zxwU*v$-U#azxs-eXT@C-|1VwKIwhhRaM=(LO@`O69)$?J|{=5 zxoMS9RJdG(>%-Gwhs5XK_VTo>x6(<6`R2dEBxKqWAZf1ih8cHu`$7jU$??ob*x&=+lvuA3nSnwOGpAaPBq6myWeYizlo@MUH>MoU@e<8$s~O-P5L@m1VKf zeEGvpnXSi{-&ip(sCRaBT-8RbN22I(Fi?1eFX;yUCi#-QPN7(Y?wH$P-`10{r8o+p zG>Y^2)%Jfs8DzM0Uwiy|Tcb;7CnQ=uBx~52cgUpfdrQ~wgalcMvNPc#tuNmkdXJJr z+PdcK;OiAlO)^6lD%~Xw9uz+8yCPXKKXLFL-?Jn8zMOqAaC#TJ#)CZRouAapF{y)8 zq161zTD}7}&SWfTOgXU!rAf~Z@5(+;5t$>Gql1g1g}!L>fM05bc3VTOj#IMm2Uq<+ zr-JqO?g7Uij_ZeS-pu3Xi^>;tA_%NnQP03~@vb-r)ns{%E%Lfcr?wN6Y;Pl0k9DK%$)XF{hW$F?y9)~!^Va zJV5}RMq{t{e4+E)g5A&U)=pI!w4JA4DydIhLWELVPl(7IsiGfnjtB`$bts6mpk%=u zm_BK-ttGvGQ(2&I*Sri-fK*1!xn2+g90a&?u$u{FnJ$ zQkdN(-6_Z1RveoeWg|b|$pf<%fG|sLXof9u>iqfXD8E`Svdw#3pef*y`d`b}eAO={ zMjQ$cpAtS-MHW41NO-b)Hk%6paYfpAuf_Xi52FT3%`o_x8f!tJL(E!I0&1iP$@!5h zo&}VDynBg)A&_RxMy@tvEfkq8jfke}E>vgmc8N)k_bE}+|J~4pJ>@|~_^aIx(+$cw z5O)YwEV0mNJAtl2YT`i;_9%piiQgFb6xz)+hLj2P3W_^WRsu z->zRzMY7Wki`&=Pp?S#PKZ%}85=@0&egfNB2`afIOWU@A!kN+8D)T*sS^xUUU*A9u zC4L6XUJZQ2v5Voxs*f(P#bQDEl4|OQ3o2aV8wKxgQs7o`}86KD`paR3~FX2qM5hPsq*RDGe+(GBxO>{|nu%dw4!3_Rp&LzUW=rUl1`oz3<%DSku=e!a zOK0ucz1v|$%3Up}>k^@`3G3>383%REZf#~qdxm3SVf}hnZi|BTCL?P+3p+ZT?+u{1 z_{%FZeSb5aK`Jh`jjXAAyL@Ya9sXN4rle_(7Jpwl3eNIW2I(}S-2Nm%~vH^+y>K#-$3>&@9-h6CmwP||Lvwh5Z|j&DD0U41rl4!^8b8z*m*}GqN7u| z^Z#9R=$}SygHN5ucc2dYkMK~e^mCC)QD_GeC}Q6j2Lzx=3%42%GXH=DaZsPIJ4U*eVk< zGonD$e;J(8pSXzgOTvPX$LkxTQ!*FA1yT^Wuazk$SJs}CmhR)j@-Lj3d$BF{2^~rj}Eg7c2f`nbn%VY#S87 z%|~z&xGRW*NhC{Qb1DdNgs0lIPgt}lc+ZMJ))AgB(P3ReTQ;w{|9P-ec*IgjSzTdx z!mH!!`l6Bp*Jl^H7PSkHq&3cOE;sPV4i5HQx-vRTAz{S_clG{s#X{W*53fHy^$-4L zM%_cu0~&QxMEs&QdSyD_uP-e_3&3Ewf5w<=iQ5GIgI`oMsT2YoH?Uwz<%+JMalrLaFXGv&17tFx!AscYj5#K7}MRhJ$!lA%#r(C^)r zerENYBpKS;+Qp7O3N~qB3qW_l>Cpz$Ld+5i*z*oXMVXz+_Jc+My|Dv>@v>q{|3r8v zY@QzhCIDOmRA~uCN#hnQ9w-?#)INXy9M*bn7(=ijmpywX0ZHFJ9NT0|gfaw2tUWlq zjrW~Z7U9sAkEzT+xgUMCLX`OpCT#>B6h5FIfcWnBfHs2>2!a)6-=B!eL@ zCvHCY+wv_RK_Ny_)=>+)wldM)r}Tq5FT3hxuyO9!f)v}}!`6x^WCu|}EnlL`;AIFq zT6CVx)}b(>_Sm$sZyx+x1`&|s*nLJI0mXCV}Zm2JkxHJU$}Erl4s5G?)Z z=xF;~(>Wuj$NG(sL3Jd=G)4qCXxEcA~r=gWVkVh(_&u=C)tH+TNHVFW&?I2peauS-raDIxqG4nO=)K z5^k0#&wf4S`PC9-#n-&g(Ji~m{`3aEEz8&zm(4T2Wy`)v=6n@huSv(Qu==#Rm*XCq%5 z8XDM-cdl|h)6)}f#p|#$-U_HiM;BodHw3vf z-v{Gum>rAz#FAa?(D-m@D-VS=XFtHZs8gjiV?>+K1uy{68kvmt5xOW7k-i;7cqa`` zT*CmHf{i+I9s2vE`t|?u@qBW^HsE zNM#Ag;>BB0QBgNhdj!x}2dl;&Dw23DzeR8uSXC*U8^Md=i>A&CJ^<7t2fxS(&#eO) z9?%cZ%bq(_%Y%o4D?)|1n0Q330Bo9_$?sh=Q5RmOxP%Xou7u+zUG7L;J)qFG$Xvpx z8xuRhM(|5fSWRy(f~O7jFPi2;F52b(;{#pPOS&KyK8XEc!h{KTA#s8wU6_1^c!jLg zs{B~(0E@@tPD9urR~ZT3Y!p(Vw$b3!Uz~l#pMUGPzlK`0~;$%#e)R-Cd9zwZ-@k) z{+~PIAm=~GH-Qpkge=iW&N^fkRI^&a43Mbpu*5a8U5FlnW82E=?hKOzYEWh8t+UG zhbt~wQdnL(n};M?s6QZUvc3rx53l4z$f@M=ujD7-GU;;dMOV}wrgF4Gb80qg3-rAV zhma!;aS_jtH1}LhhLGyM-<@l$XyXi{kTiQEHJ-H$n;mh?jd*u$H&%a~l%6=4BAPU36 zJ1AniE5$&3sDVN20B?9k2@m2wxha5S%@z#T};rxh^5##TNOS0iJ`VfxTF6oJ| z5m=}9pAV^@qpWO>8<7Xgyzbn(#Z(^ccR_6@bUaWO4_5t9ypQ}@!O#N`}fm~2x+XStO9JViUB}+ zOJe|o6g>yhd2lhJ2os!@(wJ|Cel2_zL}ClNJM2whygHCfXeXl~$P}`+ zt8Kf_M??@T-%Q(1#ER5*BNQQkrTQ3fkOx*lm9MI*Dq$dq=0}GpG*4E$YT!SpsjoQmEIunmTkihMbR8Ns^zN&PbmC)tVqE_j-*u zfb^_`AK!%I6o!>jn6#fz{=0pCzll+XON^s)CnU6*1CKo@z`-mYMhn0z^yZ&x!5$D)mA)XLc0cLA*qaM4~Ro|z(sR@okg2;gTW)DRA zE)35k@du{7(6BKQ3((#Khxf5vNw9wthm!~-kdt3W|CgHqYBX;yb|F~3IC1#XUJ*`X z7_pq13L;YRZIQ?od(zxIycMiz63fjpUa?{#gIo>__%85AfV|DySjA5^B>e+ZIT)UA zpU$Eq+0-=hPLfESsy<(xYAwcVVAsk>O2RgVRbx&;m#`ie(a!;Fh;xrSa^i$MBLq_s zGM+Vqle9b%wp(yIOZfCcz@+BFW|W7HVt+SmZr6b1g%c>#i9o}eN(!rn`_DIFrWdvW zM>Nx8X7cj#_JTJOTXX^o!@A<7Jkh!YqD&OD<;9FCPkrK79&7r(U1`_)p&1SW#yNKf}qpkTz<3$Kd(#>2c9A zr2p@mtfXiLPbMhF3) z^tqmUAnKC%oG_LSj;Rm(myS_}j;=6N!dwDQhKylzNJjKluf>5O5!X@o&)ZZjOx|ov zd4TyfrrCAdHlyM#6su#LGaA3-4~B>hfsoNWT@ky;74ud~^azaB@UrALofBukuwZmC z<%ML;`Ph|etIQV&ZMd`ux`0t~jK{=-9MWF+TkVw};UB|v%w_!HAt7bOJ|Ad?!|HYG zU{Ltk@-7FB;GQq3Z$B03}1FEKD)T{q#RzzM51 zi`NWGL2RQqY4H&(PJuN<58Gre?nJcu_REVii1CE{noyAio-g@9Sq-Vx9-4};qyitg z7xWk00#sbaadQ)tlSsU18Sc~ceuA(QOA-+$wDsTh6I8X<@|gezk09Era;Qp$pg--@ zMpJ{tJRBZZP}8x;i9i_((snMIZ9&ot%mOg2#Fg;w5gyt1={~LI77G@s|(#_V4;gZ z$o3ouBTW@2ASZV4m&K>Bcaj8L*Dx$IrL)aoa)EX!)I*xl*?=!" + "
" ] }, "metadata": {}, @@ -181,28 +181,33 @@ "\n", "\n", " # Finally, let's plot the solutions\n", - " f, axarr = plt.subplots(4, sharex=True,figsize=(20, 30))\n", + " f, axarr = plt.subplots(4, sharex=True,figsize=(9, 9))\n", " #f.figure(figsize=(20, 6), dpi=500)\n", "\n", - " axarr[0].plot(values['x'], values['p'], linewidth=1.5, color='b')\n", - " axarr[0].plot(xLBM,rhoLBM*temperatureLBM,'.',label='P')\n", - " axarr[0].set_ylabel('pressure')\n", + " markEveryNumber=20\n", + "\n", + " axarr[0].plot(values['x'], values['p'], linewidth=1.5, color='k')\n", + " axarr[0].plot(xLBM,rhoLBM*temperatureLBM,'.',label='P',color='r',markevery=markEveryNumber)\n", + " axarr[0].set_ylabel('Pressure')\n", "\n", - " axarr[1].plot(values['x'], values['rho'], linewidth=1.5, color='r')\n", - " axarr[1].plot(xLBM,rhoLBM,'.',label='rho')\n", - " axarr[1].set_ylabel('density')\n", + " axarr[1].plot(values['x'], values['rho'], linewidth=1.5, color='k')\n", + " axarr[1].plot(xLBM,rhoLBM,'.',label='rho',color='r',markevery=markEveryNumber)\n", + " axarr[1].set_ylabel('Density')\n", "\n", - " axarr[2].plot(values['x'], values['u'], linewidth=1.5, color='g')\n", - " axarr[2].plot(xLBM,velocityLBM,'.',label='ux')\n", - " axarr[2].set_ylabel('velocity')\n", + " axarr[2].plot(values['x'], values['u'], linewidth=1.5, color='k')\n", + " axarr[2].plot(xLBM,velocityLBM,'.',label='ux',color='r',markevery=markEveryNumber)\n", + " axarr[2].set_ylabel('Velocity')\n", "\n", - " axarr[3].plot(values['x'], values['p']/values['rho'], linewidth=1.5, color='m')\n", - " axarr[3].plot(xLBM,temperatureLBM,'.',label='T')\n", - " axarr[3].set_ylabel('temperature')\n", + " axarr[3].plot(values['x'], values['p']/values['rho'], linewidth=1.5, color='k')\n", + " axarr[3].plot(xLBM,temperatureLBM,'.',label='T',color='r',markevery=markEveryNumber)\n", + " axarr[3].set_ylabel('Temperature')\n", "\n", + " plt.xlabel(\"x\")\n", + " #plt.legend(loc='best')\n", "\n", - " plt.suptitle('Shocktube results at t={0}\\ndust fraction = {1}, gamma={2}'\\\n", - " .format(t, dustFrac, gamma))\n", + " #plt.suptitle('Shocktube results at t={0}\\ndust fraction = {1}, gamma={2}'\\\n", + " # .format(t, dustFrac, gamma))\n", + " plt.suptitle('Shocktube results at t='+str(t)+r', $\\gamma$'+'='+str(gamma))\n", " plt.show()\n", "\n", "sodAnalytical(maxTime)"