source: trunk/src/get_field.m @ 693

Last change on this file since 693 was 693, checked in by sommeria, 11 years ago

bug repaired geometry_calib/Detect grid, + corrections in get_field

File size: 39.0 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%
[674]6% OUPUT:
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%
[204]18%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
19%  Copyright Joel Sommeria, 2008, LEGI / CNRS-UJF-INPG, sommeria@coriolis-legi.org.
20%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
21%     This file is part of the toolbox UVMAT.
22%
23%     UVMAT is free software; you can redistribute it and/or modify
24%     it under the terms of the GNU General Public License as published by
25%     the Free Software Foundation; either version 2 of the License, or
26%     (at your option) any later version.
27%
28%     UVMAT is distributed in the hope that it will be useful,
29%     but WITHOUT ANY WARRANTY; without even the implied warranty of
30%     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
31%     GNU General Public License (file UVMAT/COPYING.txt) for more details.
32%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
33
34function varargout = get_field(varargin)
35
[693]36% Last Modified by GUIDE v2.5 23-Oct-2013 23:44:17
[204]37
38% Begin initialization code - DO NOT EDIT
[581]39gui_Singleton = 1;
[204]40gui_State = struct('gui_Name',       mfilename, ...
41                   'gui_Singleton',  gui_Singleton, ...
42                   'gui_OpeningFcn', @get_field_OpeningFcn, ...
43                   'gui_OutputFcn',  @get_field_OutputFcn, ...
44                   'gui_LayoutFcn',  [] , ...
45                   'gui_Callback',   []);
46if nargin && ischar(varargin{1})&& ~isempty(regexp(varargin{1},'_Callback','once'))
47    gui_State.gui_Callback = str2func(varargin{1});
48end
49
50if nargout
51    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
52else
53    gui_mainfcn(gui_State, varargin{:});
54end
55% End initialization code - DO NOT EDIT
[674]56     
[204]57%------------------------------------------------------------------------
58% --- Executes just before get_field is made visible.
[648]59%------------------------------------------------------------------------
[581]60function get_field_OpeningFcn(hObject, eventdata, handles,filename,ParamIn)
[204]61
[648]62%% GUI settings
[581]63handles.output = 'Cancel';
[204]64guidata(hObject, handles);
[434]65set(hObject,'WindowButtonDownFcn',{'mouse_down'}) % allows mouse action with right button (zoom for uicontrol display)
[204]66
[648]67%% enter input data
[674]68if ischar(filename)% input file name
69    set(handles.inputfile,'String',filename)% fill the input file name
70    Field=nc2struct(filename,[]);% reads the  field structure, without the variables
71else
72    'bad input to get_field'
73end
74if ~exist('ParamIn','var')
75    ParamIn=[];
76end
[648]77if isfield(Field,'Txt')
[674]78    msgbox_uvmat('ERROR',['get_field/nc2struct/' Field.Txt])% display error message for input file reading
[648]79    return
[227]80end
[204]81if ~isfield(Field,'ListVarName')
82    return
83end
[582]84
[648]85%% look at singletons and variables with a single dimension
86Field.Display=Field;
87Field.Check0D=zeros(size(Field.ListVarName));% =1 for arrays with a single value
88NbVar=numel(Field.VarDimName);
89for ilist=1:NbVar
[582]90    if ischar(Field.VarDimName{ilist})
91        Field.VarDimName{ilist}={Field.VarDimName{ilist}}; %transform string into cell
92    end
[648]93    NbDim=numel(Field.VarDimName{ilist});
[582]94    check_singleton=false(1,NbDim);
95    for idim=1:NbDim
[648]96        dim_index=strcmp(Field.VarDimName{ilist}{idim},Field.ListDimName);%index in the list of dimensions
97        check_singleton(idim)=isequal(Field.DimValue(dim_index),1);%check_singleton=1 for singleton
[582]98    end
[693]99    Field.Check0D(ilist)=(isequal(check_singleton,ones(1,NbDim)))||(~isequal(Field.VarType(ilist),4)&&~isequal(Field.VarType(ilist),5));% =1 if the variable reduces to a single value
[648]100    if ~Field.Check0D(ilist)
101    Field.Display.VarDimName{ilist}=Field.VarDimName{ilist}(~check_singleton);% eliminate singletons in the list of variable dimensions
[582]102    end
[648]103    %Field.NbDim(ilist)=numel(Field.VarDimNameNonSingleton{ilist});%nbre of array dimensions after elimination of singletons
[582]104end
[669]105if ~isfield(Field,'VarAttribute')
106    Field.VarAttribute={};
107end
[648]108if numel(Field.VarAttribute)<NbVar% complement VarAttribute by blanjs if neded
109    Field.VarAttribute(numel(Field.VarAttribute)+1:NbVar)=cell(1,NbVar-numel(Field.VarAttribute));
[582]110end
[672]111% Field.Display = list of variables and corresponding properties obtained after removal of singletons
[648]112Field.Display.ListVarName=Field.ListVarName(~Field.Check0D);
113Field.Display.VarAttribute=Field.VarAttribute(~Field.Check0D);
114Field.Display.VarDimName=Field.Display.VarDimName(~Field.Check0D);
115Field.Display.DimValue=Field.DimValue(Field.DimValue~=1);
116Field.Display.ListDimName=Field.ListDimName(Field.DimValue~=1);
[581]117
[648]118%% analyse the input field cells
119[CellInfo,NbDim,errormsg]=find_field_cells(Field.Display);
120if ~isempty(errormsg)
121    msgbox_uvmat('ERROR',['get_field / Field_input / find_field_cells: ' errormsg])
122    return
123end
124[Field.MaxDim,imax]=max(NbDim);
125
[581]126%% set time mode
[646]127ListSwitchVarIndexTime={'file index'};% default setting: the time is the file index
[582]128% look at global attributes with numerical values
[674]129check_numvalue=false(1,numel(Field.ListGlobalAttribute));
[646]130for ilist=1:numel(Field.ListGlobalAttribute)
131    Value=Field.(Field.ListGlobalAttribute{ilist});
132    check_numvalue(ilist)=isnumeric(Value);
133end
[648]134Field.Display.ListGlobalAttribute=Field.ListGlobalAttribute(check_numvalue);% select the attributes with float numerical value
135if ~isempty(Field.Display.ListGlobalAttribute)
[646]136    ListSwitchVarIndexTime=[ListSwitchVarIndexTime; {'attribute'}];% the time can be chosen as a global attribute
[582]137end
138if Field.MaxDim>=2
[693]139    ListSwitchVarIndexTime=[ListSwitchVarIndexTime;{'variable'};{'matrix index'}];% the time can be chosen as a dim index
[582]140end
[674]141
142%% select the Time attribute from input
143if isfield(ParamIn,'TimeAttrName')
144    time_index=find(strcmp(ParamIn.TimeAttrName,Field.Display.ListGlobalAttribute),1);
145else
[693]146    time_index=find(~cellfun('isempty',regexp(Field.Display.ListGlobalAttribute,'Time')),1);% look for global attribute containing name 'Time'
[674]147end
148if ~isempty(time_index)
[648]149    set(handles.SwitchVarIndexTime,'Value',2);
[674]150    set(handles.TimeName,'UserData',time_index)
[581]151else
[648]152    set(handles.SwitchVarIndexTime,'Value',1);
[581]153end
[646]154set(handles.SwitchVarIndexTime,'String',ListSwitchVarIndexTime)
[582]155set(handles.get_field,'UserData',Field);% record the finput field structure
[648]156SwitchVarIndexTime_Callback([], [], handles)
[581]157
[582]158%% set vector menu (priority) if detected or scalar menu for space dim >=2, or usual (x,y) plot for 1D fields
[648]159set(handles.vector_x,'String',Field.Display.ListVarName)% fill the menu of x vector components
160set(handles.vector_y,'String',Field.Display.ListVarName)% fill the menu of y vector components
[672]161set(handles.vector_z,'String',[{''} Field.Display.ListVarName])% fill the menu of y vector components
162set(handles.vec_color,'String',[{''} Field.Display.ListVarName])% fill the menu of y vector components
[648]163set(handles.scalar,'Value',1)% fill the menu of y vector components
[672]164set(handles.scalar,'String',Field.Display.ListVarName)% fill the menu of y vector components
[648]165set(handles.ordinate,'Value',1)% fill the menu of y vector components
[672]166set(handles.ordinate,'String',Field.Display.ListVarName)% fill the menu of y vector components
[654]167if isfield(Field,'Conventions')&& strcmp(Field.Conventions,'uvmat/civdata')
[672]168    set(handles.FieldOption,'String',{'1D plot';'scalar';'vectors';'civdata...'})
[654]169else
[672]170    set(handles.FieldOption,'String',{'1D plot';'scalar';'vectors'})
[654]171end
[582]172if Field.MaxDim>=2 % case of 2D (or 3D) fields
[674]173    check_vec_input=0;
174    if isfield(ParamIn,'vector_x')&& isfield(ParamIn,'vector_y')
175        ichoice_x=find(strcmp(ParamIn.vector_x,Field.Display.ListVarName),1);
176        ichoice_y=find(strcmp(ParamIn.vector_y,Field.Display.ListVarName),1);
177        if ~isempty(ichoice_x)&&~isempty(ichoice_y)
178            set(handles.vector_x,'UserData',ichoice_x)
179            set(handles.vector_y,'UserData',ichoice_y)
180            check_vec_input=1;
181        end
182    end
183    if ~check_vec_input && isfield(CellInfo{imax},'VarIndex_vector_x') &&  isfield(CellInfo{imax},'VarIndex_vector_y')
184        set(handles.vector_x,'UserData',CellInfo{imax}.VarIndex_vector_x(1))
185        set(handles.vector_y,'UserData',CellInfo{imax}.VarIndex_vector_y(1))
186        check_vec_input=1;
187    end
188    if check_vec_input
[646]189        set(handles.FieldOption,'Value',3)% set vector selection option
[204]190    else
[674]191        set(handles.FieldOption,'Value',2)% set scalar selection option
[204]192    end
[582]193else % case of 1D fields
[644]194    set(handles.FieldOption,'Value',1)
[204]195end
[582]196
[648]197%% fill the general list of dimensions, variables, attributes
198if isfield(Field,'ListDimName')&&~isempty(Field.ListDimName)
199    Tabcell(:,1)=Field.ListDimName;
200    for iline=1:length(Field.ListDimName)
201        Tabcell{iline,2}=num2str(Field.DimValue(iline));
202    end
203    Tabchar=cell2tab(Tabcell,' = ');
204    set(handles.dimensions,'String',Tabchar)
205end
206set(handles.variables,'Value',1)
207set(handles.variables,'String',[{'*'} Field.ListVarName])
208variables_Callback(handles.variables,[], handles)% list the global attributes
[644]209
[648]210%% fill menus for coordinates and time
211FieldOption_Callback(handles.variables,[], handles)% list the global attributes
212
[672]213%% Make choices of coordinates from input
214if isfield(CellInfo{imax},'CoordIndex')
215    CoordIndex=CellInfo{imax}.CoordIndex;
216    if numel(CoordIndex)==2
217        YName=Field.ListVarName{CoordIndex(1)};
218        XName=Field.ListVarName{CoordIndex(2)};
219        ListCoord=get(handles.Coord_x,'String');
220        XIndex=find(strcmp(XName,ListCoord));
221        if ~isempty(XIndex)
222            set(handles.Coord_x,'Value',XIndex)
223        end
224        YIndex=find(strcmp(YName,ListCoord));
225        if ~isempty(YIndex)
226            set(handles.Coord_y,'Value',YIndex)
227        end
228    end
229end
[648]230
231%% put the GUI on the lower right of the sceen
232set(hObject,'Unit','pixel')
233pos_view_field=get(hObject,'Position');
234set(0,'Unit','pixels')
235ScreenSize=get(0,'ScreenSize');
236pos_view_field(1)=ScreenSize(1)+ScreenSize(3)-pos_view_field(3);
237pos_view_field(2)=ScreenSize(2);
238set(hObject,'Position',pos_view_field)
239set(handles.get_field,'WindowStyle','modal')% Make the GUI modal
240drawnow
241uiwait(handles.get_field);
242
243return
244
245%% set z coordinate menu if relevant
246if Field.MaxDim>=3
247    set(handles.vector_z,'Visible','on')
248    set(handles.vector_z,'String',Field.ListVarName(~Field.Check0D))
249    set(handles.Coord_z,'Visible','on')
250    set(handles.SwitchVarIndexZ,'Visible','on')
251    set(handles.Z_title,'Visible','on')
252else
253    set(handles.vector_z,'Visible','off')
254    set(handles.Coord_z,'Visible','off')
255    set(handles.Z_title,'Visible','off')
[581]256end
[204]257
[674]258%% make selections according to ParamIn
259if isfield(ParamIn,'vector_x') && isfield(ParamIn,'vector_y')
260   
261end
262
[204]263%------------------------------------------------------------------------
[672]264% --- Outputs from this function are returned to the command line.
[648]265%------------------------------------------------------------------------
[582]266function varargout = get_field_OutputFcn(hObject, eventdata, handles)
[672]267
[582]268varargout{1} = handles.output;
269delete(handles.get_field)
270
[648]271%------------------------------------------------------------------------
272% --- Executes when user attempts to close get_field.
273%------------------------------------------------------------------------
274function get_field_CloseRequestFcn(hObject, eventdata, handles)
[582]275
[648]276if isequal(get(handles.get_field, 'waitstatus'), 'waiting')
277    % The GUI is still in UIWAIT, us UIRESUME
278    uiresume(handles.get_field);
279else
280    % The GUI is no longer waiting, just close it
281    delete(handles.get_field);
282end
283
[672]284%------------------------------------------------------------------------
[648]285% --- Executes on button press in OK.
286%------------------------------------------------------------------------
287function OK_Callback(hObject, eventdata, handles)
288
289handles.output=read_GUI(handles.get_field);
290guidata(hObject, handles);% Update handles structure
291uiresume(handles.get_field);
292drawnow
293
[582]294% -----------------------------------------------------------------------
295% --- Activated by selection in the list of variables
[672]296% ----------------------------------------------------------------------
[582]297function variables_Callback(hObject, eventdata, 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');
325    var_select=list_var{index};
326    set(handles.attributes_txt,'String', ['attributes of ' var_select])
327    if isfield(Field,'VarAttribute')&& length(Field.VarAttribute)>=index-1
[672]328        %         nbline=0;
[582]329        VarAttr=Field.VarAttribute{index-1};
330        if isstruct(VarAttr)
331            attr_list=fieldnames(VarAttr);
332            for iline=1:length(attr_list)
333                Tabcell{iline,1}=attr_list{iline};
[674]334                val=VarAttr.(attr_list{iline}) ;
[582]335                if ischar(val);
336                    Tabcell{iline,2}=val;
337                else
[672]338                    Tabcell{iline,2}=num2str(val);
[582]339                end
340            end
341        end
342    end
[672]343   
[582]344end
345if ~isempty(Tabcell)
346    Tabchar=cell2tab(Tabcell,'=');
[672]347    %     Tabchar=[{''};Tabchar];
[582]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
356        dim_indices=1:length(Field.ListDimName);
357        set(handles.dimensions_txt,'String', 'dimensions')
358    else
359        DimCell=Field.VarDimName{index-1};
360        if ischar(DimCell)
361            DimCell={DimCell};
[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
369        set(handles.dimensions_txt,'String', ['dimensions of ' var_select])
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
381%------------------------------------------------------------------------
[648]382% --- Executes on selection change in FieldOption.
[582]383%------------------------------------------------------------------------
[648]384function FieldOption_Callback(hObject, eventdata, handles)
385
[672]386Field=get(handles.get_field,'UserData');
[648]387FieldList=get(handles.FieldOption,'String');
388FieldOption=FieldList{get(handles.FieldOption,'Value')};
389switch FieldOption
[672]390   
[648]391    case '1D plot'
[654]392        set(handles.Coordinates,'Visible','on')
[648]393        set(handles.PanelOrdinate,'Visible','on')
394        pos=get(handles.PanelOrdinate,'Position');
395        pos(1)=2;
396        pos_coord=get(handles.Coordinates,'Position');
397        pos(2)=pos_coord(2)-pos(4)-2;
398        set(handles.PanelOrdinate,'Position',pos)
399        set(handles.PanelScalar,'Visible','off')
400        set(handles.PanelVectors,'Visible','off')
401        set(handles.Coord_y,'Visible','off')
402        set(handles.Y_title,'Visible','off')
403        set(handles.Coord_z,'Visible','off')
404        set(handles.Z_title,'Visible','off')
405        ordinate_Callback(hObject, eventdata, handles)
[672]406       
[648]407    case 'scalar'
[654]408        set(handles.Coordinates,'Visible','on')
[648]409        set(handles.PanelOrdinate,'Visible','off')
410        set(handles.PanelScalar,'Visible','on')
411        set(handles.PanelVectors,'Visible','off')
412        pos=get(handles.PanelScalar,'Position');
413        pos(1)=2;
414        pos_coord=get(handles.Coordinates,'Position');
415        pos(2)=pos_coord(2)-pos(4)-2;
416        set(handles.PanelScalar,'Position',pos)
417        set(handles.Coord_y,'Visible','on')
418        set(handles.Y_title,'Visible','on')
[672]419        %default scalar selection
420        test_coord=zeros(size(Field.Display.VarDimName)); %=1 when variable #ilist is eligible as structured coordiante
421        for ilist=1:numel(Field.Display.VarDimName)
422            if isfield(Field.Display,'VarAttribute') && numel(Field.Display.VarAttribute)>=ilist && isfield(Field.Display.VarAttribute{ilist},'Role')
423                Role=Field.Display.VarAttribute{ilist}.Role;
424                if strcmp(Role,'coord_x')||strcmp(Role,'coord_y')
425                    test_coord(ilist)=1;
426                end
427            end
428            dimnames=Field.Display.VarDimName{ilist}; %list of dimensions for variable #ilist
429            if numel(dimnames)==1 && strcmp(dimnames{1},Field.Display.ListVarName{ilist})%dimension variable
430                test_coord(ilist)=1;
431            end
432        end
433        scalar_index=find(~test_coord,1);%get the first variable not a coordiante
434        if isempty(scalar_index)
435            set(handles.scalar,'Value',1)
436        else
437            set(handles.scalar,'Value',scalar_index)
438        end       
[648]439        scalar_Callback(hObject, eventdata, handles)
[674]440             
[648]441    case 'vectors'
[674]442        set(handles.PanelVectors,'Visible','on')
[654]443        set(handles.Coordinates,'Visible','on')
[648]444        set(handles.PanelOrdinate,'Visible','off')
445        set(handles.PanelScalar,'Visible','off')
446        pos=get(handles.PanelVectors,'Position');
447        pos(1)=2;
448        pos_coord=get(handles.Coordinates,'Position');
449        pos(2)=pos_coord(2)-pos(4)-2;
450        set(handles.PanelVectors,'Position',pos)
451        set(handles.Coord_y,'Visible','on')
[672]452        set(handles.Y_title,'Visible','on')
453        %default vector selection
[674]454        vector_x_value=get(handles.vector_x,'UserData');
455        vector_y_value=get(handles.vector_y,'UserData');
456        if ~isempty(vector_x_value)&&~isempty(vector_y_value)
457            set(handles.vector_x,'Value',vector_x_value)
458            set(handles.vector_y,'Value',vector_y_value)
459        else
460            test_coord=zeros(size(Field.Display.VarDimName)); %=1 when variable #ilist is eligible as structured coordinate
461            for ilist=1:numel(Field.Display.VarDimName)
462                if isfield(Field.Display,'VarAttribute') && numel(Field.Display.VarAttribute)>=ilist && isfield(Field.Display.VarAttribute{ilist},'Role')
463                    Role=Field.Display.VarAttribute{ilist}.Role;
464                    if strcmp(Role,'coord_x')||strcmp(Role,'coord_y')
465                        test_coord(ilist)=1;
466                    end
467                end
468                dimnames=Field.Display.VarDimName{ilist}; %list of dimensions for variable #ilist
469                if numel(dimnames)==1 && strcmp(dimnames{1},Field.Display.ListVarName{ilist})%dimension variable
[672]470                    test_coord(ilist)=1;
471                end
472            end
[674]473            vector_index=find(~test_coord,2);%get the two first variables not a coordinate
474            if isempty(vector_index)
475                set(handles.vector_x,'Value',1)
476                set(handles.vector_y,'Value',2)
477            else
478                set(handles.vector_x,'Value',vector_index(1))
479                set(handles.vector_y,'Value',vector_index(2))
[672]480            end
481        end
[648]482        vector_Callback(handles)
[672]483       
[654]484    case 'civdata...'
485        set(handles.PanelOrdinate,'Visible','off')
486        set(handles.PanelScalar,'Visible','off')
487        set(handles.PanelVectors,'Visible','off')
488        set(handles.Coordinates,'Visible','off')
[582]489end
490
491%------------------------------------------------------------------------
[204]492function ordinate_Callback(hObject, eventdata, handles)
493%------------------------------------------------------------------------
[582]494Field=get(handles.get_field,'UserData');
[648]495y_index=get(handles.ordinate,'Value');
496y_menu=get(handles.ordinate,'String');
497YName=y_menu{y_index};
498
499%% set list of possible coordinates
500test_component=zeros(size(Field.Display.VarDimName));%=1 when variable #ilist is eligible as unstructured coordinate
501test_coord=zeros(size(Field.Display.VarDimName)); %=1 when variable #ilist is eligible as structured coordiante
502ListCoord={''};
503dim_var=Field.Display.VarDimName{y_index};%list of dimensions of the selected variable
504
505for ilist=1:numel(Field.Display.VarDimName)
506    dimnames=Field.Display.VarDimName{ilist}; %list of dimensions for variable #ilist
507    if isequal(dimnames,dim_var)
508        test_component(ilist)=1;
509    elseif numel(dimnames)==1 && ~isempty(find(strcmp(dimnames{1},dim_var)))%variable ilist is a 1D array which can be coordinate variable
510        test_coord(ilist)=1;
[582]511    end
512end
[648]513var_component=find(test_component);% list of variable indices elligible as unstructured coordinates
514var_coord=find(test_coord);% % list of variable indices elligible as structured coordinates
515ListCoord=Field.Display.ListVarName([var_component var_coord]);
[204]516
[648]517%% set default coord selection
518if numel(find(test_coord))>3
[693]519     SwitchVarIndexTime=get(handles.SwitchVarIndexTime,'String');
520    if numel(SwitchVarIndexTime)<3
521        SwitchVarIndexTime=[SwitchVarIndexTime;'matrix_index'];
522        set(handles.SwitchVarIndexTime,'String',SwitchVarIndexTime)
523    end
[648]524    set(handles.SwitchVarIndexTime,'Value',3)% the last dim must be considered as time
[693]525    SwitchVarIndexTime_Callback([], [], handles)
[582]526end
[648]527if numel(var_component)<2
528    if numel(test_coord)<2
529        ListCoord={''};
530    else
531        set(handles.Coord_x,'Value',2)
532        set(handles.Coord_y,'Value',1)
533    end
[582]534else
[648]535    coord_val=1;
536    for ilist=1:numel(var_component)
537        ivar=var_component(ilist);
538        if isfield(Field.Display,'VarAttribute') && numel(Field.Display.VarAttribute)>=ivar && isfield(Field.Display.VarAttribute{ivar},'Role')
539            Role=Field.Display.VarAttribute{ivar}.Role;
540            if strcmp(Role,'coord_x')
541                coord_val=ilist;
542            end
543        end
544    end
545    set(handles.Coord_x,'Value',coord_val)
[582]546end
[648]547set(handles.Coord_x,'String',ListCoord)
[582]548
549
[648]550%% set list of time coordinates
551menu=get(handles.SwitchVarIndexTime,'String');
552TimeOption=menu{get(handles.SwitchVarIndexTime,'Value')};
553switch TimeOption
554    case 'variable'
555        if numel(find(test_coord))<3
556            ListTime={''};
557        else
558            ListTime=Field.Display.ListVarName(find(test_coord,end));
559        end
560        set(handles.TimeName,'Value',1)
561        set(handles.TimeName,'String',ListTime)
[693]562    case 'matrix index'
[648]563        if numel(find(test_coord))<3
564            ListTime={''};
565        else
566            ListTime=Field.Display.VarDimName{find(test_coord,end)};
567        end
568        set(handles.TimeName,'Value',1)
569        set(handles.TimeName,'String',ListTime)
570end 
571update_field(handles,YName)
[674]572         
[582]573%------------------------------------------------------------------------
[204]574% --- Executes on selection change in scalar menu.
[648]575%------------------------------------------------------------------------
[204]576function scalar_Callback(hObject, eventdata, handles)
[648]577
[586]578Field=get(handles.get_field,'UserData');
[648]579scalar_index=get(handles.scalar,'Value');
580scalar_menu=get(handles.scalar,'String');
581ScalarName=scalar_menu{scalar_index};
[204]582
[648]583%% set list of possible coordinates
584test_component=zeros(size(Field.Display.VarDimName));%=1 when variable #ilist is eligible as unstructured coordinate
585test_coord=zeros(size(Field.Display.VarDimName)); %=1 when variable #ilist is eligible as structured coordiante
586dim_var=Field.Display.VarDimName{scalar_index};%list of dimensions of the selected variable
[693]587CheckDimensionX=get(handles.CheckDimensionX,'Value');
588CheckDimensionY=get(handles.CheckDimensionY,'Value');
589if ~CheckDimensionX  || ~CheckDimensionY
590    %look for coordinate variables among the other variables
591    for ilist=1:numel(Field.Display.VarDimName)
592        dimnames=Field.Display.VarDimName{ilist}; %list of dimensions for variable #ilist
593        if isequal(dimnames,dim_var)
594            test_component(ilist)=1;% the listed variable has the same dimension as the selected scalar-> possibly chosen as unstructured coordinate
595        elseif numel(dimnames)==1 && ~isempty(find(strcmp(dimnames{1},dim_var), 1))%variable ilist is a 1D array which can be coordinate variable
596            test_coord(ilist)=1;
597        end
[648]598    end
[227]599end
[672]600var_component=find(test_component);% list of variable indices elligible as unstructured coordinates
[648]601var_coord=find(test_coord);% % list of variable indices elligible as structured coordinates
602ListCoord=Field.Display.ListVarName([var_component var_coord]);
[227]603
[648]604%% set default coord selection
605if numel(find(test_coord))>3
[693]606    SwitchVarIndexTime=get(handles.SwitchVarIndexTime,'String');
607    if numel(SwitchVarIndexTime)<3
608        SwitchVarIndexTime=[SwitchVarIndexTime;'matrix_index'];
609        set(handles.SwitchVarIndexTime,'String',SwitchVarIndexTime)
610    end
[648]611    set(handles.SwitchVarIndexTime,'Value',3)% the last dim must be considered as time
[693]612    SwitchVarIndexTime_Callback([], [], handles)
[648]613end
[672]614% if numel(var_component)<2
615%     if numel(test_coord)<2
616%         ListCoord={''};
617%     else
618%         set(handles.Coord_x,'Value',2)
619%         set(handles.Coord_y,'Value',1)
620%     end
621% else
622coord_val=[0 0];
623% look for labelled unstructured coordinates
624for ilist=1:numel(var_component)
625    ivar=var_component(ilist);
626    if isfield(Field.Display,'VarAttribute') && numel(Field.Display.VarAttribute)>=ivar && isfield(Field.Display.VarAttribute{ivar},'Role')
627        Role=Field.Display.VarAttribute{ivar}.Role;
628        if strcmp(Role,'coord_x')
629            coord_val(1)=ilist;
630        elseif strcmp(Role,'coord_y')
631            coord_val(2)=ilist;
[227]632        end
633    end
634end
[672]635if numel(find(coord_val))<2
636    if numel(var_coord)>=2
637        coord_val=[numel(var_component)+2 numel(var_component)+1];
638    else
639        coord_val=[1 2];
640    end
641end
[693]642if  CheckDimensionX
643    set(handles.Coord_x,'Value',2)
644    set(handles.Coord_x,'String',dim_var')
645else
646    set(handles.Coord_x,'Value',coord_val(1))
647    set(handles.Coord_x,'String',ListCoord)
648end
649if  CheckDimensionY
650    set(handles.Coord_y,'Value',1)
651    set(handles.Coord_y,'String',dim_var')
652else
653    set(handles.Coord_y,'Value',coord_val(2))
654    set(handles.Coord_y,'String',ListCoord)
655end
[227]656
[648]657%% set list of time coordinates
658menu=get(handles.SwitchVarIndexTime,'String');
659TimeOption=menu{get(handles.SwitchVarIndexTime,'Value')};
660switch TimeOption
661    case 'variable'
662        if numel(find(test_coord))<3
663            ListTime={''};
664        else
665            ListTime=Field.Display.ListVarName(find(test_coord,end));
666        end
667        set(handles.TimeName,'Value',1)
668        set(handles.TimeName,'String',ListTime)
669    case 'dim index'
670        if numel(find(test_coord))<3
671            ListTime={''};
672        else
673            ListTime=Field.Display.VarDimName{find(test_coord,end)};
674        end
675        set(handles.TimeName,'Value',1)
676        set(handles.TimeName,'String',ListTime)
677end 
678update_field(handles,ScalarName)
[582]679
[648]680% --- Executes on button press in check_rgb.
681function check_rgb_Callback(hObject, eventdata, handles)
[582]682
683
684%------------------------------------------------------------------------
[204]685% --- Executes on selection change in vector_x.
[672]686%------------------------------------------------------------------------
[204]687function vector_x_Callback(hObject, eventdata, handles)
[672]688
[648]689vector_x_menu=get(handles.vector_x,'String');
690vector_x_index=get(handles.vector_x,'Value');
691vector_x=vector_x_menu{vector_x_index};
692vector_Callback(handles)
693update_field(handles,vector_x)
[227]694
[204]695%------------------------------------------------------------------------
[648]696% --- Executes on selection change in vector_x.
[204]697function vector_y_Callback(hObject, eventdata, handles)
698%------------------------------------------------------------------------
[648]699vector_y_menu=get(handles.vector_x,'String');
700vector_y_index=get(handles.vector_x,'Value');
701vector_y=vector_y_menu{vector_y_index};
702vector_Callback(handles)
703update_field(handles,vector_y)
[204]704
705%------------------------------------------------------------------------
706% --- Executes on selection change in vector_z.
707function vector_z_Callback(hObject, eventdata, handles)
708%------------------------------------------------------------------------
[648]709vector_z_menu=get(handles.vector_z,'String');
710vector_z_index=get(handles.vector_z,'Value');
711vector_z=vector_z_menu{vector_z_index};
712vector_Callback(handles)
713update_field(handles,vector_z)
[204]714
715%------------------------------------------------------------------------
716% --- Executes on selection change in vec_color.
717function vec_color_Callback(hObject, eventdata, handles)
[227]718%------------------------------------------------------------------------
[204]719index=get(handles.vec_color,'Value');
720string=get(handles.vec_color,'String');
721VarName=string{index};
[648]722vector_Callback(handles)
723update_field(handles,VarName)
[204]724
[227]725%------------------------------------------------------------------------
[648]726% --- Executes on selection change in vector_x or vector_y
727function vector_Callback( handles)
[227]728%------------------------------------------------------------------------
[648]729Field=get(handles.get_field,'UserData');
730vector_x_index=get(handles.vector_x,'Value');
731vector_y_index=get(handles.vector_y,'Value');
732vec_color_index=get(handles.vec_color,'Value');
[204]733
[648]734%% set list of possible coordinates
735test_component=zeros(size(Field.Display.VarDimName));%=1 when variable #ilist is eligible as unstructured coordinate
736test_coord=zeros(size(Field.Display.VarDimName)); %=1 when variable #ilist is eligible as structured coordiante
737check_consistent=1;%check that the selected vector components (and possibly color var) have the same dimensiosn
738ListCoord={''};
739dim_var=Field.Display.VarDimName{vector_x_index};%list of dimensions of the selected variable
740if ~isequal(dim_var,Field.Display.VarDimName{vector_y_index})
741    check_consistent=0;
742elseif vec_color_index~=1 && ~isequal(dim_var,Field.Display.VarDimName{vec_color_index})
743    check_consistent=0;
[204]744end
[672]745% the two vector components have consistent dimensions
[648]746if check_consistent
747    for ilist=1:numel(Field.Display.VarDimName)
748        dimnames=Field.Display.VarDimName{ilist}; %list of dimensions for variable #ilist
749        if isequal(dimnames,dim_var)
750            test_component(ilist)=1;
751        elseif numel(dimnames)==1 && ~isempty(find(strcmp(dimnames{1},dim_var)))%variable ilist is a 1D array which can be coordinate variable
752            test_coord(ilist)=1;
[204]753        end
754    end
[675]755    var_component=find(test_component);% list of variable indices elligible as unstructured coordinates
[648]756    var_coord=find(test_coord);% % list of variable indices elligible as structured coordinates
757    ListCoord=Field.Display.ListVarName([var_component var_coord]);
758   
759    %% set default coord selection
760    if numel(find(test_coord))>3
761        set(handles.SwitchVarIndexTime,'Value',3)% the last dim must be considered as time
[204]762    end
[648]763    if numel(var_component)<2
764        if numel(test_coord)<2
765            ListCoord={''};
766        else
767            set(handles.Coord_x,'Value',2)
768            set(handles.Coord_y,'Value',1)
[204]769        end
770    else
[672]771        coord_val=[0 0];
[648]772        for ilist=1:numel(var_component)
773            ivar=var_component(ilist);
774            if isfield(Field.Display,'VarAttribute') && numel(Field.Display.VarAttribute)>=ivar && isfield(Field.Display.VarAttribute{ivar},'Role')
775                Role=Field.Display.VarAttribute{ivar}.Role;
776                if strcmp(Role,'coord_x')
777                    coord_val(1)=ilist;
778                elseif strcmp(Role,'coord_y')
779                    coord_val(2)=ilist;
780                end
[526]781            end
782        end
[675]783        if isempty(find(coord_val))
[672]784            coord_val=var_coord;% case of dimension coordinates
785        end
786        if numel(find(coord_val))<2
[675]787            %coord_val=[numel(var_component)+2 numel(var_component)+1];
788            coord_val=[1 2];
[672]789        end
[648]790        set(handles.Coord_x,'Value',coord_val(1))
791        set(handles.Coord_y,'Value',coord_val(2))
[526]792    end
[204]793end
[648]794set(handles.Coord_y,'String',ListCoord)
795set(handles.Coord_x,'String',ListCoord)
[204]796
797
[648]798%% set list of time coordinates
799menu=get(handles.SwitchVarIndexTime,'String');
800TimeOption=menu{get(handles.SwitchVarIndexTime,'Value')};
801switch TimeOption
802    case 'variable'
803        if numel(find(test_coord))<3
804            ListTime={''};
805        else
806            ListTime=Field.Display.ListVarName(find(test_coord,end));
[204]807        end
[648]808        set(handles.TimeName,'Value',1)
809        set(handles.TimeName,'String',ListTime)
810    case 'dim index'
811        if numel(find(test_coord))<3
812            ListTime={''};
813        else
814            ListTime=Field.Display.VarDimName{find(test_coord,end)};
[204]815        end
[648]816        set(handles.TimeName,'Value',1)
817        set(handles.TimeName,'String',ListTime)
818end 
[204]819
[227]820%------------------------------------------------------------------------
[648]821% --- Executes on selection change in SwitchVarIndexX.
[227]822%------------------------------------------------------------------------
[648]823function SwitchVarIndexX_Callback(hObject, eventdata, handles)
[227]824
825%------------------------------------------------------------------------
[648]826% --- Executes on selection change in Coord_x.
[672]827%------------------------------------------------------------------------
[648]828function Coord_x_Callback(hObject, eventdata, handles)
[672]829
[648]830index=get(handles.Coord_x,'Value');
831string=get(handles.Coord_x,'String');
832VarName=string{index};
833update_field(handles,VarName)
[227]834
835%------------------------------------------------------------------------
[648]836% --- Executes on selection change in Coord_y.
[672]837%------------------------------------------------------------------------
[648]838function Coord_y_Callback(hObject, eventdata, handles)
[672]839
[648]840index=get(handles.Coord_y,'Value');
841string=get(handles.Coord_y,'String');
842VarName=string{index};
843update_field(handles,VarName)
[227]844
[646]845%------------------------------------------------------------------------
[648]846% --- Executes on selection change in Coord_z.
[672]847%------------------------------------------------------------------------
[648]848function Coord_z_Callback(hObject, eventdata, handles)
[672]849
[648]850index=get(handles.Coord_z,'Value');
851string=get(handles.Coord_z,'String');
852VarName=string{index};
853update_field(handles,VarName)
[581]854
[646]855%------------------------------------------------------------------------
[581]856% --- Executes on selection change in SwitchVarIndexTime.
[646]857%------------------------------------------------------------------------
[581]858function SwitchVarIndexTime_Callback(hObject, eventdata, handles)
[672]859
[648]860Field=get(handles.get_field,'UserData');
[581]861menu=get(handles.SwitchVarIndexTime,'String');
862option=menu{get(handles.SwitchVarIndexTime,'Value')};
[648]863
[581]864switch option
865    case 'file index'
[648]866        set(handles.TimeName, 'Visible','off')% the time is taken as the file index
[581]867    case 'attribute'
[648]868        set(handles.TimeName, 'Visible','on')% timeName menu represents the available attributes
[674]869        time_index=get(handles.TimeName,'UserData');    %select the input data
870        if isempty(time_index)
871            PreviousList=get(handles.TimeName, 'String');
872            if ~isempty(PreviousList)
873                PreviousAttr=PreviousList{get(handles.TimeName, 'Value')};
874                index=find(strcmp(PreviousAttr,Field.Display.ListGlobalAttributes),1);
875            end
[581]876        end
[674]877        if isempty(time_index)
878            time_index=find(~cellfun('isempty',regexp(Field.Display.ListGlobalAttribute,'Time')),1);% index of the attributes containing the string 'Time'
879        end     
[581]880        if ~isempty(time_index)
[674]881            set(handles.TimeName,'Value',time_index)
[581]882        else
[648]883            set(handles.TimeName,'Value',1)
[581]884        end
[648]885        set(handles.TimeName, 'String',Field.Display.ListGlobalAttribute)
[674]886
[648]887    case 'variable'% TimeName menu represents the available variables
888        set(handles.TimeName, 'Visible','on')
889        TimeVarName=Field.Display.SingleVarName;% slist of variables with a single dimension (candidate for time)
890        List=get(handles.TimeName,'String');
891        option=List{get(handles.TimeName,'Value')};
[582]892        ind=find(strcmp(option,TimeVarName));
[581]893        if isempty(ind)
[648]894            set(handles.TimeName, 'Value',1);
[581]895        else
[648]896            set(handles.TimeName, 'Value',ind);
[581]897        end
[648]898        set(handles.TimeName, 'String',TimeVarName)
[693]899    case 'matrix_index'% TimeName menu represents the available dimensions
900        set(handles.TimeName, 'Visible','on')     
901        set(handles.TimeName, 'Value',1);
902        set(handles.TimeName, 'String',Field.Display.ListDimName)
[581]903end
904
[648]905%-----------------------------------------------------------------------
906function update_field(handles,VarName)
907%-----------------------------------------------------------------------
908Field=get(handles.get_field,'UserData');
909index=name2index(VarName,Field.ListVarName);
910if ~isempty(index)
911    set(handles.variables,'Value',index+1)
912    variables_Callback(handles.variables, [], handles)
[644]913end
[648]914
[672]915%------------------------------------------------------------------------
916% --- give index numbers of the strings str in the list ListvarName
917% -----------------------------------------------------------------------
[648]918function VarIndex_y=name2index(cell_str,ListVarName)
[672]919
[648]920VarIndex_y=[];
921if ischar(cell_str)
922    for ivar=1:length(ListVarName)
923        varlist=ListVarName{ivar};
924        if isequal(varlist,cell_str)
925            VarIndex_y= ivar;
926            break
927        end
928    end
929elseif iscell(cell_str)
930    for isel=1:length(cell_str)
931        varsel=cell_str{isel};
932        for ivar=1:length(ListVarName)
933            varlist=ListVarName{ivar};
934            if isequal(varlist,varsel)
935                VarIndex_y=[VarIndex_y ivar];
936            end
937        end
938    end
939end
[693]940
941% --- Executes on button press in CheckDimensionY.
942function CheckDimensionX_Callback(hObject, eventdata, handles)
943FieldList=get(handles.FieldOption,'String');
944FieldOption=FieldList{get(handles.FieldOption,'Value')};
945switch FieldOption
946    case '1D plot'
947       
948    case 'scalar'
949       scalar_Callback(hObject, eventdata, handles)
950    case 'vectors'
951end
952
953% --- Executes on button press in CheckDimensionY.
954function CheckDimensionY_Callback(hObject, eventdata, handles)
955FieldList=get(handles.FieldOption,'String');
956FieldOption=FieldList{get(handles.FieldOption,'Value')};
957switch FieldOption
958    case '1D plot'
959       
960    case 'scalar'
961       scalar_Callback(hObject, eventdata, handles)
962    case 'vectors'
963end
964% hObject    handle to CheckDimensionY (see GCBO)
965% eventdata  reserved - to be defined in a future version of MATLAB
966% handles    structure with handles and user data (see GUIDATA)
967
968% Hint: get(hObject,'Value') returns toggle state of CheckDimensionY
Note: See TracBrowser for help on using the repository browser.