## filename of output file
bin.proper.dat

## FILTER TYPE selection
3	/* 0) switched off; 1) run. avg.; 2) run. min/max; 3) min. a/max. e; 4) representative surface; 5) FMFT */

## Write real*8 binary output instead of real*4? [T|F]
F

## sampling time-step should be an integer multiple of dtfilter!
## (refer to filter.in file) - osculating elements are calculated
## only every dtfilter (thus it is the maximum sampling rate)
-1	/* sampling time-step [day], set to -1 for equal to dtfilter */

## options for running-average filter
5.0 	/* width of running window [Myr] */
0.1  	/* output time-step [Myr] */

## options for running-minimum/maximum filter of the given elements,
## the output of others is for the respective time of min./max. -
## ie. far from resonances, where a, e oscillations are NOT coupled,
## it produces misleading results; the ABOVE running window options apply also!
1	/* orbital element to be monitored: 1 a, 2 e, 3 inc, 4 node, 5 peri */
T	/* T minimum, F maximum */

## options for INDIVIDUAL minimum/maximum filter of a, e, i
## (the above running window options apply also!)
T	/* semimajor axis - T minimum, F maximum */
T	/* eccentricty */
T	/* inclination */

## options for representative surface filter (so called pseudo-proper,
## or resonant elements); the above running window options apply also!
## a critical angle of an arbitrary mean motion resonance is:
##   sigma = (p+q)/q*lambda_PL - p/q*lambda_TP - varpi_TP
## a representative surface is defined by conditions:
##   sigma ~ sigma_crit, d sigma/d t > or < 0 and eventually also
##   varpi-varpi_PL ~ varpi_crit, Omega-Omega_PL ~ Omega_crit
-2	/* identification of the planet (PL id) */
1	/* p */
1	/* q */
0.	/* critical value of sigma [deg] */
5.	/* max. value of abs(sigma-sigma_crit) [deg] */
0.	/* critical value of varpi-varpi_PL */
5.	/* maximal difference from ^^^, set to -1 to switch the check off! */
180.	/* critical value of Omega-Omega_PL */
-1.	/* max. diff. from ^, set to -1 to switch the check off! */
-1	/* signum of d sigma/d t [-1|1], set to 0 to switch it off! */

## options for FMFT filter
2	/* the variant of frequency analysis, see fmft.c */   
10	/* number of computed frequencies */
8192	/* number of used data, must be a power of 2! */ 
T	/* write `fmft.out' file with TPid, time, g and s frequencies (in arcsec/yr) */

-1000.	/* minimum frequency (arcsec/yr), frequencies smaller than MIN_F are ignored */ 
1000.  	/* maximum frequency */ 
0.02	/* absolute value of the smallest allowed frequency */
0.03	/* allowed difference between calculated and plan. freq., see below */
7	/* number of (known, forced) planetary frequencies to be discarded */

  4.25749319	/* list of plan. freq.: g5 */
 28.24552984	/* g6 */
  3.08675577	/* g7 */
  0.67255084	/* g8 */
-26.34496354	/* s6 */
 -2.99266093	/* s7 */
 -0.69251386	/* s8 */

-1000.	/* minimum g frequency (arcsec/yr), otherwise, it is considered planetary */
1000.	/* maximum g */
-1000.	/* minimum s */
1000.	/* maximum s */
F	/* add vectorially all frequencies matching the above criteria */
T	/* drop intervals, where an interrupt (a longer time-step) occurs */
F	/* calculate also proper elements of planets */
F	/* write all frequencies, amplitudes and phases to `propgs.out' */

