Add docs

conda_package/docs/api.rst

@@ -225,6 +225,8 @@ Ocean Tools
    depth.compute_depth
    depth.compute_zmid
 
+   barotropic_streamfunction.compute_barotropic_streamfunction
+
 .. currentmodule:: mpas_tools.ocean.inject_bathymetry
 
 .. autosummary::

conda_package/docs/index.rst

@@ -50,6 +50,7 @@ analyzing simulations, and in other MPAS-related workflows.
    ocean/coastline_alteration
    ocean/moc
    ocean/depth
+   ocean/streamfunction
    ocean/visualization
 
 .. toctree::

conda_package/docs/ocean/streamfunction.rst

@@ -0,0 +1,14 @@
+.. _ocean_streamfunction:
+
+Computing streamfunctions
+=========================
+
+Computing the barotropic streamfunction
+---------------------------------------
+
+The function :py:func:`mpas_tools.ocean.barotropic_streamfunction.compute_barotropic_streamfunction()`
+computes the barotproic streamfunction at vertices on the MPAS-Ocean grid.
+The function takes a dataset containing an MPAS-Ocean mesh and another with
+``normalVelocity`` and ``layerThickness`` variables (possibly with a
+``timeMonthly_avg_`` prefix).  The streamfunction is computed only over the
+range of (positive-down) depths provided and at the given time index.