Changeset 6473

Timestamp:

11/02/10 11:40:50 (14 years ago)

Author:

jschierm

Message:

Modified kml_part_elems.m to show actual elements in partition rather than entire group.

Location:

issm/trunk/src/m/kml

Files:

• kml_part_edges.m (modified) (1 diff)
• kml_part_elems.m (modified) (4 diffs)
• kml_partitions.m (modified) (1 diff)

issm/trunk/src/m/kml/kml_part_edges.m

@@ -128 +128 @@
 
         [icol,irow]=find(epart'==k);
-        if isempty(icol)
+        if isempty(irow)
             continue;
         end

issm/trunk/src/m/kml/kml_part_elems.m

@@ -123 +123 @@
     for k=1:md.npart
 
-%  for each partition, find all the included elements and determine the
-%  perimeter (including those shared by another partition)
+%  for each partition, find all the included elements
 
         [icol,irow]=find(epart'==k);
-        if isempty(icol)
+        if isempty(irow)
             continue;
         end
         irow=unique(irow);
         elem=md.elements(irow,:);
-        nodecon=nodeconnectivity(elem,md.numberofgrids);
-        [edgeper,elemper,iloop]=edgeperimeter(elem,nodecon);
-        iloop(end+1)=size(edgeper,1)+1;
 
 %  determine the data to be used for the colors (if any)
@@ -155 +151 @@
 
         kplace=kml_placemark();
-        if (length(iloop)-1 > 1)
-            kplace.name      =sprintf('Partition %d (%d loops)',k,length(iloop)-1);
-        else
-            kplace.name      =sprintf('Partition %d',k);
-        end
+        kplace.name      =sprintf('Partition %d (%d elements)',k,size(elem,1));
         kplace.visibility=1;
         if exist('pdata','var')
@@ -176 +168 @@
 
         kmgeom=kml_multigeometry();
-        kmgeom.geometry  ={repmat(kml_polygon(),1,length(iloop)-1)};
-
-%  loop over each loop of the perimeter for the given partition
-
-        for i=1:length(iloop)-1
+        kmgeom.geometry  ={repmat(kml_polygon(),1,size(elem,1))};
+
+%  loop over each element for the given partition
+
+        for i=1:size(elem,1)
             kpoly=kml_polygon();
             kpoly.extrude   =1;
@@ -186 +178 @@
 
             kring=kml_linearring();
-            kring.coords    =zeros(iloop(i+1)-iloop(i)+1,3);
-
-%  loop over the element edges on the loop of the partition
-
-            for j=iloop(i):iloop(i+1)-1
-                [lat,long]=mapxy(md.x(edgeper(j,1)),md.y(edgeper(j,1)),'s');
-                kring.coords(j-iloop(i)+1,:)=[long lat alt];
+            kring.coords    =zeros(size(elem,2)+1,3);
+
+%  loop over the element nodes
+
+            for j=1:size(elem,2)
+                [lat,long]=mapxy(md.x(elem(i,j)),md.y(elem(i,j)),'s');
+                kring.coords(j,:)=[long lat alt];
             end
             kring.coords(end,:)=kring.coords(1,:);

issm/trunk/src/m/kml/kml_partitions.m

@@ -129 +129 @@
 
         [icol,irow]=find(epart'==k);
-        if isempty(icol)
+        if isempty(irow)
             continue;
         end