// app3.ini
// Phase volume during Kepler oscillator evolution
// simultaneously computed with the Direct Midpoint Method (yellow dots)
// and the Runge Kutta method of second order (white dots)
// New controls added 2009-08-03.
//////////////////////////////////////////////////////////////////////////
data
//////////////////////////////////////////////////////////////////////////
R E=0.0
// Initial value of the 'excentric anomaly'
//R eps=0.05
R eps=0.10
// numerical excentricity, e.g. 0.17
R nPerRev=50
// computed points per revolution (needs not to be integer)
// e.g. 32
//Z nTot=200
//Z nTot=1500
//Z nTot=16000
Z nTot=15879
// total number of computed steps (per system)
// e.g. 3600
Z nK=4000
// closed curve is made of so many points
R fac=1.1
//R fac=1.1
// factor by which the axes of the display area are larger than the
// exact orbit, e.g. 2.2.
// Needs to be considerably larger than 1, if the Runge Kutta dots
// are to stay within the frame after a few revolutions
R xrs=0.
// x-shift of the graphical window in fractions of its width
R yrs=-0.3
// y-shift of the graphical window in fractions of its height
R acc=1e-10
// accuracy used in solving Keplers equation by iteration.
R tWaitFinal1=6
// time in s for inspecting the final phase space curve
R tWaitFinal2=6
// time in s for inspecting the phase area development curve
// This is a new feature. Phase area is virtualy constant for DALF
// and linearly decreasing for ADALF
R tWaitFrame=0
// time for viewing the individual frame
Z method=4
// 1 Euler, 2 Runge-Kutta, 3 DMI, 4 ALF, 5 DALF, 6 LF, 7 ADALF, 8 AALF, 9 HAM
// DALF was added later. Earlier it was accessed as
// Z method=2
// B stepRed_=1
Z displayPeriod=1
// after as many steps are computed, they get transfered to
// the display screen. Since displaying takes time, having values
// >1 here lets the program run much faster
// e.g. 50
//Zs filingLimits=0 0
Zs filingLimits=15800 15879
// lowest and highest frame number to be written to file
B fileOnShow=1
// only if this is true, filing of screen frames will happen.
// The last frame will be then filed anyway, other frames only
// if their frame number is in the range given by the
// filingLimits.
R facX=1.
R facV=1.
// controls the size of the phase space subset the boundary of which
// will be visualized.
W imageName=kneading
B backEvl=1
// back evolution (corresponds to interaction picture).
// If this is true the computed position is evolved back
// by the reversible exact dynamics. Otherwise we have
// the computed dynamics followed by an operation of centering.
W movieName=alf