Heat flux boundary condition

Hi all,

I think it can be really useful to implement the input of the constant heat flux boundary condition, either in W or W/m2, since conductvity will vary through the simulation and it will facilitate its input.

Thanks!