def Fe(mw_atmos, mh=1):
    """Defines properties for Fe as condensible"""
    if mh != 1:
        raise Exception(
            'Alert: No M/H Dependence in Fe Routine. Consult your local theorist to determine next steps.'
            )
    gas_mw = 55.845
    gas_mmr = 0.0013
    rho_p = 7.875
    return gas_mw, gas_mmr, rho_p
