KRC for Mars
From krc
(Difference between revisions)
(→Command Line Example) |
(→Running in the Command Line) |
||
| (29 intermediate revisions by 2 users not shown) | |||
| Line 1: | Line 1: | ||
| − | = | + | ==Notes== |
| − | + | Basalt is the default material for Mars (Mat1 = "basalt"), and T_user = 220 (temperature at which the inertia is defined). | |
| − | + | ||
| − | + | ||
| − | + | ||
| + | By default, PTOTAL = 545 Pa. | ||
| − | + | =Running in the Command Line= | |
| − | = Command Line | + | |
KRC is run within Davinci. | KRC is run within Davinci. | ||
| Line 14: | Line 11: | ||
> davinci | > davinci | ||
| − | For Mars, it can be run with as little input as a single latitude, e.g. 25 | + | For Mars, it can be run with as little input as a single latitude, e.g. 25°N: |
| − | + | OUT = krc(lat=25.) | |
Model output is saved in structure 'OUT'. The decimal is required for parameters to be initialized as a floating point number. All other model inputs are retrieved from lookup tables (e.g., longitude, surface albedo, thermal inertia etc.). | Model output is saved in structure 'OUT'. The decimal is required for parameters to be initialized as a floating point number. All other model inputs are retrieved from lookup tables (e.g., longitude, surface albedo, thermal inertia etc.). | ||
| Line 24: | Line 21: | ||
tsurf: 96x1x360 array of double, bsq format [276,480 bytes] | tsurf: 96x1x360 array of double, bsq format [276,480 bytes] | ||
| − | Additional fields can be prescribed within the parentheses when calling krc | + | Additional fields can be prescribed within the parentheses when calling krc. E.g., the same latitude but for a longitude of 120° and surface albedo of 0.3: |
| − | + | OUT = krc(lat=25.,lon=120.,ALBEDO=0.3) | |
| − | If a particular season provided as solar longitude (in units of degrees) is desired, the annual dimension (e.g., 360) is removed. E.g., for Ls = 90 | + | If a particular season provided as solar longitude (in units of degrees) is desired, the annual dimension (e.g., 360) is removed. E.g., for Ls = 90° (northern summer solstice): |
| − | + | OUT = krc(lat=25.,ls=90.) | |
| − | + | Alternatively, the model can run for a specific Gregorian Date, (currently ranging from 1990-Jan-01 to 2040-Jan-01, format: ????-Mmm-DD, | |
| − | + | with Mmm:Jan, Feb, Mar, Apr, May, Jun, Jul, Aug, Sep, Oct, Nov, Dec | |
| − | + | ||
| − | + | OUT = krc(lat=12.,GD="2010-Jan-05") | |
| − | + | ||
| − | + | or a specific Julian Date: | |
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| + | OUT = krc(lat=12.,JD=2455201) | ||
| − | + | Note: the possibility to specify the date with GD is only currently available for Mars, the Moon, Bennu, and Europa. | |
| − | + | Now 'OUT.tsurf' has the dimensions of: | |
| + | |||
| + | tsurf: 96x1x1 array of double, bsq format [768 bytes] | ||
| − | + | One can plot the diurnal temperature series against local true solar time (LTST) with: | |
| − | + | plot(OUT.tsurf,xaxis=OUT.time,"25N,Ls=90",w=4 ,color=1) | |
| + | labelxy("LTST","Temperature (K)") | ||
| − | + | [[Image:mars_temp_1.png|border|800px]] | |
| − | + | ||
| − | + | ||
| − | + | Alternatively, one can prescribe local time and output will be provided for a full Mars year, e.g., for a local time of 3 AM: | |
| − | + | OUT = krc(lat=25.,hour=3.) | |
| − | + | ||
| − | + | 'OUT.tsurf' will now be sampled roughly once per degree of solar longitude, e.g.,: | |
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | tsurf: 1x1x360 array of double, bsq format [2,880 bytes] | |
| − | + | ||
| − | + | ||
| − | + | And can be plotted against time (with the addition of 3PM local time as 'OUT_2') with: | |
| − | + | ||
| − | [ | + | plot(OUT.tsurf[12,1,],w=4,color=4,"Dawn",OUT.tsurf[60,1,],w=4,color=1,"Noon") |
| + | labelxy("Solar Longitude","Temperature (K)") | ||
| − | + | [[Image:mars_temp_2.png|border|800px]] | |
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | == Table of Input Parameters == | |
| − | * | + | Other common fields that can be prescribed are included in the table below (*NOTE fields are case sensitive*): |
| − | + | [table of parameters (include example ranges?)] | |
| − | + | ||
| − | + | ||
| − | [ | + | |
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | {| class="wikitable" | |
| − | KRC | + | |+Input Parameters |
| + | |- | ||
| + | |'''Parameter''' | ||
| + | |'''KRC Syntax''' | ||
| + | |'''Range''' | ||
| + | |'''Units''' | ||
| + | |- | ||
| + | |Latitude | ||
| + | |lat | ||
| + | | -90–90 | ||
| + | |- | ||
| + | |Longitude (°E) | ||
| + | |lon | ||
| + | |0–360 | ||
| + | |- | ||
| + | |Albedo | ||
| + | |ALBEDO | ||
| + | |0–1 | ||
| + | |- | ||
| + | |Thermal Inertia | ||
| + | |INERTIA | ||
| + | |20–2000 | ||
| + | |- | ||
| + | |Elevation | ||
| + | |ELEV | ||
| + | |- | ||
| + | |Local True Solar Time | ||
| + | |hour | ||
| + | |0–24 | ||
| + | |- | ||
| + | |Solar Longitude | ||
| + | |ls | ||
| + | |0–360 | ||
| + | |- | ||
| + | |Local True Solar Time | ||
| + | |hour | ||
| + | |} | ||
| − | + | =Deriving Thermal Inertia and the Simulated One Point Mode= | |
| + | In order to derive thermal inertia values for a single temperature measurement one can use the simulated one-point mode in KRC. *''available only for Mars at present'' | ||
| − | + | Please refer to the [[Simulated One Point Mode]] page for more information. | |
Latest revision as of 16:16, 9 April 2020
Contents |
[edit] Notes
Basalt is the default material for Mars (Mat1 = "basalt"), and T_user = 220 (temperature at which the inertia is defined).
By default, PTOTAL = 545 Pa.
[edit] Running in the Command Line
KRC is run within Davinci.
> davinci
For Mars, it can be run with as little input as a single latitude, e.g. 25°N:
OUT = krc(lat=25.)
Model output is saved in structure 'OUT'. The decimal is required for parameters to be initialized as a floating point number. All other model inputs are retrieved from lookup tables (e.g., longitude, surface albedo, thermal inertia etc.).
By default, the output is stored into multidimensional arrays sampled at 96 values per sol and 360 values per Mars year. Hence the structure element surface temperature ('OUT.tsurf') appears as:
tsurf: 96x1x360 array of double, bsq format [276,480 bytes]
Additional fields can be prescribed within the parentheses when calling krc. E.g., the same latitude but for a longitude of 120° and surface albedo of 0.3:
OUT = krc(lat=25.,lon=120.,ALBEDO=0.3)
If a particular season provided as solar longitude (in units of degrees) is desired, the annual dimension (e.g., 360) is removed. E.g., for Ls = 90° (northern summer solstice):
OUT = krc(lat=25.,ls=90.)
Alternatively, the model can run for a specific Gregorian Date, (currently ranging from 1990-Jan-01 to 2040-Jan-01, format: ????-Mmm-DD, with Mmm:Jan, Feb, Mar, Apr, May, Jun, Jul, Aug, Sep, Oct, Nov, Dec
OUT = krc(lat=12.,GD="2010-Jan-05")
or a specific Julian Date:
OUT = krc(lat=12.,JD=2455201)
Note: the possibility to specify the date with GD is only currently available for Mars, the Moon, Bennu, and Europa.
Now 'OUT.tsurf' has the dimensions of:
tsurf: 96x1x1 array of double, bsq format [768 bytes]
One can plot the diurnal temperature series against local true solar time (LTST) with:
plot(OUT.tsurf,xaxis=OUT.time,"25N,Ls=90",w=4 ,color=1)
labelxy("LTST","Temperature (K)")
Alternatively, one can prescribe local time and output will be provided for a full Mars year, e.g., for a local time of 3 AM:
OUT = krc(lat=25.,hour=3.)
'OUT.tsurf' will now be sampled roughly once per degree of solar longitude, e.g.,:
tsurf: 1x1x360 array of double, bsq format [2,880 bytes]
And can be plotted against time (with the addition of 3PM local time as 'OUT_2') with:
plot(OUT.tsurf[12,1,],w=4,color=4,"Dawn",OUT.tsurf[60,1,],w=4,color=1,"Noon")
labelxy("Solar Longitude","Temperature (K)")
[edit] Table of Input Parameters
Other common fields that can be prescribed are included in the table below (*NOTE fields are case sensitive*): [table of parameters (include example ranges?)]
| Parameter | KRC Syntax | Range | Units |
| Latitude | lat | -90–90 | |
| Longitude (°E) | lon | 0–360 | |
| Albedo | ALBEDO | 0–1 | |
| Thermal Inertia | INERTIA | 20–2000 | |
| Elevation | ELEV | ||
| Local True Solar Time | hour | 0–24 | |
| Solar Longitude | ls | 0–360 | |
| Local True Solar Time | hour |
[edit] Deriving Thermal Inertia and the Simulated One Point Mode
In order to derive thermal inertia values for a single temperature measurement one can use the simulated one-point mode in KRC. *available only for Mars at present
Please refer to the Simulated One Point Mode page for more information.
