mirror of
https://github.com/HDFGroup/hdf5.git
synced 2025-01-06 14:56:51 +08:00
89fbe00dec
* commit '54957d37f5aa73912763dbb6e308555e863c43f4': Commit copyright header change for src/H5PLpkg.c which was added after running script to make changes. Add new files in release_docs to MANIFEST. Cimmit changes to Makefile.in(s) and H5PL.c that resulted from running autogen.sh. Merge pull request #407 in HDFFV/hdf5 from ~LRKNOX/hdf5_lrk:hdf5_1_10_1 to hdf5_1_10_1 Change copyright headers to replace url referring to file to be removed and replace it with new url for COPYING file.
302 lines
8.3 KiB
Fortran
302 lines
8.3 KiB
Fortran
! * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
|
|
! Copyright by The HDF Group. *
|
|
! Copyright by the Board of Trustees of the University of Illinois. *
|
|
! All rights reserved. *
|
|
! *
|
|
! This file is part of HDF5. The full HDF5 copyright notice, including *
|
|
! terms governing use, modification, and redistribution, is contained in *
|
|
! the COPYING file, which can be found at the root of the source code *
|
|
! distribution tree, or in https://support.hdfgroup.org/ftp/HDF5/releases. *
|
|
! If you do not have access to either file, you may request a copy from *
|
|
! help@hdfgroup.org. *
|
|
! * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
|
|
!
|
|
!
|
|
! This program creates two files, copy1.h5, and copy2.h5.
|
|
! In copy1.h5, it creates a 3x4 dataset called 'Copy1',
|
|
! and write 0's to this dataset.
|
|
! In copy2.h5, it create a 3x4 dataset called 'Copy2',
|
|
! and write 1's to this dataset.
|
|
! It closes both files, reopens both files, selects two
|
|
! points in copy1.h5 and writes values to them. Then it
|
|
! uses an H5Scopy to write the same selection to copy2.h5.
|
|
! Program reopens the files, and reads and prints the contents of
|
|
! the two datasets.
|
|
!
|
|
|
|
PROGRAM SELECTEXAMPLE
|
|
|
|
USE HDF5 ! This module contains all necessary modules
|
|
|
|
IMPLICIT NONE
|
|
|
|
CHARACTER(LEN=8), PARAMETER :: filename1 = "copy1.h5" ! File name
|
|
CHARACTER(LEN=8), PARAMETER :: filename2 = "copy2.h5" !
|
|
CHARACTER(LEN=5), PARAMETER :: dsetname1 = "Copy1" ! Dataset name
|
|
CHARACTER(LEN=5), PARAMETER :: dsetname2 = "Copy2" !
|
|
|
|
INTEGER, PARAMETER :: RANK = 2 ! Dataset rank
|
|
|
|
INTEGER(SIZE_T), PARAMETER :: NUMP = 2 ! Number of points selected
|
|
|
|
INTEGER(HID_T) :: file1_id ! File1 identifier
|
|
INTEGER(HID_T) :: file2_id ! File2 identifier
|
|
INTEGER(HID_T) :: dset1_id ! Dataset1 identifier
|
|
INTEGER(HID_T) :: dset2_id ! Dataset2 identifier
|
|
INTEGER(HID_T) :: dataspace1 ! Dataspace identifier
|
|
INTEGER(HID_T) :: dataspace2 ! Dataspace identifier
|
|
INTEGER(HID_T) :: memspace ! memspace identifier
|
|
|
|
INTEGER(HSIZE_T), DIMENSION(1) :: dimsm = (/2/)
|
|
! Memory dataspace dimensions
|
|
INTEGER(HSIZE_T), DIMENSION(2) :: dimsf = (/3,4/)
|
|
! File dataspace dimensions
|
|
INTEGER(HSIZE_T), DIMENSION(RANK,NUMP) :: coord ! Elements coordinates
|
|
! in the file
|
|
|
|
INTEGER, DIMENSION(3,4) :: buf1, buf2, bufnew ! Data buffers
|
|
INTEGER, DIMENSION(2) :: val = (/53, 59/) ! Values to write
|
|
|
|
INTEGER :: memrank = 1 ! Rank of the dataset in memory
|
|
|
|
INTEGER :: i, j
|
|
|
|
INTEGER :: error ! Error flag
|
|
INTEGER(HSIZE_T), DIMENSION(2) :: data_dims
|
|
|
|
|
|
!
|
|
! Create two files containing identical datasets. Write 0's to one
|
|
! and 1's to the other.
|
|
!
|
|
|
|
!
|
|
! Data initialization.
|
|
!
|
|
do i = 1, 3
|
|
do j = 1, 4
|
|
buf1(i,j) = 0;
|
|
end do
|
|
end do
|
|
|
|
do i = 1, 3
|
|
do j = 1, 4
|
|
buf2(i,j) = 1;
|
|
end do
|
|
end do
|
|
|
|
!
|
|
! Initialize FORTRAN interface.
|
|
!
|
|
CALL h5open_f(error)
|
|
|
|
!
|
|
! Create file1, file2 using default properties.
|
|
!
|
|
CALL h5fcreate_f(filename1, H5F_ACC_TRUNC_F, file1_id, error)
|
|
|
|
CALL h5fcreate_f(filename2, H5F_ACC_TRUNC_F, file2_id, error)
|
|
|
|
!
|
|
! Create the data space for the datasets.
|
|
!
|
|
CALL h5screate_simple_f(RANK, dimsf, dataspace1, error)
|
|
|
|
CALL h5screate_simple_f(RANK, dimsf, dataspace2, error)
|
|
|
|
!
|
|
! Create the datasets with default properties.
|
|
!
|
|
CALL h5dcreate_f(file1_id, dsetname1, H5T_NATIVE_INTEGER, dataspace1, &
|
|
dset1_id, error)
|
|
|
|
CALL h5dcreate_f(file2_id, dsetname2, H5T_NATIVE_INTEGER, dataspace2, &
|
|
dset2_id, error)
|
|
|
|
!
|
|
! Write the datasets.
|
|
!
|
|
data_dims(1) = 3
|
|
data_dims(2) = 4
|
|
CALL h5dwrite_f(dset1_id, H5T_NATIVE_INTEGER, buf1, data_dims, error)
|
|
|
|
CALL h5dwrite_f(dset2_id, H5T_NATIVE_INTEGER, buf2, data_dims, error)
|
|
|
|
!
|
|
! Close the dataspace for the datasets.
|
|
!
|
|
CALL h5sclose_f(dataspace1, error)
|
|
|
|
CALL h5sclose_f(dataspace2, error)
|
|
|
|
!
|
|
! Close the datasets.
|
|
!
|
|
CALL h5dclose_f(dset1_id, error)
|
|
|
|
CALL h5dclose_f(dset2_id, error)
|
|
|
|
!
|
|
! Close the files.
|
|
!
|
|
CALL h5fclose_f(file1_id, error)
|
|
|
|
CALL h5fclose_f(file2_id, error)
|
|
|
|
!
|
|
! Open the two files. Select two points in one file, write values to
|
|
! those point locations, then do H5Scopy and write the values to the
|
|
! other file. Close files.
|
|
!
|
|
|
|
!
|
|
! Open the files.
|
|
!
|
|
CALL h5fopen_f (filename1, H5F_ACC_RDWR_F, file1_id, error)
|
|
|
|
CALL h5fopen_f (filename2, H5F_ACC_RDWR_F, file2_id, error)
|
|
|
|
!
|
|
! Open the datasets.
|
|
!
|
|
CALL h5dopen_f(file1_id, dsetname1, dset1_id, error)
|
|
|
|
CALL h5dopen_f(file2_id, dsetname2, dset2_id, error)
|
|
|
|
!
|
|
! Get dataset1's dataspace identifier.
|
|
!
|
|
CALL h5dget_space_f(dset1_id, dataspace1, error)
|
|
|
|
!
|
|
! Create memory dataspace.
|
|
!
|
|
CALL h5screate_simple_f(memrank, dimsm, memspace, error)
|
|
|
|
!
|
|
! Set the selected point positions. Because Fortran array index starts
|
|
! from 1, so add one to the actual select points in C.
|
|
!
|
|
coord(1,1) = 1
|
|
coord(2,1) = 2
|
|
coord(1,2) = 1
|
|
coord(2,2) = 4
|
|
|
|
!
|
|
! Select the elements in file space.
|
|
!
|
|
CALL h5sselect_elements_f(dataspace1, H5S_SELECT_SET_F, RANK, NUMP,&
|
|
coord, error)
|
|
|
|
!
|
|
! Write value into the selected points in dataset1.
|
|
!
|
|
data_dims(1) = 2
|
|
CALL H5dwrite_f(dset1_id, H5T_NATIVE_INTEGER, val, data_dims, error, &
|
|
mem_space_id=memspace, file_space_id=dataspace1)
|
|
|
|
!
|
|
! Copy the daspace1 into dataspace2.
|
|
!
|
|
CALL h5scopy_f(dataspace1, dataspace2, error)
|
|
|
|
!
|
|
! Write value into the selected points in dataset2.
|
|
!
|
|
data_dims(1) = 2
|
|
CALL H5dwrite_f(dset2_id, H5T_NATIVE_INTEGER, val, data_dims, error, &
|
|
mem_space_id=memspace, file_space_id=dataspace2)
|
|
|
|
!
|
|
! Close the dataspace for the datasets.
|
|
!
|
|
CALL h5sclose_f(dataspace1, error)
|
|
|
|
CALL h5sclose_f(dataspace2, error)
|
|
|
|
!
|
|
! Close the memoryspace.
|
|
!
|
|
CALL h5sclose_f(memspace, error)
|
|
|
|
!
|
|
! Close the datasets.
|
|
!
|
|
CALL h5dclose_f(dset1_id, error)
|
|
|
|
CALL h5dclose_f(dset2_id, error)
|
|
|
|
!
|
|
! Close the files.
|
|
!
|
|
CALL h5fclose_f(file1_id, error)
|
|
|
|
CALL h5fclose_f(file2_id, error)
|
|
|
|
!
|
|
! Open both files and print the contents of the datasets.
|
|
!
|
|
|
|
!
|
|
! Open the files.
|
|
!
|
|
CALL h5fopen_f (filename1, H5F_ACC_RDWR_F, file1_id, error)
|
|
|
|
CALL h5fopen_f (filename2, H5F_ACC_RDWR_F, file2_id, error)
|
|
|
|
!
|
|
! Open the datasets.
|
|
!
|
|
CALL h5dopen_f(file1_id, dsetname1, dset1_id, error)
|
|
|
|
CALL h5dopen_f(file2_id, dsetname2, dset2_id, error)
|
|
|
|
!
|
|
! Read dataset from the first file.
|
|
!
|
|
data_dims(1) = 3
|
|
data_dims(2) = 4
|
|
CALL h5dread_f(dset1_id, H5T_NATIVE_INTEGER, bufnew, data_dims, error)
|
|
|
|
!
|
|
! Display the data read from dataset "Copy1"
|
|
!
|
|
write(*,*) "The data in dataset Copy1 is: "
|
|
do i = 1, 3
|
|
print *, (bufnew(i,j), j = 1,4)
|
|
end do
|
|
|
|
!
|
|
! Read dataset from the second file.
|
|
!
|
|
CALL h5dread_f(dset2_id, H5T_NATIVE_INTEGER, bufnew, data_dims, error)
|
|
|
|
!
|
|
! Display the data read from dataset "Copy2"
|
|
!
|
|
write(*,*) "The data in dataset Copy2 is: "
|
|
do i = 1, 3
|
|
print *, (bufnew(i,j), j = 1,4)
|
|
end do
|
|
|
|
!
|
|
! Close datasets.
|
|
!
|
|
CALL h5dclose_f(dset1_id, error)
|
|
|
|
CALL h5dclose_f(dset2_id, error)
|
|
|
|
!
|
|
! Close files.
|
|
!
|
|
CALL h5fclose_f(file1_id, error)
|
|
|
|
CALL h5fclose_f(file2_id, error)
|
|
|
|
!
|
|
! Close FORTRAN interface.
|
|
!
|
|
CALL h5close_f(error)
|
|
|
|
END PROGRAM SELECTEXAMPLE
|