tadd deglaciation timing to plots - 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 dd05f749f8a4c49781ac9dbdacaf68fe2fa15969
 (DIR) parent 22d416e2361d75c329c09cc9afd280b47f341b32
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Tue, 10 Nov 2015 12:02:15 +0100
       
       add deglaciation timing to plots
       
       Diffstat:
         M matlab/generate_plots.m             |       9 ++++++++-
       
       1 file changed, 8 insertions(+), 1 deletion(-)
       ---
 (DIR) diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
       t@@ -376,6 +376,9 @@ xs_cross = (i_cross-1)*dt;
        xs_cross = [0;xs_cross];
        %title([fn,'.  minrad = ',num2str(minrad),' Ma'],'interp','none')
        
       +% deglaciation timing
       +plot([t_degla_med, t_degla_med], get(gca,'YLim'), 'k--');
       +
        xs_cross(2) = 0.011;
        
        % Exposure
       t@@ -398,7 +401,8 @@ ylabel('Exposure [%]')
        axis tight
        xlim([min(xs), max(xs)])
        ylim([-20, 120])
       -
       +% deglaciation timing
       +plot([t_degla_med, t_degla_med], get(gca,'YLim'), 'k--');
        
        % Erosion rate
        axh(2)=subplot(3,1,3);
       t@@ -422,6 +426,9 @@ xlim([min(xs), max(xs)])
        ylims = get(gca,'YLim');
        ylim([ylims(1)*0.95, ylims(2)*1.05])
        
       +% deglaciation timing
       +plot([t_degla_med, t_degla_med], get(gca,'YLim'), 'k--');
       +
        hold on
        d18Oth = midvalue;
        %ErateInt=1e-6; ErateGla=1e-7;