def mdot(snu, nu=10, Te=1e4, vinf=2000, mue=1.3, d=1):
    """
    [snu] = mJy
    nu = 10 GHz (frequency)
    Te = 10^4 K (ionized gas temp)
    vinf = 2000 km/s (maximum wind speed)
    mue = 1.3 (mean atomic weight per free electron)
    d = distance (kpc)

    return: mdot in Msun/yr
    """
    mdot = ( snu / 7.26 * (nu/10.)**-0.6 * (Te/1e4)**-0.1 * d**2 * (mue*vinf/100)**(4/3.) ) **(3/4.) * 1e-6
    return mdot