Hertz subsurface stress is now properly computed for r=0 (i.e. on the center line axis).