def draw_orbits(number):
    """ Semi-major axis is fixed at 100 au and long. of asc. node is fixed at 0 deg.
    Written by Logan Pearce, 2019
    """
    import astropy.units as u
    import numpy as np
    sma = 100.0 * u.au
    sma = np.array(np.linspace(sma, sma, number))
    ecc = np.random.uniform(0.0, 1.0, number)
    cosi = np.random.uniform(-1.0, 1.0, number)
    inc = np.degrees(np.arccos(cosi))
    argp = np.random.uniform(0.0, 360.0, number)
    lon = np.degrees(0.0)
    lon = np.array([lon] * number)
    orbit_fraction = np.random.uniform(0.0, 1.0, number)
    return sma, ecc, inc, argp, lon, orbit_fraction
