1 | #include "slalib.h"
|
---|
2 | #include "slamac.h"
|
---|
3 | void slaAfin ( char *string, int *iptr, float *a, int *j )
|
---|
4 | /*
|
---|
5 | ** - - - - - - - -
|
---|
6 | ** s l a A f i n
|
---|
7 | ** - - - - - - - -
|
---|
8 | **
|
---|
9 | ** Sexagesimal character string to angle.
|
---|
10 | **
|
---|
11 | ** (single precision)
|
---|
12 | **
|
---|
13 | ** Given:
|
---|
14 | ** string c*(*) string containing deg, arcmin, arcsec fields
|
---|
15 | ** iptr int where to start decode (1st = 1)
|
---|
16 | **
|
---|
17 | ** Returned:
|
---|
18 | ** iptr int advanced past the decoded angle
|
---|
19 | ** a float angle in radians
|
---|
20 | ** j int status: 0 = OK
|
---|
21 | ** +1 = default, A unchanged
|
---|
22 | ** -1 = bad degrees )
|
---|
23 | ** -2 = bad arcminutes ) (note 3)
|
---|
24 | ** -3 = bad arcseconds )
|
---|
25 | **
|
---|
26 | ** Example:
|
---|
27 | **
|
---|
28 | ** argument before after
|
---|
29 | **
|
---|
30 | ** string '-57 17 44.806 12 34 56.7' unchanged
|
---|
31 | ** iptr 1 16 (points to 12...)
|
---|
32 | **
|
---|
33 | ** a ? -1.00000f
|
---|
34 | ** j ? 0
|
---|
35 | **
|
---|
36 | ** A further call to slaAfin, without adjustment of iptr, will
|
---|
37 | ** decode the second angle, 12deg 34min 56.7sec.
|
---|
38 | **
|
---|
39 | ** Notes:
|
---|
40 | **
|
---|
41 | ** 1) The first three "fields" in string are degrees, arcminutes,
|
---|
42 | ** arcseconds, separated by spaces or commas. The degrees field
|
---|
43 | ** may be signed, but not the others. The decoding is carried
|
---|
44 | ** out by the dfltin routine and is free-format.
|
---|
45 | **
|
---|
46 | ** 2) Successive fields may be absent, defaulting to zero. For
|
---|
47 | ** zero status, the only combinations allowed are degrees alone,
|
---|
48 | ** degrees and arcminutes, and all three fields present. If all
|
---|
49 | ** three fields are omitted, a status of +1 is returned and a is
|
---|
50 | ** unchanged. In all other cases a is changed.
|
---|
51 | **
|
---|
52 | ** 3) Range checking:
|
---|
53 | ** The degrees field is not range checked. However, it is
|
---|
54 | ** expected to be integral unless the other two fields are absent.
|
---|
55 | ** The arcminutes field is expected to be 0-59, and integral if
|
---|
56 | ** the arcseconds field is present. If the arcseconds field
|
---|
57 | ** is absent, the arcminutes is expected to be 0-59.9999...
|
---|
58 | ** The arcseconds field is expected to be 0-59.9999...
|
---|
59 | **
|
---|
60 | ** 4) Decoding continues even when a check has failed. Under these
|
---|
61 | ** circumstances the field takes the supplied value, defaulting
|
---|
62 | ** to zero, and the result a is computed and returned.
|
---|
63 | **
|
---|
64 | ** 5) Further fields after the three expected ones are not treated
|
---|
65 | ** as an error. The pointer iptr is left in the correct state
|
---|
66 | ** for further decoding with the present routine or with slaDfltin
|
---|
67 | ** etc. See the example, above.
|
---|
68 | **
|
---|
69 | ** 6) If string contains hours, minutes, seconds instead of degrees
|
---|
70 | ** etc, or if the required units are turns (or days) instead of
|
---|
71 | ** radians, the result a should be multiplied as follows:
|
---|
72 | ** for to obtain multiply
|
---|
73 | ** string a in a by
|
---|
74 | ** d ' " radians 1 = 1.0f
|
---|
75 | ** d ' " turns 1/2pi = 0.1591549430918953358f
|
---|
76 | ** h m s radians 15 = 15.0f
|
---|
77 | ** h m s days 15/2pi = 2.3873241463784300365f
|
---|
78 | **
|
---|
79 | ** Called: slaDfltin
|
---|
80 | **
|
---|
81 | ** Last revision: 16 November 1993
|
---|
82 | **
|
---|
83 | ** Copyright P.T.Wallace. All rights reserved.
|
---|
84 | */
|
---|
85 | {
|
---|
86 | double ad;
|
---|
87 |
|
---|
88 | /* Call the double precision version */
|
---|
89 | slaDafin ( string, iptr, &ad, j );
|
---|
90 | if ( *j <= 0 ) *a = (float) ad;
|
---|
91 | }
|
---|