KRC for Exoplanets
| Line 9: | Line 9: | ||
Where | Where | ||
| − | :e: Eccentricity | + | :e: Eccentricity (Default=0) |
| − | :a: | + | :a: Semi-Major Axis in AU (Default=1) |
| − | :i: | + | :i: Inclination of mean orbit to ecliptic in degrees (Default=0) |
| − | :node: | + | :node: Longitude of the asceding node in degrees (Default=0) |
| − | :per: | + | :per: Argument of perihelion in degrees (Default=0) |
| − | :m: | + | :m: Mean Anomaly at epoch in degrees (Default=0) |
| − | :rot_per: | + | :rot_per: siderial rotation period in hours (Default=23.9345) |
| − | :polera: | + | :polera: Right Ascension of the pole in degrees (Default=0) |
| − | :merid: | + | :merid: Prime meridian at epoch in degrees (Default=0) |
| − | :poledec: | + | :poledec: Declination of the pole in degrees (Default=0) |
| − | :period: | + | :period: Sideral orbital period in Earth days (Default=365.256) |
| − | :name: | + | :name: Body name (default ="None") |
| + | |||
| + | |||
| + | printf("\nGenerate a generic PORB structure for use with the davinci porb function\n") | ||
| + | printf("This is formatted in the minor body style for porb and permits the following values described below\n\n") | ||
| + | printf("\tname = body name (Default=\"None\"\n") | ||
| + | printf("\tepoch = epoch in Julian Date (Default = 2451545.0, year 2000)\n") | ||
| + | printf("\ta = Semi-Major Axis in AU (Default=1)\n") | ||
| + | printf("\te = Eccentricity (Default=0)\n") | ||
| + | printf("\ti = Inclination of mean orbit to ecliptic in degrees (Default=0)\n") | ||
| + | printf("\tnode = Longitude of the asceding node in degrees (Default=0)\n") | ||
| + | printf("\tperi = Argument of perihelion in degrees (Default=0)\n") | ||
| + | printf("\tm = Mean Anomoly at epoch in degrees (Default=0)\n") | ||
| + | printf("\tpolera = Right Ascention of the pole in degrees (Default=0)\n") | ||
| + | printf("\tpoledec = Declination of the pole in degrees (Default=0)\n") | ||
| + | printf("\tmerid = prime meridian at epoch in degrees (Default=0\n") | ||
| + | printf("\trot_per = siderial rotation period in hours (Default=23.9345)\n") | ||
| + | printf("\tperiod = siderial orbital period in days (Default=365.256)\n") | ||
Revision as of 20:33, 15 January 2019
A user can define the orbital parameters of any body rotating around any star. PORB generates the rotation matrix accordingly.
Solar System-like
First, run porb to define the orbital parameters of the body in a structure:
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=poledec, merid=7.7, period=4332.589, name="IdealJupiterTrojan")
Where
- 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")
printf("\nGenerate a generic PORB structure for use with the davinci porb function\n")
printf("This is formatted in the minor body style for porb and permits the following values described below\n\n")
printf("\tname = body name (Default=\"None\"\n")
printf("\tepoch = epoch in Julian Date (Default = 2451545.0, year 2000)\n")
printf("\ta = Semi-Major Axis in AU (Default=1)\n")
printf("\te = Eccentricity (Default=0)\n")
printf("\ti = Inclination of mean orbit to ecliptic in degrees (Default=0)\n")
printf("\tnode = Longitude of the asceding node in degrees (Default=0)\n")
printf("\tperi = Argument of perihelion in degrees (Default=0)\n")
printf("\tm = Mean Anomoly at epoch in degrees (Default=0)\n")
printf("\tpolera = Right Ascention of the pole in degrees (Default=0)\n")
printf("\tpoledec = Declination of the pole in degrees (Default=0)\n")
printf("\tmerid = prime meridian at epoch in degrees (Default=0\n")
printf("\trot_per = siderial rotation period in hours (Default=23.9345)\n")
printf("\tperiod = siderial orbital period in days (Default=365.256)\n")
poledec=0; 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=poledec, merid=7.7, period=4332.589, name="IdealJupiterTrojan")
out=krc(lat=0.,lon=0,INERTIA=20.,ALBEDO=.05,LKofT="F",body=tmp,bodytype="minor") #
PLOTS AND EXPLANATION NEEDED HERE
