Mesoscale height and temperature fields can be extracted from the observed wind field by making use of the full divergence equation. Mass changes associated with irrotational ageostrophic motions are retained for a nearly complete description of the height field. Above the boundary layer, in the absence of friction, the divergence equation includes terms composed of the components of the wind and a Laplacian of the geopotential height field. Once the mass field is determined, the thermal structure is obtained through application of the hypsometric equation. In this paper an error analysis of this divergence method is undertaken to estimate the potential magnitude of errors associated with random errors in the wind data. Previous applications of the divergence method have been refined in the following ways. (i) The domain over which the method is applied is expanded to encompass the entire STORM-FEST domain. (ii) Wind data from 23 profiler and 38 rawinsonde sites are combined in the analysis. (iii) Observed profiler and rawinsonde data are interpolated to grid points through a modified objective analysis, and (iv) the variation in elevation of the profiler sites is taken into account. The results of the application of the divergence method to the combined wind data from profiler and rawinsonde sites show good agreement between the retrieved heights and temperatures and the observed values at rawinsonde sites. Standard deviations of the difference between the retrieved and observed data lie well within the precision of the rawinsonde instruments. The difference field shows features whose magnitude is significantly larger than the errors predicted by the error analysis, and these features are systematic rather than random in nature, suggesting that the retrieved fields are able to resolve mesoscale signatures not fully captured by the rawinsonde data alone. The divergence method is also applied solely to the profiler data to demonstrate the potential of the divergence method to provide mass and thermal fields on a routine basis at synoptic times when operational rawinsonde data are not available. A comparison of the heights derived from the profiler winds with those independently measured by rawinsondes indicates that valuable information on the evolution of atmospheric height and temperature fields can be retrieved between conventional rawinsonde release times through application of the divergence method. The implications of the results for applications of the method in weather analysis and in numerical weather prediction are discussed.