AstronomicalAlgorithms/DateFromYearDay.c at main · efalk/AstronomicalAlgorithms · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
/**************************************************************************
  AstronomicalAlgorithms

  A portable ANSI C implementation of some of the algorithms published in

        Astronomical Algorithms
        by Jean Meeus
        2nd edition (December 1998)
        Willmann-Bell
        ISBN: 0943396638

  by Christophe DAVID (christophe.david@christophedavid.org)

  You may use parts of this source code as long as

       - you mention clearly that its latest version can be obtained
         free of charge at

         http://www.christophedavid.org/

       AND

      - you send me a free copy of whatever you make using this code.

     Comments and suggestions welcome.

**************************************************************************/

/*!
@file
@brief       calculates a date from the day in the year
@author      Christophe DAVID \n
             christophe.david@christophedavid.org \n
             http://www.christophedavid.org
@since       01/07/1999
@version     1.0
@date        05/04/2001
@bug         no known bug

@param double  YearDay
@param short   year
@param double *month
@param double *day

@return
- 0 : completed successfully
- 1 : error in parameter 1
- 2 : error in parameter 2
- 3 : error in parameter 3
- 4 : error in parameter 4

@if logger
@image html http://www.mot.be/cgi-bin/logger.cgi?DateFromYearDay.c
@endif

*/

/* see page 66 */

#ifdef ASTROALGO
   #include <math.h>
   #include <stdio.h>
   #include "AstronomicalAlgorithms.h"

   #if _WIN32
       __declspec(dllexport) short __stdcall
   #endif
#else
   short
#endif

int
ShDateFromYearDay(double doYD, short shYear, double *pdoM, double *pdoD)
   {
   short shReturnValue = (short) 0;
   short shIsLeap      = (short) 0;

   if (pdoM != NULL)
      {
      *pdoM = (double) -1;
      }

   if (pdoD != NULL)
      {
      *pdoD = (double) -1;
      }

   if (ShIsLeapYear(&shIsLeap, shYear,
                         (short) (shYear < 1582 ? JULIAN : GREGORIAN)) != 0)
      {
      shReturnValue = (short) 2;
      }

   else if (   (doYD < 0)
            || (doYD > (shIsLeap == 0 ? 365 : 366))
           )
      {
      shReturnValue = (short) 1;
      }

   else if (pdoM == NULL)
      {
      shReturnValue = (short) 3;
      }

   else if (pdoD == NULL)
      {
      shReturnValue = (short) 4;
      }

   else
      {
      if (doYD < 32)
         {
         *pdoM = (double) 1;
         }
      else
         {
         *pdoM = floor(  (   (  (double) 9
                              * (  (double) (shIsLeap == 1 ? 1 : 2)
                                 + doYD
                                )
                             )
                           / (double) 275
                         )
                       + (double) (0.98)
                       );

         }

      *pdoD = (   (  doYD
                   - floor( (  ((double) 275 * (*pdoM))
                             / (double) 9
                            )
                          )
                   )

                +  (  (double) (shIsLeap == 1 ? 1 : 2)
                    * floor(  (*pdoM + (double) 9)
                            / (double) 12
                           )
                   )
                + (double) 30
              );

      }

   return shReturnValue;
   }