tgenerate output files with ascii data - cosmo - front and backend for Markov-Chain Monte Carlo inversion of cosmogenic nuclide concentrations

git clone git://src.adamsgaard.dk/cosmo

Log

Files

Refs

README

LICENSE

---

commit 5f1e1cff49491cbb75276a6ed084c916bcbe0aa4

parent 64ca77a04780213f339c748f3f89e0221dbafe90

Author: Anders Damsgaard 

Date:   Mon,  7 Dec 2015 15:44:19 +0100

generate output files with ascii data

Diffstat:
  M matlab/generate_plots.m             |      10 ++++++++++

1 file changed, 10 insertions(+), 0 deletions(-)
---

diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m

t@@ -338,12 +338,16 @@ for i1 = 1:M % for each model parameter
     % save median values for later
     if i1 == 1
         epsilon_int_med = med;
+        epsilon_int_data = data;
     elseif i1 == 2
         epsilon_gla_med = med;
+        epsilon_gla_data = data;
     elseif i1 == 3
         t_degla_med = med;
+        t_degla_data = data;
     elseif i1 == 4
         record_threshold_med = med;
+        record_threshold_data = data;
     else
         error('Unknown parameter');
     end
t@@ -1090,3 +1094,9 @@ html = ['\n' ...
 fileID = fopen(filename,'w');
 fprintf(fileID, html);
 fclose(fileID);
+
+%% Save general data
+dlmwrite([save_file, '-eps_int.txt'], epsilon_int_data);
+dlmwrite([save_file, '-eps_gla.txt'], epsilon_gla_data);
+dlmwrite([save_file, '-t_degla.txt'], t_degla_data);
+dlmwrite([save_file, '-d18O_threshold.txt'], record_threshold_data);