Reading tiled HDF_EOS grid datasets does not read all data

[ajain]

I found that there is an issue with reading HDF_EOS datasets in tiled manner using GDAL.

I am attaching a smaple dataset and a sample C++ program that highlights the issue. For this particular sample, what I found was that the block size returned by GDAL is 1200x833, and when I specify this block size to be read using IReadBlock API on the band, only the first 300x300 = 90000 values are read. The rest of the buffer is still in uninitialized state.

If I turn off tiled access (by setting bReadTile to FALSE), the dataset is read properly, because the tile size returned is 1200x1, and the data buffer has valid values.

[rouault]

Apparently related 1.8-esri commit that fixes/workaround the issue : r23057