surfcalc
|
Contents: Citation, Description, Arguments, Usage, Examples, Images, Sub-Functions, Related Functions, Source Supported Under Version: >1.00
Citation Bandfield, J. L., and C. S. Edwards (2008), Derivation of martian surface slope characteristics from directional thermal infrared radiometry, Icarus, 193(1), 139-157.
Description Creates a surface brightness temperature given surface and atmospheric parameters. This is used to estimate surface kinetic temperature rather than simply relying on the highest measured brightness temperature.
Arguments and Return Values Arguments: Surface emissivity, atmospheric opacity, pressure-temperature profile, surface pressure, measured radiance, emission angle. Return Value: Temperature as a function of wavelength corrected for both atmospheric opacity and surface emissivity.
Usage Syntax: surfcalc(surf_em, tau, pt, p_surf, radiance, eangle,[wn]) surf_em: surface emissivity (m x n x 1 array) tau: atmospheric column opacity (m x n x 1 array) pt: pressure-temperature profile (m x 21 x 1 array using pressure levels from pt) p_surf: surface pressure (m x 1 x 1 array, in mbar) radiance: measured radiance (m x n x 1 array, in W cm-2 sr-1 cm) eangle: emission angle of simulated observation (m x 1 x 1 array, degrees) wn: wavenumber of single point calculation (optional value, used only if surf_em and tau are single point arrays) struct: return a structure with individual transmitted and radiance components - 0 (default) returns radiance only, 1 returns individual components surfcalc only works on 1, 143, or 286 point spectra
Examples dv> surf_em 1x143x1 array of float, bsq format [572 bytes] dv> tau 1x143x1 array of float, bsq format [572 bytes] dv> pt 1x21x1 array of double, bsq format [168 bytes] 0.000000000 0.000000000 0.000000000 0.000000000 235.0474968 228.5224953 222.0349960 215.6124992 209.5499954 203.7799988 198.5674934 194.0149994 189.8474922 186.3049965 183.1299934 180.2024956 177.6199989 175.1824989 172.8649979 170.7399940 168.6849976 dv> p_surf 7.236700058 dv> rad 1x143x1 array of double, bsq format [1,144 bytes] dv> eangle 0.7899999823 dv> out = surfcalc(surf_em, tau,pt, p_surf, rad, eangle) 1x143x1 array of double, bsq format [1,144 bytes]
Images Plot showing the output from the example above. In a perfect world the temperature would be constant at all wavelengths.
|
DavinciWiki Mini-Nav Bar Contents
Contact Developers
All other topics
Recent Library Changes Created On: 12-31-2007 |
