Get climate data

Get data from PAVICS THREDDS catalog

Usage

data_PAVICS(
  OPENDAP_url =
    "https://pavics.ouranos.ca/twitcher/ows/proxy/thredds/dodsC/datasets/simulations/RCM-CMIP6/CORDEX/NAM-12/day/NAM-12_NorESM2-MM_historical_r1i1p1f1_OURANOS_CRCM5_v1-r1_day_19500101-20141231.ncml",
  first_date = "1950-01-01",
  end_date = "1951-01-1",
  variable = "tas",
  geom = data_planning_areas(),
  maskGeom = TRUE
)

Arguments

OPENDAP_url

character string for an OPENDAP url from the THREDDS catalog. Defaults to a RCM-CMIP6 daily simulation

first_date

character string for the first date to be included in the data subset. Defaults to "1950-01-01"

end_date

character string for the last date to be included in the data subset. Defaults to "1950-01-01"

variable

character string for the variable to query from the data. Defaults to "tas" for Near-Surface Air Temperature. See the OPENDAP url for more options

geom

sf data.frame used to limit the spatial extent of the data subset (via bounding box). Defaults to the Scotian Shelf and Bay of Fundy planning area

maskGeom

logical defaults to TRUE to turn values outside the geom to NA instead of returning all the values that are downloaded as is.

Value

matrix

Examples

if (FALSE) { # \dontrun{
require(MarConsNetData)
require(dplyr)
# get the Scotian Shelf polygon
ssbof <- data_planning_areas()

# get CMIP6 data dor Near-Surface Air Temperature (tas)
climate <- data_PAVICS(OPENDAP_url = "https://pavics.ouranos.ca/twitcher/ows/proxy/thredds/dodsC/datasets/simulations/RCM-CMIP6/CORDEX/NAM-12/day/NAM-12_NorESM2-MM_historical_r1i1p1f1_OURANOS_CRCM5_v1-r1_day_19500101-20141231.ncml",
                       first_date = "1950-01-01",
                       end_date = "1950-03-1",
                       variable = "tas",
                       geom = ssbof,
                       maskGeom = TRUE)

# plot a single time point
climate |>
  slice("time",1) |>
  plot()

# calculate and plot monthly means
monthlyClimate <- aggregate(climate,
                            by = "1 month",
                            FUN = mean,
                            na.rm = TRUE)

plot(monthlyClimate)

} # }