mirror of
https://github.com/HDFGroup/hdf5.git
synced 2024-11-21 01:04:10 +08:00
8055378bce
Adding the Fortran interface to the HDF5 library Description: Fortran is now a subdirectory of the HDF5 library tree. Platforms tested: Solaris and IRIX (O2K)
283 lines
7.5 KiB
Fortran
283 lines
7.5 KiB
Fortran
!
|
|
! 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(HSSIZE_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
|
|
LOGICAL :: status
|
|
|
|
|
|
!
|
|
! 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 predefined datatypes
|
|
!
|
|
CALL h5init_types_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.
|
|
!
|
|
CALL h5dwrite_f(dset1_id, H5T_NATIVE_INTEGER, buf1, error)
|
|
|
|
CALL h5dwrite_f(dset2_id, H5T_NATIVE_INTEGER, buf2, 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.
|
|
!
|
|
CALL H5dwrite_f(dset1_id, H5T_NATIVE_INTEGER, val, 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.
|
|
!
|
|
CALL H5dwrite_f(dset2_id, H5T_NATIVE_INTEGER, val, 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.
|
|
!
|
|
CALL h5dread_f(dset1_id, H5T_NATIVE_INTEGER, bufnew, 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, 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 predefined datatypes.
|
|
!
|
|
CALL h5close_types_f(error)
|
|
|
|
END PROGRAM SELECTEXAMPLE
|