R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Milankovitch graph for a given astronomical configuration

Milankovitch

R Documentation

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 radiants

Value

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

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)