From 319a5766cab5ae0be53704e3d792d2af8533091d Mon Sep 17 00:00:00 2001 From: caes Date: Sun, 14 Aug 2016 01:41:32 -0400 Subject: [PATCH] =?UTF-8?q?tophat=20transfer=20function=20basically=20done?= =?UTF-8?q?=20but=20for=20some=20reason=20proper=20output=20requires=20V(x?= =?UTF-8?q?=E2=81=BB=C2=B9)=20to=20be=20multiplied=20by=20=CE=B9=C2=B2?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- plot.sh | 11 +++++----- pres/img/tophat_freqdomain.pdf | Bin 32808 -> 32732 bytes pres/img/tophat_timedomain.pdf | Bin 15119 -> 15268 bytes scripts/tophat_fft.pl | 38 +++++++++++++++++---------------- 4 files changed, 25 insertions(+), 24 deletions(-) diff --git a/plot.sh b/plot.sh index 9b4b55c..688cdf5 100755 --- a/plot.sh +++ b/plot.sh @@ -10,13 +10,10 @@ echo_bands=$(ls analyses/*CM*|sed 's|[^≺]*≺_\(.\{5\}\).*|\1|') # echo Using list of reverberated bands: $echo_bands -echo Propagating tables. -scripts/propagate_tables.sh > /dev/null - case $1 in "PSD"|"psd"|"PSDs"|"PSDS"|"psds") - echo "Producing PSD atlas." gnuplot_file=psd_atlas.gp + scripts/propagate_tables.sh gnuplot_input=$(cat scripts/templates/${gnuplot_file}|perl -pe 's|\n|␤|g') for tabfile in analyses/tables/PSD_*${errtype}*.tab; do @@ -33,9 +30,9 @@ case $1 in ;; "lags"|"lag"|"delay"|"delays") - echo "Producing time delay atlas." gnuplot_file=timelag_atlas.gp gnuplot_input=$(cat scripts/templates/${gnuplot_file}|perl -pe 's|\n|␤|g') + scripts/propagate_tables.sh for tabfile in analyses/tables/timelag_*${errtype}*.tab; do ref_band_extracted=$(basename $tabfile|sed 's|timelag_\([^≺]*\)[_ ]≺[_ ][^≺_ ]*[_ ]{[^_ ]*}.tab|\1|') @@ -51,7 +48,9 @@ case $1 in ;; "tophat") - gnuplot_file="tophat_w_fft.gp" + scripts/tophat_fft.pl + gnuplot scripts/templates/tophat_freqdomain.gp + gnuplot scripts/templates/tophat_timedomain.gp ;; *) diff --git a/pres/img/tophat_freqdomain.pdf b/pres/img/tophat_freqdomain.pdf index dbff5dbbcd5ee90743b2c455fb98b3d2f4909815..8041846e2d905ee4cc1d084aee617d08f957a210 100644 GIT binary patch delta 18680 zcmV)KK)S!EfCAkA0gz39olB25&5_mj^(lHeEZkJyGq;F~gv4UNyU{F0U=YX##t5PL z_H<>Ob0V@HkGo%4*4r7+zbda46&V@vtAD)ikLQ2?@!x&?%YV2({@ee#?$7q||GuRE z{_Fq$U3{QghZzx=PiyMOoflaJs3*XwV7SMGoD`~PwMo8NnX|7u^-=Unq?(E9wky#zn$ zi!QC@v)t(?*P_;c-apINf4x=+eo(7l%38m_mbKoW^Xr08s|kKmt1eZo-(HJa+t)4D z3bjJ;gIfJk*82UmsI`7oR%ckF5d5H4U8*&DT>IlS+Vf`_ckcDqHJX3Z8m;v+-Ep-* z@RPpiQjXyF*P_-k_rQ6r5d5H4zm&Cpe=Tb*vBvVF+R9OX(_g&?wPdLw)Rcd$!TeIz z`2Dr4l?sDvGXGfDWY;xVlgG6`UX$zoO!+Ub$@H7nWO{3`^kWT{e$imxU&~tW&(UI~ zyZnW9@bSqP@Pq)@ohSKvwDcsqXR%QRVknqDJZoEY0;(4YORq z*up%n{K>k1%=L5R+*|*el<~Q~UZQ3FvN}J0%i^3b%{70jxaM1LeSalt3|edcL}|^a zv;O)OR>OE{tRI_PTJtAdf45Gy;vrpPD>UBrq7mb-!%^<~x1nKF3#Me_4UQakH&=5e4nD|7bJ9&d`$n=CUu9 z5}gZw`jFte3i(40I*Q*_*cLnhtWbjQD&!AQA&cKt*vH;}dxa8wS0R6h3R(QF!msvuRj6&9qMC$vHnIs)0fp!T5A1E%d~%X{okK=`pFZJ za{tUNbVh>jMvp&4Z&>_p^is?F9T#ep+MjfB;QubU$Rn`zpZrE zZ9jxcCHb~e{*sll{B5P{J)i5R^SRFXTtA!7eD~LUy5zem^-I*L^0$?KJv9BSRFZFh zE9Ea)Da+qhn%^teCEr$>&y}n4_m$4qE6OF`s2lxIb)x)R6;96g=gfCW-)Nj4DV#h? z?948GJ9^Kr*o)`yZ%1wYGHQ~1TiJRI8_VBSTGH4%lweIr@~t_+U(%jn`P)iMdRY`M z`L@!M4x2*dZ!7(>6wgzpUo3h1tJu z(4VFJ^`zWu^#0M)=c=jR{g?Io`~Uu@KYaY_Km73@{^{fI{>wl7)BpMVXEE=q_l?$< z;eXdY=g0cYqqe@EfWP_uAO7^e|JR@Y;g28x_RoL#$N&3>|M{PB^}aZNjG%&#yO4ul=i^N2Py{=PdoT=ZwyhicM|)8nNvJso2y$f^D~d zjgIMqNAkFCA4zUoFu(q|w&32o#7DA~rY$Zef{eC zwMAtYt>4L=u?Y<*cpoo#w;##w_m^bYMfJJp+7jDS`>STD>D7k#OTv~nd+ja`@t5Ga z{7W5A^`K$x&%N3oNr%3es;hstT`st$!7qjFNe z%GTq`NzE0iv2xviN%qozy{i^#=9;P+2Ck{fHI=#5u4+nys^6>jg2*~h)S`wS6txI8 zVNru@l~H9}k6@UDtr@$7F&Wrm6*Yvryv!As_iDVh$d+85xGHNcOw3v4LX>fnslD>r zbzF~#ZCMyTshuy0jb3Jh3YQtrtscp9)wM;oLcT;DKf=hA8}%oDw&?5a$s3_XCHale zqQb?OFhV6^Jkp-BN5XAl;#2x1vGKSyc_rgM(sCu^M)0TPQg#6&Sc=W$t>%6$$7>7N z+F~<#bA=abxEl_kTq>U9o8wu`n7ptrFIQOS_YpLaXnJ#JFVwz|zwb~3W{HrOqWv{kJuY5_K#w@91 zo=;7Vw&~O*4P9&zugOQ}akXAgnZJJ;G}QyEYe|2XW#cr}g~diupBG!qkL2bx)9S>^ zskgOiQ2EvUid$0m6>q>bI33lN%%cAPD(?v#oQ-Ph{1I4xtzUxO)A97r%DW2mtlBg4 zSWu0>#~snLnyssYvs1IZ;Md?oQEk0{3D~o(EWS)5`+JXFSFeiQu-j~1Te25fXRr2R zUwYZp!dqa8J+64WMMKeA)YYwOEz;ZQOJKKTuPSfJUNzn_J;snp&?-z|ZMZ24`vy(j z>UM+1MOZt3TGiMsYyWCoh1Kli)zncNq}9|>drLop3+48dpakmR4cIL#e7^ST+$*YS zrY2UDUW3&hfn`76C$a31J&#N8#A*-lBS3BYnZF(hdbriZk`-G`u`InZn<=7I*}-P& z>k(}A(r%syV7m{lw~Rcm%4o z$X-UP@G1QgJY!loPR(9&OK<9C#`Y>4r>k@Ay8BPt4r|hPGew(#-NcQL?4sJL5Yiak_pNP;UL#vzF z{|I)gy4gTWcGuTFITmc7w0rD$UCAcv2?LclxIM8OwsgPBqa8(b3S-}GL|2k#@0p## z*!Rp%VeEUhPI+<94GZ=(-<*q%$NA>acY87pLEdpxn2K)C>AJPW>c!L&4a9D`A(r=l zINq9>x(k?^e;`qPO0gT zx%8UorbHeIipyHfDK*c=AxPF}&hs+GH)X!}=VnLoY0s27-*<XQ?yTE$B3-;vCDdjm#T|i=jVo_mO0mGDk)*B2`v^u0=HkdN`*9Sv9Y>*oim4t2 z#vx{U6xNTK9fbxfW^xpB6&{lPlA^&m3R}S0G_%T0&{5fpgtzm6s*+7QDjUMdgE`+& zfJMlus!+%>T}P$eNrBzUZP}rJKtvAy?!mTXFO-Yj!XDnSZ5VljMq?v5H;tz6%?iHi zzJjSc`cYv04MkX~oQrPNd2!n^ZuO3i`Zy7tA(@$XTxLk4JZO1ygE|K7d(O;wnaVgp za8r2<+RmJqt>T_`p=72x3b1iGP1uCXa2}IZsS671WKNT*JU6AsG&>c4lQL-ps|+jk z5$`8x1*3W_n(~|^Y;lUF1zVgbuy+{j>|wfg$t^z?o4~1BFH3yD$|^D-=N(M(dB1G) zd9w|8gyzi+h6`*p$Qqyb%U+*1TY;&YH(PpPRqTQmZ1`!~ z@RAX~B{1_PFoTwpIn8)~m0^#di(S4TfE|s;Zf3N~%iIc}-x*s0^cs7BQ(~>c4A0C~ zu?>5iSS)N(#hYRlwh1d>7PJZTJ&X1a=(jDFLty)6Z5U=@t6YFqWwry&;u@$bTu;eX zw)PjZa#!O5?m#dQm926CTg4>;P0N3N6b}D}QP0ItU)z(}w3Fa}0%z)p{gWpMjaI8uXqKP;cYMyiZfaybpO z)63gFUh;g;!6~MHU}|55!J2HR3o^G+7X;@%YfiKi1#2V<-L@B8Q4@Owdf8!T`6dg- zTM(G_=i4sWl!xFs1X{P299{rXfGq{q*gE%VFxPEgc0pnd<(gaCxm7j`?SugP72m7b zrsp~(wt>Vt#4hU)aC?Vhiv`ny+R34HhhSp!w%8`@#VF8!*6!LR?r^atv2=W|hvj(+ z?2fhu?J2uXf!zX&&>rlgjFH3=gudZ$J(A3!#2SPKJOnn9Sb|W?hXC7Flx)6Hz*Yx~ z67RB8u;RvXmFpN38or1?S^_}LI z#TE-}BiCMk<%TQ>7PpzU94gzyT2e^7TbHX{nDOmf^h#mrq+S92NPA6iy+9?w! z3a^%_so^%@kXR+0o6b`v&O~RcU4X+ZGqP}O(Rr=*4B zF`P%80;>a~J#gZdu=QO!!sEj~*la%hGaRnruT!vhCOU4jy3hl>h*;R;u5&ny)giEP z*U=m%39U52v~#u~XgQqSVIl3F+Q2WK|6%@@&hF4L|BGjWE(LgoKH@S4qt|c#Oj$hTQRMg(O3?rrcW7vN1Cui%AHZNf^~wYss(WuhalUrD+Y!K zy@zB?KZtzUqN8Ca0&6i8ZNDZu$|Slgmi>WLA8VAf-oR2m~#w& zR_-G)Cmq+8oMl`}a+bNExtnnp_f9NI;6~)isys$^0la!xT9uB>(by`SesguTqd?Cl zSC+l6<+v<6F6%l4#cJTnYV4K^%+<+)c@>+jti~>JX*G7q1?Uo&R%4%Buy8ic1=o)P zJH@5d={{At1Xt*_1=aGyAcG@d0OtF!nI>_ zZ@N+o!cM(l{y;9Z!L#`wKmfp+7bGAlIGaQ#g7xcrtf|6Kq- zHx^*J*%=q0Ggy)7=KQ$8`m*1r347bhe=@%%7Gb*SfptNTZsaXMG=fV;bGLFq`0*)- ze&&+o9mTfmk)+^jE(mXqWeW>`!Lm$=ZrmW)BSnAoBSBYO(4&KT3#<|sj2H*<2+(pY z%#`SNw`OH_ghjZ(u?|y)p%x|95v)Z-iXQ^_qp=iIHm5va%Siz+ScxgK8!mu*ISNjz(dP*Sm@2B zpt^tjt3XGx3vA%C7GZ@W*~@Swb5Y7=>w*>?<^?hB^r-K0bf7;1v^of)2fd_r9VIrK z3$V#0?}F%1e*{IrHWwxCQ!j|24@I<@3sCVvF0ie5))G|2P!OyTlT%U;CQ<4LJ@E*b?Jg=4FCxeH)-aj-o= z&;>Z6Z7f@@ELt}?w1YN^3t*(|k-fHrr@$D^K`Vu0J=oPCh;fOJU?LVWNTP9+2789` z)VVQ+Z|d9_ulY3DH9qlE2gxotjH5JVf=NK@m`^+dBBMY?Qo=TWUj=&^@`;~%8Mih}4dZImtKJ1McJnDSc2fd$1z~k1 z-@#(auh7mduo_&D-9|#j$Q8;*P#;>2J9<|G@;1&48zlfoN7sgR2FMBUrC0Kg7Tt|uhz{1O~ z;MkJ+rsvYrn_eqZZ+fmh51(Z*SVg(8#VkJA4sH@MKz~ZZ)?K+ExJ(Me_AqZkjKtU>!2`)%rD{}EQSy>XK!yw9y4R-;8 zHdtEI3-z#wvflhIz`LId&|D0MfI$O|EyA1L1?#Yia+j&aCE9sBjuqHpX{lHaqbPSa z+y#eOlsnsh9R!D2l)DTq7o28M?rb<$rJ|C*{-l|u;$aeH1t}w1vw(i%CQ)YCvKMVs zxCDn{4kWJuZeRxWEqE zvHZA&4gwstWBDNs?i5(NvHr@9qxOzRyWDIy)?c|965xUmxLuG%cM7UCtD%hO3Ke4{ zLn-Mfu*Fz@CE2?t3XVwcEm-pupv72zC7HLG178vyvr|{V}HS z<{##NNav$K4>X58(*eB&mr3AN$1rDlpo74C!yM?$(*YN>OMdiny4B0et^NoOuX;JX z>ei;ytzMLq#i1@oL4I}&8r~SE#CmBDetEgsi!!?K+?T>WQZtldlZnT^yq)Z&u%CG9 zOK`HIY<7`+>ZP!sF0g*?LUgp3=xB!k9%&bUoR0QVxtRz8c-}E^?s>mhZgCG;-k$d& z&pVdi%kf^S_Dx)ZZ_j%P5yGXsJ?~8DwD;^~659cCIX3o_tFLj78) zd3)YV!#0de5rvL{$3bwhy`m48XoPhufE;w^%r-$r9K};>@zWDZ)babJKJmGL(hJChKCD(j(0-M z(e)_76+PF@kx1~#HM!dg6*Wz&YS-#7rS+;oS8&_d`E>dAQddS zRBv1Zrb}V)T7rp_3eDRDYao>>2|rDsy>P;ND3x0`Li4h5C!Q*;%pFIIDQUlS1d#Ti50ubx7Cwa{P(@BkBC3wM78-XA5RkB+u)Aw08 zB(Y#g^_$U-E?|V1-K8rw%X)LAo)1NpdZDqAdO6=2&uWH011rC`kf?n@f{iqrptM!= z)6-8))h6&(=hR)1KRy}Wn7X+#Cd26x(lEL2BpP$_UG;o;X6Jr?>O|#{>@-SRT4~tR zF^jH2d^0BiDS|Vr7HG#wm&zk~+QDso0?KtFf<9TM+{e@u?cBGqxMO)hc27Hp-ObF{ z!Qj+1OgD_7CT;~V7>bz2GDS1!CP>K08r<=wX+tyk22MM^a_d%fwe-eCBzX2Rfq87N5VTVLBXgd7HL?@`Me2P69VuSuggM+ z0jiQtytar?E3yR>+M`q+K)C{u%I@fT^U9!vIr!<`wwSWC8S<#mbVM#HUi+>JqDDCL z-2zsK5rBk$0Nw2*_IUgqAt4@5N0MPbEhAzTp{xm2)E)>-#VcflB2rALmJIWH6BRv~ z@0Np@iAfBrQqowuL)sK$otQwx@Fiu^B<7PCdM`1vfj3J`TOgw_`+I>n#i`i+1j}eT zR1i-*L-smxqxnpqIdMqJ`lN-Cld^lo}+iBlPt2o}(L3P{7tUqpl&4yEtl z$ayFe21LW5q!I%f4tw+fGaO9t5Xw+#E+QFy|(BDGgPuc20%^``rsLg z$R(bCVW*+&Kts->c!rut;S32{B+QWJ^A%}^=pdS&d zb4&^x@Y0p^t_w3&>Kg-9UJ+~L4>G-)07A_y{nf5bNg&kgi~eJ0GWJ)RXAu{MEw$rY zsaOC9or>6rMPV}s6U;argDG!Jv4ua?2_ya;Dc-~upgsph>fm&0;EpnO)ok3ssg$9A z^0MiP8JvQ*f-*RDC`Qbn?0#Yfr?c?(S>fbK3owI0e^=zlOrb)66`aC)04q2f3hiSG z0{7ZP04cbYJ-Pvp*A!R+9xYe*KnhN&3~Yp5-P66~_fQlFd36f%==#M7k*&Fg5KKu< zo<~rpa!QyBe|JvC_VCPO49)>I0b72SkM=B(<&!5&J3?fJ`TWwQ*0Z5?2U21 zMEisQHW+4H_6y}8S{V(rrT+F8Si9=*=6klihi?j z6(Kty42D-`D|(4yqMI2F0uMBDwV`S+aRx!V-9M_fZ@{2dRHIjHtJ+&QYVaa~mY_jt zUV;W?Vzm~?vG~rIzk1CPc!i{YT!RfNhnBEGyo!k$q$$;U$PjHfo%Yu z;RS2B<8dk`#Gtb@aqNY^3C3an)|tHU2u$5_B?yEA+}r&o(2zDLWDIG4g95%7Z4lN2 z*@A+cbXb2Pf;di)$TW{BS-b>cFJ#UUC#Xgxk`UM^F$UPEO!h%T5k!Ha3(WjrGjaHG zqX2{e#t}5Zz=0DqVCR(qCK&i?MSL_}(9G~4+pb8Z;};cMM1uZcm+%JlrO&|%bjt*^ zz&ZSm*wGCPLJRsGtiUOMVlY^NbD#j--uP>?7i6*Y9_ii&0*rKTfB*x^8!pwj9{>w5 zsIK1WVfr)9;~88)h9zAB13I%A9&N?@01W6pB60zBE($K-<)@eg1G*vw9MBah z*nrD%bPD_CXbpYF79cW(c{g1$4xn3l;o9%jznO6Lv%u4vH?60X>pa)z1;6i8atDBR_-LyXmuw}qfU!JY^T)IN`( zK-kfw1SVnBhzXSHCMHmLA)*3B1cFOn&bNUnbQTXufzajT1WK2K6DTGPIf1v{7n>C5 z);W(i7D;~rTLKDTU|)8qcCisV{dpX>?H1F+Z(zWn*^LW-=zGpV(5dJhAOlV+1CRly zlmWnibE2>t0F&qv70`W%W5|iZ8=BqggEnL}CS9{P_(5@!59lKY6Ng}%0SiZ_A?bj= z&oMlp6e;T8h=IbF0ukt%mW05@J^fUC4GCh2HYqE~FOVk@fg*EGB)LfSk`XA~Ym7iS zBj6uatd$UfkDYso4=C5fRQd`CV4kmnMk`;1|DG*>n}mZaM}^pa8+@fA)zQRX0P#IdyT6Y z$sfLo2P5|Q?{xN}0m8YW$)dT^;0FG$;P2BBOxT~3JC)SIwKJ8Dq8g8M6aoKy4KL=S zaM)sg2kCq{WNQfc=ZHw9Z!{FK#)~ePFES#3(>0!gKz&3u{qvwcupz|iBRcEm>IE80 ztUjWx$m(NZ5jbb~rns2|r!Uam5)pRE1&d&85cd`pGKqWf*SoLw>&#dhQn8?nMKpsh zc3fHJ>yx9A`Ti`r1@n5A~z3tRt2wAPC?;+ z-X`fD_EtxMStc-vmT*7pr!o<29BT0LEaoAohn0>Y^{{fM!yYQtfJCy&vv1seAtJ#jc9yM@x`w=k|>d2Mxwvc=bL+ zJ#BhO)RWF4Q4c!nB2Z5{i$p!>ECTg^q_YUr)2>KBJstRypI2BwFqdv{gAb}YKpI5dlf07nM+yTIRnT*N+3 z-4t(LLH2NL=Lo&7qrQWG=S;;Vgd&P7>MYUkri1IugrP%C@H?ukz`7$CNwB-YdpM7LkZc%zVnhZ6iTY*hA5g`jnR4h~~%W<#kQN>9i-$8_Z*_lM9r=Vvi8{v#mkoNtvC zYpvT6lYkqRkljLLUHxnau+BkGMX4pYUlBn-w~!eJHeJlLAe>Hl6t_pDvO&kOz1Vb~ z)^O_}pKz-1@hAc6n0Zb{o$|`>O00m`c2s!*TDM}(wlV6Iv3={5RR}MC2iM%Z0$U(9 z0{b5hWIc_;!bm#LJ-&mQ6{Aim!(`Mc71+KVN-cXiqmF61MARt>_Xe-rnRdaFB%|)~ z6j(y)l;Dc*pla*033Rs_)TtlUA{`%dAp2o`ib;2?O`Th(h9X$UntB~KbNc&Xjsi*N zTA4(;+u?_uP8Rnc8`;E3iez8n+iTR(DFfj$-x{rnbRE+Q1jgQ`lF-g zTp8MHbc9wxF*)xJR-Du87GRMljLsRFlPc9YaWL*Uz!Q{BWvedPNdQOp=*V!BhF>I_4qZgamFVjciq5%k!qHJvB_zj@ z(4^$Jrb@8UuBwDEw8l$ji$-X`(GgRWf_f#J&Ou{I$EFsuMi`w_zV!7D=g1MsrJDw; z6==FE0<{t&VL8};`o2d}B?ygC9ZqF9V09_9914%{BYmum1D;>#0|D!#X927exFBVo z6gbh2WC2WT4rBo$>zX^EM6gR*Jp2NjCLvP|E>$w|S-7cMU;}$!;9~-nskp@~z#UZe z<%ByWg(evb&aXQ&HJxzB@pOZs>n6uBy(8{4ey?vyZh=*Qqv4Xpz^2t;5G3($I9N0p zSiow926XV<&7WhnieiVw+EoGhEa-!OyM$M8u`XVLwWtOwc0jkl)=sOn#g!^QB}dl3 zXM$GCfZ#|sA}<`K$H=h8ORz8-ibrsQNO?N*5-cxK%Qddcq&i?zuwRX<(0jr&FR}09 zB2_e{bL(J#b~v|=FY8-0J@cQiUobm{p1`iRwQ326#en67!(!3$^fCr5PkR-@@_J}E z-qwU6c2X|b@DX0jBQiHiYt9vK%z=jpvnbY4dR{upcj6|i|z_d^cu z;0=$%>dnZrH*QxA^yzpn@FI79-*KEBKzbbd8P2zV1=6dLcHm?^XW-epp9LL%QvmE8 zyQUoe!AtP+XV;X%&#oyziMx`dw0?F?5lh}R8U*SZ3<5y}&0cdOK*{~=8H1cyFjW}f zd`o%kmr2e%bF%ncpZXWUngn(ev})1So^%W7*Q8tE{qL=Waw4bbrdR6vqqFxMfpzl? zA4DO41y+$)b<;qw8uJZ%o4~1?D%bIQP`QrZYY`1~#}EL6@j(SARY`k6b#8kBf5U|} zb^G1BKut=4kd)3~fsN*anSh(?%d)({o@N;Z+4)0@w}5tLmjomd{txOo9dJLX+N@YI z2lXNTV5w4A4O5lwwJ<`ZM>?utfBK|BPIT&juolD{f&h30C7A{NcYGhzoLGEd{-d@? zDk6n6a={Z8R z#g!|mg4zOQaDiRBW&sO$#epk83bt!jVh4TCltnPeyajQisR;TgzkJIm#}&#p%94y%SG-fR9OQ^g`h%*Y*vPRI?9lnCeCFA`@}JWI2B;0DtCgn zYJnSGQ&(eRNnjarl~&rQ=$F`rU6qnJqJBB61Lbgh8Q6#3F;8BxMzu&^bPytHZ1Q@c zLV^%caF5qTiP@Lzq#(#sRm3ftvk|+0D&Ro?`_y9cjIk1DpB@AU(_qQ>ybY%wBnEp zD#I2WD8LTnb2XwylpTn*CPPrq&?|8J5t!2C1yGQI&A6f|nFqXVt`G+`BNk8nvHUh*wZ{G7;IO7{g8hJVvU3G3`}tqDEAeP>5FTipuUFKbaumv%65} z>$N4%LPK#a5)H*^TVab%TqO5@=&F{C-J91iaY?nDdv-%gUQyrav$LEZ8I|4{sz9rz z+|+Lya;0c|0?qCh`jybh!orFcX#r14KK^o_yZ9oL(jk`TF22ZQB-FIr)ObZVJwpDs z!0yUcrm$+LS!Urd*Rd}tB`%l3jgh#N%=qKv6s}`- zvT(XN`h~O2juy`MB3g1?0pn^rz1artZR#ab$;xN>R2dzAG;9M$z zq`?1MD&A1+r2*lh`J)nc1^@-i83vUr0QUN?XwtpKhhUra6F8HM`9HXX==*z(R5zUe>`ob?up}+tY+&}>;Y^J{IK6S?m55eF)4ql;uB3jt$c1HFKfVp2N z*N6C*P8*bj2R6Q=I9vYh_@rha4wlo^qxR@j-CVMr=NO&ZU6C>DN9V)cbQ8wL zpqZU79X?DREW?&}+T{XmLUXa*ommH})ijr>Err*3yBS`CHjMZ}QLGIo5AcGXq07J=|}2qu*V+N6WY##=%oWZr`QP9i9$puTEF+_R&IKcNq%?ha=s z`#0t6WdC$%sk@(tmSsilaoD1Sc)IFXG-rk2q$R+VFvUl;kv%Jqci$1}@ zId%cHvY%Hr-tq8P0^w=(*Y|g{>Md(DKs{wPSX#fEmDLY3TYd?iF}Wh?6JW`Bx(Pu& z?1*|(DG=gFO1Sh!eV=azfz3sIA7HwGM?wG&_`^>Cs2={RLr=Cst!dn@4&<-HHuD8wF*(qNsTfS=vJHNVSiK|nsQ>}p(*g6 zL37mnG7nt98Io6!-N@1oMK?JNh$_Dzo1I_wf&{j*6n3@Lh|_=|p>(PVhXGN+EpiLu zFEUJzSg@L37Ddg$q97=hCWY#SKtNC^4o;Y`rtCR?ZAgcC7VP=hSC-k;BD1LsV1~XoRo`&!MALPn)CCb6C^@o~7LB6yW24wY;p7Dp#V9%xRPmm0G=U z1*?QqRM##p&kwP?|8QQPs$C$Zomd<_&rbUaM8=Dyb2t%~I=WvkVbzTaw>s;n8^+!j4TPDUZ!p(6EkIof{=-sGH|ei z$U=q~UxZ?&v~*Dceo(+fWV5EEz8e?w5nwcJ*V)A?4kj3Q&GKG4~c7%fXswQn#&QHDe*l+C}toi0*aY#4h^Z8387+_b1~C@ z&8!iEF#}Z*9L$GbhGj9C^WV9pAqrFcGNLfk-Bdls<8Hb$K*hJdh)Wt~x~Xl6!%R2H z#K;MHz$|FzW0T~yp>Hy7N)2M0$+;LDDb6n zZAHG^e3i5wTf>2tjR~{-abtcDThRf37eK>1Q*F4O5I(7dSK=-e%Sqg&qCrTzT;gFa zSsU!dQqcGwo_ybKY7l$rxh@QRx2NDqx&Wa@=02&Hp3ZVBlN6a7x6(ozibrY+ahDZ^ zCM!7~0n}SiKjYz%lP<+UyF+-BNS8|0ymgR+{^@g7$wI+|x}=^G>e8W}gMLkaoz_bh) z%fO#@S2|&8`$b(k+$qWpqn_A=pT_+~gmUEtN@4vSm{~Pue9mcM^5Q!3BLwiyv_@ zfxo$fbEnup3D7tikbq$ZWE%oPrc$*K&M~V&5^v(!Fz5!toV4}G!gPKoV3<>QpbS{m z)N_PLgDGeQq6B6^VOQzUNaNyUW{25fI{_v41mjnnsAux_${D&3|z z>+35A|1zjOKQDvwCH+!ad48aUA&6&F7bMV+sKrl`(3v8?18VV;pQTIa_a_Tc9szP7 zeK(IHNqsi9Md8PYbIbgHe3vroiXw%vgn0fSUet4%AeCAQ13K9c`2De7==jIqME#*T zmDy}IZlx;qf^~^j!ul%Xg;}3yB#z5Ob}_Xd{L{`Kfp=Mo8qQ)*fGbK~=fAfrsO@rQ zraLN=Q%6)L3j2vAU*i%clPSUFO&1K5OcZx(#2VaA#I5Qt6pe3xl61MzOm;!MOb)tu znW%k9v$OgoGT~DsEK^h?z%mEHi8tXVPeqqW@OHtX!k2BFpiFqCUKoUD3VifLQJBjUNazSIw7$lIA)4EXKsY4lw!C5 zoSA4y%?JuLr8;$gWCMbxzl*Rvo>E)b@=m1{gQ)4Rt7=~3sl;M{H3MO^^B!zR=RL4+ zz}3Y2i%?AmkOHbH-7x&YLs z?`6o^+^|U-8o)Lm3Es|)?0UyfTB+IJXupz_!Q0$KH`GfX#Vu5bkeiODBy&@cH&Hj0 ztZdlLYXgiLe`6hxIV|j^FqWilDm_`$%{vk|9Rvs9rg%mT0DX495LU#0wR{qXUsfYC z*x$T478)>rn;4mD3`|zcvsWr%LN=v@TZ@EDB9Cx&JCR3i15=%NB*G>N2*u*WcQ?VB zCk|l5I?{Dzl>M`OO+44lL%q<|;gPUcU~F>w6`aMf7@r~D)sdST35M7M-OxfHm!O-D zot@{ox}clgo1~L3Hc1CvY|^9)8pEWcDz-@{QT%s*xYF6IMc+A_wO~g&xN|D-iloh( zk*C6LN+lD!DV5yuEq14H+JX{fi6~B4mOaE#fy?PuCb68@OkB?8)TbH9<&@s+XAq3A zoU}tbgP)~FiRILmgJ8}gPG@i?i^^H(^CFtl;SfM`y4?yir_mxhF?k%tf(Q!Puw7s$ z4y3?;PF**g@rddnb>ci~x2Vz^%#l^9HAYsMQCN^&nNS>9Dg$1B#N9vKP6s^o2)BlS z!BJ>0<0qWOsSG8q1OTwr|`zrW+RVA>L=2l zxGF|4x#3$DxpmOz!8I{tq_g!ecb90M0%Z<=YA9AXrQa^cd5S#)oTql7eUc0zRwpXt zjLuiClyIKXdc=7~)FzGqPnxF_;fV8e3hVk*j4NQCK@*DQ35~<=4>>VzL>Hc5o+w}| zi{_5EgU>jgm#f6U@w7o0hYRqJ5#wnW2;x~B(aCsKYceUGk@I;(>x<$U!JlvvvgN{m zc(UbDYZBury$OtGWs0Xbo^V1Akm;~A2=SC29i&kZ8mbdHkvhZHogmMOg83-f@>iy) z7@mbs&QM+`Td%y(t_o(1Ariw=JAZx+Oeb{0csiI=M|>}ereF2y zLy=BNl{y@4@Pon9-(-ABe-rX4bDfY+AzLR{e^6Xxd`hR=`IZ2#0uiy7(aCLp#83O& z0;&^qNPtgxSU8o7uORktk(i<7d3ehIDB@>;!Y)zUtlKlV(uw?>B1RkX(_J<|e!8LW_JtaOSiL0qj#g{)@Cv zxrPAybe1hn7~NL{*r(HN0sEwX^Jun0oBydkaajC#Ry3iRHx2WOw9bG&L*plKrHAPjwBSa$v_` zS&I8ds;9GdiS^X>lUPrgYxqKl;{oT9&HhKiJmaPmI8V3Pfb;ZyD$ag)UG7ve+=b{s;plL2+uItxC`?c6PC0B-Y_A;U?r_3@1Aaf~ZW&_1=t=AX>$K zUEH2m-r-2D7ebU=#+wI4Rc&CvIVO6=?K|G)vOBR|Q)pOTzwGCKVMMqwx-2IKf1X?W zrGys_LD=gb98zxfm*=;j1(}v7$3Xcl9D__snWIO`aZI(C-l9RsXLy-dqynLcMe6)W z_Zg@m2g^r1(!v5b%p44P5gLEaL{dzDAxpyr&_=M_F$%r~wQ+M52dn=;qRvCZt$GFZ zhG`6d)Jmm^6;TO){NONDs?6j}(|Ra8IFGe$yV4lD=0a~I3m@7zptT)H;XEGlXi0^3S|NHr3Bs?9WBY|@ONplWYT z=&8PV56)>+ghT0_sCF>hU6syL5faDd0NNSV;MDdT9iP-&4ka~J+fXuZwChD-OgIwa zm~bGJq$+l5muduhmC81$S?3Y^(ubEEYpOU}P*Ja_E*(mOR;@GAA6(R+E|p=bH-@K} zR^FqZ`cl+?lBQNcv+H!R_%1 zi}!@5mZPQI0gjd;8aP@>*|3W@J=6I^7@5ISqc0eLG*$S5iBfd{W~M)YLn)3Lhf=Qyr1{3a}7ARZOgVM_dS(3d{=!tT6TsE2^64(jyS545rd+7@}Lw}QS&RCSDJ1;!H#veGBQQLXfe za9AsUzA(m7t@KH7R25soJAy;Qfz^j`m?=S^s)InQA_URO7O0?&of2#nODjt7|`(cIyCtl=mFrd{n5_p-Ce6LdcJM$V@9Y_m~{nl%rr_Ky|fo0jI%n zduKo1XZGW7g0ME;k`RxZw+!$%#?LCDS`5~RX*~!=gs^rNR?LsOV{X*#+gGkOZC^Qm zYD!lWP8EkL9JWP^i`asDS+JoaUtuZo2H4wjj26wrqE=0OJFTcy*zXp#3f&@vHDp+l z$lAC|1d=u6R*}f+;5nWZvCx9a2pa1ZQmqnp!sO81wAWm4GDOl1BQx+^ouNbl>`F8; zh4gR2SdogA_*o>U)@id$YF&cSi9sZPo7E|H+B@SlpH7}ca-M_)FhWVWOR_@QWMamq z>Y$^{1*SL#I|Bcdb624|$W|emRq88cY7-~1f(wZf6eNcRd6llu2xgT&GcK!RBVlk* zp%KeEoMo7)FgX_bI7m_}p)+r_(qTalleLT0LMW?KCkq{OTu=p?gh*HF5y!26mJ{xf z3RNbDt#Fg+09GiMAV8Qgn`Nqn16imui<&tPb(Tbw;7Bq(Pbr>+v`VEB)armsQYAQW zWJcJ-{Q5?TD6LTdjle3)Ey+&|FYTR>kFkzjr3g{dTU5`qxQ|B-N zn2c6eYh58d5|mS|b$Jmna71xYmGuC&_#l7PWtw|96N($h8>3yqsBCh{HP>MmmsX81 z2OS(%%ydfTan|IYg1^&636`#x`;B8Tt-WBh|!7a*YEuGd+Pn` zeLnB|KHu-3@AE$I3m>Zd*;(mecx%?})F)hITJ;K};d-*|?-U1FPV;3A0qbs#irt0o zj(lIHSOe%ZsEj)(8JHa(hS%E$Bp;x!9Is~s*3rV(08z;1(j9z(MqycP(8^HZWv_zU`?BI@bgv*ucOQ=5KmXCdaqx2CK+Mt=|0uwEO~Cl;34#3ge`mz_ z?=7|$zORe=A%@BAH2RH5KWZ&#b<-3^f%~Yt8tm5U*8NB^M(y{<*+gK%YPi>Q60 zW85&p>KF~BSSlPhs_ILW^ML``cW}O5&2owfc!&fGeHg9c7S=mp{3t!leiehcA`-;| z*?`->5;Q{bo$YNDg>{cy53z)N-DzS&rDEK1psx30uBLLpjq*K^ciqPqwFs2)-$iTI>&Z!nd88xrxgm zz6p_B|=ITz!8s>8` z*~#V+esx*l?qqgx`$Z0^A6y}cvCp5JD9F)YR5Y9|;uPG?m^kGwg*=2V{i+#{nZ9Xo+Yj=wJ)SYXF`kSA5NxZn9%lFaHSHRaO+?5)NM{R-5?C`d^)*wsg8e5^W>0`b8$XIj z4=Tx<<&}+UnH)l@WuA(wN@C|iGjy)OXruH-*!zVYPZD+o^i3_+oWHnPNB%s(GyWo8 zzY2^4J==GOt7WDR^7VO|;CJW%7r$hfmj0n$+no8lD4`@a&n1vvTwsawRaGX|E1GF>2xRBFR5pm!P3jHbfsc~( z@DxdoC7tBLvJ(^v#6T6eoK|hhJo-hKH+ci%t>gAG!;7^jE!!{w6;B6KJ#N|V=*muy z1kgedoFm+M8Tmnw-Kzv7%loMOFjL{EB6j1W{<9|e0*6R9CeQXr*&(33!RfO$7vPkU zgp5s(JClSpTsuYY@R-H|wN(zv7=m4hvd5k3ikZ zCvnA@+nU@HX)%Me>>602n%u&0@ENtH+kBVRty}6aT@G72z$Q1pqi7`5iT!cso9Av; z^A5-8Ep1rJQ*O3>!rlmmQm+#60Zx`*D?Y-D_(M0pJ1stt3Mq^OMGL${omH)lurEPk zzu5|NG-X7NbE66eGxs`wzg7>sVF3bG9m+hQSsLdWI@(=L{aV>CVni5&++1p_ei81v zl1*eqv0D&m*Iypvh9kj#azezxM{DrI{x?BcFI-td%(?F1IJH+UT%i0WjCOk9_yBQ= z9$63Tt-N3I$Mp+Hki~)Oe%a*2Apb77_n9!}c%LqP4E(IV)?ejdycJ-|864jicYTgYCk~A%Wom zSzuRnS+uQsz{f4S~o&;R`&|Ld=R`_K2U|NP&s`>Xx> z|31?H?eqWt_4@b!^y~h*fBna+e9d3~`>%ibmG+G_|W`mg_g@Yk!1ui5giG{#r0 z>(^g@)%4X`lHj*r|L3oNcm3TbYx_!dhDXci+2>#9v)07F{pVl*-W#H*86LI9tgFX;3u`}QPukG zv#7OwUa?lF6@nks>W{M4@1I4j^;21$VU0rYgIe{d*68u2eta6to8e6S!;=ZHI^UMR*stf;WMZuOAVo> z{9_H~kFv(^pJlC77+jP2$GRrFuECl-p8fHfT=!SXe}7G;KeQ&(TZ5$^Yq0b?4fg%B zto8mHEmpeAKdF-ss><)4M2+`bd)+_QUia^`*Z0p_txFonDqTO-UH%}d{QgPQNIik2 zxqhl)mMa*4TbRd_KUtT#evO=a>z_#(U+c4pmi4>U`Rk7?&iT?@^QVeyzV+7kPol=4 zwdPNh){HvqA3tF=jF-mxvB@>w`s({9QDfTw<273O2Q}(J*7)s{R^wyHY$F4U5asV{ z)Pr1v$CE!&j_a$;KGL(fwy%-K`n_fO>knzjeyXZ}{zFaGUyAD6=d8B3k@aI6>vtyB zPc5wa(osLQuj*Szeg7nCe3?N%wt{|V1pU+ok{yEu`C|(x7(lilk0<}&f*j?m_P@Iz z|Bx9OEvp|HR)<~n6SL}L2$gu2Lizie{6Wwb$Kehn#oxXCP1jE{vih)#qTO?3myPgD8Y9X@()oVi{Dk)$3A|0g%W&MA^#8+ zviMzvpZA_Ul;FDx`G=^G#qTPd`Ri`6LJ7)$KmG#vcc_oW`Um|?AFHLb)cQxuw10K| z-(Pq7$rq4v|H>_NMuP7~kAH~Xu=w5RrIz;|7iyH`+ba2&tdr$$D;@XN+o4iPzO9sh z$x2!Nw$fR*{SYdZi zjmG(r!pWn=&g{~+qxbw3d-44H+fiG;8#PJ3t!%x9jpgquEotl%O0Xs*`PQ7^lc59` ze^`p=7U+_1E9GB?^W-jnTWPeo%00gF8vBmnv3CFVm*1`C&%*4V7xZT-f4-D^joyFt z^tEcLcmLh`{qukM>tBBTPk;ICpa1)>fBzr<`M>|Se|i@4etN%XeGLD*_Bns8uRLn& z`wRFtfBwr~|DXT+Z~y$;um9z5fBBdHfA=r{*MGp%`{DdG)|z*-cc_&6`ZTA``B}&< zzU)uwtTWHgPSb1u^z*3npYWZffApQnGgsd?YVL~@j`8DQVX1icj&S!u7?b?F<=M#GI z5$sM)EXik!z#|h?+MY2cyO}BmYN?&=dH+~yWwWtes->Abi(ZD$#O-hXBly&EL9$=L zN^Ax{%Yt=`JuUNAC3Zgj2w2G;f0XNyurJemd$$J1x;Yc4_#=4!(?^hP3f6Adr(71M zaDCb}ljdw$yAmhhvtcBACB`e8%`_9c@)20OGHd-vp8xg{SgmYmh9oJ_|JQ2D`zmn?ONIEwrB2BX1DE2x@J|C+c5i-QRA~lKBd|CUqA9` z|7m>wp`P0Ja?Pf$+SmROu&Z|DnvJU`+eb5-#vF^M`It|>gkx^oDi_MZAIyD^G!q+G}TUZpyPt_RH8JQ#4b$sZ{&aCf7$XV{J~1RWV)t zDHnDVOG?_S_?-cY#qJ)keu{0;k6`fu8-J;{(rp|j=dg5MJFw+@f7LCu+g9C@k}o#u z9|_uRV?3n}3pCZzz}l$I^~DFQQA;C>#VCDQuoxGZZ4z8f8n-^s59BexDrbHb&36|**coouWH-nBS3|=q?grmndDRVO8g$H zX~~t;&%&2B6!A)SvpoNL+QiAK1+Qc;3tx%ps+F%~H!EK{vF&^$yP^Eg`vuxSc=4Cj zFI)MGJ4$V6RMz}$vt?{rqJ*_XZwyJcVYGPCB{a;kmEj|5fN4!*F9Rny3;s9*MxEp2Oc zZB4a`rq4d3#eJ@shF*o!*X*Ns1lG@XvzeT}hP7X)M4E-_ti3fhTaUV?eFOzJz&5lk zVg9&THd{@re{wz>7Juc!g@b5mJ8gS0tgTfq2yi>xQ0(>qhNfk-kweq)Ey!qY#kAfd zK|^uJ%`Lh;nURd8!Ld0YmM)>H}?Hce>AT1^&Gu!8`mSj@O5iH$FB4B z7`x85bM(e8pmBR!*{{Rb^pT+bF4&{pQP^+m8olSgXzV6tvh!POw#BTjvBSql+vQWX zw`=_H?R7I8OtS0|j_9&Q*)^nlmt3-T+F~UQ#l2Z9SwC%Wd7m;z6N4EAf0ldx!eA-WNme9TLXpNBbif1+}ye*j_=h-_JRX!xfG&!HZDR5CGBO7q zEWG+Irx@pWJ%JpnHQ0jSfWZG-J@^!G$&05UwMf)F%pXM?A92Byc8h zMVO{^ugBOYd2JqJ(&>STo^n`lf6}&#gwtcX%RM4v(W~MD_(YnG)1sG!Lup^WJ5AuF zol^e;K2NFpu}mdT>ic;Ij|zOYQ<{D5RqGy>d9z6-xH6o$?sI_;cdB3x6?}WswCI%L zhKKi>Qr0ukgMxvbX0*vmV=muV+-+KN3_L#HT}CGoTEwA;3wX-O4Hq!Ie>&y>j$R!^ zK^y|K2}5rlVR!q?%vRBqeSE_12AY|T;@D?4$^~c?CnIC-q%3D)vpD$Bz&8rp&4YjM zjYmrGBiK96lE#}*#Zqq*7NfM=gepD+<`sIp(T||gE5xz4bFurLW~HCV1#mubGN6nX z^bmpy+&g7X<)b{UZ=4G_e_d`lKxU&04Bx?Qbga1m?ZW&tI@&lLAEkT+|JH1Dym2xb zhsy;RF3t|jQLNGqXBd3W?hKnakLs>t!Y-9}Y}%>J1Q#T5k8y77RJC|NGCRFxH!GRu zxN?EnGjgcG2kj2EIq3azLEbOr+~Hh1Y#TJmd8Wcn;72Iv{X>8=e^zBB6DKD2ZSX)u z!Q`n6ot->MsGOT!kD^(V4Q1~q+sV#Pw$lYXr=j!HY@;Z^>4w8G*-i=;+u2Y=Gn*}m zg=TVPFSgPJXeF2WVk=#+7rh_rOywM9p8FOXg;qiUEwN8}7X)aheH`mXJMnGcdp`s` z8`|JKU+Z<2c~-Qsf5V4%3RXM}uEfH%OE}p4LSdgSP{4@_D-O%`5b%_Ep?sbr;Q(IR z>Q~~yQe~X-umBgfj3CK81KAQZ2X8Fm{DTF!ayV;Vkl^r#6}U3n83fi&tiY9dGz)@l zam%V$hzWAzkR7qRpz zXV=T%1+y(-&%iM#>*dPk8FX_7Za+6?;O29C25!AWuuqv& zXGL?{xrqZef89~GHgjTh{YHlowh6a+unV}!V=CP(9%%4WuouLQQVUiO^yDcpJAj}5 zfN{YuotsaTz`EdG1mI3eL=jfH#M>DjrFX7noZBBW|$; z^G+|pf`xWLzzfXV8j@syNdQyP*2)%GEiN#ZDQ&bZfuSb>dQ&kriMC@yv0EkED%oYL zqcyeye?18TtD@g*1$sh34%e)o61&kc6q(&{ff<&gp@=>V3{~b7O^nQTWW^%-FflHf z^W%bTy)hB^5ujIGGMV#3L9tui&p&KlSg@6qwSFv4=MdmFmFoxY`Y5QlC$??zk>D1W z3(D?2b#uH0R*ehnqBu7cO?720pTgT{3hFNye z>{k$2gIXm=Vv?_ zg1`ki*$=M2i3#&E`i1p3zbweXi&h2Kf1e^Ng1~Cv@|&0%Y^aX_lf&gV*)s~3M?icE z?3D`;On_2hfW$6XJOtMtB9@|{c>u*OB3{BJh&Uc`%e}fg1 z$9H@R>;bH`M}j_J=_N&`xWIx)u<%l14_wefK*%mzB`!di$v~6PUKJS-a;s730+wP* z3`46to)iS=dtonAhM_!?-Q+`nu#=nZ*KKoyJ(}r1qrWh6rdOTg4`r%H5a5EyPMXZLVjF;Hgf^C z>nPbr(765-6{NC_JQrVMqqRO8Tg}zi*lHJG?78|BNaO;GuWVX=jZ21te{O+1C1Iyl z!D3ao_PVX~C}129MyxZ$#s$p-f1V=ZdT2GK&VU#$zs^XXAh0MLEWZ)|Z9Xvt3Ngdy z5Lgrtm!f$;4*}ZDZ8^c?+qD$24|lD;!B)Eft!BH-i{)hmD{=Kr1}gvPz&Na9= z=q^ChxdiD3J|sI92Ln})n3t%Ie z0{9tRaQGP>31gUTSaZ@Pe~<+A@H1pYIeXC@ zGi;QdIa#nHE(cbp=pZ=U4w+76m*n7K$eG>&7vNaULFFkLVK0vWt`3Z<&edMFCU9Jj zsB06ThZyH6{QD9Ya{WUsDT+o8gBI`6x(1S24<|-$f9#eEco?@<_3%&jk*m?8 z^ChrR;o4(x@gcGLxcvH8dljt?uD#AKb?Vd3$s24I`s8IBN3Oq4uL%X?3LcC6N(k5j z_%twTgMJhjSU(rP>A`MsHS&Cpf}LYgz!{cnZ}1Fz^P?U(#B%MuQV_NPZVe3ED<|O* z*j6T|Mb=^5e+4;sF|y9$UX~Z*f&@MXuE%LU2d+oYdvF098;h-)aL@#a)ynlZ>1Uuo zX%e4*(!rqO-uUIU$O~QaGIAie{1!Xq0`!SqjO?iw7|O+NaY-&rkC)**8EluCa6y1R zVF8w#)8ztYZ`jb48$IEMT|Aw`My}lKiVM&ctif`#e{n<=T@;@dz$RtWGSy#y)xZf7ly~Yz!B)xI31&gTK5T{Dn4#m$4sn zcPuh8xEvc4;x53fU?Cm5W9i(sxH}el{ki(yKK}Cd@t3!czu;)>n_=nPs)E4k#rk_W z`Ac;2Tf0Ws8q3?sU*1ms@^XD ze=olYLgHa5Z!dqD>=~Ek%gc|n;8l2|S!76HRl<(o5`4S;i@W_Xc9YqVE|}-;`1bM_ z*%e%a(bIoHVtW`0d%MEF;g;p1d^`Hf;?{*gt49J4GgsxJhnb7;?dmTJdoY(E?TDAe zcgXtk_Sak4l{RuUdiyU(Y#YWhtT!()e}Vwse;c{_cK27A5nP7#?e4FKyT8hqP!Zn$ zw&O$|@}T?t_n%0Bby(jn|N3(ISHLjdlHlX7Z>N8KIsL0U{jZ~pby(j%|N8d%*UWYw zg6QL~nLTlV^#lu(>EMT8J6&z8;HGH!Pw?SlA=XSEt_wy;M}5jW2>p=26^(TWe*n-S zNLCBhVGWM|>kwca7CN6*Zo$c4-(LS(bbG;77~TFa!D6FMuYZ--|Jo(CnM+dl01d@r z?|T<9ZWUR!q3lk76P!N(szhHdz{)1`Y*GCh00Y^6SHMziY;9n1rtj%OmP7T|UkPqv z2!8trp1a7z0N#}CBg_a^N*Q2{e@q3LvdvPado*QRk>%cny$ctb1rOW)Q{-@^Y#O>S zMl|6JR2eOmiuv*?x9tSexpKx*u~o!${6}J+e`SbVf-z!EDCP)gbr*m;?z>>SK^u&E zOvReh7}1xiIX1hUFA_eA&aX6Kx%H0!j@0nZK-z}}{Z!Z)SD5N{jE|&0f3WR%WV`{~ zCZ*#iAVdd0i_}av4|?_!Qu4)7dU%V^E^2H zW1oz{3hZ-rwS`OQG0~z9G8^0n3Ac^|>n<(Fs-e@!isWrnLZq^*xdkB8q!S_`M|!H; zpk4`(a^C^-g}_+kS3SYZ2x@g0*NtSoj>2w{kPwIqb-2|He?PtFhwiw{0GqlC^CX)D z5$a@7j)-O7;lc<~RE+ms8T7p)ke~2ClQ+7d*!Bi@8Pfn5X-5xAW~7|hZNL@s1e`r= z&FO^-Z?d3iMv;lt!T>J8OMokMfXJ;2V5AfG#`Mq?@-h;9 zhxYe5q^JX2f2fpia0^2obSjsp4JZCui=+Dr9uTU}%m(cfLKmhx@#Ap-t~f}E>DF0| zZrBEj{-z-TS-xum3qQLDN-brA9~;=W#EePyn>Fn7M2CoA z0lzR=XypZw^+&K0-2A6oxo=<;`VK`%;Ui(9kYNhLf3`(TVY-#TVxU6b$4DxCBsd_IXN9TZ4}VON8d zoX#{!e=-b}OGt(xQ||K1s1(6asa|#w6jHm1G?c0&(r|@GpgRfKyP5LAhJk75*(qQe zdJ+h4MRpZM$p%=7g*Yr)w`3E_HAy<5Tzn6m&=qJY@HfDuB#zL*RKOAHXiglVE7ao0 zQoOV>plFvGgN1+rY>W)TQYekVCr%QK;U<^Ze}kTgx!HrRaSJ&yN&gwl6ekN>1fEnf z;DeroMqj4BNB_Yqz64jEOWKcqMkP>fzhY3{1|SUFi1Ie;0Uf-ydS>$%)R~fnjXM~X zx_Q+ebxp{hUB+EOr^D@)ZjtB`6kb8e)Tqd(8!MRZB3?>;pCAYqd?e@q=z{4MwhjIy ze{2L7M5Yj>0dU?8I{}D-{!?Ux3c`T8h!Au|gER8Qqjk$}9A$t{?G=8w2|s^;1qrqz8z!`aeP+TN~ z-3BwpA_X9DYE=Oc=&+%l0yPMHK-Xy<>6D`SkP+wzLNWqTIM5Z22ee7Q0w2&`CpsW` zY*{D;G%0_@We!^46`*zRoPXj19U7-Be^$age;#y+4|KqqxLoLVd4b1zfAte#$0iyW zTYv*P_NwzZCp!ahz(AM+9q+NpNh>gm(LXAvBO1NIN)8&(7v!f!qy6Rth1^>33T>UhZED` zQa}YxCkO6@?{ltS_vBzsf8Uuef&v98FixN!6~GB}LO`Dc0tg(qm6K=UckSFsk047M zN~A00oL^_mFef+==Z~pp=)fT-P^Zbx-JfFutTqCqHYHG#PD-HWoRmNn6utuib61Q0 zp0O+_GF4Si93htC1C~AFgaZnhJY*;o{rwTyAU2>(-H~+XuSXEPe?ZDv4zS`OD*}iz z8WA?hm-v7(GmH;NL@@Dy`sIzmi&U)p5iU|~*^Lq?y+bxeLotVhAm9N5MRza^TR+q> zWOT4gGCDHfiG&Ms)+!S&=y$+5&7Zc50tv{jfQ|pnVgCeddY0z3r9o1BO%kIPQY>as3L~ zhZPKFA2&Wi`?T>5^7x(V$0G~cxWrWf!tZh@3K@PrPV&KUf1Jd_;M@p^0h472h~cC_ zo`4!(3j75LPIm;8MbZTbet{fz9~mFy_XVWat(a154O1iwRmaduuhL@Kb|ES1A5~iA z7sD0-@KcZ>34RI+B(cvmCy0F>K1XDqdy#NArdTo}`gja8K%cA5Ne;7g=LYmCIEsWm z#WRu6=cxB7f9g7qvorJHM&MUbn(m?i=QgymoPjM4<`XP^DnS&(R~S=K^p5Fa8=~?x z-j3y29aTOR1Hu--eb!N63L=K7(%yB%-;vEHK3^BE57ckHAZEm6-k>{_n{d8?0`wvq z(CAo0?;DfL=z5#Y&W-%f$Ri^Vt??KsupdEVjelT!f0k9(gOiSAKH1tN^U3O(B1Ut< zyLePxFCimKFCmW#3y;3RrFTHN^5jSopO8@i?-5uzbQ;B=4N^H>(4ae~-q1~Vvep63 z$K@_(%}L!0-kYf(J8|G0a82wy%}DlrQXg{kCk$C+WLHDFGmvZ9tIOFb7gbxZHYZCX8*N+%`@{Dc| zJ{N3Hhi}mKG6Tzm+mj`_HZ2q)uQd`#$@4#tz5${S+7g4K(MH&c^(NY;$9&d zM%X?5T!&Hi_){E9rYn(6bb6sM2!$Y$57CjCe~-iI@8~-O6|&0$39!=w39!@R2CxZ? z8dMkEZ&vviP_<;a2_WF@fU>oAmhdwYOCf@6K~a`Hhn*%cn^bu9z$O){TMUE@V4iMm zh}{i`0AT-_g(Z_%OugX0;GnxcCUS3}0QJk&NZXhV($4@;A1@iqY*FR)GnYm^k(5iG ze?39{?0Uv|ku8!S#oyyJHELm~$0gsNmZyF^py3aI0Mr9;qHqud*yAfZPQy0n{~BS> z>3BKNex0w!{&ei}H%h0AB?ym|hAEI20jKe|oM6Lw=4>dfBzYWG08x_Fze(O;ef?J(|>v z4MiT}`Eks*YMe^t&^uq6QR#X23!TR+fTf%~Y6(xV!$>H_PFTE-v_00MNO)=jhu-KN zG4NdwJlz_If~OGl`#zPVKzIuOz7^?C63|V16eZlGI>e5~sD&*PbEhC>V(vT)f0>*+ z#SxQprzl}^?il$?%$+zL!rTeELCl?K8l>70Zb7V_NEO7`$um&MvFo8708NlX7qASZ z*one$V~aMd0uXkhB;4}sa_F)wSqxn^L$>I39pr_^vZ>LoH%2uog57Xx9=s4DlNDE^ zqS#H1SYZ^XBj=!KMU=BqpK5ZUf1`qpE(C{82~4-4K)rUeMO17j#O@9?z`FAqpw_HS zHT+JU79Kw}L4gFjNMx*{ALxCgB@&tUf+bceK~7tx z79?;kVhWqGEu9Saq6;?Khj`!li*f;AC%u52`KmXbS$YK62~t9UoiBA}f5eONad(Gx zfI-J~`wEQ^UT&z9T`OI%J6_JQlV|*d_;Ppn18jfHx4gBzIWZ32<%S|~AC7Z_ubGQy zaF^CVfG0ckM%x2fsFB}sbNk~E4cx*Kx4}?Ar*~>YDFVM{?bY>vTt+;$9JU)(t zW40|CjU@JrLJM20@Ny27f7Aw9^ge8mZ1hB~qqmtUoZjd&h4Wr7n!wH9{zIX-Z|C~W zsd7#m)P{bX)zIyv#*Txt~yA`yZMOc&xtap{k{mb@nI zwgS-=V&-l>!pyV9N3b9Bi59kRr_Sw5!QUcwZtww^5rhdb3vH?be@nk%VI7H3Dds#S z?bB#~qCcPTC+*X?Q6Pnm9I^Bj@_7iHbNW_#)3a31;Ff+Lga9r6RMV}VmxJUiZrvT;iOroICo^Z~ z=q3KP>a7~EF*-_b2G_I@cG8ZiVxP56ilEa_#C4wI=13*L`|Cj^NWRyzy^)ral)9eu zz8)io1h~rx{UYEldwY3u5QC5T`fS09fkq`Uvz@d$IMSVWe;31$1i+IP1p%+%-}G;D zz?*WyEB_pakjx68HH#`Bykx0ONT=q)X|H_h)TBmL3-hTLA0Y~a9qA!$IWhI!aeBP9 zWJkI%8rUCFKRFsC=}FP{S3)nzdah@?;`BKgotl&#=~ShIFFi70iQ22aJ#?PXyel~Q zuqz{?CdVtKe>4qTDR_bOuaffvGu5ERcV8EpQ=^CIimNu@|{$* z{2(TE?bEmm5-8i7&m0jwk6@Yqi(ka)Z3f`)3go%8e@?%Ph<^zM_+8k{MFTEzlV9X~ zas1_?W-}LjHLBq7)i@teh9&ucu0PLEN-+bI_lu#G4lblS$Etlmg4zB-R}T43?&PkM z0=agG2#mfRK?DijI$YcW3OY243Hd@`-KS%xSekDDSiM}}Vh#ofm7L;2^aZ1o5JJ(S z(&+{Qe}nKqJ;{&GCk`k27PJuYAtYgO^p~PU*e^KVrz(XK50JvG*jU&;3u)IZU%YCh z-$~tMbgC+RIi7M3PNytaeM9jGW!-{9(Hge!a;Frrg;KTT7Rsv43S3#Ll<}~}Jv==9 zz!(P8mJ~zj-ke{Ktqj7DbZ_>%fW>!&y8tl_f92-N#r5tbjv=k#{0hJ&*HC(xSi{bO z++Cs`t`HjY2uP1bm&yyXnt64g2I7D!C7jF4BbIfUfL4dQ1VIgi!k6Dq`^KLXX~&%i~4udq`0X?4DRHaYlgtSD?2%BROB^M1l;nBr02KK$1#bVdLSActt4DP}Kct4$S#Abmw!@l;YO=bGA<02oSpWPKG>LmDxA7!ToWk)+TDBl46Ny8(H zjMP<-(CQsdu;)zKW z$YxFwNZ71N`uFXW2A$CFh=36vjd$EJY!UvLOk;oP>4-cf5Pd| zou0^x0@Si-40Cved}g;0k-&B$=E4oUETgdK6;O;~L{c&G6$(nYP&TE+WiK1i5suJ` z-t2McwIzS1_ro**3-&rRr@maWHa)!rJ{?C{j$^qnlPZUHe`42rfhG$I z=Z0s=t!y(vlUgi=v&~Kv1=?UD3N}2>L(%?57GTj4;pQXU;TE4ccQBGBD+=>IL8e7J znA&~rV4jrAc2{~^_P|#j0h}v&ZfPV@R)HhWEv-~V^Hjb8TiR%6nKM;&lim0V>~VtI ziB83C@FyzNi`Qj#5vT1E9x0sl+@A(~ zDG`#{&s|AeO?B(dmV}_1Mwr1uX4<;zmnm{G-ZrTouwW-$1Z?x9i-2vO)CZiadEBIt zw^`Yq zX{P7p`Ub>N88lUn>6WU&?FZBkWg2q#xbao%m1 zPH%OOP)P-?L6budS(X{5K30@FX+$< zEm$D6>0SgQ-E`{N4KT9n!*>0K>Zg%6hrP|r zmqEK9$j$2kFA}KvY}zC9ZFUzhxD?3SbP}Uwe`w;bVc?AR#q#r4L-S`)%|}48B;c?j zeuAL66GJcOXO7)AkTVApoCwYw&ZV@k$tjEz?))kEClFy%^A3*Ye>pHyP=g1UsT?*EGw1r{v`pn^BxEQm^^CksMIsZJ zsZe7AGZ7C;T;{q3PF)I}#2BtgzyMyk6m(Q6fZq&Z{P*p@*FkY^{v17ye~3-% zaw?WuZ~aZhfh$CO7CnL#6Hzvm=IE47lt@KkVxb6|N^P`=vUH2jDCDNH!G+vROfzug zrqFM|ZYps{*iGdy3cHDq!=al>%@K4nt4E|638Ep+rXxX-?GI?&7Ll{bGsPpBcoF76 zgz%hB`4kh82u?r|xii7*Efa=ZPfPQ+mUR-xICNq? zH56MQL7D$|Dw(!moz9&*i(fOEVi~~GDtm}T&k1+|jx{omE&@Frky1!FL7@deilICm z2?c-P6xAOj&p-hvM8cx5=UDp~+-jZ#3pbf3{KL!T3IFvbPz7Mhf17!i5pj7sY@JDRIZX`?G0u{LY&c6@C=~RSzd@?c zxxS)Y0Q?LIsZ~)cfA}P;!cwTc>W~xu=AX2XqYogc+nxY|dX9db3VvZbCfsUl#-u-* z%@97F!!dPfN6KhONfAL*RKoi5i9Z5F3t~?Tvxsh!2RR+u_6MQ4DF5^jiB)=2cMOju=4=Te9O@3lzlGi2cNEnxLX2X&F;gc0QS+(vM<_s)aq$nVl<&_6gL*_JP#JS?**e zqO|X1!Z`tme^Jh|BS!H-2FNH8p@fVI<2iJaV4DdUl?E19RA&ouMQ>(51h6Pd_FAcA zK}E%G6jW3+Q@EYNiNu(qkAOqG1bk$B9ZpT^ssAv|60DR*9fAdKn0yX_yfJ)I_`r?g zU_-)ai>TCo>ySbs<&EK!J=TwtPfxUz#cwg}XaJ@5fAeQxMIEi-C4;K=f&l88bWPo8 zU=qn5CEc&Tf~q%qG=%(-J?g5rD*`&{mAGpVHk20$WO4GIfIZbM0Glb=6+x52T{*qf z@=JYmYEaIdcs~@Vi}5z`tT;RTA#@&Pz$KN}pIA~?);+{blSxWteThlxycHlxU43_j zU$D)#e}fPF+j1}iO*A59J9VniqR6i1;?j>hk?IHSkXO0UO#$Xj(2>caJ zoqak&1Xuv!ii13;W|Q6sH4935kFFN5M+3pu_I#pDCm91A@+qPgDWQ}5hG3+`x!imN zOk8$g5q`kMy#Q};2JWIVX1A6pW3HS~hfsmu-_|Q(Xk)#0%Jlt7EAA#K_1WYh)F{bRFQ3C89 ze^^rg@r?~8y>u)Cm(*dsxB>|zOZl}2pV%jz{ z4|qfV@~W41dk3E|o~9a(UrL$U2LUGa6QN;9?Lau_aaS83bH~qz!%XH&)Zo#;kYfWy zQDOtT(=%}k;vax}@G+`RO+wqcRnL_Gf6B72_6~N4*19MWY_XU~TnM@p;9?LYf{S73 zaLG;#Wmq@l?jY3!>tT@AgsTP8O`W|oNP0j4hf@yx6pX3z&qXRhPv(JV0@>6L zyn}Za?8RhwZsTaJz-j{AR53JUoBC;ZCC z$TIaT7U4_9Z<_vN4V8a^`)*sJ0`#>S&9OjPREbtkopJ4X*k0h;&il9V;QI1p-1WP>p`#a;h3eBu~|YmT(TMAuIyl zl>PFb2f~u6ii%^u^n*c}?`Sa~MvPNA5heqL)lCHkH&qK$*g>R6D6mkcBC|HNi~Vi_ zVmH2{dEl>Lr=5ira$hc?PG$T~1?dK!=^bKT1UPlNYe%DG{o>~t7J09of0!-yv(vZ1 zesEhhaQ$x@)j0Fj&g5Z+VB$J&`Q@Lgc zHFYH8ZN4H0l$!a?0<3p*eaRF~Y66V$z%V`RckHDB50340^6yV;Xs1?9uQAaa@d|dM zWFA@kshdZPpb&{2u4T8Eu$6{iw?uT9n@zTLy48g3)7g_*L^gC2f8pxxlgD=Ih7^Ls zK&Q+!UXp5rLmohajXiZ%WqFxY*d^Vgzet~9w~I`yhVp4(EwN^tY1-`yz|pF(^tc%! z+|wz)tA+~#B)0%{DO3z+wgz+SqGZirZY|{3%irkvNtHaa=gq61NYarX2Ykkk=Ka-yO-Z?m1&+Gpocx&)c1 zms@MtXqIz{JD3i5n#LTbSt%OGQWP20@tplM@0LDHliWzJ{qEQQyhD1iy z^$p0VG6fmRA=z`=pEZ!hB z7WdrAdL%g2G0*E|J;H$^x2A(sr5Qw0ood%G#NwN4+_+QyB=At#V>o7vgKOn7=dMAY2XoR9N{|xO) zye}iZF)2l4@uldB!&%{EWhq*!#YwD5%_MS2&0Oxb-AKTYGU6mw4qcWRizX7|=4{pN z$MUpp(pbC0R_&DzTRqHL(pOFYSc}+;6Tj*d%@Q(fe@S0;q-MdYQW7}vs{;l2c{yl> z_*f2FiB>S+s)rdY{;EKog%^b!xZ-v($`iJBQM*@kFK+;=;78=I68%|EBl#kJ)mfn8 zi*s1T7k7ed3rt(n1ec@iNMLoEsD3z*Bdee1S!yM%=*#)5&Khk;fzD|T^NT(X{;K~X zn*glBf9A@K<_OnDB5tE)FV0D$`r?!{YGO`FqxvEd@S(F-;aH)BsU*7i zJLn;4YMhb|3X4;s`-zXZg;V3CbZTWzN>}0va<+(gbscNoJHGJDoFm@VHSBP> zjdL(~xE5-ax>MSm@3R{ z^qEO>i;=;x4lai4;i^Iva-7A-Xe^pG=SZ&?EAHkG!Mb){?4}E(t$%YutdzA97hxp` ze-P^$dr`l0pq06})NTv2cnVJd@c<>3o&yXy!LU-@1j9Os?xtEK;q8{w!`ZPh2pVYE zqiB_4%o(z*A0Z?Pb_gvC4@grh5i~*~TsjJNY;es@8mv=C6LdYvf{6|hg76aiAcP!4 zzw`vOW6={}jFN1Z=xHa}F2PsmC{vFge^b-pduPkxS{3-2RM?wFPNkN2rN3m&qmN!jSPe4rds8Un4}0gj!MOC!NGWq%vI*hp+U--f5RA6 zoG?{X03If;22>UVCn?LxuN*d|&^5LS2WPG@1@+y0B@=bkY$oce*%mA1aM$>6rRk#0 z9u$Z+zsRo_(XUF-OAH%iVa04cWtwR;)D0KF4aPDuwwa>P7`w#m{v&OaL1Z6mGxqF} zeV#OSvNOg~vec^>Swfl=Mz(}X48l{EkS!%t>f!BuKE0>zIrp4-Crc}*CT4ZR4J{iNb%T-xX!G@{+w*ZoZ3+tyo<*7_L7JImaMgy# zK#G0@7^LO;>sCB>E;ahq&g?U&BQ>M@pj*5N4eS-zwGurWblk?;j5PMxF|#)9W$N5X zW@RN*Ho#I2m|*PkWH#4BoyVjOZ_*Xs5MU5es9IrtN1!Q*ttMT^Utmk7xw+;dO|X6n zU_-f?zQ>iSsP}3?SNK{~uJEZEe@ikfUY4lT!1%+;9=`AnxHy#gH?mecoQbU;=T>=T z7KWQG`!W7L|FK$e{G*VwHO^nW&^1vQ^?C(f<*x%>?3!ui^*-0;GXVOv|GRNC*#g!< z^`tJ`ARJlp{(2}nibOMO1Gm|{(R_S|B7P~JlJLzrma5K9xXNrTl1c_~TJ4PLUe#Hi zO38|MnsjLD)D!W#5%qIN%iu|uSfthJk^QU!mMeZ zE`6zPwT#q%B{DeRgUoiH@F!$mq;zMt_rh_KfzBHqYYPq5v3+N?{=iXG6(ndy?z zJ|z|V&?;%fN-)@&>fu{FU%_wISzek>EEOgq>U_OF-j)wj-`cudUpjpAOv+*{7)Acmp@`O#bbk)Oi4^a zV(B$z4JAN&LkU}jxiu5_6L|w8lAw4*mnHHl1d5-oly#J`Kt2i21Pk7fjowejW$nrz z<(_CeAZMFM$KE8?Tg3D1Vgki5Rj<5@Y-D|lY_z;P4nj7q3pjDJ`e>Ysm?Cn9b^`YU z?>+%l_YUC&FY%+Q+MP*k?rK!z+XjHF7s@=A+ra+-sosyxSanwmbTY1Yw-)Yh zkd&dj4BJk%1SHlFv1W4FZb!Eu(t4u{mf zjKh@%%>(C|T!kc-NWnog9j4zCV+eoW)mJRo*f@<1a#)cpVuU z1cIKWfL=&aVgsn!#+({LOTeXYk={E9j)Y(dh}3quYLn-b;v$X|^WaXc%#kkb zcB!j9>&`O<&U*T7EoE}2g3n6PA5Bw1Ewm_&GsX#H@SE)wUvNm6Z&xc+zQsA46B}Zs z(eo4x43#nCkMFg{q+Jct1UBQ{XL#9Z=|*p zY%Fg?r`XpyCPcI?h)E>%t{)eNuhcsx>Jz^)QN#oGjTMD5vfU2kAQ*iybVK%i+`-QK z1i5V*v(vKLv?Pu4ZQQ_RFu~YSg>>NX^sg_eEttCuKTuP3_s&CS$zPMyA@pwFfIetH z=x)>ylKp4_ey95y6|fARg*U_by_ABkAJ6SqMpxkTh}@Ox*3?uL_eP39fqr{&$`1@1 zw$rVv`lL^zbb-m+K~HcNZMJSTUs;xnB9FrQ8j}&%?Y4R>SoVd>F8`=RjtZB4?dMUv zGMZBwf#L`~6BP&F_rA>$t!Htt_k?1-sktr%+#gNDEdAl!p67$O=B1|L<$oH6S#iUt z4yD)fyiEg3_mdsS)~J3p;2s9)vlJob8A@nfZAvXghauMJISu0{1dws*xc-<6R~RNx+Maf$H@7ImnP$g#gHH z2W7Ikr+eEL9jOJZbzg6vs5T1uJUs_lzE0hC8@w3)RS_8Iz8>Zn=D+TvO(v$_rI>UI zYO6QJZPS4z`4YYh(vyS1-9Czq^%;^d@2%`=M~`irwhXT<{_$m^X?OFEMW}|W!h%ZP zOxb)P{B1-?I%H6V<(3`%UZD1}+1_~zuBScubPNzc;0uShx6oi7cCProqN5iyBM9w8 zT3xyC+2QB{VvW3>xY7|ymmOG)jnKU<*L%dwR8x3So>craPH*_ZYb2ZN)A*2$mKs^n z`>65Px#$yXH9X6SF}oXly??PA_wXmy|B)+=xCAk`aH6<>2`t_Hxu*RCy!UZu59BiK z+Q)cAI?RV9U>he%Y{Gsb9)bFwAOaDPe+@K%(;+gdt}4@U_}_h4vh1`7a;M7kLSheJb@mCu?RLjDV} C`^8-V diff --git a/pres/img/tophat_timedomain.pdf b/pres/img/tophat_timedomain.pdf index 2808dedc7c84dacfc5bdcc80877b57f34aeee6b7..969bfba061274af0b0f85525de11ea78899e79a0 100644 GIT binary patch delta 1628 zcmZXRe>f9(9LIZzT}xsK6d%IN1vTN(M0D6`2{6K!bA z&%#D3%HvY&$0k(DHCDDuPGm0Uadgk!b5Fj1eV+H{{rP-epU=18xahS=UA~KVa*{6 zCx(1$vU#r`-SE{)wcB{yNmjuA8ytxdC8M7Bs#jU&*`goiOphNA1-I$pI`uESd_+9V zY4&i)OdT84Ts-=vYn_actezU0BHmlsT-sck_?l6Bo>qWpybPey@YQA1SQG@^i@zy? zN?c2Rugaok<|=(PxA0#Q%OkoBxFQ^_5r51(^_gv>hh-H6>co;lC52C}*Wjh|&JZV+ zcWYQB(OBccc;_8;cT=AV!PSwWo>1vL3O^76XCiSU{r2!{(#-go|0*DrVhCpSTDM&srRt- zW#K})H_kKgT1u_6+B$A@oZVN%y#(qtT*>)RW|DllM!DmHWUy4Cn80UyMdaYZC&Mni zpn%KG9vu!6O1nKzsy@}UovLwGCQi04?NhK5CU7A8i{bX> z4)!-du?~m3iAO7*E9t_prFgldDnrCN*F{Mxifv2*Evq*HXdWUaDd`0m&sNs!RkZkY zlUav<8V=VJAr-;0nm?p%UbqeBSQ<28D@dwXMaafQ5CB>O-IM2Gzzgf2cS zwD(5xlzpaosRzV8+PJdT2C%@NydopGtK#X1GRd%z4gCgq;DrqFA7J}*4yzKo%i-I? z?_`r`2akT`0(+&ON_)!=E|5{8ppAoFK6kPG77OK=%$cJgx+zayus^vWGA`kjUcfj~ zM_RZ9;@X;do|uwvoc!&zaoSXBOB#Pn-aj&;jP3bqR#S8>GRo0E@tCy%R#p_e(2eF~ z%90{qdcvoqNWOF6*pP0(#@T%1!p{uN_+j|YW7%x)&ITcq@iSm|ZXkvH4$ygX%SqA& z-s{_n>-q?>-vOo#VEWBDx<0Ul!=EC2JQAHZu;J6&4|=V~u<)@h%D)1{zW6X^W?HRorrTwvO>hxRYuY-Qw4!GzE zqB2!VUXKBho?5NXjT%dDy6-bxHhRbLb4*pz?alG|o^nk47kfRm-b&*C6+_Y-91Stx zT@jzyw@r{e7wkuGnp_mxH4aaQJr`435TP7HRHq*moKniH0Kg<38Aptna*Uf<*Qgm`$5CmJGWwMmjLI;`q{&Un zt}8U1ThTa@Th5lmI+Mvc((cr=`#k&b{`Eee_xrxz@0)i%F^o$F1bpi=ss;+8nJ>@~ z@;1udzbxL#kgmuW9~$r&G;hJ0T@CE5(dArPIp<59phXlxS*4a=jvkG-HpZ?5ZKzy* zvg-Eq(zc2AVVPd}`DS|Ok)@L6#>WwIdplxl$MY_^{F%}-JX3X|i#Wd$I!p1(Oto|< zF{7p?;w8idmHCW|K!s5u)8tnt*^-6}xsj1=*+#DyS9y}z^4%@)*SDuX-A`-6yPcK) z8;s{q@my|BiY3Ep{j}SUSe{0wauBDGb2rnJjRC@YWYr*+m8wTHh z@eR7;72I5lk#FUm^BQKK;2qWASio`5Vdov}X1lBIbWTG6E**%@F-5XdUaeFKy~!vP z+;h6xxyP8e>$*ITqPo(K?;IuP(ax`*2;R=Ofr|(Yv6!Ti~&f@N`7tsc2qEOjrJCG`%2i0v=eb(=bVKdZXvBwd`&OAgV|FRo z*7!<4a@ojuDK}Y1FFvkfv(1o7M3e4dn?KHOlwqf zbH=m9EA?HQ*l#!4VoW@fmESZ95+Qz&z)!<u=p8h2PcbkUIP0i-Hgf!P z+jO|Ew#L-KQGV%k44QDIWwN-1I>_!vmj_xKfup^2W=pgC-v*AB>CfsaUKG?cj=Cz; zczb+sP@luc`wc$SN0#FA+&^0EPVu`BAGY}t2DR1K^OnI-LN#kscKq@o%PF2wG5MZL z*=tDLPkNc9#6MhpCIYapha_1h!E{O5nD+r0&gCRJNa`id>1!5%PsO_4Py$^Q2jXsB zjxdi~@Fm?U8fY!1H9H9)hBX@<4=Gr|V~R%dp%i~CBZ{_RAp^_2XHE~6@}&6f0~Aj+ zO`ayroCT0c}4Wp1IzeFwZDl zu*w*2HT%r!))?hV-K}sTT*;~se@OfI_otO5jX-D|nGcWdvkz_~%sQlRc>QH_7=hp4D!TYY&dekkX_H>== zez4U?^^d+-(erO#B}-ZFBRhe~l_1#T5P$>f>#w#v$>cxTFgGOYMg8_LXYkbuaA2l! ztm)0zGW{|%D6II1ZLF-H67xdD>13Q&kW1`xlinvals($xmin,$xmax); @@ -46,25 +45,28 @@ foreach my $pars (@tophat_list) { my $num_elements = nelem($U); say "Found $num_elements elements."; - if ($num_elements % 2 == 0) { - # even number of bins - $f = $U->xlinvals( - -(${num_elements}/2-1)/${num_elements}/${xres}, 1/2/${xres} - )->rotate(-(${num_elements}/2 -1)); - } - else { - #odd number of bins - $f = $U->xlinvals( - -(${num_elements}/2-0.5)/${num_elements}/${xres}, - (${num_elements}/2-0.5)/${num_elements}/${xres} - )->rotate(-(${num_elements}-1)/2); - } + $f = ($num_elements % 2 == 0) ? + $U->xlinvals( + -(${num_elements}/2-1)/${num_elements}/${xres}, + 1/2/${xres} + )->rotate(-(${num_elements}/2 -1)) + : + $U->xlinvals( + -(${num_elements}/2-0.5)/${num_elements}/${xres}, + (${num_elements}/2-0.5)/${num_elements}/${xres} + )->rotate(-(${num_elements}-1)/2); + # Output frequency-domain tophat + wcols $f, $U, $V, "fft${tophat_count}.tab"; # Calculate x offset - my $φdiff = atan2($V,$U); + # This currently multiplies the imaginary component by -1, + # and I really need to figure out why this is necessary for + # proper output. + my $φdiff = atan2(-$V,$U); my $offset = $φdiff/($_2π*$f); + # Output frequency-domain time delay for given tophat wcols $f,$offset,"analyses/tables/tophat_fft${tophat_count}.tab"; } @@ -77,5 +79,5 @@ sub tophat { push @vals, ($x_coord >= ($mean - $halfwidth) && $x_coord <= ($mean + $halfwidth)) ? 1/$width : 0; } - return @vals; + return pdl(@vals); } \ No newline at end of file