nloading
                f = ku * (ovlp - res);
            else
                f = 0;
            end
            
            % Total tangential force vector (against deformation and motion)
            this.total_force = -f * dir;
        end
    end
end
##### SOURCE END #####
--></body></html>