From 9ee07718a66e064ecf4ca2cd941d707aaa2cd9ff Mon Sep 17 00:00:00 2001 From: Kate Hedstrom Date: Wed, 1 Sep 2021 17:18:17 -0800 Subject: [PATCH] More on background mixing. --- docs/images/background_varying.png | Bin 0 -> 51342 bytes .../vertical/_Internal_tides.dox | 10 ++-- .../vertical/_V_diffusivity.dox | 47 +++++++++++++++--- 3 files changed, 44 insertions(+), 13 deletions(-) create mode 100644 docs/images/background_varying.png diff --git a/docs/images/background_varying.png b/docs/images/background_varying.png new file mode 100644 index 0000000000000000000000000000000000000000..44a65175a049779d84a89e99bc71a90aee71e194 GIT binary patch literal 51342 zcmXtA2RxVU_x@(@ot>2x5<L*(00(d3t~U z&#Sl3M?Bx>x$kqIbDitD&J(4jp-6~Fi-$s?2$hxObWkWXL-<3Cg9YEITF*;{zpyM- z73EMD$bWgQCCGPhU6hO-qfnG=$e&BfI;>mpO>B2%b$RS%T(nE*RH0tNa2Y0)vYfP@ z*Y~w-Z$G{Lw--_?AKILaaB*=>2vkXJNVQ6bR57Xa>h!1nxM`2-@2#LU{0?KuGkEmy zVP(^#@%6#$ToVN;l?}#nn9T1fdr(G70vQ=hMd_BpA{pa=Zm=8o%o5 z>+@f}#4@g=r0?&v7wlSgoBuUh*mV56yynrAFukN7bud$2OA9^SEhjk(3+9iKCu^Zr z^b$T45yXsrd#mHy101<`D)p^?e#mgg#0h4?34Sk1&HQa9MNW`5Y`5ih_H@Y%jXKQROJws&^sCI`*i_T< z+S&rN&KPP-zA4MfqHnaFIEn9n*2RhFKJ&>|GACs2$+A-^PNEG_!dEFKbEv1q3k#;b z8ndz88LLtptWu2D%@^|fKx=1ht(DiI-seuA1eqc7l!?+_@B!00!;VL~*VO<@CPiuhIdAJs5C^?5lhdH15*J=yv9;|0x z;_=_V8&~thBqfoRPW+aBEU~mOz3kSx>dqo1#Yhx^9?TSI^MGh(b~4z4`tU5^^Tb4_ zWldgbDGmy2xc=Eie0VrkU1MWmQ4t1v@&m1eE9kBt-494Z;)_+&7hXJk-Fa))FT_h+ z@)rRg3JHw#dmc0l0f!bCA=2+k(Y>yVGZ0hqnns#9w_esLd939_O (?#UmQ`7o-5 z^~TwcQH0nb54^dK#1GcKh3t+E8r7JcQmhr3tYFNPg+ z!fk1VmTl~D_3n@lp-`O}T{NPp!G^?|>>CO<9;M>$_)X3lSWVcna0!^6FQy=pQh zEsgTQ!HhZmjieKu_J(Dx?E06F!&CjHJ4xt8o-YP0h0K}nm_PaHP5mLL=q(8j_g^+w za~j8XMN=P@dJG}MkaAsh6|5nNY#svaT?yWcD? z{POz8Wai_zUdtIg!eSlV_(fK2vM4!{g#O!Es=~8!1C#u;amW)g*DhhHns_5k(eK<)(g1{!p2GY4RIa+c|np?zFUH zC*MiDkMHVzlQM?wG~RL3z%n=KKE5;wnd7>76VbL8JZ)lpmFJ{~0+;wQ=YcYP`XBC==EHOu>FM(L|p za0FC-EzLxgH~#R!&%EWVvBg;=mJ646oN6FUO-+3+?!A*os5jl_e;b0yb|iwj?NrCLTX{F^jQAt_NrjsJ-U zsv}GEao&VOz7{7|Z|ptEQbutv65qdTof&pHQG&-I>ur;E*?Xh;p-(SM^V^-9{!h-Z$ zxnKJ}NdA02Z2KYOCbmj(?YIj{n#5cJ~fg9<>x4}O1rjg#= z6>o+_B@DPg;i!AFrPZ5nTM183(F+yge!twT%7S#3YclJ)fb2&s`D#s2Vwc3YgD%wE7+wFJQqwjBKMEiEnmHTsrb z%pHf2kTBb)O71UYp9lsYoZXd?{sQ?JtS5gTz0DiOtaa7rygw%HcgV-!zo!k=ay}TJ zwxAfPnoXyxsmwFyhJ%ccp~|&NG>!*Brdb4O0c?%!JH;}i`h=8Kr7h0I z#+PZ;&Y~@-KY#gRIgrE#Cx#;Jx(QZD=QF+x6$Z__T$gP>JinB!F?pcJArs%Z+XAJp z8p0yKx?0QYXNH|v3p-Bm7UV*l%a<>AaaXmB&DXb82b}p>^hA<&!Ah}RF*qg1kS?vN zqJk?iPrD6W51!rJ+&U^$Ou$QKS0@WNU1Y>hR4VFkZ{&a&zKdewF>4Wd)Smv{Ke#im z+Po9(QXuLMg!PNLeV;N5DC9mw^dd}#mHP8gRH^(V?-KY>3Zpj)ctva65kyzIx{0M?|GaJ{nlybJxAZpY>E`kN zAJ>H*T=A_AbUtc=D;~FF=|#y@(|9_+RGMrrr{DMZ6Td$EL8F58<8BMZZ3U8<;rEjj zMj>%@w=Bnsl%Mliput@Q9v}NSFahoxLd+ButKHGH_fZ%Kq`czD>}!J+GOnQpV=9O+esz!h%*R85*8S|ap) ziwRNyPCVvbmA{q|q*bT+L>(np{q!P0En7@zxt9vlz{rRO42`I_NL1KTYBEBS9ky>a-W`i%jB{F7IM=wkK&}jqn8Sx#|pXB9YMV5hT^<_Jxss6 z7XaZ5hef0Jt}TLMw760Y-@35@U<={M@sc4kUO36p;Bato;CuAF&M9>tZr641m-d5n z7Whz&67rxsD7adLPFc_4_F}#s7o`5{Hl7fTOkznlqZF2#Z^5tA6Gc!grZ{i_Vyiaf z;pFC)y?5_k!mb7zBRUGnicin}+&8Djns=dgZa)*w@%~98c3we5I)X2q`$r zgoOedn?J+Q;HYxRn_t6jsNHMr^wbZDDmW>~^`b9xSNk2gIex9aBqSt+yYQNultaqk z$j*L?1#2OI4C4}qW{w;WrR?3imr)&-HRh*>+xAXQeNp5bSn`-CfV5yNtP-t@<3`;)T8Iu@Kez-XjI&S2S+Io0}Hh;hoEF`K_(3x`|zU1OP=w zKYe2AO)VmS^DXw~h&Z`o0)x89jcknw3E#i&t3s&2nPyMbO^(Q*=cO+65u0r%M3ga@ z-Qk4EW|`vN0x2qNL|4R8v~qdU$=b;fwA^*E{D-tR%J9X{Px-R<0$%6mU)uIEH9i-~ zJ}iIzEJhe=Q@t_2)3k`+b}v=-eUdbJ5qgCuPlOZ_7zp8Lw5Cu$i;Cu~9vmFZoBUHe z*>=#n2_SJU)45Grgg))~mUk@@f@xu@Bvtr^eR z&c$u>g`TLAd-9SKc1y1JC~nHSPuc;DK;po%Co*)P>wV%?Z5^FptNwUc#C;z(Sk}>5 zk2sreDcr^je=@c@xF6R@$0=N8*T{+_q!Tuu?+SyHv99Ge{%3yy5>yXdpxDoF;SMvE z_|6C*<=?*@cuX1yxD8(S^hA*Zk}0Ou-I=eE62SZzt&h`n>{_-CklxfS%I*i zo7+9&dwTAdl9HlS^=;h)s?b9F#Yy}2pZgcIn3$NjZ(qZT3Oh|l_A@@EpS+E3m#=XPww{r{$*qY+o7#{22 zQ?*d~`C3{P1fy4p7{qbr^cGCK^NqtTw2Y0T)iZ^aUE`s*`Su}SHo^rkjmqK!R#g8xK3Fp6ArbcpNM{?ahQAz<$OU5ZxV|@dJiTZzVBYu^Lb(LYG~-GIShR>XLwI;4Rc z8H<-nkhdNo`KfWngYL6~kYO5Tzw2|+$F^Nf$_y1|e9k|_{kA%>4Vyi#oB8kGT|F8| z=J>bRE#A9%`sL(U*wQ2}=SM4xfJE**etZ+UD;y|uTgR6TOoJ8{EGtYJi5UI1F8%Oc zxuTq>o+D24^st)|xvah53yixwK8XtP^ifa#d_f}RaI2em`}|}x#DW@vV=hsVp|gI{ z4FAoy>5^F|L0VS;mF(h3=@5r%Syv>x7kO1N{!Ym_f=gblUi+Kf)g^-}ncmo0!G&^+ zd-SuhH-?JC@oQpx4ldox56cR7AyT)x@Y!l8@WSBkwM|W@%ebFeP*YsKJpauu7akq< z$ZPOGJwpIf_oZ4VoJWGjsq1Xluk-MzJzwQ`n&5S*KaMUUmD>OTPC9um%F0o!hXT^# z4Cx=Euw)C~oG4T%Y2Ck1+-jE4)zwARr>Hgsjo@<aGbuo3Ir;N(<#8^e2d5id9psFRc$6*o=FZoX>cO8CN^ozP|D!9fFNLXBOw?}niC z-jW&Vi8%xxv>)_D*lQlX5erKm%~xpSu=e-&=e*6hun`&}i^Z-ab3jIkH=t)ymk?*? zGw!`&C8ZyaIG?uNLXSJWZEwezcwJCPD8I6j@aX6Wfj{i*V%bRQ?guQeuEzo=gS8rK zHP)23E=k_N4n0`h>{YoHM@)hBrAazh`5`VjVKt?tvxB99WI!F9JUkGxs$1r`!KhX< zVmSY7f7lY?PWa_8(x#JBe!CiArZ3|Th;hHUdOon2;q#|8DW+}wh)3Fbp>5uq;mMOH zk#a1Yyu9I~_G^6kB$)cTy5zvk2+>QE^k|U+2a%Q26y3jQ@92n>=|K%P-)ccx!ZnZE zGWjIN#>S=Fkn|y$VI-=#O+?Af9!r_C3>PPu zgtXPrtcr?f2&e}s=}-i)EG;iz3S_OekIH1ikrAdTW{1Dp#)TvYb$u(l{t z7!!0e-MB#vA2*-?ffjiYB5o+VYI)&2-2G$2gEgEIkcCqt9bQ;LH(5SstzlYfPx`OU z0>=#q0J;T3j*Tyx4cY?_(%^1qJWA~aqX8AOP!phu)EDXM>L%Cfl7`Z^&l1B|NGB?$ zH|dj|96Bp6X)O=;L^CTuVg2t@Q#B7Kd8y+RpgtJ#gliWnDCZoisSzogFPC<>HP18w zwA$XEysvsi@q>Y!8CC&za3>S26FrflOw#C*kr_~5uz26UeY}g zxL-HPs8lT-oNO<|+*oLsazaOyv0&oz6sY#Ri_|UA2%B(7mAeb)T=L@d@s2RbJpdOZ zcwt&bMz~=zqxrG%k;AqjqxNzpCUj^)cU1ThSVU;ZLz(sXqbOiK*n>nRxmRU)OjmK+ z|J`iuH%G&Ihr3|##mos2aDKcBe87WdC^M!OznCNf^*64-Ik?JE?(ZWPy;S>YvD`_ z;mpy`pE25JHT|QRA}i`1c9b<~!}^AFesZRBd(T3Kk!TK~PK8q1*!`H6GHiPjjvUk> zIAuZ%kfL#Uf|-&aPEm92VwX^1VPPe)avV9=QdPo&gVE#)wV`ER-9?Bpz?u!7)5u-8^r{NKATUl8-Clc>gLHZvp$zo(4)@=dzX2_GqfXmZ7Ljk3t-zegJ%pr0BDepurepPB6AZQ-#FZ+Ioe)UmG13=ZRImNRS%CiI}e%TFkEo*==#+ za0?)Yl=Vm@x)!&Y$f7#tU4VS`Y~JlN=j(8&FKBJe7KdYN-8wE!Pu{C83e?`pQEY5N z%76qm=}or?8D>BxIKil5ohW-jl{i^43%sg@uC`MP5BthRqPAN10!opEVu*%mf;ygO32t34_lsX48P9 z<)k>(#xs_n)IvEHQa=FKhj1i1_+pnPbT}!_USNk(n?^QYk`J1=d!H2J63Dy!l49WBkO{-7qa|#z~QrA{cNk7Ap=5cZS##E;GMvdEmMp1V2wt zpOiK)CMPo#jfUF(;K75g$AHZ-Q3LWSBm#LCE zA>*WBq$zp2zVr|BApOa-@EYgnBIttqrz)Fzq@D{scH-K%)(I|2p_8uU{_Em=ny|g!Ut9^ zO79KESmjgmCOfpKVb{rtVD#wkSQ#_0ijf;aqEslXQB+otC%{4n>CIt!6n`RyetZn~+xKJ_#ndzf+(D^Xy{+6JIe;25 zU@gN+rc8nm?B7EDI9mN~<)6BW6D4N3kH+8HqW3?aK`@x*a~L?CG%`P|e=!j+?BHcp4G%u>5s*z(Lt&$5yl_zU_4Uxv7YkCl z(PKw+2fg_y#%)!P^Dpunr`8dM4Hz9Zi@U(t0*bP&*dXKgA|x5vO*$BGsQ>#hk~5Vr zLv2okm-(*j%UDVf8JuZuInk%_-s4(6jOkx|F80Jhp*dj3k9BZ;pm=p6mvuptg;-5@ogH+m{mvQ+YEC+k)e@McTVoB}5 z!$3>Y`hDHsw2}M^pn2KZVvS0)#cQs(wr= z(?C_;w;0g)rUp;@Tose6(StpnH1j2F_S{UvD&E;~n+E<0Ht-eekA>}~J+W>+ucPJk znt-!oq8wHr`xE+Ip_8!sR%=t2qG;_$#`5>XMI_do4^m3fvu6z@2N}JVK?;RIQ|oD# zXK-A@1w>ornxJDh_1~X>hUmG_LrQ3?aT5ANR&t**BuJf$1ng$c&S6U(Uee>bn~;7J zL;4Oh#n_B~fz6w3bWp@aHd}V{M(yFFfoixv2kKmZvZiE|iPKw56qy`7>U1~v0s%9? z^Jg;^$c=n0)h;ZxBRh+p-Go*;e8wrgf&W0ECXE5er)y44*pyXAm znoW7dcR$?pI-HL{7$#uuIiyXaW}hyn8Z1z?AIhh<7dG<^sH$_}Sh#qIu|)rKGJP z1Rvb+*L9E6o%Hq;+VE z*G_X1z#Y|7>$s_{x_3`EhLehw46aezcU;xkzml^u3h6dy#W%G;dQBJT_mc5Whk6#k z)GNDBex)O+kELN@+pA=gw13e(C1pr7sF#lbav7+*H%+~lQ=d!v`*b<(p*_4p1Up=s8k{n58vIxREae9;$?-miwh0TE9q*4l zlO`r6aAZIL;K=n021s(Pa18)Vms?iF1RbWDKEnk9|0W9iJ`VE9TiX^>i&PPpxg8{5y z%-EO?ay*;n!|8l@I1iu~fakvUM2?FB53}w3_~E;u$Q1wYfA&?_oORSyd}O4hFNwq_ zr{N%-HYS8Ma9OakPaI}}as!%o?b{nj9?IqeVd2krU%h@!1Tc;$qEHqO>!=+#AVhly zkz>4ABys#3S=OPcukY=EqoMl+pzee4_AnTQuCd%9kU1=Jak}felI>Xlgwx>oT?2!N z^!ouc2Q!|GAVmQo^JzILmkDQ1gvhyN8w+759^mq&dUAcsyfYLkYof(wcf-lq8TjJR z87IwOBdwzP78dzS(l5`U(hD?@&+>D@VvFq5F0J{tnT_!1s~nQ0bxnnkZ0@H#T@WITUZGmM?P{ei4}HD`v^$ zI2bNFOUkiyw_Y23E7v8|0iu|9qh(hEQHKm$yj-z2QzAKk_H=k0k-co{Nl#CYK+p(? z@Dcrst8R8w-F!L)dT*(2Ty{%9i{g^(sJJQc4^!sLO#xUWq zB!kdmU;mxAH`&j0v;6&F)U}s>e{94V6hX^?z7A+@wh=9^QSJNp6aN177$|fXD+pQf zK3#Wy+C5gR>N<42(R}~+1mFeSpv^K~y`5T{9npB2n^*zu?TSxkmQ-4r*GpElmgSs= zH=whIzTt!I%{`y92Qkhuy=XxHz(vbub?#rqVaZ&}ZnFG5({uC6QfNm~3`Lks5W=xL z3R`XdXidTaQUfev%enZ&CjkhZdrof&eStfWT+Qw*A z^PJ|WSHo?DcAly0N@ zkmwcxSM~;D;%@;#$MNC(_E2Iy2-t#xg6{ABa^s)@O@q+2z-8`VyPP_+78*H?>!i^`j_k>}>1ICHgq?rHLh?{DWA%RU0{Rl< zF!h`s?Lxy|+^G<2wd2!j|LBptID;mR9F_bu6z&*B25$_hcb7WMGKY-6pNLh*P`gJ03O8go9x!56mmOKZm=U zTo*fPd4RwsVw6NyU=D!^ngKiHAF@P~1uV#En(cbpE)!q@swCU#%O(Q+1^^y-vMUZQ z|MeoI`Nb0sHLwUg#&tMIPTcd6;Gx8a=!WtMH@&aZ5wEXU_Z%t85>-rMemWxp7EA_jr1(->2w1k1N zsEuEub|5n>!5!J9mhcfomKh{R^YNGJa5=#I2o!Sxn)nGPJ_}v~0wv4!E0e?xckn@3 z!I=&?SqJ?xDUCuj{?D8|ImV^l-do*(&A*5W=kkk z$?+fB)y54D1Q=M{9ghxr6qp5IGub)k^;~rZ3e|WYG2w z&!2_u9@hLFwnn{z>vt)=fvP}x@{>0oUa<_UIc@S^y|dFlAQOirA+rQjs}?7PPSWW3 zDxIYim1k>1eZ6)!uru`wucYMkHSl@7E3|8WDKjYPq)x7?@kKF|FUgFRwMo5ok7L5o ziR-Z*F{-WpJwzmvvjWZq(wN>NE6fn-)kBI9c`SuB4oK>-u{uPJKl;oDPLhUgfzJgd z$r>Fr>dMON`CT!i_9?<;i5+5uiK)Qo)}JLY=EkpZ+?D9~X!hr~gg?!ld%RT8tVKg% z4YV-=hagy>h`N%geOM1#|MTel+CMfMUenpPXrAHSuXo?GL25`f=-j5h>-A&L8NeLl z_)UPvu=)-y`S_J2+1@U<&$7%u>^=NlE&bxe(tVuoI>ZpXqjSaF02uofk>x~c`BJm*rP(TGkM5t>;M(_&X@Z6TaZ>x1 z90^ve+V*7HHAfj83Hv-_rLpxA+EPpZk zRKu1*&-xuxKOmqj#SRcqFa|J6vS|O`l0yIX?2A>o9+wDxd}#1U3MGiNNW2l?2C+rw zy7Q9a1v7agS9$obJVi9)f4?sZ7Xq2>)#c|O*<|=5>+As5aAaE!m|DBlk{t)VraX`P@^YQ!{sptgzPVqMX z3`-0gIG`?p`*mvBR><=F;hmKCsOefQ2$+ndzg@g=k(}*hvJfxtpO8UU;^ZMuVp8g&@;DAA=@MVwCE7#Pap!# z>4{+?x4SqPr>Ct=+IF&W`_Qi~R3uXE?rRRH)3G3^{W~U|=1QWrww7p0(i@^Gui{zy zx*}71qViP%Ridm}?2Od|e$fK?*R^gW1%< zz;=LXvBO@_wjwu*%?l9Qe}NP^9VU*0fdbyMX~yZ7yV&#LpRwN6JuG9uyn_m~Pz7fN z;j}l#0Kvj`z}6IwKL2trwBy3N8kAFji{Nx8Z;*zlMNI50D-dLGs;j>_>S#Py0}kmw z)r=LF2l!PanKMOZ^{jhbrN9FUN2ct98ke{pCj*Z&u?U*G}!R3Oe|12f9oKvF%m1-Z6#G$LYp>#I!HgZi+ z0Emit7eg|`jAO}E70mL9Xb3=+i8~}cfm3(*iS+7Z3+yDkKd4dZadBWf>rHZy6DsIKSG+~CHBMKao^!k^FK3(J_28+oCgQJ?Yc%94fYSx!z>OQ4 z0)iQblG?zghWz~bGnR%Q5Y1~xQWWC+1zRi;REC9Z5$F*Av5}w#CAAo=HfX|`s*?o) zQBO|~)H9{#MnD{hkR24>bm;NcH0e;eQCs=7iu0RqfJSI(Y01ox*wlng13NCa_j)fM z?ps2&Mt?-2LK|{wL|DbNFVCSugI72ll-yvfCgyNzEBI~r_$&>(;Eur_M{NoteQ*Yv zYtLTAv1EEK_<#%nHX)r%xYx2Z(zToA(6}Sy(;(M-X+V)g=jE?E4#t7((>w!vIJS0O zw;=vsy#=)#6ehMk2{L1VBe%q|HMoesFbYR{kP-z+3#I0X2(lNXB)cqEtiA+BVsMkONvhiA6sWQ4 zT*XRveEngYm3qcg+1J57qux^qx!j3Nv^wkph3LI&x4q+AFS`An}(P!4Zdkl0fn+i{JBEEjDcR z5dtQ#Xfe9ENhBV<-}P|l{&_IM%)(=U`~cSjoP=yMgWr}k;zLd%xfKw=e2DQP zx%h=#02_$*i8S-Y#l>DgY`e^LTtYw+I07#JPlyEw7MrFUS*IxpnF=5TpU>W}fFq}W z4ijX&1TbVYD6W?8KUHEstY*Wa81aPO+6(46kcMsH;;j9G!#f11ZGMX`46w6a26mAf z+`SQCn{71+nV7ve_w?Ht;P3+Ff#3N@OX%d&y+TF{;CP!h-uGUj! zuSG?2gYCVpSt+~#F<8=9zk~qC4d||Go)5l&IM_2hjge!4+&d@Y)L~m{DWg88$Fzxz zrhwnta(kh7j4V5g=BDjjIGs}ms9upSBGcvBm|$9CKn4roUs%0LymGKdc0pj_;(fbB z<~ZASg*Tm>xpV`Ji}ElKkj7(#Ryg`~z6I=bA-QL3cF43s+v(1b8DY*$eJiIx+1K4q z6o!yTfh8vB&a0@Xm~AsLji-2Q0l9y#WSzWTO%rWqoQ7bw5_@Q6*8+r^u@44J=As8mlK)6LJ1-&L+Uk8$=JQNZyhr(K0iOFDPKqGcX8Q6oD+qi>b>>wpwka z&0El;BWjtr;{Jre(gl|XkmsjQfnZaEQXSB2psK10CIC1PNY&CWFS{Vc7~K@+Ba339 zzV-)D$LTA{G>N2{;v#9f2^L*J!U@}vnTm?8mU%MO`~YCoCTgt}cknXSQrm;_HY z+uKbxJL?{dyTGYEJUi!AeF0ux;OB4|VSGwPQKgLt!WcF#MC58Xkz_70&`Gd|NrUox z5$|ci*a*@+@*KcbwSjr6Kq~77G6D!5q}=LPsa*Gq0g8clIAm;WOmeD^0rPDFtq>jn z!9+QEH>a0+Iaxsta<119pv?9AD-ECo8yg!D@ZYBA6fA8fmOLS8m0!3(R0*(1Qi8v+ zFXfspS*YDVpV#^*(QBZf{f<{lz$P+dmJ35-z;2Zc30b*DaT?K+%L0Y-X>1HFcO(Kf z>6iPvfJ_Q7o`TL+0PeR=bH3#y38XuBWOU)$Tx#<_Va?4-5gW*a>;D&bLR4d<5uVHzHx0O^!)b*@uA~AUMQNx4Qq?;=#A7j2`9Xx9u;X z%3lBf;oE_U4ZlLQz!o1 zK?_PAK-MOh%>Wn~iYO5KX~jq_}qQ$n-d2m_)`S4|%L@`=nKAP^dR@E1qF?~MP!43y+ea5AGN zK~8x^#hO6-<$pB9*Ee>=h%FslWQDt{%Dg<=Ee~AiUi2*kLI9=?ag(92&ZW z0_xVCA0V7((ygDxV_4`Lo#6MU??>|JAqeTYmw#|;FGZ~+f$}lm9bpAV{%i8jIRhLX zJ_LfGC~tQu@cxnAQZD(8JJ5ePwH?{gi@ITqNL`4-AgtqnD(it+RR6@EePk>Jtj}R! zuzmK4K%^d>6cx`X6%XiH?KvOwvPOhn`+av$k0%mJhJ^RJ$q2AUzkbEZ9L9?19^K4Y zfdleDM;0-ycd0n-iqUcRuJrQPR|zLc%#b5e+T@?%K`bF^0_0O9Pq$=2aAl9AuSxWm zWRyJqbrk{@eDI{WL9+$VR4(K2;-Epez$~mw{}oRJHX({X_^D%wWHfoldSZn9kKJL$ zbDr_yP!W5coiZc<=3l_$QtW4vK~~iBiJ6*U0T>CubP_w!FX_rB@L&c3D38%zNb^cA z#%K@o9-BSkU;yVb;yecB|8BWH2+PX9fx_sqZU{x2Erd+(Jrw_a{4lb!WEWb6v&hPIZHH5cHXQ5x%8P?Ty631GrS6si0O z-wq==CHM5@jBzl~M7(!c0b=BjuO9C=V@d}CyD{NVj5O8Z1SRz`S{sCe`p`&6^fUMp z0At((Aj@Uu$)v^JuTlF44?2A}TO%q>8W%p^-t2*oCBrEE&qGdm!z>mc4(PBIrfzH8GKl7W-7#TsKMqWl zv^Q-(AkGJ591YrFHLzhrr_BMZb;6hsiJT`Oc^GYAf)QNWxRGF{ptp=f{5zdv63VM`eK(WFfu%l07g-ldZm%>aB6r61^#l56>QD8_15&#SZJOL7&C4UBvP#nx<&^-n0hzvT< zgA@WliIa<~3lT)7-N)5IzfMsJQlhiGJ_?3Q7!$liY6#rW*Y>AUK%^kfFD3J0848wf zYxh9x0LlZASl}|Tz-jdDjPDp38R`EV1>vvp+2j=p#vxt00a>rHE4hB(8F;pfj=AafRwM6t`**=1_v8$m z{wA!B8bduJ7^+3)0JRzTj(9uuHm!^(i6UfGR9?ImPFZPHF?***8PxM>4q$7l00prg z|KHkL0dV{Hf$W1AB(P-_f&6SHCL{lUFTk~>%&$wXlNogpC9N@9ZE^f`gw}52DW-bv z=`;6u6?BAO@H}sQbHIk(=lPf6=w!^OO@X2vxC>bYU#{evb+*m`!eB&6t4vR)?;@y) zQ8A7dGky`K!p=i|QxG2Q%pzoIk)_}DN?(`zRKd6zOUuq6H&_?Nyp(Gnm?v_7Uh4a! zW*udz)K)P$v3;9dQ{ssC$n_nU>tnHY2{_kat`8kwO?f>6(-&kI0{nu|ks&7us6;^3 zSUuIy8$ZDH1A8WLOSqX*KxPOiL*_$N{%adA7DY_mN97NH9JEbV>K4ypq}lVph|15D zA4B&Wj^2Y_tab%@_D>5<8NT_=5YMc)tbHq=cgwoNS}f0^vT+pA<>t$ zA5&OZxa-e8Vn3javFTvH4ZR<*ZAd!=4(VcdQV>W0Ac(px;UMGW4;rxU`czmxVim6zFDL+*-R^cin!RK4k(Eo4sOPwc!a zz#@B+N5R0Pf(8zJ#x{(+Tf*xb+!Mk7-E;*bFgh9!W{oOB2-SEK(lykR{w4~O;LO2t z(KR^F=4kbAf>E_|B;fB&Mk-@7sf!=DKU^}5nSON96qV5RyG zr)JN?4SBG_Ew+Hzu-uq_dk7sJn8IBFTFpZi&gcF?98w|fWyWwr(X3*cYBHFK30^+X z&_Ib5BF*I>N#d`)mSC6T&5oJ_2nBqY91^s)RJ{}RORm~}yz9w>?jJQWGc#*k7Zs4c z&TAoxBTuIchR(P5#x^JxvZ z!+vvzU=Z>j=Zzcv@L5HCm&ndBq?cG8kQ>4~G9Yt900WJeG%5}qI}lI^WkhoZ138eq zq2NAw#Rte}DpZjIG2p^PF0ED#$B!uxZa{%o0DnJlj5E6kEBO0kU@cGPvn*Y~EW{iz zY>QCpGG;v9FR%q_F%koU*fd;FtV=#LkpiL%o(_VvdUx+a0nTq|NE4s@-vz)53JrV7 z__cL?dwT|RDrW_mU63uGP>~!F>6{p}nM|N3XhuXLNI`3+*CqZe&Zl6tnsn{SJBWS^_pigMp8&fSl zzzGh7rG+$!Na~??u@b2k(;9CQ+P{4y*Hh#T2@V;zhv2YioWkbWG~Mykzk64@Om6|^ z67{H(nIZ$8a77C1@Ky-iotMk_7y@&j7BG-$*x*P5nT8ITZL-ULLdn~*1!G4$S}-d} z2b#4EP|7f7nxXY2Bvjw0m0F@KS!pS!l(NbfW*FlcB%)(@-{A;nUVr!z8stPI%XEfd z1BnhWU6Tsf@_>wA1O0v0XNjr6|4)~#f zIGS&4Kf6wVCEFF7DbD?%rTSXWT{lS{md2_RWT*llEtq@|_IOzI#jq_h3y7E(wio+X z3!vY;0q+?gdJ=teby_(FTa`B!`6R?X2}r}ExmA|$1Q~=My3uZG%xeu}`E<7u89I^I zTL4hk;>3WFb08<_=-6&xCuT+-X!B9`z%PS75dw&rPNmqbW0QY7z@SX;5Ju$^0F4nd zK0W0J%0dQt)L1|@F2E8zs2$i92k#YG7suH3hB21v9&$2I66mOqc_&!nXzE{@efq`7t`ithWS$btJYdnj zrWObvhs<;@z${8)U$+u#Fw)(^k@M|SS?0WX6!K2YNT4VKJV5Zi1f;Z3-?SwFpK6_W z^OHaNJCr~9EKg_P;KCG1C|IIT;N*?G+ogveQi!7sgO+;{l4Ihl@4$as2g8_wsY@s3 zm^CjwkrkM}Xkyw2n>tvqu+nURQr`wwsVSK*cnRTre|`6)JmpR7dv%fc?v0sdgyVv! ztk^XcAVAL?4u`GFd-skOd5gursgW9G4wRn0K6_6<%%bJ1?7t7_VopvJ$<9J#@?jVz z9E383!bH7huj}?@e7ArWhj%8xxdj=eX{~msc87k#K?tU;(Zel}2~QAF-g&|Oe%^pQ z1dk0m;R8o&BVj9U2gj%9B9?PvX{BG3w|@RyQv)VkP)GrDFEKFz*r(aGL;Ybt%4#f zn8_S^1q5Xw4N3?$F$nvrEQz6zr&k*i&m=Yw5k)$-r3b`Hq!B}2nE=w271RQw3^e(+ z@+V-Kh29v9q3B>nbXx$PiogTpvmU10>;XGl5K+aD-s781x&<6HqwApUA@84vh4Cn2 z01NC;DA{|llHSAo6M|^ad23Pd`Wy!2pFkdY-FEB%4l<1)HXD4fFM7^Ap@KdG>Jwq( zl=?#Se_Ek6HdyMb+K#!9lMO;G^fol8+FSOf1OhhHSXB4@_-%(@gxzcQmxZbA&zLmdFjrxX@DU(D6mk86`=$?G>rgAu)OxdU#YMzrMB?BM zFb05I5!SU-_?k2c4Y|>oaAFZVM`WiR{Zc}s@-m>lbmHLh=xp{dwZubn1Mi+V1ovAEw7Zm06o|#QP_LU7 z@E$z4KAyTTe&h>P1esn1By$D0FL-YDTJ9|Q^Tt>z>z*J0}YL2?v43(MEd_&P_2(Qu7V{! zXQg#R}q2pLxz0T^~sXZxWAN&2Mfd`II8PXk;B@E8J@RT-#!S*YT2+H#JL zJYZgrLs<5VjEvU5o*K6sc|=60i6!<(5RVjWp(QLY0%l>Fo|Kvz2ZrLHTEqX72R>Vt z;S+I?k)pt5C;s=V6}5mh4jiE%bZkDGK?)%)BOXe`EY`I8GBa!NThs2QeZ-}5u(u)2 zC8X|Z09N!OE~p!=dt}gL-);f-HgXA&uqyPqp(8Y=kSy9L;rVH*s~Z9g-o7w2_#h}w z{_)+f`37U-ZjT?6l9FBu>1>?jHi$Q{x91EF2_a=rTreZVudO*E5?r&#dKY93WGWwsL9mnn(grDJmntXH*v}r?W_ZDZ z{*+krOg!-W^#_BNH66BZ%GL`x+`o!F0umii7NW8Mf{KA;JO}c3UkI^eB%M>!B_vVs zHq8a=U)bx6M8(GTz~LZ-XPBF0?gHNQ5)4(FRGp`Qybq(su}w`P$U84Ix+mHcl2}2n z#lXUP)onTW?70V8Q&kE-p#0wxiGz%-A3msA+uAOqnK%W_I7Q4j8Tf9)J1vf}lOo^W zh5dv=39n8==BA;|zskgP=S+)=oP6Cb`Kc88QF=-$|xf=%@T@GSs@{N?~+**S>ZAwiRh9B zDwX&Ay5Hyh<2{b&c%I|AZ!XvMyT0GgIM4HQexNzDe8$bf!7=CL+k5Qdqx!cJxEI^2 zXX$b)r?$eoWH58}kw_7dEn7wr>zV;W-MM!!q_*~;fWfzWkjrCkprNTr)|`XmPVZFr z7L&h=3(Z6kT6>Xg)$DXbH);h^fH2U&cqt?$L$_{=dpZhmoW6$U=EXm21Ju+$jM3^a z8BxAOf7KaZ(r!0AbNc;x_pJ%9FQ!vvq$MP-NwtPVM4Sl`6HXcC&zEObZQ)Iy{w|-`{g% zihK?W)kkwb-?m#1OBnb#s2w~=V`F1e?jI8JqKrO9R!k?x4)hm23?`aUj~G~5vej11 zM1Jf0@WIv(EZM5ntH%r;>S7KAzYNv8#)S*JyjF2={PL_lO5p0`V1et;o*f_t2#?az zVB_~h$;79|(!#Ir-MdYvi|%V)rLfM_F3sJR+qH|lbqr~T+`4uA%NYAkyoggRFj3-Q z<>n?s0^y_d^fF$aZn>aDaUH$3Mbm-v8NJhH<{59|uE&VqMBS%sVZqL}FD3N#Ymnf!eZYL&&16?uloSn^5DCOmCiRf#OPI45X zNcrHwgIUgx#lD3M{I16!m6{rZz?DC4pbAJm_M)u~y`ke*3`AC45;xGw%<1_srORLQ zRIKdvxw5S0CuIy$)IM)5EhZX435gq5q(ee7j`c&kVt~mEm!hyawYbSr0jc16TA-Jp z9>`rLXrK(p)spw6o}$mhZ%zq!Bqt{;6&6K1Pt80WHJy7}RcGk3G>$%Z@A;1`*RNk^ zU%Pg^z-pDSuy7ueiFQ~P;=SU}5&ypGi#yv6EcunN#mE+zRvo&arlwXV{dU0g;`PMI z#*UyTZ!OIflH=;r+ZJ)*VFn>2FV80$bz?CkB*gE4wZwCKJTxLL{z?$x!q#oaeEH-z z$SbN>nYij%Dy36v9I_5y=QMhG@!7tOY`6r%`}Yf=3|9zT{(ahKy6y)6Y$2J$8SZdE zaB!Bg3e#@qIg0m9@tZGRzixygo3th=sHMT60IaPgO&9RNdg~n-v;{>qrY*il>mtQ( zlA&UKL&HDE3}hCaWB=cHc075QE8Y+V!n0djtOaKrFTVthsb<`H@E z8G34JEG84(ye);M@z;*nE5|{|DuKRX5p9T2Ce$D_u0Qd!0Dsg zul!m)Q$EH9<*sh9B$Yjpvy`Z}p{`cuRwt9vlf8LXaf_PiT_@i?y#yzTd$M-gp|^3I zk)uBC%BC3Gzw@=Xm%z9?J3l`xeEPCPr?@yc=#65L0WsOb0$qEGxO|V9gS73%!fR5f zej&1#asP0r8uT6`+1(C5q;IUR!U$j2e0AC|6P6(|O7ULFEtsht1-Q6*7)vcabaRsp!1zTrd z+T+JIaP9J@mh{fzt4#rLgRq?<{=_O$B-1K!D$O1=8A6IH(8^PjZ<3{b(osAKZ{AFp z?Dr_dsPl>?JXi70FnvM_M8YEwfhQDJ!KwoTnutHQMVLA|Ivh;8W;Q6ZcWceCa4hZN z2JE#u>SaJk$f~@IGa71YHR-IMK2v8tye;5CSU~~8vT-jv1of)8u|aJ7Qy1pqpIZ|z zC0ExSq%ZOqRri{C0YHYEn|lNkB3{f}DYYqGpKd>Rz-+)x1bPjPjR|w?$kk57r&6KG!Qe@!Gs+`Gq1lTbDF)#;%t!v72*m!72J37`BMp77(Jz+$g`U@K~%d{(Owo6Db z;U>5`)*j^>0=9}&&9Z6#;3v;s{wGou?r4NacRrz*S-q3B=G!3g^OBY(11|7;y#+L; z9;ph1q&FC4lJha8l@awtL0F5EGha^;VCy@bSxY##ZGZO1eb95j z`OuLghN=^p$I@rPAB{lu87S3yPr>V4ReiX2%jM~uG9qG0Pu6!%khbgV7)( zw+VL_j^ccSTnY#b6q1t4ZoH_jrq-Q^4r3GyjE94nnt{{qX5xAY#nb0+j+MKuX`fkY zP4fAeId~p8AXyv2dEyTlXD;MFue8 zsaUF%^L1gN`Ju9ZJl*yAk|UyAq%d)!6f_ zS@wQUWyr=2@z`=<9#x*fneSR?HUxY>%P_25>E$i+z%$naG5Kuxa0*ZZOTu&7g^mHOc)HF0~2&i0fo`?y1IvJclIkm@?rd)Pm7Ke!^ z##>apz03LXk#*WcH>0t|7~;7NjUB(ur&Tb2`zMr4s zk~}+yLEcD|k_mDZV&hND1`n|4rii}>@)wkpWQA03_Y;KuB{by@C80;tyOd!ig0hOx zr+^nHN?&>oB@7O!BiF85hZMKN8g&M}kylqx@WY3K=5CE^3Cckf%q8->ykO8skXB*@ zhW;YqiGOS1Es$Suh{{-U$t%9*`OnxlA`|cs>7d~vSSK+laUhG~Cb6DI?Na{!rsSYc zqq$4fH{BC*8Rg}?sDxpX{U>aTurOnGB8wS$W!O|+LmFXE^83c*Bo`;&K84fUu=J;W zNl;dnlR!Ue1{atKl-Bg+@`8dnq<23t3}<#X?yU4G@ftlU)0#ESgB9lqm+f zrD-|;4)K9tN$?vfI5#AZVtO;KEOPr^k11~a&{uqfQ*cGN^pqRD2JAYKFV>PeYTdf# zcP7HZHnme)ks?C-_N^z+#q-N>82T%wo}M0J8A?c7z&xAf1fUmgC_ZZWQ1e+kTU+betF%+IkEB{10#8uqtHv4}a%h}f z?njVuVwiV6dh{s2T%3a=fr~2Zj#kCltNrngi=Zx;_csOp&)qyQO+Fbf3EbOC=AVozYeniGKh8jd zM>i9td>=2E(7Bv(F*jY|V+NSkDA)W3Oy+t{{YtCPW$FE`PiguE%i-q$98K@;o1Nz$ ziq(~}JKLQwS6sQ1WgRznt&`klQPH3)5~$f7lQq+i4QSZi4ENNnQ1IyUPqQe}^M&H8 zuf}%|5FO$w0Axc%L_|zty@)) z)A9sGgj3GTd44vjUi!C)p<#z$|H9IoNc+K@`l~&_+sPm>=I0$oHSMmO0OJzTs)POu zMjHk0^PSisH3z5Bz6gRW_A z9&ajeH6sgeiymTjU@Onxy?fT1lZF=?R<=JtanMxkczVGnauM1!5=I{VR~wR8fWtsT z135s>mX6-6h=m}oS8-)}zC!9?nlQ##smA87z`=3xx)>!^ccqjv`wr_`n(71Jw#@xc z3jo~Br1F9U0YqTue8@s0C~T9SsF!opHy709%P65!-q!8PND(J%J3n2RT5C}#g}`i0 zSXb1g#2>xx>@HdJJDMLpepCS{A}Z?REKZ*5k#y#;t&Tq>!s4Y`p+DL9?^pX3VFiUa z-#df4-0G0n2@49+p#Eer#xwJ!)ltaP=2rJeZUEGYd|3E(bJ?4>dX7kapri@7$HCFz z__i9aOeZg+w(5@+$F{WJnP_~zHgAFmY=4RN-VOV|td^v{KJl&T*I~-}7XU!uyaZAb z=O+FNYM6;`o~=|J$E;58Lg#(yI{$E3h=Tdq(%&zNK9+te67M99y|*s;IGS{vGsZWK z$jY&PeWsDyGSHIiDKDNXM3FNXFq%JGdsMD%)b&E&F_j@bCgB#Q)c9S94>j)w>b3E9 zVfxu4ASH~W%ciiT(D}GHxaH<9%9@3lbL&_qzNxloY~Aqc>i8usV%vb>RyIh4gj!|c zTBVyH-y{6eHw%^Np~wJ)^q-51-yd>vbnH_=n@GH!fByQV9Ak%)l8k;6EnK^FbeH7| zoBVIO=+5za1f4-&5FdWOGju$xn5%y^GqYf;EDVjC-`JFi}bmmpg`V4`MA}gGz)>3WV0iE=zP39I%=AX-OI&D z4}^NNsOX-Y_u*$-2_awinz_XCqeHbrGrq8Uz$BN3Mo^d5v7<**M7>--y6(N5l(Yur zM11^b?pofD2;NRBR<1k+u4fe+o6!Jgn)LC(!A^P{l|Kc;$S??%5n?2ybu&=u7FRjjz zj-|NH5+&ndKPh~j2Ed^M0vOUIS6yVoL98eZq@;YJEGl|*4DXAz=iPPoHf6GcF|&^-7_lRL#qO_?#@^C)MZNp20xkU@B&-vmRf-Fiq` z;a+&$H{XO4jh9c{HVgGhN)a^Ar0YN)0yU3Zh>A;i#}0N>0RSim$GS&<9Y1=MPAl?u zWMo}}ikjLbNdOYKLqse>XaooVZ$Vr!0w-9!Is;8WK!D)p&GZ2KdF3yNqugo^Dhj?H zI_*3R4+)^(Fr8;a#S663<@_Qks=NM`c*J>LXy9@_9P5)0F83>$_=oUc9lag^DH+6N zC^D}c1kgt#xM$D01OI+*E?L1yA$DZZk`|nZb%qcx{nbB*sOQpZvd(VLV+l&s9f)Qe*>%_!YAPr ztjh>#6r7&u!weP8@;(t%WrLZcXY>NksKiLM0x`Eduu1pCrYU^-+3(*~j~?Y*Bh10! zabHABcfIO{-0Iv&kdp&n!M%+@`z*E%O;d@js-yARmO>VE^tW!Ykde1TdV0-($Nd#% z%2ji^Uu(UZStFTiKiJ-`r&Lf~(yT8bi&o?&x>SKf_CVxHTn8KW+*jXW?DToP8tsXW`kBI)d`6)8~T^qz6di+($+E<8qKr zO5$W;S=}Ynr>>af&RQI_l^ul?paIowO{MFZbLY0vJhXatS$BdFcenp2r4ooW!4aQ7 ze?H|W>1;N6!=1IVXy_2wnHreZa}x0fh*3V?pg(!KEPrqgOhD_7NQu%dUD8nci zHltY{9Ua{$jEGWJB9oiChZ4z}cy>d(#M76Bw{G8Nhg6i2%zu!2j*JGNry-IBqDDX$ zMRO38A+hwtMNh>-9{gMK$1qMjb>W1FujJq+BcFl&Yt=bGZ(;4g_%a<)I zhe6Q$u`vYXCVoV%nz#E_NcA#jKrwkpN2eZ7WhKybpbs8LtH2hZ7Xe)yuiaDBeA672 zRx+kBWMpJkF*CofR>YMvL*WpV$qIo1xK2F_)wa2>@yArDFI!vvfk~k$#Q>@A-Wwqq zE-SLgvd$8wadwZ}{ z>!HJkG2T3(j9hLxz$j}yCp0|hZ&x;ej*CPIWN1d@`KcbN5PeL9*ahl0wVa-%*%+#m zlN0HSfxU5Za#m8gaPa18JPQ%U5w5suNUo;OOTj)2=;gOvj)agw6<-ah zb%0iT2;D|1$iw3V;k5DEwQCLv_W*eSb@;+;q8*`n?p*FMyW?qos=iz=?k?g}%F4>_ z<^GU-)ezv#0g%7Yq|qP;V_vSuAugbxG_d^zPGnBw$~ z$Z+}KN@;7_S{Rx4T)#~=ivS)q-Sr^1xMu*)SqLQF-Ab^L5>iy$2oG$A*(;&uF#aeb zZ`>ztPk`hW*VwyvVXPpGnDQJ5+n|d5G$BGUfPc`UNiW6@-*WoCmYLid?>5x{(s>fU z7@`>ml6b7Ti z_zzrGB!CaGOG1ZJ=2Io2m7slh=#Hs{RhC#*OPPmy(!88jWcT0RXlEzsgBJ%3FSh(+ zSLZm>7Jp4d3uO19xB(=;?w(Tps>z~n<2Tm)u&BJLA7ghI7Sqs~2lqRjbF(rFv-R4Q zcwKO>^jMeCL$~JAa!;Ak`RbO^;&p$1Qz*VIz9j?H7dWcVtw-2}KqCd9iqqy@#l-Zy zE^u{ER#&xz$XTw;X}OB?3MGy*rDf8TL?vo;EM=l&rv53XEX`a$vLqRwLkbp(MM*Vdp_6@+KSTnThDs7yep^qk)NW=0IG420AlK71gY z8HAPqSP{(;@Ggm3iyNOuMBI#`1DwgTXnfi|^hRcZTNStpD<$_rg9oxyaq)8z0d@3` z9fPAeoFR!k7BwjuZ9`q9A6u3VWM<*rl>;!sC96dQ<>f1TnkAG>K$AK4<~Cv!P9Ia| zp!Qm2aSa&RHiiVQUS|XJ{#?=|HRoCm@Fi3nBsF#~v<`y{oH|dNRWd5s`%AaQHFmms zGj%q$YJ)*C-*zq_E{mUU zr01NEhWG`^1sN$_(qP|V^n648+b*i}XTRms=?&_dTcXU4we1=G9gk&cm`U6-d}+Z# z%A%4RD)#d}p8;RqMh*Y&?~K<^OzzwO^^zHApg`F+u6@9O8xD#-#GeOauRoKPer*C; zmOCUM2G>&{kcz)2Qy~zS2N8nViEM%HA3rV-g6G*^s>?F`gLey)t8)X!!=;MhWY~!2 zD!yjViPvbwJW5`FYB|*w)Ml{(lN``&K#&2O({8~|(>D+%674HsF`d-BA&C@mvMC3f z*@)92k?5B2{m=jyIi2@R0JP&ZST|bced%_Mci)#T)k;7qJbLm)lht6vt&8HR}yT!<;r~>b=He`?y8jh(85f7&8+J@Y|l}pY_Osh(K@2!?*;p8y*@|Mf_ zCf;Fg_LV(wbipDSUa&wp+wLDj6a*efM@Bx;dkqJWa(oAjnHRFAW+D81N(QS$zzV}8r0N6@N`V)xo z;q{C;CJ&!Y{047jJ9Arx?kDTPiWNGFSC56mdFdm|wqLhw@Hxz0Im z9+F$B!?fDi>eB9`$BxP1aTD}qd8h~8hJ@I`d&bToJ>3VwbX8jP2Jj;`!;D)n5GFV*_u z(-Hi9W@!6mh+~rQOZF629u%lJb;m^IBKd&rrbz)X{YoJ z6@_-~UwgJ~V~bGywu($@5)BPluX+%Wq`xk9x}^~N#Nbfm2ksctEU~oe%huLV!E2W7 zp5&Jf$r*}B=odiVCaCLb9x43=M@PfyIlyp#=gkKFzP}b=DxG+OkAW{2`Vmz~LGqUC z%?-$V!C&m&EplTHixz+KEghfT0#&o(~6tcjF2O$xb338iKolw3Q zPWX6;?ApZ*l0PT6Dq~0!IUYwpRj13J98P0aT zzv?Z<9}W7fO)jkc2B7l>gqW#-kT+fW!9j>wqC3W`mZQ~Mh9aZNhW44KeJ%zFmkb0K zossUq#BW0*qaX-$AieUbHU^t51Q!3_zkje&*3VjbOyR_3xCTHLq6r~HaN0Wa`jd{5 zfsT4Sc2CXqiu#>MdWNJ5P`&C}4vvcG(j#vCLiW?G7Lpm|&ip+1q(qH?G)T0VGJZc@ zvn{H~UZ>#y7K4z0PrebdCX7tcBN|m*dTvMKaO%>ZpKj?{ycVGeiC}#qMNWbTKDWj< zhl!1?2|{&j0?0KgvtZf}!jeowF|zQE9G7hEC9OVw9yo%o@!k2fz}y=lFZNXA9QFqg zikg))@ZdFx$03@dEv!L-fhUoR2~0>_L-8CU9`%6m1`7Ee4r>uUzK*T)0RI33+Kkj* z8|w76`LAQcrGM6QQpIins8%5IC?Q(U_y&J0Ed2J7MKizYBTl-;>==&j?UOROjkaEeaCaR3#x)7+C{bP8llAC0;u z3R#8h>oeh0n7|PMAV@3m9-nK|1?&y!NP1 z!QG&v>zljTUfX9qS-=qM-(icqyg!aG=nOvahZfZaY!>dxb7%TN%0X725W0@Q!)2?k zflMb#z(nl10{@wyjfEEx77MOH1A62*vBLxA9aB zbdMebHc9?68Pg`;we1NFWCHSZa9Qh}DtaH_bv$lgW?LJN1%AY!w}{ zGSn%6Jm1XZZzXNW{x^;T$L;uSsnnS#gTrU&jI0hMqa-KD%G}mNz*duz@)E6;_%BOt z-Aa}j0UO$c2xs00#z*lZiOO(V2OmiDr!w*d2n8AYP7#j}FNY*6cUb?*z@Yl_WGtqa zuUcg_NWQEoa@^>16fHZNLja9*0If=z!=O^(a&E)JeKK!yBP9eHvy+%V5inBIq(t0F$Bque%5i#TU37W z0%J3GxEZB3;qVba2X?D_hyFolc?tnZsq@$epZAHV(U89hdgKrE$N=iadVGiPrhY$l zj<%lK>bc>E5L816`bvPx5$#a50vuQayu*0<77kBERlTykFLGgUR?ZGCvzS~G>GrO_k#phs+uq((Pzs>P&Hd4x`*G68tP%5Y z8H&NuT3~gXWn?(e*T#>a++ZpmU z9i(5?DWdDhk>;7xri_zRP6v7k`n|c!mg*?D$ZuqzhXHO+)&iKCK4_kbK%XXKeMeKg zMn$P+{-F_oam<;ak_p907XU8w&08VEFT{9%x8uYLUk8Rjd@YDta}1q{CKRshDp**e zp2j2|+GSP%GP(nxM?QSY?=gMzO)c(4n!Dd!@taWYX@F8&#t_KZnyFGp<@biR_I6G1 zSqR9Ggl8}wWaHxGxRZ{gbPD0u!Cnv*Kuk;!9;C(BQx>C? zoxynG1O06`g^65QixW3;vN(Vi&pVyqNTm4weFZ!iVd*i^T-Hoc$7!&Wi*tLuyVn|U zl7lQ3Bh_MJ*(bZ~#ksQ00G?@SX^{ZDn2kTzK;oxIZa%6%nEHY8MstGkfXANN-YFxy z@piysL8)2~rFKTRG)f+5gH}NDF*QGT?=Fp9QNKP6HGpsJww5u$dlDNEpyA!?DoBRA zA=fQQI6i|N69N~V@tzJ8P?&cp$w(Q)&K`^nxDH_`7+0dIuMF(`D4}FH^mR1RnB({C z%MU&tYI+MihR2T+*&U!06ul#Zq!l$9>SrrlnH7T7sC>K@N^)QI5d_@Q(8M}{ z?EPv3r@FanG>MT1kOMRwYXB#CI4<@MQ3dFYHnExvC45jfAj9Jj$U|fF+#h_MQYV|I z@nZ8;^&zlwWP>PBGrY6eC*!{OAKv^7<0pA}&A^sBo?W>DQvCcs6elv)90`(kCb9|& z8&LW{y|)Qhke5(Bw4r_O--Q;ybN>38KwLYmjYP!VDE7dL}432FhHTvHy0m_@=j4d>@)0 zf`_8r0F~dj{CuV`2T+eaa{#@2gzC^u0Ld>{UCW#AcNcl(@ZnqX$$IPzh$9(VPpj|_ zN!@k2NC*fEqyM!~MJ^75GTHlw(bu_d<@z8$%J2GRk~vTKhcbmGwRgJbVA`L7gc?nJ zeb`ckOb*^{=-+R=ns!wNHOHM=;qEoiaCO|jc$~Lb?2cTC*E-xoEO%0l6d|hPmBBzH zvz?DW&N#b${@d}m7eXqV)||*$ku?*4#Skqu#L0-BLP|;&H4!x?I2jhs$)$hKJ8x3z z5x{e$eU(a>G+ZxKaFDuP9R*JJTTS7`RdoTSi|rugK?`8xT%l&{ln}0luG7~S)Z(#z zw0jIR1euUe2sdUY`IF~Y%naXyuBUC7H(}@I4gh%K;eA_(9Fy<#R~)HH=>x*c_w()+ znqn~%uK=tHg9^qK@;)_nb)jwBatHhckW9^zx+5X%jj-t_7PGED+^KU5$yD23MSKCHRz; znfSkflJ1B90TTP{)}5_Skwu}!JPVpS_$d$-3TgKq_&4U0OEQiu zZBl5<1T+i5x5uH8_(Nkh^pJDe;EMIxnFbn3!n0O^dUzsEE)?AS@F_@Xp7mDV-GvmOXO%_hVlP z)fKVq&sp1jf#1>asd3Xu1l<*h$uS1 zCTx59)5aXzPl289MgP?1pc^TBTJoka`AMSig-iZ$p|$g)qHsqRdIT9Bcfy11jF zqd?V&CVJ#}<-xdoQa(e)T-o?{ShxQl6uSr}!ho=#ZvG>U?87K!RR}0^5aVbeVPTCJ zbCL0wI*vwtc`OXzdjZDO$-h0b+xD^I44RTA2xXz>@<*}o9z=o$*bFibO|UW|R$1-> z=59;W7=WJE#Fzl|!c#m$k%dyHhz@I^Ci(PTM;A*JiAaQQn9W4&z8k0xDw;N{>}XKe z-KJ_l!DI_FHir^{52kFMsG2}#7o@*nnR9{g{k!N^Y$<0SIYCv|N-9ox0HXE4e-jyKM~XUBna>ODT486J5F|K_h>4N_ z7YJ7jKZCnA_k__o;%)zrC-$2L-Y(iZyC*5=nKUBN;!nbxkG(8-ACaC=SFx)X^sF8! zey*)>(z=SS@<91c1n-QC?}V z5z|#->0pn%`{>w$bRCvx4WQJoyh?m2|mU8_3Nk3J1!ZL8Vi>m zWUp>n-0Rb1ypZ?r%><1pqi-&fr6if`0~`V@!74O-rL!3iFE3sdwOB@J!KMAy6>{^~Skq(YWZqAm zFQQy)`ChV4|89PA0ROno{g6%k{M7E$M+x_^m{!JX#ys2VnE)sva%2Vu_oo@3dz*a% z3#p8P;cL-rMc#(Ov@Cmq;eWfKh49@WQdpOhQKd=0^*gL-h-?1ChYxj?hWV+nZLyvY zIc9=EaJ#lm#O?bqG$T z8@ZkHjG)~|D(1_v@a!GIvINa!|p?)F-4?t|tDsg2OW0eE5-!sE!Q&x8Q)9zC0b?C8$lR<4YR&YnKRLzhahVt zW5tBpM(uUX2W`&$UWo(i)~yQ$o=hAPfP7w%7o@naFe`OZhH471b~=x| zaf`Dl*Jlbr%}6{rF`D+r@|X?=j|I_@AyiFq?#o?-L@Xn3Hd8<CTYq>!BrEgOmm%)zYHJ+dQuECt(NuwNS zSp4epE@m-VrWk!b%X|UI2&b0zFfkXE$^Iz zz_bY`4dO6}sli4P$4bQ?Z&bZe3N7zZ14dS7;WSuxC7b%CxFboNdWZu%m#3WP@x!C{@6xo!g{wfV-w_M*{S|BV9 zQ7!;-4Tr1*(|#$rxojW=pLgjUpdLSaY`sl2L{1>O{mB4hU$H$G(ZQg90XRxDdAdLUUQbzge#`&mJ#dCi&+dKDjGv>Lt!Ad0jWq@}_r z&O;Cu7`(vcBM|bTLm_;0#j2*ye!X^?C38e*bIE|>gpP!=76$TAaS_=M(L$3rfcK7+ zL5fAxO&l9He%^J>BE2iK{yK3>!@itQ+`mI5|}+RQk3(s zbeI28*@Kvtrg_AwOnLd4@@CwW^=Zlbhm_`uAW1|rYGp_+*YoHBr2@-`^B#j05!bAG zy|)592d9n4-*6K0o^O@dZ3nL(>RaVVzzIIV_+FF=;Bq9 zmttLnKO|i!RgNd5H{Hey58EKtk?IERA;46k+lP7^LU{QwCmCylV4*`GsDMNwhD-$dLe_PfJ;>|aTvj%jWTf_;#9$X`%dupgSLmL@LU=FHQG>9_ynMG ztAiLWCD#sv)dkpn_9{_n!1;VNRMP;c3;-?vH{%ggLj>F&L~KKJpV(Q>dQ3=2Xv$%` zMS>hu`pqgaN5^#bx!vW;%E}6G%xc>N;|G8bAfuN#&hHQ*z1Tp?Uo|zsUAu~V+VEajOuE-XRd%lI6av>FZEaG!`_kin=O$MVA6TLRh8hFl6Zz-8vc#EEqSg<1~lU7+H)c`;Vo*w?z1vVtXuLu;$ zpzxRPO1{k+9{jc7RhKZ&jvV3eQ-+siYKb=fdm#-vCRkbF_SK@+ynXlXhtf1JR49X5 zzh<|$5>XEvC%flXK{<6xBehw66`tmg08%@)L`wqaUJd+U2QQ{Uy2p2jXe~Q%084@i zc8|R`aP1cxVBX3CYIHk}7*y$i{)x|wH9D!O#uMhDx8>L7BPJ0d`}0thF{>d;gTD*r zNo6p}gB}ARk{9S~Qv!No%RB@nE5EOSrbW;;aKu`8-uHP{c#DmmrvzK1g+xUei4_lY z2M=DsVm%P@6`1bEdu)&C4D^}mmLOdK8IizrofnA0i#5PFZN^)o0iV*$PMcKA7;<07m=#5ID5`fI1Era3YH@CFP;?}~xn4DQAF42rog`^|W z(s*4)sc$bIcHbW`t!h7AG>WoW6LlCF>7MhnD>#gEJKXuNB4!g_WSfDmcv`WFnK6xn zF}LxPfsY?PEbEcTb%cFvvies8yq-&X*QfZ8td9OM^Q{RcD_b?Wd_f_QGY`ID?t}d0 zKjl7G97hi!S}*p(Y#oom{QT?Pup6+=z+bZ@S}(Hx4+G~tPi7FXG)|v>Y#RPzh`QlF zCIL|m3-4ACCP7!*nCFtw_2Z7kVkaQ+C5z@TltWC7svU@SG7Wc~xQsh5fkYej`+kJ* z_t|4taa7O(1Ox@W$Cs>WZ=cB*PZ(-422~%Wb${^b$$CaIyo`Me^76i$gwXC?Pq) zcpvj$M-EuZ@SuQaZ*%~NsU0uNdVSNf70P|)91zH@vCdO>wL;NKAfbMlM{Mz}O^fG9 ziVcu4zyG^P3!^ykUF`+67;8i1D~Wn*w%Qzep+@3l0nNk+G@U`uo;~}{RDWXzG4_L`6$Y;_`NK4c(4ou4Tk_CPM8W3E8ArQD7ejf&T!_zH>zPTQGz-)iV87GAS4?v0kW?pt_26LBp?T_9yWjMXfoY$8i>R!nG^v~0j5nAc{bPT zSvb&{+MywSn2fDQ>3DWLj58NI)W^sU0^BlynM7{+U!hH;Bv2V(Sbb)*?^`sC0HaqV zPQb0t?bkOunp(IUeT6cM=K7<7G_#{3v2m&xNMroPx_n;ji{2+PLsb&&%< z>@zJ*jJ*;Q6Q51~yzv81ASV5l9P2 z!M+wXUrfL<0!(xD@K`%={A9wARn}|A*rHX?(qUCDv1$o@pC!c6iE~E$%wgyfc23Mm zlKUV4eJGlt?4Vh>a;40n;TZfVaizm#`g^zOsVU$Aq#4Ie>*ms}^F*aEJ3D*O({q0N zWk%WiEDj}nY;0_<7cO8o8|R8>2q(`-!H*knNE2dfYI0N**G7K+O}_fRRqL?2V+=bV zfQo_{WM6Sq^#Y!sd&Ndn)aOGDKVwlbtYfyYHxfYz1+fDrfeg{z0Y?ZLm=3*XOGOioPyv_pFguel8<|CA{BFblYIMkd_Y%s z_v=+vRndHBkhAE8*V@SI}7t~<1)h%85* z7^Zft!fvxQYu32QDW9j)N$q%u)*rg0FlM3Wjz&gCWc_!O^*#j!;zEoqj{WEDODJNT z0r2L=OEP7PHbxKw6%^Z84xO4PdjD_`I1;Q{U_p<6M8NDpm=cIL8doZ1x1{oxzZR>` zwzsyD9)Gfg8TwhkV{J(X-hKW{X83T&v3<`rZTZd8O_!;%%L9_aKYHJ}BMevpO%pFF zHB+f!pok+gGvT<{yMES?+7W8sjv8Mt8Kvt*MLPi^KHxUNGhDN2e*^L+>IGUcF|op< zB)obvSRACA=HLhjNBAdWUSR9L((Y_(7}Tr#Qc+QXJB*Q#Vv6q1JB&|PCxMQf{rVM5 zcg)Wk>KYnS*6-iH|Gm8CEt$Rh_6cFW8lW_KX_wi`_l#_Z#VzjH z@G+`&JY_;FB-!Q`JbZ25H{2k|$2RlwtKcB9f;j81%V7C6Nb!h@4Ufub=f)(+J#66* z0Q3V(eM-wCcCehl3=1fp+JS+pe!ml)iv=c9I2uOv_e%%*E^B@wWlhL@&;{WKZ~nE6 zyGx!4%8$HBz1(?MuTEQ9B;To_L1%kS0IY~lj$*nC4vv`fK|%veMYdSLiG`3HKx2$n zV9g|48}e-$nc+=e4FRLg@Yd1MbVMi6UZnUX02A=;y?gi8V$cm5aYC-e$H(VG6E7aE zjh*#8M&LmRv<9khT}XJi2L5m3UP)ks5J8QRF#_ZXW-F_AqDn`uyMcYWclfbrnV~p0 z1phQ3vjb4AXQjD+&hgU{Bthi}_Dmo|tiS_@HWx2bSy`!}9`%rZaA;_QMOvDn^IFH| zobJZZbzfmG&H*0?kXXt%l~5(2=fGCX>qKb_ptabi@0JWxoXYSq5-UN= zyoY5>EjUQQ*w(msb#azcmod@NSLNlRG-@ZQP2LN7UntnJMtpsD(k4so4KFggvnA>! zwiadwg-}cAqPEr8<}C<4I?hpDJ@38Ay#o6_hXgog>WAOCBwG}eEtqZH`gu;kX@wES zuXyQ);?H94p&k(tFp?TwmQ29O8k7vE^+rMA14D6jb0fBDBR=|F)&{6^rChAs@=Gb=)j!A8B@++0^Lue=4>{?-3Y+~NveOHEC+ zFfS_pd~L~_151WRM@DF@tgNb}R#uo-^?teZKP>=Tj!`ILBrm z2xzoXPS;Hk&}pA5E02Z72B)XIs;a)RaRqarN;Vh8$A8<}Qv7*9kw^l#A~u>0i7g-+ z7|CxVEmhalY-)Awih8;dQ4SM8U<&J+n^z*6c}71n^&ZAA;@6SY5F%bbU55}(jKM1} zeERmnK9Z7C%emDje5}DTK|h&!O}Wsi*r$`7mWNI?$IfitCBLzGx7m(PIv$_4XBS?y z*~;#g%>CfL4>|=hs|>GiF4%AvW(BF4u3wt*dR7midOo}ZFvi%t%{@?ys=h6D47Sk3 z-ig%whyYlxSp2J{mj^kT(AcQwnx}fQzSgbY4%^*6!_9Pm#3sxk?S=;~%3b18Q^j#p zB|$W{wE#UCarHETaGET%=-7E?jh)b-Kb(zlkI7U8YLc6dPbsQr^i97kOhq{vB=31?;mT65^FZmz?hJ#J*ClgU1ZT`RjST#5}waypFcS|uHxdv!l3fXs+y-~P;#IyaP{;oRk*Y2 zpwl0$iNk(j)11<*@L=aNi1&5KJGnMeft7mRYfw~S?Gx$`YTfs)dq=RZ0b%}zN1jrD z)DM2+t-nVNG3BRAjg~bZ=tXE#$Xn?8~9jZ?Vjxna!CMwQfQ?YIbY9 z`0f{#lASqSRMTX&RXoo9(Ed+>+pbVsOvJ@}c|=t%7B`)v*ti&ntSmP6EexwK4*r&Z z>0R0vHJvH{(htAxSmrs8o!M?<%(qk2kJ;09yPI2S@Agac^V|bx;?~xd-s1Cp=6r0g z;#5tk*)_XGYp39^e7>&N7cw)GmCfp-Dvw{$dh;{YNKo3M;RIj)*>)8jidmUv#El;U z0kZkMi}@v2d}}TX)A%7sn#nO&>A#*S=RZfb~XbOvr<`m(NH#nI035sJBjIRn-aT zOZ9G{(4p8x>!=*N&EoUC011o;wW%noYy@`uvR)|jGB-D$r8wqYlVKz^$mCl~`8(IC zB!Uo=AegyY^dCJXXS?RBA8YiD!=kDXC=*9%u+sg851;-J@J?(dlqm``xA_DmwqETO zQN^1(kE_fy&;2VlyVheAVBveab2K09|M3tllsl_UJ9CVr%IeE{BZ(6%9fq|zZ-1kNu ziz2@K&Y$Z!PdKMGZwuR8RL0mb|9~oUEq`dwvzWXx|PP(Fd7fET`x}8 zJ$_2K)zVozBk2*BXyRRqJGD~>QB)+KX-@=NUTPM>)JRh-jTyusz-w@9&a>Pvw>9+L zM}g0VCQONuc%nal|Gs&@o+>Sxn_T$t_k3CX^XIe#f`@yX;kP$pCyl8-(kgX)Si!oy ze%jvaCHIX|Ei{%`N3NY#*m}ZBZGgjA<#Gyo;~h;S(7Ji4mZ}XDNIXA>JNVq1Jxcz7 zR)-E<7QN8ng0y>^-uIprq`EQh>5%)}?D;}tsv@QGruW@>-QN<-Nn#>}8FHTwM&EPZ ze>;f#fbyUB!+g5~8a7iRcF|WVNlP6rR(jONFfdQmD)IV&?TLT`%!6+?tUb9P93eo9 zO0^?K&`q6ege65F*+;xwWd6T@fSt1cRh?v|GWpN{_YYF)4~8mvxjxDYJ@(P|eQ)7* zv#P?^E>3aI_dnl0o9Wyy+{f9k-(29|Z1^^!N=amm(aNBSi&=B`omg*0{^5$yJ{5hp zyVN@|o$=0D4^{EMcbG4Pk7fLuvi5kD=W*Cf!7=B^o!XS9c85Ebna$Vu6yIz8*c%vV zl4g^{rt$Z;fWeImq$VIwdEn8g6qNt_v)tLT-8L*jVu|Aeyd%we1UUg6JYe3J?>t~C+~l;twSV>M)n74Tx*h=k^`n^$7>>iglLLAQ zwMk*)f^tl6TZ19UwsY^E5+n>JY=YwA%Ulz9uSu14ETIz@iDw|5OC1U(xIQzW&cp}? z88jmFij9r1(x9jLHFTqF6YWW5g_{FaPTGl~0m1ggYY`C<(Qi0?e8KGyUsNznBa@RM zpjjcBLgrE4=)x>aa@SEhvKGh>g&;}9Ys2|_mFe)B0AonG zh3-};>Lt;vA&mm3*we?Zy|8X7A5%^!YKSatl2RhMapT79Alh;FeZ-`SpbRj|#eNX8RFN`N!Uy1L*xj8$q;h!-(YKNNb*2Kjow~b?a!#~aW$idL6^lXRl*S$QU zff&7;$6DC}eDtn*=X8+WfVh?xpDndQC3FG~+dX<9U{|kNWlhtELgWR!xyCM94kk@p zeBAY^4^#84V6E@R#;ya%OmL}eIIw`qIQbbL|JK;vrtj6|Uw}8?hk1jfl$5Pv`_7{N zSKf(dT&lk^P7hV>$;!ycmtS-8iU&vgcPS;Mge{_?vPfDy0s?IaqfPiUeOGDe``-PV z-(zHKEJj(e45jTzoai5_^6Dq=2nNn7elyhP}sFv#3mq2gS`o8|y$JnxXxhh6~sC`T5iM23`#e#5FcI7h&ox9|f(O z?@22O+?iA8H>4m_UctzS1zh4E$HyBX{O*IB+)FgdIsd1tD-WbHZ=VNmt1_WtoU5vSeEL&ZGan{;Y7l z?{9gYbIx;)M`>AE@|p1PA{f@iYS$*WdAOd__IdYY(d!&5iAlVL;whuYq0NFoABm{L zI6lF~u+7ba70Jlw&nD2pQo02we}WRwg)eOO!v>;34smhUoUUzZoI0?Q?b&ShMUo4p^2T9Nx)MkE;aQWaPw-kh)`Wn;Y#*l zK{KZymx-)fRB^-8(=)bypG&Ma_&PF>Mx+$h{O_?vGA(AXJ($f4Mp6O`*%IJ_Ls0hv zBB1F&*g+GrF@hf$y^6!U1{bx(YhrcU+0vHo663IUbT=gLH9 z0?Y)wZzw9x`_3rxQ8kwNf6-WUHtn@7DJfN#E=&E zVz>~;9BwT z*jU)3$>;-@!?5o@%=DQ3D(!#wTe^!O1Z*$C%P8owlYTEBhSfqNnyqT63Y7dVhy4H) z*@{Ee5$7{T%8H6qs>Jb6H7G)#Vfj7>5sX6yFO>~2MaO^ABhT+Mu19CQSC_@@qMMxC z?rN{Bb%2iX7Y1W5MvaA|t1jY~n}qFaA%wzbqmye~_TuhYuB8==1fl~s@tANrvytQ= zzIgFd3<gkkMBzWp*7?WzV()0!!36p;MK$IeR)8S66eUZdXi_umFxtqnNA zXSBym)WJv3;PIG2iv9@=cpXTC6}X}D_>A5=@?K9tB<3nferBDB=_Edh=}vY|dHLlW zxJzPaI|^@obP6{L3ro%eG~dn4%-~%WC*zbVmQv7~1(}Cpo4pO}i z{8oOG@21j5%y~7Y8$w>TF>*r3gkTGIUVjwG%bFfm9 za#%-EMD-K|{ZI$Cu9nVfzOfD4S8T@frl(?ix;0}U6bZ87aNnJ8V`GCOl(`|;#8(gv z4bDkZ2xb6xnRYBwRV9x(eE7vH{E;k$htmrLAMiA?j!rvi)S(cAD0;MtN6s1Ho-`o} z=_AH#@MPM1{N+zK3zqqgDW|5T)i~U!egE<=Av&<2b3-MtK{=tSt?gHmCy;RxU>!VS z2+p%{-HkizJLI_Htio>Z!2~WD9LXA-42IIiZ&e}8E@NR%hD{tX$}oASvWPgHk@0at zZkm^)ru<5;5DTq1@S4tl9_U;-Oh=giM zeQt3M79WoLhd7U300m>Wt=ofAdGxxG2@)Bi5W1;#w|lbzjN$y7nDSJV2N{ zgkpX&03!-??EzBj?(4(A7^Uct9SSKffKB8N5gSALRgC;ynMm#8o|v38hmrs~7naiC z20U?q;a;NQCRx`HSSc2@8=0s)ZD=p?{_|Tx>uKN zqX$t4An`{G5rIVg#xV#@LNL4K>w8x<6Ac?KC{S^z9BUHignjb(aKAQF{Ep04cfB;5 z)b1A!^Q$F8Q&JX#4#n`Ud-1EVJ3U!^*Lzpus3X)A_cHl9xx$nfJ9~+tp^SKxS*J~2 zvk3s0mX^}=B;)nhgK-_$0Yr@Hy*Y#d6h=D#fSu^xbutE7&+v5A)YV6kKH>}-9H&Sn zNdo|ZQB@12BMnjzjEt;4N8cOpR_GoWU}$Ox#sjMZx$$L1 za$$NVIrb#x()H?@hr?QG++AEkf(0Hf+aX?1p|8K6^bbgn<#;F6?*jzMjsDA;F_WaB zq;gnur`j#lNH)T1qHZ}YvD$ncjObk(%r{0m)z!TAhsV%-;QM6l0Z;Rl0RaJAna1DZ z{6yz*6pO98Pq-!|Uh;T8{XsU@?iTb`Xuu1PcAHQFEe@&qEPj}alG7nh+;9vpEXF37 z>Xn;OZ4X^ilya!zf@MzPLHFgiOUHxu?U{7@UAnHD%x-b)kUbkYA}cGKZ`Se0(e^B3 zog64$eAi;NJu$teN0`Shb>@6<)Y-S{yvv$@tK50q_kb`hDJiiC6?RdCY6{D1ph{## zQ9hA+Fj+%GBQgA>zMf&w7^Tx<7plywc7>yQ%dP$(cfXHN;~0e>kwitCch}SwwM3${ ziFZBn;D9lJXPU4#tp9InhP>@O2Rn^xIuz&Hw9TsDs zYYLt6Xjj!-%ufE+rcSgiE=q9oHqkO14<(gHmYDDE!T}-rKnO>QlLqFEjPED*OH?B%+eEW7DTU%Ro zuUCNhAxO!F@SuJS7Ag;MaHr}0;rD%b)uP-VVi-5T)B)WbtF12`11rP?6x7 zcMT0K*WiBuvjdD>05*6!C`kqc;AWAgqw)N5X|dH^rx?f{DxW^J8ZV3Ii7EtDH%B>v zfe&@?g_wLmI-ENU{lucs9Yz67u!>!}cW)JIR+#>T)05T(* zZ1;LGF#-YITVBl-i_TKpMik$OQc?tE2ijXJFq=|{6wxHE+VS)-HUhJi9)+i+KXqWm4@UbBzR(aah^JF zo7n`048kNeTv2L#C3OtOO&{k7EEbw?f^ozd5H@GI{Y~B}ZrzL;vC zQqVh=ZRuqTwfqQxlh{Swy}hBk@B2Bv(*#vqeQDJJ@Jg(03`_Mo5|(f{?6z@v({Jl@ z>@STr} zwIe70DS^u->Fv$0aGu;#npcp4;XJsu5s_1q>K^jXe^XPLSe8h5x4rLYt-BF72GuIi zi~P1&c*e@V$W)&^)0QhLh6R2ft$Vg_rDNcBN5}eU_=psuNM_pEtrv=N2H|rPE$a}( z)neNdl3dsht-wVb!D+Q&e`shZRL!>Ui5D*x;W8}5Be)GpxTPZq`;*$z{5=o;e$$?8 z`fXR#2)=N+wsss58h$wj zWhlfN#<+K&=%=m^KpURVY-%w$66hv{S>6Zh<727-plE>>rtDYPoY3q%nx=4+sm35l z(AXlf(k0Z*s00~_YmfE8T26|LW5C2ixlPWNPGv{AY2Xlrq^$2(1I}c}s`P$xGeqgX z7ZO|&%#$1r+o2m0uGXIk4zN{TRy*9*w_j}XEu1tK2uOwD6_P9hs!n@2KPCbq^N!^)!ADJbW+AEqur51Ed8I6&EP)7&4qKD{Afw-9QbuScTI(Au~0r!{bR-~ zDLe<=FCx|XA<3&kK-xs!AJdzhi#MQkfx|r}m;mwdo;Y!0i-Snm3Y5YS$DccfM}S(| zgWUk&iV}sLwLrvGHZ%-xXGkgy{ergx+ziLPuAW{s>;%+ot}R2EjR%e)Ir*fg-NVBI z#XGEpis!PRfOd3b)JmqNrF9Pu^5xjK_2d)-9r~G>16Bj1%BcyTxOq(dA?{76$n z8f9zyi<}BtxHHb?My!UTvCl@I=JqU!foH}y9`fc||F*EF$*II=z^7#QN>K;XJOOsr z)Pj(Mxx>60{2p92KI~Sr@i@P~`(c){%q_Qk8PVZ;yZtWs^_3o!sQ2!WZ!(Q*;)`0c zV`;q9jgwNhGs30p22672$uBvnyYrP+-1ViGg||Pr>KP~7Yq3`5F=EFdMH#3~Trwmw zs7$c?VNucT^C<6r7J8@HQo`ctFCz9_+9+U=6U+Pw*>kRmxTYD}a1@r0x4*c1wPVe3 zn*kGcNmO0?^VY+I<%d53W5|1e`0KV8zL(4U`eM}Yr}w3wuqu)NnT9H}`KEeOeAOCzYLbx`0ThpU!Ve#Ys9c~{4a`w3m@Ef>&uDSMcVv86@)O2HL(`xVs z%IsZEUz1#ul5QPZ;Ls1hYrydWXO%lx_3ykenwVfN>^L&(v1SB26=8>+poM(#?{ zdyvl+f9o!6p(?#nPs&yMc;9uZ=&Ox6{K5K=hI+|SqwaP_$U#XZ!#Ut9)KDA&eL}6a zNFvR1n~=S7DCluet&sifIbx^}2?II-KrdG+DmpE-=w~ixR6oOAhTcCIrMDyta9uD2 zFhG5mB-9`6s^lgINM^O2uuVv%*z_8#dUQO914cD6Vjv_Tk>vL@coQgmV7_Qp(Ikl! zMECnTJvb!;E!~+K7q+N1NNj2^ktVJpwUb=ih-47LP92?LYh*#BZzboj*RcdE^cyf|GMyAx5BKs#WpE2n7P=08{AzpkbZ?!D4Ti|o6fi`u=@+Tgw1X$dX(gTdhS(8nAo8P0$gyvyq@BQ$j92r_zbiJ6bM3#J4>8=8%I?bq3u>0uVxr~z zyH0{R_KX~RN6pNg5lQxHi9NAlBm}uuCMo6Wkekn zs5TriIFDOeA ze1G)HE7pGfyhdj5>b}E8Y5Nc8v)JrbwT`Wf{6OUr6FWZ{`?%HkCC7(5B=CHU(P=ds zb_f>eC`vE0o;$l7vP1b+R|9G;f!Hlv6g5l3JeaIyW}B*+ULD}M8dB{ul`kJB z|KrVTmQq7ibE{5C$j}zykcq{XxmpnK!s8$t?_9brvfo4Xqvh9SfeR>$9;+?J?Qp0@ z+M*MPOd*Fy6xlzgm@!vNE%!6mu_NArD@6WmNGp){YLSA=03R+Jl)9SBa*sL<*r_tq zX8BH2o{hqMLb3*xz#XXW(QZJsG#3Lg;jk2ci{zuhqioX?D1c+6_x!r0 zii%esDnGg4S});hJ!d`YLkjB~y`)_Wv<=(dJlsqeand30ANSmV42Rs!p6$Hq-{OCIfH^jJnTi%w_2e;VCaqop(V zNe}kkdB$@<^H$qn@mB#fo?Rb+`?^g6pOyy@ME9SWUS z8@o_FLgposR33L$qx*L5JHMY~u_=$qNJ+p%^L9-5RH6ovTDTPEJ*|HS>Qlzh)YR15 zfMh#f4T;wZ9v)c-uz#Xz;Fk5O=R@LI^}QK5)hWgZc4z-yTFNU*y3U3u8he$-JQHON zV-Vr(`x1jj7OHFrdBp<)1A|^3>hq=8C(H?{_OKjorGC z3#bCwFOZG4j?QcMpbnNf=ANfUy6MfGGiN2VaZ<{42E6sG_xl_fbbww5&jEppeB`7T zy?6~B+tPJVh(Kuobr%3RO3pcnXPsrAD-|H#H2)GrVS^}^AWxtE+Va2!STvbV;?l~m z*xrHui4Cfh6uJ7tNoLZau(V;D6wJz+{>V?n=k)J5d+`?7>jYK6HK~!nc@O}6W zRVPu>_-ZXsn0dbm<~a5PFO0%yu4#o1mnwX~Kp_x9h=lVcka|nkInHnR<%=K21jv~? zn$HYt1cG7mBE)CgmpYF#Z>9Mz5floNCZVS$2Dcet8X8#l?khcyFlLQr8yY)ew@H|d zz&5|zy&jY#iPl65DFnb_@;-qmfJER7y3xqDZL)ts=?f;`9H=BP0jbm<2XK{w(!pW{ z1Ad*ecWw8n2G>=Xb)MC5VRu5NjKLC9gfM`rB#aSytGu>wlSS#6CWk&<3{Wd`sW#|y zh$QPUMo1z@9v&V*jnWF3c_8T|@Sg-1EhGVu87>-;db^(oJ(m2Vgisgkb}9q^U!d{c b?)$n{*_#dfkAH09!XFDW8`E1R&LRH;_DSR> literal 0 HcmV?d00001 diff --git a/src/parameterizations/vertical/_Internal_tides.dox b/src/parameterizations/vertical/_Internal_tides.dox index bf6e16ea5c..a07663d4a1 100644 --- a/src/parameterizations/vertical/_Internal_tides.dox +++ b/src/parameterizations/vertical/_Internal_tides.dox @@ -4,7 +4,7 @@ Two parameterizations of vertical mixing due to internal tides are available with the option INT_TIDE_DISSIPATION. The first is that of \cite st_laurent2002 while the second is that of \cite polzin2009. Choose between them with the INT_TIDE_PROFILE option. There are other relevant -paramters which can be seen in MOM_parameter_doc.all once the main tidal +parameters which can be seen in MOM_parameter_doc.all once the main tidal dissipation switch is turned on. \section section_st_laurent St Laurent et al. @@ -69,7 +69,7 @@ case the maximum of all the contributions is used. The vertical diffusion profile of \cite polzin2009 is a WKB-stretched algebraic decay profile. It is based on a radiation balance equation, -which links the dissipation profile associtated with internal breaking to +which links the dissipation profile associated with internal breaking to the finescale internal wave shear producing that dissipation. The vertical profile of internal-tide driven energy dissipation can then vary in time and space, and evolve in a changing climate (\cite melet2012). \cite melet2012 @@ -135,9 +135,9 @@ at the ocean floor, so that in both formulations: \int_{0}^{H} \epsilon (z) dz = \frac{qE}{\rho} . \f] -Whereas \cite polzin2009 assumed tthat the total dissipation was locally in balance with the +Whereas \cite polzin2009 assumed that the total dissipation was locally in balance with the barotropic to baroclinic energy conversion rate \f$(q=1)\f$, here we use the \cite simmons2004 value -of \f$q=1/3\f$ to retain as much consistency as passible between both parameterizations. +of \f$q=1/3\f$ to retain as much consistency as possible between both parameterizations. \subsection subsection_vertical_decay_scale Vertical decay-scale reformulation @@ -221,7 +221,7 @@ the implementation in MOM6, it is required that you provide an estimate of the TKE loss due to the Lee waves which is then applied with either the St. Laurent or the Polzin vertical profile. -IS THERE A SCRIPT to produce this somewhere or what??? +\todo Is there a script to produce this somewhere or what??? */ diff --git a/src/parameterizations/vertical/_V_diffusivity.dox b/src/parameterizations/vertical/_V_diffusivity.dox index 1d79f58997..3e0ca7d4e1 100644 --- a/src/parameterizations/vertical/_V_diffusivity.dox +++ b/src/parameterizations/vertical/_V_diffusivity.dox @@ -3,8 +3,8 @@ Sets the interior vertical diffusion of scalars due to the following processes: -# Shear-driven mixing: two options, \cite jackson2008 and KPP interior; --# Background mixing via CVMix (Bryan-Lewis profile) or the scheme described by - \cite harrison2008. +-# Background mixing via CVMix (Bryan-Lewis profile), the scheme described by + \cite harrison2008, or that in \cite danabasoglu2012. -# Double-diffusion, old method and new method via CVMix; -# Tidal mixing: many options available, see \ref Internal_Tidal_Mixing. @@ -50,11 +50,31 @@ parameterization of \cite large1994 is as follows, where the diffusivity \f$\kap is given by \f[ - \kappa = \kappa_0 \left[ 1 - \min \left( 1, \frac{\mbox{Ri}}{\mbox{Ri}_c} \right) ^2 \right] ^3 , + \kappa = \kappa_0 \left[ 1 - \min \left( 1, \frac{\mbox{Ri}}{\mbox{Ri}_c} \right) ^2 \right] ^3 ,\ \f] with \f$\kappa_0 = 5 \times 10^{-3}\, \mbox{m}^2 \,\mbox{s}^{-1}\f$ and \f$\mbox{Ri}_c = 0.7\f$. +One can instead select the \cite pacanowski1981 scheme within CVMix. Unlike +the \cite large1994 scheme, they propose that the\ vertical shear +viscosity \f$\nu_{\mbox{shear}}\nf$ be different from the vertical shear +diffusivity \f$\kappa_{\mbox{shear}}\f$. For gravitationally stable +profiles (i.e., \f$N^2 > 0\f$), they chose + +\f[ + \nu_{\mbox{shear}} = \frac{\nu_0}{(1 + a \mbox{Ri})^n} +\f] + +\f[ + \kappa_{\mbox{shear}} = \frac{\nu_0}{(1 + a \mbox{Ri})^{n+1}} +\f] + +where \f$\nu_0\f$, \f$a\f$ and \f$n\f$ are adjustable parameters. Common settings are \f$a = 5\f$ +and \f$n = 2\f$. + +For both CVMix shear mixing schemes, the mixing coefficients are set to +a large value for gravitationally unstable profiles. + \subsection subsection_kappa_shear Shear-driven mixing in Jackson While the above parameterization works well enough in the equatorial @@ -117,10 +137,10 @@ that the TKE reaches a quasi-steady state faster than the flow is evolving and faster than it can be affected by mean-flow advection so that \f$DQ/Dt = 0\f$. Since this parameterization is meant to be used in climate models with low horizontal resolution and large time steps compared to the -mixing time scales, this is a reasonable assumtion. The most tenuous +mixing time scales, this is a reasonable assumption. The most tenuous assumption is in the form of the dissipation \f$\epsilon = Q(C_N N + c_S S)\f$ (where \f$c_N\f$ and \f$c_S\f$ are to be determined), -which is assumed to be dependent on the buoyancy frequeny (through loss +which is assumed to be dependent on the buoyancy frequency (through loss of energy to internal waves) and the velocity shear (through the energy cascade to smaller scales). @@ -138,7 +158,7 @@ diffusivity, the second term as a source, and the final two as sinks. This equation with \eqref{eq:Jackson_11} are simple enough to solve quickly using an iterative technique. -We also need boundary contitions for \eqref{eq:Jackson_10} +We also need boundary conditions for \eqref{eq:Jackson_10} and \eqref{eq:Jackson_11}. For the turbulent diffusivity we use \f$\kappa = 0\f$ since our diffusivity is numerically defined on layer interfaces. This ensures that there is no turbulent flux across @@ -189,7 +209,7 @@ The background vertical mixing in \cite bryan1979 is of the form: \kappa = C_1 + C_2 \mbox{atan} [ C_3 ( |z| - C_4 )] \f] -where the contants are runtime parameters as shown here: +where the constants are runtime parameters as shown here: @@ -227,7 +247,10 @@ the diffusivity is where \f$H_t = 2500\, \mbox{m}\f$, \f$\delta_t = 222\, \mbox{m}\f$, and \f$\kappa_d\f$ is the deep ocean diffusivity of \f$10^{-4}\, \mbox{m}^2 -\, \mbox{s}^{-1}\f$. +\, \mbox{s}^{-1}\f$. Note that this is the vertical structure described +in \cite harrison2008, but that isn't what is in the code. Instead, the surface +value is propagated down, with the assumption that the tidal mixing parameterization +will provide the deep mixing: \ref Internal_Tidal_Mixing. There is also a "new" Henyey version, taking into account the effect of stratification on TKE dissipation, @@ -248,6 +271,14 @@ The original version concentrates buoyancy work in regions of strong stratificat \subsection subsection_danabasoglu_back Danabasoglu background mixing +The shape of the \cite danabasoglu background mixing has a uniform background value, with a dip +at the equator and a bump at \f$\pm 30^{\circ}$ degrees latitude. The form is shown in this figure + +\image html background_varying.png "Form of the vertically uniform background mixing in \cite danabasoglu2012. The values are symmetric about the equator." +\imagelatex{background_varying.png,Form of the vertically uniform background mixing in \cite danabasoglu2012. The values are symmetric about the equator.,\includegraphics[width=\textwidth\,height=\textheight/2\,keepaspectratio=true]} + +Some parameters of this curve are set in the input file, some are hard-coded in calculate_bkgnd_mixing. + \section section_Double_Diff Double Diffusion */
Bryan Lewis parameters