Philip Chmielowiec

For documentation, I have already set up a notebook for a user guide section that you can fill out. `uxarray/docs/user-guide/zonal-mean.ipynb` You should also include any functions in both our User...

You can use one of my work in progress sections as a reference if you'd like. https://uxarray--711.org.readthedocs.build/en/711/user-guide/topological-aggregations.html

> I attempted to run this on one of our sample meshes, but received an error. > > ```python > base_path = "../../test/meshfiles/ugrid/outCSne30/" > grid_path = base_path + "outCSne30.ug" >...

> > ```python > > base_path = "../../test/meshfiles/ugrid/outCSne30/" > > grid_path = base_path + "outCSne30.ug" > > data_path = base_path + "outCSne30_vortex.nc" > > uxds = ux.open_dataset(grid_path, data_path) > >...

In the process of review this now, I'll also add a simple ASV benchmark that we can run through this PR.

Here's an example I was able to generate  ```Python import matplotlib.pyplot as plt from matplotlib import gridspec import uxarray as ux base_path = "../../test/meshfiles/ugrid/outCSne30/" grid_path = base_path + "outCSne30.ug"...

> I noticed that the grid.node_x,grid.node_y,grid.node_z when reading in the MPAS data, are all unnormalized, is that what we expect? Yeah, we parse these directly from the MPAS grid file...

Can you try out something like this? ```Python edge_node_connectivity = uxgrid.edge_node_connectivity.values node_lon = uxgrid.node_lon.values node_lat = uxgrid.node_lat.values # Longitudes of the nodes that make up each edge in radians: (n_edge,...

@hongyuchen1030 I didn't consider the ordering. I think I can update the implementation I provided to take that into account. Looking into it!

I think this should work. ```Python def build_face_edges_lonlat_rad(face_node_conn, n_nodes_per_face, n_face, n_max_face_edges, node_lon, node_lat): face_node_conn_shift = np.roll(face_node_conn, -1, axis=1) for i, final_node_idx in enumerate(n_nodes_per_face): face_node_conn_shift[i, final_node_idx-1] = face_node_conn[i, 0] face_edge_conn =...