Index: /issm/trunk/src/c/objects/Gauss/GaussPenta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Gauss/GaussPenta.cpp	(revision 5674)
+++ /issm/trunk/src/c/objects/Gauss/GaussPenta.cpp	(revision 5675)
@@ -195,4 +195,22 @@
 }
 /*}}}*/
+/*FUNCTION GaussPenta::GaussFaceTria{{{1*/
+void GaussPenta::GaussFaceTria(int index1, int index2, int index3, int order){
+
+	/*in debugging mode: check that the default constructor has been called*/
+	ISSMASSERT(numgauss==-1);
+
+	/*Basal Tria*/
+	if(index1==0 && index2==1 && index3==2){
+		GaussLegendreTria(&numgauss,&coords1,&coords2,&coords3,&weights,order);
+		coords4=(double*)xmalloc(numgauss*sizeof(double));
+		for(int i=0;i<numgauss;i++) coords4[i]=-1.0;
+	}
+	else{
+		ISSMERROR("Tria not supported yet");
+	}
+
+}
+/*}}}*/
 /*FUNCTION GaussPenta::begin{{{1*/
 int GaussPenta::begin(void){
Index: /issm/trunk/src/c/objects/Gauss/GaussPenta.h
===================================================================
--- /issm/trunk/src/c/objects/Gauss/GaussPenta.h	(revision 5674)
+++ /issm/trunk/src/c/objects/Gauss/GaussPenta.h	(revision 5675)
@@ -42,4 +42,5 @@
 		void GaussPoint(int ig);
 		void GaussVertex(int iv);
+		void GaussFaceTria(int index1, int index2, int index3, int order);
 		void GaussCenter(void);
 		void SynchronizeGaussTria(GaussTria* gauss_tria);