Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 4845)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 4846)
@@ -2151,7 +2151,4 @@
 	GetDofList(&doflist[0],&dummy);
 
-	/*Get node data: */
-	GetVerticesCoordinates(&xyz_list[0][0], nodes, numvertices);
-
 	/*Use the dof list to index into the solution vector: */
 	for(i=0;i<numdof;i++){
@@ -2174,4 +2171,7 @@
 	penta=this;
 	for(;;){
+
+		/*Get node data: */
+		GetVerticesCoordinates(&xyz_list[0][0],penta->nodes,numvertices);
 
 		/*Now Compute vel*/
@@ -2186,5 +2186,5 @@
 
 		/*Now compute pressure*/
-		inputs->GetParameterValues(&surface[0],&gauss[0][0],6,SurfaceEnum);
+		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]);
@@ -2448,4 +2448,5 @@
 	 * status, otherwise, we'll wipe them off: */
 	this->inputs->ChangeEnum(VzEnum,VzOldEnum);
+	this->inputs->ChangeEnum(PressureEnum,PressureOldEnum);
 
 	/*Add vz and vel as inputs to the penta element: */
Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 4845)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 4846)
@@ -5792,6 +5792,6 @@
 
 	int          dummy;
-	Input*       VzInput=NULL;
-	double*      VzPtr=NULL;
+	Input*       vz_input=NULL;
+	double*      vz_ptr=NULL;
 	
 	/*Get dof list: */
@@ -5810,11 +5810,11 @@
 
 	/*Get Vz*/
-	VzInput=inputs->GetInput(VzEnum);
-	if (VzInput){
-		if (VzInput->Enum()!=TriaVertexInputEnum){
-			ISSMERROR("Cannot compute Vel as Vz is of type %s",EnumAsString(VzInput->Enum()));
+	vz_input=inputs->GetInput(VzEnum);
+	if (vz_input){
+		if (vz_input->Enum()!=TriaVertexInputEnum){
+			ISSMERROR("Cannot compute Vel as Vz is of type %s",EnumAsString(vz_input->Enum()));
 		}
-		VzInput->GetValuesPtr(&VzPtr,&dummy);
-		for(i=0;i<numvertices;i++) vz[i]=VzPtr[i];
+		vz_input->GetValuesPtr(&vz_ptr,&dummy);
+		for(i=0;i<numvertices;i++) vz[i]=vz_ptr[i];
 	}
 	else{
@@ -5870,6 +5870,6 @@
 
 	int          dummy;
-	Input*       VzInput=NULL;
-	double*      VzPtr=NULL;
+	Input*       vz_input=NULL;
+	double*      vz_ptr=NULL;
 	
 	/*Get dof list: */
@@ -5888,11 +5888,11 @@
 
 	/*Get Vz*/
-	VzInput=inputs->GetInput(VzEnum);
-	if (VzInput){
-		if (VzInput->Enum()!=TriaVertexInputEnum){
-			ISSMERROR("Cannot compute Vel as Vz is of type %s",EnumAsString(VzInput->Enum()));
+	vz_input=inputs->GetInput(VzEnum);
+	if (vz_input){
+		if (vz_input->Enum()!=TriaVertexInputEnum){
+			ISSMERROR("Cannot compute Vel as Vz is of type %s",EnumAsString(vz_input->Enum()));
 		}
-		VzInput->GetValuesPtr(&VzPtr,&dummy);
-		for(i=0;i<numvertices;i++) vz[i]=VzPtr[i];
+		vz_input->GetValuesPtr(&vz_ptr,&dummy);
+		for(i=0;i<numvertices;i++) vz[i]=vz_ptr[i];
 	}
 	else{