tupdate xlims to new 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 87bff390226ab8f488960b38784affa6c73ca87c
(DIR) parent 5e3be79e8b3cb07aa8ea4fc4ae044f419668b3ea
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Wed, 18 Nov 2015 16:05:49 +0100
update xlims to new units
Diffstat:
M matlab/generate_plots.m | 30 ++++++++++++++++++++++--------
1 file changed, 22 insertions(+), 8 deletions(-)
---
(DIR) diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
t@@ -250,9 +250,16 @@ for i1 = 1:M % for each model parameter
end
%set (gca,'xtick',[1e-7:1e-3]);
- switch mDistr(i1,:)
- case 'uniform', set(gca,'xscale','lin','xlim',mminmax(i1,:))
- case 'logunif', set(gca,'xscale','log','xlim',mminmax(i1,:))
+ if i1 == 1 || i1 == 2 % shift axes for new units
+ switch mDistr(i1,:)
+ case 'uniform', set(gca,'xscale','lin','xlim',mminmax(i1,:)*1000.)
+ case 'logunif', set(gca,'xscale','log','xlim',mminmax(i1,:)*1000.)
+ end
+ else
+ switch mDistr(i1,:)
+ case 'uniform', set(gca,'xscale','lin','xlim',mminmax(i1,:))
+ case 'logunif', set(gca,'xscale','log','xlim',mminmax(i1,:))
+ end
end
end
end
t@@ -297,7 +304,7 @@ for i1 = 1:M % for each model parameter
elseif i1 == 4
record_threshold_med = med;
else
- error('Unknown parametr');
+ error('Unknown parameter');
end
%ylims = get(gca,'YLim');
%text(med, ylims(1) + (ylims(2) - ylims(1))*0.9, ...
t@@ -315,7 +322,7 @@ for i1 = 1:M % for each model parameter
if i1 == 1
%xlabel(['Interglacial erosion rate [mm/yr], median = ' ...
- xlabel(['Interglacial erosion rate [mm/yr], median = ' ...
+ xlabel(['Interglacial erosion rate [m/Myr], median = ' ...
num2str(med, 4) ' m/Myr'])
%num2str(med, 4) ' mm/yr'])
text(0.02,0.98,'a', 'Units', ...
t@@ -346,9 +353,16 @@ for i1 = 1:M % for each model parameter
end
%set (gca,'xtick',[1e-7:1e-3]);
- switch mDistr(i1,:)
- case 'uniform', set(gca,'xscale','lin','xlim',mminmax(i1,:))
- case 'logunif', set(gca,'xscale','log','xlim',mminmax(i1,:))
+ if i1 == 1 || i1 == 2 % shift axes for new units
+ switch mDistr(i1,:)
+ case 'uniform', set(gca,'xscale','lin','xlim',mminmax(i1,:)*1000.)
+ case 'logunif', set(gca,'xscale','log','xlim',mminmax(i1,:)*1000.)
+ end
+ else
+ switch mDistr(i1,:)
+ case 'uniform', set(gca,'xscale','lin','xlim',mminmax(i1,:))
+ case 'logunif', set(gca,'xscale','log','xlim',mminmax(i1,:))
+ end
end
end