nasapower (version 1.0.1)

get_power: Get POWER Data and Return a Tidy Data Frame

Description

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.

Usage

get_power(community, pars, temporal_average, lonlat = NULL,
  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 <code>community</code>

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

For regional coverage

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

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_v8_methodology.pdf

Examples

# NOT RUN {
# 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")
                                      )
# }
# NOT RUN {
# }