1 /*
2 *+
3 *  Name:
4 *     palNutc
5 
6 *  Purpose:
7 *     Calculate nutation longitude & obliquoty components
8 
9 *  Language:
10 *     Starlink ANSI C
11 
12 *  Type of Module:
13 *     Library routine
14 
15 *  Invocation:
16 *     void palNutc( double date, double * dpsi, double *deps, double *eps0 );
17 
18 *  Arguments:
19 *     date = double (Given)
20 *        TT as modified Julian date (JD-2400000.5)
21 *     dpsi = double * (Returned)
22 *        Nutation in longitude
23 *     deps = double * (Returned)
24 *        Nutation in obliquity
25 *     eps0 = double * (Returned)
26 *        Mean obliquity.
27 
28 *  Description:
29 *     Calculates the longitude * obliquity components and mean obliquity
30 *     using the SOFA/ERFA library.
31 
32 *  Authors:
33 *     TIMJ: Tim Jenness (JAC, Hawaii)
34 *     {enter_new_authors_here}
35 
36 *  Notes:
37 *     - Calls eraObl06 and eraNut06a and therefore uses the IAU 206
38 *       precession/nutation model.
39 *     - Note the change from SLA/F regarding the date. TT is used
40 *       rather than TDB.
41 
42 *  History:
43 *     2012-03-05 (TIMJ):
44 *        Initial version
45 *        Adapted with permission from the Fortran SLALIB library.
46 *     {enter_further_changes_here}
47 
48 *  Copyright:
49 *     Copyright (C) 2012 Science and Technology Facilities Council.
50 *     All Rights Reserved.
51 
52 *  Licence:
53 *     This program is free software; you can redistribute it and/or
54 *     modify it under the terms of the GNU General Public License as
55 *     published by the Free Software Foundation; either version 3 of
56 *     the License, or (at your option) any later version.
57 *
58 *     This program is distributed in the hope that it will be
59 *     useful, but WITHOUT ANY WARRANTY; without even the implied
60 *     warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
61 *     PURPOSE. See the GNU General Public License for more details.
62 *
63 *     You should have received a copy of the GNU General Public License
64 *     along with this program; if not, write to the Free Software
65 *     Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
66 *     MA 02110-1301, USA.
67 
68 *  Bugs:
69 *     {note_any_bugs_here}
70 *-
71 */
72 
73 #include "pal.h"
74 #include "palmac.h"
75 #include "pal1sofa.h"
76 
palNutc(double date,double * dpsi,double * deps,double * eps0)77 void palNutc( double date, double * dpsi, double *deps, double *eps0 ) {
78   eraNut06a( PAL__MJD0, date, dpsi, deps );
79   *eps0 = eraObl06( PAL__MJD0, date );
80 }
81