orbital_dt function

Convert to data.table code