Description
This example shows: How to create a new BUFR message by cloning an existing message.Source code
bufr_clone.c
/* * 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: clone * * Description: How to create a new BUFR message 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; codes_handle *source_handle = NULL; const void *buffer = NULL; size_t size = 0; int err = 0; if (argc != 3) { usage(argv[0]); return 1; } in = fopen(argv[1],"r"); out = fopen(argv[2],"w"); if (!in || !out) { perror("ERROR: unable to open files"); fclose(out); fclose(in); return 1; } /* loop over the messages in the source bufr and clone them */ while ((source_handle = bufr_new_from_file(0,in,&err))!=NULL) { /* 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 */ CODES_CHECK(codes_set_long(clone_handle, "centre", 250),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) { perror(argv[1]); 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; }
bufr_clone.f90
! 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. ! ! ! ! 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) do while (iret/=CODES_END_OF_FILE) ! 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',250) ! write cloned messages to a file call codes_write(ibufr_out,outfile) ! relase the clone's handle call codes_release(ibufr_out) ! relase the sources 's handle call codes_release(ibufr_in) ! next message from source call codes_bufr_new_from_file(infile,ibufr_in,iret) end do call codes_close_file(infile) call codes_close_file(outfile) end program bufr_clone