core
/
eccodes
mirror of https://github.com/ecmwf/eccodes.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

GRIB-814: grib_nearest on Lambert Conformal grid picks wrong point

This commit is contained in:
Shahram Najm 2015-08-11 14:42:07 +01:00
parent 03bdf4f77c
commit f3bfa78c22
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
src/grib_nearest_class_lambert_conformal.c
View File

@ -186,7 +186,8 @@ static int find(grib_nearest* nearest, grib_handle* h,
int ilat=0, ilon=0;
int idx_upper=0, idx_lower=0;
double lat1=0, lat2=0; /* inlat will be between these */
double dist=0, lat_delta=0.1;
double dist=0;
const double LAT_DELTA = 10.0; /* in degrees */
if (grib_is_missing(h,self->Ni,&ret)) {
grib_context_log(h->context, GRIB_LOG_DEBUG,"Key '%s' is missing", self->Ni);
@ -234,7 +235,7 @@ static int find(grib_nearest* nearest, grib_handle* h,
i = 0;
while(grib_iterator_next(iter,&lat,&lon,&the_value))
{
if (lat > lat2+lat_delta || lat < lat1-lat_delta) {
if (lat > lat2+LAT_DELTA || lat < lat1-LAT_DELTA) {
/* Ignore latitudes too far from our point */
}
else {