Curvilinear coodinates
Add curvilinear coordinate grid as an auxiliary variable