Changeset 6785


Ignore:
Timestamp:
12/22/10 09:29:16 (15 years ago)
Author:
Eric.Larour
Message:

bypassing mapxy if lat long are found

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/kml/kml_mesh_elem.m

    r6461 r6785  
    139139    kring.coords    =zeros(size(md.elements,2)+1,3);
    140140
    141     for j=1:size(md.elements,2)
    142         [lat,long]=mapxy(md.x(md.elements(i,j)),md.y(md.elements(i,j)),'s');
    143         kring.coords(j,:)=[long lat alt];
    144     end
     141        if ~isempty(md.lat) & ~isempty(md.long),
     142                for j=1:size(md.elements,2)
     143                        lat=md.lat(md.elements(i,j)); long=md.long(md.elements(i,j));
     144                        kring.coords(j,:)=[long lat alt];
     145                end
     146        else
     147                for j=1:size(md.elements,2)
     148                        [lat,long]=mapxy(md.x(md.elements(i,j)),md.y(md.elements(i,j)),'s');
     149                        kring.coords(j,:)=[long lat alt];
     150                end
     151        end
    145152    kring.coords(end,:)=kring.coords(1,:);
    146153
Note: See TracChangeset for help on using the changeset viewer.