Procedures to set constants and read data files. More...

Functions/Subroutines
subroutine	set_thesky_constants ()
	Set the initial values of the variables and constants used in this package.

subroutine	thesky_readdata ()
	Read input data files.

subroutine	read_deltat ()
	Read deltat.dat with historical values for DeltaT.

subroutine	readconstellations ()
	Read constellation names: abbreviation, Latin, genitive, Dutch and English.

subroutine	readnutation ()
	Read nutation input files.

subroutine	readplanetdata ()
	Reads VSOP planet data from planets.dat.

subroutine	readpluto ()
	Read periodic terms for the position of Pluto.

subroutine	readplanetelements ()
	Read the data needed to compute the orbital elements of the planets.

subroutine	readmoondata ()
	Read ELP-82B periodic terms from moondata.dat.

subroutine	readmoondata_la ()
	Read low-accuracy moon-position data from Meeus (moon_la.dat)

subroutine	elp_mpp02_initialise_and_read_files (mode, ierr)
	Initialise ELP-MPP02 data and read the data files.

subroutine	elp_mpp02_initialise (mode)
	Initialization of the constants and parameters used for the evaluation of the ELP/MPP02 series.

subroutine	elp_mpp02_read_files (ierr)
	Read the six data files containing the ELP/MPP02 series.

real(double) function	elp_dms2rad (deg, min, sec)
	Function for the conversion: sexagesimal degrees -> radians.

subroutine	readasteroidelements ()
	Read orbital-element data for the asteroids.

subroutine	readcometelements ()
	Read orbital-element data for the comets.

subroutine	read_bsc
	Read the Bright Star Catalogue.

Detailed Description

Procedures to set constants and read data files.

Function/Subroutine Documentation

◆ elp_dms2rad()

real(double) function thesky_data::elp_dms2rad	(	integer, intent(in)	deg,
		integer, intent(in)	min,
		real(double), intent(in)	sec )

Function for the conversion: sexagesimal degrees -> radians.

Parameters

deg	Degrees
min	Minutes or arc
sec	Seconds of arc

Return values

elp_dms2rad The angle in radians

Definition at line 1129 of file data.f90.

References elp_dms2rad().

Referenced by elp_dms2rad(), and elp_mpp02_initialise().

Here is the call graph for this function:

◆ elp_mpp02_initialise()

subroutine thesky_data::elp_mpp02_initialise ( integer, intent(in) mode )

Initialization of the constants and parameters used for the evaluation of the ELP/MPP02 series.

Parameters

mode	Index of the corrections to the constants: 0: LLR observations for 1970-2000, 1: DE405 ephemeris for 1950-2060

Note

Returns: Set of the constants of ELPMPP02 solution in module TheSky_elp_mpp02_constants.
Remarks: The nominal values of some constants have to be corrected. There are two sets of corrections, one of which can be chosen using the parameter 'mode' (used in elp_mpp02_initialise()):
- mode=0, the constants are fitted to LLR observations provided from 1970 to 2001; it is the default value;
- mode=1, the constants are fitted to DE405 ephemeris over one century (1950-2060); the lunar angles W1, W2, W3 receive also additive corrections to the secular coefficients.
Moon constants: nu : mean motion of the Moon (W1(1,1)) (Nu) g : half coefficient of sin(F) in latitude (Gamma) e : half coefficient of sin(l) in longitude (E) np : mean motion of EMB (eart(1)) (n') ep : eccentricity of EMB (e')

p is the precession rate and t is the time

See also

ELPdoc: Lunar solution ELP, version ELP/MPP02, Jean Chapront and Gerard Francou, October 2002

Definition at line 786 of file data.f90.

References elp_dms2rad().

Referenced by elp_mpp02_initialise_and_read_files().

Here is the call graph for this function:

◆ elp_mpp02_initialise_and_read_files()

subroutine thesky_data::elp_mpp02_initialise_and_read_files	(	integer, intent(in)	mode,
		integer, intent(out)	ierr )

Initialise ELP-MPP02 data and read the data files.

Parameters

mode	Index of the corrections to the constants: 0: LLR observations for 1970-2000, 1: DE405 ephemeris for 1950-2060 (input)
ierr	File error index: ierr=0: no error, ierr=1: file error (output)

Definition at line 735 of file data.f90.

References elp_mpp02_initialise(), and elp_mpp02_read_files().

Referenced by thesky_moon::elp_mpp02_xyz().

Here is the call graph for this function:

◆ elp_mpp02_read_files()

subroutine thesky_data::elp_mpp02_read_files ( integer, intent(out) ierr )

Read the six data files containing the ELP/MPP02 series.

Parameters

ierr	File error index: ierr=0: no error, ierr=1: file error (output)

Note

module elp_mpp02_constants: Set of the constants of ELP/MPP02 solution (input)
module elp_mpp02_series: Series of the ELP/MPP02 solution (output)

See also

Data files: ftp://cyrano-se.obspm.fr/pub/2_lunar_solutions/2_elpmpp02/

Definition at line 999 of file data.f90.

References thesky_constants::theskydatadir.

Referenced by elp_mpp02_initialise_and_read_files().

◆ read_bsc()

subroutine thesky_data::read_bsc

Read the Bright Star Catalogue.

Definition at line 1299 of file data.f90.

References thesky_bscdata::bsc_abbr, thesky_bscdata::bsc_bv, thesky_bscdata::bsc_dec, thesky_bscdata::bsc_mult, thesky_bscdata::bsc_name, thesky_bscdata::bsc_par, thesky_bscdata::bsc_pma, thesky_bscdata::bsc_pmd, thesky_bscdata::bsc_ra, thesky_bscdata::bsc_ri, thesky_bscdata::bsc_rv, thesky_bscdata::bsc_sao, thesky_bscdata::bsc_sptype, thesky_bscdata::bsc_ub, thesky_bscdata::bsc_var, thesky_bscdata::bsc_vm, thesky_bscdata::bsc_vm_indx, thesky_bscdata::n_bsc, and thesky_constants::theskydatadir.

◆ read_deltat()

subroutine thesky_data::read_deltat

Read deltat.dat with historical values for DeltaT.

Note

See also: http://hemel.waarnemen.com/Computing/deltat.html

Definition at line 128 of file data.f90.

References thesky_constants::deltat_0, thesky_constants::deltat_accel, thesky_constants::deltat_change, thesky_constants::deltat_maxyr, thesky_constants::deltat_minyr, thesky_constants::deltat_n, thesky_constants::deltat_nmax, thesky_constants::deltat_values, thesky_constants::deltat_years, and thesky_constants::theskydatadir.

Referenced by thesky_readdata().

◆ readasteroidelements()

subroutine thesky_data::readasteroidelements

Read orbital-element data for the asteroids.

Note

data are passed via the module planetdata, in asterElems(*,1:9):
Epoch (JD), a, e, i, omega, Omega, M, H, G, for J2000.0

See also

asteroids.dat: http://sf.net/projects/libthesky/files/asteroids.dat.bz2 (selection: H<15, a<100 - 25% of all bodies)
original: http://ssd.jpl.nasa.gov/dat/ELEMENTS.NUMBR.gz

Definition at line 1155 of file data.f90.

References thesky_planetdata::asterelems, thesky_planetdata::asternames, thesky_planetdata::nasteroids, and thesky_constants::theskydatadir.

Referenced by thesky_readdata().

◆ readcometelements()

subroutine thesky_data::readcometelements

Read orbital-element data for the comets.

Note

data are passed via the module comet_data, in cometElems(*,1:9)
colums: 1: Epoch (JD), 2: q, 3: e, 4: i, 5: omega, 6: Omega, 7: Tp, for J2000.0, 8: H, 9: G

See also

comets.dat: http://ssd.jpl.nasa.gov/dat/ELEMENTS.COMET
comets_mpc.dat: http://www.minorplanetcenter.net/iau/Ephemerides/Comets/Soft00Cmt.txt

Definition at line 1210 of file data.f90.

References thesky_cometdata::cometdatfile, thesky_cometdata::cometdiedatp, thesky_cometdata::cometelems, thesky_cometdata::cometnames, thesky_cometdata::ncomets, thesky_cometdata::ncometsmax, and thesky_constants::theskydatadir.

Referenced by thesky_readdata().

◆ readconstellations()

subroutine thesky_data::readconstellations

Read constellation names: abbreviation, Latin, genitive, Dutch and English.

Definition at line 187 of file data.f90.

References thesky_stardata::conabr, thesky_stardata::conid, thesky_stardata::conidabr, thesky_stardata::coniddecl, thesky_stardata::conidral, thesky_stardata::conidrau, thesky_stardata::enconnames, thesky_stardata::genconnames, thesky_stardata::latconnames, thesky_stardata::nconid, thesky_stardata::nconstel, thesky_stardata::nlconnames, and thesky_constants::theskydatadir.

Referenced by thesky_readdata().

◆ readmoondata()

subroutine thesky_data::readmoondata

Read ELP-82B periodic terms from moondata.dat.

Note: Constants are defined in the module moondata

See also

Chapront-Touzé & Chapront, A&A, 124, 50 (1983)
Chapront-Touzé & Chapront, A&A, 190, 342 (1988)
ftp://cdsarc.u-strasbg.fr/pub/cats/VI/79/

Definition at line 415 of file data.f90.

◆ readmoondata_la()

subroutine thesky_data::readmoondata_la

Read low-accuracy moon-position data from Meeus (moon_la.dat)

Definition at line 693 of file data.f90.

References thesky_planetdata::moonla_arg, thesky_planetdata::moonla_lrb, and thesky_constants::theskydatadir.

Referenced by thesky_readdata().

◆ readnutation()

subroutine thesky_data::readnutation

Read nutation input files.

Definition at line 242 of file data.f90.

References thesky_constants::nutationdat, and thesky_constants::theskydatadir.

Referenced by thesky_readdata().

◆ readplanetdata()

subroutine thesky_data::readplanetdata

Reads VSOP planet data from planets.dat.

Definition at line 276 of file data.f90.

References thesky_constants::theskydatadir, thesky_planetdata::vsopdat, thesky_planetdata::vsopnblk, thesky_planetdata::vsopnls, and thesky_planetdata::vsoptruncs.

Referenced by thesky_readdata().

◆ readplanetelements()

subroutine thesky_data::readplanetelements

Read the data needed to compute the orbital elements of the planets.

Note

use planetelements() to compute the orbital elements
L, a, e, i, Omega, pi for Mer-Nep, mean equinox of date and J2000.0 (=2x6x8)
Terms in one row are a_0, a_1, a_2, a_3 and should be used as: Sum(i=0,3) a_i t^i, with t the time in centuries since 2000.0
Angles are still expressed in degrees

Definition at line 366 of file data.f90.

References thesky_planetdata::plelemdata, and thesky_constants::theskydatadir.

Referenced by thesky_readdata().