timprove presentation of the histograms figure 4 - 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 331b9b208c103e35308172c1c6c7019c1fba34f4
(DIR) parent 3496a2f8ba02021438f2077ffffca7740be3746e
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Thu, 22 Oct 2015 14:51:35 +0200
improve presentation of the histograms figure 4
Diffstat:
M index.php | 12 ++++++++++--
M matlab/generate_plots.m | 11 ++++++++++-
2 files changed, 20 insertions(+), 3 deletions(-)
---
(DIR) diff --git a/index.php b/index.php
t@@ -68,7 +68,7 @@ if (isset($_GET['wait_id']) && !empty($_GET['wait_id'])) {
<div class="container">
<br><br>
<!-- page header -->
- <h2 class="header center orange-text">Inversion results</h2>
+ <h2 class="header center orange-text">MCMC Inversion results</h2>
<div class="row">
t@@ -158,7 +158,15 @@ if (isset($_GET['wait_id']) && !empty($_GET['wait_id'])) {
<span class="card-title">Fourth figure</span>
</div>
<div class="card-content">
- <p>This is the fourth figure, and it shows this and that.</p>
+ <p>The histograms show the distribution of (a)
+ interglacial erosion rate, (b) glacial erosion rate, (c)
+ timing of last deglaciation, and (d)
+ d<sup>18</sup>O<sub>threshold</sub> levels that provide
+ the best fit to the supplied TCN concentrations. The
+ <!--fraction indicates the number of simulations included in-->
+ vertical axis indicates the number of simulations included
+ in each bin out of the 10,000 simulations that followed
+ the MCMC burn-in phase.</p>
</div>
<div class="card-action">
<a href="output/<?php
(DIR) diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
t@@ -189,12 +189,21 @@ for i1 = 1:M
xlabel(fixed_stuff.mname{i1})
if i1 == 1
xlabel('Interglacial erosion rate [mm/yr]')
+ text(0.02,0.98,'a', Units', ...
+ 'Normalized', 'VerticalAlignment', 'Top')
elseif i1 == 2
xlabel('Glacial erosion rate [mm/yr]')
+ text(0.02,0.98,'b', 'Units', ...
+ 'Normalized', 'VerticalAlignment', 'Top')
elseif i1 == 3
xlabel('Timing of last deglaciation [yr]')
+ text(0.02,0.98,'c', 'Units', ...
+ 'Normalized', 'VerticalAlignment', 'Top')
elseif i1 == 4
- xlabel('$\delta^{18} O_\text{threshold}$','Interpreter','LaTeX')
+ xlabel('$\delta^{18} O_\mathrm{threshold}$', ...
+ 'Interpreter', 'LaTeX')
+ text(0.02,0.98,'d','Units', ...
+ 'Normalized', 'VerticalAlignment', 'Top')
end
%set (gca,'xtick',[1e-7:1e-3]);