mirror of https://github.com/ecmwf/eccodes.git
140 lines
3.4 KiB
Modula-2
140 lines
3.4 KiB
Modula-2
# (C) Copyright 2005- ECMWF.
|
|
|
|
# moved here to allow different bitsPerValue in second order packing
|
|
unsigned[1] bitsPerValue : dump;
|
|
alias numberOfBitsContainingEachPackedValue = bitsPerValue;
|
|
alias accuracy = bitsPerValue;
|
|
|
|
# For grib1 -> grib2
|
|
#constant dataRepresentationTemplateNumber = 51;
|
|
|
|
constant PUnset = -32767;
|
|
|
|
# N is a pointer to the start of the packed data (i.e. gives octet number)
|
|
unsigned[2] N : read_only,dump;
|
|
signed[2] P = PUnset;
|
|
|
|
unsigned[1] JS=0 : dump;
|
|
unsigned[1] KS=0 : dump;
|
|
unsigned[1] MS=0 : dump;
|
|
|
|
alias subSetJ=JS;
|
|
alias subSetK=KS;
|
|
alias subSetM=MS;
|
|
|
|
constant GRIBEXShBugPresent = 1;
|
|
if (gribex_mode_on()) {
|
|
transient computeLaplacianOperator=0 : hidden;
|
|
} else {
|
|
transient computeLaplacianOperator=1 : hidden;
|
|
}
|
|
|
|
meta data.laplacianOperator scale(P,oneConstant,grib1divider,truncateLaplacian) : dump;
|
|
meta laplacianOperatorIsSet evaluate(P != PUnset && !computeLaplacianOperator);
|
|
|
|
if (localUsePresent) {
|
|
if (changed(localDefinitionNumber)) {
|
|
transient TS = 0;
|
|
meta TScalc spectral_truncation(JS,KS,MS,TS) : read_only,hidden;
|
|
meta Nassigned octet_number(N,4*TScalc) : hidden;
|
|
}
|
|
}
|
|
|
|
position offsetBeforeData;
|
|
meta values data_g1complex_packing(
|
|
section4Length,
|
|
offsetBeforeData,
|
|
offsetSection4,
|
|
unitsFactor,
|
|
unitsBias,
|
|
changingPrecision,
|
|
numberOfCodedValues,
|
|
bitsPerValue,
|
|
referenceValue,
|
|
binaryScaleFactor,
|
|
decimalScaleFactor,
|
|
optimizeScaleFactor,
|
|
|
|
GRIBEXShBugPresent,
|
|
ieeeFloats,
|
|
|
|
laplacianOperatorIsSet,
|
|
laplacianOperator,
|
|
subSetJ,
|
|
subSetK,
|
|
subSetM,
|
|
pentagonalResolutionParameterJ,
|
|
pentagonalResolutionParameterK,
|
|
pentagonalResolutionParameterM,
|
|
|
|
halfByte,
|
|
N,packingType,spectral_ieee,precision
|
|
) : dump;
|
|
|
|
meta data.packedValues data_sh_packed(
|
|
section4Length,
|
|
offsetBeforeData,
|
|
offsetSection4,
|
|
|
|
unitsFactor,
|
|
unitsBias,
|
|
changingPrecision,
|
|
numberOfCodedValues,
|
|
bitsPerValue,
|
|
referenceValue,
|
|
binaryScaleFactor,
|
|
decimalScaleFactor,
|
|
optimizeScaleFactor,
|
|
|
|
GRIBEXShBugPresent,
|
|
ieeeFloats,
|
|
|
|
laplacianOperatorIsSet,
|
|
laplacianOperator,
|
|
subSetJ,
|
|
subSetK,
|
|
subSetM,
|
|
|
|
pentagonalResolutionParameterJ,
|
|
pentagonalResolutionParameterK,
|
|
pentagonalResolutionParameterM
|
|
) : read_only;
|
|
|
|
meta data.unpackedValues data_sh_unpacked(
|
|
section4Length,
|
|
offsetBeforeData,
|
|
offsetSection4,
|
|
|
|
unitsFactor,
|
|
unitsBias,
|
|
changingPrecision,
|
|
numberOfCodedValues,
|
|
bitsPerValue,
|
|
referenceValue,
|
|
binaryScaleFactor,
|
|
decimalScaleFactor,
|
|
optimizeScaleFactor,
|
|
|
|
GRIBEXShBugPresent,
|
|
ieeeFloats,
|
|
|
|
laplacianOperatorIsSet,
|
|
laplacianOperator,
|
|
subSetJ,
|
|
subSetK,
|
|
subSetM,
|
|
|
|
pentagonalResolutionParameterJ,
|
|
pentagonalResolutionParameterK,
|
|
pentagonalResolutionParameterM
|
|
) : read_only;
|
|
|
|
meta packingError simple_packing_error(bitsPerValue,binaryScaleFactor,decimalScaleFactor,referenceValue,ibm) : no_copy;
|
|
meta unpackedError simple_packing_error(zero,binaryScaleFactor,decimalScaleFactor,referenceValue,ibm) : no_copy;
|
|
|
|
# nearest sh(values,radius,J,K,M);
|
|
|
|
meta numberOfCodedValues g1number_of_coded_values_sh_complex(bitsPerValue,offsetBeforeData,offsetAfterData,halfByte,numberOfValues,subSetJ,subSetK,subSetM) : dump;
|
|
|
|
template statistics "common/statistics_spectral.def";
|