tshow table of input values - cosmo - front and backend for Markov-Chain Monte Carlo inversion of cosmogenic nuclide concentrations
HTML git clone git://src.adamsgaard.dk/cosmo
DIR Log
DIR Files
DIR Refs
DIR README
DIR LICENSE
---
DIR commit 49eb5cb6f34627f32882fe2d6c44a1d3d9d0700b
DIR parent ad3acb3166b42f0aadb967638f8ba59b2827fce9
HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Tue, 1 Dec 2015 12:11:41 +0100
show table of input values
Diffstat:
M index.php | 20 ++++++++++++++++++++
M matlab/generate_plots.m | 257 +++++++++++++------------------
2 files changed, 124 insertions(+), 153 deletions(-)
---
DIR diff --git a/index.php b/index.php
t@@ -290,6 +290,26 @@ if (isset($_GET['wait_id']) && !empty($_GET['wait_id'])) {
</div>
</div>
+ <div class="row">
+ <div class="col s12">
+ <div class="card">
+ <div class="card-image">
+ <?php
+ include('output/' . $_GET['results_id'] . '_Walks-input.html');
+ ?>
+ </div>
+ <div class="card-content">
+ <h5 class="blue-text">
+ Input parameters</h5>
+ <p>This table summarizes the supplied input parameter
+ values.</p>
+
+ </div>
+ </div>
+ </div>
+ </div>
+
+
<?php
} else {
?>
DIR diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
t@@ -825,178 +825,129 @@ html = [...
' <thead>\n'...
' <tr>\n'...
' <th data-field="param">Parameter</th>\n'...
- ' <th data-field="param">Percentile</th>\n'];
-
-for i=1:Nwalkers
- html = [html, ...
- ' <th data-field="w1">Walker ', num2str(i), '</th>\n'];
-end
-
-% epsilon_int
-html = [html, ...
- ' <th data-field="avg">Average</th>\n'...
+ ' <th data-field="val">Value</th>\n'...
' </tr>\n'...
' </thead>\n'...
' <tbody>\n'...
' <tr>\n'...
- ' <td> </td>\n'...
- ' <td align="center">25%%</td>\n'];
-for i=1:Nwalkers
- html = [html, ' <td>',...
- num2str(epsilon_int_25(i),3),'</td>\n'];
-end
-
-html = [html, ' <td>',...
- num2str(sum(epsilon_int_25)/Nwalkers,3),'</td>\n'...
+ ' <td>Sample ID</td>\n'...
+ ' <td>' sample_id '</td>\n'...
+ ' </tr>\n'...
+ ' <td>Name</td>\n'...
+ ' <td>' name '</td>\n'...
+ ' </tr>\n'...
+ ' <td>Email</td>\n'...
+ ' <td>' email '</td>\n'...
' </tr>\n'...
' <tr>\n'...
- ' <td align="center">ε<sub>int</sub> [m/Myr]</td>\n'...
- ' <td align="center">50%%</td>\n'];
-
-for i=1:Nwalkers
- html = [html, ' <td>',...
- num2str(epsilon_int_50(i),3),'</td>\n'];
-end
-
-
-html = [html, ' <td>',...
- num2str(sum(epsilon_int_50)/Nwalkers,3),'</td>\n'...
+ ' <td>Latitude</td>\n'...
+ ' <td>' num2str(lat) '</td>\n'...
' </tr>\n'...
- ' <tr style="border-bottom:1px solid #D0D0D0">\n'...
- ' <td> </td>\n'...
- ' <td align="center">75%%</td>\n'];
-
-for i=1:Nwalkers
- html = [html, ' <td>',...
- num2str(epsilon_int_75(i),3),'</td>\n'];
-end
-
-html = [html, ' <td>',...
- num2str(sum(epsilon_int_75)/Nwalkers,3),'</td>\n'...
- ' </tr>\n'];
-
-
-% epsilon_gla
-html = [html, ...
' <tr>\n'...
- ' <td> </td>\n'...
- ' <td align="center">25%%</td>\n'];
-for i=1:Nwalkers
- html = [html, ' <td>',...
- num2str(epsilon_gla_25(i),3),'</td>\n'];
-end
-
-html = [html, ' <td>',...
- num2str(sum(epsilon_gla_25)/Nwalkers,3),'</td>\n'...
+ ' <td>Longitude/td>\n'...
+ ' <td>' num2str(long) '</td>\n'...
' </tr>\n'...
' <tr>\n'...
- ' <td align="center">ε<sub>gla</sub> [m/Myr]</td>\n'...
- ' <td align="center">50%%</td>\n'];
-
-for i=1:Nwalkers
- html = [html, ' <td>',...
- num2str(epsilon_gla_50(i),3),'</td>\n'];
-end
-
-
-html = [html, ' <td>',...
- num2str(sum(epsilon_gla_50)/Nwalkers,3),'</td>\n'...
+ ' <td><sup>10</sup>Be concentration</td>\n'...
+ ' <td>' num2str(be_conc) ' atoms/g</td>\n'...
' </tr>\n'...
- ' <tr style="border-bottom:1px solid #D0D0D0">\n'...
- ' <td> </td>\n'...
- ' <td align="center">75%%</td>\n'];
-
-for i=1:Nwalkers
- html = [html, ' <td>',...
- num2str(epsilon_gla_75(i),3),'</td>\n'];
-end
-
-html = [html, ' <td>',...
- num2str(sum(epsilon_gla_75)/Nwalkers,3),'</td>\n'...
- ' </tr>\n'];
-
-
-% record_threshold
-html = [html, ...
' <tr>\n'...
- ' <td> </td>\n'...
- ' <td align="center">25%%</td>\n'];
-for i=1:Nwalkers
- html = [html, ' <td>',...
- num2str(record_threshold_25(i),3),'</td>\n'];
-end
-
-html = [html, ' <td>',...
- num2str(sum(record_threshold_25)/Nwalkers,3),'</td>\n'...
+ ' <td><sup>10</sup>Al concentration</td>\n'...
+ ' <td>' num2str(al_conc) ' atoms/g</td>\n'...
' </tr>\n'...
' <tr>\n'...
- ' <td align="center">δ<sup>18</sup>O<sub>threshold</sub> [‰]</td>\n'...
- ' <td align="center">50%%</td>\n'];
-
-for i=1:Nwalkers
- html = [html, ' <td>',...
- num2str(record_threshold_50(i),3),'</td>\n'];
-end
-
-
-html = [html, ' <td>',...
- num2str(sum(record_threshold_50)/Nwalkers,3),'</td>\n'...
+ ' <td><sup>10</sup>C concentration</td>\n'...
+ ' <td>' num2str(c_conc) ' atoms/g</td>\n'...
' </tr>\n'...
- ' <tr style="border-bottom:1px solid #D0D0D0">\n'...
- ' <td> </td>\n'...
- ' <td align="center">75%%</td>\n'];
-
-for i=1:Nwalkers
- html = [html, ' <td>',...
- num2str(record_threshold_75(i),3),'</td>\n'];
-end
-
-html = [html, ' <td>',...
- num2str(sum(record_threshold_75)/Nwalkers,3),'</td>\n'...
- ' </tr>\n'];
-
-% E
-html = [html, ...
' <tr>\n'...
- ' <td> </td>\n'...
- ' <td align="center">25%%</td>\n'];
-for i=1:Nwalkers
- html = [html, ' <td>',...
- num2str(E_25(i),3),'</td>\n'];
-end
-
-html = [html, ' <td>',...
- num2str(sum(E_25)/Nwalkers,3),'</td>\n'...
+ ' <td><sup>10</sup>Ne concentration</td>\n'...
+ ' <td>' num2str(ne_conc) ' atoms/g</td>\n'...
' </tr>\n'...
' <tr>\n'...
- ' <td align="center">E [m]</td>\n'...
- ' <td align="center">50%%</td>\n'];
-
-for i=1:Nwalkers
- html = [html, ' <td>',...
- num2str(E_50(i),3),'</td>\n'];
-end
-
-
-html = [html, ' <td>',...
- num2str(sum(E_50)/Nwalkers,3),'</td>\n'...
+ ' <td><sup>10</sup>Be conc. uncertainty</td>\n'...
+ ' <td>' num2str(be_uncer) ' %%</td>\n'...
' </tr>\n'...
- ' <tr style="border-bottom:1px solid #D0D0D0">\n'...
- ' <td> </td>\n'...
- ' <td align="center">75%%</td>\n'];
-
-for i=1:Nwalkers
- html = [html, ' <td>',...
- num2str(E_75(i),3),'</td>\n'];
-end
-
-html = [html, ' <td>',...
- num2str(sum(E_75)/Nwalkers,3),'</td>\n'...
- ' </tr>\n'];
-
-
-% footer
-html = [html, ' </tbody>\n'...
+ ' <tr>\n'...
+ ' <td><sup>10</sup>Al conc. uncertainty</td>\n'...
+ ' <td>' num2str(al_uncer) ' %%</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td><sup>10</sup>C conc. uncertainty</td>\n'...
+ ' <td>' num2str(c_uncer) ' %%</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td><sup>10</sup>Ne conc. uncertainty</td>\n'...
+ ' <td>' num2str(ne_uncer) ' %%</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td>Observation depth</td>\n'...
+ ' <td>' num2str(zobs) ' m</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td><sup>10</sup>Be production (spallation)</td>\n'...
+ ' <td>' num2str(be_prod_spall) ' atoms/g/yr</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td><sup>10</sup>Al production (spallation)</td>\n'...
+ ' <td>' num2str(al_prod_spall) ' atoms/g/yr</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td><sup>10</sup>C production (spallation)</td>\n'...
+ ' <td>' num2str(c_prod_spall) ' atoms/g/yr</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td><sup>10</sup>Ne production (spallation)</td>\n'...
+ ' <td>' num2str(ne_prod_spall) ' atoms/g/yr</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td><sup>10</sup>Be production (muons)</td>\n'...
+ ' <td>' num2str(be_prod_muons) ' atoms/g/yr</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td><sup>10</sup>Al production (muons)</td>\n'...
+ ' <td>' num2str(al_prod_muons) ' atoms/g/yr</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td><sup>10</sup>C production (muons)</td>\n'...
+ ' <td>' num2str(c_prod_muons) ' atoms/g/yr</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td><sup>10</sup>Ne production (muons)</td>\n'...
+ ' <td>' num2str(ne_prod_muons) ' atoms/g/yr</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td>Rock density</td>\n'...
+ ' <td>' num2str(rho) ' kg/m<sup>3</sup></td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td>ε<sub>gla</sub></td>\n'...
+ ' <td>' num2str(epsilon_gla_min) ...
+ ' to ' num2str(epsilon_gla_max) ' m/Myr</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td>ε<sub>int</sub></td>\n'...
+ ' <td>' num2str(epsilon_int_min) ...
+ ' to ' num2str(epsilon_int_max) ' m/Myr</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td><i>t</i><sub>degla</sub></td>\n'...
+ ' <td>' num2str(t_degla_min) ...
+ ' to ' num2str(t_degla_max) ' yr</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td>Climate record</td>\n'...
+ ' <td>' record '</td>\n'...
+ ' </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'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td>MCMC walkers</td>\n'...
+ ' <td>' num2str(nwalkers) '</td>\n'...
+ ' </tr>\n'...
+ ' </tbody>\n'...
'</table>\n'...
];
fileID = fopen(filename,'w');