def compute_primary_air_flow(Firepower):
    """Apply agenbroad method for estimating the intake flow rate
    Args:
        Firepower (double): user-defined firepower for the analysis
    Returns:
        m_dot_primary_air (double): mass flow rate of naturally entrained air (kg/s).
        mol_dot_primary_air (double): molar flow rate of naturally entrained air (mol/s).
    """
    cp_average = 1.044# kJ/kg-K
    Q_dot = Firepower # kJ/s (kW)
    T_amb = 300 # K
    T_h = 1274 # K

    m_dot_primary_air = (Q_dot)/(cp_average*(T_h - T_amb)) #kg/s

    mw_air = 0.018 #kg/mol

    mol_dot_primary_air = m_dot_primary_air/mw_air

    print("mass flow rate of air:")
    print(m_dot_primary_air)

    print("molar flow rate of primary air:")
    print(mol_dot_primary_air)

    return mol_dot_primary_air, m_dot_primary_air