@@ -0,0 +1,1633 @@
+from __future__ import division%0Afrom __future__ import print_function%0Aimport numpy as np%0A'''%0AUse delta-nu, nu-max, and temp to estimate oscillation amplitudes.%0ABased on Corsaro et al. 2013.%0A'''%0A# Read in text file containing useful info%0Ainfile = '../../RGEB_seismicinfo.txt'%0A%0A# Define columns%0AKICcol = 0%0Anumaxcol = 1; numaxerrcol = 2%0Adnucol = 3; dnuerrcol = 4%0Ateffcol = 5; tefferrcol = 6%0Aampcol = 7%0A%0A# Read in data%0AKICs = np.loadtxt(infile, comments='#', usecols=(KICcol,), unpack=True)%0Anumaxes, numaxerrs = np.loadtxt(infile, comments='#', usecols=(numaxcol,numaxerrcol), unpack=True)%0Adnus, dnuerrs = np.loadtxt(infile, comments='#', usecols=(dnucol,dnuerrcol), unpack=True)%0Ateffs, tefferrs = np.loadtxt(infile, comments='#', usecols=(teffcol,tefferrcol), unpack=True)%0Aamps = np.loadtxt(infile, comments='#', usecols=(ampcol,), unpack=True)%0A%0A#for star in KICs:%0A#    print(star)%0A%0Adef ampguess(numax = 3100, dnu = 135, teff = 5777):%0A    '''%0A    Calculate amplitude given numax, dnu, and teff%0A    '''%0A    # Solar values for reference%0A    amp_s = 3.6%0A    numax_s = 3100%0A    dnu_s = 135%0A    teff_s = 5777%0A    # Constants from Corsaro et al. 2013 paper (M_%7B4,beta%7D specifically)%0A    s = 0.602; s_err = 0.008; r = 5.87; r_err = 0.14; t = 1.31; t_err = 0.02%0A    ampfrac = (numax/numax_s)**(2*s-3*t) * (dnu/dnu_s)**(4*t-4*s) * (teff/teff_s)**(5*s-1.5*t-r+0.2)%0A    amp = ampfrac*amp_s%0A    return amp%0A    %0A%0Afor star, numax, dnu, teff, amp in zip(KICs, numaxes, dnus, teffs, amps):%0A    ampcalc = ampguess(numax, dnu, teff)%0A    print('%7B0%7D: %7B1:.1f%7D predicted, %7B2:.1f%7D observed (%7B3:.0f%7D percent)'.format(int(star), ampcalc, amp, amp/ampcalc*100))%0A

use Corsaro et al 2013 to naively predict oscillation amplitude