krc()
From krc
Function and Definition
Run KRC through the Davinci Interface: This function takes several inputs to run the KRC model, and output surface temperatures.
Input Arguments
- This Wiki has divided the possible Inputs for Davinci into three parts:
- A Parameter Glossary of the actual input parameters you type into Davinci (i.e. N5).
- A Variable Glossary of the meaning of the input parameters you type into Davinci (i.e. Model Run-Time).
- A Category Glossary of the fields or subjects in which fall the Parameters you type into Davinci (i.e. Time).
- See the Input Parameter Glossary or the KRC planetary surface temperatures: Helplist (PDF) for a full list of all possible input parameters
Output
- Davinci functions as an analytical tool for data. When you run KRC in davinci, the output is contained in a series of structures and elements.
- Elements are simply values (strings, integers, arrays, etc.) generated by davinci
- Structures contain multiple elements (For more see Structures on the DavinciWiki)
When KRC is run in Davinci, the output is contained within a 14 element structure.
- See the Simple Runs Section of the Beginner Tutorial for more.
Examples
The most basic run you can perform in davinci is by only designating a latitude, the only required parameter in davinci.
Start davinci as normal.
>davinci dv>
Then type krc(lat=12) or any latitude you wish. The following will result.
dv> krc(lat=12)
Found Default File: $DV_SCRIPT_FILES/krc_support/porb_defaults/Mars_Mars.porb.hdf
./nfs/software/davinci_install/share/davinci//library/script_files/krc_support/ti_map2ppd_v4.vicar: VICAR bsq image: 720x360x1, 32 bits
./nfs/software/davinci_install/share/davinci//library/script_files/krc_support/albedo_2ppd.vicar: VICAR bsq image: 720x360x1, 32 bits
./nfs/software/davinci_install/share/davinci//library/script_files/krc_support/mola_2ppd.vicar: VICAR bsq image: 720x360x1, 32 bits
Running KRC model with:
Output:
Every ~1.000 Ls (1.9083 JD) for 360 seasons starting at 0.10 Ls
For 96 times per day, every 15 local minutes
Running model WITH NO temperature prediction
because DELJUL < 3.0779
For a Surface With:
29 Layers with uniform material properties and 384 daily model time steps
and T-dependant material properties
and an insulating bottom edge condition
With Observation Parameters:
LATITUDE: 12.00 LONGITUDE: 0.00
ELEVATION: -1.478km^ SLOPE/AZIMUTH: 0.00/0.00deg
ALBEDO: 0.236^ OPACITY: 0.300
HOUR: All LS: All
Upper Material:
INERTIA: 52.00
CONDUCTIVITY: 0.0026*
DENSITY: 1770.00*
SPECIFIC HEAT: 595.55
POROSITY: 0.40
* = derived value unless overwritten by user
^ = from TES 2ppd map unless overwritten by user
& = ajusted for max authorized value
struct, 14 elements
tsurf: 96x1x360 array of double, bsq format [276,480 bytes]
tbol: 96x1x360 array of double, bsq format [276,480 bytes]
tatm: 96x1x360 array of double, bsq format [276,480 bytes]
down_ir: 96x1x360 array of double, bsq format [276,480 bytes]
down_vis: 96x1x360 array of double, bsq format [276,480 bytes]
time: 96x1x1 array of float, bsq format [384 bytes]
ls: 1x1x360 array of double, bsq format [2,880 bytes]
lat: 12.00000000
elev: -1.477849960
layer: struct, 8 elements
thickness: 28x2x1 array of double, bsq format [448 bytes]
center: 28x2x1 array of double, bsq format [448 bytes]
top: 28x2x1 array of double, bsq format [448 bytes]
center_mass: 28x1x1 array of double, bsq format [224 bytes]
mass_burden: 28x1x1 array of double, bsq format [224 bytes]
thermal_scales: 28x1x1 array of double, bsq format [224 bytes]
tmax: 28x1x360 array of double, bsq format [80,640 bytes]
tmin: 28x1x360 array of double, bsq format [80,640 bytes]
anc: struct, 11 elements
krccom: struct, 5 elements...
JDate: 1x1x360 array of double, bsq format [2,880 bytes]
ref_pressure: 1x1x360 array of double, bsq format [2,880 bytes]
taud: 1x1x360 array of double, bsq format [2,880 bytes]
total_frost: 1x1x360 array of double, bsq format [2,880 bytes]
avg_heat_flow: 1x1x360 array of double, bsq format [2,880 bytes]
frost_alb: 1x1x360 array of double, bsq format [2,880 bytes]
frost: 1x1x360 array of double, bsq format [2,880 bytes]
tatm_predict: 1x1x360 array of double, bsq format [2,880 bytes]
delta_t_rms: 1x1x360 array of double, bsq format [2,880 bytes]
converge_days: 1x1x360 array of double, bsq format [2,880 bytes]
version: "v3.2.1"
alb: 0.2363921404
body: "Mars"
If you set the function krc() equal to a variable, you can access the elements contained within the output.
dv> test = krc(lat=12)
To access the elements simply type the variable name followed by a period, then the element you wish to examine. Be careful to keep the dimensions of the element consistent.
To see basic information about the element:
dv> test.tsurf[,,] 96x1x360 array of double, bsq format [276,480 bytes]
Or:
dv> test.tsurf 96x1x360 array of double, bsq format [276,480 bytes]
To see all values for a specific season:
dv> test.tsurf[,,105] 96x1x1 array of double, bsq format [768 bytes] 171.5856935 171.1550320 170.7370510 170.3293537 169.9308028 169.5406723 169.1583461 168.7833589 168.4152575 168.0536924 167.6983010 167.3488043 167.0048976 166.6663525 166.3329117 166.0043847 165.6805461 165.3612358 165.0462598 164.7354822 164.4287279 164.1258806 164.8784105 169.1850639 175.8574576 184.1348010 193.0209648 201.8500108 210.2748614 218.1462596 225.4199940 232.1023999 238.2214758 243.8125454 248.9114546 253.5516366 257.7629989 261.5717402 265.0005054 268.0687358 270.7930092 273.1874139 275.2638321 277.0322298 278.5008521 279.6764380 280.5643349 281.1686471 281.4922865 281.5370772 281.3037515 280.7920112 280.0004782 278.9267210 277.5671609 275.9170552 273.9703619 271.7196781 269.1560659 266.2689434 263.0458944 259.4725498 255.5324779 251.2072449 246.4767807 241.3204027 235.7192896 229.6618763 223.1557308 216.2532801 209.1076064 202.0783448 195.7509294 190.7815122 188.2797335 186.4750573 185.0143727 183.7710532 182.6794830 181.7005257 180.8088599 179.9870621 179.2225887 178.5060773 177.8303253 177.1896796 176.5795809 175.9963471 175.4369088 174.8987397 174.3796661 173.8778849 173.3918072 172.9201071 172.4615814 172.0152291
To see a value for a specific time of day during a specific season:
dv> test.tsurf[13,,105] 167.0048976
- See the Interpreting the Output Section of the Beginner Tutorial
