{
"cells": [
{
"cell_type": "markdown",
"id": "14b764a5-2679-43e1-9f37-04f781fa01d1",
"metadata": {
"tags": []
},
"source": [
"# Barotropic Streamfunction in MOM5/MOM6\n",
"\n",
"The barotropic streamfunction ($\\psi$) is obtained from the integration of the velocity field starting from a physical boundary at which we know the transport is zero. The difference between to streamlines is a measure of the transport between them. \n",
"\n",
"There are different ways to calculate it depending on your choice of boundary for the integration. This notebook calculates it integrating in the meridional space, starting from the Antarctic continent using the zonal velocity field:\n",
"\n",
"$$\n",
"\\psi = \\int_{y_{\\rm Antarctica}}^{y} U \\, \\mathrm{d}y ,\n",
"$$\n",
"\n",
"where $U = \\int u \\, \\mathrm{dz}$ is the depth-integrated $u$-velocity."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "577aa03f-c391-4009-a759-5dd361c83c09",
"metadata": {},
"outputs": [],
"source": [
"import cosima_cookbook as cc\n",
"import numpy as np\n",
"\n",
"import cf_xarray as cfxr\n",
"import cf_xarray.units\n",
"import pint_xarray\n",
"from pint import application_registry as ureg\n",
"\n",
"import matplotlib.path as mpath\n",
"import matplotlib.pyplot as plt\n",
"import cartopy.crs as ccrs\n",
"import cartopy.feature as cfeature\n",
"import cmocean\n",
"\n",
"import dask.distributed as dask\n",
"import warnings\n",
"warnings.simplefilter(action='ignore', category=UserWarning)"
]
},
{
"cell_type": "markdown",
"id": "cf8f3cd7-1add-41bb-91c8-12fae6376e07",
"metadata": {},
"source": [
"Initialise a dask client"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "2f25a245-9737-4444-8c53-0b91ec9a5c06",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"