it is extraordinary but true that the galaxy of the local group whose global stellar populations are least well observed is our own milky way .
deeply immersed within the optically opaque , dusty lanes of our galaxy s spiral arms , astronomers have been frustrated for centuries in their attempts to map the stellar populations of the milky way .
a complete census of our galaxy for every member of even one class of star would have seemed like an impossible goal even a decade ago .
major advances in instrumentation are transforming this daunting task from near - impossibility to increasingly likely .
wide - field , high - resolution , sensitive surveys , particularly in the near - infrared and x - ray parts of the electromagnetic spectrum ( where the galaxy is relatively transparent ) , are key to locating and characterizing all members of one or more classes of stellar object .
the goal of the project described in this paper is to detect and spectrographically characterize at least 90% of the wolf - rayet stars in the milky way within ten years
.    a set of well - defined tests of stellar evolution theory will follow from detections of complete samples of wolf - rayet and other related stars .
for example , the radial abundance gradient across our galaxy and the increase of the wr / o number - ratio with increasing z suggests that many more wr stars will be found in the inner parts of the milky way than in the outer regions .
in addition , our previous hst survey of the hii regions in the sciii galaxy ngc 2403 @xcite suggested that the distribution of wolf - rayet and red supergiant stars ( rsg ) is a sensitive diagnostic of the recent star - forming history of these large complexes : young cores of o and wr stars are surrounded by older halos containing rsg .
theory predicts that the number - ratio wr / o increases with increasing metallicity ; thus , relatively fewer wr stars form at lower z. we will also be able to determine if superclusters , dominated by wolf - rayet stars , are common in the milky way .
finally , we note that wr stars are predicted to end their lives as supernovae , and in rare cases as gamma ray bursts .
the wolf - rayet stars in the milky way may be abundant enough for one to erupt as a type ib or ic supernova within a few generations .
this comes from the assumption that the mw contains @xmath06000 wr stars , each lasting @xmath1 yrs .
the clear identification of a wr star as the progenitor of one of these eruptions would be a dramatic confirmation of a key prediction of stellar evolution theory .
wolf - rayet ( wr ) stars are massive ( with initial masses greater than @xmath2 at @xmath3 ) stars with strong winds ( @xmath4 ) displaying the heavier elements created by what are normally internal nuclear processes .
distinctive spectra with strong , broad emission lines of helium , and either nitrogen ( wn ) or carbon ( wc ) are the defining observational characteristics of wr stars . as they have relatively short lifetimes ( about @xmath5 years ) , wr stars are excellent tracers of star formation , and they are also believed to be type ib or ic supernova progenitors , because they have removed their outer h - rich layers ( wn ) or even he - rich layers ( wc / o ) .
about 300 wr stars have been previously identified in the milky way @xcite , with distribution models predicting @xmath01000 - 6500 total expected ( @xcite ; @xcite ) .
optical narrow band surveys have been severely limited by interstellar extinction @xcite , and so the natural solution is to turn to the near infrared .
emission line magnitudes of 40 known wr stars are presented in appendix a. a model of the milky way , predicting the numbers and distributions of wr stars visible in the k band is presented in appendix b.    in an initial attempt at a narrowband near - infrared survey @xcite had limited success , while @xcite were somewhat more successful with their color - based selection of objects from the 2mass @xcite and glimpse @xcite surveys .    utilizing a new narrow - band survey , described in @xmath6 [ observations ] along with the candidate selection methods ,
we have found 41 new wr stars : 15 wn and 26 wc .
spectrographic follow up and data reduction is described in @xmath6 [ specobsred ] .
our resulting k - band spectra , line measurements , and subtype classification are presented in @xmath6 [ results ] .
distances to the new wr stars are calculated and their distribution within the galaxy considered in @xmath6 [ distribution ] .
our conclusions are summarized in @xmath6 [ conclusion ] .
the most reliable optical technique to detect individual wolf - rayet stars in crowded fields consists of subtracting a normalized continuum image from an image obtained with a narrowband filter centered on the heii 4686  line .
this works well for individual wr stars with equivalent widths of heii(4686 ) @xmath7   or larger , and even for dense , unresolved clusters that include a very small fraction of wr stars ( @xcite ) . unfortunately
, dust extinction makes this technique infeasible for the large majority of galactic wr stars . only in the near - infrared
can we hope to detect the wr stars farther than about 5 kpc .
a near - infrared survey of the plane of the galaxy was carried out under the umbrella of the smarts consortium ( see http://www.astro.yale.edu/smarts/ ) .
the imaging data were taken over approximately 200 nights in 2005 - 2006 on the cerro tololo inter - american observatory ( ctio ) 1.5 meter telescope , using the universit de montral s cpapir camera .
images are 35@xmath8 on a side , with a plate scale of 1.03@xmath9 per pixel , and cover 1@xmath10 above and below the galactic plane from galactic longitude _ _ l__=-90@xmath10 to _ _ l__=60@xmath10 .
each of the 1200 fields was imaged in the j and k bands , and through a selection of narrow bands designed for identifying wr stars .    motivated by the wr infrared spectra in figer et al . (
1997 ) , we purchased a filter set ( table [ filters ] ) which targeted the emission line features he i 2.062 @xmath11 m , c iv 2.081 @xmath11 m , h i br@xmath12 2.169 @xmath11 m , and he ii 2.192 @xmath11 m .
in addition , two narrow - band continuum filters were used which were selected to be relatively devoid of emission lines , one blueward ( 2.033 @xmath11 m ) and one redward ( 2.255 @xmath11 m ) of the emission line filters .
these were then used to linearly interpolate a continuum magnitude at each of the emission line bands , so we could calculate the difference in measured and interpolated magnitudes , @xmath13 , indicative of emission or absorption in that band . in this paper , we are working from a catalogue of the calculated @xmath13 values for all of the stars contained in the survey area which had been processed by the time of our spectroscopic follow - up .
this consisted of about 75% of the survey area , and in general we excluded those areas most crowded with stars , in particular the approximately eight degrees in longitude closest to the galactic center .
the more than @xmath14 science and dome flat images of the survey require a customized , streamlined pipeline to reduce this large amount of data .
the pipeline was constructed ( by jg and dz ) in idl and uses the 2mass catalogue extensively as a reference for both astrometry and photometry .
each of the @xmath15 fields has been imaged in each of the @xmath16 filters at @xmath17 dither positions with separations of @xmath18 .
dome flats were used to flatten each of the images for sensitivity and chip defects .
these @xmath17 dither positions were median combined together without shifting to remove most of the stellar sources and to create a sky image .
the sky image was subtracted from each image .
sources for each field / filter were matched to 2mass to determine the world coordinate system ( wcs ) .
once a wcs was fit , another iterative process was performed to minimize the residuals to 2mass and determine the best geometric distortion solution . the idl procedure warp_tri applied the geometric distortion solution .
the idl procedure hastom aligned all the images taken in a dither series .
these images were combined to create a final deep exposure .
sources were identified as wr candidates through emission in the narrowband filters .
this made it necessary to have not only the magnitude of a source in the narrowband filter , but also the magnitude of the continuum at that wavelength .
as a result , each narrowband filter had to be examined concurrently with the cont1 and cont2 filters .
sources on the final deep exposures were detected using the idl procedure find .
aperture photometry was then carried out using the idl procedure aper with a 2-pixel aperture and a sky annulus from 10 to 20 pixels .
objects matched to the 2mass catalog ( at least 100 in each field ) were used to determine a zero point for each filter in each deep image . an object with a flat spectrum through our filters
has the same magnitude in all 6 filters .
sources were considered matched across filters if their positions were consistent within 0.5 pixel .
these matched sources with 2mass - calibrated magnitudes were then used to construct emission - magnitude diagrams ( emds ) .
an emd was constructed for each narrowband filter with the wavelength - interpolated continuum magnitude vs continuum - subtracted narrowband magnitude .
the stars scatter around a continuum - subtracted narrowband magnitude of zero . in each magnitude bin
we calculate the standard deviation sigma of the continuum - subtracted narrowband magnitude .
objects that are 5@xmath19 or more in the negative direction from the locus of stars are considered candidates .
we determined the offset to convert instrumental magnitudes to apparent magnitudes using the 2mass ks band catalogue .
the images were divided into an 8x8 grid , with each of the 64 areas having an individually determined offset , to compensate for an observed color dependence ( probably due to variable reddening ) across the field .
the idl procedure aper determined the magnitude for the source at the coordinates given by 2mass .
once the offsets to convert to apparent magnitude were calculated , find , an idl procedure , identified sources in an image that were a given deviation above the background .
aper found the instrumental magnitude for the sources , which were then converted to apparent magnitudes .
the cont1 filter was taken as a reference image and hastrom aligned all filters from a field .
the sources from cont1 were then matched to the sources from the cont2 image and sources within 0.5 pixels were kept as matches .
the matched sources were then compared to the sources of a narrowband filter .
the result was a list of sources found within 0.5 pixels of each other in both of the continuum filters and the narrowband filter .
a linear fit was found between the cont1 and cont2 filter magnitudes of each source .
then , using the central wavelength for the narrowband filter , an interpolated continuum value at the narrowband wavelength was determined .
the magnitude of the source in the narrowband filter was then subtracted from the interpolated continuum value , giving the negative emission magnitude for the source in that narrowband filter . the emission magnitude for a source
was also estimated by subtracting the cont1 and the cont2 magnitude from the narrowband magnitude , resulting in a total of three estimates for the emission magnitude of a source .
an emd of continuum magnitude vs. emission magnitude was created for each narrowband filter and the standard deviation was determined for the sources in bins of 1 magnitude .
sources that had emission magnitudes of @xmath20 or greater from the center of the emd were marked as wr candidates .
1 arcmin @xmath21 1 arcmin finder charts were produced for each wr candidate , showing the candidate in all filters of the survey and also in xdss red and xdss infrared images .
these are presented in appendix c. candidates with emission magnitudes that were similar to those of the known wr stars covered in the survey were selected for spectrographic follow up .
this refined list of candidates was blinked by eye to remove any candidates that did not resemble stars .      in this initial , exploratory phase of the survey we used two techniques when selecting targets for spectroscopic follow - up .
we began by selecting targets with such powerful emission lines that the star appeared brighter in the narrow - band images when compared to the continuum images even when examined by eye .
this corresponds to a minimum 0.5 to 1 magnitude difference in brightness between the narrow- and broad - band images .
we initially selected candidates displaying a brightening of at least 0.5 magnitudes in at least one narrow band image relative to the continuum .
this resulted in the detection of 34 new planetary nebulae ( which will be reported elsewhere ) whose very strong , sharp emission lines , he i 2.058 @xmath11 m and br@xmath12 2.166 @xmath11 m , fell within our he i and br@xmath12 narrow - band imaged fields .
a few of the planetary nebulae were slightly resolved on the br@xmath12 narrow - band images .
no new wr stars have yet been found this way .    our second , much more successful , technique relied on using known wr stars to calibrate our selection of targets .
forty known wr stars were selected within the survey area , and patterns were found in their narrow band @xmath13 values which distinguished broad subtypes of wr stars .
( see appendix [ earlywork ] ) .
wc stars generally showed strong ( -0.8 magnitudes or less ) emission excess in the c iv filter , and slightly weaker emission ( between -0.4 and -0.8 magnitudes ) in the he i filter , due to the blue side of the c iv line extending into the range of the he i narrow - band filter .
early wn stars generally showed moderate emission in the br@xmath12 filter and the he ii filter , but slight absorption ( 0.1 magnitudes ) in both he i and c iv .    using these criteria ,
173 candidate targets were selected which appeared at least @xmath20 brighter in their narrow - band filter than did other stars in the field and which also fit the criteria suggested by the known wrs .
though we are reliably detecting stars to magnitude 14 - 14.5 in all of the filters ( by judicious choice of exposure times ) , during this exploratory stage candidates were selected to have emission - band magnitudes brighter than k = 11.5 .
this is because the initial spectrographic follow - up is being done with a 1.5 m telescope ( see below ) .
( thousands of uncrowded @xmath20 candidates as faint as faint as k = 14.5 will be the subjects of future papers ) .
there are strong selection effects for those wr subtypes which were used to determine the selection criteria , and as a result , no wc9 or late wn type ( @xmath22wn6 ) stars have yet been found .
improvements now underway in survey image reduction will permit discovery of the less strongly distinguished subtypes .
the spectrographic follow - up data were taken between 28 april and 6 june 2007 , with the near - infrared ( 0.8 - 2.5 @xmath11 m ) simon spectrograph @xcite of the universit de montral mounted on the 1.5 meter telescope at ctio .
simon has a scale of 0.46 " /pixel on the ctio 1.5 meter telescope .
targets were observed in the k band with a resolving power of r @xmath0 1500 .
each target was observed 5 times , with a nod to move the target along the slit between each observation .
total integration times ranged from 10 to 30 minutes per candidate .
all data were reduced using iraf routines .
images were dark subtracted and flat fielded to remove any instrument signature ; spectra were then extracted using the apall task .
this task also provided sky subtraction by fitting the background on either side of the object along the slit .
the five exposures were scaled and median - combined .
standard stars , which were observed periodically throughout the night , were similarly reduced , and corrections made for br@xmath12 absorption in the telluric standard .
object spectra were then divided by the temporally closest standard - star spectra to remove , as best as possible , the atmospheric absorption features , particularly those at 2.008 @xmath11 m and 2.059 @xmath11 m , and a fainter feature at 2.199 @xmath11 m .
wavelength calibrations were done using the atmospheric oh emission observed off - target by the spectrograph during each object observation .
lines were first identified using the identify task and were compared with the coordinate list ohlines.dat included therein ( @xcite ) .
the wavelength solution was then refitted to each observed target using the reidentify task .
in general , the root - mean - square error of the residuals of wavelength fits was less than 1  , and in most cases between 0.1 and 0.3  .
these dispersion solutions were each applied to the corresponding target object , and then corrected for the intrinsic heliocentric motion using the rvcorrect task to identify the heliocentric velocity and the dopcorrect task to doppler shift the wavelength scale .
from our target list of 173 candidates we have discovered 41 new wolf - rayet stars : 15 of type wn and 26 of type wc .
right ascension and declination , as well as j , h and k@xmath23 magnitudes , were obtained from 2mass , and are listed in table [ tbl - pos ] .
all 2mass objects were then referenced in nomad @xcite to obtain b , v , and r magnitudes , when available , which are also included in table [ tbl - pos ] .
spectra are grouped by assigned spectral type , and are presented in figures [ wn4 ] through [ wc82 ] .      the iraf task splot was used to to fit the continuum and then optimize and deblend gaussian functions to fit the observed emission lines in all spectra .
this gave measurements of line centers , equivalent widths ( ews ) and the full - widths - at - half - maximum ( fwhms ) .
a number of errors contribute to reduce the accuracy of these measurements .
blending of many lines creates the largest emission features , and introduces errors into line center measurements up to 10  .
it also skews the shape of the feature to be a poor fit to a gaussian .
in general , the rms error of the residuals for the fits was about 10% of the continuum value . additionally , the continuum was quite difficult to determine precisely due to the abundance of emission lines throughout the spectrum , especially in wc stars , resulting in errors in the equivalent width measurements .
a number of our wn detections were rather faint , with peak fluxes only twice the continuum level , resulting in lower signal - to - noise ratios for these spectra .
finally , ground - based observatories must peer through the murky atmosphere , so that the strong c iv 2.08 @xmath11 m line in wc stars falls on the edge of the equally strong atmospheric absorption feature at 2.06 @xmath11 m .
division by a standard star removes most of this feature ; however , changes in atmospheric conditions between object and standard - star observations result in residual features on the blue side of the emission line .      the strong emission lines visible in the spectra , while easily distinguished as belonging to either type wc or wn , are the result of overlapping blends of emission lines of various elements . in the optical ,
wr stars are categorized by looking at ratios of equivalent widths of various nitrogen species for wn , carbon and oxygen species for wc , and he ii and he i for both . however , the heavy blending of lines present in the k band makes it more difficult to find either isolated spectral lines or distinguishing ratios of blended lines @xcite .
spectral subtypes were assigned following @xcite .
these are presented in tables [ tbl - eqws - wn ] and [ tbl - eqws - wc ] , along with the ratio of equivalent widths used for categorization .
the ratio @xmath24 was used to categorize the wn stars , and the ratio @xmath25 was used to categorize the wc stars . also following @xcite , wn stars with fwhm(he ii 2.189 )
@xmath26 @xmath27 were classified as broad / strong , and the letter b appended to the subtype designation .
subtypes thus assigned are expected to be accurate to within one subtype .
the 40 known wr stars were similarly assigned subtypes , which agreed within one subtype of their published spectral classifications .
those wc stars with especially broad , heavily blended c iv and c iii lines , as presented in figure [ wce ] and which match our observed spectra of wr19 ( wc4 ) and spectra observed in @xcite of wr146(wc4 ) , wr143(wc5 ) and wr150(wc5 ) , are classified more generally as wce .
blended lines complicate the calculations of accurate radial velocities , as there is no longer any fixed line - center with which to compare theoretical and actual wavelengths .
the purest line in our k band spectra , he ii at 21891  , allows measurement of radial velocity with respect to the motion of the sun in principal ; however , this line is extremely weak or not discernable in most of the wc stars .
line - center measurements are further complicated by the difficulty of accurately fitting gaussian curves to wr emission lines , as described in section [ spectra ] , resulting in errors of 5 - 10   in determination of peak wavelengths .
this corresponds to errors on the order of 50 - 150 km / s . in general
, the direction of motion with respect to the sun follows the clockwise rotation of the galaxy ; however , as noted , the error on these measurements is large .
measured radial velocities ( based on the measured line centers , which may be shifted , depending on details of the line formation mechanism ) for wns and , when possible , wcs , are included in tables [ tbl - eqws - wn ] and [ tbl - eqws - wc ] .
distances to all new wr stars were estimated from the 2mass j , h and ks@xmath23 color excesses , using the method of spectroscopic parallax described in @xcite .
intrinsic @xmath28 and @xmath29 colors specific to wr subtype were taken from @xcite and used to calculate e@xmath30 and e@xmath31 .
extinction ratios taken from @xcite then allow two calculations for @xmath32 :    @xmath33    and @xmath34    the average of these values , @xmath35 , was used to calculate the distance modulus , taking the apparent k@xmath23 magnitude from 2mass and the subtype - specific absolute magnitude , m@xmath36 , from @xcite .
derived k@xmath23 band extinctions and galactocentric distances , r@xmath37 , are listed in table [ tbl - dist ] .
calculations for r@xmath37 assume the iau standard solar galactocentric distance r@xmath38=8.0 kpc . , and the known galactic wr stars are on the same scale in figure 2 . because these calculations are based on the inherent absolute magnitude and @xmath28 and @xmath29 colors for each subtype , the errors are highly dependent on the accuracy of these measurements .
the measured scatter about the adopted color values is approximately 0.02 magnitudes ( table a1 in @xcite ) , which is negligible in a@xmath36 compared to the average scatter of 0.4 magnitudes from the adopted m@xmath36 values .
the redundant calculations of @xmath32 also provide some indication of the reliability of our measurement .
the two values are generally in agreement to within 0.2 magnitudes , especially for the wns , however they can differ by as much as 0.66 magnitudes for some of the wcs , indicating that more accurate subtype specific colors and absolute magnitudes are needed .
these uncertainties give typical errors on the order of 25% in our distance measurements , though they may range as high as 40% in some cases .
the distances to new stars were constrained by the limiting observable magnitude of @xmath3911.5 . using the overall average extinction @xmath40
, we can calculate the typical measurable distances by subtype . as all discovered wn stars were early types
, we can distinguish the faintest observed , strong and weak - lined wn as having average distances of 8.0 and 9.4 kpc , respectively , while observations of the faintest wc stars yield typical distances of 8.9 kpc .    in figure [ galaxyplot ] our new wr stars ( in bold ) have been over - plotted with the previously - known wrs onto the plane of the galaxy , with distances to the known stars taken from the 7th catalogue of galactic wolf - rayet stars @xcite
. the galactic center is labeled , and circles of radius 4 , 8 , and 12 kpc are plotted .
the new wr stars largely follow the distribution pattern established by the known stars , though we can see that we are beginning to push out to larger heliocentric distances .
we also find new stars without optical counterparts within a few kpc of our sun , reinforcing the necessity of wr surveys in the near infrared .
@xcite , along with the more recent reanalysis in @xcite , maintain that wr stars trace the spiral structure of the galaxy .
one arm may be seen along roughly the 8 kpc radius , and an inner arm can perhaps begin to be traced along the inner 4kpc radius
. however , the distance error bars are not trivial , so that firm conclusions about the utility of wr stars as spiral tracers should not yet be drawn .
we have discovered 41 new galactic wr stars , 15 of type wn and 26 of type wc , using a new , near - infrared narrow - band survey of the galactic plane . the reduced extinction from dust and gas in the near infrared makes this the optimal method for future discovery of the thousands of undetected galactic wr stars .
of the 254 total candidates observed spectrographically , 75 proved to be emission line objects .
all of the emission line objects that were not wr stars ( 34 objects ) were planetary nebulae ( pn ) . as the key goal of this survey is the detection of new galactic wr stars
, we amended our selection criteria to eliminate likely pn .
our modified selection criteria yielded 173 wr star candidates which were observed spectrographically : 41 proved to be new wr stars .
with such a 23% detection rate , we have barely scratched the surface of the wealth of new wr stars expected to be discovered within our survey area with the available data .
an initially fairly simple sky - subtraction methodology resulted in relatively scattered color - magnitude diagrams , raising our cut for emission objects to 5@xmath19 .
it also meant that most of our non - detections were erroneously selected objects with featureless spectra . improved sky subtraction ( using entire nights of data , median - filtered in each filter as skyflats ) will allow us to lower this limit to 3@xmath19 and will improve the detection rate of emission - line objects .
we expect this survey to yield thousands of additional discoveries in the coming years .
our survey limits will be pushed fainter by the use of a larger infrared telescope for spectroscopic follow - up .
as we increase the number of known stars , we will also increase the statistical significance of distribution plots , and subtype abundances , allowing us to learn more about our galaxy s structure and composition .
the galactic center is expected to prove an especially rich area for discovery , but it is still largely terra incognita as the crowding of stars there is very high .
the vast majority of galactic wolf - rayet stars remain to be discovered , but we now have a proven technique to continue the search .
benjamin , r.a . ,
churchwell , e. , babler , b.l . ,
bania , t. m. , clemens , d. p. , cohen , m. , dickey , j. m. , indebetouw , r. , jackson , j. m. , kobulnicky , h. a. , lazarian , a. , marston , a. p. , mathis , j. s. , meade , m. r. , seager , s. , stolovy , s. r. , watson , c. , whitney , b. a. , wolff , m. j. , & wolfire , m. g. 2003 , , 115 , 953 buckalew , b.a . ,
kobulnicky , h.a .
, dufour , r.j .
2005 , , 157 , 30 clark , j.s . &
negueruela , i. 2002 , , 396 , l25 crowther , p.a . ,
hadfield , l.j .
, clark , j.s . , negueruela , i. & vacca , w.d .
2006 , , 372 , 1407 conti , p.  s. , & vacca , w.  d.  1990 , , 100 , 431 doyon , r. , nadeau , d. & vallee , ph .
2000 , asp conf .
195 , 548 drissen , l. , moffat , a.f.j . &
shara , m. m. 1993 , , 105 , 1400 drissen , l. , roy , j - r , moffat , a.f.j . & shara , m. m. 1999 , , 117 , 1249 drimmel , r. , & spergel , d. n. 2001 , , 556 , 181 figer , d.f . ,
mclean , i.s . & najarro , f. 1997 , , 486 , 420 hadfield , l.j . , crowther , p.a . ,
schild , h. & schmutz , w. 2005 , , 439 , 265 hadfield , l.j . ,
van dyk , s.d . ,
morris , p.w . ,
smith , j.d . ,
marston , a.p . ,
peterson , d.e .
2007 , , 376 , 248 homeier , n.l . ,
blum , r.d . ,
pasquali , a. , conti , p.s . & damineli , a. 2003 , , 408 , 153 indebetouw , r. , et al .  2005 , , 619 , 931 massey , p. , conti , p.s . ,
moffat , a.f.j . & shara , m.m .
1987 , , 99 , 816 meynet , g. , & maeder , a. 2005 , , 429 , 581 moffat , a.f.j .
, drissen , l.  shara , m.m .
1994 , , 436 , 183 moffat , a.f.j . & shara , m.m .
1987 , , 320 , 266 shara , m.m .
, moffat , a.f.j . ,
smith , l.f . ,
niemela , v.s . , potter , m. & lamontagne , r. 1999 , , 118 , 390 shorlin , s.l . ,
turner , d.g . , &
pedreros , m.h .
2004 , , 116 , 170 skrutskie , m. f. , cutri , r. m. , stiening , r. , weinberg , m. d. , schneider , s. , carpenter , j. m. , beichman , c. , capps , r. , chester , t. , elias , j. , huchra , j. , liebert , j. , lonsdale , c. , monet , d. g. , price , s. , seitzer , p. , jarrett , t. , kirkpatrick , j. d. , gizis , j. e. , howard , e. , evans , t. , fowler , j. , fullmer , l. , hurt , r. , light , r. , kopan , e. l. , marsh , k. a. , mccallon , h. l. , tam , r. , van dyk , s. , & wheelock , s. 2006 , , 131 , 1163 smartt , s.j . , & rolleston , w.r.j .
1997 , , 481 , 47 steed , a.j . &
baker , d.j .
1979 , , 18 , 3386 tutukov , a.v . , & cherepashchuk , a.m. 2003 , astronomy reports , 47 , 386 vanbeveren , d. , van bever , j. & belkus , h. 2007 , , 662 , 107 van der hucht , k.a .
2001 , new astronomy review , 45 , 135 van der hucht , k.a .
2006 , , 458 , 453 zacharias , n. , monet , d.g . , levine , s.e . , urban , s.e . , gaume , r. , & wycoff , g.l .
2005 , vizier catalog i/297    mms , jg and dz acknowledge with gratitude hilary lipsitz , whose ongoing support has been essential to the success of this program .
afjm , rd , ld are grateful to nserc ( canada ) and fqrnt ( quebec ) for financial aid .
we also thank the american museum of natural history for essential funding , and a careful referee for excellent suggestions .
most of the infrared imaging was expertly carried out by claudio aguilera , alberto miranda and alberto pasten .
* appendix a *
in tables 6 and 7 we give a summary of the emission line @xmath13 values measured for the known wrs found in the survey area .
all survey stars were compared with finder charts for the known wr stars to ensure that photometric measurements were made for the correct star .
all stars were found by eye to be present in the survey images , though some had not been identified photometrically due to nebulosity and occasional out - of - focus images .
the @xmath13 magnitude for each narrow - band filter gives the difference in magnitude from the interpolated continuum magnitude .
negative values of @xmath13 indicate emission , positive values indicate absorption for lines in the narrow - band emission - line filter , assuming there are no other problems .
several fields were imaged multiple times , and on more than one night .
this allowed us to estimate the accuracy of the emission line magnitudes in tables 6 and 7 as @xmath410.03 magnitudes .
* appendix b *
following shara et al . ( 1999 ) , we assume that the star distribution expressed in cylindrical galactic coordinates ( i.e. @xmath42 ) follows that of the interstellar dust and adopt an axisymmetric exponential disk formulation with an outwards flair taken from drimmel & spergel ( 2001 ) :    @xmath43    with , + @xmath44    a linear flair is added by imposing the following functional form to @xmath45 :    @xmath46 h_o   &   r(r , l ) \leq r_f \end{array }   \right.\ ] ]    in these equations , @xmath47 is the local stellar density , @xmath48 pc and @xmath49 pc are respectively the solar galactocentric distance and distance from the galactic plane .
the scale length @xmath50 and the constants @xmath51 , @xmath52 and @xmath53 which define the scale height @xmath45 are taken from the fit to the fir emission of the interstellar dust of drimmel & spergel ( 2001 ) and are respectively @xmath54 pc , @xmath55 pc , @xmath56 and @xmath57 pc .    although the assumption that the star distribution follows the dust is justified for o stars , it must be modified for wr stars to account for the observed metallicity dependence of their galactic distribution . using the results of maeder & meynet ( 1994 ) for the metallicity dependence of the wr / o number ratio and the radial metallicity distribution of the galactic disk from smartt & rolleston ( 1997 ) , shara et al .
( 1999 ) found that the wr star density can be expressed as :    @xmath58    where , + @xmath59    knowing the local surface density of wr stars from van der hucht ( 2001 ) :    @xmath60    and applying this constraint to @xmath61 we obtain :    @xmath62    finally , the conversion of the wr star density from a spatial ( @xmath61 ) to a magnitude ( @xmath63 ) dependence is done in two steps .
first , we use the conservation of stars in both systems to write :    @xmath64    then , the radial distance @xmath65 is converted to the magnitude @xmath66 using the inverse square - law of light attenuation accounting for interstellar extinction :    @xmath67    where the extinction @xmath68 allowing for a spherical hole in the interstellar dust at the center of the galaxy is ( see drimmel & spergel 20001 ) :    @xmath69 a_k(0.5r_o , l , z)e^{- \left ( \frac{r(r , l ) - 0.5r_o}{2500}\right)^2 }   & r(r , l ) < 0.5r_o \end{array } \right.\ ] ]    with @xmath70 mag / pc ( mathis 1990 ) and the intrinsic magnitude of wr stars in the @xmath71-band @xmath72 mag ( van der hucht 2001 ) .
these equations have been solved numerically to determine the star density @xmath73 to an accuracy better than @xmath74 using a monte - carlo method with a sobol quasi - random number generator in @xmath75 . by integrating over the whole galaxy the model predicts @xmath76 wr stars which is very close to the @xmath77 predicted by van de hucht ( 2001 ) .
further validation can be carried out by applying this model to the @xmath78-band where most wr stars have been detected so far .
it is then possible to compare the observed to the predicted number of wr stars .
to do so , we adopt @xmath79 mag ( van der hucht 2001 ) and @xmath80 mag / pc . according to the model ,
the number of wr stars observable up to a magnitude of v = 15 , 12 and 10 are respectively 153 , 79 and 42 .
these numbers are very close to the respective observed numbers reported by van der hucht ( 2001 ) of 159 , 80 and 36 .
it is again apparent how important it is to continue the search for new wr stars in the infrared .
figure [ con95 ] shows the contour plot representing the necessary target k magnitude to detect @xmath81 of all wr stars along a line of sight .
figure [ numwr ] shows the number of wr stars expected per cpapir field and figure [ numwr2 ] presents the same information for @xmath82 .
figure 6 shows the cumulative number of expected galactic wr stars as a function of k magnitude . at k= 11 , 12 , 13 and 14 ,
the expected numbers of galactic wr stars are 1200 , 2500 , 4200 and 5400 .
* appendix c *
we present in figures 7a through 7 g the finder charts for the 41 new wolf - rayet stars discovered in our survey .
l r r r r c c c r r r 668_4 & 10 16 26.22 & -57 28 05.7 & 283.26 & -0.64 & - & - & 17.46 & 11.73 & 10.38 & 9.53 + 740_21 & 11 16 03.53 & -61 26 58.3 & 291.79 & -0.66 & 16.76 & 16.82 & 15.97 & 13.16 & 11.57 & 10.47 + 740_16 & 11 19 42.96 & -61 27 12.4 & 292.20 & -0.51 & 13.67 & 12.35 & 13.03 & 11.50 & 10.99 & 10.33 + 768_6 & 11 46 06.66 & -62 47 12.7 & 295.54 & -0.86 & 15.56 & 17.17 & 14.11 & 12.15 & 11.36 & 10.78 + 772_17 & 11 50 04.24 & -62 52 15.4 & 296.00 & -0.83 & 18.23 & - & 16.47 & 12.57 & 11.58 & 10.94 + 776_3 & 11 55 52.11 & -62 45 02.2 & 296.62 & -0.56 & 18.38 & - & 14.94 & 11.65 & 10.76 & 10.14 + 791_12c & 12 13 28.29 & -62 41 42.9 & 298.59 & -0.14 & - & - & - & 14.05 & 12.67 & 11.53 + 808_14 & 12 28 41.91 & -63 25 46.1 & 300.39 & -0.67 & 20.28 & - & 17.82 & 12.78 & 11.83 & 10.94 + 808_23 & 12 28 50.99 & -63 17 00.2 & 300.40 & -0.52 & - & - & 18.61 & 13.56 & 12.55 & 11.84 + 807_16 & 12 30 03.86 & -62 50 17.1 & 300.50 & -0.07 & - & - & 19.20 & 13.16 & 12.05 & 11.09 + 816_10 & 12 38 18.78 & -63 24 19.7 & 301.46 & -0.57 & - & - & 18.43 & 12.93 & 11.76 & 10.96 + 832_25 & 12 55 44.26 & -63 35 50.0 & 303.41 & -0.73 & - & - & - & 12.48 & 11.12 & 10.22 + 856_13c & 13 03 11.08 & -63 42 16.2 & 304.23 & -0.86 & - & - & 18.72 & 13.13 & 11.92 & 11.05 + 839_12 & 13 04 50.08 & -63 04 40.2 & 304.45 & -0.25 & - & - & - & 15.69 & 14.00 & 12.54 + 845_34 & 13 12 21.30 & -62 40 12.5 & 305.33 & 0.10 & - & - & 16.24 & 10.75 & 9.57 & 8.77 + 845_35 & 13 12 27.67 & -62 44 22.0 & 305.34 & 0.03 & - & - & 17.95 & 13.16 & 11.82 & 10.71 + 847_8 & 13 12 45.35 & -63 05 52.0 & 305.34 & -0.33 & - & - & - & 13.06 & 11.34 & 10.26 + 853_9 & 13 22 16.08 & -62 30 57.4 & 306.48 & 0.14 & - & - & - & 14.77 & 12.69 & 11.55 + 858_26 & 13 28 15.87 & -62 06 23.5 & 307.22 & 0.46 & - & - & 18.70 & 12.53 & 11.40 & 10.68 + 883_18 & 13 52 02.36 & -62 26 46.0 & 309.88 & -0.39 & - & - & - & 14.59 & 12.41 & 10.99 + 885_11 & 13 54 13.45 & -61 50 01.8 & 310.27 & 0.14 & - & - & - & 13.60 & 12.07 & 10.86 + 897_5 & 14 10 10.01 & -61 15 25.5 & 312.25 & 0.18 & - & - & - & 15.30 & 12.43 & 10.60 + 903_15c & 14 12 36.54 & -61 45 32.7 & 312.38 & -0.38 & - & - & - & 13.58 & 11.33 & 9.71 + 907_18 & 14 16 27.37 & -61 17 56.2 & 312.96 & -0.09 & - & - & - & 13.60 & 11.39 & 10.01 + 956_25 & 15 01 30.11 & -59 16 12.0 & 318.88 & -0.49 & - & - & - & 13.84 & 11.83 & 10.64 + 979_11 & 15 20 35.91 & -57 27 11.9 & 321.95 & -0.20 & - & - & - & 13.88 & 12.15 & 11.06 + 1011_24 & 15 43 04.68 & -55 11 12.3 & 325.81 & -0.13 & 19.93 & 17.31 & 16.06 & 10.96 & 9.94 & 9.06 + 1053_27 & 16 11 43.70 & -51 10 16.6 & 331.67 & 0.17 & - & - & 18.53 & 10.26 & 8.97 & 8.10 + 1059_34 & 16 14 37.23 & -51 26 26.3 & 331.81 & -0.34 & - & - & - & 15.05 & 12.79 & 11.54 + 1081_21 & 16 24 58.86 & -48 56 52.4 & 334.75 & 0.27 & - & - & - & 13.28 & 11.76 & 10.73 + 1093_34 & 16 31 29.23 & -47 56 16.4 & 336.22 & 0.19 & - & - & - & 15.56 & 12.86 & 11.32 + 1093_33 & 16 31 49.06 & -47 56 04.4 & 336.26 & 0.15 & - & - & - & 15.11 & 12.86 & 11.47 + 1093_53 & 16 32 12.98 & -47 50 35.8 & 336.37 & 0.17 & - & - & - & 15.00 & 12.70 & 11.34 + 1096_22 & 16 35 23.31 & -48 09 18.0 & 336.51 & -0.43 & - & - & - & 14.90 & 12.77 & 11.43 + 1222_15 & 17 22 40.74 & -35 04 52.9 & 352.20 & 0.74 & - & - & - & 15.11 & 12.05 & 10.33 + 1385_24 & 18 13 42.47 & -17 28 12.2 & 13.15 & 0.13 & - & - & 20.02 & 11.21 & 9.70 & 8.57 + 1425_47 & 18 23 03.42 & -13 10 00.4 & 18.01 & 0.18 & 15.23 & 14.49 & 14.41 & 10.34 & 9.28 & 8.27 + 1462_54 & 18 29 33.84 & -08 39 02.1 & 22.75 & 0.87 & - & - & - & 14.07 & 12.75 & 11.93 + 1509_29 & 18 41 48.45 & -04 00 12.9 & 28.27 & 0.31 & - & - & - & 15.62 & 13.32 & 11.99 + 1613_21 & 19 06 36.53 & + 07 29 52.4 & 41.33 & 0.07 & - & - & - & 14.24 & 12.65 & 11.61 + 1671_5 & 19 20 40.38 & + 13 50 35.2 & 48.55 & -0.05 & - & - & - & 13.57 & 11.80 & 10.76 +    l r r r r r r r r r r l 668_4 & -8 & 141 & -11 & 121 & -31 & 147 & -91 & 153 & 2.9 & 55 & wn5b + 740_21 & -14 & 208 & -20 & 200 & -39 & 165 & -203 & 242 & 5.2 & -82 & wn4b + 768_6 & -6 & 38 & -19 & 126 & -42 & 120 & -101 & 106 & 2.4 & 41 & wn5 + 772_17 & -3 & 38 & -12 & 113 & -53 & 126 & -95 & 121 & 1.8 & - & wn6 + 776_3 & - & - & -19 & 109 & -60 & 93 & -72 & 94 & 1.2 & - & wn6 + 808_23 & - & - & - & - & -46 & 101 & -141 & 126 & 3.1 & 68 & wn5 + 816_10 & - & - & -29 & 227 & -48 & 140 & -137 & 136 & 2.9 & -14 & wn5b + 847_8 & -8 & 70 & -47 & 142 & -65 & 119 & -114 & 120 & 1.8 & -96 & wn6 + 853_9 & - & - & -29 & 95 & -36 & 66 & -92 & 97 & 2.6 & 27 & wn6 + 858_26 & - & - & -26 & 96 & -38 & 74 & -56 & 78 & 1.5 & 55 & wn6 + 907_18 & -16 & 171 & -19 & 140 & -70 & 175 & -150 & 153 & 2.1 & -68 & wn5b + 956_25 & -11 & 202 & -35 & 222 & -31 & 163 & -180 & 366 & 5.8 & -315 & wn4b + 979_11 & - & - & -57 & 260 & -62 & 207 & -212 & 220 & 3.4 & -14 & wn4b + 1093_53 & - & - & -52 & 219 & -60 & 185 & -145 & 182 & 2.4 & -68 & wn5b + 1462_54 & - & - & -59 & 202 & -72 & 164 & -166 & 126 & 2.3 & 68 & wn5
+    l r r r r r r r r r r l 740_16 & -987 & 255 & -226 & 330 & - & - & -50 & 241 & 4.4 & - & wce + 791_12c & -1590 & 282 & -498 & 443 & - & - & -75 & 228 & 3.2 & - & wce + 808_14 & -948 & 306 & -392 & 477 & - & - & -140 & 623 & 2.4 & - & wce + 807_16 & -985 & 178 & -206 & 147 & -31 & 103 & -74 & 116 & 4.8 & -55 & wc7 + 832_25 & -454 & 236 & -73 & 236 & -37 & 351 & -15 & 121 & 6.2 & - & wc5 - 6 + 839_12 & -1658 & 230 & -244 & 195 & - & - & -71 & 157 & 6.8 & - & wc5 - 6 + 845_34 & -374 & 271 & -131 & 284 & - & - & -35 & 384 & 2.9 & - & wc8 + 845_35 & -1789 & 249 & -412 & 258 & -42 & 323 & -70 & 195 & 4.3 & 0 & wc7 + 856_13c & -712 & 191 & -138 & 180 & -14 & 86 & -61 & 157 & 5.2 & -164 & wc5 - 6 + 883_18 & -1093 & 301 & -373 & 453 & -10 & 103 & -77 & 320 & 2.9 & - & wce + 885_11 & -1489 & 210 & -279 & 180 & - & - & -60 & 174 & 5.3 & - & wc5 - 6 + 897_5 & -525 & 205 & -188 & 165 & -56 & 354 & -35 & 126 & 2.8 & - & wc8 + 903_15c & -165 & 186 & -66 & 123 & -13 & 125 & -14 & 125 & 2.5 & - & wc8 + 1011_24 & -884 & 244 & -275 & 226 & - & - & -83 & 360 & 3.2 & - & wc8 + 1053_27 & -303 & 222 & -81 & 213 & -59 & 223 & -84 & 202 & 3.7 & - & wc8 + 1059_34 & -729 & 210 & -242 & 181 & -28 & 366 & -44 & 199 & 3.0 & - & wc8 + 1081_21 & -419 & 219 & -214 & 174 & -75 & 185 & -60 & 143 & 2.0 & -137 & wc8 + 1093_33 & -369 & 199 & -155 & 161 & -69 & 161 & -69 & 132 & 2.4 & -219 & wc8 + 1093_34 & -475 & 226 & -178 & 222 & -62 & 212 & -84 & 187 & 2.7 & -260 & wc8 + 1096_22 & -571 & 202 & -230 & 156 & -71 & 169 & -76 & 141 & 2.5 & -82 & wc8 + 1222_15 & -551 & 197 & -171 & 154 & - & - & -45 & 230 & 3.2 & - & wc8 + 1385_24 & -537 & 220 & -245 & 196 & -63 & 361 & -33 & 138 & 2.2 & - & wc8 + 1425_47 & -1593 & 238 & -300 & 296 & -48 & 250 & -75 & 179 & 5.3 & - & wc5 - 6 + 1509_29 & -1578 & 240 & -355 & 224 & - & - & -110 & 273 & 4.4 & - & wc7 + 1613_21 & -659 & 365 & -245 & 518 & - & - & - & - & 2.7 & - & wce + 1671_5 & -335 & 215 & -94 & 210 & -49 & 337 & -20 & 151 & 3.6 & - & wc8 +    l l c c r r r r r r 668_4 & wn5b & 1.23 & 1.06 & 1.14 & 9.53 & -4.77 & 13.16 & 4.28 & 8.60 + 740_21 & wn4b & 1.55 & 1.51 & 1.53 & 10.47 & -4.77 & 13.71 & 5.51 & 8.24 + 740_16 & wce & 0.37 & 0.15 & 0.26 & 10.33 & -4.59 & 14.66 & 8.56 & 9.52 + 768_6 & wn5 & 0.80 & 0.76 & 0.78 & 10.78 & -4.41 & 14.41 & 7.62 & 8.63 + 772_17 & wn6 & 0.97 & 0.87 & 0.92 & 10.94 & -4.41 & 14.43 & 7.68 & 8.60 + 776_3 & wn6 & 0.89 & 0.84 & 0.86 & 10.14 & -4.41 & 13.69 & 5.46 & 7.78 + 791_12c & wce & 1.27 & 1.02 & 1.15 & 11.53 & -4.59 & 14.97 & 9.88 & 9.46 + 808_14 & wce & 0.82 & 0.56 & 0.69 & 10.94 & -4.59 & 14.84 & 9.29 & 8.87 + 808_23 & wn5 & 1.03 & 1.00 & 1.02 & 11.84 & -4.41 & 15.23 & 11.14 & 10.02 + 807_16 & wc7 & 0.97 & 0.69 & 0.83 & 11.09 & -4.59 & 14.85 & 9.33 & 8.88 + 816_10 & wn5b & 1.07 & 0.96 & 1.02 & 10.96 & -4.77 & 14.71 & 8.76 & 8.44 + 832_25 & wc5 - 6 & 1.10 & 0.58 & 0.84 & 10.22 & -4.59 & 13.97 & 6.22 & 7.26 + 856_13c & wc5 - 6 & 0.98 & 0.53 & 0.75 & 11.05 & -4.59 & 14.89 & 9.49 & 8.46 + 839_12 & wc5 - 6 & 1.70 & 1.60 & 1.65 & 12.54 & -4.59 & 15.48 & 12.48 & 10.39 + 845_34 & wc8 & 1.04 & 0.76 & 0.90 & 8.77 & -4.65 & 12.52 & 3.19 & 7.15 + 845_35 & wc7 & 1.23 & 0.96 & 1.10 & 10.71 & -4.59 & 14.20 & 6.93 & 7.22 + 847_8 & wn6 & 1.76 & 1.67 & 1.71 & 10.26 & -4.41 & 12.96 & 3.90 & 7.01 + 853_9 & wn6 & 2.04 & 1.78 & 1.91 & 11.55 & -4.41 & 14.05 & 6.46 & 6.98 + 858_26 & wn6 & 1.12 & 1.02 & 1.07 & 10.68 & -4.41 & 14.02 & 6.37 & 6.88 + 883_18 & wce & 2.00 & 1.53 & 1.76 & 10.99 & -4.59 & 13.82 & 5.80 & 6.53 + 885_11 & wc5 - 6 & 1.42 & 1.15 & 1.28 & 10.86 & -4.59 & 14.17 & 6.81 & 6.62 + 897_5 & wc8 & 2.86 & 2.64 & 2.75 & 10.60 & -4.65 & 12.50 & 3.16 & 6.79 + 903_15c & wc8 & 2.30 & 2.26 & 2.28 & 9.71 & -4.65 & 12.08 & 2.61 & 7.01 + 907_18 & wn5b & 2.16 & 2.02 & 2.09 & 10.01 & -4.77 & 12.69 & 3.45 & 6.65 + 956_25 & wn4b & 1.90 & 1.67 & 1.79 & 10.64 & -4.77 & 13.62 & 5.31 & 5.70 + 979_11 & wn4b & 1.64 & 1.49 & 1.57 & 11.06 & -4.77 & 14.26 & 7.12 & 5.26 + 1011_24 & wc8 & 0.98 & 0.91 & 0.95 & 9.06 & -4.65 & 12.76 & 3.57 & 5.90 + 1053_27 & wc8 & 1.16 & 0.89 & 1.03 & 8.10 & -4.65 & 11.72 & 2.21 & 6.64 + 1059_34 & wc8 & 2.06 & 1.58 & 1.82 & 11.54 & -4.65 & 14.37 & 7.47 & 4.01 + 1081_21 & wc8 & 1.42 & 1.18 & 1.30 & 10.73 & -4.65 & 14.08 & 6.54 & 3.80 + 1093_34 & wc8 & 2.55 & 2.11 & 2.33 & 11.32 & -4.65 & 13.64 & 5.34 & 4.21 + 1093_33 & wc8 & 2.15 & 1.84 & 1.99 & 11.47 & -4.65 & 14.13 & 6.69 & 3.59 + 1093_53 & wn5b & 2.20 & 1.98 & 2.09 & 11.34 & -4.77 & 14.02 & 6.36 & 3.70 + 1096_22 & wc8 & 2.04 & 1.75 & 1.89 & 11.43 & -4.65 & 14.19 & 6.88 & 3.51 + 1222_15 & wc8 & 2.91 & 2.44 & 2.68 & 10.33 & -4.65 & 12.30 & 2.89 & 5.65 + 1385_24 & wc8 & 1.48 & 1.37 & 1.42 & 8.57 & -4.65 & 11.80 & 2.29 & 6.29 + 1425_47 & wc5 - 6 & 0.97 & 0.78 & 0.88 & 8.27 & -4.59 & 11.98 & 2.49 & 6.18 + 1462_54 & wn5 & 1.31 & 1.20 & 1.26 & 11.93 & -4.41 & 15.08 & 10.39 & 4.16 + 1509_29 & wc7 & 2.02 & 1.36 & 1.69 & 11.99 & -4.59 & 14.89 & 9.50 & 4.50 + 1613_21 & wce & 1.35 & 0.84 & 1.09 & 11.61 & -4.59 & 15.11 & 10.51 & 6.97 + 1671_5 & wc8 & 1.59 & 1.20 & 1.40 & 10.76 & -4.65 & 14.01 & 6.34 & 6.41
+     l l r r r r 48c & wn3h+wc4 & 0.01 & 0.02 & -0.37 & -0.20 + 18 & wn4 & -0.07 & 0.06 & -0.13 & -0.32 + 35b & wn4 & 0.10 & 0.12 & -0.18 & -0.48 + 44a & wn4 & 0.19 & -0.22 & -0.42 & -0.19 + 45b & wn4 & 0.11 & - & -0.19 & 0.11 + 62a & wn4 & -0.03 & 0.01 & -0.08 & -0.13 + 31 & wn4+o8v & -0.01 & 0.03 & -0.12 & -0.21 + 51 & wn4+ob ? & -0.08 & -0.01 & -0.15 & -0.38 + 38a & wn5 & 0.04 & 0.02 & -0.19 & -0.25 + 42c & wn5 & 0.09 & 0.12 & -0.16 & -0.35 + 42d & wn5 & 0.03 & 0.06 & -0.09 & -0.29 + 45a & wn5 & 0.08 & 0.05 & -0.17 & -0.41 + 36 & wn5 - 6+ob ? & 0.00 & 0.08 & -0.25 & -0.44 + 111c & wn6 & 0.08 & 0.12 & -0.21 & -0.38 + 35a & wn6h & 0.00 & 0.01 & -0.15 & -0.13 + 21a & wn6+o / a & 0.09 & - & -0.05 & 0.02 + 28 & wn6(h)+ob ?
& 0.06 & 0.06 & -0.21 & -0.22 + 47 & wn6+o5v & -0.10 & 0.01 & -0.18 & -0.19 + 55 & wn7 & -0.08 & 0.03 & -0.23 & -0.17 + 74 & wn7 & -0.06 & 0.00 & -0.17 & -0.16 + 84 & wn7 & -0.07 & 0.03 & -0.16 & -0.22 + 111d & wn7 ? &
0.03 & 0.09 & -0.04 & -0.09 + 19a & wn7:(h ) & -0.06 & 0.00 & -0.20 & -0.10 + 26 & wn7/wce & 0.03 & -0.59 & -0.16 & -0.50 + 107 & wn8 & -0.23 & -0.07 & -0.22 & -0.02 + 47b & wn9h & 0.08 & 0.03 & -0.26 & -0.06 + 108 & wn9h+ob & -0.11 & -0.02 & -0.01 & 0.12 +     l l r r r r 38 & wc4 & -1.06 & -1.56 & - & -0.20 + 47c & wc5 & -0.16 & -0.42 & -0.05 & -0.08 + 32 & wc5+ob ? & -0.41 & -1.77 & -0.13 & -0.37 + 41 & wc5+ob ? & -1.06 & -2.52 & -0.40 & -79.86 + 23 & wc6 & -0.46 & - & -0.12 & -0.28 + 27 & wc6+a & -0.71 & -1.36 & -0.13 & -0.31 + 31c & wc6+ob & -1.04 & -1.51 & -0.19 & -0.28 + 38b & wc7+ob & -0.77 & -1.14 & -0.08 & -0.11 + 39 & wc7+ob ? & -0.31 & -0.76 & -0.13 & -0.17 + 42 & wc7+o7v & -0.36 & - & -0.11 & -0.12 + 50 & wc7+ob & -0.76 & -1.42 & -0.02 & -0.15 + 111b & wc9d & 0.12 & 0.18 & -0.11 & -0.19 + 48b & wc9d & -0.03 & -0.06 & - & -0.24 +