It would be nice to have the ability to specify a spatially varying heat flux boundary condition, or in other words a heat flux profile or distribution. This could be done by defining a function of heat flux as a function of position on a surface with appropriate limits. As of right now, the only way to accomplish this is to break the surface into descrete surfaces and apply average heat fluxes to each region. This can be especially painful when you have a steep slope in the distribution of heat flux vs. position function.
This would be useful for applications where the heat flux applied to a surface is non-uniform, such as concentrated solar power applications where you might have a hot spot in the center and a profile much like a gaussian distribution.
