<div class="gmail_quote">On Sat, Nov 5, 2011 at 14:34, Dharmendar Reddy <span dir="ltr">&lt;<a href="mailto:dharmareddy84@gmail.com">dharmareddy84@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
 I am writing a 3D FEM code for learning purpose (experimenting with object oriented concepts in Fortran 2003). Can some one tell me (pseudo code) how to implement a non homogenous Neumann boundary condition. You can also point me to a book. I am using tetrahedral elements for solving  Poisson equation. I am using a 4 point Gauss quadrature, and linear basis functions for cell volume integrals. I am confused on how to do the integration on the facet. </blockquote>
</div><br><div>You need to compute an integral over the face. It appears in the weak form. Any book on finite element methods should cover this.</div>