KRC for Exoplanets

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