AstronomicalAlgorithms/RhoCosPhiTwo.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
/**************************************************************************
  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       conversion
@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 Phi
@param double Height

@return
converted value

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

#include <math.h>
#include <stddef.h>

#include "AstronomicalAlgorithms.h"

/* see page 82 */

double DoRhoCosPhiTwo(double doPhi, double doHeight)
   {
   double doReturnValue = (double) 0;

   double doU =  atan((((double) 1 - DO_EARTH_FLATTENING)
                       * tan(DEGREES2RADIAN(doPhi))));

   #if 0
   (void) fprintf(stdout, "\n%s line %ld doU = %8.6f", __FILE__, __LINE__,
                  RADIAN2DEGREES(doU));
   #endif

   doReturnValue = (   cos(doU)
                    +  (  (  doHeight
                           / (double) DO_EARTH_EQUATORIAL_RADIUS
                          )
                        * cos(doU)
                       )
                    );

   return doReturnValue;
   }