From 7ee4031fe13d4f73b964b1cfd9f7db41f0026f7c Mon Sep 17 00:00:00 2001 From: Thibault Barnouin Date: Tue, 15 Mar 2022 15:50:32 +0100 Subject: [PATCH] add option to leave WCS unchanged when overplotting --- plots/IC5063_x3nl030/18GHz_overplot.png | Bin 442983 -> 445472 bytes .../IC5063_x3nl030/18GHz_overplot_forced.png | Bin 0 -> 439462 bytes plots/IC5063_x3nl030/24GHz_overplot.png | Bin 351340 -> 352073 bytes .../IC5063_x3nl030/24GHz_overplot_forced.png | Bin 0 -> 355040 bytes src/lib/plots.py | 21 +++++++++++++++--- 5 files changed, 18 insertions(+), 3 deletions(-) create mode 100644 plots/IC5063_x3nl030/18GHz_overplot_forced.png create mode 100644 plots/IC5063_x3nl030/24GHz_overplot_forced.png diff --git a/plots/IC5063_x3nl030/18GHz_overplot.png b/plots/IC5063_x3nl030/18GHz_overplot.png index 89b2b3a635108890d38055b63a391ff141c01de4..7b95a91eeea8d5043295a3371f130919aa89552f 100644 GIT binary patch literal 445472 zcmeGEcQ}`Q_y>$@3zdqHM46c>MJghDW$%$t_TD2cl2CS5LS>V(XHi6Dlo7I`?Cias z*L8ot$M5+4@%{h#=W!hOeILi2KA-El-tYH$zSem?QdW}QM|PBqgoI?Dtc-*z2?_ZR z5|X{6q`UEdNR8h6;XfCgCAFQ^>~1={89SPhC>lH4TiZEXTijrBHFI>bu(LhOdybc% zo5|eS+1^Q*kI&}+_Y1ssj<@)J8qZeXL-yLs=s1y(9Iqz+*=F19B0{o_ghWa zIkh=BPbvz3U3)0%Zg4Zo=G{c{(s1cw>HM4F_*>l(XPe?pPLkduIr)G7Bp(P!Pujus zfBxA?PHTnN?EmKPkw|9vw5 zZ#=-{1IiK-I~5cZsE+sc_WJMt(qCRtAvq;4A@LxoX7&2lIKgZ4n~9>nz0Jl4Lqok* zW?7?n^!)Gc*psYVU}^8@crPX<#@y==Bct7j@}mHmrdhkS?jk$8k*A*nR+g6oD$k0D zh~QH>=(@c=^$8z582Yy{QIxOiVPl+N6F!poq1)wjbt7h} z2zc^D+R~D9Vq)U@ty_Dvw6xT7^yxG*)Wv3gf7-{uYr(|MzGvxQPtE;32c{=F&$6(y zGjnl~ySlpCyShrXWoW3Bx!MQ`3DF9xfBf`m;9pM}p7Q%lt=y+!VcQQMJ{-ZJvOOdu z#NNe4^8W69{X;{h#0>-l$XxC0?0k07jW>$dAh5Mn{pv$Cj^mfy@9a5nynpmZs*Hky z1g^^3#%3%z$j6To+S=Oh>+AQPvuc+ZtPM7{ zu`x)SdHL$q2@Vd@q@<*P;NVYjf_B$!Y=Uom3jFZx+df?JRs3T%dG^eiLr)k*Yk#C( zn_K<*`STUr??8Y5oik<)FN|yaalw%VM#M4ut|%iziu-yP&1ZD-{P`pJA#VL!?74x_ z(Y;(;T=zdl#>71GTKae5ZMl!g!+n{rx0+_Uwt^)F5*iZ$B3}ML`^U3AN_ez$XZdHMRUVh1v_NKTFKp1XUqvbc|lc++mJ&mG`3 zEDw0_U}wRYMwa%$V0uBRa5lx}bTxK&PtWRB^>TMPY%fl&oGcA%M~jY}gxa%1l^%1u zd3kwd6_q_mj9z{Jc9i+{Db!pZ9(^ zE;cs#_QvW+UzMoTlj8o~orhI@f!|DCz34jezB%ct*(1q@ z*mG9X_3Vip5jr|LCpb9+1%)ZGZ&}1B`~Up$!)`k&u-_-~9V;vAK>vsPKg-pB5{H`-r62BR$ls^)9ed@p!*7X%?)uN4?>u>O zpgl|HK5qI|)4PHK$y~#V(@uXC1nheChi^yLH#AglZmi`Q*X)kJ=%zFfe=RUP{BW$G zT?ke@#cithXaf6XpH(WXeeItw(GGKe6?g90lj+lFkgm$=?%`pK(|42?`2cy1lNQWL@LGPybbTKri{~@whWDalpBP*mZPuf8*+t zAFclV7WZ&vZrJKqWkJEYKpMV~L#%Qyq=Jun3|cSFClOdf2Z$@2?MM|5;?mNb>&Er(H(A8|Nxl$GA#r1Omd@uW?dbiF@|!gjz!P+3Y`x=>=A9{NY z+ucKsUk;2pd+Vftfcm7cgM)*R>jZ5G!$toJucZv1lto;-O1ADf9_C%gMQ;b;(A*gP z{%U%5wmOiSOQhI|0Vm%iGCrQd$f&RUqlc)dDE11SFn5@Ws_LQ2<=^DiKi_P_N#IOO z3=6x5y>Pwh`DNm-<2k=_N%j;Yc7NjS)xB7URv`xi%AF4$K4iXh>Db!hWX({0cuUgN zJqob`dw1;F&mBQ85D+l35$!vX>$A4_Hpe|ND#|}UU!dH5W;>3skoWTGBfR{Gh-T8KRpp*}m)VwwAAaVJ@heqeVG&a_@XX~k?Dluef=ZV*n_K}g3UyVC&bNa?d zQhdJYBxVo|7LOdX>z(Wl2pkJyo?eHhHoakFtMwI92nQgsDEN(l#0z z8b6$oixdn;4*HOsIdg_s)W8N5kpTSaNp5ZmKEv|&U!wV_kMG~Ve|LbU6N)`{^SL+f zGdFZXB<*@i{jQ0u>@2YAXe_&1lGa&ZMT4bI)ja8Q8bp&uCt5fq=$fItIWX9Qk4`b;ew_4QMHzt)z)-1v%n}B zpN-Xp!#*^lTix$&9Vm6348+|W-?RUiJKy0^>^PedO4F24PEoXA^R^6XJfQ!+pZNR~ zsbKr9O`qU#&cpQdkI{9KY`P1D9DYCGu8YN9SYGP$CDs7dngq+t712gYtZO7T2)Y(m z1a-Tq3;NBo%U1BU@b~X~sEozM<#Ud(a}*qT89B*1zS?2sT)=o#&5kaZ!NOw;C1)%Vm^0H zW7H`_J$?Iu0|z3wbq=8iqCFHhk@39sPZD`~`GmnsXf3rN z43E*IUp3|B=VQO_L5HG{mX_8=J^B_W=)SprLHz2~NZv2#I)RD4n{=%5Q6GmBeRG_$ zBtN0*^Sg|xmwU{y<#J!VNROqh73}k7_VNYRNq&pnn+v4IP<<4!g zva$qd!Ip|(>woC)U;X6^d9G2F)T2YJ^V4MWgE(I3GcQeNRq}3pYJOvURXtsG7ygp_ zbpN({_wGqFl{pRX7k2xHvu#AhrKvEkt*WNRZ{L3<9{?P^ozNl4X9Z zjT{Y5yd%fJEcj7CJ`P8!y?6F_;8Fa3^56n~C83MuyOg;2y^%5*#;CoA*=U7#cAC%X zCa0vdWN2``5|lwH-4F2KpP9)O79Kt`(!#7?b(`VBg$ra&Ub91cg>kU(9FxWU$b>&l zjgLRZ=E;nwO;>#-m!@05*${DRhn}9^vuDqkPMz9|okxn}?d+WYy%gp9UUW3ga}l4w z**{;#srU{a^eNulSTbsH8LKv2*%+vKIP$CHEKx2svvsL0`y?bJ)R&(rT$7N{jXP%* zWqPtAGLl+kV_q$Y;o=cI9xYu_tl5VTAB0@SsL%&f9M7+!kP=<&XkCd(-Ts!Zdeofi zSFT*SQ){E6t!+NPjxz9=_lQNT0wX2mi#eu}l$4aNrk=NE_%Cu4+tfI>puP@JhZfrcV@vXXgJ9=tO zno3G@p-uPp=>EqcA+_j)d1haZG&Tkl(ToG#aLYON6+5I%)X!rve!ew(oaX)Y8UN|C zXTuBlPPzc1%k&P3`I+>Vdpr}gd!AIobduTaPXo0|Cw8>j%fq8Te*7R#`mqa+60=jd z$+|j5-ypI3yOZ#^&29nDHkWZAM)rGa#T2|s~Pz$P*-lQT4j%0#l1gHP0bo3 zV`pr9*JJLlF_tNrdXH^y`6FPof#sP&vOjMBetp0$&T|>p%)Wt7`+&P)IXtRYVClzK zwSLH;)cI~gLW1cA_E7u7^ZlcDkPSw)Kl zbWYK;DB(Bi!&u48qN1W`FMH;5Pi7`hF6=m!=(9-Y`tMh|vgTV;+4+AxdfBTbhY!lK zeiy{IR&8xw)Ge~*5q{`_qH%k3=#8 z@$~o-#Up7`E>F1=Y-bRwQH=9Vw10rH0Jb1`RKBMtuG&xaUdCzFe@7#R-7HD=A7
!`ES-0|$J?{Ry(@vjHNTL^e;ayb zqZixY8)jJIXx5sorj&C1}}TaE9#yY(EWy#O|~3$0@~bXY}k4go(upi+Uy zeZ<-xWMs^pj>i>VCwe-HXnx?_NJ|Q>pl!TqJUUbjqwo4%Tv+^ON>@u|;p3S3&5hngJ%r9hwgsFi4WhM= zD)8ylOv>mRl&BnWF{)2k>dDz>0N3RwT7M02faOt?5xb3l~o(k1=WWso`Y z$@n1y+z(YxytSR(lZ~~d6sO_FFqfkO=68ro$JW(;cQd%U4<~K6ot<5_rr7FeYT_Ha zJ`uk&^z`(^Q!}XzQqif>&bvWC;MJw6C)8Y;{pjkJOOr!G4-ziAC8M<-H7MpmlfxzQ z18fEzzkJkdaYAL;OF@C+2#+p{oF_BjPEEamN2p?r&mzJh!gR5%-JR0!qUt0dTg16dZ|sg3qQuZ>=mGDuX)t z#(TvP4R#l}uB!Gn0#_0^_~zGVNq@Okw$|6zPhWpeQXgqfoV-Xbtz`6k%c5bxdQ|J{ z)erW%FHfH!7%}MLiL+t`KNuWtAnjQxRRwqk7CUBA^)ed6WcwYPOzY`D;_!>tfaJLkW!6T5p%t#r&PoTq9>$Bv~ zFMIg|*ftz|{Q(?6W?;*&Vm+QMIXU^X^N4bb!W0`Dn|7u5h0T>=#>hwWEBhFPsnCPM z6|DQ-7~lU+&rUXrE1+N$(T|X^v$PBZ1|N8TcZVhW_hGFfTZ4)2q66$o@mH|v@+}8_ znqMo@t7mE|s`D;c_oz-5ShkV<`}gnd;#HiXspy|Swf(j(Uc6{kelLQ5WwNO64NAR= z$1Ax=`-L$L@0m|Dy4fS9H*fkMWR~hjg<4+hwllWW&G+AC*IRxnb{`eV?X3;BSC8gP z6ciQnEZaEqtvglci|nd4GzYaI8l3js^!{%nfmUm0ORTPgq&Mf8NEakdVeg!{>TiB| zM=Y2TXbEo&Q-RN#YaFq|T3=W_oR^%H-!8<1iEF+L4LQ zcVU9|O3zMzAZZ!7Ioy~K#JeFY1v3s&`isOy61 z+2z-;uJfC63JU7WRRskFcTcwT^pv^rdo5b4cwkG||NC_zSgRoH1{4|s01`U|u*=yt z_YH^qmDQ1yC#@O8rUE2*92SWcv|pUi3BC~{;v@X_?OUn_v`0Us1YwKTkP|U+aYTPj zo=r>RKp!fa-a&{hBv-H!d5(YNg6PifMopIdxBIicD%Vp@pg-sb=dK8phI;a#t4cu7 zY0Idx={|qScT+#K`Gs{fB+7y&GCmQP7){vWUrAg0g(J< zEdgoMfPl|0#Kpy@e}5uZSpnwVM;ScR;F~-(G7@m!ri+#%_8v;w{xTj^(f)l5qMXngpG?g0YVWPuPTStkJrIN-%@LH zK!U>RRl#S^UJ)R+LnV$+^d;B1#YJ96!95QiJm_zT=GexQ>^6X(JVGhfjdVW z^ZmWN_u#>nwQs+gw^1HcX*mpp#OVYbH9S0gTQa}l%B`uX&UfrzzXdyWTtrb&iFx&k z6*>{K9%G91lUj%@#Lo1a#m2Y}_?ZR-e0O1&*wc8~Jb-`g0?SBihsrPCzpF-dUU{4h z%$ou(lKfZ{qE`g#fC1-Si{=+%Z;FaKN^CA*zaIcV7B@9zQB_sFZe{ftx0_<|<28p? zPJHj`w}u9BfV+QvRZm7n%wN8KUERdfK2l9^UGYj_VE7w9hU7`k03c& zRaK%OjWeh*BzGP^Cb#U!34W~@cXTq(G>6cr@OAF9L$bXS zGoK#QNTYLLA$EZ4eRzNHXf&5rjtpCu!n4!cuqYk!i6{u37FxF0aFUjm8R9yqNX@V0 zX@a|tUUcPlo$Ae2cI|(^&D6{+eIZU4$ScjdG5*3))&+WUD=Vww;^F|?kfxqqT3waP zyQ9?PTM7>ilyztMFVYYN+h=Hl8Cjf9KczO4Oml8@l%4c($ z3UjbXQ$z$6ATV<-50F{}ufb8EBuPNqg=Nmb=ctg(kmi!I-e}!!>WLaMYB8|%bF-nH zEg$_~x(%K{$;kKwn_pW;=L1Xw95N7d;F*EOg*5zTd&CY4BWoz19Y`7yaR(xo+OxxkCWV=pw?Ioc)4lEm)y|It}d_Y&O zdU{xcXLq5^;WpS85!Md%irpRQkDooG{Q#bs)P!PoM4KJBn2Md`1pbqo*IVKw{k9Kf zp}tq^`t=v4Mt$CM!40L@(-B!^fR&wt+zNg(B~KVGYPN3v%rntw@6I!MnD!jn70&c4 zwMz|^P7qr=(hXc>D;TFL`=C%%gWa7hE-8_%jaf5q&x#@7l#WT6ir$z=xu-Vq;MzX7hkLb&^xzYlUe5($n>?J}ytcLLW00IEn z`Sa&Fy9){nj{qFALr9eT+|Zyp$KacZ8vC}q5YJGyP}J?8enDbph)G?@UP4{V%F0r& zWZ$uyf>J47D9Oc0wT_;^-0qEN(d-aLB18dbT)#Xjq>vHstNjiU5`~Rf)Y#0R*qXq?N}J9KWzM`A+j2c zau+0-$UZgvB*7J_z9zcCf);AklPlcajbcJdTpTwKPh{<2e?OB#%o$1d#fccDsr$bS z%iT|WO}MC)rkNevu#^w;qZ&%q15k#Hp2}8NF{5rZn-Ryx^WMw$DhANY?r(1H5~BpP zOkJ=F`+{;T!5_pZYM5CF{J0CnM#T$#P+domlT+?A)VMUQLTjDuCsjdr@87>aQUNXc znV8ABOt2IE=hbA_&!41dj7NaEQ;RxM_xel?z6@b> z>1mQpXfmH`b}pNmn&LxU75l><_3G6t!kID#y%kqd3X$^R2!H;Zmr#V$3@beE`Uu+o z(*9ZQ@lL#I>{!j0Nn&{xJb*b1#;frBp_v>m82I-uN9CC4#(a2MFpe#!PbJ$_WNY4o z^H1*8+G%XGo7}kZAh&W=4lbQWS7k{_3GsymA$uvlDG%0;;e>{A6ap4PUk^6Bd6NuX z;84X8#l=ZJqsoZ&EN8ksbMy051w~MaBNuG@DtYL>fR<;vaIe0wjyu7gKEtii?RBA1 zGV-!B-MH$WI+tD%p&aJdUsftJ9Ofn_tc?Mt!@9HF7~zv11YP5NnVLF2ESVwJ<{$wN z;WRYJ?M(&%D9qft`7&_exZdRAQImq=QUR|-7HYsf8h9eQ+3cPr68+jl&3AJpbz}u{ z@b`?Vb&;)gCzb9#->nKMoAuT4+~htbDJiMDE-wo^OGEa^F-RX5vcHmgexE~c#M7rt zsMZ>^qPM;8c6{me!BO7}lOUbogbAySRa@VDy|>`n8h7Pe+=C$YT^k+6k9cx2GBO-* za?xwaA&ko2wSnU(2Rk4r9k8#2Tq9Fv0pl5k(Kfr;nsW*N7ssL3og6I4Nx|1=t zdrDoB!wogqJ^V{KX1@vc$t-%IwO7ouyEweUIpEy0Ln?^&M8+)}8y*;8qpXe9H@?=- zqpgJ9-(R)4nyED7z0jr=wj<5aKt}7;3o%HiMv=1fD~yxfMI9a_mq0s4I&u@GHvWnL zShE9iSo})Osh`I)QjM8@IYz+EY~Pw)22@=TFbIqixk0mFx(dLSra-H*0K^t3CHR&q z(scFdQ_3WC3f-zHm>)-%bd()s<>UfQZN)%-aySErdPdUfYd(AsTONEOn(rw}-qm~4)RYv= zpoQZ+tRdM<2U@f4&++lcj5pU;cCnHC z5fM;)ZB%TZ)Ih(h*XLPW969+B?&h&t>I;k8{0g3UF@}vYir(g)pX}Dk6uETyvLA4~ z#rjG)7YBz-n$Ek|QX7R9KiKTTNUpN{NWAQ&*_z7sC>w29Yqd*&E^eeXJ>0xZRXN#M zMecQ9rB9jEFdfN*CnBpd<2wthuVXU^+Vw^R*`P}J?^iTwxNU^ms96?_)v7-)P}0dM znm%bg)mx#_zJ@a(xq$41vr^t)DdrJF@>^XJwVW?b?xR)RA%dYP{g5 zx;k?7rl;}O#d{uN8K%(7#G$-9};idPUA5l?>Nlj%RYtPP%-h3Ig6b$tl z-brM8J-nKk!8*x5hK7dr&dw2MFL^$Mi0*%Yk(ZOuqT6kNYD^$yCObk&tKj9y3W4(W z0xcDqQgd1_Y~Yr152u#5;bXH)OAk|1Q`@_nf#8X=@-yEeTyvqI01`I~E9+=$ZKcmz zbH0TdpsCU8v^!fHQ@*r9_Dn>R(~GMeL#>k3)_$gBbIEg_l%Q$w_DusDrU|{DO)<9r z`zv|fwkrF9?cg&(;3y4;uHCqCQrLOq4zQBBey;Di0{!Hxdm(lo2rk!EuPISf8ob_NZPBJ+oaBnV{b^pKN(wPbfsCabq$CvL&hNu1 zx~i!e)jnAcj7<)uu@26oKp+p4e<-2+gm$j;yIBe=4CUn^p@IXCwiMrXbUbsjG5!&v zHz#>`6b~FcT7|Pj1Ai=e7Mu4oNegW>5M=#^9x743ty+dHuv&xXt;YC+ksh9&mJAP* zY*3NaHLGs>M4;U?o1FS;avjno#qmqpN$|1sW;dEYVgVgcF8$sDGrhp`<{T&13ML+* zX~I9wa_9qNTZW29C=bDlUmAvnhOQ5N4kZ#Z@IjqzrwkAPsfOlsuFLnPnVB(o3GF#c zs<-dnnFUKz{~K&|`~7h*VI0CT`kt`~QqYWwl3la|kL@aUGv%Q}$to3(Y51TP9tj8v zYRxb*D6sB41hb0FOm&UktnvBf%z3-9LiFwsARUf~a}coMvRc*|XWd@;LqiCDDRUGg z+C{bupv$ilU-UYmGSWmwBJqR5BLHP(1v0dR8IjGwHf*zH=H*p-TaI$t<x;+$mz_LM?BuTA`SFfxsm1vY>&p%;SPzqwwgRWC34%blM^V>0M#}}@ zL!t^w0@s_-d3`Q$9#suGdX5YU56cdZ)%jD{(o9`Yd--<|ZD7R3-_I`zQb|ByplaEE zcw)3)_c4$_l)Vf}@({4a{JjN4hmoIo%neaIRxdsWT?L6$N)*XfZZEPk+99d-28dpM zzy9Lr4{19e{CkgDsuFc%B0LP7oQSN$l$3W+h8{QIF>rc$zchzyMTm>3%C1ZQ4A3-_ z@Dz_0xq%Q%NK11?#KgwlnE06|{3$&X+>LM#;68W|M2{eegu%OK&%^RKk(a#1qF##> z5Xk&s>B>cMr@K~xGmd@E>$o1!AYuiw7$!|yLx zi;jb)g3yK^T!(5Sj z3T<*3mJ1Q_A;ieuLqbHK{)FF=*wKLsLVL!P6bcM-b`ODa;A7B{H93IbzPtI2&^SzQ z-68@VsCLPN->|ACSWKJ%P6b7mI7TAm<}%iL7uRy^*s<%7z#+su@beo0Ob?^c^6~ao zE_Jq8U!Q9dUf5U8FoK8qOg!e@yR)$Z7WY8#Po6tRO~l<0gZguh zQwz}XxK`W%S%^D=eDC+@(O$w7K%Y9v$|^bk3&;67Of@0Lp@7E3TPn45`vsM7?025%C&a3P71A9_m90e;?xr8GM4;s4-4z@dXs=T#5_ z$?nclN%2pcK=6ZeT(;+#28yQ2MI5f05~pDz83C?U)7yJICntx)sPYnFPeboYk_~6m zt>b%$BWVKpHF?&qa^-ci3I9j_Y0zTHdLr6UO=IMpv@rfuZiw;h&7g+32hjeA6v{+b zf&AYjw*2$e_-2B4V|{yL;Bca(x0!wOE1`z=@=6?{j_h4e~$ z#?(sUNdoVE{`xgy-HS+1A}a%_jv2xbR5PO9aSZTcsnFAZB25B&E7x_4IC;Rr-z8Fx zp0j*kQ?moHK6gY!@#ItNQ`!gXp0uD^5D;%tByjT}QlMzdq*$aKB;Xyp@SU3O zO2m_e_f%Y3dh)^rI_Tz7_V#Z$LJ$H3C#S)FFZ-IS+>)Ajn=|(ME#UOzzCyOm)p0(R zPAK^2EC)Y_GLy+pAmsuF^Au(3LX$gdiE|5N>9_m}r-MhL2X;F#}h zJ-hCwGq@`iJ;1#73?3UuO#lYU%UdsH!&#!Ae??>hP=tBqISJ+N1IM#mkDv9$O_q;_h<$0}oSXe4pKu)V>OIInXHDgeXHS*MQ%!g zWA6gmt`~xxo~udW7`o~M%+E8`hksL}V6v`UK(i%7}ssT1F&`)ZT3^AWQp1{q&kgHl#3RM zK3P!V;^wBMy}cYtS~c$f0k&BEk}(hl%z~0M2B@Y3kjE|SPCkTrDR-LFVyL(H_Pfh2>7+9Nl9ozZ?bj z@S@tTDhyfeU^;f$Q+>Px4)?`^@gs{OykG&T z8*#0-dErRpUIJ>RYIWsua!J;8Sd`j9RKw3BUN!O zP|iKw4U}v>^(7RfUG1U5eJYEFW~Qbb-ih!bPGw{^E}aINwz#y?ovEFdjDUN>(iThx zQtXk)OO=gZzfvKt!M;-FVsB4Ln9QDxW9&+Ox^h;U`dV)L6h8Z2hgp0HaDx!sH+UEIfx&(qn zFJ7?Me>f=IAhXfiTJewDMP%{c&szxD`k`^CjA~mVq1Bu$NnSry1BKPpJ2FF-;cA7~ zzwgMaen2aF`wzJ30c1t_$`D|8QV>uLj=9;{Lxe2X*^d3z*9(^}}yLu+uAc^>%tTee0eHO2Mz zxSq54F0qU}$pZxQw~K~@G7J#W$hj*x6{#&=Pkrepuw+k}iHzR>61MW*ZK?Fvt5fEF5El&=qR(jo( zNIE<^nyz`-j7Aaah&&X@=z*rYja z)2G~Nj4~ItP=#!6L1j1bSZO&LdOEiET3&NyVA?QAM-sx1x`0$*HTpt zBcQ#?7XE^I**>$i)nC$BoChcqk(2tB-a3Tb*0wZ1&z(=NWB050gvU%rP0iJSkPuCs zE#$w*Xr`SQ(d6Y0_q`mdqD*Nr@zTh@d@FpF8_KiV#-6O*2H~Sli~tjdbBwBT9Jat> z4;Qo*JD4OGssT`#2JN!Z+gh6vZJGP_?b|Lf_qcf2_GxWvvz@mR>lfn`m6VPo?2wX{ zCI@1qfga7ej*zk3VI^#pQ!siP4n1U}@Sxy!iRZ^eLbhQYqPLPZ1r|S~%m`*47;=9b zDG|14pjUJ1DR#KZq4qi`Ef@zgL!;`&)NsH-W|qTuoNb7E${w5}gtItiRO%P)m8zWk zT@G^u88fmkZ7Oy+iJXG)Dh{tga*~L^PCiJ4_C^X#gr)#rd7xknF>38uuk%P{M7jVe zN3fCp6(M7e>BRwUoyxcPq23Nq9w(jvRufQJG{&+rbj`LXKVRL$F#XF)V(Uy;4%sRN zyJGCi_R2*Z*RZ`xkoA6C5GkRvI#ba~Yli^kzN6uVFKcD;OJ+{<^Q+9h1O9y0o_8uv zs2`>Oe`<_d^_^aGyQKca;@ z$(}D~F%c}f;2eDa;4UNu5B6NIAR7OeT*6cl9I5wKBX}T`)ipJ{0K1HF@t8ZhN^RTt z)#q7kEEz2ImXQM58jl~4vK}BJYvr;-3(blM(Wvyr_PTQxY<=?yf=`z{^YGDO_P{Y@ zqLo1;Kd)5@58nf44XevxmI_i{TcAWl;Yt0xgoISqi*0;9ms=~8qr=0QFZ7>Sz*9O% z03v{PQUaSG0`Yp`la3A5%`-l(u7&--eGazS z$BspaRACe&dAI{=@RY+5f0U6eBuogSz16?8{%FbO`rI`YdVnL2d6wt*cKd7r{T2`d zGn8G0He8k|FX*p2q8X8bDkYbKs^X4izESL_^5*A>7CFQvK)~{V)h)~pwB_VlGwOhR zM!w6Wd3hQq#?g_5@9$4IMrfix*=_;T&S^~nD0%Kc` z`6{1UOLt~1PN_yI*kQ$7PgXFMSI{yStV$X;Hou%$(!dfH*u89MXb6kH%)Zp}?$)(} z!AD4u=vRof?PGT&t{dBhZA?^1NIQFV7sFs%mX2BqJuB?O6P=Ajmet4OXuaq~^?$Bc z!^yh_1_qJ`%}5A`ugx|kjr~X6r%$rCLy+O%iGQK-lq=>q&T`T14bR}!E%=8Llg=oY zsU1{wbf@V0&=MlrLt`!4T3S?uoc;z4?MsfmgqqL%<-!b5u)x1=;8P^#N}h&?r#NlJ zo|)&LhQTk%wHKbXZsy0xq^l1DS5Y^rG}^;UpKZ2j8Cm?=0x_4EGJ{6dnvpqC)>2@l zm0+u_)t-}?nb~!AcyD!gC*7+vB}8Q^ih8x_go&9i(LN(73G~}MIPwUxBvYQ`J*JhTFJ^6>C9VS{iGnKF5==h%Ee5MbDvUs? zRdmF~w-YJAL`DV%9!{B@(9;j{omkVipoXRCQMDxW@Tp=1SuF|0^}3UR@UMi z4>$L8hmkLFW|P*TBNXE+@4DbjJw~>}?+iYJ=&>7jws|jh-h#?fkme27`e=a(;X#Q3 z%2PM&XApkrREDG&k}{7QC#n3=e2)+o4aZr2{tu|sEfDo9mW-0bcX`i!6C^WhX@La3 zQQEr}71!$43#mbLefqxCBy3DlFc^(Y9~YnXnL0vS>uN$2n8-x3aW2ix{Se+=V*RW5 zpr)p#<>ijAUnzi?!q?-#kcrr6!5U!bdt|01{`f6C$H?Y&aMR>j9E$N@i^M=sxw}I_ zVqdnNAkk98U z9PiWB(@Q~WgvcQ{AhEM^*Dhih6H-L6&)Ru(UOxom3Ee1c0dt?zEz*n@`EbOa&_qHQ zR_KmazV4G==~&^tqWYf=g9D2E&gI;%U0p|!7dW_|Q6vb$!tc&^w@67z--9x>0KV53 zVz&ojOIXYwb=F9c+LvWib)3-HFf(U-LDs)sLD&mmtz>4;pFdwlyztx@BT4Leo4_W~E&U{JCMJ@tjitm|L{F;f(9-!VMNt{oy`y3 zXouv>ivO}i(7H2QO4iF}w<@{_t_r3*_@?uFg!K`s7KB((wo?(r z1OP<6GoKyYc|*oIR-0>G=FNz`&Tv*MLwzri3yEA?Tl+#6Rr}v(^S^?wR0_JaQiMG& z|NIBX4`_O=85;6Q+{oHESLoQHuRS|$d{10QN@_RS^9z-W1aD54WczXWU}!rAI2|Ut&SAeu50+vq1N@%_ z?!jF`QGzy)lvmr;2tTf^T?-f>@v6uAd6&HH>C2d1>>mE3-FhaJsTJ4kf71egg&$7SXA} z_;rLm)sZ9qNXtNy)y_K!Ci&9J+^@Q>ZV&G0NF*PoHO)>&#S>Xk&N9yh%cKRqhnOg| z?yEFV7XglH`No@pRF$kPd7*6&XDMmMene&nhyQ3txNspEDQR;t!v&Xrzv%WI6=+Ti zc}p&bE-D4WFJPI8U7h?E(_NOW{=CSG-6f;un}x8K`{ItcM|z~fblo3le!pJp+}Tm9CeD|j^6Vv(n4nhX;Z|p3?bi5LHpv&BdM_) zOE4yrUBugvxss@#>N0FyoShX%oFzX)?w0VI@OPqAR$38 z*bx@ppvgyl+~X0}j0`!;@U9SeQ&$Hk0}z~ieJwFZg`q^2Vb75WQbI?Ac`MT@wdslw zlGgQ;#p;Yx@@F=nk!H;7bBFPl{BaofO}KE4&Ak`zF(F8~@0d_ChaDx43}#r7mgM1F zT3qyp0j_C-cBQ3$1TjRIf3VM;odq;=4YTK`Vvm0rbF2(Sf*4W%V8)fD+^I4MPsF&7 zz+H6+^7bRmCjcF3InIQP+v`h<67w?x^#x3|LW?V4?$zLpDZ&l6nPFwOE?;U=t_np} zf3)Zp+xwZ@o7G+&YYmH}(VK6TeE-2IB^T2Fhi5G4+>R-;QAn z=h7Kxn5rBGrTJ1W5gYBV6=ITSD%+aZFdcHn_P-ZqBoH^{uVZ=g`^-a(8owd~g;hR;*FM9^;{tQtv4_Z~h znEW@nseu6%j(2vMDUhPV-btljwPcIAck_Yw;shNc>r@S2zkX%uR*?qXBwP}|Gh&qX zqu*JxySNl!0jVrk0hBll`nOL)S5n}6XVEZO?MEJVzaB@52uK^H5gBDk=g=2v-xPci zwrts~0aVS5+aMV=Ek=ein|A~QI?9f(y)=+f+5sTp_ZFe~;O@{L0L15V)ZX1xtnYi- zI-!=)nq#1+j%|E4xVuGi(9{}?DqPXyzGjWlhb$zBh=Ku zr)sxJ(!}3mCKoDgI&v6bBtIr5E1bnL;B*(ePR0)LD0xXg0SI#Gr{8%30h&_>-RsK3O9n9e)V8lX&BDv6qV?Myf03;-t_@)#b4|8- zRK!M%aJNPN;X+JIQ6D*rkJJ>oIH+5-25bTjM6H8 z3{RcG97q@dOJ*H%k!03-maYla#W7+! zj16J>6x8e+zQFcj_iM(!=0xM3Apiw51Y2mXg34K4W;TW2BU%3$`0FD~N|r>-s&a86 z)7LtBFR=ARvbw0B;{OcqwPk9By;GCj_^ZKVW%QJ=AP?Cl5p%n{_aEcv@`u*48~Cfp za8=>jHEDhQ`0qN^@S~`o=YF|$om91Nf!Clkt-{*nHq4!{vkB&NaSgu43Wv=aB2<=r zF=*RvQ=f$V#L4QUPfj*RHX*lBsW_^Sz3*Nr1Mx=tJC@*bNj{%prRbaZ*Pq-BSc7ey{+mI528&6Mn8mw6-G) zpWP(p_$RRyFQItx{-dzjyvs2IszsMyb?JOt-zW%&3_&F&8#%?VN!C$6VUM%~E1~AY z*o)Kj(ZZ}R+@a=-oMU4hxd-40`@Ovv7B$P>HQXCyf{`B>l^JbHc<$!ypD+h7TE0Dq z-`?I{Io32qZqteH{q(R=t&*i=Ee0mW%N0-WOz-#(T|V13sL2eF3~>B$LmzGnDD*d+ zx?Sy00t3a78FK&|-=3(|{B>oj#=TU()WeMSOYApmvizaRu2u2h~Uf1QiS zgad<5i|fRXpNkd9>=W6x+J%>u9Vl0k`Uy&Lg=FKGSD@TUnvM0|;aeytLTrO^+_Tp- zwLqMR7Xzf6pT!a;G>nfBi3C`)7z>{UfPVp{GWjQ@M)9ZnGe2x};puTYW7O8(!6Eb&kB!n_4U4Y)Qv@x5^+C+_ z6(dcB4uCNbgNf|_>h_@+NRzzf+ICjZhd0-0~PqXLi@#>Q<+x2PLgJx!C z!v3hvkUWh)7|vekhih}d7I`Gz2`JN#>^<^Y9N9#82yX`9h&UzRC`e@!t+lh)Fg^qC zIpNsG3L`4m=ERA|b7kvAwuGx^T>8tOO3?YHA~Iz$|d7;^)f zef7?dD;${+VCL+QboofmS8ChSTm9U9Y?o@Pt0kt2U;!!ZH3xLN7aSa1-EasEjBuFD z3rClN-*)!)MGxx1H6wg*>|;vS2_ipmrew`PRj~V%-Um5@hTDC+CcVg4yRrSwZEpHv zLXnuz9r6EyF(h-B6Do9wRlwg%hBs~AUC2YcDg!T?*$;Ch&;0Kjge{1d5`hDSXZxB3 zjncXUpZ$dTDc*?<(;BPjhCfOk?8aLs(vOuDBe$RS6>6`UZ^qEIkKaCjCSJ58&{apo z4&dbp*M%j*ZdIxq_AAI8L`OqRt_~!dx*f#)lH}NKl;br z*DiBq(J{ceAED?Yo-jOZN_hl$bC}o`qLE?j5de4C4IU^knR@vlA%S$1vbW znQb0^@jUDmb|SEaWZM7V>j+@JQQ`Ft*g%b$gA)@$XkaZ%s-4=x)3<|Czql|^a4bvd zT*xx4&~=>6z)bgu^Zn=X$Yw8|2p%7K%-ff_PEj(2W|5`2b~?S|#IRTLq@9EsK(!>+ zn5vdy|1qIEh*Uk+%(~qH${}&sf_QmAL8<(v)kIg0aONratO&*sFSuD4YpVgwqNJjF z0PIIheR?#>!{Em#7v-^I`u%o&w2g5Q2X}76HvWX0B?5g2F${t+CsTX2CQY7WL3AGQ zQGC!-wd#P)ldXL#dH|M#buq--b?oq>Bm~8^iHS-?eq=CIjz{qZUQa$UF7FS=BfK_% z{2cE`bC>9?oE4zR5pTyVJ5}~RKmIh{HX=Uu{ngD^yc0A9rcFzCA3l6I6y%3zx2wGy zbGiRb(e2t_O?}?xCWDH`7V?u}P1xhqk+M7+S$oqZj-}t8HwMClE`0@@h#1o*UVDRB zlo))K+DG$uY`d9!Shm%(B9IRpmMjqrNvj|GT-boA91~37D8$RMuK)h{2tG6qMpI|- zK9TI3Y&@1UTEQ!Vk-qJ|iqGP>de@#~s>DUlvH9jxdden$cCCwd{+`u<;%baCO$rGWDeH7JHY%yVp0pE)ervkEo=Kh_$Q*s1#y|Uiln0d+n0M!F<<&6Xd?zt6aY*4N zhHLf_M&}Z3kcDp#`do`b0O}_a84nstu``ZD21CF%t4O_;7}!9BMKN+XGQSM6B{6AP z`aQ8SB+)jbo&CRW)I8WXe z8kn54d3MK{=^R%D#_Gz;%S|nOu|ef-*@%D`vnj;~5HG~iFLA`g#%|B~k-HGe)J}x< z;{6R!fTDl#QO{RIuj;ALPo3g~hQsRNGzO}4+N73ZhXf1KIp0Mh{_BBTR@Q5V7pgk^ z*MrT!eEhf*4u>Yk7o?TS3zu?Wvjy1RBS3R)$0lCJMr25p+3Fgmx-yKaM1b?TbhmMx zKd<@6bZ4!7Z-kMi<`K;-?YqRZBbHaY*#5-+$d z7^Etr@z&RzuzRk6f1d`{{#E97t(L;c{k<){e|0%wn}4S- z-F@IIx#LTfc6CLnU)uYC!~by4sXnUcGW?~cVKnyl_r2dFy=?u&h)kvsQgh6ld_NIVk#JmMFo1pT2}j#T)XbQ#6P#2B}xsgRiwe7vD*9QJ84-Cy9nwpe2*_-+Q~ zr3`Hy{P{U+!Z<-K+)2PEmf7xUdyx<}q3APCZ<1rIQ}>hQ4Bhw|b2jQlS_T7j>`1O!D(a~pCDD0n|HXPyRNE}f&R-?6Y`>OJ+x17~RAb*3 z6nuixAyLOL#)&B8;#A)gU}%}fXS&T_i&Ycr@WdrE)zZe%A^`qJ5qQ+kbbVZkGjjvE zxMY@c)Y(BGCjn*k@2KmOSLr%hqL<9JxQ53$g=>FZS$_WgxE0au8_gE4Lo&a8y%HdJrnut!rf8*x)_(2nJuh6;(CMXBbDT+K zDOcg5tnkd_)a@_MR%UGk`Q5$s?QcX+zM>A z^PQALjBJWeavG`mf3jx@8JjsR$r8cAM$k+(7u2RXE9CF_ zaCSwfo_mU06`Gize<_cnvMxh#K6{Y^Qf?7oB^Wg)ndP=p!pOy6U%2U7hx#iRy2Nu1_@aOWw z&L<{BBe;){B>wMJ*7CkiM0i6fj`k@+7Klg9Ppm-ta;d7$OpioHdO6YzoV z`}%SQCMj|&d~voy>h3f>QX4w7PUy@R(CTVZo3M(wg5=(yY4T9SbS=P!#b|`f%FAy> z8c#G_udV~1QRoxB^KNqB(bj2u09V96CPmE@WC;P|!ZjC~)7Jr3mBw8Ky%EH)E6^QP zP~@7fDcJlDh&xT&WX9tZ*jhTXmS|RYWPN#fnEPChuF7OMvQTEx4T0RCBCYcEMLL&l zsA628u#odhEZ%LarrtvqZ9trS+>TQ}1+ZC$>l0q~C-eakhLCut4_s0RwPq>~N=bKkvv)Pus(o7dGAHf40t|aYw(D@5<&c zpcm1XJy29B3^oXMJtE?2nUQFFgMKScUjU~Qqy46CXRWu5XrXHwIYKWR1Bnz1^ucef zdsd^mz6$f2P!PbNrxryor4?x`JZHW%P_jP$s+0iFEu9yDhN<@%k5qu>$Yx(yq6#8&L} z^KChL;)N{S04o#fl!BjPLH+7DJ)PCYGQo^%y0_rYI@^b%g_Wxphu@D6h6f_nNH(^ zk`vwAGKPD)x_0dIK6ZQUE}xon&FTBlLX}-AjKr8(mzVcEmgmZhF9 z+7)@YOKS>`JrTbNu916(JLn3#>DD139 zeQ05|n#1xH>I#D@ug*q9Mv{D=J})p3obPC@>2#uEe)0PC8cqeM?@gIy{iqx<>+|Kil~BWzJmICIZk5hDFG zOPKgaQ zA9Q$1T(LSSf=K?OEVfE0y?xc%E2BJn^uY`S#fG?8nn@bz3L@hdezH@c)a=Bm?{@7^ z@3pKOI5k^?M1PCsjO50aP8HkTkr#3&bl@rFr<0-JzacKXms>lW$3G=yA39~*xNA8r zR5hD&j<}qt@y|3d*@G2Gs;5{cf=65WDMTALiYUlB?j78$fsQB`(MGPUMcf8B#`UsX zcG(Nv#?3m`-$}wImL6`B!}xNZw82BYCl8VY^5}W58+J@BKX$xJK3gFzI3A7wL%;VY zz6KA6LnkHg)2B~dc_#Un`dhNt(dS+3+qZw7F8dzyMy|kbAK`Zb-ix4;>_TtF2nq-e zmu+}?c{%bLz`HgPe$}~Fo$!6Tw{y|`L=@tr*lz8&AmQ*RlJnMCAfqR zHlS2J-*RCleu`;tl37DdlnVp=O_+R{ttCFKB@bwbh2Vs^xmH=~7lg8Z_eRxiN%$M4v`Lc zdKB7X2Fi7)5g2riB}Cm1jO)aqnfN*g2OrQBWqW`mN#<4PA_!U-Ar#rqmJ%vcTp_DT z_~Pc$;q7+QI55Z=uFbM$O_=Je%F@%w@vEbc8V{8{J!NQk6-V*8PbX%s4UEC5gd0Z? zPwLJHFfvI4;FwVTA79Jp=rZBQQ}97Xee1Sukl;kCihG=>c2y0N!a=;6?eOw-f(x*C z!Pe>8(+3z`Y}eBG0+c2yEwd=edaR#6ftuX@&QQoy4aBCm!TG#VSPQdK%g)hI>AFF6 zw6htWn*er2X>hynN7+==(&$v=-pZN@5SfJv>h9jk9ozTq+jE;ofewRT-Hq(g-o?o* z@I66SS66QD8@Ma1$M!$LOB)cbI;fE-GaW%y54SWO{V9h z=*XdrF@#Jno)v4)Fs7Oxv}R+tFN4P^7aOjFYLc2!61*Vc!XEp%PVb)~-*CFox#5Oa zP+vlJw|>W#%Imvax3bZJ_R33Pb`NbQ3<%#5Y70B^zC}0WlW6+G>V3-x!8rDnJ9zN6 zr_CW}93i=2rpb1o@|1vfz|8GNGi77#7H%gr`I-7x%+5?w+&|e7%A9frLfg4U4X9r# z5~g-1eLeZ--WMC*+pFRs?bdnwRO!Ua({dGeSNA&T)0UaSDM1MshAfL_5Up?gsP~W< zbzQB?wGBPV zSy{Oh^Wh!(;F~H_b8~ZO%QpbHLw>H}I1BSL;&fJjjTxF;g5HVxW1BiYF49lUu7tI!!ns92Z_%`RLZ zeOXQOHx4MBqMm~%r9k00e6QRDXBVXdF&SNz%kf(}3W|%X&-s}}Q{SppCg{L{ zN=`u!{q#_?pzUZ1Z5NY|#uE?O&b*=iAiFaMH{$YpbQIxwHAOpftnBRkb%w;GqGjKK zd2pVY;FqglZ*qctT%lyGQ&%ZloXHQLi6+3REWP$rBvc$lxSY6Jm%BEu5@(88aeFko zfzfeK6l1Z)S1=B!?4=M(H7u?e>+w4brurm^A;Q54G(rhy(DSOnI}u@$vu1ea7ZUDt^T*XSVRlDX*2EL-zw63~mR zds%(GmheuZiRy2^DO z0^_4{LcV^ibv>~kPo)7V@-kISDu!Rs5euv<`+RIz9Sup}^bHcQ{UJ6J7p9TNpN4cS zbO?2D7BnCDDXJ^kJVJ8e=ci` zk_pFLU}WS|%q%#^H%{c! z4SWBD9GYk1=-Ee?0AcY|d}F0vC9z-@NI=8cHFB*rPsV7Xn6R%F0HK%K;)u0^n*`1Xh*mAV1FF`t~z22pfofQx^dI=;^c@sBdl*pb=lGi6qvbmSy*O~ z3XctQ7#f42k3K_7brlQoeDyKqXGJd3G4YwBHo1xSin4ZEr zSh^5WR|qfcuOg2HRmNl>_vo>4ja?f3VHCi7Hkxv$CAi1ar7NCl(^E&KpYKY{}>s4hLW@7ZWKfuAKu-*R16LC$;iw88| zFzB70y)_w??0c|4LwEsBBYFK4-2^rgNoCZU!3_?oFA;H>x};o(Lm0t>GiB=HOdanX z1sGvKvPrqyXAH6saWENz%%3~X(9m$3eTA4EzlthAz||$qt?lg``z9mfhg^@(I~(Nm z0;0Hzg;AmAh?|6ZAlrdLCkOd}QQyJ%GV;{DSj?+2m_a}}V2flGwpA$VP;RPr%Tv-2nH%@>-Yy$nsu!@-al;U_JjN4q75XO_BEI@NGJ$m=rXE$ArsO2R#<`>ng5T4>kPGcx+5 zrx<`n;6QCSdGGSBUAyFNy5Yg0C$*pM3vOpcPAgD-qh*+Jm_o*%k2to zyM?NnZE>L%?X=k(o4{m>^%F^Fu@JukJ@xBS_1vt+3W z*ldb}U&ysI6NIXTo*^k|X*`V=A%nXMpRu;a6f62N9*BP%o0@FyXgTLvem_(Xs$9V% z@ituqQopGScQw(lXus=@yh@d4U&f1tix=O1qYiQev2_1r-d4ny<)gE(2cPW!xbemZ z7++UJ9%rxMevS<51899CE~>rP?;sF+D8q5n8d{VPXygsS`!#P)vD}S)a+k7aakJ+j z>trF782oLN>4ReV@(81IGe?(AI`rMW>!I#fesJ1p+k61p!^F!Eq^t88cM^Soe53iP zr`C+mteYsmK4n+3!7;rM;%Gt;3>_luCpdeVr7Gb&))oN)VP zN4xG^(=!xup0LJblv{-W(z#M%yx$-tu1K*mDR@NJf}tyt6&8=TBlY6Pj}>Un5+f2U z?bs$iHLr4XWvsw8^Jm!f>W2$zSHoJvhBD>9VFJrQ|P>xY|6^v1o#kAtKHW&x|O;dEQu zc~{5HB51!?ClTrfSW5t3%@`(rLuhevi3{)Z;LXQ5PQ;vK1)cX{a)S*vF(U{SpgHM2 zn+`z?0$qcX8xBUYq#ud|&FVSz6)--+Ikff;&d89!wLKp>k`=bevv1|PV>1E)SOSc! z{0FR6C&2_K!uz`oSzv>8vwnJb!o?L_{B|2Qz(l?Zv*ZvlZLkM@MLOVOJQqfLtYnJ} zeE9HCGnsQE$3(N*Hbfnnizb#5nAaq%i0OWw->jD%NMzcL1qBKnO-&wa*01MY?^9)Cj%lGDT*Jr`L#ESqD(OYtW} zXW~~Iko(rpP)J-N#H_iHB8Mx% z!~*H7dLU5(mE?>gRVv_&^hE`F8Kpu!)k>Klh~$pGmV&+T-I`a=Sa3VGhbE%P(R^lB zqoMB{ZjqzlLiB6C9l?jZB;#qa=yq{JO~93CwLcd!&>Zv+b^^jO3g;a(ni zTOy~;LCP%WiiXvD=oKQA9g!Iygmh?4A&_fA-jBA~)`D7_lPMmWvjw0PCAAU#1DAjR zuq?13PdJ`QK=q*5yha?Jl?Bs1X^6-xjvqgc$GaHxzds!uQtMyNr)Y4T#9$<_YPeHh zqI{nGo+hpB{1N9KW)f`I>gYs