@@ -20,4 +20,12 @@ To read in and work with the HDF4 files, it is helpful to create a dedicated pyt
* ethz: plots of cloud-radiative heating for a tropical cyclone
* pythonenv: create dedicated python environment
* core: functions to read-in data and compute heating rates per granule
* solar: functions to derive daily-mean insolation, adapted from climlab of Brian Rose (see the readme file in the solar directory)
To generate the binned heating rates, run ``make_binned_heating.py`` followed by ``postprocess_binned_heating.py``. The first script generates heating rates for each year, the second merges the year and takes care of the time and vertical axes, plus adds some metadata.
``how_to_best_compute_radheating.ipynb`` illustrates that heating rates should be diagnosed by using the flux divergence divided by the pressure level thickness - this is the approach implemented in the core function.
``plot_zonaltimemean_crh.ipynb`` does a time-mean zonal-mean plot of cloud-radiative heating - this is a sanity check.
