{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Advanced Xarray\n", "\n", "\n", "*This material is adapted from the [Earth and Environmental Data Science](https://earth-env-data-science.github.io/intro.html), from Ryan Abernathey (Columbia University)*.\n", "\n", "\n", "In this notebook, we cover some more advanced aspects of `Xarray`.\n", "\n", "## Groupby\n", "\n", "`Xarray` copies Pandas' very useful groupby functionality, enabling the \"split / apply / combine\" workflow on xarray DataArrays and Datasets. In the first part of the notebook, we will learn to use groupby by analyzing sea-surface temperature data." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "import xarray as xr\n", "\n", "from matplotlib import pyplot as plt\n", "#%config InlineBackend.figure_format = 'retina'\n", "plt.ion() # To trigger the interactive inline mode\n", "\n", "plt.rcParams['figure.figsize'] = (6,5)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "First we load a dataset. We will use the [NOAA Extended Reconstructed Sea Surface Temperature (ERSST) v5](https://www.ncdc.noaa.gov/data-access/marineocean-data/extended-reconstructed-sea-surface-temperature-ersst-v5) product, a widely used and trusted gridded compilation of of historical data going back to 1854.\n", "\n", "Since the data is provided via an [OPeNDAP](https://en.wikipedia.org/wiki/OPeNDAP) server, we can load it directly without downloading anything:" ] }, { "cell_type": "code", "execution_count": 41, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.Dataset>\n", "Dimensions: (lat: 89, lon: 180, time: 708)\n", "Coordinates:\n", " * lat (lat) float32 88.0 86.0 84.0 82.0 80.0 ... -82.0 -84.0 -86.0 -88.0\n", " * lon (lon) float32 0.0 2.0 4.0 6.0 8.0 ... 350.0 352.0 354.0 356.0 358.0\n", " * time (time) datetime64[ns] 1960-01-01 1960-02-01 ... 2018-12-01\n", "Data variables:\n", " sst (time, lat, lon) float32 -1.8 -1.8 -1.8 -1.8 ... nan nan nan nan\n", "Attributes: (12/38)\n", " climatology: Climatology is based on 1971-2000 SST, X...\n", " description: In situ data: ICOADS2.5 before 2007 and ...\n", " keywords_vocabulary: NASA Global Change Master Directory (GCM...\n", " keywords: Earth Science > Oceans > Ocean Temperatu...\n", " instrument: Conventional thermometers\n", " source_comment: SSTs were observed by conventional therm...\n", " ... ...\n", " license: No constraints on data access or use\n", " comment: SSTs were observed by conventional therm...\n", " summary: ERSST.v5 is developed based on v4 after ...\n", " dataset_title: NOAA Extended Reconstructed SST V5\n", " data_modified: 2021-03-07\n", " DODS_EXTRA.Unlimited_Dimension: time
array([ 88., 86., 84., 82., 80., 78., 76., 74., 72., 70., 68., 66.,\n", " 64., 62., 60., 58., 56., 54., 52., 50., 48., 46., 44., 42.,\n", " 40., 38., 36., 34., 32., 30., 28., 26., 24., 22., 20., 18.,\n", " 16., 14., 12., 10., 8., 6., 4., 2., 0., -2., -4., -6.,\n", " -8., -10., -12., -14., -16., -18., -20., -22., -24., -26., -28., -30.,\n", " -32., -34., -36., -38., -40., -42., -44., -46., -48., -50., -52., -54.,\n", " -56., -58., -60., -62., -64., -66., -68., -70., -72., -74., -76., -78.,\n", " -80., -82., -84., -86., -88.], dtype=float32)
array([ 0., 2., 4., 6., 8., 10., 12., 14., 16., 18., 20., 22.,\n", " 24., 26., 28., 30., 32., 34., 36., 38., 40., 42., 44., 46.,\n", " 48., 50., 52., 54., 56., 58., 60., 62., 64., 66., 68., 70.,\n", " 72., 74., 76., 78., 80., 82., 84., 86., 88., 90., 92., 94.,\n", " 96., 98., 100., 102., 104., 106., 108., 110., 112., 114., 116., 118.,\n", " 120., 122., 124., 126., 128., 130., 132., 134., 136., 138., 140., 142.,\n", " 144., 146., 148., 150., 152., 154., 156., 158., 160., 162., 164., 166.,\n", " 168., 170., 172., 174., 176., 178., 180., 182., 184., 186., 188., 190.,\n", " 192., 194., 196., 198., 200., 202., 204., 206., 208., 210., 212., 214.,\n", " 216., 218., 220., 222., 224., 226., 228., 230., 232., 234., 236., 238.,\n", " 240., 242., 244., 246., 248., 250., 252., 254., 256., 258., 260., 262.,\n", " 264., 266., 268., 270., 272., 274., 276., 278., 280., 282., 284., 286.,\n", " 288., 290., 292., 294., 296., 298., 300., 302., 304., 306., 308., 310.,\n", " 312., 314., 316., 318., 320., 322., 324., 326., 328., 330., 332., 334.,\n", " 336., 338., 340., 342., 344., 346., 348., 350., 352., 354., 356., 358.],\n", " dtype=float32)
array(['1960-01-01T00:00:00.000000000', '1960-02-01T00:00:00.000000000',\n", " '1960-03-01T00:00:00.000000000', ..., '2018-10-01T00:00:00.000000000',\n", " '2018-11-01T00:00:00.000000000', '2018-12-01T00:00:00.000000000'],\n", " dtype='datetime64[ns]')
array([[[-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan]],\n", "\n", " [[-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan]],\n", "\n", " [[-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " ...,\n", "...\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan]],\n", "\n", " [[-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan]],\n", "\n", " [[-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan]]], dtype=float32)
<xarray.DataArray 'time' (time: 708)>\n", "array(['1960-01-01T00:00:00.000000000', '1960-02-01T00:00:00.000000000',\n", " '1960-03-01T00:00:00.000000000', ..., '2018-10-01T00:00:00.000000000',\n", " '2018-11-01T00:00:00.000000000', '2018-12-01T00:00:00.000000000'],\n", " dtype='datetime64[ns]')\n", "Coordinates:\n", " * time (time) datetime64[ns] 1960-01-01 1960-02-01 ... 2018-12-01\n", "Attributes:\n", " long_name: Time\n", " delta_t: 0000-01-00 00:00:00\n", " avg_period: 0000-01-00 00:00:00\n", " prev_avg_period: 0000-00-07 00:00:00\n", " standard_name: time\n", " axis: T\n", " actual_range: [19723. 80750.]\n", " _ChunkSizes: 1
array(['1960-01-01T00:00:00.000000000', '1960-02-01T00:00:00.000000000',\n", " '1960-03-01T00:00:00.000000000', ..., '2018-10-01T00:00:00.000000000',\n", " '2018-11-01T00:00:00.000000000', '2018-12-01T00:00:00.000000000'],\n", " dtype='datetime64[ns]')
array(['1960-01-01T00:00:00.000000000', '1960-02-01T00:00:00.000000000',\n", " '1960-03-01T00:00:00.000000000', ..., '2018-10-01T00:00:00.000000000',\n", " '2018-11-01T00:00:00.000000000', '2018-12-01T00:00:00.000000000'],\n", " dtype='datetime64[ns]')
<xarray.DataArray 'month' (time: 708)>\n", "array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6,\n", " 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,\n", " 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4,\n", " 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9,\n", " 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2,\n", " 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7,\n", " 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,\n", " 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6,\n", " 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,\n", " 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4,\n", "...\n", " 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7,\n", " 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,\n", " 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6,\n", " 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,\n", " 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4,\n", " 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9,\n", " 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2,\n", " 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7,\n", " 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,\n", " 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12])\n", "Coordinates:\n", " * time (time) datetime64[ns] 1960-01-01 1960-02-01 ... 2018-12-01
array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6,\n", " 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,\n", " 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4,\n", " 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9,\n", " 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2,\n", " 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7,\n", " 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,\n", " 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6,\n", " 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,\n", " 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4,\n", "...\n", " 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7,\n", " 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,\n", " 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6,\n", " 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,\n", " 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4,\n", " 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9,\n", " 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2,\n", " 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7,\n", " 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,\n", " 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12])
array(['1960-01-01T00:00:00.000000000', '1960-02-01T00:00:00.000000000',\n", " '1960-03-01T00:00:00.000000000', ..., '2018-10-01T00:00:00.000000000',\n", " '2018-11-01T00:00:00.000000000', '2018-12-01T00:00:00.000000000'],\n", " dtype='datetime64[ns]')
<xarray.Dataset>\n", "Dimensions: (lat: 89, lon: 180, time: 59)\n", "Coordinates:\n", " * lat (lat) float32 88.0 86.0 84.0 82.0 80.0 ... -82.0 -84.0 -86.0 -88.0\n", " * lon (lon) float32 0.0 2.0 4.0 6.0 8.0 ... 350.0 352.0 354.0 356.0 358.0\n", " * time (time) datetime64[ns] 1960-01-01 1961-01-01 ... 2018-01-01\n", "Data variables:\n", " sst (time, lat, lon) float32 -1.8 -1.8 -1.8 -1.8 ... nan nan nan nan\n", "Attributes: (12/38)\n", " climatology: Climatology is based on 1971-2000 SST, X...\n", " description: In situ data: ICOADS2.5 before 2007 and ...\n", " keywords_vocabulary: NASA Global Change Master Directory (GCM...\n", " keywords: Earth Science > Oceans > Ocean Temperatu...\n", " instrument: Conventional thermometers\n", " source_comment: SSTs were observed by conventional therm...\n", " ... ...\n", " license: No constraints on data access or use\n", " comment: SSTs were observed by conventional therm...\n", " summary: ERSST.v5 is developed based on v4 after ...\n", " dataset_title: NOAA Extended Reconstructed SST V5\n", " data_modified: 2021-03-07\n", " DODS_EXTRA.Unlimited_Dimension: time
array([ 88., 86., 84., 82., 80., 78., 76., 74., 72., 70., 68., 66.,\n", " 64., 62., 60., 58., 56., 54., 52., 50., 48., 46., 44., 42.,\n", " 40., 38., 36., 34., 32., 30., 28., 26., 24., 22., 20., 18.,\n", " 16., 14., 12., 10., 8., 6., 4., 2., 0., -2., -4., -6.,\n", " -8., -10., -12., -14., -16., -18., -20., -22., -24., -26., -28., -30.,\n", " -32., -34., -36., -38., -40., -42., -44., -46., -48., -50., -52., -54.,\n", " -56., -58., -60., -62., -64., -66., -68., -70., -72., -74., -76., -78.,\n", " -80., -82., -84., -86., -88.], dtype=float32)
array([ 0., 2., 4., 6., 8., 10., 12., 14., 16., 18., 20., 22.,\n", " 24., 26., 28., 30., 32., 34., 36., 38., 40., 42., 44., 46.,\n", " 48., 50., 52., 54., 56., 58., 60., 62., 64., 66., 68., 70.,\n", " 72., 74., 76., 78., 80., 82., 84., 86., 88., 90., 92., 94.,\n", " 96., 98., 100., 102., 104., 106., 108., 110., 112., 114., 116., 118.,\n", " 120., 122., 124., 126., 128., 130., 132., 134., 136., 138., 140., 142.,\n", " 144., 146., 148., 150., 152., 154., 156., 158., 160., 162., 164., 166.,\n", " 168., 170., 172., 174., 176., 178., 180., 182., 184., 186., 188., 190.,\n", " 192., 194., 196., 198., 200., 202., 204., 206., 208., 210., 212., 214.,\n", " 216., 218., 220., 222., 224., 226., 228., 230., 232., 234., 236., 238.,\n", " 240., 242., 244., 246., 248., 250., 252., 254., 256., 258., 260., 262.,\n", " 264., 266., 268., 270., 272., 274., 276., 278., 280., 282., 284., 286.,\n", " 288., 290., 292., 294., 296., 298., 300., 302., 304., 306., 308., 310.,\n", " 312., 314., 316., 318., 320., 322., 324., 326., 328., 330., 332., 334.,\n", " 336., 338., 340., 342., 344., 346., 348., 350., 352., 354., 356., 358.],\n", " dtype=float32)
array(['1960-01-01T00:00:00.000000000', '1961-01-01T00:00:00.000000000',\n", " '1962-01-01T00:00:00.000000000', '1963-01-01T00:00:00.000000000',\n", " '1964-01-01T00:00:00.000000000', '1965-01-01T00:00:00.000000000',\n", " '1966-01-01T00:00:00.000000000', '1967-01-01T00:00:00.000000000',\n", " '1968-01-01T00:00:00.000000000', '1969-01-01T00:00:00.000000000',\n", " '1970-01-01T00:00:00.000000000', '1971-01-01T00:00:00.000000000',\n", " '1972-01-01T00:00:00.000000000', '1973-01-01T00:00:00.000000000',\n", " '1974-01-01T00:00:00.000000000', '1975-01-01T00:00:00.000000000',\n", " '1976-01-01T00:00:00.000000000', '1977-01-01T00:00:00.000000000',\n", " '1978-01-01T00:00:00.000000000', '1979-01-01T00:00:00.000000000',\n", " '1980-01-01T00:00:00.000000000', '1981-01-01T00:00:00.000000000',\n", " '1982-01-01T00:00:00.000000000', '1983-01-01T00:00:00.000000000',\n", " '1984-01-01T00:00:00.000000000', '1985-01-01T00:00:00.000000000',\n", " '1986-01-01T00:00:00.000000000', '1987-01-01T00:00:00.000000000',\n", " '1988-01-01T00:00:00.000000000', '1989-01-01T00:00:00.000000000',\n", " '1990-01-01T00:00:00.000000000', '1991-01-01T00:00:00.000000000',\n", " '1992-01-01T00:00:00.000000000', '1993-01-01T00:00:00.000000000',\n", " '1994-01-01T00:00:00.000000000', '1995-01-01T00:00:00.000000000',\n", " '1996-01-01T00:00:00.000000000', '1997-01-01T00:00:00.000000000',\n", " '1998-01-01T00:00:00.000000000', '1999-01-01T00:00:00.000000000',\n", " '2000-01-01T00:00:00.000000000', '2001-01-01T00:00:00.000000000',\n", " '2002-01-01T00:00:00.000000000', '2003-01-01T00:00:00.000000000',\n", " '2004-01-01T00:00:00.000000000', '2005-01-01T00:00:00.000000000',\n", " '2006-01-01T00:00:00.000000000', '2007-01-01T00:00:00.000000000',\n", " '2008-01-01T00:00:00.000000000', '2009-01-01T00:00:00.000000000',\n", " '2010-01-01T00:00:00.000000000', '2011-01-01T00:00:00.000000000',\n", " '2012-01-01T00:00:00.000000000', '2013-01-01T00:00:00.000000000',\n", " '2014-01-01T00:00:00.000000000', '2015-01-01T00:00:00.000000000',\n", " '2016-01-01T00:00:00.000000000', '2017-01-01T00:00:00.000000000',\n", " '2018-01-01T00:00:00.000000000'], dtype='datetime64[ns]')
array([[[-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan]],\n", "\n", " [[-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan]],\n", "\n", " [[-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " ...,\n", "...\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan]],\n", "\n", " [[-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan]],\n", "\n", " [[-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " [-1.8, -1.8, -1.8, ..., -1.8, -1.8, -1.8],\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan]]], dtype=float32)
<xarray.Dataset>\n", "Dimensions: (month: 12)\n", "Coordinates:\n", " * month (month) int64 1 2 3 4 5 6 7 8 9 10 11 12\n", "Data variables:\n", " sst (month) float32 13.66 13.77 13.76 13.68 ... 13.98 13.69 13.51 13.53
array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12])
array([13.659641, 13.768647, 13.76488 , 13.684034, 13.642146, 13.713043,\n", " 13.921847, 14.093956, 13.982147, 13.691116, 13.506494, 13.529454],\n", " dtype=float32)
<xarray.Dataset>\n", "Dimensions: (lat: 89, lon: 180, month: 12)\n", "Coordinates:\n", " * lat (lat) float32 88.0 86.0 84.0 82.0 80.0 ... -82.0 -84.0 -86.0 -88.0\n", " * lon (lon) float32 0.0 2.0 4.0 6.0 8.0 ... 350.0 352.0 354.0 356.0 358.0\n", " * month (month) int64 1 2 3 4 5 6 7 8 9 10 11 12\n", "Data variables:\n", " sst (month, lat, lon) float32 -1.8 -1.8 -1.8 -1.8 ... nan nan nan nan
array([ 88., 86., 84., 82., 80., 78., 76., 74., 72., 70., 68., 66.,\n", " 64., 62., 60., 58., 56., 54., 52., 50., 48., 46., 44., 42.,\n", " 40., 38., 36., 34., 32., 30., 28., 26., 24., 22., 20., 18.,\n", " 16., 14., 12., 10., 8., 6., 4., 2., 0., -2., -4., -6.,\n", " -8., -10., -12., -14., -16., -18., -20., -22., -24., -26., -28., -30.,\n", " -32., -34., -36., -38., -40., -42., -44., -46., -48., -50., -52., -54.,\n", " -56., -58., -60., -62., -64., -66., -68., -70., -72., -74., -76., -78.,\n", " -80., -82., -84., -86., -88.], dtype=float32)
array([ 0., 2., 4., 6., 8., 10., 12., 14., 16., 18., 20., 22.,\n", " 24., 26., 28., 30., 32., 34., 36., 38., 40., 42., 44., 46.,\n", " 48., 50., 52., 54., 56., 58., 60., 62., 64., 66., 68., 70.,\n", " 72., 74., 76., 78., 80., 82., 84., 86., 88., 90., 92., 94.,\n", " 96., 98., 100., 102., 104., 106., 108., 110., 112., 114., 116., 118.,\n", " 120., 122., 124., 126., 128., 130., 132., 134., 136., 138., 140., 142.,\n", " 144., 146., 148., 150., 152., 154., 156., 158., 160., 162., 164., 166.,\n", " 168., 170., 172., 174., 176., 178., 180., 182., 184., 186., 188., 190.,\n", " 192., 194., 196., 198., 200., 202., 204., 206., 208., 210., 212., 214.,\n", " 216., 218., 220., 222., 224., 226., 228., 230., 232., 234., 236., 238.,\n", " 240., 242., 244., 246., 248., 250., 252., 254., 256., 258., 260., 262.,\n", " 264., 266., 268., 270., 272., 274., 276., 278., 280., 282., 284., 286.,\n", " 288., 290., 292., 294., 296., 298., 300., 302., 304., 306., 308., 310.,\n", " 312., 314., 316., 318., 320., 322., 324., 326., 328., 330., 332., 334.,\n", " 336., 338., 340., 342., 344., 346., 348., 350., 352., 354., 356., 358.],\n", " dtype=float32)
array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12])
array([[[-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", " [-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", " [-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", " [-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", " [-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", "...\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[-1.7995025, -1.7995973, -1.7998415, ..., -1.7997988,\n", " -1.7996519, -1.7995045],\n", " [-1.7995876, -1.7997634, -1.8000009, ..., -1.8000009,\n", " -1.7998358, -1.7996247],\n", " [-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]]], dtype=float32)
<xarray.Dataset>\n", "Dimensions: (lat: 89, lon: 180, month: 12)\n", "Coordinates:\n", " * lat (lat) float32 88.0 86.0 84.0 82.0 80.0 ... -82.0 -84.0 -86.0 -88.0\n", " * lon (lon) float32 0.0 2.0 4.0 6.0 8.0 ... 350.0 352.0 354.0 356.0 358.0\n", " * month (month) int64 1 2 3 4 5 6 7 8 9 10 11 12\n", "Data variables:\n", " sst (month, lat, lon) float32 -1.8 -1.8 -1.8 -1.8 ... nan nan nan nan
array([ 88., 86., 84., 82., 80., 78., 76., 74., 72., 70., 68., 66.,\n", " 64., 62., 60., 58., 56., 54., 52., 50., 48., 46., 44., 42.,\n", " 40., 38., 36., 34., 32., 30., 28., 26., 24., 22., 20., 18.,\n", " 16., 14., 12., 10., 8., 6., 4., 2., 0., -2., -4., -6.,\n", " -8., -10., -12., -14., -16., -18., -20., -22., -24., -26., -28., -30.,\n", " -32., -34., -36., -38., -40., -42., -44., -46., -48., -50., -52., -54.,\n", " -56., -58., -60., -62., -64., -66., -68., -70., -72., -74., -76., -78.,\n", " -80., -82., -84., -86., -88.], dtype=float32)
array([ 0., 2., 4., 6., 8., 10., 12., 14., 16., 18., 20., 22.,\n", " 24., 26., 28., 30., 32., 34., 36., 38., 40., 42., 44., 46.,\n", " 48., 50., 52., 54., 56., 58., 60., 62., 64., 66., 68., 70.,\n", " 72., 74., 76., 78., 80., 82., 84., 86., 88., 90., 92., 94.,\n", " 96., 98., 100., 102., 104., 106., 108., 110., 112., 114., 116., 118.,\n", " 120., 122., 124., 126., 128., 130., 132., 134., 136., 138., 140., 142.,\n", " 144., 146., 148., 150., 152., 154., 156., 158., 160., 162., 164., 166.,\n", " 168., 170., 172., 174., 176., 178., 180., 182., 184., 186., 188., 190.,\n", " 192., 194., 196., 198., 200., 202., 204., 206., 208., 210., 212., 214.,\n", " 216., 218., 220., 222., 224., 226., 228., 230., 232., 234., 236., 238.,\n", " 240., 242., 244., 246., 248., 250., 252., 254., 256., 258., 260., 262.,\n", " 264., 266., 268., 270., 272., 274., 276., 278., 280., 282., 284., 286.,\n", " 288., 290., 292., 294., 296., 298., 300., 302., 304., 306., 308., 310.,\n", " 312., 314., 316., 318., 320., 322., 324., 326., 328., 330., 332., 334.,\n", " 336., 338., 340., 342., 344., 346., 348., 350., 352., 354., 356., 358.],\n", " dtype=float32)
array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12])
array([[[-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", " [-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", " [-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", " [-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", " [-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", "...\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[-1.7995025, -1.7995973, -1.7998415, ..., -1.7997988,\n", " -1.7996519, -1.7995045],\n", " [-1.7995876, -1.7997634, -1.8000009, ..., -1.8000009,\n", " -1.7998358, -1.7996247],\n", " [-1.8000009, -1.8000009, -1.8000009, ..., -1.8000009,\n", " -1.8000009, -1.8000009],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]]], dtype=float32)
<xarray.Dataset>\n", "Dimensions: (lat: 89, lon: 180, time: 708)\n", "Coordinates:\n", " * lat (lat) float32 88.0 86.0 84.0 82.0 80.0 ... -82.0 -84.0 -86.0 -88.0\n", " * lon (lon) float32 0.0 2.0 4.0 6.0 8.0 ... 350.0 352.0 354.0 356.0 358.0\n", " * time (time) datetime64[ns] 1960-01-01 1960-02-01 ... 2018-12-01\n", "Data variables:\n", " sst (time, lat, lon) float32 9.537e-07 9.537e-07 9.537e-07 ... nan nan
array([ 88., 86., 84., 82., 80., 78., 76., 74., 72., 70., 68., 66.,\n", " 64., 62., 60., 58., 56., 54., 52., 50., 48., 46., 44., 42.,\n", " 40., 38., 36., 34., 32., 30., 28., 26., 24., 22., 20., 18.,\n", " 16., 14., 12., 10., 8., 6., 4., 2., 0., -2., -4., -6.,\n", " -8., -10., -12., -14., -16., -18., -20., -22., -24., -26., -28., -30.,\n", " -32., -34., -36., -38., -40., -42., -44., -46., -48., -50., -52., -54.,\n", " -56., -58., -60., -62., -64., -66., -68., -70., -72., -74., -76., -78.,\n", " -80., -82., -84., -86., -88.], dtype=float32)
array([ 0., 2., 4., 6., 8., 10., 12., 14., 16., 18., 20., 22.,\n", " 24., 26., 28., 30., 32., 34., 36., 38., 40., 42., 44., 46.,\n", " 48., 50., 52., 54., 56., 58., 60., 62., 64., 66., 68., 70.,\n", " 72., 74., 76., 78., 80., 82., 84., 86., 88., 90., 92., 94.,\n", " 96., 98., 100., 102., 104., 106., 108., 110., 112., 114., 116., 118.,\n", " 120., 122., 124., 126., 128., 130., 132., 134., 136., 138., 140., 142.,\n", " 144., 146., 148., 150., 152., 154., 156., 158., 160., 162., 164., 166.,\n", " 168., 170., 172., 174., 176., 178., 180., 182., 184., 186., 188., 190.,\n", " 192., 194., 196., 198., 200., 202., 204., 206., 208., 210., 212., 214.,\n", " 216., 218., 220., 222., 224., 226., 228., 230., 232., 234., 236., 238.,\n", " 240., 242., 244., 246., 248., 250., 252., 254., 256., 258., 260., 262.,\n", " 264., 266., 268., 270., 272., 274., 276., 278., 280., 282., 284., 286.,\n", " 288., 290., 292., 294., 296., 298., 300., 302., 304., 306., 308., 310.,\n", " 312., 314., 316., 318., 320., 322., 324., 326., 328., 330., 332., 334.,\n", " 336., 338., 340., 342., 344., 346., 348., 350., 352., 354., 356., 358.],\n", " dtype=float32)
array(['1960-01-01T00:00:00.000000000', '1960-02-01T00:00:00.000000000',\n", " '1960-03-01T00:00:00.000000000', ..., '2018-10-01T00:00:00.000000000',\n", " '2018-11-01T00:00:00.000000000', '2018-12-01T00:00:00.000000000'],\n", " dtype='datetime64[ns]')
array([[[ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", " [ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", " [ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", " ...,\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan]],\n", "\n", " [[ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", " [ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", " [ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", "...\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan]],\n", "\n", " [[-4.9746037e-04, -4.0268898e-04, -1.5842915e-04, ...,\n", " -2.0110607e-04, -3.4809113e-04, -4.9543381e-04],\n", " [-4.1234493e-04, -2.3651123e-04, 9.5367432e-07, ...,\n", " 9.5367432e-07, -1.6415119e-04, -3.7527084e-04],\n", " [ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", " ...,\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan]]], dtype=float32)
<xarray.Dataset>\n", "Dimensions: (lat: 89, lon: 180, time: 708)\n", "Coordinates:\n", " * lat (lat) float32 88.0 86.0 84.0 82.0 80.0 ... -82.0 -84.0 -86.0 -88.0\n", " * lon (lon) float32 0.0 2.0 4.0 6.0 8.0 ... 350.0 352.0 354.0 356.0 358.0\n", " * time (time) datetime64[ns] 1960-01-01 1960-02-01 ... 2018-12-01\n", " month (time) int64 1 2 3 4 5 6 7 8 9 10 11 ... 2 3 4 5 6 7 8 9 10 11 12\n", "Data variables:\n", " sst (time, lat, lon) float32 9.537e-07 9.537e-07 9.537e-07 ... nan nan
array([ 88., 86., 84., 82., 80., 78., 76., 74., 72., 70., 68., 66.,\n", " 64., 62., 60., 58., 56., 54., 52., 50., 48., 46., 44., 42.,\n", " 40., 38., 36., 34., 32., 30., 28., 26., 24., 22., 20., 18.,\n", " 16., 14., 12., 10., 8., 6., 4., 2., 0., -2., -4., -6.,\n", " -8., -10., -12., -14., -16., -18., -20., -22., -24., -26., -28., -30.,\n", " -32., -34., -36., -38., -40., -42., -44., -46., -48., -50., -52., -54.,\n", " -56., -58., -60., -62., -64., -66., -68., -70., -72., -74., -76., -78.,\n", " -80., -82., -84., -86., -88.], dtype=float32)
array([ 0., 2., 4., 6., 8., 10., 12., 14., 16., 18., 20., 22.,\n", " 24., 26., 28., 30., 32., 34., 36., 38., 40., 42., 44., 46.,\n", " 48., 50., 52., 54., 56., 58., 60., 62., 64., 66., 68., 70.,\n", " 72., 74., 76., 78., 80., 82., 84., 86., 88., 90., 92., 94.,\n", " 96., 98., 100., 102., 104., 106., 108., 110., 112., 114., 116., 118.,\n", " 120., 122., 124., 126., 128., 130., 132., 134., 136., 138., 140., 142.,\n", " 144., 146., 148., 150., 152., 154., 156., 158., 160., 162., 164., 166.,\n", " 168., 170., 172., 174., 176., 178., 180., 182., 184., 186., 188., 190.,\n", " 192., 194., 196., 198., 200., 202., 204., 206., 208., 210., 212., 214.,\n", " 216., 218., 220., 222., 224., 226., 228., 230., 232., 234., 236., 238.,\n", " 240., 242., 244., 246., 248., 250., 252., 254., 256., 258., 260., 262.,\n", " 264., 266., 268., 270., 272., 274., 276., 278., 280., 282., 284., 286.,\n", " 288., 290., 292., 294., 296., 298., 300., 302., 304., 306., 308., 310.,\n", " 312., 314., 316., 318., 320., 322., 324., 326., 328., 330., 332., 334.,\n", " 336., 338., 340., 342., 344., 346., 348., 350., 352., 354., 356., 358.],\n", " dtype=float32)
array(['1960-01-01T00:00:00.000000000', '1960-02-01T00:00:00.000000000',\n", " '1960-03-01T00:00:00.000000000', ..., '2018-10-01T00:00:00.000000000',\n", " '2018-11-01T00:00:00.000000000', '2018-12-01T00:00:00.000000000'],\n", " dtype='datetime64[ns]')
array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6,\n", " 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,\n", " 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4,\n", " 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9,\n", " 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2,\n", " 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7,\n", " 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,\n", " 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6,\n", " 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,\n", " 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4,\n", "...\n", " 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7,\n", " 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,\n", " 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6,\n", " 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,\n", " 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4,\n", " 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9,\n", " 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2,\n", " 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7,\n", " 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,\n", " 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12])
array([[[ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", " [ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", " [ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", " ...,\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan]],\n", "\n", " [[ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", " [ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", " [ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", "...\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan]],\n", "\n", " [[-4.9746037e-04, -4.0268898e-04, -1.5842915e-04, ...,\n", " -2.0110607e-04, -3.4809113e-04, -4.9543381e-04],\n", " [-4.1234493e-04, -2.3651123e-04, 9.5367432e-07, ...,\n", " 9.5367432e-07, -1.6415119e-04, -3.7527084e-04],\n", " [ 9.5367432e-07, 9.5367432e-07, 9.5367432e-07, ...,\n", " 9.5367432e-07, 9.5367432e-07, 9.5367432e-07],\n", " ...,\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan]]], dtype=float32)
<xarray.Dataset>\n", "Dimensions: (lat: 89, lon: 180, time: 13)\n", "Coordinates:\n", " * time (time) datetime64[ns] 1960-12-31 1965-12-31 ... 2020-12-31\n", " * lat (lat) float32 88.0 86.0 84.0 82.0 80.0 ... -82.0 -84.0 -86.0 -88.0\n", " * lon (lon) float32 0.0 2.0 4.0 6.0 8.0 ... 350.0 352.0 354.0 356.0 358.0\n", "Data variables:\n", " sst (time, lat, lon) float32 -0.0005707 -0.0005493 ... nan nan
array(['1960-12-31T00:00:00.000000000', '1965-12-31T00:00:00.000000000',\n", " '1970-12-31T00:00:00.000000000', '1975-12-31T00:00:00.000000000',\n", " '1980-12-31T00:00:00.000000000', '1985-12-31T00:00:00.000000000',\n", " '1990-12-31T00:00:00.000000000', '1995-12-31T00:00:00.000000000',\n", " '2000-12-31T00:00:00.000000000', '2005-12-31T00:00:00.000000000',\n", " '2010-12-31T00:00:00.000000000', '2015-12-31T00:00:00.000000000',\n", " '2020-12-31T00:00:00.000000000'], dtype='datetime64[ns]')
array([ 88., 86., 84., 82., 80., 78., 76., 74., 72., 70., 68., 66.,\n", " 64., 62., 60., 58., 56., 54., 52., 50., 48., 46., 44., 42.,\n", " 40., 38., 36., 34., 32., 30., 28., 26., 24., 22., 20., 18.,\n", " 16., 14., 12., 10., 8., 6., 4., 2., 0., -2., -4., -6.,\n", " -8., -10., -12., -14., -16., -18., -20., -22., -24., -26., -28., -30.,\n", " -32., -34., -36., -38., -40., -42., -44., -46., -48., -50., -52., -54.,\n", " -56., -58., -60., -62., -64., -66., -68., -70., -72., -74., -76., -78.,\n", " -80., -82., -84., -86., -88.], dtype=float32)
array([ 0., 2., 4., 6., 8., 10., 12., 14., 16., 18., 20., 22.,\n", " 24., 26., 28., 30., 32., 34., 36., 38., 40., 42., 44., 46.,\n", " 48., 50., 52., 54., 56., 58., 60., 62., 64., 66., 68., 70.,\n", " 72., 74., 76., 78., 80., 82., 84., 86., 88., 90., 92., 94.,\n", " 96., 98., 100., 102., 104., 106., 108., 110., 112., 114., 116., 118.,\n", " 120., 122., 124., 126., 128., 130., 132., 134., 136., 138., 140., 142.,\n", " 144., 146., 148., 150., 152., 154., 156., 158., 160., 162., 164., 166.,\n", " 168., 170., 172., 174., 176., 178., 180., 182., 184., 186., 188., 190.,\n", " 192., 194., 196., 198., 200., 202., 204., 206., 208., 210., 212., 214.,\n", " 216., 218., 220., 222., 224., 226., 228., 230., 232., 234., 236., 238.,\n", " 240., 242., 244., 246., 248., 250., 252., 254., 256., 258., 260., 262.,\n", " 264., 266., 268., 270., 272., 274., 276., 278., 280., 282., 284., 286.,\n", " 288., 290., 292., 294., 296., 298., 300., 302., 304., 306., 308., 310.,\n", " 312., 314., 316., 318., 320., 322., 324., 326., 328., 330., 332., 334.,\n", " 336., 338., 340., 342., 344., 346., 348., 350., 352., 354., 356., 358.],\n", " dtype=float32)
array([[[-5.70684671e-04, -5.49256802e-04, -5.38418710e-04, ...,\n", " -4.15881485e-04, -5.20139933e-04, -5.83092391e-04],\n", " [-5.15530526e-04, -5.89827716e-04, -6.42200292e-04, ...,\n", " -1.39196709e-04, -3.24239343e-04, -4.26093728e-04],\n", " [-1.40964985e-05, -1.79181501e-04, -2.83718109e-04, ...,\n", " -1.24106809e-04, -7.90158883e-05, -3.50972005e-05],\n", " ...,\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan]],\n", "\n", " [[-5.70684671e-04, -5.49256802e-04, -5.38418710e-04, ...,\n", " -4.15881485e-04, -5.20139933e-04, -5.83092391e-04],\n", " [-5.15530526e-04, -5.89827716e-04, -6.42200292e-04, ...,\n", " -1.39196709e-04, -3.24239343e-04, -4.26093728e-04],\n", " [-1.40964985e-05, -1.79181501e-04, -2.83718109e-04, ...,\n", " -1.24106809e-04, -7.90158883e-05, -3.50972005e-05],\n", "...\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan]],\n", "\n", " [[ 1.06710559e-02, 1.02710230e-02, 1.00688301e-02, ...,\n", " 7.78120104e-03, 9.72755440e-03, 1.09023023e-02],\n", " [ 9.64192580e-03, 1.10286605e-02, 1.20063955e-02, ...,\n", " 2.61675660e-03, 6.07069675e-03, 7.97217432e-03],\n", " [ 2.81734596e-04, 3.36324517e-03, 5.31468447e-03, ...,\n", " 2.31115357e-03, 1.49352022e-03, 6.73671544e-04],\n", " ...,\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan],\n", " [ nan, nan, nan, ...,\n", " nan, nan, nan]]],\n", " dtype=float32)
<xarray.Dataset>\n", "Dimensions: (lat: 89, lon: 180, time: 708)\n", "Coordinates:\n", " * lat (lat) float32 88.0 86.0 84.0 82.0 80.0 ... -82.0 -84.0 -86.0 -88.0\n", " * lon (lon) float32 0.0 2.0 4.0 6.0 8.0 ... 350.0 352.0 354.0 356.0 358.0\n", " * time (time) datetime64[ns] 1960-01-01 1960-02-01 ... 2018-12-01\n", " month (time) int64 1 2 3 4 5 6 7 8 9 10 11 ... 2 3 4 5 6 7 8 9 10 11 12\n", "Data variables:\n", " sst (time, lat, lon) float64 nan nan nan nan nan ... nan nan nan nan
array([ 88., 86., 84., 82., 80., 78., 76., 74., 72., 70., 68., 66.,\n", " 64., 62., 60., 58., 56., 54., 52., 50., 48., 46., 44., 42.,\n", " 40., 38., 36., 34., 32., 30., 28., 26., 24., 22., 20., 18.,\n", " 16., 14., 12., 10., 8., 6., 4., 2., 0., -2., -4., -6.,\n", " -8., -10., -12., -14., -16., -18., -20., -22., -24., -26., -28., -30.,\n", " -32., -34., -36., -38., -40., -42., -44., -46., -48., -50., -52., -54.,\n", " -56., -58., -60., -62., -64., -66., -68., -70., -72., -74., -76., -78.,\n", " -80., -82., -84., -86., -88.], dtype=float32)
array([ 0., 2., 4., 6., 8., 10., 12., 14., 16., 18., 20., 22.,\n", " 24., 26., 28., 30., 32., 34., 36., 38., 40., 42., 44., 46.,\n", " 48., 50., 52., 54., 56., 58., 60., 62., 64., 66., 68., 70.,\n", " 72., 74., 76., 78., 80., 82., 84., 86., 88., 90., 92., 94.,\n", " 96., 98., 100., 102., 104., 106., 108., 110., 112., 114., 116., 118.,\n", " 120., 122., 124., 126., 128., 130., 132., 134., 136., 138., 140., 142.,\n", " 144., 146., 148., 150., 152., 154., 156., 158., 160., 162., 164., 166.,\n", " 168., 170., 172., 174., 176., 178., 180., 182., 184., 186., 188., 190.,\n", " 192., 194., 196., 198., 200., 202., 204., 206., 208., 210., 212., 214.,\n", " 216., 218., 220., 222., 224., 226., 228., 230., 232., 234., 236., 238.,\n", " 240., 242., 244., 246., 248., 250., 252., 254., 256., 258., 260., 262.,\n", " 264., 266., 268., 270., 272., 274., 276., 278., 280., 282., 284., 286.,\n", " 288., 290., 292., 294., 296., 298., 300., 302., 304., 306., 308., 310.,\n", " 312., 314., 316., 318., 320., 322., 324., 326., 328., 330., 332., 334.,\n", " 336., 338., 340., 342., 344., 346., 348., 350., 352., 354., 356., 358.],\n", " dtype=float32)
array(['1960-01-01T00:00:00.000000000', '1960-02-01T00:00:00.000000000',\n", " '1960-03-01T00:00:00.000000000', ..., '2018-10-01T00:00:00.000000000',\n", " '2018-11-01T00:00:00.000000000', '2018-12-01T00:00:00.000000000'],\n", " dtype='datetime64[ns]')
array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6,\n", " 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,\n", " 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4,\n", " 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9,\n", " 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2,\n", " 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7,\n", " 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,\n", " 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6,\n", " 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,\n", " 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4,\n", "...\n", " 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7,\n", " 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,\n", " 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6,\n", " 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,\n", " 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4,\n", " 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9,\n", " 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2,\n", " 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7,\n", " 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,\n", " 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5,\n", " 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,\n", " 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3,\n", " 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8,\n", " 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1,\n", " 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12])
array([[[nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " ...,\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan]],\n", "\n", " [[nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " ...,\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan]],\n", "\n", " [[nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " ...,\n", "...\n", " ...,\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan]],\n", "\n", " [[nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " ...,\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan]],\n", "\n", " [[nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " ...,\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan]]])
<xarray.Dataset>\n", "Dimensions: (X: 360, Y: 180, Z: 33)\n", "Coordinates:\n", " * Z (Z) float32 0.0 10.0 20.0 30.0 50.0 ... 4e+03 4.5e+03 5e+03 5.5e+03\n", " * X (X) float32 0.5 1.5 2.5 3.5 4.5 ... 355.5 356.5 357.5 358.5 359.5\n", " * Y (Y) float32 -89.5 -88.5 -87.5 -86.5 -85.5 ... 86.5 87.5 88.5 89.5\n", "Data variables:\n", " basin (Z, Y, X) float32 ...\n", "Attributes:\n", " Conventions: IRIDL
array([ 0., 10., 20., 30., 50., 75., 100., 125., 150., 200.,\n", " 250., 300., 400., 500., 600., 700., 800., 900., 1000., 1100.,\n", " 1200., 1300., 1400., 1500., 1750., 2000., 2500., 3000., 3500., 4000.,\n", " 4500., 5000., 5500.], dtype=float32)
array([ 0.5, 1.5, 2.5, ..., 357.5, 358.5, 359.5], dtype=float32)
array([-89.5, -88.5, -87.5, -86.5, -85.5, -84.5, -83.5, -82.5, -81.5, -80.5,\n", " -79.5, -78.5, -77.5, -76.5, -75.5, -74.5, -73.5, -72.5, -71.5, -70.5,\n", " -69.5, -68.5, -67.5, -66.5, -65.5, -64.5, -63.5, -62.5, -61.5, -60.5,\n", " -59.5, -58.5, -57.5, -56.5, -55.5, -54.5, -53.5, -52.5, -51.5, -50.5,\n", " -49.5, -48.5, -47.5, -46.5, -45.5, -44.5, -43.5, -42.5, -41.5, -40.5,\n", " -39.5, -38.5, -37.5, -36.5, -35.5, -34.5, -33.5, -32.5, -31.5, -30.5,\n", " -29.5, -28.5, -27.5, -26.5, -25.5, -24.5, -23.5, -22.5, -21.5, -20.5,\n", " -19.5, -18.5, -17.5, -16.5, -15.5, -14.5, -13.5, -12.5, -11.5, -10.5,\n", " -9.5, -8.5, -7.5, -6.5, -5.5, -4.5, -3.5, -2.5, -1.5, -0.5,\n", " 0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5,\n", " 10.5, 11.5, 12.5, 13.5, 14.5, 15.5, 16.5, 17.5, 18.5, 19.5,\n", " 20.5, 21.5, 22.5, 23.5, 24.5, 25.5, 26.5, 27.5, 28.5, 29.5,\n", " 30.5, 31.5, 32.5, 33.5, 34.5, 35.5, 36.5, 37.5, 38.5, 39.5,\n", " 40.5, 41.5, 42.5, 43.5, 44.5, 45.5, 46.5, 47.5, 48.5, 49.5,\n", " 50.5, 51.5, 52.5, 53.5, 54.5, 55.5, 56.5, 57.5, 58.5, 59.5,\n", " 60.5, 61.5, 62.5, 63.5, 64.5, 65.5, 66.5, 67.5, 68.5, 69.5,\n", " 70.5, 71.5, 72.5, 73.5, 74.5, 75.5, 76.5, 77.5, 78.5, 79.5,\n", " 80.5, 81.5, 82.5, 83.5, 84.5, 85.5, 86.5, 87.5, 88.5, 89.5],\n", " dtype=float32)
[2138400 values with dtype=float32]
<xarray.Dataset>\n", "Dimensions: (Z: 33, lat: 180, lon: 360)\n", "Coordinates:\n", " * Z (Z) float32 0.0 10.0 20.0 30.0 50.0 ... 4e+03 4.5e+03 5e+03 5.5e+03\n", " * lon (lon) float32 0.5 1.5 2.5 3.5 4.5 ... 355.5 356.5 357.5 358.5 359.5\n", " * lat (lat) float32 -89.5 -88.5 -87.5 -86.5 -85.5 ... 86.5 87.5 88.5 89.5\n", "Data variables:\n", " basin (Z, lat, lon) float32 ...\n", "Attributes:\n", " Conventions: IRIDL
array([ 0., 10., 20., 30., 50., 75., 100., 125., 150., 200.,\n", " 250., 300., 400., 500., 600., 700., 800., 900., 1000., 1100.,\n", " 1200., 1300., 1400., 1500., 1750., 2000., 2500., 3000., 3500., 4000.,\n", " 4500., 5000., 5500.], dtype=float32)
array([ 0.5, 1.5, 2.5, ..., 357.5, 358.5, 359.5], dtype=float32)
array([-89.5, -88.5, -87.5, -86.5, -85.5, -84.5, -83.5, -82.5, -81.5, -80.5,\n", " -79.5, -78.5, -77.5, -76.5, -75.5, -74.5, -73.5, -72.5, -71.5, -70.5,\n", " -69.5, -68.5, -67.5, -66.5, -65.5, -64.5, -63.5, -62.5, -61.5, -60.5,\n", " -59.5, -58.5, -57.5, -56.5, -55.5, -54.5, -53.5, -52.5, -51.5, -50.5,\n", " -49.5, -48.5, -47.5, -46.5, -45.5, -44.5, -43.5, -42.5, -41.5, -40.5,\n", " -39.5, -38.5, -37.5, -36.5, -35.5, -34.5, -33.5, -32.5, -31.5, -30.5,\n", " -29.5, -28.5, -27.5, -26.5, -25.5, -24.5, -23.5, -22.5, -21.5, -20.5,\n", " -19.5, -18.5, -17.5, -16.5, -15.5, -14.5, -13.5, -12.5, -11.5, -10.5,\n", " -9.5, -8.5, -7.5, -6.5, -5.5, -4.5, -3.5, -2.5, -1.5, -0.5,\n", " 0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5,\n", " 10.5, 11.5, 12.5, 13.5, 14.5, 15.5, 16.5, 17.5, 18.5, 19.5,\n", " 20.5, 21.5, 22.5, 23.5, 24.5, 25.5, 26.5, 27.5, 28.5, 29.5,\n", " 30.5, 31.5, 32.5, 33.5, 34.5, 35.5, 36.5, 37.5, 38.5, 39.5,\n", " 40.5, 41.5, 42.5, 43.5, 44.5, 45.5, 46.5, 47.5, 48.5, 49.5,\n", " 50.5, 51.5, 52.5, 53.5, 54.5, 55.5, 56.5, 57.5, 58.5, 59.5,\n", " 60.5, 61.5, 62.5, 63.5, 64.5, 65.5, 66.5, 67.5, 68.5, 69.5,\n", " 70.5, 71.5, 72.5, 73.5, 74.5, 75.5, 76.5, 77.5, 78.5, 79.5,\n", " 80.5, 81.5, 82.5, 83.5, 84.5, 85.5, 86.5, 87.5, 88.5, 89.5],\n", " dtype=float32)
[2138400 values with dtype=float32]
<xarray.DataArray 'basin' (lat: 180, lon: 360)>\n", "array([[nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " ...,\n", " [11., 11., 11., ..., 11., 11., 11.],\n", " [11., 11., 11., ..., 11., 11., 11.],\n", " [11., 11., 11., ..., 11., 11., 11.]], dtype=float32)\n", "Coordinates:\n", " Z float32 0.0\n", " * lon (lon) float32 0.5 1.5 2.5 3.5 4.5 ... 355.5 356.5 357.5 358.5 359.5\n", " * lat (lat) float32 -89.5 -88.5 -87.5 -86.5 -85.5 ... 86.5 87.5 88.5 89.5\n", "Attributes:\n", " long_name: basin code\n", " units: ids\n", " scale_max: 58\n", " CLIST: Atlantic Ocean\\nPacific Ocean \\nIndian Ocean\\nMediterranean S...\n", " valid_min: 1\n", " valid_max: 58\n", " scale_min: 1
array([[nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " [nan, nan, nan, ..., nan, nan, nan],\n", " ...,\n", " [11., 11., 11., ..., 11., 11., 11.],\n", " [11., 11., 11., ..., 11., 11., 11.],\n", " [11., 11., 11., ..., 11., 11., 11.]], dtype=float32)
array(0., dtype=float32)
array([ 0.5, 1.5, 2.5, ..., 357.5, 358.5, 359.5], dtype=float32)
array([-89.5, -88.5, -87.5, -86.5, -85.5, -84.5, -83.5, -82.5, -81.5, -80.5,\n", " -79.5, -78.5, -77.5, -76.5, -75.5, -74.5, -73.5, -72.5, -71.5, -70.5,\n", " -69.5, -68.5, -67.5, -66.5, -65.5, -64.5, -63.5, -62.5, -61.5, -60.5,\n", " -59.5, -58.5, -57.5, -56.5, -55.5, -54.5, -53.5, -52.5, -51.5, -50.5,\n", " -49.5, -48.5, -47.5, -46.5, -45.5, -44.5, -43.5, -42.5, -41.5, -40.5,\n", " -39.5, -38.5, -37.5, -36.5, -35.5, -34.5, -33.5, -32.5, -31.5, -30.5,\n", " -29.5, -28.5, -27.5, -26.5, -25.5, -24.5, -23.5, -22.5, -21.5, -20.5,\n", " -19.5, -18.5, -17.5, -16.5, -15.5, -14.5, -13.5, -12.5, -11.5, -10.5,\n", " -9.5, -8.5, -7.5, -6.5, -5.5, -4.5, -3.5, -2.5, -1.5, -0.5,\n", " 0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5,\n", " 10.5, 11.5, 12.5, 13.5, 14.5, 15.5, 16.5, 17.5, 18.5, 19.5,\n", " 20.5, 21.5, 22.5, 23.5, 24.5, 25.5, 26.5, 27.5, 28.5, 29.5,\n", " 30.5, 31.5, 32.5, 33.5, 34.5, 35.5, 36.5, 37.5, 38.5, 39.5,\n", " 40.5, 41.5, 42.5, 43.5, 44.5, 45.5, 46.5, 47.5, 48.5, 49.5,\n", " 50.5, 51.5, 52.5, 53.5, 54.5, 55.5, 56.5, 57.5, 58.5, 59.5,\n", " 60.5, 61.5, 62.5, 63.5, 64.5, 65.5, 66.5, 67.5, 68.5, 69.5,\n", " 70.5, 71.5, 72.5, 73.5, 74.5, 75.5, 76.5, 77.5, 78.5, 79.5,\n", " 80.5, 81.5, 82.5, 83.5, 84.5, 85.5, 86.5, 87.5, 88.5, 89.5],\n", " dtype=float32)
<xarray.DataArray 'sst' (time: 708, basin: 14)>\n", "array([[-1.8 , -1.8 , 23.455315 , ..., -1.8 ,\n", " 3.3971915 , 24.182198 ],\n", " [-1.8 , -1.8 , 23.722523 , ..., -1.8 ,\n", " 0.03573781, 24.59657 ],\n", " [-1.8 , -1.8 , 24.601315 , ..., -1.8 ,\n", " -0.26487017, 26.234186 ],\n", " ...,\n", " [ 0.6758096 , 6.5041866 , 29.279472 , ..., 10.920225 ,\n", " 15.955018 , 29.419764 ],\n", " [-0.7937473 , 3.0715103 , 27.608444 , ..., 5.4078865 ,\n", " 10.673691 , 27.755806 ],\n", " [-1.8 , -0.06061229, 25.881489 , ..., 0.5253577 ,\n", " 7.2676916 , 26.16315 ]], dtype=float32)\n", "Coordinates:\n", " * time (time) datetime64[ns] 1960-01-01 1960-02-01 ... 2018-12-01\n", " Z float32 0.0\n", " * basin (basin) float64 1.0 2.0 3.0 4.0 5.0 ... 10.0 11.0 12.0 53.0 56.0\n", "Attributes:\n", " long_name: Monthly Means of Sea Surface Temperature\n", " units: degC\n", " var_desc: Sea Surface Temperature\n", " level_desc: Surface\n", " statistic: Mean\n", " dataset: NOAA Extended Reconstructed SST V5\n", " parent_stat: Individual Values\n", " actual_range: [-1.8 42.32636]\n", " valid_range: [-1.8 45. ]\n", " _ChunkSizes: [ 1 89 180]
array([[-1.8 , -1.8 , 23.455315 , ..., -1.8 ,\n", " 3.3971915 , 24.182198 ],\n", " [-1.8 , -1.8 , 23.722523 , ..., -1.8 ,\n", " 0.03573781, 24.59657 ],\n", " [-1.8 , -1.8 , 24.601315 , ..., -1.8 ,\n", " -0.26487017, 26.234186 ],\n", " ...,\n", " [ 0.6758096 , 6.5041866 , 29.279472 , ..., 10.920225 ,\n", " 15.955018 , 29.419764 ],\n", " [-0.7937473 , 3.0715103 , 27.608444 , ..., 5.4078865 ,\n", " 10.673691 , 27.755806 ],\n", " [-1.8 , -0.06061229, 25.881489 , ..., 0.5253577 ,\n", " 7.2676916 , 26.16315 ]], dtype=float32)
array(['1960-01-01T00:00:00.000000000', '1960-02-01T00:00:00.000000000',\n", " '1960-03-01T00:00:00.000000000', ..., '2018-10-01T00:00:00.000000000',\n", " '2018-11-01T00:00:00.000000000', '2018-12-01T00:00:00.000000000'],\n", " dtype='datetime64[ns]')
array(0., dtype=float32)
array([ 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11., 12., 53., 56.])
<xarray.DataArray 'sst' (month: 12, basin: 14)>\n", "array([[18.931244 , 20.88606 , 22.068876 , 15.7806225 , 2.9009485 ,\n", " 8.286119 , 25.862597 , 21.052422 , -1.7988946 , 2.5364153 ,\n", " -1.2113837 , 5.99166 , 6.69535 , 26.615591 ],\n", " [18.998413 , 20.938007 , 22.466644 , 15.171556 , 1.6216687 ,\n", " 7.475156 , 25.31563 , 20.009794 , -1.8000002 , 2.86857 ,\n", " -1.2622244 , 4.981395 , 5.473524 , 26.889442 ],\n", " [19.013428 , 20.972063 , 22.608067 , 15.258743 , 1.4993656 ,\n", " 7.793447 , 25.649214 , 20.699417 , -1.8000002 , 2.6315544 ,\n", " -1.2847868 , 5.1045027 , 5.522258 , 28.037973 ],\n", " [18.989248 , 20.968798 , 22.381023 , 16.224312 , 3.0491576 ,\n", " 10.017782 , 26.88303 , 23.203306 , -1.7996598 , 2.094409 ,\n", " -1.2763202 , 6.663208 , 9.4189415 , 29.374727 ],\n", " [19.068415 , 21.018923 , 21.714952 , 18.613165 , 6.4635086 ,\n", " 14.915698 , 28.640142 , 26.750372 , -1.7467502 , 1.571553 ,\n", " -1.1642531 , 9.768522 , 15.095652 , 29.972607 ],\n", " [19.352161 , 21.17309 , 20.790792 , 21.979132 , 11.794082 ,\n", " 20.411827 , 29.786642 , 29.794815 , -0.96558964, 1.1652179 ,\n", " -0.8448692 , 14.277155 , 20.739882 , 29.463205 ],\n", " [19.71806 , 21.479403 , 20.011818 , 24.690279 , 16.125217 ,\n", " 23.58809 , 30.706587 , 31.815105 , 2.8862731 , 0.8692938 ,\n", " -0.30307108, 18.744436 , 24.158962 , 28.86987 ],\n", " [19.959143 , 21.723642 , 19.603079 , 25.818441 , 17.186897 ,\n", " 24.277603 , 31.248007 , 32.692772 , 5.3934736 , 0.68534446,\n", " 0.16548847, 21.910128 , 24.577526 , 28.628706 ],\n", " [19.816757 , 21.665964 , 19.646137 , 24.763577 , 14.3625765 ,\n", " 21.748938 , 31.174229 , 32.241535 , 4.6783195 , 0.6313895 ,\n", " 0.03611656, 20.365383 , 21.501966 , 28.853634 ],\n", " [19.467312 , 21.347132 , 20.05091 , 22.60426 , 10.456942 ,\n", " 17.874405 , 30.591843 , 30.352285 , 2.3623 , 0.72350246,\n", " -0.554342 , 16.454807 , 16.766832 , 29.040428 ],\n", " [19.136904 , 21.059431 , 20.721912 , 19.869085 , 7.2571416 ,\n", " 13.941229 , 28.884726 , 27.106995 , 0.03359437, 1.0484893 ,\n", " -0.97754955, 12.12734 , 12.716839 , 28.491669 ],\n", " [18.968773 , 20.906147 , 21.460379 , 17.377382 , 4.863486 ,\n", " 10.682313 , 27.186651 , 23.71751 , -1.7128265 , 1.7405492 ,\n", " -1.1221472 , 8.642129 , 9.399436 , 27.350985 ]],\n", " dtype=float32)\n", "Coordinates:\n", " * month (month) int64 1 2 3 4 5 6 7 8 9 10 11 12\n", " Z float32 0.0\n", " * basin (basin) float64 1.0 2.0 3.0 4.0 5.0 ... 10.0 11.0 12.0 53.0 56.0
array([[18.931244 , 20.88606 , 22.068876 , 15.7806225 , 2.9009485 ,\n", " 8.286119 , 25.862597 , 21.052422 , -1.7988946 , 2.5364153 ,\n", " -1.2113837 , 5.99166 , 6.69535 , 26.615591 ],\n", " [18.998413 , 20.938007 , 22.466644 , 15.171556 , 1.6216687 ,\n", " 7.475156 , 25.31563 , 20.009794 , -1.8000002 , 2.86857 ,\n", " -1.2622244 , 4.981395 , 5.473524 , 26.889442 ],\n", " [19.013428 , 20.972063 , 22.608067 , 15.258743 , 1.4993656 ,\n", " 7.793447 , 25.649214 , 20.699417 , -1.8000002 , 2.6315544 ,\n", " -1.2847868 , 5.1045027 , 5.522258 , 28.037973 ],\n", " [18.989248 , 20.968798 , 22.381023 , 16.224312 , 3.0491576 ,\n", " 10.017782 , 26.88303 , 23.203306 , -1.7996598 , 2.094409 ,\n", " -1.2763202 , 6.663208 , 9.4189415 , 29.374727 ],\n", " [19.068415 , 21.018923 , 21.714952 , 18.613165 , 6.4635086 ,\n", " 14.915698 , 28.640142 , 26.750372 , -1.7467502 , 1.571553 ,\n", " -1.1642531 , 9.768522 , 15.095652 , 29.972607 ],\n", " [19.352161 , 21.17309 , 20.790792 , 21.979132 , 11.794082 ,\n", " 20.411827 , 29.786642 , 29.794815 , -0.96558964, 1.1652179 ,\n", " -0.8448692 , 14.277155 , 20.739882 , 29.463205 ],\n", " [19.71806 , 21.479403 , 20.011818 , 24.690279 , 16.125217 ,\n", " 23.58809 , 30.706587 , 31.815105 , 2.8862731 , 0.8692938 ,\n", " -0.30307108, 18.744436 , 24.158962 , 28.86987 ],\n", " [19.959143 , 21.723642 , 19.603079 , 25.818441 , 17.186897 ,\n", " 24.277603 , 31.248007 , 32.692772 , 5.3934736 , 0.68534446,\n", " 0.16548847, 21.910128 , 24.577526 , 28.628706 ],\n", " [19.816757 , 21.665964 , 19.646137 , 24.763577 , 14.3625765 ,\n", " 21.748938 , 31.174229 , 32.241535 , 4.6783195 , 0.6313895 ,\n", " 0.03611656, 20.365383 , 21.501966 , 28.853634 ],\n", " [19.467312 , 21.347132 , 20.05091 , 22.60426 , 10.456942 ,\n", " 17.874405 , 30.591843 , 30.352285 , 2.3623 , 0.72350246,\n", " -0.554342 , 16.454807 , 16.766832 , 29.040428 ],\n", " [19.136904 , 21.059431 , 20.721912 , 19.869085 , 7.2571416 ,\n", " 13.941229 , 28.884726 , 27.106995 , 0.03359437, 1.0484893 ,\n", " -0.97754955, 12.12734 , 12.716839 , 28.491669 ],\n", " [18.968773 , 20.906147 , 21.460379 , 17.377382 , 4.863486 ,\n", " 10.682313 , 27.186651 , 23.71751 , -1.7128265 , 1.7405492 ,\n", " -1.1221472 , 8.642129 , 9.399436 , 27.350985 ]],\n", " dtype=float32)
array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12])
array(0., dtype=float32)
array([ 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11., 12., 53., 56.])