From db4fc4af3c4a15fbd8c6fa62e5aa39aa2ad31818 Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Wed, 19 May 2010 18:24:17 +0000
Subject: [PATCH] Test version with new isostasy routine and temperature calc
---
src/build_system_wsmp.f90 | 22 +++++++++++++++++++++-
1 file changed, 21 insertions(+), 1 deletion(-)
diff --git a/src/build_system_wsmp.f90 b/src/build_system_wsmp.f90
index a4820636..f06305ad 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)
--
GitLab