Changeset 6312

Timestamp:

10/15/10 09:32:48 (15 years ago)

Author:

jschierm

Message:

Intermediate version of kml_mesh_write.m with overlapping partitions.

Location:

issm/trunk/src/m/kml

Files:

• edgeadjacency.m (modified) (3 diffs)
• edgeperimeter.m (added)
• kml_mesh_write.m (modified) (2 diffs)

issm/trunk/src/m/kml/edgeadjacency.m

@@ -2 +2 @@
 %  create an edge adjacency table for the elements in the model.
 %
-%  [edgeadj,edgeuns,elemuns]=edgeadjacency(elem,nodecon)
+%  [edgeadj]=edgeadjacency(elem,nodecon)
 %
 %  where the required input is:
@@ -11 +11 @@
 %    edgeadj       (numeric, edge adjacency array (elems x edges))
 %
-%  the optional output is:
-%    edgeuns       (numeric, unshared edge list (edgeuns x 2))
-%    elemuns       (numeric, unshared edge element list (edgeuns x 1))
-%
-function [edgeadj,edgeuns,elemuns]=edgeadjacency(elem,nodecon)
+function [edgeadj]=edgeadjacency(elem,nodecon)
 
 if ~nargin
@@ -49 +45 @@
 end
 
-%%  create the unshared edge list
-
-if (nargout >= 2)
-    [icol,irow]=find(edgeadj'==0);
-    edgeuns=zeros(length(irow),2);
-    if (nargout >= 3)
-        elemuns=zeros(length(irow),1);
-    end
-
-%  loop over the edges
-
-    for i=1:length(irow)
-        edgeuns(i,1)=elem(irow(i),icol(i));
-        edgeuns(i,2)=elem(irow(i),mod(icol(i),size(elem,2))+1);
-        if (nargout >= 3)
-            elemuns(i)=irow(i);
-        end
-    end
 end
-
-end

issm/trunk/src/m/kml/kml_mesh_write.m

@@ -295 +295 @@
 if (~exist('prtplt','var') || strncmpi(prtplt,'on' ,2) || strncmpi(prtplt,'y',1)) && ...
     md.npart
-    [edgeadj,edgeuns,elemuns]=edgeadjacency(md.elements,md.nodeconnectivity);
+    [edgeadj]=edgeadjacency(md.elements,md.nodeconnectivity);
+    [icol,irow]=find(edgeadj'==0);
+    edgeuns=zeros(length(irow),2);
+    for i=1:length(irow)
+        edgeuns(i,1)=md.elements(irow(i),icol(i));
+        edgeuns(i,2)=md.elements(irow(i),mod(icol(i),size(md.elements,2))+1);
+    end
     fprintf(fid,'    <Folder>\n');
     fprintf(fid,'      <name>Unshared Edges</name>\n');
@@ -327 +333 @@
 end
 
+%  write folder for partitions
+
+fprintf(fid,'    <Folder>\n');
+fprintf(fid,'      <name>Partitions</name>\n');
+fprintf(fid,'      <visibility>1</visibility>\n');
+fprintf(fid,'      <description>Partitions=%d, Nodes=%d</description>\n',...
+    md.npart,md.numberofgrids);
+
+%  write each partition loop as a polygon
+
+disp(['Writing ' num2str(md.npart) ' partitions as KML polygons.']);
+epart=md.part(md.elements)+1;
+for k=1:md.npart
+    [icol,irow]=find(epart'==k);
+    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;
+    for i=1:length(iloop)-1
+        fprintf(fid,'      <Placemark>\n');
+        if (length(iloop)-1 > 1)
+            fprintf(fid,'        <name>Partition %d, Loop %d</name>\n',k,i);
+        else
+            fprintf(fid,'        <name>Partition %d</name>\n',k);
+        end
+        fprintf(fid,'        <visibility>1</visibility>\n');
+    %     if exist('edata','var')
+    %         fprintf(fid,'        <description>Element data: %g</description>\n',edata(i));
+    %         imap = fix((edata(i)-cmin)/(cmax-cmin)*size(cmap,1))+1;
+    %         if     (imap >= 1) && (imap <= size(cmap,1))
+    %             fprintf(fid,'        <styleUrl>#MatlabColor%d</styleUrl>\n',imap);
+    %         elseif (edata(i) == cmax)
+    %             fprintf(fid,'        <styleUrl>#MatlabColor%d</styleUrl>\n',size(cmap,1));
+    %         else
+    %             fprintf(fid,'        <styleUrl>#BlackLineEmptyPoly</styleUrl>\n');
+    %         end
+    %     else
+            fprintf(fid,'        <styleUrl>#BlackLineRandomPoly</styleUrl>\n');
+    %     end
+        fprintf(fid,'        <Polygon>\n');
+        fprintf(fid,'          <extrude>1</extrude>\n');
+        fprintf(fid,'          <altitudeMode>relativeToGround</altitudeMode>\n');
+        fprintf(fid,'          <outerBoundaryIs>\n');
+        fprintf(fid,'            <LinearRing>\n');
+        fprintf(fid,'              <coordinates>\n');
+        for j=iloop(i):iloop(i+1)-1
+            [lat(j),long(j)]=mapxy(md.x(edgeper(j,1)),md.y(edgeper(j,1)),'s');
+            fprintf(fid,'                %0.16g,%0.16g,%0.16g\n',long(j),lat(j),alt);
+        end
+        fprintf(fid,'                %0.16g,%0.16g,%0.16g\n',long(iloop(i)),lat(iloop(i)),alt);
+
+        fprintf(fid,'              </coordinates>\n');
+        fprintf(fid,'            </LinearRing>\n');
+        fprintf(fid,'          </outerBoundaryIs>\n');
+        fprintf(fid,'        </Polygon>\n');
+        fprintf(fid,'      </Placemark>\n');
+    end
+end
+fprintf(fid,'    </Folder>\n');
+
 %  write trailer data