Update the documentaiton

Author: xylar

conda_package/docs/ocean/streamfunction.rst

@@ -7,8 +7,18 @@ Computing the barotropic streamfunction
 ---------------------------------------
 
 The function :py:func:`mpas_tools.ocean.compute_barotropic_streamfunction()`
-computes the barotproic streamfunction at vertices on the MPAS-Ocean grid.
+computes the barotropic 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
+``timeMonthly_avg_`` prefix). The streamfunction is computed only over the
 range of (positive-down) depths provided and at the given time index.
+
+Optionally, the Gent-McWilliams bolus velocity (``normalGMBolusVelocity``) and
+the submesoscale velocity (``normalMLEvelocity``) can be included in the
+vertically integrated velocity calculation by setting the corresponding
+arguments ``include_bolus`` and ``include_submesoscale`` to ``True``.
+
+For large meshes, performance and memory usage can be improved by specifying
+the ``horiz_chunk`` argument to control the number of edges processed at once,
+and by providing a ``tmp_dir`` argument to use a temporary directory for
+intermediate files.