Re-added debug output; Fixed incorrect density values (wrong sign)

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