DFLUX subroutine

Hi all, I am looking for a DFLUX subroutine that incorporates some sort of control algorithm. I am trying to apply a heat flux through an embedded heater layer to a set temperature.