tadd extra parameters to plotting function - 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 762780b1231bd44f056fc0d47cd4203463d5d3b8
(DIR) parent 49eb5cb6f34627f32882fe2d6c44a1d3d9d0700b
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Tue, 1 Dec 2015 12:28:16 +0100
add extra parameters to plotting function
Diffstat:
M matlab/file_scanner_mcmc_starter.m | 17 +++++++++++++++--
M matlab/generate_plots.m | 25 +++++++++++++++++++------
2 files changed, 34 insertions(+), 8 deletions(-)
---
(DIR) diff --git a/matlab/file_scanner_mcmc_starter.m b/matlab/file_scanner_mcmc_starter.m
t@@ -115,8 +115,21 @@ while 1
% generate plots
%CompareWalks2(Ss, save_file)
- generate_plots(Ss, record, matlab_scripts_folder, ...
- save_file, graphics_formats, show_figures);
+ generate_plots(Ss, matlab_scripts_folder, ...
+ save_file, graphics_formats, show_figures, ...
+ sample_id, name, email, lat, long, ...
+ be_conc, al_conc, c_conc, ne_conc, ...
+ be_uncer, al_uncer, c_uncer, ne_uncer, ...
+ zobs, ...
+ be_prod_spall, al_prod_spall, c_prod_spall, ne_prod_spall, ...
+ be_prod_muons, al_prod_muons, c_prod_muons, ne_prod_muons, ...
+ rock_density, ...
+ epsilon_gla_min, epsilon_gla_max, ...
+ epsilon_int_min, epsilon_int_max, ...
+ t_degla_min, t_degla_max, ...
+ record, ...
+ record_threshold_min, record_threshold_max, ...
+ nwalkers);
% close all figures
close all;
(DIR) diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
t@@ -1,5 +1,18 @@
-function generate_plots(Ss, record, matlab_scripts_folder,...
- save_file, formats, show_figures)
+function generate_plots(Ss, matlab_scripts_folder,...
+ save_file, formats, show_figures, ...
+ sample_id, name, email, lat, long, ...
+ be_conc, al_conc, c_conc, ne_conc, ...
+ be_uncer, al_uncer, c_uncer, ne_uncer, ...
+ zobs, ...
+ be_prod_spall, al_prod_spall, c_prod_spall, ne_prod_spall, ...
+ be_prod_muons, al_prod_muons, c_prod_muons, ne_prod_muons, ...
+ rock_density, ...
+ epsilon_gla_min, epsilon_gla_max, ...
+ epsilon_int_min, epsilon_int_max, ...
+ t_degla_min, t_degla_max, ...
+ record, ...
+ record_threshold_min, record_threshold_max, ...
+ nwalkers)
%% Copied from m_pakke2014maj11/CompareWalks2.m
% Generates and saves all relevant figures
t@@ -540,7 +553,7 @@ for iwalk = 1:Nwalkers
dXxfine = 0.01; Xxfinemax = 50; 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));
+ %pcolor(-tsfine,Xxfine,sqrt(smoothgrid));
%pcolor(tsfine,Xxfine,sqrt(histgrid)); %<<<BHJ: v?lge mellem glattet og r? farvel?gning
hold on
%plot a selection of depth histories
t@@ -917,7 +930,7 @@ html = [...
' </tr>\n'...
' <tr>\n'...
' <td>Rock density</td>\n'...
- ' <td>' num2str(rho) ' kg/m<sup>3</sup></td>\n'...
+ ' <td>' num2str(rock_density) ' kg/m<sup>3</sup></td>\n'...
' </tr>\n'...
' <tr>\n'...
' <td>ε<sub>gla</sub></td>\n'...
t@@ -940,8 +953,8 @@ html = [...
' </tr>\n'...
' <tr>\n'...
' <td>δ<sup>18</sup>O<sub>threshold</sub></td>\n'...
- ' <td>' num2str(record_treshold_min) ...
- ' to ' num2str(record_treshold_max) ' ‰</td>\n'...
+ ' <td>' num2str(record_threshold_min) ...
+ ' to ' num2str(record_threshold_max) ' ‰</td>\n'...
' </tr>\n'...
' <tr>\n'...
' <td>MCMC walkers</td>\n'...