From 10c9a2dc8ef5a6fd952e6ea08f9ab0ec35476d0a Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Wed, 8 Jun 2011 13:37:12 +0000
Subject: [PATCH] Re-added debug output; Fixed incorrect density values (wrong
sign)
---
src/build_system_wsmp.f90 | 24 ++++++++++++------------
1 file changed, 12 insertions(+), 12 deletions(-)
diff --git a/src/build_system_wsmp.f90 b/src/build_system_wsmp.f90
index 6996c279..404f271f 100644
--- a/src/build_system_wsmp.f90
+++ b/src/build_system_wsmp.f90
@@ -194,19 +194,19 @@ do ileaves=1,osolve%nleaves
iloc=ifind_loc_wsmp(i1loc,i1,ja,ia,n_iproc,nz_loc)
elem_dz=(osolve%z(osolve%icon(5,ileaves))-osolve%z(osolve%icon(1,ileaves)))*params%vex
if (k1 < 5) then
- !write (*,*) '****************'
- !write (*,*) 'k1: ',k1
- !write (*,*) 'avals before for proc ',iproc,': ',avals(iloc)
- !avals(iloc)=avals(iloc)-(params%dt*weightel(ileaves))/(2.d0*elem_dz)
- !write (*,*) 'avals after for proc ',iproc,': ',avals(iloc)
- !write (*,*) '****************'
+ write (*,*) '****************'
+ write (*,*) 'k1: ',k1
+ write (*,*) 'avals before for proc ',iproc,': ',avals(iloc)
+ avals(iloc)=avals(iloc)-(params%dt*(-weightel(ileaves)))/(2.d0*elem_dz)
+ write (*,*) 'avals after for proc ',iproc,': ',avals(iloc)
+ write (*,*) '****************'
else
- !write (*,*) '****************'
- !write (*,*) 'k1: ',k1
- !write (*,*) 'avals before for proc ',iproc,': ',avals(iloc)
- !avals(iloc)=avals(iloc)+(params%dt*weightel(ileaves))/(2.d0*elem_dz)
- !write (*,*) 'avals after for proc ',iproc,': ',avals(iloc)
- !write (*,*) '****************'
+ write (*,*) '****************'
+ write (*,*) 'k1: ',k1
+ write (*,*) 'avals before for proc ',iproc,': ',avals(iloc)
+ avals(iloc)=avals(iloc)+(params%dt*(-weightel(ileaves)))/(2.d0*elem_dz)
+ write (*,*) 'avals after for proc ',iproc,': ',avals(iloc)
+ write (*,*) '****************'
endif
endif
endif
--
GitLab