KRC for Mars
Contents |
General advice about using KRC on Mars
When using KRC on Mars the best practice is to:
- Stay away from poles
- Be aware of what season it is
- Be aware of the opacity settings
Command Line Examples
KRC is run within Davinci.
> davinci
For Mars, it can be run with as little input as a single latitude, e.g. 25°N:
dv> 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:
dv> 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):
dv> OUT = krc(lat = 25., ls = 90.)
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:
dv> plot(OUT.tsurf , xaxis = OUT.time , "25N, Ls=90" , w = 4 , color = 1)
dv> 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:
dv> 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:
dv> plot(OUT.tsurf , xaxis = OUT.ls , "LTST=3" , w = 4 , color = 3,
OUT_2.tsurf , xaxis = OUT_2.ls , "LTST=15" , w = 4 , color = 8)
dv> labelxy("Solar Longitude" , "Temperature (K)")
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 | |
| 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 |
KRC Fortran Input File Example
0 0 / KOLD: season to start with; KEEP: continue saving data in same disk file
Version 222 default values. 19 latitudes with mean Mars elevations
ALBEDO EMISS INERTIA COND2 DENS2 PERIOD SPEC_HEAT DENSITY
.25 1.00 200.0 2.77 928.0 1.0275 647. 1600.
CABR AMW [ABRPHA PTOTAL FANON TATM TDEEP SpHeat2
0.11 43.5 -0.00 546.0 .055 200. 180.0 1711.
TAUD DUSTA TAURAT TWILI ARC2 [ARC3 SLOPE SLOAZI
0.3 .90 0.5 0.0 0.5 -0.00 0.0 90.
TFROST CFROST AFROST FEMIS AF1 AF2 FROEXT [FD32
146.0 589944. .65 0.95 0.54 0.0009 50. 0.0
RLAY FLAY CONVF DEPTH DRSET DDT GGT DTMAX
1.2000 .1800 2.0000 0.0 0.0 .0020 0.1 0.1
DJUL DELJUL SOLARDEC DAU LsubS SOLCON GRAV AtmCp
-1222.69 17.174822 00.0 1.465 .0 1368. 3.727 735.9
ConUp0 ConUp1 ConUp2 ConUp3 ConLo0 ConLo1 ConLo2 ConLo3
0.038640 -0.002145 0.002347 -0.000750 2.766722 -1.298966 0.629224 -0.527291
SphUp0 SphUp1 SphUp2 SphUp3 SphLo0 SphLo1 SphLo2 SphLo3
646.6275 246.6678 -49.8216 7.9520 1710.648 721.8740 57.44873 24.37532
N1 N2 N3 N4 N5 N24 IB IC
20 384 15 19 120 48 0 9
NRSET NMHA NRUN JDISK IDOWN FlxP14 FlxP15 KPREF
3 24 0 81 0 45 65 1
K4OUT JBARE Notif IDISK2 end
52 0 20 0 0
LP1 LP2 LP3 LP4 LP5 LP6 LPGLOB LVFA LVFT LkofT
F T F F F F F F F T
LPORB LKEY LSC spare LOCAL Prt76 LPTAVE Prt78 Prt79 L_ONE
T F F F F T F T F F
Latitudes: in 10F7.2 _____7 _____7 _____7 _____7 _____7 _____7 _____7
-87.50 -80.00 -70.00 -60.00 -50.00 -40.00 -30.00 -20.00 -10.00 0.00
10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 87.50 -0.00
_____7 _____7 _____7 Elevations: in 10F7.2 ____7 _____7 _____7 _____7
3.51 2.01 1.39 1.22 0.38 0.48 1.17 1.67 1.26 0.17
-0.94 -1.28 -1.99 -2.51 -3.52 -4.08 -4.51 -4.38 -2.57 -0.00
2013 Jul 24 11:28:09=RUNTIME. IPLAN AND TC= 104.0 0.10000 Mars:Mars
104.0000 0.1000000 0.8644665 0.3226901E-01 -1.281586
0.9340198E-01 1.523712 0.4090926 0.000000 0.9229373
5.544402 0.000000 0.000000 686.9929 3397.977
24.62296 0.000000 -1.240317 0.000000 0.000000
0.000000 0.3244965 0.8559126 0.4026359 -0.9458869
0.2936298 0.1381285 0.000000 -0.4256703 0.9048783
8 0 0 'master222.t52' / Disk file name for Run 1
0/
3 10 1 'LkofT' / Temperature-dependant conductivity
0/
0/ ======================= end of run
If LkofT set to T, then
Upper material: weakly cemented particulates:
Grain: k: BasicRocks_Zoth88 2:4
Cp: Chlorite_Bert07_Fe=0.89 6:4.89
Cement: k: Limestone Zoth88 2:1
Cp: Sphene, which has relatively strong T dependence 5:0
Cement fraction 1.e-8
Yields c_0 of 0.050087, This adjusted to 0.038640 to agree with I=200 at 220 K
lower material: H2O Ice
k: koftop: 48
fit to A+B/T fit 2.766722 -1.298966 0.629224 -0.527291 <k H2O:ice3sources
Cp: koftop: @ 49,491,5,33 yields:
3.95779 > 1710.648 721.8740 57.44873 24.37532 <SpH H2O:Ice_3sources
