mirror of https://github.com/ecmwf/eccodes.git
79 lines
2.6 KiB
Fortran
79 lines
2.6 KiB
Fortran
! (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.
|
|
!
|
|
!
|
|
! Description: How to use codes_grib_find_nearest and codes_get_element
|
|
!
|
|
!
|
|
!
|
|
program find
|
|
use eccodes
|
|
implicit none
|
|
integer :: npoints
|
|
integer :: infile
|
|
integer :: igrib, ios, i
|
|
real(8), dimension(:), allocatable :: lats, lons
|
|
real(8), dimension(:), allocatable :: nearest_lats, nearest_lons
|
|
real(8), dimension(:), allocatable :: distances, values, lsm_values
|
|
integer(kind=kindOfInt), dimension(:), allocatable :: indexes
|
|
|
|
! initialization
|
|
open (unit=1, file="../../data/list_points", form="formatted", action="read")
|
|
read (unit=1, fmt=*) npoints
|
|
allocate (lats(npoints))
|
|
allocate (lons(npoints))
|
|
allocate (nearest_lats(npoints))
|
|
allocate (nearest_lons(npoints))
|
|
allocate (distances(npoints))
|
|
allocate (lsm_values(npoints))
|
|
allocate (values(npoints))
|
|
allocate (indexes(npoints))
|
|
do i = 1, npoints
|
|
read (unit=1, fmt=*, iostat=ios) lats(i), lons(i)
|
|
if (ios /= 0) then
|
|
npoints = i - 1
|
|
exit
|
|
end if
|
|
end do
|
|
close (unit=1)
|
|
call codes_open_file(infile, &
|
|
'../../data/reduced_gaussian_lsm.grib1', 'r')
|
|
|
|
! A new grib message is loaded from file
|
|
! igrib is the grib id to be used in subsequent calls
|
|
call codes_grib_new_from_file(infile, igrib)
|
|
|
|
call codes_grib_find_nearest(igrib, .true., lats, lons, nearest_lats, nearest_lons, lsm_values, distances, indexes)
|
|
call codes_release(igrib)
|
|
|
|
call codes_close_file(infile)
|
|
|
|
! will apply it to another GRIB
|
|
call codes_open_file(infile, &
|
|
'../../data/reduced_gaussian_pressure_level.grib1', 'r')
|
|
call codes_grib_new_from_file(infile, igrib)
|
|
|
|
call codes_get_element(igrib, "values", indexes, values)
|
|
call codes_release(igrib)
|
|
call codes_close_file(infile)
|
|
|
|
do i = 1, npoints
|
|
print *, lats(i), lons(i), nearest_lats(i), nearest_lons(i), distances(i), lsm_values(i), values(i)
|
|
end do
|
|
|
|
deallocate (lats)
|
|
deallocate (lons)
|
|
deallocate (nearest_lats)
|
|
deallocate (nearest_lons)
|
|
deallocate (distances)
|
|
deallocate (lsm_values)
|
|
deallocate (values)
|
|
deallocate (indexes)
|
|
|
|
end program find
|