示例#1
0
def test_ldc():
    filt_type = 6  #  V band
    temp = 6400.

    ldc = PyLDC(DB_NAME, filt_type)
    coeffs = ldc.coefficients(temp)

    assert type(coeffs) == np.ndarray
    assert coeffs.size == 4
示例#2
0
def get_catalogue_data(planet_name):
    results = exodb.queryForColumns([
        'planet.name', 'system.period', 'system.epoch',
        'system.a', 'system.i', 'star.radius', 'planet.radius',
        'star.mass', 'star.teff'],
        where='replace(planet.name, " ", "") = "%s"' %
            planet_name.replace(" ", ""))

    if len(results):
        results = results[0]

        m = PyModel()
        m.period = float(results[1])
        m.epoch = float(results[2])
        m.a = float(results[3])
        m.i = float(results[4])
        m.rs = float(results[5])
        m.rp = float(results[6])
        m.mstar = float(results[7])
        m.teff = float(results[8])

        ldc = PyLDC(
                os.path.join(
                    os.path.expanduser("~"),
                    "work",
                    "ReferenceDatabase",
                    "Reference.db"),
                6
                )

        coeffs = ldc.coefficients(m.teff)

        m.c1 = coeffs[0]
        m.c2 = coeffs[1]
        m.c3 = coeffs[2]
        m.c4 = coeffs[3]

        return m
    else:
        raise RuntimeError("Cannot create planet model")