Change heat transfer time step in solver to account for thermal_subcy…

…cle_steps > 1

src/CabanaPD_Solver.hpp

@@ -241,7 +241,8 @@ class SolverElastic
             {
                 if ( step % thermal_subcycle_steps == 0 )
                     computeHeatTransfer( *heat_transfer, *particles,
-                                         neigh_iter_tag{}, dt );
+                                         neigh_iter_tag{},
+                                         thermal_subcycle_steps * dt );
             }
 
             // Add non-force boundary condition.