$heatiterative is the command to ask the program solve thermal equation and Poisson equation self-consistenly.