Index: /issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp	(revision 322)
+++ /issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp	(revision 323)
@@ -54,5 +54,5 @@
 	#endif
 
-		if ((int)model->gridonhutter[i]){
+		if (!(int)model->gridonhutter[i]){
 	
 			spc_sid=count;
Index: /issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateElementsNodesAndMaterialsDiagnosticHutter.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateElementsNodesAndMaterialsDiagnosticHutter.cpp	(revision 322)
+++ /issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateElementsNodesAndMaterialsDiagnosticHutter.cpp	(revision 323)
@@ -231,5 +231,5 @@
 					matice_mid=i+1; //same as the material id from the geom2 elements.
 					matice_B=model->B[i];	
-					matice_n=(double)model->n[i];
+					matice_n=(double)model->n[1]; //n defined on elements not grids, so take the first value everywhere
 				
 					/*Create matice ubeam its constructor:*/
@@ -332,13 +332,4 @@
 		node=new Node(node_id,node_partitionborder,node_numdofs,node_x,node_onbed,node_onsurface,node_upper_node_id);
 
-		/*set single point constraints.: */
-		if (strcmp(model->meshtype,"3d")==0){
-			/*On a 3d mesh, we may have collapsed elements, hence dead grids. Freeze them out: */
-			if (model->deadgrids[i]){
-				for(k=1;k<=node_numdofs;k++){
-					node->FreezeDof(k);
-				}
-			}
-		}
 		/*Add node to nodes dataset: */
 		nodes->AddObject(node);