ETsolradRasters {envirem}R Documentation

Extraterrestrial Solar Radiation


Generate monthly extraterrestrial solar radiation rasters.


ETsolradRasters(rasterTemplate, year, outputDir = NULL, ...)



any rasterLayer that can be used to extract extent, resolution, projection, etc.


The year solar radiation should be calculated for. See details.


destination directory for rasters, can be NULL


additional arguments passed to writeRaster


Given the latitude values of the cells found in the raster template and the year, monthly extraterrestrial solar radiation can be calculated, using the palinsol R package. year = 0 corresponds to 1950. Although the year can take on any value, it should match the time period of the other rasters that will be used for generating ENVIREM variables. Suggestions would be year = 40 for the present, year = -6000 for the mid Holocene, and year = -21500 for the LGM.

If you are having problems with this function and the rasterTemplate is not in long/lat, try with an unprojected long/lat raster.


If outputDir = NULL, a SpatRaster is returned. Otherwise, rasters are written to disk in the designated directory, and nothing is returned. Naming of the layers uses the tag specified via assignNames.


Pascal Title


J. Laskar et al., A long-term numerical solution for the insolation quantities of the Earth, Astron. Astroph., 428, 261-285 2004.


# Find example rasters
rasterFiles <- list.files(system.file('extdata', package='envirem'), full.names=TRUE)
env <- rast(rasterFiles)

# set aside a template raster
template <- env[[1]]

# generate solar radiation for the present
solrad <- ETsolradRasters(template, year = 40, outputDir = NULL)

[Package envirem version 3.0 Index]