Changeset 6785
- Timestamp:
- 12/22/10 09:29:16 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/kml/kml_mesh_elem.m
r6461 r6785 139 139 kring.coords =zeros(size(md.elements,2)+1,3); 140 140 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 145 152 kring.coords(end,:)=kring.coords(1,:); 146 153
Note:
See TracChangeset
for help on using the changeset viewer.