source: trunk/src/get_field.m @ 1082

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

rodrigues added

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