textend stair plot to the end of range - cosmo - front and backend for Markov-Chain Monte Carlo inversion of cosmogenic nuclide concentrations
 (HTM) git clone git://src.adamsgaard.dk/cosmo
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) README
 (DIR) LICENSE
       ---
 (DIR) commit 69cb202539fc5cea6b0be5380913234ba851c00f
 (DIR) parent 07d35441b424f7a559c66603cd19d45ab4616c04
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Tue, 10 Nov 2015 11:53:54 +0100
       
       extend stair plot to the end of range
       
       Diffstat:
         M matlab/generate_plots.m             |      22 +++++++++++++---------
       
       1 file changed, 13 insertions(+), 9 deletions(-)
       ---
 (DIR) diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
       t@@ -378,16 +378,17 @@ xs_cross = [0;xs_cross];
        
        xs_cross(2) = 0.011;
        
       +% Exposure
        axh(2)=subplot(3,1,2);
        %stairs(xs_cross,(1+-1*(-1).^(1:length(xs_cross)))/2,'b','linewidth',1.5);
       -stairs(xs_cross, ...
       -    (1+-1*(-1).^(1:length(xs_cross)))/2 * 100, ...
       +exposure = (1+-1*(-1).^(1:length(xs_cross)))/2 * 100;
       +stairs(xs_cross, exposure, ...
            'b','linewidth',1.5);
        hold on
       -start1 = [xs_cross(end),2.7];
       -start2 = [1,1];
       -%plot(start1,start2,'b','linewidth',1.5);
       +start1 = [xs_cross(end), xs(end)];
       +start2 = [exposure(end), exposure(end)];
        plot(start1,start2,'b','linewidth',1.5);
       +%plot(start1,start2,'b','linewidth',1.5);
        %title('Exposure. 0 = glaciated, 1 = not glaciated')
        %axis([-0.1,2.7,-0.5,1.5])
        %axis([-0.05,0.3,-0.5,1.5])
       t@@ -398,14 +399,17 @@ axis tight
        xlim([min(xs), max(xs)])
        ylim([-20, 120])
        
       +
       +% Erosion rate
        axh(2)=subplot(3,1,3);
        %stairs(xs_cross,(1+-1*(-1).^(1:length(xs_cross)))/2 ,'r','linewidth',1.5);
       -stairs(xs_cross, ...
       -    epsilon_gla_med + ...
       +erosionrate = epsilon_gla_med + ...
            -epsilon_gla_med*(1+-1*(-1).^(1:length(xs_cross)))/2 + ...
       -    epsilon_int_med*(1+-1*(-1).^(1:length(xs_cross)))/2, ...
       -    'r','linewidth',1.5);
       +    epsilon_int_med*(1+-1*(-1).^(1:length(xs_cross)))/2;
       +stairs(xs_cross, erosionrate, 'r','linewidth',1.5);
        hold on
       +start1 = [xs_cross(end), xs(end)];
       +start2 = [erosionrate(end), erosionrate(end)];
        plot(start1,start2,'r','linewidth',1.5);
        %title('Exposure. 0 = glaciated, 1 = not glaciated')
        %axis([-0.1,2.7,-1,2])