From 78de49faeab9cdf8f46f9aa6dcb4c57b02b81c5c Mon Sep 17 00:00:00 2001 From: Bill Ladwig Date: Mon, 19 Mar 2018 11:10:49 -0600 Subject: [PATCH] Updated CI test file for new CTT changes. Updated docs for raw ctt routine. --- src/wrf/computation.py | 16 ++++++++++------ test/ci_tests/ci_result_file.nc | Bin 2148317 -> 2148537 bytes 2 files changed, 10 insertions(+), 6 deletions(-) diff --git a/src/wrf/computation.py b/src/wrf/computation.py index 164d905..641b164 100644 --- a/src/wrf/computation.py +++ b/src/wrf/computation.py @@ -1111,12 +1111,16 @@ def ctt(pres_hpa, tkel, qv, qcld, height, terrain, qice=None, True. Default is :data:`wrf.default_fill(numpy.float64)`. - opt_thresh (:obj:`float`, optional): The required amount of optical - depth (looking from top down) required to trigger a cloud top - temperature calculation. Values less than this threshold will be - treated as no cloud areas. For thin cirrus, a value of .1 might be - appropriate. For large CB clouds, 1000.0 might be appropriate. - Default is 1.0. + opt_thresh (:obj:`float`, optional): The amount of optical + depth (integrated from top down) required to trigger a cloud top + temperature calculation. The cloud top temperature is calculated at + the vertical level where this threshold is met. Vertical columns + with less than this threshold will be treated as cloud free areas. + In general, the larger the value is for this + threshold, the lower the altitude will be for the cloud top + temperature calculation, and therefore higher cloud top + temperature values. Default is 1.0, which should be sufficient for + most users. meta (:obj:`bool`): Set to False to disable metadata and return :class:`numpy.ndarray` instead of diff --git a/test/ci_tests/ci_result_file.nc b/test/ci_tests/ci_result_file.nc index 2cf9ffb9fafc9c006f3fcf9369eb985a5f187729..e15d100a4b5f5d5a4ec7e9526024de31d652a014 100644 GIT binary patch delta 11726 zcmb_i2{cvT_kSY~&ofWYT&55Tk&J~5nVL!FF*BJ`9wMPa$|aPckO+yFM21A>B4h}e zDl-kr^ncHLr(fUiTL1rl{nz^6wYc}Q_de&IyZ1hO?|04_j=v1AOuP(FkU&&6mrDA> z5iUhZS?~@NJTA;P$qz9ApG`7G?nZe$i9C};K=Sa#^yzwUoK5$gZpZwYs=+Ondnd%I5AWP zG1!3Y&Ot+3bZB`fei`9}Qu2g$CE0L$6J1vjh8C!*ry`ZaXGTa~ID(c2dQVa^;g&U$ zG8_~mwQr1zqr0oU>-pdaTFPkPD9X!WFr>Tf`41RCWui7fx_f!f(gP|ppn7>Yo3{We zMIA60+5+kIfXeI5P${Cyp-bV3Ox$TRRK#90$08V5N~2rxq>^yn;P=NgLq zSv*PSq&%RkNdbLm7>mY0lGaFel-%C65RT;1ir+?yOw1+5e0Q3(84VNOCwtGgtrllo+@_y(#1DXRBMJw;MJsKp7E$)wn zk3y4J(s1hdgJ2G_TSYSE4R+pj4Uiv5o&9Wb*ioe zR~PrKtqP$nA2?ij?P5HMxiThaHN_EoT=ZH~5Sq*JwkdneYi|P)QOhh0;mSCB%&43} zrrcf5QT|pXv$ZrBWKeh3OHV$1 zG9f?pw*mT-vry1&7C8Z%^m;{!Yub_pLanF=c4<}d2IJa^*)A5PgW5LtrOlogFci+p z*swCG>am*~eWv*qWfYTbVB2S|2hHHnf`lb8+nbu7Qlv+f#acv+ckzjjBZ_0O#&0E~ zcS$zxjW|-uhg9fFZ)@i6y521!B4zVbGch|`GFv23QrK6d(@YYbYD15b9=+L8fGUcW zHZC<~yfGzr!_nMa=NzBhV}l)rnER5ZJH+)hm0O#Y7h~-uMf161lJ*|@q-;`<-qK8< z?;shipjKFlWsxI^Cwv;-s%hq`&*7>QKE>>!FQ#nn6h?E6MTgG5jwd&Py(d18BU$wH zK?esGLz$Dox(^&eP{!JEQ?+(~RDGC@akhwg7V5?o z>=@ekW-N-ZCu)kD^>=KZ=ruEm&#Dn)Xv^b!Tn;DK3AreN+>6wyzbm4T84CV;N>^6BUx-drdsld(>?mANiZJ9UH$k1J} z2_kwir}kJfM0_$odOJAiu)fHv-LVSodf5ldp1~YBf(`8Cu=ZK5=D2KIWyH7$r;yA1 z=^ID`&9t2^0!_7->syt@;Ec~|V+jgTx>|MK*BW_r9OUT4Zn6hkqI4`wgZc$g9S=1X zCw7jC7)$WoVyXxep`$mfj|dY#qU%8;EFo7+e@)SlMO7chDrWhbo3kN{Q(D(P7TZ=g zdz1bZjkAbpr#_x7e?5iXJ&aSQC4V!aq99vGwB=UBBW0asQ|TLtE_aXVggX#i!-Y{S z7N$`=lXQsw=uqLH2o^nEvCK%hysln8F^Q-3i>&SWg|&6Cx*mE~-57*|gj@$JpV3u% zxSL*>{J~FhGy>cwl7sB#@{ghec=Qi?=$&K_!qVp2XuHV8Yk8KO$(*#2HNDv^)%|66 zR51#nd~{ceW5h&|#2)S`(<*NQqV-)#f{^Q-i1MAazL#G|J`^z*pfSnywsN%P6Arzb zWFEu%w!+4nB`4e+J+7AslM(W|k|fH_WpjMY-Zoo=Zb#jTlfpK8b*@BLYe|$anwMPl z8@=hEkGNc}=~Ev@rxI^_S>oXF;5(3>4W9d!jzv@kT%L9fei9)nDIHr$^1p_z#Skds z+&zq~zAjo>Kh6|(ns(HZ({{h6I4cT`q#+)19lhHeprFq3l08G<-Z#$a`JYTv!( zn&nWLy{I|TpHn8w_+`pinsAJp+Xd09d<|QAEi457qfmPQ0nb&jRl8WyzGi4bCx8=`_V`6WVsre+-X?zR{)3VNoPJZ5`9_y;=&>)6}2 z2IFMgcgN3NIP_UOdxvFcx9^Kj*(&!^z39yx!{|3q5}nG=U^aQRUG-Bri{T?RkJ%?~ z?hBbwdR*sknVK*VpC>klD|daj__GSj8D>!t_LE}L9BWHcPZYab z4V|XZUj?_Ol*QWv*01s$&W*lsZ%BU7a;@pxoUvQCVIe(VoFK3Jor(&-=fCHU%6wio zJ?j!MT^?=j?x5XM_A03_uhEKtNm0mcyl`MJ1G&?>O7QuYbI+3Qu<}+o9zv-U;pJw=k~>34R%9pG}YzQ>Tl~ zf0i`VIvug_bIlL_=&3pXw8i&WJu1bLIKQ`0zA9$O6z@^Ufda$o+M zYoZ(+bbtRfUpM|To23^M#^D?uDuiA`e}(g;ms&iV^Y=7#`i|Qg^|!q_HnyR0wg1

cTa_oN3;*Un3;4Xwh{>)*33b=I_PeV685ejT0Ud;c)QFnZwBnaT1e{5GKniu&#j1R4f@BU}oZ zvAVA0QBe`F%M8c(d&Bz2aMiPuFV6MYPrMj;T;uWXt)X*T@8GYqN1}?mJIat&SAT`- z4qj2S_B~#CVssrTTh^!SB_{P@7fiyl%gBe0#byaisiiY z`tF&$<83pF7pI7CS>_hJFL~|ru)5H+`d-Oy?G?Yij8AnSo_~GckeUDPh6A5kE3#5Xy6)@WE}5mmA-6)0-cUm#;gtG!^=vDpuEPi;;?O^W^=za=hSJ#*@4-=h%C&kX_YEpbywPQm?8?Y*u%7j|7% zbT4-~V7I5N#%y^gyF-{s#purcweglKUy?2prX!^_3WFBfJ!S`@@dpUwH7h>7ieIaT zzj=B#TCHW3^`!@txARnn7JEDW3~Ku1=4{fGt1dUN>K?Lq{bx#L;+M%{&$Gt6`3L-3 zMi=&0yIeH>)_&yY&j)eIs05jpg4rdX4l~VjReOz>|D1TbTP;HTJKxH+&sTqrypyn* zsaYl9(hTF8?nxC5O7;4Ayy`ec+f;2a`5`~CZlB+x#fVl&jkRBudh+d$dkkD#@#3nO z`O2T;ZLvN@vVnM&{Ja$(Zi!2XlzS=n#oO5E?vLSOblV(*gOTH43Wh5P$`a_(`a zw|{f$pAKhvS>ccqg>>5|+RI%%H>FbJ-Z9fo%bKiaqn^1>~$Z5f&SQ=z$6Z7h&iC3e&FC(1- z`~(tr5Go6-WPaaG>@*O}sI>THA8f8Yx%5b0Hkr5}pW1n`I)7mz_>AxC@5z0*n6L6p zP9_5jo9d^&p0-{Y>zKE{v(_Paar)@%cXF4$r622_D<1yv^84t0iLEH&X7i7q!{-rF z>)e&+l&0&{&z-FFcAJoU#B^A$+^;Lc)@kAWlF9Q$0v9p)%Y53o+@e{Z!Jj45xB7ac z)1y7g4`__`wS^@v?MY;uKi~gzrSOzUs6VC;VuK9X*N~-kR zP8a!$<8){UdBo|8J}=n;PxJV--E>OV?p8bxkSe}CbEc|BN~5d$@o>+S+mphO?|nKG zZ``V^@e0>Eyy-vmdyK^$@qF#6BBB`jeo5y|QQvoqjN!4llodGJ#{7Z6nUZ;-?jKKQ ze5yQ7NuHr?*)1NqnlW%ZarN!o!69pZ28NE zu+-OJwq5O4Uivcd4OCb9=M*uhJ|W8by1wZ9u2?x)b(V#uPFr~-ebYVZ(qd0~PP^xa z=X}9s0|UnCRp+jBc5U8$UMTNe-9Qtl=P|OW8~M2Wy4_xnmmJ6a`flQ;IC*Rvw#wy$ z4qB!(H8;^-#1r_^HmmF>d%^;D%cRn%t7L?ZFGbHMRf|tlEYf{mJ9BtOU_t&$IU~c% z**?Q1-r24#i>kQ*QNh-Z)i(`v>x_qtr}Bay&3DV_0V`I{e=HzlgRcno4zN1eU9P6O-W6^{Jyns z{8f3@b?xy8n#y23Uq4@84YfBy`v-_amfu7VNXf5!=@vJ+tBQWp+DiL2^Ruvq$z`Dz zlWOleb@83jc3S3y`8QxG~d13D_`D+arvdm#K~x=O}LVhaY;>4@%w;y zgZ`wZ5%!Su+T z)~6;xJFmj^!{LsJk92I_b|u`Ybu(1PTQ=}AEI+!btMKN-)klPM6Hi~i#6V512p7EW z?(YSKiks!}Uw^{{VG-7AYORkr!Cv0_U&yF(Mq z0?%66TC~?eLT<-kvw|=7Wiv)u5g$HeY2Hn9cO{`{{o-{z@}kac1MddMEiB4qUz#HE z4%hCTBEPP(?e^3ou={u^n3w3qCZx1S-HAPUAyT%h#`UCTCuUTZqfNIc7Vn?@G~sSC zV+27AFHi}qve<8?DIrkQ7Ip7P+N1F&tD5iMw&XNkR@P?EHn(e*L?3qNx-~~^yWXKV!t*LcO)#n97{at-(t~uB7=b2 zioB8%OK(<@?710nham51)^#MuQ`dWL@*X|@MT<5^4fpNVn^WS)IJrL4vP57{9Hz(L zxzv5nGwu;bn^=yhajV=3w@e3kyyc1DbvCQ(#qZv-cggIz!;I4^*U+oc+CnC zwdj)1lBwb7ZfjmF;kxVY9fq3ut(&oZO^uCBTv<2z&bd6yjXHEp_=&YCe~p-^5dMXM zkn>!3HlyC>)q~pW){7?=O~QGv(^>SnjAx~8VZFW-6lCr`{#(Xs>ZjW!g_FVWX9@K? zpIzz2<2Nh(-3Dun=9+7dedBfsQU2{#ZV;w#$ez3@D>?tQZ_tyVnfO5=Wcca~aUl|~ z(abgO&2M)6=-$Th=e^&# za2!fyg26LuiujKR*}aW9@{I8p5AOBdagM25eJ%e)27h|148L-5qkVXq5nYU!hr7a+ z@B+J@9Xe8yrK{h*I1)U(X$L>$alYtPoO-nyS)d-SYIqNF`F6BC4eOh#L`tfP#TR90SAK=h~?|X8k6A zr9~&0X?a$-X)#^ktnzKFmsF+yam4Puil;Ph;@a$@9oC#F)(Y&n{ zq<6$jXurJOlht2wTo|iU(Y+R-=i>IX(=aacE37s(8Zr;Z<42Z2Q=QhC_C~)Bc_%woQA5x#npU?C=6CE+1f`oGT|U zx#4FSD}V=DaMd(;90XmqokhTucm>P1w?TSR<*$fS*_ymlgfOusWT$C3eCa zwfSj6gU&U4#-G4c+!vb)oA|t`EAdf0bO58hz5SGcXpYUJ>g zh2!ic#P$Q6x{eH%l_^vh?9W<0mH%fIpNjvi-&66769XUQeH zJ>ERqbVTU_RM(1rjIzz;TBR&lFYnu4nBI>-c0#2?7&i5@rr1qNe_+v{+B@t@CuI`= z;GWih;wLVifu!A%lMqt}I-338t<7cErzs2mKRVEX2-i`_s}p?); zDUGX(20|fYM`(y<`*2%y#FCFRJ$y(T8rnP(rYA!G~7-*~S1|&tOM5MoUK#nV9+K=Jp0h=pulo*wZ_7CuIJih^j_hZ|JpH`|A1jf`8(F9QvCQmEz~)Hd8{GiV&qQk4SXOK!WPph zS^(IN18iqeD6Vb1hbdFw0TuPQ|(G;rMMcNuK<1h_Vfpa{TXYsP!}~J1W)4 zO`Z}o@XBnIfjehlMNs1@b|c%{Q(GBW3FK2}B1On{j)@|4=^jZ)Fex`;dr~7!7|~Fo zHS#=i2OZ0HoI@@nbciLH*jOF}XCqlPn{-5#SJ)&@#KY<6&H*g3bAXX>1|g|MG~I{0 z&yVP`B0nJLiK}OE{s@G=9~u6=(Px>?O?h*z8%~7h)F6-eKc_4T24k@({se8f;kfX9 zZj{*wNA03&18wi*YipmL56QUWgb-p)kf}S46Jgv8xwzv5_rG~V-QY#-rpjXi9DUEb z+S|G>wE&h3l|@Hp+539=`Ke7pweC1R_!Knkj>97kk3zrPal8btG0G%^qxMppU|l^) zbK2Y5-7&z?T~I}6n+!+EQ`xiu6#n4@fTKX=AcL%fGt~f9kxE4b1@kBY;yx;ombB4p zP1@E)a1j7Y3CO#`aL@v?Mw? zIRQ-eroS;A34gcno&<4&1j*O(8<#4oOGR=}Z1&0 zZa%m`Lj=qH>M&DvmKaI@##YCNj zO{+}Ard8&bOh`L^mVX6)M1uij!ol(w6_`+T4T5vB2rv$FQ$m4&q)C$=QJRzG9ty8E z3F0E<8=UvH$G~7u<+8OYGpEs#; zw}_Zx=(iI`vK+7gS^IxmfE)|{#^iGW3ytTeDZAR4h1_(|uro)I8t_QU$H8Dz9cIqk z&mfY6paOz8zz-OPqoL1V5J!;+*d)Z~%F&6Uabvb48;%n~pX1nvbi0e)p80fR78I{R zKt7pRRg}+nY#wp#9LGDjNP0fVY?SwxzB{Wzej<;Yd||KR7%5>S+Y3htBd}-u6n99n zIi*Z?8X}t~hZ_gdB7uXfI5?L0>jhp$Bf6x6gILnR8KaEkVCTvsaj?;`LbH9uCclVV z;6F!fsM4So2aQ234lt<20S1}{?**C#7-$w?pjj#=z(6$sIt0}Kz(Dn&15gbBKBxu& z26ZIBz*L|;uoMLoU@%}RAOK7SFetD91_c(tz~{jWf$0Drm=0iIx@}B)KMXiM5CEnG zlK_|!@IjOT7^F9Vr>6cH;{POD>IgqTZRFEUq$|iOy1tm3SYOO@xd0&?NIYzT|A9sH zxj`DvcuqP{TVTnD&meziJo6tCZkzPL881PuFgnT7Voj7J%eu`%yO5*-iG7=-=z^Dk zFF{r=crHAsH9)weRM{}tn>SFr3*HJrDlisY@SJ#1L4Xeea)K(5petSmLAsl`;&t($ zasUN6kV6oN#U#EC@@+MqIB$a7U?hgdbL`cK7GKAc{&wIceNX<5WYW#F);Kwtn zjZ@yobM5QLkHokJ{u~4!$eTdP)o&%YiyFj~W(8#Mn*TWBl_!+>n!f`j&2++nTorIa zSDW}t5&qoJfoA@QKM|PxLyO!Ah^yqESbt;kv<+SyzYond^B+fQzUGIrUh_B8DJYZZ zCq{O?gg9DA^kF;;g{~+^qIaxyg+O&CNO36a7LGRBt(8B|tyKVnIQ*ID(JgQki%55Y z>P7@QP_I_I2ce`@9mh<-{y72#R{=}igvBflOOBG)(MvwH5^UeaXye28@7j_5Wm;3GuJl3U&}JNSy7f zJPyZCE{u*Nk0^t;ldJuXSkY8+RcJ<3Qx~G+MDDN_qAXtMTM0F>!)J*senN1%KbaZW z@Be!(2;?8oKgd5+O!%{!0|Y<<0vO~UplpzT00#L7$bto8UZn17AireCkPQ|h7lVK#M*eFMsb8~ XEisxXq=rFcHs_HZFwpBJA+G-dFRTi{ delta 11922 zcmbt)2{@GB`?uK`Ok=IFgq*9i$hpZuF-<4fi`pOitR)l)8BugP{3Zv{>wk%PW zgruTuAyJZd#?0y0?|Z%P^?%>%eXfgVKKFf}^PJ~i&V4`UoToR#VPo;(umo}VSXc;RLbWmRhzgAKdk`;R7hDRrjz|X^2lk^e$kjG0jg_4K} zZobk%G#%xdIfvc8Z_zFfx7M{MD0{sFIH#UFJ8!Q2+jGVl}DEYSAjy&+E8zE$AoVs ziAX8)&G}GYUy{ZK_x8Q-&BH^bqC%x$q);jNXftXpo7@F7>57+y5|-i4$fkeJrX|@C zd`Rvqa9vB}@-c@~q$)iGFAQFRf?kt#O;oi-P)CK_LF*mu=Ir6_=zcy}hJiXB7?S#U zIy!`h{rM;~pt8^!AUwRiQkeji6;QoBT`zS2Dpf0VbTA*Tc0gk%OK$gp_VIJ}_rKuh zT+jn07$DBT>_~G42ox4-<7F#DUD`nY^-e}mh%Ew#l+J&B9=ei=ISIY9MFb*-s!>p* zA669$n@TgL4j})S*I-nE9fDL~hg3>~v*T~72p}qUprA#6Aw0>hihqXjrX32(Ou;n4 z+k~M$dQ=JWg(toeYUE`jp#omfgNOQ{QO?NtH3f6Xpc_R(o!p>t_CwDZQM|}&P3LEl z3?*@poxBhu1ZP4OD=oTHeUPRJx}7fnjVpcnc(-{#`9RS9wn`86DPQOrsL$UIjWCm+ zvZZ?a;r7Vj2koFT3@R1rX}cqcVp^zv2TGS{8e_{6ydKEm(sk21AgbKFs6)@m_-VX$ znKt1=UF6v^4Ydam7*Xr9gKS5%ZgeXdm6_}1i@y3+TsKiSQBtMe$4#{30$L3 zx~4E2pC%*j`2ajblZky#bXdzWUk~ej$z?fF^oYCco(>T`CS>mlPGx_2CR2KZKk~{H zJW7|<($4)R=YE^}IKN>Q#Ij8S)Qd&oEsn$B9e0>wWTXPj1~>TF^F(jsMG7AIIPLlB zdby?(O<_ct1ZsR$TkOPBdSKf3l8>=p9wN`%z2O-k&XQFDcTDOAPzHRNnL$sT+Mw(&4~p zzG6z z)RC!szsL)`fi5A(P$n4J%v_c4Aag@99Ok~R&${P+M#}+0RHKcMP0h*UYe!TfkKb2q zz&S8UYfHY;cv2J{li2x^uxly#1$&K9^Gwq)tKqTmT&^|^m0NDGqWxoW=?E7{uKsk^ti}C*I8M+bd4EQ6`@gjoIC(nJ63|_wEEU6U$n7!x8-b$cxDI7RIaL z-y87>A}aRU@~mH(l3w0{p)__9P=l!t&>13hjJ#5gn7GM4X3m`b8uTCS0<}?9je#O3 zdbLy8^EdA>oxRMZ(V4$xSXq#*Cfb=CRw%Er+$VYCrd!4-jVnYa_bWS*Y?_8w`9Ctk zUquD)3<_gAp}9Nr>Z3fFLB8GMk6xY^$+ zsJ@q>w~;B3T&lS~Sp(;~DUx+qM1n~%GTiu)uv|o2=Ao*6GJZZ>;@ThZxhikHw64%y zH(?4$5jrAZiMet_jwqGh@K}ydL_9V3-5lEIZHv>HVV%T}_DfDo)25t*Z0jrOFLaWl zOHko)gdGkN1g%)48|=b?bjikNm$R7!k-VbX4an0LS1)(txCIvIvlmR8%bPo0|PZD=@6V^x>j=C5i9(jNe^wx-@4wB*E02 zc#X*?kuwYCgO&dnE7|0GSz83n#EHqk+SIB}y^AnDxxkoYdZ=WtVnkT+`*&(ymIuXO~VcYD4x z)HJATCT5jgYwrwEQyAtQRI)Sn8qu(x-Q&Ng&QtesPewNK4p9`Z>+P$1OC+e@S2Es+ zuj}ep!bZQhZ(fd^RK~US&SHc1_Qt(d43Xht(u`3lj;8Y8*sX+NElvAVABTv(zu#W= zu#(Gc*6J7YXgfR|1tUqDd{x*PLD#D~!@pxnk}U~28Ji~5A>U@F23-_jWKwWi+tE9f zN#E`Gq9uoEX9Mwoxv?|$kai3Af*t3#EAL&CosfbGC%%d0{# zSG*weRzm)T%yLx~y61?vs6&CO>H~~kNZ{wf+wlhMv(x3ft3Hn##9zO73n9jFBX-qe zpsP>Rd+%w^f9@X}cAoNYcBq_w9Io4x9xqlRckss^vF8D&^qn`2(_$HZKKE2OGvacI zw_Q(Z>J!0i*mZyi<4ohMSbX_ls$N=W!*ls9_k?HJ<7h#SjT+nKeF~SI1IR<0sDaosHIVH=lT(lD9av;qH9) zS;Ebyp;krPgCoyuJVjmMEBXcMJUa6 z&6OwjP*;louLikOMETdGm8^rXQ0`hKGWYebXV*iIKkNQ(n|*Kgs4VZ46_*z_erPLr zao5eti|Rdv=LDX=8W7_Rice1ss2}YK3Am{9wqZPQ^lE(?>+h@~&FOn_Z`Q{rjy?K) zn!})%RIBWVACvS5Vc7X_tZ#`~Nta%96jS?E%O>j5E<~(??k3T!cQ&BmXEs9V>4eE_ z>%LD?LGped4m0FP{p|Iqa$m@be|p3I32ve(R%O?2oYoD2^aCGzy9)KIbEB@yua*zj zH-GNdz9%*Ewx<96`u6|>&mmTRmQ|wf)8oI&UHmi#9yJqRD-;agUx=&!EHE`Hxl?h` z1)D)($sLyQKe0|QbAy$kJ@X*)KJOJ4QZ~2m1EVb9*M7s537<|-3AIZ)Y#MB zK8)pE`VwEgY77-xA(@oX<*Ks$9y` zJjSNGn;=>DnAmCL{pFvc@r`n!SC#W3&}5JP4XK<^%bPw=XDpKFO9vJsy91uOp6MG2 zbqgP{d0NRkr>j(by`0!(xn}A+3v0bp8^+{2r}Qbb_g&ZX+}bwm3ZHm)pg`Zw zDUa~A$lxphdiK~p}~^*&1mgQ3<)>X==HE0PZZLA+~1tJ9Nyup9sFfb zVx`2MSA)5i&Lz${(mR-V&-7{SbkBufK^JBdgPAXxdCF(Xi+#2n-I3!`GX8xw#sB)3 zgT#u>VsE{~iqXzK6NR{J_cMtRN!6~SCRPf!QUjER`i4c1uJKND*OeXd&e-(*aq(8# z=C2rQL4}u92A#!io8zWGO4CjTY)GUu?Om82mmsdyp6XxHxEapXQ$f796g_h_KDIo) zW@TtaDbU39g|hW)j}vdc>KQGGeeb=M z)^7K{xd(R_^+5krUw^W=+ds%C@OicI29?SzBFJ z9I_fxTin5qRWz~LK~jK{9>4wW%ypNlULkUl$nlAcm%QI?W~<&J{;F8qE!b)Pjq7ynT7Z7x%dRzR#SO_~-_PE&*;413 z9tHVG_H{WV_KstB#aVf?q}QJBx>86=llomTXU-58d3tX4`q>IUhr<&KV=EIC$`Ylc zPozcz%qJ=i|1_;x_94!k&(>J1!2yw+_jqkH{rrNNf96bZ zz#HVTwY5WQVN8JrdPSr?1?BIeOOh6bixm6{4hy{6RN-jru<&-U3BB037lEB1#<84q zhcYMF=f%^l*s5jn7`O(qqukh@J}}jcLT`T4;7vR4$28tQvM@gs^WAI=^2|qDec*V! z`iXVHQEhE-;e5t?fhU3(ynN@M4XrrSq24II)!vz!;P&TD8KtRrQ|87(*B%yCb#Lr* z8M^ni-EnEr!ic5yOH`&6-F z;O6VKrcoV(fT4E^)1iu>d6(=rFS=f6TSRhph<>?YL`73%nZ zJsLB~3w;>RyM^4?#r@px92;a8u$u-w+LxRk?8W53j4 zmg^~HnJVgadZWL)ynaZ{e7JOY=0O_KYgIU|=C-%h_ZntPYsIZEJ9aGm$?v>thOR9e&^_0_c`-eQPNdrZjd0(yIrzN?mQ?ls3$k@@*U5(p{4N|xkPO#bxZsu{hWcU zxyw?B0Y!}E=L*wvKhF`r{6qvbBy(Taj2MD&yLMEhM~n6kTvP8!HGdH8ox*lidY_8+ z!44H1PK(8chyB68Yx;hO$It`EBN{1)LniWIgi7=WYgDEkaE15PlQdbpv&@r^5oPG` z2&CTbq}yg)*`^UiNcZfeNY$I)>2tW_iuh7GOmd&P6>F^<&)?-l{ld`wE)DOyyII5 zLgI$}r@r>R=iSmjV>E3i%kx69yUS$yN-RoI$24upJBYt-P|N=~sl`rR9;s1}6(ZJe zWo`u=_qWWwdicz8PrloX=!{wRzR2A2=j|uadLze{VqU(icIbTBa!yGuaG)rw4AUxI zoKtEcYVMIZsOitOT5-BV5vJJ?o_Zp`eo6dUN@Pk^BWteS^-2NDb>gI+^FNr>xID--g`f zed}jD%s>3*(1;{XI=>0+zOpqH|m&1sV#pA){_jnq1l6@+8FyAM$`Thcwb!nJQ`{m$fWM~@<+ z(HUMn@;cQuTa92wzR%aMTv-!*?)9d;@7hpgH?I0XZh(GfPr(-*z46%h*0-3uh!q25MY-b z@&8_9_=;OZLnGKY9p0LfkcfY6kRc>*vg1v1J2&qzYh>z&4}z@G?{b%snNL=<@jnId zXI7C#&)!J7*w+#Jo@e%H=r$i(i2bIaYiJ{I#>UhLA3u{xGG%JGRp95C{ZP^C9~B=r zCv%^^_wHFQ|FH{VwI`%8HMF^Z$QnG`Ao2QR?pWE!7=?iLXAE3#p9S0%d=T!6u-xB( z5HxHFf$g1bg2Y5jZfkzgiF@qFRb8#wk~3@Z=~rvrFxQ5I?OfcYGU@MIFMc!BL$ZBP z*&-7Gs@{bXM=(gENNgMoTg;@7(x~0%5z@tfF9DO7OHead1Xtg-&?YlwNqD^-Rfp(k zM{lnvBUjZ78=ywT&Gu?AobwQ6%~#HXb$fALtcSXS?nUc&`_|$k zFF?g$fgNxQXyO2`cT$@JJTJbYOWvR`gJyeBQMiA%7Vv3dpf0rc_o4z3*C!aMY7z5l z@rL%-!tA3~$huqUqlb`)l5zUKR@X_f6ZC_8hzJzJ-$Xllh71@&n5=;RtizLjvoWy3 z;IaZytQ-s^Ru1&*N(AE=nqoGpMJS8(i0f}KROQGIp+%5Z2~VC;_D;myWD`za(mbKA zN~g(9dnK}go>@yV!jK#^0>UT28@x6Qi34Ijfu+OL{UI0+FeDfQ;=Ge=F0_9bjU(3I0Zi(uKOG(XyyM;sK;)zK zgF5f#eBS=SJwT`K7Ey%H3HnlZU|>j^jmg9STsn{x{4+hT7O+HUEW`!dfHpTsas*9~20I)u9e8(u!CM0?%*PqM9Xc>LWH^Ra;C|Zuzlj-n?TtqiA6AL_ zNcOgmV~j0uME5Ksb-IWuoMR*v&M|2s5cmzIJhCrhqP}e~k-lv(=OrQPGckWTFxW^! z9FN9QUtl2>BSI{E=%`K`5ePyscK~k$_>UvfCKO+MULb3a57v4G{Xrv45!&>a)QjNm44%Tc9GA9_0@bvqe_($`X(+JcN zrh8yTbSkj@g#Vrs$a{!SThIkCq&Q7O7!U72TTgcI|a?R1A%=s0d#<~ z-+6aOdk;KB^u+SQ?LR>2o>)E$rccz^h9PBuBb@_r^AjQj^M z#XJATjAX187_uA8Q)EcF@RRF{G)E1`Q!Sk7ixol2FF7_+=6^2es}~*(P5WY5k!lP} zqtslb<%bP~p`m&|tU5AR>cIobXjMtt{#Xu9#N}ixHJWX2VM(^PSX9jst?Mi_x6Jv$ zLdyBUTF-)rFhC<#*pG1J#Lg;AbN_@Uu5DAWq1z zQ$~TrL`m#tCrRw*kY+%Hs&PN0SY5W=rvC|9xS6ETv z1n!FOXDlwEngb4%4WbW09Bz0E8DJ8MIsP^YWwQAjQ%(x(sqVQ#F?**Ha6IY@($T#U zU?B$?`J#XG-*Xj6C1CCf^FE*rEp_s9LdP!)aY23VcrjQKwC;}YMgGV_h)^uZ2~D2G zA40_6#Xxp0@0`UWLgVo?^IIk0NmdD*PXiH0IJjsuE*uw$3rDyb zjc6z({KqA-@(E3h5(88-BPgcif>=KZctT$v614xQlya!zZxvFK>fe}h{vo8|n;c^t^^bm65=gw#VJ zbaWX+i*W3RoEM!3+d2$bbNYgiaaI_FVNR1=EIj3Wq9*y8Vbiy8Vbdt^lqP zK-zVN=R0O6C`OoRa@5(BjtqyQI7;X+s3c*J*ts)d@K{&N-%(Al>; z0v}1>F@THP?ugBuM_d-)Nv&eI{xu>%2U9WVV!j;DWu?TFBJPasEv zcI0AZ!<|P8PQFW=<&$Gk+J@t9q0@C#SB-a-H95^4y)4*sL296xc&^oPenNV|- zVjI6QV#X5-33c+fGwwe~X71nOAEneJoKP!4kN}0W^QSPTa8bX)Sg1+eu!H|UNr0kZ z$^*dEPWb@+8&hTuqB{Et>H{Y__?6*Kt&m_Re-HAaheij5Y=VICux$?kVtv7nLed4> z@Kb69v8dxO_({iK2r!`G=cY;8uLO)Shz3_YMYJ0^g?7w0V9O$m_Tj3`N=?J>$iMp+mRO0KOG0~)0`=HD8nIHa0#&{8YRs_T|K$} zkm{|`kf9&e1@hrzL!+=JMm{n~(wY(Ql7WxxiXssNNhE@h0}igF zOyb=q^neL&N`gi&3teOd!`MZB?5We~!Qeb#B@PtUWXATD7vd@Thbuv@qiJ;N&A2Q1 zF8LRjIe-=B9Nr61edB~1o)r=&3P~Y8DF{(LQixK|&n7MJr+TEdf>0U!ayZ1QC=~Xm z3K1fU{H>my|H8s#+$(oQlRUCc9GIrawC(FldJPx#>7ZN_B~R@>hpEqp>Xn61$dAfX zaujzoqJF;d9}_}x!ZmH61%Dw<>RC&Yu&q!h4#rBF_ZNaOQjVW~4G|LfbNm$88{|Fk z>LBj{4DuepAn$?pAnyST(jLGd?a?ri4&*)1A;^0GgS-bi0C^AaLEZxxJC$0{KuG`yuPkNg7-Ww#4){V&*U(}w^6