pro m67_props,ra,dec,v,bmv,teff,logg,logz,ext ; This routine returns the observed stellar properties of all stars ; in M67 that have spectroscopic teff values from Latham and V, B-V ; from Montgomery et al. ; log(g) is inferred from V, teff, and assumed distance modulus = 9.69 ; and assumed reddening = .04, A_V = .12 ;constants vsun=4.82 ; solar V abs mag tsun=5777. ; solar teff dm=9.69 ; M67 dist modulus, fro Sarajadini et al. extin=.12 ; A_V bc=[-2.48,-1.62,-1.25,-1.02,-.61,-.5,-.42,-.31,-.28,-.21,-.2,-.18,$ -.16,-.14,-.11,-.09,-.15,-.2,-.3,-.8,-1.46,-2.35] tbc=[3380.,3540.,3690.,4050.,4390.,4660.,4830.,5150.,5310.,5560.,5790.,$ 5940.,6250.,6650.,7000.,7300.,8180.,9000.,9790.,11400.,15200.,20900.] bv=[1.63,1.6,1.56,1.5,1.16,1.0,.91,.81,.74,.68,.63,.58,.52,.44,.35,.3,.15,$ .05,-.02,-.11,-.17,-.24] mss=[.30,.40,.46,.56,.67,.70,.75,.79,.85,.92,.99,1.05,1.15,1.4,1.5,1.6,2.0,$ 2.55,2.9,3.8,5.9,8.5] ; read Latham's table filin='~/costi_linux/tomthop/phot/m67.sumt' rd_specdat,filin,ral,decl,teffl,logzl,loggl,chi2l ral=ral*15. nstarl=n_elements(ral) ; read photometry table from Montgomery et al. rd_m67_montgomery,filmont,ram,decm,vm,bmvm nstarm=n_elements(ram) ; convert, precess coords to J2000.0 ram=ram*15. precess,ram,decm,1950.0,2000.0 ; make teff from B-V teffm=interpol(tbc,bv,bmvm) ; make radii for Montgomery stars from V magnitudes and B-V colors bcm=interpol(bc,bv,bmvm) radm=(tsun/teffm)^2*10.^(-(vm-vsun+bcm-dm-extin)/5.) ; find positional matches between Latham's and Montgomery's stars; ; associate V, B-V, and radius with matched stars rag=dblarr(nstarl) decg=dblarr(nstarl) tefflg=fltarr(nstarl) teffmg=fltarr(nstarl) loggg=fltarr(nstarl) logzg=fltarr(nstarl)-0.1 massg=fltarr(nstarl) vg=fltarr(nstarl) bmvg=fltarr(nstarl) radg=fltarr(nstarl) for i=0,nstarl-1 do begin dra=(ram-ral(i))*3600. ddec=(decm-decl(i))*3600. dist=sqrt(dra^2+ddec^2) s=where(dist le 10.,ns) ; capture radius is 3 arcsec if(ns eq 1) then begin rag(i)=ral(i) decg(i)=decl(i) vg(i)=vm(s) bmvg(i)=bmvm(s) tefflg(i)=teffl(i) teffmg(i)=teffm(s) radg(i)=radm(s) endif endfor s1=where(rag ne 0.,ns1) if(ns1 gt 0) then begin rag=rag(s1) decg=decg(s1) tefflg=tefflg(s1) teffmg=teffmg(s1) loggg=loggg(s1) logzg=logzg(s1) massg=massg(s1) vg=vg(s1) bmvg=bmvg(s1) radg=radg(s1) endif ; make masses using mass-teff relation for average of Latham & Montgomery ; teff, except set mass=1.25 for stars bright enough to be above turnoff tavg=(tefflg+teffmg)/2. massg=interpol(mss,tbc,tavg) s=where(vg le 13.,ns) if(ns gt 0) then begin massg(s)=1.25 endif loggg=alog10(massg)-2.*alog10(radg)+4.44 ra=rag dec=decg v=vg bmv=bmvg teff=tavg logg=loggg logz=logzg ext=fltarr(ns1)+extin end