Subset and download a specific marine product from Copernicus
Source:R/cms_download_subset.r
cms_download_subset.Rd
Subset and download a specific marine product from Copernicus. You need to register an account at https://data.marine.copernicus.eu before you can use this function.
Arguments
- username
Your Copernicus marine user name. Can be provided as
options(CopernicusMarine_uid = "my_user_name")
, or as argument here.- password
Your Copernicus marine password. Can be provided as
options(CopernicusMarine_pwd = "my_password")
, or as argument here.- destination
File or path where the requested file will be downloaded to.
- product
An identifier (type
character
) of the desired Copernicus marine product. Can be obtained withcms_products_list
.- layer
The name of a desired layer within a product (type
character
). Can be obtained withcms_product_details
.- variable
The name of a desired variable in a specific layer of a product (type
character
). Can be obtained withcms_product_details
.- region
Specification of the bounding box as a
vector
ofnumeric
s WGS84 lat and lon coordinates. Should be in the order of: xmin, ymin, xmax, ymax.- timerange
A
vector
with two elements (lower and upper value) for a requested time range. Thevector
should be coercible toPOSIXct
.- verticalrange
A
vector
with two elements (minimum and maximum) numerical values for the depth of the vertical layers (if any). Note that values below the sea surface needs to be specified as negative values.- overwrite
A
logical
value. WhenFALSE
(default), files at thedestination
won't be overwritten when the exist. Instead an error will be thrown if this is the case. When set toTRUE
, existing files will be overwritten.
Value
Returns a logical
value invisibly indicating whether the requested file was
successfully stored at the destination
.
Examples
if (FALSE) { # \dontrun{
destination <- tempfile("copernicus", fileext = ".nc")
## Assuming that Copernicus account details are provided as `options`
cms_download_subset(
destination = destination,
product = "GLOBAL_ANALYSISFORECAST_PHY_001_024",
layer = "cmems_mod_glo_phy-cur_anfc_0.083deg_P1D-m",
variable = "sea_water_velocity",
region = c(-1, 50, 10, 55),
timerange = c("2021-01-01 UTC", "2021-01-02 UTC"),
verticalrange = c(0, -2)
)
mydata <- stars::read_stars(destination)
plot(mydata["vo"])
} # }