/usr/share/doc/libgrib-api-dev/examples/count_messages.py is in libgrib-api-dev 1.25.0-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 | #
# Copyright 2005-2017 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_pf_ecmwf.grib2'
VERBOSE = 1 # verbose error reporting
def example():
f = open(INPUT)
mcount = grib_count_in_file(f)
gid_list = [grib_new_from_file(f) for i in range(mcount)]
f.close()
keys = [
'Ni',
'Nj',
'latitudeOfFirstGridPointInDegrees',
'longitudeOfFirstGridPointInDegrees',
'latitudeOfLastGridPointInDegrees',
'longitudeOfLastGridPointInDegrees',
'jDirectionIncrementInDegrees',
'iDirectionIncrementInDegrees',
]
for i in range(mcount):
gid = gid_list[i]
print ("processing message number", i + 1)
for key in keys:
print ('%s=%g' % (key, grib_get(gid, key)))
print ('There are %d, average is %g, min is %g, max is %g' % (
grib_get_size(gid, 'values'),
grib_get(gid, 'average'),
grib_get(gid, 'min'),
grib_get(gid, 'max')
))
print ('-' * 100)
grib_release(gid)
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())
|