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

ECC-445: Fix case when lon1 > lon2

This commit is contained in:
Shahram Najm 2018-06-20 17:41:44 +01:00
parent 0a8da3ca2f
commit def3e2f382
1 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/grib_gaussian_reduced.c
View File

@ -262,7 +262,7 @@ static void gaussian_reduced_row(
//auto Nw = (w / inc).integralPart();
Fraction_value_type Nw = fraction_integralPart( fraction_operator_divide(w, inc) );
Fraction_type Nw_inc = fraction_operator_multiply_n_Frac(Nw, inc);
//if (Nw * inc < w) { //??
//if (Nw * inc < w) {
if (fraction_operator_less_than(Nw_inc, w)) {
Nw += 1;
}
@ -394,11 +394,15 @@ void grib_get_reduced_row(long pl, double lon_first, double lon_last, long* npoi
void grib_get_reduced_row2(long pl, double lon_first, double lon_last, long* npoints, long* ilon_first, long* ilon_last )
{
long long Ni_globe = pl;
Fraction_type west = fraction_construct_from_double(lon_first);
Fraction_type east = fraction_construct_from_double(lon_last);
Fraction_type west;
Fraction_type east;
long long the_count;
double the_lon1, the_lon2;
/*printf("Using gaussian_reduced_row...\n");*/
while (lon_last < lon_first) lon_last += 360;
west = fraction_construct_from_double(lon_first);
east = fraction_construct_from_double(lon_last);
gaussian_reduced_row(
Ni_globe, // plj
west, // lon_first