Changeset 180

Timestamp:

Jan 18, 2011, 10:50:40 PM (13 years ago)

Author:

sommeria

Message:

rationalisation of uvmat, introduction of the new function read_field, links with get_field, several bug repairs

Location:

trunk/src

Files:

• calc_field.m (modified) (1 diff)
• civ.m (modified) (3 diffs)
• editxml.m (modified) (2 diffs)
• get_field.m (modified) (20 diffs)
• proj_field.m (modified) (1 diff)
• read_civxdata.m (modified) (3 diffs)
• set_object.m (modified) (3 diffs)
• sub_field.m (modified) (1 diff)
• update_obj.m (modified) (1 diff)
• uvmat.m (modified) (52 diffs)
• view_field.m (modified) (1 diff)

trunk/src/calc_field.m

@@ -67 +67 @@
     units_cell={};
     for ilist=1:length(FieldName)
-        [VarName,Value,Role,units]=feval(FieldName{ilist},DataIn);%calculate field with appropriate function named FieldName{ilist}
-        ListVarName=[ListVarName VarName];
-        ValueList=[ValueList Value];
-        RoleList=[RoleList Role];
-        units_cell=[units_cell units];
+        if ~isempty(FieldName{ilist})
+            [VarName,Value,Role,units]=feval(FieldName{ilist},DataIn);%calculate field with appropriate function named FieldName{ilist}
+            ListVarName=[ListVarName VarName];
+            ValueList=[ValueList Value];
+            RoleList=[RoleList Role];
+            units_cell=[units_cell units];
+        end
     end
     %erase previous data (except coordinates)

trunk/src/civ.m

@@ -1598 +1598 @@
 box_test(5)=get(handles.FIX2,'Value');
 box_test(6)=get(handles.PATCH2,'Value');
-index=find(box_test==1);
-if isempty(index)
+index_first=find(box_test==1,1);
+if isempty(index_first)
     msgbox_uvmat('ERROR','no selected operation')
     return
 end
-index_first=min(index);
-index_last=max(index);
+index_last=find(box_test==1,1,'last');
 box_used=box_test(index_first : index_last);
 [box_missing,ind_missing]=min(box_used);
-if isequal(box_missing,0)
+if isequal(box_missing,0); %there is a missing step in the sequence of operations
     msgbox_uvmat('ERROR',['missing' cell2mat(operations(ind_missing))]);
     return
@@ -1640 +1639 @@
 %% reinitialise status callback 
 if isfield(handles,'status')
-set(handles.status,'Value',0);%suppress status display
-status_Callback(hObject, eventdata, handles)
+    set(handles.status,'Value',0);%suppress status display
+    status_Callback(hObject, eventdata, handles)
 end
 
 %% set the list of files and check them
 display('checking the files...')
-%compare=get(handles.compare,'Value');%test for usual PIV (compare=1) or displacement (=2) or stereo PIV (=3)
 [filecell,num1_civ1,num2_civ1,num_a_civ1,num_b_civ1,num1_civ2,num2_civ2,num_a_civ2,num_b_civ2,nom_type_nc]=...
     set_civ_filenames(handles,box_test);
@@ -4623 +4621 @@
 box_test(4)=get(handles.CIV2,'Value');
 box_test(5)=get(handles.FIX2,'Value');
-box_test(6)=get(handles.PATCH2,'Value')
-find(box_test)
-option_civ=max(find(box_test))
+box_test(6)=get(handles.PATCH2,'Value');
+option_civ=find(box_test,1,'last');%last selected option (non-zero index of box_test)
 filecell=get(handles.civ,'UserData');
 if ~isfield(filecell,'nc')

trunk/src/editxml.m

@@ -330 +330 @@
     head_name=get(t,1,'name');
     %Path to shemas:
-    path_uvmat=which('editxml');% check the path detected for source file uvmat
-    path_UVMAT=fileparts(path_uvmat); %path to UVMAT
-    xmlparam=fullfile(path_UVMAT,'PARAM.xml');
+%     path_uvmat=which('editxml');% check the path detected for source file uvmat
+%     path_UVMAT=fileparts(path_uvmat); %path to UVMAT
+%     xmlparam=fullfile(path_UVMAT,'PARAM.xml');
+    xmlparam='PARAM.xml'; %will find PARAM.xml whose path is set in priority
     if exist(xmlparam,'file')
         tparam=xmltree(xmlparam);
@@ -350 +351 @@
             '*.xsd',  '.xsd files '; ...
             '*.*',  'All Files (*.*)'}, ...
-            'Pick a .xsd schema' ,schemapath); %file browser
+            'Pick a .xsd schema' ,schemafile); %file browser
             if ischar(PathName) && ischar(FileName) && exist(fullfile(PathName,FileName),'file')
                 DataIn.Schema=fullfile(PathName,FileName);

trunk/src/get_field.m

@@ -45 +45 @@
 %------------------------------------------------------------------------
 % --- Executes just before get_field is made visible.
-function get_field_OpeningFcn(hObject, eventdata, handles,filename,Field,haxes)
-%------------------------------------------------------------------------
-global nb_builtin
-
-%% look at the existing figues in the work space
+function get_field_OpeningFcn(hObject, eventdata, handles,filename,multiple)
+%------------------------------------------------------------------------
+global nb_builtin % nbre of functions to include by default in the menu of  functions called by RUN
+
+%% look at the existing figures in the work space
 browse_fig(handles.list_fig)
 
@@ -61 +61 @@
 set(hObject,'WindowButtonUpFcn',{@mouse_up_gui,handles})
 
-%% prepare the list of builtin fcts and set their paths
-menu_str={'PLOT';'FFT';'filter_band';'histogram'}; %list of functions included by default in 'get_field.m'
+%% prepare the list of RUN fcts and set their paths
+% functions included by default in 'get_field.m
+menu_str={'PLOT';'FFT';'filter_band';'histogram'}; 
 nb_builtin=numel(menu_str);
 path_uvmat=fileparts(which('uvmat'));%path of the function 'uvmat'
@@ -78 +79 @@
 rmpath(fullfile(path_uvmat,'get_field'))
 dir_perso=prefdir;
+
+% look for functions previously used (names and paths saved in the personal file uvmat_perso.mat):
 profil_perso=fullfile(dir_perso,'uvmat_perso.mat');
 if exist(profil_perso,'file')
-    % menu={'RUN';'raw2phys';'histogram';'FFT';'peaklocking'};
       h=load (profil_perso);
      if isfield(h,'get_field_fct') && iscell(h.get_field_fct)
@@ -99 +101 @@
      end
 end
+
 menu_str=menu_str(testexist==1);%=menu_str(testexist~=0)
 fct_handle=fct_handle(testexist==1);
@@ -110 +113 @@
 if exist('filename','var') && ischar(filename) %transfer input file name in slave mode
     set(handles.inputfile,'String',filename)% prefill the input file name 
-    set(handles.inputfile,'Enable','off')% desactivate the input file edit box   
-  %  set(handles.list_fig,'Value',2)% plotting axes =uvmat selected
-  %  set(handles.list_fig,'Visible','off')% 
-    set(handles.RUN,'String','REFRESH')% passive mode, get_field used to define the field for uvamt
-    set(handles.MenuOpen,'Visible','off')
-    set(handles.MenuExport,'Visible','off')
-    set(handles.MenuHelp,'Visible','off')
-    inputfile_Callback(hObject, eventdata, handles)
+    Field=nc2struct(filename);% reads the whole field
+    if isfield(Field,'Txt')
+        msgbox_uvmat('ERROR',Field.Txt)
+    else
+        set(handles.figure1,'UserData',Field);
+        Field_input(eventdata,handles,Field);
+    end
 else  %master mode
     set(handles.inputfile,'String','')   
-    set(handles.RUN,'String','RUN')%
+%     set(handles.RUN,'String','RUN')%
     % load the list of previously browsed files for the upper bar menu Open 
     dir_perso=prefdir;
@@ -143 +145 @@
     end
 end
-%transfer input field  in slave mode
-if exist('Field','var') && isstruct(Field)
-        Field_input(eventdata,handles,Field)
-%         if exist('haxes','var')
-%             'TESTget'
-%             get(haxes,'Tag')
-%             Field.PlotAxes=haxes;
-%         end
-    set(hObject,'UserData',Field);
-end
-
+
+%% remove already opened get_field GUI with name get_field
+if ~(exist('multiple','var') && isequal(multiple,1)) %set single occurrence
+    hget_field=findobj(allchild(0),'Name','get_field'); %hget_field(1)= new GUI
+    if length(hget_field)>1
+        delete(hget_field(2))
+    end
+else
+    set(hObject,'name','get_field_1')
+end
 
 
@@ -163 +164 @@
 
 %------------------------------------------------------------------------
-% --- Executes on button press in browse.
-function browse_Callback(hObject, eventdata, handles)
-%------------------------------------------------------------------------
-
-%------------------------------------------------------------------------
+% --- Executes when a new input file name is introduced.
 function inputfile_Callback(hObject, eventdata, handles)
 %------------------------------------------------------------------------
@@ -178 +175 @@
 Field_input(eventdata,handles,Field);
 end
-
-%------------------------------------------------------------------------
+huvmat=findobj(allchild(0),'tag','uvmat');
+if ~isempty(huvmat)
+    delete(huvmat)%delete uvmat for plot reinitialisation 
+end
+
+%------------------------------------------------------------------------
+% --- update the display when a new field is introduced.
 function Field_input(eventdata,handles,Field)
 %------------------------------------------------------------------------
@@ -256 +258 @@
                 set(handles.coord_y_scalar,'Value',VarType{imax}.coord_y+1)
             end
+            if ~isempty(VarType{imax}.coord_z) 
+                set(handles.coord_z_scalar,'Value',VarType{imax}.coord_z+1)
+            end
             if ~isempty(VarType{imax}.coord)
-                set(handles.coord_y_scalar,'Value',VarType{imax}.coord(1)+1)
-                if numel(VarType{imax}.coord)>=2
-                    set(handles.coord_x_scalar,'Value',VarType{imax}.coord(2)+1)
+                if numel(VarType{imax}.coord)>=maxdim-2
+                    set(handles.coord_z_scalar,'Value',VarType{imax}.coord(maxdim-2)+1)
                 end
+                if numel(VarType{imax}.coord)>=maxdim-1
+                    set(handles.coord_y_scalar,'Value',VarType{imax}.coord(maxdim-1)+1)
+                end
+                if numel(VarType{imax}.coord)>=maxdim
+                    set(handles.coord_x_scalar,'Value',VarType{imax}.coord(maxdim)+1)
+                end     
             end
         end
@@ -272 +282 @@
 function ordinate_Callback(hObject, eventdata, handles)
 %------------------------------------------------------------------------
-%update_field(hObject, eventdata, handles)
-% A REVOIR
 hselect_field=get(handles.inputfile,'parent');
 Field=get(hselect_field,'UserData');
-% xindex=get(handles.abscissa,'Value');
 list=get(handles.ordinate,'String');
 yindex=get(handles.ordinate,'Value');
@@ -298 +305 @@
 coord_x_index=coord_x_index(coord_x_index~=0);
 set(handles.abscissa,'String',[{''}; (Field.ListVarName(coord_x_index))'; (Field.ListVarName(VarIndex))'])
-% Field.VarIndex.y=yindex;
-% set(hselect_field,'UserData',Field);
-%update_UserData(handles)
 
 %------------------------------------------------------------------------
@@ -794 +798 @@
 %---------------------------------------------------------
 % --- Executes on button press in RUN.
+
 function RUN_Callback(hObject, eventdata, handles)
 %---------------------------------------------------------
@@ -799 +804 @@
 index=get(handles.list_fig,'value');
 figstring=figcell{index};
+
+% plot requested in uvmat
 if isequal(figstring,'uvmat')
+    inputfile=get(handles.inputfile,'String');
     huvmat=findobj(allchild(0),'tag','uvmat');
-    hhuvmat=guidata(huvmat);
-    uvmat('run0_Callback',hObject,eventdata,hhuvmat); % display field in uvmat
+    if isempty(huvmat)
+        inputfile=get(handles.inputfile,'String');
+        uvmat(inputfile)
+    else
+        hhuvmat=guidata(huvmat);
+        uvmat('run0_Callback',hObject,eventdata,hhuvmat); % display field in uvmat
+    end
+   
+% other kind of plot
 else
-    huvmat=findobj(allchild(0),'tag','uvmat');
-    delete(huvmat)% 
     index=get(handles.ACTION,'Value');
     list_func=get(handles.ACTION,'UserData');
@@ -1077 +1090 @@
 set(handles.attributes,'String',Tabchar);
 
-% list_var=get(handles.dimensions,'String');
-% val=get(handles.dimensions,'Value');
-
 % update dimensions;
 if isfield(Field,'ListDimName')
@@ -1109 +1119 @@
 end  
 
+%------------------------------------------------------------------------
 % --- Executes on button press in check_1Dplot.
 function check_1Dplot_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 val=get(handles.check_1Dplot,'Value');
 if isequal(val,0)
     set(handles.Panel1Dplot,'Visible','off')
-%      set(handles.scalar,'Visible','off')
-%     set(handles.ordinate,'Max',2.0)%allow multiple ordinate input option
-%     if isequal(get(handles.check_vector,'Value'),0);
-%         set(handles.coord_z_vectors_scalar,'Visible','off')
-%     end
 else
     set(handles.Panel1Dplot,'Visible','on')
-%     set(handles.scalar,'Visible','on')
-%     val=get(handles.ordinate,'Value');
-%     val=val(1);
-%     set(handles.ordinate,'Value',val);%suppress multiple ordinates
-%     set(handles.ordinate,'Max',1.0);%suppress multiple ordinate input option
-%       set(handles.coord_z_vectors_scalar,'Visible','on')
-end
-
+end
+
+%------------------------------------------------------------------------
 % --- Executes on button press in check_scalar.
 function check_scalar_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 val=get(handles.check_scalar,'Value');
 if isequal(val,0)
@@ -1138 +1141 @@
 end
 
-%---------------------------
+%------------------------------------------------------------------------
 % --- Executes on button press in check_vector.
 function check_vector_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 val=get(handles.check_vector,'Value');
 if isequal(val,0)
@@ -1148 +1152 @@
 end
 
-%-----------------------------
+%------------------------------------------------------------------------
 function mouse_up_gui(ggg,eventdata,handles)
+%------------------------------------------------------------------------
 if isequal(get(ggg,'SelectionType'),'alt') 
     message='';  
@@ -1164 +1169 @@
 end
 
-%---------------------------------------------
+%------------------------------------------------------------------------
 % --- Executes on selection change in ACTION.
+%------------------------------------------------------------------------
 function ACTION_Callback(hObject, eventdata, handles)
 global nb_builtin
@@ -1175 +1181 @@
 % add a new function to the menu
 if isequal(ACTION,'more...')
-%     pathfct=fileparts(path_get_field);
-%     browse_name=fullfile(path_get_field,'FIELD_FCT');
-%     if length(list_path)>nb_builtin
-%         browse_name=list_path{end};% initialize browser with  the path of the last introduced function
-%     end
     [FileName, PathName] = uigetfile( ...
        {'*.m', ' (*.m)';
@@ -1300 +1301 @@
      end
 end
-list=['view_field';list];
+list=['uvmat';list];
 set(menu_handle,'Value',1)
 set(menu_handle,'String',list)

trunk/src/proj_field.m

@@ -1377 +1377 @@
                     Xbound(1)=Coord{NbDim}(1)+DXinit*(min_indx-1);
                 end 
-                iz_str=':'; %default
                 if NbDim==3
                     DimCell(1)=[]; %suppress z variable

trunk/src/read_civxdata.m

@@ -61 +61 @@
 [var,role,units,vel_type_out_cell]=varcivx_generator(FieldNames,VelType);%determine the names of constants and variables to read
 [Field,vardetect,ichoice]=nc2struct(filename,var);
-if isfield(Field,'Txt')
-    return % error in file reading
-end
+% if isfield(Field,'Txt')
+%     return % error in file reading
+% end
 if vardetect(1)==0
      Field.Txt=[ 'requested field not available in ' filename '/' VelType];
-     return
 end
 var_ind=find(vardetect);
@@ -89 +88 @@
 Field.NbCoord=Field.nb_coord;
 Field.NbDim=Field.nb_dim;
-
-%determine the appropriate constant for time and dt for the PIV pair
-test_civ1=isequal(VelTypeOut,'civ1')||isequal(VelTypeOut,'interp1')||isequal(VelTypeOut,'filter1');
-test_civ2=isequal(VelTypeOut,'civ2')||isequal(VelTypeOut,'interp2')||isequal(VelTypeOut,'filter2');
-if test_civ1
-    Field.Time=double(Field.absolut_time_T0);
-    Field.dt=double(Field.dt);
-elseif test_civ2
-    Field.Time=double(Field.absolut_time_T0_2);
-    Field.dt=double(Field.dt2);
-else
-    Field.Txt='the input file is not civx';
-    display(Field.Txt)
-end
 
 % CivStage
@@ -119 +104 @@
 end 
 
+%determine the appropriate constant for time and dt for the PIV pair
+test_civ1=isequal(VelTypeOut,'civ1')||isequal(VelTypeOut,'interp1')||isequal(VelTypeOut,'filter1');
+test_civ2=isequal(VelTypeOut,'civ2')||isequal(VelTypeOut,'interp2')||isequal(VelTypeOut,'filter2');
+Field.Time=0; %default
+if test_civ1
+    if isfield(Field,'absolut_time_T0')
+        Field.Time=double(Field.absolut_time_T0);
+        Field.dt=double(Field.dt);
+    else
+       Field.Txt='the input file is not civx'; 
+       Field.CivStage=0;
+       Field.dt=0;
+    end
+elseif test_civ2
+    Field.Time=double(Field.absolut_time_T0_2);
+    Field.dt=double(Field.dt2);
+else
+    Field.Txt='the input file is not civx';
+    Field.CivStage=0;
+    Field.dt=0;
+end
+
+
+
 %% rescale fields to pixel coordinates
 if isfield(Field,'pixcmx')

trunk/src/set_object.m

@@ -426 +426 @@
         set(handles.XObject,'TooltipString',['XObject:  x coordinate of the axis origin for the ' ObjectStyle])
         set(handles.YObject,'TooltipString',['YObject:  y coordinate of the axis origin for the ' ObjectStyle])
-        set(handles.ZMax,'TooltipString',['ZMax: range of projection normal to the plane'])
+        set(handles.ZMax,'TooltipString','ZMax: range of projection normal to the plane')
         if test3D
             set(handles.Theta,'Visible','on')
@@ -733 +733 @@
 if isequal(get(hhuvmat.list_object_2,'Visible'),'on')
     IndexObj_2=get(hhuvmat.list_object_2,'Value');
-    List2=get(hhuvmat.list_object_2,'String')
+    List2=get(hhuvmat.list_object_2,'String');
     if IndexObj_2==length(List2)
         IndexObj_2=[];% '...' selected
@@ -803 +803 @@
     ObjectData.DisplayHandle_uvmat=hhuvmat.axes3;
     ObjectData.DisplayHandle_view_field=[];
-elseif ~isfield(UvData.Object{IndexObj},'DisplayHandle_uvmat')% save the previous object graph handles
+elseif isfield(UvData.Object{IndexObj},'DisplayHandle_uvmat')% save the previous object graph handles
     ObjectData.DisplayHandle_uvmat=UvData.Object{IndexObj}.DisplayHandle_uvmat;
     ObjectData.DisplayHandle_view_field=UvData.Object{IndexObj}.DisplayHandle_view_field;

trunk/src/sub_field.m

@@ -289 +289 @@
     SubData.ListVarName=[SubData.ListVarName {XName_1_1} {YName_1_1} {AName_1_1}];
     DimCell=[DimCell Field_1.VarDimName(ivar_C_1)]; %(TODO: check for dimension names)
+    if testX_1
+        for icell=1:numel(DimCell)
+            if isequal(DimCell{icell}{1},SubData.VarDimName{1}{1})
+                DimCell{icell}{1}=[DimCell{icell}{1} '_1'];
+            end
+        end
+    end
     SubData.VarDimName=[SubData.VarDimName DimCell];
     if isfield(Field_1,'VarAttribute')

trunk/src/update_obj.m

@@ -31 +31 @@
             hobject=Object_out{iobj}.DisplayHandle_uvmat;%graphic handle of object #iobj in the uvmat plot
         end
-        Object_out{iobj}.DisplayHandle_uvmat=plot_object(Object_out{iobj},Object_out{IndexObj_1},hobject,'m');%update the object representation
+        Object_out{iobj}.DisplayHandle_uvmat=plot_object(Object_out{iobj},Object_out{IndexObj_1},hobject,'m');%update the object representation of Object_out{iobj} on Object_out{IndexObj_1}
     end
 % plot view_field

trunk/src/uvmat.m

@@ -212 +212 @@
 set(hObject,'Units','Normalized')
 movegui(hObject,'center')
-UvData.PosColorbar=[0.805 0.022 0.019 0.445];
-UvData.SetObjectOrigin=[-0.05 -0.03]; %position for set_object
-UvData.SetObjectSize=[0.3 0.7];
-UvData.CalOrigin=[0.95 -0.03];%position for geometry_calib (TO IMPROVE)
-UvData.CalSize=[0.28 1];
+UvData.OpenParam.PosColorbar=[0.805 0.022 0.019 0.445];
+UvData.OpenParam.SetObjectOrigin=[-0.05 -0.03]; %position for set_object
+UvData.OpenParam.SetObjectSize=[0.3 0.7];
+UvData.OpenParam.CalOrigin=[0.95 -0.03];%position for geometry_calib (TO IMPROVE)
+UvData.OpenParam.CalSize=[0.28 1];
 
 %functions for the mouse and keyboard
@@ -235 +235 @@
     set(handles.list_object_1,'String',{'1-PLANE'});
 end
+set(handles.Fields,'Value',1)
+set(handles.Fields,'string',{''})%TODO:  PUT IN THE GUI
 
 %TRANSFORM menu: builtin fcts
@@ -305 +307 @@
 set(handles.uvmat,'UserData',UvData)
 %initiates menu of vector colors 
-list_menu=calc_field;
-%list_menu=[{'ima_cor'};{'black'};{'white'};list_menu(3:end)];
-set(handles.col_vec,'String',list_menu)
+% list_menu=calc_field;
+% %list_menu=[{'ima_cor'};{'black'};{'white'};list_menu(3:end)];
+% set(handles.col_vec,'String',list_menu)
 
 %check the path and date of modification of all functions in uvmat
@@ -346 +348 @@
         display_file_name(hObject, eventdata, handles,inputfile)
         testinputfield=1;
-%     else
-%         UvData.TestInputFile=0;
-    end
-    if ~isempty(Field)
-        set(handles.Fields,'Value',1)
-        set(handles.Fields,'String',{'get_field...'})    
-        testinputfield=1;
-        
-        % set the colorbar position on the interface:
-       % UvData.PosColorbar=[0.805 0.022 0.019 0.445];
-    elseif ischar(input)
-        scan_i_Callback(handles.scan_i, eventdata, handles);        
-    end
+    end
+%     if ~isempty(Field)
+%         menu_str=update_menu(handles.Fields,'get_field...');
+% %         set(handles.Fields,'Value',1)
+% %         set(handles.Fields,'String',{'get_field...'})    
+%         testinputfield=1;
+%     elseif ischar(input)
+%         scan_i_Callback(handles.scan_i, eventdata, handles);    
+%     end
 else
    if ishandle(handles.UVMAT_title)
@@ -391 +389 @@
             msgbox_uvmat('ERROR',errormsg)
         end
-    else
-        update_rootinfo(hObject,eventdata,handles);
     end
 end
@@ -874 +870 @@
 
 % set default options in menu 'Fields'
+
 if testima
-   set(handles.Fields,'Value',1) % set menu to 'image'
-   set(handles.Fields,'String',{'image';'get_field...';'velocity';'vort';'div';'more...'})
+%   set(handles.Fields,'Value',1) % set menu to 'image'
+%    set(handles.Fields,'String',{'image';'get_field...';'velocity';'vort';'div';'more...'})
 elseif isequal(FileExt,'.nc')||isequal(FileExt,'.cdf')
    Data=nc2struct(FileName,'ListGlobalAttribute','absolut_time_T0','civ');
-   col_vec=get(handles.col_vec,'String');
-   if ~isempty(Data.absolut_time_T0)&& ~isequal(Data.civ,0)
-       set(handles.Fields,'String',{'image';'get_field...';'velocity';'vort';'div';'more...'})
-       set(handles.Fields,'Value',3) % set menu to 'velocity'
-       col_vec{1}='ima_cor';
+   %col_vec=get(handles.col_vec,'String');
+   if ~isempty(Data.absolut_time_T0)&& ~isequal(Data.civ,0)%if the new input is Civx
+       FieldList=calc_field;
+       set(handles.Fields,'String',[{'image'};FieldList;{'get_field...'}]);%standard menu for civx data
+       set(handles.Fields,'Value',2) % set menu to 'velocity'
+       col_vec=FieldList;
+       col_vec(1)=[];%remove 'velocity' option for vector color (must be a scalar)
    else %general netcdf file (not civx)
        set(handles.Fields,'Value',1) % set menu to 'get_field...
        set(handles.Fields,'String',{'get_field...'})
-       col_vec{1}='get_field...';
+       col_vec={'get_field...'};
+       hget_field=findobj('Name','get_field');
+       if ~isempty(hget_field)%delete any existing get_field GUI for reinitialisation withthe new file series
+           hhget_field=guidata(hget_field);
+           if ~strcmp(get(hhget_field.inputfile,'String'),FileName)%delete any existing get_field GUI with file name different than the input
+                delete(hget_field)
+           end
+       end
    end
    set(handles.col_vec,'String',col_vec)
@@ -963 +969 @@
 %     search the files, recognize their type according to their name and fill the rootfile input windows
 function MenuBrowse_1_Callback(hObject, eventdata, handles)
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
 % huvmat=get(handles.run0,'parent');
 UvData=get(handles.uvmat,'UserData');
@@ -1031 +1037 @@
 % --- Open again as second field the file whose name has been recorded in MenuFile_1
 function MenuFile_1_1_Callback(hObject, eventdata, handles)
+% -----------------------------------------------------------------------
 fileinput_1=get(handles.MenuFile_1_1,'Label');
 display_file_name_1(hObject,eventdata,handles,fileinput_1)
@@ -1037 +1044 @@
 % --- Open again as second field the file whose name has been recorded in MenuFile_2
 function MenuFile_2_1_Callback(hObject, eventdata, handles)
+% -----------------------------------------------------------------------
 fileinput_1=get(handles.MenuFile_2_1,'Label');
 display_file_name_1(hObject,eventdata,handles,fileinput_1)
@@ -1043 +1051 @@
 % --- Open again as second field the file whose name has been recorded in MenuFile_3
 function MenuFile_3_1_Callback(hObject, eventdata, handles)
+% -----------------------------------------------------------------------
 fileinput_1=get(handles.MenuFile_3_1,'Label');
 display_file_name_1(hObject,eventdata,handles,fileinput_1)
@@ -1049 +1058 @@
 % --- Open again as second field the file whose name has been recorded in MenuFile_4
 function MenuFile_4_1_Callback(hObject, eventdata, handles)
+% -----------------------------------------------------------------------
 fileinput_1=get(handles.MenuFile_4_1,'Label');
 display_file_name_1(hObject,eventdata,handles,fileinput_1)
@@ -1055 +1065 @@
 % --- Open again as second field the file whose name has been recorded in MenuFile_5
 function MenuFile_5_1_Callback(hObject, eventdata, handles)
+% -----------------------------------------------------------------------
 fileinput_1=get(handles.MenuFile_5_1,'Label');
 display_file_name_1(hObject,eventdata,handles,fileinput_1)
@@ -1160 +1171 @@
 % default choice of fields
 if isequal(ext_test,'.image')
-   set(handles.Fields_1,'String',{'';'image';'get_field...';'velocity';'vort';'div';'more...'})
-   set(handles.Fields_1,'Value',2) % set menu to 'image'
+%    set(handles.Fields_1,'String',{'';'image';'get_field...';'velocity';'vort';'div';'more...'})
+%    set(handles.Fields_1,'Value',2) % set menu to 'image'
 elseif strcmp(FileExt_1,'.nc')||strcmp(FileExt_1,'.cdf')
    Data=nc2struct(fileinput_1,[]);
-   if isfield(Data,'absolut_time_T0')
-       set(handles.Fields_1,'String',{'';'image';'get_field...';'velocity';'vort';'div';'more...'})
-       set(handles.Fields_1,'Value',4) % set menu to 'velocity'
-   else   
-       set(handles.Fields_1,'Value',2) % set menu to 'get_field...'
-       set(handles.Fields_1,'String',{'';'get_field...'});
-   end
+%    if isfield(Data,'absolut_time_T0')
+%        set(handles.Fields_1,'String',{'';'image';'get_field...';'velocity';'vort';'div';'more...'})
+%        set(handles.Fields_1,'Value',4) % set menu to 'velocity'
+%    else   
+%        set(handles.Fields_1,'Value',2) % set menu to 'get_field...'
+%        set(handles.Fields_1,'String',{'';'get_field...'});
+%    end
 end  
 set(handles.SubField,'Visible','on')
@@ -1179 +1190 @@
 % --- Called by action in RootPath_1 edit box
 function RootPath_1_Callback(hObject,eventdata,handles)
+% -----------------------------------------------------------------------
 update_rootinfo_1(hObject,eventdata,handles)
 
@@ -1184 +1196 @@
 % --- Called by action in RootFile_1 edit box
 function RootFile_1_Callback(hObject, eventdata, handles)
+% -----------------------------------------------------------------------
 update_rootinfo_1(hObject,eventdata,handles)
-%-----------------------------------------------------------------------
 
 %------------------------------------------------------------------------
@@ -1191 +1203 @@
 function FileIndex_1_Callback(hObject, eventdata, handles)
 %------------------------------------------------------------------------
-
 run0_Callback(hObject, eventdata, handles)
 
@@ -1198 +1209 @@
 %     calibration from an xml file, then refresh current plots
 function update_rootinfo_1(hObject,eventdata,handles) %A REVOIR
-
+% -----------------------------------------------------------------------
 set(handles.RootPath_1,'BackgroundColor',[1 1 0])% indicate active program by yellow color
 drawnow
@@ -1472 +1483 @@
 set(handles.masklevel,'Value',z)
 
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
 % --- Executes on button press in view_xml.
 function view_xml_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 [FileName,RootPath,FileBase,FileIndices,FileExt]=read_file_boxes(handles);
 option=get(handles.view_xml,'String');
@@ -1482 +1494 @@
 end
 
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
 % --- Executes on button press in mask_test.
 function mask_test_Callback(hObject, eventdata, handles)
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
 %case of view mask selection
 if isequal(get(handles.mask_test,'Value'),1)
@@ -2061 +2073 @@
 %------------------------------------------------------------------------
 
-%initialisation
-% errormsg=[]; % default error message
+%% initialisation
 abstime=[];
 abstime_1=[];
@@ -2070 +2081 @@
 end
 UvData=get(handles.uvmat,'UserData');
-
-
 if ishandle(handles.UVMAT_title) %remove title panel on uvmat
     delete(handles.UVMAT_title)
@@ -2078 +2087 @@
 %% determine the main input file information for action
 FileType=[];%default
+if ~exist(filename,'file')
+    errormsg=['input file ' filename ' does not exist'];
+    return
+end
+%Ext=get(handles.FileExt,'String');
+NomType=get(handles.FileIndex,'UserData');
+%update the z position index
+nbslice_str=get(handles.nb_slice,'String');
+if isequal(nbslice_str,'volume')
+    z_index=num_j1;
+    set(handles.z_index,'String',num2str(z_index))
+else
+    nbslice=str2num(nbslice_str);
+    z_index=mod(num_i1-1,nbslice)+1;
+    set(handles.z_index,'String',num2str(z_index))
+end
+% refresh menu for save_mask if relevant
+masknumber=get(handles.masklevel,'String');
+if length(masknumber)>=z_index
+    set(handles.masklevel,'Value',z_index)
+end
+
+%% read the first input field if a filename has been introduced
 if ~isempty(filename)
-    if ~exist(filename,'file')
-        errormsg=['input file ' filename ' does not exist'];
-        return
-    end
+    ObjectName=filename;
+    FieldName=[];%default
+    VelType=[];%default
     Ext=get(handles.FileExt,'String');
-    NomType=get(handles.FileIndex,'UserData');
-    %update the z position index
-    nbslice_str=get(handles.nb_slice,'String');
-   % z_index=1;%default
-    if isequal(nbslice_str,'volume')
-        z_index=num_j1;
-        set(handles.z_index,'String',num2str(z_index))
-    else
-        nbslice=str2num(nbslice_str);
-        z_index=mod(num_i1-1,nbslice)+1;
-        set(handles.z_index,'String',num2str(z_index))
-    end
-    % refresh menu for save_mask if relevant
-    masknumber=get(handles.masklevel,'String');
-    if length(masknumber)>=z_index
-        set(handles.masklevel,'Value',z_index)
-    end
-
-    % determine the input file type
-    if isequal(Ext,'.nc')||isequal(Ext,'.cdf')
+    if strcmp(Ext,'.nc')||strcmp(Ext,'.cdf')
         FileType='netcdf';
+        list_fields=get(handles.Fields,'String');% list menu fields
+        index_fields=get(handles.Fields,'Value');% selected string index
+        FieldName= list_fields{index_fields}; % selected field
+        if ~strcmp(FieldName,'get_field...')% read the field names on the interface get_field...
+           VelType=setfield(handles);
+        end
+        if strcmp(FieldName,'velocity')
+            list_code=get(handles.color_code,'String');% list menu fields
+            index_code=get(handles.color_code,'Value');% selected string index
+            if  ~strcmp(list_code{index_code},'black') &&  ~strcmp(list_code{index_code},'white')
+                list_code=get(handles.col_vec,'String');% list menu fields
+                index_code=get(handles.col_vec,'Value');% selected string index
+                ParamIn.ColorVar= list_code{index_code}; % selected field
+            end
+        end
     elseif isfield(UvData,'MovieObject')
+        ObjectName=UvData.MovieObject;
         FileType='movie';
-        FieldName='image';
     elseif isequal(lower(Ext),'.avi')
         FileType='avi';
-        FieldName='image';
     elseif isequal(lower(Ext),'.vol')
         FileType='vol';
-        FieldName='image';
+        if isfield(UvData.XmlData,'Npy') && isfield(UvData.XmlData,'Npx')
+            ParamIn.Npy=UvData.XmlData.Npy;
+            ParamIn.Npx=UvData.XmlData.Npx;
+        else
+            errormsg='Npx and Npy need to be defined in the xml file for volume images .vol';
+            return
+        end
     else 
        form=imformats(Ext(2:end));
@@ -2122 +2154 @@
                FileType='image';
            end
-           FieldName='image';
        end
     end
-else
-    FileType='netcdf';
-    FieldName='get_field...';
-end
-VelType=[];%default
-if isequal(FileType,'netcdf')
-    list_fields=get(handles.Fields,'String');% list menu fields
-    index_fields=get(handles.Fields,'Value');% selected string index
-    FieldName= list_fields{index_fields}; % selected field
-    if ~isequal(FieldName,'get_field...')% read the field names on the interface get_field...
-       VelType=setfield(handles);
-    end
-end
-
-%% choose a second field if Subfield option is 'on'
+    ParamIn.FieldName=FieldName;
+    ParamIn.VelType=VelType;
+    ParamIn.GUIName='get_field';
+    [Field{1},ParamOut,errormsg] = read_field(ObjectName,FileType,ParamIn);
+    if isfield(ParamOut,'Npx')&& isfield(ParamOut,'Npy')
+        set(handles.npx,'String',num2str(ParamOut.Npx));% display image size on the interface
+        set(handles.npy,'String',num2str(ParamOut.Npy));
+    end
+end
+
+%% choose a second field if Subfield option is 'on', and if the field has changed
+VelType_1=[];%default
 FieldName_1=[];
-scal_color=[];
-VelType_1=setfield_1(handles);
-% sub_value=get(handles.SubField,'Value');
-FileType_1='none';%default
+ParamOut_1=[];
 if ~isempty(filename_1)
-    % test for a constant second field (comparison with a fixed field)
-    NomType_1=get(handles.FileIndex_1,'UserData');
-    Ext_1=get(handles.FileExt_1,'String');
-    % determine the input file type
-    if isequal(Ext_1,'.nc')||isequal(Ext_1,'.cdf')
-        FileType_1='netcdf';
-    elseif isfield(UvData,'MovieObject_1')
-        FileType_1='movie';
-        FieldName_1='image';
-    elseif isequal(lower(Ext_1),'.avi')
-        FileType='avi';
-        FieldName_1='image';
-    elseif isequal(lower(Ext_1),'.vol')
-        FileType_1='vol';
-        FieldName_1='image';
-    elseif length(Ext_1)>1
-       form=imformats(Ext_1(2:end));
-       if ~isempty(form)% if the extension corresponds to an image format recognized by Matlab
-           if isequal(NomType_1,'*');
-               FileType_1='multimage';
-           else
-               FileType_1='image';
-           end
-           FieldName_1='image';
-       end
-    end
-    if ~isequal(FieldName_1,'image')
-        list_fields=get(handles.Fields_1,'String');% list menu fields
-        index_fields=get(handles.Fields_1,'Value');% selected string index
-        FieldName_1= list_fields{index_fields}; % selected field
-        if isequal(VelType_1,'*')% free veltype choice
-            VelType_1=[];
-        elseif isequal(VelType_1,'"')% veltype the same as for the first field
-            if isempty(VelType)
-                VelType_1=[];
-            else
-                VelType_1=VelType;
-            end
-        end
-    end
-    test_keepdata_1=0;% test for keeping the previous stored data if the input files are unchanged
-    if ~isequal(NomType_1,'*')%in case of a series of files (not avi movie)
-        if isfield(UvData,'filename_1')&& isfield(UvData,'VelType_1') && isfield(UvData,'FieldName_1')
-            test_keepdata_1= isequal(filename_1,UvData.filename_1)&&...
-                isequal(VelType_1,UvData.filename_1) && isequal(FieldName_1,UvData.FieldName_1);
-        end
-    end
-    if test_keepdata_1
-        Field{2}=UvData.Field_1;   
-    elseif ~exist(filename_1,'file')
+    if ~exist(filename_1,'file')
         errormsg=['second file ' filename_1 ' does not exist'];
         return
-    end
-end
-
-%% read the input field(s)
-%read images
-if ~isempty(filename) && isequal(FieldName,'image')
-    Npz=1;%default
-     switch FileType
-        case 'movie'
-            try
-                A=read(UvData.MovieObject,num_i1);
-            catch
-                errormsg=lasterr;
+    else
+        Name=filename_1;
+        FieldName_1=[];%default
+        VelType_1=[];%default
+        Ext_1=get(handles.FileExt_1,'String');
+        NomType_1=get(handles.FileIndex_1,'UserData');
+        if isequal(Ext_1,'.nc')||isequal(Ext_1,'.cdf')
+            FileType_1='netcdf';
+        elseif isfield(UvData,'MovieObject_1')
+            Name=UvData.MovieObject_1;
+            FileType_1='movie';
+        elseif isequal(lower(Ext_1),'.avi')
+            FileType_1='avi';
+        elseif isequal(lower(Ext_1),'.vol')
+            FileType_1='vol';
+            if isfield(UvData.XmlData_1,'Npy') && isfield(UvData.XmlData_1,'Npx')
+                ParamIn.Npy=UvData.XmlData_1.Npy;
+                ParamIn.Npx=UvData.XmlData_1.Npx;
+            else
+                errormsg='Npx and Npy need to be defined in the xml file for volume images .vol';
                 return
             end
-        case 'avi'
-            try
-                mov=aviread(filename,num_i1);
-            catch
-                errormsg=lasterr;
-                return
+        else 
+           form=imformats(Ext_1(2:end));
+           if ~isempty(form)% if the extension corresponds to an image format recognized by Matlab
+               if isequal(NomType_1,'*');
+                   FileType_1='multimage';
+               else
+                   FileType_1='image';
+               end
+           end
+        end
+        if strcmp(FileType_1,'netcdf')
+            list_fields=get(handles.Fields_1,'String');% list menu fields
+            index_fields=get(handles.Fields_1,'Value');% selected string index
+            FieldName_1= list_fields{index_fields}; % selected field
+            if ~isequal(FieldName_1,'get_field...')% read the field names on the interface get_field...
+                VelType_1=setfield_1(handles);
             end
-            A=frame2im(mov(1));
-        case 'vol'
-            A=imread(filename);
-            if isfield(UvData.XmlData,'Npy')
-                Npz=size(A,1)/UvData.XmlData.Npy;
-                A=reshape(A',UvData.XmlData.Npx,UvData.XmlData.Npy,Npz);
-                A=permute(A,[3 2 1]);
+            if strcmp(VelType_1,'*')% free veltype choice
+                VelType_1=[];
+            elseif strcmp(VelType_1,'"')% veltype the same as for the first field
+                if isempty(VelType)
+                    VelType_1=[];
+                else
+                    VelType_1=VelType;
+                end
             end
-        case 'multimage'
-            A=imread(filename,num_i1);
-        case 'image'
-            A=imread(filename);
-    end 
-    npxy=size(A);
-    set(handles.npx,'String',num2str(npxy(2)));% display image size on the interface
-    set(handles.npy,'String',num2str(npxy(1)));
-    Rangx=[0.5 npxy(2)-0.5]; % coordinates of the first and last pixel centers
-    Rangy=[npxy(1)-0.5 0.5]; %
-    Field{1}.NbDim=2;%default
-    Field{1}.AName='image';
-    Field{1}.ListVarName={'AY','AX','A'}; % 
-    if ndims(A)==3 
-        if Npz==1;%color
-             Field{1}.VarDimName={'AY','AX',{'AY','AX','rgb'}}; %
-             Field{1}.AY=[npxy(1)-0.5 0.5];
-             Field{1}.AX=[0.5 npxy(2)-0.5]; % coordinates of the first and last pixel centers
+            if strcmp(FieldName_1,'velocity')
+                list_code=get(handles.color_code,'String');% list menu fields
+                index_code=get(handles.color_code,'Value');% selected string index
+                if  ~strcmp(list_code{index_code},'black') &&  ~strcmp(list_code{index_code},'white')
+                    list_code=get(handles.col_vec,'String');% list menu fields
+                    index_code=get(handles.col_vec,'Value');% selected string index
+                    ParamIn.ColorVar= list_code{index_code}; % selected field
+                end
+            end
+        end
+        test_keepdata_1=0;% test for keeping the previous stored data if the input files are unchanged
+        if ~isequal(NomType_1,'*')%in case of a series of files (not avi movie)
+            if isfield(UvData,'filename_1')&& isfield(UvData,'VelType_1') && isfield(UvData,'FieldName_1')
+                test_keepdata_1= strcmp(filename_1,UvData.filename_1) && strcmp(VelType_1,UvData.VelType_1) && strcmp(FieldName_1,UvData.FieldName_1);
+            end
+        end
+        if test_keepdata_1
+            Field{2}=UvData.Field_1;
         else
-            Field{1}.NbDim=3;
-            Field{1}.ListVarName=['AZ' Field{1}.ListVarName];
-            Field{1}.VarDimName={'AZ','AY','AX',{'AZ','AY','AX'}}; 
-            Field{1}.AZ=[npxy(1)-0.5 0.5];
-            Field{1}.AY=[npxy(2)-0.5 0.5];
-            Field{1}.AX=[0.5 npxy(3)-0.5]; % coordinates of the first and last pixel centers
-        end
-    else
-        Field{1}.VarDimName={'AY','AX',{'AY','AX'}}; %
-        Field{1}.AY=[npxy(1)-0.5 0.5];
-       Field{1}.AX=[0.5 npxy(2)-0.5]; % coordinates of the first and last pixel centers
-    end
-    Field{1}.A=A;
-    Field{1}.CoordUnit='pixel'; %used for mouse_motion
-end
-
-%read a second image
-if  ~isempty(filename_1) && ~test_keepdata_1 && isequal(FieldName_1,'image')
-    switch FileType_1
-        case 'movie'
-            A=read(UvData.MovieObject_1,num_i1);
-        case 'avi'
-            mov=aviread(filename_1,num_i1);
-            A=frame2im(mov(1));
-        case 'vol'
-            A=imread(filename_1);
-        case 'multimage'
-            A=imread(filename_1,num_i1);
-        case 'image'
-            A=imread(filename_1);
-        case 'netcdf'
-        otherwise
-            errormsg=['unknown input file type ' filename_1];
-    end 
-    npxy=size(A);
-    set(handles.npx,'String',num2str(npxy(2)));% display image size on the interface
-    set(handles.npy,'String',num2str(npxy(1)));
-    Rangx=[0.5 npxy(2)-0.5]; % coordinates of the first and last pixel centers
-    Rangy=[npxy(1)-0.5 0.5]; %
-    Field{2}.AName='image';
-    Field{2}.ListVarName={'AY','AX','A'}; % 
-    if size(A,3)==3;%color
-        Field{2}.VarDimName={'AY','AX',{'AY','AX','rgb'}}; %
-    else
-        Field{2}.VarDimName={'AY','AX',{'AY','AX'}}; %
-    end
-    Field{2}.AY=Rangy;
-    Field{2}.AX=Rangx;
-    Field{2}.A=A;
-    Field{2}.CoordType='px'; %used for mouse_motion
-    Field{2}.CoordUnit='px'; %used for move_mou
-end
-
-%read ncfile(s)
-CivStage_1=0;%default
-VelType_out_1=[];
-InputField={FieldName};
-InputField_1={FieldName_1};
-if (~isempty(filename)&& isequal(FileType,'netcdf')) || (~isempty(filename_1)&& isequal(FileType_1,'netcdf')) ;
-    %read the velocity field(s) from netcdf rootfile(s)
-    list_code=get(handles.col_vec,'String');% list menu fields
-    index_code=get(handles.col_vec,'Value');% selected string index
-    scal_color= list_code{index_code(1)}; % selected field
-    if isequal(FieldName,'velocity')&& ~isequal(scal_color,'black') && ~isequal(scal_color,'white')
-        InputField=[InputField scal_color];
-    end
-    if isequal(FieldName_1,'velocity') && ~isequal(scal_color,'black') && ~isequal(scal_color,'white')
-        InputField_1=[InputField_1 scal_color];
-    end
-    if isequal(FileType,'netcdf')  %read the first nc field
-        if isequal(FieldName,'get_field...')% read the field names on the interface get_field.
-            hget_field=findobj(allchild(0),'Name','get_field');%find the get_field... GUI
-            if isempty(hget_field)
-                hget_field= get_field(filename);%open the get_field GUI   
-              %  set(hget_field,'name','get_field')
-            elseif UvData.NewSeries% refresh the fet_field GUI for a new series
-                delete(hget_field)
-                hget_field= get_field(filename);%open the get_field GUI 
-                set(hget_field,'name','get_field')%rename  get_field GUI for the 'slave' mode
-            end
-            hhget_field=guidata(hget_field);
-            set(hhget_field.list_fig,'Value',1)
-            set(hhget_field.list_fig,'String',{'uvmat'})
-            funct_list=get(hhget_field.ACTION,'UserData');
-            funct_index=get(hhget_field.ACTION,'Value');
-            funct=funct_list{funct_index};%select  the current action in get_field, e;g. PLOT
-            Field{1}=funct(hget_field); %%activate the current action selected in get_field, e;g.read the names of the variables to plot 
-            CivStage=0;
-            VelType_out=[];    
-        else
-            [Field{1},VelType_out]=read_civxdata(filename,InputField,VelType);
-            if isfield(Field{1},'Txt')
-                errormsg=Field{1}.Txt;
-                return
-            end
-            CivStage=Field{1}.CivStage;
-            %UvData.NbDim=Field{1}.NbDim;
-        end
-    end
-    if ~isempty(filename_1) && ~test_keepdata_1 && isequal(FileType_1,'netcdf') %read the second file
-        if isequal(FieldName_1,'get_field...')% read the field names on the interface get_field.
-            hget_field_1=findobj(allchild(0),'Name','get_field_1');%find the get_field... GUI
-            if isempty(hget_field_1)
-                 hget_field_1= get_field(filename_1);%open the get_field GUI
-                 set(hget_field_1,'name','get_field_1')
-            elseif UvData.NewSeries% refresh the fet_field GUI for a new series
-                delete(hget_field_1)
-                hget_field_1= get_field(filename);%open the get_field GUI 
-                set(hget_field_1,'name','get_field_1')%rename  get_field GUI for the 'slave' mode
-            end
-            hhget_field_1=guidata(hget_field_1);%handles of GUI elements in get_field
-            set(hhget_field_1.list_fig,'Value',1)
-            set(hhget_field_1.list_fig,'String',{'uvmat'})
-            funct_list=get(hhget_field_1.ACTION,'UserData');
-            funct_index=get(hhget_field_1.ACTION,'Value');
-            funct=funct_list{funct_index};
-            Field{2}=funct(hget_field_1); %read the names of the variables to plot in the get_field GUI
-        else
-            [Field{2},VelType_out_1]=read_civxdata(filename_1,InputField_1,VelType_1);
-            CivStage_1=Field{2}.CivStage;
-        end
-        if ~isequal(FileType,'netcdf')
-            VelType_out=VelType_out_1;
-        end
-    end
-end
-
-%store the second field for possible latter use
-if numel(Field)==2
-    UvData.Field_1=Field{2};
+            ParamIn.FieldName=FieldName_1;
+            ParamIn.VelType=VelType_1;
+            ParamIn.GUIName='get_field_1';
+            [Field{2},ParamOut_1,errormsg] = read_field(Name,FileType_1,ParamIn);         
+            UvData.Field_1=Field{2}; %store the second field for possible use at next RUN
+        end
+    end
+end
+
+%% update uvmat interface
+if isfield(ParamOut,'Npx')
+    set(handles.npx,'String',num2str(ParamOut.Npx));% display image size on the interface
+    set(handles.npy,'String',num2str(ParamOut.Npy));
+elseif isfield(ParamOut_1,'Npx')
+    set(handles.npx,'String',num2str(ParamOut_1.Npx));% display image size on the interface
+    set(handles.npy,'String',num2str(ParamOut_1.Npy));
 end
 
 %% update the display buttons for the first velocity type (first menuline)
 veltype_handles=[handles.civ1 handles.interp1 handles.filter1 handles.civ2 handles.interp2 handles.filter2];
-if ~isequal(FileType,'netcdf')
-    reset_vel_type(veltype_handles)
-elseif isempty(VelType) && ~isequal(FieldName,'get_field...')
-    set_veltype_display(veltype_handles,CivStage)%update the display of available velocity types for the first field
-    if isempty(VelType_out)
+if ~isequal(FileType,'netcdf')|| isequal(FieldName,'get_field...')
+%     reset_vel_type(veltype_handles)
+    set(veltype_handles,'Visible','off')
+else% if isempty(ParamOut.VelType) && ~isequal(FieldName,'get_field...')
+    set_veltype_display(veltype_handles,ParamOut.CivStage)%update the display of available velocity types for the first field
+    if isempty(ParamOut.VelType)
         reset_vel_type(veltype_handles)
     else
-        handle1=eval(['handles.' VelType_out]);
+        handle1=eval(['handles.' ParamOut.VelType]);
         reset_vel_type(veltype_handles,handle1)
     end
 end
+field_index=strcmp(ParamOut.FieldName,ParamOut.FieldList);
+set(handles.Fields,'String',ParamOut.FieldList); %update the field menu
+set(handles.Fields,'Value',find(field_index,1))
 
 %% update the display buttons for the second velocity type (second menuline)
-veltype_handles_1=[handles.civ1_1 handles.interp1_1 handles.filter1_1 handles.civ2_1 handles.interp2_1 handles.filter2_1];
-if ~isequal(FileType_1,'netcdf')
-    reset_vel_type(veltype_handles_1)
-elseif isempty(VelType_1) && ~isequal(FieldName_1,'get_field...')
-    set_veltype_display(veltype_handles_1,CivStage_1)%update the display of available velocity types for the first field
-    if isempty(VelType_out_1)
-        reset_vel_type(veltype_handles_1)
-    else
-        handle1=eval(['handles.' VelType_out_1 '_1']);
-        reset_vel_type(veltype_handles_1,handle1)
-    end
+if ~isempty(filename_1)
+    veltype_handles_1=[handles.civ1_1 handles.interp1_1 handles.filter1_1 handles.civ2_1 handles.interp2_1 handles.filter2_1];
+    if ~isequal(FileType_1,'netcdf')|| isequal(FieldName_1,'get_field...')
+        set(veltype_handles_1,'Visible','off')
+        %reset_vel_type(veltype_handles_1)
+    else %if isempty(VelType_1) && ~isequal(FieldName_1,'get_field...')
+        set_veltype_display(veltype_handles_1,ParamOut_1.CivStage)%update the display of available velocity types for the first field
+        if isempty(ParamOut_1.VelType)
+            reset_vel_type(veltype_handles_1)
+        else
+            handle1=eval(['handles.' ParamOut_1.VelType '_1']);
+            reset_vel_type(veltype_handles_1,handle1)
+        end
+    end
+    field_index=strcmp(ParamOut_1.FieldName,ParamOut_1.FieldList);
+    set(handles.Fields_1,'String',ParamOut_1.FieldList); %update the field menu
+    set(handles.Fields_1,'Value',find(field_index,1))
 end
 
@@ -2432 +2322 @@
 end           
 
-%multislice case
-% if ~isempty(filename) &&(~isfield(UvData,'NbDim') || isequal(UvData.NbDim,2))&&...%2D case
-%       isfield(UvData,'XmlData') && isfield(UvData.XmlData,'GeometryCalib')&& isfield(UvData.XmlData.GeometryCalib,'SliceCoord')
-%        siz=size(UvData.XmlData.GeometryCalib.SliceCoord);
-% end
-
-%store the current open names, fields and vel types in uvmat interface 
-UvData.filename=filename;
+%% store the current open names, fields and vel types in uvmat interface 
 UvData.filename_1=filename_1;
-UvData.VelType=VelType;
-UvData.VelType_1=VelType_1;
-UvData.FieldName=FieldName;
-UvData.FieldName_1=FieldName_1;
-if ~isempty(scal_color)
-    UvData.CName=scal_color;
-end
-
-%% coordinate transform or user fct
+UvData.VelType_1=[];%default
+UvData.FieldName_1=[];
+if isfield(ParamOut_1,VelType)
+    UvData.VelType_1=ParamOut_1.VelType;
+end
+if isfield(ParamOut_1,FieldName)
+    UvData.FieldName_1=ParamOut_1.FieldName;
+end
+
+%% apply coordinate transform or other user fct
 XmlData=[];%default
 if isfield(UvData,'XmlData')%use geometry calib recorded from the ImaDoc xml file as first priority
@@ -2480 +2364 @@
 
 %% calculate scalar
-if isequal(FileType,'netcdf') && ~isequal(FieldName,'get_field...')%
-    Field{1}=calc_field(InputField,Field{1});
-end
-if length(Field)==2 && ~test_keepdata_1 && isequal(FileType_1,'netcdf') && ~isequal(FieldName_1,'get_field...')
-    Field{2}=calc_field(InputField_1,Field{2});
+if isequal(FileType,'netcdf') && ~isequal(ParamOut.CivStage,0)%&&~isempty(FieldName)%
+    Field{1}=calc_field([{ParamOut.FieldName} {ParamOut.ColorVar}],Field{1});
+end
+if length(Field)==2 && ~test_keepdata_1 && isequal(FileType_1,'netcdf') && ~isequal(ParamOut_1.FieldName,'get_field...')%&&~isempty(FieldName_1)
+    Field{2}=calc_field([{ParamOut_1.FieldName} {ParamOut_1.ColorVar}],Field{2});
 end
 
@@ -2494 +2378 @@
 end
 
-%% test 3D , default projection menuplane and typical mesh (needed to menuopen set_object)
+%% get bounds and mesh (needed for mouse action and to open set_object)
 test_x=0;
 test_z=0;% test for unstructured z coordinate
 [UvData.Field,errormsg]=check_field_structure(UvData.Field);
 if ~isempty(errormsg)
-    errormsg=['error in uvmat/run0_Callback/check_field_structure: ' errormsg];
+    errormsg=['error in uvmat/refresh_field/check_field_structure: ' errormsg];
     return
 end
-[CellVarIndex,NbDim,VarType]=find_field_indices(UvData.Field);
+[CellVarIndex,NbDim,VarType,errormsg]=find_field_indices(UvData.Field);
+if ~isempty(errormsg)
+    errormsg=['error in uvmat/refresh_field/find_field_indices: ' errormsg];
+    return
+end
 [NbDim,imax]=max(NbDim);
 if ~isempty(VarType{imax}.coord_x)  && ~isempty(VarType{imax}.coord_y)    %unstructured coordinates
@@ -2558 +2446 @@
 UvData.Field.YMin=YMin;
 
-%% create a default projection menuplane
-UvData.Object{1}.Style='plane';%main plotting plane
-UvData.Object{1}.ProjMode='projection';%main plotting plane
-UvData.Object{1}.DisplayHandle_uvmat=[]; %plane not visible in uvmat
-
 
 %% 3D case (menuvolume)
 if NbDim==3% && UvData.NewSeries
     test_set_object=1;
-    hset_object=findobj(allchild(0),'tag','set_object');
+    hset_object=findobj(allchild(0),'tag','set_object');% look for the set_object GUI
     ZBounds(1)=UvData.Field.ZMin; %minimum for the Z slider
     ZBounds(2)=UvData.Field.ZMax;%maximum for the Z slider
-    if ~isempty(hset_object)
-        hhset_object=guidata(hset_object);
-        ZBounds_old(1)=get(hhset_object.z_slider,'Min');
-        ZBounds_old(2)=get(hhset_object.z_slider,'Max');
-        if isequal(ZBounds_old,ZBounds)
-            test_set_object=0;
-        else
-            delete(test_set_object);
-        end
-    end
-    if test_set_object% reinitiate set_object
-        delete_object(1);
+    if ~isempty(hset_object) %if set_object is detected
+%         hhset_object=guidata(hset_object);
+% %         ZBounds_old(1)=get(hhset_object.z_slider,'Min');
+% %         ZBounds_old(2)=get(hhset_object.z_slider,'Max');
+% %         if isequal(ZBounds_old,ZBounds)
+%             test_set_object=0;% do not refresh the GUI set_object
+%         else
+            delete(hset_object);% delete the GUI set_object if it does not fit 
+%         end
+    end
+    if test_set_object% reinitiate the GUI set_object
+        delete_object(1);% delete the current projection object in the list UvData.Object, delete its graphic representations and update the list displayed in handles.list_object_1 and 2
         UvData.Object{1}.Style='plane';%main plotting plane
         UvData.Object{1}.ProjMode='projection';%main plotting plane
@@ -2597 +2480 @@
         set_object(UvData.Object{1},PlotHandles,ZBounds);
         set(handles.list_object_1,'Value',1);
+        set(handles.list_object_1,'String',{'1-PLANE'});
+        set(handles.edit_object,'Value',1)% put the plane in edit mode to enable the z cursor
+        edit_object_Callback([],[], handles)
     end
     %multilevel case (single menuplane in a 3D space)
@@ -2615 +2501 @@
         UvData.Object{1}.ZObject=UvData.ZIndex;
     end
+else
+    % create a default projection menuplane
+    UvData.Object{1}.Style='plane';%main plotting plane
+    UvData.Object{1}.ProjMode='projection';%main plotting plane
+    UvData.Object{1}.DisplayHandle_uvmat=[]; %plane not visible in uvmat
+    set(handles.list_object_1,'Value',1);
+    set(handles.list_object_1,'String',{'1-PLANE'});
 end
 if ~isfield(UvData.Object{1},'plotaxes')
-    UvData.Object{1}.plotaxes=handles.axes3;%default plotting axis
-    set(handles.list_object_1,'Value',1);
-    set(handles.list_object_1,'String',{'1-PLANE'});  
+    UvData.Object{1}.plotaxes=handles.axes3;%default plotting axis 
 end
 
@@ -2640 +2531 @@
 PlotParam{1}=read_plot_param(handles);%read plotting parameters on the uvmat interfac
 keeplim(1)=get(handles.FixedLimits,'Value');% test for fixed graph limits
-PosColorbar{1}=UvData.PosColorbar;%prescribe the colorbar position on the uvmat interface
+PosColorbar{1}=UvData.OpenParam.PosColorbar;%prescribe the colorbar position on the uvmat interface
 
 % second projection object (view_field display)
@@ -2729 +2620 @@
             hget_field=findobj(allchild(0),'name','get_field');
             if isempty(hget_field)
-                get_field([],ObjectData)% the projected field cannot be automatically plotted: use get_field to specify the variablesdelete(hget_field)
-            else
-                errormsg='The field defined by get_field cannot be plotted';
-                return
+                get_field(filename)% the projected field cannot be automatically plotted: use get_field to specify the variablesdelete(hget_field)
             end
+            errormsg='The field defined by get_field cannot be plotted';
+            return
         end
     end
@@ -2853 +2743 @@
         set(handles.Dt_txt,'String',['Dt=' num2str(1000*dt,3) '  m' UvData.TimeUnit] )
     end
-end
-
-%% update the input file name in the get_field GUI
-if isequal(FieldName,'get_field...')
-    set(hhget_field.inputfile,'String',filename)
-    Tabchar={''};%default
-    Tabcell=[];
-    if isfield(Field{1},'ListGlobalAttribute')&& ~isempty(Field{1}.ListGlobalAttribute)
-        for iline=1:length(Field{1}.ListGlobalAttribute)
-            Tabcell{iline,1}=Field{1}.ListGlobalAttribute{iline};
-            if isfield(Field{1}, Field{1}.ListGlobalAttribute{iline})
-                eval(['val=Field{1}.' Field{1}.ListGlobalAttribute{iline} ';'])
-                if ischar(val);
-                    Tabcell{iline,2}=val;
-                else
-                    Tabcell{iline,2}=num2str(val);
-                end
-            end
-        end
-        if ~isempty(Tabcell)
-            Tabchar=cell2tab(Tabcell,'=');
-            Tabchar=[{''};Tabchar];
-        end
-    end
-    set(hhget_field.attributes,'String',Tabchar);%update list of global attributes in get_field 
-end
-if isequal(FieldName_1,'get_field...')
-    set(hhget_field_1.inputfile,'String',filename_1)
-    Tabchar={''};%default
-    Tabcell=[];
-    if isfield(Field{2},'ListGlobalAttribute')&& ~isempty(Field{2}.ListGlobalAttribute)
-        for iline=1:length(Field{2}.ListGlobalAttribute)
-            Tabcell{iline,1}=Field{2}.ListGlobalAttribute{iline};
-            if isfield(Field{2}, Field{2}.ListGlobalAttribute{iline})
-                eval(['val=Field{2}.' Field{2}.ListGlobalAttribute{iline} ';'])
-                if ischar(val);
-                    Tabcell{iline,2}=val;
-                else
-                    Tabcell{iline,2}=num2str(val);
-                end
-            end
-        end
-        if ~isempty(Tabcell)
-            Tabchar=cell2tab(Tabcell,'=');
-            Tabchar=[{''};Tabchar];
-        end
-    end
-    set(hhget_field_1.attributes,'String',Tabchar);%update list of global attributes in get_field 
 end
 
@@ -3713 +3555 @@
     [hset_object,UvData.sethandles]=set_object(data,PlotHandles);% call the set_object interface with action on haxes,
                                                       % associate the set_object interface handle to the plotting axes
-    if isfield(UvData,'SetObjectOrigin')                                                
+    if isfield(UvData.OpenParam,'SetObjectOrigin')                                                
     pos_uvmat=get(huvmat,'Position');
-    pos_set_object(1:2)=UvData.SetObjectOrigin + pos_uvmat(1:2);
-    pos_set_object(3:4)=UvData.SetObjectSize .* pos_uvmat(3:4);  
+    pos_set_object(1:2)=UvData.OpenParam.SetObjectOrigin + pos_uvmat(1:2);
+    pos_set_object(3:4)=UvData.OpenParam.SetObjectSize .* pos_uvmat(3:4);  
     set(hset_object,'Position',pos_set_object)
     end
@@ -4365 +4207 @@
     UvData.MouseAction='none';
     set(handles.edit_object,'BackgroundColor',[0.7,0.7,0.7])   
+    hset_object=findobj(allchild(0),'tag','set_object');% look for the set_object GUI
+    if ~isempty(hset_object)
+        delete(hset_object)% delete the current GUI set_object
+    end
 end
 set(handles.uvmat,'UserData',UvData);
@@ -4410 +4256 @@
 function update_object(handles,IndexObj,option,ObjectName)
 %------------------------------------------------------------------------
-UvData=get(handles.uvmat,'UserData');%read UvData properties stored on the uvmat interface 
+UvData=get(handles.uvmat,'UserData');%read UvData properties stored on the uvmat interface
 if numel(UvData.Object)<IndexObj;
     return
@@ -4420 +4266 @@
 end
 ZBounds=0; % default
-    if isfield(UvData,'ZMin') && isfield(UvData,'ZMax')
-        ZBounds(1)=UvData.ZMin; %minimum for the Z slider
-        ZBounds(2)=UvData.ZMax;%maximum for the Z slider
-    end
+if isfield(UvData.Field,'ZMin') && isfield(UvData.Field,'ZMax')
+    ZBounds(1)=UvData.Field.ZMin; %minimum for the Z slider
+    ZBounds(2)=UvData.Field.ZMax;%maximum for the Z slider
+end
 hset_object=findobj(allchild(0),'tag','set_object');
 if ~isempty(hset_object)
     delete(hset_object)% delete existing version of set_object
+end
+edit_test=get(handles.edit_object,'Value');
+if edit_test
+    ObjectData.enable_plot=1;
+else
+        if isfield(ObjectData,'enable_plot')
+        ObjectData=rmfield(ObjectData,'enable_plot');
+    end
 end
 hset_object=set_object(ObjectData,[],ZBounds);% call the set_object interface,
@@ -4434 +4288 @@
 if option==1%length(UvData.Object)>= IndexObj && isfield(UvData.Object{IndexObj},'plotaxes')&& ishandle(UvData.Object{IndexObj}.plotaxes)
     PlotHandles=handles;
-%     UvData.Object{IndexObj}.plotaxes=handles.axes3;
 else
     hview_field=findobj(allchild(0),'tag','view_field');
@@ -4441 +4294 @@
     end
     PlotHandles=guidata(hview_field);
-%     hview_field=view_field(ProjData);
-%     hhviewfield=guidata(hviewfield);
-   % UvData.Object{IndexObj}.plotaxes=hhviewfield.axes3;
 end
 if ~isempty(ProjData)
-plot_field(ProjData,PlotHandles.axes3,PlotHandles);
+    plot_field(ProjData,PlotHandles.axes3,PlotHandles);
 end
 set(handles.uvmat,'UserData',UvData)
@@ -4457 +4307 @@
         end
     elseif isequal(get(hother(iobj),'Type'),'image')
-           Acolor=get(hother(iobj),'CData');
-           Acolor(:,:,1)=zeros(size(Acolor,1),size(Acolor,2));
-           set(hother(iobj),'CData',Acolor);
+        Acolor=get(hother(iobj),'CData');
+        Acolor(:,:,1)=zeros(size(Acolor,1),size(Acolor,2));
+        set(hother(iobj),'CData',Acolor);
     else
-         set(hother(iobj),'Color','b')
+        set(hother(iobj),'Color','b')
     end
     set(hother(iobj),'Selected','off')
@@ -4467 +4317 @@
 hother=findobj('Tag','DeformPoint');
 set(hother,'Color','b');
-set(hother,'Selected','off')    
+set(hother,'Selected','off')
 if isfield(ObjectData,'DisplayHandle_uvmat')
-%     for iview=1:length(ObjectData.HandlesDisplay)
-        if ishandle(ObjectData.DisplayHandle_uvmat)
-            uistack(ObjectData.DisplayHandle_uvmat,'top')
-            linetype=get(ObjectData.DisplayHandle_uvmat,'Type');
-            if isequal(linetype,'line')
-                set(ObjectData.DisplayHandle_uvmat,'Color','m'); %set the selected object to magenta color
-            elseif isequal(linetype,'rectangle')
-                set(ObjectData.DisplayHandle_uvmat,'EdgeColor','m'); %set the selected object to magenta color
-            elseif isequal(linetype,'patch')
-                set(ObjectData.DisplayHandle_uvmat,'FaceColor','m'); %set the selected object to magenta color
-            end 
-            SubObjectData=get(ObjectData.DisplayHandle_uvmat,'UserData');
-            if isfield(SubObjectData,'SubObject') & ishandle(SubObjectData.SubObject)
-                uistack(SubObjectData.SubObject,'top')
-                for iobj=1:length(SubObjectData.SubObject)
-                    hsub=SubObjectData.SubObject(iobj);
-                    if isequal(get(hsub,'Type'),'rectangle')
-                        set(hsub,'EdgeColor','m'); %set the selected object to magenta color
-                    elseif isequal(get(hsub,'Type'),'image')
-                       Acolor=get(hsub,'CData');
-                       Acolor(:,:,1)=Acolor(:,:,3);
-                       set(hsub,'CData',Acolor);
-                    else
-                        set(hsub,'Color','m')
-                    end
+    if ishandle(ObjectData.DisplayHandle_uvmat)
+        uistack(ObjectData.DisplayHandle_uvmat,'top')
+        linetype=get(ObjectData.DisplayHandle_uvmat,'Type');
+        if isequal(linetype,'line')
+            set(ObjectData.DisplayHandle_uvmat,'Color','m'); %set the selected object to magenta color
+        elseif isequal(linetype,'rectangle')
+            set(ObjectData.DisplayHandle_uvmat,'EdgeColor','m'); %set the selected object to magenta color
+        elseif isequal(linetype,'patch')
+            set(ObjectData.DisplayHandle_uvmat,'FaceColor','m'); %set the selected object to magenta color
+        end
+        SubObjectData=get(ObjectData.DisplayHandle_uvmat,'UserData');
+        if isfield(SubObjectData,'SubObject') & ishandle(SubObjectData.SubObject)
+            uistack(SubObjectData.SubObject,'top')
+            for iobj=1:length(SubObjectData.SubObject)
+                hsub=SubObjectData.SubObject(iobj);
+                if isequal(get(hsub,'Type'),'rectangle')
+                    set(hsub,'EdgeColor','m'); %set the selected object to magenta color
+                elseif isequal(get(hsub,'Type'),'image')
+                    Acolor=get(hsub,'CData');
+                    Acolor(:,:,1)=Acolor(:,:,3);
+                    set(hsub,'CData',Acolor);
+                else
+                    set(hsub,'Color','m')
                 end
             end
-            if isfield(SubObjectData,'DeformPoint') & ishandle(SubObjectData.DeformPoint)
-                set(SubObjectData.DeformPoint,'Color','m')
-            end
-        end
-%     end
+        end
+        if isfield(SubObjectData,'DeformPoint') & ishandle(SubObjectData.DeformPoint)
+            set(SubObjectData.DeformPoint,'Color','m')
+        end
+    end
+    %     end
 end
 pause(0.1)
@@ -4545 +4394 @@
 end
 
-
+%------------------------------------------------------------------------
 % --------------------------------------------------------------------
 function MenuExportMovie_Callback(hObject, eventdata, handles)
@@ -4631 +4480 @@
 msgbox_uvmat('CONFIRMATION',{['movie ' aviname ' created '];['with ' num2str(imax) ' frames']})
 
-
-% ------------------------------------------------------------------
+%------------------------------------------------------------------------
 function MenuCalib_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 
 UvData=get(handles.uvmat,'UserData');%read UvData properties stored on the uvmat interface 
@@ -4651 +4500 @@
 [FileName,RootPath,FileBase,FileIndices,FileExt,SubDir]=read_file_boxes(handles);
 set(handles.view_xml,'Backgroundcolor',[1 1 0])%indicate the reading of the current xml file by geometry_calib
-if isfield(UvData,'CalOrigin')
+if isfield(UvData.OpenParam,'CalOrigin')
     pos_uvmat=get(handles.uvmat,'Position');
-    pos_cal(1)=pos_uvmat(1)+UvData.CalOrigin(1)*pos_uvmat(3);
-    pos_cal(2)=pos_uvmat(2)+UvData.CalOrigin(2)*pos_uvmat(4);
-    pos_cal(3:4)=UvData.CalSize .* pos_uvmat(3:4);
+    pos_cal(1)=pos_uvmat(1)+UvData.OpenParam.CalOrigin(1)*pos_uvmat(3);
+    pos_cal(2)=pos_uvmat(2)+UvData.OpenParam.CalOrigin(2)*pos_uvmat(4);
+    pos_cal(3:4)=UvData.OpenParam.CalSize .* pos_uvmat(3:4);
 end
 geometry_calib(FileName,pos_cal);% call the geometry_calib interface    
@@ -4661 +4510 @@
 set(handles.MenuCalib,'checked','on')% indicate that MenuCalib is activated, test used by mouse action
 
-
-% ------------------------------------------------------------------
+%------------------------------------------------------------------------
 function MenuMask_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 UvData=get(handles.uvmat,'UserData');%read UvData properties stored on the uvmat interface 
 ListObj=UvData.Object;
@@ -4766 +4615 @@
 end
 
-%set(handles.uvmat,'UserData',UvData);
-
-% ------------------------------------------------------------------
+%------------------------------------------------------------------------
 %-- open the GUI set_grid.fig to create grid
 function MenuGrid_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 %UvData=get(handles.uvmat,'UserData');%read UvData properties stored on the uvmat interface 
 
@@ -4786 +4634 @@
 %set(handles.uvmat,'UserData',UvData);
 
-%----------------------------------------------------------------
+%------------------------------------------------------------------------
 % open the GUI 'series'
-%----------------------------------------------------------------
 function MenuSeries_Callback(hObject, eventdata, handles)
-%-------------------------------------------------------------------
+%------------------------------------------------------------------------
 series; %first display of the GUI to fill waiting time
 [param.FileName]=read_file_boxes(handles);
@@ -4835 +4682 @@
 series(param); %run the series interface
 
-% ------------------------------------------------------------------
+%------------------------------------------------------------------------
 % -- open the GUI civ.fig for civx (PIV)
-% ------------------------------------------------------------------
 function MenuPIV_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
  
 [FileName,RootPath,filebase,FileIndices,ext,SubDir]=read_file_boxes(handles);
@@ -4874 +4721 @@
 param.ImaExt=ext;
 civ(param);% interface de civ(not in the uvmat file)
-   
-% ------------------------------------------------------------------
+
+%------------------------------------------------------------------------
 function MenuTools_Callback(hObject, eventdata, handles)
-
-% ------------------------------------------------------------------
+%------------------------------------------------------------------------
+
+%------------------------------------------------------------------------
 function MenuEditObject_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 set(handles.edit_object,'Value',1)
 edit_Callback(hObject, eventdata, handles)
 
-%--------------------------------------------------------------------------
+%------------------------------------------------------------------------
 function enable_transform(handles,state)
+%------------------------------------------------------------------------
 set(handles.transform_fct,'Visible',state)
 set(handles.TRANSFORM_txt,'Visible',state)    
@@ -4897 +4747 @@
 function MenuEditVectors_Callback(hObject, eventdata, handles)
 %------------------------------------------------------------------------
-% set(handles.MenuEditVectors,'checked','on')
 set(handles.edit_vect,'Visible','on')
 set(handles.edit_vect,'Value',1)
 edit_vect_Callback(hObject, eventdata, handles)
 
-% --------------------------------------------------------------------
+% -----------------------------------------------------------------------
 function Menupoints_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 data.Style='points';
 data.ProjMode='projection';%default
 create_object(data,handles)
 
-% --------------------------------------------------------------------
+% -----------------------------------------------------------------------
 function Menuline_Callback(hObject, eventdata, handles)
-% set(handles.create,'Visible','on')
-% set(handles.create,'Value',1)
-% LINE_Callback(hObject,eventdata,handles)
+%------------------------------------------------------------------------
 data.Style='line';
 data.ProjMode='projection';%default
@@ -4954 +4802 @@
 create_object(data,handles)
 
-% ------------------------------------------------------------------
+%------------------------------------------------------------------------
 function Menuplane_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 data.Style='plane';
 data.ProjMode='projection';%default
@@ -4961 +4810 @@
 create_object(data,handles)
 
-% ------------------------------------------------------------------
+%------------------------------------------------------------------------
 function Menuvolume_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 data.Style='volume';
 data.ProjMode='interp';%default
@@ -5109 +4959 @@
 end
 
-
+%------------------------------------------------------------------------
 % --- Executes on button press in delete_object.
 function delete_object_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 IndexObj=get(handles.list_object_2,'Value');
-    if IndexObj>1 
-        delete_object(IndexObj)
-    end
+if IndexObj>1 
+    delete_object(IndexObj)
+end
     
 

trunk/src/view_field.m

@@ -74 +74 @@
     hhuvmat=guidata(huvmat);
     set(hhuvmat.list_object_2,'Visible','on')
+    set(hhuvmat.view_field_title,'Visible','on')
     % handles_mouse.create=hhuvmat.create;
     handles_mouse.edit=hhuvmat.edit_object;