On 04/08/2014 04:01 PM, seismo_phil wrote:

Hi,

I have another question relating to the thermal diffusion problem I am working on. There is a background heat production distributed throughout the volume, but the heat production for any given cell will depend on which region the cell is in, where region is determined by cell group. I.e., all cells of group 1 are defined to comprise region 1, and they have a heat production that is different from that in region 2.

So far the examples I've seen have a load that depends only on coordinate. Is there any example in which the load depends on the region?

You can do it in the exactly same way as you define the diffusion coefficient. The 'load.f' is a material parameter just like 'coef.val', so use a function similar to get_coef().

Best regards, r.