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