diff --git a/src/build_system_wsmp.f90 b/src/build_system_wsmp.f90
index 6996c279f3b143990883caa83efc0053f1d383a3..404f271f5ddc9f5e633bf33edde52a1b77959187 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