montu 0.5.2

Montu Python: astronomical ephemerides for the ancient world

Montu Python /mnṯw ꜥꜣpp(y)/

Astronomical ephemerides for the ancient world

MontuPython (transileterated mnṯw ꜥꜣpp(y)) is a Python package intended to compute astronomical ephemerides in the ancient world, thousands of years before present. It was initially designed to compute ephemerides for the ancient Egypt, but it can also be used to study astronomical phenomena in other sites of interest for cultural astronomy (archeoastronomy).

Download and install

Describe here how the package can be downloaded and install it in different arquitectures.

If you are using PyPI installation it's as simple as:

pip install montu

You can also test the unstable version of the package with:

pip install -i https://test.pypi.org/simple/ montu

Quick start

In this section you should provide the most simple instructions to use your package.

For instance:

import montu
print(montu.version())

Code examples

For a fully-fledged working example see examples/montunctions.ipynb.

In order to properly use MontuPython you need to import the package and load the required data:

import montu
from montu import *
Montu.load_kernels()
allstars=Stars()

It may take a long time to download the planetary kernels, so be patient!

Get an historical date in all calendars and scales

One of the most interesting and basic functionalities of MontuPython is to convert date among different type of calendars and astronomical scales. You may taste these functionalities using:

mtime = MonTime('bce2501-01-01 12:00:00')

whose output will be:

Montu Time Object:
--------------------------
General:
    Calendar: proleptic
    Is bce: True
    Components UTC: [-1, 2500, 1, 1, 12, 0, 0, 0]
Uniform scales:
    Terrestrial time:
        tt: -142006202700.32
        jtd: 807954.6909685184
    UTC time:
        et: -142006262400.0
        jed: 807953.9999999999
    Delta-t = TT - UTC = 59699.68000000001
Strings:
    Date in SPICE format: 2501 B.C. 01-01 12:00:00.000000
    Date in proleptic calendar: -2500-01-01 12:00:00.0
    Date in proleptic calendar (jul.year): 2214-01-01 12:00:00.000000
    Date in mixed calendar: -2500-01-22 12:00:00
    Date in mixed calendar (jul.year): 2212-01-22 12:00:00
Objects:
    Date in datetime64 format: -2500-01-01T12:00:00.000
    Date in datetime format proleptic: 2500-01-01 12:00:00
    Date in datetime format proleptic (julian year): 2214-01-01 12:00:00
    Date in datetime format mixed: 2500-01-22 12:00:00
    Date in datetime format mixed (julian year): 2212-01-22 12:00:00
    Date in PyPlanet Epoch: 807953.9999999999
    Date in PyEphem Epoch: -2501/1/22 12:00:00
    Date in AstroPy Time: 807954.6909685184
Astronomical properties at Epoch:
    True obliquity of ecliptic: 23:58:33.587
    True nutation longitude: 00:00:10.214
    Greenwhich Meridian Sidereal Time: 18:40:25.323
Hash: 3649516878321859143

Notice that the date in Gregorian proleptic will be bce2501-01-01 but in mixed calendar will be bce 2501-01-22.

Position of mars at an historical date and place

First we need to prepare some basic objects, the Earth (where the observer is), the location on Earth, the time of observation and the object to be observed:

earth = PlanetaryBody('Earth')
tebas = ObservingSite(planet=earth,lon=33,lat=24,height=0)
mtime = MonTime('bce2501-01-01 12:00:00')
mars = PlanetaryBody('Mars')

Now we can ask for the position of Mars:

mars.calculate_sky_position(mtime,tebas)

The result will be:

Computing position of body 'mars' at epoch: jtd = 807954.6909685184 
Updating orientation of site (old time 2000-01-01 11:58:56.126200, new time 2501 B.C. 01-01 12:00:00.000000)
Method 'SPICE':
	Position Epoch:  -2500-01-01 12:00:00.0 807953.9999999999
	Coordinates @ J2000: 
		Equatorial: 12:31:48.754 01:37:12.184
		Ecliptic: 186:39:46.949 04:38:36.308
	Coordinates @ Epoch : 
		Equatorial: 08:32:9.796 24:06:28.555
		Ecliptic: 124:21:21.542 04:39:5.339
	Observing conditions: 
		Distance to site [au]:  0.660450348841685
		Distance to sun [au]:  1.6261149729985933
		Solar elongation [deg]:  157:49:18.876
		Phase angle [deg]:  13:21:31.981
		Magnitude:  -1.1
	Other properties: 
		Local true sidereal time:  20:52:25.323
		Hour angle @ Epoch:  12:20:15.527
		Local coordinates @ Epoch:  06:11:24.275 -41:38:31.094

There are different methods to calculate the position of a planets with MontuPython. You can try all of them:

mars.calculate_sky_position(mtime,tebas,method='all')

NOTE: Some methods depend on the availability of a network connection.

The result of the previous command will be:

Computing position of body 'mars' at epoch: jtd = 807954.6909685184 
Method 'Horizons':
	Position Epoch:  -2500-01-01 12:00:00.0 807953.9999999999
	Coordinates @ J2000: 
		Equatorial: 12:31:49.147 01:37:6.708
		Ecliptic: 186:39:54.558 04:38:33.609
	Coordinates @ Epoch : 
		Equatorial: 08:32:9.283 24:06:29.556
		Ecliptic: 124:21:14.472 04:39:4.619
	Observing conditions: 
		Distance to site [au]:  0.66052182424896
		Distance to sun [au]:  1.626124866723
		Solar elongation [deg]:  157:47:51.000
		Phase angle [deg]:  13:22:14.880
		Magnitude:  -1.1
	Other properties: 
		Local true sidereal time:  20:52:25.323
		Hour angle @ Epoch:  12:20:16.054
		Local coordinates @ Epoch:  06:11:33.727 -41:38:29.317
Method 'VSOP87':
	Position Epoch:  -2500-01-01 12:00:00.0 807953.9999999999
	Coordinates @ J2000: 
		Equatorial: 12:31:48.360 01:37:33.109
		Ecliptic: 186:39:33.205 04:38:53.193
	Coordinates @ Epoch : 
		Equatorial: 08:32:9.718 24:06:40.920
		Ecliptic: 124:21:17.390 04:39:17.040
	Observing conditions: 
		Distance to site [au]:  0.6604883074760437
		Distance to sun [au]:  1.626107096672058
		Solar elongation [deg]:  157:49:6.094
		Phase angle [deg]:  13:21:48.632
		Magnitude:  -1.13
	Other properties: 
		Local true sidereal time:  20:52:25.323
		Hour angle @ Epoch:  12:20:15.604
		Local coordinates @ Epoch:  06:11:23.906 -41:38:18.685
Method 'SPICE':
	Position Epoch:  -2500-01-01 12:00:00.0 807953.9999999999
	Coordinates @ J2000: 
		Equatorial: 12:31:48.754 01:37:12.184
		Ecliptic: 186:39:46.949 04:38:36.308
	Coordinates @ Epoch : 
		Equatorial: 08:32:9.796 24:06:28.555
		Ecliptic: 124:21:21.542 04:39:5.339
	Observing conditions: 
		Distance to site [au]:  0.660450348841685
		Distance to sun [au]:  1.6261149729985933
		Solar elongation [deg]:  157:49:18.876
		Phase angle [deg]:  13:21:31.981
		Magnitude:  -1.1
	Other properties: 
		Local true sidereal time:  20:52:25.323
		Hour angle @ Epoch:  12:20:15.527
		Local coordinates @ Epoch:  06:11:24.275 -41:38:31.094

Create a map of the sky

Now we want to create a map of the sky surrounding Aldebaran. For that purpose we get first the information about Aldebaran:

aldebaran = allstars.get_stars(ProperName='Aldebaran')

Now we select all stars between magnitude -1 and 5 around Aldebaran in a radius of 10 degrees:

hyades = allstars.get_stars_area(RA=aldebaran.data.RA,Dec=aldebaran.data.Dec,radius=10,Mag=[-1,5])

Plot the selected stars:

fig,ax = hyades.plot_stars(pad=0.0,labels=False,figargs=dict(figsize=(8,8)))

The resulting figure will be:

What's new

Versions 0.5.*:

• New properties (distances, phase angle, magnituded) added to PlanetaryBody class method.
• Solved the DEBUG problem.
• First fully-fledge working script, yay!
• MonTime class is now able to produce dates in all relevant calendars and formats.
• Verified ephemerides calculations for planets.
• Major improvements of all functionalities.

Version 0.1.*:

• First classes created and tested with study case.
• A proper identifying image of the project has been found.
• The project is started!

---

This package has been designed and written by Jorge I. Zuluaga with the historical advise of egyptologist Francisco "Tito" Vivas (C) 2023