krc one

From DavinciWiki

Jump to: navigation, search

Contents: Description, Arguments, Usage, Examples, Sub-Functions, Related Functions, Source

Supported Under Version: >1.00

Current .dvrc Version: 1.01

Description

NOTE: KRC is not included with standard davinci distributions.

Determines a surface thermal inertia from an input temperature using the KRC one point thermal model.

Output can be a single value or a structure containing ancillary information.

Arguments and Return Values

Returns a single value (ti) or a structure containing:

ti_scale: 1x100x1 array thermal inertia values used in the model
tk: 1x100x1 array       surface kinetic temperatures corresponding to ti_scale
tbol: 1x100x1 array     bolometric temperatures corresponding to ti_scale
tes_ti: 170.8125000     thermal inertia for input lat/lon from 2ppd TES map
alb: 0.1864313632       albedo used in the model
elev: -1.410959959      elevation used in the model
ti: 411.8508911         thermal inertia interpolated using ti_scale and tk or tbol

Usage

Syntax: ti = krc_one(temp [, lat = FLOAT] [, lon = FLOAT] [, elev = FLOAT] [, alb = FLOAT] [, tau = FLOAT] [, slope = FLOAT] [, azimuth = FLOAT] [, hour = FLOAT] [, bol = BOOLEAN] [, ls = FLOAT] [, struct = BOOLEAN])

temp    = surface kinetic or bolometric temperature
lat     = latitude (default=0)
lon     = east longitude (default=0)
elev    = elevation km (default from lat,lon 2ppd MOLA)
alb     = albedo (default from lat,lon 2ppd TES)
tau     = dust visible opacity (default=0.30)
slope   = surface slope in degrees from horizontal (default=0)
azim    = slope azimuth in degrees east from north (default=0)
hour    = hour of day for data (default=2)
bol     = use bolometric rather than kinetic temperature (default=0)
ls      = Lsubs (default=0)
struct  = Output ancillary info in a structure (default=0)

Examples

dv> a=krc_one(280,hour=14,lat=0,ls=300,struct=1)
./usr/share/davinci/krc/krc_bin/ti_map2ppd_v3.vicar: VICAR bsq image: 720x360x1, 32 bits
./usr/share/davinci/krc/krc_bin/albedo_2ppd.vicar: VICAR bsq image: 720x360x1, 32 bits
./usr/share/davinci/krc/krc_bin/mola_2ppd.vicar: VICAR bsq image: 720x360x1, 32 bits

TES Inertia: 170.812500
TES Albedo: 0.186431
Elevation (km): -1.410960 

 ?* Input file name or / for default =
 krc.inp                                                                         
  ?* Print file name or / for default =
 krc.prt                                                                         
  TDISK:2, WRITE, BUT NO FILE OPEN
  TDISK:2, WRITE, BUT NO FILE OPEN
  TDISK:2, WRITE, BUT NO FILE OPEN
  TDISK:2, WRITE, BUT NO FILE OPEN
  TDISK:2, WRITE, BUT NO FILE OPEN...



...Apparent file size: 12x100x1
Read ASCII file: 12x100x1
struct, 7 elements
    ti_scale: 1x100x1 array of float, bsq format [400 bytes]
    tk: 1x100x1 array of float, bsq format [400 bytes]
    tbol: 1x100x1 array of float, bsq format [400 bytes]
    tes_ti: 170.8125000    
    alb: 0.1864313632    
    elev: -1.410959959    
    ti: 411.8508911

DavinciWiki Mini-Nav Bar

Contents