First, download and extract the sec4code.tgz archive.

To reproduce the actual MCMC chains and posterior distributions for the results in section 4, modify the "sec4run.m" MATLAB script for your location and then run that script.

After that has concluded, modify the "plotres.R" R script for your location and run that script to generate the final PDF.

NOTE that the full set of results will consume about 75GB of disk storage.