crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 06:05:58 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

successfully wrote example dataset 2D array to HDF5 output file

Browse source
This commit is contained in:
Bryan Herman 2012-01-05 14:27:13 -08:00
parent f032d4b9ef
commit 67fc7ce085
1 changed files with 30 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

31
src/cmfd_utils.F90
View file

@ -519,9 +519,23 @@ contains
character(LEN=8), parameter :: filename = "filef.h5" ! File name
character(LEN=4), parameter :: grpname = "cmfd" ! Group name
character(LEN=9), parameter :: dsetname = "cmfd/dset" ! dataset name
integer(HID_T) :: file_id ! File identifier
integer(HID_T) :: group_id ! Group identifieer
integer(HID_T) :: group_id ! Group identifier
integer(HID_T) :: dataspace_id ! Data space identifier
integer(HID_T) :: dataset_id ! Dataset identifier
integer :: error ! Error flag
integer :: i,j
integer, dimension(3,3) :: dset_data
integer(HSIZE_T), dimension(2) :: dims = (/3,3/)
integer :: rank = 2
! Create example dataset
do i = 1,3
do j = 1,3
dset_data(i,j) = j
end do
end do
! Initialize FORTRAN interface.
call h5open_f(error)
@ -532,6 +546,21 @@ contains
! Create the CMFD group
call h5gcreate_f(file_id, grpname, group_id, error)
! Create the data space for the dataset
call h5screate_simple_f(rank, dims, dataspace_id, error)
! Create the data set in CMFD group
call h5dcreate_f(file_id, dsetname, H5T_NATIVE_INTEGER,dataspace_id,dataset_id,error)
! Write the dataset
call h5dwrite_f(dataset_id, H5T_NATIVE_INTEGER, dset_data, dims, error)
! Close the data space
call h5sclose_f(dataspace_id, error)
! Close the dataset
call h5dclose_f(dataset_id, error)
! Close the CMFD group
call h5gclose_f(group_id, error)