diff --git a/docs/pangeo/xarray_introduction.ipynb b/docs/pangeo/xarray_introduction.ipynb old mode 100755 new mode 100644 index a47e588..8984f62 --- a/docs/pangeo/xarray_introduction.ipynb +++ b/docs/pangeo/xarray_introduction.ipynb @@ -64,15 +64,37 @@ "source": [ "## Context\n", "\n", - "We will be using the [Pangeo](https://pangeo.io/) open-source software stack for visualizing the near-surface temperature and computing time averaged values (such as seasonal mean and other statistics).\n", + "We will be using the [Pangeo](https://pangeo.io/) open-source software stack for visualizing the Particle Matter < 2.5 μm and computing time averaged values (such as daily mean and other statistics).\n", "\n", "### Data\n", "\n", - "In this episode, we will use [CMIP6](https://www.wcrp-climate.org/wgcm-cmip/wgcm-cmip6) data.\n", + "We will be using data from [Copernicus Atmosphere Monitoring Service](https://ads.atmosphere.copernicus.eu/)\n", + "and more precisely PM2.5 ([Particle Matter < 2.5 μm](https://en.wikipedia.org/wiki/Particulates#Size,_shape_and_solubility_matter)) 4 days forecast from December, 22 2021.\n", "\n", - "This dataset can be discovered through the [CMIP6 online catalog](https://pangeo-data.github.io/pangeo-cmip6-cloud/) or from [ESGF](https://esgf.llnl.gov/).\n", - "\n", - "The same dataset can also be downloaded from [Zenodo](https://zenodo.org/): [Near-surface Temperature from CMIP6 NCAR CESM2 historical monthly dataset for CLIVAR CMIP6 Bootcamp](https://zenodo.org/record/7181714/files/CMIP6_NCAR_CESM2_historical_amon_gn.nc)." + "The dataset can be downloaded from [Zenodo](https://zenodo.org/): [PM2.5 4 days forecast from December, 22 2020 retrieved from Copernicus Monitoring Service](https://zenodo.org/records/5805953)." + ] + }, + { + "cell_type": "markdown", + "id": "6d69e8e9-5e58-460b-bb67-be18fd5cfd08", + "metadata": {}, + "source": [ + "
\n", + "\n", + "Comment: Remark\n", + "This tutorial uses data on a regular latitude-longitude grid. More complex and irregular grids are not discussed in this tutorial. In addition,\n", + "this tutorial is not meant to cover all the different possibilities offered by Xarrays but shows functionalities we find useful for day to day\n", + "analysis.
\n", + "
\n", + "" ] }, { @@ -87,6 +109,7 @@ "- xarray {cite:ps}`a-xarray-hoyer2017` with [`netCDF4`](https://pypi.org/project/h5netcdf/) and [`h5netcdf`](https://pypi.org/project/h5netcdf/) engines\n", "- pooch {cite:ps}`a-pooch-Uieda2020`\n", "- numpy {cite:ps}`a-numpy-harris2020`\n", + "- cmcrameri {cite:ps}`a-cmcrameri-crameri2018`\n", "\n", "Please install these packages if they are not already available in your Python environment (see [Setup page](https://pangeo-data.github.io/foss4g-2022/before/setup.html)).\n", "\n", @@ -98,6 +121,36 @@ { "cell_type": "code", "execution_count": 1, + "id": "dfb32910-f236-4d50-88e0-70871bba51e6", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Requirement already satisfied: cmcrameri in /srv/conda/envs/notebook/lib/python3.11/site-packages (1.9)\n", + "Requirement already satisfied: matplotlib in /srv/conda/envs/notebook/lib/python3.11/site-packages (from cmcrameri) (3.8.2)\n", + "Requirement already satisfied: numpy in /srv/conda/envs/notebook/lib/python3.11/site-packages (from cmcrameri) (1.26.2)\n", + "Requirement already satisfied: packaging in /srv/conda/envs/notebook/lib/python3.11/site-packages (from cmcrameri) (23.2)\n", + "Requirement already satisfied: contourpy>=1.0.1 in /srv/conda/envs/notebook/lib/python3.11/site-packages (from matplotlib->cmcrameri) (1.2.0)\n", + "Requirement already satisfied: cycler>=0.10 in /srv/conda/envs/notebook/lib/python3.11/site-packages (from matplotlib->cmcrameri) (0.12.1)\n", + "Requirement already satisfied: fonttools>=4.22.0 in /srv/conda/envs/notebook/lib/python3.11/site-packages (from matplotlib->cmcrameri) (4.46.0)\n", + "Requirement already satisfied: kiwisolver>=1.3.1 in /srv/conda/envs/notebook/lib/python3.11/site-packages (from matplotlib->cmcrameri) (1.4.5)\n", + "Requirement already satisfied: pillow>=8 in /srv/conda/envs/notebook/lib/python3.11/site-packages (from matplotlib->cmcrameri) (10.1.0)\n", + "Requirement already satisfied: pyparsing>=2.3.1 in /srv/conda/envs/notebook/lib/python3.11/site-packages (from matplotlib->cmcrameri) (3.1.1)\n", + "Requirement already satisfied: python-dateutil>=2.7 in /srv/conda/envs/notebook/lib/python3.11/site-packages (from matplotlib->cmcrameri) (2.8.2)\n", + "Requirement already satisfied: six>=1.5 in /srv/conda/envs/notebook/lib/python3.11/site-packages (from python-dateutil>=2.7->matplotlib->cmcrameri) (1.16.0)\n", + "Note: you may need to restart the kernel to use updated packages.\n" + ] + } + ], + "source": [ + "pip install cmcrameri" + ] + }, + { + "cell_type": "code", + "execution_count": 2, "id": "51c76db6-f5aa-40eb-97f1-2df7cde50e83", "metadata": {}, "outputs": [], @@ -119,7 +172,7 @@ }, { "cell_type": "code", - "execution_count": 2, + "execution_count": 3, "id": "08ba3f1d-71ea-42ec-904f-dc2452640590", "metadata": {}, "outputs": [], @@ -129,7 +182,7 @@ }, { "cell_type": "code", - "execution_count": 3, + "execution_count": 4, "id": "f247c3cd-7f19-43c8-b721-6b159f401379", "metadata": { "tags": [ @@ -138,9 +191,9 @@ }, "outputs": [], "source": [ - "tas_file = pooch.retrieve(\n", - " url=\"https://zenodo.org/record/7181714/files/CMIP6_NCAR_CESM2_historical_amon_gn.nc\",\n", - " known_hash=\"md5:5f86251e5bc5ef9b86a3a86cd06a536b\",\n", + "cams_file = pooch.retrieve(\n", + " url=\"https://zenodo.org/record/5805953/files/CAMS-PM2_5-20211222.netcdf\",\n", + " known_hash=\"md5:c4a6bb0a5a5640fc8de2ae6f377932fc\",\n", " path=f\".\",\n", ")" ] @@ -155,12 +208,12 @@ }, { "cell_type": "code", - "execution_count": 4, + "execution_count": 5, "id": "6974bdbf-f874-4b0d-bb9b-13a0161af3b0", "metadata": {}, "outputs": [], "source": [ - "tas_ds = xr.open_dataset(tas_file)" + "cams = xr.open_dataset(cams_file)" ] }, { @@ -179,13 +232,13 @@ "metadata": {}, "source": [ ":::{tip}\n", - "If you get an error with the previous command, first check the location of the input file some_hash-CMIP6_NCAR_CESM2_historical_amon_gn.nc: it should have been downloaded in the same directory as your Jupyter Notebook.\n", + "If you get an error with the previous command, first check the location of the input file some_hash-CAMS-PM2_5-20211222.netcdf: it should have been downloaded in the same directory as your Jupyter Notebook.\n", ":::" ] }, { "cell_type": "code", - "execution_count": 5, + "execution_count": 6, "id": "098a2be0", "metadata": {}, "outputs": [ @@ -453,6 +506,11 @@ " grid-column: 4;\n", "}\n", "\n", + ".xr-index-preview {\n", + " grid-column: 2 / 5;\n", + " color: var(--xr-font-color2);\n", + "}\n", + "\n", ".xr-var-name,\n", ".xr-var-dims,\n", ".xr-var-dtype,\n", @@ -474,14 +532,16 @@ "}\n", "\n", ".xr-var-attrs,\n", - ".xr-var-data {\n", + ".xr-var-data,\n", + ".xr-index-data {\n", " display: none;\n", " background-color: var(--xr-background-color) !important;\n", " padding-bottom: 5px !important;\n", "}\n", "\n", ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n", - ".xr-var-data-in:checked ~ .xr-var-data {\n", + ".xr-var-data-in:checked ~ .xr-var-data,\n", + ".xr-index-data-in:checked ~ .xr-index-data {\n", " display: block;\n", "}\n", "\n", @@ -491,13 +551,16 @@ "\n", ".xr-var-name span,\n", ".xr-var-data,\n", + ".xr-index-name div,\n", + ".xr-index-data,\n", ".xr-attrs {\n", " padding-left: 25px !important;\n", "}\n", "\n", ".xr-attrs,\n", ".xr-var-attrs,\n", - ".xr-var-data {\n", + ".xr-var-data,\n", + ".xr-index-data {\n", " grid-column: 1 / -1;\n", "}\n", "\n", @@ -535,7 +598,8 @@ "}\n", "\n", ".xr-icon-database,\n", - ".xr-icon-file-text2 {\n", + ".xr-icon-file-text2,\n", + ".xr-no-icon {\n", " display: inline-block;\n", " vertical-align: middle;\n", " width: 1em;\n", @@ -545,134 +609,125 @@ " fill: currentColor;\n", "}\n", "Agenda\n", + "In this tutorial, we will cover:
\n", + "\n", + "
\n", + "- Analysis
\n", + "\n", + "
\n", + "- Import Python packages
\n", + "
<xarray.Dataset>\n", - "Dimensions: (lat: 192, nbnd: 2, lon: 288, member_id: 1, time: 1980)\n", + "Dimensions: (longitude: 700, latitude: 400, level: 1, time: 97)\n", "Coordinates:\n", - " * lat (lat) float64 -90.0 -89.06 -88.12 -87.17 ... 88.12 89.06 90.0\n", - " lat_bnds (lat, nbnd) float32 -90.0 -89.53 -89.53 ... 89.53 89.53 90.0\n", - " * lon (lon) float64 0.0 1.25 2.5 3.75 5.0 ... 355.0 356.2 357.5 358.8\n", - " lon_bnds (lon, nbnd) float32 -0.625 0.625 0.625 ... 358.1 358.1 359.4\n", - " * time (time) object 1850-01-15 12:00:00 ... 2014-12-15 12:00:00\n", - " time_bnds (time, nbnd) object 1850-01-01 00:00:00 ... 2015-01-01 00:00:00\n", - " * member_id (member_id) object 'r1i1p1f1'\n", - "Dimensions without coordinates: nbnd\n", + " * longitude (longitude) float32 335.0 335.1 335.2 ... 44.75 44.85 44.95\n", + " * latitude (latitude) float32 69.95 69.85 69.75 69.65 ... 30.25 30.15 30.05\n", + " * level (level) float32 0.0\n", + " * time (time) timedelta64[ns] 00:00:00 01:00:00 ... 4 days 00:00:00\n", "Data variables:\n", - " tas (member_id, time, lat, lon) float32 ...\n", - "Attributes: (12/50)\n", - " Conventions: CF-1.7 CMIP-6.2\n", - " activity_id: CMIP\n", - " branch_method: standard\n", - " branch_time_in_child: 674885.0\n", - " branch_time_in_parent: 219000.0\n", - " case_id: 15\n", - " ... ...\n", - " variant_label: r1i1p1f1\n", - " status: 2019-10-25;created;by nhn2@columbia.edu\n", - " netcdf_tracking_ids: hdl:21.14100/d9a7225a-49c3-4470-b7ab-a8180926f839\n", - " version_id: v20190308\n", - " intake_esm_varname: tas\n", - " intake_esm_dataset_key: CMIP.NCAR.CESM2.historical.Amon.gn" + " pm2p5_conc (time, level, latitude, longitude) float32 ...\n", + "Attributes:\n", + " title: PM25 Air Pollutant FORECAST at the Surface\n", + " institution: Data produced by Meteo France\n", + " source: Data from ENSEMBLE model\n", + " history: Model ENSEMBLE FORECAST\n", + " FORECAST: Europe, 20211222+[0H_96H]\n", + " summary: ENSEMBLE model hourly FORECAST of PM25 concentration at the...\n", + " project: MACC-RAQ (http://macc-raq.gmes-atmosphere.eu) " ], "text/plain": [ "
<xarray.DataArray 'time' (time: 1980)>\n", - "array([cftime.DatetimeNoLeap(1850, 1, 15, 12, 0, 0, 0, has_year_zero=True),\n", - " cftime.DatetimeNoLeap(1850, 2, 14, 0, 0, 0, 0, has_year_zero=True),\n", - " cftime.DatetimeNoLeap(1850, 3, 15, 12, 0, 0, 0, has_year_zero=True),\n", - " ...,\n", - " cftime.DatetimeNoLeap(2014, 10, 15, 12, 0, 0, 0, has_year_zero=True),\n", - " cftime.DatetimeNoLeap(2014, 11, 15, 0, 0, 0, 0, has_year_zero=True),\n", - " cftime.DatetimeNoLeap(2014, 12, 15, 12, 0, 0, 0, has_year_zero=True)],\n", - " dtype=object)\n", + "" ], "text/plain": [ - "<xarray.DataArray 'time' (time: 97)>\n", + "array([ 0, 3600000000000, 7200000000000, 10800000000000,\n", + " 14400000000000, 18000000000000, 21600000000000, 25200000000000,\n", + " 28800000000000, 32400000000000, 36000000000000, 39600000000000,\n", + " 43200000000000, 46800000000000, 50400000000000, 54000000000000,\n", + " 57600000000000, 61200000000000, 64800000000000, 68400000000000,\n", + " 72000000000000, 75600000000000, 79200000000000, 82800000000000,\n", + " 86400000000000, 90000000000000, 93600000000000, 97200000000000,\n", + " 100800000000000, 104400000000000, 108000000000000, 111600000000000,\n", + " 115200000000000, 118800000000000, 122400000000000, 126000000000000,\n", + " 129600000000000, 133200000000000, 136800000000000, 140400000000000,\n", + " 144000000000000, 147600000000000, 151200000000000, 154800000000000,\n", + " 158400000000000, 162000000000000, 165600000000000, 169200000000000,\n", + " 172800000000000, 176400000000000, 180000000000000, 183600000000000,\n", + " 187200000000000, 190800000000000, 194400000000000, 198000000000000,\n", + " 201600000000000, 205200000000000, 208800000000000, 212400000000000,\n", + " 216000000000000, 219600000000000, 223200000000000, 226800000000000,\n", + " 230400000000000, 234000000000000, 237600000000000, 241200000000000,\n", + " 244800000000000, 248400000000000, 252000000000000, 255600000000000,\n", + " 259200000000000, 262800000000000, 266400000000000, 270000000000000,\n", + " 273600000000000, 277200000000000, 280800000000000, 284400000000000,\n", + " 288000000000000, 291600000000000, 295200000000000, 298800000000000,\n", + " 302400000000000, 306000000000000, 309600000000000, 313200000000000,\n", + " 316800000000000, 320400000000000, 324000000000000, 327600000000000,\n", + " 331200000000000, 334800000000000, 338400000000000, 342000000000000,\n", + " 345600000000000], dtype='timedelta64[ns]')\n", "Coordinates:\n", - " * time (time) object 1850-01-15 12:00:00 ... 2014-12-15 12:00:00\n", + " * time (time) timedelta64[ns] 00:00:00 01:00:00 ... 4 days 00:00:00\n", "Attributes:\n", - " axis: T\n", - " bounds: time_bnds\n", - " standard_name: time\n", - " title: time\n", - " type: double" + " long_name: FORECAST time from 20211222
<xarray.DataArray 'tas' (member_id: 1, time: 1980, lat: 192, lon: 288)>\n", - "[109486080 values with dtype=float32]\n", + "" ], "text/plain": [ - "<xarray.DataArray 'pm2p5_conc' (time: 97, level: 1, latitude: 400,\n", + " longitude: 700)>\n", + "[27160000 values with dtype=float32]\n", "Coordinates:\n", - " * lat (lat) float64 -90.0 -89.06 -88.12 -87.17 ... 88.12 89.06 90.0\n", - " * lon (lon) float64 0.0 1.25 2.5 3.75 5.0 ... 355.0 356.2 357.5 358.8\n", - " * time (time) object 1850-01-15 12:00:00 ... 2014-12-15 12:00:00\n", - " * member_id (member_id) object 'r1i1p1f1'\n", + " * longitude (longitude) float32 335.0 335.1 335.2 335.4 ... 44.75 44.85 44.95\n", + " * latitude (latitude) float32 69.95 69.85 69.75 69.65 ... 30.25 30.15 30.05\n", + " * level (level) float32 0.0\n", + " * time (time) timedelta64[ns] 00:00:00 01:00:00 ... 4 days 00:00:00\n", "Attributes:\n", - " cell_measures: area: areacella\n", - " cell_methods: area: time: mean\n", - " comment: near-surface (usually, 2 meter) air temperature\n", - " description: near-surface (usually, 2 meter) air temperature\n", - " frequency: mon\n", - " id: tas\n", - " long_name: Near-Surface Air Temperature\n", - " mipTable: Amon\n", - " out_name: tas\n", - " prov: Amon ((isd.003))\n", - " realm: atmos\n", - " standard_name: air_temperature\n", - " time: time\n", - " time_label: time-mean\n", - " time_title: Temporal mean\n", - " title: Near-Surface Air Temperature\n", - " type: real\n", - " units: K\n", - " variable_id: tas" + " species: PM2.5 Aerosol\n", + " units: µg/m3\n", + " value: hourly values\n", + " standard_name: mass_concentration_of_pm2p5_ambient_aerosol_in_air
<xarray.DataArray 'tas' (member_id: 1, time: 1980, lat: 192, lon: 288)>\n", - "[109486080 values with dtype=float32]\n", + "" ], "text/plain": [ - "<xarray.DataArray 'pm2p5_conc' (time: 97, level: 1, latitude: 400,\n", + " longitude: 700)>\n", + "[27160000 values with dtype=float32]\n", "Coordinates:\n", - " * lat (lat) float64 -90.0 -89.06 -88.12 -87.17 ... 88.12 89.06 90.0\n", - " * lon (lon) float64 0.0 1.25 2.5 3.75 5.0 ... 355.0 356.2 357.5 358.8\n", - " * time (time) object 1850-01-15 12:00:00 ... 2014-12-15 12:00:00\n", - " * member_id (member_id) object 'r1i1p1f1'\n", + " * longitude (longitude) float32 335.0 335.1 335.2 335.4 ... 44.75 44.85 44.95\n", + " * latitude (latitude) float32 69.95 69.85 69.75 69.65 ... 30.25 30.15 30.05\n", + " * level (level) float32 0.0\n", + " * time (time) timedelta64[ns] 00:00:00 01:00:00 ... 4 days 00:00:00\n", "Attributes:\n", - " cell_measures: area: areacella\n", - " cell_methods: area: time: mean\n", - " comment: near-surface (usually, 2 meter) air temperature\n", - " description: near-surface (usually, 2 meter) air temperature\n", - " frequency: mon\n", - " id: tas\n", - " long_name: Near-Surface Air Temperature\n", - " mipTable: Amon\n", - " out_name: tas\n", - " prov: Amon ((isd.003))\n", - " realm: atmos\n", - " standard_name: air_temperature\n", - " time: time\n", - " time_label: time-mean\n", - " time_title: Temporal mean\n", - " title: Near-Surface Air Temperature\n", - " type: real\n", - " units: K\n", - " variable_id: tas" + " species: PM2.5 Aerosol\n", + " units: µg/m3\n", + " value: hourly values\n", + " standard_name: mass_concentration_of_pm2p5_ambient_aerosol_in_air
<xarray.Dataset>\n", - "Dimensions: (latitude: 192, nbnd: 2, longitude: 288, member_id: 1, time: 1980)\n", + "Dimensions: (lon: 700, lat: 400, level: 1, time: 97)\n", "Coordinates:\n", - " * latitude (latitude) float64 -90.0 -89.06 -88.12 ... 88.12 89.06 90.0\n", - " lat_bnds (latitude, nbnd) float32 -90.0 -89.53 -89.53 ... 89.53 89.53 90.0\n", - " * longitude (longitude) float64 0.0 1.25 2.5 3.75 ... 355.0 356.2 357.5 358.8\n", - " lon_bnds (longitude, nbnd) float32 -0.625 0.625 0.625 ... 358.1 359.4\n", - " * time (time) object 1850-01-15 12:00:00 ... 2014-12-15 12:00:00\n", - " time_bnds (time, nbnd) object 1850-01-01 00:00:00 ... 2015-01-01 00:00:00\n", - " * member_id (member_id) object 'r1i1p1f1'\n", - "Dimensions without coordinates: nbnd\n", + " * lon (lon) float32 335.0 335.1 335.2 335.4 ... 44.75 44.85 44.95\n", + " * lat (lat) float32 69.95 69.85 69.75 69.65 ... 30.25 30.15 30.05\n", + " * level (level) float32 0.0\n", + " * time (time) timedelta64[ns] 00:00:00 01:00:00 ... 4 days 00:00:00\n", "Data variables:\n", - " tas (member_id, time, latitude, longitude) float32 245.3 ... 245.0\n", - "Attributes: (12/50)\n", - " Conventions: CF-1.7 CMIP-6.2\n", - " activity_id: CMIP\n", - " branch_method: standard\n", - " branch_time_in_child: 674885.0\n", - " branch_time_in_parent: 219000.0\n", - " case_id: 15\n", - " ... ...\n", - " variant_label: r1i1p1f1\n", - " status: 2019-10-25;created;by nhn2@columbia.edu\n", - " netcdf_tracking_ids: hdl:21.14100/d9a7225a-49c3-4470-b7ab-a8180926f839\n", - " version_id: v20190308\n", - " intake_esm_varname: tas\n", - " intake_esm_dataset_key: CMIP.NCAR.CESM2.historical.Amon.gn
PandasIndex(Index([335.04998779296875, 335.1499938964844, 335.25,\n", + " 335.3500061035156, 335.45001220703125, 335.54998779296875,\n", + " 335.6499938964844, 335.75, 335.8500061035156,\n", + " 335.95001220703125,\n", + " ...\n", + " 44.04999923706055, 44.150001525878906, 44.25,\n", + " 44.349998474121094, 44.45000076293945, 44.54999923706055,\n", + " 44.650001525878906, 44.75, 44.849998474121094,\n", + " 44.95000076293945],\n", + " dtype='float32', name='lon', length=700))
PandasIndex(Index([ 69.94999694824219, 69.8499984741211, 69.75,\n", + " 69.6500015258789, 69.55000305175781, 69.44999694824219,\n", + " 69.3499984741211, 69.25, 69.1500015258789,\n", + " 69.05000305175781,\n", + " ...\n", + " 30.950000762939453, 30.850000381469727, 30.75,\n", + " 30.649999618530273, 30.549999237060547, 30.450000762939453,\n", + " 30.350000381469727, 30.25, 30.149999618530273,\n", + " 30.049999237060547],\n", + " dtype='float32', name='lat', length=400))
PandasIndex(Index([0.0], dtype='float32', name='level'))
PandasIndex(TimedeltaIndex(['0 days 00:00:00', '0 days 01:00:00', '0 days 02:00:00',\n", + " '0 days 03:00:00', '0 days 04:00:00', '0 days 05:00:00',\n", + " '0 days 06:00:00', '0 days 07:00:00', '0 days 08:00:00',\n", + " '0 days 09:00:00', '0 days 10:00:00', '0 days 11:00:00',\n", + " '0 days 12:00:00', '0 days 13:00:00', '0 days 14:00:00',\n", + " '0 days 15:00:00', '0 days 16:00:00', '0 days 17:00:00',\n", + " '0 days 18:00:00', '0 days 19:00:00', '0 days 20:00:00',\n", + " '0 days 21:00:00', '0 days 22:00:00', '0 days 23:00:00',\n", + " '1 days 00:00:00', '1 days 01:00:00', '1 days 02:00:00',\n", + " '1 days 03:00:00', '1 days 04:00:00', '1 days 05:00:00',\n", + " '1 days 06:00:00', '1 days 07:00:00', '1 days 08:00:00',\n", + " '1 days 09:00:00', '1 days 10:00:00', '1 days 11:00:00',\n", + " '1 days 12:00:00', '1 days 13:00:00', '1 days 14:00:00',\n", + " '1 days 15:00:00', '1 days 16:00:00', '1 days 17:00:00',\n", + " '1 days 18:00:00', '1 days 19:00:00', '1 days 20:00:00',\n", + " '1 days 21:00:00', '1 days 22:00:00', '1 days 23:00:00',\n", + " '2 days 00:00:00', '2 days 01:00:00', '2 days 02:00:00',\n", + " '2 days 03:00:00', '2 days 04:00:00', '2 days 05:00:00',\n", + " '2 days 06:00:00', '2 days 07:00:00', '2 days 08:00:00',\n", + " '2 days 09:00:00', '2 days 10:00:00', '2 days 11:00:00',\n", + " '2 days 12:00:00', '2 days 13:00:00', '2 days 14:00:00',\n", + " '2 days 15:00:00', '2 days 16:00:00', '2 days 17:00:00',\n", + " '2 days 18:00:00', '2 days 19:00:00', '2 days 20:00:00',\n", + " '2 days 21:00:00', '2 days 22:00:00', '2 days 23:00:00',\n", + " '3 days 00:00:00', '3 days 01:00:00', '3 days 02:00:00',\n", + " '3 days 03:00:00', '3 days 04:00:00', '3 days 05:00:00',\n", + " '3 days 06:00:00', '3 days 07:00:00', '3 days 08:00:00',\n", + " '3 days 09:00:00', '3 days 10:00:00', '3 days 11:00:00',\n", + " '3 days 12:00:00', '3 days 13:00:00', '3 days 14:00:00',\n", + " '3 days 15:00:00', '3 days 16:00:00', '3 days 17:00:00',\n", + " '3 days 18:00:00', '3 days 19:00:00', '3 days 20:00:00',\n", + " '3 days 21:00:00', '3 days 22:00:00', '3 days 23:00:00',\n", + " '4 days 00:00:00'],\n", + " dtype='timedelta64[ns]', name='time', freq=None))
<xarray.Dataset>\n", - "Dimensions: (lat: 192, nbnd: 2, lon: 288, member_id: 1, time: 1980)\n", + "" ], "text/plain": [ - "<xarray.DataArray 'pm2p5_conc' ()>\n", + "array(2.396931, dtype=float32)\n", "Coordinates:\n", - " * lat (lat) float64 -90.0 -89.06 -88.12 -87.17 ... 88.12 89.06 90.0\n", - " lat_bnds (lat, nbnd) float32 -90.0 -89.53 -89.53 ... 89.53 89.53 90.0\n", - " * lon (lon) float64 0.0 1.25 2.5 3.75 5.0 ... 355.0 356.2 357.5 358.8\n", - " lon_bnds (lon, nbnd) float32 -0.625 0.625 0.625 ... 358.1 358.1 359.4\n", - " * time (time) object 1850-01-15 12:00:00 ... 2014-12-15 12:00:00\n", - " time_bnds (time, nbnd) object 1850-01-01 00:00:00 ... 2015-01-01 00:00:00\n", - " * member_id (member_id) object 'r1i1p1f1'\n", - "Dimensions without coordinates: nbnd\n", - "Data variables:\n", - " tas (member_id, time, lat, lon) float32 245.3 245.3 ... 245.0 245.0\n", - "Attributes: (12/50)\n", - " Conventions: CF-1.7 CMIP-6.2\n", - " activity_id: CMIP\n", - " branch_method: standard\n", - " branch_time_in_child: 674885.0\n", - " branch_time_in_parent: 219000.0\n", - " case_id: 15\n", - " ... ...\n", - " variant_label: r1i1p1f1\n", - " status: 2019-10-25;created;by nhn2@columbia.edu\n", - " netcdf_tracking_ids: hdl:21.14100/d9a7225a-49c3-4470-b7ab-a8180926f839\n", - " version_id: v20190308\n", - " intake_esm_varname: tas\n", - " intake_esm_dataset_key: CMIP.NCAR.CESM2.historical.Amon.gn" + " lat float32 59.95\n", + " level float32 0.0\n", + " time timedelta64[ns] 00:00:00\n", + " lon float32 -14.95\n", + "Attributes:\n", + " species: PM2.5 Aerosol\n", + " units: µg/m3\n", + " value: hourly values\n", + " standard_name: mass_concentration_of_pm2p5_ambient_aerosol_in_air
<xarray.DataArray 'tas' ()>\n", - "array(300.38446, dtype=float32)\n", + "" ], "text/plain": [ - "<xarray.DataArray 'pm2p5_conc' (lat: 10, lon: 10)>\n", + "array([[2.396931, 2.40307 , 2.406154, 2.408797, 2.415234, 2.418361, 2.414907,\n", + " 2.408982, 2.398707, 2.388873],\n", + " [2.175199, 2.175696, 2.172883, 2.173877, 2.177146, 2.17976 , 2.219366,\n", + " 2.337629, 2.281311, 2.306678],\n", + " [2.153712, 2.11665 , 2.083326, 2.175028, 2.317506, 2.509566, 2.59864 ,\n", + " 2.742937, 2.716505, 2.692815],\n", + " [2.673005, 2.624574, 2.577082, 2.55173 , 2.527671, 2.507051, 2.495241,\n", + " 2.491348, 2.499505, 2.51255 ],\n", + " [3.034373, 3.064997, 3.083713, 3.099743, 3.114579, 3.122835, 3.130595,\n", + " 3.134872, 3.066902, 3.080857],\n", + " [2.827477, 2.889905, 2.940652, 3.008623, 3.070142, 3.122423, 3.16949 ,\n", + " 3.204718, 3.228834, 3.248346],\n", + " [2.72496 , 2.726921, 2.744543, 2.776133, 2.823825, 2.872682, 2.926299,\n", + " 3.074007, 3.489817, 3.465331],\n", + " [2.830262, 2.818837, 2.809443, 2.803788, 2.796654, 3.25818 , 3.697267,\n", + " 3.924114, 3.825889, 3.679787],\n", + " [4.127571, 3.362629, 3.57602 , 3.699882, 3.674302, 3.445621, 3.212406,\n", + " 3.258507, 3.308671, 3.361635],\n", + " [4.588827, 4.447983, 4.203201, 3.700336, 3.73332 , 3.784365, 3.837968,\n", + " 3.880246, 3.91303 , 3.947676]], dtype=float32)\n", "Coordinates:\n", - " lat float64 4.241\n", - " lon float64 125.0\n", - " time object 1850-01-15 12:00:00\n", - " member_id <U8 'r1i1p1f1'\n", + " * lat (lat) float32 59.95 59.85 59.75 59.65 ... 59.35 59.25 59.15 59.05\n", + " level float32 0.0\n", + " time timedelta64[ns] 00:00:00\n", + " * lon (lon) float32 -14.95 -14.85 -14.75 -14.65 ... -14.25 -14.15 -14.05\n", "Attributes:\n", - " cell_measures: area: areacella\n", - " cell_methods: area: time: mean\n", - " comment: near-surface (usually, 2 meter) air temperature\n", - " description: near-surface (usually, 2 meter) air temperature\n", - " frequency: mon\n", - " id: tas\n", - " long_name: Near-Surface Air Temperature\n", - " mipTable: Amon\n", - " out_name: tas\n", - " prov: Amon ((isd.003))\n", - " realm: atmos\n", - " standard_name: air_temperature\n", - " time: time\n", - " time_label: time-mean\n", - " time_title: Temporal mean\n", - " title: Near-Surface Air Temperature\n", - " type: real\n", - " units: K\n", - " variable_id: tas" + " species: PM2.5 Aerosol\n", + " units: µg/m3\n", + " value: hourly values\n", + " standard_name: mass_concentration_of_pm2p5_ambient_aerosol_in_air
<xarray.DataArray 'tas' (time: 5, lat: 10, lon: 10)>\n", - "array([[[300.38446, 300.08688, ..., 300.71 , 300.7822 ],\n", - " [300.2441 , 299.95798, ..., 300.53357, 300.62396],\n", - " ...,\n", - " [297.77292, 298.71323, ..., 299.71173, 299.8677 ],\n", - " [297.73315, 298.60205, ..., 299.6976 , 299.812 ]],\n", - "\n", - " [[299.65933, 299.26968, ..., 300.4081 , 300.55255],\n", - " [299.5267 , 299.12592, ..., 300.40277, 300.5284 ],\n", - " ...,\n", - " [297.3739 , 298.0038 , ..., 299.71512, 299.80475],\n", - " [297.3448 , 297.89917, ..., 299.4545 , 299.56235]],\n", - "\n", - " ...,\n", - "\n", - " [[300.5706 , 300.41907, ..., 300.87616, 300.82022],\n", - " [300.4315 , 300.26505, ..., 300.84323, 300.8234 ],\n", - " ...,\n", - " [299.8353 , 299.99384, ..., 300.41165, 300.41684],\n", - " [299.85242, 299.9631 , ..., 300.28253, 300.27316]],\n", - "\n", - " [[300.81683, 300.7148 , ..., 301.1554 , 301.1403 ],\n", - " [300.5332 , 300.56613, ..., 300.99164, 301.0214 ],\n", - " ...,\n", - " [300.82352, 300.53986, ..., 300.4276 , 300.4221 ],\n", - " [300.87964, 300.58102, ..., 300.36655, 300.33893]]], dtype=float32)\n", + "" ], "text/plain": [ - "<xarray.DataArray 'pm2p5_conc' ()>\n", + "array(2.396931, dtype=float32)\n", "Coordinates:\n", - " * lat (lat) float64 4.241 5.183 6.126 7.068 ... 9.895 10.84 11.78 12.72\n", - " * lon (lon) float64 125.0 126.2 127.5 128.8 ... 132.5 133.8 135.0 136.2\n", - " * time (time) object 1850-01-15 12:00:00 ... 1850-05-15 12:00:00\n", - " member_id <U8 'r1i1p1f1'\n", + " lat float32 59.95\n", + " level float32 0.0\n", + " time timedelta64[ns] 00:00:00\n", + " lon float32 -14.95\n", "Attributes:\n", - " cell_measures: area: areacella\n", - " cell_methods: area: time: mean\n", - " comment: near-surface (usually, 2 meter) air temperature\n", - " description: near-surface (usually, 2 meter) air temperature\n", - " frequency: mon\n", - " id: tas\n", - " long_name: Near-Surface Air Temperature\n", - " mipTable: Amon\n", - " out_name: tas\n", - " prov: Amon ((isd.003))\n", - " realm: atmos\n", - " standard_name: air_temperature\n", - " time: time\n", - " time_label: time-mean\n", - " time_title: Temporal mean\n", - " title: Near-Surface Air Temperature\n", - " type: real\n", - " units: K\n", - " variable_id: tas" + " species: PM2.5 Aerosol\n", + " units: µg/m3\n", + " value: hourly values\n", + " standard_name: mass_concentration_of_pm2p5_ambient_aerosol_in_air
<xarray.DataArray 'tas' ()>\n", - "array(300.38446, dtype=float32)\n", + "" ], "text/plain": [ - "<xarray.DataArray 'pm2p5_conc' (level: 1, lat: 400, lon: 700)>\n", + "array([[[0.500001, 0.500001, ..., 1.419116, 1.425483],\n", + " [0.500001, 0.500001, ..., 1.238894, 1.257112],\n", + " ...,\n", + " [2.416178, 2.482116, ..., 9.897312, 9.944606],\n", + " [2.588073, 2.540324, ..., 9.613549, 9.609912]]], dtype=float32)\n", "Coordinates:\n", - " lat float64 4.241\n", - " lon float64 125.0\n", - " time object 1850-01-15 12:00:00\n", - " member_id <U8 'r1i1p1f1'\n", + " * lat (lat) float32 69.95 69.85 69.75 69.65 ... 30.35 30.25 30.15 30.05\n", + " * level (level) float32 0.0\n", + " time timedelta64[ns] 2 days\n", + " * lon (lon) float32 -24.95 -24.85 -24.75 -24.65 ... 44.75 44.85 44.95\n", "Attributes:\n", - " cell_measures: area: areacella\n", - " cell_methods: area: time: mean\n", - " comment: near-surface (usually, 2 meter) air temperature\n", - " description: near-surface (usually, 2 meter) air temperature\n", - " frequency: mon\n", - " id: tas\n", - " long_name: Near-Surface Air Temperature\n", - " mipTable: Amon\n", - " out_name: tas\n", - " prov: Amon ((isd.003))\n", - " realm: atmos\n", - " standard_name: air_temperature\n", - " time: time\n", - " time_label: time-mean\n", - " time_title: Temporal mean\n", - " title: Near-Surface Air Temperature\n", - " type: real\n", - " units: K\n", - " variable_id: tas" + " species: PM2.5 Aerosol\n", + " units: µg/m3\n", + " value: hourly values\n", + " standard_name: mass_concentration_of_pm2p5_ambient_aerosol_in_air
<xarray.DataArray 'tas' (member_id: 1, time: 1, lat: 192, lon: 288)>\n", - "array([[[[245.71487, ..., 245.71481],\n", - " ...,\n", - " [246.21568, ..., 246.21373]]]], dtype=float32)\n", + "" ], "text/plain": [ - "<xarray.Dataset>\n", + "Dimensions: (level: 1, time: 97)\n", "Coordinates:\n", - " * lat (lat) float64 -90.0 -89.06 -88.12 -87.17 ... 88.12 89.06 90.0\n", - " * lon (lon) float64 0.0 1.25 2.5 3.75 5.0 ... 355.0 356.2 357.5 358.8\n", - " * time (time) object 2010-01-15 12:00:00\n", - " * member_id (member_id) object 'r1i1p1f1'\n", + " lat float32 46.35\n", + " * level (level) float32 0.0\n", + " * time (time) timedelta64[ns] 00:00:00 01:00:00 ... 4 days 00:00:00\n", + " lon float32 8.75\n", + "Data variables:\n", + " pm2p5_conc (time, level) float32 5.469 5.207 5.003 ... 4.19 4.119 4.017\n", "Attributes:\n", - " cell_measures: area: areacella\n", - " cell_methods: area: time: mean\n", - " comment: near-surface (usually, 2 meter) air temperature\n", - " description: near-surface (usually, 2 meter) air temperature\n", - " frequency: mon\n", - " id: tas\n", - " long_name: Near-Surface Air Temperature\n", - " mipTable: Amon\n", - " out_name: tas\n", - " prov: Amon ((isd.003))\n", - " realm: atmos\n", - " standard_name: air_temperature\n", - " time: time\n", - " time_label: time-mean\n", - " time_title: Temporal mean\n", - " title: Near-Surface Air Temperature\n", - " type: real\n", - " units: K\n", - " variable_id: tas" + " title: PM25 Air Pollutant FORECAST at the Surface\n", + " institution: Data produced by Meteo France\n", + " source: Data from ENSEMBLE model\n", + " history: Model ENSEMBLE FORECAST\n", + " FORECAST: Europe, 20211222+[0H_96H]\n", + " summary: ENSEMBLE model hourly FORECAST of PM25 concentration at the...\n", + " project: MACC-RAQ (http://macc-raq.gmes-atmosphere.eu)
<xarray.Dataset>\n", - "Dimensions: (nbnd: 2, member_id: 1, time: 1980)\n", + "Dimensions: (level: 1, time: 97)\n", "Coordinates:\n", - " lat float64 46.65\n", - " lat_bnds (nbnd) float32 46.18 47.12\n", - " lon float64 8.75\n", - " lon_bnds (nbnd) float32 8.125 9.375\n", - " * time (time) object 1850-01-15 12:00:00 ... 2014-12-15 12:00:00\n", - " time_bnds (time, nbnd) object 1850-01-01 00:00:00 ... 2015-01-01 00:00:00\n", - " * member_id (member_id) object 'r1i1p1f1'\n", - "Dimensions without coordinates: nbnd\n", + " lat float32 59.95\n", + " * level (level) float32 0.0\n", + " * time (time) timedelta64[ns] 00:00:00 01:00:00 ... 4 days 00:00:00\n", + " lon float32 -14.95\n", "Data variables:\n", - " tas (member_id, time) float32 272.7 273.8 274.2 ... 278.1 275.4 272.2\n", - "Attributes: (12/50)\n", - " Conventions: CF-1.7 CMIP-6.2\n", - " activity_id: CMIP\n", - " branch_method: standard\n", - " branch_time_in_child: 674885.0\n", - " branch_time_in_parent: 219000.0\n", - " case_id: 15\n", - " ... ...\n", - " variant_label: r1i1p1f1\n", - " status: 2019-10-25;created;by nhn2@columbia.edu\n", - " netcdf_tracking_ids: hdl:21.14100/d9a7225a-49c3-4470-b7ab-a8180926f839\n", - " version_id: v20190308\n", - " intake_esm_varname: tas\n", - " intake_esm_dataset_key: CMIP.NCAR.CESM2.historical.Amon.gn" + " pm2p5_conc (time, level) float32 2.397 2.281 2.223 ... 2.175 2.192 2.274\n", + "Attributes:\n", + " title: PM25 Air Pollutant FORECAST at the Surface\n", + " institution: Data produced by Meteo France\n", + " source: Data from ENSEMBLE model\n", + " history: Model ENSEMBLE FORECAST\n", + " FORECAST: Europe, 20211222+[0H_96H]\n", + " summary: ENSEMBLE model hourly FORECAST of PM25 concentration at the...\n", + " project: MACC-RAQ (http://macc-raq.gmes-atmosphere.eu) " ], "text/plain": [ "
<xarray.Dataset>\n", - "Dimensions: (nbnd: 2, member_id: 1, time: 1980)\n", - "Coordinates:\n", - " lat float64 4.241\n", - " lat_bnds (nbnd) float32 3.77 4.712\n", - " lon float64 125.0\n", - " lon_bnds (nbnd) float32 124.4 125.6\n", - " * time (time) object 1850-01-15 12:00:00 ... 2014-12-15 12:00:00\n", - " time_bnds (time, nbnd) object 1850-01-01 00:00:00 ... 2015-01-01 00:00:00\n", - " * member_id (member_id) object 'r1i1p1f1'\n", - "Dimensions without coordinates: nbnd\n", - "Data variables:\n", - " tas (member_id, time) float32 300.4 299.7 299.9 ... 302.3 302.0 301.8\n", - "Attributes: (12/50)\n", - " Conventions: CF-1.7 CMIP-6.2\n", - " activity_id: CMIP\n", - " branch_method: standard\n", - " branch_time_in_child: 674885.0\n", - " branch_time_in_parent: 219000.0\n", - " case_id: 15\n", - " ... ...\n", - " variant_label: r1i1p1f1\n", - " status: 2019-10-25;created;by nhn2@columbia.edu\n", - " netcdf_tracking_ids: hdl:21.14100/d9a7225a-49c3-4470-b7ab-a8180926f839\n", - " version_id: v20190308\n", - " intake_esm_varname: tas\n", - " intake_esm_dataset_key: CMIP.NCAR.CESM2.historical.Amon.gn
array([69.95, 69.85, 69.75, 69.65, 69.55, 69.45, 69.35, 69.25, 69.15, 69.05,\n", + " 68.95, 68.85, 68.75, 68.65, 68.55, 68.45, 68.35, 68.25, 68.15, 68.05,\n", + " 67.95, 67.85, 67.75, 67.65, 67.55, 67.45, 67.35, 67.25, 67.15, 67.05,\n", + " 66.95, 66.85, 66.75, 66.65, 66.55, 66.45, 66.35, 66.25, 66.15, 66.05,\n", + " 65.95, 65.85, 65.75, 65.65, 65.55, 65.45, 65.35, 65.25, 65.15, 65.05,\n", + " 64.95, 64.85, 64.75, 64.65, 64.55, 64.45, 64.35, 64.25, 64.15, 64.05,\n", + " 63.95, 63.85, 63.75, 63.65, 63.55, 63.45, 63.35, 63.25, 63.15, 63.05,\n", + " 62.95, 62.85, 62.75, 62.65, 62.55, 62.45, 62.35, 62.25, 62.15, 62.05,\n", + " 61.95, 61.85, 61.75, 61.65, 61.55, 61.45, 61.35, 61.25, 61.15, 61.05,\n", + " 60.95, 60.85, 60.75, 60.65, 60.55, 60.45, 60.35, 60.25, 60.15, 60.05,\n", + " 59.95, 59.85, 59.75, 59.65, 59.55, 59.45, 59.35, 59.25, 59.15, 59.05,\n", + " 58.95, 58.85, 58.75, 58.65, 58.55, 58.45, 58.35, 58.25, 58.15, 58.05,\n", + " 57.95, 57.85, 57.75, 57.65, 57.55, 57.45, 57.35, 57.25, 57.15, 57.05,\n", + " 56.95, 56.85, 56.75, 56.65, 56.55, 56.45, 56.35, 56.25, 56.15, 56.05,\n", + " 55.95, 55.85, 55.75, 55.65, 55.55, 55.45, 55.35, 55.25, 55.15, 55.05,\n", + " 54.95, 54.85, 54.75, 54.65, 54.55], dtype=float32)
array([0.], dtype=float32)
array([ 0, 3600000000000, 7200000000000, 10800000000000,\n", + " 14400000000000, 18000000000000, 21600000000000, 25200000000000,\n", + " 28800000000000, 32400000000000, 36000000000000, 39600000000000,\n", + " 43200000000000, 46800000000000, 50400000000000, 54000000000000,\n", + " 57600000000000, 61200000000000, 64800000000000, 68400000000000,\n", + " 72000000000000, 75600000000000, 79200000000000, 82800000000000,\n", + " 86400000000000, 90000000000000, 93600000000000, 97200000000000,\n", + " 100800000000000, 104400000000000, 108000000000000, 111600000000000,\n", + " 115200000000000, 118800000000000, 122400000000000, 126000000000000,\n", + " 129600000000000, 133200000000000, 136800000000000, 140400000000000,\n", + " 144000000000000, 147600000000000, 151200000000000, 154800000000000,\n", + " 158400000000000, 162000000000000, 165600000000000, 169200000000000,\n", + " 172800000000000, 176400000000000, 180000000000000, 183600000000000,\n", + " 187200000000000, 190800000000000, 194400000000000, 198000000000000,\n", + " 201600000000000, 205200000000000, 208800000000000, 212400000000000,\n", + " 216000000000000, 219600000000000, 223200000000000, 226800000000000,\n", + " 230400000000000, 234000000000000, 237600000000000, 241200000000000,\n", + " 244800000000000, 248400000000000, 252000000000000, 255600000000000,\n", + " 259200000000000, 262800000000000, 266400000000000, 270000000000000,\n", + " 273600000000000, 277200000000000, 280800000000000, 284400000000000,\n", + " 288000000000000, 291600000000000, 295200000000000, 298800000000000,\n", + " 302400000000000, 306000000000000, 309600000000000, 313200000000000,\n", + " 316800000000000, 320400000000000, 324000000000000, 327600000000000,\n", + " 331200000000000, 334800000000000, 338400000000000, 342000000000000,\n", + " 345600000000000], dtype='timedelta64[ns]')
array([-2.450012, -2.349976, -2.25 , ..., 42.25 , 42.350006, 42.449997],\n", + " dtype=float32)
PandasIndex(Index([ 69.94999694824219, 69.8499984741211, 69.75,\n", + " 69.6500015258789, 69.55000305175781, 69.44999694824219,\n", + " 69.3499984741211, 69.25, 69.1500015258789,\n", + " 69.05000305175781,\n", + " ...\n", + " 55.45000076293945, 55.349998474121094, 55.25,\n", + " 55.150001525878906, 55.04999923706055, 54.95000076293945,\n", + " 54.849998474121094, 54.75, 54.650001525878906,\n", + " 54.54999923706055],\n", + " dtype='float32', name='lat', length=155))
PandasIndex(Index([0.0], dtype='float32', name='level'))
PandasIndex(TimedeltaIndex(['0 days 00:00:00', '0 days 01:00:00', '0 days 02:00:00',\n", + " '0 days 03:00:00', '0 days 04:00:00', '0 days 05:00:00',\n", + " '0 days 06:00:00', '0 days 07:00:00', '0 days 08:00:00',\n", + " '0 days 09:00:00', '0 days 10:00:00', '0 days 11:00:00',\n", + " '0 days 12:00:00', '0 days 13:00:00', '0 days 14:00:00',\n", + " '0 days 15:00:00', '0 days 16:00:00', '0 days 17:00:00',\n", + " '0 days 18:00:00', '0 days 19:00:00', '0 days 20:00:00',\n", + " '0 days 21:00:00', '0 days 22:00:00', '0 days 23:00:00',\n", + " '1 days 00:00:00', '1 days 01:00:00', '1 days 02:00:00',\n", + " '1 days 03:00:00', '1 days 04:00:00', '1 days 05:00:00',\n", + " '1 days 06:00:00', '1 days 07:00:00', '1 days 08:00:00',\n", + " '1 days 09:00:00', '1 days 10:00:00', '1 days 11:00:00',\n", + " '1 days 12:00:00', '1 days 13:00:00', '1 days 14:00:00',\n", + " '1 days 15:00:00', '1 days 16:00:00', '1 days 17:00:00',\n", + " '1 days 18:00:00', '1 days 19:00:00', '1 days 20:00:00',\n", + " '1 days 21:00:00', '1 days 22:00:00', '1 days 23:00:00',\n", + " '2 days 00:00:00', '2 days 01:00:00', '2 days 02:00:00',\n", + " '2 days 03:00:00', '2 days 04:00:00', '2 days 05:00:00',\n", + " '2 days 06:00:00', '2 days 07:00:00', '2 days 08:00:00',\n", + " '2 days 09:00:00', '2 days 10:00:00', '2 days 11:00:00',\n", + " '2 days 12:00:00', '2 days 13:00:00', '2 days 14:00:00',\n", + " '2 days 15:00:00', '2 days 16:00:00', '2 days 17:00:00',\n", + " '2 days 18:00:00', '2 days 19:00:00', '2 days 20:00:00',\n", + " '2 days 21:00:00', '2 days 22:00:00', '2 days 23:00:00',\n", + " '3 days 00:00:00', '3 days 01:00:00', '3 days 02:00:00',\n", + " '3 days 03:00:00', '3 days 04:00:00', '3 days 05:00:00',\n", + " '3 days 06:00:00', '3 days 07:00:00', '3 days 08:00:00',\n", + " '3 days 09:00:00', '3 days 10:00:00', '3 days 11:00:00',\n", + " '3 days 12:00:00', '3 days 13:00:00', '3 days 14:00:00',\n", + " '3 days 15:00:00', '3 days 16:00:00', '3 days 17:00:00',\n", + " '3 days 18:00:00', '3 days 19:00:00', '3 days 20:00:00',\n", + " '3 days 21:00:00', '3 days 22:00:00', '3 days 23:00:00',\n", + " '4 days 00:00:00'],\n", + " dtype='timedelta64[ns]', name='time', freq=None))
PandasIndex(Index([ -2.45001220703125, -2.3499755859375, -2.25,\n", + " -2.1500244140625, -2.04998779296875, -1.95001220703125,\n", + " -1.8499755859375, -1.75, -1.6500244140625,\n", + " -1.54998779296875,\n", + " ...\n", + " 41.55000305175781, 41.649993896484375, 41.75,\n", + " 41.850006103515625, 41.94999694824219, 42.05000305175781,\n", + " 42.149993896484375, 42.25, 42.350006103515625,\n", + " 42.44999694824219],\n", + " dtype='float32', name='lon', length=450))
<xarray.DataArray 'tas' (member_id: 1, time: 1980, lat: 18, lon: 35)>\n", - "array([[[[277.14154, ..., 255.24866],\n", + "<xarray.DataArray 'pm2p5_conc' (time: 97, level: 1, lat: 155, lon: 450)>\n", + "array([[[[ 855.6641 , 852.48096, 833.32465, ..., 879.8937 ,\n", + " 902.90106, 865.896 ],\n", + " [ 885.7343 , 891.02075, 878.97 , ..., 971.65314,\n", + " 1061.4657 , 1015.5647 ],\n", + " [ 893.89136, 828.7204 , 800.25604, ..., 1032.3335 ,\n", + " 990.852 , 982.55286],\n", " ...,\n", - " [271.17374, ..., 259.73807]],\n", - "\n", - " ...,\n", - "\n", - " [[283.0737 , ..., 269.6892 ],\n", + " [10043.194 , 10672.565 , 11238.228 , ..., 7507.282 ,\n", + " 7414.5137 , 7444.7827 ],\n", + " [13826.692 , 14108.494 , 14563.739 , ..., 7519.0063 ,\n", + " 7214.311 , 7085.4043 ],\n", + " [14878.069 , 14873.55 , 14858.358 , ..., 7069.7866 ,\n", + " 7237.475 , 7302.5605 ]]],\n", + "\n", + "\n", + " [[[ 868.89185, 899.5873 , 930.3964 , ..., 936.80554,\n", + " 839.1485 , 772.7412 ],\n", + " [ 858.319 , 891.28815, 934.2618 , ..., 772.7412 ,\n", + " 733.67554, 716.23883],\n", + " [ 873.09827, 925.96265, 940.6993 , ..., 729.8386 ,\n", + "...\n", + " 6790.51 , 6748.986 ],\n", + " [ 1623.8978 , 1855.4778 , 2173.332 , ..., 5969.662 ,\n", + " 6049.47 , 6131.5386 ],\n", + " [ 2051.7156 , 2147.7383 , 1941.084 , ..., 5981.9263 ,\n", + " 5846.8096 , 5826.9146 ]]],\n", + "\n", + "\n", + " [[[ 1362.0598 , 1371.5527 , 1349.8243 , ..., 2480.518 ,\n", + " 2496.761 , 2526.6606 ],\n", + " [ 1398.2976 , 1351.3591 , 1315.5477 , ..., 2294.1143 ,\n", + " 2327.9363 , 2370.2136 ],\n", + " [ 1287.1544 , 1304.108 , 1268.4387 , ..., 2281.4312 ,\n", + " 2239.3528 , 2212.153 ],\n", " ...,\n", - " [272.89108, ..., 271.8542 ]]]], dtype=float32)\n", + " [ 1829.8173 , 1892.3165 , 1997.889 , ..., 6986.432 ,\n", + " 6914.4253 , 6935.685 ],\n", + " [ 1650.3838 , 1853.0875 , 2025.124 , ..., 6504.3003 ,\n", + " 6350.027 , 6176.655 ],\n", + " [ 1954.8159 , 2019.9797 , 1916.7521 , ..., 5647.4995 ,\n", + " 5836.5464 , 6040.9556 ]]]], dtype=float32)\n", "Coordinates:\n", - " * lat (lat) float64 55.13 56.07 57.02 57.96 ... 68.32 69.27 70.21 71.15\n", - " * lon (lon) float64 0.0 1.25 2.5 3.75 5.0 ... 38.75 40.0 41.25 42.5\n", - " * time (time) object 1850-01-15 12:00:00 ... 2014-12-15 12:00:00\n", - " * member_id (member_id) object 'r1i1p1f1'\n", - "Attributes:\n", - " cell_measures: area: areacella\n", - " cell_methods: area: time: mean\n", - " comment: near-surface (usually, 2 meter) air temperature\n", - " description: near-surface (usually, 2 meter) air temperature\n", - " frequency: mon\n", - " id: tas\n", - " long_name: Near-Surface Air Temperature\n", - " mipTable: Amon\n", - " out_name: tas\n", - " prov: Amon ((isd.003))\n", - " realm: atmos\n", - " standard_name: air_temperature\n", - " time: time\n", - " time_label: time-mean\n", - " time_title: Temporal mean\n", - " title: Near-Surface Air Temperature\n", - " type: real\n", - " units: K\n", - " variable_id: tas
array([69.95, 69.85, 69.75, 69.65, 69.55, 69.45, 69.35, 69.25, 69.15, 69.05,\n", + " 68.95, 68.85, 68.75, 68.65, 68.55, 68.45, 68.35, 68.25, 68.15, 68.05,\n", + " 67.95, 67.85, 67.75, 67.65, 67.55, 67.45, 67.35, 67.25, 67.15, 67.05,\n", + " 66.95, 66.85, 66.75, 66.65, 66.55, 66.45, 66.35, 66.25, 66.15, 66.05,\n", + " 65.95, 65.85, 65.75, 65.65, 65.55, 65.45, 65.35, 65.25, 65.15, 65.05,\n", + " 64.95, 64.85, 64.75, 64.65, 64.55, 64.45, 64.35, 64.25, 64.15, 64.05,\n", + " 63.95, 63.85, 63.75, 63.65, 63.55, 63.45, 63.35, 63.25, 63.15, 63.05,\n", + " 62.95, 62.85, 62.75, 62.65, 62.55, 62.45, 62.35, 62.25, 62.15, 62.05,\n", + " 61.95, 61.85, 61.75, 61.65, 61.55, 61.45, 61.35, 61.25, 61.15, 61.05,\n", + " 60.95, 60.85, 60.75, 60.65, 60.55, 60.45, 60.35, 60.25, 60.15, 60.05,\n", + " 59.95, 59.85, 59.75, 59.65, 59.55, 59.45, 59.35, 59.25, 59.15, 59.05,\n", + " 58.95, 58.85, 58.75, 58.65, 58.55, 58.45, 58.35, 58.25, 58.15, 58.05,\n", + " 57.95, 57.85, 57.75, 57.65, 57.55, 57.45, 57.35, 57.25, 57.15, 57.05,\n", + " 56.95, 56.85, 56.75, 56.65, 56.55, 56.45, 56.35, 56.25, 56.15, 56.05,\n", + " 55.95, 55.85, 55.75, 55.65, 55.55, 55.45, 55.35, 55.25, 55.15, 55.05,\n", + " 54.95, 54.85, 54.75, 54.65, 54.55], dtype=float32)
array([0.], dtype=float32)
array([ 0, 3600000000000, 7200000000000, 10800000000000,\n", + " 14400000000000, 18000000000000, 21600000000000, 25200000000000,\n", + " 28800000000000, 32400000000000, 36000000000000, 39600000000000,\n", + " 43200000000000, 46800000000000, 50400000000000, 54000000000000,\n", + " 57600000000000, 61200000000000, 64800000000000, 68400000000000,\n", + " 72000000000000, 75600000000000, 79200000000000, 82800000000000,\n", + " 86400000000000, 90000000000000, 93600000000000, 97200000000000,\n", + " 100800000000000, 104400000000000, 108000000000000, 111600000000000,\n", + " 115200000000000, 118800000000000, 122400000000000, 126000000000000,\n", + " 129600000000000, 133200000000000, 136800000000000, 140400000000000,\n", + " 144000000000000, 147600000000000, 151200000000000, 154800000000000,\n", + " 158400000000000, 162000000000000, 165600000000000, 169200000000000,\n", + " 172800000000000, 176400000000000, 180000000000000, 183600000000000,\n", + " 187200000000000, 190800000000000, 194400000000000, 198000000000000,\n", + " 201600000000000, 205200000000000, 208800000000000, 212400000000000,\n", + " 216000000000000, 219600000000000, 223200000000000, 226800000000000,\n", + " 230400000000000, 234000000000000, 237600000000000, 241200000000000,\n", + " 244800000000000, 248400000000000, 252000000000000, 255600000000000,\n", + " 259200000000000, 262800000000000, 266400000000000, 270000000000000,\n", + " 273600000000000, 277200000000000, 280800000000000, 284400000000000,\n", + " 288000000000000, 291600000000000, 295200000000000, 298800000000000,\n", + " 302400000000000, 306000000000000, 309600000000000, 313200000000000,\n", + " 316800000000000, 320400000000000, 324000000000000, 327600000000000,\n", + " 331200000000000, 334800000000000, 338400000000000, 342000000000000,\n", + " 345600000000000], dtype='timedelta64[ns]')
array([-2.450012, -2.349976, -2.25 , ..., 42.25 , 42.350006, 42.449997],\n", + " dtype=float32)
PandasIndex(Index([ 69.94999694824219, 69.8499984741211, 69.75,\n", + " 69.6500015258789, 69.55000305175781, 69.44999694824219,\n", + " 69.3499984741211, 69.25, 69.1500015258789,\n", + " 69.05000305175781,\n", + " ...\n", + " 55.45000076293945, 55.349998474121094, 55.25,\n", + " 55.150001525878906, 55.04999923706055, 54.95000076293945,\n", + " 54.849998474121094, 54.75, 54.650001525878906,\n", + " 54.54999923706055],\n", + " dtype='float32', name='lat', length=155))
PandasIndex(Index([0.0], dtype='float32', name='level'))
PandasIndex(TimedeltaIndex(['0 days 00:00:00', '0 days 01:00:00', '0 days 02:00:00',\n", + " '0 days 03:00:00', '0 days 04:00:00', '0 days 05:00:00',\n", + " '0 days 06:00:00', '0 days 07:00:00', '0 days 08:00:00',\n", + " '0 days 09:00:00', '0 days 10:00:00', '0 days 11:00:00',\n", + " '0 days 12:00:00', '0 days 13:00:00', '0 days 14:00:00',\n", + " '0 days 15:00:00', '0 days 16:00:00', '0 days 17:00:00',\n", + " '0 days 18:00:00', '0 days 19:00:00', '0 days 20:00:00',\n", + " '0 days 21:00:00', '0 days 22:00:00', '0 days 23:00:00',\n", + " '1 days 00:00:00', '1 days 01:00:00', '1 days 02:00:00',\n", + " '1 days 03:00:00', '1 days 04:00:00', '1 days 05:00:00',\n", + " '1 days 06:00:00', '1 days 07:00:00', '1 days 08:00:00',\n", + " '1 days 09:00:00', '1 days 10:00:00', '1 days 11:00:00',\n", + " '1 days 12:00:00', '1 days 13:00:00', '1 days 14:00:00',\n", + " '1 days 15:00:00', '1 days 16:00:00', '1 days 17:00:00',\n", + " '1 days 18:00:00', '1 days 19:00:00', '1 days 20:00:00',\n", + " '1 days 21:00:00', '1 days 22:00:00', '1 days 23:00:00',\n", + " '2 days 00:00:00', '2 days 01:00:00', '2 days 02:00:00',\n", + " '2 days 03:00:00', '2 days 04:00:00', '2 days 05:00:00',\n", + " '2 days 06:00:00', '2 days 07:00:00', '2 days 08:00:00',\n", + " '2 days 09:00:00', '2 days 10:00:00', '2 days 11:00:00',\n", + " '2 days 12:00:00', '2 days 13:00:00', '2 days 14:00:00',\n", + " '2 days 15:00:00', '2 days 16:00:00', '2 days 17:00:00',\n", + " '2 days 18:00:00', '2 days 19:00:00', '2 days 20:00:00',\n", + " '2 days 21:00:00', '2 days 22:00:00', '2 days 23:00:00',\n", + " '3 days 00:00:00', '3 days 01:00:00', '3 days 02:00:00',\n", + " '3 days 03:00:00', '3 days 04:00:00', '3 days 05:00:00',\n", + " '3 days 06:00:00', '3 days 07:00:00', '3 days 08:00:00',\n", + " '3 days 09:00:00', '3 days 10:00:00', '3 days 11:00:00',\n", + " '3 days 12:00:00', '3 days 13:00:00', '3 days 14:00:00',\n", + " '3 days 15:00:00', '3 days 16:00:00', '3 days 17:00:00',\n", + " '3 days 18:00:00', '3 days 19:00:00', '3 days 20:00:00',\n", + " '3 days 21:00:00', '3 days 22:00:00', '3 days 23:00:00',\n", + " '4 days 00:00:00'],\n", + " dtype='timedelta64[ns]', name='time', freq=None))
PandasIndex(Index([ -2.45001220703125, -2.3499755859375, -2.25,\n", + " -2.1500244140625, -2.04998779296875, -1.95001220703125,\n", + " -1.8499755859375, -1.75, -1.6500244140625,\n", + " -1.54998779296875,\n", + " ...\n", + " 41.55000305175781, 41.649993896484375, 41.75,\n", + " 41.850006103515625, 41.94999694824219, 42.05000305175781,\n", + " 42.149993896484375, 42.25, 42.350006103515625,\n", + " 42.44999694824219],\n", + " dtype='float32', name='lon', length=450))
<xarray.DataArray 'tas' (member_id: 1, time: 1980, lat: 18, lon: 35)>\n", - "array([[[[ 3.9915466 , 3.9372864 , 3.63208 , ..., -17.947784 ,\n", - " -18.09903 , -17.901337 ],\n", - " [ 3.9906616 , 3.8585205 , 3.4658813 , ..., -19.287094 ,\n", - " -19.638443 , -19.815659 ],\n", - " [ 3.8772888 , 3.693756 , 3.3085938 , ..., -20.736603 ,\n", - " -21.195053 , -21.357864 ],\n", + "<xarray.DataArray 'pm2p5_conc' (time: 97, level: 1, lat: 155, lon: 450)>\n", + "array([[[[ 855.6641 , 852.48096, 833.32465, ..., 879.8937 ,\n", + " 902.90106, 865.896 ],\n", + " [ 885.7343 , 891.02075, 878.97 , ..., 971.65314,\n", + " 1061.4657 , 1015.5647 ],\n", + " [ 893.89136, 828.7204 , 800.25604, ..., 1032.3335 ,\n", + " 990.852 , 982.55286],\n", " ...,\n", - " [ -0.11868286, 0.1071167 , 0.3340149 , ..., -11.844849 ,\n", - " -12.702545 , -14.013733 ],\n", - " [ -0.9946594 , -0.6482544 , -0.4326172 , ..., -10.813141 ,\n", - " -12.086548 , -13.261871 ],\n", - " [ -1.9762573 , -1.6532898 , -1.2767944 , ..., -10.560425 ,\n", - " -11.828857 , -13.411926 ]],\n", - "\n", - " [[ 6.5706787 , 6.261444 , 5.76593 , ..., -10.512054 ,\n", - " -10.819702 , -11.010132 ],\n", - " [ 6.6134033 , 6.326996 , 5.8705444 , ..., -10.703583 ,\n", - " -10.812073 , -11.228058 ],\n", - " [ 6.63974 , 6.3713074 , 5.99823 , ..., -10.811035 ,\n", - " -11.002502 , -11.324585 ],\n", + " [10043.194 , 10672.565 , 11238.228 , ..., 7507.282 ,\n", + " 7414.5137 , 7444.7827 ],\n", + " [13826.692 , 14108.494 , 14563.739 , ..., 7519.0063 ,\n", + " 7214.311 , 7085.4043 ],\n", + " [14878.069 , 14873.55 , 14858.358 , ..., 7069.7866 ,\n", + " 7237.475 , 7302.5605 ]]],\n", + "\n", + "\n", + " [[[ 868.89185, 899.5873 , 930.3964 , ..., 936.80554,\n", + " 839.1485 , 772.7412 ],\n", + " [ 858.319 , 891.28815, 934.2618 , ..., 772.7412 ,\n", + " 733.67554, 716.23883],\n", + " [ 873.09827, 925.96265, 940.6993 , ..., 729.8386 ,\n", "...\n", - " [ 4.5212708 , 4.881195 , 5.349457 , ..., -0.38009644,\n", - " -0.14691162, -0.15356445],\n", - " [ 3.3636475 , 3.8493652 , 4.25 , ..., 0.430542 ,\n", - " 0.32699585, 0.23458862],\n", - " [ 2.291626 , 2.7035522 , 3.230133 , ..., 0.6605835 ,\n", - " 0.5307617 , 0.37503052]],\n", - "\n", - " [[ 9.923706 , 10.213959 , 10.246948 , ..., -2.7084656 ,\n", - " -3.117157 , -3.460785 ],\n", - " [ 10.048218 , 10.249817 , 10.259094 , ..., -2.855835 ,\n", - " -3.1342468 , -3.5550232 ],\n", - " [ 9.972595 , 10.138733 , 10.1493225 , ..., -3.176941 ,\n", - " -3.5891113 , -3.8943787 ],\n", + " 6790.51 , 6748.986 ],\n", + " [ 1623.8978 , 1855.4778 , 2173.332 , ..., 5969.662 ,\n", + " 6049.47 , 6131.5386 ],\n", + " [ 2051.7156 , 2147.7383 , 1941.084 , ..., 5981.9263 ,\n", + " 5846.8096 , 5826.9146 ]]],\n", + "\n", + "\n", + " [[[ 1362.0598 , 1371.5527 , 1349.8243 , ..., 2480.518 ,\n", + " 2496.761 , 2526.6606 ],\n", + " [ 1398.2976 , 1351.3591 , 1315.5477 , ..., 2294.1143 ,\n", + " 2327.9363 , 2370.2136 ],\n", + " [ 1287.1544 , 1304.108 , 1268.4387 , ..., 2281.4312 ,\n", + " 2239.3528 , 2212.153 ],\n", " ...,\n", - " [ 2.2316895 , 2.5167236 , 2.8789368 , ..., -1.0142822 ,\n", - " -1.0447388 , -1.3068542 ],\n", - " [ 0.93966675, 1.3753662 , 1.7021484 , ..., -0.4562683 ,\n", - " -0.82281494, -1.1673279 ],\n", - " [ -0.25891113, 0.1281128 , 0.6257019 , ..., -0.47140503,\n", - " -0.86395264, -1.2958069 ]]]], dtype=float32)\n", + " [ 1829.8173 , 1892.3165 , 1997.889 , ..., 6986.432 ,\n", + " 6914.4253 , 6935.685 ],\n", + " [ 1650.3838 , 1853.0875 , 2025.124 , ..., 6504.3003 ,\n", + " 6350.027 , 6176.655 ],\n", + " [ 1954.8159 , 2019.9797 , 1916.7521 , ..., 5647.4995 ,\n", + " 5836.5464 , 6040.9556 ]]]], dtype=float32)\n", "Coordinates:\n", - " * lat (lat) float64 55.13 56.07 57.02 57.96 ... 68.32 69.27 70.21 71.15\n", - " * lon (lon) float64 0.0 1.25 2.5 3.75 5.0 ... 38.75 40.0 41.25 42.5\n", - " * time (time) object 1850-01-15 12:00:00 ... 2014-12-15 12:00:00\n", - " * member_id (member_id) object 'r1i1p1f1'
array([69.95, 69.85, 69.75, 69.65, 69.55, 69.45, 69.35, 69.25, 69.15, 69.05,\n", + " 68.95, 68.85, 68.75, 68.65, 68.55, 68.45, 68.35, 68.25, 68.15, 68.05,\n", + " 67.95, 67.85, 67.75, 67.65, 67.55, 67.45, 67.35, 67.25, 67.15, 67.05,\n", + " 66.95, 66.85, 66.75, 66.65, 66.55, 66.45, 66.35, 66.25, 66.15, 66.05,\n", + " 65.95, 65.85, 65.75, 65.65, 65.55, 65.45, 65.35, 65.25, 65.15, 65.05,\n", + " 64.95, 64.85, 64.75, 64.65, 64.55, 64.45, 64.35, 64.25, 64.15, 64.05,\n", + " 63.95, 63.85, 63.75, 63.65, 63.55, 63.45, 63.35, 63.25, 63.15, 63.05,\n", + " 62.95, 62.85, 62.75, 62.65, 62.55, 62.45, 62.35, 62.25, 62.15, 62.05,\n", + " 61.95, 61.85, 61.75, 61.65, 61.55, 61.45, 61.35, 61.25, 61.15, 61.05,\n", + " 60.95, 60.85, 60.75, 60.65, 60.55, 60.45, 60.35, 60.25, 60.15, 60.05,\n", + " 59.95, 59.85, 59.75, 59.65, 59.55, 59.45, 59.35, 59.25, 59.15, 59.05,\n", + " 58.95, 58.85, 58.75, 58.65, 58.55, 58.45, 58.35, 58.25, 58.15, 58.05,\n", + " 57.95, 57.85, 57.75, 57.65, 57.55, 57.45, 57.35, 57.25, 57.15, 57.05,\n", + " 56.95, 56.85, 56.75, 56.65, 56.55, 56.45, 56.35, 56.25, 56.15, 56.05,\n", + " 55.95, 55.85, 55.75, 55.65, 55.55, 55.45, 55.35, 55.25, 55.15, 55.05,\n", + " 54.95, 54.85, 54.75, 54.65, 54.55], dtype=float32)
array([0.], dtype=float32)
array([ 0, 3600000000000, 7200000000000, 10800000000000,\n", + " 14400000000000, 18000000000000, 21600000000000, 25200000000000,\n", + " 28800000000000, 32400000000000, 36000000000000, 39600000000000,\n", + " 43200000000000, 46800000000000, 50400000000000, 54000000000000,\n", + " 57600000000000, 61200000000000, 64800000000000, 68400000000000,\n", + " 72000000000000, 75600000000000, 79200000000000, 82800000000000,\n", + " 86400000000000, 90000000000000, 93600000000000, 97200000000000,\n", + " 100800000000000, 104400000000000, 108000000000000, 111600000000000,\n", + " 115200000000000, 118800000000000, 122400000000000, 126000000000000,\n", + " 129600000000000, 133200000000000, 136800000000000, 140400000000000,\n", + " 144000000000000, 147600000000000, 151200000000000, 154800000000000,\n", + " 158400000000000, 162000000000000, 165600000000000, 169200000000000,\n", + " 172800000000000, 176400000000000, 180000000000000, 183600000000000,\n", + " 187200000000000, 190800000000000, 194400000000000, 198000000000000,\n", + " 201600000000000, 205200000000000, 208800000000000, 212400000000000,\n", + " 216000000000000, 219600000000000, 223200000000000, 226800000000000,\n", + " 230400000000000, 234000000000000, 237600000000000, 241200000000000,\n", + " 244800000000000, 248400000000000, 252000000000000, 255600000000000,\n", + " 259200000000000, 262800000000000, 266400000000000, 270000000000000,\n", + " 273600000000000, 277200000000000, 280800000000000, 284400000000000,\n", + " 288000000000000, 291600000000000, 295200000000000, 298800000000000,\n", + " 302400000000000, 306000000000000, 309600000000000, 313200000000000,\n", + " 316800000000000, 320400000000000, 324000000000000, 327600000000000,\n", + " 331200000000000, 334800000000000, 338400000000000, 342000000000000,\n", + " 345600000000000], dtype='timedelta64[ns]')
array([-2.450012, -2.349976, -2.25 , ..., 42.25 , 42.350006, 42.449997],\n", + " dtype=float32)
PandasIndex(Index([ 69.94999694824219, 69.8499984741211, 69.75,\n", + " 69.6500015258789, 69.55000305175781, 69.44999694824219,\n", + " 69.3499984741211, 69.25, 69.1500015258789,\n", + " 69.05000305175781,\n", + " ...\n", + " 55.45000076293945, 55.349998474121094, 55.25,\n", + " 55.150001525878906, 55.04999923706055, 54.95000076293945,\n", + " 54.849998474121094, 54.75, 54.650001525878906,\n", + " 54.54999923706055],\n", + " dtype='float32', name='lat', length=155))
PandasIndex(Index([0.0], dtype='float32', name='level'))
PandasIndex(TimedeltaIndex(['0 days 00:00:00', '0 days 01:00:00', '0 days 02:00:00',\n", + " '0 days 03:00:00', '0 days 04:00:00', '0 days 05:00:00',\n", + " '0 days 06:00:00', '0 days 07:00:00', '0 days 08:00:00',\n", + " '0 days 09:00:00', '0 days 10:00:00', '0 days 11:00:00',\n", + " '0 days 12:00:00', '0 days 13:00:00', '0 days 14:00:00',\n", + " '0 days 15:00:00', '0 days 16:00:00', '0 days 17:00:00',\n", + " '0 days 18:00:00', '0 days 19:00:00', '0 days 20:00:00',\n", + " '0 days 21:00:00', '0 days 22:00:00', '0 days 23:00:00',\n", + " '1 days 00:00:00', '1 days 01:00:00', '1 days 02:00:00',\n", + " '1 days 03:00:00', '1 days 04:00:00', '1 days 05:00:00',\n", + " '1 days 06:00:00', '1 days 07:00:00', '1 days 08:00:00',\n", + " '1 days 09:00:00', '1 days 10:00:00', '1 days 11:00:00',\n", + " '1 days 12:00:00', '1 days 13:00:00', '1 days 14:00:00',\n", + " '1 days 15:00:00', '1 days 16:00:00', '1 days 17:00:00',\n", + " '1 days 18:00:00', '1 days 19:00:00', '1 days 20:00:00',\n", + " '1 days 21:00:00', '1 days 22:00:00', '1 days 23:00:00',\n", + " '2 days 00:00:00', '2 days 01:00:00', '2 days 02:00:00',\n", + " '2 days 03:00:00', '2 days 04:00:00', '2 days 05:00:00',\n", + " '2 days 06:00:00', '2 days 07:00:00', '2 days 08:00:00',\n", + " '2 days 09:00:00', '2 days 10:00:00', '2 days 11:00:00',\n", + " '2 days 12:00:00', '2 days 13:00:00', '2 days 14:00:00',\n", + " '2 days 15:00:00', '2 days 16:00:00', '2 days 17:00:00',\n", + " '2 days 18:00:00', '2 days 19:00:00', '2 days 20:00:00',\n", + " '2 days 21:00:00', '2 days 22:00:00', '2 days 23:00:00',\n", + " '3 days 00:00:00', '3 days 01:00:00', '3 days 02:00:00',\n", + " '3 days 03:00:00', '3 days 04:00:00', '3 days 05:00:00',\n", + " '3 days 06:00:00', '3 days 07:00:00', '3 days 08:00:00',\n", + " '3 days 09:00:00', '3 days 10:00:00', '3 days 11:00:00',\n", + " '3 days 12:00:00', '3 days 13:00:00', '3 days 14:00:00',\n", + " '3 days 15:00:00', '3 days 16:00:00', '3 days 17:00:00',\n", + " '3 days 18:00:00', '3 days 19:00:00', '3 days 20:00:00',\n", + " '3 days 21:00:00', '3 days 22:00:00', '3 days 23:00:00',\n", + " '4 days 00:00:00'],\n", + " dtype='timedelta64[ns]', name='time', freq=None))
PandasIndex(Index([ -2.45001220703125, -2.3499755859375, -2.25,\n", + " -2.1500244140625, -2.04998779296875, -1.95001220703125,\n", + " -1.8499755859375, -1.75, -1.6500244140625,\n", + " -1.54998779296875,\n", + " ...\n", + " 41.55000305175781, 41.649993896484375, 41.75,\n", + " 41.850006103515625, 41.94999694824219, 42.05000305175781,\n", + " 42.149993896484375, 42.25, 42.350006103515625,\n", + " 42.44999694824219],\n", + " dtype='float32', name='lon', length=450))
<xarray.DataArray 'tas' (member_id: 1, time: 1980, lat: 18, lon: 35)>\n", - "array([[[[ 3.9915466 , 3.9372864 , 3.63208 , ..., -17.947784 ,\n", - " -18.09903 , -17.901337 ],\n", - " [ 3.9906616 , 3.8585205 , 3.4658813 , ..., -19.287094 ,\n", - " -19.638443 , -19.815659 ],\n", - " [ 3.8772888 , 3.693756 , 3.3085938 , ..., -20.736603 ,\n", - " -21.195053 , -21.357864 ],\n", + "<xarray.DataArray 'pm2p5_conc' (time: 97, level: 1, lat: 155, lon: 450)>\n", + "array([[[[ 855.6641 , 852.48096, 833.32465, ..., 879.8937 ,\n", + " 902.90106, 865.896 ],\n", + " [ 885.7343 , 891.02075, 878.97 , ..., 971.65314,\n", + " 1061.4657 , 1015.5647 ],\n", + " [ 893.89136, 828.7204 , 800.25604, ..., 1032.3335 ,\n", + " 990.852 , 982.55286],\n", " ...,\n", - " [ -0.11868286, 0.1071167 , 0.3340149 , ..., -11.844849 ,\n", - " -12.702545 , -14.013733 ],\n", - " [ -0.9946594 , -0.6482544 , -0.4326172 , ..., -10.813141 ,\n", - " -12.086548 , -13.261871 ],\n", - " [ -1.9762573 , -1.6532898 , -1.2767944 , ..., -10.560425 ,\n", - " -11.828857 , -13.411926 ]],\n", - "\n", - " [[ 6.5706787 , 6.261444 , 5.76593 , ..., -10.512054 ,\n", - " -10.819702 , -11.010132 ],\n", - " [ 6.6134033 , 6.326996 , 5.8705444 , ..., -10.703583 ,\n", - " -10.812073 , -11.228058 ],\n", - " [ 6.63974 , 6.3713074 , 5.99823 , ..., -10.811035 ,\n", - " -11.002502 , -11.324585 ],\n", + " [10043.194 , 10672.565 , 11238.228 , ..., 7507.282 ,\n", + " 7414.5137 , 7444.7827 ],\n", + " [13826.692 , 14108.494 , 14563.739 , ..., 7519.0063 ,\n", + " 7214.311 , 7085.4043 ],\n", + " [14878.069 , 14873.55 , 14858.358 , ..., 7069.7866 ,\n", + " 7237.475 , 7302.5605 ]]],\n", + "\n", + "\n", + " [[[ 868.89185, 899.5873 , 930.3964 , ..., 936.80554,\n", + " 839.1485 , 772.7412 ],\n", + " [ 858.319 , 891.28815, 934.2618 , ..., 772.7412 ,\n", + " 733.67554, 716.23883],\n", + " [ 873.09827, 925.96265, 940.6993 , ..., 729.8386 ,\n", "...\n", - " [ 4.5212708 , 4.881195 , 5.349457 , ..., -0.38009644,\n", - " -0.14691162, -0.15356445],\n", - " [ 3.3636475 , 3.8493652 , 4.25 , ..., 0.430542 ,\n", - " 0.32699585, 0.23458862],\n", - " [ 2.291626 , 2.7035522 , 3.230133 , ..., 0.6605835 ,\n", - " 0.5307617 , 0.37503052]],\n", - "\n", - " [[ 9.923706 , 10.213959 , 10.246948 , ..., -2.7084656 ,\n", - " -3.117157 , -3.460785 ],\n", - " [ 10.048218 , 10.249817 , 10.259094 , ..., -2.855835 ,\n", - " -3.1342468 , -3.5550232 ],\n", - " [ 9.972595 , 10.138733 , 10.1493225 , ..., -3.176941 ,\n", - " -3.5891113 , -3.8943787 ],\n", + " 6790.51 , 6748.986 ],\n", + " [ 1623.8978 , 1855.4778 , 2173.332 , ..., 5969.662 ,\n", + " 6049.47 , 6131.5386 ],\n", + " [ 2051.7156 , 2147.7383 , 1941.084 , ..., 5981.9263 ,\n", + " 5846.8096 , 5826.9146 ]]],\n", + "\n", + "\n", + " [[[ 1362.0598 , 1371.5527 , 1349.8243 , ..., 2480.518 ,\n", + " 2496.761 , 2526.6606 ],\n", + " [ 1398.2976 , 1351.3591 , 1315.5477 , ..., 2294.1143 ,\n", + " 2327.9363 , 2370.2136 ],\n", + " [ 1287.1544 , 1304.108 , 1268.4387 , ..., 2281.4312 ,\n", + " 2239.3528 , 2212.153 ],\n", " ...,\n", - " [ 2.2316895 , 2.5167236 , 2.8789368 , ..., -1.0142822 ,\n", - " -1.0447388 , -1.3068542 ],\n", - " [ 0.93966675, 1.3753662 , 1.7021484 , ..., -0.4562683 ,\n", - " -0.82281494, -1.1673279 ],\n", - " [ -0.25891113, 0.1281128 , 0.6257019 , ..., -0.47140503,\n", - " -0.86395264, -1.2958069 ]]]], dtype=float32)\n", + " [ 1829.8173 , 1892.3165 , 1997.889 , ..., 6986.432 ,\n", + " 6914.4253 , 6935.685 ],\n", + " [ 1650.3838 , 1853.0875 , 2025.124 , ..., 6504.3003 ,\n", + " 6350.027 , 6176.655 ],\n", + " [ 1954.8159 , 2019.9797 , 1916.7521 , ..., 5647.4995 ,\n", + " 5836.5464 , 6040.9556 ]]]], dtype=float32)\n", "Coordinates:\n", - " * lat (lat) float64 55.13 56.07 57.02 57.96 ... 68.32 69.27 70.21 71.15\n", - " * lon (lon) float64 0.0 1.25 2.5 3.75 5.0 ... 38.75 40.0 41.25 42.5\n", - " * time (time) object 1850-01-15 12:00:00 ... 2014-12-15 12:00:00\n", - " * member_id (member_id) object 'r1i1p1f1'\n", + " * lat (lat) float32 69.95 69.85 69.75 69.65 ... 54.85 54.75 54.65 54.55\n", + " * level (level) float32 0.0\n", + " * time (time) timedelta64[ns] 00:00:00 01:00:00 ... 4 days 00:00:00\n", + " * lon (lon) float32 -2.45 -2.35 -2.25 -2.15 ... 42.15 42.25 42.35 42.45\n", "Attributes:\n", - " cell_measures: area: areacella\n", - " cell_methods: area: time: mean\n", - " comment: near-surface (usually, 2 meter) air temperature\n", - " description: near-surface (usually, 2 meter) air temperature\n", - " frequency: mon\n", - " id: tas\n", - " long_name: Near-Surface Air Temperature\n", - " mipTable: Amon\n", - " out_name: tas\n", - " prov: Amon ((isd.003))\n", - " realm: atmos\n", - " standard_name: air_temperature\n", - " time: time\n", - " time_label: time-mean\n", - " time_title: Temporal mean\n", - " title: Near-Surface Air Temperature\n", - " type: real\n", - " units: K\n", - " variable_id: tas
array([69.95, 69.85, 69.75, 69.65, 69.55, 69.45, 69.35, 69.25, 69.15, 69.05,\n", + " 68.95, 68.85, 68.75, 68.65, 68.55, 68.45, 68.35, 68.25, 68.15, 68.05,\n", + " 67.95, 67.85, 67.75, 67.65, 67.55, 67.45, 67.35, 67.25, 67.15, 67.05,\n", + " 66.95, 66.85, 66.75, 66.65, 66.55, 66.45, 66.35, 66.25, 66.15, 66.05,\n", + " 65.95, 65.85, 65.75, 65.65, 65.55, 65.45, 65.35, 65.25, 65.15, 65.05,\n", + " 64.95, 64.85, 64.75, 64.65, 64.55, 64.45, 64.35, 64.25, 64.15, 64.05,\n", + " 63.95, 63.85, 63.75, 63.65, 63.55, 63.45, 63.35, 63.25, 63.15, 63.05,\n", + " 62.95, 62.85, 62.75, 62.65, 62.55, 62.45, 62.35, 62.25, 62.15, 62.05,\n", + " 61.95, 61.85, 61.75, 61.65, 61.55, 61.45, 61.35, 61.25, 61.15, 61.05,\n", + " 60.95, 60.85, 60.75, 60.65, 60.55, 60.45, 60.35, 60.25, 60.15, 60.05,\n", + " 59.95, 59.85, 59.75, 59.65, 59.55, 59.45, 59.35, 59.25, 59.15, 59.05,\n", + " 58.95, 58.85, 58.75, 58.65, 58.55, 58.45, 58.35, 58.25, 58.15, 58.05,\n", + " 57.95, 57.85, 57.75, 57.65, 57.55, 57.45, 57.35, 57.25, 57.15, 57.05,\n", + " 56.95, 56.85, 56.75, 56.65, 56.55, 56.45, 56.35, 56.25, 56.15, 56.05,\n", + " 55.95, 55.85, 55.75, 55.65, 55.55, 55.45, 55.35, 55.25, 55.15, 55.05,\n", + " 54.95, 54.85, 54.75, 54.65, 54.55], dtype=float32)
array([0.], dtype=float32)
array([ 0, 3600000000000, 7200000000000, 10800000000000,\n", + " 14400000000000, 18000000000000, 21600000000000, 25200000000000,\n", + " 28800000000000, 32400000000000, 36000000000000, 39600000000000,\n", + " 43200000000000, 46800000000000, 50400000000000, 54000000000000,\n", + " 57600000000000, 61200000000000, 64800000000000, 68400000000000,\n", + " 72000000000000, 75600000000000, 79200000000000, 82800000000000,\n", + " 86400000000000, 90000000000000, 93600000000000, 97200000000000,\n", + " 100800000000000, 104400000000000, 108000000000000, 111600000000000,\n", + " 115200000000000, 118800000000000, 122400000000000, 126000000000000,\n", + " 129600000000000, 133200000000000, 136800000000000, 140400000000000,\n", + " 144000000000000, 147600000000000, 151200000000000, 154800000000000,\n", + " 158400000000000, 162000000000000, 165600000000000, 169200000000000,\n", + " 172800000000000, 176400000000000, 180000000000000, 183600000000000,\n", + " 187200000000000, 190800000000000, 194400000000000, 198000000000000,\n", + " 201600000000000, 205200000000000, 208800000000000, 212400000000000,\n", + " 216000000000000, 219600000000000, 223200000000000, 226800000000000,\n", + " 230400000000000, 234000000000000, 237600000000000, 241200000000000,\n", + " 244800000000000, 248400000000000, 252000000000000, 255600000000000,\n", + " 259200000000000, 262800000000000, 266400000000000, 270000000000000,\n", + " 273600000000000, 277200000000000, 280800000000000, 284400000000000,\n", + " 288000000000000, 291600000000000, 295200000000000, 298800000000000,\n", + " 302400000000000, 306000000000000, 309600000000000, 313200000000000,\n", + " 316800000000000, 320400000000000, 324000000000000, 327600000000000,\n", + " 331200000000000, 334800000000000, 338400000000000, 342000000000000,\n", + " 345600000000000], dtype='timedelta64[ns]')
array([-2.450012, -2.349976, -2.25 , ..., 42.25 , 42.350006, 42.449997],\n", + " dtype=float32)
PandasIndex(Index([ 69.94999694824219, 69.8499984741211, 69.75,\n", + " 69.6500015258789, 69.55000305175781, 69.44999694824219,\n", + " 69.3499984741211, 69.25, 69.1500015258789,\n", + " 69.05000305175781,\n", + " ...\n", + " 55.45000076293945, 55.349998474121094, 55.25,\n", + " 55.150001525878906, 55.04999923706055, 54.95000076293945,\n", + " 54.849998474121094, 54.75, 54.650001525878906,\n", + " 54.54999923706055],\n", + " dtype='float32', name='lat', length=155))
PandasIndex(Index([0.0], dtype='float32', name='level'))
PandasIndex(TimedeltaIndex(['0 days 00:00:00', '0 days 01:00:00', '0 days 02:00:00',\n", + " '0 days 03:00:00', '0 days 04:00:00', '0 days 05:00:00',\n", + " '0 days 06:00:00', '0 days 07:00:00', '0 days 08:00:00',\n", + " '0 days 09:00:00', '0 days 10:00:00', '0 days 11:00:00',\n", + " '0 days 12:00:00', '0 days 13:00:00', '0 days 14:00:00',\n", + " '0 days 15:00:00', '0 days 16:00:00', '0 days 17:00:00',\n", + " '0 days 18:00:00', '0 days 19:00:00', '0 days 20:00:00',\n", + " '0 days 21:00:00', '0 days 22:00:00', '0 days 23:00:00',\n", + " '1 days 00:00:00', '1 days 01:00:00', '1 days 02:00:00',\n", + " '1 days 03:00:00', '1 days 04:00:00', '1 days 05:00:00',\n", + " '1 days 06:00:00', '1 days 07:00:00', '1 days 08:00:00',\n", + " '1 days 09:00:00', '1 days 10:00:00', '1 days 11:00:00',\n", + " '1 days 12:00:00', '1 days 13:00:00', '1 days 14:00:00',\n", + " '1 days 15:00:00', '1 days 16:00:00', '1 days 17:00:00',\n", + " '1 days 18:00:00', '1 days 19:00:00', '1 days 20:00:00',\n", + " '1 days 21:00:00', '1 days 22:00:00', '1 days 23:00:00',\n", + " '2 days 00:00:00', '2 days 01:00:00', '2 days 02:00:00',\n", + " '2 days 03:00:00', '2 days 04:00:00', '2 days 05:00:00',\n", + " '2 days 06:00:00', '2 days 07:00:00', '2 days 08:00:00',\n", + " '2 days 09:00:00', '2 days 10:00:00', '2 days 11:00:00',\n", + " '2 days 12:00:00', '2 days 13:00:00', '2 days 14:00:00',\n", + " '2 days 15:00:00', '2 days 16:00:00', '2 days 17:00:00',\n", + " '2 days 18:00:00', '2 days 19:00:00', '2 days 20:00:00',\n", + " '2 days 21:00:00', '2 days 22:00:00', '2 days 23:00:00',\n", + " '3 days 00:00:00', '3 days 01:00:00', '3 days 02:00:00',\n", + " '3 days 03:00:00', '3 days 04:00:00', '3 days 05:00:00',\n", + " '3 days 06:00:00', '3 days 07:00:00', '3 days 08:00:00',\n", + " '3 days 09:00:00', '3 days 10:00:00', '3 days 11:00:00',\n", + " '3 days 12:00:00', '3 days 13:00:00', '3 days 14:00:00',\n", + " '3 days 15:00:00', '3 days 16:00:00', '3 days 17:00:00',\n", + " '3 days 18:00:00', '3 days 19:00:00', '3 days 20:00:00',\n", + " '3 days 21:00:00', '3 days 22:00:00', '3 days 23:00:00',\n", + " '4 days 00:00:00'],\n", + " dtype='timedelta64[ns]', name='time', freq=None))
PandasIndex(Index([ -2.45001220703125, -2.3499755859375, -2.25,\n", + " -2.1500244140625, -2.04998779296875, -1.95001220703125,\n", + " -1.8499755859375, -1.75, -1.6500244140625,\n", + " -1.54998779296875,\n", + " ...\n", + " 41.55000305175781, 41.649993896484375, 41.75,\n", + " 41.850006103515625, 41.94999694824219, 42.05000305175781,\n", + " 42.149993896484375, 42.25, 42.350006103515625,\n", + " 42.44999694824219],\n", + " dtype='float32', name='lon', length=450))
<xarray.DataArray 'tas' (member_id: 1, time: 1980, lat: 18, lon: 35)>\n", - "array([[[[ 3.9915466 , 3.9372864 , 3.63208 , ..., -17.947784 ,\n", - " -18.09903 , -17.901337 ],\n", - " [ 3.9906616 , 3.8585205 , 3.4658813 , ..., -19.287094 ,\n", - " -19.638443 , -19.815659 ],\n", - " [ 3.8772888 , 3.693756 , 3.3085938 , ..., -20.736603 ,\n", - " -21.195053 , -21.357864 ],\n", - " ...,\n", - " [ -0.11868286, 0.1071167 , 0.3340149 , ..., -11.844849 ,\n", - " -12.702545 , -14.013733 ],\n", - " [ -0.9946594 , -0.6482544 , -0.4326172 , ..., -10.813141 ,\n", - " -12.086548 , -13.261871 ],\n", - " [ -1.9762573 , -1.6532898 , -1.2767944 , ..., -10.560425 ,\n", - " -11.828857 , -13.411926 ]],\n", - "\n", - " [[ 6.5706787 , 6.261444 , 5.76593 , ..., -10.512054 ,\n", - " -10.819702 , -11.010132 ],\n", - " [ 6.6134033 , 6.326996 , 5.8705444 , ..., -10.703583 ,\n", - " -10.812073 , -11.228058 ],\n", - " [ 6.63974 , 6.3713074 , 5.99823 , ..., -10.811035 ,\n", - " -11.002502 , -11.324585 ],\n", - "...\n", - " [ 4.5212708 , 4.881195 , 5.349457 , ..., -0.38009644,\n", - " -0.14691162, -0.15356445],\n", - " [ 3.3636475 , 3.8493652 , 4.25 , ..., 0.430542 ,\n", - " 0.32699585, 0.23458862],\n", - " [ 2.291626 , 2.7035522 , 3.230133 , ..., 0.6605835 ,\n", - " 0.5307617 , 0.37503052]],\n", - "\n", - " [[ 9.923706 , 10.213959 , 10.246948 , ..., -2.7084656 ,\n", - " -3.117157 , -3.460785 ],\n", - " [ 10.048218 , 10.249817 , 10.259094 , ..., -2.855835 ,\n", - " -3.1342468 , -3.5550232 ],\n", - " [ 9.972595 , 10.138733 , 10.1493225 , ..., -3.176941 ,\n", - " -3.5891113 , -3.8943787 ],\n", - " ...,\n", - " [ 2.2316895 , 2.5167236 , 2.8789368 , ..., -1.0142822 ,\n", - " -1.0447388 , -1.3068542 ],\n", - " [ 0.93966675, 1.3753662 , 1.7021484 , ..., -0.4562683 ,\n", - " -0.82281494, -1.1673279 ],\n", - " [ -0.25891113, 0.1281128 , 0.6257019 , ..., -0.47140503,\n", - " -0.86395264, -1.2958069 ]]]], dtype=float32)\n", - "Coordinates:\n", - " * lat (lat) float64 55.13 56.07 57.02 57.96 ... 68.32 69.27 70.21 71.15\n", - " * lon (lon) float64 0.0 1.25 2.5 3.75 5.0 ... 38.75 40.0 41.25 42.5\n", - " * time (time) object 1850-01-15 12:00:00 ... 2014-12-15 12:00:00\n", - " * member_id (member_id) object 'r1i1p1f1'" + "
<xarray.DataArray 'pm2p5_conc' ()>\n", + "array(375.96038818)" ], "text/plain": [ - "
<xarray.DataArray 'tas' ()>\n", - "array(-26.84191895)" + "
<xarray.DataArray 'pm2p5_conc' (level: 1, lat: 155, lon: 450)>\n", + "array([[[ 2121.0813, 2219.4773, 2217.5305, ..., 3652.913 ,\n", + " 3473.7568, 3577.1978],\n", + " [ 2229.2544, 2233.2195, 2240.6362, ..., 4046.014 ,\n", + " 3602.7063, 3208.9517],\n", + " [ 2284.292 , 2295.2058, 2309.7578, ..., 4741.6777,\n", + " 4632.1265, 4161.8325],\n", + " ...,\n", + " [13408.843 , 13095.977 , 14038.214 , ..., 18508.65 ,\n", + " 18356.64 , 17824.658 ],\n", + " [17083.7 , 16497.771 , 16912.217 , ..., 19211.908 ,\n", + " 19491.762 , 18477.52 ],\n", + " [19680.336 , 19178.678 , 18560.22 , ..., 18558.336 ,\n", + " 18246.791 , 18232.215 ]]], dtype=float32)\n", + "Coordinates:\n", + " * lat (lat) float32 69.95 69.85 69.75 69.65 ... 54.85 54.75 54.65 54.55\n", + " * level (level) float32 0.0\n", + " * lon (lon) float32 -2.45 -2.35 -2.25 -2.15 ... 42.15 42.25 42.35 42.45" ], "text/plain": [ - "
<xarray.DataArray 'tas' (member_id: 1, lat: 18, lon: 35)>\n", - "array([[[17.913177 , 18.16388 , 18.275024 , 18.313934 , 18.328522 ,\n", - " 18.410065 , 18.476685 , 19.338104 , 20.345673 , 19.73639 ,\n", - " 19.681091 , 20.134735 , 19.119598 , 19.005066 , 19.041931 ,\n", - " 19.106995 , 19.019653 , 22.693756 , 24.226562 , 24.650146 ,\n", - " 24.417786 , 23.973938 , 23.803497 , 22.828125 , 22.704803 ,\n", - " 22.352753 , 22.175903 , 22.134857 , 23.450043 , 24.378357 ,\n", - " 25.166748 , 25.610962 , 24.972687 , 25.00238 , 24.528198 ],\n", - " [17.241577 , 17.338257 , 17.492157 , 17.585663 , 17.65509 ,\n", - " 17.713348 , 17.98352 , 20.150513 , 21.748383 , 21.886383 ,\n", - " 21.556824 , 20.893677 , 20.363037 , 18.943329 , 18.613983 ,\n", - " 18.611572 , 18.620941 , 19.760834 , 23.221924 , 23.647827 ,\n", - " 23.514404 , 23.093231 , 23.336456 , 22.503937 , 22.339935 ,\n", - " 22.338959 , 22.166626 , 22.34134 , 23.142212 , 23.4458 ,\n", - " 24.20691 , 24.140076 , 24.62323 , 24.352112 , 24.227814 ],\n", - " [16.807587 , 16.770233 , 16.859833 , 16.93045 , 16.981567 ,\n", - " 17.039246 , 17.455383 , 18.455078 , 19.253387 , 20.740967 ,\n", - " 21.366882 , 21.06958 , 20.757294 , 20.318634 , 18.286224 ,\n", - " 18.358154 , 18.349121 , 19.311615 , 21.430908 , 22.219482 ,\n", - " 22.33014 , 22.559174 , 22.694366 , 22.51947 , 22.080627 ,\n", - " 22.51355 , 22.614471 , 22.480713 , 22.66214 , 23.090118 ,\n", + "<xarray.DataArray 'pm2p5_conc' (time: 97, level: 1, lat: 155, lon: 450)>\n", + "array([[[[ 855.6641 , 852.48096, 833.32465, ..., 879.8937 ,\n", + " 902.90106, 865.896 ],\n", + " [ 885.7343 , 891.02075, 878.97 , ..., 971.65314,\n", + " 1061.4657 , 1015.5647 ],\n", + " [ 893.89136, 828.7204 , 800.25604, ..., 1032.3335 ,\n", + " 990.852 , 982.55286],\n", + " ...,\n", + " [10043.194 , 10672.565 , 11238.228 , ..., 7507.282 ,\n", + " 7414.5137 , 7444.7827 ],\n", + " [13826.692 , 14108.494 , 14563.739 , ..., 7519.0063 ,\n", + " 7214.311 , 7085.4043 ],\n", + " [14878.069 , 14873.55 , 14858.358 , ..., 7069.7866 ,\n", + " 7237.475 , 7302.5605 ]]],\n", + "\n", + "\n", + " [[[ 868.89185, 899.5873 , 930.3964 , ..., 936.80554,\n", + " 839.1485 , 772.7412 ],\n", + " [ 858.319 , 891.28815, 934.2618 , ..., 772.7412 ,\n", + " 733.67554, 716.23883],\n", + " [ 873.09827, 925.96265, 940.6993 , ..., 729.8386 ,\n", "...\n", - " 12.900909 , 13.222961 , 13.52597 , 14.403259 , 15.553009 ,\n", - " 16.8042 , 16.932892 , 16.954376 , 17.034729 , 16.96762 ,\n", - " 16.96637 , 17.100616 , 17.330872 , 17.345306 , 17.114044 ,\n", - " 16.488129 , 14.975281 , 12.001068 , 13.143677 , 10.707916 ,\n", - " 9.694275 , 9.177216 , 9.027313 , 8.803345 , 8.669434 ],\n", - " [ 9.139069 , 9.590576 , 9.904144 , 10.247589 , 10.48056 ,\n", - " 10.866425 , 11.181763 , 11.512421 , 11.786804 , 12.199493 ,\n", - " 12.608215 , 12.951721 , 13.20874 , 13.357086 , 13.322723 ,\n", - " 14.709045 , 13.948395 , 15.049042 , 16.412964 , 17.884033 ,\n", - " 17.783173 , 17.81784 , 17.365814 , 17.036255 , 16.102875 ,\n", - " 12.359375 , 9.921783 , 9.611908 , 9.409424 , 9.122467 ,\n", - " 8.811737 , 8.50293 , 8.18045 , 7.9012146, 7.5783386],\n", - " [ 8.040253 , 8.353821 , 8.763489 , 9.011475 , 9.314575 ,\n", - " 9.718872 , 10.178741 , 10.651154 , 11.184814 , 11.681763 ,\n", - " 12.085083 , 12.411102 , 12.637634 , 12.749329 , 12.723297 ,\n", - " 12.770172 , 12.895966 , 12.92099 , 12.873657 , 13.131958 ,\n", - " 14.395782 , 12.900848 , 12.762329 , 11.248535 , 10.335876 ,\n", - " 9.764557 , 9.28595 , 8.989227 , 8.640686 , 8.288513 ,\n", - " 7.905426 , 7.5679626, 7.2664795, 7.0660706, 6.9654236]]],\n", - " dtype=float32)\n", + " 6790.51 , 6748.986 ],\n", + " [ 1623.8978 , 1855.4778 , 2173.332 , ..., 5969.662 ,\n", + " 6049.47 , 6131.5386 ],\n", + " [ 2051.7156 , 2147.7383 , 1941.084 , ..., 5981.9263 ,\n", + " 5846.8096 , 5826.9146 ]]],\n", + "\n", + "\n", + " [[[ 1362.0598 , 1371.5527 , 1349.8243 , ..., 2480.518 ,\n", + " 2496.761 , 2526.6606 ],\n", + " [ 1398.2976 , 1351.3591 , 1315.5477 , ..., 2294.1143 ,\n", + " 2327.9363 , 2370.2136 ],\n", + " [ 1287.1544 , 1304.108 , 1268.4387 , ..., 2281.4312 ,\n", + " 2239.3528 , 2212.153 ],\n", + " ...,\n", + " [ 1829.8173 , 1892.3165 , 1997.889 , ..., 6986.432 ,\n", + " 6914.4253 , 6935.685 ],\n", + " [ 1650.3838 , 1853.0875 , 2025.124 , ..., 6504.3003 ,\n", + " 6350.027 , 6176.655 ],\n", + " [ 1954.8159 , 2019.9797 , 1916.7521 , ..., 5647.4995 ,\n", + " 5836.5464 , 6040.9556 ]]]], dtype=float32)\n", "Coordinates:\n", - " * lat (lat) float64 55.13 56.07 57.02 57.96 ... 68.32 69.27 70.21 71.15\n", - " * lon (lon) float64 0.0 1.25 2.5 3.75 5.0 ... 38.75 40.0 41.25 42.5\n", - " * member_id (member_id) object 'r1i1p1f1'
array([69.95, 69.85, 69.75, 69.65, 69.55, 69.45, 69.35, 69.25, 69.15, 69.05,\n", + " 68.95, 68.85, 68.75, 68.65, 68.55, 68.45, 68.35, 68.25, 68.15, 68.05,\n", + " 67.95, 67.85, 67.75, 67.65, 67.55, 67.45, 67.35, 67.25, 67.15, 67.05,\n", + " 66.95, 66.85, 66.75, 66.65, 66.55, 66.45, 66.35, 66.25, 66.15, 66.05,\n", + " 65.95, 65.85, 65.75, 65.65, 65.55, 65.45, 65.35, 65.25, 65.15, 65.05,\n", + " 64.95, 64.85, 64.75, 64.65, 64.55, 64.45, 64.35, 64.25, 64.15, 64.05,\n", + " 63.95, 63.85, 63.75, 63.65, 63.55, 63.45, 63.35, 63.25, 63.15, 63.05,\n", + " 62.95, 62.85, 62.75, 62.65, 62.55, 62.45, 62.35, 62.25, 62.15, 62.05,\n", + " 61.95, 61.85, 61.75, 61.65, 61.55, 61.45, 61.35, 61.25, 61.15, 61.05,\n", + " 60.95, 60.85, 60.75, 60.65, 60.55, 60.45, 60.35, 60.25, 60.15, 60.05,\n", + " 59.95, 59.85, 59.75, 59.65, 59.55, 59.45, 59.35, 59.25, 59.15, 59.05,\n", + " 58.95, 58.85, 58.75, 58.65, 58.55, 58.45, 58.35, 58.25, 58.15, 58.05,\n", + " 57.95, 57.85, 57.75, 57.65, 57.55, 57.45, 57.35, 57.25, 57.15, 57.05,\n", + " 56.95, 56.85, 56.75, 56.65, 56.55, 56.45, 56.35, 56.25, 56.15, 56.05,\n", + " 55.95, 55.85, 55.75, 55.65, 55.55, 55.45, 55.35, 55.25, 55.15, 55.05,\n", + " 54.95, 54.85, 54.75, 54.65, 54.55], dtype=float32)
array([0.], dtype=float32)
array([-2.450012, -2.349976, -2.25 , ..., 42.25 , 42.350006, 42.449997],\n", + " dtype=float32)
array(['2021-12-22T00:00:00.000000000', '2021-12-22T01:00:00.000000000',\n", + " '2021-12-22T02:00:00.000000000', '2021-12-22T03:00:00.000000000',\n", + " '2021-12-22T04:00:00.000000000', '2021-12-22T05:00:00.000000000',\n", + " '2021-12-22T06:00:00.000000000', '2021-12-22T07:00:00.000000000',\n", + " '2021-12-22T08:00:00.000000000', '2021-12-22T09:00:00.000000000',\n", + " '2021-12-22T10:00:00.000000000', '2021-12-22T11:00:00.000000000',\n", + " '2021-12-22T12:00:00.000000000', '2021-12-22T13:00:00.000000000',\n", + " '2021-12-22T14:00:00.000000000', '2021-12-22T15:00:00.000000000',\n", + " '2021-12-22T16:00:00.000000000', '2021-12-22T17:00:00.000000000',\n", + " '2021-12-22T18:00:00.000000000', '2021-12-22T19:00:00.000000000',\n", + " '2021-12-22T20:00:00.000000000', '2021-12-22T21:00:00.000000000',\n", + " '2021-12-22T22:00:00.000000000', '2021-12-22T23:00:00.000000000',\n", + " '2021-12-23T00:00:00.000000000', '2021-12-23T01:00:00.000000000',\n", + " '2021-12-23T02:00:00.000000000', '2021-12-23T03:00:00.000000000',\n", + " '2021-12-23T04:00:00.000000000', '2021-12-23T05:00:00.000000000',\n", + " '2021-12-23T06:00:00.000000000', '2021-12-23T07:00:00.000000000',\n", + " '2021-12-23T08:00:00.000000000', '2021-12-23T09:00:00.000000000',\n", + " '2021-12-23T10:00:00.000000000', '2021-12-23T11:00:00.000000000',\n", + " '2021-12-23T12:00:00.000000000', '2021-12-23T13:00:00.000000000',\n", + " '2021-12-23T14:00:00.000000000', '2021-12-23T15:00:00.000000000',\n", + " '2021-12-23T16:00:00.000000000', '2021-12-23T17:00:00.000000000',\n", + " '2021-12-23T18:00:00.000000000', '2021-12-23T19:00:00.000000000',\n", + " '2021-12-23T20:00:00.000000000', '2021-12-23T21:00:00.000000000',\n", + " '2021-12-23T22:00:00.000000000', '2021-12-23T23:00:00.000000000',\n", + " '2021-12-24T00:00:00.000000000', '2021-12-24T01:00:00.000000000',\n", + " '2021-12-24T02:00:00.000000000', '2021-12-24T03:00:00.000000000',\n", + " '2021-12-24T04:00:00.000000000', '2021-12-24T05:00:00.000000000',\n", + " '2021-12-24T06:00:00.000000000', '2021-12-24T07:00:00.000000000',\n", + " '2021-12-24T08:00:00.000000000', '2021-12-24T09:00:00.000000000',\n", + " '2021-12-24T10:00:00.000000000', '2021-12-24T11:00:00.000000000',\n", + " '2021-12-24T12:00:00.000000000', '2021-12-24T13:00:00.000000000',\n", + " '2021-12-24T14:00:00.000000000', '2021-12-24T15:00:00.000000000',\n", + " '2021-12-24T16:00:00.000000000', '2021-12-24T17:00:00.000000000',\n", + " '2021-12-24T18:00:00.000000000', '2021-12-24T19:00:00.000000000',\n", + " '2021-12-24T20:00:00.000000000', '2021-12-24T21:00:00.000000000',\n", + " '2021-12-24T22:00:00.000000000', '2021-12-24T23:00:00.000000000',\n", + " '2021-12-25T00:00:00.000000000', '2021-12-25T01:00:00.000000000',\n", + " '2021-12-25T02:00:00.000000000', '2021-12-25T03:00:00.000000000',\n", + " '2021-12-25T04:00:00.000000000', '2021-12-25T05:00:00.000000000',\n", + " '2021-12-25T06:00:00.000000000', '2021-12-25T07:00:00.000000000',\n", + " '2021-12-25T08:00:00.000000000', '2021-12-25T09:00:00.000000000',\n", + " '2021-12-25T10:00:00.000000000', '2021-12-25T11:00:00.000000000',\n", + " '2021-12-25T12:00:00.000000000', '2021-12-25T13:00:00.000000000',\n", + " '2021-12-25T14:00:00.000000000', '2021-12-25T15:00:00.000000000',\n", + " '2021-12-25T16:00:00.000000000', '2021-12-25T17:00:00.000000000',\n", + " '2021-12-25T18:00:00.000000000', '2021-12-25T19:00:00.000000000',\n", + " '2021-12-25T20:00:00.000000000', '2021-12-25T21:00:00.000000000',\n", + " '2021-12-25T22:00:00.000000000', '2021-12-25T23:00:00.000000000',\n", + " '2021-12-26T00:00:00.000000000'], dtype='datetime64[ns]')
PandasIndex(Index([ 69.94999694824219, 69.8499984741211, 69.75,\n", + " 69.6500015258789, 69.55000305175781, 69.44999694824219,\n", + " 69.3499984741211, 69.25, 69.1500015258789,\n", + " 69.05000305175781,\n", + " ...\n", + " 55.45000076293945, 55.349998474121094, 55.25,\n", + " 55.150001525878906, 55.04999923706055, 54.95000076293945,\n", + " 54.849998474121094, 54.75, 54.650001525878906,\n", + " 54.54999923706055],\n", + " dtype='float32', name='lat', length=155))
PandasIndex(Index([0.0], dtype='float32', name='level'))
PandasIndex(Index([ -2.45001220703125, -2.3499755859375, -2.25,\n", + " -2.1500244140625, -2.04998779296875, -1.95001220703125,\n", + " -1.8499755859375, -1.75, -1.6500244140625,\n", + " -1.54998779296875,\n", + " ...\n", + " 41.55000305175781, 41.649993896484375, 41.75,\n", + " 41.850006103515625, 41.94999694824219, 42.05000305175781,\n", + " 42.149993896484375, 42.25, 42.350006103515625,\n", + " 42.44999694824219],\n", + " dtype='float32', name='lon', length=450))
PandasIndex(DatetimeIndex(['2021-12-22 00:00:00', '2021-12-22 01:00:00',\n", + " '2021-12-22 02:00:00', '2021-12-22 03:00:00',\n", + " '2021-12-22 04:00:00', '2021-12-22 05:00:00',\n", + " '2021-12-22 06:00:00', '2021-12-22 07:00:00',\n", + " '2021-12-22 08:00:00', '2021-12-22 09:00:00',\n", + " '2021-12-22 10:00:00', '2021-12-22 11:00:00',\n", + " '2021-12-22 12:00:00', '2021-12-22 13:00:00',\n", + " '2021-12-22 14:00:00', '2021-12-22 15:00:00',\n", + " '2021-12-22 16:00:00', '2021-12-22 17:00:00',\n", + " '2021-12-22 18:00:00', '2021-12-22 19:00:00',\n", + " '2021-12-22 20:00:00', '2021-12-22 21:00:00',\n", + " '2021-12-22 22:00:00', '2021-12-22 23:00:00',\n", + " '2021-12-23 00:00:00', '2021-12-23 01:00:00',\n", + " '2021-12-23 02:00:00', '2021-12-23 03:00:00',\n", + " '2021-12-23 04:00:00', '2021-12-23 05:00:00',\n", + " '2021-12-23 06:00:00', '2021-12-23 07:00:00',\n", + " '2021-12-23 08:00:00', '2021-12-23 09:00:00',\n", + " '2021-12-23 10:00:00', '2021-12-23 11:00:00',\n", + " '2021-12-23 12:00:00', '2021-12-23 13:00:00',\n", + " '2021-12-23 14:00:00', '2021-12-23 15:00:00',\n", + " '2021-12-23 16:00:00', '2021-12-23 17:00:00',\n", + " '2021-12-23 18:00:00', '2021-12-23 19:00:00',\n", + " '2021-12-23 20:00:00', '2021-12-23 21:00:00',\n", + " '2021-12-23 22:00:00', '2021-12-23 23:00:00',\n", + " '2021-12-24 00:00:00', '2021-12-24 01:00:00',\n", + " '2021-12-24 02:00:00', '2021-12-24 03:00:00',\n", + " '2021-12-24 04:00:00', '2021-12-24 05:00:00',\n", + " '2021-12-24 06:00:00', '2021-12-24 07:00:00',\n", + " '2021-12-24 08:00:00', '2021-12-24 09:00:00',\n", + " '2021-12-24 10:00:00', '2021-12-24 11:00:00',\n", + " '2021-12-24 12:00:00', '2021-12-24 13:00:00',\n", + " '2021-12-24 14:00:00', '2021-12-24 15:00:00',\n", + " '2021-12-24 16:00:00', '2021-12-24 17:00:00',\n", + " '2021-12-24 18:00:00', '2021-12-24 19:00:00',\n", + " '2021-12-24 20:00:00', '2021-12-24 21:00:00',\n", + " '2021-12-24 22:00:00', '2021-12-24 23:00:00',\n", + " '2021-12-25 00:00:00', '2021-12-25 01:00:00',\n", + " '2021-12-25 02:00:00', '2021-12-25 03:00:00',\n", + " '2021-12-25 04:00:00', '2021-12-25 05:00:00',\n", + " '2021-12-25 06:00:00', '2021-12-25 07:00:00',\n", + " '2021-12-25 08:00:00', '2021-12-25 09:00:00',\n", + " '2021-12-25 10:00:00', '2021-12-25 11:00:00',\n", + " '2021-12-25 12:00:00', '2021-12-25 13:00:00',\n", + " '2021-12-25 14:00:00', '2021-12-25 15:00:00',\n", + " '2021-12-25 16:00:00', '2021-12-25 17:00:00',\n", + " '2021-12-25 18:00:00', '2021-12-25 19:00:00',\n", + " '2021-12-25 20:00:00', '2021-12-25 21:00:00',\n", + " '2021-12-25 22:00:00', '2021-12-25 23:00:00',\n", + " '2021-12-26 00:00:00'],\n", + " dtype='datetime64[ns]', name='time', freq=None))
<xarray.DataArray 'tas' (member_id: 1, year: 165, lat: 18, lon: 35)>\n", - "array([[[[10.190738 , 10.236267 , 10.129463 , ..., 4.487551 ,\n", - " 4.577349 , 4.4684777 ],\n", - " [10.055203 , 10.016124 , 9.869583 , ..., 3.5996563 ,\n", - " 3.602525 , 3.5518277 ],\n", - " [ 9.951228 , 9.84331 , 9.695267 , ..., 2.8245113 ,\n", - " 2.760255 , 2.7445972 ],\n", + "<xarray.DataArray 'pm2p5_conc' (day: 5, level: 1, lat: 155, lon: 450)>\n", + "array([[[[ 1169.239 , 1166.6467 , 1160.973 , ..., 1441.5961 ,\n", + " 1446.9415 , 1449.0151 ],\n", + " [ 1173.6224 , 1171.9554 , 1171.7555 , ..., 1428.3876 ,\n", + " 1443.7915 , 1446.5171 ],\n", + " [ 1181.1288 , 1189.5563 , 1205.5222 , ..., 1435.0128 ,\n", + " 1436.0138 , 1443.3612 ],\n", " ...,\n", - " [ 2.131218 , 2.480194 , 2.8605297 , ..., -4.5990014 ,\n", - " -4.9325995 , -5.297928 ],\n", - " [ 0.8806508 , 1.3956375 , 1.7613932 , ..., -4.671717 ,\n", - " -5.130834 , -5.525271 ],\n", - " [-0.27815247, 0.1952184 , 0.7522685 , ..., -4.968567 ,\n", - " -5.4446716 , -5.893354 ]],\n", - "\n", - " [[10.707288 , 10.862836 , 10.859578 , ..., 7.673986 ,\n", - " 7.586838 , 7.378306 ],\n", - " [10.575684 , 10.646342 , 10.60096 , ..., 7.274389 ,\n", - " 7.101011 , 6.786873 ],\n", - " [10.451709 , 10.4361925 , 10.362701 , ..., 6.5675125 ,\n", - " 6.483808 , 6.3353195 ],\n", + " [ 9959.0625 , 9706.599 , 10212.829 , ..., 8834.37 ,\n", + " 8843.07 , 8738.081 ],\n", + " [11996.876 , 11674.386 , 11926.778 , ..., 8548.44 ,\n", + " 8362.801 , 8002.9683 ],\n", + " [13141.722 , 12680.837 , 12554.7 , ..., 7965.3184 ,\n", + " 7982.112 , 7990.6177 ]]],\n", + "\n", + "\n", + " [[[ 1686.1406 , 1708.1188 , 1722.3102 , ..., 1365.6635 ,\n", + " 1362.4203 , 1382.3403 ],\n", + " [ 1709.118 , 1726.6967 , 1740.0448 , ..., 1310.9231 ,\n", + " 1319.0963 , 1326.693 ],\n", + " [ 1732.6241 , 1738.5929 , 1748.2808 , ..., 1290.9888 ,\n", "...\n", - " [ 3.899109 , 4.278895 , 4.683767 , ..., 0.22347514,\n", - " 0.07540639, -0.13505554],\n", - " [ 2.6409404 , 3.1836903 , 3.5581539 , ..., 0.25100455,\n", - " -0.05962372, -0.3530426 ],\n", - " [ 1.38296 , 1.857727 , 2.4163895 , ..., 0.07100169,\n", - " -0.25650534, -0.61615247]],\n", - "\n", - " [[11.859332 , 12.051389 , 12.075607 , ..., 8.140752 ,\n", - " 8.121837 , 7.935308 ],\n", - " [11.677655 , 11.791171 , 11.787402 , ..., 7.7124076 ,\n", - " 7.5722404 , 7.4204307 ],\n", - " [11.52195 , 11.550827 , 11.534116 , ..., 7.284935 ,\n", - " 7.149447 , 7.030052 ],\n", + " 4905.723 , 4991.6763 ],\n", + " [ 1478.167 , 1535.4425 , 1637.3916 , ..., 4755.4307 ,\n", + " 4620.1465 , 4505.0903 ],\n", + " [ 1717.3419 , 1729.2365 , 1662.1465 , ..., 4393.5034 ,\n", + " 4420.8423 , 4419.639 ]]],\n", + "\n", + "\n", + " [[[ 1362.0598 , 1371.5527 , 1349.8243 , ..., 2480.518 ,\n", + " 2496.761 , 2526.6606 ],\n", + " [ 1398.2976 , 1351.3591 , 1315.5477 , ..., 2294.1143 ,\n", + " 2327.9363 , 2370.2136 ],\n", + " [ 1287.1544 , 1304.108 , 1268.4387 , ..., 2281.4312 ,\n", + " 2239.3528 , 2212.153 ],\n", " ...,\n", - " [ 4.1399918 , 4.518181 , 4.9494348 , ..., 1.4126383 ,\n", - " 1.2889252 , 1.1015829 ],\n", - " [ 2.8445637 , 3.3857346 , 3.787501 , ..., 1.4600092 ,\n", - " 1.1696981 , 0.8996175 ],\n", - " [ 1.5644354 , 2.0351918 , 2.6172917 , ..., 1.269165 ,\n", - " 0.9806595 , 0.6546936 ]]]], dtype=float32)\n", + " [ 1829.8173 , 1892.3165 , 1997.889 , ..., 6986.432 ,\n", + " 6914.4253 , 6935.685 ],\n", + " [ 1650.3838 , 1853.0875 , 2025.124 , ..., 6504.3003 ,\n", + " 6350.027 , 6176.655 ],\n", + " [ 1954.8159 , 2019.9797 , 1916.7521 , ..., 5647.4995 ,\n", + " 5836.5464 , 6040.9556 ]]]], dtype=float32)\n", "Coordinates:\n", - " * lat (lat) float64 55.13 56.07 57.02 57.96 ... 68.32 69.27 70.21 71.15\n", - " * lon (lon) float64 0.0 1.25 2.5 3.75 5.0 ... 38.75 40.0 41.25 42.5\n", - " * member_id (member_id) object 'r1i1p1f1'\n", - " * year (year) int64 1850 1851 1852 1853 1854 ... 2011 2012 2013 2014
array([69.95, 69.85, 69.75, 69.65, 69.55, 69.45, 69.35, 69.25, 69.15, 69.05,\n", + " 68.95, 68.85, 68.75, 68.65, 68.55, 68.45, 68.35, 68.25, 68.15, 68.05,\n", + " 67.95, 67.85, 67.75, 67.65, 67.55, 67.45, 67.35, 67.25, 67.15, 67.05,\n", + " 66.95, 66.85, 66.75, 66.65, 66.55, 66.45, 66.35, 66.25, 66.15, 66.05,\n", + " 65.95, 65.85, 65.75, 65.65, 65.55, 65.45, 65.35, 65.25, 65.15, 65.05,\n", + " 64.95, 64.85, 64.75, 64.65, 64.55, 64.45, 64.35, 64.25, 64.15, 64.05,\n", + " 63.95, 63.85, 63.75, 63.65, 63.55, 63.45, 63.35, 63.25, 63.15, 63.05,\n", + " 62.95, 62.85, 62.75, 62.65, 62.55, 62.45, 62.35, 62.25, 62.15, 62.05,\n", + " 61.95, 61.85, 61.75, 61.65, 61.55, 61.45, 61.35, 61.25, 61.15, 61.05,\n", + " 60.95, 60.85, 60.75, 60.65, 60.55, 60.45, 60.35, 60.25, 60.15, 60.05,\n", + " 59.95, 59.85, 59.75, 59.65, 59.55, 59.45, 59.35, 59.25, 59.15, 59.05,\n", + " 58.95, 58.85, 58.75, 58.65, 58.55, 58.45, 58.35, 58.25, 58.15, 58.05,\n", + " 57.95, 57.85, 57.75, 57.65, 57.55, 57.45, 57.35, 57.25, 57.15, 57.05,\n", + " 56.95, 56.85, 56.75, 56.65, 56.55, 56.45, 56.35, 56.25, 56.15, 56.05,\n", + " 55.95, 55.85, 55.75, 55.65, 55.55, 55.45, 55.35, 55.25, 55.15, 55.05,\n", + " 54.95, 54.85, 54.75, 54.65, 54.55], dtype=float32)
array([0.], dtype=float32)
array([-2.450012, -2.349976, -2.25 , ..., 42.25 , 42.350006, 42.449997],\n", + " dtype=float32)
array([22, 23, 24, 25, 26])
PandasIndex(Index([ 69.94999694824219, 69.8499984741211, 69.75,\n", + " 69.6500015258789, 69.55000305175781, 69.44999694824219,\n", + " 69.3499984741211, 69.25, 69.1500015258789,\n", + " 69.05000305175781,\n", + " ...\n", + " 55.45000076293945, 55.349998474121094, 55.25,\n", + " 55.150001525878906, 55.04999923706055, 54.95000076293945,\n", + " 54.849998474121094, 54.75, 54.650001525878906,\n", + " 54.54999923706055],\n", + " dtype='float32', name='lat', length=155))
PandasIndex(Index([0.0], dtype='float32', name='level'))
PandasIndex(Index([ -2.45001220703125, -2.3499755859375, -2.25,\n", + " -2.1500244140625, -2.04998779296875, -1.95001220703125,\n", + " -1.8499755859375, -1.75, -1.6500244140625,\n", + " -1.54998779296875,\n", + " ...\n", + " 41.55000305175781, 41.649993896484375, 41.75,\n", + " 41.850006103515625, 41.94999694824219, 42.05000305175781,\n", + " 42.149993896484375, 42.25, 42.350006103515625,\n", + " 42.44999694824219],\n", + " dtype='float32', name='lon', length=450))
PandasIndex(Index([22, 23, 24, 25, 26], dtype='int64', name='day'))
<xarray.DataArray 'year' (year: 165)>\n", - "array([1850, 1851, 1852, 1853, 1854, 1855, 1856, 1857, 1858, 1859, 1860, 1861,\n", - " 1862, 1863, 1864, 1865, 1866, 1867, 1868, 1869, 1870, 1871, 1872, 1873,\n", - " 1874, 1875, 1876, 1877, 1878, 1879, 1880, 1881, 1882, 1883, 1884, 1885,\n", - " 1886, 1887, 1888, 1889, 1890, 1891, 1892, 1893, 1894, 1895, 1896, 1897,\n", - " 1898, 1899, 1900, 1901, 1902, 1903, 1904, 1905, 1906, 1907, 1908, 1909,\n", - " 1910, 1911, 1912, 1913, 1914, 1915, 1916, 1917, 1918, 1919, 1920, 1921,\n", - " 1922, 1923, 1924, 1925, 1926, 1927, 1928, 1929, 1930, 1931, 1932, 1933,\n", - " 1934, 1935, 1936, 1937, 1938, 1939, 1940, 1941, 1942, 1943, 1944, 1945,\n", - " 1946, 1947, 1948, 1949, 1950, 1951, 1952, 1953, 1954, 1955, 1956, 1957,\n", - " 1958, 1959, 1960, 1961, 1962, 1963, 1964, 1965, 1966, 1967, 1968, 1969,\n", - " 1970, 1971, 1972, 1973, 1974, 1975, 1976, 1977, 1978, 1979, 1980, 1981,\n", - " 1982, 1983, 1984, 1985, 1986, 1987, 1988, 1989, 1990, 1991, 1992, 1993,\n", - " 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005,\n", - " 2006, 2007, 2008, 2009, 2010, 2011, 2012, 2013, 2014])\n", + "" ], "text/plain": [ - "<xarray.DataArray 'day' (day: 5)>\n", + "array([22, 23, 24, 25, 26])\n", "Coordinates:\n", - " * year (year) int64 1850 1851 1852 1853 1854 ... 2010 2011 2012 2013 2014" + " * day (day) int64 22 23 24 25 26\n", + "Attributes:\n", + " long_name: FORECAST time from 20211222