krc()

From krc
Revision as of 18:18, 30 July 2017 by Aaron (Talk | contribs)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

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 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
Personal tools