diff --git a/matlab/profiler.m b/matlab/profiler.m
index e578665159ac924839c9c4da56c09a0eb703807c..f5a6821b1e42accd1fe0a44c4ecd3aa7c8bde2fd 100644
--- a/matlab/profiler.m
+++ b/matlab/profiler.m
@@ -10,29 +10,30 @@ adv_field_Tc = h5read(filename,'/profiler/Tc_adv_field');
 poisson_Tc   = h5read(filename,'/profiler/Tc_poisson');
 Sapj_Tc      = h5read(filename,'/profiler/Tc_Sapj');
 diag_Tc      = h5read(filename,'/profiler/Tc_diag');
+comm_Tc      = h5read(filename,'/profiler/Tc_comm');
 checkfield_Tc= h5read(filename,'/profiler/Tc_checkfield');
 step_Tc      = h5read(filename,'/profiler/Tc_step');
 Ts0D         = h5read(filename,'/profiler/time');
 
-missing_Tc   = step_Tc - rhs_Tc - adv_field_Tc -...
+missing_Tc   = step_Tc - rhs_Tc - adv_field_Tc - comm_Tc -...
                poisson_Tc - Sapj_Tc -diag_Tc -checkfield_Tc;
 total_Tc     = step_Tc;
 
-TIME_PER_FCT = [diff(rhs_Tc); diff(adv_field_Tc); diff(poisson_Tc);...
+TIME_PER_FCT = [diff(rhs_Tc); diff(adv_field_Tc); diff(poisson_Tc); diff(comm_Tc);...
     diff(Sapj_Tc); diff(checkfield_Tc); diff(diag_Tc); diff(missing_Tc)];
-TIME_PER_FCT = reshape(TIME_PER_FCT,[numel(TIME_PER_FCT)/7,7]);
+TIME_PER_FCT = reshape(TIME_PER_FCT,[numel(TIME_PER_FCT)/8,8]);
 
 TIME_PER_STEP = sum(TIME_PER_FCT,2);
 TIME_PER_CPU  = trapz(Ts0D(2:end),TIME_PER_STEP);
 
 %% Plots
-TIME_PER_FCT = [diff(rhs_Tc); diff(adv_field_Tc); diff(poisson_Tc);...
+TIME_PER_FCT = [diff(rhs_Tc); diff(adv_field_Tc); diff(poisson_Tc); diff(comm_Tc);...
     diff(Sapj_Tc); diff(checkfield_Tc); diff(diag_Tc); diff(missing_Tc)];
-TIME_PER_FCT = reshape(TIME_PER_FCT,[numel(TIME_PER_FCT)/7,7]);
+TIME_PER_FCT = reshape(TIME_PER_FCT,[numel(TIME_PER_FCT)/8,8]);
 fig = figure;
 
 p1 = area(Ts0D(2:end),TIME_PER_FCT,'LineStyle','none');
-legend('Compute RHS','Adv. fields','Poisson','Sapj','Check+Sym','Diag','Missing')
+legend('Compute RHS','Adv. fields','Poisson', 'comm', 'Sapj','Check+Sym','Diag','Missing')
 xlabel('Sim. Time [$\rho_s/c_s$]'); ylabel('Step Comp. Time [s]')
 xlim([Ts0D(2),Ts0D(end)]);
 title(sprintf('Proc. #1, total sim. time  ~%.0f [h]',CPUTIME/3600))