We consider the problem of the generation of adaptive grids, where, for a given domain Ω, we have to compute a partition that depends on a given map f providing the desired local mesh size. We propose a numerical method for the solution of this problem in the case of planar quadrilateral grids....