Index: /issm/trunk/src/m/classes/@node/display.m
===================================================================
--- /issm/trunk/src/m/classes/@node/display.m	(revision 1041)
+++ /issm/trunk/src/m/classes/@node/display.m	(revision 1042)
@@ -10,4 +10,5 @@
 disp(sprintf('   z: %g',node.z));
 disp(sprintf('   s: %g',node.s));
+disp(sprintf('   surface: %g',node.surface));
 disp(sprintf('   doflist: %i %i %i %i %i %i',node.doflist(1),node.doflist(2),node.doflist(3),node.doflist(4),node.doflist(5),node.doflist(6)));
 disp(sprintf('   gridset: %s',node.gridset));
Index: /issm/trunk/src/m/classes/@node/node.m
===================================================================
--- /issm/trunk/src/m/classes/@node/node.m	(revision 1041)
+++ /issm/trunk/src/m/classes/@node/node.m	(revision 1042)
@@ -13,4 +13,5 @@
 	node.z=NaN;
 	node.s=NaN;
+	node.surface=NaN;
 	node.doflist=zeros(6,1); %6 degrees of freedom available per node.
 	node.gridset='';
Index: /issm/trunk/src/m/solutions/cielo/diagnostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic_core.m	(revision 1041)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic_core.m	(revision 1042)
@@ -40,5 +40,5 @@
 
 	displaystring(debug,'\n%s',['computing pressure according to MacAyeal...']);
-	p_g=ComputePressure(m_dhu.elements,m_dhu.nodes,m_dhu.loads,m_dhu.materials,m_dhu.parameters);
+	p_g=ComputePressure(m_dhu.elements,m_dhu.nodes,m_dhu.loads,m_dhu.materials,m_dhu.parameters,inputs);
 
 	displaystring(debug,'\n%s',['update boundary conditions for macyeal pattyn using hutter results...']);
@@ -56,5 +56,5 @@
 
 	displaystring(debug,'\n%s',['computing pressure according to MacAyeal...']);
-	p_g=ComputePressure(m_dh.elements,m_dh.nodes,m_dh.loads,m_dh.materials,m_dh.parameters);
+	p_g=ComputePressure(m_dh.elements,m_dh.nodes,m_dh.loads,m_dh.materials,m_dh.parameters,inputs);
 end
 	
@@ -74,5 +74,5 @@
 
 	displaystring(debug,'\n%s',['computing pressure according to Pattyn...']);
-	p_g=ComputePressure(m_dh.elements,m_dh.nodes,m_dh.loads,m_dh.materials,m_dh.parameters);
+	p_g=ComputePressure(m_dh.elements,m_dh.nodes,m_dh.loads,m_dh.materials,m_dh.parameters,inputs);
 	
 	if isstokes,
Index: /issm/trunk/src/mex/ComputePressure/ComputePressure.cpp
===================================================================
--- /issm/trunk/src/mex/ComputePressure/ComputePressure.cpp	(revision 1041)
+++ /issm/trunk/src/mex/ComputePressure/ComputePressure.cpp	(revision 1042)
@@ -15,4 +15,5 @@
 	DataSet* loads=NULL;
 	DataSet* materials=NULL;
+	ParameterInputs* inputs=NULL;
 	int      numberofnodes;
 
@@ -34,4 +35,11 @@
 
 	/*!Generate internal degree of freedom numbers: */
+	/*Fetch inputs: */
+	inputs=new ParameterInputs;
+	inputs->Init(INPUTS);
+
+	UpdateFromInputsx(elements,nodes,loads, materials,inputs);
+
+	/*!Generate internal degree of freedom numbers: */
 	ComputePressurex(&p_g, elements,nodes,loads,materials,numberofnodes);
 
@@ -44,4 +52,5 @@
 	delete materials;
 	delete loads;
+	delete inputs;
 	VecFree(&p_g);
 
@@ -52,5 +61,5 @@
 void ComputePressureUsage(void) {
 	_printf_("\n");
-	_printf_("   usage: [p_g] = %s(elements, nodes, loads, materials, params);\n",__FUNCT__);
+	_printf_("   usage: [p_g] = %s(elements, nodes, loads, materials, params,inputs);\n",__FUNCT__);
 	_printf_("\n");
 }
Index: /issm/trunk/src/mex/ComputePressure/ComputePressure.h
===================================================================
--- /issm/trunk/src/mex/ComputePressure/ComputePressure.h	(revision 1041)
+++ /issm/trunk/src/mex/ComputePressure/ComputePressure.h	(revision 1042)
@@ -22,4 +22,5 @@
 #define MATERIALS (mxArray*)prhs[3]
 #define PARAMS (mxArray*)prhs[4]
+#define INPUTS (mxArray*)prhs[5]
 
 /* serial output macros: */
@@ -30,5 +31,5 @@
 #define NLHS  1
 #undef NRHS
-#define NRHS  5
+#define NRHS  6