span> (~isempty(drv.particles) &amp;&amp;<span class="keyword">...</span>
               (drv.type == drv.THERMAL || drv.type == drv.THERMO_MECHANICAL))
                tavg = mean([drv.particles.temperature]);
                tdev = std([drv.particles.temperature]);
                tmin = min([drv.particles.temperature]);
                tmax = max([drv.particles.temperature]);
                fprintf(<span class="string">'Average temperature....: %.3e\n'</span>,tavg);
                fprintf(<span class="