Replaced rhoMax with max over all mpi procs in torus_2d/problem.c