In the initialisation step, since we are computing only the densities, there is no need to keep the tend,rho and gpart communications; they can all be skipped.
This should be a trivial change but I have made enough mistakes lately that another pair of eyes on this would be welcomed.