@wenqing Here I've implemented the OGS5 approach to map physical to natural coordinates. But as far as I tested, calculated natural coordinates are wrong. Is the approach really working in OGS5? As you can see at the bottom of this function, I also implemented Newton-Raphson approach, which works well but of course is expensive.