Page History

# 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.

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

from __future__ import absolute_import

import sys
import traceback

from eccodes import *

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


def example():
    # open bufr file
    f = open(INPUT, 'rb')

    cnt = 0

    # loop for the messages in the file
    while 1:
        }
# get handle for message
    
    bufr    free(values);
= codes_bufr_new_from_file(f)
        if bufr is None:
         /* delete handle */break

        codes_handle_delete(h);print("message: %s" % cnt)

        
# we need to instruct ecCodes to expand cnt++;
all the   }descriptors
    
    fclose(in);
    return 0;
}

/*
 * 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: bufr_expanded
 *
 * Description: how to read all the expanded data values from BUFR messages.
 *
 */

#include "eccodes.h"

static void usage(const char* prog) {
    printf("usage: %s infile\n",prog);
    exit(1);
}

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

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

    double *values = NULL;
    size_t values_len=0;
    int i, err=0;
    int cnt=0;

    if (argc!=2) usage(argv[0]);

    filename=argv[1];

    /* open bufr file */
    in=fopen(filename,"rb");
    if (!in) {
        printf("ERROR: unable to open file %s\n", filename);
        return 1;
    }

    /* loop over

# 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_1.bufr'
VERBOSE=1 # verbose error reporting
  
def example():
    
    # open bufr file
    f = open(INPUT)

    cnt=0    
    
    # loop for the messages in the bufr file */
    while 1:
        # get handle for message((h = codes_handle_new_from_file(NULL, in, PRODUCT_BUFR, &err)) != NULL || err != CODES_SUCCESS)
    {
        if gid(h = codes_bufr_new_from_file(f)= NULL) {
        if gid is None printf("Error: break

unable to create handle     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 all the expanded data values
        #-----------------------------------for message %d\n",cnt);
        key='numericValues'
    cnt++;
    
        # get sizecontinue;
        num=codes_get_size(gid,key)}

        print/* get 'the  size of %sthe is: %s' % (key,num)values array*/
        
  CODES_CHECK(codes_get_size(h, "numericValues", &values_len),0);
      # get values
printf("  number of expanded values: %lu\n",  values=codes_get_array(gid,key)_len);

        for/* iallocate in xrange(len(values)):
   array for data values */
        values print "   %d %.10e" % (i+1,values[i])= (double*)malloc(values_len*sizeof(double));

        cnt+=1

     /* get the expanded data values*/
   # delete handle
   CODES_CHECK(codes_get_double_array(h, "numericValues",    codes_release(gid)

values, &values_len),0);

     # close the file
    f.close()

def main():
   for(i = 0; i < values_len; i++)
        try:{
         example()
   printf(" except CodesInternalError,err: %.10e\n",values[i]);
        if VERBOSE:}

         free(values);

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

   print >>sys.stderr,err.msg

    cnt++;
    return 1}

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