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

ECC-1795: Geoiterator for polar stereographic: Keys orientationOfTheGridInDegrees and LaDInDegrees decoded as integers

This commit is contained in:
shahramn 2024-04-04 17:00:05 +01:00
parent 157dbfb192
commit 49c4ca3cbf
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/grib_iterator_class_polar_stereographic.cc
View File

@ -123,7 +123,8 @@ static int init(grib_iterator* iter, grib_handle* h, grib_arguments* args)
double *lats, *lons; /* arrays for latitudes and longitudes */ double *lats, *lons; /* arrays for latitudes and longitudes */
double lonFirstInDegrees, latFirstInDegrees, radius; double lonFirstInDegrees, latFirstInDegrees, radius;
double x, y, Dx, Dy; double x, y, Dx, Dy;
long nx, ny, centralLongitudeInDegrees, centralLatitudeInDegrees; long nx, ny;
double centralLongitudeInDegrees, centralLatitudeInDegrees;
long alternativeRowScanning, iScansNegatively, i, j; long alternativeRowScanning, iScansNegatively, i, j;
long jScansPositively, jPointsAreConsecutive, southPoleOnPlane; long jScansPositively, jPointsAreConsecutive, southPoleOnPlane;
double centralLongitude, centralLatitude; /* in radians */ double centralLongitude, centralLatitude; /* in radians */
@ -173,9 +174,9 @@ static int init(grib_iterator* iter, grib_handle* h, grib_arguments* args)
return ret; return ret;
if ((ret = grib_get_long_internal(h, s_southPoleOnPlane, &southPoleOnPlane)) != GRIB_SUCCESS) if ((ret = grib_get_long_internal(h, s_southPoleOnPlane, &southPoleOnPlane)) != GRIB_SUCCESS)
return ret; return ret;
if ((ret = grib_get_long_internal(h, s_centralLongitude, &centralLongitudeInDegrees)) != GRIB_SUCCESS) if ((ret = grib_get_double_internal(h, s_centralLongitude, &centralLongitudeInDegrees)) != GRIB_SUCCESS)
return ret; return ret;
if ((ret = grib_get_long_internal(h, s_centralLatitude, &centralLatitudeInDegrees)) != GRIB_SUCCESS) if ((ret = grib_get_double_internal(h, s_centralLatitude, &centralLatitudeInDegrees)) != GRIB_SUCCESS)
return ret; return ret;
if ((ret = grib_get_double_internal(h, s_Dx, &Dx)) != GRIB_SUCCESS) if ((ret = grib_get_double_internal(h, s_Dx, &Dx)) != GRIB_SUCCESS)
return ret; return ret;