Versions Compared

Key

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

...

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
 */
#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[])
{
    FILE *in = NULL;
    FILE *out = NULL;

    /* message handle. Required in all the eccodes calls acting on a message.*/:: err,i,iret
  integer  codes_handle *source_handle = NULL;

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

    if (argc != 3) {
   :: nx, ny
  integer  usage(argv[0]);
        return 1;
    }

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

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

    /* 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:: couldibufr_out
 not create handle for message");
 character(len=2)           fclose(out);
        fclose(in);
        return 1;
   :: }step

  double  /* 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;
        }

        /*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! E.g. we change the bufrHeaderCentre */

    call    CODES_CHECK(codes_set_long(cloneibufr_handleout, "'bufrHeaderCentre"', 222),0);
 
    ! write cloned  /* get the coded message inmessages to a buffer */file
        CODES_CHECK(call codes_get_messagewrite(cloneibufr_handle,&buffer,&size),0);

out,outfile)
    
    /*! writerelase the buffer to a file */
        if(fwrite(buffer,1,size,out) != size) {
            perror("ERROR: could not write message to file");
            return 1;
        }

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

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

    fclose(out);
    fclose(in);

    return 0;
}
Tabs Page
titleFortran 90
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
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
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          clone_id = codes_clone(gid)

        # this is the place where you may wish to modify the clone
       :: err,i,iret codes_set(clone_id, 'bufrHeaderCentre', centre)

  integer      # write the cloned message to a file
        codes_write(clone_id, fout)

        # relase the clone's handle
     :: nx, ny  codes_release(clone_id)

  integer  # release the source's handle
    codes_release(gid)

    fin.close()
    fout.close()


def main():
    try:
        example()
    except CodesInternalError as  :: infile,outfileerr:
  integer      if VERBOSE:
            traceback.print_exc(file=sys.stderr)
        else:
            :: ibufr_inprint >>sys.stderr, err.msg

  integer      return 1

if __name__ == "__main__":
    sys.exit(main())
Tabs Page
titleC
Code Block
languagecpp
titlebufr_clone.c
linenumbersfalse
/*
 * Copyright 2005-2015 ECMWF.
 *
 * This software is licensed under the terms of the Apache Licence Version 2.0
 * which can be obtained at :: 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
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):
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"

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

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

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

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

        #/* relasewrite the clone's handlebuffer to a file */
        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:
    /* relase 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;
}