Added debug output

src/build_system_wsmp.f90

@@ -193,10 +193,20 @@ do ileaves=1,osolve%nleaves
                     if (kfix((osolve%icon(k1,ileaves)-1)*ndof+idof1) < 1) then
                       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                        
+                      if (k1 < 5) then
+                        write (*,*) '****************'
+                        write (*,*) 'k1: ',k1
+                        write (*,*) 'avals before for proc ',iproc,': ',avals
                         avals(iloc)=avals(iloc)-(params%dt*weightel(ileaves))/(2.d0*elem_dz)
+                        write (*,*) 'avals after for proc ',iproc,': ',avals
+                        write (*,*) '****************'
                       else
+                        write (*,*) '****************'
+                        write (*,*) 'k1: ',k1
+                        write (*,*) 'avals before for proc ',iproc,': ',avals
                         avals(iloc)=avals(iloc)+(params%dt*weightel(ileaves))/(2.d0*elem_dz)
+                        write (*,*) 'avals after for proc ',iproc,': ',avals
+                        write (*,*) '****************'
                       endif
                     endif
                   endif