Implement staggered grid production from node information only

I am interested in implementing a staggered grid construction function from node information only. For example, in ORCA12, only latitudes, longitudes, and cartesian co-ordinates are provided. Together these give the nodes. With the information of whether it is a B- or C-grid, we could code a function like: ```Julia function restore_staggered_grid_from_nodes(lambda, phi, e1t, e2t; format = "c-grid") ``` I think this probably belongs in Oceananigans and somewhere in Grids. I wanted to check what people think about this, and proposed APIs, before I implement in a PR. @glwagner @simone-silvestri @navidcy and others?