Switched to using mpi_init_thread to initialize MPI with thread support

src/DOUAR.f90

@@ -25,7 +25,7 @@ include 'mpif.h'
 integer i,j,is,iter,istep,niter,nedge,nedgen,ileaves,matel
 integer ntriangle,ndof,ndoft,material0,nsurfacen
 integer numarg,iargc
-integer ierr,nproc,iproc,k,nrem
+integer ierr,nproc,iproc,k,nrem,provided
 integer current_level
 integer ie_ov,ie_loc,ie_level,iter_nl
 integer err,iunit,cltemp
@@ -109,7 +109,9 @@ nleaves_old_mem2=0
 !------------------------------------------------------------------------------|
 !------------------------------------------------------------------------------|
 
-call mpi_init(ierr)
+provided=0
+!call mpi_init(ierr)
+call mpi_init_thread(MPI_THREAD_MULTIPLE,provided,ierr)
 call mpi_comm_size (mpi_comm_world,nproc,ierr)
 call mpi_comm_rank (mpi_comm_world,iproc,ierr)