Page History

Description

This example shows: how to create a new BUFR

...

message by cloning an existing message.

Source code

Tabs Container

direction	horizontal

Tabs Page

title	CFortran 90

Code Block

language	cppnone
title	bufr_clone.cf90
linenumbers	false

/*
! * 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
 */
#include <stdio.h>
#include "eccodes.h"

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

int main(int argc, char *argv[])
{
     :: i,iret
  integer      FILE *in = NULL;:: infile,outfile
  integer  FILE *out = NULL;
   :: ibufr_in
  integer  /* 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");
        fclose(out);
        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 centre */:: 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

title	Python

Code Block

language	python
title	bufr_clone.py
linenumbers	false

# 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
        CODES_CHECK(clone_id = codes_set_long(clone_handle, "centre", 222),0);clone(bufr)

        /* get# this is the coded message in a buffer */ place where you may wish to modify the clone
        CODES_CHECK(codes_get_messageset(clone_handle,&buffer,&size),0);
       id, 'bufrHeaderCentre', centre)

        /*# write the cloned buffermessage to a file */
        if(fwrite(buffer,1,size,out) != size) {
codes_write(clone_id, fout)

        # release the clone's handle
   perror("ERROR: could not write message to file"); codes_release(clone_id)

    # release the source's handle
    return 1;codes_release(bufr)

    fin.close()
    }
 fout.close()


def main():
    try:
   
     example()
   /* relaseexcept theCodesInternalError clone's handle */as err:
        codes_handle_delete(clone_handle);if VERBOSE:
    }
        traceback.print_exc(file=sys.stderr)
    /* release the source's handle */ else:
    codes_handle_delete(source_handle);

    fclose(out);
    fclose(in);

sys.stderr.write(err.msg + '\n')

        return 0;
}1


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

Tabs Page

title	Fortran 90C

Code Block

language	nonecpp
title	bufr_clone.f90c
linenumbers	false

!/*
 * 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 */

/*
 * C implementation: bufr_clone
! *
! * Description: how to create a new BUFR messagemessages by cloning 
! *              an existing message.
!
!
!
program bufr_clone
  use eccodes
  implicit none
  integer                                       :: err,i,iret
  integer                                       :: nx, ny
  integer       *
 */
#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 :: infile,outfileoutput */
  integer  if (!in || !out) {
        perror("ERROR: unable to open files");
        if (out) fclose(out);
        if   :: ibufr_in(in) fclose(in);
  integer      return 1;
    }

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

    if (source_handle == :: ibufr_outNULL) {
  character(len=2)      perror("ERROR: could not create handle                    :: 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 centre 
    call codes_set(ibufr_out,'centre',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

Tabs Page

title	Python

Code Block

language	python
title	bufr_clone.py
linenumbers	false

# 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):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");
 clone the message
        clone_id = codes_clone(gid) return 1;
        }

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

        CODES_CHECK(codes_set_long(clone_id,'centre',centre)
 handle, "bufrHeaderCentre", 222),0);

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

   # relase the clone's handle
 /* write the buffer to a  codes_release(clone_id)file */
  
    # release the source's handle
    codes_release(gid)
 if(fwrite(buffer, 1, size, out) != size) {
        
    fin.close()
    fout.close()

def main():
perror("ERROR: could not write message to file");
     try:
       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;
}

Space shortcuts

Page tree

Versions Compared

Old Version 2

New Version Current

Key

Description

Source code