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
|
# Resolution and component flags
flags[1] resolutionAndComponentFlags 'grib1/7.table' : edition_specific,no_copy ;
# Note our flagbit numbers run from 7 to 0, while WMO convention uses 1 to 8
# (most significant to least significant)
flagbit ijDirectionIncrementGiven(resolutionAndComponentFlags,7) = 1 ;
# For grib 1 to 2
alias iDirectionIncrementGiven = ijDirectionIncrementGiven;
alias jDirectionIncrementGiven = ijDirectionIncrementGiven;
alias DiGiven = ijDirectionIncrementGiven;
alias DjGiven = ijDirectionIncrementGiven;
flagbit earthIsOblate(resolutionAndComponentFlags,6) : dump;
if (earthIsOblate) {
# Earth assumed oblate spheroidal with size as determined by IAU in 1965
transient earthMajorAxis = 6378160.0;
transient earthMinorAxis = 6356775.0;
alias earthMajorAxisInMetres=earthMajorAxis;
alias earthMinorAxisInMetres=earthMinorAxis;
}
flagbit resolutionAndComponentFlags3(resolutionAndComponentFlags,5) = 0: read_only;
flagbit resolutionAndComponentFlags4(resolutionAndComponentFlags,4) = 0: read_only;
flagbit uvRelativeToGrid(resolutionAndComponentFlags,3) : dump;
flagbit resolutionAndComponentFlags6(resolutionAndComponentFlags,2) = 0: read_only;
flagbit resolutionAndComponentFlags7(resolutionAndComponentFlags,1) = 0: read_only;
flagbit resolutionAndComponentFlags8(resolutionAndComponentFlags,0) = 0: read_only;
|
