Description Usage Arguments Details Value Examples
View source: R/get_ensemble_temp.R
Function wraps get_ensemble_climate_data
and returns
precipitation by basin or country in mm. Output is the 10th 50th and 90th percentile
for all gcm's for the a1 and b2 scenarios.
1 |
locator |
A vector of either watershed basin ID's from http://data.worldbank.org/sites/default/files/climate_data_api_basins.pdf It can be just a single basin id, or a vector of ids. ids should be strings. |
type |
the type of data to retrieve, must be "mavg" for monthly averages, "annualavg" for annual averages, "manom" for monthly anomaly, and "annualanom" for annual anomaly. |
start |
the start year to gather data for. |
end |
the end year to gather data to. |
start and end year can be any years, but all years will be coerced into periods outlined by the API (http://data.worldbank.org/developers/climate-data-api) anomaly periods are only valid for future scenarios and based on a reference period of 1969 - 1999, see API for full details.
a dataframe with precipitation predictions in mm for all scenarios, gcms, for each time period.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | ## Not run:
# Get data for 2 basins, annual average precipitation for all valid time periods
# then subset them, and plot
temp_dat <- get_ensemble_temp(locator=c(2,231), type="annualavg", start=1900, end=3000)
temp_dat <- subset(temp_dat,temp_dat$scenario!="b1")
temp_dat$uniqueGroup <- paste(temp_dat$percentile,temp_dat$locator,sep="-")
ggplot(temp_dat, aes(x=fromYear, y=data, group=uniqueGroup,
colour=as.factor(locator), linetype=as.factor(percentile))) +
geom_path()
### Get data for 2 countries with monthly precipitation values
temp_dat <- get_ensemble_temp(locator = c("USA","BRA"), type = "mavg", start = 2020, end = 2030)
temp_dat <- subset(temp_dat,temp_dat$scenario!="b1")
temp_dat$uniqueGroup <- paste(temp_dat$percentile, temp_dat$locator,sep="-")
ggplot(temp_dat, aes(x=as.factor(month), y=data, group=uniqueGroup, colour=locator)) +
geom_path()
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.