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
|
program sds_vrs_coordvar
implicit none
C
C Parameter declaration.
C
character*7 FILE_NAME
parameter (FILE_NAME = 'SDS.hdf')
integer DFACC_READ, DFNT_INT32
parameter (DFACC_READ = 1,
+ DFNT_INT32 = 24)
integer MAX_VAR_DIMS
parameter (MAX_VAR_DIMS = 32)
C
C Function declaration.
C
integer sfstart, sfselect, sfiscvar, sffinfo, sfginfo
integer sfendacc, sfend
C
C**** Variable declaration *******************************************
C
integer sd_id, sds_id, sds_index, status
integer rank, data_type
integer n_datasets, n_file_attrs, n_attrs
integer dim_sizes(MAX_VAR_DIMS)
character*256 sds_name
C
C**** End of variable declaration ************************************
C
C
C Open the file and initialize the SD interface.
C
sd_id = sfstart(FILE_NAME, DFACC_READ)
C
C Obtain information about the file.
C
status = sffinfo(sd_id, n_datasets, n_file_attrs)
C
C Get information about each SDS in the file.
C Check whether it is a coordinate variable, then display retrieved
C information.
C Output displayed:
C
C SDS array with the name SDStemplate
C Coordinate variable with the name X_Axis
C Coordinate variable with the name Y_Axis
C
do 10 sds_index = 0, n_datasets-1
sds_id = sfselect(sd_id, sds_index)
status = sfginfo(sds_id, sds_name, rank, dim_sizes,
+ data_type, n_attrs)
status = sfiscvar(sds_id)
if (status .eq. 1) then
write(*,*) "Coordinate variable with the name ",
+ sds_name(1:6)
else
write(*,*) "SDS array with the name ",
+ sds_name(1:11)
endif
C
C Terminate access to the data set.
C
status = sfendacc(sds_id)
10 continue
C
C Terminate access to the SD interface and close the file.
C
status = sfend(sd_id)
end
|
