From 876adfffeef164e42ec68fb1663b40d6ef48df02 Mon Sep 17 00:00:00 2001 From: Steve Martell Date: Thu, 15 Jan 2015 12:23:22 -0800 Subject: [PATCH] Stash --- Rsrc/R/plot-sizetransition.R | 19 +++++++++++++++++++ docs/manuscript/GLBAM.pdf | Bin 226245 -> 225914 bytes docs/manuscript/GLBAM.tex | 5 +++-- examples/bbrkc/bbrkc.ctl | 12 ++++++------ src/gmacs.tpl | 22 +++++++++++++--------- 5 files changed, 41 insertions(+), 17 deletions(-) diff --git a/Rsrc/R/plot-sizetransition.R b/Rsrc/R/plot-sizetransition.R index d9e6f354..32058ff9 100644 --- a/Rsrc/R/plot-sizetransition.R +++ b/Rsrc/R/plot-sizetransition.R @@ -20,3 +20,22 @@ plot_sizetransition <- function(replist){ pSizeTransition <- p return(pSizeTransition) } + +# THis routine includs the molting probability +plot_growthtransition <- function(replist){ + A <- replist + df <- data.frame(stm = A$size_transition %*% A$P) + colnames(df) <- A$mid_points + nrow <- dim(A$size_transition)[2] + df$sex <- c(rep(1,length=nrow),rep(2,length=nrow)) + df$col <- A$mid_points + mdf <- melt(df,id=c("sex","col")) + + p <- ggplot(mdf) + p <- p + geom_point(aes(variable,col,size=value),alpha=0.4,col="red") + p <- p + scale_size_area(max_size=10) + p <- p + labs(x="Post-molt carapace width",y="Pre-molt carapace width",size="Probability") + p <- p + facet_wrap(~sex) + pGrowthTransition <- p + return(pGrowthTransition) +} diff --git a/docs/manuscript/GLBAM.pdf b/docs/manuscript/GLBAM.pdf index 236dfe532762e07e389d0af778449659bc717a17..ff733abb8956901109d05632c53cea4f036f5bb1 100644 GIT binary patch delta 23074 zcmV(=MDPg4X{WC0x>p|VFW0DjaNyt+cpxu`&aNtRhbM4?sIcy64#{0m&v;& zIY1PigfbTe0PRwK{q#nP@{=ob5rIad(filRqH~l*AKztv^ZoJV-48$Kb)>Siu2mJi zyhdG}7MY6bs!JPHM=z)7xA?;$0a|=ElW}^mf0+{#JS-C*+|q}|twik)dF+jU_oJB( zs*VGb$N7-J?qmkuzs~OP$ICBK4iD8ul&CW8GI{TxZt+3)Mo%7jK*c7?%d{wKg(|jb zo3&9=RB72MzVV->w*~~m3!k(LW*A6#lILlzPrbVlNqNFU-Rp$>!aD?u|_wK9x{Uoo!oCfsDQl(m!@Nv|= zKO|*S#8$ow)QS98bi6E^g!g06(8afZrm~U?h+>G8;ravc1UoM7uzkLYU zNBtT1#`AW;^f@GXEH~y?{%St{YEN?T@1Z%a`I0i_h$GSoKaH zdG(mMo{1Dl_gw3fOtll_p$%SkynxsARUhqO*XqK&t=0Aq@e^}*{OGi+9ArYliY!h zm@EC^2M{l(i#Yy?5!7v!ik*RMwP81r z0Fj;S?AmX@QHNH=h?!-B-e;6m1!x)QZG{7+gZ0YJ!lC)##L<9rq1v2bvfdUfgvY-T z!zb1{KA>K$WYD!e6T#`pGm9SPUI4_$XcDaQ!9s=6NcRK|y|=JvyZE>Wga=j{nKdtJ zNl+($;QR_*h0Zq}kxIh{bq|rqDBWPUI)nGjiBBKw*c!oJ@2pW8`bNrbje zkm!&|rIg*yVC>5X3kWZPH-xpo#7YeV?>wY`EENiH@B{!vJptMc!Y@S%=TK-N4)oEf>=7fGe&h&j*q^I5Jb4 z!ZFzn`Z&-KST}-1!1;(o4V#bO6^VdnJ<_~ouz^Qcya5^+|3vi_edy2C9dKa7fCJip z3E%*U07u$({T6UU0fz5T+ykLG)^Sb!lmWil6v8V!RCBQw zHKQEQJxsVRZiOY||AB>2lBzk1mOH}Q4@Nus6N2{uLlK?n-yU&YWHHQV0ejfY?nG^* zsATMbc5<=bhatGw%an(H>+WFP8*ZMNR41OHqKXw&9l(P-L(TssGa@ojIN+;)`GdhW zD4ktGb1{$y?GU(eX}LF=GeVWT@{kh3BVf@Lt7z~u{F{De5>2-_fh@=5%|s!+99GUi zBk}h?jCe!(aBuNHOwYvo95L;>3dz6MXGM5Uvks4`I68yz>$p>G)AM`Y&PPfVg!y%@ zyXToI@-%BoDA77Cpb*Kd{DJ2QrXzv zbZFINF49j?Fr0>PqaLPC7yCClqP%~Ql+4(+J()$$mQpKG*;q9;ww7^BCXtB?inZ8Z z&5_gGtXB^O!85bX9rCl=g07`_jkZ45V1U>NK>Ku=O>7yVHNFCwHYSgWT(3WKl{xd? zgtS&a--hMow$mI3fmO4ogzu`x?;gYI)>ns$ao%4SpAEsZxZt^TGvslPa}9s9<+jUT zZ#Sz?IrF0iC6JlsfXB?{|5|f}?wRfD8F{pB$J!%CR zxVBZ_lFWaiBiu!c;iHRb0tAywX&?MHN-h=vU-HUCc+puG zt_Pim)(v z@QK3@kyZo)2Is_DU%8&6`DmPM_NL5Rnxl_;i^!e30u&__aMQ_KkqdgJd5D;L7%rX@ ziW#O|O$A%UM<)q#3?>&`wO+vb=~I#_$gM5GUa1P8Qa;m9^AkrxWQo4sz*4+QxfTO8vk)xkGh;X z(dDA~2$!QZ#{xFgM!<|I*SMl>0jPv`?c9txUpH*4N%A;QmyWqxU+o??huigC;0cGa z@wNG5y;yDbhu!AG!yw@TkBrEuVZs&IlpiKsg)pm`Ad{o}%42sS%XG4jJI+LJm6d8Mewyf@kN{RCs<@j4-}i{?zP;)X*^7wC3Jdg$ldMN77+bMAVa2`Z!b}6SdEN?2}HFLThc)G z+0cxk(S>G1OdJ^m8J(}25>Lte;F%CGo*zQOfeGjhh*ena0Hy}T$6EOwU=YHs$XO`> zMQL;=I zmLo8CQPd8mW*vY8sjUFfQlV8QsQ%+c_jzZ@nwgV|ggjf*d*;xJZhC7>kG9vIfk1K; zfn(koKR4s|LMe5EbJAx|ZjgivnKJM-VNBUeH0G5pyq)QLG$4%L!Nirqn zk;h@i8e)G+T$J2qur=i|bE$BWs7iNxO;@F_6IJPJ!rdIRR)3?KSx>5&;e_tvyP@2U z5$4=cxqa;T=;+R9D!n?arfd7XX&I_zr#Fe+eM8q%JnlK%(K;Gs62;VTg9WQrA!l`I zty(k)LBLxiAJX`N=1w?{R$(1bRgg$Lh2z3h(dmB&j{5Lh?&4s*puO%QVYLg=>%g1w za=Ujh_|;?LsAj-dj)_>ec)F%v3&S(wpCD`FUDE+^o!HuIl2CB4v4R5)2HB-!LvgG6 zaciCvxA9qFGHxwJ8PjnK)8p2Tgstp_En=FlPw07rs=^=(ZbS6>2eebT%=eN~P;R8?LRy z#~jI1mY3k-$8>|rJ2*!9<_=D)*p7uS!dZVSGSwKFA4x~B7EXs+BC_<9$nr;d@I@lB zw0joR5s{@|B(lUb3tmRnT+Eng-AD;n z+Po7T`0U(KHM6VxyVa(+@0x;5=@An#u?@l|Zljyfw3ik^z1r#vVq+WK%XH zmnS(e*4nQP|CXH23u?Bwb&a-1F(S02%Od(ebQGbKa9#pt;1&Ri;l}1SsCW%gQj-Fm z^klYxxN#N$Z$vp9Bvqu`HpQKJtlR_5Odk@}#{+#o!YGEecm* zTCe_jT#u`BIwKh8IS8rhuy4%A#Y;Q6dI=}*XqWbK6KVHId--F+%V}D|ou|*ClU}YS zdHJ#8QeHmklw9J-aCfGEBS=8o$uW^*O)7m}ER^U{{I(#!jQ77q`>@y*hsA$(%FhiX z_Xq-Z5N~u?vz-7q!fV`?a@Y1`>s*W2F8B25M7Yv8+RN;aNeSfTStL^6r(J-L=e(e= z;)YH!A+3A!0`mHV*jJrCuqd~=P05>MmtOZn`_@_eQ8)RjEB!GME}V-WLxPi+{82sN z^{`P>4zZ7L0v2MOzyTBEfC+yM3<>!K(NmIi4+m6pqxn3Mc#RZ+A36JMrw<*a{Hh(e zA?*c9Q_+2~U%4`5hV(K{X&3aad!;8g7IZ}8I~m1=fW&}9-3I$UfIDx^U_SAHkcN8> zn0l-o4h6@o3+GR@LB5pXx(8JXK_j=NHZFEdOP>r|w2rdO>B+w{>%1I`wGGZV@!bi( z-3vEfI7Ul#M+CaRk^4VOLW-BOCoC5#0Wy(lZ)$r-*0wT_F8tlYipxRuLf*!wL2s~&J2g&a4B`c z6{p0G0{%EUyHrVzisKORL5{Vk!y_SP$Jt5J>bSbd1vi6>LTM)CVzh_RMNPqpqOAD9 z^JITQj51XxD03(-p)B^HPD4iJIb<0iWl2$jVp z0ERG{EV!>O84gV@X~OgDQVudxbXfob&{Hh|LM{^q36V@qBPc673^hwZnO%+M+!-=(Glxb4W9Y=3EN~?Q#)8C6w7Jxm=n>}Lc!2MNyIZEG=`?2 zCV+-UD3Gg!@xsooN`S4Q*9xoQ0v#y2rl9D7DlQVbfYKzQHzdbau+A1wbpH1Vsd}xIP`XSF=yqL}3|6_4}vHHoqc;Rr}70s;= z|M_L<*v-FLU#>Ure_5O@-Td3-zc0>~@7}z^<@%qCid3C?=K8NFU%qq)wPk-p3kS5d ztAAfxxO}y~TA2?`b2@k(povqK8R*#%%*#1wg=KQ-)+vV!qo&31C{J3<@ zXY=6Y%h~+xW_|X4xka1vx8J^Yr!(CA^6meSo)2kjN?+P<>r1h}+w*|u`Tzel|9=nm zf1dvb+FwT5Uu|?D#h&?hx2}JmxxV1ti^q7!J$1%7M<3;ZvEhNH9_wa%O^uEW>@^s3 z-nqyq&)9CRo*m;0cMS5Z0lCDpV>~;?Eyrl8AtilwjHw;t5SsBoQ~%4sM?A4|9dfn< z91h;)=5N*?)w#~9rbidf3ri)$0uUD(})&?uobb0OnZFkLMOopo;w$1lH zU#xzfK|`DK<)+i!{EzcL&A&OJ=8GXY!`LW@m(Xn7c-w0uB2S7kW>490eC_rI{fF)1 zhs&ksuWY~_h@L9G|J#4xT_%^A-0JG`^7Pk4Fis9bB9Z{3rXew^)fdgT#BR-dX~o|z zw%g@q<+%L{oU2K=G)#+^h8un_vloC4cj7!-`?yg}EDZ6U`th;MLYX&E9OJ74Ry?qdOwKe zzpfS*$$Xd+NhX9YEGtd!Wt-x@*oJ>8zUtWNco+g!@)6-lVGpu*9o2I#_RhfvEyuA3 zMT}SM#RvIh8}xrUda-vEQ9WkBkUgf4)2@H{7_N&;J1%E0$**O=Id)vm-gob~94_*X*n>sogLqv@`@F8i$MCuv1r{CD zabLt1E(9XbD0BDUOwT>t_VszAH|RFFQV!8!r8_#TPt#!?&Qkg!OCiNNX;V4-tV07c~#=u_E z5C(7oOVVHkp;4mXnyJK5qFC4&ZazYkp)XCN05F(!#sW1)l+Y^WF`|S8Hz0)>BZ}a! zC4JTy(TmG=wny&7j@xk=scY3cvHOYfg-G`UjGfPUQudu%gm$#3Ex9YDD`Qt`SEjDa zU0Hv+vK^6~Cu+-S_kG}LDo5#Opr{Qd1r8-L>c)MR_%&x!j7(q zEnkC-m8MK=EU^lfU>PR@7efACN zdy>Fa?9?UJ#%#y;kZ;)mKFiKiIih)pZIMQ*SwK{W3<8y0ET_8D5!h$#AnvktKi5%)j|zNxEm-_9n%D>e3wf&e9yisHK_Om*{=X z?QlJ$-M$PKZ=~4vbE^oW#?OCAe)^jz{f(2p8KV7`CmyEmdRN`~QR?ndA?_#;cNC~Q z@`Z=^BTU z_@Mk5G?J0t+iN?Jhpf$^_GU4)UE7*Rub%$}2gyaP3T19&b98cLVQqs~V24*=0k>CR z0@Ml*H8lz^Ol59obZ9alIWRDn5kvtM1u-x%HZ+&vkpU=wwgpfe*cL5{Lm)Uo5~Oi= zcXxMp8fc^&cXtgC90EZ?@Zj$5E`bg%2_a~3e`IFv%-s8b)qB-dbg!*voxS!csL9nd znZ>Oj7C;#Y*o~Q;g-rk;sjR`y3t(g8WMO0DK&7VE0=YQ=|1zUe>i}I`K@hOOzXBv( zfaY$`HYsy|x92!z2pFK??f_ut1hDf6u=5JAu>m;P*!cew2yqbrNSS+ptN_X^00jsb z=!!}$332jv0omBPJ-7MKPXL`IJ%F8`pO@*6bbz=c&;?{^4hATjyV(LApF3KbI{-8x zmLQ;;_kWe36S8%4a}r=>_4M>)F?V!jfwvM)X*b3+Zcy12Rlve_%Isw6d1}ps; zzy$c~ZUF2o?0=X0EBkjL5crQ|b4yEzqmwz<8w9oiSc4pZ096?!7B??9CV)BE>UW^I zgDd2JIp5sF9OPhb@f`5S$;|;W;_3kNXM%sxbG38g#Im0R*_YxC8xt|Lyp1gv!niumV}S0W5$vATa7b)t}8k>p$`H@LfP&07JHC?Xd&c zet-V{W&A8JD+t)Z`=9VX=F6%kE~}%b%J?VazkL!C5HEl)Gbbm2nS+}Rz|X}E;N{|f z2l)SYmB0D_hx;E^HFMBkXJh-Ppghqzhla;rs-XL8W9R|@-9ZKNEM*{o z?jPkFv2n9mKL26=f7bsWDgS@6|6S$(ru_dskc_*7!yj(CKm7k6x49$8!TYa(XC=G4 zJ#U3F&c^mXzUO_jw10k2T%Tw4hYR?;PXFCf8f*!%`n`7?+&lnt z7Z-DHRJLb%ad2}3eA%BD)C%bJhq3{zEMSP+a|+;@ogg7ksYl`Y*@>V0HT!eCFu>59D}0 zx6dE&zi!#m-Nohk_4#8Do=5P1AN0JehYVh3BYu zbu;LheV1K2+;1_EQ|YVThpf9?il>hBVy$dT(_M%`$sT-vw$vj>wkE5$-TB>{CTYy< zwxKRf;}1@K{3%{HK=Oi^SxfAj--ENCPOv?E3v8DHb)K_3ABI{1?zLxszpPi?K-t&v zh{av?Z&f^sX!m8a%<K~kYjuh zjk=^jPG_s)JzdHc#G9S@aGK*u1rEYY9I-Ai-n zV^8OaO@Z`4v`ZgfgM>b|KF0pK+MH;rEVLjG3-PO3CT}O^YIe_mw{#J+SW}c=?QKLj z@q!qRAwt)2i!xfuTX$1BMJ38|aGR7nUuL`wWIIwUPvgk7s`WP0ul#N7o%wU*a0m|2 z@slomD<4?JeWR{Lt|K(U4{z|&dTF+G7xpi>g7+>NytZs-+t;}w8%X3m)11&p-wEYP z@hzU{bLtT6d1MWL&3t~4CuO-0zlTh(Mrqd?E6emXGTpE554Gxc#)Y0TY1J)|Sq!0= zt+9-1StGsSxy7clx+Rv^;D-?tqAy##FMMPQBktCWOTG%KXZ&DeJ_*%0qDU4xag+_p zp${aU+pNQUUD0bQ4ugz4;8i6JA~9l-wmJAUTYvODdo)^q&Gow=yrWG`Rr~wFBdC|O zW2stq-T;!e#^gRI8XaCD_#Ly z`}H@Q@1z}S2Q(1Zso3U#;9oJqKBdNFmGG+z4rfV5gJ=&$$w|I3vcp%SF*lgP9?Q%k zBk25iI79b;I00hK$g!4KnZ^^w3a{4*$KH0-3%}q(%{yK41f6^#l|wJS$l}0P>4`MI zkhZcuj*N{da3#!HkXl2*#hPOF^)^XWZ^cX5b|l6m3-U~q#gGo%cq!1amxJo}TAr3Hj0A7R^YeEHJNsi!JT6{}5`ImJ=T5euaR-y3h6#P>%Lh}ajr)vogs zKL4114>3&`T6OpFIVMM5>L9|{$`$w!*by>*kZ49sg3a+sQ0-o-i8K3D%+F0O9P#xm z@wA&Uo_&wDSZ*$_HK&u2veXd@#*{Rf;gfv`W?2oh!KC9)NjY}0H!`XsXlQ3zw8k|} zN3(2cV^C-Vzs4f&%PrF{J+Vm(vf>?=M;MrYNrgotZK)-m#k%@hy=#M^+d)Te^@!Bz|dd)>PU94nixkoiaI~DuX)l>ax z&-CQAi(nF*0tiabgWfHs_LM#h&bC6@E}d4U-YD+}d!*cf%T@*#4s;_Ut=gUO>k*lM z=i@4%1Y4Ny-eE)86@g>Pqd`oh;c+GM)g(t%y)yQKFRwI6xEG;v6?;a0iH8t`8|15f zioT(sU50|eP~G=ziepj{TAAf&b9MmFt%cCz^d0+aB;GH5^UV359=qWm8&s{c zCFpD5wPz%)rTJLu(B*a6eAl*yWamGB9fhaty=cIdsJV)n2$#UsAabI#ZCgdc#@O%L zQZ!wUyy#Lszg+gymWo-+EMJIam^fq@SU%~QRr7M}%&m4I-G8+i%QwZUH>ySSinl`} zOoG)N@1^Ip9!Kz1fbOVHYz?3_)a1%m3bg;LzNZWLo=i)y#Bi^Wy_0ldTU@?>Q?UqF zo&B`>4q+|CqD|)lerz73v3cAu_!v+)HogNN$Nxl0`v%FVxT!Q31B&J}xnGf_z|CD= zjb!5k&zl%*BP&-KH6=+G{C0Wn4q6}1!kjNBv_XV2L4+C2*CwT~wrpUrR2L{a??s?W z({=oW3Z#NQf1FENr>Vb#Kv|6xM*8&6^b?e_zY_taG#Qy2y+Q zmsP~4(Smn5_kq^#yF3vRQ#2K57qZ87H?!T_$pW1l=lDo3rdwG^tD;@L9ZoGc8gmAH z7tY8gEKOg=K(Cg*Y*=W2BkWMmu&B8Tb9OW2+5TSt9sWB9b=`Syh6f-$ARC4QbjO~p z1Rr2_s-$pd)By>Qv{SS_Tn)%#Dzl(=Zl*GrxLhrcP!WvHc)%VcUWoDV zGC+LLBNT>pU`jc?8IVb6Y{CjcpwJ~NJh3|~MuntkLzjfe7A$vv4n8*A=#g}=vkF~z zcs1@~>ODzjeP~ubiquoursx(a%~;KcXO#hs4>uw!dvcrV&xAOCOz9xIv4jw;b%onm zv}(+>O_}ATms3po8xtmyVZRHD@%q@UNxRJ9{oxqz(*|E81zfMSab34 zu>O>K@>ropa=veW^_RTAQg=eMreSC0t*v6#`d)Gg7oAc+$x+IFWo$a`2bSUr3IblX z4-V;eI#W@tv$a6^J0`Eu9u?lezOOnNUOzBbn^O|Q7lVggR#h-->>&0C;!%yQ4Vh%K z0MU^tp&Jc0f5TPY?74j6z^YB7mU>TMaewNWC+Cb--lq0{CWkC8!d&3ku>5y2Rp$%+ z(A4G16|-=(Ye9^eEQiTGVU~Wu^_Yx)#P@*R{x2*#MXNrr6JrWKh0fc*rM%tz_C+u8ITMF}xxl31J4-ixd$O>RjoSN<6)8=p-t-3zSOhN`8D2(-8g|-rQBNm-jireG z1%L*;jB%L!TIr*Nn=vVH=UQypaV(9rn=U-yjp_$u@T5w_xbb;>hd~lpkY!4lE8^~- zfxy4QZ@XE_ubzRE`}~IgKtZT6d96RiOYelN)w^zgR!}_+;(bCw}F+4=wjKoYCHz}b1?pO`BpBJ}vG9jPXW6LPph zkbd)s=FPI(1BW$l{aUAR;Wbk82756od~is?GW90u2i+v6C>&olCseBiDzPkzw*m|U zstIg=KNNt_lH{h4-X4`4^*8n#<@*tOA5elES}*mzMRPVv#+Eyk%kJ@rJeaqkCROsl zqFLrQ$ne2U>7t_bIxpiL%4H2Hawa4Nsk#I6OOVRk5D(}Y)MtxuP~XjRH%6>!m7jH6 zPu(=1XIYL*Me@tHAxNOi6VNo#vejY zpg25rIoUfGC2-ij;~{L}7i~GlG<8I2DLP>R6_c~w`c^AMDsm}JK`%VeGdx7-rYg06 z@u>507>QGQc6$v6y2kior4pHU!P@MyjUSZIPGhEqOJ!CbR7CfOMrbWXNlf7z4HrDn zZ!A5v1Sfp#CiJ2jK{gDs_+&?+ubTG2v8_l^pw{-24Esv$dNG5hU^0Xu98d~YD@W_Q z!?#sVZB0#14J#Fy^Jdveq|22*g%M_d6t_2?>r1A5Om+q?b9Q7eFMvOUo5oj=}P(4ojuZ7#Z-ipCI{Oi}|W_9VexKI5Nc& z;d5mp(m6bec)dvfLOQC(tAfz1fMNcQga}E^Hs$V`){F1_VNi@cL_`HSymY+b8KKv* zs@hz3R+a5_R(vB;#Jfgno1mVia?UM*bhXCxVU>K|^eZ?HoAR=ON{w3X8VBdJ*R5pz z7jeg7q6_n1<~QMdA+wvLTaPGz0xc_!ku)!sGQy|b_dlbPJDG~+_E>u81u?yt2(XI= zt>uvw4GDdUTEt}{AQ>}0)%MFg7m=t8#X8iK4)t6Q5&nuuW%snEeMHOV{TQ?(h1i$x z5RYVMmTydgkD=$K)rdOYcI9y2g8=(vAOgC?AEi?frkjBEeD4br$3Ay|ZvILCQ>eQd zs3L%SKUFD_%CF~!^KL86TX_9^xFaEd;oid_Sy~kysWp1UW)oJP@CT-K+>bN88IG;b ztj_$6_VpXJ2f#OUR-QEsqUx3rU8m%@)N`Q?;_od+>;}*~LhtU(uW)$UD4~cAoRIpr$}lg}VRC-b0Z#pN#3>pOUi_4QhdpDWJYR@W(%4}A zR2+$Fm&CH#a2$5Ki00Ino#3^;Z=@ouU&Hn;}C zD`7KjU^SM$#GJ4X6<6+kQL_|0fN-yFj~)IpkY6c%0?d4ze(VKd%ci(*w3_<0)dlzpL%i`K%afxS*{q zcu}Ut!~SBsr7(ubBR(v3Rk6txXIgPd~QzDN!BQh)6G1Y?BZz?lZgWHGi^X@A$~WJh2L=wV9Z^_%wbr z!$>22;K?iJNg|#HOCI;HcjKG(rVKI?*9l)~G&~?X610EVz^Xla#~$+qSXp85Idsnd zLkhS}CL%%J_2HYyp?gIQ09Fa!20W&_s3wWOs^YPKp7YVc%^OxifKUZCLCoGE9b7zY zgstcw^(%krj4{mXs_F`i+UjHEpSzt}NHi#cXm52uZZa@J zndqH%DU&5V+fEeYhXJWqgVF2fMEO|(uQsYOe+dPHCPsF>q|EA_)jh3CoYFdE%Dy(= zlu^=uDa4D>!uwD63Li<3lAGP1(jOLfT0aG=vO*j`c|@OK+C|sbf(535KXQxa9HB>* zjc}%%VRO;JlG#<07~Z?FqGW@blo&HjFM>4n1>y<=KtDgq*-h!h50K{5)&9CIaFB^4 zW9-$j7r38b=FMtO>B$q+&l!SamAEt$tc}WlZn`#e+$`}Beq-*GxpO&GER9wNsu5HH zzD{h-lc3T&j7<#KN2Ekmw!qa46F<6&l_BI~xn7NpXuYRClfi;=+*+eUQuHj{$o%OG2;F>h#L8a682XQmbo(C+)pdjA zOS;4z>sjl}?>OeA-Z_~)0Yw>(yFK#0(wnKdopSSOpR7ro;YZZyv-9_1eptBm2*mXB zb1Yx0y}6yTwK^ADE4$>sx-mqiDqSn^!eu{RXu_PgNSC^I_O@QOnC&l)Sl`iQ`R>Q?-Fe;lvMzqLALV%yTBj+ zR;H38RO4F&f3pH&el$dXm3H|qV+$#wemquG~{T8cs2NHeQZA=`psNqCa%&VXCF zG>CX1l4_#(;(ViDGgmPaM+g~+TmiibE&a5m@`a>VxnD##R^8D9tG#$m_m(m=q(Bxt zk&QEEr{Q=!Pi#vw@=aOB2xdH0eV%SW6EnN+7@k@L?WB7s&Gq6Za~lJq4h83ZWm2vOk++lA+rX#;yd)o&dUl0AF`Jh-fe;Vwfi4n)u z|5g)as`+AgRoJ9PjdQ}eDp=;t7&HM#Mfm#0p2V-iXn@RG;DB|#eksG4@uFg8(KmTq zeCLI7Dkjqy0g-OXTsH+JjOrN5B^{y5vL0qI?Wc%WHBF>{MJ*{8Y>`1jaS{P&DfRNz z8&{ShyyWD2(Vj`56C|btBEGmG%ckmYJ!$R=>(Dr&aXhamGV9ehhSh9j_cHjfHs#_k zyY79e_dE?f#s$HC5|Jv6ty9=a>5M6f<5fx}#pj%Im5bt~p`(Hm$uG3mJJCy8%yV7V ziwCLpqqNw62u!0!j3Y@|!vd8VZ8PpP_7*XkPhnQlQ6Z4i11)&r9h?H&Hyh+(a_fhu z4!@vp+=9zSM2`!XH$~?XK9enhc9xUR4A(Aar%y@7-Tqy8;xLtTuEiP2}E3o7~ z2sA%)i5@_lP0a_lHq^nh5p5jQS*|o7AXHB2QcO1|Iy0Qr&(qKa-N}&nd7b>pb~|x~ zPhO`F?bjHbf7xs(h5|&S()iEyH@{3(DjA=M1_1 zXa*r9Jf0YmB5kDyDVDITr9gak7Pp7S(ZNrHwW6}_^n$UI$g?jztSfETrL2l@AD@&y z;JM`^+@>@rT`WlCpTMoNTmrwtfIk=U zQQ_4!*o<%cw9RJ^GWQXXc)~y{!XW3mG^Q|&;Hst`YB`h$_kYEsqN@t)EuHw|YPf#s zJ&}N4-RqZZdL~*AvUxJW!U(Qd_+RsXlyc3cB#eusq-ez(k&n*etM=@KpkLL>%*w0s zGc*J=M(<4%nZYm;tcY!dpCty@$w}5H`Oo^fA~c3QcTE$)^|gzt+3%uI@guArE<2qW z=0e5qQpn&?a`6#RM-)+-_-w~AvS5*yU%MM~ocDb~N!t6OvARhq|K5^O(#;@$zRasQ zahHL@jFPNN9^>pwPj37$Z7Ac5r9>H9S3}-U88AM1EWH_5PseEebs#ziy7B>151IbW z6qCM>=_9&2&g1CSt@Bx`D~?IrGTdX=R6myFKEAEHzj33&dF4%_7Lz4gesi^+^hPkN zog#(ST8F@|hC6-?=BNKe`fU4ueNrK@26fKn#_T5^K-{$YD+@)2oqV8Gs73j)g*n!@yCsF)xr1YWEdGR*4jG{* z)1L2UZ~T%;L=c!>U7*var;Zz=v=UCTpP{JekkZc^ z_=!tA*$-K|aNCZp-kvyYdfHCb#^i_LQPCeFlOau zm8Me5(-p$}32zL4lzuXb5O?!fDfaNTFbnM~c0nTbKH7e5wa09A0)<&_snM(jh6+K) zo8$6rsQ&$xhQO9Uw8Hx;fWovd#^dF1PcO~QRdCr9hm_=TuDBpWVF zL=Kz=)qv^aQw_6j{L!8r%<4j3X;Rs@cG*EKhJvYtNL#UgD*4AU4Ei-w-7Qipw?jBQ z!h0S!@6}}M6DydEtY2T|3{1mAYWX(b_BJSlQG8%7Ekc@jIX3Z&*CG7;YAS7ebRM5R z&Rwar4ks_Gc04RTvg0`CS0$Xo{evEoM)oZtyM(A($ZP&FZ1@K*o3>=$?XFFplL|$v zW?QL!)I;KbTCEG@=~%+IHItF48sPo!dP&6xKQ=kfBw}C#SRPrtIQx@{=5Bb0?BOyQ zsmB{Ud#uxFq*VI3>(B#87mykV_sFd%q5(E^xw^yNSycRI#S5#ROr?Sclj0< zXdVfEI<6D4bi3r_d~9LPg}pv5%RqhDa5&*DtJLCTm>=Puo|Sr0j0N>)2-<~v%m3i- z40D(qtUBK41KD|qy3=m=!;|in?Cbi7p_>C&vikP1pU-kUx(T++q-yQl$guQ~CSjI> zF7)MpvHpE8Wli7$kt$9!)v}i24sqW}O}jo3OFNmy+X)=4AI$W24Aj{XwF3*TDzDJ_ ztNAcgl8Xn-*R4VcHLFtFxWnsQd4Nb#(kdG5L^cTyz+EPcx`WK&#j|TLAdmSh>xgL4 zJ48oRX3or9;JaDZjk@+zH~r%77#2I+2^)cbtR%A3jv8<2Aipybt0{uEd83>5bZ}Gb zXH#zaWVO+>RApfz6)*LpgU6}P`>v~rvD*|wRk6=ap2sGu;lBFV6%|UNKeR%WZAG0W zrC86k#46*@p|?jn_L^V@{acPqlF$l#Zh<5lowgN<+I>~Aj?FQ}U1bx%ylH$*N=5U3 z;}dV0kvV-R%+5DjHe-8(20t?kFVo2RQKXC&a%78cCxLA2_ts8UlNYiKCQp%_Ps^dc zqD3G;Mf_y=1v5R$D7$HDeo2JX)pTefE-o3n`*&wKuT1IKwq#2L1nz9;$%Yl+D(c0j zn&5g%7?EW-U!L5%b|@4+lHP0qZ|}u_RY|!QO(_mjlB;O+ zzR~R-mv4q6BMB_ly*FF+_|aI-j9oa>sJLBN3FX?PQ$qz8xuQYY!9jEgS@qAcotY$D z2gIGR+V3CZYeV=ZXk7J@A!;ykiOiyw76pQ13_Yo9DUF|kxHzE*`)f;ZQI{=$vaZ)g z)aceAxUb9#`p$zjNzHSIOMWyw!%BfZXR~odIv^nscw~Act#~A->&DzI6~}xT{W6$5 z1wHL?HaRl=^sUh9F(z`U9hEORMi{q~+-DT@zGj-uE&oy3tXECj*&089NhzQfW@Lt} za3KIR!JwQb@tZ_RAKr^uRE^JnWOM;<<%u~iMx}*5dwhaooK=7Ae)H-vAp{e@%B+jU zbN^ktG6MGuhHmCD`n+19NE{DV_M)eT$R`bTsB95)7K^~rcn^AAtr=-`Bq5E)BSSE@ zk9l$S!ptv<6nCvrRxc~40-IvmSA`M2p;*5%)d9>MOL092T<;RaokeDUzqQufV+QhU ztpAuqVSsL*-+uLeT6{WzrrINo_WFi^=ZDK;kPq<$zTo` z!W?}SsGBq%ro$QlL?PWt>Gl%#Gcd&&!zt=Sp?4-VMLvS3N9vCL1>vx9cvT`P3B{kB zH!h+e^39V(@=Ip=&ghzd#5Xkd9lpwP0lcMMmh<{@>!(8LD8!x}V-dKR;cIKtg`8DW zoDH-2EGE}-a5)$14FRob+*s0V$ILq!#BS)FUw3XWW=c$`?l4A z*(P>pT`v|j5jeAToAZQk+KIs*)KBs4Fr50=Dn3tOwwr}E`gNg~m6?CyDP2A80fdp| zNVKPDBq({5J}#gqlJjEE7Kf0kePgsoz_L_WBxhcux2~~E{FQtxD^oHy zoHmY`)jqX>rdv^YWvmt8++K6y85S2%LHRnnWJ+fz!%^cswH{XqmhvmHuL`mgLh{5s05dFyA8XOz z)$v%S&?ch~zN))cxo_l|*Hy4JY}>^Ncr9bh{bRJ;P9x62xm-`e94kXb{V4m*c!?L;N9|HTUMbWP;_>k4HRBG6s7??Agu_1i?Ea5VUWn zBlr%bYq2o{KTVCiO$2`hFV%M>^TPDLG;-=yRyLjf@$8XnJZ+>fpmfa15}G4QyArxS zImkm%3eRAQYhr-u8Bxb#7S64g3#PRAeb<59-Z!Jycg@fFuKbujOt|ax&i>we614So zQ=U+i!UvQ$R(+&B1h(w{?>F|g68bSi+^5a(eFwKcFv{XCM3@{d=)#qp+(Av2M;%(V zGm@f|kMHeJ!!IgQH7&#^ay*R`r8nopaOJ}7!$AKj$&lkoj%_9kvSw5|tBgT0mZZRV znfk9rz)8`jv=y8z;{X}@o+drEUBO%wVFlq@Hc7G^Mq}_Dm6$N0B&WbiH_2sB3iDni zUr(9_`n*-V=?W2XoiIIFu8`XxXg%IUUR+<`y|Xxlq>@s)_V=^1Jjfx9T-7Z$0}175 ziN7Gp*gINfm#7G9Vp&Pwh$Z;NRa`8v2?&zsEplch=HynWKYUVw zEzL!r9JnILtZ5qjMyONhr@8lN^zM3N5sYUg0T+?${(WyD<-p6qL&XVIUVi>hzsHyO zWZrwil zu9vVcpLlmRc|uvF3uz$tBOm2ie(*oJYWPK><}mPWG=2@C%eWZ)2!ZZVbO6y%e7$qY zIR^R3&4B1VcP+{ek1UC8EHL;%FqD52iF_U5U{!r5fl3Z6Kh!oGdmT7nh*%&J`qp_w z=SMMwW+5ok0YNMf(DW%Vkh;H$zF6VL!P30x4-Z7#CE8#=6&(7CSYq<0y^J6pcqFP& z1kjJ?Hb{@Z%k**$W*e9X2iMLFKL=FR%gc6QC>f0AaI6(i2q9a@YWm=_`f0aLS*|NjuH;`@tU~M&k(A>B z$q5R2j%Fw~Lq^n+hQcY1CbtVpRoq)Ood>;Dx0=i!x-ix_i(DpwD!}EpKMqX2&-A>+ z$Khq225>(hL3J;@*?$(V-8+YeMX3l*%NixZ5xM6pGG*u9@vd=#PY#3;9S*mZLzQ3N z%18g`8n`gko*UY*<<_0d;7Ar|U`zT@5q-K%RmNRfeY-PYq!HKN(|es zt0WRBJopnRY;S>E>L7y~EQW^~s^l5{ZVQ18-cKpS4pHiZwZ2fWHu)?HrzW-+Cni5i zst)DF3XXJ7o0cD1stKcma1nk~k$kpry4L~$k9IUXer3ZVz)`y2M0K&kdZ9TrAbO6V zm`Bd*>6g+RP-54^^=l!jG%y7Pi?yq{@UxGtIyg>=l^8a{zn&Ft$1*g zMf>yb-V>bZ;flHN?~5i=dlaYc7W9ho)>6dl29V||kE|i#YL%<^xD2~J4bC9sDl#i9 zm=4PJ-Il>4&WtBm(QC)Dq2ndx5{5|g zj$J);-#ixm06eAo?YA~u@NoXwx3))HJrmaka1h9o?5H1(ISn59c-_OI0=MM}!4Gm~ z=Ln6YI8=#M&d-RFz{NhxG&H`op7eaU(Uv@9%Nmk0U+KziQ3x+2jE8M)y18K1Vr29D znq{Fix+&AkFqNVtKqzHT$JlP$oc8gcD{J#rzh0oh=V9#*R$a}nZOJ61x94;vw8L!y zh#Lsn;!&fUzOzkCq4`=ktk{ki1=`DuE>+;XSD?+wvluTw)qFxQ*TTyFpRwtRFTBbBUWk) z%3=0z*zwIMS>WF)_86x@1oQtLbe1Jy|8IqbnV6Wl-SVf*BYig}OT6rs&P;UxdUdDM z0yASY!-}@g*Tqt8Jk8Qn|1xEo8adbV1i2nv_UzGJWVoJfl;oS19_%4!_$b9l!JW*S zl(0F{KuzKfAR^nMOwwA$iEeWBSU#M99;qyIS$bq<*!j}-xgu`%GR;Zn*>W>|IfWjE zb`!u+%|fIkk@Ew)O|HR(4olV5a}R0<^6IZEggLryVlw{zqAGWhVf8JIRnCc8juty{ zUr9_lb*NcG``@1CJ4##?lsqJ5`f#EaI@VOoDm$0~K0(Te^$(s>c`-Z|!vPv;({wFX z5S|sH4&XkoA8mrX8FE?MR>zrJ(9%*bcc~X2Zo{Cp;!JdD9xa2ZP^D8j-YW;^@WaAZ zf!rb6m@8ypAUMUzJk2Hq4H>_*T8d_p7sJU20ttrk+X(QD(mr zz}%cco9-5+4;XH8-5bGi`pRmSjO6`QVv(;G0CZ%fV*dm}$gWGX_hu8pTl`{NT{ zEx;4{#q(8bYGF=TiJ3tR)yr8m=!*VR>@+W@Z7e4EOHC5pHJMWn_cu99t~};nWk$SN zpl%vzY7E5(r~vl9s~P)}aNnQlXaQm6n;c|$#yJax-RDmIG+)~BRS}6zlM;wesi8gS z)hjFgt}bWpwBgVU7+P)d9ox4Zl0K>40rHlZO4avxXFq~OOPB@)q~r70@0W~(K7~_s zh84}fPtQP~$TU$@E-xvl(WTQ=%;=#;X1s11Pua}LaQiH4eEiMF_;U+6cb@HXMqrvj z=Lbg^bdIprfuCZ`X(e<|D_JeKw=^=4BZ4PEmi?CedRD2}3l0HWW z0q;>(_6iwKaq!mM%d%!F=qX*daLt<2p)=Y&wrIyI{mUo;t0+U)YwklTGn=eWp*>{m zGPo>8X`|FWTIv^3Iv@&LL_wfST_FM^KJ|2d=YxxI_Y@N>DZS%k>w9m4N@}g<-$LXg zq2U-feDx;k7q8Pb?6G3lDv{Mr+}JQDwbj6GW>(<8SW;PWulQG&z1lyS0eqLqrZE2?8W&&|(E8BIN-Kx9I z6^H8KxePCW6>9fh=N5sS#c{nb-POhqqEadcEHfu{WkDX9)#bO^+dnx!u&8}*){Xp! zmaf!G69;|qD>*&QSa4R`lA`z%w=(#V0XysPVw@k zd09z&bpL4pcRzSlJa=QsBT)cRn?(J|8N>MH3R*v7E*^LZ;8*~see)PNNii|=)0w2* zS8Ac>L5F?3UiJ8Yn*~Fb9Rdf=5hD|FaAFQ5lSX3{zvK?7Yw!f!LiUHto3}6Tw+4^4 zFDr&`pPJ_;By&p6ZqAQsB8g5t`SLEj_xoWuk5+Gul81~hEd?*VdApY+)ro%p{BhiN zobR^h88xT~`=$VtKDLJ)y@Vg)Qvy#$2tiToRTIe+m4eidNhFgTRId$tiO7f|r|}zjA$^exEG`9z{<23rQg%;G}>u zy~a}4k=v`6_+Hef_p9;cxrKJV{G@h6?`E5pq)$YRRdCHN$mA7W#F3iv(b{g2^!gDD zhIg-N3jt?7vRhP#6eOwJ&prwMktH!S>_%@ zpqv;vpm^f`R0jy*+b=kq84Djex_@t?ksLEmyLE{<5(F0pdpqmDwHR~0=W`;H?qB_S zas6t8>(ik<9mrw-x^oyZk|z+~Iet-*y~ko=50JHOh&>$bZSY#Ln%Vg7c(!&i_F1ot z6I0CFc&hZniJTJ;RA%Fj=+?#Kum}~N_1B&Fvr@hI-NBj#<@>vH$y08<2~RHA)w^(! zpI=(tz#Ls3FQm`vaL4i2NB1O@_fVE*6X8{hjkEp`^Bgt7OYlwb<*`6aNs+~I>s9+C zu&yuR3EVwwoSAssIi=G;z}<40{mAq2hIm)4*q2(qG4;-&kgrnc&9Th}F)Qss!LEvAwT7FL-o#@xDZo z&0XlbhxRNl*qyoqiUN zC*h-1Q<%Rbs7jfKp-o~@KKI=>Iy%%9&g^n(5|2)EFjTKRE(JRIh)7UK`#U&(mkcq} zo2Kjaq9F&2shI1He00^Di&FE zV05&tedbM;LpBOlcG-Ic0`5S?^wt50>gT z;~xct{6*`Fban5BK4FpN1ORk%)pX&%s~nZB$(O(=uHhYh3466ggV%x@1rk;@t3ZEa z!}z+twt?d)>=3jQ>T;H~(?H!>*}>F11e>_OvF;iwU7aj_yK#B@nM)N3eFqncj3uV!57H5cvS=h*J0t;p^Joy`}C8h1;O=*#M5`SP6Z>u_A{ME&MKO?vXl zUvJg*LweuL_VT`vcA!yRyk6&kl=yG6Q^jW39gH=c|5(AIsLmkUU9SGu?cXJg1K(pm zNS6g}ugW{$W^q{S`wpEU+Sf_$P4mHzc%MW_NRTVThwsf7}Bo(igz(pb+H>R z_LiRW-&I0u=Y`rnqT7}1(f!{PHEykwM&%*sU!0=8xWsM(@cz!EHzJp+9o`|dCw^a# ztCvszQ-!MRCkd>WrJI8r2I(`tDq8H^}9h)w-i%B;mMT=uCLG>6-zL*_&XB^-ce}9bg7z3B5@s?7y``pFvC4n~#f| zwitKx%Id-XUFF^7=asmrLEN~FyIV+(2C3n;%(0WF2Uu6@s(N#R`FNSO8#*jBdq1#{ zHQ}qiv3!n5-fwCGUx5G+y^&FFX6K!GbWP37WSXLRmVK#5R>M0NR8c8udY}6Yo`!T; zj!@Apk!%NQjV@^x#@dvmaI-=@-LQ2_oC=h$((fveeJ!afP@kj`h2ZK!?KQury&aT0 zA*s0tBW4YK?eJ2nN!Ehw#0HI5!dSTxih)TqN-6mgGK_EiaeBii0HC z2=iMX3dIL=cglC7!>RJfFih3AhK7Q-*iDdI1U^;xBt2Yt2PJ4+s(VBB1&^pY3VH4a z57zAl23baVRY?IuZBkB{wlu;%A7uH$fu0%;1eKQ&hgPyJ-DhZrW=Ci5M-xgR4{#0t zA0JtDjJ}?ze$4wntkSKr`76YX?gTYw*0`n`ry|NTCAKKHGC#()!WqDw|0En>9sWPw z4gW6y>j?srdxiBn%M8_hpE%2jyhbmJq=7iB^4aWWp?>9UdOq}c(m*4kit`y{C0kXH zAMinc-SLL6NO_M~qZnKhe9c!j0bR-;<$t)JoR4IX6cAT)YU+fIW|TeHYs`b&-2O!e zYVieU7bo0^{ko3D;zm6TC^<-7R4^v{>{>!{nJTF!7afb$Xgc3|s*x@6%Kioq!6`8Y V<#{7gV)6+K3u4gI%WEiL{0}o(==J~r delta 23362 zcmV({K+?bZ-=}MqNJ!`RIA>qiNj^aE8+BS7aM|^-ty|*)H{?8EdN_T(_<+0awZlV5}Rf`SY$KWP_$jg79` zv^}EH$22P|*7m*W%!mh0U@IJRe6GkvEtTcq*%Af3ddMF>ltq+h;K4RWe-tIijgk_E zsugFP;H7jIzIR^To=@@`oZEncS*p~LC43xp2PT(I5nCBG4rT0v4%P>M+v$PnuOfqp z0b}48cj(#ABEWB<(qO0{ZDh6{ojRSF_UyplU-97VfuRTq`*g7zDFBkr1d!)fwJTaNmsjuSgQ zI6+&mV>?QkrXfclQJIKI@rga)na=As?9P|0Knuk`y6JS(ge`&zKBe=47kq*9-~i7Q z((aaj-1KKxpis$wiUzX+Q>Hy`)d04wjNO{JqRZkp3Gd$DDlS?2X-jk*oA}4Y4n0*^ z)z-0}nQAA9LmRxTcmkK{t3KJmuGxi~L5>)(xrX%S=o?=61Uhik!7$GS_Z9?!UPPlO zaJ*_XH-MMYN8bL4;w}S+Sh7jhoU#|N3O%OH_UB3>pv+D=s9`~zWYDgI96BNP_ReVKy?EumCV4a7~>_bOKb zf$W44OA207)bS#9L^!@J*>neGBJl3@SLB3?4m=}E&+{W)w-R0yQs9cBgz}0t5J^T; zI?Bo~CvX0Lzic>yG?c(IJ9W-bJwzFHk~2^eQ|=}JgYza;c~uhnsWw$LHWL8TMY%nD zk;lP-r7-q@Z0Mh@7s=DMs;CzMJ5-&jqn!Ko78?syV)>!8zwnZK9GHx^b)!%f)IT(E z?;!wZrK51GEZqajz`$mMD9`kH?_eAY^LD%rBIfFU}hVX=SbKi`&vk+RfRHZp@bNbM3s8qhEC2g!A=lqqg*E~ z2lXMael{8#5HGunIR1eVRBe@toq=$*;ie-2B0JgHwcvoGup>m!EF1Jb!>lSm%RqN4 z?06l2%vW|64%G)+js~0y>4yC$lJz!WA$0XJ;`8mjAZ}2tR?>;vj)~y(pjX40~ z{c94e^2S1m(Mb6O4t;B3actw$A`l)}X=K(ssU<+2fb%ny724l)L@IS3R6ax^qjZDi z>I~j9J?}oU<6iHqUdso}Scj~ZFu@ral{JEYu#F(Buy#S zVf}T@z^uNh%GZWc-5%SmB30w06MRzv7w|aPZ?(?`6Aw4LBl!XE%VOB+-D%&Upwq z**(~NGVs~TK;!=jW8?|QU#dYM(8G&=U=6OfJ_8rj-@q7{yW$kZ*piS%+HLS5vz5bs z8Aqsq!dI@n9m&I&ySt-`!8PFWJmx+VDSu>DX3}%|_E8KyVX;9m@^zX;dF1fZ6equY zd=vdnorJpkc!hHyG{-!ysh%>xSDQk3gqvzk)}m&V(62X}^?|4W)gWT5AOQS*S0ZBRP9gc@TY585GcQ;GeJb|OVEo$dR9pM}mY4I10tI1y?d$e9Q$?O;O$jMlr+FoRL-ISG|31C> z8z~w(gtP1oLK6cvI5?9rAzlK~lCzN^$N_&}!day=v~m*dk9+#e%Ki0CG#iQ>+yvu% zw56nrTjs2Jfi0aHe$S-a@pc=FVHl9-7{_5Rq*$H^a^Rrd1f5U%>K!pE^?OQ)Vp7XS zvxTu1s(d1=e1Z;a&>*-4m~apYbrpz$Ps2Bd4y~HZMfxcUhSLyk)Wg*2V*f@*lox-J zk{R2!C$q@eQfehC8>_~~)-sOCBr-A@@GH2eKkk;zw+pxUccADcLuxj>{@LkpT-D6nY z`sz?I&im`)vmuxk7d)44hCJ?Zu3>+++;;iv?Pm2UXMXqP$o=1oFLvvF@dZw>DQ4RK zPG{*tsb~Y1#Qvilx}N&OmW+GFjEx}vyoi@npdf|v)}0zCog#s86;YN~PE}g(31F$(EEB3{4-h9?pMPlKD?` zgu7@ld~`8QfM9Yd?StP&$;ATTOJ11>FFNbO^`H~9hP9T{HdHbD1aGMbUG1r~w>U() zaxFiVNkF1~v2Ihc|9H`T-q8-}s8XkE7qGX-O)!Y0E$93%4GH*0k+yW?s1XG~+9d03 zdxf=}pxw}bYBsk3OTjSp;-`N?oQNlIoh0ZSK5_UV(u!cf;G9_NE7x;0AB~gE-jsPu zbM#Sf5xH|$fTDx~ZaR4@azW2D4-r!j!^Lw#F~hX0sbH)4=p;do!Q_Ig)(coaeM(XV zSytMuXrU8;#a$LurIG?+h~p%NOo@3~hmCV;2rnJZ>WBIO5=U5r*c*SAZCmzBrhE)3 zV>{c{gBzFxHyCVM^$5p!N+Nwu<3CR8QI}IEx?D6L;c~R*Sipwb2$(VD8dtO}0G05r zotrV|>xOMLNgfC4(lK}ItKGxqaJ#+>JmF9_zBYfX7pu+wu-kli7$jWakr5d+Ot=D@ z^23Cy5N0)#Tz08IZ_9rjG)PHWEUMVXIR#8=gOsDrQ?wf=HG&e@Ap;y;$U!GQ!&aF| z@a&wL3eWF~5yp4RpPK!U8X9Js*4+E5+@a9NIDwDzt=@>T!wqWiqFFwM+K6hF_k^tL!mvTk}~ zOOH0kz5qYuD1OKWKZ4azEZ3(sjj2iyNr%*(O?Cv)$=e#GYkl0^2Zr;xI=(C0C`fi5 zc8{B~#ggPX$|DQINE_lgT$IdYu(9MZ@2GGNs7iNxP1k>vuM;)pYr@?evsQnjidRpn zc;STJ;kzg^BWa169LzWh5~8Kx%-)Y4GQx>{%~S?;SPR$od($#h%T8|+yZeT&hIrgD zxTAG6$|Opm;r0mDoI=j(!dSIP4CZ-nkvd3B0_oOFI80Vy9Z=PONE?OY!c@`e2M+4+ zTkhguy`X=+?jm8e3)1VroAF|~cdYl-W8$b%{wv2stXn)?(`|*}8Szh$weha$fVfU< z?KN4taj>z11MLCXWn4pXtNL+kz9eqrlfGo!T8b5>;})jJtsMzl*$Z35G+&>u^9EIg z(fI$)X-x^9;SgfCzu9~$AoXw`on|CK6%6AB20?!n%W5S$pc)gBN|?gZ1tW_1)I_%CSqb6giYK=H=*S&ErNQr)fdFZxV(Er z+zRR)y_e@OFxHr_4gZ#$&I@X`xbcj(M=>Hal*@l2`ag6Op_Fji0cPL~AH{HE^ZQS{ zhA62?flhidTR_|!3xGGG91fByQf`~#&OBD`fo7%;iR$Bl?zJh6xuvtRRDKG-mmJNb9dVW-TDOc8V7!F&N_1wSTi1Dr22{A**ls{wByP7kVWTss zutk62B1`MlKacBibOStp& zDRR=w)g&)JR$R)ry5gwU zPWiEbrM|h3fQtsNGY@L5=5!>aSK7a@p5=VQP9Wp6_y!whn3jDO| z>+zfy^hMgxDJG%5&Qp?f4+m6pANf3yc#RZ+A36KQP9J|d zO8G@Oa6{S)l%}HFUcYi>$PDRaoYF4nUH3{)ZY=1C#&8i*@`rr3AyDNJwJKnXmr?ja7TU;%N z(jspHTI(t+}j*8oM%YygnlHt(gk|xw=mvWGqV#opzfSzgzAaa>7Nswe}Izd_4VX9dQ z%ItE?AYzv_REMc%#h^8~oI{#{H0N3{X_t!#8llYA$mI&Ae`WMjZII+jz;gFVl^yA3vp1RvH_DOR--T)V+Ttzx@;7NUaU4*jo7=CFoqay1~({L@5RPw zMy%QBfJ%a)1iSGBO(L70uzebCk!rnG>)qxhMrTS+e^?D@AM)(Qi`o4B-xuc>tDoJA z7Y@%|(Y*Tb?_ZaW-Td42<$Ck}*Tvb=&A(gz>*8$r?#(MauK%&PSh&^Mc_I z{It0IxOB~DQ}OcUZ2oq$K6}61qR;u;?_Rsp8D4(<@_&fWhx9e2F75aAwaDM?X~6UP z|NmP5{|)kgUjGNmUq;DaZFV8Wp7?jKuIpT%fAH?ZW3=O*He-~dkMcm<@IY6ORkOXW zMq38<8Z2_&smLhJ*lw<#72_*c4Dzf2xx}+#JS)a6#b}x#C4E+msTJdpFyn!){+Er9 zcw*%`q-+N$9JI^L->g5Zf!v*l!Fq>DU2Lsgqma$S=5H52FK1Au75d+PbRtGLo4;Fr ze_UT}&Xyl-*ROf==jHju;???7JH>7f;hJhJ zhN~a9&G$cDtbUn+L!0yErsLiG5A#3HzdfPmiy=9~+$e~b&}`gz(`F*_q$p$dh#W!p zjoX{_@3xB{FPENwZxil7^jPu3U;pMZf4R)$R#%sor=KIiJUI-JNCFJWLu6E|Pqo}4 zyS476C4alvZkL;tS0 z+g`NoMcZDqPN{=fu)fslQcS+uezs(vF^QzrYkEKhsjq^EVTe#P+pE@6eFU#!70<;s z)UiJ5{UDbAyjqwg^I=LPnP9r$R+`+)HpP9h4gXSn)v?R*Fa)gRBg~UR9%N@7)pIU( z=HP>t2OF$UST%{vJg^ud%|k1@EYtn4QgwdkRt)Mf3rX)1+pm# z96?A_;iqgXNCvgFPPlF|L=9iXnt=fEe^p=MRT&13A|De;f&pW!ToPa4&x{v#PRI zCa$|-V9`My_esog!4M&gGIjsW_}t@dpPx5;18;*XWW)!W!QH|`Disv}4TTL>K4n`pGxbb6Uf1hOB_%W%^@Q0&l zr68gkhc^cHl7=uK7O*4@QV<$NidZw1IEoaGouTF^@OjPP^{|?WAh_jELP43W<+u;wZ>E2a^esX*2JP7${Lj zJ_WB`kxkguH8J-!7_riniH#*z!4fQsi4Y4R|4nAi5EIcbQ*=rsNE|rZJwl)CJ{RJ~ z*ra~yHS;l{2z*4gtBf5K0h2W6(NHGTlw5s@a*f2YKU-#Ff03jqqvry7SPG2fCKsqk zZMLio8u9RU(S*4t2~@>yUBWkJyS|5f%P#O)cAm-+#Y1dP8mVSMpaN$?NQRN?7%rCh zxQya3B(NY12f`UC4cVPe2ubXIF8Pw^h%DhmzGl5f!nUg@G8)N*YQm)=%mS1Yz1{NU z3Tat3j@7O7e~5j4yF@z_Lf@u=PGZoo7d9OszUaph{aE8xes}1;JxRP1J3RGl$1VW3 z&b|u*s0#wP3j(+cD!8NLE~wxxsNgQB;4Y}(JY~?p6uey^PnHS<3r1lD5tqQ+_)_x& zN_0@X&>YzHcBf0T$CnVYu$}qA++@4cEoO@Ew=-uKf7_)$yIgF3;Sl(z3Bh1N7es11 zFhVdr6psaFf8SoPmf;DIR6@15ApR$d8LF3Dj=9Z6$_Sv6iyVAKW@S6ZVn2frOUn3| zF`ZhtSbOHgYU_V!Z z+F~Kwf7(Il4@+uCJbWB(Fk&YAg@vWvS?y;uupZGsA%F>k0soXy9HKUg*Fa)LWD9>V z63CP>ob8DNfds5v`{_rnjVY0&V`XAqS7g{T-=g7Wmm=>Supc?|Z+=*k?%B4zN%5cR znj_zF%`uE}&D0*E_c^!2^$>UaG;rQXvF+zpe-T2BUo!aVZ=&=!PWobq_FJBK7`y9T zcIQX2yB7*^7Xon?0(DnDsJs6{-Svj4_CpFI95!3f3XNf+3wUXklHxH>dAZubj9w!V zW=~^r#3@e_)pJPKIF!T(>U+eqU zti{Z;``Pc_dqG96qRu2{27L>XhC-Z}*qB-Qff7n;Y&<|#Rt{!XR(1dtl?K?!2K2WW zK&1t8a0Ejk{C^loIDkx?o@A0HPJd5!N>B(;-q{Ao#sOsG=4a#KXJrMlv$FF2#}Mkk z50o@<0h}Z_cK};hIz6TN&rcw&DIJiFkB^7(w>wbG7UTdnHGu$? zOq?u1woe&NO>BVbP*X6-$^E}V&4oB zk=e=3i4kZ5G5cj`V&e#X@_#pRF#+3{ynQnGUAYNRT1*va@XQUJfJ~p-?#}Yp%UVOAt`N_^fjJmrX8wzWnX?^> zCIoEn43d-l+vG_E_*Z5DasqO(va<5=O#_Sk@|zdrx`GI$y^stg881hv6^A|LDZUp>9A=CJqiD z6FV0xkdKoK$ivA8^nd>EDF5*P5BJ|v6%+8^wXyyiP!3`a1@iqx;i-53q2coPRM7st zF?7KH&Y%o^8f6fW_TS0ZXXRowefq=p|2+S{UH<>Z{#TU$+vNXmLDJ4PHow)hzxDqg zwTUg*#{F-Dr%85pdfEyl=+kaM{?}9+^w%aUfy}_pw*RY@bAK{<+72;@h0Q;01UpKD z-9Tn4U?)?{zeeakukaM!1`GkIKpnxqZVix$jg|F(bWi(cYW;Li9G~9mw+i&MPXC=! z3StU1`?YuMT--nt2L}^(0PE9uv2$?&J=vZX)C}bId$NHn%n+#4lMC>Po;T1O>Hzrl zqTF0S7O`KV-+zdY2gqXbFX98TnEr!)RmWoX2V?`Xfc}8&Ko;{qAP11e;t$9PWU>4Q zKFuYIrMsQw)4BQA;z+SfGEvr&8?xfKRy{{(w)3 z9RGmaKo+Mz;8O>j|Ap*Np*()W|LVG_vxCFaHTu05Pk&GHAN>3Dfk18`Q^4XJ)KnnI zswSxQx=M`5mFe4<&@|OI?G!pD&qaqe=Ua5-_jHvR0c#GIV(&k9V=iw=(O!rykv(|s zf2uGuQb=TRWZ)QDK$~D*>I|(TXbqPq|FR(`@sXfJalH~E{ zJ&Qj?Vm2`kv*$g?XBr%Y0~X}TX$`;m^>|v1DSvb!OCDUVM^L+gLfVG6ul&4`P)RU*V4@@jrC@9c2i#{=xRT)$dl>|cj)1*7uUhkLEl?bnHG63vy_*EiTFh&m8TtJ zC4aqZJ6YH`(g@{8YYQRvI{|dt0Kx01IcW{0&AV}}!eS*^xQ$m|=O^8DW!e%=Pou~+ zs@`lQUwKZab(p543P^xwIpciXg_ zYF*+w^hPbl%kpj`F&Jc=KDSCKFo{01l*Pf1sk)=@>3>7v zDvsX;;B76cD_b-AkCxn|)MHW;&pW;IC@o?|%|9ta7SPi-O+F9+KQZcbXRt%7{9&kC z4ZFrDywO(&qd?#&FIiQ9*jz>26`2UCy@+dJNtAJmbwK^&{`{N|7k_);9@|H>2|3afGu2?!K>p<#!Ekh2ozOE*K+frc zEBIudRQ7q%MH)Mva#yIyg_N23acE>jo+Cl#tmG;R4(2$Mr@P^M)n?qpEn8v?GC$W? z8FVS%bu|9AolJn&Oa4y7wdl!_f&$(G1`^VZGQ$Ka#iE! zeMC3hSd)=sj;`e6iLwOPuj%=ehNdXss&ie$aiT6h_H78$nHMEN47|-IW5UIkYF~nX zwn4Hv66bJatppOIoSccpl!Mg0UwjSrEy*J=EP!pn_T)`Ap&BnYQm8UoFKnF`ty$|Obbnru3pI)vSaJ67I3`D4Xd^`5%;L}XZ3`GVh&3iA!D9a) zpmHzS$dP_3>g6OGjQDbjc*4m5*SbqnG%Jh8oWo9EN%9BnsxcAE}u zyjV`ra1Lu)YFFq{RekSGbM{6~vk)f6&IiBvxnKO;`>v9Q{;6hY%cb4&IFZtxzf0mB zq;$D|_CPx{)U4GWuYV4aX(pG&rhQu8QQSvRm3(0PRYR zglld|wtPq5EA|kIaD#lcN6|Cj_l-WUKTtb^Rbf~%L?g8fwRR~;W@7kkSz|WvIQgsf zH4@K!&kR%Uho7DBKkHwcr;F27!)s1Tm`m|8*FKliX7ya%9Dk6R`EV4Rxbv(YN4)wf zVl-GBM~%?#m1WBc5*GSi$EJeOTIfZG()s10m!@RIYHHbRB>m_iec$3q*OZEzZF^Rg z1L@w2jY!^cmN!EhgfDp7)PlrWoN>`yuivoyU-@VcX+>58n*$B6EG5BvKkB+VKpA8j z0>yee`E2c^vwvG+a_tI*II3)?Rd)!h0dHHhF5rh}z-k-E_5DA6@`p#h!bkBvzM>&Q z(l2T($wFU3wHwJp!JVloW5JC3ub1@l@(kkL>VQ- z2LJU|C+kx2tFYAbDvKCme4v@|9}0w*xqw@k;P5Q}73y^u!*0;SpvGr-%C z)V40ZYJb=!=Y()a4V`;Sx%_osQrlnyvz2bak7%5@*fzMA<_>t^<4p)n;6i?Bz-Tb{+N2O6vIa)d~X;gsNS&>haZ$*t~o=B4Qo z_5wOl?TUOFWzDj!U@C#3h_j{mUe(F zR_>2otP(>g_s3#5VDl5pM}I)`5$ksegkk9$S4wUIrVMtj>x6&_vCp z1%E-ZS<|lv1@$*?NZQy~1h3oN8opt?d6Y=YZc;i5eWSQV(J5S#vXTqWA`KoHY(Q3W zbMG8reExu~21pNBQuCfQ1BJgj-D1|)aK=r+`)%;VehMJvT^7h|Oq=YOJ_kc!QeXITH8xs$FnK}g?1C8MA`vGLTM?w}qM z|5*b)T9}AlyF~}pM8cOyitr!6B_Fg1o3Sqy9*Q_AWBhlHMW$`XQb^yDh5Fn9UNE|k zuL&8}KaPB*mjLCZ8NJF9cJ@v|;9KUi+$iBwO~KB3y5YMh2-Su!b;fzGk@!X!<;dHsg*zb5~*pOtq1_`ACR|5wLzM# z9cLGY?a68fFq@?mO{0kBr+@Ez9mBdS4_YctXbkA?QqEK*vR*IS3we``;%C!*spBq^ zxn4ZH*sfH1k4xynw6$bdDF-T?Vj@C@_is!V5vkKcduLN7qeqcBDj`7G>6=@ORO*Cy zKwGamRfr9UpWogy~X*y3c9gz&>lWRc`N14vM(qoZ`yMMC6el0>to({#S zuk1vy_hkIJl4X;R>Ryc}IpGr4%u7z$(wS$@J{9@x#ggt@cL8!%78&g*efl&{u(KlfYU<-N~JhG#}7mYiR$Y!K4aC7TP zigkr%ov3c3KN-RRHaHQpnY}swFuKyNrLobjez`n##yCBWbbqm8UkG7RVQZsF9Q!Q3 z0)oV+)H}VMc2sEgY_6Bv@Zc@pH_*`M@hcdV1K+QS;I~&3QaUc_*l_^bPkUqRwvv$n+#BUkr$ZeGsxE#6C$DBqSRU4 zeDR$xXbF7>5m82qH)i(A?FPq7FFQSfvL}q8^XExwGp?^~wq?B z-uw8R+|EcOtIO2ojUVH)Q6H;t@M;cO;eg%`PU=grQEi>&Ls2NeDvo)j{y6A%4%MzFJ;rTq)2&SXovkaznAKtN z7PT+qET(0BJv6ab*uAmNK=&F9uZYFCj@eLxhB0a#D5li?ta`z}58+vhmrpVI z^ZmA1batcY@TI)V1taiYva@!I(rDK0Q-8F&2TT7~O4q(*sduFjxH=o?HZm~%``3A? z$v&Hw?GXlqpc@~o7LrWkMRBn{YIN0ynUmBz6QNq?EP zt{Oc>u`oWjax6O~4uLqcVBMM60!f8aUoBDTRX&x^ z=BwSlNEEpJ*C&So`Qjx#x-Xlf(+mgBJM~bF*TAIWOL3Rb<1L2IgrpbCmT`!0?^8Rh z)%Tgxzk0~QJhBKSwHTT>cr+9kqbHF*aOae9#}Us!B#wL6I`NFU6Z;v6YlW`V>K~A8 z@msUkF>B7^*&^ma73FU~22OitCqi1JLt^9{AGU=Loy)6%u!`^&kYVjP6@Ll56=j#T z%mN!HcUW3^4|o6=Xb^6x}x;Jqihg^t8Y$&K$%=??ST%^&?=vp{V> zxP+f!ScTWsK=>y>yIF66uv==m$vjt@&}AXCoa-#?Y;0&3`}2#FN&P*p(xolQ{s#B7c5qEKn1c-gs?n zyHV^SL}X%@`t@?4ND8$UTrHpsdKufCBTo6|Ff!I>5AhX1=`D_Wkl4{xq%;99^Yuy; zHJhnO$g6aG@`z<|A*AZ_pZ%MQ2KW#2q5MOLGQ5!6I!XJ)&3&YMrF@*y8Px;%x0qvT z9}P9qr|Z#rtjvE*fPYqXC@+fIXxNSvUjYqZ@*2}^62(eum@cQoN3$}ghQ5;oxaHpA z%eKE{M6Bp0h@dNIpxw*XQ`PpDEA9}ptz)S*xnrM^jJGp>1c}fece>hvExqKsx_{IK2F##pNMKQi%M2Fv_c=W8< zKVs6(VSjklB!z)=Wve$&`pytAI~q1cbsr+AW>92KE-QNIC_mmYD$#q{LDr7IS@`8= zf^+VT)Qt0m+$8D3+|mm&g^-_sU6MS z9@Dnj-HSCjQC5Jwy6N_PZ2Z6ft$a_8P=)u_e}CHygi)XekaGAgZ3!ye>8&c4rQVa4 zT!=tjPcp7nCfkHzj(L>mOo3az)D3wi{N7OE+4*{}dX_>ewjeSHxqRs=u;jz0@;phm zQm?Rfq^hk8W^2*3_U)^{fIJzKN&@$M96b zsDCG&14)j@j}c8gkeM{oLFRDVPWbCntLg7i5zU=z?t7S9cT#Qicn`8f(#g0+T$d*3 z%(aGwvfa)occScqepS}MysQYH(d;_2pLwZaE#hl`b!fHFU4Rh3tit>lX$i1 z_mP?NAF!;|Eu+U=8&fh!XcfO{|lvTE8+C<{>BF33rVHpCB>D5b{P1m^M~zcO^N?uPsFp zj^Mh5k(sX$=~c0k-Am)aT9k>Qb=-SY?YQc>j0ixy#6y)En#ZvelNl0YM=BMIi_SS@ zE9S&X0*3@f6P{_VwLdTZWRm5uR)5q_xfiCvif=@e!; z832Wz9%#S|eZ|hRBw8mAl3hDIwfV6`!l>aYJM}<|EkxnA_T9+GGAGiWVDw zMV-$Rj?fbd7Ytv-m9Xt>{SwnYR@b(1jq>wY3FRUmPInr{PJ_R9`?-6eMSpXxi+I42 zTW#2shwTC87D@T^Y{K-WnWW^EdX`0Tr^a~9Ei~%R+ADZ^+L>z~;lW9NRmuBDb&+zS zFl~wazETZ|!%x}bIy$SLN~B3s5;RR$3x@bKuTF#-nK2UVRDPfVkbLCe&)iR7CVQG zm#ol3PO6sGKg}YFAB}SBP?cLX5wzv|QxAQhM;(qsq6beV+Tge;B7e@KQ9A~F%2RL) z-3ia=Uc>t+%tL8q;S6{=INUK`c)1_ppIQ)}vF<)0-f$L`G~IB>zQ2=FnA2lG^e7Q5 zr)WlS+=S2zM3kSYsj2xeVrIQ*^vvIG@vp+BgM$90u7rOd{wYgC)v>P|f) zNt|8UHB8CL(~5_mAa98Iuz$5pR{0|XMSDiynfEo#iFH=rk1eP9DF>_`7Ofa=%{QM{ zd$&Gv$W#K60$j=rDUp|o&P`YB3T55l%$OrB5Q4YGGFXcx^M4YwL9o6nIT-4DFC@3! z2XE|wtGjf?I}Tg13JTPMW?fWA3^FoPA<3GU^kp)IDFbAJ+A(pW4+-BFf=k#}g< zBqx)!iQ@Zvnxqn6lJD_^R)rMAwu^{W}|GCdG?x(Gs}%2R9Y(H4j5xA&o*u98E9hE#zv$4GBv2<#xQUzsM?xq z^*9nbICF#+&#uCBv!49pZlrp?RHA)%~Ex^tD@dYCy83w^A46QArnu%n_ zNu90Z09_6a)RlHKP&!TAe+pGVDoNHL4>c?X}* z+kfzqVYw1h_}~X^A-$YI%WFg(^t(5iB?PeG62%E3krB!qc=pEE7m~tcMQidpX`*z57x z!o{p3L_g&AD(?=)*0c zMrXW`;2JNuJk#rJrAo!p%~ZcX$A8pbH#p=E|x^waUex5%J>HNw_=gLo{3kr{tHu8${dWV^W18Qu}Lw{Ly6z6CI zfZ#i^d7WK;@l)Jm?Wc~>e>&gs_08RdGVCM>FWgu-PF>qpVh~Zok`PhKqJLW?I4|7l z(zO@c3r{oI3vLE1;LID_Ny-z$jveUJylf`RSgrbQ#57zX-oG2}ky53C7ta(@x_BGU zN1hHqtaRSt#b1rU(w^A`FSIPi3n)r)0d99^)#uhj6?D|b%Wjxp{E|gF!l-=9y2o9n zjZmXxgHt+&IF$(~!fX4PmVeHgFL0NlUyj8%{;<%W$fec%m8@pAGZ3xGxioUA=oXj*9(m(DJ z9Y~MhdPKp?3s^JF$A->Oqakfe{-~eyj80bxXQIBd?m%(Wi&;rU+dOFe2z#g+`dPZ8 zqyv-BRN)GKmJ9Rk>VM!k_71s$6i@tJRLU1<;X|fq-o6=aF6aE-*b9Wn9r}sF@nUyy z{IaS9qE)U#>&XOMM9Sku(92Wg?W{WJc=4z5JYyEjvRv}(0?{H|56oi166N)DrMV& zRt(TeB%%(BmI|~fC{^vnuqkxdy~b_sTn)dX#}E^r@2MW{3+3T_h%LCCA!Ydrey%-c z9X&ar@S>|w1;2Wh6c)?ICB+uYmOmM=ys6ZjemPlIy?c@L0-y<}Se7u8z;6sJ0iJdS zA?6*rLFC+LmVYB?IZ?}JUUpI5>L<1Ecy4fR0u58pM&0^i`@isHU+a_}MWh1aMp&h9 zXN5=emHnHzlFx#<&jo`eV}L4l#Yd$|HM(~pF+iV5 z$sxCGkb(`Ci@>nXJJx+v&JW1>9PQpZDlA7DIE$G)oqyG?+w6&U!Y1V&1DN|_8lSjP znul&u> z?&OQyB2@#i4lQM}E$_qfpO~o~$o{FJM)3xRyvNDoQ9y*d(UG1Ijb1aXm4fcu_zj8W zv@Y6N{D<5L``I(xbIJ3HQ5&TUaX$K9pCqvLO;u!9f#2ZXcFz&hGg|t@h(P7d`jDMR zwpX|idR9TCDj3?nYQo#V>NkXUQPUl-F6KeS$`6bZ&)A!uV>Kj#NUS8k;$4m@RITEb zJ_m%W7b0+YjWQKPG;NH|bUbP1#2HBZPI9gZlWrw%qVHTkRU8VlCD5Z|5+b1}5a(WNyLQUC4EIrdMK+_`6W&-HSVM2bVG97zHn z3BfO<$U)t4gX_+#=4S^sOEg|`Ld4l(;vN}`u2?cZdJEU_>H+tqxI?Ap=+RgqQ=5o~ zLe&SS0;^E&)CG2Bq&({qw-$z^atmK%DEVFlMe+i_F~XInJ!Pw&w*zgeyH$%cCb)wN zt!LtedZL0)q%K<3@h`S*6m5z$`AOe4Ei&@$#k6ws3NabvXiJREHL84%zqjtos2&^- zW^C@t+cYW#)04{}F3Md(6s$S({B!Snp9GC1vK~O;Un|4Cvtzo+jO7iKY7zLC zyGBl(dc4O$7hdk}5*`Kx+$C&ffKMpZxumwIgEPBAQ`3#d>V!d2698AA)n5eaDApw}0QG zxnoE1K}D^GyF=VE#JJet zH*go&%#gp)2U!r9-gRt`27(?)@`fQyNb;XbAf;PewKA zw&VCpgLcL&krq^gMAVd;h~&|O$n!s+HpV{NIvdK@zE+uhl^Ean(B8o*$rWfhKs}&L zsT>gdWe31MMZmNqtx0Vi@G4r;swlvsL13wRfl$79V%VWQS5~QL&^B47=KHevLwEJ? zK;+ci+P(jIMOfYDT!iLf?VCoQ8tN*;{-s@OG8fYVM^+8Y|Hu{15<7>T5t(B&1;%plkQ?%f3SN4WmLa&13_5& zC^IUUZBt0(i#f!0Adno4|Gr3+hN>5)LR!iSb0BRmMx#nDB(z-Y3+yiq#3A9EG~!9< zEmIaxVAOK68jwU5w#1np-Xi$0ta>a+}uV`nK z9Sx!&6)Y}v@w=eF&?Ux|n#8?st($=T z=ksq*uM(VcM(~4^6TilpbEAOFq35*iy7g~QtJg!DHdh+m26N7@63_du+i?bKCpXUW zsC8$sODjpOV3Bt+Y&Trt8|2@qB#HBVXZ_*+|IR&L>#kp~tDw`bj?kNdT^9jy35rsw zC*S)x=*b2R_^R+3&Ev~slM#F*%w)*;htmWHSS?^#2B1ae707?Ym7mklYXpCI+3~CW z#tkAv+~e6Fa}l6W9Y+Z8e*XJ!v94p$Mc}9&>*NMJF=6GpRTB7Lmcs)4fxE<>h9}le zJh!ZG2|WAt5CF0jhvJl4h@RUY{kp?$(B6ONYWO4*%()MsdCeMn>|R4F=-fb6|D_K#`~G|^aZ)4&2a$t(p~ z(}j~+7ycRr?P^aae-T-)gg@pitPAPZ<{bN0{ikud=_qS;XXpIB*&Ar?U%>|>ltR{R z6;kR+i|p^!JfAx8rIEPs`^-efC3g^3wrwz z2@)t`@u)A-jPwBlQ0%U=O&s$BvC}PD$foI_%9&R@s!UjvRbYj~R{~%Ppy-?}@geZP zTDnl-mivN2;s-)e1%(8oMnK**q-9F*ouQ=W&tjJve8vNX2Nq0Wa)Om*wM3dpD^no) zW-K{y=ok2k!LA(Gl!6y*F^7?si__(gBaQn)S2oI)+TmRJxNRX;d)Tv z*%|h!K_I*|*a^l}l?zwAmut^zuIQ$))G zx-eYY6ZXkLQ!1aR^Vf!Sk8_gGD!W8?&e>zWUc6Is9$%j^c%xP~DRpP18n#2O@e-x_ zY3gJ$;nHE_RP6R{eDZG+Ledi)9P_)ZYC+v#y{=zR!04X|rGIBSrT>vQNBT)#bZ0cK zoO9mkTYr&wXgFna6rRpr%D+Dz=;Yqg>)4=aNf|hYA`)f$QY{*XQNFu|$uKeyvAtJc}54B4(DRY%TG3ApJg zom?E8w`bkXpDsOD6BsAPw#cXt`?xTCNMHl=h{bu7RVJPGvejAUpc zD(DtO3jMuX4s#OGTMm>`naiE*JTz=lzB|ZQ)()S#rYn0aVv0vhd|O}xEJ#t>XUn5D zd%u@PmuqSD2Z|@_W)2V|I^jVW<}I`x^~PMs(upJxD!P>?R%_U=%w!@MqZB~cfPScQ zT>*WFXro3nFCn37d5oT||0Mwt*}0oohCCl@H;T3DR=l?{F$2Oy^Ei5E0}1J$7^=Gz zyd==K)Yk|Xm+Ic>T2igpkrLe}DeK=v@2;lj=L-1=kSp~nP~o0m4DWsgve1tiTLPoj z>LU{|2ZL*+Y%w(!n@^jPabkOC8>gXsV0NAlb+t=P^{0zUh8Y^>che9J5-DFqWwSOf zQ)Y_c#L;N=`;=)$`tfOALY_`J|3I$yLpv|h;+7#xq!*Ei`R<&Ou6>W;t>k~j{)~bx z6DYPQD-aoikfRa!5`2qP<&wci2rS!WA1o!`di)j?ih=6)NDM}wt>~~rB?W+>N2ABKZshm1!>!X^tt+o(!QFmW{u zSX6~b6Y`0Ke6r#Lt}#WAjqr4%a!a#93yM9-H0Pp`x;`aRo{#aEvfh{uOih@ZZz{E5 zxsD?xG4NwdJFhB%8H1g#HE^rQuA7gV6JY~e(5tXy*8|TOlNm(Uj$dtINkNwL6V)S# z6&n6%`9#lHC|A&4X$(}*kTx~=i9;7aZ5+k4EwB0`HcdNg3-vj7G?cAOK7YbZU+1^v z)esw`hT*#n=(qK3+6RbNKo}b8P^qY8Ht~~ zBjrXP86c?gpN3wLpEBorul5OMDxIC7`^U11kHdP)_XMH^yZo3js2;pUzm%t>QLPyo z`1pIS%he3BD>xLZfD3Y42zx+9#++V5MQVBpqrc%orVM48+|!aJ3UT9wsrFeR=hn#k zt2v47W{ZukTc}NZU>hB+q!+GU#lo<)b{jAc4ynvjff*)84gnJj`mY3Y#s~uu+wcXz84@WOjE z-fClly$!8LJYqQ;>wc=L5%^n4O)1`_O$kA(QY=={t{P(}{*l-+8GU1$E$S+Cy;|cZ z8pDWM8)Ir$Hq>=cW2CI_{5Wpo9iw@qEe@b^ZTS>Q;6PZfVAtI>-eX2^GKfu-6PjL# z2h??`hhU~`5$ps6Q03XwTB%}{xHzW>L%X+(s_qlkKP8_F?R+gu=k>Yz9LBBSc zlPo_!#w;dI-5S#)=~vjn4h5HzrigPgci)$F2pOo8Vz9nVDivuEPlMR8>PZ!AnRxvT z%47>Yl{8h03U{3#eqFOl`rLPxYcV802RZZ;`)QvGvNtdl(`L;9f>05#Trfsd!pY3m zdj~}$W~Gj8ppyc7FQx29g;@uF9MrFJ#15WkL&J6DubD+cOwi9EP=ejM-ad*t#H#Zh zr~|~DTfGl*_8S!Wfjpm$fwI35i5Q^ZKAidEGaa$e_|Z&F#DvnpZfLEK}SI2 z3~6}9Q}w5ru#}+AN}{1=YFoUk2i@zE>-~CWHNr*kE)(pq+rp1qG%p!9(*zfMxwCI# z((ZqC1f}j;g1*YKf`rJaRYi4~J$JJZ)*mPRt#^!>2Cbc6YA0?2*HO{PljFij;gm6P zv6liNIl0Rhi{CPF^Cm)>wu9xZ<(XhOBj;msCGZiTM+v-0P1>gW{lH&z3F|jQcqz?Y z-cxYEH28UOeMC4nmV4ReU%>r_K&K~-1oZAjadX50x-pW-X?(H#!sS=y2QX=SO?nG^ zvxT=q=eo|;Z=8#_-yFMiAH`1OU-#dxn9oHqe#?Fk%9kgng*?}H5z-fGX{z}O6WClBn3699VhQW5*z;%{zcfV#o-nN4Q>$<5a3olopT;psk_(@4x zNy$}A9V?*jzl{TqGIns{BC?kO=}T?IUp~gIUIQ%!FBprhBeR#64_&0O3Vs?{L#{lx zVUcz3b|#|weM>gz!?@r>U;0Gc{66y z+cj&y>zvY7>9@6ZHX_(|-o3i48tHC#MxgUXqQYseI4CI6%QQ!I_76eyK9Il)-?|#N z!Md`N$A!QDR9Dib$cY{KT>|4=a20vO0FK5o z!89Z&@+b>l@4aRFbu>d4@J%M3y6%o;Z>pjKA3Q7bg~3t62onb43Qt)>g6rlfAsb?uI~_h}VqD`+-?tZX z*y>tXOp)waXVS1VZ>@Y65vqk-7ch)r$Z=QpkhlswV4XHY%&`hldmSM4JzsmJ{p^Im z<>%x1V$_mgIF7n8l@$XrhN|?Jl!wGD52JLzABT=j^3X-IXt8`pb7>NQ0h>G$KgI=u z^=93ubeDxQ{&AYTO17=6f36C7R1K`f)u~#fC^*-WF_6<2S!{{w=x?gx>F6fsm)P5) z^<1J-lfP2_&YC95%999ZmJYbRbO&XOdx;N`8{-^PEbOYm1N?+Up}E6Q;orA+$KSY) zmj`Z7-P+owMpkv!0`6~@V{TWkpUyzr=Fb(GtLNR%PYm5lXL2TsM%7&F;l@R$O(v{4 z|GqUHw>=Y2mCv3yHLm+t_86_NuJ`zsKltXqDe1R#!iiU=94Ho=R@!{^Rt<`pI$w!5 zGITZv`tDwwts2(?pp`_1%d1BeX4Nsa4#um0k08g{x2x3%;7SMTtg-W;wgyTc|7+G} zw*dNM&*xR$y2abE>OhW$3f?NV6IJnl_nk*aRqL5>si}Fx>JlIOEKz^huuj{v z*hORN_e$pg6lV_pwzm1{S256-T~B>yOIO5oO2@g`-m&c%X%4JDY0K&9>J3t5rA2$y#H#a{YKd&G=mk<*d7ZVExr>e7s zjG4P7EwhXW4;T0UBWd!s2aEtT$V&E+G7mpDux$`Z zX{%#w2lw_&53lu*MIExSg@;6hz=InYJ?)}%xCJ03TwYzaoiA1si*6I+e!jAWR=oJX zIVQ1vv4Vb%RW90YZWB2c`L~RyE1+v*hziWVsN~j`wW2!XujRl(y?aM z@-%H)4i2YQROsE;^00=q{rFge;(&sZGI?zHt6PAh1Gv+e0kS*- z6|OtzCh*IK&;bh5A4L}^1|{;dcINVH!n5$r`$NhDb=$Uc%rL!XX%VGXRT1p*K=Dc3 z@$ZrUWWM`dol_+O1EoEi?JZdw)R9Xk_4vpAVG>>K9QnKPBCe zura>{#yMbfEvZ8eea{9sufF(2fipcbzfD6oYKI@3*S{s9pq7w9^Yn*%B8OFmlu^+2 zX5G8R8ZhtzM|SCc=I&}kbp4bWBvn4Cxc$Y*CoxSo`#)^P>y(2YHswEL|6k0K%+ecv zfgZQlCB@wSat-j{=02@YtN_}1kbmQBU~F`>AbDB-Uz*MS2k3vp%MyUFi>>3JzPl(2 zJ&Pai|50L=Qtd1jQ~Kh&ZL95ZlUBk#*Wd%*{>5qRo2mYKF*Q?owC|f+%*QhCojJc_ zHl?MQl(sH5W^?tHRZ`1n^6o!NT{U&jWEXowgBRM1b7FGUIU|dYEr$n3>3QP@@!SG4 z WTF + init_number_vector log_fbar(1,nfleet,f_phz); ///> Male mean fishing mortality + init_vector_vector log_fdev(1,nfleet,1,nFparams,f_phz); ///> Male f devs + init_number_vector log_foff(1,nfleet,foff_phz); ///> Female F offset to Male F + init_vector_vector log_fdov(1,nfleet,1,nYparams,foff_phz); ///> Female F offset to Male F // Recruitment deviation parameters - init_bounded_dev_vector rec_ini(1,nclass,-7.0,7.0,rdv_phz); ///> initial size devs - init_bounded_dev_vector rec_dev(syr+1,nyr,-7.0,7.0,rdv_phz); ///> recruitment deviations + init_bounded_dev_vector rec_ini(1,nclass,-7.0,7.0,rdv_phz); ///> initial size devs + init_bounded_dev_vector rec_dev(syr+1,nyr,-7.0,7.0,rdv_phz);///> recruitment deviations // Time-varying natural mortality rate devs. init_bounded_dev_vector m_dev(1,nMdev,-3.0,3.0,Mdev_phz); @@ -1010,6 +1010,7 @@ FUNCTION calc_fishing_mortality if(fhit(i,k)) { log_ftmp = log_fbar(k) + log_fdev(k,ik++); + if(yhit(i,k)) { log_ftmp += (h-1) * (log_foff(k) + log_fdov(k,yk++)); @@ -1021,7 +1022,7 @@ FUNCTION calc_fishing_mortality sel = exp(log_slx_capture(k)(h)(i)); ret = exp(log_slx_retaind(k)(h)(i)) * slx_nret(h,k); tmp = elem_prod(sel,ret + (1.0 - ret) * lambda); - + F(h)(i) += ft(k,h,i) * tmp; } } @@ -1138,7 +1139,7 @@ FUNCTION calc_size_transition_matrix At(l)(l,nclass) = At(l)(l,nclass) / sum(At(l)); } size_transition(h) = At; - + } @@ -1224,6 +1225,8 @@ FUNCTION calc_natural_mortality * @details \f$ S = exp(-Z) \f$ * @return NULL * + * ISSUE, for some reason the diagonal of S goes to NAN if linear growth model is used. + * Due to F. * */ FUNCTION calc_total_mortality @@ -1241,7 +1244,7 @@ FUNCTION calc_total_mortality S(h)(i)(l,l) = mfexp(-Z(h)(i)(l)); } } - + //COUT(F(h)); } @@ -2369,6 +2372,7 @@ REPORT_SECTION REPORT(pMoltInc); } REPORT(survey_q); + REPORT(P);