A numerical algorithm to compute steady-state ground level concentration from elevated sources by means of a K-model which takes into account the spatial variability of wind and diffusivity and neglects horizontal diffusion is discussed. The boundary value problem to be treated, also for a point source, is always reduced to a two dimensional one and it is solved on an optimized grid. In this way the proposed method is made computationally comparable with the classical Gaussian plume model.