; spcalflux.pro - integrated fluxes (in Jy km/s), linewidths, and Vhel ; for a set of calibrators (spectral line and flux) ; Returns: ; answer, which is a 3-element array: ; answer[0] = line flux (Jy km/s) ; answer[1] = line-width at 50% of peak level (km/s) ; answer[2] = mean velocity (km/s)] ; These are all set to -1 if the source is not found. ; ; Example: ; answer = spcalflux('UGC0317') ; ; Note: Each source is repeated four times, in case they are observed ; repeatidly on the same night. function spcalflux, source rec = {specline, name:"blank", Sdv:0.0, W20:0.0, W50:0.0, Vhel:0.0 } spcal = replicate({specline}, 76) i=0 ;; spectral line calibrators from O'Neil, 2004, AJ, 128, 2080. spcal[i].name='UGC0317' & spcal[i].Vhel=5361.0 & spcal[i].Sdv=1.1 & spcal[i].W20= 37.0 & spcal[i].W50= 33.0 & i+=1 spcal[i].name='UGC0317a' & spcal[i].Vhel=5361.0 & spcal[i].Sdv=1.1 & spcal[i].W20= 37.0 & spcal[i].W50= 33.0 & i+=1 spcal[i].name='UGC0317b' & spcal[i].Vhel=5361.0 & spcal[i].Sdv=1.1 & spcal[i].W20= 37.0 & spcal[i].W50= 33.0 & i+=1 spcal[i].name='UGC0317c' & spcal[i].Vhel=5361.0 & spcal[i].Sdv=1.1 & spcal[i].W20= 37.0 & spcal[i].W50= 33.0 & i+=1 spcal[i].name='UGC2432' & spcal[i].Vhel= 758.0 & spcal[i].Sdv=5.9 & spcal[i].W20=114.0 & spcal[i].W50=105.0 & i+=1 spcal[i].name='UGC2432a' & spcal[i].Vhel= 758.0 & spcal[i].Sdv=5.9 & spcal[i].W20=114.0 & spcal[i].W50=105.0 & i+=1 spcal[i].name='UGC2432b' & spcal[i].Vhel= 758.0 & spcal[i].Sdv=5.9 & spcal[i].W20=114.0 & spcal[i].W50=105.0 & i+=1 spcal[i].name='UGC2432c' & spcal[i].Vhel= 758.0 & spcal[i].Sdv=5.9 & spcal[i].W20=114.0 & spcal[i].W50=105.0 & i+=1 spcal[i].name='UGC3755' & spcal[i].Vhel= 318.0 & spcal[i].Sdv=8.2 & spcal[i].W20= 67.0 & spcal[i].W50= 77.0 & i+=1 spcal[i].name='UGC3755a' & spcal[i].Vhel= 318.0 & spcal[i].Sdv=8.2 & spcal[i].W20= 67.0 & spcal[i].W50= 77.0 & i+=1 spcal[i].name='UGC3755b' & spcal[i].Vhel= 318.0 & spcal[i].Sdv=8.2 & spcal[i].W20= 67.0 & spcal[i].W50= 77.0 & i+=1 spcal[i].name='UGC3755c' & spcal[i].Vhel= 318.0 & spcal[i].Sdv=8.2 & spcal[i].W20= 67.0 & spcal[i].W50= 77.0 & i+=1 ;; CHANGED THE FLUX. O'NEIL PUBLISHES Sdv=7.9; I MEASURE 9.75+/-0.25 WITH GBT spcal[i].name='UGC6421' & spcal[i].Vhel=1502.0 & spcal[i].Sdv=9.75 & spcal[i].W20=175.0 & spcal[i].W50=169.0 & i+=1 spcal[i].name='UGC6421a' & spcal[i].Vhel=1502.0 & spcal[i].Sdv=9.75 & spcal[i].W20=175.0 & spcal[i].W50=169.0 & i+=1 spcal[i].name='UGC6421b' & spcal[i].Vhel=1502.0 & spcal[i].Sdv=9.75 & spcal[i].W20=175.0 & spcal[i].W50=169.0 & i+=1 spcal[i].name='UGC6421c' & spcal[i].Vhel=1502.0 & spcal[i].Sdv=9.75 & spcal[i].W20=175.0 & spcal[i].W50=169.0 & i+=1 ;; CHANGED THE FLUX. O'NEIL PUBLISHES Sdv=7.9; I MEASURE 7.4 WITH GBT spcal[i].name='UGC9901' & spcal[i].Vhel=3160.0 & spcal[i].Sdv=7.4 & spcal[i].W20=246.0 & spcal[i].W50=237.0 & i+=1 spcal[i].name='UGC9901a' & spcal[i].Vhel=3160.0 & spcal[i].Sdv=7.4 & spcal[i].W20=246.0 & spcal[i].W50=237.0 & i+=1 spcal[i].name='UGC9901b' & spcal[i].Vhel=3160.0 & spcal[i].Sdv=7.4 & spcal[i].W20=246.0 & spcal[i].W50=237.0 & i+=1 spcal[i].name='UGC9901c' & spcal[i].Vhel=3160.0 & spcal[i].Sdv=7.4 & spcal[i].W20=246.0 & spcal[i].W50=237.0 & i+=1 spcal[i].name='UGC11461' & spcal[i].Vhel=3122.0 & spcal[i].Sdv=6.6 & spcal[i].W20=322.0 & spcal[i].W50=298.0 & i+=1 spcal[i].name='UGC11461a' & spcal[i].Vhel=3122.0 & spcal[i].Sdv=6.6 & spcal[i].W20=322.0 & spcal[i].W50=298.0 & i+=1 spcal[i].name='UGC11461b' & spcal[i].Vhel=3122.0 & spcal[i].Sdv=6.6 & spcal[i].W20=322.0 & spcal[i].W50=298.0 & i+=1 spcal[i].name='UGC11461c' & spcal[i].Vhel=3122.0 & spcal[i].Sdv=6.6 & spcal[i].W20=322.0 & spcal[i].W50=298.0 & i+=1 ;; spectral line calibrators from Hogg et al. 2007, AJ, 134, 1046 spcal[i].name='UGC2947' & spcal[i].Vhel=863.0 & spcal[i].Sdv=39.3 & spcal[i].W20=178.0 & spcal[i].W50=178.0 & i+=1 spcal[i].name='UGC2947a' & spcal[i].Vhel=863.0 & spcal[i].Sdv=39.3 & spcal[i].W20=178.0 & spcal[i].W50=178.0 & i+=1 spcal[i].name='UGC2947b' & spcal[i].Vhel=863.0 & spcal[i].Sdv=39.3 & spcal[i].W20=178.0 & spcal[i].W50=178.0 & i+=1 spcal[i].name='UGC2947c' & spcal[i].Vhel=863.0 & spcal[i].Sdv=39.3 & spcal[i].W20=178.0 & spcal[i].W50=178.0 & i+=1 spcal[i].name='UGC4165' & spcal[i].Vhel=514.0 & spcal[i].Sdv=32.7 & spcal[i].W20=102.0 & spcal[i].W50=102.0 & i+=1 spcal[i].name='UGC4165a' & spcal[i].Vhel=514.0 & spcal[i].Sdv=32.7 & spcal[i].W20=102.0 & spcal[i].W50=102.0 & i+=1 spcal[i].name='UGC4165b' & spcal[i].Vhel=514.0 & spcal[i].Sdv=32.7 & spcal[i].W20=102.0 & spcal[i].W50=102.0 & i+=1 spcal[i].name='UGC4165c' & spcal[i].Vhel=514.0 & spcal[i].Sdv=32.7 & spcal[i].W20=102.0 & spcal[i].W50=102.0 & i+=1 spcal[i].name='UGC4543' & spcal[i].Vhel=1960.0 & spcal[i].Sdv=29.1 & spcal[i].W20=112.0 & spcal[i].W50=112.0 & i+=1 spcal[i].name='UGC4543a' & spcal[i].Vhel=1960.0 & spcal[i].Sdv=29.1 & spcal[i].W20=112.0 & spcal[i].W50=112.0 & i+=1 spcal[i].name='UGC4543b' & spcal[i].Vhel=1960.0 & spcal[i].Sdv=29.1 & spcal[i].W20=112.0 & spcal[i].W50=112.0 & i+=1 spcal[i].name='UGC4543c' & spcal[i].Vhel=1960.0 & spcal[i].Sdv=29.1 & spcal[i].W20=112.0 & spcal[i].W50=112.0 & i+=1 spcal[i].name='UGC7047' & spcal[i].Vhel=210.0 & spcal[i].Sdv=35.1 & spcal[i].W20=52.0 & spcal[i].W50=52.0 & i+=1 spcal[i].name='UGC7047a' & spcal[i].Vhel=210.0 & spcal[i].Sdv=35.1 & spcal[i].W20=52.0 & spcal[i].W50=52.0 & i+=1 spcal[i].name='UGC7047b' & spcal[i].Vhel=210.0 & spcal[i].Sdv=35.1 & spcal[i].W20=52.0 & spcal[i].W50=52.0 & i+=1 spcal[i].name='UGC7047c' & spcal[i].Vhel=210.0 & spcal[i].Sdv=35.1 & spcal[i].W20=52.0 & spcal[i].W50=52.0 & i+=1 spcal[i].name='UGC7985' & spcal[i].Vhel=652.0 & spcal[i].Sdv=51.8 & spcal[i].W20=165.0 & spcal[i].W50=165.0 & i+=1 spcal[i].name='UGC7985a' & spcal[i].Vhel=652.0 & spcal[i].Sdv=51.8 & spcal[i].W20=165.0 & spcal[i].W50=165.0 & i+=1 spcal[i].name='UGC7985b' & spcal[i].Vhel=652.0 & spcal[i].Sdv=51.8 & spcal[i].W20=165.0 & spcal[i].W50=165.0 & i+=1 spcal[i].name='UGC7985c' & spcal[i].Vhel=652.0 & spcal[i].Sdv=51.8 & spcal[i].W20=165.0 & spcal[i].W50=165.0 & i+=1 spcal[i].name='UGC10445' & spcal[i].Vhel=963.0 & spcal[i].Sdv=27.9 & spcal[i].W20=146.0 & spcal[i].W50=146.0 & i+=1 spcal[i].name='UGC10445a' & spcal[i].Vhel=963.0 & spcal[i].Sdv=27.9 & spcal[i].W20=146.0 & spcal[i].W50=146.0 & i+=1 spcal[i].name='UGC10445b' & spcal[i].Vhel=963.0 & spcal[i].Sdv=27.9 & spcal[i].W20=146.0 & spcal[i].W50=146.0 & i+=1 spcal[i].name='UGC10445c' & spcal[i].Vhel=963.0 & spcal[i].Sdv=27.9 & spcal[i].W20=146.0 & spcal[i].W50=146.0 & i+=1 spcal[i].name='UGC10521' & spcal[i].Vhel=852.0 & spcal[i].Sdv=30.4 & spcal[i].W20=225.0 & spcal[i].W50=225.0 & i+=1 spcal[i].name='UGC10521a' & spcal[i].Vhel=852.0 & spcal[i].Sdv=30.4 & spcal[i].W20=225.0 & spcal[i].W50=225.0 & i+=1 spcal[i].name='UGC10521b' & spcal[i].Vhel=852.0 & spcal[i].Sdv=30.4 & spcal[i].W20=225.0 & spcal[i].W50=225.0 & i+=1 spcal[i].name='UGC10521c' & spcal[i].Vhel=852.0 & spcal[i].Sdv=30.4 & spcal[i].W20=225.0 & spcal[i].W50=225.0 & i+=1 spcal[i].name='UGC11218' & spcal[i].Vhel=1484.0 & spcal[i].Sdv=34.7 & spcal[i].W20=329.0 & spcal[i].W50=329.0 & i+=1 spcal[i].name='UGC11218a' & spcal[i].Vhel=1484.0 & spcal[i].Sdv=34.7 & spcal[i].W20=329.0 & spcal[i].W50=329.0 & i+=1 spcal[i].name='UGC11218b' & spcal[i].Vhel=1484.0 & spcal[i].Sdv=34.7 & spcal[i].W20=329.0 & spcal[i].W50=329.0 & i+=1 spcal[i].name='UGC11218c' & spcal[i].Vhel=1484.0 & spcal[i].Sdv=34.7 & spcal[i].W20=329.0 & spcal[i].W50=329.0 & i+=1 ;; for flux calibrators observed in spectral line mode, take flux at ;; freq=1420.4 from ottflux.pro spcal[i].name='3C48Line' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C48',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C48LineA' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C48',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C48LineB' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C48',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C48LineC' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C48',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C161Line' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C161',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C161LineA' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C161',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C161LineB' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C161',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C161LineC' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C161',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C295Line' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C295',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C295LineA' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C295',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C295LineB' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C295',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C295LineC' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C295',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C286Line' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C286',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C286LineA' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C286',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C286LineB' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C286',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C286LineC' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C286',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C249.1Line' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C249.1',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C249.1LineA' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C249.1',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C249.1LineB' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C249.1',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 spcal[i].name='3C249.1LineC' & spcal[i].Vhel=1.0 & spcal[i].Sdv=ottflux('3C249.1',1420.4) & spcal[i].W20=1.0 & spcal[i].W50=1.0 & i+=1 answer=[-1.0,-1.0,-1.0] index = -1 for i=0,(n_elements(spcal)-1) do begin if source eq spcal[i].name then begin index = i break endif end if index eq -1 then begin print,'That source is not found. Valid sources are:' print, spcal.name return, answer endif answer[0]=spcal[index].Sdv answer[1]=spcal[index].W50 answer[2]=spcal[index].Vhel W20=spcal[index].W20 return, answer end