dummy functions for shear stress and q
nx and ny don't really need to be of type double
mac() now accepts an additional PointGrille* argument
mac() is now of type void
added some boilerplate and console spam on stdout