FUNCTION bb_1600_phot,t_v_photnum h=6.62606957e-34 c=float(299792458) k=double(1.3806488e-23) t_v_photnum=fltarr(36,2) resp_struct=aia_get_response(/uv,/dn,/full) ;instrument response function structure inst_resp=resp_struct.a1600_full.effarea ;instrument response at 1600A, units of cm^2 DN phot^-1 sa_a=resp_struct.a1600_full.platescale ;size of pixel times solid angle for detector resp=fltarr(1001) j=0 for t=4000,7500,100 do begin spect_tot=bb_specf(t) ;calls function to find BB spectrum, puts into units of W cm^-2 A^-1 st^-1 spect_inc=spect_tot*0 spect_inc(*,1)=spect_tot(*,1)*sa_a*1e-14 ;multiplies, finds amt incident on sensor spect_inc(*,0)=spect_tot(*,0) spect_phot=spect_inc*0 ;set up array for photon spectrum for i=0,1000 do begin spect_phot(i,0)=i+1000 ;makes 0 index (wavelength) 1000-2000 spect_phot(i,1)=spect_inc(i,1)/(h*c/((i+1000)*1e-10)) ;finds # of photons at each wavelength if i lt 200 then begin ;lines up indices of response function (0th entry is 1000A, instead of 1200A) resp(i)=0 endif else begin resp(i)=inst_resp(i-200) endelse endfor spect_act=spect_phot(*,1) spect_act=spect_act*resp t_v_photnum(j,0)=t t_v_photnum(j,1)=int_tabulated(spect_phot(*,0),spect_act)/2.1 j=j+1 endfor return,t_v_photnum end