Changeset 25058

Timestamp:

06/18/20 00:04:23 (5 years ago)

Author:

Eric.Larour

Message:

CHG: added latlon plot for 3Dsurface meshes.

Location:

issm/trunk-jpl/src/m/plot

Files:

• plot_coastlines.m (modified) (1 diff)
• processmesh.m (modified) (1 diff)

issm/trunk-jpl/src/m/plot/plot_coastlines.m

@@ -9892 +9892 @@
         set(p,'LineWidth',getfieldvalue(options,'coast_linewidth',1));
 else
-        [x,y,z]=AboveGround(coastlat,coastlon,mesh.r(1),1000);
-        hold on, p=plot3(x,y,z,'k-'); 
+        if ~strcmpi(getfieldvalue(options,'coord','xy'),'latlon'),
+                [x,y,z]=AboveGround(coastlat,coastlon,mesh.r(1),1000);
+                hold on, p=plot3(x,y,z,'k-'); 
+        else
+                hold on, p=plot(coastlon,coastlat,'k-'); 
+        end
         set(p,'Color',getfieldvalue(options,'coast_color','k'));
         set(p,'LineWidth',getfieldvalue(options,'coast_linewidth',1));

issm/trunk-jpl/src/m/plot/processmesh.m

@@ -24 +24 @@
 if strcmp(domaintype(md.mesh),'3Dsurface'),
         [x y z elements is2d isplanet] = processmesh(md.mesh,options);
+        if strcmpi(getfieldvalue(options,'coord','xy'),'latlon'),
+                x=md.mesh.long;
+                y=md.mesh.lat;
+                elements=delaunay(x,y);
+                z=md.mesh.lat; z(:)=0;
+        end
         return;
 end