sunstardb short example.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "# Load the necessary modules\n",
    "import sunstardb.database\n",
    "import astropy.table\n",
    "import os\n",
    "from matplotlib import pyplot as plt\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Connect to the database\n",
    "db = sunstardb.database.SunStarDB()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Fetch the cycle periods from the Baliunas1995 dataset\n",
    "cycs = db.fetch_data_table('Baliunas1995', ['Pcyc'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Fetch the rotation periods from the Donahue1996 dataset\n",
    "rots = db.fetch_data_table('Donahue1996', ['Prot'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<i>Table length=23</i>\n",
       "<table id=\"table4513717160-782451\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>idx</th><th>star</th><th>Pcyc</th><th>Prot</th></tr></thead>\n",
       "<tr><td>0</td><td>* 107 Psc</td><td>9.6</td><td>35.2</td></tr>\n",
       "<tr><td>1</td><td>* 11 LMi</td><td>7.9</td><td>18.6</td></tr>\n",
       "<tr><td>2</td><td>* 12 Oph</td><td>17.4</td><td>21.07</td></tr>\n",
       "<tr><td>3</td><td>* 15 Sge</td><td>2.6</td><td>13.94</td></tr>\n",
       "<tr><td>4</td><td>* 36 Oph B</td><td>5.7</td><td>21.11</td></tr>\n",
       "<tr><td>5</td><td>* 61 Cyg A</td><td>7.3</td><td>35.37</td></tr>\n",
       "<tr><td>6</td><td>* 61 Cyg B</td><td>11.7</td><td>37.84</td></tr>\n",
       "<tr><td>7</td><td>* 9 Cet</td><td>9.1</td><td>7.78</td></tr>\n",
       "<tr><td>8</td><td>* bet Com</td><td>16.6</td><td>12.35</td></tr>\n",
       "<tr><td>9</td><td>* kap01 Cet</td><td>5.6</td><td>9.24</td></tr>\n",
       "<tr><td>10</td><td>HD 115404</td><td>12.4</td><td>18.47</td></tr>\n",
       "<tr><td>11</td><td>HD 152391</td><td>10.9</td><td>11.43</td></tr>\n",
       "<tr><td>12</td><td>HD 160346</td><td>7.0</td><td>36.4</td></tr>\n",
       "<tr><td>13</td><td>HD 16160</td><td>13.2</td><td>48.0</td></tr>\n",
       "<tr><td>14</td><td>HD 190007</td><td>13.7</td><td>28.95</td></tr>\n",
       "<tr><td>15</td><td>HD 4628</td><td>8.37</td><td>38.5</td></tr>\n",
       "<tr><td>16</td><td>HD 78366</td><td>12.2</td><td>9.67</td></tr>\n",
       "<tr><td>17</td><td>HR 1321</td><td>7.8</td><td>7.15</td></tr>\n",
       "<tr><td>18</td><td>HR 3750</td><td>8.17</td><td>40.2</td></tr>\n",
       "<tr><td>19</td><td>HR 6806</td><td>15.8</td><td>42.4</td></tr>\n",
       "<tr><td>20</td><td>Sun</td><td>10.0</td><td>26.09</td></tr>\n",
       "<tr><td>21</td><td>V* HN Peg</td><td>6.2</td><td>4.86</td></tr>\n",
       "<tr><td>22</td><td>V* V2213 Oph</td><td>7.4</td><td>7.78</td></tr>\n",
       "</table><style>table.dataTable {clear: both; width: auto !important; margin: 0 !important;}\n",
       ".dataTables_info, .dataTables_length, .dataTables_filter, .dataTables_paginate{\n",
       "display: inline-block; margin-right: 1em; }\n",
       ".paginate_button { margin-right: 5px; }\n",
       "</style>\n",
       "<script>\n",
       "\n",
       "var astropy_sort_num = function(a, b) {\n",
       "    var a_num = parseFloat(a);\n",
       "    var b_num = parseFloat(b);\n",
       "\n",
       "    if (isNaN(a_num) && isNaN(b_num))\n",
       "        return ((a < b) ? -1 : ((a > b) ? 1 : 0));\n",
       "    else if (!isNaN(a_num) && !isNaN(b_num))\n",
       "        return ((a_num < b_num) ? -1 : ((a_num > b_num) ? 1 : 0));\n",
       "    else\n",
       "        return isNaN(a_num) ? -1 : 1;\n",
       "}\n",
       "\n",
       "require.config({paths: {\n",
       "    datatables: 'https://cdn.datatables.net/1.10.12/js/jquery.dataTables.min'\n",
       "}});\n",
       "require([\"datatables\"], function(){\n",
       "    console.log(\"$('#table4513717160-782451').dataTable()\");\n",
       "    \n",
       "jQuery.extend( jQuery.fn.dataTableExt.oSort, {\n",
       "    \"optionalnum-asc\": astropy_sort_num,\n",
       "    \"optionalnum-desc\": function (a,b) { return -astropy_sort_num(a, b); }\n",
       "});\n",
       "\n",
       "    $('#table4513717160-782451').dataTable({\n",
       "        order: [],\n",
       "        pageLength: 5,\n",
       "        lengthMenu: [[5, 10, 25, 50, 100, 500, 1000, -1], [5, 10, 25, 50, 100, 500, 1000, 'All']],\n",
       "        pagingType: \"full_numbers\",\n",
       "        columnDefs: [{targets: [0, 2, 3], type: \"optionalnum\"}]\n",
       "    });\n",
       "});\n",
       "</script>\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Join the cycle and rotation data on the star name and display the results\n",
    "cycrots = astropy.table.join(cycs, rots, keys=['star'])\n",
    "cycrots.show_in_notebook(display_length=5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5,1,'Cycle Period vs Rotation Period')"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAEXCAYAAACDJtDpAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAG/ZJREFUeJzt3XucXWV97/HPdxIUBhAICYqEmQGpFEQuEhTxUopiwVK0HmihAwctNi9pq0BtEcx5ofTVaRUVpKen53QUDmAGUKioh5cWcqCAVS4mETCASj1mAiaQcL9EuSS/88d6drKzM/s2s/Ze+/J9v177NXuvtWat396z9m+e9TzPeh5FBGZmeRkoOgAz6y1OKmaWKycVM8uVk4qZ5cpJxcxy5aRiZrlyUsmJpJWS3tsBcXxP0mnT/N2OeA+tNJPPp10kDUl6XtKsafzukZIeaUVcjerLpCLpTyQtTX+4NelEe2dBsRwpaWOK5TlJP5P0kenuLyKOjYgr8owxLzN9r5JulfTRJrb/rKTF5cta9flIulzSS+m9PSlpiaTfns6+ImJVROwQERvyjrMd+i6pSPor4MvA3wOvBYaAfwY+UGBYqyNiB+A1wKeAr0jav5kdKNMNf8/y93o22Xvdt+CY8nJhem/zgbXA5c3uQNLsvINqu4jomwewE/A8cGKV9a8D1gO7li07FFgHbJNe/xnwIPAc8ADwlrR8JfDe9HwAOBf4BfAE8A1gTpVjHgk8UrFsHXBCen448EPgaeBe4Miy7W4FxoAfAL8G9knLPloWx38DJslO8iuBncp+/9S07glgUfl7qIjncOBRYFbZsj8E7kvP3wosBZ4FHgMuauK9ri3/ewBHAD8Cnkk/j0jLx4ANwG/S3/Cf0vJLgIfTsZcB70rLjwFeAl5O299b9pnV/XyAESCA04BVwOPAohrn1uXA35W9/n3g+XrnQ9lxTk/Hub1s2ey0zeuB7wBPAv8J/FnZcbZLx36K7Hz8m8rPuO3fs6K/6G19s9mJ9krpj1Vlm+8CZ5S9vhj47+n5icCvgMMAkX2Jh9O6TV9I4CzgTrL/WK8G/gW4usrxNn3R0sn3h+mLsC+wRzoJ35/WHZ1ezyv7gqwC3gTMBrap+NL8aToJ9wZ2AL4JfC2t2z992d6dYrwofTZbJZW0/S+Ao8teXwucm57fAZyanu8AHN7gez0e2AgckpbNSV+OU9P7OTm93rXs/X60Yp+nALum7T9Jlvy2Tes+Cyyu2L7Rz2eE7Iv9FbIv7kHAi8B+Vd7b5aSkkvZ1FfD9eudD2XGuBLZPxyotKyWV28hK09sCB5P903lPWvc54Pvps9sTWIGTSluTyijwaJ1t/hj4QXo+K52kb02vbwTOrPJ7K9mcVB4s/dHT693JEsVWySx90TaSlUSeBO4BTkrrPlU6ycu2vxE4LTZ/Qf62Yn35l+Zm4M/L1u1bigM4H7imbN32ZP/ZqyWVvwMuS893BF5gc0K9HbgAmFvnsy1/ry+SlTzOKlt/KnB3xe/cAXy48r3VOMZTwEHp+WepnVRqfT6lL/b8svV3l/42Uxz3crJS1NPpnPkO8IZ650PZcfYuW19aNpssUWwAdixb/w/A5en5/wOOKVu3kIKTSjdcg+fpCWBunevWbwP7S9qbrGTwTETcndbtSfYfu55h4HpJT0t6muyk2kBWhzOV1RGxc0TMiYiDI+Kasv2cWNpP2tc7yU7KkodrxPF6sqJ9ySTZifratG7T70bEC2SfTzVXAR+S9GrgQ8DyiCjt+3TgjcBPJf1I0nE19rM6InYmq1P5R+CoGvGWYt6j2s4kfVLSg5KeSZ/PTsDcGscvV+vzKXm07Pl6slJINV9Mf8fXRcTxEVE6Vxo5H6r9HV8PPBkRz1XEuUfZ+ocr1hWq35LKHWT/TT5YbYOI+A3ZNe8o2X/Or5Wtfhh4QwPHeRg4Np1gpce2EfGrJuN9mKykUr6f7SPic+Uh1/j91WQndMkQ2SXOY8AasiQJgKRBssuIKUXEA2Qn7LHAn5AlmdK6hyLiZGA34PPAdZK2r/XGIuJFspLYmyWV/h6V8ZZiLn1uW7xXSe9K+/gjYJeUrJ4huzTdavsp1Pp88tTI+VAt1tXAHEk7VsRZ+t0t/o5pXaH6KqlExDNkxf7/IemDkgYlbSPpWEkXlm16JfBhsmv+8ibJrwJ/LenQ1Nqyj6TKLwHA/wLGSuskzZM0ndalxcAfSPo9SbMkbZuaZec3+PtXA2dL2kvSDmQtXl+PiFeA64DjJL1T0quAv6X++XAV8AmyephrSwslnSJpXkSULm0g+09cU0S8BHyJ7G8CWX3WG1OT/2xJf0xW93NDWv8YWf1HyY5kSWAdMFvS+WQlIMq2H6nRKlbr88nTtM+HiHiYrKL+H9Lf/0CykuFE2uQbwHmSdknnxcdzjr1pfZVUACLiIuCvyGr915H9F/lL4Ftl2/yA7Np/eUSsLFt+LVkrxFVkrT/fIqsgq3QJ2TX1TZKeI6uke9s0Yn2YrKn702Wx/g2N/90uIytp3Q78kqyU9vG07/uBv0jvZQ1ZXUS9TlNXk9WL3BIRj5ctPwa4X9LzZO/9pFTiazTGIUl/EBFPAMeRVbg+AZwDHFd2rEuAEyQ9JekfyeqXvgf8nKwU9Ru2vBQoJb4nJC2vcuwpP5+czfR8OJmsnmU1cD3wmYhYktZdQPbefwncxJYl60IoVe5YBUm3AFdFxFeLjsWsmzipTEHSYcASYM+KCjIzq6PvLn/qkXQF8H/JmjqdUMya5JKKmeWqpSUVSZdJWitpRdmygyXdKemedFPfW1sZg5m1V0tLKpLeTdYV/MqIOCAtuwm4OCK+J+n9wDkRcWS9fc2dOzdGRkZaFquZTW3ZsmWPR8S8Rrdv6R2REXG7pJHKxWzuS7ATWTNZXSMjIyxdujS/4MysIZKa6qVbxG3WZwE3Svoi2eXXEdU2lLSQ7F4GhoYK7yhoZg0oovXnDODsiNiTbDyNS6ttGBHjEbEgIhbMm9dw6cvMClREUjmN7BZzyHo8uqLWrIcUkVRWA7+Tnh8FPFRADGbWIi2tU5FUuldkbhqM9zNkI6ddkoYf+A2pzsTMekNLSyoRcXJE7B4R20TE/Ii4NCL+IyIOjYiDIuJtEbGslTHY1CYmJhgZGWFgYICRkREmJibq/5JZA7p/kF1r2sTEBAsXLmT9+vUATE5OsnBhVmAcHR0tMjTrAb73pw8tWrRoU0IpWb9+PYsWLSooIuslTip9aNWqVU0tN2uGk0ofqtaR0B0MLQ9OKn1obGyMwcHBLZYNDg4yNjZWUETWS5xU+tDo6Cjj4+MMDw8jieHhYcbHx11Ja7nomvFUFixYEL6h0Kz9JC2LiAWNbu+SipnlyknFzHLlpGJmuXJSaQF3gbd+5m76OXMXeOt3LqnkzF3grd85qeTMXeCt3zmp5Mxd4K3fOankzF3grd85qeTMXeCt37mbvpnV1FHd9Kea9jQt/7ikn0m6X9KFrYzBzNqr1Zc/lwPHlC+Q9LvAB4ADI+JNwBdbHIOZtVGrB76+HXiyYvEZwOci4sW0zdpWxmBm7VVERe0bgXdJukvSbZIOq7ahpIWSlkpaum7dujaGaNZ+vXJ7RxHd9GcDuwCHA4cB35C0d0xRYxwR48A4ZBW1bY3SrI166faOIkoqjwDfjMzdwEZgbgFxmHWMXrq9o4ik8i2y6U6R9EbgVcDjBcRh1jF66faOVjcpXw3cAewr6RFJpwOXAXunZuZrgNOmuvQx6ye9dHtHS+tUIuLkKqtOaeVxzbrN2NjYFnUq0L23d7ibvlkH6KXbO9xN38xq6qhu+mbWf5xUzCxXTipmlisnFTPLlZOKmeXKScXMcuWkYma5clIxs1w5qZhZrpxUzCxXTipmlisnFTPLlZOKmeXKScXMcuWkYma5KmSGwrTuryWFJA96bdZD2j5DIYCkPYGjge4b1dfMaipihkKAi4FzgO4Yds7MGtb2OhVJxwO/ioh7231sM2u9ts5QKGkQWAS8r8HtFwILoTunKjDrR+0uqbwB2Au4V9JKYD6wXNLrpto4IsYjYkFELJg3b14bw7Qi9Mpcwv2urSWViPgJsFvpdUosCyLCMxT2uV6aS7jfFTFDodlWemku4X5X1AyFpfUjrTy+dY9emku437lHrXWEXppLuN85qVhHGBsbY3BwcItl3TqXcL9zUrGO0EtzCXe6VreyeS5lsz5S2coGWYmwVgL3XMpmVlU7WtmcVMz6SDta2ZxUzPpIO1rZnFTM+kg7WtmcVMz6SDta2dz6Y2Y1ufXHzArlpGJmuXJSMbNcOan0CA9wZJ2irYM0WWt4gCPrJC6p9AAPcGSdxEmlAZ1+aeEBjqyTOKnUUbq0mJycJCI2XVp0UmLxAEfWSdo+7amkL0j6qaT7JF0vaedWxjBT3XBp4QGOrJMUMe3pEuCAiDgQ+DlwXotjmJFuuLTwAEfWSVo98PXtkkYqlt1U9vJO4IRWxjBTQ0NDTE5OTrm8k4yOjjqJWEcouk7lT4HvVVspaaGkpZKWrlu3ro1hbeZLC7PmFJZUJC0CXgGq1nh2wgyFvrSwSp3eGli4iGjpAxgBVlQsO41skrHBRvdz6KGHhlnRFi9eHIODgwFsegwODsbixYuLDq1lgKXRxHe+7SUVSccAnwKOj4j19bY36yTd0BpYtJoVtZLmNLCPjRHxdJXfvxo4Epgr6RHgM2StPa8GlkgCuDMiPtZM0GZF6YbWwKLVa/1ZnR6qsc0sYMqmkJh62tNLGwvNrPN0S2tgkepd/jwYEXtHxF7VHsAT7QjUrBO4NbC+eknl7Q3so5FtzHqCWwPrqztGraQB4L6IOKA9IU3NY9SaFSP3MWojYiNwryRfNJpZXY12098duF/S3cALpYURcXxLojKzrtVoUrmgpVGYWc9otPPbm8nqVW4rf7QysG7ibttmmzVaUnkd8CNJy4HLgBujXg1vn/D4sGZbaniGQmXdX98HfARYAHwDuDQiftG68Dbr1NafkZGRKTtDDQ8Ps3LlyvYHZJazls1QmEomj6bHK8AuwHWSLmw6yh7ibttmW2ooqUj6hKRlwIXAD4A3R8QZwKHAf2lhfB3P48OabanRkspc4EMR8XsRcW1EvAyb+rAc17LouoC7bZttqaGkEhHnR8TWFQfZugfzDam7uNu22ZZqVtRKWh4Rb6m5gwa2yUOnVtSa9bpmK2rrNSnvJ+m+WscDdmr0YGbW++olld9uYB8b8gjEzHpDzaRSrR7FzKyaoqfoMLMeU8S0p3MkLZH0UPq5SytjsGL4fqj+1Wjnt70kbVv2ervKmQeruJytpz09F7g5In4LuDm9th7SDZPaW+s0WlK5FthY9npDWlZTRNwOPFmx+APAFen5FcAHG4zBuoSnsehvjSaV2RHxUulFev6qaR7ztRGxJu1nDbBbtQ07YdpTa57vh+pvjSaVdZI2jfIm6QPA460JabPogGlPrXm+H6q/NZpUzgA+LWmVpFVkMwwunOYxH5O0O0D6uXaa+7EO5fuh+lujSeWXEXE4sD/wpog4YgbjqHyHbC5l0s9vT3M/1qG67X4ot1TlrJEJl4FfAl8A9mtmombgamAN8DLwCHA6sCtZq89D6eecRvblCdqtFfpxwvVm0eQE7Q2N/CZpR+AkslHfBsiGlLwmIp7NMb/V5BsKrRU8cl99LRn5LSKei4ivRMQRwDlkE62vkXSFpH2mGatZ4dxSlb9GO7/NknS8pOuBS4AvAXsD/wf4bgvjM2spt1Tlr9GK2ofIOq19ISIOiYiLIuKxiLgO+LfWhWfWWm6pyl+jSeXAiDg9In5YuSIiPpFzTG3lmv/+1m0tVV2hkdpcsu70O5e93gW4rJka4Zk+WtH645p/s/posvWnmZLK02WJ6CngkDyTWxF8j4pZ/hpNKgPlQxRImkPjsxt2LNf8m+Wv0cTwJeCHkq4ju0z4I6Dra7KGhoam7KPgmn+z6Wu0n8qVZJOGPQasI5sD6Gul9d060JJr/s3y1/AlTEQ8ADxQZfXNQMun6chbqYZ/0aJFrFq1iqGhIcbGxlzzbzYDDU/QXnMn0o8joqUVt+6mb1aMlk3QXsfMM5OZ9QSPpm9mucorqSin/Zh1FffI3lrNito0gv7HgH2AnwCXRsQrU2z6nhbEZtbRSrMGlDpQlmYNAPq6sr/eBO1fJxtg6fvAscBkRJzZpti24Ipa6zT9MhZL3hW1+0fEKRHxL8AJwLtmFJ1ZD2mkR3Y/Xh7VSyovl55Uuewx61v1xmLp10nV6iWVgyQ9mx7PAQeWnkua0VCSks6WdL+kFZKuLp8B0awb1OuR3a83rNZMKhExKyJekx47RsTssuevme5BJe0BfAJYEBEHALPIxsC1GejHonaR6o3F0q83rBZ5p/FsYDtJLwODwOoCY+l6bokoxujoaNXPt19vWC2k81tE/Ar4IrCKbAqPZyLipsrtWjXtaS/+R+/XonYn69sbVpsZ0SmvB9nIcbcA84BtgG8Bp9T6nbxGfuvV0d4kbfGeSg9JRYfW1xYvXhzDw8MhKYaHh7vyPKMV8/7kTdKJwDERcXp6/V+BwyPiz6v9Tl79VHq1b0Gvvi8rXlE3FDZrFXC4pEFJIuuR+2BbDtyjlWd9W9S2jlNUncpdwHXAcrLu/wPAeDuO3avzvHhUeOsUhVz+TEdelz+VrSQAkogIhoeHPUiTWYVuufwpTPl/dNicUIC+6fFoW+vFFsHCNFOrW+SjFfP+DA8PT9liMjw8nPuxrHP1aotgXuiG1p/paMVdygMDA0z1/iWxcePGXI9lncstZ7X58qcJvVppa83p1RbBovR1UnEzrEH1fyIDAwOuY5mGvk4qboY1mPqfC8CGDRv6asiCvPR1nYpZycTExKb5nwYGBtiwYcNW2/RrHUuzdSpOKmYVXIG/JVfUms2QK/BnxknFrIIr8GfGScWsgivwZ8Z1KmZWk+tUzKxQTipmlisnFTPLVc8kFd+6btYZipyiIzeensKscxRWUpG0s6TrJP1U0oOS3j7dfXl6CrPOUWRJ5RLg3yLiBEmvIptQbFp867pZ5yikpCLpNcC7gUsBIuKliHh6uvtzt2qzzlHU5c/ewDrgf0v6saSvStp+ujtzt2qzzlFUUpkNvAX4nxFxCPACcG7lRo1Oe+pu1Wado6ik8gjwSGTz/0A2B9BbKjeKiPGIWBARC+bNm1dzh6Ojo6xcuZKNGzeycuVKJ5Q63ARvrVLUZGKPAg9L2jcteg/wQBGx1NKrX7xSE/zk5KRHNrP8NTP0fp4P4GBgKXAf2QTtu9TavhVTdNTSy9M2eGoSawaeoiMfvTxtg0c2s2b4LuWc9HLfFzfBWys5qVTRy188N8FbKzmpVNHLXzw3wVsruU6lhvJpG4aGhhgbG/MXz/qOp+gws1y5otbMCuWkYma5clIxs1w5qZhZrpxUzCxXTipmlisnFTPLlZOKmeXKScXMcuWkYma5clIxs1w5qZhZrpxUzCxXhSYVSbPSvD83FBmHmeWn6JLKmcCDBcdgZjkqcoL2+cDvA18tKgYzy1+RJZUvA+cAVYdvb3SGQjPrHEVN0H4csDYiltXaLpqYodDMOkNRJZV3AMdLWglcAxwlaXFBsZhZjoqa9vS8iJgfESPAScAtEXFKEbGYWb6Kbv0xsx4zu+gAIuJW4NaCwzCznLikYma5clIxs1w5qZhZrpxUzCxXTipmlisnFTPLlZOKmeXKScXMcuWkYma5clIxs1w5qRRgYmKCkZERBgYGGBkZYWJiouiQzHJT+L0//WZiYoKFCxeyfv16ACYnJ1m4cCEAo6OjRYZmlguXVNps0aJFmxJKyfr161m0aFFBEZnly0mlzVatWtXUcrNu46TSZkNDQ00tN+s2TiptNjY2xuDg4BbLBgcHGRsbKygis3w5qbTZ6Ogo4+PjDA8PI4nh4WHGx8ddSWs9QxFRdAwNWbBgQSxdurToMMz6jqRlEbGg0e2LmqJjT0n/LulBSfdLOrOIOMwsf0X1U3kF+GRELJe0I7BM0pKIeKCgeMwsJ0VN0bEmIpan58+Rzae8RxGxmFm+Cq+olTQCHALcNcU6T3tq1mUKTSqSdgD+FTgrIp6tXO9pT826T2GtP5K2AW4AboyIixrYfh0w2fLA6psLPF50ENPk2IvR7bFvHxEN/1cvJKlIEnAF8GREnNX2AGZA0tJmmtc6iWMvRr/FXuQE7aeSTcx+T3q8v6BYzCxHhTQpR8R/ACri2GbWWoW3/nSh8aIDmAHHXoy+ir1ruumbWXdwScXMcuWkYma5clKpQdJlktZKWlG2bI6kJZIeSj93KTLGaqrdtNkN8UvaVtLdku5NsV+Qlu8l6a4U+9clvaroWKciaZakH0u6Ib3uirgBJK2U9JPUIrs0LWvqnHFSqe1y4JiKZecCN0fEbwE3p9edqHTT5n7A4cBfSNqf7oj/ReCoiDgIOBg4RtLhwOeBi1PsTwGnFxhjLWeS3c9W0i1xl/xuRBxc1j+luXMmIvyo8QBGgBVlr38G7J6e7w78rOgYG3wf3waO7rb4gUFgOfA2sl6ps9Pyt5P1xi48xop456cv3lFkPcbVDXGXxb8SmFuxrKlzxiWV5r02ItZAdrc1sFvB8dRVcdNmV8SfLiHuAdYCS4BfAE9HxCtpk0fozDvbvwycA2xMr3elO+IuCeAmScskLUzLmjpnPO9Pj6u8aTO7Q6LzRcQG4GBJOwPXA/tNtVl7o6pN0nHA2ohYJunI0uIpNu2ouCu8IyJWS9oNWCLpp83uwCWV5j0maXeA9HNtwfFUlW7a/FdgIiK+mRZ3TfwAEfE0cCtZvdDOkkr/COcDq4uKq4p3AMdLWglcQ3YJ9GU6P+5NImJ1+rmWLJm/lSbPGSeV5n0HOC09P42srqLjpJs2LwUejC3vAu/4+CXNSyUUJG0HvJes4vPfgRPSZh0Xe0ScFxHzI2IEOAm4JSJG6fC4SyRtn0ZiRNL2wPuAFTR7zhRdMdTJD+BqYA3wMtm18Olk18g3Aw+ln3OKjrNK7O8kK2bfB9yTHu/vhviBA4Efp9hXAOen5XsDdwP/CVwLvLroWGu8hyOBG7op7hTnvelxP7AoLW/qnHE3fTPLlS9/zCxXTipmlisnFTPLlZOKmeXKScXMcuWkYma5clKxrUjakG59XyHpWkmDOezz0w0c7/VTrPuwpH9Kz8+WtKr02jqTk4pN5deR3fp+APAS8LF6v6BMrfOpalIpO17N7usRcTFwfr1YrFhOKlbP94F9plohaSQNAvXPZMMT7Cnp5DTIzwpJn0/bfQ7YLpVGJuodUNJHJP1c0m1k99NYF3FSsarSTXDHAj+psdm+wJURcQjZ7QyfJ7uR7mDgMEkfjIhz2VwaGa1zzN2BC8iSydHA/jN/J9ZOTio2le3SWCZLgVVkNyZWMxkRd6bnhwG3RsS6yMYPmQDe3eSx31a2j5eArzf5+1Ywj6diU/l1RBzc4LYvlD3Pa7AW35DWxVxSsTzdBfyOpLmSZgEnA7eldS+n8V0a2ceRknZN25/YolitRVxSsdxExBpJ55GNHyLguxFRGntjHLhP0vJa9SppH58F7iAbdmI5MKu1kVuePPSBFU7S8xGxQ4PbfhhYEBF/2dqobLp8+WOd4Nlqnd/KSTobOA94tj1h2XS4pGJ1SSqN/FXpPRHxRLvjsc7mpGJmufLlj5nlyknFzHLlpGJmuXJSMbNc/X8FBCqzCAuRrQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 288x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Plot cycle period versus rotation period\n",
    "plt.figure(figsize=(4,4))\n",
    "plt.plot(cycrots['Prot'], cycrots['Pcyc'], 'ko')\n",
    "plt.xlabel('P_rot [d]')\n",
    "plt.ylabel('P_cyc [yr]')\n",
    "plt.title(\"Cycle Period vs Rotation Period\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}