tadd citation info, remove legacy plotting routines - cosmo - front and backend for Markov-Chain Monte Carlo inversion of cosmogenic nuclide concentrations
HTML git clone git://src.adamsgaard.dk/cosmo
DIR Log
DIR Files
DIR Refs
DIR README
DIR LICENSE
---
DIR commit ea8c83fcdb73d972e5ff2b880876c61372cbecb5
DIR parent a93793bedaa193aae5c7d675b9512d372e07377c
HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Thu, 19 Nov 2015 16:31:28 +0100
add citation info, remove legacy plotting routines
Diffstat:
M matlab/generate_plots.m | 165 +------------------------------
M pages/methods.html | 34 +++++++++++++++++++++++++++++++
2 files changed, 36 insertions(+), 163 deletions(-)
---
DIR diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
t@@ -509,169 +509,8 @@ linkaxes(axh,'x')
% stairs(-tStarts,relExpos+2,'r')
-% %% Exhumation history from InspectDepthConcTracks_medTrueDepthIndsat.m
-% disp('Exhumation history')
-% %[fn,pn]=uigetfile('*.mat')
-% %load([pn,fn])
-% DoConfLevelTrim = 1;
-%
-% fh = [fh;figure('visible', show_figures)];
-%
-% Nwalks = length(Ss);
-% for iwalk = 1:Nwalks
-% % iwalk=input(['What iwalk?[1..',num2str(length(Ss)),']']),
-% subplot(Nwalks,1,iwalk)
-% lump_MetHas = Ss{iwalk}.lump_MetHas;
-% fixed_stuff = Ss{iwalk}.fs;
-% % if ~isempty(lump_MetHas.zsss{1})
-% % % % % zsss = lump_MetHas.zsss;
-% % % % % tss = lump_MetHas.tss;
-% % % % % dtfine = -1000; tfinemax = -20e5; dzfine = 0.1; zfinemax = 20,
-% % % % % iz = 1;
-% % % % % DepthExposureTimeDens(zsss,tss,dtfine,tfinemax,dzfine,zfinemax,fixed_stuff,iz)
-% % % % % title('z(time)')
-% % % % % axis ij
-% % % % % hold on
-%
-% %making quantiles to plot ontop of exhumation densities:
-% tss = lump_MetHas.tss;
-% dtfine = -1000; tfinemax = -20e5;
-% Xxsss = lump_MetHas.zsss;
-% dXxfine = 0.01; Xxfinemax = 10; iz=1; %Xx may be depth or nucleide concentration or ...
-% dzfine = dXxfine; zfinemax = Xxfinemax;
-% [smoothgrid,histgrid,tsfine,Xxfine]=XxTimeDens(Xxsss,tss,dtfine,tfinemax,dXxfine,Xxfinemax,fixed_stuff,iz);
-% pcolor(tsfine,Xxfine,sqrt(smoothgrid));
-% hold on
-% %plot a selection of depth histories
-% for i=1:ceil(length(tss)/50):length(tss)
-% plot(tss{i},Xxsss{i},'.-w')
-% end
-% %Compute and plot quantiles
-% N = sum(histgrid(:,1));
-% fractions = [0.25,0.5,0.75]; %quartiles
-% tsfine = 0:dtfine:tfinemax; %bin boundaries
-% Ntfine = length(tsfine);
-% zsfine = 0:dzfine:zfinemax; %bin boundaries
-% % quants{iwalk} = GetHistgridQuantiles(histgrid,N,fractions,tsfine,zsfine);
-% % plot(tsfine,quants{iwalk}(1,:),'r','linewidth',3)
-% % plot(tsfine,quants{iwalk}(2,:),'k','linewidth',3)
-% % plot(tsfine,quants{iwalk}(3,:),'b','linewidth',3)
-% quants(:,:,iwalk) = GetHistgridQuantiles(histgrid,N,fractions,tsfine,zsfine);
-%
-% grid on; shading flat; axis tight; set(gca,'fontsize',8); hold on
-% lh(1)=plot(tsfine,quants(1,:,iwalk),'r','linewidth',2);
-% lh(2)=plot(tsfine,quants(2,:,iwalk),'k','linewidth',2);
-% lh(3)=plot(tsfine,quants(3,:,iwalk),'g','linewidth',2);
-% legend(lh,'25%','median','75%','location','northwest')
-% axis ij
-% %track_handle=AddTrueModelDepthTrack(Ss{iwalk}.fs,'*-m'); %<<< BHJ: added 2014 dec 09
-% end
-% colormap(1-copper)
-% % subplot(2,2,1);
-
-% the figure below is really computationally intensive
-% fh = [fh;figure('visible', show_figures)];
-% for iwalk = 1:Nwalks
-% plot(tsfine,quants(1,:,iwalk),'r','linewidth',3); hold on
-% plot(tsfine,quants(2,:,iwalk),'k','linewidth',3)
-% plot(tsfine,quants(3,:,iwalk),'g','linewidth',3)
-% end
-% %axis([-2e6 0 0 4])
-% axis ij
-
-% load handel; sound(y(1:5000))
-% iwalk=input(['What iwalk?[1..',num2str(length(Ss)),']']);
-% lump_MetHas = Ss{iwalk}.lump_MetHas;
-% fixed_stuff = Ss{iwalk}.fs;
-% % if ~isempty(lump_MetHas.zsss{1})
-% fh = [fh;figure('visible', show_figures)];
-% zsss = lump_MetHas.zsss;
-% tss = lump_MetHas.tss;
-% dtfine = -1000; tfinemax = -20e5; dzfine = 0.1; zfinemax = 20,
-% iz = 1;
-% DepthExposureTimeDens(zsss,tss,dtfine,tfinemax,dzfine,zfinemax,fixed_stuff,iz)
-% title('z(time)')
-% axis ij
-% hold on
-%
-% %making quantiles to plot ontop of exhumation densities:
-% Xxsss = lump_MetHas.zsss;
-% dXxfine = 0.1; Xxfinemax = 10;
-% [smoothgrid,histgrid,tsfine,Xxfine]=XxTimeDens(Xxsss,tss,dtfine,tfinemax,dXxfine,Xxfinemax,fixed_stuff,iz);
-% %plot a selection of depth histories
-% for i=1:ceil(length(tss)/50):length(tss)
-% plot(tss{i},zsss{i},'.-w')
-% end
-%
-% %Compute and plot quantiles
-% N = sum(histgrid(:,1));
-% fractions = [0.25,0.5,0.75]; %quartiles
-% tsfine = 0:dtfine:tfinemax; %bin boundaries
-% Ntfine = length(tsfine);
-% zsfine = 0:dzfine:zfinemax; %bin boundaries
-% quants = GetHistgridQuantiles(histgrid,N,fractions,tsfine,zsfine);
-% plot(tsfine,quants(1,:),'r','linewidth',3)
-% plot(tsfine,quants(2,:),'k','linewidth',3)
-% plot(tsfine,quants(3,:),'g','linewidth',3)
-
-% fh = [fh;figure('visible', show_figures)];
-% zsss = lump_MetHas.zsss;
-% tss = lump_MetHas.ExposureTimeSinceNows;
-% dtExposurefine = 1000; tExposurefinemax = 20e5; dzfine = 0.1; zfinemax = 20,
-% iz = 1;
-% DepthExposureTimeDens(zsss,tss,dtExposurefine,tExposurefinemax,dzfine,zfinemax,fixed_stuff,iz)
-% title('z(ExposureTime)')
-% % end
-% axis ij
-% hold on
-% for i=1:ceil(length(tss)/50):length(tss)
-% plot(tss{i},zsss{i},'.-w')
-% end
-%
-% figure %five subplots:
-% tss = lump_MetHas.tss;
-% dtfine = -1000; tfinemax = -20e5;
-% iz = 1;
-%
-% axh(1)=subplot(5,1,1) %depth track
-% Xxsss = lump_MetHas.zsss;
-% dXxfine = 0.1; Xxfinemax = 10;
-% [smoothgrid,histgrid,tsfine,Xxfine]=XxTimeDens(Xxsss,tss,dtfine,tfinemax,dXxfine,Xxfinemax,fixed_stuff,iz);
-% pcolor(tsfine,Xxfine,sqrt(smoothgrid));
-% % axis(sqrt(3)*[-1,1,-1,1])
-% % set(gca,'xtick',[-1.5:0.5:1.5],'ytick',[-1.5:0.5:1.5])
-% grid on; shading flat; axis tight; set(gca,'fontsize',8); hold on
-% lh(1)=plot(tsfine,quants(1,:),'r','linewidth',2)
-% lh(2)=plot(tsfine,quants(2,:),'k','linewidth',2)
-% lh(3)=plot(tsfine,quants(3,:),'g','linewidth',2)
-% legend(lh,'25%','median','75%','location','northwest')
-% axis ij
-%
-% Xxsss{2}=lump_MetHas.c14Csss;
-% Xxsss{3}=lump_MetHas.c10Besss;
-% Xxsss{4}=lump_MetHas.c26Alsss;
-% Xxsss{5}=lump_MetHas.c21Nesss;
-%
-% for ic=2:5
-% axh(ic)=subplot(5,1,ic) %10Be track
-% dXxfine = 0.01*Xxsss{ic}{end}(1,end); Xxfinemax = 1.3*Xxsss{ic}{end}(1,end);
-% [smoothgrid,histgrid,tsfine,Xxfine]=XxTimeDens(Xxsss{ic},tss,dtfine,tfinemax,dXxfine,Xxfinemax,fixed_stuff,iz);
-% pcolor(tsfine,Xxfine,sqrt(smoothgrid));
-% % axis(sqrt(3)*[-1,1,-1,1])
-% % set(gca,'xtick',[-1.5:0.5:1.5],'ytick',[-1.5:0.5:1.5])
-% grid on; shading flat; axis tight; set(gca,'fontsize',8,'tickdir','out');
-% hold on
-% end
-% for ic=1:4
-% set(axh(ic),'xticklabel',[])
-% end
-% subplot(5,1,1); title(['Depth track:',fn],'interpreter','none'); set(gca,'layer','top')
-% subplot(5,1,2); title('^{14}C'); set(gca,'layer','top')
-% subplot(5,1,3); title('^{10}Be'); set(gca,'layer','top')
-% subplot(5,1,4); title('^{26}Al'); set(gca,'layer','top')
-% subplot(5,1,5); title('^{21}Ne'); set(gca,'layer','top')
-% linkaxes(axh,'x')
-% colormap(flipud(gray))
+%% Exhumation history
+
%% position figure windows at certain coordinates on the screen
DIR diff --git a/pages/methods.html b/pages/methods.html
t@@ -24,5 +24,39 @@
flexible.
</p>
+ <p>For a full description of the methods and application possibilities,
+ see the open-access publication by <a
+ href="http://www.sciencedirect.com/science/article/pii/S1871101415300558">Knudsen
+ et al. (2015)</a>.
+ </p>
+
+ <p class="flow-text">If you use the results generated by this tool in a
+ scientific publication, please acknowledge this fact by citing:<br><br>
+ <blockquote>
+ Knudsen, M.F., Egholm, D.L., Jacobsen, B.H., Larsen, N.K., Jansen,
+ J.D., Andersen, J.L., Linge, H.C., 2015. A multi-nuclide approach to
+ constrain landscape evolution and past erosion rates in previously
+ glaciated terrains. Quaternary Geochronology 30, 100-113,
+ doi:10.1016/j.quageo.2015.08.004.
+ </blockquote>
+ You may use the following BibTeX entry:
+ <div class="row">
+ <pre><code class="language-markup col s12">
+ @article{Knudsen2015100,
+ title = "A multi-nuclide approach to constrain landscape evolution and past erosion rates in previously glaciated terrains ",
+ journal = "Quaternary Geochronology ",
+ volume = "30, Part A",
+ number = "",
+ pages = "100 - 113",
+ year = "2015",
+ note = "",
+ issn = "1871-1014",
+ doi = "http://dx.doi.org/10.1016/j.quageo.2015.08.004",
+ }
+ </code></pre>
+ </div>
+
+
+
</div>
</div>