tcall inversion 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 64900976c865f661a74355d368768e4479cd317a
(DIR) parent b78530b6b2d91c57ce162c5f2512139c2fab7a93
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Wed, 26 Aug 2015 13:09:01 +0200
call inversion function
Diffstat:
M matlab/file_scanner_mcmc_starter.m | 22 +++++++++++++++++++---
1 file changed, 19 insertions(+), 3 deletions(-)
---
(DIR) diff --git a/matlab/file_scanner_mcmc_starter.m b/matlab/file_scanner_mcmc_starter.m
t@@ -18,10 +18,15 @@ prefix = 'cosmo_';
% information
archivefolder = '/Users/ad/tmp/cosmo-archive';
+% folder containing matlab scripts to path
+matlab_scripts_folder = 'm_pakke2014maj11/';
%% general settings
debug = true; % show debugging output to matlab console
+%% initialization
+addpath(matlab_scripts_folder);
+
%% main loop
while 1
t@@ -55,14 +60,25 @@ while 1
record_threshold_min, record_threshold_max] ...
= import_php_file(infile, 1, 1); % only read first line
+ % run inversion
+ 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, ...
+ rock_density, ...
+ epsilon_gla_min, epsilon_gla_max, ...
+ epsilon_int_min, epsilon_int_max, ...
+ t_degla, t_degla_uncer, ...
+ record, ...
+ record_threshold_min, record_threshold_max);
+
% delete or archive the file so it is not processed again
%delete(infile)
%movefile(infile, archivefolder);
-
-
+
+ %keyboard
end
-
% for debugging purposes; ends loop after first iteration
break