Changeset 842

Timestamp:

06/08/09 12:47:39 (15 years ago)

Author:

seroussi

Message:

optionnal subdivision parameter

File:

• issm/trunk/src/m/utils/DataProcessing/addtrack.m (modified) (5 diffs)

issm/trunk/src/m/utils/DataProcessing/addtrack.m

@@ -1 @@
-function [x_m2 y_m2 values2]=addtrack(x_m1,y_m1,values1,track_coord,track_values,distance,exponent),
+function [x_m2 y_m2 values2]=addtrack(x_m1,y_m1,values1,track_coord,track_values,distance,exponent,varargin),
 %ADDTRACK - modify a map to take values of a track
 %
@@ -15 +15 @@
 %
 %   Usage:
-%      [x_m2 y_m2 values1 values2]=addtrack(x_m1,y_m1,values1,track_coord,track_values,distance,exposant)
+%      [x_m2 y_m2 values1 values2]=addtrack(x_m1,y_m1,values1,track_coord,track_values,distance,exposant[,subdivision])
 %
 %   Example:
 %      [xnew ynew valuesnew]=addtrack(x_m,y_m,thickness,'trackcoord.exp',thickness_track,1000,2)
+%      [xnew ynew valuesnew]=addtrack(x_m,y_m,thickness,'trackcoord.exp',thickness_track,1000,2,3)
 
 %Create a new picture precise enough to be modified by the tracks
@@ -52 +53 @@
 
 %Calculate the multiplicate factor for the new values:
-%mult=round((x_m1(2)-x_m1(1))/(1*dist_av));
-mult=3;
+if nargin==7,
+        mult=round((x_m1(2)-x_m1(1))/(1*dist_av));
+else
+        mult=varargin{1};
+end
 
 %Plug the values in the new multiplied matrix
@@ -102 +106 @@
 %Loop over the points of the track
 nel=size(track_values,1);
-fprintf('%s','      track processing progress:  0.00 %');
+fprintf('%s','      track processing progress:   0.00 %');
 for i=1:nel;
         if mod(i,1000)==0,
                 fprintf('\b\b\b\b\b\b\b')
-                fprintf('%4.2f%s',i/nel*100,' %');
+                fprintf('%5.2f%s',i/nel*100,' %');
         end
 
@@ -136 +140 @@
 if nel>1000,
         fprintf('\b\b\b\b\b\b\b\b')
-        fprintf('%4.1f%s\n',100,' %');
+        fprintf('%4.2f%s\n',100,' %');
 end