KRC for Europa
(→Known Issues) |
|||
Line 99: | Line 99: | ||
:Date: ??? | :Date: ??? | ||
− | dv> OUT = krc(lat = 0., INERTIA = 045., N1 = 32, body = | + | dv> OUT = krc(lat = 0., INERTIA = 045., N1 = 32, body = "Jupiter,Europa", bodytype = "minor", N24 = 96, Eclipse = "T", Eclipser="Jupiter_Jupiter", Ecl_Cent_Hr = 12., Bias = 0., Eclipse_Style = 1., Date = 5000.) |
[[Image:Europa_F03.png|500px]] Example of Europa diurnal curve with Jupiter Eclipse centered at Noon | [[Image:Europa_F03.png|500px]] Example of Europa diurnal curve with Jupiter Eclipse centered at Noon |
Revision as of 21:59, 8 January 2019
Notes
Some parameters are set to realistic values by default, including:
- Mat1 = "H2O" which assignes water ice Cp, Density, and Conductivity properties
- Mat2 = "H2O" which assignes water ice Cp, Density, and Conductivity properties
- ALBEDO = 0.67
- PTOTAL = 0 which eliminates the atmosphere (ELEV, TAU, etc. unused)
When running the DaVinci interface, default Mars atmospheric values might be printed on the screen but not effectively used (PTOTAL = 0.)
Command Line Examples
Basic surface temperatures for Europa:
dv> OUT = krc(lat = 0, INERTIA = 45., body="Jupiter,Europa", bodytype = "minor", ALBEDO = .55, LKofT = "F")
Example of simple diurnal temperature curve
labelxy("LTST","Temperature [K]") plot(OUT.tsurf[,1,1],xaxis=OUT.time,"45 Kieffer, No Flux",w=2,color=2)
To include the visible and thermal infrared flux from Jupiter, two approaches are possible:
1: Use default built-in parameters, set PFlux = "T", and Lon_Hr
- Lon_Hr [0-24] is the surface longitude relative to the sub-Jupiter point, expressed in hours.
- When Lon_Hr < 6, the surface point is on the Antijovian hemisphere => No flux contributed.
- When 6 < Lon_Hr < 18, the surface point is on the subjoin hemisphere => Flux is contributed (max at Lon_Hr = 12.).
- When Lon_Hr > 18, the surface point is on the Antijovian hemisphere => No flux contributed.
dv> OUT = krc(lat = 0., INERTIA = 70., body = "Jupiter,Europa", bodytype = "minor", ALBEDO = 0.55, PFlux = "T", Lon_Hr = 12., LKofT = "F")
2: Provide values for all the necessary input parameters, and set PFlux = "T" (Default provided for common bodies)
- BT_Avg : Average Brightness Temperature [K]
- BT_Min : Min Brightness Temperature, if diurnal cycle [K]
- BT_Max : Max Brightness Temperature [K]
- Dis_AU : Distance from Sun in AU
- Geom_alb : Geometric Albedo [1]
- Mut_Period : Mutual Period [?]
- Orb_Radius : Orbiting Radius [km]
- Radius : Radius of the Obiting body [km]
- Lon_Hr : Longitude Hour of the surface point (see above)
dv> OUT = krc(lat = 0., INERTIA = 70., body = "Jupiter,Europa", bodytype = "minor", ALBEDO = 0.55, PFlux = "T", BT_Avg = 127., BT_Min = 127., BT_Max = 127., Dis_AU = 5.203, Geom_alb = 0.52, Mut_Period = 3.55, Orb_Radius = 670900, Radius = 670900, Lon_Hr = 12., LKofT = "F")
3: Provide Visible and IR flux tables vs. LTST, and set PFlux = "T"
- The interface fits sin functions through the table values, and extracts parameters required by KRC (amplitude, phase, etc.). For Europa, this is not the preferred option.
- IR: A 2 x n x 1 array with IR flux (1st col.) vs. LTST (2nd col.)
- Vis: A 2 x n x 1 array with Vis flux (1st col.) vs. LTST (2nd col.)
dv> OUT = krc(lat = 0., INERTIA = 70., body = "Jupiter,Europa", bodytype = "minor", ALBEDO = 0.55, PFlux = "T", Lon_Hr = 12., LKofT = "F")
Comparison between simple diurnal temperature curves (with vs. without Jupiter flux)
dv > OUT_2 = krc(lat = 0, INERTIA = 70., body = "Jupiter,Europa", bodytype = "minor", ALBEDO = .55, LKofT="F") dv > OUT_3 = krc(lat = 0., INERTIA = 70., body = "Jupiter,Europa", bodytype = "minor", ALBEDO = 0.55, PFlux = "T", Lon_Hr = 12., LKofT = "F") dv > labelxy("LTST","Temperature [K]") dv > plot(OUT_2.tsurf[,1,1], xaxis = OUT_2.time, "70 Kieffer, No Flux", w = 2, color = 2, OUT_3.tsurf[,1,1], xaxis = OUT_3.time, "70 Kieffer, With Flux", w = 2, color = 3)
To include an Eclipse by Jupiter, set Eclipse = "T" and specify the following parameters:
- Eclipse = "T" forces an eclipse (Default = "F")
- body = "Jupiter,Europa"
- Eclipser = "Jupiter_Jupiter" Eclipser name, for Example "Jupiter_Jupiter" or "Mars_Phobos"
- Eclipse_Style = 1 because eclipses are assumed to occur daily; see dedicated Eclipse Section [Build Link Here]
- Ecl_Cent_Hr: Eclipse central hour [subjovian point => =12.; Antijovian point => =0.]
- Bias = 0.0: Eclipse Bias (0 => perfect alignement; 1 => partial eclipse); see dedicated Eclipse Section [Build Link Here]
- Date: ???
dv> OUT = krc(lat = 0., INERTIA = 045., N1 = 32, body = "Jupiter,Europa", bodytype = "minor", N24 = 96, Eclipse = "T", Eclipser="Jupiter_Jupiter", Ecl_Cent_Hr = 12., Bias = 0., Eclipse_Style = 1., Date = 5000.)
Example of Europa diurnal curve with Jupiter Eclipse centered at Noon
labelxy("LTST","Temperature [K]","45 Kieffer, Eclipse centered at Noon") plot(OUT.tsurf[,1,1],xaxis=OUT.time,"Noon Eclipse")
Known Issues
1: The default Jupiter_Europa.porb.hdf file has an issue:
EUROPA = load($DV_SCRIPT_FILES + "/krc_support/porb_defaults/Jupiter_Europa.porb.hdf") EUROPA.planet_flux.Radius = 3121.60
==> 3121.60 IS THE DIAMETER, NOT THE RADIUS #REGENERATE .HDF WIHT RADIUS VALUE (ONLY EUROPA SEEMS IMPACTED)
2: Simulated One Point Mode Never Tested