Changeset 24882

Timestamp:

05/22/20 09:25:00 (5 years ago)

Author:

Eric.Larour

Message:

CHG: extended plot coastlines to 2D plots too.

File:

• issm/trunk-jpl/src/m/plot/plot_coastlines.m (modified) (1 diff)

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

@@ -9881 +9881 @@
 end
 
-[x,y,z]=AboveGround(coastlat,coastlon,mesh.r(1),1000);
-hold on, p=plot3(x,y,z,'k-'); 
-set(p,'Color',getfieldvalue(options,'coast_color','k'));
-set(p,'LineWidth',getfieldvalue(options,'coast_linewidth',1));
+%check on the type of mesh:
+if strcmpi(class(mesh),'mesh2d'),
+        %project: 
+        xl=getfieldvalue(options,'xlim',xlim);
+        yl=getfieldvalue(options,'ylim',ylim);
+        [x,y]=gdaltransform(coastlon,coastlat,'EPSG:4326',mesh.proj);
+        hold on, p=plot(x,y,'k-');
+        xlim(xl); ylim(yl); 
+        set(p,'Color',getfieldvalue(options,'coast_color','k'));
+        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-'); 
+        set(p,'Color',getfieldvalue(options,'coast_color','k'));
+        set(p,'LineWidth',getfieldvalue(options,'coast_linewidth',1));
+end