function [ FIGURE ] = real_plot_1D( data, options )
nplots = 0;
if data.Nx > 1; nplots = nplots + 1; end;
if data.Ny > 1; nplots = nplots + 1; end;
if data.Nz > 1; nplots = nplots + 1; end;


options.PLAN      = 'xy';
options.COMP      = options.COMPZ;
options.POLARPLOT = 0;
options.INTERP    = 0;

toplot = process_field(data,options);
t = data.Ts3D; frames = toplot.FRAMES;

colors = jet(numel(frames));

switch options.NAME
    case '\Gamma_x'
    FIGURE.fig = figure; FIGURE.FIGNAME = ['transport_1D_avg_',data.PARAMS]; 
    yname = '\Gamma';
    case 'v_y'
    FIGURE.fig = figure; FIGURE.FIGNAME = ['ZF_1D_avg_',data.PARAMS]; 
    yname = 'v_y';
    case 's_y'
    FIGURE.fig = figure; FIGURE.FIGNAME = ['ZS_1D_avg_',data.PARAMS]; 
    yname = 's_y';
    case '\phi'
    FIGURE.fig = figure; FIGURE.FIGNAME = ['phi_1D_avg_',data.PARAMS]; 
    yname = '\phi';
    case 'n_i'
    FIGURE.fig = figure; FIGURE.FIGNAME = ['ni_1D_avg_',data.PARAMS]; 
    yname = 'n_i';
    case 'n_e'
    FIGURE.fig = figure; FIGURE.FIGNAME = ['ne_1D_avg_',data.PARAMS]; 
    yname = 'n_e';
end

switch options.COMPX
    case 'avg'
        compx  = @(x) mean(x,2);
        ynamex = ['$\langle ',yname,'\rangle_x$'];
    otherwise
        compx  = @(x) x(:,1);
        ynamex = ['$',yname,'(y,x=0)$'];
end
    
switch options.COMPY
    case 'avg'
        compy  = @(x) mean(x,1);
        ynamey = ['$\langle ',yname,'\rangle_y$'];
    otherwise
        compy  = @(x) x(1,:);
        ynamey = ['$',yname,'(x,y=0)$'];
end
    

set(gcf, 'Position',  [20 50 1000 500])
subplot(1,nplots,1)

    X     = data.x;
    xname = '$x$';
    switch options.COMPT
        case 'avg'
            Y    = compy(mean(toplot.FIELD(:,:,:),3));
            Y    = squeeze(Y);
            if options.NORM
                Y    = Y./max(abs(Y));
            end    
            plot(X,Y,'DisplayName',['$t\in[',num2str(t(frames(1))),',',num2str(t(frames(end))),']$'],...
            'Color',colors(1,:)); hold on;
            legend('show')
        otherwise
            for it = 1:numel(toplot.FRAMES)
                Y    = compy(toplot.FIELD(:,:,it));
                Y    = squeeze(Y);
                if options.NORM
                    Y    = Y./max(abs(Y));
                end
                plot(X,Y,'DisplayName',['$t=',num2str(t(frames(it))),'$'],...
                    'Color',colors(it,:)); hold on;
            end
            legend('show','Location','eastoutside')
    end
    grid on
    xlabel(xname); ylabel(ynamey)
    xlim([min(X),max(X)]);
    
subplot(1,nplots,2)

    X     = data.y;
    xname = '$y$';
    switch options.COMPT
        case 'avg'
            Y    = compx(mean(toplot.FIELD(:,:,:),3));
            Y    = squeeze(Y);
            if options.NORM
                Y    = Y./max(abs(Y));
            end    
            plot(X,Y,'DisplayName',['$t\in[',num2str(t(frames(1))),',',num2str(t(frames(end))),']$'],...
            'Color',colors(1,:)); hold on;
            legend('show')
        otherwise
            for it = 1:numel(toplot.FRAMES)
                Y    = compx(toplot.FIELD(:,:,it));
                Y    = squeeze(Y);
                if options.NORM
                    Y    = Y./max(abs(Y));
                end
                plot(X,Y,'DisplayName',['$t=',num2str(t(frames(it))),'$'],...
                    'Color',colors(it,:)); hold on;
            end
            legend('show','Location','eastoutside')
    end

    grid on
%     title('HeLaZ $k_y$ transport spectrum'); 
    legend('show','Location','eastoutside');
    xlabel(xname); ylabel(ynamex)
        xlim([min(X),max(X)]);
        
%% Z plot
if data.Nz > 1
    options.PLAN      = 'yz';
    options.COMP      = options.COMPX;
    options.POLARPLOT = 0;
    options.INTERP    = 0;

    toplot = process_field(data,options);
    t = data.Ts3D; frames = toplot.FRAMES;
 
switch options.COMPZ
    case 'avg'
        compz  = @(x) mean(x,1);
        ynamez = ['$\langle ',yname,'\rangle_y$'];
    otherwise
        compz  = @(x) x(1,:);
        ynamez = ['$',yname,'(z,y=0)$'];
end

subplot(1,nplots,3)

    X     = data.z;
    xname = '$z$';
    switch options.COMPT
        case 'avg'
            Y    = compz(mean(toplot.FIELD(:,:,:),3));
            Y    = squeeze(Y);
            if options.NORM
                Y    = Y./max(abs(Y));
            end    
            plot(X,Y,'DisplayName',['$t\in[',num2str(t(frames(1))),',',num2str(t(frames(end))),']$'],...
            'Color',colors(1,:)); hold on;
            legend('show')
        otherwise
            for it = 1:numel(toplot.FRAMES)
                Y    = compx(toplot.FIELD(:,:,it));
                Y    = squeeze(Y);
                if options.NORM
                    Y    = Y./max(abs(Y));
                end
                plot(X,Y,'DisplayName',['$t=',num2str(t(frames(it))),'$'],...
                    'Color',colors(it,:)); hold on;
            end
            legend('show','Location','eastoutside')
    end

    grid on
%     title('HeLaZ $k_y$ transport spectrum'); 
    legend('show','Location','eastoutside');
    xlabel(xname); ylabel(ynamez)
        xlim([min(X),max(X)]);

end
end
