Index: /issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp	(revision 559)
+++ /issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp	(revision 560)
@@ -83,43 +83,4 @@
 	xfree((void**)&model->gridonhutter);
 
-	//deal with mpcs for 2d-3d mesh transitions
-	
-	/*Fetch data: */
-	ModelFetchData((void**)&model->penalties,&model->numpenalties,NULL,model_handle,"penalties","Matrix","Mat");
-
-	if (strcmp(model->meshtype,"3d")==0){
-	
-		if (model->numpenalties){
-			for (i=0;i<model->numpenalties;i++){
-				for (j=0;j<model->numlayers-1;j++){
-
-					//constrain first dof
-					rgb_id=count;
-					rgb_dof=1;
-					rgb_nodeid1=(int)*(model->penalties+(model->numlayers-1)*i+j);
-					rgb_nodeid1=(int)*(model->penalties+(model->numlayers-1)*i+j+1);
-					
-					rgb = new Rgb(rgb_id,rgb_nodeid1,rgb_nodeid2,rgb_dof);
-					constraints->AddObject(rgb);
-					count++;
-
-					//constrain second dof
-					rgb_id=count;
-					rgb_dof=2;
-					rgb_nodeid1=(int)*(model->penalties+(model->numlayers-1)*i+j);
-					rgb_nodeid1=(int)*(model->penalties+(model->numlayers-1)*i+j+1);
-					
-					rgb = new Rgb(rgb_id,rgb_nodeid1,rgb_nodeid2,rgb_dof);
-					constraints->AddObject(rgb);
-					count++;
-	
-				}
-			}
-		}
-	}
-
-	/*Free data: */
-	xfree((void**)&model->penalties);
-
 	/*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these 
 	 * datasets, it will not be redone: */