itemp = self.macstemp_typ
for i=0L,nen-1 do begin
    qty[*,*,i] = qty[*,*,i]*data_temperature[i,itemp]
    err[*,*,i] = err[*,*,i]*data_temperature[i,itemp]
endfor
   
data_temperature[*,itemp] = en/(kb*data_temperature[*,itemp])
