We present an improved adaptive mesh process that allows the accurate control of the numerical solution of interest derived from the solution of the partial differential equation. In the cases of two-dimensional studies, such an adaptive meshing is applied to compute phenomenon involving high field gradients in near-field (electric intensity, Poynting's vector, optical forces,...). We show, that this improved scheme permits to decrease drastically the computational time and the memory requirements.