Sorry that comment doesn't make sense! I'll try harder:

Firstly I do now understand why you construct new latitudes. I forgot that the ones from file aren't a monotonic axis but a mapping of latitude values for *every* grid point.

What I mean to suggest is that if you could come up with a way to construct new latitude values that does not depend on `jDirectionIncrementInDegrees` (this is not defined for a gaussian grid due to variable latitudinal spacing) then you could make this work for quasi-gaussian grids as well as quasi-regular girds, which would be nice I think.