TTV model output

[macgrac]

Dear Stefano,

The plot(k, type='ttv') function produces great figures, but I need to export the full TTV model (sketched with a red line in the plot) to the ascii file. In the 'Transit timings' panel, one can export model values only for times of observations. Is there any command to access the full model data?

Thanks, Gracjan (Systemic 2.12)

[stefano-meschiari]

Hi Gracjan,

You can get all the transits in a given time interval by using the kintegrate function:

# Assuming you have already loaded the planets in the "k" kernel and are interested in transits
# within the span of observations

i <- kintegrate(k, times=seq(from=k$trange[1], to=k$trange[2], by=min(k[, 'period'])), transits=TRUE)

The "i" now contains the results of the integration (done with the currently-selected integration method). You can check out the contents of the object by typing str(i).

You can print the transits times using

i$transits[[1]]
# or, with more digits:
sprintf("%.10f", i$transits[[1]])

which is a vector containing the times of primary transits (/not/ corrected for light travel). The number between double square brackets is the number of the planet transiting. You can save those times to file:

write(sprintf("%.10f", i$transits[[1]]), file="mytransits.txt")

I'll add a menu item/button to achieve this in a future release (together with the ability to print secondary eclipses as well).

Hope this helps! -Stefano

[macgrac]

Dear Stefano,

Many thanks! It works great!

Best, Gracjan

WiadomoϾ napisana przez Stefano Meschiari w dniu 6 lut 2014, o godz. 18:13:

Hi Gracjan,

You can get all the transits in a given time interval by using the kintegrate function:

Assuming you have already loaded the planets in the "k" kernel and are interested in transits

within the span of observations

i <- kintegrate(k, times=seq(from=k$trange[1], to=k$trange[2], by=min(k[, 'period'])), transits=TRU E) The "i" now contains the results of the integration (done with the currently-selected integration method). You can print the transits times, using

i$transits[[1]]

or, with more digits:

sprintf("%.10f", i$transits[[1]]) which is a vector containing the times of primary transits (/not/ corrected for light travel). The number between double square brackets is the number of the planet transiting. You can save those times to file:

write(sprintf("%.10f", i$transits[[1]]), file="mytransits.txt") I'll add a menu item/button to achieve this in a future release (together with the ability to print secondary eclipses as well).

Hope this helps! -Stefano

— Reply to this email directly or view it on GitHub.

[macgrac]

Dear Stefano,

Many thanks for the new version of Systemic. The precession rate is a great feature. I noticed however that it works in the Keplerian approach only, not with dynamical integrators.

Best, Gracjan