# Missing semi-major axis column in MPCORB table of DP0.3?

**URL:** <https://www.rubin.community/t/missing-semi-major-axis-column-in-mpcorb-table-of-dp0-3/7801>\
**Category:** Support\
**Tags:** solarsystem, dp0, dp03\
**Created:** [August 2, 2023, 11:09pm UTC](https://www.rubin.community/t/missing-semi-major-axis-column-in-mpcorb-table-of-dp0-3/7801 "2023-08-02T23:09:25Z")\
**Posts on this page:** 4\
**Page:** 1

<div class="post-metadata">

**Author:** ![talister](https://sea2.discourse-cdn.com/flex002/user_avatar/www.rubin.community/talister/32/875_2.png) [@talister](https://www.rubin.community/u/talister)\
**Post date:** [August 2, 2023, 11:09pm UTC](https://www.rubin.community/t/missing-semi-major-axis-column-in-mpcorb-table-of-dp0-3/7801/1 "2023-08-02T23:09:25Z")

</div>

Hi folks, thanks so much for making the simulated solar system objects and DP0.3 available. As a first time user of the RSP, I tried the tutorial to use the `MPCORB` table (dp03\_catalogs\_10yr.MPCORB) through the TAP service. Using the suggested cut of H\<20 produced me a table of the 50,000 (limit) objects as expected. However it seemed that the semi-major (‘a’) column which I would expect and which should be there according to Table 4/page 31 of the DPDD, is missing and is also missing in the schema browser. Can do the remedial orbital elements math and reconstruct it from the other parameters but this seems like it should be there… ?  
Ta muchly,  
Tim

---

<div class="post-metadata">

**Author:** ![MelissaGraham](https://sea2.discourse-cdn.com/flex002/user_avatar/www.rubin.community/melissagraham/32/2549_2.png) [@MelissaGraham](https://www.rubin.community/u/MelissaGraham)\
**Post date:** [August 2, 2023, 11:24pm UTC](https://www.rubin.community/t/missing-semi-major-axis-column-in-mpcorb-table-of-dp0-3/7801/2 "2023-08-02T23:24:00Z")

</div>

Hi Tim, it’s true that the DP0.3 data set did not come with a pre-computed semi-major axis (`a`), and that users have to derive this for themselves using the orbit’s ellipticiy (`e`) and perihelion distance (`q`), as `a` = `q` /(1 - `e`).

Your comment made me realize that the beginner-level introductory DP0.3 tutorials for the Notebook and Portal Aspect don’t cover this conversion. (They’re the only tutorials available so far, and don’t happen to use semi-major axis; but it is covered in a few of the intermediate-level tutorials, coming soon).

Furthermore our [DP0.3 data products definitions document](https://dp0-3.lsst.io/data-products-dp0-3/index.html) doesn’t mention this conversion either – I’ll get it added. Thank you!

---

<div class="post-metadata">

**Author:** ![talister](https://sea2.discourse-cdn.com/flex002/user_avatar/www.rubin.community/talister/32/875_2.png) [@talister](https://www.rubin.community/u/talister)\
**Post date:** [August 2, 2023, 11:54pm UTC](https://www.rubin.community/t/missing-semi-major-axis-column-in-mpcorb-table-of-dp0-3/7801/3 "2023-08-02T23:54:26Z")

</div>

Thank you for the quick response @MelissaGraham I added another column to the table with the formula and that’s working for me now and I could produce the classical _a_ vs _e_ plot:

 ![chart_a-au-e](https://us1.discourse-cdn.com/flex002/uploads/rubin/original/2X/7/7fa81bced95c9da9ade36e3a43f9c3710d5b7fd3.png)  
Now I need to figure out how to overlay the NEO q\<1.3 curve…

The need to compute _a_ was indeed not mentioned in the tutorials or the dp0.3 definitions document as I checked there first in case I was missing something; I also checked the “unpopulated fields” section

---

<div class="post-metadata">

**Author:** ![MelissaGraham](https://sea2.discourse-cdn.com/flex002/user_avatar/www.rubin.community/melissagraham/32/2549_2.png) [@MelissaGraham](https://www.rubin.community/u/MelissaGraham)\
**Post date:** [August 3, 2023, 1:24am UTC](https://www.rubin.community/t/missing-semi-major-axis-column-in-mpcorb-table-of-dp0-3/7801/4 "2023-08-03T01:24:03Z")

</div>

If you’re working in the Portal aspect, I can think of a way to hack plot a q\<1.3 curve – but it’s better to use the Notebook aspect for making plots like that.

> **The hack.**
>
> The hack is to create a new column of values equal to 1 - 1.3/a, and overplot a new trace that is those values vs. a. It will appear as a q\<1.3 curve as a line so long as the trace style is “connected points” and the “`q_LT_1p3`” column is sorted. Note that the “Plot a Function” tool is on the long-term list of functionality to be added.
> 
> ![Screen Shot 2023-08-03 at 6.06.17 PM](https://us1.discourse-cdn.com/flex002/uploads/rubin/original/2X/1/114cb47582e6e3fc81fac563c376af611d06fc64.png)
