Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

Description

This example shows: how to create a new BUFR

...

message by cloning an existing message.

Source code

Tabs Container
directionhorizontal


Tabs Page
titleCFortran 90


Code Block
languagecppnone
titlebufr_clone.cf90
linenumbersfalse
/*
 *! 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 FORTRAN 90 implementation: bufr_clone
 *!
! * Description: how to create a new BUFR messagesmessage by cloning
! *              an existing message.
!
!
program *bufr_clone
 */
 use eccodes
  implicit none
  integer       :: i,iret
  integer       :: infile,outfile
  integer       :: ibufr_in
  integer       :: ibufr_out

  ! Open source file
  call codes_open_file(infile,'../../data/bufr/syno_multi.bufr','r')

  ! Open target file
  call codes_open_file(outfile,'bufr_clone_test_f.clone.bufr','w')

  ! The first bufr message is loaded from file,
  ! ibufr is the bufr id to be used in subsequent calls
  call codes_bufr_new_from_file(infile,ibufr_in,iret)

  ! Create several clones of this message and alter them
  ! in different ways
  do i=1,3

    ! Clone the current handle
    call codes_clone(ibufr_in, ibufr_out)

    ! This is the place where you may wish to modify the clone
    ! E.g. we change the bufrHeaderCentre
    call codes_set(ibufr_out,'bufrHeaderCentre',222)

    ! Write cloned messages to a file
    call codes_write(ibufr_out,outfile)

    ! Release the clone's handle
    call codes_release(ibufr_out)
  end do

  ! Release the original handle
  call codes_release(ibufr_in)

  call codes_close_file(infile)
  call codes_close_file(outfile)

end program bufr_clone



Tabs Page
titlePython


Code Block
languagepython
titlebufr_clone.py
linenumbersfalse
# 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_clone
#
# Description: how to create a new BUFR message by cloning
# an existing message.
#

from __future__ import absolute_import
import sys
import traceback

from eccodes import *

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


def example():
    # open BUFR file
    fin = open(INPUT, 'rb')

    # open output BUFR file
    fout = open(OUTPUT, 'wb')

    # get handle for message
    bufr = codes_bufr_new_from_file(fin)

    # create several clones of this message and alter them
    # in different ways
    for centre in range(0, 3):
        # clone the message
        clone_id = codes_clone(bufr)

        # this is the place where you may wish to modify the clone
        codes_set(clone_id, 'bufrHeaderCentre', centre)
#include <stdio.h>
#include "eccodes.h"

void usage(char *app)
{
    fprintf(stderr,"Usage is: %s input_file ouput_file\n", app);
}

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

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

    const void *buffer = NULL;
    size_t size = 0;
    int i, err = 0;

    if (argc != 3) {
        usage(argv[0]);
        return 1;
    }

    in = fopen(argv[1],"r");
    out = fopen(argv[2],"w");

    /* open input and output */
    if (!in || !out) {
        perror("ERROR: unable to open files");
        if (out) fclose(out);
        if (in) fclose(in);
        return 1;
    }

    /* create a handle for the first message */
    source_handle = codes_handle_new_from_file(NULL,in,PRODUCT_BUFR,&err);

    if (source_handle == NULL) {
        perror("ERROR: could not create handle for message");
        fclose(out);
        fclose(in);
        return 1;
    }

    /* create several clones of this message and alter them
       in different ways */

    for(i=0; i < 3; i++) {

        /* clone the current handle */
        codes_handle *clone_handle = codes_handle_clone(source_handle);

        if (clone_handle == NULL) {
            perror("ERROR: could not clone field");
            return 1;
        }

        /* This is the place where you may wish to modify the clone
           E.g. we change the bufrHeaderCentre */

        CODES_CHECK(codes_set_long(clone_handle, "bufrHeaderCentre", 222),0);

        /* get the coded message in a buffer */
        CODES_CHECK(codes_get_message(clone_handle,&buffer,&size),0);

        /* write the buffer to a file */
        if(fwrite(buffer,1,size,out) != size) {
        # write   perror("ERROR: could not writethe cloned message to a file");
        codes_write(clone_id, fout)

   return 1;
    # release the clone's }handle

        /* relase the clone's handle */
        codes_handle_deletecodes_release(clone_handle);
    }id)

    /*# release the source's handle */
    codes_handle_delete(source_handle);release(bufr)

    fclosefin.close(out);
    fclosefout.close(in);


def main():
   return 0;
}
Tabs Page
titleFortran 90
Code Block
languagenone
titlebufr_clone.f90
linenumbersfalse
! 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.
!
!  FORTRAN 90 implementation: bufr_clone
!
!  Description: how to create a new BUFR message by cloning 
!               an existing message.
!
!
!
program bufr_clone
  use eccodes
  implicit none
  integer                                       :: err,i,iret
  integer                                       :: nx, ny
  integer                                       :: infile,outfile
  integer                                       :: ibufr_in
  integer                                       :: ibufr_out
  character(len=2)                              :: step
  double precision, dimension(:,:), allocatable :: field2D
  
  !open source file 
  call codes_open_file(infile,'../../data/bufr/syno_multi.bufr','r')
  
  ! open target file
  call codes_open_file(outfile,'bufr_clone_test_f.clone.bufr','w')
  
  ! the first bufr message is loaded from file
  ! ibufr is the bufr id to be used in subsequent calls
  call codes_bufr_new_from_file(infile,ibufr_in,iret)

  ! create several clones of this message and alter them
  ! in different ways

  do i=1,3
    
    ! clone the current handle
    call codes_clone(ibufr_in, ibufr_out)

    ! This is the place where you may wish to modify the clone 
    ! E.g. we change the bufrHeaderCentre 
    call codes_set(ibufr_out,'bufrHeaderCentre',222)
 
    ! write cloned messages to a file
    call codes_write(ibufr_out,outfile)
    
    ! relase the clone's handle
    call codes_release(ibufr_out)
  end do  
  
! relase the sources 's handle
  call codes_release(ibufr_in)
  
  call codes_close_file(infile)
  call codes_close_file(outfile)

end program bufr_clone
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
titleC


Code Block
languagecpp
titlebufr_clone.c
linenumbersfalse
/*
 * 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_clone
 *
 * Description: how to create new BUFR messages by cloning
 *              an existing message.
 *
 */
#include <stdio.h>
#include "eccodes.h"

static void usage(const char *app)
{
    fprintf(stderr,"Usage is: %s input_file ouput_file\n", app);
}

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

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

    const void *buffer = NULL;
    size_t size = 0;
    int i, err = 0;

    if (argc != 3) {
        usage(argv[0]);
        return 1;
    }

    in = fopen(argv[1],"rb");
    out = fopen(argv[2],"wb");

    /* open input and output */
    if (!in || !out) {
        perror("ERROR: unable to open files");
        if (out) fclose(out);
        if (in) fclose(in);
        return 1;
    }

    /* create a handle for the first message */
    source_handle = codes_handle_new_from_file(NULL, in, PRODUCT_BUFR, &err);

    if (source_handle == NULL) {
        perror("ERROR: could not create handle for message");
        fclose(out);
        fclose(in);
        return 1;
    }

    /* create several clones of this message and alter them
       in different ways */

    for(i=0; i < 3; i++) {

        /* clone the current handle */
        codes_handle *clone_handle = codes_handle_clone(source_handle);

        if (clone_handle == NULL) {
            perror("ERROR: could not clone field");
Tabs Page
titlePython
Code Block
languagepython
titlebufr_clone.py
linenumbersfalse
# 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_clone
#
# Description: how to create a new BUFR message by cloning
# an existing message.
#
#

import traceback
import sys

from eccodes import *

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

def example():

    # open bufr file
    fin = open(INPUT)

    # open otput bufr file
    fout = open(OUTPUT,'w')

    # get handle for message
    gid = codes_bufr_new_from_file(fin)

    # create several clones of this message and alter them
    # in different ways

    for centre in range(0,3):

        # clone the message
 return 1;
      clone_id = codes_clone(gid) }

        #/* thisThis is the place where you may wish to modify the clone
           E.g. we change the bufrHeaderCentre */

        CODES_CHECK(codes_setset_long(clone_idhandle,'bufrHeaderCentre',centre) "bufrHeaderCentre", 222),0);

        #/* writeget the clonedcoded message toin a filebuffer */
        CODES_CHECK(codes_get_writemessage(clone_id,fout)handle, &buffer, &size),0);

        /* write the #buffer relaseto thea clone's handlefile */
        codes_release(clone_id)

    # release the source's handle
    codes_release(gid)

    fin.close()
    fout.close()

def main():
    try:
if(fwrite(buffer, 1, size, out) != size) {
            perror("ERROR: could not write message to file");
            return example()
1;
       except CodesInternalError,err: }

        if VERBOSE:
    /* release the clone's handle */
        traceback.print_exc(file=sys.stderr)codes_handle_delete(clone_handle);
        else:}

    /* release the source's handle */
   print >>sys.stderr,err.msg codes_handle_delete(source_handle);

    fclose(out);
    return 1fclose(in);

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