Rasmus Handberg

Astronomer, Software Developer...

Timeseries Analysis Tools

The programs are written in Fortran 90 and uses OpenMP for parallisation

.

Installation instructions

Using existing installation at SAC, Aarhus University

If you are and have access to the computer systems here, the easiest way to use the programs is simply to use the pre-compiled binaries available on the systems.

/usr/users/f041080/tidsserie/bin/machine

where machine is the system you are running on.

If you wish the program to use all available CPU cores available on your machine (and why wouldn't you?), you need to edit the file ~/.bashrc or ~/.cshrc (depending on which shell you use) and add the following line:

setenv OMP_NUM_THREADS 8 (csh, tcsh)
export OMP_NUM_THREADS=8 (bash)

The entire setup for the computer "ast", with 8 CPU cores, would for example look like this (in bash):

export PATH=$PATH:/usr/users/f041080/tidsserie/bin/ast
export OMP_NUM_THREADS=8

If you have any questions, feel free to ask me.

Compiling it yourself

First of all, you should send me an e-mail. Then I will grant you access to the SVN repository where the code is stored, and you should then be able to check out the code from the repo. using the following command:

svn checkout https://xp-dev.com/svn/TimeseriesTools .
make; make install

This will compile the program (by default using gfortran) and put the compiled binaries into the "bin" directory. If you wish to use a different compiler, you can simply modify the file "makefile.preample" where you can specify which compiler (and settings) to use.

Documentation

Input files should be ASCII files with three columns: time, measurement and statistical error.

spec

Weighted Least Squares Spectrum

Input:
  spec [options] [inputfile] [outputfile]

Options:
  -p         Power spectrum (default)
  -a         Amplitude spectrum
  -pd        Power density spectrum
  -tsec      Treat times as seconds (default)
  -tday      Treat times as days
  -kplrraw   Input file is a Kepler data file (use raw columns)
  -kplrpdc   Input file is a Kepler data file (use PDC columns)
  -kplrwg    Input file is a Kepler data file (use WG columns)
  -auto      Automatic. No interaction required (experimental)
  -quiet     Print nothing to screen
  -version   Print program version information

Written by
  Rasmus Handberg
  Department of Physics and Astronomy
  Aarhus University

window

WINDOW
Calculate the Spectral window-function
of a timeseries.

Input:
   window [options] [inputfile] [outputfile]

Options:
   -tday      Treat times as days (default)
   -tsec      Treat times as seconds
   -kplrraw   Input file is a Kepler data file (use raw columns)
   -kplrpdc   Input file is a Kepler data file (use PDC columns)
   -kplrwg    Input file is a Kepler data file (use WG columns)
   -quiet     Print nothing to screen
   -version   Print program version information

Written by
  Rasmus Handberg
  Department of Physics and Astronomy
  Aarhus University

clean

CLEAN
Iterative Sine Wawe Fitting

Input:
   clean [options] [inputfile] [outputfile]

Options:
   -tday      Treat times as days (default)
   -tsec      Treat times as seconds
   -version   Print program version information

Written by
  Rasmus Handberg
  Department of Physics and Astronomy
  Aarhus University

lclean

LCLEAN
Iterative Sine Wawe Fitting

Description:
  Cleans a timeseries until given limit.

Input:
  lclean [options] [inputfile] [outputfile]

Options:
  -p         Power spectrum (default)
  -a         Amplitude spectrum
  -pd        Power density spectrum
  -tday      Treat times as days (default)
  -tsec      Treat times as seconds
  -version   Print program version information

Written by
  Rasmus Handberg
  Department of Physics and Astronomy
  Aarhus University

sclean

SCLEAN
Simultaneous Iterative Sine Wave Fitting

Description:
  Cleans a timeseries for coherent signals,
  taking into account the influences of the
  frequencies on each other.

Input:
  sclean [options] [inputfile] [outputfile]

Options:
  -tday      Treat times as days (default).
  -tsec      Treat times as seconds.
  -version   Print program version information.

Written by
  Rasmus Handberg
  Department of Physics and Astronomy
  Aarhus University

lowpass

LOWPASS
Lowpass filtering using Weighted Least Squares

Input:
  lowpass [options] [inputfile] [outputfile]

Options:
  -tday      Treat times as days (default).
  -tsec      Treat times as seconds.
  -version   Print program version information.

Written by
  Rasmus Handberg
  Department of Physics and Astronomy
  Aarhus University

highpass

HIGHPASS
Highpass filtering using Weighted Least Squares

Input:
  highpass [options] [inputfile] [outputfile]

Options:
  -tday      Treat times as days (default)
  -tsec      Treat times as seconds
  -version   Print program version information

Written by
  Rasmus Handberg
  Department of Physics and Astronomy
  Aarhus University

bandpass

bandstop