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

1
2
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

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