KRC for Phobos

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?

OUT_6 = krc(lat=12.,INERTIA=50.,Mat1="basalt",Mat2="H2O",TI_CO=1400.,body=body,INERTIA2=200.,thick=0.1,LKofT="T",PFlux="T",Lon_Hr=8.)

labelxy("LTST","Temperature [K]","12N, 70 Kieffer, With Flux")

plot(OUT_5.tsurf[,1,1],"No Flux",OUT_6.tsurf[,1,1],"Flux",Xaxis=OUT_5.time,"XXX",w=2,color=2) #

Eclipse = "T" #forces an eclipse (Default = "F")

Eclipser = "Mars_Mars" #Eclipser name, for Example "Jupiter_Jupiter" or "Mars_Phobos"

Ecl_Cent_Hr = 12. #Eclipse central hour [subjovian point => =12.; Antijovian point => =0.]

Bias = 0.0 #Eclipse Bias (0 => perfect alignement; 1 => partial eclipse)

Eclipse_Style = 1 #==1 > daily ; ==2 > rare

Date = 5000.

test = krc_eclipse(Eclipser=Eclipser,Ecl_Cent_Hr=Ecl_Cent_Hr,Bias=Bias,Eclipse_Style=Eclipse_Style,Date=Date,body=body)

1. 14 1 5.2026 71492. 0.6711D6 1560.8 3.551 0.01 6000. 12. 2 77 77 / Europa #HHK EXAMPLE
2. "14 1 5.20 69911. 670900.00 3121.60 3.55 0.00 0.00 12.00 2 / Eclipse" #DV INTERFACE

      #Eclipse Style   #Mutual center-of-mass orbit radius       #debug code
        #Distance to SUN AU      #EB(Europa) Radius, km XXX should be 1560.8
               #OB (Jupiter) Radius in km XXX Should be 71492.     #N/A
                                          #Mutual solar synodic period, days
                                                #Eclipse Bias         #N/A
                                                           #Eclipse central hour (Should be calculated from subjovian latitude by default)

OUT_7 =

krc(lat=12.,INERTIA=50.,Mat1="basalt",Mat2="H2O",TI_CO=1400.,body=body,INERTIA2=200.,thick=0.1,LKofT="T",PFlux="T",Lon_Hr=8.,TUN8=0,N24=96*4,Eclipse=Eclipse,Eclipser=Eclipser,Ecl_Cent_Hr=Ecl_Cent_Hr,Bias=Bias,Eclipse_Style=Eclipse_Style,Date=Date,body=body,bodytype="minor")

plot(OUT_7.tsurf[,1,1],xaxis=OUT_7.time) #Daily Eclipse

body = "Mars"

Eclipse = "T" #forces an eclipse (Default = "F")

Eclipser = "Mars_Phobos" #Eclipser name, for Example "Jupiter_Jupiter" or "Mars_Phobos"

Ecl_Cent_Hr = 12. #Eclipse central hour [subjovian point => =12.; Antijovian point => =0.]

Bias = 0.0 #Eclipse Bias (0 => perfect alignement; 1 => partial eclipse)

Eclipse_Style = 6 #==1 > daily ; ==2+ > rare

Date = 6000.

test = krc_eclipse(Eclipser=Eclipser,body=body,Ecl_Cent_Hr=Ecl_Cent_Hr,Bias=Bias,Eclipse_Style=Eclipse_Style,Date=Date)

OUT_8 = krc(TUN8="F",Date=Date,ls=359,lat=0.,INERTIA=50.,Mat1="basalt",body=body,LKofT="F",Eclipse=Eclipse,Eclipser=Eclipser,Ecl_Cent_Hr=Ecl_Cent_Hr,Bias=Bias,Eclipse_Style=Eclipse_Style,N24=96*100)

OUT = OUT_8.tsurf

plot(OUT[dim(OUT)[1]/2-16:dim(OUT)[1]/2+16,,1],xaxis=OUT_8.time[dim(OUT)[1]/2-16:dim(OUT)[1]/2+16,,1],ignore=-999.) #Daily Eclipse