de.kah2.zodiac.nova4jmt

Class EllipticMotion

java.lang.Object

de.kah2.zodiac.nova4jmt.EllipticMotion

public class EllipticMotion
extends Object

Constructor Summary

Constructors

Constructor and Description
EllipticMotion()

Method Summary

Modifier and Type	Method and Description
static double	ln_get_ell_body_earth_dist(double JD, LnEllOrbit orbit) double ln_get_ell_body_earth_dist(double JD, LnEllOrbit orbit) \param JD Julian day.
static double	ln_get_ell_body_elong(double JD, LnEllOrbit orbit) double ln_get_ell_body_elong(double JD, LnEllOrbit orbit); \param JD Julian day \param orbit Orbital parameters \return Elongation to the Sun.
static void	ln_get_ell_body_equ_coords(double JD, LnEllOrbit orbit, LnEquPosn posn) void ln_get_ell_body_equ_coords(double JD, LnEllOrbit orbit, LnEquPosn posn) \param JD Julian Day.
static int	ln_get_ell_body_next_rst_horizon_future(double JD, LnLnlatPosn observer, LnEllOrbit orbit, double horizon, int day_limit, LnRstTime rst) double ln_get_ell_body_next_rst_horizon(double JD, LnLnlatPosn observer, LnEllOrbit orbit, double horizon, LnRstTime rst); \param JD Julian day \param observer Observers position \param orbit Orbital parameters \param horizon Horizon height \param day_limit Maximal number of days that will be searched for next rise and set \param rst Pointer to store Rise, Set and Transit time in JD \return 0 for success, else 1 for circumpolar (above the horizon), -1 for circumpolar (bellow the horizon) Calculate the time of next rise, set and transit (crosses the local meridian at upper culmination) time of a body with an elliptic orbit for the given Julian day.
static int	ln_get_ell_body_next_rst_horizon(double JD, LnLnlatPosn observer, LnEllOrbit orbit, double horizon, LnRstTime rst) double ln_get_ell_body_next_rst_horizon(double JD, LnLnlatPosn observer, LnEllOrbit orbit, double horizon, LnRstTime rst); \param JD Julian day \param observer Observers position \param orbit Orbital parameters \param horizon Horizon height \param rst Pointer to store Rise, Set and Transit time in JD \return 0 for success, else 1 for circumpolar (above the horizon), -1 for circumpolar (bellow the horizon) Calculate the time of next rise, set and transit (crosses the local meridian at upper culmination) time of a body with an elliptic orbit for the given Julian day.
static int	ln_get_ell_body_next_rst(double JD, LnLnlatPosn observer, LnEllOrbit orbit, LnRstTime rst) double ln_get_ell_body_next_rst(double JD, LnLnlatPosn observer, LnEllOrbit orbit, LnRstTime rst); \param JD Julian day \param observer Observers position \param orbit Orbital parameters \param rst Pointer to store Rise, Set and Transit time in JD \return 0 for success, else 1 for circumpolar (above the horizon), -1 for circumpolar (bellow the horizon) Calculate the time of next rise, set and transit (crosses the local meridian at upper culmination) time of a body with an elliptic orbit for the given Julian day.
static double	ln_get_ell_body_phase_angle(double JD, LnEllOrbit orbit) double ln_get_ell_body_phase_angle(double JD, LnEllOrbit orbit); \param JD Julian day \param orbit Orbital parameters \return Phase angle.
static int	ln_get_ell_body_rst_horizon(double JD, LnLnlatPosn observer, LnEllOrbit orbit, double horizon, LnRstTime rst) double ln_get_ell_body_rst_horizon(double JD, LnLnlatPosn observer, LnEllOrbit orbit, double horizon, LnRstTime rst); \param JD Julian day \param observer Observers position \param orbit Orbital parameters \param horizon Horizon height \param rst Pointer to store Rise, Set and Transit time in JD \return 0 for success, else 1 for circumpolar (above the horizon), -1 for circumpolar (bellow the horizon) Calculate the time the rise, set and transit (crosses the local meridian at upper culmination) time of a body with an elliptic orbit for the given Julian day.
static int	ln_get_ell_body_rst(double JD, LnLnlatPosn observer, LnEllOrbit orbit, LnRstTime rst) double ln_get_ell_body_rst(double JD, LnLnlatPosn observer, LnEllOrbit orbit, LnRstTime rst); \param JD Julian day \param observer Observers position \param orbit Orbital parameters \param rst Pointer to store Rise, Set and Transit time in JD \return 0 for success, else 1 for circumpolar (above the horizon), -1 for circumpolar (bellow the horizon) Calculate the time the rise, set and transit (crosses the local meridian at upper culmination) time of a body with an elliptic orbit for the given Julian day.
static double	ln_get_ell_body_solar_dist(double JD, LnEllOrbit orbit) double ln_get_ell_body_solar_dist(double JD, LnEllOrbit orbit) \param JD Julian Day.
static void	ln_get_ell_geo_rect_posn(LnEllOrbit orbit, double JD, LnRectPosn posn) void ln_get_ell_geo_rect_posn(LnEllOrbit orbit, double JD, LnRectPosn posn); \param orbit Orbital parameters of object.
static void	ln_get_ell_helio_rect_posn(LnEllOrbit orbit, double JD, LnRectPosn posn) void ln_get_ell_helio_rect_posn(LnEllOrbit orbit, double JD, LnRectPosn posn); \param orbit Orbital parameters of object.
static double	ln_get_ell_last_perihelion(double epoch_JD, double M, double n) !
static double	ln_get_ell_mean_anomaly(double n, double delta_JD) double ln_get_ell_mean_anomaly (double n, double delta_JD); \param n Mean motion (degrees/day) \param delta_JD Time since perihelion \return Mean anomaly (degrees) Calculate the mean anomaly.
static double	ln_get_ell_mean_motion(double a) double ln_get_ell_mean_motion (double a); \param a Semi major diameter in AU \return Mean daily motion (degrees/day) Calculate the mean daily motion (degrees/day).
static double	ln_get_ell_orbit_avel(LnEllOrbit orbit) double ln_get_ell_orbit_avel(LnEllOrbit orbit); \param orbit Orbital parameters \return Orbital velocity in km/s.
static double	ln_get_ell_orbit_len(LnEllOrbit orbit) double ln_get_ell_orbit_len(LnEllOrbit orbit); \param orbit Orbital parameters \return Orbital length in AU Calculate the orbital length in AU.
static double	ln_get_ell_orbit_pvel(LnEllOrbit orbit) double ln_get_ell_orbit_pvel(LnEllOrbit orbit); \param orbit Orbital parameters \return Orbital velocity in km/s.
static double	ln_get_ell_orbit_vel(double JD, LnEllOrbit orbit) double ln_get_ell_orbit_vel(double JD, LnEllOrbit orbit); \param JD Julian day.
static double	ln_get_ell_radius_vector(double a, double e, double E) equ 30.2
static double	ln_get_ell_smajor_diam(double e, double q) double ln_get_ell_smajor_diam (double e, double q); \param e Eccentricity \param q Perihelion distance in AU \return Semi-major diameter in AU Calculate the semi major diameter.
static double	ln_get_ell_sminor_diam(double e, double a) double ln_get_ell_sminor_diam (double e, double a); \param e Eccentricity \param a Semi-Major diameter in AU \return Semi-minor diameter in AU Calculate the semi minor diameter.
static double	ln_get_ell_true_anomaly(double e, double E) equ 30.1
static double	ln_solve_kepler(double e, double M) double ln_solve_kepler (double E, double M); \param E Orbital eccentricity \param M Mean anomaly \return Eccentric anomaly Calculate the eccentric anomaly.
static double	sgn(double x) the BASIC SGN() function for doubles

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Constructor Detail

EllipticMotion

public EllipticMotion()

Method Detail

sgn

public static double sgn(double x)

the BASIC SGN() function for doubles

ln_solve_kepler

public static double ln_solve_kepler(double e,
                                     double M)

double ln_solve_kepler (double E, double M); \param E Orbital eccentricity \param M Mean anomaly \return Eccentric anomaly Calculate the eccentric anomaly. This method was devised by Roger Sinnott. (Sky and Telescope, Vol 70, pg 159)

ln_get_ell_mean_anomaly

public static double ln_get_ell_mean_anomaly(double n,
                                             double delta_JD)

double ln_get_ell_mean_anomaly (double n, double delta_JD); \param n Mean motion (degrees/day) \param delta_JD Time since perihelion \return Mean anomaly (degrees) Calculate the mean anomaly.

ln_get_ell_true_anomaly

public static double ln_get_ell_true_anomaly(double e,
                                             double E)

equ 30.1

ln_get_ell_radius_vector

public static double ln_get_ell_radius_vector(double a,
                                              double e,
                                              double E)

equ 30.2

ln_get_ell_smajor_diam

public static double ln_get_ell_smajor_diam(double e,
                                            double q)

double ln_get_ell_smajor_diam (double e, double q); \param e Eccentricity \param q Perihelion distance in AU \return Semi-major diameter in AU Calculate the semi major diameter.

ln_get_ell_sminor_diam

public static double ln_get_ell_sminor_diam(double e,
                                            double a)

double ln_get_ell_sminor_diam (double e, double a); \param e Eccentricity \param a Semi-Major diameter in AU \return Semi-minor diameter in AU Calculate the semi minor diameter.

ln_get_ell_mean_motion

public static double ln_get_ell_mean_motion(double a)

double ln_get_ell_mean_motion (double a); \param a Semi major diameter in AU \return Mean daily motion (degrees/day) Calculate the mean daily motion (degrees/day).

ln_get_ell_helio_rect_posn

public static void ln_get_ell_helio_rect_posn(LnEllOrbit orbit,
                                              double JD,
                                              LnRectPosn posn)

void ln_get_ell_helio_rect_posn(LnEllOrbit orbit, double JD, LnRectPosn posn); \param orbit Orbital parameters of object. \param JD Julian day \param posn Position pointer to store objects position Calculate the objects rectangular heliocentric position given it's orbital elements for the given julian day.

ln_get_ell_geo_rect_posn

public static void ln_get_ell_geo_rect_posn(LnEllOrbit orbit,
                                            double JD,
                                            LnRectPosn posn)

void ln_get_ell_geo_rect_posn(LnEllOrbit orbit, double JD, LnRectPosn posn); \param orbit Orbital parameters of object. \param JD Julian day \param posn Position pointer to store objects position Calculate the objects rectangular geocentric position given it's orbital elements for the given julian day.

ln_get_ell_body_equ_coords

public static void ln_get_ell_body_equ_coords(double JD,
                                              LnEllOrbit orbit,
                                              LnEquPosn posn)

void ln_get_ell_body_equ_coords(double JD, LnEllOrbit orbit, LnEquPosn posn) \param JD Julian Day. \param orbit Orbital parameters. \param posn Pointer to hold asteroid position. Calculate a bodies equatorial coordinates for the given julian day.

ln_get_ell_orbit_len

public static double ln_get_ell_orbit_len(LnEllOrbit orbit)

double ln_get_ell_orbit_len(LnEllOrbit orbit); \param orbit Orbital parameters \return Orbital length in AU Calculate the orbital length in AU. Accuracy: - 0.001% for e < 0.88 - 0.01% for e < 0.95 - 1% for e = 0.9997 - 3% for e = 1

ln_get_ell_orbit_vel

public static double ln_get_ell_orbit_vel(double JD,
                                          LnEllOrbit orbit)

double ln_get_ell_orbit_vel(double JD, LnEllOrbit orbit); \param JD Julian day. \param orbit Orbital parameters \return Orbital velocity in km/s. Calculate orbital velocity in km/s for the given julian day.

ln_get_ell_orbit_pvel

public static double ln_get_ell_orbit_pvel(LnEllOrbit orbit)

double ln_get_ell_orbit_pvel(LnEllOrbit orbit); \param orbit Orbital parameters \return Orbital velocity in km/s. Calculate orbital velocity at perihelion in km/s.

ln_get_ell_orbit_avel

public static double ln_get_ell_orbit_avel(LnEllOrbit orbit)

double ln_get_ell_orbit_avel(LnEllOrbit orbit); \param orbit Orbital parameters \return Orbital velocity in km/s. Calculate the orbital velocity at aphelion in km/s.

ln_get_ell_body_solar_dist

public static double ln_get_ell_body_solar_dist(double JD,
                                                LnEllOrbit orbit)

double ln_get_ell_body_solar_dist(double JD, LnEllOrbit orbit) \param JD Julian Day. \param orbit Orbital parameters \return The distance in AU between the Sun and the body. Calculate the distance between a body and the Sun.

ln_get_ell_body_earth_dist

public static double ln_get_ell_body_earth_dist(double JD,
                                                LnEllOrbit orbit)

double ln_get_ell_body_earth_dist(double JD, LnEllOrbit orbit) \param JD Julian day. \param orbit Orbital parameters \returns Distance in AU Calculate the distance between an body and the Earth for the given julian day.

ln_get_ell_body_phase_angle

public static double ln_get_ell_body_phase_angle(double JD,
                                                 LnEllOrbit orbit)

double ln_get_ell_body_phase_angle(double JD, LnEllOrbit orbit); \param JD Julian day \param orbit Orbital parameters \return Phase angle. Calculate the phase angle of the body. The angle Sun - body - Earth.

ln_get_ell_body_elong

public static double ln_get_ell_body_elong(double JD,
                                           LnEllOrbit orbit)

double ln_get_ell_body_elong(double JD, LnEllOrbit orbit); \param JD Julian day \param orbit Orbital parameters \return Elongation to the Sun. Calculate the bodies elongation to the Sun..

ln_get_ell_body_rst

public static int ln_get_ell_body_rst(double JD,
                                      LnLnlatPosn observer,
                                      LnEllOrbit orbit,
                                      LnRstTime rst)

double ln_get_ell_body_rst(double JD, LnLnlatPosn observer, LnEllOrbit orbit, LnRstTime rst); \param JD Julian day \param observer Observers position \param orbit Orbital parameters \param rst Pointer to store Rise, Set and Transit time in JD \return 0 for success, else 1 for circumpolar (above the horizon), -1 for circumpolar (bellow the horizon) Calculate the time the rise, set and transit (crosses the local meridian at upper culmination) time of a body with an elliptic orbit for the given Julian day. Note: this functions returns 1 if the body is circumpolar, that is it remains the whole day above the horizon. Returns -1 when it remains the whole day below the horizon.

ln_get_ell_body_rst_horizon

public static int ln_get_ell_body_rst_horizon(double JD,
                                              LnLnlatPosn observer,
                                              LnEllOrbit orbit,
                                              double horizon,
                                              LnRstTime rst)

double ln_get_ell_body_rst_horizon(double JD, LnLnlatPosn observer, LnEllOrbit orbit, double horizon, LnRstTime rst); \param JD Julian day \param observer Observers position \param orbit Orbital parameters \param horizon Horizon height \param rst Pointer to store Rise, Set and Transit time in JD \return 0 for success, else 1 for circumpolar (above the horizon), -1 for circumpolar (bellow the horizon) Calculate the time the rise, set and transit (crosses the local meridian at upper culmination) time of a body with an elliptic orbit for the given Julian day. Note: this functions returns 1 if the body is circumpolar, that is it remains the whole day above the horizon. Returns -1 when it remains the whole day below the horizon.

ln_get_ell_body_next_rst

public static int ln_get_ell_body_next_rst(double JD,
                                           LnLnlatPosn observer,
                                           LnEllOrbit orbit,
                                           LnRstTime rst)

double ln_get_ell_body_next_rst(double JD, LnLnlatPosn observer, LnEllOrbit orbit, LnRstTime rst); \param JD Julian day \param observer Observers position \param orbit Orbital parameters \param rst Pointer to store Rise, Set and Transit time in JD \return 0 for success, else 1 for circumpolar (above the horizon), -1 for circumpolar (bellow the horizon) Calculate the time of next rise, set and transit (crosses the local meridian at upper culmination) time of a body with an elliptic orbit for the given Julian day. This function guarantee, that rise, set and transit will be in <JD, JD+1> range. Note: this functions returns 1 if the body is circumpolar, that is it remains the whole day above the horizon. Returns -1 when it remains the whole day below the horizon.

ln_get_ell_body_next_rst_horizon

public static int ln_get_ell_body_next_rst_horizon(double JD,
                                                   LnLnlatPosn observer,
                                                   LnEllOrbit orbit,
                                                   double horizon,
                                                   LnRstTime rst)

double ln_get_ell_body_next_rst_horizon(double JD, LnLnlatPosn observer, LnEllOrbit orbit, double horizon, LnRstTime rst); \param JD Julian day \param observer Observers position \param orbit Orbital parameters \param horizon Horizon height \param rst Pointer to store Rise, Set and Transit time in JD \return 0 for success, else 1 for circumpolar (above the horizon), -1 for circumpolar (bellow the horizon) Calculate the time of next rise, set and transit (crosses the local meridian at upper culmination) time of a body with an elliptic orbit for the given Julian day. This function guarantee, that rise, set and transit will be in <JD, JD+1> range. Note: this functions returns 1 if the body is circumpolar, that is it remains the whole day above the horizon. Returns -1 when it remains the whole day below the horizon.

ln_get_ell_body_next_rst_horizon_future

public static int ln_get_ell_body_next_rst_horizon_future(double JD,
                                                          LnLnlatPosn observer,
                                                          LnEllOrbit orbit,
                                                          double horizon,
                                                          int day_limit,
                                                          LnRstTime rst)

double ln_get_ell_body_next_rst_horizon(double JD, LnLnlatPosn observer, LnEllOrbit orbit, double horizon, LnRstTime rst); \param JD Julian day \param observer Observers position \param orbit Orbital parameters \param horizon Horizon height \param day_limit Maximal number of days that will be searched for next rise and set \param rst Pointer to store Rise, Set and Transit time in JD \return 0 for success, else 1 for circumpolar (above the horizon), -1 for circumpolar (bellow the horizon) Calculate the time of next rise, set and transit (crosses the local meridian at upper culmination) time of a body with an elliptic orbit for the given Julian day. This function guarantee, that rise, set and transit will be in <JD, JD + day_limit> range. Note: this functions returns 1 if the body is circumpolar, that is it remains the whole day above the horizon. Returns -1 when it remains the whole day below the horizon.

ln_get_ell_last_perihelion

public static double ln_get_ell_last_perihelion(double epoch_JD,
                                                double M,
                                                double n)

!\fn double ln_get_ell_last_perihelion (double epoch_JD, double M, double n); \param epoch_JD Julian day of epoch \param M Mean anomaly \param n daily motion in degrees Calculate the julian day of the last perihelion.