examples/H5Part/...

- bugfixes
This commit is contained in:
2015-09-09 18:25:27 +02:00
parent 106be817d9
commit 4d9d3dc879
3 changed files with 5 additions and 6 deletions
+1 -1
View File
@@ -27,7 +27,7 @@ program read_setviewf
! initialize MPI & H5hut
comm = MPI_COMM_WORLD
call mpi_init (mpi_error)
call mpi_init (mpi_ierror)
call mpi_comm_size (comm, comm_size, mpi_ierror)
call mpi_comm_rank (comm, comm_rank, mpi_ierror)
call h5_abort_on_error ()
-1
View File
@@ -71,4 +71,3 @@ main (
H5CloseFile(file);
return MPI_Finalize();
}
+4 -4
View File
@@ -17,7 +17,7 @@ program write_setview
character (len=*), parameter :: fname = "example_setview.h5"
! H5hut verbosity level
integer*8, parameter :: h5_verbosity = H5_VERBOSE_DEFAULT
integer*8, parameter :: h5_verbosity = H5_DEBUG_ALL
! we are going to write multiple consecutive blocks
integer*8, parameter :: num_blocks = 4;
@@ -44,16 +44,16 @@ program write_setview
! with H5PartSetview(). Otherwise we have to define the total number
! of particles with H5PartSetNumParticles().
offset = comm_rank * num_blocks * num_particles_per_block+1
h5_ierror = h5pt_setview (file, offset, offset + num_blocks*num_particles_per_block)
h5_ierror = h5pt_setview (file, offset, offset + num_blocks*num_particles_per_block - 1)
! write multiple consecutive blocks
allocate (data (num_particles_per_block))
do i = 1, num_blocks
! create fake data
do j = 1, num_particles_per_block
data (i) = int((j-1) + i*num_particles_per_block + offset)
data (j) = int((j-1) + (i-1)*num_particles_per_block + offset - 1)
end do
h5_ierror = h5pt_setview (file, offset + i*num_particles_per_block, (i+1)*num_particles_per_block)
h5_ierror = h5pt_setview (file, offset + (i-1)*num_particles_per_block, offset - 1 + i*num_particles_per_block)
! write data
h5_ierror = h5pt_writedata_i4 (file, "data", data)
end do