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
|
/*
This program is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 2 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program; if not, write to the Free Software
Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
Copyright 2004 Liam Girdwood
*/
#include <libnova/transform.h>
#include <libnova/julian_day.h>
#include <libnova/utility.h>
#include <stdio.h>
/*
* Do some RA/DEC <--> ALT/AZ conversions for Alnilam
* in Leon (Spain) on 25/4/2004
*/
int main (int argc, const char *argv[])
{
struct lnh_equ_posn hobject, hequ;
struct lnh_lnlat_posn hobserver;
struct ln_equ_posn object, equ;
struct ln_hrz_posn hrz;
struct lnh_hrz_posn hhrz;
struct ln_lnlat_posn observer;
struct ln_date date;
double JD;
/*
* observers position
* longitude is measured positively eastwards
* i.e. Long 5d36m30W (Leon, Spain) = 354d24m30
* Lat for Leon = Lat 42d35m40 N
*/
hobserver.lng.degrees = -5;
hobserver.lng.minutes = 36;
hobserver.lng.seconds = 30.0;
hobserver.lat.degrees = 42;
hobserver.lat.minutes = 35;
hobserver.lat.seconds = 40.0;
/* Alnilam */
hobject.ra.hours = 5;
hobject.ra.minutes = 36;
hobject.ra.seconds = 27.0;
hobject.dec.neg = 1;
hobject.dec.degrees = 1;
hobject.dec.minutes = 12;
hobject.dec.seconds = 0.0;
/* UT date and time */
date.years = 2004;
date.months = 4;
date.days = 25;
date.hours = 12;
date.minutes = 18;
date.seconds = 49.0;
JD = ln_get_julian_day(&date);
ln_hequ_to_equ(&hobject, &object);
ln_hlnlat_to_lnlat(&hobserver, &observer);
ln_get_hrz_from_equ(&object, &observer, JD, &hrz);
fprintf(stdout, "(Alnilam) Equ to Horiz ALT %f\n", hrz.alt);
fprintf(stdout, "(Alnilam) Equ to Horiz AZ %f\n", hrz.az);
ln_hrz_to_hhrz(&hrz, &hhrz);
fprintf(stdout, "ALT %d:%d:%f AZ %d:%d:%f\n",
hhrz.alt.degrees, hhrz.alt.minutes, hhrz.alt.seconds,
hhrz.az.degrees, hhrz.az.minutes, hhrz.az.seconds);
ln_get_equ_from_hrz (&hrz, &observer, JD, &equ);
fprintf(stdout, "(Alnilam) Horiz to Equ RA %f\n", equ.ra);
fprintf(stdout, "(Alnilam) Horiz to Equ DEC %f\n", equ.dec);
ln_equ_to_hequ(&equ, &hequ);
fprintf(stdout, "RA %d:%d:%f DEC %d:%d:%f\n",
hequ.ra.hours, hequ.ra.minutes, hequ.ra.seconds,
hequ.dec.degrees, hequ.dec.minutes, hequ.dec.seconds);
return 0;
}
|
