Last active
September 30, 2021 22:10
-
-
Save bamford/a55afccbb68d0ed12b1dc5ebc5af37cd to your computer and use it in GitHub Desktop.
Eclipsing Binaries
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| #!/bin/python | |
| from astropy import time, coordinates as coord, units as u, table as tab | |
| from astropy.table import Table | |
| from astropy.time import Time, TimeDelta | |
| import numpy as np | |
| path = 'eclipses' | |
| t = Table.read('eclipsing_binaries.fits') | |
| t.write('{}/binaries.html'.format(path), overwrite=True, | |
| htmldict={'cssfiles': ['./eclipseTable.css'], | |
| 'table_class': 'eclipseTable'}) | |
| # Now calculate times of next minima, and select objects suitable for observation... | |
| t['coord'] = coord.SkyCoord(t['ra'], t['dec'], unit=('hour', 'deg')) | |
| # compute next minima | |
| Mnow = Time.now().jd | |
| t['epoch'] = (np.ceil((Mnow - t['M0']) / t['P'])).astype(np.int) | |
| t['Mnext'] = t['M0'] + t['epoch'] * t['P'] | |
| # compute local utc of next minima (correct from heliocentric time) | |
| nottingham = coord.EarthLocation(lon=-1.19201, lat=+52.94166, height=35.0) | |
| times_hjd = time.Time(t['Mnext'], format='jd', scale='utc', location=nottingham) | |
| ltt_helio = times_hjd.light_travel_time(t['coord'], 'heliocentric') | |
| times_jd = times_hjd - ltt_helio | |
| t['Mnext_UTC'] = times_jd.iso | |
| # This is a worse-case-scenario error, reality may not be so bad: | |
| t['Mnext_err_min'] = ((t['M0_err'] + t['epoch'] * t['P_err']) * 24 * 60).round(1) | |
| # compute sky position | |
| altaz = t['coord'].transform_to(coord.AltAz(obstime=times_jd, location=nottingham)) | |
| t['alt'], t['az'] = altaz.alt.round(1), altaz.az.round(1) | |
| # Get all EA/B with primary eclipse between 5pm and 10pm UTC tonight with duration less than 3 hours, altitude above 30 deg and brighter than V=12... | |
| Tearly = Time(Time.now().value.replace(hour=17, minute=0, second = 0)) | |
| Tlate = Time(Time.now().value.replace(hour=22, minute=0, second = 0)) | |
| sel = t['minimum_type'] != 'SEC' | |
| sel &= t['Mnext'] > Tearly.jd | |
| sel &= t['Mnext'] < Tlate.jd | |
| sel &= [('EA' in x)|('EB' in x) for x in t['var_type']] | |
| sel &= (t['D_day'] > 0.01) & (t['D_day'] < 3/24.) | |
| sel &= t['alt'] > 30 | |
| sel &= t['V'] < 12 | |
| useful_columns = ['constellation', 'id', 'var_type', 'V', 'Mnext_UTC', | |
| 'Mnext_err_min', 'P', 'D_min', 'alt', 'az', 'ra', 'dec', 'url'] | |
| tout = t[sel][useful_columns] | |
| # Use the url to see the O-C diagram. A positive O-C means the minimum is later than predicted. Use the left axis to find the O-C in days. 1 hour = 0.04 days. | |
| fn = '{}/tonight'.format(path) | |
| tout.write(fn+'.fits', overwrite=True) | |
| tout.write(fn+'.html', overwrite=True, | |
| htmldict={'cssfiles': ['./eclipseTable.css'], | |
| 'table_class': 'eclipseTable'}) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| { | |
| "cells": [ | |
| { | |
| "metadata": {}, | |
| "cell_type": "markdown", | |
| "source": "This notebook parses the tables prepared by J.M. Kreiner from the [Mt Suhora Astronomical Observatory](http://www.as.up.krakow.pl/ephem/) into a convenient format. The resulting combined FITS table is [available here](https://www.dropbox.com/s/jrmhnk6vgezq63f/eclipsing_binaries.fits?dl=0). The end of this notebook also provides an example of calculating times of next minima and selecting suitable targets for observation. Note that these are linear ephemerides: O-C offsets need to be considered." | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "from astropy import time, coordinates as coord, units as u, table as tab\nfrom astropy.table import Table\nfrom astropy.time import Time, TimeDelta\nfrom uncertainties.core import parse_error_in_parentheses", | |
| "execution_count": 1, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "trusted": true, | |
| "collapsed": true | |
| }, | |
| "cell_type": "code", | |
| "source": "# Process first table\n!wget -qr http://www.as.up.krakow.pl/ephem/EPHEM.TXT", | |
| "execution_count": 2, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true, | |
| "hide_input": false | |
| }, | |
| "cell_type": "code", | |
| "source": "h1 = 'con id min M0 X period Y A R S E C V G P years date D GCVS'\n#h2 = '+ -+ -+ -+ -+- + -+- + -+ -+ -+ -+ -+ -+ -+ -+ -+ -+ -+ -'\nh2 = '+ -+ -+ -+ -+ -+ -+ -+ -+ -+ -+ -+ -+ -+ -+ -+ -+ -'\n# AND RT ALL 2452500.3511(1) 0.62892857 (5) 132 96 36 37 95 0 0 0 1997-2013 17.Mar.14 0.170 10012", | |
| "execution_count": 3, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "names = [\n 'constellation',\n 'id',\n 'minimum_type',\n 'M0_str',\n 'P_str',\n 'n_min_all',\n 'n_min_pri',\n 'n_min_sec',\n 'n_min_e',\n 'n_min_ccd',\n 'n_min_v',\n 'n_min_pg',\n 'n_min_p',\n 'years',\n 'date',\n 'D',\n 'GCVS']\nh2 = np.array(list(h2))", | |
| "execution_count": 4, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "starts = tuple(np.where(np.array(list(h2)) == '+')[0])\nends = tuple(np.where(np.array(list(h2)) == '-')[0] - 1)", | |
| "execution_count": 5, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "t1 = Table.read('EPHEM.TXT', format='ascii.fixed_width_no_header',\n names=names, col_starts=starts,\n col_ends=ends, data_start=2)", | |
| "execution_count": 6, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "M0 = np.array([parse_error_in_parentheses(x.replace(' ', '').replace('(*)', '(100)').replace('()', '(100)'))\n for x in t1['M0_str']])\nt1['M0'], t1['M0_err'] = M0.T\nt1.remove_column('M0_str')\n\nP = np.array([parse_error_in_parentheses(x.replace(' ', '').replace('(*)', '(100)').replace('()', '(100)'))\n for x in t1['P_str']])\nt1['P'], t1['P_err'] = P.T\nt1.remove_column('P_str')", | |
| "execution_count": 7, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "trusted": true, | |
| "collapsed": true | |
| }, | |
| "cell_type": "code", | |
| "source": "# Process second table\n!wget -qr http://www.as.up.krakow.pl/ephem/allstars-cat.txt", | |
| "execution_count": 8, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "t2 = Table(names=('id', 'constellation', 'lastyear', 'V', 'spec_type', 'var_type', 'minimum_type',\n 'ra_h', 'ra_m', 'ra_s', 'dec_d', 'dec_m', 'dec_s', 'D', 'P', 'M0'),\n dtype=('U5', 'U3', 'i4', 'f4', 'U9', 'U11', 'U3', 'i2', 'i2', 'f4', 'i2', 'i2', 'f4', 'f4',\n 'f4', 'f8'))", | |
| "execution_count": 9, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "with open('allstars-cat.txt') as f:\n start = False\n for line in f:\n if start:\n line = line.replace('| ', '?').replace('|', '').replace(' : ', ' ')\n ls = line.split()\n if odd:\n #star, const, year of last minimum observations, Vmax, Vmin, spectral type, variability type, elements for ALL/PRI or SEC minima\n #RT And 2013 8.97 |F8V+K1 EA/RS for ALL\n try:\n id, constellation, lastyear, V, spec_type, var_type, _, minimum_type = ls\n except:\n print(line)\n odd = False\n else:\n #RA(hh mm ss.s), DEC(dd '' \"\"), D(in days; form GCVS), d (days), Period, Mo(HJD), phase of secondary\n #23 11 10 53 1 33 2000.0 0.11 0.0 0.6289285 2452500.3511 0.5\n try:\n ra_h, ra_m, ra_s, dec_d, dec_m, dec_s, _, D, _, P, M0, _ = ls\n except:\n print(line)\n t2.add_row((id.upper(), constellation.upper(), lastyear, V, spec_type, var_type, minimum_type,\n ra_h, ra_m, ra_s, dec_d, dec_m, dec_s, D, P, M0))\n odd = True\n elif line.startswith('^^^^^'):\n start = True\n odd = True", | |
| "execution_count": 10, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "# Make everything uppercase\nfor t in (t1, t2):\n for col in t.colnames:\n if t[col].dtype.kind == 'U':\n t[col] = [x.upper() for x in t[col]]", | |
| "execution_count": 11, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "# Remove perfect duplicate columns\nfor col in t1.colnames:\n if col in t2.colnames:\n t2.replace_column(col, t2[col].astype(t1[col].dtype))\n if (np.array(t1[col]) == (np.array(t2[col]))).all():\n t2.remove_column(col)", | |
| "execution_count": 12, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "# combine tables\nt = tab.hstack((t1, t2))", | |
| "execution_count": 13, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "t['ra'] = ['{:02d}:{:02d}:{:02.0f}'.format(x['ra_h'], x['ra_m'], x['ra_s']) for x in t]\nt['dec'] = ['{:02d}:{:02d}:{:02.0f}'.format(x['dec_d'], x['dec_m'], x['dec_s']) for x in t]\nfor x in ('ra_h', 'ra_m', 'ra_s', 'dec_d', 'dec_m', 'dec_s'):\n t.remove_column(x)", | |
| "execution_count": 14, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "t.rename_column('D_2', 'D_day')\nt.rename_column('D_1', 'D_period')\nt.rename_column('P_1', 'P')\nt.remove_column('P_2')", | |
| "execution_count": 15, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "t['D_min'] = (t['D_day'] * 24 * 60).round(1)", | |
| "execution_count": 16, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "t['url'] = ['http://www.as.up.krakow.pl/minicalc/' + x + y + '.HTM' for x, y in zip(t['constellation'], t['id'])]", | |
| "execution_count": 17, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "trusted": true, | |
| "collapsed": true | |
| }, | |
| "cell_type": "code", | |
| "source": "t.write('eclipsing_binaries.fits', overwrite=True)", | |
| "execution_count": 18, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true, | |
| "hide_input": false | |
| }, | |
| "cell_type": "code", | |
| "source": "# Now calculate times of next minima, and select objects suitable for observation...\nt['coord'] = coord.SkyCoord(t['ra'], t['dec'], unit=('hour', 'deg'))", | |
| "execution_count": 19, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "collapsed": true, | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "# compute next minima\nMnow = Time.now().jd\nt['epoch'] = (np.ceil((Mnow - t['M0']) / t['P'])).astype(np.int)\nt['Mnext'] = t['M0'] + t['epoch'] * t['P']", | |
| "execution_count": 20, | |
| "outputs": [] | |
| }, | |
| { | |
| "metadata": { | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "# compute local utc of next minima (correct from heliocentric time)\nnottingham = coord.EarthLocation(lon=-1.19201, lat=+52.94166, height=35.0)\ntimes_hjd = time.Time(t['Mnext'], format='jd', scale='utc', location=nottingham)\nltt_helio = times_hjd.light_travel_time(t['coord'], 'heliocentric')\ntimes_jd = times_hjd - ltt_helio\nt['Mnext_UTC'] = times_jd.iso\n# This is a worse-case-scenario error, reality may not be so bad:\nt['Mnext_err_min'] = ((t['M0_err'] + t['epoch'] * t['P_err']) * 24 * 60).round(1)", | |
| "execution_count": 21, | |
| "outputs": [ | |
| { | |
| "output_type": "stream", | |
| "text": "WARNING: Tried to get polar motions for times after IERS data is valid. Defaulting to polar motion from the 50-yr mean for those. This may affect precision at the 10s of arcsec level [astropy.coordinates.builtin_frames.utils]\n", | |
| "name": "stderr" | |
| } | |
| ] | |
| }, | |
| { | |
| "metadata": { | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "# compute sky position\naltaz = t['coord'].transform_to(coord.AltAz(obstime=times_jd, location=nottingham))\nt['alt'], t['az'] = altaz.alt.round(1), altaz.az.round(1)", | |
| "execution_count": 22, | |
| "outputs": [ | |
| { | |
| "output_type": "stream", | |
| "text": "WARNING: Tried to get polar motions for times after IERS data is valid. Defaulting to polar motion from the 50-yr mean for those. This may affect precision at the 10s of arcsec level [astropy.coordinates.builtin_frames.utils]\n", | |
| "name": "stderr" | |
| } | |
| ] | |
| }, | |
| { | |
| "metadata": {}, | |
| "cell_type": "markdown", | |
| "source": "Get all EA/B with primary eclipse between 6pm and 10pm UTC (7pm-11pm BST) on 2017-10-27 with duration less than 3 hours, altitude above 30 deg and brighter than V=12..." | |
| }, | |
| { | |
| "metadata": { | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "Tearly = Time('2017-10-27T18:00:00')\nTlate = Time('2017-10-27T22:00:00')\nsel = t['minimum_type'] != 'SEC'\nsel &= t['Mnext'] > Tearly.jd\nsel &= t['Mnext'] < Tlate.jd\nsel &= [('EA' in x)|('EB' in x) for x in t['var_type']]\nsel &= (t['D_day'] > 0.01) & (t['D_day'] < 3/24.)\nsel &= t['alt'] > 30 * u.deg\nsel &= t['V'] < 12\nuseful_columns = ['constellation', 'id', 'var_type', 'V', 'Mnext_UTC',\n 'Mnext_err_min', 'P', 'D_min', 'alt', 'az', 'ra', 'dec', 'url']\nt[sel][useful_columns]", | |
| "execution_count": 23, | |
| "outputs": [ | |
| { | |
| "output_type": "execute_result", | |
| "execution_count": 23, | |
| "data": { | |
| "text/plain": "<Table length=5>\nconstellation id ... dec url \n ... \n str3 str5 ... str9 str48 \n------------- ---- ... -------- ---------------------------------------------\n CAS CC ... 59:33:49 http://www.as.up.krakow.pl/minicalc/CASCC.HTM\n CEP WW ... 69:51:40 http://www.as.up.krakow.pl/minicalc/CEPWW.HTM\n CEP XZ ... 67:09:03 http://www.as.up.krakow.pl/minicalc/CEPXZ.HTM\n CEP AI ... 56:55:02 http://www.as.up.krakow.pl/minicalc/CEPAI.HTM\n CYG KV ... 36:47:37 http://www.as.up.krakow.pl/minicalc/CYGKV.HTM", | |
| "text/html": "<Table length=5>\n<table id=\"table4605045504\" class=\"table-striped table-bordered table-condensed\">\n<thead><tr><th>constellation</th><th>id</th><th>var_type</th><th>V</th><th>Mnext_UTC</th><th>Mnext_err_min</th><th>P</th><th>D_min</th><th>alt</th><th>az</th><th>ra</th><th>dec</th><th>url</th></tr></thead>\n<thead><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>deg</th><th>deg</th><th></th><th></th><th></th></tr></thead>\n<thead><tr><th>str3</th><th>str5</th><th>str10</th><th>float32</th><th>str23</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>str8</th><th>str9</th><th>str48</th></tr></thead>\n<tr><td>CAS</td><td>CC</td><td>EB/DM</td><td>7.06</td><td>2017-10-27 20:50:00.598</td><td>45.4</td><td>3.366308</td><td>43.2</td><td>56.7</td><td>53.7</td><td>03:14:05</td><td>59:33:49</td><td>http://www.as.up.krakow.pl/minicalc/CASCC.HTM</td></tr>\n<tr><td>CEP</td><td>WW</td><td>EA/SD:</td><td>11.1</td><td>2017-10-27 21:19:08.369</td><td>4.9</td><td>4.600848</td><td>72.0</td><td>70.6</td><td>339.0</td><td>22:18:27</td><td>69:51:40</td><td>http://www.as.up.krakow.pl/minicalc/CEPWW.HTM</td></tr>\n<tr><td>CEP</td><td>XZ</td><td>EB/DM:</td><td>8.0</td><td>2017-10-27 20:29:27.118</td><td>4.4</td><td>5.09725</td><td>72.0</td><td>75.5</td><td>353.3</td><td>22:32:25</td><td>67:09:03</td><td>http://www.as.up.krakow.pl/minicalc/CEPXZ.HTM</td></tr>\n<tr><td>CEP</td><td>AI</td><td>EB/DM</td><td>9.18</td><td>2017-10-27 18:41:04.954</td><td>17.1</td><td>4.225334</td><td>57.6</td><td>82.3</td><td>53.5</td><td>21:46:22</td><td>56:55:02</td><td>http://www.as.up.krakow.pl/minicalc/CEPAI.HTM</td></tr>\n<tr><td>CYG</td><td>KV</td><td>EB/SD</td><td>11.5</td><td>2017-10-27 18:39:01.694</td><td>4.0</td><td>2.8390041</td><td>43.2</td><td>72.2</td><td>209.4</td><td>20:15:38</td><td>36:47:37</td><td>http://www.as.up.krakow.pl/minicalc/CYGKV.HTM</td></tr>\n</table>" | |
| }, | |
| "metadata": {} | |
| } | |
| ] | |
| }, | |
| { | |
| "metadata": {}, | |
| "cell_type": "markdown", | |
| "source": "Use the url to see the O-C diagram. A positive O-C means the minimum is later than predicted. Use the left axis to find the O-C in days. 1 hour = 0.04 days." | |
| }, | |
| { | |
| "metadata": {}, | |
| "cell_type": "markdown", | |
| "source": "All the columns, in case they provide additional insight..." | |
| }, | |
| { | |
| "metadata": { | |
| "trusted": true | |
| }, | |
| "cell_type": "code", | |
| "source": "t[sel]", | |
| "execution_count": 24, | |
| "outputs": [ | |
| { | |
| "output_type": "execute_result", | |
| "execution_count": 24, | |
| "data": { | |
| "text/plain": "<Table length=5>\nconstellation id minimum_type n_min_all ... Mnext_err_min alt az \n ... deg deg \n str3 str5 str3 int64 ... float64 float64 float64\n------------- ---- ------------ --------- ... ------------- ------- -------\n CAS CC ALL 14 ... 45.4 56.7 53.7\n CEP WW ALL 33 ... 4.9 70.6 339.0\n CEP XZ ALL 13 ... 4.4 75.5 353.3\n CEP AI ALL 8 ... 17.1 82.3 53.5\n CYG KV PRI 17 ... 4.0 72.2 209.4", | |
| "text/html": "<Table length=5>\n<table id=\"table4604947256\" class=\"table-striped table-bordered table-condensed\">\n<thead><tr><th>constellation</th><th>id</th><th>minimum_type</th><th>n_min_all</th><th>n_min_pri</th><th>n_min_sec</th><th>n_min_e</th><th>n_min_ccd</th><th>n_min_v</th><th>n_min_pg</th><th>n_min_p</th><th>years</th><th>date</th><th>D_period</th><th>GCVS</th><th>M0</th><th>M0_err</th><th>P</th><th>P_err</th><th>lastyear</th><th>V</th><th>spec_type</th><th>var_type</th><th>D_day</th><th>ra</th><th>dec</th><th>D_min</th><th>url</th><th>coord</th><th>epoch</th><th>Mnext</th><th>Mnext_UTC</th><th>Mnext_err_min</th><th>alt</th><th>az</th></tr></thead>\n<thead><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>deg,deg</th><th></th><th></th><th></th><th></th><th>deg</th><th>deg</th></tr></thead>\n<thead><tr><th>str3</th><th>str5</th><th>str3</th><th>int64</th><th>int64</th><th>int64</th><th>int64</th><th>int64</th><th>int64</th><th>int64</th><th>int64</th><th>str9</th><th>str9</th><th>float64</th><th>int64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>int32</th><th>float32</th><th>str8</th><th>str10</th><th>float64</th><th>str8</th><th>str9</th><th>float64</th><th>str48</th><th>object</th><th>int64</th><th>float64</th><th>str23</th><th>float64</th><th>float64</th><th>float64</th></tr></thead>\n<tr><td>CAS</td><td>CC</td><td>ALL</td><td>14</td><td>10</td><td>4</td><td>7</td><td>7</td><td>0</td><td>0</td><td>0</td><td>1972-2010</td><td>30.JUN.11</td><td>0.01</td><td>80104</td><td>2452503.33</td><td>0.02</td><td>3.366308</td><td>7e-06</td><td>2010</td><td>7.06</td><td>O9IV+O9I</td><td>EB/DM</td><td>0.0299999993294</td><td>03:14:05</td><td>59:33:49</td><td>43.2</td><td>http://www.as.up.krakow.pl/minicalc/CASCC.HTM</td><td>48.52083333333333,59.56361111111111</td><td>1649</td><td>2458054.37189</td><td>2017-10-27 20:50:00.598</td><td>45.4</td><td>56.7</td><td>53.7</td></tr>\n<tr><td>CEP</td><td>WW</td><td>ALL</td><td>33</td><td>18</td><td>15</td><td>9</td><td>19</td><td>5</td><td>0</td><td>0</td><td>1994-2007</td><td>18.FEB.09</td><td>0.01</td><td>45</td><td>2452501.167</td><td>0.001</td><td>4.600848</td><td>2e-06</td><td>2007</td><td>11.1</td><td>G3</td><td>EA/SD:</td><td>0.0500000007451</td><td>22:18:27</td><td>69:51:40</td><td>72.0</td><td>http://www.as.up.krakow.pl/minicalc/CEPWW.HTM</td><td>334.61249999999995,69.8611111111111</td><td>1207</td><td>2458054.39054</td><td>2017-10-27 21:19:08.369</td><td>4.9</td><td>70.6</td><td>339.0</td></tr>\n<tr><td>CEP</td><td>XZ</td><td>ALL</td><td>13</td><td>12</td><td>1</td><td>9</td><td>4</td><td>0</td><td>0</td><td>0</td><td>1971-2012</td><td>17.DEC.13</td><td>0.01</td><td>51</td><td>2452503.451</td><td>0.002</td><td>5.09725</td><td>1e-06</td><td>2012</td><td>8.0</td><td>O9.5V</td><td>EB/DM:</td><td>0.0500000007451</td><td>22:32:25</td><td>67:09:03</td><td>72.0</td><td>http://www.as.up.krakow.pl/minicalc/CEPXZ.HTM</td><td>338.10416666666663,67.15083333333334</td><td>1089</td><td>2458054.35625</td><td>2017-10-27 20:29:27.118</td><td>4.4</td><td>75.5</td><td>353.3</td></tr>\n<tr><td>CEP</td><td>AI</td><td>ALL</td><td>8</td><td>4</td><td>4</td><td>0</td><td>8</td><td>0</td><td>0</td><td>0</td><td>1991-2009</td><td>24.NOV.14</td><td>0.01</td><td>63</td><td>2452502.192</td><td>0.004</td><td>4.225334</td><td>6e-06</td><td>2009</td><td>9.18</td><td>B0.5V:P</td><td>EB/DM</td><td>0.0399999991059</td><td>21:46:22</td><td>56:55:02</td><td>57.6</td><td>http://www.as.up.krakow.pl/minicalc/CEPAI.HTM</td><td>326.5916666666666,56.91722222222222</td><td>1314</td><td>2458054.28088</td><td>2017-10-27 18:41:04.954</td><td>17.1</td><td>82.3</td><td>53.5</td></tr>\n<tr><td>CYG</td><td>KV</td><td>PRI</td><td>17</td><td>17</td><td>0</td><td>4</td><td>13</td><td>0</td><td>0</td><td>0</td><td>1990-2013</td><td>21.NOV.13</td><td>0.01</td><td>10255</td><td>2452501.186</td><td>0.001</td><td>2.8390041</td><td>9e-07</td><td>2013</td><td>11.5</td><td>B0</td><td>EB/SD</td><td>0.0299999993294</td><td>20:15:38</td><td>36:47:37</td><td>43.2</td><td>http://www.as.up.krakow.pl/minicalc/CYGKV.HTM</td><td>303.9083333333333,36.79361111111111</td><td>1956</td><td>2458054.27802</td><td>2017-10-27 18:39:01.694</td><td>4.0</td><td>72.2</td><td>209.4</td></tr>\n</table>" | |
| }, | |
| "metadata": {} | |
| } | |
| ] | |
| } | |
| ], | |
| "metadata": { | |
| "_draft": { | |
| "nbviewer_url": "https://gist.github.com/a55afccbb68d0ed12b1dc5ebc5af37cd" | |
| }, | |
| "gist": { | |
| "id": "a55afccbb68d0ed12b1dc5ebc5af37cd", | |
| "data": { | |
| "description": "Eclipsing Binaries ", | |
| "public": true | |
| } | |
| }, | |
| "hide_input": false, | |
| "kernelspec": { | |
| "name": "python3", | |
| "display_name": "Python [default]", | |
| "language": "python" | |
| }, | |
| "language_info": { | |
| "name": "python", | |
| "version": "3.6.1", | |
| "mimetype": "text/x-python", | |
| "codemirror_mode": { | |
| "name": "ipython", | |
| "version": 3 | |
| }, | |
| "pygments_lexer": "ipython3", | |
| "nbconvert_exporter": "python", | |
| "file_extension": ".py" | |
| }, | |
| "latex_envs": { | |
| "eqNumInitial": 1, | |
| "eqLabelWithNumbers": true, | |
| "current_citInitial": 1, | |
| "cite_by": "apalike", | |
| "bibliofile": "biblio.bib", | |
| "LaTeX_envs_menu_present": true, | |
| "labels_anchors": false, | |
| "latex_user_defs": false, | |
| "user_envs_cfg": false, | |
| "report_style_numbering": false, | |
| "autocomplete": true, | |
| "hotkeys": { | |
| "equation": "Ctrl-E", | |
| "itemize": "Ctrl-I" | |
| } | |
| } | |
| }, | |
| "nbformat": 4, | |
| "nbformat_minor": 1 | |
| } |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
eclipsing_binaries.pygenerates the nightly table of eclipses here.