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
|
% /*-----------------------------------------------------------------*/
% /*!
% \file mexmultiple.m
% \brief Example of usage of the multiple file access functions
% with the Octave/Matlab interface.
% This example reads the constant EMRAT, AU, GM_Mer and print their values.
% It computes for a date
% the geocentric moon coordinates,
% the value TT-TDB
% the heliocentric coordinates of Mars.
%
% \author M. Gastineau
% Astronomie et Systemes Dynamiques, IMCCE, CNRS, Observatoire de Paris.
%
% Copyright, 2018-2023, CNRS
% email of the author : Mickael.Gastineau@obspm.fr
%
% History:
% */
% /*-----------------------------------------------------------------*/
%
%
% /*-----------------------------------------------------------------*/
% /* License of this file :
% This file is 'triple-licensed', you have to choose one of the three licenses
% below to apply on this file.
%
% CeCILL-C
% The CeCILL-C license is close to the GNU LGPL.
% ( http://www.cecill.info/licences/Licence_CeCILL-C_V1-en.html )
%
% or CeCILL-B
% The CeCILL-B license is close to the BSD.
% (http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.txt)
%
% or CeCILL v2.1
% The CeCILL license is compatible with the GNU GPL.
% ( http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.html )
%
%
% This library is governed by the CeCILL-C, CeCILL-B or the CeCILL license under
% French law and abiding by the rules of distribution of free software.
% You can use, modify and/ or redistribute the software under the terms
% of the CeCILL-C,CeCILL-B or CeCILL license as circulated by CEA, CNRS and INRIA
% at the following URL 'http://www.cecill.info'.
%
% As a counterpart to the access to the source code and rights to copy,
% modify and redistribute granted by the license, users are provided only
% with a limited warranty and the software's author, the holder of the
% economic rights, and the successive licensors have only limited
% liability.
%
% In this respect, the user's attention is drawn to the risks associated
% with loading, using, modifying and/or developing or reproducing the
% software by the user in light of its specific status of free software,
% that may mean that it is complicated to manipulate, and that also
% therefore means that it is reserved for developers and experienced
% professionals having in-depth computer knowledge. Users are therefore
% encouraged to load and test the software's suitability as regards their
% requirements in conditions enabling the security of their systems and/or
% data to be ensured and, more generally, to use and operate it in the
% same conditions as regards security.
%
% The fact that you are presently reading this means that you have had
% knowledge of the CeCILL-C,CeCILL-B or CeCILL license and that you accept its terms.
% */
% /*-----------------------------------------------------------------*/
% /*-----------------------------------------------------------------*/
% /* main program */
% /*-----------------------------------------------------------------*/
jd0=2442457
dt=0.5E0
% open the ephemeris file
peph = CalcephBin.open('example1.dat');
disp('The ephemeris is already opened\n')
% print the values of AU, EMRAT and GM_Mer
AU = peph.getconstant('AU')
EMRAT = peph.getconstant('EMRAT')
GM_Mer = peph.getconstant('GM_Mer')
% compute and print the coordinates
% the geocentric moon coordinates
disp('geocentric coordinates of the Moon in AU and AU/day')
[ found, moon ] = peph.getidbyname("Moon",0)
[ found, earth ] = peph.getidbyname("Earth",0)
PV = peph.compute(jd0, dt, moon, earth)
disp('geocentric coordinates of the Moon in km and km/s')
PV = peph.compute_unit(jd0, dt, NaifId.MOON, NaifId.EARTH,
Constants.UNIT_KM+Constants.UNIT_SEC+Constants.USE_NAIFID)
% the value TT-TDB
PV = peph.compute(jd0, dt, 16, 0);
printf('TT-TDB = %f\n',PV(1));
% the heliocentric coordinates of Mars
disp('heliocentric coordinates of Mars')
PV = peph.compute_unit(jd0, dt, NaifId.MARS_BARYCENTER, NaifId.SUN,
Constants.UNIT_AU+Constants.UNIT_DAY+Constants.USE_NAIFID)
% print the whole list of the constants
disp('list of constants');
for j=1:peph.getconstantcount()
[ nameconstant, valueconstant ] = peph.getconstantindex(j);
printf('%s\t= %e\n',nameconstant, valueconstant)
end
% close the ephemeris file
peph.close();
disp('The ephemeris is already closed');
|
