$solvethermalddcc enable the thermal solver to use the heat source from the DDCC solver. The heat source in each elemenent will be
heatsource=E→⋅J