Added the Jacobian back into the isothermal wall bc and fixed the heat flux calculation