function [b_temp]=radtot(freq,radiance);
c1=1.1911E-8;
c2=1.439;
b_temp=c2.*freq./(log(1+(c1.*freq.*freq.*freq)./radiance));
