Get POWER global meteorology and surface solar energy climatology data and return a tidy data frame. All options offered by the official POWER API are supported.
get_power(community, pars, temporal_average, lonlat = NULL, dates = NULL)
community | A character vector providing community name: "AG", "SB" or "SSE". See argument details for more. |
---|---|
pars | A character vector of solar, meteorological or climatology
parameters to download. See |
temporal_average | Temporal average for data being queried, supported values are DAILY, INTERANNUAL and CLIMATOLOGY. See argument details for more. |
lonlat | A numeric vector of geographic coordinates for a cell or region
entered as x, y coordinates. Not used when |
dates | A character vector of start and end dates in that order, |
A data frame of POWER data including location, dates (not including CLIMATOLOGY) and requested parameters. A header of metadata is included.
The associated metadata are not saved if the data are exported to a file format other than a native R data format, e.g., .Rdata, .rda or .rds.
community
There are three valid values, one must be supplied. This will affect the units of the parameter and the temporal display of time series data.
Provides access to the Agroclimatology Archive, which contains industry-friendly parameters formatted for input to crop models.
Provides access to the Sustainable Buildings Archive, which contains industry-friendly parameters for the buildings community to include parameters in multi-year monthly averages.
Provides access to the Renewable Energy Archive, which contains parameters specifically tailored to assist in the design of solar and wind powered renewable energy systems.
temporal_average
There are three valid values.
The daily average of pars
by day, month and year.
The monthly average of pars
by year.
The monthly average of pars
at the surface of the
earth for a given month, averaged for that month over the 30-year period
(Jan. 1984 - Dec. 2013).
lonlat
To get a specific cell, 1/2 x 1/2 degree, supply a
length-2 numeric vector giving the decimal degree longitude and latitude in
that order for data to download,
e.g., lonlat = c(-89.5, -179.5)
.
To get a region, supply a length-4 numeric
vector as lower left (lon, lat) and upper right (lon, lat) coordinates,
e.g., lonlat = c(xmin, ymin, xmax, ymax)
in that order for a given
region, e.g., a bounding box for the southwestern corner of Australia:
lonlat = c(112.5, -55.5, 115.5, -50.5)
.
*Maximum area processed is 4.5 x 4.5 degrees (100 points).
dates
If one date only is provided, it will
be treated as both the start date and the end date and only a single day's
values will be returned, e.g., dates = "1983-01-01"
. When
temporal_average
is set to "INTERANNUAL", use only two year values
(YYYY), e.g. dates = c(1983, 2010)
. This argument should not be used
when temporal_average
is set to "CLIMATOLOGY".
https://power.larc.nasa.gov/documents/POWER_Data_v8_methodology.pdf
# Fetch daily "AG" community temperature, relative humidity and precipitation # for January 1 1985 for Kingsthorpe, Queensland, Australia daily_ag <- get_power(community = "AG", lonlat = c(151.81, -27.48), pars = c("RH2M", "T2M", "PRECTOT"), dates = "1985-01-01", temporal_average = "DAILY" ) # Fetch global AG climatology for temperature, relative humidity and # precipitation climatology_ag <- get_power(community = "AG", pars = c("RH2M", "T2M", "PRECTOT"), temporal_average = "CLIMATOLOGY" ) # Fetch interannual solar cooking parameters for a given region interannual_sse <- get_power(community = "SSE", lonlat = c(112.5, -55.5, 115.5, -50.5), dates = c("1984", "1985"), temporal_average = "INTERANNUAL", pars = c("CLRSKY_SFC_SW_DWN", "ALLSKY_SFC_SW_DWN") )