Abstract
A numerical method for the computation of radiative flux divergence is given, which is particularly suited for micrometeorological use. By avoiding derivatives of the emissivity for very short path lengths, this method is relatively insensitive to uncertainties in emissivity values over that range. The method permits the use of different emissivities without the need of lengthy tabulations while offering all the advantages of Brooks' method.