From 26efbb82982566d2a4d92c2c2a621849f53796da Mon Sep 17 00:00:00 2001 From: caes Date: Fri, 4 Aug 2017 00:26:00 -0400 Subject: [PATCH] moved magdziaz sed --- .../magdziarz1997}/components.csv | 0 .../magdziarz1997}/magdziarz1997.tab | 0 src/sed/magdziarz/.RData | Bin 33181 -> 0 bytes src/sed/magdziarz/.Rhistory | 146 ------------------ 4 files changed, 146 deletions(-) rename {src/sed/magdziarz => sed/magdziarz1997}/components.csv (100%) rename {reference => sed/magdziarz1997}/magdziarz1997.tab (100%) delete mode 100644 src/sed/magdziarz/.RData delete mode 100644 src/sed/magdziarz/.Rhistory diff --git a/src/sed/magdziarz/components.csv b/sed/magdziarz1997/components.csv similarity index 100% rename from src/sed/magdziarz/components.csv rename to sed/magdziarz1997/components.csv diff --git a/reference/magdziarz1997.tab b/sed/magdziarz1997/magdziarz1997.tab similarity index 100% rename from reference/magdziarz1997.tab rename to sed/magdziarz1997/magdziarz1997.tab diff --git a/src/sed/magdziarz/.RData b/src/sed/magdziarz/.RData deleted file mode 100644 index 9ff24376295a26a430b5a330db2971505af4a92d..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 33181 zcmV(wK}zJv-rxPb2S!Co zi;qx?5DNSXd{|VSYaW_e$6i!na${Exgm}d)k z+CXh_ocnEV38*iaPfD{6g2pAXU$tVFKqIT8OX__rXj)(F5nDpMq{` za+vALFVOweRC|L^40OLX&X{(S!M!iN9N47ZptpKs5KXiXdiyShUtYC?zVBWXeM15m zTx^K$i7o>am8yY%4GrJRLhorw4{!PEv0T$l(52sd=%dLU3QNO2`q5mxh_aZh+rf}M+I9>XpH*jvYya#V1F{VT6+1#ez(JmVA)_IV5(1uv31uD$@L zHlxIYTTI|gq@~5vo(j&A^^L@AE#SQN$%cI>7hKUll!iW;fve>AHVLe4aCLe+@}Py zBMJWDt3krJeGs5l`u;45C>DD8I81igWSxe{lP(WVg|a{->%E-I_!$soF#CZq!2qI?)#5zQ%tO?$ z_vkk7FNhxMosPzS2(hbgYf|pgLEMQkcQUbLh;#lH&r;U|ajp$H*K+Y7&WkMX{Nrqh zdq(<|F=-a!<@R&?ZG9l#D2<2c>?w#ZJ{&PmybcM1i`i-HpCQ5I^!vmtHb_Y5JrsXD z1&O&oey*MfhQ!Aa^`%{RAhE%r;PG1}NV;0V%0I9QNwP!SO!q_~NpoS7&i4T%#oTq> z9G!-wO@`5m&y$c`=7Gm~X$g{_i7MP)V1wj0>qEQBjgVpxHPyS;11XMk{^bqtAZ0y> z!Xx=Lq-=BLF?l9KYJq^m2i1B=U0;$4rHp|z^0mV#)H6u4vJ#y+pn|lQgHlrcwvc|g zuR}K>0MgM6ip}(bkbd*(VoQGpr0ZYRkkNez>2KMm#zrC_Lr1IA;^A+|SZGYpH5P(Q z^wP$NKM7s`1q( z$k9=HV%x9@IR{)eWomMei^sRea;pOJaQfb-zMqG@=s7$d=S9eSt|%AAeHQXBdN715 zNI||b`I3BW3go*tvigsNLw+1St=`%&S1qx>-Q^Kt2plHTf@4#9ciuHJ3 zyjUTD;u^K()_b3!R8wc1*^B{7Ct7k-Cv>1poI)coLk%jJx8e-WKy6_m-jLBC)RySqcDL1nhwdlosA|RF zVal=ucg7ApO#k(gCHF2o4*$Bq7dQ!zU+Zi5FAqTdIy*blrZCiRJ!YCtPlE=7<_OUW z8E7nR^$6#`1kVwfsH&Ta&{Y0Ig|F{9G-uAdcS>-F7nff7TwPa!mXn-y-`3rrC4#Fh zlFbrYdV8+DQ)PqJ=~3>;@lI&-@XDT07KXMI85V5saOkk8Q)J4ghz=Bo0(^0>E8^fx^DTLEd6jl?S z9^hC1g0(9`OdQV~VJ)Wn?jepQtOw=Ikks?S230DtXAUlGaFxmqpPGZ8ydLkAcB0^? z%j?IYS|spu{Dughi!=QE**Q1zLJ>AU7{ADSb_KS=8eM}OGGQxQ{&qllB5X4W3%y@A zg`L!=*KR1X!|si!^9pHtu=}uLDpk}P_F`^^Hh(+`zgF87E=6?0Z+w|6xhEQ6|Hgo4 z*S06@|M*i74!FoYCS|MP5VJRvXYz2E)@;J;B_n~Qs@MK_i$VgcYFAe&wpIeG#<8XK zcdWz-^mjK(NRT*TiFd=5Dplgd;Y7rbhm{i8#`@R7g%u>Q(@rzc+BHg?lzb*2`N>k^ zl>Dxo`Y5vmUU!_;8_s$OLU|X#`Rf@HM0ybf*yAT8NW`i}7Le-_Wc1FAGK3`(%TzBtsFL`oW zLaOMfx_2dHNr`tIWThqKy)kDUcUwZi<8S?MxqcNlx@20Z9u}f4O8q|xP-|8RTVxjvxFILEp10;iG=y35em$nBrL@1rSspk zN?4l5n>p0dNm!|2!o|}P)-r772LdLal9$(_Ay6gcNxaM~1VQpOX{l-m za&F1DQz?XCoQ{a`rAH9FP_=)x>H&l>Mpm>(heOCPdvTt~MF?B)8QFW-2w@vv^%Z67 zAR;soci2D!A{M`XxzVixksR?_0dD6Z(ilHYZSyxo8FueJdQuEgIXuI}v|JEXHsJEo z{Q#nK<5Y8>NkdEp8)0&=2*lL?{K+qQ6=KKyo$HMQATE2lQm{}N;?_1juf2Q$@tQPm zIt`T}es6T2R74mOo{a=w+h~DAqF(pWr}mIom9CuQa|)99hT8K!(Lz$M`J%6}I3%AL z7x+Yd6_Pjf@ZPI_h7_SaGm?k|Nb$m*%#2%tl;Gmd48_er&0!=Z~<*!{~-t&x@XY_>i~zy86zgL&$H66k-i+g#3xv9UkU6$lugAe}7c~ z3WQ~o7V+$$pzcTSRWm6@5lKw9} z6NVzF|0oc035s+r5JT~MP-J))S=_z@MebEea$n$ZPR4M@!DF$7*t*#42QbaC< za#~Py)t*LK@d^~-Cpb5JYlXs(-u{9#olw|vNoZ<`3JU8e#r0YApir-i_^0_BD3s0l zg>8XA;RPM^g@!;VXnOA{=lmQBVoZf>4^p5&Makb54U{1- z!*zDC5(n~(BXZQv`9Pl0v!v%UvXFP&e)A;azID}TLhiRe{UGGNi$_-kdq8f))zO~J z7|7XYdBXN14RU7EGOdQ~Am;-SQZ>B`ISNFZQTR_G=el#J`ZHO`ktk_Iy(fVj;<9r! z6P%EJXhZONFbA@Cq_^*$?}6-TKJ%}sm5|-UI)8`aJ!E^Ja3vxqA=~kbtc_|WWLqB< zNP}#V-+3Ruw?H;urQ-IRSp2D4u}VkwFII>@A|j~C z?pO;z#(=_yv&$S9#(rC|#es~$x2nnnoRBfZDoe-I3K{SD87@9+f{e@%?CA+&$cRG6 z1Nt0H%J}4B0kfb{LvADm2SM#z||3H*`A~Zd8OamqG@P zP^WQ)FJzn&?4C38h4j5ovgGlhkiN4rSf} zzSAZMiHG!;Z*?N15+J?BxE_t%!f@yv&aE{_Z%uqu5_t*I6SyoD)qg^|yH|Vez5%4) zx=&u9CenF&tPm*uRP4>T9($d`R2=fUU8v4rwbf zL{2*0koNsc+Pftz43}^^Meae`54oLG^>s*l7Oy|MF$-x88b7f!n;@;cCSIjb9MZgQ zo!XA;gftJ-7cAPRkY?0JTFYJqsfSDBgs0ab^>a82pAs*m26)L7c(y>Qx|QS=QGZB1 zD0eU%sBMaS2q49LAf)CYA*7HFkf$&1Kyuw`nG^FWBwyMoe`78UNrxf2UJpkg zar3Nz3#$qwjPG6%=c<5&t~iBn%{dVNs}i@|)EeT489VPYNkDvR3{y^z7{qf`&TKaL zK-^aoj-YmVh*Q@3GG8hMvA-T!Fbcnh*o8j>AH>vrw|(5+1~K~MrmQ%pAo@FVoM~$< zL=QCLq%Cbi^l8W5u%A2-W#+d&Qqluanx}o^PEtbTtmmq@$v8xexDM<&>_CKm=b}JU z9z^h@DZ#-{2>%eB$}#j2!bLOG{ly0$T<|V?Zp$wS`~AFh^3E=VP0ZB2nz#XBgK-%v z8nqBs`PxpV5wkvAL*12a6v8+??<+pwgV6Oa!_jZKA#@G>O4q6#LR0GuiNsYPG>5=q z>fJX8eOQ;4I6(!WEjAY$&sIQa7XyV}^B{zNNy_BB5C>tk(Y>=`0T6y|zxQ^^JBXl@ zXPkbi1(Eo^YylDB5Lp^%%EBoFQCWPwL`9qsolfmA6uky9Ifnu)jExZcpc#FpryJt# zpTtQ;Tp(^^jza#s0mO&TcG)IBgZMQf>dnWBknrmM&BKTNkjVez(twQ)Bz}4Lf#y^( zBt^R_m_E&gWJzw;MxpPJe5lGhC_V|PCkt-=Xv>7uNY*4JECkZ-D&g4n%46DXj}WG7 zhx9JIxRqUN$f%^*zbwxOnX^jI6+Me!x zM*`W=x1*lL;6iq+&fybHj1EBcSj-)iLY9n&@TY(x$Rs>{x#{Nuq<^+5XbFD}X^!+l zkq_h`C1dj$`-C4!t2istjOhRYB?w)W*f(GB7lDr7b8aO23KFlxE0lSQ@)-4=Xu)*<0U-zp9i%(4B;sOF- zy7!3qbYu}2+r(Fsy*GsWdi2FIF1TQzd`qztHv;ahFJ73Gl)`IjT{P{*_%N4bxs*@;TvRKl=wL%JKBzLd`0cXt`WRV^5h7Dh`@X7bPeGvEbzG&pV5_i5`5Vnju8&u z1YgyFd*|w6!1qqC&L=BX@H^q*@DF_sbn`kg4hQN;-*`0g65JcLRca=5;f?^hoYeKI>FyFh+ zbJ9`}yeaTzYQ7ah5aSe`55f@2`-|bjvLS@J+n-__sfN%}qAkPM3K06h_)tW92Erbg z1-eoiWBP{+A6awsAZ#kAV+*N=@VQ1jYPLrZ!NjWPG@1<&iCDV~l2s6qLvue({XIkm zdsST}xCD`rKcbcdZ6I>{VyWGF9Yj&rl5-oZKvc8r&q&Ru5Va^|aYZW67(H4hD>$Ww>dPX^~ZqVGc-+PC?ZzZ}HL#x=-?5kOq&1)dyE zZ-{&8rtJ7+6ykdEGi!d{f_VI--R&|-h`$#odHzra;`Q-;RNJ;7-pUN0`TTbb6?p_a z`5<0i?dZxwyxc1PC_OF2OQW|Rsun}McuSs5?i$3?vtbm#C5S&?YiTWH1o5ZMTcR=T zF|M8XaoM+Yh^xFV68;#|uLqWpDNfvj*ooq|{ZjR0E({RU^n*b7 zDHB9DDp_>V2t#zxcalYZVu3B8syg zjYj)HI4gzwOZ=A*R*;!DI_L=@n^uRopPL})+|z4Y7YiT|ty>*37YhM=XJ0La(u4oL zQ3k(*BlzW;XO+B^1;1?KjK1v06>M|{?pxBv0n{Cz#&?{mx=xp*sqQ=5HSXT z9RabNtd%F&%3^0Po>l;B9|^h<%6c$A@Nq7A%L8VsUEaeE8(@ZgC26(U8t#ktHY<2> zgW(deJRTqb!;$5T4gCo)xL8Pjt8`=Oo3$sNEnS7XIElg) zO&p+szCT^TzXs}k&kDTTdq6$C7@O0-9aI}cmT}%%fvUfb{=tR`+|F0Yi|plw+urrD z)Guq`hLpue)LuKtVyo#asIY)cE;b7pZyrdOCz89)bb-`6VL2j>yCA8gN|F2S6kLDD zaDxfE0j{TH^C|Sy!*%YGUv|ywAR$j#(8zoUs4FyW`$3B!;^smxu)P67!QKH=**`$w zInHXqX%FCEYWeo7Z3%b<^To|>H3Cn6Z$b`5Ca`UOv+%!w3#=4gHe#8)z*v5%k0J6Y z&<}eXJ{~Us`bSC5w3m5-hW2M;HYyw_X@4m{Fue~Hj&l{W_z7^%ux5H1Zxl#Bk>}#x z4Fi&<=MQt}Er5iZ)~V`oEfAgi`D3r-BoJICTuYIig;NjCK0QxN1gCgpTuo$R0f&u} z`GzPFU~{nYlTnny3A4G9wW|flL53qMPI4OZYo?o(*q0FbNr=OFwj~T%<0XnKsNh6? zusu2PPF)$9+PwTNm*FNdab1xjgrXJs;&{tVNzNJR{lW6^O~*^*L%R%4<)Q%c{(#;V zo$v!`ro_Lfv%iisY0BW3#g8FXU-1mwma34-f>P$DIdY^jaL?XJwiwBjJt$vU7eO)@ z-gX9X)FT;q?dK{yY>|la`N7}mG)P#Eq(Q%rBjPUVeqrrhIO2X`3g1C%9IbsLo^&6FMpl0MpW>nOKHdH5oMbEk0gj7A}wAH- zldSI%#HCHiqbUs$AnzjeBmIooo+MhE+7cCDvRpYdJ{sr z?Qz9`K@lPP%xqfvV-Y#o)LU?&TmpSy!BAd|VnpxWnJZZL-9xWeJbqhwc@aH-pw-LC zWrv=|SKC;07e^1f)Szu&aic$7IjMgZiyPf{LOWt5R1V$6q=)M;L62_7c_F;g%ZP5f zb4$y}&Kv#co*A3O4?}dF&eOBEAJ(D^8y~5P?(?AwQP|IDpLwA3ER$L@nZ?ivMZf%M ztBKGF{7I`5#P#S{QUeBgxiECV@9_@b8zJa`I4P3$2n)2g5!sw*mIB&-jg$JMl^EK3 zx2fA~1BEt}ZrQ8i!$<4r7Jl|g*g~s5lYHi9`3SA>o&4_Q{Sma(Rbg|vS|2nzFkL&g z?kZY{m4`vPUkc5>wl#fSWd+S9*HTi{#g3*5$@nycYm6o~#3@g)+d~s&*ljlj_@hrG zwc8Dp`k}FK+~YLTT2Q}Fzv?8I2uE%CSInJzHj7$6z^0eUDn%^@X#3@Si$l#?t~SUu zpF_>e>^*KRpg?^I@1I(VphS%xD(27d#G!^tcNI4j%u#~}4y~1hDyZHhJo?3H1oe*f z3%wlUchp;-n29obVN{og&n?YlAJj|1stb#sYEX@z0xzy)cA%c}%vryl8$>-!Jg7@$ z=R{T9T{*b-_!p}DH7+sZW)`ZDMgHBjD_W?$2#pyHo13WI+r}TQQu9!0YtcUrUv;2T zcDWAFs1{U$sji8O85=5g7U!HPcPA?L6?4Wj(W|I%C-ppUlUY<~`?q(896_i+f-bD} ziB^=KO&hl3i*1xo!kdVvR49}iDFGP?iWy}V;2mgHIEB(=E6x49cLF7wI)-0Q(%v^0`xuUI_Rx1-tocK}i_{>ESQUrY{#8#;v8+ed`;7CV zVv!8Jbi4uSVxDD_ZojU)60=}=ksUy@DyH@F`h0X1gP1%peXqXpSxmG&a4e}eM2zi5 zS}cCx88NEm7QPWqDKX+9+W05C3S!vPvQ%fJ%|v&Yo}tVUP0{a@<;p!#?V__c(k^Ac z#S{IaT{rV1wpR2L>dq?%Ahb)Vw+4AF}gshnMVb!$vCPgeCg?R7cP46-S^p#(|MltBEZ zIs+=v1iX7a@wdE1W9ixo3iR_t!&bfZbN7ZtgTs_|+mu{H1Cpm*cH`Pa{Z4=1bmUO!NBoS)Tu=7{O_MUqV}I?SJoE9MeS103mBh&BWhLr+HB9WNYw24 ziUvCYxv24-MRnfJPonp~vM-+aN-1i|F$1kdz!}Y-2SxfY8Y{1PW_DWIBx230( zDey%#7<`!T(o2e}($Fv)T#WxOellQ5Z=LkUu(y<%b`ZSH391JyPk^^w2a$eUBlz8t zBml})@K@EcD$DKq*H8Y}PyTTNT$(jcNM^8&(3xYq5 zTl{JDT<|}OO?TaP9=u;g5|dra2e0JYlkz+mUi?Asthojr-_CTKl0F3Yr(GOZH&wyS zll_(YlXKukq&KsCNCmF1Cu<)kN`veDnaa6B5pedE~$Bu-9Hr$c7jv03w;S-tE+XRVjGj9egub&FXT0IXY74PiF&vk*(bzwXf zv;i1imaeRQtq1zD5yz84Ab2K6fx(oF@AKxHDD<}~LRD1T+A z5g(9*TP4*MJ59-;#PV=dY@ZNrwv1>p8@vPQ8x%4GYu`W;kNRRT69HU1QAqpMiw7j! z)hbH)2?2R=;@g&J5}=HXs{QT;gSh+lrI4;@5N*bhDYmTuk(e2OiBIJq@Tzph^|ruahgx~Er~quYp1R$Bdka``8R_P4NWfLs1j#KzDWJ#hsm}IH z2D;aH)L$uTfU-1$ZZEU}DAn}J9+KPx(i>W3y3uVw^6frX2V)12eEoXv)mL6x zekBTc%wqMgFLDCjB|bwYob#ojOEPlN(C_b3`V42y~iv@{zHDt+BAj zRLIE3+=XA6Cy`-2weDx5mPmh-o6jBNCrEE4RHPVk+*)?0kU|g%KTpDg6!1R!rFqbSWK$|gee4TDGBg@| z3#9{)^l^66{hoXz;VFk{!);w8{t4yA?G7a*BJF$p#AYuNaXtC>y5$xUhUT-pgUUt% z_~wT-R80~8yu}c2AxcHa3M0ukMDcN^QwYBsBK^+n)SRCaB5_$elaXr}0hzD@TrD?5 z1Q)NpsN*>zwAkDf!Fdl6v}y}{8=8mkvqd=}XXp{$6k7&5J~4#z_T#SASK0`h3bmUN zz7fK@U&YJLfrT)>A8yq~pF|i2)VtOm_99oksj(=?Gfrsi!Z0$R*=i;)K3utHiSZSNDRf|hLH1JDX1|nM$TPH+sY(sMMy6cs4y3H zAZMQAIZ7}eCa`}y6kG3o5+){ML>4>DHJjW`jmvb(P6I_hY(^))thnaUeV6T(mEyu!^H`ziu- zYVn}~WB5CC9B0$@!*D-z$ie5NAnXFPpN1(z(l0Bt%U$m~*X}4{=4|7gll|prhyyBRU@5fCA@HrRI^iPsVHKrAh9SknqT#SLJ=#Pe~$v&XVewV zVKi(1E^`Nc)oRu`b1npZIZIfQ;a3QnY)0|pV>B0$QTOA)Y)NkQ&&DtN;vmWF$Z=8X8eR+3&J60Ff z+T#dzF;7v?)D5B^i&&x_3v<3))V__XrMaJB_p=vO{#c$i*ijo*e3OFP1lUpeg}9cP zUo}zLB93eo^_Hk~2C0g%SHDon1w~d1V-u)^>;Ugj6h0~rM?t=wq!<;I7kz7Xa1#}t z!b8n$eGV1!rBq!j$psZ;yXsaR{{$5vC0T9Wa0=z8a`RhpONF|N zISUVelKiJmHve$lx}k*;!+t}w&&C3bI(#DhityxsiIOt5{i*xlJZQpFwv+|VUvG^GDO?BFOFWvBiCW;AH{5^TL<(G2o(0ly zl7Sn~DJawxV8ws2HFMvPR zrlTIUEcko9eUPx(1OZp-XyrDWAs~w&(w4vt0^Y>b;LNr|;L%ez_`DEU_4RV~m#YxC zedW}rN(6#Vc#k+0|A62>kHcMs;I)+6IwCa)IjK|G=YJhSJ;P@%u3ACp6QV>D_lFR6 zq45suwd)Xe^|cN|$0-P><>=%sy9h_0=JCvWLU?t;Eqd>h5dM{Bd_hYJBFdIz1X7u7#j7baP`p&9zY8^t{p&EX2&r^t}CH)RWIuIXbdY#Ly5fY9b-+KHG5|p1Xh-Pv? zg0*?%o;uX&djQpi=Z@4niF zlt7d|cJ2j8iQ1xO;o5-I6L%cXQJ;gmyvEkVdGe|9NT=(#}7a=W?BfG>)Reh9N~rQ>$WqNx=?j26npo zoTiZGAgqoxtqy6=?AgxZmG1R#d z((0R2K1?4#+S|+PMei#hZDUUarfeW>(40*?0%q|*^J@PMU8Jq>koCW_XX316TAg=dq)izb_ecK012U zUx(X6+AE^x6~-AD{&}DjL-)&0Jy94&qb)j))+cV5T4MI|Dql~(p%l_;{~Qp~lC6kr z)f6Dj)91IIK`f+Mk^d$M6^1mOm()-5mLW~u;I2_cF{FtZFpQrXg|y>&RoayX4)&o< zNV_n)Kp9XCX?W9n#J{#7b?1$L+S(^b-QIJg^@4CQ|a-g5;XZlV1i2A=w=(>&sdLB)b)qpZUZ9$vU3w z7sC-qmU>XVr}_huj~@3)?T2Jiwl~$YcOe;@Z-SYs2$KFhK6(a{hI7s498DpqEA*T7 ziX~x)Q-2PL)AvbbmFXbyi5j8gtx3$hG5IQ)Y6~Rl z_b}Me^+F<=RYNRS1QKcMlx>5kAdw`2#`F9QNa)295xDvj66!UVoVSS(qW= z?utpC93iGZvymMWE`o%Uyk{1qKSKQBB2UDWEW|Hr(%m{m2l3e(nbtKt5TEgkogf_- z;sb^JFijZZ#a2@W?VdrraQQQLB0h-Q=m@)0wFz;iuT2DBVERwN+2l*r!w{<{F)VC! z31UW&w-k%a5JR;qNLct9qO}s?Oz;hea!+$m;n;vk-wVQOYzz>Q`*y1BL>h!&xb&ex zPZPpsDv75H`XEG8VdaO86a-zARq7@gguoe-%L&5~;GfL#-puDB_|sP?Z3JS*t02J| zQFCnYCEuCwI>8Cv`aBel{DjqY{T3=f(MuN$- zBuL9(oGS6iBH-oGyYUp_Y3CK!{-kd^x zfol(bagn}w4-)9}cE+EE0F6iY<%Xjbh@sk8xTqXKxPI#W>1qw&*}ntu+dF6KTY$W@H-DPJ)rOSd|JscmsX^iypG%~5O(Jgi z={kH_8i?Vek-i70pfI*k6K|)AthUl#X1 zs;ZX1BaVXvm7?zUL;oBL%2!iiyQYZ;WoA~_8Z*C&x=pL|6H!t`iJT^Ax;PYsqHnb> z{Y;LHB9g(0?WX@GzBPQ~<@>m&;?o~a?Udp7i}!L(-$TU%@IOO>kLK9*8@ih6e35(@_nFII|+Z>;%JJ~I^;@b$-+ zNoWw~I&(nlsktN0;wF#k&hQkc`|+mP#92|CT8M1v=llin^OnyPCt7X9N$oYL28mO| z&lq!)zmmEvj(0_HkX%Di9J_z5xh;N8?4UgVe7<46*v{s)_-=9=v7agxPOqMrimmG_ z51e%C5L;PMH@!U1EVghjME>MOZL#U0jqu&3=VFsUX47Ef8h)KWmo31its>Sm=UWtTLlF1IejvJd;`JY!FP{* z_dzIn_J-j4DhRdaRK+rmgRuNSpeykQ5RuU)*$pcPk$_NYu7f8a%Gb>N@jeNNixl{Q z>llbjHg~L_*#T6?1&;o}F+iIsDHq0G1oXfIfz_rvfJlx<-+JW&64W(jhKpk$aU)wT zB+3RPRO;^P4b*|e*7BQ&qdRb|TwjQ`&u|P|(2ViMn|f6w7Zb^C^>o;%~K&&(Af3Qp6$2=kV*Gl;suwto0ozp9;{S z^AG@K!aP+@y-rZ}by72aI0ee(pT`K)U&8IXrrLx0G;rH5)9$OdD5#K7GLc4QfC?X} zGxdElP)Xp=Xr9S#2Cpc}wcugDb& z_k6fNKkHS5d*h``Df9P1k2yLztDYP5=^`bFW0OH&=&P9UV+GL9di8qSj}7#T|NOxj z7?47my|M=wTt%;!53+!P;f0-x8ZW>g<#Q1`B_0@j7A+v+Nd|*e*;H=LKE5RB|Y2G-*+@scYsbC~hds5*P9JYWrs+O&s+#?ONB8LKNq zr~SdWFu1bhbQ&0!PU%$cJAv`o^8kgs6)^r1HSPUE7L338aX#(l0pnk4PWh(kU?S!u zO(&}eCNka&jifi-o$)p(NX2F!% z=MAC912AQo^SRn}2TWD(`)VkdfvG{Jr#*2rn3_&Be#;aCQ}2soR~QArw5W~Zxz=4U zU3h7q;YtjqKa5q+|4s$d&5HP#Fexy*^u4_;`XrcLv8<%&V*xYd(gkeKQ(y){`K7h^ zU}ied$9+Z@%sgsu=+kUu1u9`~*v$gpO}4gS(a4Da>M;&jYIkA1*HOGO#9Ms-e&d1?&3s(aL<-pF%t4x;H8tlB48!Fh(fc>fZ4}ATmV1KqYYf1MZ*b`?keGv@+``n9> z*GyT#zBXia)8#BUNCSfw!%1+Ev3$c4PzeqN?|u^FVuQn{5&i0AQ*b1SH1n~^1V?$r z-d8vr;AlP0x@?#Pj()XCp_7x~I6=;=rr8NjWbgDXZ`FVk=QXTeAq1Q<>asE_=D}&; z9Lh@?!`$<)B(Y(&fyu<-6r9tTVHY#u# zEgkG;*a269*3Ihcj^L_2=P6&H53V1BrVdac;CA+`e)94>xP^*PUb!0wZq4CDIL>#$ zt#5^1v%MbNSspP*@a2H}EfmWwH7{_#GipWks0BQ@_R}4;&w+>P?rmdwTkxdtqEyo> z2hYoI^e5Yx!4u7?d&)c?JnK(zson|$&ljl`Lnj}C=e7fCT=N-tRi2AzM7fy1}XoM_pEAZiQw zE7&~{q{&NuH6qvaa%!QoiN+k=o|B1`({HxTmSU{F9U6GHk-CH!{d{(V0A-{+J6mwZ0iI~J=`nUdP})75LO#0T$DpC6?VX>*809-j}9{c!=$50)W{ZSguT+8Cmx1rcT0R)}`l zkQ*7@glPYmZ??MtF)q}@GAjtgzF!J{6Z9Hl$D=NVt*?nYM$>(f1N@=kOtm;FPd}gA-)E?jhUPZG@}?POHuQYRGKui9skc$5a6KTOqc)m}vK{h9)#OJ| z8<3}_Fs@GM1bNCT{ZCy>AWtCu?5KzxhP-+|b#x(*Z(+ixRuA&%1hB;^YazF1-1l}z z2jn)AXqr$OLT=zo(hvLNUOJnSXoC$S{mP9VdQG2=h;X^It+^kXm z?nDMT@^yQSw)Bv5;a2IF!(Whd68ptFt_zU8r$MzsTLIbQx@l8JF_2wQFjPC}3fc0r zeH^;%kbV6$`o;_!WTW+Ta7n}<`y3mq{4HO|J}WB`qPYQCbJjW^8=E0(dZ<#7RtmBT zi<2I(Dq#9;#)(s7IFJ=VWoPF)3Rxkf=PhQ>LsmenB5kD>WO=Ald|$#WyQWW0+2euOEM_}3)`GcQ9mwK9@E#rQAOB-ovfMOeuI0YmRcnRJ&rAfNnca}t z6O%xx*8rInE$KPqB9K|~_F7uhEM!JiYs%XtLuSxm;^tBXWIDI5I*4J$6*KFGUR6KH zG;J1b$MT0vwWs_8O~a6ReYtYsxfo=kJ>p)_xIiY~doG=sUdZJ5S!GKl0GTY6tB?5# zAd|`B9hLDr$h`Dn=w`nfWS$myIWyD@8NU{WTSu}WW1M|uQyg8?Z64f$O9-xmyW8Lr z2rj|h-QC^YWpHqvD~cb#6f_g>tW8^y`^Vfr3H zmZD_F%)&V7Z*XD0hNbTQMiag+z&*Av_~D7Y%ypci)QQuSvuD}~k=CSp>2v>O0$S1G zqKLYNq-q<6oncM|sA-Yf5NgGMTYamz!p5He?kBGAF_apiO3GWXryu8r6AH1DL zwV(vO!coH^CmfYtU?%Ah_5#HY(1Nu#$VDaHy9dortTXX6x@N&iy05KfLr@0lWa21k zd5?8+Af-JA4{4e&5}LMaU}=^W^gL_T@i6&IO8OyJVLf0sO^EB>iz>S?1WvVNB5 zjMCGbu+01LvXLz!ZO!0Nh2;ffnHYC_q&Xw174rVPeas?#3zTvybwXK)u^u%adnKJd zrp9~Z!YIC}vlnz53|-p3-}#-4nV{wF{kGZ7U%dTC>yU;ZT(y>-ze+?V0bE4dTA+ii zJkjenuD~0BeW6)GK02Au-Q>8&#s-swoPFIattZB2NQNVJ%`d%<+m*fogo>2)(l>X? zAjOzfCOBsn8+<7v#3^e_2ymtS8*}pC@Z+iuX~s$5Q^x&8|2RiP8Q{w}ELDUx&ihP9 zq%fggkJ)68K*S$!hFmKqxP}Q?EB;MF^km`1TY4QtBk`w`x-G(hW9mm_%RZKBRkNk} zH<=+CTlAO6+^AGLE;_E=bjs}{(jX<( zBI?}+Y0OGlYlkwF{w zh%(7LH+owB@7AMiXhRVYub%$z;n+XhiLO6{SF81$R342kx^J%i+2y|Z=s0Yyi+|m%8?avwKZn{<{ z)5g294K&DeJ!=3yX6MHvmDKPR_tWiOeD4Z44Dun1tNeH@nIeNyYG>H{0L;(;WJ`dc#8jW*P4hMgi zZ5x7;@$&tba@b6jRQ;g2)VlO-T)Jq^@`hM+Ayo62Dp@wsh5@JOdngPgRf>TJ#~T}_ z2Gu%K?@8!tVu$E> zsX%*Q{AYycuyN8QCS=t1VdZ7rgVI34MP3Cm!F02sH zjpbAlCpt@Y$LBO+VB1jmh$wmiX@2V_Ln?UVk*UAifazGPyOg2m2)5*qbYFA5z{sT| zjF;q~Frr6yPE=7u>CpFBpL;Ff!ED{3kO^*y5qGGemN_>ikJK0qApwk#As&(sLJucN z$DkEB|MkThe$>eff0=+6f4i0Lx7-yJuPH6(Hvbk&lb$2bY3v_F^C46-^W#q-w#f*# z7PYS!QtOnq&an{bb84G>$-9FX_{?sd{REv1)(y}KJK~xs%KMhPh~XtYgJv$NZtO`} zGR2Ct~%s3B`30~JsD<~;JK@)TtvT%X3-+Q zdx(u)uHMNro0HyQ=M_ZOGGwzrMdO%XP}Iel$u+|R_U)ovs1$sxr!zUEUVwIN`!l@k z=5-(<+0AKTuw%@Z<{!U6{(<}ROf0#sUW488RF3 zlG)(l`6?B0*xU+~p+GCtbQbJF0^rPJ%Y0v+c?sEk=OV;@@Ib*}-H~B6=?uxq;*L5M zDvPu5pN`$ZvXDUWw(h2l)kyHZ;+&T66$nI-GH%IfB@#bYs1#zxmJ@wigR;rab0UW$ zcu3xVXowhj&}Eb06&&7u!HV>#;}Z8QXAFKykQVQ9E%`;n|jZwsu zl$EPpopXc?9k-SA48c6&R2@%WMtH=Bki7|cd$@z8&(nhqFFHmkN6!G87aRs#7sVru zH+DyDl0FfxJ6y-)-QLo4!DGZ_`SqG*6Jp4b*Vi})J$lim=q4{&9QgS$il4(~zl57* zcCP9y>By<|W_^Ujy~*LS)W+;ns0**`-H%Y#5$63Qs#J|-0{WTk_OQO2!4&wQOQ88h z7_K4w0rYW*0}beuZoQ$L+}n$DocdAbAp`3Y*SlPW5@ofv_DaDt+hfGXe{!TQDV6u; zIWI{!8N+{Mb=INF*lrupm`Sw`}v1viD*9V zCmZHzjkpUHjcR!a(L8}#X=T=rNCYJfU5ADd-4F(8BhHC7l#Z9%+lbeGfcSuE*S6*! z+-}O6#xPq&V0$aKeovoiTzmUov`6SX>3JQd*1j|F^Towtl1l}m%X9>2LvSn+PLsea z|74`@4g)0=jnA=z%wPDkXr^I!+x!c7j)KUbQqyj05nUo)4pz=7in5MzTdcDdx9bv+ z-@qhL!^jw6j7>DYZ!F@`c@oiOdbsnxIyE4-COnZ9WhGsr2#J)KH8?qr2aKQM)A z?F{H6#mjwncynDz5K;GEz*Qb_HwsTwtSajPt0YOEm8H^17wq`~zIlQ1ZKT?`HDfHX#V{9YvtJETk{p5f~^QIqkq%c?0g z^o~J#k7pAqaUAn4f%Ib(RvRq)$*LjJ-;b5+)4$Xr{`k5}q%7g>qIkeCp6QY@;dWh{ zynPuqHb-rbs}%|Q?d+cxQkCFP$1L@b@Ex-!wv9<{RK%#SdIlO6Pnw*^+lWCFDk&t{ zQIMoYlr3;#sz8UB5_!4fE#fG|WZ2t-4gjT$z>LaO-yu99FrD6>ewxs@{O6Bj9I|qF z*Q$F1ESaq-X{!4)^nzgNkI@?}QfO_evCw)9wqBjN^@Nu2DR-6h>0JpjS>-@8#YH6QIPY8#nSUjCqQ%dhhP2MuziR0VtcxFN`dCsvI{blfW*U5`PnYi@7lq=tF?HNR zP3TaiJ`_H9q=GRore{SHxTqYP;z4hYe2Yqy$#MuZ_TCLWC>A?lH>YHH5u;0S)0XNd zh#j8lH$n@85>blyUo=K}ERpebztsdS!W28I-U)-q5H?teTrFri`>8;m8IkZ;8{i5H9xgCHFQ-j4Y!hTHet2*LgZz95S^?@wgYmePj5&$&6rAOGDeB&W6JI#f z*${Yz4n)2N)q_2%#ipTEV)b}+%1`<_y8ntIa{@d4720A%8CL?l%C}ILkgk60*Y|{L z4UO^812oGy6KgT7&I7eJf{gNZY89B%-$&lUij1}<}?%#hN!3Cw$zVmHW~%Pf)P(rq=Vs_L+n;SLC&!^)xisS3YN5f z54aIyLNA%so{Qw|L+3n_e@RPGalyR3tiw?N%YZ zFsK>sW6VF)GzRyyk0KVNz@34sl|i(zM=y@`zt{hPYrMxYoW7Cj3pvotWFmv{Bdh!m zjV^=5V#{>ZUK84w%8aK};5~Ox2}Cb+$gUWQAtB28~c|}UUiTea6Yd# z3r(Um7DwEgOnx48o-IL?Zjnu<9mh?OUXR*XsgvfH(khmw!T$Vv?mv?XEN?r8oe&0j zUEf6(+d!vtVxI$H=GrfO@+EP3c*dR^c4sBj6gf# zWX`Um?GcuZQHzJEl18lbpC7ZzKj^U~V9TNJU-(S-3U!z5b*%e6b=aNf(ROK5`LY+A7$ z6?FGjdLVFwA^p+dez~$FAq^_}iLnNCpPXD znL)BZa>1m-w-kCe0kHC&?9Fgu?)|=}WH83{&=35+P()M1%V|Znbw`%_JOfos5IPNM zd`CR-asm1NbSBl+jET8a;wqNi=yQCHjnwaT9CuhmT0o=%r`Q>sM%T59|M4MSIMj!0 z2>UCto-%7oH3}o`{U$M6QB1hN*T3^`Y!K%KpSX`;ge{qp5agQJYKqrQ~kN7!op_YC-E#OU^)>5 zNebFIhzE7=1ikM82u4glZjx(E7vG=8yZD3rAtq7KQ9wVo9)~Fs;Pd563(&7I5BfA7 zv?-l^H>xnhY4o-S@V6;=+5^P9SH9gOb3;JyH_7$~rXTmoleM3>=7^&oFXIW|6$Ll1 z&&)5T-iH;1wEfpT0A8-+`xyYV{nmA`{kcQZdVm*g&vZLKFRQ-eol)G zpd7Drz!~A|IGO)vVDo@6koS63;qr5W9pcyV<jj z-(U9?o7Wo4q>$uk%yTVCJx~FcEo?hx=i$jOazT95fdLe-u2-qp{~b*e&FEFZ@=!Gp zAV1s%m@5r?`Wt+`7Fz%C)4ll53g#r^FPT+$3D> z7fG!+g&=Jz?`eE{3Y9rGZePfHCcJEG{CYp2$@6G)AnFKWN*L{`{Wsz=VM?vjI%|uQod1Q+MX5Do6dog5K4Dmem0rOwpN7D)HYpOL}HN zZ!)|j+3t=27{1>n1MCO9hGqe!%5D`K!r>M^OKYHM7&I`npS3F?R(GXO%NycBWv?}! zp0OGkfTA6c1*g#nS2|`vZ(RT7I5lqTnx6Id*d5CgK4E+0JZRxsm<1>dcr{X1Vm%Ei z8dh47k+SG%V2od5!lCRR*B8qa9WDzXsAaRT1N>nVwItuX1-_8SD(C@=c;E{? zSo;peGlWw}+v4Z>VTZG$#k%tZU~l9Nlk2S61IR&)|7^Tr%fnUIPp4!8629-8v{gav zYx)#7uepBd5WJjH10*2)Ebx8DLZ%F>2X4sI^tU)`!5nn`aUiNvyOay(&@sn<1#?(^ ztDW5~#riuvbU!_Mc8gdh0yZUE)LEfTfVi2z*K-mSU#6_UMymxavG=K_3;8DC##y9J z^1__($Eh1J--@iepfV?t!WEMw*`{yM1unrK03?9-!bypcb&3+P1ePB<=ckBYXeuT8 zTVRC!|J)T9eh--7dbF6(7B-8i0ICXp*>Ca~U)8V@hRT!R%l!&<^XJ|ceE9`E-nOp2 z=mle#R))aO|2H(o&Tt(>wTd)9|6=Ug0w0E_kv>$P^}ii@>K=mkrT!K0X-~QX^R?hM z^OKz_n$IQyA3b+wVESQdXu&et-{<$10>mPXj@F*@LaoSvdzq1KfIBwJhmU`@!(Sv8 zUgw-C3blYLh{B}ztv-~mwr!rY8{8a_Mh!H9553M1>_(7!>tH%sVo(?Jt#d`2kZ7)a2{rlIyfCmiIJFO*UgPwvLY}iv(MXPi4R2%h^B&B~$rH#P z*=BDTOE{CZiMPueO1QAOVtJ`Aj9OTM(Sk{C@}b!y+^EUkIGJ`njCm3V!9G_HDU~yM zyZaK~JgjNb!hiV+7<4d^Air7*V+QDUj8>WO+D!p_< z5c<5QVN}#0Yw^F0p;V-4CFno_x0 z3=tw?jtSflJ1$dsW2WzfuDDcACMJv__3;6HZ z#Fm6nB+gfA+?QNYDDCL9F_fGqI5lVml?Z2W6}iDL1VO(xkJh@ZP9O+ zq6RrBY=*n-YEvIY0Z!A|kkFYb?_<3<_MPZyJBsd+8KxFQlL#AmO_#wQg$ShC7#7R< zui`2-e(Qu%(5ya6KgyiFF@GZ#K(*k5JuwrkTA^ z6^BzhY13PX&9IOF?UPk{Z#sUJepHoNYqFv-;VS<}nInJauQAw1O;#WpsxTfsqy< zmU;`{XBbvT>5M62V=0uBN9mOG9t_|39dL6Ajrgl3Tn8H#E0lB1*RGv?beMSF_NXBO zRt6;O%`~6z<5!`Mx1SmRb!^mm1=?>x&Fgf&2}S*dE0pM=+smqsrI|p`@3A-u*^HER zsy(=nXa{+27W3_eN?YwR98|54R|lPxUyNX=)AJF86$ARItM$JALxwY-EKdLq7gj6vGp#v-PBY24D-MY3Ip9 zQ7shBVLNjPA1Xs%d9c~N77!dBxyU}Ka8nEEF~0x)r)+^`@uuJ&A#Z&n|z9{rhl+7)$~8kRM4mXby_&bO>vBiGioME zZoV(-PrHcsPX%E}{Zh4k%J@_N+EtoPh zlv~|&5{u!KHt$V(C4KxhJhAn-oQSM0e&i%R7xC~nyxGl*gl;}t@r|-hoFepfpVLPPeXrM`hYI8AKKfOb((Zm z#o%;4cK-^k$8_%fcV)J7W~EU_k)pnR-tD~m)kia)ZL|x=XW7x-r{_IzF^|7z*aNep z{>Q<1-+=^b^gCE%Q{T=Va-a~|r}LU!Hs%iL(}8FM6W8wha-O{*A9ow0EALKs_JZPJ zU*Dbn+WU&XNWV8e+V_ZWf2OHw_x7FBJ#`g*Jm~tC*`QR%60!CX+Wz$D2FEnSHUwAj zBp2!mU(#PE9)x6|w_MeV9o7^=&Dc51WV8z!~qp2G03 zWA;EWKgX(JT8)_7{uk{;k5+I0mht$(rfJok(qtldV~e3#JHZc-$d zjKzvZWC<@=3PI)jvmYg7B_x!gPB?**NxS;hlj8d_6%XqyzV~Dp%K>Si-r%G2zi~Xk z>^#QiPjVT+##3!LOlcQLyurT<-|vKDaba~U?or=U$O7m%kEtF}amX;{B7Sl8YtHsc z%W+YlOPyOTK%X`)w}2h2OFa|wKA4>H7Gj}nx)|&>x}4K{IqgaDlI~pseINLGvVY`% zT5=tW8i&K$`~4Cw&aXtiYEAecH?b>D-_P;4`|(DrYDDG{PnB&3a`(wZ_uOp(`E|@@ zR|kNi;lDZ@eao$=N3hb}RWy02`Y^@#P?ewzM7K`jxMdq^K-=({L0$b3iMz#qkV^Yx zMXREWQ`HkTs{yzog{bISBTz$O?JImOm>=AyZp=DsrJ-A|b~5Vj7{`b#K({V!bb1mr zh4*gA&tG+4Sj+lk`P3!4=bP9ABo?aLzOy6WV&)y@9VqqNy691x=ylm>{DQ?P`|fd? zmvI<4`0@uXNI$`|Xu#ufX|P)7)SzZ^WcsahgMD%FRM`{lsC9eKA&YjrcC1U~b@N1NH`tb22B&sfl6zj-Ijvxu#`Y0XXx8RrH^_J&?#h<#jNS#PK;{lqso z7i4TSb@gac{bhxXb$o+sSG27Fb6ybH4GskJ890`TKhrI^MUpJ}0tXuFd>W@}Z`qN? zVjIg;Ien%n9~l`M;w^*-_(g^Zp2)YTapL3^ZR602Ngee7h|Obg1j|sv^trj=x{&;> zZiznDZGX~k9vpf8l5tF2I>={ud*F}FaDVD~*Om{wYUdnhVi#m!^YAiQ50TiEtvUb{D@D(ji^vXGDqn7ciq2mGDwr-7GwdGWE$>`y%6_{bTTB@w?d zfxqZ@sZHBLSXM4TloiH~gjMS^7f+n&a00ITQ}uhYyx(5}&$$ys@VIew`~H=U#8=M6 z3kYI)eIT7rVcdTiNR9@MrfswxOzCKtfUSZ{Jw)jud`D-nXf*zK94q|M0L^CSCo&@b zB`vkWpKF!UwbYI6!qgM_sB}REaqh3#irt2?W;|W$3zhFGxF7`J`d!*ntha9uQc_O9hTR<_RbvgcKWKLP&>=M<07P8mM_!||MO{mOe&701|=Bw z2<`C3QaO+}A;x-Z@>*{Oyu-xS5}I9lo2|P`Nm*wU*oVGxF5f-QTAl8`Bm27RJgUpv z-=psnI3HIrLpTUQ3iS=k^T4Zg-_purXHkCt9hR`x{Bth*Fev?x1iT|?+1|%=i*w1w zrV_V$zpt;Jw;m~MY1cXels2NadlkINz%Dw`dMSW_jg^?bnlz9=-{O-VT$(7)h`mDvl zKjSiZfj8<#^jhhT&C&F*!nPIn3}0n4jU2uH`f7WcMz59Kq*@Dw*$yYNrTgmgVmdGI z7crIYv=$fnl$Hk(odaIZmrXZQ2%mE`E zS1Z4t)`3i^KpNI0MRr+Oq&q(MPDQV{LI@nFfLnpUz|nz#R}4@U@&glDON#ZfI?EYB zwuNEY%9aW)Vx+Q1?@zn#Tu7C$bR2HA75`{g$qwf21UWf@Gd2ZkFJ?Hs+g8)jeCHjIKom&AH|q6h^J zWOB&k%FJNoi+Xer1?W{?mXX|)#Qknwp3$01tJl&eJ+{$$U(s2>uW>`u;Rgd3Rf56q zL@;f0T3lxj#QJH=Rw^VZ+`jvM1!j`z+iBV!Np2wTqKfn4(eJG#>u&^AEs}(Evsgkh zw;8w|qc^=wYz~I54Kj^VQBiAm1lCLBDAQ(jg30ruLy{KiwO*5Z);<7xR-={*1V3-Z ze1G?qnz|pj4PI_7;4Vg)x+lBa$<~Ztr(cOm4(qOvmxNa{m%2(K(u77ichji`w$A9? z@j=^`PJU7{H>Sh>k&ziY`Rl)G#vnxnPL{6^U6fVU5v;OD0?rxQnpT2&-5x!AyftT)#@epNN%(qp+KLzM0`ORq^uNy+<ma5H*KGth4QyEmfU!>Ua|EU!fhW1ABj^9EPbU+nxnD4f3bCX2=H2 zKhw1x`VRdO~^yd+wSfJ;yH$a9~VT|J1vWBbhfevT0z5#yeT@TiM(dNQAH2AK~pCzj*W| z{1XmT{rW(Ed1D_7xJN740S3qXFAVt~IQ73U=zn1K|H82Uf#?4V1MlNWm;HV6hP}J< zG?9yegbc3 zpN2j9pBZ0wUt241n-CtYXEf^cn-|$_YxorG4LZW{j1MoC(-&&)X6=}>`TShNvT4&Z zNyMMt23C0mn?WI!}2vDg{@ zPfn+gCPbBusUrqTX#=_s^DTEck&No4wmX`3_$y6uN1mcpmgmy&`-3ad(y`2!+68X( zxvEX47XU%VW{GWoixK>VDEcuj2fV~g5HT8qhC%O31kE|iS9uP{*z|BFNL>=M+$hES zf;pPO%hb@meX=Wy8vUwIb`>%;1+>rYt6Tj73T zi8B}xODOEY}vbAb93#{oDw9ob(y zK6G#wTpgOB^pW{vZ{5ESgZ}zk&VI}uA~%^#rjZ1}8iWa2iVrfvmfDIOa3@pC{R9Cj zC94hS>VF9wnthYg`ktQtueiTUy}BdIUFABVBb@5Pn!+SResl>eDI%=h*p8@&oAIP3 z&@ms3wjOIfdDIJpm1uSR-87Wp!;IJUzR&UEfBGKEE3wn`Hso%QW_zH6SqGSk!23~P zN<7#+c`{wV!j_Tc7;mqP7D(&9$Q>;S~ubnLlD*X`}_t^a+OQJ`4k#i~C59)HUE z{>CK+UGctBMP3F&E@-?~8GI4^a2pb8#IK@Q9Baj8oDF4`6vEx8#P`v38~XfsMm@Vf z&th)z)S7ZdKZW#?k$1oY&wh-dKlEai(b~Gh!r=q2LPFc{LtkxRY3m2y2mcd`%&}|P z&^v3lEURzhgAvy0!&7MNGLx9sQP}hhtuI2v<%7S{5gK*Dl{AQ%EznKu-`C21R7FZH zrTZbzuX5q`uRq{y|B55NYOo{y{^+_5wdduMb&~@a`Vm$=p1|5>ye=5`= zRiZ2|*k`l+Bq(Lmj2K@RFgO*pc4oEYk{foEN7ro)cju}Rnq96hxGb9blb8OSz~mp= zRO5yLB?sLUHZCwgc`k1Ko+>cs*TB}c*n);ckQ%?;r3JmX+1lKb1eFsp^p)>5k}$LT zU3EP8X8GgK^WefQJfme9{$Fj)0Ex9ElqvsR;RA-Y{v4Cg*IW1+mSD@CVpj?l!APBW9E?Ya>7#B*L`!l zg8}T}+3{-2e#+8xy9k1i@un5|9Y9Pt>x6h2HL|`(bZz7|lU4v9KND7kp@wwV)ta+; z6*5GF{iXwnZh$fi-J&+{%f$vnt~WGDv7& z|7FNZABDbKf(@%anw|&ae_lul>J8L`>D0!bB!n0&pP2nCNkc;!kzIZGl`{z5u5-WL ze>3wN1tUr=4cyINytHDMs@8V{b%&0T`n6#ebP^Zp+iZ^RMhMe)55^4vp#`O`g_$WZQ5Fm zMEy({WFjDv|0dU=9mbD~jXm!v->cK_LJR0Mmx8XxEtw;vcTb3XI%rI9u-H$f9jW8q zWRt$Q0e{%#9Y{sW44jT;4{fywBW^lm#yVv%cqOlhBc52zzkrZeb7+|#gIkdv_*Ts?Dab@Ng` z?xN&q5vM}Dm@19JhBzq}t#;ooV!$iY*_Fb&;s$<%645@-P8uXRJ^VzC#jsfMX}qz2 zM?s=|AL6?L1U6KyNs^eu*Y9$Rwgr=rVTnKbjsM0cCAxPX^A)p5w0c_`wYpsV*zg2@ zJO`10%dbHp4E30Gr@iTXh_X^!^OJnP>KkTqJg-oi>L^Uo(bx<-Lx>AF(1!2YbgKST z{$_F9%mY}wqtvtLVcb6Tb_0A0nAcYK?Sg0ZWr+OHX9g2SAOYi}XV?G!P0k|OPFyR$ z7yaI~_ARgE_v|5G!Ac@YSd{}+eCc}%Z4=nfG}x6$$)e8rZ@q(EbVmQ>pP4mi2+%EN zu(TWnA-9D4rO~51E%YJ!o1i4a9oSBKVTH924S|+6eRzRg4$S@#%iD?0OiD(RG#oj1 zHc1u4U=VcP-xTp@Gjpp1cV0+r zHCVU<^j|&zFCGsrd|YZ`1!eBk zAqeO)9p7w+pQ;_C9gtbHadp9*QORPpf#C<0d^ej;|NY0`vmpggy**NYdL@vXIwlvQ zC#35YA`Y*AFfCoD4SWTub2P1?mU1}B@=rrg(j;uxm94vTKXzCA>s;e|% zo`e&be;M6ZK<*>$%{`e(XTZI_?00JQHbwbO3x&aL`Y2+pl$LJv_s%n2lb~M z`~wD=nBvh(Gkrh;9@2HIjVNCWspLPd3JDVik0q_*cAw@iMKqoz(9!?RZQ2|3dX0YF z*DU_b&yfqz3oFrjuP0giBMZtc1B=4V`0<6_@?ch!Q4+plOTgYPlTG5KSsCCux58>o zvff^1uz#L=lBJpyEU7XQ$QTBswVuA;H96bu!^V8zOEw!+1IJ8zzBzB1NuPR7%nn_< zD=}r+VxvSVdiSN=1_yd=h!6Af{DkHk`t$#SL}V{$DHgzXJc=SRsoudWx$F&Js*-eV9=R}<$8Y!+xSP` zwcwkN!Q2~Pbq?Ej=eDD%-~5H~LSvMJM0x*cCb{0SiL+RWL-IQ57wnXe+AC6wOuBEB%A@0M5_{+zB$0OQh+ZuvDay;BnXx5hz!I5h9#2h#J<&75fA-8 zyS?J?PgcECti~y*Mc;K`;K>46Li*UeRW7km$~ShVD$&9eM&DHZ>F@3`Bt?MMzih^+ zppTI~FxITL;d;K^v+JLr2ue+F-Xi-L>3)_)ay>?(?rgYRTx`ej@U-`y!6 z>MUphwfVc_zT9Tf6ef# z#9v;s-bG2Vuj-pj!xPXTkw*` zbCyB#gcUt1`!wHbTeoDK1p4I)O_itk>i!m<-t$t{XYM%25nF-!^Nbiw|>;F}|+pQIqu zg029R@SS^HpPC;_K@$-nL{SbHx{(1I+nuXby zr1Tj0K5mHx%(lB+%%;&!7q;XzEF4AN4tvj|XXCC7AP|Bo0qp%YR6i7ItzeV8!YQJk z6cKE1eA8jpFn<5gFzpZfV6?xc45~Dh<$o&jPg`|N?Qf8rafS-oF`M5r!+E`kVymLN zF|Bp)Gm9+Z{<^`2BrXj~YX|4G zI@xOcQ==oybE}$sb;YMbUwRAy1V2aGe#P%}F~F?lmYjej*P(i8PNMP~?&keoFU+X# zb+{%SX2N_LtSd|>QW5HwE!(S{A;LlCkxeLQ11C{$$=U$E0iGf8xSh4HOIdcQyD}O5 zv2|9DQnu~|A9^#U5EPpybdBdRovHwrVgJPy{U)f+W+fFUOR~IL#Sl9%NK3Ht~IPdS*9Xuh|Xf=%lrHvU7Xm&!>HY^ zfZFg%3QlV;G&;LqA*WV0!kNnd&xDf(@jV7#iFNPr=1HaAF+TnHfOagqX^^Jm$5^XQ z76)}o;2_;AJ6bXx8MnWeBmEqnk@l=Vr@B9Eps}`do_%*e#q3_*2%1=Z|0!KT=Oez0 zGl4P%@}0sc}OadpBU5R}duSaGp#l)U z)AyWZvHbGyagnz$ML|Se)9$ETl?BEtiN8{UFX%{`kt}u`ZWWIcNTJFKZ+zd2LT~0h zN0?sxO~4hebMTs}M581rlQpovj6=BqG_Jmr}gOSB}c-rH*1o#|#vWDR4{abYSVQy6eT@&XJ9_o|K z4jr8c&rbsdu2W=IwJjI6=MnVJzv4oPdaEDdvOxz&Yh=1n2lGbQ^R9 zV*W<7z*5gHn6%k`5P#gKhb{q3VM#lW1UDHj!T5Ok9L!GH2$1=$9y;7W>h*ch-r)Pu zlqT^mk`|qQ`Hx?IqFs8giB7?Xha3O5hQzgrEBi;tAE+hbSD^#eqC4#oamC>PzqT=_ zFP>IBk|rcFq*x!wu~ztBi>Fj#8+QtVhoz`8@{*x;sr30)k#R)ivLp+nUNOZBjq9(k zyhSord=FFe{ScIx7J6Ebpawz=fvBr<=*@sp=kYi6pfFg3^6OmEdQqB3icVMv-im+% zVNCGOpRt;=TjRg?rZxYKcfnr%dHpkaoBcP5CM@`IhQMbQ<<5$>8*bKa9>vukdzJsD zQKq9YpkqQMH+QHje;ji2uoP=)F%UtBJM#9ibM9{Q`fNV*pX=+-Mbr<+ZI=~4*Ulvo z5C_w13Hfq>yt%Nr@edF=C zcFS;q#AwIyDtdK0cJg9_;TyL14*xoLUvnPDV)Ib%jXmfDH9q#?YrbBFQl9UYF07>M zqt0t5&a8fR$NveU0A2qi*w>xu+nIgrhkb2O_ddSEFYK#Utk4P*d0=0@u<)vW;XV73 z?KL;^i|Xx*gv^pNqo3I4TlLOrzou+o;hk*}Yd7CMKWp+QV3BH{EcVaSZPt5OyB~5M z9@BpY$O*~J&~x!vJW^8Or#J^>Wa@dq4vFS=EsPJ(Oex9GbJT-qW`SwOdb~171v5wm Y(kabS{xoXM)02aw8$^ZZW diff --git a/src/sed/magdziarz/.Rhistory b/src/sed/magdziarz/.Rhistory deleted file mode 100644 index cd1272f..0000000 --- a/src/sed/magdziarz/.Rhistory +++ /dev/null @@ -1,146 +0,0 @@ -data=read.csv("components.csv",header=T) -plot(data,log="xy",ylim=c(.0003,.1),xlim=c(.001,1000)) -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -data$test_sum = data$blackbody_sum_1 + data$comptonization + data$excess -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -data$test_error = data$total - data$test_sum -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -warnings() -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -head(data) -mean(data$total - data$test_error) -mean((data$total - data$test_error)^2) -mean(sqrt((data$total - data$test_error)^2)) -data$test_sum = data$blackbody_sum_2 + data$comptonization + data$excess -mean(sqrt((data$total - data$test_error)^2)) -mean(sqrt((data$total - data$test_sum)^2)) -data$test_sum = data$blackbody_sum_1 + data$comptonization + data$excess -mean(sqrt((data$total - data$test_sum)^2)) -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -data$test_error = data$total - data$test_sum -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -data$test_sum = data$blackbody_sum_2 + data$comptonization + data$excess -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -data$test_error = data$total - data$test_sum -mean(sqrt((data$total - data$test_sum)^2)) -rms -?rms -library(rms) -data$test_error = data$total - data$test_sum -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -data$test_error_sq = (data$total - data$test_sum)^2 -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -data = read.csv("components.csv",header=t) -data = read.csv("components.csv",header=T) -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -data$test_sum_1 = data$blackbody_sum_1 + data$comptonization + data$excess -data$err_1_sq = (data$total - data$test_sum)^2 -data$test_sum_2 = data$blackbody_sum_2 + data$comptonization + data$excess -data$err_2_sq = (data$total - data$test_sum_2)^2 -data$err_1_sq = (data$total - data$test_sum_1)^2 -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -warnings() -data$err_1_sq = (data$total - data$test_sum)^2 -data$err_1_sq = (data$total - data$test_sum_1)^2 -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -attach(data) -plot(err_1_sq ~ energy..keV.) -plot(err_2_sq ~ energy..keV.) -head(err_1_sq) -head(err_2_sq) -plot(err_2_sq ~ energy..keV.,log="x") -plot(err_1_sq ~ energy..keV.,log="x") -plot(total ~ energy..keV.,log="x") -contour() -?read.table -TOTL1216 = read.table("../../../TOTL1216.tab",header=T) -head(TOTL1216) -contour(TOTL1216) -man contour() -?contour -contour(TOTL1216$Hden,TOTL1216$Phi.H.,TOTL1216$Eq_Width.A.) -?contour -contour(TOTL1216$Phi.H.,TOTL1216$Hden,TOTL1216$Eq_Width.A.) -TOTL1216 -contour(TOTL1216$Phi.H.,TOTL1216$Hden,TOTL1216$Eq_Width.A.) -contour(y = TOTL1216$Phi.H.,x = TOTL1216$Hden,z = TOTL1216$Eq_Width.A.) -?contour -?contour -contour(seq(1:4),seq(5:8),seq(3:6)) -?contour -?op -?op -?opar -?contour -filled.contour -?filled.contour -filled.contour(seq(1:4),seq(5:8),seq(3:6)) -contour(y = TOTL1216$Phi.H.,x = TOTL1216$Hden,z = TOTL1216$Eq_Width.A.) -TOTL1216$Hden -TOTL1216$Phi.H. -?sort -TOTL1216[order(TOTL1216$Phi.H.,TOTL1216$Hden)] -TOTL1216[order(Phi.H.,Hden)] -TOTL1216[order(TOTL1216$Phi.H.,Hden)] -TOTL1216[order(TOTL1216$Phi.H.,Hden),] -TOTL1216[order(TOTL1216$Phi.H.,TOTL1216$Hden),] -TOTL1216[order(TOTL1216$Hden,TOTL1216$Phi.H.),] -library(ggplot) -library(ggplot2) -library(ggplot2) -ggplot(Statsrep, aes(x=X, y=Y, z=Zn)) + - geom_density2d() -ggplot(Statsrep, aes(x=TOTL1216$Hden, y=TOTL1216$Phi.H., z=TOTL1216$Eq_Width.A.)) + -ggplot(Statsrep, aes(x=TOTL1216$Hden, y=TOTL1216$Phi.H., z=TOTL1216$Eq_Width.A.)) + geom_density2d() -library(stats) -ggplot(Statsrep, aes(x=TOTL1216$Hden, y=TOTL1216$Phi.H., z=TOTL1216$Eq_Width.A.)) + geom_density2d() -?ggplot -ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + geom_density2d() -TOTL1216$Eq_Width.A. -ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + geom_density2d() -ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + geom_density2d() -ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) -?ggplot -?interp -ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + geom_contour() -ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + geom_contour(alpha=0.5) -ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + geom_contour(alpha=1) -ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + - geom_tile() + - coord_equal() + - geom_contour(color = "white", alpha = 0.5) + - s - theme_bw() -ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + - geom_tile() + - coord_equal() + - geom_contour(color = "white", alpha = 0.5) + - theme_bw() -ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + - geom_tile("TOTL 1216") + - coord_equal() + - geom_contour(color = "white", alpha = 0.5) + - theme_bw() -ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + - geom_tile() + - coord_equal() + - geom_contour(color = "white", alpha = 0.5) + - theme_bw() -?ggplot -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -plot(err_1_sq ~ energy..keV.) -plot(err_1_sq ~ energy..keV.,log="x") -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -plot(err_2_sq ~ energy..keV.,log="x") -plot(total ~ energy..keV.,log="x") -mean(sqrt((data$total - data$test_sum)^2)) -mean(sqrt((data$total - data$test_sum_1)/data$total)) -mean((data$total - data$test_sum_1)/data$total) -mean((data$total - data$test_sum_2)/data$total) -plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000)) -ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + - geom_tile() + - coord_equal() + - geom_contour(color = "white", alpha = 0.5) + - theme_bw() -q()