<div dir="ltr">Hi,<div><br></div><div>If I use GDAL to read GRIB files like these:</div><div><a href="ftp://hydro1.sci.gsfc.nasa.gov/data/s4pa/NLDAS/NLDAS_NOAH0125_H.002/2015/065/NLDAS_NOAH0125_H.A20150306.1300.002.grb">ftp://hydro1.sci.gsfc.nasa.gov/data/s4pa/NLDAS/NLDAS_NOAH0125_H.002/2015/065/NLDAS_NOAH0125_H.A20150306.1300.002.grb</a><br></div><div><br></div><div>the CRS I get back is:</div><div><div>GEOGCS["Coordinate System imported from GRIB file",</div><div>    DATUM["unknown",</div><div>        SPHEROID["Sphere",6371200,0]],</div><div>    PRIMEM["Greenwich",0],</div><div>    UNIT["degree",0.0174532925199433]]</div></div><div><br></div><div>Plotting these files shows that they are shifted 1.5 pixels north, but setting CRS to EPSG:4326 that uses SPHEROID["WGS 84",6378137,298.257223563 instead fixes the issue.</div><div><br></div><div>I'm trying to understand whether the source files incorrectly specifies that it's a sphere, or whether the GDAL driver is wrong, but I cannot find a GRIB reader that would show me this information about the file. What's the best way to debug this?</div><div><br></div><div>Thanks,</div><div>Simon</div></div>