I am working with ROMS ocean model outputs, which have curvilinear grids. I would like to be able to save a variable (subset) from those outputs to a NetCDF file in a way that it is recognised to have a curvilinear grid when that file is accessed again (e.g. from R, with the stars package, but I may also want to use tools outside of R, e.g. gdalwarp). Is there a way to construct the NetCDF file appropriately?