Correctly handle phase interpolation across branch cuts

[Former user]

Hi,

I have a question regarding the spatial interpolation of constituent phases.

I was browsing through the code looking for clues about how phase discontinuities (0 deg = 360 deg) are handled. Does the interpolation routine take care of branch cuts? Or will a point between two fes grid points 0 and 360 receive a value of 180?

Thanks a lot for your help!

Regards, Stefan Riha https://sriha.net

[Frederic Briol]

The code handles the discontinuity of the grids. To interpolate a value between 359 and 0 we use the last line and the first line of the grids.