treturn data from inversion, try plotting results - 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 18f347b2531e003200eb970faa5c465cb5d0cfa8 (DIR) parent 24e0cd3ec43d514f0a9d7bf66d80f9660f4c8824 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk> Date: Wed, 26 Aug 2015 15:17:55 +0200 return data from inversion, try plotting results Diffstat: M matlab/file_scanner_mcmc_starter.m | 5 ++++- M matlab/mcmc_inversion.m | 15 ++++++++++++--- 2 files changed, 16 insertions(+), 4 deletions(-) --- (DIR) diff --git a/matlab/file_scanner_mcmc_starter.m b/matlab/file_scanner_mcmc_starter.m t@@ -61,7 +61,7 @@ while 1 = import_php_file(infile, 1, 1); % only read first line % run inversion - mcmc_inversion(matlab_scripts_folder, debug, ... + [Ss, save_file] = mcmc_inversion(matlab_scripts_folder, debug, ... be_conc, al_conc, c_conc, ne_conc, ... be_uncer, al_uncer, c_uncer, ne_uncer, ... be_prod, al_prod, c_prod, ne_prod, ... t@@ -72,6 +72,9 @@ while 1 record, ... record_threshold_min, record_threshold_max); + % generate plots + CompareWalks2(Ss, save_file) + % delete or archive the file so it is not processed again %delete(infile) %movefile(infile, archivefolder); (DIR) diff --git a/matlab/mcmc_inversion.m b/matlab/mcmc_inversion.m t@@ -1,4 +1,4 @@ -function [] = mcmc_inversion(matlab_scripts_folder, debug, ... +function [Ss, save_file] = mcmc_inversion(matlab_scripts_folder, debug, ... be_conc, al_conc, c_conc, ne_conc, ... be_uncer, al_uncer, c_uncer, ne_uncer, ... be_prod, al_prod, c_prod, ne_prod, ... t@@ -98,8 +98,17 @@ switch fs.g_case fs.tStarts = NaN; %load or compute fixed times of more or less glaciated periods fs.relExpos = NaN; %load or compute degree of exposure in periods case 'd18OTimes' - % fs.d18Ofn = 'lisiecki_triinterp_2p6Ma_5ky.mat'; - fs.d18O_filename = 'lisiecki_triinterp_2p6Ma_30ky.mat'; % zachos_triinterp_2p6Ma + %fs.d18Ofn = 'lisiecki_triinterp_2p6Ma_5ky.mat'; + %fs.d18O_filename = 'lisiecki_triinterp_2p6Ma_30ky.mat'; % zachos_triinterp_2p6Ma + if strcmp(record, 'rec_5kyr') + fs.d18O_filename = 'lisiecki_triinterp_2p6Ma_5ky.mat'; % zachos_triinterp_2p6Ma + elseif strcmp(record, 'rec_20kyr') + fs.d18O_filename = 'lisiecki_triinterp_2p6Ma_20ky.mat'; % zachos_triinterp_2p6Ma + elseif strcmp(record, 'rec_30kyr') + fs.d18O_filename = 'lisiecki_triinterp_2p6Ma_30ky.mat'; % zachos_triinterp_2p6Ma + else + error(['record ' record ' not understood']); + end fs.tStarts = NaN; %load or compute fixed times of more or less glaciated periods fs.relExpos = NaN; %load or compute degree of exposure in periods end