Set up communications
subroutine comms_setup
!! Set up communications
implicit none
#ifdef MPI
integer :: ierr
call mpi_init(ierr)
if (ierr .ne. 0) call io_error('MPI initialisation error')
#endif
call comms_setup_vars
end subroutine comms_setup