Milankovitch {palinsol}

Milankovitch graph for a given astronomical configuration

Description

Computes the distrubition in latitude and longitude of incoming solar radiation, known as a Milankovitch graph, with possibility of plotting with a dedicated plot function

Usage

Milankovitch(
  orbit,
  S0 = 1365,
  lat = seq(-pi/2, pi/2, l = 73),
  long = seq(0, 2 * pi, l = 145),
  deg = TRUE
)

Arguments

orbit

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

S0

Total solar irradiance

lat

latitudes, passed as an array

long

true solar longitudes, passed as an array

deg

If true : the axes of the Milankovitch object are expressed in degrees. Inputs are always in radians

Value

A object of Milankovitch class, which may be plotted using the regular plot function

Note

The polar night option may not be bullet-proof for exotic obliquities

Author(s)

Michel Crucifix, U. catholique de Louvain, Belgium.

References

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

Examples

orbit <- c(eps=0.409214, ecc=0.01672393, varpi=4.92251)
M <- Milankovitch(orbit)
plot(M, plot=contour)
plot(M, plot=contour, month=FALSE)