diff --git a/src/build_system_wsmp.f90 b/src/build_system_wsmp.f90
index a4820636d9ac5de53cb24fe38f8b5f783a0bb6d3..f06305adf8beaa76fd862864f40565650c24291d 100644
--- a/src/build_system_wsmp.f90
+++ b/src/build_system_wsmp.f90
@@ -215,9 +215,29 @@ if (ndof==3) then
    call mpi_allreduce (is_plastic,osolve%is_plastic,osolve%nleaves,mpi_logical,mpi_lor,mpi_comm_world,ierr)
    call mpi_allreduce (vbounded,vbounded2,osolve%nleaves,mpi_logical,mpi_lor,mpi_comm_world,ierr)
    if (iproc.eq.0) write (*,'(a,i10)') shift//'viscosity capped in', count(vbounded2)
-   call mpi_allreduce (weightel,weightel_glob,osolve%nleaves,mpi_double_precision,mpi_min,mpi_comm_world,ierr)
+   !call mpi_allreduce (weightel,weightel_glob,osolve%nleaves,mpi_double_precision,mpi_min,mpi_comm_world,ierr)
+   !print *,'ndof: ',ndof
+   !print *,'weightel_glob min: ',minval(weightel_glob)
+   !print *,'weightel_glob max: ',maxval(weightel_glob)
+   !print *,'weightel min: ',minval(weightel)
+   !print *,'weightel max: ',maxval(weightel)
 end if
 
+!weightel_glob=0.d0
+!print *,'1. ndof: ',ndof
+!print *,'1. weightel_glob min: ',minval(weightel_glob)
+!print *,'1. weightel_glob max: ',maxval(weightel_glob)
+!print *,'1. weightel min: ',minval(weightel)
+!print *,'1. weightel max: ',maxval(weightel)
+
+call mpi_allreduce (weightel,weightel_glob,osolve%nleaves,mpi_double_precision,mpi_min,mpi_comm_world,ierr)
+
+!print *,'2. ndof: ',ndof
+!print *,'2. weightel_glob min: ',minval(weightel_glob)
+!print *,'2. weightel_glob max: ',maxval(weightel_glob)
+!print *,'2. weightel min: ',minval(weightel)
+!print *,'2. weightel max: ',maxval(weightel)
+
 call heap (threadinfo,'is_plastic','build_system',size(is_plastic),'bool',-1)
 deallocate (is_plastic)                              
 call heap (threadinfo,'vbounded','build_system',size(vbounded),'bool',-1)