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
|
program write_data
implicit none
C
C Parameter declaration.
C
character*7 FILE_NAME
character*11 SDS_NAME
integer X_LENGTH, Y_LENGTH, RANK
parameter (FILE_NAME = 'SDS.hdf',
+ SDS_NAME = 'SDStemplate',
+ X_LENGTH = 5,
+ Y_LENGTH = 16,
+ RANK = 2)
integer DFACC_WRITE, DFNT_INT32
parameter (DFACC_WRITE = 2,
+ DFNT_INT32 = 24)
C
C Function declaration.
C
integer sfstart, sfselect, sfwdata, sfendacc, sfend
C
C**** Variable declaration *******************************************
C
integer sd_id, sds_id, sds_index, status
integer start(2), edges(2), stride(2)
integer i, j
integer data(X_LENGTH, Y_LENGTH)
C
C**** End of variable declaration ************************************
C
C
C Data set data initialization.
C
do 20 j = 1, Y_LENGTH
do 10 i = 1, X_LENGTH
data(i, j) = i + j - 1
10 continue
20 continue
C
C Open the file and initialize the SD interface.
C
sd_id = sfstart(FILE_NAME, DFACC_WRITE)
C
C Attach to the first data set.
C
sds_index = 0
sds_id = sfselect(sd_id, sds_index)
C
C Define the location and size of the data to be written
C to the data set. Note that setting values of the array stride to 1
C specifies the contiguous writing of data.
C
start(1) = 0
start(2) = 0
edges(1) = X_LENGTH
edges(2) = Y_LENGTH
stride(1) = 1
stride(2) = 1
C
C Write the stored data to the data set named in SDS_NAME.
C Note that the routine sfwdata is used instead of sfwcdata
C to write the numeric data.
C
status = sfwdata(sds_id, start, stride, edges, data)
C
C Terminate access to the data set.
C
status = sfendacc(sds_id)
C
C Terminate access to the SD interface and close the file.
C
status = sfend(sd_id)
end
|
