PaleoTracer/func_plot_FLU_FLUPYR_temporal.m at master · ghosh17/PaleoTracer · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
function [] = func_plot_FLU_FLUPYR_temporal(FLU_FLUPYR, IsotopicValue, numsamples, t)

    c = linspace(t(1),t(numsamples),length(t));
    hold on
    grid on;

    yyaxis left;
    figure('Name', 'FLU/FLUPYR temporal');
    fire_freequency_plot = scatter(t, FLU_FLUPYR,[],c);
    %set(gca,'yscale','log')
    xlabel('Time (Ma)');
    ylabel('FLU/(FLU+PYR)');

    set(gca, 'XDir','reverse')

    buffer = (t(1)-t(numsamples))/numsamples;
    line([t(1)+buffer, t(numsamples)-buffer],[0.4,0.4],'Color','black','LineStyle','--');
    line([t(1)+buffer, t(numsamples)-buffer],[0.5,0.5],'Color','black','LineStyle','--');

    %Secondary plot showing Carbon isotopic ratios
    yyaxis right;
    isotope_plot = scatter(t,IsotopicValue,'*','r');
    ylabel('δ13C (‰)');
    ylim([-35, -10]);
    legend([fire_freequency_plot,isotope_plot], {'Fire freequency' , 'δ13C of Bulk Organic Matter'})


    %%Save plot file
    %Location to save Plots
    PlotFolder = strcat(pwd,'/Plots/');
    filename = 'Combustion_Character_Flu_Pyr_Temporal.jpg';
    %Corner case handling
    %if no Plot folder exists
    if ~exist(PlotFolder, 'dir')
        %create a figure folder
        mkdir(PlotFolder);
    end
    %Save all figures to this folder
    saveas(gcf, fullfile(PlotFolder, filename));

    hold off

end