KRC for Comets
From krc
(Difference between revisions)
(Blanked the page) |
|||
Line 1: | Line 1: | ||
+ | Orbital parameters have already been processed for some comets: | ||
+ | 19P-Borrelly | ||
+ | |||
+ | 22P-Kopff | ||
+ | |||
+ | Churyumov-Gerasimenko | ||
+ | |||
+ | Giacobini-Zinner | ||
+ | |||
+ | Halley | ||
+ | |||
+ | Tempel_1 | ||
+ | |||
+ | Tempel_2 | ||
+ | |||
+ | Wild_2 | ||
+ | |||
+ | For these bodies, directly run KRC using the body's name as follow: | ||
+ | |||
+ | : body = "minor,BODY_NAME" | ||
+ | |||
+ | For example: | ||
+ | |||
+ | out = krc(lat = 0., INERTIA = 50., T_user = 220., LkofT = "T", body = "minor_Wild_2", ls = 0.) | ||
+ | |||
+ | HERE PLOT SHOWING RESULTS | ||
+ | |||
+ | |||
+ | |||
+ | For other bodies -not in the list above-, the user needs to generate the rotational matrix. First, run PORB to create a structure that can be ingested to generate the orbital rotation matrix. The following arguments are needed: | ||
+ | |||
+ | :e: Eccentricity (Default=0) | ||
+ | |||
+ | :a: Semi-Major Axis in AU (Default=1) | ||
+ | |||
+ | :i: Inclination of mean orbit to ecliptic in degrees (Default=0) | ||
+ | |||
+ | :node: Longitude of the asceding node in degrees (Default=0) | ||
+ | |||
+ | :per: Argument of perihelion in degrees (Default=0) | ||
+ | |||
+ | :m: Mean Anomaly at epoch in degrees (Default=0) | ||
+ | |||
+ | :rot_per: siderial rotation period in hours (Default=23.9345) | ||
+ | |||
+ | :polera: Right Ascension of the pole in degrees (Default=0) | ||
+ | |||
+ | :merid: Prime meridian at epoch in degrees (Default=0) | ||
+ | |||
+ | :poledec: Declination of the pole in degrees (Default=0) | ||
+ | |||
+ | :period: Sideral orbital period in Earth days (Default=365.256) | ||
+ | |||
+ | :name: Body name (default ="None") | ||
+ | |||
+ | |||
+ | Then run generic_porb() with these arguments, for example: | ||
+ | |||
+ | tmp=generic_porb(e = 0, a = 1., i = 1.3, node = 100.4, peri = 293.9, m = 79.6, rot_per = 200., polera = 273.8, poledec = 0., merid = 7.7, period = 4332.5, name = "IdealJupiterTrojan_like_Comet") | ||
+ | |||
+ | |||
+ | Then run KRC using this new body as the "body": | ||
+ | |||
+ | out = krc( lat = 0., lon = 0, INERTIA = 20., ALBEDO =.05, LKofT = "F", body = tmp, bodytype = "minor") | ||
+ | |||
+ | This case predict surface temperature for a new body named "IdealJupiterTrojan". | ||
+ | |||
+ | HERE PLOT SHOWING RESULTS |
Revision as of 21:12, 15 January 2019
Orbital parameters have already been processed for some comets:
19P-Borrelly
22P-Kopff
Churyumov-Gerasimenko
Giacobini-Zinner
Halley
Tempel_1
Tempel_2
Wild_2
For these bodies, directly run KRC using the body's name as follow:
- body = "minor,BODY_NAME"
For example:
out = krc(lat = 0., INERTIA = 50., T_user = 220., LkofT = "T", body = "minor_Wild_2", ls = 0.)
HERE PLOT SHOWING RESULTS
For other bodies -not in the list above-, the user needs to generate the rotational matrix. First, run PORB to create a structure that can be ingested to generate the orbital rotation matrix. The following arguments are needed:
- e: Eccentricity (Default=0)
- a: Semi-Major Axis in AU (Default=1)
- i: Inclination of mean orbit to ecliptic in degrees (Default=0)
- node: Longitude of the asceding node in degrees (Default=0)
- per: Argument of perihelion in degrees (Default=0)
- m: Mean Anomaly at epoch in degrees (Default=0)
- rot_per: siderial rotation period in hours (Default=23.9345)
- polera: Right Ascension of the pole in degrees (Default=0)
- merid: Prime meridian at epoch in degrees (Default=0)
- poledec: Declination of the pole in degrees (Default=0)
- period: Sideral orbital period in Earth days (Default=365.256)
- name: Body name (default ="None")
Then run generic_porb() with these arguments, for example:
tmp=generic_porb(e = 0, a = 1., i = 1.3, node = 100.4, peri = 293.9, m = 79.6, rot_per = 200., polera = 273.8, poledec = 0., merid = 7.7, period = 4332.5, name = "IdealJupiterTrojan_like_Comet")
Then run KRC using this new body as the "body":
out = krc( lat = 0., lon = 0, INERTIA = 20., ALBEDO =.05, LKofT = "F", body = tmp, bodytype = "minor")
This case predict surface temperature for a new body named "IdealJupiterTrojan".
HERE PLOT SHOWING RESULTS