nasapower (version 1.1.3)

get_power: Get NASA POWER data

Description

Get POWER global meteorology and surface solar energy climatology data and return a tidy data frame tibble. All options offered by the official POWER API are supported.

Usage

get_power(community, pars, temporal_average, lonlat, dates = NULL)

Value

A data frame of POWER data including location, dates (not including “CLIMATOLOGY”) and requested parameters. A header of metadata is included.

Argument details for “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.

AG

Provides access to the Agroclimatology Archive, which contains industry-friendly parameters formatted for input to crop models.

SB

Provides access to the Sustainable Buildings Archive, which contains industry-friendly parameters for the buildings community to include parameters in multi-year monthly averages.

SSE

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.

Argument details for <code>temporal_average</code>

There are three valid values.

DAILY

The daily average of pars by day, month and year.

INTERANNUAL

The monthly average of pars by year.

CLIMATOLOGY

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).

Argument details for <code>lonlat</code>

For a single point

To get a specific cell, 1/2 x 1/2 degree, supply a length-two numeric vector giving the decimal degree longitude and latitude in that order for data to download, e.g., lonlat = c(-89.5, -179.5).

For regional coverage

To get a region, supply a length-four 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).

For global coverage

To get global coverage for CLIMATOLOGY, supply “GLOBAL” while also specifying “CLIMATOLOGY” for the temporal_average.

Argument details for <code>dates</code>

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”.

References

https://power.larc.nasa.gov/documents/POWER_Data_v9_methodology.pdf https://power.larc.nasa.gov

Examples

# NOT RUN {
# Fetch daily "AG" community temperature, relative
# humidity and precipitation for January 1 1985
# for Kingsthorpe, Queensland, Australia
ag_d <- get_power(
  community = "AG",
  lonlat = c(151.81, -27.48),
  pars = c("RH2M", "T2M", "PRECTOT"),
  dates = "1985-01-01",
  temporal_average = "DAILY"
)

# Fetch single point climatology for air temperature
ag_c_point <- get_power(
  community = "AG",
  pars = "T2M",
  c(151.81, -27.48),
  temporal_average = "CLIMATOLOGY"
)

ag_c_point

# Fetch global AG climatology for air temperature
ag_c_global <- get_power(
  community = "AG",
  pars = "T2M",
  lonlat = "GLOBAL",
  temporal_average = "CLIMATOLOGY"
)

ag_c_global

# Fetch interannual solar cooking parameters
# for a given region
sse_i <- 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")
)

sse_i

# }
# NOT RUN {
# }