← Back to branch summary

~stellarium/stellarium/trunk

« back to all changes in this revision

Viewing changes to src/core/planetsephems/stellplanet.c

  • Committer: georg-zotti
  • Date: 2016-01-16 17:11:08 UTC
  • mfrom: (7721.1.41 socis2015-de430)
  • Revision ID: georg.zotti@univie.ac.at-20160116171108-rng3giwntia88mnb
Merged in SoCiS2015 project to get planetary positions from JPL DE430 and DE431 ephemerides

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/core/planetsephems/stellplanet.c
1
 
/*
2
 
Copyright (C) 2003 Fabien Chereau
3
 
 
4
 
This program is free software; you can redistribute it and/or modify
5
 
it under the terms of the GNU Library General Public License as published by
6
 
the Free Software Foundation; either version 2 of the License, or
7
 
(at your option) any later version.
8
 
 
9
 
This program is distributed in the hope that it will be useful,
10
 
but WITHOUT ANY WARRANTY; without even the implied warranty of
11
 
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12
 
GNU General Public License for more details.
13
 
 
14
 
You should have received a copy of the GNU General Public License
15
 
along with this program; if not, write to the Free Software
16
 
Foundation, Inc., 51 Franklin Street, Suite 500, Boston, MA  02110-1335, USA.
17
 
*/
18
 
 
19
 
#include "vsop87.h"
20
 
#include "elp82b.h"
21
 
#include "marssat.h"
22
 
#include "l1.h"
23
 
#include "tass17.h"
24
 
#include "gust86.h"
25
 
 
26
 
/* Chapter 31 Pg 206-207 Equ 31.1 31.2 , 31.3 using VSOP 87
27
 
 * Calculate planets rectangular heliocentric ecliptical coordinates
28
 
 * for given julian day. Values are in AU.
29
 
 * params : Julian day, rect coords */
30
 
void get_pluto_helio_coords(double jd, double * X, double * Y, double * Z);
31
 
 
32
 
void get_pluto_helio_coordsv(double jd,double xyz[3], void* unused)
33
 
  {get_pluto_helio_coords(jd, &xyz[0], &xyz[1], &xyz[2]);}
34
 
 
35
 
 
36
 
/* Return 0 of course... */
37
 
void get_sun_helio_coordsv(double jd,double xyz[3], void* unused)
38
 
  {xyz[0]=0.; xyz[1]=0.; xyz[2]=0.;}
39
 
 
40
 
void get_mercury_helio_coordsv(double jd,double xyz[3], void* unused)
41
 
  {GetVsop87Coor(jd,VSOP87_MERCURY,xyz);}
42
 
void get_venus_helio_coordsv(double jd,double xyz[3], void* unused)
43
 
  {GetVsop87Coor(jd,VSOP87_VENUS,xyz);}
44
 
 
45
 
void get_earth_helio_coordsv(const double jd,double xyz[3]) {
46
 
  double moon[3];
47
 
  GetVsop87Coor(jd,VSOP87_EMB,xyz);
48
 
  GetElp82bCoor(jd,moon);
49
 
    /* Earth != EMB:
50
 
       0.0121505677733761 = mu_m/(1+mu_m),
51
 
       mu_m = mass(moon)/mass(earth) = 0.01230002 */
52
 
  xyz[0] -= 0.0121505677733761 * moon[0];
53
 
  xyz[1] -= 0.0121505677733761 * moon[1];
54
 
  xyz[2] -= 0.0121505677733761 * moon[2];
55
 
}
56
 
 
57
 
void get_mars_helio_coordsv(double jd,double xyz[3], void* unused)
58
 
  {GetVsop87Coor(jd,VSOP87_MARS,xyz);}
59
 
void get_jupiter_helio_coordsv(double jd,double xyz[3], void* unused)
60
 
  {GetVsop87Coor(jd,VSOP87_JUPITER,xyz);}
61
 
void get_saturn_helio_coordsv(double jd,double xyz[3], void* unused)
62
 
  {GetVsop87Coor(jd,VSOP87_SATURN,xyz);}
63
 
void get_uranus_helio_coordsv(double jd,double xyz[3], void* unused)
64
 
  {GetVsop87Coor(jd,VSOP87_URANUS,xyz);}
65
 
void get_neptune_helio_coordsv(double jd,double xyz[3], void* unused)
66
 
  {GetVsop87Coor(jd,VSOP87_NEPTUNE,xyz);}
67
 
 
68
 
void get_mercury_helio_osculating_coords(double jd0,double jd,double xyz[3])
69
 
  {GetVsop87OsculatingCoor(jd0,jd,VSOP87_MERCURY,xyz);}
70
 
void get_venus_helio_osculating_coords(double jd0,double jd,double xyz[3])
71
 
  {GetVsop87OsculatingCoor(jd0,jd,VSOP87_VENUS,xyz);}
72
 
void get_earth_helio_osculating_coords(double jd0,double jd,double xyz[3])
73
 
  {GetVsop87OsculatingCoor(jd0,jd,VSOP87_EMB,xyz);}
74
 
void get_mars_helio_osculating_coords(double jd0,double jd,double xyz[3])
75
 
  {GetVsop87OsculatingCoor(jd0,jd,VSOP87_MARS,xyz);}
76
 
void get_jupiter_helio_osculating_coords(double jd0,double jd,double xyz[3])
77
 
  {GetVsop87OsculatingCoor(jd0,jd,VSOP87_JUPITER,xyz);}
78
 
void get_saturn_helio_osculating_coords(double jd0,double jd,double xyz[3])
79
 
  {GetVsop87OsculatingCoor(jd0,jd,VSOP87_SATURN,xyz);}
80
 
void get_uranus_helio_osculating_coords(double jd0,double jd,double xyz[3])
81
 
  {GetVsop87OsculatingCoor(jd0,jd,VSOP87_URANUS,xyz);}
82
 
void get_neptune_helio_osculating_coords(double jd0,double jd,double xyz[3])
83
 
  {GetVsop87OsculatingCoor(jd0,jd,VSOP87_NEPTUNE,xyz);}
84
 
 
85
 
/* Calculate the rectangular geocentric lunar coordinates to the inertial mean
86
 
 * ecliptic and equinox of J2000.
87
 
 * The geocentric coordinates returned are in units of AU.
88
 
 * This function is based upon the Lunar Solution ELP2000-82B by
89
 
 * Michelle Chapront-Touze and Jean Chapront of the Bureau des Longitudes,
90
 
 * Paris. ELP 2000-82B theory
91
 
 * param jd Julian day, rect pos */
92
 
void get_lunar_parent_coordsv(double jd,double xyz[3], void* unused)
93
 
  {GetElp82bCoor(jd,xyz);}
94
 
 
95
 
void get_phobos_parent_coordsv(double jd,double xyz[3], void* unused)
96
 
  {GetMarsSatCoor(jd,MARS_SAT_PHOBOS,xyz);}
97
 
void get_deimos_parent_coordsv(double jd,double xyz[3], void* unused)
98
 
  {GetMarsSatCoor(jd,MARS_SAT_DEIMOS,xyz);}
99
 
 
100
 
void get_io_parent_coordsv(double jd,double xyz[3], void* unused)
101
 
  {GetL1Coor(jd,L1_IO,xyz);}
102
 
void get_europa_parent_coordsv(double jd,double xyz[3], void* unused)
103
 
  {GetL1Coor(jd,L1_EUROPA,xyz);}
104
 
void get_ganymede_parent_coordsv(double jd,double xyz[3], void* unused)
105
 
  {GetL1Coor(jd,L1_GANYMEDE,xyz);}
106
 
void get_callisto_parent_coordsv(double jd,double xyz[3], void* unused)
107
 
  {GetL1Coor(jd,L1_CALLISTO,xyz);}
108
 
 
109
 
void get_mimas_parent_coordsv(double jd,double xyz[3], void* unused)
110
 
  {GetTass17Coor(jd,TASS17_MIMAS,xyz);}
111
 
void get_enceladus_parent_coordsv(double jd,double xyz[3], void* unused)
112
 
  {GetTass17Coor(jd,TASS17_ENCELADUS,xyz);}
113
 
void get_tethys_parent_coordsv(double jd,double xyz[3], void* unused)
114
 
  {GetTass17Coor(jd,TASS17_TETHYS,xyz);}
115
 
void get_dione_parent_coordsv(double jd,double xyz[3], void* unused)
116
 
  {GetTass17Coor(jd,TASS17_DIONE,xyz);}
117
 
void get_rhea_parent_coordsv(double jd,double xyz[3], void* unused)
118
 
  {GetTass17Coor(jd,TASS17_RHEA,xyz);}
119
 
void get_titan_parent_coordsv(double jd,double xyz[3], void* unused)
120
 
  {GetTass17Coor(jd,TASS17_TITAN,xyz);}
121
 
void get_hyperion_parent_coordsv(double jd,double xyz[3], void* unused)
122
 
  {GetTass17Coor(jd,TASS17_HYPERION,xyz);}
123
 
void get_iapetus_parent_coordsv(double jd,double xyz[3], void* unused)
124
 
  {GetTass17Coor(jd,TASS17_IAPETUS,xyz);}
125
 
 
126
 
void get_miranda_parent_coordsv(double jd,double xyz[3], void* unused)
127
 
  {GetGust86Coor(jd,GUST86_MIRANDA,xyz);}
128
 
void get_ariel_parent_coordsv(double jd,double xyz[3], void* unused)
129
 
  {GetGust86Coor(jd,GUST86_ARIEL,xyz);}
130
 
void get_umbriel_parent_coordsv(double jd,double xyz[3], void* unused)
131
 
  {GetGust86Coor(jd,GUST86_UMBRIEL,xyz);}
132
 
void get_titania_parent_coordsv(double jd,double xyz[3], void* unused)
133
 
  {GetGust86Coor(jd,GUST86_TITANIA,xyz);}
134
 
void get_oberon_parent_coordsv(double jd,double xyz[3], void* unused)
135
 
  {GetGust86Coor(jd,GUST86_OBERON,xyz);}
136