Sublimating surfaces are out of equilibrium. It has been proposed that sublimation can be compensated by an impinging atomic flux to obtain equilibrium. This work concerns the effect of such an impinging flux on the stability of surfaces in various situations. For this purpose we combine Kinetic Monte Carlo Simulations with analytical developments based on the Burton-Cabrera-Frank (BCF) classical theory. We show that a perfect compensation of the sublimation is possible for vicinal surfaces but not when 2D islands are present on a surface. We thus study the effect of an impinging flux on the dynamic of a 2D island on a surface. We show that the 2D island area generally varies with time t as −tα. In absence of any impinging flux the value of the exponent α enables to identify the main mechanism at work (diffusion limited or attachment-detachment limited). On the contrary, in presence of an impinging flux the value of the exponent α is not enough to identify the main mechanism limiting the area change.