00001 ! Copyright 2005-2007 ECMWF 00002 ! 00003 ! Licensed under the GNU Lesser General Public License which 00004 ! incorporates the terms and conditions of version 3 of the GNU 00005 ! General Public License. 00006 ! See LICENSE and gpl-3.0.txt for details. 00007 ! 00008 ! 00009 ! Description: how to set a bitmap in a grib message 00010 ! 00011 ! 00012 ! Author: Enrico Fucile 00013 ! 00014 ! 00015 program set_bitmap 00016 use grib_api 00017 implicit none 00018 integer :: infile,outfile 00019 integer :: igrib, iret 00020 integer :: numberOfValues 00021 real, dimension(:), allocatable :: values 00022 real :: missingValue 00023 logical :: grib1Example 00024 00025 grib1Example=.true. 00026 00027 if (grib1Example) then 00028 ! GRIB 1 example 00029 call grib_open_file(infile,'../../data/regular_latlon_surface.grib1','r') 00030 else 00031 ! GRIB 2 example 00032 call grib_open_file(infile,'../../data/regular_latlon_surface.grib2','r') 00033 end if 00034 00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,00035 >00041 call grib_set(igrib_out,