|
21 | 21 | "</ul>\n",
|
22 | 22 | "</div>\n",
|
23 | 23 | "\n",
|
24 |
| - "The datasets are stored on the [PAVICS THREDDS server](https://boreas.ouranos.ca/twitcher/ows/proxy/thredds/catalog/birdhouse/eccc/RDRS_v21/catalog.html). You'll find data at both hourly and daily frequencies. Hourly data are stored in individual files for each month and variable. Daily data are stored in files holding three years of data for each variable. \n", |
| 24 | + "The datasets are stored on the [PAVICS THREDDS server](https://pavics.ouranos.ca/twitcher/ows/proxy/thredds/catalog/birdhouse/eccc/RDRS_v21/catalog.html). You'll find data at both hourly and daily frequencies. Hourly data are stored in individual files for each month and variable. Daily data are stored in files holding three years of data for each variable. \n", |
25 | 25 | "\n",
|
26 |
| - "For convenience, we provide [*NcML aggregations*](#ncml), that is, single links for the [hourly](https://boreas.ouranos.ca/twitcher/ows/proxy/thredds/catalog/datasets/reanalyses/catalog.html?dataset=datasets/reanalyses/1hr_RDRSv2.1_NAM.ncml) and [daily](https://boreas.ouranos.ca/twitcher/ows/proxy/thredds/catalog/datasets/reanalyses/catalog.html?dataset=datasets/reanalyses/day_RDRSv2.1_NAM.ncml) datasets, combining multiple netCDF files (~500 files per variable for hourly data) into a single virtual dataset. The following demonstration uses those NcML aggregations. \n", |
| 26 | + "For convenience, we provide [*NcML aggregations*](#ncml), that is, single links for the [hourly](https://pavics.ouranos.ca/twitcher/ows/proxy/thredds/catalog/datasets/reanalyses/catalog.html?dataset=datasets/reanalyses/1hr_RDRSv2.1_NAM.ncml) and [daily](https://pavics.ouranos.ca/twitcher/ows/proxy/thredds/catalog/datasets/reanalyses/catalog.html?dataset=datasets/reanalyses/day_RDRSv2.1_NAM.ncml) datasets, combining multiple netCDF files (~500 files per variable for hourly data) into a single virtual dataset. The following demonstration uses those NcML aggregations. \n", |
27 | 27 | "\n",
|
28 | 28 | "## Data access\n",
|
29 | 29 | "\n",
|
|
2950 | 2950 | "source": [
|
2951 | 2951 | "import xarray as xr\n",
|
2952 | 2952 | "\n",
|
2953 |
| - "url = \"https://boreas.ouranos.ca/twitcher/ows/proxy/thredds/dodsC/datasets/reanalyses/day_RDRSv2.1_NAM.ncml\"\n", |
| 2953 | + "url = \"https://pavics.ouranos.ca/twitcher/ows/proxy/thredds/dodsC/datasets/reanalyses/day_RDRSv2.1_NAM.ncml\"\n", |
2954 | 2954 | "\n",
|
2955 | 2955 | "# Open dataset. For hourly data, we suggest using the following chunking pattern: dict(time=720, rlon=50, rlat=50)\n",
|
2956 | 2956 | "ds = xr.open_dataset(url, chunks=dict(time=1461, rlon=50, rlat=50))\n",
|
|
2966 | 2966 | "\n",
|
2967 | 2967 | "We then subset the data for a zone and period of interest, namely eastern Canada for the period 1981-2010, using [clisops subset_bbox](https://clisops.readthedocs.io/en/latest/api.html#clisops.core.subset.subset_bbox) function. We also demonstrate how to select a single grid cell according to its longitude and latitude using [clisops subset_gridpoint](https://clisops.readthedocs.io/en/latest/api.html#clisops.core.subset.subset_gridpoint). Note that we can't simply use xarray's `sel` method here because dimensions are in rotated pole coordinates. \n",
|
2968 | 2968 | "\n",
|
2969 |
| - "Other, and more advanced subsetting capabilities are also available. For more information and examples, consult those [subsetting examples](https://boreas.ouranos.ca/climate_analysis.html#b). " |
| 2969 | + "Other, and more advanced subsetting capabilities are also available. For more information and examples, consult those [subsetting examples](https://pavics.ouranos.ca/climate_analysis.html#b). " |
2970 | 2970 | ]
|
2971 | 2971 | },
|
2972 | 2972 | {
|
|
0 commit comments