#' Create a grid of sites
#'
#' Flexibly create a grid of sites for ensemble trajectory model runs.
#'
#' @param lat,lon Latitude and longitude values in decimal degrees for the
#' center point of the grid.
#' @param range The latitude and longitude range about the `grid_ref`.
#' @param division The division distances across the latitude and longitude
#' ranges.
#' @export
create_grid <- function(lat = 49.263,
lon = -123.250,
range = c(5, 5),
division = c(0.5, 0.5)) {
n_s_points <- range[1] / division[1]
w_e_points <- range[2] / division[2]
lat_vec <-
sort(
(lat + range[1]/2) - seq(0, range[1], division[1]),
decreasing = TRUE
)
lon_vec <- sort((lon + range[2]/2) - seq(0, range[2], division[2]))
for (i in seq_along(lat_vec)) {
for (j in seq_along(lon_vec)) {
if (i == 1 & j == 1) {
coords <- vector("list", 2)
names(coords) <- c("lat", "lon")
}
coords$lat[length(coords$lat) + 1] <- lat_vec[i]
coords$lon[length(coords$lon) + 1] <- lon_vec[j]
}
}
coords
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.