label='$%s~\mathrm{Gyr}$~$\mathrm{M\&H~(2008)}$' %ages[i],
zorder=0)
xs, ys = iso_calc(pars3, ages[i])
xs -= .45
pl.plot(xs, ys, color='k', linestyle='-', linewidth=lw,
label = '$%s~\mathrm{Gyr}$~$\mathrm{Angus~\emph{et~al.}~(2015)}$'
% ages[i], zorder=0)
xs, ys1 = iso_calc(pars3-pars3_err, ages[i])
xs -= .45
xs, ys2 = iso_calc(pars3+pars3_err, ages[i])
xs -= .45
# pl.fill_between(xs, ys1, ys2, facecolor='0.5', alpha=0.3, edgecolor='None',
# zorder=0)
from prob_uncerts import sigmas
samples = sigmas("ACHF")
for j in range(np.shape(samples)[1]):
par_samp = np.zeros(4)
par_samp[:3] = samples[:, j]
par_samp[3] = 0.45
xs, ys = iso_calc(par_samp, ages[i])
xs -= 0.45
pl.plot(xs, ys, color="k", zorder=0, alpha=.04)
pl.xlabel("$\mathrm{B-V-}~c$")
pl.ylabel("$\mathrm{P_{rot} (days)}$")
pl.xlim(10**-3, 1.4)
pl.ylim(1, 10**2)
pl.legend(loc='upper left')
pl.loglog()
pl.savefig("/Users/angusr/Python/Gyro/gyro_paper/p_vs_bv%s.pdf" % i)
pl.clf()
sun = a==4.568
bv -= 0.45
print bv_err[l], "\n"
print t_err[l]
print logg_err[l]
print p_err[l], "\n"
pl.errorbar(bv[l], p[l], xerr=bv_err[l], yerr=p_err[l], color='k', \
fmt='o', mec='k', capsize=0, markersize=5, ecolor='.6')
pl.errorbar(bv[l2], p[l2], xerr=bv_err[l2], yerr=p_err[l2], color='k', \
fmt='o', mec='.7', capsize=0, markersize=5, ecolor='.6')
pl.errorbar(bv[sun], p[sun], xerr=bv_err[sun], yerr=p_err[sun], color='r', \
fmt='o', mec='r', capsize=0, markersize=10, ecolor='.6')
from prob_uncerts import sigmas
samples = sigmas(fname)
for i in range(np.shape(samples)[1]):
par_samp = np.zeros(4)
par_samp[:3] = samples[:, i]
par_samp[3] = 0.45
xs, ys = iso_calc(par_samp, age)
xs -= 0.45
pl.plot(xs, ys, color=c, zorder=0, alpha=.1)
# Add Isochrones
xs, ys = iso_calc(pars, age)
xs -= 0.4
c = "k"
pl.plot(xs, ys, color=c, linestyle='-.', linewidth=lw, \
label='$%s~\mathrm{Gyr}$~$\mathrm{(Barnes~2007)}$' %age, zorder=0)
xs, ys = iso_calc(pars2, age)
