Index: /issm/trunk/src/c/objects/Beam.cpp
===================================================================
--- /issm/trunk/src/c/objects/Beam.cpp	(revision 321)
+++ /issm/trunk/src/c/objects/Beam.cpp	(revision 322)
@@ -362,11 +362,8 @@
 		GetJacobianDeterminant(&Jdet, &z_list[0],gauss_coord);
 		
-		//Get nodal functions
-		GetNodalFunctions(l1l2, gauss_coord);
-
 		for(j=0;j<NDOF2;j++){
 			pe_g_gaussian[NDOF2+j]=constant_part*pow((s[0]-z_g)/B,n)*slope[j]*Jdet*gauss_weight;
 		}
-
+		
 		//add pe_gaussian vector to pe:
 		for(j=0;j<numdofs;j++){
@@ -387,4 +384,11 @@
 		pe_g[1]+=vb;
 	}
+
+	/*Add pe_g to global vector pg: */
+	VecSetValues(pg,numdofs,doflist,(const double*)pe_g,ADD_VALUES);
+
+	cleanup_and_return: 
+	xfree((void**)&segment_gauss_coord);
+	xfree((void**)&gauss_weights);
 }