Changeset 7342 for issm/trunk/src/m/kml/kml_mesh_elem.m

Timestamp:

02/04/11 14:43:13 (14 years ago)

Author:

jschierm

Message:

kml: Converted all utilities to use lat/long fields in model rather than x/y, converting if necessary.

File:

• issm/trunk/src/m/kml/kml_mesh_elem.m (modified) (4 diffs)

issm/trunk/src/m/kml/kml_mesh_elem.m

@@ -11 +11 @@
 %
 %  and the optional input is:
+%    latsgn        (numeric, +1/-1 for north/south latitude)
 %    data          (numeric, element or nodal results data)
 %    alt           (numeric, altitude for polygons, default 10000)
@@ -43 +44 @@
     if ischar(varargin{iarg})
         if ~isempty(strmatch(varargin{iarg},...
-                {'data','alt',...
+                {'latsgn','data','alt',...
                  'cmin','cmax','cmap'},...
                 'exact'))
@@ -58 +59 @@
     end
     iarg=iarg+2;
+end
+
+if isempty(md.lat)  || ((numel(md.lat) == 1)  && isnan(md.lat)) || ...
+   isempty(md.long) || ((numel(md.long) == 1) && isnan(md.long))
+    if     ~exist('latsgn')
+        error(['Missing ''latsgn'' parameter to calculate missing lat/long data.']);
+    elseif (abs(latsgn) ~= 1)
+        error(['Incorrect latsgn=' num2str(latsgn) ' parameter to calculate missing lat/long data.']);
+    else
+        display('Converting x/y data to lat/long data.');
+        [md.lat,md.long]=xy2ll(md.x,md.y,latsgn);
+    end
 end
 
@@ -143 +156 @@
 
     for j=1:size(md.elements,2)
-        [lat,long]=xy2ll(md.x(md.elements(i,j)),md.y(md.elements(i,j)),-1);
-        kring.coords(j,:)=[long lat alt];
+        kring.coords(j,:)=[md.long(md.elements(i,j)) md.lat(md.elements(i,j)) alt];
     end
     kring.coords(end,:)=kring.coords(1,:);