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
|
program definegdflds
integer status, gddeffld, gdsetfill, gdwrmeta
integer gddetach, gdclose
integer*4 gdfid, gdid1, gdid2, gdopen, gdattach
real*4 fillval1, fillval2
integer DFACC_RDWR
parameter (DFACC_RDWR=3)
integer DFNT_FLOAT32
parameter (DFNT_FLOAT32=5)
integer DFNT_FLOAT64
parameter (DFNT_FLOAT64=6)
integer HDFE_NOMERGE
parameter (HDFE_NOMERGE=0)
integer HDFE_AUTOMERGE
parameter (HDFE_AUTOMERGE=1)
fillval1=-7.0
fillval2=-9999.0
gdfid = gdopen("GridFile_created_with_hadeos_sample_file_write"//
1"r_of_HDFEOS2_version_219_or_higher_release.hdf", DFACC_RDWR)
gdid1 = gdattach(gdfid, "UTMGrid")
status = gddeffld(gdid1, "Pollution", "XDim,YDim,Time",
1 DFNT_FLOAT32, HDFE_NOMERGE)
status = gddeffld(gdid1, "Vegetation", "XDim,YDim",
1 DFNT_FLOAT32, HDFE_NOMERGE)
status = gdwrmeta(gdid1, "Extern", "XDim,YDim",
1 DFNT_FLOAT32)
status = gdsetfill(gdid1, "Pollution", fillval1)
1
gdid2 = gdattach(gdfid, "PolarGrid")
status = gddeffld(gdid2, "Temperature", "XDim,YDim",
1 DFNT_FLOAT32, HDFE_AUTOMERGE)
status = gddeffld(gdid2, "Pressure", "XDim,YDim",
1 DFNT_FLOAT32, HDFE_AUTOMERGE)
status = gddeffld(gdid2, "Soil Dryness", "XDim,YDim",
1 DFNT_FLOAT32, HDFE_NOMERGE)
status = gddeffld(gdid2, "Spectra", "XDim,YDim,Bands",
1 DFNT_FLOAT64, HDFE_AUTOMERGE)
status = gdsetfill(gdid2, "Pressure", fillval2)
status = gddetach(gdid1)
status = gddetach(gdid2)
status = gdclose(gdfid)
stop
end
|
