Changeset 7342 for issm/trunk/src/m/kml/kml_unsh_edges.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_unsh_edges.m (modified) (4 diffs)

issm/trunk/src/m/kml/kml_unsh_edges.m

@@ -11 +11 @@
 %
 %  and the optional input is:
+%    latsgn        (numeric, +1/-1 for north/south latitude)
 %    alt           (numeric, altitude for polygons, default 10000)
 %    prtplt        (char, 'off'/'no' for partition edge plot)
@@ -40 +41 @@
     if ischar(varargin{iarg})
         if ~isempty(strmatch(varargin{iarg},...
-                {'alt','prtplt'},...
+                {'latsgn','alt','prtplt'},...
                 'exact'))
             eval([varargin{iarg} '=varargin{iarg+1};']);
@@ -51 +52 @@
     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
 
@@ -91 +104 @@
 
         for j=1:2
-            [lat,long]=xy2ll(md.x(edgeuns(i,j)),md.y(edgeuns(i,j)),-1);
-            kline.coords(j,:)=[long lat alt];
+            kline.coords(j,:)=[md.long(edgeuns(i,j)) md.lat(edgeuns(i,j)) alt];
         end