CPML is now applied to update of the surface function.
This could be necessary because these functions overwrite the horizontal components sxx and szz on the surface.
Tests with sofi2D showd that if CPML is not applied again for the horizontal surface points, it could lead to
instabilities.