Description
This example shows: How to use keys_iterator to get all the available keys in a GRIB message.
Source code
grib_keys_iterator.f90
! (C) Copyright 2005- 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. ! ! ! Description: ! How to use keys_iterator to get all the available ! keys in a GRIB message. ! ! ! program keys_iterator use eccodes implicit none character(len=20) :: name_space integer :: kiter,ifile,igrib,iret character(len=256) :: key character(len=256) :: value character(len=512) :: all1 integer :: grib_count call codes_open_file(ifile, & '../../data/regular_latlon_surface.grib1','r') ! Loop on all the messages in a file. call codes_grib_new_from_file(ifile,igrib, iret) grib_count=0 do while (iret /= CODES_END_OF_FILE) grib_count=grib_count+1 write(*,*) '-- GRIB N. ',grib_count,' --' ! valid name_spaces are ls and mars name_space='ls' call codes_keys_iterator_new(igrib,kiter,name_space) do call codes_keys_iterator_next(kiter, iret) if (iret .ne. CODES_SUCCESS) exit !terminate the loop call codes_keys_iterator_get_name(kiter,key) call codes_get(igrib,trim(key),value) all1=trim(key)// ' = ' // trim(value) write(*,*) trim(all1) end do call codes_keys_iterator_delete(kiter) call codes_release(igrib) call codes_grib_new_from_file(ifile,igrib, iret) end do call codes_close_file(ifile) end program keys_iterator
grib_keys_iterator.py
# # (C) Copyright 2005- 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. # import traceback import sys from eccodes import * INPUT = '../../data/reduced_latlon_surface.grib1' VERBOSE = 1 # verbose error reporting def example(): f = open(INPUT, 'rb') while 1: gid = codes_grib_new_from_file(f) if gid is None: break iterid = codes_keys_iterator_new(gid, 'ls') # Different types of keys can be skipped # codes_skip_computed(iterid) # codes_skip_coded(iterid) # codes_skip_edition_specific(iterid) # codes_skip_duplicates(iterid) # codes_skip_read_only(iterid) # codes_skip_function(iterid) while codes_keys_iterator_next(iterid): keyname = codes_keys_iterator_get_name(iterid) keyval = codes_get_string(gid, keyname) print("%s = %s" % (keyname, keyval)) codes_keys_iterator_delete(iterid) codes_release(gid) f.close() def main(): try: example() except CodesInternalError as err: if VERBOSE: traceback.print_exc(file=sys.stderr) else: sys.stderr.write(err.msg + '\n') return 1 if __name__ == "__main__": sys.exit(main())
grib_keys_iterator.c
/* * (C) Copyright 2005- 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. */ /* * C Implementation: grib_keys_iterator * * Description: * Example on how to use keys_iterator functions and the * codes_keys_iterator structure to get all the available * keys in a GRIB message. * */ #include <assert.h> #include <stdlib.h> #include <stdio.h> #include <unistd.h> #include <string.h> #include "eccodes.h" #define MAX_VAL_LEN 1024 static void usage(char* progname); int main(int argc, char *argv[]) { /* To skip read only and computed keys unsigned long key_iterator_filter_flags=CODES_KEYS_ITERATOR_SKIP_READ_ONLY | CODES_KEYS_ITERATOR_SKIP_COMPUTED; */ unsigned long key_iterator_filter_flags = CODES_KEYS_ITERATOR_ALL_KEYS | CODES_KEYS_ITERATOR_SKIP_DUPLICATES; /* Choose a namespace. E.g. "ls", "time", "parameter", "geography", "statistics" */ const char* name_space="ls"; /* name_space=NULL to get all the keys */ /* char* name_space=0; */ FILE* f = NULL; codes_handle* h=NULL; int err=0; int msg_count=0; char value[MAX_VAL_LEN]; size_t vlen=MAX_VAL_LEN; if (argc != 2) usage(argv[0]); f = fopen(argv[1],"r"); if(!f) { perror(argv[1]); exit(1); } while((h = codes_handle_new_from_file(0,f,PRODUCT_GRIB,&err)) != NULL) { codes_keys_iterator* kiter=NULL; msg_count++; printf("-- GRIB N. %d --\n",msg_count); if(!h) { printf("ERROR: Unable to create grib handle\n"); exit(1); } kiter=codes_keys_iterator_new(h,key_iterator_filter_flags,name_space); if (!kiter) { printf("ERROR: Unable to create keys iterator\n"); exit(1); } while(codes_keys_iterator_next(kiter)) { const char* name = codes_keys_iterator_get_name(kiter); vlen=MAX_VAL_LEN; bzero(value,vlen); CODES_CHECK(codes_get_string(h,name,value,&vlen),name); printf("%s = %s\n",name,value); /* Alternative way of getting the string value */ CODES_CHECK(codes_keys_iterator_get_string(kiter, value, &vlen),0); } codes_keys_iterator_delete(kiter); codes_handle_delete(h); } fclose(f); return 0; } static void usage(char* progname) { printf("\nUsage: %s grib_file\n",progname); exit(1); }