Skip to content

Instantly share code, notes, and snippets.

@RuthAngus
Created April 7, 2019 12:31
Show Gist options
  • Select an option

  • Save RuthAngus/73b82c3813ee6c446bd800e079dcd635 to your computer and use it in GitHub Desktop.

Select an option

Save RuthAngus/73b82c3813ee6c446bd800e079dcd635 to your computer and use it in GitHub Desktop.
The rotation period distribution of hot stars.
import pandas as pd
import matplotlib.pyplot as plt
import numpy as np
url = "https://arxiv.org/src/1402.5694v2/anc/Table_1_Periodic.txt"
mc = pd.read_csv(url)
plt.hist(np.log10(mc.Prot.values[mc.Teff.values > 6250]), 50)
plt.xlabel("log10(Rotation period [days])")
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment