Add procedure to abort if we're in an omp parallel region

Some set up routines are known not to be thread-safe.  If called in a
parallel region, race conditions may lead to invalid code.  To guard
against accidental usage in a threaded region, add an abort to known
thread-unsafe functions.  The new subroutine
abort_if_in_parallel_region wraps this check and should be called at
the beginning of thread-unsafe procedures.