Made a timed version of region_mod, so that all subcalls of interest can be timed.
For some reason, I had not seen that only half of the routines had been parallelized. Now all of them are.
Merge branch 'performance'
Add output features
Add new utility program