Created a function that branches calls to DOPAIR1 so that correct version is called. This decision is based upon whether MPI is running, vectorisation and orientation of the neighbouring cells.