palinsol (version 0.92)

Insol: Computes incoming solar radiation (insolation)

Description

Computes incoming solar radiation (insolation) for a given astronomical configuration, true solar longitude and latitude

Usage

Insol (orbit,long=pi/2, lat=65*pi/180,S0=1365)

Arguments

orbit

Output from a solution, such as ber78, ber90 or la04

long

true solar longitude

lat

latitude

S0

Total solar irradiance

Value

• Daily-mean insolation (assuming fixed astronomical parameters during a true solar day)

Details

True solar longitude is measured in radiants: ll{ pi/2 for summer solstice pi for autumn equinox 3 * pi/2 for winter solstice 0 for spring equinox } It may be obtained for a given day in the year using the function day2l.

References

Berger, A. L. (1978). Long-term variations of daily insolation and Quaternary climatic changes, J. Atmos. Sci., 35, 2362-2367.

Examples

## make a little wrapper, with all default values

insolation <- function(times, astrosol=ber78,...)
  sapply(times, function(tt) Insol(orbit=astrosol(tt)))

tts <- seq(from = -400e3, to = 0, by = 1e3)
isl <- insolation(tts, ber78)
plot(tts, isl, typ='l')