source: trunk/src/get_field.m @ 1086

Last change on this file since 1086 was 1084, checked in by sommeria, 4 years ago

various updates

File size: 47.1 KB
Line 
1%'get_field': display variables and attributes from a Netcdf file, and OK selected fields
2%------------------------------------------------------------------------
3% GetFieldData=get_field(FileName,ParamIn)
4% associated with the GUI get_field.fig
5%
6% OUTPUT:
7% GetFieldData: structure containing the information on the selected
8%      fields, obtained by applying the fct red_GUI to the GUI get_field
9%   .FieldOption='vectors': variables are used for vector plot
10%                  'scalar': variables are used for scalar plot,
11%                  '1Dplot': variables are used for usual x-y plot,
12%                  'civdata...': go back to automatic reading of civ data
13%   .PanelVectors: sub-structure variables used as vector components
14%   .PanelScalar:
15% INPUT:
16% FileName: name (including path) of the netcdf file to open
17% ParmIn: structure containing parameters for preselecting menus:
18%   .Title: set the title of the GUI get_field
19%   .SwitchVarIndexTime='file index','variable' or 'matrix index': select the default option for 'time'
20%   .TimeAttrName: preselect the name of a global attribute for time
21%   .SeriesInput=1 if get_field is called by the GUI series,=0 otherwise (plot options provided in the latter case)
22%   .Coord_x,.Coord_y,.Coord_z, names of the variables used as the three coordinates
23%   .scalar : set the default choise of the scale variable
24%   .vector_x, .vector_y : set the default choise for the variables used for the x and y vector components
25
26%=======================================================================
27% Copyright 2008-2020, LEGI UMR 5519 / CNRS UGA G-INP, Grenoble, France
28%   http://www.legi.grenoble-inp.fr
29%   Joel.Sommeria - Joel.Sommeria (A) legi.cnrs.fr
30%
31%     This file is part of the toolbox UVMAT.
32%
33%     UVMAT is free software; you can redistribute it and/or modify
34%     it under the terms of the GNU General Public License as published
35%     by the Free Software Foundation; either version 2 of the license,
36%     or (at your option) any later version.
37%
38%     UVMAT is distributed in the hope that it will be useful,
39%     but WITHOUT ANY WARRANTY; without even the implied warranty of
40%     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
41%     GNU General Public License (see LICENSE.txt) for more details.
42%=======================================================================
43
44function varargout = get_field(varargin)
45
46% Last Modified by GUIDE v2.5 18-Feb-2015 23:42:12
47
48% Begin initialization code - DO NOT EDIT
49gui_Singleton = 1;
50gui_State = struct('gui_Name',       mfilename, ...
51                   'gui_Singleton',  gui_Singleton, ...
52                   'gui_OpeningFcn', @get_field_OpeningFcn, ...
53                   'gui_OutputFcn',  @get_field_OutputFcn, ...
54                   'gui_LayoutFcn',  [] , ...
55                   'gui_Callback',   []);
56if nargin && ischar(varargin{1})&& ~isempty(regexp(varargin{1},'_Callback','once'))
57    gui_State.gui_Callback = str2func(varargin{1});
58end
59
60if nargout
61    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
62else
63    gui_mainfcn(gui_State, varargin{:});
64end
65% End initialization code - DO NOT EDIT
66     
67%------------------------------------------------------------------------
68% --- Executes just before get_field is made visible.
69%------------------------------------------------------------------------
70function get_field_OpeningFcn(hObject, eventdata, handles,filename,ParamIn)
71
72%% GUI settings
73handles.output = 'Cancel';
74guidata(hObject, handles);
75set(hObject,'WindowButtonDownFcn',{'mouse_down'}) % allows mouse action with right button (zoom for uicontrol display)
76set(hObject,'CloseRequestFcn',{@closefcn,handles})
77
78%% enter input data
79if ischar(filename) % input file name
80    set(handles.inputfile,'String',filename)% fill the input file name
81    [Field,tild,tild,errormsg]=nc2struct(filename,[]);% reads the  field structure, without the variables
82else
83    msgbox_uvmat('ERROR','get_field requires a file name as input')% display error message for input file reading
84    return
85end
86if ~isempty(errormsg)
87    msgbox_uvmat('ERROR',['get_field/nc2struct/' errormsg])% display error message for input file reading
88    return
89end
90if ~isfield(Field,'ListVarName')
91    msgbox_uvmat('ERROR',['no variable found in ' filename])% display error message for input file reading
92    return
93end
94if ~exist('ParamIn','var')
95    ParamIn.Coord_z='';
96end
97
98%% look at singletons and variables with a single dimension
99Field.Display=Field;
100Field.Check0D=zeros(size(Field.ListVarName));% =1 for arrays with a single value
101NbVar=numel(Field.VarDimName);%nbre of variables in the input data
102for ilist=1:NbVar
103    if ischar(Field.VarDimName{ilist})
104        Field.VarDimName{ilist}={Field.VarDimName{ilist}}; %transform string into cell
105    end
106    NbDim=numel(Field.VarDimName{ilist});
107    check_singleton=false(1,NbDim);%  check singleton, false by default
108    for idim=1:NbDim
109        dim_index=strcmp(Field.VarDimName{ilist}{idim},Field.ListDimName);%index in the list of dimensions
110        check_singleton(idim)=isequal(Field.DimValue(dim_index),1);%check_singleton=1 for singleton
111    end
112    Field.Check0D(ilist)=(isequal(check_singleton,ones(1,NbDim)))||(~isequal(Field.VarType(ilist),4)&&~isequal(Field.VarType(ilist),5)&&~isequal(Field.VarType(ilist),6));% =1 if the variable reduces to a single value
113    if ~Field.Check0D(ilist)
114    Field.Display.VarDimName{ilist}=Field.VarDimName{ilist}(~check_singleton);% eliminate singletons in the list of variable dimensions
115    end
116end
117if ~isfield(Field,'VarAttribute')
118    Field.VarAttribute={};
119end
120if numel(Field.VarAttribute)<NbVar% complement VarAttribute by blanjs if neded
121    Field.VarAttribute(numel(Field.VarAttribute)+1:NbVar)=cell(1,NbVar-numel(Field.VarAttribute));
122end
123% Field.Display = list of variables and corresponding properties obtained after removal of variables with a single value and singleton dimensions
124Field.Display.ListVarName=Field.ListVarName(~Field.Check0D); %list of variables available for plots, after eliminating variables with a single value
125Field.Display.VarAttribute=Field.VarAttribute(~Field.Check0D);
126Field.Display.VarDimName=Field.Display.VarDimName(~Field.Check0D);
127Field.Display.ListDimName=Field.ListDimName(Field.DimValue~=1);% list of non singleton dimension names
128Field.Display.DimValue=Field.DimValue(Field.DimValue~=1);% corresponding list of non singleton dimension values
129
130
131%% analyse the input field cells
132[CellInfo,NbDim,errormsg]=find_field_cells(Field.Display);
133if ~isempty(errormsg)
134    msgbox_uvmat('ERROR',['get_field / Field_input / find_field_cells: ' errormsg])
135    return
136end
137if isempty(CellInfo)
138    [Field.MaxDim,imax]=max(cellfun(@numel,Field.Display.VarDimName));% maximum number of dimensions for the input fields
139    ListDim=Field.Display.VarDimName{imax};
140    check_cellinfo=false;
141else
142    [Field.MaxDim,imax]=max(NbDim);% maximum number of dimensions for the input fields identified by attributes
143    check_cellinfo=true;
144end
145
146%% set time mode
147ListSwitchVarIndexTime={'file index'};% default setting: the time is the file index
148% look at global attributes with numerical values
149check_numvalue=false(1,numel(Field.ListGlobalAttribute));
150for ilist=1:numel(Field.ListGlobalAttribute)
151    Value=Field.(Field.ListGlobalAttribute{ilist});
152    check_numvalue(ilist)=isnumeric(Value);
153end
154Field.Display.ListGlobalAttribute=Field.ListGlobalAttribute(check_numvalue);% select the attributes with float numerical value
155if ~isempty(Field.Display.ListGlobalAttribute)
156    ListSwitchVarIndexTime=[ListSwitchVarIndexTime; {'attribute'}];% the time can be chosen as a global attribute
157end
158
159Check_index=0;
160if Field.MaxDim>=2
161    ListSwitchVarIndexTime=[ListSwitchVarIndexTime;{'variable'};{'matrix index'}];% the time can be chosen as a dim index
162else
163    for ilist=1:numel(Field.Display.VarDimName)
164        NbComponent=numel(Field.Display.VarDimName{ilist});
165        if NbComponent>=2% multicomponent matrices without coordinate variables (thus not considered in the fct find_field_cell)
166            ListSwitchVarIndexTime=[ListSwitchVarIndexTime;{'matrix index'}];% the time can be chosen as a dim index
167            Check_index=1;
168            break
169        end
170    end
171end
172
173%% select the Time attribute from input
174if Field.MaxDim >2
175    variable_index=find(strcmp('variable',ListSwitchVarIndexTime),1);
176    set(handles.SwitchVarIndexTime,'Value',variable_index);
177else
178    if isfield(ParamIn,'TimeAttrName')&& ~isempty(ParamIn.TimeAttrName)
179        time_index=find(strcmp(ParamIn.TimeAttrName,Field.Display.ListGlobalAttribute),1);
180    else
181        time_index=find(strcmp('Time',Field.Display.ListGlobalAttribute));% look for global attribute containing name 'Time'
182    end
183    if isempty(time_index)
184        set(handles.SwitchVarIndexTime,'Value',1);
185    else
186        set(handles.SwitchVarIndexTime,'Value',2);
187        set(handles.TimeName,'UserData',time_index)
188    end
189end
190set(handles.SwitchVarIndexTime,'String',ListSwitchVarIndexTime)
191set(handles.SwitchVarIndexTime,'UserData',ListSwitchVarIndexTime); % keep string in memory for check3D
192set(handles.get_field,'UserData',Field);% record the finput field structure
193SwitchVarIndexTime_Callback([], [], handles)
194
195%% set vector menu (priority) if detected or scalar menu for space dim >=2, or usual (x,y) plot for 1D fields
196set(handles.vector_x,'String',Field.Display.ListVarName)% fill the menu of x vector components
197set(handles.vector_y,'String',Field.Display.ListVarName)% fill the menu of y vector components
198set(handles.vector_z,'String',[{''} Field.Display.ListVarName])% fill the menu of y vector components
199set(handles.vec_color,'String',[{''} Field.Display.ListVarName])% fill the menu of y vector components
200set(handles.scalar,'Value',1)% fill the menu of y vector components
201set(handles.scalar,'String',Field.Display.ListVarName)% fill the menu for scalar
202%set(handles.ordinate,'Value',1)% fill the menu of y vector components
203%set(handles.ordinate,'String',Field.Display.ListVarName)% fill the menu of y coordinate for 1D plots
204checkseries=0;
205if isfield(ParamIn,'SeriesInput') && ParamIn.SeriesInput% case of call by series
206    set(handles.FieldOption,'value',1)
207    if isfield(Field,'Conventions')&& strcmp(Field.Conventions,'uvmat/civdata')
208    set(handles.FieldOption,'String',{'scalar';'vectors';'civdata...'})
209    else
210       set(handles.FieldOption,'String',{'scalar';'vectors'})
211    end
212    checkseries=1;
213    set(handles.scalar,'Max',2)
214elseif isfield(Field,'Conventions')&& strcmp(Field.Conventions,'uvmat/civdata')
215    set(handles.FieldOption,'String',{'1D plot';'scalar';'vectors';'civdata...'})% provides the possibility to come back to civdata
216    set(handles.scalar,'Max',1)
217else
218    set(handles.FieldOption,'String',{'1D plot';'scalar';'vectors'})
219    set(handles.scalar,'Max',1)
220end
221
222%% set default field options
223checknbdim=cellfun('size',Field.Display.VarDimName,2);
224% if max(checknbdim)<=1
225%     Field.MaxDim=1;% only 1D fields, considered as a time series by default
226% end
227if Field.MaxDim>=2 && ~checkseries% case of 2D (or 3D) fields
228    check_vec_input=0;
229    % case of vector initially selected from uvmat input
230    if isfield(ParamIn,'vector_x')&& isfield(ParamIn,'vector_y')
231        ichoice_x=find(strcmp(ParamIn.vector_x,Field.Display.ListVarName),1);
232        ichoice_y=find(strcmp(ParamIn.vector_y,Field.Display.ListVarName),1);
233        if ~isempty(ichoice_x)&&~isempty(ichoice_y)
234            set(handles.vector_x,'UserData',ichoice_x)
235            set(handles.vector_y,'UserData',ichoice_y)
236            check_vec_input=1;
237        end
238    end
239    % otherwise select vectors marked as attributes in the input field
240    if check_cellinfo && ~check_vec_input && isfield(CellInfo{imax},'VarIndex_vector_x') &&  isfield(CellInfo{imax},'VarIndex_vector_y')
241        set(handles.vector_x,'UserData',CellInfo{imax}.VarIndex_vector_x(1))
242        set(handles.vector_y,'UserData',CellInfo{imax}.VarIndex_vector_y(1))
243        check_vec_input=1;
244    end
245    if check_vec_input
246        set(handles.FieldOption,'Value',3)% set vector selection option
247    else     
248        set(handles.FieldOption,'Value',2)% set scalar selection option
249    end
250else % case of 1D fields
251    set(handles.FieldOption,'Value',1)
252end
253
254%% fill the general list of dimensions, variables, attributes
255if isfield(Field,'ListDimName')&&~isempty(Field.ListDimName)
256    Tabcell(:,1)=Field.ListDimName;
257    for iline=1:length(Field.ListDimName)
258        Tabcell{iline,2}=num2str(Field.DimValue(iline));
259    end
260    Tabchar=cell2tab(Tabcell,' = ');
261    set(handles.dimensions,'String',Tabchar)
262end
263
264%% fill menus for coordinates and time
265FieldOption_Callback(handles.variables,[], handles)% list the global attributes
266
267%% put the GUI on the lower right of the sceen
268set(hObject,'Unit','pixels')
269%pos_view_field=get(hObject,'Position');
270set(0,'Unit','pixels')
271ScreenSize=get(0,'ScreenSize');
272pos_view_field(3:4)=[955 648];
273pos_view_field(1)=ScreenSize(1)+ScreenSize(3)-pos_view_field(3);
274pos_view_field(2)=ScreenSize(2);
275set(hObject,'Position',pos_view_field)
276set(handles.get_field,'WindowStyle','modal')% Make the GUI modal
277if isfield(ParamIn,'Title')
278    set(hObject,'Name',ParamIn.Title)
279end
280
281%% set z coordinate menu if relevant
282if Field.MaxDim>=3 && prod(Field.DimValue)<10^8 && ~(isfield(ParamIn,'Coord_z') && isempty(ParamIn.Coord_z)) % 3D field (with memory content smaller than 400 Mo)
283    set(handles.Check3D,'Value',1)
284else
285    set(handles.Check3D,'Value',0)
286end
287Check3D_Callback(hObject, eventdata, handles)
288set(handles.variables,'Value',1)
289set(handles.variables,'String',[{'*'} Field.ListVarName])
290variables_Callback(handles.variables,[], handles)% list the global attributes
291drawnow
292uiwait(handles.get_field);
293
294% -----------------------------------------------------------------------
295% --- Activated by selection in the list of variables
296% ----------------------------------------------------------------------
297function variables_Callback(hObject, VarName, handles)
298
299Tabchar={''};%default
300Tabcell=[];
301hselect_field=get(handles.variables,'parent');
302Field=get(handles.get_field,'UserData');
303index=get(handles.variables,'Value');%index in the list 'variables'
304
305%% list global TimeAttribute names and values if index=1 (blank TimeVariable display) is selected
306if isequal(index,1)
307    set(handles.attributes_txt,'String','global attributes')
308    if isfield(Field,'ListGlobalAttribute') && ~isempty(Field.ListGlobalAttribute)
309        for iline=1:length(Field.ListGlobalAttribute)
310            Tabcell{iline,1}=Field.ListGlobalAttribute{iline};
311            if isfield(Field, Field.ListGlobalAttribute{iline})
312                val=Field.(Field.ListGlobalAttribute{iline});
313                if ischar(val);% attribute value is char string
314                    Tabcell{iline,2}=val;
315                elseif size(val,1)==1 %attribute value is a number or matlab vector
316                    Tabcell{iline,2}=num2str(val);
317                end
318            end
319        end
320        Tabchar=cell2tab(Tabcell,'=');
321    end
322    %% list Attribute names and values associated to the Variable # index-1
323else
324    list_var=get(handles.variables,'String');
325    if index>numel(list_var)
326        return
327    end
328    VarName=list_var{index};
329    set(handles.attributes_txt,'String', ['attributes of ' VarName])
330    if isfield(Field,'VarAttribute')&& length(Field.VarAttribute)>=index-1
331        VarAttr=Field.VarAttribute{index-1};
332        if isstruct(VarAttr)
333            attr_list=fieldnames(VarAttr);
334            for iline=1:length(attr_list)
335                Tabcell{iline,1}=attr_list{iline};
336                val=VarAttr.(attr_list{iline}) ;
337                if ischar(val);
338                    Tabcell{iline,2}=val;
339                else
340                    Tabcell{iline,2}=num2str(val);
341                end
342            end
343        end
344    end
345end
346if ~isempty(Tabcell)
347    Tabchar=cell2tab(Tabcell,'=');
348end
349set(handles.attributes,'Value',1);% select the first item
350set(handles.attributes,'String',Tabchar);
351
352%% update dimensions;
353if isfield(Field,'ListDimName')
354    Tabdim={};%default
355    if isequal(index,1)%list all dimensions if '*' is selected as the variable
356        dim_indices=1:length(Field.ListDimName);
357        set(handles.dimensions_txt,'String', 'dimensions')
358    else   % a specific variable has been selected
359        DimCell=Field.VarDimName{index-1};
360        if ischar(DimCell)
361            DimCell={DimCell};% transform into a cell for a single dimension defined by a char string
362        end
363        dim_indices=[];
364        for idim=1:length(DimCell)
365            dim_index=strcmp(DimCell{idim},Field.ListDimName);%vector with size of Field.ListDimName, =0
366            dim_index=find(dim_index,1);
367            dim_indices=[dim_indices dim_index];
368        end
369        set(handles.dimensions_txt,'String', ['dimensions of ' VarName])
370    end
371    for iline=1:length(dim_indices)
372        Tabdim{iline,1}=Field.ListDimName{dim_indices(iline)};
373        Tabdim{iline,2}=num2str(Field.DimValue(dim_indices(iline)));
374    end
375    Tabchar=cell2tab(Tabdim,' = ');
376    Tabchar=[{''} ;Tabchar];
377    set(handles.dimensions,'Value',1)
378    set(handles.dimensions,'String',Tabchar)
379end
380
381%% propose a plot by default if variables_Callback has not been already called by FieldOption_Callback (VarName is not a char string)
382if ~ischar(VarName) && ~isequal(index,1)
383    if numel(DimCell)==1
384        set(handles.FieldOption,'Value',1)%propose 1D plot
385    else
386        set(handles.FieldOption,'Value',2)%propose scalar plot
387    end
388    if numel(DimCell)<=2
389        set(handles.Check3D,'Value',0)
390    else
391        set(handles.Check3D,'Value',1)
392    end
393    FieldOption_Callback(hObject, VarName, handles)
394end
395
396%------------------------------------------------------------------------
397% --- Executes on selection change in FieldOption.
398%------------------------------------------------------------------------
399function FieldOption_Callback(hObject, VarName, handles)
400
401Field=get(handles.get_field,'UserData');
402FieldList=get(handles.FieldOption,'String');
403FieldOption=FieldList{get(handles.FieldOption,'Value')};
404switch FieldOption
405    case '1D plot'
406        set(handles.Coordinates,'Visible','on')
407        %set(handles.PanelOrdinate,'Visible','on')
408        %pos=get(handles.PanelOrdinate,'Position');
409%         pos(1)=2;
410%         pos_coord=get(handles.Coordinates,'Position');
411%         pos(2)=pos_coord(2)-pos(4)-2;
412        %set(handles.PanelOrdinate,'Position',pos)
413        set(handles.PanelScalar,'Visible','off')
414        set(handles.PanelVectors,'Visible','off')
415        set(handles.Coord_y,'Visible','on')
416        set(handles.Coord_y,'Max',2)%allow multiple selection
417        set(handles.Y_title,'Visible','on')
418        set(handles.Coord_z,'Visible','off')
419        set(handles.Z_title,'Visible','off')
420        set(handles.Coord_x,'String',Field.Display.ListVarName')
421        Coord_x_Callback(hObject, VarName, handles)
422        %set(handles.Coord_y,'String',Field.Display.ListVarName')
423        %Coord_x_Callback(hObject, VarName, handles)       
424    case {'scalar'}
425        set(handles.Coordinates,'Visible','on')
426        %set(handles.PanelOrdinate,'Visible','off')
427        set(handles.PanelScalar,'Visible','on')
428        set(handles.PanelVectors,'Visible','off')
429        pos=get(handles.PanelScalar,'Position');
430        pos(1)=2;
431        pos_coord=get(handles.Coordinates,'Position');
432        pos(2)=pos_coord(2)-pos(4)-2;
433        set(handles.PanelScalar,'Position',pos)
434        set(handles.Coord_y,'Visible','on')
435        set(handles.Y_title,'Visible','on')     
436        if ~ischar(VarName)     
437            %default scalar selection
438            test_coord=zeros(size(Field.Display.VarDimName)); %=1 when variable #ilist is eligible as structured coordiante
439            for ilist=1:numel(Field.Display.VarDimName)
440                if isfield(Field.Display,'VarAttribute') && numel(Field.Display.VarAttribute)>=ilist && isfield(Field.Display.VarAttribute{ilist},'Role')
441                    Role=Field.Display.VarAttribute{ilist}.Role;
442                    if strcmp(Role,'coord_x')||strcmp(Role,'coord_y')
443                        test_coord(ilist)=1;
444                    end
445                end
446                dimnames=Field.Display.VarDimName{ilist}; %list of dimensions for variable #ilist
447                if numel(dimnames)==1 && strcmp(dimnames{1},Field.Display.ListVarName{ilist})%dimension variable
448                    test_coord(ilist)=1;
449                end
450            end
451            scalar_index=find(~test_coord,1);%get the first variable not a coordinate
452            if isempty(scalar_index)
453                set(handles.scalar,'Value',1)
454            else
455                set(handles.scalar,'Value',scalar_index)
456            end
457        end
458        scalar_Callback(hObject,VarName, handles)       
459    case 'vectors'
460        set(handles.PanelVectors,'Visible','on')
461        set(handles.Coordinates,'Visible','on')
462        %set(handles.PanelOrdinate,'Visible','off')
463        set(handles.PanelScalar,'Visible','off')
464        pos=get(handles.PanelVectors,'Position');
465        pos(1)=2;
466        pos_coord=get(handles.Coordinates,'Position');
467        pos(2)=pos_coord(2)-pos(4)-2;
468        set(handles.PanelVectors,'Position',pos)
469        set(handles.Coord_y,'Visible','on')
470        set(handles.Y_title,'Visible','on')
471        %default vector selection
472        vector_x_value=get(handles.vector_x,'UserData');
473        vector_y_value=get(handles.vector_y,'UserData');
474        if ~isempty(vector_x_value)&&~isempty(vector_y_value)
475            set(handles.vector_x,'Value',vector_x_value)
476            set(handles.vector_y,'Value',vector_y_value)
477        else
478            test_coord=zeros(size(Field.Display.VarDimName)); %=1 when variable #ilist is eligible as structured coordinate
479            for ilist=1:numel(Field.Display.VarDimName)
480                if isfield(Field.Display,'VarAttribute') && numel(Field.Display.VarAttribute)>=ilist && isfield(Field.Display.VarAttribute{ilist},'Role')
481                    Role=Field.Display.VarAttribute{ilist}.Role;
482                    if strcmp(Role,'coord_x')||strcmp(Role,'coord_y')
483                        test_coord(ilist)=1;
484                    end
485                end
486                dimnames=Field.Display.VarDimName{ilist}; %list of dimensions for variable #ilist
487                if numel(dimnames)==1 && strcmp(dimnames{1},Field.Display.ListVarName{ilist})%dimension variable
488                    test_coord(ilist)=1;
489                end
490            end
491            vector_index=find(~test_coord,2);%get the two first variables not a coordinate
492            if isempty(vector_index)
493                set(handles.vector_x,'Value',1)
494                set(handles.vector_y,'Value',2)
495            else
496                set(handles.vector_x,'Value',vector_index(1))
497                set(handles.vector_y,'Value',vector_index(2))
498            end
499        end
500        vector_Callback(handles)     
501    case 'civdata...'
502        %set(handles.PanelOrdinate,'Visible','off')
503        set(handles.PanelScalar,'Visible','off')
504        set(handles.PanelVectors,'Visible','off')
505        set(handles.Coordinates,'Visible','off')
506end
507
508
509function set_coord_y_options(handles,VarName)
510%------------------------------------------------------------------------
511Field=get(handles.get_field,'UserData');
512VarIndex=find(strcmp(VarName,Field.Display.ListVarName),1);
513DimCell=Field.Display.VarDimName{VarIndex};
514% y_index=get(handles.Coord_y,'Value');
515% y_menu=get(handles.Coord_y,'String');
516% if isempty(y_menu)
517%     return
518% else
519% YName=y_menu{y_index};
520% end
521
522%% set list of possible coordinates
523% test_component=zeros(size(Field.Display.VarDimName));%=1 when variable #ilist is eligible as unstructured coordinate
524test_coord=zeros(size(Field.Display.VarDimName)); %=1 when variable #ilist is eligible as structured coordiante
525% ListCoord={''};
526% dim_var=Field.Display.VarDimName{y_index};%list of dimensions of the selected variable
527
528for ilist=1:numel(Field.Display.VarDimName)
529    dimnames=Field.Display.VarDimName{ilist}; %list of dimensions for variable #ilist
530    if isequal(dimnames,DimCell)||isequal(dimnames(1:end-1),DimCell)||isequal(dimnames(2:end),DimCell)
531        test_coord(ilist)=1;
532    end
533end
534ListCoord=Field.Display.ListVarName(find(test_coord));
535set(handles.Coord_y,'String',ListCoord)
536val_y=1;
537if strcmp(VarName,ListCoord{1})&& numel(ListCoord)>=2
538    val_y=2;
539end
540set(handles.Coord_y,'Value',val_y)
541
542%% set default coord selection
543% if numel(find(test_coord))>3
544%      SwitchVarIndexTime=get(handles.SwitchVarIndexTime,'String');
545%     if numel(SwitchVarIndexTime)<3
546%         SwitchVarIndexTime=[SwitchVarIndexTime;'matrix_index'];
547%         set(handles.SwitchVarIndexTime,'String',SwitchVarIndexTime)
548%     end
549%     set(handles.SwitchVarIndexTime,'Value',3)% the last dim must be considered as time
550%     SwitchVarIndexTime_Callback([], [], handles)
551% end
552% if numel(var_component)<2
553%     if numel(test_coord)<2
554%         ListCoord={''};
555%     else
556%         set(handles.Coord_x,'Value',2)
557%         set(handles.Coord_y,'Value',1)
558%     end
559% else
560%     coord_val=1;
561%     for ilist=1:numel(var_component)
562%         ivar=var_component(ilist);
563%         if isfield(Field.Display,'VarAttribute') && numel(Field.Display.VarAttribute)>=ivar && isfield(Field.Display.VarAttribute{ivar},'Role')
564%             Role=Field.Display.VarAttribute{ivar}.Role;
565%             if strcmp(Role,'coord_x')
566%                 coord_val=ilist;
567%             end
568%         end
569%     end
570%     set(handles.Coord_x,'Value',coord_val+1)
571% end
572% set(handles.Coord_x,'String',[{''}; ListCoord])
573
574
575% %% set list of time coordinates
576% menu=get(handles.SwitchVarIndexTime,'String');
577% TimeOption=menu{get(handles.SwitchVarIndexTime,'Value')};
578% switch TimeOption
579%     case 'variable'
580%         if numel(find(test_coord))<3
581%             ListTime={''};
582%         else
583%             ListTime=Field.Display.ListVarName(find(test_coord,end));
584%         end
585%         set(handles.TimeName,'Value',1)
586%         set(handles.TimeName,'String',ListTime)
587%     case 'matrix index'
588%         if numel(find(test_coord))<3
589%             ListTime={''};
590%         else
591%             ListTime=Field.Display.VarDimName{find(test_coord,end)};
592%         end
593%         set(handles.TimeName,'Value',1)
594%         set(handles.TimeName,'String',ListTime)
595% end 
596% if ~ischar(DimCell)
597% update_field(handles,YName)
598% end
599         
600%------------------------------------------------------------------------
601% --- Executes on selection change in scalar menu.
602%------------------------------------------------------------------------
603function scalar_Callback(hObject, VarName, handles)
604
605Field=get(handles.get_field,'UserData');% get the input field info stored in UserData of the GUI
606scalar_menu=get(handles.scalar,'String');% read the menu for scalar selection
607if ischar(VarName)% case of a call with input variable
608    ScalarName=VarName;
609    scalar_index=find(strcmp(VarName,scalar_menu));
610    set(handles.scalar,'Value',scalar_index)% select the input variable field in the menu
611else % no input variable, the variable ScalarName is selected from the menu
612    scalar_index=get(handles.scalar,'Value');
613    ScalarName=scalar_menu{scalar_index};
614end
615
616%% set list of possible coordinates
617test_component=zeros(size(Field.Display.VarDimName));%=1 when variable #ilist is eligible as unstructured coordinate
618test_coord=zeros(size(Field.Display.VarDimName)); %=1 when variable #ilist is eligible as structured coordiante
619dim_var=Field.Display.VarDimName{scalar_index};%list of dimensions of the selected variable
620%if ~get(handles.CheckDimensionX,'Value')
621%look for coordinate variables among the other variables
622for ilist=1:numel(Field.Display.VarDimName)
623    dimnames=Field.Display.VarDimName{ilist}; %list of dimensions for variable #ilist
624    if isequal(dimnames,dim_var)
625        test_component(ilist)=1;% the listed variable has the same dimension as the selected scalar-> possibly chosen as unstructured coordinate
626    elseif numel(dimnames)==1 && ~isempty(find(strcmp(dimnames{1},dim_var), 1))%variable ilist is a 1D array which can be coordinate variable
627        test_coord(ilist)=1;
628    end
629end
630%end
631var_component=find(test_component);% list of variable indices elligible as unstructured coordinates
632var_coord=find(test_coord);% % list of variable indices elligible as gridded coordinates
633var_coord(var_coord==scalar_index)=[];
634var_component(var_component==scalar_index)=[];
635ListCoord=Field.Display.ListVarName([var_coord var_component]);
636coord_val=zeros(size(ListCoord));
637
638%% set default selection for grid coordinates
639if numel(var_coord)>=2
640    coord_val(1)=var_coord(end);
641    coord_val(2)=var_coord(end-1);
642    if numel(var_coord)>=3
643        coord_val(3)=var_coord(end-2);
644    end
645end
646% if numel(find(test_coord))>3
647%     SwitchVarIndexTime=get(handles.SwitchVarIndexTime,'String');
648%     if numel(SwitchVarIndexTime)<3
649%         SwitchVarIndexTime=[SwitchVarIndexTime;'matrix_index'];
650%         set(handles.SwitchVarIndexTime,'String',SwitchVarIndexTime)
651%     end
652%     set(handles.SwitchVarIndexTime,'Value',3)% the last dim must be considered as time
653%     SwitchVarIndexTime_Callback([], [], handles)
654% end
655
656%% default selection for labelled unstructured coordinates
657for ilist=1:numel(var_component)
658    ivar=var_component(ilist);
659    if isfield(Field.Display,'VarAttribute') && numel(Field.Display.VarAttribute)>=ivar && isfield(Field.Display.VarAttribute{ivar},'Role')
660        Role=Field.Display.VarAttribute{ivar}.Role;
661        if strcmp(Role,'coord_x')
662            coord_val(1)=ilist;
663        elseif strcmp(Role,'coord_y')
664            coord_val(2)=ilist;
665        elseif strcmp(Role,'coord_z')
666            coord_val(3)=ilist;
667        end
668    end
669end
670if numel(find(coord_val))<2 % no predefiend components
671    if numel(var_coord)>=3
672        coord_val(3)=3;
673    end
674    coord_val([1 2])=[1 2];
675end
676
677%% set menu and default selection for coordinates
678set(handles.Coord_x,'Value',coord_val(1))
679set(handles.Coord_x,'String',ListCoord)
680set(handles.Coord_y,'Value',coord_val(2))
681set(handles.Coord_y,'String',ListCoord)
682if numel(find(coord_val))>=3
683    set(handles.Coord_z,'Value',coord_val(3))
684    set(handles.Coord_z,'String',ListCoord)
685    set(handles.Coord_z,'Visible','on')
686    set(handles.Check3D,'Value', 1)
687end
688
689%% set list of time coordinates
690menu=get(handles.SwitchVarIndexTime,'String');
691TimeOption=menu{get(handles.SwitchVarIndexTime,'Value')};
692switch TimeOption
693    case 'variable'
694        if numel(find(test_coord))<3
695            ListTime={''};
696        else
697            ListTime=Field.Display.ListVarName(find(test_coord,end));
698        end
699        set(handles.TimeName,'Value',1)
700        set(handles.TimeName,'String',ListTime)
701    case 'dim index'
702        if numel(find(test_coord))<3
703            ListTime={''};
704        else
705            ListTime=Field.Display.VarDimName{find(test_coord,end)};
706        end
707        set(handles.TimeName,'Value',1)
708        set(handles.TimeName,'String',ListTime)
709end
710if ~ischar(VarName)
711    update_field(handles,ScalarName)
712end
713
714% --- Executes on button press in check_rgb.
715function check_rgb_Callback(hObject, eventdata, handles)
716
717
718%------------------------------------------------------------------------
719% --- Executes on selection change in vector_x.
720%------------------------------------------------------------------------
721function vector_x_Callback(hObject, DimCell, handles)
722
723vector_x_menu=get(handles.vector_x,'String');
724vector_x_index=get(handles.vector_x,'Value');
725vector_x=vector_x_menu{vector_x_index};
726vector_Callback(handles)
727if ~ischar(DimCell)
728update_field(handles,vector_x)
729end
730
731%------------------------------------------------------------------------
732% --- Executes on selection change in vector_x.
733%------------------------------------------------------------------------
734function vector_y_Callback(hObject, DimCell, handles)
735
736vector_y_menu=get(handles.vector_x,'String');
737vector_y_index=get(handles.vector_x,'Value');
738vector_y=vector_y_menu{vector_y_index};
739vector_Callback(handles)
740if ~ischar(DimCell)
741update_field(handles,vector_y)
742end
743
744%------------------------------------------------------------------------
745% --- Executes on selection change in vector_z.
746function vector_z_Callback(hObject, DimCell, handles)
747%------------------------------------------------------------------------
748vector_z_menu=get(handles.vector_z,'String');
749vector_z_index=get(handles.vector_z,'Value');
750vector_z=vector_z_menu{vector_z_index};
751vector_Callback(handles)
752if ~ischar(DimCell)
753update_field(handles,vector_z)
754end
755%------------------------------------------------------------------------
756% --- Executes on selection change in vec_color.
757function vec_color_Callback(hObject, DimCell, handles)
758%------------------------------------------------------------------------
759index=get(handles.vec_color,'Value');
760string=get(handles.vec_color,'String');
761VarName=string{index};
762vector_Callback(handles)
763if ~ischar(DimCell)
764update_field(handles,VarName)
765end
766%------------------------------------------------------------------------
767% --- Executes on selection change in vector_x or vector_y
768function vector_Callback( handles)
769%------------------------------------------------------------------------
770Field=get(handles.get_field,'UserData');
771vector_x_index=get(handles.vector_x,'Value');
772vector_y_index=get(handles.vector_y,'Value');
773vec_color_index=get(handles.vec_color,'Value');
774
775%% set list of possible coordinates
776test_component=zeros(size(Field.Display.VarDimName));%=1 when variable #ilist is eligible as unstructured coordinate
777test_coord=zeros(size(Field.Display.VarDimName)); %=1 when variable #ilist is eligible as structured coordinate
778check_consistent=1;%check that the selected vector components (and possibly color var) have the same dimensiosn
779ListCoord={''};
780dim_var=Field.Display.VarDimName{vector_x_index};%list of dimensions of the selected variable
781if ~isequal(dim_var,Field.Display.VarDimName{vector_y_index})
782    check_consistent=0;
783elseif vec_color_index~=1 && ~isequal(dim_var,Field.Display.VarDimName{vec_color_index})
784    check_consistent=0;
785end
786% the two vector components have consistent dimensions
787if check_consistent
788    for ilist=1:numel(Field.Display.VarDimName)
789        dimnames=Field.Display.VarDimName{ilist}; %list of dimensions for variable #ilist
790        if isequal(dimnames,dim_var)
791            test_component(ilist)=1;
792        elseif numel(dimnames)==1 && ~isempty(find(strcmp(dimnames{1},dim_var)))%variable ilist is a 1D array which can be coordinate variable
793            test_coord(ilist)=1;
794        end
795    end
796    var_component=find(test_component);% list of variable indices elligible as unstructured coordinates
797    var_coord=find(test_coord);% % list of variable indices elligible as structured coordinates
798    var_component(var_component==vector_x_index|var_component==vector_y_index)=[];
799    var_coord(var_coord==vector_x_index|var_coord==vector_y_index)=[];% remove vector components form te possible list of coordinates
800    ListCoord=Field.Display.ListVarName([var_coord var_component]);
801   
802    %% set default coord selection
803    if numel(find(test_coord))>3
804        set(handles.SwitchVarIndexTime,'Value',3)% the last dim must be considered as time
805    end
806    if numel(var_component)<2 %unstructured coordinates excluded
807        if numel(find(test_coord))<2
808            ListCoord={''};
809        else
810            if numel(find(test_coord))>=3
811                set(handles.Coord_x,'Value',3)
812                set(handles.Coord_y,'Value',2)
813                set(handles.Coord_z,'Value',1)
814            else
815                set(handles.Coord_x,'Value',2)
816                set(handles.Coord_y,'Value',1)
817            end
818        end
819    else
820        coord_val=[0 0 0];
821        for ilist=1:numel(var_component)
822            ivar=var_component(ilist);
823            if isfield(Field.Display,'VarAttribute') && numel(Field.Display.VarAttribute)>=ivar && isfield(Field.Display.VarAttribute{ivar},'Role')
824                Role=Field.Display.VarAttribute{ivar}.Role;
825                if strcmp(Role,'coord_x')
826                    coord_val(1)=ilist;
827                elseif strcmp(Role,'coord_y')
828                    coord_val(2)=ilist;
829                elseif strcmp(Role,'coord_z')
830                    coord_val(3)=ilist;
831                end
832            end
833        end
834        if isempty(find(coord_val))
835            coord_val=var_coord;% case of dimension coordinates
836        end
837        if numel(find(coord_val))<2
838            coord_val=[1 2 3];
839        end
840        set(handles.Coord_x,'Value',coord_val(end))
841        set(handles.Coord_y,'Value',coord_val(end-1))
842        if numel(coord_val)>=3
843            set(handles.Coord_z,'Value',coord_val(end-2))
844        end
845    end
846end
847set(handles.Coord_z,'String',ListCoord)
848set(handles.Coord_y,'String',ListCoord)
849set(handles.Coord_x,'String',ListCoord)
850
851
852%% set list of time coordinates
853menu=get(handles.SwitchVarIndexTime,'String');
854TimeOption=menu{get(handles.SwitchVarIndexTime,'Value')};
855switch TimeOption
856    case 'variable'
857        if numel(find(test_coord))<3
858            ListTime={''};
859        else
860            ListTime=Field.Display.ListVarName(find(test_coord,end));
861        end
862        set(handles.TimeName,'Value',1)
863        set(handles.TimeName,'String',ListTime)
864    case 'dim index'
865        if numel(find(test_coord))<3
866            ListTime={''};
867        else
868            ListTime=Field.Display.VarDimName{find(test_coord,end)};
869        end
870        set(handles.TimeName,'Value',1)
871        set(handles.TimeName,'String',ListTime)
872end 
873
874%------------------------------------------------------------------------
875% --- Executes on selection change in SwitchVarIndexX.
876%------------------------------------------------------------------------
877function SwitchVarIndexX_Callback(hObject, eventdata, handles)
878
879%------------------------------------------------------------------------
880% --- Executes on selection change in Coord_x.
881%------------------------------------------------------------------------
882function Coord_x_Callback(hObject, DimCell, handles)
883DimCell
884index=get(handles.Coord_x,'Value');
885string=get(handles.Coord_x,'String');
886VarName=string{index};
887if ~ischar(DimCell)
888    update_field(handles,VarName)
889end
890if isequal(get(handles.FieldOption,'Value'),1)
891set_coord_y_options(handles,VarName)
892end
893
894%------------------------------------------------------------------------
895% --- Executes on selection change in Coord_y.
896%------------------------------------------------------------------------
897function Coord_y_Callback(hObject, DimCell, handles)
898
899index=get(handles.Coord_y,'Value');
900string=get(handles.Coord_y,'String');
901VarName=string{index};
902
903if ~ischar(DimCell)
904update_field(handles,VarName)
905end
906
907%------------------------------------------------------------------------
908% --- Executes on selection change in Coord_z.
909%------------------------------------------------------------------------
910function Coord_z_Callback(hObject, DimCell, handles)
911
912index=get(handles.Coord_z,'Value');
913string=get(handles.Coord_z,'String');
914VarName=string{index};
915if ~ischar(DimCell)
916update_field(handles,VarName)
917end
918
919%------------------------------------------------------------------------
920% --- Executes on selection change in SwitchVarIndexTime.
921%------------------------------------------------------------------------
922
923function SwitchVarIndexTime_Callback(hObject, eventdata, handles)
924
925Field=get(handles.get_field,'UserData');
926menu=get(handles.SwitchVarIndexTime,'String');
927option=menu{get(handles.SwitchVarIndexTime,'Value')};
928
929switch option
930    case 'file index'
931        set(handles.TimeName, 'Visible','off')% the time is taken as the file index
932    case 'attribute'
933        set(handles.TimeName, 'Visible','on')% timeName menu represents the available attributes
934        time_index=get(handles.TimeName,'UserData');    %select the input data
935        if isempty(time_index)
936            PreviousList=get(handles.TimeName, 'String');
937            if ~isempty(PreviousList)
938                PreviousAttr=PreviousList{get(handles.TimeName, 'Value')};
939                index=find(strcmp(PreviousAttr,Field.Display.ListGlobalAttribute),1);
940            end
941        end
942        if isempty(time_index)
943            time_index=find(~cellfun('isempty',regexp(Field.Display.ListGlobalAttribute,'Time')),1);% index of the attributes containing the string 'Time'
944        end     
945        if ~isempty(time_index)
946            set(handles.TimeName,'Value',time_index)
947        else
948            set(handles.TimeName,'Value',1)
949        end
950        set(handles.TimeName, 'String',Field.Display.ListGlobalAttribute)
951
952    case 'variable'% TimeName menu represents the available variables
953        set(handles.TimeName, 'Visible','on')
954        VarNbDim=cellfun('length',Field.Display.VarDimName); % check the nbre of dimensions of each input variable
955        TimeVarName=Field.Display.ListVarName(VarNbDim==1);% list of variables with a single dimension (candidate for time)
956        List=get(handles.TimeName,'String');% list of names on the menu for time
957        if isempty(List)
958            ind=1;
959        else
960            option=List{get(handles.TimeName,'Value')};% previous selected option
961            ind=find(strcmp(option,TimeVarName)); %check whether the previous selection is available in the newlist
962            if isempty(ind)
963                ind=1;
964            end
965        end
966        if ~isempty(TimeVarName)
967            set(handles.TimeName, 'Value',ind);% select first value in the menu if the option is not found
968            set(handles.TimeName, 'String',TimeVarName)% update the menu for time name
969        end
970    case 'matrix index'% TimeName menu represents the available dimensions
971        set(handles.TimeName, 'Visible','on')     
972        set(handles.TimeName, 'Value',1);
973        set(handles.TimeName, 'String',Field.Display.ListDimName)
974end
975TimeName_Callback(hObject, [], handles)
976
977%-----------------------------------------------------------------------
978% update the display of the variable 'VarName' and its dimensions in the list of variables
979function update_field(handles,VarName)
980%-----------------------------------------------------------------------
981Field=get(handles.get_field,'UserData');
982index=name2index(VarName,Field.ListVarName);
983if ~isempty(index)
984    set(handles.variables,'Value',index+1)
985    variables_Callback(handles.variables, VarName, handles)
986end
987
988%------------------------------------------------------------------------
989% --- give index numbers of the strings str in the list ListvarName
990% -----------------------------------------------------------------------
991function VarIndex_y=name2index(cell_str,ListVarName)
992
993VarIndex_y=[];
994if ischar(cell_str)
995    VarIndex_y=find(strcmp(cell_str,ListVarName),1);
996elseif iscell(cell_str)
997    for isel=1:length(cell_str)
998        varsel=cell_str{isel};
999        for ivar=1:length(ListVarName)
1000            varlist=ListVarName{ivar};
1001            if isequal(varlist,varsel)
1002                VarIndex_y=[VarIndex_y ivar];
1003            end
1004        end
1005    end
1006end
1007
1008
1009
1010% % --- Executes on button press in CheckDimensionY.
1011% function CheckDimensionY_Callback(hObject, eventdata, handles)
1012% FieldList=get(handles.FieldOption,'String');
1013% FieldOption=FieldList{get(handles.FieldOption,'Value')};
1014% switch FieldOption
1015%     case '1D plot'
1016%         
1017%     case {'scalar','pick variables'}
1018%        scalar_Callback(hObject, eventdata, handles)
1019%     case 'vectors'
1020% end
1021%
1022%
1023% % --- Executes on button press in CheckDimensionZ.
1024% function CheckDimensionZ_Callback(hObject, eventdata, handles)
1025% FieldList=get(handles.FieldOption,'String');
1026% FieldOption=FieldList{get(handles.FieldOption,'Value')};
1027% switch FieldOption
1028%     case '1D plot'
1029%         
1030%     case 'scalar'
1031%        scalar_Callback(hObject, eventdata, handles)
1032%     case 'vectors'
1033% end
1034
1035% --- Executes on selection change in TimeName.
1036function TimeName_Callback(hObject, eventdata, handles)
1037Field=get(handles.get_field,'UserData');
1038index=get(handles.SwitchVarIndexTime,'Value');
1039MenuIndex=get(handles.TimeName,'Value');
1040string=get(handles.TimeName,'String');
1041TimeName='';%default
1042if ~isempty(string)&&iscell(string)
1043TimeName=string{MenuIndex};
1044end
1045switch index
1046    case 1
1047        set(handles.num_TimeDimension,'String','')
1048        set(handles.TimeUnit,'String','index')
1049    case 2
1050        set(handles.num_TimeDimension,'String','')
1051        attr_index=find(strcmpi([TimeName 'Unit'],Field.ListGlobalAttribute));% look for time unit
1052        if ~isempty(attr_index)
1053            AttrName=Field.ListGlobalAttribute{attr_index};
1054            set(handles.TimeUnit,'String',Field.(AttrName))
1055        else
1056            set(handles.TimeUnit,'String','')
1057        end
1058    case {3 ,4}
1059        if index==3  % TimeName is used to chose a variable
1060            VarIndex=name2index(TimeName,Field.ListVarName);
1061            DimName=Field.VarDimName{VarIndex};
1062            DimIndex=name2index(DimName,Field.ListDimName);
1063            DimValue=Field.DimValue(DimIndex);
1064            set(handles.num_TimeDimension,'String',num2str(DimValue))
1065            unit='';
1066            if isfield(Field,'VarAttribute')&& isfield(Field.VarAttribute{VarIndex},'Unit')
1067                unit=Field.VarAttribute{VarIndex}.Unit;
1068            end
1069            set(handles.TimeUnit,'String',unit)
1070            update_field(handles,TimeName)
1071        elseif index==4% TimeName is used to chose a dimension
1072            DimName=string{MenuIndex};
1073            DimIndex=name2index(DimName,Field.ListDimName);
1074            DimValue=Field.DimValue(DimIndex);
1075            set(handles.num_TimeDimension,'String',num2str(DimValue))
1076            set(handles.TimeUnit,'String','index')
1077        end
1078end
1079
1080%-----------------------------------------------------------------------
1081% --- Executes on button press in Check3D.
1082%-----------------------------------------------------------------------
1083function Check3D_Callback(hObject, eventdata, handles)
1084if get(handles.Check3D,'Value')% 3D fields
1085    status='on';
1086else% fields studied as 2D
1087    status='off';
1088end
1089
1090set(handles.Coord_z,'Visible',status)
1091% set(handles.CheckDimensionZ,'Visible',status)
1092set(handles.Z_title,'Visible',status)
1093set(handles.vector_z,'Visible',status)
1094set(handles.W_title,'Visible',status)
1095Field=get(handles.get_field,'UserData');
1096if strcmp(status,'on')% ask for 3D input       
1097    if Field.MaxDim>3% for 4D fields, propose to use the fourth variable as time
1098        %set(handles.Time,'Visible','on')
1099        menu=get(handles.SwitchVarIndexTime,'String');
1100        val=find(strcmp('variable',menu));
1101        if ~isempty(val)
1102            set(handles.SwitchVarIndexTime,'Value',val)
1103        end
1104    else
1105        set(handles.SwitchVarIndexTime,'Value',1)
1106        set(handles.SwitchVarIndexTime,'String',{'file index';'attribute'})
1107    end
1108else
1109   set(handles.SwitchVarIndexTime,'String',get(handles.SwitchVarIndexTime,'UserData'))
1110   if Field.MaxDim >=3
1111       var_index=find(strcmp('variable',get(handles.SwitchVarIndexTime,'UserData')));
1112       set(handles.SwitchVarIndexTime,'Value',var_index)
1113   end
1114end
1115SwitchVarIndexTime_Callback(handles.SwitchVarIndexTime,[], handles)
1116
1117%------------------------------------------------------------------------
1118% --- Executes on button press in OK.
1119%------------------------------------------------------------------------
1120function OK_Callback(hObject, eventdata, handles)
1121handles.output=read_GUI(handles.get_field);
1122guidata(hObject, handles);% Update handles structure
1123uiresume(handles.get_field);
1124drawnow
1125% this function then activate get_field_OutputFcn
1126
1127%------------------------------------------------------------------------
1128% --- Executes when the GUI is closed by the mouse on upper right corner.
1129%------------------------------------------------------------------------
1130function closefcn(hObject, eventdata, handles)
1131handles.output=[];
1132guidata(hObject, handles);% Update handles structure
1133uiresume(handles.get_field);
1134drawnow
1135
1136%------------------------------------------------------------------------
1137% --- Outputs from this function are returned to the command line.
1138%------------------------------------------------------------------------
1139function varargout = get_field_OutputFcn(hObject, eventdata, handles)
1140
1141varargout{1} =handles.output;
1142delete(handles.get_field)
Note: See TracBrowser for help on using the repository browser.