message: python/test: generalize set up of global array a bit
  So this would now work with other domain decompositions, too.
diff: @@ -27,14 +27,11 @@ count = [Nx, Ny]
  start = [rank * Nx, 0]
  shape = [size * Nx, Ny]
  -temperatures = np.empty(count, dtype=np.int)
  +temperatures = np.zeros(count, dtype=np.int)
  for i in range(0, Nx):
  - iGlobal = start[0] + i
  -
  for j in range(0, Ny):
  - value = iGlobal * shape[1] + j
  - temperatures[i, j] = value
  + temperatures[i, j] = (start[0] + i) * shape[1] + (j + start[1])
  # print(temperatures)
  # ADIOS2 read
  