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 49eb5cb6f34627f32882fe2d6c44a1d3d9d0700b
parent ad3acb3166b42f0aadb967638f8ba59b2827fce9
Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date:   Tue,  1 Dec 2015 12:11:41 +0100

show table of input values

Diffstat:

Mindex.php | 20++++++++++++++++++++


Mmatlab/generate_plots.m | 257++++++++++++++++++++++++++++++++-----------------------------------------------


2 files changed, 124 insertions(+), 153 deletions(-)

diff --git a/index.php b/index.php
@@ -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 {
 ?>
diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
@@ -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>&nbsp;</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">&epsilon;<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>&nbsp;</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>&nbsp;</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">&epsilon;<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>&nbsp;</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>&nbsp;</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">&delta;<sup>18</sup>O<sub>threshold</sub> [&permil;]</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>&nbsp;</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>&nbsp;</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>&nbsp;</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>&epsilon;<sub>gla</sub></td>\n'...
+    '            <td>' num2str(epsilon_gla_min) ...
+    ' to ' num2str(epsilon_gla_max) ' m/Myr</td>\n'...
+    '        </tr>\n'...
+    '        <tr>\n'...
+    '            <td>&epsilon;<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>&delta;<sup>18</sup>O<sub>threshold</sub></td>\n'...
+    '            <td>' num2str(record_treshold_min) ...
+    ' to ' num2str(record_treshold_max) ' &permil;</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');