KRC for Phobos
Line 31: | Line 31: | ||
HERE Figure T vs LTST with and Without Mars shine on Phobos | HERE Figure T vs LTST with and Without Mars shine on Phobos | ||
+ | |||
====Mars Eclipses==== | ====Mars Eclipses==== |
Revision as of 10:00, 16 January 2019
Contents |
Basic Case
Phobos models require to specify body = "Mars,Phobos"
out = krc(lat = 12. ,INERTIA = 50. ,body = "Mars,Phobos", LKofT = "T")
Note that the screen may print:
OPACITY = 0.3 as well as some other atmospheric properties values. Because PTOTAL = 0., KRC runs adequately - without an atmosphere -, but the interface prints default (Mars) atmosphere values. This will be fixed in future iterations.
Fluxes from Mars
Generally, Phobos modeling requires the addition of the visible and IR fluxes from Mars:
out = krc(lat = 12. ,INERTIA = 50. ,body = "Mars,Phobos", LKofT = "T", PFlux = "T",Lon_Hr = 8.)
...if standard flux values are acceptable. Alternatively, the user can provide Visible and IR fluxes versus local time (krc_planetary_flux_table() generates the input parameters):
HERE Figure IR and VIS vs. local time
In this case, the flux from Mars on Phobos can be formalized like this:
out = krc(lat = 12. ,INERTIA = 50. ,body = "Mars,Phobos", LKofT = "T", PFlux = "T",Lon_Hr = 8., IR = IR, Vis = Vis)
See below a comparison of the predicted surface temperature with and without Mars IR and visible shine on Phobos:
HERE Figure T vs LTST with and Without Mars shine on Phobos
Mars Eclipses
To include the eclipse by Mars using standard Phobos/Mars orbital properties:
- Eclipse = "T" #forces an eclipse (Default = "F")
- Eclipser = "Mars_Mars" #Eclipser name
- Ecl_Cent_Hr = 12. #Eclipse central hour
- Bias = 0.0 #Eclipse Bias (0 => perfect alignement; 1 => partial eclipse)
- Eclipse_Style = 1
- Date = 5000.
out = krc( lat = 12. ,INERTIA = 50. ,body = "Mars,Phobos", LKofT = "T", PFlux = "T",Lon_Hr = 8., IR = IR, Vis = Vis, Eclipse = Eclipse, Eclipser = Eclipser, Ecl_Cent_Hr = Ecl_Cent_Hr, Bias = Bias, Eclipse_Style = Eclipse_Style,Date = Date, bodytype = "minor")
!!! Do we need bodytype = "minor" for Phobos?
When working with Phobos eclipses, it may be beneficial to increase the time resolution of the output by increasing N24
out = krc( lat = 12. ,INERTIA = 50. ,body = "Mars,Phobos", LKofT = "T", PFlux = "T",Lon_Hr = 8., IR = IR, Vis = Vis, Eclipse = Eclipse, Eclipser = Eclipser, Ecl_Cent_Hr = Ecl_Cent_Hr, Bias = Bias, Eclipse_Style = Eclipse_Style,Date = Date, bodytype = "minor", N24 = 96 * 100)