Wednesday, December 17, 2008
Making SCC files for Celestia
SSC files are fairly simple text files. Here's the content of the 164P-Christensen.ssc file
# 164P C/2008 X4 C/2003 K2 Christensen
EllipticalOrbit # elements for epoch 2008 11 30
Epoch 2454820.4610 # 2008-Nov-30.0000000
Radius 10.45 # Hn=12.50 Hn=14.10-2.5*log(Albedo)-5*log(Radius) Albedo=0.04
Orientation [ 90 0 0 1 ] # random value
RotationPeriod 15.613880 # random value
Obliquity 77.387310 # random value
EquatorAscendingNode 239.190649 # random value
First off you can ignore everything under InfoURL, they are made up figures as newly discovered comets will not have this data. All you really need is the orbital elements, so where do you get these? From the Minor Planet Center, where else?
Now, you can get them in three ways, if the comet is really recent, you can look up the Minor Planet Electronic Circular. Otherwise, you can use either the Orbital elements of comets page, or you can use the Minor Planet Comet and Ephemeris Service. As a demonstration, I'll use the Ephemeris service. After clicking on the ephemeris service link, scroll down and type your comet in the box (use P/2008 X4 for this demonstration), then scroll down the Ephemeris Options below, ignoring all except the "Also display elements for epoch" box. Type the current date in there using year-month-day format 2008-12-18 (see image above, click image to enlarge). Ignore everything else and click on the Get Ephemeris button.
Down the bottom of the resulting page you will see:
Epoch (TT/JD) T (TT/JD) Peri. Node Incl. e q a Epoch (TT)Which look almost nothing like
2454818.50000 2454820.4610294 345.76463 93.89262 10.21562 0.8316725 0.5348100 3.17720 2008/12/18.00
Epoch 2454820.4610 # 2008-Nov-30.0000000
So Epoch = Epoch, that's clear UPDATE: No it's not, Epoch actually means Julian date of perihelion (that's why I had to use the second of the two dates). Except for P/2008 X4 I had to use the second of the two Epoch dates to get the comet in the right place.
Peri = ArgOfPericentre, Node = Ascending node, Incl = Inclination(so far so good), e? what the heck is e? e= eccentricity (hey obvious), q = PericenterDistance (only Celestia uses Pericenter, everyone else uses Perihelion), ignore a and use 0.0 for mean anomaly.
Where is the Period you ask? Good question, it's not there. But see the link up the top?
Perturbed ephemeris below is based on elements from MPEC 2008-X81.Click on the hyperlink, scroll down and just before the residuals you will find P 5.66. This is the period. This won't work for all comets, the old, established comets don't have that link but new ones will (obviously the name of the link will be different for other comets).
Alternatively you can use the JPL Horizons system (a PDF tutorial on using the system is here) In the JPL Horizons system Epoch=Epoch W = ArgOfPericentre, OM = Ascending node, IN = Inclination, EC= eccentricity, QR = PericenterDistance, PER= period and MA= mean anomoly (when given). However, the JPL Horizons system currently gives the WRONG orbital elements for P/2008 X4, it uses the old P/2003 K2 elements, which is why I didn't go there in the first place. If you are lucky your comet will be listed in Wikipedia (here is a list of all currently known period comets, but without C/2008 X4).
For known comets that aren't currently included in the Celestia main download, go to the Celestia Motherload comets page. There is a full catalog of periodic comets there (minus P/2008 X4 of course), but if you download that into extras and turn the orbits on, you won't see a thing for comet orbits. It may be best if you store this file somewhere else and extract comets as you need them.
Selden Ball has a list of instructions for converting MPC elements to Celestia format which is mots helpful (and also helps with other formats like the JPL Horizons format). He also has a neat list of Celestia resources which you might like.
Happy SSC file making (you can now do this for new asteroids and exoplanets as well)
Thought I'd drop you a line and say thank you for such a great blog. My girls and I have been really enjoying it. We got a great photo of the moon the other night thanks to your information.
Merry Christmas Ian
I'm pleased I've been ale to help you enjoy our great night skies! Don't forget December 29th, with the occultation of Jupiter. Hope you have a Merry and Safe Christmas and New year down at the farm!
Links to this post: