itgenerate output files with ascii data - cosmo - front and backend for Markov-Chain Monte Carlo inversion of cosmogenic nuclide concentrations Err mx1.adamsgaard.dk 70 hgit clone git://src.adamsgaard.dk/cosmo URL:git://src.adamsgaard.dk/cosmo mx1.adamsgaard.dk 70 1Log /src/log.gph mx1.adamsgaard.dk 70 1Files /src/files.gph mx1.adamsgaard.dk 70 1Refs /src/refs.gph mx1.adamsgaard.dk 70 1README /src/file/README.md.gph mx1.adamsgaard.dk 70 1LICENSE /src/file/LICENSE.gph mx1.adamsgaard.dk 70 i--- Err mx1.adamsgaard.dk 70 1commit 5f1e1cff49491cbb75276a6ed084c916bcbe0aa4 /src/commit/5f1e1cff49491cbb75276a6ed084c916bcbe0aa4.gph mx1.adamsgaard.dk 70 1parent 64ca77a04780213f339c748f3f89e0221dbafe90 /src/commit/64ca77a04780213f339c748f3f89e0221dbafe90.gph mx1.adamsgaard.dk 70 hAuthor: Anders Damsgaard URL:mailto:anders.damsgaard@geo.au.dk mx1.adamsgaard.dk 70 iDate: Mon, 7 Dec 2015 15:44:19 +0100 Err mx1.adamsgaard.dk 70 i Err mx1.adamsgaard.dk 70 igenerate output files with ascii data Err mx1.adamsgaard.dk 70 i Err mx1.adamsgaard.dk 70 iDiffstat: Err mx1.adamsgaard.dk 70 i M matlab/generate_plots.m | 10 ++++++++++ Err mx1.adamsgaard.dk 70 i Err mx1.adamsgaard.dk 70 i1 file changed, 10 insertions(+), 0 deletions(-) Err mx1.adamsgaard.dk 70 i--- Err mx1.adamsgaard.dk 70 1diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m /src/file/matlab/generate_plots.m.gph mx1.adamsgaard.dk 70 it@@ -338,12 +338,16 @@ for i1 = 1:M % for each model parameter Err mx1.adamsgaard.dk 70 i % save median values for later Err mx1.adamsgaard.dk 70 i if i1 == 1 Err mx1.adamsgaard.dk 70 i epsilon_int_med = med; Err mx1.adamsgaard.dk 70 i+ epsilon_int_data = data; Err mx1.adamsgaard.dk 70 i elseif i1 == 2 Err mx1.adamsgaard.dk 70 i epsilon_gla_med = med; Err mx1.adamsgaard.dk 70 i+ epsilon_gla_data = data; Err mx1.adamsgaard.dk 70 i elseif i1 == 3 Err mx1.adamsgaard.dk 70 i t_degla_med = med; Err mx1.adamsgaard.dk 70 i+ t_degla_data = data; Err mx1.adamsgaard.dk 70 i elseif i1 == 4 Err mx1.adamsgaard.dk 70 i record_threshold_med = med; Err mx1.adamsgaard.dk 70 i+ record_threshold_data = data; Err mx1.adamsgaard.dk 70 i else Err mx1.adamsgaard.dk 70 i error('Unknown parameter'); Err mx1.adamsgaard.dk 70 i end Err mx1.adamsgaard.dk 70 it@@ -1090,3 +1094,9 @@ html = ['\n' ... Err mx1.adamsgaard.dk 70 i fileID = fopen(filename,'w'); Err mx1.adamsgaard.dk 70 i fprintf(fileID, html); Err mx1.adamsgaard.dk 70 i fclose(fileID); Err mx1.adamsgaard.dk 70 i+ Err mx1.adamsgaard.dk 70 i+%% Save general data Err mx1.adamsgaard.dk 70 i+dlmwrite([save_file, '-eps_int.txt'], epsilon_int_data); Err mx1.adamsgaard.dk 70 i+dlmwrite([save_file, '-eps_gla.txt'], epsilon_gla_data); Err mx1.adamsgaard.dk 70 i+dlmwrite([save_file, '-t_degla.txt'], t_degla_data); Err mx1.adamsgaard.dk 70 i+dlmwrite([save_file, '-d18O_threshold.txt'], record_threshold_data); Err mx1.adamsgaard.dk 70 .