pmsafe

erfa.pmsafe(ra1, dec1, pmr1, pmd1, px1, rv1, ep1a, ep1b, ep2a, ep2b)[source]

Star proper motion: update star catalog data for space motion, with special handling to handle the zero parallax case.

Parameters:
ra1double array
dec1double array
pmr1double array
pmd1double array
px1double array
rv1double array
ep1adouble array
ep1bdouble array
ep2adouble array
ep2bdouble array
Returns:
ra2double array
dec2double array
pmr2double array
pmd2double array
px2double array
rv2double array

Notes

Wraps ERFA function eraPmsafe. The ERFA documentation is:

- - - - - - - - - -
 e r a P m s a f e
- - - - - - - - - -

Star proper motion:  update star catalog data for space motion, with
special handling to handle the zero parallax case.

Given:
   ra1    double     right ascension (radians), before
   dec1   double     declination (radians), before
   pmr1   double     RA proper motion (radians/year), before
   pmd1   double     Dec proper motion (radians/year), before
   px1    double     parallax (arcseconds), before
   rv1    double     radial velocity (km/s, +ve = receding), before
   ep1a   double     "before" epoch, part A (Note 1)
   ep1b   double     "before" epoch, part B (Note 1)
   ep2a   double     "after" epoch, part A (Note 1)
   ep2b   double     "after" epoch, part B (Note 1)

Returned:
   ra2    double     right ascension (radians), after
   dec2   double     declination (radians), after
   pmr2   double     RA proper motion (radians/year), after
   pmd2   double     Dec proper motion (radians/year), after
   px2    double     parallax (arcseconds), after
   rv2    double     radial velocity (km/s, +ve = receding), after

Returned (function value):
          int        status:
                        -1 = system error (should not occur)
                         0 = no warnings or errors
                         1 = distance overridden (Note 6)
                         2 = excessive velocity (Note 7)
                         4 = solution didn't converge (Note 8)
                      else = binary logical OR of the above warnings

Notes:

1) The starting and ending TDB dates ep1a+ep1b and ep2a+ep2b are
   Julian Dates, apportioned in any convenient way between the two
   parts (A and B).  For example, JD(TDB)=2450123.7 could be
   expressed in any of these ways, among others:

          epNa            epNb

       2450123.7           0.0       (JD method)
       2451545.0       -1421.3       (J2000 method)
       2400000.5       50123.2       (MJD method)
       2450123.5           0.2       (date & time method)

   The JD method is the most natural and convenient to use in cases
   where the loss of several decimal digits of resolution is
   acceptable.  The J2000 method is best matched to the way the
   argument is handled internally and will deliver the optimum
   resolution.  The MJD method and the date & time methods are both
   good compromises between resolution and convenience.

2) In accordance with normal star-catalog conventions, the object's
   right ascension and declination are freed from the effects of
   secular aberration.  The frame, which is aligned to the catalog
   equator and equinox, is Lorentzian and centered on the SSB.

   The proper motions are the rate of change of the right ascension
   and declination at the catalog epoch and are in radians per TDB
   Julian year.

   The parallax and radial velocity are in the same frame.

3) Care is needed with units.  The star coordinates are in radians
   and the proper motions in radians per Julian year, but the
   parallax is in arcseconds.

4) The RA proper motion is in terms of coordinate angle, not true
   angle.  If the catalog uses arcseconds for both RA and Dec proper
   motions, the RA proper motion will need to be divided by cos(Dec)
   before use.

5) Straight-line motion at constant speed, in the inertial frame, is
   assumed.

6) An extremely small (or zero or negative) parallax is overridden
   to ensure that the object is at a finite but very large distance,
   but not so large that the proper motion is equivalent to a large
   but safe speed (about 0.1c using the chosen constant).  A warning
   status of 1 is added to the status if this action has been taken.

7) If the space velocity is a significant fraction of c (see the
   constant VMAX in the function eraStarpv), it is arbitrarily set
   to zero.  When this action occurs, 2 is added to the status.

8) The relativistic adjustment carried out in the eraStarpv function
   involves an iterative calculation.  If the process fails to
   converge within a set number of iterations, 4 is added to the
   status.

Called:
   eraSeps      angle between two points
   eraStarpm    update star catalog data for space motion

This revision:   2023 April 7

Copyright (C) 2013-2023, NumFOCUS Foundation.
Derived, with permission, from the SOFA library.  See notes at end of file.