def data_for_scipy(rad=10., vel=21000, den=3000., stre=1e5, theta=45., dt=0.05):
    """
    Create the data necessary for the scipy solver

    Parameters
    ----------

    rad : float, optional
        The radius of the asteroid in meters. Default value is 10 meters.

    vel : float, optional
        The entery speed of the asteroid in meters/second.
        Default value is 21000 meters/second.

    den : float, optional
        The density of the asteroid in kg/m^3. Default value is 3000 kg/m^3.

    stre : float, optional
        The strength of the asteroid (i.e., the ram pressure above which
        fragmentation and spreading occurs) in N/m^2 (Pa). Default value is 10^5 N/m^2.

    theta : float, optional
        The initial trajectory angle of the asteroid to the horizontal
        By default, input is in degrees. Default value is 45 degrees.

    dt : float, optional
        The output timestep, in s. Default value is 0.05 seconds.

    Returns
    -------

    dict
        A dictionary with all the values needed for the scipy solution,
        including initial conditions and constants
    """
    input_data = {'radius': rad,
                  'velocity': vel,
                  'density': den,
                  'strength': stre,
                  'angle': theta,
                  'init_altitude':100000,
                  'dt': dt,
                  'radians': False}

    #data within the Planet class, can be customized
    setting = {'atmos_func':'constant', 'atmos_filename':None,
               'Cd':1., 'Ch':0.1, 'Q':1e7, 'Cl':1e-3,
               'alpha':0.3, 'Rp':6371e3,
               'g':9.81, 'H':8000., 'rho0':1.2}

    #integrate
    alldata = {**setting, **input_data}

    return alldata