c_pl_tx.m

function out=c_pl_tx(varargin)
%C_PL_TX   Plot data from all four Cluster spacecraft in the same plot
%
% C_PL_TX(x1,x2,x3,x4,[column],[linestyle],[dt1 dt2 dt3 dt4])
% C_PL_TX('x?',[column],[linestyle],[dt1 dt2 dt3 dt4])
%	plot variables x1,x2,x3,x4 with time shift dt1...dt4
%	time is 1st column, default plot 2nd column
% C_PL_TX(...,'sc_list',sc_list) specify list of spacecraft to plot
% H=C_PL_TX(..) return handle to plot
% C_PL_TX(AX,...) plot in the specified axis
%
%   column - gives which column to plot. All columns will be plotted
%            in separate panels if set to empty string or ommited.
%   linestyle - string or cell (size 4) in format accepted by plot.
%            Usefull to set line style and marker (but not color).
%   dt1 dt2 dt3 dt4 - timeshifts array
%
%   Example:
%      c_pl_tx('irf_abs(B?)')
%      % plot 3 components + magnitude of B1:4.
%      c_pl_tx('B?',3:4)
%      % plot 3th and 4th components of B1:4.
%      c_pl_tx('B?',3:4,[0 2 3 .5])
%      % plot 3th and 4th components of B1:4 with timeshifts
%      c_pl_tx('B?','',[0 2 3 .5])
%      % plot all components of B1:4 with timeshifts
%      c_pl_tx('B?','.-')
%      % plot all components of B1:4 using '.-' (line with dot markers)
%      c_pl_tx('B?','',[0 2 3 .5],{'.-','*','+','-'})
%      % plot all components of B1:4 using timeshifts and individual
%      % linestyles for each sc
%
% See also IRF_PLOT, PLOT
%
% $Id$

[ax,args,nargs] = axescheck(varargin{:});
if isempty(ax), % if empty axes
    ax=gca;
end
hcf=get(ax,'parent'); % get figure handle

if nargs == 0, % show help if no input parameters
    help c_pl_tx;
    return
end

% Defaults
sc_list=1:4;	% default plot all s/c data
column = [];	% column to plot 
delta_t = [];	% time shift
line_style = {};% line styles

% Check which are input variables
if ischar(args{1})
    % Variables defined in form 'B?'
    getVariablesFromCaller = true;
    variableNameInCaller   = args{1};
    args = args(2:end);
else
    % Variables given as 4 input paramters
    if length(args)<4, error('use c_pl_tx(x1,x2,x3,x4) or c_pl_tx(''x?'')'), end
    c_eval('x? = args{?};'); % assign x1,x2..x4
    args = args(5:end);
    getVariablesFromCaller = false;
end

% Check if column to plot is specified as input
if ~isempty(args)
    if isnumeric(args{1})
        column = args{1};
        args(1) = [];
    elseif isempty(args{1}) % empty string means default matrix size
        args(1) = [];
    end
end

% check for other input parameters
while ~isempty(args)
    if ischar(args{1})
        if strcmp(args{1},'sc_list')
            args(1) = [];
            sc_list=args{1};
			if isempty(sc_list),
				irf_log('fcal','sc_list empty');
				return;
			end
        else
            % assume that argument defines Linestyle
            if isempty(line_style), c_eval('line_style(?)={args{1}};')
            else irf_log('fcal','L_STYLE is already set')
            end
        end
    elseif iscell(args{1}) && length(args{1})==4
        % Individual linestyles for each sc
        if isempty(line_style), line_style = args{1};
        else irf_log('fcal','L_STYLE is already set')
        end
    elseif iscell(args{1})
        % Individual linestyles for each sc
        irf_log('fcal','L_STYLE must be a cell with 4 elements')
    elseif isnumeric(args{1}) && length(args{1})==4
        % dt1..dt4
        if isempty(delta_t), delta_t = args{1};
        else irf_log('fcal','DELTA_T is already set')
        end
    else
        irf_log('fcal',['ignoring input argument: ' args{1}])
    end
	args(1) = [];
end

% Get variable values from caller if needed
if getVariablesFromCaller,
    for cl_id=sc_list,
        ttt = evalin('caller',irf_ssub(variableNameInCaller,cl_id),'[]'); 
        c_eval('x? =ttt;',cl_id); clear ttt
    end
end

% If column empty, check which columns to plot
if isempty(column)
	for ic=sc_list
		eval(['nCol = size(x' num2str(ic) ',2);']);
		if ~isempty(nCol) && nCol > 0
			if nCol >= 2,
				column = 2:nCol;
			end
			break;
		end
	end
end
if isempty(column)
    irf_log('fcal','all inputs are empty')
    return
end

% define Cluster colors
cluster_colors={'[0 0 0]';'[1 0 0]';'[0 0.5 0]';'[0 0 1]'};
l_style=cell(1,4);
if isempty(line_style),
    for ic=1:4, l_style(ic)={['''color'','  cluster_colors{ic}]};end
else
    for ic=1:4, l_style(ic)={['''' line_style{ic} ''',''color'','  cluster_colors{ic}]};end
end

% t_start_epoch is saved in figures 'userdata' variable
% check first if it exist otherwise set it if variables
% use isdat epoch or put to zero without changing 'userdata'
ud=get(hcf,'userdata');
t_start_epoch = double(0); % default
if isfield(ud,'t_start_epoch'),
	t_start_epoch = double(ud.t_start_epoch);
else
	variablesUseEpochTime = false;
	c_eval(['variablesUseEpochTime = or(variablesUseEpochTime,' ...
		'(~isempty(x?) && x?(1,1)>1e8) );'],sc_list);
	if variablesUseEpochTime 
		% Set start_epoch if time is in isdat epoch,
		tt = [];
		c_eval('if ~isempty(x?), tt=[tt; x?(1,1)]; end')
		t_start_epoch = double(min(tt)); clear tt
		ud.t_start_epoch = t_start_epoch; set(hcf,'userdata',ud);
		irf_log('proc',['user_data.t_start_epoch is set to ' ...
			irf_time(t_start_epoch,'iso')]);
	end
end
if isempty(delta_t), delta_t = [0 0 0 0]; end
c_eval('ts?=t_start_epoch+delta_t(?);')

% check which spacecraft data are available
sc_list_with_data=[];
c_eval('if ~isempty(x?), sc_list_with_data=[sc_list_with_data ?];end',sc_list);

% if more than one column reset figure 
if length(column) > 1 && numel(ax) ~= numel(column)
	ax = irf_plot(length(column),'reset'); 
end

for j=1:length(column)
    for jj=sc_list_with_data
        c_eval(['if ~isempty(x?),'...
            'hl=irf_plot(ax(j), [x?(:,1)-delta_t(?) x?(:,column(j))],' l_style{jj} ');'...
            'hold(ax(j),''on'');set(hl,''Tag'',''C?''); end, '],jj);
    end
    hold(ax(j),'off');
    irf_zoom(ax(j),'y'); % optimize Y zoom to skip labels at top and bottom
    grid(ax(j),'on');
end
%ud = get(hcf,'userdata'); ud.subplot_handles = ax; set(hcf,'userdata',ud);
irf_timeaxis(ax);
irf_figmenu;


if nargout > 0, out = ax; end