sll_k_hdf5_ser_read_attribute.F90

Go to the documentation of this file.

integer(hid_t)   ::   loc_id
integer(hid_t)   ::  attr_id
integer(hsize_t) :: dims(1)
type(h5o_info_t) :: object_info
 
dims = int([1], hsize_t)
 
! Inquire about type of object in target path: group or dataset?
call h5oget_info_by_name_f(handle%file_id, objpath, object_info, error)
sll_assert_always(error == 0)
 
! Open existing group or dataset
select case (object_info%type)
case (h5o_type_group)
   call h5gopen_f(handle%file_id, objpath, loc_id, error)
   sll_assert_always(error == 0)
case (h5o_type_dataset)
   call h5dopen_f(handle%file_id, objpath, loc_id, error)
   sll_assert_always(error == 0)
end select
 
! Open existing scalar attribute
call h5aopen_f(loc_id, attrname, attr_id, error)
sll_assert_always(error == 0)
 
! Read attribute
call h5aread_f(attr_id, datatype, attrvalue, dims, error)
sll_assert_always(error == 0)
 
! Close attribute
call h5aclose_f(attr_id, error)
sll_assert_always(error == 0)
 
! Close existing group or dataset
select case (object_info%type)
case (h5o_type_group)
   call h5gclose_f(loc_id, error)
   sll_assert_always(error == 0)
case (h5o_type_dataset)
   call h5dclose_f(loc_id, error)
   sll_assert_always(error == 0)
end select