$usenonRAugRasheat

For the thermal solver, the major heat source is the $\vec{E} \cdot \vec{J}$. However, if we consider the nonradiative recombination as the heat source, we need to use the command of $usenonRAugRasheat to activate this function.