From 28f5b9faf7d8840b0d120705ad1ea092d1cba9db Mon Sep 17 00:00:00 2001 From: James Turner Date: Sat, 3 Jul 2010 16:30:48 -0400 Subject: [PATCH] JT: added PyFITS plugin, tests and data --- scikits/image/data/multi.fits | 988 ++++++++++++++++++++++ scikits/image/data/simple.fits | Bin 0 -> 138240 bytes scikits/image/io/_plugins/fits_plugin.ini | 4 + scikits/image/io/_plugins/fits_plugin.py | 149 ++++ scikits/image/io/tests/test_fits.py | 71 ++ 5 files changed, 1212 insertions(+) create mode 100644 scikits/image/data/multi.fits create mode 100644 scikits/image/data/simple.fits create mode 100644 scikits/image/io/_plugins/fits_plugin.ini create mode 100644 scikits/image/io/_plugins/fits_plugin.py create mode 100644 scikits/image/io/tests/test_fits.py diff --git a/scikits/image/data/multi.fits b/scikits/image/data/multi.fits new file mode 100644 index 00000000..daa20229 --- /dev/null +++ b/scikits/image/data/multi.fits @@ -0,0 +1,988 @@ +SIMPLE = T / Fits standard BITPIX = 16 / Bits per pixel NAXIS = 0 / Number of axes EXTEND = T / File may contain extensions ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator DATE = '2001-09-10' / Date FITS file was generated IRAF-TLM= '2010-07-02T19:33:03' / Time of last modification COMMENT FITS (Flexible Image Transport System) format defined in Astronomy andCOMMENT Astrophysics Supplement Series v44/p363, v44/p371, v73/p359, v73/p365.COMMENT Contact the NASA Science Office of Standards and Technology for the COMMENT FITS Definition document #100 and other FITS information. INSTRUME= 'GMOS-N ' / Instrument used to acquire data OBJECT = 'NGC1068 ' / Object Name OBSTYPE = 'SCI ' / Observation type GEMPRGID= 'GN-ENG20010910' / Gemini programme ID OBSID = 'GN-ENG20010910-15' / Observation ID / Data label DATALAB = 'GN-ENG20010910-15-001' / DHS data label OBSERVER= 'Jorgensen' / Observer OBSERVAT= 'Gemini-North' / Name of telescope (Gemini-North|Gemini-South) TELESCOP= 'Gemini-North' / Gemini-North PARALLAX= 0. / Parallax of Target RADVEL = 0. / Heliocentric Radial Velocity EPOCH = 2000. / Epoch for Target coordinates EQUINOX = 2000. / Equinox of coordinate system TRKEQUIN= 2000. / Tracking equinox SSA = 'Walther ' / SSA RA = 40.66966667 / Right Ascension DEC = -0.01316667 / Declination of Target ELEVATIO= 0. / Current Elevation AZIMUTH = -135.173831944444 / Current Azimuth CRPA = 0. / Current Cass Rotator Position Angle HA = '+01:14:25.87' / Telescope hour angle LT = '05:00:13.7' / Local time at start of observation TRKFRAME= 'FK5 ' / Tracking co-ordinate DECTRACK= 0. / Differential tracking rate Dec TRKEPOCH= 0. / Differential tracking reference epoch RATRACK = 0. / Differential tracking rate RA FRAME = 'FK5 ' / Target coordinate system PMDEC = 0. / Proper Motion in Declination PMRA = 0. / Proper Motion in RA WAVELENG= 4750. / Effective Target Wavelength RAWIQ = '80-percentile' / Raw Image Quality RAWCC = '50-percentile' / Raw Cloud Cover RAWWV = 'Any ' / Raw Water Vapour/Transparency RAWBG = '80-percentile' / Raw Background RAWPIREQ= 'NO ' / PI Requirements Met RAWGEMQA= 'USABLE ' / Gemini Quality Assessment CGUIDMOD= 'Basic ' / Carousel Mode UT = '15:00:14.2' / UT at observation start CONID = 'GMOS v8.0' / Detector controller ID DETECTOR= 'GMOS + Red1' / Detector name DEWAR = 'Cryostat A' / Dewar name DEW-TEMP= -121.220657276995 / Dewar temperature (Celsius) DETSIZE = '[1:6144,1:4644]' / Detector size NCCDS = 3 / Number of CCD chips NAMPS = 1 / Number of amplifiers SHUTTER = 'OPEN ' / Shutter state during observation AMPINTEG= 5000 / Amplifier integration time TIMESYS = 'UT1 ' / Time system used DATE-OBS= '2001-09-10' / UT Date of observation (YYYY-MM-DD) TIME-OBS= '15:01:11.5' / Time of observation GMOSTHT = 'GOOD ' / Shutter health NSUBEXP = 1 / Number of sub exposures UTSTART = '15:01:11.5' / UT at observation start UTEND = '15:02:11.5' / UT at observation end EXPTIME = 60.0029079914093 / Exposure time in seconds ELAPSED = 60.0029079914093 / Elapsed observation time in seconds DARKTIME= 96.9185329675674 / Dark current integration in seconds INPORT = 1 / Number of ISS port where GMOS was located GMOSCC = 'GMOS-N (V2-8)' / GMOS components controller s/w MASKID = -99 / Mask/IFU barcode MASKNAME= 'IFU ' / Mask name MASKTYP = 0 / Mask/IFU type (0=none/-1=IFU/1=mask) MASKLOC = 1 / Mask/IFU location (-1=unknown/0=FP/1=cassette) FILTER1 = 'open1-6 ' / Filter 1 name FILTID1 = 20000018 / Filter 1 barcode FILTER2 = 'g_G0301 ' / Filter 2 name FILTID2 = 20000018 / Filter 2 barcode GRATING = 'MIRROR ' / Grating name GRATID = 30000008 / Grating barcode GRWLEN = 0. / Grating central wavelength (nm) GRORDER = 0 / Grating order GRTILT = 65. / Grating tilt angle (degrees) DTAX = 103.2 / Detector translation X position (microns) DTAY = 56.64 / Detector translation Y position (microns) DTAZ = 2842. / Detector translation Z position (microns) DTAZST = 0. / Focus at observation start (microns) DTAZEN = 0. / Focus at observation end (microns) DTAZME = 0. / Mean focus during observation (microns) DTMODE = 'FOLLOW-XY' / Detector translation stage mode ADCMODE = ' ' / ADC mode GMOSDC = 'GMOS DC ' / GMOS detector controller s/w DETTYPE = 'SDSU II CCD' / Detector array type DETID = 'EEV9273-16-03EEV9273-20-04EEV9273-20-03' / Chip IDs EXPOSURE= 60. / Requested exposure time in seconds ADCUSED = 1 / ADC used? (1=yes/0=no) DETNROI = 1 / No. regions of interest ADCENPST= 0. / Start entrance prism angle of ADC ADCENPEN= 0. / End entrance prism angle of ADC ADCENPME= 0. / Mean entrance prism angle of ADC ADCEXPST= 0. / Start exit prism angle of ADC ADCEXPEN= 0. / End exit prism angle of ADC ADCEXPME= 0. / Mean exit prism angle of ADC ADCWLEN1= 0. / Lower wavelength for ADC calculation ADCWLEN2= 0. / Upper wavelength for ADC calculation DETRO1X = 1 / ROI 1 X start DETRO1XS= 3072 / ROI 1 X size DETRO1Y = 1 / ROI 1 Y start DETRO1YS= 2304 / ROI 1 Y size AIRMASS = 0. / Airmass averaged throughout the exposure AMSTART = 0. / Airmass at start of exposure AMEND = 0. / Airmass at end of exposure NSCIEXT = 3 / Number of science extensions PIXSCALE= 0.1454 / Pixel scale in Y in arcsec/pixel NEXTEND = 3 / Number of extensions OBSMODE = 'IMAGE ' / Observing mode (IMAGE|IFU|MOS|LONGSLIT) GEM-TLM = '2010-07-02T20:32:41' / UT Last modification with GEMINI GPREPARE= '2010-07-02T20:32:41' / UT Time stamp for GPREPARE END XTENSION= 'IMAGE ' / Image extension BITPIX = 16 / Bits per pixel NAXIS = 2 / Number of axes NAXIS1 = 132 / Axis length NAXIS2 = 288 / Axis length PCOUNT = 0 / No 'random' parameters GCOUNT = 1 / Only one group BSCALE = 1.000000E0 / REAL = TAPE*BSCALE + BZERO BZERO = 3.276800E4 / EXTNAME = 'SCI ' / Extension name EXTVER = 1 / Extension version INHERIT = F / Inherits global header ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator DATE = '2010-07-02T19:55:22' / Date FITS file was generated IRAF-TLM= '2010-07-02T19:59:26' / Time of last modification CTYPE1 = 'RA---TAN' / R.A. in tangent plane projection CRPIX1 = 196.375233233945 / Ref pix of axis 1 CRVAL1 = 306.871725219222 / RA at Ref pix in decimal degrees CTYPE2 = 'DEC--TAN' / DEC. in tangent plane projection CRPIX2 = 143.509119025456 / Ref pix of axis 2 CRVAL2 = 37.3733845984976 / DEC at Ref pix in decimal degrees CD1_1 = -2.9266657155837E-4 / WCS matrix element 1 1 CD1_2 = 1.36727521936194E-4 / WCS matrix element 1 2 CD2_1 = 1.36706238220991E-4 / WCS matrix element 2 1 CD2_2 = 2.93036096725254E-4 / WCS matrix element 2 2 RADECSYS= 'FK5 ' / R.A./DEC. coordinate system reference EQUINOX = 2000. / Equinox of coordinate system MJD-OBS = 52162.6261539214 / MJD of start of obseration CCDNAME = 'EEV 9273-16-03' / CCD chip name(s) AMPNAME = 'EEV 9273-16-03, right' / Amplifier name(s) GAIN = 2.1 / Amplifier gain RDNOISE = 3.5 / Readout noise CCDSIZE = '[1:2048,1:4644]' / CCD chip size CCDSUM = '2 2 ' / CCD sum CCDSEC = '[1:2048,1:4608]' / CCD section(s) DETSEC = '[1:2048,1:4608]' / Detector section(s) DATASEC = '[33:1056,1:2304]' / Data section(s) BIASSEC = '[1:32,1:2304]' / Bias section(s) FRMNAME = 'dhsFrame' / Name assigned to a frame of data. FRAMEID = '0 ' / Frame Id DATATYP = ' ' / Type of data frame WCSDIM = 2 LTV1 = 0.4375 LTV2 = 0.4375 LTM1_1 = 0.125 LTM2_2 = 0.125 WAT0_001= 'system=image' WAT1_001= 'wtype=tan axtype=ra' WAT2_001= 'wtype=tan axtype=dec' END <r7^Ņ߅ .8A; #],/kʅ<]yÆņȆˆΆ҆Ն؆چ߆ + !#%())+.0149<=<;=AA: ?b\Ʌ>^z†ĆƆɆ̆Ά҆Ն؆ۆކ ##!#%&)+,+.12249>?@@ADE?#_Ä*><]zĆȆˆ͆φцԆ؆ۆކ + !('&'()+./0256559?CDEEGIF +1~Nd:Zw†ƆˆΆІцӆֆچކ + $))(*+,.02479;;:?BDGIIKMH F!O)VwÆņƆʆΆ҆Ԇֆچ݆ #&()**-1578::;<=@BDEGHJMNH N8 mĆ Gr†ņdžɆ̆Іӆ׆ۆކ + "$')+++-158:;===>BDEFIKMPPJ T˄N')`Ćdžʆ͆Іӆֆۆ߆ + "$'*,.//1369;=@BBCEFGILPSSSL Zلd؅@@tÆdžˆ͆φӆֆچ߆  "$'+.0256899;=@DEFHIJKMPTVUSN ^vU QņˆΆφцՆ؆݆ !$'+.0137:<<<;=ADFGJLMORTVWUTN  bfdž]Ædž̆φІӆ؆܆ +  !!#&),.02469;<;;?BDFHKMPTVXYYXVO f sӆ%eÆƆɆ̆φ҆׆ۆ !$%&')*,.02468;<;=BCEHKLPTVXZ\\]ZP + + + i|چ,iņɆˆ͆Іֆچކ !#&()*++,.0258;=>?ADEGKNPTWWY\]_`]S + + + +!kކ/jÆdž̆ΆІՆچ܆߆ !#%(*+,,-/1259;>@ABDFGJMQTWXY[]_aa_V + + + + +!k/kÆņʆφцՆچ݆ކ  "$'*-.--/2469<>@@@CGHJMPTWXY[]^`ab`X + + + + !k 1l†ĆƆˆцԆ؆܆ކ !"#&)+./.0258;=>ABABEILNQTWXY[]^^`bcbZ + + + !k!3n†ņdžˆцֆنۆކ !&&&()*.0/257;?@ADEEGJNPRUWXXZ]___acee] + + + + + "k"4oĆɆ͆҆ֆنۆ߆ + $*,++*,.0048:=@BEHIILPSTUWYXX\_``acegh_ + + + + + "k"4oņʆφԆֆ؆ۆކ  &,..--./026:=@BDGJKMOSXYXYZYZ]`aabdgii_ + + + + #k"4qÆdž̆цԆ׆ن܆߆  %*.00000036:?DGIJLOPRV\][Z[[\_`bccehig_ + + + + #k#6t†ƆɆ͆цՆن܆ކ +!$).122102359AHLMMPTUVY]]\\\]^`acddfhhfa + #j#8uņɆˆΆ҆ֆچ݆߆ "%'*.12101346;CIPQPRVXZ\^^]]^_`bdeeeghhgd + #k#7sƆʆ̆φӆֆن܆߆ "!!"$'*,-.0212446:?DIPRRSWY\^_____`acfiigghijf + #k$5qƆɆ̆ІԆ׆ن܆߆ "#$&'*/10/145667:>AEIORTVWZ^`abba`abeillkklmmh + #k&7tņȆ̆цԆ׆چ݆ #%&(+0210269999=@CFKOTXZZ]`bdfedcceilmmmnpqok + #k';z†ĆȆ͆Іӆ׆چކ #%'),/00137:;<=?BEIMRX^aba`aeikjgfilmmlmpsrql + #j'?}ÆƆʆφ҆Ն؆ۆ߆  #&)+...03569<>@ACGKQX`enpg`_cjonjilmlllmqssrn + $j'@}Ædž̆цԆ׆چ݆!&*-///24579yĆɆΆ҆׆ۆކ #(,/13447;=?BEHKPUYZ\\]]\]_bfjmopomnqstuvyyr + %m*=yņʆΆ҆؆܆ #(,/1346:=?BEILOSWXY[\\[[]`bfilsxvqqstuxz||u + %m +=zņɆ͆ӆ׆ۆ !%)+.1579<>ADIMPRTVXZ[\\]]_acdfjs{ytqruw{~y + %m ,>|ņʆφԆ׆ۆ + $')*-26:<>@DILPSTUWY[\\]__`cccejqvvsrtx{~~ + &m .@~Ædž̆цԆ؆݆ #&()+/48;>@DHLORUVXYZ[\^_abccbchlpstrtw{} &m /AņʆφӆՆچ߆ !%()+.25:=@DHLNPSVXYYY[\_abddcbejnprstw{} &m 0Cdž̆цֆن݆  #&(*-147;?DILNPSTWYYYZ[\`ceggfgimopqsv{~ &m 0EȆ̆ц؆݆߆ !$''(),/36:=BINOORUVXYZ\]\]bfilmmnopqrsvz~ &m 0EɆΆӆنކ !'+--,,/369BDGILOQSVWXZ]acdegjlmoruvwz~ + 'l1KƆφԆ׆ۆ  ')(*.114;BFKNPQSVWXYZ]`dggikmnopsttv| + &k1KƆІԆ׆ۆ !""&+/26EGIJMSWWUTV[_`_`begilmloruwx{} + &j3OʆԆچ #15/,.26:AFHGIMSVTTTX]a`^_beghjlnqsuwz|~ + &j3N͆׆݆ =moG1159@DFHKMORTVY]_^`cfggimpt{ + + + + %j> _Ɔ؆ +#'*.0368:<>ACEFILORTUY\^_behjklorv~ + + + + %i@ aȆن $(,/147:;@ABCCEINPRUWZ]`bceiosvyz|~ + + + + + + + $i CfɆ܆ #(,/268:=ACCBBBBFLOOQUX[`cedcipvy|}} + + + + + + + + $h Egˆ݆  $(*.367:>CGHFCBDIMOPRUY^adghgkpuz} + + + + + + + $i Fgˆކ "%')-2469CGIJKMPU]a]Z]_aeilpvȈ=ÇȇƂ + + + + + + + %i …Kmφ +%*/24556:=@CEHJKMPTWZYZ^abdhjnv+ԏˋUVLJˇǂ + + + + + + %i ÅLm҆ +  '.465568<@CEFHKLNQTVY[\_bdehkow݉k:6ƇLJ‚ + + + + + + + + + + %i ÅLnԆ  $*398668ADFINSZ`]YWWXX[^afjmpqrtw|ÇȇӇ݇ۇ҇Ȃ + + + + + + + + + + + + + + + + $h ąN!vۆ !&*-257:?CEHJMPUXXVUUVY[_cgknqrtwz‡Ƈ̇ч҇ЇƂ + + + + + + + + + + + + + $i ŅO!wކ "%(,1358=BFJKKMQTVUUUWZ]aehknprvy|ÇLJʇχ҇Їł + + + + + + + + + + + + + $i ŅO!w߆  $&)-1237CGIIIKMQUWWXZ\`efhlnprv{ÇLJʇ̇ЇӇӇ˂ + + + + + + + + + + + + + %i ƅQ$y߆ !&)*+/37<@DGHHIKMQTVX[]^begkoqrux}ÇƇɇ̇χӇԇЂ + + + + + + + %i DžS†&zކ "&(*,/39=@CFGHJLNPRTX\^`cehlprtx{‡ćƇʇ·чՇׇӂ + + + %i ȅTĆ'{ "%'+.15;>@CEFHKNPQRTWY\_beilosvy}ćƇɇ·҇ԇևهԂ + + + %i ɅTĆ(| !"%+036;>ACEFGKOPQSTTVZ^bgklouyz}ÇŇLJ̇҇ՇՇׇهՂ + %h ȅTĆ(}  #&+036:>BCDEFJNPQSTTVZ^cimnqw{|~‡ćŇɇχՇ؇؇ڇ܇ق + %h DžTņ(~ !%)-147:>ACDFHJMOQRSUX\`djnqvz}ÇŇȇ̇҇؇ۇ݇߇܂ + %h DžUƆ)~ "',/269;>@BEIKLNPQRTX[_bfjoty|ćLJ̇чՇ؇݇ + &h ȅUņ)~ + $)-036:=?@CFILNORTTVZ]`cglqw{~ćɇχԇׇهއ &i ȅTņ)} +  &*-147:=@DFGJMNPRUWZ]_`chmsxz~ćɇЇՇ؇܇ &i DžTņ*~  &*.147:=BGGHLOOPRTW[]^`dhnswx|ÇƇʇЇԇه߇ &h DžTņ*  %+/136:>DGHJOPPQRSVYY\aeimrvx{Ňȇ̇ЇՇ܇ &h ƅTƆ*  %+/026;?DFHLOPQRSTVWY]bfimrvxz}Ňʇχׇ҇އ  %g +ƅTƆ+ !'+./38@CDGKNRWZ]`bejnqsuxz{}ćɇ͇Їԇهއ  %e…QÆ) "&,5:88:=@BDFJLOUZ\^acfjmprtuw{~‡Ƈʇ̇Їԇه߇   %eąSņ+ +"%)0579;>BDFJLNQW[^`cehiknrsrtz‡LJɇ͇҇Շڇ  %fȅXɆ.  "$).37:<@BDHLNORVZ_bdehhjnqstv{ÇLJɇ·Ӈ؇އ  &g ΅aц2† !#&*.159<@CEINPRTVY^bdfhhjmoswy|ŇȇˇЇՇۇ   &g Յkۆ:Ɔ  $(+-0368=BEINPRTUY]`cgihikotwy}ćȇ·ԇ؇݇   &g܅vFʆ + %(*-1367;AEIMOQSUX[_cfhijkpuwy}‡ȇчׇڇއ   &hㅁP͆ !$'*.258:>DGILNORUWZ_begiklquwy|‡ʇчևهއ   'i녌XІ  #&*/259=BGJKMNORUX\`cegjknruwy|~ÇʇχӇ؇އ +  'j󅘆`ӆ +"&,0149?DHKLMNPSWZ]`cfhikosuwy{ćLJˇ҇ڇ  (ki؆ "(-136;AFJLNOPRUZafedfgilosuy{}ÇŇƇˇՇ݇  (ks܆ #).26:>CGJMOPRTYglfgjmoruy|ćŇȇχ؇߇   (i%{ #).49<@DGJLMNPT]ӈ{hfilmqtw{‡Ňʇ҇ه߇   &d-  %*069<@CEIJJLOS]/ׇihjlnrux}Ƈ̇҇ه #\셲4Æ !&-47:=@ADGHILPT[rˇwiiklnsw{LJ͇҇؇ !Sۅ9Ȇ  &/7:<>@ACEGINRVZ`mvqjgijlpty}‡ɇχӇڇ  LڄŅ=͆ !(19<=?@ACEGKOTXZ[^`bcehknrw|~ćˇчׇ܇ !Gʄ@҆ + #)18;>??@DFHLPTXYZ]_`adhmqux{}}ҡ҇ه߇  #HÄCֆ !%*069<=?ADFILOSUWZ^abbehntwyz{}Ƈ͇ԇ܇  &L„F؆53!&,1478;>@BFJKMPSVZ^acdfjotwyz{~ȇЇׇ݇ )QDŽG؆/jd."'-3777:=?AEIIKNRW[]_aehlpsvxz|‡ʇЇׇއ  +V΄G؆ *\X-!$(.4877:=?AEGHKOSX[\^aeimpruxz}ćˇЇև݇ +  -ZքD׆')"$).1468;=@BDEIMPTXY[_bfjmoqtvz}LJׇ̇҇݇   /^܄Bֆ"%)-/269<>@BCFJMPTWY\adgjkmpruy~ćɇχӇׇ܇  0b…@Ն"&*-036:=?@ACGKNQUXZ^acfijlpruy}LJˇЇԇ؇݇ !! 0fх<ӆ  #&*.258<>@AADGKNQVZ\_acfhimpsvy}ŇȇˇЇՇڇއ && /k六9φ "$&+/36:>?@ACEHJMQVZ]_bdfgjmqtxz~Ňȇˇчׇ܇އ߇"+) .o6Ɇ  #$'+/37:=>?AEGIKLOUY\_bdehknqux{͇ЇƇ‡Ƈˇчׇ܇އ߇&,$ -q!0ņ "%(*-269<>?ACFIKMPUX[^cefhlorux|هۇɇćɇχև݇ %( +o") !%(),058<@@ABCFJNRVX[_bdehlnqtw|~ʇ̇Çćȇ·ԇۇ!"  )l"y !&(*,/48>CECABEINSWZ^accegjmqux{}‡Ňȇ̇Їׇ    'ks  &)*,/38@JLFBCGJNRV[bhheefhkoswz}ćƇɇ͇Շއ +  &jk݆  &)*-047>JMEBEHJLQV\dllfcdgjlpuz~ÇŇLJ̇ևއ + &jb؆ + $(+.1359ADAAFHILQV\affb`cfghmtyÇŇLJ͇Շ݇ &jYԆ + !&+.2457:=>@DGJMQVZ^__^^beehmsy~ÇƇʇЇև܇  'kꅅPφ !&*/3457:?BFIMPQRUX[_abcglsz}ƇˇЇԇه܇݇ + &jӅdՆ:Ɇ #&),/479;<=>ADGKMPSVXZ]`behlrx|ȇ·҇Շ؇ڇ܇  + &i҅bԆ8Ȇ "%(*.2479:<=?BEHKOSVWY[^cfimsy}ˇՇۇ܇ۇڇ݇ + + &i҅bԆ7Ɇ "&(*-03568:=?AEHLORTWYZ]afkpuz}ɇׇ߇܇  + + + + &j҅bՆ9ˆ $')+-/2446:>?BEIMPQSWZ\^agnswz|~ćЇ܇އ +  + + + + + + + + + + + + + &jхbֆ;̆ !'))+-/0136:=?ADHMPQSXZ\_dipuxyz}ɇ҇ه܇݇  + + + + + + + + + + + + + &jхaԆ;͆ + (.,*,-...159;=?@FLQQTWZ\`ejouwy{~Ƈ̇҇ևه߇ + + + + + + + + + + + + + + &jЅ`ӆ:͆ "-40++,,+,059:<=?EKOQSW[_bfjotvx|Ňʇ·҇և߇ + + + + + + + + + + + + + + + + + &iυ_҆9̆ %-1.*)*++,059:;=AFJMQTX^cfhjnsux{‡Ňʇ·чՇ߇ + + + + + + + + + + + + + + + + + + + + + &iυ^҆8ʆ%(('&')+,-0478:>CGJNQUZ`dgiknqtvzćɇχчՇ߇ + + + + + + + + + + + + + + + + + %iυ_ц6Ȇ!%$"!"%(+-/0357:?CGKNRW]acfhjmpruy}ɇ·чՇ݇ + + + + + + + + + + + + %hυ_ц5Ȇ! #&)-/0369=@BEIOTY^`cfghknquw{LJ̇χԇ܇ + + + + + + + + + + + %g υ_І4Ȇ"$(+-037;>?@BHNTY]`cefgjlptvzŇʇ·Ӈۇ + + + + + + + + + $g ΅_І4dž  "'*,/37;>>?BHMSX\_befgiknruyćȇ͇Ӈۇ  + + + + + + + + +  $g ΅^φ4dž  #&),/37;=>@DHMRVY]adgiijlptx}Çȇ͇ԇ݇  + + + + + + +  $g ͅ]φ4ņ + !$&(,047;<>BEIMRVWZ_dgiijlnrv{‡Ƈׇ̇ + + + + + + + + +  $f ͅ]І5ņ !!#%(-048;=@CEHMQUWY^befhjkmqvzćʇև + + + + + + $f ͅ^ц6ņ !#%)-/27@ADGINRTUV[^_`beimqvz~؈/ + + + + + $g ͅT*† +$)-37:<>@BDEGLQTTVY\^adgilquzԈ   + + + $g ʅO%  &+047:<>@ACEGJNQTUX[^bfgglqty}LJۇއه܇އ߇$ + + + + $g ʅQ( + !',0358;=?@AEHJKPTVWZ]adegkosw{}ȇ͇чևڇ܇އ8: + + + + $g ̅U†, "(,.0258;=>AEIKMQSUWY]`cehlorvy{}ćȇ͇Շڇ܇އ?9 + + + + $g ΅Ydž- + #(+,-0247:=BFIJMOQTWY]_aeiloqtvxzzyÇƇ̇ԇڇۇ܇އ % + + + + + $g҅]Ɇ- #'((+./148CHJMPORXZ[]_adgjkkkmrux{}~LJЇԇԇԇׇ؇ڇ܇߇߂ + $fƅS%zކ !$&()*,159>DHIOSQQUWY\^`cehjkmptxz|}Ƈ͇ЇЇ҇ևׇه܇߇܂ + + $fŅS†&{ކ  &()*+.26;@DFHLONPRV[]^abcegjnswz{|}Ňʇ͇̇ЇӇԇׇڇއق + + + $eŅS†&{߆ + %()*,/38?CCDFHKNPS[ggbaa`beiouxz{{|‡LJˇ̇·Ї҇Շ؇܇ق + + + + $eŅS†%x݆  $'(),/39@BABDGKOQTe|yg`__`chpuxyzz}‡ȇ͇̇·чԇׇڇއ߇ڂ + + + + $eŅT†$x܆ +  %)''+/39=>?ACGKMNScywf_^^_ciosvwxzŇʇ̇·чԇ؇ۇ݇߇ۂ + + + $eŅT†$yކ !),(',05<=<=@CGIIKPYcd_\\]_dkoruwwz‡Ƈʇ͇ЇՇه܇݇އڂ + + + $fŅS${ +  )+''+.3:<::@BEFHJLMQTVX^beinprw}‡ȇˇˇʇˇ·ч˂ + $g ɅXdž)~߆ +!#$'+.0357;=>@CEFHILQSTX]behlosx}ŇLJȇʇ͇·χȂ + #g ʅXƆ(} !$)-/0247:=>@BCEEGKOQTY^cgilpuz}‡ćLJ̇χЇχɂ + #f ʅYƆ(~  %*-/1258;<>@ABCEHKNQUZ^dgjmrw{ÇLJˇ·Їч˂ + "e ̅Ydž)߆ + !&),/1369;<>??@BEHJMRVY]beiotw{~ÇLJʇ̇χч̂ + "e ͅ[Ȇ*݆ "&),/1479;<==>@BEHJNSVX[_ciptw{~ÇƇɇ͇χʂ + "e ΅\Ɇ*~ۆ #&),/1589;==<>@BDGKPUVWY]bjpsv|ҡ·Ȃ + "f ΅]ʆ*}ن #&(+.1467;@@??ABDGLRWXYZ]binrv{桷Ȃ + !e ΅\Ɇ*}؆  $'),0235;AC@?BDFINSWY[\^afkptx|ćɇˇł + !e υ\Ȇ){؆"%'*-.16><=ADIPUX[\`aa`bdhkquvwxyz|‡‡ +  e ΅]ʆ*|׆ +!$')-48:;;=ADGMSVY[^``abfilrvxxxwx{~ +  eυ]Ɇ){׆ "%(,14689BCDGKPRQQSUX[`cfilnpruy{zz|}~ +  fЅ]Ȇ'wӆ +  #&)-39>ACEHKPSRPQRV[_bfijloqswxxz||}~ +  fЅ[Ɔ%w҆ "&*.4:>ABDEGLPPOPRVZ]aeggilnqstvxz{||~ + + fЅ[ņ$vцކ "&*/5;<>@ABDHKMNPTY[]adddfilopqsvxy{z| + eЅ\Ɔ$uІކ  #&)-4889;>@BEIKLPY^\\_aabdgjnopruwyyyz} + + fЅ\Ɔ$tφ݆ !#%'+03568BFGFJPTUVX[_bdgiklmopruvxz{} + + +  gυZĆ"r̆ۆ߆ "#$')+.1347:=ADEEHKORSUY^beijiikmnpruy{|}} + + + +  gυZĆ!rˆنކ !##$'(*,/1248ACDEGLPSW]`begggilnprstxz} + + + + + + fυYÆ odžԆچކ + #$&&&()+-/16:=?@ACGJMRY]^`bbbeimprsttux{}~ + + + + + + f΅Xmņцֆۆ #''&')+-/048;=>?CFGJQWYZ^`_`cglprsttuwz|} + + + + + + f΅XlÆІֆۆ߆  $%&')+-./369;<=ADEIOSTX^a``beimprsuwxyz} + + + + + + fͅWk†Іֆچ݆ +!$&)+,-/3689:ADFJOWbiga^`dgjlnqvxxxz~ + + + + + e̅UjΆӆՆن߆  "&)*+,05989;>@CEIP]msi_]`ceiloruvvxz + + + + + e ʅSj†φӆֆچކ + $(***.3789;>@CFIQ_rwi__abcgkostuvx{ + + + + + e ɅQg†ІՆ؆چކ + "'***-2579;=AEIMQZgkd`bbbbfinrsuwz| + + + + + d ɅReφՆ׆چކ !$(*,034679?CGLRUX`ktupllmmnsxyyz|}~~~ + + + + +eɅPdʆ҆؆ކ #).23348::;=@CGLQWdwunmlknsxyzz{|}~ +eɅOdɆц׆܆ "(,/2224679<>ADHOWgxnkjknrvyyz{||~ +e ȅPcȆІՆن݆ + + $'+/0/02489CLU_rzjfghjmquwxyyxz|}~ +eɅOaɆφ҆ֆچކ "$&'')-2788:=BKSZalunedfgilortwyxwwy||}~ +f˅Oadž͆цԆ׆ۆ߆  "#%%'+0468:=BHQXYZ_``bdfhjloqtwwwvxz{{|}~ + iЅPbƆ̆φ҆Ն؆܆ + !##%),/49<>ADKRSSUY\`bdgijlnqsuvuvyzz{{{|}y + i҅RbƆʆ͆ІԆ׆ۆކ !#%&(,39=?ABFJKMPUY]`bfghiknpqrsuxyyyxy{zs + hυRaņɆ̆Ά҆Նن܆߆  #$$$*29=?@BFGGIMQVZ^adeggilnmnpsuvwwvxzxp + g˅P_†Ɔʆ̆φ҆ֆچ݆߆ +  "##$)08=?@CGIIJKNSX\`bcdfhkkjloqsuutvyyun  + gɅN]ĆȆʆ͆ІԆ؆ۆ݆  "#%)/7=?ADGJKKKMQV[^`bdfiiiilopqssruxxtl  gDžL\ÆƆɆ̆φ҆ֆن܆ !$&(.5;@CEEGIJKLOTX[^adfhhhikmopppqsuusl f ÅH YĆȆ̆φӆֆ؆ۆ߆ +  $&(-39?BCCEHJKKNRUX[_bdeefgikmmlmprqrrk dAW†ƆʆΆӆֆ׆چ݆ "$&+17<>?BEGJJJLPSUY]`bbcdegjkjiloonooh b4SÆdž̆цԆֆن݆ +  "%*0599;@DFHIIJMPRUY\_`abceijhhkmlllmg `$J†ƆʆΆцԆ׆ۆ߆ + "&+0357:>BDFGHIKNPRTWZ]^_adhhgikkjjjjh ]|:~†ƆȆˆφ҆Նن݆ + #'+.0259?@BEGJLLLNQUXZ^dgeccdfggfea RӄgL +QÆņȆ̆φцԆ؆ކ +!$$&*,.0369;<>?AEHJKLMORVZ_dfda``cdefc\ KÄNŅ05jÆdžˆΆцӆ׆܆ !#"$)+,-0357:<>@CFGIKLMQTX\`a`_]]_acec\ +C4tΆMuƆɆ͆Іӆֆۆކ  "&()*-02479?AABGJKLMLMOQQQTZ^_``Z \Ʉ=f7Yy†džɆ̆ΆцԆ؆چ܆ !$(+-/0368:=>?BFIJJJIJLOQQRV[\\\U < +xޅ:7[{Ædžʆ͆І҆Ն؆ۆކ "%')+.1236:;ADFDCDHNTUTSTUWVN/ze̅V6ZyÆɆ̆Άцֆ؆چ݆  !"%')+.26778;>ABCBCFLTYYUSTVUN +=߅!W6XwĆɆ̆ІԆֆن܆ކ !$&)-1455689;=??@BIS[ZTPQSSK F"V 3WxĆȆ͆ІԆ؆چۆ݆  $',/233466579:FMNLKKLMJ + J#T 3WxÆɆΆ҆Ն׆ۆކ !%''(+-.///024446>FJKJIKNK + J!S 2UsņʆΆ҆Նن݆ + """%')*+,-.121028AGIHHLQP + J S 0StĆʆІԆ׆ۆ߆  "$&'()*,.//.05>FGGHLRV + + J!߅Q߆1Vwdžφӆֆن݆  ""#%'(*,--./2;CFEGKQY + K"ޅPކ +3Ww͆҆Ԇ׆چކ  "$')*+-//16<@CEINS + L#ޅN߆ +2VwņΆ҆Ԇ׆چ܆݆ކ +  "$'()-./025;@CFGE + M#݅O߆ +2WyĆΆцԆ׆؆ن܆ކ  !"$'+,,-/28>@BCFILOQOJ + &`Ȅ7*KiDžхمޅ߅߅߅ XTENSION= 'IMAGE ' / Image extension BITPIX = 16 / Bits per pixel NAXIS = 2 / Number of axes NAXIS1 = 132 / Axis length NAXIS2 = 288 / Axis length PCOUNT = 0 / No 'random' parameters GCOUNT = 1 / Only one group BSCALE = 1.000000E0 / REAL = TAPE*BSCALE + BZERO BZERO = 3.276800E4 / EXTNAME = 'SCI ' / Extension name EXTVER = 2 / Extension version INHERIT = F / Inherits global header ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator DATE = '2010-07-02T19:55:22' / Date FITS file was generated IRAF-TLM= '2010-07-02T19:59:26' / Time of last modification CTYPE1 = 'RA---TAN' / R.A. in tangent plane projection CRPIX1 = 65.7568549170406 / Ref pix of axis 1 CRVAL1 = 306.871725219222 / RA at Ref pix in decimal degrees CTYPE2 = 'DEC--TAN' / DEC. in tangent plane projection CRPIX2 = 143.564898774826 / Ref pix of axis 2 CRVAL2 = 37.3733845984976 / DEC at Ref pix in decimal degrees CD1_1 = -2.9278046959101E-4 / WCS matrix element 1 1 CD1_2 = 1.37000324401069E-4 / WCS matrix element 1 2 CD2_1 = 1.36969765153489E-4 / WCS matrix element 2 1 CD2_2 = 2.93069759666565E-4 / WCS matrix element 2 2 RADECSYS= 'FK5 ' / R.A./DEC. coordinate system reference EQUINOX = 2000. / Equinox of coordinate system MJD-OBS = 52162.6261539214 / MJD of start of obseration CCDNAME = 'EEV 9273-20-04' / CCD chip name(s) AMPNAME = 'EEV 9273-20-04, right' / Amplifier name(s) GAIN = 2.337 / Amplifier gain RDNOISE = 3.3 / Readout noise CCDSIZE = '[1:2048,1:4644]' / CCD chip size CCDSUM = '2 2 ' / CCD sum CCDSEC = '[1:2048,1:4608]' / CCD section(s) DETSEC = '[2049:4096,1:4608]' / Detector section(s) DATASEC = '[33:1056,1:2304]' / Data section(s) BIASSEC = '[1:32,1:2304]' / Bias section(s) FRMNAME = 'dhsFrame' / Name assigned to a frame of data. FRAMEID = '1 ' / Frame Id DATATYP = ' ' / Type of data frame WCSDIM = 2 LTV1 = 0.4375 LTV2 = 0.4375 LTM1_1 = 0.125 LTM2_2 = 0.125 WAT0_001= 'system=image' WAT1_001= 'wtype=tan axtype=ra' WAT2_001= 'wtype=tan axtype=dec' END wybiƅDžDžʅЅ؅߅ׅυŅ}ywxywuuwwutvwurdy{pŅ…………ąƅȅɅ˅ͅ΅ЅӅՅ؅܅݅ۅۅۅۅۅۅممڅ܅݅݅݅ۅՅυʅDžDžȅʅ˅̅ͅυЅЅЅх҅҅Ѕ΅ͅͅͅ˅ʅɅȅDžƅŅąÅÅz|{؅܅߅ + +߅߅݅ۅڅمׅօՅԅӅхυ̅˅ʅʅʅɅDžDžɅʅʅDžz|ȅ  !  ""#%()&#!  +߂z|Ն !#(*(&'('')+13.*)*,---.-,,,,+++*+--.142012.*)+,.24321/..//.-+((*)('&%$" z|܆ +  #$#$$&*,--.-,,-4EL=1-.022343221/---.0457:<:8::633468<==;:86775444101000/.-+)('%%$!  z| !"#%'''''')+.10//./9R[E3/024467776410012479;>?=;<=:89:;>ABA@><<<:87876443443110/.-+*'%$"!z| "%&')*+++*)(*-0//0/07GL>3113568:;:9644567889==>?@CFFDCBAA?=;::987679975543110-*)(&$#!z| +!&)++,,----+)*,../00049:6333468:<<<;:99:<<:99<@CB?=?@@@BDFIIGFEDCA?>=<;:98:>?=;:85453/,,.-)&#!!!!"""!  z|  !$)-...../0/.,---./012355555579:<>??ABA>?@?<;;>ACCAABCBBDGIJJJIGFDCBA@?>=;;<:8862.043.*'$$#$&%$$""##"#$#"""!z| + #$$%',//0122243100/./23345678889:;<>ABEINLFCB@???ACEDDEGGEDFILKKKJIGFFEEDCC@>>>@DFDA?><:8425:71-+)''()''&%&)(''('&%$#"z|  "&())+.012356799763224677779;<<===<=@DGJNSQLHEDDCCDGIHHJKKIIJLNNMMLKJJIHHHHECBAACGHEDCB?<977;>:3/.,++++)))*,-+*+,+)(''&"z| #&)*+,-0245689;8679:;<<<>@@AAA@?@CGJLNQROLJIIGGHJLMMOPPNNNOPQQQONMLKJJKJGFEDEGIJIHGDB?===>>;631/.//-,+,/0/--//.-++,+&z|""$'+-...03689::;>BBB@=:;<>?@@ACDDDDDDCCFILNOQRQPPONLKKLOQRSTUSSRRRTUUTROMLKLMLJHGGHJLOOMHEDCCC@>=<:64223210/133101210/..00*z|!%&(,/023347:<<===>ACCCA@??ABCDFFGHGHIHGGILOPRSTSSUUTQPPQSTTUWXWWUUUVWXXWSPMLMONLKJJLNRYWPJFFHGEB?>?=:7676653356666665431143. z|$(*,0356788:=>???@@ACEFFFDCCDEGIJKKKLLKJKMOQRTUVWWXZYVTTUVVWXYYYYXXXYZ[[YWSPNOPPONNNOQX`]SMJJKJFCAAA?=<<<;::989;;;:9998754550%z|'+-035899:;=@AAABCCCFGIIIFDEGIKLMNNMOPNNOQQRSUWYZZ[\\ZYYYYZ[[\[\\\[[[[[ZYXVSQQRTTRQPQSW\ZTQNNOMIFEDCBA@BB@@@???@@?<;;;:976652&z|*/13579::;=?ACDEEEEFIKJIHGGHJLMNOOOPQSSRSTSTUWY[\]^^^^]\\]]^^^^___^]]\[[ZYWUUVW\]YTSTTUUTSRPPPOMJHHGEDDFGFEEDDDEEB@??=<;:8773(z|.3567:<<==>@BEGHHGGJLMKHHIKLLMNPPPPRSTUVVVUVXZ[]_`aabb`^^_``_``abb`_]]]]][XVXZ]bd_YWWVTSSSRPOOOPNLKKHGHKMLKIIIIIHFDDCB?><<;96,z|3679:=?@@??ADFIJJIJLNOMJILNOONOQRQRSTTWXXXWXZ\]_bccdedb``abbaaabccb_^^_`_][Z[^`bda]\[XTQSSRQPOQRPNNMKKMOPONNNMLLKIHGGEB@??=;80z|" 7::;???@BCCCCDEGHIKOPPQTTRQSTSRSTUVWXXXWWY[\ZZ[^aceghhgfdcdeeeedeedca````bbcddeggffdb`_^[XWVWVVWWVTRSSSTVXWUTTSPPOOOPPNKIGECA@>3z|%>CDDDDDDDDEFGHHJMQRRTVWVUUUTUUUVXYZZZYY\^^\[]_bdegjjihgeefghgfffdcbaaabceggfhjkjjgdca_^^][ZXYZZYXWVWWXYZZXWVTQQQPQSTTROLIEDDB8z|&AEFGGGGHGFGHIJKMORTVXYYYXWWWWWVWY[[\\[]_aa`_`beggilmllkigghiihhhffeccddehjjikmoomjgfcbbbb`^\\\\\[ZZ[[[[\[ZYWVUTSRSU[a`XPKGGGE>z|(BEFHIIJKJIIJKMNOPRVZ\[ZZ[[ZZYYXX[\\]^^_`abbcdfhjkklmnopmkiijjjkkjigffghikmmnprsrokjgeddeeda_____^^^__^^^][ZXWXWVUUYdro_SMJIIICz|*AEGHIJMOMLLKLNPQQSX\^]\]^^]\\[ZZ\^^^`abbcdefhjlnnnmopqrqnlkkjknonkjijjkmopqsuvvsoljhffffgfbabbbabbccbba`_^\ZZZYXWW[hvraUOMKLMFz|,AEGIJMPQPNMMNPQRSVZ^``__``^]^]]]]_``bdeefgikmnoqrrssrrrrponmlnqrqnlllmoqrrtuwxxtqmjhhhhhigdccddeghhfedaa`^]]]\[YXX[bhf\UQOOQQJz|-BFHILNPRRPPPPRSTVY\`cbaba_^^__^^^`bbdffghjlopqqstvxxvttsqpppqstsrpnooprssuvwyzyvsoljihijkjgeefhjlmkigdaa`]\___\YXZ[\]]ZVTSTVVOz|-DGHJLNPRSSSRSUVWY\^beeddb``a`````bcdfhiklmorsstuwx{{xvvtsrqsvxwusrqrrsttuwxyzzywuqnkjijmnljihjloqpnlhecaa_^_`a_\\]\[Z[[YWWXYYQz|.GIIKLNRTUUUTUWY[\^adfeeedcddccbbbdefhknppqtuvwxxyz{{zxwvutsvyzywvvvuuvvvvwy{||{yvromklnqqnmmmnruusrolifecba`abb``_]\\\\[ZZ\]]Sz|/IKKLMORTVWVVWZ\]^`bdeefffffeffeddffhkrxyvvxyz{{zz{|}|zxwwwxxz|{zzzyxwyyxvvz}~~|yurpnoqsttrqqrsvyywuromkifdcbcddcba`_`__]\]_dg]z|1IKKMOQSUVWWYZ\^_`bccdefgiigfghgffghjox|z}}}~~}~~|zzz{{{|~}|{zz||{xx|}zvrmlqvwwwvuvvxz{{zxtrpnlifdefggedccccba```bksjz|1HKLOQTVWXXY[]_abddccdhiikkihiiiiijlnsz~}}~}~~}|}}{}}|vnmtyz{{zyzz{}}}|{wtrqpmjhiiihgfeeedddccbdmvlz|2JNOQSVXZ[[\]`bdefeccgjlmmmmmllmmllnrw{|}~~yxz|~}}~~~|yvuutqnnnlkkjiiihgghgfeflpez|3NRRRTWY\]^^_befffeefhjnqqpppppqpooqty}|}}{yyywtssrqqponnmlkllkjijmmaz|4"RUTSUXZ\_`abdgghhhhhikpttsrrrstsqsuw{~~}|||zxwwwzzvttsqpoooomklmk^z|6%UXWWXZ\^acdefghjllkjjlrvvutttuuuuvx{~}{zy|{yxvtrrqqqommlh^z|8'X[[[]_`aceggfgjmoonmmotwvvuuuuvwxy{~~||~}|zwusttsrqonlh^z|9'X^^_accdeghhhiloqrrqqruwvvvvvvwyz{}{wuuvwusqpnlgXz|:(Y_abdfgghhijjkmoruvuttuwvwwwxz{||~{wvvxywtrpmkfOz|:)[acehkkjjjjjjlnpsvwvuuvwwyz{}~}zyyzxvtronmiWz|;+]cehkoolkkklmnprtvuuuuwyz{}}|||yvtqooolbz|=-_dgjnsspmmnopqrsuvuuuvxz|}~}{wsppppofz|>/`egkpturpqrrsttuvxwvwxy{~~|wtrrsrofz|@1bfilqsssssstuvwwxyyxyz{|}zwvvutpez|A3dikmpqqrsstvxz{{zzyyz{|~†}|{yxvrfz|B4fjlnopprstuwy{}}||{{|~~~}zujz|D6glnprqqsuwxxz|~~~~}wmz|G8hmqstsstwzzz{~zpz|H9iortuttvy||{|}~}qz|H:jpsvwwwwy{||||~~pz|Ipuwy{{{{{|}~}jz|L@rwxz{{{}~~~{kQz|NCty{||||~oY=z|PEv{~~~~mY<̂z|RHx}nV8ƅpz|SJ|̆҆~jQ5e8 z|TL~Æ܆̆p[LIP[fv}gJ*䅿`.τz|VO††φԆņtR/2Li|eG#Յ[+ńlDz|YRÆÆÆĆĆĆƆ†††††ÆƆÆÆzX0؅Յ;_yv^A˅sFe@!z|ZT†ņƆdždždžƆƆƆņņņĆĆĆĆņƆƆņƆÆÆÆ}]5ՅrNBJ`}ȅ75;I`}ńz|b^†ĆņņĆÆÆÆĆdžɆɆˆ͆ІӆՆ׆׆ֆԆӆԆԆԆԆԆӆ҆ц҆ӆӆ҆цІІ҆ӆ҆цˆWޅTτ^1 䃻fB+-B_@_oodWOKNWg~Ą<`wz|daƆÆĆņƆdžƆƆƆdžʆ̆͆φ҆Ԇֆن܆܆نֆֆ׆׆׆׆׆ֆԆՆՆՆՆԆԆӆՆ׆׆Ն҆Æf#ۅ<鄙O؃w^E,#5Oy7Q\]_es„:fӅz|ebņɆȆņƆÆņƆdžɆʆɆɆ̆φІ҆Ԇֆ؆چ݆ކ݆چنچۆۆۆۆۆن؆نننن؆؆؆چ܆܆چӆ<셗;܄}!σ[;.**$*Q#IcwƄ6^݅7Dz|fcĆʆ͆̆ɆdžĆÆÆĆdžɆʆ̆Ά͆ΆцԆՆֆ؆ۆ݆߆߆݆܆݆ކކކ߆߆ކކކކ݆݆܆݆ކކކӆnX\ $"!;r%[ 6`օ#BZm{}z|geÆɆ͆φІΆˆȆƆdžʆ̆ΆцӆӆԆ؆ننۆ݆ӆ_(9V܂ȂĂԂ),-+%9qCԅRׅ?\uz|hg†ņȆˆφӆӆцΆ̆͆ΆІ҆Ն׆؆ۆ݆ކ߆ֆW􅄅 !܂Ȃ "-32-$*H̄t̅'}Å>[sz}ih†ĆņƆdžˆІӆՆԆӆ҆ӆԆֆ؆نۆކ܆Z􅁅Ђڂ)13/''7[G sՆ%Wtz}jhĆņdžȆɆˆІԆֆ׆׆ֆ׆نچۆ݆ކ熺iԂς *-+&&=joڅEdz}kiÆĆƆȆʆˆ͆ІՆ؆نچچچ݆߆߆Ն%957Ƃ&)'+G̈́,q؆6~z}ml†ĆƆdžʆ̆͆φцԆ׆نۆކކ߆ 􆲆O߅eee˂؂ +(5>KmZ&@z}nmÆƆdždžɆ̆ΆφцԆֆ׆؆ۆ߆  ֆs"ۂ‚‚ւ,Hax:FDz}mm†ĆņdžɆˆ̆̆͆φцӆՆ׆نۆ܆ + + + (>ȄN˃MЂǂɂԂ=fЃC΅n +Mz}moĆƆȆɆˆ͆Іццц҆Ն׆؆ن܆߆   #ņVt) -Z9\~ۅ`&_ÆÆz}pq†ƆɆʆ̆ΆІӆՆԆԆֆ؆چۆ݆ +   (#DτX僈PAOlÃ+[҄CwM{ĆdžƆ†z}rrĆɆˆ̆ΆІӆֆ؆؆نچۆ܆߆  +  !+0\&[ +׃ƃЃ@rׅ,Qs߆0V{ÆȆʆdžĆz}rtÆņĆņȆˆ͆φцԆ؆ۆ܆݆ކކކ + !$,2"솝E慇-ۄm^czʄ-^օ/Jc{ņʆˆˆɆƆƆz}svĆƆdžɆ̆͆φӆ׆نۆ߆ #&&%$%&&&&%%%%%%$" !#%&'*.)܆?텡]'/X7Umdž̆ΆΆ͆ˆɆȆņÆz}uyÆdžɆʆ͆φІӆ׆ۆ݆߆ +  ! !##$),,--.///.--,+***)''')***+/3.؆P܅:bʆцӆ҆҆ԆՆԆӆ҆φ̆ˆȆņĆÆz}v{†Ɔʆ̆ΆцӆՆ؆ކ + !$&(''(**+,*+034677667766410/00///100005=?/ ކ{S8,&!(@`ц݆߆܆ۆن؆׆Ԇцφ͆ʆdžņÆz}vz††Ædžˆ͆цՆ؆ۆ߆ + !$'*,/11124565349>?@@?>=>?@@=9657766788767;AB;+ֆֆކ܆چ؆Ն҆цΆ̆ʆdžÆz}vy†ÆĆņƆʆ͆φӆنކ + #%),/2469:<<=>?>=<>DIKLKIGFGIJIFC@??>==>?@@?@CDB=6-"ކۆ؆ֆԆцІΆ͆Ɇņ†z}vyÆĆĆņȆˆ͆Іӆ؆݆ #'*.37;=>@CDGHHGGFFGJNRUVVUSQRSSQOMKJHEEFFEFHIJKHC@<83-&"  ݆ۆنՆӆӆ҆φˆȆƆÆz}x|†ĆņƆȆˆΆцԆ؆܆ "!"'-159=ADEFHLOQQOMNOQSWZ\_`aba_^^\ZXVTSQOPQONOQRRQMIFC@;730.,& ߆܆؆׆׆ՆцΆˆdžÆz}|ÆƆȆʆ͆φ҆Ն׆ۆ߆ "$()(+18<@DILNMMPTXYXUSUX[_dghjklmllkigda_^]\[[[YYZ[ZXWSPMKHC>9531+# ߆ۆچنֆ҆φ̆ȆņÆz}ÆņȆˆφ҆Ԇ׆ن܆ކ +$*,-.04:?CHNUXXWWY[_aa^\\`ejnrtvxxwwwurqpmllkifedcdeeca_[WUSNHC>854/(" +߆݆چ׆ԆІ̆ɆȆdžĆz}ƆɆ͆ІӆՆ؆ۆކ +"*/246:>AEHMW]`abbbcgijkkjkosv{}{{{z|~|uomklnoomjgc_]YSNHB<861+'" ߆݆چֆцΆ̆ˆɆņz}ÆɆΆ҆Նֆ؆ۆ߆ $+26:=@DGJMQW^adgjllnqqsy}|zzz|xttvxxxvsokgd_ZUOIB<950+&"  ݆نԆцφ͆ʆdžz}ņ̆цֆنچ݆  &-8>@BFINQRV\addfkorvy{{{xuokga\WQIB>;50+'%" ߆ۆ׆Ԇ҆ІΆʆz}ȆφՆن݆߆ +!'-5@HJLPTX[]`eijkmqux~}xurjb\WQKFA;61-+(&! ކۆ؆ֆԆ҆Ά‚z}ˆ҆ֆ܆ #(+/5;CLQU[`dgikorstvz|~ukc^ZWQIA=9531-'" +ކ܆نֆՆ҆Ȃz}…ΆԆ؆ކ &/469?DLSY`hmrvxxz}~‡·‡yphda\SJEA=;83-($ ܆ن׆ֆ̂z}ÅІ׆܆ +#*2:>@CJRZaelsy‡LJƇƇчه̇pje]TNIEB>93.)# ߆܆چ׆΂z}ąՆۆ +&/7?DGKPW_gnrw}ćЇׇ݇݇χʇ͇Ӈ҇χ҇և̇ʇˇɇɇɇƇLJ׈4?҇smg`WQMJE?:4.(" +݆نЂz}ƅن߆ &/9AHMQW]cjqx}Ƈ̇χ·ʇȇ͇Շ؇߇ ݇҇·Їч͇̇ˇŇƇч؇ڇهևևׇׇԇчχχ·ˇŇɈBKyc{tqk`YURMG@93-(# ߆܆҂z}ƅ܆%/7@GOV]dintzÇׇ̇'&҇Їڇ߇ևʇÇćˇ҇ه߇߇ڇчƇÇ׈Xs{|ztjc^ZTMF?93-(# ߆Ԃz}Dž݆ &08@GOWbjpw~ȇӇ߇0/'%'#ׇهۇ݇ۇчԈ{|tmhb[TMF@93-'  ӂz}Džކ +%/7?GPYdow}Ň̇ԇ߇ +.AQZ^K)):;8:=>3'B5 + هڇć|vrjb[TNF>82+% ӂz}Dž߆ +$.6=EP]it~χۇ߇߇BeֈdKELQUY_c\G53>JJDHS@ + + ۇ҇ɇ{tkc\TLD>70*" +Ղz}Ʌ !-7>EN]n}Çч߇ +Bp։"q͊2_ӈi`glq~{icaZWbwY1   ܇χ}skd[RKD<6/& قz}Ʌ*5?FLWgy·݇#3:04ZƈX܊.Fψow}spbP9   ܇ŇzskcZSKC<4,%! +z}ʅ '2>HNT]j{χ /Rbhlsx|ԉF+|gوxMY׈ڈڈш#w>!*,#)'؇~wnd[RJB:3,)'" z}̅ +%/:FO[hqsżއ.Ga~LjBljKӈr>JΈ&K`XHCN_fW3[:/,,7:& #/1)҇wnd\TKC;51,'! z}̅#-7CNWcvч %FdĈ̈ۈ=zۉʉ6ۈZw-BGEE_܊߉T{I;=?:, 060#чvng_UKB<72,&z}ͅ *4>KV^i}ׇ&CdˈوΈȈ̈وHwWوщ'g}r~!q$߉͉ ^(҈]TF4)  .1* ӇÇzsi_TJC=61+# z}΅‡ &3>FR^fnӇ*Eb}ψ߈ԈˈĈĈۈ11ƈӉ,̉}<"芛Q -:$򉱉k!ӈwX9*%$.,$ۇχʇƇ~sg]SKC<6/'! z}ЅLJ"-;GOWcluχ ,AYm|шΈɈ܈ވΈՉuΊ0- +77{QӊK2@VGa)戦sK1(%!!,51&"$ ݇ˇ}qg]SKD;2*$ z}хˇ +)4BOV\gr~Շ):GWgqԈՈEщ,*!,OxGE8T[&ωv;创jF0(%#,96)$+,! +ӇŇ{qf[RJ@6.'  z}Ӆχ &2>LX`gq|և '05BYp͈шˉKɉ;qċh8ߊ\U>͉2∦vO4%#/1,1<;/*'Շȇɇ̇ˇŇzodYPF;2+$z}Ӆ҇+8GWdlu͇ *NzՈ݈ɈʈB׉1fڋ-siN퉸p-󈻈P/ &/@OMFD=/͇݇̇Ӈׇԇɇyl`UKB80)! z}ԅՇ$1>Pesz‡އ $Jwڈ و҈'_|rr։Dyދ*cwNJ7ωN!vL2!".CSRQUK:+҇·ׇևȇtg\RI@7.&! z}օه#4BVm{҇ !! *Gi͈҈҈؈߈ۈ׈ 1EGH^DhJn􊉊S/Ή?׈|^C1,2@JIKOF92#ۇчه҇ć~pe[QG>4,(%z}؅܇4EWjyȇ$10'*;PkȈ͈ɈȈψшΈ[iamy}s]LB?ETp.GIRQ<]ibF݉( шqQC@@BCA>4-2- هއχÇyodYND:2-)#z}څĆ'9IXixׇ0L\[Yaq҈݈߈݈3}v[$ É݊SGtmt8hzNE!*ՈwbQB?C>4(#.3؇̇wj_ULA7/*$z}څƆ%2?M\m}̇'KvĈ &9N[mYŠŊwF?FLOYaZK0`يފ4DҊNJL)2*"Ȉz_GBE>6.(,2#އՇ̇~od\TI>3+' z}څƆ)6CP_p؇ "=fՈ؈ĈɈHpЊ=^H=VpS\fLIhm$N3'VagTA$爪fPLICDG@76+ ܇ՇʇwkbZRF90+$ z}څƆ,:FRasˇ&&'4Qyڈ B܉܉؉ ҉ԉ5.'z}څdž1?HTdvه)-+3EaĈA4  ̉ϊ_9Ќbы;Dec)͊}9@5]JcW/xN6׈xiehr{u_J2݇χ~qf\OC90) z}ۅȆ4BKWg|̇ ...?Wr܈KNr:"É͊cԋS*׏X3֌Ӌ_B_e@$^0!m'ӉcS>̈oP- ##! +ׇȇznbTH<2*!z}ۅɆ6EO\nه&/3MXfx҇!+B_t(Xz 'ߊ 1976Fo5ߍ؎C'7j<ËYD̊cmFp=Ԋ슟;։ՉQ6=)툸xvvfQ>5 E,!抵fo V̉ъՉtPSJ؈zcE;BC/чʇLJ|l^QE8,z}߅ц@P]l‡"5Qfov{ՉAcN-) RxFʍ+oώЍ׍[ˋD势s@%ӊnɉȊ6ɉzmc>\>@D-؇͇͇ȇscUG9. z}ޅц@O\nɇ#9Qi~# Չĉ݉#! +eŠ̊\4 +]"~\ѓRď`Z\T97B>O#銄+މԊ$DɉpX#ማn_nmME;߇߇܇هˇyhYK=1 z}߅ԆAN[ȯ&&z}؇J\m߇(hщS҉tVٜݏpZJBdЊĊy|`ωN +͈XFLK(ۇ +·reWG2z}؇I[lEV}{pgX;8dO + +#nیQ'a^2H‘2̫2O;D(̌ƌXc#ъ⊳r)74􉵉ljω^㈤hMYgD҇vhXG4z}هL^oHMhg\LA=<9:EQ_tuknĊ fNtc5 ˙NN.K~UƵ]bOLы^#Њҋ +芴^=R8쉴 ߉}%\[d<߇އևxgUE2z~ڇM_oهE2SU@$$JbeanҊk (L?aWy _7p:ߑHdWŋZߊʊ&*Z3+xv<;׈mN?ڇև‡xeRC1z~هK[k҇EӉ 6@%6NSF8@RWYg~,iԊ!%5ZȌL'4]a Ѻ5.L3Y)o9Ïj= ɋCϊ݋;FŊ`kъ]׈sE' ׇćuaOA1z~ׇGVf|͇C7C?609JX]j:mʊ*=9@qDϛ0 l@}1jd5Y⌊ݍ+{GH͊MƉԊ<劔vÈoI.36'ևÇq\M@0z~ՇCSdyʇ@xɈ݈݈؈݈ %43.'3RbY[;k,UIC+t2şjWvLm-,ƭ2:ΎӍ<@印_ `Z-g{XHs􈪈oWNK8.DYVE5(ԇlYLA.z}ՆBRcwLJ߇8b|ψ*3+%$4PZOS/^ȋNE8W΍@񛺝9Ʀuǵ UݷX<,;dJ ׎tŊRce"ʉ&8, ʉV׈~unxseiiZF=1Ӈ~hVMC.z}Ԇ@Qbu‡ׇ&=DQf|78&+;EEUJ{38+HVǐ̖מנMħ( CgǓVˎ}ԋ!.~҉ 'rhmvd?),+ ӇuaRKA-z}ӆ=N_r·և߇(Lt߉$;8 *39@[ފ +1jĊf@ +: X?ަ$鮜sBW!xϒ)+Ŋnˊ:׉Ɋ + Љ҉ ۉhꈈWS]paЇjZPK@+z~߅ц:JZnćʇԇHtш42%6:7;ZӉ ]NJ+ދzțȞWϡ졀}5t'ېT Ίt?Vtω֊܉뉟dJKTlڈ` ˇzfZRK@.z~߅І7FVjÇ̇܇Dh| +&+0BID=PÉ"eŊъڊ9_Je=2KlͣN'b֗ !񊌊k ՉkӉΉn툉[QTZlڈS Ƈvh`WMA/z~ޅΆ4CSf{ч;Td|Έ(!!-=OTKLm%pҊ݊ tښb Ξ&H斟+ALE ->ԋ%kMϊJ ≥h\ۉOڈgaeinrwq?$ +܇ulcYMA.z}݅̆1AQcv؇ ':L\sÈ$!'>TWQ^Ί[|jUmÊȊ>ԝ(u1;M^k^%N0U‹2`Ut8 쉡haljˆmntrmcUMB1(&҇xodXLA1z}܅Ȇ.?O`qׇ5AFRau-ITR]Š6F0'PӋb唳):󚳚[t旦[Ԛ vl鐓Xw[\W^gx.:A@Dgˉ(c֊ Vk!Ϗԓϖ VSMȔL DŘP󎢍ˌ ފ؊ŊvhmX1󉘉.Èy`rxT)ЇćwiYMB0z}ۅƆ*7BQdv҇/MVTW[iԈ 1:>;>aމB}wj۔uɓJ|̔i Ȓ͎vŒL;Y}w7ىu#Ȉ{{[% ݇̇{obUJ@-z}مÆ&2=K^n{̇:HMMKTo +..1;Xډ+'ˊ}r܋*݌m{Џ]Б$xєh D id.pU;x{ QшtmruM݇ʇrcXOG<+z}؅".8EVeqLJއ*:BB=DZ{ш'3Qt׊c҉҉"[vꋂDm3:@ߋ+8c)ꊜK ߉ƉЉЉnЈ2%,Xz݊K*tOSjrbB߇ʇ{pke_ZVQME:)z}Յև!%)/9ER_l{· 9IJCES`ilk.*&5Pi~3 +# 0RreljV ';>;AIVwi.(Y#Έbzk>؇‡qiga\XURMD8&z}Յԇ#*5AO\huՇ&42)+@YjpjqЉ!$/Fiʊ?ߊɊeb%P{k쉬R#2$:#‹>C?vCAQ숎Ȉb/ ˇqkgb[WTQKA4z}ԅӇ )3>JVam{ȇ.Pgne`x&;]zJ܊Ԋ>틈4"A_@lZZ8ЉpJ9',CGYfيۊȋW#}nbщY鈎ɈM" އχ{snic]VRNG=0z}ԅ҇(09CMWbp~ч;Xc]T^ !0E`ijxA􊽊?#׋s8+;'ꊭnȉgv^]-F0쉏)|w\8ЇysplcZRNH@7,z}ӅЇ $+3;DNYguևDZVLKdމ$((;HHXyˉ2{ՊԊvьЋ}Y_}r+؊ɊΊXʼnΉىHC }|zG g,:؈bTQI<*Їwtrk^RLGA91(z}҅̇$,4=GSan{ɇ݇ 4QQH@Hf-'%,=hߊ(OȊЊ͋YB0L~Y NJ?鉸ˉȉ:ۈˆ͉*3ՈiL<51,!|vtpfXMGC=6/$z}хʇ &.7AMYerՇ$݉ƉÉtX=%܈ƈ}z[A5+ }N vi`ZTNJFA8-!z}Ѕ† ")0:DLWgx͇܇#/2$3Ojۈ @lϊW|s=܉։aC)ˈ~rv}rbJ1$ ևƇqe[TOKGC=3(z}υ +#+6?GP^l|‡Ӈ%( 4Nbrֈ܈ڈ߈J͊8A"쉺ĉ߉ˉ]*׈ԈֈӈĈxou}xyug]YQ@*  Ƈwk`UMHEA<6,$z}΅ (29AJUcpȇ؇-DRZf~ȈLjĈΈ 9v؉։vs׉ӉÉb$Ԉ̈ˆuggy~}pXF=9. Їwnf[OGB>:4.&z}ͅ +%-3;CLYesˇڇ  +&=KMN_yψ%Rt}vfRK`jH܈Ljˆticfu}vlfW>)χxohaWLC<73-'! +z}ͅ "(.4;DOZfw͇هއ %:A?GZm|҈ڈ$8<:5''CbtseTB&܈Ljyoga_[XUQJ@1 ߇هχ{umd]UK@71-'"z}̅ +#(-3HR`p|LJ̇χЇ̇LJ·އڇ·!3>J`{xlmĈɈ̈Έňy`SJ@<=7& ڇՇӇч̇ȇćunie`\XSMF<2)  +z}˅ #(08?DPakvȇهۇӇۇ #4BVgkdYUav{\F>80+%ՇЇчЇ͇Ƈrf_[YUQLF?6-$ z}Ʌ !*39=GS]iwʇԇׇԇև +1HVWQHEMZfo}{|uY?0+& ܇ۇۇՇЇ·͇ɇ‡zl`YVUQLGA:2)! ܂z}ȅ $+39AJT^juŇ͇͇̇̇чۇ-=B?:9<@K^xzlmoh^Z[]abZK9)#܇؇ԇӇ҇ч·ɇƇƇŇ|zwoe]XTQLGB;3,%܂z}Dž #)19BKT^iu‡ÇćŇ͇ևه؇݇*010-&,B[f\TXXOE><;974." +؇χȇƇȇˇʇƇ}vroke`[VQJD?:3,&  ۂz}Dž !)08@IS^is|ȇχ҇Їчڇ &'  #'*3>?5,))'# ! އԇ̇ćÇŇÇzuolgb^ZTOIC=61*$ ւz}Dž "(.6@KT]fov}ćƇɇʇȇ·ׇڇ݇   ؇чʇÇ~{vpkgb]XRKEA<5.(# ӂz}ƅކ +"(0;DJQZbipx~Çʇˇɇч܇ڇׇчʇć~~}{unic\VQJC>83-'! ҂z}Ņۆ $,5<@HOW^emtx{~̇ۇ܇؇ۇއӇLJʇχ͇·ччӇׇهևч̇Ňzz||xrkd\TNGC>80+&  ߆҂z}Ņ؆݆ '.38@GLS[binps{͇ՇևևڇهӇЇ҇؇ڇՇ͇ćÇƇŇÇ{xwvtojc[SLE?=91)$  ߆тz}ŅՆچچ݆ !&*18>BHOW]adipw}ćɇʇ̇̇Ƈȇ··Ƈ{vsqokgb^YSLE?;82*" ކςz}Ņ҆ֆ׆چކ߆ $+05:?EMSVZ_cfktLJŇ~|xqligea[VSPKE@<83+$ ܆΂z}ąΆӆֆچ܆݆#',26962-&  ۆ͂z}Å̆҆ֆ؆چ݆ $*/5;@CFHHJP[ht~wvxxvslgc^ZVTROKC;50.*& ߆چ΂z}̆ӆՆՆ׆ۆކ߆߆ $).37:<>@@CLXep{~wtrqnkfa\VRONLJE>6/*'$ +߆݆؆˂z}ʆ҆ӆ҆Ԇ؆ۆ܆݆݆܆݆ +#&*.137;;81+%! ߆݆ۆֆɂz}džΆΆΆцԆ׆چچننۆކ +!%(,1678=DLU^gmqtuvxzzz|~zupje`[XUOJHHEA>;61+%! +݆ۆ؆ӆȂz}ņɆɆˆΆ҆Ն׆ֆՆֆچކ +!&+/159>BHPY`eikmopopqrsw{~~|||||{{zvpkgc_XSOKHFEC>:63.(#  ݆چ׆ՆІĂz}†ņdžˆ͆І҆ӆ҆҆Ն؆܆  $',159<@FNTY]abdefgfefjmosxywxzzy{{yvutsqqpqplfa^\XQKGDBBA<72.+&! +߆ۆ׆Ԇ҆͆z}Ćdžʆ̆ΆφφφцԆֆچކ #'+.27:>CHLRVWY\__\Z\^afklmortrpooonljhfeedc`\XVTPJE@==>;60+(# ݆نՆцΆɆz}ÆƆȆʆˆ̆͆φцӆՆ؆܆ކ +  #%'+/259=AFJLOSVVTSSVartiehkkifdccba`^\[ZYWUQNMKGC>:99961,($ ߆܆׆ӆΆˆƆz}†ĆņdžȆɆ̆І҆ӆՆ؆ۆ܆ކ +  "%(+-/37;>AEILLLLMSgl_``_]\ZYYYXWUSSRPOMIFCA>;855431-(%! ߆܆نՆц͆ɆĆz}†ÆÆņņdžˆІ҆ӆՆنچۆۆކ  #%(+.036:?BDDDFM^z|dWVUSQPPPQRQOMMLJHFEB?;87520/.,*'#  ކ܆چ׆ԆІ̆ȆÆz}†ÆĆƆʆ͆ІӆՆ׆؆؆چ݆ !#%'*-059::;=BJWZSNMMLHGHHIJJHGEDA?>=;85210-*))%"  +܆چچ؆Նӆφ̆ɆÆz}ÆÆÆƆʆ͆ц҆ӆՆֆ؆܆ކ߆  !$'+.01248>DHHGGHGDBA@ABBA?=;9875420-,+(%$$  ݆ن؆׆ՆӆцΆˆȆÆz}†Æ†ĆȆˆΆІц҆ӆՆنۆۆ݆߆ + !$')+-27;?AABBA?>=;::99860,./.--+)'$"  ކچ׆ՆԆӆцΆˆȆņz}}††ĆȆˆ̆ΆφІІӆֆ؆نچ݆ !$&).259::;;9998642220-$#''''%"  ݆ۆ؆ՆӆцІΆˆdžņ†z}zĆȆʆ̆͆͆ΆΆІ҆Ն׆چ܆݆݆݆ކ !$(,.1321122110.,,+)(!"""! +ۆ؆ֆӆцІΆˆȆņ†z}uyÆņȆˆ̆ʆʆˆ̆φ҆Ն׆؆چچچۆކ !%'(**))**)**)''%$$! چֆӆ҆І͆ˆɆdžÆz}ru†ĆƆdžƆdžȆʆ͆ц҆҆Ԇֆ׆׆ن܆݆ކ  !"#####$$#"  ݆نֆ҆φ͆ʆɆȆņ†z}sv†ÆĆņƆȆ͆φΆΆцӆԆՆ؆چۆކ  ߆݆݆݆ۆن׆ԆІ͆ˆɆdžņÆz}ux†ÆĆdžˆ͆̆͆φІ҆Ԇֆֆن݆߆ +  +݆چن؆ֆԆӆ҆φ̆ˆȆņ†z}uy†ņȆɆʆ̆ΆφцӆӆԆ׆چ܆ކ    +߆ۆ؆ֆԆ҆цІφ͆̆ʆdžÆz}ux†Æņdžʆ̆ΆφІцԆֆنچچچۆ݆ކ߆߆߆ +  +ކ܆ن׆ԆцІφΆ̆ˆˆɆņz}ux†ÆƆɆˆ͆ΆІ҆Ԇ׆؆׆ֆ׆؆نچۆ݆߆ ݆چ؆ֆԆцΆ͆͆̆ʆɆɆdžĆz}uw†ĆƆȆʆ̆͆ΆІӆՆՆՆՆԆՆֆ؆چۆ݆߆߆چֆԆ҆цφ̆̆ˆɆȆdždžƆÆz}suÆĆƆɆʆʆˆ̆ΆІ҆ӆӆӆԆԆӆԆ׆چ܆ކ߆ކކ݆ۆ׆ӆІΆΆΆ̆ˆɆƆņņƆƆ†z}rs†ĆdžȆȆȆɆʆ̆φІІцццІцӆֆچ݆߆߆݆ۆن׆ӆІΆ̆̆̆ˆɆdžƆÆÆÆz}qs†ÆĆƆdždžȆɆˆˆˆ͆ΆΆΆφ҆Ԇ׆ۆކ߆߆߆݆چ׆Ն҆φ͆ʆʆʆȆƆĆƆ†z}pr†ņņņņƆdždžȆˆˆ̆͆φцӆֆن܆ކކކ߆߆ކކ݆ۆ؆ՆӆцφˆȆȆȆņ†ÆÆƆz}npÆƆ†ÆÆņdžȆȆʆˆ͆ΆцԆ؆چ܆܆܆ކކކކ݆݆݆݆ކކ݆܆ۆنՆ҆цφ͆ɆdžƆņƆz}ln†ÆņƆƆȆɆʆ͆І҆ֆنچچۆ܆݆ކ߆ކ܆ۆۆچچۆۆچۆކ݆ۆچچن؆ՆцΆ͆ˆȆĆƆ†z}kkÆÆÆņƆdžɆ̆φІԆ׆ننچچچۆކ߆ކ݆݆݆܆چنننننننچۆ݆݆݆ۆ؆׆ֆֆ׆Նφ̆̆ʆdž†z}ii†Æņdžʆ͆φІӆՆֆ׆؆نچۆ܆݆ކކ݆݆ކ߆ކ݆ۆچچن؆ֆֆ؆ن؆׆نۆ܆܆ۆۆن׆Ԇ҆҆ӆц͆̆ˆɆƆ†z}ig†ÆņȆˆ͆φІц҆ԆԆՆ؆چچچچچچچۆ܆݆߆ކ߆ކކ߆߆܆ن؆׆ֆԆՆ׆׆ֆֆ؆ۆچچن׆ֆԆ҆φΆΆ͆ˆʆɆȆņz}gf†ÆdžɆʆˆΆІІІцІцԆֆֆ׆ֆֆ׆׆؆ن܆߆݆ۆ܆݆ކ߆߆ކކ݆چنچކ܆؆ֆԆӆӆԆՆՆԆՆ׆؆؆؆ֆ҆цІφ̆ʆɆȆdžƆņÆz}fdÆņdžȆȆʆ̆ΆφΆ͆̆͆φцц҆҆ԆֆֆՆֆن݆݆ن؆نچچۆ܆܆ۆۆۆۆنֆՆ׆ކنֆԆ҆цц҆҆ӆӆԆՆֆֆՆ҆φΆ͆ˆɆȆdžņĆƆz}fdÆĆĆņƆdžɆˆ̆̆̆ʆɆˆ͆͆͆͆ΆцՆՆԆՆ׆׆׆ֆֆֆֆֆֆֆ׆׆׆׆ֆՆՆӆԆچ݆نՆӆ҆ІІцццц҆ӆӆԆԆ҆φ͆ˆɆdžƆƆņÆƆz}fdņ†ĆȆȆɆʆʆʆɆdždžɆʆʆʆɆˆΆІІ҆ԆՆԆӆԆԆԆӆ҆҆҆ԆԆԆ҆҆ӆӆц҆ԆՆ҆ІІφΆІцццц҆҆ццІφ͆ˆȆdžņƆ†z}dbĆȆφˆɆɆɆɆdžƆņƆdždžȆȆɆɆʆˆΆцццц҆ццІІІІццццццІφІІφΆΆ͆̆͆φІφφІцІφ͆ˆˆˆȆņņņ†z}b`ˆц̆ȆdžȆȆƆņĆĆĆņdždždždždžɆ̆͆ΆΆΆ͆͆ΆφφΆφΆΆІӆ҆ІΆΆ͆̆ˆˆˆˆ̆͆Ά͆͆͆Άφφ͆ʆȆȆdžņ††z|`]†Ȇ̆ɆƆņĆņņņĆĆĆĆņņņņdžɆˆˆˆˆʆʆˆ͆φΆΆ͆ˆˆΆццΆ͆͆ˆʆʆʆʆˆ͆Ά͆͆͆͆͆͆͆ˆɆȆƆÆz|_[ÆņƆƆ†ÆƆÆĆĆĆĆĆĆƆȆɆʆʆɆȆȆɆˆ͆͆̆ˆɆɆˆ͆͆ˆˆ̆ˆˆʆʆʆ̆͆̆̆͆Ά͆ˆʆɆȆȆdžƆz|]YƆ†ÆÆĆĆĆĆƆdžȆɆɆɆȆdždžȆȆȆȆȆdždžɆʆɆdždžɆʆˆʆʆ̆͆̆ˆ̆Ά͆ˆʆɆȆdžņĆƆz|[VÆÆĆņƆƆƆƆƆƆdžņĆņƆÆÆÆņdždžņĆĆņƆȆɆʆ̆̆ˆˆ̆͆ˆɆɆȆdžņ†z|ZT†ÆĆĆÆƆ†ÆƆ††ņņ††Æņdžʆˆʆˆ̆̆ʆȆƆƆņƆz|ZT†ÆƆȆȆɆˆ̆ˆɆdžņĆÆz|[T†ņƆƆȆɆɆȆȆƆņƆz}[SÆņņƆƆņņņĆÆÆz}ZPÆĆĆƆ†z}YM}††z}UJz~Ć}z|PEuz~Ć҆Նʆyz|L@pw{}}Ć/,Æzz|I/`filnquxzz|~{qz|>._egilosvwxz}}zoz|=-^cehjmrttuwyz{|~}znz|<,]ceghknpqsttuwy{}~}~~||{{zoz|;+\bdfghjmopppruwxyyxyz{}~|zyxwwpz|;+\aceeegjklllnqtvutuvwx{}~~~}zyxvtsmz|<-]`bdccdegghhikoqpqsuuvz|~~~~|{ywvtrpgz|=,\`bbaaaacdeeffikloqrruy|~}zxxwvtrqnbz|<+\_``_^^_abcdedegjmnnptx{}~~~}}~~}zvututrpok`z|;)Y\]]]\\]_`bcdddehjklnquz|}~~}|}}|}|zvtsrrqpnmi^z|9&UYYZZYZ\]_`acccdfhiklnquxz}}{zyz}~{z{|{{yurqpopomlkg[z|7$SWWVVVWY[]^^`abcdfgiklnorvxxwwwwz{|~~{xwxxxwuronmmmljjieZz|5!PTTTTTUXZ\\\^`abccdfiiiknqrrtvwwwwxz}~~~~~~~~{xvuttsrpnmlkkigghhdXz|1MQRRSSTWY[[[]^_`abbcddegiknoqtvvvuvvy|~}~~~~~~~~~}|}}|{|}|{yvtsrqomlkkjihfdeffbWz|/KPPQRSTWXZ[[\\\]`aaaa`acdgklmpstsrstvx{||}~~~}|{z||zz|||}~~}}}}}}~~~}||~|zxxxwxwvsqpoljihhhggededc_Tz|.JOOPQRSUWYZZZZ[\^```_^^_`cgiijnoooprtuvwy{}|zzzzyywuwz||}}{yxyz{|}~~||}~~}}||{{{{{|}~}|||{zzxwvvutuvuronliggffffedeeda\Pz|0KNNNOQRSUXXXXZ[\]^__^]\\\`deefhjklmorstuwyyvtuuuutsrsuvxyz{zyxxxz~|zz{}}}}{zyyyxxxxxy|~}{{||ywutsstsrstspmkhfffededcdddb`\Oz|2KNMMMOQQSUUUVY[\]^]]][ZZZ\^_`bdefhjknprtwwvspoonoppopqqrtwyzywuux~|yxxyyyxwuuvwwuuvvxz{{{||zwutrqqqpoopomkhfeedbbbabbba`^[Pz|0IMLKMOQQRSRSUWY[^^[XXXXXXXYZ\]_``dghjmqsttspmkihilmllnnnoruutttuw}zwwvvuutrqrtvustuuvwxy{|zwutsrponlkklljhfedca__`aa_^^][XNz|-FJJILOPQRQQRTVWYZ[XUUUUUUUVVXZ]]^acegkoqppomlkheeghhikllmopoopstvx{}yuuvvutrpopprttssttttvwyxvutrponmkihhiigeddc`^]^`a`]\\YVTIz|*CGGHJKMOQPQRSSTTUUTTUUSQRTVVVXZ[\_abeilmllklopkeccdefhijlmmlmnqrrsssrqrtsrqpnnnnoqrssssrrsvwusuurmjjigeeggfeba`_^]]^`_][ZXVSODz|(?DFGHHIKOQRRQQPPPQQRTTRPQTUUVVVWY\^_aeghhhgm|vgccccdfghjkkklmopponnmmnoonlllmlkklnqrrqrrsttsrstoigffcbcddcb_\[[\\\]][XVUTSPK@z|% ;ADFGGFHLPQPPOMMMMMOPPPPQQSTTTTTVXZ\]`abbbck|hbbbbceffffghjkmlkjjjjjkkjhfgiiiiiilnoopqqrqponnnjeccca`aa`_^[XVWZZZ[ZVSRQPOLH=z|" 8>BDEDCFJNNMMMLLLKJJJLMNNNOQRQQSTVXZZZ[\]]]bqzoda_^`bddccccegiigfgggggggedefedegghjkllnoponlkihgd`^_``__^][YVTTUVVVWTQOONMLIE;z| 6;>?@ABEHIIJJKLKJJHGHIIJJKMNOOPQSTVWWWXXXXYZ^a`_^\\_bdcaabcdffeefffedccb`bild`bddeghijlmljihfdca^[[\]__][YWTRSSRRRRROMMMLKIFB6z|49:;<>@CDDEFGHJIHHIHHGFEGIKMNOPPQRSRRTTUUVWWXXYZZZ\_cgeaabdeecbcddcb`^]]\^gjb^_``bdfghjigeccb`_\[YXX[]][XURPQRRPONNNMKJJJHEC?2z|16779;>???ACDFFEEGGGGFEEFHIKMOOOOOONNPQRUWWWVUUVWY\^dmkb_`cdc`__`__][YYXXY]^]]]]]`bdeefeca_^]\ZXYXWVXYXVSQOOQQOMMKJKMKHHGDB@=1z|,33468:;;=?@ACCBCDDDEEDFFGGIKMNNNMMLLLMORUVWUSSTVZ\]bmka]]_a`^]\[ZYXVVUVVWWXZ[\\]_`aaaba`^\[ZXVUVVTTUUSQONOQQNKKJHGHKJGFC@?>;0z|(.023568:;<<>@??@A@BCCDFFFGHJKLLLLMLKJJKMPSUUTTTW\]Z[ab]ZZ\]\[ZYWUTSRRRSTVVVWXY[\]^^^^]]\[YXVTSRRPOQRRPNLNPPOKIHGECCEEDB@>>=8/z|%*-/01357899:;;;<=>@ABDDDFHHIJKKKLMKIHGGHKPSTTSTW[YTTXYWVVXYXVUTSRPPQPOOQTTSTUWYZZ[[[ZXWVVVTQOOOMKKMNNNMLMMLIFDB@>;;=?>=<;::6.z|#(*+,-/0234667667:<=?ABCCEGGGIJIIJJIGFDDEGLPQQPQSSQNOQRQQRSTRQONMNNOOMLLOPQQQQTVVVVVVUTRQPONMLKJIGGHIJKJIGC?:41.,)),02100//0.&z|#%&&&'(*,.011113579:<=>?ACCCDFEEEDCB@>>@CGKLMLLLKJKKKLLLMNNMLKIHIKLKJJKLLLLMOQRQONLKKJIHGEDDB@??><;;;:73-# + + +  z| "#$%')**,.//036:977898898653457;?ABCDDEFGFEEFFGGFFGFCBBDDCCEEEECCDEDA=:730.-,+)'%# +܅υƅ…ą…ƅʅĂz|˅    %+/246:>AA?>=<;:98777642110010/.,*(" +ۅօх̅ȅƅ…~oc\XSLFCBCIRWVTTUX[Xz|w؅υ҅҅҅҅҅Ӆԅօ؅ۅ݅ޅ$+1442.*&" ߅Ѕytuvqib\VPLGA;6/$݄؄Ԅ΄Ƅ˄Єфτ΄Є҄Մтz|iy………ÅąŅŅƅDžʅ̅ͅЅׅ +ׅυȅwfRB81) ބք΄Ƅ~m\MEB>6.)&%'.7;;98:?D@z{YGmnmllnopqsuwz|}}ɅЅօڅۅօυƅtg\PE<4.*&! ل„}naWRRSQG<3,%  σz{Kn>@?>>@BABDFHLNNNPSUVXYYY[\\]_^\[YWUTSSTW\dlu}|rj`SC0 Մʄ{lX?.&؃σ˃˃˃ƃyqeVG8,&$  !)'XTENSION= 'IMAGE ' / Image extension BITPIX = 16 / Bits per pixel NAXIS = 2 / Number of axes NAXIS1 = 132 / Axis length NAXIS2 = 288 / Axis length PCOUNT = 0 / No 'random' parameters GCOUNT = 1 / Only one group BSCALE = 1.000000E0 / REAL = TAPE*BSCALE + BZERO BZERO = 3.276800E4 / EXTNAME = 'SCI ' / Extension name EXTVER = 3 / Extension version INHERIT = F / Inherits global header ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator DATE = '2010-07-02T19:55:23' / Date FITS file was generated IRAF-TLM= '2010-07-02T19:59:26' / Time of last modification CTYPE1 = 'RA---TAN' / R.A. in tangent plane projection CRPIX1 = -66.6459346493144 / Ref pix of axis 1 CRVAL1 = 306.871725219222 / RA at Ref pix in decimal degrees CTYPE2 = 'DEC--TAN' / DEC. in tangent plane projection CRPIX2 = 143.15396209091 / Ref pix of axis 2 CRVAL2 = 37.3733845984976 / DEC at Ref pix in decimal degrees CD1_1 = -2.9402995676385E-4 / WCS matrix element 1 1 CD1_2 = 1.36706049992071E-4 / WCS matrix element 1 2 CD2_1 = 1.37399803289694E-4 / WCS matrix element 2 1 CD2_2 = 2.93132453649777E-4 / WCS matrix element 2 2 RADECSYS= 'FK5 ' / R.A./DEC. coordinate system reference EQUINOX = 2000. / Equinox of coordinate system MJD-OBS = 52162.6261539214 / MJD of start of obseration CCDNAME = 'EEV 9273-20-03' / CCD chip name(s) AMPNAME = 'EEV 9273-20-03, left' / Amplifier name(s) GAIN = 2.3 / Amplifier gain RDNOISE = 3.0 / Readout noise CCDSIZE = '[1:2048,1:4644]' / CCD chip size CCDSUM = '2 2 ' / CCD sum CCDSEC = '[1:2048,1:4608]' / CCD section(s) DETSEC = '[4097:6144,1:4608]' / Detector section(s) DATASEC = '[1:1024,1:2304]' / Data section(s) BIASSEC = '[1025:1056,1:2304]' / Bias section(s) FRMNAME = 'dhsFrame' / Name assigned to a frame of data. FRAMEID = '2 ' / Frame Id DATATYP = ' ' / Type of data frame WCSDIM = 2 LTV1 = 0.4375 LTV2 = 0.4375 LTM1_1 = 0.125 LTM2_2 = 0.125 WAT0_001= 'system=image' WAT1_001= 'wtype=tan axtype=ra' WAT2_001= 'wtype=tan axtype=dec' END q|wvusolkkiededa^[[[YUU\b^UOKD8)ȃiG"ׂ` ā{xuoeWF0݃wR+ւVǁԄބ܄ڄل؄ׄՄф΄̈́̄ʄƄĄÄ„xeM1Ѓ]3قXɁۄׄ؄ل؄քՄ҄Є̈́ʄDŽńĄ}aC$ރb6 +ۂYɁ ߄ۄلׄքԄф΄̄Ʉ„lN. 惽d8 ۂYɁ'*&" +ބ܄لքԄӄτƄrT3烾e8 +ڂY ʁ,.*'%#! +݄ۄلׄԄτȄwW4郿e7 ق[ ˁ./+*('&# ބڄׄՄ҄ʄxX6郾e7 ۂ] ʁ!00-,,+*'#  + ݄ۄڄ؄Ԅ̈́ÄyY7ꃿd7 +܂[ +ʁ!2310/0.)%#  +߄ބ݄ڄՄτĄzZ8ꃾe8 ܂Z +ʁ$4654320,(%#  ߄ۄ؄фĄ{[8郿e8 ۂZ ˁ'7987642.+(&$"  + ބڄ҄Ä{[7ꃿe8 +ڂ[ ˁ(9<::9641-+)'%#! ߄ۄфÄ{[7郾e7 +ڂ\ ˁ*;>>=;8630.-+(&$#"!  لфĄzZ7郾d6 ڂ] ́,=@@?=:84210.,)(&%#! ݄؄фÄzZ7郿d7 +ۂ]́/?CB@><964320/-+)'$"! +߄܄؄Є„yY7g9 ܂_΁2CGEBA?<:86541/.,(%$#  ݄؄ЄĄyY7g: ܂_΁7HKIFEC@><987420-)'&%"  ܄ׄфĄxY8h: ۂ_ρ<:9641-+*)'$"  ބۄׄЄ„xZ8g9 ۂaρ@NQPONLIDA><:852/.,*)&$! +ބڄՄ̈́zZ7g9 ܂a΁?NRTRPRPJDA><:9631.,*&""!  ܄؄ӄ̄yZ7f9 ۂ_́?PUWUTXXOGC@><;9630-+$ "#"!  +߄ۄׄӄ̄yY6郾a2҂W ȁATYZXWXWPIEB@><:7431.&!$&%$"!  +߄܄؄ӄ˄wX3 +݃}L낼EFX\[ZYWUPKHEB?><86651)#'*)&$$#  +߄܄ׄфʄnF胷T$łd)遷K\_]\ZXUQNKHEB@>;9962-)+.-*('%#!  ݄؄ӄ̄m@߃Oe7́N]__]\ZWSPMKHDB@=<:6310021.,+(&$#" ܄Ԅ˄vO#ǃk=܂|Q(ҁR_aaa^\YVSPMKHECA>;755335420.+(&%%" +܄ӄĄfE ̓tI삺Y.Vceeda^\YUQOMKHFD@=;976665420-*)('$  +؄ʄhI'ヽkA삼W)݁Zikieb`^ZVRPOMKIFB??=;;:86531/-+)'$  ܄˄gH'ჿvP(ׂ}LЁ_ppkfcb`]YUSQONKHDBB@??=;86431/,)&$" + ބτcE$݃qL&؂a5 āeuslgeedb]YURQPMJHFEEC@>=;8642/,(%%$! + τa@كiDт]7Ёfutokihhd`[XUSRPMKJJIFB@?=;9641.+('&$! ҄~]<Ӄb<ƂzT/ ҁfuurnkjhda^ZWTSQPONMLHDB@?@?;741-)()'# ф}\9σ\4 傽nI%΁Ł́ссāfuvtqnkgda_[WTSSSSQOMJGDBEMMD;50,)((% ЄyX5ǃ}V-܂e?ʁŁԁ灺euxvsplhec`\YVVWWVTSPLHEDM^`N=4.+)(%  ΄rQ/ +僾rK$ӂZ5ځǁǁׁ&5@GH: +ȁewzxuroljgb]ZZZ[ZXWVRNIFFO^aN<3/.+&̈́mJ'ރg@ȂwO)ԁāȁׁ1Kcvp+ԁgy|zwusqole`_^^][YXWSOKHHKPPF:42.' DŽhE Ճ^5 傼lE ЁÁɁׁ1LgȂӂׂՂʂHށi{}|zxwvtnifdba_][YVSOMJIHGE?82,!Ąb>̓}U+ق_:΁Łāʁׁ0Khق ς_䁦l{~~}|zyvqmkhedb`\YWTQNKJHE@8/#„[7ÃsJ ͂{S. ܁ρʁ΁؁/Kgۂ0GXac^P6p遧o~~}{yvsplhfec_[YWSOKIF?4%ބyV0 +⃺jAnH$ځҁҁځ-Ieڂ3Qm] {끨r~~{uokhgda^]ZTNHA7'܄tO)ۃ^6䂶b<ہׁ܁0Kdۂ4QnÃ؃ڃŃy쁨ryrmjhd`^]XPF9(քmG!ӃT)ՂX1߁*Gkڂ6TqɃ'-)僽$큨szrmjgc_[VM>*Ԅe>ǃuJȂtN* (D`Ȃ7Ut̓!;SenhX>ă%tyqljfaYN?+Є_6僼h>炻hB!&A^}ڂEbyЃ +'A[srM"ƃ$큩uxqmh_SA)̈́X0܃]1ڂ\9&@\z؂8Wwԃ+Gb{ȄɄzP#ƃ#쁩wtkbVC+ DŽ{R(ԃ|P$͂yS2&?Zxׂ7Wwփ/Lf˄Մڄ݄݄ׄĄ}R$ƃ"쁩yqcS@' +„tK ɃrE낿nK. $#ᄻkA샿e8 ߂dD+([vƄۄ  ߄ʄX*Ƀ"|聧fiW;̈́}Oe6ۂpZLIO]r͂4Wyރ @`˄  ߄ʄX+˃$}聨91ƄqG냸X*ӂq_WXcv͂3Ux"Aaӄ #" ̄Y,̓&~遨ㄿi;烸Oʂvieky˂5Xz&Ffք#(*++*'#  ΄[.΃)遨a4փPǂ}vwɂ 0V{ă&Hi„ۄ(.01110/,*'$  Є^/΃(遨>&̃l@ʂʂ /Sw +-MmȄ &.3678754331/-,(# ф_0̓$|聨Ճc7̂˂ -Qv -Pq΄#-5;>>===<:875310/-)%  +҄_0̓!z聨lS/ тς*Ot -QtԄ,7>CEGGEB@??@><96300.+(%! ӄ_0̓!z聨قʂɂق +Mr -Quׄ$4AIMNNNNMKHECBBBA>96320.*'%! +ׄ`0̃ y灨ނւӂւ0Ru .Rvڄ(:ITZ\[XVTRQPNKIGFEEB>:8520.+(%! ބ_0ʃx灨؂1Sw 1Uxބ/APZaeffd`][XUSQONLJHGFC?<9620.-*'# ᄸ^/ʃ"|聨2Rwà 3Z5JZdkopomkjfba_[WUSRQOMKIGDA>;62/.-,)&! ݄_0˃)ꁩKsƒ7[~Ʉ 7L^ltwyyywtpnkhfea]ZXVTTROLJHFDA<730/.-+'# +ل`1̓(ꁩ9`˄ *BSbox~}zwspnkjhd`][XWWUQMJIHEA=9620//,($  +لb3Ѓ#z灩*\ӄ.G]oz|yvsqomjfb`^\ZYVSOLJIEA>;7410/-($! ڄc4҃!w恩΄7Ocu~{xurqnkgdba_[XVTQMKIEB>;8520.+'$" ۄc4҃!v恩3Sk~{ywtrolifeb^[XVTQNLJHD@<9531.*&$"  +ۄc4҃ u偨Wz|{yvspnljhc^[YWTROLKJGB=9632/+(&$!  ۄd5҃ u偨}{ywtrpqokf`^\XUSPMLKGA=:8530,)(%#"!  ݄d5҃t䁨~{xvusstrmhdb_YVTQPOKFB?<:730-+)'%#!  ݄f6Ӄ t䁨|ywvtssqligfa[WUSQOKFDA=:730.,*)&$" +ބf7Ӄ!u䁨~{yxusroljhfb]YVTQNJGEA>:730/-+*'%" +ބh9Ӄ!u䁨~{zwuspmigec_\XTPMJGDA=:731/-+)'%! j:ԃ!t䁨~|{xuqmifeca^ZUPMKGCA>;9520.,)(%! +„k:փ sと}yvrnjfdca^[WROKGCB?><941/,)(%! l; ׃ rと}yurplhdb`][YUQLGECA@>:630-*'$"  +m= ׃rと}yusqnkgda][ZVRMIGDB@=:751.*&$"! +m? ׃rと…~ywusqnieb_\ZUPMKIFC@=:853/+'$!  „n?؃q⁨DžÅ|{zwtpkheb]YTPNMKHEB>;8651+(%"  +ńp@كpၨɅŅ…}}{wqmjgb^ZURQOMJGD@=9761,)'$  ȄrAڃoၧȅƅÅ|{xsokfa][WTROMKIFB>;972-+'#  ʄsBۃoDžDžą…|zwtplfa_\YUROLKJHD@<:61.+&"! ̄uD܃oʅ˅Dž…{zvsplfb`_[WSPMLKIF@<951.*&$! ΄uDۃo̅υʅÅ~|yvsokfca`\XUROMKIE@<963.*'%" ΄vEۃnͅυ˅ƅÅ}zwtpmieca]YWTQNKHC@=;83.+)&"  ΄wG܃nͅͅʅȅƅą}{yvrpmifc`]ZUQNJFB?><83/-+&$" +΄yH݃n̅̅ͅ˅ɅƅÅ~|{wtrnkhfc`\WSNIEA?>;732/+'%#  τyGރn̅ЅЅ΅˅ɅŅ~{xusoljgeb\WSOIFCA>;7430,(%#"  ЄyGރm҅҅ͅЅ΅˅ȅŅą…~|xvspnkheb]WRNKIFC?<8530,)'&$!҄yH߃mЅӅӅ҅хЅͅ˅ȅŅÅÅ}yvsqnlhe`\WROMJGD@=952/-+*(%! +ӄzI߃l߁ą҅҅҅ӅӅԅ҅΅Ʌƅą…zusqnlhaZYXTPMJGDA>:61..,*(%! ӄzI߃k߁ŅӅԅӅԅօׅԅ΅ɅƅŅ…~zvtqokf[RUWTPKHGEA>;71..-*'$"҄zIރkށŅօ؅օօ؅مՅ΅ʅȅŅÅzwurokdWOTURNJHGD@=;9520-*'%! ҄zH߃jށƅم܅څڅڅمՅυ̅ʅDžą|yurokbVSVTPMKIGC@>=;851.*(%! ӄ{H߃jށͅ߅݅܅܅مօх΅˅DžŅ{wsojb[ZYUQNLKHDBA@>:51.,)%! Ԅ|I߃jށԅޅ݅ۅمׅԅЅ̅ȅƅÅ}xsojfb^[WTQNLHECBA?;620-)&" Մ|I߃kށم݅ۅمׅՅӅυ˅ȅŅ~xrnjhd`]YVTRNJGFCA?;741-)&" Մ}Jjށ܅݅م؅օӅυ˅ɅƅÅ|vqmjgca^[XXWSOLHDA?<840,)%!քLボj݁݅߅ۅمօхͅ˅ʅȅą~yuqmifdb_\ZZZWSOIEB@<840,)$ քM䃝i݁ޅ߅܅څՅх΅̅ʅɅƅÅ~yvsnheca_][ZZXTPKGDB>951-($! ׄN僝h݁߅߅܅مՅ҅Ѕ΅̅Ʌȅą}{xsnidb`_^\ZZXUPLHEC?:62.*&$  ؄P僞h܁߅߅ۅ؅Յ҅хυͅʅDžÅ~{wsokeb`_`_^\YUQMIEB?:62/+(&! +لP僝g܁݅ۅ؅օԅхυͅʅƅ…|xtplgdbabca^ZVROKGC>:62.+)&" ڄQ䃜f܁߅݅مօԅ҅υ̅ɅŅ…~yuplhfdcccb_[WSPNJD>:51.+(%! ۄR䃛eہ߅ۅօԅ҅ЅͅɅƅązupmjgfdccb_[XTPNKF@;61-+'$! +܄Q䃜dہ܅ׅԅӅхͅʅȅŅÅ~ytqomigecba_[WTPNLGB=82.+'$! ܄Q僝dځۅׅՅԅ҅Ѕͅ˅ȅą…|xtssqmiebaa^ZWTRPMIC?940,(%! ۄP惞eہڅׅׅօԅԅ҅υ̅Džą…}ywuusokgdba^[WUTRPJD?:62.*%! ݄Q烞eہۅممׅׅօԅ҅΅ʅDžƅÅ|yxvsoljhea_\XVTSPIC>;841-'! ߄S烝dځޅ܅څمم؅ԅ҅Ѕͅ˅ȅÅ}zvsolkieb_\YVSQMFA>;842.(# T 胞cځ݅ۅۅ܅څօ҅хυͅɅÅ|wspmkhdba^[WSOJEA?;842.*&# U!郟cځ݅܅݅ޅ܅؅Յ҅ЅͅȅÅ}yurnjgedca]XSOJEB>;741/+'$ ᄶT!ꃠcځ߅݅܅څօӅхͅɅƅ…~{xuplheedb]XTPKFB>:8530,'# ㄷU!ꃟcځ܅ۅڅօӅхЅ΅ʅą~{ywtpkgdc`[VSQLGC?<:852-'# 愻V!郟cځ߅݅܅څօӅ҅҅Ӆ΅Dž~zwwwunhdb^ZVSPKFC@?=:62.*&"  +넿Y"ꃠcځޅ݅ۅ؅ԅ҅ӅՅхɅ~zwuutmgca^ZVSOIECA?=:62.,($""" „[$샢cځ߅ޅ݅مԅӅӅӅυʅÅ{wsqnieb`]YVTOHDCA>;:61.+'%##% ]&cځޅ݅څօՅӅυ˅Ʌƅ{vrolhdb_\XWVPHDC@<9852/*&%$%(! ]&cځ߅݅ۅم؅ׅӅ΅ʅʅȅązurokifc_[XVTOIFD@<9641.+(&$&(^(cځ ޅۅ؅؅ׅօ҅υ̅ͅɅŅ}yvrnkjhc_[WUSOKHEA><83/-+)'%#!a-􃦃bځ +߅݅ۅم؅ׅԅххЅͅɅŅ}yvrmjigc`]YWUQMJEA@>:4/-*(%"kC bځ ߅ޅ܅څ؅ԅ҅҅Ѕͅʅƅ…~zwsnkifdb_][XSNJEA?=:50-*&#„\Ճ%cځ + +߅܅مՅӅ҅Ѕ΅˅DžÅ{xtpligec`_]YSOJEA=;:50-+'$Ą-dځ ޅڅׅԅӅ҅Ѕ̅ȅŅ…|xspmjheb`_]YTPMGA=:940.,)%!Ą*dځ ߅ۅ؅օԅӅх΅˅ȅą|wsqnlifa_^]ZVTRLD>:730.,)&# Ƅ|c.ƃ%cځ +܅؅ՅՅԅӅЅ̅Ʌƅ…}xtromkgaZ[][XVTNE?:620-*(&$  ɄwK%cځ ޅمօՅՅԅх̅ȅƅ…}yvspmlhaXY\[XVRKE@<731.*)'$! ΄~QÃ'cځ +߅܅څׅՅՅӅυ˅ȅą|yvspnlje^\\YWVQKEA>9420-+(&# Մ^&ʃ)cف  +߅݅مօՅԅх΅ʅą}yvromllid_\XWVSNHC?;852/-*(%"܄l2у,cف  +݅مօՅԅ҅΅ȅ…~{wsommmlf`\YWWTPJEB?;731.+)&$ńw<؃.cف  ޅمׅՅԅх̅ƅ|zuqponmhb]ZXWTOIEDA=731/,*(%΄Eރ/cف  ܅مօԅх˅ƅ}{xtrqnkhc]YVURLGEDA<8420.+(& ׄM0b؁   ߅ۅׅӅЅ˅ƅ…~|yurqmjgb\XTQNJFEB?<9631/,(&  V1a؁##߅ۅׅ҅υ˅ƅ…{wsqoljfa[VROKHFD@><97630,(% ̄_3bف&&  ߅ۅׅӅυ˅ƅzurpnkhd_ZUROKIGC@=;::840-)&!քh5d݁()# ߅܅؅ԅх̅ȅ…{vrpnieb^ZVRONKHC@><<<951.+(" ބp7i䁲)*%  ݅څօӅ΅Ʌą|wsplfb`^ZVTRPMIFCA@?=;741/.)"儷s9n큹)+'"ޅ܅؅Ӆ΅ɅDžŅ~{wrmheb`^ZXVRNMLJHFCA?><99;=92,$t;t)-*% ܅ׅ҅ͅʅȅŅ}zvojfdb`][XUOLLLMLHECBCEGO\bXI;72*DŽv:Ճ+iꁺ&:;4.(! +܅؅Յ҅΅Ʌƅ|yurnkjhfc_\ZYWTQNMKHECA=951.+$ لh/̓&d⁳':;4.'  ۅօӅЅͅʅƅ|xurnkjhfc_^\ZXVSQOKGDA><93.,)"ф[#ă"aށ(<<4-& ߅مՅх΅̅Ʌą|xvsolkigea_]ZYWTQOKGDB?=82.+' ̄|M_ہ)=<4.' + ݅؅ӅЅͅʅƅ{yvsqnljheb_\ZYVROMJGEB?;5/-)$˄p?]ځ*<;50( +܅؅ԅхͅȅą|yvsqomkheb`]YVTQNKHED@<72.+'!„c0\ف*:94/(  +܅؅օ҅ͅȅÅ}{vspomkheb`]YVSPLHFDB>952/+& ۄW#냝Z؁,::5.' +مօׅօ҅̅DžÅ{wromkjgeb`]ZVSOKGDC@>952/+&!ׄ~Mვ Y؁.=<70' +݅Ѕ҅ͅӅЅ˅ƅ…{wsomjgfeda]ZVRNKGC@?=:51.*'#ՄxFك Xׁ)<<70& +ޅӅυхх΅Ʌƅ}yurnjgggeb]YTPMKFA>=<951-*'#քvA ҃ Xׁ%894-$ ߅مօӅυ˅ȅą}xuqmkigda\WSOMJD@><:751-*&"ׄyAσY؁'75/'! +ۅ؅ׅ҅̅ȅƅ…~zwtqnjfb_ZVSPMIDA?<9741.+%!؄|B΃Yف,94*$߅ׅׅօхɅŅÅ~{xvtpjfb^ZWTQNIEB?<96320+& ׄx?̓Y؁,83*#݅օՅԅͅƅ…~{xvtpkgc_[WTRNIEA>;85320+& ׄt<΃Xׁ,84+$ ߅؅ӅυɅÅ|xvtpkgc`\XURLHEB>;8642/*%!ׄt<ЃWׁ,83,% م҅̅Dž}yvrokfa^[WUQIDDB?<9752-)% Մu>ЃWց,73,% ޅ؅Ӆͅƅ~|xtpmid`^[XVPEACB@=:752-($Ԅu?ЃWց*51*# ܅؅ӅͅDž~|zwsojea^\[YVQHDDB@?<962-&!Մv?ЃWց%2.'! ߅ۅׅ҅̅ƅ~|{wtqnhb_]\[YVQKHEBA@>:51+%!Մu?уXց-+$ +߅څօхʅą~{zyuqolfa_^\[ZWPKIFCA@=950*%!Մt@ ҃ Xׁ*(" ݅څׅӅʅÅ~{yxspnjda`_\ZYVQMJGC?=:850+%!ԄvB ӃXׁ('  ߅څׅօӅʅÅ|zvsrojeba_\YWURPLHC?;8752,% ԄwB ԃWց)' ܅ׅӅ҅υʅą}zvturlgdb`]ZWUSPLGC?;8543,#ՄvA ԃWց)' ޅمՅхυ˅Dž…|ywuwtnida_]ZWVRNJFB?;8431*$ Ԅu@ ԃWՁ&$܅؅ԅЅͅɅŅ{xuttrlgc`^[YXUPLIFB>:7420+&!҄t@ ԃ WՁ! ܅ׅӅЅ˅Džązvsrpnjfb_\ZZXTOLJGC>:7420,'! ӄu? +ӃVՁ ߅ۅׅӅυʅDžązusqnkheb`][[XSMJIGC>:7420+&! ӄu? +ӃVՁ + ޅۅׅӅͅȅŅÅyvtqmhedca_^[WRNJHGC>:631.)# Ԅv@ ӃUՁ݅مׅхʅŅÅzvsqnifedb`_[VROLHEA=9410-'" Մu? ԃTԁ ߅ۅׅԅЅɅą…}wsqpmigda`^[VSPLGC?<720/-(" Ԅt? ԃSԁ݅مօ҅ͅȅą…ysqpnkgb``^YUSQLF@=:51/.-)# Ԅt@ ԃRӁޅۅ؅ՅЅ˅DžÅysqomieb``^XSQPKE?;730.-,)$ ԄvB ԃPҁ܅څׅӅ΅Ʌƅąxsqolhfc`^\WQNMID?9531.-,)#ԄxC ӃNс߅ۅ؅Յυ˅ȅDžƅ…|wtrpmjgc_\ZVQNKGC?9542/-+($" ԄxB ҃LЁޅمօх̅ɅȅȅDžŅ…{xuronkfa][YVSQMGB>9752/,+)&# ҄u@ ҃Jρ܅؅Յх΅̅ɅDžDžŅÅ|ytpnmjd_\ZWTSRME@<9752/-,+($ τs? +ԃH́߅ۅׅԅӅ҅υ˅ȅƅą}xspnkgb^[XTRRPJC>;97410..,)&!τs@كF́݅مօԅԅԅυ˅ɅƅÅ{uromiea]ZURPPMGA><:7410/.,(%"Єs@܃Eˁ݅؅օօօӅ΅ʅɅŅ…~xsqnkhd_\XTQONKFA><;9520/.*'$ ҄t? ׃Dʁ܅؅ׅׅօ҅΅˅ȅą{vspligb^[XTPNKJFA><;962/.-*%"҄u@ ԃCʁޅ܅م؅ׅԅЅυ̅ȅÅ{wrnjhe`][WSOLIGD@=<;:62/.,)$ ҄u@ ӃBɁׅۅمׅՅ҅Ѕυ̅Dž…|wrmifb^\YVRNJHEB>==<9630.+'" ӄu@ ӃBɁ΅ۅ؅ׅՅ҅Ѕͅʅƅ~{vplifb^ZWTRNKHFB>==<8531-)%! ҄uA ҃BɁ҅܅م؅Յ҅υ˅ȅŅ|xtpmjfb^YVUROLIGC@=<:6320,)$  фuA +уBɁԅ߅܅څ؅Յх̅ȅDžązurolieb^ZWURNKHFEA><9520/,)%  +Єv@ ЃBʁ҅ޅۅمօ҅ͅɅƅąyspmjgda]ZXUQNJFEEB?=;730.+)%  фu?уCʁ҅ޅۅׅӅυʅƅą~xrnkheb_\ZWTPLGDDCA>=;851,)'#фt>уCʁ҅߅څՅхͅȅŅ…}xrmjgda^\ZWTOJEBCB?==;851+('#Єs>҃Dˁ҅ޅڅԅЅ̅Džą{wsnjgc`_\ZWSOKHFEDA@><840-*'$τs>҃Eˁԅ߅ۅ؅Ӆυ˅DžÅ}wurmifc`_]ZVROMJHEDCB@=720-*(% ΄r> ӃF́Ӆ݅مօхͅʅƅ…}ytrokhfc`^\YVSPNKGECBA@<61.+('$̈́q> +ԃǴӅޅڅ؅ՅЅ˅ɅŅ{wrpmihfc`^[YVTRNJGECB@>941.*'%"̄r> +ԃǴӅ߅ۅ؅ՅӅυ˅ȅą~{vrnjhhfb_][XWUSNJGFCA?;731.*&# ̈́s> +ԃǴ҅݅څׅՅӅ΅ʅŅ|xurnjhgeb_]ZYXURMIHFC?=:52/,)&" ΄s? +ԃǴӅ݅مօӅυ˅ƅ…{usqnkigeb`][ZXTQMJHEA><951-+)&! ΄s@ ԃF́Յޅ؅Ӆ΅ʅDžąytrolihgdba^\ZWTQMJGC@>;961-+)%  τs@ ԃǴ҅ޅׅЅɅƅą…~xtqmjgfedcb`]YVTPMIFC@=:840-+'" τs@ +ӃH́̅څӅͅDžŅÅ{wuqmigecccc`][XSOKIGD@<:73/-+&! τuA ԃH́ƅׅׅӅ΅ʅDžą{xwtpmjhecbbb`][XSNJHGD@=;72/-*&  ΄uB ԃÍÅ҅ӅЅ̅ʅDžÅ}ywuroljhfeca`^[ZWSOKIHGB<852/,)% ̈́uB ՃÍЅӅЅ΅̅ȅÅ~{xvtspmjigfda_\YWURPLIJJC3+/0.+(# ̈́tA փÍυӅ҅Ѕͅȅą}zxvutrnljhfc`^[XVTQNKHIIA-").,)%#!̈́tA փJ́υ҅хυ˅ƅ…~zwvvuromjgeb`][WVTPLIGEC>3+,-*&$$" ΄s@ փJ΁…Ѕхυ̅DžÅ~zwuutromjgec_\YVTSOKHEA=:61/,)'%%# +τs@ փK΁΅υ˅ȅą|zwutsrqmifeb^[XURQMJFB>;840.,)'&%# +΄s@ ՃK΁̅̅ȅŅÅ~{zwuttsqmifda^[XTQNJGC@=962/.,)'%%# +̄r? ԃL΁ȅȅŅą…}zzxvutspmieba_[XTQMHEB@<741/.+(%$$# ʄr? ӃLρŅƅŅ|{zxwuromid``^YVTPLGDB?;74420+'%##! ʄr@ ԃLρ…ƅą~}|yupmkhda^[WURMJGEC@;66640,(&#! Ȅq? ՃMρÅƅą~~~ztnjhfda]ZWUPKIFDCA=9751.+)'# +DŽp> ՃMρąƅÅ}}~}zslhgec_[YXTOJHEDDDA=94/,*)(# DŽp= +ՃNЁƅƅ…~|zxqkhgeb^ZWWTNIGECCDB>93/,*)(" DŽq> փNρŅŅ~{xsnjhfd`]YVTQLHFEB@@?<62/-)'$ Ƅq? ׃Nρ~zuokigdb^ZXURNJGECA??>:63/+(%! ńp? ؃M΁|xsnjhfc`\XVTQMIGEB@?><974/*&#  ńo? ؃M΁~|zuqnkgda^ZVUROLIGEA?><:853.)%" ńp?كNρ}{wqnljfb_]XUROLJIGEB@>;9741.*&#  +ƄqAڃOρ}zxtoljgc`^[WROMJHHHEB?><9632/+'# ƄrCۃOρ~zxvrnkiea_]ZVRNLJHGGDA?><8422/+'# DŽsCۃNρ}ywuqnkhda_\YVRNMLJIFCA?=:6311.*'# ɄtDۃNρ~{xusomjfda_\YUQONLKIFCA>;7420/,)&" ɄuDكPЁ|zxvspmkheca_]YTQPNKJHEC@<8530.-+(%" ȄtD؃Pс{xvtqnljheca`]YTRQMJIGDB>;841.,+*'$! ȄtB؃Pρ|xusqnkigfcb`\WUSOKJIFDA=:841-+*(&"  DŽsBك P΁|xvtqmjihfdb_[XVRLIIHEC@=9620.+*(%! DŽrBك P΁~zxwtpmjihfc`^\YVQLIGFEC@<830/-+*(%  +DŽrBكPρ|yxwtqnkihea^_^YTPLIGFED@;620.+*)'$ ƄrBڃPρ~{yxvtspligd`^^\WRNKHFEDB?9410.+*)&#  ńrBۃQЁ|{yvtsrnjihea^\YUPLIGFCA@=8410.,)(%#! +ńqAۃ Rс}zyxvtqokihgeb_\YVPKHGFB?=;742/-,)'%$! ƄqAۃ +Sс~zxwuusqmigfeca^]\XQKHFD@=;9752.,+)'&# ńpA܃ Sс~{xwtsrqoljhfda^]\ZWPKHEB><:8762.,)''%  ÄoA܃ RЁ|zxvsqpnmlkigda][ZWSOKGC@=;98752/+(('$ ބpAۃ RЁ~}{yvsqpnmmligea]ZWTPMJGC?<:98531.+*)($  ݄p@ۃ RЁ|xtsrponmkifc`\YVSOKIFB>;:96420.-,*'# ܄p? ڃ SЁ|xusqqqpmkida^ZWURMJGC@=;:85430...)%# ܄o> ؃ TЁ~}{yusqqqoljgb^[XVTQMIFB?=;976430/.+&$# ܄m> ك TЁ}{zywutsssrnkiea\YWUSPLIEB@>:76531/.,($#! ܄k= ڃ TЁ~{xvusrstvuqmjfc_ZWVTQMJGDBA>:76420.,)&$" ܄l= ۃ TЁ|yvtsrrrtuuqkgc_\XUTROKGECA@>;97530-*'%#! ڄn?ڃ Uс}zxvusrqqqrspid`\XUSROLIFCA?><:88630,(&#! +؄n@ڃ Uс{yxwvtspoonnlgb^YVTROLJHEB?>=;98763/*'$" ׄm@ڃ Uс}zxwwvtqnlmmkgda\WTSRNLJFCA@>=<:8652/*&#  քnB܃ Vс~{xwwvurokjkkhda^ZVTSPMKIFB@?>><:8753/*$  ׄvL߃ Vс|xvwwvsqmihhgeb_[WTSQNLJHFC@>>=<98752-(" لĄ^&⃐ Uс~{wwwwurolhgfdb`\YWSQOMJHFEDA?=<;87530+&" ߄Єr,߃RЁ{~zxxxvtqmkhfec`]ZXWTQNKIFEEDB@=<:862/+'$! ܄ʄu"ЃyMρ{|zyywuroligedb^[XVVTQMKHFEDCB@=:974/+&#"  ՄRcG΁||{yxvspnlhdba_]ZWUTRPLJGEDCBA>;9751-($"  ԄZփKڂÁ{~}zxvroljfa^]\[YVTRPMKIFEDCA?<9753/+'$! ҄^(탰u2ɂ9ˁx~~|ywupljgc^\[ZYWUTQNLJGFDCB@=;9631-)&$! ӄl7ȃX0ȁw}|zxurnjgc`][ZYWVUSQMKIGEDCA><:852/+&#"  +҄zK߃r=%ā{~|zyvrolhea^\ZXWVUSROLKJGEDC@=:8530-(%"  фY)􃾃T 傋~}~~}zxwtoljfb][YWWVUTRPMKJIFCBA>;8520-)&$" +ф^5҃g4ǂn|}|||{xvtpljhd_XVWWWTSRPNMKJIEBA?<851.,)&%$! Ԅ_9xE₥J恬r~}{zzyvtqmigeb^YWWWWSPOMLMKIGDA@>:62/+)'%$" քa;郺S!򂿂{ʁk}}{zyyxtpnjfec`^[XWWVSPMLLLKHEB@><73/-*('&#  + ؄c>_-͂D쁳h{}{zywvuqmjfcb`][YWVUTSQNLKJIFB?><950-*(''&%# ؄e@ăf6ׂXÁgz}{yxurrpnkgc`_]ZXXWTSRRQONLJGD?<;852.*)'%%&)'  ؄„gBŃj; ݂`΁gx{ywuqmllkhda_\YXWXVTSRPNNMKIGC=:752/,)(&$"%*)  + ل„hCŃl=߂dҁfvywurokihfdb`]ZWVWXWVTRPNMKHGFB=8420-*('%#"#&# +݄ńiDƃk>eӁdtwurpnkifca`_\YVUVXXXVRONKGECA?<730/,)'$""!  +߄ɄjEŃk>ႫeӁ`rsponljgda_][YXWUUUWYVPMKHEC@<==710.*'$#"! ߄ɄlGǃk>ႪdӁ^oomlkifeca^ZXWWVSRQSUSOLIECA>;;:4/-+'$##! +ʄnI ȃl>₪bҁ]mmkjgdddb_\YVUTSQQPOONLIFC@=;8740-*(%#"" ߄˄nJ!Ƀn@オbс]klkhdbbb`]ZWUTQOPQOMLJHFCA=9642/-*'&%#! +ɄmH!ʃnA䂫bЁYhkjfb``^\ZWVTRONOPMJHFEB@>:631.,)'&&%" ʄnH ˃oA䂬cсTeiiea_][YWUTQNLLLKIGEB@=<:730.+)'&%%#  +̄nH!ʃpBカdсRbfec`\ZYVSRQNKJIGFEEE@;98731.,)'%%$#! ΄oI!̓oA䂫cҁQ_baa^[YWSPOOLKJGDCBDF?96530.-*'%##"  +΄nI"̓rC傭dсM\^^]\ZXUQONMLKJFC@>?@;752/-,+(&$"!  τpJ"̓sE傭eҁHY\[ZYXVSPMLKKJHEB>;97642/-+)(&$"! ЄrM$΃rD悮eӁFWYXWVUTQNLJJIGED@<96321/-+('&%#! +фrM&σsE悭eӁDTWVTSRPMLKJIGDCB=9641/.,*(%$$$#  ҄sM%σtF悮dҁ@PTSRQOMJIHHGDA@>9641/.,)'$   +ӄtN%΃tF炭dс=<:641/.,)'%" +ՄsM%σtF炭dс6EHIJIFB>;:85444320-,+)'%#  +ԄtM%ЃuG肮cЁ /35640)#!&))'&$"! +΄rM%уvG邯dЁބ߄ ӄjF ΃vH邯eЁ{yuv{Ԅ܄؄Ԅ҄τʄ„tW6ǃsFꂰeЁ-@Yu„nV<܃kA邱eЁ`mqokc[X\bkwރ4IZht}yrhZJ5ƃ]7傰fЁȂӂւԂ΂ƂĂׂ̂ +!=]ԃ(9FQWYWUSQPNI?3$˃pO+߂fЁEMOMG?9:ALYi{ʂ ,KgЃ%(*('%#σcC"߂sӁ \ No newline at end of file diff --git a/scikits/image/data/simple.fits b/scikits/image/data/simple.fits new file mode 100644 index 0000000000000000000000000000000000000000..c54dcb1f71a6b31a1e9d965051443909b74c25bd GIT binary patch literal 138240 zcmeFaceq^TdG5VN*zt)^n)AiUIdPJ+lh|>iNJ4ZlsDeNRq%C?ALJd?v5<>LeJEC_G z9R!$ey0J|)c8ptW$K5IJiNkTeH0Qg2>$e})Y$0qXSFY>(CoJ~NXlBitwch(F_fy_? z&73@S*C`VR2DTem&wqmhn+!}^xMKOh@)fg}96Wp3!Ry&)owvqK9^7^Ep6fr}sLgY{ zaUAckxyuF)TexcOqIJIg#(z#5yXWMYd3(L*H#*0gcEsX?apkAUz@!$BMRNxh z&pv8k&eA0-W-nYaFn868xl5KWT)Jev7e9T*ro@O9_B!Ib{=o`505u_THocE`;G@_zr=sz$4}UG^7v`f`;Pbew>I0f zJ)R;kcFI8Z9mvjuvt~>jm^N$bxQR1f-}UtulLseH1>p4@FIU_iZ^p#kW=)(qIB~+j z?DOIq@I@b&!-Cybr3%ZB47kDa;U#KG|!*ZHjb zjDBo1eca3(Z{%j1j2yGcmRs~4t?->VFk!+bQ>SjS*Iw&>f8Br280$MZ-l$Et9=+Ap zqsK&Ax#oJlpGam+6xRaxG(H+v*h7d+fbq`}i;CH+uBeqsDAKlH&}!-n}PI zn>@JFz=X*&2ggnu-&!NTv%be!u?75>;|*l+w#fhM9&gsP$-7USF*71FBwK!C{f{%d z?-w0# zaCdfJpFC$^Jl^2$b&BhZ)%yN9e(J6pO&R+JJlxgejh{Ys*ReAu&&)Nl=lYJ)y*HS% z+bEa34!-@zW0zkvY_A^>gSJr^+s`@P>-+vPFUC#YXwsCiJFdqb7al8KF<3%rW1rw+IqGRzZ<7fQW9&f@fzwmfCj$itEC;j5{n>1tcv>hk@ zG8dcc&GBZAA3J5@8}WVJJ`Wv_-E!c8$YKl*8o0jy24_s3%KNdWV-H(45R(!!IBL|u z;If5_2Nt%Ja(Orv^ejAjpG~*kY}CjN+sFMnBi?u5af3B}!;d$8<=kb<=geLbn!9Y- z&;ggPShn!6_OP49acwqc%z6$yY0~&F{CW}}f~Hd=?hyy@58^n=%rx8C!Mwf_(P_xkZhZpry=HG0_P{>P8E$M~5O zkmxt+M@}sM@<00Dl)>HOwbysNksFT&;Q!!%IbM4fZ#=*M?D3|KI$-_r@5qh+Q|C8z z^Z_};H=f_Xf9iO9j2#?#z{ruKVqcaYwS2|g#oHaW;Hc$^_!bQV`j;NX9B=6Iw3E06a}U&w#&_}?S&-y`tfBk=$82(-z`|H5-<{;b{J zJM_%I-jnUmntJEGdByg=p?Cc3{rNXred_1$Yu?{{u=zl{y+4odZ?_LL8)n-m+a}HE zY@_nLasJpw=3kpOn>Sl#+alZM&1Ts)&(^)(qS>n1y4|*DKV=WLG5O;=?8&xeyWcj= zw#^RBxMur2Z-~%*eKDGqsuA?AT0b#^#*1Z_kfk@v1j%$9-hRdy^K44TBTg`}=TL5=tJ#zwkvK zx~;?BS4zTrnm^6az#SO?dmj1M2ZtnJlZ^J}`5Tq<0D1nz-_dd`M4nL-CC86FxE8Ul zM?R}j0kv&gRK^CA9rCI$EGb}5u-Gckbp#~`so5rDHnv5Qe{J1L0eFCexbX;T;NET9 zx1+7B)!Fbn{hJcSagmySDd2Z%PvxBB*QTU!79~#{x$xgP_)l$iYi2dO=W%KzVq&fw z{5uKYS~rW1sd?(&M!`cpeSfR#U=H3L!v3p!$36VKLm#aNXYIgyns?_Iywhnuwm?b4 z0&J3V*(M{;xnKuIwa=U8{2T$UzpzmjM9I-)1k?c^9c%V;d>EAv1&i%-^<(oWhGJj) zxVOKiH2CNd{Qa$3ggfm6ZH~%uJ+kO%0h9o|LNAyu-#cxt`4YDb9A9A8jES{O4%l?|Lg zg=e+jWBFZm9Hc!)zZ9@X?YrgcIEj|>G&bWpp&j#y?XQTwzoOjK_X~e5OP#Oymj%#& zfKDfQyb}N6_J;Fi!*Xt#KSpEgoWsP{15AtN&&)PGvanM#Dd&hk8aL#XxB|Y9r`t68 zrUwB3@&KT(PcQt@UojT_vJJdf7NMRCeV*UnB8*QFdv#y*u>x2C{XTlHN74r%4|;iR zYmF-CDgpH=7$a>^aV*g&aqj0{ih9 zZ!~V_JYr*!FgESJDEfVTOZgYx>GZ$QE#Fc2`#x9V<}d4I+{bZ%t!M%8RpWCeJICrz z3-;6V%$W_x-|^}2$L6a4urVb7>^tl3t98-uaPRCwzb_E|4+Py&Iqq=lKOBF47e9dK z_iSVW`T$tXj_vs3z10W(L1#U=8Z9*VmnUVcc*U^R} zWyOp5g{^ap9a;-8CCKwPZJ7S+0r)OQ-706I#qH34BmU}l*@ChL;M^}EB^!l*XYJAd zjiTq+kJ|M8GQiUQ;{95qVgWG2Z8OfhM2;uMitW&zx2SNnA}y|12}V(^=M^k`Y@Rzj zJA3Y!GZKtDQiB&S@5tEKxR*Z>4{^mBY?q@?&L~gGBl5|H1jw_nC;SGFKFeH{pAO8c< z&KLZA{6$~>rw@<`g8RDoV`sn~+mE)pfqT7z1!w`b%bBwozZG)m-?$dtPW#1E8@zqw zWpZpEn6in8-b1!j`583T{Bw5Z5*MH>o63bl80^ctx4IRXbT&F6~NLm&fud| zDk&VXqOP(4`XA26^AP2_I{(EZ5?a#8YdwIRhEF^~?W5Ei1;PTb0PCKaN8y`tK<&R) zc>VT𝔢T0YJOZ@5jI2>zDS+{rMN%Pv6h8uSyd11u`3bwd3Hq!@C*%9rLKSEnlEb zK&M4#aHGN11 z$^z(vlmILRXgS_Suf;Mt8LO)+nN$YLmUYrl_-kc2F25q`*hl=s*jx=`vvbaB$}s$u zP4$6q!U17`#Q-*rPHq@v^#35Lm6rkdSU3+yTt5K=EX`NS^y#sez&Ya z<$y*;|INj0(c-4wZIkQMQWMYXaZC0YJN=_$uw$-W>32WTQwG#N(09ZpW~BCP<}~|9 zVX_Tr-S@d=ntXtjTOcs6>*NbejQK(0mYx%@&V^(9 z+SpRbsl5vQF8a{AA?w(o9dqXD24jO2n3|2Vnlz*b*i6okBW#=V*}Uc4XdO6|4=AJy zeKGf6zAp5AR{#B$_^bK#TI9uDSs(%|576N+_Id=kpqV&gui`peWrUFZ^4NN{VtsfD zB%gQ(9Z&;^flH@58Q_?VIe9@?{4UL4vsbf6ki=$`6$g#l7=6Cs2=gsQc(;6yOxv_} z69djiEvxS;3;KN-Vd1ZI6tZf(zhYhi%OU>V?{m&XAE3!zlVX8(&)M&ub9ZFM=e?bB z{T&BH0w!du5ktb|fTLUV`|%&H|E2N$-oGqENrQ4grl+>b26r5wU;mK;?0D;suTay!PNAfrTLwEu;_+OPg|6-tTt`iOoi<%K`z!DpNW z_=CRMJU(Z9K(n|xteKxPfWL5#My6N*a{%J+H#j>v3JjpbU;AJDuW*0g>#<SO6pEikLCfr@YBhab8JsE6;ob7*Ha$op36omcDlB9bj47;tQmbS5|G_rfx_f>HzAX#n zywo+HyJ|q$emqh6hR#lu^{+GKO1W0|ELr53B?FTpd51P9W?RySKXwocpjRsG#{y_I z#9t4f|L-ioaQ*MM`y%Z-px=Ggd4hTc?r4AMeLbuDM*Mi5QLUZt@JI9E{vwZ-;5Umo zfqzwUu;Azdww}Bv)^1uJwcv8PIbn5sWrWOxV*}M29`zGie==sf4Q+UW{lW$Ijtop{ z$E##aOFOGQ9)8ta>z>(@ZTI%qi3ju`dRmxD`JajkxT<1*;%j`#>$3hzwp`TrwGG%4 zrAte+V`O$tv$8p{!&28azHXYaOdCa~Tw%`g1X z_D=hel79YI$n&iKSO1IigE!Kr{!{%it3_o1s|WY1NLTxhRWx$f4^Isudp7$7H?p1} zfOf*^=Qj(3MXvm z`gLBZp^33OCfkl7koH4j730{Rfu<)-e5SG|De8bo6dSMR6qV79l zzCI}^)Oqe_faZVuwV!H(zgJ$ra?mdioh?w;WqwLP*@D6! z3BdM~XDIy3?icH`gQKgsan14__%tIjF{5klg(}Cctpz?fR(DY_JS4tw@1Qm#f9T`h z!RFw+%O@bZdyv|zl@Bx>eA&?XdC}V0k)*w|#{t>)4zj%7r+p8NhZVvxV2T`nM&6;) zx<`B5SvktAoW&W<(Rpt|&KfMsW{=HTi@!`x%TQ9F{&#!xs+c|=z$5ATWrnan?0@-m zeZhxXTe3rb3Fb#m0RCT+Fgzbr48SaqQGl^{>3doKe(c5Hclwd%HTW0r*YbCnU*tjb z#ovrir~7C>5x;tmzlYm{xA7VBjO1e39bFgSr3{C)38$%6dtpU>fPpuS(^0jXG&Bk!BN(0Jo~YzDG}Ef{P! zqz%6?Gru`2qfRzX9R@oj7?wIQ@0Ks`mjLl!*UlU18AT{fj!AhF^nlF{5BLT-J?fxz z0-^09ZchlgAS{^qAoEoC)GiBDQFuT8>tf&k?~d_x{(oKkhjYKeUdvyaj}#aOp!;%T zy(T$UDjAhs!TK7TX~%hgM)Z0A7Jsbvq890cg7rb!_RD8z^W30_2M}HCx*ANE43e@} z@ST%==H`0`WUmABYF_>pB&udxoMY{oZ*`d-I&2DmHGlW)K@3O~xl1qsx6)lKHhhgy zFD=(H=nIPVBQbui)QB`Jknge8J7_=JQ1~NlNIQta-6KUX=utb$5Rx zMyb)#dmKH#k-zvW0m@a?n8=c9H|2fc&)5}{u_8o{#QH(ifyx?J{9n3U2g%TLD_k&(HMvY9d$vn*IwtU4BVR}*&KJcIjSA?D|AQI2#}wj*<_WAc1hi@y>Bvw$_A z^&L+@l4iHIX5aR#@BnyUt!;6^vfc97!deT!RU7Tg8cNgUb)Eh1B!W*Wt}Fbj!UTIN z<|GRQfBY9VzzCS{z&MCP`aQtL>(u{1x_PbhRnR{q*e;0- zz!K-=(?yY#rFnNjK09>ic~RcMcdTmBN1hJPCvbyrBqx*2h{e%CUoD)8P zOc2=t*onC)EG_uz6U@MMx%jdT^^Udy_F$Hh9Dp%^)d!D)chSF-+Zo^LKN{5a0@)2_ zUKGYk%z_QDDvV_+3Fxu_a6t55oIBj>mG7(byw>icgrEe`&sbe%8d) zeyZ_Acbp~XJ}V=8V5@Uj+`aQJIJq%B91rbAqvr<86QcjiS|p7HK@2%54gfzSOMqv` zCajEhFVAb@23S5)f`vFN_G48(J2v{iYRIE3Zuh{S91y%gf4>%c?9kpr_+wQL$oqS? zqYU3h*QgbMt=DBW*2o=RYd`?K;{CJ{TFk!AFS)&g?q>P*M~n^Up_$1m;K}Q&HGoD-wqGZ z#|s)MgFM@`=swmJZjH_n%b|@2XM@$_*}>}Y$iR}yddtWf`;p!+n2Q`+7FlUm1TgKI(trZ-0+=U)BdN?&Ox})w3>N zyMM+PEiLbTNHlp_a5^EHdU%WB+{lVt9`y0?XypF+!`mMkR9Cg~a&(@LZ%)g0QqaX_ zEX;<4ENV9-1ut-HJjzM=J2tP6$hIJVGhp=<|p0NB`woe}7jP7QkMh z>IYpn{_F5B^!p?K&f^#E_59}e+u>=Rzw6z?M`Gj96YVDXcCvVTMyIgivx6z4eo*^* ze$ayRuM9fa{rSO}7>#T&x;Q&}4_Ck1TceNAg7>+;4%lk_Ekpb%4_y}!)(F9h-5u~m`#0+<>+LckD(TQ`(#DCO~?=RGS z?)-mQ{dFyXwtjfM7YUI4>Hn$uRG!!RuDzJ#!i;~A_~Nx(CE0!KKDoGq+CPv*|6z5D zf&}`1T(Ddox!5~`-dwi2Qojk8hPoC2*HF(^xjrl$+<%vxCz?ZU z9vch(Sn+wSmKF0OQ*%nQy16*$EY7F&aGBM`*1;Lf`Psw;4-%7x$$jAlJ{&Rx=>3N@caC!EF4d4Y9 zx7xq5wF=^YXnxBqW0zyxGw0Ip3(#DVug7VT)8+ZiMY&%5b4iL8p=?I&U9n}iAHJaD z0N_u}bV{?P_s<)AV8m}UPv+{0RLFf30T^BJeY>c^q4I$JJipM^+Sk8$*6-8DUDZA9 zK&SJtJy~DJ`{exC{txAsCq@eBUzYRZtH|}L0Wg0M#GjK_rOFDAhxxHVH2nc;7=B0Z?J^f9~UowM&tnturjaDjjpfE zXK;&C^U3Kg3pl&A0x$vmz|tHO|35#kI0`-i9zo^6H~`y2y^%~Ic|&-Za-gqAm(l25 zvd%np2a-y^5^<`t zmSIlx(OGGSvod>enIvjH@9cL2j`&zG(VEE zDDT&i_-^qDV#SrQ?k5ICGDJtj8k`FVfDZuN{(w`@0}aG`z;ek?AJV#vm{Sfl`dm9 zdk{bwmVmGDb3V1Irw^z+F_>%(= zas7S&>rN)t<$lEgI(v^a@UJfFqrR)>?|M1)`^==G&I49(y*=9VgEvv5M`NwglK(O9 zhW29<=)GCmoD~c3YV&k+O)y-Me=p3)6PJ^TgAME%d0Eh+V+R1yg!Mu8BBn4CxFndd z-*M6S(}FR_IxI+^9F}llb9s)xrn$fQNOMQ5!%@+7IU<}7F0p(_7T_o35RnA3%lIz3 z2;B5dofQB9n9ic!$PI3K%WEUTU*#p%!=&?K)*mKCmlP{#U{@V|3{+tPf3Y*5?LSIo{H>@}xF z2kVFRw_Z%O*USLYfYm04L*bpZLlGEm6-jF!Jq^Ggc)^&Rx} z{e?c(U+u5HUO6iGlkK$9t2LvCn`|{U-^>pcR+$TYV87sWP*B8{mOP+q2WRAX4;k7$ zujdAXCE1&I;C)C8{XfJLa0hCEbB1ssf1?B#Re~NF;X{M)v5~K{!U6Ay9sP9kxvO2}BxlzVoLY{$@)5-&5Gl&{sfNzF>KkjeL z_{jb{>)&tr&Gc8Lr_1?J@v6!P8Bk(##Tvut3^J5~ej@`c&p zP4;i~aCSzBv66w&txJuc?19}v*b=a}LWB;h#l(p)18n;-tqr&|Eb8&dz@5?2Rj~`i z07u0;oS7K%oOlcDEV6|J8DFwTaRPJEcpUnSwcOQ(LiZ|>Cm9HPP`El`7jjZv8J^+kJ??;+xc z_Z9cE(noin**>)&{LO&LZP0R=p80hc6g5Eem{>h*16GU?Budb0VA0&T;+c_ldO9ly z!@j%cJ@ovbj5v84JTzHkct85i98g#WdW>8kC03S?56)LL4@LrR2^pl;D8S=LlWy<2-cGR^u7Du;N0s{e?0NkILp_`cpi zi69EZ-j@Z?{+I3VJioS|YM=O9>Ee32yoT4s42<&RNUZv=wUfPIIpxOclC@U&aY+Dn z1gXKUz@qHiL8piytk{ydr3wr}vD7eEPKTtje4+z(ykAyg)A_#v>y^SOmIK7UanIG`Ij2ZF(2Lo|*zV zXCg}^QHE;-#Fe|E4|M_NwvivfmUVoN>(fivcV=$~nl>pTp%3gB72R{%f2z<{f%E{S z|KMJzmlaU^L0sPFtBUs2f3rRIcw+rK+YkQs_FL_R1Hezv1MIjtYPcHMleZ^Rb3ioo zoLJdac@N3OdKn+9`Pu=HHgjq=L`#4_K_bkgz%11}JIJgoXbJ9vmC%pH<$ z^dI>m-()V)EEaJfI&0<@D{r(fLk53lyx5f}mz2{n%ZIfQ$m)@2A-b}IM*QJZ;P2c} zJ%BMFIe?D)cc}ZS&VC<)Le$_bDF@weZ>s;m011wj9eDj5ehBmjBqeivpaH?VH3Bp5Bo56S9S zdnPBnFz0o8vO&}aj|;}9oaW@J$xy@p zuo%=Ad7oM)8jpWa;;W)uRuA@oY#1l12}Xie5sm7RfbC-Cowwxm>H%iuJ@6-+=(;CH z0GquY?|0U}EPe65LSOAy>y-e}7k_4zq5IC}E8VZvw_vCECL%w2JkUKVJjCtJr-QB7 zS-B`Y{^_Ko)`3rxnzN?;9m#00hg+B2(F+5dU zayMr0u>SPd5b;}AwdWF@abtAQ4{|CX54Maa{D$$x%vzC-iv}rc04zB5A2ZvBY*O|*g5m3@pW~(3VD0L&kt5t0{Pl{A5o3h~C_gCr%9D%_ z8Ki3rTTd(i{#J6Z_T&TX9>-Gg84|WCq_Q+8S8Nt*pCSB}XKRDy*~$_W_BjWQA=) z-ggWOe{`mDie3LG9pr_7=lxm38E#(|pm+fI7w0qj6aVV*EA-|4)!$3BPc54{i@WBm z=49lo!_jF=1y>*7Y#n43JB`$8r)Nhy?ZP7uh&JjuWQ=+Z@GtGhQ|bqJg$!6#8g2!s z2r^f|T#&sy#Q$aq%@~=-@?C$eN2SkOAoO{AKfv>BL_S0P^Mi)&24TvVZlz<9-$MQRgSy z4^P`Q*GP4DZbpG@3cVFrL^wEc`J7;@HPGjmorZZW&bWaIdsq`zCy9o!=wxx^28@X8`yFtn0S<5jBR4#a}E#UXfC1?+vhq5Yu{)*HipKd=FUa%$cHFTD*z?;E>E$bY>!KVary_@c^CmC%(7NuLUUNOZSIkU)}w+ zj)yC(!}jfgca2i>k!UoX+oyybot6<;nb*|tiMZ7^Dt+LO{)06b$4*_wfHz_Ne$dZv@OcKV$?6vp+o2 zd}T64pKku9`G1>#ZoboeG(AK&HP;6HoAY>S7y(w`v~Wgz020TMowq_J8Pt(VWSR3b z`|aEv`@hI#xM>0U2ZI4W#_@gJVO&O1JwO2;@kVt^|lQZ0U9sJAYm(R!ccYYt_ zEALlve>bbEzr)u#X;}T5&$CZP58MvTXk8NjXe32v#iEP=5$|$kF8C<;i@%f`m9Bz1cq=9Z@QRCNO z07wDSWxj}5{@7^y6U{ez|D*TL^mKeSSt7XLEkXaLd~#`|ZovW z^q-N16|mtg&C1$=EzWA8pJ21qQM@!*nq6ab$>5M@qjzsbiz!Fpk7y4zDcc9bcmlhL?3Z1AjjiW72FB1_H_86^wJm6{&}BkLc`Dy;gi86e|#`@%wazn`0G!IG22?>PoDZRDDBm^THlE5ck60XV?Ik zgMWE`R&zA#W6n?C?|gnT&^7ys`CQl?JsXUKYRR&UFd3c;GOEnCI3-eWSUVon!N>wt zHdg(_;xbamg1I_l1+6!BLqx$JEXmr^)q6tZtjooUE3%Oo4G0|5PP4HUnA4T>^>gt>ZgAbgn4EL0I#UvG7p$O4*Uh(pjS438AkU~Z1`#U zwpko%cW@P}2Xp`&8C-8o)OTO<*0)CL@Hj_D##hG&+!PtOBcH9!^Odmy7f15X%DEkx zV_@U3LE_I2l{GH*8_qT`Ab@`}_o7*v-SJKJr&@o}{NS1&kuo8s6n;h+lVRInLMm6r} z{F}_Yb35=2{L$H?uOJ_f^biGN4~~lb+?iKJusUzcta}J;@MU-%7wGx$?8(Gw4vohS}W+Rz!j7zH^0Y27s0Te*hoE z3m^w*m>Uu?BdG72$LSe+#(@rSk4csk9#|_CU?a>F5dZ&aRsei=9{;rYlfKu!$0Hjx zLS({PVO5MOKZrM|sxQ@^55>kHanvBG<>IgP;>bG8=kQ!Tv#~B|UJ7G6x%Zoq1kM-x zXXXpphBq)GmAB&Rolnc$0XebRed2j*2iO?775rN5raq^JlXGS5OlFQ#n{mD;H5e>E zwqNaMmN>JGR>WqU7rnhRJnsB_i%*dl`@E0{RwX22 zw)+`xLM9x4Fe9=JsQLKp$T70X|l)79?l?a1hw(>Jow1K?Zm{`eOn0BWmr>l63k$*HuG*`bGv z?X+M|wf5X-_c?iWW?oaTLl*c9{)TTjE*PDbJ?I9%BGO7#@4TG#1?}Ev<};WpQVl1d z2b6!oagWG1FNnTBkht&@$p<~2=cmFCZpe0LB;mvH2AAf1$bDcR^?CTY^6s!XI~?d| z!rB-AhxIWgW(m=5^n~gWJtTMpa>>SgmH#8EtTi#QFwQq7TV!0fVbI?un%p1d73=G& zyw&AV5r6T1bieTDjwH&ViE}21u$^U3j_rbs%Yzx6tD^;{Jv$u5u&<8M%UFH@Z`+o#-GFg-~d>Fvw}Yn*_qjR1v5V> zxEVFz3*Zb#MeENG`j^LB$yqLq%wC*pfXknleOQYXD?k?z>;dj@Opbp^xc|Mu{O%U< zkB3D*7SyQ|d@L5>kx0VH`6gKf^Vr1emE$p|hXv5riM{ckai62Yn7GkBOwSY^gl@=+ z`-rEhD&kkkARtrXZwyd~)4NI5SL@!-`Pa1o+I&1dn>nAV_7Htc$yLyKO|I8fI$f0o zePl$B%&1)&^w(w_9tsLeTH9ZJ(o}3^!Jy7YPM>S6?ie|(-Ev$do}XOMocL3iCKUr$ zc&J(4cz=8W>wK+@&f*8O0HDuouk%|(&kFwZ@HqDqKM!*_Cfa^!@V_E_p33eOv0Rtu z876;vY=j*iuo~ox_>p;p)C6yc1Q7SHPE_!GJD1>*;QoosU!kx3mi!jEFIST4Vl?#I zm3w4?wZ5YrfA~Rb0cwQWkH1AD=<4!|wdtdvH~x^8lgX&i5HO&d^5LwW`qAJ(54;Q?+lO64+pU3F zy?1QPkXs*=*(wp(03ty%dSr#{F0iIUUfGNwOn~zNA7=fW8QD84zmLT~FKAwz$1{Tc zDXj&dbHmkzsP9}9FLiB@XRpiKoE};IOJW6>M@S6={Ex|XQPZQkcU~(gS7v{@K(0xY z`-A4MgZwM``$F@z{M{94aRm$XA3v@BYumX$kt+w7MJERU(}VwuKl^DJl>~f;S)*NA z{OKkzX0ACv)h#f-MIY-={=Y1M&)}{c)T`kV0Q#)x;u?-ao#nsz%8wV>Zdv*0hIA)^InF8qlU zsT837WDGhs0Q#_k$`M($g&&*{jI{q}K0M|L0l%N-j^)Y zGw)W#NA$`0Y>~6uEoX%-@jet(`DLJV$t;aEq{>Z)6PfcvS}RNzjn#Wj%D67eyYw{5 z^uV8u{x!7Pb(hSMV#UN?nI=~P{@A$cA_RSEH24Cu#&mH|6CfIJB_Y>`!&aXb9c4WL ztiU;Wz9=3Y3jjB`AUI!=f03JX&JVr+H^ggQ8|!~{9fj8V-1Gqx!F&Hsq-`tVm0M7}tHt45;xuIey8=de?- zXQe_S74j0yJ=iZ+eQ7X0KG%Oj#@;+TIFRKx4}?{);^zz{tUu`E;ne@i?~o&@=s=r5 zM=PAb91yu-w1_-IWg5sfEeXoZ)1ZIbj5|_jy^on*r?%L`1kMSHFhD9a_=9sJ3ztL& zz>l>dkBMElG>EN<1Y8%bzb?q%p6HIe?iKkA2_WOJs?~Pp2w@$l`4A;s*89~FAJ|}n z5$E^*vUh54u=lRih#$}Km*tG@elp_)^Re#>q{I^Z!QD9%75|w9HWGmOd;e#nJ8Up9 z0OO+{bjXUo&I-{&>`QCe(8UI#&7FZ_x1c%`l(2=691BR)GC_B;@ay* z=9JuNZLHr6{-O`|*g4k#C59wNZavU3!lGdbMcq!u;ae-dmGfyfp~Y zWrvMsGw)A4!iqoS16cv$oM6t3|I;G_M22|$8-nvC!H-zrghXPQ`z1v3YG~%~=PxY?qosb&QxAM&kGeY=(vVc|n zF#<#fytj;Nq||C4pkLh}odhrjFgz>=>|cf;kCr<-Q%(;;_E-F={%v@kO@|Sa@4qzTcWm^KtQmDrG!?X{_z{^D{@SRW zBf-pgqHcf$pxc!fCIVFZc_sp=Zb7|2*)M$0WxahyJTv05-lnm|Fgc@iED+wn>hOx3 z6LkSH1E;nd>pz_z)UQpXc3-v|V*{8Ud`#rv%AkE$kiRmPea%p8cUz={+`t*J0Z7@A zt^Yqe-y*YddHCOjx%Ru8zw2!j|G!OdQtv~(4bm;Z8be3qJm<8$rTTt||A^p;_ssve zcAy<}f2o+><~yg!Nz)<&S7dja1DT#Df_7QU9kW8~&%?FWig$fqLbFI_A%c zPWJf`^EEtEQHId^YN?Cu9_^$_TMWD_I))FsR|f*_r>0_PcT2=sz>w z;4pZA;sMnqK`dMNYYj>Qumsi_=rg#v_t*wcZ*Wv^RqQamcd}6WPVooM0)(Hs^6ydY z`NG|Z1+Hr2^iRf5Kb~)s`D1Ok3tO$fEx2A9{O?Z=V{NRSc=3)%0;13%n za*#vsopI&g#uJQH)fCJjR{jspjs%#+AC7)S|9X8jzqH>TZ_xK%Cbjll<>Tq_TM~Ui zn~0|A@WFeY9nGNc0nMl5#V!xfuNoFime1bec-4Oe}h~>+l!_~xJZUy4b?#1G(?P>!Sl*o%r8~ucONo19I(jtG8(J zs}{t%;biK%3NSpMUS4?n<-zyS=6k*0YJL>-SP$%&;JPX@eR=S{w#D_XpnqrNZfU@G5d?Cil9M`K}GIfE__u;&UoZ->UqSb3Jw7~3ip#rKSA zOZIEkh~fBKb7Gwym>VPAd3`cK&fXvsNyXo)pP51MzbZTrtGzb8KVOdr_(<&gS-Bck zo2RD7=dk-*V*A&$k9S75uM2B9HOHbdj0CWb2vOw~k)+!rC+vN1a$XPTxXdYLEiwDh zlv-FIYjuD>wr;oR6t)`sPfZ1nh{O{c5+AE6#Qu&CdIajk((LPM5{#2|Gv*0%Bv&DJ zrbzL>&iaFXVPAP4EkD{{wx3m8>GGm~$n2Ce)CY6k^D?^Vl-)FFzq4b>qS6z4Y_8ib zPwwA_w|54otQD_kuQ-~3bQsAZ8D@gD12RV0CN{xNLA76A>B@hs%s~&YoUc~6}an|$w$vD{C_S3e%?E{sJ#Cit^*6uAL5d^;UpAbV@{pIGsj{1$Zw zB!#?>k>GW)4{KBNc_hbrBvNogY@>6;?IY0xkORP8gRHYPsn%-!iHP8ZvSj_gnIbG? zS$W24M#hi5(M48GAt%Tfd#_SgjdZrS{onm)57+$;f1-Zz_b%;Z`?34gpCEQvW9GNh zg+uQg^>kNNf&rWy4f$+V)VUytP$5U!xl^N23BE~o%PXu2x{s}M-T>TIUr_iX4|qx| zUFd)1ea!t(g{3>t^*_;FYJpVfSo4XMKw;(=wDFmFV6y(lM*>cZw!>FeM*=<^?SCK$ zUzYvI2AVfKFW;ko?V`MIhYvBph4~I^i^4C+F=GQyiVYx^xI4%BXk_V`;D29^M=Zoj z(^Wkp4y+siQM>scA~jYBX66)nLFPy1kFHZAk^_M{V+Cg>4+I|E*%186BQmzG8tmL} z@%MgUrSF~nMO30L_(=8-U&-Hrkt%z@tOVPcK^Mju`AzEl zWb`?s5FGL0&i9F!}SYhsh>`hl$-E#)7XZ341{@7SXsq_p> z2-4LfgEe%tcu)Q;N0^y0HFoUQe(*jdWdE1}Xw9$iE@|jz{+01B=I6@&&h5bha8`IE zJUUom0nCI^2SIDFg(v3q>SzNq;?aL99`F?8)_e!Yfv4c}>k z;MtKXdOFGYQt@HN&cf&y_MG__E$pJRdcBhXqipV{R(8PR%Z?b0TICsSaYVrbUD-9q zA%{fvxFY|m{J8=Q$oK2J|9T_-WPHi_81tL)tMkHMWBXxrn;E0)eXdy86+`F`s@_7=nRi1P6Ibe&LXhOJ#4 zt>a#um$YYfRk-~ZQ-k|d@-X*D8*gv%zbxlVj~7*6GDOsUs12Q-_vsR0{tgU){sB4* z?IK_gVtsytd@&N>3d3G8U8t8v!k&o4y%hiV>HPBN!z-@Ko@fW$US`-ZcE*JJlNfvHr*Aeb}4r=_~=ZzrWX~ ztMn7|dsO5KcA@=8Q(cRVb(-mhqh}Jny&@KY3K?9SE`OtX&K!x=m&pUkPjhS}-YlTk zi3H$X>?L*<82yCR^}!*``9jtZeI<5a zZJW`%D(|A>^yPAn|m7 z-0%i^4zgZ?5-~3XP*am>N5|ChsJ92%()SDutwuatsWem zS7;=0Kb0OhIsHz~OXIwiV$T%bT2ZR_%LA=s5EsEN;rMp?kxAxw)RveRLfxNP*>-$U z-#tH=&~3%`Y!H7bXkVT8Zj9akc)amP+PI8(es!c3Ek7yp3R9HB!C;9LU=PgGwU&%! zVAeNTKCpu^(9cX=kSsUzns|?O1nDs2@0v&&*~D)(|I+)z-n$}uH|A5SMEDh0BQ-L- zq8SlZc3~{Eu%KV7fH4l{_VkF#0#;vyVa(_gu@PK?PrAOzof5Nkv7i3GL%g#9*neku zknio*@gwGE{bpAgG$KL5uuXU$)@bFv3#_6_e4G%Go?vYcd_MU459{3!EoD zimGSHz5X9LXM^o=jG7TfOozXxo`@Bo@5hxn=-VTY!@N7@ zZ2oPD-l^u&*>zs4iO$o9uTU2{J*XZZEbR!R%Y&I8%=jd`O@{BxJfi8y8din=Zn6N} zKjheqz$G~cx_!w6-xZ1bTLsK<`&`PRtIWdWgixis*GWma|iXY(1L5co-E` zUl-n!tf-6_FAD3&Mi4JyLyV2c0Co94YCm*qkoRW;f6(vZfBk;t`^@&Udb9h}l=&$QrT5MLE_r$^5+3{HQm* z_stQ%-dh#+v9`5KYw{^wWUMA=j{_Ndv;W$|(uUG6dWP*6GA?3&JAa78h?k6(@eGw; zMV?^RbT->9#{zp}k!A}vIr&KM=*gC$}k-o*N&{n#f~cC+huL9kyHd|08C zv!MIjDlJ!7u{_$ZZ#RPlS1DPqNQdvU)iJ2a>z$U#DIE%=tU1MdG&T{9TE;u CTo(p-Tadp!@7|j~sx{(oS8s84 zR&Jz64eL)Q7rIRii#cOQX58^(SODgGPz@xSV6||5gT8L6fvgyRQXBCSJzksNzA*@s zFS;&g^s&qly`Xo0Z@YAWJsSq{V&cMUvp;iF7sfi06|*V=qUJ*?^258%j;yn0C_Pgk zOdSnBTJaNc82SYY&XhA#%37j$)^W+(uojaXfH@!)`=S4OeOdp=p*bDfwi-XN{hmEu z=sDh6iIWEwL#?XqoSY@=nplyBjUtio{fb*SOOW@Gh=!;~7E$tx)v%^ZE)P3z7bzUC z<`Yo0IW8OZ27LC#vG4ap*Y67&w`a_;@nmzajc(qacOMBC;C(u}nTx@CP0YwUH8_*4 zCSQoY!|Qp5<#BfqR;!V{vU1eYd;)A9F+=$aZTi4$KGEa9;bX^Z?(`+8ORF zbYH4}-^%uh*rO*C4L%>4pfU>QBrb*P2{tx~S{&2oxfF1s1`<&6~N)Fh6wV%5;kcal{I5RubkT3d< zs9rHAj;bniPAW5Lol4Fs=En2NfmVv+7Q>Rs^5@f#aFW? z=eeyW-qmJxumQJc%$apXMR!dk;-To<)8PUSB(sZruwoY+9Q>L0L)1v6pV?v71LzEb z%P>;|3q+*B9R#RnAOW}LIMf5z1lxP_$jS!9i)ZGX_|@7x_=rL)4mywBsV1X zxAsHy4+j)~Fb97)kgaJ>W-9MF0v{&cn6s(^S<`Z=+J#Zme+Ny?d>8S z#SdVW@Nmef5vVi#B?hWpVGWpY3M0C$aG^$wy$*S%-_oa1Ay6fVxZAH2T0g9OxS!AO#wUp4bv{ok_1UHmKC zV{f;aek6s>)q7dT&DFe$p>uZF9qLQ)0@k)?)(5%oMbQ;_zxPMPXIR^`&IR^m|IKBW zgy`?B=rQK5-hu>>*;^jXqnbdUCf#~y{pyUdKAs9AG2K>o2b z8A-F?um$%8fB80c!PyYp3C`X+_`lsmt`uqhXY*ltU!9%Y7e@Sz6I{)yvjJ84!!y~H zK#WNSh|b@Ga~;%KxEeAhbfZ|Sgqu>YGs>l3sd78KRvNe}*8p)=x2yM%L1L4go!X8$ zcj>u2xQhMS#COlep57NUS$zfj`eHJze;$M%kEX2&3O7an9?833Ykry@Jl1Mpy{5m- z9`y8Fo{df)cm!Ep?h8W>09(L(!iBME^mkJwVx~BK4)i!u5pw+mA_V+by4DiudrJM2zk>{B?7EvrQC-)j#6bq04huzr=YK%Q2sISN|ug|X?luyN( zu{562%ep3}F_`$PeH_ESN_+smoD6_7fVKSA_|5f@73AJP*ncB=Q z*ddd9Me-J8oLE7HP6#41?&pkuvRl2fQph$Dq*V-anOqf{7OZ?fM3eO7?EI`}^k?)OyK$1U0G(#SD=!*n6SV2_9_%xV2US`GeY zMv1EzMjEL#pO`)P)Y#8H8Y}=dwECf7VC37VjM5JxJ0zAg0${CZSN|Ln&tKV|LZ44q z8M0f&)%Cma1N3_l1u(-vOwnPZ4ZIFHSbQ|PPi~&<+<`e?Dq2W@TV*$_XlVCk1ckpA z%)2WPZ*y->>$0FvZ#J{2;oQ%(Gpat7`83~c{;~O4vg_Z8U4J!ye;w`qS=QrwEqeZ= zX!)0-`5#Lb?v>^{y|?t<)q8jEz3Ir^w)ZE^-^T)cILBh1F0+xz<=)hyhb{4dFo z@s+R+t<1cP7g3s|iv7)^)jVH2z(InsBS&QCY==J#fR!KMe{w(W z{{Z8wOs_T{^x1~5-;E~#f9H9@^gJ_tIy{v5zI0b%tGO$e>|MTEabV?(h!(7A+B;Cu z0jz+i-^w2Ihp{EtS?}MAoxM6Ec~`9a}t* zG5;PB3HaCe`XA?=Z>0a^(~*qNH{b4!=>0+OPqMwg_rd7@-^Xe^9lP;NFkybztvUYP zu@&SE>G!@ac7Wb4DwuXP;H&K-XO#)~8(x5}e|Z49e}9hoTx)lj@4Ys<&q^?lb~^)p%sVbwQJbMXdWY5qBuf7{+xy+7-{Gjl|KnD;)_ z)`PB(y}vqA@K7}Osi5>iI(EO4{U40(;i1h0GAE4q7S9bMg~7v0i3pe%eo1Wp1F`fU z&+ok)sd_L+B-$k3e_#C2vpJKOat)8<-xub4$Q1L289P>;HU_M20C@g^d7s>r{1BE8 z?pS~FW4tm4fTv(C?3^Jhq93tqkQ`9|$}j3Tx-a90^KBbkwh!j4;tulo0y#j3Kb=3! z^E20jwV~gd73}H$qbfs>J*#`$uOgS?@iOfa_ZO zOg@NvLsKQh12C7!S?$^Yst-g1)cf%@pUB_iiA3%Q2e>uTi~K&=N-vSe6Sw!g0(S~cN;~|)Onuu0Z75_@vV#{`4{Uf%8%SY2I$>C?flWpZ{E)4buAC@ z=RVC_=9!uOg@3n_yZ+w_5cz0&eyqqbiwpd;1&j(S5IHN4cm^)h<)Nt-*zX4iq*7@f zkRwn*hWWEvB#etbEZE*<@vG!mo()<*h_(M-^IGq>dz#`% zO3s8;Am9L9_Q@EK3OTyttS9ala9Cu7-f3#YH6O-G8omH;0LQ{6Gxyq=Qs!3}Su-%Z zMCZKwtL#sfXPw729Pt2P?=uY0jvr_E+xNvxPcwaB1>4|!i0Fv%2J_ok17>(S%L~b% z&l~1fmFTKe%LdI0R5q~g*U5dgyc-)`v}9T?3o`cvEu#1Dhu!@&c7H-JJi0fh_kr9K z_#Yw@-wjevW%Ob5ti^w2^y|XtEE(c6a)dh~4KKG^`NOO{_Ul2J&TV9iHQH_o`WNQ) zqu~hO4PSUMw&F*Tsh`A#UmDb155`VHqBSZ{v-9*q+6cVFt`6V?`Ny!KgC8oxuN~9P_073O4nWL=@-bGe04_=NfrN zf<^;Wjp)L1W(9I#C4fl5PRoA$)qni|a0yTnkb}CDhv;)}8GF9z`=q~=<05zdlGQw{ z4-tE<#Y0`5m8Z-On2D&45h94n3t$2Ixu5rNW?bRM#N+gG+?vtiu8d!gPyT)^^B?zi z=&zg1d(eF%`u=G4d^o|T z6IvEYH~#I>#Y$@O!icB+gAh10&mf@4Z-MOb{=NFf>+Qix#u7Qi4FD)de4*sr#*Iy-|g2rqj zIswNfw{uzxw8n~06s z*jGk_sSD7_4Wnl65$qc8j2*!?z$WavDAr}KCw*AUV~1}ILd?|rc5J|F;hHdmJK|61 z?_+KN`j3nq6Y0VOP?t~w=pH$<`Njt4|HhyH@4x=7$KG;tW*>7m0J76}WV=0Q$qGQ9 zi4+r$J`-#3Xlwv0!*O3p`oUQZ3;Ykv*}x5t$}f?D@;*^`cjt^7S)uQ0yF8$BzxG#a zRqNg#6WZv`zCXSHW^MdW-q+QfN(#JtulR$$(Z8|3>pWDg$y_un!_2O+30{HTeD8jV zu2+|*K1$ytH&OpVpYI|y)DP*hqJxu)?hQfqk)ZJT#PolVi2Z}%d++Q0UhgffZT?8c zm+bznLH~vzWZyA07`pM%Q+oBB*G(@s^ZdC}7c*3AP5{zD^#@jeTr9(#kv2N^z7^^E zd~6Q0_HKz3Uetcmo=qwNWV*~vlL?~gwIauPu=nF%x%};~y=}~hy<@BHieF)#2rQ7b zBp;4F;m*Q;5&eIvtw_BXDWGylJ{TD!ngUnu6-g}xOS3#<0{62z!ZT5T=)?Do2CxRy z&Oo0rMl&hy4})vcKW@(~ra*DZv=VVGSBC zN>12&eJH8UXsx_ndiGx9FUcZji{&e%Shaitam`mEh^XKcbfK?Jm^A5opBL)Ufv=-i;w%SsG> zoZeBQLNu3Hl1?G)1@juve%27s^E0!Mjwj?{S&nd4j?2sg`$3oi!kpquTWO}w@KVtK zbT;N7t&L=|eh4d#(kaS%SyllV59f#58P-K~&UmLnMvk006;uU?^Ucbb8-x>*6{0)K zRbTD=mHWf~3V%92Y7_rX|Fr;>`(ZvJ_ZYX^AMJ;Ia7V{^t<|zi7v|T=fOZOr5NG6* zbF%`ByrTWqEx0!h>jmKVU4!MxSYUMj_rvtJ$_gs)OWp1*$=QD~>}hSrl^$MJZeqiI~w;L1&4vTmwXo9=ddaB2)(L%ol-JJMzR{=jdM`Xn>t+pK;a0@EJ zvIb_hljouWz??DWmSPFutoEPNT}J)z>_`Li4mh`4TiJXr67WjA!1K}nJMun^lDiEv z&yYJpndK`Vk0!C{H)Tn7)3e4No`-t~Rh89sQ?S2e-7DuR{;>b5{13<9RoyD`)BacP z7ya+>H~uGLLwB8jxqoaB9Y638-fcwn`WkZZCZz5T2J4F25Tm5vSSHFP1=p)%$N$c zB?Ba?r@qR)OR6h~{7UI1-?fvxtNd5*ztz7nt>H21zS(9T;HX<%6uy^I(WMIINRry`NjlHiQ z`aYVVkn(Dyw*@*wmI%2z!Wcm;A|5{7qC}Sdg0G1TZcXdNDr1ZfuGli}&vuS42VcO3A z!TQ_#TlgFMo8=#lzjJ-6@;e;=&i|ADW1azC9E%Rd%cq^FS8VCdvw zZXZ?#t*)*R@&$AwlldVUaD7{PNKXm!=8DW+-Ia5CGWonuwmzRq)F&bXPedYEPl^mQ z?DXpV!U-86YJiLvQh?Sti-vL3=ewpQl}DtNduV|_qbWx)-Zu7gRt>$MRQ{d!Lxr#K zcaCSL{r%{Z^)>pZ`cM4t9UaO4@>}K}c{J~UZTGy*0(thFjtoN?BT zfdR46*G0zU_KevJxkEscxbF`W@$J?-Hgl(bJ6%3>?cJJjlf}WVxYG{zTf8Gk<2C5A zp~Hf%9y$ex-*3o%mqpLKD+i38yEEZ0Zi$UxrvK%!C(IFN#t(I8m>W@>_wi)i2N@&l zI&v2Pvxj8=$uQwLU;@tjCQ}3_W&MXc^855C(U<-0)4TG!X-xjX-w$uIzhJzf!Y6tgi$$U5=N zr?@OS_(-&YIEmR*_XIU6W#s#RoNWA$Q;Yd(^qMu-xSI-hR)tqz7bMYTR$(O1!d*A9 zr99eo%W5CY2BJTk?2jw6!~W2JxZaYJB=5avq=Vs z+&qXA`_Tu)dqje)aL!CL;!d~*9_g|C2D~4yUp+w1NyPr&m|wmr`>&2<(d9!Y9sLD! z+7Tt=@2Xa{dZ~AuxisIoBKR|R(8?*?8s2Ahh1?%2Ke6r$-e2}_=jWzv#7F;URRD?n zz+dJ^y@$CX-BvTe=>B9r0`f+OG3!1#Ps8~U$N`c7??WmM%4_peAWjAVFGbwvz1|Pa z2%elVxN9iOc70wwl2Ljok>3|$F~5|2`6r?)_XjEFTCn;S9n|Di(0$OqD!$c5%(Jf;0Ijk#g*dW zio}Xoe0o99NM`<$5ndMSa9+O0$^%Ft`N)T2VYm+n=(85kbFKd0lXHP#5*@Ir_{EW| zrS0eu4cSLUo=3m$tXDFPE5k70Nvut@BA>$d;fbsyYF3r_n$Az~9~b<~`{U2I3dZ7J z_TLNu9f0Bw_hV1uzS{Kv@_+8}LVlXO6q0AJruW72-sxa&EQsf~qfO4B*QZ~IyXSc) zb>`DKXA^7rSp5F0!TsxD|DO(PdLp`cb3O-mJUBTU?jMH!GrHiZW|h*+f_8H<;5+wGRBT z0!IW9){CGo82qUsaQ_!N5!M9RH902T#}CBPJf5pzT`_+3sr(u;MQ=Kp4C0tmB1>o! zu^sEflHDh2c7$EOj-FHIrIYtMCMXc=g93FUEwa|1jtf?j9^Z#QnS6Tu3jfaflLbWg z`2&BrA5ot(2JQW{{!6D3U8&R|;rR4^Qn95gl6n{`ds~w*hhdKhY(ebRetI01WTadv z^4yFAew?+kJ|5J+leH2?#P-uE_eA_RSTc8;96mEV=(L2(BQ>|>dsqgzIs6V5KnEu) zKcIIoE}{Xd@KoX+%sb2tqRWqdvMRVK7U3Uk$J-C$&*qu2)?2~!FiW}epbE7?1 zWOP|!le?@iiBDN92ZB9- z{evy*Leq!`K9X@KE>ylgpBmg}g287}a(398TOy6vg1g&uq7#w(g1wRm;Qs7GUIIVG)jyVf zU;t$G7)RDI<(|7;CEOzi0(=}FmaVNIPtg7-kE=Rf3QcpEBYIi%#YDu zRe;F;sQIAJhP6QR_5HOP;mn*3ohsOPywjYV2h~sOA#^IjT;ZgkT;0qdk35*ogfTA2 zXsyYJuL;{CPF|C7ekNGInyCM+y@B5UN}l&iLFaE$hy8M_!o!h)$Mc6bKwqfqktu+o zyM_z*bhtWWN~YzZSoe=bcIe}w+k(}XUTsnS^BnW1=^6W%^b>q19Uec*;};?Uug2Ou zop13w%m(>Lq~mk#@7^1Ux;u8}hMc8$@FYvXd=Tdc()rBlV05c-Uo#%jeDjCeDAobR zzuX(tAC9Dv|05Q-zqK$(3oPNTFh}eORUaZfatf>#O5DdN)6c@jU24hv?GsB(Ee3|r zf49ua{C6w9ia!zn`>%XIYx`{yYe3Ey{-@n1OI-MiKGk2d|5I~TPo@6yP|n6aA68!B zoLTF|bvc;J$V>+0Lz_gdm;NxK82W>>Xe*<$*G6-Rm9B2(AJ6~gApXzk-2Xw_>HGKT zzf9EP)^7H@^^B`(l2Dxa{$K`;oljn=j;OtoDC@WPlkVtQW&7FLssS z1*jWW-k0mSKIcX466Q$d&-EzCK$6vfbJCmoRL+t4s83}iVV^{>^3cL)#d`^=i>I*V5A8SJlmmE1AH`{CDgRJS+EC^}f-u0rr1%_!}km?d9O^g1&MMqzZR%1Ky&Q}C`x~=Kg#Zt5G0{3mXE+dJD zrNa;X|9G_e%kB2n9O)mENq%SV*K%xn$3NKHzPD4ZWXs6G+ah0o6}It{bQFD~^#kzq z=fgn1oqyqZ;3W@7LX9F>vjq!Gu7?>M%+8>Pg#JOY3-DCXl&io<$u-`bPhbV)dax%S z&Tr6N@IXH08>|N9{jjMlI)da+T<^q{#py$4RbAu*-VfU+7elU1TT5mY31CeqSN*Hp z5BQh$*Z-sY#{O6U^8WCCvDg1Q=WB<=5ra7%)xT-Y(OJjhh**HrBb!7%!G ziLFA*$yjlY_9nYNG+M)4XXo>>mJ+#g`tI+IP5g9p@O#Ply*+yQw&3)4WPc(Jbag+P zaYF*&% z$}f4PyJA7^X)&eV07oIZp!4wFT*pUpj;uO}zo!TIw$^_k0o-eV7?Qgri9dcFcD=ln zF60jD05>?O)hsJ<X#ezwY;dW_R2V4=@~ma=%^HpZKry{!@bOw8WCE z0I_%E^p1=+`EUAD=#5@9#Q2Fq%~O)IBtqP)l}0!UYr~N#cjg40)2^s{LG0wM(bgAY zf4`H?{vn=HUHmD{cxY$Z}^tp)xCLdRk`UNTST8(|2a3RXlKYg+u^npA9=1*+W#duIPo8DjRwWzC&A?|pl}m2;&lh3FCP zPYl2s1ayBsk*gryp*lr}0XZP-h4++1yIIwW_{#k6;`|!3-t7@a_YcZ=qE}>k>BS<( zLJG*HEy!!K!0LbDPX>3J7XQxs8~v3BuQ!}5&=odJLZ*a1wA5WYu^r1~c4@WWz?8XEI*?EUv*f&V3_zawb;evtTC{=OGh zNCZ#5jz}2~On#2BrOHkR7}XiNzhKeG$)l|=`C8cBU&ph5JvD@Jy+eBE_HOIF*n7Hn zT<}@mJ3G>_xOaH(q~7td13Skb>>GLbQ11^iFX$`zK9R_0!_w}~Zz5M@6kJ#6%$zw_ z&&KRO=BF|5n{|aQZa&gGc*Ka_;jtj(D$stqQ?dKZ3j9)7#nZVCW*3vQWDk58l}P$O z&=1hZc9KDx6KhGnk9*6}qiCg-HR72+W5pMyO{~3tzI%N9!@jN0oswEDl^^vVEVqt^ zSMC@6uY0=oukGK>_!586CkCkX-dT-iFlR-M#$FF&*S{%8F3&*Jre6nXecJTT0g*>%jZCl|ju z<4ufChK76hp?Amwt31H1KOc_&!$gs9O}{|+WNgWV%oaJLcUkYI-X*d5%X;_pUhDm+ zx3+g))@P85=d@*PASlHPob2K=`ld;=aA=Z`RzJB1&Y;cgFR!rua z6<{xDZjY^5)%#zOMDh{Df94UeVD$ffA|C(xeCM9bk=Y; zdAgUWI1%$POVm3=xL%U$YZ+BhVW)$b`GC{_rbb_>{F%c?`wM@V-^lcT(E)}PApRdS zx68XOt{Jx*Rzsh*63BIh;MOOvfnm#Y`d1Qm=tRZ(oB*^^G zkLEaGjhPSj+rja7BS9}WU(4$;(eQ=6Gke$cuI*jdd#LxXBlaFKH4^Y#=IUJAySR6K z-n}8xa6lwwhx8b}Be5Sm@{2i#uSEu#LH>Ajcujtt`X6h>TP>tdl~tzcE~bLXoz0)g zxsWNOFM;eKJcSsK2#}bc6{lTO`sAED{V6wQfBOHNSw=41HD_4YjGhzoOfW9;{q$OD z{asgwSPN@`&0sd5D@M(X1t|QT_lu7=`?Eg$>GtosKKdoV+K;P0f&YO)p8kl_Ba_}; z5B%xirhnfm6`qJZAXxzV`j{Kc`5zvuNrl|5E!J_MI{HfV^n20PSo~<;5952km$>im zgXV{WFLP_DzFpq_QL`nZPd=AyAXyf2hs@;rZ1ertkKYT&+vfk{>ps9`FROEbpPAW2 zj4gI^O*Dy7v3IW~RzPBpXkv|`0&1+-P-3ir6&poRiZm&LbWo74A_7Jf4L0nGy(UqU z8@=yu-t+mNeK2|EIcLwFJ$v@de|@XH>s@Qj?mN5L^89==zU02~{_mcBYkq_I7W2QH zUu*v7^X=!4n_px0>-mP+m9qqfzDfrL=*lIyU5TMr2V1?xL-TuU9}2nY7-owb6=fFCqkZ}HpcqDB)@|d4fVl} z5&2(~tNYdH>VE{4`E7EGujl)3B=Wy9C*$DdGwD1I_j_{nU$cvvFj+>qls7 zqK7|;wm)R{=tTIoEq*2wKempSzhu7k{LuM_a(nlD+xac$cc0%YC*jWXznwo}zQyc& zv&&}Z=7hX5(cm7lyNKeb43lq3LhPC$SR*L%?@t#`s^) zL8n8>_~m_=Mz8=;{=Rv?Fl!RGC&-qV*097C(*%0@!J?=(kB&oIqp<)|8>WL8h!+Okb_*%YiHt?C8j=l0X;Dep0 zNz>(0Dmw5L)sN%%3&C?(F-svm*IjXNL#fxL~%=?B+R7tc=cKoy2G7`?KP! z_RqT)fp*qArWk%uH%N}DL9fgACKjkFlK=a1pQVBd@pq@>BpsD|ep0)D2VmNqX&0gj z%t2Yr?z!I@V<(-Otx&XI%s24$;;Nk>4Q(0H@5G-@d>h`s@|XKf@*ff4r~uRgGzT0A zAPflm%aq}##!l_MJllINqyO>~x_x6z{t=VF157W#l%oGKfL$Vg@w-}Vb7hkmz7$LS z(?omUORj!#r2olS@=FtQ;ra0vd*=GUQ{S6w3Wrw-6g|8zdhcp}F#7rJ4Zn`}xKa8< zo;Z7MJo*2fT_{rqrej&{y+48;(ZSfFvWi`K(imz{vIExzBs4zsCT@@srHOi%}$@+E%HBa ze);_Ci?zj#7dKnne{t~Q^2LFR9TtyV+;G0{{DJdZCJuO8a)bAzcjvrJh`3A6q^doP zUiJ38yfgfrBO-tLKka)lHO|&gF+l5?`p(P%HuVFtkKg0m9v2((nVf7BM?R99EYnO% zmBABpMfiHtBfUB-2-XMxrBlAb&Yv0EYr>v+w2)LY;nY?zE86Xk^*^?LjG6BP`y1JR zJwSImk^hJQM(wwpKbSxKuoEE7ephW!Z+b!A>5;i2AIa6gLbUa$lX4-!aDUp5t@^CU zA4Z5RAy-r>dtEfB>PL-RLX zX|v!&&nHBSF3NBGYPLwGcb_0T+eY{QD6f|JPqSOh_Klx^)qMZ?S@TcKzqa_J#jO{& zT|9hoNaTOf;)KO37k@rKB2mCyX8#?0;6=$2&x|G5Z}y0sr}Lw!7v-H?od0$0UKl%L zYIk|SI)SkOpqOeHtqt}4z9~^FTUADyDo<;Q7skGQDX062e7!uk&*d!03H^U$J{Lql zj%%+pb}uBbwLf(rSphT%4hSdEgYe{Nq!>`9B`;w6&1W06|7qUuYVJSuzcygl{}KJc z{vHrpFQ0hGQu>|!@Q6sCS692=C09X4gS8+IX7hH7bVVrYy}VTAkFg@IdFk>zV@R3o z!J^C9VIn3rT{l_Ry_2cU69axH+13g1mAYb0lUGnP{AMzdEAvms(Ybk7`b@?C@HjXC3~0B= z9xp_9xG9p}7X}bUGW&FVPJ(Fl$YrK`S)$b066Lx{J}q?n<4gIQ4r*Qi9^$S!S6Jnv z%f!S9n6UjyxMrQP2nm!8z+Er_rrF```r9&6Y}fy={;mB~`py2!0BZ$clK=g3j-Hp# z^YWaiR)AE09uc`hC7&Mq!z-CmArCewk&eDO?+qiRIe$0u%KZ5KMgZl?`K|J|4_(IB zzsmfcdj^R;Fj4zF_EGPgdWcwh*X0S=EjEz{5dYh_D>m|1-yE-cT=evdi5vej2>xFu z&wETJ1N>34h+iZYeAeuH^G7GnJ8b^p`A6oTU;N49af@dxUbnb#?E!0hukE__gta>_ z{$+8G#Dpi!-!Q*@wzPjKsK6&?FU#bYA4F5Xm+awcsYhHdn8fFj5gwoOw^#n^>+`yv)`)?U&{Kn+k`c(!Md$ zS5KE)n=&N=6lJsiw#;@WHJ-)emsw!hIB=gao4Gr`D*tx@U`Y^07uc=3VM|Q4doAe@)(Nr<9rHRQa3bBQ|jBRF6*CwL|GG zxh@}$?0%a1%e`ZB-P!v--{sQf;@h-jk+IkNAkMply_l@|F`_G<|YR*l9Bi$$%z)NzK z|9yV%*#AAE|36OLxAWpFYj;}v=w_d~+7VaZe)BW_;Iz#iy0+Ka=8Ly2UbuM7{0s9} z&$pZ3DEa@6u}AX_-%eF<%gmU1Kr)H%M^E(?>)E5}HUy~3*!L}$t@JVT?1nQ_hWBsX zUxoV3dFOfoWGcUjojWX1#=*0#lA)Xz8_^1}Nf2;BD6l%Z%6})%msg7spnfL1W@TIJ z?ViM1f3r5l0F{4hziK*VepLRA{RaA9k1*_i=YF^H$L=3O!NQ6^LPQ)32G@G z?X8UevckP%b`1KrQ=%3;NlY1&-LQVojps9$?<=vkqW|*~Egz8UNBVs1@Bs9m_wTfC zl@>9wDxhubSUYFub$*jR(Yq(#yH+Z^+b5EL&up*RDYNg-?=XMK{OrYs#cdWZUtGDi z?PmA7+KvD4H`o5wwf=Nv$ez|Ux=nccFA825H8>O61u`fLGU2)d*=|8g%0`Ay;uD)hp;v^NT(P zM|)j+vH1a;U9^<{0gI<4F4$%95Az?+?;ri&CFkb>L6Nsj&ES=>Dpy}((jvekBi&By z(kJEwumL(lb&XW}vCQGqY7;60x(nZ!$Mk-knKN){Fs`4b!u9i5E%~IWvnEblnD5>g zeQI`xM_(3=o(i$Qf0$laH5eYIgLi&ryXo2mYU1-$`&0EF(x>a|wO{pjQ~7u1do`cm zZ{645GUX4Z(-}Q5T^(@dmhF&t3W01dpbW5af?1B81_u`u{n-d>L-xC)4VNX`yE6Fd zm!n&sjU7~7lU?!rovN!&>U&)|B1SV9`gjZ5E$N zKKQo9ixUZ+J-G5u1KB5$z0p1r2uzmi?86W{-SLhF`zE|wsm!btfkNhr;O|+p6lM|9o z^L2Q1gVr8-gs1IvAM$r1Slc%x#`;t=^JhWmf03x_nz7h-oc(nE!udVsH%vZ%K)m#Cl@zJoOj6DAFq94vy=YdqgQ)!y#Ik~zg|0f?eEuSi*GK@P6Y6k`7>re zkN&?avBQ&s1$-)4;2t>%e-`8U#R9Tdhi#*m}TgqR}P< zy&})jRRp(o&UrAF0|#!e*cD8`ePSu>gXy$KD2g1~HPi*vKH$#OwtM5@i2!5*ei8uaC;`B`WUgp!+NG5Ur`3{Q>^C!(dnkxMNO+~U-Vy7r#6{nmC|yVl}ciwhS|4n}mB*|#F~4~7}^_|$_xIy)kD zA$uJ zh=h+w9CKt|^Pi$wb^(+@-Arr8rHoxkEDJRrGPZTkR)qz`rm2qa9>s5 zVf%~xlK+_hvnv1gebD}HAN6#E?d{~hbb`6T zt?f>c8>^Qdm^}F3W3R83D|<@pxcrH)=l6B@V04((RxPQP*{k%v*a;azIbg>nEC1P2 z`qxUG=bD-CeT~^S=DP*Ud21rUKS=g|qhJ8%1=IQAe9QSIi~B78A(-#+Yk#u#+O><< z&Rg4g?J;XlUE6x`AB)o$?_X>;KQ(p1Yt1(c0`#KTji)CLc*1Opt6jnMmP3NuGWM6^rpzkM0UzCR~N7a{rc z`K_JS8nnNcmT$Y1zH0nq^BkxEu1F_B@(bOgy;Ur@+;1R1IG^?PmH)8+Wba?AfA>nQ zxYHhL`K$YtzP`fB|GCRN&i-=mlsD>3aP>$wt5tvOH7B5Z`|ucd%Bgu{PQT8`;}iS+ zCJ4>-BmWcPE$xdI8OZ5W^mSgtRQP&&RJO>$k^cL0TI5@11f3Kd_1{zV**toFo6Ph7 z(T30E>mjM+KX&%V$=Uxq2;fVj#h;&D5ex91#pa8@%FRmMu@KLD=UpF;?rv*j0Q>KK+rX)r+(aBVy<8t;ybY&cMZkv2h1$P9QBub2@-Sf%EB&xmhscaUJ&zX2+a_kw2flNq$vP92E-?#P29{q}uS`v>=j_iybF^N;l(eE(tl$Fvsy z-`pWtU)gBeE#Hy6lYQt^XdX#V0_nFu1X?Ip{O8E+pyWB{r1oJJ`{rSJTpF*se;$XX zs+*XYqC>Zx_o9DOTK8SzDesTP$N4lZ3C4CxZrxOSM#SFIk zU$F*vNd|EK>}#{1&2Kn={QSuIR}%$1FjnA@#n0CMV(n&&uP;tr9KLwRV!Oq4=9kR( zn?GrO{rLRXg$cMtx==P>)+E237)AH~xq0W9I+#7>o~Q-!N}^Qk4_kza3cf5?2_GJQ zgX{RF{En$M&|rHA-;xudPA)@!X+Gg0vs=t>F~8017I}4g_2QhSgQE>L=?=E9jL$Sr zTeq5x!Ts4^J0phuHph|H}M$2J%9i`m^HrH`f-Y<{QNL3c?ccmC`* zsQ}+7=-$@zePRP1n+m`-i3Cm#^7rP&-!GmO`9Ehd4@R`d{GPKfruzH9{N9$ctzrXi zm^1L3*p)9PKU4oVCtM7I&Cv~l^q%Jz)#?Z&fBUcDv*y{^7-=IgqzaD!m;Lbk&MTt( zy_#lxpE&#W+8fr6m|rtd21fmha3jPYyPa9B z+9Y>%92v4{%O?Z@`*LjR*CM4a#HyOGidV&>+op%z;VP4{oXV6wkK=h%_vxGH>s>h! z7smp8HDX9!djV|!I4qxS%#bZ7?tQbf*KWCXkJ-P* zzr8Himyb0OcKhWO9;^;keKXihTEzb7wwO&P{+RaJT`%O{%x_?S&HcLXakcidsr#a>pgQOJ{sSAWjw$(@%d-Yo)#QH z$Nfb)Sz(Z+o_O!X@8`~MGv8`{&3N_WW}nWz15yus#r)vZgD;$)Fh4$+-wEM>JSv>v zQ-U3QZ1&Df8oggmz@JC=|9hFCek$kGygoc_7RnCsgCl=aWgDr=_IHomaRJl^UY<{n z`?B-W_S+VDMl$>h^M22c{B3^`KXq;_RP2p8#UGd*6P>Utj>F@_SgG2@^8lv+<)Z$w_o0;|7?G|KU983_9yd``OEutlsEbl+u8G0 zbh>?>*l1<`VtBii<=wfcX4RO#4>jf6-p@?Zx};^o3^3;%-s%IXOe$Nk#1z$=4nltAA{scSs)nbSAlPSk{A|8_RIr`JLk( z&WvOa3Kp@56;(adXZdy5RP{zlce+F#9|cwfbSvm*b%a*az@_I++H{KMVj}>KxL??*3YR}K9fta>?rrnc8_5BamduXC;Iw+@x z5tRL95A4YNJMC-4exyINf9QX$f8{Up8*^XU_oLfmV7`s(od9Q`R=yrycUMsd(dU!% zj$5gg7hw9q2YouRAC?^ScE4D`lXE3Mov2Mdb$p~{W2Q*TRWc!$_cjx~9sTsn%v~4> ze63zqw;kTyo{`1rv8kVre}(~Hl>E%TReg9TCKJ4N?D$K91iWi@bQszt(+4e{O^M-D&@{zxE&7AL=_cbajG1FYwr$ zIh9{*v~v7<`Y@0&lXL_5Ox@6pr=Tyk{J2iSJ~~&?T!nLE=gsj_&+CNb@+0&}n%->= z4?IYej{~dUrZ}s!z{@%h`fc@8>h%8nvVw3_uCHnRdq)Zv1Q~c>r2ebfrP;jx%IN+Z z!Va)AYNy%R^Xtv;mpuQo@dh`^?fN+l+sv;s`(&_##|JmEQREZx)u-q6b$*#Za#-Zu zeV#b|So??O@0%>ACcWQM_D9Au*zLPVt^<3B$txb*C+9~0>f!l3-7sx-l-gM+ zf+R0%Vp&=AZrD|O!R;ZD$42k$R$=>b{#yAf^Ji@ND&u#={sZ@E)KAwKs~xqLkny$_Br`Pu9R^ZSPF zef!veZE|}^s`@uih3K!+5Acat{huTv{6u&d-;4$LL3qFK4{!7yiT8h<=%lm%&FE9Z zH=~7inxvt9?4g16-7`T!wUw*K+>WkafU^nB{dYBEGiJ-0`Xa_}e+T_%f5rr7GhFopkDFNqWqX*V*Zth4^K}0T?jJ3+7rVV2 zvTE`ly*-!^W!Sk|SOfa)bdc-06KDBn*Mlt)A6mu~7ezW^a`^oIvDyDVJ0+a2&F5Rj z?*D1(LmQU7yR~GCw+*iQrql)A9TwIrX8T6+yXTR+29Mm3zLWn>Md9m-{5}@_$IWQJ z0o&gD^d5WPyESvtXBxaD9cuUi7F=2fkKVJ*8 z|5m|qznw?siS2$82FBHb7(F&Je^}~6FHKkIlV{fp%Ja=U=f}Y$f0;beR4DUiOcfBb zb*@|^z|qmAca@JgGVdJ|wtK_*E9_mpF6~{2_nfT%$oz`?RrhbZVT;)RJH`6re0Jtr zXZ|(%8`&TKQ`VQdul@tEr`V%i=iS=#jQB9QteekF4@^j5t1)8q`-t6JCoQMz4LJ=T z$o009QC_7F-fUbJSBxy~WAW|pdsXD!U4Ehoxu43qx=?$_4~iC^m)mi96}AI^>}`?Q z=dvAi=g9uhSoo*TZV-R4NS)vgvm1uZ_toWW@#`c+sV{Vlvo;c_%tQ zbncrVcWC5q+87U@3gM)wQ<+rWpPS!mCcu~qi|gjyV=3dS?w)h3zf-p~Eh;0XyF`nb z>G9Zpqx?x-@4TG5ogL1<>irRnO zek1y`>8EUB6TkU>d|wfYY6Lu})jo*-vvQ^z)93@CqhcJ>ePCgFy>xI_V)6k@;WtN8 zAI_D&Bzk5;`a!v-XDqVt_qE3X4+!Q-&z+7m0j$>8sTS$lu?3 zOX9Wj(gASiSos&n%kMb5Sx}-s45H8K-!K{8*K-Pfs^dI2Q=l%*-`A%BgVgzhRcP`X zwvKMSBl8tHDI)C5LH*qxwIv3_8w*6n0tptQs$5{I!75)BttWT-)43j<>4iDV8WcUj z^KOUlz})my>+z}ntM_`1`d%kM)kkK~X+Kc^A^%QwD*D^ZANEK8$^W*?=7V-U;=CAAb z(@$&**vYxc{<>4tYzcG4?3=5Xan|(qce*r-q)m}lXKg<&_o2%61oOQ(z#@g#{9E-Iv;N@zkRO|U?;SmF z_BYA@b{o;ZR{rezYPa8zzs$3_M?Jvf^DfFlaz1&cPJq|x^UTPfr6V6b>rTD72=C7o zJuz34{KY=@Xm^t?_OEOjcl)%B@;`oQ{f}5`Kikgt%OP}Fb;4G?Cv*wUh)-qev%WgB zf%Jvp2%Dv=lM1Kivytm}H~ctW{|52t<}Q9MoTKZfqVv_r`{G!i&*vokC=sEZuV#FS z^XR(1Vv;wVK^KpiGPZV#Ph^Ao*YrO@dqg#uAi7laB-=1x!mO#GIN-QVZ_HvFA3L2#AHjzx9v;7HI+{6ix*{JvfDY{+Oa#BRdugbbA@x%I1~ z=db?4v*TYnzZ(kZ{pwvaYs^eoh~De-enl!W5AR4Mp);Ak@ZR>w`^Dyv z>xrZIf1Q2L%gsr^kLYxtCoQE~&P{hMK-u*8Lf`P{pB0=usX>V)alJx+fgt6-)F`D4C7 zW%y*4uM^qrzp=;ds9YPJlI%51%~i3Hl{cj`C&#*;kVoY2Jn3tf)7QPIo7OS(Bzkw+2{Fvc5c9fK9U$H<~C6KcE4&J79Kp zwi7?DHa^K;wN&({t7qDc3UyUkAG20FU$H$lWqha#bj~Z}w{u>K{5AJ0_B&PmCjXE7 zQRde@m)^g9Ua&v0Uwepo{}BZ|Vj1(vIdD6nKlJ~J`5bT_Yo&#{oAwQlU9Bc}J5uwyokOu;@D9K&u~Rzu=dFsaL?d(rYlsI_4fnBsm+{w{8D z(FGgB4~s3yM7CJ@pC!IiaWL)Mt_Ra%&dUAs;(tz!l#Yuuy<>4G>FeC*1I&}MH@;Jd zO&_&0Mr~T8fu%t5{fjib-Jq)6x_?X#g$EvySC$>_w3I%4R|U12PwpfCSJ5L?(Y0II z+-7HXx=;H*R9?#W>yAHs&(8UtV!!tKIswD)!}|vEBm2KienS4me^h?T``aq7Me3D5 zrjZ*is5p?_Z?$)aWfc0foY@yGKQ)BL)C!xdOzDs*9TS<-O!tn>6%+rFqpdCVnx>(c z(ar)KAMH9JzpLV>zK=_*zp@B`o_hHe+*Ah~HmwL7x8lf1@gvdGV{)G4nYyuA7PbUS z?9B8AZWjOl>1AGbe5BKvZggAcz3zPIZLPn~EnSW4@0LESXe{om%-y*=-j-!4O@U2%|~R?=O~)^yZz#Qx3zK97A1`l;=|P4fRS_p8j0!TTG0|N4KnrTl;2>j1zv#{)8d z`Z=Yxp=oXk+;5Q9KX-zy5V=Id%(Op|HHI)T5b+R(LNg&#nuPrxnshL6VM5( zs-&*7$-A94j3LMB*cS*#>{gu*Mz=AvO@w->>WCKk@!L=c#_*MmbXiJNFF>2>u{u!A29|E$998fV@0<3l*TJ@-!L?!(}neMk9aBddxp)xhi?Fq~-)>j;WDXhs?)& zLw@_XXy_Lb$A33=-7DAx+zp^Av*(5z@ZH#PXrrp0xvQ<(ks_NlCN+!c?FIR8KEcJ& zeVCw{^apd($!(UW4f6O5HnZdB>V(H^Ju093tehOa&6KB8VqcEWy|?9g`XpdCqldTi z@>zbkZ*^KeKVdqqexS45sBy_rtQ=c>iMFQ5O?AgB5{zY~Pz^SO3M=DJ|I zbe081>VEB~GxZYNxP7A-b}!A}S10XU4DUoI(7Phj zU(EhKeOuR#9U*hEoEdT0B?m^MPmb(P$w{z%?fl%cLyrH@S*$W4=5q%5BhNEO@(a=a zbCs^GJngiV)o^@yoRlBA5GMFeqQknIVErtF@JM@EUCu* z;z0-x{(fs|PC(h5)j5p%PSa=yNKXk~5XR5nL_(0i`)`c>s{Qr*b_K#$Or7-IRDg8TxMnQCvq6vMJ(U zBK=N;KP$gQ|BuT}$EL0haf!JcrrZ5f{=#vw7GKQY=Ea;en;uw$gYsPUGH2g&yD^;a zSM?PU@ci<*V(R4I{f%-i2)jN&e4pa(p2_wPw(r)-+jRQr^qf%tDefct%3s{)-xj}< z|ET^$|5kUrlX&NdH}``>jK0rqC)jcMzD)v_hpFt-n)6~~znaR*9cQ-;7w>|6_SY`; z7bZTUot+}bee!Ix6!1CiY`_G;s_vXuu|W|)7%(p{!chB{3$p)^X#FJo5U&6ue*1ar)+B|tkUj6ksrE-3CQ>>`&%5XL4zG{kS zo=u1Kw5Hw)kctn_@7C>lBT8t+zxzEq<-f82;C!r){Bb>|6Cn0?D>^vw{wsg|er0{O z?!?Em)2MibDqs7(;epSNJj7pclH>B(?bP{N<}KfAcGLL&4@CxM6}>9*+$~>Wp02@b zqmw5uk$+onO>j9j*RL!PnHB7&^@3u+s(+ z4hTPfTwdACRXvigeR?`+ z+GXhknE6BdWqza6llO1kmj`IpKNvst{~~|bf2F_0Ch|A&sZ(BB`N8>iBmXi^WiB{v zqu*aVVj9((a~0ubW^kU9sPD6>p#LWR-~2ZHj4zAy-y2I0|7-@J`x`cF-od-0FQ?=l z&L4lkPIn)I3cne67_{m^vTNC&xLLg%2OaL*89W#t$41_=If$Fr+CTSIaqL+(T@yAv zCWo{5GV~MkJ9={1_UgJG4P6>?5lE6#z*opW+Q)#=yI&s14JaSrs7U`iv3M8eT(Kv( zVJg1PooVr%lPu#e?r$Ojj%PD@GKJ=g0Gj{ty|Qin_Hu^aVE)B}n;W+3Q?|Dae{w$_ zU?4w3{;Q~8Ykym8xOJY@tq-H;+oT8CzkH9*ivR}uK(7h*51-5akDWvqZu(T9CWq!) zVsk)Q?CY_kTgLv1NblV_B`lvjj-?)Y=h~85cj55{CTWQb+gr>3W7L=TN&Y-;IkL8n zy)gC%?g(`~EsWhC#s9GH?XsfFc5j{%iM%5c^i2Kv*n03gY;&2N$=))49@~ayr~sYE zLe%H`#(V~v+f6T>UwBBKEhc~`I&l}I)AX-XO)}fttM8q6-F^^}@viypN?yL->F8qE zc1PC&3^rexUah*<@7Ma*ORy7A|7pKl`M$JY>_1I;BmG8zL-U9BvjMlx&zyi(e#L&3 zzj&YgJ4r?!;3SOwJ*`-XWgE-bMdyi(;qoj)8ay>;^aPe%fO8m~j2?OVn&HfQk@k>}xg zeHD1TQ4yz!{dRx08w>t3HbU|O6K?mOJm0h}91ZW^zHcW|v}bB1Z0O9KH~7PU&Tf>< zXP;aVe3v)mld9*uJU-y?Jip$rXdmq_)06ktE$akye^(_7YcXqW>Ypy#dsKX@mt)G$ z`WbK4{>r|%?}YhI)_;qo{JZl@*TSb(@p{!jj{rXFzGTlwpb!Gd{K zzK%M;KJm4eri%Mtk+w+WtVqVxK6|6+sQNul8Ow8eY}8SCzfee3sP2#KyaCxF`=iUp z=B9lkZJAEJce}MapZBD&gnkiw^SxwzZ~;3E&doW}$LhxB%kH0(6Y#oxC;t`wsPzo4 zS35ZB`&&;b=G(mA_1pkJMuZQQ?K|o|x;^aeq5rFWp5=R%^)qVzBl^E}^uO3| zdEXQIH|9N*d)TUve6-wubt~6Grx)Eb`x54>=GI*oonZIcoDCQ)7Kl!ma^Y?0BVBw* zkrrN^3M%Qy`=CLcquS{{o$~1{>=xfE_bd8a`%nIQYj?}5+9Po0#%vqHBih zKA#izvseJ{9}7SaQ1{4_zMSuAxjCt1UgT^bJ~G&`C$m}&@ok5%>>Jk)VxU@J?n^iP z4E^tv&rbO6cGrpXIpzKX{cGGed4AHr)kf>TzUz0?{;K`#K9bU5{4WA@!rqwY;ro9umS5x#na9f$iO_smo!c?_?$iyRN?-VIb5f4ZZ*^WD zgi#cMbI={mqDr2J|H9E%wet7Po5c^(1=wvPpA7H(tfk~}JWO*bb3qS3-Np35^%5_& zcE|UZ>Dkz?EbdnPT5EOfA?@W`RQ@CTYm89+uk|0e9}dr^*iX@)RqHSFCw({hUuS*| z-gkNb<^HtqPiJwbdf@+I9m4@^^{)JJd3aQ{_Cw>%b;Z4VX#?#AaN^}T&&ze(Epk06 zw(`_y$$R2?%|ZTnY%z&JQ%P9|s5z;o0QElXX5<0z7W=kU$-7^QX6~GGV-~>C$-Pg^ zYdJf%3b97|h?|SCG^Q-W1MG5nfBp`XK_p@thl-JNhqwDrk)+Cituwk1OmkJez-Kaj zL_JGw)9zjsLI^K!^0qwFDS=VgbGLouKV)C*kM|V)4VHiTv$DO$d@{OrT$|}A-zwXy z*y=u&P^pIncJ zq`BRO2IZ424}HM#FN>BwH`}l7mv8sTrS_6%EraBDuXZ^=^2YFwod`xwpRqVgH z--!F1f^i1a{>A^(#3z;iYQOxC{0GOA&Ub@LKN;Jtq*-&i6=Plza zIB>v$r|DHSBVRrT4ymJ0iBjG(`!+tXR zqws~!$vKtB+sxOzLAMgj2zJ!YHB~;nAI_6%F#7==^X|yptOR=F1~nH4%bvvbFkjMd z?>70%>D1ofE$Tgf-`N?{RCHXK|H4ib^B=n%Rx`h00XB{Oi~9`u-zxFn#QkaZC-<-R zPx9B>yL~~*+)Vt@(>(tk+O>aeTVkeOhFYw&RM7JK^hxth&aC<(# zoGW~EzMFJi9Ek~ZSRPf8(;KWp#0FRXIws%`uJDfe=X8n7%!{O>7>2ka-Vp@Q26sD# zQ(6qyj9-M$D;K%^X;4GT^8-XAn(FvCu|^$G`tpBj%t!~ z?#7du1kxJQ9(f1$H`cz>T-CltDk8v%Zl1>?yP9~w`@{xcl-ioSTt)5B zyzalmGT8QEJ^;3aN&t^)n-XT_u1m>w^D*p{2cS{S0r_l`IOPZ!rnq9djYO8Fb?y^8 zb5bc7xq86GQ$h-^&uku6U#7dtmJ8IfWGHC+P-?PM*Q%7IW|b zwRrMc1L6nPE;!=eu^g z~tM`mVP53qgLhLR3 zgQ#O|;9`p8?dy_P!Es)%C^sMVhxDcH|x>vE>RX3EMPgQnnRz!2!ai3&Qv zWPi=e;0gMsY5`;*p%%cHqM3 z{jrfa58gT6wr=5@>tDe{w0HLvc`rT#`{rlPh8`woL2r_7ei=XA7WLJ4^Zz0s_c zOi;a-{GmR@du4twHu+trde!DP0_==0lXi;xH1->L|5WuW@-un=O@03q``tQGKcrmekl4PKJZCn97OY^PIovBy;GCtfqTf$fE|?v*P; zV(LG|SJ+Pz+t19k|6wxmKME_~?9SUo{@;rQkgefM;sRh3kv5 z#|9OPsuRi&U~X@UJn^@V&nbdH$WuiH@&YlOX&j%)8MNJ(hw!e=C!x(|xXJz0EcLH9 zl6~RwwDW4L8zjmcW3{JNrPU>PCb50JKELmb@&LnIdDq+L3bN=%UwV*b4{S3~BfuUBLZxi`t73{X?udfi3y;(pxOLxQbel9^q%&V-%6&Wybq zCKmi8f&g6X|^S(#S<5iOPgn721kfH!eidTi+;nQLGtk_>)BPIxZxt(z>pdzy?erOQKIp&9 ze`x;{`w#u^*OC7Z|KGmP&iwIy)T1kZbBsFE=XrT=kBhD9K6f0<;_2<|J}}o%r=F{O zLLzRATntnPq;`EoL)T8>ilc-xJtWoV&imKs>PVar)-AWk(^b8s%H246`h9o zn{?JW4V~A6Nkm$A*0G@*1C;9-P0)T>aN(KKt%rY&|gj zm)Myf&Tq00=KtYho3nLjwDXOT4TSeykytH2Pk=t4_E77cF*lXAJ!k0+^bBB~RU^Bp z!ei}8-90kK+}tVOza#I=RuA)K#HZ~ugEo1^?%^1dty*vF@A5R|!W7RZbN+hL|H{8; z#WWRfM0D-ohOVmn>Z5CY#bk{#{)!!nCvANc`SJYi@?T~8i{nP8|5Wvx^nalKjs0$$ zh;WN|fx-Re{U`bB{OK&m!8+r^4w zJv8UbJb*@lWL_P7LF7m`iww&NJDZ-LkSU75Oku=XHKPEhK&*g`^V#{6kKM%`7zsbL`sC%S7Y(L8b@oxobRDX8Jm8?gn{hi;5=YvmS z-ltmNG;cwr&=k1NcHMF%?`yqZ#eb6j`n-?NU(NfRfZ+uO{$txu^Z!o#6#cj7M3mB4 zhOZI}ib2WW465yOJzkuvq#i0>!q8WV?F3J3ARS$oB)->~*?deD_MJHu-%rH<*v#a4 ze{L_zq|Li!BFE;*$KIRkFXt+rwr{?dhsgrowNd}cpRUtB99lNOcNoI$u6C`vah8?9 zvLc(EmLA=_1+&ZXi2N@yJtF!3_i`rp&0m&5_JoVK=*BctQ@$*s(o z=|q&Fg74=U?eFq>^k4js@zDri)K^T_eQvI}%s}NI-%?kY8QYKEB z6!($#Wc>&BFAA*u@39g8o7i8Xzw*w?ITnlXRZqx!Vu$p8>IQ8M_$5n=Dj#acGbB%p zr5!q+Gqd^D+10jlW^X(u^LXx+NboO$37r&qsIFrabf>${ zDU*N1-Mep$07<_5X_8izZER*wuG~OeAYP;4@JeS`tw`TUSy-J_uuEQd&z!N7VtEdW z^vzplA3B@GeR1YMZb(qYOO6SHz-03rv0>-l+1J<%zcxcpcP*52CqnF|jAnnXrc)pa zD8??|XYHK>nI!r9od1~XZ~EUhOP=rO@GaM;x<9gi@}E3^HNVm)`6>6m)rMRCPX5!( z?{UHZ*>^~P{_Udu56iyRhtIAT9A=;VwtWcuE^nA%X8o(Vv){VGl zDZn9h3!6W|5AX+-V=Vw5Ryo#+o6vz(j*0%RrRC=v%FQFqN|pbhxz%l!M*P1^{yr(= z9<(CeR~x)uXRJ{JzF%L94d6@6O{v`beT;tM;H}1gL5T9oQ!p;iEx+t_V}`=NtYb^EV!rSa6Toqx0R7(RA~k zD}Sgui$Ly>eNoxYk&|0}d+SjDoaekco&ipyW^31wSwv0?{TUuw{cZ3t#q=**p0?fc zs`y#Gb&^dd^RVmeOS3&QsJs{QZ{Aty@=tG??_7n&r?r}vHe=zrz! z1kkJJ=dYN4x^rYNO2GHfhf#dSjHWT+XL$bdJy#h&gXuL{etw|xzts`}xK-+dcZvMf z|5xk3X3SOoFvyXGl&J!T>h}3o2A{Y*tuM$G>6S)UN0$DY$p7eQ{-=YmoET4_XAa5& z8^*%gHD^HA?m6)ex~jh!-uBZorSp|xP@NNf?8MA+aMS{KkKDy$vNhhHzOe%3*U)wS z0nH6*J-sz?-N8W;tn_sD~vn!kmHQB54yfU^#HoV!69T%KD?uTw-+R1q-}kCUk5wyJ;#f z$H!Nwt(dW|`T_@rxQy94ov+l~dsy?v>zzDG7pwp67N`|(7Uq)_GaMRB*8Dwbv)p9m z1XQEthw>VB$J8;%5DVlCc6uXQ)8FIsS(m;J={q5f_~eA0*{G7*b5Osp`dhR`#3u*r z9K9Fj`E=i0*VRh@nC;bl9ZtZM`K^!rrxVbqZ=8Xih^hWJx_@l=9uuF%kIwIl zQxX1U{J@=OH%ndc_+`9kjvn-=yObcUP&wD0_KV+(x=rzXB~9zf)-47=PW-RPySPlH zL{e-(R#1rtM-VSRQ%X}bCUIN zmhY|8M~AhnPzqs)p z@c=MDlRs2{2lJ<1NmlAg?GX8!^1EZ+ITTWbSWlQrFkUwcz>c!ZuKz1y>G3SoO!#sP z8ht$V6i}i=gAM7Y!pk}@^}z3E64!lZe-(DtC3(~gWt0Bo63|(7KvvPKjF}MQZ?XB+ z!(RRNS(vJB8e1d5W&%z?tA>1s>_J=x9jv@tS@eur?^Xchd(H=|;!MdhhL7kxk5$C4 z-EqSUbW);>PPM;#v0M4+Be3+=*qtD@$B8lptY_m9D^Wk$*XonL*llS4==I#R%V)jC zZ|+|^aO-9N_p1Eaf9Ow{A3VGhpySu4RrQsRJY#7AbcaBR%%N8Y;G1OuGPD=vQ;5n< zfYNuRa%vlzxUe%^TZz$iWfli!)<;qY{BGuR{&BW9E>a0RBOX94_tlX(bl&_xyj-(pN$WL<2N{7;yo6-tgnIl(ErALQ{+#^)ql*tDH9y}uTC_{9|!P(u{94{M*sX# zb59e;TnC?DG_zgw-VTzHC#nR~7QFy8eCK>e14PdJquDUxXWz9~2UcYB_Xnc)`h`A} z6JREjeu^JP`oByU#RnpH{=~LGlhv>Q_>s;|$FZ#q5aEg5Xs?_6s~QP~C0{P{Uq00k9nY-Lv?2XySP7#u{RaVlq#UeBh!}?b*qI(z5){7jNKuwxAt8cz)vy zbaRVGU@Fa`MWw}wNneiLnZDaBJ-@pAy_eRWyw}UWnSc9yApVu~zdhrUg~C7?$CxU{9kU`XF5dw%3E{wKe;B{no8tJBfU&A4JC9(a$i z@qQgd{`_SO*xI}vjoot+;0Jm^oQC!TcP=!YrMpE3CakSjmK7jh{;0S=Y#)lmp4Dow z|5*QcUgaC_8Y|^>J+r5PMKGVQs9gJxSPl_n=eO7HKOE7Cvv+u1e41^RS+!vRaX-zh-m2AU7DZ{k4MIEEHR zl71{~W0ePdZ}DeYVq!gV^6qhvF~pn~_6HZ$bJ+{I!7Fk$WOpygE8>V@t;knk49x)G zTCikzfW3W<_c1eANw&Y0e_zY|Y-cAvSM7eQ{H^!#`;8!8l)qvwpnOGS0%ft^wRqp^ zM5oG|yehMP9c<6md`8?qY5HUXCh0dS+*JNh;39w`z;2i#f4*EUi}h_E5$-9L$mkRi zEyx0U_W4{_oW*S>RSyiMH>nGk=S`6mPWEZx;Fg#R2{3q$l{D$wJa=hW~$Je3%{W|i$&VZ3o$N-xG z$m`naBkN;Cs>ipQSeVg&;yWId;9v9y=fGxSRkLBsM?>Xo@&izPP~LhoHI`7t-NWXU832468sLi9T`QZAIetrL@@~`Hv>i>}c zI00h@j2W;HfX;hYho;LSw-@9q-X}a_>!n{ZEr9)hNFLLxh#@4#acb1D#nrN3ECDxM zk@{=%5#ATQ|61a_FC+svA=V!s{E`hF^XUh!)|=SXw%D zUzMNi9b2NZq`Sb(h3;SO{9svY(fR7X^J2~%M%v!dJGp%kNf|0zuv_G>MrGo_YoeX3 zSY;|(gyf3n`Nrxw`@;J>0ra_9s0xKizaEzhl(}EDtTAa%is(~Z4Ly+c;Z~UBkJ~$T zIFWcW^-6uR10EfGk5$Rn`wyuPUohl9#R1MvGr(4X%>aJjQvUcK^8eAX$nAPAjifFTj;#+mP*XR1`o;^C&-%fS2yU)(OD>i&9bGmOj z+bUeJ8)gH`ckUD{Cc25M9}x{)o#jX0SG|5^HZ1=%&wt2$ zd?o))<=*dEe?P7MJ6ZC{W&jYNA^#ooPSrh%VB?>TJp*b{7=w0&m>2QzynETANEc>Y z%*ZYSzM%O~=OExu>Il%Y%@ga0t%lCsknc`7eq5PfEn6`|9nH zKmM?pp8S93!N~qalbso%M*#nzyCKpag1<-7MFh0`Ep;7ca5XVs)K~ zIq^y0h2*NMvnayggYTTbV+xch(0eZ95&0sFL0qCzDSO;$qtAfmc7g`K==_L|S$(!) zbO3gL>jQIjr)dtW@*i41bidWURog$G1&OMw9UY@%<;{pN>9qg#bX3af6)(b!8o3c!)Ag z6(cr4&j3b(yg-zvPC@I7^vV6qdca#X!N$B$(Zz0gyzD=FIiOa}H{b^4qVYT`SoS?0 zncsU?{+9U>MG?Ix(;^Rhd2hXdB z@1_>M5kRB7>HkUlJtqGf7Jvts5CA(y>;Gkb>en`qsQTD3qeqndN5^pc#kb8f-~*$N z!JIG|V`tl7w@rM}5p(UHD`<+w+1b)^jadGNhYz}CcJKbf?5?r@mnJ^IAGaG$Z*a8r&*mJ609r5BhfemSjkQmE`RC%sytL}D`Q1Dcriv*+qRU-! ze&JdA%i)z6X=Xm?Y7@WcV?H9Xx6ib@!S~O*f^dj%Og9y|=qWHmMnw-60234=%a*Z# z%l;XYe24Tm?dV#SzZk&xI((*_Z}Rp-%f};={N0-@!maWOwBL;f7=3^;fL8WJajlYw z;~pO^H)+y_2peXd0BBJs!~2a+u2$WU4fuTOo=&~bU96KwQH!Cca(eT!PYX-%n_+)! zHrqUY;I7%mcE@aBydw7BKA!jI85kfYJee(aOupAepd(1UrxJjZX^%E^$4d?{L!OeX z_4353Uk!R;E-bbKW`jvM@>e%q@FF2d0UiiX$p*BayWRhzqqp2XRsG@hE1z)!?w-io z^}lCicZZdnZ^(6${|(m@?EdKqyh%L3mASTh zf%Oc+^Jg2cY$?Iw< zgSGS9)tBJQ&W$M|ru(w}Z;LLQ2>-!&B(ns~3T27R2);C3i$BTcz@O#apP0X44j$>l zw$90K9gvfSc>vq#{GUz%mlY2Fhq$jfpMJj9bw;x2Q~+ zVs6(oCR2CJd4*l0!5@evG7+|Ofb|uhp6@RRW9fQ%w>QeP)94?#v!5#eR?(`?Z-(>_sSm9mp8sAOcAVX9{)pKXxh`AAeou)1 zB=`0AD*5UEa({aYy8mY&05}EmxypZR{DYmkwestIiu&E|mz!+iLAjL=pf1!ok)yM@ z+%{EC6>}#*j_d@$y?6nAJ;%j%n#piMYCu=U^Zz`j;F{?Qc-ZW2sRXMK=m-`4v%Ti) z+9rq>fDa)1Qm2LAnzSb$6A`Eiy1m9Q!l2oSR zu}gWIDXL;|TCmBF=k}!~Z~dw1WSlMg95CHYgE3*+MAL(FMqigxN%o}g+?53`<0pHl z0Nrm#p9+6tz9D~dZ>-ly?XNfNHv84u)fab~T|d8f=Z*4TFY)QSO1ZN4sF-h({}lh# z{*yoZzh&M}5kPZa2BPVmZPGAX+7A9 zC^xVST-q+@f=8DD%=Wjs#RtIX)$6-cK%cOxkjzjfv3EQVj)HApN9DbID1PapoYv2! zhv}ZlTei&H-!qo6N#_7|3f}(tO)*4c3~0|dbF8Dd54%fOM?3lPdBy~%RlOJK_jc!; zk7v%Flzqm399y|%?05BllK6D58Nuoa=>7Ko`~Rx)ZwB~)ydxE$Zv3RNuA3A5l)U4= ziM&a^5r8PrDd@~dyK3Zu<&XDt>lGW^cBr>K2{@UMR#DpFvH3sFF6JAh1NeH;{L5oE z&EoxN5TWbL*5Xx9j`ck@I%OuSZChmn;suxw#^b&vvcmk}0dV~w0Jc+<4FqX<-qP2z zKjmbRK8dqQFulry?`lVbEVZ5D&e_;Y-${bU#`>Ehc3Nb9NstSh4!)F=VV;yt=5hiX zBHkDaH@HTyFz=yV6yoIWZfJ#Hwmh8 z2kPD1DkG9jcjGTX?g%}9uD}~k+g|?$0jHEO04Rck-vGuCOMyy zpLN?FtXD>0`jR>x`SU7x-Y<{Ns@-Gzm7OW#gRs}~w=zTOaug?9Zf4{5h+oZJec|$i zcIMN*`8x+L&sz5Z+Uj=1^85G3-rEam%8aQ&_Q4;WpY0b9B=2wb)LJ+UpgWt%UN7HP zk!C0NmfJT){6qeO(_Ng$6#qNttkA%_EwjJMe(ZD@Ct=9AHh+9KP@n1hT9qQpQwA^) z;BFln^&Hs>1o+|c0N(Y`e_urblLhDmLw>-mGX%T8h<$xoqHA~vw#Cci#m`U9ckNjI z@6WE5INJOz^8A zcbmtkWO*jm@5pT+9FUl)7u#e5u_PN~vYY?ioC7O1j#rbvIi)t|gU?@iBZB;a%v z0T>w|4{-lT7zdmmFHeP^YzH7_%HRA>{$#InOw7=&G2IwiZm=$KQ>5D4)ZP||0xKx4 zcwTJvn>TzUe*QY?20byB;H2n}u76Y34v7`DU9lTt&kw?KVSK|E^Z%)NOqPyQ4)fQW z1K}-ucjt}LF!!Tfg6u-sp=1xsZJ(!lyO^xiz4m6Z3Rqe@EbWN7!XgZl$n1V_E<58; zXNt*BuaERGTzNA1pv{l_=l)Ccx{Wl(At#>SS>$v05mR&}e1L;kJok^I=+&5?7K4bPuAEZZhJd_bh`M*3xq zmJI+2s{ilq0JSWgO_=1JXV%F!Dxl&@lRJ+NwN0g0I9v zd@J+SH&4Il#W@-FE35gN8w#n#(Zkb$^QkU56YUJZ-5A*cHX<}=@V&&eDqLlGcgt9( z_I9>y_sg0m%+{FvWp*`W7mLBoK1#eIh$G$^h66SSw5N66(eG=cFO0{Y#`g03b*}a5 zXdu%cbs*NOO{ZYrtkeke!O&;4iU{ zk6P+KUJCiUjmc5v0+2ub@AiPsv99+QUzbth%CrJvs-$YXT0HGJBtG7Z-jfnb9T-Wf z6YQIq;Fq)OC&&JE@~rPBGBiV|^SW>YRe4o*Oyhd<(sH~vkG?W`uIAr4&Z5Dx`$(<3 z=th5^9s|__lj;3Lr<3|cvq@DU92EIxd1l3zVT`o?m65eMXuCxE{Jrn$Pehk$zc6uG z{?x@Meunc+llvVM=f{_oa^sdKcmiDiQI0PS~x0hSBZsJs1gvT`-oS4W57tnXb$ z)zH8*V&!dte|_xkX~_kyoz9TUgRxwk^C1sBIA^B2BKC;fDt)_M4$o;hFIMFp8!fZ^ z4W@(nOP|&2%o>LdVg2lqGX|wsO{T>#0xW<0Wfn!H5*~wgAho8^F%k5%rr=W1?1V$8kKU4?H5Nr?~Gj81GC{z4&o!be^2cF$HNHxWF9lMRlf*y<-nW|df8di zcpSKZ2PZT8M4op<=4zc+{4|3oBQf8gQP}DOxLZQi=wY zRRk5pE1~aCSzZ;JI?8^NRTSH$E5qVSB)|?mpq=)Ax_Re4qsuHRBcRMTA`Qgi$POffd_mFCFU(uvY5@Uc}S2Z3y zU@@MmzNyCeAWw-#h;ih~K3#i0aKy^vYkm(al`--7V{$|LH}&t1tgo5>Wci2mCyaOX z1PobD@~@m%&%%)NkpF5l;AfNkM+Lwv01=&tzPyfpI$fZh{A1rJ7Qm;@97kz1wPKzw#gGKHD$1o9w<= ze>L8()SZJN|DK7}@5X&!#|ao;`?IUsf7ZVops@+qCb-pE0Aun8{dWRW@?(|HD`o;X(vuz@8bD2L7{kPli zlEeU)1vC6)Y=BLycE=nZSvmpjKdrZJbWHkpE2-ve%)Ldz}nCFItTg zK~WDWQAq*=Hq-m5iA!_gbR7?cJeu)%ua~ zSN=4=v-{}2+pqzZ-&N&5orfWP5+5gE({nJ%fAs%M`Y!{l1!yM-n_djJsHF_H^M zF+cm?lfpLe`%Zu)-8!wcWVYAsPD4@z|%O{`Z6VX0!XJHuU~T z%|t=_xGxGKd}-T##k z+fd*HhyqBS{CR(qiPeHx08wIjAu35G35X**PeNZo`{Z_9T1s}*JTWmhbRMVnb+OE! zO()phQxCEspxfefi+?1wmkT2Q%W@LVjr7k-z37{H?A_5=vnP+*@U1Wh&R_cc{bDt2 z@iaqF1xG}=cTR>5AjmHCpmWCdO6;Ix#1t4cqUYu$?3{a|_-Q`3os3R5dtYgrkl-vh z7lXqy_HmXavFi1wp1l!u7wZ{5zm|WXJ(J~MjrIrDJITA6KJ>jZuZEMld$;Cds?LB1D)Vob(ydz?Vl#ro^8V-PPM;GP8cs4@c^kBnJ3okigI6I~~LqWRBQ} zf)dc_{Wg3k6UDa*O28hxbLkxpTH1P>e-F*iF?-Z>-@cT;P88n2zADutm;rtUex6RO z!OL^v`+1!1&c2cnv1%mm^R<@OUA$AgH${BS??$w@s`vHz!}AaAujOa^>jj3@zpAYt zdOpQ>>*YUrf}RrdfAxY-RfduO-*d?VtNiVPDVo<=AN^l!hXrg;2x*gWbshH08n6V# zIX%Z2qBVXSo+y$u-_E|qH$~UWj5bH&snLB`9dnQmP$87LV+vtl9heC4U!v<*#si!m zjXyWm{D|oLskt4V^I@j9ng3tT{D3o--T_69#{MEmA7t=D`S$Ld$_YxU)PQ^sG7FPc;D zUk@;}erWq`ayEJ%`bhPA{6Egc_}b^(A^QEWXt=+=Y5~SXk!ApDIwPhR0gP9Qrxw+E33QwwsR7QX$1Vt-PL|3F9xjm!3 zAB^6AJtyGfk-V5dwSZNyW!rpdyC+N;HSzyBvF9+`b_vMe)rfFE+k*tn(N$dR_yw`~ zGLTL<7+hz0ZEvnbyB;flF#vn-q}STXh>Gk^x<7KgX8G0p5$mtVc-j2QfAD!LRXV@w z`-l9;*WnE+@!{pi*P-{T^1n^;(B0!R56UO1|M!j{&qMwL18N4y0+4^B`tB#H<>&P) ze>*A0T>23;c1r>I*NX5W+eX^>cF&GR+n2mgqy!hQcNr7M)d<>GH|8O|B$|3uWPWmV z{?g=r$L8I?XPLpN&D%Sn9@MSTdOCN??{AlLR6Y-Em}f7W7uyrxw>h6o-ZU2s0Q}@m z=xhHq>6h(Xj=K}k__qG6dst*cG#@V(GipTpuzmR6ko}1FDt~%6MSHFA48OnL+pk7@ z&O~*1$iH97xsUf847p$R>#bui92Bp0{nZpPs3%*%|$sSbG%(HU(m!;#lrzCuXQU^9MUc{wGE2 zKN~jCPZRGS5Zf)gm#y>UYQSo~t;d=`V!sZKHEF90c4CjmI->;(wNqf8naA3-0=>b; zchX09J8ur*^ckApHX7p3=tB2zZ@H8_J{(;CuYQN}bUNF}@M`mG`-k?Ah;QV3!{={$ z1}f*@FaM$Y#Q~?Y;7kn7pU%r|^Ii{}%@Z*_d?U_S6sQOQ>6-;718iqsEj#)6cdfL;q#DbiZf7X?9gcFO+KDl;uzHUVZ+k?9%%w-kLoBRdao^A9??(?Qi@? zru7X&rjteJ`RL=7{p1NI%W$WBn(gxb@48X`wf}>91^ zlqXdroN|?26F^k&aC+3?oPlCZnAWXdl!>M;<0kKD_l;efov=@*YuWnl%Wd`PoaDwm zPJ?@m^;d0vJ$^T_4veQc9xa=+zMl84d1TXQpTsw{0*wQQ)cyRbdOc3TB=>O+ngRBI zKNbDskKAS{|5XgQvw&fg^uHLulmk}w;&}C3Qz@InglT! z-Hv<;PeQx&9qEDQJ!BR=W%P5qSaJI|&FIJCRZHrO9oV71GK_MyJJlBq6U>Whed}RNz<31c9^)dTHNY0uxYhxmGqYRBQn6TwooQONUQV%B>{*bH`ir6vXRdkv zl(#gB72Q1`w{o4kCri9X^DFzJJk9i4MOW7u@+bRA`dEBAFJM6Qe_*?-dOu`exsU%V zxnTz?_hI43*UG=1V3Ks7JN@t3ga2C#FyVg?YA3k7uJZaS|9WcoMQo3bEZ&<~{0-4J zGb?Hxhp%)3R^?v{T12rq6K$&(?e&HeRQ@;y{KkkSNUn7#RW|5*d6aBOaRt0x7_!Qn zS&Ps)ab>Yc6B< z{;dLt0VfoI4poa>!6ya3`Eq9Uou63Rb^=nYe;n46O&IpSUwfL02f5lq{(i&Nu5Hl| zL+2Vz`U%bK87SfcaVzT9IDucUKW{f3TjjNRC|(`uLL5xm?(4xYYpZx9 zf4O*WyI$J!-227#YxP~TA^l>zu4`?6<-aO@(LO6M;rm1O#rm9sA@P3gBSW&|D`{5m z`~M{W%6V1(!~fSC_<84lH!Ja;I|sJXCx0>EkiRR&Vovh^+vwZw+2V4|U@|X_{AEI} zTjQi+gZyS=f}Vp$hm8V;H>pJGFPa0;+lS_LXqwEcu_S4X_>*2k6MF{aYdkq?P|M^u ziq}3O8r}_yJb)9|{MCD?71}nR15fwKOPhhEH9iU7&xWaUj_jWP^ZCQ;Pd;C!&*qP8 zuks(V-&Fd=_E+UUG=E5be69TbblCkN{Z*OwB=|`szxuz&Ch4!u;CWiTvxgYYY0_lt+<2uU=iQbuKHi9)P~}XQs@M zKQyrsS!1$F)_Z7GVMv{3vmR46L5f&x5J!9ybp{b<oTqwU+4uiSy;|SPuBkNa0W z0rbD`lkQjkodY%OKW3-Mjm50|S%9IXWF#?6YsfAeLJK-#e0qTROUPd^|ccXF`XK( zS7!Cv&7FD%YC)SfIZfo;{7wW`rp)AXp@5A5{YE7}~`srPq&o4b^!G~{0k(5T7vqW|@tG^XB~j?y=M4Yo+sde#nj zyw`JFRcFAhvUy~rG%-ndhen1I!>EzPc(r~vwtz>N)uH#`dg_ zSLivc|M36F_!_r($7~n>KSlelk7#exJ5Q=Ydf``Uf~uj>BoR^q=Q|4ED2pMX^f z59zNW0jodXGzOStzUuoM|4k7fJ4@<}s_j>u-lRyhaYp${@%Sy5st~oim zk$!KzlI%mzM>SD?RnIW}M%kqLqKmgQtE^lrnO5Y;*Q@kw_*%=&UX1tc{;)Z#@+WQQ z#my&jLZ(yF->>u=^YQkp_}-}L)_0P8Z>)Z^yH#|5lKynBr(wu`M0}Oy|1AF@>+y_D zt-_GMukG{g31IujS?GQ5x7MHMT15aV|6;(C-F9^vsTZ@+C921Jpf%)=H#p=tG@t$q zJ!)1rCrEAdI{&CF!x1EuRJaPIuYo*}Xxu}b+( zzW2$*G^3vD+2gN|PAyqptNE4vi1zyZDj^+WrI)cs~pO8wPp_E!$Tr~#+~4f)elHgnP+`p6P&y_A!v_`bQ{IabY` zDlZbH^ObC^bG6wSU^Qwt#_4fhinC6b+c;I@`Qz0pC0SX|%8!ZMG{?PTxeZ+Q^?QQcH z?iKyCwSD#0Q?OqCSC#(I`0*33u&MlgqEWA&EC3Hs>)R_;?>H>Lu*hbWij3?JcxaH8 zM@6R9*6KJ}RRb%9aW+{0dfDn5S<-r6-SdAX>kRe`RsPOQ{fbyqWWN69w>qC_vIbtQ zXV2g7EA)MdE&ExzWpB#wn%&ptyE>DlbZ>b6_42RxU!U_8=_U0cYZjnZp%#940@5ER zVSJ_icgr=t`#SlP`l|dZ)gj|4(yJAi9>4X91!?@}vylE(<==Urz4Fru81f&s`~DkI zFg2sCVns~|Y3%Kplg83=dfe;piVyjdL%nw;SLyeD|EIl`9?AQ+ zw0O2>i5*5}PTJ(}YkSjs7RLDyOPE;Izp1uRd*3?lr1!1m)Zh1-3~yfxFf@Ovo$vrd z`os2@&qvOa{D;)p0{172P;Kx3mH3$M=tT7K%B}W($gt8T+4YZ9>Z`t=O<1)9_Ww5j zYvgxT`L9lQXa?R`xmj{ggZO#NoRgx4PJpbjb)(wVX6QxIht{zV{YMtn!Epjee6j$; zviAh^oQ(g6_aE}Dw0i=`za8qe4LuKKHS}GmMgC3Rd->YMb6dw9@?Vwxs>L7jcm0O% zCw*36ME)c88(%B=%6$0$=>$w(f1CpdLeEEk=Ko2;)4fTe9_u-n9;@{qavxv)+)sX= z{GIv{0rdLefziB?Ll^sLuH1^vgCa54th|BBzgbtkh6uZee5GB>JMLE-YvEVD0Lc$~ z;A_uBwZ7UvS%~#A9_T!ol6X&o|BGw1t|3N&SF#o4?%AW7TXb)d{}k_4=B=9y*;nr5 zKcqj9AJSi+0}k0wUrBse{c5`JakQqLK9jtsd#fj6lK55cP4@q)@^_N!{Tl&R`m6HC z0-5ajKv5>?%b*{WXZ+Rj+~^*y{70%6rn}O1s}zx~*&0t9iY8jz+J}ihH8{W_u&wtJSXs81i40{p9`0?Dw64q5YNm zu=~EN^rth>d-6SJ3_EPnd$J|}Rf{m${Yrk4@i-~J??m*mN`F=U?Ex4y;9k!u17F3_ znll#>R4=>#@|unAaWY~v+#kwuQrco^dwAD6|2{ivjMg0exv%{-))32EdMwG ztC~M_dR6k3f6qld#L)OM9wvOZRN9`8{CR=;|8Z}M0@mApavt)Zw0=mx z^1s_A@*fiIiJ0_%oRIaB9?xF?=<0Vp*Ny#a>-wp6fRqF18nMCy)hmVu2WSnps42-l zC~^`5$Qj0o84+JkK%+sAvAK=I#D*gmtNu6dn=HF`(1>n20X+$x(^KN->rX&!{d59` z+@Rz)t)=`S;UTlJ%a=p=-|@1JL>+W-zmU_mId{6j<&s zkHNCj>3Vx-fG_Aa<8eY*6pz(T3_C!+mGOw<$2qKT@U!11`Ege2@rT{7?AOcOcdPR5 z{Ym;G=WkBDYX67yeXZT!^#5T2rWjzo^&eV4$$wMp?;QB)@Do>+KOaE8{+~2@NPP8? zA^%MyLErma_4||d^Zo2UO&;=Z?A;uY9k0a4324=n{N;f(YeWFm=*pjL8}ZkVv)|L1 zn4aT#lMSep8^3wgo{Gi<<7d_Uktft+c&zVUOgHJgN2~vn#K$LYP1lw1$ngEm{Kg5O z_2UFo|ECOKeH2hHAod&5C(Usi*1zBTZqj|y?unS9!eQZuo=;yVSvL~!ePe(r66hHi z+TX90e^I|$x_;{60V+w?kfpBl>i@}It&9hl7+12u5l2^2jp%9zM17Su`%k_jK4248 zrSJYY0RzLM&m#(O2AT~F>tEY|sV)xfqdgPe{i^rx_ml2B1GV)&dwE=yzpLOL=}#wM zNWNdM8VR!b!|qr9lk~^EN#f%SRMwN!EB#^N$Ja^X>+enSADGa1WXOAbzv61* z@qpNdIpQzYsKDbo_S?Gu=Y4&PG6TpP_m1q$1}TBjZ?2oO4t3Qq{Pln+fVis#I|qn- z)FAv@n*k>NeYR^raXyh=i@4|Yc|F{FrR_4n{ViiWZ24eEoOhQC^aywBJBNRoXY|`U zI62VN{dK+0c_3zim=Csw)>TPQ zC;-^Q0@lSn>!SSowXYA}!-d*Fnt1|dh-O4M%#aEG?kBpyBL2C5k*eP-01a6 zE;V&TKR$r!Nvgo8y0>PRSOpuA4^NQ76Uv2<^iyNj&H}Rl?>1M6H9U|Na?4iO^s??n z_*{3RysjfmUa&fM7%)yi4=^5ZKYR0_8ld|39PG;nXlMB63iJZMZ<&L!`M$O$xo@8@R`{VW_!{EKuM;MN&`uXoWJ z{-X!_KzWcKG;Vk-2Q(^hOrQs390>Y3znKTx@dwRjf5F@Basx< 0: + ext_list.append((filename, n)) + hdulist.close() + + return io.ImageCollection(ext_list, load_func=FITSFactory, + conserve_memory=conserve_memory) + + + imread(fname, as_grey=as_grey, dtype=dtype, + plugin=plugin, **plugin_args) + + +def FITSFactory(image_ext): + """Load an image extension from a FITS file and return a NumPy array + + Parameters + ---------- + + image_ext : tuple + FITS extension to load, in the format (filename, ext_num). + The (extname, extver) format is unsupported, since this function + is not called directly by the user and imread_collection() does + the work of figuring out which extensions need loading. + + """ + + # Expect a length-2 tuple with a filename as the first element: + if not isinstance(image_ext, tuple): + raise TypeError("Expected a tuple") + + if len(image_ext) != 2: + raise ValueError("Expected a tuple of length 2") + + filename = image_ext[0] + extnum = image_ext[1] + + if type(filename) is not str or type(extnum) is not int: + raise ValueError("Expected a (filename, extension) tuple") + + hdulist = pyfits.open(filename) + + data = hdulist[extnum].data + + hdulist.close() + + if data is None: + raise RuntimeError("Extension %d of %s has no data" % + (extnum, filename)) + + return data + diff --git a/scikits/image/io/tests/test_fits.py b/scikits/image/io/tests/test_fits.py new file mode 100644 index 00000000..2cfeae47 --- /dev/null +++ b/scikits/image/io/tests/test_fits.py @@ -0,0 +1,71 @@ +import os.path +import numpy as np +import scikits.image.io as io +from scikits.image import data_dir +import scikits.image.io._plugins.fits_plugin as fplug + + +import_success = True + +try: + import pyfits +except ImportError: + import_success = False + + +def test_fits_plugin_import(): + # Make sure we get an import exception if PyFITS isn't there + # (not sure how useful this is, but it ensures there isn't some other + # error when trying to load the plugin) + try: + io.use_plugin('fits') + except ImportError: + assert import_success == False + else: + assert import_success == True + +def test_imread_MEF(): + if import_success: + io.use_plugin('fits') + testfile = os.path.join(data_dir, 'multi.fits') + img = io.imread(testfile) + assert np.all(img==pyfits.getdata(testfile, 1)) + +def test_imread_simple(): + if import_success: + io.use_plugin('fits') + testfile = os.path.join(data_dir, 'simple.fits') + img = io.imread(testfile) + assert np.all(img==pyfits.getdata(testfile, 0)) + +def test_imread_collection_single_MEF(): + if import_success: + io.use_plugin('fits') + testfile = os.path.join(data_dir, 'multi.fits') + ic1 = io.imread_collection(testfile) + ic2 = io.ImageCollection([(testfile, 1), (testfile, 2), (testfile, 3)], + load_func=fplug.FITSFactory) + assert _same_ImageCollection(ic1, ic2) + +def test_imread_collection_MEF_and_simple(): + if import_success: + io.use_plugin('fits') + testfile1 = os.path.join(data_dir, 'multi.fits') + testfile2 = os.path.join(data_dir, 'simple.fits') + ic1 = io.imread_collection([testfile1, testfile2]) + ic2 = io.ImageCollection([(testfile1, 1), (testfile1, 2), + (testfile1, 3), (testfile2, 0)], + load_func=fplug.FITSFactory) + assert _same_ImageCollection(ic1, ic2) + +def _same_ImageCollection(collection1, collection2): + """Ancillary function to compare two ImageCollection objects, checking + that their constituent arrays are equal. + """ + if len(collection1) != len(collection2): + return False + for ext1, ext2 in zip(collection1, collection2): + if not np.all(ext1 == ext2): + return False + return True +