Description
This example shows: How to use keys_iterator to get all the available keys in a GRIB message.
Source code
Tabs Container |
---|
|
Tabs Page |
---|
|
Code Block |
---|
language | none |
---|
title | grib_keys_iterator.f90 |
---|
linenumbers | false |
---|
| ! (C) 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.
!
!
! 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
|
|
Tabs Page |
---|
|
Code Block |
---|
language | python |
---|
title | grib_keys_iterator.py |
---|
linenumbers | false |
---|
| #
# (C) 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 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(iteridgid, 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())
|
|
Tabs Page |
---|
|
Code Block |
---|
language | cpp |
---|
title | grib_keys_iterator.c |
---|
linenumbers | false |
---|
| /*
* (C) 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.
*/
/*
* 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);
}
|
|
|
...