KRC for Phobos

Revision as of 21:24, 15 January 2019

General advice about using KRC on Phobos

body = "Mars,Phobos"

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

!!!OPACITY = 0.3 !!!! CORRECTION MADE: if PTOTAL == 0., force TAUD = 0.

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