Page History

# Copyright 2005-2015 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.

#
# Python implementation: bufr_get_keys
#
# Description: how to read values of different type of keys from BUFR messages.
#
#

import traceback
import sys

from eccodes import *

INPUT = '../../data/bufr/syno_multi.bufr'
VERBOSE = 1  # verbose error reporting


def example():

    # open bufr file
    f = open(INPUT)

    cnt = 0

    # loop for the messages in the file
    while 1:
        # get handle for message
        gid = codes_bufr_new_from_file(f)
        printf("  numericValues:\n");if gid is None:
        for(i = 0; i < values_len; i++) break

        {
print "message: %s" % cnt

        printf("# we  %.10e\n",values[i]);
        }

need to instruct ecCodes to expand all the descriptors
        /* free allocated arrays */# i.e. unpack the data values
        free(descriptors);codes_set(gid, 'unpack', 1)

        free(values);# ---------------------------------------------
        free(typicalDate);

        /* delete handle */
# get values for keys holding a single value
        #   codes_handle_delete(h);---------------------------------------------
        # Native type integer
        key = 'blockNumber'

        cnt++;try:
    }

       fclose(in);
 print '  return 0;
}

/*
 * Copyright 2005-2015 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: bufr_get_keys
 *
 * Description: how to read values of different type of keys from BUFR messages.
 *
 */

#include "eccodes.h"

int main(int argc,char* argv[])
{
    FILE* in = NULL;

    /* message handle. Required in all the eccodes calls acting on a message.*/
    codes_handle* h=NULL;

    long *descriptors=NULL;
    double *values = NULL;
    char* typicalDate= NULL; 
    long longVal;
    double doubleVal;
    size_t values_len=0, desc_len=0, len=0;
    int i, err=0;
    int cnt=0;
    const char* infile = "../../data/bufr/syno_multi.bufr";

    in=fopen(infile,"r");
    if (!in) {
        printf("ERROR: unable to open file %s\n", infile);
        return 1;
    }

    /* loop over the messages in the bufr file */
    while ((h = codes_handle_new_from_file(NULL,in,PRODUCT_BUFR,&err)) != NULL || err != CODES_SUCCESS)
    {
        if (h == NULL) {

# Copyright 2005-2015 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.

#
# Python implementation: bufr_get_keys
#
# Description: how to read values of different type of keys from BUFR messages.
#
#

import traceback
import sys

from eccodes import *

INPUT='../../data/bufr/syno_multi.bufr'
VERBOSE=1 # verbose error reporting

def example():

    # open bufr file
    f = open(INPUT)

    cnt=0

    # loop for the messages in the file
    while 1:
        # get handle for message
        gid = codes_bufr_new_from_file(f)
        if gid is None: break

        print "message: %s" % cnt

        # we need to instruct ecCodes to expand all the descriptors
        # i.e. unpack the data values
        codes_set(gid,'unpack',1)

        #----------------------------------------------
        # get values for keys holding a single value
        #----------------------------------------------
        #Native type integer
        key='blockNumber'
        
        try:
            print '  %s: %s' % (key,codes_get(gid,key))
   printf("Error: unable to create handle for message %d\n",cnt);
     except CodesInternalError,err:
      cnt++;
      print 'Error with key="%s" : %s' % (key,err.msg)continue;
        }

        #Native type integerprintf("message: %d\n",cnt);

        key='stationNumber'
        try:/* we need to instruct ecCodes to expand the descriptors 
          i.e. unpack the data values */
  print '  %s: %s' % CODES_CHECK(key,codes_get(gid,key))
        except CodesInternalError,err:set_long(h,"unpack",1),0);

        /* read  and print 'Errorsome with key="%s" : %s' % (key,err.msg) 

data values */ 

        /* #Nativelong typevalue float*/
        key='airTemperatureAt2M'CODES_CHECK(codes_get_long(h,"blockNumber",&longVal),0);
        try:
printf("  blockNumber: %ld\n",longVal);

        /* long value */
   print '  %s: %s' % (key, CODES_CHECK(codes_get(gid,key))_long(h,"stationNumber",&longVal),0);
        except CodesInternalError,err:printf("  stationNumber: %ld\n",longVal);

        /* double value */
 print 'Error with key="%s" : %s' % (key,err.msg) 

 CODES_CHECK(codes_get_double(h,"airTemperatureAt2M",&doubleVal),0);
        printf(" #Native type stringairTemperatureAt2M: %f\n",doubleVal);

        key='typicalDate'
    /* ---- string value  -----------------*/

    try:
    /* get the size and allocate memory*/
  print '  %s: %s' % CODES_CHECK(key,codes_get_length(gidh,key))
        except CodesInternalError,err: "typicalDate", &len), 0);
            print 'Error with key="%s" : %s' % (key,err.msg) typicalDate = (char*)malloc(len*sizeof(char));

        #---------------------------------/* get the values*/
        # get values for an arraycodes_get_string(h, "typicalDate", typicalDate, &len);
        #---------------------------------
        # Native type integer
        key='bufrdcExpandedDescriptors'printf("  typicalDate: %s\n", typicalDate);

        /* ---- array of long ----------------*/

        #/* get the size and allocate memory*/
        num=CODES_CHECK(codes_get_size(gid,key)
    h,"bufrdcExpandedDescriptors",&desc_len),0);
    print  '  size of %s is: %s'descriptors %= (key,num)long*)malloc(desc_len*sizeof(long));

        #/* get the values */
        values=CODES_CHECK(codes_get_array(gid,key)long_array(h,"bufrdcExpandedDescriptors",descriptors,&desc_len),0);
        printf("  bufrdcExpandedDescriptors:\n");
        for(i = 0; i in< xrange(len(values)):desc_len; i++)
        {
    print    "   %d %06dprintf(" % (i+1,values  %ld\n",descriptors[i])
;
        # Native type float
        key='numericValues'}

        /* ---- array of double ---------------*/

        #/* get the size and allocate memory*/
        num=CODES_CHECK(codes_get_size(gid,key)
    h,"numericValues",&values_len),0);
    print  '  size of %s is: %s'values %= (key,num)double*)malloc(values_len*sizeof(double));

        #/* get the values*/
        values=CODES_CHECK(codes_get_double_array(gid,key)h,"numericValues",values,&values_len),0);
        forprintf(" i in xrange(len(values)): numericValues:\n");
        for(i = 0; i print "   %d %.10e" % (i+1,values[i])

< values_len; i++)
        cnt+=1

{
         # delete handle
 printf("   %.10e\n",values[i]);
    codes_release(gid)

    #}

 close the file
    f.close()

def main():
    try: /* free allocated arrays */
        examplefree(descriptors);
      except CodesInternalError,err:  free(values);
        if VERBOSE:free(typicalDate);

        /* delete   traceback.print_exc(file=sys.stderr)handle */
        else:codes_handle_delete(h);

        cnt++;
    print >>sys.stderr,err.msg}

    fclose(in);
    return 1

if __name__ == "__main__":
    sys.exit(main())0;
}