Index: /issm/trunk/src/c/objects/Pengrid.cpp
===================================================================
--- /issm/trunk/src/c/objects/Pengrid.cpp	(revision 467)
+++ /issm/trunk/src/c/objects/Pengrid.cpp	(revision 468)
@@ -233,5 +233,5 @@
 #undef __FUNCT__ 
 #define __FUNCT__ "Pengrid::PenaltyCreatePVector"
-void  Pengrid::PenaltyCreatePVector(Vec pg,void* inputs,double kmax,int analysis_type){
+void  Pengrid::PenaltyCreatePVector(Vec pg,void* inputs,double kmax,int analysis_type,int sub_analysis_type){
 
 	/*No penalty applied, do nothing: */
Index: /issm/trunk/src/c/objects/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Penta.cpp	(revision 467)
+++ /issm/trunk/src/c/objects/Penta.cpp	(revision 468)
@@ -1112,5 +1112,5 @@
 		 * and compute Du*/
 		tria=(Tria*)SpawnTria(0,1,2); //grids 0, 1 and 2 make the new tria (lower face).
-		tria->Du(du_g,u_g,u_g_obs,inputs,analysis_type);
+		tria->Du(du_g,u_g,u_g_obs,inputs,analysis_type,sub_analysis_type);
 		delete tria;
 		return;
@@ -1119,5 +1119,5 @@
 		
 		tria=(Tria*)SpawnTria(3,4,5); //grids 3, 4 and 5 make the new tria (upper face).
-		tria->Du(du_g,u_g,u_g_obs,inputs,analysis_type);
+		tria->Du(du_g,u_g,u_g_obs,inputs,analysis_type,sub_analysis_type);
 		delete tria;
 		return;
@@ -1182,5 +1182,5 @@
 		 * and compute Misfit*/
 		tria=(Tria*)SpawnTria(0,1,2); //grids 0, 1 and 2 make the new tria (lower face).
-		J=tria->Misfit( velocity,obs_velocity,inputs,analysis_type);
+		J=tria->Misfit( velocity,obs_velocity,inputs,analysis_type,sub_analysis_type);
 		delete tria;
 		return J;
@@ -1189,5 +1189,5 @@
 
 		tria=(Tria*)SpawnTria(3,4,5); //grids 3, 4 and 5 make the new tria (upper face).
-		J=tria->Misfit( velocity,obs_velocity,inputs,analysis_type);
+		J=tria->Misfit( velocity,obs_velocity,inputs,analysis_type,sub_analysis_type);
 		delete tria;
 		return J;