Changeset 26219

Timestamp:

04/29/21 11:50:05 (4 years ago)

Author:

Eric.Larour

Message:

CHG: new contours rendering routine.

File:

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

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

@@ -289 +289 @@
                 shpdisp(filenamei,1,stylei,linewidthi,getfieldvalue(options,'unit',1));
         end
+end
+if exist(options,'contours'),
+
+        hold on;
+        contours=getfieldvalue(options,'contours');
+        style=getfieldvalue(options,'contourstyle',{'-'});
+        linewidth=getfieldvalue(options,'linewidth',{1});
+        color=getfieldvalue(options,'contourcolor',{'r'});
+        contourheight=getfieldvalue(options,'contourheight',1); 
+
+        radius=md.solidearth.planetradius;
+        ratio=1+(contourheight*1000/radius);
+
+        
+        if ~isa(contours,'cell'),
+                contours={contours};
+        end
+        nc=length(contours);
+        if ~isa(style,'cell'), error('contour style should be a cell array'); end
+        if ~isa(linewidth,'cell'), error('contour line width should be a cell array'); end
+        if ~isa(color,'cell'), error('contour color should be a cell array'); end
+
+        for i=1:length(contours),
+                ci=contours{i};
+                if length(style)==1, sti=style{1}; else sti=style{i}; end
+                if length(color)==1, coli=color{1}; else coli=color{i}; end
+                if length(linewidth)==1, li=linewidth{1}; else li=linewidth{i}; end
+
+                for j=1:length(ci),
+                        cijx=ci(j).x*ratio;
+                        cijy=ci(j).y*ratio;
+                        cijz=ci(j).z*ratio;
+
+                        plot3(cijx,cijy,cijz,'LineWidth',li,'LineStyle',sti,'Color',coli);
+                end
+        end
+
 end