Changeset 594

Timestamp:

Mar 26, 2013, 9:56:10 PM (12 years ago)

Author:

sommeria

Message:

bugs corrected. One step further for using compiled fcts and cluster with series

Location:

trunk/src

Files:

• civ_matlab.m (modified) (1 diff)
• mouse_down.m (modified) (1 diff)
• read_civdata.m (modified) (2 diffs)
• series.m (modified) (19 diffs)
• series/aver_stat.m (modified) (1 diff)
• series/check_data_files.m (modified) (6 diffs)

trunk/src/civ_matlab.m

@@ -105 +105 @@
     Data.VarAttribute{5}.Role='warnflag';
     
-    if 'PIV volume'
+    if strcmp(Param.ListCompareMode, 'PIV volume')
         Data.ListVarName=[Data.ListVarName 'Civ1_Z'];
         Data.Civ1_X=[];Data.Civ1_Y=[];Data.Civ1_Z=[];

trunk/src/mouse_down.m

@@ -384 +384 @@
     if isequal(get(h_edit_append,'Value'),1) && ~isempty(haxes)
         h_ListCoord=hh_geometry_calib.ListCoord; %findobj(h_geometry_calib,'Tag','ListCoord');
-        coord_value=get(hhuvmat.transform_fct,'Value');% set uvmat to pixel coordinates, run it again if not
+        coord_value=get(hhuvmat.TransformName,'Value');% set uvmat to pixel coordinates, run it again if not
         if ~(isequal(coord_value,1)||isequal(coord_value,3)); %active only with no transform or px (no phys)
-            set(hhuvmat.transform_fct,'Value',1)
-            uvmat('transform_fct_Callback',hObject,eventdata,hhuvmat); %file input with xml reading  in uvmat
+            set(hhuvmat.TransformName,'Value',1)
+            uvmat('TransformName_Callback',hObject,eventdata,hhuvmat); %file input with xml reading  in uvmat
             set(hhuvmat.CheckFixedLimits,'Value',0)% put FixedLimits option to 'off'
             set(hhuvmat.CheckFixedLimits,'BackgroundColor',[0.7 0.7 0.7])

trunk/src/read_civdata.m

@@ -96 +96 @@
             Field.SubDomain=Field.Patch1_SubDomain;
             Field.ListGlobalAttribute=[Field.ListGlobalAttribute {'SubDomain'}];
-        end     
+        end 
+        if isfield(Field,'Civ1_Dt')
         Field.Dt=Field.Civ1_Dt;
+        end
+        if isfield(Field,'Civ1_Time')
         Field.Time=Field.Civ1_Time;
+        end
     case{'civ2','filter2'}
         if isfield(Field,'Patch2_SubDomain')
@@ -104 +108 @@
             Field.ListGlobalAttribute=[Field.ListGlobalAttribute {'SubDomain'}];
         end
+        if isfield(Field,'Civ2_Dt')
         Field.Dt=Field.Civ2_Dt;
+        end
+        if isfield(Field,'Civ2_Time')
         Field.Time=Field.Civ2_Time;
+        end
 end
 Field.ListGlobalAttribute=[Field.ListGlobalAttribute {'Dt','Time'}];

trunk/src/series.m

@@ -186 +186 @@
 
 %% Adjust the GUI according to the binaries available in PARAM.xml
-path_uvmat=fileparts(which('uvmat')); %path to civ
-addpath (path_uvmat) ; %add the path to civ, (useful in case of change of working directory after civ has been s opened in the working directory)
-errormsg=[];%default error message
-xmlfile='PARAM.xml';
-if exist(xmlfile,'file')
-    try
-        t=xmltree(xmlfile);
-        sparam=convert(t);
-    catch ME
-        errormsg={' Unable to read the file PARAM.xml defining the  binaries:';ME.message};
-    end
-else
-    errormsg=[xmlfile ' not found: path to binaries undefined'];
-end
-if ~isempty(errormsg)
-    msgbox_uvmat('WARNING',errormsg);
-end
-test_batch=0;%default: ,no batch mode available
-if isfield(sparam,'BatchParam') && isfield(sparam.BatchParam,'BatchMode')
-    test_batch=strcmp(sparam.BatchParam.BatchMode,'sge'); %sge is currently the only implemented batch mod
-end
-RUNVal=get(handles.RunMode,'Value');
-if test_batch==0
-   if RUNVal>2
-       set(handles.RunMode,'Value',1)
-   end
-   set(handles.RunMode,'String',{'local';'background'})
-else
-    set(handles.RunMode,'String',{'local';'background';'cluster'})
-end
+% path_uvmat=fileparts(which('uvmat')); %path to civ
+% addpath (path_uvmat) ; %add the path to civ, (useful in case of change of working directory after civ has been s opened in the working directory)
+% errormsg=[];%default error message
+% xmlfile='PARAM.xml';
+% if exist(xmlfile,'file')
+%     try
+%         t=xmltree(xmlfile);
+%         sparam=convert(t);
+%     catch ME
+%         errormsg={' Unable to read the file PARAM.xml defining the  binaries:';ME.message};
+%     end
+% else
+%     errormsg=[xmlfile ' not found: path to binaries undefined'];
+% end
+% if ~isempty(errormsg)
+%     msgbox_uvmat('WARNING',errormsg);
+% end
+% test_batch=0;%default: ,no batch mode available
+% if isfield(sparam,'BatchParam') && isfield(sparam.BatchParam,'BatchMode')
+%     test_batch=strcmp(sparam.BatchParam.BatchMode,'sge'); %sge is currently the only implemented batch mod
+% end
+% RUNVal=get(handles.RunMode,'Value');
+% if test_batch==0
+%    if RUNVal>2
+%        set(handles.RunMode,'Value',1)
+%    end
+%    set(handles.RunMode,'String',{'local';'background'})
+% else
+%     set(handles.RunMode,'String',{'local';'background';'cluster'})
+% end
 
 %% introduce the input file name(s) if defined from input Param
@@ -1308 +1308 @@
 set(handles.RUN,'BackgroundColor',[0.831 0.816 0.784])
 drawnow
-[h_fun,Series,filexml,errormsg]=prepare_jobs(handles);
+[Series,filexml,errormsg]=prepare_jobs(handles);
 if ~isempty(errormsg)
     msgbox_uvmat('ERROR',errormsg)
@@ -1320 +1320 @@
 ActionList=get(handles.ActionName,'String');
 ActionName=ActionList{get(handles.ActionName,'Value')};
-
+path_series=fileparts(which('series'));
+
+% create the Action fct handle if RunMode option = 'local'
+if strcmp(RunMode,'local')
+    if ~isequal(ActionPath,path_series)
+        eval(['spath=which(''' ActionName ''');']) %spath = current path of the selected function ACTION
+        if ~exist(ActionPath,'dir')
+            msgbox_uvmat('ERROR',['The prescribed function path ' ActionPath ' does not exist']);
+            return
+        end
+        if ~isequal(spath,ActionPath)
+            addpath(ActionPath)% add the prescribed path if not the current one
+        end
+    end
+    eval(['h_fun=@' ActionName ';'])%create a function handle for ACTION
+    if ~isequal(ActionPath,path_series)
+        rmpath(ActionPath)% add the prescribed path if not the current one
+    end
+end
+
+%% Get RunTime code from the file PARAM.xml (needed to run compiled functions)
+errormsg='';%default error message
+xmlfile=fullfile(path_series,'PARAM.xml');
+test_batch=0;%default: ,no batch mode available
+if ~exist(xmlfile,'file')
+    [success,message]=copyfile(fullfile(path_series,'PARAM.xml.default'),xmlfile);
+end
+RunTime='';
+if strcmp(ActionExt,'.sh')
+    if exist(xmlfile,'file')
+        s=xml2struct(xmlfile);
+        if strcmp(RunMode,'cluster') && isfield(s,'BatchParam')
+            if isfield(s.BatchParam,'RunTime')
+                RunTime=s.BatchParam.RunTime;
+            end
+            if isfield(s.BatchParam,'NbCore')
+                NbCore=s.BatchParam.NbCore;
+            end
+        elseif (strcmp(RunMode,'background')||strcmp(RunMode,'local')) && isfield(s,'RunParam')
+            if isfield(s.RunParam,'RunTime')
+                RunTime=s.RunParam.RunTime;
+            end
+            if isfield(s.RunParam,'NbCore')
+                NbCore=s.RunParam.NbCore;
+            end
+        end
+    end
+    if isempty(RunTime) && strcmp(RunMode,'cluster')
+        msgbox_uvmat('ERROR','RunTime name not found in PARAM.xml, compiled version .sh cannot run on cluster')
+        return
+    end
+    Series.RunTime=RunTime;
+end
+   
 %% set nbre of processes
-if isfield(Series,'NbSlice')&&~isempty(Series.NbSlice)
-    NbProcess=Series.NbSlice;
-else
+if ~isfield(Series,'NbSlice')
+    Series.NbSlice=[];
+end
+if isempty(Series.NbSlice)
     switch RunMode
         case {'local','background'}
@@ -1338 +1392 @@
     end
     NbProcess=NbCore;% choose one process per core
-end
+else
+    NbProcess=Series.NbSlice;% the nbre of run processes is equal to the number of slices
+end
+
+%% read index ranges
+first_i=1;
+last_i=1;
+incr_i=1;
+first_j=1;
+last_j=1;
+incr_j=1;
+if isfield(Series.IndexRange,'first_i')
+    first_i=Series.IndexRange.first_i;
+    incr_i=Series.IndexRange.incr_i;
+    last_i=Series.IndexRange.last_i;
+end
+if isfield(Series.IndexRange,'first_j')
+    first_j=Series.IndexRange.first_j;
+    last_j=Series.IndexRange.last_j;
+end
+if last_i < first_i || last_j < first_j , msgbox_uvmat('ERROR','last field number must be larger than the first one'),...
+    set(handles.RUN, 'Enable','On'), set(handles.RUN,'BackgroundColor',[1 0 0]),return
+else
+    BlockLength=ceil(numel(first_i:incr_i:last_i)/NbProcess);
+end
+
+%% main processing
 for iprocess=1:NbProcess% TOD0
 %     split the input files
@@ -1346 +1426 @@
 switch RunMode
     case 'local'
-        switch ActionExt
-            case '.m'
-                for iprocess=1:NbProcess
-                    if ~isempty(filexml)
-                        t=struct2xml(Series);
-                        t=set(t,1,'name','Series');
-                        save(t,filexml);
-                    end
-                    Series.RUNHandle=handles.RUN;
-                    Series.WaitbarHandle=handles.Waitbar;
-                    Series=h_fun(Series);
+        Series.RUNHandle=handles.RUN;
+        Series.WaitbarHandle=handles.Waitbar;
+        for iprocess=1:NbProcess
+            if isempty(Series.NbSlice)
+                Series.IndexRange.first_i=first_i+(iprocess-1)*BlockLength;
+                Series.IndexRange.last_i=first_i+(iprocess)*BlockLength-1;
+            else
+                Series.IndexRange.first_i= first_i+iprocess-1;
+                Series.IndexRange.incr_i=incr_i*Series.NbSlice;
+            end
+            if ~isempty(filexml)
+                t=struct2xml(Series);
+                t=set(t,1,'name','Series');
+                if isequal(Series.IndexRange.last_i,Series.IndexRange.first_i)
+                    term_i=num2str(Series.IndexRange.first_i);
+                else
+                    term_i=[num2str(Series.IndexRange.first_i) '-' num2str(Series.IndexRange.last_i)];
                 end
-            case '.sh'
-                for iprocess=1:NbProcess
-                    CivmBin=fullfile(ActionPath,[ActionName '.sh']); %path to the source directory of uvmat
+                if isequal(Series.IndexRange.last_j,Series.IndexRange.first_j)
+                    term_j=num2str(Series.IndexRange.first_j);
+                else
+                    term_j=[num2str(Series.IndexRange.first_j) '-' num2str(Series.IndexRange.last_j)];
+                end
+                save(t,[filexml '_' term_i '_' term_j '.xml']);
+            end
+            switch ActionExt
+                case '.m'
+                    h_fun(Series);
+                case '.sh'
                     switch computer
-                        case {'PCWIN','PCWIN64'}
+                        case {'PCWIN','PCWIN64'} %Windows system
                             filename=regexprep(filename,'\\','\\\\');% add '\' so that '\' are left as characters
                             % TODO launch command in DOS
-                        case {'GLNX86','GLNXA64','MACI64'}
-                            cmd=['#!/bin/bash \n '...
-                                '#$ -cwd \n '...
-                                'hostname && date \n '...
-                                'umask 002 \n'...
-                                CivmBin ' ' Param.xml.RunTime ' ' regexprep(filename,'(.+)([/\\])(.+$)','$1$20_XML$2$3.xml') ' ' Param.OutputFile '.nc'];%allow writting access to created files for user group
+                        case {'GLNX86','GLNXA64','MACI64'}%Linux  system
+%                             cmd=['#!/bin/bash \n '...
+%                                 '#$ -cwd \n '...
+%                                 'hostname && date \n '...
+%                                 'umask 002 \n'...
+%                                 fullfile(ActionPath,[ActionName '.sh']) ' ' Series.RunTime ' ' filexml '_' term_i '_' term_j '.xml'];%allow writting access to created files for user group
+                             system([fullfile(ActionPath,[ActionName '.sh']) ' ' Series.RunTime ' ' filexml '_' term_i '_' term_j '.xml']);
                     end
-                end
+            end
         end
     case 'background'
@@ -1384 +1479 @@
             t=set(t,1,'name','Series');
             save(t,filexml);
-            path_uvmat=fileparts(which('uvmat'));
             
             filename_bat=regexprep(filexml,'.xml$','.bat');% create executable file to run program in background
@@ -1401 +1495 @@
                         '. /etc/sysprofile \n'...
                         'matlab -nodisplay -nosplash -nojvm -logfile ''' filelog ''' <<END_MATLAB \n'...
-                        'addpath(''' path_uvmat '''); \n'...
+                        'addpath(''' path_series '''); \n'...
                         'addpath(''' Series.Action.ActionPath '''); \n'...
                         '' Series.Action.ActionName  '( ''' filexml '''); \n'...
@@ -1413 +1507 @@
                 case {'PCWIN','PCWIN64'}
                     text_matlabscript=['matlab -automation -logfile ' regexprep(filelog,'\\','\\\\')...
-                        ' -r "addpath(''' regexprep(path_uvmat,'\\','\\\\') ''');'...
+                        ' -r "addpath(''' regexprep(path_series,'\\','\\\\') ''');'...
                         'addpath(''' regexprep(Series.Action.ActionPath,'\\','\\\\') ''');'...
                         '' Series.Action.ActionName  '( ''' regexprep(filexml,'\\','\\\\') ''');exit"'];
@@ -1577 +1671 @@
 % --- Main launch command, called by RUN and BATCH
 
-function [h_fun,Series,filexml,errormsg]=prepare_jobs(handles,run)
+function [Series,filexml,errormsg]=prepare_jobs(handles)
 %INPUT: 
 % handles: handles of graphic objects on the GUI series
@@ -1584 +1678 @@
 
 %------------------------------------------------------------------------
-h_fun=[];
 filexml='';
 errormsg='';
-if ~exist('run','var')
-    run=1;
-end
+
 %% Read parameters from series
 Series=read_GUI(handles.series);
@@ -1596 +1687 @@
 end
 
-%% read index ranges
-first_i=1;
-last_i=1;
-incr_i=1;
-first_j=1;
-last_j=1;
-incr_j=1;
-if isfield(Series.IndexRange,'first_i')
-    first_i=Series.IndexRange.first_i;
-    incr_i=Series.IndexRange.incr_i;
-    last_i=Series.IndexRange.last_i;
-end
-if isfield(Series.IndexRange,'first_j')
-    first_j=Series.IndexRange.first_j;
-    incr_j=Series.IndexRange.incr_j;
-    last_j=Series.IndexRange.last_j;
-end
-
-%% read input file parameters and set menus
-menu_coord_state=get(handles.TransformName,'Visible');
-if isequal(menu_coord_state,'on')
-    menu_index=get(handles.TransformName,'Value');
-    transform_list=get(handles.TransformName,'UserData');
-    Series.FieldTransform.TransformHandle=transform_list{menu_index};% transform function handles
-end
-
-if last_i < first_i || last_j < first_j , msgbox_uvmat('ERROR','last field number must be larger than the first one'),...
-    set(handles.RUN, 'Enable','On'), set(handles.RUN,'BackgroundColor',[1 0 0]),return,end;
 
 %% projection object
-if isfield(Series,'CheckObject')
-    if Series.CheckObject
-        hset_object=findobj(allchild(0),'tag','set_object');
-        Series.ProjObject=read_GUI(hset_object);
-        CheckObject_Callback([], [], handles)
-    end
-else
-    Series.CheckObject=0;
-end
+% if isfield(Series,'CheckObject')
+%     if Series.CheckObject
+%         hset_object=findobj(allchild(0),'tag','set_object');
+%         Series.ProjObject=read_GUI(hset_object);
+%         CheckObject_Callback([], [], handles)
+%     end
+% else
+%     Series.CheckObject=0;
+% end
 
 %% get_field GUI
@@ -1644 +1707 @@
 end
 
-if ~run
-    return
-end
-
 %% defining the ActionName function handle
-list_action=get(handles.ActionName,'String');% list menu action
-index=get(handles.ActionName,'Value');
-action= list_action{index}; % selected string
-%Series.hseries=handles.series; % handles to the series GUI
-path_series=which('series');
-ActionPathList=get(handles.ActionName,'UserData');
-ActionPath=ActionPathList{index}; %path stored for the function ACTION
-if ~isequal(ActionPath,path_series)
-    eval(['spath=which(''' action ''');']) %spath = current path of the selected function ACTION
-    if ~exist(ActionPath,'dir')
-        errormsg=['The prescribed function path ' ActionPath ' does not exist'];
-        return
-    end
-    if ~isequal(spath,ActionPath)
-        addpath(ActionPath)% add the prescribed path if not the current one
-    end
-end
-eval(['h_fun=@' action ';'])%create a function handle for ACTION
-if ~isequal(ActionPath,path_series)
-        rmpath(ActionPath)% add the prescribed path if not the current one    
-end
+% list_action=get(handles.ActionName,'String');% list menu action
+% index=get(handles.ActionName,'Value');
+% action= list_action{index}; % selected string
+% %Series.hseries=handles.series; % handles to the series GUI
+% path_series=which('series');
+% ActionPathList=get(handles.ActionName,'UserData');
+% ActionPath=ActionPathList{index}; %path stored for the function ACTION
+% if ~isequal(ActionPath,path_series)
+%     eval(['spath=which(''' action ''');']) %spath = current path of the selected function ACTION
+%     if ~exist(ActionPath,'dir')
+%         errormsg=['The prescribed function path ' ActionPath ' does not exist'];
+%         return
+%     end
+%     if ~isequal(spath,ActionPath)
+%         addpath(ActionPath)% add the prescribed path if not the current one
+%     end
+% end
+% eval(['h_fun=@' action ';'])%create a function handle for ACTION
+% if ~isequal(ActionPath,path_series)
+%         rmpath(ActionPath)% add the prescribed path if not the current one    
+% end
 
 %% create the output data directory and write in it the xml file from the GUI config
@@ -1712 +1771 @@
         end
     end
-    filexml=fullfile(OutputDir,[Series.InputTable{1,3} '.xml']);% name of the parameter xml file set in this directory
+    filexml=fullfile(OutputDir,Series.InputTable{1,3});% name of the parameter xml file set in this directory
 end
 %removes redondant information
@@ -1816 +1875 @@
 
 %% check the current ActionPath to the selected function
-PathName=ActionPathList{ActionIndex};%current recorded path
-set(handles.ActionPath,'String',PathName); %show the path to the senlected function
+ActionPath=ActionPathList{ActionIndex};%current recorded path
+set(handles.ActionPath,'String',ActionPath); %show the path to the senlected function
 
 %% reinitialise the waitbar
@@ -1829 +1888 @@
 set(handles.Mask,'Visible','off')
 
-%% run the current action function to prepare the input GUI
-[h_fun,Series,tild,errormsg]=prepare_jobs(handles);
+%% create the function handle for Action
+path_series=which('series');
+if ~isequal(ActionPath,path_series)
+    eval(['spath=which(''' ActionName ''');']) %spath = current path of the selected function ACTION
+    if ~exist(ActionPath,'dir')
+        errormsg=['The prescribed function path ' ActionPath ' does not exist'];
+        return
+    end
+    if ~isequal(spath,ActionPath)
+        addpath(ActionPath)% add the prescribed path if not the current one
+    end
+end
+eval(['h_fun=@' ActionName ';'])%create a function handle for ACTION
+if ~isequal(ActionPath,path_series)
+        rmpath(ActionPath)% add the prescribed path if not the current one    
+end
+
+%% prepare the input param
+[Series,tild,errormsg]=prepare_jobs(handles);
 if ~isempty(errormsg)
     msgbox_uvmat('ERROR',errormsg)
@@ -1877 +1953 @@
         end
     end
+else
+% check index ranges
+first_i=1;last_i=1;first_j=1;last_j=1;
+if isfield(Series.IndexRange,'first_i')
+    first_i=Series.IndexRange.first_i;
+    last_i=Series.IndexRange.last_i;
+end
+if isfield(Series.IndexRange,'first_j')
+    first_j=Series.IndexRange.first_j;
+    last_j=Series.IndexRange.last_j;
+end
+if last_i < first_i || last_j < first_j , msgbox_uvmat('ERROR','last field number must be larger than the first one'),...
+    set(handles.RUN, 'Enable','On'), set(handles.RUN,'BackgroundColor',[1 0 0]),return,end;
 end
 
@@ -2248 +2337 @@
 function MenuExportConfig_Callback(hObject, eventdata, handles)
 global Series
-[tild,Series,errormsg]=prepare_jobs(handles,0);
+[tild,Series,errormsg]=prepare_jobs(handles);
 % Series=read_GUI(handles.series);
 
@@ -2492 +2581 @@
 ActionName=ActionList{get(handles.ActionName,'Value')};
 if strcmp(ActionExt,'.sh')
-    ActionFullName=fullfile(get(handles.ActionPath,'String'),[ActionName ActionExt]);
+    ActionFullName=fullfile(get(handles.ActionPath,'String'),[ActionName '.sh']);
     if ~exist(ActionFullName,'file')
         answer=msgbox_uvmat('INPUT_Y-N','compiled version has not been created: compile now?');
@@ -2502 +2591 @@
         end
     end
+    currentdir=pwd;
+    cd(get(handles.ActionPath,'String'))
+    sh_file=dir([ActionName '.sh']);
+    m_file=dir([ActionName '.m']);
+    if isfield(m_file,'datenum') && m_file.datenum>sh_file.datenum
+        set(handles.ActionExt,'BackgroundColor',[1 1 0])
+        answer=msgbox_uvmat('INPUT_Y-N',[ActionName '.sh needs to be updated: recompile now?']);
+        if strcmp(answer,'Yes')
+            compile(ActionName)
+        end
+        set(handles.ActionExt,'BackgroundColor',[1 1 1])
+    end
+    cd(currentdir)
 end
 

trunk/src/series/aver_stat.m

@@ -62 +62 @@
 
 %%%%%%%%%%%%  STANDARD PART  %%%%%%%%%%%%
+ParamOut=[];
 %% read input parameters from an xml file if input is a file name (batch mode)
 checkrun=1;

trunk/src/series/check_data_files.m

@@ -34 +34 @@
 function ParamOut=check_data_files(Param)
 
-% %% set the input elements needed on the GUI series when the action is selected in the menu ActionName
-% if ~exist('Param','var') % case with no input parameter 
-%     ParamOut={'NbViewMax';'';...% max nbre of input file series (default='' , no limitation)
-%         'AllowInputSort';'off';...% allow alphabetic sorting of the list of input files (options 'off'/'on', 'off' by default)
-%         'NbSlice';'on'; ...%nbre of slices ('off' by default)
-%         'VelType';'off';...% menu for selecting the velocity type (options 'off'/'one'/'two',  'off' by default)
-%         'FieldName';'off';...% menu for selecting the field (s) in the input file(options 'off'/'one'/'two', 'off' by default)
-%         'FieldTransform'; 'off';...%can use a transform function
-%         'ProjObject';'off';...%can use projection object(option 'off'/'on',
-%         'Mask';'off';...%can use mask option   (option 'off'/'on', 'off' by default)
-%         'OutputDirExt';'';...%set the output dir extension
-%                ''};
-%         return
-% end
 %% input preparation mode (no RUN)
 if isstruct(Param) && isequal(Param.Action.RUN,0)
-    ParamOut.AllowInputSort='off';...% allow alphabetic sorting of the list of input file SubDir (options 'off'/'on', 'off' by default)
-    ParamOut.WholeIndexRange='off';...% prescribes the file index ranges from min to max (options 'off'/'on', 'off' by default)
-    ParamOut.NbSlice='on'; ...%nbre of slices ('off' by default)
-    ParamOut.VelType='off';...% menu for selecting the velocity type (options 'off'/'one'/'two',  'off' by default)
-    ParamOut.FieldName='off';...% menu for selecting the field (s) in the input file(options 'off'/'one'/'two', 'off' by default)
-    ParamOut.FieldTransform = 'off';...%can use a transform function
-    ParamOut.ProjObject='off';...%can use projection object(option 'off'/'on',
-    ParamOut.Mask='off';...%can use mask option   (option 'off'/'on', 'off' by default)
+    ParamOut.AllowInputSort='off';% allow alphabetic sorting of the list of input file SubDir (options 'off'/'on', 'off' by default)
+    ParamOut.WholeIndexRange='off';% prescribes the file index ranges from min to max (options 'off'/'on', 'off' by default)
+    ParamOut.NbSlice='on';%nbre of slices ('off' by default)
+    ParamOut.VelType='off';% menu for selecting the velocity type (options 'off'/'one'/'two',  'off' by default)
+    ParamOut.FieldName='off';% menu for selecting the field (s) in the input file(options 'off'/'one'/'two', 'off' by default)
+    ParamOut.FieldTransform = 'off';%can use a transform function
+    ParamOut.ProjObject='off';%can use projection object(option 'off'/'on',
+    ParamOut.Mask='off';%can use mask option   (option 'off'/'on', 'off' by default)
     ParamOut.OutputDirExt='';%set the output dir extension (blank=no output dir)
 return
 end
 %%%%%%%%%%%%  STANDARD PART  %%%%%%%%%%%%
-%% select different modes,  RUN, parameter input, BATCH
-% BATCH  case: read the xml file for batch case
+
+%% read input parameters from an xml file if input is a file name (batch mode)
+checkrun=1;
 if ischar(Param)
-        Param=xml2struct(Param);
-        checkrun=0;
-% RUN case: parameters introduced as the input structure Param
-else
-    if isfield(Param,'Specific')&& strcmp(Param.Specific,'?')
-        checkrun=1;% will only search interactive input parameters (preparation of BATCH mode)
-    else
-        checkrun=2; % indicate the RUN option is used
-    end
-    hseries=guidata(Param.hseries);%handles of the GUI series
+    Param=xml2struct(Param);% read Param as input file (batch case)
+    checkrun=0;
 end
-ParamOut=Param; %default output
-%OutputDir=[Param.OutputSubDir Param.OutputDirExt];NO OUTPUT FILE
-    
+
 %% root input file(s) and type
 RootPath=Param.InputTable(:,1);
@@ -126 +103 @@
         Tabchar{4}=FileInfo.ImageType;
         Tabchar{5}=['  compression' FileInfo.VideoCompression];
-        Tabchar{6}=[ 'quality ' num2str(FileInfo.Quality)];   
+        Tabchar{6}=[ 'quality ' num2str(FileInfo.Quality)];
     else
         Tabchar={};
@@ -135 +112 @@
             datnum=zeros(1,nbfield_j);
             for ifile=1:nbfield_i
-                stopstate=get(hseries.RUN,'BusyAction');
+                if checkrun
+                    stopstate=get(Param.RUNHandle,'BusyAction');
+                    update_waitbar(Param.WaitbarHandle,ifile/nbfield_i)
+                else
+                    stopstate='queue';
+                end
                 if isequal(stopstate,'queue')% enable STOP command
-                    update_waitbar(hseries.Waitbar,ifile/nbfield_i)         
                     file=filecell{iview,index_slice(ifile)};
                     [Path,Name,ext]=fileparts(file);
@@ -148 +129 @@
                             datnum(ifile)=datfile.datenum;
                             filefound(ifile)={datfile.name};
-                        end                     
+                        end
                         lastfield='';
                         [FileType{iview},FileInfo,Object]=get_file_type(file);
                         if strcmp(FileType{iview},'civx')||strcmp(FileType{iview},'civdata')
                             if isfield(FileInfo,'CivStage')
-                            liststage={'civ1','fix1','patch1','civ2','fix2','patch2'};
-                            lastfield=liststage{FileInfo.CivStage};
+                                liststage={'civ1','fix1','patch1','civ2','fix2','patch2'};
+                                lastfield=liststage{FileInfo.CivStage};
                             end
                         end
-                        lastfield=[FileType{iview} ', ' lastfield];                   
+                        lastfield=[FileType{iview} ', ' lastfield];
                     end
                     Tabchar(1,i_slice)={['slice #' num2str(i_slice)]};
@@ -168 +149 @@
                 message=['no set of ' num2str(NbSlice) ' (NbSlices) files found'];
             else
-                 message='no file found';
+                message='no file found';
             end
         else
@@ -177 +158 @@
             message={['oldest modification:  ' filefound{ind} ' : ' datestr(first)];...
                 ['latest modification:  ' filefound{indlast} ' : ' datestr(last)]};
-        end 
+        end
         if ~isempty(Tabchar)
-          Tabchar=reshape(Tabchar,NbSlice*(nbfield_i+1),1);
+            Tabchar=reshape(Tabchar,NbSlice*(nbfield_i+1),1);
         end
     end