Index: /issm/trunk/src/m/classes/@acceleratedtriaelem/CreateKMatrix.m
===================================================================
--- /issm/trunk/src/m/classes/@acceleratedtriaelem/CreateKMatrix.m	(revision 38)
+++ /issm/trunk/src/m/classes/@acceleratedtriaelem/CreateKMatrix.m	(revision 39)
@@ -166,5 +166,5 @@
 end
 
-Drag_operator=sparse2(6*element.nods,6*element.nods);
+Drag_operator=sparse(6*element.nods,6*element.nods);
 
 %Create basal stiffness
@@ -228,5 +228,5 @@
 	end  
 
-	Drag_operator=sparse2(row_location_drag,col_location_drag,drag_operator_value,6*element.nods,6*element.nods);
+	Drag_operator=sparse(row_location_drag,col_location_drag,drag_operator_value,6*element.nods,6*element.nods);
 end
 
@@ -234,3 +234,3 @@
 Ke=globalmatrix;
 Ke.nrows=numdof*element.nods;
-Ke.matrix=Drag_operator+sparse2(row_location,col_location,Ke_terms,6*element.nods,6*element.nods);
+Ke.matrix=Drag_operator+sparse(row_location,col_location,Ke_terms,6*element.nods,6*element.nods);
Index: /issm/trunk/src/m/classes/@acceleratedtriaelem/CreatePVector.m
===================================================================
--- /issm/trunk/src/m/classes/@acceleratedtriaelem/CreatePVector.m	(revision 38)
+++ /issm/trunk/src/m/classes/@acceleratedtriaelem/CreatePVector.m	(revision 39)
@@ -69,6 +69,6 @@
 end
 
-Rhs=full([sparse2(row_rhs,ones(element.nel*27,1),Rhs_x,element.nods,1)
-   sparse2(row_rhs,ones(element.nel*27,1),Rhs_y,element.nods,1)]);
+Rhs=full([sparse(row_rhs,ones(element.nel*27,1),Rhs_x,element.nods,1)
+   sparse(row_rhs,ones(element.nel*27,1),Rhs_y,element.nods,1)]);
 
 %Index icefront is skipped: 
Index: /issm/trunk/src/m/classes/public/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 38)
+++ /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 39)
@@ -371,5 +371,5 @@
 	%SPARSITY BETWEEN 0 AND 1
 	if ( (md.sparsity<=0) | (md.sparsity>1)),
-		disp(['model ' md.name ' sparsity should be inside the ]0 1] range']);
+		disp(['model ' md.name ' sparsity should be inside the [0 1] range']);
 		bool=0; return;
 	end
@@ -400,8 +400,10 @@
 
 	%NP
-	if md.np==0,
-		disp(['model ' md.name ' has a =0 value in field ' md.np ' !']);
-	elseif md.np<0,
-		disp(['model ' md.name ' has a negative value in field ' md.np ' !']);
+	if ~strcmpi(md.cluster,'none'),
+		if md.np==0,
+			disp(['model ' md.name ' has a =0 value in field np !']);
+		elseif md.np<0,
+			disp(['model ' md.name ' has a negative value in field np !']);
+		end
 	end
 
Index: /issm/trunk/src/m/classes/public/solve.m
===================================================================
--- /issm/trunk/src/m/classes/public/solve.m	(revision 38)
+++ /issm/trunk/src/m/classes/public/solve.m	(revision 39)
@@ -47,5 +47,4 @@
 %Use package to set solution namespace
 usenamespace(package);
-
 
 %Check model is self-consistent