//Run4bMigBS with dissipation >> Model of paper GV-EPJInitConfig file for 4B problem with dissipatiob 
//VERSION 3
//Press Z during run to set dissipative forces to zero
//Press H during run to set dissipative force amplitude (C[i]) to half

workfile = k51
integration accuracy = 1.0e-15
masses (m1,m2,m3,m0) = 0.0000066  0.0000126   0.0000239   1.0

//dissipation params index 0 is free param, index 1...3 means planet
//db[i] params are devided by mass[i] to give C[i] in EOM of EPJ-paper, i.e. C[i]::db[i]/m[i]
db[i]  = 0.0000  0.0e-10  0.37e-10  1.0e-10
da[i]  = 0.0000  1.00     0.99     0.99

//Initial orbital elements of 3 planets (normed a of kepler-51 : c 1.5274463   d 2.0246619 
a=  1.0   	1.7     	2.3
e=  0.01       0.008           0.01
i=  0.0   	0.0       	0.0
om= 0.0  	0.0      	0.0
Om= 0.0  	0.0      	0.0
M=  0.0  	0.0        	0.0

//time parameters
Iterations N=500000;
Iteration time step DT = 6.2831853
starting step dt=0.001;

Screen Out modulo = 200
Data file out modulo = 50
Time-Scale factor for output = 0.0001233 //Kyears for kepler-51
Target resonance (kepler-51)
Tb/Tc (p12/q12)  = 1 2
Tc/Td (p23/q23)  = 2 3
=========== END ==========================