{ "cells": [ { "cell_type": "markdown", "id": "0702ec6c-f370-4cad-bb32-bb80d4335bdb", "metadata": {}, "source": [ "# Event" ] }, { "cell_type": "code", "execution_count": 1, "id": "4ce4b199-5572-4e33-83b6-593ce09e57a4", "metadata": {}, "outputs": [], "source": [ "import xarray\n", "import numpy\n", "import pandas\n", "import climtas\n", "import xesmf\n", "import dask.array" ] }, { "cell_type": "markdown", "id": "bfc85e52-2dc0-44f7-b271-8e59361fbc6e", "metadata": {}, "source": [ "We have a Dask dataset, and we'd like to identify periods where the value is above some threshold" ] }, { "cell_type": "code", "execution_count": 2, "id": "7c1dcb98-fbaf-403e-a388-523b73fc6821", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.DataArray 'temperature' (time: 1095, lat: 50, lon: 100)>\n", "dask.array<random_sample, shape=(1095, 50, 100), dtype=float64, chunksize=(90, 25, 25), chunktype=numpy.ndarray>\n", "Coordinates:\n", " * time (time) datetime64[ns] 2001-01-01 2001-01-02 ... 2003-12-31\n", " * lat (lat) float64 -90.0 -86.33 -82.65 -78.98 ... 78.98 82.65 86.33 90.0\n", " * lon (lon) float64 -180.0 -176.4 -172.8 -169.2 ... 169.2 172.8 176.4
\n",
"
| \n",
"\n", "\n", " | \n", "
array(['2001-01-01T00:00:00.000000000', '2001-01-02T00:00:00.000000000',\n", " '2001-01-03T00:00:00.000000000', ..., '2003-12-29T00:00:00.000000000',\n", " '2003-12-30T00:00:00.000000000', '2003-12-31T00:00:00.000000000'],\n", " dtype='datetime64[ns]')
array([-90. , -86.326531, -82.653061, -78.979592, -75.306122, -71.632653,\n", " -67.959184, -64.285714, -60.612245, -56.938776, -53.265306, -49.591837,\n", " -45.918367, -42.244898, -38.571429, -34.897959, -31.22449 , -27.55102 ,\n", " -23.877551, -20.204082, -16.530612, -12.857143, -9.183673, -5.510204,\n", " -1.836735, 1.836735, 5.510204, 9.183673, 12.857143, 16.530612,\n", " 20.204082, 23.877551, 27.55102 , 31.22449 , 34.897959, 38.571429,\n", " 42.244898, 45.918367, 49.591837, 53.265306, 56.938776, 60.612245,\n", " 64.285714, 67.959184, 71.632653, 75.306122, 78.979592, 82.653061,\n", " 86.326531, 90. ])
array([-180. , -176.4, -172.8, -169.2, -165.6, -162. , -158.4, -154.8, -151.2,\n", " -147.6, -144. , -140.4, -136.8, -133.2, -129.6, -126. , -122.4, -118.8,\n", " -115.2, -111.6, -108. , -104.4, -100.8, -97.2, -93.6, -90. , -86.4,\n", " -82.8, -79.2, -75.6, -72. , -68.4, -64.8, -61.2, -57.6, -54. ,\n", " -50.4, -46.8, -43.2, -39.6, -36. , -32.4, -28.8, -25.2, -21.6,\n", " -18. , -14.4, -10.8, -7.2, -3.6, 0. , 3.6, 7.2, 10.8,\n", " 14.4, 18. , 21.6, 25.2, 28.8, 32.4, 36. , 39.6, 43.2,\n", " 46.8, 50.4, 54. , 57.6, 61.2, 64.8, 68.4, 72. , 75.6,\n", " 79.2, 82.8, 86.4, 90. , 93.6, 97.2, 100.8, 104.4, 108. ,\n", " 111.6, 115.2, 118.8, 122.4, 126. , 129.6, 133.2, 136.8, 140.4,\n", " 144. , 147.6, 151.2, 154.8, 158.4, 162. , 165.6, 169.2, 172.8,\n", " 176.4])
\n", " | time | \n", "lat | \n", "lon | \n", "event_duration | \n", "
---|---|---|---|---|
0 | \n", "2 | \n", "23 | \n", "12 | \n", "10 | \n", "
1 | \n", "3 | \n", "1 | \n", "12 | \n", "11 | \n", "
2 | \n", "3 | \n", "9 | \n", "21 | \n", "11 | \n", "
3 | \n", "3 | \n", "20 | \n", "6 | \n", "11 | \n", "
4 | \n", "6 | \n", "0 | \n", "6 | \n", "10 | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
2589 | \n", "1074 | \n", "31 | \n", "82 | \n", "11 | \n", "
2590 | \n", "1079 | \n", "36 | \n", "79 | \n", "11 | \n", "
2591 | \n", "1081 | \n", "34 | \n", "98 | \n", "12 | \n", "
2592 | \n", "1084 | \n", "25 | \n", "79 | \n", "10 | \n", "
2593 | \n", "1083 | \n", "35 | \n", "78 | \n", "11 | \n", "
2594 rows × 4 columns
\n", "\n", " | time | \n", "lat | \n", "lon | \n", "event_duration | \n", "
---|---|---|---|---|
0 | \n", "2001-01-03 | \n", "-5.510204 | \n", "-136.8 | \n", "10 days | \n", "
1 | \n", "2001-01-04 | \n", "-86.326531 | \n", "-136.8 | \n", "11 days | \n", "
2 | \n", "2001-01-04 | \n", "-56.938776 | \n", "-104.4 | \n", "11 days | \n", "
3 | \n", "2001-01-04 | \n", "-16.530612 | \n", "-158.4 | \n", "11 days | \n", "
4 | \n", "2001-01-07 | \n", "-90.000000 | \n", "-158.4 | \n", "10 days | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
2589 | \n", "2003-12-11 | \n", "23.877551 | \n", "115.2 | \n", "11 days | \n", "
2590 | \n", "2003-12-16 | \n", "42.244898 | \n", "104.4 | \n", "11 days | \n", "
2591 | \n", "2003-12-18 | \n", "34.897959 | \n", "172.8 | \n", "12 days | \n", "
2592 | \n", "2003-12-21 | \n", "1.836735 | \n", "104.4 | \n", "10 days | \n", "
2593 | \n", "2003-12-20 | \n", "38.571429 | \n", "100.8 | \n", "11 days | \n", "
2594 rows × 4 columns
\n", "\n", " | sum | \n", "mean | \n", "
---|---|---|
0 | \n", "7.720646191765621 | \n", "0.7720646191765621 | \n", "
1 | \n", "7.545168908782278 | \n", "0.6859244462529344 | \n", "
2 | \n", "8.776404953284827 | \n", "0.7978549957531661 | \n", "
3 | \n", "8.820629283522708 | \n", "0.8018753894111552 | \n", "
4 | \n", "7.350620818505183 | \n", "0.7350620818505182 | \n", "
... | \n", "... | \n", "... | \n", "
2589 | \n", "8.710790884253576 | \n", "0.7918900803866887 | \n", "
2590 | \n", "6.73266081786626 | \n", "0.6120600743514782 | \n", "
2591 | \n", "8.868592214995974 | \n", "0.7390493512496645 | \n", "
2592 | \n", "7.790211594277864 | \n", "0.7790211594277864 | \n", "
2593 | \n", "8.669867314918553 | \n", "0.7881697559016867 | \n", "
2594 rows × 2 columns
\n", "\n", " | time | \n", "lat | \n", "lon | \n", "event_duration | \n", "sum | \n", "mean | \n", "
---|---|---|---|---|---|---|
0 | \n", "2001-01-03 | \n", "-5.510204 | \n", "-136.8 | \n", "10 days | \n", "7.720646191765621 | \n", "0.7720646191765621 | \n", "
1 | \n", "2001-01-04 | \n", "-86.326531 | \n", "-136.8 | \n", "11 days | \n", "7.545168908782278 | \n", "0.6859244462529344 | \n", "
2 | \n", "2001-01-04 | \n", "-56.938776 | \n", "-104.4 | \n", "11 days | \n", "8.776404953284827 | \n", "0.7978549957531661 | \n", "
3 | \n", "2001-01-04 | \n", "-16.530612 | \n", "-158.4 | \n", "11 days | \n", "8.820629283522708 | \n", "0.8018753894111552 | \n", "
4 | \n", "2001-01-07 | \n", "-90.000000 | \n", "-158.4 | \n", "10 days | \n", "7.350620818505183 | \n", "0.7350620818505182 | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
2589 | \n", "2003-12-11 | \n", "23.877551 | \n", "115.2 | \n", "11 days | \n", "8.710790884253576 | \n", "0.7918900803866887 | \n", "
2590 | \n", "2003-12-16 | \n", "42.244898 | \n", "104.4 | \n", "11 days | \n", "6.73266081786626 | \n", "0.6120600743514782 | \n", "
2591 | \n", "2003-12-18 | \n", "34.897959 | \n", "172.8 | \n", "12 days | \n", "8.868592214995974 | \n", "0.7390493512496645 | \n", "
2592 | \n", "2003-12-21 | \n", "1.836735 | \n", "104.4 | \n", "10 days | \n", "7.790211594277864 | \n", "0.7790211594277864 | \n", "
2593 | \n", "2003-12-20 | \n", "38.571429 | \n", "100.8 | \n", "11 days | \n", "8.669867314918553 | \n", "0.7881697559016867 | \n", "
2594 rows × 6 columns
\n", "