diff options
Diffstat (limited to 'doc/html/Tutor/examples/selectele.f90')
| -rw-r--r-- | doc/html/Tutor/examples/selectele.f90 | 282 |
1 files changed, 0 insertions, 282 deletions
diff --git a/doc/html/Tutor/examples/selectele.f90 b/doc/html/Tutor/examples/selectele.f90 deleted file mode 100644 index 8727bd9..0000000 --- a/doc/html/Tutor/examples/selectele.f90 +++ /dev/null @@ -1,282 +0,0 @@ -! -! 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 - 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 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. - ! - 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_f(error) - - END PROGRAM SELECTEXAMPLE |