{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Using masks and computing weighted average" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This example is based from xarray example http://xarray.pydata.org/en/stable/examples/area_weighted_temperature.html" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Import python packages" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "xr.set_options(display_style='html')\n", "import intake\n", "import matplotlib.pyplot as plt\n", "import cartopy.crs as ccrs\n", "import numpy as np\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
pangeo-cmip6 catalog with 7674 dataset(s) from 514818 asset(s):
\n", " | unique | \n", "
---|---|
activity_id | \n", "18 | \n", "
institution_id | \n", "36 | \n", "
source_id | \n", "88 | \n", "
experiment_id | \n", "170 | \n", "
member_id | \n", "657 | \n", "
table_id | \n", "37 | \n", "
variable_id | \n", "700 | \n", "
grid_label | \n", "10 | \n", "
zstore | \n", "514818 | \n", "
dcpp_init_year | \n", "60 | \n", "
version | \n", "736 | \n", "
derived_variable_id | \n", "0 | \n", "
\n", " | activity_id | \n", "institution_id | \n", "source_id | \n", "experiment_id | \n", "member_id | \n", "table_id | \n", "variable_id | \n", "grid_label | \n", "zstore | \n", "dcpp_init_year | \n", "version | \n", "
---|---|---|---|---|---|---|---|---|---|---|---|
0 | \n", "CMIP | \n", "NCC | \n", "NorESM2-LM | \n", "historical | \n", "r1i1p1f1 | \n", "Amon | \n", "tas | \n", "gn | \n", "gs://cmip6/CMIP6/CMIP/NCC/NorESM2-LM/historica... | \n", "NaN | \n", "20190815 | \n", "
<xarray.Dataset> Size: 110MB\n", "Dimensions: (lat: 96, bnds: 2, lon: 144, member_id: 1,\n", " dcpp_init_year: 1, time: 1980)\n", "Coordinates:\n", " height float64 8B ...\n", " * lat (lat) float64 768B -90.0 -88.11 -86.21 ... 86.21 88.11 90.0\n", " lat_bnds (lat, bnds) float64 2kB dask.array<chunksize=(96, 2), meta=np.ndarray>\n", " * lon (lon) float64 1kB 0.0 2.5 5.0 7.5 ... 352.5 355.0 357.5\n", " lon_bnds (lon, bnds) float64 2kB dask.array<chunksize=(144, 2), meta=np.ndarray>\n", " * time (time) object 16kB 1850-01-16 12:00:00 ... 2014-12-16 12:...\n", " time_bnds (time, bnds) object 32kB dask.array<chunksize=(1980, 2), meta=np.ndarray>\n", " * member_id (member_id) object 8B 'r1i1p1f1'\n", " * dcpp_init_year (dcpp_init_year) float64 8B nan\n", "Dimensions without coordinates: bnds\n", "Data variables:\n", " tas (member_id, dcpp_init_year, time, lat, lon) float32 109MB dask.array<chunksize=(1, 1, 990, 96, 144), meta=np.ndarray>\n", "Attributes: (12/65)\n", " Conventions: CF-1.7 CMIP-6.2\n", " activity_id: CMIP\n", " branch_method: Hybrid-restart from year 1600-01-01 of ...\n", " branch_time: 0.0\n", " branch_time_in_child: 0.0\n", " branch_time_in_parent: 430335.0\n", " ... ...\n", " intake_esm_attrs:variable_id: tas\n", " intake_esm_attrs:grid_label: gn\n", " intake_esm_attrs:zstore: gs://cmip6/CMIP6/CMIP/NCC/NorESM2-LM/hi...\n", " intake_esm_attrs:version: 20190815\n", " intake_esm_attrs:_data_format_: zarr\n", " intake_esm_dataset_key: CMIP.NCC.NorESM2-LM.historical.Amon.gn