Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 6933)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 6934)
@@ -1752,5 +1752,4 @@
 ElementVector* Penta::CreatePVectorCouplingMacAyealStokes(void){
 
-	_error_("coupling MacAyeal Stokes not implemented yet");
 	/*compute all load vectors for this element*/
 	ElementVector* pe1=CreatePVectorCouplingMacAyealStokesViscous();
@@ -1794,5 +1793,5 @@
 	Input* vy_input=inputs->GetInput(VyEnum);               _assert_(vy_input);
 	Input* vz_input=inputs->GetInput(VzEnum);               _assert_(vz_input);
-	Input* vzpattyn_input=inputs->GetInput(VzMacAyealEnum);   _assert_(vzpattyn_input);
+	Input* vzmacayeal_input=inputs->GetInput(VzMacAyealEnum);   _assert_(vzmacayeal_input);
 
 	/* Start  looping on the number of gaussian points: */
@@ -1806,5 +1805,5 @@
 		GetNodalFunctionsP1Derivatives(&dh1dh6[0][0],&xyz_list[0][0], gauss);
 
-		vzpattyn_input->GetParameterDerivativeValue(&dw[0],&xyz_list[0][0],gauss);
+		vzmacayeal_input->GetParameterDerivativeValue(&dw[0],&xyz_list[0][0],gauss);
 
 		this->GetStrainRate3d(&epsilon[0],&xyz_list[0][0],gauss,vx_input,vy_input,vz_input);
@@ -1859,5 +1858,5 @@
 	Input* vy_input=inputs->GetInput(VyEnum);               _assert_(vy_input);
 	Input* vz_input=inputs->GetInput(VzEnum);               _assert_(vz_input);
-	Input* vzpattyn_input=inputs->GetInput(VzMacAyealEnum);   _assert_(vzpattyn_input);
+	Input* vzmacayeal_input=inputs->GetInput(VzMacAyealEnum);   _assert_(vzmacayeal_input);
 
 	for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
@@ -1875,6 +1874,6 @@
 		GetNodalFunctionsP1(l1l6, gauss);
 
-		vzpattyn_input->GetParameterValue(&w, gauss);
-		vzpattyn_input->GetParameterDerivativeValue(&dw[0],&xyz_list[0][0],gauss);
+		vzmacayeal_input->GetParameterValue(&w, gauss);
+		vzmacayeal_input->GetParameterDerivativeValue(&dw[0],&xyz_list[0][0],gauss);
 
 		BedNormal(&bed_normal[0],xyz_list_tria);
@@ -2460,5 +2459,5 @@
 	Input* vy_input=inputs->GetInput(VyEnum); _assert_(vy_input);
 	Input* vzstokes_input=NULL;
-	if(approximation==PattynStokesApproximationEnum){
+	if(approximation==PattynStokesApproximationEnum || approximation==MacAyealStokesApproximationEnum){
 		vzstokes_input=inputs->GetInput(VzStokesEnum); _assert_(vzstokes_input);
 	}
@@ -2475,5 +2474,5 @@
 		vx_input->GetParameterDerivativeValue(&du[0],&xyz_list[0][0],gauss);
 		vy_input->GetParameterDerivativeValue(&dv[0],&xyz_list[0][0],gauss);
-		if(approximation==PattynStokesApproximationEnum){
+		if(approximation==PattynStokesApproximationEnum || approximation==MacAyealStokesApproximationEnum){
 			vzstokes_input->GetParameterDerivativeValue(&dw[0],&xyz_list[0][0],gauss);
 			dwdz=dw[2];
Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 6933)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 6934)
@@ -1644,5 +1644,5 @@
 	Input* vy_input=inputs->GetInput(VyEnum);               _assert_(vy_input);
 	Input* vzstokes_input=NULL;
-	if(approximation==PattynStokesApproximationEnum){
+	if(approximation==PattynStokesApproximationEnum || approximation==MacAyealStokesApproximationEnum){
 		vzstokes_input=inputs->GetInput(VzStokesEnum);       _assert_(vzstokes_input);
 	}
@@ -1658,5 +1658,5 @@
 		vx_input->GetParameterValue(&vx, gauss);
 		vy_input->GetParameterValue(&vy, gauss);
-		if(approximation==PattynStokesApproximationEnum){
+		if(approximation==PattynStokesApproximationEnum || approximation==MacAyealStokesApproximationEnum){
 			vzstokes_input->GetParameterValue(&vz, gauss);
 		}