GitLab

Provide all essential building blocks for the toy example.
Still missing: proper step-by-step description.

A consitent way to use file name patterns is essential to support the addition
of further future examples.

Offset variable tapers are defined be picking times and offsets in a
seismogram display. It is quite likely that some of these picks, in particular
for near-source traces, are not within the time window of seismic sample data.
In such cases the four-point taper complaines and aborts. The offset variable
taper now ensures reasonable values to be passed to the four-point taper in
any case, by applying lower and upper bounds to the time values.

- provide tapers for muting
- provide muting operation
- provide arrival time picks
- provide refracted wave travel times and model
- provide plots for all of them

- calculate seismograms with reflectivity method
- add gaussian noise

The example will be based on synthetic data. Provide refmet config files.

Conflicts:
	src/green/gremlin1/Makefile

copy literally from TFtexte/pjt/distel/abb/synt/2lay

Use proper libtsioxx output operators.

soutifuamp.gpt was accidentally deleted by target simulclean

HERE documents are cannot be defined within make rules. A work-around are
multi-line make variables exported to the environment and printed to file
within the shell. This way explanatory texts are easier to maintain and are
even easier to read within the Makefile

- provide detailed documentation
- add noise to synthetic test data
- provide control plots

Multiple calls to ts::rnd::dugauss() resulted in identical sequences of random
numbers.

Issue regarding random number generator in libtsxx
==================================================

Location: src/libs/libtsxx/random.cc

The function ts::rnd::dugauss() allocates a GSL random number generator and
immediately seeds it. The seeding operation is required, since the default
seed is a constant, which would result in identical sequences. Function
ts::rnd::dugauss() uses the current system time as returned by function time()
as a seed value. Unfortunately, if ts::rnd::dugauss() is called several times
during the same second, this also will result in identical sequences.

Currently ts::rnd::dugauss() is only used by src/ts/wf/noisymize.cc and by the
filter GaussianNoise in libtsxx.

The appropriate solution appears to be:
- encapsulate the GSL randon number generator in a class object, just like is
  provided in tfxx::numeric::RNGgaussian
- use a static variable of this class type in ts::rnd::dugauss()
  this variable will be initialized only once upon the first call to the
  function
- the destructor is guaranteed to be called after termination of the program,
  such that the destructor is able to properly call gsl_rng_free

The filter GaussianNoise is provided in tidofi through keyword noi.
It adds white Gaussian Noise of a given rms amplitude (specified through
filter parameter) to each time series.

The filter is implemented in libtsxx. tidofi is additionally linked against
libgsl to provide random number generator.