Added a program to calculate properties at the surface.

This program is exceedingly inefficient and slow, but it does
determine the surface velocity of the slab and the mass flux at
the surface.