% Planck's spectral hemispherical blackbody emission
function Eb = Eblkl(lrow,Tcol)
  C1=0.59552197e8; C2=14387.69;
  lT=bsxfun(@times,max(Tcol(:),0),lrow);
  Eb=bsxfun(@rdivide,2*pi*C1./(lrow.^5),exp(C2./lT)-1.0);
end
