En este artículo, se presenta un proceso h-adaptativo que define una malla óptima para calcular la solución de problemas de contorno parabólicos usando el método de Diferencias Miméticas. La estimación del error, en la variable espacial, se hace a partir de la versión discreta del operador gradiente. La experimentación numérica evidencia el buen comportamiento del procedimiento.