Index: /issm/trunk/src/c/Container/Inputs.cpp
===================================================================
--- /issm/trunk/src/c/Container/Inputs.cpp	(revision 5659)
+++ /issm/trunk/src/c/Container/Inputs.cpp	(revision 5660)
@@ -209,60 +209,4 @@
 }
 /*}}}*/
-/*FUNCTION Inputs::GetParameterValues(double* values,double* gauss_pointers, int numgauss,int enum_type){{{1*/
-void Inputs::GetParameterValues(double* values,double* gauss_pointers, int numgauss,int enum_type){
-
-	vector<Object*>::iterator object;
-	Input* input=NULL;
-	bool   found=false;
-
-	/*Go through inputs and check whether any input with the same name is already in: */
-	for ( object=objects.begin() ; object < objects.end(); object++ ){
-
-		input=(Input*)(*object); 
-		if (input->EnumType()==enum_type){
-			found=true;
-			break;
-		}
-	}
-
-	if (!found){
-		/*we could not find an input with the correct enum type. No defaults values were provided, 
-		 * error out: */
-		ISSMERROR("could not find input with enum type %i (%s)",enum_type,EnumToString(enum_type));
-	}
-
-	/*Ok, we have an input if we made it here, request the input to return the values: */
-	input->GetParameterValues(values,gauss_pointers,numgauss);
-
-}
-/*}}}*/
-/*FUNCTION Inputs::GetParameterValues(double* values,double* gauss_pointers, int numgauss,int enum_type,double* defaultvalues){{{1*/
-void Inputs::GetParameterValues(double* values,double* gauss_pointers, int numgauss,int enum_type,double* defaultvalues){
-
-	int i;
-	vector<Object*>::iterator object;
-	Input* input=NULL;
-	bool   found=false;
-
-	/*Go through inputs and check whether any input with the same name is already in: */
-	for ( object=objects.begin() ; object < objects.end(); object++ ){
-
-		input=(Input*)(*object); 
-		if (input->EnumType()==enum_type){
-			found=true;
-			break;
-		}
-	}
-
-	if (!found){
-		/*we could not find an input with the correct enum type. Return the default values: */
-		for(i=0;i<numgauss;i++) values[i]=defaultvalues[i];
-	}
-	else{
-		input->GetParameterValues(values,gauss_pointers,numgauss);
-	}
-
-}
-/*}}}*/
 /*FUNCTION Inputs::GetParameterAverage{{{1*/
 void Inputs::GetParameterAverage(double* pvalue,int enum_type){
Index: /issm/trunk/src/c/Container/Inputs.h
===================================================================
--- /issm/trunk/src/c/Container/Inputs.h	(revision 5659)
+++ /issm/trunk/src/c/Container/Inputs.h	(revision 5660)
@@ -50,6 +50,4 @@
 		void GetParameterValue(double* pvalue,GaussTria* gauss,int enum_type);
 		void GetParameterValue(double* pvalue,double* gauss,int enum_type,double defaultvalue);
-		void GetParameterValues(double* values,double* gauss_pointers, int numgauss,int enum_type);
-		void GetParameterValues(double* values,double* gauss_pointers, int numgauss,int enum_type,double* defaultvalues);
 		void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, double* gauss,int enum_type);
 		/*}}}*/
Index: /issm/trunk/src/c/Makefile.am
===================================================================
--- /issm/trunk/src/c/Makefile.am	(revision 5659)
+++ /issm/trunk/src/c/Makefile.am	(revision 5660)
@@ -480,6 +480,4 @@
 					./modules/ComputeBasalStressx/ComputeBasalStressx.h\
 					./modules/ComputeBasalStressx/ComputeBasalStressx.cpp\
-					./modules/ComputePressurex/ComputePressurex.h\
-					./modules/ComputePressurex/ComputePressurex.cpp\
 					./modules/ComputeStrainRatex/ComputeStrainRatex.h\
 					./modules/ComputeStrainRatex/ComputeStrainRatex.cpp\
@@ -1034,6 +1032,4 @@
 					./modules/ComputeBasalStressx/ComputeBasalStressx.h\
 					./modules/ComputeBasalStressx/ComputeBasalStressx.cpp\
-					./modules/ComputePressurex/ComputePressurex.h\
-					./modules/ComputePressurex/ComputePressurex.cpp\
 					./modules/ComputeStrainRatex/ComputeStrainRatex.h\
 					./modules/ComputeStrainRatex/ComputeStrainRatex.cpp\
Index: /issm/trunk/src/c/modules/modules.h
===================================================================
--- /issm/trunk/src/c/modules/modules.h	(revision 5659)
+++ /issm/trunk/src/c/modules/modules.h	(revision 5660)
@@ -15,5 +15,4 @@
 #include "./Chacox/Chacox.h"
 #include "./ComputeBasalStressx/ComputeBasalStressx.h"
-#include "./ComputePressurex/ComputePressurex.h"
 #include "./ComputeStrainRatex/ComputeStrainRatex.h"
 #include "./ConfigureObjectsx/ConfigureObjectsx.h"
Index: /issm/trunk/src/c/objects/Elements/Element.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/Element.h	(revision 5659)
+++ /issm/trunk/src/c/objects/Elements/Element.h	(revision 5660)
@@ -51,5 +51,4 @@
 		virtual void   InputDepthAverageAtBase(int enum_type,int average_enum_type,int object_enum)=0;
 		virtual void   ComputeBasalStress(Vec sigma_b)=0;
-		virtual void   ComputePressure(Vec p_g)=0;
 		virtual void   ComputeStrainRate(Vec eps)=0;
 		virtual double MassFlux(double* segment,bool process_units)=0;
Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5659)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5660)
@@ -611,50 +611,4 @@
 	/*Add value to output*/
 	VecSetValue(sigma_b,id-1,(const double)value,INSERT_VALUES);
-}
-/*}}}*/
-/*FUNCTION Penta::ComputePressure {{{1*/
-void  Penta::ComputePressure(Vec pg){
-
-	int i;
-	const int numgrids=6;
-	int    doflist[numgrids];
-	double pressure[numgrids];
-	double rho_ice,g;
-	double surface[numgrids];
-	double xyz_list[numgrids][3];
-	double gauss[numgrids][4]={{1,0,0,-1},{0,1,0,-1},{0,0,1,-1},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
-
-	/*inputs: */
-	bool onwater;
-	Input* surface_input=NULL;
-
-	/*retrieve inputs :*/
-	inputs->GetParameterValue(&onwater,ElementOnWaterEnum);
-
-	/*If on water, skip: */
-	if(onwater)return;
-
-	/*Get node data: */
-	GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
-
-	/*pressure is lithostatic: */
-	//md.pressure=md.rho_ice*md.g*(md.surface-md.z); a la matlab
-
-	/*Get dof list on which we will plug the pressure values: */
-	GetDofList1(&doflist[0]);
-
-	/*Retrieve all inputs we will be needing: */
-	surface_input->GetParameterValues(&surface[0],&gauss[0][0],6);
-
-	/*pressure is lithostatic: */
-	rho_ice=matpar->GetRhoIce();
-	g=matpar->GetG();
-	for(i=0;i<numgrids;i++){
-		pressure[i]=rho_ice*g*(surface[i]-xyz_list[i][2]);
-	}
-
-	/*plug local pressure values into global pressure vector: */
-	VecSetValues(pg,numgrids,doflist,(const double*)pressure,INSERT_VALUES);
-
 }
 /*}}}*/
@@ -4483,5 +4437,5 @@
 	const int  numvertices = 6;
 	double     value[numvertices];
-	double     gauss[numvertices][4]={{1,0,0,-1},{0,1,0,-1},{0,0,1,-1},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
+	GaussPenta *gauss              = NULL;
 
 	/*Recover input*/
@@ -4493,9 +4447,12 @@
 
 	/* Start looping on the number of vertices: */
+	gauss=new GaussPenta();
 	for (int iv=0;iv<numvertices;iv++){
-		input->GetParameterValue(&pvalue[iv],&gauss[iv][0]);
+		gauss->GaussVertex(iv);
+		input->GetParameterValue(&pvalue[iv],gauss);
 	}
 
 	/*clean-up*/
+	delete gauss;
 }
 /*}}}*/
@@ -4506,5 +4463,5 @@
 	const int  numvertices = 6;
 	double     value[numvertices];
-	double     gauss[numvertices][4]={{1,0,0,-1},{0,1,0,-1},{0,0,1,-1},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
+	GaussPenta *gauss              = NULL;
 
 	/*Recover input*/
@@ -4515,11 +4472,17 @@
 
 	/* Start looping on the number of vertices: */
-	if (input)
-	 for (int iv=0;iv<numvertices;iv++) input->GetParameterValue(&pvalue[iv],&gauss[iv][0]);
-	else
-	 for (int iv=0;iv<numvertices;iv++) pvalue[iv]=defaultvalue;
+	if (input){
+		gauss=new GaussPenta();
+		for (int iv=0;iv<numvertices;iv++){
+			gauss->GaussVertex(iv);
+			input->GetParameterValue(&pvalue[iv],gauss);
+		}
+	}
+	else{
+		for (int iv=0;iv<numvertices;iv++) pvalue[iv]=defaultvalue;
+	}
 
 	/*clean-up*/
-	//delete gauss;
+	delete gauss;
 }
 /*}}}*/
@@ -5293,8 +5256,6 @@
 
 		/*Now compute pressure*/
-		penta->inputs->GetParameterValues(&surface[0],&gauss[0][0],6,SurfaceEnum);
-		for(i=0;i<numvertices;i++){
-			pressure[i]=rho_ice*g*(surface[i]-xyz_list[i][2]);
-		}
+		GetParameterListOnVertices(&surface[0],SurfaceEnum);
+		for(i=0;i<numvertices;i++) pressure[i]=rho_ice*g*(surface[i]-xyz_list[i][2]);
 
 		/*Now, we have to move the previous Vx and Vy inputs  to old 
@@ -5399,9 +5360,7 @@
 	rho_ice=matpar->GetRhoIce();
 	g=matpar->GetG();
-	inputs->GetParameterValues(&surface[0],&gauss[0][0],6,SurfaceEnum);
-
-	for(i=0;i<numvertices;i++){
-		pressure[i]=rho_ice*g*(surface[i]-xyz_list[i][2]);
-	}
+	GetParameterListOnVertices(&surface[0],SurfaceEnum);
+	for(i=0;i<numvertices;i++) pressure[i]=rho_ice*g*(surface[i]-xyz_list[i][2]);
+
 	/*Now, we have to move the previous Vx and Vy inputs  to old 
 	 * status, otherwise, we'll wipe them off: */
@@ -5482,9 +5441,7 @@
 	rho_ice=matpar->GetRhoIce();
 	g=matpar->GetG();
-	inputs->GetParameterValues(&surface[0],&gauss[0][0],6,SurfaceEnum);
-
-	for(i=0;i<numvertices;i++){
-		pressure[i]=rho_ice*g*(surface[i]-xyz_list[i][2]);
-	}
+	GetParameterListOnVertices(&surface[0],SurfaceEnum);
+	for(i=0;i<numvertices;i++) pressure[i]=rho_ice*g*(surface[i]-xyz_list[i][2]);
+
 	/*Now, we have to move the previous Vx and Vy inputs  to old 
 	 * status, otherwise, we'll wipe them off: */
@@ -5565,9 +5522,7 @@
 	rho_ice=matpar->GetRhoIce();
 	g=matpar->GetG();
-	inputs->GetParameterValues(&surface[0],&gauss[0][0],6,SurfaceEnum);
-	
-	for(i=0;i<numvertices;i++){
-		pressure[i]=rho_ice*g*(surface[i]-xyz_list[i][2]);
-	}
+	GetParameterListOnVertices(&surface[0],SurfaceEnum);
+	for(i=0;i<numvertices;i++) pressure[i]=rho_ice*g*(surface[i]-xyz_list[i][2]);
+
 	/*Now, we have to move the previous Vx and Vy inputs  to old 
 	 * status, otherwise, we'll wipe them off: */
@@ -5653,9 +5608,7 @@
 	rho_ice=matpar->GetRhoIce();
 	g=matpar->GetG();
-	inputs->GetParameterValues(&surface[0],&gauss[0][0],6,SurfaceEnum);
-
-	for(i=0;i<numvertices;i++){
-		pressure[i]=rho_ice*g*(surface[i]-xyz_list[i][2]);
-	}
+	GetParameterListOnVertices(&surface[0],SurfaceEnum);
+	for(i=0;i<numvertices;i++) pressure[i]=rho_ice*g*(surface[i]-xyz_list[i][2]);
+
 	/*Now, we have to move the previous Vz inputs to old 
 	 * status, otherwise, we'll wipe them off: */
Index: /issm/trunk/src/c/objects/Elements/Penta.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.h	(revision 5659)
+++ /issm/trunk/src/c/objects/Elements/Penta.h	(revision 5660)
@@ -71,5 +71,4 @@
 		void   AverageOntoPartition(Vec partition_contributions,Vec partition_areas,double* vertex_response,double* qmu_part);
 		void   ComputeBasalStress(Vec sigma_b);
-		void   ComputePressure(Vec p_g);
 		void   ComputeStrainRate(Vec eps);
 		void   Configure(Elements* elements,Loads* loads,DataSet* nodes,Materials* materials,Parameters* parameters);
Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 5659)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 5660)
@@ -575,28 +575,4 @@
 }
 /*}}}*/
-/*FUNCTION Tria::ComputePressure {{{1*/
-void  Tria::ComputePressure(Vec pg){
-
-	const int numvertices= 3;
-	int       doflist[numvertices];
-	double    pressure[numvertices];
-	double    thickness[numvertices];
-	double    rho_ice,g;
-	
-	/*Get dof list on which we will plug the pressure values: */
-	GetDofList1(&doflist[0]);
-
-	/*pressure is lithostatic: */
-	rho_ice=matpar->GetRhoIce();
-	g=matpar->GetG();
-	GetParameterListOnVertices(&thickness[0],ThicknessEnum);
-
-	for(int i=0;i<numvertices;i++) pressure[i]=rho_ice*g*thickness[i];
-
-	/*plug local pressure values into global pressure vector: */
-	VecSetValues(pg,numvertices,doflist,(const double*)pressure,INSERT_VALUES);
-
-}
-/*}}}*/
 /*FUNCTION Tria::ComputeStrainRate {{{1*/
 void  Tria::ComputeStrainRate(Vec eps){
@@ -5614,5 +5590,5 @@
 	/* Start looping on the number of vertices: */
 	gauss=new GaussTria();
-	for (int iv=0;iv<3;iv++){
+	for (int iv=0;iv<numvertices;iv++){
 		gauss->GaussVertex(iv);
 		input->GetParameterValue(&pvalue[iv],gauss);
Index: /issm/trunk/src/c/objects/Elements/Tria.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.h	(revision 5659)
+++ /issm/trunk/src/c/objects/Elements/Tria.h	(revision 5660)
@@ -68,5 +68,4 @@
 		void   AverageOntoPartition(Vec partition_contributions,Vec partition_areas,double* vertex_response,double* qmu_part);
 		void   ComputeBasalStress(Vec sigma_b);
-		void   ComputePressure(Vec p_g);
 		void   ComputeStrainRate(Vec eps);
 		void   Configure(Elements* elements,Loads* loads,DataSet* nodes,Materials* materials,Parameters* parameters);
Index: /issm/trunk/src/c/objects/Inputs/BoolInput.cpp
===================================================================
--- /issm/trunk/src/c/objects/Inputs/BoolInput.cpp	(revision 5659)
+++ /issm/trunk/src/c/objects/Inputs/BoolInput.cpp	(revision 5660)
@@ -175,7 +175,4 @@
 void BoolInput::GetParameterValue(double* pvalue,GaussPenta* gauss){ISSMERROR(" not supported yet!");}
 /*}}}*/
-/*FUNCTION BoolInput::GetParameterValues{{{1*/
-void BoolInput::GetParameterValues(double* values,double* gauss_pointers, int numgauss){ISSMERROR(" not supported yet!");}
-/*}}}*/
 /*FUNCTION BoolInput::GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, double* gauss){{{1*/
 void BoolInput::GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, double* gauss){ISSMERROR(" not supported yet!");}
Index: /issm/trunk/src/c/objects/Inputs/BoolInput.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/BoolInput.h	(revision 5659)
+++ /issm/trunk/src/c/objects/Inputs/BoolInput.h	(revision 5660)
@@ -50,5 +50,4 @@
 		void GetParameterValue(double* pvalue,GaussTria* gauss);
 		void GetParameterValue(double* pvalue,GaussPenta* gauss);
-		void GetParameterValues(double* values,double* gauss_pointers, int numgauss);
 		void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, double* gauss);
 		void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, GaussTria* gauss);
Index: /issm/trunk/src/c/objects/Inputs/DoubleInput.cpp
===================================================================
--- /issm/trunk/src/c/objects/Inputs/DoubleInput.cpp	(revision 5659)
+++ /issm/trunk/src/c/objects/Inputs/DoubleInput.cpp	(revision 5660)
@@ -188,7 +188,4 @@
 void DoubleInput::GetParameterValue(double* pvalue,GaussPenta* gauss){ISSMERROR(" not supported yet!");}
 /*}}}*/
-/*FUNCTION DoubleInput::GetParameterValues{{{1*/
-void DoubleInput::GetParameterValues(double* values,double* gauss_pointers, int numgauss){ISSMERROR(" not supported yet!");}
-/*}}}*/
 /*FUNCTION DoubleInput::GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, double* gauss){{{1*/
 void DoubleInput::GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, double* gauss){ISSMERROR(" not supported yet!");}
Index: /issm/trunk/src/c/objects/Inputs/DoubleInput.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/DoubleInput.h	(revision 5659)
+++ /issm/trunk/src/c/objects/Inputs/DoubleInput.h	(revision 5660)
@@ -49,5 +49,4 @@
 		void GetParameterValue(double* pvalue,GaussTria* gauss);
 		void GetParameterValue(double* pvalue,GaussPenta* gauss);
-		void GetParameterValues(double* values,double* gauss_pointers, int numgauss);
 		void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, double* gauss);
 		void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, GaussTria* gauss);
Index: /issm/trunk/src/c/objects/Inputs/Input.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/Input.h	(revision 5659)
+++ /issm/trunk/src/c/objects/Inputs/Input.h	(revision 5660)
@@ -28,5 +28,4 @@
 		virtual void GetParameterValue(double* pvalue,GaussTria* gauss)=0;
 		virtual void GetParameterValue(double* pvalue,GaussPenta* gauss)=0;
-		virtual void GetParameterValues(double* values,double* gauss_pointers, int numgauss)=0;
 		virtual void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, double* gauss)=0;
 		virtual void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, GaussTria* gauss)=0;
Index: /issm/trunk/src/c/objects/Inputs/IntInput.cpp
===================================================================
--- /issm/trunk/src/c/objects/Inputs/IntInput.cpp	(revision 5659)
+++ /issm/trunk/src/c/objects/Inputs/IntInput.cpp	(revision 5660)
@@ -176,7 +176,4 @@
 void IntInput::GetParameterValue(double* pvalue,GaussPenta* gauss){ISSMERROR(" not supported yet!");}
 /*}}}*/
-/*FUNCTION IntInput::GetParameterValues{{{1*/
-void IntInput::GetParameterValues(double* values,double* gauss_pointers, int numgauss){ISSMERROR(" not supported yet!");}
-/*}}}*/
 /*FUNCTION IntInput::GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, double* gauss){{{1*/
 void IntInput::GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, double* gauss){ISSMERROR(" not supported yet!");}
Index: /issm/trunk/src/c/objects/Inputs/IntInput.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/IntInput.h	(revision 5659)
+++ /issm/trunk/src/c/objects/Inputs/IntInput.h	(revision 5660)
@@ -50,5 +50,4 @@
 		void GetParameterValue(double* pvalue,GaussTria* gauss);
 		void GetParameterValue(double* pvalue,GaussPenta* gauss);
-		void GetParameterValues(double* values,double* gauss_pointers, int numgauss);
 		void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, double* gauss);
 		void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, GaussTria* gauss);
Index: /issm/trunk/src/c/objects/Inputs/PentaVertexInput.cpp
===================================================================
--- /issm/trunk/src/c/objects/Inputs/PentaVertexInput.cpp	(revision 5659)
+++ /issm/trunk/src/c/objects/Inputs/PentaVertexInput.cpp	(revision 5660)
@@ -192,21 +192,4 @@
 }
 /*}}}*/
-/*FUNCTION PentaVertexInput::GetParameterValues{{{1*/
-void PentaVertexInput::GetParameterValues(double* values,double* gauss_pointers, int numgauss){
-	/*It is assumed that output values has been correctly allocated*/
-
-	int i,j;
-	double gauss[4];
-
-	for (i=0;i<numgauss;i++){
-
-		/*Get current Gauss point coordinates*/
-		for (j=0;j<4;j++) gauss[j]=gauss_pointers[i*4+j];
-
-		/*Assign parameter value*/
-		GetParameterValue(&values[i],&gauss[0]);
-	}
-}
-/*}}}*/
 /*FUNCTION PentaVertexInput::GetParameterDerivativeValue(double* p, double* xyz_list, double* gauss){{{1*/
 void PentaVertexInput::GetParameterDerivativeValue(double* p, double* xyz_list, double* gauss){
Index: /issm/trunk/src/c/objects/Inputs/PentaVertexInput.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/PentaVertexInput.h	(revision 5659)
+++ /issm/trunk/src/c/objects/Inputs/PentaVertexInput.h	(revision 5660)
@@ -50,5 +50,4 @@
 		void GetParameterValue(double* pvalue,GaussTria* gauss){ISSMERROR("not implemented yet");};
 		void GetParameterValue(double* pvalue,GaussPenta* gauss);
-		void GetParameterValues(double* values,double* gauss_pointers, int numgauss);
 		void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, double* gauss);
 		void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, GaussTria* gauss){ISSMERROR("not implemented yet");};
Index: /issm/trunk/src/c/objects/Inputs/TriaVertexInput.cpp
===================================================================
--- /issm/trunk/src/c/objects/Inputs/TriaVertexInput.cpp	(revision 5659)
+++ /issm/trunk/src/c/objects/Inputs/TriaVertexInput.cpp	(revision 5660)
@@ -181,22 +181,4 @@
 }
 /*}}}*/
-/*FUNCTION TriaVertexInput::GetParameterValues{{{1*/
-void TriaVertexInput::GetParameterValues(double* values,double* gauss_pointers, int numgauss){
-	/*It is assumed that output values has been correctly allocated*/
-
-	int i,j;
-	double gauss[3];
-
-	for (i=0;i<numgauss;i++){
-
-		/*Get current Gauss point coordinates*/
-		for (j=0;j<3;j++) gauss[j]=gauss_pointers[i*3+j];
-
-		/*Assign parameter value*/
-		GetParameterValue(&values[i],&gauss[0]);
-	}
-
-}
-/*}}}*/
 /*FUNCTION TriaVertexInput::GetParameterDerivativeValue(double* p, double* xyz_list, double* gauss){{{1*/
 void TriaVertexInput::GetParameterDerivativeValue(double* p, double* xyz_list, double* gauss){
Index: /issm/trunk/src/c/objects/Inputs/TriaVertexInput.h
===================================================================
--- /issm/trunk/src/c/objects/Inputs/TriaVertexInput.h	(revision 5659)
+++ /issm/trunk/src/c/objects/Inputs/TriaVertexInput.h	(revision 5660)
@@ -50,5 +50,4 @@
 		void GetParameterValue(double* pvalue,GaussTria* gauss);
 		void GetParameterValue(double* pvalue,GaussPenta* gauss){ISSMERROR("not implemented yet");};
-		void GetParameterValues(double* values,double* gauss_pointers, int numgauss);
 		void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, double* gauss);
 		void GetParameterDerivativeValue(double* derivativevalues, double* xyz_list, GaussTria* gauss);
Index: /issm/trunk/src/mex/Makefile.am
===================================================================
--- /issm/trunk/src/mex/Makefile.am	(revision 5659)
+++ /issm/trunk/src/mex/Makefile.am	(revision 5660)
@@ -12,5 +12,4 @@
 				Chaco\
 				ComputeBasalStress\
-				ComputePressure\
 				ConfigureObjects \
 				ControlOptimization\