Calculation of heating rates. Output is only provided at altitudes specified by zout. To get heating rate profiles a number of altitudes must thus be specified by zout. Heating rates is the change of temperature with time in units of K/day. For spectral calculations the default output is a matrix: 0.0 zout1 zout2 … \Verbatimlambda1 heating_rates … lambda2 . . . . . For integrated calculations (output sum or output integrate) the default output is in two columns with column 1 being the altitude and column 2 the heating rates. The output of heating_rate can also be specified with the output_user option. Note that heating rates are only well-behaved up to altitudes for which the respective correlated-k options are valid. E.g. about 60 km for fu and about 80 km for kato, kato2, kato2.96, and lowtran. Attention: For spectral calculations, the extraterrestrial spectrum is assumed to be in mW/(m2 nm).
Two different methods are implemented to calculate the heating rate, which can be selected
with an optional keyword:
heating_rate [method]
where method may be either layer_cd (heating rates are derived from centered
differences of the flux, this is the default method) or local (heating rates are derived
from the actinic flux). Attention: heating_rate local introduces new levels into the profile which
slightly affects the model output.
There is also a third method called layer_fd, which means that heating rates are derived
from forward differences of the flux over one layer. Please be aware: If using layer_fd, the output is
not representativ for a level, but for the layer from the z-level of the line in the output file,
where it is written, up to next output level above!