Test version with new isostasy routine and temperature calc

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)