![]() |
SuperNOVAS v1.5
The NOVAS C library, made better
|
Macros | |
#define | _GNU_SOURCE |
for strcasecmp() | |
Functions | |
short | make_cat_entry (const char *restrict name, const char *restrict catalog, long cat_num, double ra, double dec, double pm_ra, double pm_dec, double parallax, double rad_vel, cat_entry *source) |
int | make_cat_object (const cat_entry *star, object *source) |
int | make_cat_object_sys (const cat_entry *star, const char *restrict system, object *source) |
int | make_ephem_object (const char *name, long num, object *body) |
short | make_object (enum novas_object_type type, long number, const char *name, const cat_entry *star, object *source) |
int | make_orbital_object (const char *name, long num, const novas_orbital *orbit, object *body) |
int | make_planet (enum novas_planet num, object *restrict planet) |
int | make_redshifted_cat_entry (const char *name, double ra, double dec, double z, cat_entry *source) |
int | make_redshifted_object (const char *name, double ra, double dec, double z, object *source) |
int | make_redshifted_object_sys (const char *name, double ra, double dec, const char *restrict system, double z, object *source) |
void | novas_case_sensitive (int value) |
double | novas_helio_dist (double jd_tdb, const object *restrict source, double *restrict rate) |
int | novas_init_cat_entry (cat_entry *restrict source, const char *restrict name, double ra, double dec) |
enum novas_planet | novas_planet_for_name (const char *restrict name) |
int | novas_set_catalog (cat_entry *restrict source, const char *restrict catalog, long num) |
int | novas_set_distance (cat_entry *source, double parsecs) |
int | novas_set_lsr_vel (cat_entry *source, double epoch, double v_kms) |
int | novas_set_parallax (cat_entry *source, double mas) |
int | novas_set_proper_motion (cat_entry *source, double pm_ra, double pm_dec) |
int | novas_set_redshift (cat_entry *source, double z) |
int | novas_set_ssb_vel (cat_entry *source, double v_kms) |
double | novas_solar_power (double jd_tdb, const object *restrict source) |
int | proper_motion (double jd_tdb_in, const double *pos, const double *restrict vel, double jd_tdb_out, double *out) |
int | starvectors (const cat_entry *restrict star, double *restrict pos, double *restrict motion) |
short | transform_cat (enum novas_transform_type option, double jd_tt_in, const cat_entry *in, double jd_tt_out, const char *out_id, cat_entry *out) |
int | transform_hip (const cat_entry *hipparcos, cat_entry *hip_2000) |
Functions that define an observing target.
short make_cat_entry | ( | const char *restrict | name, |
const char *restrict | catalog, | ||
long | cat_num, | ||
double | ra, | ||
double | dec, | ||
double | pm_ra, | ||
double | pm_dec, | ||
double | parallax, | ||
double | rad_vel, | ||
cat_entry * | source ) |
Fully populates the data structure for a catalog source, such as a star, with all parameters provided at once.
Alternatively, you may use novas_init_cat_entry()
to initialize just with the name and R.A./Dec coordinates, and then add further information step-by-step as needed, icluding using alternative parameters (SSB vs. LSR velocity, vs. redshift; parallax vs. distance). The latter approach provides more flexibility, and will result in more readable code, which is also easier to debug.
name | Object name (less than SIZE_OF_OBJ_NAME in length). It may be NULL if not relevant. | |
catalog | Catalog identifier or epoch (less than SIZE_OF_CAT_NAME in length). E.g. 'HIP' for Hipparcos, 'TY2' for Tycho-2; or 'ICRS', 'B1950', 'J2000'. It may be NULL if not relevant. | |
cat_num | Object number in the catalog. It is not used internally, so you may set it to anything that is meaningful to you, or set it to 0 as a default. | |
ra | [h] Right ascension of the object. | |
dec | [deg] Declination of the object. | |
pm_ra | [mas/yr] Proper motion in right ascension. | |
pm_dec | [mas/yr] Proper motion in declination. | |
parallax | [mas] Parallax. | |
rad_vel | [km/s] Radial velocity relative to the Solar-System Barycenter (SSB). To convert velocities defined against the Local Standard of Rest (LSR), you may use novas_lsr_to_ssb_vel() to convert appropriately. Or, to convert from a redshift value, you might use novas_z2v() . | |
[out] | source | Pointer to data structure to populate. |
References novas_init_cat_entry(), novas_set_catalog(), novas_set_parallax(), novas_set_proper_motion(), novas_set_ssb_vel(), and rad_vel().
Populates and object data structure with the data for a catalog source. The astrometric parameters must be defined with ICRS. To create objects in other reference systems, use make_cat_object_sys()
instead.
star | Pointer to structure to populate with the catalog data for a celestial object located outside the solar system, specified with ICRS astrometric parameters. | |
[out] | source | Pointer to the celestial object data structure to be populated. |
References make_object(), NOVAS_CATALOG_OBJECT, cat_entry::starname, and cat_entry::starnumber.
Populates and object data structure with the data for a catalog source for a given system of catalog coordinates.
star | Pointer to structure to populate with the catalog data for a celestial object located outside the solar system. | |
system | Input catalog coordinate system epoch, e.g. "ICRS", "B1950.0", "J2000.0", "FK4", "FK5", or "HIP". In general, any Besselian or Julian year epoch can be used by year (e.g. "B1933.193" or "J2022.033"), or else the fixed value listed. If 'B' or 'J' is ommitted in front of the epoch year, then Besselian epochs are assumed prior to 1984.0. (See ‘novas_epoch() for more). | |
[out] | source | Pointer to the celestial object data structure to be populated with the corresponding ICRS catalog coordinates, after appying proper-motion and precession corrections as appropriate. |
References make_cat_object(), and object::star.
int make_ephem_object | ( | const char * | name, |
long | num, | ||
object * | body ) |
Sets a celestial object to be a Solar-system ephemeris body. Typically this would be used to define minor planets, asteroids, comets and planetary satellites.
name | Name of object. By default converted to upper-case, unless novas_case_sensitive() was called with a non-zero argument. Max. SIZE_OF_OBJ_NAME long, including termination. If the ephemeris provider uses names, then the name should match those of the ephemeris provider – otherwise it is not important. | |
num | Solar-system body ID number (e.g. NAIF). The number should match the needs of the ephemeris provider used with NOVAS. (If the ephemeris provider is by name and not ID number, then the number here is not important). | |
[out] | body | Pointer to structure to populate. |
References make_object(), and NOVAS_EPHEM_OBJECT.
short make_object | ( | enum novas_object_type | type, |
long | number, | ||
const char * | name, | ||
const cat_entry * | star, | ||
object * | source ) |
make_cat_object()
, make_redshifted_object()
, make_planet()
, make_ephem_object()
, or make_orbital_object()
.Populates an object data structure using the parameters provided. By default (for compatibility with NOVAS C) source names are converted to upper-case internally. You can however enable case-sensitive processing by calling novas_case_sensitive() before.
NOTES:
orbit
field (added in v1.2) with zeroes to remain ABI compatible with versions <1.2. type | The type of object. NOVAS_PLANET (0), NOVAS_EPHEM_OBJECT (1) or NOVAS_CATALOG_OBJECT (2), or NOVAS_ORBITAL_OBJECT (3). | |
number | The novas ID number (for solar-system bodies only, otherwise ignored) | |
name | The name of the object (case insensitive). It should be shorter than SIZE_OF_OBJ_NAME or else an error will be returned. The name is converted to upper internally unless novas_case_sensitive() was called before to change that. | |
star | Pointer to structure to populate with the catalog data for a celestial object located outside the solar system. Used only if type is NOVAS_CATALOG_OBJECT, otherwise ignored and can be NULL. | |
[out] | source | Pointer to the celestial object data structure to be populated. |
References object::name, NOVAS_CATALOG_OBJECT, NOVAS_OBJECT_TYPES, NOVAS_PLANET, NOVAS_PLANETS, object::number, object::star, and object::type.
int make_orbital_object | ( | const char * | name, |
long | num, | ||
const novas_orbital * | orbit, | ||
object * | body ) |
Sets a celestial object to be a Solar-system orbital body. Typically this would be used to define minor planets, asteroids, comets, or even planetary satellites.
name | Name of object. It may be NULL if not relevant. | |
num | Solar-system body ID number (e.g. NAIF). It is not required and can be set e.g. to -1 if not relevant to the caller. | |
orbit | The orbital parameters to adopt. The data will be copied, not referenced. | |
[out] | body | Pointer to structure to populate. |
References make_object(), NOVAS_ORBITAL_OBJECT, and object::orbit.
int make_planet | ( | enum novas_planet | num, |
object *restrict | planet ) |
Sets a celestial object to be a major planet, or the Sun, Moon, Solar-system Barycenter, etc.
num | Planet ID number (NOVAS convention) | |
[out] | planet | Pointer to structure to populate. |
References make_object(), NOVAS_PLANET, NOVAS_PLANET_NAMES_INIT, and NOVAS_PLANETS.
int make_redshifted_cat_entry | ( | const char * | name, |
double | ra, | ||
double | dec, | ||
double | z, | ||
cat_entry * | source ) |
Populates a celestial object data structure with the parameters for a redhifted catalog source, such as a distant quasar or galaxy. It is similar to make_cat_object()
except that it takes a Doppler-shift (z) instead of radial velocity and it assumes no parallax and no proper motion (appropriately for a distant redshifted source). The catalog name is set to EXT
to indicate an extragalactic source, and the catalog number defaults to 0. The user may change these default field values as appropriate afterwards, if necessary.
name | Object name (less than SIZE_OF_OBJ_NAME in length). It may be NULL. | |
ra | [h] Right ascension of the object (hours). | |
dec | [deg] Declination of the object (degrees). | |
z | Redhift value (λobs / λrest - 1 = frest / fobs - 1). | |
[out] | source | Pointer to structure to populate. |
References novas_init_cat_entry(), novas_set_catalog(), and novas_set_redshift().
int make_redshifted_object | ( | const char * | name, |
double | ra, | ||
double | dec, | ||
double | z, | ||
object * | source ) |
Populates a celestial object data structure with the ICRS astrometric parameters for a redhifted catalog source, such as a distant quasar or galaxy. To create redshifted objects with in other reference systems, use make_redshifted_object_sys()
instead.
It is similar to make_cat_object()
except that it takes a Doppler-shift (z) instead of radial velocity and it assumes no parallax and no proper motion (appropriately for a distant redshifted source). The catalog name is set to EXT
to indicate an extragalactic source, and the catalog number defaults to 0. The user may change these default field values as appropriate afterwards, if necessary.
name | Object name (less than SIZE_OF_OBJ_NAME in length). It may be NULL. | |
ra | [h] ICRS Right ascension of the object (hours). | |
dec | [deg] ICRS Declination of the object (degrees). | |
z | Redhift value (λobs / λrest - 1 = frest / fobs - 1). | |
[out] | source | Pointer to structure to populate. |
References make_cat_object(), and make_redshifted_cat_entry().
int make_redshifted_object_sys | ( | const char * | name, |
double | ra, | ||
double | dec, | ||
const char *restrict | system, | ||
double | z, | ||
object * | source ) |
Populates a celestial object data structure with the parameters for a redhifted catalog source, such as a distant quasar or galaxy, for a given system of catalog coordinates.
name | Object name (less than SIZE_OF_OBJ_NAME in length). It may be NULL. | |
ra | [h] ICRS Right ascension of the object (hours). | |
dec | [deg] ICRS Declination of the object (degrees). | |
system | Input catalog coordinate system epoch, e.g. "ICRS", "B1950.0", "J2000.0", "FK4", "FK5", or "HIP". In general, any Besselian or Julian year epoch can be used by year (e.g. "B1933.193" or "J2022.033"), or else the fixed value listed. If 'B' or 'J' is ommitted in front of the epoch year, then Besselian epochs are assumed prior to 1984.0. | |
z | Redhift value (λobs / λrest - 1 = frest / fobs - 1). | |
[out] | source | Pointer to the celestial object data structure to be populated with the corresponding ICRS catalog coordinates. |
References make_redshifted_object(), and object::star.
void novas_case_sensitive | ( | int | value | ) |
Enables or disables case-sensitive processing of the object name. The effect is not retroactive. The setting will only affect the celestial objects that are defined after the call. Note, that catalog names, set via make_cat_entry() are always case sensitive regardless of this setting.
value | (boolean) TRUE (non-zero) to enable case-sensitive object names, or else FALSE (0) to convert names to upper case only (NOVAS C compatible behavior). |
double novas_helio_dist | ( | double | jd_tdb, |
const object *restrict | source, | ||
double *restrict | rate ) |
Returns a Solar-system body's distance from the Sun, and optionally also the rate of recession. It may be useful, e.g. to calculate the body's heating from the Sun.
jd_tdb | [day] Barycentric Dynamical Time (TDB) based Julian date. You may want to use a time that is antedated to when the observed light originated from the source. | |
source | Observed Solar-system source | |
[out] | rate | [AU/day] (optional) Returned rate of recession from Sun |
References ephemeris(), NOVAS_CATALOG_OBJECT, NOVAS_HELIOCENTER, NOVAS_REDUCED_ACCURACY, and novas_vlen().
int novas_init_cat_entry | ( | cat_entry *restrict | source, |
const char *restrict | name, | ||
double | ra, | ||
double | dec ) |
Initializes a catalog entry, such as a star or a distant galaxy, with a name and coordinates. All other fields of the cat_entry
structure will be initialized with zeroes. After the initialization, you may set further properties, such as:
The typical process for setting up a catalog source is captured by the following steps:
cat_entry
with source name and RA / Dec coordinates using novas_init_cat_entry()
.novas_set_ssb_vel()
, ornovas_set_lsr_vel()
, ornovas_set_redshift()
.novas_set_proper_motion()
.novas_set_parallax()
, ornovas_set_distance()
.novas_set_catalog()
.After the initialization (Step 1), order does not matter.
[out] | source | Output structure to populate. |
name | Object name (less than SIZE_OF_OBJ_NAME in length). It may be NULL if not relevant. | |
ra | [h] Right ascension of the object. | |
dec | [deg] Declination of the object. |
References SIZE_OF_OBJ_NAME.
enum novas_planet novas_planet_for_name | ( | const char *restrict | name | ) |
Returns the NOVAS planet ID for a given name (case insensitive), or -1 if no match is found.
name | The planet name, or that for the "Sun", "Moon" or "SSB" (case insensitive). The spelled out "Solar System Barycenter" is also recognized with either spaces, hyphens ('-') or underscores ('_') separating the case insensitive words. |
References NOVAS_PLANET_NAMES_INIT, NOVAS_PLANETS, and NOVAS_SSB.
int novas_set_catalog | ( | cat_entry *restrict | source, |
const char *restrict | catalog, | ||
long | num ) |
Sets the optional catalog information for a sidereal source. SuperNOVAS does not use the catalog information internally. It is entirely for the convenience of the user to set these values if it is useful to them.
[out] | source | Output structure to populate with the parameters. |
catalog | Catalog identifier or epoch (less than SIZE_OF_CAT_NAME in length). E.g. 'HIP' for Hipparcos, 'TY2' for Tycho-2; or 'ICRS', 'B1950', 'J2000'. It may be NULL if not relevant. | |
num | Object's number in catalog. |
References SIZE_OF_CAT_NAME.
int novas_set_distance | ( | cat_entry * | source, |
double | parsecs ) |
Sets the distance for a catalog source.
[out] | source | Output structure to populate with the parameters. |
parsecs | [pc] distance |
References novas_set_parallax().
int novas_set_lsr_vel | ( | cat_entry * | source, |
double | epoch, | ||
double | v_kms ) |
Sets a radial velocity for a catalog source, defined w.r.t. the Local Standard of Rest (LSR).
[out] | source | Output structure to populate with the parameters. |
v_kms | [km/s] Radial velocity of source w.r.t. the Local Standard of Rest (LSR). | |
epoch | [yr] Coordinate epoch. |
References cat_entry::dec, novas_lsr_to_ssb_vel(), novas_set_ssb_vel(), and cat_entry::ra.
int novas_set_parallax | ( | cat_entry * | source, |
double | mas ) |
Sets the parallax (a measure of distance) for a catalog source.
[out] | source | Output structure to populate with the parameters. |
mas | [mas] Parallax |
References cat_entry::parallax.
int novas_set_proper_motion | ( | cat_entry * | source, |
double | pm_ra, | ||
double | pm_dec ) |
Sets the proper-motion for a (Galactic) catalog source.
[out] | source | Output structure to populate with the parameters. |
pm_ra | [mas/yr] Proper motion in right ascension. | |
pm_dec | [mas/yr] Proper motion in declination. |
References cat_entry::promodec, and cat_entry::promora.
int novas_set_redshift | ( | cat_entry * | source, |
double | z ) |
Sets a redhift for a catalog source, as a relativistic measure of velocity.
[out] | source | Output structure to populate with the parameters. |
z | [-1:inf] The redshift measure z, defined as (1 + z) = sqrt((1 + v/_c_) / (1 - v/_c_)) |
References novas_set_ssb_vel(), and novas_z2v().
int novas_set_ssb_vel | ( | cat_entry * | source, |
double | v_kms ) |
Sets a radial velocity for a catalog source, defined w.r.t. the Solar-System Barycenter (SSB).
[out] | source | Output structure to populate with the parameters. |
v_kms | [km/s] Radial velocity of source w.r.t. the Solar-System Barycenter (SSB). |
References NOVAS_C, and cat_entry::radialvelocity.
double novas_solar_power | ( | double | jd_tdb, |
const object *restrict | source ) |
Returns the typical incident Solar power on a Solar-system body at the time of observation.
jd_tdb | [day] Barycentric Dynamical Time (TDB) based Julian date. You may want to use a time that is antedated to when the observed light originated ( was reflected) from the source. |
source | Observed Solar-system source |
References novas_helio_dist(), and NOVAS_SOLAR_CONSTANT.
int proper_motion | ( | double | jd_tdb_in, |
const double * | pos, | ||
const double *restrict | vel, | ||
double | jd_tdb_out, | ||
double * | out ) |
Applies proper motion, including foreshortening effects, to a star's position.
REFERENCES:
jd_tdb_in | [day] Barycentric Dynamical Time (TDB) based Julian date of the first epoch. | |
pos | [AU] Position vector at first epoch. | |
vel | [AU/day] Velocity vector at first epoch. | |
jd_tdb_out | [day] Barycentric Dynamical Time (TDB) based Julian date of the second epoch. | |
[out] | out | Position vector at second epoch. It can be the same vector as the input. |
int starvectors | ( | const cat_entry *restrict | star, |
double *restrict | pos, | ||
double *restrict | motion ) |
Converts angular quantities for stars to vectors.
NOTES:
REFERENCES:
star | Pointer to catalog entry structure containing ICRS catalog | |
[out] | pos | [AU] Position vector, equatorial rectangular coordinates, It may be NULL if not required. |
[out] | motion | [AU/day] Perceived motion of star, in equatorial rectangular coordinates. It must be distinct from the pos output vector, and may be NULL if not required. Note, that it is suitable only for calculating the apparent 3D location of the star at a different time, and should not be used as a measure of physical velocity, e.g. for spectroscopic radial velocity determination. |
References NOVAS_KMS, novas_los_to_xyz(), and radec2vector().
short transform_cat | ( | enum novas_transform_type | option, |
double | jd_tt_in, | ||
const cat_entry * | in, | ||
double | jd_tt_out, | ||
const char * | out_id, | ||
cat_entry * | out ) |
Transform a star's catalog quantities for a change the coordinate system and/or the date for which the positions are calculated. Also used to rotate catalog quantities on the dynamical equator and equinox of J2000.0 to the ICRS or vice versa.
'date_incat' and 'date_newcat' may be specified either as a Julian date (e.g., 2433282.5 or NOVAS_JD_B1950) or a fractional Julian year and fraction (e.g., 1950.0). Values less than 10000 are assumed to be years. You can also use the supplied constants NOVAS_JD_J2000 or NOVAS_JD_B1950. The date arguments are ignored for the ICRS frame conversion options.
If 'option' is PROPER_MOTION (1), input data can be in any reference system. If 'option' is PRECESSION (2) or CHANGE_EPOCH (3), input data is assume to be in the dynamical system of 'date_incat' and produces output in the dynamical system of 'date_outcat'. If 'option' is CHANGE_J2000_TO_ICRS (4), the input data should be in the J2000.0 dynamical frame. And if 'option' is CHANGE_ICRS_TO_J2000 (5), the input data must be in the ICRS, and the output will be in the J2000 dynamical frame.
This function cannot be properly used to bring data from old star catalogs into the modern system, because old catalogs were compiled using a set of constants that are incompatible with modern values. In particular, it should not be used for catalogs whose positions and proper motions were derived by assuming a precession constant significantly different from the value implicit in function precession().
option | Type of transformation | |
jd_tt_in | [day|yr] Terrestrial Time (TT) based Julian date, or year, of input catalog data. Not used if option is CHANGE_J2000_TO_ICRS (4) or CHANGE_ICRS_TO_J2000 (5). | |
in | An entry from the input catalog, with units as given in the struct definition | |
jd_tt_out | [day|yr] Terrestrial Time (TT) based Julian date, or year, of output catalog data. Not used if option is CHANGE_J2000_TO_ICRS (4) or CHANGE_ICRS_TO_J2000 (5). | |
out_id | Catalog identifier (0 terminated). It may also be NULL in which case the catalog name is inherited from the input. | |
[out] | out | The transformed catalog entry, with units as given in the struct definition. It may be the same as the input. |
References cat_entry::catalog, CHANGE_EPOCH, CHANGE_ICRS_TO_J2000, CHANGE_J2000_TO_ICRS, cat_entry::dec, frame_tie(), ICRS_TO_J2000, J2000_TO_ICRS, NOVAS_JD_J2000, NOVAS_KMS, novas_los_to_xyz(), novas_vlen(), novas_xyz_to_los(), cat_entry::parallax, PRECESSION, precession(), cat_entry::promodec, cat_entry::promora, PROPER_MOTION, cat_entry::ra, radec2vector(), cat_entry::radialvelocity, SIZE_OF_CAT_NAME, SIZE_OF_OBJ_NAME, cat_entry::starname, and cat_entry::starnumber.
Convert Hipparcos catalog data at epoch J1991.25 to epoch J2000.0, for use within NOVAS. To be used only for Hipparcos or Tycho stars with linear space motion. Both input and output data is in the ICRS.
hipparcos | An entry from the Hipparcos catalog, at epoch J1991.25, with 'ra' in degrees(!) as per Hipparcos catalog units. | |
[out] | hip_2000 | The transformed input entry, at epoch J2000.0, with 'ra' in hours(!) as per the NOVAS convention. It may be the same as the input. |
References cat_entry::catalog, NOVAS_JD_HIP, cat_entry::ra, and transform_cat().