For quasi-regular grids the latitudes are fixed right? So I'm assuming you do this so that the latitudes are ascending when put into the interpolation routine?

The reason I'm interested is that if you just used the latitudes you extract in Line 636 this regularisation should work nicely for quasi-gaussian grids too...