tput titles below image, change step 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 07d35441b424f7a559c66603cd19d45ab4616c04
(DIR) parent f4b79e447518f809e2106f464ecc90211f8640cd
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Tue, 10 Nov 2015 11:39:18 +0100
put titles below image, change step plots
Diffstat:
M index.php | 17 +++++++++++++----
M matlab/generate_plots.m | 38 ++++++++++++++++++++++++-------
2 files changed, 43 insertions(+), 12 deletions(-)
---
(DIR) diff --git a/index.php b/index.php
t@@ -77,10 +77,12 @@ if (isset($_GET['wait_id']) && !empty($_GET['wait_id'])) {
<div class="card-image">
<img src="output/<?php
echo($_GET['results_id']); ?>_Walks-6.png">
- <span class="card-title blue-text text-darken-2">
- Generalized model parameter values</span>
+ <!--<span class="card-title blue-text text-darken-2">
+ Generalized model parameter values</span>-->
</div>
<div class="card-content">
+ <h4 class="blue-text">
+ Generalized model parameter values</h4>
<p>The histograms show the distribution of (a)
interglacial erosion rate, (b) glacial erosion rate, (c)
timing of last deglaciation, and (d)
t@@ -117,10 +119,14 @@ if (isset($_GET['wait_id']) && !empty($_GET['wait_id'])) {
<div class="card-image">
<img src="output/<?php
echo($_GET['results_id']); ?>_Walks-5.png">
+ <!--
<span class="card-title blue-text text-darken-2">
- Model parameter values per inversion walker</span>
+ Model parameter values per inversion walker</span>
+ -->
</div>
<div class="card-content">
+ <h4 class="blue-text">
+ Model parameter values per inversion walker</h4>
<p>The histograms show the distribution of (a)
interglacial erosion rate, (b) glacial erosion rate, (c)
timing of last deglaciation, and (d)
t@@ -153,10 +159,13 @@ if (isset($_GET['wait_id']) && !empty($_GET['wait_id'])) {
<div class="card-image">
<img src="output/<?php
echo($_GET['results_id']); ?>_Walks-7.png">
+ <!--
<span class="card-title blue-text text-darken-2">
- Erosion history</span>
+ Erosion history</span>-->
</div>
<div class="card-content">
+ <h4 class="blue-text">
+ Erosion history</h4>
<p>
(a) Two-stage glacial-interglacial model based on the
selected climate record, based on the median
(DIR) diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
t@@ -359,13 +359,16 @@ colneg = [1,0,0];
midvalue = record_threshold_med;
axh(1)=subplot(3,1,1);
[~,~,~,i_cross]=fill_red_blue(xs,ys,colpos,colneg,midvalue);
+text(0.02,0.98,'a', 'Units', 'Normalized', 'VerticalAlignment', 'Top')
xlabel('Age BP [Ma]')
ylabel('\delta^{18}O')
axis tight
-axis([-0.1,2.7,3.0,5.2])
+xlim([min(xs), max(xs)])
+%axis([-0.1,2.7,3.0,5.2])
%axis([-0.05,0.3,3.0,5.2])
-xlim([-0.1,2.7])
+%xlim([-0.1,2.7])
%xlim([-0.05,0.3])
+
hold on
%plot(A(:,1)/1000,A(:,2),'.-m')
axis ij
t@@ -376,29 +379,48 @@ xs_cross = [0;xs_cross];
xs_cross(2) = 0.011;
axh(2)=subplot(3,1,2);
-stairs(xs_cross,(1+-1*(-1).^(1:length(xs_cross)))/2,'b','linewidth',1.5);
+%stairs(xs_cross,(1+-1*(-1).^(1:length(xs_cross)))/2,'b','linewidth',1.5);
+stairs(xs_cross, ...
+ (1+-1*(-1).^(1:length(xs_cross)))/2 * 100, ...
+ 'b','linewidth',1.5);
hold on
start1 = [xs_cross(end),2.7];
start2 = [1,1];
+%plot(start1,start2,'b','linewidth',1.5);
plot(start1,start2,'b','linewidth',1.5);
%title('Exposure. 0 = glaciated, 1 = not glaciated')
-axis([-0.1,2.7,-0.5,1.5])
+%axis([-0.1,2.7,-0.5,1.5])
%axis([-0.05,0.3,-0.5,1.5])
+text(0.02,0.98,'b', 'Units', 'Normalized', 'VerticalAlignment', 'Top')
xlabel('Age BP [Ma]')
+ylabel('Exposure [%]')
+axis tight
+xlim([min(xs), max(xs)])
+ylim([-20, 120])
axh(2)=subplot(3,1,3);
-stairs(xs_cross,(1+-1*(-1).^(1:length(xs_cross)))/2,'r','linewidth',1.5);
+%stairs(xs_cross,(1+-1*(-1).^(1:length(xs_cross)))/2 ,'r','linewidth',1.5);
+stairs(xs_cross, ...
+ epsilon_gla_med + ...
+ -epsilon_gla_med*(1+-1*(-1).^(1:length(xs_cross)))/2 + ...
+ epsilon_int_med*(1+-1*(-1).^(1:length(xs_cross)))/2, ...
+ 'r','linewidth',1.5);
hold on
plot(start1,start2,'r','linewidth',1.5);
%title('Exposure. 0 = glaciated, 1 = not glaciated')
-axis([-0.1,2.7,-1,2])
+%axis([-0.1,2.7,-1,2])
%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]')
+axis tight
+xlim([min(xs), max(xs)])
hold on
d18Oth = midvalue;
-ErateInt=1e-6; ErateGla=1e-7;
+%ErateInt=1e-6; ErateGla=1e-7;
+ErateInt = epsilon_int_med;
+ErateGla = epsilon_gla_med;
[~,~] = ExtractHistory2(ti,d18O_triang,d18Oth,ErateInt,ErateGla);
linkaxes(axh,'x')