{
"cells": [
{
"cell_type": "markdown",
"metadata": {
"colab_type": "text",
"id": "JkUviKhSYLMF",
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"## Climatology and Anomaly in Geoscience\n",
"An example of how Python is used to explore global sea surface temperature (SST).\n",
"* Wenchang Yang (wenchang@princeton.edu)\n",
"* Department of Geosciences, Princeton University\n",
"* Junior Colloquium, Nov. 09, 2020\n"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"## What's covered so far\n",
"1. Python basics: `number`, `string`, `list`, `function`, `module`, `package`, ...\n",
"2. Scientific computation: `numpy`, `scipy`\n",
"3. Data visualization: `matplotlib`\n",
"4. High-level and user-friendly packages: `pandas` and `xarray`"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"## Today's plan\n",
"\n",
"1. A short introduction to `xarray`.\n",
"2. SST dataset.\n",
"3. Scientific questions to keep in mind.\n",
"4. Step-by-step data analysis using Python.\n"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"## What's `xarray` able to do? \n",
"\n",
"
\n",
"\n",
"1. Open/save datasets (single/multiple, local/remote): `open_dataset`, `open_mfdataset`.\n",
"2. Data selection: `sel`, `isel`.\n",
"3. Computation: `mean`, `std`, `max`, `min`, `differentiate`, `integrate`, ...\n",
"4. Split-apply-combine: `groupby`.\n",
"5. Plot: `plot`, `plot.line`, `plot.contourf`, `plot.hist` ...\n",
"\n",
"http://xarray.pydata.org/"
]
},
{
"cell_type": "markdown",
"metadata": {
"colab_type": "text",
"id": "TOPbnilz9XD5",
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"## SST data\n",
"* ERSST version 5: global **monthly** SST.\n",
"* $2^\\circ$ longitude $\\times$ $2^\\circ$ latitude \n",
"* From Columbia University [data library](http://iridl.ldeo.columbia.edu): http://iridl.ldeo.columbia.edu/SOURCES/.NOAA/.NCDC/.ERSST/.version5/.sst/T/%28Jan%201979%29/%28Dec%202018%29/RANGE/X//lon/renameGRID/Y//lat/renameGRID/T//time/renameGRID/time/(days%20since%201979-01-01)/streamgridunitconvert%5Bzlev%5Daverage\n",
"* It covers 1854-present, but we focus on **1979-2018** today.\n",
"* Downloaded and available on the Adroit server:
`/home/wenchang/JC2020/ersst5_1979-2018.nc`"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"## Some scientific questions to keep in mind\n",
"1. What does global SST pattern look like?\n",
"2. Is SST getting warmer over the recent decades?\n",
"2. How do El Nino/La Nina vary during this period?\n",
" \n",
"\n",
"https://blog.weatherops.com/hubfs/blog-files/elnino-vs-lanina-noaa.jpg"
]
},
{
"cell_type": "markdown",
"metadata": {
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"## Start analysis"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"ExecuteTime": {
"end_time": "2019-10-21T16:43:26.868932Z",
"start_time": "2019-10-21T16:43:25.222395Z"
},
"colab": {},
"colab_type": "code",
"id": "2nyvaFZnYYtm",
"slideshow": {
"slide_type": "fragment"
}
},
"outputs": [],
"source": [
"# xarray is the core package we are going to use\n",
"import xarray as xr\n",
"import matplotlib.pyplot as plt # we also use pyplot directly in some cases\n",
"import os"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"ExecuteTime": {
"end_time": "2019-10-21T16:44:03.579148Z",
"start_time": "2019-10-21T16:44:03.564588Z"
},
"slideshow": {
"slide_type": "fragment"
}
},
"outputs": [],
"source": [
"# some configurations on the default figure output\n",
"%config InlineBackend.figure_format ='retina'\n",
"plt.rcParams['figure.dpi'] = 120"
]
},
{
"cell_type": "markdown",
"metadata": {
"colab_type": "text",
"id": "azoylog2-GzB",
"slideshow": {
"slide_type": "slide"
}
},
"source": [
"## Open the SST data file\n",
"Use `xr.open_dataset` to open the SST data file.\n"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"ExecuteTime": {
"end_time": "2019-10-21T16:44:56.636403Z",
"start_time": "2019-10-21T16:44:53.628388Z"
},
"colab": {
"base_uri": "https://localhost:8080/",
"height": 244
},
"colab_type": "code",
"executionInfo": {
"elapsed": 835,
"status": "ok",
"timestamp": 1568235456809,
"user": {
"displayName": "Wenchang Yang",
"photoUrl": "https://lh3.googleusercontent.com/a-/AAuE7mAzexWcchBOvHGw9u_Nm2D-vWc4ApTqQ4uLX1i-=s64",
"userId": "02317458745209383076"
},
"user_tz": 240
},
"id": "bik3mA4tYefp",
"outputId": "c95cc7cf-0866-43b3-ca86-44fa0d0320f8",
"slideshow": {
"slide_type": "fragment"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"file to be opened: ersst5_1979-2018.nc\n"
]
}
],
"source": [
"# Am I running the notebook on Adroit or not?\n",
"if os.uname().nodename.startswith('adroit'):\n",
" ifile = '/home/wenchang/JC2020/ersst5_1979-2018.nc'\n",
"else:\n",
" ifile = 'ersst5_1979-2018.nc'\n",
"print('file to be opened:',ifile)\n",
"ds = xr.open_dataset(ifile)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"slideshow": {
"slide_type": "subslide"
}
},
"outputs": [
{
"data": {
"text/html": [
"
<xarray.Dataset>\n",
"Dimensions: (lat: 89, lon: 180, time: 480)\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] 1979-01-16T12:00:00 ... 2018-12-16T12:00:00\n",
"Data variables:\n",
" sst (time, lat, lon) float32 ...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(['1979-01-16T12:00:00.000000000', '1979-02-15T00:00:00.000000000',\n",
" '1979-03-16T12:00:00.000000000', ..., '2018-10-16T12:00:00.000000000',\n",
" '2018-11-16T00:00:00.000000000', '2018-12-16T12:00:00.000000000'],\n",
" dtype='datetime64[ns]')[7689600 values with dtype=float32]
<xarray.DataArray 'sst' (time: 480, lat: 89, lon: 180)>\n",
"[7689600 values with dtype=float32]\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] 1979-01-16T12:00:00 ... 2018-12-16T12:00:00\n",
"Attributes:\n",
" pointwidth: 1.0\n",
" long_name: Extended reconstructed sea surface temperature\n",
" units: Celsius_scale\n",
" standard_name: sea_surface_temperature\n",
" expires: 1604275200\n",
" file_missing_value: -999.0\n",
" history: Averaged over zlev=0.0 meters minimum 0.0% data present[7689600 values with dtype=float32]
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(['1979-01-16T12:00:00.000000000', '1979-02-15T00:00:00.000000000',\n",
" '1979-03-16T12:00:00.000000000', ..., '2018-10-16T12:00:00.000000000',\n",
" '2018-11-16T00:00:00.000000000', '2018-12-16T12:00:00.000000000'],\n",
" dtype='datetime64[ns]')<xarray.DataArray 'sst' (lat: 89, lon: 180)>\n",
"array([[ nan, nan, nan, ..., nan, nan,\n",
" nan],\n",
" [ nan, nan, nan, ..., nan, nan,\n",
" nan],\n",
" [ nan, nan, nan, ..., nan, nan,\n",
" nan],\n",
" ...,\n",
" [-1.7999686, -1.799723 , -1.7995721, ..., -1.799789 , -1.7998605,\n",
" -1.7999344],\n",
" [-1.7991427, -1.7990254, -1.7989469, ..., -1.7997802, -1.7994727,\n",
" -1.7992976],\n",
" [-1.7990803, -1.7990966, -1.7991083, ..., -1.7993528, -1.7991874,\n",
" -1.7990723]], dtype=float32)\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.0array([[ nan, nan, nan, ..., nan, nan,\n",
" nan],\n",
" [ nan, nan, nan, ..., nan, nan,\n",
" nan],\n",
" [ nan, nan, nan, ..., nan, nan,\n",
" nan],\n",
" ...,\n",
" [-1.7999686, -1.799723 , -1.7995721, ..., -1.799789 , -1.7998605,\n",
" -1.7999344],\n",
" [-1.7991427, -1.7990254, -1.7989469, ..., -1.7997802, -1.7994727,\n",
" -1.7992976],\n",
" [-1.7990803, -1.7990966, -1.7991083, ..., -1.7993528, -1.7991874,\n",
" -1.7990723]], dtype=float32)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)<xarray.DataArray 'sst' (lat: 89, lon: 180)>\n",
"array([[ 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",
" [9.4175339e-05, 1.0762215e-03, 1.6802549e-03, ..., 7.9834461e-04,\n",
" 5.2642822e-04, 2.3126602e-04],\n",
" [3.3977032e-03, 3.8667917e-03, 4.1812658e-03, ..., 8.4733963e-04,\n",
" 2.0782948e-03, 2.7782917e-03],\n",
" [3.6476851e-03, 3.5823584e-03, 3.5355091e-03, ..., 2.5581121e-03,\n",
" 3.2191277e-03, 3.6793947e-03]], dtype=float32)\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",
"Attributes:\n",
" long_name: SST change from 1979-1988 to 2009-2018\n",
" units: $^\\circ$Carray([[ 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",
" [9.4175339e-05, 1.0762215e-03, 1.6802549e-03, ..., 7.9834461e-04,\n",
" 5.2642822e-04, 2.3126602e-04],\n",
" [3.3977032e-03, 3.8667917e-03, 4.1812658e-03, ..., 8.4733963e-04,\n",
" 2.0782948e-03, 2.7782917e-03],\n",
" [3.6476851e-03, 3.5823584e-03, 3.5355091e-03, ..., 2.5581121e-03,\n",
" 3.2191277e-03, 3.6793947e-03]], dtype=float32)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)<xarray.DataArray 'sst' (month: 12, lat: 89, lon: 180)>\n",
"array([[[ 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.7999998, -1.7999998, -1.7999998, ..., -1.7999998,\n",
" -1.7999998, -1.7999998],\n",
" [-1.7999998, -1.7999998, -1.7999998, ..., -1.7999998,\n",
" -1.7999998, -1.7999998],\n",
" [-1.7999998, -1.7999998, -1.7999998, ..., -1.7999998,\n",
" -1.7999998, -1.7999998]],\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.7999998, -1.7967713, -1.7949593, ..., -1.7999998,\n",
" -1.7999998, -1.7999998],\n",
" [-1.7904856, -1.7887875, -1.7874559, ..., -1.7974577,\n",
" -1.7940279, -1.7922757],\n",
" [-1.7907732, -1.7902114, -1.7898309, ..., -1.794199 ,\n",
" -1.7924855, -1.7912834]],\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.7999998, -1.7999998, -1.7999998, ..., -1.7999998,\n",
" -1.7999998, -1.7999998],\n",
" [-1.7993208, -1.7996114, -1.7999998, ..., -1.7999998,\n",
" -1.7997364, -1.7993886],\n",
" [-1.7992064, -1.7993428, -1.7997385, ..., -1.7996912,\n",
" -1.799461 , -1.7992189]]], dtype=float32)\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 12array([[[ 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.7999998, -1.7999998, -1.7999998, ..., -1.7999998,\n",
" -1.7999998, -1.7999998],\n",
" [-1.7999998, -1.7999998, -1.7999998, ..., -1.7999998,\n",
" -1.7999998, -1.7999998],\n",
" [-1.7999998, -1.7999998, -1.7999998, ..., -1.7999998,\n",
" -1.7999998, -1.7999998]],\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.7999998, -1.7967713, -1.7949593, ..., -1.7999998,\n",
" -1.7999998, -1.7999998],\n",
" [-1.7904856, -1.7887875, -1.7874559, ..., -1.7974577,\n",
" -1.7940279, -1.7922757],\n",
" [-1.7907732, -1.7902114, -1.7898309, ..., -1.794199 ,\n",
" -1.7924855, -1.7912834]],\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.7999998, -1.7999998, -1.7999998, ..., -1.7999998,\n",
" -1.7999998, -1.7999998],\n",
" [-1.7993208, -1.7996114, -1.7999998, ..., -1.7999998,\n",
" -1.7997364, -1.7993886],\n",
" [-1.7992064, -1.7993428, -1.7997385, ..., -1.7996912,\n",
" -1.799461 , -1.7992189]]], dtype=float32)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])
<xarray.DataArray 'sst' (time: 480, lat: 89, lon: 180)>\n",
"array([[[ 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",
" [-1.1920929e-07, -1.1920929e-07, -1.1920929e-07, ...,\n",
" -1.1920929e-07, -1.1920929e-07, -1.1920929e-07],\n",
" [-1.1920929e-07, -1.1920929e-07, -1.1920929e-07, ...,\n",
" -1.1920929e-07, -1.1920929e-07, -1.1920929e-07],\n",
" [-1.1920929e-07, -1.1920929e-07, -1.1920929e-07, ...,\n",
" -1.1920929e-07, -1.1920929e-07, -1.1920929e-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",
" [-1.1920929e-07, -3.2286644e-03, -5.0406456e-03, ...,\n",
" -1.1920929e-07, -1.1920929e-07, -1.1920929e-07],\n",
" [-9.5143318e-03, -1.1212468e-02, -1.2544036e-02, ...,\n",
" -2.5422573e-03, -5.9720278e-03, -7.7242851e-03],\n",
" [-9.2267990e-03, -9.7885132e-03, -1.0169029e-02, ...,\n",
" -5.8009624e-03, -7.5144768e-03, -8.7165833e-03]],\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",
" [-1.1920929e-07, -1.1920929e-07, -1.1920929e-07, ...,\n",
" -1.1920929e-07, -1.1920929e-07, -1.1920929e-07],\n",
" [-6.7913532e-04, -3.8850307e-04, -1.1920929e-07, ...,\n",
" -1.1920929e-07, -2.6357174e-04, -6.1130524e-04],\n",
" [-7.9357624e-04, -6.5720081e-04, -2.6142597e-04, ...,\n",
" -3.0875206e-04, -5.3894520e-04, -7.8105927e-04]]], dtype=float32)\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] 1979-01-16T12:00:00 ... 2018-12-16T12:00:00\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 12array([[[ 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",
" [-1.1920929e-07, -1.1920929e-07, -1.1920929e-07, ...,\n",
" -1.1920929e-07, -1.1920929e-07, -1.1920929e-07],\n",
" [-1.1920929e-07, -1.1920929e-07, -1.1920929e-07, ...,\n",
" -1.1920929e-07, -1.1920929e-07, -1.1920929e-07],\n",
" [-1.1920929e-07, -1.1920929e-07, -1.1920929e-07, ...,\n",
" -1.1920929e-07, -1.1920929e-07, -1.1920929e-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",
" [-1.1920929e-07, -3.2286644e-03, -5.0406456e-03, ...,\n",
" -1.1920929e-07, -1.1920929e-07, -1.1920929e-07],\n",
" [-9.5143318e-03, -1.1212468e-02, -1.2544036e-02, ...,\n",
" -2.5422573e-03, -5.9720278e-03, -7.7242851e-03],\n",
" [-9.2267990e-03, -9.7885132e-03, -1.0169029e-02, ...,\n",
" -5.8009624e-03, -7.5144768e-03, -8.7165833e-03]],\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",
" [-1.1920929e-07, -1.1920929e-07, -1.1920929e-07, ...,\n",
" -1.1920929e-07, -1.1920929e-07, -1.1920929e-07],\n",
" [-6.7913532e-04, -3.8850307e-04, -1.1920929e-07, ...,\n",
" -1.1920929e-07, -2.6357174e-04, -6.1130524e-04],\n",
" [-7.9357624e-04, -6.5720081e-04, -2.6142597e-04, ...,\n",
" -3.0875206e-04, -5.3894520e-04, -7.8105927e-04]]], dtype=float32)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(['1979-01-16T12:00:00.000000000', '1979-02-15T00:00:00.000000000',\n",
" '1979-03-16T12:00:00.000000000', ..., '2018-10-16T12:00:00.000000000',\n",
" '2018-11-16T00:00:00.000000000', '2018-12-16T12: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",
" 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])<xarray.DataArray 'sst' (time: 480)>\n",
"array([-1.04900680e-01, -1.82981327e-01, 3.79271880e-02, 1.02511741e-01,\n",
" -1.56986818e-01, -2.13878661e-01, -3.84980589e-01, -9.16783959e-02,\n",
" 2.33983338e-01, 1.16299964e-01, 3.04566145e-01, 4.96357203e-01,\n",
" 5.66177428e-01, 2.53595203e-01, 7.83962235e-02, 2.09114216e-02,\n",
" 1.65222645e-01, 2.91420966e-01, 2.07556952e-02, -3.65579039e-01,\n",
" -3.15712154e-01, -2.53466249e-01, -8.54855329e-02, 2.84164567e-02,\n",
" -3.33171606e-01, -6.02614760e-01, -5.78435421e-01, -4.41379726e-01,\n",
" -4.88844067e-01, -3.81764174e-01, -5.59118807e-01, -5.22406876e-01,\n",
" -2.35358313e-01, -2.97763467e-01, -3.94616067e-01, -2.15333968e-01,\n",
" 1.56619325e-01, -1.26914948e-01, 1.85399741e-01, 2.95923799e-01,\n",
" 5.35673141e-01, 6.08795702e-01, 4.48002100e-01, 7.32922375e-01,\n",
" 1.44802785e+00, 1.97870481e+00, 1.92627609e+00, 2.17196774e+00,\n",
" 2.36158586e+00, 1.96175349e+00, 1.42730427e+00, 1.03574896e+00,\n",
" 9.88953352e-01, 6.29534662e-01, -1.97548345e-02, -2.45766759e-01,\n",
" -5.52735507e-01, -1.14381588e+00, -1.26753747e+00, -1.04204881e+00,\n",
" -6.42956257e-01, -1.56747371e-01, -4.61359173e-01, -5.92787027e-01,\n",
" -6.51417553e-01, -8.22249234e-01, -4.28176045e-01, -2.39528924e-01,\n",
" -3.80766124e-01, -6.75847292e-01, -1.16851592e+00, -1.35727656e+00,\n",
" -1.12902832e+00, -6.89546764e-01, -7.34117508e-01, -1.08534360e+00,\n",
" -9.41993594e-01, -8.35659027e-01, -6.60370827e-01, -5.49420297e-01,\n",
"...\n",
" 3.38472098e-01, 2.61887074e-01, 1.68515474e-01, -2.75554210e-01,\n",
" -4.97084290e-01, -4.83319312e-01, -2.04805434e-01, -1.70178846e-01,\n",
" -3.27110559e-01, -3.59871030e-01, -3.20950657e-01, -3.80554438e-01,\n",
" -2.30542362e-01, -2.09253475e-01, -1.38686135e-01, -2.06214979e-01,\n",
" -4.76016909e-01, -5.97726822e-01, -2.38656268e-01, 1.58438459e-01,\n",
" 3.86419833e-01, 2.90219992e-01, 8.36944557e-04, -3.28716263e-02,\n",
" 2.31316730e-01, 4.24140334e-01, 7.49024451e-01, 6.81629837e-01,\n",
" 5.40025949e-01, 4.66240525e-01, 5.33091247e-01, 7.78923273e-01,\n",
" 9.95902777e-01, 1.24814522e+00, 1.52608073e+00, 1.92614651e+00,\n",
" 2.15197659e+00, 2.35080457e+00, 2.71775723e+00, 2.63437271e+00,\n",
" 2.58447456e+00, 2.27873635e+00, 1.67325842e+00, 9.90715742e-01,\n",
" 3.94458622e-01, 5.84407821e-02, -3.88212562e-01, -5.55807173e-01,\n",
" -6.04830086e-01, -7.31975555e-01, -7.34857500e-01, -5.19446015e-01,\n",
" -3.96952122e-01, -5.30423373e-02, 8.23768601e-02, 2.86082298e-01,\n",
" 4.33291554e-01, 4.06300545e-01, 3.33628446e-01, -1.46758303e-01,\n",
" -4.73812044e-01, -5.67844987e-01, -9.50382292e-01, -1.01649618e+00,\n",
" -9.64271903e-01, -7.61266887e-01, -7.58066595e-01, -4.35995132e-01,\n",
" -1.26900762e-01, 1.22105628e-01, 2.05426410e-01, 1.05248615e-01,\n",
" 4.28478599e-01, 8.99333298e-01, 9.17508781e-01, 8.60169709e-01],\n",
" dtype=float32)\n",
"Coordinates:\n",
" * time (time) datetime64[ns] 1979-01-16T12:00:00 ... 2018-12-16T12:00:00\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",
"Attributes:\n",
" long_name: Nino3.4 indexarray([-1.04900680e-01, -1.82981327e-01, 3.79271880e-02, 1.02511741e-01,\n",
" -1.56986818e-01, -2.13878661e-01, -3.84980589e-01, -9.16783959e-02,\n",
" 2.33983338e-01, 1.16299964e-01, 3.04566145e-01, 4.96357203e-01,\n",
" 5.66177428e-01, 2.53595203e-01, 7.83962235e-02, 2.09114216e-02,\n",
" 1.65222645e-01, 2.91420966e-01, 2.07556952e-02, -3.65579039e-01,\n",
" -3.15712154e-01, -2.53466249e-01, -8.54855329e-02, 2.84164567e-02,\n",
" -3.33171606e-01, -6.02614760e-01, -5.78435421e-01, -4.41379726e-01,\n",
" -4.88844067e-01, -3.81764174e-01, -5.59118807e-01, -5.22406876e-01,\n",
" -2.35358313e-01, -2.97763467e-01, -3.94616067e-01, -2.15333968e-01,\n",
" 1.56619325e-01, -1.26914948e-01, 1.85399741e-01, 2.95923799e-01,\n",
" 5.35673141e-01, 6.08795702e-01, 4.48002100e-01, 7.32922375e-01,\n",
" 1.44802785e+00, 1.97870481e+00, 1.92627609e+00, 2.17196774e+00,\n",
" 2.36158586e+00, 1.96175349e+00, 1.42730427e+00, 1.03574896e+00,\n",
" 9.88953352e-01, 6.29534662e-01, -1.97548345e-02, -2.45766759e-01,\n",
" -5.52735507e-01, -1.14381588e+00, -1.26753747e+00, -1.04204881e+00,\n",
" -6.42956257e-01, -1.56747371e-01, -4.61359173e-01, -5.92787027e-01,\n",
" -6.51417553e-01, -8.22249234e-01, -4.28176045e-01, -2.39528924e-01,\n",
" -3.80766124e-01, -6.75847292e-01, -1.16851592e+00, -1.35727656e+00,\n",
" -1.12902832e+00, -6.89546764e-01, -7.34117508e-01, -1.08534360e+00,\n",
" -9.41993594e-01, -8.35659027e-01, -6.60370827e-01, -5.49420297e-01,\n",
"...\n",
" 3.38472098e-01, 2.61887074e-01, 1.68515474e-01, -2.75554210e-01,\n",
" -4.97084290e-01, -4.83319312e-01, -2.04805434e-01, -1.70178846e-01,\n",
" -3.27110559e-01, -3.59871030e-01, -3.20950657e-01, -3.80554438e-01,\n",
" -2.30542362e-01, -2.09253475e-01, -1.38686135e-01, -2.06214979e-01,\n",
" -4.76016909e-01, -5.97726822e-01, -2.38656268e-01, 1.58438459e-01,\n",
" 3.86419833e-01, 2.90219992e-01, 8.36944557e-04, -3.28716263e-02,\n",
" 2.31316730e-01, 4.24140334e-01, 7.49024451e-01, 6.81629837e-01,\n",
" 5.40025949e-01, 4.66240525e-01, 5.33091247e-01, 7.78923273e-01,\n",
" 9.95902777e-01, 1.24814522e+00, 1.52608073e+00, 1.92614651e+00,\n",
" 2.15197659e+00, 2.35080457e+00, 2.71775723e+00, 2.63437271e+00,\n",
" 2.58447456e+00, 2.27873635e+00, 1.67325842e+00, 9.90715742e-01,\n",
" 3.94458622e-01, 5.84407821e-02, -3.88212562e-01, -5.55807173e-01,\n",
" -6.04830086e-01, -7.31975555e-01, -7.34857500e-01, -5.19446015e-01,\n",
" -3.96952122e-01, -5.30423373e-02, 8.23768601e-02, 2.86082298e-01,\n",
" 4.33291554e-01, 4.06300545e-01, 3.33628446e-01, -1.46758303e-01,\n",
" -4.73812044e-01, -5.67844987e-01, -9.50382292e-01, -1.01649618e+00,\n",
" -9.64271903e-01, -7.61266887e-01, -7.58066595e-01, -4.35995132e-01,\n",
" -1.26900762e-01, 1.22105628e-01, 2.05426410e-01, 1.05248615e-01,\n",
" 4.28478599e-01, 8.99333298e-01, 9.17508781e-01, 8.60169709e-01],\n",
" dtype=float32)array(['1979-01-16T12:00:00.000000000', '1979-02-15T00:00:00.000000000',\n",
" '1979-03-16T12:00:00.000000000', ..., '2018-10-16T12:00:00.000000000',\n",
" '2018-11-16T00:00:00.000000000', '2018-12-16T12: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",
" 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])