act.retrievals.calculate_pbl_heffter(ds, temperature='tdry', pressure='pres', height='alt', smooth_height=3, base=5.0)[source]

Function for calculating the PBL height from a radiosonde profile using the Heffter technique. There are differences from the ARM VAP at times due to different averaging schemes. Larger differences do occur at times and are unknown as to the cause but it is being investigated and is potential a code issue with the VAP.

  • ds (xarray Dataset) – Dataset housing radiosonde profile for calculations

  • temperature (str) – The name of the temperature field.

  • pressure (str) – The name of the pressure field.

  • height (str) – The name of the height field

  • smooth_height (int) – Number of points to do a moving average on sounding height data to reduce noise

  • base (int) – Interval for pressure gridding. In testing, 5 mb was found to produce results with the lowest RMS


ds (xarray Dataset) – xarray dataset with results stored in pblht_liu_liang variable


Heffter JL. 1980. “Transport Layer Depth Calculations.” Second Joint Conference on

Applications of Air Pollution Meteorology, New Orleans, Louisiana.

Sivaraman, C., S. McFarlane, E. Chapman, M. Jensen, T. Toto, S. Liu, and M. Fischer.

“Planetary boundary layer (PBL) height value added product (VAP): Radiosonde retrievals.” Department of Energy Office of Science Atmospheric Radiation Measurement (ARM) Program (United States) (2013).