;Compute some data from the BASE observations restore,'BASEidl.dat' bins[0]=3.5 dt=fltarr(11)*0-999. for i=0,10 do begin srch=1 bb=n_elements(bins)-1 nn=0 while srch and bb ge 0 do begin nn0=nn nn=total(npb[i,bb:*]*3.0) ;3 um bin widths if nn gt 3.0 then srch=0 else bb=bb-1 endwhile if srch eq 0 then begin if bb lt n_elements(bins)-1 then dt[i]=(bins[bb+1]+(bins[bb]-bins[bb+1])*(3.0-nn0)/(nn-nn0)) else dt[i]=bins[bb] endif endfor end