Hello, my 2 cents : apparently, the boundary conditions are not applied as you want (from the file edge_elts.pre, I see that zero is assigned to all the edge elements on the boundary). Now, I don't know how to apply these BC... (and I'm interested) O.C.