KRC for Phobos

From krc
(Difference between revisions)
Jump to: navigation, search
Line 1: Line 1:
 
== General advice about using KRC on Phobos ==
 
== General advice about using KRC on Phobos ==
  
−
#
 
−
###PHOBOS#################################################################################
 
−
#
 
 
body  = "Mars,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")
 
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.
+
 
 +
 
 +
!!!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.)
 
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")
 
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)            #
 
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")  
 
Eclipse      = "T"                                                                      #forces an eclipse (Default = "F")  
Line 18: Line 25:
 
Date          = 5000.
 
Date          = 5000.
 
test          = krc_eclipse(Eclipser=Eclipser,Ecl_Cent_Hr=Ecl_Cent_Hr,Bias=Bias,Eclipse_Style=Eclipse_Style,Date=Date,body=body)
 
test          = krc_eclipse(Eclipser=Eclipser,Ecl_Cent_Hr=Ecl_Cent_Hr,Bias=Bias,Eclipse_Style=Eclipse_Style,Date=Date,body=body)
 +
 +
 
#  14 1 5.2026 71492.  0.6711D6  1560.8  3.551 0.01 6000. 12.  2 77 77 / Europa        #HHK EXAMPLE
 
#  14 1 5.2026 71492.  0.6711D6  1560.8  3.551 0.01 6000. 12.  2 77 77 / Europa        #HHK EXAMPLE
 
#  "14 1 5.20  69911.  670900.00 3121.60  3.55  0.00 0.00  12.00 2      / Eclipse"      #DV INTERFACE
 
#  "14 1 5.20  69911.  670900.00 3121.60  3.55  0.00 0.00  12.00 2      / Eclipse"      #DV INTERFACE
Line 26: Line 35:
 
                                                 #Eclipse Bias        #N/A
 
                                                 #Eclipse Bias        #N/A
 
                                                             #Eclipse central hour (Should be calculated from subjovian latitude by default)
 
                                                             #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")
+
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
 
plot(OUT_7.tsurf[,1,1],xaxis=OUT_7.time)                                                  #Daily Eclipse
  
Line 39: Line 51:
 
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_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
 
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
 
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

Revision as of 21:23, 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

Retrieved from "https://krc.mars.asu.edu/index.php?title=KRC_for_Phobos&oldid=2365"
Personal tools