def carbon_abundance(Z):
    """
    Set the carbon abundance. 
    We assume Asplund et al 2009 abundance at Zsun and that Ac scales linearly with Z. Z in solar units
    """
    Ac = 2.7e-4 
    return Ac*Z