comsol_plot_sequence.m

function p = comsol_plot_sequence(varargin)
%function p = comsol_plot_sequence(varargin)
%
% This function was mainly copied from plot_sequence and adopted to special
% COMSOL usage. So please see the explanation there...
%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%function p = plot_sequence(sim_data,params[,callbackfigure, cbhandle])
% plotting a sequence of data slices on polygonal
% 2d grid (constructed from params if empty) and providing a slider
% for changing between the data
% slices.
%
% A new figure is opened and the handle returned in p.  If further parameters
% are set, the call is assumed to stem from a callback-function
%
% Every column of data is interpreted as one discrete function
% dof vector forwarded to the params.plot() function.
%
% parameters:
%  varargin: usually called with 3 arguments:
%            @code plot_sequence(data, grid, params) @endcode
%   - data - data vector to be plotted
%   - grid - the underlying grid
%   - params - plotting parameters
%   .
%            Alternatively, there can be 5 arguments:
%            @code plot_sequence(data, grid, params, callbackfigure, cbhandle) @endcode
%            This function is set as a callback function for the time slider
%            uicontrol. The additional arguments are:
%   - callbackfigure - handle to the figure
%   - cbhandle       - handle to the object calling the callback function. This
%                      is usually the time slider.
%
% return values:
%  p : figure handle of plot
%
% required fields in params:
%  - 'plot' -- pointer to the plot-function performing the plotting of a
%              single slice, e.g.  plot_element_data(), plot_vertex_data(),
%              fv_plot(), ldg_plot().
%
% optional field of params:
%  - 'title' -- string indicating the title of the newly opened figure
%  - 'clim'  -- 2-vector giving the range of the colormap. If this is set,
%               identical range is used for all slices. Default is the min and
%               max of all slices.
%  - 'clim_tight' -- if this flag is set, the colorbar is set tightly to the
%                    range of every single data slice. Clearly only one of
%                    'clim' or 'clim_tight' should be set.
%
% see also the chosen 'params.plot_function' for its further params-options
%
% Bernard Haasdonk 9.5.2007
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Oliver Zeeb, 2013



% open figure
if nargin < 3 % no callback, open new figure

  %  keyboard;
  f = figure;
  di = 0.1;
  textwidth= 0.15;
  textheight = 0.05;
  txt = uicontrol(f,'Style','text','Units','normalized',...
                  'Position',[di,di,textwidth,textheight],...
                  'String','Slice 1','Tag','text1');
  slider = uicontrol(f,'Style','slider','Units','normalized',...
                     'Position',[di+textwidth+di, di ,1-3*di-textwidth,...
                    textheight],'Tag','slider1',...
                    'Callback','comsol_plot_sequence([],[],gcbf,gcbo)');  
  ax = axes('Units','normalized',...
            'Position',[di,2*di+textheight,1-2*di,1-3*di - textheight],...
            'Tag','axes1');

  ud = [];
  ud.sim_data = varargin{1};
  ud.params = varargin{2};

  set(f,'Userdata',ud);

  sl = findobj(f,'Tag','slider1');
  %set(sl,'min',0);
  %set(sl,'max',ud.params.nt);
  %set(sl,'sliderstep',[1/ud.params.nt,1]);
  set(sl,'min',1);
  if size(ud.sim_data.U,2)>1
    set(sl,'max',size(ud.sim_data.U,2));
    set(sl,'sliderstep',[1/(size(ud.sim_data.U,2)-1),1]);
  else
    set(sl,'max',size(ud.sim_data.U,2)+eps);
    set(sl,'visible','off');
    %  set(sl,'sliderstep',[0,0]);
  end;
  set(sl,'value',1);
  th = findobj(f,'Tag','text1');
  set(th,'String',('data slice 1') );
  replot(f,sl,1);
  %  set(f,'Resize','on');  
  p = f;
end;

if nargin >=3 % callback-function
  cbf = varargin{3};
  cbo = varargin{4};
  th = findobj(gcbf,'Tag','text1');
  v = round(get(gcbo,'value'));
  set(gcbo,'value',v)
  set(th,'String',['data slice ',num2str(v)]);
  replot(gcbf,gcbo,v);
  p = gcbf;
end;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function replot(fh,oh,v)
ud = get(fh,'Userdata');


ud.sim_data.comsol_model.result('pg1').set('solnum',v);

ax = findobj(fh,'Tag','axes1');
cla(ax);

mphplot(ud.sim_data.comsol_model,'pg1','rangenum',1)