From 1895a84fd16698b394da07edd6d5a0c731f7d715 Mon Sep 17 00:00:00 2001 From: daw538 Date: Wed, 30 Jun 2021 09:51:04 +0100 Subject: [PATCH 1/4] Vertical coordinate docs --- docs/source/modules/index.rst | 3 +- docs/source/modules/model_levels.png | Bin 0 -> 72777 bytes docs/source/modules/vert_coordinate.rst | 94 ++++++++++++++++++++++++ 3 files changed, 96 insertions(+), 1 deletion(-) create mode 100644 docs/source/modules/model_levels.png create mode 100644 docs/source/modules/vert_coordinate.rst diff --git a/docs/source/modules/index.rst b/docs/source/modules/index.rst index a53a154cc..f44b75efe 100644 --- a/docs/source/modules/index.rst +++ b/docs/source/modules/index.rst @@ -14,4 +14,5 @@ Components of Isca idealised_moist_phys two_stream_gray_rad convection_simple_betts_miller - socrates \ No newline at end of file + socrates + vert_coordinate \ No newline at end of file diff --git a/docs/source/modules/model_levels.png b/docs/source/modules/model_levels.png new file mode 100644 index 0000000000000000000000000000000000000000..a0909fb68923868da87c0f18e0a521f96e7e2aaf GIT binary patch literal 72777 zcmc$_WmJ^W+c!FNmw*T&F-Vt$(o#c6cekV@#I@JQdyaMRJ$B@6YF5x4m& zNkpRB-@E_BcIt=I?uQ7~eyxb=)>iM`$Ioo-RI619d$Uou_dj_bT;J;~c#C88LGH$G zgm0Qd<8LU(1ikm_r}gw{pz#0sx|K~}`iK2LFX28mFnVp+e|>F)N@_Su{r4+h{{PQ! zZG54_l&}>i3=D6ff<=4Rv6|J-79b@Z^nA9j0*87LM=qv7SSNU%s1#(npm zC~-kKkCn&Q+#&(40tz+8P#t7Pg^}Q{ zVEa?}=kf_E5u!RkmUMSoFauXuFeTDyFbUAai|9p-2!lJ(#4#VE(y9`oDvJ}&S(ZC6 zVSeKgCwx>>*0)pp+wVnNRy$wHuakN*)KxeQCKluqc3?Iu<~PwmnqW&q+hQ6_deJi@ zG;!;6%?EyQ!Ua@u>x}x2EE>%5rFr*o8ced9Bzc^`DlNib|1$4cX3Rh@s_B6*RVn zq!m(sdR;$a|EG{sq@p19b|MG`5);-+blJK*N)a**p}*(vYAr}vZ4x!oll1yBv}~u} z&bAQot0^m5e+5WqkPyc|$QaTLni#P^7laX#chb}t0*TtF>0tA>hGf)^&A;Appo36A z9wE(`&p$Bx`|1uqMKT!hVwi()AnP36-)&q7_d=66b!&!4ym#1Do}mBf&fvWq)B|IXv09u`P+Z3k@{#2%R+Cwqx$NNiKbd`BN-QXjK7XMGnCCtQOpU^tR> z9ox-o5vt-6`o}Nic2{wb8hwZ9hZb{$h&aZYDgGKKZ+HFkU26ELFOQ`gK?ec zkO;Un7m)>FET_N7<)Pt{0*n?$g0ep`+u>&_PFRPefa-3J62?T{MGm>v-MS+=E}JH6 zp%zO}j9n8qsyVB-5KGvmhNuqS@%l36a>LBuz~0u9=m`YOCb1(aj4q`FdVBtXZJ?NN zH#JhwVt&ICCmgAbdWw|I4IvOG%)}0y*d+dUet)Kk6l~dg&rv0Hpx<>hx4Oegkdfv1fLVyb zl8}Q=!zRh0?8r?W>WVyUY4-Q_U^FSJX?8X)yj891BZ)@kcsn(y+72`|NBzWSnpNF0 zO8@6pxWu-H#b=XHkl5&QEt>abU2(g(#0g!yhPGH7l_#Ex`H2(CLH;6j84;qQKn6|n zi6Asp4tqD)>o?QN7$&l?xbdb!F-T~bwK>&%CCUpD)b7V@KI%2UWsqdiR5E7#1dNZp zsIpWYh+>sDdHvl2X%bK-o{7bmodrN{;RAng40Bcox z+B*jB;5E$0iC7SV)Xu?%3S>xSH2?XR0Wl;vx8Vp8d$Vz<2;}||YEEpIHFn@(DGlZv zXk;fLxR88VfNEGBR3+6XYaJeIT%YyO)G5KQiC^{(T`Yo~bMoDM5 zuP_5m(%l>^(8Mn{?>_&z5Qf1zNMR&cEBO9*E9jC-e+7|b2e5yujK>901k&&_^;#Gq z7qManmUoG#uy$ok#-mXdUT4~sfTzRhB${g%?t`L)oWhDJ(OmYSvgeyce^k(M6&i(GJ4e`yFFZ#R>G07x!*)xtm1@mT_WASumow? z5#&4a3=Yi`yIFQ??hslN#?xtCx2n*eqcnWiLB<}Qp7O8&NQ=cU&PNMx)1%T4M!Bja zZ?-SWzWZ4%Uc3$9dA)!t)ZVH5U;9XPAeHfA*=~O0pvLD3Xm@CQEe(*{u7hnONK2O2 zYT3azZb))>0BwQE!}*90!Ny(f@SktdWz$YYa;z zB~n#e*wlo6X=q-YEV<|2!#>B;i$;x6% zNJtP85e3Z5==ApX(!6|03f_v5eXgRa8W0|iT~}9^U0RBRf`XEgn%cdvK$(=3R8>{w zA0A$0NlL#qNd}(n)uZ4Lq`^V|`j=JC^e;SK#H3(GP~k1KLgigqCgKlYi@Pa3HunBXeAVunv-w54ixaMI!6vC5ViS ze7L_Tba{2P$#k&PE@@yuIh4fy*?DsaaVg)we}D1fg@0ru9`DPSIRyou;^^fiSQ4Q; zJj7vPVesf^850wl%-0{oV`2!~yf5%qS67*|D=6*k>{ObaEu^ZetHrO5H&tuQ2XlQ% zM1{X;QxuT&zm|eIt1XN=;mLSe8N+i3I=XQI2TaUM(D6=iqWmfA*y?K6Otsj2Ls!H5GJa|ybcHPZ zHZL@5YX~Vz)i*oK?@1=m2A0noAxxFawAr zfPetRvrE&@Qw*UT^E)Dx{xQX>7&n63?(JXlmY@w5s2DDL^TJ96^ERx`*w$ zk2_V`Ki<`iR>0>m=pcam_>7m3ew-7AuSU7H$BVC)U>p;n0R1UB7Afi#G3m{cWH z<&xOK)bizW5YI|oJ?PUXbY5Ow#4~*M?AhP#ZE<5`>g((47reY0By@Ci;DW`qwfWpP z?{CBqY$ZHqZPB`*kXSDH&+3PlFO8DZu&r&XIhne4h#y_%vBx7ZGy7L?#3b};iMzQq zuuza65qe1=o7j!WWv5YO=b^NDqG%L<3%63hCD`lq?CH6@zFw$u-5nURhM`V0aXGFc zo7fHy%{9F5iXbnumV!afs@?=~j>prkg`z;B-?(YurA&^f-4kbY8k9JuU zaz}5lN4viaN;GkB(7HA=yY!vxqxhFNLU9+;O$NBQj8(c!A(a|VW7`^A8q&%s!n|<_tuCAN2^EMFt10V9w zN|p@zA2c;)>LuhqPg|HDIxFBl9868kuo%=b=-1}u`U|?vcr8}Zver!&qIjP>N*){# zuo+<1*FRIJddwdcG06)b{LGi#kYO;XCX5UyQq}R)-y69 zV&unnONai=;aR^E7VgyDAmd_!n{|g~95{ zka?FCQ|$3WY+KtyXZ(Z1&&kQy%*-fIrqAi=nHQMX*4CU|v7vNyM5Lr40|T-!X=$i{ z06983x`13Nw}nc90y(%A34Q&ib8~YO6BGUg1q@TaeyO&&+7uNRi)(6Xn!eXkRTUqM zqepaQO1UNju4FQSbIx&T6w;~j%ar&W8Ck8*wm$S{8v;uutt! z%xioRaCqpr(CmWIXutR^kyTrzM1$FN`zM-V>j$afBzDyjjp9lVt1NK^T-Ma>Z96pk z*|-N|dik39{$;-tZ&FgU-8nS71U>R%kgsL3brH@;kcIW!!&OhPZRxGw?mTmr; zn21^Oi#}GTSkmUG(MPJ)N7>YpXx2QugKC>-XlP`Yl$5C0AqD`bg9HvkJyR(Sb#=+@ zpZPfz75Hv$Zhcj)EiHnGQJ<58#LmucZDS*@u1*X}iOI0#becLUI&xEK_U^{Nl(iFQ z{kv1KM%c2H%LxHK4h~LGX(?-9P!M9S*46PMicMY~Cy_<-o3iK~Xd%M{>^ZcwxAyfm zJ$pS7)g-f3BANl#BN``ZgP9kZa&qcGG*Vb5i( zwm*glcPh=^3{%Uyy4G0_k`ggXguCTosalGXrz)Uxa&j(o`cb+&EOnN}tnv}yy!d-09wZO*I9fyF zFl-6_OvVZD0m}RN2B|`&jENjcj08^O`+1RIYBE)$PIvgsKz3eFo{*m2!eC9!nv%K# z?y%luK{FK1%N(T1%=?~Xab{7-^i&c zjXK+D{AWTRV4#SyFFJ;m-m$nJBBKOY4Euk)@%Ou56|8?hoMbL9&lM!MN|@ z@*ID1NnHWh5R!_ptfGC1}VKt)I|Yn*nZ=y4iSRGFf1-^Z4=OPdn3<>;Hb0HhG>f zfFk&jD}}>rJppNGXeiRHWf)3eK5?U-hzRgRVhDC&j2e!A1C}j#`X=`iO*-WBm2)appdG%f;MRi-_M!wQ1=&Q74Bv3j4&;}(ZleM+A9sU`8 z+U9#lvA@3$x}w6CkL(MC9&gw`xi7Z*?%1e_8R^M@)I2y6eDZ*(9cAwdOjP0DEa-LI0> z!Yp42u)HSE$4&JKWul$MzI2M}#}p_GR?m))6OeHl_5Lc=K3eHUf%qpTlA?->nVCVs z)&l@K4I{Cj=c-R8$mDH5wKcY(O%?O|nNu z6r{vr7ebYwu}sKtFr{Mw(ER(?S}G#ngcf2nJ=*0KtD8G8GB4b4DFATn}{diO3aDKA2U@GjL6ng2fuwv57qB&z` zV|F~)bHkmAOF$*u6-y@ruvGM1z3tn97}~!-@?^&F3n-g2Hx9~xjG5a{p$yQbyU^3zv%LE8EuNDlC{G$zn#)5lik+C1yz2e%W7Z#@c^yyRhEFAPY zl4h>*@^Vuf8(n_Z%~tmX6M`pQ-L4{DRorRizFT7%Jr#znD}VlcTK2ucL$sfr-{kxkZ7*zMW0@e`pBY=EMA6wc_3nVY^v!KEJ4W!#wW;8rUp z0VD$x6Vt!4k`ti$mF;aTWo2a#4-YXJ8JQmfG0C74_Kb||dJi4xF0fvoZD_Js^cA_< zXhm$KvbD7fr0}tX-|4byTP8ncWo_W5;_s6C*Aq%PGF?c3@fe?KghN2Y>+pW(x6|*> zS0@6p?!Vcz@W<0vea4rTd_-_QbH1}k7&P}d-oP457YK@HRL)gM<$m(Q1cHuFnbJt@ zEsd^HLShBct%Dz&eaACd5yN-Z)V8v4{Eln#i=}}nG_+P;oCDw|v++!knR;8Yh{#AT z=MCA@{ly!}b#hE4{C{tZht5b#B-#9*1V3>GCzsl?Sb4nd6+BN7+;DAuJ*TiRnU1fZErDJC zRsXIOn8(Hy$Ic;IQX!(ZZ~OnuqywxMgh62Tt3>k!A7B4W^;=L;(dpr?=FhFJ{!TG( zZQz-(metns{H`$kEabj#@;!naJfJe?Cz{YzSDDVc?;Rb2f?buJzRfo`FgT{TxVXpJ zKT^cN;d5A88U4a{cz!M|Cx^|VRaW#@13>~UhLb|TX!?eSSiZhK(p-;W+fuFIR`*$D zedwwkSGoR)iklwA){73`44%9vPaOaF>(zTOEssaD=O%NWIDRPCODC7oPGB1)DxLar)!)|Fpcyokc7|JSPRwm_l9ax z-)EQ6LujI$)7`b(?qVw$xH%XAUO(Es-Rm5dslY7Da_TpJQk?MM@c40h7?eH)rP)mV<+wo4;=;g*~cthFJpxF*k>tym!A$Tt8d+l7KAO=Z)efP4#NeEVf?0TqYemAfx<#_ia!eMP+4zvwy!N z5p%u%ohFbAx`29RF_aLHmNs|6MFEXJ;!|Q!|D&%gYa(A+z+e@}^>Xuq^L1u29o?0F zVfkoA2sY`=T2}f92lMKxPgwzjLu11qmxBDGRwEsHZNlIJ$T{)^H1{}=W(xg98vOeo zmk7juHxGh`pCxY(#L^{bgXcVqw)rdJVF)F4LtOD#jKPybvH53iMTJpYVF@m_S!d_< z*=%673**h|)%CGC<^jN|PN#oW+4Sl#fB*iiQ)7-r@$%i*h3+4!HqD2JCLz3h^t|km zlT>HBgVcyZRB8-2}F+$tn2`u9Q!9<{>1TLrU<1Q$Hh1g z4{M2OpD-;HL0xHS#kq9-ebkUV8Q<(nvT?$`UM6PIRT`vw-rc?SyFHRQ@?>g0@mu~5 zDKOR^x|jeK@_}ze@mN1jvTg(-qI&{Vq$rakrb}ZBPP;rY0%6ZP1*%&FAXTc$R4B zmoHDj&>$os30?t%wZC5yXfTqJbTU@TnYB8Kig^G2{lk9i50oLZKVa;ZWCPkF=Uzx1 zFTvWZ3^_%141oVLFZX@}yPCazI6zZ@ZHV+r*ESJ99oEwP`0-pw$TvMye{sQM&bawN zk7K}`duD$AZFe9V5NvvZh&j>XX1j_1g|9Z(ppnD4e*~91Re>vS7GQrGOk{xCOy0kb z1a-)YVxpAwXG-EYwR0#zTv{6amBUg5Abe)KGu0^qE)2ZRwfuwT1`c?5R2=C%q3rTN z%b=(gy*DX6ilMbb!^K)P-TzD`kYf2A7eCi_$)Fe&26Gfysx^8rfgkHB#DMvAvH#4d z-76%S(-^V#oJMV307t3j%aKl%>yzP8@B*fUr_%aiPuWEVo4|^C<&{wY>gT{Y((D?uuUpN4NCTqQ6&BJTwdv=Ohe3cX9 zKtGTqVoMV-(ReM96c?A9(nbpq1o$_H7Z;(0g^U0&o4k7$4jN1jxHjvBCN2tIYsjc@ z!LxgO=C(L_SR|^rN2cPS8(r?)73H}foj(*x9X4bAm@zPrKTwD^;0|-T0!^!I#k7)$ zZYH!$>a#()z-+ASbN-Llyh$8}*g$|}GiW00ygmN>Gxy|>7Q=d>bG*l4a|jeG5(Yjc zkaOp0;s+=@N`d35<4ul1po#CDng}P@8wjHDoXvkUlOBf zFICI+>c8DM3>WVCUd5r5I1j2p;|CpgXShZ(MV{8%PJd6~dV4(L&U%qyT>p-gO4K(3 zv}`kg5e3*=kBi6TlgvZB)_TAFYLAi4S3@YseSK)kS+D!R6FmYO>GpI|*&1x0T!0kL zzx>|8iEADC!&m8+gltEU%?u%f^$`N%TDyAr1GwfXBtb4GiF3rz>uZN`xYQ_wOqn&Y0eL-sWs*GQ_I_ z3BI|hX(9xN%o?Ci-nZXz|1Gxi{i-z10t+oGFaP=7bcJHO_eFt#%b$tvAdJDm!7zZH zIt`?_xi8&kUgVZ)pDgSF@-Wrp^!D9EHqv6J-$OhHLc1){uB5uV-l+uG^k@}dTpT9! z;qZ9A@wS`^s-~{KeSddT_wEKbee>Zpi;yBWZe$fxS+JwG;_SeToe z+}~dDJ8xh{Q3}LOy1nt_O1g2`@<7ekD284{GeQx2Y;X_@QTAYdtbjG+wiwFRtA7VZ zKR4+4!Qcv5k=G(3B0B88r+quw-p2-BGOA+?9R1^3U&k`KI&5Qy=2cpjgc+0|^-GH5 zP;T1GW>dXRDxH=Dl&-aP0RpRJj&!*1-`bN9I^hW-M3%fV(_HbT^;t+W@yw zkHw>}R8v7G34205pXWzF+Gud2KXL}Rn{RM2`vNf{v?gf7mKxYiYbQk23zUcYA0 z`T@yd`M4~InDO|ce>scs6X@ANE=x^LD1exJ4@ARxVL6i0w!`nvp*;5DgI!0*gH*r* z;~Y@TjEsyB>;gzr$S5eJsj7fA{IA_PGJm z7<(_?8&6CdI&>MlrPfBFKopx>sKb}$aS>%=fyOUXOR5Bxz8oEo|Bc6B+^#y<0G0Cc zZ*dGiKTdCvlmA@q+8^b#Gy})f>hOKYU?zlG7L4j&DzPIQy)WAfT#TiOm zUx`@n7g+uiPMU68W2M}2U#u=1V&TfmDrGzhwsBGTj7iiy(wWx`?#~}J2nz8-q@^+O z@uh*S0t7F#-Q~`#JEHhFomB7*W3ES~Q$$#p>Q{!Ak2kM!b8`_^sYaQOYHXk)kO}kU z5+#5lORD<<==M24Dl#-QOi4@oHZ?^8>;>cRKhHO?FjkwJohp7eyEN~Mcpu;VtY>xC zO3!GXoM`34`1Y%~M{D19nxFHq)jd7)Z()|lvRvX!=-gt)@{eqqh9q6Rn!Yw?3oVwV zj!oe&NJzIv$?U@i@6P`1@@*dwfSp{n>DCucdg}Wpej(f-)k72~ylsH)w6@SHo|XO` zfCbPjWIajA^_{|hHwHf^J&|$r^~Y+ls|+@XP24~`ilBIj4W?LsnQpD~5*-l`_0mg0 zIp4psyg*ADZ0fPG{5~uHyK;Q>%l-MDJ|H8cq7Sz`uk066T7RjtY&rnpR6q2i4@6nH zLrH6FcP(6j0l3oXxdXMiU%!53`YAOuDLW*jb&RilKBn`*P@}^j#y9QbO`GQlD5F-K zE_KKWWM)vlqn}+ztj+|<`_0S@p{Lg)5pQXYMvltbS~|Vfq~y^Ncwo@1afK$U5ufw( z9~FjeW%*JGF=DQ($g0U4Z54igxp^|RZ7*sT;-jcG&tPrKAviwiNslY%>pzdn_f9u| z{8DLa`F-bi#ikdLlL&A)yyj;P1&=a&y|2@~S-hw_ZimnJ45mfzX^xL*kIyzOpt1C@+_rOCunVT4B^Kr(UR(pXGZU1|$ocnJQD)qctfxMa9vFm812jaamSbIT{SxBFoXVsZKvXuiHy5 zYC1ov{wN|T8Fe5Koc;X~1iTMLpsuyCzf82gcD)}R8EGszIy&0e*nk0LJQs*Z6F>9i zZ9t3SeLu&W_QDi7lAPxgXr6z$~4rvr(rI@#1~P@8}O8 z=J6@>z4=$Q^W>2cjcb+GW+KCu-r=WwkJ^0|Nyz%|6JW51d-?E-4$?5n%o&t!l?x3!G*N;3pLTz@r8L&l`Yn0BcwJ{Foie$XFPN z()j#QZjNL#?Hh%W9}p@5=U`ksyzRr)o;E0h^sG)nff}$Q%Z)oT!H;bKgmk^Vbi&5M z$^j;n*Zr*kuwKJJ?w}YbIi~ykBBAXR3lU zp!Kd>W6$Qm{65>8PX=)hd5{Ds1cJ+tW|u7{0fDQ}#*LWF?>tF1>uJ){b&3<@!qRw6 zbrlrvto!!QcBVC2KfM13@ZNs6W>^@^?@sFZ!B=1Go=?m|N4mg7{TYl!3?$Zf8#_BU zHe22oQJn%gIWYHw}^yCDRSNj9zQWP8KcueVW`U7Cu z;^K4dNggUKUh(s%fZ~!wlskiBS`6^uU)s&Y6ewpY0lIX5iMgLmlgne}anSDf*(gT z4x8J(JU^ocD+2@)2HzI^~m_0&c84vV&c*N zo7blOJd1^oUjUFtsMp25BH-7BzIWF@K+6OdEVLDM?2;)uo<}e{Q)7XXnvyc%b_MU7 zhB45EM_ZM068&iNgy}aNIuluaxWpz=8bot^>fdvd|6gQaDst^XF^~N5V`=RQgB*a) z)Br;)wwtXnH!%^*5OOb|mrt4~*KcInR$qi>if(;BSQN6LsIdbP1GSz{c-t4te%kml z9bQtrrh#Tzw$l}{)zw_@mpcq+tA@;BkrYHBA=N(@%AZz_Q&73NxX{_!+5**vg4*Tm z9scl71X}Oy4F)%xTq~iRRg=#ul2ItFD;Fd!@ofZ#Ek@nDF~28&LRd z&Ua_OR#qy1ygAQDd@{JVzd(UM1vd~M5P*dD^i{;|-JMIO!@W<^=0y`2wy9}pKT9-A zI}BSyp&z>8(OWlVhF^pTGBaURRLd3k#?e4rIIBstg#1R)G5nt$d=zTS3nCu>;SN^g zRB6JV-1QF2nV|W;F*eQu84P(~w-gdHDop_9jS%Uu?kCuIm-m?>@uEd#ES=8*Nltrd ztF@2EPlVxob$NKl>o7KS-1%@EObXkT@wgY}JAV0&iCwSN{p5qTb_4Q9C?Kf7|EH(m zvBU+4FbwqRVyp4Y|IkW6BAndbBe-wH##vrW=VUaM6Jd_V&~KQhmCC+61w~qEshDZA z*BQ$l;QUj-ZuYu4f4?te;l{N2Oh(}R{M_8yS{4{a4}WbzIEC3T0)Azo2d@s-iX&EK?Ad`*2?80r;!G_g|8?FUh~jU|5w7v@pG#%w$J;P(2_V~ z%TOO>7U=PVOtsRRU{0D^7;3S`wRn`Gr|&2tXu`-y$Hheils8#fSuya?!89Z!CJrhp zVuH!Y&;g+c@I@7EZJF5`3sZZ0e9$m~-kJ@VVpT)Kr~3M{WJRj1gkEBR1G<}EjE?Xt zYEl6y7Y!I0(trY01Nsee1xKI!A(<&qVagX=-C(h+T(*=DPOz4iR#zCYmj{(jKtHvcKet^}Xv%OO!wiPgVLA7WyQV03%JgMPWg zskC12M7M>_j*76SUACg2t8GA%K0P}-np74A#x*vRafdJ>Ao%v}+b2(-ZVhzmwzzIr z@|lV$`FeQXtqm;a#So4;85TdP@@eMYTrN`W%99b!FE3`E=n&a@Q=ze4x*_|)$*CJe zM^@L@gTO`Sfa-4U6moJ(|8{u3N6CEIU^g2G{QK3lzRyK;Sz^)*>&#tyyZ-XY{Lhxc zNd=A%T3TBX>bOq3*GsTzD%+G?dAClCnht|*Kr`+EZXzP=*Y}w$4-}x3P$d}ZEQC4ZNBOH_3Ia@ROWMoOj#q4F#?2yNxc9UtPpSzAt1L0q7Gdk zxA$%tN_)N9lOEh|A)ndlyE#Nj%-W``qC!}pkcw<&Wwkq3kFz)5sBQ0alf;_3APdF@5hZ0=%1c%b4t#)k5xET*G@zxW z1)5x1Qo@4B(g5Nn%aNrWc-7|X0>ry~V19I1KK+t9c%|p&+C-OMKFZ^qA=EEWW!$#n z-Yi34a__T`t5wv6tttj_-^k6qTmK zkpSR)wi*TNnyvMr7K@twjwlVlMpuD1z95h&+qABohwexBJ!184eT|>dK0{D=KDSu_ zvI4*yBKCq%II{fiQl z{`hhG;C-7Sgj{QNaz{Wl5NUQ$70oRDK_0r+kx3+ zHXvBjfXCu>bxhyxdzZG<<|zplabtV?4Uhv)0l_;yImz?)hhPJ^Q&sgM9H0P@s*sbD z%LVel+JP-OS-<2badBY-&Uhjbk8FjWcmmLQ>O79jxL>@;zB<`b4;N0RC%UN`T7y<1+K0SR-KtO;%gw3t3lb|#QMn@w-vjYEBVy38{5^y&d z$;il@ot$W`thRh`NKsNp8gH){MPu@05_^VvJ`r8rxNw+UotNu}MNzJ98vCw3B*zo~ z@JDL1zjn3QXzzI4nrG0Q>yu?UvO+n-O|Ckgd5hoD@+jl===^3!Ew^q(&(Sh=G02W^ z8W+5JKL-SaO$5LYu;ArWQ-e!K-)Aq~TpTch>`}hYBPUP(Ec|J1@8I z1_u)=W|SNn0xU+YtD{r(V5~v%;}af?n>-&7sZji9MSc(WyL_yP&XW3h=r{Brdq=dP z#>|Jgz3uvMikX{>%e0U6Rnkd7!PgiNVFG4P6pE=ickW;7^Wxz@>P+#2e?j}60?q^r zNV+2c5tWD+Z2K7fN}iagyHmv5IVMx#T!b-?^PT5NCqKUwC59OqR)YruhC7c;$yS1 z+*rpuCM(%H)Ol;z=9ks^k&DaDnr7`neV^y!6*u)$s*3c(Z?n3d6)sy%JkKc8MW~%} z5*NsHJ4cZgP$uxY;*q0>iHHW@{NRfaZ5d&B&LkQr`KcE)*1guO_DZ#!|&2ugj^g3lq{qiLq7*Un*i05-Hi9p3DdfOibtVmHU$Abmcv9_}g zv)B|_UiXS`&fefQxnI8A45O}Lo|R8x%K_k#gd${W4#v+2B| z&7=P0X^FVlSaBe)d|IkK{A8a@7||{w(sl|CUh4<; z(Jaxd0>v-V;Hl;VgK`4s&&7ahPJzfF2+uwOxzTWhqyb0|2oDwzFjf!Xq(+o$!F;GW zIOKK|j|vgbBW`*e#*!EZg6jy~2?4s_#4~_YWrcFB3In~7@`M)u@estc0@kqixT#YSZJDMu8Z@=>f(SnASEEFQ26U3^@sg8 z>uLPqbE{*MpCOnyJ(^j5dLZb<+duL)D;x$z;6GJRB8W}Mpq!chQ1@lJvUXhO)rLJ9 z2M0$RlqcT$&MS9HL$BGH{#UtvGCNaL6gEx;=u8Ns4(My!V20%YCj-&K5p3xGa-G1w z%edn?6l!D1;I{oUgh45NfH|pVfgP9yML^O<)W3y={Q}KgTWWG=*P#o@lz4i2{>YIG zwcZ$r)u=EKC;`!b&^G?|Y+QK8rjLX_EVk|{%k++OMYrYk`R69}A7~qO5`u2Dv9cmY z%%DI6LVtiqK4u3o;mN|@M7c!CrJNjUI}i-c!zN{!czxAN3@gjF%twkJG_P~sl;h^( zD-1n0=r}k!-jqi;cxMO8#x~Vw=L~bzL2!=(kVWGG3jFU-42OjQ&!$cyBqStW_k9D9 z70v}A-3c(=J39yt8gnM(D$xC@Aa*tPuTUm({zn2aJF@28q@KL=u;>D-)|Q?CM7uR+r%vv z0e{**T^R=ddwX@F(O{>Gjf<;ssWkNbkzR+70NAJ;+I(<;6M*g5S>$tddDl9%K>fMVC>j4tx;x=?#_5Zt+m)dpFTNF9ZqB>Cw4dqf89i8Hj&IE3 z7Y`RpS?x&-;B=;U`;Pm8K4uUOj>$w{Zw|xM3z){YcSr7>k8K$kd=pdM0|NuC!3`pA z2lV*-V)eorXe#yny}_2%cnyXCA<;Mx94(W00=M zrcwljWI;V&^}_eJsVPkH%kfccEKH3VAwz?_BoiBNf|V3qF6gi^AhHWUP3^ag7%>2;H2p?lx4*xH{*w7(CVu2ctJzGk^zBCuMDgqQQ> zDWAx}a2~T||1Q_xs!&*JAJ#I=br9T1fkA1#y}jSfRM8@|j0T70`}XU-$9jWqXxYfw zYDvT2&>_huBKFzgfZ%^jJU z27b^-2XER6$6E?>9Y{_%w*a$F=Cq=xc@9g(|1(6|(D3Mp&_@CHfKK}nrKJT(4WYh&<{OR!Q1<0T$!0M`aC+yG9(0Up;^(ysrzl-&`_6)c5QZcVc~Kn zjM#Ir-Ft{vf9H4g&!^OtyXK2c-MZ}XnU)c)w%p`FYL!97H*aV`ItI~_-2P34fx}5f zV3wI%S;>I;!pgrl-G8|M%-!ebPaUpNG$3{Z!jNa|dnJ_`xq4=6fz1kSao)rRvLp*n z2{<`(A^S53V_}vN3WC;P=h4eFcitawh9BgCnFS#C&)wPDlL|jsOc1vl+GSnF1Rn@w z<~%?gVuALmm??baeKOBp7q=64CWO)zT(-IOfBJ*pfHcp4y57g@gvV?;VqkmHgI!DU z%I=9F5LDJDt>_sZ{VMIpzVS{PjEKPfPhSQ1BCu66K19tDKWIIXEq>Zh?N-t1wkl-3yZ}7@q>d5G`Wn*K-FhrXOBd1KCr+3 zjD1>oNP<(n_K^jT%2jE}#xtwYA#wU{k58*CBSYXNIT5Spv5ikCt|^)bTSU16b$|jx z9u64?qUja`XT|mO?~k@Za*idEdSTy-<1~=6iKZ4UKsaN7P8BI-2qG92h{(_YF;~$0 zoJ~qv+P4{};A_2b@d%VIv7-+Q%ec>RWot_Yq{X1%2qdsZh!7)StFgjE&6@rK%l{q%&e+G) z?_2y|;7iAkk!V2S1;OEAJWBo;AbcZCJzhW*fP;Nr$gCR`HT3;?vgbY>6j}ttt5Yxs zS?F5x@@D^KxX_B0LV|x_@>7-9(_xByoSX5;(1%!kqv&DmCWn@C9 zrlyL($zxs+i?`=hPN!*c2}NKQa17J+{?^le;e4MSs7C+;sJ#+)->*`pCK3T9gm6It z5Tf1^q={!nxya2+Z2nMbWbwJEwwC*a<+d6XIrCuxH=?><+Ww9xD=)VJ9^4T6!%a45 zwyiPKK=V94Jj`a*u9yIojbI2<%`R3oH8oIhB&w#qJ|8fnA3${323hQ^%&2{J>*f2r zrBNWWsDlYGgBAc^e&Dep5y$|TE+?~%m zXcQ_v%M|`Lo4SqxM^rMe={MBZ{k2vg(8m1l5FF{@=H0pQKB**OFqj{*kdI1O-?Hw}b>4FSQM-&$f% z0glEmh#mTb=*|H)#PYFYIvg&m6yBnG%Rw<34U!)}0hhN1=VW;8W?s}ZG%$dD2&_Fu z5TxCi%&QFJ#}*V$|7d)Id5;L%>+0ztHoi-XU~li?1|B!Li05%Ia5oSx+s8Mpt~p?{ zBETFXG4^-kP2a>7kTC}r7Sh3~Bq=T~u7u1?C2#_V*Q5(lqTk42#Q{eD{av#Ca*NdF zP-32ksA%r5(lEF+kR1^Ikd`^0!`I^iprKWk&2AW%FHN66QfCRvT6HF4c7kr?zpg!Y@#`_5zPJ=*( z^!1xvSw8J#R&UFs5_xuZbA06uT3BLEp}h1s?`CXl#|g29>k62Q{(XV)k9fSRF3PM@A{?D7fk^C@t+dk0m1e@dH(z;*iM3sB);>>&bOg1 zT_28!Tb4bUnVHF;wr)6en!@_?VM~f522=(qVM?@J4-YIq;p5w~gq*8XHOBhD+oMO1 z;_njD>;L(A$$(t#-m@niU57oUIzn6yM;IYfOZA+q)$NE~p#z z3(Qe8CV4CPkjD0=PN*nzB`WSpROM-D*ETWPiT^qj*Ka1k-BSIl(rnJ%5rE_S0i!dn zuLs8O6n$4%gSM`(uaCr-g{9&L>oa32UNTbCa(`vQ z88ECRd|wy?+f7ZrBY5S|C4+WJ7a~l3$BqR76DHi@a7}`Jvv+hP%+ioy9f2j+3dlH) zC-{PDAmfl&3z-#&VpdhPm@1>#+MuenRZBPsontDG*V$iAj`gI46IoYpwq%tlPmw%9v5b|?y1)tGN!mGA4 zoN8>uoR&pcR4FJ(8Xl1D*`xN{DAn%|999tQ^`js1o*&tP2S){?08Z;(COho8N8p$;b;TuMS`^Pgl89)epH)`HaPckJYgUdP$Ex~9HtWP3J*JRjDI zb0k?|&u;Nj7m%>Aw7RZvo)+bR3XJaw_y{4shbN95(|~3V!vN#%-PHf{C-h?koSBVKlj-6sV6Dq@P3`m6YrsCe-ihtv7dMC8dIQ6ggkv7tElCX(Z~ZkmhQ6 z`ueeiZB9;e%5Z)4p<7?Z-*Lx(;(CgYr-N>qx?SFK<^mv(wS%Ltd`#~O8z)f$Dd?zg5)kPr8!4} z$Qbzf8XwCtfPEy48AU}!9)tOIY72iVrpx)|<@XQ*XX3M##4s7W7k&cOCVDvg`T1s3 zlPsF@(YiWs4pvGY9-fY0zp6Af$W2Rkn-4N|3jo}St4F2QS;_EH>d+&RUzBA!lsY1+}guT^10h3vBme^cVdqSNt5 z(nrrS@L4UpS@kHTS!zQSM-F~l6)P4M~Q}AUF ztQR4S1}IipTN@^;n-y4IuI^82_3g<~M7LC-0lTRc8F(@|Z~4!uu8%i%$sV*&568Q6 z*p2sr)34iT;nHxGyIr|Llp=agxxn@xFM}S4>(nHOx275Li@q<|ruX3XrA%#T>-P|L z0{oirE%(*4S}j~AdG!vrJ~OV*PtZYmr{*42C-b_Mt$!&C%vSD=kFGkvUIjJVbQy1t zcE)$@(lilCt=Y6T_2qnF5V;h`qK@jJLd=LSu1L(pPk-;(qURiFJ$Gt0BkBz&YchX zi#(t)reYV|IiS{$9s(2ddZUH;KizOuG~Fmt*IY9 z3OiuDjd*iWyJwFzJn$`xqTBcKSs$(t*!!jb{LH};%4Z6S(l+3Uxz!HnUE($7(vmq zA?rv;fp3mp4{aC>C#Y`>xCc=YG64FHwZG^J>rM4=BHKqqoO00mkiD!>n5j&Vbw1rhAaG z5g17;;nX3>#f?`MssX}_4Yc+)?p=?}NLtwaty~#M2!DRL{{|H2^g-{W)11x*PgG50 z9ZTkPOj7$$i&)y2yT%ke_sR0O~tb$yWL{KnCLzkM7 zrMma+^y4l&8SzWqP6C4gBy=P*WPv*Z#YjBa+~~r!M6C@U+!}hdJDw+W%zpS^`L2V$ z5;rx@;{pa5s9g8$TkPM-$t%srTXK>Q2};Zk(GepRq%G%}nANWi)1Y9-K;B060h*xv z!J6mzL6L%g@O$5Ln)|D3w+G0hDez$l|o z)-Y(>znLCoVYDp!wX}ujx~}{9;2*xwl$PT)9naJmOyKVzlL;Li2xT?&Jn00ddeHpC zy#Vg6HWT`g@t^uCMd~qN;&lb$VoZehhiejbmnR-`0|O1~N2SfJT1U^)(s=o3S*N~m zQ<_SGQ3a3<#oxaRo~%BruSHi*>gf@U`N@*eiDPZA_UsqjPyD`enRUh2hWs~`po@tK z3L@99;!Az#4tq$2nunLC0O(J`(9#;sarJs7*=5|4b?fV4f7;TQKggkf8w|fN^fN_q ziM*32gM>n_Dwwvs`!=1dt8ljN77fc6Uk#6+I`hBMI&Nk!SZwQaIDM8z?e+g3Vjj5r z+rIva#X-xaPTD6%oAv=iwiE~sk-pfto1m^S9**{W_)t8z!pFDsrIM0TEo5;y=0#N^ zcR2S(+!4TI04~tCP5k8%WlSZw$D@tjdic~23=lJ68+s9jZ`6}LA0Oj2sr}sf(0<)J zqtMl1rmvE~!L=Ucp?-e{9xdN}W@{63^*gWW^Dj=nX8+tE*b^3M+{C>Yyz?uoj)ZN95`(IfaJBLt5MfF;3~ zdvK30#ps6N7j3d<1JK{6)Hb(dDx8JB@_U7fUszId9rStl701S(9M!D3W2oyU5zV3a zfKcoasqG(}j$mFG?%YvMxw)z2VAOgtdqdv0IF>eTU*O}F$v7PphPtDMvJt8e8UbAHxSz^n z`1Mg~Q_X0cC=)zn#9>EtP9_EmBa?QM6dv9xtDW{PwL>rM&<^O0bbKX40g!kk4S?^`vEHygZrCt*cYkgttgiHf)uDQkV_C;iyO|=R zqq9)Y0Jr^s++eaA)p2E`S%3op5PLzsG(9YF+uPdnaj3ul2*|#jQ+;aR_G&yM_%H~* zV|~0^M%)Jp&G*H-do;WXU9$<$6B0rDfwUw0>T1)RC)y3c5gd(8Imxg%osC|_=$sx+ zYhUHQ9Sqqxi{A`aq8@5?ZCCNw_j9cXaTnlgQmQ930hDr@vS9wlE5mkR+|WwFYO_mu5_Bpsf>#dXE>OXMCfHodu}BidGJ z-(q|7_Ue~#?5Z(77C=hmp5*7{rMVA&8CE_#lPIa#d6DcmG%Lc7a5U( zS+{lMzDT(3A%e!=FwCZ}`}^0vzbd$`^t>{S??QxPThKNpCc{5ZPqdXA^0Rm3@*H(+ zOY7y3eH!%rDGvNmz^UoD@#orFpFbzC5*Xv_TTvb}4m-|$EF2oLC1jCeW(EHhZ{j5| zGc$Vu+ZnjUID{WPfA);v=MUSr1Zv;mJPMJg(McWi-6_h0u3x@Bu5V2n&__#TOhQ&s zR`huI{|S(-XI-U@jeYV-`A%;Ro&HA86-Gwp&5Txx)^8OD-k`V4uMU#ak;aGBdKb^v zvi2Hh?`%KCmw~IEKzabA3x#rtTMC?e8-np0AK(r!hOhpG0h>ob$~dtC7qss@6wLg0 zVPwY|ea$D|XJZy6+C!r{haN@|rUj1nT*25nAr+MrD1bWkSUHc+mX+B$IW-3GD8KNQ z!Q^g&G>~&>0YaDVr!+0Ga5vFFR-JE$re@Y7`=nE)&(a;S22Cfi=7##@$)c#;lD7a@ z=z#fe^s=O!A7Lsgntj`oKRkG6#S_pkTY%5k$WNO7xmiYwbtm;l>3@8QE>(W4S%4{|~({cg}qH@cF;{erpd{<)6nRgR6j3 z#ar#&H4)e8?#J2bQIR%ac}J&(CT>9HV5zxlx7?n>S2F=tfd)7K~n_) zJtjFh*{JgCh<=G{uL@9kZmMEiTgH`T1N2O1$5E?0-#fOwvC$#u@80Qqdo?wQ^a287 zd~e@+YjbnWb-1nErsi8o}?kBBgU59H@7BxHY%^w-9bZf>Xm5>fYtcI zwVE}GZOj|<4B~frn!m@^i~TY>(5$;KnN4vJ!d{j5E?HO8rl6pe2vAscvRxf-`ZV*1Cho#8FOWZ*X&05j8w7bBs*~+2t_)^70Q6MIX?M^ zVoitk)F{^~=6`MoSo(h7oTQvQaa&))`eH-Y|4+ER-tbV$yuoL_<1#yA+SMI$9=F9_ zYH9|oS3W4*CGIVJoDUT_eCP=G>C>-x!AyO3W_Hg+DJ@m<6cq{d_t#vSnpN=8n5Tr8 zG96KBF(6vsJA4lzsN@Ur3u{WgwWm*HF#KV``2_-#3H0&6{V&Vj8L0uVFsD~kr~wmy z0ueL|mxAWdfI9L%_MsoFck>jWp?`|(7bB!U9N<5;5=l+}6#9UBphS}8uY?eADCnAf z5aMvF`_0`Jl7C1*zvb5^^@c{Y!IPsX9M65LKr-&&U~6MjUFmxMyr$o(5PJ%=ibR*P zfdm{&0VBnOcV3qq+G7?ypLbhu|Nd5lFw`|Rif&cB>k7LwGWsZ1UD};5SMjRu_y9U= z)r1+hvNL>lJQLOq*Y`p}v-oD8^%9+ro7IC}>6I0x*A@*Ee<$zrQ@_2`kR2*^{h{a! z8JWYU#(5@&G&wVY|kZnZ*Kv{1B+SHDEJ{rD2%(wua3+JECd z(W_NH5JT+GY~AVt|CXXxZB^B+oa0Hlf(+bG7Eiq}NPH``*4X>;7$o|(moLBgSFV5Z zDVUQM#5>NVJ0Wr|lfyC*-vA-nA;A5UNpiCehgMXYn$|ZGbi?W&*ikNjGMw|C6`&UwLsogODomiqRQlxL==*H3y)!JXgKP~theX3;ih zY7m3KsIu242itS(hNVkVN5{YGMXryvuA7@%gbV18T%+A3oEr2-GN3TWVzA9rQxisf zxA(`tL?oNkPoH6W5`4woJsr=3upcHX-zc!snye&8cB02CkRNSX$982kuYS+gMDK~V zjF97Ihf*EQJ)iaK{pHs8J2)U)KJm}7;la@Gu$CPxaFy>J_tEai(xGJ|jE3XmP08`h z9LGp??#4zJ$y=De4HvunvJ8Kb{!|olI3h z?*C}ho0kUXew}#}+fyEpN?f*xo#Wk7e|o9m$W`er72jvaH`?#rgbehPjXmX{x=k9?vDRYgdN4P?0z(V$R!1 z`jUQ@6(N1#f09-5dM6_z=42(GSX_2mfWnaKA)8-VxTWUa$Xav_Qk<-;X=PL_LMy=Y zl%HKDYOGz1kB^Wa5P1mmFMG=M_r@ZQb!}di$zE>hLF$U6ZDBv>)e@H zU}5gKrXi3C5IPQa5Fc)3Jq&Qp9eFpe{~J33aPadDAA*=)DgGeDz9-m>f~obL1Ln2i zXHzgG5>Y#tI31z9wC|~uPpz+ub4qRwLXP7(gwntm{T;sVvaW08D(Gleb`_N`BrEEP z{L8U$%A#jDbLW-P^(TKx%s*62`yP+_R;Vnb+j-=~aZ#7*X3bBhYk@Y!q1;1lsg+VA z=eoWgKP5i>=K8{Af_(O(uC8WXf^gsTk&LHk>qDul%QF}L7?i7I$vWviU0sbgf6n|q z$0*o^j!W+ibgP8m1vqz0`~l4(wv^-1Dc`>xR zzBM!?pi)N}S=+ECSKU~eGK5xINJeI{I-)7*=(|h${XFMX@FKS4+coU$eczn_XGG_^ zLQGffy^s9Y|2>;i?krZP_{lqcp#zW#eCj#cZqbrHLdOY%T6b12?++EIpVdy zzC4re5u{`w`trg(s)wc57b0F~Ta*~-nd}7i+7h(36#nN&=Cf>lpRu{#05mmEKwZWl zF;A=9uVK_@J(1&P{LTVbRL%95%a(q^&-j zq@)}^uJvFo$oi$fuM#3QeT0v0Zf>?^>+=u@5fU>e2Nm{e9m536z@HmD3=K$i z;Qk}iULraIHX`X=f^r;q82h0=$fCT`rE7;o@TPU_zmB;{lz$eo!>ty9{&C_yS6t#@R*;U zUg4Ycki(9#)ffS&7-Ca%#pO8_KC$&P@sIV?-?z1~)r@&<7I(i)<+I64?H>*C>Hh#{ z8k|WmLWqp)#7|>shVC-P>92MH%FVXT*N>QIAAMBIE}UL#b_ z>cD=j@qmZ4ghL$am)3vd-|wsBXv*Jh$S;Q7=^@~kFQ*eB_xTc?85TxO%!aAE#9{Mz zY84Pbpj=_|qv1+QV<4}ef7gAb+=V@7H3GrBNuQlE8z0qynGwBk=D+zkvkMmHG|F_yJw1WZmoTIf;ssoM zBg(^)OlDUd+k_-CHUW@HUmY_Lh7*jKEg*W*3ZHZ0{2T*(>u>`qene2SYMkL#R@QR=|O?CKRj35nNunRE(QXZL!|17g^9Mk@aia(I_-g< zEllQ>=@G7gm0-dpGI9D<8X-B4eI*w|AS_6qsGdVumL}B9fBxu!0sAMDP&}%ZLb+e; z#KZ2sPN~IofdOw81_pYy`(B=^MpMnniaG)U0?<+4xLhV(Hze8F5QLVXA zoO*&Y)#R&ta~|viN`G)skq^_4;}n4gZLC6pxd-1mI5_-Y;0{4)BLWTOuo97|Jv}`$ z3k$(GB{-nLu+Uewx3fW7p^4+Wmpa(n*TJEpzq%qdfK1kvY`C6Obt8Tv2XtiOVMjLb z0kcSGO2d>W@8vy5A?Hcg@oqbb8?qsycvQKG$Ti5-)hM1kwO)_f14@3M%xm-Q=4Dat z^z?9Z*8Dp6o-rrm)@2rfzOSvB5{Q@cy4a|OY+?gskr4^LAg~Ayu6%1V60$+cwQ~g& zzP>F#zV1+!b=+36w)7@VkM+~BtPXp-;IJ2`46jv?cz?W3S>5YTo?}^Kqp!aX)C-;~ zgiRh`;-j< zTv#;|H;MD(6PeMPw~RG+vvNgGhJ-AODRKUMB&+k0iebjHj@dgTG;^&o$a{YEu*uVt zfVnw@%-GBdoGNkY_Kl6z_Qo~uJP3x;SG&ZirRb+?(#??%cN&a(j>s3!r0c2(S#OP%~ z0|E6E0(>K4A3h9+Y6~|Dud*^cN#nXFPG}*51-iR?@$t9+7Zt^eR5nOF5GQ;DWFAV+ z1zTGRsNE1w#s(dd*8e!)V7#!E;E<8@lNunUqJr%Y2HEt;+sw1wp7!*?=oDfemjPW7 zoEeHgCGe{H=-Jtd`yy)P{t>5u>=8-^@X?{T1IRQ&NT61EYkO4MwTF{iKedWiU|Xet zz%vZNN;u0gS^Lqc+(JT>4aD!d{}!}{FMjFZ|I?OxH!X-N{d|dx_QCUX=)gqJePFD;TmjIDR#sNl{Da$(&Sxxd@*(iG za#%Tg_Y4w}016Xnn`ozRqBS~#0Gy5FK=;8wy9N?JKdtULemo-v9+gj!qxJv&8wG0~ z5zLQWE`(_qL${T^JvESOg7QTg5<5FP!MQ^4h=@KKR#sMG`2;X|=&bhhQmf}ss?iUx z_-Pdr!pU>zXSoXB0P^|ZyK;!0TPX}pUK&EPcEjyYAij8w_%q4Euxsjeb& zKuO47C*_dMpZi%L=J~hdCJ>SV&xQxOYb%#bI3-9J;Dvo-Q)7^2t)*0cW7jnS{V%hJ z9VPs0=#|cWWDl*bZY^}RB(`9<@8GKf7jB@ZMM=9${7ct%776<PKpGL~B{Y(0&oTFQIh7+4yhw!ONE# zs4O^{5g7)T||iK9@Q;6PQ{P|R>KY-6L$Gt zU(%j0-wNsNDl*75*(U`o&3m_f=_yi?+lg+?|now`Xm! ze{^wDvOtGnd37xl!9>LEu=wgMT+VwluQ{S!+lez!fhOa=;KwUe_U&{Z{5Ys-$U+!K zN1x;}?9A8;+UPTkoJ2A?E{1g6%|+87NU&80Bx8CGAv%Jx_jX7KpMnBA$}61CniyD# z{A@%Kqjlg+mh{mPE9^f;Qrk22v3DCYlajJBa>)i1jXm}=Z+Gb}i$)7L2q!VY%ShPP zGNOHklk zfV!K%#>YFZ57U^ZC|9=)H_STzF6j7hL-z0F-xl{FU2iVcJK6#3H#Z~MEk#DR}=ma8y!pXosQ2^#XC6y+P%xZiOMzZQMRnvmja6(e71M4mus2xA~| z`#Sx4x1E}rKhVPLW% zsfC5-6*>PTF}sxy3Y?0bbf`98FwaOz7v2fJ%V2ibvV2?HkF!S(Y=+IPK17BDs0NM? zP^f_ozeVV$(EW@d|4|rAmzG(UCTSZVQ?)&d89hk&`SI`ckBw1-%jcDU*gJtmN6hk zqE=N=;m-E@Tam*enCHCgCjF;_>1%J<=lm9fq{9+wqfhqz9q(jOJ!_iwBbG-6SP{PP}^GLIwsCi`KwdAtg-@DuTIszd^uX;h~!4%=7VUG={Rbk%^wn4lF&b2fL0ykD1mr<110SrPQxf7<(F%YLYJ1Z%s{{6&?$h@B9F?^9jmX_lb5f?wnHp z!9vaF5aS=S!;5Z2b!hGkXP8eQp{3`D4n(bN58nF`M~- ztqX`C<2mIDC&se(jlpX_b2=QH=t2ec-GqgxYp_&6kVo~)%E@HnMXhaRLbLR_D_!Gk z-DEchi5s+nxb;|oYkxuMrky@GFx{-@7bnnFnW(TpKH;LSPPUDKI|>mIN%B{=VK694 zo#>P*_+(L3JA$k0k+j1tRIzkI zkP4d6{lW62%b3qsQ7+aa;u zStuT>7as@`spsA0;MbQgKjUEx6cH9qWZxau%ETn+V0|lOx+i_&gZoM`c|C-p6~oTI z^Bg)Pye&(qpNDTR7TA9Qd;f9{J4Y1jrw%j~Uz51zp&nZ3+FIh}8U<&8&Y2+^A2wJa zOh35o<>ca$1URQvvLU%3#Q1!t+zqy|ACVik=C{+|kAo-`@01n%En(5oT3eeNM8-Wp zobnILgeT9B#EOim$acemq2e@O{r*>{kIXtt1n|<($TZNxz0%Ry$(pchev9%=>-z2? zz2@!`r#B8QB+m>;y;qoZj}O)MlZmEvZpW;I6YmaFDG{%N{Cf)5a z@`vuPJG2D^u*rb1pkb}Yg5g2}f@2GuIf$n8j`OdI7an${r5=>Uuo#qP?LTk+Au#{c ziwY5eZ78BaSe0}2+Evr-kO)IOj^P+bh6=XVK(zA#=VzX4k~wwL98aRYeVfFtVZ^`) z{Uk7kP&7LI5JnMBq^VKJXtrD<*_`^t&d{Srlsi!R`6-{z%;(`#+=WbBT+v81x*ZzI zFD<*vSB&Q611Blu3afHnrbel12QtZ`ZS%crg-?p8Z>U-5LI%Nr3R z6-Fk)Xh8HpxYguee~Dyse?ZE?L2UT|H(&kH`L{t;6TOJy)j%vr@?4w<6|N4cyzoo< zU&jX~W4w%6zYEjRFWk$^V|e>|s2U$FJ^N+kSc14D03G*>HT0%!uOk=j>ie=HX!P{< znaX$_Oj%w@$DN7X_1ge#-G->Cip#umRa8`n63DIOjirF88c4(kl!3VO+sG;SzXdiB zK0!8`2~`l!^k~w9FTEXWOPyNd{O(>}HGnM$fiMLHMVeM5+cO}4$e;+q7heUUntTJ{ zFQ+{?NaRR)_c`9*6IiA>l+=9Tb7tjwC3YQzBE65s69M37afaf z|68mKPNPYqF`ZEkd&Wy;fW_1G*!yXPVw&po`M!;nRUqtP0sZH~hcw=}NE_~ccwZzH z0U#0%4(=}=DO@vU@|_QP9Ue;W?*FU!TIc1Nzk8^hlzpmHE>O{vh?(Ux*c&nHaO`5n zMEA#1CS~wb0NW|2EcJ zu3EabRP!6l%1U(@D5ic<{`erbTws3_K`5~ML{H00#fkNGlAYMt*mwhPBB!c9C%`-+ z)}L4x6cG)rJlw)o7cP*Bh`d}k-c!fQeB;Jxz23+Ugq0mZBhdo$DN@Fzx8z73Msb*t zlY45MExa=kFhrdl2s64f8=KSc_wQ@JPWM`iB*)$~g`&A{`c8 z#6KJCEE?F`9Afn34LfW_AFo$@yy0|eJ#33ukAu(=dSFVhKS_+PAr9|u6G}TDbm$f> zNO^CRk|J+w%lp@@=v`c5!xI=$=+dusU)X+QfO2tRJjh)?CGB+#lS|gq3|wESz9R5L z5mh^>Iq696O>)EG&$NHTqd%)YcwkDM=Kfe__ud#@`xz^!(Ff4I>+x*7u_=l>r^{~w zLIn5m8^DIR>=^K_e@0iGTrPU)V*W&x>XzG`2_kcZ53%8pc7;VnX<}&#jx@1t3Uny~ ziG&rbI^U*4LCY93AbUe5zU%6bzY|#5vgM$c!%OWPR(pi;wt0z>=2e;08$0iJH0a{-09=TFJ-{rRIN;}~ip{2AkgJ8p7fx)eCuZs6@LB>Dl>h(Ox_xZG-T{!LcC z5zET4?XWiWU`-FZ8RnC#`~5&H4DvaUO`?3XcmL{ zp`%W{*kkorMyuq`e|^rs>0X}N+KR5vk^M;A#uJ8;E|0XMQ2ghmMdXDk`OFbAYT)Ar zPevjlBQ=OTaWXQiix)!(*cnQ5O~7O5^HeeS3NP>rU~NZ~H60M*F>*{M#`H@6PC9 zm4FmK=Ac_0*H!(MurD&(_^|@3oaZa=hKajvUwX=$d&53Q-fU?hdz)G7?%|=1EFlT+ z87tAV`I_T45SAUJ?=+_)raCBPDmpv&WK~r&ad6W=zQZ-Uyy<6Ne4TY^x>pC~muYe& zN?Nu)#6ESvX6(PYGP=%iudc(= zf6n$~^P*qOL{NQhR*azb#G_+%s*Rwe4w3|zl$K)20C%0zHQ}W;(u%Gmjz&eV%r@!=K7Z|iLKH->&N=-jZuN z=S*i>BrD>mYrh+G7ghh}kln_u)OStfblxpmH=9#W=?q@7iF_at5jb1F_j3EXHH`~i zH_xfAaA1X2NWMlg*7o)F^*mzUW9>p@$+asslr=39p|>6b>cj350{O*uz7(U%6UT?> z=no{xyF))~k8BGCuc^kepr1Jxrt9K}f5!iE59d1Q&-RdZ|ku21zf3*XYQV zF1KVeh%<1}n7sIOP-5=(GSaYE)6+FpJ{~jxJ|9w?wRQ+^NWOlegFypp=;}MX<5O;g8J6B*(dX+ell>t%$HK3 zmef(YXFe(E-q-+{l#Y(cyMlx9w`fTqGp0UhNLIJ-=ZNPN=k;^vGMXQVdEZzNW^YW| zZpOuBuH;E)-^pzJ<<4g`O@x(E<-+aa^ zStWatolV@aa&ST%v@DClY;FS|C^i&XR330XBvnM}e$_Uz`v!V~G$=&!-n^*yhXj$BNuI}eo^qrUbZ}D6TRn3?U$oGbQBb}ywE9A~Sdwe3<500ojm#|xSP;CgUmpM{D?i_^m*w<>M~1F8_AQ#v zclNDhh)u3}63%jSWu6J|_~pZ{^vH6QOHm9}UR+FUFA!_*Dh^J4aXJ{vA5d6v8Ytnazk6E&+vpUyhC&l>|Kru3zq#jw>sJ?c zFa5Oa_~6XJ$hAY%@(#{mMQtrL&ciVI@4`!zKU-rgj?o;24%!`ocgX4?2Lk>1gZm(M zvb^X@x7J$BPLbJB``P+8a`HbE6%i_hXBH(gz}iBv`9S@r&XA(@=*h)OabxWR0q6yf z(oEs2FhHCH+~tGlV|Y$B-ixV>tFF~?{Lj{nIH^QAPf=I1dT%Bu6?V?;YRP(44dhh43$-m1W71$B7k)t+H3~i+E1p@~hIt zUura%S%OL!xJj^5lr2lIvZ;wM&1=IWk*angTCtr=k!;TDhop_S5v#bk+LVaz3~xyo zX_ve1$)@)zP@xbTsSJ?hgscfhO!F(f71+Ymx4Fc#HC<+`B2ZGLMDm>EB$9$g%+D%?Qu4~msy%N;Vy z;HI}A;?38|>+)GR8{5&Xa+moAjEK*)pl?tl!h>Nt62DaCg>8br=;(Hyt|NA!=;(x6 zP{6B~m#?lDCjc$F`?qg-@Ya-i zTnqj^@srqwGvzrIm00F?!6KU;&aaApi_TUZ61RgxiBudp@0m0YzOdUAr?7A|Nh7CnCW&T_#mzu5Timu3+r6U&6`{n|dqP3;z?n4{vyr6gS z&a8I{-ozW23Ej50-xTMcja3o<@URLRRn_(H?+zr1oc@%4W68~D;c~^7Z6+^vf*~aZ zr3wHVZRZQfRX*8aQS;zjG+nW{tHY!!!6^^I?*V~?NcpED$iB|+n_SWe8@d2+ot z8<)Ex|MrKUg-UnqRf2JY21ruwB7x5C-o1M>pPK~v=H4(e z;g19gMR0Wd+EIb%`eBikaP3jTY;A?xo|=IH``!~IfjQkZb{j6Ug{Zd7!gp9s7Hh-m1>Xje%2+1OO-Q?FWU38IBjfi^ekzParQkl6-Zk1}8}} z&r)^l;Q6R)jz0;!*>c~_l1Go&kqUCFnnNb-WuYvb%sa|e{L);VqHv#B*W4uk%ds)E z{rHLFZrJ)qx91SW-o0l#d9GcRBfq6tRbE{Dx1P7TLH@e*+Tumqyt~3oW}z%Bs?Q=* zQn-m*0@W@ZClTgmV;7S%qn8Q-0#ev`i*ozW;`F0Agw{LLvfQ5#NHAY1O3{@R6bK30 zAb}IT@7LjiZ}pdI*H%W=4}5@tj6ipwPU;^SSxmom>{ELa`{57O50HFDh0tzNgnuQD zjbZI`+F*?&99XZC<(L%M*p@ACT;gK;@k?vAmoXxuyE{6=+MGkt;7v)%+OOBW?{m7l zyTxRzMXW`H!i1>5B}5iA$UxH^q_8eQ2R0rG(|mGC%GKsuk;0qQVbA=-dB(d zVIKCtay9*c@GsH9z5dhO$CHJBudkP_n(kvO(^db=aqsGxhvkx3O=60>M7AMhKLNfF zyU2#?Q|A_2`l}sE-E6 zdGzA%ZF?EAV|zqc8kPtpDSFXEFj=LfT<2(AWhWv+(_3~T=eaQvk@@Q5m{(O#PfDn) z%Z%HQ?l#pM+68ob65lyP(YDuhBldlY!_V2S7UC8OWIi3B(bPtGoE zBkM;~L|}4+P78Lbh+X0*v0$-$BWYP?J^RbVnaXO7{Dqadqs9Odt*!H}y);z|z5R(+ zuACC{9hS`oC93fdZ zHr56KzL39qC6}44qvGYoy;E3IL|k0^>9Lxv*!|U2RRbl9mXUyT=^In^rw>{Pd7z=f zOo!zy<9$}#r(b*uHp(*fUOG5B+IVDmfb6x|Abe90o!6zR2CVl4x*pLDv!@p~6ZraQ znyVVM2|2rFQ1*$QcxOqb6{9s*7dNNqD6OjL`IoBi@mF3WGX9%N5gGqZJvRxynGun} zx%H{x>X9jZ+-uoIL^ReJ@7_RU)JV-~c(IDH=I75{SowhCqls^=?$Kcph?n^;H6>wM4}~JzC{2wt z-?nbf#ypl=`PKpBeQWIOtA?tc@6UZjnCV$aw}GDOcRYicBqleaX%pscLxh77k!x72 zz^i1QpBlMK$`2RfhjONRU0MjibuK!V#D+iR=z_ zWhMLGLh)uEWpz7SrjzZ0zh14OhULt|9c4>? zK${364%Id)OUM3=^J5{NoubQH9=hlUv3D#Qte^*e0Dc=E(21bW1d#1(`}NU%bgyv- zEfk228!>#82aKMCwFVuQmX2zT+(`OR$~Nulh^{P_J9oWNUWQlg(DKSWZQ-v^6!~^R z0wRHw)^+;d zG(0tT@Ub9Exr7x5hv>1>R!SnxI~<@}$%9|lCXBEN60We@8($w&jE(nW@^+|;U6*p@ zMjp>f`iU648tIsbaZRimf)?=%tVu*+uD&LBT zKAcLL7!V%UPM$e5%ezIyEm=FIscB_kr4ts@(p2hgp9C8@W zmifZDvZ=cM@fxF0MC=n{{l*R9Ux8{Iyx2(psq0;(1`^P0Tz>WO;WZ>Iv?Zh} zC!ft)Ub1Yp2N0Q*torbyAC|-t3rDmfcB0VjFIx}{+lP!Z*FQ5-)R~I z_Z`gDQvVIOF7-)`$gYDzZ8>74m$YIe-To#fzR=Xvjl1yjG~bc?hq6reUu=qh)$>@E zIqCSLeLEg%MMjY7j5Hn!$+1khn5h-HPtojY)Qsfh-+laBamVFwruP8i#oddGXm0RP z4tuodqs8f6&#M!^vIES3j)Ma8lGawghWiF~jKKgmi^w?okkBhHeXeMKeOTSett?b* z&vbWgC_N&Rzx~UXL?If#wS_hJ#ZMI&G7nst+s9Po(RyD0?&6(0N%lQShqCqOJ=#Tc zZum$#^{PxEcSF`YZzl4}sJ7uw?00Wmo*je*kQT~p!e|Wr+1ziouP0CDyZNjw28F0X z>fwHD+?<^5sk>HiggjkpxOQ*(;B#Z< zu#%EUV~giT8F9({q0*w~%MZ%*xwu%fn>-0K5Skl~t10KN{g^4F=W}1rE?Y>ibIf^( zRZDW?zp}~j1zRA|zqT+84_*e&e;3wckZ)`K)2}t?Vi%vRC`O;Pw5)sL(tB`Ems95_ z?qz;(3f8R|TG-NIJ-;%{ej2Rz}t_v_6XWYrVA=~wl7FJ34 z6~1NbH<#|!;Ptzn1A`5qn>U4pf)c!ef8{;a?}fZwB6hImn3st@s<<*b2afr4f=KAM z=Pd6U5J|tHzg1+W1&n0ZkBkD26FJqKv*Eoi&UI}wU-?ObbI7d5ri?2mO z&6p3;`O)R~S%YNPpnUptj52v$tvt4bgAdAQANB{fJ*&iih&FXuqj=SLnIA zK4KNBu(P!6N>TaDYoiT^%|DIu)N`*d;zB3X!%MWgyzh(d8Jr)w(9@9c@32&+;azT% zsrG1$PK5IHtfu7$URtlARcue&+mxc*0160v{9eSd6ZIRt?mg6REX+EBH}=5!QaR## zm3yyaIJ8>eVj+50mK{k`V6{_GlqKhsscK!9zj;jh{L^j9^-a_*@4{1^A;^v-1wpO)Xf9(Iv9Njm0RO-Opvy+rp zRy>NVJSPO6vMl`ZYd5WGZ-0Pc9_xm-;T~pjOy0ALie9(wgd;O^G?*hBXXme2R1y4^ zgH7NCX+%}E&fRB7H^e7;Vph)m?P(T*`Al+6!uP@itACTssr#)?Qmwy4+vUTsrVdTL zIQ4{qCZprYQ8bd)LzB_baV!|jKBJ$rO;4}VQSvcb_SV%#df@`Wlwv3W#iotNSO;8y zZU?X5>c64p)KthN{1t58*@nCA?QLkp%=j4jDS6my`$#bUy|S-qes~a7vuUVk~s%E&*_V`c^6U($~UyZ zH6_)t?>LsKY~z%wvVLWnJMP9mtADZkJFwua`t7CXs#k>|f3iKZGyVK)i^qKnVn|AP z+DwSyp<_V~WHo9Vkv9RNFcZ_yzscXbuDq@@G`G*qK9i$sx4$W41oKZD0!|?MkVJ39 zqy2VOYk7pXkXvHA5NQ(W3$kujt)hS72^)*#5(oB5NZiByyUys`b+;0RaZ>Va_x6(R zFy!*5z8W|3v!e5CzFu8y`02s<(R{HJkHjx~%*@UbiGsuqPkhPWp@wBRcJIiM497h}&I_3FazvEU{ zu@H>&C_0O*qm%z1OWz&OWxxL~O|u9gQbI_{C`71iQAT#irm_p!H0-Ph*`tK)j50GL zGojqsN|cPU`8}`m{rzzs=W)*Ibo2RK*Zci?J?DuLIEjT@6JmKTxU0L=k*+zu~R0*nHgb;v$ig1W_NrZ=!(*5e|f5i|tU; zaM{(OqO8Kk3%F(Z`1puYdHit7D|i3>qlK(yV_|mm(9zakUG55c<9|;cIB=jAd!z0< zR{>}_OfV(TdGY5k?w6E|{P@;Y?e$)V6mD{8=O=K&0fk8l(UbL6j!n@!8-qLs%nb#S z4y6BnXKs9W%9mr5%(pxz)bcgqZ=6u?-@muy%vx&A*w{vgws#?jY}9E@SeI9G1i?^x zlHMKvwZ{o)=ji^Iy3-tMAPa3h)U~F+Jcp~IvMPBrs%<>V-@bK%6%g8rt2p2>3wAuV zp;1$d+%zg2^XGk5RwjbqnFE1Tz!^d{$rMPq)I1U%JWJ!Y^o@~rF+Ox{OK8VyYAHrM zTvD4tbh4v>OrmLx z4-q3F1FL6w>$xs-;i&ls*B!s1r_wpR%2W;fCzvO2WBQTYq1gD^J%ZOXwCAcN6A;L#GXh zU_3tXM7wexeHpeK618SVIiZ$*U#lW1E322jl`u1Le!~dHgFlb08~m<6wr*n+lZEL+ zHTOA()K+|?UTNjvW8V);7({>qTE2L?xUUCcF+~PcEd=n~baXPv;6QL;U}9F5M1tV% z<92q;DoO28ZL!v{MBJ>ao7fKI+6%N9o&i6sa5T3qn={ho2W#Afeuc3fZ}>0qr2Ddo zPSyFZZ*JclYY4=$r{ZiY?i^0yym_C@FprYv9i-$Rj}o~!CgpIm;XHb7Qd&J=@w^sd=`nxp_d)SH-+yd6d%-7 zN8-9eMdBa();%|>`Gzw+R$ce=<7)v)>j+GCW$1>Kl-r%pxtN0<;E#WPt$RFF;=i{- zmV<&xlAb7FMY?mqQGmP%H+hOD+`HFhO`WKh7Shnr%F}NSpgFF9#drHt9>w)nBj|F2 zhYF+Ti+BIrk6hsO7zjOFvnur8&=9d*L-aL^e|sZANt>dT!}a|V3xZhAbEOk*wX0tb z%gcRmi#QqkQk0UFT+FFHgvo^J)~%B)-+D4LFTa)EA?8y5{$c(m>-Qe_l?zR`Xo+I> z0l+<=>6!-yp1~*Di(uv+uR6b{rp;yo8#UXUEA}Xk2U37L`w6l5qcQN}w>NhLc8owd z@e#+ESc#j^E9Ys#QFTJ~d96(wL#6Xinb3c~9K08H@rUktVcEGaEF=T>x75M`4fw6l{xWqKi?5`o>hWj-YiIoN%5sLFG~v){Pu|+T(Y<4Mt7|xFcat$9)PG4?Qu#5*BgLo* zyXhYGwe6*BVNM_CS{otWs{4uu#>Ut^1j4`9c21uF`fbA| zjU2TP?Z{E(oJM@AxX>-}N`|nCGZ5?6x!WratI1oNoPS>H3K{z3aFYX?GKCUyU|^tO zPkC5a(Xe&zrO&ySp)&&-r3bf)Q$NPH+D*;zC${5aWnCkGT!`C-kh$TS@sEpR z!ps%yS>RyW-F0W!&{-=zLfplcNU*~}jb<>h>gsB?ZIa8E#ayKQV#w(^Tl-|qRlcPD$sn^8eYMMX7;*e0m3*eW9? zKX3wIme#8%y0-uKx2H2Ush84(572W4H~c7!k#KV)Le7Ei35bY%hr)V(SUYda>S;5= z2e8;j@Eb|uxBCdnpKKM#jmfyOpTUWdnJ-2!RH6{7yzfB4(BkDg`ywxIwk$?FcHc11?XTyO&`~w`elr~>LWTEY)^C?JreU7!%kJ$Z zp-P;}9~^d)`8?hqyMBoO!^6)FGJ>~`dQDR?w`>ZSnuxvr!Y30iEOmu(2=05lB^ocU z9V6~TJZQvR4){Q+{Gg5+hxC4GhO<7)&qMHRGs0p-)GNT9H(%2eVYYnD9E0#_-psjBR=vx*M`1S!hfZG@j=5v*_fk!DbX(s% zJ!!f#;}E}@uZivH=c?hmW3}c-=S{&_35P#q3DBuN7tl2dA9esF8sdD2PQ&RPDduEB zWNE@Jgz2rmvGF1HU-fE_)y`WTHc_JI&N*(Gw~&xMiy?-$3P*DgyNzM2;xb%q&GOGz zLqbwiSqm&P|K)9qre7&Xi@tZL zobT4%%i;fP#!QFd;GW3cqipqKPwU}o&HcHZF1O^c60cTkKhY!47dzTYL zF!y&I|CWA&SP*@8?iEBNE@0c`oL~2A=ql%fI3zK;16*pxK97K2!T^yZ;daTFboUC; z76pAB7O%B)v1eo>b~QvCpATXE)(j|V8YmO4-mKVjGp|FDG5RR}EEy$b2I2wm90u%I zU8w6$A3{wF3@^rZcS0{qliW!-JKmkAuZii|e7vdKWz8AccRa8a_0O$~)OL4Aw8fql zsMHmjz0O-?+4iW)e`kWoo%%EfQ=4ZALbgAyWh?OM%N;%AG<5~H9pN3ahCY;-gt1d+ zTut6;yK*v1^(n71%jPiVaKa4$(GZ(n8VmEK*>Z)9~`AqZLXeMTak)yii5iYLJ7kz;^kyBjsL6EnRh$& zqW7~F$9qj)*!om}3A@4!W!;+*(MtGhypUoh>h_lC_Mmu*kw1$H|@d zTQLtYNqV??yn37#a9C}yh-|?O6zoUD54kN&D28NCKXHy*+o7N?9z{}Z{mw?1+0Wk} zR=%Ui$_BiF6@#(1WbdC370m||WgVpL!0kzbd2&=sco+=r5>%SsrGK|H4SBJ|&C+MHFh z`}RkB=5xr;{}QNAYWE1b*j6%d*BW&ub=NmOTx6y~%Mjx^$BCSkT3M8wSt{derMKw3whTyne+ZybRvferQO670`v_n)AlU>^ zCTP|o4f7L^TG}p5!C3%8$d>znBI|yGD3tCYaK4+CmXDu4?R)&CR{Y++Q@@0=4j>VQ z<-+}FXj)8=aRnpKOT;Dut;xVJLqM!pTcdb1vWy%;YE(`flFn4;+!=IiXPco~bZ3N2 z)Tall;`RUimk7ne0ebqsZBqJ{l=B*Yk4I)OdulLAj3`Xzo+pPHKv@5F%#928GKIyhzrKaD| zkHASFNpK4WLoa#;q>){3^{1Gp@}acB1ERj5`EC)-Qs8wtZ?TA-du{mC*rGXX*YE;# zcp%qyK>dzJzkr7i8($pm^Zs5m!_4|Z6%C^`c$o)z@Q*@)GjoA5oKrhBM8R& zdKf5=b8ryL2Zm_`OgMyW-ou6163ugv;Mts}+Ob*Vn|2I8arFrJpW|}eQsg3bN{g0( zv}1zC5>i|KL(m98|g{rT1@nR=V-1U`s@z&|KW*!tb&8{C;} zI`O@^y)Q?%#F+LkANA8#uZcmJFn+%;{86?SE1g@C@XbQ6s0~h6Sg8I4=M+^2f^2AkDkgO)RC;tBZ+uhrHIa9rHqt7sTxMo`tSI*i; zivyv)XMK(f{#o3)Nhd4zH!ybEZ)mOoKaLv?$@(I3DJFm>P(4d|eJ`ON4(7OEy%5}# zo5qbxb`ew5Gm13)BtVO!QF;k%7p_=^O+7Bs3m15d-dIPC^b(3wtb#h| z&jjZh;3=IK(|`Uv`S*dAfR+{o3>C7tn4XC`q)Bx49};GMlPJ*&$wPUG>0TYZqB}pX zHHEM$*RajHZ*vmwD^hE7F>DP_edNJr4p@leg%W3vL4t^X(h-Dbn!;smQ0);1%?qN! zdimRenjfYGJiVQ8E^;(fXcfn6B4JFsbB8ryQX{9X?|PWE)yqG=!~SyVqp&}Kk~xHv znW!#>wGqNC<#rzt#V7sv|8NszBMZnzL?Q%0S%~xNeOK3ky;zkXd3SoW^1Q;~HhvO>h*1LmDr4>fnvahNvrX zFK^m*p!%eBVT0(Fy=n68T{QPKpI6E;e=<~A1=*xjmy<4-PBwedJb(V&{=CzpvzbZ) zt+hU(NF7GMfc<%|=fdARH9Eha~KD5C|TmQVckEfjR_3*onjSCE)OF zgphj3J~(oL^zl z?UxAN(MniAL{_a!W_c%f5bZ$1NwRC)>6T4u z_|+Bxj&Uh59C}bh_`_{(C0uY`p;73G;u;pS=(9NsepeuFH&n)iN@e}uDac5PO9hkL z0;(yACmQ?~aWoTo?TAf@K7C{DtDv0Tw_j4*0Ts3KEWyxJ5tbA!EiL#l!f|FadZTQo z=Z~P(we~E+4LyPm?R5EW4Z3-?WYLHKIvx`%YD_YuA|kx!>{O9_|MFUDeZTkPxdyzc zq6j#D9n6qwY1>gSSge>NA+|Wx9+CUcJya7|swKWa{16YZEoc;(DLP`h{Wm+rBKoB9Fm#xue8!V zY(9#Nf!+51B%!7NdTFLU-n!YL(6iSOkXw*kdZ4b_t?i~$o^moW4bhV5q#%PRMN||` z9@4l*6cW#$1z|>zPAHCf2=_xcEKZJ7?TVPxAfl1N#kTFB63QHEkjro2b9+#7O&E+7 zF}LH97(m3}#;m!_12e|!ygaS4SCioIH^yy)6I|Y}>LqtX$o+looftq$zn zy_=ZkVB>K^M>|9Zp}HxeqN2t~WP`GN8c|%3RWK|Dm%A>jk`ARDjaa?(i5Hn2O1Kns zpa_}$ds%Hl^zkF5KN>H)ENfyRi-Vj168~fj@P9xI8yFan5H~0BN#Od1Mj{ALLA7AQ8G5n6QQR2+5GfAhNR0^ApXVj)TuypZK|KrD=2ia zTjdOs@Ussj?JzRR4nb3k@6O1{vjP$AE7QG>P7V=SomQ}l_Co)Jw?YMIjJc^o zfW6lfR{u}g4tMW{3WO`3`6K0R_|o|DBZZL!VQ$&E@dW>&>(YkC{DK8vNdDDcCy`a7 z=K2$#Ax`4G6PZPP-F;M$I)V72K(cILC7U-novdO{)D&12DQ@p3Du$9Fo;#~T*xdx` zWT6>T1f@_92k=h+dGOLAZno-hYSGi)In2(il8l zYJZGif`Qhm^L9ilnvGT;2udTw#{`>=6~CaY%o5};(Xp{5;IXWFzVOZ+(;`B-1dnhc z204#sOy}=w4*jYw+{or8M_sY2yE`#ZBqtlSn<&3M5F8L7ou_pWF+Lg#nmK)Jr`)(o z?hwD_{dbX{{DvOyfB3WaU4`al1VJrKuBXNF2GwI}M|v3@=3xt*K!p^%vJq1y)t<*X zo}(kl0SDYegg}}XnbE{GWos1Tt|puf2%rsf`abaEcSnbirRdjkfi}Bu4g7cX_K1px zZ+Eahdcrv0#@zAkaYDjjyL-P(&@KpCb8{eusRQOfI`#PYmp1(%G-Bk<&OQ2!$I zTLxIzK`(E2_wINvqZ$$G5j-9$$*wwTbeWcyzSfvDI|uk<=qg5Rq|C=LKQm$EVbIZLZ< zVN#?08tyos*78$E#pu@sdmObc&)j( zxB#tU(Go7cy_$KOq@~z_A;{5lv>XvW= zP=eQJFSh@EMn=N}3!d!KYIIt;I(5XYP;hgn3iG*)Tb)V742iwwrdrwQ*+KuF3L*D@ zzs&)YJC1DWX?7E{uk!QQvm3rpBvtN8u-~&5t?Ug+on`A)mBs`%R7-8dI2DbN-79zsnmH=kJ5&ZHGf6IdLGF)3aCN!Y1Y|6SH+Z4?wGsp+jTL12TIS;7wTDDdFTF_}+IO&IsC z!Cz&8fE(xtgoL9@vNei-)bHNP*4SGbOkujT)MCFT3|xPcbI*ZfKZH+AVhD>rzL$#a z!k9yj){s(cm}-ue!}6GH51&3)6_$#ZeKs(Mbzs#fEG~Wz*b%ARB;f+4au}uLfg5mZ zzh)uY9fpO!rTQ-&j6VGllw!78`X20i5sj9b(a*Bzg6|aed?9C(@-kqcVSC>iCFqrJ zKvPkxmfk?NKoQO8MO!TvCBFaJ?gykVXiNh}gVJt8#!PpbTSuUC}M|mN-k8}w7 z+U4iWkpN+mwRNj--nzaShwyD~8&W2ro_B6e#5YfuclepyFcQccIEJE^r$Wnj8&)p)>Ro0Z?>vk&1@7QS~KV~v?iXEeAR z9ox?9DJqDp@tcTZZK;4Ky8wq6e4e7K?w)${4Lq zwlN}f1W#++nVY0|E*ZvYWPS7WQ#DRBsXlK3Fj|PfwfBCzE6hwU@ixYx+uf z-`8ZmJ&vORR{KBy{x!9>1`?5?By(}3?Jc*qTyL4opz-9-Luac(l$dm9#gQz-Q=yLh76iq|~yUNZ;mAfvnxp|+} z9mnO!&O~vp-}cFwTzZP+JIt9!fCSis^Ac2h7i+zbjEra?)x?kXnN;McwBPH+@5fB# zfSETzm+`=~@;1}#a3{~nlUQ&dCtF!Q63XU>>%nRA1O2-odJOLZ;E*7j1vY%Z((?Pj z3DxA&QE{=_xrx6bInTJf9%4Ip#kJ+)#Yl6{{a9PZVE-aYK7enu$vsu#DuNVmodKwy zO#z|Y6Sl4g6D>J3ZBX$*DEGyQdsQN~3Jm;#d`nS$B^SSZDJ`&hI%M}r!L@Kdn^Mye znM-XsF z>|)MY`@vU&(Lzq(97D6*@#Cu>U*CwDd6saP_y49@OtS4CKQf^)9WPHk1j@J`aN*Y> zO^%F;BC21)eUFzH!S@3%>#6p1e_k^FR+hqACioGw(=%|Jp;Z7>qsLmXSsRQQEjp*Q zuKeKB6aS-Iu@r>b1{DbGdBKw2o7{i zc}>t1JpFTU{}DQ}w87%>JJeKEnE*hLa2Senv_aiz>SL5An^+_^GNl+k^e4<@>zBRW z<#qqQ=FCtWjH*qQXE=8NjopvHvZYFIX=Xt zaT`=R2SN*qp50oUztck%i;t;{UcaXOl&eEyXLOdC3YQq-=!A51bt8lro2x5y8j?1}Z|}m-pN&Eb0PY)>Z>7O5FI3gi?2K6M{Ti#wQgNJ_G{Ut zuLI?3qqeBG zKcHenjc5{Pi7v?V=C=6-j*9nI zHr%=?yf*a=jb24FhS2ftixcz|x^a`bizW4>=B}O{&GfRXpBxnqO?Yv_vE7QjHp(nq zACqSH&c0FiaN*al*xtW&A=fQxv5~Sf)Kw7|RT%90X6To2=^=-5e~Pbx@oB^85O9Qo zt)E|UGG)Af%y~_>O7Fw%V#~{YuV`C-f8k_GoqO(#F0nXaYg%Yja5FLib31I zi&mz>G3R~1I0AQ-gMBtb{}eopP0ahPcAhO|28D)XYG5ExdT0pln`~SB)O!%w8ECOu zGi6P!J_w!p^9$}V7zbMBv3Cdvo55S~&SN!_$ba}Q>BTelXKic@2WAOky}>h?*f1w9 z!>Vxs1HyPk8qp*TY2G#iokQ%#z}NGoNk-Z*O+gE`%q&eB+$dn8one{c`ydyglIbJpG{kS-eQMs}Db zL;yN^dSlbmJcmO=WZsJ2s2qso?0S7e{LP)Oa={FzL>a$dpnVYIy;8G%vBI&3yUa23 zc5nZaJMqOuk+Oj~w{EXGB1+UL|BimD`Ow?FJNJHgGgxfOL*JL9Y`yV;KEi;MXMX07 zitEasV|WLTp^q7Gav+qDK<72?FaJ)ox6e4Ge3-F;!X9&dFkA^uO-;~SKcf8Nlx-^i zQuiQC^B5HxpwE^GM2tqt(~|_*4k-Y8y57e{{W4MTm+4WQHE^M@JX-Soby=y32=zuMZYe z{FPCIt7{JxOAA&IqNvs7h$P)D$E#O0_wFyk908vTc%86!@{K`)hVf^EQ;9}eKgetd zl$7wO`uksIXUVo}jx}?n^4;O{d83Y5hn4k@7p3C;r*WT%&|U;$$8-1_gzqWA_Ka&xT!d*b`CWpl0ex&Vm zkX3^blL`|bLEnP_eK?dHep4Ck=(O*$uChG|BeH?v9tI3@eRlrc{p(hH(c$hMjQgVp z7J38<5p=PnaGU}&HgtZUG;RskKR!4D3ot){eEl&{c>IFoPAWJdt0sHHKf44*g>cvEWp5wVXf| zW14Hmkq1DK6VD)lt->Iw27ez8vdee@vahB3u?Ts4ck>}DC-TI1q)`*H2Vm=wMtL#3 zz;20;63R^gFC0aILF}Ji2hopjo(oDVaGX>EFjx<%I{^`^rwe%Dp1WYY&#LrhC5er! z+I{aab43%9Eo#uRAX`WB@6;=%fot3gkXLc`zm>#P;4!pVG6AWMTI=O4Jsu;g8|j^p5yp~~L5WtJ zek6#j`bAWybK{_8(xPAG3L;&9Yy9*fr6?>$5Duh5PQVy7~>?H z|1U3ckMsE1)cx1#Fsu-D#v(?Tf+4KU0#*Ty()A4c>Rk*B z>=93Ws?(#Q{^xsufd%?Ip*O2CZFQlh7WR0;@ys_9%zklcL`1o}yGDWWMVzgtSn4=} zkk}P*_Lf!{izFZC<{m1^7KG-?*dvS(9Lpngo#+SZfI_bWC?0|Iww?cW$ZjPh9O%yRCFetZ3k3yBcybao$om-22$p z;`mtvc*)Dl6Y5)RPZi$Vl6XiiWvZ)zP6H@=8MF+LzaQTZSCj-NI*d06%v=yDK^`(; ziAGnAnCqXT2$~Du`#G29=ci~MJi{g!0dTkJmf5o2_gwcw=0SRAZ}1D%o9bDV7jNFg zff<4wQe8GL>vI4NN87UqTz4DKcmCZ;fP_E-?UTi0&3=6mFrL*J8XAHPE}V~%$z=RXj(9uw3J8SAxae3Mx*d2`jhH{s!CV;)~-#RBqddWFVhg zsp9H8zU~w4!#i+<;6j1z`=Gma^*(GvnUJ_BW(oZjxvVZ_RpZgl8{~(t9TV(021<;b z46A*TbVwhwW7kLIf89Id~yAj%Ypq2Af{{5D4{`;PgnzWO8;H#qJ zQGWnV=!3g;D)*(mgKt#gX?+2>-o&~Js>ej=(%d&ic>-Ix|E&g8KdyEsWDf12W>LU} z-kLVa#C2q=W&TyS6^@BoSdWn8hY$b)KY`{E5`(dyKR*zH1K_wY{7Rpd+Yog%Gdm6F=&>cxfI-KA#7N>Ki=Kvq{tRDmsZV6!T8B~KMu2bcrs1TJt>Plq|bb=;ruNvKE6{-NLC}> z%R;B`s^F?V)NHOFq@kDnui6uK=zT;=CO5B=CfN~*14(PCmmRB4>S<0)VC0l51~nTt z4^otXn(@*bgT?q!<+k(v0z^uUU(C|Si?*!B2S?~Ox!in@1f_X= zc&CF6Y6tO}VQCBu2vC^4vC0nZZVv{Vt4?{)fNixXE-H^AHx_-H_`=7*Dg%LMgC3gT zsd2CocqJavbU8&a4v}TjKE>nr0TO{M`U+sK!8IPQTPKz!x%E|=Ce6NnH4RU-N~u+h{*yo zT_HhB;`|_9E@(L7ysJDIIp*iTz#}P}r zxLkv6@J3GCm}+Za<cTcj^FU%4)qc6~EKI8fQ{if+3| z7VaR*4a=ak2i{Yhl}DF_Ae18jTar7E*X;bj(+Yuhy>V`NCgRinfkRJpOZ8N+j6#^U zYtNoBlqFyUM9waetqMOLrWrzAfiUJK;A4oLAV!r4c3Eze>)=_)#1C3n8EZrlD(B3; zy)0E2rSgM`sXyFSNxIU;rBYDol?a(C@nuDNXFj}ZG7A1$6K%a9 zto##MgR*&ML~a{Hst2fxD5O%u?f#_qep(vGPFpjMg)sMj2k7V@0;t6AX$)uW4>n24 z)29t`Gy}2gI4w@8j5f&M>c7v-jMJMQ2{VL})MR+QgEjR+d;H&v#UtNSx~)*g(b(1& z3}wXs1qp=ehRDjw(l9Xv%B5e#%ie8ujc&I}$@$liH^jzIR^5+~PG<|0CHu zgT{2_&l5>bNIGzC{%@eeEC5%B30p`&IG@7S#{3luM4q&7Z!Jq<5#sUCm!-)YR!ww8 zg$O($xEM}AHy}0TzCntsAQTYO9|RRdCprylV2k_*V7bFLVbvbFO{1d z^J*xSVAVUB_6zqHDE943!Sxw(SR@S7Ky$V4o4UD%!2R)Iv_#tm9`;Fm7DNgtj6jIQ z4FpeXXlTd*#Ar(53S3Z_9iEUamDCvxwRf$Hoch*zfE8e&;uG!2dX^mSv^Q!cFNM!l zyB%M+edLJARwN1bGr&Zqc>Y6}&j#UoQDze7Rq6=U&d;bgKxs^@rg#&uz6I;f((~zr6ACA|0tRmN|DCyX!qNVysHG}Q zx`dfexj8h-=gqZa2xqA!daDpZXx~ho8gnHn1o49yIW2JqvYi{KtO&IIhpuG{!37!* zr*L=LeLIhk3Zy*_O1rY%aEsQPgzF|gst^p6UX{ucWXJIxKoE<{hKPZ;3ty$xq9H%j z)7Nh{dsMQ+4-eXvdkP6mN^ydm_*l39#(9Lw3BGbT8zNOde;<-f$PAIL7XoXs=Sy1q z=}3jehHn09lt|+aY5wv>FJr;aPe%W2rLd)~Z5Sk%c-5pv+cqUYOhOx}ghPWi*7xbffD(~&i$o`C8 zVc?R4q-8*)%VPYr!GXpgJD$Mfx-56jUepf)zK&mw-OxfIfwjlNf6U@x$=tzZerbHu zgwjU9e<=9S!*3eZ5}L31wr`aNOVKP5Vn*UFH(MF+{$|l0a92jgb)GP3X!|xQaCsw%5vvdCY{V<5WP%#nX)))KEjDj%pLILouiuFMf} zE^1R=dL2EbRv@st8CF!zlKuoKi99mAR*lf-gCA6xEiWnZi{xyaJzgo$4gE3pbU`Fu zBE{_3%KZll~2hT{N4Q2rS=N zH@P)m9%v3@p^1vxv2&!6f=PdlN25U&5yTkzOdwQ^=F^SDY8|1LL+UiX=W1^z(U?Ps z)@z1_QoKCBJ$w*wta0k&ryTu5HvJ6aPsbWT`|#p5GKLxko6p4i{X?b7<^TANK9pX* z^)ga)&$(MDFiq!e+sQ=rZs6XT`UG$2L z%zh^+Dba&KpbVgEfCBw_f4ey6R_;NqB~yq?iNg%Ay7;-;zAx`P5*|a+=fJYYy3@c1 zuUZejK3Hoyu@%vJMV{NKUHWR*nM`G69}4wHpL?UYihs;#;Zh$!O%qa#Xx^IoNVrXu z|NC)_g1lie7cWW)RzWZSyz1(gf@UAdJ=c5^|La`&dA_0XtxHk)pG9=lntvTEjx5xs zYi*I04b(`AFqhtUAL`fVJ7JNH?|N2TxGFR!t@I+u5$R-^J=u{+dotG?T zOn40s7?3-`E+ppA=CV*T@8VQy+2u9)LB!+j_R&h`h?>J9O+B{p=jd=cm%klWjU;+;(<_N42xfX082iDEZH zUoFk@hM!8Z2g2m#BKahqqtYi9pdd`=%Bl!2qS$!J9nZjX30J~&Urua@JR2La-QjB| ztkyQ3p514$zNDoQMiu;+=CI>$Y%ULzo%n$}zx48O<@JuZ5#x$C{$`IvY@^q+vWf}{ zOP}sOy-!x@ZGg>xde!~EzNU5UaQs!_A@!PT9G%Gy3R`vn>Vt2G=M>^Cmtu2G_ga`J#C` zD?vp0+F_bw-ZdmhTgALZjXW3R3MgOx)vzT=y7Oim16fwu2}$fYkEvNn@Y0t#KZig3 z`->}Ny!6m$5||YRxJEmJ^-wg(eZZ5Ejl60A)T`&i3gwQAxVx`XrJwS3!Jon5QjXC4 z*_E=wh(O->wnGrO!?|pYWNn<&$}5|liBZF^+?)og^tp-~iZ@igjQ?G``|!ilOtpOz zW%l;rw=0v5tgI=;7@Uk)jyn7R_zNL+zSaMl6EiRF?H=&cD~lKH3wCvGt|?yUjNd{z zR1;9`g5u(F+^Z993EHpsvYz*epn-dVIF>=C`r`&{6*VmMJG`H*2R!y`!z!iG_PEiX z*I#{13QbD-X4jm!&4_W@)J#r3lvx&EtfZk@4_D`FL;9fib^auT)f#_QAwoXd-`699 zpB?qvxEZlm>pQ>CH@Jv5^6sfq+1GOzFm%&wIzTX;CS+a@2nCu@*ll&OIh4JnyYHqr zCBh*SKM=$~KUn-JdlKgCNZIh%Sk1+)HRVs|fmScNZ+ltwDRpd?*=6=DuLwl*t$uuN za55z#jXVCTq<|!X1>IpgHO4Ij0|Cds*Bf+LdQ-TSLcW=RxQD}via-qzj)A3T`}4^e zHlK{GN$>V&eP4c7zHQr$svX@s+S%-!1&}mt+_`n~;_pKXYhyqoh(HA(!;IsGF)3rd zIXPDX_EI+p9mDR3NtAF5;HWzyZaw>jlVyZcnS`iS0;f?Z{78yO>Ic~vj7*Ctyt_mM zBoXigP>6u)@nRz%^7r_H!AY$aGdg)CoTIiG3K4#I`P&}lXg@5vR%YAo1g;96=@(cM z@QnT*7ymK0{~%ZG!leC1`IT=U6M}p3$r^FqSYlnV3+bOD;zPi=ngO{ZpeX1paCZ5l zVRz$JmE3Wr2}&joA@RF5l)qZwdBmEj4xpJJMK0ojo!WEx(Fj1^GTR#jhKAFP)8A!* zs+Lgy~*snSX5zCNv=k2jXxWf1K|NOAEA@L?Oi0tDTHIRe&lZ{?8}n zpPg;(abfx;x31=u_+9A!8C@_(sz7Jp!4oh*kd{~6oSVDa_p^`CZgZZy?eqGETLLEq zkZO=ADh%Y7o#c5|Gg{KGtAGl8M_$4tj2#5or8uqQ?G*l&iJRm(smYSVU0M^P#}KN) zK3X-L6uh3<^*_S1J|0L``k3lV@d(WoQ&5|6ABSzu*yKi5||^sgC}*W%-j`0Y!oE+go1gjzdk@o&R zrmX|;Ss?@WBPJQF%A8OOrlh1itHFqxh2Bge(&VQ2+`$tW$4=d~LEt3_VsoGmD7Dl) zqD`l34Aj231!=eMu?pZh6hy9Y=Uhf!-o{y>DB2%`#YL`zRr>Lt@6E-Wo?Ae4D6d1l z)7__BS~dfFtZ*3L^JveWk=b_H|M@Fw@lG|Cp5ayE8nr znbe&-<~5}gf9F(gEUa@|Z7wf=krc+#Jf?6Q=whe0?3gWR&$+S62Sl zW5dU{H9&F*28}1cn+VDji`i9z?uF$VauY%=P565tnRY zFdi+$B|ARK)W+u-o>)rv=C&Gb1e4_|W4H5wZE3TPk zY}`{xxD>;O@bWPHbcB+N%s#!3HQHkpc~vCXpv=Gg8BL23G})Di)*d*$y`6n12mvuc zY~W7An*Rb@AIP?3Y^hDpK4Jck1SaOt^hJ-<;LssNax%ok5F#R&!BYZb3};Z_2#Sc_ zr5NpjEeD9812}h~A&+}@orzIF)_ig5*(O2Z;Rv4qW`_m(DX<=bHix)V&Yu13GfzsR z$ksdQ+#jkt@3`-;ykaX1*N84iipcfRgNHHq&;dySY|WAB*xdQxw9wL1m|IP8FcLT} zreLls&!n*}OWB_N`Rmlh>uoWSjjo-OR+!reD;p$<4Ke@H(k9q{|2p;WnKWKgbT`)a z^ek9iX6?L9~B|LYh&^jkNhmai0B6BG>@K($){z=~$S9iE;)s zp`dnr-7ZTdk#{`b{wDkT_u1E9P-8c=xU{*3ieTHr=Lm8Ln$(M{KN_n|tQ58{`r||% zhj^#3=51cyu8$7Cx@lf-V>~0;Y>7;37<~yGmf*)A?0~oHGJv<%@MG^gbH^Y&7=4|_ zJtjFmwnNXoU65fbN9eHlv4J!;Y;ehpV`2EwbQ1Cpy^-OWm@{$@I0TKccB!Vu% z2AM+0262hU+Sp5O$NB$?IJMSzJKUS7?2DlxST}}eggwuY^+10y_@>?(Z}0yetR5l$ z2Ed}vMC^`WJ&FwFH>zKmuea37#)I=7@D?uXCYuGrD&8F-Zw{Vk1}j7Cg$N8H#J@J2 z#u5<$3O(a}Lqifxsh}GKn>K)I_+b`fGz~$g6#`EbU{$(eCqp=@k#kG|WX)qwFSu*xG3VQiPfj+zR7*z@ZsdQewd1I$87*eB7%8fD-; z_g~b^?cOZ)9m>WI%$TsTt5J2vxZN1jGDB4wrODNQen# z03{`57aY}f7{bMhUurkhfszF}MgkDk3{#KUHfi=N zBOi_nhL{#sd@tVwT4KT%TFTIjD~AJ-+fr|)MMQN3SJu{gbnbr-+js*_un^@A>vBT` z2SSCgAM*n71UL^bI1&s9UKF7jNE!B}IWEU> z&q7ti&{R^9=6vVq`?S9Ewk&ifxgODUbda65&+fAUER75}4LoOOFrdG{ueFd;R`#rF zc7WzMjSJ+l1ek^|kHeUX;vlnsKL6E?P2%3f{gdt$NPFp|QiLCdkYcp%VFcSe(Rp|{ z6pt_Aj>E*E3q`}y&xeBdh3{}(%|b9y`<>wL1-noUT&dUO+Oc`MnH{7c&vnDWQ-88GsHjc z0i=!)L$NRS8yr2(F5_qQevh?}a;>d;0@c6eWkQ2izuC*qt0bsaS0@nv4$7nvT=vd0 zpMRAtR$ir)fWioesQvd2*48(-gxtIT)e=Jezx_tUJW>5cPCc)vsmh&wUvjPJGOT&P zmQnc>0o!~x#GOR#0@eZVzjF)2%=tm=Yt4;8tz+GF7oCI9cr^b1iu>+=uG_zTxkOS# z2qil!q!dXZyfZ>ZDivi*MrD?nkUdH?BwIsD6cLfVDM{HRQW@Fvdz_!^b3eZKe{lE1 zbv;}cdcVebp3mp;JjR+md!n81%T3&J+|0~eig%6J6Oqjq4>(T8XH!G>k0v!r^H;9= z20hjQ^c;j0NlFYfORC|P6J~NwzL4G$o+bKwHE3^A@8&t%K6_=I0U(O-u#27VFw&Rx z^d#aygatn<#VBjNW!+?BQzQ|ZfS4whq01B}b>*egd6Fdta!>U$eT`UgvGuFH=SW$A z=q~KGA=#P@p@<3?%W*}Yyd)Z*#1ZEu& z9+LBfXm0e!gfImJ0-}0y2xGK$0|f>kc%t)Pl>0I?w7c(Kyx<%(t+D{K2>?eJuY;Gi zHg8wFqp{do5YIb`;tWqR9794_H@|i4$zJ8azsvE_V5#Fr^GM3%rrn*qy_tq6-)?~U z=!8{uav_K!;P)|vEwO$H%P2vMV!;k+PY6&_dcz8 zrBX*ubMZSF+U^1Yzsz^<*8Q<_$yU7OPHHS*xaZJ?-3DHd=cxyQ{2V=qp=zc4Z?zOL z{KR{YT={BLF4eJnGK`8hE-ia(<j(QMV9%3kD2mw4>`3uv!{cjE z!9f5+2eMO&bB@J1`-bkRk-Pisxm<;HqL#ky9sWdeEB?ct*FX^n!9DVv^P{qVQi1vPL0zWz_|uU!RoPM zY~*>Mkd@ZniA|Q@?(;L^-wl2KE$Cl(;=Hc7Z3#t~PQWT>+|^BV_|-%nLz;q)mOn?N zp}AlfD0JJ@&_J>QkbjqqUuRmm+Qvrx(D!>NKXUv!SYs?+>?byVWyObunr4HM_v4HoV%8Q zHd`YnhNEe^gmrqq|JKIWBcC2`M)qkGZC=*t7R-bpC4vv_%LY4nPic}Ck1?)pW1?_J z+fxy4`DP+*ydxp5b3<5{<5uFtvWI>q1Y9BC+tHUu@oc{HVo z#g&V5)z=?iwoJ+xY!F{)oBpXG6%ZbNhR!I9MOcA>A=;<@;n~K9vxS0kCU10DxoJ1} zhOnh^jn8B~I33Ajx6!j$zPgQSpOTcCo{6gC6@AvU_dgh+sanvbEkNZ$vf5VJG@d0= zhg1Bqqo8bEN?{fKk48Veq@;wfKR{f^O`mvcyqWF9^J+npyeGpdTz*0ExgD{Xi?k z9ERu5x4!?BN^c|Me)SqTWl;-@Eo2#YXD2tmF;yfcVyH6@Nm-|VpaWG1h?qZmT9T$J z=f%}-Rw05IPaA)%?p*gEGSti!o{25_2Lt^G&FM&DMsv2D} zPary1WY25_m8zT{;ir9_6uQ0Jw{ItIcz7Pzz*$!-n3-crPvx$QVLv{oc+iARtrT8Q zqyj@x>Cl-^xRmN~lHBM-lOT|XUH_^|+|x(ANec1z>o^=bee0F<%itMZ*uqtMwWJUC4k)kIx!hbfcY4@N)t}W4+dRw zz!z`xZJhtDU@Xi%ILlqk1|WR|-Cq~VYJdrmtrg#S)9 zJ9}$xjBNJNHYmXys+@vhM?|>inRbC9xD1ofES5RhX`=93l}ue0nfuX zuzLQX>d6fzMVo5z#Ds&U?h z3S6e%K@a{V?e`0lmLQ*y7|b&EkMf>uf4bk*;VOt#SD+Xhm!jg&Tw4`%m>kEGu8-Tu z+M9uCTc4WNHU5%T^W zfBg2z6R@uvOkVhPxB9BuHv?}*e7^+@y}%-iLj&r4QnJv;bcaT)m`N{STd~RrZY_*( zoWrdY*BO5;10^{1(x>k2ItkbNKA1n;SS78s z_0bv#sJLg07ysIjQ6{R2{=va=KpB~k0f#O7QgJTs+jlH`JAs(7bf!kzBY~&YUb!cO zse}7t5#mS+2F>HgBf1JtR=M#gG3_LsQK8i~ zr1UDYlB>sDSBtBuE0xh>fiy{oE>J2E_0ofcgc|VhafgM$=1bLzi@==D3`RI$Dw=Uo z@~$=+Z(pve=k5OS(GLEqL1RppOQ=`x^3LKKQ9Y_z&y8<}OP4y9)#H;B2E_s;8)K$Y zYnryb;tKCPta5(Th(l%_B&=Xu3u@9Gf4bgrjG-5}20^SrL&4uUANBB&TIBYajS=R; z*Wkp0j&Gzn=|g+XosscVXnnbXVPWrP)mn8Qr_zk2j$IFfr;b;`CMG5mS%-V7h8p)#u~wRLNi9ji8|Zx}Yt&kU9(sU-Uzdlr98 z`}M_O>n?`F+VLRP5&i~tZ|Id2UZ42huVyRrwQvg+XgNgOjhiz3F(Yz(D3t9Vb8~9y ze?(8e-igyh=R$;g$w4xTgpge}`|Jn71+D_;j@{VXLgCM_z%7pi2@!}9Kx#y49$58a z!G5xigkgsl8rVSN2yTGv5h7d;m30z9CWTIm@toz0o>wJU z)y~631B}!oWB?EE@56DW%fYDMh+>52a*4Mh!ozNjalTqviE#m*Y_mqDJOvV$upikDU!rg84_AR0sc$} z-n+~>;1D6-49`a5<6_~tQzC6!f~oO-L`4w?P}pM=bqPJ4aR|RAhp4{#mg!jWjV~1pOxrpAU}O~&oRsHN@Pa9N5YTo;5a+N7Cqc*bp6?JFxUC7wAJB=SW7g9DX#U*)g_^KX(c+OmZ{^NwV!{fTT% zTkdkz@^7b82VdLYvX7I}qGeWnK3qaAC3PzGf{a;cmdNhFl)~YTJgwcX9)BypTAe=W zxv5uPPHxx5FT2)nk*vpE2dGU8m0}E)?9QY2jtKA%7HbB?DD1~pm<015FoTa$pRu!Z zUtk?7b9g>aMYZ*bhN`k*qng1>e$qE$q6bXqCv|7y4Tw!hOw#9Ho%9FXU|z&V2=p*r z`*P^2;15cg=CGP7aah9vq=vR2KZ{%Z?jEfjrn*cK+oA_4>No4@c^fVShb6*(_TQmY z`N4|M`&`e3aS7e?;!TKiaA=LTbCCSDFvHd_ukMPu)S!Q^XJQ2R`5!t3*=HNZ{yUjE zS>$j5O1{~(^AEJ*2Cw8jtzXMvQKB>hQ`nh8^=PKzkN0b1JKk7FO1@uSi-Jz>kWa*` zna=r;9#M^iYZ}}X%A~CBc6y3ITF628irqY?u&#kug;wN*&eTR_>VfBHzptcg^V7#3 zko^>q{@ce~0MmJ5A%s+0VoMb-Z4Tz0Buogs<5f51apoh8Pc}(4py=tD|VoAgG+r@61Rpk62DQk~TRx6h5AFZT8L+ z&6TB?O}1Ue*Lw=xgnAS98#InwXe-Yz?>sXzP|rypk)gx(M78ctj_(5V&-JE$f_u)d zs6!a>5;GGF>>w6fxkBIv+?0}6Ha@KnWp2HX&x4B^%3t!@fWZm+ew5T^PRC}YoSN^o z72kko5r;mQtbV+-;4Mw+QP9mjR+b8O>aeQiWY@*1v?NP6g@%P$V48>VSpWznA>*E{ z?|r<9xiWg2{3|~gPdq&j_#68NDtCA})n zP!1WJsOC+EV&mh!p7Z#PV1Dnq)UOIwMkApn<>sLb41)wmMbL1?Y$q9 z61brE25tG=^UvzM#$PXRCCHcxXC$#lKx0M#2vh>*h-xX=!mu&*2-97Wp28QjYuD?#voE3j+je6)b8JzBO z=*+C^-LV_*Fi((`I^Na|_s#)I?>*#lJ$gl!YpK&iS8{MG7m%u9)qx zkP{2vhN;?Nv+!xmeG(6OJpr~^Q7x%Pp;X!oUU0b#_abwd%ZLibZb70YkzasKe|&0c z5X%}f`wE*ucKFYdnP{$!9x{K_(q}7;6lm7ct7YnN+lplAsW`OYwFmvCIt|~V?twZw0GBlm{7{#?I{$%|3?BfvSLp25=DG-;GkSe8TVWB> zRM#)<;GFpdx=8=KpwLhTJl`SA#c)z>;^0t6Ix!i+!>mFOEMEYG57A2?x8*1{dlIk; zhS{Q$@CjI0Dkq%z?In>rilzfSJrG z3CrRQkuD{P1q9CoB~F@u36n8Fu*ity5~3IesTn(`sZ5!9%LXy|Pxobg#t#4vj>QCt z`UMC;CmBqmtlmtGBsrSOaBakjw!jMjq4Ts`R_n4fB>I+k379RWyWqhHaZA< zI$!q%Odmhze}mB0Mfp!Vf}|nIni&q;9>alNrx7D}7zaGa$3%1pHFa!eW+umv)avT$ z_-DsD>yP12*tHm=Pv4-EmXeY}t3ecdxOgH~#a2sFI|w@uzQ5Aa((^Mtzq$|Fst{g~f8VIrw$fSUH_n~)B9q8S$MwLRzfnO1_za)zeNx1%wdHz$>;Lw6`IOPLI(|B_> zpqT!F9gA2a;f0IB#0W*Z8um#IbYd9X{C^AWlz0erPVw1W4KF`yUhFV~kt^2G{P8Ql z?#m`$TNnzb$rx?dqRl?1EMmIipP>jPC7EQ7ewJ97N75KD*#nHS_AF_EIhM+4peuev zJ_KyN`q?)3k!8PSMaYhPEOBq zH2Qnn*A|Gy{(Qr9Sqm!>6~cg9*v}|Ic*X|+qQ%5U6B<>sf1`16=fq(sCC$mFqVnNC zcRC7W4??p9FOG;aF%rRlhQkvG-BpqH3hr~qc2oHHScbjJCuh;aW0Z=pt;2waeU8Z2 z@q}XJL?oiyVXOk6&?2RzWL)-@<)RsbX`XR7L5LqJ^x+1>t(gn<&3G3k{{9WESr4gN z1r$4{`0!m2m!-6LMzLM>VKih70C^?Vn|k@ne#G-My4ce#eh5l#Y)%ii8Df``yliP) zkez=3=1+LN3ALHz>|C9)L8-c_@X_>NfA}q!7vsL`-$40-{GprJE3jv+cCoPg5b4JM z*l(OnJakyOz?nWQ98dxiVM*b>!RK0T1#GT>GX<)9`-b(RKQ843)1KWis|8>O4`%=s z$=eMMEP+W$e7FINp8oo+r9O78OD3w;O!51)cbH^FT%smV3JSawAhx%4c?(lj`UkK? zQMEbY2$|!YvHi%49Fk&9W2B|F^_cS68=n#p2~i(D=p$*qga5FzYF%7h;0spKN2>HK z;@onZwsBY7b^5%sIDtwQ`YG?Xc;|44yLQb_~MSBVi~{yXRzu1iOh8-CQ9SK~T@%v4;C6XN$hqnK#7N-eyh- zs*9^w;OxJ@CrwK@P@*cb%4ka7jeVhB_)*;;1(s1Hv~Cj%3yX^DXFgjMQ1?LH#h^E_ zC#h7#);f^IuKr0|zpd$Gmik1gr?b%lieGQ12E@f}K?RNM$(D5Ot-E(TG=0lxdK06h z9Q$wk0Od|6bx*G>?=H3dC#9UH>FvaNNPI72(*YZl$iHUCdNS}%!g;6_%`H>j#jeiO z^#+e)EJ9D`QG0peWB@cuBtcjPb8b)GFIzLSu~nduA+iP{L%YPuA!3+PN{R^Lg)bxZyc6Dq}E0`NDMYRJGViPjoxpC zk9$q|6Bu*kyJ(59KwPdcI5b2-4}qc)4bi^DJ3=4%G>DBCb}T}10P&fr(2b$XvcEk3 z2e>4p-bGXtS}+=PD=Cm+LGCtk;`Hv(>z%WGV`EKNR1@Rl%`Ka~&50)f8b!>+?n9_3 zx|Fd)+oNju;;#UP{>y3EQZhGAzv=yT3RAy_&xd|&0!amDGFXx8P&1{1uBqr5=Q1}> zq|x9{f{{M@VSQM#Iz`sb$oA^J{hQab#n9dnY%Y(vv1yM^fk328=8J3>dQ$))O2eh#^VLb^ip$9zRmcTU%FQ-Zc2ebCTx4gZE|p z>3qMv=H}-o#`{!odiug~1QfA|Q^!E-B{eN#p|BQ13o7>dH*2oGlbW6_9z?n`HVGGH za1Q^a)J|}$G2>p@M&30rOxCmSjQ>&?vxy+u>DucLUR9d?v&$pL!nQIZ_PyE6#U%C( zJd&jq6``%Xug|Ud^@VRCuHcGF*V{>-r$#Sx{ReAG?B0-`4XvZc!jvc$K?Y_g@*Xb~ zso*?!eG(FgL!?-?;} zzg=8u(ZBl$E+}E%*iNj@d zM~vC@XsBr2wa$58Un-10>agvKnieadh5KmfxVnk-1H@ARl6H_TM7MHpdvrp|mfhfX z_Vd>wA&22K0ZPt9J#TJJ)ocTAcid9wwGv!J$*$^zecHh*;^Kg41738#$sTv1gkRsb z-)t3J3p%^6G9&$?-%bqb0*|w$mzC`(cFj_M;QVWy_PuJ;SWU)1jR}+c54eRKI^9#; z8k-R6G319!OUEE`akDW_@^K4L2OWJr$F8IMU8roXB#Tjid{bpv?tA*A|iy0fL zFK3(7U~E#;ovu0E5%_{9fM`EqPL0i(NRTk=Y|*~sLzFumZ(YQ}AKL;4l6hALmvg)~ zfyO9Yh8P4NGHZ5_FjmqiKi2qBE2dLe6*tn7QrtGh6sK^a8C1yBS>=KYKa ziJ4C}lF-v`c35qssgmt2iK}U|kpZodl|a`^OSoVkoG!Ga8OLF=Bkh@KOhBxRwe|_~A+QJd$1SnrV#L1?3{946ak+DUD@S3LIW|7Q zoX{$E{M@1|y=HA(oLEEXmUHF;--k>6)xJuAM8~d01)lLy#I)e_A?^-vbFD#<>gK(d zJ4TFE4G;+}A~2Ddp+gQbo~bETx1uXQR6jfP9FlXQs*T+oVtnpVb?CCflAH5oxTL(y z$w82c2gI(dC{4jG;URi`w0k6A^*>|lKf~jk`2_{iuD?zq@>UDTvf0J^qis(!7@Ss} zh{GQ?g=kz2y1@$Y*Ksx8OaOR&W`6~@al%L{=Y{UWb(ZVSwF;grI{8QTf}r2c`}e=R zxK8u#v4K|F(2#abbsKV6B~d!_(B$LAeP&wnkhGTQs<9A{U=R;Ww-i9TA7KAMKvy)E z$?&4W)aL=XZXUa?SUrC9DO*Zx18>pqYIT4%B`Y!ZcJ1riHV!IWo7W$v3{=&<7B~WH zy#htAx5ow^B%H5d>2QS?Q@N6ehO)Am+{8<-mCd^gxWX|rkOl3m5+^JUk>Hn3pZZD| zaGR#4rwIcPVxprLCAt))_Dpf-BPS8OsUcKZ*oTPW9V!8cq9E1EklQho0|-~ z5M0oq*W2#Nq1K(QW`<)sj8sHaD?&prXw<~8*c5q*-|c8VMt`Qnz^VyV5uWTNkPy+c zssY;k*Bx8%yoNOQiaU_FOK{x+5(jDreO4~R)27mIr>T1`cD(;;{GjR6g}!NH<4Vcb z_4VflD5qYmrK`JwJYLEt92qJ_pd}9TvA$&$UD@`Q>=7_MA{_|~8HjQPCS*8*-o7^q zVj>AX_>wSvpo!wxx>dN+LGZT*3|u;Jz7lV8fMx~o3qUm`49+x>po0GeMw;rRX;V)t zA3Z+Pd_C#_AsGVmh@A-u@=E=*V7Yo`spVEW$S|Nkr;-6I~)eGAbhRhxwTSn1;Yt{)#`@t+70~^n8Ml zt8wKd8?dEC+Nx9jz<~jHY=B~ejesB98QF4B)QuiJ;XEe-7YlrbPK21DFgbeks9B{Q zxN1Jfn-@_`;n#?xwUNHPCeN5#fPwkMH-Z3A7C+ zfWL8apjaS^LqZo*3ETi6mydGGKvwzZxCBg3us;CHb;kyYF_O^EonfR^1}}u;{QQ6C z&je|L=>&yBA;umQ6RiW@c>{)~ri$HjUWx{tXmkJ(AT5|fSNMb-4YON4K51oVe?$lR zRP?3b+Y9}dsYdv1R`uY)9Ipy!@4L^uZCgE1a>DrrST#2{S3WF**z{ohy_NCGB&r3_ zTIKKoyC0{sD9>FVRB{*6hqt5MYX&vrUzYku-^M~+I&QPRG!e}Y<=w%N%Cc52&3xOP z2aeNmdrJw*S^2P^Ri@*1ZvLNtuxxz0&+31Dy&^Y!{mlRA*GoL#88c`--)t Date: Wed, 30 Jun 2021 09:54:37 +0100 Subject: [PATCH 2/4] Removed some draft text --- docs/source/modules/vert_coordinate.rst | 4 ---- 1 file changed, 4 deletions(-) diff --git a/docs/source/modules/vert_coordinate.rst b/docs/source/modules/vert_coordinate.rst index 337bc41d6..d111dc590 100644 --- a/docs/source/modules/vert_coordinate.rst +++ b/docs/source/modules/vert_coordinate.rst @@ -62,10 +62,6 @@ Even Sigma ^^^^^^^^^^ This option will create an evenly spaced set of levels from sigma level 0 (top) to 1 (surface), scaled by the surface pressure set in ``spectral_dynamics``. To use this scheme, set ``'vert_coord_option': 'even_sigma'`` and also provide a value for ``num_levels``. Increasing the number of levels will decrease the thickness of each model layer. -For example, if surf_res = 1 then you get even spacing in height. If surf_res = 0 then you get a height depending on zeta**exponent. -For surf_res in between, you get a mix of the two. -scale_heights sets the model top height, and exponent determines how heavily biased your level spacings are towards the troposphere. Larger exponent values are more tropospherically biased. - Uneven Sigma ^^^^^^^^^^^^ The use of evenly-spaced model levels is often not the best option for modelling atmospheres, since it can lead to low resolution near the surface and in the troposphere, where the atmosphere is often the most dynamic). In addition, poor resolution near the model top can in certain configurations be unstable, since breaking gravity waves deposit significant amounts of energy into the upper layers. From e48391e34ad6e9e3514c731b9e5c1cfeda621863 Mon Sep 17 00:00:00 2001 From: daw538 Date: Mon, 2 Aug 2021 17:45:41 +0100 Subject: [PATCH 3/4] Updated following comments --- docs/source/modules/index.rst | 1 + docs/source/modules/vert_coordinate.rst | 16 ++++++++++------ 2 files changed, 11 insertions(+), 6 deletions(-) diff --git a/docs/source/modules/index.rst b/docs/source/modules/index.rst index f44b75efe..061819d2a 100644 --- a/docs/source/modules/index.rst +++ b/docs/source/modules/index.rst @@ -10,6 +10,7 @@ Components of Isca introduction dynamics physics + constants output idealised_moist_phys two_stream_gray_rad diff --git a/docs/source/modules/vert_coordinate.rst b/docs/source/modules/vert_coordinate.rst index d111dc590..e74004c59 100644 --- a/docs/source/modules/vert_coordinate.rst +++ b/docs/source/modules/vert_coordinate.rst @@ -5,25 +5,29 @@ Vertical Coordinates Summary ------- -Climate models split the atmosphere into a number of different levels in order to calculate various quantities (e.g. temperature :math:`T`, zonal wind :math:`u`) throughout the depth of the model. There are a number of ways one may wish to divide the atmosphere; this is the purpose of the ``vert_coordinate`` module. Isca provides four different methods for prescribing model levels: +Climate models split the atmosphere into a number of different levels in order to calculate various quantities (e.g. temperature :math:`T`, zonal wind :math:`u`) throughout the depth of the model. There are a number of ways one may wish to divide the atmosphere; this is the purpose of the ``vert_coordinate`` module. Isca provides four different methods for prescribing model levels, which will be detailed later: 1. Manually prescribed levels 2. Even Sigma levels 3. Uneven Sigma levels 4. Hybrid Scheme -which we will elaborate on later, but first a note about model levels in climate models. +Glancing at any of the output files produced by Isca, one will notice that there are two different pressure levels listed. ``p_half`` represents the upper/lower boundaries of each model layer whilst ``p_full`` is the pressure height of the mid-point of each layer. Each of Isca's four level schemes produce an array of ``p_half`` values, from which the corresponding ``p_full`` are calculated. + Pressure or Sigma Levels?? ^^^^^^^^^^^^^^^^^^^^^^^^^^ -Whilst at first glance one might think the best way to divide up the atmosphere is using pressure levels (i.e. 1000-900hPa, 900-800hPa etc.), it is far more common in climate models to use what are known as sigma (:math:`\sigma`) levels, which are scaled pressure levels. Doing so is advantageous when topography is included in the model, as the different model levels form continuous functions; there are no discontinuities where the terrain intercepts the pressure levels. The figure below compares how conventional pressure levels intercept a mountain in the model, whereas the two schemes using :math:`\sigma`-levels follow the curvature of the terrain. All three models have 10 levels, but note how the distribution changes when using the uneven scheme. +Whilst at first glance one might think the best way to divide up the atmosphere is using pressure levels (i.e. 1000-900hPa, 900-800hPa etc.), it is far more common in climate models to use what are known as sigma (:math:`\sigma`) levels, which are scaled pressure levels. Doing so is advantageous when topography is included in the model, as the different model levels form continuous functions; there are no discontinuities where the terrain intercepts the pressure levels. + +The figure below compares a conventional pressure level based scheme against two different :math:`\sigma`-level schemes (even/uneven), with black lines and red dashes representing ``p_half`` and ``p_full`` values respectively. Note how the pressure levels intercept the mountain, whereas the two schemes using :math:`\sigma`-levels follow the curvature of the terrain. All three models have 10 levels. .. figure:: model_levels.png :scale: 100 -.. note:: Isca models will usually contain two pressure variables in their output files. ``p_half`` represents the boundaries of each model layer (black lines above) whereas ``p_full`` is the pressure height of the mid-point of each layer, denoted by the dashed red lines above. +.. note:: The model levels in any output files will follow the scheme selected in the namelist. If you wish to interpret these results in a pressure coordinate basis, it is necessary to interpolate back to standard pressure levels. If your model has no topography then :math:`\sigma`-levels are equivalent to pressure levels, so this is not required. *(Interpolation documentation will follow in future)* + Namelist options ---------------- @@ -64,13 +68,13 @@ This option will create an evenly spaced set of levels from sigma level 0 (top) Uneven Sigma ^^^^^^^^^^^^ -The use of evenly-spaced model levels is often not the best option for modelling atmospheres, since it can lead to low resolution near the surface and in the troposphere, where the atmosphere is often the most dynamic). In addition, poor resolution near the model top can in certain configurations be unstable, since breaking gravity waves deposit significant amounts of energy into the upper layers. +The use of evenly-spaced model levels is often not the best option for atmospheric modelling, since the vertical structure of planetary atmospheres does not tend to fall into regions of even thickness. For example if one splits an Earth-like (surface pressure 1000hPa) model into 10 even regions, the entire stratosphere, which contains the Ozone Layer, is represented by only two levels. In addition, poor resolution near the model top can in certain configurations be unstable, since breaking gravity waves deposit significant amounts of energy into the upper layers; leading to violation of the `CFL criterion `_. An improvement one can make is to use an uneven sigma scheme, where the model levels are not equally spaced, selected with ``'vert_coord_option': 'uneven_sigma'``. In addition to a number of model levels, there are two additional arguments in the namelist that control the spacing of the levels. ``scale_heights`` sets the overall depth of the model, whilst ``surf_res`` changes the distribution of the levels and ``exponent`` controls the bias of these levels towards the surface. When ``surf_res=1`` levels are evenly spaced, and when ``surf_res=0`` model levels follow a function set by the exponent (with mixed behaviour at intermediate values). Setting ``exponent`` to higher values increases the amount of levels in the troposphere. Hybrid Scheme ^^^^^^^^^^^^^ -In models with particularly high top levels (for example models including the stratosphere or of other planets), running with a conventional sigma level scheme can be undesirable since the model can become quite unstable; calculations in the uppermost layers are more stable if the layers are flat. To overcome this we can specify a hybrid model level scheme, which combines :math:`\sigma`-levels in the lower atmosphere with :math:`p`-levels in the upper atmosphere (and a transition region in between). +In models with particularly high top levels (for example models including the stratosphere and/or of other planets), running with a conventional sigma level scheme can be undesirable since the model can become unstable as a result of the surface topography being 'imprinted' on all model levels. Calculations in the uppermost layers are more stable if the layers are flat. To overcome this we can specify a hybrid model level scheme, which combines :math:`\sigma`-levels in the lower atmosphere with :math:`p`-levels in the upper atmosphere (and a transition region in between). To use this scheme we set ``'vert_coord_option': 'hybrid'``, along with any namelist controls we previously used for the uneven sigma scheme, but we must also include two additional parameters to the namelist, ``p_press`` and ``p_sigma``. These two variables are values between 0 and 1, where 0 is the top :math:`\sigma`-level and 1 the bottom (surface) :math:`\sigma`-level, that represent the boundaries of the transition region between pure-pressure and pure-sigma levels (where ``p_press`` < ``p_sigma``). From b840e1ca9caca81dbf9346e4d0883a37e33d34ec Mon Sep 17 00:00:00 2001 From: daw538 <40756463+daw538@users.noreply.github.com> Date: Fri, 10 Dec 2021 10:44:20 +0000 Subject: [PATCH 4/4] Update vert_coordinate.rst --- docs/source/modules/vert_coordinate.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/docs/source/modules/vert_coordinate.rst b/docs/source/modules/vert_coordinate.rst index e74004c59..478b723bf 100644 --- a/docs/source/modules/vert_coordinate.rst +++ b/docs/source/modules/vert_coordinate.rst @@ -91,4 +91,4 @@ References Authors ------- -This documentation was written by Daniel Williams, peer reviewed by Stephen Thomson +This documentation was written by Daniel Williams, peer reviewed by Stephen Thomson and quality controlled by Neil Lewis.