Index: /issm/trunk/src/m/solutions/cielo/modelsize.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/modelsize.m	(revision 808)
+++ /issm/trunk/src/m/solutions/cielo/modelsize.m	(revision 809)
@@ -1,11 +1,20 @@
 function dof=modelsize(m_dh,m_dv,m_ds,m_dhu,m_sl);
+%DOF - return the maximum number of degrees of freedom of the model
+%
+%   Usage:
+%      dof=modelsize(m_dh,m_dv,m_ds,m_dhu,m_sl)
 
-
+dof=0;
+if ~isempty(m_ds.nodesets),
+	dof=m_ds.nodesets.fsize; %biggest dof number
+end
 if ~isempty(m_dh.nodesets),
-	dof=m_dh.nodesets.fsize; %biggest dof number
-else
-	if m_dhu.parameters.ishutter,
+	if dof<m_dh.nodesets.fsize,
+		dof=m_dh.nodesets.fsize; %biggest dof number
+	end
+end
+if ~isempty(m_dhu.nodesets),
+	if dof<m_dhu.nodesets.fsize,
 		dof=m_dhu.nodesets.fsize;
 	end
 end
-