tfix old labels reflecting internal units - 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 58818076edeed2691b677e3854cac1f09abd84b8
 (DIR) parent 56a526072365d711c52d0350ce0e90670cacaa02
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu,  3 Dec 2015 12:41:14 +0100
       
       fix old labels reflecting internal units
       
       Diffstat:
         M matlab/generate_plots.m             |      14 +++++++-------
       
       1 file changed, 7 insertions(+), 7 deletions(-)
       ---
 (DIR) diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
       t@@ -249,7 +249,7 @@ for i1 = 1:M % for each model parameter
            
            if i1 == 1
                title(['MCMC walker ' num2str(iwalk)])
       -        %xlabel('Interglacial erosion rate [mm/yr]')
       +        %xlabel('Interglacial erosion rate [m/yr]')
                %xlabel('Interglacial erosion rate [m/Myr]')
                xlabel('\epsilon_{int} [m/Myr]')
                text(0.02,0.98,'a', 'Units', ...
       t@@ -259,7 +259,7 @@ for i1 = 1:M % for each model parameter
                epsilon_int_75(iwalk) = prctile(Ss{iwalk}.ms(i1,:)*1.0e6, 75);
                
            elseif i1 == 2
       -        %xlabel('Glacial erosion rate [mm/yr]')
       +        %xlabel('Glacial erosion rate [m/yr]')
                %xlabel('Glacial erosion rate [m/Myr]')
                xlabel('\epsilon_{gla} [m/Myr]')
                text(0.02,0.98,'b', 'Units', ...
       t@@ -362,19 +362,19 @@ for i1 = 1:M % for each model parameter
            hold off
            
            if i1 == 1
       -        %xlabel(['Interglacial erosion rate [mm/yr], median = ' ...
       +        %xlabel(['Interglacial erosion rate [m/yr], median = ' ...
                %xlabel(['Interglacial erosion rate [m/Myr], median = ' ...
                xlabel(['\epsilon_{int} [m/Myr], median = ' ...
                    num2str(med, 4) ' m/Myr'])
       -            %num2str(med, 4) ' mm/yr'])
       +            %num2str(med, 4) ' m/yr'])
                text(0.02,0.98,'a', 'Units', ...
                    'Normalized', 'VerticalAlignment', 'Top')
            elseif i1 == 2
       -        %xlabel(['Glacial erosion rate [mm/yr], median = ' ...
       +        %xlabel(['Glacial erosion rate [m/yr], median = ' ...
                %xlabel(['Glacial erosion rate [m/Myr], median = ' ...
                xlabel(['\epsilon_{gla} [m/Myr], median = ' ...
                    num2str(med, 4) ' m/Myr'])
       -            %num2str(med, 4) ' mm/yr'])
       +            %num2str(med, 4) ' m/yr'])
                text(0.02,0.98,'b', 'Units', ...
                    'Normalized', 'VerticalAlignment', 'Top')
            elseif i1 == 3
       t@@ -500,7 +500,7 @@ plot(start1,start2,'r','linewidth',1.5);
        %axis([-0.05,0.3,-1,2])
        text(0.02,0.98,'c', 'Units', 'Normalized', 'VerticalAlignment', 'Top')
        xlabel('Age BP [Ma]')
       -%ylabel('Erosion rate [mm/yr]')
       +%ylabel('Erosion rate [m/yr]')
        ylabel('Erosion rate [m/Myr]')
        axis tight
        xlim([min(xs), max(xs)])