KRC for Comets

From krc
(Difference between revisions)
Jump to: navigation, search
(Created page with "== Online Tool == == KRC Davinci Interface Template == == Command Line Example == == KRC Fortran Input File Example == == Common Problems ==")
 
 
(30 intermediate revisions by 2 users not shown)
Line 1: Line 1:
== Online Tool ==
+
Careful!!! For most comets, the rotation period is not known sufficiently well to be reported. In this case, KRC uses a rotation period of 24h by default.
  
== KRC Davinci Interface Template ==
+
Also comets here are treated with the physical properties as for an asteroid!
  
== Command Line Example ==
+
If the rotation period is known by the user, set it with PERIOD.
  
== KRC Fortran Input File Example ==
+
For example:
  
== Common Problems ==
+
out = krc(lat=0.,INERTIA=50.,T_user=220.,LKofT="T",body="Halley",ls=0.)
 +
 
 +
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 ascending 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)
 +
 
 +
This case predict surface temperature for a new body named "IdealJupiterTrojan_like_Comet".
 +
 
 +
[[Image:Made_up_Comet.png|800px]]

Latest revision as of 13:51, 20 August 2020

Careful!!! For most comets, the rotation period is not known sufficiently well to be reported. In this case, KRC uses a rotation period of 24h by default.

Also comets here are treated with the physical properties as for an asteroid!

If the rotation period is known by the user, set it with PERIOD.

For example:

out = krc(lat=0.,INERTIA=50.,T_user=220.,LKofT="T",body="Halley",ls=0.)

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 ascending 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)

This case predict surface temperature for a new body named "IdealJupiterTrojan_like_Comet".

Made up Comet.png

Personal tools