Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15752)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15753)
@@ -9891,6 +9891,7 @@
 
 	const int    numdofm=NDOF2*NUMVERTICES;
-	const int    numdofs=NDOF4*NUMVERTICES;
+	const int    numdofs=NDOF3*NUMVERTICES;
 	const int    numdof2d=NDOF2*NUMVERTICES2D;
+	const int    numdofpressure=NDOF1*NUMVERTICES;
 
 	int     i;
@@ -9898,4 +9899,5 @@
 	IssmDouble  SSA_values[numdofm];
 	IssmDouble  FS_values[numdofs];
+	IssmDouble  Pressure_values[numdofs];
 	IssmDouble  vx[NUMVERTICES];
 	IssmDouble  vy[NUMVERTICES];
@@ -9908,4 +9910,5 @@
 	int*    doflistm        = NULL;
 	int*    doflists        = NULL;
+	int*    doflistpressure = NULL;
 	Penta   *penta          = NULL;
 
@@ -9917,4 +9920,5 @@
 	penta->GetDofList(&doflistm,SSAApproximationEnum,GsetEnum);
 	GetDofList(&doflists,FSvelocityEnum,GsetEnum);
+	GetDofListPressure(&doflistpressure,GsetEnum);
 	this->parameters->FindParam(&FSreconditioning,DiagnosticFSreconditioningEnum);
 
@@ -9927,7 +9931,6 @@
 		SSA_values[i+numdof2d]=solution[doflistm[i]];
 	}
-	for(i=0;i<numdofs;i++){
-		FS_values[i]=solution[doflists[i]];
-	}
+	for(i=0;i<numdofs;i++)FS_values[i]=solution[doflists[i]];
+	for(i=0;i<numdofpressure;i++) Pressure_values[i]=solution[doflistpressure[i]];
 
 	/*Transform solution in Cartesian Space*/
@@ -9937,8 +9940,8 @@
 	/*Ok, we have vx and vy in values, fill in vx and vy arrays: */
 	for(i=0;i<NUMVERTICES;i++){
-		vx[i]=FS_values[i*NDOF4+0]+SSA_values[i*NDOF2+0];
-		vy[i]=FS_values[i*NDOF4+1]+SSA_values[i*NDOF2+1];
-		vzFS[i]=FS_values[i*NDOF4+2];
-		pressure[i]=FS_values[i*NDOF4+3]*FSreconditioning;
+		vx[i]=FS_values[i*NDOF3+0]+SSA_values[i*NDOF2+0];
+		vy[i]=FS_values[i*NDOF3+1]+SSA_values[i*NDOF2+1];
+		vzFS[i]=FS_values[i*NDOF3+2];
+		pressure[i]=Pressure_values[i*NDOF1]*FSreconditioning;
 
 		/*Check solution*/
@@ -10193,5 +10196,5 @@
 		vy[i]=FS_values[i*NDOF3+1]+HO_values[i*NDOF2+1];
 		vzFS[i]=FS_values[i*NDOF3+2];
-		pressure[i]=FS_values[i*NDOF1]*FSreconditioning;
+		pressure[i]=Pressure_values[i*NDOF1]*FSreconditioning;
 
 		/*Check solution*/