file='Rm1_o16.dat'
run='XX'
t1=30.
t1b=100.
t2=1e9