Index: /issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.cpp	(revision 17478)
+++ /issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.cpp	(revision 17479)
@@ -1199,5 +1199,5 @@
 }/*}}}*/
 /*FUNCTION GaussLegendreTetra{{{*/
-void GaussLegendreTetra( int* pngaus, IssmPDouble** pl1, IssmPDouble** pl2, IssmPDouble** pl3, IssmPDouble** pl4, IssmPDouble** pwgt, int iord ) {
+void GaussLegendreTetra( int* pngaus, IssmDouble** pl1, IssmDouble** pl2, IssmDouble** pl3, IssmDouble** pl4, IssmDouble** pwgt, int iord ) {
 	/* Gauss quadrature points for the tetrahedron.
 
@@ -1402,9 +1402,9 @@
 		*pngaus=np[iord-1];
 
-		*pl1  =xNew<IssmPDouble>(*pngaus);
-		*pl2  =xNew<IssmPDouble>(*pngaus);
-		*pl3  =xNew<IssmPDouble>(*pngaus);
-		*pl4  =xNew<IssmPDouble>(*pngaus);
-		*pwgt =xNew<IssmPDouble>(*pngaus);
+		*pl1  =xNew<IssmDouble>(*pngaus);
+		*pl2  =xNew<IssmDouble>(*pngaus);
+		*pl3  =xNew<IssmDouble>(*pngaus);
+		*pl4  =xNew<IssmDouble>(*pngaus);
+		*pwgt =xNew<IssmDouble>(*pngaus);
 
 		for (i=0; i<*pngaus; i++) {
@@ -1422,9 +1422,9 @@
 		*pngaus=nigaus*nigaus*nigaus;
 
-		*pl1  =xNew<IssmPDouble>(*pngaus);
-		*pl2  =xNew<IssmPDouble>(*pngaus);
-		*pl3  =xNew<IssmPDouble>(*pngaus);
-		*pl4  =xNew<IssmPDouble>(*pngaus);
-		*pwgt =xNew<IssmPDouble>(*pngaus);
+		*pl1  =xNew<IssmDouble>(*pngaus);
+		*pl2  =xNew<IssmDouble>(*pngaus);
+		*pl3  =xNew<IssmDouble>(*pngaus);
+		*pl4  =xNew<IssmDouble>(*pngaus);
+		*pwgt =xNew<IssmDouble>(*pngaus);
 
 		/*  get the gauss points in each direction  */
Index: /issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.h
===================================================================
--- /issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.h	(revision 17478)
+++ /issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.h	(revision 17479)
@@ -12,5 +12,5 @@
 void GaussLegendreTria(int* pngaus, IssmDouble** pl1, IssmDouble** pl2, IssmDouble** pl3, IssmDouble** pwgt, int iord);
 #define MAX_TETRA_SYM_ORD    6
-void GaussLegendreTetra(int* pngaus, IssmPDouble** pl1, IssmPDouble** pl2, IssmPDouble** pl3, IssmPDouble** pl4, IssmPDouble** pwgt, int iord);
+void GaussLegendreTetra(int* pngaus, IssmDouble** pl1, IssmDouble** pl2, IssmDouble** pl3, IssmDouble** pl4, IssmDouble** pwgt, int iord);
 #define MAX_LINE_GLOB_PTS    5
 void GaussLobatto(IssmPDouble** pxgaus, IssmPDouble** pxwgt, int ngaus);