1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
|
!************************************************************
!
! This example shows how to write a C string to a dataset
! and read it back as a Fortran string.
! The program first writes the C
! strings to a dataset with a dataspace of DIM0, then closes the file.
! Next, it reopens the file, reads back the data into a Fortran
! fixed character string, and outputs it to the screen.
!
! This file is intended for use with HDF5 Library version 1.8
! with --enable-fortran2003
!
! ************************************************************/
PROGRAM main
USE hdf5
USE ISO_C_BINDING
IMPLICIT NONE
CHARACTER(LEN=21), PARAMETER :: filename = "h5ex_t_stringC_F03.h5"
CHARACTER(LEN=3) , PARAMETER :: dataset = "DS1"
INTEGER , PARAMETER :: dim0 = 4
INTEGER(SIZE_T) , PARAMETER :: sdim = 7
INTEGER(HID_T) :: file, filetype, memtype, space, dset ! Handles
INTEGER :: hdferr
INTEGER(HSIZE_T), DIMENSION(1:1) :: dims = (/dim0/)
INTEGER(HSIZE_T), DIMENSION(1:1) :: maxdims
CHARACTER(LEN=sdim), DIMENSION(1:dim0), TARGET :: &
wdata = (/"Parting", "is such", "sweet ", "sorrow."/)
CHARACTER(LEN=sdim), DIMENSION(:), ALLOCATABLE, TARGET :: rdata
INTEGER :: i
INTEGER(SIZE_T) :: size
TYPE(c_ptr) :: f_ptr
!
! Initialize FORTRAN interface.
!
CALL h5open_f(hdferr)
!
! Create a new file using the default properties.
!
CALL h5fcreate_f(filename, H5F_ACC_TRUNC_F, file, hdferr)
!
! Create file and memory datatypes. For this example we will save
! the strings as C strings
!
! Include the NULL TERMINATION of string in C (i.e. add +1 to the length)
CALL H5Tcopy_f(H5T_C_S1, filetype, hdferr)
CALL H5Tset_size_f(filetype, sdim+1, hdferr)
CALL H5Tcopy_f( H5T_FORTRAN_S1, memtype, hdferr)
CALL H5Tset_size_f(memtype, sdim, hdferr)
!
! Create dataspace.
!
CALL h5screate_simple_f(1, dims, space, hdferr)
!
! Create the dataset and write the string data to it.
!
CALL h5dcreate_f(file, dataset, filetype, space, dset, hdferr)
f_ptr = C_LOC(wdata(1)(1:1))
CALL H5Dwrite_f(dset, memtype, f_ptr, hdferr)
!
! Close and release resources.
!
CALL h5dclose_f(dset , hdferr)
CALL h5sclose_f(space, hdferr)
CALL H5Tclose_f(filetype, hdferr)
CALL H5Tclose_f(memtype, hdferr)
CALL h5fclose_f(file , hdferr)
!
! Now we begin the read section of this example.
!
! Open file and dataset.
!
CALL h5fopen_f(filename, H5F_ACC_RDONLY_F, file, hdferr)
CALL h5dopen_f(file, dataset, dset, hdferr)
!
! Get the datatype and its size.
!
CALL H5Dget_type_f(dset, filetype, hdferr)
CALL H5Tget_size_f(filetype, size, hdferr)
! Make sure the declared length is large enough,
! the C string contains the null character.
IF(size.GT.sdim+1)THEN
PRINT*,'ERROR: Character LEN is to small'
STOP
ENDIF
!
! Get dataspace.
!
CALL H5Dget_space_f(dset, space, hdferr)
CALL H5Sget_simple_extent_dims_f(space, dims, maxdims, hdferr)
ALLOCATE(rdata(1:dims(1)))
!
! Create the memory datatype.
!
CALL H5Tcopy_f(H5T_FORTRAN_S1, memtype, hdferr)
CALL H5Tset_size_f(memtype, sdim, hdferr)
!
! Read the data.
!
f_ptr = C_LOC(rdata(1)(1:1))
CALL H5Dread_f(dset, memtype, f_ptr, hdferr, space)
!
! Output the data to the screen.
!
DO i = 1, dims(1)
WRITE(*,'(A,"(",I0,"): ", A)') DATASET, i, rdata(i)
END DO
!
! Close and release resources.
!
CALL H5Dclose_f(dset, hdferr)
CALL H5Sclose_f(space, hdferr)
CALL H5Tclose_f(filetype, hdferr)
CALL H5Tclose_f(memtype, hdferr)
CALL H5Fclose_f(file, hdferr)
DEALLOCATE(rdata)
END PROGRAM main
|
