Changeset 6299

Timestamp:

10/14/10 15:02:17 (15 years ago)

Author:

jschierm

Message:

Addition of unshared edges to kml_mesh_write.m.

Location:

issm/trunk/src/m/kml

Files:

• edgeadjacency.m (added)
• kml_mesh_write.m (modified) (4 diffs)

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

@@ -12 +12 @@
 %
 %  and the optional input is:
-%    data          (numeric, results data)
+%    data          (numeric, element or nodal results data)
 %    alt           (numeric, altitude for polygons, default 10000)
 %    lwidth        (numeric, line width in pixels, default 1)
 %    popac         (numeric, polygon opacity, default 0.50)
-%    cmin          (numeric, minimum of color scale)
-%    cmax          (numeric, maximum of color scale)
+%    cmin          (numeric, minimum of color map)
+%    cmax          (numeric, maximum of color map)
 %    cmap          (char or numeric, colormap definition)
+%    prtplt        (char, 'off'/'no' for partition segment plot)
 %
 function []=kml_mesh_write(varargin)
@@ -256 +257 @@
 %  write folder for partition segments
 
-if md.npart
+if (~exist('prtplt','var') || strncmpi(prtplt,'on' ,2) || strncmpi(prtplt,'y',1)) && ...
+    md.npart
     [xseg,yseg]=flagedges(md.elements,md.x,md.y,md.part);
     fprintf(fid,'    <Folder>\n');
@@ -264 +266 @@
         md.npart,size(xseg,1));
 
-%  write each element as a polygon
+%  write each segment as a linestring
 
     disp(['Writing ' num2str(size(xseg,1)) ' partition segments as KML linestrings.']);
@@ -289 +291 @@
 end
 
+%  write folder for partition edges
+
+if (~exist('prtplt','var') || strncmpi(prtplt,'on' ,2) || strncmpi(prtplt,'y',1)) && ...
+    md.npart
+    [edgeadj,edgeuns,elemuns]=edgeadjacency(md.elements,md.nodeconnectivity);
+    fprintf(fid,'    <Folder>\n');
+    fprintf(fid,'      <name>Partition Edges</name>\n');
+    fprintf(fid,'      <visibility>1</visibility>\n');
+    fprintf(fid,'      <description>Partitions=%d, Edges=%d</description>\n',...
+        md.npart,size(edgeuns,1));
+
+%  write each partition edge as a linestring
+
+    disp(['Writing ' num2str(size(edgeuns,1)) ' partition edges as KML linestrings.']);
+    for i=1:size(edgeuns,1)
+        fprintf(fid,'      <Placemark>\n');
+        fprintf(fid,'        <name>Edge %d</name>\n',i);
+        fprintf(fid,'        <visibility>1</visibility>\n');
+        fprintf(fid,'        <styleUrl>#RedLineRedPoly</styleUrl>\n');
+        fprintf(fid,'        <LineString>\n');
+        fprintf(fid,'          <extrude>1</extrude>\n');
+        fprintf(fid,'          <tessellate>1</tessellate>\n');
+        fprintf(fid,'          <altitudeMode>relativeToGround</altitudeMode>\n');
+        fprintf(fid,'          <coordinates>\n');
+        for j=1:2
+            [lat(j),long(j)]=mapxy(md.x(edgeuns(i,j)),md.y(edgeuns(i,j)),'s');
+            fprintf(fid,'            %0.16g,%0.16g,%0.16g\n',long(j),lat(j),alt);
+        end
+
+        fprintf(fid,'          </coordinates>\n');
+        fprintf(fid,'        </LineString>\n');
+        fprintf(fid,'      </Placemark>\n');
+    end
+    fprintf(fid,'    </Folder>\n');
+end
+
 %  write trailer data