Solar astronomy using Meeus' algorithm
The exact julian day computation is adapted from ode::julianDay() and tested againts test cases in Redas and Andreas (2008, table A4.1) for validity up to year 4712 BC.
calendar_change julian_day(time) julian_day_fast(time) julian_century(time) geom_mean_lon_sun(x) geom_mean_anom_sun(x) eccent_earth_orbit(x) sun_eq_of_ctr(x, anom) sun_rad_vector(eccent, anom) sun_app_lon(x, lon) mean_obliq_eclip(x) obliq_corr(x, eclip) sun_rt_ascen(app.lon, obliq.corr) sun_decline(app.lon, obliq.corr) var_y(obliq.corr) eq_of_time(mean.lon, eccent.earth, anom.mean, var.y) ha_sunrise(lat, decline, nag = 0) solar_noon(lon, eq.of.time) sunrise(noon, ha.sunrise) sunset(noon, ha.sunrise) sunlight_duration(ha.sunrise, unit.out = "hours") solar_datetime(time, lat, lon, eq.of.time) solar_tod(time, lat, lon, eq.of.time) hour_angle(solar.time) zenith_angle(lat, hour.angle, decline) elevation_angle(lat, hour.angle, decline) atm_refraction_approx(elevation.angle) azimuth_angle(lat, hour.angle, zenith.angle, decline)
time |
dateTime |
x |
numeric Julian century |
anom |
numeric Solar anomaly in degrees |
eccent |
numeric Eccentricity of Earth orbit |
eclip |
numeric Ecliptic |
app.lon, obliq.corr, mean.lon, nag, decline |
numeric Angles in degrees |
lat, lon |
numeric Geographic coordinates in degrees |
eq.of.time, ha.sunrise, noon |
numeric |
zenith.angle, elevation.angle, hour.angle |
numeric Angles in degrees |
An object of class POSIXct
(inherits from POSIXt
) of length 1.
Low level functions based on NOAA's Excel worksheet
datetime
numeric
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.