We developed a new way to model how planetary surfaces rise and sink as the deep interior slowly flows. Existing approaches are either costly or unstable. Our method represents the surface smoothly within a fixed grid, which avoids artificial air layers and numerical problems. Tests show it matches established results while running faster and working in more realistic settings, such as loaded surfaces and global models. This makes simulations of surface evolution more reliable and accessible.