Index: /issm/trunk/src/m/kml/nodeconnectivity.m
===================================================================
--- /issm/trunk/src/m/kml/nodeconnectivity.m	(revision 6346)
+++ /issm/trunk/src/m/kml/nodeconnectivity.m	(revision 6346)
@@ -0,0 +1,58 @@
+%
+%  create a node connectivity table for the elements in the model.
+%
+%  [nodecon]=edgeadjacency(elem,ngrids,mxepg)
+%
+%  where the required input is:
+%    elem          (numeric, element connectivity array (elems x nodes))
+%
+%  and the required output is:
+%    nodecon       (numeric, node connectivity array (ngrids x mxepg+1))
+%
+%  the optional input is:
+%    ngrids        (numeric, number of grids)
+%    mxepg         (numeric, max elements per grid)
+%
+function [nodecon]=nodeconnectivity(elem,ngrids,mxepg)
+
+if ~nargin
+    help nodeconnectivity
+    return
+end
+
+if ~exist('ngrids','var') || isempty(ngrids)
+    ngrids=max(max(elem));
+end
+if ~exist('mxepg','var') || isempty(mxepg)
+    mxepg=25;
+end
+
+%%  create the node connectivity array
+
+nodecon=zeros(ngrids,mxepg+1);
+
+%  loop over the elements
+
+for i=1:size(elem,1)
+
+%  loop over the nodes for each element
+
+    for j=1:size(elem,2)
+        if elem(i,j)
+            nodecon(elem(i,j),nodecon(elem(i,j),end)+1)=i;
+            nodecon(elem(i,j),end)=nodecon(elem(i,j),end)+1;
+        end
+    end
+end
+
+%%  sort the node connectivity array
+
+%  loop over the nodes
+
+for i=1:size(nodecon,1)
+    if (nodecon(i,end) > 1)
+        nodecon(i,1:nodecon(i,end))=sort(nodecon(i,1:nodecon(i,end)));
+    end
+end
+
+end