thesky_asteroids Module Reference

Procedures for asteroids. More...

Functions/Subroutines
subroutine	asteroid_elements (tjm, as, hcr, om1, nu, cosi, sini, coso, sino)
	Calculate orbital elements for asteroids.
subroutine	asteroid_lbr (t, as, l, b, r)
	Calculate heliocentric ecliptical coordinates l,b,r for asteroid as at time t (Julian Millennia after 2000.0 in DT)
subroutine	asteroid_eq (t, as, l0, b0, r0, ra, dec, delta)
	Calculate geocentric(?) equatorial coordinates ra, dec, delta for asteroid as at time t (Jul. Millennia DT after 2000)
real(double) function	asteroid_magn (as, delta, r, pa)
	Calculate asteroid magnitude for asteroid AS.
subroutine	kepler (ecc, mea, eca)
	Solve Kepler's equation.

Detailed Description

Procedures for asteroids.

Function/Subroutine Documentation

asteroid_elements()

subroutine thesky_asteroids::asteroid_elements	(	real(double), intent(in)	tjm,
		integer, intent(in)	as,
		real(double), intent(out)	hcr,
		real(double), intent(out)	om1,
		real(double), intent(out)	nu,
		real(double), intent(out)	cosi,
		real(double), intent(out)	sini,
		real(double), intent(out)	coso,
		real(double), intent(out)	sino )

Calculate orbital elements for asteroids.

Parameters

tjm	Dynamical time in Julian Millennia after 2000.0
as	Asteroid ID
hcr	Heliocentric distance (AU) (output)
om1	Argument of perihelion (rad) (output)
nu	True anomaly (rad) (output)
cosi	cos(i): cosine of inclination (output)
sini	sin(i): sine of inclination (output)
coso	cos(Om): cosine of longitude of ascending node (output)
sino	sin(Om): sine of longitude of ascending node (output)

See also

Meeus, Astronomical Algorithms, 1998, Ch.33

Note

Used by asteroid_lbr() and asteroid_eq()

Definition at line 50 of file asteroids.f90.

References thesky_planetdata::asterelems, kepler(), and thesky_coordinates::precess_orb().

Referenced by asteroid_eq(), and asteroid_lbr().

Here is the call graph for this function:

asteroid_eq()

subroutine thesky_asteroids::asteroid_eq	(	real(double), intent(in)	t,
		integer, intent(in)	as,
		real(double), intent(in)	l0,
		real(double), intent(in)	b0,
		real(double), intent(in)	r0,
		real(double), intent(out)	ra,
		real(double), intent(out)	dec,
		real(double), intent(out)	delta )

Calculate geocentric(?) equatorial coordinates ra, dec, delta for asteroid as at time t (Jul. Millennia DT after 2000)

Parameters

t	Dynamical time in Julian Millennia after 2000.0
as	Asteroid ID
l0	Heliocentric longitude of Earth
b0	Heliocentric latitude of Earth
r0	Heliocentric distance of Earth
ra	Geocentric right ascension (rad) (output)
dec	Geocentric declination (rad) (output)
delta	Geocentric distance (AU) (output)

See also

Meeus, Astronomical Algorithms, 1998, Ch.33

Definition at line 164 of file asteroids.f90.

References asteroid_elements(), and thesky_coordinates::calcsunxyz().

Here is the call graph for this function:

asteroid_lbr()

subroutine thesky_asteroids::asteroid_lbr	(	real(double), intent(in)	t,
		integer, intent(in)	as,
		real(double), intent(out)	l,
		real(double), intent(out)	b,
		real(double), intent(out)	r )

Calculate heliocentric ecliptical coordinates l,b,r for asteroid as at time t (Julian Millennia after 2000.0 in DT)

Parameters

t	Dynamical time in Julian Millennia after 2000.0
as	Asteroid ID
l	Heliocentric ecliptical longitude (rad) (output)
b	Heliocentric ecliptical latitude (rad) (output)
r	Heliocentric distance (AU) (output)

See also

Meeus, Astronomical Algorithms, 1998, Ch.33

Definition at line 115 of file asteroids.f90.

References asteroid_elements().

Referenced by thesky_planets::planet_position().

Here is the call graph for this function:

asteroid_magn()

real(double) function thesky_asteroids::asteroid_magn	(	integer, intent(in)	as,
		real(double), intent(in)	delta,
		real(double), intent(in)	r,
		real(double), intent(in)	pa )

Calculate asteroid magnitude for asteroid AS.

Parameters

as	Asteroid number
delta	Distance to the Earth
r	Distance to the Sun
pa	Phase angle

Return values

asteroid_magn

Magnitude of the asteroid

See also

Meeus, Astronomical Algorithms, 1998, Ch.33

Definition at line 243 of file asteroids.f90.

References thesky_planetdata::asterelems, and asteroid_magn().

Referenced by asteroid_magn(), and thesky_planets::planet_position().

Here is the call graph for this function:

kepler()

subroutine thesky_asteroids::kepler	(	real(double), intent(in)	ecc,
		real(double), intent(in)	mea,
		real(double), intent(out)	eca )

Solve Kepler's equation.

Parameters

ecc	Eccentricity
mea	Mean anomaly
eca	Eccentric anomaly (output)

Solve Kepler's equation for given eccentricity ecc, mean anomaly mea and eccentric anomaly eca

See also

Meeus, Astronomical Algorithms, 1998, p.206

Definition at line 276 of file asteroids.f90.

Referenced by asteroid_elements().