/usr/share/doc/libhdf4-doc/examples/fortran/VD_get_vdata_info.f is in libhdf4-doc 4.2.13-2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 | program vdata_info
implicit none
C
C Parameter declaration
C
character*18 FILE_NAME
integer DFACC_READ, FULL_INTERLACE
integer FIELD_SIZE
C
parameter (FILE_NAME = 'General_Vdatas.hdf',
+ DFACC_READ = 1,
+ FULL_INTERLACE = 0,
+ FIELD_SIZE = 80)
C
C Function declaration
C
integer hopen, hclose
integer vfstart, vsfatch, vsfgid, vsfinq,
+ vsfisat, vsfdtch, vfend
C
C**** Variable declaration *******************************************
C
integer status
integer file_id, vdata_id, vdata_ref
integer n_records, interlace_mode, vdata_size
character*64 vdata_name
character*80 fieldname_list
C
C**** End of variable declaration ************************************
C
C
C Open the HDF file for reading.
C
file_id = hopen(FILE_NAME, DFACC_READ, 0)
C
C Initialize the VS interface.
C
status = vfstart(file_id)
C
C Set the reference number to -1 to start the search from the beginning
C of the file.
C
vdata_ref = -1
10 continue
C
C Use vsfgid to obtain each vdata by its reference number then
C attach to the vdata and get information. The loop terminates
C when the last vdata is reached.
C
vdata_ref = vsfgid(file_id, vdata_ref)
if (vdata_ref .eq. -1) goto 100
C
C Attach to the current vdata for reading.
C
vdata_id = vsfatch(file_id, vdata_ref, 'r')
C
C Test whether the current vdata is not a storage for an attribute,
C then obtain and display its information.
if (vsfisat(vdata_id) .ne. 1) then
C Initialize buffers before getting values back.
vdata_name = ' '
fieldname_list = ' '
n_records = -1
vdata_size = -1
status = vsfinq(vdata_id, n_records, interlace_mode,
+ fieldname_list, vdata_size, vdata_name)
if (status .eq. 0) then
write(*,*) 'Vdata: ', vdata_name
write(*,*) 'contains ', n_records, ' records'
if (interlace_mode .eq. 0) then
write(*,*) 'Interlace mode: FULL'
else
write(*,*) 'Interlace mode: NONE'
endif
write(*,*) 'Fields: ', fieldname_list(1:30)
write(*,*) 'Vdata record size in bytes :', vdata_size
write(*,*)
endif
endif
C
C Detach from the current vdata.
C
status = vsfdtch(vdata_id)
goto 10
100 continue
C
C Terminate access to the vdata and to the VS interface, and
C close the HDF file.
C
status = vsfdtch(vdata_id)
status = vfend(file_id)
status = hclose(file_id)
end
|