An array of nine square heated elements mounted in a square three-by-three pattern with no gap between the elements on a large vertical adiabatic surface with natural convective flow over the elements has been considered. Each of the elements has a uniform heat flux over its surface, the heat fluxes over eight of the elements being the same and the heat flux over the ninth element being higher than that over the other eight elements. The basic aim of the study was to determine the effect the position of the higher heat flux element on the mean temperatures of the other eight elements. The situation considered is an approximate model of situations that can arise in electronic cooling. The flow has been assumed to be steady and laminar and it has been assumed that the fluid properties are constant except for the density change with temperature which gives rise to the buoyancy forces, this having been treated by using the Boussinesq approach. The solution has been obtained by numerically solving the full three-dimensional form of the governing equations, these equations being written in terms of dimensionless variables using the commercial cfd code FLUENT. The solution has the heat flux Rayleigh number, the Prandtl number, the ratio of the heat flux over the high heat flux element to the heat flux over the other eight elements, and the position of the high heat flux element as parameters. Because of the application that motivated this work results have only been obtained for Pr = 0.7. Results have been obtained for a wide range of values of the other input parameters and the effect of these parameter values on the mean surface temperatures of each of the elements has been studied.