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
|
#
# Copyright 2005-2018 ECMWF.
#
# This software is licensed under the terms of the Apache Licence Version 2.0
# which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
#
# In applying this licence, ECMWF does not waive the privileges and immunities granted to it by
# virtue of its status as an intergovernmental organisation nor does it submit to any jurisdiction.
#
from __future__ import print_function
import traceback
import sys
from gribapi import *
INPUT = '../../data/tigge/tigge_ecmf_pl_t.grib'
OUTPUT = 'out.set_missing.grib'
VERBOSE = 1 # verbose error reporting
def example():
fin = open(INPUT)
fout = open(OUTPUT, 'w')
gid = grib_new_from_file(fin)
grib_set_long(gid, "scaledValueOfFirstFixedSurface", 15)
grib_set_long(gid, "scaleFactorOfFirstFixedSurface", 1)
level = grib_get_double(gid, "level")
assert(level == 1.5)
# set type of level to surface
grib_set(gid, 'typeOfFirstFixedSurface', 'sfc')
grib_set_missing(gid, 'scaleFactorOfFirstFixedSurface')
grib_set_missing(gid, 'scaledValueOfFirstFixedSurface')
grib_write(gid, fout)
grib_release(gid)
fin.close()
fout.close()
def main():
try:
example()
except GribInternalError as err:
if VERBOSE:
traceback.print_exc(file=sys.stderr)
else:
print (err.msg, file=sys.stderr)
return 1
if __name__ == "__main__":
sys.exit(main())
|