{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Manipulation of CMIP6 model data using Pangeo catalog" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Import python packages" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "xr.set_options(display_style='html')\n", "import intake\n", "import cftime\n", "import matplotlib.pyplot as plt\n", "import cartopy.crs as ccrs\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Open CMIP6 online catalog" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
pangeo-cmip6 catalog with 7693 dataset(s) from 518834 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", "709 | \n", "
grid_label | \n", "10 | \n", "
zstore | \n", "518834 | \n", "
dcpp_init_year | \n", "60 | \n", "
version | \n", "718 | \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", "NCAR | \n", "CESM2 | \n", "historical | \n", "r1i1p1f1 | \n", "Amon | \n", "tas | \n", "gn | \n", "gs://cmip6/CMIP6/CMIP/NCAR/CESM2/historical/r1... | \n", "NaN | \n", "20190308 | \n", "
<xarray.Dataset>\n", "Dimensions: (lat: 192, lon: 288, member_id: 1, nbnd: 2, time: 1980)\n", "Coordinates:\n", " * lat (lat) float64 -90.0 -89.06 -88.12 -87.17 ... 88.12 89.06 90.0\n", " lat_bnds (lat, nbnd) float32 dask.array<chunksize=(192, 2), meta=np.ndarray>\n", " * lon (lon) float64 0.0 1.25 2.5 3.75 5.0 ... 355.0 356.2 357.5 358.8\n", " lon_bnds (lon, nbnd) float32 dask.array<chunksize=(288, 2), meta=np.ndarray>\n", " * time (time) object 1850-01-15 12:00:00 ... 2014-12-15 12:00:00\n", " time_bnds (time, nbnd) object dask.array<chunksize=(1980, 2), meta=np.ndarray>\n", " * member_id (member_id) <U8 'r1i1p1f1'\n", "Dimensions without coordinates: nbnd\n", "Data variables:\n", " tas (member_id, time, lat, lon) float32 dask.array<chunksize=(1, 600, 192, 288), meta=np.ndarray>\n", "Attributes: (12/48)\n", " Conventions: CF-1.7 CMIP-6.2\n", " activity_id: CMIP\n", " branch_method: standard\n", " branch_time_in_child: 674885.0\n", " branch_time_in_parent: 219000.0\n", " case_id: 15\n", " ... ...\n", " variable_id: tas\n", " variant_info: CMIP6 20th century experiments (1850-2014) with ...\n", " variant_label: r1i1p1f1\n", " status: 2019-10-25;created;by nhn2@columbia.edu\n", " intake_esm_varname: ['tas']\n", " intake_esm_dataset_key: CMIP.NCAR.CESM2.historical.Amon.gn
array([-90. , -89.057592, -88.115183, -87.172775, -86.230366, -85.287958,\n", " -84.34555 , -83.403141, -82.460733, -81.518325, -80.575916, -79.633508,\n", " -78.691099, -77.748691, -76.806283, -75.863874, -74.921466, -73.979058,\n", " -73.036649, -72.094241, -71.151832, -70.209424, -69.267016, -68.324607,\n", " -67.382199, -66.439791, -65.497382, -64.554974, -63.612565, -62.670157,\n", " -61.727749, -60.78534 , -59.842932, -58.900524, -57.958115, -57.015707,\n", " -56.073298, -55.13089 , -54.188482, -53.246073, -52.303665, -51.361257,\n", " -50.418848, -49.47644 , -48.534031, -47.591623, -46.649215, -45.706806,\n", " -44.764398, -43.82199 , -42.879581, -41.937173, -40.994764, -40.052356,\n", " -39.109948, -38.167539, -37.225131, -36.282723, -35.340314, -34.397906,\n", " -33.455497, -32.513089, -31.570681, -30.628272, -29.685864, -28.743455,\n", " -27.801047, -26.858639, -25.91623 , -24.973822, -24.031414, -23.089005,\n", " -22.146597, -21.204188, -20.26178 , -19.319372, -18.376963, -17.434555,\n", " -16.492147, -15.549738, -14.60733 , -13.664921, -12.722513, -11.780105,\n", " -10.837696, -9.895288, -8.95288 , -8.010471, -7.068063, -6.125654,\n", " -5.183246, -4.240838, -3.298429, -2.356021, -1.413613, -0.471204,\n", " 0.471204, 1.413613, 2.356021, 3.298429, 4.240838, 5.183246,\n", " 6.125654, 7.068063, 8.010471, 8.95288 , 9.895288, 10.837696,\n", " 11.780105, 12.722513, 13.664921, 14.60733 , 15.549738, 16.492147,\n", " 17.434555, 18.376963, 19.319372, 20.26178 , 21.204188, 22.146597,\n", " 23.089005, 24.031414, 24.973822, 25.91623 , 26.858639, 27.801047,\n", " 28.743455, 29.685864, 30.628272, 31.570681, 32.513089, 33.455497,\n", " 34.397906, 35.340314, 36.282723, 37.225131, 38.167539, 39.109948,\n", " 40.052356, 40.994764, 41.937173, 42.879581, 43.82199 , 44.764398,\n", " 45.706806, 46.649215, 47.591623, 48.534031, 49.47644 , 50.418848,\n", " 51.361257, 52.303665, 53.246073, 54.188482, 55.13089 , 56.073298,\n", " 57.015707, 57.958115, 58.900524, 59.842932, 60.78534 , 61.727749,\n", " 62.670157, 63.612565, 64.554974, 65.497382, 66.439791, 67.382199,\n", " 68.324607, 69.267016, 70.209424, 71.151832, 72.094241, 73.036649,\n", " 73.979058, 74.921466, 75.863874, 76.806283, 77.748691, 78.691099,\n", " 79.633508, 80.575916, 81.518325, 82.460733, 83.403141, 84.34555 ,\n", " 85.287958, 86.230366, 87.172775, 88.115183, 89.057592, 90. ])
\n",
"
| \n",
"\n", "\n", " | \n", "
array([ 0. , 1.25, 2.5 , ..., 356.25, 357.5 , 358.75])
\n",
"
| \n",
"\n", "\n", " | \n", "
array([cftime.DatetimeNoLeap(1850, 1, 15, 12, 0, 0, 0),\n", " cftime.DatetimeNoLeap(1850, 2, 14, 0, 0, 0, 0),\n", " cftime.DatetimeNoLeap(1850, 3, 15, 12, 0, 0, 0), ...,\n", " cftime.DatetimeNoLeap(2014, 10, 15, 12, 0, 0, 0),\n", " cftime.DatetimeNoLeap(2014, 11, 15, 0, 0, 0, 0),\n", " cftime.DatetimeNoLeap(2014, 12, 15, 12, 0, 0, 0)], dtype=object)
\n",
"
| \n",
"\n", "\n", " | \n", "
array(['r1i1p1f1'], dtype='<U8')
\n",
"
| \n",
"\n", "\n", " | \n", "