radcalc
|
Contents: Citation, Description, Arguments, Usage, Examples, Images, Sub-Functions, Related Functions, Source Supported Under Version: >1.00
Citation Bandfield, J. L., and M. D. Smith (2003), Multiple emission angle surface-atmosphere separations of Thermal Emission Spectrometer data, Icarus, 161, 47-65.
Description Creates a simulated radiance spectrum given surface and atmospheric parameters. This is used for various TES atmospheric correction routines.
Arguments and Return Values Arguments: Surface emissivity, atmospheric opacity, pressure-temperature profile, surface pressure, surface temperature, emission angle. Return Value: Total radiance (in W cm-2 sr-1 cm) as measured from above the atmosphere. Optionally, transmitted and emitted components of radiance can be returned in a structure separately.
Usage Syntax: radcalc(surf_em,tau,pt,p_surf,T_surf,eangle,[wn], [struct = binary]) 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) T_surf: surface kinetic temperature (m x 1 x 1 array, in degrees Kelvin) 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 radcalc 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> T_surf
279.6499938
dv> eangle
0.7899999823
dv> out = radcalc(surf_em, tau,pt, p_surf, T_surf, eangle, struct = 1)
struct, 3 elements
rad: 1x143x1 array of double, bsq format [1,144 bytes]
trans: 1x143x1 array of double, bsq format [1,144 bytes]
em: 1x143x1 array of double, bsq format [1,144 bytes]
1x143x1 array of double, bsq format [1,144 bytes]
Images Plot showing the elements of the output structure example above:
|
DavinciWiki Mini-Nav Bar Contents
Contact Developers
All other topics
Recent Library Changes Created On: 12-31-2007 |
