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? I cannot.

Due to this example, the dimension of the command line buffer in Bongo has been extended to 1000 lines.

The situation is the following:

  1. 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).
  2. 73 manually extracted (and corrected from 'comma' to 'dot'-separated decimal digits) files (name-of-proxy.dat). The two colums of each file are 'tab'-separated.
  3. I need to plot all the datasets into a single frame (graph).
  4. I don't like (I hate!) to write 8-10 lines of Bongo commands for any file to be plotted (580-730 lines and complex file names!).
  5. 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 ...

My solution:

  1. 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
  2. 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)
  3. marcott-all.bon is executed to produce the plot (pdf), (png)
  4. 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)
  5. Running marcott2-all.bon produces the .ntb files

Page Written: March 17, 2013        Last Updated: