FUNCTION inter,wave,dn,temp if wave eq 1600 then begin data=bb_1600_phot() ;calls program to get temp/intrument reading data endif else begin data=bb_1700_phot() endelse data_x=data(*,1) data_y=data(*,0) n=n_elements(dn) temp=fltarr(n) temp(*)=spline(data_x,data_y,dn(*)) return,temp end