Testing: Increase coverage: F90 grib_get_data

examples/F90/CMakeLists.txt

@@ -25,6 +25,7 @@ if( HAVE_BUILD_TOOLS )
           bufr_copy_message
           grib_get_keys
           grib_get_data
+          grib_get_data_real4
           grib_get_pl
           grib_get_pv
           grib_keys_iterator

examples/F90/grib_get_data_real4.f90

@@ -0,0 +1,76 @@
+! (C) 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.
+!
+!
+!
+program get_data_real4
+   use eccodes
+   implicit none
+   integer            :: ifile
+   integer            :: iret, i
+   real(kind=4), dimension(:), allocatable     :: lats, lons, values
+   integer, dimension(:), allocatable          :: bitmap
+   integer(4)        :: numberOfPoints
+   logical            :: is_missing_value
+   integer            :: count1 = 0, count2 = 0, bitmapPresent = 0, bmp_len = 0
+   integer            :: igrib  ! message identifier
+
+   ifile = 5
+
+   call codes_open_file(ifile, '../../data/reduced_latlon_surface.grib2', 'r')
+
+   ! loop on all the messages in the file
+   do while (.true.)
+      call codes_grib_new_from_file(ifile, igrib, iret)
+      if (iret == CODES_END_OF_FILE) exit
+
+      count1 = count1 + 1
+      print *, "===== Message #", count1
+      call codes_get(igrib, 'numberOfPoints', numberOfPoints)
+      call codes_get(igrib, 'bitmapPresent', bitmapPresent)
+
+      allocate (lats(numberOfPoints))
+      allocate (lons(numberOfPoints))
+      allocate (values(numberOfPoints))
+      if (bitmapPresent == 1) then
+         ! get the bitmap
+         call codes_get_size(igrib, 'bitmap', bmp_len)
+         allocate (bitmap(bmp_len))
+         call codes_get(igrib, 'bitmap', bitmap)
+      end if
+
+      call codes_grib_get_data(igrib, lats, lons, values)
+
+      do i = 1, numberOfPoints
+         ! consult bitmap to see if the i'th value is missing
+         is_missing_value = .false.
+         if (bitmapPresent == 1 .and. bitmap(i) == 0) then
+            is_missing_value = .true.
+         end if
+         ! only print non-missing values
+         if (.not. is_missing_value) then
+            print *, lats(i), lons(i), values(i)
+            count2 = count2 + 1
+         end if
+      end do
+      print *, 'count of non-missing values=', count2
+      if (count2 /= 214661) then
+         call codes_check(-2, 'incorrect number of missing', '')
+      end if
+
+      deallocate (lats)
+      deallocate (lons)
+      deallocate (values)
+
+      call codes_release(igrib)
+
+   end do
+
+   call codes_close_file(ifile)
+
+end program

examples/F90/grib_get_data_real4.sh

@@ -0,0 +1,12 @@
+#!/bin/sh
+# (C) 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.
+
+. ./include.ctest.sh
+
+${examples_dir}/eccodes_f_grib_get_data_real4 > /dev/null