bug in mpi code, output now tested