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 ...
lambda1    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 that the output is not a level property, but representativ for the layer from the z-level of the line in the output file, where is written up to next output level above!





Arve Kylling 2010-03-10