Computes caloric summer insolation for a given astronomical configuration and latitude.

1 |

`orbit` |
Output from a solution, such as |

`lat` |
latitude |

`...` |
Other arguments passed to Insol |

The caloric summer is a notion introduced by M. Milankovitch. It is defined as the halve of the tropical year during for which daily mean insolation are greater than all days of the other halves. The algorithm is an original algorithm by M. Crucifix, but consistent with earlier definitions and algorithms by A. Berger (see examples). Do not confuse this Berger (1978) reference with the Berger (1978), J. Atm. Sci. of the astronomical solution.

Time-integrated insolation in kJ/m2 during the caloric summer.

Michel Crucifix, U. catholique de Louvain, Belgium.

Berger (1978) Long-term variations of caloric insolation resulting from the earth's orbital elements, Quaternary Research, 9, 139 - 167.

1 2 3 4 5 6 7 8 9 10 11 12 13 | ```
## reproduces Table 2 of Berger 1978
lat <- seq(90, 0, -10) * pi/180. ## angles in radiants.
orbit_1 = ber78(0)
orbit_2 = orbit_1
orbit_2 ['eps'] = orbit_2['eps'] + 1*pi/180.
T <- sapply(lat, function(x) c(lat = x * 180/pi,
calins(orbit_2, lat=x, S0=1365) / (4.18 * 1e1)
- calins(orbit_1, lat=x, S0=1365) / (4.18 * 1e1) ) )
data.frame(t(T))
# there are still some differences, of the order of 0.3 %, that are probably related to
# the slightly different methods.
# 41.8 is the factor from cal/cm2 to kJ/m2
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.