Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 10525)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 10526)
@@ -5547,8 +5547,8 @@
 
 	/*Constants*/
-	const int    numdofp=NDOF2*NUMVERTICES;
-	const int    numdofs=NDOF4*NUMVERTICES;
-	const int    numdoftotal=(NDOF2+NDOF4)*NUMVERTICES;
-	int          i,j;
+	const int numdofp     = NDOF2 *NUMVERTICES;
+	const int numdofs     = NDOF4 *NUMVERTICES;
+	const int numdoftotal = (NDOF2+NDOF4) *NUMVERTICES;
+	int       i,j;
 
 	for(i=0;i<numdofs;i++) for(j=0;j<NUMVERTICES;j++){
@@ -6469,4 +6469,7 @@
 	}
 
+	/*Transform coordinate system*/
+	TransformLoadVectorCoord(pe,nodes,NUMVERTICES,XYZPEnum);
+
 	/*Clean up and return*/
 	delete gauss;
@@ -6539,4 +6542,7 @@
 		}
 	}
+
+	/*Transform coordinate system*/
+	TransformLoadVectorCoord(pe,nodes,NUMVERTICES,XYZPEnum);
 
 	/*Clean up and return*/
@@ -7632,4 +7638,8 @@
 	for(i=0;i<numdofs;i++) stokes_values[i]=solution[doflists[i]];
 
+	/*Transform solution in Cartesian Space*/
+	TransformSolutionCoord(&pattyn_values[0],this->nodes,NUMVERTICES,XYEnum);
+	TransformSolutionCoord(&stokes_values[0],this->nodes,NUMVERTICES,XYZPEnum);
+
 	/*Ok, we have vx and vy in values, fill in vx and vy arrays: */
 	for(i=0;i<NUMVERTICES;i++){