Funkcja ta umożliwia automatyczne zagęszczenie siatki ES na powierzchniach. Zagęszczenie siatki jest stopniowe. Na każdym kroku siatka ES jest tworzona na podstawie porównania błędów wyników w poprzednim kroku obliczeń. Błąd numeryczny jest obliczany na podstawie wyników dla elementów powierzchniowych i jest oparty na sformułowaniu energetycznym Zienkiewicza-Zhu.
Ocena błędu jest przeprowadzana dla liniowej analizy statycznej. Wybieramy przypadek obciążenia (lub kombinację obciążeń), dla którego wygenerowana jest siatka ES. Siatka ES jest następnie wykorzystywana do wszystkich obliczeń.