diff --git a/Examples/Along_Isobath_Average.ipynb b/Examples/Along_Isobath_Average.ipynb
new file mode 100644
index 00000000..277b6b2d
--- /dev/null
+++ b/Examples/Along_Isobath_Average.ipynb
@@ -0,0 +1,1700 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "id": "723a3b22-c5f9-4a34-b278-d36e29e0c296",
+ "metadata": {},
+ "source": [
+ "# Along-isobath average\n",
+ "\n",
+ "This code shows how to do an along-isobath averaging for scalar properties (and vector properties not rotated to along- and across-contour) in ACCESS-OM2-01.\n",
+ "\n",
+ "This code is relevant for Antarctic shelf analyses where a zonal average over a section of cross-shelf regions is needed.\n",
+ "\n",
+ "### Some caveats: \n",
+ "\n",
+ "1) If there are deep valleys or troughs around the coastline, those will be averaged with deeper isobaths further offshore. Similarly, if there are seamounts off-shore, these will be averaged with shallower regions on the continental slope. Careful adjustment of the mask is necessary to ensure your results are interpretable. \n",
+ "\n",
+ "Please also inspect the ht_coarse and hu_coarse variables derived from the variables of interest. \n",
+ "\n",
+ "2) More caveats should be added here! \n"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "id": "a476d924-8b6d-4c03-801b-386e07944e2f",
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "# Load necessary modules\n",
+ "import cosima_cookbook as cc\n",
+ "import matplotlib.pyplot as plt\n",
+ "import numpy as np\n",
+ "import xarray as xr\n",
+ "from xhistogram.xarray import histogram\n",
+ "from dask.distributed import Client"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "id": "e3cc23c9-8622-4a9a-9b5f-2313a8553955",
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "