KRC for Exoplanets
From krc
(Difference between revisions)
| Line 1: | Line 1: | ||
A user can define the orbital parameters of any body rotating around any star. PORB generates the rotation matrix accordingly. | A user can define the orbital parameters of any body rotating around any star. PORB generates the rotation matrix accordingly. | ||
| − | + | First, run PORB to create a structure that can be ingested to generate the orbital rotation matrix. The following arguments are needed: | |
| − | + | ||
| − | First, run PORB the orbital | + | |
:e: Eccentricity (Default=0) | :e: Eccentricity (Default=0) | ||
| Line 29: | Line 27: | ||
:name: Body name (default ="None") | :name: Body name (default ="None") | ||
| − | + | Then run generic_porb() with these arguments, for example: | |
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| + | tmp=generic_porb(e = 0, a = 1., i = 1.30439695, node = 100.47390909, peri = 293.923, m = 79.668, rot_per = 200., polera = 273.85, poledec = 0., merid = 7.7, period = 4332.589, name = "IdealJupiterTrojan") | ||
| + | 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 above assumes the Sun as the central star. Any star can be defined by setting SOLCON (solar constant defined as the total stellar irradiance at 1AU in W/m^2. In the solar system, SOLCON ~ 1,361 W/m2. | |
| − | out=krc(lat=0.,lon=0,INERTIA=20.,ALBEDO=.05,LKofT="F",body=tmp,bodytype="minor") | + | out = krc( lat = 0., lon = 0, INERTIA = 20., ALBEDO =.05, LKofT = "F", body = tmp, bodytype="minor", SOLCON = 2722.) |
| − | + | The case above assumes a start outputting twice as much energy as the Sun. | |
Revision as of 20:51, 15 January 2019
A user can define the orbital parameters of any body rotating around any star. PORB generates the rotation matrix accordingly.
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.30439695, node = 100.47390909, peri = 293.923, m = 79.668, rot_per = 200., polera = 273.85, poledec = 0., merid = 7.7, period = 4332.589, name = "IdealJupiterTrojan")
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 above assumes the Sun as the central star. Any star can be defined by setting SOLCON (solar constant defined as the total stellar irradiance at 1AU in W/m^2. In the solar system, SOLCON ~ 1,361 W/m2.
out = krc( lat = 0., lon = 0, INERTIA = 20., ALBEDO =.05, LKofT = "F", body = tmp, bodytype="minor", SOLCON = 2722.)
The case above assumes a start outputting twice as much energy as the Sun.
