def qO2_Calc(qH2, rho):
    """
    Calculate qO2.

    :param qH2: molar flow of hydrogen [kmol.s^(-1)]
    :type qH2 : float
    :param rho: hydrogen-oxygen flow rate
    :type rho : float
    :return: qO2 [kmol.s^(-1)] as float
    """
    try:
        result = (qH2 / rho)
        return result
    except (TypeError, ZeroDivisionError):
        print(
            "[Error] qO2 Calculation Error (qH2:%s, rho:%s)" %
            (str(qH2), str(rho)))