Changeset 13836

Timestamp:

10/26/12 15:53:47 (12 years ago)

Author:

Eric.Larour

Message:

NEW: curvedarrow routine + changes to add more options.

Location:

issm/trunk-jpl/src/m/plot

Files:

• applyoptions.m (modified) (1 diff)
• colormaps/getcolormap.m (modified) (1 diff)
• curvedarrow.m (added)
• plot_riftrelvel.m (modified) (4 diffs)

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

@@ -265 +265 @@
 if exist(options,'northarrow'),
         northarrow(getfieldvalue(options,'northarrow'));
+end
+
+%curved arrow
+if exist(options,'curvedarrow'),
+        curvedoptions=getfieldvalue(options,'curvedarrow');
+        curvedarrow(curvedoptions{:});
 end
 

issm/trunk-jpl/src/m/plot/colormaps/getcolormap.m

@@ -21 +21 @@
 
 if strcmpi(map,'Ala'),
-        map = jet;
+        map = jet(256);
         map = map(128:end,:);
 elseif strcmpi(map,'redblue'),

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

@@ -15 +15 @@
 end
 options=addfielddefault(options,'scaling',2);
+
+%markersize: 
+markersize=getfieldvalue(options,'markersize',1);
 
 %recover vx and vy:
@@ -37 +40 @@
 subplot(nlines,ncols,index); 
 hold on
-
-%plot mesh boundaries
-for i=1:size(md.mesh.segments,1),
-        plot(x(md.mesh.segments(i,1:2)),y(md.mesh.segments(i,1:2)),'k-');
-end
 
 %plot rifts vel
@@ -81 +79 @@
                 %if penetration is negative, plot in black, positive, plot in red;: ie: if rift is closing, black, if rift is opening, red.
                 if(penetration>0),
-                        p2=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'*'); set(p2,'Color',[140 140 140]/255);
+                        p2=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'*','MarkerSize',markersize); set(p2,'Color',[140 140 140]/255);
                         isp2=1;
                 else
-                        p1=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'k*');
+                        p1=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'k*','MarkerSize',markersize);
                         isp1=1;
                 end
@@ -90 +88 @@
 
         %point out the tips
-        h2=plot(x(md.rifts.riftstruct(i).tips(1)),y(md.rifts.riftstruct(i).tips(1)),'g*');
-        plot(x(md.rifts.riftstruct(i).tips(2)),y(md.rifts.riftstruct(i).tips(2)),'g*');
+        h2=plot(x(md.rifts.riftstruct(i).tips(1)),y(md.rifts.riftstruct(i).tips(1)),'g*','MarkerSize',markersize);
+        plot(x(md.rifts.riftstruct(i).tips(2)),y(md.rifts.riftstruct(i).tips(2)),'g*','MarkerSize',markersize);
         segments=md.rifts.riftstruct(i).segments(:,1:2);
 end