def globalOutgassing(Fmod_out, Q, m):
    """
    This function will calculate the outgassing flux (equation S9).

    Inputs:
        Fmod_out - the modern Earth's outgassing rate [mol C yr-1]
        Q        - pore space heat flow relative to modern Earth [dimensionless]
        m        - scaling parameter [dimensionless]

    Returns:
        F_out - the global outgassing flux [mol C yr-1]
    """

    F_out = Fmod_out*Q**m
    
    return F_out