Commit 35785c6d authored by Uwe Schulzweida's avatar Uwe Schulzweida
Browse files

grib reduced_gg: gribapi gridType detector doesn't like lonIncr

parent 14ad5c33
......@@ -3,6 +3,10 @@
* using CGRIBEX library version 1.9.1
* Version 1.9.5 released
2018-07-09 Uwe Schulzweida
* grib reduced_gg: gribapi gridType detector doesn't like lonIncr
2018-07-06 Uwe Schulzweida
* set_gridtype: ignore attribute CDI_grid_type for curvilinear and unstructured grids
......
......@@ -1653,7 +1653,7 @@ void cgribexDefGridRegular(int *isec2, double *fsec2, int gridID, int gridtype,
{
ISEC2_FirstLon = 0;
ISEC2_LastLon = (int)lround(1000*(360.-360./(nlat*2)));
ISEC2_LonIncr = (int)lround(1000*360./(nlat*2));
//ISEC2_LonIncr = (int)lround(1000*360./(nlat*2)); // gribapi gridType detector doesn't like lonIncr
}
else
{
......
......@@ -1993,7 +1993,7 @@ void gribapiDefGridRegular(grib_handle *gh, int gridID, int gridtype, bool gridI
if ( gridtype != GRID_GAUSSIAN_REDUCED ) GRIB_CHECK(my_grib_set_long(gh, "Ni", nlon), 0);
GRIB_CHECK(my_grib_set_double(gh, "longitudeOfFirstGridPointInDegrees", xfirsto), 0);
GRIB_CHECK(my_grib_set_double(gh, "longitudeOfLastGridPointInDegrees", xlasto), 0);
GRIB_CHECK(my_grib_set_double(gh, "iDirectionIncrementInDegrees", xinc), 0);
if ( gridtype != GRID_GAUSSIAN_REDUCED ) GRIB_CHECK(my_grib_set_double(gh, "iDirectionIncrementInDegrees", xinc), 0);
GRIB_CHECK(my_grib_set_long(gh, "Nj", (long)nlat), 0);
GRIB_CHECK(my_grib_set_double(gh, "latitudeOfFirstGridPointInDegrees", yfirst), 0);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment