A Large Example
As an example of how Bongo can manage a large number of data files, mixing
fortran code and Bongo commands, an application to 73 datasets from
Supplementary Material (SM) of Marcott et al.,2013
(see also the WUWT site).
By using the proxy data in the SM, Marcott et
al. derive the presence of a large "spike" after industrial revolution (a
novel "Hockey Stick").
Can you understand as such dispersed and often opposite data can give the
spike shown here or here?
Due to this example, the dimension of the
command line buffer in Bongo has been extended to 1000 lines.
The situation is the following:
- A .xlsx file is available, containing 73 datasets (and much more material) with
proxies for global temperatures, from 20000 to 0 years before present (B.P.,
present set to 1950).
- 73 manually extracted (and corrected from 'comma' to 'dot'-separated
decimal digits) files (name-of-proxy.dat). The two colums of each file are
- I need to plot all the datasets into a single frame (graph).
- I don't like (I hate!) to write 8-10 lines of Bongo commands for any file to
(580-730 lines and complex file names!).
- Bongo can read tab-separated data, but my fortran code for power
spectrum does not. So, I need to write datasets as blank-separated-columns
files. Of course the Bongo command POWER can manage tab-separated, but ...
- Prepared a file (list.txt) with all
file-names.ext of the datasets by the shell commands:
|ls -l *.dat>dum.txt ||!list directory to the file dum.txt
|cut -c 48-69 dum.txt >list.txt || !extract columns 48 to 69
from dum.txt to list.txt
- Written a fortran code to prepare the plotting of all datasets, changing
colors. This code reads list.txt and produces a separate command
file, to be inserted in the main Bongo command file (marcott-all.bon)
- marcott-all.bon is executed to produce the plot (pdf), (png)
- Written a fortran code to manage
datafiles as blank-separated-columns output files (extension .ntb, i.e. no-tab).
This code writes a separate Bongo command
file, to be inserted in the main Bongo command file (marcott2-all.bon)
- Running marcott2-all.bon produces the .ntb
Page Written: March 17,