From 1c18cc35f4497a977c0f87fa8d6f484812105ed7 Mon Sep 17 00:00:00 2001 From: "Edgar P. Burkhart" Date: Tue, 5 Apr 2022 10:15:48 +0200 Subject: [PATCH] Reflection test + orbitals (report) --- report/chapters/swash.tex | 27 +++++++++++++++++++++++++++ report/fig/orbitals.pdf | Bin 0 -> 9913 bytes report/fig/r_test.pdf | Bin 0 -> 21530 bytes 3 files changed, 27 insertions(+) create mode 100644 report/fig/orbitals.pdf create mode 100644 report/fig/r_test.pdf diff --git a/report/chapters/swash.tex b/report/chapters/swash.tex index 0aa7a04..79dbec6 100644 --- a/report/chapters/swash.tex +++ b/report/chapters/swash.tex @@ -124,5 +124,32 @@ Model runs with 4 layers. \autoref{fig:rests4lay}. \caption{Results with real timeseries, 4 layers.}\label{fig:rests4lay} \end{figure} +\section{Reflection coefficient verification} + +A small Python script has been written to generate theoretical wave height and homogenized velocity for a combination of +incident and reflected waves. + +The method was taken from \parencite{huntley1999use}. Incident waves are modelled by white noise, reflected waves are +incident waves shifted and multiplied by the reflection coefficient. Water level is the sum of incident waves and +reflected waves, velocity is the difference of reflected waves and incident waves. Additionnal noise is added to the +water level and velocity. + +Results are displayed in \autoref{fig:r_test}. + +\begin{figure} + \centering + \includegraphics{r_test.pdf} + \caption{Reflection coefficient testing (puv method from \cite{huntley1999use}).}\label{fig:r_test} +\end{figure} + +\section{Plotting orbitals from buoy measurements} +\autoref{fig:orbitals}. Orbital for the large wave have been plotted in the average motion plane of the buoy. + +\begin{figure} + \centering + \includegraphics{orbitals.pdf} + \caption{2Dv buoy trajectory for wave event of 20170228.}\label{fig:orbitals} +\end{figure} + %\subsection{2D Model} %Working on 2D model which might work with overtopping. diff --git a/report/fig/orbitals.pdf b/report/fig/orbitals.pdf new file mode 100644 index 0000000000000000000000000000000000000000..6c81c8e22ed94e59fa2278b4008eb0a2fcadd3d8 GIT binary patch literal 9913 zcmb_?2{={V*SJU-ijtCqYfj?M+$&_Bm3f}IxaQ$X6e^j@n3>2tlQLzd2+5QwbEae{ zWp2R#9MbE3d%xfR_k7QHJCh z3Lzl~RN2H5=K?_hJ*p6>y1OF|0+lc!nAkd4LPSI;aQ5aPBIdvVsT%|tHG&BN2SM!( zl6J5sK#+Ua$|eK?&e9_YaC-7DeD z@h0L9Hvk?$3=}_%4+BA>(0u$51AsvjX8?I$d{vx_gR8R{pbt>s-}(SmbU+y(;eG8% z<81+XAy8>sKp~PiGY4}(e~LJJOM(>yzNb`xZ5M(w&cu%5#-q`MwvhB@C*O@V+MTz$ zMVw198CFAwOVTh!bbfTiO2JaAUCCr|PnMP>%{4Wy^<`(V1aOjmQDsa#O;0~?DtyqP zc5~f(XSKb&ZfALA>T#V&O~v+VduY~;Tg`P#A&$WRD~w+^A9{}`5?`sg5_P7-=i0Re zHNLfxO>TXWt&O?W=AMM{ z(eO$5#X5Qxr9|VBm3&-%S*MFa>HYSPv-;xarp`J-?NC!995eYTOpHC5B6&LIkt>&q z%dW3Fb$B+KKYNMo7kSELUY7SVgrZsBB#Nq5OEtK=b>l5>#jH@44ws-!=4ABIyIEnK zEN?`rTIdDv%W>?bRFNIlv{!g^(aV@-p}`4J*KFZou~NZ-6E}r*g{N3E95{fXn@%Nh8RN3lAClZrIwGGodb{uO+JQR;Ywzv)9-UoGaQ#;fiYmxVc}{Zrma)x*EOey%x{pD(B@GPUG1 z#4^EX-8ge=?Y%84o2*al8zFSbbTxBJIr|fgab79?<42o`fgOhQSMubl{9Zi`4)tnL zvnh#(9J|hXo~tak1>dO5nXox8`NlW(lrNI?s#aQTy2`XpxN>gBkXShTF&!Pf)CHOa z)^YOZ=1Yt`t9Se4%;>A04;_6*EyT+lrQ9wog`;81>XH#QCf|v@A8vx6eJbPYZ6~~G z9q70HL1kuEic<1b=gGTI#F?yW;ud8eG7KUfxKS-HW;HMKc5{i`=T$M{z(4dWcCKs@ zHp#k#NI0Tkh<6U@Kgp4by(QUr{%&<$Dq$JDeexN($VnpKjI;ni>IT(hoOqKT}x?P`izCrn$*;-1T)#TmiQP5uv0u z^-8+z<f`PMsiWR&c<$mpzD)3>*R;~IH?^6SQh#E`qBGquk3{UgnM0^)1B+I_P)(buSX(IDLCGt8~@uo>ae^;q)T9!|VR(cCdoW=K7Q`^Li$No=We zC>EN1TbKhYLX7JZZ}|+~tI&9QZ0G@1CCdv*C_6hzU#gGyWesHl;+A8^2O? zGo(tJ?y-DYTp6Qu7ba@Onlit}+}5V=q&g8rk^QRcsoo_z0|@H9k~&ex8$B-?=1;G} zvsIDFBhRW?RWR{SsNWU7aV^TOoFLpy(~r+7EOUxSp;A$&=qSA}%!nM4zo_oT-oE^m zU-bJ@`wE4-NMMH>S)X#W%yOK*=|cB+#m|HHIKA?OpOqdd&TOO4wjTXxt=jXM2|uuW zf?53Wl7D4~bdpBxK0aihcvf6m-D#4^Y1Z%Rv>I*F zlcB>^Haeg5BDrt$*zjLF&yzx{dMEKg1_%1WrFZ_aA-piVy z&<)+J0y^3S{+%xT1rnU!DVV{#Bu5vq(Oqw4mlL(Bu3zwsT$%IwR2n>BV?dDDv3<}- z+V6q&c!9WEYJ@KIXl;x$!yT`?m?w+%GP#k?-w?#c?Vql#f2Kr4khjv)%6i(mhGwi9 zoKXL=rruN0PmQ*5ef|uMT@5)R_&I96#ICij{c^njq>!iLBTr)7kD*IsmMOwz-I6MN zTt}WV>C?QSyXAG`p zHTLxAT)ljq{a9Cez0^vD!&pE5B56lm7-aZ*mHs74{CdXIAr9KlGb5HnH zZ){l%&sw!4$MKIZ%5Gj-LdU1W6!ZRt>p3wsuKGCo0`dFc<%-M`z=5~AOaN9v~LpDBhE-el6 z9&^vBeL_#SSogSq81}{2``ipEP>%}+913blvTjowo?gfMsfW^8Wao_t^*OENQ3^*+ASS0> z3un%zm1soa6JyTwhQD=B_f5B0fA~$brC>WJC?ikR`;qB6M=^9y3{2JMmqy%=v!O3SV4H&YY8 z?3v45vAc<$%ei1VFVC{^dF|R%O)~v`sncOECZ3!6*c+Yb?!aitQ=M~zTqCR>RvNnh zoX&+gg@@dt{jmkofUEovWz}d+k2Q(csO_1i_z_nI5-C?-cXzGUbG*aaN414CiA`+8 zj<(0%hZrq)s1IMAmqBmGtNSrNnI<&|dROay+>%+oVoH0#PgUytYq-RTj^QnG;xf6B zdM6_OD(@@`axMHD)j6kryP_}p(~l*zHl`oUeLn4qQ>+`eo6ht<8AEBZmdrcsu^1z~ zI^%QCHv}R=f-9A(i|Oh3()YDos7khZX?1K8GV|C=kU(;Fy5>c@5fgK@c$^1C3d5VW zgm8wpGaM~*HmBy6^QT5Xc>g49=7qRHi{2i8#Q1qdO~^}dtKR$NSDo9TAGYF;a0ZXJ zgj5cw5pJz)uYMbzYOkADY5uw%%60Q_io%oE0^{2Qb5}yVQ**aIY<5n4TDrQkyuM1Z zeAwyu0QGn^olht55nXlV=de-j_REbcUhCV}#&>3#rylHJYXv>mR&oX%^^c#fGKzFj zxEYmdWV-I8!LyLTbcE^MrBB|HeG{n@3|rfe=5z1*sdwLSdBHHW(7Q%lT$jpcpqcYR z&0hV3k>&u>|0nizrhug9+biRbzWS{T`qtDhKDlu0!Vo%`4>tr^m-3C;;qo1=cFXg) z-jf?#;xP8gp!IR=Hj_)b@}nOAX@5LSBe+*NF{It8;mdJ-1+Hc89LJx!suRIUe&<3& zUvBBIXf5ZVz1zM@uXJK+lif!pHP2lc(wXn)vZ~#pKtHhzO4F!;xo3W7_@3vw6EqvV z{CMS-*bfgnn%;LCDWcyFuf-UpeP7ExWM%6h<3YMZ@6a2S;Q^aU&*fjhE&M9!aa-f=6sz2#@vscYuAA(U4bb}H-ITA zdP(JTE79ul+fpylN)~M-)z>AwYiJzY$1Z-Q?~IAe&N|K9Y^%jlR_9N{F|<8=f80A; z%aA3x*3ZPRUJAbI9Dc+v(ncV z-2GPd{+J0tNogxZax62_Z}Co%QN(Nk!9Z#%zH+`)Vdr_&*$ynKN}}8I#z?K1R>$Ip zhy{X5s#ib-+Pf33R<+FhRB<#4C7fT~5=UstpGc@|?~kEbY|dBC{yOC05Vvt&{ajbcsis1QAB4+; z2`>wqV-*JFjC#yyTuP}-Zz}6W3KnvIO(WkBY`?iGU|LkY!&bZzy3-O=`uMwyHS;Bd zOClW`c|ALa1fsa0V5;(0p6yrCayS15ChAaFw7`BQb2m*0cIpA?zkNh7zxdx|^pIFI zm8*gN=l0&c*6L_$!#_>#&L$1p^u<-_~-& zLOr3(TYP%DaqaSS_?}`IWtJ%si5Dfa+jF~9HpQcTkA_|Sl0x~y z6T3YHz1=WSbo}{*FhZ<|{R6l7D|;3Ccp9;jCbQwEv(uFy>%16O`Mj`!AhmV~%Ahnm zze-+0DsCs+NIWk}B)J%=UJw>y$$qNt^r)rjdE?pcXNL;u)g_~EE;a|`bYurrOE!vJ zQ~G`(T9H@tmAx9h(a*@B-dhCHi#e#}$mK4L;x0<`&Rv*?wCVKSBHc<-i$c zWD|m>!&m`poDwr7^L|*z^<*Jq#q$or9poI+k|1=cGqfWvQ*I-C9yA6 zRE5${8ygWKiNV|+sh5h*#aBxs8%$1|?B96VLUL})x%vbu%Z9Z?n>C@-R5JL!t&4EY z;eN%XCT#!9Q!i|n68lr0A(0nZ?haM@SucoQFRop2EPH$f`7zb*UyOe+W&NMBv{=;N zOez)JaIzjkc6di!jGTgE=o)8)|N2hw3B)bqS*p+XNGZ$QrI|^iI1Nv0AC)&7Q4STY znQ-Q99U-N5e7V-KlAHUW7FZ>xYECP{jHkGvI&T$fdpU+PUQX*VZ*V*n9}E2y|3yM)^C(jW5~4gS`l`wl zp-OU>{hbTmi7e2WGtfJ|I%%A%AU$POSSPM+JVGI(i&5$7$3so8o8*!WLo2Htvmkw4Hbm6UUqDTCupBduUjGjX*=e4oZg;}=Z++(BBD#o0`@ z%X+7x(isb;=m%kAZ7?;{l!+qJ-Vdg@yPs#K&`(|jd{T)pcQ%15twRWOC4B z$R(}nneaTTm*YtI-EmTXFP;HA?P*bUh@6=+I_3ITk5#63PO;=~BKnJtS>#jy4B{U7 z*fasBw=5sP^`2NXcp@Ap)|#*1$Z?8JvCo9flIltvg0I9ySO79BtsP?$X0Le;NfUzd zm1XhGO-S!JEJhPlC>_l2O#34P8cOif2+oEty71Z=bt++Ox)!zbm!~6{44nA5=%D%M zSJ(JS=v46e%c8RNTlUmzF3%#yZz#-moavdm+1|SRLfyD2VZdILzi)E9I$-KY(39EO zZ$C+ljP(NkMTcN*=D%@f1pMzIAFYC=B4IhUxN&-0lsyP~<67#4_U*d%pl(g327SMh zNS3%|)rf4C88s1)IKD@fBLi=&QM1t{ye_w&msLXx9@Yn11o$-IW1k<27dj*{(4DO! zLC4h9o1u20XymZS{c?NiQM1xJF=Z^Xfe_OR^q0iibEHn4vFqeYN7orXU>J+}G}@w; z@hbZ4og;@dG3N(94t$9czQ);2IF0DdL${qJKPAN~^erT_@6>`W>Meq%_tnXKYm(~A zTd=D{M5mnEh|M3%M;1>vIy_x1KBDz~`_N+<_2@tK14#Xse#A{&aJxl5sJN1filiow zB+lAI(^bvH-UTQYTA4V5Wj>$+_p53rWe=3Wz$zYCUX;XJSm1yf8(37+gP<@Fs0Ys3 z0fL4>paeIdML?hyK+O$;MnRx1_!|&3KLm=yTUrso>fByA@;5yor~pyW*@Xbk3GfF6 znfm`lO#np_=VIoJcO*DCL*TnLOMnc}tz<%Q24)26f_yL-xTb)=-&a7uRW$iE69SEf zK(CtE0j0e^6mTH!g4`na1Mi-`wD9HxD?l@Q)lq=|{WlJfjeYkMb1?&JkXQs3EWJ9) z?uieG2l5R*;YA<;tpFuPOBVA@={su^dp|#hZ!QTiW74Eh?GZ*yR&o69t$< zjexV$1pD{24+?TXfl%$={2dSnln)d~55f;oDJ2g=Z(0x&25 z076&{L_h$zM^|{BQ&y;NH9_U|uxvB7($%<4_PZ&|&~(f$c6JKs^3kC4enB|E|$M1p+((ilqSO zL<3`Xm9ozq*zS57I1>uIfq^Ol*#<=eTL5= zTmhW#yTQJ>_FKSY6#Ie-dUBu1{{1lk#KMY_5gq-1-v=9hXOVR^#W{AdjWdy0MCK32UQ~h z96TZ59qc8+^EwwK1raa=0)~R2;V?K7f#885xL`0YZVKow92#zj-P^jJSM|=*>mcSrj+c$SL1CGBO9OYA`^Sq(AzQfXn{_4~+zt=s)^k2-qL=&`AD2co7f~__GfUgp+^J6F`7V z1^oS`7Yqh?YxqYW8U_1<7c>Ugdj5gO4?_X21pfY_$BzKtRsPY(5Bq~Zek>4T{(*;u zV-L(ra5ll);+%KSn912&0P7j_OJHrQIyivtynC&N_rL+pu=kZ>;!N0e4=fT5oK)kG Jx+YEW{{Y8a)AaxV literal 0 HcmV?d00001 diff --git a/report/fig/r_test.pdf b/report/fig/r_test.pdf new file mode 100644 index 0000000000000000000000000000000000000000..7dc67f3c2583a811a03754a6ccb75cf2b2cab577 GIT binary patch literal 21530 zcma%jV|Zmvw{53mbZqX}9b?Bv2OV^5+w9o3ZQEAIM#r{oCpZ0l-?{g^=e*B--u+|O z+Uu!RHRdRdT6@+cl@Su5W%$YrLt3^1C~bjZ0MG-hbWLHnxBzrAIz|R|0LG6yasav? zKdlV_bb>neI_6eJ03IF~155os8QK5a1H#VsqVo1S_67jve?KH*WoZv!`uDYrj=jBs zttEiv&rcXS8Cxqoc>{ZZ+DAP?A|K@&INJm0#4SGR6!=dn_@7i9paGzhH`dWNu`~j( z{`+y|@4BW2diDUee--}!;4u7m1Nbd1t?cap4D|m-AozjC!1Ch;@ZZ0s4D?NO1gxAt zF#O5E%u4^29l*rQ@|6{!@c~fB_9OFunJ;HxXXRk4_rb>>ga0ode<=EIWc(@kzt|Ho zG5@d^Kqq4UL5Prno|XOwf0724M)t-4hJTa#0o%^r))T_Vy#+{Sk88+V9i(hx1I!$2*p{=aEFe&JI@- z;#G{dL$E<$(jP>f;6vX?V2I`s*GI4 zd?q}(#2xE=zBYa%opUaaViPdIz6W-2CEOW@9_T%YNRZZMxo2|EJ}UBeJ`hSi-B~|k zz)2W~T)bsIEBLaOvAt#Qq`zYakW$<2yB0qrA`s*Ly07Y12uLQ3yJ*3i6Dk^itt5O= zO>X`rE~=VHh#zyV6D$*wOq6ihLTl^mk$ct|pqy)J{&HSIin%S;5y?;?rynNMkJG2>U2aX%|)`_gPD!WBto<8pcfYq;iJ+8H{C8BAq;W=w! z1-*4w=I-4?;CbAc$h*9cZoGmJq)mcNkG4oLOW=f0;TR_+W|(q2q1--akaZJN%Jx-% zf~Mc-pcA87Mp{^G<@TQsP|z1Rq>?>4SVcX3{(84%nJZT_>~824gdOO-Tqt=9THe z6mYv3#b+muA3xvzW{oP3OV`WYFyK%*NQMq`-d7%yVN-~l2 zZo=&TSBXd|`^McjxM$_DjgnP`0KMT}*gMDC%HQT(6R@j-hFZj_u}03vmNjI$t=LF! zFoD%`{I$|9L=Wd{8tn<6=E%&9<(|!tTa!WPH4X1?8fq!wA~hIY1z}Tc!SOj0*3w*; zO2P4$cR0%R`(k1y?N-j*N!lfRGb=cZ>YX;LCF=R;4(FU;)E$dxH|AQ8G;Z)6Hdv-v zw{;@CGF-yo&RLVo>ZZl^K7qq?wj&=n$brEx&&yh-;McDm|9H%w&mNd`K($~s6Q*3g~~K_zmf2C6i{$? zB(xvR$shNldMKLzYzC=)tUnaAl2WPQ$;i9Mkut13D1s!g?;0W_jh%v)(PB4|Jg#8xY*n zUpA_CT_!{dLhI-yS^YuEt7K+qh6seA0_B;SO6Um*P z$Pv|BmL8Ev*|@Non2^;EeB?;oFX$d7g9?LCpjR6Mq&ga_)*mITckcy}FbxVZ9cwC< z;Q*Th$g~s+1?1!-7IIZ1=49#?G3GiJ3Bak*6ovHOIz^C9;~y3!1&81&q;HN3 zJ0H;-l<_^ZJSdFbx?9Orh(D>}%mDAj+ox`sqPyJJ$2@+87{+IQa{oGlCJ0;a;E1Ig z@S<&mLr*{hg*2!=KImf@!y=*2|Dt8dI?#0uAPuJPFT}g7Vi@42VdlB0(p%_B%K<~p z8>yp3c0}E|_-t7W%mI`|@Y|mS`8+6n;jb4#?Il;P#gmCDNbT~)Qb9bo4N(Avy=wi^ z-X@H$pV&`stIv<>y*a%#;rqfd&}50X!|!pjUyGftQ~M~~1K?AGQ3`Jm1e zl~%&jkQ)pDCzd|Yizc{fe{6&yY#!`5<|-*gmbz^JsRF|XL7%-Wz$~sf;;ZvmeTLX%OyrR6 zGwAciMDX=5CH?nx8t3wyo1viB^jri|%<3L%UF3XAF8d{2i4x!_i&S0eZA?>9kN(QO zFu_Qa|0X=KhnE-}ICrDzwNs#21UzPqlNggg97ESklpzhMhTbS1MNCexPO}y&w5%F; za(Wt3c5fLV?D3P32jcZQUQCETgR50)6bE37dt6^h4CAVA*UNblp(db4RT877E<6K$ z+CR#BEe14g0pXmCJwce>i1*o4G!E(n5{qoTc)*@8q$kdHz?_f(v=eGB?e#(#B%E4& z_fTMvikW3Q!=A9-*K+tLkXWvDck_u}SLfL)iRkrGBkCsPW&ol$&!P7Rz<|~)itE8> zJ>9g@*V9~4Sw)WXy2OP!=KM^kpOmwQ z73&KgS*gf;DC3BjGc&~NjdbXYdGD$(#6dx?@laKXW1H zW~gy*#2OF?v^RD{^yT5n{4`KPUqsW&Z|sNy0gtmw+e89$sKiU?Bvgk3Eo7)oph&9i z0i|0vqPO%Al$-}QM%mq>F@m0T@}|;B5uR~`xN)CcZ**_y-Me)nGgnWz8mZF<_Dzp{ zq#^m{%!%L@TQerM?fj`HKwwM^#szgFCfsoLQoUg)edU96(_jv3T}Nj=qovS5hk*H? zps9}>gTRzqaC2qW5%^%V3m|T16*q19#xM@k|tK*Dy33GmK?&GvF zFe=x7Ye(o{2MUyk;w!0SIxdi2moeY3m>jbFKC27MIEEWgFRPR#UaF)p$-_Ly&z&voN)?gWrM4V-g8?y5dauO7sV~o2d&iEGQbLu0^;~zW#Uj>CkuKB3j z*FS|Lg`-$ZxP*2YE4_&&x&VupwIoBB;|I1n2__xsT%N5&heDaSfN5BWo{C2yE`ya- z1u7#@%3XKh;f(D6<>nmSE|^(ppGtLOXiq++D|SXtIvf~CZh~jK3dbcbg-D{dA?~V> zeB~Pv9-yY`&Uvsvv#MHmbbKHt>^K5~tAgp$fmR?ygwtq1drArPEB&3 zWDzsf_4iRlEJ-T7g_0gobR}{$k$4SSmL--qwhI0{pY4XoAlWtl^_Iw!#MNz2Pk5A0 z1`pUYe1tvLg3ExgfhxRa3lV#YcVj>uXoUGIOCxsoujuupV&{I@jln(0ptz=z(R~g9 zPZ9wpnkptpqCC!eY`cK5tDj5P*3MA@*hObBT*m3y-^)8V5SHea+%RQ@#7_I!8IcIX z+nkVtfJmyI%R2Q7Eo`xG-cXc+eu_HMe&!*4eWQGSG;r-#!BFLi+TB=*0jC;TjShf3 z%iM?0HJ|fdQ`KovaYgdfmL1zekyKnIp+wX%_eAe6>UBwx%U!3)KePTObvf;TDH9?} zHB2<{KF-R+sx_e=;TLDUy^B{T7r#sYHi8 zAi&y~`MsL}$fu{rrmS%5WZY|~-tX!oVs@C~X9XW#hAEI};+7TthOin}j>4jPFeqhy zI8$`-1Qui-tQp_R<7IqVB1z6W1?co@h$55rT`2?Bo7#3wbW5+cYHkeB0XW39AbE%B zi;PXCA>}$kG)dfS2AKvjboq_AjwKBqZ7Tth5$y4BJYk~l4DE=bz@5d122&Wm0s6Ma zLP&7?G2~3sS3xU#A^VNF&+W-X3fqg>Y^oI}`B!OcG$pz6z3Kp*oJYj~L903Tr0!+5 zxclu=s6y3K2bXQ5z{bQLoK~_%v4J+N?J6j80+2O8YFb%eTjtN{lMXo=7d;8MQFgtx z$a7m-7sxylE*v?>qsl>^p!EYs9CwCTtGK-yi=M@MiFBl>8b+2d4I4bhF%x|E=$gCH z{J39{@r18oCwd^3^yf>!3a(F<)*(s(yA*gE5U%P|At}>KS+I`+0n10iC*zHtN2{AYaTosMK-L!)kUhxj+VgU(a<=iLwXi!yS^AAtIP^_6U0-dY=R+~>=mg~-*6rkZR(@5xSb~=v) zA8AwlyZd42eX}`aEw|he-lmL7H_+-4r8n8U480MyE?e4X#iUqZ!Dy z^-7gOP=S5A>*aY-o_V+gEBU$DnZ9hLEdJc%!}f5hocJh)55g zY{8~dID7@_=ErfqTQfVUjUTXW{Q>#nx(yLuSSQs)vzSVS=2FJ9nXsw3Gf>j2p8x`+ z5ozNKRqEX<-s7VwJxU(!SP&+((qBzT!yH-mF7m1vgW=5u58}2`Du88VZ-Mg%l1$W_ z*3RLtPU(GSkgUTHHspvnQooT%Dk4PQstKsZvrkUuz7FK$FB<0^4qizb$8 z^ehd%sMOa)9q^?A7ZzB_IS>P}Z1;rB1H5~Dre}6`_;A&7o*}wHZ6?9z%9!YJ2|!cS zq|)chBF$-MIVzG@OHi`O9}f=o9s?MS*v-24aE8-1DFo|&ch`|WF5E*J9eE#@;l(&Y z;a1dsGj|w|-K}+2`1y9=!yVd^q$OIT_e|~S+B8n3xHBJk{IGWy9yab^f4kJ+Jlkx! z;Fj+yV3nQ?^NWM>ayOX+eSR4K@iV7g?IF5jjh;rbDrz;Qk`Xt1{o?1->I!51sk-Tm zTMw{jWaG#$%hPjq(yla>r;l-^dXd7SY}8z%1)IoKA%DBpjp1o&*QKK+e-z7V&q&!nRst zKeQ>&q1c5kCYF!2TKX8^K-3R?zc9+T104nT41K)4iYu*cDh$^f_nB9i_=ZwLc$Fs2>iJLslW~RAS1e*0=^$~4?5(5fi^Nmi5#VdkHzM(!plTRc-%t>us?bwG$h3*MgO8yuKMie#Brx8t zZ9e!bhz3ja(tL5Wk3nVoU07lzVQ;SAr;SWXgNK7L&{@c}UUPFAfdjD!ze>5z&v; z?Y5V`EFTOXu3kO=w0zqmfhfPeL zFJVMPP>dU!O{5l|IJUO615H+*aCna7E5Eqo9e1o;L`2M9TGBVLxV&LIO4$9G!9v5# zrTG0BF>?OEIrY`z=*paF=rOLT0kI$ksRFy0Iz#&n`=}i~>jpHH7)JXTG@Y1Cck?Gy z$r>B4FO=PRizs#x$l5WLy2F%lj!DeRKrP%P=6iyr1RL#`YLs;60Rw$lC*fl%NU}wq zStbGa*tz4ul}2j9K_w)DZ=-K2ad;WXc4_)xy}m$5d!gCM%H2>%JQWQgOoLTiyu$%Z z4F|(tQ?k+OMse^QH8u<~GeHroP|T3RyU=5hd#X0xNs;Ay;*wf`2C-oAO(?D*fLGLQ zUHoLr_un<+Am4bD0+);%?89ntHoo7d}Hg3wkuM-BFvW?>TI@1U=}S}(|ww(W6)9sM+72b%~hk?8?rg1 zzbjUhP5 z#OOz{4g&6((+{12IsbU8ajFw>|7CQk-%a=~TYp3rNj75N5Mj~h zNyaHGaK1S)W{|prgfy%kqt#IbT8~nELr`87b0nu}c9UaUmz0nsk=ma!+TWO;1{@De z1s_^6!MCdhE)J|3l?5Ju9$o}z6D6_JGEi5lumBB%zYNWvqRfsSq8jzoZvj^^znli? z9aQ(Uc`;VH8bZhHoxl9qdRIKsQ1=;JM4Sdz2OP~}qYs3xQXG$Lle zuAO+ZXj$vG5Q!Ah7G{|DOjc|i1kIrFHKC3t<1T4>;hr&Ob!!`~V^l0yQ_%f{b8O3_ z>&@gSN{-I(Fw)0p{PrhbSDj+$<(W;tG^Zd3X3WnxadU46soE<<@sbPGTgAyfJ0z%t z?@@P06zyShCDVTuVnGQjv3Tg9sR?QW z(dZSwe;LEWQyIh03h)R=0w*@uiW39M5i81reogDnHUJ+uB(L)?Q*?ds1-25Gat~;z zVAP~I!V=4~@9By_3YoBI43;bv+HEWbPy@#YJx2SYX0bniA&vJhBU@I((ZyAnC3jKa zkQy{SiZ6c5V=66~+J_odJ>h||XFrY7m-~_M+iikU5v0T%q?EtR9q^cR43r&eAy7qM zNRh2c0H&Fh_mxbiRLzxucgfqz<<=+b8~HN~s3tdwmaudd&&s!p)SnzwO&#fk<|=^L zN%EeWtH`-J+_v=`H9ANMnU@EJYAZnEEjakpJI|U0tH?Rvj&2Gj%$RGsjTC{D^{Za& zc6bT$gmvZZ82yeK3{rL<@y?pl9m3aSkS)*InPFfoRF!WsdmXmR+aoD|a%tyNDzQ%o zw;!+G+FLB5n5+-E&8Q2o$}*F#;SS(Dph{3&mD$XrUiX{uGq4R8R!Aly*I?o-=px3A zYry15;%ep=Q$o@Fr;Wyp!1G6wrEbKdztgyR$n;Hk5vM9s;( znzLjc9D{r)4Ge|9ovTGN=bEPZVpTxx>yzPpO~ORajrXa>Dz{f*I^Id5J#f9l=xqIujB z@=G?crBI)Y+FiY-Sj+b}cSuLy{!p85-fJ-*c*S%oyQ;sRL3B zTZBfedYg(9%fdr(*`S3{nJfzQ;CMtvNH)7FEoJR6@V##FkWs{lYP>1m#@G>Se{P_% zsPE`WMNgdZ25LRMtL_7*=nM|iv}~UIH>okU984k)8kY?CX+@uC%?tCYwUE~;dj0bKcRjE5oK>Ke9+lOi7`m@I z!CTbh=7B%V>kZM*9d&H#ao*}cF0>gG*0B2#(Bxa&>;=2Bb^Gy)vUWWRbzOdF)%){@ zP2tEFHc@vm>qyh9pX-UJK%vJ3ha94{9T1+tL}fQAT{y}z7_8P#!S^xpO<$S@>YAph z(jtXmMZjvywtfcdo1J?uJcnVS*u$V7`tdeVoP^!DTBm=mJbIDf_RXd$DRTo)PN%hS zHjF968egm;#ub))0L|1JwS`Vhd$ee1HJt< zbc~&oIXw&u^7(?8M!LW509dh;^h+kd1{i2x>ECwdijJFLLSdLFxu-qa^=abRJFHY= zIk$%(Fl_O5ZjP*0uphmlVS+L(?D8vH?r7Ax!q4Gx4vg9o^&kR}P~~!k(seizde(XA z7KvT2YGJJODo*v$@##WaHAtW0s?(ZQ4aCuYo_$FCoxXm{{6JqGMNLk?G2*CKV70jr z6zkQ+GYdSJ+Q^g6yP$;#LHXHT$cZl?sWE_2^+HZzpQTmI{=l**)^mC{Cnv+&sV?fFa|C7De%q8Kok>fI7 z_Q)rAd0Vmm9lOBT?K@nE5VSk+q~Vdm2}ml!(l%52wEOBSM0Tgq_9Z7#Coj|=a(hM< zc$A0Cuj)X3p`oMulDJvIZ)^Hyt?Hj;e~Y(d;sKkLTfnA`7%)bh_HS^j_H|9h8yw7* z3BM~wn*742uf5K5y_8_aV_MY5nI85^PH&t=mGb(&R&RG|?YcEJ6)UYi!5tNm5ogdT zf0kv`nA=?I$4Mg9u2Fx!r%l#suq-jNmn2^6<+$Y@ee@Z>FSL;`+QFZ@-~m6SMjW9~ zvO`I6c{8br3w@RA1WpETBXV;MUkX*CZuA~+O2Q4TziaZx`=os>+)lC_bjC+dbyON}#z&PkgG5k+-@PV}v*Oap>7u}xOkJnpJZM!X6i&I{V1o0| zxH)xmnk`26h@@(e=qdC5(IrB&%}s)F=JaQb72XuZ3-=>>dFViR< z7t8bS1>;b|w*oWMDNHCFPsQSB;I}g`dV^z*@{f#EgyIa33e=bpWS7`u>U$l7{Omj*OZF;DOCGtQg~Z#dsG|c)B^qfMTOc>AZWT=<<-Ns=V(AdTo=m zPA9Pp*rS&6wrD*b3@JsR9xQUMNb8r^mijzzit_{l+Vl$Qmpz=;hmVpKMex^>IX}z% z&hUlx#E_ATKjc(NCx0Z{v#{YTX&rIn zc%sCy$#H-o2*da5;|H9+rkwb)+t@mkUhpmXq5wbd7N;vfa!V@ie&WHFk4v%4;v%jb za%OANUpu@v_i&YDx7#F6vBE2}&4IwK)GN2_*#J(p6tgoyqc4 zf@yZMA)gvmbwf*Mp~9n%hqX}Razl7_A1zp{NnQK0LtSwG%oNj0k>5HBDIcC$oL=jHR-(} zM}rj|LZX-+y=VF|#lsa#ln%yDw)xgaxFfj}DmETqqb5!^cLSL_%^G#BG|l;HTlnsF z41JS^c^Y@c4wjM_C(;!I=;@POO|S|m6cgYQiy~!z&!&-X6{ai}#rbbCy0|AfInnI@p z=^0^`s~xDqXwGkF=X~SQoEoDGh_%7SI6|(8s_Xc3!!<)hHlz)lgN9)^ypH^@V7bX6it?Jon;c1q&7mkZom^M5xBk|B^Bht%^Sv}f;28-z1UIoetnY;NdEtg_6 zW=CkN#x}}x>K>1>7DkKlr=~wFlI;_A?bHvazA9DIxRP)%Q5RvA*_?}shQ1{Cka`uf(^eQH#q=A=q9m`HCOvtBH4PXHRP z`jYh4Gtf=ook?u5IRn$ifxtWI<*sZdID|K4u7}tjckwwsiwc|5sQN-}Z>;h8R5FTp z%5W;J%hGW?*F{gn@tL&@5reu@`S~`kBBN`J#g{iHQ5~EX?m+@n!!Z&033due-s2bd zo@M`TbcfW94Q*7M?f~Vws>Hb&U$r>v(3tqf7vz>$ z^?WL$)l+Y)r`NDTJP35Cl$ObTeFD(27s4FQS?%~I531=%#jv9;|6v4%UY(3XTq-!%jSO?^l!#B`RM z?X$j8cAr@Wg&C&_aMWaBSWCXMXtP;s4oWlkq-NiPl)6?k)qI!OjFSk5b=l~rN!#!K zfU$ye`#dxj`I$z^(;4$&5SAY)MJmH+y!~XR4Rklf?d#HGzjtucky+H1ZK#Ucp+v?C zpeJ~bE6C+)-9PJ1!OhSpjLS01*19JsK^!UM8kU0l1g~S8&Y16JKXLB_nY0e)`W0jp zso?3-Sp;LSq(bT<*Y}*}+^vRFWXDzM)Vg&~s6eNVwa-l=WU#ZQD%89S1?$xUhntzF zkQPfuSc0s3ymF>5P6xXU?YvK~ZDniQTgIhzWNQP-@zN1lQnTfz0%!{j>-(CM`p$7&LNJ{f@#m$H9w8XLX7TZyu9LCFwkqWVi}rIpR};VUOd#E zVX$bkzA9FDB)OeINPxd?Ua&|e|h~DnJjEB2kp8EI#ja@M? z)4FeSs480&z(!}1*n-bKjVj`2NdS#NYTX9DPs1dYx8pNWg47m67W_UD?lp&~>dRgB zO+rh4jUX-5V>0Hh#h9wj*{9iz0BJUwX4QiIzkH<7DeN&vI$Q#t}Fb9jKd^!Dux66#6dPi{5p9Tuz?@pZ(jlGsOupoq4k8tO>v z&bEG`%!J^dv#7Pl#>+ZLoILx9I6*(l4F-T`bbAW8YB0s)ZATHM4|HFE0}JIE<3DsY z9k{ieZ1&=N*wTMn@3FjoxLQVgvqQ1)8t9Kwdunlj?T(i64_6l1cEneHk|r(1hFT1N zkyhVV-2C<%^NEj2lWMnld}1iWlPDI^;gxGot@HcKO>MH`XyOaWO>3vk3wY;2GThCB zd`54ZR>uW7=i{^4OFnF`prv4!Ox;xdiP9SW*1C|*W@QzN#(m{|?oQ`pYi497Y#&?q zxhK!N^zOpNCYwk7HNyMF5FNSvwUA^xb7LkLuA~)@iAQiIOQV4^PKx@a#09v^b0j%w zP4}|tzSGQCR>6U^Ta{CZmpn&O9xEx9TaUQ5%7^A^PpPk+?d3mzirQ_r4aK5=Ptn>U zAG5%7KZNjTCm=mH^EeAh_r)-8EC1;(dffqRF0tb^`>2MQeC(19flezp=AMN&F;C%n zce;_q>fCtzMsWNYrH|k@c{>pJi+t?**72jB`)%Z?(RCO%wez=zwLUT@ukigUXS!g*A>Z6rEAf2^4gPi$BTknko3Qb*(x-7e#~9kL(4#1~KEUbJ;0I(%jo~>tlF~Q5@e*Hqs>; zUzwt=pV&j6+TuoFZ^G6-$0iaXOf~+#y=*i{y|D3kBi!ke%Qk*hsYHHgG<(2#Kb8$F zWv?`@8*4Pppc@p3d4w>60yBRFlX-xBySPd8Rk5Xi&Rb+-*$`HAzPsIfn_pk`cy85t zz1*+vlj`{{=aMO3bIsZ!BG{Ud*rgEe!#jOOLxKaS7;@go11L-ZfdX$A%sddC0qE&rF}_Mf|I9}e07mGh6s_doHmW@G!~nEikK zt69Q-ed(fy0)3bK1;y?!b5~r#^ZT`vR>o-;90VHR7uwriVZG`(RC119_~w-j-d=nTcE$-tTt!XahDOF^al1< zdrIrvB|kg2^qy3hl%~J|^Lnwe^j%59BYHA%^8@ig)^TIo?X;DWA8V3k?fqM+woa|X zW=jd(?SaEj?At$prBr6otVqeiW;nk@@K35Ha!ZXzRH#ggXZM(_N5U5nJPnl+qfWuUCMkP|8Sp zG%%4&B6*9#F~n7Fla^;ycz1>hTv~1&`*s6L1?2$aE(k&PCk3q0C>(fcOuRa5dJ$d1 z3=nCD_66LKWWeWc7(Y84gA+wD^57Dph$E&5$_y&kOlNeO(?M1>@)+(Nz^0sB@KAnc z+%T`ja2yPW;S?ARa`T2hiWK`utGHYkRvhCIet{#Adl}bYF2{ks)OZUuk z>qf05ed+iCd;HOdw7PL$AX>})b}%cXRQS@t$vuh*6`fowo@%;>19JQ-<~}BJqLC-o zy8ink( z2^Ec%Ol3a^B=|!{Nnd2ccDn=Lw&}cn%~ySC+T>>BuOiu_6zI2<#Y#cuN7C7kMk~sb z$yDiHmpwbaX9PF33d@Jp!@Ys70~fFm>-~Jh_etm^?nil4gb@i!C)%o!F0S_eR1HW4 z)(@cs*OP&u;_;%eHlbdg?@~_`36ivm{g(1rS}$>7Gj8_aRK?6!aaWTHwUe;Q!9S7j zfG=VX`*Y@m0UVZdSHZgkOnr>VEVH~ zITnA@rvLr$%KvW~VPpPFlh)EV$k5BsGtYCf2_| zl8WSq2W5Z)J^$&S$p_$bo={XkSbo=%Sa6(%#UJuAmxb^q4Rj8~*_>KNSQlnRGulqm zlB3imO!{t}NiUBZPU&=|7wKgaXG&9O?5rsx4p^U!E9?Wz}?+crO}n5&cDY zf)-X_-aO`!(HqhYGfC88+a!Nr_0GKV?&Z_CIkWiR81&~f(myfie@;A7vf?42^l+{F zXit2Enz3z#g2{sW@AJd({y;G3%9z-MBxP6x20sEC?uTeXJGNIGkDSF;6_=Qy#rLFr zB6-;dxuUPS0vbU$6WK@?+}UQxD%I_mK$_fSyADVK|0Xb|Y!E$namMDNy?8wm|Jq4~!+#-l0ieL% zfAB-Yz$|?yiw(Sc4@YEl)82&c4!40!J6E2 zY4S%}!nBh~w73H0N*?mhcX%VRsG^@!dg=GDde6jMe}nSRF|&VylI^bxG81Ft$55(b z20;!$Cxh8Dgn*$8Q0MC?4sI6siXjHqI}B0_=T%3V<`ZP(l?1Rr>;s?+RVN2y3E`!I zM2~EI1N}t0j6U=?c>bJQ`zLr9SegDZkc6W+kwHHcf~#~k@dGg^&x9c(CF{MF-h#PY z*R(Rgw-v=_X2~b`z$dVQbQS!(9Lz}X<&f2q1$hU|!Ami3+c-|ny}8dM-u}{<)s7mr z90ku1-X$&^(G$2Q_vizw*e+nbXoT^p-cU zQa0~FeWIXaq9QQ@-I(5L$OpY54_qqc<|TmPxceJh|EWY{V*D%lEV8W@ee`fAEZ`XQ z01wi!ZGc#Gqh~7TQnA={8-ao*W4`)VhaoBzovVC7hmZod48wvK_}i%MB>Y9f7WugZ zS4{MeP1ZNO>J*xLYgpEfPbtWDu$T8x7t~&|T|!xwm-h%QD#(1;s4u7oDWwU@x}U#< z`j;=FZ7Hr&DmRrrIFd9&LfbN^9LU29hfp?lxna>PMysYB+p564Te^S7*70XQsSOp{ z_}K+2N%E%eWj?{9c5-wi%<34xemLnB>XkZ10>UrPUdq#?^K9C>x1UTw&JLWb%wg#G zBE2Uex2DuQFp(z1yQ48;P+VG zm0z#yzEsLbl1E!wb!=&-6~luf;l`x25o;({e?>+zVd0mW=u)NiYHief!>J=~f#bUi z?zKsgi^%RYv(BHO+a;G@_#2A;sT5`Z9|JS2_0NVv$hANl%!BNXM>ei|t@8^asbv%S zDz`cwbTXwSQ;#4-1{_U*hK=`40zIs{s7yrqNv_Up#}QcgZGqhNIYL~@rmtVlR^X5=%ZEUS{J9+&C)l}%#N>FbuU&h!HqpjwI@ zuiKPnJV3}6kAZr>=>Lmn`&3+iFiRKjTiD%az>gjHMc?S!jPe4`E8t7q#BSb%O8H#% zZ7@v~KL}ow4VSF<*J82xbKJC{tY{Yc*lvPopUhw>DRU;}d~Z=73IFjjE{28LsZWq- zHb5J88ED%qY~SU(h|Gb43{ORS|5FRPkPj)rK z8}6<4i?NA&h9c`yzTB@a#~n@tfr6ew6-}nFmtCi56JtO9QIf+^@TM-Gz(NEd&^80N zTY)ff0mYXP&ZKubI|X&oPX{o&J7X%j!cR@%o6I2EO8O=~y94Op#&JH(ZxP)JFE5Zt zhn4+*)4cyw?y)fZWo$CB*8N-VLB@y#@FheC3xp72h<)lE-SvUq^HPD}!_HxctRM@d z_ND&vX|wK2(%)45k98d@{a@-9jshZo=sGaFb{}#vb)WmPk@{ax@5PXFqB)13nQ>dt zrG{Nljd>2MMY3Ba%?@b3#R&aEkHM2E)C=OI5R(+sixF$U_a#nIuGL6N{y|kMw^?V> z*WJ0OIy)#7ZXNy-dC^$?SRpHn0fYr>o1Z%Y{!3N=s(qaOiX{3E3|SFwtWD&Wjh#+ET#v@9j-(B z7U%Z&hL~7k5qRwzhTVYRZ*D4p#@xBy#_&--|INzMQU3|Cbn^~1H$b-tR)?^0Ufj*@Br!P&tjK6j zP%*YoMQO|Po9LJRF5RKTbW8|DgD%!pRr+Ounyd@GjI1RuE1%+m_4a&qcFF1S20~-! zX7@ME@=rx1+h4mM6BDcahdiVIj9gF<0JA5O@Gp7B6BIuJjxqvLJ4!R^O+ySDm4N$I z5LS=gLf1S`+fWxmUj@JbXkl#n{V6F01@v#w{8O>W^s(#y%S)w?y)^j8JUcmS86End3fd9 zSc~*)BlPmR*_&11GkdxQK4-hzES23sr?RmcJ=SLlV@^<0#&0<2!y9^NF$dGBnxyG_O&B@LfTA=_xdDdd6ghE)l_9lNq)7I{gFqs<+n&KMvgj|E3}Su}c1{ zTKU(E=+%QNSi$@WZx2-NZ%Y2hnw0r}HZjohLkJze!ngyxESJw0KQ!q}=a=F|R8Sjx zxj0VMAC4n|{y*wZ>>ummmyp?ZoDBSXp@&M5Pvru{irN$jXU~%vD;rsizy{^=yLa+- zY=(Szd-4lm2Ob-$D|ZL+BJq{ch_i*d82fq;uoK45stnaU5=^=xo&D0N=V%9%!R2sP z*-u06vwzaoKm@fXCxpP@DF&l^@Ly($hP&Tk{!j)?o#KjHhXu<-dpUY9iNOI~X&hw< zZkQ+-P!5)Ggje{zNn?_LeJ%Dk>kpIFS#JwrY@Hlbhia-SdcS#aPUHG%7!O^PZ|tw< z27FskA+crA zf2>mf$4W>mD{c$Q57jzBedA|XxUQE-9k}y84>b%4@g@DfJ3a?rQQfMwJ;Z!m*+ zqEX46%s*77-b;#tAvvo@P;>Uy_(MTz0h+hp%Dtl^lbv4EN@i6%<3(r0j}M ztn1Yh;nuDucHLQGe;j>!+jVT{s{6;sfxvl7KGxZ-^|qkx=djBC{l^zjEp4@+zv;AcFH#MY4hz_>(jF@)T7oU_2p> zmI$kQ_1=jMg!{n&T?#B3{;^_V6m*LT@;sK#4}HIuc)B_$pWmN==?wv*$c{UIFnoKBW7Ufu#`8Ol-Og~bYDQFw#E4}omPbfAcpY=1p#ntZg zT2K z$M#QjnvcLS{|PoDZ21wrZ(x_`83Fz7%b-CDw~&Tm4WVYZk27ngI#W4%ny5B!LO1OLu>p%$y5rqN7k^3MSD)-BXy-$RdZ~)M5rUB4!(;OaL$ajMxAB& z{9p=vx-r%{*4-h=<>t6OsqWiH7$VjXg8cnopz-79M$k7W2wIUKA`gBf%Fz+JF!8 zIq2gs5Th1>=qXHoV6q;~Ok1?W@UO2YmX`CtJ zgQkioHkZa0Aok#N5fYCj%1G*WxPgF?g~f{5L#ZG+htmyubM-Aumh^&%R9OKc{Jy35 z^n;10M2a(k6Gxr6P$;bv>FXGTJg6S%2+;@V2U8GT!!PI)BM^U}PxK(Z()k4aVGNW@ zp235@(F0eWYY!VjR{FpT5DyBzXfCmcFgR*H($$%s2}6j%X$!D^Z<^1m+j? zURgCv z?_)OYW{;^7aYFi-Jz+ZBCo!qID8Ud(`#{S2byx7l$aLtRj`zYTw}?GfxULoGF%cR>O;CbPmCO zfkSXw;IKKTI6UPrC9jh#P#nE`f-@%2|0g)Dm$HD4=9|!&vgH1q;TVACI26>F<3fqV zIioXSNoR2goYUD%!}+H3+JakuigP9`=%9jyEY2gT&BvBVpfg+s2VCF`lEIuVW4z;a z(!I&brC#97b3ZWNB4f?Q`eHe+=K*6zm+P?LmwOyO=5%?s(0!p9#dr-e7lxsAhV$7% zzM@?C8vMQ(&D%DsYPV^mhUdrqwnhq;SdiY-HA94i(%b4!lnkF literal 0 HcmV?d00001