We explore the dynamics of cosmological models with two coupled dark components with energy densities ρA and ρB and constant equation of state (EoS) parameters wA and wB. We assume that the coupling is of the form Q = H q(ρA, ρB), so that the dynamics of the two components turns out to be scale independent, i.e. does not depend explicitly on the Hubble scalar H. With this assumption, we focus on the general linear coupling q = qo + qA ρA + qB ρB, which may be seen as arising from any q(ρA, ρB) at late time and leads in general to an effective cosmological constant. In the second part of the paper we consider observational constraints on the form of the coupling from SN Ia data, assuming that one of the components is cold dark matter (CDM), i.e. wB = 0, while for the other the EoS parameter can either have a standard (wA > −1) or phantom (wA < −1) value. We find that the constant part of the coupling function is unconstrained by SN Ia data and, among typical linear coupling functions, the one proportional to the dark energy density ρA is preferred in the strong coupling regime, |qA| > 1. Models with phantom wA favor a positive coupling function, increasing ρA. In models with standard wA, not only a negative coupling function is allowed, transferring energy to CDM, but the uncoupled sub-case falls at the border of the likelihood.PACS numbers: 98.80.Jk; 95.35.+d; 95.36.+x