Question

我想在Fortran中使用HDF5 Tables，但我遇到了一些麻烦。一个例子很有用，但只有C个例子是provided。

我的部分问题是如何处理所需的偏移量和大小。使用gfortran，我可以使用sizeof和loc，但这些是gfortran特定的扩展。我看到Fortran的新版本有c_loc和c_sizeof这可能会有所帮助，但也许它们只能处理c变量？

我也不知道chunk_size应该是什么。

FWIW，下面的丑陋代码使用loc和sizeof。它编译并运行但是给出了错误（并且没有在hdf5文件中添加任何内容）：

 size:                    8
 offsets:                    0                    4
 types:           0           0
 initialized
 file open
HDF5-DIAG: Error detected in HDF5 (1.8.8) thread 140636908803840:
  #000: ../../../src/H5Tcompound.c line 370 in H5Tinsert(): not a datatype
    major: Invalid arguments to routine
    minor: Inappropriate type
 table created

所以偏移量可能有意义，但类型并不是真的。

任何帮助将不胜感激。感谢。

module tdef
  type blarg
    integer :: arg
    real    :: blah
  end type
end

PROGRAM H5_TABLE

  use tdef
  USE HDF5
  use h5lt
  use h5tb

  IMPLICIT NONE

  INTEGER(HID_T) :: file_id
  INTEGER     ::   error
  integer(HSIZE_T) :: nfields, nrecords
  integer(SIZE_T)  :: type_size
  type(blarg) :: test
  character(len=4), dimension(2) :: field_names = (/' arg', 'blah'/)
  integer(SIZE_T), dimension(2) :: field_offset
  integer(HID_T), dimension(2) :: field_types
  integer(HSIZE_T) ::  chunk_size = 1
  integer :: compress = 0

  nfields  = 2
  nrecords = 2

  type_size = sizeof(test)
  print *, "size:", type_size

  field_offset(1) = loc(test%arg) - loc(test)
  field_offset(2) = loc(test%blah) - loc(test)
  print *, "offsets:", field_offset

  field_types(1) = H5T_NATIVE_INTEGER
  field_types(2) = H5T_NATIVE_REAL
  print *, "types:", field_types

  CALL h5open_f(error)
  print *, "initialized"

  CALL h5fcreate_f("cpt.h5", H5F_ACC_TRUNC_F, file_id, error)
  print *, "file open"

  call h5tbmake_table_f('Table Title', file_id, "Steps", nfields, &
            nrecords, type_size, field_names, field_offset, &
            field_types, chunk_size, compress, error) 
  print *, "table created"

  CALL h5fclose_f(file_id, error)

  CALL h5close_f(error)

END PROGRAM H5_TABLE

Answer 1

事实证明这是一个非常简单的问题。

使用H5T_NATIVE_INTEGER和H5T_NATIVE_REAL后调用h5open_f。但是，在调用h5open_f之前，不会为这两个值分配适当的值。调用h5open_f首先解决了我的问题。

如何在Fortran中使用HDF5表

1 个答案: