size = comm.Get_size() ## Writing print "\n>>> Writing ... (rank = %d)\n" % rank config = "config_mpi.xml" if len(sys.argv) > 1: config = sys.argv[1] ad.init(config, comm) fd = ad.open("temperature", "adios_test_mpi.bp", "w", comm) NX = 10 groupsize = 4 + 4 + 4 + 8 * 1 * NX t = np.array(range(NX), dtype=np.float64) + rank*NX ad.set_group_size(fd, groupsize) ad.write_int(fd, "NX", NX) ad.write_int(fd, "rank", rank) ad.write_int(fd, "size", size) ad.write(fd, "temperature", t) ad.close(fd) ad.finalize() ## Reading if rank == 0: print "\n>>> Reading ...\n" f = ad.file("adios_test_mpi.bp", comm=MPI.COMM_SELF) f.printself()