Changeset 874 for trunk/src/series/time_series.m

Timestamp:

Feb 19, 2015, 2:17:59 AM (9 years ago)

Author:

sommeria

Message:

merge_proj corrected + various bugs

File:

• trunk/src/series/time_series.m (modified) (6 diffs)

trunk/src/series/time_series.m

@@ -92 +92 @@
         end
     end
+    % test for subtraction
+       if size(Param.InputTable,1)==2
+            answer=msgbox_uvmat('INPUT_Y-N','substract the two input file series (Yes) or concatene them (No)');
+            ParamOut.ActionInput.CheckSub=str2double(answer);
+        end
     % check the existence of the first and last file in the series
      first_j=[];
@@ -161 +166 @@
 CheckNc=cell(1,nbview);
 frame_index=cell(1,nbview);
+
 for iview=1:nbview
     if ~exist(filecell{iview,1}','file')
@@ -284 +290 @@
     VarMesh=Param.ActionInput.VarMesh;
     end
+end
+CheckSub=0;%default
+if isfield(Parma,'ActionInput') && isfield(Param.ActionInput,'CheckSub')
+CheckSub=ParamOut.ActionInput.CheckSub;
 end
 %%%%%%%%%%%%%%%% loop on field indices %%%%%%%%%%%%%%%%
@@ -312 +322 @@
         % coordinate transform (or other user defined transform)
         if ~isempty(transform_fct)
-            switch nargin(transform_fct)
-                case 4
-                    if length(Data)==2
-                        Field=transform_fct(Data{1},XmlData{1},Data{2},XmlData{2});
-                    else
+            if CheckSub
+                switch nargin(transform_fct)
+                    case 4
+                        if length(Data)==2
+                            Field=transform_fct(Data{1},XmlData{1},Data{2},XmlData{2});
+                        else
+                            Field=transform_fct(Data{1},XmlData{1});
+                        end
+                    case 3
+                        if length(Data)==2
+                            Field=transform_fct(Data{1},XmlData{1},Data{2});
+                        else
+                            Field=transform_fct(Data{1},XmlData{1});
+                        end
+                    case 2
                         Field=transform_fct(Data{1},XmlData{1});
-                    end
-                case 3
-                    if length(Data)==2
-                        Field=transform_fct(Data{1},XmlData{1},Data{2});
-                    else
-                        Field=transform_fct(Data{1},XmlData{1});
-                    end
-                case 2
-                    Field=transform_fct(Data{1},XmlData{1});
-                case 1
-                    Field=transform_fct(Data{1});
+                    case 1
+                        Field=transform_fct(Data{1});
+                end
+            else
+                if nargin(transform_fct)>=2
+                    for iview=1:nbview
+                        Data{iview}=transform_fct(Data{iview},XmlData{iview});
+                    end
+                else
+                    for iview=1:nbview
+                        Data{iview}=transform_fct(Data{iview});
+                    end
+                end
             end
         end
@@ -344 +366 @@
             end
         end
-%         nbfile=nbfile+1;
+        %         nbfile=nbfile+1;
         
         % initiate the time series at the first iteration
@@ -441 +463 @@
                 VarVal=Field.(VarName);
                 if testsum(ivar)==2% test for recorded variable
-                    if isempty(errormsg)                     
+                    if isempty(errormsg)
                         VarVal=shiftdim(VarVal,-1); %shift dimension
                         DataOut.(VarName)(index,:,:)=VarVal;%concanete the current field to the time series